CN-121995451-A - Formation parameter determination method, apparatus, device, storage medium and program product

CN121995451ACN 121995451 ACN121995451 ACN 121995451ACN-121995451-A

Abstract

The application provides a formation parameter determination method, a device, equipment, a storage medium and a program product. Relates to the technical field of petroleum exploration seismic data acquisition and processing. The method comprises the steps of obtaining an initial stratum model, determining a target stratum inversion function and an iteration residual vector based on the initial stratum model and seismic exploration data, wherein the target stratum inversion function is determined by a Scholte wave dispersion curve, the iteration residual vector is determined by a Scholte wave deflection pilot dispersion curve, determining a stratum parameter updating formula based on the target stratum inversion function and the iteration residual vector, updating and iterating the initial stratum model based on the stratum parameter updating formula and a preset horizon layering rule to determine a target stratum model, and determining stratum parameters based on the target stratum model, wherein the stratum parameters comprise longitudinal wave velocity, transverse wave velocity and stratum medium density. The method effectively improves the accuracy and the reliability of the Scholte wave technology in seismic exploration.

Inventors

WANG YANJIAO
CHEN HAIFENG
WANG KEBIN
SUN PENGYUAN
LIU XIAOBO
ZHANG LINA

Assignees

中国石油天然气集团有限公司
中国石油集团东方地球物理勘探有限责任公司
中油油气勘探软件国家工程研究中心有限公司

Dates

Publication Date: 20260508
Application Date: 20241105

Claims (10)

1. A method of determining formation parameters, comprising: Acquiring an initial stratum model, wherein the initial stratum model is constructed based on known seismic exploration data, and comprises N stratum layers uniformly distributed along the depth direction, and N is a positive integer; Determining a target stratum inversion function and an iteration residual vector based on the initial stratum model and the seismic exploration data, wherein the target stratum inversion function is determined by a Scholte wave dispersion curve, and the iteration residual vector is determined by a Scholte wave deviation pilot dispersion curve; determining a stratum parameter updating formula based on the target stratum inversion function and the iteration residual vector; updating and iterating the initial stratum model based on the stratum parameter updating formula and a preset horizon layering rule to determine a target stratum model, wherein the target stratum model comprises M stratum horizons distributed along the depth direction, M is a positive integer, and M is smaller than N; formation parameters including compressional wave velocity, shear wave velocity, and formation media density are determined based on the target formation model.
2. The method of claim 1, wherein the determining a target formation inversion function and an iterative residual vector based on the initial formation model and the seismic survey data comprises: Determining an initial Scholte wave dispersion curve and an initial Scholte wave deflection frequency dispersion curve based on the initial stratum model; Determining a measured Scholte wave dispersion curve based on the seismic exploration data, and establishing a target stratum inversion function based on the measured Scholte wave dispersion curve and the initial Scholte wave dispersion curve; And determining an iteration residual vector based on the target stratum inversion function and the initial Scholte wave partial pilot scattering curve.
3. The method of claim 2, wherein the determining an initial Scholte waveguide dispersion curve based on the initial formation model comprises: Establishing a Helmholtz equation of any stratum layer based on the initial stratum model; Deriving the Helmholtz equation with respect to formation parameters to obtain a partial derivative equation; Solving a general solution of the partial derivative differential equation; determining displacement stress partial guide vectors of any stratum layer based on the general solution; establishing a recurrence relation based on the displacement stress deviation vector to obtain a recurrence relation between the displacement stress deviation vector of the first stratum layer and the displacement stress deviation vector of the Nth stratum layer; determining a Scholte wave bias pilot dispersion equation based on the recurrence relation; And solving a Scholte wave deviation pilot frequency dispersion equation to obtain the initial Scholte wave deviation pilot frequency dispersion curve.
4. The method of claim 1, wherein the target formation inversion function F (m) is expressed as: wherein f c is an initial Scholte wave dispersion curve, an To measure the Scholte wave dispersion curve.
5. The method of claim 1, wherein the iterative residual vector d k is represented as d k ＝-(G k TG k +α k I) -1 G k T f k , wherein d k is an iterative residual vector of a kth iteration, G k is a wave velocity bias matrix of the kth iteration, and the wave velocity bias matrix is obtained from the initial Scholte wave bias pilot dispersion curve, I is an identity matrix, a k is a positive real constant, and f k is a calculated value of a kth iteration of a target formation inversion function; The stratum parameter updating formula m k+1 is expressed as m k+1 ＝m k +λ k d k , wherein m k is a stratum parameter vector matrix of the kth iteration, m k+1 is a stratum parameter vector matrix of the (k+1) th iteration, and lambda k is the step length calculated by the kth iteration.
6. The method of claim 1, wherein updating the initial stratigraphic model based on the stratigraphic parameter updating formula and a preset horizon layering rule to determine a target stratigraphic model comprises: acquiring stratum parameters updated by the stratum parameter updating formula; Determining whether a formation parameter difference between two adjacent formation layers in the initial formation model is smaller than a first preset threshold value based on the updated formation parameters; And in response to the difference in formation parameters between two adjacent formation horizons being less than the first preset threshold, merging the two adjacent formation horizons to determine a target formation model.
7. A formation parameter determination apparatus, comprising: The acquisition module is used for acquiring an initial stratum model, the initial stratum model is constructed based on known seismic exploration data, the initial stratum model comprises N stratum layers which are uniformly distributed along the depth direction, and N is a positive integer; the determining module is used for determining a target stratum inversion function and an iteration residual vector based on the initial stratum model and the seismic exploration data, wherein the target stratum inversion function is determined by a Scholte wave dispersion curve, and the iteration residual vector is determined by a Scholte wave deviation pilot dispersion curve; the determining module is further configured to determine a formation parameter update formula based on the target formation inversion function and the iteration residual vector; The model updating module is used for updating and iterating the initial stratum model based on the stratum parameter updating formula and a preset horizon layering rule to determine a target stratum model, wherein the target stratum model comprises M stratum horizons distributed along the depth direction, M is a positive integer, and M is smaller than N; And the stratum parameter determining module is used for determining stratum parameters based on the target stratum model, wherein the stratum parameters comprise longitudinal wave velocity, transverse wave velocity and stratum medium density.
8. An electronic device comprising a processor and a memory communicatively coupled to the processor; The memory stores computer-executable instructions; the processor executes computer-executable instructions stored in the memory to implement the method of any one of claims 1-6.
9. A computer readable storage medium having stored therein computer executable instructions which when executed by a processor are adapted to carry out the method of any one of claims 1-6.
10. A computer program product comprising a computer program which, when executed by a processor, implements the method of any of claims 1-6.

Description

Formation parameter determination method, apparatus, device, storage medium and program product Technical Field The application relates to the technical field of petroleum exploration seismic data acquisition and processing, in particular to a stratum parameter determining method, a stratum parameter determining device, stratum parameter determining equipment, a storage medium and a program product. Background Currently, in the field of marine geological exploration, the exploration of the sedimentary structure of a submarine stratum is one of the main targets of marine exploration and geological exploration. Traditional geological exploration methods often rely on ground observation and drilling sampling, and the methods have the problems of high cost, long period, limited coverage range and the like. In recent years, with the development of marine geophysical technology, a stratum exploration method based on submarine Scholte waves gradually becomes a research hotspot. However, in the conventional geological research using Scholte waves, the number of stratum layers and the thickness of the stratum are required to be known, but the actual data of the number of stratum layers and the thickness of the stratum are not known, so that the actual data of the number of stratum layers and the thickness of the stratum are required to be set manually and empirically, and therefore, the geological research using Scholte waves is greatly limited, and the actual situation of geology is difficult to accurately reflect. Disclosure of Invention The application provides a stratum parameter determining method, a stratum parameter determining device, stratum parameter determining equipment, a stratum parameter determining storage medium and a stratum parameter determining program product, which are used for solving the problems that geological research by utilizing Scholte waves is greatly limited and real situations of geology are difficult to accurately reflect in the prior art. In a first aspect, the present application provides a method of determining a formation parameter, comprising: Acquiring an initial stratum model, wherein the initial stratum model is constructed based on known seismic exploration data, and comprises N stratum layers uniformly distributed along the depth direction, and N is a positive integer; Determining a target stratum inversion function and an iteration residual vector based on the initial stratum model and the seismic exploration data, wherein the target stratum inversion function is determined by a Scholte wave dispersion curve, and the iteration residual vector is determined by a Scholte wave deviation pilot dispersion curve; determining a stratum parameter updating formula based on the target stratum inversion function and the iteration residual vector; updating and iterating the initial stratum model based on the stratum parameter updating formula and a preset horizon layering rule to determine a target stratum model, wherein the target stratum model comprises M stratum horizons distributed along the depth direction, M is a positive integer, and M is smaller than N; formation parameters including compressional wave velocity, shear wave velocity, and formation media density are determined based on the target formation model. In one possible design, the determining a target formation inversion function and an iterative residual vector based on the initial formation model and the seismic survey data includes: Determining an initial Scholte wave dispersion curve and an initial Scholte wave deflection frequency dispersion curve based on the initial stratum model; Determining a measured Scholte wave dispersion curve based on the seismic exploration data, and establishing a target stratum inversion function based on the measured Scholte wave dispersion curve and the initial Scholte wave dispersion curve; And determining an iteration residual vector based on the target stratum inversion function and the initial Scholte wave partial pilot scattering curve. In one possible design, the determining an initial Scholte waveguide dispersion curve based on the initial formation model includes: Establishing a Helmholtz equation of any stratum layer based on the initial stratum model; Deriving the Helmholtz equation with respect to formation parameters to obtain a partial derivative equation; Solving a general solution of the partial derivative differential equation; determining displacement stress partial guide vectors of any stratum layer based on the general solution; establishing a recurrence relation based on the displacement stress deviation vector to obtain a recurrence relation between the displacement stress deviation vector of the first stratum layer and the displacement stress deviation vector of the Nth stratum layer; determining a Scholte wave bias pilot dispersion equation based on the recurrence relation; And solving a Scholte wave deviation pilot frequency dispersion equation to obtain the