Structure-guided gravity inversion for layered density modeling with an application in the Chezhen Depression, Bohai Bay Basin

Geophysics ◽  
2021 ◽  
pp. 1-54
Author(s):  
Jie Liu ◽  
Jianzhong Zhang
Keyword(s):  
Spline Interpolation ◽  
A Priori ◽  
Structural Similarity ◽  
Gravity Inversion ◽  
Inversion Method ◽  
Bohai Bay ◽  
Bohai Bay Basin ◽  
Lower Paleozoic ◽  
Gas Targets ◽  
Density Modeling

Gravity inversion, as a static potential field inversion, has inherent ambiguity with low vertical resolution. In order to reduce the nonuniqueness of inversion, it is necessary to impose the apriori constraints derived by other geophysical inversion, drilling or geological modeling. Based on the a priori normalized gradients derived from seismic imaging or reference models, a structure-guided gravity inversion method with a few known point constraints is developed for mapping density with multiple layers. The cubic B-spline interpolation is used to parameterize the forward modeling calculation of the gravity response to smooth density fields. A recently proposed summative gradient is used to maximize the structural similarity between the a priori and inverted models. We first demonstrate the methodology, followed by a synthetic fault model example to confirm its validity. Monte Carlo tests and uncertainty tests further illustrate the stability and practicality of the method. This method is easy to implement, and consequently produces an interpretable density model with geological consistency. Finally, we apply this method to the density modeling of the Chezhen Depression in the Bohai Bay Basin. Our work determines the distribution of deep Lower Paleozoic carbonate rocks and Archean buried hills with high-density characteristics. Our results are consistent with the existing formation mechanism of the “upper source-lower reservoir” type oil-gas targets.

Three‐dimensional gravity inversion using simulated annealing: Constraints on the diapiric roots of allochthonous salt structures

Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1438-1449 ◽  
Author(s):  
Seiichi Nagihara ◽  
Stuart A. Hall
Keyword(s):  
Simulated Annealing ◽  
Gulf Of Mexico ◽  
Gravity Data ◽  
A Priori ◽  
Gravity Inversion ◽  
Inversion Method ◽  
A Priori Information ◽  
Final Density ◽  
Smoothing Effects ◽  
Priori Information

In the northern continental slope of the Gulf of Mexico, large oil and gas reservoirs are often found beneath sheetlike, allochthonous salt structures that are laterally extensive. Some of these salt structures retain their diapiric feeders or roots beneath them. These hidden roots are difficult to image seismically. In this study, we develop a method to locate and constrain the geometry of such roots through 3‐D inverse modeling of the gravity anomalies observed over the salt structures. This inversion method utilizes a priori information such as the upper surface topography of the salt, which can be delineated by a limited coverage of 2‐D seismic data; the sediment compaction curve in the region; and the continuity of the salt body. The inversion computation is based on the simulated annealing (SA) global optimization algorithm. The SA‐based gravity inversion has some advantages over the approach based on damped least‐squares inversion. It is computationally efficient, can solve underdetermined inverse problems, can more easily implement complex a priori information, and does not introduce smoothing effects in the final density structure model. We test this inversion method using synthetic gravity data for a type of salt geometry that is common among the allochthonous salt structures in the Gulf of Mexico and show that it is highly effective in constraining the diapiric root. We also show that carrying out multiple inversion runs helps reduce the uncertainty in the final density model.

Gravity inversion of a discontinuous relief stabilized by weighted smoothness constraints on depth

Geophysics ◽  
1999 ◽  
Vol 64 (5) ◽  
pp. 1429-1437 ◽  
Author(s):  
Valéria C. F. Barbosa ◽  
João B. C. Silva ◽  
Walter E. Medeiros
Keyword(s):  
A Priori ◽  
Bouguer Anomaly ◽  
Initial Solution ◽  
Maximum Depth ◽  
Normal Faults ◽  
Gravity Inversion ◽  
Inversion Method ◽  
Solution Stability ◽  
Flat Bottom ◽  
Northern Portion

We present a new stable gravity inversion method applied to the mapping of an interface separating two homogeneous media. In contrast with previous similar methods, it does not impose an overall smoothness on the estimated interface to stabilize the solution. The density contrast between the media is assumed to be known. The interpretation model for the upper medium consists of rectangular juxtaposed prisms whose thicknesses represent the depths to the interface and are the parameters to be estimated. The true interface is assumed to be flat everywhere except at faults. To incorporate this attribute into the estimated relief, we developed an iterative process in which three kinds of constraints are imposed on parameters: (1) proximity between values of adjacent parameters, (2) lower and upper bounds to parameters, and (3) proximity between the values of parameters and fixed numerical values. Starting with an initial solution which presents an overall smooth relief, the method enhances initially estimated geometric features of the interface; that is, flat areas will tend to become flatter and steep areas will tend to become steeper. This is accomplished by weighting the constraints, which requires proximity between adjacent parameters. The weights are initialized with values equal to unity and are updated automatically to enhance any discrepancy between adjacent depths that have been detected at the initial solution. Constraints 2 and 3 are used both to compensate for the decrease in solution stability caused by the introduction of small weights and to reinforce flatness at the basin bottom. Constraint 2 imposes that any depth be nonnegative and smaller than an a priori known maximum depth value, whereas constraint 3 imposes that all depths be closest to a value greater than the maximum depth. The trade‐off between these conflicting constraints is attained with a final relief presenting flat bottom and steep borders. The method was tested with a synthetic gravity anomaly produced by a simulated sedimentary cratonic extensional basin whose basement consists of steep edges and a flat bottom. The results showed an improvement in the resolution of the relief, leading to a reliable mapping both of the sharp discontinuities at the borders and of the lateral extent of the base of the basin. Additionally, the method produced excellent estimates for the average dip angles of the basin edges (presumably controlled by normal faults), indicating, in this way, its potential in interpreting data produced by this kind of basin. The method was applied to the Bouguer anomaly from the northern portion of Steptoe Valley, Nevada, delineating an isolated basin with a wider, flat base and relatively straight borders as compared with the estimate imposing overall smoothness on the relief.

Three-dimensional gravity inversion using graph theory to delineate the skeleton of homogeneous sources

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. G53-G66 ◽  
Author(s):  
Rodrigo Bijani ◽  
Cosme F. Ponte-Neto ◽  
Dionisio U. Carlos ◽  
Fernando J. S. Silva Dias
Keyword(s):  
Genetic Algorithm ◽  
Graph Theory ◽  
Minimum Spanning Tree ◽  
Gravity Data ◽  
A Priori ◽  
Three Dimensional ◽  
Real Data ◽  
Gravity Inversion ◽  
Inversion Method ◽  
Dimensional Gravity

We developed a new strategy, based on graph theory concepts, to invert gravity data using an ensemble of simple point masses. Our method consisted of a genetic algorithm with elitism to generate a set of possible solutions. Each estimate was associated to a graph to solve the minimum spanning tree (MST) problem. To produce unique and stable estimates, we restricted the position of the point masses by minimizing the statistical variance of the distances of an MST jointly with the data-misfit function during the iterations of the genetic algorithm. Hence, the 3D spatial distribution of the point masses identified the skeleton of homogeneous gravity sources. In addition, our method also gave an estimation of the anomalous mass of the source. So, together with the anomalous mass, the skeleton could aid other 3D methods with promising geometric a priori parameters. Several tests with different values of regularizing parameter were made to bespeak this new regularizing strategy. The inversion results applied to noise-corrupted synthetic gravity data revealed that, regardless of promising starting models, the estimated distribution of point masses and the anomalous mass offered valuable information about the homogeneous sources in the subsurface. Tests on real data from a portion of Quadrilátero Ferrífero, Minas Gerais state, Brazil, were performed for complementary analysis of the proposed inversion method.

Sign in / Sign up

Export Citation Format

Share Document