电偶源频率电磁测深二维层状模型参数化反演

本文在大地电磁测深松弛（ｒａｐｉｄ）反演理论方法的基础上，探讨了电偶源频率电磁测深（赤道装置）阻抗视电阻率的二维反演问题。由二维层状模型反演目标函数的近似简化出发。将全剖面的二维反演约化为沿剖面各测点的逐一测点一维反演，从而构制了一个简易可行的二维层状模型参数化反演方案。反演试验表明，这一反演方法是实用和有效的     

RRI方法在EH-4数据解释中的应用

RRI(Rapid Relaxation inversion)方法是一种新的MT(Magnetotelluric)反演方法.文章介绍了RRI方法的基本理论,并应用该方法对不同的反演模式和不同的初始模型所得的反演结果分别进行对比,结果表明,二维RRI方法具有快速、稳定,反演效果好等特点.并应用二维RRI方法对内蒙古某矿山EH-4实测数据进行反演,结果表明该方法快速、实用,反演结果符合地质规律,便于地质解释.     

二维视电阻率断面的快速最小二乘反演

快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方     

点源二维地电模型反演方法的研究

为了提高电法勘探数据解释的能力,寻找较好的二维反演方法,本文在研究改进的阻尼最小二乘法在二维反演中应用的同时,成功地引进和发展了Box-Kanemasu反演法。在微机上实现了电法勘探数据的二维反演解释。理论和实际地电模型的二维反演算例表明,该方法的理论是正确的,方法是切实可行的。     

大地电磁测深反演技术有效性对比试验

设计了一个简单的三维模型:三维低阻异常体和高阻异常体位于一维层状介质模型中,以"十"字型和"米"字型观测剖面方式作为模拟方式,开展模型一维、二维、三维反演技术有效性对比试验,其中一维反演计算采用自适应正则化(ARIA)反演,二维反演计算采用非线性共轭梯度(NLCG)反演,三维反演计算采用REBOCC三维反演,将剖面下方不同的反演结果与原始模型进行对比分析,发现:在该模型条件下,一维和二维反演都能得到反映模型真实信息的结果,REBOCC三维反演更偏向重建原始模型的宏观电性特征;在一维反演结果中,对于深部的电性结构特征,通常TE极化模式的反演结果好于TM极化模式,TE/TM几何平均值反演结果介于前两者之间;在二维反演结果中,通常TM极化模式的反演结果好于TE极化模式,TE/TM联合模式反演结果与TM极化模型相当,甚至更好;模型REBOCC三维反演相比一维、二维反演更易受反演结果多解性影响,REBOCC三维反演结果偏向重建原始模型的宏观电性特征。     

基于二维小波分解的MT二维反演方法

对于大地电磁反演,期望得到快速稳定的反演结果并且能够反映清晰的、未发生畸变的地质界面.为了解决这一问题,提出了基于二维小波分解的MT二维反演技术,对小波分解得到的背景场采用光滑反演,细节部分采用聚焦反演.通过模型试验,证明了该方法的有效性,同时对某区的实际资料进行了处理,表明了该方法的适用性和效果.     

直流电阻率测深二维自适应正则化反演

直流电阻率测深二维反演中,正则化参数的选取影响反演结果分辨率及反演过程稳定性。利用主动约束平衡正则化因子,进行直流电阻率光滑约束最小二乘二维自适应反演,改善直流电阻率测深二维反演的分辨率与稳定性。在反演迭代过程中,正则化因子根据模型参数的空间展布函数进行自适应计算、正则化参数的自适应计算。模拟数据反演结果验证了该方法的有效性与可行性,反演结果能准确地反映地下模型的真实电性结构。     

软土地基沉降的二维非线性粘弹性反演与预测

基于Merchant粘弹性模型和简化双曲线模型，建立软土地基沉降的二维非线性粘弹性Biot固结有限元计算公式；根据最小二乘法建立目标函数，采用Simplex直接优化算法反演确定参数；为反演能获得在物理意义上合理且数学上唯一的参数，提出从简单模型到复杂模型的分步反演分析方法。实例分析表明，该方法是可行的。     

二维直流电阻率测深曲线的快速反演

本文采用作者提出的一维直流电阻率测深的直接反演公式,结合二维正演和渐近迭代反演,推出了一种二维直流电阻率测深的快速反演方法,经实际应用,效果很好。     

二连—东乌旗地区固定翼三频航电数据反演方法应用对比研究

固定翼频率域航电测量效率高,空间分辨率高。以往我国三频航电实测数据已获得良好的解释成果,但均以定性解释和一维反演为主,二维反演应用较少。本文将快速近似反演、一维反演、二维反演引入到二连—东乌旗重点航电异常实测数据的解释中,并将三种方法的应用效果进行了对比。快速近似反演以视电阻率计算为基础,一维反演采用Brent方法,二维反演采用数据空间Occam反演,采用Rodi法求解雅克比矩阵。结果显示,快速近似反演可快速有效地提供大数据量地下电性分布信息,一维反演可提供稳定的、具一定深度的反演结果,二维反演则具有明显的横向、纵向分辨率优势。     

利用一维正则化反演进行大地电磁测深数据拟二维反演解释

大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。     

On the application of combined geoelectric weighted inversion in environmental exploration

Geophysical surveying methods are of great importance in environmental exploration. Inversion-based data processing methods are applied for the determination of geometrical and physical parameters of the target model. The use of this geoelectric inversion method is advantageous in environmental research where highly reliable information with large spatial resolution is required. The 2D combined geoelectric inversion (CGI) method performs more accurate parameter estimation than conventional 1D single inversion methods by efficiently decreasing the number of unknowns of the inverse problem (single means that data sets of individual vertical electric sounding stations are inverted separately). The quality improvement in parameter space is demonstrated by comparing the traditional 1D inversion procedure with a 2D series expansion-based inversion technique. The CGI method was further developed by weighting individual direct current geoelectric data sets automatically in order to improve inversion results. The new algorithm was named combined geoelectric weighted inversion, which extracts the solution by a special weighted least squares technique. It is shown that the new inversion methodology is applicable to resolve near-surface structures such as rapidly varying layer boundaries, laterally inhomogeneous formations and pinch-outs.     

二维地层结构和速度的走时反演

介绍了在倾斜均匀各向同性层状介质情况下,根据地震反射走时反演二维地层结构和地层速度的地震射线层析成像方法——走时反演。正演模型射线追踪是根据Fermat原理实现的,即求解满足该原理的非线性方程组得到射线与界面的交点,进而计算相应的走时。反演是先假设一初始模型,用最优化方法使射线追踪走时与观察走时的残差极小。最后计算了分辨矩阵和信息密度矩阵,以评价反演结果。对有噪情形也进行了反演。     

用STRATA进行地震数据的反演

STRATA为一套交互2D/3D模型和反演软件,其功能是将叠后的地震道转换成波阻抗地震道。除了地震反演外,STRATA软件还包括其它一些功能如地震处理,子波提取和井编辑等。所有任务可在多个连接窗口中完成和显示,可以建立地质模型,分析地震数据和交互研究反演结果。在软件中提供了多种反演方法如基于模型反演,稀疏脉冲反演和人工智能神经网络反演。STRATA还提供了分析分步反演结果的工具。     

大地电磁面积性资料和稀疏测线资料的三维反演解释

对理论模型合成的稀疏测线数据和面积性数据以及实测单一测线数据进行了三维反演。通过对比分析合成稀疏测线数据和面积性数据的三维反演结果发现,使用4个张量阻抗元素作为反演数据的稀疏测线数据的三维反演结果可以较好地反映三维目标体,虽然效果可能不如面积性数据的三维反演结果。实测单一测线数据的三维反演结果也较好地揭示测线周围地下基岩界面的三维起伏情况。合成数据和实测资料的三维反演结果表明:可以采用三维反演解释的方法对大地电磁稀疏测线资料进行反演解释。同时,从合成数据的三维反演结果分析可以看出,对大地电磁稀疏测线资料进行三维反演,应尽可能反演所有张量数据;对大地电磁面积性资料进行三维反演,可以只使用两个非对角元素数据,但使用所有张量数据的效果更好。     

ZTEM二维非线性共轭梯度反演研究

ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。     

基于乘积型目标函数的电阻率法三维反演

为了减小地球物理反演的多解性,通常采用累加型目标函数,即数据拟合差项加上正则化项。基于累加型目标函数的反演存在如何优选正则化因子的问题,这个过程通常需要做大量的反演计算。本研究系统分析和实现了基于乘积型目标函数(数据拟合差项和正则化项相乘)的电阻率法三维反演,乘积型目标函数的反演不存在优选正则化因子的问题。三维正演采用非结构化有限单元法,三维反演采用有限内存拟牛顿方法。使用理论模型合成数据进行了三维反演,检验了基于乘积型目标函数的电阻率法三维反演的可行性与有效性。反演算例结果表明：基于乘积型目标函数的反演方法能够可靠恢复异常体的电阻率值、形态和位置。     

基于L—BFGS算法的全波形反演Q值方法研究

估计地层品质因子Q对描述地下构造分布和油气预测都有着重要的意义。考虑到地震波在介质中传播时的黏滞吸收作用,笔者选择2D黏滞性声波方程全波形反演,将观测波场与正演模拟计算得到的波场构建一个目标函数,引入L--BFGS算法反演出地层品质因子。将该方法应用于异常体模型和Marmousi模型测试,结果均表明反演的模型与理论模型的相对误差较小,计算精度高,反演效果好。     

一维层状介质大地电磁模拟退火反演法

大地电磁模拟退火反演法是一种最优化的非线性反演方法,与传统的线性反演方法相比该方法具有：（１）不依赖于初始模型的选择;（２）能寻找全局最小点而不陷入局部极小;（３）在反演过程中不用计算雅可比偏数矩阵等优点;通过对各种类型的大地电磁测深理论曲线试算,结果表明模拟退火法能准确地自动反演地电参数（地层电阻率,厚度）最后对实际资料进行了处理,取得了较好的效果。     

隧道勘探AMT数据二维非线性共轭梯度反演的关键参数探讨

研究了二维非线性共轭梯度反演对于铁路隧道音频大地电磁(AMT)勘探的适应性。通过模型反演实验,确定控制反演进程的关键参数选取方法,以便于在不损失分辨率的前提下避免假异常。实验结果表明,在类似的地电条件下:①采用TE+TM模式联合反演;②使用浅部细分网格以更好地拟合地形;③不使用1000 Hz以上易受干扰的数据;④使用较小的正则化因子;⑤选用中、浅层电阻率的均匀空间作为初始模型;以上这些手段能提升反演结果的准确性,更适应铁路隧道勘探解释需要。