Airborne geophysical data often suffer from corrugations or line level errors, which can be adjusted or removed by using levelling procedure. Levelling is a critical step in airborne geophysical data processing and interpretation. Microlevelling routine can be applied to removing the remaining line level after tie line leveling. Based on practical aeromagnetic data levelling, the authors deal in detail with the principle of microlevelling and key steps in microlevelling procedure. For the purpose of microlevelling data, a directional high pass filter perpendicular to the flight line direction is first employed to produce a decorrugation noise grid. The noise grid is then extracted as new channel flight data. Amplitude limiting and low pass filtering can be applied to the noise channel so as to remove the residual geological signal and leave only the component of line level drift, which is then subtracted from the original data to produce the final microlevelled data. This paper has also discussed some key technical problems in airborne geophysical data leveling, especially in microlevelling. The discussion in this paper may be useful in practical data processing.

 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.

 It is difficult to exactly detect small objects in complex shallow subsurface. As the energy of backscattered signal from the target is low, backscattered signals from target and ground surface are overlapped. In addition, clutters created by ground surface and movement strongly depress the signaltoclutter ratio when shallowly buried small objects are detected by portable ground penetrating radar (GPR). An algorithm for detection of shallowly buried objects based on principal component analysis (PCA) is thus proposed in this paper. Via PCA decomposition, current Ascan data are projected onto the projecting direction of background data. A set test function is compared with adaptive threshold to decide if current Ascan data are from an object. Detection of shallowly buried small objects can be achieved in combination with background data dynamic updating. The data tested in sand, laterite, clay and lawn were processed, and the results show that shallowly buried objects can be detected using algorithm based on PCA.

In Shanghai area, borehole injection piles are used in large quantities for construction foundation in soft soil strata. The reinforcement cage for each injection pile consists of several parts linked together through welding. After the completion of the borehole injection pile construction, the inspection of the length of the reinforcement cage and the examination of welding quality constitute a very difficult problem. This paper recounts the technology and effects of applying high precision foundation pile diagnostic device and geological radar to the inspection of the length of reinforcement cage for the borehole injection pile.

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.

Based on studying the distribution of trace elements in the dioritic weathered crust of Zigui area, the authors have divided the trace elements into three types. The first type consists of inactive elements which include Sc, V, Cr, Co, Ni, Cu, Zn, Zr, Hf, Nb and Ta, U, Th, Mo, W. During the weathering process, they can be effectively preserved. Due to the influence of redox condition variation, U, Cr and Cu sometimes show local enrichment, without obvious regularity. The second type is composed of active elements such as Ca, Rb, Cs, Sr, Ba, Pb, Ga, Gd, and Tl, which are gradually leached from profiles along with the weathering. The third type is represented by REE, which are leached from the upper section downward, and enrichment is observed in the mid-lower section, with the leaching degree of HREE higher than that of LREE. Eu gradually changes from positive anomalies in parent rocks into negative anomalies on the surface due to the weathering leaching of plagioclase in parent rocks under the oxidation environment. Under the surface oxidation condition, Ce easily produces tetravalence oxides (cerianite) and shows obvious surface enrichment, and the variation of the redox conditions leads to the undulate variation of Ce.

:The double current supply, inline and combined pole-dipole/dipole-pole array for time domain induced polarization method is an unsymmetrical arrangement sounding technology composed of forward pole-dipole and reverse dipole-pole array. This paper shows the sounding principle with graphical representation. For collecting resistivity and chargeability raw data distribution features in the pseudosection with forward pole-dipole and reverse dipole-pole array, it is similar to MT in that MT raw data obtain "static state effect", and resistivity and chargeability raw data also obtain "partial topography " and "electromagnetism noise" effect and assumes "false belt anomaly" features. Forward pole-dipole and reverse dipole-pole array obtains resistivity and chargeability data and can respectively process two-dimension inversion, but only the combined pole-dipole/dipole-pole array can get the best geological effect, indicating that the strongpoint of the electrodes configuration array is that it can obtain massive collection data, strong induced polarization signal, larger penetration depth and highly surveying accuracy. With the electrical property parameters for resistivity and chargeability two-dimension inversion imaging techniques in the geoelectrical section, we can accurately determine the spatial distribution of the electric anomaly body and can also provide accurate geophysical basis for drilling verification of the electric anomaly.

As an important method in potential field data processing and interpretation, Euler deconvolution has aroused widespread attention among researchers because it can provide automatic or semi-automatic estimates of source locations and depths under the condition of less priori information. This paper described the principles of Euler deconvolution, summarized its improvement as well as related problems that geophysical experts both in China and abroad have been working on in recent decades, analyzed the choice of structural index and the stability of Euler solution and the others. The future research emphasis and developing direction of Euler deconvolution are pointed out in this paper, which are of important significance for practical use of Euler deconvolution in potential fields.

Dajingpo Airborne Geophysical Test Site is the first airborne geophysical(magnetic, gravity and electromagnetic)integrated test site in China, at the sponsor of 863 project. The selection of integrated test site is an important and difficulty work, because it need to consider of all the conditions for the test of magnetic, gravity and electromagnetic systems. Therefore, the selection of integrated test site need to consider the basic principles as follow: the interfere is small in the test site and around; there is diversity in topography; there is affluent in geology, and possess of iron points (or iron deposits) are best; magnetic, gravity and electromagnetic abnormities are obvious., and so on. According to the basic principles and methods of "Select a large region first, and the site wherein ", after the synthetical analysis and contrast in topography, geology and geophysical materials and check in field, Siziwangqi Test Site is filtered from all the county, then do more detailed work to Siziwangqi Test Site, Dajingpo Test Site is determined in the Dajingpo area at the NE corner of the Siziwangqi Test Site.

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.

It is convenient to calculate derivative in wave number space using the inbuilt FFT function in Matlab. This paper describes a matlab algorithm for derivative calculation of 1D gravity and magnetic anomaly in wave number space, gives code of program and deals with some programming skills for precision improvement. Through model test and data analysis, it is found that the precision of the FFT method is much higher than that of the Fourier Series method in vertical derivative calculation. In horizontal derivative calculation, however, the precisions of the two methods are almost the same. During kainite exploration in a certain basin, satisfactory result was obtained in the processing of gravity profile data using this method.

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.

Exemplified by the utilization of geophysical data in search for a concealed iron deposit and based on the deduction and interpretation of the gravitational, magnetic and electrical cross-shaped profile data in a certain area, this paper sums up the train of thought and working methods for reducing the ambiguity in the geophysical anomaly interpretation. It is held that various kinds of constraint data should be accurately and completely collected during the field work and data acquisition so as to reduce the ambiguity, and that the one-by-one judgment and exclusive method seems to be an effective way for this purpose.

Starting with algorithm designing, training image,and practical application,the authors analyzed multiple-point statistics research trends both in China and abroad.On such a basis,the core principles of the main four MPS algorithms,their applicable ranges and advantages as well as disadvantages were analyzed,so as to forecast the trend of reservoir modeling.Multiple-point geostatistics is an international forefront research tool in stochastic modeling;nevertheless,as the algorithm is not yet mature,it should be further improved.On the basis of previous researches,the existing problems of the multiple-point geostatistics that need to be modified and the direction of the processing of training images are proposed,such as suitable processing of training image,choice of similarity methods,choice of standardized data,smoothing,coupling among algorithms,and expansion of the simulation extent for fractured-vuggy carbonate reservoir.In order to improve the usage of multiple-point geostatistics for reservoir stochastic modeling,we should spare no efforts to increase the simulation area,improve simulation accuracy,expand the scope of simulation of reservoir types,save simulation time and improve simulation efficiency.

The Laowan gold ore belt is an important component part of the Tongbai-Dabieshan (northern slope) gold-silver metallogenic belt, and the Laowan gold deposit is one of the major deposits in this belt and also one of the important large-size metamorphic clastic rock gold deposits in China.Based on an analysis of geological,geophysical and geochemical characteristics of the Laowan gold deposit,the authors established the integrated ore-prospecting criteria and ore-prospecting model to guide the ore exploration and target-locating work.The ore-prospecting criteria and ore-prospecting model are of great reference value in further work.

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.

 According to electric differences between the bauxite layer and the wall rock, the authors applied such electric

techniques as intermediate gradient method, IP sounding method and highdensity resistivity method to the prospecting for concealed

bauxite deposits in western Henan. The application of the model combining geophysical exploration with drilling has reduced the

exploration cost and speeded up the exploration work in the ore district.

Usually, city road consists of the surface layer and the subbase layer, and each layer contains several layers of different media with different relative dielectric constants. So using one relative dielectric constant to analyze the GPR data of all depths in city road detection is not accurate. This paper deduced the calculation formula of the equivalent dielectric constant of multi-layer underground media. Then the authors took two-layer media model and three-layer media model as examples to calculate the equivalent dielectric constant and the relation between the depth and the travel time. At last, theoretical modeling data of two-layer media were used to check the formula, and the results prove that the calculation formula of equivalent dielectric constant of multi-layer underground media is much better for practical GPR data analysis when there are multi-layer underground media.

This paper discusses methods of estimating airborne gravity anomalies. Based on the Kalman smoothing techniques, the authors developed the software for estimating airborne gravity anomalies. Research shows that, if this software is used to estimate airborne gravity anomalies, the standard deviation of the discrepancy is within 1 mGal in comparison with the data processing results by the imported GT-1A system software.

Break points or negative values often appear in the data of charging rate during the performance of the multi

electrode induced polarization method in the field, which causes great difficulty in data interpretation. Based on an analysis of

the causes of break points and negative values, this paper gives prominence to the useful anomaly by using interpolation of cubic

spline to process the data. Threepoints slurring filtering is conducted to reduce random interference effectively.

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.

The present paper makes a brief description of the progress of several main electric methods which havebeen developed quite rspidly since 1980's.They include induced polarization method,frequency spectrum IPmethod,trandient electromagnetic method, controlled-source audio-frequency magnetotelluric method and GPR.

In this paper, a simple method for normalization of geochemical data is described. It can be used to reduce the differences between different analytical batches and normalize the data from different areas which have different background values. The method lies in using the ratio of the original analytical value to mean or median value in the batch multiplied by a reference value. The transfer value is given the name apparent content. Two examples are given in this paper to illustrate the details of adjusting the differences.

A soil gas radon survey was performed on a large scale to determine the distribution of radon in soil of Zhuhai City in Guangdong Province by means of a portable radon monitor of a semiconductor alpha spectroscopy. The survey sampled 469 sites covering an area of more than 100 km². The average of soil radon concentration in the soil depth of 0.6 m is 55.94 ± 58.54 kBq/m³ in Zhuhai urban area, whereas the concentration is 7.14±8.75, 37.64±25.92, and 151.25±196.23 kBq/m³ in the Quaternary sediments, the mixtures of sediments and weathered grain of granite, and the weathered granite in Doumen District, respectively. The high radon potential areas are located within biotitic granites and new industrial districts, as indicated by the strong correlation between the radioactivity level and geological lithology. The mean value of soil gas radon concentration in Zhuhai urban area (ZUA) is about ten times as high as that in Guangzhou, Quanzhou and Jinjing City. The results show that Zhuhai area has higher radon potential, and hence protective measures against radon should be taken into account.

The Yindongpo gold deposit is a large-size gold deposit hosted in metamorphosed clastic rocks. Based on a geochemical study of the Yindongpo ore deposit, this paper puts forward some useful indicators for further prospecting. They include complexity of the element association with an obvious Au, Ag, Pb, Zn, As, Cd anomaly, an apparent element concentration center with clear element zonation; w(As)/w(Mo)≥4, w(As)+w(Cd)+w(Au)/w(Ni)+w(Co)+w(Mo)≥10, K_k (enriched Clark value) ≥0.5, C_v (Variation coefficient)≥1, regular distribution of striped anomaly as well as stratigraphic and structural control. On the basis of these indicators, ore prognosis was conducted in the periphery of the Yindongpo deposit, covering the area between W4 and W26 exploration lines. The result has been verified by drilling, with new reserves of 11.945 t metals of gold obtained.

The problems encountered in the application of tilt derivative method to aeromagnetic data processing have been analyzed in this paper. Influenced by the shallow-source interference field, the shape of the aeromagnetic data anomaly processed by the tilt derivative method is complex. An improved method that combines the upward continuation with tilt derivative is therefore put forward in this paper. With practical examples, this paper deals with the problems that are worthy of attention in the application of tilt derivative method to aeromagnetic data processing. The improved method can fully manifest the advantages of the aeromagnetic data in the aspect of geological data interpretation.

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.

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.

This paper deals with three gravity satellites, namely CHAMP, GRACE and GOCE. Satellite gravity data are widely used in such fields as the recovery of the earth's gravity field, the land water storage, the study of the geoid change, the elastic thickness, the crustal thickness, the gravity variance before and after earthquake and the crustal movement. Lots of successful and valuable research results were obtained. It is considered that the satellite gravity survey technology is by for superior to traditional gravity survey and will surely have larger application space.

 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.

From the viewpoint of chemical composition categorization and structure classification of rocks, an effective method was proposed to identify the lithology of volcanic rocks by using logging data. On the one hand, the conventional logging data could be obtained by core wafer identification. Thus, after processing the data with Support Vector Machines (SVM) method of statistical theory, we could get the lithologic type of the volcanic rocks, which are classified according to the chemical composition of rocks. On the other hand, the volcanic rocks can be classified as volcanic lava, pyroclastic lava and pyroclastic rock according to the rock structure. Typical formation micro-resistivity imaging logging (FMI) image mode can be concluded by establishing the corresponding relationship between FMI images and lithology of volcanic rocks with different structures. As a result, the lithologic type of the volcanic rock classified by rock structure can be determined. Finally, by combining these two kinds of lithology, the ultimate rock lithology can be determined, too. In this paper, the authors presented a novel method to identify the lithology of volcanic rocks by combining SVM processed logging data and FMI image mode.

In recent years, with the implementation of three-dimensional deep geological mapping and shale gas bearing basin prospective evaluation, magnetotelluric method, as one of geophysical methods, has become more and more important. There are a finite number of challenging problems in magnetotelluric data processing and interpretation such as impedance tensor rotation, TE and TM mode identification, static distortion correction, and selection of TE and TM for 2D inversion. Based on magnetotelluric module in WinGlink software, the authors adopted five impedance tensor rotation methods and three curve smoothing methods to process the magnetotelluric data gathered in Jarud Banner, Inner Mongolia. Three stations with different qualities in one profile were chosen to process. The results of five rotation methods and three smoothing methods were analyzed and assessed. Polar diagram analysis of three stations showed that the geological models under themselves were three-dimensional. The comparative study is critical for three-dimensional geological models, and is especially much instructive to the processing and interpretation of magnetotelluric data by using WingGlink software.

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.

This paper used staggered-grid high-order finite-difference scheme to simulate the propagation of source wavefield and reverse time propagation of receiver wavefield to finish elastic prestack reverse time migration in combination with energy flux density vector correlation imaging condition. The results of imaging for a two-layer model and a metal mine model prove that the scheme described in this paper can image the complex structure with steep reflector and eliminate artifacts in reverse time migration for elastic wave.

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.

In order to evaluate environmental effects of the chemical speciation of cadmium, the authors studied soils with high content of cadmium in Guiyang, and investigated the crop, water pollution and endermic diseases. Some conclusions have been reached: Cadmium is existent mainly in exchangeable fraction, residual phase, and Fe-Mn oxide phase, subordinately in weak organic combinational state and carbonate phase, and very lowly in water soluble phase. In vertical profile of soil,cadmium content decreases with increasing soil depth, cadmium is more concentrated in surface soil, and the distribution of chemical speciation of cadmium in deep soil is similar to that of surface soil. The water has not been polluted by cadmium in abnormal areas, and cadmium in such crops as rice, pepper and tobacco has exceed the state standard. In addition, local disease symptoms induced by cadmium has appeared in some areas. Therefore, the crops and human health are adversely affected by the high soil background values of cadmium in Guiyang, and the potential threat is serious.

 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.

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.

The feasibility tests of varied geophysical methods determined the physical character differences between the shell bank and the surrounding media. Data obtained from the Rayleigh wave survey and the reflection wave optimum offset survey were comparatively verified to detect the mode of occurrence of the shell bank in the working area and infer the boundary, buried depth and thickness of the shell bank. These data are very valuable for the protection of the shell bank in the planning of the new urban area of the Gangdong District.

Sedimentary exhalative deposits (SEDEX) widely distributed in the world constitute an important type of ore deposits. They are major sources of Zn and Pb and possess more than 50% and 60% of world Zn and Pb reserves, respectively. In this paper, based on consulting various research achievements and vast literatures concerning this type of ore deposits, the authors provide an overview of several aspects of the sedimentary exhalative deposits, such as their tectonic settings, geological and geochemical characteristics, oreforming model, progress of exploration, and prospecting characteristics, and point out several directions in the further study of  SEDEX deposits.

The faults-enclosed oil/gas accumulation is controlled by faults, and the radioactive anomaly model above the oil/gas accumulation has its unique characteristics. Based on a discussion on the mechanism governing the formation of radioactive anomaly, this paper puts forward a radioactive anomaly pattern above the faults-enclosed oil/gas accumulation. An analysis of the radioactive anomaly characteristics of the radiometric profile above an oil accumulation in Yanji basin shows that the radioactive anomaly pattern above the faults-enclosed oil/gas accumulation does exist.

Using resistivity，radioactivity and density logging curve，the authors analyzed the sedimentary environment of coal seams in some coalfields of Guizhou Province. According to numerical values and shapes of logging curves，the authors hold that coal seams of the exploration area consist of a set of alternate marine and continental sediments，and the coal seams were mainly formed in a peat bog environment of the delta plain. The construction of the sedimentary mode for the coal seam exploration area can guide the exploration of coal mines and help the study of coal accumulation characteristics and the forecast of coal thickness as well as the comparison of coal seams.

 In order to evaluate the performance index of the airborne gravity system, researchers usually use the method for calculating the internal accord accuracy of repeated flight lines. The authors have designed a method for precise calculation of the corresponding point of repeat lines, in which the gravity value is calculated by using interpolation. High precision internal accord accuracy could be achieved by using RMS calculation for gravity value of corresponding points on repeat lines. 

The preliminary study of the geochemical features in Buqingshan area of East Kunlun in Qinghai Province was based on the 1∶50000 stream sediment survey. The study included element distribution characteristics, element anomalies, element correlation and element anomaly assemblage. The result shows that the study area might be a prospective area for Au, Cu, Co, Ni and Sb. It is inferred that ore deposits in Buqingshan area are probably porphyry Cu-Au polymetallic deposits related to granite and granite porphyry, whereas Co-Ni polymetallic deposits are related to the basic and ultrabasic rocks. The Mnite area is expected to be a prospective area for porphyry Cu-Au polymetallic deposits, the Delisitan area is likely to have Co-Ni polymetallic deposits associated with basic and ultrabasic rocks, and the Gerizhuoduo area is expected to be a prospective area for porphyry Cu-Au polymetallic deposits.

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.

During the exploration of the underground medium structure, the gradiometer as well as its new deployment styles for various applications has become an important research topic of world concern. This paper describes various kinds of gradiometers based on torsion balance and difference algorithm and deals with the principle of each instrument. In addition, the absolute gradiometers based on the laser and atom interference are also described. In the last part of this paper, the design of the measurement light path is presented. The results obtained by the authors can serve as the technologic reference for developing gradiometer.

Exemplified by the fathoming data of a survey line using direct current, the sectional drawings of apparent resistivity isograms were plotted by manual and Surfer software respectively. It is concluded that the plotting of the isograms by correct utilization of Surfer software is of great importance. The key points for an ideal effect include the presentation of the x and y value before mesh generation, the choice of the interpolation methods and the correct interpolation. If the plotting is incorrect, the explanation of the data can be directly affected and even a wrong conclusion will be reached.

In large scale gravity prospecting, when total accuracy of Bouguer gravity anomaly is increased to 0.025×10^-5m/s², the height mean square error of the measuring point must be 0.050 m. Using EGM2008 model and GPS leveling data to fit local quasi-geodic surface, the authors transformed the altitude of the measuring point GPS, thus reaching the accuracy. By using GPS leveling data of a certain area, the authors studied the feasibility of the transforming method. The result shows that the method is completely applicable to large scale gravity prospecting.

This paper proposes PML absorbing boundary based on Cerjan attenuation function based on the study of Cerjan absorbing boundary and PML absorbing boundary.Wave field modeling calculations by using two absorbing boundary conditions show that PML absorbing boundary based on Cerjan attenuation function can perform absorption much more clearly and absorb the boundary reflection from various angles.Besides,its operation is easy,the calculation time is short, and it is suitable to seismic wave numerical modeling in two-dimensional situation.

Tilt-angle derivatives method requires no source property and structural index and can quickly estimate magnetic source boundary location and depth distribution, thus having good prospect in rapid interpretation of large-area magnetic anomalies. In this paper, the authors present some algorithm improvements based on Salem's tilt-angle derivatives method. Finite difference method is used in space domain to convert horizontal derivatives of magnetic anomalies, and ISVD algorithm is employed to convert high-level vertical derivatives of magnetic anomalies. A method is also presented which combines the filtering generalized tilt-angle total horizontal derivatives of magnetic anomalies with screening the effectiveness of source depth to remove bad solutions. Theoretical model test shows that this method is simple and quick, and hence its effect is obvious.