Multi-parameter full-waveform inversion method and device
By employing a multi-parameter full-waveform inversion method and utilizing various parameterization modes of the velocity-density and modulus-density equations, the reliability issues of velocity, density, and modulus inversion in complex regions were resolved, thereby improving modeling accuracy and resolution.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-14
- Publication Date
- 2026-04-14
AI Technical Summary
Existing multi-parameter seismic inversion methods lack reliability in complex regions, especially in the full waveform inversion of velocity, density, and modulus, where effective methods are lacking.
A multi-parameter full waveform inversion method is adopted, which realizes the full waveform inversion of velocity, density and modulus through various parameterization modes of velocity-density and modulus-density equations, including the step-by-step establishment and inversion of velocity model, density model and modulus model.
It improves the modeling accuracy of full waveform inversion of velocity, density, and modulus, can clearly characterize the model, has higher vertical and horizontal resolution, and the results are closer to the real model.
Smart Images
Figure CN121857044A_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of seismic data processing and interpretation technology, specifically to a multi-parameter full waveform inversion method and apparatus. Background Technology
[0002] The development of multi-parameter seismic inversion has facilitated reservoir monitoring. With the advancement of seismic acquisition, processing, and interpretation technologies, multi-parameter seismic inversion has continuously progressed from theoretical methods to practical applications, and researchers have proposed various processing and interpretation methods and strategies accordingly. Traditional multi-parameter inversion methods have certain feasibility, but due to the inherent limitations of ray theory assumptions, the reliability of inversion results cannot be guaranteed in complex regions. With the development of full-waveform inversion, multi-parameter full-waveform inversion has also gradually developed, achieving certain results in both theoretical and practical applications, but it also faces greater challenges. Velocity, density, and modulus are key elastic parameters in seismic interpretation, and currently, there are no relevant full-waveform inversion methods for these three parameters.
[0003] Based on this technical background, the present invention studies a multi-parameter full waveform inversion method and apparatus. Summary of the Invention
[0004] To address the shortcomings of existing technologies, this invention provides a multi-parameter full waveform inversion method and apparatus. This method is based on the multi-parameter full waveform inversion technology of velocity, density, and modulus under multiple parameterization modes. It utilizes multiple parameterization modes to achieve full waveform inversion of the three parameters of velocity, density, and modulus, thereby improving the modeling accuracy of velocity, density, and modulus full waveform inversion.
[0005] To achieve the above objectives, a first aspect of the present invention provides a multi-parameter full waveform inversion method, comprising:
[0006] The velocity model is obtained by performing a full waveform inversion of velocity using the velocity-density equation;
[0007] Using the velocity model as input, the density model is obtained by performing full waveform inversion of density using the velocity-density equation;
[0008] Convert the velocity and density in the velocity and density models into moduli;
[0009] Using the density model and modulus as input, the modulus model is obtained by performing full waveform inversion of the modulus using the modulus-density equation;
[0010] Using the density model and modulus model as input, the final density model is obtained by performing full waveform inversion of density using the modulus-density equation.
[0011] A second aspect of the present invention provides a multi-parameter full waveform inversion apparatus, comprising:
[0012] The velocity model building module is used to obtain the velocity model by performing full waveform inversion of velocity using the velocity-density equation.
[0013] The density model establishment module is used to obtain the density model by performing full waveform inversion of density using the velocity-density equation, with the velocity model as input.
[0014] The modulus acquisition module is used to convert the velocity and density in the velocity model and density model into modulus.
[0015] The modulus model establishment module is used to obtain the modulus model by performing full waveform inversion of the modulus using the modulus-density equation, taking the density model and modulus as inputs.
[0016] The final density model building module is used to obtain the final density model by performing full waveform inversion of density using the density model and modulus model as inputs and the modulus-density equation.
[0017] A third aspect of the present invention provides an electronic device, the electronic device comprising:
[0018] Memory, which stores executable instructions;
[0019] A processor that executes the executable instructions in the memory to implement the multi-parameter full waveform inversion method described in the first aspect.
[0020] A fourth aspect of the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the multi-parameter full waveform inversion method described in the first aspect.
[0021] The beneficial effects of this invention include:
[0022] The multi-parameter full waveform inversion method proposed in this invention is based on the technology of multi-parameter full waveform inversion of velocity, density and modulus under multiple parameterization modes. It utilizes multiple parameterization modes to realize the full waveform inversion of the three parameters of velocity, density and modulus, thereby improving the modeling accuracy of the full waveform inversion of velocity, density and modulus.
[0023] Other features and advantages of the present invention will be described in detail in the following detailed description section. Attached Figure Description
[0024] The above and other objects, features and advantages of the present invention will become more apparent from the more detailed description of exemplary embodiments of the invention in conjunction with the accompanying drawings.
[0025] Figure 1 This is a flowchart illustrating the multi-parameter full waveform inversion method proposed in this invention.
[0026] Figure 2 This is a flowchart illustrating a specific implementation of the multi-parameter full waveform inversion method proposed in this invention.
[0027] Figure 3 This is a schematic diagram of the modulus-density parameter mode in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention.
[0028] Figure 4 This is a schematic diagram of the velocity-density parameter mode in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention.
[0029] Figure 5 This is a real velocity diagram in a specific implementation of the multi-parameter full waveform inversion method proposed in this invention.
[0030] Figure 6 This is a schematic diagram of the actual density in a specific implementation of the multi-parameter full waveform inversion method proposed in this invention.
[0031] Figure 7 This is a schematic diagram of the actual modulus in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention.
[0032] Figure 8 This is a schematic diagram of the initial velocity in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention.
[0033] Figure 9 This is a schematic diagram of the initial density in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention.
[0034] Figure 10 This is a schematic diagram of the initial modulus in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention.
[0035] Figure 11 This is a schematic diagram of the inversion speed in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention.
[0036] Figure 12 This is a schematic diagram of the inversion density in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention.
[0037] Figure 13 This is a schematic diagram of the inverted modulus in a specific embodiment of the multi-parameter full waveform inversion method proposed in this invention. Detailed Implementation
[0038] Preferred embodiments of the invention will now be described in more detail. While preferred embodiments of the invention are described below, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein.
[0039] This invention provides a multi-parameter full waveform inversion method, such as... Figure 1 As shown, it includes:
[0040] The velocity model is obtained by performing a full waveform inversion of velocity using the velocity-density equation;
[0041] Using the velocity model as input, the density model is obtained by performing full waveform inversion of density using the velocity-density equation;
[0042] Convert velocity and density in velocity and density models into modulus;
[0043] Using density model and modulus as input, modulus model is obtained by full waveform inversion of modulus using modulus-density equation;
[0044] Using the density model and modulus model as input, the final density model is obtained by performing full waveform inversion of density using the modulus-density equation.
[0045] In this invention, based on the technique of multi-parameter full waveform inversion of velocity, density, and modulus under multiple parameterization modes, the full waveform inversion of the three parameters of velocity, density, and modulus is achieved by using multiple parameterization modes, thereby improving the modeling accuracy of the full waveform inversion of velocity, density, and modulus.
[0046] According to the present invention, the expression for the velocity model is:
[0047]
[0048] Where ω is the angular frequency, ρ is the density, v is the velocity, G0(r,ω;x) is the Green's function solution of the background medium, P0(x,ω) is the background wave field, and Υ v (x) represents the velocity model.
[0049] According to the present invention, the expression for the density model is:
[0050]
[0051] Where θ is the angle between the incident wave and the reflected wave, Υ vρ-ρ (x) represents the density model.
[0052] According to the present invention, the formula used to convert velocity and density in the velocity model and density model into modulus is as follows:
[0053] κ=ρv 2 .
[0054] According to the present invention, the expression for the modulus model is:
[0055]
[0056] Among them, Υ κ (x) is the modulus model.
[0057] Preferably, the expression for the final density model is:
[0058]
[0059] Among them, Υ ρ (x) represents the final density model.
[0060] According to the present invention, the method further includes:
[0061] The final velocity, density, and modulus are output through the velocity model, final density model, and modulus model.
[0062] The present invention will be described in more detail below through embodiments.
[0063] Example 1:
[0064] like Figure 2 As shown, this embodiment proposes a multi-parameter full waveform inversion method, which realizes a fused parameterized mode multi-parameter full waveform inversion through velocity-density and modulus-density equations, thereby achieving velocity, modulus, and density inversion;
[0065] The basic principle of this method is that, under the modulus-density model, the contributions of large-angle and small-angle parameters to seismic data are roughly equal, i.e., the coupling degree is relatively strong (e.g., Figure 3 As shown: From top to bottom, the angle ranges from 0 degrees to 180 degrees. The top is 0 degrees, which is a small angle, and the bottom is 180 degrees, which is a large angle. This angle is the angle between the incident wave and the reflected wave. Figure 4 Similarly, in Figure 3 The solid line represents the modulus, and the dashed line represents the density. Figure 4 (The solid line represents density, and the dashed line represents velocity), making them difficult to distinguish; the relative coupling is weaker at medium angles; in the velocity-density mode, the coupling is stronger at small angles and weaker at large angles; therefore, the coupling characteristics of different parameter modes can be used to perform multi-parameter full waveform inversion;
[0066] The implementation steps of this method are as follows:
[0067] Step 1: Perform full velocity waveform inversion using the velocity-density equation:
[0068]
[0069] Where ω is the angular frequency, G0(r,ω;x) is the Green's function solution of the background medium, and P0(x,ω) is the background wave field;
[0070] Step 2: Using the model obtained in Step 1 as the velocity model, perform density full waveform inversion using the velocity-density equation:
[0071]
[0072] Here, θ represents the angle between the incident wave and the reflected wave;
[0073] Step 3: Combine the velocity and density obtained in Step 1 and Step 2 using the formula κ = ρv 2 Convert to modulus;
[0074] Step 4: Using the density model obtained in Step 2 and the modulus model obtained in Step 3 as inputs, perform full waveform inversion of the modulus using the modulus-density equation:
[0075]
[0076] Step 5: Using the density model obtained in Step 2 and the modulus model obtained in Step 4 as inputs, perform full waveform inversion of density using the modulus-density equation:
[0077]
[0078] This embodiment uses a two-dimensional model for testing; Figure 5 , Figure 6 , Figure 7 These are the actual velocity, density, and modulus, respectively. Figure 8 , Figure 9 , Figure 10 These are the initial velocity, density, and modulus, respectively. Figure 11 , Figure 12 , Figure 13 These are the velocity, density, and modulus inverted in the method of this invention, respectively. The inversion results show that this invention can clearly characterize the velocity, density, and modulus models. Compared with the initial model, it has higher vertical and horizontal resolution, is closer to the real model, and has obvious effects, stability, and practicality.
[0079] Example 2:
[0080] This embodiment provides a multi-parameter full waveform inversion method, such as... Figure 1 As shown, it includes:
[0081] The velocity model is obtained by performing a full waveform inversion of velocity using the velocity-density equation;
[0082] Using the velocity model as input, the density model is obtained by performing full waveform inversion of density using the velocity-density equation;
[0083] Convert velocity and density in velocity and density models into modulus;
[0084] Using density model and modulus as input, modulus model is obtained by full waveform inversion of modulus using modulus-density equation;
[0085] Using the density model and modulus model as inputs, the final density model is obtained by performing full waveform inversion of density using the modulus-density equation.
[0086] In this embodiment, the expression for the velocity model is:
[0087]
[0088] Where ω is the angular frequency, ρ is the density, v is the velocity, G0(r,ω;x) is the Green's function solution of the background medium, P0(x,ω) is the background wave field, and Υ v (x) represents the velocity model;
[0089] In this embodiment, the expression for the density model is:
[0090]
[0091] Where θ is the angle between the incident wave and the reflected wave, Υ vρ-ρ (x) represents the density model;
[0092] In this embodiment, the formula used to convert velocity and density in the velocity model and density model into modulus is as follows:
[0093] κ=ρv 2 ;
[0094] In this embodiment, the expression for the modulus model is:
[0095]
[0096] Among them, Υ κ (x) represents the modulus model;
[0097] In this embodiment, the expression for the final density model is:
[0098]
[0099] Among them, Υ ρ (x) represents the final density model;
[0100] In this embodiment, the method further includes:
[0101] The final velocity, density, and modulus are output through the velocity model, final density model, and modulus model.
[0102] Example 3:
[0103] This embodiment provides a multi-parameter full waveform inversion device, including:
[0104] The velocity model building module is used to obtain the velocity model by performing full waveform inversion of velocity using the velocity-density equation.
[0105] The density model building module is used to obtain the density model by performing full waveform inversion of density using the velocity-density equation, with the velocity model as input.
[0106] The modulus acquisition module is used to convert velocity and density in the velocity and density models into modulus.
[0107] The modulus model building module is used to obtain the modulus model by taking the density model and modulus as input and using the modulus-density equation to perform full waveform inversion of the modulus.
[0108] The final density model building module is used to obtain the final density model by taking the density model and the modulus model as input and using the modulus-density equation to perform full waveform inversion of density.
[0109] In this embodiment, the expression for the velocity model is:
[0110]
[0111] Where ω is the angular frequency, ρ is the density, v is the velocity, G0(r,ω;x) is the Green's function solution of the background medium, P0(x,ω) is the background wave field, and Υ v (x) represents the velocity model;
[0112] In this embodiment, the expression for the density model is:
[0113]
[0114] Where θ is the angle between the incident wave and the reflected wave, Υ vρ-ρ (x) represents the density model;
[0115] In this embodiment, the formula used to convert velocity and density in the velocity model and density model into modulus is as follows:
[0116] κ=ρv 2 ;
[0117] In this embodiment, the expression for the modulus model is:
[0118]
[0119] Among them, Υ κ (x) represents the modulus model;
[0120] In this embodiment, the expression for the final density model is:
[0121]
[0122] Where, γ ρ (x) represents the final density model;
[0123] In this embodiment, the method further includes:
[0124] The final velocity, density, and modulus are output through the velocity model, final density model, and modulus model.
[0125] Example 4:
[0126] This invention provides an electronic device including a memory and a processor.
[0127] Memory, which stores executable instructions;
[0128] The processor executes executable instructions in memory to implement a multi-parameter full waveform inversion method.
[0129] This memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.
[0130] The processor may be a central processing unit (CPU) or other form of processing unit with data processing and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the invention, the processor is used to execute computer-readable instructions stored in the memory.
[0131] Those skilled in the art will understand that, in order to solve the technical problem of how to achieve a good user experience, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included within the protection scope of this invention.
[0132] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.
[0133] Example 5:
[0134] This invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements a multi-parameter full waveform inversion method.
[0135] A computer-readable storage medium according to embodiments of the present invention stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of the present invention are performed.
[0136] The aforementioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).
[0137] The multi-parameter full waveform inversion method proposed in the embodiments of the present invention is based on the multi-parameter full waveform inversion technology of velocity, density and modulus under multiple parameterization modes. It utilizes multiple parameterization modes to realize the full waveform inversion of the three parameters of velocity, density and modulus, thereby improving the modeling accuracy of the full waveform inversion of velocity, density and modulus.
[0138] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.
Claims
1. A multi-parameter full waveform inversion method, characterized in that, include: The velocity model is obtained by performing a full waveform inversion of velocity using the velocity-density equation; Using the velocity model as input, the density model is obtained by performing full waveform inversion of density using the velocity-density equation; Convert the velocity and density in the velocity and density models into moduli; Using the density model and modulus as input, the modulus model is obtained by performing full waveform inversion of the modulus using the modulus-density equation; Using the density model and modulus model as input, the final density model is obtained by performing full waveform inversion of density using the modulus-density equation.
2. The method according to claim 1, characterized in that, The expression for the velocity model is: Where ω is the angular frequency, ρ is the density, v is the velocity, G0(r,ω;x) is the Green's function solution of the background medium, P0(x,ω) is the background wave field, and γ v (x) represents the velocity model.
3. The method according to claim 2, characterized in that, The expression for the density model is: Where θ is the angle between the incident wave and the reflected wave, γ vρ-ρ (x) represents the density model.
4. The method according to claim 3, characterized in that, The formula used to convert velocity and density in the velocity and density models into modulus is: k=ρv 2 。 5. The method according to claim 4, characterized in that, The expression for the modulus model is: Where, γ κ (x) is the modulus model.
6. The method according to claim 5, characterized in that, The expression for the final density model is: Where, γ ρ (x) represents the final density model.
7. The method according to claim 1, characterized in that, Also includes: The final velocity, density, and modulus are output through the velocity model, final density model, and modulus model.
8. A multi-parameter full waveform inversion device, characterized in that, include: The velocity model building module is used to obtain the velocity model by performing full waveform inversion of velocity using the velocity-density equation. The density model establishment module is used to obtain the density model by performing full waveform inversion of density using the velocity-density equation, with the velocity model as input. The modulus acquisition module is used to convert the velocity and density in the velocity model and density model into modulus. The modulus model establishment module is used to obtain the modulus model by performing full waveform inversion of the modulus using the modulus-density equation, taking the density model and modulus as inputs. The final density model building module is used to obtain the final density model by performing full waveform inversion of density using the density model and modulus model as inputs and the modulus-density equation.
9. An electronic device, characterized in that, The electronic device includes: Memory, which stores executable instructions; A processor that executes the executable instructions in the memory to implement the multi-parameter full waveform inversion method according to any one of claims 1-7.
10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the multi-parameter full waveform inversion method according to any one of claims 1-7.