The Southwestern Tianshan has a very good potential of gold and non-ferrous mineral resources. 1:50 000 regional geological and mineral survey provided abundant information for the geological and ore-forming background in Bileti area. The regional geochemical characteristics were identified by carrying out stream sediment survey. The authors' work has greatly enriched large scale regional geological and geochemical data. Based on the above work, 3 ore-forming prospective areas have been delineated. Calculation using geochemical quantitative prediction methods (areal productivity method) shows that the potential resource of Au is about 309 kg, and the potential resource of Cu is about 38 168 t in Kalatieke ore-forming prospective area, suggesting a very good mineral prospecting potential in the study area.

The authors conducted 1:50 000 rock structure geochemical survey in the Yiliu orefield of northern Guangdong. Structural fissure filling material, mineralized altered rock, ore veins and limonitization as the sampling media, ore-inducing geochemical anomalies could be effectively delineated, and the delineated anomalies were basically consistent with the stream sediment anomalies in spite of the existence of some differences in anomaly size and intensity. The results show that rock structure anomalies are obviously more consistent with the known ore deposits (ore spots) than stream sediment anomalies. Due to the limitation of the sampling conditions, this method is only suitable for areas with well exposed basement rock and areas dominated by limestone. This method has reference value for geochemical exploration in strongly weathered, thickly covered areas dominated by limestone in southern China.

Deep prospecting has always been difficult and ineffective. Logging as an important means of prospecting for metal ores has been used for many years, but its deep prospecting application is still in the exploration stage, and its effect also needs evaluation. In this paper, the authors took the mining borehole log data from the Fengyan-Xiekeng ore district of Meixian Town in Youxi County of Fujian Province as the research object by analysis of the tri-component magnetic survey, resistivity logging,and IP logging method application effect of deep prospecting in this area. Through the comparison of core observation and analysis, the identification effect of mineralization section and various logging response characteristics were analyzed comprehensively. It is thought that among several logging methods, the resistivity logging is relatively effective in search for lead-zinc deposits, whereas magnetic survey is effective for lead-zinc ore block with high ore grade. On the whole, the two kinds of logging methods are respectively effective to a certain extent, but their effects in search for deep molybdenum ares are not ideal, with no corresponding relationship. Based on a comprehensive analysis of the logging application for deep ore prospecting in this area, the authors hold that, for deep prospecting, researchers should analyze well logging application conditions. In the process of mineral exploration, a more effective means of logging detection should be used to find the corresponding relationship with mineralization, which can help reduce the blindness of regional prospecting work. In the borehole suitable for the exploration conditions, well logging is somewhat effective in the detection and prediction of deep orebody, determination of the attitude and the size of the orebody and locating of its spatial location. Generally, the logging method can play a supplementary role and can improve prospecting success rate.

The Sigou rock mass is a Yanshanian acid rock mass in the middle of the Taihang Mountain and has the same formation period with the Mapeng rock mass. Lots of geological and geophysical work has been carried out in the past years, and porphyry crytptoexplosive breccia type silver and polymetallic sulfide deposits have been found. However, due to the uncertainty of the spatial distribution of the rock mass, the deep exploration prospect has been doubtful. In this paper, the author employed Audio Magnetotelluric (AMT) to divide the boundary of the rock mass based on obvious electrical differences from the surrounding rock, and drew the spatial distribution characteristics. It is believed that the electrical structure of the Sigou rock mass is mainly characterized by "Y" type low resistivity, with the upper part being exposed veins, and the lower part being hidden rock mass. The bifurcation compounding part of the "concave" low-resistance zone is favorable for the formation of cryptoexplosive breccia polymetallic sulfide deposits, whereas the high and low resistance gradient zone (inner zone) is the contact part between the rock mass and the surrounding rock, which is conducive to the formation of alteration belt type polymetallicl sulfide deposits.

In this paper, the authors utilized the high-precision ground magnetic survey and Controlled Source Audio Frequency Magnetotelluric (CSAMT) technique to study the Yuejinshan iron deposit in Qinghai Province. The integrated geological and geophysical survey in this area shows that it is vital to design appropriate strategies and adopt suitable methods for ore explorations in areas of complex topography and geological setting. The choice of strategies and methods has a direct and significant influence on the final survey results. The studies also indicate that a combination of high-precision ground magnetic survey and CSAMT technique can produce highly reliable survey results for this area. The authors made an integrated evaluation for the mine exploration prospect in Yuejinshan area based on an analysis of local geological settings and geophysical anomalies. The geological-geophysical prospecting model developed in this study can serve as a valuable prototype for mine exploration in areas with similar geological evolution settings.

In view of the serious water shortage at the southern edge of Tenggli Desert, the authors carried out conventional resistivity sounding. The resistivity sounding curves are interpreted qualitatively, quantitatively and in detailed by comprehensive geological, hydrogeological and electrical characteristics. It is concluded that the Neogene strata with resistivity of about 25 Ω·m are sandstone, sand gravel is relatively thick and mudstone is relatively thin, suggesting that this horizon is the target horizon in search for cranny pore water of clastic rock in this area. The drilling results show that the conventional resistivity sounding method, which is portable and economical, can be used to find groundwater in the arid desert area of Alxa with severe water shortage.

In order to development the Aero Geophysical Survey and technology application in the key metallogenic zone ,China Aero Geophysical Survey and Remote Sensing Center for Land and Resourc conducted 1:50 000 aeromagnetic surveys in Dunhuang-Akesai from 2015~2016. They organized ground verification and collected magnetic parameters . In the process, they collected a series of Nb-Ta deposits information and date such as YuShishan ore. According to the classification of different types of deposits, the researcher compared this geological datas with aeromagnetic ΔT characteristics. The result shows these Nb-Ta deposits has typical magnetic field zonation indication on the aeromagnetic ΔT anomaly map. And this research has important guiding significance in expanding the known scale of Nb-Ta ore and in finding new rare metal deposits.

Through the study of ore deposit geology of uranium deposits in Jiling pluton of Longshoushan metallogenic belt, the authors hold that, although the ore-bearing sandstones are different, different types of uranium mineralization have some similarity in that the ore deposits all underwent such evolution processes as magmatic differentiation, enrichment at the edge of rock mass, and hydrothermal activity and all exhibit the uranium element differentiation, migration and enrichment mechanism. Deposits of the geophysical field can be compared with each other. Through in depth study and a comparison of physical characteristics of three typical types of uranium deposits (siliceous veins, cryptoexplosive breccia type, alkaline complex type), the authors found that uranium deposits are generally located in the transition region between strong and weak magnetic anomalies with some favor to the weak magnetic anomaly side; in the resistivity profile, uranium deposits are mostly located in high and low resistance variation gradient belt or in the middle and lower resistance zone; in the resistivity plan, uranium deposits are located in the middle and lower resistance transition zone with some favor to the low resistivity side. The area of the electric and magnetic intensity anomaly transition zone with some favor to the weak anomaly side can be used as a new geophysical prospecting mark.

The resource extent of gold is diminishing with thirty years of continuous mining in Dashui gold concentration area of south Gansu Province. In order to obtain a new breakthrough in prospecting in the periphery and deep of the Dashui gold concentration area, 1:50 000 high-precision magnetic measurement of the ground was implemented. With the help of 1:50 000 high accuracy ground magnetic measurement and 2.5D human-computer interaction inversion, the author delineated the geomagnetic field anomalies in the study area and the precision section is explained. The author analyzed the characteristics of magnetic anomalies and the relationship between magnetic anomalies and gold deposits and explained the reasons for the abnormal formation. In addition, the author also summarized the geophysical prospecting signs in the study area. Finally, with the help of the stream sediment survey and geological mapping, the author predict the potential ore-prospecting area and newly discovered gold resources of 10 tons. It shows that magnetic survey is effective in prospecting in the Dashui gold concentration area and the distribution range of acid intrusive reflected by magnetic anomalies is favorable area for exploration of gold deposits.

Wutaishan-Hengshan region is one of the impotant copper and gold polymetallogenic districts in Shanxi Province, theYanshanian rock mass is closely related to ore deposits.In the past, the main object of the study was strong aeromagnetic anomalies,there is not enough attention to weak anomalies.In this area, by use of residual anomaly method of moving average, weak anomalies are extracted from 1:50 000 aeromagnetic data, and combines with geological, gravityand geochemical exploration, the paper infers the presence of hidden rock mass by outdoor inspection of abnormities, and the prospecting targets and delineated for future geological prospecting work.

As a main and important method in regional geological survey, geochemical survey can not only offer theoretical guidance or inference based geological evidence but also provide direct evidence for mineral resources. The technology was implemented in two periods of specification guidance, i.e., DZ/T 0011-91, DZ/T 0011-2015. In this study, the authors find that there exist some differences in the key links of sample collection, analysis, test and data interpretation in different periods and different working processes by different workers, and hence the results might fail to fully reflect the geochemical information. The authors believe that the real and effective sampling (sample representation), accurate data testing (reliability of results), and scientific interpretation of prediction constitute the key to geochemical survey work.

According to aeromagnetic and gravitational data,the boundary and scope of Sichuan basin were determined based on the analyses of the characteristics of gravity and magnetic fields.And on the basis of compiling depth map of the metamorphic basement and structuring zoning map,the research has been 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 to 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 near-northeast and northwest-west construction,which has a pattern of four depression and three uplifts.All these new cognitions and conclusions provide a reference for further oil and gas exploration in the basin.

A research area in Nzega region of Tanzania is covering so severe by fourth system that conventional geological prospecting methods are not available. In order to understand the stratigraphic or lithological boundary of the study area and to construct the production environment, it provides guidance function and effective prospecting information for the selected target area, carries out 1:10 000 the ground high-precision magnetic survey, delineate multiple linear magnetic anomaly zone of beaded displays for NW-SE and NE-SW. Based on the magnetic characteristics of gabbro, It is speculated that the magnetic anomaly belt is caused by intrusive gabbro dyke, thus clarifying the distribution characteristics of the linear structure in the study area. Combination of regional structure and distribution of gold deposits (points), that research the shear zone for NW-SE is the most close structure with gold mineralization relationship, to look for tectonic altered rock type and shear zone type gold deposit;After the surface engineering of the shear tectonic belt (SZ6) for NW-SE was revealed, gold mineralized alteration zone was discovered.

This paper introduces the refineddata processing of magnetization vector inversion and Maxwell plate inversion, and takes the skarn type lead-zinc mine in Galaiaoyi area of Heilongjiang province as an example, summed up a set of data processing, analysis and interpretation of the idea that the magnetization vector inversion with the apparent resistivity profile can analyze the metallogenic environment, analysis and extraction of ore caused by abnormal, and then use Maxwell has a plate-like inversion of the abnormal lot to simulate the distribution of the orebody (low-resistance) in the ground, thus guiding the drilling and laying, and the test results in the known mining area are better.

With the change of season, the physical parameters in the process of freezing and thawing in active layer of seasonal frozen soil region change significantly. Taking the seasonal frozen soil in Northeast China as an example, the authors used the Gaussian distribution rough surfaces to simulate the rough freezing and melting layers, established the random media model which can accurately described heterogeneity of the active layers, and carried out forward modeling. The results show that the depth of freezing and melting layers changes with the seasons, together with the change of permittivity and conductivity. The scattered waves in radar profile are very developed because of the heterogeneous active layer and the undulating freezing and melting layers. With the change of time, the greater the fluctuation of the melting layer, the stronger the scattered wave energy in the radar profile, the harder the reflection of the melting and the freezing layers. At the same time, it is proved that the application of GPR to monitor the seasonal variation, frozen depth and melting depth of seasonal frozen soil is a practical method.

As a device form of transient electromagnetic method, the ground-borehole TEM is being studied by more and more experts both in China and abroad. It has such advantages as small electromagnetic interference,strong signal for its receiving probe in the drilling or downhole tunnel, and closeness to the target anomaly. In this paper, the author used loop source as the excitation source, and established the geological model using finite-difference time-domain method (FDTD) with plate-like body and rectangular low-resistance body. The author calculated the transient response of the rectangular loop source ground-borehole TEM in the half space, and studied its characteristics and laws from the perspective of ground-borehole observation and ground-tunnel observation. The results show that, when multiple anomalies exist simultaneously, the ground-borehole TEM can distinguish different anomalies from the anomalous body and can accurately locate the depth and the horizontal position of different anomalies. The results can provide reference for the quantitative explanation of the depth and location of the anomaly with ground-borehole TEM.

Tomography is one of the most important velocity building methods.Travel time tomography in image domain, implemented with migration, is widely used in velocity model building currently. The authors first introduced the method to output the azimuth-reflection angle gathers, and then, under the assumption of the first-order Born and Rytov approximation of wave equation,started with the imaging condition of Gaussian Beam Migration to derive the linear relation between traveltime perturbation and velocity perturbation in the image domain,with which the authors constructed the explicit expression of kernel function for the wave equation traveltime tomography and established the traveltime tomography equation.The key to computing the kernel is how to compute the Green function in the background model.Making use of the Gaussian beam propagation operator to compute the kernel function can be flexible and efficient.Together with the implementation of Gaussian beam propagation operator in migration,the authors truly realized the integrated technological process of velocity building and migration.Numerical tests and field data application demonstrate that the Gaussian-beam-propagator based traveltime tomography in image domain is effective.

Studying the basic characteristics of near-surface scattering waves and determining the location of near-surface heterogeneities help seismic imaging,survey planning,and avoiding geological disasters.In order to study the wave fields characteristics of near-surface heterogeneous bodies,the authors used the high-order finite difference numerical simulation technique and perturbation theory method to simulate the wave field records of shallow surface scatters and analyze the basic characteristics of near-surface scattering waves.The near surface heterogeneity was located by introducing of reverse time migration imaging technology with the scattering wave field as an extrapolated wave field.The numerical results indicate that surface-wave scatterings are usually stronger than those body-wave scatterings and that forward scatterings are also stronger than backward scatterings.Moreover,near-surface scatters can be precisely imaged by the elastic wave reverse time migration.The near-surface scatters can be regarded as a secondary source,which increases the near-surface illumination.Using the scattered wave field can improve the accuracy of near-surface velocity inversion and seismic imaging.

In order to solve the problem of low prediction resolution for fan delta front thin sand body,the authors,considering basic geological characteristics of Nv59 fault-block in the Shenusi oilfield,utilized the advantages and disadvantages of the constrained sparse spike inversion and geostatistics inversion method,predicted sand body of Ek21 of this block,and evaluated the effectiveness of the geological effect of the calculated result.According to different well densities and different underground geological conditions,the authors chose adaptive inversion result.Compared with conventional methods,geostatistics inversion based on characteristic curve structuring is more suitable for single body prediction in development areas with dense wells.In areas where the wells are scarce,constrained sparse spike inversion is more suitable for prediction of thick sand group.This technique has guided the development of Ek21 oil sands at Nv59 fault block effectively.

The matching pursuit algorithm has higher time-frequency resolution and local adaptive capability and obtains precise position both in time and frequency domain.It is difficult to apply it in practical work due to the low-frequency shadow and low computation efficiency.This paper proposes a method to improve the calculation precision and operation efficiency of Ricker wavelet matching pursuit algorithm.The accuracy of frequency localization and time resolution of time-frequency spectrum is improved by optimizing iteration scheme which improves the calculating method of atom main frequency and compressed atom's sidelobe.The efficiency of algorithm is greatly improved by adaptive acquisition of atomic parameters.By using model experiment and actual seismic data,the effect of improved algorithm is proved.And in some thin reservoir area,the distribution of sandstone is characterized well by the time-frequency slice calculated by improved algorithm.It would help the sedimentary microfacies characterization.

In view of the problem of viscous imaging and the advantage of least-squares reverse time migration (LSRTM),a 3D viscoacoustic LSRTM method is proposed by using the viscoacoustic wave equation of the General Standard Linear Solid (GSLS),which is based on the adjoint operator of the 3D viscoacoustic wave equation and the frame of LSRTM.The polar encoding technology is used to dramatically decrease the computational cost and memory storage of the 3D viscoacoustic LSRTM,which can make it possible to be practical.In the numerical examples,the authors firstly verify the advantage of the proposed method over the traditional acoustic LSRTM using the Marmousi model.Finally,the test of the 3D flat model and salt dome model with Q proves the accuracy of the proposed 3D viscoacoustic LSRTM.

Since the difference of physical properties between formation oil and formation water is small,the physical properties change of sand and mudstone affects seismic reflection,and there exists noise of seismic data,it is very difficult to identify oil-bearing sand bodies using seismic data in that it has strong multiple solutions,and hence the successful application cases are very few.Although it is difficult for single hydrocarbon-bearing sands to generate observable difference,hydrocarbon is accumulated in underground space,and seismic data response from multi-hydrocarbon zone could be identified easily.In this paper,the authors propose peak energy sum hydrocarbon detection method.The method first utilizes the variable factor generalized S transform to achieve frequency division,and then analyzes seismic frequency range which highlights the difference between oil and water and determines energy threshold;on such a basis,the peak energy of the hydrocarbon sensitive frequency section is calculated.For high abundance oil-bearing areas,many sets of hydrocarbon reservoirs are developed longitudinally,and hydrocarbon information can be further highlighted by accumulation of peak energy in the longitudinal direction.Therefore,the peak energy sum hydrocarbon detection factor can effectively predict the high abundance of hydrocarbon distribution.The model data test proves that this method is effective.This method was applied in several oilfields in Bohai,and the results of hydrocarbon detection are in good agreement with the hydrocarbon distribution areas obtained by drilling.This method has guided many successful drilling wells and obtained good effects.

Har Hu area is a part of the frozen soil region of the Tibetan Plateau.The results of geological survey show that there exist favorable accumulation conditions and prospecting potential for gas hydrate in this area.Therefore,high-precision shallow seismic reflection experiments were carried out in this area,with the selection of KZ-28 large vehicle-borne Vibroseis.By analyzing the seismic geologic conditions of the test area in combination with the depth of the shallow seismic exploration and the characteristics of the reflection waves of the target layer,the authors studied the number of the vibration sources,the scanning frequency,the scanning length,vibration number and the driving amplitude before the seismic profile acquisition,and the seismic records obtained by different single excitation parameters were compared and analyzed.Under the prerequisite of satisfying the detection requirements,the seismic exploration can be determined for the purpose of maximizing the seismic record resolution and SNR.The optimal excitation parameters of the Vibroseis are 1 vibrator,2~3 sweeps,the sweep frequency is 8~100 Hz,the sweep length is 12 s,and the drive level is 75%.

The dolomite reservoir,with huge prospecting potential,is widespread in Sichuan basin.Due to various reasons,the prediction of dolomiet reservoir is extremely difficult.In order to guide reservoir prediction,this paper introduces forward modeling to study the seismic response of dolomite reservoir.The seismic response of Permian Qixia dolomite reservoir is not clear.In order to clarify the seismic response characteristics,the authors constructed different forward models in combination with the situation of drilling.The analysis of forward modeling results shows that,on the one hand,reservoir seismic response is weak wave reflection when the underlying strata are high-speed formation and,on the other hand,the response characteristics of the reservoir are influenced by the underlying strata,featured by trough reflection,but still show the weak wave reflection after eliminating the influence of the underlying strata of low speed formation.Meanwhile,the authors combined seismic attributes and high resolution seismic inversion to predict dolomite reservoir.The forward modeling provides a basis for studying seismic response characteristics of dolomite reservoirs and also serves as a means for further study of the dolomite reservoir.

When people evaluate and explain geochemical anomalies, they often focus on "high, large and overall" geochemical anomalies and pay little attention to anomalies with small areas, weak intensity and low superposition. Based on geochemical characteristics of a large number of large and superlarge ore deposits in combination with geochemical practice, the authors hold that there are several kinds of weak anomalies which deserve attention, especially the geochemical anomalies in the porphyry deposits related to tube like rocks, which are usually small in size and low in intensity. However, the size of the deposits related to these small and weak anomalies may be large or even superlarge, resulting in an abnormal phenomenon of "small anomaly corresponds to large ore deposit". This kind of anomalies is easily ignored in the interpretation of evaluation. In addition, regular single point anomalies and small and weak anomalies in the coverage area can represent certain geological significance. These anomalies deserve attention in the interpretation and evaluation of geochemical anomalies.

This paper presents a method for the determination of ultra-trace platinum and ultra-trace palladium in rock ore by graphite furnace atomic absorption spectrometry after extraction from phenylthiourea-octanol. The conditions for extraction of platinum and palladium, the interference of co-existing ions and the optimum instrumentation conditions were studied. The liquid phase extraction was conducted after the sample was dissolved by aqua regia, and then the sample test was measured in the organic phase. The detection limit of platinum and that of palladium were 0.5×10 ^-9 and 0.1×10 ^-9. The relative standard deviations of the method were 6.0% and 4.7%, respectively, for eight samples with platinum 14.7×10 ^-9 and palladium15.2×10 ^-9 from the repeated sampling and dissolution. The relative standard deviations are 6.0% and 4.7% respectively, showing relatively good analytical quality.

In view of the fact that a certain copper mine in Chile is with great terrain elevation difference, sharp vertical cutting and low efficiency in running way of conventional symmetrical schlumberger units, this paper brings forward the proposal of utilization of complex unit combining symmetrical four-electrode and central gradient, that is to say, small polar distance is provided with symmetrical four-electrode detection device whereas large polar distance is provided with central gradient detection device. In the same sounding profile of the mine, the comparison of the sounding inversion results by symmetrical schlumberger unit and complex unit (symmetrical schlumberger + central gradient) indicated that the inversion results of the complex unit equals to the sounding effects of symmetrical schlumberger unit, but the complex unit is more efficient and more suitable for construction in complicated topography areas.

The surface and underground structure of the plateau are complex and changeable, and the static correction problem is distinct. Through years of exploration, the authors gradually formed a series of static correction technologies for plateau salty surface construction: the initial model of complex surface; extremely low SNR standard layer static correction technology; high steep structure comprehensive static correction technology; static correction, surface gun domain iterative detection of saline desert flow; the original surface refraction static correction technology. This series of static correction technologies can effectively solve the problem of static correction in the Western Plateau under the condition of saline land surface. It has played an active role in improving the signal-to-noise ratio and the image quality of the surface seismic data.

Map drawing is one of the important means for researching geological science with multidisciplinary and multipurpose.It has drawn up a solution that the graphics exchange is carried out between the different kinds of mapping software used by SVG.Structures and attributes of graphic entities are described in comparison with MapGIS.Conversion and comparison between various color spaces and processing between SVG and MapGis are described in detail.The related codes of Python are partly presented.Practice shows that the conversion method has the advantages of higher precision,strong commonality and higher efficiency.

Due to historical reasons, the present coordinate systems of the geological, geophysical and geochemical exploration results are basically WGS84, BJ54 or XA80 systems; nevertheless, according to the NASG requirements on the overall use of CGCS2000 coordinates, in the future the coordinates for all kinds of results should be CGCS2000 coordinates, which causes inconvenient situation in using the results and in comprehensive research work. As a professional geographic information system software, ArcGIS has a wide range of applications in various walks. In this paper, the authors briefly introduced the ArcGIS built-in coordinate system, studied and deduced the calculation formula of the MOLODENSKY coordinate conversion method in ArcGIS, put forward the method to get the conversion and precision evaluation of transformation parameters between different ellipsoids, and cited practical examples for verification. On such a basis, the specific ideas and points for attention were analyzed for each coordinate system in the conversion of ArcGIS software to CGCS2000.

This study focuses on the C and O isotope signature for the rocks from the Permian?Triassic boundary in Guanba area of Qianxi, Guizhou Province. The result reveals that δ13C has a distinct negative drift,and the δ18O drops abruptly. Three layers of tuff,recognized at the boundary,have recorded the volcanic events that occurred in late Permian to early Triassic in this area . The volcanic events led to the abrupt change of sedimentary environments and ultimately triggered the mass extinction during the end of Permian.

To understand the nutritional quality, safe quality and growing environment of Chimonobambusa utilis in Nanchuan, Chongqing, the authors analyzed the basic nutritional components, mineral elements and heavy metal elements of Chimonobambusa utilis bamboo shoots, and evaluated the grade of nutrient elements, environmental quality, and Se and Ge in the surrounding soil. The results show that the concentrations of moisture, protein, total sugar, K, Zn and Cu of Chimonobambusa utilis bamboo shoots were high. The bamboo shoots contain 17 kinds of amino acids, and essential amino acids account for 39.9%. Chimonobambusa utilis bamboo shoots are Se-rich food, with the content of Se 0.006×10 ^-6~0.042×10 ^-6 and Se-rich rate 73%. The bamboo shoots are not Ge-rich food, with the content of Ge being 0.0002×10 ^-6~0.0018×10 ^-6. The soil nutrients in the surrounding soil are abundant. Except for K and CaO, the nutrient grades are all relatively rich or rich and the content of Se in soil is high. Some Cd elements in the soil exceed the standard which do not affect the safety quality of the bamboo shoots, and other heavy metals are better in environmental quality.