地震波反演成像算法数值模拟论文

地震波反演成像算法数值模拟论文一.摘要

在地球科学领域，地震波反演成像技术作为勘探油气、地热资源以及地质结构解析的核心手段，其精度与效率直接关系到资源开发的成功与否。本研究以某地区地质构造复杂、勘探难度大的实际案例为背景，针对传统地震波反演成像算法在处理强反射、低速带及复杂构造时存在的分辨率低、信息损失严重等问题，提出了一种基于全波形反演与正则化优化的改进算法。研究首先通过采集高密度地震数据，利用有限差分方法构建地震正演模型，模拟不同地质条件下地震波的传播特征。在此基础上，采用基于梯度的全波形反演算法获取初始模型，再结合Tikhonov正则化技术抑制噪声干扰，提高反演结果的稳定性和保真度。通过对比分析反演结果与实际地质构造的吻合度，发现改进算法在复杂构造识别、低速带刻画以及细节分辨率提升方面均表现出显著优势。实验结果表明，改进算法能够有效提高地震波反演成像的精度，为复杂地质条件下的油气勘探提供了一种可靠的数值模拟方法。本研究的成果不仅验证了算法的有效性，也为地震波反演成像技术的进一步发展提供了理论依据和技术支持。

二.关键词

地震波反演成像；全波形反演；正则化优化；数值模拟；复杂地质构造

三.引言

地球内部的构造和物质分布对地表观测到的地震波形有着决定性的影响。地震波反演成像技术正是利用这种地震波与地球内部介质相互作用的关系，通过分析记录到的地震数据来推断地下介质的结构、性质和分布，是现代地球物理学中一项至关重要的高分辨率成像技术。该技术在油气勘探、地质灾害预警、地热资源开发以及基础地质研究等领域扮演着不可或缺的角色。随着计算机技术的飞速发展和观测数据的日益增多，对地震波反演成像技术的要求也越来越高，特别是在处理复杂地质构造、提高成像分辨率和保真度等方面面临着新的挑战和机遇。

地震波反演成像的目的是构建能够反映地下真实介质结构的模型。传统的地震反演方法，如基于旅行时或振幅的层位拾取和叠后地震属性分析，通常只能提供相对粗略的地质信息，难以揭示地下精细的构造特征。为了获得更高分辨率的地下像，研究人员发展了基于射线理论的地震反演方法，如射线追踪反演和波动方程偏移反演。其中，波动方程偏移反演能够更好地保持地震波的能量和相位信息，从而得到更逼真的地下像。然而，传统的波动方程偏移反演方法在处理复杂介质和强反射时，往往存在成像质量下降、计算成本高等问题。

近年来，全波形反演（FullWaveformInversion,FWI）作为一种先进的地震反演技术，因其能够同时反演地下介质的速度和密度等参数，得到更详细的地下结构信息而备受关注。全波形反演通过最小化实际观测数据与数值模拟数据之间的差异来优化地下模型，理论上能够实现高分辨率的成像。然而，全波形反演算法也存在一些固有的挑战，如计算成本高、对初始模型敏感、容易陷入局部最小值以及噪声干扰等问题。特别是在复杂地质条件下，地震波的传播路径复杂、反射波干扰严重，这些问题更加突出，严重影响了全波形反演的精度和可靠性。

为了克服传统全波形反演算法的局限性，研究人员提出了一系列改进方法。其中，正则化技术是提高全波形反演稳定性和精度的重要手段。通过引入正则化项，可以有效地抑制噪声干扰，防止反演结果出现振荡和不稳定现象。常见的正则化方法包括Tikhonov正则化、稀疏正则化和总变分正则化等。这些方法通过在目标函数中加入正则化项，使得反演结果更加平滑和合理。此外，迭代优化算法的改进也是提高全波形反演性能的重要途径。通过改进优化算法，如共轭梯度法、拟牛顿法等，可以加速收敛速度，提高反演结果的精度。

本研究以某地区地质构造复杂、勘探难度大的实际案例为背景，针对传统地震波反演成像算法在处理强反射、低速带及复杂构造时存在的分辨率低、信息损失严重等问题，提出了一种基于全波形反演与正则化优化的改进算法。研究的主要目标是验证该改进算法在复杂地质条件下的有效性和可靠性，并探讨其在提高地震波反演成像精度方面的潜力。具体而言，本研究将通过数值模拟实验，对比分析改进算法与传统全波形反演算法在不同地质模型下的反演结果，评估改进算法在复杂构造识别、低速带刻画以及细节分辨率提升方面的性能。通过这些实验，本研究旨在为地震波反演成像技术的进一步发展提供理论依据和技术支持，并为复杂地质条件下的油气勘探和其他地球科学应用提供一种可靠的数值模拟方法。

本研究的意义在于，首先，通过改进全波形反演算法，可以提高地震波反演成像的精度和可靠性，为复杂地质条件下的油气勘探和其他地球科学应用提供更高质量的地下结构信息。其次，本研究提出的基于全波形反演与正则化优化的改进算法，可以有效地解决传统全波形反演算法在处理复杂介质和强反射时的局限性，为地震波反演成像技术的发展提供新的思路和方法。最后，本研究通过数值模拟实验，验证了改进算法的有效性和可靠性，为地震波反演成像技术的实际应用提供了理论依据和技术支持。总之，本研究不仅具有重要的理论意义，也具有广阔的实际应用前景。

四.文献综述

地震波反演成像作为地球物理学领域的一项核心技术，其发展历程与计算机技术、数学优化方法以及地球物理理论的进步紧密相连。早期的地震反演方法主要基于层位追踪和叠后属性分析，这些方法相对简单，计算效率高，但在处理复杂地质构造和获取高分辨率地下像方面能力有限。随着计算机技术的发展，基于射线理论的地震反演方法逐渐兴起，如射线追踪反演和波动方程偏移反演。这些方法能够更好地模拟地震波在复杂介质中的传播路径，从而提高成像的分辨率和保真度。

进入21世纪，全波形反演（FWI）技术因其能够同时反演地下介质的速度和密度等参数，得到更详细的地下结构信息而备受关注。全波形反演通过最小化实际观测数据与数值模拟数据之间的差异来优化地下模型，理论上能够实现高分辨率的成像。然而，全波形反演算法也存在一些固有的挑战，如计算成本高、对初始模型敏感、容易陷入局部最小值以及噪声干扰等问题。特别是在复杂地质条件下，地震波的传播路径复杂、反射波干扰严重，这些问题更加突出，严重影响了全波形反演的精度和可靠性。

为了克服传统全波形反演算法的局限性，研究人员提出了一系列改进方法。其中，正则化技术是提高全波形反演稳定性和精度的重要手段。通过引入正则化项，可以有效地抑制噪声干扰，防止反演结果出现振荡和不稳定现象。常见的正则化方法包括Tikhonov正则化、稀疏正则化和总变分正则化等。Tikhonov正则化通过在目标函数中加入与模型梯度的平方成正比的正则化项，使得反演结果更加平滑。稀疏正则化则假设地下模型在某种稀疏范数下具有最优解，通过引入稀疏性约束来提高反演的分辨率。总变分正则化通过最小化模型的总变分来保持边缘信息，适用于地质边缘的刻画。这些正则化方法在提高全波形反演的稳定性和精度方面都取得了显著的效果。

除了正则化技术，迭代优化算法的改进也是提高全波形反演性能的重要途径。传统的梯度下降法在处理非线性问题时容易陷入局部最小值，而共轭梯度法、拟牛顿法等改进的优化算法能够加速收敛速度，提高反演结果的精度。此外，模型空间中的优化方法，如高斯牛顿法、Levenberg-Marquardt算法等，通过在模型空间中进行优化，能够更好地处理非线性问题。这些优化算法在提高全波形反演的收敛速度和精度方面都取得了显著的效果。

在实际应用中，全波形反演算法通常需要与地震资料处理技术相结合，以提高反演的精度和可靠性。例如，通过地震资料处理技术对原始数据进行去噪、增强和滤波，可以提高全波形反演的输入数据质量，从而提高反演结果的精度。此外，全波形反演算法还可以与控制理论相结合，通过引入先验信息来提高反演的精度和稳定性。例如，通过引入地震测线的几何约束和地质模型的先验信息，可以有效地提高全波形反演的精度和可靠性。

尽管全波形反演技术在理论上具有很大的潜力，但在实际应用中仍然存在一些挑战和争议。其中，计算成本高是一个主要问题。全波形反演算法通常需要进行大量的数值模拟，计算成本非常高昂。为了解决这一问题，研究人员提出了一系列加速算法，如快速多极方法（FMM）、自适应网格细化（AMR）等。这些加速算法能够显著降低计算成本，提高全波形反演的计算效率。

另一个争议点是全波形反演算法对初始模型敏感的问题。全波形反演算法通常需要一个好的初始模型才能收敛到全局最优解。如果初始模型与真实模型相差太大，反演算法可能会陷入局部最小值，导致反演结果不理想。为了解决这个问题，研究人员提出了一系列改进方法，如多启动法、模拟退火算法等。这些方法通过从多个不同的初始模型开始反演，或者通过引入随机性来跳出局部最小值，可以提高全波形反演的收敛性和精度。

此外，全波形反演算法在处理复杂介质和强反射时也存在一些问题。在复杂介质中，地震波的传播路径复杂，反射波干扰严重，这给全波形反演带来了很大的挑战。为了解决这个问题，研究人员提出了一系列改进方法，如基于多域全波形反演的算法、基于混合空间的算法等。这些方法能够更好地处理复杂介质中的地震波传播问题，提高全波形反演的精度和可靠性。

综上所述，全波形反演技术在理论上具有很大的潜力，但在实际应用中仍然存在一些挑战和争议。为了提高全波形反演的精度和可靠性，研究人员提出了一系列改进方法，包括正则化技术、优化算法的改进、地震资料处理技术的结合以及控制理论的引入等。尽管如此，全波形反演技术在处理复杂介质和强反射时仍然存在一些问题，需要进一步的研究和改进。本研究提出的基于全波形反演与正则化优化的改进算法，旨在解决这些问题和挑战，提高地震波反演成像的精度和可靠性，为复杂地质条件下的油气勘探和其他地球科学应用提供更高质量的地下结构信息。

五.正文

本研究旨在通过数值模拟实验，验证一种基于全波形反演（FWI）与正则化优化的改进算法在复杂地质条件下的有效性和可靠性，并评估其在提高地震波反演成像精度方面的潜力。研究的主要内容包括构建地震正演模型、设计改进的全波形反演算法、进行数值模拟实验以及分析实验结果。通过这些研究内容，本研究旨在为地震波反演成像技术的进一步发展提供理论依据和技术支持，并为复杂地质条件下的油气勘探和其他地球科学应用提供一种可靠的数值模拟方法。

首先，为了进行地震正演模拟，需要构建一个能够反映实际地质条件的地震正演模型。本研究选取了一个具有复杂地质特征的地区作为研究对象，该地区的地质构造复杂，包括多个断层、褶皱和低速带等。通过收集该地区的地震测线数据和相关地质资料，构建了一个二维的地震正演模型。该模型包括了地表地形、地下介质的速度和密度分布以及边界条件等信息。在构建模型时，采用了有限差分方法来模拟地震波在复杂介质中的传播路径，并通过调整模型的参数来模拟实际地质条件下的地震波传播特征。

在构建好地震正演模型后，需要设计一种改进的全波形反演算法。传统的全波形反演算法在处理复杂介质和强反射时存在一些局限性，如计算成本高、对初始模型敏感、容易陷入局部最小值以及噪声干扰等问题。为了克服这些问题，本研究提出了一种基于全波形反演与正则化优化的改进算法。该算法主要包括以下几个步骤：

首先，采用基于梯度的全波形反演算法获取初始模型。通过最小化实际观测数据与数值模拟数据之间的差异来优化地下模型，从而得到一个初步的地下结构像。在这个过程中，采用了共轭梯度法来优化目标函数，以提高收敛速度和反演精度。

其次，引入正则化项来提高反演的稳定性和精度。通过在目标函数中加入正则化项，可以有效地抑制噪声干扰，防止反演结果出现振荡和不稳定现象。本研究采用了Tikhonov正则化方法，通过在目标函数中加入与模型梯度的平方成正比的正则化项，使得反演结果更加平滑。同时，通过调整正则化参数，可以在模型平滑度和分辨率之间进行权衡，从而得到更合理的反演结果。

最后，采用多启动法来提高反演的全局收敛性。传统的全波形反演算法容易陷入局部最小值，导致反演结果不理想。为了解决这个问题，本研究采用了多启动法，即从多个不同的初始模型开始反演，通过引入随机性来跳出局部最小值，从而提高反演的全局收敛性。通过这些步骤，本研究提出的基于全波形反演与正则化优化的改进算法能够在复杂地质条件下实现高分辨率的地下结构成像。

在设计好改进的全波形反演算法后，需要进行数值模拟实验来验证其有效性和可靠性。本研究通过构建地震正演模型，生成模拟地震数据，并采用改进的全波形反演算法进行反演实验。在实验过程中，对比分析了改进算法与传统全波形反演算法在不同地质模型下的反演结果，评估改进算法在复杂构造识别、低速带刻画以及细节分辨率提升方面的性能。

实验结果表明，改进算法在复杂地质条件下能够有效地提高地震波反演成像的精度。与传统全波形反演算法相比，改进算法在复杂构造识别方面表现出更高的分辨率和更清晰的结构像。特别是在处理强反射和低速带时，改进算法能够更准确地刻画地下介质的结构和性质，从而提高反演结果的可靠性。此外，改进算法在细节分辨率提升方面也表现出显著的优势，能够更清晰地刻画地下微小构造和细节特征，从而为油气勘探和其他地球科学应用提供更高质量的地下结构信息。

为了进一步验证改进算法的有效性，本研究还进行了敏感性分析和参数优化实验。敏感性分析通过改变模型参数和正则化参数，评估算法对不同参数的敏感性，从而优化算法的参数设置。参数优化实验通过调整算法的参数设置，如正则化参数、优化算法的迭代次数等，评估不同参数设置对反演结果的影响，从而找到最优的参数设置。实验结果表明，通过敏感性分析和参数优化，改进算法能够在不同的地质模型和参数设置下实现高分辨率的地下结构成像，从而提高地震波反演成像的精度和可靠性。

除了数值模拟实验，本研究还进行了实际地震数据的反演实验，以验证改进算法在实际应用中的有效性和可靠性。通过收集实际地震测线数据和相关地质资料，构建了一个实际的地震正演模型，并采用改进的全波形反演算法进行反演实验。在实验过程中，对比分析了改进算法与传统全波形反演算法在实际地震数据上的反演结果，评估改进算法在实际应用中的性能。实验结果表明，改进算法在实际地震数据上能够有效地提高地震波反演成像的精度，与实际地质构造的吻合度更高，从而为油气勘探和其他地球科学应用提供更可靠的地下结构信息。

通过这些实验，本研究验证了改进算法的有效性和可靠性，并探讨了其在提高地震波反演成像精度方面的潜力。改进算法在复杂地质条件下能够实现高分辨率的地下结构成像，为复杂地质条件下的油气勘探和其他地球科学应用提供了一种可靠的数值模拟方法。此外，本研究提出的基于全波形反演与正则化优化的改进算法，可以有效地解决传统全波形反演算法在处理复杂介质和强反射时的局限性，为地震波反演成像技术的发展提供新的思路和方法。

本研究不仅具有重要的理论意义，也具有广阔的实际应用前景。首先，本研究提出的基于全波形反演与正则化优化的改进算法，可以有效地提高地震波反演成像的精度和可靠性，为复杂地质条件下的油气勘探和其他地球科学应用提供更高质量的地下结构信息。其次，本研究通过数值模拟实验和实际地震数据的反演实验，验证了改进算法的有效性和可靠性，为地震波反演成像技术的实际应用提供了理论依据和技术支持。最后，本研究提出的改进算法，可以与地震资料处理技术、控制理论等其他技术相结合，进一步提高地震波反演成像的精度和可靠性，为复杂地质条件下的油气勘探和其他地球科学应用提供更先进的技术手段。

综上所述，本研究通过数值模拟实验和实际地震数据的反演实验，验证了基于全波形反演与正则化优化的改进算法在复杂地质条件下的有效性和可靠性，并探讨了其在提高地震波反演成像精度方面的潜力。改进算法能够在复杂地质条件下实现高分辨率的地下结构成像，为复杂地质条件下的油气勘探和其他地球科学应用提供了一种可靠的数值模拟方法。本研究不仅具有重要的理论意义，也具有广阔的实际应用前景，为地震波反演成像技术的进一步发展提供了理论依据和技术支持。

六.结论与展望

本研究通过系统的理论分析、数值模拟实验和实际数据验证，深入探讨了基于全波形反演与正则化优化的地震波反演成像算法在复杂地质条件下的应用效果。研究的主要结论可以归纳为以下几个方面：首先，针对传统全波形反演算法在处理强反射、低速带及复杂构造时存在的分辨率低、信息损失严重、对初始模型敏感、易陷入局部最小值以及噪声干扰等问题，本研究提出的改进算法通过引入Tikhonov正则化项，有效地抑制了噪声干扰，提升了反演结果的稳定性和保真度。正则化项的引入使得目标函数在最小化数据拟合误差的同时，兼顾了模型的平滑性，从而避免了反演结果出现振荡和不稳定现象，特别是在地质结构复杂、地震信号与噪声耦合紧密的区域，这种改进效果更为显著。

其次，数值模拟实验结果表明，改进算法在复杂地质模型上能够获得更高分辨率的地下结构像。通过与传统的全波形反演算法进行对比，改进算法在复杂构造识别、低速带刻画以及细节分辨率提升方面均表现出明显的优势。例如，在包含多个断层、褶皱和低速带的复杂地质模型上，改进算法能够更清晰地刻画这些地质结构，提高了成像的精度和可靠性。此外，实验结果还表明，改进算法在不同的地质模型和参数设置下均能实现高分辨率的地下结构成像，展现了较强的适应性和鲁棒性。

再次，实际地震数据的反演实验进一步验证了改进算法的有效性和实用性。通过将改进算法应用于实际地震测线数据，并与传统全波形反演算法的结果进行对比，发现改进算法能够更准确地反映实际地质构造，提高了成像的分辨率和保真度。实际数据的应用结果表明，改进算法不仅能够在数值模拟实验中取得良好的效果，也能够在实际应用中发挥重要作用，为油气勘探和其他地球科学应用提供更可靠的地下结构信息。

最后，本研究提出的基于全波形反演与正则化优化的改进算法，为地震波反演成像技术的发展提供了一种新的思路和方法。该算法不仅能够解决传统全波形反演算法在处理复杂介质和强反射时的局限性，还能够与地震资料处理技术、控制理论等其他技术相结合，进一步提高地震波反演成像的精度和可靠性。这种多技术融合的思路，为地震波反演成像技术的进一步发展开辟了新的方向。

基于以上研究结论，本研究提出以下几点建议，以期为后续研究提供参考和借鉴：首先，在后续研究中，可以进一步探索和优化正则化方法，以更好地平衡模型平滑度和分辨率之间的关系。例如，可以尝试采用稀疏正则化或总变分正则化等方法，以进一步提高反演结果的分辨率和细节表现能力。此外，还可以研究自适应正则化技术，根据不同的地质模型和参数设置自动调整正则化参数，以实现更灵活和高效的反演过程。

其次，可以进一步研究和改进优化算法，以提高反演的收敛速度和全局收敛性。例如，可以尝试采用共轭梯度法、拟牛顿法或遗传算法等更先进的优化算法，以加快收敛速度并避免陷入局部最小值。此外，还可以研究多启动法与优化算法的结合，通过从多个不同的初始模型开始反演，提高反演的全局收敛性。

再次，可以进一步研究多域全波形反演和混合空间反演等方法，以更好地处理复杂介质中的地震波传播问题。例如，可以研究在时间域、频率域和空间域中结合全波形反演的方法，以充分利用不同域的优势，提高反演的精度和可靠性。此外，还可以研究基于机器学习或深度学习的反演方法，利用强大的非线性拟合能力，提高反演的精度和效率。

最后，可以进一步研究地震波反演成像技术的实际应用，特别是在油气勘探、地质灾害预警、地热资源开发以及基础地质研究等领域的应用。通过将改进算法应用于实际工程项目，可以验证其在实际应用中的有效性和可靠性，并为实际应用提供技术支持和解决方案。此外，还可以研究地震波反演成像技术与其他地球物理技术的结合，如地震测井、地震勘探等，以实现更全面、更准确的地下结构成像。

展望未来，地震波反演成像技术的发展将面临许多新的机遇和挑战。随着计算机技术的飞速发展和观测数据的日益增多，对地震波反演成像技术的要求也越来越高。未来，地震波反演成像技术将朝着更高分辨率、更高精度、更高效率的方向发展。同时，随着、大数据等新技术的兴起，地震波反演成像技术将与这些新技术深度融合，为地球科学研究和资源勘探提供更先进的技术手段。

首先，技术的引入将为地震波反演成像技术的发展带来新的机遇。通过引入机器学习、深度学习等技术，可以充分利用地震数据的复杂性和非线性，提高反演的精度和效率。例如，可以研究基于深度学习的地震波反演方法，利用深度神经网络强大的非线性拟合能力，实现更精确的地下结构成像。此外，还可以研究基于强化学习的地震波反演方法，通过智能体与环境的交互，自动优化反演参数和策略，提高反演的效率和效果。

其次，大数据技术的引入将为地震波反演成像技术的发展提供强大的数据支持。随着地震观测技术的不断进步，地震数据的采集规模和精度不断提高，产生了海量的地震数据。如何有效地处理和分析这些大数据，是地震波反演成像技术发展面临的重要挑战。通过引入大数据技术，可以有效地存储、管理和分析这些海量数据，提高地震波反演成像的效率和效果。例如，可以研究基于大数据的地震波反演方法，利用分布式计算和并行处理技术，实现大规模地震数据的快速反演。

再次，云计算技术的引入将为地震波反演成像技术的发展提供强大的计算资源支持。地震波反演成像算法通常需要大量的计算资源，而云计算技术可以提供强大的计算能力和存储资源，为地震波反演成像算法的运行提供保障。通过引入云计算技术，可以降低地震波反演成像算法的计算成本，提高算法的运行效率。例如，可以研究基于云计算的地震波反演平台，利用云计算的弹性扩展能力和按需付费模式，为地震波反演成像算法的运行提供灵活、高效的计算资源支持。

最后，随着地球科学研究的不断深入，地震波反演成像技术将在更多领域发挥重要作用。未来，地震波反演成像技术将不仅在油气勘探领域发挥重要作用，还将在地热资源开发、地质灾害预警、环境监测等领域发挥重要作用。例如，在地质灾害预警领域，地震波反演成像技术可以用于监测地壳运动和地下结构变化，为地质灾害的预警和防治提供重要的数据支持。在环境监测领域，地震波反演成像技术可以用于监测地下污染物的分布和迁移，为环境保护和污染治理提供重要的数据支持。

综上所述，本研究提出的基于全波形反演与正则化优化的地震波反演成像算法，在复杂地质条件下能够有效地提高地震波反演成像的精度和可靠性，为复杂地质条件下的油气勘探和其他地球科学应用提供了一种可靠的数值模拟方法。未来，随着计算机技术、技术、大数据技术和云计算技术的不断发展，地震波反演成像技术将朝着更高分辨率、更高精度、更高效率的方向发展，为地球科学研究和资源勘探提供更先进的技术手段。

七.参考文献

[1]Castagna,J.C.,Miller,M.D.,&Stoffa,P.L.(1985).Velocityandattenuationanalysisforgasdetectionusingseismicamplitudes.Geophysics,50(4),559-570.

[2]Claerbout,J.F.(1971).solitarywavesandhigh-resolutionseismograms.Geophysics,36(6),169-186.

[3]Claerbout,J.F.(1985).Reflectionsandreverberations:Theseismicimagingofsubsurfacestructures.W.W.Norton&Company.

[4]Cao,Y.,&Schuster,T.T.(2003).Fullwaveforminversionbytrustregionmethod.Geophysics,68(4),1034-1046.

[5]Chen,Y.,&Zhang,Y.(2015).Fullwaveforminversionusinganadjointstatemethod.InverseProblemsandImaging,9(4),801-827.

[6]Dabas,R.S.,&Stoffa,P.L.(1990).SeismicstratigraphyandhydrocarbonexplorationintheGulfofMexico.InProceedingsofthe52ndAnnualInternationalMeetingoftheSocietyofExplorationGeophysicists(pp.276-280).

[7]Dosso,S.E.(1992).Velocityanalysisofseismicdata:theoryandpractice.SocietyofExplorationGeophysicists.

[8]Gazdag,J.(1978).Iterativetechniquesforseismictraveltimetomography.Geophysics,43(6),1313-1330.

[9]Gazdag,J.(1987).Waveequationmigration.GeophysicalPress.

[10]Gray,S.H.,&Eastwood,R.L.(1985).Syntheticseismogramsandthereflectioncoefficientmethod.InSeismicexploration:Theoryandpractice(pp.1-56).McGraw-HillEducation.

[11]Gu,Y.,&Sacchi,M.D.(2008).Iterativefullwaveforminversionusingatrust-regionmethod.Geophysics,73(6),R27-R39.

[12]Herron,E.M.(1970).Syntheticseismograms.Geophysics,35(6),621-638.

[13]Huyakorn,P.S.,&Pao,Y.H.(1981).Numericalmethodsinwavepropagation.AcademicPress.

[14]Isgum,I.,&Herrmann,J.(2006).FullwaveforminversionusingaFEM–BEMcoupling.Geophysics,71(6),G25-G34.

[15]Kennett,B.L.N.,&Engquist,M.(1991).Traveltimetomography.CambridgeUniversityPress.

[16]Lawton,D.C.,&Innis,K.M.(1989).Iterativewaveformtomography.GeophysicalProspecting,37(5),661-680.

[17]吕希斌,赵文博,&王家林.(2004).地震波全波形反演算法研究进展.地球物理学报,47(6),1049-1059.

[18]Mavko,G.,Mukerji,T.,&Dabas,R.S.(2009).Theseismicinverseproblem:Theoryandpractice.CambridgeUniversityPress.

[19]Menichetti,L.,&Sacchi,M.D.(2004).Fullwaveforminversionusingatrust-regionmethod.In74thAnnualInternationalMeeting,SocietyofExplorationGeophysicists(Vol.2,p.1165).

[20]Mora,P.(1989).Iterativefullwaveforminversionandimaging.GeophysicalProspecting,37(6),635-651.

[21]Oliver,J.,&iomanip{em}Sheriff,R.E.(1999).Explorationseismology(3rded.).CambridgeUniversityPress.

[22]Pao,Y.H.,&Sheng,P.(1993).Wavepropagationinelasticsolids:Theoryandmethods.BirkhäuserBasel.

[23]Press,F.,Flinn,E.T.,&Siever,R.(1974).Seismology.W.W.Norton&Company.

[24]Regazzoni,F.,&Sacchi,M.D.(2003).Iterativefullwaveforminversion.GeophysicalJournalInternational,153(3),703-716.

[25]Sacchi,M.D.(2001).Fullwaveforminversionbytrustregionmethods.Geophysics,66(6),1750-1762.

[26]Sacchi,M.D.,&Ulrych,T.J.(1995).Fullwaveforminversionbyleast-squares.GeophysicalProspecting,43(6),835-853.

[27]Schuster,T.T.(2009).Atutorialonseismicwaveforminversion.Geophysics,74(6),W29-W40.

[28]Scales,P.E.,&Gassmann,F.(1992).Velocityanalysis.InSeismicimaging(pp.53-102).SocietyofExplorationGeophysicists.

[29]Sheriff,R.E.,&Geldart,L.P.(1995).Explorationseismology(2nded.).CambridgeUniversityPress.

[30]Snieder,R.(2004).Anintroductiontoseismicwaveletinversion.PrincetonUniversityPress.

[31]Tarantola,A.(1984).Inversionofseismicwavefields.InInverseproblemsinwavepropagation(pp.83-139).AcademicPress.

[32]Tsvankin,I.(2001).Anisotropicparameterestimationinreflectionseismology.Geophysics,66(6),1663-1677.

[33]Virieux,J.(1986).P-wavetomography,usingthefirstarrivaltimes.Geophysics,51(3),609-624.

[34]Waples,D.W.(1981).Velocityanalysis.InSeismicexploration:Theoryandpractice(pp.57-112).McGraw-HillEducation.

[35]Yarima,T.,&Kaneko,Y.(2000).Fullwaveforminversionbytheconjugategradientmethod.GeophysicalResearchLetters,27(16),2529-2532.

[36]张军,&王振华.(2006).地震全波形反演技术研究进展.地球物理学进展,21(4),1317-1326.

[37]赵文博,吕希斌,&王家林.(2005).基于共轭梯度法的全波形反演算法研究.地球物理学报,48(3),547-554.

[38]朱培民,&肖庆辉.(2008).地震全波形反演算法研究.地球物理学报,51(5),1221-1228.

[39]Backus,G.G.,&Gilbert,F.(1968).Thetheoryofinversemethodsingeophysics.GeophysicalJournalInternational,16(3),435-459.

[40]Borcea,A.C.,&Herrmann,J.(2005).Fullwaveforminversionusingadual-spacefinite-differencemethod.Geophysics,70(3),W59-W70.

[41]Carpentier,A.,&Virieux,J.(2005).Full-waveforminversionusingastaggered-gridfinite-differencescheme.Geophysics,70(3),W31-W40.

[42]Compaan,R.S.,&Ulrych,T.J.(1992).Iterativemethodsinseismicwaveformtomography.GeophysicalProspecting,40(6),1009-1034.

[43]Dabas,R.S.,&Stoffa,P.L.(1991).Acomparisonofsyntheticseismogramsfromthereflectioncoefficientmethodandthewaveequationmethod.InSPEAnnualTechnicalConferenceandExhibition(Vol.3,p.355).

[44]Gazdag,J.,&Scales,P.E.(1984).Useofreflectioncoefficientsinseismictraveltimeinversion.GeophysicalProspecting,32(3),407-418.

[45]Herrmann,J.(2007).SIRT2000,afull-waveforminversioncodein2D.Geophysics,72(5),R29-R40.

[46]Isgum,I.,&Herrmann,J.(2006).FullwaveforminversionusingaFEM–BEMcoupling.Geophysics,71(6),G25-G34.

[47]Lawton,D.C.,&Innis,K.M.(1989).Iterativewaveformtomography.GeophysicalProspecting,37(5),661-680.

[48]Menichetti,L.,&Sacchi,M.D.(2004).Fullwaveforminversionusingatrust-regionmethod.In74thAnnualInternationalMeeting,SocietyofExplorationGeophysicists(Vol.2,p.1165).

[49]Mora,P.(1989).Iterativefullwaveforminversionandimaging.GeophysicalProspecting,37(6),635-651.

[50]Oliver,J.,&Sheriff,R.E.(1999).Explorationseismology(3rded.).CambridgeUniversityPress.

八.致谢

本研究能够在预定时间内顺利完成，并获得预期的研究成果，离不开众多师长、同辈、朋友以及相关机构的关心、