【摘 要】
:
有限差分方法(Finite Difference Method,FDM)作为一种快速而精确的数值方法,在波动方程正演数值模拟领域得到了最为广泛的应用,然而,当在不规则或者地表起伏不平的区域中模拟时,基于规则网格的有限差分法对波动方程离散时会产生阶梯状近似,影响模拟的精度,给有限差分方法的求解带来困难。通过求解椭圆型偏微分方程生成的贴体网格为这种复杂区域中的有限差分波场模拟问题提供了一种有效的手段。
论文部分内容阅读
有限差分方法(Finite Difference Method,FDM)作为一种快速而精确的数值方法,在波动方程正演数值模拟领域得到了最为广泛的应用,然而,当在不规则或者地表起伏不平的区域中模拟时,基于规则网格的有限差分法对波动方程离散时会产生阶梯状近似,影响模拟的精度,给有限差分方法的求解带来困难。通过求解椭圆型偏微分方程生成的贴体网格为这种复杂区域中的有限差分波场模拟问题提供了一种有效的手段。借助贴体网格以及链式法则,将空间导数的计算从不规则的物理区域转换到规则的计算域,在计算域中继续使用有限差分方法。相比其他的不规则网格剖分方法,如插值法、映射法以及非结构网格,贴体网格在普遍性、精度和稳定性方面都有优势。本文在复杂区域综合应用了贴体网格、声波方程正演模拟、逆时偏移(Reverse Time Migration,RTM)以及完全匹配层(Perfectly Matched Layer,PML)等技术。首先通过求解椭圆偏微分方程法生成贴体网格,根据链式法则代入顺序不同,将二阶位移形式的二维声波方程改写到曲线网格中,得到了两种形式的贴体坐标二维声波方程,第二种形式的波动方程相比第一种形式而言,形式对称、更紧凑。为了压制截断边界造成的人工边界反射,推导了两种形式的声波方程对应的PML方程。在复杂起伏地表区域中,分部求和(Summation-by-Parts,SBP)有限差分方法的使用能够确保曲线网格中非均匀介质数值模拟的稳定性。SBP有限差分法特别适合求解含变系数的导数项。因此,对两种形式的波动方程及PML针对性地采用了不同的差分方法:第一种形式中,导数项的系数可预先求得,采用中心差分方法进行数值离散;第二种形式中含变系数,因此采用显式的二阶精度SBP有限差分方法进行离散,并采用Fourier谱分析方法讨论了这种离散格式的稳定性。与中心有限差分方法对比,分部求和有限差分方法的稳定性更高。四阶精度的有限差分方法在降低存储需求和提高效率方面有很大的优势。在二阶精度中心差分方法和SBP有限差分方法的基础上,将两种形式的贴体网格声波方程和PML中空间导数项的离散均扩展到四阶精度。采用Fourier谱分析方法研究了该四阶SBP差分离散格式的稳定性,得到了离散方程的稳定性条件,并验证了四阶精度SBP有限差分方法比中心有限差分方法的稳定性更高。同时,在同样尺寸的模型中比较了二阶和四阶精度SBP有限差分方法的计算时间及存储需求,验证了四阶方法的高精度、低频散特性,显示了四阶方法在提高运算效率、降低存储方面的能力。本论文将贴体网格的方法分别应用于地面地震和VSP(Vertical Seismic Profiling,垂直地震剖面)资料的逆时偏移研究,得到了曲线坐标系中各自对应的成像数据体,为解决实际地震采集中观测系统不规则、地表起伏等情况下的数据处理难题提供一种有效的技术手段。
其他文献
近地表低速层的高吸收衰减特征严重影响地震资料的分辨率。虽然现有的地震处理方法和软件能够对近地表低速层的衰减效应进行处理,但大部分处理结果只是部分消除地层衰减效应,这主要是由于缺乏对近地表低速层衰减特征更全面、准确、定量的研究与分析。利用地震物理模拟方法参数可控性高、模拟结果真实性强、数据资料信噪比高的特点,可以为全面、准确、定量的研究近地表低速层衰减特征提供高质量的地震物理模拟数据。无论是野外勘探
段塞流是油气混输过程中油气输送通道路径改变而引起的一种特殊现象,导致流体压力波动并进而和结构耦合。对于深水立管系统,这种动力耦合严重恶化了系统的动力响应,是立管的疲劳破坏的主要因素之一。由于段塞流形成过程流体瞬时压力波动及其传播的机理和规律没能突破,段塞流和立管系统的耦合作用的研究更加困难。本博士论文就是在这种背景下开展的相关理论和数值模拟工作,核心进展体现在:1、通过深入研究段塞流的形成及与立管
平潭海峡公铁大桥起于长乐,终于平潭,全长16.323千米,将人屿岛、长屿岛、小练岛、大练岛串联到一起。桥址所在处风大、浪高、水深、流急,被人们视为"建桥禁区"。2013年11月动工建设,2020年10月、12月,公路段、铁路段分别通车运营,结束了平潭岛不通铁路的历史,也拿下了世界最长跨海峡公铁大桥的纪录。
岩浆侵入活动影响或改造围岩,从而形成岩浆岩及接触变质岩油气储层。受岩浆侵入活动影响的变质围岩储层逐渐成为我国东部油气勘探新的储量增长点。深入探讨岩浆侵入对沉积围岩的影响机理具有重要意义。本论文以苏北盆地高邮凹陷的油田实例和松辽盆地南部的团山子采石场露头实例来探讨辉绿岩侵入对围岩储层的影响作用。基于高邮凹陷北斜坡地区阜宁组岩心、地震及测井等资料,研究了该地区辉绿岩及其接触变质围岩的储集空间和储层物性
本文面向大数据背景下随机退化设备剩余寿命预测的现实需求,结合随机退化设备监测大数据特点及剩余寿命预测不确定性量化这一核心问题,深入分析了机理模型与数据混合驱动的剩余寿命预测技术、基于机器学习的剩余寿命预测技术、统计数据驱动的剩余寿命预测技术以及机器学习和统计数据驱动相结合的剩余寿命预测技术的基本研究思想和发展动态,剖析了当前研究存在的局限性和共性难题.针对存在的局限性和共性难题,以多源传感监测大数
全波形反演利用观测地震数据与模拟地震数据的残差来构建目标函数,通过求取目标函数的极小值问题来获得地下介质参数,其结果分辨率高,近年来成为一个研究热点。由于单一的纵波速度信息已无法满足地震勘探成像日益复杂化和精细化的要求,本文针对时间域弹性波全波形反演展开了相关研究。本文基于时间域一阶速度-应力弹性波动方程,根据拉格朗日伴随方法理论推导出了伴随方程和相应的模型参数梯度表达式。在正演模拟时为了防止人工
传统的有限差分方法离散波动方程时分别对时间和空间的偏导数进行近似,通过使用较长的差分算子来提高空间近似的精度,抑制空间频散,然而只使用二阶差分算子近似时间偏导数,使得传统的有限差分方法虽然具有空间任意偶数阶精度,但时间上的精度只有二阶。当激发源的频率较高或使用的时间步长较大时,时间二阶差分方法会引入较大的时间频散误差,进而导致地震波形的畸变。为了提高传统有限差分方法的模拟精度,本文研究了时-空域有
根据大庆油田制定的高含水后期油层分类标准,二类油层与主力油层对比不论在纵向上还是在平面上,沉积环境变化都较大,砂体发育规模明显变小,层数增多,单层厚度变薄、渗透率变低、粘土矿物含量变大,平面及纵向非均质更严重。大庆油田主力油层已大部分被聚驱动用,三元复合驱(ASP)的可应用空间主要集中于剩余地质储量相对较大的二类油层。室内研究及矿场试验表明,ASP复合驱可在水驱的基础上提高采收率约20%。复合体系
目前,刚性摆臂式多通道高精度通径检测器广泛应用于检测管道变径、变形、内表面缺陷以及腐蚀等特征。但相对国外,国内对于刚性摆臂式多通道高精度通径检测器的研究起步较晚,直到2012年才第一次实现了刚性摆臂式多通道通径检测器的自主研发,这就造成了多通道通径检测器方面的理论研究相对较少。本文基于以上背景开展了对刚性摆臂式多通道通径检测器相关关键技术的研究。具体内容包括如下:(1)研制了适用于?355mm天然
本文基于国家自然科学基金项目“深水钻井隔水管-导(套)管力学特性研究与水下井口稳定性分析(51574261)”,研究了考虑地震作用的深水水下井口-浅层管柱的稳定性。在浅层管柱组合研究时,提出一种考虑隔水管安全余量和浅层气安全余量的深水浅层钻井液安全密度窗口计算方法,得出使用双梯度钻井(Dual-Gradient Drilling,DGD)的方法可解决窄钻井液安全密度窗口的问题,DGD技术可以减少1