Quantitative analysis method and device for carbonate rock bench evolution

Abstract

The invention discloses a quantitative analysis method and device for carbonate rock bench evolution, and the method comprises the steps: analyzing outcrop data and logging data, building a sedimentary sequence, and carrying out the absolute time calibration of the sedimentary sequence; the sea level lifting cycle is analyzed, the relative sea level variation is obtained, and the relative sea level change rate is calculated; the tectonic settlement amount is obtained through a back-stripping method, and the tectonic settlement rate is calculated; respectively utilizing the relative sea level variation and the tectonic settlement to calculate the carbonate rock yield; the sum of the tectonic settlement rate and the relative sea level change rate is calculated, and the capacity space consumption rate is obtained; and quantitatively analyzing the bench evolution of the carbonate rock by using the difference between the carbonate rock yield and the tolerable space consumption rate. According to the method, the carbonate rock bench evolution can be quantitatively analyzed, a basis is provided for improving reservoir prediction precision, the exploration cost is reduced, and the exploration period is shortened.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas exploration, and in particular to a quantitative analysis method and device for carbonate platform evolution. Background Art

[0002] This section is intended to provide a background or context to the embodiments of the invention that are recited in the claims. No statement herein is admitted to be prior art by virtue of its inclusion in this section.

[0003] The evolution of carbonate platforms is influenced by a variety of geological factors, with climate, sea level, and tectonic activity being key. Climate and sea level fluctuations are controlled by astronomical cycles. Current research indicates that cyclostratigraphy reflects periodic changes in the sedimentary environment, while tectonic activity primarily controls changes in platform morphology and sedimentary space.

[0004] However, previous studies on carbonate platform evolution have mostly focused on qualitative descriptions, focusing primarily on the influence of various controlling factors. This qualitative approach struggles to accurately reveal the complex relationships and specific mechanisms underlying carbonate platform evolution, limiting our in-depth understanding and accurate prediction of carbonate platform evolution. Summary of the Invention

[0005] An embodiment of the present invention provides a quantitative analysis method for carbonate platform evolution, which is used to quantitatively analyze carbonate platform evolution and provide a basis for improving reservoir prediction accuracy, thereby reducing exploration costs and shortening exploration cycles. The method includes:

[0006] By analyzing outcrop data and well logging data, the sedimentary cycle units are identified, the sedimentary sequence is established, and the absolute time calibration of the sedimentary sequence is performed;

[0007] Identify and analyze sea-level cycles in sedimentary sequences, obtain relative sea-level changes at each absolute time point, and calculate relative sea-level change rates;

[0008] The tectonic settlement amount at each absolute time calibration point is obtained by backstripping method, and the tectonic settlement rate at each absolute time calibration point is calculated;

[0009] The carbonate rock yield at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point.

[0010] The sum of the tectonic subsidence rate and the relative sea level change rate at each absolute time calibration point was calculated to obtain the accommodation space consumption rate at each absolute time calibration point. The difference between the carbonate rock yield and the accommodation space consumption rate at each absolute time calibration point was used to quantitatively analyze the evolution of the carbonate platform.

[0011] The present invention also provides a quantitative analysis device for carbonate platform evolution, which is used to quantitatively analyze carbonate platform evolution and provide a basis for improving reservoir prediction accuracy, thereby reducing exploration costs and shortening exploration cycles. The device includes:

[0012] Sedimentary cycle analysis module is used to: identify sedimentary cycle units, establish sedimentary sequences, and perform absolute time calibration on sedimentary sequences by analyzing outcrop data and well logging data;

[0013] The relative sea level analysis module is used to identify and analyze sea level cycles in sedimentary sequences, obtain relative sea level changes at each absolute time point, and calculate relative sea level change rates.

[0014] The structural settlement analysis module is used to obtain the structural settlement amount at each absolute time calibration point through the backstripping method and calculate the structural settlement rate at each absolute time calibration point;

[0015] The carbonate rock yield analysis module is used to calculate the carbonate rock yield at each absolute time calibration point using the relative sea level change and tectonic subsidence at each absolute time calibration point;

[0016] The carbonate platform evolution analysis module is used to calculate the sum of the tectonic subsidence rate and the relative sea level change rate at each absolute time calibration to obtain the accommodation space consumption rate at each absolute time calibration; and use the difference between the carbonate rock yield and the accommodation space consumption rate at each absolute time calibration to quantitatively analyze the evolution of the carbonate platform.

[0017] An embodiment of the present invention also provides a computer device, including 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 quantitative analysis method for the evolution of carbonate platforms is implemented.

[0018] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the computer program implements the above-mentioned quantitative analysis method for carbonate platform evolution.

[0019] 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, it implements the above-mentioned quantitative analysis method for carbonate platform evolution.

[0020] Compared with the technical solutions for qualitatively describing the evolution of carbonate platforms in the existing technology, this method establishes a sedimentary sequence by analyzing outcrop data and logging data, and performs absolute time calibration on the sedimentary sequence; takes absolute time as the benchmark, integrates multi-source data, and improves the analysis accuracy of carbonate platform evolution; obtains relative sea level change through sedimentary cycle analysis, and calculates the relative sea level change rate; uses the backstripping method to obtain the tectonic subsidence, and calculates the tectonic subsidence rate; uses the relative sea level change and tectonic subsidence to calculate the carbonate rock yield at each absolute time calibration point; calculates the sum of the tectonic subsidence rate and the sea level change rate at each absolute time calibration point, and obtains the accommodative space consumption rate at each absolute time calibration point; uses the difference between the carbonate rock yield and the accommodative space consumption rate at each absolute time calibration point as a mathematical model for analyzing carbonate platform evolution, which can realize quantitative analysis of carbonate platform evolution, provide a basis for improving reservoir prediction accuracy, and thus reduce exploration costs and shorten exploration cycles. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] 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 the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative work. In the drawings:

[0022] Figure 1 This is a flow chart of a quantitative analysis method for carbonate platform evolution in an embodiment of the present invention;

[0023] Figure 2 This is an example diagram of a seismic profile across the Jiangliangping Trough in an embodiment of the present invention;

[0024] Figure 3 This is a diagram of the growth pattern of carbonate platform in an embodiment of the present invention;

[0025] Figure 4 A diagram showing the growth process of a carbonate platform restored in an embodiment of the present invention;

[0026] Figure 5 Schematic diagram of a quantitative analysis device for carbonate platform evolution according to an embodiment of the present invention;

[0027] Figure 6 This is another schematic diagram of a quantitative analysis device for carbonate platform evolution according to an embodiment of the present invention. DETAILED DESCRIPTION

[0028] To make the purpose, technical solutions and advantages of the embodiments of the present invention more clear, the embodiments of the present invention are further described in detail below with reference to 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.

[0029] In order to overcome the shortcomings of the existing technology, the quantitative analysis method of carbonate platform evolution proposed in the embodiment of the present invention mainly improves the analysis accuracy of carbonate platform evolution by constructing a mathematical model for calculating the vertical growth rate of the platform. Figure 1 This is a flow chart of a quantitative analysis method for carbonate platform evolution in an embodiment of the present invention.

[0030] like Figure 1 As shown, the method includes the following steps:

[0031] Step 101: Identify sedimentary cycle units by analyzing outcrop data and well logging data, establish a sedimentary sequence, and perform absolute time calibration on the sedimentary sequence;

[0032] Step 102: Identify sea level rise and fall cycles in the sedimentary sequence, analyze the sea level rise and fall cycles, obtain relative sea level changes at each absolute time calibration point, and calculate the relative sea level change rate;

[0033] Step 103: Obtain the tectonic settlement amount at each absolute time calibration point by backstripping method, and calculate the tectonic settlement rate at each absolute time calibration point;

[0034] Step 104: Calculate the carbonate rock yield at each absolute time calibration point using the relative sea level change and tectonic subsidence at each absolute time calibration point;

[0035] Step 105: Calculate the sum of the tectonic subsidence rate and the relative sea level change rate at each absolute time calibration point to obtain the accommodation space consumption rate at each absolute time calibration point; use the difference between the carbonate rock yield and the accommodation space consumption rate at each absolute time calibration point to quantitatively analyze the evolution of the carbonate platform.

[0036] Compared with the technical solutions for qualitatively describing the evolution of carbonate platforms in the existing technology, this method establishes a sedimentary sequence by analyzing outcrop data and well logging data, and performs absolute time calibration on the sedimentary sequence; takes absolute time as the benchmark, integrates multi-source data, and improves the analysis accuracy of carbonate platform evolution; obtains relative sea level change through sedimentary cycle analysis, and calculates the relative sea level change rate; inverts the tectonic subsidence history through the backstripping method, obtains the tectonic subsidence, and calculates the tectonic subsidence rate; calculates the carbonate rock yield at each absolute time calibration point using the relative sea level change and the tectonic subsidence; calculates the sum of the tectonic subsidence rate and the sea level change rate at each absolute time calibration point, and obtains the accommodative space consumption rate at each absolute time calibration point; uses the difference between the carbonate rock yield and the accommodative space consumption rate at each absolute time calibration point as a mathematical model for analyzing carbonate platform evolution, which can realize quantitative analysis of carbonate platform evolution, provide a basis for improving reservoir prediction accuracy, and thus reduce exploration costs and shorten exploration cycles.

[0037] When constructing mathematical models to analyze the evolution of carbonate platforms, time dimension information is essential. To enable the mathematical models to more accurately reflect geological processes, a continuous time scale was established for the sedimentary sequence.

[0038] In the embodiment of the present invention, outcrop data and well logging data are analyzed to identify sedimentary cycle units, establish a sedimentary sequence, and perform absolute time calibration on the sedimentary sequence.

[0039] Among them, outcrop data mainly include: field observation of lithology, sedimentary structure (such as bedding, reef structure), and stratigraphic contact relationship (such as unconformity surface); well logging data: through gamma, resistivity, acoustic wave and other curves, identify lithology changes and cycle boundaries (such as low natural gamma values ​​may indicate pure carbonate rocks).

[0040] In one embodiment, sedimentary cycle units are identified by analyzing outcrop data and well logging data; the sedimentary cycle units are arranged in chronological order to obtain a sedimentary sequence; the characteristics of the well logging data and the continuity of the seismic reflection layer are used to divide the cyclostratigraphy in the sedimentary sequence; astronomical information is extracted from the well logging data through spectral analysis; the absolute age data of the cyclostratigraphy is obtained using isotope chronology and biostratigraphy; the absolute age data of the cyclostratigraphy is correlated with the astronomical information to calibrate the absolute time of the cyclostratigraphy.

[0041] Specifically, cyclostratigraphic delineation requires combining the morphological and numerical characteristics of well logging data. Seismic data contains information about subsurface structures (such as faults, uplifts, and stratigraphic boundary reflections), covering regional-scale tectonic patterns. Cyclostratigraphic boundaries are primarily indicated by reflection structures such as onlap and truncation, which demonstrate the continuity of seismic reflectors.

[0042] On the basis of the division of cyclostratigraphy, cyclostratigraphy analysis is further carried out, mainly by performing spectral analysis on the logging data, identifying the astronomical information in the logging data, and combining isotope chronology and biostratigraphy to constrain the astronomical information, thereby realizing the calibration of the absolute time of the cyclostratigraphy.

[0043] For example, the absolute time of the cyclostratigraphic formation can be calculated according to formula (1):

[0044]

[0045] Where T is the total time span of the cyclostratigraphic formation, H i is the thickness of the ith sedimentary cycle unit in the cyclolith (in one embodiment, the average thickness of the sedimentary cycle unit calculated from the measured outcrop data is 3 m), R i is the sedimentation rate of the i-th sedimentary cycle unit in the cyclostratigraphy (the sedimentation rate of limestone is 2 cm / kyr, and the sedimentation rate of marl is 1 cm / kyr), and n is the number of sedimentary cycle units in the cyclostratigraphy.

[0046] It should be noted that in the embodiments of the present invention, data of different resolutions (such as outcrop data, seismic data, logging data, analytical and testing data, etc.) are integrated to resolve differences in temporal and spatial scales (such as outcrop high resolution but local, seismic low resolution but regional coverage).

[0047] The evolution of carbonate platforms is governed by both sedimentary and tectonic factors. Cyclic characteristics reflect the periodic changes in the sedimentary environment, such as sea level rise and fall and climate change, while tectonic characteristics demonstrate the influence of crustal movements on the platform's morphology and stability. Analysis of cyclical characteristics can identify regular changes in the sedimentary sequence and reconstruct paleoenvironmental and sedimentary history. Analysis of tectonic characteristics reveals the tectonic setting and evolutionary dynamics of the platform. By combining these two characteristics, we can more accurately reconstruct the evolution of carbonate platforms, providing a scientific basis for related geological research and resource exploration.

[0048] In an embodiment of the present invention, sea level rise and fall cycles are identified in a sedimentary sequence, analyzed, relative sea level changes at each absolute time calibration are obtained, and the relative sea level change rate is calculated.

[0049] In the embodiment of the present invention, the tectonic settlement amount at each absolute time calibration point is obtained by the backstripping method, and the tectonic settlement rate at each absolute time calibration point is calculated.

[0050] By using seismic data to identify fault characteristics (fault vertical throw, fault activity time span) and uplift characteristics (paleo-uplift location), tectonic activities can be divided into different types based on the fault characteristics and uplift characteristics of carbonate platforms, and different types of tectonic activities can be analyzed using different tectonic parameters.

[0051] For example, in a carbonate platform region, the embodiment of the present invention uses the tectonic subsidence rate as a tectonic parameter for quantitatively analyzing tectonic deformation activities (such as basin subsidence, frontal depression, and subsidence of superimposed basins) and calculates it according to formula (2):

[0052]

[0053] where S(t) is the tectonic subsidence at the absolute time calibration point t (unit: m), H(t) is the formation thickness at the absolute time calibration point t (obtained by inversion of well logging data / seismic data), φ(t) is the porosity decay function with depth, and H0 is the initial sedimentary thickness.

[0054] It is also possible to reflect the activity intensity of the fault and the characteristics of local structural deformation by calculating the fault activity rate based on the local characteristics of the carbonate platform. Therefore, in the process of tectonic activity analysis, the vertical displacement rate of the fault can also be used as another structural parameter for quantitative analysis of structural deformation activities (such as strike-slip faults, normal faults, and negative faults). The calculation method is shown in formula (3):

[0055]

[0056] Where Vf is the vertical displacement rate of the fault, ΔD is the vertical throw of the fault identified by the seismic profile, and Δt is the time span of fault activity.

[0057] Figure 2 This is an example diagram of a seismic profile across the Jiangliangping Trough in an embodiment of the present invention. Figure 2 As shown in the figure, the yellow dotted line at the position of well Longg2 indicates the position of the platform edge, the red arrow at the position of well longg10 indicates the direction of platform migration, the cyan dotted line marks the seismic reflection line of the fourth section of the Triassic Feixianguan Formation during the T1f4 period, the red dotted line marks the seismic reflection line of the first section of the Triassic Feixianguan Formation during the T1f1 period, the green dotted line marks the seismic reflection line of the Permian Changxing Formation during the P2ch period, and the blue dotted line marks the seismic reflection line of the Permian Liangshan Formation during the P2I period. Faults can cause seismic reflection lines to be interrupted or dislocated at the fault, so that the vertical fault throw of the fault can be identified by the seismic profile.

[0058] Before constructing a mathematical model for quantitative analysis of carbonate platform evolution, the spatiotemporal correlation (such as Pearson coefficient, cross-spectral analysis, etc.) of statistical cycle parameters (such as the thickness of sedimentary cycle units and the frequency of sedimentary cycle units) and structural parameters (such as the vertical displacement rate of faults and the structural subsidence rate) is first performed to further improve the accuracy of the quantitative analysis of carbonate platform evolution in the embodiment of the present invention.

[0059] Taking the Pearson coefficient as an example, the correlation coefficient calculation is shown in formula (4):

[0060]

[0061] Where R is the correlation coefficient (range [-1, 1], measuring the correlation between tectonic activity and cycle parameters), x i is the ith structural parameter (such as structural subsidence rate), y i is the i-th cycle parameter (such as thickness, frequency).

[0062] Through spatiotemporal correlation analysis, the tectonic subsidence rate is selected in the embodiment of the present invention to perform quantitative analysis on the evolution of carbonate platforms.

[0063] In the embodiment of the present invention, the relative sea level change and the tectonic subsidence at each absolute time calibration point are respectively used to calculate the carbonate rock yield at each absolute time calibration point.

[0064] In one embodiment, the actual water depth at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point; the carbonate rock yield at each water depth is calculated to obtain the carbonate rock yield at each absolute time calibration point; wherein the distribution of the carbonate rock salt yield in the water depth direction obeys a Gaussian distribution.

[0065] For example, the carbonate rock yield at each water depth can be calculated according to formula (5):

[0066]

[0067] In the above formula, P(z) is the carbonate rock yield at water depth z, P max is the optimal yield of carbonate rock, z0 is the water depth of the optimal yield of carbonate rock, k is the attenuation coefficient, and e is a natural constant.

[0068] In one embodiment, the actual water depth at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point according to formula (6):

[0069] D(t)=D0+S(t)-ΔL(t) (6)

[0070] In the above formula, D(t) is the actual water depth at the absolute time calibration point t, D0 is the initial water depth, S(t) is the tectonic subsidence at the absolute time calibration point t, and ΔL(t) is the relative sea level change at the absolute time calibration point t.

[0071] The vertical growth rate of a carbonate platform is a key indicator of its vertical development and changes, reflecting the combined impact of sediment accumulation and tectonic activity on the platform. This embodiment of the present invention calculates the vertical growth rate of carbonate platforms by analyzing the relationship between cyclic and tectonic characteristics and carbonate yield, thereby quantitatively analyzing carbonate platform evolution.

[0072] In an embodiment of the present invention, the sum of the tectonic subsidence rate and the relative sea level change rate at each absolute time calibration is calculated to obtain the accommodation space consumption rate at each absolute time calibration; the difference between the carbonate rock yield and the accommodation space consumption rate at each absolute time calibration is used to quantitatively analyze the evolution of the carbonate platform.

[0073] For example, the mathematical model for quantitative analysis of carbonate platform evolution is shown in formula (7):

[0074]

[0075] In the above formula, x is the platform position, t is the absolute time calibration point, h(x,t) is the vertical thickness of the platform at position x and absolute time calibration point t (unit: m), Pcarb(x,t) is the carbonate rock yield (unit: m / kyr, related to light, water temperature, and nutrients), Stect(x,t) is the tectonic subsidence (unit: m, calculated by the backstripping method), and SL(t) is the relative sea level change (unit: m).

[0076] In addition to the vertical growth rate of the platform, the lateral migration rate of the platform margin is also an important indicator in the analysis of carbonate platform evolution. It can provide key information about platform expansion or contraction, sedimentary environment changes, platform stability, paleogeographic reconstruction, and reservoir distribution, and analyze the migration of the platform margin.

[0077] In one embodiment, the method further includes: calculating the water depth change using the relative sea level change and the tectonic subsidence; interpreting the seismic profile of the seismic data to obtain the platform front slope; obtaining the platform edge migration distance based on the water depth change and the platform front slope, and calculating the lateral migration rate of the platform edge.

[0078] For example, the lateral migration rate of the platform edge can be calculated according to formula (8), which is as follows:

[0079]

[0080] In the above formula, V m is the lateral migration rate of the platform edge, Δt is the time period, ΔD is the change in water depth within Δt, and θ is the slope of the platform front (obtained from the measured seismic profile).

[0081] Figure 3This is a diagram of the growth pattern of carbonate platform in an embodiment of the present invention. Figure 3 In the figure, TST is the transgressive system tract, HST is the highstand system tract, and FSST & LST are the forced regressive system tract & lowstand system tract.

[0082] In one embodiment, when the difference between the carbonate rock yield and the accommodation space consumption rate is greater than zero, the terrace at the absolute time calibration grows upward; when the difference between the carbonate rock yield and the accommodation space consumption rate is equal to zero, the terrace at the absolute time calibration stops growing; when the difference between the carbonate rock yield and the accommodation space consumption rate is less than zero, the terrace at the absolute time calibration is submerged.

[0083] For example, if the platform edge Pcarb = 0.5 m / kyr, the tectonic subsidence rate Sea level rise rate but This indicates that the platform is vertically accreting at a rate of 0.2 mm per year.

[0084] Figure 4 Figure 2 is a diagram of the growth process of the carbonate platform restored in the embodiment of the present invention. Figure 4 As shown, the growth process of the carbonate platform from profile position A to profile position B (see red marks in the figure) is shown: from the deposition of the late Permian P2, the end of the late Permian P2, the end of the late Triassic T3, the end of the Jurassic J to the present period. Figure 4 The middle geological ages are as follows: J1z represents the Lower Jurassic, T3x1 represents the Upper Triassic, T2l represents the Middle Triassic, T1j1 represents the Lower Triassic, T1f1 represents the Lower Triassic, F2l represents the Middle Devonian, P1l represents the Lower Permian, O3l represents the Upper Ordovician, O1n represents the Lower Ordovician, ∈1l represents the Lower Cambrian, ∈1q represents the Lower Cambrian, Z2n represents the Upper Sinian, and Z1n represents the Lower Sinian.

[0085] In this embodiment, by integrating multi-source data, coupled cyclical and tectonic features, the relative control of tectonic subsidence and sea-level change on platform evolution is clarified. This significantly improves the accuracy of quantitative analysis, overcoming the inability to quantify rates and pathways in traditional seismic facies research on carbonate facies migration. It also provides a basis for influencing oil and gas exploration by improving the accuracy of later reservoir predictions.

[0086] The present invention also provides a device for quantitatively analyzing carbonate platform evolution, as described in the following embodiments. Because the principles underlying the device are similar to those of the quantitative analysis method for carbonate platform evolution, the implementation of the device can be referenced to the implementation of the quantitative analysis method for carbonate platform evolution, and any repetitions will not be repeated.

[0087] Figure 5 Schematic diagram of a quantitative analysis device for carbonate platform evolution according to an embodiment of the present invention. Figure 5 As shown, the device includes:

[0088] Sedimentary cycle analysis module 501 is used to: identify sedimentary cycle units by analyzing outcrop data and well logging data, establish sedimentary sequences, and perform absolute time calibration on the sedimentary sequences;

[0089] The relative sea level analysis module 502 is used to identify sea level rise and fall cycles in the sedimentary sequence, analyze the sea level rise and fall cycles, obtain the relative sea level change at each absolute time calibration point, and calculate the relative sea level change rate;

[0090] The structural settlement analysis module 503 is used to obtain the structural settlement amount at each absolute time calibration point by backstripping method and calculate the structural settlement rate at each absolute time calibration point;

[0091] The carbonate rock yield analysis module 504 is used to calculate the carbonate rock yield at each absolute time calibration point using the relative sea level change and tectonic subsidence at each absolute time calibration point;

[0092] The carbonate platform evolution analysis module 505 is used to calculate the sum of the tectonic subsidence rate and the relative sea level change rate at each absolute time calibration to obtain the accommodation space consumption rate at each absolute time calibration; and use the difference between the carbonate rock yield and the accommodation space consumption rate at each absolute time calibration to quantitatively analyze the evolution of the carbonate platform.

[0093] In one embodiment, the sedimentation cycle analysis module 501 is specifically configured to:

[0094] Identify sedimentary cycle units by analyzing outcrop data and well logging data; arrange sedimentary cycle units in chronological order to obtain a sedimentary sequence; and divide cyclostratigraphy within the sedimentary sequence using the characteristics of well logging data and the continuity of seismic reflection layers.

[0095] Astronomical information is extracted from well logging data through spectral analysis; absolute age data of cyclostratigraphy are obtained using isotope chronology and biostratigraphy; and the absolute age data of cyclostratigraphy are correlated with astronomical information to calibrate the absolute time of cyclostratigraphy.

[0096] In one embodiment, the carbonate rock yield analysis module 503 is specifically configured to:

[0097] The actual water depth at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point.

[0098] The carbonate rock yield at each water depth is calculated to obtain the carbonate rock yield at each absolute time calibration point; among which, the distribution of carbonate rock salt yield in the water depth direction obeys Gaussian distribution.

[0099] In one embodiment, the carbonate rock yield analysis module 503 is specifically configured to:

[0100] The actual water depth at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point according to the following formula:

[0101] D(t)=D0+S(t)-ΔL(t)

[0102] where D(t) is the actual water depth at the absolute time calibration point t, D0 is the initial water depth, S(t) is the tectonic subsidence at the absolute time calibration point t, and ΔL(t) is the relative sea level change at the absolute time calibration point t.

[0103] Figure 6 FIG. 1 is another schematic diagram of a quantitative analysis device for carbonate platform evolution according to an embodiment of the present invention. Figure 6 As shown, the device further includes a terrace edge migration analysis module 601, which is used to:

[0104] Calculate the change in water depth using the relative sea level change and tectonic subsidence;

[0105] The seismic data is interpreted as a seismic profile to obtain the platform front slope; based on the water depth change and the platform front slope, the platform edge migration distance is obtained and the lateral migration rate of the platform edge is calculated.

[0106] In one embodiment, the carbonate platform evolution analysis module 505 is specifically configured to:

[0107] When the difference between the carbonate rock yield and the accommodation space consumption rate is greater than zero, the value of the upward growth of the platform at the absolute time calibration is calculated using the difference between the carbonate rock yield and the accommodation space consumption rate at the absolute time calibration;

[0108] When the difference between the carbonate production rate and the rate of consumption of the available space is equal to zero, the platform at the absolute time calibration stops growing;

[0109] When the difference between the carbonate rock yield and the accommodation space consumption rate is less than zero, the difference between the carbonate rock yield and the accommodation space consumption rate at the absolute time calibration is used to calculate the value of the platform being flooded at the absolute time calibration.

[0110] An embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned quantitative analysis method for carbonate platform evolution when executing the computer program.

[0111] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program, which implements the above-mentioned quantitative analysis method for carbonate platform evolution when executed by a processor.

[0112] 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, it implements the above-mentioned quantitative analysis method of carbonate platform evolution.

[0113] Compared with the technical solutions in the prior art for qualitatively describing the evolution of carbonate platforms, the embodiments of the present invention integrate multi-source data, perform spatiotemporal correlation analysis on structural parameters and cycle parameters, select strongly correlated cycle features and structural features for coupling, clarify the relative control effects of structural subsidence and sea level changes on platform evolution, and construct a mathematical model for analyzing the evolution of carbonate platforms, thereby significantly improving the accuracy of quantitative analysis, providing a basis for improving the accuracy of reservoir prediction, and thus reducing exploration costs and shortening exploration cycles.

[0114] It will be understood by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Thus, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Furthermore, 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 magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0115] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, 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 generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0116] 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 an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0117] These computer program instructions can 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, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0118] 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

1. A quantitative analysis method for carbonate platform evolution, characterized in that: include: By analyzing outcrop data and well logging data, the sedimentary cycle units are identified, the sedimentary sequence is established, and the absolute time calibration of the sedimentary sequence is performed; Identify and analyze sea-level cycles in sedimentary sequences, obtain relative sea-level changes at each absolute time point, and calculate relative sea-level change rates; The tectonic settlement amount at each absolute time calibration point is obtained by backstripping method, and the tectonic settlement rate at each absolute time calibration point is calculated; The carbonate rock yield at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point. The sum of the tectonic subsidence rate and the relative sea level change rate at each absolute time calibration point was calculated to obtain the accommodation space consumption rate at each absolute time calibration point. The difference between the carbonate rock yield and the accommodation space consumption rate at each absolute time calibration point was used to quantitatively analyze the evolution of the carbonate platform.

2. The method according to claim 1, wherein By analyzing outcrop data and well logging data, the sedimentary cycle units are identified, the sedimentary sequence is established, and the absolute time calibration of the sedimentary sequence is performed, including: Identify sedimentary cycle units by analyzing outcrop data and well logging data; arrange sedimentary cycle units in chronological order to obtain a sedimentary sequence; and divide cyclostratigraphy within the sedimentary sequence using the characteristics of well logging data and the continuity of seismic reflection layers. Astronomical information is extracted from well logging data through spectral analysis; absolute age data of cyclostratigraphy are obtained using isotope chronology and biostratigraphy; and the absolute age data of cyclostratigraphy are correlated with astronomical information to calibrate the absolute time of cyclostratigraphy.

3. The method according to claim 1, wherein The relative sea level change and tectonic subsidence at each absolute time calibration point are used to calculate the carbonate rock yield at each absolute time calibration point, including: The actual water depth at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point. The carbonate rock yield at each water depth is calculated to obtain the carbonate rock yield at each absolute time calibration point; among which, the distribution of carbonate rock salt yield in the water depth direction obeys Gaussian distribution.

4. The method according to claim 3, wherein The actual water depth at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point according to the following formula: D(t)=D0+S(t)-ΔL(t) where D(t) is the actual water depth at the absolute time calibration point t, D0 is the initial water depth, S(t) is the tectonic subsidence at the absolute time calibration point t, and ΔL(t) is the relative sea level change at the absolute time calibration point t.

5. The method according to claim 1, wherein Also includes: Calculate the change in water depth using the relative sea level change and tectonic subsidence; Seismic data were interpreted as seismic profiles to obtain the platform front slope; Based on the change in water depth and the slope of the platform front, the platform edge migration distance is obtained and the lateral migration rate of the platform edge is calculated.

6. The method according to claim 1, wherein The difference between the carbonate production rate and the accommodation space consumption rate at each absolute time calibration point is used to quantitatively analyze the evolution of carbonate platforms, including: When the difference between the carbonate rock yield and the accommodation space consumption rate is greater than zero, the value of the upward growth of the platform at the absolute time calibration is calculated using the difference between the carbonate rock yield and the accommodation space consumption rate at the absolute time calibration; When the difference between the carbonate production rate and the rate of consumption of the available space is equal to zero, the platform at the absolute time calibration stops growing; When the difference between the carbonate rock yield and the accommodation space consumption rate is less than zero, the difference between the carbonate rock yield and the accommodation space consumption rate at the absolute time calibration is used to calculate the value of the platform being flooded at the absolute time calibration.

7. A quantitative analysis device for carbonate platform evolution, characterized in that: include: Sedimentary cycle analysis module is used to: identify sedimentary cycle units, establish sedimentary sequences, and perform absolute time calibration on sedimentary sequences by analyzing outcrop data and well logging data; The relative sea level analysis module is used to identify and analyze sea level cycles in sedimentary sequences, obtain relative sea level changes at each absolute time point, and calculate relative sea level change rates. The structural settlement analysis module is used to obtain the structural settlement amount at each absolute time calibration point through the backstripping method and calculate the structural settlement rate at each absolute time calibration point; The carbonate rock yield analysis module is used to calculate the carbonate rock yield at each absolute time calibration point using the relative sea level change and tectonic subsidence at each absolute time calibration point; The carbonate platform evolution analysis module is used to calculate the sum of the tectonic subsidence rate and the relative sea level change rate at each absolute time calibration to obtain the accommodation space consumption rate at each absolute time calibration; and use the difference between the carbonate rock yield and the accommodation space consumption rate at each absolute time calibration to quantitatively analyze the evolution of the carbonate platform.

8. The device according to claim 7, wherein Sedimentary cycle analysis module, specifically used for: Identify sedimentary cycle units by analyzing outcrop data and well logging data; arrange sedimentary cycle units in chronological order to obtain a sedimentary sequence; and divide cyclostratigraphy within the sedimentary sequence using the characteristics of well logging data and the continuity of seismic reflection layers. Astronomical information is extracted from well logging data through spectral analysis; absolute age data of cyclostratigraphy are obtained using isotope chronology and biostratigraphy; and the absolute age data of cyclostratigraphy are correlated with astronomical information to calibrate the absolute time of cyclostratigraphy.

9. The device according to claim 7, wherein Carbonate rock yield analysis module, specifically used for: The actual water depth at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point. The carbonate rock yield at each water depth is calculated to obtain the carbonate rock yield at each absolute time calibration point; among which, the distribution of carbonate rock salt yield in the water depth direction obeys Gaussian distribution.

10. The device according to claim 9, wherein Carbonate rock yield analysis module, specifically used for: The actual water depth at each absolute time calibration point is calculated using the relative sea level change and tectonic subsidence at each absolute time calibration point according to the following formula: D(t)=D0+S(t)-ΔL(t) where D(t) is the actual water depth at the absolute time calibration point t, D0 is the initial water depth, S(t) is the tectonic subsidence at the absolute time calibration point t, and ΔL(t) is the relative sea level change at the absolute time calibration point t.

11. The device according to claim 7, wherein Also includes a platform edge migration analysis module for: Calculate the change in water depth using the relative sea level change and tectonic subsidence; The seismic data is interpreted as a seismic profile to obtain the platform front slope; based on the water depth change and the platform front slope, the platform edge migration distance is obtained and the lateral migration rate of the platform edge is calculated.

12. The device according to claim 7, wherein Carbonate platform evolution analysis module is specifically used for: When the difference between the carbonate rock yield and the accommodation space consumption rate is greater than zero, the value of the upward growth of the platform at the absolute time calibration is calculated using the difference between the carbonate rock yield and the accommodation space consumption rate at the absolute time calibration; When the difference between the carbonate production rate and the rate of consumption of the available space is equal to zero, the platform at the absolute time calibration stops growing; When the difference between the carbonate rock yield and the accommodation space consumption rate is less than zero, the difference between the carbonate rock yield and the accommodation space consumption rate at the absolute time calibration is used to calculate the value of the platform being flooded at the absolute time calibration.

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

14. 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 6 is implemented.

15. 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 6 is implemented.

Images