A geological survey method for recovering protolith types of metamorphic rocks by using boron isotope fractionation characteristics

Abstract

The application provides a geological survey method for restoring protolith types of metamorphic rocks by using boron isotope fractionation characteristics. The method establishes a benchmark database containing boron isotopes and contents of key inactive elements of various unmetamorphic protoliths. A double-control correction model of boron isotope fractionation is established, quantitative function relationships between distribution coefficients and pressure and between fractionation coefficients and temperature are respectively constructed by using high-temperature and high-pressure experimental data, and a fluid activity factor is introduced to form an inversion calculation equation. The original boron isotope composition is obtained by using the model, and then the original boron isotope composition is analyzed in cooperation with the key element index to construct a two-dimensional discrimination diagram, and finally, the protolith type is determined by comparing the discrimination area in the benchmark database. The application effectively corrects the metamorphic fractionation effect, and significantly improves the accuracy and reliability of the protolith restoration of high-grade metamorphic rocks by combining the double-geochemical index constraint. The method can be matched with portable field analysis equipment to realize rapid discrimination in the field.

Description

Technical Field

[0001] This invention relates to the field of geological survey technology, and in particular to a geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics. Background Technology

[0002] Protolith reconstruction of high-grade metamorphic rocks, such as granulite facies rocks, has long been a challenging problem in geology. Traditionally widely used geochemical discrimination diagrams become significantly less reliable under high-grade metamorphic conditions due to the substantial migration of numerous elements during the metamorphic process. Therefore, finding stable geochemical indicators that can still effectively trace protolith origins in strongly metamorphic alterations has become crucial in this field.

[0003] Boron and its stable isotope δ 11 Boron (B), as a sensitive geochemical tracer, exhibits significant differences in protoliths of different origins. However, boron also possesses high fluidity during metamorphism, readily migrating and undergoing isotopic fractionation in open tectonic-thermal events. This leads to substantial bias in protolith identification using directly measured boron isotopic compositions of metamorphic rocks, as the signal is a superposition of protolith signals and the effects of multiple metamorphic fractionation phases.

[0004] Currently, there is a lack of systematic and effective technical means to quantitatively isolate the alteration of boron isotope composition by metamorphic processes, thereby extracting original boron isotope information that can represent the protolith type. Existing research has neither established a sufficient protolith boron isotope reference database, nor does it have a correction model that can simultaneously simulate the effects of pressure conditions and fluid activity on fractionation, nor has it formed a complete methodological system integrating correction and multivariate discrimination. This technological gap severely limits the practical application of boron isotope tracing technology in the protolith reconstruction of advanced metamorphic terranes.

[0005] Therefore, there is an urgent need to develop a new geological survey method that can systematically correct the alteration of boron isotope composition by metamorphism and synergistically utilize stable geochemical indicators to achieve reliable restoration of the protolith types of advanced metamorphic rocks. Summary of the Invention

[0006] To address the shortcomings of the existing technologies, this invention proposes a geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics.

[0007] The present invention specifically provides the following technical solution:

[0008] A geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, the method comprising the following steps:

[0009] S1. Collect and determine the boron content, boron isotope composition and key inactive element content of various known unmetamorphosed protoliths, and establish a benchmark database including marine sedimentary rocks, oceanic crust basalt, island arc magmatic rocks and terrigenous clastic rocks;

[0010] S2. Establish a dual-control calibration model for boron isotope fractionation. The dual-control calibration model is used to simulate the distribution behavior of boron between the mineral phase and the fluid phase and the boron isotope fractionation effect under metamorphic conditions. The dual-control calibration model uses the pressure conditions and fluid activity during the metamorphic process as the core control variables.

[0011] S3. Obtain the measured whole-rock boron isotope composition data and key inactive element content data of the metamorphic rock sample to be tested, and determine the temperature and pressure conditions of the main metamorphic stages experienced by the sample.

[0012] S4. Based on the determined temperature and pressure conditions, the measured whole-rock boron isotope composition data are fractionated and corrected using the dual-control correction model, and the original boron isotope composition of the sample before metamorphism is obtained by inversion calculation.

[0013] S5. The original boron isotope composition obtained from the inversion calculation is analyzed in conjunction with the content data of the key inactive elements to form a comprehensive discrimination index;

[0014] S6. Compare the comprehensive discrimination index with the benchmark database, and determine the protolith type of the metamorphic rock sample to be tested according to the preset dual-index discrimination rule.

[0015] Optionally, in step S1, establishing the benchmark database specifically includes:

[0016] Marine sedimentary rocks, oceanic basalts, island arc igneous rocks, and terrigenous clastic rocks with clear geological backgrounds and no significant metamorphic alteration were selected as standard samples.

[0017] The whole-rock boron content and whole-rock boron isotope δ¹⁸O of the standard sample were determined. 11 Composition B, and for the aforementioned marine sedimentary rock standard samples, the single-mineral boron isotope δ¹⁸ of representative boron-bearing minerals in the samples were additionally measured. 11 B value;

[0018] The measured boron content data and whole-rock boron isotope δ¹⁸O will be used to... 11 B value, boron isotope δ of single mineral 11 The B-value is correlated with the content data of key inactive elements corresponding to the standard sample and geological background information to construct a standardized benchmark database.

[0019] Optionally, in step S2, establishing the dual-control correction model includes establishing a quantitative functional relationship between the distribution coefficient and the pressure, wherein the quantitative functional relationship is determined in the following manner:

[0020] Based on high-temperature and high-pressure experiments simulating geological metamorphic processes, data on the distribution coefficients of boron between the target metamorphic minerals and coexisting fluid phases under different pressure conditions were obtained. A functional relationship between the distribution coefficient D and the pressure P was derived through nonlinear regression fitting. Specifically, the functional relationship is as follows:

[0021]

[0022] in, , , For the fitting parameters, For pressure.

[0023] Optionally, establishing the dual-control correction model further includes establishing a quantitative functional relationship between the boron isotope fractionation coefficient and temperature, and constructing an inversion calculation equation, wherein: the quantitative functional relationship between the fractionation coefficient and temperature is determined in the following way:

[0024] Based on high-temperature and high-pressure experiments simulating geological metamorphic processes, boron isotope equilibrium fractionation values ​​between target metamorphic minerals and fluid phases under different temperature conditions were obtained. A nonlinear regression fitting was used to derive a functional relationship between the fractionation value Δ and temperature T. Specifically, the functional relationship is as follows:

[0025]

[0026] Where A, B, and C are fitting parameters, and T is the absolute temperature;

[0027] The inversion calculation equation is:

[0028]

[0029] in, It has the original boron isotopic composition. To measure the whole-rock boron isotope composition, As a fluid activity factor, For allocation coefficients, This is the fractionation value.

[0030] Optionally, in S3, the acquisition of measured whole-rock boron isotope composition data is performed by using a multi-receiver inductively coupled plasma mass spectrometer to measure the powder pellets or fused glass targets prepared from the metamorphic rock sample to be tested.

[0031] The determination of the temperature and pressure conditions of the main metamorphic stages experienced by the sample is specifically based on the petrographic observation results of the metamorphic rock sample to be tested, and quantitative calculation is performed using a symbiotic mineral geothermobarometer; wherein, the symbiotic minerals include typical symbiotic assemblages of target boron-rich minerals, and the temperature and pressure conditions are calculated based on the chemical composition of the characteristic minerals in the symbiotic assemblage.

[0032] Optionally, determining the temperature and pressure conditions for the main metamorphic stages experienced by the sample specifically includes:

[0033] Based on petrographic observations of the metamorphic rock samples to be tested, specific mineral pairs associated with the target boron-rich minerals are identified; wherein, the specific mineral pairs include at least one of the following: a garnet-biotite mineral pair for temperature calculation, or a garnet-orthopyroxene-plagioclase-quartz mineral assemblage for pressure calculation; the chemical composition of the specific mineral pairs is determined by electron probe microanalysis.

[0034] The chemical composition data is input into the corresponding geological thermo-barometer calculation formula to calculate the temperature and pressure values ​​respectively.

[0035] Optionally, the dual-control correction model is invoked to perform fractionation correction and inversion calculations, specifically including:

[0036] The values ​​of metamorphic temperature and metamorphic pressure determined in S3 are input into the dual-control correction model to calculate the corresponding boron isotope fractionation coefficient Δ and mineral-fluid partition coefficient D, respectively.

[0037] Based on the petrographic analysis results of the metamorphic rock sample to be tested, the fluid activity factor in the dual-control calibration model was set. The value;

[0038] The measured whole-rock boron isotope composition obtained in S3 The calculated Δ and D, and the set The inversion calculation equation is input and solved to output the original boron isotope composition. .

[0039] Optionally, in S5, the collaborative analysis and formation of comprehensive discrimination indicators is specifically achieved by constructing a "boron isotope-key element" collaborative discrimination system, including:

[0040] The original boron isotope composition obtained from the inversion calculation is used as the core isotope index;

[0041] From the obtained key inactive element content data, at least one element ratio that is significant in distinguishing the original rock type is calculated as a key element index. The element ratio is composed of elements that are chemically stable and have strong migration inertness during advanced metamorphism. The key inactive element is at least two of Ti, Nb, Ta, Zr, Hf, and Th.

[0042] A synergistic correspondence is established between the core isotope index and the key element index. A two-dimensional discriminant diagram is constructed with the core isotope index as the horizontal axis and the key element index as the vertical axis. The comprehensive discriminant index corresponds to the coordinate point in the two-dimensional discriminant diagram.

[0043] Optionally, in step S6, the step of judging according to a preset dual-index rule specifically includes:

[0044] In the two-dimensional discrimination diagram, a specific discrimination region is defined for the standard sample data of each type of protolith in the benchmark database. The boundary of the discrimination region is determined by the statistical distribution range of the core isotope index and the key element index of the corresponding type of protolith.

[0045] If the coordinates of the metamorphic rock sample to be tested fall within the discrimination area of ​​a certain type of protolith, then the protolith of the sample is determined to belong to the corresponding category.

[0046] Optionally, determining the protolith type of the metamorphic rock sample to be tested specifically includes:

[0047] Calculate the Euclidean distance from the coordinates of the metamorphic rock sample to be tested to the center point of each type of protolith discrimination region in the benchmark database;

[0048] The metamorphic rock sample to be tested is identified as the protolith type corresponding to the discrimination region with the smallest Euclidean distance.

[0049] This invention offers the following beneficial technical effects: It provides a geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics. By establishing a boron isotope reference database, this invention constructs a dual-control calibration model with pressure conditions and fluid activity as core control variables, and establishes a "boron isotope-key element" collaborative discrimination system. This significantly improves the accuracy and reliability of protolith reconstruction for advanced metamorphic rocks, and has important value for regional geological surveys and resource assessment. Attached Figure Description

[0050] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0051] Figure 1 This is a schematic diagram of a geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, provided as an embodiment of the present invention.

[0052] Figure 2 A conceptual model diagram of a modified system provided for embodiments of the present invention.

[0053] Figure 3 The core calculation relationship diagram of the boron isotope fractionation dual-control correction model provided for the embodiments of the present invention.

[0054] Figure 4 A diagram illustrating the synergistic discrimination of boron isotopes and key elements provided for embodiments of the present invention. Detailed Implementation

[0055] The embodiments of this application will now be described in detail with reference to the accompanying drawings.

[0056] The following specific examples illustrate the implementation of this application. Those skilled in the art can easily understand other advantages and effects of this application from the content disclosed in this specification. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of them. This application can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of this application. It should be noted that, in the absence of conflict, the following embodiments and features in the embodiments can be combined with each other. Based on the embodiments in this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0057] It should be noted that various aspects of embodiments within the scope of the appended claims are described below. It will be apparent that the aspects described herein can be embodied in a wide variety of forms, and the structure and / or function of any markings described herein are merely illustrative. Based on this application, those skilled in the art will understand that one aspect described herein can be implemented independently of any other aspect, and two or more of these aspects can be combined in various ways. For example, any number and aspects set forth herein can be used to implement the apparatus and / or practical methods.

[0058] It should also be noted that the illustrations provided in the following embodiments are only schematic representations of the basic concept of this application.

[0059] Additionally, specific details are provided in the following description to facilitate a thorough understanding of the examples. However, those skilled in the art will understand that practice can be carried out without these marked details.

[0060] The purpose of this invention is to provide a geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, establishing a quantitative correlation model between boron isotopes and protolith type, proposing a pressure-fluid dual-control model for boron isotope fractionation during metamorphism, and realizing real-time protolith identification in the field.

[0061] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0062] All formula calculations involved in this invention are presented after parameter normalization and dimension removal.

[0063] This example uses the protolith restoration of boron-rich, high-grade metamorphic rocks in the Larsmann Hills region of East Antarctica. (Refer to...) Figure 1 This paper illustrates a geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, according to an embodiment of this application. The method includes the following steps: S1. Collecting and measuring the boron content, boron isotope composition, and key inactive element content of various known unmetamorphosed protoliths, and establishing a benchmark database including marine sedimentary rocks, oceanic basalt, island arc magmatic rocks, and terrigenous clastic rocks.

[0064] In S1, establishing the benchmark database specifically includes: selecting marine sedimentary rocks, oceanic basalts, island arc igneous rocks, and terrigenous clastic rocks with clear geological backgrounds and no significant metamorphic alteration as standard samples; and determining the whole-rock boron content and whole-rock boron isotope δ¹⁸O of the standard samples. 11 Composition B, and for the aforementioned marine sedimentary rock standard samples, the single-mineral boron isotope δ¹⁸ of representative boron-bearing minerals in the samples were additionally measured. 11 B value; the measured boron content data and whole-rock boron isotope δ¹⁸O values. 11 B value, boron isotope δ of single mineral 11 The B-value is correlated with the content data of key inactive elements corresponding to the standard sample and geological background information to construct a standardized benchmark database.

[0065] Specifically, typical protoliths with clear geological backgrounds and no significant metamorphic alteration were selected as standard samples. These included marine sedimentary rocks such as organic-rich shale and mudstone known to have formed in stable continental margin environments, or boron-bearing strata within evaporite systems; oceanic basalts collected from modern mid-ocean ridges or well-preserved fragments of ancient oceanic crust, such as basalt within ophiolite suites; island arc igneous rocks including typical island arc basalts and andesites; and terrigenous clastic rocks such as quartz sandstone with high compositional maturity. Detailed geological background information, including field occurrence, geographical location, and stratigraphic unit, was recorded for each type of sample. Systematic testing was performed on the collected standard samples. All standard samples were pulverized to below 200 mesh, prepared using alkali fusion or acid dissolution methods, and their whole-rock boron content and the content of key inactive elements, including but not limited to Ti, Nb, Ta, Zr, Hf, and T, were determined using inductively coupled plasma mass spectrometry (ICP-MS). Simultaneously, the whole-rock boron isotopic composition δ¹⁸O was determined using a multi-receiver inductively coupled plasma mass spectrometer (MC-ICP-MS). 11 B. Single-mineral supplementation analysis was performed, especially for standard samples of marine sedimentary rocks. Since boron may be enriched in specific minerals such as tourmaline, representative boron-bearing mineral grains within these samples needed to be selected under a microscope. In-situ micro-area boron isotope determination of these single minerals was performed using laser ablation multi-collector inductively coupled plasma mass spectrometry (LA-MC-ICP-MS) to obtain more refined isotopic information that may represent the initial sedimentary fluid. All the above test data, including whole-rock boron content and whole-rock δ¹⁸O₂, were then analyzed. 11 B value, δ of boron-loaded minerals 11 The B-value, the content of key inactive elements, and the corresponding geological background information of the sample, including but not limited to rock type, origin, and tectonic environment, are correlated and entered into a computer database management system to construct a structured and standardized benchmark database. This database forms the basis for subsequent quantitative comparisons in discrimination.

[0066] S2. Establish a dual-control calibration model for boron isotope fractionation. This model simulates the distribution behavior of boron between the mineral and fluid phases and the boron isotope fractionation effect under metamorphic conditions. The model uses pressure conditions and fluid activity during the metamorphic process as the core control variables. Figure 2 As shown, the conceptual model of the metamorphic system illustrates the coexistence of mineral particles and intergranular fluids at the microscale. Pressure P and fluid activity X are two key external factors controlling boron isotope fractionation.

[0067] First, for the target metamorphic minerals in the study area, namely tourmaline and columnar quartz, commonly found in the boron-rich rocks of the Larsmann Hills, a dataset of boron partition coefficients between these minerals and coexisting fluid phases, such as chloride-rich fluids, was collected or obtained through high-temperature and high-pressure experiments at a specific temperature (800°C) and under different pressure conditions, including 0.5 GPa, 1.0 GPa, 1.5 GPa, and 2.0 GPa. The partition coefficient D is defined as D=C 矿物 / C 流体 Where C is the boron concentration. Then, a nonlinear regression fitting is performed on this dataset, with pressure P as the independent variable, to establish a predictive equation for the partition coefficient. This equation describes how pressure changes systematically alter the partitioning behavior of boron between the solid and liquid phases. The fitted relationship is:

[0068]

[0069] in, , , For the fitting parameters, The pressure is the factor. This fitted relationship function can reflect the rapid adjustment of the distribution coefficient due to changes in fluid density in the initial stage of pressure increase, as well as the linear effect of the compression effect of mineral crystal structure under high pressure.

[0070] Secondly, for the same mineral-fluid system, a dataset of fractionation values ​​Δ (in units of ‰) of boron isotopes at equilibrium between the two phases was collected or obtained experimentally under different temperature conditions, including 600℃, 700℃, 800℃, and 900℃. This dataset was then fitted to establish a quantitative relationship between the fractionation value Δ and the absolute temperature T. The fitting equation describing the complex fractionation behavior is as follows:

[0071]

[0072] Where A, B, and C are the fitting parameters.

[0073] Finally, combining the distribution coefficient D and the fractionation value Δ, and introducing a key variable, namely the fluid activity factor,... Its value ranges from 0 to 1. =0 represents a completely closed system. =1 represents an extremely open system with complete fluid exchange, and the core inversion calculation equation is constructed as follows:

[0074]

[0075] This equation is the core of the calculation of the dual-control correction model, and it can convert measured values... Correct back to original value Among them, pressure P and temperature T determine the values ​​of D and Δ, while the factor This quantifies the impact of fluid activity on the fractionation path, thus achieving a dual-control simulation of the altered fractionation effect. For example... Figure 3 The figure shows two core calculation relationships of the dual-control correction model for boron isotope fractionation. The left figure shows the relationship between the distribution coefficient D and the pressure P. D changes non-linearly with increasing pressure, reflecting the response of boron's distribution behavior between minerals and fluids to pressure. The right figure shows the relationship between the fractionation value Δ and the temperature T. It shows that Δ decreases with increasing temperature, reflecting the response of the degree of boron isotope fractionation to temperature.

[0076] S3. Obtain the measured whole-rock boron isotope composition data and key inactive element content data of the metamorphic rock sample to be tested, and determine the temperature and pressure conditions of the main metamorphic stages experienced by the sample.

[0077] The boron-rich metamorphic rock sample to be tested, corresponding to tourmaline-bearing quartz gneiss, was pulverized and ground to a particle size of less than 200 mesh to prepare a uniform powder. A portion of the powder was then used in an alkaline fusion method, including co-melting with sodium hydroxide, to prepare a uniform glass target. The boron isotope composition δ¹⁸O₅ of the prepared target was analyzed using a multi-receiver inductively coupled plasma mass spectrometer (MC-ICP-MS). 11 Accurate determination of the B value. Simultaneously, the content of key inactive elements in the same sample powder, including but not limited to Ti, Nb, Ta, Zr, Hf, and Th, was determined using inductively coupled plasma mass spectrometry (ICP-MS). These data will be used for subsequent correction and discrimination.

[0078] Probe slides were prepared for the test samples, and detailed petrographic observations were conducted under a microscope. The focus was on identifying specific mineral assemblages that were in equilibrium with the target boron-rich mineral, tourmaline, and suitable for geothermobaric calculations. For example, if tourmaline was observed to be in stable contact with garnet (Grt) and biotite (Bt) without twinning, reaction edges, or other imbalance phenomena, the Grt-Bt mineral pair could be identified as the thermometer's calculation object. If further coexistence with garnet (Grt), orthopyroxene (Opx), plagioclase (Pl), and quartz (Qtz) was found, this could be used as the pressure gauge's calculation combination. Using an electron probe microanalyzer (EPMA), precise in-situ major elemental chemical composition analysis was performed on the identified coexisting mineral pairs: the Gr core and adjacent Bt, and the Gr core and adjacent Opx and Pl, ensuring that the measurement points were in a homogeneous equilibrium region.

[0079] Substituting the chemical composition of Grt and Bt, including the contents of elements such as Fe, Mg, and Mn, into the specific calculation formula of the garnet-biotite Grt-Bt, Fe-Mg exchange thermometer, the primary metamorphic temperature experienced by the sample was calculated, yielding T = 800 ± 30℃. Substituting the chemical composition of Grt, Opx, Pl, and Qtz—the Ca, Mg, and Fe contents in Grt, the Al content in Opx, and the An molecular percentage in Pl—into the formula of the garnet-orthopyroxene-plagioclase-quartz Grt-Opx-Pl-Qtz geobaric pressure gauge, the corresponding pressure value was calculated, yielding P = 1.2 ± 0.1 GPa.

[0080] Finally, the calculated temperature T≈800℃ and pressure P≈1.2GPa were used as representative temperature and pressure conditions for the main metamorphic phase of the sample and input into the subsequent correction model.

[0081] S4. Based on the determined temperature and pressure conditions, the measured whole-rock boron isotope composition data is fractionated and corrected using the dual-control calibration model, and the original boron isotope composition of the sample before metamorphism is obtained through inversion calculation. The fractionation correction and inversion calculation using the dual-control calibration model specifically includes: inputting the metamorphic temperature and pressure values ​​determined in S3 into the dual-control calibration model to calculate the corresponding boron isotope fractionation coefficient Δ and mineral-fluid partition coefficient D; and setting the fluid activity factor in the dual-control calibration model based on the petrographic analysis results of the metamorphic rock sample to be tested. The numerical value; the measured whole-rock boron isotope composition obtained in S3. The calculated Δ and D, and the set The inversion calculation equation is input and solved to output the original boron isotope composition. .

[0082] Specifically, the representative metamorphic temperature and pressure conditions determined above are used as inputs: temperature T≈800℃, i.e. 1073 K; pressure P≈1.2GPa.

[0083] Calculate the distribution coefficient D: Substitute the pressure value P = 1.2 GPa into the established distribution coefficient-pressure function relationship. In the middle. Assuming that the parameters for the tourmaline-chlorinated fluid system have been obtained through previous experimental fitting as D0=5.0, α=0.8, β=0.1, then we can calculate: D=5.0×e 0.8×1.2 +0.1×1.2≈5.0×0.383+0.12≈1.915+0.12≈2.035

[0084] Calculate the fractionation coefficient Δ: Substitute the temperature value T=1073K into the established fractionation coefficient-temperature function relationship. In the middle. Assume the parameters fitted for the same system are A=5000, B= 2.5, C=15, then we can calculate: Δ=5000 / 1073+( 2.5)×ln(1073)+15≈4.66+( 2.5)×6.978+15≈4.66 17.445 + 15 ≈ 2.215‰

[0085] Setting fluid activity factor The setting of this factor should be based on a comprehensive evaluation of geological evidence obtained from petrographic analysis of the samples.

[0086] Observations revealed that the contact boundaries between tourmaline and minerals such as garnet and biotite in the tested samples were straight and clear, without obvious dissolution or reaction rim structures, indicating that the minerals were basically in equilibrium during the main metamorphic stage. However, a small number of late-stage quartz-carbonate veinlets were found in the rock, cutting through the early foliation, suggesting the existence of localized and limited fluid activity after the metamorphic peak.

[0087] Based on the above evidence, it can be concluded that the system was close to equilibrium during the main metamorphic stage, but not absolutely closed, with limited external fluid intervention or local reequilibrium of internal fluids. Therefore, the fluid activity factor... Set to a small value between 0 and 1. Based on empirical evaluation, set... =0.2, used to characterize the system state of "mainly in near-closed equilibrium, accompanied by slight fluid activity".

[0088] Perform inversion calculations to solve for the original boron isotope composition: The obtained measured whole-rock boron isotope values ​​are then used to determine the composition. The set values ​​are +5.0‰, and the calculated Δ = 2.215‰ and D = 2.035. =0.2 and substitute it into the core inversion equation: Perform numerical calculations: =5.0+2.215×(1 0.2)×2.035(1 2.035) = 5.0 + 2.215 × 0.8 × 2.035( 1.035) = 5.0 + 1.772 × ( 0.5086) = 5.0 0.901≈4.099‰

[0089] The original whole-rock boron isotope composition before metamorphism of the boron-rich metamorphic rock sample from the Larsmann Hills was obtained by calculation using a dual-control calibration model. It is approximately +4.1‰. This value has largely eliminated the influence of boron isotope fractionation caused by metamorphism at approximately 800℃ and 1.2 GPa, providing a geochemical basis closer to its initial origin for subsequent protolith identification.

[0090] S5. The original boron isotope composition obtained from the inversion calculation is synergistically analyzed with the content data of the key inactive elements to form a comprehensive discrimination index. In S5, the synergistic analysis and formation of the comprehensive discrimination index is specifically achieved by constructing a "boron isotope-key element" synergistic discrimination system, including: using the original boron isotope composition obtained from the inversion calculation as the core isotope index; calculating at least one element ratio that is significant for discriminating the protolith type from the acquired key inactive element content data as the key element index, wherein the element ratio consists of elements with stable chemical properties and strong migration inertness during advanced metamorphism; the key inactive elements are at least two of Ti, Nb, Ta, Zr, Hf, and Th; establishing a synergistic correspondence between the core isotope index and the key element index, constructing a two-dimensional discrimination diagram with the core isotope index as the horizontal axis and the key element index as the vertical axis, and the comprehensive discrimination index corresponding to the coordinate point in the two-dimensional discrimination diagram.

[0091] Specifically, the original boron isotope composition output directly from the inversion calculation is used. =+4.1‰ was used as the core isotope index. This index carries the boron isotope signal, after correction, that best reflects the characteristics of the original source region. From the obtained data on the content of key inactive elements in the samples to be tested, elements that are chemically extremely stable during advanced metamorphism and hardly migrate with fluid activity were selected, and their ratios were calculated.

[0092] For the regional tectonic setting identification studied in this embodiment, Th and Nb elements were selected. Th is strongly enriched in felsic crust, while Nb is relatively enriched in mantle-derived rocks, and both exhibit high inactivity under granulite facies metamorphic conditions. Therefore, the ratio Th / Nb is an effective indicator for determining whether the protolith was added with ancient continental crustal material. Assuming that ICP-MS analysis shows that the sample has a Th content of 12.5 ppm and an Nb content of 8.0 ppm, the key element index is calculated as: Th / Nb = 12.5 / 8.0 = 1.56. This ratio of 1.56 is the quantified key element index.

[0093] Constructing a two-dimensional discriminant diagram and determining the comprehensive discriminant coordinate points: based on the original boron isotope composition. A two-dimensional Cartesian coordinate system is established with the x-axis as the horizontal axis and the key element ratio Th / Nb as the y-axis. This system provides a discriminant diagram. The two index values ​​X= in this embodiment are used to represent the sample under test. =+4.1‰, Y=Th / Nb=1.56 is calibrated in this diagram. This unique coordinate point (+4.1, 1.56) is the comprehensive discrimination index corresponding to this sample. This point integrates the sample's corrected boron isotope characteristics and elemental ratio characteristics unaffected by metamorphism, forming the core data pair for subsequent comparison with the protolith database. For example... Figure 4 The figure illustrates the specific process of determining the type of protolith using the "boron isotope-key element" synergistic discrimination system. The figure also shows the details of Euclidean distance discrimination: calculating the distance from the sample point to the center point of each type of protolith region.

[0094] S6. Compare the comprehensive discrimination index with the benchmark database, and determine the protolith type of the metamorphic rock sample to be tested according to the preset dual-index discrimination rule. In S6, the preset dual-index discrimination rule specifically means: in the two-dimensional discrimination diagram, a dedicated discrimination region is defined for the standard sample data of each type of protolith in the benchmark database. The boundary of the discrimination region is determined by the statistical distribution range of the core isotope index and the key element index of the corresponding type of protolith. If the coordinate point of the metamorphic rock sample to be tested falls within the discrimination region of a certain type of protolith, then the protolith of the sample is determined to belong to the corresponding category.

[0095] The determination of the protolith type of the metamorphic rock sample to be tested specifically includes: calculating the Euclidean distance from the coordinate point of the metamorphic rock sample to be tested to the center point of the various protolith discrimination regions in the benchmark database; and classifying the metamorphic rock sample to be tested as the protolith type corresponding to the discrimination region with the smallest Euclidean distance.

[0096] Specifically, before performing the discrimination, a dedicated discrimination region needs to be pre-defined for each type of protolith in the benchmark database on the constructed two-dimensional discrimination diagram. From the benchmark database, the original boron isotope composition δ¹⁸, corresponding to two key indicators of all standard samples of each type of protolith, including marine sedimentary rocks and oceanic basalt, should be used after method calibration. 11 The ratio of B to the selected key element Th / Nb.

[0097] Statistical analysis was performed on the two sets of data for each type of protolith, including but not limited to calculating the mean and standard deviation. Taking "marine sedimentary rocks" as an example, assuming that the δ0.05 of all its standard samples... 11 The B value is mainly distributed between -10‰ and +10‰, and the Th / Nb ratio is mainly distributed in the range >5. Based on this, a polygonal or elliptical region roughly covering this densely distributed data area can be drawn on the discriminant diagram as the discriminant region for "marine sedimentary rocks". Similarly, for "oceanic basalt", its δ... 11B is approximately +5‰ to +10‰, Th / Nb < 0.2, and other types such as "terrigenous clastic rocks" are delineated into their respective regions. These regions serve as the direct basis for subsequent comparisons.

[0098] Plot the coordinates (+4.1‰, 1.56) of the sample to be tested onto the diagram with the pre-defined discrimination area, and make a discrimination based on the following rules: observe which preset discrimination area the coordinate point falls into. For example, if the point happens to fall into the preset discrimination area for "island arc igneous rock", then it is preliminarily determined that its protolith may be an island arc igneous rock.

[0099] To improve the objectivity and accuracy of the discrimination, especially when the coordinate point is close to the regional boundary, a secondary confirmation is performed by calculating the Euclidean distance. The Euclidean distance from the coordinate point to the mean point of each data index corresponding to the statistical center point of each type of protolith discrimination region in the database is calculated. The calculated center point coordinates for the "Island Arc Igneous Rock" category in the benchmark database are (+5.0‰, 1.2), and the center point coordinates for the "Oceanic Crust Basalt" category are (+7.0‰, 0.1). The distances from the measured point (+4.1, 1.56) to these two centers are calculated: Distance to the center of "Island Arc Igneous Rock": Distance to the center of the "oceanic basalt": Based on the rule of "identifying the sample as the protolith type corresponding to the region with the smallest Euclidean distance", since 0.969 < 3.247, the protolith type of the boron-rich metamorphic rock sample to be tested in the Larsmann Hills was finally determined to be island arc magmatic rock.

[0100] The solution of this invention is not limited to a laboratory environment. To meet the needs of rapid, real-time field operations in geological surveys, a dedicated portable boron isotope rapid on-site analyzer has also been developed. This integrated on-site analysis equipment typically includes at least the following modules: Portable boron isotope analysis module: The core is a miniaturized, interference-resistant multi-receiver inductively coupled plasma mass spectrometer (MC-ICP-MS). This module is used to directly perform boron isotope delta analysis on rock powder pellets or simply prepared target samples on-site. 11Rapid, semi-quantitative, or quantitative determination of boron (B). Key element rapid analysis module: Integrates a portable X-ray fluorescence spectrometer. This module is used for rapid, non-destructive, or minimally destructive on-site determination of the content of key inactive elements in rock samples, including but not limited to Ti, Nb, Zr, and Th, to calculate key element ratios such as Th / Nb. Data processing and model calculation unit: Built-in high-performance embedded computer system, pre-stored with: a benchmark database and dual-control calibration models including the functional relationship between distribution coefficient and pressure, the functional relationship between fractionation coefficient and temperature, and inversion calculation equations. The established thermo-barotropic calculation program integrates common geological thermo-barotropic calculation methods. An established "boron isotope-key element" co-discrimination program and discrimination diagrams. Human-computer interaction and display module: Includes a touch screen, data interface, and result report printing output unit.

[0101] This integrated field analysis equipment eliminates the reliance on large, fixed laboratories, compressing the process of sample delivery, batch testing, offline calculations, and manual interpretation—which previously took weeks or even months—into hours in the field. This significantly improves the efficiency and responsiveness of geological surveys, especially in remote areas and extreme environments.

[0102] Based on the same inventive concept, according to another aspect of the present invention, embodiments of the present invention also provide a computer device, comprising:

[0103] At least one processor; and

[0104] The memory stores a computer program that can run on a processor. When the processor executes the program, it performs the steps of a geological survey method described above for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics.

[0105] Based on the same inventive concept, according to another aspect of the present invention, embodiments of the present invention also provide a computer-readable storage medium storing computer program instructions, which, when executed by a processor, perform the steps of the geological survey method described above for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics.

[0106] Finally, it should be noted that those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. The storage medium can be a magnetic disk, optical disk, read-only memory (ROM), or random access memory (RAM), etc. The above computer program embodiments can achieve the same or similar effects as any of the corresponding foregoing method embodiments.

[0107] Furthermore, typically, the devices and equipment disclosed in the embodiments of this invention can be various electronic terminal devices, such as mobile phones, personal digital assistants (PDAs), tablet computers (PADs), smart TVs, etc., or they can be large terminal devices, such as servers. Therefore, the scope of protection disclosed in the embodiments of this invention should not be limited to a specific type of device or equipment. The client disclosed in the embodiments of this invention can be applied to any of the above-mentioned electronic terminal devices in the form of electronic hardware, computer software, or a combination of both.

[0108] Furthermore, the method disclosed in the embodiments of the present invention can also be implemented as a computer program executed by a CPU, which may be stored in a computer-readable storage medium. When the computer program is executed by the CPU, it performs the functions defined in the method disclosed in the embodiments of the present invention.

[0109] Furthermore, the above-described method steps and system units can also be implemented using a controller and a computer-readable storage medium for storing a computer program that enables the controller to perform the functions of the above-described steps or units.

[0110] Furthermore, it should be understood that the computer-readable storage medium (e.g., memory) described herein can be volatile memory or non-volatile memory, or may include both volatile and non-volatile memory. By way of example, and not limitation, non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM), which can act as external cache memory. By way of example, and not limitation, RAM may be available in various forms, such as synchronous RAM (DRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), and direct Rambus RAM (DRRAM). The storage devices disclosed herein are intended to include, but are not limited to, these and other suitable types of memory.

[0111] The above are exemplary embodiments disclosed in this invention. However, it should be noted that various changes and modifications can be made without departing from the scope of the embodiments of this invention as defined by the claims. The functions, steps, and / or actions of the methods according to the disclosed embodiments described herein do not need to be performed in any marked order. Furthermore, although the elements disclosed in the embodiments of this invention may be described or claimed individually, they may be understood as multiple unless explicitly limited to a singular.

[0112] In this specification, the same or similar parts between the various embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments. In particular, the descriptions of the embodiments described later are relatively simple, and relevant parts can be referred to the descriptions of the foregoing embodiments.

[0113] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, characterized in that, The method includes the following steps: S1. Collect and determine the boron content, boron isotope composition and key inactive element content of various known unmetamorphosed protoliths, and establish a benchmark database including marine sedimentary rocks, oceanic crust basalt, island arc magmatic rocks and terrigenous clastic rocks; S2. Establish a dual-control calibration model for boron isotope fractionation. The dual-control calibration model is used to simulate the distribution behavior of boron between the mineral phase and the fluid phase and the boron isotope fractionation effect under metamorphic conditions. The dual-control calibration model uses the pressure conditions and fluid activity during the metamorphic process as the core control variables. S3. Obtain the measured whole-rock boron isotope composition data and key inactive element content data of the metamorphic rock sample to be tested, and determine the temperature and pressure conditions of the main metamorphic stages experienced by the sample. S4. Based on the determined temperature and pressure conditions, the measured whole-rock boron isotope composition data are fractionated and corrected using the dual-control correction model, and the original boron isotope composition of the sample before metamorphism is obtained by inversion calculation. S5. The original boron isotope composition obtained from the inversion calculation is analyzed in conjunction with the content data of the key inactive elements to form a comprehensive discrimination index; S6. Compare the comprehensive discrimination index with the benchmark database, and determine the protolith type of the metamorphic rock sample to be tested according to the preset dual-index discrimination rule; In S2, establishing the dual-control correction model includes establishing a quantitative functional relationship between the distribution coefficient and the pressure, which is determined in the following way: Based on high-temperature and high-pressure experiments simulating geological metamorphic processes, data on the distribution coefficients of boron between the target metamorphic minerals and coexisting fluid phases under different pressure conditions were obtained. A functional relationship between the distribution coefficient D and the pressure P was derived through nonlinear regression fitting. Specifically, the functional relationship is as follows: ； in, , , For the fitting parameters, For pressure; The establishment of the dual-control correction model also includes establishing a quantitative functional relationship between the boron isotope fractionation coefficient and temperature, and constructing an inversion calculation equation, wherein: the quantitative functional relationship between the fractionation coefficient and temperature is determined in the following way: Based on high-temperature and high-pressure experiments simulating geological metamorphic processes, boron isotope equilibrium fractionation values ​​between target metamorphic minerals and fluid phases under different temperature conditions were obtained. A nonlinear regression fitting was used to derive a functional relationship between the fractionation value Δ and temperature T. Specifically, the functional relationship is as follows: ； Where A, B, and C are fitting parameters, and T is the absolute temperature; The inversion calculation equation is: ； in, It has the original boron isotopic composition. To measure the whole-rock boron isotope composition, As a fluid activity factor, For allocation coefficients, The fractionation value; In S5, the collaborative analysis and formation of comprehensive discrimination indicators are specifically achieved by constructing a "boron isotope-key element" collaborative discrimination system, including: The original boron isotope composition obtained from the inversion calculation is used as the core isotope index; From the obtained key inactive element content data, at least one element ratio that is significant in distinguishing the original rock type is calculated as a key element index. The element ratio is composed of elements that are chemically stable and have strong migration inertness during advanced metamorphism. The key inactive element is at least two of Ti, Nb, Ta, Zr, Hf, and Th. A synergistic correspondence is established between the core isotope index and the key element index. A two-dimensional discriminant diagram is constructed with the core isotope index as the horizontal axis and the key element index as the vertical axis. The comprehensive discriminant index corresponds to the coordinate point in the two-dimensional discriminant diagram.

2. The geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics according to claim 1, characterized in that, In S1, establishing the benchmark database specifically includes: Marine sedimentary rocks, oceanic basalts, island arc igneous rocks, and terrigenous clastic rocks with clear geological backgrounds and no significant metamorphic alteration were selected as standard samples. The whole-rock boron content and whole-rock boron isotope δ¹⁸O of the standard sample were determined. 11 Composition B, and for the aforementioned marine sedimentary rock standard samples, the single-mineral boron isotope δ¹⁸ of representative boron-bearing minerals in the samples were additionally measured. 11 B value; The measured boron content data and whole-rock boron isotope δ¹⁸O will be used to... 11 B value, boron isotope δ of single mineral 11 The B-value is correlated with the content data of key inactive elements corresponding to the standard sample and geological background information to construct a standardized benchmark database.

3. The geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics according to claim 1, characterized in that, In S3, the measured whole-rock boron isotope composition data are obtained by using a multi-receiver inductively coupled plasma mass spectrometer to measure the powder pellets or fused glass targets prepared from the metamorphic rock sample to be tested. The determination of the temperature and pressure conditions of the main metamorphic stages experienced by the sample is specifically based on the petrographic observation results of the metamorphic rock sample to be tested, and quantitative calculation is performed using a symbiotic mineral geothermobarometer; wherein, the symbiotic minerals include typical symbiotic assemblages of target boron-rich minerals, and the temperature and pressure conditions are calculated based on the chemical composition of the characteristic minerals in the symbiotic assemblage.

4. A geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, as described in claim 3, is characterized in that... The temperature and pressure conditions for determining the main metamorphic stages experienced by the sample specifically include: Based on petrographic observations of the metamorphic rock samples to be tested, specific mineral pairs associated with the target boron-rich minerals are identified; wherein, the specific mineral pairs include at least one of the following: a garnet-biotite mineral pair for temperature calculation, or a garnet-orthopyroxene-plagioclase-quartz mineral assemblage for pressure calculation; the chemical composition of the specific mineral pairs is determined by electron probe microanalysis. Input the chemical composition data into the corresponding geothermal barometer calculation formula to calculate the temperature and pressure values ​​respectively.

5. A geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, as described in claim 4, is characterized in that... The dual-control correction model is invoked to perform fractionation correction and inversion calculations, specifically including: The values ​​of metamorphic temperature and metamorphic pressure determined in S3 are input into the dual-control correction model to calculate the corresponding boron isotope fractionation coefficient Δ and mineral-fluid partition coefficient D, respectively. Based on the petrographic analysis results of the metamorphic rock sample to be tested, the fluid activity factor in the dual-control calibration model was set. The value; The measured whole-rock boron isotope composition obtained in S3 The calculated Δ and D, and the set The inversion calculation equation is input and solved to output the original boron isotope composition. .

6. A geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, as described in claim 1, is characterized in that... In step S6, the step of judging according to the preset dual-index rule specifically includes: In the two-dimensional discrimination diagram, a specific discrimination region is defined for the standard sample data of each type of protolith in the benchmark database. The boundary of the discrimination region is determined by the statistical distribution range of the core isotope index and the key element index of the corresponding type of protolith. If the coordinates of the metamorphic rock sample to be tested fall within the discrimination area of ​​a certain type of protolith, then the protolith of the sample is determined to belong to the corresponding category.

7. A geological survey method for reconstructing the protolith type of metamorphic rocks using boron isotope fractionation characteristics, as described in claim 6, is characterized in that... Determining the protolith type of the metamorphic rock sample to be tested specifically includes: Calculate the Euclidean distance from the coordinates of the metamorphic rock sample to be tested to the center point of each type of protolith discrimination region in the benchmark database; The metamorphic rock sample to be tested is identified as the protolith type corresponding to the discrimination region with the smallest Euclidean distance.

Images