电(磁)法技术在地下管线探测中的应用

电(磁)法主要有充电法、高密度电法、甚低频法、电磁感应法和电磁波法(探地雷达)等。在地下管线探测中,电(磁)法是最主要的探测方法。本文介绍了电(磁)法探测地下管线的原理和特点,并列举了各种电(磁)法技术在地下管线探测中的应用。笔者对电(磁)法技术在地下管线探测中的应用前景进行了展望,指出电(磁)法技术在地下管线探测中拥有巨大的应用空间。

基于数字功放的任意波形电(磁)法信号发送技术研究

针对桥式逆变发送机系统只能发送跳变的矩形波信号且接收回路易受电磁耦合尖峰干扰以及谐波污染等问题,设计了基于数字功放的任意波形电(磁)法信号发送方案,通过对低通滤波后的矩形波放大输出,显著降低发送机输出矩形波边沿变化率,压制人工源信号中高次谐波电流在接收机输入回路中引起的电磁耦合干扰,或直接对正弦波或组合正弦波信号放大输出以压制电磁耦合干扰、消除谐波污染.本文分析了电磁耦合产生机理,模拟计算了发送低通滤波处理后的矩形波对接收回路耦合尖峰的压制效果,搭建了基于数字功放的人工源电(磁)法信号发送装置,对包括低通滤波后的组合矩形波信号在内的任意波形信号进行放大输出测试,并以双频信号为例进行了感应耦合压制实验,取得了良好的实验效果,具有较重要的理论和实际应用价值.

杭州地区地下空间调查地球物理方法技术研究

随着中国经济的快速发展和城市化进程的加快,有限的土地资源和城市发展之间的矛盾越来越突出,城市地下空间的安全、合理利用和地质环境保护具有重要的战略意义。为了解决G20、亚运会以及数字经济为杭州市快速发展带来的人口快速增长与土地资源有限的瓶颈问题,针对以杭州为代表的南方丘陵地区地下空间精细探测需求,开展了弹性波法、电法与电磁法等多种地球物理方法的可行性研究。结果表明:不同勘探方法在探测深度、分辨率以及勘探效率上具有明显的差异性,需要根据地下地质状况以及地表条件选择合适的勘探方法进行探测,这对类似丘陵地区城市地下空间开发及利用具有指导意义和参考价值。

Electrical Characteristics of Tangjiawan Landslide in Lixian, Sichuan

A wide range of terrain features and landforms,which are exemplified by intricate geological formations and diverse rock compositions,are found in the western mountainous regions of China.These areas frequently encounter geological disasters.As one of the natural disasters,landslides lead to considerable loss of human life and property.Considering mitigation of the losses caused by landslide disasters,a necessary measure for disaster prevention and mitigation involves conducting detailed investigations and monitoring of landslides,which is also the cornerstone of landslide warning.This study compares and analyzes the feasibility of the magnetotelluric detection method for landslides using the results of engineering geological surveys and landslide monitoring.The study aims to address the scientific problem of the validity of using magnetotelluric methods to detect landslide development processes.The Tangjiawan landslide signal on the left side of the K94+000~K94+145 section of the Wenma Expressway is analyzed by employing engineering geological survey,magnetotelluric detection,landslide monitoring,landslide analysis,and other methods.Analysis results provide the static electrical characteristics of lithology,structure,and groundwater,as well as the dynamic electrical characteristics of landslide development.This study focuses on analyzing the relationship between the methods of magnetotelluric detection and engineering geological surveys and the results of landslide monitoring.The workflow and methods for data collection,processing,inversion,interpretation,and analysis using the magnetotelluric method to detect the dynamic development process of landslides are presented in the conclusion.Preliminary conclusions indicate a strong correlation between the dynamic changes in magnetotelluric wave impedance with the surface displacement of landslides and the dynamic changes in groundwater.The use of the magnetotelluric method for landslide detection and monitoring is a feasible example.The research results can offer certain technical references for the detection and monitoring of landslides using magnetotelluric methods and also provide references and guidance for the selection of diversified landslide monitoring methods in the future.

Electrical structure identification of deep shale gas reservoir in complex structural area using wide field electromagnetic method

To fully exploit the technical advantages of the large-depth and high-precision artificial source electromagnetic method in the complex structure area of southern Sichuan and compensate for the shortcomings of the conventional electromagnetic method in exploration depth,precision,and accuracy,the large-depth and high-precision wide field electromagnetic method is applied to the complex structure test area of the Luochang syncline and Yuhe nose anticline in the southern Sichuan.The advantages of the wide field electromagnetic method in detecting deep,low-resistivity thin layers are demonstrated.First,on the basis of the analysis of physical property data,a geological–geoelectric model is established in the test area,and the wide field electromagnetic method is numerically simulated to analyze and evaluate the response characteristics of deep thin shale gas layers on wide field electromagnetic curves.Second,a wide field electromagnetic test is conducted in the complex structure area of southern Sichuan.After data processing and inversion imaging,apparent resistivity logging data are used for calibration to develop an apparent resistivity interpretation model suitable for the test area.On the basis of the results,the characteristics of the electrical structure change in the shallow longitudinal formation of 6 km are implemented,and the transverse electrical distribution characteristics of the deep shale gas layer are delineated.In the prediction area near the well,the subsequent data verification shows that the apparent resistivity obtained using the inversion of the wide field electromagnetic method is consistent with the trend of apparent resistivity revealed by logging,which proves that this method can effectively identify the weak response characteristics of deep shale gas formations in complex structural areas.This experiment,it is shown shows that the wide field electromagnetic method with a large depth and high precision can effectively characterize the electrical characteristics of deep,low-resistivity thin layers in complex structural areas,and a new set of low-cost evaluation technologies for shale gas target layers based on the wide field electromagnetic method is explored.

Cone-shaped source characteristics and inductance effect of transient electromagnetic method

Small multi-turn coil devices are used with the transient electromagnetic method （TEM） in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and for detecting shallow geological targets in environmental and engineering fields. However, the equipment involved has strong mutual inductance coupling, which causes a lengthy turn-off time and a deep “blind zone”. This study proposes a new transmitter device with a conical-shape source and derives the radius formula of each coil and the mutual inductance coefficient of the cone. According to primary field characteristics, results of the two fields created, calculation of the conical-shaped source in a uniform medium using theoretical analysis, and a comparison of the inductance of the new device with that of the multi-turn coil, show that inductance of the multi-turn coil is nine times greater than that of the conical source with the same equivalent magnetic moment of 926.1 A·m2. This indicates that the new source leads to a much shallower “blind zone.” Furthermore, increasing the bottom radius and turn of the cone creates a larger mutual inductance but increasing the cone height results in a lower mutual inductance. Using the superposition principle, the primary and secondary magnetic fields for a conical source in a homogeneous medium are calculated; results indicate that the magnetic behavior of the cone is the same as that of the multi-turn coils, but the transient responses of the secondary field and the total field are more stronger than those of the multi-turn coils. To study the transient response characteristics using a cone-shaped source in a layered earth, a numerical filtering algorithm is then developed using the fast Hankel transform and the improved cosine transform, again using the superposition principle. During development, an average apparent resistivity inverted from the induced electromotive force using each coil is defined to represent the comprehensive resistivity of the conical source. To verify the forward calculation method, the transient responses of H type models and KH type models are calculated, and data are inverted using a “smoke ring” inversion. The results of inversion have good agreement with original models and show that the forward calculation method is effective. The results of this study provide an option for solving the problem of a deep “blind zone” and also provide a theoretical indicator for further research.

Three-dimensional tensor controlled-source electromagnetic modeling based on the vector finite-element method

Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.

A seafloor electromagnetic receiver for marine magnetotellurics and marine controlled-source electromagnetic sounding

In planning and executing marine controlled-source electromagnetic methods, seafloor electromagnetic receivers must overcome the problems of noise, clock drift, and power consumption. To design a receiver that performs well and overcomes the abovementioned problems, we performed forward modeling of the E-field abnormal response and established the receiver＇s characteristics. We describe the design optimization and the properties of each component, that is, low-noise induction coil sensor, low-noise Ag/AgCI electrode, low-noise chopper amplifier, digital temperature-compensated crystal oscillator module, acoustic telemetry modem, and burn wire system. Finally, we discuss the results of onshore and offshore field tests to show the effectiveness of the developed seafloor electromagnetic receiver and its performance： typical E-field noise of 0.12 nV/m/rt（Hz） at 0.5 Hz, dynamic range higher than 120 dB, clock drift lower than 1 ms/day, and continuous operation of at least 21 days.

Singularity-free Green's function for EM sources embedded in a stratified medium

We present a method to unify the calculation of Green＇s functions for an electromagnetic（EM） transmitting source embedded in a homogeneous stratified medium.A virtual interface parallel to layer interfaces is introduced through the source location.The potentials for Green＇s function are derived by decomposing the partial wave solutions to Helmholtz＇s equations into upward and downward within boundaries.The amplitudes of the potentials in each stratum are obtained recursively from the initial amplitudes at the source level.The initial amplitudes are derived by coupling with the transmitting sources and following the discontinuity of the tangential electric and magnetic fields at the source interface.Only the initial terms are related to the transmitting sources and thus need to be modified for different transmitters,whereas the kernel connected with the stratified media stays unchanged.Hence,the present method can be easily applied to EM transmitting sources with little modification.The application of the proposed method to the marine controlled-source electromagnetic method（MCSEM） demonstrates its simplicity and flexibility.

Field testing of the surface electromagnetic prospecting system

To test the performance of the Chinese whole-surface electromagnetic prospecting （SEP） system, system integrations, instrument performances, and large-scale production viabilities in Liaoning province and Inner Mongolia were measured via extensive field tests. Resultant electric fields, magnetic fields, apparent resistivities, impedance phases, and inversion profiles compared favorably with results of commercial equipment from other countries. The inversion results agreed well with the geologic information from boreholes. Field tests showed that the SEP system is stable, reliable, lightweight, and easy to operate, making it suitable and ready for real-field exploration.

Numerical analysis of multicomponent responses of surface-hole transient electromagnetic method

We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.

Forward modeling of marine DC resistivity method for a layered anisotropic earth

Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.

Microstructures and formation mechanism of headstand pyrocarbon cones developed by electromagnetic-field-assisted CVD

Novel headstand pyrocarbon cones （HPCs） with hollow structure were developed on the surfaces of pyrocarbon layers of the carbon/carbon （C/C） composites at 650-750 °C by the electromagnetic-field-assisted chemical vapor deposition in the absence of catalysts. The fine microstructures of the HPCs were characterized by high-resolution transmission electron microscopy. The results show that the textural features of the HPCs directly transfer from turbostratic structure in roots to a well-ordered high texture in stems. And the degree of high texture ordering decreases gradually from the stem to the tail of the HPCs. The formation mechanism of the HPCs was inferred as the comprehensive effect of polarization induction on electromagnetic fields and particle-filler property under disruptive discharge.

PC-based artif icial neural network inversion for airborne time-domain electromagnetic data

Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.

Research on 3D marine electromagnetic interferometry with synthetic sources for suppressing the airwave interference

In order to suppress the airwave noise in marine controlled-source electromagnetic （CSEM） data, we propose a 3D deconvolution （3DD） interferometry method with a synthetic aperture source and obtain the relative anomaly coefficient （RAC） of the EM field reflection responses to show the degree for suppressing the airwave. We analyze the potential of the proposed method for suppressing the airwave, and compare the proposed method with traditional methods in their effectiveness. A method to select synthetic source length is derived and the effect of the water depth on RAC is examined via numerical simulations. The results suggest that 3DD interferometry method with a synthetic source can effectively suppress the airwave and enhance the potential of marine CSEM to hydrocarbon exploration.

The analysis on IP signals in TEM response based on SVD

During transient electromagnetic method （TEM） exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon （a voltage change from positive to negative） caused by the influence of the induced polarization （IP） effect, which affects the depth and precision of the TEM detection. The conventional inversion method is inefficient because it is difficult to process the data. In this paper, the Cole-Cole model is adopted to analyze the effect of Dc resistivity, chargeability, time constant, and frequency exponent on the TEM response in an homogeneous half space model. Singular Value Decomposition （SVD） is used to invert the measured TEM data, and the Dc resistivity, chargeability, time constant and frequency exponent were extracted from the measured TEM data in the mine area. The extracted parameters are used for interpreting the detection result as a supplement. This reveals why the TEM data acquired in the area has a low resolution. It was found that the DC resistivity and time constant do not significantly change the results, however, the chargeability and frequency exponent have a significant effect. Because of these influences, the SVD method is more accurate than the conventional method in the apparent resistivity profile. The area of the copper mine is confined accurately based on the SVD inverted data. The conclusion has been verified by drill and is identical to the practical geological situation.

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 drill- hole 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.

Finite element solution based on fast numerical technique for large-scale electromagnetic computation

A numerical technique of the target-region locating （TRL） solver in conjunction with the wave-front method is presented for the application of the finite element method （FEM） for 3-D electromagnetic computation. First, the principle of TRL technique is described. Then, the availability of TRL solver for nonlinear application is particularly discussed demonstrating that this solver can be easily used while still remaining great efficiency. The implementation on how to apply this technique in FEM based on magnetic vector potential （MVP） is also introduced. Finally, a numerical example of 3-D magnetostatic modeling using the TRL solver and FEMLAB is given. It shows that a huge computer resource can be saved by employing the new solver.

EVALUATION OF ELECTROMAGNETIC SCATTERING OF TWO DIMENSIONAL CONDUCTING OBJECTS--TE CASE

A hybrid technique is developed for the evaluation of two dimensional electromagnetic scattering from electrically large conducting bodies with cracks on their surfaces (TE case). The edge based finite element method (FEM) is employed to compute the scattering from the cracks. Physical optics (PO) and physical theory of diffraction (PTD) are utilized to evaluate the scattering from the large bodies with the cracks filled with perfect conductors. These two methods are combined by an efficient coupling scheme. Some of numerical results are presented. It is shown that the hybrid technique has some advantages over other methods in regard to saving computer memory units and CPU time.

Static characteristics of a novel flux-switching permanent magnet linear motor

A novel flux-switching permanent magnet linear motor（FSPMLM） is proposed for linear direct driving machine tools.First,the two-and three-dimensional topological configuration of the proposed motor is presented;the basic operational principle of the FSPMLM is introduced;and the magnetic fields at the two typical conditions of no-load are analyzed.Secondly,the FSPMLM is analyzed by the two-dimensional finite element method（FEM） to investigate the static electromagnetic characteristics such as flux-linkage,back EMF（electromotive force） and inductance performances.The cogging forces of two kinds of FSPMLMs with different shaped cores are analyzed and compared,and the results show that the cogging force is significantly reduced by using the E-shaped cores.Additionally,based on the co-energy method,the thrust equation is derived and verified by the simulation results obtained by the FEM.Finally,an experimental prototype is used to test the characteristics under open circuit and load conditions.The simulation and experimental results indicate that the proposed motor has advantages of a sinusoidal back-EMF waveform,a small cogging effect and a high thrust density,and it is suitable for the application of linear direct driving machine tools.