This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zer...This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.展开更多
The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)h...Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)has covered whole south China.展开更多
The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity i...The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.展开更多
Based on three-dimensional joint finite element, this paper discusses the theory and methodology about inversionof geodetic data. The FEM and inversion formula is given in detail; also a related code is developed. By ...Based on three-dimensional joint finite element, this paper discusses the theory and methodology about inversionof geodetic data. The FEM and inversion formula is given in detail; also a related code is developed. By use of theGreen's function about 3-D FEM, we invert geodetic measurementS of coseismic deformation of the 1989 Ms=7. 1Loma Prieta earthquake to datermine itS source mechanism. The result indicates that the slip on the fault plane isvery heterogeneous. The maximum slip and shear stress are located about 10 kin to northwest of the eathquakesource, the stress drop is about more than 1 MPa.展开更多
Intermediate acid-complex rock masses with low-density characteristics are the most important prospecting sign in the Beiya area, of western Yunnan province, and provide a physical basis for good gravity exploration. ...Intermediate acid-complex rock masses with low-density characteristics are the most important prospecting sign in the Beiya area, of western Yunnan province, and provide a physical basis for good gravity exploration. It is usually difficult to obtaining solutions in connection with actual geological situations due to the ambiguity of the conventional gravity-processing results and lack of deep constraints. Thus, the three-dimensional (3D) inversion technology is considered as the main channel for reducing the number of solutions and improving the vertical resolution at the current stage. The current study starts from a model test and performs nonlinear 3D density-difference inversion called “model likelihood exploration”, which performs 3D inversion imaging and inversion of the known model while considering the topographic effects. The inversion results are highly consistent with those of the known models. Simultaneously, we consider the Beiya gold mine in Yunnan as an example. The nonlinear 3D densitydifference inversion technology, which is restricted by geological information, is explored to obtain the 3D density body structure below 5 km in the mine area, and the 3D structure of the deep and concealed rock masses are obtained using the density constraints of the intermediate-acid-complex rock masses. The results are well consistent with the surface geological masses and drilling-controlled deep geological masses. The model test and examples both show that the 3D density-difference nonlinear inversion technology can reduce inversion ambiguity, improve resolution, optimize the inversion results, and realize “transparency” in deeply concealed rock masses in ore-concentrated areas,which is useful in guiding the deep ore prospecting.展开更多
Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and ...Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.展开更多
Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotrop...Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic （TDAEM） method is still confined to one dimensional （1D） cases, and the corresponding three-dimensional （3D） numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell＇s equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.展开更多
To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is perf...To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM ＂footprint＂ concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.展开更多
Due to the induced polarization(IP)eff ect,the sign reversal often occurs in timedomain airborne electromagnetic(AEM)data.The inversions that do not consider IP eff ect cannot recover the true umderground electrical s...Due to the induced polarization(IP)eff ect,the sign reversal often occurs in timedomain airborne electromagnetic(AEM)data.The inversions that do not consider IP eff ect cannot recover the true umderground electrical structures.In view of the fact that there are many parameters of airborne induced polarization data in time domain,and the sensitivity diff erence between parameters is large,which brings challenges to the stability and accuracy of the inversion.In this paper,we propose an inversion mehtod for time-domain AEM data with IP effect based on the Pearson correlation constraints.This method uses the Pearson correlation coeffi cient in statistics to characterize the correlation between the resistivity and the chargeability and constructs the Pearson correlation constraints for inverting the objective function to reduce the non uniqueness of inversion.To verify the eff ectiveness of this method,we perform both Occam’s inversion and Pearson correlation constrained inversion on the synthetic data.The experiments show that the Pearson correlation constrained inverison is more accurate and stable than the Occam’s inversion.Finally,we carried out the inversion to a survey dataset with and without IP eff ect.The results show that the data misfit and the continuity of the inverted section are greatly improved when the IP eff ect is considered.展开更多
To reduce the dependence of EM inversion on the choice of initial model and to obtain the global minimum, we apply transdimensional Bayesian inversion to time-domain airborne electromagnetic data. The transdimensional...To reduce the dependence of EM inversion on the choice of initial model and to obtain the global minimum, we apply transdimensional Bayesian inversion to time-domain airborne electromagnetic data. The transdimensional Bayesian inversion uses the Monte Carlo method to search the model space and yields models that simultaneously satisfy the acceptance probability and data fitting requirements. Finally, we obtain the probability distribution and uncertainty of the model parameters as well as the maximum probability. Because it is difficult to know the height of the transmitting source during flight, we consider a fixed and a variable flight height. Furthermore, we introduce weights into the prior probability density function of the resistivity and adjust the constraint strength in the inversion model by changing the weighing coefficients. This effectively solves the problem of unsatisfactory inversion results in the middle high-resistivity layer. We validate the proposed method by inverting synthetic data with 3% Gaussian noise and field survey data.展开更多
In this paper, we analyze lithospheric density distribution of China and surrounding regions on the basis of 30＇ × 30＇ gravity data and 1°×1 °P-wave velodty data, Firstly, we used the empirical e...In this paper, we analyze lithospheric density distribution of China and surrounding regions on the basis of 30＇ × 30＇ gravity data and 1°×1 °P-wave velodty data, Firstly, we used the empirical equation be- tween the density and the P-wave velocity difference as the base of the initial model of the Asian lithospheric density. Secondly, we calculated the gravity anomaly, caused by the Moho discontinuity and the sedimentary layer discontinuity, by the Parker formula. Thirdly, the gravity anomaly of the spherical harmonics with 2 40 order for the anomalous body below the lithosphere is calculated based on the model of EGM96. Finally, by using Algebra Reconstruction Techniques （ART）, the inversion of 30＇ ~ 30＇ residual lithospheric Bouguer gravity anomaly caused by the lithosphere yields a rather detailed struc- tural model. The results show that the lithospheric density distribution of China and surrounding regions has a certain connection with the tectonic structure. The density is relatively high in the Philippine Sea plate, Japan Sea, the Indian plate, the Kazakhstan shield and the Western Siberia plain, whereas the Tibetan Plateau has low-density characteristics. The minimum value of density lies in the north of Philippines, in the Taiwan province and in the Ryukyu island arc.展开更多
The Beishan orogen,located in the central segment of the Tianshan–Solonker suture within the southern Central Asian Orogenic Belt(CAOB),is crucial for understanding the accretionary processes and continental growth i...The Beishan orogen,located in the central segment of the Tianshan–Solonker suture within the southern Central Asian Orogenic Belt(CAOB),is crucial for understanding the accretionary processes and continental growth in Central Asia.This orogen developed through the episodic amalgamation and accretion of continental margin arcs,island arcs,ophiolites,and accretionary wedges,undergoing a complex process of accretion and evolution.Since the Phanerozoic,the Beishan orogen has experienced multiple phases of magmatic and collision events.The intricate distribution of magmatic arc rocks has obscured the complete basement traces,and the spatial superposition of multiple magmatic arc phases has complicated the study of its evolutionary history.展开更多
A novel method based on the relevance vector machine（RVM） for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.T...A novel method based on the relevance vector machine（RVM） for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine（SVM） based approach.展开更多
We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular...We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.展开更多
There is promising indication of helium-rich natural gas in Weihe Basin.However,the lack of a detailed investigation on the origin and the spatial distribution of helium source rocks(mainly Yanshanian granites)limits ...There is promising indication of helium-rich natural gas in Weihe Basin.However,the lack of a detailed investigation on the origin and the spatial distribution of helium source rocks(mainly Yanshanian granites)limits the helium potentiality evaluation in Weihe Basin(WB).We performed three-dimensional(3 D)inversions of magnetic data in Weihe Basin and its adjacent areas to figure out the crustal thermal structure and the temporal-spatial distribution of deep granite in the basin.Based on this,we have proposed a geological model of helium accumulation and predicted the potential area of helium distribution.The results of 3 D magnetic inversion indicate that the granites in the deep Weihe Basin are mostly located at the central and southern parts of the basin,which are connected spatially with the granites in the North Qinling Orogenic Belt.These granites were all derived from largescale intra-continental orogeny in the Qinling area during the mid-Mesozoic,providing a good material basis for crust-derived helium in the basin.The local uplift of the Curie surface makes the thermal fluid more actively,which contributes to helium accumulation.The faults developed in the WB are the migration pathway of crust-derived helium and the upward migration of the mantle-derived helium.The wells with high percentage helium are mostly located near the Weihe fault and the areas on the south of it.The Wugong-Xi’an-Lantian area in the central and the Lintong-Weinan-Tongguan area in the eastern Weihe Basin are the most promising helium distribution areas.Furthermore,the region from the north of Taibai Mountain to Baoji City in the western Weihe Basin may also be another potential area of Helium resource.展开更多
Local CO2 laser treatment has proved to be an effective method to prevent the 351-nm laser-induced damage sites in a fused silica surface from exponentially growing, which is responsible for limiting the lifetime of o...Local CO2 laser treatment has proved to be an effective method to prevent the 351-nm laser-induced damage sites in a fused silica surface from exponentially growing, which is responsible for limiting the lifetime of optics in high fluence laser systems. However, the CO2 laser induced ablation crater is often surrounded by a raised rim at the edge, which can also result in the intensification of transmitted ultraviolet light that may damage the downstream optics. In this work, the three-dimensional finite-difference time-domain method is developed to simulate the distribution of electrical field intensity in the vicinity of the CO2 laser mitigated damage site located in the exit subsurface of fused silica. The simulated results show that the repaired damage sites with raised rims cause more notable modulation to the incident laser than those without rims.Speeifically,we present a theoretical model of using dimpled patterning to control the rim structure around the edge of repaired damage sites to avoid damage to downstream optics. The calculated results accord well with previous experimental results and the underlying physical mechanism is analysed in detail展开更多
The 2,026 earthquake events registered by the Sichuan regional digital seismic network and mobile seismic array after the April 20th, 2013 Lushan earthquake and 28,188 pieces of data were selected to determine direct ...The 2,026 earthquake events registered by the Sichuan regional digital seismic network and mobile seismic array after the April 20th, 2013 Lushan earthquake and 28,188 pieces of data were selected to determine direct P waves arrival times. We applied the tomographic method to inverse the characteristics of the velocity structure for the three-dimensional （3D） P wave in the mid-upper crust of the seismic source region of the Lushan earthquake. The imaging results were combined with the apparent magnetization inversion and magnetotelluric （MT） sounding retest data to comprehensively study the causes of the deep seismogenic environment in the southern section of the Longmenshan fault zone and explore the formation of the Lushan earthquake. Research has shown that there are obvious differences in velocity structure and magnetic distribution between the southern and northern sections of the Longmenshan fault zone. The epicenter of the Lushan earthquake is located near the boundary of the high and low-velocity anomalies and favorable for a high-velocity section. Moreover, at the epicenter of the Lushan earthquake located on the magnetic dome boundary of Ya＇an, the development of high velocity and magnetic solid medium favors the accumulation and release of strain energy. Low- velocity anomalies are distributed underneath the are of seismogenic origin, The inversion results of the MT retest data after the April 20th Lushan earthquake also indicate that there a high-conductor anomaly occurs under the area of seismogenic origin of the Lushan earthquake, Therefore, we speculated that the presence of a high-conductivity anomaly and low-velocity anomaly underneath the seismogenic body of the Lushan earthquake could be related to the existence of fluids. The role of fluids caused the weakening of the seismogenic layer inside the mid-upper crust and resulted in a seismogenic fault that was prone to rupture and pIayed a triggering role in the Lushan earthquake.展开更多
One of the main factors of laser induced damage is the modulation to incident laser which is caused by the defect in the subsurface of the fused silica. In this work, the repaired damage site irradiated by CO2 laser i...One of the main factors of laser induced damage is the modulation to incident laser which is caused by the defect in the subsurface of the fused silica. In this work, the repaired damage site irradiated by CO2 laser is simplified to a Gaussian rotation according to the corresponding experimental results. Then, the three-dimensional finite-difference time-domain method is employed to simulate the electric field intensity distribution in the vicinity of this kind of defect in fused silica front subsurface. The simulated results show that the modulation is notable, the Emax is about 2.6 times the irradiated electric field intensity in the fused silica with the damage site （the width is 1.5 μm and depth is 2.3 μm） though the damage site is repaired by CO2 laser. The phenomenon and the theoretical result of the annular laser enhancement existed on the rear surface are first verified effectively, which agrees well with the corresponding experimental results. The relations between the maximal electric field intensity in fused silica with defect depth and width are given respectively. Meanwhile, the corresponding physical mechanism is analysed theoretically in detail.展开更多
Technical development oriented on the detailed seismic imaging and velocity model building coupled with rapid increase of the computing power available nowadays make it possible to process large volumes of seismic dat...Technical development oriented on the detailed seismic imaging and velocity model building coupled with rapid increase of the computing power available nowadays make it possible to process large volumes of seismic data using numerically intensive approaches based the wavefield propagation.展开更多
基金supported by the Development of airborne gravity gradiometer(No.2017YFC0601601)open subject of Key Laboratory of Petroleum Resources Research,Institute of Geology and Geophysics,Chinese Academy of Sciences(No.KLOR2018-8)
文摘This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
基金co-supported by the China Geological Survey Project(DD20190012 and DD20160082)
文摘Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)has covered whole south China.
基金supported by the Key Laboratory of Geological Survey and Evaluation of Ministry of Education (China University of Geosciences)(No. GLAB2020ZR13)
文摘The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.
文摘Based on three-dimensional joint finite element, this paper discusses the theory and methodology about inversionof geodetic data. The FEM and inversion formula is given in detail; also a related code is developed. By use of theGreen's function about 3-D FEM, we invert geodetic measurementS of coseismic deformation of the 1989 Ms=7. 1Loma Prieta earthquake to datermine itS source mechanism. The result indicates that the slip on the fault plane isvery heterogeneous. The maximum slip and shear stress are located about 10 kin to northwest of the eathquakesource, the stress drop is about more than 1 MPa.
基金The authors would like to thank the China Geological Survey (DD20190033)National Natural Science Foundation (41804144) for the financial support,Yunnan Gold and Mineral Group Co.,Ltd. for providing the original geological information,and the reviewers for providing valuable comments.
文摘Intermediate acid-complex rock masses with low-density characteristics are the most important prospecting sign in the Beiya area, of western Yunnan province, and provide a physical basis for good gravity exploration. It is usually difficult to obtaining solutions in connection with actual geological situations due to the ambiguity of the conventional gravity-processing results and lack of deep constraints. Thus, the three-dimensional (3D) inversion technology is considered as the main channel for reducing the number of solutions and improving the vertical resolution at the current stage. The current study starts from a model test and performs nonlinear 3D density-difference inversion called “model likelihood exploration”, which performs 3D inversion imaging and inversion of the known model while considering the topographic effects. The inversion results are highly consistent with those of the known models. Simultaneously, we consider the Beiya gold mine in Yunnan as an example. The nonlinear 3D densitydifference inversion technology, which is restricted by geological information, is explored to obtain the 3D density body structure below 5 km in the mine area, and the 3D structure of the deep and concealed rock masses are obtained using the density constraints of the intermediate-acid-complex rock masses. The results are well consistent with the surface geological masses and drilling-controlled deep geological masses. The model test and examples both show that the 3D density-difference nonlinear inversion technology can reduce inversion ambiguity, improve resolution, optimize the inversion results, and realize “transparency” in deeply concealed rock masses in ore-concentrated areas,which is useful in guiding the deep ore prospecting.
基金supported by the National Natural Science Foundation of China (Grant No. 40974039)High-Tech Research and Development Program of China (Grant No.2006AA06205)Leading Strategic Project of Science and Technology, Chinese Academy of Sciences (XDA08020500)
文摘Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.
基金financially supported by National Nonprofit institute Research Grant of IGGE(Nos.AS2017J06,AS2017Y04,and AS2016J10)Survey on coastal area for airborne magnetic method of UNV in Jiangsu(No.DD20160151-03)+3 种基金Key National Research Project of China(No.2017YFC0601900)Key Program of National Natural Science Foundation of China(No.41530320)Natural Science Foundation(No.41274121)China Natural Science Foundation for Young Scientists(No.41404093)
文摘Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic （TDAEM） method is still confined to one dimensional （1D） cases, and the corresponding three-dimensional （3D） numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell＇s equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.
基金supported by the Key Natural Science Foundation(No.41530320)Natural Science Foundation(No.41274121)+1 种基金Natural Science Foundation for young scientist(No.41404093)the Projects on the Development of the Key Equipment of Chinese Academy of Science(No.ZDYZ2012-1-03)
文摘To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM ＂footprint＂ concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.
基金This paper was fi nancially supported by the National Natural Science Foundation of China(Nos.42030806,41774125,41904104,41804098)the Pioneer Project of Chinese Academy of Sciences(No.XDA14020102).
文摘Due to the induced polarization(IP)eff ect,the sign reversal often occurs in timedomain airborne electromagnetic(AEM)data.The inversions that do not consider IP eff ect cannot recover the true umderground electrical structures.In view of the fact that there are many parameters of airborne induced polarization data in time domain,and the sensitivity diff erence between parameters is large,which brings challenges to the stability and accuracy of the inversion.In this paper,we propose an inversion mehtod for time-domain AEM data with IP effect based on the Pearson correlation constraints.This method uses the Pearson correlation coeffi cient in statistics to characterize the correlation between the resistivity and the chargeability and constructs the Pearson correlation constraints for inverting the objective function to reduce the non uniqueness of inversion.To verify the eff ectiveness of this method,we perform both Occam’s inversion and Pearson correlation constrained inversion on the synthetic data.The experiments show that the Pearson correlation constrained inverison is more accurate and stable than the Occam’s inversion.Finally,we carried out the inversion to a survey dataset with and without IP eff ect.The results show that the data misfit and the continuity of the inverted section are greatly improved when the IP eff ect is considered.
基金This paper was financially supported by the Key National Research Project of China （Nos. 2017YFC0601900 and 2016YFC0303100）, and the Key Program of National Natural Science Foundation of China （No. 41530320） and Surface Project （No. 41774125）.
文摘To reduce the dependence of EM inversion on the choice of initial model and to obtain the global minimum, we apply transdimensional Bayesian inversion to time-domain airborne electromagnetic data. The transdimensional Bayesian inversion uses the Monte Carlo method to search the model space and yields models that simultaneously satisfy the acceptance probability and data fitting requirements. Finally, we obtain the probability distribution and uncertainty of the model parameters as well as the maximum probability. Because it is difficult to know the height of the transmitting source during flight, we consider a fixed and a variable flight height. Furthermore, we introduce weights into the prior probability density function of the resistivity and adjust the constraint strength in the inversion model by changing the weighing coefficients. This effectively solves the problem of unsatisfactory inversion results in the middle high-resistivity layer. We validate the proposed method by inverting synthetic data with 3% Gaussian noise and field survey data.
基金supported by Project SinoProbe-02: Experiment and Integration of Deep Probe Techniques in ChinaNational Natural Science Foundation of China (NSFC, Grant No. 40874067)the Research Fund for the Doctoral Program of Higher Education (Grant No. 20070491520)
文摘In this paper, we analyze lithospheric density distribution of China and surrounding regions on the basis of 30＇ × 30＇ gravity data and 1°×1 °P-wave velodty data, Firstly, we used the empirical equation be- tween the density and the P-wave velocity difference as the base of the initial model of the Asian lithospheric density. Secondly, we calculated the gravity anomaly, caused by the Moho discontinuity and the sedimentary layer discontinuity, by the Parker formula. Thirdly, the gravity anomaly of the spherical harmonics with 2 40 order for the anomalous body below the lithosphere is calculated based on the model of EGM96. Finally, by using Algebra Reconstruction Techniques （ART）, the inversion of 30＇ ~ 30＇ residual lithospheric Bouguer gravity anomaly caused by the lithosphere yields a rather detailed struc- tural model. The results show that the lithospheric density distribution of China and surrounding regions has a certain connection with the tectonic structure. The density is relatively high in the Philippine Sea plate, Japan Sea, the Indian plate, the Kazakhstan shield and the Western Siberia plain, whereas the Tibetan Plateau has low-density characteristics. The minimum value of density lies in the north of Philippines, in the Taiwan province and in the Ryukyu island arc.
基金supported by the China Geological Survey(Grant No.DD20230254)。
文摘The Beishan orogen,located in the central segment of the Tianshan–Solonker suture within the southern Central Asian Orogenic Belt(CAOB),is crucial for understanding the accretionary processes and continental growth in Central Asia.This orogen developed through the episodic amalgamation and accretion of continental margin arcs,island arcs,ophiolites,and accretionary wedges,undergoing a complex process of accretion and evolution.Since the Phanerozoic,the Beishan orogen has experienced multiple phases of magmatic and collision events.The intricate distribution of magmatic arc rocks has obscured the complete basement traces,and the spatial superposition of multiple magmatic arc phases has complicated the study of its evolutionary history.
基金Project supported by the National Natural Science Foundation of China (Grant No. 61071022)the Graduate Student Research and Innovation Program of Jiangsu Province,China (Grant No. CXZZ11-0381)
文摘A novel method based on the relevance vector machine（RVM） for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine（SVM） based approach.
基金supported by the Young Scientists Fund of the National Natural Science Foundation of China(No.41304082)the China Postdoctoral Science Foundation(No.2016M590731)+2 种基金the Young Scientists Fund of the Natural Science Foundation of Hebei Province(No.D2014403011)the Program for Young Excellent Talents of Higher Education Institutions of Hebei Province(No.BJ2016046)the Geological survey project of China Geological Survey(No.1212011121197)
文摘We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.
基金supported by the National Natural Science Foundation of China(No.41904115)the Natural Science Basic Research Plan in Shaanxi Province of China(No.2018JQ4034)。
文摘There is promising indication of helium-rich natural gas in Weihe Basin.However,the lack of a detailed investigation on the origin and the spatial distribution of helium source rocks(mainly Yanshanian granites)limits the helium potentiality evaluation in Weihe Basin(WB).We performed three-dimensional(3 D)inversions of magnetic data in Weihe Basin and its adjacent areas to figure out the crustal thermal structure and the temporal-spatial distribution of deep granite in the basin.Based on this,we have proposed a geological model of helium accumulation and predicted the potential area of helium distribution.The results of 3 D magnetic inversion indicate that the granites in the deep Weihe Basin are mostly located at the central and southern parts of the basin,which are connected spatially with the granites in the North Qinling Orogenic Belt.These granites were all derived from largescale intra-continental orogeny in the Qinling area during the mid-Mesozoic,providing a good material basis for crust-derived helium in the basin.The local uplift of the Curie surface makes the thermal fluid more actively,which contributes to helium accumulation.The faults developed in the WB are the migration pathway of crust-derived helium and the upward migration of the mantle-derived helium.The wells with high percentage helium are mostly located near the Weihe fault and the areas on the south of it.The Wugong-Xi’an-Lantian area in the central and the Lintong-Weinan-Tongguan area in the eastern Weihe Basin are the most promising helium distribution areas.Furthermore,the region from the north of Taibai Mountain to Baoji City in the western Weihe Basin may also be another potential area of Helium resource.
基金supported by the Fundamental Research Funds for the Central Universities (Grant No. ZYGX2010J045)the National Natural Science Fundation of China and the China Academy of Engineering Physics United Foundation (NSAF) (Grant No. 11076008)the Foundation for Young Scholars of University of Electronic Science and Technology of China (Grant No. L08010401JX0806)
文摘Local CO2 laser treatment has proved to be an effective method to prevent the 351-nm laser-induced damage sites in a fused silica surface from exponentially growing, which is responsible for limiting the lifetime of optics in high fluence laser systems. However, the CO2 laser induced ablation crater is often surrounded by a raised rim at the edge, which can also result in the intensification of transmitted ultraviolet light that may damage the downstream optics. In this work, the three-dimensional finite-difference time-domain method is developed to simulate the distribution of electrical field intensity in the vicinity of the CO2 laser mitigated damage site located in the exit subsurface of fused silica. The simulated results show that the repaired damage sites with raised rims cause more notable modulation to the incident laser than those without rims.Speeifically,we present a theoretical model of using dimpled patterning to control the rim structure around the edge of repaired damage sites to avoid damage to downstream optics. The calculated results accord well with previous experimental results and the underlying physical mechanism is analysed in detail
基金supported by China earthquake scientific array exploration-northern section of North South seismic belt (20130811)National Natural Science Foundation of China (41474057)Science for earthquake Resllience of China Earthquake Administration (XH15040Y)
文摘The 2,026 earthquake events registered by the Sichuan regional digital seismic network and mobile seismic array after the April 20th, 2013 Lushan earthquake and 28,188 pieces of data were selected to determine direct P waves arrival times. We applied the tomographic method to inverse the characteristics of the velocity structure for the three-dimensional （3D） P wave in the mid-upper crust of the seismic source region of the Lushan earthquake. The imaging results were combined with the apparent magnetization inversion and magnetotelluric （MT） sounding retest data to comprehensively study the causes of the deep seismogenic environment in the southern section of the Longmenshan fault zone and explore the formation of the Lushan earthquake. Research has shown that there are obvious differences in velocity structure and magnetic distribution between the southern and northern sections of the Longmenshan fault zone. The epicenter of the Lushan earthquake is located near the boundary of the high and low-velocity anomalies and favorable for a high-velocity section. Moreover, at the epicenter of the Lushan earthquake located on the magnetic dome boundary of Ya＇an, the development of high velocity and magnetic solid medium favors the accumulation and release of strain energy. Low- velocity anomalies are distributed underneath the are of seismogenic origin, The inversion results of the MT retest data after the April 20th Lushan earthquake also indicate that there a high-conductor anomaly occurs under the area of seismogenic origin of the Lushan earthquake, Therefore, we speculated that the presence of a high-conductivity anomaly and low-velocity anomaly underneath the seismogenic body of the Lushan earthquake could be related to the existence of fluids. The role of fluids caused the weakening of the seismogenic layer inside the mid-upper crust and resulted in a seismogenic fault that was prone to rupture and pIayed a triggering role in the Lushan earthquake.
基金Project supported by the National High Technology Research and Development Program of China (Grant No.2008AA8040508)the Youth Foundation of University of Science and Technology of China (Grant Nos.L08010401JX0834 and L08010401JX0806)the Fundamental Research Funds for the Central Universities of China
文摘One of the main factors of laser induced damage is the modulation to incident laser which is caused by the defect in the subsurface of the fused silica. In this work, the repaired damage site irradiated by CO2 laser is simplified to a Gaussian rotation according to the corresponding experimental results. Then, the three-dimensional finite-difference time-domain method is employed to simulate the electric field intensity distribution in the vicinity of this kind of defect in fused silica front subsurface. The simulated results show that the modulation is notable, the Emax is about 2.6 times the irradiated electric field intensity in the fused silica with the damage site （the width is 1.5 μm and depth is 2.3 μm） though the damage site is repaired by CO2 laser. The phenomenon and the theoretical result of the annular laser enhancement existed on the rear surface are first verified effectively, which agrees well with the corresponding experimental results. The relations between the maximal electric field intensity in fused silica with defect depth and width are given respectively. Meanwhile, the corresponding physical mechanism is analysed theoretically in detail.
基金by the SEISCOPE consortium(http://seiscope2.osug.fr),sponsored by AKERBP,CGG,CHEVRON,EQUINOR,EXXON-MOBIL,JGI,SHELL,SINOPEC,SISPROBE,and TOTALthe Polish National Science Center,(Grant No.2019/33/B/ST10/01014)
文摘Technical development oriented on the detailed seismic imaging and velocity model building coupled with rapid increase of the computing power available nowadays make it possible to process large volumes of seismic data using numerically intensive approaches based the wavefield propagation.
