Gravity gradient data three-dimensional inversion method based on partial smoothness constraints

Abstract

The invention relates to a gravity gradient data three-dimensional inversion method based on partial smooth constraints. The method comprises the following steps of firstly, calculating second-order partial derivatives of the gravitational potential of a cuboid unit in x, y and z directions of a three-dimensional space to obtain a forward modeling formula of a cuboid element, giving a model, and performing forward modeling calculation according to the forward modeling formula to obtain observation data; manually assigning values to the roughness constraint matrix according to certain prior information to obtain a roughness matrix, and overcoming a skin effect by adopting a weighting mode of integral sensitivity; and finally, writing a target equation containing parameters such as an integral sensitivity matrix, a roughness matrix, regularization parameters, observation data and the like, and solving the target equation by adopting a conjugate gradient method to obtain an inversion result. The technology is suitable for three-dimensional gravity and gravity gradient data inversion under certain geological or geophysical prior conditions.

Description

technical field [0001] The invention relates to the inversion problem of three-dimensional gravity gradient data in the field of geophysical exploration, in particular to a three-dimensional inversion method of gravity gradient data based on partially smooth constraints. Background technique [0002] How to quickly and accurately determine the location and physical properties of geological bodies is one of the most important issues in prospecting geological exploration, and the physical property inversion technology of gravity or gravity gradient is the key technical link to solve this problem. In actual work, due to the limited effective data collected, the measurement accuracy is greatly disturbed by various factors, especially the equivalence of field sources makes the inversion of 3D physical properties have serious multi-solutions, so the inversion problem is generally an underdetermined problem . The most common way to solve this problem is to add prior information pr...

AI Technical Summary

Problems solved by technology

In actual work, due to the limited effective data collected, the measurement accuracy is greatly disturbed by various factors, especially the equivalence of field sources makes the 3D physical property inversion have serious multiple solutions, so the inversion problem is generally an underdetermined problem

Smooth inversion generally adopts global smoothing, and the obtained inversion results have blurred boundaries and expanded range, and the model resolution is low

Images