CN-116338782-B - Crack zone plane feature detection method, device, equipment and medium

CN116338782BCN 116338782 BCN116338782 BCN 116338782BCN-116338782-B

Abstract

The method comprises the steps of obtaining three-dimensional seismic data corresponding to a plurality of space positions in a fracture zone, determining inclination angles and azimuth angles of geological textures in target space positions and coordinates of the space positions corresponding to the geological textures according to the three-dimensional seismic data, determining analysis time windows corresponding to the geological textures, calculating curved surface coefficients of the geological textures according to the analysis time windows, guiding the curved surface coefficients into a fracture enhancement curvature function to obtain curvature values of the geological textures, calculating projections of the geological textures corresponding to the space positions and the ground in parallel according to the curvature values of the geological textures, carrying out projection recombination on all the space positions to obtain fracture zone images, and enhancing description effects on the geological texture characteristics through the method.

Inventors

PENG DA
XU XIANG

Assignees

中国石油天然气股份有限公司

Dates

Publication Date: 20260512
Application Date: 20211223

Claims (6)

1. The crack band plane characteristic detection method is characterized by comprising the following steps of: acquiring three-dimensional seismic data corresponding to a plurality of spatial positions in a fracture zone; determining the inclination angle and the azimuth angle of the geological texture at the target space position according to the three-dimensional seismic data; determining an analysis time window of the geological texture according to the inclination angle and the azimuth angle of the geological texture and the coordinates of the corresponding target space position of the geological texture; calculating the surface coefficient of the geological texture according to the analysis time window; the curvature value of the geological texture is obtained by bringing the curvature coefficient of the geological texture into a crack enhancement curvature function, wherein the crack enhancement curvature function is a first curvature function representing the anticline characteristic of the geological texture And a second curvature function characterizing syncline features of the geologic texture The expression form is as follows: ; Calculating the projection of the geological texture corresponding to the target space position in the direction parallel to the ground according to the curvature value of the geological texture; By the method, projections of geological textures corresponding to all spatial positions in the fracture zone are obtained, projection recombination is carried out, and a fracture zone image is obtained; the determining the analysis time window of the geological texture according to the inclination angle and the azimuth angle of the geological texture and the coordinates of the target space position corresponding to the geological texture further comprises: acquiring the inclination angle and the azimuth angle of the geological texture at the target space position; According to the formula Constructing the analysis time window of the geological texture; Wherein, to the right of the equation And Inclination and azimuth angles of the geological texture at the target spatial position, respectively, to the right of the equation 、 And Coordinates of the target space position, respectively, to the left of the equation Expressed as an analysis time window function, wherein the analysis time window function is used for constructing the analysis time window; The calculating the surface coefficient of the geological texture according to the analysis time window further comprises: According to the formula 、、、 And Calculating a plurality of surface coefficients a, b, c, d and e of the geological texture; wherein a is an inclination angle direction differential coefficient, b is an azimuth angle direction differential coefficient, c is a differential mean coefficient, d is an inclination angle of the geological texture, and e is an azimuth angle of the geological texture; P is the analysis time window function Q is the analysis time window function Is a direction angle of (a).
2. The method of claim 1, wherein determining the dip and azimuth of the geologic texture at the target spatial location from the three-dimensional seismic data comprises: three gradient vector bodies which are in one-to-one correspondence with three directions of geological textures at the target space position are obtained from the three-dimensional seismic data; performing vector-parallel operation on the three gradient vector bodies to obtain a gradient structure tensor matrix of the geological texture; carrying out spectrum decomposition on the gradient structure tensor matrix to obtain three eigenvectors corresponding to the geological texture and three eigenvalues corresponding to the three eigenvectors one by one; selecting a feature vector corresponding to the largest feature value in the three feature values; and determining the inclination angle and the azimuth angle of the geological texture according to the selected feature vector.
3. The method of claim 2, wherein before acquiring three gradient vector volumes corresponding one-to-one to three directions of a geological texture at a target spatial location from the three-dimensional seismic data, comprising: And filtering the three-dimensional seismic data through a three-dimensional Gaussian smoothing filter.
4. A crack strip plane feature detection device, comprising: The data acquisition unit is used for acquiring three-dimensional seismic data corresponding to a plurality of spatial positions in the fracture zone; The parameter acquisition unit is used for determining the inclination angle and the azimuth angle of the geological texture at the target space position according to the three-dimensional seismic data; The time window creation unit is used for determining an analysis time window of the geological texture according to the inclination angle and the azimuth angle of the geological texture and the coordinates of the target space position corresponding to the geological texture; The curved surface coefficient calculating unit is used for calculating the curved surface coefficient of the geological texture according to the analysis time window; A curvature value calculation unit, configured to bring the curvature coefficient of the geological texture into a fracture enhancement curvature function to obtain a curvature value of the geological texture, where the fracture enhancement curvature function is a first curvature function characterizing a anticline feature of the geological texture And a second curvature function characterizing syncline features of the geologic texture The expression form is as follows: ; the projection calculation unit is used for calculating the projection of the geological texture corresponding to the target space position in the direction parallel to the ground according to the curvature value of the geological texture; the image reorganization unit is used for acquiring projections of geological textures corresponding to all the spatial positions in the fracture zone through the units, and carrying out projection reorganization to obtain fracture zone images; the time window creation unit is specifically configured to: acquiring the inclination angle and the azimuth angle of the geological texture at the target space position; According to the formula Constructing the analysis time window of the geological texture; Wherein, to the right of the equation And Inclination and azimuth angles of the geological texture at the target spatial position, respectively, to the right of the equation 、 And Coordinates of the target space position, respectively, to the left of the equation Expressed as an analysis time window function, wherein the analysis time window function is used for constructing the analysis time window; The surface coefficient calculating unit is specifically configured to: According to the formula 、、、 And Calculating a plurality of surface coefficients a, b, c, d and e of the geological texture; wherein a is an inclination angle direction differential coefficient, b is an azimuth angle direction differential coefficient, c is a differential mean coefficient, d is an inclination angle of the geological texture, and e is an azimuth angle of the geological texture; P is the analysis time window function Q is the analysis time window function Is a direction angle of (a).
5. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the crack zone planar feature detection method according to any one of claims 1-3 when executing the computer program.
6. A computer readable storage medium, characterized in that the computer readable storage medium stores a computer program, which when executed by a processor implements the crack strip plane feature detection method according to any one of claims 1-3.

Description

Crack zone plane feature detection method, device, equipment and medium Technical Field The invention relates to the technical field of petroleum exploration, and can be used in the field of geological curvature of oil and gas reservoirs, in particular to a method, a device, equipment and a medium for detecting planar characteristics of a fracture zone. Background More than half of the global oil and gas production is statistically distributed in fractured reservoirs, however, because of the unique low pore, heterogeneity and anisotropy of carbonate reservoirs, the difficulty of oil and gas exploration is greatly increased, and the fractures are the channels of the main oil and gas reservoirs and migration of such reservoirs, so that accurately predicting the fractures in such reservoirs is equivalent to finding oil and gas directly in such reservoirs. At present, a scheme for detecting cracks in a crack zone through an image is required to design the size of an analysis time window through manual experience, and the geological texture in three-dimensional seismic data is vertically and horizontally scanned according to the analysis time window, and due to the limitation of the size and scanning direction of the analysis time window, only the curved surface coefficient of the geological texture at the central position in the analysis time window can be selected, and the characteristics of all the geological textures contained in the analysis time window are represented, so that the data scanning is incomplete, and the real geological appearance is not consistent with the characteristics fed back by the image. Disclosure of Invention In view of the foregoing problems in the prior art, an object of the present disclosure is to provide a method, an apparatus, a device, and a medium for detecting a planar feature of a fractured zone, so as to solve the problem in the prior art that data scanning is incomplete, and a real geological appearance does not coincide with a feature fed back by an image. In order to solve the technical problems, the specific technical scheme is as follows: in one aspect, provided herein is a fracture zone planar feature detection method comprising: acquiring three-dimensional seismic data corresponding to a plurality of spatial positions in a fracture zone; determining the inclination angle and the azimuth angle of the geological texture at the target space position according to the three-dimensional seismic data; determining an analysis time window of the geological texture according to the inclination angle and the azimuth angle of the geological texture and the coordinates of the corresponding target space position of the geological texture; calculating the surface coefficient of the geological texture according to the analysis time window; Bringing the surface coefficient of the geological texture into a crack enhancement curvature function to obtain a curvature value of the geological texture; Calculating the projection of the geological texture corresponding to the target space position in the direction parallel to the ground according to the curvature value of the geological texture; By the method, projections of geological textures corresponding to all spatial positions in the fracture zone are obtained, projection recombination is carried out, and a fracture zone image is obtained. As one embodiment herein, the fracture enhancement curvature function is derived from a first curvature function and a second curvature function, further comprising: Performing mean square error calculation on the first curvature function and the second curvature function to obtain the crack enhancement curvature function; Wherein the first curvature function characterizes anticline features of the geological texture; The second curvature function characterizes syncline features of the geological texture. As one embodiment, the determining the analysis time window of the geological texture according to the inclination angle and the azimuth angle of the geological texture and the coordinates of the spatial position of the target corresponding to the geological texture further comprises: acquiring the inclination angle and the azimuth angle of the geological texture at the target space position; Constructing the analysis time window of the geological texture according to the formula g (t, p, q) =g (t-px-qy); wherein p and q on the right of the equation are the inclination angle and azimuth angle of the geological texture at the target space position, respectively, x and y on the right of the equation are the coordinates of the target space position, respectively, and g (t, p, q) is expressed as an analysis time window function, wherein the analysis time window function is used for constructing the analysis time window. As one embodiment herein, the calculating the surface coefficient of the geological texture according to the analysis time window further includes: According to the formula D=p (x, y, t