用户名: 密码: 验证码:
复杂地表模型菲涅尔体层析反演静校正方法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
在各种油气勘探方法中,地震勘探一直是最重要的技术方法之一,并在现代石油勘探中发挥着重要作用。我国在平原地区的油气勘探开发已经比较深入,目前陆地地震勘探的主战场大多已转向我国西北、西南等地区。这些地区大多以山地、沙漠、黄土塬等复杂地表条件为主。静校正、叠前去噪、偏移成像等是此类复杂地区地震资料处理的重点和难点问题。尤其是静校正是地震资料处理中首先要完成的工序,静校正处理效果的好坏不仅直接影响后续的处理工序,而且有可能影响最终的成像效果。
     本文针对复杂地表模型下静校正技术的重点和难点,深入研究了初至波规律及利用初至波进行层析反演得到速度模型的方法,并应用于实际资料的层析反演与静校正,取得较好效果。
     论文首先研究分析了国内外层析反演的各种理论、方法与技术,然后对二、三维射线追踪算法进行了研究,改进了旅行时线性插值算法,提高了射线走时正演的精度;根据菲涅尔体的定义提出并实现了快速计算菲涅尔体的筛选法;然后改进反投影算法,并将其用于基于菲涅尔体的层析反演,在保证精度的情况下具有较高的计算效率;通过二维、三维几种典型理论模型的试算验证了本文提出的正演、反演算法的正确性与有效性,并应用于实际资料处理,取得较好的应用效果。
     综上所述,本文实现了基于菲涅尔体的复杂地表速度模型的正反演,能完成复杂地表地震资料的层析反演成像与静校正。本文方法有良好的应用前景。
In all kinds of oil and gas exploration method, seismic exploration has been one of the most important techniques in modern petroleum exploration, and plays an important role. As it has been deeply used in the plain area, seismic exploration for oil and gas exploration and development has shifted it’s main battlefield to northwest and southwest areas, where the main surface are mountain,desert and loess plateau. Static correction, pre-stack denoise and migration imaging are key and most difficult seismic data processing problems to be solved in these complicated surface areas. Static correction is especially the first and most important seismic data processing task because the effect of static correction will not only affect the following processing task and may affect the final imaging.
     Aiming at the main and deficullt problems at these complex surface model, this paper studied the theory of first arrival wave and used to tomographic inversion model reversing method, and applied to the practical data tomographic inversion and static correction, obtained good effect.
     Based on the research of the analysis of tomographic inversion theory, method and technology from home and abroad, the paper has studied the traditional ray tracing methods of 2D and 3D space, and has improved LTI method to increase the precition of ray tracing. Proposed and implemented a rapid sieve method to calculate Fresnel volume, then improved the BPT method and applied it to the tomographic inversion based on Fresnel volume, got a high computational efficiency while ensuring the accuracy. The forward and inverse methods proposed by the paper have been proved to be correct and effective by the typical 2D and 3D theoretical models, and have got some good effect when applying to actual seismic data.
     In summary, this paper realized the forwarding and inversion methods of the velocity model based on Fresnel volume, and is able to perform tomographic inversion and static correction to complex surface seismic data. The methods have good application prospect.
引文
[1] COX M. Static Corrections for Seismic Reflection Surveys:反向地震勘探静校正技术[M].柯本喜,李培明,译.石油工业出版社, 2004.
    [2]顾庆雷.近地表模型约束层析静校正技术研究与应用[D].中国海洋大学地质工程, 2009.
    [3]陈世军,张建中.初至波射线层析成像在复杂区静校正中的应用[J].石油物探, 2006(01):454-457.
    [4]冯泽元,李培明,唐海忠,等.利用层析反演技术解决山地复杂区静校正问题[J].石油物探, 2005(03):636-641.
    [5]韩晓丽,杨长春,麻三怀,等.复杂山区初至波层析反演静校正[J].地球物理学进展, 2008(02):475-483.
    [6]井西利.复杂介质的层析建模及模型反演[J].石油地球物理勘探, 2010(01):66-71.
    [7]陈爱萍,梁波,邹文,等.初至波地震层析技术及其在四川复杂地区的应用[J].世界地质, 2006(04):440-444.
    [8]李晓斌,常锁亮,杨晓东.复杂地表条件剖分折射法静校正处理方法探索[J].中国煤田地质, 2006(02):54-56.
    [9]王振华,袁明生,阎玉魁,等.复杂地表条件下的静校正方法[J].石油地球物理勘探, 2003(05):101-106.
    [10]酆少英,何瑞珠,王巍,等.复杂地表条件下的静校正方法研究[J].河南石油, 2002(03):73-81.
    [11]段洪有,曾庆才,李琛.复杂地表条件下折射静校正技术的应用[J].物探与化探, 2005(02):142-145.
    [12]宋常洲,高建华.复杂地区静校正方法研究[J].石油地球物理勘探, 2007:32-35.
    [13]李红.复杂地区静校正技术应用研究[J].安徽地质, 2007(01):598-602.
    [14]王克斌.复杂地表条件下初至折射波静校正方法研究[D].成都理工大学地球探测与信息技术, 2004.
    [15]邹强.山地静校正若干问题研究[D].成都理工大学地球探测与信息技术, 2004.
    [16]成谷.地震反射走时层析理论与应用研究[D].同济大学固体地球物理学, 2004.
    [17]于富文.中国西部地区静校正方法对比研究[D].中国地质大学(北京)地球探测与信息技术, 2006.
    [18]高磊.地震初至波自动拾取和折射波表层调查方法研究[D].成都理工大学地球探测与信息技术, 2008.
    [19]何晓冬,夏忠谋,张国保.模型法静校正在西部地区的应用,中国安徽合肥, 2009[C].
    [20]夏忠谋,何晓冬.模型法静校正在准噶尔盆地腹部地区的应用[J].中国西部油气地质, 2005(2):203-206.
    [21]罗英伟,段卫星,徐维秀,等.几种静校正方法的研究与比较[J].油气地球物理, 2010(1):34-39.
    [22]华勇.折射层析静校正方法及其应用效果[J].勘探地球物理进展, 2002(5):52-54.
    [23] Taner M. Turban, Matsuoka Toshifumi, Baysal Edip,等.折射地震波成象,美国新奥尔良, 1992[C].
    [24]赵德亨.折射波广义互换法在城市地质调查中的应用[D].吉林大学, 2005.
    [25]华勇.折射层析静校正方法及其应用效果[J].勘探地球物理进展, 2002(5):52-54.
    [26]刘洪雷.复杂地区的折射波静校正应用研究[D].中国地质大学(北京), 2006.
    [27]王克斌,王顺根.利用扩展广义互换折射波静校正方法解决MX地区资料的野外静校正闭合差[J].石油物探, 2001(2):126-130.
    [28]李颜贵.山西复杂地表区煤田地震勘探折射静校正方法应用研究[D].太原理工大学, 2006.
    [29]皮金云,段云卿,康平.综合折射静校正方法在柴达木盆地的应用[J].地球科学(中国地质大学学报), 2007(4):487-490.
    [30]罗有春,雷宛,王怀坤,等.层析反演静校正技术及其应用[J].内蒙古石油化工, 2009(9):72-74.
    [31]张正峰,韩立国.变换坐标系下的转换波波动方程延拓静校正,中国安徽合肥, 2009[C].
    [32]王守东.复杂地表波动方程反演延拓静校正[J].石油地球物理勘探, 2005(1):31-34.
    [33]崔兴福,徐凌,陈立康.复杂近地表波动方程波场延拓静校正[J].石油勘探与开发, 2006(1):80-82.
    [34]方伍宝,李满树,孙爱萍,等.基于Born近似的波动方程静校正技术[J].石油物探, 2004(1):26-29.
    [35]刘素芹.基于应用网格环境的复杂地表波动方程基准面静校正研究[D].中国石油大学, 2008.
    [36]李玮.地震信号差分法初至折射波自动剩余静校正[D].西安理工大学, 2008.
    [37]李录明.动校前初至波自动剩余静校正[J].石油地球物理勘探, 1989(2):179-186.
    [38]高磊,潘树林,周熙襄,等.利用转换波初至进行转换波剩余静校正[J].大庆石油地质与开发, 2008(1):130-133.
    [39]刘连升.约束初至拾取与初至波剩余静校正[J].石油地球物理勘探, 1998(5):604-610.
    [40]李继光.不同复杂近地表校正技术分析与应用[J].石油天然气学报, 2010(1):223-227.
    [41]张继国,刘连升.复杂区初至层析反演静校正[J].石油地球物理勘探, 2006(4):383-385.
    [42]顾庆雷.近地表模型约束层析静校正技术研究与应用[D].中国海洋大学, 2009.
    [43]李振春,王希萍,韩文功.地震数据处理中的相位校正技术综述[J].地球物理学进展, 2008(3):768-774.
    [44]唐建侯,张金山.高效模拟退火剩余静校正[J].石油地球物理勘探, 1994(3):382-387.
    [45]潘文勇.基于改进模拟退火算法的剩余静校正及程序实现[J].物探与化探, 2010(4):528-531.
    [46]田连玉.基于模拟退火和遗传算法的剩余静校正方法研究[D].中国石油大学, 2007.
    [47]井西利.模拟退火方法的改进及其在地震剩余静校正中的应用[J].西安石油学院学报(自然科学版), 2003(4):19-22.
    [48]郑鸿明,彭立,李生杰.模拟退火静校正[J].新疆石油地质, 2001(1):32-34.
    [49]曲凤山,黄刘生,同利文,等.模拟退火剩余静校算法(ARSCSA)及其并行实现[J].计算机工程, 2001(5):119-122.
    [50]秦亚玲.沙漠地区地震资料处理技术研究[D].中国地质大学(北京), 2006.
    [51] Normark E.剩余静校正估算:模拟退火温度状态标定法的应用[J].石油物探译丛, 1994(4):76.
    [52]于富文.中国西部地区静校正方法对比研究[D].中国地质大学(北京), 2006.
    [53]马德堂.弹性波场数值模拟及井间地震初至波旅行时层析成像[D].长安大学地质工程, 2006.
    [54]罗有春,雷宛,王怀坤,等.层析反演静校正技术及其应用[J].内蒙古石油化工, 2009(9):72-74.
    [55]陈会芬.层析技术在静校正中的应用[J].油气地球物理, 2003(2):40-43.
    [56]张树林,朱介寿,贺振华.井间地震弯曲射线层析成像[J].石油地球物理勘探, 1992(2):270-278.
    [57]王宏伟.菲涅耳体地震层析[J].石油地球物理勘探, 1996(03):309-326.
    [58]陆基孟. [M].石油大学出版社, 1993.
    [59]陈仲侯,付唯一.浅层地震勘探[M].成都:成都地质学院, 1986.
    [60]杨文采.地球物理反演和地震层析成像[M].北京:地质出版社, 1989.
    [61]徐涛,徐果明,高尔根,等.三维复杂介质的块状建模和试射射线追踪[J].地球物理学报, 2004(06).
    [62]徐涛,徐果明,高尔根,等.三维试射射线追踪子三角形法[J].石油地球物理勘探, 2005(04):391-399.
    [63]宋葵著蒋选艺刘贤功.复杂构造模型正演模拟[M]. 2004.
    [64] Koketsu K., Sekine S.,雷建设,等.具有间断面的球形地球三维地震射线追踪的伪弯曲法[J].世界地震译丛, 2000(05):50-58.
    [65]陈会芬.层析技术在静校正中的应用[J].油气地球物理, 2003(02):397-399.
    [66]陈国金,曹辉,吴永栓,等.最短路径层析成像技术在井间地震中的应用[J].石油物探, 2004(04):327-330.
    [67]卞爱飞,於文辉.三维最短路径法射线追踪及改进[J].天然气工业, 2006(05):43-45.
    [68]段心标,金维浚.井间层析成像混合最短射线追踪旅行时反演[J].工程地球物理学报, 2007(03):226-230.
    [69]杨薇.基于弯曲射线的跨孔层析成像算法研究[D].吉林大学地球探测与信息技术, 2008.
    [70]王兆国.东北地区地震综合信息联合成像与震源机制反演[D].吉林大学固体地球物理学, 2010.
    [71]梅胜全.三维表层速度结构层析反演方法研究[D].成都:成都理工大学信息工程学院, 2010.
    [72] I N, K Y. A numerical experiment on nonlinear image reconsruction from first-arrival times for two-dimensional island ard structure[J]. J. Phys, Earth, 1986(34):195-201.
    [73] MOSER. Shortest path calculation of seismic rays[J]. Geophysics, 1991,56(1):59-67.
    [74] L K, M K. 3-D network ray tracing[J]. Geophys, 1994,116(4):726-738.
    [75] Asakawa Eiichi,蒋录全.用于地震射线追踪的旅行时线性插值法[J].石油物探译丛, 1993(04):673-677.
    [76] ASAKAWA, EIICHI, KAWANAKA. Seismic ray tracing using linear traveltime interpolation[J]. Geophysics, 1993,5(2):326-333.
    [77] Asakawa Eiichi,石河.利用线性旅行时内插的地震射线追踪方法[J].石油物探译丛, 1994(01):107-113.
    [78] Asakawa Eiichi, Kawanaka Taku,严又生.旅行时线性内插地震射线追踪[J].国外油气勘探, 1995(01):59-67.
    [79]赵改善,郝守玲,杨尔皓,等.基于旅行时线性插值的地震射线追踪算法[J].石油物探, 1998(02):14-24.
    [80]聂建新,杨慧珠.地震波旅行时二次/线性联合插值法[J].清华大学学报(自然科学版), 2003(11):1495-1498.
    [81]彭直兴,张芬,周熙襄,等.地震波旅行时非线性/线性联合插值法[J].成都理工大学学报(自然科学版), 2005(03):322-324.
    [82]彭直兴,张春红,周熙襄,等.基于非线性插值的地震波旅行时计算[J].物探化探计算技术, 2006(01):10-13.
    [83] VIDALE, E J. Finite difference calculation of traveltime[J]. Bull. Seis. Soc. Am, 1988(78):2062-2076.
    [84] PODVIN P, LAECOMTE. Finite difference computation of travel-times in very contrasted velocity model: a massively parallel approach and its associated tools[J]. Geophys J. Int, 1991,105:271-284.
    [85] F Q, K. O, Y. L, et al. Finite-difference solution of the eikonal equation along expanding wavefronts[J]. Geophysics, 1992(57):478-487.
    [86]张霖斌,许云,乌达巴拉,等.地震反射波波前法射线追踪[J].石油地球物理勘探, 1995(01):50-55.
    [87]张霖斌,姚振兴,纪晨.地震初至波走时的有限差分计算[J].地球物理学进展, 1996(04):47-52.
    [88]赵连锋,朱介寿,曹俊兴,等.有序波前重建法的射线追踪[J].地球物理学报, 2003(03).
    [89] A S J, M. P. Three dimensional traveltimes computation using the Fast Marching method[J]. Geophysics, 1999,64(2):516-523.
    [90] R K, A. S J. A Fast marching methods on triangulated domains[J]. Proc Nat Acad Sci, 1998(95):8341-8435.
    [91]黄联捷,李幼铭,吴如山.用于图像重建的波前法射线追踪[J].地球物理学报, 1992(02).
    [92]李录明,罗省贤.复杂三维表层模型层析反演与静校正[J].石油地球物理勘探, 2003,38(6):636-641.
    [93]李录明,罗省贤.三维波场延拓复杂表层模型校正方法及应用[J].地球物理学进展, 2005,20(4):1024-1034.
    [94]罗省贤,李录明.地面地震初至波层析反演复杂表层速度结构方法[J].地球科学进展, 2004,19:30-34.
    [95] P S, S. F. Huygens wavefront tracing: A robust alternative to ray tracing[J]. Geophys.,Expanded Abstracts, 1998(29):1961-1964.
    [96] HAGEDOORN J G. A process of seismic reflection interpretation[J]. Geophysisical Prospecting, 1954(2):85-127.
    [97] STEPHAN HUSEN E K. Local earthquake tomography between rays and waves:fat ray tomography[J]. Physics of the earth and Planetary, 2001(123):129-149.
    [98] A H F K E. 3D crustal structure fromlocal earthquake tomography around the Gulf of Arta(Ionianregion,NW Greece)[J]. Tectomophysics, 1999(304):201-218.
    [99] V S, J D J. Seismic tomography with local earthquakes in Costa Rica[J]. Tectonophysics, 2000(304):201-218.
    [100] F H, E. K. Investigating efects of 3D ray tracing methods in local earthquake tomography[J]. Physics of the Earth and Planetary Interiors,123:103-114.
    [101] AL M R G N. Finite-Frequency Tomography Reveals a Variety of Plumes in the Mantle [J]. 2004,303:338-343.
    [102]朱丹.波前快速推进法三维走时计算技术[J].大庆石油地质与开发, 2001,20(5):65-66.
    [103]张建中,杨国辉,林文,等. Fresnel层析成像并行算法研究[J].计算机研究与发展, 2007(10):1661-1666.
    [104]朱生旺.近地表速度菲涅耳带初至层析反演[J].新疆石油地质, 2007(06):752-755.
    [105]刘玉柱,董良国,李培明.初至波菲涅尔体地震层析成像[J].地球物理学报, 2009,52(9):231-2320.
    [106]赵妍.层析成像静校正方法研究[D].中国石油大学地球探测与信息技术, 2009.
    [107]陈誌敏.龙格-库塔法及其Mathematica实现[J].武汉工程职业技术学院学报, 2006(2):72-76.
    [108]李夏云,陈传淼.用龙格-库塔法求解非线性方程组[J].数学理论与应用, 2008(2):62-65.
    [109]吴律.层析基础及其在井间地震中的应用[M].北京:石油工业出版社, 1997.
    [110] VIDALE J. Finite-difference calculation of traveltime in three dimensions[J]. Geophysics, 1990(55):521-526.
    [111] DIJKSTRA E W. A note on tow problems in connection with graphs[J]. Number Math, 1959(1):269-271.
    [112]喻维秋.初至波表层模型层析反演及其应用[D].成都理工大学固体地球物理学, 2007.
    [113] FICHER R, LEES J M. Shortest path ray tracing with sparse graphs[J]. Geophysics, 1993,58(7):987-996.
    [114] VAN AVENDONK, HARDING A J. Shortest path and ray bending method for traveltime and raypath calculations[J]. Geophysics, 2001,66(2):648-653.
    [115]张建中,陈世军,余大祥.最短路径射线追踪方法及其改进[J].地球物理学进展, 2003(01):146-150.
    [116]张建中,陈世军,徐初伟.动态网络最短路径射线追踪[J].油气地球物理, 2003(04):13-18.
    [117]张建中,陈世军,徐初伟.动态网络最短路径射线追踪[J].地球物理学报, 2004(05):900-905.
    [118]杨慧珠,聂建新,谭桂华.反射波地震层析成像算法研究[J].辽宁工程技术大学学报, 2004(02):159-161.
    [119]彭直兴.地震波初至旅行时射线追踪正演方法研究[D].成都理工大学地球探测与信息技术, 2005.
    [120] MCMECHAN M A. Seismic tomography in boreholes[J]. Geophysics, 1983,74(2):326-333.
    [121]刘伊克著常旭.地震正反演与成像[M]. 2001.
    [122]姜国博.初至波表层模型层析正反演及层析静校正研究[D].中国地质大学(北京)地球探测与信息技术, 2007.
    [123]吕子强.初至波地震层析成像研究[D].中国地震局兰州地震研究所固体地球物理, 2009.
    [124] R W P, H W M. Resolution limits in ray tomograph due to wave behavior-numerical experiments[J]. Geophysics, 1993,58(5):727-735.
    [125] NOLET G. Seismic tomography with applications in global seismology and exploration geophysics[J]. D.Reidel Publishing Co., 1987.
    [126] YOMOGIDA K. Fresnel zone inversion for heterogeneities in the earth[J]. Pageoph, 1992,138(3):391-406.
    [127]俞寿朋.高分辨率地震勘探[M].北京:石油工业出版社, 1993.
    [128]杨文采.应用地震层析成像[M].北京:地质出版社, 1993.
    [129] H Z, A M J, A L C, et al. Constraining the magnitude of velocity perturbations in traveltime tomography[J]. J Seism Exploration, 1993(2):365-380.
    [130] K A, W L. Determination of three-dimensional anomalies under a seismic array using first P arrival times from local earthquakes.[J]. J Geophys Res, 1976,81:4381-4399.
    [131] W C R, R C P. A tomographic analysis of mantle heterogeneities from body wavetravel-times[J]. EOS trans, 1983(62):776.
    [132] M D A. Mapping the lower mantle-determination of lateral heterogeneity in P velocity up to degree and order 6[J]. J Geophys Res, 1984(89):5929-5952.
    [133] JINSONG L, FUTIAN L, JUN L. the optimization and parallelization of sequential orthogonal triangularization algorithms used in seismic tomography[J]. Process in Geophysics, 2005,20(5):911-915.
    [134] V C, P S E. Fresnel volume ray tracing[J]. Geophysics, 1992,57(7):902-915.
    [135]杨国辉.菲涅耳体旅行时层析成像方法及应用研究[D].厦门大学信号与信息处理, 2009.
    [136]刘清林.初始模型和约束条件对层析反演结果的影响[J].石油物探, 1996(02):41-46.
    [137]常旭,卢孟夏,刘伊克.地震层析成像反演中3种广义解的误差分析与评价[J].地球物理学报, 1999(05):695-701.
    [138]刘玉柱,董良国.初至波层析影响因素分析[J].石油地球物理勘探, 2007(05):544-553.
    [139]景月红.地震初至波走时层析成像与近地表速度建模[D].长安大学地球探测与信息技术, 2009.
    [140]伊尔马兹渥.地震数据处理[M].北京:石油工业出版社, 1994.

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

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

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