CN-121980114-A - Multi-metal ore resource quantity estimation method and system thereof
Abstract
The invention relates to the technical field of mineral resource quantity calculation, in particular to a multi-metal mineral resource quantity estimation method and a system thereof, which uniformly restrict the spatial correlation of multi-metal grades in adjacent calculation units through a collaborative kriging method, so that the consistent relation of different metals in trend and trend change directions is brought into the same judgment system, the space splitting risk generated during the multi-metal calculation is reduced, the contribution ordering of the multi-metal grades in the same calculation unit is carried out through the Skerman grade correlation coefficient, the relative change relation between the metal grades is converted into a comparable grade sequence, the unstable factors of low-contribution metals in numerical conversion are avoided, the abrupt metal items are identified and removed through the spatial correlation constraint, the interference of the local abnormal grades on the whole result is reduced, and the formed multi-metal resource quantity result is enhanced in the aspects of spatial consistency, numerical comparability and overall reliability.
Inventors
- YANG HONGWEI
- LIU WEIJUN
- XU JIATIAN
- GUO ZIHAO
Assignees
- 四川省金属地质调查研究所
Dates
- Publication Date
- 20260505
- Application Date
- 20260123
Claims (10)
- 1. A method for estimating the amount of a polymetallic ore resource, comprising the steps of: s1, reading a coordinate positioning unit based on boundary coordinates, drilling section grades, unit volumes and unit densities recorded by a ore body model file, matching drilling sections to unit indexes, and binding a volume density writing sequence to obtain an ore body calculation unit space sequence; S2, based on the ore body calculation unit space sequence, extracting multi-metal grades of adjacent units, comparing trend and trend change directions and judging continuous range intervals according to space correlation constraint requirements of a collaborative Kriging method, removing abrupt metal items, registering symbiotic pairing and forming a multi-metal space correlation constraint structure; S3, based on the multi-metal space association constraint structure, a screening unit allows metal items and reads corresponding grade values, a retention item is selected after contribution sorting is calculated based on the Squalman grade correlation coefficient, and the retention item is written into a unit record after low contribution items are removed, so that a combined metal grade value set is generated; S4, checking grades and constraints based on the multi-metal space association constraint structure and the combined metal grade value set, selecting reserved metals, eliminating violating items, comparing adjacent difference values and smoothing abnormal units to obtain an effective mineralization calculation unit set; And S5, reading the unit volume and the unit density based on the effective mineralization calculation unit set, completing conversion of the paired metal grades, accumulating the unit conversion values according to the metal types, merging and summarizing according to the spatial index, and outputting a multi-metal resource quantity numerical list.
- 2. The method for estimating the multi-metal ore resource quantity according to claim 1, wherein the ore body calculation unit space sequence comprises a unit index, a unit space coordinate, a drilling section grade, a unit volume and a unit density, the multi-metal space association constraint structure comprises a metal change direction consistency mark, a change amplitude continuity mark and a symbiotic metal pairing record, the combined metal grade value set comprises a calculation unit index, a reserved metal mark and a corresponding metal grade value, the effective mineralization calculation unit set comprises a calculation unit index, a reserved metal mark and a corresponding metal grade value, and the multi-metal resource quantity value list comprises a metal category mark, a corresponding resource quantity value and a space index merging result.
- 3. The method for estimating the amount of polymetallic ore resources according to claim 1, wherein the specific steps of obtaining the spatial sequence of ore body calculation units are as follows: Based on boundary coordinates, drilling section grades, unit volumes and unit densities recorded by the ore body model file, reading the boundary coordinates, checking and closing to generate unit indexes, writing the drilling coordinates into the corresponding indexes, and obtaining a calculation unit positioning corresponding relation table; And reading the unit index and writing the grade of the drilling section based on the calculation unit positioning corresponding relation table, binding the unit volume and the unit density according to the index, writing the same record, and arranging the record items according to the index sequence to obtain the ore body calculation unit space sequence.
- 4. The method for estimating the amount of polymetallic mineral resources according to claim 1, wherein the specific steps of forming the polymetallic spatial association constraint structure are: Based on the ore body calculation unit space sequence, reading adjacent calculation unit records according to a space index sequence, extracting each metal grade value in a unit one by one, classifying and sorting metal categories, and arranging similar metal grades of adjacent units into continuous records to obtain an adjacent unit metal grade arrangement sequence; Based on the adjacent unit metal grade arrangement sequence, according to the calculation requirement of the collaborative kriging method on the multi-metal space correlation constraint, comparing the trend change direction of the same metal in the adjacent units and the trend change direction, calculating the adjacent grade value difference interval, and marking the metal items beyond the interval to obtain a metal change consistency marking set; and based on the metal change consistency marking set, carrying out item-by-item rejection on marked metal items, carrying out pairing registration on unmarked metal items according to the index of the computing unit, collecting and storing metal pairing results in the same unit, and establishing a multi-metal space association constraint structure.
- 5. The method for estimating the multi-metal ore resource amount according to claim 1, wherein the collaborative kriging method firstly calls a plurality of metal grade values participating in analysis in the range of adjacent calculation units, performs joint investigation on the change relation of each metal grade based on the same spatial position, performs constraint check on the change direction consistency, change amplitude continuity and grade difference distribution condition of the metal grade of the same type in the adjacent units according to the established spatial association relation between the metals, and is used as the calculation basis for the subsequent metal item reservation and rejection judgment.
- 6. The method for estimating the quantity of multi-metal ore resources according to claim 1, characterized in that the specific step of generating the set of combined metal grade values is: Based on the multi-metal space association constraint structure, reading a metal identifier corresponding to the index of the computing unit, judging that the identifier state allows participation in the current unit processing, selecting allowed identifier metal, reading a metal corresponding grade value, writing the metal corresponding grade value into a unit temporary storage area, and establishing a unit metal grade reading set; Based on the unit metal grade reading set, reading the grade values of all metals in the same calculation unit, carrying out grade correlation check on the metal grade sequence by adopting a Sherman grade correlation coefficient, outputting a sequencing basis, comparing the grade value magnitude relation to generate grade sequence numbers, marking the metals with the front sequence as a reserved state and marking the metals with the rear sequence as an exclusion state, and obtaining a unit metal accepting and rejecting sequencing record; Based on the unit metal sorting records, reading the grade values corresponding to the reserved state metal, writing the grade values into a formal recording area of the computing unit, removing the storage contents corresponding to the excluded state metal, and updating the unit metal recording sequence to generate a combined metal grade value set.
- 7. The method for estimating the resource quantity of the multi-metal ore according to claim 1, wherein the spearman grade correlation coefficient is characterized in that each metal grade value in the same calculation unit is read in a unit metal grade reading set, each metal grade value sequence is generated into a grade sequence number according to the value, the grade sequence number is used as a grade representation of the metal grade sequence, grade sequence numbers corresponding to different metals are paired item by item, the grade sequence number difference after pairing is summarized, a grade correlation check value is output, and finally the grade correlation check value is used as a sorting basis to be written into a unit metal sorting record.
- 8. The method for estimating the amount of polymetallic mineral resources according to claim 1, wherein the specific step of obtaining the set of effective mineralization calculating units is: Reading the metal grade records and checking the marks based on the multi-metal space association constraint structure and the combined metal grade value set, judging the reserved marks and clearing the violating items, collecting the rest records and writing the rest records into a list, and generating a constraint checking reserved metal record set; And (3) based on the constraint check, reserving a metal record set, reading similar metal grade values of adjacent units, calculating a difference value, comparing a threshold value, adjusting an overrun unit value, screening a stable unit index, and collecting records to obtain an effective mineralization calculation unit set.
- 9. The method for estimating the amount of polymetallic ore resources according to claim 1, wherein the specific step of outputting the numerical list of the amount of polymetallic ore resources is: Based on the effective mineralization calculation unit set, reading the volume value of each calculation unit, reading the corresponding density value, reading the metal grade value of the same unit, performing pairing calculation on the volume value and the density value, and writing the calculation result and the metal grade into the unit record in a item-by-item conversion way to generate a unit metal resource amount conversion record; and reading metal resource quantity values converted by each calculation unit based on the unit metal resource quantity conversion record, accumulating the resource quantity values item by item according to metal categories, merging and sorting accumulated results according to space indexes, and forming and outputting a multi-metal resource quantity value list.
- 10. A multi-metal ore resource quantity estimation system, characterized in that the multi-metal ore resource quantity estimation method according to any one of claims 1 to 9, the system comprising: The unit construction module is used for matching the drilling section with the spatial index based on boundary coordinates, drilling section grades, unit volumes and unit densities recorded by the ore body model file, and binding volume and density information to obtain an ore body calculation unit spatial sequence arranged according to the spatial index; The association constraint module is used for reading multi-metal grade information of adjacent units based on the ore body calculation unit space sequence, checking continuous intervals of trend and trend change directions and amplitude, eliminating abrupt metal items, registering symbiotic pairing and forming a multi-metal space association constraint structure; The grade combination module is used for screening metal items allowed to participate in processing in the calculation unit based on the multi-metal space association constraint structure, reading corresponding grade values, carrying out contribution sorting on the metal grades and completing choosing and sorting, and generating a combined metal grade value set of the calculation unit; The effective screening module is used for carrying out consistency check on the metal grade records and eliminating violations based on the multi-metal space association constraint structure and the combined metal grade value set, and carrying out gentle processing on the difference values of adjacent units to obtain an effective mineralization calculation unit set; And the resource summarizing module is used for reading the volume and the density of the computing units based on the effective mineralization computing unit set, completing conversion of paired metal grades, accumulating the conversion results according to metal types, merging according to spatial indexes and outputting a multi-metal resource quantity numerical list.
Description
Multi-metal ore resource quantity estimation method and system thereof Technical Field The invention relates to the technical field of mineral resource quantity calculation, in particular to a multi-metal mineral resource quantity estimation method and a system thereof. Background The technical field of mineral resource amount calculation particularly relates to a technology for carrying out mathematical modeling and numerical calculation on spatial distribution and grade parameters of various metal mineral bodies and occurrence boundaries based on computer digital processing in mining and geological engineering, and the technical field takes geometric modeling, grade assignment, volume calculation and resource amount summarization of the mineral bodies as cores, and realizes quantitative estimation on resource amounts of different metal elements in the multi-metal mineral deposit by carrying out calculation processing on drilling data, block model parameters and constraint conditions, thereby being mainly used for mineral exploration evaluation, mine design and mining planning stages. The method aims to realize collaborative estimation of the multi-metal resource quantity in the same ore body model, avoid boundary inconsistency and result deviation caused by separate calculation of single metal, and can obtain a plurality of metal resource quantity calculation results under the constraint condition of the unified ore body by the method, thereby improving the stability of the resource quantity estimation results in terms of space consistency and numerical value comparability. In the prior art, the calculation of the multi-metal ore resource quantity is mostly carried out by adopting an operation mode of single-metal independent modeling and calculation, although the statistics of the multi-metal resource quantity can be completed in practical application, the lack of joint constraint on spatial distribution trend of different metals leads to the lack of consistency control on the change direction and the change amplitude of different metal grades in the same calculation unit and adjacent blocks, when the distribution of drilling holes is uneven and the local mineralization condition is complex, the single-metal calculation result is easily affected by local abnormal grades and is directly reflected into the collection of the resource quantities, the fluctuation of the result is caused, the processing of the multi-metal grades in the unit in the prior art usually lacks systematic contribution distinction, different metals often participate in the calculation in a parallel mode, so that the superposition influence on the result is caused by the low-correlation and low-stability metals, in the multi-metal symbiotic ore deposit, when the certain metal grade is obviously jumped and is not effectively constrained, the problem of spatial concentration and discrete imbalance can occur to the result of the resource quantity, the judgment on the whole scale and the grade structure of the ore body is affected, and the application stability of the resource quantity in mine design and mining planning is reduced. Disclosure of Invention The invention aims to solve the defects in the prior art, and provides a multi-metal ore resource quantity estimation method and a system thereof. In order to achieve the purpose, the invention adopts the following technical scheme that the multi-metal ore resource quantity estimation method comprises the following steps: s1, reading a coordinate positioning unit based on boundary coordinates, drilling section grades, unit volumes and unit densities recorded by a ore body model file, matching drilling sections to unit indexes, and binding a volume density writing sequence to obtain an ore body calculation unit space sequence; S2, based on the ore body calculation unit space sequence, extracting multi-metal grades of adjacent units, comparing trend and trend change directions and judging continuous range intervals according to space correlation constraint requirements of a collaborative Kriging method, removing abrupt metal items, registering symbiotic pairing and forming a multi-metal space correlation constraint structure; S3, based on the multi-metal space association constraint structure, a screening unit allows metal items and reads corresponding grade values, a retention item is selected after contribution sorting is calculated based on the Squalman grade correlation coefficient, and the retention item is written into a unit record after low contribution items are removed, so that a combined metal grade value set is generated; S4, checking grades and constraints based on the multi-metal space association constraint structure and the combined metal grade value set, selecting reserved metals, eliminating violating items, comparing adjacent difference values and smoothing abnormal units to obtain an effective mineralization calculation unit set; And S5, reading the unit volume and the u