Along-curved-surface gravitational field inversion method and device based on target layer

Abstract

The invention relates to the technical field of geophysical exploration, in particular to an along-curved-surface gravitational field inversion method and device based on a target layer, and the method comprises the steps: taking a target layer curved surface of a target region as a constraint, peeling off a shallow-layer residual gravity anomaly field in a Bouguer gravity anomaly field of the target region layer by layer, and subtracting a region gravity anomaly field below the target layer to obtain a gravity anomaly field; obtaining a residual gravity anomaly field of the target layer curved surface; and in combination with the fluctuation degree of the curved surface of the target layer, performing inversion on the residual gravity anomaly field of the curved surface of the target layer by using deconvolution to obtain a density superposition field of the curved surface of the target layer. According to the method, a new along-curved-surface gravity field inversion method is established, depth parameters of a target layer curved surface are considered during inversion calculation, interference signal influences are eliminated, a density field signal along the curved surface is obtained, a density superposition field of the target layer curved surface is obtained, and an oil-gas-water enrichment area is depicted for follow-up density change interface prediction and oil-gas-water enrichment area depicting. And a basis is provided for providing reference suggestions for well position deployment and production regulation and control.

Description

Technical Field

[0001] The invention relates to the technical field of geophysical detection, and is a method and device for inverting a gravity field along a curved surface based on a target layer. Background Art

[0002] Microgravity detection technology is one of the most important means in geophysical exploration. It can be used to accurately and reliably predict the density change interface of gas-water, oil and gas, and define the relative enrichment range of oil and gas in oil and gas reservoirs. It is of great significance for adjusting well site deployment and ensuring the normal and stable production of oil and gas wells.

[0003] Although the existing microgravity detection methods have made great progress, the current gravity inversion is still an underdetermined problem due to the lack of longitudinal parameter constraints, and the multi-solution phenomenon is serious. The density inversion technology can reduce the uncontrollability caused by depth variables and transform the three-dimensional inversion into a two-dimensional inversion, but the current density inversion technology can only achieve density inversion at a certain depth, and cannot eliminate the influence of interference signals when facing a target layer with complex and changeable structural fluctuations. Summary of the invention

[0004] The present invention provides a method and device for inversion of the gravity field along a curved surface based on a target layer, which overcomes the shortcomings of the above-mentioned prior art and can effectively solve the problem that the existing density inversion method can only realize density inversion at one depth and cannot eliminate the influence of interference signals when facing a target layer with complex and changeable structural fluctuations.

[0005] One of the technical solutions of the present invention is achieved by the following measures: a method for inverting the gravity field along a curved surface based on a target layer, comprising:

[0006] Taking the target layer surface of the target area as a constraint, the shallow residual gravity anomaly field is peeled off layer by layer in the Bouguer gravity anomaly field of the target area, and the regional gravity anomaly field below the target layer is subtracted to obtain the residual gravity anomaly field of the target layer surface;

[0007] Combined with the undulation of the target layer surface, the residual gravity anomaly field of the target layer surface is inverted by using deconvolution to obtain the density stacking field of the target layer surface.

[0008] The following are further optimizations and / or improvements to the above technical solutions:

[0009] The above-mentioned obtaining of the Bouguer gravity anomaly field of the target area includes:

[0010] Obtain the original observed gravity anomaly field of the target area;

[0011] The original observed gravity anomaly field is corrected to obtain the Bouguer gravity anomaly field of the target area, where the corrections include solid tide correction, zero drift correction, mesosphere correction, latitude correction and terrain correction.

[0012] The above method takes the target layer surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface, including:

[0013] Expand the Bouguer gravity anomaly field in the target area;

[0014] The multi-scale surface fitting method is used to fit each sampling point in the Bouguer gravity anomaly field to obtain the corresponding separation scale.

[0015] The depth information of the target layer surface is converted into a corresponding scale, and the separation scale of each sampling point is updated based on the converted scale to obtain the abnormal field of the depth of the target layer surface;

[0016] The shallow residual gravity anomaly field is peeled off layer by layer, and the regional gravity anomaly field below the target layer is subtracted to obtain the residual gravity anomaly field of the target layer surface.

[0017] In the above, combined with the undulation of the target layer surface, the residual gravity anomaly field of the target layer surface is inverted by deconvolution to obtain the density stacking field of the target layer surface. The inversion model is as follows:

[0018]

[0019] Among them, h i,j is the layer depth of point (i, j) on the target layer surface.

[0020] The second technical solution of the present invention is achieved by the following measures: a device for inverting the gravity field along the curved surface based on the target layer, comprising:

[0021] The first processing unit, taking the target layer surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface;

[0022] The second processing unit, in combination with the undulation degree of the target layer surface, uses deconvolution to invert the residual gravity anomaly field of the target layer surface to obtain the density superposition field of the target layer surface.

[0023] The following are further optimizations and / or improvements to the above technical solutions:

[0024] The first processing unit comprises:

[0025] The acquisition module acquires the Bouguer gravity anomaly field of the target area, including:

[0026] Obtain the original observed gravity anomaly field of the target area;

[0027] The original observed gravity anomaly field is corrected to obtain the Bouguer gravity anomaly field of the target area, where the correction includes solid tide correction, zero drift correction, mesosphere correction, latitude correction and terrain correction;

[0028] The separation module takes the target layer surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface, including:

[0029] Expand the Bouguer gravity anomaly field in the target area;

[0030] The multi-scale surface fitting method is used to fit each sampling point in the Bouguer gravity anomaly field to obtain the corresponding separation scale.

[0031] The depth information of the target layer surface is converted into a corresponding scale, and the separation scale of each sampling point is updated based on the converted scale to obtain the abnormal field of the depth of the target layer surface;

[0032] The shallow residual gravity anomaly field is peeled off layer by layer, and the regional gravity anomaly field below the target layer is subtracted to obtain the residual gravity anomaly field of the target layer surface.

[0033] The second processing unit comprises:

[0034] The inversion model building module uses deconvolution and depth parameters to build an inversion model, where the inversion model is as follows:

[0035]

[0036] Among them, h i,j is the layer depth of point (i, j) on the target layer surface;

[0037] The inversion module, combined with the undulation of the target layer surface, uses the inversion model to invert the residual gravity anomaly field of the target layer surface to obtain the density stacking field of the target layer surface.

[0038] The present invention establishes a new method for inverting the gravity field along the curved surface. When performing the inversion calculation, the depth parameter of the target layer curved surface is considered to obtain the density field signal along the curved surface, and the density superposition field of the target layer curved surface is obtained, which provides a basis for subsequent density change interface prediction, delineation of oil, gas and water enrichment areas, and reference suggestions for well site deployment and production regulation. In addition, the present invention uses the target layer curved surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer curved surface. The generated separation result has good directional selectivity and weaker depth aliasing, effectively presents interference signals, and the separation result is conducive to identifying abnormal bodies, providing key information for obtaining the density superposition field of the target layer curved surface in the later stage. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Attached Figure 1 It is a schematic diagram of the method flow of the present invention.

[0040] Attached Figure 2 It is a schematic flow chart of the method for separating the residual gravity anomaly field in the present invention.

[0041] Attached Figure 3 This is a schematic diagram of the top structure of a section of the Toutunhe Formation in Example 3 of the present invention.

[0042] Attached Figure 4 This is a schematic diagram of the superimposed density field at the top of a section of the Toutunhe Formation in Example 3 of the present invention.

[0043] Attached Figure 5 This is a local detail diagram of the layer density inversion field of the Dixi 121 well area in Example 3 of the present invention.

[0044] Attached Figure 6 This is the logging curve diagram of the DXHW1212 well in Example 3 of the present invention.

[0045] Attached Figure 7 This is the logging curve diagram of the DXHW1215 well in Example 3 of the present invention.

[0046] Attached Figure 8 The figure is a schematic diagram of the structure of a device of the present invention.

[0047] Attached Fig. 9 This is a schematic diagram of another device structure of the present invention. DETAILED DESCRIPTION

[0048] The present invention is not limited by the following embodiments, and specific implementation methods can be determined based on the technical solution of the present invention and actual conditions.

[0049] The present invention will be further described below in conjunction with embodiments and drawings:

[0050] Embodiment 1: As attached Figure 1 As shown, the embodiment of the present invention discloses a method for inverting a gravity field along a curved surface based on a target layer, comprising:

[0051] Step S110, taking the target layer surface of the target area as a constraint, peeling off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracting the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface;

[0052] Step S120, combining the undulation degree of the target layer surface, using deconvolution to invert the residual gravity anomaly field of the target layer surface, to obtain the density superposition field of the target layer surface.

[0053] The present invention establishes a new method for inverting the gravity field along the curved surface. When performing the inversion calculation, the depth parameter of the target layer curved surface is considered to obtain the density field signal along the curved surface, and the density superposition field of the target layer curved surface is obtained, which provides a basis for subsequent density change interface prediction, delineation of oil, gas and water enrichment areas, and reference suggestions for well site deployment and production regulation. In addition, the present invention uses the target layer curved surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer curved surface. The generated separation result has good directional selectivity and weaker depth aliasing, effectively presents interference signals, and the separation result is conducive to identifying abnormal bodies, providing key information for obtaining the density superposition field of the target layer curved surface in the later stage.

[0054] In the above step S120, the residual gravity anomaly field of the target layer surface is inverted by using deconvolution in combination with the undulation of the target layer surface to obtain the density stacking field of the target layer surface. The inversion of the residual gravity anomaly field of the target layer surface by using deconvolution is calculated by using an inversion model. The inversion model is as follows:

[0055]

[0056] Among them, h i,j is the layer depth of point (i, j) on the target layer surface.

[0057] In actual geological conditions, it is very rare for the strata to be nearly horizontal, and the target layer structure often has certain fluctuations or is even complex and changeable. Therefore, even if the horizontal density stacking field depth is located near the target layer, there is a large amount of non-target layer interference information in the field. The present invention updates the depth parameter in the layer density inversion method in real time during the calculation inversion, changing the original depth constant z into a depth variable z i,jThe gravitational acceleration is the convolution of density and Green's function, so it can be simplified to: Δg(x,y,z)=G·r(x,y,z)*σ(x,y,z). Assuming that in a very thin curved layer, the density does not change with depth, Δg can be obtained from the Fourier transform. h (x,y)=G·r h (x,y)*σ h (x, y), the layer gravity value is the convolution of the layer density and the Green function, and the calculation formula for the layer density inversion is: σ h (x,y)=G·r h -1 (x,y)*Δg h (x,y), obtained by Fourier transform When the layer thickness is very small, the Green's function changes approximately unchanged, which can be approximated by the above inversion model h i,j is the layer depth of point (i, j) on the target surface.

[0058] Embodiment 2: As attached Figure 1 , 2 As shown, the embodiment of the present invention discloses a method for inverting the gravity field along the curved surface based on the target layer, which is a further optimization of the above embodiment, wherein the curved surface of the target layer in the target area is used as a constraint, the shallow residual gravity anomaly field is peeled off layer by layer in the Bouguer gravity anomaly field of the target area, and the regional gravity anomaly field below the target layer is subtracted to obtain the residual gravity anomaly field of the curved surface of the target layer, including:

[0059] Step S210, obtaining the Bouguer gravity anomaly field of the target area, includes:

[0060] Step S211, obtaining the original observed gravity anomaly field of the target area;

[0061] Step S212, correcting the original observed gravity anomaly field to obtain the Bouguer gravity anomaly field of the target area, wherein the correction includes solid tide correction, zero drift correction, intermediate layer correction, latitude correction and terrain correction;

[0062] Since the original observed gravity anomaly field is affected by factors related to time and spatial position, it is necessary to calculate various correction items to eliminate interference before further analysis and research. Therefore, the time-related correction items of the original observed gravity anomaly field include but are not limited to solid tide correction and zero drift correction, and the spatial position-related correction items include but are not limited to intermediate layer correction, terrain correction and dimension correction.

[0063] Step S220, taking the target layer surface of the target area as a constraint, peeling off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracting the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface, as shown in the attached figure. Figure 2 As shown, including:

[0064] Step S221, performing edge expansion processing on the Bouguer gravity anomaly field in the target area;

[0065] Step S222, using a multi-scale surface fitting method to fit each sampling point in the Bouguer gravity anomaly field to obtain a corresponding separation scale;

[0066] When performing anomaly separation for a certain sampling point (i, j) as mentioned above, it is necessary to fit and calculate a total of 25 sampling points from (is*2, js*2) to (i+s*2, j+s*2), where s is the scale. After fitting, it is found that the scale is in multiples of the depth. In the conventional multi-scale surface fitting method, the scale of each layer of anomaly separation is a constant. After optimization, the new layer separation technology can be converted into the corresponding scale according to the depth of a specific surface, and the separation scale of each sampling point is updated to obtain the anomaly field of a specific surface depth;

[0067] Step S223, converting the depth information of the target layer surface into a corresponding scale, and updating the separation scale of each sampling point in combination with the converted scale to obtain an abnormal field of the depth of the target layer surface;

[0068] Step S224, peel off the shallow residual gravity anomaly field layer by layer, and subtract the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface.

[0069] Embodiment 3, attached Figure 3 Based on the top structural depth signal of the first section of the Toutunhe Formation shown in the figure, the gravity field inversion calculation along the curved surface is performed. In order to concentrate the inversion signal on the target layer, the depth parameter selection needs to be updated in real time during the calculation using the method proposed by the present invention, and the original depth constant z is changed to the depth variable z. i,j , calculating the layer density value at a specific point will greatly reduce the computational workload and improve computational efficiency.

[0070] In Fourier transform:

[0071]

[0072] in, and is the depth of the upper and lower interfaces of point (i, j) on the target layer surface. When the layer thickness is very small, the Green function changes approximately unchanged and can be approximated as:

[0073]

[0074] h i,j is the layer depth of point (i, j) on the target surface.

[0075] The stacked density field at the top of the first section of the Toutunhe Formation is obtained by inversion calculation based on the above formula, as shown in the attached figure. Figure 4 As shown, its physical meaning is the superimposed density effect within the extremely thin range of the surface.

[0076] A total of 8 vertical wells and 2 horizontal wells in the study area were used for verification, including Dixi 123, Dixi 128, Dixi 2, Dixi 125, Dixi 149, Dixi 121, Dixi 129, Dixi 146, and DXHW1211 and DXHW1212. According to the well logging data, the average density of the gas-bearing layer was read and correlated with the layer density inversion results. The average density of the well logging was correlated with the layer density field reading, and the linear trend was basically consistent, indicating that the layer density inversion results have high accuracy and the consistency is 90%. The specific verification well consistency statistics are shown in Table 1:

[0077] Table 1 Statistical table of the consistency of the verification wells

[0078]

[0079] Taking Dixi 129 as an example, the local details of the layer density stacking field of Dixi 121 well area obtained by the method disclosed in the present invention are shown in the attached figure. Figure 5 As shown in the figure, the logging curves of DXHW1212 well and DXHW1215 well in Dixi 121 well are shown in the attached figure. Figure 6 , 7 As shown. Figure 5 It can be found that in the density superposition field along the curved surface, the entire section of the DXHW1212 well is in a low-density area, the gas abundance in the sand body is large, and the gas production capacity is strong. There is a local high-density area in the middle section of the DXHW1215 well, with low gas abundance and average production capacity. Figure 6 and attached Figure 7 According to the logging interpretation results, DXHW1212 and DXHW1215 wells are interpreted as gas-bearing layers throughout the entire section, and the hydrocarbon content in the middle of the well is relatively small, which is consistent with the prediction results. Therefore, the method disclosed in the present invention is effective and meets the accuracy requirements.

[0080] Embodiment 4: As attached Figure 8 As shown, the embodiment of the present invention discloses a device for inverting a gravity field along a curved surface based on a target layer, comprising:

[0081] The first processing unit, taking the target layer surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface;

[0082] The second processing unit, in combination with the undulation degree of the target layer surface, uses deconvolution to invert the residual gravity anomaly field of the target layer surface to obtain the density superposition field of the target layer surface.

[0083] Embodiment 5: As attached Fig. 9 As shown, the embodiment of the present invention discloses a device for inverting a gravity field along a curved surface based on a target layer, comprising:

[0084] The first processing unit comprises:

[0085] The acquisition module acquires the Bouguer gravity anomaly field of the target area, including:

[0086] Obtain the original observed gravity anomaly field of the target area;

[0087] The original observed gravity anomaly field is corrected to obtain the Bouguer gravity anomaly field of the target area, where the correction includes solid tide correction, zero drift correction, mesosphere correction, latitude correction and terrain correction;

[0088] The separation module takes the target layer surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface, including:

[0089] Expand the Bouguer gravity anomaly field in the target area;

[0090] The multi-scale surface fitting method is used to fit each sampling point in the Bouguer gravity anomaly field to obtain the corresponding separation scale.

[0091] The depth information of the target layer surface is converted into a corresponding scale, and the separation scale of each sampling point is updated based on the converted scale to obtain the abnormal field of the depth of the target layer surface;

[0092] The shallow residual gravity anomaly field is peeled off layer by layer, and the regional gravity anomaly field below the target layer is subtracted to obtain the residual gravity anomaly field of the target layer surface.

[0093] The second processing unit comprises:

[0094] The inversion model building module uses deconvolution and depth parameters to build an inversion model, where the inversion model is as follows:

[0095]

[0096] Among them, h i,j is the layer depth of point (i, j) on the target layer surface;

[0097] The inversion module, combined with the undulation of the target layer surface, uses the inversion model to invert the residual gravity anomaly field of the target layer surface to obtain the density stacking field of the target layer surface.

[0098] Embodiment 6: The embodiment of the present invention discloses a storage medium, on which is stored a computer program that can be read by a computer, and the computer program is configured to execute a method for inverting a gravity field along a curved surface based on a target layer when running.

[0099] The above storage medium may include, but is not limited to: a USB flash drive, a read-only memory, a mobile hard disk, a magnetic disk or an optical disk, and other media that can store computer programs.

[0100] Embodiment 7: The embodiment of the present invention discloses an electronic device, including a processor and a memory, wherein the memory stores a computer program, and the computer program is loaded and executed by the processor to implement a method for inverting a gravity field along a curved surface based on a target layer.

[0101] The processor may be a central processing unit (CPU), a general purpose processor, a digital signal processor (DSP), an ASIC, an FPGA or other programmable logic device, a transistor logic device, a hardware component or any combination thereof. It may implement or execute various exemplary logic blocks, modules and circuits described in conjunction with the disclosure of this application. It may also be a combination that implements computing functions, such as a combination of one or more microprocessors, a combination of a DSP and a microprocessor, and the like. The memory may include, but is not limited to, various media that can store computer programs, such as a USB flash drive, a read-only memory, a mobile hard disk, a magnetic disk or an optical disk.

[0102] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of complete hardware embodiments, complete software embodiments, or embodiments in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code. The scheme in the embodiments of the present application can be implemented in various computer languages, for example, object-oriented programming language Java and literal scripting language JavaScript, etc.

[0103] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0104] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0105] The above technical features constitute the best embodiment of the present invention, which has strong adaptability and best implementation effect. Non-essential technical features can be added or reduced according to actual needs to meet the requirements of different situations.

Claims

1. A method for inverting the gravity field along a curved surface based on a target layer, characterized in that: include: Taking the target layer surface of the target area as a constraint, the shallow residual gravity anomaly field is peeled off layer by layer in the Bouguer gravity anomaly field of the target area, and the regional gravity anomaly field below the target layer is subtracted to obtain the residual gravity anomaly field of the target layer surface; Combined with the undulation of the target layer surface, the residual gravity anomaly field of the target layer surface is inverted by using deconvolution to obtain the density stacking field of the target layer surface.

2. The method for inverting the gravity field along a curved surface based on the target layer according to claim 1, characterized in that: The step of obtaining the Bouguer gravity anomaly field of the target area includes: Obtain the original observed gravity anomaly field of the target area; The original observed gravity anomaly field is corrected to obtain the Bouguer gravity anomaly field of the target area, where the corrections include solid tide correction, zero drift correction, mesosphere correction, latitude correction and terrain correction.

3. The method for inverting the gravity field along a curved surface based on the target layer according to claim 1 or 2, characterized in that: The method uses the target layer surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface, including: Expand the Bouguer gravity anomaly field in the target area; The multi-scale surface fitting method is used to fit each sampling point in the Bouguer gravity anomaly field to obtain the corresponding separation scale. The depth information of the target layer surface is converted into a corresponding scale, and the separation scale of each sampling point is updated based on the converted scale to obtain the abnormal field of the depth of the target layer surface; The shallow residual gravity anomaly field is peeled off layer by layer, and the regional gravity anomaly field below the target layer is subtracted to obtain the residual gravity anomaly field of the target layer surface.

4. The method for inverting the gravity field along a curved surface based on the target layer according to claim 1 or 2, characterized in that: The residual gravity anomaly field of the target layer surface is inverted by using deconvolution in combination with the undulation of the target layer surface to obtain the density stacking field of the target layer surface, wherein the inversion model is as follows: Among them, h i,j is the layer depth of point (i, j) on the target layer surface.

5. The method for inverting the gravity field along a curved surface based on the target layer according to claim 3, characterized in that: The residual gravity anomaly field of the target layer surface is inverted by using deconvolution in combination with the undulation degree of the target layer surface to obtain the density superposition field of the target layer surface, wherein the inversion calculation formula is as follows: Among them, h i,j is the layer depth of point (i, j) on the target layer surface.

6. A device for inverting a gravity field along a curved surface based on a target layer using the method as claimed in any one of claims 1 to 5, characterized in that: include: The first processing unit, taking the target layer surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface; The second processing unit, in combination with the undulation degree of the target layer surface, uses deconvolution to invert the residual gravity anomaly field of the target layer surface to obtain the density superposition field of the target layer surface.

7. The device for inverting the gravity field along a curved surface based on the target layer according to claim 6, characterized in that: The first processing unit comprises: The acquisition module acquires the Bouguer gravity anomaly field of the target area, including: Obtain the original observed gravity anomaly field of the target area; The original observed gravity anomaly field is corrected to obtain the Bouguer gravity anomaly field of the target area, where the correction includes solid tide correction, zero drift correction, mesosphere correction, latitude correction and terrain correction; The separation module takes the target layer surface of the target area as a constraint, peels off the shallow residual gravity anomaly field layer by layer in the Bouguer gravity anomaly field of the target area, and subtracts the regional gravity anomaly field below the target layer to obtain the residual gravity anomaly field of the target layer surface, including: Expand the Bouguer gravity anomaly field in the target area; The multi-scale surface fitting method is used to fit each sampling point in the Bouguer gravity anomaly field to obtain the corresponding separation scale. The depth information of the target layer surface is converted into a corresponding scale, and the separation scale of each sampling point is updated based on the converted scale to obtain the abnormal field of the depth of the target layer surface; The shallow residual gravity anomaly field is peeled off layer by layer, and the regional gravity anomaly field below the target layer is subtracted to obtain the residual gravity anomaly field of the target layer surface.

8. The device for inverting the gravity field along a curved surface based on the target layer according to claim 6 or 7, characterized in that: The second processing unit comprises: The inversion model building module uses deconvolution and depth parameters to build an inversion model, where the inversion model is as follows: Among them, h i,j is the layer depth of point (i, j) on the target layer surface; The inversion module, combined with the undulation of the target layer surface, uses the inversion model to invert the residual gravity anomaly field of the target layer surface to obtain the density stacking field of the target layer surface.

9. A storage medium, characterized in that: The storage medium stores a computer program that can be read by a computer, and the computer program is configured to execute the method for inverting the gravity field along a curved surface based on a target layer as described in any one of claims 1 to 5 when running.

10. An electronic device, characterized in that: It comprises a processor and a memory, wherein a computer program is stored in the memory, and the computer program is loaded and executed by the processor to implement the method for inverting the gravity field along the curved surface based on the target layer as described in any one of claims 1 to 5.

Images