瞬变电磁虚拟波场二阶Born近似成像算法

依据瞬变电磁扩散场与虚拟波动场之间的数学表达式,通过波场反变换得到虚拟波场,进而借用地震勘探上成熟的成像方法定位地质界面的位置以及形态.本文首先将时间域波动方程转换为频率域波动方程,将总波动场分为背景场和散射场、总波速分为背景速度和扰动速度,然后运用Green定理得到了散射场关于地质界面速度扰动量的表达式,并使用Bor...     

瞬变电磁虚拟波场层状Green函数Born近似成像方法

瞬变电磁虚拟波场Born近似成像是通过Born近似对瞬变电磁虚拟波场方程进行线性化处理, 求取速度扰动量实现成像的方法.Born近似采用半空间Green函数, 在背景速度变化较大时成像分辨率较低.因此, 本文从虚拟波场出发, 建立辐射场与散射场满足的Helmholtz方程, 导出层状介质Green函数表达式, 将其应用于Born近似中, 得到了基于层状介质格林函数的Born近似电导率计算公式和散射场关于速度扰动量的关系式, 最终利用速度扰动量实现成像.通过对理论模型的实验计算, 验证了利用层状模型Green函数的Born近似成像算法, 能够重建地层电阻率分布, 改善Born近似虚拟波动方程逆散射反演的垂向分辨率, 获得较好的成像效果.对甘肃地区的实际资料处理结果表明, 层状介质Green函数的Born近似, 是一个精确度较高、稳定性好、有实际应用价值的成像方法.

     

瞬变电磁扩散场到虚拟波场的精细积分变换算法

现有的瞬变电磁解释与反演方法主要对电阻率参数进行解释,较难给出准确的地质构造信息.为了能够较多获得地质构造信息,可将瞬变电磁扩散场信号转换为虚拟波场信号,提高瞬变电磁法的分辨率.本文采用精细积分法,实现了瞬变电磁扩散场到虚拟波场的转换.将一个高度病态的线性方程组转换成求积分的过程,大大降低了解决病态问题的难度.积分步长...     

瞬变电磁虚拟波场高分辨成像技术研究

本文提出了从瞬变电磁场到波场的优化算法.目的是通过对现有渡场变换过程进行插值的办法提高波场变换的分辨率,从而为后期数据处理的质量打好坚实的基础.在进行波场反变换的过程中,将正则化算法用于逆变换过程的计算中,通过采用偏差原理和Newton迭代格式选出最优的正则化参数,使得反变换所得到的波场稳定、可靠.解决了积分系数过多而...     

多道瞬变电磁法(MTEM)虚拟波场提取技术

多道瞬变电磁法（MTEM）是近年来发展起来的一种新的地球物理勘查技术,其数据采集方式与地震法类似,因此,采用瞬变电磁拟地震解释方法对MTEM数据进行处理解释具有一定的优越性.研究MTEM虚拟波场有效提取方法具有重要意义.在以往奇异值分解法、预条件正则化共轭梯度法两种不同的常规提取方法研究基础上,本文提出采用相关叠加法提取MTEM虚拟波场信息.首先采用预条件正则化共轭梯度法对全时段MTEM数据进行虚拟波场提取,然后采用同种方法对划分的各时间段数据进行虚拟波场提取,最后对全时段提取结果与各个时间段提取结果进行相关性叠加,叠加结果作为最终的提取结果.实测MTEM数据虚拟波场提取结果表明,采用相关叠加法可以得到稳定、光滑的虚拟波场波形曲线,且抗干扰能力强.     

瞬变电磁虚拟波场的三维曲面延拓成像

根据瞬变电磁场扩散方程与波动方程间存在的数学对应形式,将已知瞬变电磁场数据转换为虚拟波场数据,在此基础上,提出了用克希霍夫积分法将所转换的波场从地面向地下反向延拓的方法.在延拓的数值计算中,为提高计算效率,采用三维边界元技术,把边界积分分解为诸三角单元积分的积分之和, 建立了曲面延拓方程式;通过地面测点的波场值,求出地下某一点的波场值,实现了地表为曲面的向下延拓成像计算.通过对理论模型计算和实际资料处理,证明了该方法可以增强瞬变电磁法识别地下电性分界面的能力,使瞬变电磁法对地下目标体的三维精细探测成为可能.     

虚拟波场降速条件下电性源瞬变电磁垂直磁场分量的波场反变换方法

通过波场反变换从瞬变电磁信号中提取虚拟波场是定位电性界面和刻画地下导电目标形态的有效手段.波场反变换过程中, 降低虚拟波场速度有利于改善虚拟波场对电性界面的分辨能力.然而, 波场变换核函数随时间、虚拟时间以及虚拟波场速度变化的动态范围极大, 且随虚拟波场降速而进一步扩大, 造成波场反变换问题高度不适定, 制约了虚拟波场提取的精度和效果.对此, 在虚拟波场降速的条件下, 本文以核函数在无穷区间上积分的解析解为准绳, 通过最优化与数据拟合得到了经验公式, 确定了波场反变换所需的最优积分区间; 在适应采集时间与虚拟波场速度的同时, 最大限度缩短了积分区间长度, 降低了核函数变化的动态范围, 压制了波场反变换方程的病态程度.此外, 通过引入精细积分方法求解波场反变换方程, 提高了波场反变换方程的求解精度.为验证方法的有效性与可靠性, 对正演模拟的电性源瞬变电磁垂直磁场分量及其时间导数分量进行了波场反变换; 结果回代后的拟合数据与正演数据的相对误差小于5%.在所提取的虚拟波场记录中观察到了合理的运动学特征, 如实反映了低阻目标层深度、厚度和二次场观测时间范围的变化.在计算条件相同的前提下, 通过降低虚拟波场速度, 改善了波场反变换的精度和效果.降速后提取的虚拟波场记录对低阻层上下界面的刻画更为清晰.在对实测数据的处理中, 本文的虚拟波场降速、积分区间最优化和精细积分求解方法得到了合理可信的虚拟波场提取结果, 为进一步利用虚拟波场开展低阻目标体成像提供了数据保障.

     

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

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

电性源瞬变电磁有限差分偏移成像方法

目前,瞬变电磁数据的主要处理结果是给出地下介质的（视）电阻率信息,较难给出清晰的地质界面信息.本文基于有限差分偏移成像原理,实现了电性源瞬变电磁虚拟波场的偏移成像技术.根据瞬变电磁扩散场与虚拟波场之间的数学关系,本文使用精细积分法对第一类Fredholm积分方程求解,然后消除直达波,得到只含有地层反射波的虚拟波场.在此基础上,本文从时间域波动方程出发,对上行波的一级近似方程反向外推,得到时间-空间域二维与三维的波场递推公式.通过求解地下空间的波场,并使用反褶积方法削弱虚拟波场的波形展宽效应,从而提高对地质界面的纵向分辨率,获得较为准确的地质界面信息.本文分别采用均匀半空间模型、D型模型、HK型模型以及两层含水采空区模型对本文方法进行验证,并对实测数据进行处理,偏移成像的结果与地质资料相吻合,说明瞬变电磁有限差分偏移成像方法在识别电性界面以及地层特征方面具有可行性.

     

井中磁源瞬变电磁响应特征研究

井中瞬变电磁波勘探是一个全空间地球物理场问题.采用Gaver-Stehfest逆拉氏变换方法,正演计算了瞬变信号激励下接收线圈上的电磁场响应.分析了包含井眼泥浆、套管、水泥环和地层的轴对称多层介质模型的电磁场响应特征,考察了各层介质参数对井中瞬变电磁响应的影响.不同电导率井眼泥浆的电磁场响应衰减曲线表明,井眼泥浆电导率的变化对井中瞬变电磁响应的影响极小.不同套管几何参数的电磁响应数值模拟结果显示,套管内径变化对电磁响应有较大影响;在套管厚度变化的端点处,电磁响应存在明显异常.通过对不同套管磁导率电磁响应特征的讨论,认为套管磁导率参数对电磁响应有重要影响,套管相对磁导率越高,电磁波信号越难穿过套管,响应信号幅度越低.最后,对不同水泥环参数的电磁场响应进行了数值模拟,发现对于低导地层,低阻水泥环会产生较大的测量误差;而对高导地层,水泥环厚度参数对响应信号影响可忽略不计.     

从瞬变电磁场到波场的优化算法

在波场的正变换中,保证瞬变电磁场的计算精度的同时,应用两步最优化算法,成功地控制和减少了积分系数个数和离散数字积分的采样点个数,解决了在波场反变换式中,由于积分系数过多而产生的欠定方程组的问题,同时改善了第一类算子方程的不适定性.在波场反变换中,将正则化算法用于逆变换过程的计算中,通过采用偏差原理和Newton迭代格式选出最优的正则化参数,使得反变换所得到的波场稳定、可靠.通过对数值计算结果与已知波场函数对比,证明了该方法的有效性和实用性.     

矿井瞬变电磁法超前探测合成孔径成像研究

在矿井瞬变电磁超前探测数据处理和解释中,为了更好地识别巷道前方含水体电性分界面,利用扩散场与波动场间的函数关系,将矿井瞬变电磁法超前探测中瞬变电磁波场转换为拟地震波场,实现波场转换.同时,为了解决波场转换后信号幅值偏弱,以增加信噪比利于解释,利用雷达的合成孔径方法对波场转换后的拟地震信号进行相关叠加处理.将矿井瞬变电磁法巷道超前探测的同一测点不同探测方向的数据作为孔径数据进行合成孔径成像,以突出弱异常而提高信噪比,结合应用实例验证了方法的效果.研究表明,矿井瞬变电磁法合成孔径成像方法能够提高电性界面信息的信噪比和界面几何分辨率,对矿井瞬变电磁法精细超前探测含水构造具有重要的理论意义和应用价值.     

瞬变电磁拟地震子波宽度压缩研究

导电媒质中电磁场所满足的扩散方程与无耗媒质中电磁场所满足的无阻尼波动方程之间存在着数学转换关系式,可以实现瞬变电磁法拟地震资料解释.但是转换出的虚拟波往往波形较宽,使分辨能力下降,影响了瞬变电磁拟地震方法的发展和推广.本文针对虚拟波形展宽这一现象进行了研究,认为虚拟波形展宽的原因,不是由于波在介质中传播时能量损耗所致,...     

瞬变电磁测深中数据处理偏差及修正方法研究

瞬变电磁测深方法中通常采用阶跃波激励方式的模型进行处理,但在实际应用中,为了便于重复观测,一般采用方波形式进行激励,而依然采用阶跃波模型进行数据处理和解释.这种处理方式理论上是存在比较大的偏差.为了修正这种偏差,本文通过理论计算,建议采用单脉冲电磁测深方法来处理方波激励源的电磁测深结果.通过实例分析,得知该方法可行.     

瞬变电磁法拟地震成像研究进展

瞬变电磁法是工程地球物理勘察中应用较多的一种勘探方法之一,瞬变电磁法拟地震成像方法研究是当前电磁探测理论与应用研究的热门.文中详细叙述了瞬变电磁法拟地震成像的两种不同思路,一种是建立时间—频率转换经验公式,把TEM数据转换成平面波数据,借鉴大地电磁法拟地震解释思路进行成像,另一种是通过积分变换,把扩散场变成虚拟波场,把电磁衰减信号变成拟地震子波,然后按照地震解释法进行解释,并分析了各自存在的的问题.指出下一步研究的方向.     

从瞬变电磁测深数据到平面电磁波场数据的等效转换

中心回线源装置的瞬变电磁场是一种涡流场,在地下主要以扩散形式传播.平面电磁波在地下介质中传播时具有反射、折射等电磁学特性,为了借鉴平面电磁波场测深成熟的解释方法,文中开展了从瞬变电磁涡流场测深数据到平面波场数据的快速等效转换的研究.通过大量的理论模型的正演计算、曲线对比,误差统计,对瞬变电磁场、平面波场在地下传播机制特性的分析,建立了一种由瞬变电磁测深视电阻率数据向平面波场测深视电阻率数据转换的时间-频率等效对应关系.利用这一关系式,把扩散场数据快速转换成等效平面波场测深数据. 为进一步进行瞬变电磁拟平面波成像解释打下理论基础.     

瞬变电磁合成孔径成像方法研究

瞬变电磁法虽然广泛应用于资源勘探和工程勘察中,但仍存在一些亟待解决问题.迫切需要引进一些新的技术与方法,以进一步提高解释精度.由于瞬变电磁资料可以通过数学公式转换成虚拟波;多孔径瞬变电磁物理模拟证明TEM多孔径具有相干性;相邻位置上同一地质体的反射回波具有较好的相关性.所以,对于瞬变电磁资料可以进行多孔径合成成像.本文借助于合成孔径雷达的基本思想,提出一套新的数据处理方法,即采用相关叠加技术,实现多孔径数据合成.在完成瞬变电磁虚拟波提取后,将传统的以剖面为主的处理方式发展成为以测点为中心的多孔径合成,将传统的以单点处理方式发展成为逐点推移多次覆盖的处理方法.采用相关叠加的方法来进行合成孔径,由此大大提高瞬变电磁法的分辨率.从波场的角度拓展和丰富了瞬变电磁场的内涵,使得从实测资料中提取到常规瞬变电磁法提取不到的信息,对地下目标体成像更为有利.通过对所设计模型和实测资料处理,结果表明合成孔径成像效果较好.研究成果为发展瞬变电磁成像技术,提高分辨率提供了一条新途径.     

瞬变电磁法三维拟地震成像信息提取技术

在瞬变电磁法中,为提高对地下电性界面的分辨率,可移植地震勘探中的诸多技术对界面直接成像,以克服反演中的不适定性.首先为使得波形在高低不同的电性基底所产生的虚拟反射波振幅有可比性,在将原有瞬变电磁场逆变换成虚拟波场数据的基础之上,需要求取均衡系数;接着通过对虚拟波场的离散数据求取反褶积,可基本抵消波场变换的展宽作用和大地滤波影响,以达到对电性界面精确定位的效果;然后,对波形做低通滤波处理后运用克希霍夫积分公式,将虚拟地震波场由地面向地下反向逆时偏移成像.把波场分析原理用于对瞬变电磁场的解释,就形成了瞬变电磁偏移方法.从波场角度看,它拓展和丰富了瞬变电磁场的内涵.使用本方法,可更加容易得到从前用常规解释方法处理无法提取的信息.在延拓成像中可更精确的确定界面深度,并为在瞬变电磁法中合理引入拟地震处理解释方法提供重要的理论依据和实验参照.     

三维瞬变电磁近似反演

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

从瞬变电磁扩散场到拟地震波场的全时域反变换算法

将瞬变电磁满足的扩散方程转变为波动方程,然后利用地震类成像方法实现瞬变电磁虚拟波场成像,是实现瞬变电磁三维反演的有效手段之一.为了实现由扩散场到虚拟波场的转换,文中采用预条件正则化共轭梯度法求解波场反变换问题.首先,对几种离散方式进行比较,采用条件数最小的离散方式进行离散;然后选择最优的正则化参数,并利用超松弛预条件技术对系数矩阵进行预条件处理;最后,利用共轭梯度法进行迭代求解.超松弛预条件有效降低了系数矩阵的条件数,正则化方法使得反变换得到的波场稳定、可靠,共轭梯度法能够保证计算快速收敛.将反变换结果与已知虚拟波场函数对比,证明算法稳定、可信.将文中算法结果与前人研究结果进行对比,说明方法效果.通过实测数据的波场变换处理给出了文中方法的实际应用效果.结合反变换算法,对不同参数模型进行分析,总结了虚拟波场在色散介质中的传播规律.