地面瞬变电磁法一维烟圈反演技术研究

本文首先进行了瞬变电磁一维正演理论研究,推导出了回线源在均匀水平层状介质中的电磁响应计算公式,采用快速汉克尔变换和正、余弦变换的数值滤波算法以及G-S逆拉普拉斯算法计算了垂向感应电压.区别于传统一维反演一般基于晚期视电阻率,利用二分查找法和平移算法计算得到了全区视电阻率,在全区视电阻率的基础上,根据蒋邦远提出的烟圈反演方法进行了地面瞬变电磁法一维烟圈反演技术研究,最后对理论数据及实测数据进行了一维反演解释,证明了这种烟圈反演技术有助于提高实际反演解释工作的准确性.     

不规则回线瞬变电磁法一维烟圈反演研究

瞬变电磁法中最常用矩形或圆形回线,但有时由于地形限制,只能使用不规则回线发射和接收.如果仍用现有的常规回线理论进行处理解释,得到的结果会存在很大的偏差.本文首先进行了不规则回线瞬变电磁法一维正演理论研究,基于电偶极子源的频率域响应公式,通过沿回线积分和时频转换,推导出不规则回线源在水平层状介质中的时间域响应公式.采用欧拉算法、高斯积分和快速汉克尔变换,计算不规则回线内任意一点处的磁场响应.利用改进的二分搜索法计算全区视电阻率,并在此基础上利用烟圈法反演电阻率和深度.通过四个典型地电模型的正反演计算,表明烟圈法反演能够有效反映地电模型的大致形态,可以用于不规则回线瞬变电磁数据的快速反演解释.     

三维瞬变电磁近似反演

根据烟圈理论,将发射源在地球中感应的电流线看成是随时间向下扩散的烟圈,称之为发射源象．将源象在地表的磁场和源附近测得的响应磁场相对比,确定源象深度,由此用三次样条插值法求出源象扩散速度,并与均匀半空间中烟圈向下扩散的速度进行对比,从而近似估算出随深度变化的电阻率．本文主要讨论圆形回线发射源情形,编制出相应的近似反演程序,并对一维层状理论模型和实测资料进行反演试验．     

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

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

SOTEM数据一维OCCAM反演及其应用于三维模型的效果

本文基于垂直磁场分量研究了SOTEM数据的一维OCCAM反演方法,并将其应用于理论三维数据及野外实测数据的反演.对于大部分一维模型,OCCAM反演可取得较好的反演效果,且反演结果不依赖于偏移距;噪声对SOTEM数据的OCCAM反演具有较大影响,但当信号含噪水平不超过5%时,反演结果仍具有较好的准确性;若浅层存在较厚的低阻层,OCCAM反演结果对下部地层的分辨能力下降,仅能获得具有平均效应的电阻率.将一维算法应用于SOTEM三维数据的反演,会产生较大的误差,尤其是在异常体边缘地带影响最为严重.该影响程度与异常体和背景电阻率之间的差异有关,对于大多数电性近似呈连续变化的真实大地而言,一维OCCAM反演算法仍可获得较好的效果.最后通过陕西某煤田深部富水性调查的实测SOTEM数据反演验证了本文的研究成果.

     

瞬变电磁数据L-PSO反演方法

本文针对粒子群优化算法(PSO)搜索步长短的问题,提出将Lévy flight搜索策略引入其中.利用Lévy flight兼具短距离搜索与偶尔长距离游走的特点,在不降低PSO算法收敛效率的情况下赋予其大步长游走能力,提高了PSO算法搜索全局最优的能力,极大增加了PSO算法飞出局部最优的概率.最后本文将L-PSO算法应用于实际资料的反演中,结果表明L-PSO算法适用于瞬变电磁法的一维反演.

     

地面与半航空瞬变电磁法三维联合反演

半航空瞬变电磁法(SATEM)因其具有工作效率高、探测成本低等优势, 成为了一种发展迅速的新兴地球物理探测技术.基于一维反演的SATEM数据处理方法难以满足复杂地质条件下精细解释的需求.同时, 采用单一的探测方式, 存在探测能力有限、反演多解性强等问题.本文开发了适用于地面瞬变电磁法和SATEM的三维反演及联合反演算法.正演采用基于非结构化四面体网格的矢量有限元法, 可以精确模拟地形和复杂地电结构.反演采用具有较高收敛速度的高斯-牛顿法.通过对理论模型分别进行地面瞬变电磁、半航空瞬变电磁单独反演及其联合反演, 验证了反演算法的有效性.通过对比不同装置的单独反演结果, 以及单独反演与联合反演结果, 分析了不同装置的探测能力, 并展示了联合反演可以有效压制反演多解性、提高反演结果分辨率的优势.

     

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

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

直升机航空瞬变电磁自适应正则化一维反演方法研究

本文以中心回线式直升机航空瞬变电磁法为例对时间域航空电磁资料的一维反演算法进行理论研究.由于时间域航空电磁反演问题是典型的不适定问题,多解性严重且稳定性能较差.为此,论文使用Tikhonov正则化反演方法进行一维解释,可以得到用最光滑模型、最平坦模型和最小构造模型约束条件下的反演结果,而正则化因子是根据各次迭代的数据目...     

瞬变电磁法全空间三维伪谱法模拟

本文根据均匀介质中时间域麦克斯韦方程组,应用伪谱法模拟三维全空间瞬变电磁场的扩散.将电磁场的六分量整理成统一的方程,通过初值条件替代源项将源响应问题转化为初值问题.使用快速傅里叶变换近似空间偏导数,应用切比雪夫多项式近似演化算子来求解方程.通过与全空间均匀介质瞬态线源响应的解析解的比较,验证伪谱法的有效性.计算和分析不同时刻全空间均匀介质中高阻体和低阻体的响应,说明瞬变电磁法对高阻体灵敏度不高,但能较好的分辨低阻体的特点,并反映出时间域伪谱法的优点,即可以连续地计算出电磁场扩散及与低阻体的响应过程.最后讨论伪谱法未来可能应用于矿井瞬变电磁和海底瞬变电磁等领域.     

航空瞬变电磁数据光滑拟二维反演

航空瞬变电磁法适用于浅表地层的大面积探测,由于实测数据量巨大,一维反演作为一种快速高效的方法被广泛应用.在层状地电条件下,横向约束反演能有效的改善传统单点阻尼最小二乘反演结果中界面横向不连续的情况,但其存在反演过程依赖初始模型,且在复杂层状情况下迭代不稳定的问题.本文提出一种基于OCCAM反演并在正则化项中引入横向光滑约束的光滑拟二维反演方法,在保持OCCAM方法对初始模型依赖性小优势的同时,也使得相邻测点电性横向连续.在反演过程中采用令满足迭代误差的测点不再参与后续迭代的策略,以此减少不必要的正演和雅可比计算,提高反演效率.通过对三个不同地电情况的理论模拟数据进行光滑拟二维反演,并与其他传统反演方法对比,验证了该方法能有效的提高反演稳定性,得到与理论模型高度吻合的反演结果.     

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

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

大地电磁法三维共轭梯度反演研究

Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjugate gradient inversion algorithm doesn＇ t need to compute and store the Jacobian matrix but directly updates the model from the computation of the Jacobian matrix. Requiring only one forward and four pseudo-forward modeling applications per frequency to produce the model update at each iteration, this algorithm efficiently reduces the computation of the inversion. From a trial inversion with synthetic magnetotelluric data, the validity and stability of the 3D conjugate gradient inversion algorithm is verified.     

回线源瞬变电磁法的三维正演研究

本文利用矢量有限元法对回线源激发的瞬变电磁场进行了三维数值模拟.根据瞬变电磁场扩散的“烟圈”原理,H型和K型地电模型可近似看作埋置于均匀半空间中的水平板状体地电模型,称为准H型或准K型地电模型.在准H型地电模型中,采用了3种不同的网格剖分方案;从方案1至方案3,水平低阻板状体范围不断扩大,逐渐接近于H型地电模型.中心回线瞬变电磁法的数值模拟结果表明:在12.8 ms以内,不同网格剖分方案矢量有限元法解与一维解析解相对误差均在5％以内;在19.2 ms时,三种网格剖分方案的相对误差依次为12％、5％和2％左右.再利用准H型和准K型地电模型对大回线源激发瞬变电磁场的数值结果进行了验证,结果表明:两种地电模型矢量有限元法解均满足精度要求,但准K型地电模型数值精度更高.最后对直立低阻薄板地电模型进行了数值模拟,取得了与物理模拟一致的感应电动势曲线形状.     

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

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

圆锥型场源瞬变电磁数据AWPSO算法优化反演

瞬变电磁圆锥型场源装置有效的减小了线圈间的电感,提高了小装置探测地下浅层的分辨率,但常规反演方法需给定初始模型且反演精度不高.针对瞬变电磁法反演计算问题,通过对粒子群优化算法(PSO)和神经网络算法(BP)分析研究,改进了一种基于神经网络算法Sigmoid函数的自适应加权粒子群优化(AWPSO)算法.采用标准测试函数对算法进行试算,建立多个理论层状地质模型对该算法进行理论验证,最后在地质资料已知地区开展现场实验.计算结果表明,新提出的算法具有更高的全局搜索寻优能力和收敛速度快、计算精度高,且不需要给定初始模型;实验结果显示实测数据反演结果与高密度电法探测结果吻合,证明该算法能够对瞬变电磁探测数据进行反演计算且精度较高,可以在同类型的浅层探测任务中提供参考.     

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

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

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

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

大地电磁法三维快速松弛反演

实现大地电磁法快速三维反演的关键在于找到快速计算灵敏度矩阵的方案. 本文在对大地电磁三维张量阻抗表达式进行深入分析的基础上，获得了三维快速松弛反演算法的灵敏度函数表达式，解决了三维快速松弛反演的核心问题；为了减小反演解的多解性，针对三维问题定义了最小构造函数，实现了求最小构造的三维快速松弛反演算法. 对二维棱柱体、三维棱柱体理论模型的大地电磁合成数据进行了反演试算，反演结果与理论模型相吻合. 对日本Kayabe地区实测资料进行了XY模式反演，反演得到的地电模型较好地反映了地热资源的分布. 反演试算结果表明大地电磁三维快速松弛反演具有计算速度快、结果稳定可靠等特点；使三维快速松弛反演在普通微机上得以实现.     

浅水区的瞬变电磁法：一维数值模拟结果分析

在浅水环境中,利用频率域海洋可控源电磁法探测油气藏薄层面临着挑战,这是因为空气波支配着电磁响应而它只含有少量海底电阻率结构的信息.本文研究了空气波对时间域可控源电磁(CSEM)响应的影 响. 空气波到达时间取决于海水深度.在浅水区空气波到达早,而在深海水区空气波到达晚,空气波与来自深部电阻体的电磁波出现在不同时间段.虽然在中等深度的海水区空气波和来自深部电阻体的电磁信号几乎同时到达,但与不含有油气藏的背景模型相比,仍可以看到明显的异常.浅水区勘探的好处在于可以使用表面拖曳系统.