收藏本站
《同济大学》 2004年
收藏 | 手机打开
二维码
手机客户端打开本文

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

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

知网文化
【引证文献】
中国期刊全文数据库 前3条
1 孟宪海;金颖;李吉刚;谭文磊;杨钦;;基于三角域快速行进法的地震波走时计算[J];软件;2011年11期
2 秦宁;李振春;杨晓东;;基于角道集的井约束层析速度反演[J];石油地球物理勘探;2011年05期
3 成谷;张宝金;;基于参数模型两种剖分方式的射线追踪比较[J];中山大学学报(自然科学版);2008年05期
中国重要会议论文全文数据库 前1条
1 黄光南;刘洋;;层析成像数据和射线不均匀覆盖问题的解决方法[A];中国地球物理学会第二十七届年会论文集[C];2011年
中国博士学位论文全文数据库 前1条
1 李军;复杂地表模型菲涅尔体层析反演静校正方法研究[D];成都理工大学;2011年
中国硕士学位论文全文数据库 前4条
1 孟庆利;基于大炮初至的层析静校正算法研究与实现[D];中国石油大学;2011年
2 李国磊;基于大炮初至数据的层析成像算法研究[D];中国石油大学;2010年
3 宋立芳;金属矿地震勘探成像方法技术应用研究[D];长安大学;2012年
4 王亚强;改进的阻尼最小二乘层析算法及影响因素分析[D];燕山大学;2013年
【共引文献】
中国期刊全文数据库 前10条
1 郜延红,周云轩,刘万崧;地球物理位场可视化建模初步探讨[J];长春科技大学学报;2000年02期
2 孙建国;何洋;;基于波前构建的射线追踪:一种Java实现[J];吉林大学学报(地球科学版);2007年04期
3 韩复兴;孙建国;杨昊;;基于二维三次卷积插值算法的波前构建射线追踪[J];吉林大学学报(地球科学版);2008年02期
4 张营;冯晅;刘财;赵旭;兰慧田;张瑾;卢勇旭;杨庆节;;转换点位置对转换波速度分析精度的影响[J];吉林大学学报(地球科学版);2010年S1期
5 孔雪;李振春;叶月明;黄建平;徐秀刚;;多尺度地震资料联合成像[J];吉林大学学报(地球科学版);2011年S1期
6 曹小林,洪学海,曹俊兴;面波波形反演中的模拟退火法[J];成都理工学院学报;2000年03期
7 王志,贺振华,薛彬,黄德济,王熙明,周亚同;广角反射研究及高斯射线正演模拟分析[J];成都理工学院学报;2002年01期
8 吴建平,冯锐;井间电磁波CT中的非块体分段反演算法[J];CT理论与应用研究;1993年03期
9 张兆田,邱佩璋;一种优化的图象重建外插算法[J];CT理论与应用研究;1993年04期
10 闻则刚,何璟琳;用于势场反演的特殊解法[J];CT理论与应用研究;2001年02期
中国重要会议论文全文数据库 前10条
1 刘福平;杨长春;;孔隙度敏感系数的数值计算[A];中国科学院地质与地球物理研究所二○○三学术论文汇编·第四卷(油气资源)[C];2003年
2 耿小兵;张仲山;;核钻地弹爆炸地下爆震破坏范围确定方法探讨[A];全国第六届核仪器及其应用学术会议论文集[C];2007年
3 龙士国;马天朗;;稀土超磁致伸缩换能器的研制及其在桥梁检测中的应用[A];第六届中国功能材料及其应用学术会议论文集(3)[C];2007年
4 余易豪;舒玉峰;丁国辉;李铭铭;;震波CT探测技术及其应用[A];矿山地质灾害成灾机理与防治技术研究与应用[C];2009年
5 李满树;;同时反演复杂地区速度与界面的速度建模方法[A];油气地球物理实用新技术——中国石化石油勘探开发研究院南京石油物探研究所2004年学术交流会论文集[C];2004年
6 龙士国;马天朗;;超磁震源声波CT仪的研制及其在工程无损评价中的应用[A];中国仪器仪表学会第九届青年学术会议论文集[C];2007年
7 罗水余;申秋学;;井间地震射线层析成像的分辨率讨论[A];勘探地球物理2005学术交流会论文集[C];2005年
8 夏建军;宋强弓;尹业才;宋锡熊;;复杂地表条件下的地震波场外推法模拟分析[A];中国地球物理第二十一届年会论文集[C];2005年
9 李满树;;一种利用多种地震波走时建立表层速度模型的方法[A];中国地球物理学会第22届年会论文集[C];2006年
10 徐兆涛;;单程波动方程叠前深度偏移并行计算[A];中国地球物理学会第二十三届年会论文集[C];2007年
中国博士学位论文全文数据库 前10条
1 刘俊峰;海拉尔盆地复杂储层地震识别方法研究[D];中国地质大学(北京);2010年
2 鲁晶津;地球电磁三维数值模拟的多重网格方法及其应用研究[D];中国科学技术大学;2010年
3 曹令敏;地震层析成像方法的研究及其在金属矿复杂模型模拟中的应用[D];中国地质大学(北京);2011年
4 王连山;廊固凹陷大兴砾岩体天然气藏评价预测研究[D];中国地质大学(北京);2011年
5 汪自军;基于卫星临边辐射的大气痕量气体含量反演研究[D];吉林大学;2011年
6 张向宇;工业炉温度场可视化与辐射特性参数解耦重建研究[D];华中科技大学;2011年
7 张风雪;有限频体波走时层析成像及其在华北地区的应用[D];中国地震局地球物理研究所;2011年
8 彭远新;能源消费与经济发展的多尺度分析和反演[D];南京师范大学;2011年
9 韩翀;海相碳酸盐岩地震储层预测技术及应用[D];成都理工大学;2011年
10 朱四新;琼东南盆地深水区低幅构造速度研究与PSDM应用[D];成都理工大学;2011年
中国硕士学位论文全文数据库 前10条
1 田禹;煤岩双相介质弹性波场数值模拟[D];山东科技大学;2010年
2 李要;稀土超磁致伸缩换能器的动态特性及其应用[D];湘潭大学;2010年
3 宋小超;南海均衡重力异常研究[D];中国地质大学(北京);2011年
4 郭伟;隧道掘进超前地质勘探的研究与仿真[D];沈阳工业大学;2011年
5 赵旭;转换波多次覆盖道集抽取研究[D];吉林大学;2011年
6 杨文;利用背景噪声方法研究云南强震前后的波速变化[D];中国地震局地震预测研究所;2011年
7 刘乾;茶汗通古盆地侏罗纪地层分布特征研究[D];新疆大学;2011年
8 冯红刚;隐患空区群下采场结构参数优化研究[D];中南大学;2011年
9 王冉;复杂地电模型的2.5-D直流电阻率边界单元数值模拟[D];中南大学;2011年
10 杨婷;华北地区地壳上地幔S波速度结构研究[D];中国地震局地球物理研究所;2011年
【同被引文献】
中国期刊全文数据库 前10条
1 李红;;复杂地区静校正技术应用研究[J];安徽地质;2007年01期
2 彭直兴,张芬,周熙襄,钟本善;地震波旅行时非线性/线性联合插值法[J];成都理工大学学报(自然科学版);2005年03期
3 王金峰;罗省贤;李录明;;改进的神经网络并行算法及其在地震初至拾取中的应用[J];成都理工大学学报(自然科学版);2007年03期
4 张朋,张兆田;几种CT图像重建算法的研究和比较[J];CT理论与应用研究;2001年04期
5 吴万龙,李元景,桑斌,何文俊,李玉兰;CT技术在安检领域的应用[J];CT理论与应用研究;2005年01期
6 徐明才,高景华,荣立新,刘建勋;从金属矿地震方法的试验效果探讨其应用前景[J];中国地质;2004年01期
7 皮金云;段云卿;康平;;综合折射静校正方法在柴达木盆地的应用[J];地球科学(中国地质大学学报);2007年04期
8 朱丹;波前快速推进法三维走时计算技术[J];大庆石油地质与开发;2001年05期
9 高磊;潘树林;周熙襄;史艳华;;利用转换波初至进行转换波剩余静校正[J];大庆石油地质与开发;2008年01期
10 高星;地震层析成像研究的回顾与展望[J];地球物理学进展;2000年04期
中国重要会议论文全文数据库 前2条
1 何晓冬;夏忠谋;张国保;;模型法静校正在西部地区的应用[A];中国地球物理·2009[C];2009年
2 张正峰;韩立国;;变换坐标系下的转换波波动方程延拓静校正[A];中国地球物理·2009[C];2009年
中国博士学位论文全文数据库 前8条
1 巩向博;金属矿地震高精度成像与数据处理方法研究[D];吉林大学;2011年
2 王克斌;复杂地表条件下初至折射波静校正方法研究[D];成都理工大学;2004年
3 邹强;山地静校正若干问题研究[D];成都理工大学;2004年
4 于富文;中国西部地区静校正方法对比研究[D];中国地质大学(北京);2006年
5 马德堂;弹性波场数值模拟及井间地震初至波旅行时层析成像[D];长安大学;2006年
6 高磊;地震初至波自动拾取和折射波表层调查方法研究[D];成都理工大学;2008年
7 刘素芹;基于应用网格环境的复杂地表波动方程基准面静校正研究[D];中国石油大学;2008年
8 王兆国;东北地区地震综合信息联合成像与震源机制反演[D];吉林大学;2010年
中国硕士学位论文全文数据库 前10条
1 朱龙生;多方位角逆VSP层析成像[D];长安大学;2003年
2 赵德亨;折射波广义互换法在城市地质调查中的应用[D];吉林大学;2005年
3 彭直兴;地震波初至旅行时射线追踪正演方法研究[D];成都理工大学;2005年
4 刘洪雷;复杂地区的折射波静校正应用研究[D];中国地质大学(北京);2006年
5 秦亚玲;沙漠地区地震资料处理技术研究[D];中国地质大学(北京);2006年
6 尹纪超;基于网格的球面波波前射线追踪方法研究[D];成都理工大学;2006年
7 张丽;复杂地区静校正技术研究[D];长安大学;2006年
8 姜国博;初至波表层模型层析正反演及层析静校正研究[D];中国地质大学(北京);2007年
9 喻维秋;初至波表层模型层析反演及其应用[D];成都理工大学;2007年
10 田连玉;基于模拟退火和遗传算法的剩余静校正方法研究[D];中国石油大学;2007年
【二级引证文献】
中国期刊全文数据库 前7条
1 秦宁;李振春;杨晓东;贾正虹;岳玉波;;叠前多级优化联合偏移速度建模[J];地球物理学进展;2013年01期
2 秦宁;李振春;杨晓东;张凯;王俊;;共散射点道集与角道集串级优化叠前偏移速度分析[J];吉林大学学报(地球科学版);2013年02期
3 刘润泽;于师建;田清伟;陆二男;;基于波前最小走时单元的三角网射线追踪全局算法[J];地球物理学进展;2013年02期
4 秦宁;李振春;张凯;杨晓东;;弹性矢量波成像域联合层析反演[J];地球物理学报;2013年09期
5 潘兴祥;秦宁;曲志鹏;李振春;;叠前深度偏移层析速度建模及应用[J];地球物理学进展;2013年06期
6 秦宁;李振春;杨晓东;周卿;陈富强;桑运云;;自动拾取的成像空间域走时层析速度反演[J];石油地球物理勘探;2012年03期
7 秦宁;李振春;杨晓东;岳玉波;张凯;张晴;;高斯波束角道集的旅行时层析速度分析[J];石油地球物理勘探;2013年03期
中国重要会议论文全文数据库 前1条
1 秦宁;李振春;;基于角度域共成像点道集的走时层析速度反演[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
中国硕士学位论文全文数据库 前3条
1 孟庆利;基于大炮初至的层析静校正算法研究与实现[D];中国石油大学;2011年
2 高杨;三维高密点地震数据预处理方法研究[D];南京理工大学;2013年
3 曾同生;基于角道集的层析速度反演方法研究[D];中国石油大学(华东);2012年
【相似文献】
中国期刊全文数据库 前10条
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期
中国重要会议论文全文数据库 前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-62791813
  • 010-62985026