三分量井孔地震资料波场分离与应用
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
当今勘探和开发工作面临的挑战需要高精度?高清晰度和高可信度的技术,井孔地震的研究日益受到广泛的应用,尽管目前无论垂直地震剖面(Vertical Seimic Profile,简称VSP)还是井间地震的成本较其他技术来说是比较高的?从复杂的井孔地震记录中提取纵横波场,进行有效波成像日益关键?
     本文针对井间三分量采集数据,研究纵横波分离方法,提取纵横波分量,进而研究反射波波动方程成像方法?根据高分辨率Radon变换基本原理,针对井孔地震采样率高?频带宽?主频高?道间距小等特点,对高分辨率Radon正反变换方法进行了研究,编制了相关程序,调试参数使收敛效果最佳?根据井中三分量接收系统的独特性,结合纵波和横波在不同分量上的极性特性,在Radon域用极化滤波对纵横波进行分离,并编程现实?对理论模型和实际数据进行处理,成功分离了纵横波波场?然后对提取的纵横波波场,进行井间波动方程有限差分偏移成像?从波动方程偏移基本原理出发,结合井间地震的特点,对井间剖面进行网格划分,利用射线追踪方法求取成像时间,利用有限差分方法对井间波场进行延拓,并应用成像条件,实现了有限差分逆时偏移算法?模型和实际资料处理结果表明纵横波波场成像能够真实地反应地层产状,且其成像结果具有较高的分辨率和信噪比?通过对比发现,横波成像结果要比纵波成像结果好,更充分地利用了井中的地震波场信息?
Nowdays the challenge of seismic E&D needs the techniques of high accuracy,high definition,high confidence;this gives rise to quick development of borehole seismic technique,though at present the cost of no matter VSP or crosswell is higher.How to get P/S wave from complex borehole seismic record and crosswell seismic imaging are hot research point of crosswell seismic data processing.
     Based on the acquisition of 3C borehole seismic, we research the separation of P and S waves, extract P and S wave component, and study wave equation imaging methods. Depending on principle of high resolution frequency domain radon transform, according to the characters of crosswell data such as: high sample rate, wide spectrum, high main frequency quantity, we research high resolution radon transform; corresponding processing program are developed and various parameters are tested and adjusted. Considering the specific feature of 3C borehole seismic acquisition system and the difference of P and S wave in polar nature, we utilize high resolution forward/inverse radon transform to develop the P/S wavefield separation method in ? ?pdomain, and design corresponding program.
     Processing of simulated data of theoretical model and practical crosswell model separates P/S wavefield successfully. In the basis of extracting P and S wavefield, depending on data structure of existing crosswell seismic processing system, we transplant wave equation finite-difference migration method into crosswell seismic data processing. Considering the character of crosswell seismic, we grid the crosswell section, calculate imaging time by ray-tracing, use finite-difference method to achieve crosswell wavefield reverse time continuation, and finally apply imaging condition to get finite-difference reverse-time migration profile. Processing results implies that imaging from extracted P/S wavefiled can reflect actual attitude of strata actually, and the imaging results have high resolution and high S/N.
引文
[1]杨勤伟,徐丽萍.地震勘探技术新进展.勘探地球物理进展,2002, 25(1):5-10
    [2]张胜.VSP仪器现状及发展展望.物探装备,2008, 18(2):86-93
    [3]郭建.VSP技术应用现状及发展趋势.勘探地球物理进展,2004, 27(1):1-8
    [4]曹辉.井间地震技术发展现状.勘探地球物理进展,2002, 25(6):6-10
    [5] Payton.A new tool for reservoir management.Geophysics-New Horizons in Exploration & Production.Hart Energy Publication,1999, 8(9):2-23
    [6]宋建国.井间地震技术在油田开发中的应用潜力.石油地球物理勘探,2003,37(4):377-382
    [7] Joe Wong.Crosshole seismic imaging for sulfide orebody delineation near Sudbury. Ontario,Canada.Geophysics, 2000, 65(6):1900-1907-
    [8] Deans S. R. The Radon transform and some of its applications[M]. News York.Jhon Wiley and Sons Inc,1983, 12(1):120-126
    [9] Durrani T. S..Bisset.The Radon Transform and its properties[J].Geophysics,1984, 49(11): 1180-1187
    [10] Chapman C. J. Generalized Radon Transforms and slant stacks[J].Geophysics J.Roy .Astr. S. oc.,1981, 54(4): 481-518
    [11] Thorson J. R., Claerbout J. F.. Velocity stack and slant stack stochastic inversion[J] Geophysics,1985,50(12): 2727-2741
    [12] Kostov C.Toeplitz structure in slant stack inversion[A].In 60th Annual Internat. Mtg, Soc. Expl.Geophysics,1990 , Expanded Abstrancts[C]:1618-1621
    [13] Beylkin G. Discrete radon transform[J].IEEE Trans Acoust,Speech and Signal Processing, 1987, 35(2):162-172
    [14] Moster J. D.,Mosher C. C.Suppression of multiples reflection using thr Radon transform[J]. Geophysics,1992,57(3): 386-395
    [15] Zhou B., Greenhalgh S. A.. Linear and parabolic z-p revisited[J]. Geophysics,1994,59 (8): 1133-1149
    [16] Qiang Li.High resolution hyperbolic Radon transform multiple removal. the university of Alberta,2001:51-64
    [17]吴律. ? ?p变换及应用.石油工业出版社,1993:18-39
    [18]李远钦.一种非线性Radon变换及非零偏移距VSP波场分离.石油物探,1994,33(3):33-39
    [19]曹务祥,周国安.一种新的波场分离方法.石油物探,1997,36(4):104-109
    [20]张秉铭,董敏煜,滕吉文等.固体液体分界面处P波和S波波场的分离.地球物理学进展,1999:14(4):18-27
    [21]刘喜武,刘洪,李幼铭.高分辨率Radon变换方法及其在地震信号处理中的应用.地球物理学进展,2004,19(1):8-15
    [22]曾有良,乐有喜,单启铜等.基于高分辨率Radon变换的VSP波场分离方法.石油物探, 2007, 46(2):115-120
    [23] Schneider W.A.. Integral formulation for migration in two and three dimensions. Geophysics,1978, 43(1): 49-76
    [24] Gazdag J.. Wave-equation migration with the Phase shift method.GeoPhysics,1978, 43 (7): 1342-1355
    [25] Stolt R.H..Migration by Fourier Transform. GeoPhysics,1978, 43(1): 23-41
    [26] C1aerbout J. F. Towards aunified theory of reflection mapping.GeoPhysics,1971,36 (4 ) :467-581
    [27] Claerbout J.F.,Ooherty S.M..Downward continuation of move-out corrected seism ogram. GeoPhysics, 1972, 37(5): 741-746
    [28]谢里夫,吉尔达特L.P..勘探地震学.北京:石油工业出版社,1999
    [29] Chang W. F., McMechan G. A.. Reverse migration of offset vertical seismic profiIing data using the excitation tim eimaging condition.Geophysics,1986, 51(6): 67~84
    [30] Sun R., McMechan G…A pre-stack reverse-time migration for elastic with application to synthetic offset vertical seismic profiles. Proc, Inst. Electr. And Electron.Eng,1986, 71:457-465
    [31] Dong, Zhengxin; McMechan, George A. 3~D prestack migration in anisotropic media. Geophysics,1993, 58(1): 79-90
    [32]何兵寿.矢量波场叠前逆时深度偏移方法研究.石油大学博士学位论文,2002:29-46
    [33]闫有平,顾春桥,高峻.井中地震勘探技术和仪器的最新进展.物探装备,2003,13(1):1-7
    [34]曹辉.井中地球物理技术综述.勘探地球物理进展,2004, 25(6):6-10
    [35]陆基孟.地震勘探原理.东营:石油大学出版社,1993:92-96
    [36]朱光明.垂直地震剖面方法.北京:石油工业出版社,1988:32-46
    [37] Gal’perin E.I.Vertical seismic profiling.SEG special publication,1974,12:101-110
    [38] Hardage B.A.Vertical seismic profiling part A: principles.1984:241-249
    [39]陈世军,刘洪,周建宇等.井间地震技术的现状与展望.地球物理学进展,2003, 18(3):524-529
    [40]赵军锋,孙红.VSP采集技术应用探讨.石油仪器,.2003.17(6):57-59
    [41]曹辉,郭全仕,唐全良等.井间地震资料特点分析.勘探地球物理进展,2006:29(5):312-318
    [42]张军华,吕宁,田连玉等.地震资料去噪方法综合评述.石油地球物理勘探,2006:40(增刊):121-128
    [43]左建军.垦71井区井间地震资料处理方法.勘探地球物理进展,2007, 30(2):119-123
    [44]严又生,宜明理,魏新等.三维三分量VSP数据处理方法及效果.石油地球物理勘探,2005, 40(1):18-27
    [45]田雪丰,张奋轩,曹新领.弹性波逆时偏移及成像条件.中国煤田地质,2002,14(2):58-61
    [46]戚敬华,李萍.利用? ?p变换技术实现多波波长分离.煤田地震与勘探,1998,26(5):54-57
    [47]孙显义,陈可为,许世勇.采用? ?p变换法进行纵横波波场分离.大庆石油地质与开发,2002, 21(4):76-77
    [48] Carswell A.,Tang R.,Dillistone C.,et al. A new method of Field Separation in VSP Data Processing[J]. Expanded Abstracts of 54th Annual International SEG Mtg,1984: 40-42
    [49] Orhan Yilmaz, M.Thurhan taner.Discrete plane-wave decomposition by least-mean-square-error method.Geophysics,1994, 59(6):973-982
    [50] Carswell A., Tang R.,Dillistone C.,et al. A new method of Field Separation in VSP Data Processing[J]. Expanded Abstracts of 54th Annual International SEG Mtg,1984: 40-42
    [51]陈建平.LU分解递归算法的实现.数值计算与计算机应用,2004, 31(6):141-142
    [52]戴华林.快速高分辨率抛物线拉东变换.硕士学位论文,长安大学西安,2001,6(2):37-51
    [53]陈遵德,段天友等.SVD滤波方法的改进及应用[J] .石油地球物理勘探,1994,29(6):783-792
    [54] J.Mars, N. Thirion. Wave field Separation Based On A Bland Algorithm. SEG Technical Program Expanded Abstracts,1996:1547-1550
    [55]张军华,吕宁,雷凌等.抛物线拉东变换消除多次波的应用因素分析.石油地球物理勘探,2004,39(4):395-406
    [56]巩向博,单刚义,刘冰等.压制截断效应的Radon变换在波场分离中的应用.吉林大学学报(地球科学版),2007, 37(增刊):74-77
    [57]刘喜武,刘洪,李幼铭.高分辨率Radon变换方法及其在地震信号处理中的应用.地球物理学进展,2004, 19(1):8-15
    [58]曾有良,乐有喜,单启铜等.基于高分辨率Radon变换的VSP波场分离方法.石油物探,2007, 46(2):115-120
    [59] Mauricio Dino Sacchi. Aperture Compensated Radon and Fourier Transforms. The University of British Columbia,1996: 53-99
    [60] Douglas J.Foster, James E. Gaiser. Elastic Wave Field Decomposition of offset VSP Data. SEG Technical Program Expanded Abstracts,1986: 563-565
    [61] Trad, Ulrych, Sacchi M D. Accurate interpolation with high resolution time-variant Radon transforms [J].Geophysics,2002, 67(4): 44-656
    [62] Daniel Trade, Tadeusz Ulrych, Mauricio. Sacchi. Latest views of the sparse Radon transform. Geophysics,2002, 68(1): 386-399
    [63] Kurt J. Marfurt, Robert V. Schneider, Michael C. Mueller. Pitfalls of using conventional and discrete Radon transforms on poorly sampled data. Geophysics,1996,61(5): 1467-1482
    [64] Mark Ng, Mike perz. High resolution Radon transform in the t-x domain using intelligent prioritization of the Gauss-Seidel estimation sequence. SEG Int’l Exposition and 74th Annual Meeting Denver, Colorlado,2004: 2160-2163
    [65] Peter W. Cary. The simplest discrete Radon transform. SEG Expanded Abstracts,1998: 1999-2002
    [66]陈建伟.用高分辨率双曲Radon变换实现波场分离.硕士学位论文,长安大学硕士学位论文,西安.2005,6:20-37
    [67]加尔彼林.地震勘探偏振法.石油工业出版社,1984:97-114
    [68]刘财,董世学等.极化滤波在广角地震P?S波场分离中的应用.物探化探计算技术,1995,17(2):15-18
    [69]王海军.应用三分量极化滤波法压制地滚波等椭圆偏振波.西安工程学院硕士论文,1996:21-46
    [70]黄中玉,高林等.三分量数据的偏振分析及其应用.石油物探,1996,35(2):9-16
    [71]李铜基.极化滤波.海洋技术,1996,15(4):130-137
    [72]李英康,崔作舟.分离纵波和横波的偏振旋转法.地球物理学报,1994,37(2):372-381
    [73]张玉芬,周建新.影响空间方向滤波效果的因素分析.石油与天然气地质,1999,20(3):212-215
    [74]胡天跃,张广娟等.多分量地震波波场分解研究.地球物理学报,2004,47(3):504-508
    [75]戚敬华,李萍.利用? ?p变换技术实现多波波长分离.煤田地震与勘探,1998,26(5):54-57
    [76]孙显义,陈可为,许世勇.采用? ?p变换法进行纵横波波场分离.大庆石油地质于开发,2002,21(4):76-77
    [77] Douglas J. Foster. Charles C. Mosher. Up and Downgoing Wave Field Separation for Multi-offset Single Level VSP Data.SEG Technical Program Expanded Abstracts,1988: 816-818
    [78] Xianhuai Zhu, McMechan G.A.. Acoustic modeling and migration of stacked cross-hole data,1988,53(4): 492-500
    [79] Balch A.H.,Chang,H., Hofland,G. S.et al. Use of forward and bac-scattered P and S wave converted waves in cross~borehole imaging. GeoPhysical Prospecting.1991, 39(7): 887-913.
    [80]张昕.井间地震反射波迭前深度偏移方法研究.石油大学硕士学位论文,1996:21-29
    [81]马在田.地震成像技术有限差分偏移[M] .北京:石油工业出版社,1989:61-84
    [82] Stewart A.Levin.Principle of reverse-time migration.GEOPHYSICS,1984, 49(5):581-583
    [83] Dan Loewenthal, Irshad R. Mufti. Reversed time migration in spatial frequency domain, 1983, 48(5):401-409
    [84]傅旦丹,朱宏彰.复杂介质多波地震正演模拟的动力学射线追踪法.中国海上油气(地质),1999,13(5):328-333
    [85] 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 a heavy oil steamflood.Geophysics, 1995, 60(3): 651-659
    [86]赵成斌,许云等.地震波走时和射线的有限差分计算,1995,13(4):41-46
    [87]梁非,杨慧珠等.横观各向同性介质中弹性波的吸收边界条件.固体力学学报,1994,15(4):311-318
    [88] Randall C J. Absorbing boundary conditions for the elastic wave equation. Geophysics,1988,53:611-624
    [89]马在田.高阶近似有限差分偏移.石油地震物理勘探,1982,17(1):5-15

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

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

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