ENHANCEMENT OF ELECTROMAGNETIC ANOMALIES BY A CONDUCTING OVERBURDEN II *

The influence of a partially conducting overburden/host rock on the electromagnetic response of a horizontal, tabular conducting ore body, investigated with the aid of quantitative scale model experiments, was discussed in an earlier paper (Gaur, Verma and Gupta 1971), which will be referred to as I. This paper presents the results of more comprehensive experiments subsequently carried out to study the combined effect of various geological parameters, namely: the dip and depth of burial of the ore body, its electrical conductivity and that of the zone surrounding it. These results obtained for four different transmitter-receiver configurations confirm the general enhancement of response in varying degrees, brought about by a conducting overburden. However, the transformation of the shapes of the anomaly profiles with increasing overburden conductivity is observed to depend on the dip of the ore body, being more drastic for gently dipping ones. Variations in the inphase and quadrature components as well as in the phase of the anomaly have been studied for varying depths of burial of the ore body and for a number of values of the solution conductivity. Anomaly index diagrams have been constructed with a view to predicting possible values of the geological parameters from a knowledge of the anomaly components. It is felt that the notable overburden effects are caused by a drastic redistribution and concentration of currents, mostly in a narrow loop at the top of the ore body, brought about by its galvanic contact with an extended medium of relatively poorer conductivity.     

VLF ANOMALIES FROM A PERFECTLY CONDUCTING HALF PLANE BELOW AN OVERBURDEN*

A theoretical solution to the electromagnetic problem of a perfectly conducting half plane below a conducting overburden has been obtained. The VLF anomalies have been computed for different overburden conductivity and thickness and also for different dip angles of the half plane. In the computations the contribution to the secondary magnetic field from the electric Hertz potential has been neglected. The anomaly curves which are displayed as EM 16 readings, show a fairly complicated behaviour. This is mainly due to the phase shift and attenuation of the field caused by the conductivity of the overburden and the host rock. From the anomaly curves it is possible to define the apparent depth to the top of the conductor as the distance between the peak value and the cross-over of the real component. The apparent depth is usually larger than the actual depth, but it is possible to determine the actual depth to the conductor from the relation between the peak-to-peak anomaly and the apparent depth. When the peak-to-peak anomaly is fairly large, it is also possible to make estimates of the dip angle. However, a complete set of master curves will be a necessary tool for interpretation of VLF data when there is need to obtain more accurate estimates of the half plane parameters. In a specific case the theoretical calculations are shown to be in good agreement with measured data.     

ENHANCEMENT OF ELECTROMAGNETIC ANOMALIES BY A CONDUCTING OVERBURDEN*

The effect of a conducting overburden on the electromagnetic response of sulphide ore-bodies has been studied with the help of quantitative electromagnetic model experiments. These experiments were conducted at a fixed, crystal controlled, frequency of 100 kHz using a number of transmitter-receiver configurations, though the results discussed here mostly pertain to a horizontal coplanar system. An analysis of the anomaly profiles—after accounting for the regional anomaly—indicates a general enhancement of the response. This is attributed (a) to the concentration of the current lines in the host-rock effected by the embedded target and (b) to the change in the phase and space orientation of the field vectors, brought about by the conducting overburden. The enhancement is most pronounced in the case of non-symmetrical bodies such as veins and ribbon-like structures represented by sheet models, but is also observed to a lesser degree for isometric structures. These results are expected to pave the way for a more accurate interpretation of the induction prospecting data.     

INDUCTION ANOMALY DUE TO AN ELONGATED COVERED ORE-ZONE EXCITED BY A LONG CURRENT CARRYING CABLE*

A detailed analysis of the electromagnetic anomaly due to an elongated mineralized zone of low resistivity in presence of a fixed-transmitter (a long cable carrying an, alternating current Ie^Iωt) is presented where (i) the host rock is not highly resistive and might contain some disseminated mineralization, (ii) the ore-body has an inhomogeneous conductivity, (iii) there is a contrast in magnetic properties of the ore-body and the host rocks, and (iv) source is close to the conducting system i.e., anomalous zone is lying shallow. The numerical computation of the generalized analytical expressions have been made. The selective screening behaviour of the cover is delineated and suitable frequency ranges for maximum detectability of a covered conductor have been obtained. A paradoxical decrease of the in-phase component of the response function with the increase of the core conductivity has been found for large conductivities of the cover. For moderate values of response parameter the inhomogeneity in the conductivity is found to significantly affect the induced field. The present analysis of the various features of the secondary field will aid to the existing interpretation sensitivity of the induction prospecting data for porphyry conducting ore deposits with zonal wall-rock alteration and sulphide distributions.     

深海热液硫化物矿体3D瞬变电磁正演

深海热液硫化物矿体瞬变电磁的正演是考虑深海环境的全空间条件下三维体的涡流电磁响应.采用全空间矢量有限元法模拟计算深海热液硫化物矿的三维瞬变电磁响应,对硫化物矿体采用矩形单元模型剖分,应用Galerkin法推导有限元方程,先计算频率域响应,再通过Fourier反变换将其转换至时间域,得出深海热液硫化物矿矿体的瞬变电磁响应.并用双半空间模型的解析解检验了全空间矢量有限元法模拟计算算法和程序的正确性,最后按照等比例缩小电磁物理实验原则,比对数值计算和物理实验结果论证了全空间3D模型数值的正确性.结果表明:对于海水、矿体以及围岩复杂电磁边界,应用全空间矢量有限元法模拟计算深海热液硫化物矿瞬变电磁响应异常与物理模拟结果一致,而且计算方法简单精确,异常幅值明显,边界清晰.     

基于二次场电导率分块连续变化的三维可控源电磁有限元数值模拟

从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态.     

COMPUTATION OF VLF RESPONSE OVER HALF-PLANE AND WEDGE MODELS*

A theoretical solution is presented to the problem where a VLF anomaly is generated by a conducting half-plane or a perfectly conducting wedge below a stratified overburden. The solution is obtained by the use of a scattering matrix for plane-wave eigenfunctions. VLF anomalies have been computed for different values of the conductance and dip of the half-plane. The phase of the VLF anomaly due to a conducting half-plane depends on the conductance and the distance to the half-plane. Close to the half-plane the tilt angle and ellipticity are of opposite sign for a perfect conductor, but the ellipticity will change sign for a poor conductor. The VLF anomaly for a perfectly conducting wedge is essentially determined by the position of the upper surface of the wedge, i.e. the anomaly will closely resemble the anomaly of a perfectly conducting half-plane in the same position as the upper surface of the wedge.     

单一层面高频电偶极子场的近似计算

本文对一个倾斜平面界面的高频电磁波场近似计算问题做了解答。这是钻井高频电磁波勘探方法中的基本问题。文中写出了问题的严格解,并用鞍点法对φ=0的条件,得出其近似计算式和数值结果.这个结果具有便于计算和明确的物理意义。根据计算结果,指出了特别在高导介质中,界面附近的干涉异常,必须予以充分重视。     

Analogue model measurements for a horizontal magnetic dipole embedded within a conducting medium

The validity of an analogue model method employing a horizontal magnetic dipole source situated within a conducting layer for the cases of a poorly conducting model earth and a highly conducting model ocean is studied by comparing model magnetic field measurements with theoretical calculations. The model is then used to study one example for each case; the response of a conducting cylinder simulating an ore body embedded in the earth, and the response of a conducting wedge simulating a shelving ocean.     

深海6000 m拖曳式瞬变电磁系统及其应用

针对我国西南印度洋合同区热液硫化物矿快速、有效以及便捷地质勘探装备的要求,研发了由甲板控制系统、万米光电复合缆、仪器舱拖体和天线拖体组成的深海6000 m拖曳式瞬变电磁系统.为了便于拖体布放和快速发现异常,选择重叠回线收发装置类型;采用理论数值计算确定了在拖曳高度不大于50 m的前提下,观测1～100ms窗范围内的二次场响应,可以发现近海底深海热液硫化物矿堆;另外,通过不同拖延深度海上试验,研究了拖曳深度对瞬变电磁的影响规律及仪器性能.大洋第30航次第二航段在西南印度洋脊热液区,应用该系统发现了明显的瞬变电磁异常.印证了系统的有效性.     

水平n层导电介质中时谐垂直电偶极子的电磁场

对于包含n层线性、均匀、各向同性的水平成层导电介质,本文应用汉克尔变换及其逆变换推导求得了位于顶层导电媒质中的时谐垂直电偶极子在该层产生的电磁场的表达通式．通过这些表达式能够很方便地求得包含任意多层导体半空间中的电磁场．文中最后计算了导电半体空间和二层导电介质中时谐垂直电偶极子产生的电磁场,验正了本文所得结果的正确性。     

导电率随时间变化的导电球的电磁响应问题

本文采用非齐次麦克斯韦方程组,给出了导电率随时间变化时均匀的和非均匀的导电球,在均匀的和非均匀的外激励场的作用下的电磁响应问题的解析解。并简单讨论了这类电磁响应问题在地球物理学中的可能应用。     

多通道瞬变电磁法2D有限元模拟

多通道瞬变电磁法的反演工作是在大地脉冲响应的基础上进行的,因此本文用2D有限元正演来求取大地脉冲响应,实现多道瞬变电磁法的正演模拟研究.研究工作首先在有源大地电磁法2D有限元模拟的基础上进行,原有2D有限元法的频带相对于多通道瞬变电磁法的频带是窄频的,为此采用校正法将其扩展到所需的宽度,实现了频率域宽频带有源电磁勘探方法的正演模拟;然后通过频时变换变换到时间域,得到瞬变电磁法的阶跃响应;最后通过求取阶跃响应的时间导数,得到大地脉冲响应.针对多通道瞬变电磁方法的装置特点,用有一定埋深的顺层成矿模型分析了大地脉冲响应特征,发现对有一定埋深的顺层产出的矿脉模型来说,频率域电磁场结果以及时间域大地脉冲响应结果均需移除背景场才能突出异常体的存在,然而,在实际工作中,未知的背景场限制了其实用性.为此,本文利用多通道瞬变电磁方法数据量大的特点,提出用不同偏移距的等时曲线与等时断面来展示脉冲响应模拟结果,事实表明这两种展示方式可以清晰地分辨出矿体引起的异常,且能准确定位矿体的中心在地面的投影位置,从而说明多通道瞬变电磁方法相对于频率域有源电磁方法,对有一定埋深的顺层矿床有更高的分辨能力.     

Analysis of a guided wave along a conducting structure in a borehole

A guided wave along a borehole is often observed in borehole radar measurements. These guided waves deform the antenna pattern and can cause artefacts in radar measurements. A water-filled borehole or a conducting logging cable can function as a waveguide for electromagnetic waves under some conditions. We describe the theoretical characteristics of such a guided wave in a borehole and compare them with our experiments. The measured signal discussed was obtained with a directional borehole radar. This radar uses a cylindrical conformal array antenna as receiver and is a model of a conducting structure in a borehole. The induced field around the borehole was compared with the theory. The most fundamental symmetrical and asymmetrical modes were TM₀₁ and HE₁₁, and they were identified in the measured signals using time–frequency distribution analysis and by observation of the azimuthal field distribution of the magnetic field.     

地球物理探测卫星数据分析处理技术与地震预测应用研究项目及研究进展

地震监测预警和预测预报是当前地球科学及相关学科所面临的最艰巨的问题之一,是关系到人类社会安全与国计民生的亟待攻克的科学难题.为进一步提高地震预测科学研究水平,推进地震监测预测能力建设,我国于世纪之交提出了建立地震立体观测体系的战略发展思路,并希望突破三维地球物理场获取能力瓶颈,发展地球多圈层耦合模型,通过卫星观测获取全...     

电磁监测台站观测的舒曼谐振背景变化

舒曼谐振是由闪电激发的电离层与地面之间的全球电磁震荡,在地球与电离层空腔中传播,由于共振作用,其电磁波能量明显高于其他频率电磁波能量.舒曼谐振具有稳定的频域参数和频谱结构,近几年的研究发现,地震发生前会使这些固有参数发生扰动,因此利用舒曼谐振异常来进行地震短临监测可能是一种行之有效的手段.要判断舒曼谐振地震电磁异常,了解舒曼谐振电磁场背景变化规律是前提.本文利用极低频电磁台站观测的天然电磁场数据,处理和分析了各台站观测的舒曼谐振频率电磁场功率谱随时间的变化,得到了舒曼谐振随季节和纬度的背景变化规律.并提出将舒曼谐振背景变化规律应用到地震短临监测预报中,能够为今后辨别地震监测中的舒曼谐振异常提供物理依据.     

INLINE AND BROADSIDE EM DIPOLE DIPOLE PROFILING OVER A THIN VERTICAL INFINITELY CONDUCTING VEIN *

Master diagrams for electromagnetic responses of nine dipole dipole systems are presented for various depths of a vertical infinitely conducting vein. The YY system gives the minimum anomaly for both the inline and broadside arrays. Among other inline systems, it is difficult to decide the clearcut superiority of one system over the other, whereas the XX system in broadside array gives maximum anomaly.     

A study on the discrete image method for calculation of transient electromagnetic fields in geological media

We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized: the real number of the electromagnetic field solution based on the Gaver–Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the “ground excitation–stratum measurement” method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drillhole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.     

A numerical study of the topographic effect on electromagnetic fields in a three-dimensional conductivity model

Summary The topographic effect on electromagnetic fields in a three-dimensional conductivity model is investigated using a finite-difference numerical method. The model is that of a conducting block buried in the host Earth beneath a conducting hill. The altitude effect as well as the effect of conductivity of the hill has been studied by comparing the results for the field components at the surface of the hill with those at the surface of a flat half-space. The results indicate that both altitude and the hill conductivity influence the behaviour of the electromagnetic fields. For certain traverses beyond the edge of the block, it is seen that the variations in some field components are mainly due to the topographic effect and that the contribution from the block itself is negligible.     

An introduction to electromagnetic induction in the ocean

The governing equations for the induction of electromagnetic fields in the ocean by ionospheric and oceanic sources are presented. A uniformly conducting layered model and a nonuniformly conducting thin-sheet model are discussed with reference to the interpretation of fields observed in the ocean. A procedure for the separation of the electric field continuum into parts of ionospheric and oceanic origin is presented.