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

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

成谷  
【摘要】: 速度建模是地震勘探中的一项重要技术。其中偏移速度分析和层析成像是两种很重要的手段。从理论上看,速度建模的精度应以层析成像为高,但层析成像在实践过程中受诸多因素的制约,其难点不在于理论方法上,而在于具体实现的技巧上。 本文对两种速度建模方法进行了系统的研究,并重点研究了利用反射地震资料进行走时层析成像速度建模的方法。考虑到层析成像在地震中应用的很多理论和实际上的困难,更由于反射层析中速度和反射深度的耦合性,本文确立了用偏移速度分析的结果为层析成像提供初始速度和反射深度的思想,使得层析成像能够具有一个好的起点。 本文对偏移速度分析的发展现状进行了简单的介绍,详细阐述了层析成像中存在的问题和国内外的发展现状,并针对这两种速度建模方法进行了详细的比较。本文对偏移速度分析的要求不高,偏移速度分析的结果只需给出常速模型即可,本文更对偏移速度分析提供的初始反射深度感兴趣。鉴于速度分析的迭代型和交互性,本文选取了速度较快的局部区域的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];石油物探;2000年01期
2 周辉,韩波,邱东玲,裴建新,王兆磊;地质雷达资料的偏移速度分析和叠前偏移[J];吉林大学学报(地球科学版);2005年02期
3 H.J.Tieman;偏移速度分析:考虑横向速度变化的影响[J];勘探地球物理进展;1995年05期
4 Gerhard P.Dohr;P.K.Stiller;吴晖;;偏移速度分析—应用[J];石油物探;1975年04期
5 邵幼英;;符号位统计检测偏移速度分析[J];石油地球物理勘探;1980年03期
6 Johann W.Sattlegger;吴晖;;偏移速度分析—原理[J];石油物探;1975年04期
7 PaulDocherty,周英;利用遗传算法的偏移速度分析[J];海洋石油;1999年01期
8 成谷;张宝金;;改进的RCA偏移速度分析相干谱计算方法[J];中山大学学报(自然科学版);2007年S2期
9 李录明,罗省贤;相移加校正叠前深度偏移及偏移速度分析[J];石油地球物理勘探;1997年01期
10 张敏;李振春;;偏移速度分析与建模方法综述[J];勘探地球物理进展;2007年06期
11 李导;张恩嘉;熊晓军;;基于G-D波场分解的角度域共成像点道集提取及成像[J];物探化探计算技术;2010年03期
12 邓富求;;偏移速度分析的应用效果[J];石油地球物理勘探;1982年02期
13 刘守伟;王华忠;程玖兵;马在田;;时空移动成像条件及偏移速度分析[J];地球物理学报;2008年06期
14 杨淑卿;王华忠;程玖兵;马在田;;波动方程偏移角度域共成像道集[J];油气地球物理;2004年02期
15 孙成龙,方伍宝,孔祥宁,王于静,段文超;深度-射线域偏移速度分析[J];石油物探;2005年05期
16 潘宏勋;方伍宝;;地震速度分析方法综述[J];勘探地球物理进展;2006年05期
17 周洪波;;最大熵方法在偏移速度分析中的推广及应用[J];石油物探;2007年06期
18 杨淑卿;王华忠;程玖兵;马在田;胡中标;陈建平;;共方位角波动方程偏移速度分析方法及其应用[J];油气地球物理;2003年03期
19 杨淑卿;王华忠;程玖兵;马在田;孟祥宾;胡中标;;基于地质层位约束的2D波动方程偏移速度分析方法[J];油气地球物理;2003年04期
20 刘守伟;程玖兵;王华忠;马在田;;偏移距域/角度域共成像点道集与偏移速度的关系[J];地球科学(中国地质大学学报);2007年04期
中国重要会议论文全文数据库 前10条
1 胡莲莲;王赟;芦俊;;弹性矢量波场偏移速度分析[A];中国地球物理学会第二十七届年会论文集[C];2011年
2 杨祥森;;面向剩余偏移速度分析的三维叠前深度偏移平面波合成算子[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
3 胡先岷;贺振华;黄德济;;提高共成像道集偏移速度分析的精度和计算速度[A];中国地球物理学会年刊2002——中国地球物理学会第十八届年会论文集[C];2002年
4 孙成龙;刘守伟;王华忠;方伍宝;;时移角度域偏移速度分析[A];中国地球物理学会第二十三届年会论文集[C];2007年
5 薛春玲;;叠前深度偏移速度建模技术[A];中国地球物理学会第二十四届年会论文集[C];2008年
6 巩向博;韩立国;;偏移成像中速度敏感度分析[A];中国地球物理·2009[C];2009年
7 张洪宙;杨长春;冷传波;;建立精细速度-深度模型的综合方法[A];寸丹集——庆贺刘光鼎院士工作50周年学术论文集[C];1998年
8 岳玉波;李振春;;真振幅高斯波束共角度偏移[A];中国地球物理·2009[C];2009年
9 ;工程地质与水资源研究室SCI论文首页[A];中国科学院地质与地球物理研究所2009年度学术论文汇编[C];2010年
10 杨长春;徐升;王德明;;一种共炮集偏移的走时层析成象方法[A];1994年中国地球物理学会第十届学术年会论文集[C];1994年
中国博士学位论文全文数据库 前4条
1 刘守伟;角度域共成像点道集及偏移速度分析[D];同济大学;2007年
2 陈辉;复杂介质波动方程叠前深度偏移及偏移速度分析[D];成都理工大学;2009年
3 成谷;地震反射走时层析理论与应用研究[D];同济大学;2004年
4 巩向博;金属矿地震高精度成像与数据处理方法研究[D];吉林大学;2011年
中国硕士学位论文全文数据库 前10条
1 李芳;基于弹性Kirchhoff偏移的纵横波偏移速度分析方法[D];中国石油大学;2011年
2 韩波;遗传算法在地质雷达成像中的应用研究[D];中国海洋大学;2005年
3 周振兴;偏移速度分析方法研究[D];大庆石油学院;2006年
4 侯波;弹性矢量波场偏移速度分析方法研究[D];中国石油大学;2009年
5 张春堂;分频叠前深度偏移处理方法研究[D];大庆石油学院;2007年
6 吴永国;复杂地质体波动方程地震波场模拟与偏移成像[D];成都理工大学;2007年
7 邓盾;井间地震纵横波速度分析方法研究[D];中国石油大学(华东);2012年
8 曾同生;基于角道集的层析速度反演方法研究[D];中国石油大学(华东);2012年
9 王鹏燕;基于共成像点道集的层析速度反演[D];中国石油大学;2010年
10 徐磊;地震转换波速度交互分析软件的研究与实现[D];电子科技大学;2013年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978