对大地电磁资料曲线对比反演法的一点改进

对大地电磁资料曲线对比反演法的一点改进晋光文孙洁王继军（国家地震局地质研究所，北京１０００２９）在大地电磁测深资料反演解释中，Ｂｏｓｔｉｃｋ反演法因其不需要初始模型以及实施的方便、快捷而被广泛应图１在Ｓ５２（上）和Ｓ６４（下）号测点，两种反演方法得到     

一种并行的大地电磁场非线性共轭梯度三维反演方法

本文改进并验证了大地电磁测深数据的三维反演算法和并行计算程序,程序对计算机物理内存和CPU速度及数量要求较低,使普通家用机进行三维反演计算成为可能.本文在Newman和Alumbaugh(2000)提出的三维非线性共轭梯度算法和Rodi和Mackie(2001)给出的大地电磁场二维NLCG反演预处理方法的基础上实现了大地电磁场NLCG三维反演算法,改进了的预处理方法,将反演计算对初始模型的依赖性降到最低,并且通过理论模型验证了程序的正确性,并根据日本KAYABE地区实测数据的反演结果验证了算法的实用性.     

实数编码遗传算法在大地电磁测深二维反演中的应用

作为全局非线性优化的新方法之一的遗传算法,近年来已从生物工程流行到大地电磁测深资料解释中.然而,大地电磁反演问题具有不适定性,解的非唯一性.通过结合求解不适定问题的Tikhonov正则化方法,本文采用实数编码遗传算法求解大地电磁二维反演问题.此算法在构建目标函数时引入正则化的思想,利用遗传算法求解最优化问题.常规的基于局部线性化的最优化反演方法易使解陷入局部极小值,而且严重的依赖初始模型的选择.与传统线性化的迭代反演方法相比,实数编码遗传算法能够克服传统方法的不足且能获得更好的反演结果.通过对大地电磁测深理论模型进行计算,结果表明:该算法具有收敛速度快、解的精度高和避免出现早熟等优点,可用于大地电磁资料解释.     

大地电磁阻抗张量模的拟地震解释法研究及应用

阻抗张量模具有多维信息综合和降维的优势,适用于基于严格意义上的一维地电断面情况下拟地震反演法,相比单个阻抗要素反演,可更多地凸显地质体在横向和纵向上的电性分界面,同时又具有最小构造的特性.因此,反演所得反射系数时间剖面不仅可用于判断极化模式和给定二维电阻率反演的初始模型,还可用于优选反演结果模型;在解释时可弥补常规大地电磁资料反演所得电阻率的模糊像,对划分电性层,推断隐伏断层的位置等都具有重要的实际应用意义.并将该方法用于中上扬子油气勘探大地电磁资料处理中,取得了较好的效果,方法具有一定实用价值和应用前景.     

大地电磁测深法中静态效应及其反演

大地电磁测深法的观测结果常受静态效应的影响而畸变,易将静态效应推断为陡立的深大断裂或垂向大延深的异常体,因此,静态效应校正是大地电磁资料处理必不可少的环节.在全面了解国内外静校正研究现状的基础上,用有限元二维正演程序计算了均匀半空间、均匀层状介质及较复杂的地电模型的TE,TM模式响应,得出静态效应在不同地电模型中的响应特征.并将静位移作为一个反演参数直接参与反演,在反演过程中同时对表层局部不均匀体引起的静态位移作校正.结果表明静态效应直接反演法对静态效应的压制具有较好的效果.     

大地电磁资料精细处理和二维反演解释技术研究(三)——构建二维反演初始模型的印模法

在主流的线性最优化大地电磁二维反演中,如何合理构建初始模型是一个亟待解决的问题.常用的是采用均匀半空间或一维反演结果构建初始模型,不易获得稳定可靠的反演效果.实践表明,尽管基于不同初始模型的大地电磁二维反演结果差别较大,但均较初始模型更为接近真实模型.基于这样一种认识,经过反复的理论和实践探索,我们提出构建大地电磁二维反演初始模型的印模法.印模法的基本思想是依据已有反演结果和均匀半空间模型之间的加权来确定下一步二维反演的初始模型,它一方面保留了已有反演结果中关于真实模型的宏观轮廓信息,另一方面,保证了深部电性结构的均匀性,从而满足大地电磁二维正演所要求的底边界条件.基于印模法,本文进一步提出了迭代重构的反演思想.通过多个理论模型和实测数据的反演计算,验证了上述方法可在很大程度上压制初始模型对反演结果的影响.     

瞬变电磁测深资料的ABC算法反演研究

目前瞬变电磁测深资料的反演基本局限于一维,且多数是基于最小二乘法,实际上是将非线性问题做线性化处理,再利用迭代的方法逐步逼近最优解,所以在一定程度上会丢失许多细节信息,同时极度依赖初始模型,易陷入局部最优解,而非全局最优解.为提高瞬变电磁测深法的反演解释精度,本文以瞬变电磁测深的一维正演理论为基础,同时针对野外实测资料的特点,建立非线性反演问题的目标函数,考虑到现有非线性反演算法存在的易陷入局部优化和过学习等问题,首次引入了一种当前优化领域的十分新颖且性能优良的智能优化算法——ABC(人工蜂群)算法,来求解反演问题.ABC算法对求解多极值非线性问题优势明显,它不需要了解待求解问题特殊信息,只需要对问题进行优劣的比较,通过单独人工蜂的局部寻优的行为,最终将全局最优值突显出来,具有很快的收敛速度.通过在Matlab平台下对算法参数的详细试算研究,最终确定了一套最佳反演参数,并以二、三、四层的地电模型为例进行了验证.在野外实测资料的反演解释中,经过与实地钻孔资料及其他反演方法的对比分析,表明利用ABC算法建立的瞬变电磁测深一维反演系统具有明显的理论优势和实用价值.     

复杂条件下大地电磁资料二维反演及初始模型构建研究

大地电磁资料反演方法较多,但构建二维反演初始模型常规方法主要有两种:以均匀半空间为初始模型、直接以一维反演结果为初始模型.本文利用典型实测剖面资料对这两种方法进行反复试验对比,结果表明:在地形及构造复杂地区,由于实测资料存在严重的静态效应,利用常规方法构建反演初始模型进行二维反演,难以达到理想的反演效果.基于此,本文提出了一种构建二维反演初始模型与反演新的思路方法,其核心为用于反演与构建初始模型数据为同一基础数据,但对基础数据分别进行处理,得到两个处理及校正程度不同的数据,分别用于初始模型构建与二维反演输入数据,这样得到的反演结果既能较好地保留浅部至深部有用异常信息,又能够最大程度上消除实测资料中的静态效应.本文进一步对实测剖面反演结果与已知地质信息进行对比,验证了本方法的有效性和可靠性,认为此方法为地形及地质条件复杂地区大地电磁资料处理反演提供了新的思路,具有一定的参考意义和借鉴作用.     

大地电磁测深资料的二次函数逼近非线性反演

将二次函数逼近非线性优化首次应用于大地电磁测深反演问题,该反演方法利用二次函数有唯一最小值的特点进行逼近大地电磁反演模型,从而避免了常规的迭代反演过程中陷入局部极小问题,实现了对目标函数求全局极小,较好地解决了非唯一性问题;同时该方法不用求灵敏度矩阵,且对初始模型无任何要求。通过理论模型检验、井旁MT点反演结果与测井曲线的对比及MT测线的反演电阻率深度剖面与地震测线的时间剖面对比均表明,本文方法取得较好的应用效果。     

大地电磁测深在西安凹陷基底探测中的应用

为了研究西安凹陷深部地层及基底构造特征,采用大地电磁测深十字剖面法对西安凹陷的电性结构进行定性分析,并使用非线性共轭梯度(NLCG)二维反演方法结合地震勘探结果对其定量解释,划分出西安凹陷的主要地层及埋深.大地电磁测深反演结果与地震勘探解释结果相吻合,有效反映了西安凹陷深部地层结构的真实形态,为研究西安凹陷基底结构特征提供有效的依据.     

海洋三维大地电磁的两步反演法:模型测试

大地电磁(MT)响应的反演已被广泛用于探索地球内部的电导率分布.因为大多数三维反演方法的反演结果都会受到初始模型的影响,因此在三维大地电磁反演中,使用好的初始模型非常重要.尽管均匀初始模型广泛用于三维反演,但一维的初始模型也是个不错的选择,因为一维模型较均匀模型可以更恰当地描绘地球内部的实际电导率变化.本研究提出了一种分两步进行的三维反演方法,该反演方法尤其适用于海洋大地电磁研究.该方法首先反演得到一维初始模型,进而依此反演三维电导率结构,并且这两次反演中均考虑了表层三维地形变化.该方法测试了含有不同地形变化(凹陷、平滑变化、槽形和方形高原地形)和不同三维棋盘类型异常(尖锐或平滑横向电导率变化)的三维电导率模型.文章比较了基于所提出方法获得初始模型和基于地形校正方法得到的初始模型,还对基于本研究方法得到的三维反演结果与基于不同的均匀初始模型得到的反演结果进行了比较.利用新引入的模型矢量参数的方向信息讨论了这些反演的结果,验证了该方法的性能和有效性.文章还揭示了新引入的模型向量参数可用于显示反演的收敛性,将有助于选择合适的反演模型.     

二维地电模型多参数反演方法及其在MT测深中的应用

二维地电模型电（电磁）测深法参数化反演,宜采用多参数反演方案.文中导出了多参数反演的基本数学模型,并进行了分解处理,使之变成与剖面上反演点数相等的多个相互独立的小数学模型的组合,用广义逆法逐个求解,即可得一次迭代中整个二维模型参数的修改量.大地电磁测深理论模型和实测数据反演结果表明所提出的二维多参数反演方法可行.它具有稳定的收敛性、较快的收敛速度、能适应较复杂构造的反演和计算量小等特点.     

大地电磁二维反演方法分析对比

基于地球物理正则化反演理论,详细给出了目前常用的三种大地电磁(MT)二维反演方法——奥可姆法(OCCAM)、简化基奥可姆法(REBOCC)和非线性共轭梯度法(NLCG)的理论基础.设计了两个理论地电模型,一个是高、低阻体组合模型,另一个是断裂模型,使用这三种MT二维反演方法分别对TE极化模式、TM极化模式以及TETM极化模式响应数据进行了反演,对反演结果进行对比分析,考察了不同的反演方法以及不同的极化模式对反演结果的影响.研究结果表明:三种反演方法有各自的特点,使用不同的极化模式数据反演将会得到不同的反演效果.     

大地电磁非线性共轭梯度拟三维反演

提出了非线性共轭梯度法大地电磁拟三维反演.该方法选取共轭梯度反演算法为拟三维反演的核心.在计算灵敏度（Jacobian）矩阵时,吸取近似灵敏度矩阵思想,采用一维灵敏度矩阵来代替三维灵敏度矩阵,并对非测点的灵敏度元素提出一种近似方法.在第一次反演之后,采用拟牛顿法更新灵敏度矩阵.拟三维反演法在很大程度上节省了计算时间,并且理论模型和实际资料的反演试算结果表明大地电磁拟三维反演法具有一定的实用价值.     

倾子资料三维共轭梯度反演研究

在对倾子响应和共轭梯度算法深入分析的基础上,我们实现了倾子资料三维共轭梯度反演算法.基于倾子资料的三维共轭梯度反演研究,探讨了利用倾子资料进行三维反演定量解释的方法.通过对理论模型合成数据 进行反演试算,验证了所实现的倾子资料三维共轭梯度反演算法的有效性和稳定性.该反演算法可用于对大地电磁测深和地磁测深(地震地磁台站进...     

网格剖分方式对大地电磁数据反演的影响

随着计算机技术的飞速发展,大地电磁测深数据的各种反演方法得到空前发展。对于反演过程中的网格剖分问题,前人研究有限。本文就此进行研究,确定影响网格剖分方式的因素,找到一种合适的网格剖分方法。为此,设计一个复杂模型,对正演模拟得到的数据加入随机噪声作为反演的输入数据,选择不同的网格剖分方式,使用DASOCC与NLCG反演方法,以相同参数分别进行反演,对所得结果进行对比分析。结果显示,在设定的模型及参数条件下,针对DASOCC反演方法,纵向及横向网格均不加密,反演效果比较理想,即横向采用一个测点对应一个网格的剖分方式,纵向采用表层加密,100 m以下后一网格为前一网格厚度1.1倍的剖分方式;针对NLCG反演方法,纵向与DASOCC反演方法相同,采用不加密网格的剖分方式得到的结果较好,而横向则建议采用一个测点对应一个网格,并在2个测点之间插入一个网格的剖分方式。     

尖锐边界反演在临江市六道沟可控源音频大地电磁测深数据反演中的应用

在实际大地电磁测深数据反演过程中,为了得到更加清晰的电性单元分界面,传统的Occam平滑反演方法遇到了挑战,本文讨论了sharp boundary inversion(SBI)反演的理论及其应用.通过模型试算,对比Occam和SBI的反演效果,证明了SBI反演划分层状介质的可靠性.本文对临江市六道沟地区实测可控源音频大地电磁测深(CSAMT)数据进行了全区视电阻率转换,然后对全区视电阻率进行Occam反演和SBI反演.最后对比已有地质资料,表明SBI反演能够用于CSAMT数据的处理中.     

基于有限差分正演的带地形三维大地电磁反演方法

本研究实现了一套基于有限差分（FD）方法的大地电磁测深数据带地形三维反演算法及代码.其中,在大地电磁场正演数值模拟方面,开发了起伏地形条件下基于交错网格剖分、有限差分方法的大地电磁测深三维正演代码;在满足平面波场假设的前提下,使用长方体网格剖分模拟三维起伏地形,实现了带地形三维正演计算;并设计理论模型进行试算,经试算结果与前人的有限元法计算结果对比,验证了所研发的带地形三维正演计算的正确性与可靠性.在反演方面,本研究基于非线性共轭梯度方法编写了大地电磁测深带地形三维反演代码,试验了不同的共轭梯度搜索因子β,避免了目标函数对海森矩阵（参数二次导数矩阵）的显式计算和存储,初步实现了大地电磁资料的带地形三维反演.最后,对一系列理论模型进行正演计算,利用其生成的合成数据模拟实测数据进行反演,并与现有的不带地形大地电磁测深三维反演结果比较,检验了所研发的带地形三维反演计算的可靠性与稳定性.     

ZTEM三维数据空间OCCAM反演研究

Z轴倾子电磁法(ZTEM)是一种极具前景的新型天然源频率域航空电磁法,具有勘探深度较大,工作效率高的特点.本文基于大地电磁法(MT)三维数据空间OCCAM反演算法,实现了三维ZTEM倾子资料的数据空间反演算法.将该算法应用到理论合成算例中进行反演,反映出ZTEM倾子反演对模型横向边界的约束优于垂向的特性,并且与MT阻抗反演结果进行对比,验证了所实现的倾子资料三维数据空间反演算法的有效性,表明ZTEM倾子反演在约束电阻率模型横向位置能力上更具优越性.该反演算法可用于ZTEM所整理得到的倾子数据进行三维定量反演,实现大区域的地质构造三维成像,获得地下地质模型的电阻率结构信息.     

大地电磁测深和电偶源频率电磁测深资料反演

在研究大地电磁响应函数频散关系的基础上,构制了一套滤波系数算法,以用于由一套视电阻率资料估算相应的阻抗相位.理论模型和实际大地电磁观测资料的数字试验表明,该法是行之有效的.由频散关系估算的相位值与观测的相位资料的比较,可用于检验观测资料是否满足频散关系.利用经频散关系校正的阻抗相位值,进行大地电磁阻抗的联合反演则可望获得更为可靠的结果.研究了观测频带相互衔接的电偶源频率电磁测深和大地电磁测深视电阻率的一维联合反演问题.对两个实测点两种电磁法的观测资料进行了联合反演试验,与钻井资料对比表明,所获得的电性分层参数是较为可信的.在补充了由频散关系获取的电偶源频率电磁相位资料后,对于两种电磁法的视电阻率————阻抗相位、阻抗实部视电阻率-阻抗虚部视电阻率进行了拟大地电磁反演,获得了相近的反演结果.