Three-dimensional constrained gravity inversion of Moho depth and crustal structural characteristics at Mozambique continental margin

Mozambique's continental margin in East Africa was formed during the break-off stage of the east and west Gondwana lands. Studying the geological structure and division of continent-ocean boundary(COB) in Mozambique's continental margin is considered of great significance to rebuild Gondwana land and understand its movement mode. Along these lines, in this work, the initial Moho was fit using the known Moho depth from reflection seismic profiles, and a 3D multi-point constrained gravity inversion was carried out. Thus, highaccuracy Moho depth and crustal thickness in the study area were acquired. According to the crustal structure distribution based on the inversion results, the continental crust at the narrowest position of the Mozambique Channel was detected. According to the analysis of the crustal thickness, the Mozambique ridge is generally oceanic crust and the COB of the whole Mozambique continental margin is divided.

Three-dimensional gravity inversion based on 3D U-Net++

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.

Extrapolated Tikhonov method and inversion of 3D density images of gravity data

Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.

A 3D Geological Model Constrained by Gravity and Magnetic Inversion and its Exploration Implications for the World-class Zhuxi Tungsten Deposit, South China

The Zhuxi tungsten deposit in Jiangxi Province,South China,contains a total W reserve of about 2.86 Mt at an average grade of 0.54 wt%WO3,representing the largest W deposit in the world.Numerous studies on the metallogeny of the deposit have included its timing,the ore-controlling structures and sedimentary host rocks and their implications for mineral exploration.However,the deep nappe structural style of Taqian-Fuchun metallogenic belt that hosts the W deposit,and the spatial shape and scale of deeply concealed intrusions and their sedimentary host rocks are still poorly defined,which seriously restricts the discovery of new deposits at depth and in surrounding areas of the W deposit.Modern 3 D geological modeling is an important tool for the exploration of concealed orebodies,especially in brownfield environments.There are obvious density contrast and weak magnetic contrast in the ore-controlling strata and granite at the periphery of the deposit,which lays a physical foundation for solving the 3 D spatial problems of the ore-controlling geological body in the deep part of the study area through gravity and magnetic modeling.Gravity data(1:50000)and aeromagnetic data(1:50000)from the latest geophysical surveys of 2016-2018 have been used,firstly,to carry out a potential field separation to obtain residual anomalies for gravity and magnetic interactive inversion.Then,on the basis of the analysis of the relationship between physical properties and lithology,under the constraints of surface geology and borehole data,human-computer interactive gravity and magnetic inversion for 18 cross-sections were completed.Finally,the 3 D geological model of the Zhuxi tungsten deposit and its periphery have been established through these 18 sections,and the spatial shape of the intrusions and strata with a depth of 5 km underground were obtained,initially realizing―transparency‖for ore-controlling bodies.According the analysis of the geophysical,geochemical,and geological characteristics of the Zhuxi tungsten deposit,we discern three principles for prospecting and prediction in the research area,and propose five new exploration targets in its periphery.

Optimization algorithm for rapid 3D gravity inversion

The practical application of 3D inversion of gravity data requires a lot of computation time and storage space.To solve this problem,we present an integrated optimization algorithm with the following components:(1)targeting high accuracy in the space domain and fast computation in the wavenumber domain,we design a fast 3D forward algorithm with high precision;and(2)taking advantage of the symmetry of the inversion matrix,the main calculation in gravity conjugate gradient inversion is decomposed into two forward calculations,thus optimizing the computational efficiency of 3D gravity inversion.We verify the calculation accuracy and efficiency of the optimization algorithm by testing various grid-number models through numerical simulation experiments.

3D inversion of borehole gravity data using cokriging

Borehole gravity has been used in mineral exploration recently with the advent of slim-hole gravimeters. It is logical to perform inversion to utilize the information in the newly acquired data. The inversions were carried out by using cokriging,which is a geostatistical method of estimation that minimizes the error variance by applying cross-correlation between several variables. In this study the estimated densities are derived by using boreholes gravity and known densities along the borehole. This method does not need iterative process and computes efficiently. The selection of examples demonstrates that this method has the ability to include physical property from borehole measurements in the inversion. The synthetic examples demonstrate the density variation along a borehole can be well determined without depth constraints in the inversion. The resolution of the recovered model can be further improved by including the densities along the borehole for inversion. However,this capability decreases dramatically with the increasing of distance between the anomalous body and the borehole.

3D Gravity Inversion with Correlation Image in Space Domain and Application to the Northern Sinai Peninsula

We present a 3D inversion method to recover density distribution from gravity data in space domain.Our method firstly employs 3D correlation image of the vertical gradient of gravity data as a starting model to generate a higher resolution image for inversion.The 3D density distribution is then obtained by inverting the correlation image of gravity data to fit the observed data based on classical inversion method of the steepest descent method.We also perform the effective equivalent storage and subdomain techniques in the starting model calculation,the forward modeling and the inversion procedures,which allow fast computation in space domain with reducing memory consumption but maintaining accuracy.The efficiency and stability of our method is demonstrated on two sets of synthetic data and one set of the Northern Sinai Peninsula gravity data.The inverted 3D density distributions show that high density bodies beneath Risan Aniza and low density bodies exist to the southeast of Risan Aniza at depths between 1~10 and 20 km,which may be originated from hot anomalies in the lower crust.The results show that our inversion method is useful for 3D quantitative interpretation.

Gravity inversion using the frequency characteristics of the density distribution

Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to search over the entire solution space for a more refined result. However, the inversion will be difficult with the increased parameters in the large search space and the number of computations increases exponentially. ｜n this paper, we propose a novel approach based on the frequency characteristics of the density distribution over the mesh. The purposes of our study are to reduce the parameters of three- dimensional gravity inversion and to lighten the image quality of the inversion result. The results show that the new method can expedite the inversion processing and get a better geological interpretation than tradition methods.

Inversion of gravity gradient data based on spatial gradient weighting

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.

Coastal bathymetry inversion using SAR-based altimetric gravity data:A case study over the South Sandwich Island

The global bathymetry models are usually of low accuracy over the coastline of polar areas due to the harsh climatic environment and the complex topography.Satellite altimetric gravity data can be a supplement and plays a key role in bathymetry modeling over these regions.The Synthetic Aperture Radar(SAR)altimeters in the missions like CryoSat-2 and Sentinel-3A/3B can relieve waveform contamination that existed in conventional altimeters and provide data with improved accuracy and spatial resolution.In this study,we investigate the potential application of SAR altimetric gravity data in enhancing coastal bathymetry,where the effects on local bathymetry modeling introduced from SAR altimetry data are quantified and evaluated.Furthermore,we study the effects on bathymetry modeling by using different scale factor calculation approaches,where a partition-wise scheme is implemented.The numerical experiment over the South Sandwich Islands near Antarctica suggests that using SARbased altimetric gravity data improves local coastal bathymetry modeling,compared with the model calculated without SAR altimetry data by a magnitude of 3:55 m within 10 km of offshore areas.Moreover,by using the partition-wise scheme for scale factor calculation,the quality of the coastal bathymetry model is improved by 7.34 m compared with the result derived from the traditional method.These results indicate the superiority of using SAR altimetry data in coastal bathymetry inversion.

3-D velocity structure in the central-eastern part of Qilianshan

The 3-D velocity tomography image of the central-eastern part of Qilianshan is obtained by the joint inversion of 3-D velocity structure and focal parameters based on the S-P data of micro-earthquakes recorded by the digital seismic network set up for a Sino-French cooperation program since 1996. The inversed velocity structure does primarily reflect some important features of the deep structure in the region and provide the scientific background for the further study of active tectonic structure and the calculation of earthquake parameters.

Earthquake relocation and 3-dimensional crustal structure of P-wave velocity in central-western China

A simultaneous inversion of earthquake relocation and three-dimensional crustal structure of P-wave velocity in central-western China (21癗~36癗, 98癊~112癊) were performed in this paper. The crustal P-wave velocity model and earthquake relocation for this region are obtained using Pg and Sg phase readings of 9 988 earthquakes from 1992 to 1999 recorded at 193 seismic stations within central-western China by SPHYPIT90 and SPHREL3D90 programs. A lateral inhomogeneous structure of P-wave velocity in this region was obtained. Ob-vious contrast of P-wave velocities was revealed on both sides of active fault zones. Relocated epicenters of 6 459 events show clear lineation along active faults, which indicated a close correlation between seismicity and the active faults in this region. Focal depths of 82% relocated events ranged from 0 to 20 km, which is in good agreement with that from double-difference earthquake location algorithm.

FAST APPROXIMATE INVERSION OF MAGNETOTELLURIC(FAIMT) DATA FOR A HORIZONTALLY LAYERED EARTH WITH NEAR-SURFACE INHOMOGENEITIES

This paper discusses use of approximations and the Integral Mean Value Theorem to show that 6 coefficients approximately describe the distortions of near surface inhomogeneities on the MT field of a horizontally layered earth model. When these 6 coefficients are considered together with those of the magnetic field of a horizontally layered earth model,the analytic and approximate wave impedance equations can be derived for the MT response of a horizontally layered earth model with near-surface 2-D and 3-D inhomogeneities. These approximate wave impedance equations are used with inverted MT data for 2-D and 3-D forward modelling. Although these 6 coefficients cannot be determined before inversion,initial estimates can be used. The 6 coefficients and the asistivity and thickness of each layer of a horizontally layered earth can be obtained by using published inversion methods. The 6 coefficients give important informaion (depths and resistivities) on the near-surface inhomogenelties.The authors inverted 2-D and 3-D theoretical models for Fast Approximate Inversion of Magnetotelluric (FAIMT) data for a horizontally layered earth with near-surface inhomogeneities compares favorably with traditional invrsion methods, especially for inverting regional or basin structures. This method simplifies computation and gives a reasonable 1 -D geological model with fewer nonuniquenas problems.

3-D seismic tomography for velocity and interface structure of the crust and upper mantle(theoreticalpart)

A method of three dimensional (3-D) model parameterization is presented that makes forward and inverse problems become easy. The velocity and interface structure of crust and upper mantle are described by a set of highly smoothed functions. Shooting ray tracing method is chosen to calculate the ray paths for both forward and inverse problems. The partial derivatives of traveltime with respect to parameters of the model grids are calculated analytically while rays are being traced. Because velocity and interface functions have second-order continuous partial derivatives, the geometrical shadow zones at the surface caused by scattering and focusing of ray paths can be prevented. After ray tracing, an equation consisting of matrix and vectors for inverse problem is obtained. We use singular value decomposition method with damped factor to solve the equation. A synthetic data set which consists of several in-line profiles is used to test the methods. The results show that the methods are robust. Compared with the two dimensional method, the 3-D inversion method can give the right position of interfaces and the velocity structure when the crustal model is complicated.

Reconstruction of Temperature Field in 3-D, Absorbing, Emitting, and Anisotropically Scattering Medium

The soft measurement technology of flame temperature field is an efficient method to learn the combustion status in furnace. Generally, it reconstructs the temperature field in furnace through the image of flame, which is a process to solve radiative inverse problem. In this paper, the flame of pulverized coal is considered as 3-D, absorbing, emitting, and anisotropically scattering non-gray medium. Through the study on inverse problem of radiative heat transfer, the temperature field in this kind of medium has been reconstructed. The mechanism of 3-D radiative heat transfer in a rectangular media, which is 2 m×3 m× 5 m and full of CO2, N2 and carbon particles, is studied with Monte Carlo method. The 3-D temperature field in this rectangular space is reconstructed and the influence of particles density profile is discussed.

Three-Dimensional Density Distribution and Seismic Activity along the Guxiang–Tongmai Segment of the Jiali Fault,Tibet

The Guxiang-Tongmai segment of the Jiali fault is situated northeast of the Namche Barwa Syntaxis in northeastern Tibet.It is one of the most active strike-slip faults near the syntaxis and plays a pivotal role in the examination of seismic activity within the eastern Himalayan Syntaxis.New study in the research region has yielded a 1:200000 gravity dataset covering an area 1500 km^(2).Using wavelet transform multiscale decomposition,scratch analysis techniques,and 3D gravity inversion methods,gravity anomalies,fault distributions,and density structures were determined across various scales.Through the integration of our new gravity data with other geophysical and geological information,our findings demonstrate substantial variations in the overall crustal density within the region,with the fault distribution closely linked to these density fluctuations.Disparities in stratigraphic density are important causes of variations in the capacity of geological formations to endure regional tectonic stress.Earthquakes are predominantly concentrated within the density transition zone and are primarily situated in regions of elevated density.The hanging wall stress within the Guxiang-Tongmai segment of the Jiali fault exhibits a notable concentration,marked by pronounced anisotropy,and is positioned within the density differential zone,which is prone to earthquakes.

Advanced Computational Methods of Rapid and Rigorous 3-D Inversion of Airborne Electromagnetic Data

We develop a new computational method for modeling and inverting frequency domain airborne electromagnetic(EM)data.Our method is based on the contraction integral equation method for forward EM modeling and on inversion using the localized quasi-linear(LQL)approximation followed by the rigorous inversion,if necessary.The LQL inversion serves to provide a fast image of the target.These results are checked by a rigorous update of the domain electric field,allowing a more accurate calculation of the predicted data.If the accuracy is poorer than desired,rigorous inversion follows,using the resulting conductivity distribution and electric field from LQL as a starting model.The rigorous inversion iteratively solves the field and domain equations,converting the non-linear inversion into a series of linear inversions.We test this method on synthetic and field data.The results of the inversion are very encouraging with respect to both the speed and the accuracy of the algorithm,showing this is a useful tool for airborne EM interpretation.

STUDY ON CRUSTAL THICKNESS IN TAIWAN REGION——PRELIMINARY RESULTS BY 3D GRAVITY INVERSION

Bouguer anomaly can be considered as the total sum of the attraction effects resulting from the geological bodies beneath the observational point with different scales, densities and depths. How can we extract tlae information of deep crust from Bouguer anomaly?The authors suggested using the formula of 3D positive calculation to determine the surface gravity anomaly caused by the shallow strata. After layering forward calculation, the deep residual gravity anomaly reflecting crustal structure is obtained. If we take Moho depth of 3 l km from DSS in Fujian Province as a controlling point, the crustal thickness throughout the whole region is inversed. As for the calculation and error, see Refs.

Geophysical Modeling with Satellite Gravity Data: Eigen-6C4 vs. GGM Plus

Satellite data sets are an asset in global gravity collections;their characteristics vary in coverage and resolution. New collections appear often, and the user must adapt fast to their characteristics. Their use in geophysical modeling is rapidly increasing;with this in mind we compare two of the most densely populated sets: EIGEN-6C4 and GGMplus. We characterize them in terms of their frequency histograms, Free Air anomalies, power spectrum, and simple Bouguer anomalies. The nature of the digital elevation models used for data reduction is discussed. We conclude that the GGMplus data set offers a better spatial resolution. To evaluate their effect in geophysical modelling, we chose an inland region with a prominent volcanic structure in which we perform 3D inversions of the respective Bouguer anomalies, obtaining density variations that in principle can be associated with the geologic materials and the structure of the volcanic edifice. Model results are analyzed along sections of the inverted data;we conclude that the GGMplus data set offers higher resolution in the cases analyzed.

A Two-Streamfunction-Coordinate Method for the Numerical Solution of 3-D Aero-Thermodynamics Direct and Inverse Problems

This paper describes a new method which has been developed for the solution of direct and inverse problems of 3-D compressible flows in turbomachinery.Two types of streamfunctions are proposed in the paper and the streamfunction-coordinate system is applied in numerical computations.The algorithm is applied to stator blades and the results are compared with experimental data,It is shown that the comparisons are very satis- factory.