矿井瞬变电磁法全空间视电阻率解释方法研究

从全空间中心回线源响应公式出发,推导了全空间晚期和全区视电阻率公式;通过与半空间公式比较和对实测资料的分析,研究了全空间视电阻率公式的应用效果.结果显示,全空间视电阻率为半空间视电阻率的(5/2)^2/3倍,其晚期视电阻率与真值的相对误差小于0.56%,而相应的半空间视电阻率与真值误差均大于42%.与常规解释方法相比,采用全空间全区视电阻率换算方法能更好地圈定低阻异常范围,且电阻率差异大.因此,该公式更适用于矿井瞬变电磁法资料处理与解释,并具有更好的效果.     

Feasibility of central loop TEM method for prospecting multilayer water-filled goaf

With deep mining of coal mines, prospecting multilayer water-filled goaf has become a new content that results from geophysical exploration in coalfields. The central loop transient electromagnetic (TEM) method is favorable for prospecting conductive layers because of the coupling relationship between its field structure and formation. However, the shielding effect of conductive overburden would not only require a longer observation time when prospecting the same depth but also weaken the anomalous response of underlying layers. Through direct time domain numerical simulation and horizontal layered earth forward modeling, this paper estimates the length of observation time required to prospect the target, and the distinguishable criterion of multilayer water-filled goaf is presented with observation error according to the effect of noise on observation data. The observed emf curves from Dazigou Coal Mine, Shanxi Province can distinguish multilayer water-filled goaf. In quantitative inversion interpretation of observed curves, using electric logging data as initial parameters restrains the equivalence caused by coal formation thin layers. The deduced three-layer and two-layer water-filled goafs are confirmed by the drilling hole. The result suggests that when observation time is long enough and with the anomalous situation of underlying layers being greater than the observation error, the use of the central loop TEM method to prospect a multilayer water-filled goaf is feasible.     

露天煤矿地下水电磁探测与类型判别

为准确确定露天煤矿地下水分布与类型,根据不同赋存条件电磁差异,基于磁场强度与电阻率参数进行异常划分与类型判别。以新疆某露天煤矿为例,采用高精度磁法通过磁异常强度确定烧变岩边界,并采用高密度电法和瞬变电磁法通过电阻率圈定低阻异常区,结合异常响应差异分别推断地下水分布范围与赋存类型。地面钻探验证探测成果的准确性。结果表明:露天煤矿烧变岩水与砂岩裂隙水电磁响应差异明显,基于磁场强度与电阻率的双参数综合电磁法能准确圈定其分布范围,并判别其赋存类型。     

瞬变电磁法中心回线装置资料解释方法的改进

中心回线装置是瞬变电磁勘探中最常用的装置之一.中心回线装置的视电阻率一般从回线中心点场的公式出发导出.在现场实际施工中,为提高工作效率,把发射回线中部三分之一的区域作为观测范围.通过对场分布特征的分析研究表明,与中心点相比,观测区边缘处的感应电动势数值偏离达15%~25%,这与广泛存在于华北型煤系中赋存深度为400~1000 m的陷落柱、导水小断层等引起的异常相比,已经不可忽略.大定源回线公式可以准确地表示任意场点的感应电动势,由此导出视电阻率无边缘效应影响;观测按照中心回线方式并保持在近区进行,可使感应电动势和视电阻率之间的转换简单而直接.理论和实际应用结果表明,这些措施进一步提高了瞬变电磁探测中心回线方法的探测精度.在有上覆低阻屏蔽层的情况下,对埋深500~700 m的陷落柱给出了明显的异常反映,并被井下掘进所验证.     

电性源瞬变电磁地空逆合成孔径成像

电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性.     

多辐射场源地空瞬变电磁法多分量全域视电阻率定义

本文针对多辐射场源地空瞬变电磁法理论,建立了相应的多分量全域视电阻率定义方法.分析了利用磁场强度进行多辐射场源地空系统全域视电阻率定义的优点,针对磁场强度的各个分量,提出了各自的全域视电阻率算法,实现了多分量、全时域、全空域视电阻率计算,并分析了偏移距对全域视电阻率的影响.通过调整源的相对位置及电流方向等参数,多辐射源瞬变电磁地空系统不仅可以加强不同分量信号强度,削弱随机干扰,还可更好地分辨地下异常体的位置.通过对所设计模型的处理,证实了多辐射场源地空系统多分量全域视电阻率算法的有效性,也验证了多辐射场源地空系统的优势.     

Continuous TDEM for monitoring shale hydraulic fracturing

Monitoring and delineating the spatial distribution of shale fracturing is fundamentally important to shale gas production. Standard monitoring methods, such as time-lapse seismic, cross-well seismic and micro-seismic methods, are expensive, timeconsuming, and do not show the changes in the formation with time. The resistivities of hydraulic fracturing fluid and reservoir rocks were measured. The results suggest that the injection fluid and consequently the injected reservoir are characterized by very low resistivity and high chargeability. This allows using of the controlled-source electromagnetic method (CSEM) to monitor shale gas hydraulic fracturing. Based on the geoelectrical model which was proposed according to the well-log and seismic data in the test area the change rule of the reacted electrical field was studied to account for the change of shale resistivity, and then the normalized residual resistivity method for time lapse processing was given. The time-domain electromagnetic method (TDEM) was used to continuously monitor the shale gas fracturing at the Fulin shale gas field in southern China. A high-power transmitter and multi-channel transient electromagnetic receiver array were adopted. 9 h time series of Ex component of 224 sites which were laid out on the surface and over three fracturing stages of a horizontal well at 2800 m depth was recorded. After data processing and calculation of the normalized resistivity residuals, the changes in the Ex signal were determined and a dynamic 3D image of the change in resistivity was constructed. This allows modeling the spatial distribution of the fracturing fluid. The model results suggest that TDEM is promising for monitoring hydraulic fracturing of shale.     

The geoelectrical properties of the earth around the Karymshina multidisciplinary observation site in southern Kamchatka

This paper is concerned with the Earth’s electromagnetic field, a deep geoelectrical section, and the dynamics of earth conductivity. This analysis is based on MTS curves along directions that to a first approximation go along and across southern Kamchatka. It is shown that the longitudinal and transverse MTS curves are subject to the influence of local and regional geoelectrical inhomogeneities. The coast effect was studied by 3-D numerical modeling. The patterns we have found were used for interpretation of generalized MTS curves. The resulting parameters of the geoelectrical section were refined by reducing the longitudinal curve to the standard curve of apparent resistivity. The results from this interpretation are to be refined as more geoelectrical information is forthcoming. The dynamics of lithospheric conductivity were studied from data acquired during the 2005–2008 monitoring of magnetotelluric impedance in the range of periods between several hundred and several thousand seconds. The basis for the analysis was assumed to be the transverse impedance and its phase. The latter was found to exhibit anomalous bay fluctuations that might have been related to earthquakes. A possible origin of these anomalies is discussed.     

Sensitivity of transversal apparent resistivity to the changes in electrical resistivity of the elements of a 2D geoelectrical structure

The sensitivity of transversal apparent resistivity to the changes in the electrical resistivity of elements of a two-dimensional (2D) geoelectrical structure is studied by numerical modeling. This sensitivity is found to have a series of specific features, due to which the monitoring of transversal apparent resistivity can serve as a helpful addition to the monitoring of longitudinal apparent resistivity to trace the dynamics of subsurface and shallow crustal elements of the medium. It is shown that the method previously suggested for the inversion of relative changes in apparent resistivity into relative changes in electrical resistivities of the elements of the geoelectrical structure is applicable to the transversal electrical resistivity.     

二维水平电偶极变频测深阻抗视电阻率的有限元正演计算

在二维地电断面情况下,电偶源变频测深阻抗视电阻率的有限元计算,存在着两个困难:一个是对源的处理;另一个是边界条件的确定。本文提出了解决这两个问题的方法,构成了含源波动方程的边值问题,用有限元法求出了相应的电磁场分量,从而计算出二维地电断面的阻抗视电阻率。对各种地电断面计算的结果表明,计算方法是正确的,可对实际地电断面进行有限元模拟,求得其电磁响应。     

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

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

时间域瞬变电磁法全区视电阻率的平移算法

深入分析了均匀半空间中瞬变电磁场的解析表达式,发现瞬变响应曲线随地下电导率、发射回线边长与观测时间具有平移伸缩特性,并据此提出了一种新的直接计算全区视电阻率的方法.该法与现有计算全区视电阻率的方法的不同点在于：无需迭代和求解非线性方程,速度快,精度高,算法稳定,且只有惟一解.     

核磁共振与瞬变电磁三维联合解释方法

传统核磁共振地下含水量解释多采用基于均匀半空间或层状导电模型的一维反演,分层给出地下含水信息.然而,这些方法忽略了地下复杂电阻率分布信息对结果的影响,也不能很好地反映局部三维含水构造.本文从三维电介质中核磁共振响应的正演理论出发,提出首先利用瞬变电磁数据进行基于等效导电平面法的快速电阻率成像,然后将成像结果作为核磁共振三维反演的电性模型,进行联合解释.激发磁场的分布采用有限元法直接求解,通过引入伪δ源实现电流源的加载,并强加散度条件排除了三维磁场模拟中"弱解"的影响.针对核磁共振灵敏度矩阵的病态性和数据中存在的干扰信号,提出考虑罚项的非线性拟合目标函数,利用线性化方法进行核磁共振反演.模型数据表明该方法能较准确反映地下三维含水构造,实测算例进一步证明了方法的有效性.本研究将促使核磁共振方法在岩溶、裂隙水、孤立水体等复杂水文地质条件及隧道、矿井灾害水源探测等方面得到有效应用.     

基于神经网络的视电阻率快速算法

本文从瞬变电磁均匀半空间二次磁场响应公式出发,提出了一种基于神经网络的视电阻率快速计算方法.以中心回线为例,根据瞬变响应公式的特点,简化网络结构,选用三层BP神经网络和误差训练算法,用均匀半空间样本数据进行训练,确定了收敛快、误差小的一步正割法和隐含单元数,得到基于不同采样时窗的一组网络参数.用本文方法与二分法、牛顿迭代法做模型计算比较,及最后的实验计算,说明算法的快速,准确.本文方法不依赖初始模型,避开了复杂的电磁场数值计算,实现了视电阻率的快速计算,对瞬变电磁法资料的快速解释有一定的参考价值.     

瞬变电磁法的探测深度问题

用解析分析、时域有限差分、时-频分析的方法,以地面中心回线装置和阶跃脉冲激励源为例,分析讨论了瞬变电磁测深法的勘探深度问题,以便为野外勘探工作设计提供依据,达到预期的探测目的.解析计算证实了瞬变场在地下以有限速度传播,数值模拟表示出了准静态条件下瞬变场的反射.研究结果表明,由于时间域电磁场遵循因果律,瞬变电磁法的探测深度主要由观测时间决定. 瞬变电磁场的初始传播速度与大地电阻率无关,继后在大地色散作用下,阶跃脉冲前沿逐渐变得平缓,各频率分量的传播速度与电阻率有关,在低阻地层中探测同样的深度需要较长的观测时间. 最大探测深度是在给定时间内电磁波往返地下某一深度的单程距离,最小探测深度受仪器性能的限制,但是埋藏较浅的异常体也有可能在晚时段被观测到.从时-频密度谱中可得到瞬变电磁场信号时间与频率的关系.     

地震地球物理观测台网地电数据集

我国地震地球物理观测台网中地电站网由92个地电阻率站点和111个地电场站点组成,地电数据集由地电阻率、地电场观测数据及相关产品构成,含我国地电站网在运行131个观测台站的原始观测数据、预处理数据和产品数据,并附带观测日志、基础信息等辅助信息,可为地震预测预报、地球动力学、电磁环境监测等多方面的科学研究提供强有力的数据支撑。     

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.     

CSAMT全区电阻率法数值模拟及应用探讨

常规的可控源电磁法理论在计算视电阻率公式上,多半采用其电磁场的渐近特征,难以直接反映全区视电阻率的值,及直观地显现地下介质的地质构造.文中采用水平偶极子激发的电磁场,提出了电场的全区精确表达式,直接计算出大地电阻率.利用汉克尔数值滤波算法和逆样条插值算法对水平层状电磁场进行正演计算,并与计算的卡尼亚视电阻率的对比和野外试验结果表明：该方法的结果在远区等价卡尼亚电阻率,在近区和过渡带则明显地改善了卡尼亚电阻率的非波场区场畸变,从而能更好地接近基底的真电阻率,更形象地反映了地下介质的垂向电性变化.     

Joint inversion of long‐offset and central‐loop transient electromagnetic data: Application to a mud volcano exploration in Perekishkul,Azerbaijan

Mud volcanism is commonly observed in Azerbaijan and the surrounding South Caspian Basin. This natural phenomenon is very similar to magmatic volcanoes but differs in one considerable aspect: Magmatic volcanoes are generally the result of ascending molten rock within the Earth's crust, whereas mud volcanoes are characterised by expelling mixtures of water, mud, and gas. The majority of mud volcanoes have been observed on ocean floors or in deep sedimentary basins, such as those found in Azerbaijan. Furthermore, their occurrences in Azerbaijan are generally closely associated with hydrocarbon reservoirs and are therefore of immense economic and geological interest. The broadside long‐offset transient electromagnetic method and the central‐loop transient electromagnetic method were applied to study the inner structure of such mud volcanoes and to determine the depth of a resistive geological formation that is predicted to contain the majority of the hydrocarbon reservoirs in the survey area. One‐dimensional joint inversion of central‐loop and long‐offset transient electromagnetic data was performed using the inversion schemes of Occam and Marquardt. By using the joint inversion models, a subsurface resistivity structure ranging from the surface to a depth of approximately 7 km was determined. Along a profile running perpendicular to the assumed strike direction, lateral resistivity variations could only be determined in the shallow depth range using the transient electromagnetic data. An attempt to resolve further two‐dimensional/three‐dimensional resistivity structures, representing possible mud migration paths at large depths using the long‐offset transient electromagnetic data, failed. Moreover, the joint inversion models led to ambiguous results regarding the depth and resistivity of the hydrocarbon target formation due to poor resolution at great depths (>5 km). Thus, 1D/2D modelling studies were subsequently performed to investigate the influence of the resistive terminating half‐space on the measured long‐offset transient electromagnetic data. The 1D joint inversion models were utilised as starting models for both the 1D and 2D modelling studies. The results tend to show that a resistive terminating half‐space, implying the presence of the target formation, is the favourable geological setting. Furthermore, the 2D modelling study aimed to fit all measured long‐offset transient electromagnetic Ex transients along the profile simultaneously. Consequently, 3125 2D forward calculations were necessary to determine the best‐fit resistivity model. The results are consistent with the 1D inversion, indicating that the data are best described by a resistive terminating half‐space, although the resistivity and depth cannot be determined clearly.     

Transient process and optimal design of receiver coil for small‐loop transient electromagnetics

Distributed parameters of the receiver coils greatly affect transient electromagnetic signals over short time periods, causing a delay in the signal's effective sampling time and the loss of shallow exploration information. This paper investigates the influence of transient process on apparent resistivity calculation and analyses the relations between the error of apparent resistivity and receiver coil design. We find that, under the same effective area, different radii of the receiver coils lead to different levels of impact on the estimation of the apparent resistivity. An optimization model is proposed to determine the optimal receiver coil size that gives rise to the smallest estimation error of the apparent resistivity. The relationship between the optimal radius and the effective areas is developed, which serves as a guideline for the optimal receiver coil design. The results may provide a useful means for improving the accuracy of the small loop transient electromagnetic instrumentation for shallow‐depth mapping.