粒子群-梯度算法在频率域地震波形反演中的应用

本文将优化领域应用较广的全局随机非线性粒子群算法与局部迭代梯度法相结合,构造了一种粒子群-梯度算法,并将其应用于频率域波形速度结构反演.数值实验结果表明,粒子群-梯度算法能继承梯度法快速收敛和粒子群法全局寻优的特点,适用于频率域波形反演问题,算法具有一定的抗噪能力,无论在计算精度还是在降低解的非唯一性方面,都有较明显的改善.     

频率域粘弹性声波透射波形速度反演

在用稀疏矩阵的LU分解技术对频率域粘弹性声波方程进行直接求解的基础上,根据失配函数二范数最小准则, 用预条件梯度类方法对粘弹性声波介质的速度结构进行了逐频反演. 局部非均匀介质模型和层状介质模型速度结构反演的实验结果表明,不同频率能够反映地下介质的多尺度物性结构（低频数据对应与介质物性的大尺度结构）,用低频反演结果作为高频反演的初值逼近这一顺序模式,能大大改善反演过程中解的非唯一性. 而且,在反演过程中用Hess矩阵的对角线元素来做梯度类方法的预条件算子,能够吸收了高斯牛顿法的二次收敛优势, 使得本文算法具有较快的收敛速度.      

三角网格有限元法声波与弹性波模拟频散分析

本文对声波与弹性波方程进行有限元法离散,构造有限元法频散关系的一般特征值问题,分析了时间离散格式为中心差分的三角网格有限元法声波与弹性波模拟的频散特性. 比较了三种质量矩阵即分布式质量矩阵、集中质量矩阵和混合质量矩阵对有限元法频散的影响;选取四种典型三角网格,分析了混合质量矩阵有限元（MFEM）频散的方向各向异性;数值频散、方向各向异性随插值阶数的增加逐渐减弱,当空间为三阶插值时,频散主要表现为随采样率的变化而几乎无明显方向各向异性, 其频散幅值也较小. 控制其他影响因素不变的情况下,研究了不同波速比介质中弹性波的数值频散. 最后给出了三角网格MFEM的数值耗散性.