We combine gradient data from the Macao Science Satellite-1(MSS-1),CHAllenging Minisatellite Payload(CHAMP),Swarm-A,and Swarm-C satellites to develop a 110-degree lithospheric magnetic field model.We then comprehensiv...We combine gradient data from the Macao Science Satellite-1(MSS-1),CHAllenging Minisatellite Payload(CHAMP),Swarm-A,and Swarm-C satellites to develop a 110-degree lithospheric magnetic field model.We then comprehensively evaluate the performance of the model by power spectral comparisons,correlation analyses,sensitivity matrix assessments,and comparisons with existing lithospheric field models.Results showed that using near east–west gradient data from MSS-1 significantly enhances the model correlation in the spherical harmonic degree(N) range of 45–60 while also mitigating the decline in correlation at higher degrees(N > 60).Furthermore,the unique orbital characteristics of MSS-1 enable its gradient data to provide substantial contributions to modeling in the mid-to low-latitude regions.With continued data acquisition from MSS-1 and further optimization of data processing methods,the performance of the model is expected to improve.展开更多
Aeromagnetic gradient data needs to be reduced to the pole so that it can be better applied to geological interpretation through theoretical derivation.In this paper,we conduct research on the morphological characteri...Aeromagnetic gradient data needs to be reduced to the pole so that it can be better applied to geological interpretation through theoretical derivation.In this paper,we conduct research on the morphological characteristics of the total and horizontal gradient modules before and after reduction to the pole and design models at different latitudes,with consistent and inconsistent magnetic field direction and geological body magnetization direction.We discuss how to use the total gradient module and horizontal gradient module in geological interpretation.The reduced-to-the-pole（RTP） method is required for the horizontal gradient module method but not for the total gradient module.Finally,the conclusions derived from the theoretical models are verified through analysis of real data.The position determination of a geological body using the total gradient method,gradient data,or total-field data works better without RTP,ensuring data primitive authenticity.However,the horizontal gradient module should be reduced to the pole to determine the boundary of the geological body.Finally,the theoretical model is verified by actual data analysis.Both the total and horizontal gradient methods can be applied to geological interpretation.展开更多
Joint inversion is one of the most effective methods for reducing non-uniqueness for geophysical inversion.The current joint inversion methods can be divided into the structural consistency constraint and petrophysica...Joint inversion is one of the most effective methods for reducing non-uniqueness for geophysical inversion.The current joint inversion methods can be divided into the structural consistency constraint and petrophysical consistency constraint methods,which are mutually independent.Currently,there is a need for joint inversion methods that can comprehensively consider the structural consistency constraints and petrophysical consistency constraints.This paper develops the structural similarity index(SSIM)as a new structural and petrophysical consistency constraint for the joint inversion of gravity and vertical gradient data.The SSIM constraint is in the form of a fraction,which may have analytical singularities.Therefore,converting the fractional form to the subtractive form can solve the problem of analytic singularity and finally form a modified structural consistency index of the joint inversion,which enhances the stability of the SSIM constraint applied to the joint inversion.Compared to the reconstructed results from the cross-gradient inversion,the proposed method presents good performance and stability.The SSIM algorithm is a new joint inversion method for petrophysical and structural constraints.It can promote the consistency of the recovered models from the distribution and the structure of the physical property values.Then,applications to synthetic data illustrate that the algorithm proposed in this paper can well process the synthetic data and acquire good reconstructed results.展开更多
Compared with traditional gravity measurement data,gravity gradient tensor data contain more high frequency information,which can be used to understand the earth's interior structure,mineral resources distribution...Compared with traditional gravity measurement data,gravity gradient tensor data contain more high frequency information,which can be used to understand the earth's interior structure,mineral resources distribution etc. In this study,the authors present an algorithm for inverting gravity gradiometer data to recover the three-dimensional( 3-D) distributions of density. Spatial gradient weighting was used to constrain the extent of the body horizontally and vertically. A more accurate inversion result can be obtained by combining the prior information into the weighting function and applying it in inversion. This method was tested on synthetic models and the inverted results showed that the resolution was significantly improved. Moreover,the algorithm was applied to the inversion of empirical data from a salt dome located in Texas,USA,which demonstrated the validity of the proposed method.展开更多
The cross-gradients joint inversion technique has been applied to multiple geophysical data with a significant improvement on compatibility, but its numerical implementation for practical use is rarely discussed in th...The cross-gradients joint inversion technique has been applied to multiple geophysical data with a significant improvement on compatibility, but its numerical implementation for practical use is rarely discussed in the literature. We present a MATLAB-based three-dimensional cross-gradients joint inversion program with application to gravity and magnetic data. The input and output information was examined with care to create a rational, independent design of a graphical user interface (GUI) and computing kernel. For 3D visualization and data file operations, UBC-GIF tools are invoked using a series of I/O functions. Some key issues regarding the iterative joint inversion algorithm are also discussed: for instance, the forward difference of cross gradients, and matrix pseudo inverse computation. A synthetic example is employed to illustrate the whole process. Joint and separate inversions can be performed flexibly by switching the inversion mode. The resulting density model and susceptibility model demonstrate the correctness of the proposed program.展开更多
The processing and interpretation of gravity and gradient data plays an important role in geophysics.The cross gradient joint inversion is usually used for achieving structure coupling of multiple geophysical models. ...The processing and interpretation of gravity and gradient data plays an important role in geophysics.The cross gradient joint inversion is usually used for achieving structure coupling of multiple geophysical models. In order to realize the coupling of gravity and gravity tensor data,the authors analyzed each component.The results show that different types of data contain different direction information,and derived the joint inversion based on cross gradient function and applied it to model data. The theoretical model results show that the cross gradient method can reduce the multi solution and significantly improve the resolution of the inversion.The method was also applied to inverse the gravity tensor data in Vinton salt dome,showing that this method can get higher resolution results than the separate linear inversion,and be closer to the real density from drilling data.展开更多
With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion a...With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.展开更多
This study aims to improve knowledge of the structure of southwest Cameroon based on the analysis and interpretation of gravity data derived from the SGG-UGM-2 model. A residual anomaly map was first calculated from t...This study aims to improve knowledge of the structure of southwest Cameroon based on the analysis and interpretation of gravity data derived from the SGG-UGM-2 model. A residual anomaly map was first calculated from the Bouguer anomaly map, which is strongly affected by a regional gradient. The residual anomaly map generated provides information on the variation in subsurface density, but does not provide sufficient information, hence the interest in using filtering with the aim of highlighting the structures affecting the area of south-west Cameroon. Three interpretation methods were used: vertical gradient, horizontal gradient coupled with upward continuation and Euler deconvolution. The application of these treatments enabled us to map a large number of gravimetric lineaments materializing density discontinuities. These lineaments are organized along main preferential directions: NW-SE, NNE-SSW, ENE-WSW and secondary directions: NNW-SSE, NE-SW, NS and E-W. Euler solutions indicate depths of up to 7337 m. Thanks to the results of this research, significant information has been acquired, contributing to a deeper understanding of the structural composition of the study area. The resulting structural map vividly illustrates the major tectonic events that shaped the geological framework of the study area. It also serves as a guide for prospecting subsurface resources (water and hydrocarbons). .展开更多
The distribution of sampling data influences completeness of rule base so that extrapolating missing rules is very difficult. Based on data mining, a self-learning method is developed for identifying fuzzy model and e...The distribution of sampling data influences completeness of rule base so that extrapolating missing rules is very difficult. Based on data mining, a self-learning method is developed for identifying fuzzy model and extrapolating missing rules, by means of confidence measure and the improved gradient descent method. The proposed approach can not only identify fuzzy model, update its parameters and determine optimal output fuzzy sets simultaneously, but also resolve the uncontrollable problem led by the regions that data do not cover. The simulation results show the effectiveness and accuracy of the proposed approach with the classical truck backer-upper control problem verifying.展开更多
Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse ...Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse matrix factorization, combined with nine-point finite-difference for the linear system of equations. In the waveform inversion, we use preconditioned gradient method where the preconditioner is provided by the diagonal of the approximate Hessian matrix. We successfully applied waveform inversion method from low to high frequency in two sets of Marmousi data. One is the data set generated by frequencydomain finite-difference modeling, and the other is the original Marmousi shots data set. The former result is very close to the true velocity model. In the original shots data set inversion, we replace the prior source with estimated source; the result is also acceptable, and consistent with the true model.展开更多
On the basis of Space-Wise Least Square method, three numerical methods including Cholesky de- composition, pre-conditioned conjugate gradient and Open Multi-Processing parallel algorithm are applied into the determin...On the basis of Space-Wise Least Square method, three numerical methods including Cholesky de- composition, pre-conditioned conjugate gradient and Open Multi-Processing parallel algorithm are applied into the determination of gravity field with satellite gravity gradiometry data. The results show that, Cholesky de- composition method has been unable to meet the requirements of computation efficiency when the computer hardware is limited. Pre-conditioned conjugate gradient method can improve the computation efficiency of huge matrix inversion, but it also brings a certain loss of precision. The application of Open Multi-Processing parallel algorithm could achieve a good compromise between accuracy and computation efficiency.展开更多
基金the support of the National Natural Science Foundation of China (Nos. 42250103, 41974073, and 41404053)the Macao Foundation and the preresearch project of Civil Aerospace Technologies (Nos. D020308 and D020303)funded by China’s National Space Administration, the Specialized Research Fund for State Key Laboratories。
文摘We combine gradient data from the Macao Science Satellite-1(MSS-1),CHAllenging Minisatellite Payload(CHAMP),Swarm-A,and Swarm-C satellites to develop a 110-degree lithospheric magnetic field model.We then comprehensively evaluate the performance of the model by power spectral comparisons,correlation analyses,sensitivity matrix assessments,and comparisons with existing lithospheric field models.Results showed that using near east–west gradient data from MSS-1 significantly enhances the model correlation in the spherical harmonic degree(N) range of 45–60 while also mitigating the decline in correlation at higher degrees(N > 60).Furthermore,the unique orbital characteristics of MSS-1 enable its gradient data to provide substantial contributions to modeling in the mid-to low-latitude regions.With continued data acquisition from MSS-1 and further optimization of data processing methods,the performance of the model is expected to improve.
基金surpported by the National 863 Program(Grant No.2013AA063901)
文摘Aeromagnetic gradient data needs to be reduced to the pole so that it can be better applied to geological interpretation through theoretical derivation.In this paper,we conduct research on the morphological characteristics of the total and horizontal gradient modules before and after reduction to the pole and design models at different latitudes,with consistent and inconsistent magnetic field direction and geological body magnetization direction.We discuss how to use the total gradient module and horizontal gradient module in geological interpretation.The reduced-to-the-pole（RTP） method is required for the horizontal gradient module method but not for the total gradient module.Finally,the conclusions derived from the theoretical models are verified through analysis of real data.The position determination of a geological body using the total gradient method,gradient data,or total-field data works better without RTP,ensuring data primitive authenticity.However,the horizontal gradient module should be reduced to the pole to determine the boundary of the geological body.Finally,the theoretical model is verified by actual data analysis.Both the total and horizontal gradient methods can be applied to geological interpretation.
基金supported by the National Key Research and Development Program(Grant No.2021YFA0716100)the National Key Research and Development Program of China Project(Grant No.2018YFC0603502)+1 种基金the Henan Youth Science Fund Program(Grant No.212300410105)the provincial key R&D and promotion special project of Henan Province(Grant No.222102320279).
文摘Joint inversion is one of the most effective methods for reducing non-uniqueness for geophysical inversion.The current joint inversion methods can be divided into the structural consistency constraint and petrophysical consistency constraint methods,which are mutually independent.Currently,there is a need for joint inversion methods that can comprehensively consider the structural consistency constraints and petrophysical consistency constraints.This paper develops the structural similarity index(SSIM)as a new structural and petrophysical consistency constraint for the joint inversion of gravity and vertical gradient data.The SSIM constraint is in the form of a fraction,which may have analytical singularities.Therefore,converting the fractional form to the subtractive form can solve the problem of analytic singularity and finally form a modified structural consistency index of the joint inversion,which enhances the stability of the SSIM constraint applied to the joint inversion.Compared to the reconstructed results from the cross-gradient inversion,the proposed method presents good performance and stability.The SSIM algorithm is a new joint inversion method for petrophysical and structural constraints.It can promote the consistency of the recovered models from the distribution and the structure of the physical property values.Then,applications to synthetic data illustrate that the algorithm proposed in this paper can well process the synthetic data and acquire good reconstructed results.
基金Supported by Project of Natural Science Fund of Jilin Province(No.20180101312JC)
文摘Compared with traditional gravity measurement data,gravity gradient tensor data contain more high frequency information,which can be used to understand the earth's interior structure,mineral resources distribution etc. In this study,the authors present an algorithm for inverting gravity gradiometer data to recover the three-dimensional( 3-D) distributions of density. Spatial gradient weighting was used to constrain the extent of the body horizontally and vertically. A more accurate inversion result can be obtained by combining the prior information into the weighting function and applying it in inversion. This method was tested on synthetic models and the inverted results showed that the resolution was significantly improved. Moreover,the algorithm was applied to the inversion of empirical data from a salt dome located in Texas,USA,which demonstrated the validity of the proposed method.
文摘The cross-gradients joint inversion technique has been applied to multiple geophysical data with a significant improvement on compatibility, but its numerical implementation for practical use is rarely discussed in the literature. We present a MATLAB-based three-dimensional cross-gradients joint inversion program with application to gravity and magnetic data. The input and output information was examined with care to create a rational, independent design of a graphical user interface (GUI) and computing kernel. For 3D visualization and data file operations, UBC-GIF tools are invoked using a series of I/O functions. Some key issues regarding the iterative joint inversion algorithm are also discussed: for instance, the forward difference of cross gradients, and matrix pseudo inverse computation. A synthetic example is employed to illustrate the whole process. Joint and separate inversions can be performed flexibly by switching the inversion mode. The resulting density model and susceptibility model demonstrate the correctness of the proposed program.
基金Supported by Project of National Key Research and Development Plan(No.2017YFC0601606,2017YFC0602203)National Science and Technology Major Project(No.2016ZX05027-002-03)+1 种基金National Natural Science Foundation of China(No.41604098,41404089) State Key Program of National Natural Science of China(No.41430322)
文摘The processing and interpretation of gravity and gradient data plays an important role in geophysics.The cross gradient joint inversion is usually used for achieving structure coupling of multiple geophysical models. In order to realize the coupling of gravity and gravity tensor data,the authors analyzed each component.The results show that different types of data contain different direction information,and derived the joint inversion based on cross gradient function and applied it to model data. The theoretical model results show that the cross gradient method can reduce the multi solution and significantly improve the resolution of the inversion.The method was also applied to inverse the gravity tensor data in Vinton salt dome,showing that this method can get higher resolution results than the separate linear inversion,and be closer to the real density from drilling data.
基金partially supported by the National Natural Science Foundation of China (No.41230318)
文摘With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.
文摘This study aims to improve knowledge of the structure of southwest Cameroon based on the analysis and interpretation of gravity data derived from the SGG-UGM-2 model. A residual anomaly map was first calculated from the Bouguer anomaly map, which is strongly affected by a regional gradient. The residual anomaly map generated provides information on the variation in subsurface density, but does not provide sufficient information, hence the interest in using filtering with the aim of highlighting the structures affecting the area of south-west Cameroon. Three interpretation methods were used: vertical gradient, horizontal gradient coupled with upward continuation and Euler deconvolution. The application of these treatments enabled us to map a large number of gravimetric lineaments materializing density discontinuities. These lineaments are organized along main preferential directions: NW-SE, NNE-SSW, ENE-WSW and secondary directions: NNW-SSE, NE-SW, NS and E-W. Euler solutions indicate depths of up to 7337 m. Thanks to the results of this research, significant information has been acquired, contributing to a deeper understanding of the structural composition of the study area. The resulting structural map vividly illustrates the major tectonic events that shaped the geological framework of the study area. It also serves as a guide for prospecting subsurface resources (water and hydrocarbons). .
基金This project was supported by State Science &Technology Pursuing Project (2001BA204B01) of China and Foundation forUniversity Key Teacher by the Ministry of Education of China.
文摘The distribution of sampling data influences completeness of rule base so that extrapolating missing rules is very difficult. Based on data mining, a self-learning method is developed for identifying fuzzy model and extrapolating missing rules, by means of confidence measure and the improved gradient descent method. The proposed approach can not only identify fuzzy model, update its parameters and determine optimal output fuzzy sets simultaneously, but also resolve the uncontrollable problem led by the regions that data do not cover. The simulation results show the effectiveness and accuracy of the proposed approach with the classical truck backer-upper control problem verifying.
基金Supported by the National Natural Science Foundation of China (69983005)
文摘Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse matrix factorization, combined with nine-point finite-difference for the linear system of equations. In the waveform inversion, we use preconditioned gradient method where the preconditioner is provided by the diagonal of the approximate Hessian matrix. We successfully applied waveform inversion method from low to high frequency in two sets of Marmousi data. One is the data set generated by frequencydomain finite-difference modeling, and the other is the original Marmousi shots data set. The former result is very close to the true velocity model. In the original shots data set inversion, we replace the prior source with estimated source; the result is also acceptable, and consistent with the true model.
基金supproted by the National Natural Science Foundation of China(40874012,40904003,40974016,41004007)
文摘On the basis of Space-Wise Least Square method, three numerical methods including Cholesky de- composition, pre-conditioned conjugate gradient and Open Multi-Processing parallel algorithm are applied into the determination of gravity field with satellite gravity gradiometry data. The results show that, Cholesky de- composition method has been unable to meet the requirements of computation efficiency when the computer hardware is limited. Pre-conditioned conjugate gradient method can improve the computation efficiency of huge matrix inversion, but it also brings a certain loss of precision. The application of Open Multi-Processing parallel algorithm could achieve a good compromise between accuracy and computation efficiency.
