In the light of typical coneshaped columns in the Lu'an mine of Shanxi Province, the authors

established a mathematical model for the collapse column, and used the wave equation model for

seismic ray tracing and wavefield simulation of collapse columns. The simulation results show that,

due to the special nature of the collapse columns, such waves as the normal reflected wave, the fault

point diffraction wave, the delay diffraction wave, the delay reflected wave and  the "diffraction

wave" consisting of Pwave field are formed around the collapse columns. On such a basis, the field

shot was simulated, and then a single shot record for routine processing  was generated to produce

stacking sections and migrated sections. The forward section and the actual data analysis reveal that

the fallen columns of the special wave constitute an important feature for recognizing collapse

columns; nevertheless, the resultant scale of the actual seismic data interpretation is often smaller

than the scale of actual collapse column. Some suggestions are also put forward for reference.

Low density geochemical survey in Yishenjilike mountain area led to the discovery of a huge gold geochemical block, and the verification and evaluation of gold anomalies led to the discovery of the Bogutu gold deposit. In combination with the geological characteristics of the prospecting area, the authors carried out a series of geological-geophysical-geochemical exploration work, delineated quite a few geochemical and IP anomalies, and detected the characteristics of ore-bearing structural belt. Through trenching and drilling verification, the authors found more than 40 gold orebodies, thus achieving good ore-prospecting results.

The existing geological data show that there are several buried faults in the main urban area of Chengdu. However, the specific location and distribution of these faults are still unclear, which poses great security risks to the comprehensive and scientific exploitation and utilization of underground space resources and the optimization of urban construction planning and layout in Chengdu. In view of such a situation, four geophysical methods, namely, micromotion survey, high-density electrical method, transient electromagnetic method and soil radon measurement, were used in this paper to comprehensively explore the buried Baojiangqiao fault in the work area. The integrated geophysical prospecting methods not only identified the stratigraphic structure along the survey line, but also obtained the location, property, attitude and scale of the buried Baojiangqiao fault. This work indicates that the integrated geophysical prospecting methods can achieve better results in the exploration of urban buried fault.

A review on the airborne gravity survey and airborne magnetic survey are given in this paper. The

authors hold that these techniques have entered into a new development period, as evidenced by the

emergence of such technologic indicators as airborne vector magnetometry, magnetic fulltensor LTSSQUID

and HTSSQUID gradiometry, airborne gravimetry, gravity gradiometry and geophysical survey by UAV.

This paper gives a brief review of the history of the airborne gravimetry. Based on the principle of the airborne gravimetry, the paper deals emphatically with the history, the present  research situation and the developments of the airborne scalar gravity survey system. Further development trends are also indicated.

 With powerful Matlab image processing functions, this paper has realized the 3D display of ground penetrating radar data, whose procedures are simple to prepare and easy to learn. A detailed description of the code is given, and the three-dimensional test data show that, by setting transparency, the target can be displayed in an intuitive and visual way.

The R-mode cluster analysis is a mathematic statistical method for obtaining the quantitative similarity of several elements. Its procedure includes: the conversion of the original data; the solution of the relevant coefficient ; the clustering of the result. The above operation can be realized by using the data analysis tool of Excel. This method is quite suitable for field utilization.

This paper puts forward a new method for calculating anomaly area—the polygonous approximation method,and gives corresponding algorithm and precision controlling technique.With this method,one can calculate anomaly area rapidly and precisely,thus yielding relatively satisfactory results.

In drawing the contour map of geophysical data, suitable gridding methods should be chosen according to objective

environment characteristics and characteristics of data themselves. With practical examples, this paper deals with some common

gridding methods such as inverse distance to a power, Kriging, Minimum Curvatrue, Nearest Neighbor, Polynomial Regression, Radial

Basis Function, and Triangulation/Liner Interpolation, probing into their choosing means, applicable fields and parameter

installment.

This paper has described and reviewed the pinciples,expermental set up and observational results of applying atom interferometer to measuring the earth's gravitational acceleration.This method was developed by Nobelist Steven Chu,M.Kaservich,A.Peters et al.,who demonstrated a resolution of 10^-10.The application of this technique to geophysics and the related problems are also discussed.

IP is an important method for mineral resources exploration, but it can only be applied in a few pivotal places because its difficult performance, and hence the information obtained from the survey area is very limited. The aim of this study is to find a simple and effective IP method with which we can get much valuable information form the survey area so as to improve the exploration effect. In the Yinan gold ore district, the authors applied the gradient sounding profile method, calculated 2D resistivity and obtained the IP model. The results show that the gradient sounding profile method is simple and effective.

Element abundances of the continental crust and rocks cited in the past and recent literature have been based on compilations of data from various studies. This leads to uncertainties in sample representativeness and data quality. The present study is based on systematic collection of 28 253 individual rock samples over an area of 3.3 million km² in eastern China, east of eastern longitude of 105°.The sampling involves more than 800 igneous intrusive bodies and metamorphic complexes as well as more than 500 type stratigraphic sections. From the individual rock samples, 2 718 composite samples were made and analyzed for 77 elements by 15 specific methods, dominantly XRF and INAA. Analytical quality was controlled by international and national preliminary geochemical reference materials of GSR, GAu and GPt series. Synthetic results from geological, geochemical and geophysical studies were used to construct crustal constitution model, from which element abundances of the continental crust in the North China platform, the upper crust and exposed crust in eastern China, chemical composition of igneous rock in China and of sedimentary rock and metamorphic rock in eastern China were derived.

As an important prospecting technique and an effective means for obtaining mineralization data,

geochemical exploration has been proved to be very successful in mineral exploration. This paper focuses on

commenting the present research situation and application effect of some new methods such as structural

superimposed halos method, heat released mercury method, separatory electrogeochemistry method, enzyme leach, and

geogas and selective leaching of mobile metals method. It is emphasized that any one of these methods has its

unique applicability and that, in the practical application, we should pay attention to the cooperation of

geochemical exploration, geologicalgeophysical exploration and remote sensing and depend on the study of

geological background so as to demonstrate the usefulness and effectiveness of geochemical exploration.

In this paper, recent advances in the research and development of quantum magnetometers have been described together with some suggestions on further research work.

Based on two-dimensional Fourier transform and half-wave theory, this paper has studied the seismic Rayleigh wave dispersion curve extraction in f-k domain and made this theory fit in with a program by means of Delphi7.0. It is concluded that the f-k method overcomes the shortcomings of the one-dimensional digital processing technology and makes full use of multi-channel Rayleigh wave data record.

Porphyry copper deposit is the most important copper deposit type in China. With the deepening of mineral exploitation, exploration geochemistry in mineral exploration has become increasingly important. Based on related literature, this paper sums up the exploration geochemical research results of porphyry copper deposits, which include such aspects as geochemical characteristics, exploration methods, anomaly evaluation and prospecting indicators. Exemplified by the Fujiawu copper deposit, this paper reports the latest advances in the study of porphyry copper deposits.

According to aeromagnetic and gravitational data, the boundary and range of Ordos basin were determined based on an analysis of the characteristics of gravity and magnetic fields. And on the basis of compiling depth map of the metamorphic basement and structuring zoning map, research was conducted on the basin’s basement structure, characteristics and features of depth change, structure framework, and caprock thickness. The research indicates that the metamorphic basement is composed of Proterozoic metamorphic series, and the buried depth of crystalline basement can reach 5 000 to 20 000 meters. The caprock is the layers of Ediacaran, Paleozoic, and Mesozoic. The occurrence and development of the basin are restricted by nearly NE-and NWW-trending structures, forming a pattern of four depressions and three uplifts. All the new understanding and conclusions provide a reference for further oil and gas exploration in the basin.

According to the data obtained from air-borne electromagnetic and magnetic survey and

ground two-frequency IP and magnetic survey, this paper discusses the relationship between the

features of air-borne and ground geophysical composite anomalies and the ore resource distribution.

On the basis of an integrated analysis, geophysical composite anomaly criteria in search for iron, zinc

and lead deposits were established, and 12 important ore-prospecting targets were delineated, which

provides very important clues for finding iron, zinc, lead and some other ore resources.

 Numerical integration of Hankel transforms is effective tools for EM Sounding 's forward numerical simulation, this paper made out formula derivation of numerical integration of Hankel transforms by digital filtering, and use digital coefficients to do numerical compute which was put forward by Guptasama and Singh，finally contrasted to theoretical resolve expression and analyzed this algorithm's error distribution. The results show that the calculation of this algorithm continuously approximate its theoretical solution, it has no oscillation, high precision and great practical value in the numerical simulation study.

 Vibroseis exploration, as an important method of seismic exploration, has become increasingly valued by geologists. During the field work, different geological conditions need different parameters, and the choice of suitable excitation parameters has become a very important problem. This paper mainly introduces six kinds of parameters, namely numbers of controlled seismic source, scanning bandwidth, vibration period, scanning length, scanning slopes, and vibration rate. The influence of these parameters on the quality of seismic records was studied in detail, and the simulation of the proper parameters was performed in Jiulishan area to enhance the resolution of vibroseis exploration and improve the signaltonoise ratio of seismic records.

On the basis of accurate analysis of 154 composite soil samples collected from various landscape in China and other reference information, abundance of 79 elements and composition of soils in China have been given, chemical composition of soils under the influence of matrix rock and supergenesis geochemistry condition have been studied preliminarily.

This paper has put forward the idea of utilizing Surfer to draw the initial data profile-plan. The programme compiled can realize the construction of Bln file. The profile-plans for regular net and irregular net have been constructed, and the results are satisfactory.

Migration velocity analysis and migration imaging constitute two important parts in seismic data processing. At

present, time migration has become mature, and depth migration is getting more and more perfect. The common method for time

domain migration imaging is pre-stack time migration. By adopting stacking velocity analysis along the layer, exact layer

stacking velocity can be obtained. Through dip correction, pre-stack time migration and CRP de-migration velocity analysis, the

velocity is optimized step by step, and then a desired RMS velocity field that is up to the geological rule is acquired. In

addition, the method for building the accurate migration velocity field is summed up through the study of the depth migration

method. A new seismic data processing flow in combination of Kirchhoff arithmetic based on ray tracing and wave equation

arithmetic based on wave field extrapolation is presented. As a result, migration velocity analysis and migration imaging are

considerably improved.

The implied contour level and color scheme in Surfer software fail to express the subtle difference of DEM and the effect

of color solid. In addition, the artificial setting of the levels and color values are timeconsuming and laborious, and the

results are sometimes not perfect. With the consideration of contour level and color scheme and on the basis of Surfer platform

automation technology, this paper presents a program that can automatically generate a level file with different intervals and

different colors, thus resulting in satisfactory effect and efficiency.

On the basis of mote than230 English and Russian papers or monographsissued in the last 30 years as well as nearly 50 Chinese articles published overthe past 10 years on the inversed problem in g avity and magnetic explora-tion, combined with his practice in the study of the inversed problem, the all-thor expounds the inversed methods in gravity and magnetic exploration cur-rently used both at home and abroad, briefs the readers about the researchsituation in China and foreign countries, and makes a detailed review on achi-evements gained in China in comparison with the research levels abroad.Final-ly, suggestions are put forward concerning some subjects which are worthnoticing in future research work.

This paper has summed up the progress of the ERI method over the past decade of years as well as its future development trend in the following aspects: ① A comparison of the performances of the main ERI instruments used at present shows that the ERI instruments tend to develop in the multi-channel, multi-parameter, multi-functional, high-power direction; ② ERI measurement environment has changed from surface measurement to water surface, underwater and cross-hole measurements, with the last three kinds of measurements analyzed in this paper; ③ On the basis of analyzing ERI data processing method and inverse development status, this paper describes three-dimensional and four-dimensional inversion theory of ERI with practical examples; ④ ERT applications are summed up, and several new applications are introduced. It is concluded that, with the improvement of the probing depth and observation precision as well as the diversification of the observation models, the application field of ERI will become broader and broader, and this technique will surely have wide development prospect.

This paper proposes an equivalent second order elastic wave equation to solve the problem of being unable to completely

separate the coupled P and S wave by full elastic wave equation. Through solving this equivalent wave equation by high-order

staggeredgrid finite difference scheme together with Flux Correction Technology (FCT) and separating wave fields of the

isotropic model and layered model, the authors accurately obtained the hybrid wave field and completely separated fields of pure

P wave and pure S wave fields. An analysis of the numerical results shows that the method is effective and reliable in isotropic

media, and there exists abundant energy transform information in separated pure P and pure S wave field. The result of the study

is of significance in understanding the propagating law and the elastic wave theory in the complex wave field.

The development of the gravity gradiometer is described simply in this paper, and the measuring principle and development experience of the rotating accelerometer gravity gradiometer are emphatically discussed. On the basis of researches on the existing airborne gravity gradiometer, the prospects of the airborne gravity gradiometer are also presented.

With the extensive adoption of computer technology in geology and geophysics, computer graphics becomes more and more important. For the purpose of improving quality and efficiency, several software should be combined in practical work. In this paper, some basic characteristics of Surfer and MapGIS are discussed, and their application in geological field is described from the angle of geological mapping, thus drawing forth the necessity of the conversion between Surfer and MapGIS. The general method and steps for the conversion of maps between Surfer and MapGIS are presented in detail, and several problems concerning the conversion and corresponding solutions are emphatically discussed.

The algorithm of the finite difference oforder velocitystress staggered grids has been built, which is

suitable for the perfectly matched layer (PML) absorbing boundary for elastic media. Meanwhile, the construction of the

PML absorbing boundary condition and the realization of the finitedifference algorithm are discussed in detail. Wave

field modeling calculations show that, compared with the conventional decaying exponential absorbing boundary and non

absorbing boundary, the PML absorbing boundary can perform absorption much more clearly and absorb the boundary

reflection from various angles, whose absorptivity (the ratio of absorption energy to unabsorbable energy) can reach

99.99%. The PML absorbing boundary can well eliminate the periodic folding effects, which makes the calculation of the

wave field characteristics very distinct, and the Rayleigh wave can be clearly shown in the waveform record.

 Basic principles and techniques for the shallow seismic refraction wave method and the highdensity resistivity method

are described in this paper in combination with engineering practice to explain the results of the refraction wave method and the

high-density resistivity method in railroad and highway tunnel exploration. The application effects of the combination of these two

methods in detecting positions, sizes, natures and characteristics of such harmful geological bodies as concealed faults and shatter

zones are also discussed. Some measures for improving the investigation accuracy of these two methods are put forward.

In the light of the characteristic time-frequency properties of different time-frequency analysis methods, it is

important to select an optimal time-frequency analysis method for high precision seismic exploration in the coal field on the

basis of comparing and studying these methods. Aimed at probing into STFT, Wavelet Transform, S-Transform, Choi-Williams

Distribution, Wigner-Ville Distribution and its improved methods, this paper studied advantages and disadvantages of every method

through numerical simulation in time resolution, frequency resolution and interference. The First Break Picking and Separate

Frequency Display for Collapse Column interpretation were realized by SPWVD and S-Transform, with a good result obtained.

The main minerals in the Yiliu polymetallic ore district are mainly As, Sb, Pb, Zn, Ag, W, Sn, Nb and Ta, the deposit types and the ore-controlling factors are complex and diverse. Recently, by using 1:50000 geochemical stream sediment survey, geochemical soil survey and CSAMT method in this ore district, supergene distribution of elements were found, and many integrated anomalies were delineated. A comprehensive evaluation of all kinds of anomalies shows that the F3 fault has features of ore-passage structure and ore-storage structure, and the area sandwiched between F2 and F3 fault are favorable ore-forming place, whereas the depression parts of biotite granite which have high resistance and multiple sets of fracture intersections are most favorable places for mineralization.

In the light of the pulse width of the radar wave, this paper deals with the difference and the relationship between the vertical resolution and the transverse resolution of the ground-penetrating radar, indicates the influence of the noise upon the resolution and, with practical examples, points out that the digital handling of the signal can greatly improve the resolution of the radar.

In this paper,a historical review of the gravity gradiometry is given, the application of this technique to oil exploration and some other fields is described,and its development in future is predicted.

With practical examples, this paper deals with the basic principle and technique of the Rayleigh surface wave exploration, the technology for its application to railways, the data acquisition, processing and interpretation, and the calculation method for converting the velocity of the Rayleigh surface wave to the bearing capacity. The reliability of applying this technique to railways is also pointed out.

This thesis supply a profile about performance of TerraTEM system and the  distinguishing feature of TEM. By exploring

work on two conceal iron deposit in HeNan province and analyzing. the abnormal value of the TEM system.combining with drill

drawing,The writer prove that the TEM can take a obvious effect on exploring iron mineral.

In this paper, the working principle of three satellites (CHAMP, GRACE and GOCE) has been described, and the important geopotential models and the newest geopotential models (EIGEN-CG01C, GGM02) have been introduced.

As a kind of clean energy, geothermal energy has attracted the attention of scholars all over the world in recent years. Previous geochemical exploration methods for geothermal resources are limited to the analysis of individual geochemical indices. Moreover, previous studies of microorganisms in geothermal fields mostly focus on hot spring outcrops, lacking ecological studies of geothermal resources in complex terrains. This study investigated the soil geochemistry and microbial diversity of the Bantang Hot Spring geothermal field in Chaohu, Anhui Province. Geochemical indices such as head-space gas, soil gas, acid-hydrolyzed hydrocarbons, and altered carbonate were detected in this study. Combined with the microbial high-throughput sequencing technology, this study analyzed the composition and spatial-temporal distribution of the microbial population above the geothermal fields in uplifted mountains and the relationship between these bioinformatics characteristics and the geochemical indices. The results indicate that the acid-hydrolyzed hydrocarbons on the surface of the geothermal field showed a maximum methane concentration of 43.7 μL/kg in the area between faults F2 and F3, adequately reflecting the fault location of the geothermal field.Bacillaceae, Hydrogenophilaceae, and Thermodesulfovibrionaceae in the geothermal field and the background area showed large relative abundance differences, which were 0.178%, 0.108%, and 0.060%, respectively. This result indicates that they are sensitive to geothermal resources and correspond well to geochemical indices above the known geothermal field. This study preliminarily investigated the diversity of geothermal microorganisms in the geothermal field and analyzed the corresponding relationships between microbial distribution characteristics and geochemical indexes, providing technical support for the microbiological exploration of geothermal resources.

Following China's planning for ore prospecting in overburden areas, China's geochemical exploration researchers have conducted extensive research on the fundamental theories, methods, and technologies of geochemical exploration for overburden areas in the past decade. They achieved significant advances mainly in two aspects: (1) the research on the migration mechanism, occurrence state, and anomaly formation mechanism of elements in overburden areas; (2) advances in methods and technologies, including geoelectrochemistry, active state of elements, geogas, separation of micro-fine-sized soil particles, soil thermomagnetic composition, and integrated gas survey, as well as numerous experimental demonstrations. These advances represent continuous progress in the research on the fundamental theories, methods, and technologies of geochemical exploration for overburden areas, providing new geochemical methods and technologies for ore prospecting breakthroughs in overburden areas.

The Fengtai ore concentration area is an important producing area of plumbum-zinc and gold ores in the middle of the Qinling orogenic belt. It hosts many large to super-large deposits, such as Qiandongshan-Dongtangzi, Bafangshan-Erlihe, Baguamiao, and Shuangwang deposits. With the exploitation proceeding, the reserves of these deposits have decreased significantly, and the ore prospecting in these deposits has shifted from the surface to the deep part. However, the geophysical fields in the deposits are yet to be ascertained, severely restricting research on the metallogenic regularity of the deposits and the ore prospecting and exploration in the peripheral zones. Based on the latest 1∶50,000 aeromagnetic and airborne radioactivity survey data, as well as gravity data, this study investigated the multi-source geophysical fields on the scales of the region, the ore concentration area, and deposits, aiming to summarize the distribution patterns of geophysical anomalies of different scales and provide evidence for research on metallogenic regularity and prospecting prediction. The results show that the Fengtai and the Xicheng ore concentration areas, with similar magnetic structures, are separated by the zone with strong magnetic anomalies caused by the crystalline basement of the Huicheng Basin. The first vertical derivative of gravity reveals that the Fengtai and Xicheng ore concentration areas have density structures similar to the Huicheng Basin. It can be inferred that the shallow part is a unified basin and that a large prospecting space exists between the two ore concentration areas. Many NW-trending linear magnetic anomaly zones occur in the Fengtai ore concentration area. Their locations are highly consistent with those of fault structures, and they were offset by NE-trending structures due to late transformation. The second vertical derivative of gravity reveals that many intermediate-acid intrusive stocks or veins have developed in the deep part of the Fengtai ore concentration area. Structural boundaries can be effectively identified based on the first vertical derivative, X-directional derivative, and wavelet transform of aeromagnetic data. Furthermore, ore bodies are mostly distributed in the NW direction along the tectonic belt. All these will play an important role in guiding ore prospecting and exploration in the Fengtai ore concentration area.

Geothermal resources are significant clean energy and tourism mineral resources. The Rucheng Basin, a carbonate basin in the southeastern mountainous area of Hunan Province, possesses favorable conditions for the formation of convective geothermal energy. However, the basin is enclosed on three sides by the giant Zhuguangshan rock mass, and its basement is subjected to the intrusion and destruction by the rock mass, resulting in severely deformed formations, crisscrossing faults, and significantly different eastern and western structures. The understanding of the basin's water- and heat-conducting pathways and deep reservoir structures remains elusive, thus restricting the investigation of the basin's geothermal potential. Hence, this study probed the basin's deep structures through gravity survey and audio magnetotellurics (AMT), obtaining the following insights: (1) The Rucheng Basin has developed into a bidirectional ramp structure due to east-west differentiation. The synclinorium in the east experienced compression and clockwise rotation due to the emplacement of the Yanshanian rock mass, rocks were fragmented in the core zone, and strike-slip fracture zones were found at the boundary. The faults have vertical cutting depths exceeding 4 km, widths ranging from 300~600 m, and dip angles between 80°~90°. (2) The basin's basement anticlinal axis hosts several NWW-directed concealed rock masses, with diameters from 3~4 km and buried depths from 0.5~1.5 km. Hot springs reside in the fracture zones crossing the boundaries of the concealed rock masses. (3) The basin boasts favorable conditions for the formation of convective geothermal energy. Folds, fault zones, and concealed rock masses match each other to form a unified spatial combination of heat-controlling elements, manifesting heat accumulation characterized by east-west recharge and intermediate discharge. With more thriving deep geothermal reservoirs in the east, the basin has high potential for geothermal resources.

This study aims to explore the feasibility of the carbon dioxide (CO₂) and sulfur dioxide (SO₂) gas geochemical survey method for mineral exploration in forested areas. Based on the newly designed gas rapid analysis instrument, this study conducted an experimental investigation of the method in the forested Jiapigou gold concentration area, Jilin Province. The results show that significant CO₂ and SO₂ anomalies were observed above the concealed ore bodies and structures. In the forested area, the CO₂ and SO₂ gas geochemical survey method reflected the fault structures and effectively indicated the deep concealed gold deposit. This method holds critical significance for the breakthrough of prospecting technology in China's covered areas.

With the rapid development of urban infrastructure,the demand for urban geological work is increasing,and urban geological surveys become particularly important.Urban geophysical exploration has different exploration purposes and working environments from conventional geophysical exploration.Accordingly,compared with conventional geophysical exploration methods,the geophysical exploration methods for urban geological surveys face the challenges of many interference factors,limited construction sites and time,and high requirements for exploration accuracy.The opposing-coils transient electromagnetic(OCTEM) method enjoys a strong anti-interference ability,convenient and efficient construction,and high resolution.Therefore,this study employed the OCTEM method to investigate the test profile in the urban geological survey and evaluation of Haidong City.This test profile was subjected to numerous interference sources since it crossed 11 highways and railways and passed through factories,schools,logistics parks,villages,living quarters,and rivers.Consequently,the OCTEM results agree well with the results of single-point resistivity sounding and drilling results.Therefore,the OCTEM method proposed in this study is effective for urban geological surveys.

Ji'nan possesses highly abundant geothermal resources, which are hosted by Ordovician-Cambrian karst-fissured geothermal reservoirs and Neogene-Paleogene clastic pore-fissure geothermal reservoirs. The geothermal exploration in this study focuses on the Ordovician-Cambrian karst fissured geothermal reservoirs in Daqiao Town in northern Ji'nan. Through geophysical profile measurements, this study aims to identify the distributions of strata and fault structures and the burial depths of geothermal reservoirs, infer the attitudes and spatial morphologies of fault structures associated with heat control and conduction, delineate the target area for geothermal well construction, and conduct drilling verification in the favorable underground water-rich position. Building on the collected data, this study interpreted and inferred the fault structures in the study area and comparatively analyzed the water-bearing properties by employing direct-current sounding, controlled source audio magnetotellurics, and magnetotelluric survey. A geothermal exploration and production combined well was constructed in a favorable position of the geothermal target area, manifesting a completion depth of 1 532.06 m, a static-water burial depth of 13.03 m, a wellhead water temperature of 50.1 ℃, a water yield of 132.998 m³/h, and a dropdown depth of 18.27 m.

The Eerguna block with metallogenic geological conditions is an important metallogenic area in Heilongjiang Province. Globally, China boasts the richest resource of antimony. However, the high mining intensity in recent years imposes huge challenges to this resource advantage of China. In this context, it is necessary to ascertain the geochemical characteristics of antimony in the Eerguna block. Based on the data of the 1∶250 000 stream sediment survey in the Eerguna block, this study explored the geochemical parameters of antimony in different tectonic units and the regional geochemical anomalies of this block. The results show that the study area has median and average concentrations of antimony of 0.33×10^-6 and 0.55×10^-6, respectively. The Mohe foreland basin is rich in antimony, with median and average concentrations of antimony higher than those of the study area. Furthermore, zones with high and extremely high antimony concentrations in the study area are distributed primarily in the Mohe foreland basin. Based on the 85% cumulative percentage, this study determined 66 geochemical anomalies of antimony, among which two reach the scale of geochemical provinces. Furthermore, this study identified significant geochemical anomalies of antimony in the discovered gold, antimony, and plumbum deposits or ore occurrences (mineralization points). Based on the spatial distributions of geochemical anomalies and metallogenic geological conditions of antimony, arsenic, and gold, this study delineated three metallogenic prospect areas of antimony: the Beijicun-Sanlianshan metallogenic prospect area, the Wangsushan-Daling metallogenic prospect area, and the Baikalushan-Huzhong metallogenic prospect area. In addition, the geochemical anomalies and metallogenic prospect areas for antimony, arsenic, and gold provide important areas for searching for sulfide deposits such as gold, antimony, and plumbum ones in the study area.

Despite abundant geothermal reserves of the Tangyu geothermal field in Meixian County, Shaanxi Province, long-term exploitation has decreased the water temperatures and yields of its existing geothermal wells year by year. Hence, there is an urgent need to explore new potential geothermal resources in the geothermal field. Since the known geothermal wells in the geothermal field are significantly controlled by faults, investigating the deep fault propagation holds critical significance for exploring the geothermal field’s potential geothermal resources. Due to the method limitations and the topographic influence, identifying thermal control faults through conventional geological route investigation or large-scale engineering is not applicable to the geothermal field. Therefore, a new technical method combining the penetrating soil radon measurement and the controlled source audio magnetotelluric (CSAMT) method was employed in this study to find concealed faults and delineate potential geothermal areas. Based on the measured surface soil radon concentration anomaly data and the subsurface electrical structure model derived from the CSAMT data inversion, this study inferred six new concealed faults on the basis of corroborating the known faults, predicted two potential geothermal areas, and built a conceptual model for the Tangyu geothermal field. As revealed by the results, the soil radon concentrations at concealed faults are much higher than the regional background value, and the concealed faults are located in the low-resistivity fracture zones as indicated by the apparent resistivity results based on CSAMT data inversion. Besides, the two potential geothermal areas spread from 450~750 m and 850~1 150 m on the profile, respectively, at depths of approximately 250~300 m. This study concludes that the geothermal field resides in a low-resistivity region with soil radon anomalies three times the regional background value. The results of this study provide a reference for the subsequent sustainable production and utilization of potential geothermal resources in the region.

The progress in seismic acquisition techniques characterized by wide azimuths,wide frequency bands,and high densities has greatly promoted the application of the prestack P-wave fracture inversion technique based on the azimuthal anisotropy theory.Azimuthal anisotropy-based inversion can yield the azimuths and intensities of fractures.However,different inversion techniques yield different parameters for fracture intensity characterization,resulting in inconsistent inversion results.Consequently,the azimuthal anisotropy-based inversion results of fractures are non-unique,leading to confusion about accurate results.Based on the Thomsen anisotropy theory,as well as the interrelationships between fracture models(the Hudson coin model and the Schoenberg linear sliding model),this study established the connections of anisotropic parameters between different fracture inversion techniques(VVAZ,Ruger's approximation,and Fourier series),presenting the real meanings and mathematical expressions of results from different azimuthal anisotropy-based fracture inversion techniques.Additionally,this study summarized the relationships of parameters between different inversion techniques and fracture models,further deepening the research on azimuthal anisotropy-based fracture inversion.This study lays solid theoretical and technical foundations for large-scale fracture detection based on the seismic data obtained using the seismic acquisition techniques featuring wide azimuths,wide frequency bands, and high densities.

The Laiwu area in central Shandong Province, situated in the eastern North China Craton, is a significant production area of skarn iron-rich ores. Its ore deposits occur primarily in the contact zone between the mining rock mass and the Middle Ordovician carbonate formation. Based on the latest areal gravity and magnetic survey results, this study thoroughly investigated the characteristics of gravity and magnetic anomalies along the Shijiaquan-Liujiamiao area in the western periphery of the mine rock mass. Then, this study delineated the deep prospecting target combining the characteristics of gravity and magnetic fields of the known iron deposits in the Laiwu area. Large-scale gravity and magnetic profiles were arranged in the favorable mineralization area. With the known boreholes as constraints, the gravity and magnetic anomalies were qualitatively and quantitatively interpreted using the 2.5D gravity-magnetic joint inversion technique. The interpretation results provide a basis for the location and depth of the borehole to be placed, which revealed a 15.8 m-thick iron-rich ore deposit, suggesting remarkable prospecting effects. This study holds critical indicative significance for further exploration of skarn iron ore deposits in this area.

Airborne geophysical surveys, characteristic of being green, economical, efficient, and subjected to less influence by surface factors, serve as the one of most effective means of basic surveys and scientific research on the Qinghai-Tibet Plateau. This study reviewed the progress in the airborne geophysical surveys in the Qiangtang Basin of the Qinghai-Tibet Plateau in the past thirty years, systematically summarizing the progress and geological interpretation results of comprehensive airborne geophysical surveys in the basin. Furthermore, this study presented research progress and understanding of major basic geological issues of the basin, such as the basin's boundaries, central uplift zone, basement properties, deep structures, and cover characteristics, as well as the identification of favorable structural areas for oil and gas exploration. Finally, based on an analysis of the developmental trends of airborne geophysical surveys in the Qiangtang Basin, this study proposed opinions and suggestions for key research directions in the future.

The Tiaojishan Formation in northwestern Hebei Province is home to volcanic-sedimentary strata. Due to the lack of fossil organisms, insufficient isotopic dating data, and outdated dating methods, the formation epochs of these strata have been controversial. To accurately determine their formation epochs and examine their regional geotectonic setting, this study conducted a detailed field investigation of the lithologic assemblage of the Tiaojishan Formation in the Xuanhua Basin. Petrological, petrogeochemical, and high-precision isotopic dating studies were conducted on the trachytic volcanic rocks in the upper segment of the formation, obtaining the LA-ICP-MS-based zircon U-Pb isotopic ages, which were 161.1±1.2 Ma and 162.5±1.3Ma. As indicated by the petrological and petrogeochemical characteristics, these trachytic volcanic rocks belong to the shoshonite series, exhibiting enriched light rare earth elements, relatively enriched elements including Rb, K, Th, Ce, Zr, and Hf, and relatively depleted Ba, Nb, Sr, P, and Ti. According to the comparison of principal parameters and graphical discrimination, the magma originated primarily from the melting of continental crustal materials and formed in the tectonic setting of compressional continental margin volcanic arcs. The comprehensive research reveals that the Tiaojishan Formation formed primarily during the Middle Jurassic and continued to the Late Jurassic, and the trachytic volcanic rocks in its upper segment formed in the intraplate compressional tectonic setting. The results of this study provide new data for the division and correlation of Mesozoic volcanic-sedimentary strata and the analysis of their formation environment.

Low-density geochemical mapping is characterized by high order streams to be sampled and a large coverage area and can be used to effectively trace geochemical blocks with high metal contents. Based on the 1:1,000,000 low-density geochemical mapping data of Ethiopia and the processing of the testing data of Cu in the original stream sediments, this study calculated the anomaly threshold of Cu at 37×10^-6 through iterative deletion. Then, this study delineated three geochemical blocks and two regional anomalies with 37×10^-6, 42×10^-6, 47×10^-6, 52×10^-6, 59×10^-6, and 66×10^-6 as grading intervals. It calculated the mineralization coefficient of Cu ore bodies in the study area at 0.055% by referencing the known reserves of Cu deposits in geochemical blocks with a high level of copper exploration in the same metallogenic belt. Moreover, this study estimated the Cu resources in the study area at 2,600,000 t based on a rock mass thickness of 1,000 m. By combining the analysis of metallogenic geological conditions, this study determined that the zones where geochemical blocks nos. 2, 3, and 4 are located can be considered key metallogenic prospect areas for further detailed exploration.

The performance of ray-based tomography inversion is affected by many factors,such as initial model error and low-velocity interlayer.The conventional tomography method based on first-arrival wave travel time,which constrains or smooths models,destroys the relative relationship between model parameters and rays and affects the inversion stability.By testing the performance of first-arrival wave travel time-based tomography inversion under different initial models,this study proposed a first-arrival wave travel time-based tomography inversion method with surface wave information as constraints.The process of this method is as follows:(1)Given that surface waves feature high energy and frequency dispersion in seismic data,the surface-wave frequency dispersion curves are obtained through the multi-channel analysis of surface waves;(2)Using the damped least squares method,the shallow-surface shear wave (S-wave) velocities are determined through inversion;(3)With the S-wave velocity structure as the constraint,the initial compressional wave (P-wave) model is established,and accordingly,the first-arrival wave travel time-based tomography inversion that considers regularization is achieved.This method improves the accuracy and stability of shallow structure inversion by fully utilizing the surface wave information in seismic data to counteract the inherent defects of tomography inversion.The effectiveness of this method has been verified using actual data.

The transient electromagnetic method (TEM) commonly uses the all-time apparent resistivity parameter for interpretation, which involves complex formulas and time-consuming iterative processes. Based on the characteristics of TEM data, this study employed the artificial neural network (ANN) for TEM pseudo-resistivity imaging. First, this study designed a multi-hidden-layer BP neural network and calculated a response amplitude through TEM analysis. The response amplitude, as the mapping parameter of pseudo resistivity, was used for network training. Then new data outside the training set were used to test the trained network. A homogeneous half-space and one-dimensional layered model was built to verify the correctness and adaptability of the neural network. The imaging of the three-dimensional geoelectric model was performed. As revealed by the results, the pseudo resistivity calculated based on the neural network can reflect the target anomalies of the geoelectric model, with highly accurate network imaging results. Finally, the measured data were processed using the neural network algorithm, further indicating that the neural network-based imaging can serve as a basis for data interpretation. This study verified the feasibility of the ANN in TEM imaging, thus providing a new approach for TEM imaging.

Distributed acoustic sensing (DAS) technology, one of the most advanced sound field detection technologies, can achieve distributed, long-distance, and high-precision real-time detection of the ambient vibration and sound field information interacting with optical fiber. The optical fiber exploration system of the DAS technology solves the problems of high cost and deployment difficulty of conventional geophones in complex geological environments. In recent years, the DAS technology has experienced rapid development, especially in monitoring application scenarios that require long-term and large-scale deployment. However, its systematic understanding is insufficient due to divergent research results. To further understand the research advances of the DAS technology in geophysical exploration for more effective subsequent research, this study systematically classified and summarized the development history of the DAS technology and its recent research results in geophysical exploration based on the oil and gas, marine, and environmental engineering application scenarios through literature research. This study focused on the development process of the DAS technology in different directions, the research advances in data processing, and relevant literature with specific results. Finally, this study generalized the development trend and urgent problems of the DAS acquisition system, analyzing the DAS development prospect.

Since individual geophysical exploration methods suffer the multiplicity of solutions, comprehensive geophysical prospecting has been extensively applied in deep ore prospecting presently. This study conducted the geological exploration of the Duhu porphyry copper deposit in Xinxing County using multiple geophysical methods such as high-precision magnetic survey and controlled source audio-frequency magnetotellurics (CSAMT). It was inferred that the CSAMT-derived medium-low resistivity anomalies and the low-gentle anomalies derived from the high-precision magnetic survey serve as significant prospecting indicators. Satisfactory results were achieved in follow-up verification of the anomalies. Specifically, copper, molybdenum, silver, and gold mineralized bodies with a cumulative thickness of 178.2 m were identified in a 1 000 m deep borehole, with the highest copper grade of 1.45%. The application of comprehensive geophysical prospecting holds great significance in guiding the exploration of porphyry copper deposits in western Guangdong.

Gravity survey is an important part of geodetic mapping and geophysical exploration. Different gravity specifications adopt different methods to calculate the zero drift rates in two-way reciprocal observations (i,j,j',i'). The repeated measurements of point j are subjected to static drift in the geological survey specification. By contrast, the zero drift rate is calculated through regression analysis in the petroleum specification. Different processing methods for the repeated measurement interval of point j yield significantly different zero drift rates based on the two specifications. Through theoretical formula analysis and comparison of measured data, this study analyzed the influence of g'_j- g_j {}_{} and t'_j-t_j in the repeated measurements of point j on the calculation of the zero drift rate. By comparing with the regression analysis, this study illustrated the importance of static drift in reducing the uncertainties caused by instrument performance and measurement environment. Moreover, this study proposed the applicable conditions of the two specifications.

It is difficult to explore the overlapping double-layer waterlogged goafs using the transient electromagnetic method. The reason is that upper waterlogged goafs will hinder the propagation of the electromagnetic field, thus prolonging the observation of the lower waterlogged goafs and reducing the signal-to-noise ratio. Besides, the burial depths and layer spacings of double-layer waterlogged goafs affect the signal-to-noise ratio and the observation time of transient electromagnetic signals. By building a double-layer waterlogged goaf model based on the Majiayan coal mine in Shanxi, this study analyzed the electromagnetic field propagation under layer spacings of 25 m, 50 m, 75 m, and 100 m,and calculated the observation time of waterlogged goafs with different layer spacings. Furthermore, it quantitatively characterized the differences between induced voltages in the double-layer waterlogged goafs with different layer spacings using root mean square errors. Additionally, this study proposed the identification criteria for explorable lower waterlogged goafs based on the record errors and noise levels during the observation. The results of physical simulation experiments are as follows: The differences between the induced voltages of double-layered waterlogged goafs with different layer spacings occur mainly in the late stage; the differences between induced voltages gradually decrease as the layer spacing and the burial depth of upper waterlogged goafs increases; the difference between induced voltages is close to the noise level when the layer spacing is greater than 75 m. The actual detection of the double layer waterlogged goaf with a spacing of 75 meters was conducted in Majiayan Coal Mine, and the results showed that the lower waterlogged goaf was not effectively identified.Therefore, It is difficult to effectively explore the lower waterlogged goafs when the layer spacing is greater than 75 m.

Seismic data acquired from mature industrial areas frequently contain a large amount of local strong noise with high amplitude due to the continuous operation of production equipment.However,such local strong noise can be hardly suppressed using conventional denoising methods.This study integrated dilated convolution(DC) and U-Net into a DC-UNet network for suppressing local strong noise.For the circular DC blocks at the front end of the DC-Unet network,a circularly expanded DC kernel was used to extract the features of strong noise at different scales,with the receptive field being expanded.Meanwhile,an encoder was used at the back end of the network to extract the features of strong noise and restore the details of strong noise.Subsequently,the DC-UNet network was employed to perform a nonlinear mapping from noisy data to noise.On this basis,strong noise was suppressed by subtracting the learned strong noise from the noisy data.As indicated by the experimental results of synthetic and real data obtained from the training using the PyTorch framework in the GPU environment,the DC-UNet network can effectively suppress the local strong noise and improve the signal-to-noise ratio compared with DnCNN,U-Net,and PCA-UNet networks.

Xiong'an New Area has great potential for geothermal resources. Carbonate rocks are favorable reservoirs for deep geothermal resources in this area. The integrated geophysical exploration technology is an effective way to ascertain the deep structures and the characteristic stratigraphic structure of carbonate geothermal reservoirs. Aiming at the exploration target of deep carbonate geothermal reservoirs, this study put forward a surface-line-point hierarchical and progressive geophysical exploration model. Using the high-precision gravity and aeromagnetic data, this model first investigated the distribution range of carbonate rocks, the thickness of carbonate strata, the distribution of deep-seated faults, and the fluctuation of bedrocks. Then, it analyzed the low-resistivity anomalies of geothermal reservoir strata using the magnetotelluric method. Finally, this model finely characterized the geothermal reservoir strata using two-dimensional seismic profiles and analyzed the velocity structure and regional structural characteristics of anomaly zones in the geothermal field. Based on the exploration precision and reliability of gravity, aeromagnetic, magnetotelluric, and seismic geophysical methods in the geothermal resource exploration of different stages, as well as other factors such as construction cost and efficiency, this study analyzed the economic applicability of geophysical methods in the exploration of deep karst geothermal reservoirs and suggested that the carbonate geothermal resources should be explored using the geophysical exploration technology combination of gravity, magnetic, and magnetotelluric methods.

Xiong’an New Area is one of the areas with the best conditions for the development and utilization of geothermal resources in sedimentary basins of eastern China. Wumishan Formation of Jixian System is the main thermal reservoir, with buried depth of the top boundary varies greatly, which istypical carbonate type thermal reservoir. In order to find out the stratigraphic structure and concealed faults within 6 km of the starting area of Xiong’an New Area, focus on finding out the spatial distribution and structural characteristics of the main deep thermal reservoirs, and predict the favorable geothermal areas, the time-frequency electromagnetic method is applied to the geothermal resources exploration of Xiong’an New Area. Through high-standard data acquisition, refined data processing and resistivity-constrained inversion, the pseudo-seismic imaging technology and borehole data are used to calibrate the resistivity horizon, which effectively improves the reliability of interpretation results. In this work, 8concealed faults were identified, mainly NNE strike normal faults, which controlled the salient-depression pattern and stratigraphic distribution of the study area. The study area is divided into Rongcheng salient, middle sub-depression, Niutuozhen salient and Baxian depression from west to east. The thermal reservoirs of Wumishan Formation is mainly distributed in Niutuozhen salient, Rongcheng salient and middle sub-depression, with buried depth of top interface is about 600~2600m. Vsing the inversion results of the time-frequency electromagnetic method, the three-dimensional geological model of the study area is constructed, and the Niutuozhen salient area is considered to be the optimal geothermal favorable area from the aspects of heat source, channel, reservoir, caprock and fluid, especially near the footwall of F₄.

The grounded-source transient electromagnetic (TEM) method, which enjoys the advantages of high topographic adaptability and large exploration depths, is suitable for deep resource exploration in mountainous areas. However, the TEM responses can be distorted due to topographic effects, causing great difficulties in data interpretation. This study investigated the influence patterns and correction method of topographic effects based on the three-dimensional unstructured time-domain finite element method. First, this study conducted the fine-scale description of mountains using unstructured tetrahedral grids, calculated the response of various topographic models, and analyzed the influence patterns of topographic effects. Then, it proposed a correction method for topographic effects based on the principle of the linear superposition principle of electromagnetic fields, established a geoelectric model of mountains according to the actual elevation data, and calculated the model responses through three-dimensional forward modeling. Subsequently, this study determined the topographic responses by subtracting the flat Earth model responses from the mountain model responses and then obtained the corrected TEM responses by removing the topographic responses from the total responses. The results are as follows: (1) The influence of mountains on the TEM responses is concentrated in the early stage and weakens gradually with time; (2) The topographic effects are concentrated near mountains, and their intensity depends on the distance of survey points from mountain peaks; (3) The influencing range and response amplitude of mountains are proportional to the mountain scale. In other words, a larger mountain scale corresponds to a larger influencing range and higher response amplitude; (4) Mountains with relatively high resistivity show more intense topographic influence. As shown by multiple models with simple and regular topographic anomalies, the corrected TEM responses, which match well with the responses from direct forward modeling, can effectively eliminate topographic effects to a certain extent. The research on the influence and correction method of topographic effects can be used as a reference for the processing and interpretation of TEM data of areas with complex terrain.

To ensure the high quality of data acquisition in high-precision ground magnetic surveys in the field, it is an important link to ascertain the instrument performance in a period using instrument calibration points. However, some workers fail to properly use instrument calibration points and thus fail to monitor the operation state of instruments on a day through the early and late calibration measurement of the instruments, affecting the acquisition of high-precision magnetic survey data. Based on the analysis of specific examples, this study proposed that the checking of early and late calibration measurement results of instruments should consider both the differences between these results, as well as the differences between these results and the actual magnetic field values of instrument correction points. The purpose is to ensure the high-quality field acquisition of high-precision magnetic survey data.

This study proposed a fracture prediction method based on unsupervised clustering of fracture sensitivity attributes to accurately characterize the distribution of fractures with different attitudes in tight sandstones of the Upper Triassic Xujiahe formation in the Yuanba area, northeastern Sichuan. First, the sensitivity of fractures was extracted and selected based on the optimized post-stack seismic data. Then, the convolutional neural network, a deep learning algorithm, was used to learn global massive fault and fracture databases of various types, obtaining the intensities, dip angles, and azimuths of fractures. In combination with high-precision-guided curvature attributes, an unsupervised clustering algorithm based on the Bayesian probability model was used to predict the intensities of fractures with different dip angles through dimensionality reduction using principal component analysis (PCA). The prediction results are highly consistent with both the fracture interpretation results from imaging logs and the geological results. The results of this study show that the third member of the Xujiahe Formation has more developed fractures than the second member. Fractures in the third member include both the fault-induced fractures distributed near the faults in the southeast flank of the Jiulongshan anticline and the fold-induced fractures in the areas with large formation flexures. By contrast, only fault-induced fractures near the faults occur in the second member.

The eastern Kunlun metallogenic belt, as a significant metal metallogenic belt in China, hosts extensive orogenic gold deposits and large-scale Kunlunhe, Gouli, and Wulonggou gold concentration areas. The Dachaigou gold deposit is a large-scale gold deposit newly discovered in the Wulonggou gold field in recent years. Despite its high metallogenic potential, the western extension of its ore belt has not been defined. Hence, this study conducted induced polarization (IP) sounding and wide-field electromagnetic sounding in the deposit. The results show that the known ore belt is situated in the regional gravity anomaly gradient zone, the transition zone of positive and negative weak magnetic anomalies, the edge of IP anomalies, or the electrical gradient zone. The development zone of the regional tectonic belt resides in the large-scale IP anomaly section. The regional tectonic belt is characterized by a wide range of low-resistivity anomaly zones. The IV and III alteration zones of the known ore belt are located in the opening position of the low-resistivity anomaly zone and the shallow electrical anomaly gradient zone, respectively. Based on the above understanding and the electromagnetic anomaly change patterns of several parallel profiles in the western extension segment, it was inferred that the regional ore-controlling structure extends steadily in the W-NWW direction, forming a favorable prospecting space in the western extension segment of the deposit. The results of deep geophysical exploration in the Dachaigou deposit indicate that geophysical methods manifest significant advantages in deep geological prospecting research, providing successful experience for deep prospecting in the eastern Kunlun gold deposit area.

To explore the critical factors influencing the results of microtremor-based exploration at a depth of thousands of meters,this study conducted experiments using triangular arrays based on spatial autocorrelation(SPAC) and extended SPAC(ESPAC).Focusing on factors such as array size,acquisition unit frequency,and acquisition duration,this study explored the frequency band ranges corresponding to different array sizes,the arrangement of arrays in kilometer-depth exploration for obtaining both deep and shallow data,and the improvement in deep resolution.Based on the analysis and discussion results,this study established a parameter-setting system to improve the accuracy of exploration at a depth of thousands of meters.

The Mogou fluorite deposit of Fangcheng County, residing in the Neoproterozoic Meiyaogou Formation, is a crucial part of the fluorite metallogenic belt in southern Henan Province. Based on the analysis of the geological characteristics of the Mogou fluorite deposit, this study conducted petrographic and petrogeochemical studies to define the source of ore-forming materials and probe into the genetic mechanism of the fluorite deposit. The results show that the fluorite ore body is veined or lenticular, with mineralization-related alterations composed of silicification, fluorite and sericite alterations. Its ore structures are primarily massive, followed by banded, striped, and brecciated types. The chondrite-normalized rare earth element distribution pattern of the fluorite ore body is similar to that of the Meiyaogou Formation marbles and the Yanshanian porphyritic plagiogranites, suggesting a genetic relationship. The possible genetic mechanism is that the F-enriched ore-bearing hydrothermal liquids intruded along the interbedded fault and reacted with the surrounding rocks to form the fluorite ore /mineralized body. The deposit belongs to the epithermal filling type along the interbedded structure.

Due to the lack of effective deep exploration techniques, it has become a key and difficult task to expand the scale of prospecting by obtaining the mineralization information in the deep and peripheral areas of the proven typical deposits of the Huijiabao anticline in southwestern Guizhou. To truly reflect the characteristics of deep geological anomalies, this study investigated the deep ore-bearing Longtan Formation of the Zhexiang gold deposit in the eastern Huijiabao anticline. The information on deep Au anomalies was extracted through geochemical exploration in boreholes of any azimuth according to the profiles, longitudinal profiles, and deep 3D planes of survey lines. Then this study summarized the distribution patterns of geochemical anomalies and evaluated the deep metallogenic potential of the Zhexiang gold deposit. The geochemical data obtained from deep boreholes were processed using the iterative histogram method and the sample length weighted average grade method, respectively. The results show that the Longtan Formation in the mining area had Au anomaly background values of (0.04~0.12)×10^-6 and an anomaly threshold of about 0.24×10^-6. As revealed by the anomaly contour maps generated from the processed data using the above two methods, the distribution characteristics of Au anomalies are in high agreement with those of deep ore bodies, and the distribution ranges of high Au anomalies in all profiles are highly consistent with the morphologies of the proved ore-bearing zones (bodies). Moreover, the large-scale and unclosed high Au anomalies follow the dip direction of fault F₂₀ in the first member of the Longtan Formation. The 3D distribution of Au anomalies shows that the nearly EW-directed high anomaly zone in the central mining area is consistent with the axial region of the Huijiapu anticline. As indicated by the analysis of the multi-dimensional geochemical anomalies and the geological characteristics of the study area, fault F₂₀ is the main ore transmitting and controlling fault in the mining area, and the northern and eastern deep parts around the mining area have great metallogenic and prospecting potential. In addition, two prospecting targets were delineated, needing further engineering verification. This study plays an important demonstration role in guiding the exploration of other mining areas in the Huijiabao anticline. The feasible study methods can be referenced for the study of geochemical anomalies in other mining areas.

To improve the detection level of municipal roads for rapid and effective municipal road collapse warning and rapid search for municipal pipe network distribution, the 22^nd Research Institute of China Electronics Technology Group Corporation designed and developed real-time three-dimensional ground-penetrating radar (3D GPR). With the real-time 3D GPR, which is based on the architecture of field programmable gate array and digital signal processor (FPGA&DSP), the institute achieved the design and implementation of several key technologies for the multi-channel high-speed acquisition system, enriching the road detection techniques and methods. The real-time 3D GPR enables high-speed acquisition of ten-channel radar data using the horizontally polarized antennae equipped with five transmitters and six receivers. The channels can be switched using high-speed switches, which operate in an interactive interpolation manner. The 3D GPR allows for up to 32 channels and detection speeds of above 60 km/h (channel interval: 2 cm). This is attributed to the switching of the antenna array using switches. The optimum antenna polarization design was verified by the comparison of experimental data from cavity- and metal-plate-based experimental sites. As a result, the optimal antenna scheme was determined. The measured results show that, compared to general LTD-2600 radar, the real-time 3D GPR boasts a higher acquisition speed and higher performance in terms of amplitude and phase, conducive to the improvement of road disaster detection technologies. Therefore, there is high market demand for the real-time 3D GPR.

Located in the loess tableland area,the eastern Ordos Basin features criss-cross ravines and gullies and great topographic fluctuations.Its complex surface conditions and near-surface structures caused significantly different amplitudes,frequencies,and phases in the original seismic data of different regions.These differences are unrelated to the subsurface geological information and can easily lead to misinterpretation.Therefore,there is an urgent need to develop a processing method for seismic data consistency according to the characteristics of seismic data in the study area,aiming to improve the fidelity of seismic data.Based on the summary of the current consistency processing methods,this study built a new consistency processing flow for seismic data under complex surface conditions.By innovating surface consistent amplitude compensation,this study employed the dual-domain near-surface Q-absorption compensation for the first time to broaden the frequency band and improve the amplitude and frequency consistency of seismic data collected across gullies and tablelands.In addition,this study applied the surface consistent deconvolution of Yu's wavelet to attenuate the low-frequency noise and broaden the effective signal frequency band,improving the signal-to-noise ratio,resolution,and wavelet consistency of seismic data.The actual application verifies that the method used in this study is feasible and effective and has a high popularization value.

To avoid aggregate interference on the detection target in concrete quality detection by Rayleigh waves, the abnormal signals from the aggregate and the target should be distinguished for high detection reliability. Hence, this study explored the propagation characteristics of Rayleigh waves in aggregate concrete. By building the mesostructural random aggregate models and employing the high-order staggered-grid finite difference scheme, this study examined the effects of the randomness, shape, size, and content of aggregate on the Rayleigh wave field and dispersion curves, thus obtaining the Rayleigh wave scattering, energy attenuation, distortion, and dispersion curve characteristics under different aggregate parameters. Based on the forward modeling data of different aggregate concrete models, this study quantitatively analyzed the sizes and influence ranges of aggregate-induced anomalies of the Rayleigh wave field and dispersion curves. The results are as follows. The randomness and content of aggregate could affect the energy attenuation of direct Rayleigh waves. The Rayleigh wave field was slightly influenced by the aggregate shape but significantly impacted by the aggregate size. When the aggregate size exceeded half of the dominant wavelength, Rayleigh waves would produce strong scattering, distorting its waveforms. In contrast, these aggregate parameters caused no anomalies in the dispersion curves. Therefore, data analysis using dispersion curves can avoid aggregate interference on the target in concrete quality detection by Rayleigh waves.

This study aims to ascertain the heavy metal pollution of groundwater in a typical chemical industry park in northwestern Shandong. According to the investigation and evaluation requirements for groundwater environments in chemical industry parks, this study collected 10 groundwater samples to analyze the concentrations and spatial distributions of eight heavy metals, namely Fe, Cu, Zn, Al, Cd and Pb. Furthermore, it conducted the health risk assessment of groundwater for different populations. The results indicate that: (1) except Mn and Al, all heavy metals in the groundwater of the chemical industry park showed concentrations lower than the class III water quality standard stated in the Standard for Groundwater Quality (GB/T 14848—2017), with a comprehensive pollution index ranging from 0.37 to 0.78. The monitoring points for heavy metal elements are all pollution-free; (2) the overall spatial distributions of heavy metals are roughly consistent with the distribution areas of key enterprises in the chemical industry park. This consistency indicates that external factors such as the production activities of enterprises affect the spatial distributions of heavy metals in groundwater; (3) the health risks are lower for children than for adults and lower via skin contact than via drinking. The non-carcinogenic risks under different exposure routes are acceptable for different populations. Carcinogens As and Cd have slight carcinogenic risks for adults. Local authorities should strengthen the supervision of the production activities of enterprises and thoroughly inspect enterprises with suspected pollution. Only in this way can we effectively control and then gradually reduce the risks to human health caused by groundwater pollution.

Velocity is a key parameter determining the migration imaging resolution of ground penetrating radars (GPR).The method combining minimum image entropy and migration usually estimates the medium velocity by calculating the entropy curves using the overall migration profile as a fixed window.Therefore,such a method is not applicable to non-uniformly distributed media.Moreover,for this method,a too-high or too-low test velocity will make the convergence position of hyperbolic diffracted waves go beyond the fixed window,thus reducing the estimation accuracy.This study proposed a minimum entropy method based on a velocity-controlled moving window,in which the calculation window in the migration profile is accurately controlled by the test velocity.Then,this method was combined with inverse time migration to estimate the optimal migration velocity.By automatically adjusting the position of the calculation window using the trial velocity,this method keeps the convergence position of hyperbolic diffracted waves at the center of the calculation window.In this manner,stable and accurate entropy curves can be obtained.By comparing the calculation results with those of the minimum entropy method based on a fixed window,this study verified the correctness and effectiveness of the minimum entropy method based on a velocity-controlled moving window for a typical hyperbolic diffracted wave.As revealed by numerical experiments and the tests of measured data,compared with the minimum entropy method based on a fixed window,the minimum entropy method based on a velocity-controlled moving window can keep the convergence position of hyperbolic diffracted waves accurately at the center of the calculation window,yielding more stable entropy curves,lower computational complexity,higher estimation accuracy of the migration velocity,and better imaging performance of reverse time migration.

This study aims to systematically investigate the distribution and over-limit elements of areas with high geological background values of heavy metals in soils in Yunnan Province. GIS spatial analysis was conducted based on the heavy metal content data from a province-wide 1∶200,000 stream sediment survey and the regional geological map. The analysis results were validated using the data of heavy metals in soils in Kunming, Yuxi, Zhaotong, and other regions. A total of 61 geological units were identified, with heavy metal content in soils exceeding the screening values of agricultural land, accounting for 21.09% of the total land area of Yunnan. The cultivated land in high geological background areas covers an area of approximately 2.844 1 million hectares, accounting for 7.22% of the total land area of Yunnan. The lithologies that cause over-limit heavy metals in soils primarily comprise carbonate rocks, mafic-ultramafic volcanic rocks, intermediate mafic intrusive rocks, coal-bearing clastic rocks, and clastic rocks with mafic components. The over-limit heavy metal elements in high geological background areas are dominated by Cu, Cr, Ni, and Cd. In contrast, As manifests an over-limit risk mainly in carbonate rock formations, Pb and Zn only exhibit an over-limit risk in individual strata, and Hg almost shows no over-limit risk.

Significant aeromagnetic and gravity anomalies were found in the ultra-deep coverage zone of the Pandian area in the northwestern Shandong Province. Through systematic verification of gravity and magnetic anomalies, a breakthrough in prospecting for deep skarn iron deposits has been achieved through several boreholes around the Pandian gravity-magnetic anomaly zone. However, the gravity-magnetic anomaly zone spreads extensively and its characteristics cannot directly indicate the prospecting target, thus some boreholes failed to find ores or revealed poor ore-finding conditions. Hence, to accurately delineate the ore-forming location in the gravity-magnetic anomaly zone and achieve a further breakthrough in ore prospecting, this study conducted the wide-field electromagnetic (WFEM) sounding in the verified Pandian gravity-magnetic anomaly zone. Combining the drilling verification, this study delineated the deep ore-forming interval of skarn iron deposits in the deep coverage zone and the favorable structural plane for ore-forming, suggesting the deep fine-scale prospecting targets.

The Niuxingba plumbum-zinc-gold-silver deposit (the Niuxingba deposit) occurs in the Yinkeng ore field at the intersection of the EW-trending Nanling metallogenic belt and the NNE-striking Wuyishan metallogenic belt. It is a typical hydrothermal vein deposit under the strict control of fault structures, with ore bodies occurring as vein groups. To predict the deep prospecting potential of the No. V31 major ore belt (body) of the Niuxingba deposit, this study systematically investigated the axial (vertical), transverse, and longitudinal geochemical characteristics and zoning of tectonic primary halos of the ore belt (body) based on the geological characteristics of the deposit. Then, this study determined the axial zoning sequence and the geological-geochemical model of tectonic primary halos. It was confirmed that this model was correct, reasonable, and reliable using the multivariate statistical method. The results of the geological-geochemical model show that: ① the front halo of the Niuxingba deposit has F as the indicator element; ② the near-ore halo has two types of indicator elements: the gold mineralization-related As-Au element association in the middle and upper parts of the ore body and the Pb-Zn-Ag mineralization-related Hg-Zn-Pb-Ag-Cu-Bi-Sb element association in the middle part of the ore body; ③ the rear halo has the Mo-Cr-Co-Ni element association as indicator elements; ④ the anomaly centers of the front, near-ore, and rear halos are about 200~300 m apart axially, suggesting significant coexistence of front and rear halos and the presence of inflection points of the near-ore halo broken line. Furthermore, the transformation of low-moderate-temperature Pb-Zn-Ag-Au mineralization into moderate-high-temperature Cu-Zn mineralization occurs eastward. These characteristics indicate that the No. V31 main ore belt (body) has a pitch direction of east and a large extension toward the deep part (to an elevation of -200 m at least), implying that blind ore bodies may occur. As verified by deep drilling in the eastern concealed area (the lowest ore-controlled elevation of the No. 367 survey line: -165 m), the eastern deep part of the No. V31 ore belt (body) has a high potential for increasing reserves.

Given large data volumes and low post-processing efficiency of full matrix capture-total focusing imaging,this study proposed a planar array-total focusing method(PATFM) for the imaging of planar array data.First,the wave front time of the downgoing planar array was calculated using the eikonal equation based on the characteristics of both the total focusing imaging algorithm and the planar array wave field.Then,the total focusing imaging formula was improved using the upgoing and downgoing ultrasonic propagation time based on the delay superposition principle.Finally,focusing imaging was performed on a wide range of imaging points below the planar array aiming at the derived directivity and diffusion correction coefficient of the planar array.Through Field II simulation,the PATFM was compared with three imaging methods,including phase-controlled scanning imaging,full matrix capture-based total focusing imaging,and plane wave capture-based total focusing imaging.The results show that the PATFM can be used for large-range focusing imaging of single planar array data,greatly improving the computational efficiency while obtaining corresponding accuracy.Therefore,this study provides a feasible technical means for 3D imaging of array acoustic waves.

Investigating the concentrations of natural radionuclides in soils is basic for soil pollution prevention and control. First, this study determined the background values and anomaly thresholds of radionuclides ²³⁸U, ²²⁶Ra, ²³²Th, and ⁴⁰K in soils in Guangdong Province using the geographic information system (GIS) and three methods, namely conventional statistics, exploratory data analysis (EDA), and median absolute deviation (MAD). Then, this study assessed the radioactivity levels of the soil environment using four indices: the γ radiation dose rate, the equivalent radium specific activity, the external exposure index, and the annual effective dose rate. Finally, this study assessed the health risks of the soil radioactivity using the geoaccumulation index method, the single factor index method, and the Nemerow composite index method. The results indicate that: (1) Radionuclides ²³⁸U, ²²⁶Ra, ²³²Th, and ⁴⁰K in soils of Guangdong Province have average concentrations of 79.4, 102.95, 74.59, and 541.25 (Bq·kg^-1·dry), respectively; (2) The four radionuclides have background values of 34.27、46.75、32.87 and 271.79 (Bq·kg^-1·dry), and the anomaly values of 111.35、141.805、105.12、946.26 (Bq·kg^-1·dry); (3) The radioactivity level indices for soil environment, namely D, Ra_eq, Hex, and AEDR, show averages of 121.44 ± 67.08 nGy/h, 263.48 ± 147.72 Bq/kg, 0.08~2.3, and 148.93 ± 82.26 μSv/a, respectively; (4) The radionuclide geoaccumulation index suggests clean soils, the single factor index indicates slight pollution, and the Nemerow composite index of 0.22~5.52 (average: 1.74 ± 0.97) indicates mild pollution. These three indices, which suggest roughly consistent results for the health risk assessment of radionuclides in soils, can serve as a supplement to the health risk assessment method for soil radioactivity.

Since the periphery of the Xishan rare earth deposit in Weishan County, western Shandong is mostly covered by the Quaternary strata, single geophysical exploration methods are ineffective in searching for rare earth elements in this area. To establish a geological-geophysical prospecting model for deep rare earth prospecting in this area, a combination of exploration techniques and methods suitable for the deep exploration of rare earth deposits in this area were selected through multiple geophysical exploration technology experiments on the known geological profiles. The distribution range of the underground Mesozoic Xishan alkaline complex was roughly determined through gravity survey and high-precision magnetic survey. The ore body outcrop was delineated through ground-based gamma spectrometry, and the deep ore body characteristics were revealed through drilling. Finally, this study developed an integrated geophysical exploration method including gravity-magnetic joint delineation of rock masses, radioactive positioning for ore body outcrops, and deep drilling. By using this exploration model, one super-large rare earth deposit was discovered in this area, providing a reference for the exploration of rare earth deposits in the surrounding area.

This study synthesized microtremor signals using the normal mode superposition.The dispersion curves of the synthesized microtremor signals were calculated using the spatial autocorrelation method,and the calculated results agreed well with theoretical dispersion curves.The simulation of microtremor signals involves many parameters,which are of great significance for microtremor exploration.As shown by results from numerical experiments of these parameters,the hypocentral distance and array size had a significant influence on the results.Moreover,high-order surface waves should be considered for complex strata,while small-size arrays tended to yield high-frequency dispersion information at quiet observation sites,thus improving the resolution of shallow strata.Therefore,it is necessary to consider the influences of hypocenter distribution,array size,and high-order surface waves on exploration results in microtremor exploration.

As an inevitable physical phenomenon in the application of frequency-domain electromagnetics, the static shift effect is generally suppressed or eliminated by correction. This study proposed a new approach of directly utilizing the static shift effect of natural electromagnetic methods to explore shallow electrical anomalies. The frequency selection method of telluric electricity field (FSM) is to study the variations in electrical properties of subsurface media by measuring several horizontal electric field components with different frequencies generated on the surface by the natural alternating electromagnetic field. In this study, the forward modeling of FSM data was conducted using the two-dimensional finite element method. The modeling results are shown as follows: (1) In the case of low-resistivity anomalies near the surface, the curves of horizontal electric field components along the survey line on the surface showed the same morphologies as the FSM-derived curves, with significant low-potential anomalies above the low-resistivity anomalies; (2) As the calculated frequencies increased, both the profile curves and the pseudosection map of electric field components exhibited a static shift effect, indicating that the FSM-derived anomalies were mainly caused by the static shift effect. Both the FSM application results and the drilling verification results showed that with the presence of groundwater, the FSM-derived profile curves and pseudosection map exhibited a significant static shift effect, which was consistent with the CSAMT exploration results. As indicated by theoretical and practical research, it is feasible to directly use the components of the telluric electricity field for the exploration of shallow electrical anomalies. Moreover, shallow geological exploration can be conducted by utilizing the static shift effect of the frequency domain electromagnetics.

To reduce the limitations and the multiplicity of solutions of a single geophysical inversion method, this study investigated the two-dimensional joint inversion based on the marine controlled-source electromagnetic (MCSEM) method and seismic full-waveform inversion. The MCSEM method employs the data-space Occam’s algorithm, while the seismic full-waveform inversion utilizes the gradient algorithm. By incorporating a cross-gradient function for the mutual coupling of the two types of physical property parameters, this study developed a two-dimensional joint inversion method, whose accuracy was verified using three different models. As indicated by the results, compared to a single inversion method, the MCSEM-based joint inversion yielded significantly improved inversion results, predominantly in terms of the morphology characterization of anomalous bodies, as well as the reconstruction of their structure and textures and their physical property values. Therefore, the full-waveform inversion can enhance the reliability of the MCSEM inversion results.

Huize County of Yunnan Province is situated in the karst area in southwestern China, where karst groundwater is its primary water source. To conquer local difficulties in drinking water, this study constructed a conceptual model of groundwater occurrence by fully investigating the hydrogeological conditions of the Huize area. Moreover, this study evaluated the applicability and optimal combination of geophysical methods based on the measurement results of petrophysical properties. According to the actual local needs, this study deployed a comprehensive profile combining high-density electrical resistivity tomography (HDERT) and audio magnetotellurics (AMT) in Tuogu Village, Huize County. The groundwater enrichment site was delineated relying on resistivity anomalies, effectively guiding the layout of boreholes. The boreholes achieved the maximum single-borehole water yield of 20.76 m³/d, thus effectively alleviating the local drinking water problem. The HDERT-AMT combined exploration method proves to be optimal for prospecting for groundwater in carbonate rock areas. HDERT can accurately characterize weathered layer thicknesses, bedrock boundaries, fissure evolutionary degrees, and water-bearing properties of strata, constraining groundwater recharge channels, thus counteracting AMT's defects for identification of near-surface stratigraphic structures. AMT can accurately reflect the spatial structures of fracture zones and the macrostructures of strata, limiting the boundary conditions (aquicludes) of water-bearing structures, thus making up for the defects of insufficient detection depths of HDERT in high-resistivity stratigraphic regions. HDERT and AMT, which are complementary to each other in terms of accuracy and depth, can be applied to identify and constrain the spatial occurrence conditions of groundwater migration, storage, and enrichment.

To evaluate the qualities of the marine magnetic data in the Earth magnetic anomaly grid EMAG2v3 and the marine gravity data in global gravimetric database V29, this study selected the magnetic and gravity data of the Aegir axial rift and its adjacent areas within a range of about 150 km from EMAG2v3 and V29, respectively to conduct comparative research. This study systematically collected the anomaly data of the study area from EMAG2v3 and V29 for comparison with the measured gravity and magnetic data of the study area. First, this study gridded and whitened the EMAG2v3 data, V29 data, and measured gravity and magnetic anomaly data to obtain the corresponding images. Then, this study analyzed the correlations between the EMAG2v3 data and the shipborne magnetic data and between the V29 data and the shipborne gravity data, obtaining the magnetic and gravity correlation diagrams and corresponding correlation coefficients. By comparing the correlation coefficients and differences between the two kinds of magnetic data and the two kinds of gravity data, this study conducted an integrated evaluation of magnetic the gravity data of the study area from EMAG2v3 and V29, respectively. As indicated by the results, EMAG2v3 incorporates many shipborne magnetic data, with the shipborne magnetic anomaly data showing higher quality than the data from the EMAG2v3 for areas with dense survey lines. The results also show that the shipborne gravity anomalies showed roughly the same variations as those from V29, indicating the same lateral resolution of the two types of anomaly data.

This study aims to ascertain the distribution patterns and influencing factors of fluorine in irrigation water and soils in the Gaomi area, Shandong Province, China. Through systematic sampling and testing, this study obtained data including the pH and fluoride concentration of irrigation water and the pH, organic matter, and fluorine concentration of soils. Based on these data, this study plotted the geochemical contour maps for statistical, correlation, and difference analyses. Furthermore, this study explored the influencing factors. The results are as follows: (1) The irrigation water in the study area has a fluoride concentration of 1.89 mg/L on average, which shows strong spatial variability. Zones with high fluoride concentrations in irrigation water are distributed in the northern low-flat alluvial plain, with the number of samples with fluoride concentrations greater than 2 mg/L accounting for 63.16%. There is a significant positive correlation between the fluoride concentration in the irrigation water and the pH (P < 0.01); (2) The soils have a fluoride concentration of 455×10^-6 on average, which shows an inhomogeneous spatial distribution. Zones with high or excess fluorine concentrations are distributed in the northern part of the study area. The fluorine concentrations of soils show significant positive correlations with the pH and organic matter of soils and the fluoride concentration of irrigation water (P < 0.01); (3) The fluoride concentrations in the irrigation water and soils are high in the distribution area of lime concretion black soil. The results of this study reveal the background characteristics and influencing factors of fluoride in the irrigation water and soils of the Gaomi area, providing a geochemical basis for the precise prevention and control of endemic fluorosis.

Ground-penetrating radar (GPR) data are often contaminated by random noise in the actual engineering inspection.The noise in data will reduce the signal-to-noise ratio and resolution of the data,adversely affecting the subsequent inversion and interpretation.Accordingly,this study proposed neural network-based denoising for GPR data.First,a multi-layer neural network model was constructed to integrate the data corrupted by white Gaussian noise into the noise-free data.Then,the corrupted data and their corresponding noise patches were built as training data.The weights of neurons in every layer of the model e updated using a back-propagation algorithm to minimize the model training loss.Finally,the two synthetic data and the measured radar data were input to the trained model,and the model's output was calculated using the noise characteristic weights acquired from the training.Compared with the curvelet transform,the numerical simulation test results verify the effectiveness and robustness of the method proposed in this study.Moreover,the proposed method can suppress the noise more thoroughly in areas with complex structures and weak amplitudes,and show effective signals more clearly.

Targeting the selenium resource in soils of Luoyang City, Henan Province, this study analyzed the distribution of total selenium concentration under different geological settings, soil types, land use types, acidity and alkalinity, and organic matter. Furthermore, it determined the bioavailability and speciation distribution of selenium in the soils. The results indicate that: (1) The soils in Luoyang City have total selenium concentrations of (0.03~5.67)×10^-6, with an average of 0.30×10^-6, and the area of soils with moderate and rich selenium accounts for 94.73%; (2) The average selenium concentration is the highest in soils of the Jixian Yunmengshan formation and the lowest in soils of the Cambrian Xinji Formation; (3) The average total selenium concentration is the highest in lime concretion black soil and the lowest in paddy soil; (4) Regarding land use types, the average total selenium concentration is relatively high in soils of paddy land and woodland; (5) Soils with high acidity and high organic matter concentrations are favorable for selenium enrichment; (6) In terms of soil types, the average concentration of bioavailable selenium is the highest in fluvo-aquic soils and the lowest in red clays; (7) The average concentration of bioavailable selenium is higher in alkaline soils than in neutral and acidic soils; (8) Selenium in soils primarily occurs in residue, strong organic bound, and humic acidic bound forms, with the highest average concentration of bioavailable selenium occurring in alkaline soils and red clays. The selenium resource in the soils in Luoyang City has high utilization potential, and it is feasible to conduct targeted development and planning of the selenium-rich industry according to the distribution of total selenium concentration in the study area.

The Tonglvshan ore field is the most important skarn copper-iron-gold polymetallic ore field in the southeastern Hubei ore concentration area. Since its mineralization is closely related to Tonglvshan stocks, ascertaining the deep distribution and morphologies of these stocks and their contact relationship with surrounding limestones and marbles holds great significance for deep prospecting prediction of the ore field. This study first conducted a detailed analysis of the gravity and magnetic anomaly combinations of the Tonglvshan pluton and stocks; limestones and marbles, as well as their xenoliths; and skarn ore bodies and mineralized bodies. Accordingly, this study inferred the boundary of the Tonglvshan pluton and delineated the xenoliths and mineralized bodies of limestones and marbles within the Tonglvshan pluton. Then, this study conducted the 3D geological modeling based on the results from the 3D physical property inversion and 2.5D and 3D gravity-magnetic human-computer interaction inversion, as well as the drilling and geological data. Finally, this study established a 3D geological model for the plutons and mineralized bodies at depths of 4000 m and below in the Tonglvshan ore field, providing a basis for target delineation and deep prospecting prediction.