Seismic inversion method and system based on generalized total variation regularization
A seismic inversion and total variation technology, applied in the field of seismic exploration, can solve problems such as stratigraphic step effects
- Summary
- Abstract
- Description
- Claims
- Application Information
AI Technical Summary
Problems solved by technology
Method used
Image
Examples
Embodiment 1
[0142] This embodiment provides a seismic inversion method based on generalized total variation regularization, in which the second-order generalized total variation is used as an example for illustration, combined with figure 1 The flow process of the inventive method is described:
[0143] Step 1: Enter and create the data required for the inversion
[0144] Step 1.1: Input seismic data, logging data, horizon interpretation information;
[0145] Step 1.2: Extract the seismic wavelet W according to the data input in step 1.1, and create the initial model of the parameters to be inverted;
[0146] Step 2: Establish and solve the objective function
[0147] The objective function J is:
[0148]
[0149] where W is the wavelet matrix; L is the logarithm of the parameters to be inverted; S is the observed seismic record matrix; L' is the logarithm of the initial model of the parameters to be inverted; μ is the generalized full variational regularization coefficient; η is t...
Embodiment 2
[0204] This embodiment provides a seismic inversion system based on generalized total variation regularization based on the above embodiments, mainly including:
[0205] Objective function generation module: used to generate the objective function J of the generalized full variation regularization constraint with the logarithm of the parameter to be inverted;
[0206] Objective function transformation module: perform Fourier transformation on the objective function J to obtain the frequency domain expression of the objective function;
[0207] Inversion parameter solution module: Use the alternating direction multiplier method to solve the frequency domain expression of the objective function and iteratively update the parameters until the natural logarithm L of the parameter to be inverted is L 2 The inversion result is output when the relative change of the norm is smaller than the preset threshold tol.
Embodiment 3
[0209] On the basis of Embodiment 2, the system of the present application is further described in detail:
[0210] 1. The workflow of the objective function generation module is as follows:
[0211] Seismic data, well logging data, and horizon interpretation information need to be input. At the same time, seismic wavelet W must be extracted based on the data, and an initial model of the parameters to be inverted should be created. Then, the objective function J should be created based on the above parameters:
[0212]
[0213] where W is the wavelet matrix; L is the logarithm of the parameters to be inverted; S is the observed seismic record matrix; L' is the logarithm of the initial model of the parameters to be inverted; μ is the generalized full variational regularization coefficient; η is the initial model constraint regularization coefficient; v is the auxiliary variable; TGV(L) represents the generalized total variation of L, D is the difference matrix, α 0 、α 1 ar...
PUM
Abstract
Description
Claims
Application Information
- R&D Engineer
- R&D Manager
- IP Professional
- Industry Leading Data Capabilities
- Powerful AI technology
- Patent DNA Extraction
Browse by: Latest US Patents, China's latest patents, Technical Efficacy Thesaurus, Application Domain, Technology Topic, Popular Technical Reports.
© 2024 PatSnap. All rights reserved.Legal|Privacy policy|Modern Slavery Act Transparency Statement|Sitemap|About US| Contact US: help@patsnap.com