Discussion on the application prospect of the transient electromagnetic method based on the dual launcher and mine advanced detection

The dual transmitter implements the equivalent anti-magnetic flux transient electromagnetic method, which can effectively reduce the scope of the transient electromagnetic detection blind area. However, this method is rarely reported in the detection of pipelines in urban geophysical exploration and the application of coal mines. Based on this, this paper realizes the equivalent anti-magnetic flux transient electromagnetic method based on the dual launcher. The suppression effect of this method on the blind area is analyzed by physical simulation. And the detection experiment of underground pipelines is carried out outdoors. The results show that the dual launcher can significantly reduce the turn-off time, thereby effectively reducing the impact of the blind area on the detection results, and the pipeline detection results verify the device’s effectiveness. Finally, based on the ground experimental results, the application prospect of mine advanced detection is discussed. Compared with other detection fields, the formation of blind areas is mainly caused by the equipment. If the dual launcher can be used to reduce the blind area, the accuracy of advanced detection can be improved more effectively. The above research results are of great significance for improving the detection accuracy of the underground transient electromagnetic method.

A Novel Localized Meshless Method for Solving Transient Heat Conduction Problems in Complicated Domains

This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.

3D forward modeling and response characteristics of low-resistivity overburden of the CFS-PML absorbing boundary for ground-well transient electromagnetic method

This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.

Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.

Application of geophysical methods in fine detection of urban concealed karst:A case study of Wuhan City,China

The construction of modern livable cities faces challenges in karst areas,including ground collapse and engineering problems.Wuhan,with a population of 13.74×10^(6) and approximately 1161 km^(2)of soluble rocks in the urban area of 8569.15 km^(2),predominantly consists of concealed karst areas where occasional ground collapse events occur,posing significant threats to underground engineering projects.To address these challenges,a comprehensive geological survey was conducted in Wuhan,focusing on major karstrelated issues.Geophysical methods offer advantages over drilling in detecting concealed karst areas due to their efficiency,non-destructiveness,and flexibility.This paper reviewed the karst geological characteristics in Wuhan and the geophysical exploration methods for karst,selected eight effective geophysical methods for field experimentation,evaluated their suitability,and proposed method combinations for different karst scenarios.The results show that different geophysical methods have varying applicability for karst detection in Wuhan,and combining multiple methods enhances detection effectiveness.The specific recommendations for method combinations provided in this study serve as a valuable reference for karst detection in Wuhan.

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.

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.

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.

Equivalent Method of Integrated Power Generation System of Wind, Photovoltaic and Energy Storage in Power Flow Calculation and Transient Simulation

针对工程实际开展风光储联合发电系统在潮流计算和机电暂态仿真中的等值方法研究，旨在为大容量风光储联合发电系统的并网仿真分析奠定基础。将潮流计算的等值分为单元机组和集电系统2部分来研究。单元机组等值采用根据不同控制模式选取不同节点类型的方法，针对集电系统等值提出基于损耗不变原则的方法。等值模型和详细模型的算例结果表明，潮流计算等值方法具有较好的精度。在机电暂态仿真动态等值中，基于实际工程计算的最严重工况分析原则，提出运行在满出力点的单机“倍乘”等值模型，为工程计算中的风光储联合发电系统动态等值提供了一种解决方案。

Mapping coalmine goaf using transient electromagnetic method and high density resistivity method in Ordos City, China

The research about subsurface characteristics by using transient electromagnetic method（TEM） and high density resistivity method（HDRM） were already conducted in Ordos. The objective of this research is to detect coalmine goaf areas based on rock resistivity. The data processing using wavelet transform, three point smoothing, RES2 DINV and Maxwell processing software to obtain 2D resistivity structure. The results showed that the layers with maximum resistivity values（30e33 U m on Line 1, 30e31 U m on Line 2, 32e40 U m on Line3） are founded at station 1e7, and 14e20 on Line 1,13e18 on Line 2, and 8e13 and 16e20 on Line 3 which is predicted as goaf layer, and the minimum resistivity values（20e26 U m of TEM, 45e75 U m of HDRM） at the other layers. This resistivity difference was caused by the geology and characteristics of the study area which is located close by the cleugh with rich coal, so the goaf area distinguishable with aquifer layer and coal seam. The results were also significant accidents and serious destruction of ecological environment.

Transient analysis of 1D inhomogeneous media by dynamic inhomogeneous finite element method

The dynamic inhomogeneous finite element method is studied for use in the transient analysis of one dimensional inhomogeneous media. The general formula of the inhomogeneous consistent mass matrix is established based on the shape function. In order to research the advantages of this method, it is compared with the general finite element method. A linear bar element is chosen for the discretization tests of material parameters with two fictitious distributions. And, a numerical example is solved to observe the differences in the results between these two methods. Some characteristics of the dynamic inhomogeneous finite element method that demonstrate its advantages are obtained through comparison with the general finite element method. It is found that the method can be used to solve elastic wave motion problems with a large element scale and a large number of iteration steps.

A new complex variable meshless method for transient heat conduction problems

In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.

Experiment in metal disturbance during advanced detection using a transient electromagnetic method in coal mines

Tunneling machines, or excavators, are large and good conductors and affect the reliability of data gathering and interpretation in advanced detection using transient electromagnetic methods. In our experiment, we used a coincident-loop and central loop type of configuration, where the coil plane l） vertical to and 2） parallel to the working face. A SIROTEM instrument at different locations was used to observe the transient electromagnetic responses of the excavator and to analyze the response amplitudes. The result shows that the tunneling machine affects the advanced detection data and is related to the way the coil is coupled. When the excavator is 6 m from the observatory, the interference of tunneling machine can be ignored.

Prediction Technology of Buried Water-Bearing Structures in Coal Mines Using Transient Electromagnetic Method

Buried water-conducting and water-bearing structures in front of the driving head may easily lead to water bursts in coal mines. Therefore,it is very important for the safety of production to make an accurate and timely forecast about water bursts. Based on the smoke ring effect of transient electromagnetic fields,the principle of transient electro-magnetic method used in detecting buried water-bearing structures in coal mines in advance,is discussed. Small multi-turn loop configurations used in coal mines are proposed and a field procedure of semicircular sector scanning is presented. The application of this method in one coal mine indicates that the technology has many advantages compared with others. The method is inexpensive,highly accurate and efficient. Suggestions are presented for future solutions to some remaining problems.

Alterable-element Method for Vehicle-bridge Interaction Considering the Transient Jump of Wheel

The so called "alterable-element method" (AEM) was introduced to deal with the coupling interac-tion of vehicle and sub-structure considering the actual transient jump of wheel, while the classical "contact all along" assumption based on which wheels and lower structure are always contact was abandoned. The alterable element used in this method is a conceptional element, which is used to calculate the coupling interaction of upper and lower structures and has some typical characteristics: firstly it flows along with the moving of contact point; secondly whether it is used for calculation depends on the contact state; thirdly its sizes could change according to specific problems and so on. VISUAL FORTRAN program was coded, and different moving vehicle models were presented taking into consideration the effects of random corrugation in the numerical study. The numerical solutions are favored comparing with the results obtained by alternative methods when there is no jump phenomenon existed. With abrupt irregularity, the transient jump of wheel was studied using the present method.

Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems

In this paper, the complex variable reproducing kernel particle （CVRKP） method and the finite element （FE） method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.

Transient characteristics verification method for DC transformer used in flexible HVDC system

Previous studies have proposed higher requirements for the transient characteristics of a DC transformer used in a flexible high-voltage direct current(HVDC) system to achieve faster sampling speed and meet wider bandwidth requirements of the control and protection signal, and to eventually suppress the large transient fault current. In this study, a transient characteristics verification method is proposed for transient characteristics verification of a DC transformer used in a flexible HVDC system based on resampling technology and LabVIEW measurement technology after analyzing the key technology for transient characteristics verification of a DC transformer. A laboratory experiment for the transient characteristics of a full-fiber electronic DC transformer is conducted, and experimental results show that such verification method can be employed for frequency response and step response verification of a DC transformer at 10% of the rated voltage and current, and can eventually improve the screening of a DC transformer.

Calculation of all-time apparent resistivity of large loop transient electromagnetic method with very fast simulated annealing

In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.

New nonconforming finite element method for solving transient Naiver-Stokes equations

For transient Naiver-Stokes problems, a stabilized nonconforming finite element method is presented, focusing on two pairs inf-sup unstable finite element spaces, i.e., pNC/pNC triangular and pNQ/pNQ quadrilateral finite element spaces. The semi- and full-discrete schemes of the stabilized method are studied based on the pressure projection and a variational multi-scale method. It has some attractive features： avoiding higher-order derivatives and edge-based data structures, adding a discrete velocity term only on the fine scale, being effective for high Reynolds number fluid flows, and avoiding increased computation cost. For the full-discrete scheme, it has second-order estimations of time and is unconditionally stable. The presented numerical results agree well with the theoretical results.

Distributed Lagrange Multiplier/Fictitious Domain Finite Element Method for a Transient Stokes Interface Problem with Jump Coefficients

The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.