全波形反演是一种利用地震波传播的动力学特征来获取地下介质物性参数的反演方法,可为揭示地下精细结构提供重要依据。本文以弹性波方程作为数学模型来模拟地震波传播规律并进行相应的反演方法研究。为提高计算效率与反演结果的准确性,可将近似解析离散化(NAD)算子用于频率域弹性波方程的正演模拟。本文在频率域NAD离散的基础上推导阻抗矩阵的稀疏分块结构与反演目标函数对模型参数的梯度计算公式,由此建立基于NAD算子的频率域弹性波全波形反演方法。为验证该方法的有效性,文中通过数值实验对多种典型介质模型进行反演计算,均得到了理想的反演结果。
相似文献逆时偏移成像具有高分辨率的特点,在频率域实施能提升偏移成像效率.数值求解波动方程是逆时偏移成像的关键环节,数值算法的优劣直接决定着偏移过程的计算效率与成像结果的质量.本文在频率域近似解析离散化(NAD)方法的基础上通过优化网格差分模板系数构造改进NAD方法,保证数值离散精度的同时可缓解阻抗矩阵的病态程度,以此提高频率域正演模拟的计算效率并将该方法运用于频率域声波逆时偏移.通过分别对凹陷模型、Sigsbee2B模型以及Marmousi模型进行偏移成像且与四阶标准NAD方法和四阶普通有限差分方法进行比较,改进NAD方法在压制数值频散与提高成像效率方面的优势得到了充分体现.
相似文献为提高频率域弹性波动方程数值求解的计算效率,本文引入近似解析离散化(NAD)方法将其进行数值离散并得到大型线性代数方程组.在详细分析了相应系数矩阵的稀疏分块结构与数学性质之后,本文提出采用不精确旋转分块三角预处理子加速Krylov子空间迭代方法来快速求解该线性方程组,并利用数值试验证实这种方法在弹性波场模拟方面的数值效率.通过与另外两种经典数值方法(常规有限差分方法和交错网格有限差分方法)对多种介质模型进行波场模拟、数值频散分析以及与解析解的波形对比,NAD方法显示了其在压制数值频散和提高计算效率方面的优势以及对复杂介质模型弹性波场数值模拟的有效性.
相似文献为研究南岭-武夷交汇区深部动力学过程提供深部背景资料和科学依据,本文利用远震P波接收函数H-κ叠加和共转换点(CCP)叠加两种方法获取了研究区66个宽频带流动台站及24个固定台站下方的地壳厚度、泊松比和Moho面起伏形态,揭示了扬子地块与华夏地块地壳结构及泊松比变化特征,给出了南岭和武夷之间一条莫霍凸起带的高分辨图像.
结果显示:(1)研究区内地壳厚度平均值为31.2 km,泊松比平均值为0.23,总体呈现薄地壳、低泊松比的特点.地壳厚度从西北往东南由厚变薄,与区域地壳伸展特征相一致.(2)在韶关-赣州-吉安-南昌一线存在条带状薄地壳结构,平均值为28 km,呈南西-北东向展布,对应的泊松比值略微升高.推测地壳减薄带的形成可能与来自南海方向的地幔热流上涌有关.(3)江南造山带的泊松比整体偏低,存在两处平均值小于0.21的区域.华夏地块内地壳厚度与泊松比之间存在弱的负相关,表示随着地壳厚度的增大,铁镁质的下地壳厚度在整个地壳厚度中所占的比例减小.
相似文献