Low-temperature chronology data processing method, system and device based on paleogeothermal gradient

Abstract

The present invention relates to a method, system and device for processing low-temperature age data based on paleogeothermal gradient. The method includes obtaining M samples and sample parameters based on the strata of sedimentary basins above angular unconformities; solving the burial history time of corresponding samples by using the depth-time equation for the sample parameters; taking the burial history time as the stratigraphic boundary and calculating the depth of the samples by using the depth calculation formula to obtain the burial depth of the corresponding samples, and constructing the burial history curve of the samples according to the burial depth and time of the samples; assuming N paleogeothermal gradient-time curves and constructing a geothermal gradient-time function through the N paleogeothermal gradient-time curves, and converting the burial history curve of the corresponding samples into N time-temperature curves of the samples by using the burial depth of the samples and the geothermal gradient-time function. The geothermal gradient-time curve is obtained by calculating the goodness of fit between the simulation result and the measured result through the time-temperature curve, so as to improve the accuracy of restoring the paleogeothermal gradient of sedimentary basins.

Description

Technical Field

[0001] The present invention relates to the field of geological technologies, and in particular, to a method, system and device for processing low-temperature age data based on paleogeothermal gradient. Background Art

[0002] The paleogeothermal gradient of sedimentary basins is an essential parameter in the fields of geological research, petroleum exploration, etc. The restoration of this important parameter has always been one of the focus issues in the research of related fields. At present, generally, data such as vitrinite reflectance and apatite fission track are used to restore the paleogeothermal gradient. The common problem of such methods is the use of a single type of data. However, in the actual exploration and research field, it is difficult to obtain a large amount of unified type data for a certain drilling profile; in addition, each single type of data has its own advantages and disadvantages. For example, vitrinite reflectance is not sensitive to time, and the sensitive temperature range of apatite fission track is limited to the annealing interval (about 60 - 120 °C), etc. Summary of the Invention

[0003] Embodiments of the present invention provide a method, system and device for processing low-temperature age data based on paleogeothermal gradient, which are used to solve the technical problems that the existing restoration of the paleogeothermal gradient of sedimentary basins by using a single type of vitrinite reflectance and apatite fission track is not sensitive to time and temperature, and the obtained simulation output results are inaccurate.

[0004] In order to achieve the above object, the embodiments of the present invention provide the following technical solutions:

[0005] A method for processing low-temperature age data based on paleogeothermal gradient includes the following steps:

[0006] Based on the strata of the sedimentary basin above the angular unconformity, obtain M samples of the strata of the sedimentary basin above the angular unconformity and the parameters of each sample;

[0007] Based on the parameters of each sample, use the depth-time equation to solve the burial history time corresponding to the sample;

[0008] Take the burial history time of each sample as the stratigraphic boundary, and use the depth calculation formula based on the principle of equal formation mass at different time points to calculate the burial depth corresponding to the sample, and construct the burial history curve corresponding to the sample based on the burial depth and time of the sample;

[0009] Assume N paleogeothermal gradient-time curves and construct a geothermal gradient-time function through the N paleogeothermal gradient-time curves. According to the burial depth of the sample and the geothermal gradient-time function, convert the burial history curve of the corresponding sample into N time-temperature curves of the sample;

[0010] Among them, N is a natural number greater than 500, and M is a natural number greater than 1.

[0011] Preferably, the low-temperature chronology data processing method based on paleogeothermal gradient includes:

[0012] According to M samples, M×N time-temperature curves of the sedimentary basin strata are obtained. The time-temperature curves of the same geothermal gradient-time function are extracted from the M×N time-temperature curves as a set of thermal history simulation curves, and then N sets of thermal history simulation curves are obtained;

[0013] Using the low-temperature chronology method to simulate the age of N sets of the thermal history simulation curves, and obtaining the goodness-of-fit values of each set of the thermal history simulation curves, obtaining N goodness-of-fit values;

[0014] Comparing each of the goodness-of-fit values with a goodness-of-fit threshold, screening out K sets of thermal history simulation curves corresponding to the qualified goodness-of-fit values from the N goodness-of-fit values, and obtaining the corresponding K paleogeothermal gradient-time curves according to the K sets of thermal history simulation curves;

[0015] At the same time point, taking the mean value of each data on the K paleogeothermal gradient-time curves to obtain the temperature mean value of each data; connecting the temperature mean values of all the data in chronological order into a curve, and this curve is used as the thermal history simulation result of the sedimentary basin strata above the angular unconformity;

[0016] Among them, a set of thermal history simulation curves corresponds to a paleogeothermal gradient-time curve.

[0017] Preferably, taking the mean value of each data on the K paleogeothermal gradient-time curves includes:

[0018] Dividing each of the K paleogeothermal gradient-time curves into A equal parts, taking the mean value of the temperatures at the same time point in all the K paleogeothermal gradient-time curves after dividing into A equal parts, and obtaining A temperature mean values;

[0019] Connecting the A temperature mean values in chronological order into a curve, and this curve is used as the thermal history simulation result of the sedimentary basin strata above the angular unconformity;

[0020] Among them, A is a natural number greater than 50.

[0021] Preferably, the depth-time equation is:

[0022]

[0023] In the formula, h d is the depth at the lower part of the stratum where the sample is located, t d is the time corresponding to the depth at the lower part of the stratum where the sample is located, h uis the depth of the upper part of the formation where the sample is located, t u is the time corresponding to the depth of the upper part of the formation where the sample is located, h s is the depth of the sample, t s is the burial history time corresponding to the depth of the sample, y is the porosity of the sample, is the porosity function.

[0024] Preferably, the depth calculation formula is:

[0025]

[0026] In the formula, y is the porosity of the sample, is the porosity function, h is the depth of the sample, j is the serial number of the formation boundary, i is the serial number of the time point corresponding to the formation boundary, and both i and j are natural numbers greater than 0.

[0027] Preferably, the low-temperature chronology data processing method based on paleogeothermal gradient includes: constructing a geothermal gradient-time function for N paleogeothermal gradient-time curves by using a linear interpolation function.

[0028] The present invention also provides a low-temperature chronology data processing system based on paleogeothermal gradient, including a data acquisition module, a first processing module, a second processing module, and a third processing module;

[0029] The data acquisition module is used to obtain M samples located in the formation of the sedimentary basin above the angular unconformity and the parameters of each sample based on the formation of the sedimentary basin above the angular unconformity;

[0030] The first processing module is used to solve the burial history time corresponding to the sample by using a depth-time equation based on the parameters of each sample;

[0031] The second processing module is used to use the burial history time of each sample as the formation boundary, and calculate the depth of the corresponding sample by using a depth calculation formula based on the principle of equal formation mass at different time points, obtain the burial depth of the corresponding sample, and construct a burial history curve of the corresponding sample based on the burial depth and time of the sample;

[0032] The third processing module is used to assume N paleogeothermal gradient-time curves, construct a geothermal gradient-time function through the N paleogeothermal gradient-time curves, and convert the burial history curve of the corresponding sample into N time-temperature curves of the sample according to the burial depth of the sample and the geothermal gradient-time function;

[0033] The depth-time equation is:

[0034]

[0035] In the formula, hd is the depth of the lower part of the formation where the sample is located, t d is the time corresponding to the depth of the lower part of the formation where the sample is located, h u is the depth of the upper part of the formation where the sample is located, t u is the time corresponding to the depth of the upper part of the formation where the sample is located, h s is the depth of the sample, t s is the burial history time corresponding to the sample depth, y is the porosity of the sample, is the porosity function;

[0036] The depth calculation formula is as follows:

[0037]

[0038] In the formula, y is the porosity of the sample, is the porosity function, h is the depth of the sample, j is the serial number of the formation boundary, i is the serial number of the time point corresponding to the formation boundary, and both i and j are natural numbers greater than 0;

[0039] Among them, N is a natural number greater than 500, and M is a natural number greater than 1.

[0040] Preferably, the low-temperature chronology data processing system based on the paleogeothermal gradient includes an output module. The output module is used to obtain M×N time-temperature curves of the sedimentary basin formation according to M samples, extract the time-temperature curves of the same geothermal gradient-time function from the M×N time-temperature curves as a set of thermal history simulation curves, and then obtain N sets of thermal history simulation curves; use the low-temperature chronology method to simulate the age of each set of the N sets of thermal history simulation curves to calculate the goodness-of-fit value of each set of the thermal history simulation curves, and obtain N goodness-of-fit values; compare each goodness-of-fit value with the goodness-of-fit threshold, screen out K sets of thermal history simulation curves corresponding to the qualified goodness-of-fit values from the N goodness-of-fit values, and obtain the corresponding K paleogeothermal gradient-time curves according to the K sets of thermal history simulation curves; at the same time point, calculate the mean value of each data on the K paleogeothermal gradient-time curves to obtain the temperature mean value of each data; connect the temperature mean values of all the data in chronological order into a curve, and this curve is used as the thermal history simulation result of the sedimentary basin formation above the angular unconformity; among them, a set of thermal history simulation curves corresponds to a paleogeothermal gradient-time curve.

[0041] Preferably, the output module is further used to divide each paleogeothermal gradient-time curve into A equal parts, calculate the mean value of the temperatures at the same time point in all the paleogeothermal gradient-time curves after dividing them into A equal parts to obtain A temperature mean values, and connect the A temperature mean values in chronological order into a curve, and this curve is used as the thermal history simulation result of the sedimentary basin formation above the angular unconformity; among them, A is a natural number greater than 50.

[0042] The present invention also provides a low-temperature age data processing device based on paleogeothermal gradient, including a processor and a memory;

[0043] The memory is used for storing program codes and transmitting the program codes to the processor;

[0044] The processor is used for executing the above-mentioned low-temperature age data processing method based on paleogeothermal gradient according to the instructions in the program codes.

[0045] It can be seen from the above technical solutions that the embodiments of the present invention have the following advantages: For the low-temperature age data processing method, system and device based on paleogeothermal gradient, the method includes: obtaining M samples located in the strata of the sedimentary basin above the angular unconformity and the parameters of each sample based on the strata of the sedimentary basin above the angular unconformity; solving the burial history time of the corresponding sample by using the depth-time equation based on the parameters of each sample; taking the burial history time of each sample as the stratigraphic boundary, calculating the burial depth of the corresponding sample by using the depth calculation formula of the principle of equal formation mass at different time points, and constructing the burial history curve of the corresponding sample based on the burial depth and time of the sample; assuming N paleogeothermal gradient-time curves and constructing a geothermal gradient-time function through the N paleogeothermal gradient-time curves, and converting the burial history curve of the corresponding sample into the time-temperature curve of the sample according to the burial depth of the sample and the geothermal gradient-time function. The time-temperature curve obtained by the low-temperature age data processing method based on paleogeothermal gradient avoids the phenomenon of being insensitive to time and temperature. Using this time-temperature curve as the data for simulating the thermal history of sample low-temperature chronology, the goodness of fit between the simulation result and the measured result is calculated, and finally the geothermal gradient-time curve is obtained, improving the accuracy of restoring the paleogeothermal gradient of the sedimentary basin and also improving the accuracy of the thermal history simulation result, and solving the technical problem that the existing restoration of the paleogeothermal gradient of the sedimentary basin by using a single type of vitrinite reflectance and apatite fission track is insensitive to time and temperature and the obtained simulation output result is inaccurate.

[0046] Among them, the goodness of fit between the low-temperature chronology simulation age and the measured age is calculated through the time-temperature curve based on paleogeothermal gradient (the simulation and measured results include low-temperature chronology age and fission track confined track length distribution), and screening is also carried out through the set goodness-of-fit threshold value. Finally, the paleogeothermal gradient-time curves greater than the goodness-of-fit threshold value are retained, and their mean value is calculated as the final paleogeothermal gradient-time curve, thereby improving the accuracy of restoring the paleogeothermal gradient. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can be obtained based on these drawings.

[0048] Figure 1 It is a flowchart of the steps of the low-temperature age data processing method based on paleogeothermal gradient according to the embodiment of the present invention;

[0049] Figure 2 It is a framework diagram of the low-temperature age data processing system based on paleogeothermal gradient according to the embodiment of the present invention. Detailed implementation manners

[0050] To make the objectives, features, and advantages of the present invention more obvious and understandable, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the embodiments described below are only some embodiments of the present invention, rather than all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts belong to the scope of protection of the present invention.

[0051] The embodiments of the present application provide a low-temperature age data processing method, system, and device based on paleogeothermal gradient, which are used to solve the technical problems that the existing restoration of paleogeothermal gradient in sedimentary basins using a single type of vitrinite reflectance and apatite fission track is insensitive to time and temperature, and the obtained simulation output results are inaccurate.

[0052] Embodiment 1:

[0053] Figure 1 It is a flowchart of the steps of the low-temperature age data processing method based on paleogeothermal gradient according to the embodiment of the present invention.

[0054] As Figure 1 shown, the embodiments of the present invention provide a low-temperature age data processing method based on paleogeothermal gradient, including the following steps:

[0055] S1. Based on the strata of the sedimentary basin above the angular unconformity, obtain M samples of the strata of the sedimentary basin above the angular unconformity and the parameters of each sample. Wherein, N is a natural number greater than 500, and M is a natural number greater than 1.

[0056] It should be noted that, first, within the time range of formation simulation, the Monte Carlo random search method is used to assume N paleogeothermal gradient-time curves for the formation, and the value of N is not less than 10,000. Second, for the strata of sedimentary basins above the angular unconformity, M samples are obtained by sampling at different depths and by drilling, and corresponding parameters are obtained for each sampled sample. Among them, the parameters of the samples include depth, porosity, time, etc. Among them, the Monte Carlo random search method assumes that the curve belongs to a commonly used geological testing technique, which will not be elaborated here.

[0057] In the embodiment of the present invention, the curves, samples, and parameters of the samples obtained in step S1 mainly provide data for subsequent steps.

[0058] S2. Based on the parameters of each sample, use the depth-time equation to solve the burial history time corresponding to the sample.

[0059] It should be noted that mainly the parameters of each sample obtained in step S1 are processed. Specifically, the parameters of each sample are processed to obtain the burial history time of each sample.

[0060] In the embodiment of the present invention, the depth-time equation is:

[0061]

[0062] In the formula, h d is the depth at the lower part of the formation where the sample is located, t d is the time corresponding to the depth at the lower part of the formation where the sample is located, h u is the depth at the upper part of the formation where the sample is located, t u is the time corresponding to the depth at the upper part of the formation where the sample is located, h s is the depth of the sample, t s is the burial history time corresponding to the depth of the sample, y is the porosity of the sample, is the porosity function.

[0063] It should be noted that in the process of using the depth-time equation to solve the burial history time of the corresponding sample for the parameters of each sample, the depth at the lower part of the formation where the sample is located and the depth at the upper part of the formation where the sample is located are known values. The time point is interpolated and solved based on the depth of the sample, and the assumption for its solution is that the sedimentation of the formation where the sample is located is uniform sedimentation.

[0064] S3. Take the burial history time of each sample as the formation boundary, and use the depth calculation formula based on the principle of equal formation mass at different time points to calculate the depth of the corresponding sample, obtain the burial depth of the corresponding sample, and construct the burial history curve of the corresponding sample based on the burial depth and time of the sample.

[0065] It should be noted that mainly for the burial history time of the samples obtained in step S2, the stratigraphic boundary of the strata in the sedimentary basin above the angular unconformity is re-divided. The divided strata are the samples corresponding to the burial history time. The depth calculation formula is used to calculate the depth of each sample, and the burial depth of the corresponding sample is obtained. Then, according to the burial depth and time, the burial history curve of the sample is constructed, that is, a broken line with time as the abscissa and depth as the ordinate is used as the burial history curve.

[0066] In the embodiment of the present invention, the depth calculation formula is:

[0067]

[0068] In the formula, y is the porosity of the sample, is the porosity function, h is the depth of the sample, j is the serial number of the stratigraphic boundary, i is the serial number of the time point corresponding to the stratigraphic boundary, and both i and j are natural numbers greater than 0.

[0069] It should be noted that the principle of equal formation mass at different time points means that the mass of the formation does not change with time, that is, the formation integrals of each time period of the same formation are equal, which can be represented by the following system of equations:

[0070]

[0071]

[0072] .......

[0074] By numerically solving the above system of equations, all unknowns are solved That is, the burial depth of the formation at each time point And so on. In the process of establishing and solving this system of equations, the time nodes only play the role of nodes, and the specific values of their time do not play a role. That is to say, the time values can be assigned after the solution and used to make the burial history curve.

[0075] S4. Assume N paleogeothermal gradient-time curves and construct a geothermal gradient-time function through the N paleogeothermal gradient-time curves, and convert the burial history curve of the corresponding sample into N time-temperature curves of the sample according to the burial depth of the sample and the geothermal gradient-time function.

[0076] It should be noted that mainly according to the geothermal gradient-time function (temperature = depth * geothermal gradient + surface temperature), the burial history curve is converted into the time-temperature curve of the thermal history.

[0077] In an embodiment of the present invention, during the process of processing the paleo-geothermal gradient-time curve to construct the geothermal gradient-time function, a linear interpolation function is used to process N paleo-geothermal gradient-time curves to obtain the geothermal gradient-time function.

[0078] In an embodiment of the present invention, it is assumed that the N paleo-geothermal gradient-time curves are generally curves, mainly including three types of curves: monotonically decreasing, monotonically increasing, and having undulations in the middle. The curve with undulations in the middle is caused by rapid heating events or abnormal cooling events. For example, within a certain range of geothermal gradient at a certain time, the inflection points of the geothermal gradient curve rising and falling are randomly assumed as constraints, and the constraints can be one or more. Therefore, there are undulations in the middle part of the curve. For the curve with undulations in the middle, the Monte Carlo random search method is used to assume that the paleo-geothermal gradient-time curve generally determines the random search range by the number of constraints to enclose its inflection points.

[0079] It should be noted that G is the geothermal gradient and t is the time. Points are plotted in the G-t coordinate system, and then the polyline connecting the nodes is used as the assumed paleo-geothermal gradient-time curve.

[0080] A method for processing low-temperature age data based on paleo-geothermal gradient provided by the present invention includes: obtaining M samples located in the strata of the sedimentary basin above the angular unconformity and the parameters of each sample based on the strata of the sedimentary basin above the angular unconformity; solving the burial history time of the corresponding sample by using the depth-time equation based on the parameters of each sample; taking the burial history time of each sample as the stratigraphic boundary, calculating the burial depth of the corresponding sample by using the depth calculation formula of the principle of equal formation mass at different time points, and constructing the burial history curve of the corresponding sample based on the burial depth and time of the sample; assuming N paleo-geothermal gradient-time curves and constructing the geothermal gradient-time function through the N paleo-geothermal gradient-time curves, and converting the burial history curve of the corresponding sample into the time-temperature curve of the sample according to the burial depth of the sample and the geothermal gradient-time function. The time-temperature curve obtained by the method for processing low-temperature age data based on paleo-geothermal gradient avoids the phenomenon of being insensitive to time and temperature. Using this time-temperature curve as the data for the thermal history simulation of sample low-temperature chronology improves the accuracy of the thermal history simulation results, and solves the technical problems that the existing restoration of the paleo-geothermal gradient of sedimentary basins using a single type of vitrinite reflectance and apatite fission track is insensitive to time and temperature, and the obtained simulation output results are inaccurate.

[0081] It should be noted that the method for processing low-temperature age data based on paleo-geothermal gradient obtains M samples by drilling and sampling from the strata of the sedimentary basin above the angular unconformity, which reduces the difficulty of sample sampling and improves work efficiency.

[0082] In one embodiment of the present invention, the low-temperature chronology data processing method based on paleogeothermal gradient includes:

[0083] Obtain M×N time-temperature curves of the sedimentary basin strata from M samples, extract the time-temperature curves of the same geothermal gradient-time function from the M×N time-temperature curves as a set of thermal history simulation curves, and then obtain N sets of thermal history simulation curves;

[0084] Use the low-temperature chronology method to simulate the ages of the N sets of thermal history simulation curves and calculate the goodness-of-fit values of each set of thermal history simulation curves to obtain N goodness-of-fit values;

[0085] Compare each goodness-of-fit value with the goodness-of-fit threshold, screen out K sets of thermal history simulation curves corresponding to the qualified goodness-of-fit values from the N goodness-of-fit values, and obtain the corresponding K paleogeothermal gradient-time curves according to the K sets of thermal history simulation curves;

[0086] At the same time point, calculate the mean value of each data on the K paleogeothermal gradient-time curves to obtain the temperature mean value of each data; connect the temperature mean values of all data in chronological order into a curve, and this curve is used as the thermal history simulation result of the sedimentary basin strata above the angular unconformity;

[0087] Among them, a set of thermal history simulation curves corresponds to a paleogeothermal gradient-time curve.

[0088] Among them, calculating the mean value of each data on the K paleogeothermal gradient-time curves includes:

[0089] Divide each paleogeothermal gradient-time curve into A equal parts, calculate the mean value of the temperatures at the same time points in all the paleogeothermal gradient-time curves after dividing them into A equal parts to obtain A temperature mean values;

[0090] Connect the A temperature mean values in chronological order into a curve, and this curve is used as the thermal history simulation result of the sedimentary basin strata above the angular unconformity;

[0091] Among them, K is a natural number greater than 5, and A is a natural number greater than 50.

[0092] It should be noted that N groups of thermal history simulation curves of M samples are obtained according to steps S1 to S4. The goodness-of-fit values corresponding to each group of thermal history simulation curves are calculated respectively by existing low-temperature chronology methods. Then, each group of thermal history simulation curves corresponding to the goodness-of-fit greater than the goodness-of-fit threshold and the corresponding paleogeothermal gradient-time curves are selected. Each paleogeothermal gradient-time curve is divided into A equal parts, and the average value of the temperatures at the same time points in all paleogeothermal gradient-time curves after dividing into A equal parts is obtained, resulting in A temperature average values. The A temperature average values are connected into a curve in chronological order, and this curve serves as the thermal history simulation result of the strata of the sedimentary basin above the angular unconformity, improving the accuracy of the thermal history simulation result of the strata of this sedimentary basin. In this embodiment, calculating the goodness-of-fit value by existing low-temperature chronology methods is a relatively mature technology in the field, such as the two published patent documents with publication numbers CN112685909A and CN113722917A. Among them, the goodness-of-fit threshold can be set according to requirements, and the goodness-of-fit threshold is preferably 0.05. A is preferably 100.

[0093] In the embodiment of the present invention, the thermal history simulation result of the strata of the sedimentary basin above the angular unconformity obtained by this low-temperature chronology data processing method based on the paleogeothermal gradient is compared with the test result, more intuitively showing the fitting effect. This low-temperature chronology data processing method based on the paleogeothermal gradient calculates the simulation results (i.e., the goodness-of-fit values) of M samples by screening out the paleogeothermal gradient-time curves corresponding to N groups of thermal history simulation curves; according to the calculated simulation results, N goodness-of-fit values corresponding to N groups of thermal history simulation curves are calculated, a goodness-of-fit threshold is set, the paleogeothermal gradient-time curves greater than the goodness-of-fit threshold are retained, and the average value of the paleogeothermal gradient-time curves is calculated and output as the final paleogeothermal gradient-time curve (or curve). This method is not restricted by low-temperature chronology testing methods, improving the accuracy of paleogeothermal gradient restoration. That is, the goodness-of-fit of the low-temperature chronology simulation age and the measured age is calculated through this time-temperature curve based on the paleogeothermal gradient (the simulation and measured results include low-temperature chronology age and fission track confined track length distribution), and it is also screened by the set goodness-of-fit threshold. Finally, the paleogeothermal gradient-time curves greater than the goodness-of-fit threshold are retained, and their average value is calculated as the final paleogeothermal gradient-time curve, thereby improving the accuracy of paleogeothermal gradient restoration.

[0094] Embodiment Two:

[0095] Figure 2 It is the framework diagram of the low-temperature chronology data processing system based on the paleogeothermal gradient according to the embodiment of the present invention.

[0096] As Figure 2As shown in the figure, an embodiment of the present invention further provides a low-temperature age data processing system based on paleogeothermal gradient, including a data acquisition module 10, a first processing module 20, a second processing module 30, and a third processing module 40;

[0097] The data acquisition module 10 is configured to obtain M samples located in the strata of the sedimentary basin above the angular unconformity and the parameters of each sample based on the strata of the sedimentary basin above the angular unconformity;

[0098] The first processing module 20 is configured to solve the burial history time corresponding to the sample by using the depth-time equation based on the parameters of each sample;

[0099] The second processing module 30 is configured to use the burial history time of each sample as the stratigraphic boundary, calculate the burial depth of the corresponding sample by using the depth calculation formula based on the principle of equal formation mass at different time points, and construct the burial history curve of the corresponding sample based on the burial depth and time of the sample;

[0100] The third processing module 40 is configured to assume N paleogeothermal gradient-time curves and construct a geothermal gradient-time function through the N paleogeothermal gradient-time curves, and convert the burial history curve of the corresponding sample into N time-temperature curves of the sample according to the burial depth of the sample and the geothermal gradient-time function;

[0101] The depth-time equation is:

[0102]

[0103] In the formula, h d is the depth of the lower part of the stratum where the sample is located, t d is the time corresponding to the depth of the lower part of the stratum where the sample is located, h u is the depth of the upper part of the stratum where the sample is located, t u is the time corresponding to the depth of the upper part of the stratum where the sample is located, h s is the depth of the sample, t s is the burial history time corresponding to the depth of the sample, y is the porosity of the sample, is the porosity function;

[0104] The depth calculation formula is:

[0105]

[0106] In the formula, y is the porosity of the sample, is the porosity function, h is the depth of the sample, j is the serial number of the stratigraphic boundary, i is the serial number of the time point corresponding to the stratigraphic boundary, and both i and j are natural numbers greater than 0;

[0107] Wherein, N is a natural number greater than 500, and M is a natural number greater than 1.

[0108] In an embodiment of the present invention, the low-temperature chronology data processing system based on paleogeothermal gradient includes an output module 50. The output module 50 is configured to obtain M×N time-temperature curves of the sedimentary basin strata according to M samples, extract the time-temperature curves of the same geothermal gradient-time function from the M×N time-temperature curves as a set of thermal history simulation curves, and further obtain N sets of thermal history simulation curves; use the low-temperature chronology method to simulate the ages of the N sets of thermal history simulation curves to calculate the goodness-of-fit values of each set of thermal history simulation curves, and obtain N goodness-of-fit values; compare each goodness-of-fit value with the goodness-of-fit threshold, and screen out K sets of thermal history simulation curves corresponding to the qualified goodness-of-fit values from the N goodness-of-fit values, and obtain the corresponding K paleogeothermal gradient-time curves according to the K sets of thermal history simulation curves; at the same time point, calculate the average value of each data on the K paleogeothermal gradient-time curves to obtain the temperature average value of each data; connect the temperature average values of all data in chronological order into a curve, and this curve is used as the thermal history simulation result of the sedimentary basin strata above the angular unconformity; wherein, a set of thermal history simulation curves corresponds to one paleogeothermal gradient-time curve.

[0109] In an embodiment of the present invention, the output module 50 is further configured to divide each paleogeothermal gradient-time curve into A equal parts, calculate the average value of the temperatures at the same time point in all the paleogeothermal gradient-time curves after dividing them into A equal parts to obtain A temperature average values, and connect the A temperature average values in chronological order into a curve, and this curve is used as the thermal history simulation result of the sedimentary basin strata above the angular unconformity; wherein, A is a natural number greater than 50.

[0110] It should be noted that the modules in the system of the second embodiment correspond to the steps in the method of the first embodiment. The step content of the method of the first embodiment has been elaborated in detail in the first embodiment, and the content of the modules will not be elaborated in detail in the second embodiment here.

[0111] Embodiment Three:

[0112] An embodiment of the present invention provides a low-temperature chronology data processing device based on paleogeothermal gradient, including a processor and a memory;

[0113] The memory is used to store program codes and transmit the program codes to the processor;

[0114] The processor is configured to execute the above-mentioned low-temperature chronology data processing method based on paleogeothermal gradient according to the instructions in the program codes.

[0115] It should be noted that the processor is used to execute the steps in the above embodiments of the method for processing low-temperature age data based on paleogeothermal gradient according to the instructions in the program code. Alternatively, when the processor executes the computer program, it realizes the functions of each module / unit in the above embodiments of each system / device.

[0116] Exemplarily, the computer program can be divided into one or more modules / units. One or more modules / units are stored in the memory and executed by the processor to complete the present application. One or more modules / units can be a series of computer program instruction segments capable of performing specific functions, and the instruction segments are used to describe the execution process of the computer program in the terminal device.

[0117] The terminal device can be a computing device such as a desktop computer, a notebook, a palm computer, and a cloud server. The terminal device may include, but is not limited to, a processor and a memory. Those skilled in the art can understand that this does not constitute a limitation on the terminal device, and it may include more or fewer components than shown in the figure, or combine certain components, or different components. For example, the terminal device may further include an input / output device, a network access device, a bus, etc.

[0118] The so-called processor may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0119] The memory may be an internal storage unit of the terminal device, such as the hard disk or memory of the terminal device. The memory may also be an external storage device of the terminal device, such as a plug-in hard disk equipped on the terminal device, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. Further, the memory may also include both the internal storage unit and the external storage device of the terminal device. The memory is used to store the computer program and other programs and data required by the terminal device. The memory may also be used to temporarily store the data that has been output or will be output.

[0120] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the systems, devices, and units described above can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated herein.

[0121] In several embodiments provided in the present application, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is only a logical function division, and there may be other division methods in actual implementation. For example, 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 displayed or discussed couplings or direct couplings or communication connections to each other can be indirect couplings or communication connections through some interfaces, devices, or units, and can be in electrical, mechanical, or other forms.

[0122] The units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they can be located in one place, or distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0123] In addition, in each embodiment of the present invention, the functional units can be integrated into one processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated units can be implemented in the form of hardware or in the form of software functional units.

[0124] If the above-mentioned integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in each embodiment of the present invention. The foregoing storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0125] As described above, the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for processing low-temperature age data based on paleogeothermal gradient, characterized in that It includes the following steps: Based on the strata of the sedimentary basin above the angular unconformity, M samples of the strata of the sedimentary basin above the angular unconformity and the parameters of each sample are obtained; Based on the parameters of each sample, the burial history time corresponding to the sample is solved by using the depth-time equation; Taking the burial history time of each sample as the stratigraphic boundary, the depth of the corresponding sample is calculated by using the depth calculation formula based on the principle of equal formation mass at different time points, and the burial history curve of the corresponding sample is constructed based on the burial depth and time of the sample; Assume N paleogeothermal gradient-time curves and construct a geothermal gradient-time function through the N paleogeothermal gradient-time curves, and convert the burial history curve of the corresponding sample into N time-temperature curves of the sample according to the burial depth of the sample and the geothermal gradient-time function; Wherein, N is a natural number greater than 500, and M is a natural number greater than 1; The depth-time equation is: where hd is the depth of the lower part of the formation where the sample is located, td is the time corresponding to the depth of the lower part of the formation where the sample is located, hu is the depth of the upper part of the formation where the sample is located, tu is the time corresponding to the depth of the upper part of the formation where the sample is located, hs is the depth of the sample, ts is the burial history time corresponding to the depth of the sample, y is the porosity of the sample, is the porosity function; The depth calculation formula is: where y is the porosity of the sample, is the porosity function, h is the depth of the sample, j is the serial number of the formation boundary, i is the serial number of the time point corresponding to the formation boundary, and both i and j are natural numbers greater than 0; Wherein, N is a natural number greater than 500, and M is a natural number greater than 1.

2. The low-temperature age data processing method based on paleogeothermal gradient according to claim 1, wherein It includes: According to the M samples, M×N time-temperature curves of the sedimentary basin strata are obtained. The time-temperature curves of the same geothermal gradient-time function are extracted from the M×N time-temperature curves as a set of thermal history simulation curves, and then N sets of thermal history simulation curves are obtained; Using the low-temperature chronology method to simulate the age of the N sets of thermal history simulation curves, and calculating the goodness-of-fit value of each set of thermal history simulation curves, to obtain N goodness-of-fit values; Comparing each goodness-of-fit value with the goodness-of-fit threshold, screening out K sets of thermal history simulation curves corresponding to the qualified goodness-of-fit values from the N goodness-of-fit values, and obtaining the corresponding K paleogeothermal gradient-time curves according to the K sets of thermal history simulation curves; At the same time point, the mean value of each data on the K paleogeothermal gradient-time curves is obtained to get the temperature mean value of each data; Connecting the temperature mean values of all the data in chronological order to form a curve, and this curve is used as the thermal history simulation result of the strata of the sedimentary basin above the angular unconformity; Wherein, a set of thermal history simulation curves corresponds to a paleogeothermal gradient-time curve.

3. The method for processing low-temperature chronology data based on paleogeothermal gradient according to claim 2, wherein, Obtaining the mean value of each data on the K paleogeothermal gradient-time curves includes: Dividing each paleogeothermal gradient-time curve into A equal parts, and obtaining the mean value of the temperatures at the same time point in all the paleogeothermal gradient-time curves after dividing them into A equal parts, to obtain A temperature mean values; Connecting the A temperature mean values in chronological order to form a curve, and this curve is used as the thermal history simulation result of the strata of the sedimentary basin above the angular unconformity; Wherein, A is a natural number greater than 50.

4. The low-temperature age data processing method based on paleogeothermal gradient according to claim 1, characterized in that It includes: Using a first-order linear interpolation function to construct a geothermal gradient-time function for the N paleogeothermal gradient-time curves.

5. A low-temperature chronology data processing system based on paleogeothermal gradient, characterized in that, It includes: A data acquisition module, a first processing module, a second processing module and a third processing module; The data acquisition module is used to obtain M samples of the strata of the sedimentary basin above the angular unconformity and the parameters of each sample based on the strata of the sedimentary basin above the angular unconformity; The first processing module is configured to solve the burial history time corresponding to each sample by using a depth-time equation based on the parameters of each sample; The second processing module is configured to use the burial history time of each sample as a stratigraphic boundary, and use a depth calculation formula based on the principle of equal formation mass at different time points to calculate the burial depth of the corresponding sample, and construct a burial history curve corresponding to the sample based on the burial depth and time of the sample; The third processing module is configured to assume N paleogeothermal gradient-time curves and construct a geothermal gradient-time function through the N paleogeothermal gradient-time curves, and convert the burial history curve corresponding to the sample into N time-temperature curves of the sample according to the burial depth of the sample and the geothermal gradient-time function; The depth-time equation is: where hd is the depth of the lower part of the formation where the sample is located, td is the time corresponding to the depth of the lower part of the formation where the sample is located, hu is the depth of the upper part of the formation where the sample is located, tu is the time corresponding to the depth of the upper part of the formation where the sample is located, hs is the depth of the sample, ts is the burial history time corresponding to the depth of the sample, and y is the porosity of the sample, is the porosity function; The depth calculation formula is: where y is the porosity of the sample, is the porosity function, h is the depth of the sample, j is the serial number of the formation boundary, i is the serial number of the time point corresponding to the formation boundary, and both i and j are natural numbers greater than 0; Where N is a natural number greater than 500, and M is a natural number greater than 1.

6. The low-temperature age data processing system based on the paleogeothermal gradient according to claim 5, wherein, It includes an output module. The output module is configured to obtain M×N time-temperature curves of the sedimentary basin strata according to M samples, extract the time-temperature curves of the same geothermal gradient-time function from the M×N time-temperature curves as a set of thermal history simulation curves, and then obtain N sets of thermal history simulation curves. The low-temperature chronology method is used to simulate the age of each set of the thermal history simulation curves to calculate the goodness-of-fit value of each set of the thermal history simulation curves, and N goodness-of-fit values are obtained; Compare each of the goodness-of-fit values with a goodness-of-fit threshold, and screen out K sets of thermal history simulation curves corresponding to the qualified goodness-of-fit values from the N goodness-of-fit values, and obtain the corresponding K paleogeothermal gradient-time curves according to the K sets of thermal history simulation curves; at the same time point, calculate the mean value of each data on the K paleogeothermal gradient-time curves to obtain the temperature mean value of each data; connect the temperature mean values of all the data in chronological order to form a curve, and this curve is used as the thermal history simulation result of the sedimentary basin strata above the angular unconformity; among them, a set of thermal history simulation curves corresponds to one paleogeothermal gradient-time curve.

7. The low-temperature age data processing system based on the paleogeothermal gradient according to claim 6, characterized in that The output module is further configured to divide each of the paleogeothermal gradient-time curves into A equal parts, calculate the mean value of the temperatures at the same time point in all the paleogeothermal gradient-time curves after dividing them into A equal parts to obtain A temperature mean values, and connect the A temperature mean values in chronological order to form a curve, and this curve is used as the thermal history simulation result of the sedimentary basin strata above the angular unconformity; where A is a natural number greater than 50.

8. A low-temperature age data processing device based on paleogeothermal gradient, characterized in that, It includes a processor and a memory; The memory is configured to store program codes and transmit the program codes to the processor; The processor is configured to execute the low-temperature chronology data processing method based on paleogeothermal gradient as described in any one of claims 1-4 according to the instructions in the program codes.

Images