不同模型约束的SOTEM稳定性反演

短偏移距瞬变电磁法（简称SOTEM）理论上探测精度较高,但由于地下电性结构复杂多样,目前所采用的单一模型约束反演策略具有局限性.本文提出了一种不同模型约束的SOTEM稳定性反演方案,能够快速获得具有光滑型、聚焦型、突变型结构特点的反演结果.该方案引入了光滑泛函、总变分泛函、最小梯度支撑泛函进行不同正则化,将不同稳定泛函重新构造,在统一反演框架下获得不同模型约束反演问题的解,并对正演、迭代过程、雅可比矩阵计算、正则化参数设置等多个细节进行了优化,以提高SOTEM反演的稳定性和效率.通过不同地电模型的大量样本,考察了本文所提的多种优化方案对反演的稳定性和效率的提升效果,对比了采用不同模型约束进行正则化的反演效果,从而验证了所提SOTEM反演方案的稳定性与快速有效性.

     

瞬变电磁晚期场资料的一维反演

本文研讨了由垂直磁偶源和水平电偶源所激发的瞬变电磁晚期场的垂直磁场时间偏导数及其视电阻率的性态特征。探讨了其一维反演的某些问题，研制出了对该时间域电磁响应函数的一维实用化反演程序，表明具有较好的使用效果。     

多道瞬变电磁法共中心点道集数据联合反演

多道瞬变电磁法是目前地面资源勘查领域的研究热点,开展该方法观测数据的一维反演对完成装备研发及技术推广具有积极的推动作用.针对多道瞬变电磁法特殊的源波形和数据采集方式,本文对多道瞬变电磁数据的响应提取和一维反演进行了研究.首先,采用反卷积方法从伪随机二进制序列全波形响应数据中提取大地脉冲响应.为了避免脉冲响应中空气波的干扰,提出采用反向积分算法,把脉冲响应转换成下降沿阶跃响应,并将其整理成共中点道集.最后,为充分利用各个偏移距下的响应数据对不同深度范围目标体的分辨能力,采用基于OCCAM算法的联合反演方法对共中心点道集数据进行联合反演.结果表明在不同目标体深度以及有噪声干扰的情况下,共中心点道集数据联合反演计算均可获得较好的结果.     

基于瞬变电磁矩变换的快速三维反演方法

瞬变电磁法的严格三维反演计算复杂、占用资源多,在普通计算机上难以实现.本文引入瞬变电磁矩变换的概念,提出一种快速三维反演方法.该方法基于阻性限制（resistive limit）特性,建立包含异常体的三维大地的一阶矩响应正演算法,根据不同约束条件,选择优化的最速下降法实现瞬变电磁快速三维反演.文中通过含异常体的三维大地正演一阶矩与仿真数据一阶矩的对比,验证了快速三维正演算法的有效性,之后在不同约束条件下,利用优化的最速下降法实现了对含噪声的仿真瞬变电磁数据的快速三维反演.结果表明,该方法能够在普通计算机上短时间内较为准确地反演出地下异常体的体积和位置,在瞬变电磁数据的实时解释工作中具有良好的应用前景.     

具有干预机制的CSAMT数据一维最优化反演

在CSAMT测深中,测深曲线往往受到电偶极源的影响而出现畸变,这给CSAMT数据的反演解释带来很大的困惑。文中依据CSAMT法在波区与MT理论相似的原理,提出用Bostick半定量结果作为CSAMT一维反演的初始模型,采用传统最小二乘算法,结合多种常识性物性变化特点,优化层参数并控制反演过程的迭代方向,使纯粹的数值迭代转化为按地球物理特征自动迭代,其结果可更好地满足物探解释需要。算例和实测数据反演表明,优化后的CSAMT一维反演方法,反演精度较高、计算速度较快,其反演结果的正演响应曲线和实测曲线拟合得很好。     

大地电磁一维连续介质反演的曲线对比法

根据电磁波的传播特性,把视电阻率随周期变化的曲线转化为电阻率随深度变化的曲线,并以此作为初始反演的地电模型．通过初始地电模型得到的视电阻率曲线与真实模型的视电阻率曲线的对比,对初始地电模型的电阻率值进行校正．校正后的反演模型的视电阻率曲线与真实模型的视电阻曲线的拟合程度有所提高．然后如此反复进行多次校正,获得与真实模型更为接近的反演模型,反演的拟合误差一般小于ｌ％．模型试验和实际例子表明,该方法的拟合程度优于Ｂｏｓｔｉｃｋ法．     

大地电磁一维连续介质反演的曲线对比法

根据电磁波的传播特性，把视电阻率随周期变化的曲线转化为电阻率随深度变化的曲线，并以此作为初始反演的地电模型．通过初始地电模型得到的视电阻率曲线与真实模型的视电阻率曲线的对比，对初始地电模型的电阻率值进行校正．校正后的反演模型的视电阻率曲线与真实模型的视电阻曲线的拟合程度有所提高．然后如此反复进行多次校正，获得与真实模型更为接近的反演模型，反演的拟合误差一般小于ｌ％．模型试验和实际例子表明，该方法的拟合程度优于Ｂｏｓｔｉｃｋ法．     

SOTEM响应特性分析与最佳观测区域研究

电性源短偏移距瞬变电磁法(SOTEM)是目前研究和应用较为广泛的一种人工源时间域电磁法工作装置,对深部资源地球物理精细探测具有一定的实际意义.为了深入理解方法内涵并更好地进行推广应用,本文基于电性源瞬变电磁一维正演理论,研究了SOTEM地下感应电流扩散、多分量电磁响应平面分布、多偏移距衰减等特性,然后根据上述特性研究了SOTEM的最佳观测区域.研究结果表明:电性源在地下可以产生水平和垂直两个方向的感应电流.其中,水平感应电流又分为上部水平感应电流和下部水平感应电流(又称作返回电流),水平感应电流的极大值主要集中于发射源附近并垂直向下扩散;垂直感应电流极大值沿与地面呈45°角的方向向下、向外扩散,并且具有较低的振幅和较快的扩散速度.电性源激发的六个方向的电磁场分量都具有一定的探测能力,但是考虑到地面观测的方便性和各分量的传播、分布特点,大多数情况仅利用垂直磁场分量Hz(ðB/ðt)和水平电场分量Ex.其中,Hz仅对低阻目标体敏感,且敏感区域位于赤道向区域,并集中在发射源附近;Ex既对低阻体敏感也对高阻体敏感,对低阻体的敏感区域位于赤道向区域,而对高阻体的敏感区域位于轴向区域,并且敏感区域距发射源的距离与目标体埋深和围岩电性有关.     

航空瞬变电磁法一维正反演研究

优化了航空瞬变电磁法一维正演算法,采用新的汉克尔变换系数,理论上提高了正演的精度.提出了航空瞬变电磁法一维反演算法——模型交替调整反演算法,阐述了该算法的原理和计算方法,编写出反演程序,以已知模型正演响应作为实测数据,对若干典型模型进行了反演计算,取得了较理想的反演效果,与Zohdy法相比,该方法有更高的精度.     

考虑关断时间的地面瞬变电磁三维带地形反演

起伏地形和关断时间对地面瞬变电磁响应影响严重,这给传统基于水平地表模型和理论阶跃波形的瞬变电磁数据解释技术带来很大困难.为此,本文开展考虑起伏地形和关断时间的地面瞬变电磁三维反演算法研究.正演采用基于非结构网格和后退欧拉隐式时间离散格式的时间域有限元算法,快速模拟起伏地表模型瞬变电磁响应.反演采用L-BFGS算法,减少...     

数据空间磁异常模量三维反演

强剩磁的存在通常导致了总磁化强度方向未知,进而影响了磁异常的反演和解释.磁异常模量是一种受磁化方向影响小的转换量,可以在强剩磁条件下通过反演三维磁化强度大小分布来推测场源分布状态.我们提出了一种数据空间磁异常模量反演算法来减少剩磁的影响.与标准的模型空间L2范数正则化反演方法相比,我们的方法有两个优点:一是无需搜索正则化参数(需要反复求解非线性反演问题),因而可以减少计算时间;二是反演结果更加聚焦,深度分辨率更高,我们对此进行了原因分析.通过模型和实测数据测试证明了该算法的有效性和更好的反演效果.     

重力数据3D密度成像中EXTR方法的各参数变化对反演结果的影响(英文)

Tikhonov正则化(TR)方法在重磁数据处理中发挥了重要的作用,本文在研究如何利用Tikhonov正则化方法方法解决重力数据3D反演的同时,深入讨论了可进一步提高拟合误差的Extrapolation Tikhonov正则化方法(EXTR)的原理,并就其参数选择方法及各参数对拟合误差、迭代次数及反演结果的影响进行研究。常密度及变密度组合模型试算结果表明,与TR方法相比,EXTR方法不仅可以达到解释人员设定的先验拟合误差水平,在计算时间及迭代次数相应增加的前提下有更高的拟合精度;同时其反演结果也更加紧致,进一步改善了TR反演结果的发散性;并且其反演数据范围更贴近预设模型参数范围,模型特征与预设模型密度分布吻合较好。     

三维频率域航空电磁反演研究

航空电磁数据的三维解释由于数据量大需要有高效的反演算法作为支撑.本文利用两种目前主流的数值优化技术（非线性共轭梯度和有限内存的BFGS法）实现了三维频率域航空电磁反演,并进一步比较了两种方法的有效性和运算效率.在反演过程中,为了更好地反演异常体的空间位置,模型方差矩阵中的光滑系数在反演起始阶段取值较大;当数据拟合差下降趋于平缓时,再利用较小的光滑因子约束反演过程来实现聚焦和获得精确的反演结果.理论数据反演表明这两种优化策略具有相似的内存需求,但是有限内存的BFGS技术比非线性共轭梯度法在计算时间和模型反演分辨率上具有一定的优越性,因此有限内存BFGS法更适合于求解大规模三维反演问题. 模型试验进一步表明目前主流的迭代法求解技术不适合大规模航空电磁数据反演,未来移动平台多源电磁数据快速正反演可通过引入矩阵分解技术来实现.     

利用大地电磁三维反演方法获得二维剖面附近三维电阻率结构的可行性

大地电磁野外实测数据目前大多为二维剖面数据.如何反演这些二维剖面数据获得较为接近实际地电情况的结果,是多数大地电磁工作者关心的问题.我们通过对理论模型的三维响应进行分析和对合成数据及实测资料的反演结果进行对比研究,讨论了利用三维反演的方法来获得大地电磁二维剖面附近三维电阻率结构的可行性.结果表明:可用三维反演的方法来解释二维剖面数据.对大地电磁二维剖面的张量数据进行三维反演,不仅可以沿剖面获得较好的二维断面结果,还能够得到二维反演所不能获得的剖面附近的三维电阻率结构信息.合成数据的反演算例表明:对二维剖面数据进行三维反演时,对角元素对于圈定剖面附近三维异常体的空间分布具有独特作用,应尽量反演所有的张量元素.     

模拟退火方法在三维速度模型地震波走时反演中的应用

采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。     

大地电磁三维快速松弛反演并行算法研究(英文)

为了进一步提高大地电磁三维快速松弛反演的计算效率,在深入分析大地电磁三维快速松弛反演算法的基础上,结合MPI自身的优越性,确定了并行计算的思路,实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维快速松弛反演并行程序运行结果正确,效率提高明显。此成果在普通微机上实现,推动了大地电磁三维反演技术的实用化,可为其它地球物理三维正反演研究所借鉴。