Search

Detecting and mapping subsurface utilities in urban areas is crucial for identifying defects or damages in drinking and sewage pipes that can cause leaks. These leaks make it difficult to accurately characterize the pipes due to changes in their reflective properties. This study focused on detecting leaks originating from underground pipes and distinguishing between these various types of pipes. It also aimed to create a visual fingerprint model that displays the reflection characteristics of these pipes during different leak conditions, enabling efficient maintenance and handling procedures on the pipes. To achieve this, a finite-difference time-domain (FDTD) method was used to simulate two types of pipe materials with and without leak areas to construct different scenarios. Additionally, a ground-penetrating radar (GPR) field survey was conducted using a 600 MHz antenna in a part of the El Hammam area on Egypt's northwest coast. The simulated images produced with numerical modeling were compared with the radar profiles obtained using GPR at particular locations. The numerical simulations and radar profiles demonstrated the noticeable influence of water leaks from the different pipes, wherein the reflection of saturated soil waves was interrupted due to the presence of saturated soil. Envelope and migration techniques were employed in a new application to accurately distinguish between different pipe types, specifically focusing on leak areas. The strong correlation between the real radar profile and the specific signal of a water pipe leak in the simulated models suggests that GPR is a reliable non-destructive geophysical method for detecting water pipe leaks and distinguishing between the different pipe materials in various field conditions. The simulated models, which serve as image-matching fingerprints to identify and map water pipe leaks, help us to comprehend reality better. [ABSTRACT FROM AUTHOR]

Inversion of magnetotelluric (MT) responses has been used to explore the electrical conductivity distribution of the Earth's interior. In three-dimensional (3-D) inversions, it is significant to use a good initial model, because final model obtained by most 3-D inversion methods is influenced by the initial model. Although uniform initial models are widely used in 3-D inversions, one-dimensional (1-D) initial models are alternatives, which could more appropriately represent the actual conductivity variations in the Earth's interior. This study presents a two-step 3-D inversion method, especially for marine cases. This inversion method first concentrates on obtaining a 1-D initial model and then inverts for 3-D conductivity structures with it, in both of which the 3-D topography is carefully taken into consideration. This method was tested by synthetic models of different topography variations (depression-shaped, smoothly varying, channel-shaped and square-shaped plateau topography) and of heterogeneous layers with different checkerboard-type anomalies (sharp or smooth lateral conductivity variations) embedded in 1-D model of depth-dependent conductivity. The comparisons were done about obtaining an initial model by the proposed method and that inverted from the corrected responses. The results of 3-D inversions by using the method of this study were also compared to that with different uniform initial models. Results of synthetic tests and comparisons were discussed by using directional information of newly introduced model-vector parameters. The performance and validity of this method was verified. It also revealed that some of the newly introduced model-vector parameters could be used to show the convergence of inversions and help to select inverse model. [ABSTRACT FROM AUTHOR]

Qingchengzi is an important polymetallic mineral concentration zone in China. However, long-term mining has led to the depletion of proven reserves in the Qingchengzi ore concentration zone. Therefore, it is necessary to carry out prospecting in the periphery of the existing ore bodies. Taking into account the complex surface conditions of the Qingchengzi mining area, we carried out an airborne time-domain electromagnetic survey in this region and obtained a high-precision underground resistivity model which covers the whole area. The quasi-3D inversion result shows a good agreement with the surface geological information. After analyzing the relationship between the underground geological structures and the resistivity model, we constructed a geological model with the distribution of the main strata, the rock masses and the faults. The reconstructed resistivity distribution of the entire region added new knowledge to understand the geology of the Qingchengzi survey area. [ABSTRACT FROM AUTHOR]

Fully coupled lithosphere, atmosphere, and ionosphere theory has demonstrated that extremely low-frequency electromagnetic (ELF-EM) fields present a broad application prospect in deep resource exploration, but previous studies have ignored the contribution of the Earth's curvature. This study extends the theory of ELF-EM over a stratified Earth to the case where the Earth's curvature must be taken into account, and presents an analytical solution of the ELF-EM field excited by a grounded horizontal antenna in a spherical Earth–ionosphere model, whose theoretical approach and solution method are notably different from the flat Earth–ionosphere model. Additionally, the Earth is treated as a concentric-layered sphere rather than an ideal homogeneous sphere. We aim to investigate the effects of the Earth's curvature on the surface field, so as to broaden the coverage of the ELF wave in resource exploration. The solution is mathematically accurate and physically reasonable, since it reflects the sphericity and radially stratified structure of the Earth. We first verify the correctness and reliability of the proposed method by comparing the results with FDTD in a full-space spherical model. Additionally, we then compared the spherical results with the conventional controlled-source electromagnetic method and flat Earth–ionosphere results. The results show that when the distance between the transmitter and the receiver is comparable to the Earth radius, the spherical model better reflects the resonance of the wave in the cavity, suggesting that the effect of the Earth's curvature is not negligible. Then, the numerical simulations conducted to investigate the properties of the EM fields and their sensitivities to the conductivity at depth in the Earth are discussed. Finally, the EM responses of some simple electrical conductivity structures models are modeled to illustrate their prospects in future resource exploration. [ABSTRACT FROM AUTHOR]

Skywave refers to the electromagnetic wave reflected or refracted from the ionosphere and propagated in the form of a guided wave between the ionosphere and the Earth's surface. Since the skywave can propagate over large distances, it has been widely used in long-distance communication. This paper explores and demonstrates the feasibility of skywave for deep resource and energy exploration at depths of up to 10 km. Theoretical and technical advancements were accomplished in furthering skywave applications. A new solution method based on Green's function has been developed to study skywave propagation in a fully coupled lithosphere-air-ionosphere full space model. For the first time, the model allows one to study skywave distribution characteristics in the lithosphere containing inhomogeneity such as ore deposits or oil and gas reservoirs. This model also lays a foundation for skywave data processing and interpretation. On a parallel line, we have developed a multi-channel, broadband, low-noise, portable data acquisition system suitable for receiving skywave signals. Using the skywave field excited by a high-power fixed source located in central China, actual field surveys have been carried out in some areas in China including the Biyang depression of Henan Province. The initial results appear encouraging—the interpreted resistivity models prove to be consistent with those of seismic exploration and known geological information, and the exploration cost is only ∼1/4 to 1/10 that of seismic surveys. These initial successful applications of the skywave theory lay a solid foundation for further verification of the new method. [ABSTRACT FROM AUTHOR]

A summation technique has been developed based on the continuous fractional expansion to accelerate the convergence of infinite series involving the product of Riccati-Bessel functions, which are common to electromagnetic applications. The series is transformed into a new and faster convergent sequence with a continued fraction form, and then the continued fraction approximation is used to accelerate the calculation. The well-known addition theorem formula for spherical wave function is used to verify the correctness of the algorithm. Then, some fundamental aspects of the practical application of continuous fractional expansion for Mie scattering theory and electromagnetic exploration are considered. The results of different models show that this new technique can be applied reliably, especially in the electromagnetic field excited by the vertical electric dipole (VED) source in the "earth-ionospheric" cavity. The comparison among the new technology, theWatsontransform, and the spherical harmonic series summation algorithm shows that this new technology only needs less than 120 series items which is already enough to obtain a small relative error, which greatly improves the convergence speed, and provides a new way to solve the problem. [ABSTRACT FROM AUTHOR]

The article discusses research which assessed the behavior of downhole electromagnetic telemetry (EMT) in deep subsurface prospecting and the characteristics of uphole EMT with repeater system. Topics covered include the performance of downhole EMT during sudden water influx into adjacent formation during directional drilling, the characteristics of downhole EMT in deep horizontal well, and the behaviors of the EMT system with oil based mud (OBM) and water based mud (WBM).

Geophysics offers an important means to investigate the physical processes occurring inside the earth. In particular, since the 1960s, electromagnetic (EM) methods have played important roles in mineral exploration and engineering investigation. Such investigation requires extensive data acquisition and experimental analysis based on geophysical techniques. However, high-performance geophysical equipment, particularly EM exploration equipment, has been dominated by large geophysical companies from the United States, Canada, Germany, and other European countries for decades. This has limited the development of deep exploration technology in China. Recently, we have developed a high-resolution acquisition system with a wireless control unit and a high-power transmitting system for surface EM prospecting (SEP). The new system has been tested in the high-intensity, noisy environment in Jian-sanjiang area, Heilongjiang Province. We then conducted a field survey on the western edge of Qaidam Basin, Qinghai Province. A highly conductive anomaly was found in the upper mantle below the Qinmantage Mountains, which indicates a possible northward fluid channel that runs from below the Qinmantage Mountains to the bottom of the western crust of the Qaidam Basin. Identification of this significant feature was made possible by the new SEP for its better resolution than the previous systems. Also, geophysical analysis confirmed that the thick Cenozoic sediments of the Qaidam Basin are underlain by rigid Precambrian basement rocks and are characterized by a series of folds. The resistivity profile indicates that the Qaidam Basin was formed due to the folding structures in the northern part of the Qinghai-Tibet Plateau which provided an additional evidence for the uplifting of the Qinghai-Tibet Plateau. [ABSTRACT FROM AUTHOR]

The North China Craton is the oldest continental block, and has suffered from large-scale lithospheric thinning and destruction, which in turn led to gold deposits in northern China. The decratonic gold deposits in the North China Craton became the most important gold deposits in China, and geophysical methods are key means to detect and discover gold deposits there. In this paper, based on the geological and petrophysical characteristics of the North China Craton, the geological model of the decratonic gold deposits is transformed into a geophysical model. At present, two methods of geophysical exploration of decratonic gold deposits are in use: rapid and efficient exploration on the scale of the ore concentration area, and large depth exploration on the scale of the deposit area. In detail, the airborne electromagnetic, magnetic and gravity methods are used to detect the shallow (1,500 m) anomaly area on the scale of the ore concentration area. Through the ground-controlled source electromagnetic and ground magnetotelluric methods, explorations for targets at significant depth (5,000 m) are carried out in the mining area. Then, taking the Liaodong ore concentration area as an example, geophysical methods are used to discover two prospecting areas around the Jianshanzi Fault in the Qingchengzi ore concentration area, Baiyun-Xiaotongjiapuzi deep prospecting area, and Qingchengzi deep prospecting area. Next, three prospecting areas are delineated around the Jixingou Fault in the Wulong mining area, Wulong deep prospecting area, Weishagou deep prospecting area, and Chang'an deep prospecting area. The anomalies in the ore concentration area and mining area are revealed by means of three-dimensional exploration methods, thereby providing technical support for the exploration of metal minerals such as decratonic gold deposits. [ABSTRACT FROM AUTHOR]

Transient electromagnetic method (TEM) is an important tool for mapping the subsurface resistivity. As the electromagnetic fields with different frequencies reach different depth of the Earth, the resolution of TEM method is closely related to the spectrum properties of its source signal. The pseudorandom binary sequence (PRBS) is a recently applied source signal and is thought to be superior to the conventional square wave for its adjustable frequency range and density. However, to invert the response from this optimized source, the pure secondary field needs to be extracted from the full waveform response in which parts of useful information is lost inevitably. We thus propose to directly invert the full waveform TEM response from the excitation of PRBS using a transdimensional framework. We use the reverse-jump Markov chain Monte Carlo (RJ-MCMC) sampler to get the model ensemble so that the posterior of the model can be obtained. The PRBS source is convolved with the impulse response of the Earth to fit the recorded full waveform response. The 2-D model space is discretized using a two-stage binary tree and updates in parallel, which dramatically improve the sampling efficiency. The synthetic examples compare the results from the conventional time-domain step response and the full waveform response and demonstrate that the proposed method obtains higher resolution and can better resolve the interface in the shallow and moderate depth. Furthermore, the real data example recovers the model structure that is consistent with the geologic inference, providing alternate evidence for the conjecture of the Sanshandao fault. [ABSTRACT FROM AUTHOR]

Taohemu deposit is a superlarge silver‐polymetallic deposit discovered in the Great Xing'an Range metallogenic belt of China in recent years. The deposit occurs in the terrestrial volcanic rock sequences of Late Jurassic age. The ore bodies consist mainly of rhyolitic tuffs, brecciated tuffs, and volcanic breccia bearing Ag‐Pb‐Zn mineralization, and are largely conformable with the surrounding layered volcanic and volcaniclastic rocks with a gradual transitional relationship manifesting as strata‐bounded features in general. The mineralization system appears to be generally controlled by a paleo‐volcanic chamber. It is suggested that the late volcanic hydrothermal fluids flowed within the permeable layers of rhyolitic pyroclastics consisting mainly of volcanic breccias and breccia‐bearing tuffs. This resulted in the formation of superlarge Ag‐Pb‐Zn polymetallic ore deposit metasomatism and fracture filling. The majority of the deposit occurs as shallowly dipping concealed ore bodies, as defined by sparsely distributed drill holes. Because of its short discovery history, new mineralization style and low degree of comprehensive research, the boundaries of most ore bodies have not been effectively delineated along vertical and lateral extensions. At present, it is crucial to strengthen the comprehensive studies on the deposit, conduct more detailed geological exploration, and establish the metallogenic and prospecting models. For this reason, we have designed and implemented long sections of magnetotelluric (MT) and controlled‐source audio magnetotelluric (CSAMT) sounding across the Taohemu ore field and obtained a large number of geological structures and valuable underground exploration results that provided reliable geophysical support for promoting deep resource exploration. Plain Language Summary: The geological structures associated with metallogenic system is usually covered by surface vegetation and Quaternary layer. Furthermore, surface mineralization is not relevant to deep ore bodies. To reveal the shape of ore bodies and geoelectrical structures in this study area, we applied geophysical electromagnetic methods to achieve this purpose. The results showed that the ore bodies are largely conformable with the surrounding rocks as strata‐bounded features in general tested by borehole information. The mineralization system appears to be generally controlled by a paleo‐volcanic chamber. Key Points: Underground geoelectric structures are imaged from MT and CSAMT exploration over Taohemu superlarge silver‐polymetallic deposit, Inner MongoliaOre bodies are conformable with surrounding volcanic and volcaniclastic rock with a gradual transitional relationship manifesting as strata‐bounded featuresThe mineralization system appears to be generally controlled by a paleo‐volcanic chamber [ABSTRACT FROM AUTHOR]

As an important branch of geophysical exploration method, the electromagnetic method with artificial source has advanced rapidly in the past decade. These methods are classified as airborne electromagnetic method, ground-air electromagnetic method, ground electromagnetic method, and marine electromagnetic method. Over the years, researchers in China have made significant improvement to the fundamental theory, forward modeling and inverse for series of electromagnetic detection methods. Conversely, significant progress was made in the development of corresponding equipment. The researched techniques and their developed equipment have been successfully utilized to detect underground targets as deep as 10 km. However, there is increasing need for deep resources exploration, urban subsurface study, and prediction, monitoring and detection of geological hazards. To meet the increasing need and catch up with the advanced international level of exploration technologies and developed equipment, there is urgent necessity and requirement to continue developing geophysical methods and the corresponding equipment. [ABSTRACT FROM AUTHOR]

The electromagnetic fields generated by an artificial source can be divided into: 1) the stratum wave propagating in the earth; 2) the ground wave propagating along the earth surface; 3) and the sky wave traveled in the ionosphere-earth wave guide. Traditionally, the sky wave is mainly used for communication, which receives the electromagnetic information from upper half of the space. The sky wave in this paper is used for detection and receives the electromagnetic information from the whole space. Compared with the traditional sky wave, the sky wave information we study is the whole space propagation information. [ABSTRACT FROM AUTHOR]

Previously, the multi‐channel transient electromagnetic (MTEM) system has been successfully developed and applied to hydrocarbon reservoirs exploration on land and offshore environments in France. In order to explore the mineral resources buried at the depths up to 2,000 m in China, we developed a new MTEM system. In this paper, we present the components of the system, methods, principles, the parameters of the system, and the data interpretation method. To validate the effectiveness of the methodology and the developed system, we carried out an MTEM survey to locate the controlling faults of a gold deposit in the east of China. The survey contains 30 receivers distributed over 12‐km survey. By comparing the MTEM interpretation results with the known geological background, we find that MTEM methods can provide a good estimation of the ore‐controlling faults and this method could be a potential tool for the future deep mineral exploration. [ABSTRACT FROM AUTHOR]

Mantle tectonics, such as asthenospheric upwelling and lithospheric decoupling, usually controls subsequent fluid formation and migration from deep to shallow levels in the lower lithosphere and promotes rock failure and deformation. Based on the concept whereby the fluid migration location defines the structural inheritance, we aim to relate crustal processes to mantle tectonics using fluids revealed by a new magnetotelluric (MT) array across the Ailao Shan-Red River belt. As imaged via a 3D resistivity model, two lower-lithospheric conductive anomalies are determined as consistent with the Dian-Qiong (DQ) suture and Song Da (SD) belt and interpreted to contain interconnected melt. Because potassic magmatism was enabled by partial melting in the lower lithosphere, we could infer that the DQ and SD have reworked as the major fluid migration channels. Fluid migration is considered to drive lithospheric decoupling at a low-viscosity conductive layer, which is inferred feeding by aqueous and melt fluids originating from the two channels and diffusing at depths from 15 to 20 km. Constrained by the geochronology results, this weak layer could have provided convenience to induce the entirely upper crustal translation-rotation within the time interval between potassic magmatism and strike-slip of the Red River fault. This translation-rotation process is inherited from the underlying mantle processes and may further be remotely affected by the upper crustal movement of the Tibetan Plateau, conforming with the crustal rotation observed to the north of our study region. This work provides a compelling example of the tectonic control of the mantle on inherited responses in the crust. Plain Language Summary Crustal deformation is usually controlled by mantle tectonic processes, especially the southeastern Tibetan Plateau. Some models have proposed that lithospheric shearing and lower-crustal flow have controlled the crustal deformation in the Cenozoic. However, crustal rotation models show less possibility of a channel flow and relate the crustal deformation to the remote effect of India-Asia collision. Our study investigates above hypotheses by imaging the electrical resistivity of rocks beneath the major Red River Fault area in southeastern Tibet, using magnetotelluric data with high quality in a dense array. The lithosphere of the study area is interpreted to have been divided into two horizontal systems by a fluid diffusion layer at the bottom of upper crust. Resistivity low in the lower lithosphere is inferred due to partial melting of deep lithosphere and reworking of the paleo sutures in the Late Cenozoic, which have spread at the bottom of upper crust. Rather than the channel flow, this fluid migration process sensitively reflected in our model relates the mantle tectonics to crustal rotation by providing rheological conditions. We, hence, propose an inherited structure model featured by lithospheric decoupling and upper-crustal translation-rotation that may occur before the major strike-slip event. [Display omitted] • A new 3D model reveals electrical structures on the southeastern Tibetan Plateau. • The Red River fault is discontinuous at a depth of ∼20 km and may move away from the palaeo-suture. • Palaeo-suture was reconstructed into two activated zones and reworked as channels for fluid migration. • Fluid migration has driven the lithospheric decoupling and upper crustal inheritance. • The upper crust experienced a translation-rotation process before the internal strike-slip event. [ABSTRACT FROM AUTHOR]

Three emerging techniques for surface electromagnetic (EM) exploration with high detection accuracy and large depth of investigation, namely, WEM (wireless electromagnetic method), MTEM (multi‐channel transient electromagnetic method) and SOTEM (short‐offset grounded‐wire TEM), are introduced. In detail, we proposed a newly developed sky‐wave theory in which EM wave can propagate in full space of ionosphere‐atmosphere‐lithosphere. Correspondingly, we developed the WEM recording system to measure the sky‐wave to obtain the structures down to ∼10 km depth. We then introduced MTEM as a highly effective method for mineral exploration with target depth of 3km. Using self‐developed MTEM system, we presented the results obtained from mineral exploration study carried out in China. Lastly, we proposed the SOTEM method for a greater detection depth (1.5km) than traditional loop source TEM (500m) with a small offset and high resolution. [ABSTRACT FROM AUTHOR]

Electromagnetic telemetry is bi‐directional transmission of borehole data for effective drilling, especially in unconventional wells. In this study, numerical simulations was used to investigate the possibility of applying EM telemetry for measurement while drilling (MWD) in an ultra‐deep borehole in the Tahe region of the Tarim basin, using the drill string as an antenna excited as a gaped dipole source to quantify spatial variations of the electromagnetic field from subsurface to the surface. Results confirm that the distribution of the electric field along the drill string is sensitive to vertical variations of electric conductivity in the host rock. The simulated results show a high attenuation within the conductive formations, thereby limiting the effectiveness of EM telemetry in data transfer to the surface. In addition, the changes in thickness and resistivity of the target region affects the amplitude of the desired response. [ABSTRACT FROM AUTHOR]

Abstract The irregular structures of paleokarst systems are among the major factors resulting in carbonate reservoir heterogeneity, which restricts the further exploration of paleokarst reservoirs. Using multiscale geological-geophysical (G&G) data, including thin-section, core-sample and well-log data, the caves, cave fillings and fractures in individual wells in the northern Tahe Oilfield were characterized. Guided by seismic datasets, the interwell paleokarst structures were interpreted in profile and plane view. The paleokarst system in the study area was divided into four vertical zones, and the relative genetic structures were classified in the upper three zones. The epikarst zone includes paleo-canyons, karst towers, dolines, paleo-soils, and fractures, and the water flows horizontally. The vadose zone includes vadose caves and fractures, and the water flows vertically. The epiphreatic zone includes chambers, main channels, branch channels and fractures, and the water flows primarily horizontally along the water table. For each genetic type, the total hydrocarbon production and monthly production were analyzed. Furthermore, a genetic classification model was built for the paleokarst reservoirs of the study area, which increased our understanding of the structure of the paleokarst system. The results of the structure classification and production analysis can effectively guide hydrocarbon exploration, and this approach can be applied to similar paleokarst reservoirs. Highlights • Genetic model and related types of paleokarst system in Northern Tahe oilfield. • Multiscale geological-geophysical characterizations of paleokarst genetic structure types. • Hydrocarbon production characteristics of each genetic type in paleokarst system. [ABSTRACT FROM AUTHOR]

With fast development of Chinese economy, a large amount of garbage have been produced. Especially, high level radioactive waste (HLRW) bring dangers to the security of human being. Finding a suitable site to dispose HLRW is a key management to government. A potential HLRW disposal site has been chosen in west north of China at first, however, it is necessary to understand whether there exist electrical structures inside the rock mass, whether the rock mass is intact or not. We applied the controlled source audio-frequency magnetotelluric method (CSAMT) to confirm that. Three-dimensional numerical modeling of finite difference method was performed for typical structure in the study area, and response characteristics were analyzed to confirm which actual frequency range can be chosen and what structure dimension can be detected. Then we laid out 14 parallel survey lines to carry out data acquisition. By analyzing typical sounding curves, combined with drilling information and the inversion results, we suggest that the rock is intact and suitable for building a repository, which is a good recommendation for government department. [ABSTRACT FROM AUTHOR]

This paper introduces new technique for qualitative analysis of multi-transient electromagnetic (MTEM) earth impulse response over complex geological structures. Instantaneous phase and frequency attributes were used in place of the conventional common offset section for improved qualitative interpretation of MTEM data by obtaining more detailed information from the earth impulse response. The instantaneous attributes were used to describe the lateral variation in subsurface resistivity and the visible geological structure with respect to given offsets. Instantaneous phase attribute was obtained by converting the impulse response into a complex form using the Hilbert transform. Conversely, the polynomial phase difference (PPD) estimator was favored over the center finite difference (CFD) approximation method in calculating the instantaneous frequency attribute because it is computationally efficient and has the ability to give a smooth variation of the instantaneous frequency over a common offset section. The observed results from the instantaneous attributes were in good agreement with both the subsurface model used and the apparent resistivity section obtained from the MTEM earth impulse response. Hence, this study confirms the capability of both instantaneous phase and frequency attributes as highly effective tools for MTEM qualitative analysis. [ABSTRACT FROM AUTHOR]

Semi-airborne transient electromagnetic method (SATEM) has been widely used in geophysical investigations. This method uses a grounded-wire to transmit the source signal and installs the receivers on an aircraft so that the data with a large depth of investigation can be acquired efficiently. However, the instability of aircraft when collecting the data will introduce the non-geological factors that distort the recorded data sets, which further affect the precision of the investigation result greatly. The researching of the mechanism and the correction method of these factors, including geometric effect, geomagnetic effect and mutual inductance effect between transmitter and receivers, is of great significance to the SATEM method. Thus, we proposed a thorough scheme to suppress noise generated by the non-geological factors of the SATEM configuration. We first derived the formulas for calculating the mutual inductance coupling between transmitter and receiver. The formulas of the geometric noise and geomagnetic noise that are affected by the roll, pitch, and yaw of the aircraft are further derived to analysis the influence interval of the motion noise. Then a hybrid correction scheme that uses a two-stage strategy to correct the receiver noise is put forward. The results show that the geomagnetic field and the mutual conductance contribute more noises than the geometric effect. The geomagnetic field will distort the data at all time channel, while the effect of the mutual inductance mainly focused on the early time channels. A field example demonstrated that the proposed method could well suppress the non-geological noise and improving the signal to noise ratio. • The semi-airborne transient electromagnetic method (TEM) has been widely used. • The researching of the mechanism and the correction method of these factors is of great significance to the SATEM method. • The results are verified by field measurements. [ABSTRACT FROM AUTHOR]

The Alxa area in Inner Mongolia has been selected as a possible site for geological disposal of high-level radioactive waste (HLRW). Based on results of a previous study on crustal stability, the Tamusu rock mass has been chosen as the target. To determine the geological structure of this rock mass, aeromagnetic and gravity data are collected and inverted. Three-dimensional (3D) inversion horizontal slices show that the internal density of the rock mass and the distribution of magnetic properties are not uniform, with fractures and fragmentation being present. To confirm this result, the controlled source audio-frequency magnetotelluric method (CSAMT) was applied to explore the geological structures, the typical CSAMT sounding curve was analyzed, and the response characteristics of the geological structure and surrounding rock are distinguished. The original data were processed and interpreted in combination with data from surface geology and drilling and logging data. It is found that the CSAMT results were consistent with those from 3D inversion of the gravity and magnetic data, confirming the existence of fractures and fragmentation in the exploration area. [ABSTRACT FROM AUTHOR]

Dense array Magnetotelluric sounding (MT) was applied in the study of the emergence of massive gold deposits relating to the tectonic activity responsible for the creation of the North China Craton. This study focuses on understanding the shallow crustal electrical structure characteristics (>1.5 km) of Pb Zn and Au Ag deposits in the Qingchengzi orefield, through the analysis of shallow electrical resistivity distribution at depths less than 5 km, obtained by 3D inversion of full impedance at a frequency range of 320–0.107 Hz. The results showed that the electrical characteristics of Xiaotongjiapuzi and Taoyuan deposits in the 3D resistivity model are consistent with the known 2D resistivity model, and the Jianshanzi fault plays an important role in the gold mineralization, acting as a migration channel and precipitation zone for ore-forming materials. It is inferred that the fault extends at least 5 km deep with deep mineralization prospects. The strata of Liaohe Group in the north of the orefield, the intrusive granite in the south, and the NW-trending Jianshanzi and Yushanggou faults show distinguishable electrical properties, which are either high resistance, low resistance, or high and low resistance transition zone on the depth slice map. Furthermore, the Xinling rock mass has a NE trending extension, which confirms the inference that the Jianshanzi fault extends to a depth of 5 km into the crust. • Imaging shallow crustal electrical structures through 3D inversion of MT data. • Relationship of mineralization with the electrical structure, fault, and lithology. • Guide the deep prospecting direction of polymetallic deposits in the future. [ABSTRACT FROM AUTHOR]

To identify problematic geological structures which would be encountered when driving tunnels in the high mountainous area of South West China, a joint geophysical and engineering geological study was made. The tunnels will form part of a project to divert water from the Yangtze to the Yellow Rivers. A controlled source audio-frequency magnetotelluric method (CSAMT) was carried out in 2004 between the Ma-ke and Jia-qu rivers in high, steep terrain. The paper discusses the method of data collection, processing, and analysis for the West route. The faults/fractures identified will need to be taken into account in the design and construction of the tunnel. [ABSTRACT FROM AUTHOR]

The paper derives the finite element equation for acoustic wave in time domain and presents a transparent-plus-attenuation boundary condition. Forward modeling demonstrates that the boundary condition absorbs boundary reflection wave very well. On these bases, we derive the equation satisfied by elements of Jacobi matrix used in the inversion of the physical property parameters of acoustic media. In fact, the equation is the same as that of forward modeling in form. Only the right force item is different. So with the same method of forward modeling, we can get the elements of Jacobi matrix. Because the elements are variable with time and the present inversion does not permit too many unknowns. We integrate the finite elements with the same physical property as one unknown structure unit (for example, a horizontal layer or an oblique layer, etc.) and inverse the physical property parameters of these unknown structure units instead all element's unknown parameters. The method greatly reduces calculation time and saves computer memory. Also, it improves the accuracy of the inversion results and improves the stability of the solving process. The inversion equations are solved with QR decomposition method. Model results prove that the full wave equation inversion method in time domain is effective. [ABSTRACT FROM AUTHOR]

We conducted broadband Magnetotelluric survey to obtain 2D lithospheric electrical structures of Alxa block through two NE-SE trending profiles, which revealed the electrical evidence of bidirectional subduction of Paleo-Asian Ocean with spatial difference beneath Enger Us fault belt. [Display omitted] • Lithospheric electrical structures of the Alxa block were imaged through two NW-SE BMT profiles. • There was bidirectional subduction of the Paleo Asian Ocean beneath the Enger Us fault belt. • The low resistivity bodies of the middle and upper crust in the northwest are related to large-scale horizontal nappe tectonics. • The high conductivity bodies that extend deep in the southeast are attributed to the partial melting. The Alxa Block, situated between the Tarim and North China cratons, is essential to understand the tectonics of northwest China. It has undergone several tectonic cycles, including rifting from northern Gondwana in Neoproterozoic, subduction-accretion of the Paleo-Asian Ocean during Paleozoic, and intracontinental deformation in Mesozoic. Broadband magnetotelluric data covering a frequency range of 320 Hz to ∼ 10000 s was collected along two NW-SE trending lines across major tectonic units of the Alxa region to investigate the deep structures of the Alxa Block. By applying the nonlinear conjugate gradient (NLCG) inversion approach, two-dimensional (2D) electrical resistivity models were obtained at a 100 km depth. These findings indicate that the Zhusileng-Hangwula and Zongnaishan-Shalazhashan tectonic zones have high resistivity in their lower crust and upper mantle, with an overall layered high-low–high resistivity pattern. The upper crust of the Nuru-Langshan tectonic zone is primarily characterized by high resistivity values, whereas the lower crust and upper mantle contain a large-scale low-resistivity (or conductive) zone that is likely the result of partial melting. The large-scale Early–Middle Jurassic nappe structures may be connected to sub-horizontal low-resistivity anomalies in the shallow strata, which are located up to 20 km below the Yagan and Zhusileng–Hangwula tectonic zones. Below the Engger Us fault, two low-resistivity bands are interpreted as the product of bidirectional subduction and final closure of the Paleo-Asian Ocean. The Alxa Block may have subducted southwards beneath the Ordos Block, as suggested by a clear low-resistivity anomaly between the two blocks. Furthermore, partial melting may have happened for high-conductivity bodies in the upper mantle based on the modified Archie formula's evaluation of their melting degree. [ABSTRACT FROM AUTHOR]

This paper deals with the study of multi-channel adaptive noise cancellation with a focus on its application in electromagnetic (EM) telemetry. We presented new variable step-size least mean square (LMS) techniques: regularized variable step-size LMS and regularized sigmoid variable size LMS, for electromagnetic telemetry data processing. Considering the complexity and spatial distribution of environmental noise, algorithms with multiple reference signals were used to retrieve transmitted EM signals. The feasibility of the regularized variable step size LMS algorithms with numerical simulation was analyzed and presented. The adaptive processing techniques were applied in the recovery of frequency and binary phase shift key modulated signal. The proposed multi-channel adaptive technique achieves fast convergence speed, low mean squared error and is shown to have good convergence characteristics compared to conventional methods. In addition to attaining good results from the multi-channel adaptive filter and performing the signal analysis in real-time, we implemented combined fast effective impulse noise removal techniques. The combination of median and mean filters was effective in removing a wide range of impulsive noises without distorting any other data points. Further, electromagnetic telemetry data were acquired during a drilling operation in Sichuan province, China, for real field application. Data processing workflow was designed for EM telemetry data processing based on the noise characteristics, simulation results and expected result for demodulation. To establish a comprehensive overview, a performance comparison of the acquisition array system is also provided. Conclusively, the introduced multichannel adaptive noise canceling techniques are very effective in recovering transmitted EM telemetry signals. [ABSTRACT FROM AUTHOR]

This paper deals with the study of multi-channel adaptive noise cancellation with a focus on its application in electromagnetic (EM) telemetry. We presented new variable step-size least mean square (LMS) techniques: regularized variable step-size LMS and regularized sigmoid variable size LMS, for electromagnetic telemetry data processing. Considering the complexity and spatial distribution of environmental noise, algorithms with multiple reference signals were used to retrieve transmitted EM signals. The feasibility of the regularized variable step size LMS algorithms with numerical simulation was analyzed and presented. The adaptive processing techniques were applied in the recovery of frequency and binary phase shift key modulated signal. The proposed multi-channel adaptive technique achieves fast convergence speed, low mean squared error and is shown to have good convergence characteristics compared to conventional methods. In addition to attaining good results from the multi-channel adaptive filter and performing the signal analysis in real-time, we implemented combined fast effective impulse noise removal techniques. The combination of median and mean filters was effective in removing a wide range of impulsive noises without distorting any other data points. Further, electromagnetic telemetry data were acquired during a drilling operation in Sichuan province, China, for real field application. Data processing workflow was designed for EM telemetry data processing based on the noise characteristics, simulation results and expected result for demodulation. To establish a comprehensive overview, a performance comparison of the acquisition array system is also provided. Conclusively, the introduced multichannel adaptive noise canceling techniques are very effective in recovering transmitted EM telemetry signals. [ABSTRACT FROM AUTHOR]

Equivalence causes non-uniqueness in interpretations of the grounded-wire transient electromagnetic method (TEM). Many different layered models may produce the same response curves practically. However, the conventional theory of equivalence is based on the three-layer geo-electric models while the earth may not be reducible to a simple three-layer structure. Therefore, it is important to perform research on equivalence based on more complex models to reduce the uncertainty the grounded-wire TEM interpretations. Here, we extend the theory of equivalence based on a multi-layer target in the grounded-wire TEM, consisting of SS-, HH-, and TT-type. The analysis of relative error is applied to investigate the effect of measured components and receiver locations on equivalence. The results show that equivalence of multi-model occurs when any of three conditions is satisfied: the longitudinal conductance is constant in the far area, the thickness of targets do not change, or the transverse conductance of a multilayer target remains unchanged. In order to largely resolve the effect of equivalence on data interpretation, a strategy combined with short-offset survey and joint inversion of the electric component and magnetic component is presented and verified by the field survey in Shanxi Province, China. • Equivalence causes non-uniqueness in interpretations of the grounded-wire TEM method. • The conventional theory of equivalence is based on the three-layer geo-electric models. • It is important to perform research on equivalence based on more complex models. [ABSTRACT FROM AUTHOR]

An active clamp, three-phase, six-switch power factor correction rectifier with a one-cycle controller is proposed, which can effectively suppress the reverse recovery of the reverse parallel diode of the bridge arm switch and reduce the reverse recovery loss. The main switches and auxiliary switch are both zero-voltage switches. It works in a fixed switching frequency and low power stress during the switching period. The specific working process and control circuit are analyzed in detail, and the model simulation is carried out. The experimental platform of a 2.5 KW prototype, with a complete test and verification of the soft switch technology proposed in this paper, was set up in the laboratory. [ABSTRACT FROM AUTHOR]

• The gold deposits of the Qingchengzi orefield are controlled by faults. • The Jianshanzi fault zone is the major ore-controlling structure. • A magnetotelluric study was undertaken along the The Jianshanzi fault suggesting excellent potential for new, deep-seated gold discoveries in the Qingchengzi orefield. The Liaodong Peninsula in northern China has undergone multi-stage magmatic and tectonic modifications, similar to those recorded by the Jiaodong gold province that is the largest gold-producing district in China, indicating that Liaodong Peninsula has great potential for new discoveries. However, no significant gold deposit has been discovered thus far at depths below 1.5 km, the typical depth of the gold-bearing system in the comparable Jiaodong gold province. This study focused on the Qingchengzi orefield on the northern Liaodong Peninsula. Qingchengzi recorded multiple post-magmatic hydrothermal gold-mineralizing events along the Jianshanzi fault zone, the main ore-controlling structure in the orefield. Given the potential for undiscovered, deep-seated gold mineralization along the Jianshanzi fault zone, we performed a magnetotelluric (MT) study designed to identify any deeper exploration targets. The surface electromagnetic prospecting (SEP) system, a portable instrument suitable for use in electromagnetically mountain environments, was employed for the purpose of acquiring the MT data. Shallow conductive anomalies identified in the vicinity of the Jianshanzi fault zone were tested by known drilling recorder. The survey not only helped to generate a detailed map of the deep geological architecture, but also allowed us to estimate the vertical extent of the known gold orebody to a depth of more than 2 km. [ABSTRACT FROM AUTHOR]

In traditional ground penetrating radar (GPR) field data processing, e.g. horizontal stacking and migration, the wave propagating phase velocity is commonly used, while field magnitude decay is neglected. We derived the scalar wave equation of an electric vector field or magnetic vector field component and proved that the solution of this equation is the product of a decaying exponential function factor and a solution factor, which satisfies the pure wave equation. The waveform decomposition allows the evaluation of the absorption coefficient using the amplitudes of GPR waves. A meaningful method for absorption coefficient estimation in the time domain was introduced and tests with synthesized data demonstrated that this method is simple and effective. This method can simplify the process of absorption estimation and has a practical application value. Therefore, it is not only important for GPR data processing but may also yield important information on ground electrical properties. • A new absorption coefficient β is defined. • A simple method for β estimation is introduced. • Stack and migration processing will be efficiently implemented after β is obtained. [ABSTRACT FROM AUTHOR]

The Shijiazhuang-Taiyuan (ShiTai) high-speed railway line is designed to cross the Taihang Mountains. The line will include the construction of a 27.85 km long tunnel with a maximum depth of 450 m. The Taihang Mountain Tunnel is the key to the successful construction of the railway. The geological structure of the tunnel is highly complex and the lithology changes dramatically throughout. Fault zones, karst, and groundwater ingress are the main geological engineering problems for the tunnel construction. However, conventional geophysical methods are not effective for detecting these geological structures in a region with such a complex terrain. Therefore, the controlled source audio-frequency magnetotellurics (CSAMT) method was used to effectively characterize the resistivity distribution of the underground rock. This technique is highly sensitive to low resistivity geological structures, such as fault zones, and it is effective for studying the geological engineering conditions of the tunnel. Various corrections, designed denoising, and inversion methods are used in CSAMT data processing to obtain the resistivity profile of the entire tunnel. Based on the resistivity results and existing core information, the underground lithological distribution is inferred. The location of the fault zone is determined and the distribution of dangerous geological bodies is delineated. The interpretation of the CSAMT results is verified by boreholes and related geological information. These results are important for understanding the subsurface geological engineering conditions and have provided stability, safety, and cost evaluation support for the construction of a long deeply buried tunnel in Taihang Mountains. • CSAMT method was used to detect the deeply buried long tunnel. • Weak geological structures of the tunnel were inferred based on the sounding result. • The inferred faults were verified by boreholes. [ABSTRACT FROM AUTHOR]

The Tianshan-Xingmeng Orogenic Belt (TXOB) is located in China's northern region, and is one of the country's most important Au-Pb-Zn-Ag-Mo-Cu metallogenic belts. This metallogenic belt is known to have good minerogenetic conditions and exploration potential. However, the engineering detection depths of the majority of the metal deposits are generally less than 300 m. Therefore, high precision geophysical detection techniques are necessary to achieve the required accurate detection results. Among the commonly used techniques, the Stratagem EH4 system has been found to be the better choice for the deep prospecting of the metal deposits in the TXOB. Furthermore, a good understanding of the physical characteristics of the different rock formations in the region is essential for effective and efficient geophysical prospecting. This study revealed that the resistivity values of the ore and rock formation walls were different. In previous related research, the EH4 Conductivity Image System, which is a system combining MT and CSAMT, has been successfully used for ore explorations in China. The results of the geophysical prospecting processes indicated that there was good exploration potential of the aforementioned system in the TXOB. In the current study, the hidden gold ore bodies of disseminated (altered rock-type) gold deposits within the study area were defined by resistivity values of less than 150–1000 Ωm. The concealed gold orebodies in the quartz-type gold deposits were confirmed by resistivity values of 300 Ωm. The mineralization anomalies of the epigenetic lead-zinc deposits were defined by the smaller resistivity values, which were lower than the 100–1000 Ωm range. [ABSTRACT FROM AUTHOR]

The semi-airborne transient electromagnetic method (TEM) has been widely used in resource exploration and other fields due to its high efficiency, strong adaptability to the complex environment and high resolution. However, the conventional semi-airborne electromagnetic method continues to use the survey area of the grounded-wire TEM method. Compared with the existing ground-based TEM method, the semi-airborne TEM method is more sensitive to signal intensity because observed data is easily affected by noise for receiving position in the air. Therefore, it's unreasonable to consider only signal intensity or sensitivity related to offset in survey area. To obtain more suitable area, we discuss the distribution of the transient electromagnetic field with electrical source as well as the attenuation characteristics of the signal at different flight altitude through the forward modeling. Then, we obtain the sensitive area and flight altitude applied by the error analysis method in typical models. Finally, we describe a method what is herein dubbed the fuzzy comprehensive evaluation, which greatly reflects the membership degree of the evaluation object from different single factors to each level of fuzzy subset. Therefore, we use this method to get the optimal survey area and flight altitude under different noise levels. In low noise level, the signal intensity and sensitivity account for 33.3% and 66.7% respectively in selecting the optimum area and altitude, resulting in that the optimal area and flight altitude are coincides with the sensitive area and altitude better. While in the area filled with noise, the most area of strongest signal intensity is considered the optimal area and altitude. A field example shows that the aforementioned method can be well resolved in identifying survey lines and improving data quality. • The semi-airborne transient electromagnetic method (TEM) has been widely used. • We use the method of fuzzy comprehensive evaluation to evaluate the analog observation area with low and high noise. • The results are verified by field measurements. [ABSTRACT FROM AUTHOR]

The multichannel transient electromagnetic method (MTEM), sensitive to both near-surface and deeply buried geoelectric structures, is a useful geophysical tool to map resistive targets. A typical MTEM survey consists of a large number of receiver and transmitter positions. The produced data volume is much larger than a traditional transient electromagnetic method survey, posing computational difficulties for its data inversion and interpretation. We proposed a 3-D parallel inversion algorithm for MTEM data based on nonlinear conjugate gradient method and the moving footprint technique. The electromagnetic response and Fréchet derivative are first calculated in the frequency domain with integral equation method and then cosine transformed into time domain. We not only fully parallelize the forward modelling and Fréchet derivative calculation but also use a moving footprint to speed up the inversion. The footprints for different offsets are carefully evaluated through synthetic models. The effectiveness of the developed algorithm has been demonstrated through synthetic and case studies. [ABSTRACT FROM AUTHOR]

In underwater electromagnetic communications, geophysical electromagnetic exploration, and other fields, the conduction current in the medium is much larger than the displacement current. In this case, the electromagnetic waves follow the diffusion field equation. For the simulated region truncation of the diffusion field, traditional perfectly matched layers (PMLs) for the wave equation in media that are dominated by displacement currents are not applicable. In this article, a new PML metric for the diffusion field is proposed, and its parametric form is given, which is completely different from the conventional PML formation. The parametric form is related to the frequency and conductivity of the medium and is suitable for wide frequency band and media with large conductivity ranges. The spatial distribution of characteristic parameters and the method to find the optimal characteristic parameters are introduced. The numerical experiments for 1-D, 2-D, and 3-D radiation or scattering problems verify the validity and precision of the PML through a comparison with traditional methods or analytic solutions. Through the truncation performance of different models, it is concluded that the proposed PML metric is stable and can provide a good truncation performance for the diffusion field. [ABSTRACT FROM AUTHOR]

Methods based on pseudo-random coded emission technology can effectively suppress the noise in electromagnetic exploration and improve the depth and resolution of the detection. Additionally, the synthesis of electromagnetic response signals excited by pseudo-random code sources can provide not only theoretical data for data processing but also intermediate results for evaluations of field data. Currently, pseudo-random signal simulations ignore the presence of source effects generated by band-limited sources in the Earth's impulse response. Furthermore, the simulation methods are mainly completed based on one-dimensional and two-dimensional formats and thus cannot accurately reflect the pseudo-random response of three-dimensional geoelectric structures. In this paper, 3D pseudo-random electromagnetic signals with a band-limited source were synthesized according to the propagating routine of the electromagnetic field. First, Maxwell's Equations were solved using a three-dimensional time-domain finite volume method. To avoid the singularity of the source and incorporate source effect response, the Gaussian function was introduced to simulate the band-limited impulse sources. The Earth's impulse response was obtained, which contained source effects under three-dimensional geoelectric structures. Then, ideal pseudo-random electromagnetic signals could be composed by convoluting the Earth's impulse response with a pseudo-random source waveform. Since the transmitter and receiver instruments were band-limited, a low-pass filter with a Blackman window was introduced to simulate the instrument's system response. The simulated pseudo-random signals were obtained by adding noise and the ideal pseudo-random electromagnetic signals were filtered using a low-pass filter. By comparison with signals acquired in the field, the correctness of the pseudo-random signal synthesis with source effects was successfully verified. • MTEM response is composed of band-limited source, instrument response, and noise. • The practical earth impulse response includes the ideal response and source effect. • Gaussian sources can be used to directly simulate the band-limited impulse source of MTEM. • The source effect must be considered in the pseudo-random response simulation. [ABSTRACT FROM AUTHOR]

Three-dimensional (3-D) inversion technique has become an important and practical approach for magnetotelluric (MT) data interpretation. In this article, we developed a 3-D parallelized MT inversion scheme using the edge-based finite element method and applied the developed method to the newly collected MT data in the Xinjiang Luntai area. The distorted hexahedral element is adopted to incorporate topography into the forward modeling and inversion for complicated scenarios. We use the Gauss–Newton optimization method to minimize the objective functional for MT inversion. The developed algorithm is parallelized using MPI over frequencies and parallel direct solvers when solving the forward and adjoint problems for each frequency. We compare the performance of the least-square QR (LSQR) factorization and preconditioned conjugate gradient (PCG) solvers for the model update within each Gauss–Newton iteration and found that the LSQR solver is more stable. The developed inversion algorithm is validated using several synthetic models. Finally, we applied the inversion algorithm to the subsurface resistivity imaging in the Luntai area. The recovered geoelectric model from full 3-D inversion fits well with the known geological and geophysical information. The recovered model shows a low resistivity layer which may be caused by the salt strata. Besides, the inversion results reveal the movement tectonic in this survey area within a depth of 9 km. [ABSTRACT FROM AUTHOR]

We have carried out microseismic forward modelling in orthorhombic poroelastic medium. The seismic response to isotropic (ISO) source, double couple (DC) source and compensated linear vector dipole (CLVD) source was solved by finite-difference algorithm in time domain, respectively. We used analytical solution to test our numerical algorithm and find that the match between the analytical solutions and numerical solutions is sufficiently close. We then designed several different models to calculate the seismic response and analyse the effect of medium parameters and source mechanism on the propagation of seismic waves. We observed from the modelling results that ISO source excites two kinds of waves in isotropic medium, which are fast longitudinal wave and slow longitudinal wave, while the DC and CLVD sources excite three kinds of waves in isotropic medium, including a transverse wave in addition to the fast and slow longitudinal waves. All of these three kinds of sources generate four kinds of seismic waves in orthorhombic poroelastic medium. These are two separable transverse waves in addition to the fast and slow longitudinal waves. The fluid viscosity and medium tortuosity has an effect on the propagation of the slow longitudinal wave. The slow longitudinal wave appears in propagating mode under the effect of low fluid viscosity and appears at the source location in a static mode under the effect of high fluid viscosity. The wavefield snapshot of the slow longitudinal wave has a circular shape when the tortuosity is isotropic, whereas it has an oval shape when the tortuosity is anisotropic. As far as the anisotropic parameters are concerned, the anisotropy of the fast longitudinal wave is more sensitive to the value of |$\varepsilon ({{\varepsilon _1},{\varepsilon _2}})$|⁠ , while the value of |${\rm{\delta }}({{{\rm{\delta }}_1},{{\rm{\delta }}_2},{{\rm{\delta }}_3}})$| has more effect on the anisotropic behaviour of the two separable transverse waves. [ABSTRACT FROM AUTHOR]

• Three-dimensional electrical resistivity structure beneath Ganzi Basin, eastern margin of Tibetan Plateau was obtained. • 3D electrical model reveals the shallow and deep geothermal reservoirs in Ganzi geothermal system. • The Ganzi fault provides the pathway for transport heat and the circulation of geothermal fluids. The Ganzi geothermal system located in the eastern margin of Tibetan Plateau, is found to distribute the strongly active surface manifestations, such as boiling springs, and fumaroles. The resistivity structure of high-temperature hydrothermal system beneath Ganzi Basin has been obtained from the three-dimensional (3-D) magnetotelluric inversion of 103 broadband MT stations in period range from 320 Hz to 1000s. The MT stations have an average spacing of 500 m and completely cover the Ganzi Basin. The detailed 3-D resistivity structure which the authors found has been described characteristics of the shallow and deep geothermal reservoirs as well as to correlate with the structural features present in the Ganzi Basin. The shallow geothermal reservoir is mainly composed of the sedimentary layer, with a depth of 300−800 m, while the deep geothermal reservoir mainly consists of the fault fracture zone with a depth of 1−5 km. The porosity of the geothermal reservoirs is in the range 5–10 %, estimated based on the resistivity model. The Ganzi fault provides the pathway for transport heat and the circulation of geothermal fluids. The possible heat source for the Ganzi geothermal system has been attributed to tectonic deformation heating in the crust, radiogenic heat production derived from the decay of U, Th and K, and heat from the mantle. [ABSTRACT FROM AUTHOR]

The paleokarst system is one of the main carbonate reservoirs, which can form important super-large oil fields. There are many typical paleokarst reservoirs in the Tarim Basin Ordovician strata, mainly composed of caves, vugs, and fractures. Due to the deep burial depth and strong heterogeneity, qualitative identifying the different scale fracture-vuggy reservoirs from the tight limestone around the wellbore is a real challenge in the industrial community. In this paper, machine learning methods were used to classify electrofacies. Firstly, core samples and electrical imaging logging of the paleokarst reservoirs are observed in detail and a core-electrical imaging chart was established. Secondly, conventional logging data was optimized and preprocessed for data mining, using Principal Component Analysis (PCA) algorithm and K-means algorithm. High-resolution electrical imaging logging was chosen as a constraint to recognize electrofacies, and an electrofacies-lithology database was established. Thirdly, based on the electrofacies-lithology database, Linear Discriminant Analysis (LDA) algorithm was used to build an electrofacies prediction model, which can automatically identify the electrofacies in carbonate strata, with a coincidence rate of 92.2%. Finally, the model was used to quantitatively recognize paleokarst reservoirs and their distributions. The electrofacies machine learning workflow proposed in this paper could be used in Tarim Basin and other similar paleokarst reservoirs, which can improve exploration efficiency and save economic cost. • Machine learning methods on conventional logging data can recognize electrofacies. • High-resolution cores and FMI calibrated the electrofacies recognition process. • Electrofacies recognizations can reveal heterogeneous paleokarst characteristics. [ABSTRACT FROM AUTHOR]

Abstract For controlled source audio frequency magnetotellurics (CSAMT) which supported by traditional dipole hypothesis theory, the observations point should be located in far source zone. However, during field measurements, due to limited transmitter power, the distance between transmitter and receiver is usually not too large. Sometimes, field measurements are carried out in the near source zone, and the traditional dipole hypothesis will produce large errors under such condition. So, the errors between the approximate solution of traditional dipole approximation and the exact solution of the long grounded-wire source are firstly studied in this paper. Moreover, relative error distribution between exact solution of the long grounded-wire source and electric dipole source variations with different frequencies has been analyzed. It finds that the traditional dipole hypothesis bring large errors in the near source zone and no longer be satisfied to such condition. As a further study, a new 'harmonic point charge' solution for CSAMT, instead of the traditional dipole hypothesis, is proposed in this paper, which lay a foundation to realize near source prospecting and obtain more accurate near source responses in frequency domain electromagnetic. [ABSTRACT FROM AUTHOR]

Because of its unique geographical location, the Tibetan Plateau has been widely studied in geological research. Several studies revealed the geological structure and evolutionary process of the lithosphere deep beneath the northeastern Tibetan Plateau; however, the deep structural characteristics of the Zhangye Basin, which is in the middle of the Hexi Corridor at the junction of the northeastern Tibetan Plateau and Alxa Block, have not been widely examined. To better understand the geodynamics of the Zhangye Basin, studies are needed to investigate the tectonic mechanism of the northeastern Tibetan Plateau. To determine its subsurface structure, 46 broadband magnetotelluric survey stations were deployed across the basin in 2020. The crustal electrical structure of the basin was modeled via three-dimensional nonlinear conjugate gradient inversion. The results of this inversion revealed a low-resistivity zone in the middle-lower crust of the basin, which was characterized as an intracrustal decollement, thus suggesting decoupling deformation of the crust interior. This conductive body may have formed during continuous growth of the Tibetan Plateau, which may have expanded northeastward below the Hexi Corridor and Helishan Fault. The latter may have partly underthrust southwestward in the lower part of the Hexi Corridor and Tibetan Plateau. These results regarding geothermal formation of the Zhangye Basin improve the understanding of its underground structure and may guide the future development of geothermal resources in this area. • The Zhangye Basin was modeled using three-dimensional inversion. • A conductive body in the middle and lower crust indicates intracrustal decollement. • The conductive body may have formed during the expansion of the Tibetan Plateau. • The Zhangye Basin has the potential for geothermal development. [ABSTRACT FROM AUTHOR]

Laterally constrained inversion (LCI) has been successfully applied to the inversion of dc resistivity, TEM and airborne EM data. However, it hasn't been yet applied to the interpretation of controlled-source audio-frequency magnetotelluric (CSAMT) data. In this paper, we apply the LCI method for CSAMT data inversion by preconditioning the Jacobian matrix. We apply a weighting matrix to Jacobian to balance the sensitivity of model parameters, so that the resolution with respect to different model parameters becomes more uniform. Numerical experiments confirm that this can improve the convergence of the inversion. We first invert a synthetic dataset with and without noise to investigate the effect of LCI applications to CSAMT data, for the noise free data, the results show that the LCI method can recover the true model better compared to the traditional single-station inversion; and for the noisy data, the true model is recovered even with a noise level of 8%, indicating that LCI inversions are to some extent noise insensitive. Then, we re-invert two CSAMT datasets collected respectively in a watershed and a coal mine area in Northern China and compare our results with those from previous inversions. The comparison with the previous inversion in a coal mine shows that LCI method delivers smoother layer interfaces that well correlate to seismic data, while comparison with a global searching algorithm of simulated annealing (SA) in a watershed shows that though both methods deliver very similar good results, however, LCI algorithm presented in this paper runs much faster. The inversion results for the coal mine CSAMT survey show that a conductive water-bearing zone that was not revealed by the previous inversions has been identified by the LCI. This further demonstrates that the method presented in this paper works for CSAMT data inversion. [ABSTRACT FROM AUTHOR]