TTI介质弹性波FCT有限差分数值模拟

TTI(Tilted Transversely Isotropic)各向异性是对地下岩石中广泛存在的规则发育的裂缝和层理的一种有效的弹性近似,基于TTI介质的地震波数值模拟技术是分析地震波在复杂各向异性介质中的传播机理的有效工具.同时,高精度的数值模拟算法也能为后续的逆时偏移技术提供重要的技术支撑.由于TTI介质中地震波方程的弹性参数众多且变化复杂,常规有限差分技术在解决TTI介质正演模拟问题时往往会产生严重的数值频散现象,降低了数值模拟精度.通量校正传输(FluxCorrected Transport,FCT)技术能够有效地压制由空间离散产生的数值频散.本文将FCT技术用于TTI介质中弹性波方程的交错网格高阶精度差分正演,在数值模拟过程中通过对波场进行漫射和反漫射校正实现了空间网格频散的压制.模型模拟结果表明,与常规有限差分算法相比,本文算法能够有效的压制大网格条件下的数值频散,提高模拟精度.     

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

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

地震波数值模拟中优化的通量校正传输方法

为了有效提高地震波正演数值模拟精度,提出了一种优化的通量校正传输方法,即在每一步时间递推过程中,沿坐标轴方向和对角线方向均引入通量校正参数进行波场校正处理,以高频散的均匀介质地震模拟波场为例,研究对比了不同中心网格有限差分近似阶数、传统的FCT方法和优化的FCT方法校正得到的地震波场效果.数值计算结果表明,传统的FCT方法容易产生虚假波场,而优化的FCT方法可以削弱或消除这些假波动,并指出结合高差分阶数和应用优化的FCT方法可以更好地提高数值模拟记录的信噪比.     

横向各向同性介质地震波场数值模拟研究

地震波场数值模拟是理解地震波在地下介质中的传播特点,帮助解释观测数据的有效手段,而提高计算精度和运算效率是所有波场数值模拟方法研究所追求的目标.有限差分技术是求解波动方程计算效率最高、应用最为广泛的方法之一.但传统的有限差分技术计算过程中的数值频散问题影响了该技术的计算精度与计算效率.本文通过交错网格高阶有限差分技术与通量校正传输方法（Flux|corrected transport method,FCT）相结合, 对横向各向同性介质（Transverse isotropic medium,TI）一阶速度|应力弹性波动方程组进行了数值求解研究.波场快照数值模拟结果表明,本文研究的数值模拟方法与波动方程二阶有限差分方法、交错网格四阶有限差分方法相比,在压制网格数值频散方面有明显的优势,计算精度提高,而且可以利用较大的空间步长,提高计算效率.     

用于弹性波方程数值模拟的有限差分系数确定方法

有限差分方法是波场数值模拟的一个重要方法,交错网格差分格式比规则网格差分格式稳定性更好,但方法本身都存在因网格化而形成的数值频散效应,这会降低波场模拟的精度与分辨率.为了缓解有限差分算子的数值频散效应,精确求解空间偏导数,本文把求解波动方程的线性化方法推广到用于求解弹性波方程交错网格有限差分系数;同时应用最大最小准则作为模拟退火(SA)优化算法求解差分系数的数值频散误差判定标准来求解有限差分系数.通过上述两种方法,分别利用均匀各向同性介质和复杂构造模型进行了数值正演模拟和数值频散分析,并与传统泰勒展开算法、最小二乘算法进行比较,验证了线性化方法和模拟退火方法都能有效压制数值频散,并比较了各个算法的特点.     

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

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

TTI介质的交错网格伪P波正演方法

研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果.     

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

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

基于保幅波场分离的弹性波逆时偏移方法研究（英文）

欲实现基于弹性波方程的矢量波场逆时偏移纵、横波独立成像,必须在波场延拓过程中实现纵、横波场的分离,散度和旋度算子分离的纵、横波出现振幅与相位的畸变,导致输出成像结果的振幅失真。本文提出一种在弹性波场延拓过程中实现纵、横波保幅分离的方法,在传统的弹性波方程中加入纵波压力、纵波振动速度和横波振动速度方程,实现纵横波的矢量分解,再对分解后的矢量纵波和矢量横波做标量化合成得到保幅分离的纵、横波场,对保幅分离的纵、横波场应用成像条件,然后实现矢量波场逆时偏移的保幅纵横波成像。该方法可以保证分离后纵、横波的振幅与相位不变;同时,分解后的纵波压力和纵波振动速度可用于层间反射噪音压制和横波极性校正,提高多分量地震资料联合逆时偏移的纵、横波成像质量,从而实现保幅弹性波逆时偏移的目的,为叠前深度剖面应用于叠前反演工作奠定基础。     

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

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

矢量分离纵横波场的弹性波逆时偏移

弹性波逆时偏移是当前多分量地震资料相对准确的偏移算法,它能够形成多波模式的成像剖面,从而减少纵波勘探的多解性.本文首先依据各向同性介质中矢量分离纵横波场的速度-应力方程组,利用高阶交错网格有限差分数值方法求解弹性波方程,进而构建矢量的纵波和横波波场,不同于散度和旋度算子分离纵横波场的传统方法,文中提出的矢量分离纵横波场方法保持了原始波场的振幅和相位特征.文中也提出将震源归一化的内积成像条件应用于分离后的纯纵波和横波矢量场,由此得到的转换波成像避免了传统弹性波成像方法中出现的极性反转.水平层状和复杂构造模型测试表明,文中提出的基于矢量分离纵横波场的弹性波逆时偏移方法成像精度高,转换波成像PS和SP极性无反转,所形成的多种模式纯波剖面能够准确地对复杂地下构造成像.     

双相各向同性介质弹性波高精度波场分离数值模拟方法

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场...     

基于NAD算子和FCT技术改进的Runge-Kutta方法及其波场模拟

地震波正演数值模拟中所产生的数值频散严重地影响波场模拟的精度和分辨率.基于声波方程,以8阶精度的近似解析离散算子来校正差分算子对空间偏导数进行离散,以三阶Runge-Kutta方法对时间导数进行离散,并结合通量校正传输技术以抑制数值频散,从而建立基于NAD算子和TCF技术改进的Runge-Kutta方法.波场模拟结果表明,该方法在抑制数值频散上明显优越于传统的高阶有限差分方法.同时,该方法对模拟地震波在复杂介质中的传播有着很强的适应能力.     

二维弹性随机介质中的波场特征

通过波动方程的交错网格有限差分正演,模拟了地震波在二维弹性随机介质中的传播及其自激自收时间记录.本文将理论记录剖面分割成三个不同的时间区段,在三个时间区段上分别计算并提取剖面的统计特征.运用正交试验设计和分析方法,本文全面研究了混合型弹性随机介质的五个统计特征（自相关长度、粗糙度因子、扰动标准差、纵波背景速度、纵横波速度比）及其交互作用与各种波场特征之间的关联,并得出了若干结论.     

三维弹性波方程的修正时空优化保辛数值求解方法

正演计算是反演研究的基础,为了实现基于三维弹性波方程的全波形反演成像,发展准确、高效、低数值频散的三维正演模拟方法至关重要.为此,本文将修正保辛分部龙格-库塔格式与优化有限差分算子结合,发展了用于数值求解三维弹性波方程的修正时空优化保辛方法(MTSOS).新方法使用二级龙格-库塔格式达到了三阶时间精度,且更适用于求解非均匀介质情况下的弹性波方程,数值频散误差小于同精度保辛分部龙格-库塔(SPRK)方法的误差,提高了计算精度.波场模拟结果表明,三维MTSOS方法可以精确给出数值模拟结果,能够清晰模拟地震波传播过程中产生的各种震相、有效压制数值频散.     

一种全局优化的隐式交错网格有限差分算法及其在弹性波数值模拟中的应用

在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高.     

二维弹性波频率域15点差分格式及正演模拟

为优化二维各向同性介质中弹性波频率域正演时阻抗矩阵的结构,减小正演所需内存,提高正演效率,在25点差分格式的基础上进行适当的简化,得到了二维弹性波频率域15点差分格式.利用该格式重新计算了弹性波方程中偏微分项和加速项的差分算子,减少了计算过程中的网格节点需求,构造了优化阻抗矩阵后的频率域正演矩阵方程;推导了纵波和横波相速度的频散公式,给出了不同泊松比条件下的频散曲线,得到了相速度误差控制范围±1%时每一横波波长内网格数需求.通过对比频散曲线和简单模型数值模拟时得到的波场快照、检波点处速度分量及单炮记录,验证了15点差分格式与25点差分格式相比,具有稍严格的网格间距需求、相当的计算精度、更少的计算时间和更小的阻抗矩阵带宽等特点.最后,利用复杂模型数值模拟对本方法的适应性进行了验证.     

三维VTI介质qP波方程频率空间域有限差分高阶加权算子

波动方程有限差分法是波场模拟的一个重要方法,为解决常规有限差分法存在着数值频散的问题,本文从具有垂直对称轴的三维横向各向同性(VTI)介质频率-空间域qP波动方程出发,在常规差分算子的基础上构造了适合三维VTI介质的频率空间域有限差分优化算子,然后利用最优化理论中的Gauss-Newton法求解了优化算子的系数,使差分方程的相速度与波动方程的相速度尽量吻合,从而在理论上使网格数值频散达到极小,精度对比分析及数值测试表明,有限差分优化算子具有较高的波场数值模拟精度,有效地压制了数值频散现象,为三维VTI介质频率一空间域qP波正演模拟研究提供了理论基础.     

一种基于全局优化的交错网格有限差分法

在地震波场数值模拟中, 交错网格有限差分技术得到了广泛的应用, 但是在弹性模量变化较大时, 通常会因插值而导致模拟误差增大. 旋转交错网格可以很好地克服这个缺点, 因而适合于各向异性介质正演模拟. 但是对于同样大小的网格单元, 旋转交错网格需要的步长比常规交错网格要大, 这会使梯度和散度算子的误差增大因而更易产生空间数值频散. 针对这些问题, 本文提出了旋转交错网格与紧致有限差分相结合的方法, 并基于模拟退火算法进行全局优化, 压制数值频散, 拓宽波数范围. 数值模拟结果表明, 此方法可以有效地压制数值频散, 且具有较高的模拟精度.      

基于优化差分系数的双相介质交错网格正演模拟

交错网格有限差分方法已经被广泛应用到数值模拟和地震波传播的研究中.传统交错网格有限差分方法中,一阶空间导数的高阶差分系数是通过Taylor级数展开求取的,这种表示空间导数的方法会导致数值频散的产生.本文针对时间二阶空间十阶交错网格有限差分算法,采用最小二乘法通过改变积分区间求取一系列一阶空间导数的差分系数,分析该差分系数和传统方法求取的差分系数的频散关系.选取效果最佳的最小二乘法进行数值模拟,并与传统方法相比较.数值频散分析和弹性波场模拟分析表明:介质弹性参数和离散参数相同的情况下,采用最佳积分区间的最小二乘法更能有效地压制数值频散,比Taylor级数展开法具有更高的数值模拟精度.