多学科地震预报研究中的地球物理联合反演问题

对震源区岩石破裂用扩容物理数学模型求定岩石强度和破裂的整体的临界参数（前兆）的方法已发展起来了，这个参数就是单位体积裂纹的密度根据监测到的不同性质的地球物理场的异常变化数据和联合多学科反演问题的数值解，可对这一组合前兆进行定期判定。     

综合地球物理联合反演综述

阐述了综合地球物理研究的必要性，方法学及研究思路，指出联合反演是综合地球物理重要的定量解释工具，通过回顾联合反演的研究现状与进展，指出联合反演的发展方向。     

地球物理联合反演研究综述

地球物理联合反演由于使反演问题的非唯一性得到有效限制而越来越受到人们的重视。本文概述了联合反演的发展现状及实现的方法，并讨论了其发展趋势及其局限性，指出地球物理联合反演是地球物理数据分析的理想工具，而非线性联合反演方法则是地球物理联合反演发展的方向。     

神经网络在地球物理反演中的应用

本文简要介绍神经网络与神经计算机，神经元模型，层状网络结构，误差逆传播算法及其在地球物理反演中的应用。     

地球物理联合反演研究的现状与趋势分析

本文从回顾和分析地球物理联合反演的必要性、分类方式,以及基本原理入手,通过分析和总结前人的研究思路,讨论了联合反演研究和应用中的关键性问题,归纳出联合反演研究的三个基本环节,即模型参数化、目标函数构造和反演算法选择,并对其加以讨论, 同时指出非线性联合反演方法和模型参数化研究将成为未来联合反演研究的主流和热点.     

基于交叉梯度约束的重力、磁法和大地电磁三维联合反演

为了降低单一地球物理方法反演的多解性及受噪声的影响程度,本文围绕重力、磁法和大地电磁法开展了三维联合反演的研究.重、磁采用基于对数障碍法的正则化反演算法,大地电磁使用limited-memory Broyden-Fletcher-Goldfarb-Shanno（L-BFGS）反演算法,引入交叉梯度函数实现了三种物性结构的相互耦合,最终开发出一套重磁电三维联合反演算法,并实现MPI并行加速计算.通过理论模型算例验证了算法的准确性,结果表明：不论是单棱柱体模型还是组合棱柱体模型,联合反演结果相较单独反演对于异常体的空间形态刻画以及物性数值恢复具有较好的提升;单棱柱体模型算例使得异常体的物性参数（密度、磁化率和电阻率）更加接近于真实的物性参数;组合棱柱体模型的联合反演结果不仅仅消除了围岩物性参数的假异常,而且还增强了异常体边界结构的恢复程度.

     

基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演

为了有效解决目前大地电磁和地震走时资料单方法反演结果一致性不好的问题,同时克服基于岩石不同物性参数间关系耦合约束联合反演的局限性,本文研究了基于交叉梯度耦合约束的大地电磁与地震走时资料的三维联合反演算法.以较为成熟的天然地震走时资料三维正反演和大地电磁三维正反演算法为基础,实现了具有共同的反演网格,以交叉梯度结构耦合约束,并能同时获得电阻率和速度模型的三维联合反演算法.分别利用单棱柱体模型和双棱柱体模型合成数据进行了联合反演试算.结果表明:无论是单棱柱体模型还是双棱柱体模型,联合反演结果比单独反演对异常体的空间形态都有更好的恢复,其中单棱柱体模型反演的异常体电阻率更接近于真实电阻率,双棱柱体模型的联合反演结果不仅消除了围岩的部分电阻率假异常,而且增强了对异常体深部速度结构特征的恢复程度.联合反演还能同时改善电阻率和速度反向变化异常体的单独反演结果,进一步证明交叉梯度耦合不依赖于岩石物性关系,而强调地下结构的相似性,具有更普遍的适用性.     

基于交叉梯度约束的地震初至纵波与瑞雷面波联合反演

本文提出了一种初至纵波（P波）与瑞雷面波的交叉梯度联合反演策略.通过对初至P波进行全波形反演可以获得近地表P波速度结构;通过对仅含瑞雷面波信息的地震数据转换到频率-波数域进行加窗振幅波形反演（Windowed-Amplitude Waveform Inversion,w-AWI）可获得近地表横波（S波）速度结构.在二者反演的目标函数中均加入P波速度和S波速度的交叉梯度作为正则化约束项,使得在反演过程中P波速度和S波速度相互制约,相互约束,从而实现对地震初至P波与瑞雷面波的联合反演.数值模拟结果表明交叉梯度联合反演可以提高S波速度反演分辨率,而P波速度反演结果并没有得到提高.实际资料的反演结果表明,交叉梯度联合反演能够获得更加可信的近地表速度结构.     

重震联合反演方法及其应用进展

地球物理反演是探索地下结构的最佳途径之一.地震波可以穿透到地球深部进行直接采样,是探测地球深部的主要方法.重力是结构体密度分布与地表观测点之间距离的体积积分效应,重力异常随着源深度的增加衰减很快,其对浅部结构的灵敏度明显优于地震数据.地震和重力联合反演能够相互补充和约束,提高空间分辨率,使反演结果更加稳定可靠.本文首先介绍了联系地震和重力数据的速度-密度经验关系,随后分别介绍了重震联合反演的3种常用方法—顺序反演、同步反演和交叉梯度反演,简要阐述了各种方法在国内外的应用情况.分析认为顺序反演将两类数据分开独立进行计算,原理简单,易于操作实现.但是该方法依赖于先验模型和速度-密度经验关系,可能存在分辨率较低区域(如模型边界)的误差放大效应.同步反演采用将地震和重力数据放在同一方程组中同时反演的策略,减弱了单一数据先验模型对结果的影响,但两种数据的同时运用势必引入数据权重分配问题.交叉梯度寻求不同物理参数模型在结构上的相似性,对潜在的岩石物性关系做了最少的假设,一定程度上降低了反演的非唯一性,但强制性地匹配模型的结构不一定完全符合地下介质的物性分布.因此使用交叉梯度方法反演时应注意模型的推导需要遵循客观标准,以控制模型的结构相似性和数据拟合度.最后指出重震联合反演中的速度-密度经验关系和数据的权重分配仍是值得探究的问题.     

部分区域约束下的交叉梯度多重地球物理数据联合反演

交叉梯度联合反演方法通过对多种地球物理模型实现结构耦合,在岩石物性关系不确定的情况下,既能提高反演结果的可靠性,又能减少反演的多解性,还能减少不同方法解释结果之间的矛盾.当不同的模型观测数据覆盖范围不一致时,交叉梯度联合反演通常需要取出重叠区域数据进行联合反演,并且建模时还要扩展一些模型范围.本文首先提出并实现了部分区域约束下的交叉梯度多重地球物理数据联合反演算法;接着进行了算法的模型试算;最后,我们将该反演算法用于本溪-集安深部地质调查重磁电综合地质地球物理解释中.结果表明：该算法不但能在重叠区域内很好地恢复结构相似的模型,而且在非重叠区域与重叠区域的边界处仍然可以得到平滑变化的模型;在本溪-集安10号剖面所获得的结构上相似的电阻率、密度及磁化率模型较好地反映了该区的深部地质结构,对于确定深部地质体的性质提供了有力的证据.     

基于交叉梯度约束的电阻率法和背景噪声法三维联合反演研究

电阻率法和背景噪声法是通过获得地下介质的电阻率和速度参数的分布来探究地球内部物质分布的非均匀性特征,联合反演可以有效地发挥两种方法的优势,减小反演多解性.本文采用有限内存拟牛顿（L-BFGS）算法实现了电阻率法和背景噪声法的单方法三维反演,然后基于电阻率和速度模型结构耦合约束,采用交替迭代方式实现了电阻率法和背景噪声法的三维联合反演.通过设计规则体组合模型和不规则体组合模型进行理论模型合成数据三维反演试算,结果表明：联合反演可以获得同时满足数据拟合差和结构更为相似的速度-电阻率模型;联合反演可以压制背景噪声单方法反演出现的假异常,获得更精细的速度结构;联合反演获得的电阻率模型对倾斜异常体、高阻覆层下方异常体、圈闭的高/低阻体等边界信息有明显的提升,有效克服电阻率法单方法反演的局限,提高深部电阻率的分辨率.

     

基于交叉梯度约束的电阻率法和背景噪声法三维联合反演研究

电阻率法和背景噪声法是通过获得地下介质的电阻率和速度参数的分布来探究地球内部物质分布的非均匀性特征,联合反演可以有效地发挥两种方法的优势,减小反演多解性.本文采用有限内存拟牛顿（L-BFGS）算法实现了电阻率法和背景噪声法的单方法三维反演,然后基于电阻率和速度模型结构耦合约束,采用交替迭代方式实现了电阻率法和背景噪声法的三维联合反演.通过设计规则体组合模型和不规则体组合模型进行理论模型合成数据三维反演试算,结果表明：联合反演可以获得同时满足数据拟合差和结构更为相似的速度-电阻率模型;联合反演可以压制背景噪声单方法反演出现的假异常,获得更精细的速度结构;联合反演获得的电阻率模型对倾斜异常体、高阻覆层下方异常体、圈闭的高/低阻体等边界信息有明显的提升,有效克服电阻率法单方法反演的局限,提高深部电阻率的分辨率.     

基于交叉梯度约束的时间域激发极化法二维同步反演

时间域激发极化法(Time-domain induced polarization method,简称为TDIP)已有的反演算法采用的是分步反演的思路,即先由视电阻率资料反演电阻率,固定电阻率再由视极化率资料反演极化率,这样就存在极化率结果严重依赖于电阻率反演结果的问题.为了有效解决这一问题,本文实现了TDIP二维数据空间分步反演算法,提出了基于交叉梯度约束的TDIP二维同步反演策略,实现了交叉梯度约束的电阻率和极化率二维同步反演算法.分别用电阻率和极化率结构一致和不一致的二维模型合成数据进行了分步和同步反演试算,对不同模型试算结果进行了对比分析.结果表明:对于电阻率和极化率结构一致和不一致模型,同步反演结果比分步反演结果能更好地确定异常体的空间分布范围,反演得到的电阻率和极化率值更接近真值.理论模型算例表明本文提出的同步反演算法有效解决了分步反演的问题,优于分步反演算法,具有更好的实用性.     

Two-dimensional data-space joint inversion of magnetotelluric,gravity, magnetic and seismic data with cross-gradient constraints

In recent years, joint inversion has been widely used for integrated geological interpretation. We extended a data-space joint inversion algorithm of magnetotelluric, gravity and magnetic data to include first-arrival seismic travel-time and normalized cross-gradient constraints. We describe the main features of the algorithm and apply it to synthetic data generated for hypothetical models. For the synthetic data, we find that the joint inversion with multiple parameters is superior to the joint inversion with two or even three parameters, which can reduce the multisolution of inversion results more effectively. Furthermore, data-space joint inversion involves fewer memory requirements and better calculation speeds than traditional model-space joint inversion. The normalized cross-gradient constraints can better couple model parameters of different magnitudes compared with traditional unnormalized cross-gradient constraints, resulting in higher levels of structural similarity among resistivity, density, magnetic susceptibility and velocity models.     

大地电磁、重力、磁法和地震初至波走时的交叉梯度二维联合反演研究

为了准确的探测和描绘地下复杂的地质结构,同时克服地球物理单一方法反演的多解性和单一参数反演模型的不一致性等问题,近年来基于交叉梯度联合反演的综合地球物理解释已经得到了广泛的关注和应用.本文首先研究了两种地球物理方法的交叉梯度联合反演算法,在此基础上,推导并实现了多种地球物理方法（大地电磁,重力,磁法,地震初至波走时）的多交叉梯度约束的二维联合反演算法;其次,我们设计了结构不一致模型和复杂模型,针对多物性联合反演算法的准确性和有效性进行了模拟试算,并对复杂模型的单独反演结果和联合反演结果进行了交叉梯度值和物性交会图的对比;最后,本文将成熟的卫星资料多光谱综合分析技术应用到联合反演中,将多物性参数反演模型结果图通过RGB（红-绿-蓝）模式进行合成,得到融合的RGB合成图.结果表明：通过对结构不一致模型和复杂模型的联合反演结果和单独反演结果的对比分析,可以得出联合反演得到的结果更接近真实模型,并从得到的交叉梯度值进一步证明了联合反演模型相似度高,也从物性交会图中得到联合反演的物性相关性更好的结论,反向证明了算法的正确性.最终从得到的RGB合成图像,我们可以更直观的分析反演结果,更有利于准确划分地下模型结构.