Method for quantitatively detecting mineral resource denudation preservation condition based on thermal chronology

Abstract

The invention relates to the technical field of mineral resource exploration in geology, in particular to a method for quantitatively detecting the denudation preservation condition of mineral resources based on thermal chronology. Constructing a thermal evolution history of the target ore deposit in a set time period based on thermal chronology; dividing a stable cooling time period of the target ore deposit in an evolution history in a set time period; geologic characteristic parameters and thermal chronology temperature parameters of the target ore deposit in the stable cooling time period are obtained, a thermal simulation curve is established according to the thermal chronology temperature parameters, and a second rise value and a second stripping value of the target ore deposit in the stable cooling time period are obtained according to the geologic characteristic parameters and the thermal simulation curve; and carrying out quantitative analysis on the resource preservation condition of the target ore deposit according to the first rising value and the first stripping value of the target ore deposit. According to the method provided by the invention, denudation and uplift degree analysis of the ore deposit is carried out on the basis of thermal chronology data, accurate analysis data can be obtained, and quantitative evaluation of the storage condition of the ore deposit can be better realized.

Description

Technical Field

[0001] The invention relates to the technical field of mineral resource exploration in applied geology, and in particular to a method for quantitatively detecting the denudation and preservation of mineral resources based on thermochronology. Background Art

[0002] The mineralization system includes all geological factors and mineralization dynamic processes that control the formation and preservation of ore deposits, mainly including the activation of mineralizing materials from the source area ("source"), migration ("transport"), accumulation in highly enriched forms ("aggregation"), and the role of preserving the mineralized products in the subsequent geological history ("preservation"); the end of the mineralization process "source", "transport" and "aggregation" does not represent the end point of the evolution of the mineralization system. The "preservation" after mineralization, that is, the preservation of the ore deposit, is also an important part of the mineralization system; the erosion (exposure and erosion) caused by tectonic uplift has the most significant impact on the preservation of ore deposits. If a ore deposit has experienced a large-scale tectonic uplift after its formation and has emerged on the surface, its current surface outcrop is only the remnant of the erosion, and most of it has been destroyed by erosion, then the deep exploration potential of the ore deposit is relatively small; conversely, if a ore deposit (especially a large-super-large ore deposit) is less uplifted by the tectonic uplift, its field outcrop is only the top of the ore body, then its deep exploration potential is huge.

[0003] Traditional methods such as geological mapping and structural analysis of mining areas, experimental simulation, geophysical methods (electrical and magnetic anomalies and seismic exploration), ore deposits, mineralogy, alteration petrology and geochemical methods (such as alteration models, element zoning) and mineralization models can usually only qualitatively or semi-quantitatively study the degree of ore exposure, which is not enough to accurately and quantitatively determine the specific amount of erosion. Therefore, how to obtain accurate data on the amount of ore erosion is an important technical problem that needs to be solved in the field of ore preservation research. Summary of the invention

[0004] In order to solve the problems existing in the prior art, the present invention provides a method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology. The scheme includes: constructing the thermal evolution history of the target ore deposit within a set time period based on thermochronology; dividing the stable cooling period of the target ore deposit in the evolution history within the set time period; obtaining the geological characteristic parameters and thermochronological temperature parameters of the target ore deposit in the stable cooling period, establishing a thermal simulation curve based on the thermochronological temperature parameters, and obtaining the second uplift value and the second denudation value of the target ore deposit in the stable cooling period according to the geological characteristic parameters and the thermal simulation curve; quantitatively analyzing the resource preservation of the target ore deposit according to the first uplift value and the first denudation value of the target ore deposit. The method proposed in this application analyzes the degree of denudation and uplift of the ore deposit based on thermochronological data, and can obtain accurate analysis data, thereby better realizing the quantitative evaluation of the preservation of the ore deposit.

[0005] The present invention adopts the following technical scheme, a method for quantitatively detecting the denudation and preservation of mineral resources based on thermochronology, comprising:

[0006] Construct the thermal evolution history of the target deposit within a set time period based on thermochronology;

[0007] Dividing a stable cooling period of the target ore deposit within the thermal evolution history within the set time period;

[0008] Acquiring thermochronological temperature parameters and geological characteristic parameters of the target ore deposit during the stable cooling period, and establishing a thermal simulation curve according to the thermochronological temperature parameters;

[0009] Obtaining a first uplift value and a first stripping value of the target ore deposit during a stable cooling period according to the thermal simulation curve and the geological characteristic parameters;

[0010] A quantitative analysis of the resource preservation of the target deposit is conducted based on the first uplift value and the first denudation value of the target deposit.

[0011] Furthermore, the method further comprises:

[0012] According to the geological characteristic parameters, a closed temperature method is used to obtain a second uplift value and a second stripping value of the target ore deposit during a stable cooling period;

[0013] Obtaining the thermochronological age of the target deposit based on the evolutionary history;

[0014] establishing an age-elevation diagram of the target deposit based on the thermochronological age and the elevation value of the target deposit;

[0015] modifying the second uplift value and the second denudation value based on the age-elevation diagram;

[0016] A quantitative analysis of the resource preservation of the target deposit is conducted based on the revised second uplift value and second stripping value.

[0017] Furthermore, obtaining the thermochronological temperature parameters of the target ore deposit during the stable cooling period includes: using thermochronological simulation software to read the burial temperature of the target ore deposit.

[0018] Furthermore, obtaining the first uplift value and the first stripping value of the target ore deposit during the stable cooling period includes:

[0019] The burial temperature of the target ore deposit is read according to the thermal simulation curve, and the expression of the first stripping value is constructed in combination with the geological characteristic parameters:

[0020] ΔE=(T m -T s ) / G

[0021] The expression of the first uplift value is:

[0022] ΔH=U-ΔE

[0023] In the above formula, T m represents the burial temperature of the target deposit, T s represents the ancient surface temperature, and G is the ancient geothermal gradient.

[0024] Furthermore, the thermal evolution history of the target deposit within a set time period is constructed based on thermochronology, including:

[0025] Set the starting time of the nth temperature change period experienced by the target deposit as t n , the corresponding temperature is T n ; The end time is t n+1 , corresponding to temperature T n+1 ; Wherein, n is a positive integer not less than 1;

[0026] If we take the magmatic deposit in the period [t n , t n+1 ) is ΔT, and its corresponding time period is Δt. Then the thermal evolution history of the target deposit in this period is established as: ΔT = T n -T n+1 =∫C r dt, where C r is the cooling rate within the time period Δt.

[0027] Furthermore, the geological characteristic parameters include: paleo-geothermal gradient, paleo-surface temperature, paleo-average surface elevation and present-day average surface elevation.

[0028] Furthermore, obtaining the second uplift value and the second stripping value of the target ore deposit during the stable cooling period includes:

[0029] The calculation steps for obtaining the second stripping value are:

[0030] Get the burial depth of the closed temperature range of the paleothermal chronology system:

[0031] D=(T c -T s ) / G

[0032] Among them, T c represents the closing temperature of the thermochronological system, T s represents the paleosurface temperature, G represents the paleogeothermal gradient;

[0033] The parameter d is set to represent the distance between the burial position of the target deposit and the closed temperature range of the thermochronological system. At this time, the degree of exposure is expressed as:

[0034] ΔE=D±d

[0035] The second stripping value is obtained by comprehensive analysis:

[0036] ΔE±d=(T c -T s ) / G

[0037] The steps of obtaining the second uplift value are:

[0038] The relationship between rock uplift amplitude, surface uplift amplitude and exfoliation degree is established:

[0039] H p =(U+H i )+ΔE

[0040] Among them, the rock uplift amplitude is (U+H i ), ΔE represents the degree of exfoliation, H p is the surface uplift amplitude, H i represents the paleo-average surface height, and U represents the uplift amplitude;

[0041] Let ΔH=H p -H i , we can get:

[0042] ΔH=U-ΔE

[0043] Then the expression of the second uplift value is: U±d=ΔH+(ΔE±d).

[0044] Furthermore, the second uplift value and the second denudation value are corrected based on the age-elevation diagram, specifically:

[0045] The absolute value of the elevation at which the trend line of the target mine in the age-elevation diagram turns is used as the value of the parameter d, and the corrected second uplift value and the second stripping value are obtained according to the value of d at this time.

[0046] The beneficial effects of the present invention are as follows: the technical means proposed in the present invention is based on thermochronological data, and can conduct comparative analysis on the degree of exposure and emplacement depth of the mine field, so as to better evaluate the preservation of this type of ore deposit after its formation, and provide more solid theoretical support for mineral exploration work at the mining area and regional scales, the measurement of the economic value of ore-bearing rock masses, and the study of the uplift and exposure process of magmatic deposits. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the embodiments of the present invention 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 some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0048] Figure 1 A schematic flow chart of a method for quantitatively detecting the denudation and preservation of mineral resources based on thermochronology according to an embodiment of the present invention;

[0049] Figure 2 A schematic diagram of the thermal evolution history of a sample obtained by simulation using a thermal age simulation software according to an embodiment of the present invention;

[0050] Figure 3 A schematic diagram of a cooling rate calculation principle curve according to an embodiment of the present invention;

[0051] Figure 4 A schematic diagram of the relationship between temperature, elevation and time during an uplift-stripping process according to an embodiment of the present invention;

[0052] Figure 5 A schematic diagram of a Jinchuan mining area thermochronology simulation result according to an embodiment of the present invention;

[0053] Figure 6 An age-sample elevation diagram according to an embodiment of the present invention. DETAILED DESCRIPTION

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

[0055] In the present embodiment, the term "thermochronological data" is first explained: thermochronology is an isotopic dating method that determines the cooling and exhumation age (i.e., thermochronological age) of a geological body based on the formation and accumulation of products of processes such as radioactive isotope decay or fission in rocks or minerals. Thermochronology includes (U-Th) / He dating, fission track (FT) analysis and 40 Ar / 39Ar dating and other main methods. Each thermochronological system (rock or mineral + method) has its own unique thermal diffusion (or annealing) model and closed temperature range, such as the closed temperature range of the apatite fission track (AFT) system is 120-60℃; when the ambient temperature changes and is lower than the closed temperature, the thermochronological system gradually transforms from a completely open isotope system to a completely closed system, and records the thermochronological age; that is, the thermochronological age reflects the time information of the cooling process of "the latest experience of the system closure temperature", so each thermochronological age contains the triple connotations of time, space (temperature / depth) and geology; if the basic data such as the age of multiple thermochronological systems and the length of fission tracks can be combined with the fission track annealing model and the helium diffusion model to perform thermal evolution history inversion modeling, the time-temperature thermal evolution history or the time-depth denudation process can be quantitatively restored.

[0056] A flow chart of a method for quantitatively detecting the denudation and preservation of mineral resources based on thermochronology according to an embodiment of the present invention is shown in FIG. Figure 1 As shown, including:

[0057] Based on the existing thermochronological data, combined with existing literature and other information, the thermal evolution history of the target deposit within a certain period is constructed with the help of thermochronological simulation software.

[0058] In a specific embodiment of the present invention, the thermal evolution history of a magmatic deposit (or other related geological bodies such as a rock mass, stratum, etc.) is a process in which its temperature changes within a relatively certain time range. The process changes continuously in time and temperature, and there is a unique correspondence between time and temperature. Consider that the start time of the nth temperature change period experienced by the magmatic deposit is t n , the corresponding temperature is T n ; The end time is t n+1 , corresponding to temperature T n+1 ; Where n is a positive integer not less than 1; if the magmatic deposit is taken in the period [t n , t n+1 ) is ΔT, and its corresponding time period is Δt. Then the thermal evolution history of the target deposit in this period is established as: ΔT = T n -T n+1 =∫C r dt, where C r is the cooling rate within the time period Δt, such as Figure 3 As shown, C is usually used r =ΔT / Δt to represent the average cooling rate of the deposit within the time period Δt.

[0059] Based on the thermal evolution history obtained, identify the geological events that may have led to stable cooling of the deposit and delineate the time periods when these events had an impact;

[0060] Through existing literature and other information, the geological characteristic parameters corresponding to the mineral deposits during each period of stable cooling are collected, including paleo-geothermal gradient (commonly used unit: ℃ / km, the same below), paleo-surface temperature (℃), paleo-average surface elevation (km), and current average surface elevation (km). The number of decimal places of all data are unified according to the calculation accuracy requirements.

[0061] Acquiring thermochronological temperature parameters and geological characteristic parameters of the target ore deposit during the stable cooling period, and establishing a thermal simulation curve according to the thermochronological temperature parameters;

[0062] like Figure 2 The time-temperature curve shown is the thermal evolution history given by the HeFTy software in the form of 'Pathenvelopes' based on the sample and the limit frame. If you want to obtain the burial temperature of the sample at 240Ma, draw a straight line perpendicular to the horizontal axis and intersect the figure formed by the "high-quality simulation path" at two points, and then project these two points on the vertical axis (temperature). According to the projection, it can be read that the sample temperature range corresponding to the thermochronological age of the sample is 10-47°C. When formally calculating, this temperature range should consider other factors to make necessary corrections. It should be noted that the premise of applying the thermal simulation curve method in the embodiment of the present invention is that the thermal evolution history has been accurately portrayed, and the accuracy of this method is directly related to the accuracy of the thermal evolution history.

[0063] In another embodiment of the present invention, the thermal evolution history can also be simulated by QTQt software using the MCMC (Markov Chain Monte Carlo) method. It should be noted that the software for thermal evolution simulation can also be any of the existing software such as Low-T Thermo or PQ software. The present invention aims to analyze the curve obtained after the simulation, and there is no limitation on how to simulate the curve.

[0064] Obtaining a first uplift value and a first stripping value of the target ore deposit during a stable cooling period according to the thermal simulation curve and the geological characteristic parameters;

[0065] In the embodiment of the present invention, the burial temperature of the target ore deposit is read out according to the established thermal simulation curve, and the expression of the first stripping value is constructed in combination with the geological characteristic parameters as follows:

[0066] ΔE=(T m -T s ) / G

[0067] The expression of the first uplift value is:

[0068] ΔH=U-ΔE

[0069] In the above formula, T m represents the burial temperature of the target deposit, T s represents the ancient surface temperature, and G is the ancient geothermal gradient.

[0070] In the embodiments of the present invention, other geological information (such as the emplacement depth of the ore-bearing rock mass, geometric elements, etc.) can be combined for analysis and interpretation. At the same time, text or pictures can be used to describe the degree of exposure, preservation, development potential, especially the prospects for deep mineral exploration of the target ore deposit.

[0071] In another specific embodiment of the present invention, the method for quantitatively detecting the erosion and preservation of mineral resources further comprises:

[0072] According to the geological characteristic parameters, a closed temperature method is used to obtain a second uplift value and a second stripping value of the target ore deposit during a stable cooling period;

[0073] Obtaining the thermochronological age of the target deposit based on the evolutionary history;

[0074] establishing an age-elevation diagram of the target deposit based on the thermochronological age and the elevation value of the target deposit;

[0075] modifying the second uplift value and the second denudation value based on the age-elevation diagram;

[0076] A quantitative analysis of the resource preservation of the target deposit is conducted based on the revised second uplift value and second stripping value.

[0077] Specifically, the current study discusses three types of displacements (the direction of the gravity vector is defined as the negative direction) based on the relationship between the point uplift in the region and the overall uplift: Figure 4 As shown: 1. Rock uplift, denoted as (U+H i ); 2. Exhumation, denoted as ΔE, with the same direction as (U+H i ) is opposite; 3. The surface uplift (surface uplift), denoted as H p The starting points of the two displacements except ΔE are all located on the geoid. The relationship between the above three displacements is:

[0078] H p =(U+H i )+ΔE

[0079] That is H p -H i =U+ΔE, which means that the difference between "uplift amplitude" and "exfoliation degree" is equivalent to the difference between "present surface elevation" and "ancient surface elevation" in value; among them, the rock uplift amplitude is (U+Hi ), ΔE represents the degree of exfoliation, H p is the surface uplift amplitude, H i represents the paleo-average surface height, and U represents the uplift amplitude;

[0080] Let ΔH=H p -H i , we can get:

[0081] ΔH=U-ΔE

[0082] In actual calculations, it is usually necessary to consider the amplitude of sea level rise and fall, which reflects the elevation transformation of the geoid itself. In the embodiment of the present invention, it is set to 0;

[0083] Assume that the closing temperature of the thermochronological system is T c ; The surface temperature is T s ; The paleo-geothermal gradient during this period is G. Considering that the uplift amplitude U cannot be directly limited, the burial depth D of the closed temperature range of the thermochronological system can be calculated according to the following formula:

[0084] D=(T m -T s ) / G

[0085] Since there is still a certain distance between the buried position of the ore deposit and the closed temperature range of the thermochronological system, the embodiment of the present invention introduces d, and ΔE=D±d

[0086] Therefore, the second stripping value is obtained comprehensively:

[0087] ΔE±d=(T c -T s ) / G

[0088] Substituting ΔH=U-ΔE, the expression of the second uplift value can be obtained as follows: U±d=ΔH+(ΔE±d);

[0089] In practical applications, d ≥ 0 should always be true. p <(H i -D)+U is established, which means that when the thermochronological age of the sample is recorded, the rock mass at the closed temperature range of the corresponding system has been completely eroded, and at this time ΔE=(T c -T s ) / G+d; Similarly, if there is H p >(H i -D)+U, then ΔE=(T c -T s ) / Gd.

[0090] In order to more accurately obtain the uplift amplitude and exposure degree of the rock mass, it is necessary to limit the value of d. The embodiment of the present invention makes a thermochronological age-sample elevation diagram, such as Figure 6 As shown in Figure 2, d is defined by the absolute value of the elevation corresponding to the turning point of the trend line of the sampling point. The samples above the turning point are located in the partial annealing zone before the uplift, where they have accumulated fission tracks (including early partial annealing tracks and late tracks accumulated during the uplift) and recorded thermochronological ages. The samples below the turning point only begin to accumulate fission tracks when they uplift and exceed the elevation of the bottom of the partial annealing zone. The sudden change in slope marks the beginning of a regional uplift event.

[0091] In another specific embodiment of the present invention, taking the Jinchuan deposit as an example, a simple quantitative analysis is performed on the degree of mineral resource exposure, preservation, and development potential of the Jinchuan deposit; and the good simulation curves of the results of the thermochronological simulation of the Jinchuan mining area in various literature materials are sorted into the same Tt diagram, such as Figure 5 As shown in the figure, it can be seen that after the Silurian magmatic activity, the exposure cooling stage of the Jinchuan mining area began as early as the Carboniferous, and thereafter no longer experienced obvious warming events, and its thermal evolution history can be regarded as a linear cooling process. Next, we focus on the uplift-exhumation history of the Jinchuan II mining area.

[0092] Since there is no research on the ancient geothermal gradient of the Jinchuan mining area, according to existing research, in the Late Jurassic-Early Cretaceous, the geothermal gradient of Huashen 1 well in the Huahai Basin of the Jiuquan Basin Group reached 37.5-45.0℃ / km, and that of Dan 25 well in the Jiuxi Basin was 42℃ / km. At the same time, according to the mining area data, the current geothermal gradient of the Jinchuan mining area is about 19.8°C / km. Therefore, when the embodiment of the present invention calculates the amount of erosion in each period of the Jinchuan mining area, the geothermal gradient during 320-160Ma is taken as 45°C / km; the geothermal gradient during 160-100Ma is 40°C / km; the geothermal gradient during 100-70Ma is taken as 30°C / km; the geothermal gradient from 70Ma to the present is 20°C / km; during the Carboniferous-Permian period, the fault basin in the Longshou Mountain area received marine or marine-continental transition carbonate rock-clastic rock-volcanic rock deposition, so the Longshou Mountain area should have been a marine or marine-continental transition depositional environment at that time, and the ancient average surface elevation was taken as -0.2km; the elevation of the mining area at 260Ma was taken as 0.5km; the elevation at 70Ma was 0.5km, and the current surface elevation was taken as 1.6km.

[0093] The samples currently available in the Jinchuan mining area have similar sampling elevations, and the thermochronological age-sampling elevation diagram cannot be used to limit d. In addition, there are currently no research results related to the paleogeothermal field in the Jinchuan mining area for reference. Therefore, the "thermal simulation curve method" is used in the embodiment of the present invention to calculate the first uplift value and the first denudation value of the Jinchuan mining area.

[0094] The values ​​of ΔE in each stable cooling period of Jinchuan II mining area are obtained through analysis and calculation as shown in Table 1:

[0095] Table 1. Excavation thickness of Jinchuan II mining area at different time periods

[0096]

[0097] “Exposed” in the table refers to the time when the sample sampling location (not the top of the rock mass) is exposed to the surface.

[0098] The uplift amplitude of Jinchuan II mining area in different periods is calculated as follows: ① During the Late Paleozoic (about 300-260Ma), the exposure degree of Jinchuan II mining area ΔE = 1.8km. If the elevation of the mining area at 300Ma is -0.2km and the elevation at 260Ma is 0.5km, then ΔH = 0.7km, and the uplift amplitude of the rock mass during this period is U = 2.5km; ② The exposure degree of the ore-bearing rock mass in Jinchuan II mining area during the Late Paleozoic-Late Mesozoic (about 260-70Ma) ΔE = 1.9km. If the elevation of the mining area at 70Ma is 0.5km, then ΔH = 0 and U = 1.9km during this period; ③ Since 70Ma, the exposure degree of the ore-bearing rock mass in Jinchuan II mining area ΔE = 1.9km. If the current average elevation of the mining area is taken as 1.6km, the ΔH of the rock mass at this stage is 1.1km, and the uplift amplitude is U=3.0km.

[0099] Comprehensive discussion of the uplift and denudation history of the Jinchuan II mining area since 300Ma: ① At 300Ma, the current sampling location (the middle of the ore-bearing rock mass in the II mining area) was about 5.5km from the ground. During the Late Paleozoic (about 300-260Ma), the uplift amplitude of the ore-bearing rock mass in the Jinchuan II mining area was U = 2.5km, and the degree of denudation was about 1.8km; ② At 260Ma, the current sampling location was about 3.7km from the ground. The uplift amplitude of the ore-bearing rock mass in the Jinchuan II mining area during the Late Paleozoic-Late Mesozoic (about 260-70Ma) was about 1.9km, and the degree of denudation was about 1.9km; ③ Since 70Ma, the ore-bearing rock mass in the Jinchuan II mining area has continued to uplift and exposed to the surface. The uplift amplitude of the rock mass during this stage was U = 3.0km. Until the current sampling location was exposed, the rock mass had suffered The degree of exposure is about 1.9km. Therefore, after comprehensively analyzing the current geological, geochemical and geophysical data, it is believed that the rock mass in the Jinchuan II mining area extends about 800-1200m underground. If the average value is 1000m, combined with the exposure degree of 1900m, it can be seen that the scale of the existing Jinchuan II mining area rock mass is about 34.5% of the initial rock mass; according to the mechanism of (large) mineralization of small mafic rock masses, ore-rich magma and ore pulp are mostly located at the bottom or tail of the rock mass, reflecting the characteristics of tail feather mineralization. The larger degree of exposure of the ore-bearing rock mass in the Jinchuan II mining area is beneficial to further deep prospecting to a certain extent.

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

Claims

1. A method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology, characterized in that: include: Construct the thermal evolution history of the target deposit within a set time period based on thermochronology; Dividing a stable cooling period of the target ore deposit within the thermal evolution history within the set time period; Acquiring thermochronological temperature parameters and geological characteristic parameters of the target ore deposit during the stable cooling period, and establishing a thermal simulation curve according to the thermochronological temperature parameters; Obtaining a first uplift value and a first stripping value of the target ore deposit during a stable cooling period according to the thermal simulation curve and the geological characteristic parameters; A quantitative analysis of the resource preservation of the target deposit is conducted based on the first uplift value and the first denudation value of the target deposit.

2. The method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology according to claim 1, characterized in that: The method further comprises: According to the geological characteristic parameters, a closed temperature method is used to obtain a second uplift value and a second stripping value of the target ore deposit during a stable cooling period; Obtaining the thermochronological age of the target deposit based on the evolutionary history; establishing an age-elevation diagram of the target deposit based on the thermochronological age and the elevation value of the target deposit; modifying the second uplift value and the second denudation value based on the age-elevation diagram; A quantitative analysis of the resource preservation of the target deposit is conducted based on the revised second uplift value and second stripping value.

3. The method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology according to claim 1, characterized in that: Obtaining the thermochronological temperature parameters of the target ore deposit during the stable cooling period includes: using thermochronological simulation software to read the burial temperature of the target ore deposit.

4. The method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology according to claim 1, characterized in that: Obtain the first uplift value and the first stripping value of the target deposit during the stable cooling period, including: The burial temperature of the target ore deposit is read according to the thermal simulation curve, and the expression of the first stripping value is constructed in combination with the geological characteristic parameters: ΔE=(T m -T s ) / G The expression of the first uplift value is: ΔH=U-ΔE In the above formula, T m represents the burial temperature of the target deposit, T s represents the ancient surface temperature, and G is the ancient geothermal gradient.

5. The method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology according to claim 1, characterized in that: Based on thermochronology, the thermal evolution history of the target deposit within a set time period is constructed, including: Set the starting time of the nth temperature change period experienced by the target deposit as t n , the corresponding temperature is T n ; The end time is t n+1 , corresponding to temperature T n+1 ; Wherein, n is a positive integer not less than 1; If we take the magmatic deposit in the period [t n , t n+1 ) is ΔT, and its corresponding time period is Δt. Then the thermal evolution history of the target deposit in this period is established as: ΔT = T n -T n+1 =∫C r dt, where C r is the cooling rate within the time period Δt.

6. The method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology according to claim 1, characterized in that: The geological characteristic parameters include: paleo-geothermal gradient, paleo-surface temperature, paleo-average surface elevation and present-day average surface elevation.

7. The method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology according to claim 2, characterized in that: Obtain the second uplift value and second stripping value of the target deposit during the stable cooling period, including: The calculation steps for obtaining the second stripping value are: Get the burial depth of the closed temperature range of the paleothermal chronology system: D=(T c -T s ) / G Among them, T c represents the closing temperature of the thermochronological system, T s represents the paleosurface temperature, G represents the paleogeothermal gradient; The parameter d is set to represent the distance between the burial position of the target deposit and the closed temperature range of the thermochronological system. At this time, the degree of exposure is expressed as: ΔE=D±d The second stripping value is obtained as follows: ΔE±d=(T c -T s ) / G The steps of obtaining the second uplift value are: The relationship between rock uplift amplitude, surface uplift amplitude and exfoliation degree is established: H p =(U+H i )+ΔE Among them, the rock uplift amplitude is (U+H i ), ΔE represents the degree of exfoliation, H p is the surface uplift amplitude, H i represents the paleo-average surface height, and U represents the uplift amplitude; Let ΔH=H p -H i , we can get: ΔH=U-ΔE Then the expression of the second uplift value is: U±d=ΔH+(ΔE±d).

8. The method for quantitatively detecting the erosion and preservation of mineral resources based on thermochronology according to claim 2, characterized in that: The second uplift value and the second denudation value are corrected based on the age-elevation diagram, specifically: The absolute value of the elevation at which the trend line of the target mine in the age-elevation diagram turns is used as the value of the parameter d, and the corrected first uplift value and the first stripping value are obtained according to the value of d at this time.

Images