Method and device suitable for restoring prototype basin in basin depression-augmentation stage

Abstract

The invention discloses a method and a device suitable for restoring a prototype basin in a basin depression-augmentation stage, and relates to the technical field of petroleum geology, and the method suitable for restoring the prototype basin in the basin depression-augmentation stage mainly comprises the following steps: carrying out basin boundary identification and unconformity surface analysis to determine a denudation period, a denudation point and a stratum extension trend; generating a residual thickness map; denudation amount recovery is carried out, denudation interpretation is carried out on the seismic section, and the original deposition thickness of the stratum is formed; and performing compaction correction and paleo-water depth correction to obtain the accurate original thickness of the stratum. By implementing the method and the device suitable for restoring the prototype basin in the basin depression-augmentation stage, provided by the invention, the prototype basin can be accurately restored.

Description

Technical Field

[0001] The present invention relates to the field of petroleum geology technology, and more particularly to a method and device suitable for restoring a prototype basin in a basin depression-uplift stage. Background Art

[0002] In the early days, the restoration of basin prototypes was relatively simple. Generally, the sedimentary thickness of the basin was restored vertically by the denudation thickness method. After more than 30 years of development, Tong Xiaoguang and He Dengfa (2001) proposed that the reasonable restoration of basin prototypes should include the following six aspects: (1) analysis of the tectonic background; (2) research on paleoflow direction, lithofacies, sedimentary system and filling process; (3) analysis of paleoecology, paleoclimate and paleogeographic environment; (4) analysis of tectonic subsidence; (5) analysis of basin fault system; (6) analysis of volcanic rocks and their tectonic background. However, these methods are limited to the content that needs to be restored and basically do not involve specific methods and feasibility of actual operation. Moreover, the prototype restoration methods of basins of different properties are also different, and the focus is also different. For example, the key to the restoration of foreland basins is the restoration of the thrust belt in the orogenic wedge, and the key to the restoration of strike-slip related basins is the determination of strike-slip displacement at different stages and its relationship with basin filling. Therefore, it is urgent to establish a set of scientific and feasible basin restoration methods and processes based on the characteristics of different types of basins.

[0003] The above content is only used to assist in understanding the technical solution of the present invention and does not constitute an admission that the above content is prior art. Summary of the Invention

[0004] The purpose of the present invention is to provide a method and device suitable for restoring a prototype basin in the depression-uplift stage of a basin, which can accurately restore the prototype basin.

[0005] The present invention provides a method for restoring a prototype basin in the depression-uplift stage of a basin, comprising the following steps:

[0006] S1: Based on the basin boundary identification map, perform basin boundary identification and unconformity surface analysis to determine the erosion period, erosion point, and stratigraphic extension trend;

[0007] S2: generating a residual thickness map based on the erosion period, erosion point, stratigraphic extension trend and seismic logging data;

[0008] S3: Based on the residual thickness map, perform denudation recovery and perform denudation interpretation on the seismic profile to obtain the original sedimentary thickness of the stratum;

[0009] S4: Based on the original sedimentary thickness of the stratum, compaction correction and paleo-water depth correction are performed to obtain the accurate original thickness of the stratum.

[0010] The present invention also provides a device for restoring a prototype basin in the depression-uplift stage of a basin. The device for restoring a prototype basin in the depression-uplift stage of a basin comprises the following modules:

[0011] The basin boundary identification module is configured to: perform basin boundary identification and unconformity surface analysis based on the basin boundary identification map, and determine the erosion period, erosion point, and stratigraphic extension trend;

[0012] A residual thickness map generating module is configured to generate a residual thickness map based on the erosion period, erosion point, formation extension trend and seismic logging data;

[0013] an erosion recovery module configured to: recover the erosion amount based on the residual thickness map, perform erosion interpretation on the seismic profile, and form the original sedimentary thickness of the stratum;

[0014] The correction module is configured to perform compaction correction and paleo-water depth correction according to the original sedimentary thickness of the stratum to obtain the accurate original thickness of the stratum.

[0015] The present invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-mentioned method for restoring a prototype basin in the depression-uplift stage of a basin.

[0016] The present invention also 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 program, the steps of the above-mentioned method for restoring a prototype basin in the depression-uplift stage of the basin are implemented.

[0017] The implementation of the method and device provided by the present invention for restoring a prototype basin in the depression-uplift stage has the following beneficial effects:

[0018] The present invention conducts a preliminary analysis of the prototype basin type through the basin boundary type. The fault-controlled boundary needs to be restored on the fault hanging wall, and the convex-controlled boundary needs to identify the erosion point and restore it according to the stratum thickness; identify and analyze the unconformity, determine the erosion period, erosion point and stratum extension trend; restore the erosion amount and perform compaction correction, and verify the data through the seismic stratigraphic trend method and the single well sedimentation rate method; perform compaction correction and paleowater depth correction; finally, superimpose the residual paleo-geomorphology to restore the paleo-geomorphology of the prototype basin.

[0019] The present invention starts with the seismic stratigraphic trend method and the single-well sedimentation rate method in the restoration of the prototype basin, selects the strata applicable to the corresponding methods for data calculation and mapping, and then compares and analyzes the two sets of data obtained in combination with previous research data. It is found that the results obtained by the stratigraphic trend method and the sedimentation rate method are consistent. Ultimately, the paleogeomorphological characteristics of the prototype basin can be restored more accurately, so as to more objectively evaluate the oil and gas resources in the area.

[0020] This method can use the stratigraphic trend method and the sedimentation rate method to mutually verify and constrain each other, reduce data errors, and more accurately restore the prototype basin. For the Junggar Basin studied by this method, the Permian-Jurassic strata in the Junggar Basin were characterized by frequent tectonic activity during the deposition of the Permian strata, resulting in large variations in sedimentation rates. The seismic trend method was the primary method used. The Triassic Basin was in a depression phase, with relatively stable and continuous sedimentation rates overall. Therefore, the Triassic-Jurassic strata meet the applicable conditions for the sedimentation rate method. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] The present invention will be further described below with reference to the accompanying drawings and embodiments, in which:

[0022] Figure 1 This is a flow chart of a method for restoring a prototype basin in a basin depression-uplift stage provided by the present invention;

[0023] Figure 2 This is a framework diagram for implementing the method for restoring a prototype basin in the depression-uplift stage provided by the present invention;

[0024] Figure 3 This is the basin boundary identification map in Example 1;

[0025] Figure 4 For the identification and analysis of the unconformity surface in Example 1;

[0026] Figure 5 is the residual thickness diagram in Example 1;

[0027] Figure 6 This is the P1 residual thickness map and the corresponding seismic section interpretation;

[0028] Figure 7 This is the result map after the P1 residual thickness map is restored and the corresponding seismic profile interpretation;

[0029] Figure 8 This is the restoration result diagram of the prototype basin in Example 1;

[0030] Figure 9 It is a structural block diagram of the computer device provided by the present invention. DETAILED DESCRIPTION

[0031] In order to have a clearer understanding of the technical features, purposes and effects of the present invention, specific embodiments of the present invention are now described in detail with reference to the accompanying drawings.

[0032] Figure 1 A schematic diagram of a method for restoring a prototype basin in the depression-uplift stage of a basin according to this embodiment is shown. In this embodiment, the method for restoring a prototype basin in the depression-uplift stage of a basin comprises the following steps:

[0033] S1: Based on the basin boundary identification map, perform basin boundary identification and unconformity surface analysis to determine the erosion period, erosion point and stratigraphic extension trend;

[0034] As an exemplary embodiment, in step S1, basin boundaries are identified; and a preliminary analysis of the prototype basin type is performed based on the basin boundary type.

[0035] S2: Generate a residual thickness map using seismic interpretation methods and geological tools based on the erosion period, erosion point, stratigraphic extension trend, and seismic logging data;

[0036] As an exemplary embodiment, in step S2, a residual thickness map is prepared; layer data (grid) is formed through seismic interpretation, and a residual thickness map is formed with the help of petrel software.

[0037] S3: Based on the residual thickness map, perform denudation recovery and perform denudation interpretation on the seismic profile to obtain the original sedimentary thickness of the stratum;

[0038] In an exemplary embodiment, the erosion amount recovery is performed as follows:

[0039]

[0040] Where K is the deposition rate coefficient, which is expressed as the ratio of the deposition rate of the recovery layer to the deposition rate of the reference layer. Usually, wells with stable deposition in the depression area and no erosion are selected for calculation. As the deposition rate accelerates or slows down, the K value can be greater than 1 or less than 1; ΔH A is the erosion thickness of the unconformity surface; H A 、T A are the sedimentary thickness and sedimentary time of the unconformity surface, respectively; H B is the sedimentary thickness of the stratum.

[0041] As an exemplary embodiment, in step S3, erosion recovery is performed, and erosion interpretation is performed on the seismic profile to form a new stratigraphic interface, i.e., the restored interface, which is mapped using petrel. Erosion recovery and compaction correction are performed, specifically using the stratigraphic trend method and the sedimentation rate method to recover the erosion, and the two methods are verified.

[0042] The stratigraphic trend method uses seismic data to comprehensively compare the stratigraphic thickness of uneroded areas, fill in the eroded areas, and restore the original sedimentary thickness of the eroded strata based on the regional sedimentary-tectonic background. The sedimentation rate method uses the sedimentation rate and its changing trend within the same structural layer to restore the original sedimentary thickness of the strata.

[0043] The specific calculation formula for the recovery of erosion using the sedimentation rate method is:

[0044]

[0045] Where K is the deposition rate coefficient, which is expressed as the ratio of the deposition rate of the recovery layer to the deposition rate of the reference layer (V0 / V1). Usually, wells with stable deposition in the depression area and no erosion are selected for calculation. As the deposition rate accelerates or slows down, the K value can be greater than 1 or less than 1; ΔH A is the erosion thickness of the unconformity surface, m; H A 、T A is the sedimentary thickness and sedimentary time of the unconformity surface, m, Ma; H B 、T B is the sedimentary thickness and sedimentary time of the stratum, m, Ma.

[0046] S4: Based on the original sedimentary thickness of the stratum, compaction correction and paleo-water depth correction are performed to obtain the accurate original thickness of the stratum;

[0047] In an exemplary embodiment, step S4 specifically includes: performing compaction correction and paleo-water depth correction according to the original sedimentary thickness of the stratum to obtain the accurate original thickness of the stratum, such as the formula:

[0048]

[0049] in, is the porosity of clastic rocks, is the surface porosity, z is the burial depth of the rock layer, c is the compaction factor, Z1 and Z2 are the top and bottom burial depths of the stratum after compaction, respectively. ′ 1 and Z ′ 2 is the top surface burial depth and bottom surface burial depth;

[0050] In another exemplary embodiment, compaction correction and paleo-water depth correction are performed based on the original sedimentary thickness of the stratum to obtain the accurate original thickness of the stratum, as shown in the formula:

[0051]

[0052] in, is the porosity of clastic rocks, is the surface porosity, z is the burial depth of the rock layer, k is the compaction factor, Z1 and Z2 are the top and bottom burial depths of the stratum after compaction, respectively. ′ 1 and Z ′ 2 is the top surface burial depth and bottom surface burial depth;

[0053] As an exemplary embodiment, in step S4, compaction correction and paleowater depth correction are performed; during the burial process, due to the pressure of overlying sediments and water bodies or tectonic stress, the debris particles in the sediments are arranged closely and the porosity is reduced, resulting in the formation thickness decreasing with increasing burial depth, so it is necessary to restore the compaction volume of the target layer; when the basin subsidence center is inconsistent with the deposition center, the thickness of the sediment cannot reflect the characteristics of the paleolandform, so it is necessary to correct the paleowater depth.

[0054] This embodiment provides an apparatus for restoring a prototype basin during the depression-uplift stage of a basin, comprising the following modules: a basin boundary identification module, configured to identify the basin boundary and analyze unconformity surfaces based on a basin boundary identification map, and determine the erosion period, erosion points, and stratigraphic extension trends; a residual thickness map generation module, configured to generate a residual thickness map based on the erosion period, erosion points, stratigraphic extension trends, and seismic logging data; an erosion amount recovery module, configured to recover the erosion amount based on the residual thickness map, perform erosion interpretation on seismic profiles, and determine the original sedimentary thickness of the stratum; and a correction module, configured to perform compaction correction and paleowater depth correction based on the original sedimentary thickness of the stratum to obtain an accurate original stratum thickness.

[0055] In an exemplary embodiment, the above-mentioned erosion amount recovery is performed as follows:

[0056]

[0057] Where K is the deposition rate coefficient, which is expressed as the ratio of the deposition rate of the recovery layer to the deposition rate of the reference layer. Usually, wells with stable deposition in the depression area and no erosion are selected for calculation. As the deposition rate accelerates or slows down, the K value can be greater than 1 or less than 1; ΔH A is the erosion thickness of the unconformity surface; H A 、T A are the sedimentary thickness and sedimentary time of the unconformity surface, respectively; H B is the sedimentary thickness of the stratum.

[0058] In an exemplary embodiment, the correction module is specifically configured to perform compaction correction and paleo-water depth correction according to the original sedimentary thickness of the stratum to obtain the accurate original thickness of the stratum, such as the formula:

[0059]

[0060] in, is the porosity of clastic rocks, is the surface porosity, z is the burial depth of the rock layer, c is the compaction factor, Z1 and Z2 are the top and bottom burial depths of the stratum after compaction, respectively. ′ 1 and Z ′ 2 is the top surface burial depth and bottom surface burial depth.

[0061] In an exemplary embodiment, the correction module is specifically configured to perform compaction correction and paleo-water depth correction according to the original sedimentary thickness of the stratum to obtain the accurate original thickness of the stratum, such as the formula:

[0062]

[0063] in, is the porosity of clastic rocks, is the surface porosity, z is the burial depth of the rock layer, k is the compaction factor, Z1 and Z2 are the top and bottom burial depths of the stratum after compaction, respectively. ′ 1 and Z ′ 2 is the top surface burial depth and bottom surface burial depth.

[0064] In one embodiment, a method for restoring a prototype basin during a depression-uplift phase of a basin comprises the following steps:

[0065] (1) Basin boundary identification: A preliminary analysis of the prototype basin type is conducted based on the basin boundary type. Fault-controlled boundaries need to be restored on the fault hanging wall. Convexity-controlled boundaries require identification of erosion points and restoration based on stratum thickness. Unconformities are identified and analyzed to determine the erosion period, erosion points, and stratigraphic extension trends.

[0066] (2) Residual thickness map production: The layer data (grid) is generated through seismic interpretation, and the residual thickness map is generated with the help of software such as Petrel.

[0067] (3) Restoration of erosion volume: Perform erosion interpretation on the seismic profile to form a new stratigraphic interface (restored interface), which is then mapped using petrel. Specifically, the stratigraphic trend method and sedimentation rate method are used to restore erosion volume. The stratigraphic trend method is based on seismic data, comprehensively compares the stratigraphic thickness of the uneroded area, fills the eroded part of the stratigraphic area, and restores the original sedimentary thickness of the eroded stratigraphic area based on the regional sedimentary-tectonic background. The sedimentation rate method is that the stratigraphic sedimentation rate within the same structural layer is relatively stable and continuous, and the original sedimentary thickness of the stratigraphic layer is restored by the sedimentation rate of the stratigraphic layer and the trend of the sedimentation rate change.

[0068] (4) Compaction correction and paleowater depth correction: During the burial process, the pressure of overlying sediments and water bodies or tectonic stresses causes the clastic particles in the sediments to be arranged tightly, resulting in a decrease in porosity and a decrease in stratum thickness with increasing burial depth. Therefore, it is necessary to restore the compaction volume of the target layer. When the basin subsidence center is inconsistent with the deposition center, the sediment thickness cannot reflect the characteristics of the paleogeomorphology, so paleowater depth correction is required.

[0069] Finally, the prototype basin restoration results are generated.

[0070] In some embodiments, the above-mentioned method for restoring a prototype basin in the basin depression-uplift stage may also be implemented in the following manner. Figure 2 This is a diagram illustrating an implementation framework of a method for restoring a prototype basin during a depression-uplift phase. In this embodiment, the method for restoring a prototype basin during a depression-uplift phase includes the following steps:

[0071] Step 101: Conduct a preliminary analysis of the prototype basin type through basin boundary types.

[0072] Reference Figure 3 , identify different boundary types from the figure. Growth faults control the thickening of the strata on the downthrown side, which is a fault-controlled boundary. The strata onlap onto the uplift is a convexity-controlled boundary. The fault-controlled boundary needs to be restored on the hanging wall of the fault. The convexity-controlled boundary requires identifying the erosion point and restoring it according to the stratum thickness. Figure 4 , identify and analyze unconformities, determine the erosion period, erosion points and stratigraphic extension trends.

[0073] Step 102: Generate layer data (grid) through seismic interpretation, and generate residual thickness map with the help of software such as Petrel. Figure 5 , combined with the collected seismic logging data, the residual thickness map is generated using software such as Petrel.

[0074] Step 103: Restoring the erosion volume: Using the formation trend method to restore the erosion volume, the specific operation steps are as follows: After mapping the restored interface using Petrel software, the new interface value can be directly read and subtracted from the residual thickness to obtain the erosion thickness; Figure 6 The residual thickness map of P1 and the seismic profile are the corresponding erosion interpretations, mainly the erosion of the Zhongguai uplift. Figure 7 This is the result map after restoration, which mainly illustrates the mapping situation of Petrel software;

[0075] The sedimentation rate method was used to restore the amount of erosion in the Triassic Basin. Since the Triassic Basin entered the depression stage, the overall sedimentation rate was relatively stable and continuous, which met the applicable conditions of the sedimentation rate method. The specific operation steps are as follows:

[0076]

[0077] Where ΔH A : erosion thickness of the K / J unconformity, m; K: sedimentation rate ratio of the Sangonghe and Badaowan formations; H A 、T A : Sedimentation thickness and deposition time of the K / J unconformity, m, Ma; H B 、T B : Sediment thickness and deposition time of the Sangonghe + Badaowan Formation, m, Ma; K value is the coefficient of the deposition rate, which is the deposition rate of the recovery layer divided by the deposition rate of the reference layer (V0 / V1). Usually, wells with stable deposition in the depression area and no erosion are selected for calculation. As the deposition rate accelerates or slows down, the K value can be greater than 1 or less than 1. In this embodiment, there are no wells in the depression area that are not eroded, so the K value is set to 1.

[0078] T B =194.5-178=16.5Ma (deposition time of the Sangonghe Formation);

[0079] Among them H B 、H A Determined by the selected single well and the recovered layer, T A Determined by the recovery horizon;

[0080] If you restore J 2-3 The erosion thickness is T A =166.1-145.5=20.6Ma;

[0081] The Badaowan Formation is located above the unconformity surface, with a relatively fast sedimentation rate and an incomplete structural stability. The Xishanyao Formation and its upper strata have been generally eroded. The Sangonghe Formation is relatively stable and less affected by erosion. Therefore, the Sangonghe Formation is selected as the reference layer.

[0082] Calculate J using the data from the comprehensive histogram of each well. 1-2 、J 2-3 The erosion thickness; J 1-2 The erosion thickness: H B =392m;T B =16.5Ma; T A =11.9Ma; H A =104m; therefore:

[0083]

[0084] J 2-3 The erosion thickness: H B =392m;T B =16.5Ma; T A =20.6Ma; H A=0m; so

[0085]

[0086] The sedimentation rate method is used to calculate the erosion amount of a single well. From points to lines, a well-connected erosion profile is formed. This profile is projected onto the seismic profile to constrain the point-line erosion amount of the formation trend method, ultimately forming the erosion zone.

[0087] Step 104: Compaction correction and paleo-water depth correction. The specific calculation steps for compaction correction are as follows: within a certain depth range, the porosity of clastic rock at any depth can be expressed as:

[0088] or

[0089] Where z is the depth of rock layer, m; c and k are compaction factors, 1 / m; is the surface porosity;

[0090] Assume that the initial porosity of the rock near the surface is As the burial depth increases, the top surface of the stratum after compaction is Z1, the bottom surface is Z2, and the porosity value becomes The thickness becomes Z1-Z2, then the thickness of the stratum skeleton within the depth range Z1-Z2 is H g It can be expressed as:

[0091]

[0092] Under the principle of constant stratum skeleton volume, H is set g Not affected by burial depth, that is:

[0093]

[0094] If the hole-depth function of a rock layer is known, the burial depth Z can be calculated by solving the equation using the iterative method. ′ Bottom burial depth Z at 1 o'clock ′ 2. Calculate the original thickness of the stratum and the thickness at any depth during the compaction process.

[0095] As an exemplary embodiment, in this embodiment, the initial porosity of mudstone is 50% and that of sandstone is 42%. The compaction curve of Permian-Jurassic mudstone is calculated as The compaction curve of sandstone is: The stripping thickness of multiple wells was calculated and averaged, and the compaction ratio of the Permian-Jurassic mudstone was found to be 0.55, and the compaction ratio of the sandstone was 0.42.

[0096] Referring to Table 1, different sedimentary facies have different lithologies and related physical parameters. Combined with the previous research results of trace elements and biomarkers, the following paleowater depth reference values are obtained.

[0097] Table 1: Paleowater depth parameters for different sedimentary facies

[0098]

[0099] Finally, the residual paleo-geomorphology is superimposed, the paleo-geomorphology of the prototype basin is restored, and the prototype basin restoration result map is generated, such as Figure 8 shown.

[0100] This embodiment provides a computer-readable storage medium having a computer program stored thereon. When executed by a processor, the computer program implements the steps of the above-described method for restoring a prototype basin during the depression-uplift phase. The storage medium may be a magnetic disk, an optical disk, a read-only memory (ROM), a random access memory (RAM), a flash memory, a hard disk drive (HDD), or a solid-state drive (SSD); the storage medium may also include a combination of the aforementioned types of memory.

[0101] This embodiment 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 program, the steps of the above-mentioned method for restoring a prototype basin in the depression-uplift stage of a basin are implemented.

[0102] like Figure 9As shown, the computer device 120 may include: at least one processor 121, such as a central processing unit (CPU), at least one communication interface 123, a memory 124, and at least one communication bus 122. The communication bus 122 is used to realize the connection and communication between these components. The communication interface 123 may include a display screen and a keyboard. The optional communication interface 123 may also include a standard wired interface and a wireless interface. The memory 124 may be a high-speed random access memory (RAM) or a non-volatile memory, such as at least one disk storage. The memory 124 may optionally be at least one storage device located away from the aforementioned processor 121. The memory 124 stores application programs, and the processor 121 calls the program code stored in the memory 124 to execute any of the above-mentioned method steps. The communication bus 122 may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The communication bus 122 can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 9Only one line is used to represent it, but it does not mean that there is only one bus or one type of bus. The memory 124 may include volatile memory, such as random-access memory (RAM); the memory may also include non-volatile memory, such as flash memory, hard disk drive (HDD) or solid-state drive (SSD); the memory 124 may also include a combination of the above types of memory. The processor 121 may be a central processing unit (CPU), a network processor (NP) or a combination of a CPU and an NP. The processor 121 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD) or a combination thereof. The PLD may be a complex programmable logic device (CPLD), a field-programmable gate array (FPGA), a generic array logic (GAL) or any combination thereof. Optionally, the memory 124 is further configured to store program instructions. The processor 121 may call the program instructions to implement the method for restoring a prototype basin in the depression-uplift stage of the present embodiment.

[0103] The embodiments of the present invention are described above in conjunction with the accompanying drawings, but the present invention is not limited to the above-mentioned specific implementation methods. The above-mentioned specific implementation methods are merely illustrative and not restrictive. Under the guidance of the present invention, ordinary technicians in this field can also make many forms without departing from the scope of protection of the present invention and the claims, all of which are protected by the present invention.

Claims

1. A method for restoring a prototype basin during the depression-uplift stage of a basin, characterized in that: The following steps are involved: S1: Based on the basin boundary identification map, perform basin boundary identification and unconformity surface analysis to determine the erosion period, erosion point, and stratigraphic extension trend; S2: generating a residual thickness map based on the erosion period, erosion point, stratigraphic extension trend and seismic logging data; S3: Based on the residual thickness map, perform denudation recovery and perform denudation interpretation on the seismic profile to obtain the original sedimentary thickness of the stratum; S4: Based on the original sedimentary thickness of the stratum, compaction correction and paleo-water depth correction are performed to obtain the accurate original thickness of the stratum.

2. The method for restoring a prototype basin in the depression-uplift stage according to claim 1, characterized in that: The erosion amount is recovered as follows: Where K is the deposition rate coefficient, which is expressed as the ratio of the deposition rate of the recovery layer to the deposition rate of the reference layer. Usually, wells with stable deposition in the depression area and no erosion are selected for calculation. As the deposition rate accelerates or slows down, the K value can be greater than 1 or less than 1; ΔH A is the erosion thickness of the unconformity surface; H A 、T A are the sedimentary thickness and sedimentary time of the unconformity surface, respectively; H B is the sedimentary thickness of the stratum.

3. The method for restoring a prototype basin in the depression-uplift stage according to claim 1, characterized in that: Step S4 specifically includes: performing compaction correction and paleo-water depth correction according to the original sedimentary thickness of the stratum to obtain the accurate original thickness of the stratum, as shown in the formula: in, is the porosity of clastic rocks, is the surface porosity, z is the burial depth of the rock layer, c is the compaction factor, Z1 and Z2 are the top and bottom burial depths of the stratum after compaction, respectively. ′ 1 and Z ′ 2 is the top surface burial depth and bottom surface burial depth.

4. The method for restoring a prototype basin in the depression-uplift stage according to claim 1, characterized in that: According to the original sedimentary thickness of the stratum, compaction correction and paleo-water depth correction are performed to obtain the accurate original thickness of the stratum, as shown in the formula: in, is the porosity of clastic rocks, is the surface porosity, z is the burial depth of the rock layer, k is the compaction factor, Z1 and Z2 are the top and bottom burial depths of the stratum after compaction, respectively. ′ 1 and Z ′ 2 is the top surface burial depth and bottom surface burial depth.

5. A device suitable for restoring a prototype basin during the depression-uplift stage of a basin, characterized in that: The device for restoring a prototype basin in the depression-uplift stage of a basin includes the following modules: The basin boundary identification module is configured to: perform basin boundary identification and unconformity surface analysis based on the basin boundary identification map, and determine the erosion period, erosion point, and stratigraphic extension trend; A residual thickness map generating module is configured to generate a residual thickness map based on the erosion period, erosion point, formation extension trend and seismic logging data; an erosion recovery module configured to: recover the erosion amount based on the residual thickness map, perform erosion interpretation on the seismic profile, and form the original sedimentary thickness of the stratum; The correction module is configured to perform compaction correction and paleo-water depth correction according to the original sedimentary thickness of the stratum to obtain the accurate original thickness of the stratum.

6. The device for restoring a prototype basin in the depression-uplift stage according to claim 5, characterized in that: The erosion amount is recovered as follows: Where K is the deposition rate coefficient, which is expressed as the ratio of the deposition rate of the recovery layer to the deposition rate of the reference layer. Usually, wells with stable deposition in the depression area and no erosion are selected for calculation. As the deposition rate accelerates or slows down, the K value can be greater than 1 or less than 1; ΔH A is the erosion thickness of the unconformity surface; H A 、T A are the sedimentary thickness and sedimentary time of the unconformity surface, respectively; H B is the sedimentary thickness of the stratum.

7. The device for restoring a prototype basin in the depression-uplift stage according to claim 5, characterized in that: The specific configuration of the correction module is: according to the original sedimentary thickness of the stratum, compaction correction and paleowater depth correction are performed to obtain the accurate original thickness of the stratum, as shown in the formula: in, is the porosity of clastic rocks, is the surface porosity, z is the burial depth of the rock layer, c is the compaction factor, Z1 and Z2 are the top and bottom burial depths of the stratum after compaction, respectively. ′ 1 and Z ′ 2 is the top surface burial depth and bottom surface burial depth.

8. The device for restoring a prototype basin in the depression-uplift stage according to claim 5, characterized in that: The specific configuration of the correction module is: according to the original sedimentary thickness of the stratum, compaction correction and paleowater depth correction are performed to obtain the accurate original thickness of the stratum, as shown in the formula: in, is the porosity of clastic rocks, is the surface porosity, z is the burial depth of the rock layer, k is the compaction factor, Z1 and Z2 are the top and bottom burial depths of the stratum after compaction, respectively. ′ 1 and Z ′ 2 is the top surface burial depth and bottom surface burial depth.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for restoring a prototype basin in the depression-uplift stage of a basin as described in any one of claims 1 to 4 are implemented.

10. 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 steps for restoring a prototype basin in the basin depression-uplift stage as described in any one of claims 1 to 4 are implemented.

Images