磁性条件下频率域航空电磁正演计算和成像研究

频率域航空电磁响应计算中,如何计算积分方程是频率域航空电磁法正演计算的难点,本文把电磁响应计算公式分解成汉克尔变换式和普通积分或者超几何函数的两项之和,然后分析核函数中反射系数的衰减规律,通过对比线性滤波算法和直接数值积分法的积分效果,选择了120点线性滤波算法计算汉克尔变换,采用高密度方法在确保计算精度的条件下显著提高了频率域航空电磁响应计算效率.针对强磁性区域频率域航空电磁法反演问题,在传统相位矢量图只包含电阻率和飞行高度两个参数的基础上,添加了磁导率的计算,绘制了包含磁导率、电阻率、飞机飞行高度三个参数的三维相位矢量图,计算结果表明,三维相位矢量图插值可以用于简单的一维模型,反演速度快、效率高,有助于提高磁性区域地下目标体的探测准确率,对于强磁性区域航空电磁法视电阻率填图是一种行之有效的勘探方法.     

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

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

利用数值积分提高一维模型电偶源电磁测深响应计算精度

讨论了直接利用数值积分提高电偶源电磁测深响应计算精度的方法．具体为对Hankel积分进行直接积分，结合连分式展开方法以提高积分求和的收敛速度．利用该方法对均匀半空间和层状(两层)模型的电磁测深响应进行了模拟，结果表明与常规的快速Hankel滤波方法相比，采用直接数值积分能明显提高电偶源频率测深响应计算精度．从而为获得高精度瞬变测深晚期响应提供算法基础．     

用有限元法正演京津唐地区二维构造的大地电磁测深资料

本文阐述在大地电磁测深工作中,用有限单元法,模拟计算京津唐地区二维地质构造的大地电磁响应。 文章前部分,论述了用有限单元法正演二维大地电磁测深曲线的原理及其数学公式推导。文章后部分,根据两条测线实测的大地电磁测深结果并结合其它地球物理探测资料,给出了两个二维剖面模型,并计算了两种极化的视电阻率曲线。两条实测剖面分别是横切唐山地震带的马兰峪—柏各庄剖面及横过北京地区北北东构造走向的文化营—王辛庄剖面     

基于快速拟线性积分方程法的大地电磁三维正演研究

三维大地电磁测深正演所用的时间长与占用的内存大的缺陷是如今国内外学者一直致力研究的问题.本文首先开展基于拟线性积分方程的大地电磁三维正演研究,通过压缩格林系数矩阵减少了庞大的计算量,达到了快速生成系数矩阵的目的,实现了快速拟线性近似的方法.然后,本文设计了不同背景的理论模型,对理论模型进行了快速拟线性积分方程法和积分方程法的正演模拟试算.正演模拟结果表明:采用快速拟线性积分方程法的正演视电阻率曲线与积分方程法求得的视电阻率曲线结果基本相同,表现出较高的准确性,并且快速拟线性积分方程法具有更快的计算速度,特别当地下剖分块数增大时,正演计算效率更为明显.     

地球曲率对长周期大地电磁测深法的影响

基于电磁理论,导出了层状球体大地条件下长周期大地电磁响应的计算公式,编制了计算机程序,并对相应理论模型进行正演计算。根据计算结果,分析了地球曲率对长周期大地电磁测深法视电阻率和阻抗相位的影响。理论模型的计算结果显示,随着探测周期的增大,地球曲率对长周期大地电磁响应的影响逐渐增大,视电阻率曲线会明显下降,阻抗相位曲线逐渐趋于90°;在一维情况下,对于目前的仪器测量范围可以忽略地球曲率带来的影响,但随着仪器技术的进一步发展,地球曲率对长周期大地电磁测深法的影响不可忽视。     

电偶源频率电磁测深视电阻率多参量联合反演效果检验

通过大量的数值模拟试验 ,对电偶源频率电磁测深 (FEMS)视电阻率多参量联合反演的效果进行了检验和对比 .结果表明 ,多参量联合反演具有较快的收敛速度和较高的拟合精度     

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

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

长偏移距瞬变电磁测深甚晚期响应及视电阻率的数值计算

文中讨论了提高长偏移距瞬变电磁测深甚晚期响应计算精度及相应的全区视电阻率的计算方法。在频率域 ,利用直接数值积分结合连分式展开提高低频段响应的计算精度 ,从而提高时间域甚晚期响应的计算精度。同时 ,利用连分式渐进展开 ,将均匀半空间电阻率表示为该模型长偏移距瞬变电磁测深响应的反函数 ,利用迭代方法求出适合全时间段的全区视电阻率。数值结果表明 ,文中的方法可有效地提高长偏移距瞬变电磁测深 (甚 )晚期数据的模拟与解释精度     

肥东地震台地电阻率观测场地勘选

按照《地震及前兆数字观测技术规范》和《地震台站观测环境技术要求(第二部分,电磁观测)》(GB/T 19531.2-2004),对安徽省肥东地震台地电阻率观测场地进行了勘选,采用十字电测深法勘探地下介质的电性结构分布情况,并对观测环境进行电磁骚扰测试,经过对测量结果进行计算和初步分析处理后,认为该场地适合地电阻率台站建设...     

直流电阻率测深中对转换函数曲线的直接反演法

根据三极和对称四极直流电阻率测深的转换函数曲线性质,用大地电磁测深中BOSTICK相同的方法,导出了对电阻率转换函数曲线的直接反演公式。经过对6个模型的电阻率转换函数曲线进行反演检验,该公式虽能近似反映地电断面情况,但误差较大。本文对视电阻率测深曲线的直接反演公式进行改进,得到了另外两种精度较好的直接反演公式。同时给出了对一些例子。     

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

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

瞬变电磁测深早期数据的修正

通过时间-频率转换关系,TEM数据可以转换成平面波场测深数据,从而可以对TEM资料进行拟平面波场处理解释.在对瞬变电磁视电阻率数据向平面波场测深视电阻率数据转换时,发现由于瞬变电磁使用晚期计算公式及装置问题,使测深曲线早期数据发生畸变.文中建立了视电阻率曲线进入晚期所满足的关系式,从理论上给出不同情况下瞬变电磁测深视电阻率曲线进入晚期的临界点.以瞬变电磁、大地电磁、CSAMT为例,对大量的模型进行正演计算,对计算结果进行对比分析,建立不同地表电性结构、不同时间延迟情况下,瞬变电磁早期数据误差的校正量板.     

基于全区视电阻率的线源FCSEM二维正演模拟

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.     

EDA各向异性层状介质MT正演模拟

推导了ＥＤＡ型层状各向异性介质大地电磁测深正演计算公式,对多种模型从不同的频率,测量方向计算了地表视电阻率曲线。结果表明,一般的各向异性介质的地表视电阻率值不仅与测量频率有关,而且与测量方向相对电性主轴的交角有关。     

电阻率测深的数字解释

本文主要介绍了应用积分变换的方法和采样定理将视电阻率ρ_s曲线作线性滤波,得出一新的电阻率转换函数T′曲线,然后,以层参数（各层的电阻率和厚度）算得的T用最优化数值方法在DJS-6型电子计算机上与其进行自动拟合,以达到解释电测深曲线的目的。 文中简述了戈什（Ghosh）提出的对ρ_s作线性滤波的原理,介绍了与国外不同的取样间距和滤波系数的确定以及阻尼最小二乘法和变尺度最优化法的计算框图和应用,最后附有实例和简要的讨论。     

基于遗传算法和统计学的电阻率测深二维反演研究

本文介绍了一种以遗传算法和统计学为基础的反演方法,该方法利用一维反演结果的统计学特征确定地层电性结构并设计初始模型,使用遗传算法修改模型,经过二维有限单元法正演达到拟合误差最小,从而完成电阻率测深数据的二维反演.这个方法的优点在于不需要解大型线性方程组,降低计算量,先判断地层电性结构,反演后给出用地电断面表示的直观的解释结果.理论模型试算和模型实验证明该方法效果较好.     

短偏移距电磁探测技术概述

由于人工源电磁探测效果与收发距离有一定关系,本文首先分析了收发距离为零的中心回线源瞬变电磁法（Transient Electromagnetic Method,TEM）视电阻率定义和算法特点,回溯了从圆形回线到矩形的修正式中心回线装置的发展过程;概述了从中心回线发展到接地导线源短偏移瞬变电磁（Short-Offset TEM,SOTEM）理论上的可行性;列举了将LOTEM（Long-Offset TEM）从数倍于探测深度的观测区推进到近源区后的探测实例.对于频率域人工源电磁法,介绍了起源于MT（Magneto-Telluric）的CSAMT（Controlled Source Audio-frequency Magneto-Telluric）与电磁频率测深法不同发展道路.电磁频率测深将观测区从远区推进到了中区;广域电磁法通过全区视电阻率新定义不仅提高了解释精度,而且将传统CSAMT远场观测模式推进到中区探测模式;最后指出,如能采用新的点微元假设计算方法,可有望分离自有场和辐射场,实现频域电磁近源探测.     

大地电磁测深三维正演模拟及异常特征研究

为了实现大地电磁测深三维高效快速正演算法,本文研究了大地电磁测深积分方程三维正演模拟的理论与方法.首先详细地推导了Fredholm型积分方程;然后给出积分方程的离散方法、并矢格林函数积分的求解方法以及大型矩阵方程的存储与求解方法;最后,通过理论地电模型的计算检验了方法和程序编制的正确性,并展示了几种理论地电模型的数值计算结果,分析其异常特征.算例结果表明,该正演算法在CPU主频为1.81GHz、内存为1GB的计算机上能正确地模拟多个异常体的异常特征;地表电性不均匀体会使地表视电阻率曲线产生十分严重的畸变;TE模式相对TM模式纵向分辨率比较好,而TM模式的横向分辨率较好;随着异常体埋深的增加,反映异常的分辨能力进一步降低.     

基本结构有限元算法及大地电磁测深一维连续介质正演

提出了一种新的数值方法——基本结构有限元法. 从经典的伽辽金法（Galerkin method）出发,首先定义了基本结构插值基函数,在此基础上建立了基本结构方程,然后用有限元法进行进一步的具体实现. 该算法被成功应用到大地电磁测深一维连续介质正演计算中. 采用了6种不同的单元插值基函数进行计算,并对计算结果进行了比较和分析.