Near-source noise suppression of AMT by compressive sensing and mathematical morphology filtering

In deep mineral exploration, the acquisition of audio magnetotelluric （AMT） data is severely affected by ambient noise near the observation sites; This near-field noise restricts investigation depths. Mathematical morphological filtering （MMF） proved effective in suppressing large-scale strong and variably shaped noise, typically low-frequency noise, but can not deal with pulse noise of AMT data. We combine compressive sensing and MMF. First we use MMF to suppress the large-scale strong ambient noise; second, we use the improved orthogonal match pursuit （IOMP） algorithm to remove the residual pulse noise. To remove the noise and protect the useful AMT signal, a redundant dictionary that matches with spikes and is insensitive to the useful signal is designed. Synthetic and field data from the Luzong field suggest that the proposed method suppresses the near-source noise and preserves the signal well; thus, better results are obtained that improve the output of either MMF or IOMP.

Controlled-source electromagnetic data processing based on gray system theory and robust estimation

We propose a novel method that combines gray system theory and robust M-estimation method to suppress the interference in controlled-source electromagnetic data. We estimate the standard deviation of the data using a gray model because of the weak dependence of the gray system on data distribution and size. We combine the proposed and threshold method to identify and eliminate outliers. Robust M-estimation is applied to suppress the effect of the outliers and improve the accuracy. We treat the M-estimators of the preserved data as the true data. We use our method to reject the outliers in simulated signals containing noise to verify the feasibility of our proposed method. The processed values are observed to be approximate to the expected values with high accuracy. The maximum relative error is 3.6676%, whereas the minimum is 0.0251%. In processing field data, we observe that the proposed method eliminates outliers, minimizes the root-mean-square error, and improves the reliability of controlled-source electromagnetic data in follow-up processing and interpretation.

Energy distribution and effective components analysis of 2^(n) sequence pseudo-random signal

In order to extract usable harmonics from real 2^(n) sequence pseudo-random data,a technical method is proposed.An equation for predicting the average amplitude of the main frequencies is proposed to guide the choice of signal type for different exploration tasks.By the threshold of the amplitude of the transmitted signal,a set of candidate frequencies are first selected.Then,by operating a spectrum envelope method at these candidate frequencies on received data,effective components in data are extracted.A frequency density calculation method is proposed based on a logical number summation method,to reasonably characterize the frequency density in different frequency bands.By applying this method to real data in Sichuan,China,with signal Type 13,75 effective components are extracted,including both main frequencies and harmonics.The result suggests that the number of effective frequencies in the 2^(n) sequence pseudo-random signal can be increased by extracting usable harmonics,without any additional fieldwork.

Application of pseudo-random frequency domain electromagnetic method in mining areas with strong interferences

Due to the strong electromagnetic interferences and human interference,traditional electromagnetic methods cannot obtain high quality resistivity data of mineral deposits in Chinese mines.The wide field electromagnetic method(WFEM),in which the pseudo-random signal is taken as the transmitter source,can extract high quality resistivity data in areas with sever interference by only measuring the electric field component.We use the WFEM to extract the resistivity information of the Dongguashan mine in southeast China.Compared with the audio magnetotelluric(AMT)method,and the controlled source audio-frequency magnetotelluric(CSAMT) method,the WFEM can obtain data with higher quality and simpler operations.The inversion results indicate that the WFEM can accurately identify the location of the main ore-body,which can be used for deep mine exploration in areas with strong interference.

Dynamic imaging of metallic contamination plume based on self-potential data

A dynamic imaging method for monitoring self-potential data was proposed.Based on the Darcy’s law and Archie’sformulas,a dynamic model was built as a state model to simulate the transportation of metallic ions in porous medium,and theNernst equation was used to calculate the redox potential of metallic ions for observation modeling.Then,the state model andobservation model form an extended Kalman filter cycle to perform dynamic imaging.The noise added synthetic data imaging testshows that the extended Kalman filter can effectively fuse the model evolution and observed self-potential data.The further sandboxmonitoring experiment also demonstrates that the self-potential can be used to monitor the activities of metallic ions and exactlyretrieve the dynamic process of metallic contamination.

Time-domain denoising of time–frequency electromagnetic data

Time–frequency electromagnetic data contain frequency and transient electromagnetic information and can be used to determine the apparent resistivity both in the frequency and time domains.The observation data contains three types of noise:the harmonics interference at 50 Hz,high-frequency random noise,and low-frequency noise.We use frequency-domain bandstop filtering to remove the harmonics interference noise,segmentation and extension median filtering,and fitting of fixed extremes in empirical mode decomposition to remove the high-frequency and low-frequency noise,respectively;furthermore,we base the selection of median filtering window size on the variance and skewness coefficient of the data.We first remove the harmonics interference at 50 Hz,then the high-frequency noise,and finally the low-frequency noise.We test the proposed methodology by using theory and experiments,and we find that the three types of noises are removed,the phase and amplitude information of the signal are maintained,and high-quality waveforms are obtained in the time domain.

Evidence of fluid evolution of Baoshan Cu−Pb−Zn polymetallic deposit:Constraints from in-situ sulfur isotope and trace element compositions of pyrite

In-situ LA-ICP-MS and S isotopes of pyrite from the Baoshan Cu polymetallic deposit were conducted to investigate the ore-forming process and the enrichment mechanism of elements.Three generations of pyrite (Py Ⅰ,Py Ⅱ,and Py Ⅲ) in the skarn-type ores and pyrite in the carbonate-hosted sulfide ores from central,western,and northern(C_Py,W_Py,and N_Py) mining districts are selected for comparison.Compared with Py Ⅰ and Py Ⅲ,the contents of most elements in Py Ⅱ are apparently higher.The As and Se contents are high within a wide range and are decoupled in the growth band of the C_Py.The highest As,Se,and Pb contents were found in W_Py and N_Py.These results indicate the drastic changes in the temperature and fluid mixing during the mineralization.The occurrence of fluctuation and change in temperature and f(O_(2)) was triggered by intermittent pulses of magmatic-hydrothermal fluids,mixing with meteoric water,and water-rock interactions.The sulfur isotopes of all species of pyrite indicated the magmatic source.The change in the f(O_(2)) conditions caused slight differences in the sulfur isotope compositions.Consequently,a metallogenic model was proposed to explain the ore-forming processes.

Distributed wide field electromagnetic method based on high-order 2^(n) sequence pseudo random signal

To make three-dimensional electromagnetic exploration achievable,the distributed wide field electromagnetic method(WFEM)based on the high-order 2^(n) sequence pseudo-random signal is proposed and realized.In this method,only one set of high-order pseudo-random waveforms,which contains all target frequencies,is needed.Based on high-order sequence pseudo-random signal construction algorithm,the waveform can be customized according to different exploration tasks.And the receivers are independent with each other and dynamically adjust the acquisition parameters according to different requirements.A field test in the deep iron ore of Qihe−Yucheng showed that the distributed WFEM based on high-order pseudo-random signal realizes the high-efficiency acquisition of massive electromagnetic data in quite a short time.Compared with traditional controlled-source electromagnetic methods,the distributed WFEM is much more efficient.Distributed WFEM can be applied to the large scale and high-resolution exploration for deep resources and minerals.

Occurrence characteristics and enrichment regularity of indium in pyrite: A case study of Dachang tin ore-field

To reveal the occurrence state and enrichment regularity of the dispersed element indium in pyrite, the petrology,mineralogy, geochemistry, and mineral physics were researched detailedly. The results suggest that the structure of pyrite is mainlycomposed of massive structure, disseminated structure, vein structure, reticular structure, comb structure and so on. Generally, thepyrite coexists with sphalerite, marmatite, pyrrhotite, chalcopyrite, galena, and arsenopyrite. And the texture of pyrite primarilyconsists of the metasomatic texture, solid solution texture, idiomorphic?hypidiomorphic granular texture, and disseminated texture.The content of indium in pyrite ranges from 0.491×10?6 to 65.1×10?6 with an average value of 14.38×10?6. Yet, the indium content inthe Gaofeng deposit is higher than that in the Dafulou and Tongkeng deposit, showing a particularly significant supernormalenrichment. Besides, the cadmium content in pyrite is also higher than other dispersed elements, and similarly the abnormalenrichment of cadmium in the Gaofeng deposit is also very significant. An obvious positive correlation exists between In and Cd, orTl, but a negative correlation between In and Re. It is difficult to find out a positive or negative correlation between In and Ga. Theelement zinc is of great importance to the enrichment of indium, which can possibly facilitate to the migration and crystallization of dispersed element indium.

Sources of rare earth elements REE+Y(REY)in Bayili Coal Mine from Wensu County of Xinjiang,China

By analyses on the trace elements of coals,host-rocks and wall-rocks,this study aims to trace the sources and evaluate the utilization prospects of REE+Y(REY)in coals from Bayili Coal Mine(Wensu,Xinjiang,China).The distribution patterns of REY in the coals are divided into two groups,flat-type and heavy REE-enrichment type(H-type).The REY of the former was mainly from the gneisses of the basement of the coal-bearing basin,and that of the later was partly from hydrothermal solution.The host-rocks show two types of REY patterns,middle REE-enrichment type(M-type)and H-type,which are due to the injection of REY from acid terrestrial water and hydrothermal solution,respectively.Almost all the coal samples are plotted into the promising area on the diagram of percentage of critical elements(REY_(def,rel))vs ratio of sum of critical elements to the sum of excessive elements(C_(outl))and half of the coal samples have high contents of Ga closing to the cut-off grade of Ga deposit as by-product,which indicate that the REY and Ga in Bayili Coal Mine are of utilization prospects as by-product.

Effect of nickel-plated graphite on microstructure and properties of matrix for Fe-based diamond tools

Nickel-plated graphite particles and unmodified graphite particles with different contents were added to the Fe-based diamond composites.The basic properties of those specimens were measured,including relative density,hardness,bending strength,abrasion ratio and holding force coefficient.And also,SEM,XRD and EDS were used to carry out microstructure characterization,phase analysis and element distribution of these specimens.The results show that nickel plating effectively improves the surface wettability of graphite particles.And it is determined that an element diffusion zone is formed on the transition interface between the nickel-plated graphite and the matrix materials,effectively enhancing the interfacial bonding strength.Also,the pores and cracks in the matrix generated by adding the graphite particles are reduced after nickel plating.Thus,the loss of basic properties of the specimens is restrained.But it is found the higher the graphite content is,the weaker the positive effect of nickel plating is.In addition,it is revealed that nickel plating plays a conducive part in the formation of graphite lubricants on the working surface,and nickel-plated graphites can slow down the thermal corrosion of the diamond particles inside the high-temperature sintered specimens.

New correction method for controlled-source electromagnetics source effects

A novel method for source effect correction based on integral equation method is proposed.By taking the electrical horizontal field Ex of current source as an example,the correction method is validated using both simulated data and field data.The results show that the correction method is feasible and effective for isotropic media.When the field data are processed,the correction method normalizes the sources with different geological structures,which eliminates the geological difference among sources,and retains the geological difference among receivers.The correction results are in line with the expectation in whole.

Deep metallogenic mechanism in southeastern China based on receiver function data

Crustal deformation and deep metallogenic mechanisms in southeastern(SE)China are still debated.In this study,we applied the receiver function method to measure crustal thickness and Poisson’s ratio for crustal rock using teleseismic data recorded at 207 seismic stations from China Earthquake Administration Network.The results showed that crustal thickness varied from~27 km in the eastern part to~43 km in the western part of the study area,with an average crustal thickness of 31 km.The crust is thick in the west and thin in the east.The observed Poisson’s ratio for crustal rock was relatively high in the southern Cathaysia Block(CB),with an average of 0.295,while in the Qinling—Dabie terrane,it was relatively low,with an average of 0.257.In the middle of the Yangtze craton and central east of the CB,Poisson’s ratio for crustal rock varied from 0.257 to 0.286.By comparing Poisson’s ratio of the intrusive deposits with that of igneous rocks in volcanic complexes,we deduced that the metallic mineral system might be associated with orogenic and hydrothermal deposits.These results indicated that multistage magma and mineralization in the study area might be attributed to the tectonic-magma-thermal event.The high Poisson’s ratio for crustal rock in the southeastern margin of the CB and northeastern Upper Yangtze Craton might be related to Mesozoic lower crustal mafic partial melt,which provides an important environment for various magmatic intrusions and metallogenies.

Mechanical properties and microstructure evolution of solidified copper tailings with hydantoin epoxy resin

Based on the strength and microstructure tests,the effects of the hydantoin epoxy resin content and curing time on the mechanical properties and microstructure of copper tailings specimens were studied.The results showed that the strength of the solidified specimens was increased to 20.84 MPa with 30 wt.%addition of hydantoin epoxy resin.When the specimens with 10%hydantoin epoxy resin were cured for 7 and 14 d,the strengths were 6.33 and 6.67 MPa respectively,which met the requirements as foundation filler and building materials.The microscopic tests showed that the agglomeration was enhanced and the porosities of the solidified specimens were reduced with increase in the hydantoin epoxy resin content,which could greatly enhance the strength of solidified specimens.

Ore-forming mechanism of Xiajinbao gold deposit in Pingquan, Hebei based on fluid inclusions and stable isotopes

The Xiajinbao gold deposit is located in Yong’an-Xiayingfang-Maojiagou polymetallic metallogenic belt,which is animportant metallogenic belt in North China block.In this paper,we present a detailed study on fluid inclusions and stable isotopes ofthe Xiajinbao gold deposit,Hebei Province,China,aiming at discussing the ore source,evolution of ore-forming fluid andore-forming mechanism of the deposit.The macroscopic geological characteristics,S and Pb isotopic analysis results show that thesource of ore-forming materials is mainly from granitic magma,and subordinately from country rocks.H and O isotopic compositionfeatures indicate that the ore-forming fluid is mainly derived from magmatic water.Fluid inclusion characteristics show that theore-forming fluid experienced boiling during the early mineralization stage,which led to the precipitation of gold.Fluid mixingdominated the precipitation of the ore-forming materials during the middle and late stages.The gold precipitation was caused bywater/rock reaction throughout the whole ore-forming process.

Fast forward modeling of muon transmission tomography based on model voxelization ray energy loss projection

To solve the problems associated with low resolution and high computational effort infinite time,this paper proposes a fast forward modeling method for muon energy loss transmission tomography based on a model voxelization energy loss projection algorithm.First,the energy loss equation for muon transmission tomography is derived from the Bethe–Bloch formula,and the imaging region is then dissected into several units using the model voxelization method.Thereafter,the three-dimensional(3-D)imaging model is discretized into parallel and equally spaced two-dimensional(2-D)slices using the model layering method to realize a dimensional reduction of the 3-D volume data and accelerate the forward calculation speed.Subsequently,the muon energy loss transmission tomography equation is discretized using the ray energy loss projection method to establish a set of energy loss equations for the muon penetration voxel model.Finally,the muon energy loss values at the outgoing point are obtained by solving the projection coefficient matrix of the ray length-weighted model,achieving a significant reduction in the number of muons and improving the computational efficiency.A comparison of our results with the simulation results based on the Monte Carlo method verifies the accuracy and effectiveness of the algorithm proposed in this paper.The metallic mineral identification tests show that the proposed algorithm can quickly identify high-density metallic minerals.The muon energy loss response can accurately identify the boundary of the anomalies and their spatial distribution characteristics.

Three-dimensional numerical simulation of geothermal fi eld in space-wavenumber domain

Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.

3D prestack reverse time migration of ground penetrating radar data based on the normalized correlation imaging condition

The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal and produce a high precision subsurface image when the object is buried in a complicated subsurface environment.To better handle the multi-off set GPR data,we propose a three-dimensional(3D)prestack RTM algorithm.The high-order fi nite diff erence time domian(FDTD)method,with the accuracy of eighth-order in space and second-order in time,is applied to simulate the forward and backward extrapolation electromagnetic fi elds.In addition,we use the normalized correlation imaging condition to obtain pre-stack RTM result and the Laplace fi lter to suppress the low frequency noise generated during the correlation process.The numerical test of 3D simulated GPR data demonstrated that 3D RTM image shows excellent coincidence with the true model.Compared with 2D RTM image,the 3D RTM image can more clearly and accurately refl ect the 3D spatial distribution of the target,and the resolution of the imaging results is far better.Furthermore,the application of observed GPR data further validates the eff ectiveness of the proposed 3D GPR RTM algorithm,and its fi nal image can more reliably guide the subsequent interpretation.

Groundwater resources survey of tongchuan city using the audio magnetotelluric method

The development of Tongchuan City,Shaanxi Province,located in the northwestern region of China,is restricted by water resources.The direct current resistivity and induced polarization sounding methods are typically applied in finding urban groundwater.These methods,however,are not effective due to their complicated topography and geological conditions.The application practice shows that the audio magnetotelluric(AMT)method has a large depth of exploration,high work effi ciency,and high lateral resolution.To investigate the distribution of groundwater resources,we deployed three audio-frequency magnetotelluric profiles in the city area.The impedance tensor information of AMT data is obtained using SSMT2000.AMT data dimension analysis reveals that the two-dimensional structural features of the observation area are obvious.The main structure of the observation area is about 45°northeast,as indicated by structural trend analysis.A shallow two-dimensional electrical profile of 1 km in Tongchuan City is obtained by two-dimensional nonlinear conjugate gradient inversion.Finally,combined with regional geological information,the geological structure characteristics reflected by the electrical profile were obtained along with the detailed characteristics of water-rich structures in the area.The infl uence of the structure on the groundwater distribution was analyzed,and the water-rich areas were identifi ed.This work contributes to the prospective development of Tongchuan City.