基于地震波衰减理论的地层吸收参数提取方法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
随着我国油气勘探、开发工作的不断发展,地震勘探重点由原来的构造油气藏向岩性油气藏、隐蔽油气藏转移,非均匀粘弹性介质中的地震波衰减也因此成为近年来研究的热点。尤其是在隐蔽岩性天然气藏的预测中,地层吸收参数被认为是最为敏感、可靠的依据之一。因此,研究、估算和补偿地下介质对地震波的吸收不仅对叠前地震资料处理、提高地震成像精度、合理解释AVO效应和正确反演介质物理属性等是极其重要的,而且近几年来,地层吸收参数作为一个地层固有属性在预测岩性、油气藏的位置和范围方面具有更重要的意义。
     根据目前地震波衰减理论的研究现状及发展方向,结合实际生产中天然气藏预测技术的需要,论文在对地震波衰减特征进行数值模拟分析的基础上开展了地层吸收参数提取方法的研究。首先,在研究讨论了粘弹性介质基本理论的基础上,分别采用交错网格高阶有限差分法与单程波数值模拟方法进行了粘弹性介质的地震波传播数值模拟,从不同角度分析了粘弹性介质中地震波传播特征,总结了地震波各种传播特征随着地层吸收影响的变化规律;其次,利用目前较为先进的谱分解技术(S变换)推导了地震波传播方程在S域的能量密度函数,并且在S域分析得到了地震波能量、频率等信息变化与Q值之间的对应关系,从而结合前面数值模拟得到的结论规律,在时频域提出了三种不同的地层吸收参数提取方法,并在此基础上给出了地层吸收综合参数的概念;此外,为了从视衰减中分离出地层的固有衰减属性,论文还探索性的研究了利用基于贝叶斯理论的非线性波形反演来进行Q值反演的方法,以及在AVO响应分析中加入Q值的影响因素,进行AVO-Q分析,得到了初步的认识成果。最后,将上述的各种方法在模型数据和东部某油田中深层天然气藏勘探中做了检验和应用。应用效果表明,文中讨论的地层吸收参数提取方法是可行和有效的。
With development of oil and gas exploration, more and more research has been done on the lithologic and concealed hydrocarbon reservoir rather than structural hydrocarbon reservoir. Thus studying on the seismic wave attenuation in viscoelastic medium has becoming a hotspot in recent years. Especially in the prediction of concealed lithologic natural gas reservoirs, absorption parameter was regarded as one of the most sensitive and reliable indicators. Therefore, study, evaluate and compensate for the absorption of seismic waves by underground medium is of great significance for processing prestack seismic data, improving the seismic imaging precision, explaining AVO effect logically and inversing physical attributes correctly. In recent years, earth layer absorption parameters as an intrinsical attribute of earth layer, has a significant meaning in prediction of lithology, location and range of oil and gas reservoirs.
     According to the development of seismic attenuation theories, and the technical requirements of natural gas prediction in real production, this thesis investigate the extraction of earth layer absorption parameters based on numerical simulation analysis the characteristic of seismic attenuation .First of all, based on the viscoelasticity medium theory, using two different ways the staggered-grid high-order finite difference method and one-way wave-equation numerical simulation method, seismic wave propagation of viscoelastic medium were simulated in this paper. From different aspect, we analyzed seismic wave characteristics and changes of seismic waves affected by earth layer absorption. Secondly, using a more advanced spectrum decomposition technology (S transform) inferred energy density function of the seismic wave propagation equation in the S domain, and analyzed the relationship between seismic wave energy, frequency and the Q value in the S domain. Combined with the numerical simulation we have conducted before, three different methods of earth layer absorption parameters extraction were proposed in time-frequency domain, and the concept of earth layer absorption comprehensive parameters were given in this paper. Moreover, in order to separate the intrinsic attenuation attribute from the decay, a Q value inversion method was put forwarded based on nonlinear waveform inversion of Bayes theory. Meanwhile, the Q value influencing factor was also added into the AVO response analysis to do the AVO-Q a nalysis in this study.
     Based on the above research, all the research methods and models were tested and applied to the exploration of middle-deep layer anatural gas reservoirs in an oilfield in east china. The result shows that methods of extracting attenuation parameters discussed in this paper are feasible and effective.
引文
[1]刘祖沅.岩石中地震波衰减的实验研究[J].地震科学研究,1983,6:34-40
    [2] Ricker , Norman. The form and laws of propagation of seismic wavelets[J]. Geophysics, 1952, 18:10-40
    [3]杨伟刚.基于地震波衰减属性的主参数分析研究[D].成都:成都理工大学,2007
    [4]俞寿朋.高分辨率地震勘探[M].北京:石油工业出版社,1993:35-56
    [5]陆基孟.地震勘探原理[M].山东东营:石油大学出版社,1993:261-270
    [6] ShengJie Sun. Determination of seismic attenuation from surface and downhole measurements[D]. Oklahoma:University of Oklahoma Graduate College,2000
    [7] Biot M. A .,1956, Theory of propagation of elastic waves in a fluid saturated porous solid ,II, Higher frequency range[J]. J ,Acoust. Soc. Am., v .28,p. 179-191
    [8] McDonal F J,Angona F A,Mills R L,et al. Attenuation of shear and compressional waves if Pierre shale[J]. Geophysics, 1958, 23(2):421-439.
    [9] Futterman W I. Dispersive body waves[J]. Geophys Res,1962,7(13): 5279-5291
    [10] White J E. Comuted seismic speeds and attenuation in rocks with partial saturation[J]. Geophysics, 1975, 40(2):224-232
    [11] O’Connel R J,Budiansky B. Viscoelastic properties of fluid-saturated cracked solids[J]. J. G.R., 1977, 76:2022-2034
    [12] Macko G M, Nur A. Wave attenuation in partially saturated rocks[J]. Geophysics, 1979, 42(2):161-178
    [13] Palmer I D, Treaviolia M L. Attenuation by squirt flow in undersaturated gas sands[J]. Geophysics, 1980, 45(11):1780-1792
    [14] Jones T D. Pore fluids and frequency-dependent wave propagation in rocks[J]. Geophysics, 1986,51(12):1939-1953
    [15] Murphy W F,Winkler K W,Kleinberg R L. Acoustic relaxation in sedimentary rocks:Dependence on grain contacts and fluid saturation[J]. Geophysics, 1986, 51(4): 757-766
    [16] Dvokin J,Nur A. Dynamic poroelastisity: a unified model with the squirt and the Boit mechanisns[J]. Geophysics, 1993, 58(3):524-543
    [17] Toksoz M N,Jackson D H,Timur A. Attenuation of Seismic wave in dry and saturated rocks I:Laboratory measurements[J].Geophysics, 1979, 44(4):671-680
    [18] Dutta N C,Ode H. Attenuation and dispersion of compressional waves in fluid-filled porous rocks with partial gas saturation(White model)-Part I:Biot theory[J].Geophysics, 1979, 44(11):1777-1788
    [19] Grant A G.. Fluid effect on velocity and attenuation in sandstone[J]. J. Acoust Soc.Am., 1994, 96(2):1158-1173
    [20]马昭军,刘洋.地震波衰减反演研究综述[J].地球物理学进展,2005,20(4):1074-1082
    [21] Carcione J M. Wave propagation in a linear viscoustic medium[J]. Geophys. J. Roy. Astr. Soc., 1988, Vol.93
    [22]韩文功,印兴耀等.地震技术新进展(上)[M].山东东营:中国石油大学出版社,2006:73-122
    [23] Carcione J M,Kosloff D,Kosloff R. Viscoacoustic wave propagation simulation in the earth[J].Geophysics,1988,53:769-777
    [24] Tal-EzerH,Carcione J M,Kosloff D. An accurate and efficient scheme for wave propagation in linear ciscoelastic media[J]. Geophysics, 1990, 55:1366-1379
    [25] Robertsson J O ,Blanck J O ,Symes W W. Viscoelastic finite-defference modeling[J].Geophysics,1994,59(9):1444-1456
    [26]孟凡顺,郭海燕等.复杂地质体粘弹性波正演模拟的有限差分法[J].青岛海洋大学学报,2000,30(2):315-320
    [27] Tong Xu,George A.McMechan. Efficient 3-D vciscoelastic modeling with application to near-surface land seismic data[J].Geophysics, 1998, 63(2): 601-612
    [28] Stekl,R.G.Pratt.Accurate viscoelastic modeling by freduency-domain finite differences using rotated operators[J]. Geophysics, 1998, 63(5): 811-823
    [29]Qingbo Liao,George A.McMechan*.Multifrequency viscoacoustic modeling and inversion[J]. Geophysics, 1996, 61(5): 524-522
    [30] Adam P.Koesoemadinate,George A.McMechan*. Empirical estismation of viscoelastic seismic parameters from petrophysical on sandstone[J]. Geophysics, 2001, 66(5):294-307
    [31] Upendra K.Tiwari, George A.McMechan. Effects of incomplete parameterization on inversion of full wavefield viscoelastic seismic data for petrophysical reservoir properties[R].SEG Int’l Exposition and 74th Annual Meeting *Denver, Colorado* 10-15 October 2004: 274~283
    [32]单启铜.粘弹性波动方程正演模拟与参数反演[D].山东东营:中国石油大学(华东),2007
    [33]奚先,姚姚.二维粘弹性随机介质中的波场特征分析[J].地球物理学进展,2004,19(3):608-615
    [34]刘财,张智,邵志刚等.线性粘弹体中地震波场伪谱法模拟技术[J].地球物理学进展,2005,20(3):640-644
    [35] Bath M .Spectral Analysis in Geophysics[M]. New York:Elsevier,1974.
    [36] Gladwin M T, Stacey F D. Anelastic degradation of acoustic pulses in rock[J]. Phys.Earth Planet Int., 1974,8(2):332-336.
    [37] Kjartansson E.Constant Q wave propagation and attenuation[J].J. G. R., 1979, 84(4): 4737-4748
    [38] Kan T K,Corrigan D,Huddleston P D.Attenuation measurement from VSP[R]. Presented at 51th Ann. Internat.Mtg. Soc. Expl. Geophys.1981
    [39] Quan Y L ,Harrisy J M..Seismic attenuation tomography using the frequency shift method[J].Geophysics,1997,62(3):895-905
    [40] Sun Sh J, Determination of seismic attenuation from surface and downhole measurements [D]. Oklahoma:University of Oklahoma,2000
    [41] Rainer Tonn. The determination of the seismic quality factor Q from VSP data: A comparison of different computational methods[J]. Geophysical Prospecting, 1991, 39(1):1-27
    [42] Dasgupta R, Clarkz R A. Estimation of Q from surface seismic reflection data[J]. Geophysics, 1998, 63(12):2120-2128
    [43] Hackert C L, Parra J O. Improving Q estimation from seismic reflection data using well-log-based localized spectral correction[J].Geophysics,2004,69(6):1521-1529
    [44] Zhang Ch J, Ulrychz T J.Estimation of quality factors from CMP record[J]. Geophysics, 2002, 67(5):1542-1547
    [45] Elapavuluri Pavan, Bancroft John.Estimation of Q using crosscorrelation[R]. CSEG National Convention, 2004, 1-4
    [46]周辉,渠广学,杨宝俊.用地震波频谱计算Q值的新方法[J].长春地质学院学报,1994,24(4):461-467
    [47]刘学伟,邰圣宏.一种考虑噪音干扰的地层风化层Q值反演方法[J].石油地球物理勘探,1996,31(3):367-373
    [48]王辉,常旭,刘尹克.时间域相邻道地震波衰减成像研究[J].地球物理学报,2001,44(3):396-403
    [49]王西文,杨孔庆,杨午阳等.基于小波变换拟瞬时吸收系数的计算方法和应用[J].石油地球物理勘探,2002,37(增刊):205-209
    [50]李宏兵,赵文智,曹宏等.小波尺度域含气储层地震波衰减特征[J].地球物理学报,2004,47(5):892-898
    [51] Lancaster, S. J., and M. C. Tanis. High density 3D prestack Q estimation[R].Presented at the 66th Annual Conference and Exhibition, EAGE. 2004
    [52]牛滨华,孙春岩.半空间介质与地震波传播[M].北京:石油工业出版社,2002
    [53]张静.线性粘弹介质中的频散和衰减特征研究[D].北京:中国地质大学(北京),2006
    [54] Jose M . CARCIONE. Wave Fields In Real Media: Wave Propagation In Anisotropic, Anelastic And Porous Media[J]. Seismic Exploration, 2001, 59-73
    [55] Christensen,R.M. Theory of Vioscoelasticity-An Introduction[M]. San Diego:cademic Press, 1982
    [56] Blanch,J.O., Robertsson, J.O.A.&Symes,W.W. Modeling of a constant Q:ethodology and algorithm for an efficient and optimally inexpensive viscoelastic technique[J]. Geophysics, 1995, 60:176-184
    [57] Stig Hestholm. Three-dimensional finite difference viscoelastic wave modeling including surface topography[J]. Geophys.J.Int, 1999, 139: 852-878
    [58] Robertsson J.O.A., Blanch,J.O.&Symes, W.W. Viscoelastic finite-difference modeling[J]. Geophysics, 1994, 59: 1444-1456
    [59]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419
    [60] Gazdag, J., Wave equation migration with the phase-shift[J]. Geophysics, 1978, 43:1342-1351
    [61]吴国忱.各向异性介质地震波传播与成像[M].山东东营:中国石油大学出版社,2006:87-155
    [62] Carcione J M, Helle H B. Numerical solution of the poroviscoelastic wave equation on a staggered mesh[J]. J. Comput. Phys., 1999,15(4):520-527
    [63]牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社,2005
    [64] Hokstad, K. 3-D elastic finite-difference modeling in tilted transversely isotropic media[A]. 72th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 2002, 1951-1954
    [65] Collino F and Tsogka C. Application of the perfectly matched absorbing layer modelto the linear elastodynamic problem in anisotropic heterogeneous media [J]. Geophysics, 2001, 66(1): 294-307
    [66]熊晓军.单程波动方程地震数值模拟新方法研究[D].四川成都:成都理工大学,2007
    [67] Futterman W I. Dispersive body waves[J]. Geophys Res,1962,7(13):5279-5291
    [68]赵宪生等.地震记录的Q值正演模拟与反滤波[J].石油物探,1994,33(2):42-48
    [69] Lars Hübert. Seismic Attenuation and Hybrid Attributes to reduce exploration risk– North Sea Case Study[A]. Rock Solid Images
    [70] R. G. Stockwell, L. Mansinha, and R. P. Lowe. Localization of the Complex Spectrum:The S Transform[J]. IEEE TRANSACTIONS ON SIGNAL PROCESSING, 1996, 44(4):998-1001.
    [71]邹文. S-变换时频分析技术及其在地震勘探中的应用研究[D].北京:中国地质大学,2005
    [72] JnassneA. Oality Porpeyrt of the Gaussina wnidows Pecortgrma[J]. IEEE Transactions on Acoustics,Speech nad Signal Proecssing,1991,ASSP-39(1)
    [73] Farge M . Wavelet transforms and their application to turbulence[J]. Annual Review of Fluid Mechanics, 1992, 24:395-457
    [74] Aki K, Richards PG..Quantitative Seismology, Theory and Methods[M].San Francisco: W H Freeman, 1980
    [75]印兴耀,刘杰,杨培杰.一种基于负熵的Bussgang地震盲反褶积方法[J].石油地球物理勘探,2007,42(5):499-505
    [76] Joel Walls, M. T. Taner. Novel Use of P- and S-wave Seismic Attenuation for Deep Natural Gas Exploration and Development[Z]. Rock Solid Images,2006
    [77] Scott Singleton,M. Turhan Taner, Sven Treitel. Q Estimation using Gabor-Morlet Joint Time-Frequency Techniques[A]. Rock Solid Images,2006
    [78]周静毅.基于地震剖面特征的属性分析[D].山东东营:中国石油大学(华东),2005
    [79]李庆忠.走向精确勘探的道路[M].北京:石油工业出版社,1993
    [80]陈健江.AVO三参数反演方法研究[D].山东东营:中国石油大学(华东),2007
    [81]王芳,陈生昌.解非线性反演问题的新策略[J].石油物探,1999,38{3):76-8
    [82] Tarantola A. Inverse problem theory:methods for data fitting and model parameter estimation[M]. Elsevier Science Publishers,1987
    [83] Steven R. Rutherford and Robert H. Williams. Amplitude-Versus–Offset variations in gas sands[J]. Geophysics,1989,54:680-688
    [84] Hongbo Zhang and R. James Brown. A review of AVO analysis[R]. CREWS Research Report,2001,13:357-377

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700