Gold ore magma identification method based on high-dimensional element ratio phase space differentiation

Abstract

The invention discloses a gold ore magma identification method based on high-dimensional element ratio phase space differentiation, and solves the technical problem that gold ore magma cannot be accurately identified in the prior art. According to the method, a four-dimensional feature vector is constructed according to the ratio of Ba / Zr, Nb / Y, U / Yb, V / Y and the like, the four-dimensional feature vector is projected to a ratio phase space sub-domain, and primary judgment, entropy state judgment and source region inverse deduction of the metallogenic potential magma are achieved through distribution trend, trajectory disturbance and multi-factor clustering behaviors. A multi-domain different-dimensional mapping and empowerment fusion model is introduced, initial ratio disturbance, differentiation evolution trajectory and component source coupling effect are integrated, and a set of high-order metallogenic property discrimination system for prospecting prediction is formed. Through the scheme, the method has the advantages of being simple in logic, accurate, reliable and the like, and has very high practical value and popularization value in the technical field of mineral resource prediction and magma earth discrimination.

Description

Technical Field

[0001] The present invention relates to the technical field of mineral resource prediction and magmatic earth discrimination, and in particular to a gold ore magma identification method based on the differentiation of high-dimensional element ratio phase space. Background Art

[0002] The formation mechanism of gold ore and prospecting prediction are among the frontiers of geoscience research. Research shows that the magmatic-hydrothermal process is crucial for gold mineralization, and the characteristics of the magmatic source area, melting conditions, and subsequent evolution control the enrichment and migration of gold. Traditional studies on the gold mineralization potential mostly rely on local element content, characteristics of post-formed fluids, or retrospective analysis after the discovery of ore bodies, lacking systematic identification means based on the primary information of the magmatic system, such as the public technology with the publication number CN117993578A and the name of a gold ore target prediction method, system, computer device, and medium; another example is the public technology with the publication number CN117630147A and the name of a method for locating the burial depth of Carlin-type gold ore bodies using rare earth elements and carbon isotopes; and another example is the article "The Composition of Gold and Other Elements in Mesozoic Magmatic Rocks in Jiaodong and Its Implications for Gold Mineralization".

[0003] With the development of geochemical methods, the ratios of trace elements, especially incompatible elements (such as Ba, Nb, U, Zr, Y, Yb, etc.), can serve as sensitive indicators of partial melting, fluid metasomatism, and fractional evolution in the magmatic source area. However, existing methods are often limited to single-element indicators or local ratio analysis (such as the article "The Composition of Gold and Other Elements in Mesozoic Magmatic Rocks in Jiaodong and Its Implications for Gold Mineralization"), lacking a theoretical system and technical path to systematically combine multiple groups of element ratios into a high-dimensional feature space and then discriminate the mineralization potential of the magmatic system as a whole. Especially in the gold mineralization system, the weak differentiation signals, source area enrichment status, and water content changes within the magmatic system are often difficult to effectively identify under conventional two-dimensional single-factor discrimination.

[0004] Therefore, there is an urgent need to construct a discrimination method centered on the differentiation of high-dimensional element ratio phase space that can accurately reveal the magmatic mineralization potential. Summary of the Invention

[0005] In view of the above problems, the purpose of the present invention is to provide a gold ore magma identification method based on the differentiation of high-dimensional element ratio phase space. The technical solution adopted by the present invention is as follows: A gold ore magma identification method based on the differentiation of high-dimensional element ratio phase space, comprising the following steps: Step S1, constructing an element abundance spectrum for the collected magmatic samples and performing sample spectrum screening; Step S2, generating an incompatible element ratio matrix using the screened element abundance spectrum and performing eigenvector quantization processing; Step S3: Construct a high-dimensional feature interface based on the feature tensor, preliminarily screen the metallogenic potential, and mark the preliminary discrimination label. Step S4: Qualitatively discriminate the magma water entropy flow of the magma samples with preliminary discrimination labels to obtain the entropy flow disturbance characteristics. Step S5: According to the results of the qualitative discrimination of the magma water entropy flow, perform source area enrichment factor mapping and inverse deduction of fluid enrichment degree to obtain the source area properties. Step S6: Perform multi-domain discrimination fusion and comprehensive weighted evaluation of metallogenic potential based on the preliminary discrimination label, entropy flow disturbance characteristics, and source area properties.

[0006] Furthermore, in the said Step S1, constructing an element abundance pedigree for the collected magma samples and performing sample pedigree screening includes the following steps: Collect the magma samples and measure the element contents of Ba, Zr, Nb, Y, U, Yb, V, and SiO2 to form an element abundance pedigree matrix C ; Screen and eliminate the magma samples with significant secondary transformation according to the alteration degree, mineral fidelity, and element ratio anomaly of the magma samples.

[0007] Furthermore, in the said Step S2, generating an incompatible element ratio matrix using the screened element abundance pedigree and performing feature tensor processing includes the following steps: Calculate the ratio of any group of magma samples, where: The expression of the Ba / Zr ratio is: ; The expression of the Nb / Y ratio is: ; The expression of the U / Yb ratio is: ; The expression of the V / Y ratio is: ; where represents the content of the Ba element in the magma sample; represents the content of the Zr element in the magma sample; represents the content of the Nb element in the magma sample; represents the content of the Y element in the magma sample; represents the content of the U element in the magma sample; represents the content of the V element in the magma sample; represents the content of the Yb element in the magma sample; Taking { , , , Construct a sample high-dimensional ratio feature tensor using it as a basic feature component .

[0008] Furthermore, in step S3, construct a high-dimensional feature interface based on the feature tensor, preliminarily screen the metallogenic potential, and mark the preliminary discrimination label, including the following steps: Preset the classification discrimination threshold R 1. R 2. R 3. If > R 1, > R 2, > R 3; then the magma sample is a gold metallogenic potential magma; otherwise, it is a non-mineral magma; Map the magma samples collected in step S1 to the Ba / Zr–SiO2 two-dimensional projection plane, and observe the differentiation and aggregation of gold metallogenic potential magma and non-mineral magma; Mark the preliminary discrimination label of the magma sample in the high-dimensional space according to the differentiation and aggregation of gold metallogenic potential magma and non-mineral magma.

[0009] Furthermore, in step S4, conduct a qualitative discrimination of the magma water entropy flow for the magma samples with preliminary discrimination labels to obtain the entropy flow perturbation characteristics, including the following steps: Plot the magma samples with preliminary discrimination labels in step S3 in the V / Y–Ba / Zr projection subspace and determine their distribution patterns; Analyze the water content trend and entropy flow perturbation characteristics to obtain the entropy flow perturbation characteristics and obtain the magma samples with metallogenic potential.

[0010] Furthermore, in step S5, according to the results of the qualitative discrimination of the magma water entropy flow, conduct a mapping of the source area enrichment factor and an inverse deduction of the fluid enrichment degree to obtain the source area properties, including the following steps: Plot the magma samples that have been preliminarily screened and determined to be gold metallogenic potential magma in steps S3 and S4 in the U / Yb–Nb / Y two-dimensional subspace, construct a source area enrichment mapping, and infer the source area composition type and evolution characteristics based on the sample distribution trend.

[0011] Furthermore, in step S6, conduct a multi-domain discrimination fusion and a comprehensive weighted evaluation of the metallogenic potential according to the preliminary discrimination label, entropy flow perturbation characteristics, and source area properties, including the following steps: Construct a metallogenic potential function S according to the preliminary discrimination label, entropy flow perturbation characteristics, and source area properties, and its expression is: Wherein, represents the Ba / Zr ratio , Nb / Y ratio and U / Yb ratio The metallogenic potential information characterized by the constructed high-dimensional tensor features; Indicates the magma water content trend and entropy flow perturbation information revealed based on the V / Y–Ba / Zr projection relationship; Indicates the U / Yb ratio and Nb / Y ratio The source area enrichment property and genetic type information constructed; Indicates the preliminary metallogenic potential scoring function in the high-dimensional tensor space; Indicates the magma water content trend and entropy flow perturbation mapping scoring function; Indicates the source area enrichment factor mapping scoring function; Use the ore potential function S to divide any magma sample to form the spatial distribution of the gold metallogenic potential of the magma system.

[0012] Compared with the prior art, the present invention has the following beneficial effects: (1) The present invention constructs an elemental abundance spectrum for the collected magma samples, generates an incompatible element ratio matrix using the screened elemental abundance spectrum, and performs eigen-tensorization processing, which systematically constructs the element ratios into high-dimensional eigen-tensors, comprehensively considering the source area enrichment, partial melting degree, and fluid evolution state, ensuring its accurate and reliable identification.

[0013] (2) The present invention plots the samples in the V / Y–Ba / Zr projection subspace and the U / Yb–Nb / Y two-dimensional subspace, introduces phase space differentiation analysis, and judges the metallogenic affinity of the magma system through characteristics such as the aggregation and trajectory trend of the samples in the high-dimensional space.

[0014] (3) The present invention uses the entropy flow perturbation theory to analyze the magma water content change trajectory from the V / Y–Ba / Zr projection plane, further improving the accuracy of metallogenic potential discrimination.

[0015] (4) The present invention constructs an ore potential function and performs multi-domain discrimination fusion and comprehensive weighted evaluation of metallogenic potential using preliminary discrimination labels, entropy flow perturbation characteristics, and source area properties, which is suitable for batch data processing and regional prospecting prediction.

[0016] In summary, the present invention has the advantages of simple logic, accuracy, and reliability, and has high practical value and promotion value in the technical fields of mineral resource prediction and magma geochemistry discrimination. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] To more clearly illustrate the technical solutions of the embodiments of the present invention, the following will briefly introduce the accompanying drawings required for the embodiments. It should be understood that the following drawings only show some embodiments of the present invention, and thus should not be regarded as limiting the scope of protection. For those skilled in the art, without creative efforts, other related drawings can also be obtained based on these drawings.

[0018] Figure 1 It is the logic flow chart of the present invention.

[0019] Figure 2 It is the discrimination scatter plot of gold - forming magma and ore - free magma of the present invention on the Ba / Zr–SiO2 plane.

[0020] Figure 3 It is the schematic diagram revealing the characteristics of entropy flow disturbance by the V / Y–Ba / Zr relationship of the present invention.

[0021] Figure 4 It is the schematic diagram for discriminating the enrichment degree of the source area by using U / Yb and Nb / Y of the present invention. Detailed implementation manners

[0022] To make the objectives, technical solutions and advantages of the present application clearer, the following further describes the present invention with reference to the accompanying drawings and embodiments. The implementation manners of the present invention include but are not limited to the following embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts belong to the scope of protection of the present application.

[0023] In this embodiment, the term "and / or" only describes the associated relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist simultaneously, and B exists alone.

[0024] The terms "first" and "second" etc. in the description and claims of this embodiment are used to distinguish different objects, rather than to describe a specific order of the objects. For example, the first target object and the second target object etc. are used to distinguish different target objects, rather than to describe the specific order of the target objects.

[0025] In the embodiments of the present application, words such as "exemplary" or "for example" are used to represent examples, illustrations or explanations. Any embodiment or design solution described as "exemplary" or "for example" in the embodiments of the present application should not be construed as being more preferred or having more advantages than other embodiments or design solutions. Exactly, using words such as "exemplary" or "for example" aims to present relevant concepts in a specific manner.

[0026] In the description of the embodiments of the present application, unless otherwise specified, "a plurality" means two or more. For example, a plurality of processing units means two or more processing units; a plurality of systems means two or more systems.

[0027] As Figures 1 to 4 shown, this embodiment provides a gold ore magma identification method based on the differentiation of high-dimensional element ratio phase space, which includes the following steps: The first step is to construct an element abundance spectrum for the collected magma samples and perform sample spectrum screening. Specifically: First, collect fresh or slightly altered igneous rock samples, measure the contents of elements such as Ba, Zr, Nb, Y, U, Yb, V, and SiO2 to form an element abundance spectrum matrix C. According to the sample alteration degree, mineral fidelity, and element ratio anomalies, eliminate samples with significant secondary transformation to ensure that subsequent analysis is based on the primary magma system. This embodiment takes the magma samples in a certain area as an example, and forms the basic data set shown in Table 1.

[0028] Table 1 Basic data set Sample number Category SiO2 (wt%) Ba / Zr Nb / Y U / Yb V / Y 1 Gold ore magma 58.3 3.6 1.3 1.9 0.82 2 Gold ore magma 56.7 3.2 1.1 1.7 0.79 … … … … … … … 51 Ore-free magma 59.1 1.9 0.7 0.8 1.5 52 Ore-free magma 60.5 1.6 0.6 0.9 1.7 … … … … … … … The second step is the generation of incompatible element ratio matrices and feature tensorization. Specifically: Calculate the ratio of any group of magma samples, where: The expression of the Ba / Zr ratio is: ; The expression of the Nb / Y ratio is: ; The expression of the U / Yb ratio is: ; The expression of the V / Y ratio is: ; Among them, represents the content of the Ba element in the magma sample; represents the content of the Zr element in the magma sample; represents the content of the Nb element in the magma sample; represents the content of the Y element in the magma sample; represents the content of the U element in the magma sample; represents the content of the V element in the magma sample; represents the content of the Yb element in the magma sample.

[0029] Taking { , , , Construct a sample high-dimensional ratio feature tensor as the basic feature component .

[0030] Thirdly, construct a high-dimensional feature interface based on the feature tensor, preliminarily screen the metallogenic potential, and mark the preliminary discrimination labels.

[0031] Here, set the classification and discrimination conditions. If it satisfies > R 1, > R 2, > R 3, then the sample is determined to be a magmatic rock with gold metallogenic potential, otherwise it is a magmatic rock without ore. R 1、 R 2、 R 3 are the set discrimination thresholds.

[0032] Map the samples to the Ba / Zr–SiO2 two-dimensional projection plane, draw a preliminary discrimination scatter plot, observe the differentiation and aggregation of ore-forming and non-ore samples, and record the preliminary discrimination labels of the samples in the high-dimensional space.

[0033] Fourthly, qualitatively discriminate the magmatic water entropy flow, which specifically includes the following steps: Plot the magmatic samples with preliminary discrimination labels in the V / Y–Ba / Zr projection subspace. Determine the distribution pattern. Among them, the nearly horizontal trend may represent high water content in the magma, late degassing, which is beneficial to gold mineralization; the nearly vertical trend may represent water-poor magma, early degassing, which is not conducive to gold mineralization.

[0034] Here, use the V / Y–Ba / Zr relationship to reveal the entropy flow perturbation characteristics as an important auxiliary index for the metallogenic potential.

[0035] Fifthly, map the source area enrichment factor and inversely deduce the fluid enrichment degree, which specifically includes the following steps: (501) Plot the magmatic samples determined to have metallogenic potential in the third and fourth steps in the U / Yb–Nb / Y two-dimensional subspace to construct a source area enrichment map.

[0036] (502) Determine the nature of the source area. Among them, samples with high U / Yb and high Nb / Y may indicate an enriched lithospheric mantle source; samples with low U / Yb and low Nb / Y may indicate a depleted asthenosphere source. For example, when determining the nature of the source area, based on the ratio coordinates plotted in the U / Yb–Nb / Y projection space, the following discrimination conditions are set: when the U / Yb ratio of the sample is not less than 1.5 and the Nb / Y ratio is not less than 1.0, it indicates that its source area may be an enriched lithospheric mantle and has a gold metallogenic geochemical origin; when the U / Yb ratio of the sample is not higher than 0.8 and the Nb / Y ratio is not higher than 0.6, it indicates that its source area may be a depleted asthenosphere or ordinary fore-arc mantle, and the metallogenic potential is weak; for samples between the above intervals, the ratio information in the third and fourth steps can be combined for further comprehensive evaluation.

[0037] Step 6: Conduct multi-domain discrimination fusion and comprehensive weighted assessment of metallogenic potential based on the preliminary discrimination label, entropy flow perturbation characteristics, and source area nature.

[0038] Here, construct a metallogenic potential function S based on the preliminary discrimination label, entropy flow perturbation characteristics, and source area nature. Its expression is:

[0039] Among them, represents the Ba / Zr ratio 、the Nb / Y ratio and the U / Yb ratio characterize the metallogenic potential information represented by the high-dimensional tensor features constructed; represents the magma water content trend and entropy flow perturbation information revealed based on the V / Y–Ba / Zr projection relationship; represents the U / Yb ratio and the Nb / Y ratio construct the source area enrichment nature and genetic type information.

[0040] Table 2 is the specific correspondence table of the function

[0041] Furthermore, use the metallogenic potential function S to divide any magma sample to form the spatial distribution of the gold metallogenic potential in the magma system.

[0042] Such as Figures 2 to 4 shown, Figure 2 shows the distribution of gold ore magma and non-ore magma samples on the Ba / Zr–SiO2 projection plane. From Figure 2It can be seen that the two types of samples show an obvious trend of zonal aggregation in the two-dimensional space, indicating that the constructed high-dimensional ratio feature tensor has good differentiation ability, can effectively support the preliminary identification of metallogenic potential, and reflects the discriminative visibility of the method of the present invention under low-dimensional projection. From Figure 3 It can be seen that the distribution of magmatic samples with metallogenic potential generally tends to be horizontal, while the distribution of ore-free samples is mostly vertical, reflecting the significant differences in the water content and degassing stage of magma. This trend reveals the coupling relationship between the entropy flow perturbation characteristics and metallogenic affinity, and verifies the effectiveness and interpretability of the present invention in the identification of magmatic evolution stages. Figure 4 shows the distribution pattern of samples in the two-dimensional space constructed by the ratios of U / Yb and Nb / Y. From Figure 4 It can be seen that the samples in the enriched source area and the depleted source area form two obvious trend zones, indicating that the method of the present invention can effectively infer the type of source area composition of magma and identify the lithospheric mantle or asthenospheric source based on this. This further supports the advantages of the present invention in genetic tracing and source area inference, and provides a reliable geochemical genetic judgment basis for the gold metallogenic system. Thus, Figures 2 to 4 From the three aspects of preliminary differentiation, magmatic entropy flow characteristics and source area genetic discrimination, the effectiveness, resolution and interpretability of this embodiment in the identification of gold metallogenic potential are verified, and it has good practical applicability and regional promotion prospects.

[0043] The above embodiments are only the preferred embodiments of the present invention, and do not limit the protection scope of the present invention. Any changes made by using the design principle of the present invention and non-creative labor on this basis shall fall within the protection scope of the present invention.

Claims

1. A gold ore magma identification method based on the differentiation of high-dimensional element ratio phase space, characterized in that, It includes the following steps: Step S1: Construct an elemental abundance pedigree for the collected magma samples and conduct sample pedigree screening; Step S2: Generate an incompatible element ratio matrix using the screened elemental abundance pedigree and perform eigen - tensor processing; Step S3: Construct a high - dimensional eigen - interface based on the eigen - tensor, preliminarily screen the metallogenic potential, and mark the preliminary discrimination labels; Step S4: Conduct a qualitative discrimination of the magma water - containing entropy flow for the magma samples with preliminary discrimination labels to obtain entropy flow perturbation characteristics; Step S5: According to the results of the qualitative discrimination of the magma water - containing entropy flow, conduct source - area enrichment factor mapping and inverse deduction of fluid enrichment degree to obtain the source - area properties; Step S6: Conduct multi - domain discrimination fusion and comprehensive weighted evaluation of the metallogenic potential based on the preliminary discrimination labels, entropy flow perturbation characteristics, and source - area properties.

2. The gold ore magma identification method based on the high-dimensional element ratio phase space differentiation according to claim 1, wherein In the said Step S1, constructing an elemental abundance pedigree for the collected magma samples and conducting sample pedigree screening includes the following steps: Collect the magma samples and determine the element contents of Ba, Zr, Nb, Y, U, Yb, V, and SiO2 to form an elemental abundance spectrum matrix C ; Screen and eliminate the magma samples with significant secondary transformation according to the alteration degree, mineral fidelity, and elemental ratio anomaly of the magma samples.

3. The gold ore magma identification method based on the differentiation of high-dimensional element ratio phase space according to claim 2, characterized in that, In the said Step S2, generating an incompatible element ratio matrix using the screened elemental abundance pedigree and performing eigen - tensor processing includes the following steps: Calculate the ratio of any group of magma samples, where: Ba / Zr ratio is expressed as: ; Nb / Y ratio The expression is: ; U / Yb ratio The expression is as follows: ; V / Y ratio The expression is: ; Among them, represents the content of Ba element in the magma sample; represents the content of Zr element in the magma sample; represents the content of Nb element in the magma sample; represents the content of Y element in the magma sample; represents the content of U element in the magma sample; represents the content of V element in the magma sample; represents the content of Yb element in the magma sample; Taking { , , , } as the basic feature components, construct the sample high-dimensional ratio feature tensor .

4. A gold ore magma identification method based on the differentiation of high-dimensional element ratio phase space according to claim 3, characterized in that In the said Step S3, constructing a high - dimensional eigen - interface based on the eigen - tensor, preliminarily screening the metallogenic potential, and marking the preliminary discrimination labels includes the following steps: Preset classification discrimination threshold R 1、 R 2、 R 3. If > R 1, > R 2, > R 3; then the magma sample is a magma with gold metallogenic potential; otherwise, it is a magma without ore. Map the magma samples collected in Step S1 to the Ba / Zr–SiO2 two - dimensional projection plane, and observe the differentiation and aggregation of gold - metallogenic potential magma and non - ore magma; According to the differentiation and aggregation of gold - metallogenic potential magma and non - ore magma, mark the preliminary discrimination labels of the magma samples in the high - dimensional space.

5. A gold ore magma identification method based on the phase space differentiation of high-dimensional element ratio according to claim 4, characterized in that, In the said Step S4, conducting a qualitative discrimination of the magma water - containing entropy flow for the magma samples with preliminary discrimination labels to obtain entropy flow perturbation characteristics, which includes the following steps: Plot the magma samples with preliminary discrimination labels in Step S3 on the V / Y–Ba / Zr projection subspace and determine their distribution patterns; Analyze the water - containing trend and entropy flow perturbation characteristics to obtain entropy flow perturbation characteristics and obtain the magma samples with metallogenic potential.

6. The gold ore magma identification method based on the differentiation of high-dimensional element ratio phase space according to claim 5, characterized in that In the said Step S5, according to the results of the qualitative discrimination of the magma water - containing entropy flow, conducting source - area enrichment factor mapping and inverse deduction of fluid enrichment degree to obtain the source - area properties includes the following steps: Plot the magma samples preliminarily screened and determined to be gold - metallogenic potential magma in Steps S3 and S4 on the U / Yb–Nb / Y two - dimensional subspace, construct a source - area enrichment mapping, and infer the source - area composition type and evolution characteristics based on the sample distribution trend.

7. The gold ore magma identification method based on the phase space differentiation of high-dimensional element ratios according to claim 6, characterized in that In the said Step S6, conducting multi - domain discrimination fusion and comprehensive weighted evaluation of the metallogenic potential based on the preliminary discrimination labels, entropy flow perturbation characteristics, and source - area properties includes the following steps: Construct a metallogenic potential function S according to the preliminary discrimination labels, entropy flow perturbation characteristics, and source - area properties, and its expression is: Among them, represents the Ba / Zr ratio , the Nb / Y ratio and the U / Yb ratio of the metallogenic potential information characterized by the constructed high-dimensional tensor features; represents the magma water content trend and entropy flow perturbation information revealed based on the V / Y–Ba / Zr projection relationship; represents the U / Yb ratio and the Nb / Y ratio of the source area enrichment property and genetic type information constructed; represents the preliminary metallogenic potential scoring function in the high-dimensional tensor space; represents the magma water content trend and entropy flow perturbation mapping scoring function; represents the source area enrichment factor mapping scoring function; Use the metallogenic potential function S to divide any magma sample to form the spatial distribution of the gold - metallogenic potential of the magma system.

Images