• / 5
  • 下载费用:1 下载币  

叠前地震同步反演方法及应用

关 键 词:
地震 同步 反演 方法 应用
资源描述:
第 41 卷 第 2 期 成都理工大学学报 ( 自然 科学版 ) 1 2014 年 4 月 F F 01410.3969/j. 671014.02.04 [ 文章 编号 ] 16712014) 02收 稿 日期 ] 2012基金项目 ] 中海油能源发展科技项目 ( 第一作者 ] 方中于 ( 1971 - ) , 男 , 高级工程师 , 从事地球物理处理及方法研究 , 前地震同步反演方法及应用方 中 于1, 2, 朱 江 梅1, 2, 李 勇1, 黄 丽 娜2, 李 佳 胜1( 1. 成 都 理工大学 地球物理学院 , 成都 610059; 2. 中海油能源发展物探技术研究所 , 广东 湛江 524027)[ 摘要 ] W 地区珠海组储层砂岩内部物性 、含油气性变化快 , 开展储层物性预测尤为重要 。叠前地震同步反演方法是利用不同炮检距道集数据和纵波 、横波 、密度等测井数据的反演方法 ,可同时反演得到储层岩性 、物性相关的多种参数 , 用于综合判别储层物性及含油气性 。将叠前地震同步反演方法应用于 W 地区珠海组储层的物性及含油气性预测中 , 利用获得的 P 波阻抗 、S 波阻抗及密度等地震特征参数对珠海组储层进行了综合判别 , 并对珠海组储层品质作了评价 , 优选出有利优质储层及其分布范围 , 并被钻探证实 。[ 关键词 ] 珠海组 ; 叠前 ; 同步反演 ; 物性 ; 综合判别[ 分类号 ] 4; [ 文献标志码 ] 区珠海组储层砂岩内部物性 、含 油 气性变化快 , 这些变化制约了该地区的开发 , 因此开展储层物性预测尤为重要 。叠前地震反演的目的是为了获得可靠的纵波速度 ( 波阻抗 ) 、横波速度( 波阻抗 ) 及密度等弹性参数 , 以预测储层的岩性特征和流体特征 。这个问题在国外有多位专家进行过讨论 , 1996 年对 P 波反射系数 ( RP) 、横 波 反射系数 ( RS) 和 密 度反射系数进行了反演研究[ 1], 2003 年 使用 了 演方法[ 2]。R. 人建立了叠前地震反 演 新方法[ 3, 4], 它直接反演纵 波 阻抗 ( 速度 ) 、横波阻抗( 速度 ) 及密度 , 并试图将其方法扩展到转换波数据反演方面 。叠前地震反演在国内有专家开展了一定程度研究[ 5 -8], 但 仍 处于单参数反演阶段 。由单纯叠后反演转入叠前反演 , 提高了反演结果的可信度和精度 , 并取得一些成果 。叠前地震同步反演方法是利用不同检距道集数据以及测井数据的叠前 、叠后联合反演方法 , 可同时得到纵横波速度 ( 波阻抗 ) 及密度等参数 , 这种新方法在反演过程中考虑了 vP、 度之间的关系 , 提高了储层岩性和流体的识别能力 。本文将叠前地震同步反演方法应用于 W 地区珠海组储层的物性及含油气性预测中 , 利用获得的地震特征参数对珠海组储层进行综合判别 , 取得了较好的效果 。1 储层地震特征参数叠前同步 反 演方法叠 前 地震同步反演方法是一种基于模型的叠前反演方法 , 叠前地震同步反演可获得可靠的 P 波速度 ( 、S 波 速 度 ( 和 密 度 ( ρ) 等参数 , 以达到直接预测储层的岩性特征和流体特征的目的 。1.1 叠前地震同步反演中的正演模型在零偏移距情况下 , 反射系数 R 式可写为[ 4]Ri=12Δi=12i+1- )i( 1)式 中 : 波 阻 抗 , RP, i 层 和第 i+1 层之间分界面的反射系数 。式 ( 1) 可写成矩阵形式RPN=12- 1 0 0 …0 - 1 1 …0 0 - 1 …   …PN( 2)式中 : i= i) 为 P 波阻抗的自然对数 。褶积模型可表示为下面的矩阵形式N= …   …RPN( 3)式中 : 震道的第 i 个样值 , 震子波的第 j 个样值 。联合 ( 2) 式和 ( 3) 式 , 可得叠前地震同步反演中的正演模型 , 即T =124)式 中 : W 为 式 ( 3) 中的地震子波矩阵 , D 为式 ( 2)中的导数矩阵或微分算子 。1.2 叠前地震同步反演方法J. L. G. C. 人研究的基础上 ,将 式改写为[ 9]Rθ) = 5)式 中 : 1 + - 8r = vS/3= - 0.5 2θ 为 入 射角 。式 ( 5) 中 P 波 、S 波和密度反射系数可表示为RP=12Δρ( )ρ( 6)RS=12Δρ( )ρ( 7)RD= Δρ/ρ ( 8)联 合 式 ( 4) 与式 ( 5) , 得到T( θ) =12θ) 2θ) ( θ) 9)式 中 : ρ , W( θ) 依赖 于 入射角度 。考虑到 的关系k Δ10)m Δ11)则式 ( 9) 变 为T( θ) =珘θ) θ) ( θ) 12)式 中 :珘222系数 K 和 m 可 由 汇 和 得 。式 ( 12) 可写成矩阵形式T( θ1)T( θ2)T( θN)=珘θ1) W( θ1) θ1) W( θ1) D θ1) W( θ1) θ2) W( θ2) θ2) W( θ2) D θ3) W( θ3) D  珘θN) W( θN) θN) W( θN) D θN) W( θN)LD( 13)由 式 ( 13) 可 知 , Δ为了 解 式 ( 12) , 求得 Δ首 先 建立初始化的解 。[ = [ ]T( 14)式 中 然后 , 采 用共轭梯度法即可求解 。叠前地震同步反演主要技术环节包括 : 测井岩石物理分析 , 层位标定 , 分角度子波提取 , 不同子波均衡角度部分叠加数据体之间振幅 、频率和相位差异 , 叠前初始模型建立以及叠前同步反演等 。叠前同步反演流程如图 1 所示 。2 W 地区珠海组储层地震特征参数剖面与可靠性评价图 2 是应用地震叠前同步反演方法所获 得 的珠海组储层地震特征参数剖面 , 它们分别是 叠前同步反演流程图 51·第 2 期 方 中于等 : 叠前地震同步反演方法及应用图 2 珠海组储层地震特征参数剖面 of A) 连 井 叠前 P 波阻抗预测剖面 ; ( B) 连井叠前 S 波阻抗预测剖面 ;( C) 连井叠前 S 波密度预测剖面 ; ( D) 连井叠前 vP/ 剖面阻 抗 ( 速 度 ) 剖面 、S 波阻抗剖面 、密度剖面及 vP/ ( 图 中色标标出了各参数值的变化范围 ) 。图 3 是叠前地震同步反演的 P 波阻抗 、S 波阻抗 、密度与验证井 井相应曲线对比图 。图中红色曲线和蓝色曲线分别为反演的曲线和测井曲线 。由图 3 可看出 , 反演的曲线在形态 、变化特征和分辨率等方面均可与测井曲线相对比 , 具有较好的一致性 , 其相关系数最低值为 0. 820, 最高值为 0. 940。与 9 测井参数值及变化范围对比可以发现 , 珠海组储层P 波阻抗值变化范围为 9 700 ~10 450, 表现为中高值 ; S 波阻抗值变化为 6 350 ~5 200, 表现为中高值 ; 密度 < 2. 5, 其变化范围为 2. 42 ~ 2. 45, 密度偏低 ; vP/. 70, 其 变 化范围为 1. 60 ~1. 66,vP/ 。这 种储层地震特征参数特征揭示了珠海组储层的品质 。3 W 地区珠海组储层预测图 4 是 在 井控 ( 下将反演的地震特征参数剖面 ( 图 2) 转换成珠海组岩性剖面 。对比图 4 与图 2 可看出 , 珠海组储层砂岩具有相对高的横波速度 、低密度和低vP/ 特征 , 储层砂岩发育 , 形成了良好的储集空间 。图 5 是利用图 4 与图 2 预测出的珠海组有利储层剖面 。依据在该地区砂岩越发育储层品质越好 , 在图 5 上可划出珠海组储层有利优质储层段 ,如表 1 所示 。由表 1 可知 , 珠海组第二段有利优质储层层段多于珠海组第三段的有利优质储层层段 。在 和 珠海组有利优质储层段日产天然气分别为 0. 72 ×106. 46 ×106也 为 井证实 。表 1 珠海组有利储层段表 of 上的层位时间 /利储层段范围 ( 备 注珠 海组第二段( 组 )2388 ~2418240 ~330, 390 ~535, 540 ~660,670 ~810, 1240 ~13801240 ~1380储 层 结构复杂珠海组第三段( 组 )2640 ~2682 200 ~295, 490 ~690, 1230 ~1380490 ~690储 层 结构复杂·451· 成都理工大学学报 ( 自 然 科学版 ) 第 41 卷图 3 叠前反演的纵波 、横 波 阻抗和 上纵波 、横波阻抗比较图 of 9 珠海组储层岩性预测剖面 of 珠海组有利储层预测剖面 51·第 2 期 方 中于等 : 叠前地震同步反演方法及应用4 结 束 语叠前 地震同步反演方法是一种新的反演方法 , 可同时反演多种参数 , 具有多功能的特点 , 这是对叠前反演的一种新发展 。叠前同步反演方法又是一种联合反演方法与技术 , 它充分发挥了叠前 、叠后的各自优势 , 增强了反演方法的适应性和稳定性 , 有效地提高了所反演参数的可靠性和精度 。应用这种反演方法获得的地震特征参数对W 地区珠海组储层的岩性和有利储层段进行的预测与评价 , 其预测结果与实钻结果比较吻合 , 揭示出储层的砂岩发育与储层品质之间的关系 , 为W 地区珠海组储层预测提供了技术支撑和科学依据 。[ 参 考 文 献 ][ 1] L, M. VO J] . 996, 61( 6) : 1575 -1588.[ 2] , . VO J] . 2003, 68( 1) : 185 -198.[ 3] R B. C] / /6 2004: 1 -4.[ 4] D P. of C] / / 005: 1635 -1638.[ 5] 唐 晓 花 , 成德安 , 吕金龙 . 叠前同步反演在徐家围子断陷火山岩气藏预测中的应用 [ J] . 石油物探 , 2009,45( 3) : 285 -289. H, A, L. of J] . 2009, 45 ( 3) : 285 -289. ( 6] 张 卫 卫 , 颜承志 , 庞雄 , 等 . 地震波形联合叠前反演技术在白云深水区的应用 [ J] . 石油天然气学报 , 2011,33( 12) : 71 -78. W, Z, , et J] . il 2011, 33( 12) : 71 -78. ( 7] 刘 涛 , 方 中于 , 但志伟 , 等 . 叠前同步反演在涠西南凹陷储层预测中的应用 [ J] . 石油天然气学报 , 2011, 33( 9) : 77 -81., Y, W, et of of J] . il 2011, 33( 9) : 77 -81. ( 8] 蔡 涵 鹏 , 贺振华 , 李瑞 , 等 . 两项 演研究及应用[ J] . 成都理工大学学报 : 自然科学版 , 2011, 38( 5) :581 -587. P, H, , et VO J] . ,2011, 38( 5) : 581 -587. ( 9] L, C, J, et of VO a 3-D J] . 1994, 59( 9) : 1362 -1376.on 2, 2, 10059, . 24057, It is to of to of is of on P In to of be in is to to of is on S of of by in in to of ey 51· 成都理工大学学报 ( 自 然 科学版 ) 第 41 卷
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:叠前地震同步反演方法及应用
链接地址:http://www.oilwenku.com/p-63590.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开