The inversions of complex geophysical data always solve multi-parameter, nonlinear, and multimodal optimization problems. Searching for the optimal inversion solutions is similar to the social behavior observed in swa...The inversions of complex geophysical data always solve multi-parameter, nonlinear, and multimodal optimization problems. Searching for the optimal inversion solutions is similar to the social behavior observed in swarms such as birds and ants when searching for food. In this article, first the particle swarm optimization algorithm was described in detail, and ant colony algorithm improved. Then the methods were applied to three different kinds of geophysical inversion problems: (1) a linear problem which is sensitive to noise, (2) a synchronous inversion of linear and nonlinear problems, and (3) a nonlinear problem. The results validate their feasibility and efficiency. Compared with the conventional genetic algorithm and simulated annealing, they have the advantages of higher convergence speed and accuracy. Compared with the quasi-Newton method and Levenberg-Marquardt method, they work better with the ability to overcome the locally optimal solutions.展开更多
3-D geochemical subsurface models,as constructed by spatial interpolation of drill-core assays,are valuable assets across multiple stages of the mineral industry's workflow.However,the accuracy of such models is l...3-D geochemical subsurface models,as constructed by spatial interpolation of drill-core assays,are valuable assets across multiple stages of the mineral industry's workflow.However,the accuracy of such models is limited by the spatial sparsity of the underlying drill-core,which samples only a small fraction of the subsurface.This limitation can be alleviated by integrating collocated 3-D models into the interpolation process,such as the 3-D rock property models produced by modern geophysical inversion procedures,provided that they are sufficiently resolved and correlated with the interpolation target.While standard machine learning algorithms are capable of predicting the target property given these data,incorporating spatial autocorrelation and anisotropy in these models is often not possible.We propose a Gaussian process regression model for 3-D geochemical interpolation,where custom kernels are introduced to integrate collocated 3-D rock property models while addressing the trade-off between the spatial proximity of drill-cores and the similarities in their collocated rock properties,as well as the relative degree to which each supporting 3-D model contributes to interpolation.The proposed model was evaluated for 3-D modelling of Mg content in the Kevitsa Ni-Cu-PGE deposit based on drill-core analyses and four 3-D geophysical inversion models.Incorporating the inversion models improved the regression model's likelihood(relative to a purely spatial Gaussian process regression model)when evaluated at held-out test holes,but only for moderate spatial scales(100 m).展开更多
The paper presents a short overview about the application of joint inversion in geophysics. It gives also an alternative explanation for the term of “different data sets” and discusses what types of inversion proced...The paper presents a short overview about the application of joint inversion in geophysics. It gives also an alternative explanation for the term of “different data sets” and discusses what types of inversion procedures can be considered as joint inversion. Nowadays there are no standard standpoints using the appellation joint inversion. What is joint inversion? Based on the information matrix an answer could be given for this question what could be regarded as various types of data sets that are inverted simultaneously. We would like to expand the explanation—that is professed by many researchers—of the method that regards only the simultaneous inversion of data sets based on different physical parameters as joint inversion.展开更多
In this paper we discuss the edge-preserving regularization method in the reconstruction of physical parameters from geophysical data such as seismic and ground-penetrating radar data. In the regularization method a p...In this paper we discuss the edge-preserving regularization method in the reconstruction of physical parameters from geophysical data such as seismic and ground-penetrating radar data. In the regularization method a potential function of model parameters and its corresponding functions are introduced. This method is stable and able to preserve boundaries, and protect resolution. The effect of regularization depends to a great extent on the suitable choice of regularization parameters. The influence of the edge-preserving parameters on the reconstruction results is investigated and the relationship between the regularization parameters and the error of data is described.展开更多
The Shuangjianzishan deposit in Inner Mongolia is a typical Ag-Pb-Zn deposit of the southern Great Xing’an Range.Proven reserves of Ag,Pb,and Zn in this deposit have reached the scale of super-large deposits,with fav...The Shuangjianzishan deposit in Inner Mongolia is a typical Ag-Pb-Zn deposit of the southern Great Xing’an Range.Proven reserves of Ag,Pb,and Zn in this deposit have reached the scale of super-large deposits,with favorable metallogenic conditions,strong prospecting signs,and high metallogenic potential.This paper reports a study involving integrated geophysical methods,including controlled-source audio-frequency magnetotelluric,gravity,magnetic,and shallow-seismic-reflection methods,to determine the spatial distribution of ore-controlling structures and subsurface intrusive rock for a depth range of<2000 m in the Shuangjianzishan ore district.The objective of this study is to construct a metallogenic model of the ore district and provide a scientific basis for the exploration of similar deposits in the deep and surrounding regions.We used three-dimensional inversion for controlled-source audio-frequency magnetotelluric data based on the limited memory quasi-Newton algorithm,and three-dimensional physical-property inversion for the gravity and magnetic data to obtain information about the subsurface distribution of ore-controlling structures and intrusive rocks.Under seismic reflection results,regional geology,petrophysical properties,and borehole information,the geophysical investigation shows that the Dashizhai group,which contains the main ore-bearing strata in the ore district,is distributed within a depth range of<1239 m,and is thick in the Xinglongshan ore block and the eastern part of the Shuangjianzishan ore block.The mineralization is spatially associated with a fault system characterized by NE-,NW-,and N-trending faults.The magnetic and electrical models identify large,deep bodies of intrusive rock that are inferred to have been involved in mineralization,with local shallow emplacement of some intrusions.Combining the subsurface spatial distributions of ore-bearing strata,ore-controlling faults,and intrusive rock,we propose two different metallogenic models for the Shuangjianzishan ore district,which provide a scientific basis for further prospecting in the deep regions of the ore district and surrounding areas.展开更多
Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuit...Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.展开更多
To minimize the number of solutions in 3D resistivity inversion, an inherent problem in inversion, the amount of data considered have to be large and prior constraints need to be applied. Geological and geophysical da...To minimize the number of solutions in 3D resistivity inversion, an inherent problem in inversion, the amount of data considered have to be large and prior constraints need to be applied. Geological and geophysical data regarding the extent of a geological anomaly are important prior information. We propose the use of shape constraints in 3D electrical resistivity inversion, Three weighted orthogonal vectors (a normal and two tangent vectors) were used to control the resistivity differences at the boundaries of the anomaly. The spatial shape of the anomaly and the constraints on the boundaries of the anomaly are thus established. We incorporated the spatial shape constraints in the objective function of the 3D resistivity inversion and constructed the 3D resistivity inversion equation with spatial shape constraints. Subsequently, we used numerical modeling based on prior spatial shape data to constrain the direction vectors and weights of the 3D resistivity inversion. We established a reasonable range between the direction vectors and weights, and verified the feasibility and effectiveness of using spatial shape prior constraints in reducing excessive structures and the number of solutions. We applied the prior spatially shape-constrained inversion method to locate the aquifer at the Guangzhou subway. The spatial shape constraints were taken from ground penetrating radar data. The inversion results for the location and shape of the aquifer agree well with drilling data, and the number of inversion solutions is significantly reduced.展开更多
Offshore carbon dioxide(CO_(2)) storage is an effective method for reducing greenhouse gas emissions. However, when using traditional seismic wave methods to monitor the migration of sequestration CO_(2) plumes, the c...Offshore carbon dioxide(CO_(2)) storage is an effective method for reducing greenhouse gas emissions. However, when using traditional seismic wave methods to monitor the migration of sequestration CO_(2) plumes, the characteristics of wave velocity changes tend to become insignificant beyond a certain limit. In contrast, the controllable source electromagnetic method(CSEM) remains highly sensitive to resistivity changes. By simulating different CO_(2) plume migration conditions, we established the relevant models and calculated the corresponding electric field response characteristic curves, allowing us to analyze the CSEM's ability to monitor CO_(2) plumes. We considered potential scenarios for the migration and diffusion of offshore CO_(2) storage, including various burial depths, vertical extension diffusion, lateral extension diffusion,multiple combinations of lateral intervals, and electric field components. We also obtained differences in resistivity inversion imaging obtained by CSEM to evaluate its feasibility in monitoring and to analyze all the electric field(Ex, Ey, and Ez) response characteristics. CSEM has great potential in monitoring CO_(2) plume migration in offshore saltwater reservoirs due to its high sensitivity and accuracy. Furthermore, changes in electromagnetic field response reflect the transport status of CO_(2) plumes, providing an important basis for monitoring and evaluating CO_(2)transport behavior during storage processes.展开更多
基金supported by the 973 Program(Grant No 2007CB209600)Open Fund(No.GDL0706) of the Key Laboratory of Geo-detection(China University of Geosciences,Beijing),Ministry of Education
文摘The inversions of complex geophysical data always solve multi-parameter, nonlinear, and multimodal optimization problems. Searching for the optimal inversion solutions is similar to the social behavior observed in swarms such as birds and ants when searching for food. In this article, first the particle swarm optimization algorithm was described in detail, and ant colony algorithm improved. Then the methods were applied to three different kinds of geophysical inversion problems: (1) a linear problem which is sensitive to noise, (2) a synchronous inversion of linear and nonlinear problems, and (3) a nonlinear problem. The results validate their feasibility and efficiency. Compared with the conventional genetic algorithm and simulated annealing, they have the advantages of higher convergence speed and accuracy. Compared with the quasi-Newton method and Levenberg-Marquardt method, they work better with the ability to overcome the locally optimal solutions.
基金supported by the Robert and Maude Gledden Postgraduate Research Scholarship。
文摘3-D geochemical subsurface models,as constructed by spatial interpolation of drill-core assays,are valuable assets across multiple stages of the mineral industry's workflow.However,the accuracy of such models is limited by the spatial sparsity of the underlying drill-core,which samples only a small fraction of the subsurface.This limitation can be alleviated by integrating collocated 3-D models into the interpolation process,such as the 3-D rock property models produced by modern geophysical inversion procedures,provided that they are sufficiently resolved and correlated with the interpolation target.While standard machine learning algorithms are capable of predicting the target property given these data,incorporating spatial autocorrelation and anisotropy in these models is often not possible.We propose a Gaussian process regression model for 3-D geochemical interpolation,where custom kernels are introduced to integrate collocated 3-D rock property models while addressing the trade-off between the spatial proximity of drill-cores and the similarities in their collocated rock properties,as well as the relative degree to which each supporting 3-D model contributes to interpolation.The proposed model was evaluated for 3-D modelling of Mg content in the Kevitsa Ni-Cu-PGE deposit based on drill-core analyses and four 3-D geophysical inversion models.Incorporating the inversion models improved the regression model's likelihood(relative to a purely spatial Gaussian process regression model)when evaluated at held-out test holes,but only for moderate spatial scales(100 m).
文摘The paper presents a short overview about the application of joint inversion in geophysics. It gives also an alternative explanation for the term of “different data sets” and discusses what types of inversion procedures can be considered as joint inversion. Nowadays there are no standard standpoints using the appellation joint inversion. What is joint inversion? Based on the information matrix an answer could be given for this question what could be regarded as various types of data sets that are inverted simultaneously. We would like to expand the explanation—that is professed by many researchers—of the method that regards only the simultaneous inversion of data sets based on different physical parameters as joint inversion.
基金supported in part by the National Natural Science Foundation of China under Grant-in-Aid 40574053the Program for New Century Excellent Talents in University of China (NCET-06-0602)the National 973 Key Basic Research Development Program (No.2007CB209601)
文摘In this paper we discuss the edge-preserving regularization method in the reconstruction of physical parameters from geophysical data such as seismic and ground-penetrating radar data. In the regularization method a potential function of model parameters and its corresponding functions are introduced. This method is stable and able to preserve boundaries, and protect resolution. The effect of regularization depends to a great extent on the suitable choice of regularization parameters. The influence of the edge-preserving parameters on the reconstruction results is investigated and the relationship between the regularization parameters and the error of data is described.
基金financial support from the National Key R&D Program of China(2017YFC0601305)the China Geological Survey(DD20160125,DD20160207,DD20190010)the National Natural Science Foundation of China(41504076)。
文摘The Shuangjianzishan deposit in Inner Mongolia is a typical Ag-Pb-Zn deposit of the southern Great Xing’an Range.Proven reserves of Ag,Pb,and Zn in this deposit have reached the scale of super-large deposits,with favorable metallogenic conditions,strong prospecting signs,and high metallogenic potential.This paper reports a study involving integrated geophysical methods,including controlled-source audio-frequency magnetotelluric,gravity,magnetic,and shallow-seismic-reflection methods,to determine the spatial distribution of ore-controlling structures and subsurface intrusive rock for a depth range of<2000 m in the Shuangjianzishan ore district.The objective of this study is to construct a metallogenic model of the ore district and provide a scientific basis for the exploration of similar deposits in the deep and surrounding regions.We used three-dimensional inversion for controlled-source audio-frequency magnetotelluric data based on the limited memory quasi-Newton algorithm,and three-dimensional physical-property inversion for the gravity and magnetic data to obtain information about the subsurface distribution of ore-controlling structures and intrusive rocks.Under seismic reflection results,regional geology,petrophysical properties,and borehole information,the geophysical investigation shows that the Dashizhai group,which contains the main ore-bearing strata in the ore district,is distributed within a depth range of<1239 m,and is thick in the Xinglongshan ore block and the eastern part of the Shuangjianzishan ore block.The mineralization is spatially associated with a fault system characterized by NE-,NW-,and N-trending faults.The magnetic and electrical models identify large,deep bodies of intrusive rock that are inferred to have been involved in mineralization,with local shallow emplacement of some intrusions.Combining the subsurface spatial distributions of ore-bearing strata,ore-controlling faults,and intrusive rock,we propose two different metallogenic models for the Shuangjianzishan ore district,which provide a scientific basis for further prospecting in the deep regions of the ore district and surrounding areas.
文摘Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.
基金supported by the National Program on Key Basic Research Project of China(973 Program)(No.2013CB036002,No.2014CB046901)the National Major Scientific Equipment Developed Special Project(No.51327802)+3 种基金National Natural Science Foundation of China(No.51139004,No.41102183)the Research Fund for the Doctoral Program of Higher Education of China(No.20110131120070)Natural Science Foundation of Shandong Province(No.ZR2011EEQ013)the Graduate Innovation Fund of Shandong University(No.YZC12083)
文摘To minimize the number of solutions in 3D resistivity inversion, an inherent problem in inversion, the amount of data considered have to be large and prior constraints need to be applied. Geological and geophysical data regarding the extent of a geological anomaly are important prior information. We propose the use of shape constraints in 3D electrical resistivity inversion, Three weighted orthogonal vectors (a normal and two tangent vectors) were used to control the resistivity differences at the boundaries of the anomaly. The spatial shape of the anomaly and the constraints on the boundaries of the anomaly are thus established. We incorporated the spatial shape constraints in the objective function of the 3D resistivity inversion and constructed the 3D resistivity inversion equation with spatial shape constraints. Subsequently, we used numerical modeling based on prior spatial shape data to constrain the direction vectors and weights of the 3D resistivity inversion. We established a reasonable range between the direction vectors and weights, and verified the feasibility and effectiveness of using spatial shape prior constraints in reducing excessive structures and the number of solutions. We applied the prior spatially shape-constrained inversion method to locate the aquifer at the Guangzhou subway. The spatial shape constraints were taken from ground penetrating radar data. The inversion results for the location and shape of the aquifer agree well with drilling data, and the number of inversion solutions is significantly reduced.
基金Supported by Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (2019BT02H594)Sanya Technology Innovation Special Project (2022KJCX08)。
文摘Offshore carbon dioxide(CO_(2)) storage is an effective method for reducing greenhouse gas emissions. However, when using traditional seismic wave methods to monitor the migration of sequestration CO_(2) plumes, the characteristics of wave velocity changes tend to become insignificant beyond a certain limit. In contrast, the controllable source electromagnetic method(CSEM) remains highly sensitive to resistivity changes. By simulating different CO_(2) plume migration conditions, we established the relevant models and calculated the corresponding electric field response characteristic curves, allowing us to analyze the CSEM's ability to monitor CO_(2) plumes. We considered potential scenarios for the migration and diffusion of offshore CO_(2) storage, including various burial depths, vertical extension diffusion, lateral extension diffusion,multiple combinations of lateral intervals, and electric field components. We also obtained differences in resistivity inversion imaging obtained by CSEM to evaluate its feasibility in monitoring and to analyze all the electric field(Ex, Ey, and Ez) response characteristics. CSEM has great potential in monitoring CO_(2) plume migration in offshore saltwater reservoirs due to its high sensitivity and accuracy. Furthermore, changes in electromagnetic field response reflect the transport status of CO_(2) plumes, providing an important basis for monitoring and evaluating CO_(2)transport behavior during storage processes.