三维起伏地形条件下航空瞬变电磁响应特征研究

航空瞬变电磁法以其速度快、成本低、通行性好等的优势能够有效的应用于地质地形条件复杂的地区.目前对于航空瞬变电磁法的研究主要基于平坦地形的理想情况,对于地形效应的研究相对较少,然而实际应用中地形不可避免,若忽略地形影响将对资料解释造成较大的误差,从而制约航空电磁方法的进一步发展.本文基于交错网格的时域有限差分方法对三维起伏地形条件下航空瞬变电磁进行正演模拟,在保证算法准确性的前提下给出大量模型算例.以经典地形模型为例,利用所给方法计算三维正演响应,结果显示起伏对于航空瞬变电磁数据有着显著的影响且影响主要集中在早期.而后,以实际地质资料为基础,构建起伏地形条件下包含多个异常体的三维复杂模型,计算了复杂模型的航空瞬变电磁响应,并给出三维全域视电阻率曲线,从而对地形效应的影响有了更加直观的认知.最后,通过大量模型讨论了地形的尺寸参数、电性参数、飞行轨迹与飞行高度等因素变化对于航空瞬变电磁数据的影响情况,并得出有价值的结论.

     

瞬变电磁法正反演问题研究进展

对瞬变电磁法的方法发展概况和仪器研制状况做出了综述性评价.对瞬变电磁法正反演问题的研究成果进行了系统总结.目前的数值模拟正演方法主要有一维滤波系数法,三维积分方程法,二维,三维有限差分法,2.5维有限元法等,主要的反演方法有:一维浮动薄板解释法,人机对话自动反演法,烟圈理论解释法,神经网络反演法,成像类反演等,论述了瞬变电磁法各种计算方法的特点.瞬变电磁法的正反演发展趋势主要是研究三维正反演的计算方法和目标体成像系统.     

瞬变电磁法对隐伏岩溶探测的影响因素研究

发育于地下灰岩地区中的隐伏岩溶给隧道、桥梁以及高速公路等工程的建设与运行带来了极大的安全隐患.城市中常见的路面塌陷情况,多数也是由于岩溶发育导致.利用岩溶与其周围灰岩之间的物性差异,多种地球物理探测方法对隐伏岩溶结构进行有效探测.瞬变电磁法由于其探测效率较高、勘探成本低且对低阻异常体灵敏等特点,近年来在地下隐伏岩溶探测...     

瞬变电磁法探测浅海掩埋未爆物的三维仿真研究

浅海掩埋的未爆物对近岸海洋环境保护、资源开发、工程建设及人民生活等将造成极大隐患,准确探测、识别并清除这些未爆物对于人类社会发展具有重要的现实意义.本文基于中心回线装置瞬变电磁方法,采用时间域矢量有限单元方法对浅海三维瞬变电磁响应进行了正演模拟,设计多个典型三维模型算例考察浅海瞬变电磁响应特征.不同浅海三维地电模型的仿真结果显示目标电磁响应幅值明显、异常清晰,与模型设计一致.本文研究成果初步验证了瞬变电磁方法探测浅海掩埋未爆物的可行性,对于海洋工程勘察与军事防备提供了有价值的理论依据.     

等值反磁通原理瞬变电磁法探测薄层能力

等值反磁通原理瞬变电磁法(OCTEM)的应用日益普遍.为了研究其探测能力,论文基于正演层状模型的瞬变二次磁场时间导数响应,以包含薄层的响应和不包含薄层的背景响应的相对误差作为评价指标,对比分析了不同埋深、不同厚度、不同电阻率薄层以及不同电阻率围岩等情况下的瞬变电磁响应相对异常.结果表明:OCTEM对低阻薄层较对高阻薄层更敏感,薄层与围岩电阻率差异越大,薄层厚度与埋深比越大,越容易被探测到.以响应相对异常值15%作为临界标准探讨了OCTEM方法的探测能力(薄层厚度与可探测埋深的百分比):对于低阻薄层,当与围岩电阻率比值为0.1时,OCTEM探测能力约为1.5%;对于高阻薄层,当与围岩电阻率比值为10时,OCTEM探测能力约为20%.相同情况下常规单磁源瞬变电磁法的探测能力分别为2%和25%,略低于OCTEM对薄层的分辨能力.     

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

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

地形对长偏移距瞬变电磁测深的影响研究(英文)

用基于张量格林函数的体积分方程法对三维异常体进行瞬变电磁响应的正演模拟,首先在频率域内计算电磁场分量的频率域响应,然后利用快速数字滤波技术将计算结果转换到时间域。设计和计算了水平电偶极子源激发下层状水平地层模型背景下的常见地形如山谷、山峰地形的模型,并考察分别把源和接收器放于这些地形中的瞬变电磁场响应,详细分析了这些地形对长偏移距瞬变电磁测深(LOTEM)的影响。结果表明,山谷和山峰地形对LOTEM的结果均有不同程度的影响。当电偶极子源放在山谷谷底时,地形对观测异常场的畸变非常严重;当接收器放在山谷中时,接收器处地形的影响强烈但该影响在空间和时间上只是局部的。总体来讲,不论山峰地形位于何处,其对LOTEM的影响相对较小。当地形处于发射源与接收器之间时,地形对LOTEM的影响非常小,表明在进行LOTEM勘探时,选择发射源的放置比接收器的位置更加重要,野外勘探是尽量把发射源选择在开阔的平坦位置。     

基于吸收边界条件的瞬变电磁法三维矢量有限元快速正演

瞬变电磁法的三维有限元正演通常采用齐次边界条件,为满足该边界条件,需要构建较大尺寸的模型,这降低了正演问题的求解速度.针对该问题,本文采用吸收边界条件代替齐次边界条件,以缩小模型体积,加快正演速度：首先,从时间域麦克斯韦方程组出发,推导了基于库伦规范的矢量势的微分控制方程,结合一阶吸收边界条件推导了相应的的弱形式方程;在此基础上采用一阶四面体矢量单元进行单元分析、Newmark法进行时间离散,实现了瞬变电磁法的快速三维正演.通过均匀半空间模型的解析解,H型地电断面的CR1Dmod解和相应模型有限元解的对比,验证了本文算法的正确性.均匀半空间模型分别采用吸收边界条件和齐次边界条件的正演结果对比表明：吸收边界条件确实可以提高三维正演的精度或者缩小模型尺寸、加快计算速度.

     

半航空瞬变电磁数据小波降噪方法研究

半航空瞬变电磁法野外勘探时,实测信号易受外界电磁噪声影响,导致后续反演效果不佳.常规降噪方法在处理非平稳性和相关性问题有一定局限性,因此,本文研究优化了Birge-Massart中经验系数α,提出一种基于Birge-Massart阈值策略的小波降噪法.通过正演生成理论含噪数据加入不同噪声,经本文方法降噪后信噪比(SNR...     

航空瞬变电磁法三维各向异性数值模拟和响应特征研究

随着瞬变电磁法的快速发展,三维任意各向异性介质的数值模拟成为研究热点.本文从时间域的麦克斯韦方程组出发,采用时域交错采样有限差分法,推导了时域电磁场的时间分步迭代公式,实现了任意各向异性介质的航空瞬变电磁三维正演.设计地电模型,与已有的软件计算结果进行对比,检验了三维正演算法的计算精度.设计了典型的三维各向异性地电模型,改变各向异性参数计算了航空瞬变电磁响应,分析了各向异性参数对航空瞬变电磁响应的影响.开发的正演算法为研究瞬变电磁法各向异性响应特征和三维反演提供了重要技术支撑.     

瞬变电磁系统探测海底电导率的研究

选择在海洋工程与环境地质等方面有应用前景的海洋瞬变电磁法作为研究对象．研究了在高导海水中探测海底电导率的电磁感应理论,研制了共轴水平磁偶极-偶极装置（HRHR）的瞬变探测系统和解释软件．经东海海域实地试验,取得了良好效果,证实了瞬变电磁探测的有效性和方法优势．目前,它可作为海底地质填图的主要手段．经拓展后可用于海底分层和基底深度的探测．     

瞬变电磁系统探测海底电导率的研究

选择在海洋工程与环境地质等方面有应用前景的海洋瞬变电磁法作为研究对象．研究了在高导海水中探测海底电导率的电磁感应理论,研制了共轴水平磁偶极-偶极装置（HRHR）的瞬变探测系统和解释软件．经东海海域实地试验,取得了良好效果,证实了瞬变电磁探测的有效性和方法优势．目前,它可作为海底地质填图的主要手段．经拓展后可用于海底分层和基底深度的探测．     

A FIELD OCCURRENCE OF SIGN REVERSALS WITH THE TRANSIENT ELECTROMAGNETIC METHOD*

Sign reversals have been observed with the one-loop version of the TEM method in an area near Cloncurry, Queensland. This is the only area in which such a response has been recorded in seven years of field surveys by the Australian Bureau of Mineral Resources. The geology of the area consists mainly of pyrrhotitic graphitic shale, which has resistivities as low as 0.1 Ωm and frequency effects up to 30%. The sign reversals could be due to magnetic phenomena, reflections from layers, or complex conductivity effects.     

中心回线瞬变电磁法2.5维有限单元算法

瞬变电磁法的野外资料解释基本上停留在一维水平，本文从生产实践中常用的中心回线法入手，基于电磁场本身的叠加原理，从麦克斯韦方程组出发，导出了中心回线瞬变电磁2.5维二次场（纯异常）的有限单元计算公式. 该算法采用三角形有限元网格，在尽可能拟合地下电性断面的情况下减少有限元网格的节点数和单元数；用选主元的LU分解法求解线性方程组，做到了在移动场源时只需改变右端项；在反傅氏变换中，使用新的波数选取方案，让波数m随时间t滑动. 最后给出的算例表明，该算法直接计算异常场，计算速度快、精度高.     

时间域瞬变电磁法中心方式全程视电阻率的数值计算

给出了一种时间域瞬变电磁法视电阻率的数值计算方法，利用该方法可以容易地求出中心方式的全程视电阻率. 根据中心方式磁场垂直分量时间变化率Bz/T的核函数Y′(Z)的表现特征，以参数Z把整个瞬变过程分为早期阶段（Z>1.6）、早期到晚期的转折点（Z=1.6）和晚期阶段（Z<1.6）. 首先分别得到早期视电阻率和晚期视电阻率的精确值，然后通过转折点构成一条完整的全程视电阻率曲线. 虽然磁场垂直分量Bz的核函数Y(Z)是参数Z的单值函数，但同样存在一个从早期到晚期的转折点Z=1.6，转换点两边仍然可以得到一条早期曲线和一条晚期曲线. 在数值计算中，当迭代步长ΔZ<0.005Z时，视电阻率的相对误差小于0.5髎. 理论模型和实际数据计算表明，与早期和晚期近似值比较，全程视电阻率具有更高的精度和分辨率.     

井下瞬变电磁技术在探测煤矿老空水方面的应用

介绍了全空间瞬变电磁(UWTEM)技术在煤矿井下探测老空水方面的应用,研究了井下巷道有限施工条件下探测不同赋存方位采空区的方法,证明了UWTEM技术在探测巷道掘进头前方、巷道顶底板及巷道侧帮等方向老空水的有效性,总结了井下瞬变电磁技术应用中需注意的若干关键问题。结果表明,UWTEM技术能在井下有效地探测巷道顶底板、侧帮及掘进头前方的老空水位置和范围,为井下水害预报和防治提供可靠依据     

二维地质体的瞬变电磁场响应特征

直接从时间域出发 ,应用时域有限差分方法 (FDTD)分析地下和地面的瞬变响应。由于二维情况便于图形表示 ,采取了线源二维地电模型。通过对二层均匀地电结构和含二维低阻体时地面垂直感生电动势曲线的计算 ,说明瞬变场响应有滞后现象 ,即异常地质体一旦引起瞬变响应 ,就会往后延时 ,并使响应延续较长时间。在实际的工程勘察中 ,曾用 3～ 5ms以后的延迟时间探测过较浅的异常体 ,如 2 0m至 150m深度范围内的地下老窑采空区 ,并取得了较好的探测效果 ,上述数字模拟结果为以往的实践提供了理论依据     

地形对大地电磁测深(MTS)资料的影响

地形影响是山区ＭＴＳ工作的主要问题之一，对此做了理论分析并用均匀半空间表面二维地形模型进行了数值模拟。结果表明，地形影响不仅与ＭＴＳ资料的种类、极化模式、场源电磁波周期及测点位置有关，而且还依赖于不同的地形形态和尺度。此外，还根据悬崖模型计算，给出了关于地形影响范围的估     

APPLICATION OF 3D HIGH PRECISION GRAVITY METHOD TO THE EXPLORATION OF DEEP BURIED HILL

Due to the influence of volcanic rock or weathering crust coverage, the quality of deep seismic data is often not reliable and it is difficult to identify small deep buried hill structure beneath the cover. However, the gravity exploration technology can use the remarkable density differences between the object and the overlying strata to identify this special geological reservoir. Although recently several rounds of seismic exploration have been carried out in Raoyang depression, Huabei Oilfield, North China to determine the existence of Hu 8 north deep buried hill, whether the buried hill really exists or not as well as its scale is still in controversial. In this paper, based on the information of seismic data and formation density, deep processing of 3D high precision ground gravity data has been carried out for 3D forward and inversion computation. The dimensional gravity forward calculation results show that the density of the deep anomaly body forms a relatively low gravity anomaly in the earth's surface. By using the potential processing method of vertical second derivative, and sliding filtering, the residual anomaly is separated from regional or background field, which illustrates the existence of a local high gravity anomaly at Hu 8 north area. According to the amplitude of residual gravity anomaly and formation density difference modeling, through a number of 3D forward calculations and 3D inversion of gravity data, the existence of Hu 8 north buried hill and its possible scale are basically determined. The results prove that the 3D high precision gravity method is effective in determining the deep buried hill structure in case that seismic data is not reliable.     

瞬变电磁法在场地勘探中的研究及应用

本文利用等效代换的方法建立计算瞬变电磁响应的概念及方法,通过视电阻率的计算,进行实测剖面的二维电阻率成像,并利用该成像技术进行场地勘探,获得良好的效果。