Method and device for recovering denudation thickness and thermal historical evolution process of superimposed basin

Abstract

The invention provides a method and a device for recovering the denudation thickness and the thermal historical evolution process of a superimposed basin. The method for recovering the denudation thickness and the thermal historical evolution process of the superimposed basin comprises the following steps: selecting a rock sample of a target work area; wherein the rock sample is at least one of a sandstone sample, a volcanic rock sample and a metamorphic rock sample containing apatite and zircon mineral particles; performing in-situ microcell isotope dating measurement on the rock sample to generate a dating measurement result; and determining the denudation thickness according to the ground temperature gradient of the target work area and the dating measurement result, and recovering the thermal historical evolution process. According to the method, quantitative recovery of the stratum denudation thickness caused by multi-stage tectonic cycle superposition and inversion of the dynamic evolution process of the thermal history under multi-stage tectonic activity superposition can be achieved, the defects in the prior art are effectively overcome, and the method has great significance in research and exploration evaluation of oil and gas accumulation in the complex superimposed basin.

Description

Technical Field

[0001] The present application belongs to the field of oil and gas exploration technology, in particular to the field of restoration technology of tectonic history and thermal history of sedimentary basins, and specifically to a method and device for restoring the erosion thickness and thermal history evolution process of a superimposed basin. Background Art

[0002] Superimposed basins are basins with complex structures that have undergone multiple phases of tectonic movement and are formed by the vertical superposition of multiple single-type basins (different tectonic layers). The uplift and erosion of strata caused by multiple phases of tectonic activity have a significant impact on the sedimentary history, tectonic evolution history, and thermal history of the entire basin. The thickness, age, and time interval of erosion are closely related to the formation, evolution, migration, accumulation, and preservation of oil and gas in superimposed basins. Therefore, the quantitative restoration of stratum erosion time and thickness is of great significance for the study of the thermal evolution history, burial history, hydrocarbon generation and expulsion history, and migration and accumulation history of superimposed basins, and even for oil and gas exploration and evaluation.

[0003] Most existing methods for quantitatively restoring the time and thickness of stratum erosion and the thermal evolution history of superimposed basins achieve this by establishing a fitting relationship. However, the establishment of such a fitting mathematical relationship is often non-unique, resulting in multiple solutions for the erosion thickness value and large errors. In addition, many parameter values ​​involved in the existing methods have large uncertainties, making it difficult to obtain accurate erosion thickness. In addition, some existing methods often require a large amount of drilling, logging, and core analysis data, and cannot be used in medium and low exploration areas. Existing methods also have their own applicable conditions and limitations. For example, the vitrinite reflectance method requires a large amount of test data and is significantly affected by thermal events. The stratigraphic comparison method based on seismic interpretation must be based on a small erosion range in the study area and is not applicable to areas where the erosion range covers the entire area. At the same time, existing methods cannot provide the thermal evolution process of superimposed basins under the superposition of multiple tectonic activities. Summary of the Invention

[0004] One objective of the present invention is to provide a method for recovering the erosion thickness and thermal history evolution of superimposed basins. This method can quantitatively recover the stratum erosion thickness caused by the superposition of multiple tectonic cycles and invert the dynamic thermal history evolution under the superposition of multiple tectonic activity periods. This method effectively overcomes the shortcomings of existing technologies and is of great significance for the research and exploration evaluation of oil and gas accumulation in complex superimposed basins.

[0005] Another object of the present invention is to provide a device for recovering the erosion thickness and thermal history of a superimposed basin. Another object of the present invention is to provide an electronic device comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of the method for recovering the erosion thickness and thermal history of a superimposed basin are implemented. Another object of the present invention is to provide a readable medium storing a computer program, and when the processor executes the computer program, the steps of the method for recovering the erosion thickness and thermal history of a superimposed basin are implemented.

[0006] In order to solve the technical problems in the background technology of this application, the present invention provides the following technical solutions:

[0007] In a first aspect, the present invention provides a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin, comprising:

[0008] Selecting a rock sample from the target work area; wherein the rock sample is at least one of a sandstone sample containing apatite and zircon mineral particles, a volcanic rock sample, and a metamorphic rock sample;

[0009] performing in-situ micro-area isotope dating on the rock sample to generate a dating result;

[0010] The erosion thickness is determined based on the geothermal gradient of the target work area and the dating measurement results, and the thermal history evolution process is restored.

[0011] In some embodiments of the present invention, selecting rock samples from a target work area includes:

[0012] The rock sample is selected from an overlying stratum or an underlying stratum of the eroded stratum in the target work area.

[0013] In some embodiments of the present invention, selecting rock samples from a target work area further comprises:

[0014] The rock samples are selected from the overlying strata or the underlying strata at depth intervals of 100 meters to 300 meters.

[0015] In some embodiments of the present invention, the dating measurement results include: U-Th / He age of apatite and U-Th / He age of zircon.

[0016] In some embodiments of the present invention, a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin further includes:

[0017] The current formation temperature of the target work area is generated according to the geothermal gradient and the surface temperature of the target work area.

[0018] In some embodiments of the present invention, determining the erosion thickness according to the geothermal gradient of the target work area and the dating measurement result includes:

[0019] Determine the depth corresponding to when the apatite enters the sealing temperature and the U-Th / He age when the apatite enters the sealing system according to the current formation temperature and the U-Th / He age of the apatite;

[0020] Determine the depth corresponding to when the zircon enters the closed system and the U-Th / He age when the zircon enters the closed system based on the current formation temperature and the U-Th / He age of the zircon;

[0021] The denudation thickness is determined according to the depth corresponding to when the apatite enters the closed temperature, the U-Th / He age when entering the closed system, the depth corresponding to when the zircon enters the closed temperature, and the U-Th / He age when entering the closed system.

[0022] In some embodiments of the present invention, determining the erosion thickness and restoring the thermal history evolution process based on the geothermal gradient of the target work area and the dating measurement results includes:

[0023] determining a first relationship between the U-Th / He age of the apatite and the burial depth or the measured formation temperature;

[0024] determining a second relationship between the U-Th / He age of the zircon and the depositional age of the rock sample;

[0025] The thermal history evolution process is restored according to the first relationship and the second relationship.

[0026] In a second aspect, the present invention provides a device for recovering the erosion thickness and thermal history evolution process of a superimposed basin, the device comprising:

[0027] A rock sample selection module is used to select rock samples from a target work area; wherein the rock sample is at least one of a sandstone sample containing apatite and zircon mineral particles, a volcanic rock sample, and a metamorphic rock sample;

[0028] a measurement result generating module, configured to perform in-situ micro-area isotope dating measurement on the rock sample to generate a dating measurement result;

[0029] The erosion thickness determination module is used to determine the erosion thickness according to the geothermal gradient of the target work area and the dating measurement result, and to restore the thermal history evolution process.

[0030] In some embodiments of the present invention, the rock sample selection module includes:

[0031] The first rock sample selection unit is used to select the rock sample from the overlying stratum or the underlying stratum of the eroded stratum in the target work area.

[0032] In some embodiments of the present invention, the rock sample selection module further includes:

[0033] The second rock sample selection unit is used to select the rock samples in the overlying stratum or the underlying stratum at a depth interval of 100 meters to 300 meters.

[0034] In some embodiments of the present invention, the dating measurement results include: U-Th / He age of apatite and U-Th / He age of zircon.

[0035] In some embodiments of the present invention, a device for recovering the erosion thickness and thermal history evolution process of a superimposed basin further includes:

[0036] The current formation temperature generating module is used to generate the current formation temperature of the target work area according to the geothermal gradient and the surface temperature of the target work area.

[0037] In some embodiments of the present invention, the erosion thickness determination module includes:

[0038] an apatite depth age determination unit, configured to determine the depth corresponding to when the apatite enters the sealing temperature and the U-Th / He age when the apatite enters the sealing system according to the current formation temperature and the U-Th / He age of the apatite;

[0039] a zircon depth age determination unit, configured to determine the depth corresponding to when the zircon entered the sealing temperature and the U-Th / He age when the zircon entered the sealing system based on the current formation temperature and the U-Th / He age of the zircon;

[0040] The erosion thickness determination unit is used to determine the erosion thickness according to the depth corresponding to when the apatite enters the sealing temperature, the U-Th / He age when entering the closed system, the depth corresponding to when the zircon enters the sealing temperature, and the U-Th / He age when entering the closed system.

[0041] In some embodiments of the present invention, the erosion thickness determination module further includes:

[0042] a first relationship determining unit, configured to determine a first relationship between the U-Th / He age of the apatite and the burial depth or the measured formation temperature;

[0043] a second relationship determining unit, configured to determine a second relationship between the U-Th / He age of the zircon and the sedimentary age of the rock sample;

[0044] A thermal history evolution process recovery unit is used to restore the thermal history evolution process according to the first relationship and the second relationship.

[0045] In a third aspect, the present invention provides a computer program product comprising a computer program / instruction, which, when executed by a processor, implements the steps of a method for recovering the erosion thickness and thermal history evolution process of a superimposed basin.

[0046] In a fourth aspect, the present invention provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps of a method for recovering the erosion thickness and thermal history evolution process of a superimposed basin are implemented.

[0047] In a fifth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin.

[0048] From the above description, it can be seen that the embodiment of the present invention provides a method and device for restoring the erosion thickness and thermal history evolution process of a superimposed basin. The corresponding method for restoring the erosion thickness and thermal history evolution process of a superimposed basin includes: first, selecting a rock sample from the target work area; wherein the rock sample is at least one of a sandstone sample, a volcanic rock sample, and a metamorphic rock sample containing apatite and zircon mineral particles; performing in-situ micro-area isotope dating on the rock sample to generate a dating measurement result; determining the erosion thickness based on the geothermal gradient of the target work area and the dating measurement result, and restoring the thermal history evolution process. Specifically, the present invention has the following beneficial effects:

[0049] First, the present invention can achieve quantitative recovery of stratum erosion thickness caused by the superposition of multiple tectonic cycles and inversion of the dynamic evolution process of thermal history under the superposition of multiple tectonic activities, effectively making up for the shortcomings of existing technologies;

[0050] Secondly, the present invention has the characteristics of quantitative recovery and a wide range of applications. It is less affected by objective geological data and is especially applicable to areas with low levels of geological research. The method only requires data such as well logging and core data, and the parameters involved are also few and easy to obtain.

[0051] Finally, the present invention is based on actual geological data and low-temperature thermochronological data of mineral particles. It is suitable for restoring the thickness of stratum erosion caused by the longitudinal superposition of multiple periods of structure in complex superimposed basins, and is also suitable for restoring the thickness of stratum erosion caused by tectonic uplift in single-type basins. This method has been well promoted and applied in the Sichuan Basin, Tarim Basin, Ordos Basin, Bohai Bay Basin, etc. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0053] Figure 1 Schematic diagram of a process for recovering the erosion thickness and thermal history evolution process of a superimposed basin in an embodiment of the present invention;

[0054] Figure 2 1 is a flow chart of step 100 of a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin in an embodiment of the present invention;

[0055] Figure 3 Another flow chart of step 100 of a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin in an embodiment of the present invention;

[0056] Figure 4 Another schematic flow chart of a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin in an embodiment of the present invention;

[0057] Figure 5 3 is a flow chart of step 300 of a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin in an embodiment of the present invention;

[0058] Figure 6 3 is another flow chart of step 300 of a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin in an embodiment of the present invention;

[0059] Figure 7 Schematic diagram of a process for restoring the erosion thickness and thermal history evolution process of a superimposed basin in a specific embodiment of the present invention;

[0060] Figure 8 This is a graph showing the relationship between the apatite and zircon He ages, burial depth, and current ground temperature in Well X, Yuanba area, Sichuan Basin, in a specific embodiment of the present invention;

[0061] Figure 9 A block diagram of a device for recovering the erosion thickness and thermal history evolution process of a superimposed basin in an embodiment of the present invention;

[0062] Figure 10 Schematic diagram of the structure of an electronic device in an embodiment of the present invention. DETAILED DESCRIPTION

[0063] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.

[0064] It will be understood by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Thus, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0065] It should be noted that the terms "including" and "having" and any variations thereof in the specification and claims of this application and the above-mentioned drawings are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or are inherent to these processes, methods, products or devices. The embodiments in this application and the features described in the embodiments may be combined with each other unless there is a conflict. The present application will be described in detail below with reference to the drawings and in conjunction with the embodiments.

[0066] The acquisition, storage, use, and processing of data in the technical solution of this application comply with relevant laws and regulations.

[0067] In the existing technology, there are relatively many methods for restoring the thickness of strata erosion, including stratum comparison fitting method and sedimentation-tectonic comprehensive analysis method, sandstone porosity method, sedimentation rate method, well logging curve method, paleo-geothermal method, etc. However, there is almost no mention of the restoration of the dynamic evolution of thermal history under the superposition of multiple periods of tectonic activities.

[0068] (1) Stratigraphic correlation method: For example, the patent "Method for measuring stratum erosion thickness in multi-stage differential uplift areas of superimposed basins" (patent application number: CN 201510375316.3) uses a curve fitting method to obtain the morphology, development characteristics and distribution trend of the stratum before erosion based on the development characteristics, extension trend, fault cutting relationship, fold deformation characteristics and other stratigraphic elements of the uneroded stratum interface near the unconformity interface based on drilling and fine interpretation of seismic profiles, and infers the amount of erosion in the eroded area by the stratum fitting method. The shortcomings of this method are: when the erosion area is large, the stratum thickness varies greatly in the lateral direction, and the error will be large; at the same time, the erosion thickness calculated by this method is greatly affected by the quality of seismic data, and the influence of human subjective factors cannot be ignored.

[0069] (2) Sandstone porosity method: For example, the patent "A method and device for obtaining the erosion thickness of a sedimentary basin" (patent application number: CN 201410374600.4) is based on the fitting relationship between the porosity, depth and lithology of the upper and lower adjacent strata. The fitting relationship between the porosity, depth and lithology of the upper and lower adjacent strata is compared to determine whether the lower stratum of the two upper and lower adjacent strata has been eroded. Then, the erosion thickness during sedimentation is obtained according to the compensation between the upper and lower strata interfaces. The shortcomings of this method are: in many cases, there is no necessary relationship between the porosity of the strata and the depth and lithology, and it is difficult to establish a fitting curve between them. In addition, the fitting relationship between the upper and lower adjacent strata is often inconsistent, and the erosion thickness cannot be calculated through the fitting relationship. This method assumes that the lithology of the eroded stratum is the same as that of the residual stratum, which is not true in many cases. The values ​​of many parameters in this method have great uncertainty and cannot be quantitatively restored. This method is only applicable to sandstone strata.

[0070] (3) Sedimentation rate method: For example, the patent "A method for restoring the thickness of unconformity erosion using cyclic analysis" (patent application number: CN 201710015647.5) combines regional geological data, calculates the original thickness of the stratum before compaction through decompaction correction, calculates the average stratum sedimentation rate of each sedimentary period, restores the sedimentation rate curve and performs frequency division processing, and infers the period, peak and trough corresponding to the missing curve segment; restores the complete sedimentation rate curve based on spectrum analysis, and calculates the unconformity erosion thickness by integrating the complete sedimentation rate curve. The shortcomings of this method are: the original thickness of the stratum before compaction is difficult to obtain, and the stratum sedimentation rate calculated based on this is not very reliable; the missing sedimentation curve segment often does not have periodicity, and it is impossible to perform frequency division processing to obtain the corresponding trough and peak; the parameters required by this method, such as erosion rate and absolute age of the stratum, are sometimes difficult to obtain accurate values, and it is not applicable to areas where the degree of stratigraphic research is not high.

[0071] (4) Well logging method (mudstone acoustic time difference method, porosity or rock density method, etc.): Mudstone acoustic time difference method, such as the patent "Method and device for determining erosion thickness" (patent application number: CN 202111333588.9), uses a first acoustic radar to use acoustic time difference to determine the unconformity interface at any position in the formation. Multiple sets of second acoustic radars are successively deployed outward with the position of the first acoustic radar as the reference point. The multiple sets of second acoustic radars use acoustic time difference to determine multiple depths of the unconformity interface based on the reference point, and determine the critical line of the unconformity interface. The three-dimensional simulation map of the thickness of the eroded formation is determined to obtain the formation erosion thickness. The shortcomings of this method are: this method is only applicable to areas with large erosion; this method assumes that the compaction process of mudstone is irreversible and not affected by time, and the thickness of sediments above the unconformity is less than the actual erosion thickness, but in fact some of the above conditions may not be met; this method often requires a lot of drilling and logging data and cannot be used in medium and low exploration areas. The porosity method, such as the patent "A Method for Determining the Erosion Thickness of Ancient Clastic Rock Strata" (patent application number: CN 202110893186.8), determines the unconformity interface of the stratum based on sudden changes in logging data, changes in seismic interfaces, and changes in core contact surfaces; selects rock porosity data or rock density data of the eroded stratum from the unconformity interface; determines the original depth of the eroded stratum based on the rock porosity data or rock density data, and then determines the thickness of the eroded stratum based on the drilling depth. The shortcomings of this method are that it is only applicable to areas with large amounts of erosion; in many cases, the porosity or rock density data do not necessarily change continuously with depth, making it impossible to determine the original depth of the eroded stratum based on the above data, and naturally, it is impossible to obtain the thickness of the eroded stratum.

[0072] (5) Paleo-geothermal method (paleo-geothermal gradient method, inclusion temperature measurement method, vitrinite reflectance method, apatite fission track method, etc.): For example, the patent "Method, system and equipment for calculating the angular unconformity erosion thickness of sedimentary basins" (patent application number: CN202210270337.9) includes obtaining the uniform temperature and pressure data of fluid inclusion temperature data based on fluid inclusion homogenization temperature test and fluid inclusion pressure calculation, and respectively using layer projection correction and linear fitting processing to obtain temperature gradient correction data and pressure gradient correction data; and calculating the erosion thickness of each sample in the sedimentary basin by the erosion thickness calculation formula. The shortcomings of this method are: this method is based on the fluid inclusion analysis data, but the key problem is that fluid inclusions are not always developed in sedimentary strata, so there are great defects in practical application; at the same time, due to the large uncertainty in the burial history-thermal history curve of the stratum, the maximum paleo-burial depth of the underlying stratum obtained based on the vitrinite reflectance and regional stratum comparison also has large uncertainty. For example, the patent "A Method for Determining Stratum Erosion Thickness" (patent application number: CN201810374163.4) combines regional geological data and uses vitrinite reflectance to obtain displacement pressure. Standard curves of displacement pressure and depth are established for uneroded and eroded strata, respectively. The two curves are fitted in the same coordinate system. In the fitted coordinate system, the distance difference between the intersection of the standard curve of the eroded stratum with the ordinate and the depth of the surface of the eroded stratum is the stratum erosion thickness. The shortcomings of this method are: Lower Paleozoic strata often do not contain vitrinite, so it is impossible to obtain accurate vitrinite reflectance, and naturally it is impossible to obtain accurate displacement pressure; in addition, when the thickness of the stratum deposited above the erosion surface exceeds the thickness of the eroded layer, this method cannot be used. For example, the patent "Method and Apparatus for Determining Erosion Thickness" (Patent Application No.: 201911010562.3) combines a topographic profile of the study area to obtain the temperature value corresponding to each sample point under the corresponding geological history, determine the corresponding paleotopography of the study area under the corresponding geological history, and thus determine the corresponding erosion thickness of the study area under the corresponding geological history. However, this method has shortcomings: it is only applicable to the reconstruction of erosion thickness of present-day topography; the thermal history simulation curve established by this method is based on a large number of data samples, making it unsuitable for areas with less extensive stratigraphic research. For example, the patent "Method, Apparatus, Equipment, and Storage Medium for Determining Erosion Thickness" (Patent Application No.: CN 202010465532.8) performs thermal history simulation on a rock sample to obtain a first thermal history curve. The inflection point and low-temperature segment of the first thermal history curve are determined. Based on the low-temperature segment, a trend fitting is performed on the high-temperature segment of the first thermal history curve. Based on the trend-fitted high-temperature and low-temperature segments, a second thermal history curve of the rock sample is obtained, thereby determining the erosion thickness of the rock stratum at various historical times.The shortcomings of this method are: the thermal history simulation curve established by this method is based on a large number of data samples, and is not applicable to areas where the degree of stratigraphic research is not high; the parameters such as paleotemperature and paleotemperature gradient required by this method are sometimes difficult to obtain accurate values, so this method is not applicable to areas where the degree of stratigraphic research is not high.

[0073] The paper "A New Method for Recovering Stratum Denudation Thickness Using Vitrinite Reflectance" proposes a method for estimating denudation thickness using the maximum paleotemperature method (vitrinite reflectance Ro). Specifically, the denudation thickness at the top of a structural layer is estimated using the difference in vitrinite reflectance (Ro) between the upper and lower structural layers. The Ln(Ro)-H linear regression relationship for that layer is extrapolated to Ln(0.2) to obtain an approximate paleosurface position. The difference between the obtained paleosurface position and the unconformity position is the denudation thickness. However, this method has several drawbacks: The Ro value is used to determine the maximum paleotemperature in a stratum using a variety of dynamic models or empirical formulas, each with varying scopes of applicability, making the method less universal. Furthermore, if volcanic activity or magma intrusion occurs in the strata under study, this can cause changes in the geothermal field within a certain range, leading to anomalous changes in the vitrinite reflectance of the strata. Consequently, the denudation thickness derived from this method is unreliable. The paper "Recovering Stratigraphic Erosion Thickness Based on Milankovitch Astronomical Cycles: A Case Study of the Qingshankou Formation in the X Oilfield of the Songliao Basin" uses spectral analysis, wavelet transform, and filter analysis to identify periodic variations in astronomical orbital parameters (eccentricity, slope, and precession). This method then quantitatively evaluates the thickness of sedimentary cycles driven by cyclical fluctuations in Earth's surface climate, thereby calculating stratigraphic erosion thickness. However, this method has several limitations: It requires excluding the influence of faults when selecting data, lacks predictive interpretation in less-explored areas, and the erosion thickness derived from this method is only the "minimum apparent erosion thickness" at present-day burial depth.

[0074] In summary, the existing methods have the following shortcomings:

[0075] First, most existing methods achieve quantitative recovery of erosion thickness by establishing a fitting relationship. However, the establishment of such a fitting mathematical relationship is often non-unique, resulting in multiple solutions for the erosion thickness value and large errors.

[0076] Second, many of the parameter values ​​involved in existing methods have large uncertainties, making it difficult to obtain accurate denudation thickness. In addition, existing methods often require a large amount of drilling, logging, and core analysis data, and cannot be used in areas with medium and low exploration levels.

[0077] Third, existing methods all have their own applicable conditions and limitations. For example, the vitrinite reflectance method requires a large amount of test data and is significantly affected by thermal events. The stratigraphic correlation method based on seismic interpretation must be based on a relatively small erosion area in the study area and is not applicable to areas with widespread erosion.

[0078] Fourth, existing methods cannot provide the thermal history evolution process of superimposed basins under the superposition of multiple periods of tectonic activities.

[0079] Example 1:

[0080] Based on the above reasons, the embodiment of the present invention provides a specific implementation method of a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin, see Figure 1 , specifically including the following contents:

[0081] Step 100: Selecting a rock sample from a target work area; wherein the rock sample is at least one of a sandstone sample containing apatite and zircon mineral particles, a volcanic rock sample, and a metamorphic rock sample;

[0082] Step 200: performing in-situ micro-area isotope dating measurement on the rock sample to generate a dating measurement result;

[0083] Step 300: Determine the erosion thickness based on the geothermal gradient of the target work area and the dating measurement results, and restore the thermal history evolution process.

[0084] The present invention provides a method for recovering the erosion thickness and thermal history of superimposed basins, which can achieve the quantitative recovery of the stratum erosion thickness caused by the superposition of multiple tectonic cycles and the inversion of the dynamic evolution process of thermal history under the superposition of multiple tectonic activities, effectively making up for the shortcomings of the existing technology and having important significance for the research and exploration evaluation of oil and gas accumulation in complex superimposed basins.

[0085] Example 2:

[0086] Regarding step 100, apatite and zircon provide rich information in different types of rock samples, including chronology, sedimentary environment, source characteristics, metamorphism, etc.

[0087] Apatite in sandstones originates from the source region of the sediments. It can come from eroded phosphate rocks, bioclastics (such as bones and teeth), or biochemical precipitation during marine sedimentation. The morphology and composition of apatite may reflect the source and depositional environment of the sediments; for example, marine environments may promote the formation of phosphorus-rich minerals.

[0088] Zircons in sandstones are eroded from source rocks (such as granite and gneiss). Zircons' resistance to erosion allows them to survive sedimentary transport. U-Pb dating of zircons can be used to determine the maximum depositional age of sediments—the age of the youngest zircon grains in a sedimentary layer—helping constrain the formation and evolutionary history of sedimentary basins.

[0089] Apatite is less common in volcanic rocks than in sedimentary rocks, but it can also form, particularly in phosphorus-rich magmas. Its presence may indicate a unique magma source or be a product of later hydrothermal processes. In some cases, it can also serve as a tracer mineral in volcanic sediments.

[0090] Zircon in volcanic rocks is often associated with the evolution of magma. It crystallizes early in the cooling process and is therefore often used to date volcanic rocks. Isotopic analysis of zircon can reveal the source characteristics and evolution of the magma.

[0091] Apatite in metamorphic rocks can be a remnant of original sedimentary minerals or newly formed during metamorphism. Apatite can recrystallize under high temperature and pressure, and changes in its composition and morphology can indicate the extent of metamorphism.

[0092] Zircon is relatively stable during metamorphism and can preserve information about the original magma or sedimentary age, even after significant metamorphism. The trace element and isotopic composition of zircon may change during metamorphism, thus recording changes in metamorphic temperature, pressure, and liquid phase.

[0093] Regarding step 200, in situ micro-area isotope dating is a method for measuring isotope ratios at the micron or submicron scale, which can provide age information for a single mineral grain or specific part of a rock sample. The key advantage of this technology is that it can determine the time of formation or change without destroying the sample. Commonly used in situ micro-area isotope dating techniques include laser ablation multi-collector inductively coupled plasma mass spectrometry (LA-ICP-MS), secondary ion mass spectrometry (SIMS), electron probe microanalysis (EPMA), etc.

[0094] Preferably, in situ micro-isotope dating techniques can be combined with electron microscopy to analyze the mineral microstructure and isotope analysis using SIMS or LA-ICP-MS. This can provide valuable information about the formation and evolution history of the rock, as well as subsequent geological processes.

[0095] In situ micro-isotope dating techniques are particularly useful for analyzing heat-resistant and / or corrosion-resistant minerals such as zircon, apatite, and single zircon grains, as these minerals can survive diverse geological processes in the Earth's crust. These techniques provide a more precise understanding of the formation times of these minerals, leading to a deeper understanding of Earth's history and evolution.

[0096] With respect to step 300, denudation thickness refers to the thickness of rock removed from the surface due to external forces such as weathering and erosion. Denudation causes rock below the surface to gradually rise to the surface, a process known as denudation.

[0097] Thermal history evolution refers to the thermal changes that rocks in the Earth's crust experience over time, and is usually related to regional tectonic activity and crustal movement.

[0098] In some embodiments of the present invention, see Figure 2 , step 100 includes:

[0099] Step 101: Select the rock sample from the overlying stratum or the underlying stratum of the eroded stratum in the target work area.

[0100] In some embodiments of the present invention, see Figure 3 , step 100 further includes:

[0101] Step 102: Select the rock samples from the overlying strata or the underlying strata at depth intervals of 100 meters to 300 meters.

[0102] Steps 101 and 102 include: obtaining data on the eroded strata and their overlying or underlying strata based on the lithologic characteristics, logging data, stratigraphic stratification, depth data, etc. of the drilling profile in the study area, and obtaining representative core or cuttings samples at depth intervals of 100m-300m.

[0103] Furthermore, rock samples can be coarsely crushed to 5mm-10mm, and then finely crushed according to the coarseness of the sample crystals. Coarsely crystallized samples can be crushed to approximately 40 mesh, and finely crystallized rocks can be crushed to approximately 60 mesh. The samples are then sieved through a stainless steel analytical sieve (>30 mesh). The fine fraction is then washed, electrostatically separated, and subjected to heavy liquid separation to obtain apatite, zircon, and other heavy mineral concentrates. Under a 200x stereomicroscope, the dating mineral particles are observed, and 3-5 mineral particles with relatively good crystal shape, no breakage, minimal mineral inclusions, cracks, or impurities are selected for (U-Th) / He dating analysis using an in-situ micro-area isotope dating device.

[0104] In some embodiments of the present invention, the dating measurement results include: U-Th / He age of apatite and U-Th / He age of zircon.

[0105] Specifically, based on the measurement of the concentration and distribution of Uranium (U), Thorium (Th), and Helium (He) in minerals, the accumulation of helium produced during the decay of naturally occurring radioactive elements uranium and thorium over time is used to determine the cooling age of rocks or minerals, thereby inferring their thermal history.

[0106] Apatite U-Th / He dating (AHe) is a low-temperature thermochronology technique applicable to temperatures ranging from approximately 30°C to 90°C, corresponding to depths of several thousand meters in the Earth's crust. AHe dating can be used to study the low-temperature thermal history of the near-surface, such as erosion, uplift, and other surface processes.

[0107] Uranium and thorium naturally decay within the apatite lattice, producing helium atoms. Helium is an inert gas that accumulates in minerals over time. If the mineral is heated above a certain temperature (the closure temperature), the helium escapes. By measuring the helium content in apatite, along with the uranium and thorium concentrations, the time it took for the apatite crystal to cool through the closure temperature can be calculated—the U-Th / He age. Apatite U-Th / He age determination has a wide range of applications in areas such as erosion rate measurement, tectonic uplift history, sedimentary basin evolution, and volcanic history.

[0108] Zircon U-Th / He dating (ZHe) is another low-temperature thermochronology method. Its testing temperature range is higher than that of apatite U-Th / He, approximately between 160-220°C. ZHe dating can be used to study deeper crustal cooling events.

[0109] Similar to apatite, the decay of uranium and thorium in zircon also produces helium. However, due to the tighter lattice structure of zircon, helium only begins to significantly escape at higher temperatures. Measuring the helium content in zircon, combined with the uranium and thorium concentrations, allows calculation of the zircon's U-Th / He age. Zircon U-Th / He dating is used to study deeper geological history, such as continental collision, the uplift and erosion of orogenic belts, and the cooling history of igneous rocks.

[0110] In some embodiments of the present invention, see Figure 4 A method for restoring the erosion thickness and thermal history evolution process of a superimposed basin, further comprising:

[0111] Step 400: Generate the current formation temperature of the target work area according to the geothermal gradient and the surface temperature of the target work area.

[0112] Specifically, the geothermal gradient of the study area is calculated through the measured temperature data from well logging, and then combined with the surface temperature calculation to obtain the current formation temperature of rock samples at different depths.

[0113] In some embodiments of the present invention, see Figure 5 , step 300 includes:

[0114] Step 301: Determine the depth corresponding to when the apatite enters the sealing temperature and the U-Th / He age when the apatite enters the sealing system based on the current formation temperature and the U-Th / He age of the apatite;

[0115] Step 302: Determine the depth corresponding to when the zircon enters the sealing temperature and the U-Th / He age when the zircon enters the sealing system based on the current formation temperature and the U-Th / He age of the zircon;

[0116] For steps 301 and 302, on this basis, with the depth of the rock sample and its current formation temperature as the vertical coordinate and the (U-Th) / He age of the rock sample as the horizontal coordinate, a curve is established to determine the relationship between the apatite and zircon (U-Th) / He ages and the depositional age of the sample in the rock sample. When the apatite (U-Th) / He age in the same rock sample is less than the depositional age of the sample and the zircon (U-Th) / He age is less than or equal to the depositional age of the sample, the temperature and depth corresponding to when the apatite and zircon (U-Th) / He ages in the area are close to zero can be determined based on the evolution relationship diagram of the apatite and zircon He ages with the current formation temperature. This temperature and depth are the closure temperature of apatite and zircon entering the (U-Th) / He dating system and the burial depth of the rock sample at this time, respectively.

[0117] Step 303: Determine the denudation thickness based on the depth corresponding to when the apatite enters the sealing temperature, the U-Th / He age when entering the closed system, the depth corresponding to when the zircon enters the sealing temperature, and the U-Th / He age when entering the closed system.

[0118] The calculation formula for stratum erosion thickness is as follows:

[0119] H=Vd×Td, Vd=(Da-Dz) / (Aa-Az);

[0120] Where H is the denudation thickness (m), Vd is the denudation rate (m / Ma), Td is the duration of uplift and denudation (Ma), Td is the difference between the uplift and denudation times determined by zircon (U-Th) / He dating and apatite (U-Th) / He dating for the same stratum, Da is the depth at which apatite entered the closure temperature (m), Dz is the depth at which zircon entered the closure temperature (m), Aa is the He age at which apatite entered the closure system (Ma), and Az is the He age at which zircon entered the closure system (Ma). The ratio of the difference in burial depths at which apatite and zircon entered the (U-Th) / He dating closure system to the difference in He ages between apatite and zircon in the same rock sample can be used to calculate the erosion rate of the stratum during the uplift period, and thus the uplift and denudation thickness of the stratum.

[0121] In some embodiments of the present invention, see Figure 6 , step 300 further includes:

[0122] Step 304: determining a first relationship between the U-Th / He age of the apatite and the burial depth or the measured formation temperature;

[0123] Step 305: Determine a second relationship between the U-Th / He age of the zircon and the sedimentary age of the rock sample;

[0124] Step 306: Restore the thermal history evolution process according to the first relationship and the second relationship.

[0125] In steps 304 to 306, in the relationship curve between apatite (U-Th) / He age and depth and current formation temperature (established in step 301), the apatite (U-Th) / He age in the rock sample gradually decreases with increasing burial depth or measured formation temperature, indicating that the rock sample has experienced the sealing temperature of apatite. The zircon (U-Th) / He age in the rock sample is less than or equal to the sample deposition age, indicating that the highest paleotemperature of the rock sample in the geological age corresponding to the zircon (U-Th) / He age was close to the zircon closure temperature of 170°C, and then the stratum began to uplift and erode, and the paleotemperature gradually decreased; within the geological age corresponding to the average value of the (U-Th) / He age after the zircon entered the closed system and the geological age corresponding to the average value of the He age after the apatite entered the closed system, the paleotemperature of the stratum was within the closure temperature range of apatite and zircon; from the geological age corresponding to the average value of the (U-Th) / He age after the apatite entered the closed system to the present, it is judged that the paleotemperature of the stratum was within the closure temperature range of apatite and zircon.

[0126] As can be seen from the above description, an embodiment of the present invention provides a method for recovering the erosion thickness and thermal history evolution process of a superimposed basin. This method mainly relies on low-temperature thermochronological dating of mineral particles to recover the stratum erosion thickness and thermal history dynamic evolution. The main processes involved are: (1) rock sample acquisition and selection of suitable mineral particles, and mineral particle (U-Th) / He dating analysis; (2) recovery of stratum erosion thickness; and (3) recovery of thermal history evolution process. The present invention can achieve quantitative recovery of stratum erosion thickness caused by the superposition of multiple tectonic cycles and inversion of thermal history dynamic evolution process under the superposition of multiple tectonic activities, effectively making up for the shortcomings of existing technologies and having important significance for oil and gas reservoir research and exploration evaluation in complex superimposed basins.

[0127] Example 3:

[0128] In a specific embodiment, the present invention also takes the Yuanba area of ​​the Sichuan Basin as an example to provide a specific embodiment of a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin, see Figure 7 , specifically including the following steps.

[0129] The Sichuan Basin is a typical superimposed basin. In recent years, with the deepening of oil and gas exploration, a series of large and medium-sized natural gas fields have been discovered in the Sichuan Basin, including those in Puguang, Longgang, Yuanba, Chuanxi, and Anyue, demonstrating the enormous exploration potential of deep and ultra-deep strata in superimposed basins. The Sichuan Basin has experienced multiple phases of tectonic activity, including the Caledonian, Hercynian, Indosinian, Yanshan, and Himalayan tectonic events, and has undergone multiple uplift and erosion events. The Yuanba area is a key area for natural gas exploration in northeastern Sichuan. Research on the evolution of stratigraphic erosion thickness and thermal history caused by the superposition of multiple tectonic cycles since the Mesozoic in this area is of great significance for the study of the "four histories" of the region (sedimentary history, tectonic history, hydrocarbon generation history, and oil and gas migration and accumulation history), reservoir formation mechanisms, and even oil and gas exploration and evaluation.

[0130] Step S1: Acquire regional geological, logging, drilling data and rock samples in the study area, prepare mineral particles and select suitable dating mineral particles, and conduct (U-Th) / He dating analysis using an in situ micro-area isotope dating device.

[0131] Specifically, regional geological data and rock samples are obtained, suitable mineral particles are selected, and their (U-Th) / He dating analysis is performed.

[0132] The study area is located in northeastern Sichuan. Well X is a natural gas exploration well in the area. The Mesozoic cored intervals in this well are the Xujiahe Formation Xu 2 and Xu 4 Members of the Upper Triassic, the Ziliujing Formation of the Lower Jurassic, and the Jianmenguan Formation of the Lower Cretaceous. These strata are mainly composed of siltstone, fine sandstone, argillaceous sandstone, and sand-bearing mudstone.

[0133] First, based on the lithologic characteristics, logging data, stratigraphic stratification, and depth data of the X-well drilling profile in the area, data on the eroded strata and their overlying or underlying strata were obtained.

[0134] Then, representative core or cutting samples are collected at depth intervals of approximately 100m-300m, with a sample mass of no less than 500g. The collected rock samples are coarsely crushed to 5mm-10mm, and then finely crushed according to the coarseness and fineness of the sample crystals. Coarsely crystalline samples are crushed to approximately 40 mesh, and finely crystalline rocks are crushed to approximately 60 mesh. The samples are sieved with a stainless steel analytical sieve (>30 mesh), and the fine fraction is then washed, electrostatically separated, and subjected to heavy liquid separation to obtain apatite, zircon and other heavy mineral concentrates.

[0135] Finally, under a 200x stereo microscope, observe the dating mineral grains, select 3-5 mineral grains with relatively good crystal shape, no breakage, as few mineral inclusions as possible, no cracks and impurities, and use an in-situ micro-area isotope dating device to perform (U-Th) / He dating analysis.

[0136] Step S2: Establish a relationship curve between the (U-Th) / He age of mineral particles in the rock sample and the depth and current formation temperature to determine that the area is in the main period of cooling, uplift and erosion, and use the formula to calculate the formation erosion rate and thickness during this period.

[0137] Specifically, the thickness of stratum erosion caused by the superposition of multiple tectonic cycles since the Mesozoic in Well X in the Yuanba area was quantitatively restored. The current geothermal gradient in the area was found to be 1.7°C / 100m and the average surface temperature was 24°C based on the measured temperature data of the well. The current formation temperature of rock samples at different depths was obtained by calculation. The relationship curves of apatite and zircon He ages with depth and current formation temperature were established, with the depth of the rock samples and their current formation temperature as the ordinates and the (U-Th) / He age of the rock samples as the abscissas. Figure 8 ), averaging the He age data of 3-5 apatites and 3-5 zircons from the same rock sample at different depths in Well X can respectively obtain the average He ages of apatites and zircons at different depths, based on which it can be determined that the area was in the main period of cooling, uplift and erosion.

[0138] The average value of apatite He age is 0.2Ma to 36.4Ma, indicating that the area was under intense tectonic uplift during the Paleogene-Neogene. In addition, according to the evolution relationship between apatite He age and current stratum temperature, it can be judged that the closure temperature of the apatite (U-Th) / He system in this area is about 96℃ (the black line on the far left is the apatite (U-Th) / He age curve) and the burial depth is about 4300m (the current ground temperature and depth corresponding to the apatite (U-Th) / He age curve close to 0). And the zircon in the rock sample of the second section ( Figure 8The (U-Th) / He age in the lowest triangle is younger than the stratigraphic depositional age, indicating that the sample likely reached a maximum paleotemperature close to the zircon closure temperature of around 170°C, corresponding to a burial depth of approximately 4800 m (the depth corresponding to the triangle). Using the formula for calculating stratigraphic denudation thickness (H = Vd × Td, where Vd = (Da-Dz) / (Aa-Az)), we calculated an average denudation rate of 110 m / Ma and an average cooling rate of 1.9°C / Ma, yielding a calculated uplift and denudation thickness of approximately 4000 m during this period.

[0139] At the same time, based on the relationship between the apatite (U-Th) / He age and the burial depth of the rock sample, the age-elevation method can also be used to calculate the erosion rate and erosion thickness of the stratum during this uplift period. This method can be used to verify the results of step S2. In fact, the stratum erosion thickness obtained by these two methods is consistent.

[0140] Step S3: Determine the (U-Th) / He age and closure temperature of different mineral particles in the rock samples in the study area, and invert the thermal history evolution process of the strata during different tectonic activity geological periods.

[0141] The thermal evolution of the Mesozoic formation in Well X of the Yuanba area was restored. In the relationship curve between apatite (U-Th) / He age and depth as well as the current formation temperature ( Figure 8 ), the apatite He age in the rock samples gradually decreases with the increase of burial depth or measured formation temperature, indicating that the rock samples have experienced the closure temperature of apatite, reflecting that the closure temperature of the apatite (U-Th) / He system in this area is about 96℃; Figure 8 It can be seen that the zircon He age in the Xu2 rock samples is less than or equal to the sample deposition age, indicating that the highest paleotemperature of the rock samples was close to the zircon closure temperature of 170℃ around 183Ma (average He age of zircons in the Xu2 member), and then the strata began to uplift and erode, and the paleotemperature gradually decreased; in the geological age range of 183Ma-36.4Ma (average He age of apatite in the Jianmenguan Formation), the paleotemperature of the strata was between 96℃-170℃; in the geological age range of 36.4Ma-0.2Ma, the paleotemperature of the strata was between 96℃-24℃; it can be seen that the thermal history evolution process of the Yuanba area in northeastern Sichuan since the Mesozoic has been a gradual cooling, and the cooling rate has large differences.

[0142] As can be seen from the above description, a specific embodiment of the present invention provides a method for restoring the erosion thickness and thermal history evolution process of a superimposed basin mainly by means of low-temperature thermochronological dating of mineral particles, including the following steps: acquisition of regional geological data and rock samples, preparation of mineral particles and selection of suitable dating particles, and (U-Th) / He dating analysis of mineral particles; judging the main uplift and erosion period based on the average value of the He age data of apatite and zircon in the same sample, establishing a relationship curve between the He age of apatite and zircon and the depth and the current formation temperature, and calculating the formation uplift and erosion rate and its erosion thickness by using the ratio of the burial depth difference of apatite and zircon in the same rock sample entering the (U-Th) / He closed system and the He age difference when entering the closed system; then using the average value of the He age of apatite and zircon after entering the closed system and its closed temperature to comprehensively judge the range of ancient geotemperature in different geological periods, and then the evolution process of thermal history in different geological periods can be judged. This method can quantitatively restore the thickness of stratum erosion caused by the superposition of multiple tectonic cycles and invert the dynamic evolution process of thermal history under the superposition of multiple tectonic activities, effectively making up for the shortcomings of existing technologies. It has the characteristics of quantitative recovery, little influence from objective geological data, few parameters involved and easy acquisition. It is of great significance for the study of the "four histories" of superimposed basins, reservoir formation mechanisms and even oil and gas exploration and evaluation, and has strong geological applicability.

[0143] Example 4:

[0144] Based on the same inventive concept, the embodiments of the present application also provide a device for recovering the erosion thickness and thermal history evolution process of a superimposed basin, which can be used to implement the method described in the above embodiments, such as the following embodiments. Since the principle of solving the problem by the device for recovering the erosion thickness and thermal history evolution process of a superimposed basin is similar to the method for recovering the erosion thickness and thermal history evolution process of a superimposed basin, the implementation of the device for recovering the erosion thickness and thermal history evolution process of a superimposed basin can refer to the implementation of the method for recovering the erosion thickness and thermal history evolution process of a superimposed basin, and the repetitions will not be repeated. As used below, the term "unit" or "module" can be a combination of software and / or hardware that implements a predetermined function. Although the system described in the following embodiments is preferably implemented in software, implementation in hardware, or a combination of software and hardware, is also possible and conceived.

[0145] The embodiment of the present invention provides a specific embodiment of a superimposed basin erosion thickness and thermal history evolution process recovery method, see Figure 9 A device for recovering the erosion thickness and thermal history evolution process of a superimposed basin includes:

[0146] The rock sample selection module 10 is used to select rock samples from the target work area; wherein the rock sample is at least one of a sandstone sample containing apatite and zircon mineral particles, a volcanic rock sample, and a metamorphic rock sample;

[0147] A measurement result generating module 20 is used to perform in-situ micro-area isotope dating measurement on the rock sample to generate a dating measurement result;

[0148] The erosion thickness determination module 30 is used to determine the erosion thickness according to the geothermal gradient of the target work area and the dating measurement result, and to restore the thermal history evolution process.

[0149] In some embodiments of the present invention, the rock sample selection module includes:

[0150] The first rock sample selection unit is used to select the rock sample from the overlying stratum or the underlying stratum of the eroded stratum in the target work area.

[0151] In some embodiments of the present invention, the rock sample selection module further includes:

[0152] The second rock sample selection unit is used to select the rock samples in the overlying stratum or the underlying stratum at a depth interval of 100 meters to 300 meters.

[0153] In some embodiments of the present invention, the dating measurement results include: U-Th / He age of apatite and U-Th / He age of zircon.

[0154] In some embodiments of the present invention, a device for recovering the erosion thickness and thermal history evolution process of a superimposed basin further includes:

[0155] The current formation temperature generating module is used to generate the current formation temperature of the target work area according to the geothermal gradient and the surface temperature of the target work area.

[0156] In some embodiments of the present invention, the erosion thickness determination module includes:

[0157] an apatite depth age determination unit, configured to determine the depth corresponding to when the apatite enters the sealing temperature and the U-Th / He age when the apatite enters the sealing system according to the current formation temperature and the U-Th / He age of the apatite;

[0158] a zircon depth age determination unit, configured to determine the depth corresponding to when the zircon entered the sealing temperature and the U-Th / He age when the zircon entered the sealing system based on the current formation temperature and the U-Th / He age of the zircon;

[0159] The erosion thickness determination unit is used to determine the erosion thickness according to the depth corresponding to when the apatite enters the sealing temperature, the U-Th / He age when entering the closed system, the depth corresponding to when the zircon enters the sealing temperature, and the U-Th / He age when entering the closed system.

[0160] In some embodiments of the present invention, the erosion thickness determination module further includes:

[0161] a first relationship determining unit, configured to determine a first relationship between the U-Th / He age of the apatite and the burial depth or the measured formation temperature;

[0162] a second relationship determining unit, configured to determine a second relationship between the U-Th / He age of the zircon and the sedimentary age of the rock sample;

[0163] A thermal history evolution process recovery unit is used to restore the thermal history evolution process according to the first relationship and the second relationship.

[0164] As can be seen from the above description, an embodiment of the present invention provides a device for recovering the erosion thickness and thermal history evolution process of a superimposed basin, including: a sensitive logging data selection module for selecting sensitive logging data corresponding to the target reservoir to be identified; wherein the sensitive logging data is used to characterize the lithology of the target reservoir; a sensitive logging data classification module for classifying the sensitive logging data based on the correlation between the sensitive logging data and preset calibration logging data; a target reservoir lithology identification module for identifying the lithology of the target reservoir based on the calibration logging data, the correlation coefficient between multiple sensitive logging data, and the classified sensitive logging data. Specifically, the present invention has the following beneficial effects:

[0165] Embodiment 5:

[0166] The embodiments of the present application also provide a specific implementation of an electronic device capable of implementing all steps in the method for restoring the erosion thickness and thermal history evolution process of a superimposed basin in the above embodiment, see Figure 10 , electronic equipment specifically includes the following:

[0167] Processor 1201, memory 1202, communications interface 1203, and bus 1204;

[0168] The processor 1201, the memory 1202, and the communication interface 1203 communicate with each other via the bus 1204; the communication interface 1203 is used to implement information transmission between the server device and the client device and other related devices;

[0169] The processor 1201 is used to call the computer program in the memory 1202. When the processor executes the computer program, all steps of the method for recovering the erosion thickness and thermal history evolution process of the superimposed basin in the above embodiment are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0170] Selecting a rock sample from the target work area; wherein the rock sample is at least one of a sandstone sample containing apatite and zircon mineral particles, a volcanic rock sample, and a metamorphic rock sample;

[0171] performing in-situ micro-area isotope dating on the rock sample to generate a dating result;

[0172] The erosion thickness is determined based on the geothermal gradient of the target work area and the dating measurement results, and the thermal history evolution process is restored.

[0173] Example 6:

[0174] An embodiment of the present application also provides a computer-readable storage medium that can implement all the steps of the method for restoring the erosion thickness and thermal history evolution process of the superimposed basin in the above-mentioned embodiment. The computer-readable storage medium stores a computer program, which, when executed by a processor, implements all the steps of the method for restoring the erosion thickness and thermal history evolution process of the superimposed basin in the above-mentioned embodiment.

[0175] The various embodiments in this specification are described in a progressive manner. Similar portions between the various embodiments can be referenced to each other. Each embodiment focuses on the differences between the other embodiments. In particular, the hardware + program embodiments are generally similar to the method embodiments, so their description is relatively simple. For relevant portions, refer to the description of the method embodiments.

[0176] The foregoing description of this specification describes specific embodiments. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims can be performed in an order different from that described in the embodiments and still achieve the desired results. Furthermore, the processes depicted in the accompanying drawings do not necessarily require the specific order shown or the sequential order to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0177] Although the present application provides method operation steps such as embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-creative work. The order of steps listed in the embodiments is only one way of executing the steps among many steps and does not represent the only execution order. When an actual device or client product is executed, it can be executed in the order shown in the embodiments or the drawings or in parallel (for example, in a parallel processor or multi-threaded processing environment).

[0178] For the convenience of description, the above devices are described in terms of functions divided into various modules. Of course, when implementing the embodiments of this specification, the functions of each module can be implemented in the same or multiple software and / or hardware, or the modules that implement the same function can be implemented by a combination of multiple sub-modules or sub-units, etc. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, which can be electrical, mechanical or other forms.

[0179] Those skilled in the art will also appreciate that, in addition to implementing the controller in pure computer-readable program code, it is entirely possible to implement the same functionality by logically programming the method steps in the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, embedded microcontrollers, and the like. Therefore, such a controller can be considered a hardware component, and the devices included therein for implementing various functions can also be considered structures within the hardware component. Alternatively, the devices for implementing various functions can be considered both software modules implementing the method and structures within the hardware component.

[0180] In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0181] Memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. Memory is an example of a computer-readable medium.

[0182] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between the various embodiments can be referenced across them. Each embodiment focuses on the differences from the other embodiments. In particular, since the system embodiments are generally similar to the method embodiments, their description is relatively simple. For relevant parts, reference can be made to the description of the method embodiments. Throughout this specification, reference to the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with that embodiment or example are included in at least one embodiment or example of the embodiments in this specification. In this specification, the schematic representations of these terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, those skilled in the art may combine and integrate the different embodiments or examples, and features of different embodiments or examples, described in this specification, without conflict.

[0183] The above description is merely an example of the embodiments of this specification and is not intended to limit the embodiments of this specification. For those skilled in the art, various modifications and variations of the embodiments of this specification are possible. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of this specification shall be included within the scope of the claims of the embodiments of this specification.

Claims

1. A method for restoring the erosion thickness and thermal history evolution process of a superimposed basin, characterized in that: include: Selecting a rock sample from the target work area; wherein the rock sample is at least one of a sandstone sample containing apatite and zircon mineral particles, a volcanic rock sample, and a metamorphic rock sample; performing in-situ micro-area isotope dating on the rock sample to generate a dating result; The erosion thickness is determined based on the geothermal gradient of the target work area and the dating measurement results, and the thermal history evolution process is restored.

2. The method for restoring the erosion thickness and thermal history evolution process of a superimposed basin according to claim 1, characterized in that: The rock samples of the target work area are selected, including: The rock sample is selected from an overlying stratum or an underlying stratum of the eroded stratum in the target work area.

3. The method for restoring the erosion thickness and thermal history evolution process of a superimposed basin according to claim 2, characterized in that: The step of selecting rock samples from the target work area further includes: The rock samples are selected from the overlying strata or the underlying strata at depth intervals of 100 meters to 300 meters.

4. The method for restoring the erosion thickness and thermal history evolution process of a superimposed basin according to claim 1, characterized in that: The dating measurement results include: U-Th / He age of apatite and U-Th / He age of zircon.

5. The method for restoring the erosion thickness and thermal history evolution process of a superimposed basin according to claim 4, characterized in that: Also includes: The current formation temperature of the target work area is generated according to the geothermal gradient and the surface temperature of the target work area.

6. The method for restoring the erosion thickness and thermal history evolution process of a superimposed basin according to claim 5, characterized in that: Determining the erosion thickness according to the geothermal gradient of the target work area and the dating measurement result includes: Determine the depth corresponding to when the apatite enters the sealing temperature and the U-Th / He age when the apatite enters the sealing system according to the current formation temperature and the U-Th / He age of the apatite; Determine the depth corresponding to when the zircon enters the closed system and the U-Th / He age when the zircon enters the closed system based on the current formation temperature and the U-Th / He age of the zircon; The denudation thickness is determined according to the depth corresponding to when the apatite enters the closed temperature, the U-Th / He age when entering the closed system, the depth corresponding to when the zircon enters the closed temperature, and the U-Th / He age when entering the closed system.

7. The method for restoring the erosion thickness and thermal history evolution process of a superimposed basin according to claim 5, characterized in that: Determining the erosion thickness and restoring the thermal history evolution process according to the geothermal gradient of the target work area and the dating measurement results includes: determining a first relationship between the U-Th / He age of the apatite and the burial depth or the measured formation temperature; determining a second relationship between the U-Th / He age of the zircon and the depositional age of the rock sample; The thermal history evolution process is restored according to the first relationship and the second relationship.

8. A device for restoring the erosion thickness and thermal history evolution process of a superimposed basin, characterized in that: include: A rock sample selection module is used to select rock samples from a target work area; wherein the rock sample is at least one of a sandstone sample containing apatite and zircon mineral particles, a volcanic rock sample, and a metamorphic rock sample; a measurement result generating module, configured to perform in-situ micro-area isotope dating measurement on the rock sample to generate a dating measurement result; The erosion thickness determination module is used to determine the erosion thickness according to the geothermal gradient of the target work area and the dating measurement result, and to restore the thermal history evolution process.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the steps of the method for recovering the erosion thickness and thermal history evolution process of a superimposed basin as described in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for restoring the erosion thickness and thermal history evolution process of a superimposed basin as described in any one of claims 1 to 7 are implemented.

Images