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...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.展开更多
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 com...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.展开更多
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 coalm...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.展开更多
In recent years,in order to meet the practical needs of deep edge mine detection with large depth and high precision,transient electromagnetic method(TEM)near emission source detection mode has become an international...In recent years,in order to meet the practical needs of deep edge mine detection with large depth and high precision,transient electromagnetic method(TEM)near emission source detection mode has become an international advanced method(Xue et al.,2020).展开更多
Herein,a three-dimensional(3D)inversion method in the frequency domain based on a time–frequency transformation was developed to improve the efficiency of the 3D inversion of transient electromagnetic(TEM)data.The Fo...Herein,a three-dimensional(3D)inversion method in the frequency domain based on a time–frequency transformation was developed to improve the efficiency of the 3D inversion of transient electromagnetic(TEM)data.The Fourier transform related to the electromagnetic response in the frequency and time domains becomes a sine or cosine transform under the excitation of downward-step current.We established a transformation matrix based on the digital fi ltering calculation for the sine transform,and then the frequency domain projection of the TEM data was determined from the linear transformation system using the smoothing constrained least squares inversion method,in which only the imaginary part was used to maintain the TEM data transformation equivalence in the bidirectional projection.Thus,the time-domain TEM inversion problem was indirectly and effectively solved in the frequency domain.In the 3D inversion of the transformed frequency-domain data,the limited-memory Broyden–Fletcher–Goldfarb–Shannoquasi–Newton(L-BFGS)method was used and modifi ed with a restart strategy to adjust the regularization parameter when the algorithm tended to a local minimum.Synthetic data tests showed that our domain transformation method can stably project the TEM data into the frequency domain with very high accuracy;furthe rmore,the 3D inversion of the transformed frequency-domain data is stable,can be used to recover the real resistivity model with an acceptable effi ciency.展开更多
In recent years, the anisotropic study has become a hot topic in the field of electromagnetics. Currently, inversion technologies of transient electromagnetic sounding data are mainly based on the case of an isotropic...In recent years, the anisotropic study has become a hot topic in the field of electromagnetics. Currently, inversion technologies of transient electromagnetic sounding data are mainly based on the case of an isotropic medium. However, the actual underground electrical structure tends to be complicated and anisotropic. It is often found that the isotropic inversion technologies do not lead to good results for field transient electromagnetic sounding data. We have developed an algorithm for calculating the transient electromagnetic response in a layered medium with azimuthal anisotropy. An occam inversion algorithm has also been implemented to invert the transient electromagnetic data induced by a grounded horizontal electric dipole in a layered medium with azimuthal anisotropy. Synthetic examples demonstrate the stability and validity of the inversion algorithm. Experimental results show different data for inverting have great influence on the inversion results.展开更多
Water inrush disasters poses a great threat to the safe exploitation of coal resources.To solve this problem,the transient electromagnetic method(TEM)was proposed to accurately detect the water accumulation in the goa...Water inrush disasters poses a great threat to the safe exploitation of coal resources.To solve this problem,the transient electromagnetic method(TEM)was proposed to accurately detect the water accumulation in the goaf.The electromagnetic response characteristics of diferent water-flled goaves were studied by electromagnetic feld theory,numerical simulation and feld verifcation.Through the models of 100%water accumulation,50%water accumulation,0%water accumulation,100%water accumulation with collapsed rock,50%water accumulation with collapsed rock and 0%water accumulation with collapsed rock goaf,the characteristics of induced voltage attenuation curves were studied.Meanwhile,the relationship between the attenuation voltage value and area of the transmitting coil,the depth of the goaf,the background resistivity,and the delay time were also simulated.The results illustrate that the attenuation curve of induced voltage presented a regular exponential decay form in the 0%water accumulation model but existed abnormal exaltation for voltage in water-flled model.Through the linear ftting curve,it can be seen that the abnormal intensity of the induced voltage becomes stronger as the distance between the measuring point and the center of the target decrement.Moreover,the abnormal amplitude of the induced voltage increases with the rise of the water accumulation and collapsed rock will weakly reduce the low-resistivity anomalous efect on the water-accumulated goaf.In addition,the response value of the attenuation voltage increased as the area of the transmitting coil increases,but decreased with increasing delay time and increasing background resistivity and depth of the target body.The feld detection results of the Majiliang coal mine also confrmed the theoretical analysis and the numerical simulation.展开更多
The damped least squares inversion principle is applied to the transient electromagnetic one-dimensional inversion of electrical sources,and a new model is obtained by continuously iterating the initial model,thereby ...The damped least squares inversion principle is applied to the transient electromagnetic one-dimensional inversion of electrical sources,and a new model is obtained by continuously iterating the initial model,thereby fitting the observed transient electromagnetic response,and performing one-dimensional inversion through induced electromotive force play.In this paper,in the damped least squares inversion,constraints are added to the Jacobian matrix,and simultaneous constraint equations and conventional inversion equations are solved.By weighting the constraint parameters,the difference between adjacent resistivities and layer thicknesses is minimized.Finally,K-type and H-type theoretical models were used to verify the reliability of the algorithm,and compared with the conventional transient electromagnetic damping least squares inversion.展开更多
Uncontrolled coal fires are natural disasters that may cause mineral loss and environmental damage.The traditional loop source transient electromagnetic method can effectively detect the low-resistivity region of coal...Uncontrolled coal fires are natural disasters that may cause mineral loss and environmental damage.The traditional loop source transient electromagnetic method can effectively detect the low-resistivity region of coal fires,but its detection efficiency is not so good for high-resistivity regions.In view of this limitation,a technique based on electrical source transient electromagnetics is proposed in this paper to detect high-resistivity regions in the spontaneous combustion process of coal.Considering the complex geometry of the coal fire area,an unstructured tetrahedral grid is used in this study to realize the spatial discretization of the model,and solve the electromagnetic field based on a vector finite element algorithm.Numerical analysis is used to investigate methods for detecting coal fires and the characteristics of effective anomalies are further examined to provide guidance for practical detection.展开更多
The transient electromagnetic(TEM)method becomes more urgent than ever for marine ex-ploration due to abundant resource reserves and the increasing undersea engineering construction activities,especially in the offsho...The transient electromagnetic(TEM)method becomes more urgent than ever for marine ex-ploration due to abundant resource reserves and the increasing undersea engineering construction activities,especially in the offshore exploration of mineral deposits such as Sanshandao gold mine.However,the re-search and application of TEM method in marine environment are still challenged by many problems.Such contradiction motivates our study on the coincident-loop TEM in seafloor exploration.The TEM response of coincident loops is firstly derived in the integral form,based on the potential functions in Helmholtz equ-ations for a magnetic source locating in the whole-space layered model.The frequency-domain vertical magnetic field is described as the Hankel integral with double first-order Bessel functions of first kind.Se-condly,the time-domain induced voltage is obtained by transforming the frequency-domain response through the cosine transform and then taking the derivative of time.To simultaneously solve the Hankel transform and the cosine transform,a novel algorithm is introduced by adapting the fixed-point quadra-ture and extrapolation via the Shanks transformation.Finally,a typical conductivity model for marine po-lymetallic deposit is designed to investigate the characteristic of TEM response under various conditions.Numerical results demonstrate that existence of conductive seawater causes the TEM response to increase significantly and decay slower.The air-sea reflected electromagnetic waves lead to a significantly large fake negative response(NR)in shallower seawater with depth less than 300 m.Increase in the height of loops will weaken and delay the anomaly response and shorten the observation time-window.The height of configu-ration should be no more than 100 m for shallower targets and 50 m for deeper targets,respectively.The observation time-window should cover 10-1000 ms.Increase in the radius of loops only enhances the TEM response proportionally but hardly improves the relative anomaly.The vertical resolution on the low-resistivity target approximates 20 m for the configuration considered in the study.Decreases in D.C.resistivity and chargeability cause the positive response(PR)to increase significantly and decay more ra-pidly.Meanwhile,the NR is advanced and enlarged significantly and decays slower compared with the PR.The influence of time constant is not monotony and there exists an optimal value for producing the maxi-mum NR.As the frequency parameter increases,the PR is caused to decay more rapidly without magni-tude change and the NR is advanced and decays more rapidly with significant increase in magnitude.The influence of frequency parameter is more pronounced than that of time constant.展开更多
Taking advantage of the water sensitivity of the transient electromagnetic method(TEM),this study assesses the effectiveness of the water curtain system for underground LPG storage caverns during the excavation period...Taking advantage of the water sensitivity of the transient electromagnetic method(TEM),this study assesses the effectiveness of the water curtain system for underground LPG storage caverns during the excavation period.It also detects fracture water flow during the excavation process in light of the practice of two pilot large underground LPG storage caverns in China.Comparative maps of apparent resistivity derived from TEM measurements before and after water-filling during the excavation process have been discussed to improve the quality of the water curtain system.This is the first case to apply TEM to detect the quality of the water curtain system during the construction of underground LPG storage cavern projects,and it is found to be practical,more visualized and worth popularizing.展开更多
The results of searching of equivalent for electromagnetic transient calculation in power system with the big number of transmission lines are shown. Brief information on the proper transmission line model is given. C...The results of searching of equivalent for electromagnetic transient calculation in power system with the big number of transmission lines are shown. Brief information on the proper transmission line model is given. Comparison of results of simulation in real and reduced system is presented taking into consideration three methods of identification. Some example are examined for different number of transmission lines in reduced system and consequently the recommendation are given for searching of equivalents of systems with many lines.展开更多
Calculations of the ElectroMagnetic(EM) response produced by a large horizontal loop placed over layered medium are rather complex because its integral expression contains the product of two Bessel functions and has a...Calculations of the ElectroMagnetic(EM) response produced by a large horizontal loop placed over layered medium are rather complex because its integral expression contains the product of two Bessel functions and has a divergent term. In this paper, an improved fast Hankel and Gaver-Stehfest transforms are introduced to solve the strong-oscillation and slow-decay properties of the integrand, where one Bessel function in the product is substituted by a carefully chosen polynomial of high accuracy and the other used as the digital filter coefficients in the convolution integral. Comparisons prove the validity and the efficiency of the proposed method.展开更多
In the present paper the main results of the pre-operational studies related to the energization maneuver of a Half-Wavelength Transmission Line (AC-Link), are presented. The simulations were performed in PSCAD with a...In the present paper the main results of the pre-operational studies related to the energization maneuver of a Half-Wavelength Transmission Line (AC-Link), are presented. The simulations were performed in PSCAD with actual Brazilian system data. Specifically the following studies were performed: electromagnetic transient studies of energization maneuver considering no fault and the occurrence of fault along the trunk;dynamic analysis of generator units that will energize the AC-Link;TRV on the circuit-breaker that will energize the trunk. The main conclusion is that the AC-Link energization experiment can be implemented without jeopardizing or causing time-life reduction of the equipment involved in the experiment, which comprises the generator unit that will energize the trunk, the step-up transformer, the arresters that will be kept during the experiments, as well as the circuit-breaker that will energize and then trip the AC-Link.展开更多
To realize the efficient transient simulation of a grid-connected power generation system based on multiple inverters, this paper proposes a hybrid simulation method integrating the models of electromagnetic transient...To realize the efficient transient simulation of a grid-connected power generation system based on multiple inverters, this paper proposes a hybrid simulation method integrating the models of electromagnetic transient and dynamic phasors. Based on a demonstration of the concepts and properties of dynamic phasors, the models of single-phase and three-phase inverters described by dynamic phasors are established first. Considering the numerical compatibility problem between dynamic phasors and instantaneous values, an interface scheme between dynamic phasors and instantaneous values is designed, and the efficiency and precision differences of various transformation methods are compared in detail.Finally, by utilizing MATLAB/Simulink, a hybrid simulation platform of a multi-inverter grid-connected system is built, and the efficiency and accuracy of the hybrid simulation are validated via comparison with the full electromagnetic transient simulation.展开更多
The calculation of apparent resistivity in electrical methods is relatively simple and consolidated in the literature. Furthermore, in the case of DC methods, the apparent resistivity values are for the most part, int...The calculation of apparent resistivity in electrical methods is relatively simple and consolidated in the literature. Furthermore, in the case of DC methods, the apparent resistivity values are for the most part, intuitive and enable an initial interpretation of the results. On the other hand, in TEM method the apparent resistivity values are not very intuitive and interpretation based only on the apparent resistivity values cannot be done reasonably. In this way, this paper presents a discussion about the reasons and effects of the negatives values of apparent resistivity in TEM soundings. The main objective of this paper is to clarify the meanings of these negatives values and the possible effects in TEM inversion.展开更多
In geostationary orbits and other quasi-vacuum environments,relativistic electron beams are affected by the initial emittance and space charge effects during the propagation process,resulting in beam quality degradati...In geostationary orbits and other quasi-vacuum environments,relativistic electron beams are affected by the initial emittance and space charge effects during the propagation process,resulting in beam quality degradation.Furthermore,axial energy distribution change in the beam and the axial transient electromagnetic effect caused by current changes in the head and tail regions of the beam also cause the beam to expand and affect its quality.In this study,the particle-in-cell method was used to construct a long-range propagation model of a relativistic electron beam in a vacuum environment.By calculating and simulating the axial energy distribution of the beam and the changes in the transient electromagnetic field,the axial effect during the propagation process was analyzed,and the parameter change law of the effective propagation of the beam was explored.This provided a theoretical reference for a more accurate assessment of the beam quality during propagation.展开更多
Garbage disposal has gradually become a key issue faced by environmental governance. What must be controlled in garbage disposal is the calculation of garbage volume. In this paper, according to the different landfill...Garbage disposal has gradually become a key issue faced by environmental governance. What must be controlled in garbage disposal is the calculation of garbage volume. In this paper, according to the different landfill conditions of a landfill in Anhui, including: covering soil, overlying HDPE membrane, overlying concrete, etc., combined with the physical differences between the rock mass and the garbage at the bottom of the landfill. The landfill covered with HDPE membrane adopts the transient electromagnetic method, the landfill covered with soil layer adopts the high-density electrical method, and the landfill covered with concrete adopts the combination of transient electromagnetic method and ground penetrating radar. At the same time, Combines the borehole data to determine that the resistivity value of the interface between the garbage and the soil was 29 Ω⋅m. Finally, the sections of multiple survey lines can accomplish three-dimensional stereoscopic that calculate the distribution and actual volume of the garbage. This method is used to determine the final investment cost.展开更多
Effective model reduction methods are required to deal with new challenges in active distribution network simulations that are on a large scale and have complicated structures.In the development of advanced electromag...Effective model reduction methods are required to deal with new challenges in active distribution network simulations that are on a large scale and have complicated structures.In the development of advanced electromagnetic transient simulation programs,automated model reduction plays an important role.This paper proposes an automated realization algorithm for the Krylov subspace based model reduction methods of an active distribution network with which the reduced model can be automatically established according to a given threshold of reduction error.The combined state-space nodal analysis framework is employed to apply the automated model reduction algorithm in popular EMTP-type simulation programs.Simulations are performed using PSCAD and a self-developed program to show the feasibility and validity of the proposed methods.展开更多
With wider applications of power electronic devices in modern power systems,simulation using traditional electro-mechanical and electromagnetic simulation tools suffer from low speed and imprecision.Multi-rate technol...With wider applications of power electronic devices in modern power systems,simulation using traditional electro-mechanical and electromagnetic simulation tools suffer from low speed and imprecision.Multi-rate technologies can greatly improve simulation efficiency by avoiding simulating the entire system using a small time-step.However,the drawbacks of the current synchronization mechanisms is that they introduce numerical errors and numerical instabilities in multi-rate parallel simulations.An improved multi-rate parallel technology,node splitting interface(NSI),is proposed to reduce errors and enhance simulation stability.A new synchronization mechanism is used to avoid prediction and signal delays.Theoretical analyses are carried out to prove the convergence and absolute stability of the proposed NSI algorithm.This algorithm is particularly suitable for simultaneously investigating long term dynamics of DC grids and fast transients of power electronic converters.展开更多
文摘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.
基金This work was supported by China Postdoctoral Science Foundation(No.2022M723391)the Science and Technology Innovation Project of Higher Education in Shanxi Province(No.2019L0754)+1 种基金the Central Guiding Local Science and Technology Development Fund Project(No.YDZJSX2021B021)Shanxi Province Basic Research Plan General Project(No.202203021221294).
文摘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.
基金supported by the Institute of Seismology Foundation, China Earthquake Administration (201326126)
文摘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.
基金project supported by Science and Technology Innovation Fund(Grant No.KDY2019001)Integrated Geophysical Simulation Lab of Chang’an University(Key Laboratory of Chinese Geophysical Society)
文摘In recent years,in order to meet the practical needs of deep edge mine detection with large depth and high precision,transient electromagnetic method(TEM)near emission source detection mode has become an international advanced method(Xue et al.,2020).
基金the National Key Research and Development Program of China(No.2016YFC060110403).
文摘Herein,a three-dimensional(3D)inversion method in the frequency domain based on a time–frequency transformation was developed to improve the efficiency of the 3D inversion of transient electromagnetic(TEM)data.The Fourier transform related to the electromagnetic response in the frequency and time domains becomes a sine or cosine transform under the excitation of downward-step current.We established a transformation matrix based on the digital fi ltering calculation for the sine transform,and then the frequency domain projection of the TEM data was determined from the linear transformation system using the smoothing constrained least squares inversion method,in which only the imaginary part was used to maintain the TEM data transformation equivalence in the bidirectional projection.Thus,the time-domain TEM inversion problem was indirectly and effectively solved in the frequency domain.In the 3D inversion of the transformed frequency-domain data,the limited-memory Broyden–Fletcher–Goldfarb–Shannoquasi–Newton(L-BFGS)method was used and modifi ed with a restart strategy to adjust the regularization parameter when the algorithm tended to a local minimum.Synthetic data tests showed that our domain transformation method can stably project the TEM data into the frequency domain with very high accuracy;furthe rmore,the 3D inversion of the transformed frequency-domain data is stable,can be used to recover the real resistivity model with an acceptable effi ciency.
文摘In recent years, the anisotropic study has become a hot topic in the field of electromagnetics. Currently, inversion technologies of transient electromagnetic sounding data are mainly based on the case of an isotropic medium. However, the actual underground electrical structure tends to be complicated and anisotropic. It is often found that the isotropic inversion technologies do not lead to good results for field transient electromagnetic sounding data. We have developed an algorithm for calculating the transient electromagnetic response in a layered medium with azimuthal anisotropy. An occam inversion algorithm has also been implemented to invert the transient electromagnetic data induced by a grounded horizontal electric dipole in a layered medium with azimuthal anisotropy. Synthetic examples demonstrate the stability and validity of the inversion algorithm. Experimental results show different data for inverting have great influence on the inversion results.
基金supported by the Joint Funds of National Natural Science Foundation of China and Shanxi Province(U1710258 and U1810120)Distinguished Youth Funds of National Natural Science Foundation of China(51925402)+3 种基金Ten Thousand Talent Program of China for Leading Scientists in Science,Technology and Innovation,Shanxi Science and Technology Major Project Funds(No.20201102004)Shanxi“1331 Project”Funds,Shanxi Province Key Laboratory Construction Project Funds(No.202104010910021)Shanxi-Zheda Institute of Advanced Materials and Chemical Engineering(No.2021SX-TD001,No.2021SX-TD002)National Natural Science Foundation of China(51804208).
文摘Water inrush disasters poses a great threat to the safe exploitation of coal resources.To solve this problem,the transient electromagnetic method(TEM)was proposed to accurately detect the water accumulation in the goaf.The electromagnetic response characteristics of diferent water-flled goaves were studied by electromagnetic feld theory,numerical simulation and feld verifcation.Through the models of 100%water accumulation,50%water accumulation,0%water accumulation,100%water accumulation with collapsed rock,50%water accumulation with collapsed rock and 0%water accumulation with collapsed rock goaf,the characteristics of induced voltage attenuation curves were studied.Meanwhile,the relationship between the attenuation voltage value and area of the transmitting coil,the depth of the goaf,the background resistivity,and the delay time were also simulated.The results illustrate that the attenuation curve of induced voltage presented a regular exponential decay form in the 0%water accumulation model but existed abnormal exaltation for voltage in water-flled model.Through the linear ftting curve,it can be seen that the abnormal intensity of the induced voltage becomes stronger as the distance between the measuring point and the center of the target decrement.Moreover,the abnormal amplitude of the induced voltage increases with the rise of the water accumulation and collapsed rock will weakly reduce the low-resistivity anomalous efect on the water-accumulated goaf.In addition,the response value of the attenuation voltage increased as the area of the transmitting coil increases,but decreased with increasing delay time and increasing background resistivity and depth of the target body.The feld detection results of the Majiliang coal mine also confrmed the theoretical analysis and the numerical simulation.
基金sponsored by geological Survey Project of China Geological Survey(DD20189210).
文摘The damped least squares inversion principle is applied to the transient electromagnetic one-dimensional inversion of electrical sources,and a new model is obtained by continuously iterating the initial model,thereby fitting the observed transient electromagnetic response,and performing one-dimensional inversion through induced electromotive force play.In this paper,in the damped least squares inversion,constraints are added to the Jacobian matrix,and simultaneous constraint equations and conventional inversion equations are solved.By weighting the constraint parameters,the difference between adjacent resistivities and layer thicknesses is minimized.Finally,K-type and H-type theoretical models were used to verify the reliability of the algorithm,and compared with the conventional transient electromagnetic damping least squares inversion.
基金Supported by National Key R&D Program of China(No.2018YFC0807900).
文摘Uncontrolled coal fires are natural disasters that may cause mineral loss and environmental damage.The traditional loop source transient electromagnetic method can effectively detect the low-resistivity region of coal fires,but its detection efficiency is not so good for high-resistivity regions.In view of this limitation,a technique based on electrical source transient electromagnetics is proposed in this paper to detect high-resistivity regions in the spontaneous combustion process of coal.Considering the complex geometry of the coal fire area,an unstructured tetrahedral grid is used in this study to realize the spatial discretization of the model,and solve the electromagnetic field based on a vector finite element algorithm.Numerical analysis is used to investigate methods for detecting coal fires and the characteristics of effective anomalies are further examined to provide guidance for practical detection.
基金This research was supported by the National Key R&D Program of China(No.2018YFC0603806).
文摘The transient electromagnetic(TEM)method becomes more urgent than ever for marine ex-ploration due to abundant resource reserves and the increasing undersea engineering construction activities,especially in the offshore exploration of mineral deposits such as Sanshandao gold mine.However,the re-search and application of TEM method in marine environment are still challenged by many problems.Such contradiction motivates our study on the coincident-loop TEM in seafloor exploration.The TEM response of coincident loops is firstly derived in the integral form,based on the potential functions in Helmholtz equ-ations for a magnetic source locating in the whole-space layered model.The frequency-domain vertical magnetic field is described as the Hankel integral with double first-order Bessel functions of first kind.Se-condly,the time-domain induced voltage is obtained by transforming the frequency-domain response through the cosine transform and then taking the derivative of time.To simultaneously solve the Hankel transform and the cosine transform,a novel algorithm is introduced by adapting the fixed-point quadra-ture and extrapolation via the Shanks transformation.Finally,a typical conductivity model for marine po-lymetallic deposit is designed to investigate the characteristic of TEM response under various conditions.Numerical results demonstrate that existence of conductive seawater causes the TEM response to increase significantly and decay slower.The air-sea reflected electromagnetic waves lead to a significantly large fake negative response(NR)in shallower seawater with depth less than 300 m.Increase in the height of loops will weaken and delay the anomaly response and shorten the observation time-window.The height of configu-ration should be no more than 100 m for shallower targets and 50 m for deeper targets,respectively.The observation time-window should cover 10-1000 ms.Increase in the radius of loops only enhances the TEM response proportionally but hardly improves the relative anomaly.The vertical resolution on the low-resistivity target approximates 20 m for the configuration considered in the study.Decreases in D.C.resistivity and chargeability cause the positive response(PR)to increase significantly and decay more ra-pidly.Meanwhile,the NR is advanced and enlarged significantly and decays slower compared with the PR.The influence of time constant is not monotony and there exists an optimal value for producing the maxi-mum NR.As the frequency parameter increases,the PR is caused to decay more rapidly without magni-tude change and the NR is advanced and decays more rapidly with significant increase in magnitude.The influence of frequency parameter is more pronounced than that of time constant.
基金We thank CREGC(20130503)Innovative Research Team in SWJTU–China(IRT0955)+1 种基金the Fundamental Research Funds for the Central Universities(SWJTU11ZT33)China Scholarship Council(201407000040)for their sponsorships of this project,and we also thank the consultants and engineers contacted,especially Jinzhou and Huizhou project departments of China Railway Erju Second Engineering Co.Ltd.for providing their data during this study.
文摘Taking advantage of the water sensitivity of the transient electromagnetic method(TEM),this study assesses the effectiveness of the water curtain system for underground LPG storage caverns during the excavation period.It also detects fracture water flow during the excavation process in light of the practice of two pilot large underground LPG storage caverns in China.Comparative maps of apparent resistivity derived from TEM measurements before and after water-filling during the excavation process have been discussed to improve the quality of the water curtain system.This is the first case to apply TEM to detect the quality of the water curtain system during the construction of underground LPG storage cavern projects,and it is found to be practical,more visualized and worth popularizing.
文摘The results of searching of equivalent for electromagnetic transient calculation in power system with the big number of transmission lines are shown. Brief information on the proper transmission line model is given. Comparison of results of simulation in real and reduced system is presented taking into consideration three methods of identification. Some example are examined for different number of transmission lines in reduced system and consequently the recommendation are given for searching of equivalents of systems with many lines.
基金Supported by Chinese Government SinoProbe Plan(No.SinoProbe-09-02)R&D of Key Instruments and Technologies for Deep Resources Prospecting(No.ZDYZ 2012-1-03)
文摘Calculations of the ElectroMagnetic(EM) response produced by a large horizontal loop placed over layered medium are rather complex because its integral expression contains the product of two Bessel functions and has a divergent term. In this paper, an improved fast Hankel and Gaver-Stehfest transforms are introduced to solve the strong-oscillation and slow-decay properties of the integrand, where one Bessel function in the product is substituted by a carefully chosen polynomial of high accuracy and the other used as the digital filter coefficients in the convolution integral. Comparisons prove the validity and the efficiency of the proposed method.
基金supported CNPq,CAPES and FAPESP in Brazilthe project supported by ELETROBRAS/ELETRONORTE,ELETROBRAS/CHESF and ENTE as part of ANEEL(Brazilian Electrical Energy Regula-tory Agency)strategic project 2008(proc.4500072477).
文摘In the present paper the main results of the pre-operational studies related to the energization maneuver of a Half-Wavelength Transmission Line (AC-Link), are presented. The simulations were performed in PSCAD with actual Brazilian system data. Specifically the following studies were performed: electromagnetic transient studies of energization maneuver considering no fault and the occurrence of fault along the trunk;dynamic analysis of generator units that will energize the AC-Link;TRV on the circuit-breaker that will energize the trunk. The main conclusion is that the AC-Link energization experiment can be implemented without jeopardizing or causing time-life reduction of the equipment involved in the experiment, which comprises the generator unit that will energize the trunk, the step-up transformer, the arresters that will be kept during the experiments, as well as the circuit-breaker that will energize and then trip the AC-Link.
基金supported by the State Grid Science and Technology Project (grant no. KJ2021-069)。
文摘To realize the efficient transient simulation of a grid-connected power generation system based on multiple inverters, this paper proposes a hybrid simulation method integrating the models of electromagnetic transient and dynamic phasors. Based on a demonstration of the concepts and properties of dynamic phasors, the models of single-phase and three-phase inverters described by dynamic phasors are established first. Considering the numerical compatibility problem between dynamic phasors and instantaneous values, an interface scheme between dynamic phasors and instantaneous values is designed, and the efficiency and precision differences of various transformation methods are compared in detail.Finally, by utilizing MATLAB/Simulink, a hybrid simulation platform of a multi-inverter grid-connected system is built, and the efficiency and accuracy of the hybrid simulation are validated via comparison with the full electromagnetic transient simulation.
文摘The calculation of apparent resistivity in electrical methods is relatively simple and consolidated in the literature. Furthermore, in the case of DC methods, the apparent resistivity values are for the most part, intuitive and enable an initial interpretation of the results. On the other hand, in TEM method the apparent resistivity values are not very intuitive and interpretation based only on the apparent resistivity values cannot be done reasonably. In this way, this paper presents a discussion about the reasons and effects of the negatives values of apparent resistivity in TEM soundings. The main objective of this paper is to clarify the meanings of these negatives values and the possible effects in TEM inversion.
基金National Natural Science Foundation of China(Nos.61372050,U1730247)the HighPower Microwave Key Laboratory Foundation Program(No.6142605200301)。
文摘In geostationary orbits and other quasi-vacuum environments,relativistic electron beams are affected by the initial emittance and space charge effects during the propagation process,resulting in beam quality degradation.Furthermore,axial energy distribution change in the beam and the axial transient electromagnetic effect caused by current changes in the head and tail regions of the beam also cause the beam to expand and affect its quality.In this study,the particle-in-cell method was used to construct a long-range propagation model of a relativistic electron beam in a vacuum environment.By calculating and simulating the axial energy distribution of the beam and the changes in the transient electromagnetic field,the axial effect during the propagation process was analyzed,and the parameter change law of the effective propagation of the beam was explored.This provided a theoretical reference for a more accurate assessment of the beam quality during propagation.
文摘Garbage disposal has gradually become a key issue faced by environmental governance. What must be controlled in garbage disposal is the calculation of garbage volume. In this paper, according to the different landfill conditions of a landfill in Anhui, including: covering soil, overlying HDPE membrane, overlying concrete, etc., combined with the physical differences between the rock mass and the garbage at the bottom of the landfill. The landfill covered with HDPE membrane adopts the transient electromagnetic method, the landfill covered with soil layer adopts the high-density electrical method, and the landfill covered with concrete adopts the combination of transient electromagnetic method and ground penetrating radar. At the same time, Combines the borehole data to determine that the resistivity value of the interface between the garbage and the soil was 29 Ω⋅m. Finally, the sections of multiple survey lines can accomplish three-dimensional stereoscopic that calculate the distribution and actual volume of the garbage. This method is used to determine the final investment cost.
基金supported in part by the National Key Technology Research and Development Program of China(2013BAAOlB03)in part by the National Natural Science Foundation of China(51261130473).
文摘Effective model reduction methods are required to deal with new challenges in active distribution network simulations that are on a large scale and have complicated structures.In the development of advanced electromagnetic transient simulation programs,automated model reduction plays an important role.This paper proposes an automated realization algorithm for the Krylov subspace based model reduction methods of an active distribution network with which the reduced model can be automatically established according to a given threshold of reduction error.The combined state-space nodal analysis framework is employed to apply the automated model reduction algorithm in popular EMTP-type simulation programs.Simulations are performed using PSCAD and a self-developed program to show the feasibility and validity of the proposed methods.
基金This work was supported in part by the People Programme(Marie Curie Actions)of the European Union’s Seventh Framework Programme FP7/2007-2013/under REA grant agreement(No.317221)project title MEDOW,in part supported by the Project of National Science Foundation for Distinguished Young Scholars of China(No.51407164).
文摘With wider applications of power electronic devices in modern power systems,simulation using traditional electro-mechanical and electromagnetic simulation tools suffer from low speed and imprecision.Multi-rate technologies can greatly improve simulation efficiency by avoiding simulating the entire system using a small time-step.However,the drawbacks of the current synchronization mechanisms is that they introduce numerical errors and numerical instabilities in multi-rate parallel simulations.An improved multi-rate parallel technology,node splitting interface(NSI),is proposed to reduce errors and enhance simulation stability.A new synchronization mechanism is used to avoid prediction and signal delays.Theoretical analyses are carried out to prove the convergence and absolute stability of the proposed NSI algorithm.This algorithm is particularly suitable for simultaneously investigating long term dynamics of DC grids and fast transients of power electronic converters.