Forward Modeling Research of Gravity Data Offset in North China Region

Due to the complexity and overlap effects of gravity anomalies, the gravity high point always deviates from the real position of the target. Researching the offset characteristics is helpful to determine the right location of gravity anomalies. In this paper, a series of forward models were designed to calculate the offset features under different assumptions. Different density models were established based on the measured gravity data in north China Sunhu area and a variety of methods were tried to eliminate offset effects under different conditions. The results indicate that the gravity anomalies of different density bodies can be separated effectively and the real position of target can be correctly located by using the “layer stripping method”.


Introduction
The Earth's gravity field is closely related to the density spatial distribution of material in the inner Earth.The information of gravity field can provide the basic parameter of resources exploration [1] [2] [3].So far, most of the oil and gas fields discovered in China are based on gravity method.And the application effects of gravity data heavily depend on the potential data process.The filtering effects on gravity and magnetic data operations was first studied by Nettleton (1954), who compared several different approaches in dealing with regional fields [4].S. Treitel (1967) used the Wiener filtering technique to divide the potential field [5].Wei (1984) designed a directional ellipse filter that could identify the local and regional fields with different trends [6].More recently, Guo et al.Open Journal of Yangtze Gas and Oil (2012) presented the preferential filtering method for de-noising and anomaly separation of gravity data based on Green's equivalent-layer concept and Wiener filter [7].
A lot of researches on the filtering methods have been made by the scholars both in China and abroad.However, the research on the characteristics of gravity high point offset and its correction methods is still less.In fact, it is common that the gravity high point deviates from the target's real position due to the complexity and overlap effects of gravity anomalies.Xiao et al. (1984) proposed that the location of the geological body would be affected by the regional fields and offset in the actual inversion process [8].
It is important to study the causes of the offset and find appropriate correction methods to make the geological body to its real location.In this paper, on the basis of the real data in North China Sunhu area, some different models were established, then the offset rules of gravity high point under different factors were researched in detail, and finally some methods were carried out to eliminate the offset effect.Among them, the geological modeling is a promoting method to find the shape and location of underground structures.Onur Osman, A.
Muhittin Albora et al. (2007,2016) performed Forced Neural Network (FNN) method to find the parameters of a model, and then applied it to synthetic and real data, which acquired a very good result [9] [10].

Potential Field Separation Method
In this study, two kinds of anomaly separation methods were applied-Regularization filtering and Layer stripping-to separate the anomalies.Regularization filtering method uses regularized stabilization filter factor to perform low-pass filtering on gravity and magnetic data.The filter parameters can be measured directly from the original anomaly profile or the plane contour map [11].
Suppose the measured gravity anomaly is ( ) and its spectrum is ( ) , F u v .The spectrum needs to be multiplied by a "regularization factor" ( ( ) where u and v denote angular frequencies along x and y, respectively.
λ is the maximum wavelength to be filtered out, and x λ = (line number in measuring area − 1) × line distance.
For example, the gravity anomaly component in the frequency range of can be obtained by Formulas (2) and (3): where, 01 02 λ and 02 λ are the maximum wavelengths corresponding to 01 f and 02 f , respectively.Open Journal of Yangtze Gas and Oil The thought of "Layer stripping method" is to calculate the anomalies caused by the known field sources, and then subtract them from the observed anomalies, which make the separation of remaining anomalies easier, so as to obtain the target anomalies [12].

The Characteristics of Gravity High Point Offset
In order to investigate the characteristics of the gravity anomaly offset, forward calculation has been carried out.Based on the seismic and well logging data (Table 1), a series of geological models of buried hills with different sizes, depths, quantities as well as density differences (Figure 1) are designed to calculate their gravity anomalies.In the Cartesian coordinate system, the coordinate origin is located directly above the surface of the target body.The x direction is horizontal to the right, the y direction is perpendicular to the paper, and the Z direction is perpendicular to the xoy plane.Model I (Figure 1) consists of 2 buried hills, which from left to right is named as buried hill A and B (B is target body), respectively.The density difference and coordinates center is Figure 2 and Figure 3 show that the gravity high point is shifted from the target (shown in the black box in Figure 3) due to the effect of the buried hill A.
The offset will lead to the misjudging of the real location of the target in the actual data process, so it is important to study the causes of the deviation and eliminate the offset effect.In this model, three factors-buried depth, density difference and horizontal distance of density bodies, which maybe the cause of the offset, are studied in detail.And the results are shown in Figures 4-6.It can be illustrated that, with the increase of the depth of buried hill A, the offset of the gravity high point gradually weakens (Figure 4).With the increase of the density difference of buried hill A, the offset gradually strengthens (Figure 5).And the offset gradually strengthens with the increase of the horizontal distance separation of two anomalies.But when the distance increases to a certain extent, the target's gravity high point will appear again, and correspond to the actual location of the geological body (Figure 6).

Elimination of Gravity High Point Offset
In order to eliminate the offset, the models of three buried hills (the model Ⅱ and model Ⅲ) are designed and calculated.In model II (Figure 7), the hills from left to right are buried hill A', B', C' respectively and B' is the target body.Similarly, the Cartesian coordinate system and the same notation are used in Section 3.1.The parameters used in buried hill A' are In model II, the "Layer stripping method" is used to separate the anomaly field, and the "known anomaly" is obtained by Regularization Filtering method, which has been introduced above, then subtract them (?) from the total anomalies.Figure 7 shows that the gravity response of the buried hill B' (target) and C' has been highlighted after the anomaly separation (point line).The gravity high point corresponds roughly to the actual position of geological body.
The model III is one for three buried hills with different sizes, depths and density differences.The parameters used in buried hill A" are

Figure 1 .
Figure 1.Schematic diagram of model I.

Figure 2 .
Figure 2. Schematic diagram of model I and its gravity response.

Figure 3 .
Figure 3. Anomaly contour map of model I.

Figure 4 .
Figure 4.The gravity high point offset in different buried depth of the buried hill A.

Figure 5 .
Figure 5.The gravity high point offset in different density of the buried hill A.

Figure 6 .
Figure 6.The offset of the target body gravity high point with different horizontal distance.

Figure 7 .
Figure 7.The original and processed gravity responses of model Ⅱ (upper) and the corresponding schematic geology model (lower).

..
The parameters in buried hill C" are set to be Open Journal of Yangtze Gas and Oil Model III is demonstrated in Figure8.The gravity high point of the buried hill B'' is shift to the left (full line).The buried hill C" of small size and shallow burial only shows a small protrusion on the response curve.Similarly, the "Layer Stripping method" and Regularization Filtering are used to deal with abnormal data.After this procedure is repeated twice, the offset phenomenon has been significantly weakened (point line), while the local response by buried hill C" has been effectively eliminated.A large-scale anomaly of the buried hill A" has suppressed as well.The anomaly of the target is correctly extracted out.

Figure 8 .
Figure 8.The original and processed gravity responses of model III.

Table 1 .
Logging density statistics for North China area.
20x Open Journal of Yangtze Gas and Oil