收藏本站
《中国地质大学》 2012年
收藏 | 手机打开
二维码
手机客户端打开本文

频率域全波场有限差分数值模拟及特征分析

周聪  
【摘要】:全波形反演(FWI)是一种具有挑战的数据拟合方法,它基于全波场模拟并从地震记录中提取定量信息,而频率域全波形反演是其重要的一个分支。它的核心内容是有效的频率域正演模拟方法,正演在很大程度上决定其精度和效率。由于在实际地震资料中,面波信息通常包含在地震记录里,而且面波和体波呈现不同的物理特性,包含面波的波形反演仍然是一个挑战,强能量的面波会对波形反演中的最优化策略产生额外的困难。因此,考虑自由边界的频率域全波场数值模拟是研究这一困难的基础。 频率域地震波数值模拟需要求解一个大型稀疏矩阵方程。在通常情况下,该矩阵方程具有高稀疏性、非对称性以及非正定性的特点,可以采用直接法和迭代法进行求解,二者也有各自的优点和缺点。不同方法之间的对比分析可以帮助认识和改善各自方法的缺点,提高计算效率。 本文从频率空间域弹性波波动方程出发,首先实现了两种频率域正演模拟方法——旋转九点差分法和最优化二十五点差分法,并从差分格式的特点(即稀疏矩阵特点)、网格频散、计算精度及计算效率对比分析了这两种方法的优劣。接着研究了两种求解大型稀疏矩阵的解法——共轭梯度法和LU分解法,也从计算精度、计算内存需求及计算效率分析了这两种方法的特点。然后,讨论了频率空间域弹性波正演模拟中自由边界的处理,将时间域正演中的隐式自由边界条件引入,在频率空间域模拟出了瑞雷面波。并利用Radon变换求取了面波的频率速度谱,通过与理论面波相速度曲线的对比论证了面波的正确性。最后设计了几种简单和复杂介质,分析了面波和体波的传播特点。 根据本文的研究成果,总结出如下结论: 一、旋转九点差分法和最优化二十五点差分法都属于混合网格法,都能有效应用于频率域地震波数值模拟。由于二十五点差分法用到更多周围点的信息,因此它的相速度和群速度的网格频散明显小于九点差分法的。根据本文的理论频散结果分析,可知为保证相速度误差在±1%之内,九点差分法中每一横波波长大约需要10个网格点,而二十五点法中每一横波波长只需3.3个左右。而且二十五点法的相速度频散曲线不受泊松比的影响; 二、在保证一定的相速度误差范围内,利用旋转九点差分法和最优化二十五点差分法求得的解的精度相当。虽然在相同网格数大小的情况下,九点差分法所需的计算时间和内存需求均小于二十五点差分法,但在相同大小模型的情况下,九点差分法需要将模型剖分得更细才能达到一般的精度要求; 三、共轭梯度法和LU分解法均能有效解决大型稀疏矩阵方程,其中稀疏矩阵的元素分布对LU分解法的内存需求影响很大。通过对比发现,在求解过程中,共轭梯度法所需计算内存很小,但计算时间长;而LU分解所需计算内存大,但计算时间短。并且随着模型越来越大,两者在耗费时间和所需内存的差距也越来越大; 四、最佳匹配层(PML)吸收边界条件也能较好处理频率空间域数值模拟。但由于每一个单频切片包含所有的时间信息,所以频率域的正演模拟的对吸收条件要求很高,需要将PML吸收边界设置更厚才能得到满意的效果; 五、将隐式边界条件应用于频率域数值模拟,模拟出瑞雷面波。通过从模拟记录中提取出的瑞雷面波频散曲线与理论面波相速度曲线的对比论证了本文面波模拟的正确性。并得出利用最优化二十五点差分法时,自由边界处每一横波波长需包含的网格点数超过20个才能达到较高的精度。
【学位授予单位】:中国地质大学
【学位级别】:硕士
【学位授予年份】:2012
【分类号】:P631.4

手机知网App
【参考文献】
中国期刊全文数据库 前8条
1 吴国忱;罗彩明;梁楷;;TTI介质弹性波频率-空间域有限差分数值模拟[J];吉林大学学报(地球科学版);2007年05期
2 吴国忱,王华忠;波场模拟中的数值频散分析与校正策略[J];地球物理学进展;2005年01期
3 卞爱飞;於文辉;周华伟;;频率域全波形反演方法研究进展[J];地球物理学进展;2010年03期
4 殷文;印兴耀;吴国忱;梁锴;;高精度频率域弹性波方程有限差分方法及波场模拟[J];地球物理学报;2006年02期
5 吴国忱,梁锴;VTI介质频率-空间域准P波正演模拟[J];石油地球物理勘探;2005年05期
6 梁锴;吴国忱;印兴耀;;TTI介质qP波方程频率—空间域加权平均有限差分算子[J];石油地球物理勘探;2007年05期
7 张金海;王卫民;赵连锋;姚振兴;;黏声波介质傅里叶有限差分法正演模拟[J];石油地球物理勘探;2008年02期
8 马召贵;王尚旭;宋建勇;;频率域波动方程正演中的多网格迭代算法[J];石油地球物理勘探;2010年01期
中国硕士学位论文全文数据库 前2条
1 韩冰;利用25点有限差分方法对二维弹性波正演模拟[D];中国地质大学(北京);2007年
2 戚艳平;三维VTI介质qP波正演方法研究[D];中国石油大学;2008年
【共引文献】
中国期刊全文数据库 前10条
1 时刚;高广运;张军;;评价水泥搅拌桩复合地基加固效果的新思路[J];地下空间与工程学报;2007年05期
2 张波;张仲培;郑亚东;张进江;闫淑玉;;膝褶带构造及其地球物理正演[J];北京大学学报(自然科学版);2012年01期
3 林万顺;瞬态面波法在工程勘察中的应用[J];北京水利;1999年05期
4 宋建勇;郑晓东;秦臻;苏本玉;;基于多网格的频率域全波形反演(英文)[J];Applied Geophysics;2011年04期
5 马见青;李庆春;;面波压制的TT变换法[J];吉林大学学报(地球科学版);2011年02期
6 韩利;韩立国;李翔;王德利;崔杰;;二阶声波方程频域PML边界条件及频域变网格步长并行计算[J];吉林大学学报(地球科学版);2011年04期
7 李智强;范宜仁;邓少贵;杨震;陈华;卢志远;;自适应逐次超松弛在大斜度井泥浆侵入模拟中的应用[J];测井技术;2009年06期
8 李胜平,吴杰芳,李声平;粗粒料压实干密度无损检测试验研究[J];长江科学院院报;2004年06期
9 苏向前;刘康和;;面波法与单孔检层法波速测试的工程应用[J];长江工程职业技术学院学报;2006年03期
10 韦刚,赵明阶,黄卫东,吴国雄;土石路堤模型试验及压实度的波动测试[J];重庆交通学院学报;2004年04期
中国重要会议论文全文数据库 前10条
1 冯英杰;杨长春;吴萍;;地震波有限差分模拟综述[A];中国科学院地质与地球物理研究所2007学术论文汇编(第六卷)[C];2008年
2 徐义;张剑锋;;地震波数值模拟的非规则网格PML吸收边界[A];中国科学院地质与地球物理研究所2008学术论文汇编[C];2009年
3 黄建南;黄清和;邹燕红;;强夯后地基土强度变化的瑞雷波研究[A];2010年全国工程地质学术年会暨“工程地质与海西建设”学术大会论文集[C];2010年
4 楼加丁;;面板堆石坝中的堆石体密度测试的研究[A];纪念贵州省水力发电工程学会成立20周年论文选集[C];2005年
5 赵天库;;瑞雷面波技术在碎石桩地基加固处理中的应用[A];上海市岩土工程检测中心论文集(1995—2005)[C];1995年
6 郭新毅;马力;吴国清;;浅海沉积层中界面波的传播以及掩埋目标的回波[A];中国声学学会2007年青年学术会议论文集(上)[C];2007年
7 耿光旭;李华平;刘金川;;用瑞雷波法探测水库坝堤质量[A];勘探地球物理2005学术交流会论文集[C];2005年
8 赵勇刚;刘江;;瞬态面波剖面法在滑坡探查中的应用[A];勘探地球物理2007年学术交流会论文集[C];2007年
9 郭良辉;孟小红;;地震瑞雷波速度反演求取P-SV波横波静校正量[A];中国地球物理学会第22届年会论文集[C];2006年
10 李晶;张海燕;唐湘蓉;;利用瞬态瑞雷波法调查表层结构的可行性探讨[A];中国地球物理学会第22届年会论文集[C];2006年
中国博士学位论文全文数据库 前10条
1 徐慧峰;基于弹性波的地下管线水平导向钻随钻探测理论与技术研究[D];国防科学技术大学;2010年
2 邵广周;多阶模式瑞利波频散特征与反演研究[D];长安大学;2009年
3 朱四新;琼东南盆地深水区低幅构造速度研究与PSDM应用[D];成都理工大学;2011年
4 陈山;求解波动方程的龙格—库塔型方法及其地震波传播模拟[D];清华大学;2010年
5 肖柏勋;高模式瑞雷面波及其正反演研究[D];中南大学;2000年
6 严细水;弹性地基的表面振动特性及其应用研究[D];浙江大学;2003年
7 吕苗荣;信息提取与矿山信息系统的研究开发[D];中南大学;2003年
8 裴建中;多年冻土地区路基纵向裂缝形成机理及处治对策研究[D];长安大学;2004年
9 杨建国;物探检测公路地质的综合方法研究及其工程应用[D];重庆大学;2005年
10 蔡靖;岩土隐蔽结构缺陷量化分析及小波变换应用研究[D];天津大学;2004年
中国硕士学位论文全文数据库 前10条
1 蒋爱辞;基于离散元的土石混合料动力特性分析[D];郑州大学;2010年
2 曹雪峰;微动勘察方法中的数学方法研究与应用[D];中国地质大学(北京);2011年
3 杨奎;微动面波H/V谱探测浅部地层横波速度结构的应用研究[D];中国地质大学(北京);2011年
4 彭文;长平高速公路K73+200~K73+500段路面塌陷勘察评价及治理研究[D];吉林大学;2011年
5 秦波;基于Rayleigh波理论的地基土动参数反演[D];北京交通大学;2011年
6 陈强;基于ZigBee技术的车辆交通信息采集系统研究[D];南京理工大学;2011年
7 方大为;岩溶及土洞地球物理探测的应用研究[D];中南大学;2011年
8 李培;基于奇异值分解的瑞雷面波加权阻尼最小二乘反演[D];中南大学;2011年
9 唐磊;基于改进BISQ机制的双相介质地震波场数值模拟及频散校正研究[D];中南大学;2011年
10 黄敏;铁路路基病害车载雷达探测技术研究[D];中南大学;2011年
【二级参考文献】
中国期刊全文数据库 前10条
1 吴国忱;罗彩明;梁楷;;TTI介质弹性波频率-空间域有限差分数值模拟[J];吉林大学学报(地球科学版);2007年05期
2 熊高君,贺振华,张毅祥,姜绍仁,阎贫,李显贵,张琳;三维混合延拓一步法波动方程正演模拟[J];成都理工学院学报;1999年01期
3 杨文采;非线性波动方程地震反演的方法原理及问题[J];地球物理学进展;1992年01期
4 潘冬明,杨顶辉,张慧;含煤层地质环境下地震波场的数值模拟[J];地球物理学进展;1997年01期
5 杨文采;非线性地球物理反演方法:回顾与展望[J];地球物理学进展;2002年02期
6 杨慧珠,张友生,陶果;井眼条件下弹性波传播问题的三维有限差分数值模拟[J];地球物理学进展;2003年02期
7 裴正林,牟永光;地震波传播数值模拟[J];地球物理学进展;2004年04期
8 吴国忱,王华忠;波场模拟中的数值频散分析与校正策略[J];地球物理学进展;2005年01期
9 苑书金;;叠前地震反演技术的进展及其在岩性油气藏勘探中的应用[J];地球物理学进展;2007年03期
10 吴国忱;梁锴;;VTI介质qP波方程高精度有限差分算子[J];地球物理学进展;2007年03期
中国博士学位论文全文数据库 前2条
1 骆志刚;典型结构大型线性方程组的分布式并行算法研究[D];中国人民解放军国防科学技术大学;2000年
2 赵连锋;井间地震波速与衰减联合层析成像方法研究[D];成都理工大学;2002年
中国硕士学位论文全文数据库 前1条
1 张永杰;大型稀疏方程组预处理迭代快速求解技术研究[D];西北工业大学;2006年
【相似文献】
中国期刊全文数据库 前10条
1 魏建新,孟平,赵群;固体物理模拟实验:波场观测[J];石油物探;1997年02期
2 李卫忠,张明振,王成礼,马在田,王华忠;压制面波的波场分离方法[J];石油地球物理勘探;1998年05期
3 李明,麦振洪,崔树范;近完美晶体X射线统计动力学[J];物理学报;1994年01期
4 James W.Rector;井间地震数据的多域分析和波场分离[J];勘探地球物理进展;1994年05期
5 刘建华,郭宝琦,樊庆铎;波动方程平面波深度偏移[J];哈尔滨工业大学学报;1994年01期
6 詹毅,周熙襄;小波包分析与奇异值分解(SVD)叠前去噪方法[J];石油地球物理勘探;2004年04期
7 陈朴深;全反射时的表面波特性[J];湖州师范学院学报;1989年06期
8 宫培谦,高锡铭,李幼铭,束沛镒;深源远震波形反演介质层状结构[J];地壳形变与地震;1995年04期
9 Peter S.Rowbotham;在井间地震反射处理中通过三维滤波实现波场分离[J];勘探地球物理进展;1995年02期
10 Mark Yan Schaack;西德克萨斯碳酸盐岩储层的井间高分辨率成象 第二部分:波场模拟和分析[J];勘探地球物理进展;1996年02期
中国重要会议论文全文数据库 前10条
1 周腾飞;胡天跃;;用逆频散方法压制频散面波[A];中国地球物理学会第二十七届年会论文集[C];2011年
2 雷江莉;陈文超;王伟;高静怀;;基于MCA的面波分离方法[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
3 房立华;吴建平;明跃红;王未来;;华北地区勒夫面波噪声层析成像研究[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
4 曹思远;朱艳;;基于Chirplet变换的面波去噪研究[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
5 李多;周仕勇;冯永革;陈永顺;姜明明;;鄂尔多斯地区上地幔岩石圈结构的面波层析成像结果[A];中国地球物理学会第二十七届年会论文集[C];2011年
6 赖思静;贾学明;杨建国;;路基压实度的面波检测技术[A];2004年道路工程学术交流会论文集[C];2004年
7 郑静静;印兴耀;张广智;;基于Curvelet变换的面波压制方法[A];中国地球物理·2009[C];2009年
8 李俊邦;宋建国;;利用多子波地震道分解技术压制面波[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
9 陈伟;王尚旭;;基于保角变换的角形裂纹波场研究[A];中国地球物理学会第二十七届年会论文集[C];2011年
10 李貅;薛国强;王妙月;;从瞬变电磁场到波场转换的研究进展[A];中国地球物理学会第二十七届年会论文集[C];2011年
中国重要报纸全文数据库 前10条
1 刘帅;三维垂直地震剖面与海底地震全波场联合成像研究[N];科技日报;2007年
2 高丽辉;与大山较量的勘探劲旅[N];中国石油报;2007年
3 黄琦;探测一号步入探测之旅 中星20号正式交付用户[N];中国航天报;2004年
4 ;配备电视功能的液晶台式电脑[N];大众科技报;2000年
5 本报通讯员 赵家彬 本报记者 顾春阳;吉大“找水”团队的科研神话[N];长春日报;2011年
6 本刊记者 陈欣 谷玉强;10″,可贵的攻关[N];秦皇岛日报;2010年
7 黄嵘;缘何逆势而动 QFII自圆其说不知为何[N];中国经营报;2005年
8 王跃;多波多分量:地震勘探技术正向纵深发展[N];中国石化报;2006年
9 吴国平;在“反复”中把握好机会[N];上海证券报;2007年
10 特约撰稿人 西 华;五因素支撑春季行情[N];证券日报;2003年
中国博士学位论文全文数据库 前10条
1 邓帅奇;全空间弹性波场数值模拟与逆时偏移成像方法研究[D];中国矿业大学;2012年
2 戚志鹏;瞬变电磁三维合成孔径延拓成像方法研究[D];长安大学;2011年
3 张煜;部分饱和多相介质中Rayleigh面波传播特性研究[D];中国地质大学;2012年
4 陈凌;小波束域波场的分解、传播及在地震偏移成像中的应用[D];中国地震局地球物理研究所;2002年
5 吴中;电磁波在旋磁材料中的非线性传播特性[D];上海大学;2001年
6 魏建新;正交各向异性介质波场特性的物理模拟研究[D];中国地震局地球物理研究所;2003年
7 周建宇;井间地震研究与应用[D];中国科学院研究生院(兰州地质研究所);2002年
8 张慧;非均质储层双变网格正演模拟和弹性逆时偏移方法研究[D];中国石油大学;2011年
9 唐晖;地震波动数值模拟及稳定性研究[D];中国地震局工程力学研究所;2011年
10 刘成斋;泥岩裂缝油藏地球物理特征与检测方法[D];成都理工大学;2004年
中国硕士学位论文全文数据库 前10条
1 周聪;频率域全波场有限差分数值模拟及特征分析[D];中国地质大学;2012年
2 张少梅;地震资料面波正演模拟及去噪技术研究[D];中国石油大学;2011年
3 肖扬;改善低频噪音和波场影响的偏移成像条件研究[D];中国地质大学(北京);2011年
4 李军锋;高精度有限差分波场正演模拟方法研究及应用[D];成都理工大学;2012年
5 卢秀丽;海底地震仪观测的粘弹性地层地震波场正演模拟[D];吉林大学;2004年
6 纪聪智;海相地震数据波动方程基准面校正方法研究[D];中国石油大学;2010年
7 陈高祥;地形起伏条件下的直接叠前深度偏移方法研究[D];中南大学;2012年
8 桂丽华;基追踪在面波分离与压制中的应用[D];西南交通大学;2010年
9 王治;孔隙地层震电测井波场分波分析[D];哈尔滨工业大学;2010年
10 李俊邦;基于多子波地震道分解技术消除面波[D];中国石油大学;2011年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026