Search

The borehole-ground transient electromagnetic method enhances the detection range and resolution by placing the transmitter electrode within the borehole, allowing for close-proximity excitation of subsurface targets and compensating for the shortcomings of traditional ground-based methods. Utilizing an unstructured vector finite element method combined with a second-order backward Euler variable time-stepping difference scheme and the MUMPS solver, we have achieved three-dimensional forward modeling simulation of the transient electromagnetic field for a borehole-ground electric source. Based on the validation of its accuracy, we analyzed the effectiveness of this method for dynamic monitoring of shale gas reservoir fracturing. Through forward modeling, we examined the characteristics of the radial electric field response during shale gas reservoir fracturing monitoring and evaluated the effectiveness of the borehole-ground TEM method in dynamic monitoring of shale gas reservoirs. The research results indicate that this method can meet the requirements for reservoir fracturing dynamic monitoring and has a promising application prospect.

A long wire with large current source transient electromagnetic (TEM) monitoring, with a large detection depth, low cost, safety, and environmental protection, has unique advantages in the testing and identification of unconventional reservoir fluid and the evaluation of stimulated reservoir volume. So, the TEM 3D forward modeling method has become a research hotspot. Although the finite-element method (FEM) is a type of numerical algorithm that has been widely applied in three-dimensional (3D) electromagnetic field forward modeling, the efficiency and accuracy of FEM require further improvement in order to meet the demand of fast 3D inversion. By increasing the order of the basis function and adjusting the principle of mesh discretization, the precision of the mixed-order spectral-element (SEM) result will be increased. The backward Euler scheme is an unconditionally stable technique which can ignore the impact of the scale of the time step. To achieve a better description of the nonlinear electromagnetic (EM) response of the grounded source TEM method and to optimize the efficiency and accuracy/precision of the 3D TEM forward modeling method significantly, we proposed the use of 3D TEM forward modeling based on the mixed-order SEM and the backward Euler scheme, which can obtain more accurate EM results with fewer degrees of freedom. To check its accuracy and efficiency, the 1D and 3D layered models are applied to compare the SEM results with the semi-analytical and FEM solutions. In addition, we analyzed the accuracy and efficiency of the SEM method for different types of order basis functions. Finally, we calculated the long-wire source TEM response for a practical 3D earth model of a shale gas reservoir for fracturing monitoring and tested the feasibility of the TEM method in a hydraulic fracturing monitoring area to further demonstrate the flexibility of the SEM method.

The Shanxi Graben System exhibits high seismic activity and significant destructive potential. Previous studies in this region have primarily focused on geological methods such as seismology, geological surveys, and trench excavations, with limited research in the field of geophysics. In this study, we collected two magnetotelluric profiles crossing the central segment of the Shanxi Graben System, including the fault system and the Hongdong–Linfen seismic zone. Through qualitative analysis and inversion calculations, we constructed a three-dimensional electrical model of the study area. The seismogenic environment was studied by integrating the deformation field, recent seismic geological data, and geophysical survey results of the study area and its surroundings. The Luoyunshan Piedmont Fault and the Huoshan Piedmont Fault are likely two main faults in the central segment of the Shanxi Graben System. These faults span the entire crustal scale and serve as basement faults separating the Ordos Block and the Taihangshan Block. They exhibit patterns of activity. The western side of the Shanxi Graben System, represented by the Ordos Block, exhibits a stable tectonic environment, while the eastern side, represented by the Taihangshan Block, experiences severe lithospheric destruction and thinning. The results of the magnetotelluric (MT) survey support the 1303 Hongdong earthquake as an event that occurred on the Huoshan Piedmont Fault and provide geophysical evidence for the possible dominance of the Luoyunshan Piedmont Fault in the 1695 Linfen earthquake. Multiple factors controlled the seismogenic environments of these two earthquakes. The continuous upwelling of asthenospheric material in the lower-middle crust on the eastern side of the Linfen Basin possibly leads to the regional extension of the North China Block. This, in turn, triggers inclined sliding along the Huoshan Piedmont Fault and the Luoyunshan Piedmont Fault, which may be the main controlling factors for major earthquakes in the region.

As a new green energy source, geothermal resource’s exploration, development, and utilization are an important direction in geophysical exploration at present. In this study, the actual land surface temperature was inferred based on the thermal infrared band of Landsat8 remote-sensing images, and the information about the surface anomalies and their spatial distribution was obtained through a multifactor analysis. In addition, three magnetotelluric sounding profiles were deployed in the study area, and the geo-electric sections in the study area were obtained through inversion of the measured data. Then, based on the inverse geo-electric information and the land surface temperature anomaly information, we analyzed and verified the geothermal resource genesis of the thermal anomaly area and inferred the favorable geothermal resource area in the study area. The results show that these two methods can be used to compare and analyze the possible distribution of the geothermal resources in the study area in two dimensions: the spatial distribution on the surface and the vertical distribution in the subsurface. Moreover, the results of the geothermal anomalies inferred from the thermal infrared remote sensing and the geo-electric results inferred from the magnetotelluric data are in good agreement. This study demonstrates that the integrated application of thermal infrared remote sensing and magnetotelluric technology is a promising tool for geothermal exploration.

The borehole induced polarization method has been widely used in deep mineral exploration, oil and gas resource exploration, and water resource exploration because of its high efficiency and good exploration effect. At present, the related research on the cross-well induced polarization method assumes that the underground medium is isotropic, but the electrical characteristics of the actual earth medium are anisotropic. To analyze the influence of the anisotropic characteristics on the cross-well induced polarization method, in this paper, the anisotropic forward algorithm of conductivity and polarizability in different principal axis directions based on the finite element method is studied. A three-dimensional forward simulation of the cross-well induced polarization method in anisotropic media is realized. The effectiveness and correctness of the algorithm are verified by testing and comparing complex 3-D isotropic and anisotropic models. Anisotropic geological models of the horizontal plate and inclined plate are constructed to analyze the anisotropic influences of conductivity and polarizability in different principal axis directions on the cross-well induced polarization response. The results show that the emitter sources with different depths in the well have different influences on the electrical response of the plates. Anisotropic conductivity and polarizability in horizontal plates exhibit most pronounced characteristics in the x-direction, significantly influencing the apparent polarizability curves. However, when the resistivity and polarizability are both anisotropic, the change in the z-direction is the most complicated. When the plate is inclined, the amplitude of the electrical response curve decreases to a certain extent, and the position where the amplitude appears shifts to different degrees. Notably, the response curves of the y-direction anisotropy are basically consistent with the response curves of the isotropy, regardless of the anisotropy of the conductivity and polarizability or anomalous body tilts. The results of this study improve our understanding of the influence of anisotropy on cross-well induced polarization and provide theoretical support for the interpretation of cross-well induced polarization data considering anisotropy. [ABSTRACT FROM AUTHOR]

Currently, the use of Cagniard apparent resistivity and phase is still the main method and means to process and interpret controlled-source audio-frequency magnetotelluric (CSAMT) sounding data. The CSAMT data must meet the conditions in the far zone to use the magnetotelluric (MT) sounding inversion interpretation method. The conditions in the far zone are directly related to transceiver distance; Hence, the distance measured by the CSAMT method is directly related to the reliability and credibility of the data. In this study, based on the formation wave and ground wave of the analytical expression for the electric field in a uniform half-space condition, the response laws of the formation wave and the ground wave are analyzed. After rigorous mathematical derivation, the formula for determining the optimal transceiver distance of the CSAMT is obtained, and the steps for using this formula in practical application are given. The actual field application verifies the rationality of the derived formula. This study provides a reference for the design of a CSAMT field source, and effectively guides the determination of CSAMT transceiver distance in the field. [ABSTRACT FROM AUTHOR]

In recent years, a variety of modern signal processing methods have been used to suppress the noise of Magnetotelluric (MT) signals, and have achieved certain effects on the improvement of MT data quality. However, we found that the traditional MT denoising method usually loses some useful low-frequency information in the original signal while suppressing the strong interference. Therefore, we try to apply CEEMDAN-DFA combined with endpoint detection technology to the noise suppression of MT data in order to improve the quality of MT data. Complete Ensemble Empirical Mode Decomposition with Adaptive Noise (CEEMDAN) is developed on the basis of empirical mode decomposition (EMD) and ensemble empirical mode decomposition (EEMD). CEEMDAN can well avoid the mode mixing. that often occurs in EMD, and it has the advantages of better decomposition completeness, smaller reconstruction error and lower calculation cost than EEMD. The scaling exponent α in the detrended fluctuation analysis (DFA) is used to select useful IMFs to reconstruct the signal and has a credible evaluation criterion. The results of simulation experiments show that the denoising effects obtained by using this method is obviously better than that of wavelet transform (WT) and EMDDFA. Then, in view of the fact that the whole section of data is processed in the conventional MT denoising methods, while suppressing the strong interference, it will cause signal sections without strong interference to be processed (overprocessing). In this paper, CEEMDAN-DFA is used to filter the noise signal and obtain the profile curve of the noise. Then, the starting and ending points of the interference is determined by using the endpoint detection technology, and only the signal fragments which are determined to contain strong interference in the noisy signal are filtered, and the undisturbed signal fragments are reserved. The results of the measured data after processing show that, compared with traditional methods, the proposed method can not only suppress strong interference better, but also retain more useful information in the original MT signal, and the obtained apparent resistivity-phase curve can more truly reflect the geoelectric information of the underground structure. [ABSTRACT FROM AUTHOR]

When transient electromagnetic investigation methods are carried out in the field, the measured data often contain both the induced polarization (IP) effect and the electromagnetic effect. In order to study the IP effect in the transient electromagnetic response, many researchers first calculate the electromagnetic field which considers the IP effect by replacing traditional resistivity with complex resistivity of the Cole-Cole model in the frequency domain. After the forward modeling calculation of the electromagnetic field in the frequency domain that considers the IP effect, the transient electromagnetic field in time-domain is obtained by a time-frequency transform algorithm. In this paper, the resistivity is directly replaced by the time-variant resistivity expression of the Cole-Cole model by using digital filter algorithms when we simulate the transient electromagnetic fields in time-domain. The calculated result of the Cole-Cole model in time-domain and in frequency-domain are consistent with each other, as observed in the horizontal electric field and the vertical magnetic field comparisons, which indicates the correctness of the numerical computation method adopted in this paper. The research presented herein allows us to observe the influence of the IP effect on transient electromagnetic field as well as study the mechanisms of IP directly. [ABSTRACT FROM AUTHOR]