用户名: 密码: 验证码:
井间地震纵波提取及反射波成像
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
地震勘探技术从石油勘探阶段向油气田开发阶段的渗透,使得井间地震技术得到了长足的发展。如何从复杂的井间地震记录中提取有效反射波,进行井间反射波成像是目前井间地震资料处理的研究热点。针对井间三分量采集数据,研究纵横波分离方法,提取纵波分量,进而研究反射波波动方程成像方法,对井间地震技术发展及应用具有十分重要的意义。
     本文根据高分辨率Radon变换基本原理,针对井间地震采样率高、频带宽、主频高、道间距小、单炮数据量大等特点,对高分辨率Radon正反变换方法进行了研究,设计了算法,编制了相应的处理程序,并对各参数进行了调试。根据井间地震接收系统的独特性,结合纵波和横波在不同分量上的极性关系,利用高分辨率Radon正反变换方法,研制了τ-p域纵横波波场分离方法,并编程现实。对理论模型和实际井间模型模拟数据进行处理,成功分离了纵横波波场。
     在提取纵波波场基础上,根据已有井间地震处理系统的数据结构将波动方程有限差分偏移方法移植到井间地震数据处理中。从波动方程偏移基本原理出发,结合井间地震的特点,对井间剖面进行网格划分,利用射线追踪方法求取成像时间,利用有限差分方法对井间波场进行延拓,并应用成像条件,实现了有限差分逆时偏移算法。
     模型试算和实际模拟资料处理结果表明:利用纵波波场进行成像能够更真实地反应地层产状,且其成像结果具有较高的分辨率和信噪比。
Nowdays seismic exploration technique enters the development stage of oilfiled fromoilfield prospecting stage, this gives rise to quick development of crosswell seismic tehcnique.How to get effective reflect wave from complex crosswell seismic record and crosswellseismic imaging are hot research point of crosswell seismic data processing. The acquisitionof 3C crosswell seismic, the separation of P and S waves, P waves component extracting, andwave equation imaging methods are all very important for crosswell seismic data processing.
     Based on principle of high resolution frequency domain radon transform, forward andinverse transform of high resolution Radon transform are studied and designed according tothe characters of crosswell data such as: high sample rate, wide spectrum, high mainfrequency and large single shot data quantity. Corresponding processing program aredeveloped and various parameters are tested and adjusted. Considering the specific feature ofcrosswell seismic acquisition system and the difference of P and S wave in polar nature, weutilize high resolution forward/inverse radon transform to develop the P/S wavefieldseparation method in? ? p domain, and design corresponding program. Processing ofsimulated data of theoretical model and practical crosswell model separates P/S wavefieldsuccessfully.
     In the basis of extracting P wavefield, depending on data structure of existing crosswellseismic processing system, we transplant wave equation finite-difference migration methodinto crosswell seismic data processing. Considering the character of crosswell seismic, wegrid the crosswell section, calculate imaging time by ray-tracing, use finite-difference methodto achieve crosswell wavefield reverse time continuation, and finally apply imaging conditionto get finite-difference reverse-time migration profile.
     Processing results implies that: imaging from extracted P wavefiled can reflect actualattitude of strata more clearly than image from full P/S wavefield, the imaging results havehigh resolution and high S/N.
引文
[1]宋建国.井间地震技术在油田开发中的应用潜力.石油地球物理勘探,2003.8
    [2] Payton,C. A new tool for reservoir management.Geophysics– New Horizons inExploration&Production,HartEnergyPublication,August1999:2-23
    [3] Gerhard R. Pratt§,Richard M.. Seismic waveform inversion in the frequency domain,Part 2: Fault delineation in sediments using crosshole data. Geophysics,1999,64(3):902~914
    [4]王世库,刘清林.综合利用井间地震等多种资料确定精细地质模型.石油物探,1995,34(4):10~13
    [5] Joe Wong. Crosshole seismic imaging for sulfide orebody delineation near Sudbury,Ontario,Canada.Geophysics,2000;65(6):1900~1907
    [6] Deans S. R. The Radon transform and some of its applications[M]. News York:JhonWileyandSonsInc,1983
    [7] Phinney R. A., Chowdhury K. R., Frazer L N. Transformation and analysis of recordsections[J].J.Geophysics.Res,1981,86(3):359-377
    [8] Durrani T. S., Bisset . The Radon Transform and its properties[J] . Geophysics .1984,49(11):1180-1187
    [9] Chapman C. J. Generalized Radon Transforms and slant stacks[J] . Geophysics J. Roy .Astr.S.oc.,1981,54(4):481-518
    [10] Thorson J. R., Claerbout J. F.. Velocity stack and slant stack stochastic inversion[J]Geophysics,1985,50(12):2727-2741
    [11] Kostov C.. Toeplitz structure in slant stack inversion[A]. In 60th Annual Internat.Mtg,Soc.Expl.Geophysics,1990,ExpandedAbstrancts[C]1618-1621
    [12] Moster J. D., Mosher C. C.. Suppression of multiples reflection using thr Radontransform[J].Geophysics,1992,57(3):386-395
    [13] Zhou B., Greenhalgh S. A.. Linear and parabolic z-p revisited[J]. Geophysics, 1994,59(8):1133-11149
    [14] Qiang Li.High resolution hyperbolic Radon transform multiple removal . the universityofAlberta,2001
    [15]吴律. ? ?p变换及应用.石油工业出版社,1993
    [16]曹务祥,周国安.一种新的波场分离方法.石油物探,1997,36(4)
    [17]李彦鹏,马在田.坐标拉伸后的线性拉冬变换法波场分离.石油地球物理勘探,1998,33(5)
    [18]张秉铭.固体液体分界面处P波和S波波场的分离.地球物理学进展,1999,14(4)
    [19]胡天跃.多分量地震波波场分离研究,地球物理学报,2004,47(4)
    [20]刘喜武.高分辨率Radon变换方法及其在地震信号处理中的应用.地球物理学进展,2004,19(1)
    [21]吴律.层析基础及其在井间地震中的应用.北京石油工业出版社,1997
    [22] Lazaratos, Spyros K., Rector, James W., Harris, Jerry M.. Schaack, Mark Varm,High-resolution, cross-well reflrction imaging: potential and technical difficulties,Geophysics,1993,V58(9): 1270-1280
    [23] Khalil,Ashraf A.,Stewart, Robert A.,et al. Full-waveform processing and interpretation of1 kilohertyz cross-well seismic data ,Geophysics , 1993, V58(9): 1248-1256
    [24] Smalley Nicholas, Harris J.. CDP stacking imaging of crosswell reflection. ExpandedAbstracts of 62nd SEG Mtg,1992
    [25]严又生等.井间地震测量数据的反射波成像及解释.石油地球物理勘探,2000,35(1):37~41
    [26]宋建国.井间地震偏移速度扫描的反射波超级叠加技术.石油大学学报,1999,23(1):23~26
    [27] Claerbout J.F.地震物理数据处理基础.北京石油工业出版社,1976
    [28] Samuel H., Gray, etal. Seismic migration problem and solutions. Geophysics, 2001,66(5):1622~1640
    [29]刘喜武,刘洪。波动方程地震偏移成你方法的现状与进展。地球物理学进展,2002,17(4): 582~590
    [30] Schneider W.A.. Integral formulation for migration in two and three dimensions.Geophysics,1978,43(1): 49~76
    [31] Stolt R.H..Migration by Fourier Transform. GeoPhysics,1978,43(1): 23~41
    [32] Gazdag J.. Wave-equation migration with the Phase shift method. GeoPhysics,1978,43(7): 13 42~135
    [33] C1aerbout J. F. Towards aunified theory of reflection mapping,GeoPhysics, 1971,36(4 ) :46 7~581
    [34] Claerbout J.F., Ooherty S.M.. Downward continuation of move-out corrected seismogram . GeoPhysics,1972,37(5): 741~746
    [35] Beistein N.. On imaging of reflectors in the earth. GeoPhysics,1987,52: 931~942
    [36] Hill N R. Gaussian beam migration. GeoPhysics,1990,55(9): 1416~1428
    [37] Hill N R. Pre-stack Gaussian beam depth migration. GeoPhysics,2001, 66(5): 1240~1250
    [38] Chang W.F., McMechan G. A.. 1986,Reverse migration of offset vertical seismicprofiIing data using the excitation time imaging condition. GeoPhysics,51,67~84
    [39]马在田.高阶近似有限差分偏移.石油地震物理勘探,1982,17(1):5~15
    [40]张关泉.利用低阶偏微分方程组大倾角差分偏移.地球物理学报,1986,29(3)173~282
    [41] Li Z.. Compensating finite-diference errors in 3-D migration and modeling. Geophysics,1991, 56(10): 1650~1660
    [42] Berkout A. J.. Steep dip finite-diference migration. Geophysical Prospecting, 1978,27(2):196~213
    [43]陆基孟.地震勘探原理.东营:石油大学出版社(第2版):114~117
    [44] Gazdag J., Sguazzero P..Migration of seismic data by phase shift plus interpolation.Geophysics, 1984,49(2): 124~131
    [45] Stofa P.L. et al. Split-step Fourier migration. Geophysics,1990,55(4): 410~421
    [46] Ristow D.,Ruhl T..Fourier finite-diference migration.Geophysics,1994,59(12): 1882~1893
    [47]谢里夫R.e.,吉尔达特L.P..勘探地震学.北京:石油工业出版社,1999
    [48] Chang W. F., McMechan G. A.. Reverse migration of offset vertical seismic profiIingdata using the excitation tim eimaging condition.Geophysics,1986,51: 67~84
    [49] Sun R., McMechan G…A pre-stack reverse-time migration for elastic with applicationto synthetic offset vertical seismic profiles. Proc, Inst. Electr. And Electron.Eng.,198671,457~465
    [50] Dong, Zhengxin; McMechan, George A. 3~D prestack migration in anisotropic media.Geophysics, 1993,58(1): 79~90
    [51]何兵寿.矢量波场叠前逆时深度偏移方法研究.石油大学博士学位论文,2002
    [52] Xianhuai Zhu, McMechan G.A.. Acoustic modeling and migration of stacked cross-holedata.1988,53(4): 492~500
    [53] Balch A.H.;Chang,H., Hofland,G. S.et al. Use of forward and bac-scattered P and S waveconverted waves in cross~borehole imaging. GeoPhysical Prospecting,1991,39(7):887~913.
    [54]张昕.井间地震反射波迭前深度偏移方法研究.石油大学硕士学位论文,1996
    [55] Rector J.W.. Cross-well methods: Where are we , where are we going? Geophysics ,1995,60(3): 627~630
    [56]陈世军,刘洪,周建宇,何惺华.井间地震技术的现状于展望.地球物理学进展,2003,18(3):524~529
    [57] Williams M.C..井间地震成像:这种技术的全盛时期来到了吗?邵祝华译.石油物探译丛,1997,(2): 39~46
    [58]张敏,李振春.偏移速度分析于建模方法综述.勘探地球物理进展,2007,30(6)
    [59] Carswell A., Tang R., Dillistone C., et al. A new method of Field Separation in VSP DataProcessing[J]. Expanded Abstracts of 54th Annual International SEG Mtg, 1984: 40~42
    [60] Wooil Moon,Allan Carswell, Roger Tang, et al.Radon transform wave field separationfor vertical seismic profiling data. Geophysics, 1986, 51(4): 940~947
    [61] M.D.Sacchi. Statistical and Transform Methods for Seismic Signal Processing.University of Alberta, 1999:169~207
    [62] Orhan Yilmaz, M.Thurhan taner. Discrete plane-wave decomposition by least-meansquare-error method. Geophysics, 1994, 59(6): 973~982
    [63] Kostov C.. Toeplitz structure in slant - stack inversion [J]. Expanded Abstracts of 60thAnnual Internat SEG Mtg, 1990: 1618~1621
    [64] Sacchi M., Ulrych T.. High-resolution velocity gathers and offset space reconstruction [J].Geophysics, 1995, 60(4): 1169~1177
    [65]许世勇,李彦鹏,马在田等.? ? p变换法波场分离.中国海上油气(地质),1999;13(5):334~334
    [66]戚敬华,李萍.利用? ? p变换技术实现多波波长分离.煤田地震与勘探,1998;26(5):54~57
    [67]孙显义,陈可为,许世勇.采用? ? p变换法进行纵、横波波场分离.大庆石油地质于开发,2002;21(4):76~77
    [68]戴华林.快速高分辨率抛物线拉东变换.硕士学位论文,长安大学西安:2001
    [69]张卫红,陈琳,髙志淩.肯71井区三维VSP资料波场分离方法应用研究.石油物探,2006;45(5):532~536
    [70]姚忠瑞,何惺华,左建军等.多方位Walk-away VSP处理方法[J] .石油物探,2006;54(4):380~384
    [71]陈遵德,段天友等.SVD滤波方法的改进及应用[J] .石油地球物理勘探,1994;29(6):783~792
    [72] J.Mars, N. Thirion. Wave field Separation Based On A Bland Algorithm. SEG TechnicalProgram Expanded Abstracts, 1996: 1547~1550
    [73] Jean-Luc Boelle, Guy Canadas, Pierre Hugonnet. Wave filed separation in boreholeseismic by linear Radon Decomposition. SEG Expanded Abstracts, 1998: 272~275
    [74] Trad, Ulrych, Sacchi M D. Accurate interpolation with high resolution time-variantRadon transforms [J].Geophysics, 2002, 67(4): 44~656
    [75] Daniel Trade, Tadeusz Ulrych, Mauricio. Sacchi. Latest views of the sparse Radontransform. Geophysics, 2002, 68(1): 386~399
    [76] Kurt J. Marfurt, Robert V. Schneider, Michael C. Mueller. Pitfalls of using conventionaland discrete Radon transforms on poorly sampled data. Geophysics, 1996, 61(5):1467~1482
    [77] Mark Ng, Mike perz. High resolution Radon transform in the t-x domain using“intelligent”prioritization of the Gauss-Seidel estimation sequence. SEG Int’lExposition and 74th Annual Meeting Denver, Colorlado,2004: 2160~2163
    [78] Katarina Jovanovic, Kurt.J.Marfurt. P and S waves separation from vector VSP’s by ablocky antialias Discrete Radon Transforms: Application to Vinton Dome, LA. SEG Int’lExposition and 74th Annual Meeting Denver, Colorado , 2004: 2513~2516
    [79] Peter W. Cary. The simplest discrete Radon transform. SEG Expanded Abstracts, 1998:1999~2002
    [80] Douglas J. Foster. Charles C. Mosher. Up and Downgoing Wave Field Separation forMulti-offset Single Level VSP Data. SEG Technical Program Expanded Abstracts, 1988:816~818
    [81] Mauricio Dino Sacchi. Aperture Compensated Radon and Fourier Transforms. TheUniversity of British Columbia, 1996: 53~99
    [82] Douglas J.Foster, James E. Gaiser. Elastic Wave Field Decomposition of offset VSP Data.SEG Technical Program Expanded Abstracts ,1986: 563~565
    [83] Sacchi M.D.. Statistical and Transform Methods for Seismic Signal Processing.University of Alberta, 1999:169~207
    [84]陈建伟.用高分辨率双曲Radon变换实现波场分离.硕士学位论文,长安大学硕士学位论文,西安:2005
    [85]张会星,宁书年.弹性波动方程叠前逆时偏移.中国矿业大学学报,2002.9
    [86]马在田.地震成像技术—有限差分偏移[M] .北京:石油工业出版社,1989
    [87]张秉铭,董敏煜,李承楚,等.各向异性介质中弹性波多分量联合逆时深度偏移.长春科技大学学报,2000.1
    [88]梁非,杨慧珠等.横观各向同性介质中弹性波的吸收边界条件.固体力学学报,1994,15(4):311~318
    [89] Randall C J. Absorbing boundary conditions for the elastic wave equation. Geophysics,1988,53:611~624
    [90]傅旦丹,朱宏彰.复杂介质多波地震正演模拟的动力学射线追踪法.中国海上油气(地质),1999,13(5):328~333
    [91]赵成斌,许云等.地震波走时和射线的有限差分计算.1995,13(4):41~46
    [92] Mark E.Mathisen, Paul Cunningham, Jesse Shaw, Anthony A. J. H. Justices,and N. J.Guinzy§Crosswell seismic radial survey tomograms and the3-D interpretation of aheavy oil steamflood.Geophysics, 1995,60(3): 651-659

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

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

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