中文摘要：

针对复杂介质地震波形反演问题，以流体饱和多空隙介质弹性波方程为研究对象，开展基于稀疏约束和全变差正则化方法的反演理论和反演方法研究。构造稀疏约束正则化泛函和相应的半光滑迭代法，通过大量的数值实验，探讨反演参数的小波系数稀疏性表达和基函数的选择，从而极大地减少计算量和对数据的需求，并有效地识别小的异常构造。研究非连续参数反演的全变差正则化方法和相应的半光滑方法，实现对尖、角和振荡边界的有效识别。设计全变差正则化和稀疏约束正则化混合反演策略，给出权重的选择方法，构造及分析特殊的半光滑迭代法。进一步实现全局收敛的自适应小波多尺度反演算法，完成对复杂构造的地震波全波形精细反演成像。本项目的研究成果和研究方法是全新的，同时具有普遍适用性，对于其它复杂介质地震波形反演问题乃至其它地球物理反演问题具有重要的参考价值。

结论摘要：

完全的非线性地震波形反演由于可以有效地排除偶然因素的影响，提高计算的可靠性和稳定性，并且更符合复杂的地质构造而越来越受到人们的重视。由于数据噪声的存在以及不适定性，人们必须采用正则化方法来进行求解。而传统的Tikhonov正则化是通过压制大的分量而产生光滑的近似解，对于重构阶跃性介质以及对尖、角边界和小的异常构造等不连续介质的识别十分不利，并不适合复杂构造的全波形反演。为此，我们引入先进的数学理论—稀疏约束正则化方法和全变差正则化方法，开展基于稀疏约束和全变差正则化方法的复杂介质地震波形反演理论和反演方法研究。探索此类问题的解的稀疏性，实现对复杂介质弹性波场的精细反演成像。 针对带有PML吸收边界条件的频率域二维声波方程，将经典的二阶差分格式与45度旋转差分格式加权结合，形成离散格式，然后采用CARP-CG方法并行求解方程组。方法在计算效率、存储量等方面呈现出优势，表现出可应用于更大规模的弹性波正演模拟方面的潜能，为下一步开展弹性波方程正反演研究提供了有力的技术支持。针对各向同性频率域弹性波动方程采用二阶异构有限差分法进行离散，采用稀疏存储和稀疏矩阵的LU分解直接求解离散的线性方程组，实现了弹性波方程的快速正演模拟。 针对地震波全波形反演问题，引入了非线性稀疏约束正则化方法，并采用对偶方法求解稀疏约束泛函的极小点。基于二维声波方程波形反演问题进行了数值模拟,针对不同模型对稀疏约束正则化方法进行了测试。结果表明，稀疏约束正则化方法对不连续介质模型的介质边缘具有良好的识别能力。完成了Kelvin粘弹性介质弹性波方程地震波形反演的多尺度全变差正则化方法构造，并分别针对分层模型和孔隙模型的参数反演进行了数值模拟，实现了对阶跃性介质的有效识别。为了减弱对稀疏约束项的光滑性要求，引入了贝叶斯推断，产生一组收敛于后验分布的采样点，记录了采样点的条件期望、方差、置信区间等具有统计意义的结果。在没有光滑性的要求下，贝叶斯—稀疏约束正则化方法对孔洞介质模型和分层介质模型有良好的识别能力。特别地，当减少观测数据时，方法仍能获得较好的反演结果。 整合与发展已经取得的成果，并将其应用于更为复杂的基于弹性波方程的地震勘探全波形反演是我们下一步的工作。