积水采空区瞬变电磁响应微分解释方法

瞬变电磁(TEM)视电阻率早、晚期计算公式在求取时进行了简化近似,造成部分地电信息的失真.瞬变电磁响应信号衰减过程与介质的电性特征有关,由此引入瞬变电磁响应微分衰减速度和二阶微分衰减速度两个解释参数,力图避开视电阻率计算,以衰减速度表示瞬变电磁探测结果的电性特征.首先研究了不同导电均匀半空间与低阻目标体地电模型瞬变电磁响应微分衰减速度在时间上的变化特征,在此基础上对瞬变电磁响应二阶微分衰减速度进行了分析,通过绘制TEM二阶微分衰减拟断面图实现了低阻目标体地电模型的电性层划分.研究结果表明:电性分界面与二阶微分衰减的"0"等值线相对应,可以用"0"值等值线作为电性分界面的划分依据,即瞬变电磁微分解释方法能够区分不同地电断面,可以用于瞬变电磁探测成果解释.     

从瞬变电磁响应中提取IP信息的研究(英文)

在对浙江省某铜矿进行瞬变电磁勘查时,由于感应激发极化效应的影响,造成瞬变电磁晚期测道的数据发生倒转。采用常规的瞬变电磁数据处理方法难以处理,造成晚期测道数据不可用,影响TEM的探测深度及精度。针对这个情况,采用Cole.Cole模型分析了均匀半空问模型中直流电阻率、充电率、时间常数以及频率相关系数对磁性源瞬变电磁响应的影响规律。利用奇异值分解法(The Singular Value Decomposition,简称SVD)对实测TEM数据进行分析及反演,并且从瞬变电磁响应数据中分离出勘探区瞬变电磁测量数据Cole—Cole模型的各参数,将分离出来的参数用于探测结果的辅助解释。经过分析得出充电率和频率相关系数对瞬变电磁响应的影响较大,而直流电阻率和时间常数对其影响较小的结论。通过常规方法和奇异值分解法对实测数据分别进行处理,发现后者视电阻率断面图的异常更加突出,从而较准确的圈定了铜矿体地范围,此结论得到_『钻孔的验证,与实际地质情况吻合较好。     

地-井瞬变电磁多分量响应数值分析（英文）

本文对地-井瞬变电磁法多分量响应进行计算分析。规则局部体瞬变场响应的计算方法与解释模型对于实际导电围岩模型的适用性较差,针对该问题,本文提出了一种基于地下瞬态电磁场数值模拟的计算分析方法。瞬变电磁场模拟方面,本文以时域有限差分法实现正演模拟,引入采用Gaver-Stehfest逆拉氏变换与Prony法的离散镜像法求解初始电磁场,应用透射边界条件保证迭代计算精度。通过均质半空间模型算例,证明该套方法可行。响应分析方面,设定含井旁目标体和导电围岩的地电模型,以上述方法对地下瞬态电场进行正演,以多分量观测装置为例换算感应电动势。通过对比各条件下瞬态电场与多分量响应,得出结论:地-井瞬变电磁多分量感应电动势响应反映了地下瞬态电场沿水平、垂直方向的梯度变化;响应特征取决于地层中瞬变场在不同条件下的"扩散、衰减、畸变"过程和观测位置的电磁场状态。本文的计算分析方法兼顾围岩背景场与局部体异常场,较之传统局部体瞬变场原理能够更全面的反映地质信息。     

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

井中瞬变电磁波勘探是一个全空间地球物理场问题.采用Gaver-Stehfest逆拉氏变换方法,正演计算了瞬变信号激励下接收线圈上的电磁场响应.分析了包含井眼泥浆、套管、水泥环和地层的轴对称多层介质模型的电磁场响应特征,考察了各层介质参数对井中瞬变电磁响应的影响.不同电导率井眼泥浆的电磁场响应衰减曲线表明,井眼泥浆电导率...     

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

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

瞬变电磁井间勘探方法研究

依据瞬变电磁全空间响应的传播特征,本文提出了一种井间勘探方法.在井内用线圈激发,在本井和邻近井接收瞬态响应波形.基于Maxwell方程分析了位移电流和传导电流所引起的两种传播特征,明确了在邻井接收到的波形中包含界面的电导率信息.详细研究了线圈在无限大均匀介质中激发的瞬变电磁响应波形以及电场在不同时刻的空间分布.瞬变电磁近源距响应波形幅度随源距的增大衰减很快,响应持续时间短;远源距对应的响应波形上升速度快,达到峰值之后缓慢下降.随着源距的增加,响应波形上升和下降的速度均减慢,峰值向后移动,能够看到明显的传播特征(与位移电流引起的电磁波形态不变差异很大).与此对应,瞬变电场的空间分布也随着波形幅度的上升迅速扩大,随着幅度的减小快速收缩.在电磁能量向外传播的过程中,遇到界面时会发生透射和反射,使得测量得到的波形形状会在无界面条件下波形基础上发生改变.界面信息隐含在这种波形形状改变中,需要通过瞬变电磁反射和透射模型进行提取,实现井间勘探.     

地-井瞬变电磁响应特征数值模拟分析

井中瞬变电磁法(Bore-hole transient electromagnetic method-BHTEM)是指接收线圈在钻井中观测瞬变场响应用以勘查深部矿产资源的勘探方法,其中以地-井(地面激发井中接收)组合方式研究最多、应用最广.本文应用时域有限差分法(FDTD),建立包含薄板导体的均匀半空间二维数学模型,采用线源为激发源,选用Mur吸收边界条件,对矩形回线源在半空间中产生的瞬变电磁场进行数值模拟,计算了低阻板状导体在均质半空间和有低阻覆盖层影响情况下的地-井瞬变电磁异常响应,并对响应的特征及规律进行研究分析,为研究地-井TEM提供参考.     

海底表面磁源瞬变响应建模及海水影响分析

根据电磁场理论，推导了磁偶源和接收点均位于海水中时层状海底模型的频域电磁场响应一般表达式，并通过此式，得到了海水为均匀半空间和有限海水深度两种情况下，垂直磁偶极装置、中心回线和重叠回线分别置于均匀半空间海底表面时的瞬变电磁响应（磁场和感应电压）表达式. 这些表达式将瞬变响应和海底的电导率等参数有机联系在一起，为海底瞬变电磁法的正演计算和反演解释提供了理论基础. 仿真计算表明，海水的存在不仅使得瞬变响应曲线形态发生变化，而且影响其对海底电导率的分辨能力.     

考虑关断时间的回线源激发TEM三维时域有限差分正演

从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1 μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线.     

基于航空大地电磁倾子特征的二维大地结构探测方法

倾子是航空大地电磁方法的主要参数,在地下异常体的探测中有着重要的作用和显著优势.本文简要概括了倾子测量的基本原理;基于有限单元法进行了二维大地结构的倾子建模仿真,分析了影响倾子响应特征的主要因素,包括空间性参数以及地下介质参数.通过对倾子响应特征影响因素的分析,提出了一种利用倾子及其梯度的异常位置和异常幅度,对地下二维异常体的电性特征、水平位置、埋深以及电阻率进行快速、定量识别的方法.最后,利用已知模型的倾子响应对异常体参数进行了试算,通过对试算结果和已知模型参数的对比,验证了该方法的可行性.     

用模式匹配算法研究层状各向异性倾斜地层中多分量感应测井响应

利用数值模式匹配技术,研究并建立层状各向异性倾斜地层中多分量感应测井响应的快速数值模拟算法.首先将位于井轴上的三个相互正交磁偶极子转化成关于极角θ的三个谐变分量的叠加,以便将三个正交磁偶极子电磁场的正演问题完全简化成三个谐变分量电磁场的轴对称定解问题,并给出电磁场各个谐变分量在井轴上满足的边界条件,保证电磁场在井轴附近仍然可解.然后利用数值模式匹配技术建立电磁场各个谐变分量的正演过程,得到层状各向异性地层中磁流源并矢Green函数的半解析表达式,给出计算层状各向异性倾斜地层中多分量感应测井响应的具体方法,最后通过数值计算结果证明该算法的有效性并考察几种不同情况下多分量感应测井响应特征.     

Transient electromagnetic surveys with unimodal transverse magnetic field: ideas and results

The theory behind transient electromagnetic surveys can be well described in terms of transverse magnetic and transverse electric modes. Soundings using transverse magnetic and transverse electric modes require different source configurations. In this study, we consider an alternating transverse magnetic field excitation by a circular electric dipole. The circular electric dipole transmitter is a horizontal analogue of the vertical electric dipole. Offshore surveys using circular electric dipole might represent an alternative to the conventional marine controlled‐source electromagnetic method at shallow sea and/or for exploring relatively small targets. Field acquisition is carried out by recording either electric or magnetic responses. Electric responses bear information on the 1D structure of a layered earth and successfully resolve high‐resistivity targets in marine surveys. Land‐based circular electric dipole soundings are affected by induced polarisation. On the contrary, magnetic responses are absent on the surface of a 1D earth, and as a result, they are very sensitive to any and even very small 3D conductivity perturbations. In addition, they are sensitive to induced polarisation or some other polarisation effects in the subsurface. At present, circular electric dipole transmitters and magnetic receivers are successfully used in on‐land mineral and petroleum exploration.     

Definitive superparamagnetic source identification through spatial,temporal, and amplitude analysis of airborne electromagnetic data

The aim of this paper is to add confidence to existing methods using decay shape analysis to detect superparamagnetic responses in airborne electromagnetic data. While expensive to acquire, vertical spatial gradient measurements of the electromagnetic signals can discriminate near‐surface superparamagnetic sources. This research investigated the use of horizontal spatial gradients and amplitude information as further indicators of superparamagnetic. High horizontal gradients were shown both theoretically and in field data to help discriminate superparamagnetic from deep mineral targets. Further, superparamagnetic responses have characteristically small amplitudes inconsistent with realistic mineral exploration targets at shallow depths.     

3D Modelling of the electromagnetic response of geophysical targets using the FDTD method1

A publicly available and maintained electromagnetic finite-difference time domain (FDTD) code has been applied to the forward modelling of the response of 1D, 2D and 3D geophysical targets to a vertical magnetic dipole excitation. The FDTD method is used to analyse target responses in the 1 MHz to 100MHz range, where either conduction or displacement currents may have the controlling role. The response of the geophysical target to the excitation is presented as changes in the magnetic field ellipticity. The results of the FDTD code compare favourably with previously published integral equation solutions of the response of 1D targets, and FDTD models calculated with different finite-difference cell sizes are compared to find the effect of model discretization on the solution. The discretization errors, calculated as absolute error in ellipticity, are presented for the different ground geometry models considered, and are, for the most part, below 10% of the integral equation solutions. Finally, the FDTD code is used to calculate the magnetic ellipticity response of a 2D survey and a 3D sounding of complicated geophysical targets. The response of these 2D and 3D targets are too complicated to be verified with integral equation solutions, but show the proper low- and high-frequency responses.     

全极化探地雷达正演模拟

全极化探地雷达不但能获取当前探地雷达使用的共极化信息,而且能使用交叉极化信息,进而更好地分析目标属性.本文在时域有限差分法的基础上,构建了全极化探地雷达的正演模拟方法,通过水平正交双方向同时接收获取全极化信息.标准目标体金属板和角反射器的正演模拟与实验室物理模型测试一致的共极化和交叉极化响应信息,证实了该正演模拟方法的...     

层状介质任意形状回线源瞬变电磁全区视电阻率的研究

利用时域垂直磁场求取全区视电阻率.任意形状回线的垂直磁场可以看成是多个水平电偶极子的垂直磁场累加而成;因此,不论回线是何形状,只要将其分解为水平电偶极子,就可利用水平电偶极子的垂直磁场求取任意形状回线的全区视电阻率.     

The geodynamo as a bistable oscillator

Abstract

Our intent is to provide a simple and quantitative understanding of the variability of the axial dipole component of the geomagnetic field on both short and long time scales. To this end we study the statistical properties of a prototype nonlinear mean field model. An azimuthal average is employed, so that (1) we address only the axisymmetric component of the field, and (2) the dynamo parameters have a random component that fluctuates on the (fast) eddy turnover time scale. Numerical solutions with a rapidly fluctuating α reproduce several features of the geomagnetic field: (1) a variable, dominantly dipolar field with additional fine structure due to excited overtones, and sudden reversals during which the field becomes almost quadrupolar, (2) aborted reversals and excursions, (3) intervals between reversals having a Poisson distribution. These properties are robust, and appear regardless of the type of nonlinearity and the model parameters. A technique is presented for analysing the statistical properties of dynamo models of this type. The Fokker-Planck equation for the amplitude a of the fundamental dipole mode shows that a behaves as the position of a heavily damped particle in a bistable potential ∝(1 ? a ²)², subject to random forcing. The dipole amplitude oscillates near the bottom of one well and makes occasional jumps to the other. These reversals are induced solely by the overtones. Theoretical expressions are derived for the statistical distribution of the dipole amplitude, the variance of the dipole amplitude between reversals, and the mean reversal rate. The model explains why the reversal rate increases with increasing secular variation, as observed. Moreover, the present reversal rate of the geodynamo, once per (2?3) × 10⁵ year, is shown to imply a secular variation of the axial dipole moment of ～ 15% (about the current value). The theoretical dipole amplitude distribution agrees well with the Sint-800 data.     

THE CROSS-HOLE MAGNETOMETRIC RESISTIVITY (MMR) RESPONSE OF A DISC CONDUCTOR*

The cross-hole variant of the magnetometric resistivity (MMR) method requires two bore holes in the vicinity of a conductive target. In the first, two fixed current electrodes are located, one above the other. They are linked to a low frequency current source by cables, the whole system forming a vertical current bipole. In the second, a sensitive coil measures the axial magnetic field as a function of depth. For a uniform earth, if both holes are vertical, the measured component vanishes by symmetry. However, the presence of a local conductor channels the current and causes an anomalous magnetic component which is interpreted to indicate the position, shape and relative conductance of the target. Mineral deposits are often lamellar in form. The conductive disc is the simplest bounded lamella for which MMR responses may be computed. It is excited by a single current source on its axis. The second source and the surface of the earth are assumed to be far away, a valid assumption for down-hole measurements. The numerical method introduces a new integral equation describing the interaction of current dipoles located in the plane of the disc. The equation is solved analytically for a disc of infinite radius, a layer, and the result is compared with a corresponding known boundary value solution. The computed radial current in the disc and the magnetic field generated by it are described in terms of a current channelling number. The magnitude of the computed field is of the order of one nanoTesla for a typical mining problem.     

阵列偶极源辐射场在导电地层中传播的方向特性

利用电偶极子辐射场的矢势导出了阵列偶极子在导电媒质中辐射场的一般形式,所给出的公式不仅适用于近区、感应区,而且还适用于远区.对阵列偶极子辐射场的相干特性进行了详细的讨论,由所得结果可知,对层状地层通过选择适当的方位可使接收到的电磁信号达到极大,极大值的大小与偶极源的个数成正比.本文的研究对解决地层中电磁信号衰减严重问题提供了一种可能的解决办法.     

BZ参数瞬变电磁法的优势及在沿海滩涂区的验证试验

随着我国一系列发展战略的实施,在沿海滩涂区等低阻环境下开展电磁法深部探测工作的需求越发迫切.传统瞬变电磁法(TEM)一般观测感应电动势(dBz/dt)参数,其响应函数复杂,数据解译难度大,且衰减极快,勘探深度较小,不能满足深部探测的需求.而瞬变电磁法的感应磁场(Bz)参数具有响应函数单调、衰减速度慢等优点,深部探测应用前景广阔.本文首先从理论上讨论了TEM的Bz参数相比dBz/dt参数在响应曲线形态、探测深度等方面的优越性,然后在某沿海滩涂区开展了野外验证工作.理论分析计算和野外试验结果表明,TEM的Bz参数相较于dBz/dt参数有诸多优势,基于Bz参数的瞬变电磁法在滩涂区的极低阻条件下仍具有较强的分层能力、较大的探测深度和较强的适用性.