Rate transient method is a recently-developed performance analysis tool specially designed for low-permeability or tight gas reservoirs. This method, theoretically based on pressure transient analysis, integrates mate...Rate transient method is a recently-developed performance analysis tool specially designed for low-permeability or tight gas reservoirs. This method, theoretically based on pressure transient analysis, integrates material balance principle and the concept of material balance pseudo-time proposed by Blansingame. With daily production data of gas well, it could be used to calculate OGIP, current formation pressure, permeability, skin factor, to identify complex geologic boundaries, to determine whether drainage boundary has been reached, to calculate drainage area and drainage radius for single well and to predict performance. It has been extensively employed in more than ten low-permeability gas fields. It proves that most problems in performance analysis for low permeability gas reservoirs could be solved by this method. Field practices show great economical benefits could be achieved by employing this method in gas field development.展开更多
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.展开更多
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 f...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.展开更多
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...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.展开更多
In this paper, the feasibility of measuring the gas volume fraction in a mixed gas-liquid flow by using an acoustic resonant spectroscopy (ARS) method in a transient way is studied theoretically and experimentally. Fi...In this paper, the feasibility of measuring the gas volume fraction in a mixed gas-liquid flow by using an acoustic resonant spectroscopy (ARS) method in a transient way is studied theoretically and experimentally. Firstly, the effects of sizes and locations of a single air bubble in a cylindrical cavity with two open ends on resonant frequencies are investigated numerically. Then, a transient measurement system for ARS is established, and the trends of the resonant frequencies (RFs) and resonant amplitudes (RAs) in the cylindrical cavity with gas flux inside are investigated experimentally. The measurement results by the proposed transient method are compared with those by steady-state ones and numerical ones. The numerical results show that the RFs of the cavity are highly sensitive to the volume of the single air bubble. A tiny bubble volume perturbation may cause a prominent RF shift even though the volume of the air bubble is smaller than 0.1% of that of the cavity. When the small air bubble moves, the RF shift will change and reach its maximum value as it is located at the middle of the cavity. As the gas volume fraction of the two-phase flow is low, both the RFs and RAs from the measurement results decrease dramatically with the increasing gas volume, and this decreasing trend gradually becomes even as the gas volume fraction increases further. These experimental results agree with the theoretical ones qualitatively. In addition, the transient method for ARS is more suitable for measuring the gas volume fraction with randomness and instantaneity than the steady-state one, because the latter could not reflect the random and instant characteristics of the mixed fluid due to the time consumption for frequency sweeping. This study will play a very important role in the quantitative measurement of the gas volume fraction of multiphase flows.展开更多
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 c...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.展开更多
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 use...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.展开更多
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 time...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.展开更多
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.展开更多
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 inf...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.展开更多
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.展开更多
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).展开更多
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 metho...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.展开更多
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 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.展开更多
In this paper,transient electromagnetic method was used to carry out the feasibility study on the detection and recognition of chamber blasting misfire.Firstly,an electromagnetic background field was established in th...In this paper,transient electromagnetic method was used to carry out the feasibility study on the detection and recognition of chamber blasting misfire.Firstly,an electromagnetic background field was established in the test;secondly,a benign conductor was preset in the chamber,and then the background field was eliminated after the electromagnetic field was measured;thirdly,the transient electromagnetic field was measured again after blasting;at last,the chamber blasting misfire was detected and recognized by comparing the change of eddy current field of the preset benign conductor before and after blasting.The test results showed that:When the buried depth of aluminum box target was no more than 30 m,transient electromagnetic method can clearly identify the position of the aluminum box;when the buried depth of aluminum box was more than30 m,the buried depth and position of the aluminum box was not sure due to the unknown level of secondary eddy current field generated by aluminum box.展开更多
The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this pap...The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.展开更多
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 pres...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.展开更多
文摘Rate transient method is a recently-developed performance analysis tool specially designed for low-permeability or tight gas reservoirs. This method, theoretically based on pressure transient analysis, integrates material balance principle and the concept of material balance pseudo-time proposed by Blansingame. With daily production data of gas well, it could be used to calculate OGIP, current formation pressure, permeability, skin factor, to identify complex geologic boundaries, to determine whether drainage boundary has been reached, to calculate drainage area and drainage radius for single well and to predict performance. It has been extensively employed in more than ten low-permeability gas fields. It proves that most problems in performance analysis for low permeability gas reservoirs could be solved by this method. Field practices show great economical benefits could be achieved by employing this method in gas field development.
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.41564001 and 41572185)the Natural Science Foundation of Jiangxi Province(No.20151BAB203045)
文摘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.
基金supported by the Young Scientists Fund of the National Natural Science Foundation of China(No.41304082)the China Postdoctoral Science Foundation(No.2016M590731)+2 种基金the Young Scientists Fund of the Natural Science Foundation of Hebei Province(No.D2014403011)the Program for Young Excellent Talents of Higher Education Institutions of Hebei Province(No.BJ2016046)the Geological survey project of China Geological Survey(No.1212011121197)
文摘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.
基金supported by the National Natural Science Foundation of China (Grant No. 40774099)the National High-tech Research & Devel-opment Program (Grant No. 2008AA06Z205)Special Ocean Com-monweal Research Program (Grant No. 200805005)
文摘In this paper, the feasibility of measuring the gas volume fraction in a mixed gas-liquid flow by using an acoustic resonant spectroscopy (ARS) method in a transient way is studied theoretically and experimentally. Firstly, the effects of sizes and locations of a single air bubble in a cylindrical cavity with two open ends on resonant frequencies are investigated numerically. Then, a transient measurement system for ARS is established, and the trends of the resonant frequencies (RFs) and resonant amplitudes (RAs) in the cylindrical cavity with gas flux inside are investigated experimentally. The measurement results by the proposed transient method are compared with those by steady-state ones and numerical ones. The numerical results show that the RFs of the cavity are highly sensitive to the volume of the single air bubble. A tiny bubble volume perturbation may cause a prominent RF shift even though the volume of the air bubble is smaller than 0.1% of that of the cavity. When the small air bubble moves, the RF shift will change and reach its maximum value as it is located at the middle of the cavity. As the gas volume fraction of the two-phase flow is low, both the RFs and RAs from the measurement results decrease dramatically with the increasing gas volume, and this decreasing trend gradually becomes even as the gas volume fraction increases further. These experimental results agree with the theoretical ones qualitatively. In addition, the transient method for ARS is more suitable for measuring the gas volume fraction with randomness and instantaneity than the steady-state one, because the latter could not reflect the random and instant characteristics of the mixed fluid due to the time consumption for frequency sweeping. This study will play a very important role in the quantitative measurement of the gas volume fraction of multiphase flows.
基金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.
基金support received from the National Basic Research Program of China (No2007CB209400)the National Natural Science Foundation of China (No50774085)the Young Scientists Fund of the School Science Foundation of CUMT (No2008A046)
文摘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.
基金Project 40674074 supported by the National Natural Science Foundation of China20050290501 by the Specialized Research Fund for the Doctoral Programof Higher EducationD200409 by the Scientific Research Fund for Youth of China University of Mining & Technology
文摘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.
基金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.
基金Projects(40804027,41074085) supported by the National Natural Science Foundation of ChinaProject(09JJ3048) supported by the Natural Science Foundation of Hunan Province,ChinaProject(200805331082) supported by the Research Fund for the Doctoral Program of Higher Education,China
文摘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.
基金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 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).
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Special Fund for Basic Scientific Research of Central Colleges of Chang’an University, China (Grant No. CHD2011JC080)
文摘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.
基金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.
基金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.
文摘In this paper,transient electromagnetic method was used to carry out the feasibility study on the detection and recognition of chamber blasting misfire.Firstly,an electromagnetic background field was established in the test;secondly,a benign conductor was preset in the chamber,and then the background field was eliminated after the electromagnetic field was measured;thirdly,the transient electromagnetic field was measured again after blasting;at last,the chamber blasting misfire was detected and recognized by comparing the change of eddy current field of the preset benign conductor before and after blasting.The test results showed that:When the buried depth of aluminum box target was no more than 30 m,transient electromagnetic method can clearly identify the position of the aluminum box;when the buried depth of aluminum box was more than30 m,the buried depth and position of the aluminum box was not sure due to the unknown level of secondary eddy current field generated by aluminum box.
基金support by the National Natural Science Foundation of China(grant No.51779033,51409038)the National Key Research and Development Plan(grant No.2016YFB0201001)the National Natural Science Foundation of China(grant No.51421064)
文摘The Non-uniform rational B-spline （NURBS） enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)the Innovation Fund for Graduate Student of Shanghai University of China (Grant No.SHUCX120125)
文摘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.
