We apply reweighted inversion focusing to full tensor gravity gradiometry data using message-passing interface （MPI） and compute unified device architecture （CUDA） parallel computing algorithms, and then combine M...We apply reweighted inversion focusing to full tensor gravity gradiometry data using message-passing interface （MPI） and compute unified device architecture （CUDA） parallel computing algorithms, and then combine MPI with CUDA to formulate a hybrid algorithm. Parallel computing performance metrics are introduced to analyze and compare the performance of the algorithms. We summarize the rules for the performance evaluation of parallel algorithms. We use model and real data from the Vinton salt dome to test the algorithms. We find good match between model efficiency and feasibility of parallel computing gravity gradiometry data. and real density data, and verify the high algorithms in the inversion of full tensor展开更多
The accuracy of the Earth＇s gravitational field measured from the gravity field and steady-state ocean circulation explorer（GOCE）,up to 250 degrees,influenced by the radial gravity gradient V zz and three-dimension...The accuracy of the Earth＇s gravitational field measured from the gravity field and steady-state ocean circulation explorer（GOCE）,up to 250 degrees,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij from the satellite gravity gradiometry（SGG） are contrastively demonstrated based on the analytical error model and numerical simulation,respectively.Firstly,the new analytical error model of the cumulative geoid height,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij are established,respectively.In 250 degrees,the GOCE cumulative geoid height error measured by the radial gravity gradient V zz is about 2 1/2 times higher than that measured by the three-dimensional gravity gradient V ij.Secondly,the Earth＇s gravitational field from GOCE completely up to 250 degrees is recovered using the radial gravity gradient V zz and three-dimensional gravity gradient V ij by numerical simulation,respectively.The study results show that when the measurement error of the gravity gradient is 3×10 12 /s 2,the cumulative geoid height errors using the radial gravity gradient V zz and three-dimensional gravity gradient V ij are 12.319 cm and 9.295 cm at 250 degrees,respectively.The accuracy of the cumulative geoid height using the three-dimensional gravity gradient V ij is improved by 30%-40% on average compared with that using the radial gravity gradient V zz in 250 degrees.Finally,by mutual verification of the analytical error model and numerical simulation,the orders of magnitude from the accuracies of the Earth＇s gravitational field recovery make no substantial differences based on the radial and three-dimensional gravity gradients,respectively.Therefore,it is feasible to develop in advance a radial cold-atom interferometric gradiometer with a measurement accuracy of 10 13 /s 2-10 15 /s 2 for precisely producing the next-generation GOCE Follow-On Earth gravity field model with a high spatial resolution.展开更多
The density inversion of gravity gradiometry data has attracted considerable attention;however,in large datasets,the multiplicity and low depth resolution as well as efficiency are constrained by time and computer mem...The density inversion of gravity gradiometry data has attracted considerable attention;however,in large datasets,the multiplicity and low depth resolution as well as efficiency are constrained by time and computer memory requirements.To solve these problems,we improve the reweighting focusing inversion and probability tomography inversion with joint multiple tensors and prior information constraints,and assess the inversion results,computing efficiency,and dataset size.A Message Passing Interface(MPI)-Open Multi-Processing(OpenMP)-Computed Unified Device Architecture(CUDA)multilevel hybrid parallel inversion,named Hybrinv for short,is proposed.Using model and real data from the Vinton Dome,we confirm that Hybrinv can be used to compute the density distribution.For data size of 100×100×20,the hybrid parallel algorithm is fast and based on the run time and scalability we infer that it can be used to process the large-scale data.展开更多
On the： basis of wavelet theory, we propose an outlier-detection algorithm for satellite gravity ometry by applying a wavelet-shrinkage-de-noising method to some simulation data with white noise and ers. The result S...On the： basis of wavelet theory, we propose an outlier-detection algorithm for satellite gravity ometry by applying a wavelet-shrinkage-de-noising method to some simulation data with white noise and ers. The result Shows that this novel algorithm has a 97% success rate in outlier identification and that be efficiently used for pre-processing real satellite gravity gradiometry data.展开更多
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.展开更多
At present, gravity field and steady-state ocean circulation explorer（GOCE） gravity data are always used to compute regional gravity anomaly and geoid height. In this study, the latest GOCE gravity field model data...At present, gravity field and steady-state ocean circulation explorer（GOCE） gravity data are always used to compute regional gravity anomaly and geoid height. In this study, the latest GOCE gravity field model data（from Oct. 2009 to Jul. 2010） are used to compute the gravity gradient of China's Mainland according to a rigorous recursion formula（in all the six directions）. The results show that the numerical values of the gravity gradients are larger in the T rr direction than those in the other directions. They reflect the terrain characteristics in detail and correlate with the regional tectonics; however, in the T ql and T r l directions,the numerical values are relatively smaller and the gravity gradients in the T r l direction do not reflect the terrain characteristics in detail.展开更多
With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processin...With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processing and interpretation of large-scale high-precision data, the use of the graphics processing unit process unit （GPU） and preconditioning methods are very important in the data inversion. In this paper, an improved preconditioned conjugate gradient algorithm is proposed by combining the symmetric successive over-relaxation （SSOR） technique and the incomplete Choleksy decomposition conjugate gradient algorithm （ICCG）. Since preparing the preconditioner requires extra time, a parallel implement based on GPU is proposed. The improved method is then applied in the inversion of noise- contaminated synthetic data to prove its adaptability in the inversion of 3D FTG data. Results show that the parallel SSOR-ICCG algorithm based on NVIDIA Tesla C2050 GPU achieves a speedup of approximately 25 times that of a serial program using a 2.0 GHz Central Processing Unit （CPU）. Real airbome gravity-gradiometry data from Vinton salt dome （south- west Louisiana, USA） are also considered. Good results are obtained, which verifies the efficiency and feasibility of the proposed parallel method in fast inversion of 3D FTG data.展开更多
Denoising of full-tensor gravity-gradiometer data involves detailed information from field sources, especially the data mixed with high-frequency random noise. We present a denoising method based on the translation-in...Denoising of full-tensor gravity-gradiometer data involves detailed information from field sources, especially the data mixed with high-frequency random noise. We present a denoising method based on the translation-invariant wavelet with mixed thresholding and adaptive threshold to remove the random noise and retain the data details. The novel mixed thresholding approach is devised to filter the random noise based on the energy distribution of the wavelet coefficients corresponding to the signal and random noise. The translation- invariant wavelet suppresses pseudo-Gibbs phenomena, and the mixed thresholding better separates the wavelet coefficients than traditional thresholding. Adaptive Bayesian threshold is used to process the wavelet coefficients according to the specific characteristics of the wavelet coefficients at each decomposition scale. A two-dimensional discrete wavelet transform is used to denoise gridded data for better computational efficiency. The results of denoising model and real data suggest that compared with Gaussian regional filter, the proposed method suppresses the white Gaussian noise and preserves the high-frequency information in gravity-gradiometer data. Satisfactory denoising is achieved with the translation-invariant wavelet.展开更多
Firstly, a new analytical error model of the cumulative geoid height using the three-dimensional diagonal tensors of satellite gravity gradiometry （SGG） is introduced based on the variance-covariance matrix principl...Firstly, a new analytical error model of the cumulative geoid height using the three-dimensional diagonal tensors of satellite gravity gradiometry （SGG） is introduced based on the variance-covariance matrix principle. Secondly, a study for the requirements demonstration on the next-generation GOCE Follow-On satellite gravity gradiometry system is developed using different satellite orbital altitudes and measurement accuracies of satellite gravity gradiometer by the new analytical error model of SGG. The research results show that it is preferable to design satellite orbital altitudes of 300 km–400km and choose the measurement accuracies of 10-13/s2 –10-15/s2 from satellite gravity gradiometer. Finally, the complementarity of the four-stage satellite gravity missions, including past CHAMP, current GRACE, and GOCE, and next-generation GOCE Follow-On, is contrastively demonstrated for precisely recovering the Earth’s full-frequency gravitational field with high spatial resolution.展开更多
基金supported by the Sino-Probe09(No.201011078)National High-tech R&D Program(No.863 and2014AA06A613)
文摘We apply reweighted inversion focusing to full tensor gravity gradiometry data using message-passing interface （MPI） and compute unified device architecture （CUDA） parallel computing algorithms, and then combine MPI with CUDA to formulate a hybrid algorithm. Parallel computing performance metrics are introduced to analyze and compare the performance of the algorithms. We summarize the rules for the performance evaluation of parallel algorithms. We use model and real data from the Vinton salt dome to test the algorithms. We find good match between model efficiency and feasibility of parallel computing gravity gradiometry data. and real density data, and verify the high algorithms in the inversion of full tensor
基金Project supported by the Main Direction Program of Knowledge Innovation of the Chinese Academy of Sciences for Distinguished Young Scholars (Grant No. KZCX2-EW-QN114)the National Natural Science Foundation of China for Young Scholars (GrantNos. 41004006,41131067,and 11173049)+3 种基金the Merit-Based Scientific Research Foundation of the State Ministry of Human Resources and Social Security of China for Returned Overseas Chinese Scholars (Grant No. 2011)the Open Research Fund Programof the Key Laboratory of Computational Geodynamics of the Chinese Academy of Sciences (Grant No. 2011-04)the Frontier Field Program of Knowledge Innovation of Institute of Geodesy and Geophysics of the Chinese Academy of Sciencesthe Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Grant No. PLN1113)
文摘The accuracy of the Earth＇s gravitational field measured from the gravity field and steady-state ocean circulation explorer（GOCE）,up to 250 degrees,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij from the satellite gravity gradiometry（SGG） are contrastively demonstrated based on the analytical error model and numerical simulation,respectively.Firstly,the new analytical error model of the cumulative geoid height,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij are established,respectively.In 250 degrees,the GOCE cumulative geoid height error measured by the radial gravity gradient V zz is about 2 1/2 times higher than that measured by the three-dimensional gravity gradient V ij.Secondly,the Earth＇s gravitational field from GOCE completely up to 250 degrees is recovered using the radial gravity gradient V zz and three-dimensional gravity gradient V ij by numerical simulation,respectively.The study results show that when the measurement error of the gravity gradient is 3×10 12 /s 2,the cumulative geoid height errors using the radial gravity gradient V zz and three-dimensional gravity gradient V ij are 12.319 cm and 9.295 cm at 250 degrees,respectively.The accuracy of the cumulative geoid height using the three-dimensional gravity gradient V ij is improved by 30%-40% on average compared with that using the radial gravity gradient V zz in 250 degrees.Finally,by mutual verification of the analytical error model and numerical simulation,the orders of magnitude from the accuracies of the Earth＇s gravitational field recovery make no substantial differences based on the radial and three-dimensional gravity gradients,respectively.Therefore,it is feasible to develop in advance a radial cold-atom interferometric gradiometer with a measurement accuracy of 10 13 /s 2-10 15 /s 2 for precisely producing the next-generation GOCE Follow-On Earth gravity field model with a high spatial resolution.
基金support by the China Postdoctoral Science Foundation(2017M621151)Northeastern University Postdoctoral Science Foundation(20180313)+1 种基金the Fundamental Research Funds for Central Universities(N180104020)NSFCShandong Joint Fund of the National Natural Science Foundation of China(U1806208)
文摘The density inversion of gravity gradiometry data has attracted considerable attention;however,in large datasets,the multiplicity and low depth resolution as well as efficiency are constrained by time and computer memory requirements.To solve these problems,we improve the reweighting focusing inversion and probability tomography inversion with joint multiple tensors and prior information constraints,and assess the inversion results,computing efficiency,and dataset size.A Message Passing Interface(MPI)-Open Multi-Processing(OpenMP)-Computed Unified Device Architecture(CUDA)multilevel hybrid parallel inversion,named Hybrinv for short,is proposed.Using model and real data from the Vinton Dome,we confirm that Hybrinv can be used to compute the density distribution.For data size of 100×100×20,the hybrid parallel algorithm is fast and based on the run time and scalability we infer that it can be used to process the large-scale data.
基金supported by the Director Foundation of the Institute of Seismology,China Earthquake Administration (IS201126025)The Basis Research Foundation of Key laboratory of Geospace Environment & Geodesy Ministry of Education,China (10-01-09)
文摘On the： basis of wavelet theory, we propose an outlier-detection algorithm for satellite gravity ometry by applying a wavelet-shrinkage-de-noising method to some simulation data with white noise and ers. The result Shows that this novel algorithm has a 97% success rate in outlier identification and that be efficiently used for pre-processing real satellite gravity gradiometry data.
基金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.
基金supported by Key Projects of Henan Province Department of Education Science and Technology(14B420001)
文摘At present, gravity field and steady-state ocean circulation explorer（GOCE） gravity data are always used to compute regional gravity anomaly and geoid height. In this study, the latest GOCE gravity field model data（from Oct. 2009 to Jul. 2010） are used to compute the gravity gradient of China's Mainland according to a rigorous recursion formula（in all the six directions）. The results show that the numerical values of the gravity gradients are larger in the T rr direction than those in the other directions. They reflect the terrain characteristics in detail and correlate with the regional tectonics; however, in the T ql and T r l directions,the numerical values are relatively smaller and the gravity gradients in the T r l direction do not reflect the terrain characteristics in detail.
基金the Sub-project of National Science and Technology Major Project of China(No.2016ZX05027-002-003)the National Natural Science Foundation of China(No.41404089)+1 种基金the State Key Program of National Natural Science of China(No.41430322)the National Basic Research Program of China(973 Program)(No.2015CB45300)
文摘With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processing and interpretation of large-scale high-precision data, the use of the graphics processing unit process unit （GPU） and preconditioning methods are very important in the data inversion. In this paper, an improved preconditioned conjugate gradient algorithm is proposed by combining the symmetric successive over-relaxation （SSOR） technique and the incomplete Choleksy decomposition conjugate gradient algorithm （ICCG）. Since preparing the preconditioner requires extra time, a parallel implement based on GPU is proposed. The improved method is then applied in the inversion of noise- contaminated synthetic data to prove its adaptability in the inversion of 3D FTG data. Results show that the parallel SSOR-ICCG algorithm based on NVIDIA Tesla C2050 GPU achieves a speedup of approximately 25 times that of a serial program using a 2.0 GHz Central Processing Unit （CPU）. Real airbome gravity-gradiometry data from Vinton salt dome （south- west Louisiana, USA） are also considered. Good results are obtained, which verifies the efficiency and feasibility of the proposed parallel method in fast inversion of 3D FTG data.
基金supported by the National Key Research and Development Plan Issue(Nos.2017YFC0602203 and2017YFC0601606)the National Science and Technology Major Project Task(No.2016ZX05027-002-003)+4 种基金the National Natural Science Foundation of China(Nos.41604089 and 41404089)the State Key Program of National Natural Science of China(No.41430322)the Marine/Airborne Gravimeter Research Project(No.2011YQ12004505)the State Key Laboratory of Marine Geology,Tongji University(No.MGK1610)the Basic Scientific Research Business Special Fund Project of Second Institute of Oceanography,State Oceanic Administration(No.14275-10)
文摘Denoising of full-tensor gravity-gradiometer data involves detailed information from field sources, especially the data mixed with high-frequency random noise. We present a denoising method based on the translation-invariant wavelet with mixed thresholding and adaptive threshold to remove the random noise and retain the data details. The novel mixed thresholding approach is devised to filter the random noise based on the energy distribution of the wavelet coefficients corresponding to the signal and random noise. The translation- invariant wavelet suppresses pseudo-Gibbs phenomena, and the mixed thresholding better separates the wavelet coefficients than traditional thresholding. Adaptive Bayesian threshold is used to process the wavelet coefficients according to the specific characteristics of the wavelet coefficients at each decomposition scale. A two-dimensional discrete wavelet transform is used to denoise gridded data for better computational efficiency. The results of denoising model and real data suggest that compared with Gaussian regional filter, the proposed method suppresses the white Gaussian noise and preserves the high-frequency information in gravity-gradiometer data. Satisfactory denoising is achieved with the translation-invariant wavelet.
基金Project supported by the Main Direction Program of Knowledge Innovation of Chinese Academy of Sciences for Distinguished Young Scholar (Grant No. KZCX2-EW-QN114)the National Natural Science Foundation of China for Young Scholar (Grant Nos. 41004006, 41131067, 11173049, and 41202094)+5 种基金the Merit-based Scientific Research Foundation of the State Ministry of Human Resources and Social Security of China for Returned Overseas Chinese Scholars(Grant No. 2011)the Open Research Fund Program of the Key Laboratory of Computational Geodynamics of Chinese Academy of Sciences (Grant No. 2011-04)the Open Research Fund Program of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, China (Grant No. 11-01-02)the Open Research Fund Program of the Key Laboratory of Geo-Informatics of National Administration of Surveying, Mapping and Geoinformation of China(Grant No. 201322)the Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Grant No. PLN1113)the Foundation of State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing (Grant No. PRP/open-1206)
文摘Firstly, a new analytical error model of the cumulative geoid height using the three-dimensional diagonal tensors of satellite gravity gradiometry （SGG） is introduced based on the variance-covariance matrix principle. Secondly, a study for the requirements demonstration on the next-generation GOCE Follow-On satellite gravity gradiometry system is developed using different satellite orbital altitudes and measurement accuracies of satellite gravity gradiometer by the new analytical error model of SGG. The research results show that it is preferable to design satellite orbital altitudes of 300 km–400km and choose the measurement accuracies of 10-13/s2 –10-15/s2 from satellite gravity gradiometer. Finally, the complementarity of the four-stage satellite gravity missions, including past CHAMP, current GRACE, and GOCE, and next-generation GOCE Follow-On, is contrastively demonstrated for precisely recovering the Earth’s full-frequency gravitational field with high spatial resolution.
