Methods, devices, electronic equipment and storage media for the compaction and restoration of ancient landforms and ancient structures

Abstract

This invention relates to a method, apparatus, electronic device, and storage medium for the compaction and restoration of paleomorphology and paleostructures. The embodiments pertain to the field of geological engineering, and disclose a method, apparatus, electronic device, and storage medium for the compaction and restoration of paleomorphology and paleostructures. The method includes acquiring target layer depth data, performing iterative processing using the steepest descent method on a point-by-point basis based on the target layer depth data to obtain the stratigraphic thickness, and performing planar mapping based on the stratigraphic thickness to obtain a restored true paleomorphological and paleostructural plan view. This method, applied to the compaction and restoration of paleomorphology and paleostructures, enables rapid iteration and accurate restoration of the target layer thickness. It solves the problems of complex processes and significant discrepancies in the restored stratigraphic data found in existing compaction and restoration techniques.

Description

Technical Field

[0001] This invention is applicable to technical fields such as geological engineering in oil and gas field exploration and coal field exploration, and in particular to a method, device, electronic equipment and storage medium for the compaction and restoration of ancient landforms and ancient structures. Background Technology

[0002] Compaction restoration is a technical method applied in geological engineering fields such as oil and gas field exploration and coalfield exploration. Because strata have undergone long-term sedimentary evolution, some parts have been altered by compaction, affecting the prediction of oil and gas reservoir characteristics. Therefore, compaction restoration is needed to reconstruct paleomorphology and paleotectonic structures. Existing compaction restoration technologies use well logging curves to calculate the average value of lithological compaction characteristics to obtain paleomorphological restoration data. The inventors discovered that existing compaction restoration methods fail to consider practical applications, such as the varying compaction coefficients of different rocks within a stratum. Therefore, the paleomorphology and paleotectonic structures reconstructed by existing technologies are not accurate. Summary of the Invention

[0003] The purpose of this invention is to provide a method, apparatus, electronic device, and storage medium for compaction restoration of ancient landforms and structures, which can accurately perform compaction restoration of ancient landforms and structures.

[0004] To address the aforementioned technical problems, embodiments of the present invention provide a method for the compaction and restoration of paleomorphology and paleotectonic structures, comprising: acquiring target layer data; performing layer-by-layer recursive calculations based on the target layer data; iteratively obtaining the thickness of the strata before compaction using the steepest descent method for each sample point; and performing planar mapping based on the strata thickness to obtain a restored paleomorphology and paleotectonic planar map.

[0005] Embodiments of the present invention also provide a data processing apparatus, comprising:

[0006] The acquisition module is used to acquire geological stratigraphic data such as seismic data and well data that need to be processed and calculated.

[0007] The processing module is used to perform necessary processing on the formation data;

[0008] The sending module is used to send the acquired data to other terminals;

[0009] The display module is used to display the processed data, wherein the display methods include, but are not limited to, text and images;

[0010] The processing module includes:

[0011] The profile unit is used to filter and organize well and seismic data to obtain the required seismic profile and well-connected profile.

[0012] The division unit is used to uniformly divide the geological strata after drilling. Based on the comprehensive analysis of data such as well profiles, seismic data and logging curves, the formation depth is determined.

[0013] The flattening unit is used to process the reference surface. It selects a suitable seismic reflection interface for flattening to obtain the reference surface and the erosion of the underlying strata that are closest to the paleomorphology and paleotectonic structure to be restored.

[0014] The calculation unit is used to perform calculations based on the calculation formula and the target layer data to obtain the restored paleogeographic and paleotectonic data.

[0015] Embodiments of the present invention also provide an electronic device, comprising:

[0016] At least one processor; and,

[0017] A memory communicatively connected to the at least one processor; wherein,

[0018] The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the data processing method as described in claim 1.

[0019] The present invention also provides a computer-readable storage medium storing a computer program, characterized in that the computer program, when executed by a processor, implements the above-described method for the compaction and restoration of paleomorphology and paleotectonic structures.

[0020] Compared with the prior art, the embodiments of the present invention, when performing compaction restoration of the target layer, take into account the different compaction coefficients of different rocks and derive a new iterative format using the steepest descent method. This allows the equation to iterate quickly and accurately restore the target layer, solving the problems of slow iteration speed of the integral equation in the prior art, different compaction coefficients of different rocks contained in the strata, and the lack of calculation by lithology for sub-layers, resulting in inaccurate paleomorphology and paleotectonics restored by compaction.

[0021] Furthermore, in the paleogeographic and paleotectonic compaction restoration method provided by the embodiments of the present invention, the steepest descent method is expressed as: x(k+1) => x(k) - 0.5*f / (df / dx), where f = x + ae -c2x -b,df / dx=1-a·c2e -c2·x x = d²r The porosity of a certain stratum at the Earth's surface. Let c1 be the porosity of a certain stratum underground, c2 be the compaction coefficient of the sediment, and d = (d2 - d1) / 2.

[0022] Furthermore, the paleogeographic and paleotectonic compaction restoration method provided in this invention includes target layer data such as: the current top boundary depth of the target layer, the current bottom boundary depth of the target layer, surface porosity, subsurface porosity, compaction coefficient of sediments, wave impedance, PE parameters, formation P-wave velocity, and density; and, based on the target layer data, processing yields a reference surface, seismic interface of the study area, depth-domain porosity data volume, inversion profile of PE parameters, wave impedance profile, porosity profile, and time-domain pre-stack migration profile.

[0023] In addition, the paleogeographic and paleostructural compaction restoration method provided in this embodiment of the invention, wherein the recursive calculation of the target layer point by point to obtain the thickness of the stratum before compaction further includes: obtaining the depth domain porosity profile after compaction restoration and a text file containing the x, y coordinates of the geographical location, the thickness value of the stratum after compaction restoration, and the total compaction coefficient by performing recursive calculation of the stratum point by point.

[0024] In addition, the paleogeographic and paleotectonic compaction restoration method provided in the embodiments of the present invention includes the following steps for the reference surface: processing the target layer data according to the target layer data to obtain a seismic profile; obtaining paleotectonic undulations according to the seismic profile; and selecting a seismic reflection interface to fill and complete the paleotectonic undulations to obtain the completed reference surface.

[0025] In addition, the paleogeography and paleotectonic compaction restoration method provided in the embodiments of the present invention includes the following steps: processing the completed wells in the study area according to the stratigraphic data of the target section to obtain a continuous well profile, uniformly dividing the geological strata, and performing seismic interpretation based on the divided geological strata to obtain the seismic interface of the study area.

[0026] In addition, the paleogeographic and paleostructural compaction restoration method provided in this embodiment of the invention includes the following steps: based on the target layer data, applying well logging curves and seismic data to invert and obtain a porosity data volume; and based on the porosity data volume, performing a depth domain conversion to obtain a depth domain porosity data volume. Attached Figure Description

[0027] One or more embodiments are illustrated by way of example with reference numerals in the accompanying drawings. These illustrations do not constitute a limitation on the embodiments. Elements with the same reference numerals in the drawings are denoted as similar elements. Unless otherwise stated, the figures in the drawings are not to be limited by scale.

[0028] Figure 1This is a flowchart of the paleogeographic and paleostructural compaction restoration method provided by the embodiments of the present invention;

[0029] Figure 2 This is the data processing flow provided by the embodiments of the present invention. Figure 1 ;

[0030] Figure 3 This is the data processing flow provided by the embodiments of the present invention. Figure 2 ;

[0031] Figure 4 This is the data processing flow provided by the embodiments of the present invention. Figure 3 ;

[0032] Figure 5 This is a pre-stack migration profile in the time domain, where T K T J7 T J These are the Cretaceous basal surface, the Jurassic Yan'an Formation Yan 7 basal surface, and the Jurassic basal surface, respectively. J The lower part is the Triassic Yanchang Formation;

[0033] Figure 6 This is the depth-domain porosity profile obtained from seismic inversion, where D K D J7 D J These are the Cretaceous basal surface, the Jurassic Yan'an Formation Yan 7 basal surface, and the Jurassic basal surface in the depth domain, respectively. J The lower part is the Triassic Yanchang Formation;

[0034] Figure 7 To extend the porosity profile of the depth domain after compaction recovery in Phase 7, D J It is the base of the Jurassic system;

[0035] Figure 8 This is a porosity profile of the depth domain after Cretaceous compaction recovery, Dr J This represents the Jurassic floor surface after compaction and restoration during the Cretaceous period.

[0036] Figure 9 Compaction and restoration of the ancient landform map of Phase 7;

[0037] Figure 10 Paleogeographic maps of Phase 7 after compaction and restoration;

[0038] Figure 11 Paleotectonic map of the Cretaceous period before compaction restoration;

[0039] Figure 12 The paleotectonic map of the Cretaceous period after compaction and restoration, where w1 and w2 are industrial oil wells, w3 is a water well, the black circles of the well points indicate the wellhead locations, and the endpoint of the line segment connecting them is the bottom of the Jurassic period.

[0040] Figure 13 This is a schematic diagram of the terminal structure provided in an embodiment of the present invention;

[0041] Figure 14 yes Figure 13 The processing module 1303 in the terminal device provided in the embodiment of the present invention shown

[0042] Figure 15 This is a schematic diagram of the structure of an electronic device provided by an embodiment of the present invention. Detailed Implementation

[0043] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the various embodiments of the present invention will be described in detail below with reference to the accompanying drawings. However, those skilled in the art will understand that many technical details are presented in the various embodiments of the present invention to facilitate a better understanding of this application. However, the technical solutions claimed in this application can be implemented even without these technical details and various changes and modifications based on the following embodiments. The division of the various embodiments below is for ease of description and should not constitute any limitation on the specific implementation of the present invention. The various embodiments can be combined with and referenced by each other without contradiction.

[0044] The embodiments of the present invention relate to a method for the compaction and restoration of paleomorphology and paleotectonic structures, which is applied in the process of geological engineering exploration. The process is as follows: Figure 1 As shown, it includes:

[0045] Step 101: Obtain the target layer data.

[0046] In this embodiment, the acquired data may include: the current top boundary depth of the target layer, the current bottom boundary depth of the target layer, surface porosity, subsurface porosity, compaction coefficient of sediments, wave impedance, PE parameters, formation P-wave velocity, density, and PE logging data. Of course, the above are merely specific examples; in actual use, the data and seismic data may include other content, which will not be elaborated upon here.

[0047] This embodiment does not limit the data and seismic data. In actual use, the data and seismic data can include more data required for compaction recovery, depending on the specific characteristics of the destination layer.

[0048] Step 102: Based on the target layer data, perform a steepest descent method iteration for each sample point to obtain the formation thickness.

[0049] In this embodiment, the target layer data, in addition to the current top boundary depth, current bottom boundary depth, surface porosity, subsurface porosity, sediment compaction coefficient, wave impedance, PE parameters, formation P-wave velocity, density, and PE logging data obtained in step 101, also includes the following data obtained after filtering and processing: a reference surface, seismic interface of the study area, depth-domain porosity data volume, PE parameter inversion profile, wave impedance profile, porosity profile, and time-domain pre-stack migration profile. The time-domain pre-stack migration profile is shown below. Figure 5 As shown, the time-domain pre-stack migration profile is a result of seismic data processing, and its role is that it is one of the most effective methods for imaging complex structures.

[0050] The reference surface processing procedure is as follows: Figure 2 As shown, the seismic interface processing process in the study area is as follows: Figure 3 As shown, the depth domain porosity data volume processing process is as follows: Figure 4 As shown. The specific process will be explained later, and will not be elaborated on here.

[0051] Based on the above data, recursive calculations are performed to obtain the steepest descent method iterative format. Then, the obtained steepest descent method iterative calculations are performed using this format to obtain the restored paleogeographic and paleotectonic data. The restored paleogeographic and paleotectonic data includes a compacted and restored depth-domain porosity profile and a text file containing the x, y coordinates of the geographical location, the compacted and restored thickness values ​​of the strata, and the total compaction coefficient.

[0052] The steepest descent method iterative format is x(k+1) => x(k) - 0.5*f / (df / dx), and its derivation process includes: setting assumptions based on existing technology: compaction reduces the porosity of rocks without changing the size of the skeleton particles; various secondary changes in diagenesis are not considered; compaction is an irreversible process, that is, the thickness of the strata after restoration must be greater than or equal to the thickness of the compacted strata, in order to facilitate program design and simplify the established model. Under the above assumptions, let d1 and d2 represent the current top and bottom burial depths of a certain set of clastic rock strata, respectively. After stripping the overlying strata, the top and bottom depths are d1r and d2r, respectively. The equation is derived as: f = x + ae -c2x -b = 0, differentiate the above equation, define the objective function based on the result of the differentiation, and finally obtain the gradient method iterative format. Where: x = d²r; The porosity of a certain stratum at the Earth's surface. Let c1 be the porosity of a certain stratum underground, c2 be the compaction coefficient of the sediment, and d = (d2 - d1) / 2. Different lithologies have different compaction coefficients; mudstone typically has a large compaction coefficient, generally (0.5–1.0) × 10⁻⁶. -3 m -1 The compaction coefficient of sandstone is relatively small compared to that of mudstone, mostly ranging from (0.1 to 0.5) × 10⁻⁶. -3 m -1 Carbonate rocks have a compaction coefficient between that of mudstone and sandstone. The compaction coefficient of coal-bearing strata is taken as that of mudstone; the compaction coefficient of evaporites is generally taken as zero.

[0053] Specifically, how is formation recovery performed point-by-point based on the gradient method iterative format?

[0054] When the target interval is a single lithological interval, let d1r = 0; d2r = x(1);

[0055] When the target stratigraphic segment consists of multiple lithological segments, it is determined recursively. For example:

[0056] The first layer: d1r = 0; d2r = x(1);

[0057] The second layer: d1r = x(1); d2r = x(2);

[0058] ...

[0059] For the (n-1)th layer, d1r = x(n-2); d2r = x(n-1);

[0060] The nth layer has d1r = x(n-1) and d2r = x(n).

[0061] It should be noted that the porosity profile and wave impedance profile obtained from seismic inversion, and if it includes carbonate rocks, the inversion profile with parameters such as PE, are used as inputs for compaction recovery calculation. The inversion profiles with wave impedance and PE parameters are used to determine lithology, and porosity, as a parameter for compaction recovery calculation, is also used to determine rock type.

[0062] By considering the different compaction coefficients of different rock strata, the derived equation contains more detailed information and can be applied to different types of target strata, making the technical solution provided by the embodiments of the present invention more accurate.

[0063] By utilizing the gradient method to obtain the gradient method iteration format, the equations become more concise, and the iteration speed of the technical solution provided by the embodiments of the present invention can be faster.

[0064] Of course, the above is just a simple explanation. In actual use, step 102 may involve many different data information for processing, which will not be elaborated here.

[0065] Step 103: Based on the thickness of the strata, perform planar mapping to obtain the restored landform plan.

[0066] In this embodiment, a planar mapping is performed based on the stratigraphic thickness values ​​calculated in the previous step to obtain a reconstructed paleogeographic plan. For example, as... Figure 12 The compacted and restored paleotectonic diagram of the Cretaceous period shown can be compared with, for example... Figure 11 The paleotectonic diagrams of the Cretaceous period before compaction restoration are shown for comparison. For example... Figure 9 The paleogeographic map shown above, which extends seven phases before compaction restoration, can be compared with, for example... Figure 10 The above shows the formation comparison of paleogeographic maps after compaction restoration in seven phases.

[0067] Additionally, it should be noted that the mapping also includes depth-domain porosity profiles restored after Cretaceous compaction, which can be compared with depth-domain porosity profiles before compaction restoration. The Cretaceous period was chosen because other geological studies have already established that the hydrocarbon generation-migration period in this region was the Cretaceous. The paleotectonic morphology before and after compaction can be analyzed from both the depth-domain porosity profiles and planar views.

[0068] In this embodiment, such as Figure 2 As shown, it includes:

[0069] Step 201: Based on the seismic and well data of the target section, process the seismic data of the target section to obtain a seismic profile.

[0070] In this embodiment, the seismic and well data of the target layer can be displayed in a three-dimensional manner, which is more intuitive and easier to process.

[0071] Step 202: Based on the seismic profile, obtain the paleotectonic undulations.

[0072] In this embodiment, based on the seismic profile, the sedimentary base surface, the maximum floodplain, or the peneplain is used as a reference surface for leveling. The underlying or overlying strata are then manipulated to obtain the structural morphology, which approximates the relative paleotectonic undulations of that period. After sedimentation, the target strata undergo long-term geological movement, experiencing compaction and tectonic deformation, folding, and fracturing by the overlying strata. To restore the paleomorphology and paleotectonic morphology, a suitable base surface is selected for leveling before compaction restoration, which can eliminate the influence of tectonic deformation.

[0073] Step 203: Based on the paleotectonic undulations, select the seismic reflection interface for leveling and filling to obtain the filled reference surface.

[0074] In this embodiment, by browsing and interpreting seismic data, 2-3 seismic reflection interfaces that conform to the fill-in reference surface are selected. These seismic reflection interfaces can be selected from those with relatively stable paleogeographic overlying morphology; alternatively, the maximum floodplain in the underlying strata of the paleogeographic overlying strata can be selected as the fill-in reference surface to study the erosion of the underlying strata; or, the reflection interface corresponding to the main hydrocarbon expulsion period can be selected to study the paleotectonic morphology of the main hydrocarbon expulsion period.

[0075] It should be noted that the browsing and interpretation of seismic data mentioned in this invention are not specifically limited. In practice, different interpretations can be made depending on the target layer. The selection of the reference surface to be supplemented also depends on the specific needs of the implementation process.

[0076] In this embodiment, such as Figure 3 As shown, it includes:

[0077] Step 301: Based on the seismic and well data of the target section, process the completed wells in the study area to obtain the well profile.

[0078] In this embodiment, the seismic and well data of the target layer can be displayed in a three-dimensional manner, which is more intuitive and easier to process.

[0079] It should be noted that the browsing and interpretation of seismic data mentioned in this invention are not specifically limited. In practice, different interpretations can be made depending on the target layer. The selection of the reference surface to be supplemented also depends on the specific needs of the implementation process.

[0080] Step 302: Based on the well-connected profile, perform unified geological strata division on the completed wells to obtain the divided well-connected profile.

[0081] In this embodiment, geological strata are delineated based on the acquired seismic and well data, using both seismic and strata data. Regarding strata delineation, it typically extends to the lower boundary of the sandstone, with the corresponding seismic profile showing the peaks.

[0082] Step 303: Based on the divided well profile, perform seismic interpretation to obtain the seismic interface of the study area.

[0083] In this embodiment, such as Figure 4 As shown, it includes:

[0084] Step 401: Based on the target layer data, well logging curves and seismic data are used to invert and obtain the porosity data volume.

[0085] Step 402: Perform depth domain conversion based on the porosity data volume to obtain a depth domain porosity data volume.

[0086] In this embodiment, the porosity profile before compaction recovery is as follows: Figure 9 As shown.

[0087] In this embodiment, such as Figure 13 As shown, it includes a data processing device:

[0088] The acquisition module 1301 is used to acquire data and stratigraphic information that need to be processed and calculated.

[0089] Processing module 1302 is used to perform necessary processing on the formation data;

[0090] It should be noted that the browsing and interpretation of stratigraphic data mentioned in this invention are not specifically limited. In practice, different interpretations can be made depending on the target stratigraphic segment. The selection of the reference surface to be supplemented also depends on the specific needs of the implementation process.

[0091] The sending module 1303 is used to send the acquired data to other terminals;

[0092] Display module 1304 is used to display the processed data, wherein the display method includes, but is not limited to, text, images, etc.

[0093] In this embodiment, such as Figure 14 As shown, Figure 13 The processing module in the device shown includes:

[0094] Profile unit 1401 is used to filter and organize data to obtain the required seismic profiles and well-connected profiles.

[0095] Division unit 1402 is used for the unified division of geological strata after drilling. Based on the well profile, seismic data and logging curves, the stratigraphic boundaries of the target section are determined.

[0096] The completion unit 1403 is used to process the reference surface, select a suitable seismic reflection interface for completion, and obtain the overlying reference surface that is closest to the paleomorphology and paleotectonic structure to be restored.

[0097] The calculation unit 1404 is used to perform calculations based on the calculation formula and the target layer data to obtain the restored paleogeographic and paleotectonic data.

[0098] This embodiment relates to an electronic device, such as... Figure 15 As shown, it includes:

[0099] At least one processor 1501;

[0100] and a memory 1502 communicatively connected to the at least one processor 1501;

[0101] The memory 1502 stores instructions that can be executed by the at least one processor, which are executed by the at least one processor 1501 to enable the at least one processor 901 to perform the paleomorphology and paleotectonic compaction restoration method according to the first to third embodiments of the present invention.

[0102] The memory and processor are connected via a bus, which can include any number of interconnecting buses and bridges, connecting various circuits of one or more processors and memories. The bus can also connect various other circuits, such as peripheral devices, voltage regulators, and power management circuits, which are well known in the art and will not be described further herein. The bus interface provides an interface between the bus and the transceiver. The transceiver can be a single element or multiple elements, such as multiple receivers and transmitters, providing a unit for communicating with various other devices over a transmission medium. Data processed by the processor is transmitted over the wireless medium via an antenna, which further receives data and transmits it to the processor.

[0103] The processor manages the bus and general processing, and also provides various functions, including timing, peripheral interfaces, voltage regulation, power management, and other control functions. Memory is used to store data used by the processor during operation.

[0104] This embodiment relates to a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, it implements the method embodiments described above.

[0105] That is, those skilled in the art will understand that all or part of the steps in the methods of the above embodiments can be implemented by a program instructing related hardware. This program is stored in a storage medium and includes several instructions to cause a device (which may be a microcontroller, chip, etc.) or processor to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, a portable hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0106] Those skilled in the art will understand that the above embodiments are specific examples of implementing the present invention, and in practical applications, various changes in form and detail may be made without departing from the spirit and scope of the present invention.

Claims

1. A palaeogeomorphology and palaeotectonic compaction restoration method, characterized in that, The method comprises the following steps: obtaining a plurality of target layer data and corresponding thickness data; performing steepest descent method iteration on the target layer data to obtain stratum thickness; performing plane mapping according to the stratum thickness to obtain restored palaeogeomorphology and palaeostructure plane maps; wherein the target layer data comprises: current top boundary depth of the target layer, current bottom boundary depth of the target layer, surface porosity, underground porosity, compaction coefficient of sediments, wave impedance, PE parameter, stratum P-wave velocity, density, PE logging data, and base surface, seismic interface of the study area, depth domain porosity data body, PE parameter inversion profile, wave impedance profile, porosity profile, and time domain pre-stack migration profile obtained by screening and processing the current top boundary depth of the target layer, the current bottom boundary depth of the target layer, the surface porosity, the underground porosity, the compaction coefficient of sediments, the wave impedance, the PE parameter, the stratum P-wave velocity, the density, and the PE logging data; wherein the steepest descent method iteration is as follows: x(k+1) = x(k) - 0.5*f / (df / dx), wherein, f = x+ae -c2x -b,df / dx=1-a·c2e -c2·x , x = d2r ; a = φ0 / c 2， b = [d2 - d1] + φ0[e -c1·d2 -e -c1·d1 ] / c1 + d1r + (φ0 / c2)e -c2·d1r , c1 = ln(φ0 / φ) / d, φ0 is the porosity of a certain formation at the surface, φ is the porosity of a certain formation at the underground, c2 is the compaction coefficient of the sediment, d = (d2 - d1) / 2, d1 and d2 respectively represent the current top boundary and bottom boundary depth of a set of clastic rock formation.

2. The method of claim 1, wherein, the steepest descent method iteration on the target layer data to obtain stratum thickness comprises: restored palaeogeomorphology and palaeostructure data; the restored palaeogeomorphology and palaeostructure data comprise depth domain porosity profile after compaction restoration and a text file containing geographical position x, y coordinates, stratum thickness value after compaction restoration, and total compaction coefficient.

3. The method of claim 1, wherein, the base surface is generated in the following manner: obtaining a seismic data body by seismic data processing according to newly collected original seismic data; completing well-to-seismic calibration, well-to-seismic layering after calibration, obtaining unified well layering, and then obtaining seismic interpretation horizon of the main target layer through interpretation according to the seismic data body; obtaining the structural form of the underlying or overlying stratum, i.e. approximately representing the relative palaeostructure relief at this period, by flattening the reference surface of the sedimentary base surface, the maximum flooding surface or the erosion surface, and selecting the flattened and filled surface corresponding to the seismic reflection interface as the base surface.

4. The method of claim 1, wherein, the seismic interface of the study area is generated in the following manner: obtaining a continuous well profile by calibrating a drilled well in the study area according to the target layer data; obtaining a divided continuous well profile by dividing the unified geological horizon according to the continuous well profile; obtaining the seismic interface of the study area by seismic interpretation according to the divided continuous well profile.

5. The method of claim 1, wherein, the depth domain porosity data body is generated in the following manner: obtaining a porosity data body by inversion according to the target layer data and by applying logging curves and seismic data; obtaining the depth domain porosity data body by depth domain conversion according to the porosity data body.

6. An apparatus for palaeogeomorphology and palaeotectonic compaction restoration, characterized in that, The method comprises the following steps: obtaining a plurality of target layer data and corresponding thickness data; performing steepest descent method iteration on the target layer data to obtain stratum thickness; performing plane mapping according to the stratum thickness to obtain restored palaeogeomorphology and palaeostructure plane maps; wherein the target layer data at least comprises: a current top boundary depth of the target interval, a current bottom boundary depth of the target interval, a surface porosity, a subsurface porosity, a compaction coefficient of the sediment, a wave impedance, a PE parameter, a formation P-wave velocity, a density, PE logging data, and a datum surface obtained by screening and processing the current top boundary depth of the target interval, the current bottom boundary depth of the target interval, the surface porosity, the subsurface porosity, the compaction coefficient of the sediment, the wave impedance, the PE parameter, the formation P-wave velocity, the density, and the PE logging data, a seismic interface of the study area, a depth domain porosity data volume, an inversion profile of the PE parameter, a wave impedance profile, a porosity profile, and a time domain prestack migration profile; wherein the steepest descent method format is: x(k+1) => x(k) - 0.5*f / (df / dx), wherein, f = x + ae -c2x -b, df / dx = 1 - a - c2e -c2·x , x = d2r ; a = φ0 / c 2， b = [d2 - d1] + φ0[e -c1·d2 - e -c1·d1 ] / c1 + d1r + (φ0 / c2)e -c2·d1r , c1 = ln(φ0 / φ) / d, φ0 is the porosity of a certain formation at the surface, φ is the porosity of a certain formation at the underground, c2 is the compaction coefficient of the sediment, d = (d2 - d1) / 2, d1 and d2 respectively represent the current top boundary and bottom boundary depth of a set of clastic rock formation.

7. An electronic device, comprising: comprise: at least one processor; and a memory connected with the at least one processor in communication; wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to perform the paleogeomorphology and paleostructure compaction restoration method as claimed in claim 1.

8. A computer readable storage medium storing a computer program, characterized in that, The computer program is executed by the processor to implement the paleogeomorphology and paleostructure compaction restoration method as claimed in claim 1.

Images