A method, device and related equipment for determining zircon fission track visual age

Abstract

The present application discloses a method, device and related equipment for determining the apparent age of zircon fission tracks, the method comprising: determining an annealing model of zircon fission tracks, the annealing model recording the constraint relationship between track length ratio, time and temperature; determining the track length ratio of the zircon sample in each time interval within a target simulation period based on the annealing model and a preset thermal history curve, the target simulation period being the period included in the thermal history curve; determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period; determining the apparent age of the fission track of the zircon sample based on the track density of the zircon sample in each time interval within the target simulation period. The present application calculates the apparent age of the fission track of the zircon sample by forward simulation of the fission track, which has certain practical significance in the application of low-temperature chronology.

Description

Technical Field

[0001] The present application relates to the technical field of cryogenic chronology, and more specifically, to a method and device for determining the apparent age of zircon fission tracks and related equipment. Background Art

[0002] The zircon fission track method is a low-temperature chronology method commonly used in the geological community. This method can not only provide the apparent age corresponding to the closure temperature, but also restore the thermal history of the research object. When using the zircon fission track method for thermal history simulation, the first thing to do is to establish a forward model, and then establish an inversion process based on this.

[0003] Among them, the inversion process can apply various existing low-temperature chronology thermal history inversion methods in our low-temperature thermal history simulation, such as the Monte Carlo method, etc. Since the parameters of the forward model and the annealing model more reflect the characteristics of low-temperature chronology itself, the forward model cannot simply apply various low-temperature chronology thermal history forward methods. How to establish a forward model and determine the apparent age of zircon fission tracks based on it has become a topic worth exploring. Summary of the invention

[0004] In view of this, the present application provides a method, apparatus and related equipment for determining the apparent age of zircon fission tracks, so as to determine the apparent age of zircon samples through a forward model.

[0005] To achieve the above-mentioned purpose, the first aspect of the present application provides a method for determining the apparent age of zircon fission tracks, comprising:

[0006] Determining an annealing model for zircon fission tracks, the annealing model recording the constraint relationship between track length ratio, time and temperature;

[0007] Based on the annealing model and a preset thermal history curve, determining the track length ratio of the zircon sample in each time interval within a target simulation period, wherein the target simulation period is a period included in the thermal history curve;

[0008] Based on the track length ratio of the zircon sample in each time interval within the target simulation period, the track density of the zircon sample in each time interval within the target simulation period is determined;

[0009] The fission track apparent age of the zircon sample is determined based on the track density of the zircon sample in each time interval within the target simulation period.

[0010] Preferably, the time intervals are equal-length, non-overlapping, and sequentially adjacent time periods within the target simulation period; and based on the annealing model and a preset thermal history curve, the process of determining the track length ratio of the zircon sample in each time interval within the target simulation period comprises:

[0011] Based on the thermal history curve, determining the equivalent temperature of the zircon sample in each time interval;

[0012] Based on the annealing model and the equivalent temperature of each time interval, the track length ratio of the zircon sample in each time interval is determined.

[0013] Preferably, the process of determining the equivalent temperature of the zircon sample in each time interval based on the thermal history curve comprises:

[0014] Based on the thermal history curve, determining the temperature of the first endpoint and the temperature of the second endpoint of the zircon sample at each time interval;

[0015] An average of the temperature at the first end point and the temperature at the second end point of the time interval is determined as the equivalent temperature of the time interval.

[0016] Preferably, the process of determining the track length ratio of the zircon sample in each time interval based on the annealing model and the equivalent temperature of each time interval includes:

[0017] Based on the annealing model and the equivalent temperature of the first time interval, determining the track length ratio of the zircon sample in the first time interval, wherein the first time interval is the most recent time interval;

[0018] Determining an equivalent constant time of the previous time interval relative to the next time interval based on the annealing model, the track length ratio of the previous time interval, and the equivalent temperature of the next time interval;

[0019] Determining the track length ratio of the zircon sample in the latter time interval based on the annealing model, the equivalent constant time and the equivalent temperature in the latter time interval;

[0020] The previous time interval and the next time interval are adjacent time intervals, and the time of the previous time interval is closer to the current time.

[0021] Preferably, the process of determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period includes:

[0022] The following equation is used to convert the track length ratio of the zircon sample in each time interval within the target simulation period into the track density:

[0023] ρ i =1.25×(Fr i -0.2)

[0024] Among them, Fr i The track length ratio of the zircon sample in the i-th time interval within the target simulation period, ρi Corresponding to Fr i The track density.

[0025] Preferably, the process of determining the fission track apparent age of the zircon sample based on the track density of the zircon sample at each time interval within the target simulation period includes:

[0026] The fission track apparent age of zircon samples is calculated using the following equation:

[0027] Age=k·Δt∑ i ρ i

[0028] Among them, ρ i is the ith track density, k is a preset correction factor, and Δt is the length of the time interval.

[0029] Preferably, the annealing model of the zircon fission track is:

[0030]

[0031] Wherein, c0, c1, c2, c3, α and β are preset constants, tx is time, Ty is temperature, r clr is the track length ratio.

[0032] The second aspect of the present application provides a device for determining the apparent age of zircon fission tracks, comprising:

[0033] A model determination unit, used to determine an annealing model of zircon fission tracks, wherein the annealing model records the constraint relationship between track length ratio, time and temperature;

[0034] A ratio determination unit, for determining the track length ratio of the zircon sample in each time interval within a target simulation period based on the annealing model and a preset thermal history curve, wherein the target simulation period is a period where the thermal history curve is located;

[0035] A density determination unit, for determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period;

[0036] The age determination unit is used to determine the fission track apparent age of the zircon sample based on the track density of the zircon sample in each time interval within the target simulation period.

[0037] A third aspect of the present application provides a device for determining the apparent age of zircon fission tracks, comprising: a memory and a processor;

[0038] The memory is used to store programs;

[0039] The processor is used to execute the program to implement the various steps of the above-mentioned method for determining the apparent age of zircon fission tracks.

[0040] A fourth aspect of the present application provides a storage medium having a computer program stored thereon, which, when executed by a processor, implements the various steps of the above-mentioned method for determining the apparent age of zircon fission tracks.

[0041] It can be known from the above technical solution that the present application first determines the annealing model of the zircon fission track, wherein the annealing model records the constraint relationship between the track length ratio, time and temperature. Then, based on the annealing model and the preset thermal history curve, the track length ratio of the zircon sample in each time interval within the target simulation period is determined. Wherein, the target simulation period is the period included in the thermal history curve. Next, based on the track length ratio of the zircon sample in each time interval within the target simulation period, the track density of the zircon sample in each time interval within the target simulation period is determined. Finally, based on the track density of the zircon sample in each time interval within the target simulation period, the fission track apparent age of the zircon sample is determined. The present application calculates the fission track apparent age of the zircon sample by means of forward simulation of the fission track, which has certain practical significance in the application of low-temperature chronology. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0043] Figure 1 A schematic diagram of a method for determining the apparent age of zircon fission tracks disclosed in an embodiment of the present application;

[0044] Figure 2 A schematic diagram of equally divided temperature history disclosed in the embodiments of the present application;

[0045] Figure 3 A schematic diagram of a device for determining the visual age of zircon fission tracks disclosed in an embodiment of the present application;

[0046] Figure 4 This is a schematic diagram of the zircon fission track visual age determination device disclosed in an embodiment of the present application. DETAILED DESCRIPTION

[0047] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0048] The following is a method for determining the apparent age of zircon fission tracks provided in the embodiments of the present application. Figure 1 The method for determining the apparent age of zircon fission tracks provided in the embodiment of the present application may include the following steps:

[0049] Step S101, determining the annealing model of zircon fission track.

[0050] The annealing model records the constraint relationship between track length ratio, time and temperature, and the track length ratio refers to the ratio of the track length to the original length. According to the annealing model, the track length ratio can be calculated by temperature and time.

[0051] Step S102, based on the annealing model and the preset thermal history curve, determine the track length ratio of the zircon sample in each time interval within the target simulation period.

[0052] The target simulation period is the period included in the thermal history curve, which can be the entire period covered by the thermal history curve or a certain period therein. Figure 2 The thermal history curve shown is in the period from 0 to 125 Ma, so the target simulation period can be set to 0 to 125 Ma.

[0053] Since the thermal history curve contains the temperature value of the zircon sample at each time node, on this basis, combined with the annealing model, the track length ratio of the zircon sample at each time node can be calculated. For the thermal history curve, the time points are uncountable. For the convenience of calculation, multiple time intervals can be divided for the target simulation period, and the track length ratio of the zircon sample in each time interval can be calculated.

[0054] Step S103, determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period.

[0055] For example, existing formulas can be used to convert track length ratios into track density.

[0056] Step S104, determining the fission track apparent age of the zircon sample based on the track density of the zircon sample at each time interval within the target simulation period.

[0057] Fission tracks are radiation damage caused to minerals by fission fragments with positive charge and energy within their path range. They are very small and can only be seen with a high-power electron microscope. The number of spontaneous fission tracks of uranium-238 is proportional to the uranium content and formation age of the mineral. The uranium content can be determined by neutron irradiation to produce induced fission tracks of uranium-235 in the sample. At this time, when the etching conditions of spontaneous fission tracks and induced fission tracks are exactly the same and the irradiation conditions (thermal neutron dose) are selected, the age can be calculated based on the ratio of the spontaneous fission track density to the induced fission track density of the mineral particles after etching. In some embodiments of the present application, a preset correction factor can be used to correct the track density of the zircon sample in each time interval within the target simulation period, and finally the apparent fission track age of the zircon sample is calculated.

[0058] The present application first determines the annealing model of zircon fission tracks, wherein the annealing model records the constraint relationship between track length ratio, time and temperature. Then, based on the annealing model and the preset thermal history curve, the track length ratio of the zircon sample in each time interval within the target simulation period is determined. Wherein, the target simulation period is the period included in the thermal history curve. Next, based on the track length ratio of the zircon sample in each time interval within the target simulation period, the track density of the zircon sample in each time interval within the target simulation period is determined. Finally, based on the track density of the zircon sample in each time interval within the target simulation period, the fission track apparent age of the zircon sample is determined. The present application calculates the fission track apparent age of the zircon sample by means of fission track forward simulation, which has certain practical significance in the application of low-temperature chronology.

[0059] In some embodiments of the present application, the annealing model of the zircon fission track determined in step S101 above may be:

[0060]

[0061] The above equation can be recorded as the length ratio formula of the annealing model, where c0, c1, c2, c3, α and β are preset constants. For example, c0 can be 6.24534, c1 can be 0.11977, c2 can be 314.937, c3 can be 14.2868, and α can be (60×60×24×365×10 6 ), β can be taken as 0.05721; tx is time, Ty is temperature, r clr is the track length ratio.

[0062] The above equation of the annealing model can be rewritten as follows to obtain the time formula of the annealing model:

[0063]

[0064] Among them, t clr is the time, rx is the track length ratio, and other parameters are the same as above. Based on this equation, the time t can be calculated according to the track length ratio rx and temperature Ty: clr .

[0065] In some embodiments of the present application, the above-mentioned time intervals are equal-length, non-overlapping, and sequentially adjacent time periods within the target simulation period. For example, the target simulation period can be divided into multiple time intervals using a preset step size.

[0066] Based on this, the process of determining the track length ratio of the zircon sample in each time interval within the target simulation period in step S102 based on the annealing model and the preset thermal history curve may include:

[0067] S1, based on the thermal history curve, determine the equivalent temperature of the zircon sample in each time interval.

[0068] Among them, since the thermal history curve is a two-dimensional curve in the time and temperature dimensions, for a certain time interval, its equivalent temperature depends on the temperature values ​​corresponding to the time interval. The average of the temperature values ​​corresponding to the time interval can be taken. To simplify the calculation, sampling can be performed on the time interval and then the average value can be taken.

[0069] S2, based on the annealing model and the equivalent temperature of each time interval, determining the track length ratio of the zircon sample in each time interval.

[0070] Since the annealing model describes the constraint relationship between track length ratio, time and temperature, the track length ratio of the zircon sample in each time interval can be calculated based on this constraint relationship and the equivalent temperature of each time interval.

[0071] In some embodiments of the present application, the process of determining the equivalent temperature of the zircon sample in each time interval based on the thermal history curve in S1 may include:

[0072] S11, based on the thermal history curve, determining the temperature of the first endpoint and the temperature of the second endpoint of the zircon sample in each time interval.

[0073] S12, determining an average of the temperature of the first endpoint and the temperature of the second endpoint of the time interval as the equivalent temperature of the time interval.

[0074] That is, samples are taken in each time interval to obtain the average value, and the sampling points are located at the first and last endpoints of each time interval. Figure 2The thermal history curve shown in FIG. 1 divides the entire time period of the thermal history curve into 100 time intervals with a step size of Δt, and correspondingly obtains 101 nodes. The time and temperature coordinates of each node are (t1, T1), (t2, T2), ..., (t101, T101), etc. Then, for the time interval (t1, t2), the equivalent temperature is (T1+T2) / 2, for the time interval (t2, t3), the equivalent temperature is (T2+T3) / 2, and so on.

[0075] In some embodiments of the present application, the process of determining the track length ratio of the zircon sample in each time interval based on the annealing model and the equivalent temperature of each time interval in S2 may include:

[0076] S21, determining the track length ratio of the zircon sample in the first time interval based on the annealing model and the equivalent temperature of the first time interval.

[0077] Among them, the first time interval is the time interval closest to the present. Figure 2 In the thermal history curve shown, the first time interval is (t1, t2).

[0078] S22, determining an equivalent constant time of the previous time interval relative to the next time interval based on the annealing model, the track length ratio of the previous time interval, and the equivalent temperature of the next time interval.

[0079] Among them, the previous time interval and the next time interval are adjacent time intervals, and the time of the previous time interval is closer to the current time. The previous time interval and the next time interval appear in pairs. For example, for Figure 2 For the thermal history curve shown, if the previous time interval is (t1, t2), then the next time interval is (t2, t3); if the previous time interval is (t2, t3), then the next time interval is (t3, t4), and so on.

[0080] Taking the first time interval (t1, t2) as an example, based on the time formula of the annealing model provided above, the equivalent constant time of the time interval (t1, t2) relative to the next time interval (t2, t3) can be expressed as:

[0081] t1 equivalent =t clr (Fr1,(T2+T3) / 2)

[0082] Where Fr1 is the track length ratio of the time interval (t1, t2), and (T2+T3) / 2 is the equivalent temperature of the time interval (t2, t3).

[0083] And so on, each equivalent constant time can be calculated.

[0084] S23, determining the track length ratio of the zircon sample in the latter time interval based on the annealing model, the equivalent constant time and the equivalent temperature in the latter time interval.

[0085] Taking the time interval (t2, t3) as an example, the track length ratio can be expressed as:

[0086] Fr2=r clr (Δt+t1 equivalent ,(T2+T3) / 2)

[0087] The following is a detailed description of the process of calculating the track length ratio of zircon samples in each time interval using the above method:

[0088] In the first round, since S21 has calculated the track length ratio Fr1 of the zircon sample in the first time interval (t1, t2), the first time interval (t1, t2) is taken as the previous time interval, and the equivalent constant time t1 of the previous time interval (t1, t2) relative to the next time interval (t2, t3) can be calculated according to S22. equivalent ; Then, according to S23, the track length ratio Fr2 of the next time interval (t2, t3) can be calculated.

[0089] In the second round, the time interval (t2, t3) is taken as the previous time interval. On the basis of calculating the track length ratio Fr2 of the time interval (t2, t3), the equivalent constant time t2 of the previous time interval (t2, t3) relative to the next time interval (t3, t4) can be calculated according to S22. equivalent ; Then, according to S23, the track length ratio Fr3 of the next time interval (t3, t4) can be calculated.

[0090] By analogy, the track length ratio of the zircon sample in each time interval can be calculated.

[0091] It should be noted that when the i-th round is reached, the calculated track length ratio Fr i Less than a preset value (typically 0.36 / 1.25+0.2), or the track length ratio Fr i When a complex number appears, the calculation is stopped, and the number of time steps retained is n=i-1. Within this time step, annealing occurs gradually during the temperature increase, or the track is retained during the temperature decrease.

[0092] In some embodiments of the present application, the process of determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period in step S103 may include:

[0093] The following equation is used to convert the track length ratio of the zircon sample in each time interval within the target simulation period into the track density:

[0094] ρ i =1.25×(Fr i -0.2)

[0095] Among them, Fr i The track length ratio of the zircon sample in the i-th time interval within the target simulation period, ρ i Corresponding to Fr i The track density.

[0096] In some embodiments of the present application, the process of determining the fission track apparent age of the zircon sample based on the track density of the zircon sample in each time interval within the target simulation period in step S104 may include:

[0097] The fission track apparent age of zircon samples is calculated using the following equation:

[0098] Age=k·Δt∑ i ρ i

[0099] Among them, ρ i is the ith track density, k is a preset correction factor, whose typical value is 1 / 0.995 or 1, and Δt is the length of the time interval.

[0100] The following is a description of a device for determining the apparent age of a zircon based on fission tracks provided in an embodiment of the present application. The device for determining the apparent age of a zircon based on fission tracks described below and the method for determining the apparent age of a zircon based on fission tracks described above can be referred to in correspondence with each other.

[0101] See also Figure 3 The zircon fission track visual age determination device provided in the embodiment of the present application may include:

[0102] A model determination unit 21, used to determine an annealing model of zircon fission tracks, wherein the annealing model records the constraint relationship between track length ratio, time and temperature;

[0103] A ratio determination unit 22, for determining the track length ratio of the zircon sample in each time interval within a target simulation period based on the annealing model and a preset thermal history curve, wherein the target simulation period is a period where the thermal history curve is located;

[0104] A density determination unit 23, for determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period;

[0105] The age determination unit 24 is used to determine the fission track apparent age of the zircon sample based on the track density of the zircon sample at each time interval within the target simulation period.

[0106] In some embodiments of the present application, the annealing model of the zircon fission track is:

[0107]

[0108] Wherein, c0, c1, c2, c3, α and β are preset constants, tx is time, Ty is temperature, r clr is the track length ratio.

[0109] In some embodiments of the present application, the time intervals are equal-length, non-overlapping, and sequentially adjacent time periods within the target simulation period. The process in which the ratio determination unit 22 determines the track length ratio of each time interval of the zircon sample within the target simulation period based on the annealing model and the preset thermal history curve may include:

[0110] Based on the thermal history curve, determining the equivalent temperature of the zircon sample in each time interval;

[0111] Based on the annealing model and the equivalent temperature of each time interval, the track length ratio of the zircon sample in each time interval is determined.

[0112] In some embodiments of the present application, the process of the ratio determination unit 22 determining the equivalent temperature of the zircon sample in each time interval based on the thermal history curve may include:

[0113] Based on the thermal history curve, determining the temperature of the first endpoint and the temperature of the second endpoint of the zircon sample at each time interval;

[0114] An average of the temperature at the first end point and the temperature at the second end point of the time interval is determined as the equivalent temperature of the time interval.

[0115] In some embodiments of the present application, the process in which the ratio determination unit 22 determines the track length ratio of the zircon sample in each time interval based on the annealing model and the equivalent temperature of each time interval may include:

[0116] Based on the annealing model and the equivalent temperature of the first time interval, determining the track length ratio of the zircon sample in the first time interval, wherein the first time interval is the most recent time interval;

[0117] Determining an equivalent constant time of the previous time interval relative to the next time interval based on the annealing model, the track length ratio of the previous time interval, and the equivalent temperature of the next time interval;

[0118] Determining the track length ratio of the zircon sample in the latter time interval based on the annealing model, the equivalent constant time and the equivalent temperature in the latter time interval;

[0119] The previous time interval and the next time interval are adjacent time intervals, and the time of the previous time interval is closer to the current time.

[0120] In some embodiments of the present application, the density determination unit 23 determines the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period, which may include:

[0121] The following equation is used to convert the track length ratio of the zircon sample in each time interval within the target simulation period into the track density:

[0122] ρ i =1.25×(Fr i -0.2)

[0123] Among them, Fr i The track length ratio of the zircon sample in the i-th time interval within the target simulation period, ρ i Corresponding to Fr i The track density.

[0124] In some embodiments of the present application, the process of determining the fission track apparent age of the zircon sample by the age determination unit 24 based on the track density of the zircon sample in each time interval within the target simulation period may include:

[0125] The fission track apparent age of zircon samples is calculated using the following equation:

[0126] Age=k·Δt∑ i ρ i

[0127] Among them, ρ i is the ith track density, k is a preset correction factor, and Δt is the length of the time interval.

[0128] The zircon fission track apparent age determination device provided in the embodiment of the present application can be applied to a zircon fission track apparent age determination device, such as a computer, etc. Optionally, Figure 4 The hardware structure diagram of the zircon fission track visual age determination device is shown in FIG. Figure 4 The hardware structure of the zircon fission track visual age determination device may include: at least one processor 31, at least one communication interface 32, at least one memory 33 and at least one communication bus 34.

[0129] In the embodiment of the present application, the number of the processor 31, the communication interface 32, the memory 33, and the communication bus 34 is at least one, and the processor 31, the communication interface 32, and the memory 33 communicate with each other through the communication bus 34;

[0130] The processor 31 may be a central processing unit CPU, or an application-specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present application, etc.;

[0131] The memory 33 may include a high-speed RAM memory, and may also include a non-volatile memory, such as at least one disk memory;

[0132] The memory 33 stores a program, and the processor 31 can call the program stored in the memory 33, and the program is used to:

[0133] Determining an annealing model for zircon fission tracks, the annealing model recording the constraint relationship between track length ratio, time and temperature;

[0134] Based on the annealing model and a preset thermal history curve, determining the track length ratio of the zircon sample in each time interval within a target simulation period, wherein the target simulation period is a period included in the thermal history curve;

[0135] Based on the track length ratio of the zircon sample in each time interval within the target simulation period, the track density of the zircon sample in each time interval within the target simulation period is determined;

[0136] The fission track apparent age of the zircon sample is determined based on the track density of the zircon sample in each time interval within the target simulation period.

[0137] Optionally, the detailed functions and extended functions of the program may refer to the above description.

[0138] The embodiment of the present application further provides a storage medium, which may store a program suitable for execution by a processor, wherein the program is used to:

[0139] Determining an annealing model for zircon fission tracks, the annealing model recording the constraint relationship between track length ratio, time and temperature;

[0140] Based on the annealing model and a preset thermal history curve, determining the track length ratio of the zircon sample in each time interval within a target simulation period, wherein the target simulation period is a period included in the thermal history curve;

[0141] Based on the track length ratio of the zircon sample in each time interval within the target simulation period, the track density of the zircon sample in each time interval within the target simulation period is determined;

[0142] The fission track apparent age of the zircon sample is determined based on the track density of the zircon sample in each time interval within the target simulation period.

[0143] Optionally, the detailed functions and extended functions of the program may refer to the above description.

[0144] In summary:

[0145] The present application first determines the annealing model of zircon fission tracks, wherein the annealing model records the constraint relationship between track length ratio, time and temperature. Then, based on the annealing model and the preset thermal history curve, the track length ratio of the zircon sample in each time interval within the target simulation period is determined. Wherein, the target simulation period is the period included in the thermal history curve. Next, based on the track length ratio of the zircon sample in each time interval within the target simulation period, the track density of the zircon sample in each time interval within the target simulation period is determined. Finally, based on the track density of the zircon sample in each time interval within the target simulation period, the fission track apparent age of the zircon sample is determined. The present application calculates the fission track apparent age of the zircon sample by means of fission track forward simulation, which has certain practical significance in the application of low-temperature chronology.

[0146] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0147] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The various embodiments can be combined as needed, and the same or similar parts can refer to each other.

[0148] The above description of the disclosed embodiments enables those skilled in the art to implement or use the present application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application will not be limited to the embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for determining the apparent age of zircon fission tracks, characterized in that: include: Determining an annealing model for zircon fission tracks, the annealing model recording the constraint relationship between track length ratio, time and temperature; Based on the annealing model and a preset thermal history curve, determining the track length ratio of the zircon sample in each time interval within a target simulation period, wherein the target simulation period is a period included in the thermal history curve; Based on the track length ratio of the zircon sample in each time interval within the target simulation period, the track density of the zircon sample in each time interval within the target simulation period is determined; Determine the fission track apparent age of the zircon sample based on the track density of the zircon sample at each time interval within the target simulation period; Wherein, each time interval is a time period of equal length, non-overlapping, and sequentially adjacent within the target simulation period; based on the annealing model and the preset thermal history curve, the process of determining the track length ratio of each time interval of the zircon sample within the target simulation period includes: Based on the thermal history curve, determining the equivalent temperature of the zircon sample in each time interval; Based on the annealing model and the equivalent temperature of each time interval, determining the track length ratio of the zircon sample in each time interval; The process of determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period includes: The following equation is used to convert the track length ratio of the zircon sample in each time interval within the target simulation period into the track density: ； in, The track length ratio of the zircon sample in the i-th time interval within the target simulation period, To correspond to The track density of The process of determining the fission track apparent age of the zircon sample based on the track density of the zircon sample at each time interval within the target simulation period includes: The fission track apparent age of zircon samples was calculated using the following equation: : ； Where k is the preset correction factor, is the length of the time interval; The annealing model of the zircon fission track is: ； in, , , , , and is a preset constant, For time, is the temperature, is the track length ratio.

2. The method according to claim 1, characterized in that The process of determining the equivalent temperature of the zircon sample in each time interval based on the thermal history curve includes: Based on the thermal history curve, determining the temperature of the first endpoint and the temperature of the second endpoint of the zircon sample at each time interval; An average of the temperature at the first end point and the temperature at the second end point of the time interval is determined as the equivalent temperature of the time interval.

3. The method according to claim 1, characterized in that The process of determining the track length ratio of the zircon sample in each time interval based on the annealing model and the equivalent temperature of each time interval includes: Based on the annealing model and the equivalent temperature of the first time interval, determining the track length ratio of the zircon sample in the first time interval, wherein the first time interval is the most recent time interval; Determining an equivalent constant time of the previous time interval relative to the next time interval based on the annealing model, the track length ratio of the previous time interval, and the equivalent temperature of the next time interval; Determining the track length ratio of the zircon sample in the latter time interval based on the annealing model, the equivalent constant time and the equivalent temperature in the latter time interval; The previous time interval and the next time interval are adjacent time intervals, and the time of the previous time interval is closer to the current time.

4. A zircon fission track visual age determination device, characterized in that: include: A model determination unit, used to determine an annealing model of zircon fission tracks, wherein the annealing model records the constraint relationship between track length ratio, time and temperature; The annealing model of the zircon fission track is: ； in, , , , , and is a preset constant, For time, is the temperature, is the track length ratio; A ratio determination unit, for determining the track length ratio of the zircon sample in each time interval within a target simulation period based on the annealing model and a preset thermal history curve, wherein the target simulation period is a period where the thermal history curve is located; The time intervals are equal-length, non-overlapping, and sequentially adjacent time periods within the target simulation period; based on the annealing model and the preset thermal history curve, the process of determining the track length ratio of the zircon sample in each time interval within the target simulation period includes: Based on the thermal history curve, determining the equivalent temperature of the zircon sample in each time interval; Based on the annealing model and the equivalent temperature of each time interval, determining the track length ratio of the zircon sample in each time interval; A density determination unit, for determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period; The process of determining the track density of the zircon sample in each time interval within the target simulation period based on the track length ratio of the zircon sample in each time interval within the target simulation period includes: The following equation is used to convert the track length ratio of the zircon sample in each time interval within the target simulation period into the track density: ； in, The track length ratio of the zircon sample in the i-th time interval within the target simulation period, For the corresponding The track density of An age determination unit, for determining the fission track apparent age of the zircon sample based on the track density of the zircon sample at each time interval within the target simulation period; The process of determining the fission track apparent age of the zircon sample based on the track density of the zircon sample at each time interval within the target simulation period includes: The fission track apparent age of zircon samples was calculated using the following equation: : ； Where k is the preset correction factor, is the length of the time interval.

5. A zircon fission track visual age determination device, characterized in that: include: Memory and processor; The memory is used to store programs; The processor is used to execute the program to implement each step of the method for determining the zircon fission track apparent age as described in any one of claims 1 to 3.

6. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, each step of the method for determining the apparent age of zircon fission tracks as described in any one of claims 1 to 3 is implemented.

Images