Nb-ta metal ore mineralization evaluation system based on multi-source information analysis

Abstract

The present application relates to the technical field of multi-source information analysis, in particular to a niobium-tantalum metal mine mineralization evaluation system based on multi-source information analysis, which comprises a geological structure correlation module, a mineral distribution interpretation module, a mineralization environment adaptation module, a spatial feature analysis module and a potential area determination module. In the present application, the multi-source data analysis and fusion method is introduced to effectively improve the evaluation accuracy and comprehensive ability of resource exploration. The analysis after information fusion combines the correlation and spatial pattern of different data types to ensure that the mineralization potential evaluation is more accurate. The limitations of single data source and manual analysis are overcome to improve the identification accuracy and reliability of the mineralization regularity. The potential area can be more accurately identified. The adaptability of the mineralization environment is combined to optimize the resource potential analysis, significantly improve the efficiency and accuracy of the mineralization area determination, avoid the problems of information isolation and analysis deviation in the traditional method, and enhance the scientificity and practical application value of the exploration process.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of multi-source information analysis, in particular to a niobium-tantalum metal mine mineralization evaluation system based on multi-source information analysis. BACKGROUND

[0002] The technical field of multi-source information analysis belongs to the category of information processing and comprehensive analysis, mainly studies how to extract, integrate and correlate information from different sources, different structures and different types of data, and identifies, judges and reasons the target object through unified modeling and analysis means. The core of this technical field includes multi-source data acquisition, data preprocessing, feature extraction, data fusion and correlation analysis, etc. It covers the integrated use of multi-type information such as image data, geological data, remote sensing data and sensor data, and is widely used in resource exploration, environmental monitoring, safety warning and other scenarios. This technical field emphasizes the diversity of data sources, the complementarity between data, and the depth of analysis after information fusion. The overall technical path presents the structured integration of multi-dimensional data, the unified extraction of cross-domain features, and the improvement of comprehensive judgment ability.

[0003] Among them, the traditional niobium-tantalum and other rare metal mine mineralization evaluation system refers to the process of mineral exploration. Through qualitative judgment and empirical analysis of geological exploration data, relying on single source data such as drilling profile information, mineral assemblage distribution rule and rock geochemical index, combined with artificial recognition and two-dimensional map interpretation method, the comprehensive evaluation of mineralization geological conditions is carried out. In the traditional way, geological exploration sample analysis, artificial interpretation of remote sensing images and rock and mineral experimental detection are used to identify ore body position and predict resource potential. There are problems such as single data, low information utilization rate and limited identification accuracy of mineralization regularity.

[0004] The existing technology relies on single data source and empirical judgment, and relies too much on traditional data such as drilling profile, mineral assemblage rule and geochemical index. The information acquisition in the analysis process is limited, and it is difficult to fully evaluate the potential area in complex geological environment. The traditional method relies too much on artificial interpretation of remote sensing images and sample analysis, resulting in low efficiency of data integration and limited analysis accuracy. It cannot fully tap the complementarity of multi-source information, and thus affects the accurate identification of mineralization regularity. In the process of geological exploration, too much reliance on empirical judgment and two-dimensional map interpretation method, ignoring the depth correlation and cross-domain analysis between data, making it difficult to meet the multi-dimensional information demand in complex geological environment, thus limiting the accuracy and reliability of potential area evaluation, leading to misjudgment and resource waste. SUMMARY

[0005] In order to solve the technical problems existing in the prior art, the present application provides a niobium-tantalum metal mine mineralization evaluation system based on multi-source information analysis. The technical solution is as follows:

[0006] In one aspect, a system for evaluating mineralization of niobium-tantalum metal ore based on multi-source information analysis is provided, and the system comprises:

[0007] The geological structure correlation module obtains geological structure information, including stratum distribution, fault position and fold shape, analyzes the influence degree of tectonic stress field on mineralization, and generates a tectonic stress influence index;

[0008] The mineral distribution interpretation module extracts mineral combination rules and chemical composition distribution characteristics based on the tectonic stress influence index, analyzes the corresponding relationship between mineral enrichment area and tectonic stress field, and obtains a mineral enrichment degree value;

[0009] The mineralization environment adaptation module extracts stratum lithology and structural characteristics in the mineralization environment according to the mineral enrichment degree value, analyzes the relevance of lithology combination and mineralization potential, and generates an environment adaptation index;

[0010] The spatial feature analysis module calls the environment adaptation index, extracts the spatial distribution pattern and boundary characteristics of the mineralization area, analyzes the relationship between spatial continuity and mineralization potential, and generates a spatial distribution optimization set.

[0011] As a further scheme of the present application, the tectonic stress influence index includes stress distribution uniformity, tectonic disturbance intensity and stress action directionality, the mineral enrichment degree value includes enrichment intensity level, mineral concentration coefficient and spatial aggregation distribution degree, the environment adaptation index includes lithology adaptability score, tectonic coordination degree and mineralization environment consistency, and the spatial distribution optimization set includes distribution continuity score, boundary integrity index and regional morphological regularity.

[0012] As a further scheme of the present application, the geological structure correlation module comprises:

[0013] The stratum distribution extraction submodule obtains geological structure information, extracts stratum thickness variation and lithology interface characteristics, analyzes the influence degree of stratum relief on mineralization, and generates fault activity intensity;

[0014] The fold shape analysis submodule extracts axial plane inclination and wavelength change value in fold shape according to the fault activity intensity, analyzes the action degree of fold deformation on mineral migration and enrichment, and generates fold deformation degree;

[0015] The stress field evaluation submodule extracts tectonic stress field distribution information based on the fold deformation degree, analyzes the influence range of tectonic stress concentration on mineralization, and generates a tectonic stress influence index.

[0016] As a further scheme of the present application, the mineral distribution interpretation module comprises:

[0017] The mineral assemblage extraction submodule acquires mineral assemblage sample data based on the tectonic stress influence index, decomposes the characteristics of mineral particles in the sample, determines the content proportion of mineral types, extracts the chemical composition values of the mineral particles, compares the distribution of adjacent particles in the differential mineral assemblage, identifies the paragenetic mineral assemblage block, and generates mineral type diversity;

[0018] The chemical composition analysis submodule screens key elements in the differential mineral assemblage according to the mineral type diversity, calculates the stable fluctuation interval of the elements, extracts the stable element proportion in the sample, extracts the variable element content variation section, discriminates the differential element assemblage relationship, calibrates the influence degree of the differential assemblage on the sample distribution, and generates chemical composition stability.

[0019] The enrichment zone distribution scoring submodule calls the chemical composition stability, segments the mineral enrichment zone sample unit, determines the key mineral content density in the enrichment unit, superimposes the unit and the corresponding range of the tectonic stress block, identifies the high-density enrichment block boundary, quantifies the demarcation point set of the enrichment block and the low-density area, and obtains the mineral enrichment degree value.

[0020] As a further scheme of the present application, the mineral enrichment degree value uses the formula:

[0021] ;

[0022] Wherein, represents the mineral enrichment degree value, represents the content density of the i th mineral sample, represents the distribution density of the i th mineral sample, represents the unit volume mineral quality of the i th mineral sample, represents the average value of the unit volume mineral quality of all units, represents the total number of mineral samples.

[0023] As a further scheme of the present application, the ore-forming environment adaptation module comprises:

[0024] The lithology assemblage extraction submodule extracts the stratigraphic lithology assemblage type in the ore-forming environment according to the mineral enrichment degree value, analyzes the support degree of the lithology assemblage to mineral enrichment, and generates a lithology assemblage matching degree.

[0025] The tectonic feature matching submodule calls the lithology assemblage matching degree, extracts the fault strike and fold shape in the tectonic feature, analyzes the contribution degree of the tectonic feature to mineralization, and generates a tectonic feature contribution rate.

[0026] The potential weight calculation submodule extracts a mineralization potential ranking table based on the tectonic feature contribution rate, analyzes the adaptation degree of the mineralization potential weight to the ore-forming environment, and generates an environment adaptation index.

[0027] As a further scheme of the present application, the spatial feature analysis module comprises:

[0028] The spatial continuity analysis submodule calls the environment adaptation index, collects multi-source spatial data of the mineralization area, interprets the spatial distribution unit of the target area, identifies the strength of the continuity of adjacent units, divides the continuity level blocks, analyzes the differentiated level block combination mode, and generates a spatial continuity index;

[0029] The boundary feature extraction submodule positions the boundary trend line of the mineralization area according to the spatial continuity index, measures the boundary definition distribution value range, calibrates the boundary transition band width interval, decomposes the boundary point set of the differentiated interval, analyzes the position relationship of the boundary point set corresponding to the main trunk area, and generates a boundary feature definition;

[0030] The distribution optimization scoring module sorts out the combination sequence of the spatial distribution unit based on the boundary feature definition, compares the boundary fitting state of adjacent units, filters the distribution continuity mutation area, marks the high-difference distribution unit, evaluates the influence of the high-difference unit on the overall distribution coupling degree, and generates a spatial distribution optimization set.

[0031] As a further scheme of the present application, the multi-source spatial data refers to a set of geographic spatial information from differentiated sensors or data platforms;

[0032] The spatial continuity index refers to an index of the strength of the distribution continuity between adjacent spatial units;

[0033] The boundary feature definition refers to an index of the degree of boundary form definition and structure certainty.

[0034] As a further scheme of the present application, the system further comprises a potential area determination module:

[0035] The potential area determination module analyzes the regional metallogenic potential level based on the spatial distribution optimization set, combines the original exploration data and the specification requirements, and obtains the rare metal ore metallogenic potential area evaluation result;

[0036] The rare metal ore metallogenic potential area evaluation result comprises potential level division, regional reliability score and exploration consistency label.

[0037] As a further scheme of the present application, the potential area determination module comprises:

[0038] The potential probability analysis submodule extracts the metallogenic potential index and spatial distribution features in the original exploration data based on the spatial distribution optimization set, analyzes the frequency and regularity of the metallogenic potential occurrence, and generates a potential probability interval;

[0039] The metallogenic intensity evaluation submodule calls the potential probability interval, extracts the influence degree of metallogenic events on resource reserves and quality, analyzes the distribution characteristics of metallogenic intensity, and generates a metallogenic intensity grade;

[0040] The potential score submodule analyzes the distribution law of the regional potential score based on the metallogenic intensity grade and in combination with the potential grade division standard in the specification requirement, and generates a rare metal mineralization potential region evaluation result.

[0041] The technical scheme provided by the embodiment of the present application has at least the following beneficial effects:

[0042] By introducing the multi-source data analysis and fusion method, the geological structure, mineral distribution, metallogenic environment and spatial characteristics can be systematically associated and optimized, a new idea is used to extract the metallogenic potential region, the comprehensive evaluation capability and precision of resource exploration are effectively improved, the analysis after information fusion is more profound, the correlation analysis and spatial pattern extraction of different types of data are combined, the metallogenic potential evaluation in the mineral exploration process is more accurate, the limitations of single data source and manual analysis are avoided, and the precision and reliability of the metallogenic rule recognition are improved. By using the correlation between the geological structure and the mineral enrichment degree, more effective potential regions can be identified, and the analysis of the resource potential is optimized in combination with the adaptability of the metallogenic environment, the determination efficiency and accuracy of the metallogenic region are greatly improved, the problems such as information isolation and analysis deviation in the traditional method are avoided, and the scientificity and practical application value in the exploration process are enhanced. BRIEF DESCRIPTION OF DRAWINGS

[0043] In order to more clearly illustrate the technical scheme in the embodiment of the present application, the drawings needed in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0044] Figure 1 is a schematic diagram of the niobium-tantalum metal mineralization evaluation system provided by the embodiment of the present application based on multi-source information analysis;

[0045] Figure 2 is a system framework schematic diagram of the present application;

[0046] Figure 3 is a geological structure correlation module flow chart in the present application;

[0047] Figure 4 is a mineral distribution interpretation module flow chart in the present application;

[0048] Figure 5 is a metallogenic environment adaptation module flow chart in the present application;

[0049] Figure 6 The flow chart of the spatial feature analysis module in the present application is shown in the figure.

[0050] Figure 7 The flow chart of the potential area determination module in the present application is shown in the figure. DETAILED DESCRIPTION

[0051] The technical solutions in the present application will be described below with reference to the accompanying drawings.

[0052] In the embodiments of the present application, the words such as "example", "for example" are used to represent an example, illustration or description. Any embodiment or design scheme described as "example" in the present application should not be interpreted as more preferred or more advantageous than other embodiments or design schemes. In fact, the word "example" is intended to present the concept in a specific way. In addition, in the embodiments of the present application, the meaning expressed by "and / or" can be both, or can be one of the two.

[0053] In the embodiments of the present application, "image" and "picture" can be used interchangeably at times. It should be pointed out that when the distinction is not emphasized, the meanings expressed are consistent. "Of", "corresponding" and "corresponding" can be used interchangeably at times. It should be pointed out that when the distinction is not emphasized, the meanings expressed are consistent.

[0054] In the embodiments of the present application, sometimes the subscript such as W1 can be written in the form of non-subscript such as W1. When the distinction is not emphasized, the meanings expressed are consistent.

[0055] In order to make the technical problems, technical solutions and advantages to be solved by the present application more clear, the following will be described in detail with reference to the accompanying drawings and specific embodiments.

[0056] The embodiments of the present application provide a niobium-tantalum metal ore mineralization evaluation system based on multi-source information analysis, as shown in the figure. Figures 1-2 The schematic diagram of the niobium-tantalum metal ore mineralization evaluation system based on multi-source information analysis is shown in the figure. The system comprises:

[0057] The geological structure correlation module acquires geological structure information, including stratum distribution, fault position and fold shape, analyzes the influence degree of tectonic stress field on mineralization, and generates a tectonic stress influence index;

[0058] The mineral distribution interpretation module extracts mineral combination rules and chemical composition distribution characteristics based on the tectonic stress influence index, analyzes the corresponding relationship between mineral enrichment area and tectonic stress field, and obtains a mineral enrichment degree value;

[0059] The ore-forming environment adaptation module extracts the stratum lithology and structural characteristics in the ore-forming environment according to the mineral enrichment value, analyzes the correlation between the lithology combination and the ore-forming potential, and generates an environment adaptation index;

[0060] The spatial feature analysis module calls the environment adaptation index, extracts the spatial distribution pattern and boundary characteristics of the ore-forming area, analyzes the relationship between the spatial continuity and the ore-forming potential, and generates a spatial distribution optimization set;

[0061] The potential area determination module analyzes the regional ore-forming potential grade based on the spatial distribution optimization set, in combination with the original exploration data and the specification requirements, and obtains the rare metal ore-forming potential area evaluation result;

[0062] The structural stress influence index includes stress distribution balance, structural disturbance intensity and stress action directionality, the mineral enrichment value includes enrichment intensity grade, mineral concentration coefficient and spatial aggregation distribution degree, the environment adaptation index includes lithology adaptability score, structural coordination degree and ore-forming environment consistency, the spatial distribution optimization set includes distribution continuity score, boundary integrity index and regional shape regularity, and the rare metal ore-forming potential area evaluation result includes potential grade division, regional reliability score and exploration consistency label.

[0063] Specifically, as shown in Figure 2 , 3 , the geological structure correlation module includes:

[0064] The stratum distribution extraction submodule obtains the geological structure information, extracts the stratum thickness variation and lithology interface characteristics, analyzes the influence degree of stratum fluctuation on ore-forming effect, and generates fault activity intensity;

[0065] Obtaining geological structure information, integrating seismic exploration data, drilling core data and ground geological survey data to obtain a three-dimensional geological structure model of the target area, extracting stratum thickness variation and lithology interface characteristics, the stratum thickness data revealed by drilling hole ZK01 in a certain mining area shows that the thickness of Silurian sandstone is 350 meters, and the underlying Ordovician limestone is 200 meters thick, and the lithology interface is the unconformity surface between sandstone and limestone. The unconformity surface shows erosion and truncation characteristics. Analyzing the influence of stratum fluctuation on mineralization, through the stratum isopach map and the stratum pinch-out line map, the area where the stratum thickness is less than 100 meters and the lithology interface inclination is greater than 15 degrees is determined as the strong stratum fluctuation area. The strong stratum fluctuation area causes significant hindrance and diversion effect on the migration of ore fluid, affecting the precipitation and enrichment of ore, and generates fault activity intensity. The fault activity intensity is defined as the weighted average of the number of faulted strata and the fault displacement, wherein the weight of the number of faulted strata is set to 0.6, and the weight of the displacement is set to 0.4. For example, a certain F1 fault faulted 3 main strata, and the maximum vertical displacement is 50 meters. The fault activity intensity is calculated as 3x0.6+50x0.4=21.8. The larger the fault activity intensity value, the stronger the control of the fault on mineralization.

[0066] The fold morphology analysis submodule extracts the axial plane inclination and wavelength variation value in the fold morphology according to the fault activity intensity, analyzes the degree of fold deformation on mineral migration and enrichment, and generates the fold deformation degree;

[0067] According to the fault activity intensity, relevant fold region data is extracted based on the fault activity intensity value. For example, when the fault activity intensity is greater than 15, the fold structure in the region is taken as the analysis object, the axial plane inclination and wavelength variation value in the fold morphology are extracted, and the geological profile and three-dimensional geological model of the target area are analyzed. In a certain fold structure, the axial plane inclination gradually changes from 65 degrees in the north wing to 45 degrees in the south wing, and the wavelength changes from 800 meters in the east to 600 meters in the west. The degree of fold deformation on mineral migration and enrichment is analyzed, and the area where the fold axial plane inclination is greater than 40 degrees and the wavelength is less than 700 meters is evaluated as a high fold deformation area. The fold deformation in this area causes changes in rock porosity and permeability, affecting the migration path of hydrothermal minerals. Fluids tend to accumulate in the fold core or hinge, and ore is enriched there. The fold deformation degree is calculated by the product of the axial plane inclination and the wavelength change rate, wherein the axial plane inclination is 55 degrees on average, and the wavelength change rate is (800-600) / 800=0.25. The fold deformation degree is calculated as 55x0.25=13.75. This value is used to quantify the potential impact of folds on mineral enrichment.

[0068] The stress field evaluation submodule extracts structural stress field distribution information based on the fold deformation degree, analyzes the influence range of structural stress concentration on mineralization, and generates a structural stress influence index;

[0069] Based on the fold deformation degree, the fold deformation degree value is retrieved. When the fold deformation degree is greater than 10, the structural stress field distribution information corresponding to the deformation degree is extracted. The stress tensor data of the study area is obtained through finite element simulation, including the maximum principal stress, the minimum principal stress direction and size. For example, in a certain area, the maximum principal stress is 80 MPa, the minimum principal stress is 20 MPa, and the directions are 30 degrees and 60 degrees with the fold axis, respectively. The influence range of structural stress concentration on mineralization is analyzed. According to the Mohr-Coulomb failure criterion, the rock failure probability under the current stress field is calculated. The area with a failure probability greater than 0.7 is determined as the stress concentration area. The stress concentration in this area causes rock fractures, providing channels for ore fluid circulation, thereby affecting the precipitation and distribution of ore, generating a structural stress influence index. The index is calculated by the product of the area ratio of the stress concentration area and the stress intensity difference value. The stress intensity difference value is the difference between the maximum principal stress and the minimum principal stress. For example, the stress concentration area is 10 square kilometers, the total area of the study area is 100 square kilometers, and the area ratio is 0.1. The stress intensity difference value is , and the structural stress influence index is calculated as The index reflects the comprehensive influence of structural stress on mineralization.

[0070] Specifically, as shown in Figure 2 、 4 , the mineral distribution interpretation module includes:

[0071] The mineral combination extraction submodule collects mineral combination sample data based on the structural stress influence index, decomposes the characteristics of mineral particles in the sample, determines the content proportion of mineral types, extracts the chemical composition values of mineral particles, and compares the distribution of adjacent particles in differential mineral combinations. Identify the paragenetic mineral combination block and generate mineral species diversity.

[0072] Based on the tectonic stress influence index, when the tectonic stress influence index is greater than 5, the mineral assemblage sample data is collected, the rock samples of the target area are obtained through field outcrop sampling and drilling core sampling, for example, 10 hand specimens and 50 cm drilling cores are collected in a certain mineralized alteration zone, the characteristics of the mineral particles in the sample are decomposed, the sample thin sections are observed by using an optical microscope and a scanning electron microscope, the size, shape and crystal structure characteristics of the mineral particles are obtained, for example, euhedral pyrite and granular aggregate of galena are identified, the content proportion of the mineral types is determined, the volume percentage of the main minerals in the sample is quantitatively analyzed by image analysis and X-ray diffraction (XRD) technology, for example, the content of quartz in a certain sample is 30%, the content of pyrite is 20%, and the content of galena is 15%, the chemical composition values of the mineral particles are extracted, the single mineral particle is analyzed by using an electron probe (EPMA), the content of the major elements and trace elements is obtained, for example, the content of Fe in pyrite is 46.5%, the content of S is 53.2%, and the trace Ag content is 100 ppm, the distribution of adjacent particles in the differential mineral assemblage is compared, the spatial relationship between different mineral particles is identified by using the element surface scanning map and the mineral distribution map, for example, the pyrite particles are often associated with galena particles, the boundary is straight, and the paragenetic mineral assemblage block is identified, according to the similar mineral assemblage characteristics, the regions with similar mineral assemblage proportion and spatial distribution mode are divided into paragenetic mineral assemblage blocks, for example, the pyrite-galena-quartz assemblage is divided into a block, the mineral type diversity is generated, the diversity is calculated by the weighted average value of the number of mineral types in each block and the uniformity of the content of the main minerals, for example, a certain block contains 5 kinds of minerals, the uniformity of the content of the main minerals (measured by Shannon index) is 0.8, and the mineral type diversity is 5x0.5+0.8x0.5=2.9, the weight is 0.5.

[0073] The chemical composition analysis submodule selects the key elements in the differential mineral assemblage according to the mineral type diversity, calculates the stable fluctuation interval of the elements, extracts the stable element proportion in the sample, extracts the variable element content change section, discriminates the differential element assemblage relationship, and determines the influence degree of the differential assemblage on the sample distribution, and generates the chemical composition stability.

[0074] According to the mineral species diversity, the key elements in the differential mineral combination with a mineral species diversity greater than 2.5 are screened, and the elements with a content greater than 5% in each mineral combination block are selected as the main analysis objects. For example, in the pyrite-galena-quartz combination, Fe, S, Pb, Si, and O are selected as the key elements, the element stable fluctuation interval is calculated, the normal fluctuation range of the content of each key element is determined through statistical analysis of a large amount of sample data, for example, the stable fluctuation interval of the content of Fe is set to 40% to 50%, and the fluctuation within the interval is considered to be stable, the proportion of stable elements in the sample is extracted, and the total content of the key elements within the stable fluctuation interval in the sample accounts for the proportion of the total content of all key elements, for example, the content of Fe in a sample is 45%, the content of S is 52%, the content of Pb is 12%, the content of Si is 25%, and the content of O is 30%, among which Fe, Si, and O are within the stable fluctuation interval, and the proportion is , the variable element content change section is extracted, the elements whose content exceeds the stable fluctuation interval are identified, and the percentage of the excess is recorded, for example, the content of S is 52%, which exceeds the stable interval (assuming 45%-50%) by 2%, and the content of Pb is 12%, which also exceeds the stable interval (assuming 8%-10%) by 2%, the differential element combination relationship is distinguished, the paragenesis, association or antagonism relationship between elements is identified through principal component analysis and cluster analysis, for example, Fe and S show a significant positive correlation, and Pb and Ag show an association relationship, the influence degree of the differential combination on the sample distribution is marked, the influence of the element combination relationship and the content change section on mineralization enrichment or depletion is evaluated, for example, the stable paragenetic relationship of Fe-S indicates primary sulfide mineralization, and the abnormal enrichment of Pb-Ag indicates later superimposed modification or high-grade mineralization, the chemical composition stability is calculated by weighted average of the stable element proportion and the variable element content change section, the weight of the stable element proportion is set to 0.7, and the weight of the reciprocal of the variable element content change section is set to 0.3, for example, the stable element proportion is 0.61, the total variable element change percentage is 4% (2% of S plus 2% of Pb), and the reciprocal of the variable element content change section is , the chemical composition stability is .

[0075] The enrichment zone distribution scoring sub-module calls the chemical composition stability, divides the mineral enrichment zone sample unit, measures the key mineral content density in the enrichment unit, superimposes the unit and the corresponding range of tectonic stress block, identifies the boundary of the high-density enrichment block, quantifies the demarcation point set of the enrichment block and the low-density area, and obtains the mineral enrichment degree value;

[0076] The mineral enrichment degree value is calculated by the formula:

[0077] ;

[0078] wherein, Indicates the mineral enrichment value. This represents the content density of the i-th mineral sample. This represents the distribution density of the i-th mineral sample. This represents the amount of mineral per unit volume of the i-th mineral sample. This represents the average mineral content per unit volume across all units. Represents the total number of mineral samples;

[0079] The chemical composition stability is used. When the chemical composition stability is greater than 7.5, the sample units of the mineral enrichment area are divided. Based on information such as mineralization alteration characteristics, mineral assemblage, and chemical composition stability, the study area is divided into square sample units with a side length of 50 meters. For example, a mining area is divided into 100 grid units of 50m × 50m. The content density of key minerals in the enrichment unit is measured. Through borehole core analysis, geophysical anomaly interpretation, and surface sampling results, the average content of the target rare metal mineral in each unit is obtained. For example, the average content of tantalum (Ta) in a certain unit is 300 ppm, which is converted to a content density of... (Assuming the rock density is 3) The superimposed units correspond to the tectonic stress blocks. Each sample unit is spatially superimposed with the previously generated tectonic stress blocks to determine whether each unit is located in a high-stress concentration area. For example, unit A is located in a stress concentration area, and unit B is located in a stress dispersion area. The boundaries of high-density enrichment blocks are identified, and mineral content densities exceeding a preset threshold are identified (e.g., tantalum content density greater than a certain threshold). The adjacent units of the high-density enrichment area are connected to form a high-density enrichment block, and its spatial boundary line is drawn. The set of boundary points between the enrichment block and the low-density area is quantified. Through the boundary detection algorithm, the high-density enrichment block and the surrounding low-density area (e.g., tantalum content density less than 1.5%) are identified. The specific points between these points form a clear boundary line, from which the mineral enrichment value is obtained, and the formula is used. In the calculation, the parameters in the formula... This value represents the mineral enrichment level, a comprehensive indicator used to quantify the spatial concentration of mineral deposits. Representing the The content density of each mineral sample, in units of This represents the mass of the target mineral per unit volume. Representing the Distribution density of each mineral sample, in units of This indicates the number of mineral particles per unit volume or the density at proven depths. These two parameters together reflect the actual occurrence of the mineral. Representing the The amount of mineral per unit volume of a mineral sample, in units of , which is obtained by weighing the mineral sample and dividing it by its volume, , which represents the total number of mineral samples, i.e., the total number of discrete mineral sample points participating in the calculation, , which represents the average mineral mass per unit volume of all cells, , which is the average of all values, used to measure the average mineral mass level of the entire area, , which represents the sum of all terms from 1 to ;

[0080] First part of the formula , which represents the weighted average content density, where is the content density, is the distribution density, this part calculates the average content of the mineral considering the distribution density, the higher the distribution density, the greater the contribution to the average value, reflecting the concentration of the mineral in space, the second part of the formula , which represents the standard deviation of the mineral mass per unit volume, where is the mass of the individual sample, is the average mass of all samples, is the number of samples, this part quantifies the degree of dispersion or uniformity of the mineral mass in spatial distribution, the greater the standard deviation, the greater the fluctuation of the mineral mass, the higher the enrichment or the more uneven, the formula multiplies these two parts, so that the mineral enrichment value takes into account both the average content level of the mineral and the concentration or fluctuation degree of its distribution, a certain mineral area takes 3 representative mineral sample points;

[0081] Sample 1: ;

[0082] Sample 2: ;

[0083] Sample 3: ;

[0084] , first calculate ;

[0085] Then calculate the first part :

[0086] ;

[0087] Next, calculate the second part :

[0088] ;

[0089] Finally, ;​

[0090] The formula more comprehensively reflects the spatial enrichment characteristics of the mineral by combining the content density, the distribution density and the dispersion of the mass per unit volume, and avoids the problem of ignoring the uniformity by only relying on the high content. The result shows that the mineral enrichment degree value is 107.57, indicating that the mineral in the region exists in a moderate degree of enrichment. The value is an important index for measuring the occurrence characteristics of the mineral.

[0091] Specifically, as shown in Figure 2 、 5 The ore-forming environment adaptation module includes:

[0092] The lithological combination extraction submodule extracts the stratigraphic lithological combination type in the ore-forming environment according to the mineral enrichment degree value, analyzes the support degree of the lithological combination to the mineral enrichment, and generates a lithological combination matching degree.

[0093] According to the mineral enrichment degree value, when the mineral enrichment degree value is greater than 100, the stratigraphic lithological combination type in the ore-forming environment is extracted, the rock types related to mineralization and their spatial combination relationship are identified by comprehensively interpreting the geological map, the drill columnar section and the well logging curve, for example, the granite-walling rock contact zone, the metamorphic sandstone-schist interbedding combination and the like are identified, the support degree of the lithological combination to the mineral enrichment is analyzed, the geochemical background, the physical properties (such as porosity, permeability) and the affinity with the paragenetic minerals of each lithological combination are evaluated, the favorable or unfavorable influence of the lithological combination on the enrichment of the target mineral is quantified, for example, the felsic intrusive rock has a high support degree to the enrichment of the rare metal ore due to the rich rare metal elements and the easy formation of fissures, the carbonate rock has a strong reactivity and adsorbs the ore-forming materials to generate a lithological combination matching degree. The matching degree is calculated by the quantified score of the support degree of the lithological combination to the mineral enrichment, for example, the matching degree of the granite-walling rock contact zone is scored as 0.8, the matching degree of the metamorphic sandstone-schist interbedding is scored as 0.6, and the higher the matching degree, the more favorable the lithological combination is to the mineral enrichment.

[0094] The structural feature matching submodule calls the lithological combination matching degree, extracts the fault strike and the fold shape in the structural feature, analyzes the contribution degree of the structural feature to the mineralization, and generates a structural feature contribution rate.

[0095] When the lithology combination matching degree is greater than 0.7, the strike of the fault and the fold shape in the structural feature are extracted, the three-dimensional spatial distribution direction and the dip angle of the fault and the geometric elements such as the axial direction, the plunging direction and the wing dip angle of the fold in the research region are obtained by analyzing the high-precision gravity and magnetic data, the aerial remote sensing image and the geological surveying and mapping data, the contribution degree of the structural feature to the mineralization is analyzed, the effectiveness of the fault as the ore fluid migration channel is evaluated, the control ability of the fold hinge and the wing to the ore fluid accumulation and precipitation is calculated, for example, the ore body controlled by the north-east fault zone is generally higher in grade than the north-west fault, the turning end of the fold is the preferred position of the ore body enrichment, the structural feature contribution rate is generated, the contribution rate is calculated by the weighted average of the fault ore migration capacity and the fold ore control capacity, wherein the weight of the fault ore migration capacity is set to 0.6 and the weight of the fold ore control capacity is set to 0.4, for example, the ore migration capacity score of a certain fault is 0.75, the ore control capacity score of a certain fold is 0.65, and then the structural feature contribution rate is 0.75*0.6+0.65*0.4=0.71, the higher the contribution rate, the greater the contribution of the structural feature to the mineralization.

[0096] The potential weight calculation sub-module extracts the ore-forming potential ranking table based on the structural feature contribution rate, analyzes the adaptation degree of the ore-forming potential weight to the ore-forming environment, and generates an environment adaptation index.

[0097] When the structural feature contribution rate is greater than 0.7, the ore-forming potential ranking table is extracted, a standard for classifying the ore-forming potential according to different geological, geochemical and geophysical feature combinations is established through geological expert experience, historical exploration data and regional ore-forming regularity analysis, and a potential weight distribution table is formed, for example, containing "high potential area: weight 0.9, medium potential area: weight 0.6, low potential area: weight 0.3" and the like, the adaptation degree of the ore-forming potential weight to the ore-forming environment is analyzed, the geological features of the current analysis region are compared with each index in the potential ranking table, and the corresponding potential weight is distributed according to the matching degree, for example, if a certain region has similar fault structure and lithology combination as the high potential area, it is given a potential weight of 0.9, and an environment adaptation index is generated, the index is calculated by the product of the potential weight and the analyzed geological features (such as lithology matching degree, structural feature contribution rate), for example, the potential weight of a certain region is 0.9, the lithology combination matching degree is 0.8, and the structural feature contribution rate is 0.71 (from the previous step), and then the environment adaptation index is 0.9*0.8*0.71=0.5112, which comprehensively reflects the comprehensive favorable degree of the ore-forming environment.

[0098] Specifically, as shown in Figure 2 , 6 , the spatial feature analysis module comprises:

[0099] The spatial continuity analysis submodule calls the environment adaptation index, collects multi-source spatial data in the mineralization area, interprets the spatial distribution unit of the target area, identifies the strength of the continuity of adjacent units, divides the continuity level blocks, analyzes the combination mode of the differentiated level blocks, and generates the spatial continuity index;

[0100] Multi-source spatial data refers to a set of geographic spatial information from different sensors or data platforms;

[0101] The spatial continuity index refers to an index of the strength of the continuity between adjacent spatial units;

[0102] When the environment adaptation index is greater than 0.5, the multi-source spatial data in the mineralization area is collected, the different sources of spatial data are geographically registered and fused by obtaining satellite remote sensing images (such as Sentinel 2 multispectral images), terrain elevation data (DEM), and drill hole distribution data, the spatial distribution unit of the target area is interpreted, the study area is divided into regular grid units with a side length of 200 meters based on the fused multi-source data through image segmentation and geological unit demarcation, each unit represents a spatial distribution unit, for example, in a 100 square kilometer study area, 2500 200mx200m grid units are divided, the strength of the continuity of adjacent units is identified, the difference value of the environment adaptation index between adjacent grid units is calculated, the difference value is less than 0.1, which is determined as strong continuity, the difference value is between 0.1 and 0.3, which is determined as medium continuity, and the difference value is greater than 0.3, which is determined as weak continuity, for example, the environment adaptation index of unit A is 0.55, the adjacent unit B is 0.53, the difference is 0.02, which is strong continuity, unit C is 0.48, the difference is 0.07, which is strong continuity, and unit D is 0.35, the difference is 0.20, which is medium continuity, the continuity level blocks are divided, adjacent units with similar continuity strength are aggregated to form continuity level blocks, for example, all strong continuity units are aggregated into a high continuity block, the combination mode of the differentiated level blocks is analyzed, the spatial adjacency relationship, geometric morphology and scale between different continuity level blocks are analyzed, for example, the high continuity block is distributed in a strip shape and is consistent with the strike of the fault structure, and the low continuity block is distributed in a sporadic shape, the spatial continuity index is generated, which is calculated by the weighted average of the area proportion of the strong continuity area and the average environment adaptation index of the block, the weight of the area proportion of the strong continuity area is set to 0.7, and the weight of the average environment adaptation index of the block is set to 0.3, for example, the area proportion of the strong continuity area in the total area of the study area is 0.6, the average environment adaptation index of all blocks is 0.52, and the spatial continuity index is 0.6x0.7+0.52x0.3=0.576, which reflects the uniformity and continuity of the mineralization environment in the spatial distribution.

[0103] The boundary feature extraction submodule locates the boundary trend line segment of the ore-forming area according to the spatial continuity index, measures the boundary definition value range, calibrates the boundary transition zone width interval, decomposes the differentiated interval boundary point set, analyzes the position relationship of the boundary point set corresponding to the main trunk area, and generates the boundary feature definition;

[0104] The boundary feature definition is an index of the definition degree of the regional boundary form and the structure certainty.

[0105] When the spatial continuity index is greater than 0.5, the boundary trend line segment of the ore-forming area is located, the boundary line segment representing the core of the ore-forming area is extracted by vectorizing the edge of the high continuity block, for example, a 2000-meter-long north-east boundary line is extracted on the edge of a certain long strip-shaped high continuity block, the boundary definition value range is measured, the boundary line segment is sampled at an interval of 10 meters, the standard deviation of the spatial continuity index within a radius of 50 meters around each sampling point is calculated, the standard deviation less than 0.05 is high definition, 0.05 to 0.1 is medium definition, and greater than 0.1 is low definition, the length proportion of the high definition, medium definition and low definition areas is calculated, for example, 70% of the length of the boundary line segment is high definition, 20% is medium definition, and 10% is low definition, the boundary transition zone width interval is calibrated, for each boundary sampling point, the distance from the high value to the low value of the spatial continuity index is searched outward and inward, the distance is defined as the transition zone width, and the minimum value, maximum value and average value of all transition zone widths are calculated, for example, the transition zone width interval is 50 meters to 150 meters, the differentiated interval boundary point set is decomposed, the boundary line segment is divided into different point sets according to the boundary definition and the transition zone width, for example, the point set of high definition and narrow transition zone, the point set of low definition and wide transition zone, the position relationship of the boundary point set corresponding to the main trunk area is analyzed, the spatial relationship between each boundary point set and the known ore body or main mineralization zone in the area is evaluated, for example, the high definition boundary point set is adjacent to the main trunk ore body, and the low definition boundary point set is far away, the boundary feature definition is generated, the definition is calculated by the difference between the high definition boundary length proportion and the low definition boundary length proportion, for example, the high definition boundary length proportion is 0.7, the low definition boundary length proportion is 0.1, and the boundary feature definition is 0.7-0.1=0.6, the higher the value, the more certain the boundary form of the ore-forming area.

[0106] The distribution optimization scoring submodule sorts the spatial distribution unit combination sequence based on the boundary feature definition, compares the fitting state of adjacent unit boundaries, filters the distribution continuity mutation area, marks the high difference distribution unit, evaluates the influence of the high difference unit on the overall distribution coupling degree, and generates the spatial distribution optimization set.

[0107] Based on the boundary feature clarity, when the boundary feature clarity is greater than 0.5, the spatial distribution unit combination sequence is combed, all grid units are sorted and grouped according to the similarity and spatial proximity of their geological features according to the spatial continuity index and the boundary feature clarity, for example, units with high environmental adaptation index and strong spatial continuity are classified into one category, the fitting state of adjacent unit boundaries is compared, it is checked whether the boundary line between adjacent units is smooth, continuous, and coincides with the actual geological boundary, for example, the boundary of two adjacent units coincides with a fault line by 90%, it is considered that the fitting state is good, the distribution continuity mutation area is screened, the spatial continuity index or the environmental adaptation index suddenly decreases greatly between adjacent units, for example, the boundary of a unit with an environmental adaptation index of 0.8 suddenly decreases to a unit with an environmental adaptation index of 0.2, which is a mutation area, mark the high difference distribution unit, mark those isolated units whose environmental adaptation index is significantly different from the surrounding units and do not belong to any known geological body boundary, for example, an isolated unit with an environmental adaptation index of 0.9 is surrounded by units with an environmental adaptation index of 0.3, evaluate the influence of high difference units on the overall distribution coupling degree, analyze whether the high difference unit represents an unidentified mineralization anomaly or a measurement error, quantify its interference or promotion effect on the whole area metallogenic potential evaluation, and generate a spatial distribution optimization set. The optimization set considers the comprehensive evaluation of the number of distribution continuity mutation areas and the influence degree of high difference distribution units, the optimized spatial distribution unit sequence is used for subsequent potential evaluation, for example, the spatial distribution optimization set contains all high continuity and high environmental adaptation areas, and excludes units with uncertain boundaries or high differences.

[0108] Specifically, as shown in Figure 2 , 7 , the potential area determination module includes:

[0109] The potential probability analysis submodule extracts the metallogenic potential indicators and spatial distribution features in the original exploration data based on the spatial distribution optimization set, analyzes the frequency and regularity of the occurrence of metallogenic potential, and generates a potential probability interval.

[0110] Based on the spatial distribution optimization set, the metallogenic potential indicators and spatial distribution characteristics in the original exploration data are extracted. From historical drilling data, geophysical and geochemical anomaly maps, and geological mapping data, the discovered ore point grade, ore body thickness, and anomaly intensity are obtained as metallogenic potential indicators, and the spatial position information of the indicators is extracted, for example, a drill hole reveals a tantalum (Ta) grade of 500 ppm and an ore body thickness of 10 meters, which is located at the center of a geochemical anomaly. The frequency and regularity of metallogenic potential occurrence are analyzed, and the number of high-grade mineralization points and their spatial distribution patterns under different geological backgrounds, structural positions, and lithological combinations are counted, for example, 80% of high-grade tantalum mineralization points occur near the northeast-trending faults at the contact between granite and metamorphic rock. The potential probability interval is generated. Through statistical analysis, the probability range of discovering high-grade rare metal ore in different types of geological units is determined, for example, in the granite-metamorphic rock contact zone, the potential probability interval is 0.6 to 0.8, and in the non-contact zone, the potential probability interval is 0.1 to 0.3. This interval reflects the likelihood of discovering mineral resources under specific geological conditions.

[0111] The metallogenic intensity evaluation submodule calls the potential probability interval, extracts the influence degree of the metallogenic event on resource reserves and quality, analyzes the distribution characteristics of the metallogenic intensity, and generates the metallogenic intensity grade.

[0112] When the potential probability interval is greater than 0.5, the influence degree of the metallogenic event on resource reserves and quality is extracted, the resource quantity (tons) and average grade (ppm or %) of the explored ore body are calculated, and the deposit size and ore quality formed by the metallogenic event are quantified, for example, a certain ore body has a proven resource quantity of 1 million tons and an average grade of 300 ppm tantalum. The distribution characteristics of the metallogenic intensity are analyzed, and the distribution patterns of high-grade enrichment cores, medium-grade mineralization zones, and low-grade halos are identified based on ore body shape, occurrence, and grade distribution, for example, a high-grade enrichment core is located at the intersection of faults, and the grade decreases outward. The metallogenic intensity grade is generated. According to the resource reserves, average grade, and uniformity of grade distribution, the metallogenic intensity is divided into different grades, for example, a resource quantity greater than 1 million tons and an average grade greater than 200 ppm are classified as a first-grade metallogenic intensity zone, and a resource quantity of 50-100 million tons and a grade of 100-200 ppm are classified as a second-grade metallogenic intensity zone. The higher the metallogenic intensity grade, the higher the scale and quality of mineral resources in the region.

[0113] The potential scoring submodule generates rare metal mineralization potential regional evaluation results based on the metallogenic intensity grade and the potential grade division standard in the specification requirements, analyzes the distribution pattern of regional potential scores, and generates rare metal mineralization potential regional evaluation results.

[0114] Based on the ore-forming intensity grade, combined with the potential grade division standard in the specification requirement, the determined ore-forming intensity grade is compared with the rare metal mineral resource potential evaluation standard stipulated by the state or industry, for example, referring to the rare metal mineral exploration degree grade division standard in the solid mineral exploration specification, the distribution law of the regional potential score is analyzed, the ore-forming intensity grade is combined with the potential probability interval, the final potential score of each evaluation unit is calculated through the geological model and the geostatistics method, and the regional potential score contour map is drawn, for example, the score of the high-intensity ore-forming area combined with the high-potential probability area is significantly higher than the region, the rare metal mineralization potential regional evaluation result is generated, the rare metal mineralization potential distribution map and the potential grade report in the research region are finally determined and output, and the 'extremely high potential area', 'high potential area','medium potential area' and 'low potential area' are divided, for example, a certain area is evaluated as an extremely high potential area, and it is predicted that the future exploration investment will bring significant resource discovery.

[0115] The above is only a specific embodiment of the present application, but the protection scope of the present application is not limited thereto, any person skilled in the art can easily think of changes or replacements within the technical range disclosed by the present application, which should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A system for evaluating niobium-tantalum metal mineralization based on multi-source information analysis, characterized in that, The system comprises: The geological structure correlation module obtains geological structure information, including stratum distribution, fault position and fold shape, analyzes the influence degree of tectonic stress field on mineralization, and generates a tectonic stress influence index; The mineral distribution interpretation module extracts mineral assemblage rules and chemical composition distribution characteristics based on the tectonic stress influence index, analyzes the corresponding relationship between mineral enrichment area and tectonic stress field, and obtains a mineral enrichment degree value; The mineral distribution interpretation module comprises: The mineral assemblage extraction submodule collects mineral assemblage sample data based on the tectonic stress influence index, decomposes the characteristics of mineral particles in the sample, determines the content proportion of mineral types, extracts the chemical composition values of mineral particles, compares the distribution of adjacent particles in the differential mineral assemblage, identifies the paragenetic mineral assemblage block, and generates mineral type diversity; The chemical composition analysis submodule filters key elements in the differential mineral assemblage according to the mineral type diversity, calculates the stable fluctuation interval of the elements, extracts the proportion of stable elements in the sample, extracts the content variation section of the variable elements, discriminates the differential element combination relationship, marks the influence degree of the differential combination on the sample distribution, and generates chemical composition stability; The enrichment area distribution scoring submodule calls the chemical composition stability, divides the mineral enrichment area sample unit, determines the content density of key minerals in the enrichment unit, superimposes the corresponding range of the unit and the tectonic stress block, identifies the boundary of the high-density enrichment block, quantizes the demarcation point set of the enrichment block and the low-density area, and obtains the mineral enrichment degree value; The ore-forming environment adaptation module extracts the stratum lithology and structural characteristics in the ore-forming environment according to the mineral enrichment degree value, analyzes the relevance of lithology combination and mineralization potential, and generates an environment adaptation index; The ore-forming environment adaptation module comprises: The lithology combination extraction submodule extracts the stratum lithology combination type in the ore-forming environment according to the mineral enrichment degree value, analyzes the support degree of lithology combination to mineral enrichment, and generates a lithology combination matching degree; The structural feature matching submodule calls the lithology combination matching degree, extracts the fault strike and fold shape in the structural feature, analyzes the contribution degree of the structural feature to mineralization, and generates a structural feature contribution rate; The potential weight calculation submodule extracts an ore-forming potential ranking table based on the structural feature contribution rate, analyzes the adaptation degree of the ore-forming potential weight to the ore-forming environment, and generates an environment adaptation index; The spatial feature analysis module calls the environment adaptation index, extracts the spatial distribution pattern and boundary characteristics of the ore-forming area, analyzes the relationship between spatial continuity and ore-forming potential, and generates a spatial distribution optimization set; The spatial feature analysis module comprises: The spatial continuity analysis submodule collects multi-source spatial data of the ore-forming area by calling the environment adaptation index, interprets the spatial distribution unit of the target area, identifies the strength of the continuity of adjacent units, divides the continuity level block, analyzes the differential level block combination mode, and generates a spatial continuity index; The boundary feature extraction submodule locates a metallogenic region boundary trend line segment according to the spatial continuity index, measures a boundary definition value range, calibrates a boundary transition zone width interval, decomposes a differentiated interval boundary point set, analyzes a position relationship of the boundary point set corresponding to a main trunk area, and generates a boundary feature definition; The distribution optimization scoring submodule combs a spatial distribution unit combination sequence based on the boundary feature definition, compares adjacent unit boundary fitting states, screens a distribution continuity mutation area, marks a high-difference distribution unit, evaluates an influence of the high-difference unit on overall distribution coupling degree, and generates a spatial distribution optimization set; The potential region determination module analyzes a regional metallogenic potential grade based on the spatial distribution optimization set, in combination with original exploration data and specification requirements, and obtains a rare metal ore metallogenic potential region evaluation result. The rare metal ore metallogenic potential region evaluation result includes a potential grade division, a region reliability score, and an exploration consistency label.

2. The multi-source information analysis-based mineralization evaluation system for niobium-tantalum metal mines according to claim 1, characterized in that: The tectonic stress influence index includes stress distribution uniformity, tectonic disturbance intensity, and stress action directionality, the mineral enrichment degree value includes enrichment intensity grade, mineral concentration coefficient, and spatial aggregation distribution degree, the environment adaptation index includes lithology adaptability score, tectonic coordination degree, and metallogenic environment consistency, and the spatial distribution optimization set includes distribution continuity score, boundary integrity index, and region shape regularity. 3.The system for mineralization evaluation of niobium-tantalum ore based on multi-source information analysis according to claim 1, characterized in that: The geological structure correlation module includes: A stratigraphic distribution extraction submodule obtains geological structure information, extracts stratigraphic thickness variation and lithology interface features, analyzes an influence degree of stratigraphic relief on metallogenic effect, and generates fault activity intensity; A fold shape analysis submodule extracts axial plane dip angle and wavelength variation values in fold shape based on the fault activity intensity, analyzes an influence degree of fold deformation on mineral migration and enrichment, and generates fold deformation degree; A stress field evaluation submodule extracts tectonic stress field distribution information based on the fold deformation degree, analyzes an influence range of tectonic stress concentration on metallogenic effect, and generates a tectonic stress influence index.

4. The multi-source information resolution-based mineralization evaluation system for niobium-tantalum metal mines according to claim 1, characterized in that: The mineral enrichment degree value uses a formula: ； wherein, denotes the mineral enrichment value, denotes the content density of the i-th mineral sample, denotes the distribution density of the i-th mineral sample, denotes the mineral mass per volume of the i-th mineral sample, denotes the average of the mineral mass per volume of all cells, denotes the total number of mineral samples.

5. The multi-source information resolution-based mineralization evaluation system for niobium-tantalum metal mines according to claim 1, characterized in that: The multi-source spatial data refer to geographic spatial information sets from differentiated sensors or data platforms; The spatial continuity index refers to an index of strength and weakness of distribution continuity between adjacent spatial units; The boundary feature definition refers to an index of region boundary shape definition degree and structure certainty. 6.The system for mineralization evaluation of niobium-tantalum ore based on multi-source information analysis according to claim 1, characterized in that: The potential region determination module includes: A potential probability analysis submodule extracts a metallogenic potential index and spatial distribution features in original exploration data based on the spatial distribution optimization set, analyzes a frequency and regularity of metallogenic potential occurrence, and generates a potential probability interval; A metallogenic intensity evaluation submodule calls the potential probability interval, extracts an influence degree of a metallogenic event on resource reserves and quality, analyzes distribution features of metallogenic intensity, and generates a metallogenic intensity grade; A potential scoring submodule analyzes a distribution law of regional potential scores based on the metallogenic intensity grade in combination with potential grade division standards in specification requirements, and generates a rare metal ore metallogenic potential region evaluation result.

Images