Method and apparatus for determining sequence boundaries in high-frequency seismic sequence stratigraphy

Abstract

A method for determining sequence boundaries in high-frequency seismic sequence stratigraphy. The method comprises: generating a relative geological time model on the basis of a three-dimensional seismic data volume; extracting an initial seismic-based medium-low frequency first sequence boundary data volume; extracting a seismic-based medium-low frequency first sequence boundary curve; then, performing baseline removal processing; then, obtaining a seismic-based second sequence boundary data volume by means of a first waveform difference inversion; then, extracting a seismic-based third sequence boundary curve to serve as a basic curve and fusing same with a gamma curve, a porosity curve, a carbonate rock grain curve and a carbonate rock texture curve, that reflect high-frequency sequence boundaries, and a carbonate rock stylolite frequency curve, a carbonate rock asphalt content curve and a carbonate rock dolomite content curve that reflect sequence boundaries; then, performing a second waveform difference inversion to obtain a seismic-based high-frequency third sequence boundary data volume; and then identifying high-frequency seismic sequence boundaries. Thus, high precision is achieved. The present application further relates to an apparatus for determining sequence boundaries in high-frequency seismic sequence stratigraphy.

Description

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

[0001] Related applications

[0002] This application claims priority to Chinese patent application No. 202311734898.0 filed on December 15, 2023, and cites the contents disclosed in the above patent application as part of this application. Technical Field

[0003] The present application relates to the technical field of oil and gas exploration and development in carbonate strata, and seismic division of high-frequency sea-level change cyclic sequence bodies in strata, and in particular to a method and device for determining sequence boundaries of high-frequency seismic sequence strata. Background Art

[0004] This section is intended to provide a background or context to the embodiments of the present application that are recited in the claims. No admission is made that the description herein is prior art by virtue of its inclusion in this section.

[0005] A stratigraphic unit comprises strata deposited in all genetically related depositional environments during a complete base-level cycle. Semi-cycle boundaries within a genetic sequence occur at transitions from base-level rise to fall or vice versa. In different paleogeographic settings, these transitions manifest as stratigraphic discontinuities or conformable strata recording increases or decreases in accommodation space, respectively, often forming sequence boundaries.

[0006] In sequence stratigraphy, sequences are typically divided into levels 1 to 6. Levels 1 to 3 correspond to tectonic sedimentary cycles and are low-frequency sequences (cycles), corresponding to megasequences, supersequences, and sequences, respectively. Levels 4 to 6 correspond to climatic sedimentary cycles and are high-frequency sequences (cycles), also known as Milankovitch cycles, corresponding to parasequence sets, parasequences, and rhythmic beds, respectively. Sequence division is fundamental in sequence stratigraphy, and sequence boundary identification is crucial. Common sequence boundaries include unconformities, transgressive onlaps, flooding unconformities, paleokarstification surfaces, volcanic event surfaces, and lithologic transition surfaces. Sequences of different levels correspond to different levels of sequence boundaries. For example, regional unconformities and tectonic transition surfaces often correspond to low-frequency sequences, while lithologic transition surfaces and sedimentary unconformities often correspond to high-frequency sequences.

[0007] Commonly used geological methods for determining sequence boundaries include lithologic and lithofacies variations, core observations, INPEFA curves, and wavelet-induced spectrum analysis. However, in areas primarily studied for carbonate rocks, the lithofacies include grainstone limestone, grain-bearing micrite limestone, micrite limestone (containing grains), dolomite and dolomitic limestone, and gypsum-bearing micrite limestone. Multiple thin-bedded shoals and bioclastic shoals are well developed, and multiple stages and sets of shoal bodies have developed due to multiple cycles. For highly heterogeneous carbonate formations, existing sequence boundary determination methods typically only reflect the boundaries of low-frequency sequence formations and fail to identify the higher-frequency sequence stratigraphic framework requirements. Current technical solutions are unable to effectively determine sequence boundaries in high-frequency seismic sequence formations.

[0008] Summary of the Invention

[0009] In a first aspect, embodiments of the present application provide a method for determining sequence boundaries of high-frequency seismic sequence stratigraphy, which can establish a more accurate high-frequency sequence boundary volume, thereby making the sequence boundaries of the identified high-frequency seismic sequence stratigraphy more accurate. The method includes:

[0010] Calculate the geological model grid based on the 3D seismic data volume of the target layer segment in the target area;

[0011] Generate a relative geological age model based on the geological model grid;

[0012] Process and analyze the relative geological age model to extract the initial seismic-based low- and medium-frequency first sequence boundary data volume;

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

[0014] Perform baseline removal 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.

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

[0016] 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;

[0017] The third sequence boundary curve based on seismic data is used as the base 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 tenth high-frequency sequence boundary curve based on seismic data.

[0018] 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 in the target area to obtain the high-frequency third-order sequence boundary data volume based on seismic data.

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

[0020] In a second aspect, an embodiment of the present application further provides a device for determining sequence boundaries of high-frequency seismic sequence stratigraphy, which can establish a more accurate high-frequency sequence boundary volume, thereby making the sequence boundaries of the identified high-frequency seismic sequence stratigraphy more accurate. The device includes:

[0021] A geological model grid calculation module is used to calculate the geological model grid based on the three-dimensional seismic data volume of the target layer segment in the target area;

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

[0023] Sequence thickness boundary data volume extraction module, 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;

[0024] Seismic-based sequence boundary curves are 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 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; and 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;

[0025] a fusion processing module for fusing the third-order seismic sequence boundary curve as a base 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 grain 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 the carbonate rock reflecting the sequence boundary, and a dolomite content curve of the carbonate rock reflecting the sequence boundary, to obtain a tenth-order seismic high-frequency sequence boundary curve;

[0026] 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 the seismic data to obtain the high-frequency third sequence boundary data volume based on the seismic data; based on the high-frequency third sequence boundary data volume based on the seismic data, the sequence boundaries of the high-frequency seismic sequence strata are identified.

[0027] In a third aspect, an embodiment of the present application further provides a computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the above-mentioned method for determining the sequence boundaries of high-frequency seismic sequence strata is implemented.

[0028] In a fourth aspect, an embodiment of the present application further provides a computer-readable storage medium storing a computer program, which implements the above-mentioned method for determining the sequence boundaries of high-frequency seismic sequence strata when executed by a processor.

[0029] In a fifth aspect, an embodiment of the present application further provides a computer program product, which includes a computer program. 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.

[0030] In an embodiment of the present application, 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 the 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; and from the seismic-based second sequence boundary data volume, the target layer segment at the target well point is extracted. a third seismic-based sequence boundary curve; taking the third seismic-based sequence boundary curve as the base curve, the curves are 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 tenth high-frequency seismic-based sequence boundary curve; based on the tenth high-frequency seismic-based sequence boundary curve, a second waveform difference inversion is performed on the target layer section in the target area to obtain the third high-frequency seismic-based sequence boundary data volume; based on the third high-frequency seismic-based sequence boundary data volume, the sequence boundary of the high-frequency seismic sequence stratum is identified. Compared with the existing sequence interface determination method, the embodiment of the present application extracts the seismic-based medium- and low-frequency first sequence boundary curves of the target layer at the target well point, and then 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. Finally, an inversion based on waveform differences is performed to obtain a high-frequency sequence boundary data body. Based on the high-frequency sequence boundary data body, the sequence boundary of the high-frequency seismic sequence stratum is identified. Within the high-frequency sequence boundary body, the larger the value, the closer it is to the sequence boundary. The above process integrates the gamma curve and porosity curve at the target well point, and the inverted sequence boundary characteristics 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 stratum is more accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the following briefly introduces the drawings required for the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present application. For those skilled in the art, other drawings can be obtained based on these drawings without creative work. In the drawings:

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

[0033] FIG2 is a flow chart of calculating a geological model grid in an embodiment of the present application;

[0034] FIG3 is a flow chart of generating a relative geological age model in an embodiment of the present application;

[0035] FIG4 is a flow chart of multi-curve fusion in an embodiment of the present application;

[0036] FIG5 is an example diagram of a relative geological age model calculated based on the three-dimensional seismic data volume of the region in an embodiment of the present application;

[0037] FIG6 is an example diagram of a well-connected cross section obtained by analyzing and extracting an initial seismic-based medium- and low-frequency first sequence boundary data volume in an embodiment of the present application from a relative geological age model;

[0038] FIG7 is an example diagram of a sequence boundary curve based on earthquakes in an embodiment of the present application;

[0039] 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 application;

[0040] FIG9 is a comparison diagram of gamma curves before and after processing in an embodiment of the present application;

[0041] FIG10 is a comparison diagram of the porosity curves before and after processing in an embodiment of the present application;

[0042] FIG11 is an example of the process of fusing and obtaining high-frequency sequence boundary curves in an embodiment of the present application;

[0043] FIG12 is an example of a well-connected cross section of a high-frequency sequence boundary data volume in an embodiment of the present application;

[0044] 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 application;

[0045] FIG14 is a schematic diagram of a computer device in an embodiment of the present application. DETAILED DESCRIPTION

[0046] In order to make the purpose, technical solutions and advantages of the embodiments of the present application more clear, the embodiments of the present application are further described in detail below with reference to the accompanying drawings. Here, the illustrative embodiments of the present application and their descriptions are used to explain the present application, but are not intended to limit the present application.

[0047] 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 application, comprising:

[0048] Step 101, calculating a geological model grid based on a 3D seismic data volume of a target layer segment in a target area;

[0049] Step 102, generating a relative geological age model based on the geological model grid;

[0050] Step 103: Process and analyze the relative geological age model to extract the initial seismic-based low- and medium-frequency first sequence boundary data volume;

[0051] 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;

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

[0053] Step 106, obtaining a seismic-based second sequence boundary data volume through the first waveform difference inversion based on the mid-low frequency second sequence boundary curve after the seismic-based baseline removal process;

[0054] Step 107, extracting a 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;

[0055] Step 108: The third-order seismic sequence boundary curve is used as a base curve, and is fused with a gamma curve reflecting the high-frequency sequence boundary at the target well point, a porosity curve reflecting the high-frequency sequence boundary, a carbonate rock grain 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 the carbonate rock reflecting the sequence boundary, and a dolomite content curve of the carbonate rock reflecting the sequence boundary, to obtain a tenth-order seismic high-frequency sequence boundary curve.

[0056] Step 109: performing 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 data to obtain a high-frequency third-order sequence boundary data volume based on seismic data.

[0057] Step 110 : Identify the sequence boundaries of the high-frequency seismic sequence strata based on the high-frequency third sequence boundary data volume based on the seismic.

[0058] In an embodiment of the present application, 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 the 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; and from the seismic-based second sequence boundary data volume, the target layer segment at the target well point is extracted. a third seismic-based sequence boundary curve; taking the third seismic-based sequence boundary curve as the base curve, the curves are 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 tenth high-frequency seismic-based sequence boundary curve; based on the tenth high-frequency seismic-based sequence boundary curve, a second waveform difference inversion is performed on the target layer section in the target area to obtain the third high-frequency seismic-based sequence boundary data volume; based on the third high-frequency seismic-based sequence boundary data volume, the sequence boundary of the high-frequency seismic sequence stratum is identified. Compared with the existing sequence interface determination method, the embodiment of the present application extracts the seismic-based medium- and low-frequency first sequence boundary curves of the target layer at the target well point, and then 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. Finally, an inversion based on waveform differences is performed to obtain a high-frequency sequence boundary data body. Based on the high-frequency sequence boundary data body, the sequence boundary of the high-frequency seismic sequence stratum is identified. Within the high-frequency sequence boundary body, the larger the value, the closer it is to the sequence boundary. The above process integrates the gamma curve and porosity curve at the target well point, and the inverted sequence boundary characteristics 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 stratum is more accurate.

[0059] In step 101, a geological model grid is calculated based on a 3D seismic data volume of a target layer segment in a target area;

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

[0061] Referring to FIG2 , based on the 3D seismic data volume of the target layer in the target area, the geological model grid is calculated, including:

[0062] Step 201 : Determine an initial geological model grid based on the 3D seismic data volume of the target layer segment in the target area. Two methods for determining the initial geological model grid are provided herein.

[0063] In one embodiment, determining an initial geological model grid based on a 3D seismic data volume of a target layer in a target area includes:

[0064] When seismic layer data exists in the 3D seismic data volume, the seismic layer data in the 3D seismic data volume is automatically tracked, and a geological model grid is established using the tracked seismic layer data as a constraint; wherein the seismic layer data already exists in the 3D seismic data volume;

[0065] When no seismic horizon data exists within the 3D seismic data volume, an initial geological model grid is calculated based on waveform similarity and relative distance for at least one seed point within the 3D seismic data volume within the target layer of the target area. This calculation can employ algorithms such as those based on boundary control and local mapping, which are not limited here.

[0066] Step 202: Using the initial geological model grid as the current geological model grid, the following steps are repeated until the current geological model grid and the 3D seismic data volume meet the preset conditions:

[0067] Step 2021: interactively modify the association between seismic layers in the current geological model grid; each modification will affect the links between nodes in the model grid;

[0068] Step 2022: Analyze the consistency between the modified geological model grid and the 3D seismic data volume; in specific implementation, the analysis can be performed through preview;

[0069] Step 2023: When the matching condition 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.

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

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

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

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

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

[0075] Step 303: extract multiple stratigraphic stacks from the initial relative geological age model to form a relative geological age model represented by the multiple stratigraphic stacks. Typically, there are tens of thousands of stratigraphic stacks, and a relative geological age model represented by tens of thousands of stratigraphic stacks is more accurate.

[0076] Step 103 processes and analyzes the relative geological age model to extract the initial seismic-based medium- and low-frequency first sequence boundary data volume. The sequence thickness boundary data volume is an attribute data volume that can reflect the stratigraphic sequence boundary. It can reflect the changes in the vertical and horizontal sequence boundaries and is combined with seismic data 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 seismic-geology.

[0077] 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;

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

[0079] The seismic-based low-frequency second sequence boundary curve after baseline removal is subjected to well-seismic calibration and the first waveform difference inversion to obtain the seismic-based second sequence boundary data volume.

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

[0081] In step 108, the third-order seismic sequence boundary curve is used as a base curve and is fused with a gamma curve reflecting the high-frequency sequence boundary at the target well point, a porosity curve reflecting the high-frequency sequence boundary, a carbonate rock grain 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 the carbonate rock reflecting the sequence boundary, and a dolomite content curve of the carbonate rock reflecting the sequence boundary to obtain a tenth-order seismic high-frequency sequence boundary curve.

[0082] Referring to Figure 4, the specific fusion steps include:

[0083] Step 401: Using the third-order seismic sequence boundary curve as a base curve, the curve is discretized into a plurality of columns of data points. The gamma curve reflecting the high-frequency sequence boundary is normalized and discretized into a plurality of columns of data points. The data points are then superimposed on the plurality of columns of data points corresponding to the third-order seismic sequence boundary curve at the corresponding depth point, and the resultant curves are merged into a fourth-order seismic high-frequency sequence boundary curve.

[0084] In one embodiment, the method for determining sequence boundaries of high-frequency seismic sequence strata further includes:

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

[0086] Detrending the inverse gamma curve to obtain a detrended gamma curve;

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

[0088] Step 402: Using the fourth-order seismic sequence boundary curve as a base curve, the curve is discretized into a plurality of columns of data points. The porosity curve reflecting the high-frequency sequence boundary is normalized and discretized into a plurality of columns of data points. The data points are then superimposed on the plurality of columns of data points corresponding to the fourth-order seismic sequence boundary curve at the corresponding depth point, and the resultant data is merged into a fifth-order seismic high-frequency sequence boundary curve.

[0089] In one embodiment, the method for determining sequence boundaries of high-frequency seismic sequence strata further includes:

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

[0091] Step 403: Using the fifth-order seismic sequence boundary curve as a base curve, the curve is discretized into a plurality of columns of data points. The carbonate rock grain curve reflecting the high-frequency sequence boundary is normalized and discretized into a plurality of columns of data points. The curve is then superimposed on the plurality of columns of data points corresponding to the fifth-order seismic sequence boundary curve at the corresponding depth point, and the resultant curve is merged into a sixth-order seismic high-frequency sequence boundary curve. 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 ranging from grain limestone to argillaceous limestone.

[0092] For example, multiple column data points quantified from grain 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 grains and the higher the mud content. The larger the positive value, the larger the carbonate rock grains and the higher the carbonate rock grain content.

[0093] 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 carbonate rock sequence boundary information.

[0094] Step 404: Using the sixth seismic sequence boundary curve as a base curve, the curve is discretized into a plurality of columns of data points. The carbonate rock texture curve reflecting the high-frequency sequence boundary is normalized and discretized into a plurality of columns of data points. The curve is then superimposed on the plurality of columns of data points corresponding to the sixth seismic sequence boundary curve at the corresponding depth point, and the resultant curve is merged into a seventh seismic high-frequency sequence boundary curve. The carbonate rock texture curve reflecting the high-frequency sequence boundary is a curve formed by quantizing the rock texture structure of the carbonate rock into a plurality of columns of data points.

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

[0096] Among them, the carbonate rock texture curve represents the structure of the carbonate rock 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.

[0097] Step 405: Using the seventh-order seismic sequence boundary curve as a base curve, the curve is discretized into a plurality of columns of data points. The suture frequency curve of the carbonate rock reflecting the high-frequency sequence boundary is normalized and discretized into a plurality of columns of data points. The suture frequency curve is then superimposed on the plurality of columns of data points corresponding to the seventh-order seismic sequence boundary curve at the corresponding depth point, and the resultant curve is merged into the eighth-order seismic high-frequency sequence boundary curve. The suture frequency curve of the carbonate rock reflecting the high-frequency sequence boundary is a curve formed by quantizing the suture frequency of the carbonate rock into a plurality of columns of data points.

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

[0099] 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 greater the number of sutures, the stronger the progradation characteristics caused by sea regression and the more obvious the sequence boundaries caused by coarse-grained sedimentation.

[0100] Step 406: Using the eighth-order seismic sequence boundary curve as a base curve, the asphalt content curve of the carbonate rock reflecting the sequence boundary is discretized into a plurality of columns of data points. The asphalt content curve of the carbonate rock reflecting the sequence boundary is normalized and discretized into a plurality of columns of data points. The asphalt content curve is then superimposed on the plurality of columns of data points corresponding to the eighth-order seismic sequence boundary curve at the corresponding depth point, thereby merging the curve into a ninth-order high-frequency seismic sequence boundary curve. The asphalt content curve of the carbonate rock reflecting the sequence boundary is a curve formed by quantizing the asphalt content of the carbonate rock into a plurality of columns of data points.

[0101] For example, the bitumen content of carbonate rocks is quantified on a scale of 0.1 to 0.9, with higher values ​​indicating more bitumen content;

[0102] The asphalt content curve of carbonate rocks reflects the changes in temperature, pressure, and composition during the carbonate reservoir formation process, which in turn causes changes in sequence boundaries. Therefore, integrating the asphalt content curve of carbonate rocks that reflects sequence boundaries can help accurately identify sequence boundaries.

[0103] Step 407: Using the ninth order seismic sequence boundary curve as a base curve, the curve is discretized into a plurality of columns of data points. The dolomite content curve of the carbonate rock reflecting the sequence boundary is normalized and discretized into a plurality of columns of data points. The curve is then superimposed on the plurality of columns of data points corresponding to the ninth order seismic sequence boundary curve at the corresponding depth point, and merged into a tenth order seismic high-frequency sequence boundary curve. 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 columns of data points.

[0104] For example, the dolomite to limestone content of carbonate rocks is quantified as 0.1 to 0.9, where the larger the value, the more dolomite content. Among them, the higher the dolomite content, the more exposed the carbonate rock deposition environment is, the shallower the sea water is, 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.

[0105] The final tenth high-frequency sequence boundary curve based on seismic data can reflect the information of high-frequency sequence boundaries.

[0106] In step 109, based on the high-frequency tenth-order sequence boundary curve based on seismic data, a second waveform difference inversion is performed on the target layer segment in the target area to obtain a high-frequency third-order sequence boundary data volume based on seismic data. The high-frequency sequence boundary data volume at this time is clearer and more accurate.

[0107] 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 data. 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.

[0108] The method proposed in the embodiment of the present application 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. In particular, the research on complex carbonate episodic sequence cycles, sedimentary patterns and corresponding multi-stage diagenetic evolution is weak. 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 storage control mechanism of multiple geological factors. Therefore, the method of the present application has broad application prospects.

[0109] A specific example is given below to illustrate the specific application of the method proposed in this application.

[0110] Taking a region where carbonate oil and gas reservoirs are mainly developed as an example, FIG5 is an example diagram of a relative geological age model calculated based on the three-dimensional seismic data body of the region in an embodiment of the present application, wherein (a) in FIG5 is a three-dimensional seismic data body, (b) in FIG5 is a geological model grid, and (c) in FIG5 is a relative geological age model. FIG6 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 application 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, onlap), 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.

[0111] Figure 7 shows an example of a seismic-based sequence boundary curve in an embodiment of the present application, wherein (a) in Figure 7 shows a first low- to medium-frequency seismic-based sequence boundary curve, and (b) in Figure 7 shows the curve after baseline removal from the initial sequence boundary curve. Figure 8 shows an example of a well-connected cross-section of a second seismic-based sequence boundary data volume obtained through first waveform difference inversion in an embodiment of the present application. The final extracted third seismic-based sequence boundary curve highlights the ability of sequence thickness boundary curves obtained from seismic data to identify sequence boundaries.

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

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

[0114] Figure 11 is an example of the process of fusing and obtaining a high-frequency tenth-order sequence boundary curve based on seismic data in an embodiment of the present application, wherein (a) in Figure 11 is a third-order sequence boundary curve based on seismic data, (b) in Figure 11 is a processed sequence boundary curve, (c) in Figure 11 is an initial gamma curve, (d) in Figure 11 is a gamma curve reflecting a high-frequency sequence boundary, (e) in Figure 11 is an initial porosity curve, (f) in Figure 11 is a porosity curve reflecting a high-frequency sequence boundary, and (g) in Figure 11 is a fifth-order high-frequency sequence boundary curve based on seismic data obtained by fusion. Subsequently, other curves can be further fused to finally obtain a tenth-order high-frequency sequence boundary curve based on seismic data. Figure 12 is an example of a well-connected section of a high-frequency third-order sequence boundary data body based on seismic data in an embodiment of the present application. At this time, the high-frequency third-order sequence boundary data body based on seismic data is clearer and more accurate. Finally, using the global thinking concept based on the high-frequency sequence interface data body, the sequence boundary of the high-frequency seismic sequence strata with higher accuracy can be automatically identified.

[0115] The embodiment of the present application further proposes 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 repeated here.

[0116] 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 application, comprising:

[0117] The geological model grid calculation module 1301 is used to calculate the geological model grid based on the 3D seismic data volume of the target layer segment in the target area;

[0118] The relative geological age model generation module 1302 is used to generate a relative geological age model based on the geological model grid;

[0119] The sequence thickness boundary data volume extraction module 1303 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;

[0120] Seismic-based sequence boundary curve 1304 is used to extract a 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; 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 first waveform difference inversion based on the seismic-based medium-low frequency second sequence boundary curve after baseline removal; and extract a 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;

[0121] A fusion processing module 1305 is configured to fuse the third-order seismic sequence boundary curve as a base 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 grain 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 the carbonate rock reflecting the sequence boundary, and a dolomite content curve of the carbonate rock reflecting the sequence boundary, to obtain a tenth-order seismic high-frequency sequence boundary curve.

[0122] 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 seismic sequence boundary curve to obtain a high-frequency third-order seismic sequence boundary data volume; and identify the sequence boundaries of the high-frequency seismic sequence strata based on the high-frequency third-order seismic sequence boundary data volume.

[0123] In one embodiment, the geological model grid calculation module is specifically used to:

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

[0125] 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:

[0126] Interactively modify the correlation between seismic horizons in the current geological model grid;

[0127] Analyze the consistency between the revised geological model grid and the 3D seismic data volume;

[0128] When the coincidence 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.

[0129] In one embodiment, the geological model grid calculation module is specifically used to:

[0130] When there is seismic layer data in the 3D seismic data volume, the seismic layer data in the 3D seismic data volume is automatically tracked, and a geological model grid is established with the tracked seismic layer data as a constraint;

[0131] When there is no seismic layer data in the 3D seismic data volume, an initial geological model grid is calculated based on waveform similarity and relative distance for at least one seed point in the 3D seismic data volume of the target layer segment in the target area.

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

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

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

[0135] 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.

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

[0137] The seismic-based low-frequency second sequence boundary curve after baseline removal is subjected to well-seismic calibration and the first waveform difference inversion to obtain the seismic-based second sequence boundary data volume.

[0138] In one embodiment, the fusion processing module is specifically configured to:

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

[0140] Detrending the inverse gamma curve to obtain a detrended gamma curve;

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

[0142] In one embodiment, the fusion processing module is specifically configured to:

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

[0144] In one embodiment, the fusion processing module is specifically configured to:

[0145] The third-order sequence boundary curve based on seismic data is used as the base curve and 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. The curves are then superimposed on the multiple columns of data points corresponding to the third-order sequence boundary curve based on seismic data at the corresponding depth points to form the fourth-order high-frequency sequence boundary curve based on seismic data.

[0146] The fourth-order seismic sequence boundary curve is used as the base curve and 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. The data points are then superimposed on the multiple columns of data points corresponding to the fourth-order seismic sequence boundary curve at the corresponding depth point to form the fifth-order seismic high-frequency sequence boundary curve.

[0147] The fifth-order seismic sequence boundary curve is used as a base 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. The curve is then superimposed on the multiple columns of data points corresponding to the fifth-order seismic sequence boundary curve at the corresponding depth point to form the sixth-order seismic high-frequency sequence boundary curve. 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 ranging from grain limestone to argillaceous limestone.

[0148] The sixth sequence boundary curve based on seismic data is used as a base 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. The curve is 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 form the seventh high-frequency sequence boundary curve based on seismic data. The carbonate rock texture curve reflecting the high-frequency sequence boundary is a curve formed by quantizing the rock texture structure of the carbonate rock into multiple columns of data points.

[0149] The seventh-order seismic sequence boundary curve is used as a base curve and discretized into multiple columns of data points. The suture line frequency curve of the carbonate rock reflecting the high-frequency sequence boundary is normalized and discretized into multiple columns of data points. The curve is then superimposed on the multiple columns of data points corresponding to the seventh-order seismic sequence boundary curve at the corresponding depth point, and merged into the eighth-order seismic high-frequency sequence boundary curve. The suture line frequency curve of the carbonate rock reflecting the high-frequency sequence boundary is a curve formed by quantizing the suture line frequency of the carbonate rock into multiple columns of data points.

[0150] The eighth-order sequence boundary curve based on seismic data is used as a base curve and discretized into multiple columns of data points. The asphalt content curve of the carbonate rock reflecting the sequence boundary is normalized and discretized into multiple columns of data points. The curve is then superimposed on the multiple columns of data points corresponding to the eighth-order sequence boundary curve based on seismic data at the corresponding depth point, and fused into the ninth-order high-frequency sequence boundary curve based on seismic data. The asphalt content curve of the carbonate rock reflecting the sequence boundary is a curve formed by quantifying the asphalt content of the carbonate rock into multiple columns of data points.

[0151] The ninth sequence boundary curve based on seismic data is used as a base curve and discretized into multiple columns of data points. The dolomite content curve of the 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 data at the corresponding depth points to form the tenth high-frequency sequence boundary curve based on seismic data. 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 multiple columns of data points.

[0152] In summary, in the method and device proposed in the embodiment of the present application, 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 base 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 fused to obtain a tenth seismic-based high-frequency sequence boundary curve; based on the tenth seismic-based high-frequency sequence boundary curve, a second waveform difference inversion is performed on the target layer segment in the target area to obtain a third seismic-based high-frequency sequence boundary data volume; based on the third seismic-based high-frequency sequence boundary data volume, the sequence boundary of the high-frequency seismic sequence stratum is identified. Compared with the existing sequence interface determination method, the embodiment of the present application extracts the seismic-based medium- and low-frequency first sequence boundary curves of the target layer at the target well point, and then 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. Finally, an inversion based on waveform differences is performed to obtain a high-frequency sequence boundary data body. Based on the high-frequency sequence boundary data body, the sequence boundary of the high-frequency seismic sequence stratum is identified. Within the high-frequency sequence boundary body, the larger the value, the closer it is to the sequence boundary. The above process integrates the gamma curve and porosity curve at the target well point, and the inverted sequence boundary characteristics 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 stratum is more accurate.

[0153] An embodiment of the present application also provides a computer device. Figure 14 is a schematic diagram of the computer device in an embodiment of the present application. 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.

[0154] An embodiment of the present application further provides a computer-readable storage medium storing 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.

[0155] An embodiment of the present application further provides a computer program product, which includes a computer program. 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.

[0156] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0157] The present application is described with reference to the flow chart and / or block diagram of the method, device (system), and computer program product according to the embodiment of the present application. It should be understood that each flow process and / or box in the flow chart and / or block diagram and the combination of the flow process and / or box in the flow chart and / or block diagram can be realized 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 processing machine 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 realizing the function specified in one flow chart flow or multiple flows and / or one box or multiple boxes of the block diagram.

[0158] 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 operate in a specific manner, so that the instructions stored in the computer-readable memory produce a 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.

[0159] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, so that 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.

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

Claims

1. A method for determining sequence boundaries of high-frequency seismic sequence strata, characterized in that: include: 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 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, 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.

2. The method according to claim 1, characterized in that Based on the 3D seismic data volume of the target layer in the target area, the geological model grid is calculated, including: 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.

3. The method according to claim 2, characterized in that Based on the 3D seismic data volume of the target layer in the target area, the initial geological model grid is determined, including: 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; 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.

4. The method according to claim 1, characterized in that Generate a relative geological age model based on the geological model grid, including: 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.

5. The method according to claim 1, characterized in that: 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, including: 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.

6. The method according to claim 1, characterized in that Also includes: 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.

7. The method according to claim 1, characterized in that Also includes: The porosity curve at the target well point is de-baselined to obtain a porosity curve reflecting the high-frequency sequence boundary.

8. The method according to claim 1, characterized in that: The third sequence boundary curve based on seismic is used as the basic curve, and 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 respectively fused to obtain the high-frequency tenth sequence boundary curve based on seismic, including: 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 sequence boundary curve based on earthquake is used as the basic curve, which is discretized into multiple columns of data points to reflect the high-frequency layer. After the porosity curve of the sequence boundary is normalized and discretized into multiple columns of data points, it is 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-order sequence boundary curve based on earthquake is used as the basic curve, which is discretized into a plurality of columns of data points. The asphalt content curve of carbonate rock reflecting the 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 eighth-order sequence boundary curve based on earthquake at the corresponding depth point, and merged into the ninth-order high-frequency sequence boundary curve based on earthquake; 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 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.

9. A device for determining the sequence boundary of high-frequency seismic sequence strata, characterized in that: include: 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-order sequence boundary curve based on seismic, and obtain the high-frequency third-order sequence boundary data body based on seismic; based on the high-frequency third-order sequence boundary data body based on seismic, the sequence boundary of the high-frequency seismic sequence strata is identified.

10. The device according to claim 9, characterized in that 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.

11. The device according to claim 10, characterized in that The geological model grid calculation module is specifically used for: 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; 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.

12. The device according to claim 9, characterized in that The relative geological age model generation module is specifically used for: 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.

13. The device according to claim 9, characterized in that Seismic-based sequence boundary curves are specifically used for: 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.

14. The device according to claim 9, characterized in that The fusion processing module is also used to: 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.

15. The device according to claim 9, characterized in that The fusion processing module is also used to: The porosity curve at the target well point is de-baselined to obtain a porosity curve reflecting the high-frequency sequence boundary.

16. The device according to claim 9, characterized in that 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 the base curve, which is 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, and merged into the seventh high-frequency sequence boundary curve based on seismic data. Among them, the carbonate rock texture curve reflecting the high-frequency sequence boundary is the rock texture structure quantification of carbonate rock. A curve formed for 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-order sequence boundary curve based on earthquake is used as the basic curve, which is discretized into a plurality of columns of data points. The asphalt content curve of carbonate rock reflecting the 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 eighth-order sequence boundary curve based on earthquake at the corresponding depth point, and merged into the ninth-order high-frequency sequence boundary curve based on earthquake; 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 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.

17. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 8 is implemented.

18. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 8 is implemented.

19. A computer program product, characterized in that The computer program product comprises a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 8 is implemented.

Images