Sand shale stratum micro-ancient landform depicting method and related equipment

Abstract

The invention provides a sand shale stratum micro-ancient landform depicting method, which comprises the following steps of: determining an upper datum plane according to an overlying stratum of a to-be-recovered interface, calculating impression thickness, establishing a sedimentary strata sequence framework according to a destination rock core and logging curve characteristics, extracting sand-to-ground ratio of each well section and a porosity mean value of sand shale, and calculating the thickness of the sand shale stratum micro-ancient landform. Estimating a multi-well lithology difference compaction correction value based on the sand-to-ground ratio of each well section and the porosity mean value of the sand shale, obtaining lithology difference compaction correction value distribution according to the impression thickness and the multi-well lithology difference compaction correction value, correcting the impression thickness based on the lithology difference compaction correction value distribution, obtaining corrected impression thickness distribution, and obtaining the impression thickness distribution. The corrected impression thickness distribution is converted, the elevation plane of the interface to be recovered is obtained, and a method foundation is laid for research on sedimentary micro-ancient landform recovery for constructing a relatively stable area by means of the method. The fine depicted sedimentary micro ancient landform spatial form can provide an important basis for distribution of a sediment source and a local thin mutual sand body and analysis of a time-space superposition relationship.

Description

Technical Field

[0001] The invention belongs to the field of seismic sedimentology, and specifically relates to a sandstone-mudstone stratum micro-paleomorphology characterization method and related equipment. Background Art

[0002] In areas with relatively stable structures (obvious erosion), micro-paleomorphological restoration and detailed characterization of the target layer are carried out on thin interbedded sandstones and mudstones in continental and marine-continental transitional phases, which can provide support and basis for the analysis of the provenance direction and settlement center of the target layer sediments and the analysis of the distribution characteristics of the sediments.

[0003] Under the control of the well-connected sedimentary sub-(micro) phase contrast model, tectonic erosion is not considered for the time being. Based on the top and bottom boundaries of the target sedimentary strata, the existing stamp method can only restore the overall distribution trend of the paleo-geomorphology. In fact, due to the influence of lithological changes, local strata often produce later differential compaction effects, which makes it difficult to restore and describe the micro paleo-geomorphological characteristics of the local area. In addition, due to the large vertical and lateral changes in the lithology of the thin interbedded sandstone and mudstone, the existing geostatistical interpolation method based on multi-well data is not accurate in predicting lithological characteristic parameters, and cannot meet the lateral description of micro paleo-geomorphological characteristics. Summary of the invention

[0004] The present invention provides a sandstone-mudstone stratum micro-paleomorphology characterization method and related equipment to solve the problem that the existing geological statistical interpolation method based on multi-well data has low prediction accuracy for lithological characteristic parameters and cannot meet the requirements of lateral characterization of micro-paleomorphology characteristics.

[0005] To achieve the above object, the present invention provides the following technical solutions:

[0006] A method for depicting micro-paleomorphology of sandstone and mudstone strata, comprising:

[0007] Determine the upper reference surface according to the overlying strata of the interface to be restored, and calculate the thickness of the impression;

[0008] According to the characteristics of the target core and well logging curves, a sedimentary sequence framework is established to extract the sand-to-formation ratio and the mean porosity of sandstone and mudstone in each well section;

[0009] The multi-well lithology difference compaction correction amount is estimated based on the average value of sand-to-shale ratio and sand-mudstone porosity in each well section;

[0010] The distribution of lithology difference compaction correction is obtained according to the stamp thickness and the multi-well lithology difference compaction correction;

[0011] The mold thickness is corrected based on the distribution of compaction correction amount of lithology difference to obtain the corrected mold thickness distribution, and the corrected mold thickness distribution is converted to obtain the elevation plane of the interface to be restored.

[0012] Preferably, the upper reference plane is determined according to the overlying strata on the interface to be restored, and the mold thickness is calculated as follows: a suitable stable interface is found in the overlying strata on the interface to be restored to be determined as the upper reference plane, and the coinciding axes of the interface to be restored and the upper reference plane are picked up from the seismic section through fine stratigraphic calibration, and the thickness between the interface to be restored and the upper reference plane, i.e., the mold thickness, is calculated through time-depth conversion.

[0013] Preferably, according to the characteristics of the target core and logging curves, a detailed lateral stratigraphic comparison is performed on the target layer and its upper and lower strata to establish a sedimentary sequence framework.

[0014] Preferably, the sand-to-ground ratio and the mean porosity of sandstone and mudstone of each well section are extracted as follows:

[0015] In the well section between the corresponding point of the interface to be restored in each well and the corresponding point of the upper reference surface, the sand and mudstone thickness of each well is extracted respectively, the sand-to-ground ratio in each well section is calculated, and the mean porosity of the sand and mudstone in the well section is extracted respectively.

[0016] Preferably, the multi-well lithology difference compaction correction amount is estimated based on the average value of the sand-to-formation ratio and the porosity of sand-mudstone in each well section as follows: according to the core test and logging interpretation data, a prediction model relationship between the compaction correction amount and the sand-to-formation ratio and the porosity change rate of sand-mudstone is established, and the lithology difference compaction correction amount of each well is estimated based on this relationship. The relationship is specifically:

[0017]

[0018] Among them, ΔH is the correction amount of lithology difference compaction; H1 is the thickness of the impression; α is the sand-to-ground ratio within the impression thickness section; are the mean porosities of sand and mudstone within the thickness section of the impression, respectively; are the initial porosity values ​​of sandstone and mudstone deposits, respectively.

[0019] Preferably, the distribution of the lithology difference compaction correction amount is obtained according to the mold thickness and the multi-well lithology difference compaction correction amount as follows:

[0020] The seismic frequency-fractionated bodies of the upper abdominal strata and their corresponding plane attributes are extracted in the seismic time window corresponding to the mold thickness. Taking the compaction correction of multiple well points as the target, the seismic attributes are optimized and fused for analysis based on the neural network deep learning method. The lateral prediction relationship between the compaction correction of well points and multiple seismic plane attributes is established, and the distribution of compaction correction of lithology differences is obtained using the prediction relationship.

[0021] Preferably, the stamp thickness is corrected by embedding a plurality of lithology difference compaction correction amounts into the stamp thickness based on the lithology difference compaction correction amount distribution to obtain a corrected stamp thickness distribution.

[0022] A sandstone-mudstone stratum micro-paleomorphology characterization system, comprising:

[0023] Calculation module: determine the upper reference surface according to the overlying strata of the interface to be restored, and calculate the thickness of the mold;

[0024] Extraction module: According to the characteristics of the target core and well logging curves, the sedimentary sequence framework is established to extract the sand-to-ground ratio and the mean porosity of sandstone and mudstone in each well section;

[0025] Correction amount estimation module: estimates the compaction correction amount of multi-well lithology differences based on the sand-to-ground ratio of each well section and the average porosity of sandstone and mudstone;

[0026] Correction amount distribution acquisition module: obtains the lithology difference compaction correction amount distribution according to the mold thickness and the multi-well lithology difference compaction correction amount;

[0027] Paleomorphology characterization module: Correct the mold thickness based on the distribution of lithology difference compaction correction to obtain the corrected mold thickness distribution, convert the corrected mold thickness distribution to obtain the elevation plane of the interface to be restored.

[0028] A computer device comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of a method for characterizing micro-paleogeomorphology of sandstone and mudstone strata when executing the computer program.

[0029] A computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for characterizing micro-paleomorphology of sandstone and mudstone strata are implemented.

[0030] Compared with the prior art, the present invention has the following beneficial effects: The present invention provides a method for characterizing micro-paleomorphology of sandstone and mudstone formations, specifically a method and process for characterizing sedimentary micro-paleomorphology of sandstone and mudstone formations combining well and seismic data, and focuses on proposing a single-well estimation and plane prediction method for estimating the compaction correction amount of sandstone and mudstone lithology differences. This method can be used to better restore and accurately characterize the distribution morphology of the sedimentary micro-paleomorphology of the target layer based on core, well logging, seismic and other data, laying a methodological foundation for the restoration of sedimentary micro-paleomorphology in relatively stable structural areas. The finely characterized spatial morphology of sedimentary micro-paleomorphology can provide an important basis for the analysis of the distribution and temporal and spatial superposition relationship of sediment sources and local thin intersand bodies. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] Figure 1 This is a flow chart of a method for characterizing micro-paleomorphology of sandstone and mudstone strata based on seismic multi-attribute fusion according to an embodiment of the present invention;

[0032] Figure 2 For the present invention, a seismic profile of the paleo-geomorphic interface to be restored and the upper belly datum surface is picked up;

[0033] Figure 3 The sedimentary stratigraphic framework constructed by the present invention through fine lateral comparison of the stratigraphic layers;

[0034] Figure 4 An initial impression thickness map estimated for the present invention;

[0035] Figure 5 It is a plane diagram of compaction correction amount based on earthquake multi-attribute fusion prediction of the present invention;

[0036] Figure 6 This is a morphological diagram of the paleo-geomorphic interface after restoration and depiction of the present invention;

[0037] Figure 7 This is a flow chart of a method for depicting micro-paleomorphology of sandstone and mudstone strata according to the present invention;

[0038] Figure 8 The present invention is a block diagram of a sandstone-mudstone stratum micro-paleomorphology characterization system. DETAILED DESCRIPTION

[0039] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings here can be arranged and designed in various different configurations.

[0040] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention claimed for protection, but merely represents selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0041] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, further definition and explanation thereof is not required in subsequent drawings.

[0042] In the description of the embodiments of the present invention, it should be noted that if the terms "upper", "lower", "horizontal", "inner", etc. indicate an orientation or positional relationship based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship in which the product of the invention is usually placed when in use, it is only for the convenience of describing the present invention and simplifying the description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention. In addition, the terms "first", "second", etc. are only used to distinguish the description, and cannot be understood as indicating or implying relative importance.

[0043] In addition, if the term "horizontal" appears, it does not mean that the component must be absolutely horizontal, but can be slightly tilted. For example, "horizontal" only means that its direction is more horizontal than "vertical", which does not mean that the structure must be completely horizontal, but can be slightly tilted.

[0044] In the description of the embodiments of the present invention, it is also necessary to explain that, unless otherwise clearly specified and limited, the terms "set", "install", "connect", and "connect" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal connection of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0045] In order to enable those skilled in the art to better understand the technical solution of the present invention, the present invention will be further described in detail below with reference to the accompanying drawings.

[0046] like Figure 7 As shown, the present invention provides a method for depicting micro-paleomorphology of sandstone and mudstone strata, comprising:

[0047] S101 determines the upper reference surface according to the overlying stratum of the interface to be restored, and calculates the thickness of the mold.

[0048] A suitable stable interface is found in the overlying strata of the interface to be restored and determined as the upper datum plane. Specifically, the coinciding axes of the interface to be restored and the upper datum plane are picked up from the seismic section through fine stratigraphic calibration, and the thickness between the interface to be restored and the upper datum plane, i.e., the mold thickness, is calculated through time-depth conversion.

[0049] S102 establishes a sedimentary sequence framework based on the characteristics of the target core and logging curves, and extracts the sand-to-formation ratio and the average porosity of sandstone and mudstone in each well section.

[0050] It can be divided into the following two steps:

[0051] ① Combined with the characteristics of core and logging curves, a detailed lateral stratigraphic comparison is carried out between the target layer and its upper and lower strata to build a sedimentary stratigraphic framework.

[0052] ② In the well section between the corresponding point of the interface to be restored in each well and the corresponding point of the upper reference surface, the thickness of sand and mudstone in each well is extracted respectively, the sand-to-ground ratio in each well section is calculated, and the mean porosity of sand and mudstone in the well section is extracted respectively.

[0053] S103 estimates the multi-well lithology difference compaction correction amount based on the sand-to-formation ratio of each well section and the average porosity of sandstone and mudstone.

[0054] During the normal sedimentation and diagenesis of sandstone and mudstone, the sand and mud mineral components are affected differently by compaction operations, which is mainly reflected in the difference in the pore space changes of sand and mud mineral particles before and after compaction. Generally speaking, sandstone is not easy to compact, and the pore space changes are small, while mudstone is easy to compact, and the porosity space changes greatly. Therefore, although the total mineral porosity of the rock decreases overall before and after compaction, the porosity changes in different regions are different due to the difference in sand-mud ratio, and the porosity change has a certain relationship with the compaction coefficient, that is, the compaction coefficient or compaction correction amount is affected by regional lithology differences and manifests itself as corresponding differential distribution characteristics. Based on data such as core testing or well logging interpretation, a predictive model relationship between the compaction correction coefficient or compaction correction amount and the sand-to-ground ratio and the porosity change rate of sand-mud rock can be established. Based on this relationship, the lithology difference compaction correction amount of each well can be estimated.

[0055] S104 obtains the distribution of lithology difference compaction correction amount according to the stamp thickness and the multi-well lithology difference compaction correction amount.

[0056] The seismic frequency-dividing volume of the upper stratum and its corresponding plane attributes are extracted in the seismic time window corresponding to the thickness of the mold. Taking the compaction correction of multiple wells as the target, the seismic attributes are optimized and fused based on the neural network deep learning method, and the lateral prediction relationship between the compaction correction of wells and multiple seismic plane attributes is established. The plane distribution of the compaction correction of wells is obtained using the prediction relationship.

[0057] S105 corrects the mold thickness based on the distribution of the compaction correction amount of the lithology difference to obtain the corrected mold thickness distribution, converts the corrected mold thickness distribution to obtain the elevation plane of the interface to be restored.

[0058] The mold thickness is compacted and corrected, and the lithology difference compaction correction is embedded to obtain the corrected mold thickness distribution. The upper reference plane event axis of the interface to be restored is flattened, and the mold thickness after compaction correction is converted based on the mirror relationship of the mold method to obtain the elevation plane of the interface to be restored, thereby restoring and depicting the distribution of sedimentary micro-paleomorphology.

[0059] Example:

[0060] Using the sandstone-mudstone stratum micro-paleomorphology characterization method based on seismic multi-attribute fusion analysis ( Figure 1 ), taking the Cretaceous sandstone and mudstone strata in a certain area of ​​my country as an example, the paleo-geomorphology restoration and characterization of the sandstone and mudstone strata sedimentary structure was realized, including the following steps:

[0061] (1) For the Cretaceous strata in a certain area, the geological layers are first calibrated to the corresponding seismic event axes through fine horizon calibration. At the same time, a suitable stable interface is selected from the overlying strata of the interface to be restored on the seismic section to be determined as the upper datum plane, usually an interface with good event continuity. Then, the corresponding seismic horizons are picked up along the event axes between the interface to be restored and the upper datum plane ( Figure 2 ). Then, the stratum thickness between the interface to be restored and the upper reference surface is calculated by time-depth conversion, i.e., the mold thickness H1 ( Figure 3 ).

[0062] (2) Based on the characteristics of cores and well logging curves, a detailed lateral comparison of the Cretaceous sandstone and mudstone and adjacent strata was carried out to build a sedimentary stratigraphic framework and divide the sandstone and mudstone intervals ( Figure 4 ). In the well section between the corresponding point of the interface to be restored in each well and the corresponding point of the upper reference surface, the thickness of sand and mudstone in each well is extracted respectively, and the sand-to-ground ratio α in the well section is calculated; at the same time, the mean porosity of sand and mudstone in the well section is extracted respectively and

[0063]

[0064] (3) Based on the data of core test or logging interpretation, the compaction correction amount of each well is estimated by using the prediction model relationship (1) between the compaction correction amount and the sand-to-shale ratio and the porosity change rate of sand and mudstone. The specific calculation formula is:

[0065]

[0066] In the above formula: ΔH is the lithology difference compaction correction, m; H1 is the mold thickness, m; α is the sand-to-ground ratio within the mold thickness section, %; are the mean porosities of sand and mudstone within the thickness section of the impression, %; are the initial porosity values ​​of sandstone and mudstone deposition, respectively. Here, the initial porosity of sand and mudstone is obtained based on the fitting relationship between the porosity of sand and mudstone in each well and the depth. After fitting the porosity of sand and mudstone with depth in this area, the initial porosity of sand and mudstone is 41.3% and 56.75%, respectively.

[0067] (4) Based on the seismic post-stack data, frequency division processing is performed to obtain the seismic frequency division volume, and the corresponding seismic plane attributes are extracted in the seismic time window corresponding to the mold thickness. Taking the multi-well point compaction correction as the target, the seismic attributes are optimized and fused based on the neural network method to establish the lateral prediction relationship between the well point compaction correction and various seismic plane attributes, and the interpolation trend surface of the multi-well point compaction correction is used as a constraint to obtain the plane distribution of the well point compaction correction using the prediction relationship (2) ( Figure 5 ).

[0068] ΔH xy =f(X1,X2,X3,…,X i )

[0069] In the above formula, ΔH xy is the plane distribution value of the lithology difference compaction correction, m; X i These are earthquake attributes selected after deep learning optimization.

[0070] (5) Add the lithology difference compaction correction to the impression thickness to obtain the corrected impression thickness. Flatten the datum plane phase axis, transform the impression thickness after compaction correction based on the impression method mirror relationship, obtain the elevation plane of the interface to be restored, and further restore and depict the distribution of sedimentary micro-paleo-geomorphology. From the depicted sedimentary micro-paleo-geomorphology, the paleo-geomorphology of the Cretaceous sandstone and mudstone strata in the block is relatively clear. It can be seen that the terrain in the southeast of the block is low and gradually increases toward the northwest. It can be judged that the sediment source comes from the southeast, and the sediment gradually overlaps and expands toward the northwest ( Figure 6 ).

[0071] A terminal device is provided in one embodiment of the present invention. The terminal device of this embodiment includes: a processor, a memory, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps in the above-mentioned method embodiments are implemented. Alternatively, when the processor executes the computer program, the functions of the modules / units in the above-mentioned device embodiments are implemented.

[0072] The computer program may be divided into one or more modules / units, and the one or more modules / units are stored in the memory and executed by the processor to accomplish the present invention.

[0073] The terminal device may be a computing device such as a desktop computer, a notebook, a PDA, a cloud server, etc. The terminal device may include, but is not limited to, a processor and a memory.

[0074] The processor can be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc.

[0075] The memory may be used to store the computer program and / or module, and the processor implements various functions of the terminal device by running or executing the computer program and / or module stored in the memory and calling the data stored in the memory.

[0076] If the module / unit integrated in the terminal device is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present invention implements all or part of the processes in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium. It should be noted that the content contained in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable media do not include electric carrier signals and telecommunication signals.

[0077] Although the embodiments of the present invention are described above in conjunction with the accompanying drawings, the present invention is not limited to the above-mentioned specific embodiments and application fields, and the above-mentioned specific embodiments are only illustrative and instructive, rather than restrictive. Under the guidance of the specification, a person skilled in the art can also make many forms without departing from the scope of protection of the claims of the present invention, all of which belong to the scope of protection of the present invention.

Claims

1. A method for depicting micro-paleomorphology of sandstone and mudstone strata, characterized in that: include: Determine the upper reference surface according to the overlying strata of the interface to be restored, and calculate the thickness of the impression; According to the characteristics of the target core and well logging curves, a sedimentary sequence framework is established to extract the sand-to-formation ratio and the mean porosity of sandstone and mudstone in each well section; The multi-well lithology difference compaction correction amount is estimated based on the average value of sand-to-shale ratio and sand-mudstone porosity in each well section; The distribution of lithology difference compaction correction is obtained according to the stamp thickness and the multi-well lithology difference compaction correction; The mold thickness is corrected based on the distribution of compaction correction amount of lithology difference to obtain the corrected mold thickness distribution, and the corrected mold thickness distribution is converted to obtain the elevation plane of the interface to be restored.

2. The method for depicting micro-paleomorphology of sandstone and mudstone strata according to claim 1, characterized in that: The upper datum plane is determined according to the overlying strata of the interface to be restored, and the thickness of the mold is calculated as follows: a suitable stable interface is found in the overlying strata of the interface to be restored to be determined as the upper datum plane. Specifically, the coinciding axes of the interface to be restored and the upper datum plane are picked up from the seismic section through fine stratigraphic calibration, and the thickness between the interface to be restored and the upper datum plane, i.e., the mold thickness, is calculated through time-depth conversion.

3. The method for depicting micro-paleomorphology of sandstone and mudstone strata according to claim 1, characterized in that: According to the characteristics of the target core and logging curves, a detailed lateral stratigraphic comparison is carried out between the target layer and its upper and lower strata to establish a sedimentary sequence framework.

4. The method for depicting micro-paleomorphology of sandstone and mudstone strata according to claim 1, characterized in that: The specific values ​​of the sand-to-soil ratio and the porosity of sandstone and mudstone in each well section are: In the well section between the corresponding point of the interface to be restored in each well and the corresponding point of the upper reference surface, the sand and mudstone thickness of each well is extracted respectively, the sand-to-ground ratio in each well section is calculated, and the mean porosity of the sand and mudstone in the well section is extracted respectively.

5. The method for depicting micro-paleomorphology of sandstone and mudstone strata according to claim 1, characterized in that: The compaction correction amount of lithology difference of multiple wells is estimated based on the average value of sand-to-formation ratio and sand-mudstone porosity of each well section. According to the core test and logging interpretation data, a prediction model relationship between the compaction correction amount and the sand-to-formation ratio and the porosity change rate of sand-mudstone is established. The compaction correction amount of lithology difference of each well is estimated based on this relationship. The specific relationship is: Among them, ΔH is the correction amount of lithology difference compaction; H1 is the thickness of the impression; α is the sand-to-ground ratio within the impression thickness section; are the mean porosities of sand and mudstone within the thickness section of the impression, respectively; are the initial porosity values ​​of sandstone and mudstone deposits, respectively.

6. The method for depicting micro-paleomorphology of sandstone and mudstone strata according to claim 1, characterized in that: The distribution of lithology difference compaction correction amount obtained according to the stamp thickness and multi-well lithology difference compaction correction amount is as follows: The seismic frequency-fractionated bodies of the upper abdominal strata and their corresponding plane attributes are extracted in the seismic time window corresponding to the mold thickness. Taking the compaction correction of multiple well points as the target, the seismic attributes are optimized and fused for analysis based on the neural network deep learning method. The lateral prediction relationship between the compaction correction of well points and multiple seismic plane attributes is established, and the distribution of compaction correction of lithology differences is obtained using the prediction relationship.

7. The method for depicting micro-paleomorphology of sandstone and mudstone strata according to claim 1, characterized in that: The mold thickness is corrected specifically by embedding a plurality of lithology difference compaction correction amounts into the mold thickness based on the lithology difference compaction correction amount distribution to obtain a corrected mold thickness distribution.

8. A micro-paleomorphology characterization system for sandstone and mudstone strata, characterized in that: include: Calculation module: determine the upper reference surface according to the overlying strata of the interface to be restored, and calculate the thickness of the mold; Extraction module: According to the characteristics of the target core and well logging curves, the sedimentary sequence framework is established to extract the sand-to-ground ratio and the mean porosity of sandstone and mudstone in each well section; Correction amount estimation module: estimates the compaction correction amount of multi-well lithology differences based on the sand-to-ground ratio of each well section and the average porosity of sandstone and mudstone; Correction amount distribution acquisition module: obtains the lithology difference compaction correction amount distribution according to the mold thickness and the multi-well lithology difference compaction correction amount; Paleomorphology characterization module: Correct the mold thickness based on the distribution of lithology difference compaction correction to obtain the corrected mold thickness distribution, convert the corrected mold thickness distribution to obtain the elevation plane of the interface to be restored.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the method for characterizing micro-paleomorphology of sandstone and mudstone strata as described in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of a method for characterizing micro-paleomorphology of sandstone and mudstone strata as claimed in any one of claims 1 to 7 are implemented.

Images