井间地震层析成像的现状与进展

1、2001年9月地球物理学进展第16卷第3期井间地震层析成像的现状与进展裴正林（石油大学（北京）物探重点实验室，北京，100083）摘要：综述了井间地震层析成像研究的现状，给出了小波变换域井间地震层析成像方法的最新进展，并对井间地震层析成像研究给予展望.关键词：井间地震层析成像；小波多尺度；研究进展中图分类号：P315.3+1文献标识码：A文章编号：（2001）1004-290303-0091-071井间地震层析成像的研究现状井间地震层析成像也称为井间地震CT技术，它能够提供被探测地质体的构造和岩性分布的高分率图像.井间地震CT技术是从医学CT技术发展起来的，其数学基础是Radon变换.井间地震

2、CT的研究基本始于20世纪70年代初，80年代处于对大量模型数据和少量实际数据的成像研究阶段，井间地震CT进入实用化阶段，并取得不少可喜成果，90年代以来，同时，也逐渐意识到射线CT所固有的缺点，开始研究波动方程CT.从地震波的运动学和动力学特征出发，井间地震CT方法可分为两大类：一类是基于几何光学或射线方程的方法称之为射线CT；另一类是基于波动方程的方法称之为波形CT.当非均匀体的线性尺度大于地震波长时，射线CT是适用的；而当非均匀体的线性尺度与波长相近时，衍射和散射就起主导作用了，基于射线理论的成像方法就不再适用，这时候必须用波动方程CT方法.井间地震层析成像方法主要包括两部分：正演方法和

3、反演方法.井间地震层析成像的正演方法可分为两种；一是射线追踪方法；二是波场的数值模拟方法.射线理论和射线方法是研究地震波传播理论的重要方面之一.用射线理论可以研究地下复杂构造、横向不均匀介质中的地震波传播问题.经过射线追踪，计算地震波的走时、波前和射线路径.170年代以前的各种射线追踪方法一般适合于较为简单模型的射线追踪.由于实际的介质速度变化较大（速度差大于10%），因此，需要研究复杂结构模型的射线追踪方法.收稿日期：修订日期：2001-03-15；2001-06-15.基金来源：“九五”国家科技攻关项目资助（959130602）.作者简介：裴正林，（北京）地球探测与信息技术专业博士.高级工

4、程师，现在石油大学1962年生，2000年获中国地质大学（北京）从事博士后研究.主要研究方向：信号处理，小波变换、遗传算法及神经网络应用，层析成像理论方法和地震数据处理、偏移方法等方面研究.E-mail：zhenglinpei.92地球物理学进展2001年近年来，这方面的研究取得新的进展.目前应用效果好的弯曲射线追踪方法有两类：一类是基于程函方程的有限差分方法；另一类是基于Huygens原理的最短路径射线追踪法.2，3基于程函方程的有限差分方法有有限差分方法的E-O格式，它能以较高的精度计4，5，6算较复杂模型的初至波走时.线性插值射线追踪法（LTI）是程函方程有限差分射线追7，8踪的改进.这

5、类方法是全局算法，能模拟较复杂介质，编程简单、计算速度快、精度较高.9，10基于Huggens原理的射线追踪有两种全局算法：一是基于网络理论的最短路径算法；二11，12是基于动力学的波阵面算法而且是一次计算即可得.这两种算法都能模拟复杂介质，到全部网格节点上的最短路径和最小走时，计算效率高.井间地震波场的正演数值模拟是求解变系数偏微分方程的过程.目前，解偏微分方程的13数值方法中较常用的有：有限差分法，有限元法，边界元法，伪谱法，积分法等.其中有限差分法比较容易程序实现，并且具有较快计算速度和较高计算精度，很适合在波动方程层析成像中使用.目前，井间地震层析成像2-D正演方法基本趋于成熟.井间地

6、震层析成像是非线性反问题.井间地震射线CT是通过逐次线性化反演来实现的，即可归结为求解一个大型的、稀疏的、常常是病态的线性方程组.为了在这种情况下求出满意的解来，必须从数学上把病态问题转化为良态问题求解，并把更多的先验信息引入到反演过程中.近十年来，许多研究者不断研究发展新的方法，改善了在数据不完全，有噪声和初14始模型不准情况下的成像效果.射线CT从直射线层析成像方法发展到弯曲射线层析成10，17，15，1610，15，18像方法；反演万法由最小二乘法发展到各种约束条件下的阻尼最小二乘法19以及统计法如最大熵法；观测参数数据也由单一的走时向多种参数数据（P波、S波走时、20以及多震相成像方法

7、如直振幅、信号谱等）发展.井间地震P波、S波走时层析成像方法21射线的反射十透射井间层析成像方法及弯曲射线的反射十透射井间层析成像（CCRTT）22方法；井间地震和逆VSP联合层析成像方法，一定程度上克服了井间地震透射层析成像23的多解性及有效成像区域小的缺陷.衍射CT是波动方程CT的一次近似（假定弱散射，.有两种方法：V<15%的小扰动）而Rytov近似则利用波场的相位.与射线Born法和Rytov法.Born近似是利用波场的振幅，24，25，26衍射CT提高了分辨率，可以减少由于有限观测角所造成的假像CT相比，.在一定条件约束下，在高频假定了，Born近似和Rytov近似的效果是相同

8、的；Rytov近似层析与走时层2728析是等价的或密度和压缩系数的像.声波衍射CT方法可以获得纵波速度的像，.弹性29，30，31波衍射成像方法可以获得纵横渡速度或密度和拉梅常数的像.目前衍射CT大多只32，33，34，35用于均匀的背景介质，而变背景的比较少.散射CT方法是广义Radon变换的一级近似逆.它可用于波还不均匀的背景介质.散射但在不能给出准确CT方法在己知准确的背景速度分布时可以给出分辨率很高的成像结果，36，3837背景速度场时效果很差可以采用级联算法.为了得到准确的背景速度分布，.实际应用时，衍射CT和散射CT方法都存在一些有待解决的问题，如散射波的提取、对波形产生影响的各种

9、干扰因素的消除（震源函数、介质的吸收、震源辐射、接收器的耦39合），而且这些成像方法只有当散射体与背景的速度差较小时才有效，因而使其实际应用受到限制.而波形CT方法是利用全波场信息，能够正确地反演地下的真实介质模型，因此，具有第16卷第3期裴正林：井间地震层析成像的现状与进展93潜在的科学价值.近年来，波形CT方法成为研究的热点，并有所发展.波形CT又可分声波方程CT和弹性波（全波）方程CT.声波方程是弹性波方程的近似.对2D声波方程速度反40，41，42，43演研究比较多且给出Vp，.虽然弹性波方程CT给出的图像略好于声波方程的，Vs40图像，但声波方程的计算时间约为弹性波方程计算时间的1/

10、8.由于波形反演的目标函数中存在大量的局部极小，因此，波动CT存在着收敛速度慢，对初始模型依赖性强以及易于陷入局部极小的缺陷.44，45有一种波动方程走时反演方法.该方法是利用地震波的走时差来代替波形反演中的波形拟合差进行波形反演的，提高了抗噪能力和速度成像的分辨率（相对于射线走时成像而言）加快了收敛速度.但这类方法仅利用了波形的走时信息，从某种意义上说，它减.另外，少了波形反演的局部极小个数，但这是以降低波形反演的分辨率为代价.射线CT获得是速度分布的光滑估计即速度模型的低频成分，提供低分辨率图像，而波动方程CT是估计速度模型的不均匀性和不连续性即高频成分.因此，二者联合使用则很适40合于高

11、分辨率的图像重建.井间地震层析成像问题的非线性是与介质的非均匀性联系在一起的.对于小扰动非均匀介质，常规的线性化反演方法可以得到较好的解.而对于大扰动的非均匀介质的成像问题，采用常规的线性化反演方法，反演结果不好，这是由于在此种情况下，线性化近似已经不46再有效.可以采用逐次线性化反演方法，其结果尚可满意.井间地震层析成像是非线性反问题，具有非线性性、多解性与不适定性、高维数（为达到分辨率要求，模型网格点数量很大）；另外，还存在弯曲射线、大扰动、有限频带以及有限观测角（观测数据不完备）等问题，因此，井间地震层析成像是复杂非线性反问题.井间地震层析成像反演方法可分为两类：基于算子的线性化反演方法

12、和基于模型的完全非线性反演方法.线性化或拟线性化的反演方法又称为确定性反演方法如LSQR法、梯度法、牛顿法等；完全非线性反演方法又称为随机反演方法如MonteCarlo法、模拟退火算法（SA）等.确定性反演方法一般是局部寻优的，而随机性反演方法多为全局寻优.线性化CT反演方法如正则化方法，散射法等，通过忽略高波数成分即局部极小，建立了适定的算法，换Born和Rytov近似，言之，线性化反演方法是通过忽略参数的部分非线性成分，改善了反问题的适定性，反演出有意义的解.因此，线性化CT反演方法只有当初始模型在目标函数的全局极小的邻域内（在真实模型附近）时，近似线性化的关系才成立，即只有当引人很好的背

13、景速度场先验信息约束时，线性化CT反演结果才能接近全局最优解，否则易陷入局部极小.因此，现行的井间地震层析成像方法缺点是，对初始模型的依赖性强，即图像的质量受初始模型的控制，仅适合于低对比度（速度差小）异常体成像，图像分辨率偏低以及不可缺少的先验知识等.很显然，对于复杂模型和大扰动非均匀体成像问题，现行的井间地震层析成像方法无法得到好的图像.多重网格法是目前多极值目标函数反演的有效方法之一.可以有效地改善常规一次性反演方法的性能.已有不少研究者进行了这方面的研究.地震反射波多重网格法2-D反47，48，49演，井间地震走时波形联合反演方法是解决波形反演中局部极小问题另一种方50法.尽管多尺度方

14、法一定程度上提高了线性化层析成像图像的分辨率和图像质量，但仍未94地球物理学进展2001年完全克服线性化反演方法的不足，自然无法彻底地解决非线性问题.因此，地球物理学者将随机概念引入非线性问题反演，而且取得了不小的突破，出现了许多现实的反演方法，如模拟退火法（SA）、遗传算法（GA）等.随机反演方法对目标函数不作任何数学MonteCarlo法、要求（如连续性、可导性），适合于复杂的多极值目标函数的全局优化.将SA法和GA法应用51，5552，53，54于1-D地震波形反演.应用SA法进行地震反射成像研究.应用GA法研究了5657，582-D地壳速度结构反射走时反演.采用SA法进行简单模型走时层

15、析成像研究.显见，二维井间地震层析成像的随机反演算法处于研究阶段，其应用实例也很少.对于复杂非均匀介质，现行的井间地震层析成像方法实际应用时往往图象分辨率不高、效果也不理想.22.1井间地震层析成像的进展小波多尺度地震层析成像小波分析是当前研究非线性科学的主要工具之一.小波变换也是一种多分辨变换，它能够将信号的低频成分映射到由大尺度小波基张成的粗采样子空间，把高频成分映射到由小尺度小波基张成的较细采样子空间，因此，小波变换为多重网格法提供了一种自然框架.小59，60波变换的这一特性已经应用于1-D地震波形反演.对于井间地震层析成像非线性反问61题，一种方法是通过小波分解压缩速度模型，来提高成像

16、计算速度；另外一种方法是在小62波域对线性化方程进行多尺度分解的小波多尺度地震层析成像方法.该方法将地震层析成像反问题分解为不同尺度上的一系列反问题.大尺度的反问题的目标函数局部极小少而且互相分得很开，因而用常规的反演方法可以较容易地收敛到全局极小的邻域内；将大尺度的解当作次一级尺度反问题的初始模型，再进行反演.如此类推，直到求出对应于最小尺度的原反问题的最优解.小波变换多尺度方法有效地改善了常规线性化反演方法依赖初始数模型、易陷入局部极小的不足，是复杂介质和大扰动非均匀介质成像的有效方法，另外，小波多尺度方法还能够提供不同分辨率的反演图像，有益于图像的解释.2.2井间地震多尺度遗传进化层析成

17、像针对井间地震层析成像反问题所具有的特点：模型参数维数高，且参数在空间分布上具有一定的规律性，模型可行解空间巨大等，把简单遗传进化算法（EvolutionaryProgramming，简6362称EP）应用于井间地震层析成像反演时应进行改进.（1）二维实数编码地震层析成像问题类似于二维图像的恢复问题.一个个体就是一幅图像.此时，若采用一维染色体编码则很不方便，尤其是变异操作很不直观.因此，采用二维实数编码.每个二维染色体就对应于一个二维速度模型，每个实数编码的基因就表示对应的速度网格节点上的速度.（2）模型扰动方法在遗传进化算法中，采用模型局部扰动方式作为对父代染色体的变异.二维染色体的变异不

18、同于一维染色体的点位变异，而是采用基因块变异或模型速度块扰动.模型速度块大小的选取的方法有两种：随机大小速度块和固定大小速度块.随机速度块，其大小可以从一个像元到整个模型变化.模型的速度块的速度扰动范围是自适应可变的，最大范围从模型第16卷第3期裴正林：井间地震层析成像的现状与进展95速度的下限到上限.这样一来，染色体（速度模型）的变异自动实现约束，并约束在速度参数值的选择范围内，以提高计算可靠性及计算速度.（3）多尺度逼近遗传进化算法我们将小波多尺度分析思想也引入到随机最优化过程中，将多尺度反演和进化算法反演的结合起来，称为多尺度逼近遗传进化算法（MultiscaleEvolutionary

19、ProgrammingAlgorithm，简称MEP）一是对原反问题分解为尺度由小到大的反问题序列.随.该方法的有三个特点：着尺度的增大，对应反问题的解空间范围缩小，这样以来，加速了反演搜索速度，减少了该尺度反问题的求解计算成本.二是对于大尺度反问题，可以通过小尺度来细化.三是与常规遗传进化算法所采用的固定搜索范围不同，多尺度反演过程中，通过在大尺度上解的邻域附近来扰动，既保证了有效基因不被丢失又使搜索范围是逐尺度缩小，速度基因值的分辨率提高，这样以来，逐尺度提高了解的精度和分辨率.3井间地震层析成像的发展方向由上可以看出，如何提高图像分辨率和成像效果，仍是今后提高井间地震析成像技术实用化程度

20、所研究的重点.引起图像分辨率偏低、成像效果不好的主要原因是：（1）有限观测角；（2）介质的非均匀性和各向异性；（3）观测数据含有误差；（4）成像方法本身的缺陷（如线性化、正则化和弱散射近似等）.因此，井间地震CT仍需从反演理论及方法、反演算法，野外施工与数据采集等方面进行研究、改进，进而很好地解决上述几方面的问题，才能提高地震CT的地质效果.目前，井间地震层析成像技术需要深入研究的主要内容有：（1）含噪声不完全投影数据的3-D地震层析成像方法；（2）复杂介质井间地震弹性波波场快速数值模拟方法；（3）双相介质井间地震层析成像方法；（4）各向异性介质井间地震层析成像方法；（5）3-D多波井间地震层

21、析成像方法；（6）矢量（三分量或多分量）井间地震层析成像方法；（7）四维（时间推移）井间地震层析成像方法.参考文献1CervenyV.etal.著，刘福田，等译.地震学中的射线方法M.地质出版社，1986.2VanTrierJ，J（3）：SymesMM.Upwindfinite-differencecalculationtraveltime.Geophysics，1991，56812821.3周洪波，张关泉.复杂构造区域的初至波走时计算J（4）：.地球物理学报，1994，37515520.4AsakawaE，J（l）：KanawaT.Seismicraytracingusinglineartr

22、aveltimeInterPolation.GeophysicalProspecting，1993，4199111.5赵改善，等.基于旅行时线性插值的地震射线追踪算法J（2）：.石油物探，1998，371424.96地球物理学进展2001年6张霖斌，等.有限差分法射线追踪J（6）：.石油地球物理勘探，1993，28673677.7Qinetal，Solutionoftheeikonalequationbyafinitedifferencemethod，60thSEGAnnualMeeting，1990.9Moser，T.J.，J（1）：Shortestpathcalculationofseis

23、micrays.Geophysics，1991，565967.8VidaleJ.Finite-differencecalculationoftraveltlmeinthreedimensionJ（5）：.Geophysics，1990，55521526.9MoserTJ.ShortestpathcalculationofseismicraysJ（1）：.Geophysics，1991，565967.10周兵，赵明阶.最小走时射线追踪层析方法J（2）：.物探化探计算技术，1992，14124130.11RobertL.Coultrip，Highaccuracywavefronttracingtr

24、aveltlmecalculationJ（2）：.Geophysics，1993，58284292.12刘洪，李幼铭.计算最小走时和射线路径的界面网全局方法J（6）：.地球物理学报，1995，38823831.13范祯祥，等.地震波数值模拟与偏移成像M.河南科学技术出版社，1994.14GuustNolet.冯锐，等译.地震层析技术M.地质出版社，1991.15周兵，等.一种新的地震射线层析成像计算方法J（1）：.石油物探，1994，334554.16吴律，著.层析基础及其在井间地震中的应用M.石油工业出版社，1997.17CarrionP.Dualtomographyforimagingco

25、mplexstructureJ（9）：.Geophysics，1991，5613951404.18张文生，何樵登.约束走时层析成像J（l）：.石油地球物理勘探，1997，326874.19殷军，冯锐.井间层析成像的最大滴方法J（2）：.地球物理学报，1992，35234241.20王成礼，孟凡林.井间地震纵横波走时层析成像处理方法J（4）：.石油地球物理勘探，1996，31509514.21Calmanetal.Reflection+transmissionforcrosswelltomographicimagingC.59thSEGAnnualMeeting，1989.22Wuluetal.

26、SeismiccrossholecurvedrayreflectionplustransmissiontomographyC.64thSEGAnnulMeeting，1994.23张树林，朱介寿.井间地震和逆VSP联合层折成像J（5）：.石油地球物理勘探，1993，28577583.24Wu，JR.S.andToksoz，M.N.Diffractiontomographyandmultisourceholographyappliedtoseismicimaging.Geophysics，1987，（l）：521125.25LoTW，ToksozMN，JXuSH，WuRS.Ultrasonicl

27、aboratorytestsofgeophysicaltomographicreconstruction.Geophysics，（7）：1988，53947956.26PrattRG，WorthingtonMH.Theapplicationofdiffractiontomographytocross-holeseismicdataJ（10）：.Geophysics，1988，5312841294.27吴律，朱光明.Born近似层析与Rytov近似层析对比研究J（1）：.石油地球物理勘探，1993，28116.28彭成斌，陈隅.衍射CT技术和多源全息成像技术的比较研究J（2）：.地球物理学报，1

28、990，33154162.29刁顺，等.弹性波衍射CTJ（3）：.石油物探，1994，334753.30孙豪志，范祯祥.弹性被介质中纵横波速度反演J（6）：.石油地球物理勘探，1994，29678684.31陈湛文，等.弹性介质中密度与拉梅常数、J（2）：.地球物理学报，1995，38234µ的衍射层析成像方法研究242.32PaiDM.CrossholeseismicusingverticaleigenstatesJ（7）：.Geophysics，1990，55815820.33DickensTA.Diffractiontomographyforcrosswellimagingof

29、nearlylayeredmediaJ（5）：.Geophysics，1994，59694706.34黄朕捷，吴如山.垂直非均匀背景多频背向散射层析成像J（1）：.地球物理学报，1994，3787100.35井西利，杨长春.一种非均匀背景的散射层析成像方法研究J（2）：.石油物探，1997，36714.36TarantoraA.InversionofseismicreflectiondataintheacousticapproximationJ（8）：.Geophysics，1984，4912591266.37黄联捷，杨文采.声波方程送散射反演的近似方法J（3）：.地球物理学报，1991，34

30、626634.38社剑渊，杨文采.跨孔地震层析成像的级联方法J（6）：.地球物理学报，1991，34771780.39TuraMAC，GreavesRJ，BeydounWb.Crosswellseismicreflection/diffractiontomography：Areservoircharacterizationapplica-J（3）：tion.Geophysics，1994，59351361.40PrattRG，GoultyNR.Combiningwave-equationimagingwithtraveltimetomographytoformhigh-resolutionim

31、agesfromcrossholeJ（2）：data.Geophysics，1991，56208224.41ReiterDT，J（3）：RodiW.Nonlinearwaveformtomographyappliedtocrossholeseismicdata.Geophysics，1996，61902913.42王守东，刘家琦.二维声波方程速度反演的一种方法J（6）：.地球物理学报，1995，38833839.43王大川，严忠琼，层析成像中的波形反演J（4）：.物探化探计算技术，1996，18289303.第16卷第3期裴正林：井间地震层析成像的现状与进展9744LuoY，J（5）：Schu

32、sterGT.Wave-equationtraveltimeinversion.Geophysics，1991，56645653.45樊计昌，许云.最大摘走时波形层析成像J（2）：.石油物探，1995，344652.46牛彦良，杨文采.跨孔地震CT中的逐次线性化方法J（3）：.地球物理学报，1995，38378386.47Bunks，J（5），C，Saleck，PM，Zaleski，S，Chavent，G.MultlscaleseismicwaveformInversion.Geophysics，1995，6014571473.48陈小宏，牟永光.地震波动方程反演的多重网格方法C地震出版社，.

33、中国地球物理学会年刊，1995.49钱建良，刘家琦.波动方程反问题的多尺度反演方法C地震出版社，.中国地球物理学会年刊，1995.50ZhouC，CiW，LuoY，schusteLG，HassanzadehS.AcousticwaveequationtrveltimeandwaveforminversionofcrossholeseismicJ（3）：dada.Geophysics，1995，60765774.51StoffaPL，JSenMK.Nonlinearmultlparameteroptimizationusinggeneticalgorithms：inversionofplane-

34、waveseismograms.（11）：1991，5617941810.Geophysics，52SenMK，JStoffaRL.Nonlinearone-dimentionalseismicwaveforminversionusingsimulatedannealing.Geophysics，1991，56（10）：16241638.53SenMK，JStoffaPL.Rapidsamplingofmodelspaceusinggeneticalgorithms：examplesfromseismicwaveforminversion.Geo-phys.J.Int.，1992，108：281292.54MosegaardK，JVestergaardRD.