Quantitative analysis method for transport path of source-sink system based on geological body carving technology

Abstract

This invention discloses a quantitative analysis method for transport channels in source-sink systems based on geological body sculpting technology, comprising the following steps: S1: acquiring seismic data, well logging data, well logging data, and core data; S2: geophysical comprehensive interpretation; S3: seismic attribute extraction, optimization, and fusion; S4: establishing a source-sink system transport system. The beneficial effects of this invention are: This invention is a fine analysis technology for transport channels based on geological body sculpting, combined with paleogeographic reconstruction maps, which can fuse different attributes to identify and depict transport channels, achieving separation of transport channels from background attributes and clearly depicting the spatial distribution of transport channels. Compared with existing technologies, the data used in this invention are all commonly used data, with wide applicability and easy access, lower cost, convenient and efficient operation, and higher accuracy in transport channel identification than conventional layer attribute extraction techniques.

Description

Technical Field

[0001] This invention relates to the field of oil and gas geological exploration, and in particular to a quantitative analysis method for the transport channels of source-sink systems based on geological body carving technology. Background Technology

[0002] Source-sink system analysis is a crucial step in petroleum exploration, primarily studying the source (source), transport channels (channels), and final deposition sites (sinks) of sediments. Traditional source-sink system analysis methods mainly rely on geological studies, such as field investigations and drilling data analysis; however, due to insufficient data, their analytical accuracy and reliability are somewhat limited.

[0003] In existing technologies, some methods attempt to introduce new technical means. For example, the invention patent CN104699977A, "Method for Quantitative Analysis of Source Material Using Light and Heavy Minerals," determines the source direction by screening and organizing data, normalizing light and heavy mineral data, establishing a source area discrimination formula, and quantitatively analyzing the relationship between mineral stability indices. However, this method requires sample data of light and heavy minerals from the source area and sedimentary area as a basis, making it unsuitable for well-less or few-well conditions, and the testing cost is high. Furthermore, some regions lack a systematic foundation for heavy mineral analysis. Another example is the invention patent CN115639597A, "A Fine Characterization Method of Transportation Systems Based on Seismic Attributes," which extracts the response characteristics of the transportation system through a combination of point, line, and surface analysis of three-dimensional paleogeography, two-dimensional transportation channel systems, and sedimentary facies, and quantitatively analyzes source-sink system elements to achieve fine characterization. However, this method has limited application of seismic data in source analysis, and single seismic attributes have multiple interpretations, making it difficult to reflect the integrity of the source-sink system. Overall, existing methods have technical bottlenecks in identifying transport channels and analyzing material sources. Summary of the Invention

[0004] To address the limitations of existing technologies in source analysis and the ambiguity of single seismic attributes, which make it difficult to reflect the overall source-sink system, this invention provides a quantitative analysis method for the transport channels of a source-sink system based on geological body sculpting technology.

[0005] To achieve the above-mentioned objectives, this invention provides a quantitative analysis method for the transport channels of a source-sink system based on geological body carving technology, the method comprising the following steps: S1: Acquire seismic data, well logging data, well logging data, and core data; S2: Geophysical Comprehensive Interpretation; S3: Seismic attribute extraction, optimization, and fusion; S4: Establish a source-to-consumer system for data transfer.

[0006] The S2 geophysical comprehensive interpretation includes: establishing a sequence stratigraphic framework and paleogeographic reconstruction.

[0007] The establishment of the sequence stratigraphic framework specifically involves: firstly, meticulous tracing of geological interfaces, followed by stratigraphic and fault interpretation. Stratigraphic interpretation includes tracing, closing, and densifying the interfaces of key tectonic phases, while fault interpretation includes determining the strata crossed by faults and calculating fault displacements. Based on the results of stratigraphic and fault interpretation, a sequence stratigraphic framework is constructed.

[0008] The paleogeographic restoration process involves calculating the apparent thickness of strata using well and seismic data, obtaining the true thickness of strata through dip correction, and then performing erosion recovery and compaction correction on the true thickness to obtain the original sedimentary thickness of the strata, thus achieving paleogeographic restoration. Specifically, it includes the following steps: 1) Perform stratigraphic correction on the study area. After obtaining the apparent thickness through drilling data, perform stratigraphic correction and dip correction on the target strata, and calculate the true thickness of the target strata. 2) The erosion thickness is restored by using the tectonic trend method. Combined with the lateral variation law of the stratigraphic thickness, the thickness of the uneroded stratigraphic layer and the sedimentary boundary, the plane contour map of the stratigraphic erosion thickness is obtained by interpolation or extrapolation. 3) Compaction correction was performed. Based on the relationship between the compaction coefficient and the burial depth, the compaction curve of the study area was established using the formula k=b-alnH, where: k is the compaction coefficient, a and b are constants, and H is the burial depth. For mudstone or shale, a=0.1 and b=1.46; for sandstone, a=0.08 and b=1.37. 4) Based on different depositional periods, establish paleogeographic feature models and compare them with paleogeographic models obtained from inversion simulations. If the two models are consistent within the same depositional period, the paleogeography reflected by the residual strata is confirmed to be accurate. If there are significant differences, the fault interpretation and stratigraphic interpretation are corrected, or the density of two-dimensional profile interpretation is increased to restore the paleogeography and determine the difference between the paleodepositional boundary of the depression and the present tectonic boundary.

[0009] Specifically, S3 is: S301: Extract seismic attributes based on seismic data, including amplitude, frequency, and phase attributes; S302: Optimize the extracted seismic single attributes; specifically: S3021: Analyze the correlation between seismic attributes and reservoir lithology, and select the attributes most sensitive to geological conditions; S3022: Perform correlation analysis on the selected attributes, and retain only the attribute with the best performance among the same type; S3023: To address the issues of different dimensions of earthquake attributes and multicollinearity among multiple attributes, normalization and principal component analysis are performed to optimize earthquake attributes. S303: A fusion method based on sedimentary model constraints is used to fuse the optimized single seismic attributes. During the fusion process, the influence of the attributes on the sedimentary model and their correlation and complementarity are considered. The fused seismic attributes are generated through mathematical operations and model optimization, providing an accurate data foundation for subsequent transport channel analysis.

[0010] S4 specifically involves: finely characterizing and identifying the transport channel, and constructing the source-sink system transport architecture; S4 specifically includes the following steps: S401: Based on seismic multi-attribute fusion technology, on the restored paleogeographic model, the transport path, the distribution orientation of transport channels, the combination of different transport channels and their correspondence with the source area and sedimentary area are qualitatively characterized. S402: Based on the qualitative characterization results, construct a survey profile and measure the transverse width, longitudinal depth, width-to-depth ratio, cross-sectional area, and extension distance of different transport channels on the survey profile; S403: Based on the catchment area range divided at different times, the average value of the lateral width, longitudinal depth, width-to-depth ratio, and extension distance of all transport channels within the same catchment area is taken to form the average value of the characteristic parameters of the transport channels distinguished by the catchment area.

[0011] The beneficial effects of this invention are as follows: This invention is based on a fine analysis technology for transport channels sculpted from geological bodies. Combined with paleogeographic reconstruction maps, it can integrate different attributes to identify and delineate transport channels, achieving separation of transport channels from background attributes and clearly depicting the spatial distribution of transport channels. Compared with existing technologies, the data used in this invention are all commonly used data, with wide applicability and easy availability. It is also lower in cost, more convenient and efficient in operation, and has higher accuracy in transport channel identification than conventional layer-by-layer attribute extraction techniques. Attached Figure Description

[0012] Figure 1 This is a flowchart of Embodiment 2 of the present invention.

[0013] Figure 2 This refers to the seismic sequence grid in Embodiment 3 of the present invention.

[0014] Figure 3 This is a diagram showing the results of ancient landform restoration in Embodiment 3 of the present invention.

[0015] Figure 4(a) is a diagram showing the root mean square amplitude properties of the lower sub-segment of the sand in Embodiment 3 of the present invention.

[0016] Figure 4(b) is a diagram showing the average instantaneous frequency attribute of the lower subsegment of the sand in Embodiment 3 of the present invention.

[0017] Figure 5 This is a diagram showing the effect of earthquake attribute fusion in Embodiment 3 of the present invention.

[0018] Figure 6 This is a diagram showing the result of the geological body carving and transportation channel in Embodiment 3 of the present invention.

[0019] Figure 7(a) is a cross-sectional view of the ancient valley at the Gaoqing Fault in the Boxing Depression in Embodiment 3 of the present invention.

[0020] Figure 7(b) is a cross-sectional view of the ancient valley at the Luxi Uplift in Embodiment 3 of the present invention.

[0021] Figure 8(a) is a cross-sectional view of the fracture groove at the broom-shaped fracture zone in Shicun, Boxing Depression, in Embodiment 3 of the present invention.

[0022] Figure 8(b) is a cross-sectional view of the fracture groove at the Jinjia fracture site in Embodiment 3 of the present invention.

[0023] Figure 9 This is a planar distribution diagram of the transport channels in the Boxing Depression in Embodiment 3 of the present invention.

[0024] Figure 10 Examples of ancient valley classification styles and cross-sectional views in Embodiment 2 of the present invention.

[0025] Figure 11 This is an example diagram showing the various categories and cross-sectional views of the interrupt slots in Embodiment 2 of the present invention. Detailed Implementation

[0026] To clearly illustrate the technical features of this solution, the following detailed implementation method will be used to explain the solution.

[0027] Example 1 This invention provides a quantitative analysis method for the transport channels of a source-sink system based on geological body carving technology. The method includes the following steps: S1: Acquire seismic data, well logging data, well logging data, and core data; S2: Geophysical comprehensive interpretation; wherein, the S2 geophysical comprehensive interpretation includes: establishing a sequence stratigraphic framework and paleogeographic reconstruction.

[0028] The establishment of the sequence stratigraphic framework specifically involves: firstly, meticulous tracing of geological interfaces, followed by stratigraphic and fault interpretation. Stratigraphic interpretation includes tracing, closing, and densifying the interfaces of key tectonic phases, while fault interpretation includes determining the strata crossed by faults and calculating fault displacements. Based on the results of stratigraphic and fault interpretation, a sequence stratigraphic framework is constructed.

[0029] The paleogeomorphological restoration involves calculating the apparent thickness of strata using well and seismic data, obtaining the true thickness of strata through dip correction, and then performing erosion and compaction corrections on the true thickness to obtain the original sedimentary thickness of the strata, thus achieving paleogeomorphological restoration. The paleogeomorphological restoration specifically includes the following steps: 1) Perform stratigraphic correction on the study area. After obtaining the apparent thickness through drilling data, perform stratigraphic correction and dip correction on the target strata, and calculate the true thickness of the target strata. 2) The erosion thickness is restored by using the tectonic trend method. Combined with the lateral variation law of the stratigraphic thickness, the thickness of the uneroded stratigraphic layer and the sedimentary boundary, the plane contour map of the stratigraphic erosion thickness is obtained by interpolation or extrapolation. 3) Compaction correction was performed. Based on the relationship between the compaction coefficient and the burial depth, the compaction curve of the study area was established using the formula k=b-alnH, where: k is the compaction coefficient, a and b are constants, and H is the burial depth. For mudstone or shale, a=0.1 and b=1.46; for sandstone, a=0.08 and b=1.37. 4) Based on different depositional periods, establish paleogeographic feature models and compare them with paleogeographic models obtained from inversion simulations. If the two models are consistent within the same depositional period, the paleogeography reflected by the residual strata is confirmed to be accurate. If there are significant differences, the fault interpretation and stratigraphic interpretation are corrected, or the density of two-dimensional profile interpretation is increased to restore the paleogeography and determine the difference between the paleodepositional boundary of the depression and the present tectonic boundary.

[0030] S3: Seismic attribute extraction, optimization, and fusion; specifically: S301: Extract seismic attributes based on seismic data, including parameters such as amplitude, frequency, and phase; S302: Optimize the extracted seismic single attributes; specifically: S3021: Analyze the correlation between seismic attributes and reservoir lithology, and select the attributes most sensitive to geological conditions; S3022: Perform correlation analysis on the selected attributes, and retain only the attribute with the best performance among the same type; S3023: To address the issues of different dimensions of earthquake attributes and multicollinearity among multiple attributes, normalization and principal component analysis are performed to optimize earthquake attributes. S303: A fusion method based on sedimentary model constraints is used to fuse the optimized single seismic attributes. During the fusion process, the influence of the attributes on the sedimentary model and their correlation and complementarity are considered. The fused seismic attributes are generated through mathematical operations and model optimization, providing an accurate data foundation for subsequent transport channel analysis.

[0031] S4: Establish a source-to-consumer system transportation framework; specifically, this involves: precisely characterizing and identifying transportation channels, and constructing a source-to-consumer system transportation framework; S4 specifically includes the following steps: S401: Based on seismic multi-attribute fusion technology, on the restored paleogeographic model, the transport path, the distribution orientation of transport channels, the combination of different transport channels and their correspondence with the source area and sedimentary area are qualitatively characterized. S402: Based on the qualitative characterization results, construct a survey profile and measure the transverse width, longitudinal depth, width-to-depth ratio, cross-sectional area, and extension distance of different transport channels on the survey profile; S403: Based on the catchment area range divided at different times, the average value of the lateral width, longitudinal depth, width-to-depth ratio, and extension distance of all transport channels within the same catchment area is taken to form the average value of the characteristic parameters of the transport channels distinguished by the catchment area.

[0032] Example 2 This invention provides a quantitative analysis method for the transport channels of source-sink systems based on geological body sculpting technology. It enhances the application of seismic data in the transport system of source-sink systems, establishes transport combination relationships, and quantitatively characterizes transport channels to improve the accuracy and reliability of the analysis. Furthermore, it eliminates the need for costly experimental testing and strengthens the verification of source direction. Specifically, the analysis method includes the following steps: Step 1: Comprehensive and detailed geophysical interpretation The process of detailed geophysical interpretation mainly includes two steps: the establishment of sequence stratigraphic framework and the reconstruction of paleogeography.

[0033] (1) Establishment of sequence stratigraphy framework Detailed tracing of geological interfaces is fundamental to paleogeographic characterization. The precision and accuracy of paleogeographic characterization are limited by stratigraphic tracing. Firstly, based on seismic, well logging, and core data, detailed stratigraphic and fault structural interpretations are conducted for the study area. Stratigraphic interpretation includes tracing, closing, and densifying key tectonic phase interfaces; fault structural interpretation includes identifying the strata crossed by faults and determining fault displacement. Based on this, a sequence stratigraphic framework is established.

[0034] (2) Paleomorphological restoration The paleogeomorphological restoration process involves calculating the apparent thickness of strata using well and seismic data, obtaining the true thickness of strata through dip correction, and then performing erosion recovery and compaction correction on the true thickness to obtain the original sedimentary thickness of the strata, thus achieving paleogeomorphological restoration. The basic workflow of paleogeomorphological restoration is as follows: ① Perform formation correction on the study area. The thickness obtained from drilling data is often the apparent thickness. Perform formation correction and dip correction on the target layer to determine the true thickness of the formation. ② The erosion thickness was restored using the tectonic trend method. Combined with the lateral variation law of the stratigraphic thickness, the thickness of the uneroded stratigraphic layers and the sedimentary boundary, the plane contour map of the stratigraphic erosion thickness was obtained by interpolation or extrapolation. ③ Perform compaction correction using empirical formulas: (Where k is the compaction coefficient, a and b are constants, and H is the burial depth; for mudstone and shale, a=0.1, b=1.46; for sandstone, a=0.08, b=1.37) To extend this, the compaction curve of the study area is established, and the relationship between compaction coefficient and depth is fitted. ④ Establish paleogeographic feature models for different sedimentary periods and compare them with paleogeographic models obtained from inversion simulation. If the paleogeographic morphology shows consistency within the same sedimentary period, then the paleogeography reflected by the obtained residual strata during the sedimentary period is accurate. Conversely, if there are significant differences between the two, it is necessary to make necessary corrections to the interpretation of faults and stratigraphic positions, or increase the density of two-dimensional profile interpretation to improve the accuracy and reliability of the interpretation, reconstruct the paleogeography, and determine the differences between the paleodepositional boundaries and the present tectonic boundaries of each depression.

[0035] Step 2: Seismic Multi-Attribute Fusion (1) Optimization of seismic attributes First, seismic attributes are extracted based on seismic data, well logging data, well logging data, and core data. Then, the individual seismic attributes are optimized, and the correlation between the attributes and reservoir lithology is analyzed to select the attributes most sensitive to geological conditions. Second, correlation analysis is performed among the selected attributes, and only the best-performing attribute of the same type is selected. Finally, due to the different dimensions of seismic attributes and the multicollinearity problem among multiple attributes, normalization and principal component analysis are performed before correlation analysis to optimize the attributes.

[0036] (2) Seismic attribute fusion The selected individual seismic attributes were fused using a seismic attribute fusion method based on sedimentary model constraints. During the fusion process, the influence of different attributes on the sedimentary model, as well as their correlations and complementarities, were considered. Through mathematical operations and model optimization, the fused seismic attributes were obtained. These attributes can more comprehensively reflect the subsurface sedimentary model, providing a more accurate and comprehensive data foundation for subsequent transport channel analysis.

[0037] Step 3, Geological Sculpting Techniques Based on the seismic attribute fusion, methods such as planar attribute analysis, profile analysis, and local planar attribute extraction are comprehensively utilized to carry out fine-grained sculpting of transport channels. Fine-grained sculpting of transport channels can provide accurate profile, planar, and local spatial distribution characteristics for quantitative analysis of transport channels.

[0038] Step 4: Establishment of the Yuanhui System Transportation Framework The study of the sediment transport system is crucial in the source-sink system analysis. By conducting detailed analysis of transport channels, we can clarify their types, distribution directions, combinations of different transport channels, and their correspondence with the source and sedimentary areas. Seismic data can effectively identify transport channels.

[0039] The different types, sizes, layouts, and specific combinations of transport channels play a decisive role in the direction and quantity of material transported from the source region to the basin, thus influencing the composition and development of the source-sink system. During transport, not all debris ultimately reaches the sedimentary basin; some may be deposited at the end of the transport path, while others may stall en route. Only the remaining debris continues to move with paleochannels and eventually settles within the basin. In common terrestrial rift basins, based on the type and size of sediment transport channels, they can be classified into two types: paleogullies and paleofault channels. (1) Ancient valley See Figure 10 Ancient gullies are landforms formed when the base level drops below the inflection point of the slope, causing erosion to begin on the surface. They typically develop along stratigraphic boundaries and are more common around basins or uplift peripheries. As the base level lowers, surface water begins to erode the land, leading to the formation of gullies. These gullies extend along the stratigraphic boundaries, deepening and expanding further in softer areas. The transport channels of these gullies are thin at both ends and thicker in the middle, appearing as a concave topography on seismic profiles and prone to seismic reflection axis faulting.

[0040] (2) Ancient fault groove The formation of fault channels is related to geological structures and sedimentary environments. In terrestrial lacustrine basins, crustal movement leads to the formation of faults, and the presence of faults has a significant impact on the transport and deposition of sediments. When sediments pass through fault channels, the lower elevation in the middle of the fault causes the sediments to be transported rapidly along the path of the channel and deposited at or near the end of the channel.

[0041] See Figure 11Fault-induced transport channels mainly include two types: single fault channels and double fault channels. Single fault channels typically develop along a single fault. The activity of this fault causes the tectonic block on one side to subside relative to the other, forming a sediment transport and deposition center. The subsidence of the tectonic block on one side of the fault forms a slope, and sediments move along this slope under gravity, resulting in greater sediment accumulation on the lowered side. Therefore, single fault channels are characterized by one side being steeper and thicker, and the other side being gentler and thinner, showing symmetrical graben structures and wedge-shaped filling seismic facies features on seismic profiles. Double fault channels, on the other hand, are channel structures formed under the control of two sets of relatively dipping faults. When the two sets of faults dip in opposite directions, they form a relatively low-lying area in the middle. This low-lying area becomes the main collection and transport channel for sediments. Double fault channels typically exhibit a symmetrical structure with relatively parallel edges on both sides, and parallel or subparallel seismic reflection characteristics inside. This indicates that sediments are deposited at a relatively continuous and consistent rate within the double fault channel, indicating a relatively stable depositional environment without severe external disturbance.

[0042] The basic process for establishing the Yuanhui system's data transfer system is as follows: ①Based on the fusion of multiple earthquake attributes, the transport paths, distribution orientations of transport channels, combination patterns between different transport channels, and their correspondence with source areas and sedimentary areas are depicted on the restored paleogeographic model, and the catchment area range is defined according to different periods. ② Based on the qualitative characterization results of the transport channels, a seismic profile is constructed, and the transport channels are quantitatively characterized. The quantitative characterization parameters include the width, depth, width-to-depth ratio, cross-sectional area, and extension distance of the transport channels. The quantitative characterization steps are as follows: Measuring the width, depth, width-to-depth ratio, and cross-sectional area of ​​different transport channels on the seismic profile; calculating the lateral width, longitudinal depth, width-to-depth ratio, cross-sectional area, and extension distance for each transport channel; and averaging the width, depth, width-to-depth ratio, and extension distance of all transport channels within the same catchment area, based on the catchment area divisions at different times, to form average values ​​of the transport channel characteristic parameters distinguished by the catchment area.

[0043] Example 3 This embodiment takes the Boxing Depression as an example to conduct a quantitative analysis of the source-sink system transport channels of the oilfield based on seismic multi-attribute fusion. The analysis method is the one provided in Embodiment 2. Specifically: I. Based on step 1 of Example 2, a seismic sequence framework is established, and stratigraphic correlation is performed to divide the seismic sequence framework of the Boxing Depression: The bottom of the Kongdian Formation: Tr - a regional unconformity. After being eroded by tectonic uplift in the late Mesozoic, the Kongdian Formation began to receive sediments. Top of Kong 2 section: T8- is in conformable contact with Kong 1 section, and the wave group at the interface has strong amplitude continuous reflection characteristics; the strata of Kong 2 section have a small distribution range, and the southern Kong 1 section is in unconformable contact with the underlying strata; The top of the Kongdian Formation: EK- has a local unconformable contact with the underlying ES4x strata, and a conformable or pseudo-conformable contact towards the center of the depression; the overall amplitude is medium to strong. Sand 4 Lower Top: Locally there are 3 unconformities, with medium to strong amplitude characteristics.

[0044] Based on the stratigraphic framework, paleomorphological reconstruction is completed according to step 1, such as... Figure 3 As shown.

[0045] II. Sensitive attribute optimization (as shown in Figure 4) and seismic attribute fusion (as described in step 2 of Example 2) are performed based on the method described in step 2 of Example 2. Figure 5 (As shown).

[0046] III. Following the method described in step 3 of Example 2, complete the carving of the transport channel using geological body carving technology, such as... Figure 6 As shown.

[0047] IV. Utilizing the morphological and paleogeographic features reflected in step 4 of Example 2 based on the seismic attribute fusion results, the combination relationships of transport channels were analyzed, and the transport channels were quantitatively characterized (Table 1). A statistical table of transport channel pattern types in the Dongying Depression was established (as shown in Figures 7 and 8). Figure 9 (As shown).

[0048] Table 1. Statistical Table of Transportation Channel Parameters in Boxing Depression

[0049] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A quantitative analysis method for the transport channels of a source-sink system based on geological body carving technology, characterized in that, The method includes the following steps: S1: Acquire seismic data, well logging data, well logging data, and core data; S2: Geophysical Comprehensive Interpretation; S3: Seismic attribute extraction, optimization, and fusion; S4: Establish a source-to-consumer system for data transfer.

2. The method according to claim 1, characterized in that, The S2 geophysical comprehensive interpretation includes: establishing a sequence stratigraphic framework and paleogeographic reconstruction.

3. The method according to claim 2, characterized in that, The establishment of the sequence stratigraphic framework specifically involves: firstly, meticulous tracing of geological interfaces, followed by stratigraphic and fault interpretation. Stratigraphic interpretation includes tracing, closing, and densifying the interfaces of key tectonic phases, while fault interpretation includes determining the strata crossed by faults and calculating fault displacements. Based on the results of stratigraphic and fault interpretation, a sequence stratigraphic framework is constructed.

4. The method according to claim 2, characterized in that, The paleogeographic restoration process involves calculating the apparent thickness of strata using well and seismic data, obtaining the true thickness of strata through dip correction, restoring the erosion amount and correcting for compaction on the true thickness of strata, and obtaining the original sedimentary thickness of the strata to achieve paleogeographic restoration.

5. The method according to claim 4, characterized in that, The restoration of ancient landforms specifically includes the following steps: 1) Perform stratigraphic correction on the study area. After obtaining the apparent thickness through drilling data, perform stratigraphic correction and dip correction on the target strata, and calculate the true thickness of the target strata. 2) The erosion thickness is restored by using the tectonic trend method. Combined with the lateral variation law of the stratigraphic thickness, the thickness of the uneroded stratigraphic layer and the sedimentary boundary, the plane contour map of the stratigraphic erosion thickness is obtained by interpolation or extrapolation. 3) Compaction correction was performed. Based on the relationship between the compaction coefficient and the burial depth, the compaction curve of the study area was established using the formula k=b-alnH, where: k is the compaction coefficient, a and b are constants, and H is the burial depth; 4) Based on different depositional periods, establish paleogeographic feature models and compare them with paleogeographic models obtained from inversion simulations. If the two models are consistent within the same depositional period, the paleogeography reflected by the residual strata is confirmed to be accurate. If there are significant differences, the fault interpretation and stratigraphic interpretation are corrected, or the density of two-dimensional profile interpretation is increased to restore the paleogeography and determine the difference between the paleodepositional boundary of the depression and the present tectonic boundary.

6. The method according to claim 1, characterized in that, Specifically, S3 is: S301: Extracting seismic attributes based on seismic data; S302: Optimize the extracted seismic single attributes; S303: A fusion method based on sedimentary model constraints is used to fuse the optimized single seismic attributes. During the fusion process, the influence of the attributes on the sedimentary model and their correlation and complementarity are considered. The fused seismic attributes are generated through mathematical operations and model optimization, providing an accurate data foundation for subsequent transport channel analysis.

7. The method according to claim 6, characterized in that, The preferred embodiment of S302 is as follows: S3021: Analyze the correlation between seismic attributes and reservoir lithology, and select the attributes most sensitive to geological conditions; S3022: Perform correlation analysis on the selected attributes, and retain only the attribute with the best performance among the same type; S3023: To address the issues of different dimensions of earthquake attributes and multicollinearity among multiple attributes, normalization and principal component analysis are performed to optimize earthquake attributes.

8. The method according to claim 1, characterized in that, Specifically, S4 involves: meticulously depicting and identifying the transport channel, and constructing a source-sink system transport architecture.

9. The method according to claim 8, characterized in that, S4 includes the following steps: S401: Based on seismic multi-attribute fusion technology, on the restored paleogeographic model, the transport path, the distribution orientation of transport channels, the combination of different transport channels and their correspondence with the source area and sedimentary area are qualitatively characterized. S402: Based on the qualitative characterization results, construct a survey profile and measure the transverse width, longitudinal depth, width-to-depth ratio, cross-sectional area, and extension distance of different transport channels on the survey profile; S403: Based on the catchment area range divided at different times, the average value of the lateral width, longitudinal depth, width-to-depth ratio, and extension distance of all transport channels within the same catchment area is taken to form the average value of the characteristic parameters of the transport channels distinguished by the catchment area.

Images