Abstract:SUMMARY Attention is paid to joint inversion of multiple geophysical data because of its advantages on weakening the non-uniqueness of inversion and further enhancing comprehensive interpretation. Due to the good correlation between rock velocity and density, seismic and gravity data have been widely used in joint inversion. However, the joint inversion of pre-stack seismic reflection and gravity data remains underdeveloped at the exploration scale. Without a quantitive relation between velocity and density, we develop a structure-based joint inversion using seismic reflection traveltimes, slopes and Bouguer gravity anomaly data simultaneously for building both velocity and density models. In our method, cubic B-spline interpolation is used to parametrize the common knots of velocity and density models. Incorporating seismic slopes into the joint inversion framework, we build a composite objective function which minimizes the weighted-sum of seismic/gravity data misfits, regularization and structural constraint terms. By subdividing the knot spacing, a multiscale strategy is alternative to increase the stability of inversion. First, we describe the methodology, followed by three synthetic examples to illustrate the feasibility and benefits of the method. Examination of the convergence curves via inversion suggests that the desired solution is more likely to be obtained with gentle convergence of each term, thus it can be used as an indicator for weight adjustment. Additionally, locations of scattering points and acoustic impedance can be obtained as by-products. Compared with the inversion of the respective data, the joint inversion exhibits the complementary characteristics of seismic and gravity data, improves the distribution and structural features of the resulting physical properties, especially in deep and complex tectonic situations.

What problem does this paper attempt to address?

This paper aims to solve several key problems in the joint inversion of seismic reflection data and gravity anomaly data. Specifically, the research focuses on: 1. **Reducing non - uniqueness**: The inversion of a single geophysical data usually has the problem of non - uniqueness, that is, different subsurface structures may produce similar observed data. By jointly inverting multiple geophysical data (such as seismic reflection data and gravity anomaly data), this non - uniqueness can be significantly reduced, and the reliability and accuracy of the inversion results can be improved. 2. **Enhancing comprehensive interpretation ability**: Seismic data are mainly used to infer the velocity information of strata, but are not sensitive to density changes; while Bouguer gravity anomaly data directly reflect density changes, but their inversion results are affected by factors such as non - uniqueness, resolution and depth sensitivity. Through joint inversion, the advantages of the two types of data can be fully utilized to improve the comprehensive interpretation ability of the distribution of subsurface physical properties and structural characteristics, especially in the case of deep and complex structures. 3. **Exploring the application of foreshock seismic reflection data**: Although the joint inversion of seismic and gravity data has been applied in some aspects, on the exploration scale, the joint inversion of foreshock seismic reflection data and gravity data is still in the preliminary stage. The identification and tracking of foreshock seismic reflection data are relatively difficult, especially in areas with complex geological conditions. Therefore, this research proposes a joint inversion method based on structural similarity, which uses seismic reflection travel time, slope and Bouguer gravity anomaly data to construct velocity and density models simultaneously to overcome this challenge. 4. **Developing a multi - scale strategy**: In order to increase the stability of the inversion, the research proposes a multi - scale strategy, which gradually refines the model by subdividing the node spacing. This method not only improves the stability of the inversion, but also can better capture the details of the subsurface structure. 5. **Optimizing the objective function**: The research constructs a composite objective function, including a data - fitting term, a regularization term and a structural constraint term. By adjusting the weights of each part, the influence of different data types and constraint conditions can be balanced, so as to obtain more reasonable inversion results. In summary, the main objective of this paper is to solve the technical problems in the joint inversion of foreshock seismic reflection data and gravity data by developing a new joint inversion method, and to improve the imaging quality of the subsurface structure and the description accuracy of physical parameters.

Joint inversion of seismic slopes, traveltimes and gravity anomaly data based on structural similarity

Alternating Joint Inversion of Controlled-Source Electromagnetic and Seismic Data Using the Joint Total Variation Constraint.

3-D joint inversion of seismic waveform and airborne gravity gradiometry data

Three-Dimensional Joint Inversion of MT and Gravity Data Based on Unstructured Tetrahedron Discretization

Two-Dimensional Joint Inversion of Seismic Velocity and Electrical Resistivity Using Seismic Travel Times and Full Channel Electrical Measurements Based on Alternating Cross-Gradient Structural Constraint

An efficient joint inversion strategy for 3D seismic travel time and DC resistivity data based on cross-gradient structure constraint

An Efficient Sequential Strategy for Realizing Cross-Gradient Joint Inversion: Method and Its Application to 2-D Cross Borehole Seismic Traveltime and DC Resistivity Tomography

3-D Cross-Gradient Joint Inversion of Seismic Refraction and Dc Resistivity Data

Advanced Methods of Joint Inversion of Multiphysics Data for Mineral Exploration

Joint Inversion of Gravity and Gravity Gradient Data Using Smoothed L0 Norm Regularization Algorithm with Sensitivity Matrix Compression

Joint Inversion of Gravity Gradient Tensor Data Based on L1 and L2 Norms

Binary structure constrained gravity inversion based on seismic first arrival travel time data

Joint data and model-driven simultaneous inversion of velocity and density

Time-Lapse Seismic Difference-and-Joint Prestack AVA Inversion

Joint Inversion of Gravity and Gravity Gradient Data Based on Cross-Gradient Function

A Joint Inversion Algorithm for the Establishment of High-Accuracy 3-D Marine Gravity Field

Structure-coupled joint inversion of geophysical and hydrological data

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

Cross-Gradient Joint Inversion of DC Resistivity and Gravity Gradient Data: A Multi-Disciplinary Approach for Geoscience, Heritage, and the Built Environment

Joint inversion in hydrogeophysics and near-surface geophysics

Structure constrained joint inversion of magnetotelluric data and receiver function