In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is a...In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.展开更多
Introduction The East Kunlun Orogenic Belt is located in the northeastern part of the Qinghai–Tibet Plateau(Li et al.,2007),stretching from the East Kunlun to the Elashan area in an east–west direction(Guo et al.,20...Introduction The East Kunlun Orogenic Belt is located in the northeastern part of the Qinghai–Tibet Plateau(Li et al.,2007),stretching from the East Kunlun to the Elashan area in an east–west direction(Guo et al.,2018).It is an important part of the Central Orogenic Belt(Xiong et al.,2023).It is considered one of the important gold mineralization regions in the Tethys tectonic domain(Norbu et al.,2023)and an essential potential base for mineral resources in China.Wulonggou and Gouli gold mines have been discovered successively,earning the reputation of the"Golden Belt of Qinghai Province"(Feng et al.,2004;He et al.,2023).展开更多
To fully exploit the technical advantages of the large-depth and high-precision artificial source electromagnetic method in the complex structure area of southern Sichuan and compensate for the shortcomings of the con...To fully exploit the technical advantages of the large-depth and high-precision artificial source electromagnetic method in the complex structure area of southern Sichuan and compensate for the shortcomings of the conventional electromagnetic method in exploration depth,precision,and accuracy,the large-depth and high-precision wide field electromagnetic method is applied to the complex structure test area of the Luochang syncline and Yuhe nose anticline in the southern Sichuan.The advantages of the wide field electromagnetic method in detecting deep,low-resistivity thin layers are demonstrated.First,on the basis of the analysis of physical property data,a geological–geoelectric model is established in the test area,and the wide field electromagnetic method is numerically simulated to analyze and evaluate the response characteristics of deep thin shale gas layers on wide field electromagnetic curves.Second,a wide field electromagnetic test is conducted in the complex structure area of southern Sichuan.After data processing and inversion imaging,apparent resistivity logging data are used for calibration to develop an apparent resistivity interpretation model suitable for the test area.On the basis of the results,the characteristics of the electrical structure change in the shallow longitudinal formation of 6 km are implemented,and the transverse electrical distribution characteristics of the deep shale gas layer are delineated.In the prediction area near the well,the subsequent data verification shows that the apparent resistivity obtained using the inversion of the wide field electromagnetic method is consistent with the trend of apparent resistivity revealed by logging,which proves that this method can effectively identify the weak response characteristics of deep shale gas formations in complex structural areas.This experiment,it is shown shows that the wide field electromagnetic method with a large depth and high precision can effectively characterize the electrical characteristics of deep,low-resistivity thin layers in complex structural areas,and a new set of low-cost evaluation technologies for shale gas target layers based on the wide field electromagnetic method is explored.展开更多
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.展开更多
The signs of the electric field markers in Figs.2 and 4 of the paper[Chin.Phys.B 32104211(2023)]have been corrected.These modifications do not affect the results derived in the paper.
Offshore carbon dioxide(CO_(2)) storage is an effective method for reducing greenhouse gas emissions. However, when using traditional seismic wave methods to monitor the migration of sequestration CO_(2) plumes, the c...Offshore carbon dioxide(CO_(2)) storage is an effective method for reducing greenhouse gas emissions. However, when using traditional seismic wave methods to monitor the migration of sequestration CO_(2) plumes, the characteristics of wave velocity changes tend to become insignificant beyond a certain limit. In contrast, the controllable source electromagnetic method(CSEM) remains highly sensitive to resistivity changes. By simulating different CO_(2) plume migration conditions, we established the relevant models and calculated the corresponding electric field response characteristic curves, allowing us to analyze the CSEM's ability to monitor CO_(2) plumes. We considered potential scenarios for the migration and diffusion of offshore CO_(2) storage, including various burial depths, vertical extension diffusion, lateral extension diffusion,multiple combinations of lateral intervals, and electric field components. We also obtained differences in resistivity inversion imaging obtained by CSEM to evaluate its feasibility in monitoring and to analyze all the electric field(Ex, Ey, and Ez) response characteristics. CSEM has great potential in monitoring CO_(2) plume migration in offshore saltwater reservoirs due to its high sensitivity and accuracy. Furthermore, changes in electromagnetic field response reflect the transport status of CO_(2) plumes, providing an important basis for monitoring and evaluating CO_(2)transport behavior during storage processes.展开更多
Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellit...Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellites,inverting the Earth’s three-dimensional conductivity distribution on a global scale becomes attainable.A key requirement in the global conductivity inversion is to have a forward solver with high-accuracy and efficiency.In this study,a finite volume method for global electromagnetic induction forward modeling is developed based on unstructured grids.Arbitrary polyhedral grids are supported in our algorithms to obtain high geometric adaptability.We employ a cell-centered collocated variable arrangement which allows convenient discretization for complex geometries and straightforward implementation of multigrid technique.To validate the method,we test our code with two synthetic models and compare our finite volume results with an analytical solution and a finite element numerical solution.Good agreements are observed between our solution and other results,indicating acceptable accuracy of the proposed method.展开更多
The construction of modern livable cities faces challenges in karst areas,including ground collapse and engineering problems.Wuhan,with a population of 13.74×10^(6) and approximately 1161 km^(2)of soluble rocks i...The construction of modern livable cities faces challenges in karst areas,including ground collapse and engineering problems.Wuhan,with a population of 13.74×10^(6) and approximately 1161 km^(2)of soluble rocks in the urban area of 8569.15 km^(2),predominantly consists of concealed karst areas where occasional ground collapse events occur,posing significant threats to underground engineering projects.To address these challenges,a comprehensive geological survey was conducted in Wuhan,focusing on major karstrelated issues.Geophysical methods offer advantages over drilling in detecting concealed karst areas due to their efficiency,non-destructiveness,and flexibility.This paper reviewed the karst geological characteristics in Wuhan and the geophysical exploration methods for karst,selected eight effective geophysical methods for field experimentation,evaluated their suitability,and proposed method combinations for different karst scenarios.The results show that different geophysical methods have varying applicability for karst detection in Wuhan,and combining multiple methods enhances detection effectiveness.The specific recommendations for method combinations provided in this study serve as a valuable reference for karst detection in Wuhan.展开更多
In recent years,karst construction projects in the built-up area of Wuhan(capital of Hubei Province,China)are increasing,and the karst geological disasters have aroused social concerns.The actual engineering projects ...In recent years,karst construction projects in the built-up area of Wuhan(capital of Hubei Province,China)are increasing,and the karst geological disasters have aroused social concerns.The actual engineering projects usually use shallow geophysical exploration methods to explore karst.This paper uses Spatial Auto-Correlation Method(SPAC)and electromagnetic Computerized Tomography(CT)to detect karst in urban built-up areas.Depending on the different physical properties of rock and soil,the SPAC method can better reveal the interface between soil and rock strata and the interface between soil layers.The electromagnetic CT method can identify strata according to the apparent absorption coefficient,which can better reveal the interface between soil and rock,the interface between the more intact and weathered rock.The SPAC method is mainly qualitative to measure the low-speed area,namely,the wrong geological body i.e.,karst cave,but also can detect the fracture zone or filling mode of karst cave,and at the same time,cannot use exploration holes or logging observation.The electromagnetic CT method can accurately detect the location and scale of the karst caves and has a higher accuracy detecting karst bands.In addition,exploration holes or well logging observations are also expected to be conducted,and their detection effect is greatly affected by lithology.展开更多
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.展开更多
Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as eleme...Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method (EFGM), the point interpolation method (PIM), and the radial point interpolation method (RPIM). Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled- source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG-PIM and EFG-RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM.展开更多
Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-freque...Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-frequency electromagnetic fields on the environment and which are capable of modifying and destroying our ecosystem. Particular interest was given in this article due to a massive influx of population living near high voltage lines. The analysis and simulation of the influence of low frequency electromagnetic fields on living beings in the vicinity of high voltage sources 132 kV and 220 kV in urban areas in DR Congo is the subject of our research with a view to estimating the level of exposure of humans to low frequency electromagnetic fields. To carry out our research, we used the classic method of analyzing the field produced near a high voltage line based on Maxwell’s image theory, the Maxwell-Gauss theorem and Maxwell-Ampère theorem to model and quantify low-frequency electromagnetic fields in the vicinity of high-voltage lines. The 2D FDTD numerical formulation was developed from telegraphers’ equations and allowed us to obtain models of current and voltage induced by electromagnetic fields on living beings below and near HV lines. The different simulations carried out on the proposed models illustrate the effects of the electrical and geometric parameters of the pylons on the distribution of the electromagnetic field in the vicinity of the HV lines. The results obtained were compared to the safety limits recommended by the standards.展开更多
A class of finite step iterative methods, conjugate gradients, for the solution of an operator equation, is presented on this paper to solve electromagnetic scattering. The method of generalized equivalent circuit is ...A class of finite step iterative methods, conjugate gradients, for the solution of an operator equation, is presented on this paper to solve electromagnetic scattering. The method of generalized equivalent circuit is used to model the problem and then deduce an electromagnetic equation based on the impedance operator. Four versions of the conjugate gradient method are presented and numerical results for an iris structure are given, to illustrate convergence properties of each version. Computational efficiency of these methods has been compared to the moment method.展开更多
The quality of the low frequency electromagnetic data is affected by the spike and the trend noises.Failure in removal of the spikes and the trends reduces the credibility of data explanation.Based on the analyses of ...The quality of the low frequency electromagnetic data is affected by the spike and the trend noises.Failure in removal of the spikes and the trends reduces the credibility of data explanation.Based on the analyses of the causes and characteristics of these noises,this paper presents the results of a preset statistics stacking method(PSSM)and a piecewise linear fitting method(PLFM)in de-noising the spikes and trends,respectively.The magnitudes of the spikes are either higher or lower than the normal values,which leads to distortion of the useful signal.Comparisons have been performed in removing of the spikes among the average,the statistics and the PSSM methods,and the results indicate that only the PSSM can remove the spikes successfully.On the other hand,the spectrums of the linear and nonlinear trends mainly lie in the low frequency band and can change the calculated resistivity significantly.No influence of the trends is observed when the frequency is higher than a certain threshold value.The PLSM can remove effectively both the linear and nonlinear trends with errors around 1% in the power spectrum.The proposed methods present an effective way for de-noising the spike and the trend noises in the low frequency electromagnetic data,and establish a research basis for de-noising the low frequency noises.展开更多
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.展开更多
We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method...We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized: the real number of the electromagnetic field solution based on the Gaver-Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the "ground excitation-stratum measurement method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drill- hole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.展开更多
Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is d...Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary(scattered) fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.展开更多
In an effort to reduce the shale gas exploration risks and costs, we applied the wide-field electromagnetic method (WFEM), because of its strong anti-interference capability, high resolution, ability to conduct expl...In an effort to reduce the shale gas exploration risks and costs, we applied the wide-field electromagnetic method (WFEM), because of its strong anti-interference capability, high resolution, ability to conduct exploration at large depths, and high efficiency, to the Bayan Syncline in the South Huayuan block, Hunan Province. We collected rock samples and analyzed their resistivity and induced polarization (IP) and built A series of two-dimensional models for geological conditions to investigate the applicability of WFEM to different geological structures. We also analyzed the correlation between TOC of shale and the resistivity and IP ratio to determine the threshold for identifying target formations. We used WFEM to identify the underground structures and determine the distribution, depth, and thickness of the target strata. Resistivity, IP, and total organic carbon were used to evaluate the shale gas prospects and select favorable areas (sweet spots) for exploration and development. Subsequently, drilling in these areas proved the applicability of WFEM in shale gas exploration.展开更多
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.展开更多
基金supported by the Young Scientists Fund of the National Natural Science Foundation of China(No.62102444)a Major Research Project in Higher Education Institutions in Henan Province(No.23A560015).
文摘In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.
基金supported by Qinghai Provincial Association for Science and Technology Youth Science and Technology Talent Support Project(Grant No.2023QHSKXRCTJ47)Exploration Foundation of Qinghai Province(Grant No.2023085029ky004)。
文摘Introduction The East Kunlun Orogenic Belt is located in the northeastern part of the Qinghai–Tibet Plateau(Li et al.,2007),stretching from the East Kunlun to the Elashan area in an east–west direction(Guo et al.,2018).It is an important part of the Central Orogenic Belt(Xiong et al.,2023).It is considered one of the important gold mineralization regions in the Tethys tectonic domain(Norbu et al.,2023)and an essential potential base for mineral resources in China.Wulonggou and Gouli gold mines have been discovered successively,earning the reputation of the"Golden Belt of Qinghai Province"(Feng et al.,2004;He et al.,2023).
文摘To fully exploit the technical advantages of the large-depth and high-precision artificial source electromagnetic method in the complex structure area of southern Sichuan and compensate for the shortcomings of the conventional electromagnetic method in exploration depth,precision,and accuracy,the large-depth and high-precision wide field electromagnetic method is applied to the complex structure test area of the Luochang syncline and Yuhe nose anticline in the southern Sichuan.The advantages of the wide field electromagnetic method in detecting deep,low-resistivity thin layers are demonstrated.First,on the basis of the analysis of physical property data,a geological–geoelectric model is established in the test area,and the wide field electromagnetic method is numerically simulated to analyze and evaluate the response characteristics of deep thin shale gas layers on wide field electromagnetic curves.Second,a wide field electromagnetic test is conducted in the complex structure area of southern Sichuan.After data processing and inversion imaging,apparent resistivity logging data are used for calibration to develop an apparent resistivity interpretation model suitable for the test area.On the basis of the results,the characteristics of the electrical structure change in the shallow longitudinal formation of 6 km are implemented,and the transverse electrical distribution characteristics of the deep shale gas layer are delineated.In the prediction area near the well,the subsequent data verification shows that the apparent resistivity obtained using the inversion of the wide field electromagnetic method is consistent with the trend of apparent resistivity revealed by logging,which proves that this method can effectively identify the weak response characteristics of deep shale gas formations in complex structural areas.This experiment,it is shown shows that the wide field electromagnetic method with a large depth and high precision can effectively characterize the electrical characteristics of deep,low-resistivity thin layers in complex structural areas,and a new set of low-cost evaluation technologies for shale gas target layers based on the wide field electromagnetic method is explored.
文摘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.
文摘The signs of the electric field markers in Figs.2 and 4 of the paper[Chin.Phys.B 32104211(2023)]have been corrected.These modifications do not affect the results derived in the paper.
基金Supported by Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (2019BT02H594)Sanya Technology Innovation Special Project (2022KJCX08)。
文摘Offshore carbon dioxide(CO_(2)) storage is an effective method for reducing greenhouse gas emissions. However, when using traditional seismic wave methods to monitor the migration of sequestration CO_(2) plumes, the characteristics of wave velocity changes tend to become insignificant beyond a certain limit. In contrast, the controllable source electromagnetic method(CSEM) remains highly sensitive to resistivity changes. By simulating different CO_(2) plume migration conditions, we established the relevant models and calculated the corresponding electric field response characteristic curves, allowing us to analyze the CSEM's ability to monitor CO_(2) plumes. We considered potential scenarios for the migration and diffusion of offshore CO_(2) storage, including various burial depths, vertical extension diffusion, lateral extension diffusion,multiple combinations of lateral intervals, and electric field components. We also obtained differences in resistivity inversion imaging obtained by CSEM to evaluate its feasibility in monitoring and to analyze all the electric field(Ex, Ey, and Ez) response characteristics. CSEM has great potential in monitoring CO_(2) plume migration in offshore saltwater reservoirs due to its high sensitivity and accuracy. Furthermore, changes in electromagnetic field response reflect the transport status of CO_(2) plumes, providing an important basis for monitoring and evaluating CO_(2)transport behavior during storage processes.
基金supported by the National Natural Science Foundation of China(41922027,4214200052)by the Macao Foundation+1 种基金by the Pre-research Project on Civil Aerospace Technologies No.D020308/D020303 funded by China National Space Administrationby the Macao Science and Technology Development Fund,grant No.0001/2019/A1。
文摘Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellites,inverting the Earth’s three-dimensional conductivity distribution on a global scale becomes attainable.A key requirement in the global conductivity inversion is to have a forward solver with high-accuracy and efficiency.In this study,a finite volume method for global electromagnetic induction forward modeling is developed based on unstructured grids.Arbitrary polyhedral grids are supported in our algorithms to obtain high geometric adaptability.We employ a cell-centered collocated variable arrangement which allows convenient discretization for complex geometries and straightforward implementation of multigrid technique.To validate the method,we test our code with two synthetic models and compare our finite volume results with an analytical solution and a finite element numerical solution.Good agreements are observed between our solution and other results,indicating acceptable accuracy of the proposed method.
基金jointly supported by the project of Chinese National Natural Science Foundation(42107485)National Key R&D Program(2020YFC1512400,2018YFC800804)China Geological Survey(DD20190282,DD20221734,and DD20230323)。
文摘The construction of modern livable cities faces challenges in karst areas,including ground collapse and engineering problems.Wuhan,with a population of 13.74×10^(6) and approximately 1161 km^(2)of soluble rocks in the urban area of 8569.15 km^(2),predominantly consists of concealed karst areas where occasional ground collapse events occur,posing significant threats to underground engineering projects.To address these challenges,a comprehensive geological survey was conducted in Wuhan,focusing on major karstrelated issues.Geophysical methods offer advantages over drilling in detecting concealed karst areas due to their efficiency,non-destructiveness,and flexibility.This paper reviewed the karst geological characteristics in Wuhan and the geophysical exploration methods for karst,selected eight effective geophysical methods for field experimentation,evaluated their suitability,and proposed method combinations for different karst scenarios.The results show that different geophysical methods have varying applicability for karst detection in Wuhan,and combining multiple methods enhances detection effectiveness.The specific recommendations for method combinations provided in this study serve as a valuable reference for karst detection in Wuhan.
文摘In recent years,karst construction projects in the built-up area of Wuhan(capital of Hubei Province,China)are increasing,and the karst geological disasters have aroused social concerns.The actual engineering projects usually use shallow geophysical exploration methods to explore karst.This paper uses Spatial Auto-Correlation Method(SPAC)and electromagnetic Computerized Tomography(CT)to detect karst in urban built-up areas.Depending on the different physical properties of rock and soil,the SPAC method can better reveal the interface between soil and rock strata and the interface between soil layers.The electromagnetic CT method can identify strata according to the apparent absorption coefficient,which can better reveal the interface between soil and rock,the interface between the more intact and weathered rock.The SPAC method is mainly qualitative to measure the low-speed area,namely,the wrong geological body i.e.,karst cave,but also can detect the fracture zone or filling mode of karst cave,and at the same time,cannot use exploration holes or logging observation.The electromagnetic CT method can accurately detect the location and scale of the karst caves and has a higher accuracy detecting karst bands.In addition,exploration holes or well logging observations are also expected to be conducted,and their detection effect is greatly affected by lithology.
基金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 National Nature Science Foundation of China(Grant No.40874055)the Natural Science Foundation of the Hunan Province,China(Grant No.14JJ2012)
文摘Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method (EFGM), the point interpolation method (PIM), and the radial point interpolation method (RPIM). Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled- source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG-PIM and EFG-RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM.
文摘Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-frequency electromagnetic fields on the environment and which are capable of modifying and destroying our ecosystem. Particular interest was given in this article due to a massive influx of population living near high voltage lines. The analysis and simulation of the influence of low frequency electromagnetic fields on living beings in the vicinity of high voltage sources 132 kV and 220 kV in urban areas in DR Congo is the subject of our research with a view to estimating the level of exposure of humans to low frequency electromagnetic fields. To carry out our research, we used the classic method of analyzing the field produced near a high voltage line based on Maxwell’s image theory, the Maxwell-Gauss theorem and Maxwell-Ampère theorem to model and quantify low-frequency electromagnetic fields in the vicinity of high-voltage lines. The 2D FDTD numerical formulation was developed from telegraphers’ equations and allowed us to obtain models of current and voltage induced by electromagnetic fields on living beings below and near HV lines. The different simulations carried out on the proposed models illustrate the effects of the electrical and geometric parameters of the pylons on the distribution of the electromagnetic field in the vicinity of the HV lines. The results obtained were compared to the safety limits recommended by the standards.
文摘A class of finite step iterative methods, conjugate gradients, for the solution of an operator equation, is presented on this paper to solve electromagnetic scattering. The method of generalized equivalent circuit is used to model the problem and then deduce an electromagnetic equation based on the impedance operator. Four versions of the conjugate gradient method are presented and numerical results for an iris structure are given, to illustrate convergence properties of each version. Computational efficiency of these methods has been compared to the moment method.
文摘The quality of the low frequency electromagnetic data is affected by the spike and the trend noises.Failure in removal of the spikes and the trends reduces the credibility of data explanation.Based on the analyses of the causes and characteristics of these noises,this paper presents the results of a preset statistics stacking method(PSSM)and a piecewise linear fitting method(PLFM)in de-noising the spikes and trends,respectively.The magnitudes of the spikes are either higher or lower than the normal values,which leads to distortion of the useful signal.Comparisons have been performed in removing of the spikes among the average,the statistics and the PSSM methods,and the results indicate that only the PSSM can remove the spikes successfully.On the other hand,the spectrums of the linear and nonlinear trends mainly lie in the low frequency band and can change the calculated resistivity significantly.No influence of the trends is observed when the frequency is higher than a certain threshold value.The PLSM can remove effectively both the linear and nonlinear trends with errors around 1% in the power spectrum.The proposed methods present an effective way for de-noising the spike and the trend noises in the low frequency electromagnetic data,and establish a research basis for de-noising the low frequency noises.
基金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 Young Scientists Fund of the National Natural Science Foundation of China(No.41304082)the Young Scientists Fund of the Natural Science Foundation of Hebei Province(No.D2014403011)the Geological survey project of China Geological Survey(No.12120114090201)
文摘We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized: the real number of the electromagnetic field solution based on the Gaver-Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the "ground excitation-stratum measurement method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drill- hole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.
基金supported by the Doctoral Fund Project of the Ministry of Education(No.20130061110060 class tutors)the National Natural Science Foundation of China(No.41504083)National Basic Research Program of China(973Program)(No.2013CB429805)
文摘Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary(scattered) fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.
基金financially supported by the Thirteenth Five-Year-Plan Major Project "Marine Shale Gas Exploration and Evaluation over Laifengxianfeng and Hefeng Block"(No.2016ZX05034004-004)China Huadian Engineering Co.,LTD(No.CHEC-KJ-2014-Z10)
文摘In an effort to reduce the shale gas exploration risks and costs, we applied the wide-field electromagnetic method (WFEM), because of its strong anti-interference capability, high resolution, ability to conduct exploration at large depths, and high efficiency, to the Bayan Syncline in the South Huayuan block, Hunan Province. We collected rock samples and analyzed their resistivity and induced polarization (IP) and built A series of two-dimensional models for geological conditions to investigate the applicability of WFEM to different geological structures. We also analyzed the correlation between TOC of shale and the resistivity and IP ratio to determine the threshold for identifying target formations. We used WFEM to identify the underground structures and determine the distribution, depth, and thickness of the target strata. Resistivity, IP, and total organic carbon were used to evaluate the shale gas prospects and select favorable areas (sweet spots) for exploration and development. Subsequently, drilling in these areas proved the applicability of WFEM in shale gas exploration.
基金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.