三维直流电阻率与地震走时联合反演策略研究与算法实现-固体地球物理专业论文.docx

University University of Science and Technology of China A d issertation for doctor’S deg ree Developing Effective Strategies and Numerical Algorithms for Joint 3D DC resistivity and Seismic travel time inversion Author’S Name： Ji Gao speciality： Solid Geophysics Supervisor： Prof．Haijiang Zhang r lnlShed tlme： May,2016 万方数据 摘要众所周知，不同的地球物理勘探方法对地下介质的成像各有其优缺点，而且 摘要 众所周知，不同的地球物理勘探方法对地下介质的成像各有其优缺点，而且 都存在某种程度的不确定性和多解性。引起地球物理勘探不确定性和多解性的原 因包括：有限的观测数据、环境干扰和由于仪器精度存在的观测误差、地球物理 正演方程和真实物理过程存在的偏差、地球物理等价场的存在、反演算法存在的 偏差等因素。因此为了更可靠地刻画地下介质，多种地球物理数据之间的联合反 演解释就成为一种有效的途径。 通过对目前国内外联合地球物理反演技术方法的学习研究，本文围绕着三维 直流电阻率与地震走时成像两种方法展开联合地球物理成像研究。论文的主要研 究内容和成果如下： 11对单独的直流电阻率二维和三维反演以及地震走时的二维和三维反演做 了系统的基础性研究。在直流电阻率成像方面，采用了基于伴随矩阵(Adjoint)方 法的共轭梯度反演算法，该反演方法有效的解决了电阻率法雅克比矩阵求取效率 低下及其存储巨大的问题；在地震走时成像方面，采用基于程函方程的 MSFM(Multi Stencils Fast Marching Methods)多模板行进方法求取走时场，然后采 用梯度法求取射线路径，并通过最小二乘方法求解反演矩阵得到速度结构。通过 对直流电阻率与地震走时的单独反演研究得出，地震体波走时成像与直流电阻率 成像，均面临着成像阴影区问题。对于直流电阻率成像，电场线在高阻区域分布 较少，造成高阻区域分辨率较低。对于地震走时成像，地震射线对低速区域覆盖 较差形成阴影区，造成低速区域分辨率降低。 21研究了Gallado和Meju(2003)提出的基于交叉梯度结构约束的联合地球 物理成像方法，即在要求不同的物性模型拟合各自对应数据的同时，模型之间的 结构要求一致，即两个模型之间的交叉梯度值趋于零。并通过合成理论数据测试 模型之间交叉梯度的性质，为联合反演提供理论基础。 31为了更有效地实现基于交叉梯度的结构约束的联合反演，我们提出了一 种新的基于交替结构约束的联合反演流程，即交替反演不同的数据而且在反演一 种数据时要求对应的模型与另一个模型结构一致。新的算法能够更容易地把单独 的反演系统耦合在一起，而且也更容易建立结构约束和数据拟合之间的平衡。基 于新的联合反演流程，我们测试了基于交叉梯度结构约束的二维跨孔直流电阻率 和地震走时联合成像。合成数据测试表明，我们提出的交替结构约束流程能够很 好地实现基于交叉梯度结构约束的联合成像。与单独成像结果相比，全通道电阻 率和地震走时联合成像更可靠地确定了电阻率和速度异常。 4)在交替结构约束联合反演流程的基础上，进一步提出一种新的基于模型 万方数据 摘要更新量与交叉梯度结构一致性的联合反演算法，该方法完全把基于交叉梯度的结 摘要 更新量与交叉梯度结构一致性的联合反演算法，该方法完全把基于交叉梯度的结 构约束和基于数据的模型反演部分解耦，同时可以考虑数据拟合约束和结构一致 性约束之间的平衡。新的反演策略极大的减小了实现基于交叉梯度结构约束进行 联合反演的复杂性，可以在不改变原来单独反演方法程序的基础上实现联合反演， 并取得了较好的效果。 5)利用新的基于模型更新量与结构一致性的联合反演策略，进行了直流电 阻率与地震走时三维理论数据测试，并在测试过程中详细讨论了反演参数选取的 原则与方法，同时加入了基于先验信息及根据模型分辨率自适应选择权重因子的 策略。模型数据测试表明三维联合直流电阻率与地震走时反演具有较好的改善单 独反演的效果，具有较广泛的适用性。 关键词：电阻率成像体波走时成像结构相似性交叉梯度联合反演 II 万方数据 AbslractAbstract Abslract Abstract It is well known that different geophysical exploration methods have their advantages and disadvantages to im