CN117826256A - Method and device for determining sequence boundary of high-frequency seismic sequence stratum - Google Patents

Abstract

The invention discloses a method and a device for determining an interval boundary of a high-frequency earthquake interval stratum, wherein the method comprises the following steps: according to the three-dimensional seismic data volume, a relative geologic time model is generated, an initial medium-low frequency first sequence boundary data volume based on the earthquake is extracted, a medium-low frequency first sequence boundary curve based on the earthquake is extracted, baseline removal processing is carried out, a second sequence boundary data volume based on the earthquake is obtained through first waveform difference inversion, a third sequence boundary curve based on the earthquake is extracted and used as a basic curve, and is fused with a gamma curve, a porosity curve, a carbonate particle curve and a carbonate texture curve which reflect high-frequency sequence boundaries, a carbonate suture line frequency curve, an asphalt content curve and a dolomite content curve which reflect the sequence boundaries, a Gao Pindi sequence boundary data volume based on the earthquake is obtained after second waveform difference inversion, and the high-frequency earthquake sequence boundary is identified.

Description

Method and device for determining sequence boundaries of high-frequency seismic sequence stratigraphy

Technical Field

The present invention relates to the field of carbonate stratum oil and gas exploration and development, and the technical field of seismic division of high-frequency sea level change cyclic sequence bodies, and in particular to a method and device for determining the sequence boundary of high-frequency seismic sequence strata.

Background technique

This section is intended to provide a background or context to the embodiments of the invention recited in the claims. No description herein is admitted to be prior art by inclusion in this section.

A stratigraphic unit consists of strata deposited in all genetically related depositional environments during a complete base-level cycle. Semi-cycle boundaries of a genetic sequence occur at transitions from base-level rise to fall or from base-level fall to rise. In different paleogeographic settings, these transitions are manifested as stratigraphic discontinuities or conformable strata recording increases or decreases in accommodation space, respectively, usually forming sequence boundaries.

In sequence stratigraphy, sequences are usually divided into 1 to 6 levels, among which 1 to 3 levels correspond to tectonic sedimentary cycles, belong to low-frequency sequences (cycles), and correspond to giant sequences, super sequences, and sequences respectively; 4 to 6 levels of sequences correspond to climatic sedimentary cycles, belong to high-frequency sequences (cycles), also known as Milankovitch cycles, and correspond to parasequence groups, parasequences, and rhythmic layers respectively. In sequence stratigraphy research, sequence division is the basis, and sequence interface identification is the key to sequence division. Generally, sequence interfaces include unconformity surfaces, transgression overlay surfaces, flooding unconformity surfaces, paleokarst action surfaces, volcanic event action surfaces, and lithology conversion surfaces. Different levels of sequences correspond to different levels of sequence interfaces. For example, regional unconformity surfaces and tectonic conversion surfaces often correspond to low-frequency sequences, while lithology conversion surfaces and sedimentary unconformity surfaces correspond to high-frequency sequences.

The commonly used methods for determining sequence boundaries in geology are: lithology and lithofacies changes, core observation, INPEFA curve, and wavelet variation spectrum. However, in areas where carbonate rocks are the main research object, the lithofacies include grainstone limestone, grain-containing mud-powder limestone, mud-crystal limestone (containing particles), dolomite and dolomitic limestone, gypsum-containing mud-crystal limestone, etc., and multiple sets of thin-layered shallows are developed. The bioclastic beach is relatively developed. Affected by multiple cycles, multiple stages and multiple sets of beach bodies are developed. For carbonate rock formations with strong heterogeneity, the existing sequence boundary determination methods can usually only reflect the boundaries of low-frequency sequence formations, and cannot identify the requirements of higher-frequency sequence formation frameworks. The current technical solutions cannot effectively determine the sequence boundaries of high-frequency seismic sequence formations.

Summary of the invention

In a first aspect, an embodiment of the present invention provides a method for determining a sequence boundary of a high-frequency seismic sequence stratum, which can establish a more accurate high-frequency sequence boundary volume, so that the sequence boundary of the identified high-frequency seismic sequence stratum is more accurate. The method comprises:

Calculate the geological model grid based on the three-dimensional seismic data volume of the target layer segment in the target area;

Generate a relative geological age model based on the geological model grid;

The relative geological age model is processed and analyzed to extract the initial seismic-based low- and medium-frequency first sequence boundary data volume;

Extracting the seismic-based medium-low-frequency first sequence boundary curve of the target layer segment at the target well point from the initial seismic-based medium-low-frequency first sequence boundary data volume;

Performing baseline removal processing on the first sequence boundary curve of medium and low frequency based on seismic data to obtain the second sequence boundary curve of medium and low frequency based on seismic data after baseline removal processing;

According to the low- and medium-frequency second sequence boundary curve after baseline removal based on seismic, the second sequence boundary data volume based on seismic is obtained through the first waveform difference inversion;

Extracting the seismic-based third sequence boundary curve of the target layer segment at the target well point from the seismic-based second sequence boundary data volume;

The third sequence boundary curve based on seismic is used as the basic curve, and is fused with the gamma curve reflecting the high-frequency sequence boundary at the target well point, the porosity curve reflecting the high-frequency sequence boundary, the carbonate rock grain curve reflecting the high-frequency sequence boundary, the carbonate rock texture curve reflecting the high-frequency sequence boundary, the carbonate rock suture line frequency curve reflecting the sequence boundary, the asphalt content curve of the carbonate rock reflecting the sequence boundary and the dolomite content curve of the carbonate rock reflecting the sequence boundary to obtain the high-frequency tenth sequence boundary curve based on seismic;

According to the high-frequency tenth-order sequence boundary curve based on seismic data, a second waveform difference inversion is performed on the target layer section in the target area to obtain the high-frequency third-order sequence boundary data volume based on seismic data.

The sequence boundaries of high-frequency seismic sequence stratigraphy are identified based on the high-frequency third sequence boundary data volume based on seismic.

In a second aspect, an embodiment of the present invention further provides a device for determining a sequence boundary of a high-frequency seismic sequence stratum, which can establish a more accurate high-frequency sequence boundary volume, so that the sequence boundary of the identified high-frequency seismic sequence stratum is more accurate, and the device comprises:

A geological model grid calculation module is used to calculate the geological model grid according to the three-dimensional seismic data volume of the target layer section in the target area;

A relative geological age model generation module is used to generate a relative geological age model based on a geological model grid;

The sequence thickness boundary data volume extraction module is used to process and analyze the relative geological age model and extract the initial seismic-based medium- and low-frequency first sequence boundary data volume;

Seismic-based sequence boundary curves are used to extract seismic-based medium-low frequency first sequence boundary curves of target layer segments at target well points from initial seismic-based medium-low frequency first sequence boundary data volumes; perform baseline removal on seismic-based medium-low frequency first sequence boundary curves to obtain seismic-based medium-low frequency second sequence boundary curves after baseline removal; obtain seismic-based second sequence boundary data volumes through first waveform difference inversion based on seismic-based medium-low frequency second sequence boundary curves after baseline removal; extract seismic-based third sequence boundary curves of target layer segments at target well points from seismic-based second sequence boundary data volumes;

A fusion processing module is used to fuse the third-order sequence boundary curve based on seismic as a basic curve with the gamma curve reflecting the high-frequency sequence boundary at the target well point, the porosity curve reflecting the high-frequency sequence boundary, the carbonate rock particle curve reflecting the high-frequency sequence boundary, the carbonate rock texture curve reflecting the high-frequency sequence boundary, the carbonate rock suture line frequency curve reflecting the sequence boundary, the asphalt content curve of the carbonate rock reflecting the sequence boundary and the dolomite content curve of the carbonate rock reflecting the sequence boundary, so as to obtain the high-frequency tenth-order sequence boundary curve based on seismic;

The high-frequency seismic sequence boundary identification module is used to perform a second waveform difference inversion on the target layer segment in the target area based on the high-frequency tenth sequence boundary curve based on seismic, and obtain the high-frequency third sequence boundary data body based on seismic; based on the high-frequency third sequence boundary data body based on seismic, the sequence boundary of the high-frequency seismic sequence strata is identified.

In a third aspect, an embodiment of the present invention further provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for determining the stratigraphic boundaries of high-frequency seismic sequence strata when executing the computer program.

In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned method for determining the sequence boundaries of high-frequency seismic sequence strata is implemented.

In a fifth aspect, an embodiment of the present invention further provides a computer program product, which includes a computer program, and when the computer program is executed by a processor, it implements the above-mentioned method for determining the sequence boundaries of high-frequency seismic sequence strata.

In an embodiment of the present invention, a geological model grid is calculated based on a three-dimensional seismic data volume of a target layer segment in a target area; a relative geological age model is generated based on the geological model grid; the relative geological age model is processed and analyzed to extract an initial seismic-based medium-low frequency first sequence boundary data volume; a seismic-based medium-low frequency first sequence boundary curve of the target layer segment at a target well point is extracted from the initial seismic-based medium-low frequency first sequence boundary data volume; a seismic-based medium-low frequency first sequence boundary curve is subjected to baseline removal processing to obtain a seismic-based medium-low frequency second sequence boundary curve after baseline removal; a seismic-based second sequence boundary data volume is obtained through first waveform difference inversion based on the seismic-based medium-low frequency second sequence boundary curve after baseline removal; and a seismic-based second sequence boundary data volume is extracted from the seismic-based second sequence boundary data volume. a third sequence boundary curve based on seismic; taking the third sequence boundary curve based on seismic as the basic curve, the gamma curve reflecting the high-frequency sequence boundary at the target well point, the porosity curve reflecting the high-frequency sequence boundary, the carbonate rock particle curve reflecting the high-frequency sequence boundary, the carbonate rock texture curve reflecting the high-frequency sequence boundary, the carbonate rock suture line frequency curve reflecting the sequence boundary, the asphalt content curve of the carbonate rock reflecting the sequence boundary and the dolomite content curve of the carbonate rock reflecting the sequence boundary are integrated to obtain the tenth high-frequency sequence boundary curve based on seismic; according to the tenth high-frequency sequence boundary curve based on seismic, the target layer section in the target area is subjected to a second waveform difference inversion to obtain the third high-frequency sequence boundary data body based on seismic; according to the third high-frequency sequence boundary data body based on seismic, the sequence boundary of the high-frequency seismic sequence strata is identified. Compared with the existing sequence interface determination method, after extracting the seismic-based medium- and low-frequency first sequence boundary curve of the target layer segment at the target well point, the embodiment of the present invention undergoes multiple processing, including baseline processing, inversion processing, extraction processing and fusion with multiple curves reflecting high-frequency sequence boundaries and curves reflecting sequence boundaries, to obtain an accurate high-frequency sequence boundary curve, and finally performs inversion based on waveform differences to obtain a high-frequency sequence boundary data body. According to the high-frequency sequence boundary data body, the sequence boundary of the high-frequency seismic sequence strata is identified. In the high-frequency sequence boundary body, the larger the value, the closer it is to the sequence boundary. The above process fuses the gamma curve and porosity curve at the target well point, and the inverted sequence boundary features are more obvious, thereby improving the accuracy and precision of identifying the sequence interface body using seismic-geological-well logging fusion information, so that the sequence boundary of the identified high-frequency seismic sequence strata is more accurate.

BRIEF DESCRIPTION OF THE DRAWINGS

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the prior art descriptions. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. In the drawings:

FIG1 is a flow chart of a method for determining sequence boundaries of high-frequency seismic sequence strata in an embodiment of the present invention;

FIG2 is a flow chart of calculating a geological model grid in an embodiment of the present invention;

FIG3 is a flow chart of generating a relative geological age model in an embodiment of the present invention;

FIG4 is a flow chart of multi-curve fusion in an embodiment of the present invention;

FIG5 is an example diagram of a relative geological age model calculated based on a three-dimensional seismic data volume in the region according to an embodiment of the present invention;

6 is an example diagram of a well-connected profile obtained by analyzing the relative geological age model and extracting the initial seismic-based medium- and low-frequency first sequence boundary data volume in an embodiment of the present invention;

FIG7 is an example diagram of a sequence boundary curve based on earthquakes in an embodiment of the present invention;

FIG8 is an example diagram of a well-connected cross section of a second sequence boundary data volume based on seismic data obtained by the first waveform difference inversion in an embodiment of the present invention;

FIG9 is a comparison diagram of gamma curves before and after processing according to an embodiment of the present invention;

FIG10 is a comparison diagram of porosity curves before and after processing in an embodiment of the present invention;

FIG11 is an example of a process of fusing and obtaining a high-frequency sequence boundary curve in an embodiment of the present invention;

FIG12 is an example of a well-connected section of a high-frequency sequence boundary data volume in an embodiment of the present invention;

FIG13 is a structural block diagram of a device for determining sequence boundaries of high-frequency seismic sequence strata according to an embodiment of the present invention;

FIG. 14 is a schematic diagram of a computer device in an embodiment of the present invention.

Detailed ways

To make the purpose, technical solution and advantages of the embodiments of the present invention more clear, the embodiments of the present invention are further described in detail below in conjunction with the accompanying drawings. Here, the exemplary embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

FIG1 is a flow chart of a method for determining sequence boundaries of high-frequency seismic sequence strata according to an embodiment of the present invention, comprising:

Step 101, calculating a geological model grid according to a three-dimensional seismic data volume of a target layer segment in a target area;

Step 102, generating a relative geological age model according to the geological model grid;

Step 103, processing and analyzing the relative geological age model to extract the initial seismic-based medium-low frequency first sequence boundary data volume;

Step 104, extracting the seismic-based medium-low-frequency first sequence boundary curve of the target layer segment at the target well point from the initial seismic-based medium-low-frequency first sequence boundary data volume;

Step 105, performing baseline removal processing on the medium-low frequency first sequence boundary curve based on seismic data to obtain the medium-low frequency second sequence boundary curve after baseline removal based on seismic data;

Step 106, obtaining a second sequence boundary data volume based on seismic data by performing a first waveform difference inversion based on the low- and medium-frequency second sequence boundary curve after the seismic baseline removal process;

Step 107, extracting a third-order sequence boundary curve based on seismic data of a target layer segment at a target well point from the second-order sequence boundary data volume based on seismic data;

Step 108, taking the third sequence boundary curve based on seismic as the base curve, the gamma curve reflecting the high-frequency sequence boundary at the target well point, the porosity curve reflecting the high-frequency sequence boundary, the carbonate rock particle curve reflecting the high-frequency sequence boundary, the carbonate rock texture curve reflecting the high-frequency sequence boundary, the carbonate rock suture line frequency curve reflecting the sequence boundary, the asphalt content curve of the carbonate rock reflecting the sequence boundary and the dolomite content curve of the carbonate rock reflecting the sequence boundary are merged to obtain the tenth high-frequency sequence boundary curve based on seismic;

Step 109, performing a second waveform difference inversion on the target layer section in the target area according to the high-frequency tenth-order sequence boundary curve based on seismic, and obtaining a high-frequency third-order sequence boundary data volume based on seismic;

Step 110, identifying the sequence boundaries of the high-frequency seismic sequence strata based on the high-frequency third sequence boundary data volume based on the seismic.

In an embodiment of the present invention, a geological model grid is calculated based on a three-dimensional seismic data volume of a target layer segment in a target area; a relative geological age model is generated based on the geological model grid; the relative geological age model is processed and analyzed to extract an initial seismic-based medium-low frequency first sequence boundary data volume; from the initial seismic-based medium-low frequency first sequence boundary data volume, a seismic-based medium-low frequency first sequence boundary curve of the target layer segment at a target well point is extracted; the seismic-based medium-low frequency first sequence boundary curve is subjected to baseline removal processing to obtain a seismic-based medium-low frequency second sequence boundary curve after baseline removal; based on the seismic-based medium-low frequency second sequence boundary curve after baseline removal, a seismic-based second sequence boundary data volume is obtained through first waveform difference inversion; from the seismic-based second sequence boundary data volume, the target layer segment at the target well point is extracted a third sequence boundary curve based on seismic; taking the third sequence boundary curve based on seismic as the basic curve, the gamma curve reflecting the high-frequency sequence boundary at the target well point, the porosity curve reflecting the high-frequency sequence boundary, the carbonate rock particle curve reflecting the high-frequency sequence boundary, the carbonate rock texture curve reflecting the high-frequency sequence boundary, the carbonate rock suture line frequency curve reflecting the sequence boundary, the asphalt content curve of the carbonate rock reflecting the sequence boundary and the dolomite content curve of the carbonate rock reflecting the sequence boundary are integrated to obtain the tenth high-frequency sequence boundary curve based on seismic; according to the tenth high-frequency sequence boundary curve based on seismic, the target layer section in the target area is subjected to a second waveform difference inversion to obtain the third high-frequency sequence boundary data body based on seismic; according to the third high-frequency sequence boundary data body based on seismic, the sequence boundary of the high-frequency seismic sequence strata is identified. Compared with the existing sequence interface determination method, after extracting the seismic-based medium- and low-frequency first sequence boundary curve of the target layer segment at the target well point, the embodiment of the present invention undergoes multiple processing, including baseline processing, inversion processing, extraction processing and fusion with multiple curves reflecting high-frequency sequence boundaries and curves reflecting sequence boundaries, to obtain an accurate high-frequency sequence boundary curve, and finally performs inversion based on waveform differences to obtain a high-frequency sequence boundary data body. According to the high-frequency sequence boundary data body, the sequence boundary of the high-frequency seismic sequence strata is identified. In the high-frequency sequence boundary body, the larger the value, the closer it is to the sequence boundary. The above process fuses the gamma curve and porosity curve at the target well point, and the inverted sequence boundary features are more obvious, thereby improving the accuracy and precision of identifying the sequence interface body using seismic-geological-well logging fusion information, so that the sequence boundary of the identified high-frequency seismic sequence strata is more accurate.

In step 101, a geological model grid is calculated based on a three-dimensional seismic data volume of a target layer segment in a target area;

This step combines the actual well logging and three-dimensional seismic data in the seismic data for calculation.

Referring to FIG. 2 , the geological model grid is calculated based on the 3D seismic data volume of the target layer in the target area, including:

Step 201, determine an initial geological model grid based on the three-dimensional seismic data volume of the target layer segment in the target area; two methods for determining the initial geological model grid are provided here.

In one embodiment, determining an initial geological model grid based on a three-dimensional seismic data volume of a target layer segment in a target area includes:

When there is seismic layer data in the three-dimensional seismic data volume, the seismic layer data in the three-dimensional seismic data volume is automatically tracked, and a geological model grid is established with the tracked seismic layer data as a constraint; wherein the seismic layer data is already present in the three-dimensional seismic data volume;

When there is no seismic layer data in the three-dimensional seismic data body, an initial geological model grid is calculated for at least one seed point in the three-dimensional seismic data body of the target layer segment of the target area according to waveform similarity and relative distance. When calculating specifically, a certain algorithm may be used, such as an algorithm based on boundary control-local mapping, etc., which is not limited here.

Step 202: Using the initial geological model grid as the current geological model grid, the following steps are repeatedly performed until the current geological model grid and the three-dimensional seismic data volume are consistent with the preset conditions:

Step 2021, in the current geological model grid, interactively correct the association relationship between seismic layers; each correction here will affect the links between nodes in the model grid;

Step 2022, analyzing the consistency between the modified geological model grid and the three-dimensional seismic data volume; in specific implementation, the analysis can be performed through preview;

Step 2023, when the matching condition does not meet the preset conditions, optimize the parameters of the current geological model grid, and use the optimized geological model grid as the current geological model grid.

Through the above iterative cycles, the optimal geological model grid can be obtained.

In step 102, a relative geological age model is generated based on the geological model grid;

Referring to FIG3 , in one embodiment, generating a relative geological age model based on a geological model grid includes:

Step 301, connecting and interpolating the facets of the geological model grid to obtain a processed geological model grid;

Step 302, assigning a relative geological age to each pixel in the processed geological model grid to generate an initial relative geological age model;

Step 303, extract multiple stratigraphic stacks from the initial relative geological age model to form a relative geological age model represented by multiple stratigraphic stacks. The stratigraphic stacks generally include tens of thousands of stratigraphic stacks, and the relative geological age model represented by tens of thousands of stratigraphic stacks is more accurate.

Step 103, the relative geological age model is processed and analyzed to extract the initial low- to medium-frequency first sequence boundary data body based on earthquakes; the sequence thickness boundary data body is an attribute body data that can reflect the stratigraphic sequence boundary, which can reflect the changes in the vertical and horizontal sequence boundaries and be combined with earthquakes for sequence interpretation. The larger the value, the greater the difference in the same geological age, and the more likely it is a sequence boundary. This is an interpretation method for establishing a high-frequency sequence stratigraphic framework based on the global thinking concept of earthquake-geology.

Step 104, extracting the seismic-based medium-low-frequency first sequence boundary curve of the target layer segment at the target well point from the initial seismic-based medium-low-frequency first sequence boundary data volume;

In one embodiment, extracting a seismic-based medium-low-frequency first sequence boundary curve of a target layer segment at a target well point from an initial seismic-based medium-low-frequency first sequence boundary data volume includes:

The well-seismic calibration and the first waveform difference inversion are performed on the medium-low frequency second sequence boundary curve after the seismic baseline removal to obtain the seismic second sequence boundary data volume.

The medium-low frequency first sequence boundary curve based on seismic is subjected to baseline removal processing to obtain the medium-low frequency second sequence boundary curve based on seismic after baseline removal processing. The medium-low frequency second sequence boundary curve obtained after baseline removal processing can highlight the recognition ability of the sequence thickness boundary curve obtained from the seismic data body for the sequence boundary.

In step 108, the third-order sequence boundary curve based on seismic is used as a base curve, and is fused with a gamma curve reflecting high-frequency sequence boundaries at a target well point, a porosity curve reflecting high-frequency sequence boundaries, a carbonate rock grain curve reflecting high-frequency sequence boundaries, a carbonate rock texture curve reflecting high-frequency sequence boundaries, a carbonate rock suture frequency curve reflecting sequence boundaries, an asphalt content curve of carbonate rocks reflecting sequence boundaries, and a dolomite content curve of carbonate rocks reflecting sequence boundaries to obtain a high-frequency tenth-order sequence boundary curve based on seismic;

Referring to FIG. 4 , the specific fusion steps include:

Step 401, using the third-order sequence boundary curve based on seismic as a basic curve, discretizing it into a plurality of column data points, normalizing the gamma curve reflecting the high-frequency sequence boundary, discretizing it into a plurality of column data points, superimposing it on the plurality of column data points corresponding to the third-order sequence boundary curve based on seismic at the corresponding depth point, and fusing it into a fourth-order high-frequency sequence boundary curve based on seismic;

In one embodiment, the method further comprises:

Perform reverse processing on the gamma curve at the target well point to obtain a reverse gamma curve;

Detrending the reverse gamma curve to obtain a detrended gamma curve;

The detrended gamma curve is then subjected to baseline removal to obtain a gamma curve reflecting the high-frequency sequence boundary.

Step 402, using the fourth-order sequence boundary curve based on seismic as a base curve, discretizing it into a plurality of columns of data points, normalizing the porosity curve reflecting the high-frequency sequence boundary, discretizing it into a plurality of columns of data points, and superimposing it on the plurality of columns of data points corresponding to the fourth-order sequence boundary curve based on seismic at the corresponding depth point, and merging them into a fifth-order high-frequency sequence boundary curve based on seismic;

In one embodiment, the method further comprises:

The porosity curve at the target well point is de-baselined to obtain a porosity curve reflecting the high-frequency sequence boundary.

Step 403, using the fifth sequence boundary curve based on seismic data as a base curve, discretizing it into a plurality of columns of data points, normalizing and discretizing the carbonate rock grain curve reflecting the high-frequency sequence boundary into a plurality of columns of data points, superimposing the curve onto the plurality of columns of data points corresponding to the fifth sequence boundary curve based on seismic data at the corresponding depth point, and merging them into a sixth high-frequency sequence boundary curve based on seismic data; wherein the carbonate rock grain curve reflecting the high-frequency sequence boundary is a curve formed by quantizing carbonate rock grains into a plurality of columns of data points from grain limestone to argillaceous limestone;

For example, multiple column data points quantified from granular limestone to muddy limestone can be -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6. The larger the negative value, the smaller the particles and the higher the mud content. The larger the positive value, the larger the carbonate rock particles and the higher the carbonate rock particle content.

Among them, the larger the carbonate rock particles, the stronger the energy of seawater activity. Therefore, the carbonate rock particle curve integrated with normalization can better reflect the boundary information of carbonate rock sequence.

Step 404, using the sixth sequence boundary curve based on seismic as a base curve, discretizing it into a plurality of column data points, normalizing and discretizing the carbonate rock texture curve reflecting the high-frequency sequence boundary into a plurality of column data points, superimposing it onto the plurality of column data points corresponding to the sixth sequence boundary curve based on seismic at the corresponding depth point, and merging them into the seventh high-frequency sequence boundary curve based on seismic; wherein the carbonate rock texture curve reflecting the high-frequency sequence boundary is a curve formed by quantizing the rock texture structure of carbonate rock into a plurality of column data points;

For example, the rock texture structure of carbonate rocks is quantified as integers from 1 to 15, with larger values indicating coarser texture;

Among them, the carbonate rock texture curve represents the structure of the carbonate sedimentary sequence and the characteristics of the sedimentary cycle. The fusion of the carbonate rock texture curve can more finely reflect the high-frequency sequence boundaries of the carbonate rock.

Step 405, taking the seventh-order sequence boundary curve based on seismic as the base curve, discretizing it into a plurality of column data points, normalizing and discretizing the suture frequency curve of carbonate rocks reflecting the high-frequency sequence boundary into a plurality of column data points, and superimposing it onto the plurality of column data points corresponding to the seventh-order sequence boundary curve based on seismic at the corresponding depth point, and merging them into the eighth-order high-frequency sequence boundary curve based on seismic; wherein the suture frequency curve of carbonate rocks reflecting the high-frequency sequence boundary is a curve formed by quantizing the suture frequency of carbonate rocks into a plurality of column data points;

For example, the stylolite frequency of carbonate rocks is quantified as an integer from 1 to 9, where a larger number indicates more stylolites.

Among them, the suture frequency curve of carbonate rocks reflecting high-frequency sequence boundaries reflects the changes in structure and sequence boundaries caused by diagenesis. The larger the number of sutures, the stronger the progradation characteristics caused by sea retreat, and the more obvious the sequence boundaries caused by coarse-grained deposition.

Step 406, using the eighth order sequence boundary curve based on seismic as a base curve, discretizing it into a plurality of column data points, normalizing the asphalt content curve of carbonate rock reflecting the sequence boundary, discretizing it into a plurality of column data points, and superimposing it on the plurality of column data points corresponding to the eighth order sequence boundary curve based on seismic at the corresponding depth point, and merging them into the ninth order high-frequency sequence boundary curve based on seismic; wherein the asphalt content curve of carbonate rock reflecting the sequence boundary is a curve formed by quantifying the asphalt content of carbonate rock into a plurality of column data points;

For example, the bitumen content of carbonate rocks is quantified from 0.1 to 0.9, with the larger value indicating more bitumen content;

Among them, the asphalt content curve of carbonate rocks reflects the changes in temperature, pressure, and composition during the carbonate rock accumulation process, which causes changes in sequence boundaries. Therefore, integrating the asphalt content curve of carbonate rocks that reflects sequence boundaries can help finely identify sequence boundaries.

Step 407, taking the ninth order sequence boundary curve based on seismic as a base curve, discretizing it into a plurality of columns of data points, normalizing the dolomite content curve of carbonate rock reflecting the sequence boundary, discretizing it into a plurality of columns of data points, and superimposing it on the plurality of columns of data points corresponding to the ninth order sequence boundary curve based on seismic at the corresponding depth point, and merging them into the tenth order high-frequency sequence boundary curve based on seismic; wherein the dolomite content curve of carbonate rock reflecting the sequence boundary is a curve formed by quantifying the dolomite to limestone content of carbonate rock into a plurality of columns of data points.

For example, the dolomite to limestone content of carbonate rocks is quantified as 0.1 to 0.9, and the larger the value, the more dolomite content. Among them, the higher the dolomite content indicates that the carbonate rock depositional environment is in an exposed environment, the shallower the sea water, and the carbonate rock sequence is located in the regressive environment of the upper half cycle. The dolomite content curve of the carbonate rock that reflects the sequence boundary more finely reflects the high-frequency sequence boundary.

The tenth high-frequency sequence boundary curve finally obtained based on seismic data can reflect the information of high-frequency sequence boundaries.

In step 109, according to the high-frequency tenth-order sequence boundary curve based on seismic data, a second waveform difference inversion is performed on the target layer section in the target area to obtain a high-frequency third-order sequence boundary data body based on seismic data. At this time, the high-frequency sequence boundary data body is clearer and more accurate.

In step 110, the sequence boundaries of the high-frequency seismic sequence stratigraphy are identified based on the high-frequency third sequence boundary data body based on seismic. Finally, the sequence boundaries of the high-frequency seismic sequence stratigraphy with higher accuracy can be automatically identified using the global thinking concept based on the high-frequency sequence interface data body.

The method proposed in the embodiment of the present invention can be used for carbonate sequence interfaces, which solves the technical problem that only low-cycle interfaces can be obtained in the existing process of obtaining carbonate sequence interfaces, and the sequence interface resolution is low. In areas where carbonate oil and gas reservoirs are mainly developed, the genetic mechanism of carbonate reservoir heterogeneity is a basic geological problem that restricts efficient development, especially for complex carbonate curtain sequence cycles, sedimentary patterns and corresponding multi-stage diagenetic evolution. Research is weak, and it is necessary to establish a high-frequency sedimentary cycle sequence stratigraphic framework to assist in the study of carbonate diagenetic evolution history and the coupled reservoir control mechanism of multiple geological factors. Therefore, the method of the present invention has broad application prospects.

A specific embodiment is given below to illustrate the specific application of the method proposed in the present invention.

Taking a certain area where carbonate oil and gas reservoirs are mainly developed as an example, FIG. 5 is an example diagram of a relative geological age model calculated based on a three-dimensional seismic data body in the embodiment of the present invention, wherein (a) in FIG. 5 is a three-dimensional seismic data body, (b) in FIG. 5 is a geological model grid, and (c) in FIG. 5 is a relative geological age model. FIG. 6 is an example diagram of a well-connected profile of an initial seismic-based medium- and low-frequency first sequence boundary data body analyzed and extracted in an embodiment of the present invention for the relative geological age model, showing the instantaneous change of the relative geological age of each seismic sample point, highlighting the convergence and divergence areas of the geological layer. It is sensitive to unconformity surfaces, formation termination (underlying, onlapping), erosion, compaction, and formation thickness. The initial seismic-based medium- and low-frequency first sequence boundary data body is equal to the relative isochronous geological time divided by the seismic two-way reflection time interval. If the denominator seismic two-way reflection time interval is the same, the larger the numerator relative isochronous geological time, the larger the sequence thickness boundary data body, indicating that the formation is thicker.

FIG7 is an example diagram of a sequence boundary curve based on seismic in an embodiment of the present invention, wherein (a) in FIG7 is a first sequence boundary curve of medium and low frequency based on seismic, and (b) in FIG7 is a curve after baseline removal processing of the initial sequence boundary curve. FIG8 is an example diagram of a well-connected profile of a second sequence boundary data volume based on seismic obtained by first waveform difference inversion in an embodiment of the present invention, and the third sequence boundary curve based on seismic extracted finally can highlight the recognition ability of the sequence thickness boundary curve obtained from seismic data for sequence boundaries.

FIG9 is a comparison diagram of gamma curves before and after processing in an embodiment of the present invention, wherein FIG9(a) is an initial gamma curve, and FIG9(b) is a gamma curve reflecting the high-frequency sequence boundary.

FIG10 is a comparison diagram of the porosity curves before and after processing in an embodiment of the present invention, wherein FIG10 (a) is the initial porosity curve, and FIG10 (b) is the porosity curve reflecting the high-frequency sequence boundary.

FIG11 is an example of the process of obtaining the high-frequency tenth sequence boundary curve based on seismic by fusion in an embodiment of the present invention, wherein (a) in FIG11 is the third sequence boundary curve based on seismic, (b) in FIG11 is the sequence boundary curve after processing, (c) in FIG11 is the initial gamma curve, (d) in FIG11 is the gamma curve reflecting the high-frequency sequence boundary, (e) in FIG11 is the initial porosity curve, (f) in FIG11 is the porosity curve reflecting the high-frequency sequence boundary, and (g) in FIG11 is the fifth high-frequency sequence boundary curve based on seismic obtained by fusion. Other curves can be further fused later to finally obtain the tenth high-frequency sequence boundary curve based on seismic. FIG12 is an example of a well-connected section of the high-frequency third sequence boundary data body based on seismic in an embodiment of the present invention. At this time, the high-frequency third sequence boundary data body based on seismic is clearer and more accurate. Finally, the sequence boundary of the high-frequency seismic sequence strata with higher accuracy can be automatically identified by using the global thinking concept based on the high-frequency sequence interface data body.

The embodiment of the present invention further provides a device for determining and generating sequence boundaries of high-frequency seismic sequence strata, the principle of which is similar to the method for determining and generating sequence boundaries of high-frequency seismic sequence strata, and will not be described in detail here.

FIG13 is a schematic diagram of a device for determining and generating sequence boundaries of high-frequency seismic sequence strata according to an embodiment of the present invention, comprising:

The geological model grid calculation module 1301 is used to calculate the geological model grid according to the three-dimensional seismic data volume of the target layer segment in the target area;

A relative geological age model generation module 1302 is used to generate a relative geological age model according to the geological model grid;

The sequence thickness boundary data volume extraction module 1303 is used to process and analyze the relative geological age model and extract the initial low-frequency first sequence boundary data volume based on earthquakes;

Seismic-based sequence boundary curve 1304 is used to extract a seismic-based medium-low frequency first sequence boundary curve of a target layer segment at a target well point from an initial seismic-based medium-low frequency first sequence boundary data volume; perform baseline removal processing on the seismic-based medium-low frequency first sequence boundary curve to obtain a seismic-based medium-low frequency second sequence boundary curve after baseline removal; obtain a seismic-based second sequence boundary data volume through a first waveform difference inversion based on the seismic-based medium-low frequency second sequence boundary curve after baseline removal; extract a seismic-based third sequence boundary curve of a target layer segment at a target well point from the seismic-based second sequence boundary data volume;

A fusion processing module 1305 is used to fuse the third-order sequence boundary curve based on seismic as a basic curve with a gamma curve reflecting the high-frequency sequence boundary at a target well point, a porosity curve reflecting the high-frequency sequence boundary, a carbonate rock particle curve reflecting the high-frequency sequence boundary, a carbonate rock texture curve reflecting the high-frequency sequence boundary, a carbonate rock suture frequency curve reflecting the sequence boundary, an asphalt content curve of carbonate rocks reflecting the sequence boundary, and a dolomite content curve of carbonate rocks reflecting the sequence boundary, so as to obtain a high-frequency tenth-order sequence boundary curve based on seismic;

The high-frequency seismic sequence boundary identification module 1306 is used to perform a second waveform difference inversion on the target layer segment in the target area based on the high-frequency tenth-order sequence boundary curve based on seismic, and obtain the high-frequency third-order sequence boundary data body based on seismic; and identify the sequence boundaries of the high-frequency seismic sequence strata based on the high-frequency third-order sequence boundary data body based on seismic.

In one embodiment, the geological model grid calculation module is specifically used for:

Determine the initial geological model grid based on the 3D seismic data volume of the target layer segment in the target area;

The initial geological model grid is used as the current geological model grid, and the following steps are repeated until the current geological model grid and the 3D seismic data volume meet the preset conditions:

Interactively modify the correlation between seismic layers in the current geological model grid;

Analyze the consistency between the modified geological model grid and the 3D seismic data volume;

When the matching situation does not meet the preset conditions, the parameters of the current geological model grid are optimized, and the optimized geological model grid is used as the current geological model grid.

In one embodiment, the geological model grid calculation module is specifically used for:

When the three-dimensional seismic data volume contains seismic layer data, the seismic layer data in the three-dimensional seismic data volume is automatically tracked, and a geological model grid is established with the tracked seismic layer data as a constraint;

When there is no seismic layer data in the three-dimensional seismic data volume, an initial geological model grid is calculated for at least one seed point in the three-dimensional seismic data volume of the target layer segment in the target area according to waveform similarity and relative distance.

In one embodiment, the relative geological age model generation module is specifically used to:

Connecting and interpolating the surface elements of the geological model grid to obtain a processed geological model grid;

Assigning relative geological age to each pixel in the processed geological model grid to generate an initial relative geological age model;

A plurality of stratigraphic stacks are extracted from the initial relative geological age model to form a relative geological age model represented by the plurality of stratigraphic stacks.

In one embodiment, the seismic-based sequence boundary curve is specifically used to:

The well-seismic calibration and the first waveform difference inversion are performed on the medium-low frequency second sequence boundary curve after the seismic baseline removal to obtain the seismic second sequence boundary data volume.

In one embodiment, the fusion processing module is specifically used for:

Perform reverse processing on the gamma curve at the target well point to obtain a reverse gamma curve;

Detrending the reverse gamma curve to obtain a detrended gamma curve;

The detrended gamma curve is then subjected to baseline removal to obtain a gamma curve reflecting the high-frequency sequence boundary.

In one embodiment, the fusion processing module is specifically used for:

The porosity curve at the target well point is de-baselined to obtain a porosity curve reflecting the high-frequency sequence boundary.

In one embodiment, the fusion processing module is specifically used for:

The third sequence boundary curve based on earthquake is used as the basic curve, which is discretized into multiple columns of data points. The gamma curve reflecting the high-frequency sequence boundary is normalized and discretized into multiple columns of data points, which are then superimposed on the multiple columns of data points corresponding to the third sequence boundary curve based on earthquake at the corresponding depth point, and merged into the fourth high-frequency sequence boundary curve based on earthquake;

The fourth-order sequence boundary curve based on seismic is used as the basic curve, which is discretized into multiple columns of data points. The porosity curve reflecting the high-frequency sequence boundary is normalized and discretized into multiple columns of data points, which are then superimposed on the multiple columns of data points corresponding to the fourth-order sequence boundary curve based on seismic at the corresponding depth point, and merged into the fifth-order high-frequency sequence boundary curve based on seismic.

The fifth sequence boundary curve based on seismic data is used as a basic curve and discretized into multiple columns of data points. The carbonate rock grain curve reflecting the high-frequency sequence boundary is normalized and discretized into multiple columns of data points, which are then superimposed on the multiple columns of data points corresponding to the fifth sequence boundary curve based on seismic data at the corresponding depth point to merge into the sixth high-frequency sequence boundary curve based on seismic data; wherein the carbonate rock grain curve reflecting the high-frequency sequence boundary is a curve formed by quantizing carbonate rock grains into multiple columns of data points from grain limestone to argillaceous limestone;

The sixth sequence boundary curve based on seismic data is used as a basic curve and discretized into multiple columns of data points. The carbonate rock texture curve reflecting the high-frequency sequence boundary is normalized and discretized into multiple columns of data points, which are then superimposed on the multiple columns of data points corresponding to the sixth sequence boundary curve based on seismic data at the corresponding depth point to merge into the seventh high-frequency sequence boundary curve based on seismic data; wherein the carbonate rock texture curve reflecting the high-frequency sequence boundary is a curve formed by quantizing the rock texture structure of carbonate rock into multiple columns of data points;

The seventh-order sequence boundary curve based on seismic is used as the basic curve, which is discretized into a plurality of columns of data points. The suture frequency curve of carbonate rocks reflecting the high-frequency sequence boundary is normalized and discretized into a plurality of columns of data points, which are then superimposed on the plurality of columns of data points corresponding to the seventh-order sequence boundary curve based on seismic at the corresponding depth point, and merged into the eighth-order high-frequency sequence boundary curve based on seismic; wherein the suture frequency curve of carbonate rocks reflecting the high-frequency sequence boundary is a curve formed by quantizing the suture frequency of carbonate rocks into a plurality of columns of data points;

The eighth sequence boundary curve based on seismic is used as the basic curve, which is discretized into multiple columns of data points. The asphalt content curve of carbonate rock reflecting the sequence boundary is normalized and discretized into multiple columns of data points, which are then superimposed on the multiple columns of data points corresponding to the eighth sequence boundary curve based on seismic at the corresponding depth point, and merged into the ninth high-frequency sequence boundary curve based on seismic; wherein the asphalt content curve of carbonate rock reflecting the sequence boundary is a curve formed by quantifying the asphalt content of carbonate rock into multiple columns of data points;

The ninth sequence boundary curve based on seismic is taken as the basic curve and discretized into multiple columns of data points. The dolomite content curve of carbonate rock reflecting the sequence boundary is normalized and discretized into multiple columns of data points, which are then superimposed on the multiple columns of data points corresponding to the ninth sequence boundary curve based on seismic at the corresponding depth point to merge into the tenth high-frequency sequence boundary curve based on seismic; wherein, the dolomite content curve of carbonate rock reflecting the sequence boundary is a curve formed by quantifying the dolomite to limestone content of carbonate rock into multiple columns of data points.

In summary, in the method and device proposed in the embodiment of the present invention, a geological model grid is calculated based on the three-dimensional seismic data volume of the target layer segment in the target area; a relative geological age model is generated based on the geological model grid; the relative geological age model is processed and analyzed to extract the initial seismic-based medium-low frequency first sequence boundary data volume; from the initial seismic-based medium-low frequency first sequence boundary data volume, the seismic-based medium-low frequency first sequence boundary curve of the target layer segment at the target well point is extracted; the seismic-based medium-low frequency first sequence boundary curve is subjected to baseline removal processing to obtain the seismic-based medium-low frequency second sequence boundary curve after baseline removal; based on the seismic-based medium-low frequency second sequence boundary curve after baseline removal, the seismic-based second sequence boundary data volume is obtained through the first waveform difference inversion; the target well point is extracted from the seismic-based second sequence boundary data volume. A third seismic-based sequence boundary curve of the target layer segment at the well point; taking the third seismic-based sequence boundary curve as the basic curve, the gamma curve reflecting the high-frequency sequence boundary, the porosity curve reflecting the high-frequency sequence boundary, the carbonate rock particle curve reflecting the high-frequency sequence boundary, the carbonate rock texture curve reflecting the high-frequency sequence boundary, the carbonate rock suture line frequency curve reflecting the sequence boundary, the asphalt content curve of the carbonate rock reflecting the sequence boundary and the dolomite content curve of the carbonate rock reflecting the sequence boundary at the target well point are integrated to obtain a high-frequency tenth seismic-based sequence boundary curve; according to the high-frequency tenth seismic-based sequence boundary curve, a second waveform difference inversion is performed on the target layer segment in the target area to obtain a high-frequency third sequence boundary data body based on seismic; according to the high-frequency third sequence boundary data body based on seismic, the sequence boundary of the high-frequency seismic sequence stratum is identified. Compared with the existing sequence interface determination method, after extracting the seismic-based medium- and low-frequency first sequence boundary curve of the target layer segment at the target well point, the embodiment of the present invention undergoes multiple processing, including baseline processing, inversion processing, extraction processing and fusion with multiple curves reflecting high-frequency sequence boundaries and curves reflecting sequence boundaries, to obtain an accurate high-frequency sequence boundary curve, and finally performs inversion based on waveform differences to obtain a high-frequency sequence boundary data body. According to the high-frequency sequence boundary data body, the sequence boundary of the high-frequency seismic sequence strata is identified. In the high-frequency sequence boundary body, the larger the value, the closer it is to the sequence boundary. The above process fuses the gamma curve and porosity curve at the target well point, and the inverted sequence boundary features are more obvious, thereby improving the accuracy and precision of identifying the sequence interface body using seismic-geological-well logging fusion information, so that the sequence boundary of the identified high-frequency seismic sequence strata is more accurate.

An embodiment of the present invention further provides a computer device. FIG14 is a schematic diagram of a computer device in an embodiment of the present invention. The computer device 1400 includes a memory 1410, a processor 1420, and a computer program 1430 stored in the memory 1410 and executable on the processor 1420. When the processor 1420 executes the computer program 1430, the above-mentioned method for determining the sequence boundaries of high-frequency seismic sequence strata is implemented.

An embodiment of the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned method for determining the sequence boundary of high-frequency seismic sequence strata is implemented.

An embodiment of the present invention further provides a computer program product, which includes a computer program. When the computer program is executed by a processor, the method for determining the sequence boundaries of high-frequency seismic sequence strata is implemented.

Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

The present invention is described with reference to the flowchart and/or block diagram of the method, device (system), and computer program product according to the embodiment of the present invention. It should be understood that each process and/or box in the flowchart and/or block diagram, as well as the combination of the process and/or box in the flowchart and/or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device that implements the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

These computer program instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims (19)

1. A method for determining a sequence boundary of a high frequency seismic sequence formation, comprising:

calculating a geological model grid according to the three-dimensional seismic data volume of the target interval of the target area;

generating a relative geologic time model according to the geologic model grid;

processing and analyzing the relative geologic time model, and extracting an initial medium-low frequency first sequence boundary data body based on earthquake;

Extracting a seismic-based medium-low frequency first sequence boundary curve of a target interval at a target well point from an initial seismic-based medium-low frequency first sequence boundary data volume;

performing baseline removal processing on the middle-low frequency first sequence boundary curve based on the earthquake to obtain a middle-low frequency second sequence boundary curve after the baseline removal processing based on the earthquake;

according to the medium-low frequency second sequence boundary curve after the baseline removal processing based on the earthquake, obtaining a second sequence boundary data body based on the earthquake through first waveform difference inversion;

extracting a third sub-sequence boundary curve based on the earthquake of the target interval at the target well point from the second sequence boundary data body based on the earthquake;

taking a third sequence boundary curve based on an earthquake as a basic curve, and fusing the third sequence boundary curve with a gamma curve reflecting a high-frequency sequence boundary, a porosity curve reflecting the high-frequency sequence boundary, a carbonate grain curve reflecting the high-frequency sequence boundary, a carbonate suture line frequency curve reflecting the sequence boundary, an asphalt content curve reflecting the carbonate of the sequence boundary and a dolomite content curve reflecting the carbonate of the sequence boundary at a target well point to obtain a tenth sequence boundary curve based on the earthquake;

According to a high-frequency tenth-order boundary curve based on the earthquake, carrying out second waveform difference inversion on a target interval of the target area to obtain a Gao Pindi three-order boundary data body based on the earthquake;

based on the Gao Pindi three-layer sequence boundary data body based on the earthquake, the sequence boundary of the high-frequency earthquake sequence stratum is identified.

2. The method of claim 1, wherein computing a geologic model grid from the three-dimensional seismic data volume of the interval of interest of the target area comprises:

determining an initial geological model grid according to the three-dimensional seismic data volume of the target interval of the target area;

taking the initial geological model grid as the current geological model grid, and repeatedly executing the following steps until the matching condition of the current geological model grid and the three-dimensional seismic data volume meets the preset condition:

interactively correcting the association relation among the seismic horizons in the current geological model grid;

analyzing the coincidence condition of the corrected geological model grid and the three-dimensional seismic data volume;

when the fit condition does not meet the preset condition, optimizing parameters of the current geological model grid, and taking the optimized geological model grid as the current geological model grid.

3. The method of claim 2, wherein determining an initial geologic model grid from the three-dimensional seismic data volume of the destination interval of the target area comprises:

when the three-dimensional seismic data volume has seismic horizon data, automatically tracking the seismic horizon data in the three-dimensional seismic data volume, and establishing a geological model grid by taking the tracked seismic horizon data as a constraint;

and when the three-dimensional seismic data volume does not contain the seismic horizon data, calculating an initial geological model grid according to the waveform similarity and the relative distance for at least one seed point in the three-dimensional seismic data volume of the target interval of the target area.

4. The method of claim 1, wherein generating a relative geologic time model from a geologic model grid comprises:

connecting and interpolating the surface element sheets of the geological model grid to obtain a processed geological model grid;

distributing relative geologic time to each pixel in the processed geologic model grid to generate an initial relative geologic time model;

a plurality of horizon stacks are extracted from the initial relative geologic time model to form a relative geologic time model represented by the plurality of horizon stacks.

5. The method of claim 1, wherein obtaining the seismic-based second sequence boundary data volume by first waveform difference inversion from the seismic-based de-baselined medium and low frequency second sequence boundary curve comprises:

and (3) performing well-shock calibration and first-time waveform difference inversion on the middle-low frequency second sequence boundary curve after the base line removal processing based on the earthquake to obtain a second sequence boundary data body based on the earthquake.

6. The method as recited in claim 1, further comprising:

performing reverse processing on the gamma curve at the target well point to obtain a reverse gamma curve;

performing trending treatment on the reverse gamma curve to obtain a trended gamma curve;

and (3) removing the base line of the gamma curve after trending treatment to obtain a gamma curve reflecting the boundary of the high-frequency sequence.

7. The method as recited in claim 1, further comprising:

and performing baseline removal treatment on the porosity curve at the target well point to obtain a porosity curve reflecting the boundary of the high-frequency sequence.

8. The method of claim 1, wherein the blending with the gamma curve reflecting the high frequency sequence boundary, the porosity curve reflecting the high frequency sequence boundary, the carbonate grain curve reflecting the high frequency sequence boundary, the carbonate suture line frequency curve reflecting the sequence boundary, the asphalt content curve reflecting the carbonate of the sequence boundary, and the dolomite content curve reflecting the carbonate of the sequence boundary at the target well point, respectively, to obtain the high frequency tenth sequence boundary curve based on the earthquake comprises:

Discretizing to form a plurality of column data points by taking the third-order boundary curve based on the earthquake as a basic curve, discretizing the gamma curve reflecting the high-frequency sequence boundary to form a plurality of column data points after normalization treatment, superposing the column data points on a plurality of column data points corresponding to the third-order boundary curve based on the earthquake of the corresponding depth point, and fusing the column data points into a fourth high-frequency sequence boundary curve based on the earthquake;

discretizing to form a plurality of column data points by taking a fourth time sequence boundary curve based on the earthquake as a basic curve, discretizing to form a plurality of column data points after carrying out normalization processing on a porosity curve reflecting a high-frequency sequence boundary, superposing the column data points on a plurality of column data points corresponding to the fourth time sequence boundary curve based on the earthquake of a corresponding depth point, and fusing the column data points into a fifth high-frequency sequence boundary curve based on the earthquake;

discretizing a carbonate particle curve reflecting a high-frequency sequence boundary into a plurality of column data points by taking a fifth-order sequence boundary curve based on the earthquake as a basic curve, discretizing the carbonate particle curve reflecting the high-frequency sequence boundary into a plurality of column data points after normalization treatment, and then superposing the column data points on a plurality of column data points corresponding to the fifth-order sequence boundary curve based on the earthquake of a corresponding depth point to fuse the column data points into a sixth high-frequency sequence boundary curve based on the earthquake; wherein the carbonate particle curve reflecting the high frequency sequence boundary is a curve formed by quantifying carbonate particles into a plurality of column data points from granular limestone to argillaceous limestone;

Discretizing a carbonate texture curve reflecting a high-frequency sequence boundary into a plurality of column data points by taking a sixth-order sequence boundary curve based on the earthquake as a basic curve, discretizing the carbonate texture curve reflecting the high-frequency sequence boundary into a plurality of column data points after normalization treatment, and then superposing the column data points on a plurality of column data points corresponding to the sixth-order sequence boundary curve based on the earthquake of the corresponding depth point to fuse the column data points into a seventh high-frequency sequence boundary curve based on the earthquake; wherein, the carbonate rock texture curve reflecting the high-frequency sequence boundary is a curve formed by quantifying the rock texture structure of the carbonate rock into a plurality of column data points;

discretizing a series of data points by taking a seventh layer sequence boundary curve based on the earthquake as a basic curve, discretizing a suture line frequency curve of carbonate rock reflecting a high-frequency layer sequence boundary into a plurality of series of data points after normalization treatment, superposing the series of data points on a plurality of series of data points corresponding to the seventh layer sequence boundary curve based on the earthquake of a corresponding depth point, and fusing the series of data points into an eighth high-frequency layer sequence boundary curve based on the earthquake; wherein, the suture frequency curve of the carbonate rock reflecting the high-frequency sequence boundary is a curve formed by quantifying the suture frequency of the carbonate rock into a plurality of column data points;

Discretizing to form a plurality of column data points by taking an eighth sequence boundary curve based on the earthquake as a basic curve, discretizing to form a plurality of column data points after normalizing a carbonate rock asphalt content curve reflecting the sequence boundary, and then superposing the column data points on a plurality of column data points corresponding to the eighth sequence boundary curve based on the earthquake of the corresponding depth points to fuse the column data points into a ninth high-frequency sequence boundary curve based on the earthquake; wherein, the asphalt content curve of the carbonate rock reflecting the layer sequence boundary is a curve formed by quantifying the asphalt content of the carbonate rock into a plurality of column data points;

discretizing to form a plurality of column data points by taking a ninth-order layer sequence boundary curve based on the earthquake as a basic curve, discretizing to form a plurality of column data points after carrying out normalization treatment on a dolomite content curve of carbonate rock reflecting the layer sequence boundary, superposing the column data points on a plurality of column data points corresponding to the ninth-order layer sequence boundary curve based on the earthquake of a corresponding depth point, and fusing the column data points into a tenth high-frequency layer sequence boundary curve based on the earthquake; wherein, the dolomite content curve of the carbonate rock reflecting the sequence boundary is a curve formed by quantifying the dolomite to limestone content of the carbonate rock into a plurality of column data points.

9. A high frequency seismic sequence boundary determination apparatus for a sequence of layers, comprising:

the geological model grid calculation module is used for calculating geological model grids according to the three-dimensional seismic data volume of the target interval of the target area;

the relative geologic time model generation module is used for generating a relative geologic time model according to the geologic model grid;

the layer sequence thickness boundary data body extraction module is used for processing and analyzing the relative geologic time model and extracting an initial medium-low frequency first layer sequence boundary data body based on earthquake;

the system comprises a seismic-based interval boundary curve, a first-frequency-based interval boundary curve and a second-frequency-based interval boundary curve, wherein the seismic-based interval boundary curve is used for extracting a target interval at a target well point from an initial seismic-based middle-low frequency first interval boundary data volume; performing baseline removal processing on the middle-low frequency first sequence boundary curve based on the earthquake to obtain a middle-low frequency second sequence boundary curve after the baseline removal processing based on the earthquake; according to the medium-low frequency second sequence boundary curve after the baseline removal processing based on the earthquake, obtaining a second sequence boundary data body based on the earthquake through first waveform difference inversion; extracting a third sub-sequence boundary curve based on the earthquake of the target interval at the target well point from the second sequence boundary data body based on the earthquake;

The fusion processing module is used for fusing a third sequence boundary curve based on the earthquake with a gamma curve reflecting a high-frequency sequence boundary, a porosity curve reflecting the high-frequency sequence boundary, a carbonate particle curve reflecting the high-frequency sequence boundary, a carbonate texture curve reflecting the high-frequency sequence boundary, a carbonate suture line frequency curve reflecting the sequence boundary, an asphalt content curve reflecting the carbonate of the sequence boundary and an dolomite content curve reflecting the carbonate of the sequence boundary at a target well point to obtain a tenth sequence boundary curve based on the earthquake;

the high-frequency earthquake sequence boundary identification module is used for carrying out second waveform difference inversion on a target interval of the target area according to a high-frequency tenth sequence boundary curve based on the earthquake to obtain a Gao Pindi sequence boundary data body based on the earthquake; based on the Gao Pindi three-layer sequence boundary data body based on the earthquake, the sequence boundary of the high-frequency earthquake sequence stratum is identified.

10. The apparatus of claim 9, wherein the geologic model grid computing module is configured to:

determining an initial geological model grid according to the three-dimensional seismic data volume of the target interval of the target area;

Taking the initial geological model grid as the current geological model grid, and repeatedly executing the following steps until the matching condition of the current geological model grid and the three-dimensional seismic data volume meets the preset condition:

interactively correcting the association relation among the seismic horizons in the current geological model grid;

analyzing the coincidence condition of the corrected geological model grid and the three-dimensional seismic data volume;

when the fit condition does not meet the preset condition, optimizing parameters of the current geological model grid, and taking the optimized geological model grid as the current geological model grid.

11. The apparatus of claim 10, wherein the geologic model grid computing module is configured to:

when the three-dimensional seismic data volume has seismic horizon data, automatically tracking the seismic horizon data in the three-dimensional seismic data volume, and establishing a geological model grid by taking the tracked seismic horizon data as a constraint;

and when the three-dimensional seismic data volume does not contain the seismic horizon data, calculating an initial geological model grid according to the waveform similarity and the relative distance for at least one seed point in the three-dimensional seismic data volume of the target interval of the target area.

12. The apparatus of claim 9, wherein the relative geologic time model generation module is configured to:

Connecting and interpolating the surface element sheets of the geological model grid to obtain a processed geological model grid;

distributing relative geologic time to each pixel in the processed geologic model grid to generate an initial relative geologic time model;

a plurality of horizon stacks are extracted from the initial relative geologic time model to form a relative geologic time model represented by the plurality of horizon stacks.

13. The apparatus of claim 9, wherein the seismic-based interval boundary curve is specifically configured to:

and (3) performing well-shock calibration and first-time waveform difference inversion on the middle-low frequency second sequence boundary curve after the base line removal processing based on the earthquake to obtain a second sequence boundary data body based on the earthquake.

14. The apparatus of claim 9, wherein the fusion processing module is further configured to:

performing reverse processing on the gamma curve at the target well point to obtain a reverse gamma curve;

performing trending treatment on the reverse gamma curve to obtain a trended gamma curve;

and (3) removing the base line of the gamma curve after trending treatment to obtain a gamma curve reflecting the boundary of the high-frequency sequence.

15. The apparatus of claim 9, wherein the fusion processing module is further configured to:

And performing baseline removal treatment on the porosity curve at the target well point to obtain a porosity curve reflecting the boundary of the high-frequency sequence.

16. The apparatus of claim 9, wherein the fusion processing module is specifically configured to:

discretizing to form a plurality of column data points by taking the third-order boundary curve based on the earthquake as a basic curve, discretizing the gamma curve reflecting the high-frequency sequence boundary to form a plurality of column data points after normalization treatment, superposing the column data points on a plurality of column data points corresponding to the third-order boundary curve based on the earthquake of the corresponding depth point, and fusing the column data points into a fourth high-frequency sequence boundary curve based on the earthquake;

discretizing to form a plurality of column data points by taking a fourth time sequence boundary curve based on the earthquake as a basic curve, discretizing to form a plurality of column data points after carrying out normalization processing on a porosity curve reflecting a high-frequency sequence boundary, superposing the column data points on a plurality of column data points corresponding to the fourth time sequence boundary curve based on the earthquake of a corresponding depth point, and fusing the column data points into a fifth high-frequency sequence boundary curve based on the earthquake;

discretizing a carbonate particle curve reflecting a high-frequency sequence boundary into a plurality of column data points by taking a fifth-order sequence boundary curve based on the earthquake as a basic curve, discretizing the carbonate particle curve reflecting the high-frequency sequence boundary into a plurality of column data points after normalization treatment, and then superposing the column data points on a plurality of column data points corresponding to the fifth-order sequence boundary curve based on the earthquake of a corresponding depth point to fuse the column data points into a sixth high-frequency sequence boundary curve based on the earthquake; wherein the carbonate particle curve reflecting the high frequency sequence boundary is a curve formed by quantifying carbonate particles into a plurality of column data points from granular limestone to argillaceous limestone;

Discretizing a carbonate texture curve reflecting a high-frequency sequence boundary into a plurality of column data points by taking a sixth-order sequence boundary curve based on the earthquake as a basic curve, discretizing the carbonate texture curve reflecting the high-frequency sequence boundary into a plurality of column data points after normalization treatment, and then superposing the column data points on a plurality of column data points corresponding to the sixth-order sequence boundary curve based on the earthquake of the corresponding depth point to fuse the column data points into a seventh high-frequency sequence boundary curve based on the earthquake; wherein, the carbonate rock texture curve reflecting the high-frequency sequence boundary is a curve formed by quantifying the rock texture structure of the carbonate rock into a plurality of column data points;

discretizing a series of data points by taking a seventh layer sequence boundary curve based on the earthquake as a basic curve, discretizing a suture line frequency curve of carbonate rock reflecting a high-frequency layer sequence boundary into a plurality of series of data points after normalization treatment, superposing the series of data points on a plurality of series of data points corresponding to the seventh layer sequence boundary curve based on the earthquake of a corresponding depth point, and fusing the series of data points into an eighth high-frequency layer sequence boundary curve based on the earthquake; wherein, the suture frequency curve of the carbonate rock reflecting the high-frequency sequence boundary is a curve formed by quantifying the suture frequency of the carbonate rock into a plurality of column data points;

Discretizing to form a plurality of column data points by taking an eighth sequence boundary curve based on the earthquake as a basic curve, discretizing to form a plurality of column data points after normalizing a carbonate rock asphalt content curve reflecting the sequence boundary, and then superposing the column data points on a plurality of column data points corresponding to the eighth sequence boundary curve based on the earthquake of the corresponding depth points to fuse the column data points into a ninth high-frequency sequence boundary curve based on the earthquake; wherein, the asphalt content curve of the carbonate rock reflecting the layer sequence boundary is a curve formed by quantifying the asphalt content of the carbonate rock into a plurality of column data points;

discretizing to form a plurality of column data points by taking a ninth-order layer sequence boundary curve based on the earthquake as a basic curve, discretizing to form a plurality of column data points after carrying out normalization treatment on a dolomite content curve of carbonate rock reflecting the layer sequence boundary, superposing the column data points on a plurality of column data points corresponding to the ninth-order layer sequence boundary curve based on the earthquake of a corresponding depth point, and fusing the column data points into a tenth high-frequency layer sequence boundary curve based on the earthquake; wherein, the dolomite content curve of the carbonate rock reflecting the sequence boundary is a curve formed by quantifying the dolomite to limestone content of the carbonate rock into a plurality of column data points.

17. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of claims 1 to 8 when executing the computer program.

18. 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 method of any of claims 1 to 8.

19. A computer program product, characterized in that the computer program product comprises a computer program which, when executed by a processor, implements the method of any of claims 1 to 8.