Reconstruction method and device for thermal evolution history of sedimentary basin, equipment and medium

Abstract

The invention relates to a method, device, equipment and medium for reconstructing the thermal evolution history of a sedimentary basin, and relates to the field of thermal history simulation, and the method comprises the steps: obtaining the current temperature, stratigraphic stratification and vitrinite reflectivity of the basin, and the fission track age and (U-Th) / He age of a drilling sample; determining the median age of the fission track age of the sample; a comprehensive relation graph is obtained according to the median age, the (U-Th) / He age, stratigraphic stratification and the current temperature; performing forward and reverse modeling according to the comprehensive relation graph, the basin surface temperature, the sample burying temperature and the tectonic movement condition to obtain tectonic uplift time and a temperature change result; the stratum denudation amount is calculated in combination with the tectonic uplift time and the ancient ground temperature gradient; and reconstructing the foundation heat flux evolution history by combining stratum layering, vitrinite reflectivity, current temperature and stratum denudation amount. According to the reconstruction method provided by the invention, the stratum structure deformation time, denudation amount and paleo-heat flow change conditions can be accurately revealed, and reconstruction of the foundation heat flow evolution history is realized.

Description

Technical Field

[0001] The invention relates to the field of basin thermal history simulation, and in particular to a method and device, equipment and medium for reconstructing the thermal flow evolution history of a sedimentary basin. Background Art

[0002] The restoration of thermal history of sedimentary basins is of great significance for investigating the thermal history of basins and evaluating oil and gas resources.

[0003] Currently, there are two main methods for restoring the thermal history of sedimentary basins. One is the paleothermological method, which mainly relies on technologies such as vitrinite reflectance and low-temperature thermochronology. The other is the geodynamic method, which uses the current geothermal field as a constraint and integrates regional deformation mechanisms and erosion to restore early thermal information. Due to the lack of effective constraint data, this method has low accuracy.

[0004] Commonly used paleothermometers include organic paleothermometers (such as vitrinite reflectance, asphalt reflectance, etc.) and mineral paleothermometers (such as apatite and zircon fission tracks, apatite and zircon (U-Th) / He, etc.). The latter is also called low-temperature thermochronology paleothermometer. The low-temperature thermochronology paleothermometer can not only give the corresponding temperature of the sample history period through subsequent forward and inverse simulations, but also give the time when the temperature was reached in each period. It has become an important method for quantitative reconstruction of the thermal history of sedimentary basins and restoration of stratigraphic erosion in recent years. However, in the low-temperature thermochronology paleothermometer, the fission track test results will produce the ages of multiple provenance areas. Usually, P(χ 2 ) to verify the reliability of the test age, when the chi-square test P(χ 2 )>5%, indicating that the sample has passed the chi-square test, the data has a uniform component, and the median age is used for the test age. However, when P(χ2) of all samples is 0, it indicates that there are multiple provenance areas and a more reliable comprehensive age result cannot be produced.

[0005] CN110675497A discloses a method for restoring the thermal evolution history of superimposed basins. Starting from the perspective of the impact of basin superposition and transformation on the paleo-geothermal field information of superimposed basins, according to the differences in paleo-geothermal field information recording, preservation and later superposition and transformation at different evolutionary stages of superimposed basins, with the correct geological model and a large amount of actual geological data as constraints, a multi-disciplinary and multi-paleothermal scale method was used for comprehensive research, and a new idea and method for truly restoring the thermal evolution history of superimposed basins in different evolutionary stages was proposed.

[0006] CN118130539A discloses a method for obtaining the thermal history of complex lithology in ultra-deep layers of sedimentary basins. The method includes: based on the sedimentary burial history and with the current geothermal field as a constraint, using the measured data of vitrinite reflectance or equivalent vitrinite reflectance to determine the initial paleo-heat flow history of the ultra-deep layers of the target sedimentary basin; based on the initial paleo-heat flow history of the ultra-deep layers of the target sedimentary basin, combined with the measured data of isotopic temperature of carbonate rock clusters in the ultra-deep layers of the target sedimentary basin, fitting and correcting the carbonate rock thermal history path, determining the carbonate rock thermal history of the ultra-deep layers of the target sedimentary basin; based on the carbonate rock thermal history of the ultra-deep layers of the target sedimentary basin, combined with the measured data of zircon (U-Th) / He dating of the ultra-deep layers of the target sedimentary basin, determining the sandstone thermal history of the ultra-deep layers of the target sedimentary basin through forward and inverse coupled simulation; based on the carbonate rock thermal history and sandstone thermal history of the ultra-deep layers of the target sedimentary basin, determining the geothermal gradient and geothermal heat flow evolution process of the ultra-deep layers of the target sedimentary basin.

[0007] However, when restoring the thermal history of a basin based on paleothermal scales, there are still problems such as difficulty in processing thermal history sample data of sedimentary basins and low accuracy of thermal history simulation. Summary of the invention

[0008] In view of the problems existing in the prior art, the purpose of the present invention is to provide a method and device, equipment, and medium for reconstructing the heat flow evolution history of a sedimentary basin, so as to solve the problems of difficulty in processing sedimentary basin thermal history sample data and low thermal history simulation accuracy when restoring the basin thermal history.

[0009] To achieve this object, the present invention adopts the following technical solutions:

[0010] In a first aspect, the present invention provides a method for reconstructing the heat flow evolution history of a sedimentary basin, the reconstruction method comprising:

[0011] Obtain the current temperature data, stratigraphic layer data and vitrinite reflectance data of the target basin, as well as the fission track age and (U-Th) / He age of the drilling samples;

[0012] A chi-square test was performed on the fission track ages of the samples to determine the median age;

[0013] Based on the obtained median age, (U-Th) / He age, stratigraphic stratification data and present-day temperature data, a comprehensive relationship diagram of age, depth and temperature was obtained;

[0014] Based on the comprehensive relationship diagram of age, depth and temperature, the known surface temperature of the basin, the burial temperature of the sample and the tectonic movement, forward and inverse simulations were performed to obtain the results of tectonic uplift time and temperature changes;

[0015] Combined with the obtained tectonic uplift time and paleo-geothermal gradient, the amount of stratum denudation is calculated;

[0016] The evolution history of basement heat flow is reconstructed by combining stratigraphic stratification data, vitrinite reflectance data, current temperature data and stratigraphic erosion data.

[0017] The reconstruction method provided by the present invention can effectively process low-temperature thermochronology experimental data by adopting a specifically designed reconstruction process, and at the same time combine multiple low-temperature thermochronology methods to carry out multiple paleothermometric coupling simulations, so as to more accurately reveal the deformation time of stratigraphic structures, the amount of erosion and the changes in paleoheat flow, and realize the reconstruction of the evolution history of basement heat flow.

[0018] As a preferred technical solution of the present invention, if P(χ 2 )>5%, the test age was taken as the median age.

[0019] As a preferred technical solution of the present invention, if P(χ 2 )≤5%, the fission track ages of the samples are regrouped and analyzed to obtain the age distribution probability histogram of the fission track ages of the samples, and the youngest age group is selected as the median age.

[0020] As a preferred technical solution of the present invention, the age, depth and temperature comprehensive relationship diagram includes: the relationship between the median age and the formation age and the formation temperature of the sample.

[0021] As a preferred technical solution of the present invention, the forward and inversion simulation includes: placing the sample in a partial annealing zone or a He diffusion temperature range according to the measured age, track length and Dpar value, and then transforming the tectonic movement time and the temperature of the corresponding period to obtain the degree of fit GOF. If the degree of fit GOF ≥ 0.05, the simulation result meets the requirements.

[0022] Preferably, in the forward and inverse simulations, samples with a test age < formation age are selected to carry out thermal history simulation of fission track age and (U-Th) / He age.

[0023] As a preferred technical solution of the present invention, the calculation formula of the stratum erosion amount is as follows:

[0024] H=(T1-T0) / (ΔT / ΔZ)

[0025] Where H is the amount of stratum denudation, m; T1 is the paleotemperature at the unconformity surface, °C; T0 is the surface temperature, °C; ΔT / ΔZ is the paleotemperature gradient, km / °C, and is the slope of the paleotemperature fitting line.

[0026] As a preferred technical solution of the present invention, the reconstruction of the basement heat flow evolution history includes: setting basin simulation formation parameters based on the amount of formation erosion and the stratification data in the seismic data, and using the vitrinite reflectance and the current temperature data as constraints, and adjusting the basement heat flow to obtain the simulated vitrinite reflectance value and the current temperature value to simulate and establish the formation burial history. If the current simulation result is consistent with the measured result, it indicates that the simulation result is reliable. Then, based on the formation burial history, the ancient basement heat flow of a certain period in the past is calculated, and then combined with the formation stratification data, the basement heat flow evolution history of the single well area is reconstructed.

[0027] In a second aspect, the present invention provides a device for reconstructing the heat flow evolution history of a sedimentary basin, the reconstruction device comprising:

[0028] Parameter acquisition module, used to obtain the current temperature data, stratigraphic layering data and vitrinite reflectance data of the target basin, as well as the fission track age and (U-Th) / He age of drilling samples;

[0029] The median age acquisition module is used to perform a chi-square test on the fission track age of the sample to determine the median age;

[0030] A comprehensive relationship diagram acquisition module is used to obtain a comprehensive relationship diagram of age, depth and temperature based on the obtained median age, (U-Th) / He age, combined with stratigraphic layer data and current temperature data;

[0031] The forward and inversion module is used to perform forward and inversion simulations based on the comprehensive relationship diagram of age, depth and temperature, the known surface temperature of the basin, the burial temperature of the sample and the tectonic movement to obtain the results of tectonic uplift time and temperature change;

[0032] The stratigraphic erosion acquisition module is used to calculate the stratigraphic erosion amount by combining the obtained tectonic uplift time and paleo-geothermal gradient;

[0033] The reconstruction module is used to reconstruct the evolution history of basement heat flow by combining stratigraphic stratification data, vitrinite reflectance data, current temperature data and stratigraphic erosion data.

[0034] In a third aspect, the present invention provides an electronic device, the electronic device comprising:

[0035] at least one processor; and a memory communicatively coupled to the at least one processor;

[0036] Wherein, the memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for reconstructing the thermal flow evolution history of a sedimentary basin described in the first aspect.

[0037] In a fourth aspect, the present invention provides a computer storage medium storing computer executable instructions, which, when executed by a processor, implement the method for reconstructing the thermal flow evolution history of a sedimentary basin as described in the first aspect.

[0038] Compared with the prior art solutions, the present invention has the following beneficial effects:

[0039] The reconstruction method provided by the present invention can effectively process low-temperature thermochronology experimental data by adopting a specifically designed reconstruction process, and at the same time combine multiple low-temperature thermochronology methods to carry out multiple paleothermometric coupling simulations, so as to more accurately reveal the deformation time of stratigraphic structures, the amount of erosion and the changes in paleoheat flow, and realize the reconstruction of the evolution history of basement heat flow. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 is a flow chart of a method for reconstructing the heat flow evolution history of a sedimentary basin provided by an embodiment of the present invention;

[0041] Figure 2 is a schematic diagram of a device for reconstructing the heat flow evolution history of a sedimentary basin provided by an embodiment of the present invention;

[0042] Figure 3 is a schematic diagram of an electronic device provided by an embodiment of the present invention;

[0043] Figure 4 is a flow chart of a method for reconstructing the heat flow evolution history of a sedimentary basin in Example 1 of the present invention;

[0044] Figure 5 This is a chi-square test diagram of the fission track age of the drilling sample at a depth of 1710 m in Example 1 of the present invention;

[0045] Figure 6 is a chi-square test diagram of the fission track age of the drilling sample at a depth of 1950 m in Example 1 of the present invention;

[0046] Figure 7 This is a chi-square test diagram of the fission track age of the drilling sample at a depth of 2250 m in Example 1 of the present invention;

[0047] Figure 8 This is a chi-square test diagram of the fission track age of the drilling sample at a depth of 2825 m in Example 1 of the present invention;

[0048] Fig. 9 This is a chi-square test diagram of the fission track age of the drilling sample at a depth of 2890 m in Example 1 of the present invention;

[0049] Fig.10 is a fission track age distribution histogram of the drilling sample at a depth of 1710 m in Example 1 of the present invention;

[0050] Fig.11 is a fission track age distribution histogram of the drilling sample at a depth of 1950 m in Example 1 of the present invention;

[0051] Fig.12 is a fission track age distribution histogram of the drilling sample at a depth of 2250 m in Example 1 of the present invention;

[0052] Fig.13 is a fission track age distribution histogram of the drilling sample at a depth of 2825 m in Example 1 of the present invention;

[0053] Fig.14 is a fission track age distribution histogram of the drilling sample at a depth of 2890 m in Example 1 of the present invention;

[0054] Fig.15 This is the AFT+Ahe age simulation result diagram of the drilling sample at a depth of 2250m in Example 1 of the present invention;

[0055] Fig.16 This is the AFT age simulation result diagram of the drilling sample at a depth of 2825m in Example 1 of the present invention;

[0056] Fig.17 is a diagram of the evolution history of the paleo-geothermal gradient of the borehole in the sedimentary basin in Example 1 of the present invention;

[0057] Fig.18 is a schematic diagram of the evolution of paleo-geothermal gradient in Example 1 of the present invention;

[0058] Fig.19 is a schematic diagram of the sedimentation and burial history in Example 1 of the present invention;

[0059] Fig. 20 It is a schematic diagram of the thermal flow evolution history of the substrate in a single well area in Example 1 of the present invention.

[0060] In the figure: 100-parameter acquisition module, 200-median age acquisition module, 300-comprehensive relationship diagram acquisition module, 400-forward and inversion module, 500-stratum erosion acquisition module, 600-reconstruction module;

[0061] 10 - electronic device, 11 - processor, 12 - ROM, 13 - RAM, 14 - bus, 15 - I / O interface, 16 - input unit, 17 - output unit, 18 - storage unit, 19 - communication unit.

[0062] The present invention is further described in detail below. However, the following examples are only simplified examples of the present invention and do not represent or limit the scope of protection of the present invention. The scope of protection of the present invention shall be subject to the claims. DETAILED DESCRIPTION

[0063] To better illustrate the present invention and facilitate understanding of the technical solution of the present invention, typical but non-limiting embodiments of the present invention are as follows:

[0064] This embodiment provides a method for reconstructing the heat flow evolution history of a sedimentary basin. The process is as follows: Figure 1 As shown, the reconstruction method includes:

[0065] Obtain the current temperature data, stratigraphic layer data and vitrinite reflectance data of the target basin, as well as the fission track age and (U-Th) / He age of the drilling samples;

[0066] A chi-square test was performed on the fission track ages of the samples to determine the median age;

[0067] Based on the obtained median age, (U-Th) / He age, stratigraphic stratification data and present-day temperature data, a comprehensive relationship diagram of age, depth and temperature was obtained;

[0068] Based on the comprehensive relationship diagram of age, depth and temperature, the known surface temperature of the basin, the burial temperature of the sample and the tectonic movement, forward and inverse simulations were performed to obtain the results of tectonic uplift time and temperature changes;

[0069] Combined with the obtained tectonic uplift time and paleo-geothermal gradient, the amount of stratum denudation is calculated;

[0070] The evolution history of basement heat flow is reconstructed by combining stratigraphic stratification data, vitrinite reflectance, present-day temperature data and stratigraphic erosion amount.

[0071] In the present invention, the relevant samples tested are samples at various depths obtained along the drilling depth direction. Fission track age and (U-Th) / He age tests can be performed on samples at various depths. Alternatively, the corresponding drilling depth interval can be selected according to the test advantages of the corresponding fission track age and (U-Th) / He age, and the corresponding fission track age and (U-Th) / He age tests can be performed.

[0072] Wherein, in the chi-square test, if P(χ 2 )>5%, the test age was taken as the median age.

[0073] Wherein, in the chi-square test, if P(χ 2 )≤5%, the fission track ages of the samples are regrouped and analyzed to obtain the age distribution probability histogram of the fission track ages of the samples, and the youngest age group is selected as the median age.

[0074] In the present invention, the youngest age group refers to the age group corresponding to the maximum probability in the probability histogram.

[0075] In the present invention, the fission track age, (U-Th) / He age, and age distribution probability histogram can be obtained by processing using conventional data processing software in the art, such as the IsoplotR platform.

[0076] In the present invention, the stratigraphic layering data is obtained through interpretation of oilfield seismic data, and can be obtained specifically according to conventional technical means in the art.

[0077] In the present invention, the current temperature data is obtained through drilling while drilling test, which can be obtained specifically according to conventional technical means in the field.

[0078] In the present invention, the process of obtaining a comprehensive relationship diagram of age, depth and temperature based on the obtained median age, (U-Th) / He age combined with stratigraphic stratification data and current temperature data is exemplified as follows: the lower X-axis in the image is the age range of the low-temperature thermochronology test; the left Y-axis is the stratigraphic depth, ranging from surface depth to drilling bottom depth; the right Y-axis is the stratigraphic temperature, ranging from surface temperature to drilling bottom temperature; for the apatite fission track, the track length can also be added to the upper Y-axis, and the figure can display the stratigraphic temperature corresponding to the test sample, the direct correspondence between the stratigraphic temperature and depth, and the variation characteristics of the age / apatite track length with depth.

[0079] The age, depth and temperature comprehensive relationship diagram includes: the relationship between the median age and the formation age and the formation temperature of the sample.

[0080] In the present invention, the relationship between the test age of each sample in the basin and the formation age, the formation temperature of the sample and other information can be analyzed through the comprehensive relationship diagram of age, depth and temperature. In the diagram, if the sample test age is within the formation age line, it means that the sample test age is less than the formation deposition age, and a complete annealing or He diffusion process may have occurred; if the sample test age is above the formation age line, it means that the sample test age is equivalent to the bottom layer age and may have undergone partial annealing or He diffusion; if the sample test age is outside the formation age line, it means that the sample test age is greater than the formation deposition age, and the sample may come from multiple source areas.

[0081] In the present invention, the forward and inversion simulation can be implemented using commonly used software in the field such as HeFTy software.

[0082] Among them, the forward and inverse simulation includes: placing the sample in a partial annealing zone or He diffusion temperature range according to the measured age, track length and Dpar value, and then transforming the tectonic movement time and the temperature of the corresponding period to obtain the GOF. If the GOF is ≥ 0.05, the simulation result meets the requirements.

[0083] Among them, in the forward and inverse simulations, samples with test age < formation age are selected to carry out thermal history simulation of fission track age and (U-Th) / He age.

[0084] The calculation formula of the stratum erosion amount is as follows:

[0085] H=(T1-T0) / (ΔT / ΔZ)

[0086] Where H is the amount of stratum denudation, m; T1 is the paleotemperature at the unconformity surface, °C; T0 is the surface temperature, °C; ΔT / ΔZ is the paleotemperature gradient, km / °C, and is the slope of the paleotemperature fitting line.

[0087] In the present invention, the paleo-geothermal gradient change used in the calculation of the stratum denudation amount can be obtained based on the current temperature data of the sample and the paleo-geothermal evolution, and can be obtained specifically according to the conventional acquisition process in the art.

[0088] Among them, the reconstruction of the basement heat flow evolution history includes: setting basin simulation formation parameters based on the amount of formation erosion and the stratification data in the seismic data, and using the vitrinite reflectance data and the current temperature data as constraints, and adjusting the basement heat flow to obtain the simulated vitrinite reflectance value and the current temperature value to simulate and establish the formation burial history. If the current simulation result is consistent with the measured result, it indicates that the simulation result is reliable. Then, based on the formation burial history, the ancient basement heat flow of a certain period in the past is calculated, and then combined with the formation stratification data, the basement heat flow evolution history of the single well area is reconstructed.

[0089] In the present invention, reconstruction of the evolution history of basement heat flow can be performed with the aid of commonly used simulation software in the art, such as PetroMod basin simulation software.

[0090] In the present invention, when simulating the burial history of the strata by adjusting the basement heat flow to obtain the simulated vitrinite reflectance value and the current temperature value, the basement heat flow during the tectonic activity period in the thermal history simulation can be preferably used for adjustment.

[0091] Furthermore, this embodiment provides a device for reconstructing the heat flow evolution history of a sedimentary basin, such as Figure 2 As shown, the reconstruction device comprises:

[0092] The parameter acquisition module 100 is used to obtain the current temperature data, stratigraphic layering data and vitrinite reflectance data of the target basin, and the fission track age and (U-Th) / He age of the drilling samples;

[0093] The median age acquisition module 200 is used to perform a chi-square test on the fission track age of the sample to determine the median age;

[0094] A comprehensive relationship diagram acquisition module 300 is used to obtain a comprehensive relationship diagram of age, depth and temperature based on the obtained median age, (U-Th) / He age, combined with the stratigraphic layer data and the current temperature data;

[0095] The forward and inversion module 400 is used to perform forward and inversion simulations based on the age, depth and temperature comprehensive relationship diagram, the known surface temperature of the basin, the sample burial temperature and the tectonic movement conditions to obtain the results of tectonic uplift time and temperature change;

[0096] The stratum denudation amount acquisition module 500 is used to calculate the stratum denudation amount by combining the obtained tectonic uplift time and paleo-geothermal gradient;

[0097] The reconstruction module 600 is used to reconstruct the evolution history of basement heat flow by combining stratigraphic layering data, vitrinite reflectance, current temperature data and stratigraphic erosion amount.

[0098] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0099] Further, the present invention provides an electronic device, which is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device may also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0100] like Figure 3 As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. The I / O interface 15 is also connected to the bus 14.

[0101] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0102] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the aforementioned method for reconstructing the heat flow evolution history of the sedimentary basin.

[0103] In some embodiments, the aforementioned method for reconstructing the heat flow evolution history of the sedimentary basin can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the method for reconstructing the heat flow evolution history of the aforementioned sedimentary basin described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to execute the aforementioned method for reconstructing the heat flow evolution history of the sedimentary basin in any other appropriate manner (for example, by means of firmware).

[0104] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0105] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0106] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in combination with an instruction execution system, device or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0107] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0108] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0109] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0110] The server provided in this embodiment includes: a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the method for reconstructing the thermal flow evolution history of the aforementioned sedimentary basin is implemented.

[0111] Unless otherwise specifically stated, terms such as processing, computing, calculating, determining, displaying, etc. may refer to the actions and / or processes of one or more processing or computing systems, or similar devices, which operate and convert data represented as physical (e.g., electronic) quantities within registers or memories of a processing system into other data similarly represented as physical quantities within memories, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals may be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0112] Those skilled in the art will also appreciate that the various illustrative logic blocks, modules, circuits, and algorithm steps described in conjunction with the embodiments of the present invention may all be implemented as electronic hardware, computer software, or a combination thereof. In order to clearly illustrate the interchangeability between hardware and software, various illustrative components, blocks, modules, circuits, and steps are generally described above around their functions. Whether such functions are implemented as hardware or software depends on the specific application and the design constraints imposed on the entire system. A skilled person may implement the described functions in an alternative manner for each specific application, but such implementation decisions should not be interpreted as departing from the scope of protection of the present invention.

[0113] The steps of the method or algorithm described in conjunction with the embodiments herein may be directly embodied as hardware, a software module executed by a processor, or a combination thereof. The software module may be located in a RAM memory, a flash memory, a ROM memory, an EPROM memory, an EEPROM memory, a register, a hard disk, a mobile disk, a CD-ROM, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor so that the processor can read information from the storage medium and can write information to the storage medium. Of course, the storage medium may also be an integral part of the processor. The processor and the storage medium may be located in an ASIC. The ASIC may be located in a user terminal. Of course, the processor and the storage medium may also be present in a user terminal as discrete components.

[0114] For software implementation, the techniques described in the present invention can be implemented with modules (e.g., procedures, functions, etc.) that perform the functions described in the present application. These software codes can be stored in a memory unit and executed by a processor. The memory unit can be implemented within the processor or outside the processor. In the latter case, it is coupled to the processor in a communication manner via various means, which are well known in the art.

[0115] Furthermore, in order to illustrate the reconstruction effect that can be achieved by the reconstruction method of the heat flow evolution history of a sedimentary basin provided by the present invention, an actual example is used for illustration, as follows:

[0116] Example 1

[0117] This embodiment provides a method for reconstructing the heat flow evolution history of an actual basin. The process is as follows: Figure 4 As shown, a basin in the eastern sea area is taken as the research object.

[0118] ① Collect the basin's current temperature data, stratigraphic stratification data, and vitrinite reflectance data, select borehole samples in the sedimentary basin, and conduct experimental tests using a low-temperature thermochronology test program.

[0119] ② The low-temperature thermochronological age analysis of the five fission track test results was carried out using the Isoplot experimental data processing platform, and the corresponding chi-square test P (χ 2 ) value, the result is as follows Figure 5 , Figure 6 , Figure 7 , Figure 8 and Fig. 9 As shown in the figure, it can be seen that the chi-square test value P(χ 2 ≤5%, that is, it is necessary to conduct regrouping analysis; use the IsoplotR experimental data processing platform to obtain the age distribution probability histogram and sort the youngest age group. The results are as follows Fig.10 , Fig.11 , Fig.12 , Fig.13 , Fig.14 shown.

[0120] ③ Based on the processed results, a comprehensive relationship diagram of single well age, depth and temperature is established. Fig.15 and Fig.16 As shown in the figure, it can be seen that the apatite fission track age gradually decreases with increasing depth. Among them, 3 samples are within the stratigraphic age line, indicating that they have undergone annealing. The apatite (U-Th) / He age of 6 particles is less than the stratigraphic age, indicating that He diffusion occurred in their historical period.

[0121] ④ In the above samples, the particles with test age younger than the formation age are selected to carry out thermal history simulation with the help of apatite fission track age and (U-Th) / He age using HeFTy software. For example, the apatite fission track age and (U-Th) / He age of the sample at 2250m are used for coupled simulation to jointly reveal the thermal history of the sample and obtain the time of tectonic uplift and temperature change. The results are as follows: Fig.17 As shown; further, according to the current temperature data of the sample and the paleo-geothermal evolution, the paleo-geothermal gradient change of the sample is obtained, and the results are shown as follows Fig.18 shown.

[0122] ⑤ Combined with the results of thermal history simulation, the amount of stratum erosion was calculated according to the formula. The calculation results show that the sample underwent tectonic movement from 13Ma to 7Ma, and the amount of tectonic erosion was 280-340m.

[0123] ⑥ Based on the reconstruction results of the denudation volume of thermal history simulation and the stratification data in the seismic data, the basin simulation stratigraphic parameters were set through the PetroMod software. The vitrinite reflectance data and the current temperature data were used as constraints. The basal heat flow was adjusted to obtain the simulated vitrinite reflectance value and the current temperature value to simulate the stratigraphic burial history. The results are as follows: Fig.19 shown.

[0124] ⑦ Based on the reconstruction results of the burial history, the paleobasement heat flow in each key period in the past was calculated, and combined with the stratigraphic stratification data, the evolution history of the basement heat flow in the single well area was revealed. The results are as follows: Fig. 20 shown.

[0125] The preferred embodiments of the present invention are described in detail above. However, the present invention is not limited to the specific details in the above embodiments. Within the technical concept of the present invention, a variety of simple modifications can be made to the technical solution of the present invention, and these simple modifications all belong to the protection scope of the present invention.

[0126] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any suitable manner without contradiction. In order to avoid unnecessary repetition, the present invention will not further describe various possible combinations.

[0127] In addition, various embodiments of the present invention may be arbitrarily combined, and as long as they do not violate the concept of the present invention, they should also be regarded as the contents disclosed by the present invention.

Claims

1. A method for reconstructing the heat flow evolution history of a sedimentary basin, characterized in that: The reconstruction method comprises: Obtain the current temperature data, stratigraphic layer data and vitrinite reflectance data of the target basin, as well as the fission track age and (U-Th) / He age of the drilling samples; A chi-square test was performed on the fission track ages of the samples to determine the median age; Based on the obtained median age, (U-Th) / He age, stratigraphic stratification data and present-day temperature data, a comprehensive relationship diagram of age, depth and temperature was obtained; Based on the comprehensive relationship diagram of age, depth and temperature, the known surface temperature of the basin, the burial temperature of the sample and the tectonic movement, forward and inverse simulations were performed to obtain the results of tectonic uplift time and temperature changes; Combined with the obtained tectonic uplift time and paleo-geothermal gradient, the amount of stratum denudation is calculated; The evolution history of basement heat flow is reconstructed by combining stratigraphic stratification data, vitrinite reflectance data, current temperature data and stratigraphic erosion data.

2. The method for reconstructing the heat flow evolution history of a sedimentary basin according to claim 1, characterized in that: In the chi-square test, if P(χ 2 )>5%, the test age was taken as the median age.

3. The method for reconstructing the heat flow evolution history of a sedimentary basin according to claim 1, characterized in that: In the chi-square test, if P(χ 2 )≤5%, the fission track ages of the samples are regrouped and analyzed to obtain the age distribution probability histogram of the fission track ages of the samples, and the youngest age group is selected as the median age.

4. The method for reconstructing the heat flow evolution history of a sedimentary basin according to claim 1, characterized in that: The age, depth and temperature comprehensive relationship diagram includes: the relationship between the median age and the formation age and the formation temperature where the sample is located.

5. The method for reconstructing the heat flow evolution history of a sedimentary basin according to claim 1, characterized in that: The forward and inversion simulation includes: placing the sample in a partial annealing zone or a He diffusion temperature interval according to the measured age, track length and Dpar value, and then transforming the tectonic movement time and the temperature of the corresponding period to obtain the GOF. If the GOF is ≥ 0.05, the simulation result meets the requirements; Preferably, in the forward and inverse simulations, samples with a test age < formation age are selected to carry out thermal history simulation of fission track age and (U-Th) / He age.

6. The method for reconstructing the heat flow evolution history of a sedimentary basin according to claim 1, characterized in that: The calculation formula of the stratum erosion amount is as follows: H=(T1-T0) / (ΔT / ΔZ) Where, H is the amount of stratum erosion, m; T1 is the paleotemperature at the unconformity surface, ℃; T0 is the surface temperature, ℃; ΔT / ΔZ is the paleotemperature gradient, km / ℃, and is the slope of the paleotemperature fitting line.

7. The method for reconstructing the heat flow evolution history of a sedimentary basin according to claim 1, characterized in that: The reconstruction of the basement heat flow evolution history includes: setting basin simulation formation parameters based on the amount of formation erosion and the stratification data in the seismic data, and using the vitrinite reflectance data and the current temperature data as constraints, and simulating and establishing the formation burial history by adjusting the basement heat flow to obtain the simulated vitrinite reflectance value and the current temperature value. If the current simulation result is consistent with the measured result, it indicates that the simulation result is reliable. Then, based on the formation burial history, the ancient basement heat flow of a certain period in the past is calculated, and then combined with the formation stratification data, the basement heat flow evolution history of the single well area is reconstructed.

8. A device for reconstructing the heat flow evolution history of a sedimentary basin, characterized in that: The reconstruction device comprises: Parameter acquisition module, used to obtain the current temperature data, stratigraphic layering data and vitrinite reflectance data of the target basin, as well as the fission track age and (U-Th) / He age of drilling samples; The median age acquisition module is used to perform a chi-square test on the fission track age of the sample to determine the median age; A comprehensive relationship diagram acquisition module is used to obtain a comprehensive relationship diagram of age, depth and temperature based on the obtained median age, (U-Th) / He age, combined with stratigraphic layer data and current temperature data; The forward and inversion module is used to perform forward and inversion simulations based on the comprehensive relationship diagram of age, depth and temperature, the known surface temperature of the basin, the burial temperature of the sample and the tectonic movement to obtain the results of tectonic uplift time and temperature change; The stratigraphic erosion acquisition module is used to calculate the stratigraphic erosion amount by combining the obtained tectonic uplift time and paleo-geothermal gradient; The reconstruction module is used to reconstruct the evolution history of basement heat flow by combining stratigraphic stratification data, vitrinite reflectance data, current temperature data and stratigraphic erosion data.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively coupled to the at least one processor; Wherein, the memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for reconstructing the thermal flow evolution history of a sedimentary basin according to any one of claims 1-7.

10. A computer storage medium, characterized in that: The computer storage medium stores computer executable instructions, and when the computer executable instructions are executed by a processor, the method for reconstructing the heat flow evolution history of a sedimentary basin according to any one of claims 1 to 7 is implemented.

Images