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.展开更多
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.展开更多
Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorit...Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorithms for them were described. Based on the known algorithms in DC resistivity on wave-number distribution and selection, we proposed a principle on how to choose the least wave-number concerning the central-loop transient electromagnetic method. First, observe the behavior of transformation function curve with regard to wave-number in Fourier domain. In the light of its asymptote, ascertain the coverage scope of wave-number. Compared with analytic solution, the least wave-number in Fourier domain can be derived. Furthermore, the Laplace numerical inversion algorithm which needs only a few Laplace variables in pure real domain was also introduced here. The procedure was applied to forward modeling on transient electromagnetic field of a vertical magnetic dipole over uniform half-space to demonstrate them effectiveness and general applicability.展开更多
基金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.
基金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.
基金Project(40344022) supported by the National Natural Science Foundation of China
文摘Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorithms for them were described. Based on the known algorithms in DC resistivity on wave-number distribution and selection, we proposed a principle on how to choose the least wave-number concerning the central-loop transient electromagnetic method. First, observe the behavior of transformation function curve with regard to wave-number in Fourier domain. In the light of its asymptote, ascertain the coverage scope of wave-number. Compared with analytic solution, the least wave-number in Fourier domain can be derived. Furthermore, the Laplace numerical inversion algorithm which needs only a few Laplace variables in pure real domain was also introduced here. The procedure was applied to forward modeling on transient electromagnetic field of a vertical magnetic dipole over uniform half-space to demonstrate them effectiveness and general applicability.
