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.

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.

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.

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.

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.

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.

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 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.

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.

 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.

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.

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.

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 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.

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.

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.

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.

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.

 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.

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.

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.

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, 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.

The authors made statistic calculation of stream sediment data using such means as traditional statistics, robust estimation，cumulative percentage of 85% and multifractal statistical method in Chabaqi area and compared the thresholds obtained by the four methods. The thresholds of the same elements obtained by the four methods are obviously different from each other. The thresholds calculated by the cumulative percentage statistics are lowest except for lead, and the thresholds calculated by the other three methods fail to show obvious regularity.The results are closely related to the distribution characteristics of the geochemical data and the statistical methods used. Comparing the anomaly characteristics with the geological characteristics, the authors have found that, although weak and gentle geochemical anomalies will be mostly detected by cumulative percentage of 85%, the larger anomaly areas will bring difficulties to anomaly inspection in the field. It is therefore thought that the thresholds obtained by multifractal statistical method are suitable for the study area. A comparative study of the thresholds has led the authors to believe that the reasonable method must be combined with the actual geological conditions and based on a thorough study of the geological background of the study area.