波动方程数值模拟与粘滞波形反演方法研究

《波动方程数值模拟与粘滞波形反演方法研究》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、中图分类号：P631.4单位代码：11414学号：2013315031题目波动方程数值模拟与粘滞波形反演方法研究演方法研究学科专业地质资源与地质工程研究方向地球物理勘探博士生陈汉明指导教师周辉教授二○一七年五月StudyonNumericalSimulationofWaveEquationsandViscoacousticFullWaveformInversionDissertationsubmittedtoChinaUniversityofPetroleum,BeijinginpartialfulfillmentoftherequirementsforthedegreeofDoctorof

2、Engineering(OrDoctorofNaturalScience)ByChenHanming（GeologicalResourcesandGeologicalEngineering）DissertationSupervisorProf.ZhouHuiMay,2017博士学位论文独创性声明郑重声明？？本博士学位论文是作者个人在导师的指导下，独立进行研宄工作所取得的成果。除了文中特别加以标注和致谢的地方外，论文中不包含其他个人和集体己经发表或撰写的研究成果，也不包含为获得中国石油大学或者其它单位的学位或证书所使用过的材料。对本研究做出贡献的个人和集体，均已在论文中做了

3、明确的说明并表示了谢意。作者和导师完全意识到本声明产生的法律后果并承担相应责任。°？ｒ＾－作者签名：曰期：ｎＴｉ＇丨＿（＆丨导师签名：个曰期：丨７博士学位论文版权使用授权书本学位论文作者及指导教师完全了解中国石油大学（北京）学位论文版权使用的有关规定，使用方式包括但不限于：学校有权保留并向有关部门和机构送交学位论文的复印件和电子版；允许学位论文被查阅和借阅；学校可以公布学位论文的全部或部分内容，可以采用影印、缩印或扫描等复制手段保存和汇编学位论文；可以将本学位论文的全部或部分内容编入有关数据库进行检索。本学位论文作者如需公开出版学位论文的

4、部分或全部内容，必须征得导师书一面同意，且须以中国石油大学（北京）为第署名单位。（ｆ＇作者签名：日期：１ｊ。＇ｘ丨丨：日期导师签名：了－－１摘要摘要传统的有限差分方法离散波动方程时分别对时间和空间的偏导数进行近似，通过使用较长的差分算子来提高空间近似的精度，抑制空间频散，然而只使用二阶差分算子近似时间偏导数，使得传统的有限差分方法虽然具有空间任意偶数阶精度，但时间上的精度只有二阶。当激发源的频率较高或使用的时间步长较大时，时间二阶差分方法会引入较大的时间频散误差，进而导致地震波形的畸变。为了提高传统有限差分方法的模拟精度，本文研究了时-空域有限差分方法，该方法将波动

5、方程的时间和空间离散格式视为一个整体，利用差分格式的频散关系确定差分系数，能更好地平衡时间和空间的近似精度。具体地，本文分别研究了适用于二阶声波方程和一阶弹性波动方程的空间任意偶数阶，时间四阶和六阶精度的有限差分方法，并分别在二维和三维情况下基于一般的矩形和长方体网格推导了差分系数的表达式。其中，弹性波场的数值模拟是基于一个新的速度-拟应力弹性波动方程，数值模拟该方程能获得解耦的纵波和横波波场。频散分析和数值实验证实，时-空域有限差分方法的频散误差更小，在相同的精度下比较，该方法允许使用更大的时间步长，因此提高了波场数值模拟的计算效率。速度和地震品质因子（Q）是地震勘探领域最常用的两个参数，

6、本文研究了利用粘滞全波形反演方法同时反演速度和Q的可行性。本文所使用的粘滞声波方程是最新的变分数阶常Q波动方程，该方程符合地震勘探领域最常用的常Q模型。针对目前变分数阶波动方程数值模拟面临的困难，本文提出了两种有效的数值模拟方法，能很好地处理分数阶随空间变化的问题，提高了粘滞波场数值模拟的精度与计算效率。在此基础上，本文发展了分数阶常Q波动方程的全波形反演方法，推导了表征参数的梯度公式，提出了一种新的非衰减型伴随算子，避免了波场反传过程中的数值不稳定，同时还提高了伴随波场的计算效率。为了避免子波估计不准确对反演结果造成的影响，本文发展了不依赖于震源的粘滞声波全波形反演方法。为了提高全波形反演

7、的计算效率，本文采用平面波反演策略减少反演的数据量，并将这一策略与不依赖震源的反演方法相结合，取得了理想的反演效果。数值实验证实本文发展的粘滞全波形反演方法可用于速度估计、Q估计以及双参数的同时估计。关键词：时空域有限差分；分数阶常Q波动方程；全波形反演；Q反演-II-ABSTRACTStudyonNumericalSimulationofWaveEquationsandViscoacousticFullWa