地球物理学进展  2016, Vol. 31 Issue (5): 2108-2115   PDF    
叠后地震反演面临的问题与进展
贾凌霄1,2, 王彦春3, 菅笑飞4, 吕鹏1,2, 田黔宁1,2, 谢玮3     
1. 中国地质调查局地学文献中心, 北京 100083
2. 中国地质图书馆, 北京 100083
3. 中国地质大学(北京)地球物理与信息技术学院, 北京 100083
4. 华能国际电力开发公司, 北京 100031
摘要: 叠后地震反演技术自提出以来,在储层预测和评价、油藏特征描述等方面得到了广泛的应用.叠后地震反演可以相对快速地将地震波形信息转换为具有地质意义的波阻抗等信息,指导储层的识别与刻画.因此叠后地震反演技术在油气勘探中发挥了十分重要的作用.本文首先就叠后地震反演技术的起源、分类做了简要介绍;然后对影响反演结果的因素进行了归纳,总结出叠后地震反演结果主要受地震频带宽度、子波、地质层位和井信息这4种因素的影响,并系统梳理总结了每种影响因素在叠后地震反演中发挥的作用和目前在优化反演结果中所用的方法和技术;接着分析了反演过程中构建的低频和高频分量是造成叠后地震反演多解性的主要原因,介绍了相关研究成果;最后点明了现阶段叠后地震反演技术存在的问题并展望了其发展趋势.
关键词叠后地震反演     地震频带     子波     地质层位     井信息     低频分量     高频分量     多解性    
Problems and progress in post-stack seismic inversion
JIA Ling-xiao1,2 , WANG Yan-chun3 , JIAN Xiao-fei4 , LYU Peng1,2 , TIAN Qian-ning1,2 , XIE Wei3     
1. Geoscience Documentation Center, China Geological Survey, Beijing 100083, China
2. National Geological Library of China, Beijing 100083, China
3. School of Geophysics and Information Technology, China University of Geoscience, Beijing 100083, China
4. Huaneng International Power Development CO., Beijing 100031, China
Abstract: Post-stack seismic inversion has been widely used in the reservoir prediction and evaluation, reservoir characterization since it was invented. Post-stack can relatively quickly transform seismic waveform information to geological information such as impedance, which can guide the identification and characterization of reservoirs. So post-stack seismic inversion plays a significant role in the oil and gas exploration. This paper firstly concisely introduces the origin and classification of post-stack seismic inversion; then sums up the influential factors of the inversion results, summarizes the post-stack seismic inversion results mainly influenced by seismic frequency band width, the wavelet, horizons and well information, and summarizes the effects caused by these factors in post-stack seismic inversion and the method and technology used to optimize the inversion results in recent years; the paper also analyzes the main reasons of the ambiguity in inversion are the low frequency component and the high frequency component needed in the inversion process and introduce the related research achievements; finally pointing out the current problems exiting in post-stack seismic inversion, which include:①the effects of noise is not considered in the convolution formula model, ②the generalized elastic impedance has not been wildly used in the practice, ③ambiguity problems in the process of building high frequency model components has not been solved, ④the nonlinear inversion methods currently has not wildly applied in actual large-scale productions; in addition, the paper prospects the development trend of post-stack seismic inversion.
Key words: post-stack seismic inversion     seismic frequency bandwidth     wavelet     horizon     well information     low frequency component     high frequency component     ambiguity    
0 引言

地震反演是利用观测的地震资料,加上已知的地质规律和钻井、测井资料为约束,对地下岩层空间结构和物理性质进行成像(求解)的过程,广义的地震反演包含地震数据处理解释的整个内容.

地震反演技术最早可以追溯到20世纪70年代.国外一些学者在基于测井声波阻抗正演合成地震数据的基础上提出了Seislog反演理论(Lindseth et al., 1979Becquey et al., 1979), 该方法利用速度求取波阻抗的低频分量,由地震数据求取波阻抗的高频分量(反射系数),通过两部分求和获得波阻抗反演结果.该反演理论的产生使人们开始利用波阻抗信息(更接近岩石性质)来认识和解释地震数据中的地质问题.波阻抗具有明确的物理意义,可以对储层特征进行定量和确定性的描述.李庆忠院士指出:“波阻抗反演是高分辨率地震资料处理的最终表达方式”.这说明波阻抗反演在地震反演技术中的特殊地位.

地震反演通常分为叠前和叠后反演两大类,叠后反演侧重于对储层的识别,而叠前反演侧重于对流体性质的判别(陈红等,2012;程冰洁等,2012姜岩等,2013强敏等,2015).在勘探开发初期和测井资料不全的情况下,叠后反演以其相对简单的算法和较高的稳定性,在实际储层预测工作发挥了重要的作用,是众多地球物理学家的研究重点.

1 叠后地震反演的分类

随着油气勘探对地震反演技术的需求和计算机水平的不断进步,叠后地震反演技术得到了迅速的发展,并广泛应用于实际生产中,成为储层预测中不可或缺的技术手段.叠后地震反演依据不同的标准有不同的分类,不同的分类又有不同的应用目的和研究重点.按测井资料在其中所起作用的大小可将叠后地震反演分成4类:地震直接反演、测井控制下的地震反演、测井-地震联合反演和地震控制下的测井内插外推(姚逢昌和甘利灯,2000).早期的叠后地震反演是靠线性算法实现的,现今线性反演由于其运算速度快,较为稳定等优势仍有广泛应用.20世纪80年代,随着人工神经网络等非线性学科的迅速发展,人工神经网络算法、遗传算法、模拟退火算法、小波变换算法等(陆文凯等,1996杨利强,2003孙思敏和彭仕宓,2007)逐渐应用到叠后地震反演中,因此叠后地震反演从算法上还可分为线性反演和非线性反演.自地质统计学反演被提出后(Hass and Dubrule, 1994),按反演的实现过程叠后反演还可分为确定性反演和随机反演(潘昱洁等,2011),确定性反演主要受地震数据的约束;而随机反演应用随机理论和地质统计学方法,综合井和地震数据,可产生多个等概率的反演结果.根据反演过程中处理的道数,叠后地震反演还可分为单道反演和多道反演.大多数反演算法是基于单道进行计算的,一些学者提出的多道反演方法在一定程度上可以提高反演结果的信噪比(李宏兵,1996林小竹等, 1998, 齐彦福等,2015).

目前国内比较主流的分类方法是依照实现方法将叠后地震反演分为三类:道积分、递推反演、基于模型的反演.道积分和递推反演属于直接反演法,均比较完整地保留了地震反射的基本信息(断层,产状),不存在多解性,但受地震频带宽度的限制,分辨率较低.基于模型的反演技术突破了传统意义上的地震分辨率的限制,理论上可以得到与测井资料相同的分辨率,但反演结果中测井信息和模型提供的高低频分量导致了反演的多解性,且受制于钻井的数量和井网的分布.

2 叠后地震反演影响因素 2.1 地震频带宽度

原始地震数据的分辨率达不到储层预测的要求,因此人们希望通过叠后反演方法来提高地震数据的分辨率.但任何反演方法都离不开地震数据,因此地震数据对叠后地震反演结果有着重要影响.特别是对于直接反演法,反演结果的分辨率直接依赖于地震数据,因此在此类方法中地震数据的频带直接决定了反演结果的分辨率.在基于不同频率子波的正演模型数据基础上采用有限带宽反演法、广义线性反演法和最大似然反演法进行叠后反演试验,发现高频数据反演出的砂体边界更加精确、且反演误差要小于低频数据(凌云等,2008).

在地震数据处理中人们提出了很多拓宽地震频带宽度的方法,例如预测反褶积、时变反褶积、最小平方反滤波、反Q补偿地震波传播衰减、地表一致性反褶积等方法(Robinson,1957Robinson,1967Clarke,1968Levin,1989),这些方法对提高地震数据的成像分辨率方面起到了重要作用.但实践表明,提高地震频带是以牺牲信噪比为代价的,过度的追求更宽的地震频带,得到的反演体中的高频成分有可能混杂了大量的随机噪声,从而降低反演结果的准确性,给储层预测带来误差.因此李庆忠院士提出反演前应确定地震资料的“有效频宽”,同时建议采用分频扫描的方法来检验信噪比大于1的频率带宽,在“有效带宽”内的地震资料才能用于反演(李庆忠,1998).

2.2 子波 2.2.1 子波相位

叠后地震反演主要通过井震标定中尽量使合成地震记录和实际地震记录获得较大的相关系数为前提来提取子波,子波选取的合理与否直接影响到最终的叠后反演结果.通常认为,子波的相位会对反演结果产生较大的影响,在带宽相同的条件下零相位子波的旁瓣比其他类型相位子波的要小,其能量更为集中,因此具有更高的分辩能力;在地层较薄时,相对于其他类型的子波例如最小相位子波无法分辩的相邻尖脉冲,零相位子波往往可以分辩,且零相位子波更有利于鉴别地震剖面的极性.但在实际中,地震数据经过提高分辨率处理、非线性动校拉伸处理和消除偏移频散影响等处理后,叠后地震数据仅满足混合相位子波的条件,并没有理想的方法将其处理为满足反演要求的零相位.尽管多数地球物理工作者认为地震子波的非零相位问题仍是反演结果不精确的主要影响因素,但利用实际数据通过实验,发现统计意义下的混合相位子波相对于零相位的统计子波得到的合成地震记录与井旁地震道有更高的相关系数(晓宇等,2009).事实上,大多数实际工作中提取的最终用于反演的子波为混合相位子波,证明叠后地震成像数据更满足混合相位子波的条件.

2.2.2 子波提取

波阻抗反演技术是建立在褶积模型基础理论上的(Lindseth,1978),因此子波的提取是反演中最重要的技术环节之一.一般情况下工区内都有多口井,不同的井可以提取不同的子波.在特定测线或井控区块内提取的子波往往只在限定区块内可以获得较好的反演效果,但不适用于整个工区.若分区块采用不同的子波,则有可能导致不同区块反演结果不闭合的问题,进而对储层解释带来错误,所以在一个工区进行波阻抗反演时一般只提取一个子波.传统的子波提取方法有三类:一是直接从野外观测数据中提取;二是先选取标准的子波模型(如Trapeziod、Richer、Klauder、Butterworth等),再将从地震数据中得到的频宽,主频,相位等参数值整合进来;三是在由井得到的合成记录与井旁地震道最佳匹配的条件下分别提取每口井的最佳子波,再由这些井计算出一个平均子波,这在理论上是最佳的子波提取方法,在实际中也有应用,但在正、反演试验中并没有获得理想的结果,因为简单的在时间域或频率域的相加平均运算没有考虑不同井信息的准确度和储层的横向变化,因此不能真正综合各个子波的特征.

子波的判别标准就是与井得到的反射系数褶积,通过求得与地震记录的相关性来判别子波的优劣.在工区面积不大,储层横向变化小,井网分布较为均匀的工区使用各井提取的子波,分别在每个井点与反射系数序列褶积,与该井旁地震道求相关运算,再将每口井得到的相关系数相加,总和最大者就可以认为是最佳的地震子波.但在此之前要做好井的环境校正、时深转换、标准化等前期工作(许升辉和马劲风,2003).在工区内井较多时可以利用神经网络等非线性方法将井阻抗作为神经网络的输入,将井旁道作为输出,在合成地震道与井旁地震道得到最佳拟合时完成网络训练.对每一口井都进行多级地震子波提取,所有井的子波以神经元子波的形式存在于网络中,在此过程中还可以不断加入新井,以完善子波(崔炯成等,2007).传统的地震子波提取方法对噪声的适应性较差,在低信噪比地区,还可以采用广义线性反演提取子波等方法,在传统的褶积公式中加入地震噪声和井误差的影响,可以获得较好的子波估计.

2.3 地质层位

对地震资料进行构造解释得到的地质构造层位是叠后地震反演约束模型的控制资料,是地质信息和测井信息的重要载体和框架(孟宪军等,2004).在地层比较平缓,构造相对平稳的地区,仅将顶部与底部层位作为约束层位与在其基础上增加其他约束层位相比对反演结果的影响不是很大,但储层顶底层位的准确度十分重要(Sancevero et al., 2005).有学者提出了构造建模约束地震反演技术(李占东等,2013),在常规地质解释层位的基础上,加上井-震联合小尺度构造解释,在常规地震层位之间针对小尺度地质体开展更为精细的构造解释,建立更为精细的初始模型,结果表明此反演方法可以明显提高砂泥岩薄互层条件下储层预测的精度(图 1).

图 1 同一条测线构造建模约束前后地震反演效果对比 (a)未参与构造建模过井反演剖面; (b)构造建模约束反演过井反演剖面. Figure 1 Comparison of seismic inversion results between before and after structure modeling being taken as a constraint (a) Inversion profile of not involved in structure modeling; (b) Inversion profile of involved in structure modeling.
2.4 井信息 2.4.1 井曲线

大多数叠后波阻抗反演方法都要依赖于测井信息.通常情况下,一个地区储层和非储层的声波时差有较大的差异,但由于地下储层具有非均质强、储集空间复杂、控制因素多等特点,导致在某些地区例如扇三角洲沉积区测井信息在储层和非储层的响应区分不明显.若直接使用原始的声波曲线进行反演,反演结果即使分辩率达到了要求,也无法有效识别出储层.目前解决上述问题基本上采用的是测井曲线重构的方法,重构方法又分为整合其他受环境因素影响小对储层敏感的测井曲线制作“拟声波曲线”和利用岩石物理模型重新正演出声波曲线.

目前“拟声波”法在实际应用中能取得较好效果的大约有种:一是将原始声波曲线低通滤波作为低频部分,将对岩性反映敏感的自然伽马或自然电位等曲线作为高频部分,最后将低频和高频合并得到拟声波曲线(许正龙,2002杨少虎等,2006),此方法简单快速,得到的拟声波曲线可增强对储层的区分能力,但丧失了一部分原始的高频信息.二是先找到不同岩性地层在各种测井曲线上的统计特性,再将多条测井曲线按不同的权值加入声波曲线,这样做不仅保留了声波曲线的原始频率成分,还添加了更加丰富的频率信息,但算法较为复杂(贺懿等,2008王绍忠和周红科,2009).

岩石物理正演技术在叠前反演的应用较为普遍(张永江等,2010周水生等,2012宗兆云等,2012陈怀震等,2014),但也可用于叠后反演.根据工区实际情况选择合适的岩石物理模型(Kuster-Toksz模型,Wood模型,Gassumann模型,Xu-White模型等),通过反复试验岩石骨架和流体的弹性参数,找到其与原始曲线的最佳拟合点,对声波曲线进行重构,重构的曲线相对于原始曲线可以更好地反映储层的含油气特征,并且得到的合成记录与地震记录的相关性也更高,真正实现了测井曲线的标准化.但要注意岩石物理模型的选取要适合工区的实际情况,在工区没有岩心数据时要谨慎使用此方法(黄伟传等,2007龙国清等,2008).

2.4.2 井的空间位置和井的约束

叠后地震反演一般利用测井信息对反演结果进行约束,如何更好地利用有关井资料的先验信息是叠后反演研究的重要课题.一般认为参与约束的井的数量越多,越有利于提高反演的精度,但理论研究表明在实际工作中我们要充分考虑不同反演方法的特点和井点的位置.对于确定性反演来讲,在构造与沉积背景基本相同的情况下,一般来讲只需一口井作为约束即可得到正确的相对阻抗关系,多井约束反而会破坏储层波阻抗的空间相对关系,出现波阻抗突变条带等现象,从而影响反演结果的正确性.这是因为环境因素和测井年代不同造成了井之间的波阻抗差异,另外仅通过数学插值方法难以描述地下储层的变化.图 2a为由两个非规则边界的砂体和一条河道组成的地质模型的时间切片,且标注了用于反演的两个井(1井和2井).图 2bc为分别用井1和井2单独约束下的纵波阻抗确定性反演时间切片,结果表明两个结果波阻抗的相对变化关系一致,且与初始模型吻合较好.图 2d为1井和2井同时约束下的反演结果,由于多井之间真值的差异,在通过井资料建立低频模型时在空间插值的过程中引起了假象(时间切片左下角的低值区域).因此在确定性反演时要慎重考虑是否用多口井进行约束,如果必须要用多口井进行约束,一定要对井进行环境校正或标准化处理.

图 2 初始模型与用不同数量的井约束下的确定性反演结果 Figure 2 The initial model and deterministic inversion results of using different numbers of wells to constrain

不同于确定性反演需要建立一个初始的低频背景模型,随机反演直接通过随机扰动反演波阻抗,反演过程中直接利用井的信息,故井的信息越多,先验信息就越丰富,随机反演结果的不确定性就越小,与真实情况也越接近(潘昱洁等,2011).有学者通过试验得出随着井网密度的增加,地质统计学随机反演的结果有明显的改善,当垂直井距缩小为一定距离时,丼间的反演结果基本稳定(董奇等,2013).可见不同的叠后反演方法对井约束的要求是不同的.

3 低频和高频分量的构建

叠后波阻抗反演技术是基于褶积模型的,由于受地震子波带限性质的影响,为了拓宽反演结果的频带,提高反演结果的分辨率,必须加入地质构造框架、地层沉积模式和测井资料等先验约束信息来构建低频和高频分量.低频与高频成分在地震子波通带范围之外,低频成分与高频成分与地震子波的谱是互补的(李庆忠,1993),因此从褶积模型中无法检验反演结果的正确性,这也是造成叠后地震反演多解性的主要因素之一.因此构造合理的相对精确的低频和高频分量十分重要,现阶段一般通过约束模型作为载体对地震反演结果进行高低频分量的补充.

3.1 低频分量的构建

实际生产中,低频分量主要是从测井资料与叠加速度中求取的,但从测井资料获得低频分量精度要高一些.井点处的低频分量可以根据井信息准确地得到,但在井点以外,只能靠内插和外推得到,因此井点以外低频分量的获取存在多解性,其正确与否是最终反演结果成败的关键.

早期的波阻抗反演方法只能解决单井控制的问题,低频分量采用沿地震剖面水平横向递推的方法得到,不受解释层位的约束,各道补偿的低频信息是一样的,不考虑地层产状、岩性横向变化、断层等构造的影响.80年代中期,西方地球物理公司和Sierra公司在其推出的商业软件中提出了利用构造控制,顺层递推波阻抗低频分量的方法,可应用于地层产状有横向变化和岩性横向稳定的情况.CGG公司在ROVIM软件中考虑了丼间砂层的横向变化,并提出了能有效控制横向砂层变化率的方法.有学者也提出了通过砂层的小层对比进行横向插值的方法(李庆忠,1989凌云等,2008).

实际储层往往地质条件复杂、断层和井的数量多、岩性岩相和厚度横向变化大,因此在解决实际生产问题时还应当考虑更多的因素.在构建低频分量前,要做好精细的井震标定工作,使地震资料与测井资料达到最佳的匹配,在工区存在多口井的条件下还要做好井的标准化,避免因外部环境或仪器误差导致的岩性变化假象.反演中在知道反射系数的情况下,可按下式求得第n层的波阻抗值, 公式为

(1)

式中:ρnvn为第n层的波阻抗值,Ri为第i层的反射系数,可以看出初始波阻抗值ρ0v0对反演结果至关重要,初值选择错误不仅会导致低频分量的错误,也会产生不正确的相对波阻抗序列(马劲风等,1999).因此构造低频信息时,要选择全剖面或全工区岩性稳定的层位,即ρ0v0值横向稳定的层作为横向递推使用的起始层位.

在两井间内插时,还应根据地质资料给出各井控制范围的权函数来确定每口井的控制范围.当两点外遇到岩性突变时,仅利用井点资料是难以内插外推出岩性突变信息的,可以借助速度谱或其他地质信息来确定岩性横向突变的幅度与范围.但一般的速度谱资料中的速度不具有明确的地质信息,分层能力较差,局部速度谱低频分量的水平变化需要通过加密速度谱来证实,因此要谨慎使用.另外,用测井资料得到的低频成分一定要与子波的谱互补衔接,如果子波的频率与低频成分的频率存在缺口或重叠,反演的阻抗曲线会在缺口或重叠处出现震荡,为反演结果带来误差(马劲风,2000).

3.2 高频分量的构建

由于地震和子波带限的性质,高频信息不可能通过褶积模型来做合理的检验;相反,随机噪声的频率一般较高,在噪声大于信号的高频频带里,越用褶积模型来做约束和检验,多解性越强.有学者通过理论模型证明了有效频带(信噪比较高区域)之外的地震高频信息永远是多解的.李庆忠院士提出约束叠后地震反演的前提是地震剖面中各种噪声和干扰都已作了最大的压制,建议采用分频扫描来检验信噪比大于1的频率带宽来确定“有效频宽”,根据有效频宽的频率门槛值设计高频滤波算子,用此滤波算子对测井波阻抗曲线滤波得到波阻抗曲线的高频分量用来构建高频模型.但此方法也只是一种合理的“猜测”,准确性难以检验,只能根据地质信息来作判断.此外,约束方法的正确选取也可以改善反演结果,理论模型分析表明井段的最大最小值、测井曲线范围包络、中值滤波等约束策略均能提高反演结果,其中中值滤波可以减弱噪声的干扰,合理联合运用多种约束可以达到更好的效果.

叠后波阻抗反演要想得到准确的高频信息,必须施加合理的约束条件.一般反演求得的波阻抗解应同时满足求得的解与已知井的波阻抗值方差最小和求得的合成地震道与实际地震道的方差最小.这两个条件可以用一个L2范数表示, 公式为

(2)

式中:k1k2是权系数,k1+k2=1,σs(x)是求得的波阻抗解,σw(x)是井的波阻抗值或相邻地震道的波阻抗值;fs(x)是实际地震记录,fm(x)是合成地震记录.理想情况下两个范数可以同时达到最小值,但由于受环境、测量方法和仪器本身的影响,一般情况下井的合成记录与井旁道地震道的相似度比较低.因此要保证反演结果的真实性,高频约束应为弱约束,即权系数k1一般要远小于k2.

正反演试验表明测井约束的波阻抗反演方法对地震资料有效频率之外的高频分量修正能力很弱,反演结果中的高频分量主要来源于初始模型的高频分量,但稀疏脉冲反演由于在目标函数中引入稀疏度准则,反演过程本身能够产生地震数据和初始模型都不存在的高频分量,若联合应用上述两种方法可改善反演效果(马劲风等,2000李国发等,2010),此方法在胜利油田某区块中也有较好的应用(刘成斋等,2002).

4 叠后地震反演存在的问题及发展趋势

自叠后地震反演提出的30多年来,该技术已经获得了长足的进展,新方法不断涌现,反演方法从最初的直接反演和线性反演,发展到后来的随机反演、非线性反演、深度域反演等(张静等,2010),反演的精度和分辨率也逐步提高,在油田开发,储层预测,理论研究等领域发挥了重要的作用.但叠后反演依然存在如下一些问题:

(1)褶积模型公式中没有考虑噪声的影响.但在实际地震记录中都存在噪声,虽然地震数据处理人员提出了很多压制噪声的方法,但仍然不能将噪声从地震数据中完全消除.

(2)反射系数递推公式是在建立在反射界面水平、内部介质均匀且各项同性、地震波垂直入射的假设条件基础上的,但实际地震道往往不满足这些假设条件.叠后地震资料并非自激自收的地震记录,而是共反射点道集反射振幅叠加平均的结果,因此共反射点地震资料也存在AVO问题.运用弹性阻抗不适合于共中心点叠后地震剖面.为此有学者提出了广义弹性阻抗的概念,通过理论试验发现可以较好地解决较大炮检距条件下的P波叠后资料的正反演问题(马劲风,2003王浩等,2009王浩等,2010),但目前尚没有在实际中得到广泛的应用.

(3)褶积模型对于地震资料的高频成分是没有约束力的,因为地震道是由带限的子波与反射系数序列进行褶积的结果,褶积运算就是滤波运算.对反射系数褶积后就滤掉了反射系数中的高频成分.因此想要反演出准确的高频信息,必须要有正确的高频资料作为约束,目前常用的高频资料仅为测井资料,但此约束在井点处较为可靠,在离开井点的位置是存在多解性的,这种多解性在当前的技术条件下是很难解决的.

(4)线性反演算法简单易行,运算速度快,但容易造成局部寻优,非线性反演方法大都采用随机搜索算法而不沿着某条确定的路径进行搜索,因而可以实现全局寻优.但非线性反演算法计算量较大,收敛速度慢,目前在实际中还没有得到广泛的应用(马劲风,1999马劲风等,2002张永刚,2002魏超等,2008).

总体来说,提高反演结果的精度,拓宽反演结果的频带,降低多解性,对特殊地质问题作出合理的解释,是叠后地震反演所追求的目标.相信随着地球物理技术的不断发展,叠后反演方法在理论和实际应用上定将有新的突破.

5 小结及讨论

叠后地震反演近年来被广泛应用于储层预测和理论研究中,但在实际中仍有很多问题需要注意.除了反演算法本身,地震资料、井资料、子波、层位等因素对最终的反演结果也有重要影响,对其正确合理的获取、应用和处理仍是叠后地震反演工作研究的重点.叠后地震反演的多解性主要来自地震数据之外的低频和高频信息,因此构建精确的高低频信息是反演工作成功的基础.如何减弱噪声的影响、摆脱理论假设的束缚、寻找更优的算法、减少多解性使反演结果能更加准确地反映地下的真实情况预计将是未来叠后地震反演的发展方向.

致谢 感谢审稿专家提出的修改意见和编辑部的大力支持!
参考文献
[] Becquey M, Lavergne M, Willm C .1979. Acoustic impedance logs computed from seismic traces[J]. Geophysics, 44 (9) : 1485–1501. DOI:10.1190/1.1441020
[] Chen H, Peng Z M, Wang J, et al .2011. Spectral decomposition of seismic signal based on fractional Gabor transform and its application[J]. Chinese Journal Geophysics, 54 (3) : 867–873. DOI:10.3969/j.issn.0001-5733.2011.03.028
[] Chen H Z, Yin X Y, Zhang J Q, et al .2014. Seismic inversion for fracture rock physics parameters using azimuthally anisotropic elastic impedance[J]. Chinese Journal Geophysics, 57 (10) : 3431–3441. DOI:10.6038/cjg20141029
[] Cheng B J, Xu T J, Li S G .2012. Research and application of frequency dependent AVO analysis for gas recognition[J]. Chinese Journal Geophysics, 55 (2) : 608–613. DOI:10.6038/j.issn.0001-5733.2012.02.023
[] Cui J C, Wu M, Li W X .2007. Nonlinear high-resolution impedance inversion method and its application[J]. Geophysical prospecting for petroleum, 46 (1) : 1–11.
[] Clarke G K C .1968. Time-varying deeonvolution filters[J]. Geophysics, 33 (6) : 936–944. DOI:10.1190/1.1439987
[] Dong Q, Lu S F, Zhang X J, et al. 2013. The application of construction technology for well logs to BMM area, 28(5):2579-2586, doi:10.6038/pg20130537.
[] Hass A, Dubrule O .1994. Geostatistical inversion-a sequential method of stochastic reservoir modelling constrained by seismic data[J]. First Break, 12 : 561–569.
[] He Y, Liu H S, Mao C L, et al .2008. Study on application of multi-curve sonic reconstructive technique in prediction of reservoir[J]. Oil Geophysical Prospecting, 43 (5) : 549–556. DOI:10.13810/j.cnki.issn.1000-7210.2008.05.015
[] Huang W C, Yang C C, Fan T Y, et al .2007. The application of petrophysical analysis in the reservoir prediction[J]. Progress in Geophysis, 22 (6) : 1791–1795.
[] Jiang Y, Xu L H, Zang X L, et al .2013. Prestack geostatistical inversion method and its application on the reservoir prediction of Changyuan oil field[J]. Progress in Geophysics, 28 (5) : 2579–2586. DOI:10.6038/pg20130537
[] Levin S A .1989. Surfaeoeonsistent deeonvolution[J]. Geophysies, 54 (8) : 1123–1133.
[] Li G F, Wang Y C, Xiong J L, et al .2010. Experimental analysis on seismic impedance inver-Sion[J]. Oil Geophysical Prospecting, 45 (6) : 868–872. DOI:10.13810/j.cnki.issn.1000-7210.2010.06.017
[] Li H B .1996. Generalized multi-traces linear inversion with the removal of noise[J]. Geophysical Prospecting For Petroleum, 35 (4) : 11–17.
[] Li Q Z .1989. A sand bed interpolation technique in line with the basic rule of thin layer comparision[J]. Geophysical Prospecting for Petroleum, 28 (1) : 12–21.
[] Li Q Z .1993. The Way To Obtain a Better Resolution in seismic Prospecting[M]. Beijing: Petroleum Industry Press .
[] Li Q Z .1998. On strategy of seismic restricted inversion[J]. Oil Geophysical Prospecting, 33 (4) : 423–438.
[] Li Z D, Li J, Li Y .2013. Structure modeling-constrained seismic inversion technology-a case study from Nanba area in Saertu oilfield, Daqing[J]. Oil & Gas Geology, 34 (3) : 401–405. DOI:10.11743/ogg20130318
[] Lindseth R O .1978. Synthetic sonic Logs-A process for stratigraphic interpretation[J]. Geophysics, 44 (1) : 3–26.
[] Lindseth R O .1979. Synthetic sonic logs:a process for stratigraphic interpretation[J]. Geophysics, 44 (1) : 3–26. DOI:10.1190/1.1440922
[] Lin X Z, Yang H Z, Tang L .1998. Multi-trace inversion without borehole data[J]. Oil Geophysical Prospecting, 33 (4) : 448–452.
[] Ling Y, Xi X Y, Sun D S .2008. Analysis on affecting factors of post-stack inversion and seismic attribute interpretation of thin reservoir[J]. Geophysical Prospecting For Petroleum, 47 (6) : 531–540.
[] Liu C Z, Ran J B, Yang Y S, et al .2002. Application of sparse pulse and inversion based on model in rolling exploration and development of sand body reservoir in Es4member, Wangjiagang area[J]. Oil Geophysical Prospecting, 34 (3) : 332–340. DOI:10.13810/j.cnki.issn.1000-7210.2002.04.019
[] Long G Q, Deng H W, Li Y J, et al .2008. Application of the log reconstruction technology to reservoir characterization based on forward modeling[J]. Petroleum Geology and Recovery Efficiency, 15 (1) : 66–69.
[] Lu W K, Li Y D, M Y G .1996. Seismic Inversion Using Erro-back-Propagation Neural Network[J]. Chinese Journal Geophysics, 39 (S1) : 292–301.
[] Ma J F .1999. Multiple roots in well-restricted impedance inversion[J]. Oil & Gas Geology, 20 (1) : 322–325.
[] Ma J F .2003. Forward modeling and inversion method of generalized elastic impedance in seismic exploration[J]. China Journal Geophysics, 46 (1) : 118–124.
[] Ma J F, Wang X J, Jia C H, et al .2000. Study of constraint methodology in constrained impedance inversion[J]. Geophysical Prospecting for Petroleum, 39 (2) : 52–63.
[] Ma J F, Wang X J, Jia C H, et al .2000. Study of constraint methodology in constrained impedance inversion[J]. Geophysical Prospecting for Petroleum, 39 (2) : 52–63.
[] Ma J F, Wang X J, XieY G, et al .2000. Experience and skill of constructing low-frequency components in impedance inversion[J]. Geophysical Prospecting for Petroleum, 39 (1) : 27–34.
[] Ma J F, Wang X J, Zhao S L, et al .1999. Choice of initial wave impedance value and the multisolution of seismic trace inversion[J]. Oil Geophysical Prospecting, 34 (3) : 332–340. DOI:10.13810/j.cnki.issn.1000-7210.1999.03.012
[] Ma J F, Xu S H, Wang G S, et al .2002. Problems and progress in seismic trace inversion[J]. Oil & Gas Geology, 23 (4) : 7–10.
[] Meng X J, Jin X L, Niu W M, et al .2004. Research of 3-D complex constrained model in seismic inversion[J]. Journal of the University of Petroleum, China, 28 (6) : 21–26.
[] Pan Y J, Li D, Yang K .2011. A comparison between the requirements of multi-well constrained conditions in deterministic inversion and stochastic inversion[J]. Geophysical Prospecting For Petroleum, 50 (4) : 345–349. DOI:10.3969/j.issn.1000-1441.2011.04.005
[] Qi Y F, Yin C C, Wang R, et al .2015. Multi-transient EM full-time forward modeling and inversion of m-sequences[J]. Chinese Journal Geophysics, 58 (7) : 2566–2577. DOI:10.6038/cjg20150731
[] Qiang M, Miao Q M, Zhu W M, et al .2015. Method of quantitative interpretation of seismic elastic multi-parameter cross plot on effective reservoir prediction in SG gas field[J]. Progress in Geophysics, 30 (1) : 293–299. DOI:10.6038/pg20150142
[] Robinson E A .1957. Predietive deeomposition of seismic traees[J]. Geophysics, 22 (4) : 767–778. DOI:10.1190/1.1438415
[] Robinson E A .1967. Predietive deeomposition of time series with applieation to seismic exploration[J]. Geophysics, 32 (3) : 418–484. DOI:10.1190/1.1439873
[] Sun S M, Peng S M .2007. Inversion of geotatistics based on simulated annealing algorithm[J]. Oil Geophysical Prospecting, 42 (1) : 38–43.
[] Wang H, Chen Z Q, Xiao M H .2009. Research of Condition on Seismic Elastic Impendance[J]. Journal of Yangtze University (Nat Sci Edit), 6 (3) : 164–166.
[] Wang H, Luo B, Li T .2010. Comparative analysis of elastic impedance[J]. Lithologic Reservoirs, 22 (3) : 110–113.
[] Wang S Z, Zhou H K .2009. Application of logging curve rebuilding technology in reservoir prediction at east slope of Chengdao Oilfield[J]. Progress in Exploration Geophysics, 32 (5) : 376–379.
[] Wei C, Li X F, Zhang M G .2008. The geophysical inverse method based on Quantum Monte Carlo[J]. Chinese Journal Geophysics, 51 (5) : 1494–1502.
[] Xi X Y, Ling Y, Sun D S, et al .2009. A Research on the factors to influence post-stack seismic inversion during oilfield production[J]. China Offshore Oil and Gas, 21 (5) : 308–312.
[] Xu S H, Ma J F .2003. Wavelet extraction method of wave impedance inverse with multi-well restrictions[J]. Progress in Geophysics, 18 (4) : 623–627.
[] Xu Z L .2002. Application of the curve-reconstruction technology to the lateral prediction of reservoirs[J]. Petroleum Geology & Experiment, 24 (4) : 377–380.
[] Yang S H, Huang Y S, Peng W X, et al .2006. Application of sonic reconstruction technique in inversion of reservoir[J]. Oil Geophysical Prospecting, 1 (2) : 171–176.
[] Yao F C, Gan L D .2000. Application and restriction of seismic inversion[J]. Petroleum Exploration and Development, 27 (2) : 63–56.
[] Yang L Q .2003. A review of well log constrained seismic inversion[J]. Progress in Geophysics, 18 (3) : 530–534.
[] Zhang J, Yang Q L, Wang T Q .2010. Probing for seismic inversion in depth domain[J]. Oil Geophysical Prospecting, 45 (S1) : 114–116.
[] Zhang Y G .2002. The present and future of wave impedance inversion technique[J]. Geophysical Prospecting for Petroleum, 41 (4) : 385–389.
[] Zhang Y J, Xuan Y H, Qin C G .2010. Application of petrophysical analysis for gas-bearing prediction of Panyu Gasfield[J]. Geophysical Prospecting for Petroleum, 49 (4) : 407–410. DOI:10.3969/j.issn,1000-1441.2010.04.012
[] Zhou S S, Yi W, Hao Z B, et al .2012. Experiment research and application of fluid sensitive attributes based on the pre-stack inversion[J]. Chinese Journal Geophysics (6) : 1985–1992. DOI:10.6038/j.issn.0001-5733.2012.06.020
[] Zong Z Y, Yin X Y, Zhang F, et al .2012. Reflection coefficient equation and pre-stack seismic inversion with Young's modulus and Poisson ratio[J]. Chinese Journal Geophysics, 55 (11) : 3786–3794. DOI:10.6038/j.issn.0001-5733.2012.11.025
[] 陈红, 彭真明, 王峻, 等.2011. 地震信号分数阶Gabor变换谱分解方法及应用[J]. 地球物理学报, 54 (3) : 867–873. DOI:10.3969/j.issn.0001-5733.2011.03.028
[] 陈怀震, 印兴耀, 张金强, 等.2014. 基于方位各向异性弹性阻抗的裂缝岩石物理参数反演方法研究[J]. 地球物理学报, 57 (10) : 3431–3441. DOI:10.6038/cjg20141029
[] 程冰洁, 徐天吉, 李曙光.2012. 频变AVO含气性识别技术研究与应用[J]. 地球物理学报, 55 (2) : 608–613. DOI:10.6038/j.issn.0001-5733.2012.02.023
[] 崔炯成, 吴媚, 李维新.2007. 高分辨率非线性地震波阻抗反演方法和应用[J]. 石油物探, 46 (1) : 1–11.
[] 董奇, 卢双舫, 张学娟, 等.2013. 地质统计学反演参数选取及反演结果可靠性分析[J]. 物探与化探, 37 (2) : 345–349. DOI:10.11720/j.issn.1000-8918.2013.2.29
[] 贺懿, 刘怀山, 毛传龙, 等.2008. 多曲线声波重构技术在储层预测中的应用研究[J]. 石油地球物理勘探, 43 (5) : 549–556. DOI:10.13810/j.cnki.issn.1000-7210.2008.05.015
[] 黄伟传, 杨长春, 范桃园, 等.2007. 岩石物理分析技术在储层预测中的应用[J]. 地球物理学进展, 22 (6) : 1791–1795.
[] 姜岩, 徐立恒, 张秀丽, 等.2013. 叠前地质统计学反演方法在长垣油田储层预测中的应用[J]. 地球物理学进展, 28 (5) : 2579–2586. DOI:10.6038/pg20130537
[] 李国发, 王艳仓, 熊金良, 等.2010. 地震波阻抗反演实验分析[J]. 石油地球物理勘探, 45 (6) : 868–872. DOI:10.13810/j.cnki.issn.1000-7210.2010.06.017
[] 李宏兵.1996. 具有剔除噪声功能的多道广义线性反演[J]. 石油物探, 35 (4) : 11–17.
[] 李庆忠.1989. 符合小层对比原则的砂层内插技术[J]. 石油物探, 28 (1) : 12–21.
[] 李庆忠.1993. 走向精确勘探的道路[M]. 北京: 石油工业出版社 .
[] 李庆忠.1998. 论地震约束反演的策略[J]. 石油地球物理勘探, 33 (4) : 423–438.
[] 李占东, 李吉, 李阳, 等.2013. 构造建模约束地震反演技术--以大庆萨尔图油田南八区为例[J]. 石油与天然气地质, 34 (3) : 401–405. DOI:10.11743/ogg20130318
[] 林小竹, 杨慧珠, 汤磊.1998. 无井多道反演[J]. 石油地球物理勘探, 33 (4) : 448–452.
[] 凌云, 孙德胜, 林吉祥, 等.2008. 薄储层叠后反演影响因素分析与地震属性解释研究[J]. 石油物探, 47 (6) : 531–540.
[] 刘成斋, 冉建斌, 杨永生, 等.2002. 稀疏脉冲和基于模型反演在王家岗沙四段砂岩油藏精细勘探中的应用[J]. 石油地球物理勘探, 37 (4) : 401–406. DOI:10.13810/j.cnki.issn.1000-7210.2002.04.019
[] 龙国清, 邓宏文, 李玉君, 等.2008. 基于模型正演的测井曲线重构技术在储层表征中的应用[J]. 油气地质与采收率, 15 (1) : 66–69.
[] 陆文凯, 李衍达, 牟永光.1996. 误差反传播神经网络法地震反演[J]. 地球物理学报, 39 (S1) : 292–300.
[] 马劲风.1999. 测井资料约束的波阻抗反演中的多解性问题[J]. 石油与天然气地质, 20 (1) : 322–325.
[] 马劲风.2003. 地震勘探中广义弹性阻抗的正反演[J]. 地球物理学报, 46 (1) : 118–124.
[] 马劲风, 王学军, 贾春环, 等.2000. 波阻抗约束反演中的约束方法研究[J]. 石油物探, 39 (2) : 52–63.
[] 马劲风, 王学军, 谢言光, 等.2000. 波阻抗反演中低频分量构建的经验与技巧[J]. 石油物探, 39 (1) : 27–34.
[] 马劲风, 王学军, 赵圣亮, 等.1999. 波阻抗初值选择与地震道反演的多解性[J]. 石油地球物理勘探, 34 (3) : 332–340. DOI:10.13810/j.cnki.issn.1000-7210.1999.03.012
[] 马劲风, 许升辉, 王桂水, 等.2002. 地震道反演面临的问题与进展[J]. 石油与天然气地质, 23 (4) : 7–10.
[] 孟宪军, 金翔龙, 钮学民, 等.2004. 地震反演中的三维复杂约束模型[J]. 石油大学学报(自然科学版), 28 (6) : 21–26.
[] 潘昱洁, 李大卫, 杨锴.2011. 确定性反演和随机反演对井约束条件的需求分析[J]. 石油物探, 50 (4) : 345–349. DOI:10.3969/j.issn.1000-1441.2011.04.005
[] 强敏, 苗庆梅, 朱望明, 等.2015. SG地区叠前地震弹性参数交会定量解释技术预测有效储层的方法[J]. 地球物理学进展, 30 (1) : 293–299. DOI:10.6038/pg20150142
[] 齐彦福, 殷长春, 王若, 等.2015. 多通道瞬变电磁m序列全时正演模拟与反演[J]. 地球物理学报, 58 (7) : 2566–2577. DOI:10.6038/cjg20150731
[] 孙思敏, 彭仕宓.2007. 基于模拟退火算法的地质统计学反演方法[J]. 石油地球物理勘探, 42 (1) : 38–43.
[] 王浩, 陈志强, 肖梦华.2009. 弹性阻抗对砂岩的敏感性研究[J]. 长江大学学报(自然科学版), 6 (3) : 164–166.
[] 王浩, 罗兵, 李霆.2010. 弹性阻抗对比分析[J]. 岩性油气藏, 22 (3) : 110–113.
[] 王绍忠, 周红科.2009. 测井曲线重构技术在埕岛油田东斜坡储层预测中的应用[J]. 勘探地球物理进展, 32 (5) : 376–379.
[] 魏超, 李小凡, 张美根.2008. 基于量子蒙特卡罗的地球物理反演方法[J]. 地球物理学报, 51 (5) : 1494–1502.
[] 晓宇, 凌云, 孙德胜, 等.2009. 油田开发阶段叠后地震反演影响因素研究[J]. 中国海上油气, 21 (5) : 308–312.
[] 许升辉, 马劲风.2003. 波阻抗反演中多井条件下的子波提取方法[J]. 地球物理学进展, 18 (4) : 623–627.
[] 许正龙.2002. 曲线重构技术在储层横向预测工作中的应用[J]. 石油实验地质, 24 (4) : 377–380.
[] 杨少虎, 黄玉生, 彭文绪, 等.2006. 声波重构技术在储层反演中的应用[J]. 石油地球物理勘探, 1 (2) : 171–176.
[] 姚逢昌, 甘利灯.2000. 地震反演的应用与限制[J]. 石油勘探与开发, 27 (2) : 63–56.
[] 杨利强.2003. 测井约束地震反演综述[J]. 地球物理学进展, 18 (3) : 530–534.
[] 张静, 杨勤林, 王天琦.2010. 深度域地震反演方法探索[J]. 石油地球物理勘探, 45 (S1) : 114–116.
[] 张永刚.2002. 地震波阻抗反演技术的现状和发展[J]. 石油物探, 41 (4) : 385–389.
[] 张永江, 轩义华, 秦成岗.2010. 岩石物理分析在番禺天然气区含气性研究中的应用[J]. 石油物探, 49 (4) : 407–410. DOI:10.3969/j.issn,1000-1441.2010.04.012
[] 周水生, 宜伟, 郝召兵, 等.2012. 基于叠前反演的流体敏感属性实验研究及应用[J]. 地球物理学报 (6) : 1985–1992. DOI:10.6038/j.issn.0001-5733.2012.06.020
[] 宗兆云, 印兴耀, 张峰, 等.2012. 杨氏模量和泊松比反射系数近似方程及叠前地震反演[J]. 地球物理学报, 55 (11) : 3786–3794. DOI:10.6038/j.issn.0001-5733.2012.11.025