收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

地震反射走时层析理论与应用研究

成谷  
【摘要】: 速度建模是地震勘探中的一项重要技术。其中偏移速度分析和层析成像是两种很重要的手段。从理论上看,速度建模的精度应以层析成像为高,但层析成像在实践过程中受诸多因素的制约,其难点不在于理论方法上,而在于具体实现的技巧上。 本文对两种速度建模方法进行了系统的研究,并重点研究了利用反射地震资料进行走时层析成像速度建模的方法。考虑到层析成像在地震中应用的很多理论和实际上的困难,更由于反射层析中速度和反射深度的耦合性,本文确立了用偏移速度分析的结果为层析成像提供初始速度和反射深度的思想,使得层析成像能够具有一个好的起点。 本文对偏移速度分析的发展现状进行了简单的介绍,详细阐述了层析成像中存在的问题和国内外的发展现状,并针对这两种速度建模方法进行了详细的比较。本文对偏移速度分析的要求不高,偏移速度分析的结果只需给出常速模型即可,本文更对偏移速度分析提供的初始反射深度感兴趣。鉴于速度分析的迭代型和交互性,本文选取了速度较快的局部区域的Kirchhoff偏移方法和沿层进行偏移速度分析的策略。 关于地震层析成像,本文主要介绍了几个方面的内容:层析方程组的数值解法、正则化和归一化。层析方程组的求解,经过了矩阵直接求逆和迭代解两个阶段。矩阵直接求逆的解法,只对小型的矩阵适用。迭代解法主要分为ART类方法和投影类方法。ART类方法是一种行处理方法,是针对早期计算机能力的限制而采取的方法。投影类主要有Lanczos方法、共轭梯度法(CG)、最小二乘共轭梯度法(LSCG)以及最小二乘QR分解法(LSQR)。是目前最经常使用的一类方法。在地震层析中利用正则化对欠定分量和零空间分量进行约束是非常重要的,正则化对欠定分量和零空间分量的约束主要体现为最简单解的思想:包括最平坦解、最光滑解和最小长度解思想,最平坦解思想的出发点是使解的分量间只发生平缓的变化,最光滑解的思想是解的分量间的变化是光滑的,最小长度解思想是使解矢量相对于某个先验矢量发生最小的变化。本文将正则化的具体作用方式分为加法型和乘法型两种。加法型是指将正则化方程组补在原方程组的下面的处理方式,乘法型是指将正则化方程组与原方程组相乘的处理方式。在对正则化的作用范围和作用方式进行介绍后,本文又对正则化的矩阵表述形式、压缩存储进行了阐述,并针对简单模型对导数型正则化的效果进行了分析。归一化是地震层析成像中另一个非常重要的问题。反演计算中,数据d和模型参数m的物理量纲不同,而在多参数类型同时反演时,模型参数和数据中的元本身就带有不同的物理量纲。此时Frechet微商矩阵的列向量的模相差极大,导致问题的性状很坏。因此在进行数据反演计算之前,应该把所有的模型参数和数据进行归一化化为无量纲的量。 针对反射层析中速度和深度的耦合现象,本文综合多人的研究结果得出结论如下:速度深度的耦合性在本质上是由短偏移距引起的,大偏移距的存在能够缓解这种耦合。速度深度的耦合程度受几方面因素的影响:(1)偏移距与深度的比值,偏移距与深度的比值越大,速度深度的耦合性越弱;(2)拾取误差项,拾取误差越小,速度深度的耦合性越弱;(3)界面倾角,倾角越小,速度深度的耦合性越弱:(4)扰动波长与特征波长的比值,当扰动波长大于特征波长时慢度扰动占主要地位,当扰动波长小于特征波长时,深度扰动占主要地位:当扰动波长与特征波长相等时,速度和反射深度发生耦合。特征波长取决于异常体的厚度和异常体与反射面的距离。此外本文还分析了走时误差对反射层析中利用走时扰动分辨薄层内速度扰动的稳定性的影响,并得出结论:当层厚较薄时,由于走时信息中误差的影响,利用走时扰动分辨薄层内的速度扰动是非常不稳定的。层内速度越高、偏移距越小,不稳定性越强。本文针对反射层析的横向和垂向分辨能力也进行了分析,分析证明由于反射地震勘探中以近垂直方向射线为主,因而反射层析的横向分辨能力要优于垂向分辨能力。 本文针对矩形网格参数化的形式建立了层析反演体系,在正演模拟中选用了Langan射线追踪方法,这种方法主要有三个优点:(1)网格内速度的三参数(对二维情况)表述方式在本质上照顾到了射线理论对速度平滑性的要求,方法本身有一种内在的平滑效应;(2)Langan法射线追踪具有完美的解析表达形式,且能够保证精度和效率;(3)在模型矢量选为慢度矢量差,数据矢量选为走时差的情况下,射线在网格内的长度就是反映数据变化对参数变化灵敏度的Frechet微商。考虑到均一矩形网格参数化导致的巨大的矩阵维数给计算机存储量和计算量造成的压力,在对矩形网格参数化层析反演中,本文采用了行索引的压缩存储方法,并以压缩存储的方式补入了一阶导数型正则化和二阶导数型正则化矩阵,最后针对补入正则化矩阵的压缩矩阵采用了LSQR方法迭代求解,大大节省了存储量、计算量和计算时间。 尽管在目前的地震层析成像反演中,采用最多的是矩形网格参数化的方式,但矩形网格参数化方式在很多方面不容乐观:在矩形网格参数化形式下,对地下介质的精细描述必然导致巨大的矩阵维数(各种正则化约束加入后尤甚),方程组的大型稀疏特性使得在求解过程中耗费了大量的计算机存储空间和计算时间。尽管本文针对方程组的稀疏特性设计了一种行索引的压缩存储方法以及针对压缩存储矩阵的迭代解法,大大节省了存储量、计算量和计算时间,但考虑到层析速度建模的迭代性和交互性,存储量和计算时间仍然是矩形网格参数化层析反演应用的一个限制因素。除此之外,矩形网格的均一剖分形式灵活性差,不论在构造简单区域和构造复杂区域都不能区别对待,增大了未知参数的个数,不仅在正演模拟时增加了计算量和计算时间,而且在层析反演时使问题的性态变得更坏,增加了层析反演的难度。矩形网格参数化在对界面的描述上也同样存在着缺陷,即速度模型和界面模型的不一致,由于矩形网格参数化只能用锯齿状的分布来表征界面上下的速度变化,在对界面的处理上增加了难度和引入了误差。针对矩形网格参数化方式在剖分上的不灵活给正反演方法带来的诸多不利因素,考虑到三角网格剖分的灵活性和本质上的正则化效应,本文进一步建立了可变三角网格参数化的层析反演体系,得到了很大的改善。三角网格参数化层析反演相对于矩形网格参数化层析反演无论在计算时间、计算量,还是在存储量上都占有非常大的优势,经验证两种参数化方式得到的层析结果具有可比性。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 张东;童海杰;杨艳;秦前清;徐凌;;基于波前扩展的线性走时插值射线追踪算法[J];武汉大学学报(理学版);2010年06期
2 王云安;艾东海;;基于电磁波CT的岩溶地质异常探测[J];人民黄河;2009年12期
3 黄靓,黄政宇;线性插值射线追踪的改进方法[J];湘潭大学自然科学学报;2002年04期
4 赵连锋,朱介寿,曹俊兴;并行化交错网格法地震层析成像[J];石油物探;2003年01期
5 杨慧珠,聂建新,谭桂华;反射波地震层析成像算法研究[J];辽宁工程技术大学学报;2004年02期
6 王祯祥,高尔根,徐果明,刘同庆;一种层状介质条件下射线的全路径迭代追踪法[J];煤田地质与勘探;2002年04期
7 ;Geophysics,No.1,1992论文文摘[J];勘探地球物理进展;1992年05期
8 许永忠,柴振友,宋端智;煤矿地应力异常层析成像方法的研究[J];中国煤田地质;2000年01期
9 周兵,朱介寿,S.A.Greenhalgh,C.Sinadinovski;一种新的地震射线层析成像计算方法[J];石油物探;1994年01期
10 林伯香,孙晶梅,刘清林;层析成像低速带速度反演和静校正方法[J];石油物探;2002年02期
11 宋林平;平方慢度射线追踪与非线性地震走时层析成像[J];成都理工学院学报;1994年01期
12 郑需要,张先康;地壳和上地幔三维速度与界面结构层析成像——理论部分[J];地震学报;1998年05期
13 张树林;朱介寿;贺振华;;井间地震弯曲射线层析成像[J];石油地球物理勘探;1992年02期
14 鲁传恒,顾贤明;井间地震射线层析成像[J];CT理论与应用研究;1994年03期
15 薛广盈,丁韫玉,袁志祥;渭河断陷盆地地壳速度的层析成像研究[J];地震学报;1997年03期
16 熊高君;肖又军;彭更新;贺振华;王绪本;;射线层析与波动方程联合成像[J];新疆石油地质;2010年06期
17 陈湛文,张关泉;用γ道集信息修正偏移速度模型的层析成像法[J];地球物理学进展;2004年01期
18 张霖斌,姚振兴;阻尼同时代数重建技术的应用[J];物探与化探;1997年02期
19 孔庆丰;左建军;魏国华;何惺华;;井间地震资料处理方法研究与应用[J];物探与化探;2006年06期
20 廖武林;姚运生;丁志峰;张丽芬;;三峡地区P波速度层析成像研究[J];大地测量与地球动力学;2007年03期
中国重要会议论文全文数据库 前10条
1 周连庆;赵翠萍;陈章立;;利用地方震层析成像研究龙滩水库库区三维V_P和V_P/V_S层析成像[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
2 孙晶梅;李博;林伯香;徐颖;;复杂近地表速度分析及校正[A];油气地球物理实用新技术——中国石化石油勘探开发研究院南京石油物探研究所2005年学术交流会论文集[C];2005年
3 周连庆;赵翠萍;郑现;陈章立;;三峡库区三维Vp,Vp/Vs,Qp和Qs精细结构层析成像[A];中国地球物理学会第二十七届年会论文集[C];2011年
4 李玉喜;刘四新;;跨孔层析成像四种算法的研究[A];中国地球物理·2009[C];2009年
5 曾吕明;邢达;谷怀民;杨迪武;杨思华;向良忠;;基于多元线性阵列探测器的快速微波致热声层析成像[A];第十次中国生物物理学术大会论文摘要集[C];2006年
6 裴顺平;许忠淮;汪素云;;中国大陆Sn波速度层析成像[A];中国地球物理学会年刊2002——中国地球物理学会第十八届年会论文集[C];2002年
7 高丽丽;何琳;徐果明;朱良保;;层析成像中反演参数选取的影响[A];中国地球物理学会第二十届年会论文集[C];2004年
8 房立华;吴建平;王未来;明跃红;;华北地区噪声层析成像研究[A];中国地震学会成立三十年学术研讨会论文摘要集[C];2009年
9 张宝龙;徐义贤;陈伟;蒋礼;;基于噪声的面波层析成像与观测实验[A];中国地球物理学会第二十七届年会论文集[C];2011年
10 肖化;胡广莉;何惠玲;保宗悌;;基于图象分割的人工神经网络在两相流层析成像中应用研究[A];中国图象图形科学技术新进展——第九届全国图象图形科技大会论文集[C];1998年
中国博士学位论文全文数据库 前10条
1 成谷;地震反射走时层析理论与应用研究[D];同济大学;2004年
2 王浩全;超声阵列层析成像检测方法研究[D];中北大学;2011年
3 马辉;基于MapReduce的分布式地震射线追踪方法研究[D];中国地质大学(北京);2012年
4 袁志亮;井间声波电磁波层析成像技术应用研究与软件研发[D];中国地质大学(北京);2007年
5 郝奇;VTI介质速度和各向异性参数建模研究[D];吉林大学;2010年
6 黄靓;混凝土超声波层析成像的理论方法和试验研究[D];湖南大学;2008年
7 喻振华;工程井间地震波正演模拟及层析成像技术研究[D];中南大学;2008年
8 白志明;深地震测深剖面的层析成像研究及其应用[D];中国地震局地球物理研究所;2002年
9 王振宇;土木工程的层析成像与广义反演研究[D];浙江大学;2003年
10 王兆国;东北地区地震综合信息联合成像与震源机制反演[D];吉林大学;2010年
中国硕士学位论文全文数据库 前10条
1 李国磊;基于大炮初至数据的层析成像算法研究[D];中国石油大学;2010年
2 郭彦刚;井间地震层析成像技术及应用研究[D];成都理工大学;2012年
3 周小仙;利用地震反射和折射联合层析成像探测深部构造方法研究[D];中国地质大学(北京);2010年
4 张赛民;跨孔地震层析成像研究[D];中南大学;2003年
5 燕利芳;反射波层析成像的数值模拟及在西安地裂缝中的应用研究[D];长安大学;2011年
6 李志辉;CT技术在路基病害勘查治理中的研究及应用[D];西南交通大学;2005年
7 金溪;浅层地震初至波层析成像方法技术研究[D];长安大学;2010年
8 罗文歆;井间地震层析成像研究及其在隐伏岩溶勘察中的应用[D];中南大学;2011年
9 黄靓;混凝土超声CT的数值模拟与试验研究[D];湖南大学;2003年
10 戚秀真;混凝土超声无损检测层析成像技术研究[D];长安大学;2006年
中国重要报纸全文数据库 前10条
1 华为技术有限公司 卢鹏;射线追踪技术为网络规划导航[N];通信产业报;2006年
2 黄寰;层析成像技术揭示地球奥秘[N];中国矿业报;2004年
3 本报记者 孟凡君;资源勘探的得力助手[N];中国矿业报;2008年
4 王五平;声波CT为混凝土大坝做体检[N];中国水利报;2005年
5 本报特约记者 曹菲;杨文采与他的应用地球物理反演理论[N];地质勘查导报;2006年
6 记者 刘如春 通讯员 刘文祥;地震CT技术洞悉岩溶变化[N];地质勘查导报;2008年
7 刘洪宇;英特尔Larrabee渲染架构[N];中国计算机报;2008年
8 石岩;中天山地幔精细结构研究项目[N];科技日报;2007年
9 强声;美国研制能透视救灾的无线网络系统[N];人民公安报·消防周刊;2010年
10 ;探索地球科学服务人类社会[N];科技日报;2003年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978