求解双相和黏弹性介质波传播方程的间断有限元方法及其波场模拟

间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相.     

双相各向异性研究、问题与应用前景

地球内部流体的存在和岩石各向异性是地下介质的两大表征,考虑地下流体和介质各向异性问题的双相各向异性理论是当今地震学和地球物理学理论与应用研究的前沿和难题之一。孔隙流体的存在、固体和流体之间的相互作用会弱化或硬化岩石的力学属性,上起声波或弹性声波速度的频散的振幅的衰减,并产生第二类压缩波。裂缝或裂隙的定向分布、岩层的旋性沉积、应力场的定向排列,都会引起传播速度的各向异性、横波发生分裂等重要现象,这些     

二维SH波方程的半解析解及其数值模拟

本文以波动理论为基础, 半解析化求解地震勘探中常用的SH波方程. 获得的主要结果包括: 给出了二维均匀介质中SH波方程的解析解; 利用Cagniard-de Hoop方法详细推导了二维双层介质中SH波方程的解析解, 获得了透射波的解析解表达式. 同时, 基于SH波方程的解析表达式, 给出了包含各种波(如直达波、反射波、首波以及透射波)的解析解和波形图. 对于比较复杂的积分型解析解, 利用数值积分方法给出了数值结果, 并与优化的近似解析离散化方法(ONADM)和4阶Lax-Wendroff修正方法(LWC)的数值结果进行了比较, 以验证解析解的正确性. 本文的研究成果有望在检验波动方程数值新方法的有效性、波传播理论分析等方面得到应用.     

地震波传播理论与地震“孕育”、发生和发展的深部介质和构造环境

地震灾害是人类生存、社会与经济发展所面临的重大自然灾害之一,由于地震震源和其周边介质在力源作用下可产生微破裂并逐渐形成"破裂链",而震源区介质的质点震动则以地震波动和辐射为"载体",且向四周传播.为此,可通过研究地震波传播波场的运动学和动力学响应探索地震"孕育"、发生和发展中的深部介质和构造环境,地震波动的轨迹和其发生的力源机制以及强烈地震活动的深层过程和动力成因,此乃是探索防震减灾和地震预测的核心所在,因此,研究地震波在复杂介质中的传播和响应,地球内部物质与能量的交换和其深层动力过程将必制约着地震的活动和"孕育"空间,这是防震减灾和地震预测中必须深化研究的重要领域.     

求解声波方程的辛可分Runge-Kutta方法

本文基于声波方程的哈密尔顿系统,构造了一种新的保辛数值格式,简称NSPRK方法.该方法在时间上采用二阶辛可分Runge-Kutta方法,空间上采用近似解析离散算子进行离散逼近.针对本文发展的新方法,我们给出了NSPRK方法在一维和二维情况下的稳定性条件、一维数值频散关系以及二维数值误差,并在计算效率方面与传统辛格式和四阶LWC方法进行了比较.最后,我们将本文方法应用于声波在三层各向同性介质和异常体模型中的波传播数值模拟.数值结果表明,本文发展的NSPRK方法能有效压制粗网格或具有强间断情况下数值方法所存在的数值频散,从而极大地提高了计算效率,节省了计算机内存.     

基于非结构网格求解三维D'Alembert介质中声波方程的并行加权Runge-Kutta间断有限元方法

间断有限元方法(Discontinuous Galerkin method,简称DGM)在求解地震波动方程时具有低数值频散、网格剖分灵活等优点,因此,为适应数值模拟对模拟精度和复杂地质结构的要求,本文提出一种新的加权Runge-Kutta间断有限元(weighted Runge-Kutta discontinuous ...     

各向异性动力学方程反演新算法

基于各向异性介质中的动力学波动方程组，并利用多波多分量观测数据，发展了一种适用于直接反演介质弹性参数和密度或地震波速度的动力学反演新算法．由于这种方法能充分地利用全波场信息，因此反演结果十分可靠．为了证实这种方法的有效性和适用范围，基于拟SH波方程，并选取了不同的初始模型，对地震波速度进行了反演．结果表明:即使初始模型参数的扰动量达到12.3％时，仍能获得相当好的结果．此外，由于方法构造的特殊性，使得这种方法的向量化和并行化计算非常容易实现．      

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

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

钻柱振动倒谱分析及其钻头源信号提取方法研究

随钻地震参考信号处理的关键是从众多噪音成分中提取较弱的钻头源信号.钻头随机激励源信号由于其频率范围宽,持续时间短,易受钻柱和其他机器噪音的影响.倒谱分析方法是一种非线性滤波技术,该技术不但可以将时域卷积信号转变为倒时域信号的相加,而且通过窗函数的选择,还可消除结构混响,提取源信号.本文根据钻柱模拟实验数据,利用了倒谱滤波技术进行瞬态源信号的提取,并用该方法进行SWD参考信号的处理,从而达到钻头源信号强化的目的.     

最小二乘支持向量回归滤波系统性能分析

支持向量机(Support Vector Machine: SVM)一直作为机器学习方法在统计学习理论基础上被研究和发展,本文从信号与系统的角度出发,证明了平移不变核最小二乘支持向量机(Least Squares SVM: LS-SVM)是一个线性时不变系统.以Ricker子波核为例,探讨了不同参数对最小二乘支持向量回归(Least Squares Support Vector Regression: LS-SVR)滤波器频率响应特性的影响,这些参数的不同选择相应地控制着滤波器通带上升沿的陡峭性、通带的中心频率、通带带宽以及信号能量的衰减,即滤波器长度越长通带的上升沿越陡,核参数值越大通带的中心频率越高,且通带带宽越宽,正则化参数值越小,通带带宽越窄(但通带中心频率基本保持恒定),有效信号幅度衰减越严重.合成地震记录的仿真实验结果表明,Ricker子波核LS-SVR滤波器在处理地震勘探信号的应用中,滤波性能优于径向基函数(Radial Basic Function: RBF)核LS-SVR滤波器以及小波变换滤波和Wiener滤波方法.