地震波衰减文献综述.docx

1、地震波衰减文献综述文献综述1影响地层吸收衰减的主要因素地震波在地层中传播的过程中会存在能量衰减，这种衰减会受到许多因素的影响和制约。这些因素包括：频率、压力、温度、饱和度、应变振幅以及岩石的特性等。在研究地层吸收衰减特性的过程中，了解这些因素的影响作用对于衰减问题的研究是很有帮助的。1. 频率 频率与衰减的关系目前尚未有定论。室内研究表明 Q值与频率有关，而一些对实际地震数据的研究则表明了衰减与频率无关。现有的资料表明：对不同特性的岩石，频率的影响不同。对干燥岩石，衰减与频率无关（BirchandBancroft，1938；Born，1941；PeselnickandOuterbridge，1

2、961；PanditandSavage，1973；NurandWinkler，1980；Tittmannetal.1981）；对于部分饱和或完全饱和岩石，地震波以复杂的路线传播，由流体流动类型所决定，衰减通常与频率有关。Johnston等利用与地震勘探有关的孔隙流体的粘滞系数和标准线性粘滞性模型，计算出地震频段的衰减对频率的依赖关系；OConnelland Budiansky(1977)分析了饱和碎屑岩石的弹性特征，提出了与频率相关的衰减模型，他们指出在两种特性频率条件下产生的衰减最大；White(1975)计算出了在部分饱和流体岩石的弹性波衰减，推断出 P波的衰减和频率有关，而 SH波的衰减

3、和频率无关。 2.岩性 高速的岩石，吸收性弱，而低速的岩石，吸收性强。对于大多数地区，泥岩的平均吸收性比砂岩强，砂岩的吸收比页岩和灰岩的吸收强，砂岩含有油气时，其吸收性显著增强。总之，介质弹性越好，地震波在介质中传播的能量损耗 3压力 P波、S波在所有饱和岩石中，随压力的增加Q增大（衰减减小），在高压下则保持为一稳定值（Gadneretal.1964；Gordon and Davis，1968；Walshetal.1979；Tokso.zet.al.1979）。低频时增加较快，高频时趋于一稳定值。同时在干燥岩石中随压力的增加Q增大，主要是因为增加压力能减小岩石基质中的裂缝，从而减小摩擦。 4.

4、孔隙度 同一种砂岩，孔隙度越高，Q值越小，衰减越强；对饱和流体砂岩：衰减峰的峰位随孔隙率的增加向低温方向移动，峰值增大，峰宽变窄. 总之，衰减随孔隙率的增加而增加，呈正比关系。 5.饱和度 当岩石孔隙中饱和液体较少时，衰减随饱和度的增加呈线性增大；当饱和度再增大时，衰减随饱和度的增加而增大的速度放慢，衰减随饱和度的增加变成了非线性增大，直到衰减达到极值；当饱和度还继续增加时，衰减随饱和度的增加呈非线性下降。所以，当岩石孔隙中含有部分饱和流体时，衰减才会达到最大（1997，席道瑛）。部分饱和岩石的衰减大于完全饱和岩石的衰减。Winkler（1979）指出：对于低粘度流体（例如水、油）完全饱和的岩

5、石来说，QP QS ；对于部分饱和岩石，则QP610时,衰减迅速增加（ Peselnick吸收地层中 Q 值反演规律研究andOuterbridge，1961；GordonandDavis，1968；MckavanaghandStacey，1974；BrennanandStacey，1977；JohnstonandTokso.z，1980），这个现象与摩擦滑移这一类非线性衰减机制相联系。 8.温度 在较低温度下，与温度无关；在温度高于150度时，在石英砂中，由于岩石中的热裂缝引起衰减的增加；在孔隙流体的沸点附近，衰减随温度产生剧烈的变化（Volaroich(1957)、Davis(1968)和

6、 Spancer(1981)）。 9.埋藏深度 岩石的吸收性质与埋藏深度有关，一般随埋藏深度的增加而减小。2地震波衰减机制目前国内外已经发展了多种衰减解释理论,概括的讲包括以下几种:一是把大地当成粘弹性介质,通过应用Kelvin2Voight ,Maxwell ,SL S 等模型模拟大地性质来研究大地吸收特性;二是把大地当作不均匀介质来研究,用散射理论来解释地震波能量的衰减;三是把大地当作双相介质,相关的理论包括BIOT 理论和喷射流动理论等陆基孟把吸收衰减机制归结为两个方面26 ：一是弹性理论。该理论认为，物体在外力的持续作用下，其内部结构发生变化，在外力消失后，该物体不能完全恢复其原状，存

7、在一定的剩余应变，它的存在消耗了部分弹性能量而使地震波振幅发生衰减；二是内摩擦理论。该理论认为，介质中质点在振动过程中发生相互间的摩擦作用，使部分机械能转化为热能而消耗，使地震波振幅发生衰减。 由于地层介质的复杂性，到目前为止，尽管提出了多种衰减机制，但并没有一个统一的结论。尽管不同条件下的岩石具有不同的衰减机制，但在大多数情况下,岩石颗粒表面和岩石窄裂缝之间的摩擦是主要的衰减机制。3地震资料的非弹性衰减Q补偿的研究进展最早的反Q 滤波由D. Hale (1982 年) 提出的,依据Fut terman 提出的数学模型,利用预测误差滤波方法导出迭代求出品质因子Q 值并做Q 滤波的方法,这种方法

8、运算量大, 很难用于生产. S. H.Bickel 和R. R. Natarajan20 依据Srtick (1967) 提出的数学模型来完成反Q 滤波;虽然该方法的处理效果较好,但是它是依据积分法来完成的,因此效率很低. 1991 年Hargreaves21 提出了与Stolt 偏移相仿的相移校正方法,并且用级联式常数Q 补偿的方式开时窗逐级向下计算,每次补偿Q 值的一段,这种方法计算速度较快. 裴江云,何樵登22 利用Kjar2tansson (1979) 提出的数学模型,推导出适合于处理非平稳地震记录的反Q 滤波,该方法不仅能补偿振幅的衰减,还能消除频散. Yanghua Wang (2

9、002) 24 基于波场向下延拓理论提出了一种稳定而又有效率的反Q 滤波方法. 文中假设地层是层状介质, Q 值在每一层内是一常数,在每一层内,反Q 滤波分两步: (1) 利用一个反演系统将地表的波场延拓到该层的顶部,这一步在计算中加入了一个稳定因子; (2)在这一层内做常Q 值的反Q 滤波. 2006 年25 又将这种稳定算法推广到了Q 随时间或深度连续变化的情况.4反Q滤波的原理及方法4.1反Q滤波原理Futterman（1962年）研究了大地对地震波的吸收及频散作用，得到如下的结论： （4-1）当Q为常数时，大地吸收效应是一个逐点频衰减的低通滤波。反Q滤波就是要将因Q因素而衰减损耗的地震

10、波的能量补充回来。从地震波传播的观点来看，反Q滤波可看作是地震波的逆向传播或者是偏移过程，是将地表记录到的数据反向延拓并剔除地球介质对入射地震波所产生的Q值滤波作用。4.2反Q滤波的几种典型方法（1）常Q值模型的相位反Q滤波 该方法与Stolt偏移相似，用级联式常数Q补偿的方式，开时窗逐级向下计算，每次补偿Q值的一段，相位反Q滤波的表示式 （4-2）式中： （2） 层常Q值模型的相位反Q滤波（Bano，1996） 该方法适用于每一层的Q值是常数的地质模型，能有效校正速度频散引起的相位畸变，是稳定的，但与常Q值地层模型的相位反Q滤波一样只考虑了相位补偿而忽略了振幅的影响。 （3）层常Q值模型的全

11、反Q滤波（同时作振幅、相位补偿） 为解决相位反Q滤波的不足，Wang yang hua（2002）提出了一种层常Q值模型的全反Q滤波方法。该方法对每一个独立的常Q值地层，在地表记录的地震波场首先向下连续延拓到最新地层的顶部，然后对这个层作常Q值反滤波。在覆盖层内进行延拓时，采用向上连续延拓的方法得到稳定解。在常Q值最新地层中，振幅补偿因子是旅行时和频率的2-D函数，它可以近似表示为随时间和频率变化的两个1-D函数。 该方法分两步进行： 波场延拓 在地表和最新层顶部之间的波场延拓为 （4-3）式中：两个指数项分别表示振幅补偿因子和相位补偿因子。 在覆盖层内采用波场向上连续延括的方法 （4-4）上

12、式表示为 (4-5)式中：表示向上延拓递推因子。 波场的估计表达式为 (4-6)式中：是的共轭，是用来使解稳定的正实数，是向上延拓系统的因子，即正演地层Q滤波，可表示为 (4-7)每一层内进行反Q滤波 为清楚起见，把作为简化的时间坐标的起点，则 (4-8)令 振幅补偿因子是旅行时和频率的2-D函数，它可以近似表示为随时间和频率变化的两个1-D函数 在处定义时变项 式中：是频带的中心频率， 时间间隔内的平均值定义频变项 TV (4-9)则反Q滤波可表示为 (4-10)式中：是从M（0，w）到的重采样，J是Jacobian矩阵： 该方法虽然解决了以往方法的不足，但在求解过程中的近似表示以及每一层中

13、的Q值是常数的假设条件却增加了解的不精确性。 （4）连续变化Q值模型的全反Q滤波 为使反Q滤波效果更好，Q值更接近于真实地层情况，Wang yang hua（2006）在层常Q值模型的全反Q滤波方法的基础上又提出了连续变化Q值模型的全反Q滤波。 平面波U（x，w）传播距离Vx后可表示为 (4-11)利用波数k在反Q滤波中引入大地的Q值，并用旅行时增量V代替Vx，则反Q滤波表示式为 (4-12) 式中：，是地震频带内与最高频率有关的调谐频率，是任意参考频率处的Q值。 运用上式从地表=0至时间作波场向下延拓，得到波场的稳定表达式 (4-13)其中： 对上式求解得到下面稳定式 (4-14)式中：，是

14、稳定因子。 把所有不同频率代入上式，并对所得到的平面波进行总结，得到时间域地震记录 (4-15)此式是稳定的全反Q滤波的表达式。 为提高效率，上式的求解在Gabor变换域进行。通过公式推导，最后得到补偿后的信号 u（t） (4-16)式中： 43反Q滤波的稳定性反Q滤波能够对吸收振幅和频散相位进行补偿，其中频散相位校正项是无条件稳定的；振幅补偿项是条件稳定的，可通过设置高截频的方法来实现。通过子波响应试验表明，应根据原始数据吸收频散特点选取合适模型进行反Q滤波以保证其振幅补偿和相位校正的正确性和稳定性。反Q滤波的稳定性还受限于有效信号频带范围和Q值大小。影响反Q滤波稳定性的因素主要包括：吸收衰

15、减模型，有效信号频带范围，品质因子Q. 总体而言，反Q滤波补偿的最终稳定性取决于数据的吸收衰减程度，频率成分和最大延拓深度。我们知道地震波的优势频率是分布在一个有限频带内的，如果在一个大的频段内，由选择的将优势频率的波完成大比例因子的Q补偿，对干扰波进行小比例因子Q补偿，这样既补偿了地震波的频散吸收，提高了分辨率，又达到了压制干扰的目的。5结束语本文对地层的吸收衰减，衰减介质中的地震波正演模拟以及其反Q滤波补偿的原理和方法进行了总结。目前常用的反Q滤波方法都是在频率域进行，而频率域算子较长，计算效率低。为此在时间域进行反Q滤波将大量的频率域乘法工作转化为少量的时域褶积运算，可使计算效率大大提高

16、，进一步研究时间域反Q滤波具有很大价值。参考文献 1 Toksoz M N,Johnston D H,Timur A. Attenuation of seismic waves in dry and saturated rocks:|.Laboratory measurements. Geophysics,1979;44(4):681-690 2 Johnston D H, Toksoz M N, Timur A. Attenuation of seismic waves in dry and saturated rocks:|.Mechanisms. Geophysics,1979;44(4

17、):691-711 3 席道瑛,刘爱文,刘卫.低频条件下饱和流体砂岩的衰减研究.地震学报,1995;17(4):477-481 4 席道瑛,邱文亮,程经毅等.饱和多孔岩石的衰减与孔隙率和饱和度的关系.石油地球物理勘探,1997;32(2):196-201 5 席道瑛,刘斌,谢端等.孔隙流体饱和砂岩的衰减与频率的相关性,1998;33(1):66-77 6 徐果明，孙新雷等. 温压条件下水饱和砂岩岩样纵波波速和衰减的测量与规律. 石油地球物理勘探，2002; 37(1): 39-43. 7 凌云.大地吸收衰减分析. 石油地球物理勘探,2001;36(1):1-8 8 在干燥岩石与饱和岩石中地震波

18、的衰减. 石油物探译丛,1986. 9 张树林译.计算Q值七种方法的比较.国外油气勘探,1991;3(6):90-97 10 Yanghua Wang .A stable and efficient approach of inverse Q filtering. Geophysics,2002;67(2):657-663 11 Yanghua Wang .Quantifying the effectiveness of stabilized inverse Q filtering Geophysics,2003;68(1):337-345 12 Yanghua Wang . Inverse

19、Q-filter for seismic resolution enhancement. Geophysics,2006;71(3):51-60 13 Changjun Zhang, Tadeuse J.Ulrych. Estimation of quality factorsfrom CMP records. Geophysics，2002；67(5)：1542-1547 Geophysics,1981;6(2):152-162. 14 Hargreaves and Calvert. Inverse Q filtering by Fourier transform. Geophysics,

20、1991; 56(4): 519-527. 15 李庆忠. 走向精确勘探的道路, 石油工业出版社，1993，38-39. 16 陆基孟. 地震勘探原理(下), 石油大学出版社, 1993,336. 17 安慧,黄德济,赵宪生. VSP资料的程变反Q滤波. 物探化探计算技术,1998;20(1):19-24. 18 姚振兴, 高星,李维新. 用于深度域地震剖面衰减与频散补偿的反Q滤波方法,2003； 46(2):229-233. 19 裴江云, 何樵登. 基于Kjartansson模型的反Q滤波, 1994;9(1):91-97. 20 Arnim B.Haase, Robert R. Intr

21、insic and apparent seismic attenuation in VSP data. SEG,2006: 3472-3476. 21 刘财，刘洋等. 一种频率域吸收衰减补偿方法，石油物探.2005；44（2）：116118. 22 赵宪生等. 地震记录的Q值正演模拟与滤波，石油物探.1994；33（2）；4248. 23 郑鸿明等. 粘弹性介质中地震波的吸收衰减补偿方法，新疆石油地质.2006；27（5）：556558. 24 张钋等.射线追踪方法的发展现状，地球物理学进展石油物探.2000；15（1）：3745. 25 Futterman, W. I. Dispersive

22、 body wavesJ. Journal of Geophysics Research, 1962 , 67(13): 52795291. 26 Carpenter, E. W. Absorption of elastic waves-An operator for constant Q mechanismA. In: 36th Ann. Internat . Mtg., Soc. Expl. Geophys., Expanded AbstractsC, 1966, 418430.27 Kjartansson, E. Constant Q wave propagation and atten

23、uationJ. Journal of Geophysics Research, 1979 , 82(B9): 47374748.28 Robinson, J. C. A technique for the continuous representation of dispersion in seismic dataJ. Geophysics, 1979, 44(8): 13451351. 29 Bickel, S. H., Natarajan, R. R. Plane-wave Q deconvolutionJ. Geophysics, 1985, 50(9): 14261439. 30 B

24、ickel, S. H. Constant Q inverse filtersA. In: 51st Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded AbstractsC, 1981, 172181.31 Hale, D. Q adaptive deconvolutionA. Stanford Exploration ProjectC, 1982, 30: 133158.32 Hale, D. An inverse Q filterA. Stanford Exploration ProjectC, 1981, 28: 289298. 33

25、Hargreaves, N. D. Similarity and t he inverse Q filter: some simple algorithms for inverse Q filteringJ. Geophysics, 1992, 57(2): 944947.34 马昭军,刘洋. 地震波衰减反演研究综述J. 地球物理学进展, 2005, 20(4):10741082.Ma Z Y, Liu Y. A summary of research on seismic attenuationJ. Progress in Geophysics, 2005, 20(4):10741082.3

26、5 徐平,王宝善,张尉,林建民,陈颙. 利用互相关函数求地震波衰减J. 地球物理学报, 2006 , 49(6): 17381744.Xu P, Wang B H, Zhang W, et. al. Estimating seismic attenuation using crosscorrelation functionJ. Chinese Journal of Geophysics, 2006 , 49(6): 17381744. 36 李振春,王清振. 地震波衰减机理及能量补偿研究综述J. 地球物理学进展, 2007, 22(4): 11471152.Li Z C, Wang Q Z.

27、A review of research on mechanism of seismic attenuation and energy compensationJ. Progress in Geophysics, 2007, 22(4): 11471152.37 王珺,杨长春,乔玉雷. 用稳定高效的反Q 滤波技术提高地震资料分辨率J.地球物理学进展, 2008, 23(2): 456463.Wang J, Yang C C, Qiao Y L. Improve the seismic resolution with a stable and efficient inverse Q filter

28、J. Progress in Geophysics, 2008, 23(2): 456463.38 Holberg O. Towards optimum one-way wave propagationJ. Geophysical Prospecting, 1988, 36(2): 99114. 39 Blacquire,G. Debeye, H. W. J. Wapenaar, C. P. A. and Berkhout, A. J. 3D table-driven migrationJ. Geophysical Prospecting, 37(8): 925958.40 Hale, D.

29、3-D depth migration via McClellan transformationsJ. Geophysics, 1991, 56(11): 17781785.41 Hale, D. Stable explicit depth extrapolation of seismic wavefieldsJ. G eophysics, 1991, 56(11): 17701777.42 Zhang J F, Wapenaar K. Wavefield extrapolation and prestack depth migration in anelastic inhomogeneous

30、 mediaJ. Geophysical prospecting, 2002, 50(6): 629643.43 Thorbecke, J. W., Wapenaar, K., Swinnen, G. Design of one-way wavefield extrapolation operators, using smooth functions in WLSQ optimizationJ. Geophysics, 2004, 69(4): 10371045.44 Bachrach, R., Ferber, R., et al., 2006, Two approaches for Q estimation and its application to seismic inversionJ. New Orleans 2006 Annual Meeting, SEG. 29712975. （注：可编辑下载，若有不当之处，请指正，谢谢!）