耦合PML边界条件的三维大地电磁二次场有限差分法

本文将大地电磁场分解为一次场和二次场,应用交错网格有限差分法模拟计算大地电磁二次场,并引入各向异性最佳匹配层(PML)吸收边界条件作为二次场边界条件,实现了耦合PML吸收边界条件的三维大地电磁二次场有限差分正演模拟.为了确保正演的稳定性和效率,QMR求解器和磁感应矢量散度校正技术被用于PML吸收边界条件下系数矩阵的快速求解.三维模型正演响应表明,基于二次场的三维大地电磁有限差分算法具有较高的计算精度和可靠性.通过计算分析不同PML吸收因子条件的大地电磁正演结果,显示在适当的吸收因子下,PML吸收边界条件可较大幅度的减小外边界距离,从而有效的压缩模型求解空间,最终提高三维大地电磁正演模拟的效率.     

BFO-PSO算法下的弹性波数值模拟

地震波正演模拟是地震反演与成像的基础和关键,有限差分算法广泛应用于地震波数值模拟,差分算子的精度直接影响数值模拟的质量和效率.本文提出一种BFO-PSO算法下的有限差分算子优化方法,并应用其进行弹性波数值模拟.首先,将BFO算法中的趋化、复制、驱散三个步骤引入PSO算法,形成具有更好全局搜索能力和更快收敛速度的BFO-PSO混合优化算法;之后构造包含有限差分系数的目标函数,并应用BFO-PSO混合优化算法求取最优解,获得优化的有限差分算子;最后应用此优化的有限差分算子在不同模型上进行弹性波数值模拟.根据频散曲线及数值模拟结果,可以分析得出,BFO-PSO算法优化后的有限差分算子在保证计算效率的同时,具有更高的精度,可以有效压制数值频散,提高数值模拟的精度和效率.     

基于Remez迭代算法求解差分系数的双相介质自适应有限差分正演模拟

利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性.     

塔里木地区勘探地震正演模拟研究

以塔里木前陆盆地地质背景为模型,采用三维任意差分精细积分方法并行算法实现了库车地区三维正演模拟.三维任意差分精细积分方法通过时间域采用局部积分半解析方法求得波动方程的递推算子,与常规的差分法相比计算精度有较大提高;文中计算稳定性根据实际算例进行分析,采用稳定因子约束,得到较好的稳定性;边界条件采用改进的自适应吸收边界,并通过串行程序并行化,大大减少了三维正演模拟的耗时,完成了大数据量三维正演模拟.     

裂缝型储层时空双变正演模拟研究

本文将深层裂缝型储层视为各向同性背景下的小尺度非均质体,采用时空双变正演模拟算法模拟其地震响应,并将Lanczos滤波算子推广到该算法中,在保证最优的计算效率及内存的前提下,解决了常规变网格算法长时间采样下的不稳定问题.测试结果验证了本文算法对深层微构造的模拟精度及效率.最后将本算法应用到裂缝储层正演模拟中,定量分析了不同裂缝开度、密度和倾角下的地震振幅响应规律,为裂缝型储层的检测和解释提供了可靠的理论指导.     

双相介质地震波场数值模拟的迭积微分算子及其PML边界条件

将基于计算数学中Forsyte 广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度--应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演模拟方法,来解决复杂非均匀介质模型中的波场传播问题.为了大幅衰减人工边界引起的反射,本文将完全匹配层(Perfectly Matched Layer,PML)吸收边界条件引入到所构建的方法中,以解决迭积微分算子法的边界问题.以二维波动方程为例,用迭积微分算子法实现了双相介质的地震波场正演模拟,模拟结果表明,双相介质模型较好地解释了含流体孔隙特性.同时也表明迭积微分算子法是一种非常实用、有效的数值模拟方法.     

基于CFS-CPML边界处理的LOVE面波有限差分模拟

建立近地表横波速度模型时通常需要用到面波分析的方法.Love面波是在低速层分界面附近传播的一种SH型不均匀平面波,本文利用高阶有限差分算子、合理的自由地表边界条件以及CFS-CPML吸收边界条件,获得了高精度的Love面波波场记录,并与传统的分裂式完全匹配层得到的波场记录作对比,体现了CFS-CPML吸收边界条件的优越性.在此基础上将数值模拟提取的频散特征与理论的频散特征进行对比,证明二者非常吻合,验证了Love面波有限差分模拟的精度很高,可以用来研究复杂情形下的Love面波频散特征,并分析了Love面波位移应力等一系列正演特征.     

瞬变电磁Crank-Nicolson FDTD三维正演

时域有限差分(FDTD)方法使用Yee网格剖分电磁场的空间采样,通过时间步迭代实现电磁场数值模拟,具有内存消耗低、计算简单等特点,常用于瞬变电磁三维正演.然而,常规FDTD方法的时间迭代步长Δt受Courant-Friedrich-Lewy(CFL)条件严格限制,过多的迭代次数以及过密的采样往往导致计算速度慢、累积误差不断增大.本文提出一种不受CFL条件约束的无条件稳定隐式差分算法Crank-Nicolson FDTD(CN-FDTD)用于瞬变电磁三维正演.基于Crank-Nicolson差分方法对Maxwell方程组重新离散,空间网格仍然采用Yee元胞,时间步进采用在整时间步电场、磁场同时采样的策略,建立无条件稳定FDTD格式,突破CFL条件限制.与常规FDTD交替采样相比,CN-FDTD电场、磁场同时采样的策略构成的隐式差分格式,需要求解大型稀疏矩阵方程组.通常,瞬变电磁三维正演模型中产生的矩阵阶数往往较大,需要占用大量内存和求解时间.为解决上述问题,采用Crank-Nicolson-cycle-sweep-uniform(CNCSU-FDTD)方法近似求解CN-FDTD方程,在保证求解精度的同时,计算效率大幅提高.在边界条件处理上,采用双线性变换推导了复频率参数完全匹配层(CFS-PML)吸收边界.采用均匀半空间模型、四类三层模型进行精度验证,发现CN-FDTD三维正演结果与解析解、线性数字滤波解吻合较好.之后,与接触带上的低阻复杂模型进行对比,结果显示CN-FDTD正演结果与矢量有限元、有限体积法以及FDTD计算结果吻合较好.在此基础上,研究了时间步放大对CN-FDTD计算精度的影响,发现最大时间步放大到常规FDTD的3200倍时才会在晚期出现较明显的误差.在一台CPU为Intel Core i5-7300HQ的笔记本电脑单线程计算条件下,模拟到关断后30 ms仅需要50 min.在进行并行化后,将有望实现复杂模型分钟级的三维正演,从而为三维反演提供可靠、快速的正演方法.     

波场模拟有限差分参数优化选取

有限差分方法(Finite-difference Method, FD)广泛用于地震波场数值模拟，但其存在固有的数值频散问题，影响模拟的计算效率和数值精度.本文主要研究了有限差分方法的空间数值频散误差和网格划分精度以及差分算子的关系，基于计算量最小准则，提出了最优化有限差分参数选取流程，为有限差分数值模拟参数选取提供理论指导.本文主要工作包括：(1)提出了空间数值频散正变换过程(Forward Space Dispersion Transform, FSDT)方法，该方法可以高效模拟出不同网格划分精度(波长采样点数)的带有空间数值频散的波场；(2)提出了波场空间数值频散误差衡量准则，可以定量地判断出数值模拟导致的波形频散程度，选取合适的频散误差阈值；(3)研究了给定空间数值频散误差阈值下，差分算子系数、差分算子阶数、网格划分精度与计算量之间的关系.文中基于雷米兹交换方法(Remez Exchange Method, RE)和泰勒级数展开方法(Taylor-series Expansion Method, TE)的差分系数，在空间数值频散误差阈值0.01时，数值模拟了不同差分算子阶数、网...     

波场模拟中的数值频散分析与校正策略

波动方程有限差分法正演模拟,对认识地震波传播规律、进行地震属性研究、地震资料地质解释、储层评价等,均具有重要的理论和实际意义.但有限差分法本身固有存在着数值频散问题,数值频散在正演模拟中是一种严重的干扰,会降低波场模拟的精度与分辨率.针对TI介质波场模拟的交错网格有限差分方法,本文从空间网格离散、时间网格离散和算子近似等三个方面对其产生的数值频散进行了分析,并结合其他学者的研究成果给出了TI介质波场模拟中压制数值频散的方法与策略：在已知介质频散关系时,对差分算子可实施算子校正;通过提高差分方程的阶数来提高波场模拟精度;采用流体力学中守恒式方程的通量校正传输方法来压制波场模拟中的数值频散;在实际正演模拟时,采用交错网格高阶有限差分方程,不仅在空间上采用高阶差分,而且在时间上也要采用高阶差分,否则只在单一方向上（空间或时间）提高方程的阶数对压制数值频散也不会取得理想的效果.     

二维频率空间域25点优化系数差分格式弹性波数值模拟(英文)

频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播.但是该方法的最大瓶颈是对于计算机内存的巨大需求.我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量.同时为了减少短筹分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数.为了克服边界反射,采用了最佳匹配层吸收边界条件.数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度.使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率.该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础.     

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题，会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散，本文针对VTI介质地震波数值模拟问题，构造了频率-空间域qP波波动方程高精度有限差分优化算子，根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析，并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟，通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度，有效压制了传统有限差分算子数值模拟中的数值频散现象，提高了有限差分算子精度，为VTI介质频率-空间域qP波正演模拟奠定了基础.     

探地雷达三维高阶时域有限差分法模拟研究

探地雷达数值模拟中,时域有限差分法在时间和空间上一般采用二阶精度的中心差分近似(FDTD(2,2)),其形式简单,但数值色散误差较大,在复杂模型模拟时不能很好地反映模型的精细变化.高阶时域有限差分法能很好地改善数值色散带来的误差,提高模拟精度.本文基于三维高阶时域有限差分法的基本原理实现了探地雷达正演模拟,采用单轴各向异性完全匹配层(UPML)作为吸收边界条件,可以有效地吸收外向传播的电磁波,在大大地提高计算效率的同时,也能很好地改善边界的吸收效果.分析对比正演模拟结果,通过三维高阶时域有限差分正演能获得目标体准确电磁响应信息,并能很好的提高模拟精度.     

基于多重网格有限元法的大地电磁二维正演计算研究

为提高大地电磁正演计算速度,开展了基于多重网格有限元法的大地电磁二维正演模拟计算研究.将稳定双共轭梯度算法作为多重网格法的细网格松弛迭代算法,插值算子采用完全加权算子,限制算子设计基于网格单元面积率,使多重网格法更适于求解大型复系数方程组.二维均匀半空间模型、低阻体模型和高阻体模型的大地电磁正演模拟结果表明:当计算量较小时(网格剖分数量少),多重网格法在计算效率方面并未有优势,网格剖分数量较大时,多重网格有限元算法在收敛速度方面的优势明显,多重网格有限元法的大地电磁正演精度优于一般数值算法.这为三维多重网格有限元的大地电磁正演研究奠定了基础.     

窗函数交错网格有限差分算子及其优化方法

交错网格有限差分的方法现已被广泛运用到地震波正演模拟中,而正演的精度将会直接影响到后续反演、偏移成像的精度。有限差分正演模拟研究的关键问题之一是如何有效地压制数值频散。窗函数法选取适当的窗函数去截断伪谱法的空间褶积序列,从而得到优化的有限差分算子以压制数值频散。传统窗函数法得到的有限差分算子,在低波数域内具有较高的精度,而在高波数域内,精度迅速下降。在此基础上,本文将交错网格有限差分系数的求取转化为最小二乘问题,将窗函数截断得到的交错网格有限差分算子作为迭代初值,设定误差范围确定优化区间,并采用共轭梯度法迭代求解。不失一般性,本文选取了常用的三种窗函数去截断得到有限差分算子及其优化差分算子,并将优化前后差分算子做对比验证。理论分析和数值模拟的结果表明:在窗函数的基础上,使用本文最小二乘优化方法得到的交错网格有限差分算子比传统窗函数法交错网格有限差分算子具有更高的精度,能够更好地压制数值频散。     

Lebedev网格改进差分系数TTI介质正演模拟方法研究

本文采用一种新的交错网格-Lebedev网格(LG)进行TTI介质的正演模拟研究,避免了Virieux标准交错网格(SSG)算法在处理TTI、单斜等各向异性介质时波场插值引入的数值误差,提高了模拟精度.在方法实现过程中,本文针对有限差分正演模拟面临的网格频散与边界反射两个关键性问题分别做了优化,并通过模型试算验证了它们的有效性与可行性:(1)结合最小二乘思想推导出新的频散改进差分系数(DIC),该系数比Taylor系数更能有效地压制粗网格引起的数值频散,可以节约内存,提高计算效率;(2)将分裂的多轴完全匹配层(M-PML)吸收边界条件引入到LG算法中,解决了传统PML边界条件在某些各向异性介质中的不稳定现象并且具有较好的边界吸收效果.     

瞬变电磁场有限差分与数字滤波双模型三维正演方法

瞬变电磁场数字滤波算法计算时间短,精度高,但难以实现三维正演计算;有限差分算法可以实现三维正演计算,但计算时间长,当模型网格剖分数量不足时计算精度较低.结合三维有限差分正演算法和一维数字滤波正演算法各自优势,提出双模型三维正演方法,可以减少计算时间,提高计算精度.首先,给出了瞬变电磁场三维有限差分正演算法和一维数字滤波正演算法.然后,通过理论分析和公式推导,提出了通过三维有限差分正演算法计算异常场,通过一维数字滤波算法计算背景场,然后叠加得到总场的双模型方法,并给出了具体计算公式和流程.最后,建立了经典的Newman均匀介质含低阻长方体模型,以及Commer层状介质含垂直接触带模型,分别采用双模型方法和常规三维有限差分算法进行了正演计算,对计算结果进行了对比分析.结果表明:在相同模型网格数量条件下双模型方法计算精度显著高于常规三维有限差分算法;双模型方法可以通过较少的模型网格数量取得高精度的计算结果,从而可以在保证计算精度的前提下显著提高计算效率.     

时域有限差分法在地质雷达二维正演模拟中的应用

本文从地质雷达正演原理着手,分析了差分格式中半空间步长与半时间步长的实现方法,同时通过分析数值频散的产生,进而推导出了理想频散关系和超吸收边界条件．最后文中还对比了有无边界条件的雷达正演模拟效果和精度．     

高精度频率域弹性波方程有限差分方法及波场模拟

有限差分方法是波场数值模拟的一个重要方法,但常规的有限差分法本身存在着数值频散问题,会降低波场模拟的精度与分辨率,为了克服常规差分算子的数值频散,本文采用25点优化差分算子,再根据最优化理论求取的优化系数,建立了频率空间域中弹性波波动方程的差分格式;为了消除边界反射,引入最佳匹配层,构造了各向同性介质中弹性波方程在不同边界和角点处的边界条件. 最后由弹性波波动方程和边界条件,通过频率域有限差分法,分别利用不同震源对弹性波在均匀各向同性介质、层状介质及凹陷模型中的传播过程进行了数值正演模拟,得到了单频波波场、时间切片和共炮点道集,为下一步的研究工作（如成像、反演）提供了研究基础.     

双轴各向异性介质中回线源瞬变电磁三维拟态有限体积正演算法

采用模拟离散的有限体积法实现了双轴各向异性地层回线源瞬变电磁三维正演.首先引入内积定义,采用自然边界条件,将瞬变电磁法的控制方程转化为弱形式表示.将计算区域划分为一系列的控制体积单元,采用交错网格对控制方程进行模拟有限体积空间离散,包括旋度算子离散和空间内积离散.基于斯托克斯定理的旋度积分定义公式实现旋度算子离散.中点平均实现电导率双轴各向异性的空间内积离散,从而得到离散化的控制方程.时间步迭代采用无条件稳定的欧拉后向差分格式.并通过均匀全空间中稳定电流回线源的磁场解析表达式得到回线源初始时刻的电磁场分布.为了同时保证计算精度和效率,本文采用分段等间隔的时间步迭代,利用直接法求解器PARDISO实现其快速求解.最后通过对比层状模型和各向异性半空间模型的正演计算结果,验证了本文算法的计算精度和计算效率;计算三维双轴各向异性模型的正演响应可知,水平方向电导率变化对电磁响应产生显著影响,而垂直方向的电导率变化对电磁响应几乎没有影响.产生这一现象的主要原因是回线源产生的感应电流主要是水平方向的,因此响应主要受到水平方向电导率的影响,垂直方向的电导率影响很小.