CN117057633B - Method and system for determining grade index of low-grade multi-metal ore - Google Patents
Method and system for determining grade index of low-grade multi-metal ore Download PDFInfo
- Publication number
- CN117057633B CN117057633B CN202311324867.8A CN202311324867A CN117057633B CN 117057633 B CN117057633 B CN 117057633B CN 202311324867 A CN202311324867 A CN 202311324867A CN 117057633 B CN117057633 B CN 117057633B
- Authority
- CN
- China
- Prior art keywords
- grade
- component
- ore
- comprehensive
- scheme
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 136
- 229910052751 metal Inorganic materials 0.000 title claims abstract description 42
- 239000002184 metal Substances 0.000 title claims abstract description 41
- 239000012141 concentrate Substances 0.000 claims abstract description 74
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 73
- 239000011707 mineral Substances 0.000 claims abstract description 73
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000011084 recovery Methods 0.000 claims description 34
- 239000007787 solid Substances 0.000 claims description 27
- 238000000605 extraction Methods 0.000 claims description 17
- 238000009826 distribution Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 16
- 238000000926 separation method Methods 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 14
- 238000003556 assay Methods 0.000 claims description 12
- 239000011435 rock Substances 0.000 claims description 12
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 230000015654 memory Effects 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 239000000470 constituent Substances 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 229910052797 bismuth Inorganic materials 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000002699 waste material Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method and a system for determining a grade index of low-grade multi-metal ore, comprising the following steps: s1, acquiring test data of a mineral exploration engineering sample of low-grade multi-metal ores, wherein the test data records the test grade of each component; s2, establishing a three-dimensional ore body model by using a radial basis function based on a mixed circle method and a dynamic comprehensive grade method according to the tested grades of the components, calculating to obtain a plurality of groups of resource quantity and geological grade of the ore body under different boundary grades and industrial grades, and calling each calculation result as a scheme; s3, establishing a comprehensive technical and economic analysis model for solving the total concentrate amount and the total profit of the ore body, and calculating the total concentrate amount and the total profit of the ore body; s4, taking the total concentrate amount and the total profit of the ore body as decision targets, and selecting an optimal scheme of the boundary grade and the industrial grade of the ore body by using a fuzzy comprehensive judgment method. The invention can utilize more mineral resources and reduce the loss of the mineral resources to the greatest possible extent.
Description
Technical Field
The invention relates to the technical field of low-grade multi-metal ore grade index determination, in particular to a method and a system for determining low-grade multi-metal ore grade index.
Background
The grade index of the defined ore body is an important index in the industrial index of the ore body, and directly influences the size of the ore body resource amount and the grade of the ore body. What is carried out in many mines in China is a boundary grade and industrial grade dual-grade index system. At present, common methods for determining boundary grade indexes of single metal ores include analogy methods and statistical methods, and common profit and loss balance methods for determining industrial grades. However, how to determine the boundary grade and the industrial grade of the low-grade multi-metal ore has more difficulties, and an effective method for determining the grade index of the low-grade multi-metal ore is needed.
Disclosure of Invention
The invention provides a method and a system for determining grade indexes of low-grade multi-metal ores, which are used for solving the problems in the prior art, and the technical scheme is as follows:
in one aspect, a method for determining a grade index of a low-grade multi-metal ore is provided, comprising:
s1, acquiring test data of a mineral exploration engineering sample of low-grade multi-metal ores, wherein the test data records the test grade of each component;
s2, establishing a three-dimensional ore body model by using a radial basis function based on a mixed circle method and a dynamic comprehensive grade method according to the tested grades of the components, calculating to obtain a plurality of groups of resource quantity and geological grade of the ore body under different boundary grades and industrial grades, and calling each calculation result as a scheme;
S3, establishing a comprehensive technical and economic analysis model for solving the total concentrate amount and the total profit of each scheme ore body, and calculating the total concentrate amount and the total profit of each scheme ore body;
s4, taking the total concentrate amount and the total profit of the ore body of each scheme as decision targets, and using a fuzzy comprehensive judgment method to preferably select the optimal scheme of the boundary grade and the industrial grade of the ore body.
Optionally, the mixed circle method in S2 refers to that the grade index of any one of the useful components reaches the index requirement, and is determined to meet the requirement of delineating the grade index of the ore body;
the dynamic comprehensive grade method in the S2 means that the grade index of the main component can not meet the requirement, but the grade of the co-associated component can meet the grade index requirement after being converted into the main component, and the grade index meets the requirement of the grade index of the delineating ore body; the converted grade is called as comprehensive grade; the dynamic comprehensive grade method means that the comprehensive grade changes along with the change of boundary grade and industrial grade;
the mixed circle method and the dynamic comprehensive grade method in the S2 refer to the condition that the grade index of any one useful component meets the index requirement and is determined to meet the requirement of the grade index of the defined ore body; or the grade index of the main component can not meet the requirement, but the grade of the co-associated component can meet the requirement of the grade index after being converted into the main component, and the grade index meets the requirement of the grade index of the delineating ore body.
Optionally, the S2 specifically includes:
s21, according to the assay grades of the components, counting and analyzing grade distribution rules of the components of the sample;
s22, setting a first group of boundary grades and industrial grades of each component on the basis of the grade distribution rule of each component of the sample, and carrying out sample length combination according to the first group of boundary grades, the minimum recoverable thickness and the stone-clamping rejecting thickness of each component;
s23, performing spatial interpolation on each component through RBF, and establishing a solid model of each component;
s24, dividing the solid model of each component into a plurality of small blocks, carrying out or relation graphic constraint on the block model by using the solid model of each component, and estimating the grade value of each component of each small mineral block based on a Ke-Lige method;
s25, selecting the jth mineral blocks with the grade value of any component i larger than the first group of industrial grade set by the component from the block model, calculating the resource quantity of the mineral blocks according to a formula 1, summing to obtain the resource quantity of the whole mineral body, and then calculating the geological grade of each component by using a formula 2;
(equation 1)
Wherein:-the volume of the jth ore block;
d—average density of ore;
q is the resource quantity of the whole ore body;
t is the total number of ore blocks;
(equation 2)
Wherein:-geological grade of component i;
-grade of j-th ore block of component i;
-the resource amount of the jth ore block;
q is the resource quantity of the whole ore body;
i, j-denote the constituent and the lump ore, respectively.
Optionally, the S2 further includes:
s26, dynamically calculating the comprehensive grade of the ore body, and carrying out sample length combination of the comprehensive components consisting of the main component and the co-associated components;
s27, performing spatial interpolation on the comprehensive components through RBF, and establishing a solid model of the comprehensive components;
s28, estimating the grade value of each component of each small mineral block based on a Ke-Lig method by using the or relation operation of the solid model of each component and the solid model of the comprehensive component as the constraint condition of the block model;
s29, screening j-th mineral blocks with the grade value of any component i larger than the first group of industrial grade given by the component from the block model, calculating the resource quantity of the mineral blocks according to a formula 1, summing to obtain the resource quantity of the whole mineral body, and then calculating the geological grade of each component by using a formula 2;
setting different boundary grades and industrial grades of each component, repeating the steps S22 to S29, calculating to obtain the resource quantity and geological grade of each component under each group of boundary grades and industrial grades, and calling each calculation result as a scheme.
Optionally, the step S26 of dynamically calculating the comprehensive grade of the ore body specifically includes:
s261, establishing a model for solving the mineral separation recovery rate by taking the input grade of each component as an independent variable according to mineral separation production daily report or mineral separation experimental data, wherein the model is shown in a formula 3:
(equation 3)
Wherein:beneficiation recovery of component i;
-a model for mineral separation recovery according to the grade of the beneficiation;
-grade of choice of component i;
s262, calculating the extraction grade by using a formula 4 according to the geological grade of each component under each group of boundary grade and industrial grade, and assuming that the ore of the concentrating mill is only from the mine with the boundary grade and industrial grade to be determined, considering the amount of the selected ore as the extraction ore amount and the extraction grade as the extraction grade;
(equation 4)
Wherein:-the produced grade of component i;
-geological grade of component i, obtained from equation 2;
h-rock mixing rate;
w is the grade of surrounding rock;
will beThe visual identity is +.>Substituting formula 3 to obtain->Then, the comprehensive grade is obtained by using the formula 5 and the formula 6>；
(equation 5)
(equation 6)
Wherein:-comprehensive grade;
-grade of the main component;
n-the amount of co-associated components available;
-the grade conversion factor of the co-associated component i;
-grade of co-associated component i;
mineral recovery of co-associated component i;
-the ton concentrate price of the co-associated component i;
-concentrate grade of the main component;
mineral recovery of the main component;
-the ton concentrate price of the main component;
the concentrate grade of component i is co-associated.
Optionally, the step S3 specifically includes:
s31, calculating the amount of the extracted ore by using a formula 7;
(equation 7)
Wherein:-the k-th option produces ore quantities;
-a kth scheme resource amount;
l-loss rate;
h-rock mixing rate;
s32, calculating the concentrate quantity of the j-th scheme of the i component by using a formula 8, and calculating the concentrate total quantity of each scheme;
l-loss rate;
h-rock mixing rate;
s32, calculating the concentrate quantity of the j-th scheme of the i component by using a formula 8, and calculating the concentrate total quantity of each scheme;
(equation 8)
Wherein:-i component concentrate quantity of the kth variant;
mineral recovery of the kth variant of component i,/->The method is characterized in that i component is a model for solving the ore dressing recovery rate according to the input grade;
-the amount of ore selected in the kth regime, as the amount of ore mined;
-i component k scheme grade, as produced grade;
summing the concentrate amounts of all components of the same scheme to obtain the concentrate total amount of each scheme ；
S33, calculating the income of each component of each scheme according to the unit price and concentrate quantity of each component, and the formula 9;
(equation 9)
Wherein:revenue for the kth scenario of the i component;
-unit price of component i;
-i component concentrate quantity of the kth variant;
summing the revenues of all the components of the same scheme to obtain the total revenues of each scheme；
S34, calculating the total cost of each scheme;
according to the manufacturing cost method, the operation cost of each component is calculated firstly, and then other costs are added to obtain the total cost;
s35, calculating the total profit of each scheme
The formula for calculating the total profit: total profit = total revenue-total cost-tax sales and additional.
Optionally, the S4 specifically includes:
s41, defining a membership function shown in a formula 10:
(equation 10)
Wherein:-degree of membership of the kth scheme to the p-th target, fuzzy degree of membership;
X kp -a calculation of a p-th decision target under a kth scheme;
X pmin -the minimum of the p-th decision target in all schemes;
X pmax -maximum value of p-th decision target in all schemes;
the decision targets comprise the total concentrate amount and the total profit of ore bodies;
s42, calculating the comprehensive membership degree;
after the membership degree of each decision target value of different schemes is calculated by using the membership function, the weighted average value of each decision target value is calculated, and for the kth specific scheme, the comprehensive membership degree is calculated by using a formula 11:
(equation 11)
Wherein:-the comprehensive membership of the kth scheme;
-the weight of the p-th object;
q-number of decision targets;
the scheme with the largest comprehensive membership in all schemes is the optimal scheme of the optimized boundary grade and industrial grade.
In another aspect, a system for determining a grade indicator of a low grade multi-metal ore is provided, the system comprising:
the acquisition module is used for acquiring the test data of the exploring engineering sample of the low-grade multi-metal ore, and the test data records the test grade of each component;
the first calculation module is used for establishing a three-dimensional ore body model by using a radial basis function based on a mixed circle method and a dynamic comprehensive grade method according to the assay grade of each component, calculating to obtain the resource quantity and the geological grade of the ore body under a plurality of groups of different boundary grades and industrial grades, and referring each calculation result as a scheme;
the second calculation module is used for establishing a comprehensive technical and economic analysis model for solving the total concentrate amount and the total profit of each scheme ore body and calculating the total concentrate amount and the total profit of each scheme ore body;
and the optimization module is used for taking the total concentrate amount and the total profit of the ore body of each scheme as decision targets, and optimizing the optimal scheme of the boundary grade and the industrial grade of the ore body by using a fuzzy comprehensive judgment method.
In another aspect, an electronic device is provided, the electronic device including a processor and a memory having instructions stored therein, the instructions being loaded and executed by the processor to implement the above-described method of determining a grade indicator of a low-grade multi-metal mine.
In another aspect, a computer readable storage medium having instructions stored therein is provided, the instructions being loaded and executed by a processor to implement the above-described method of determining a grade index of a low-grade multi-metal mine.
The technical scheme provided by the invention has the beneficial effects that at least:
(1) When the boundary grade and the industrial grade are determined from a plurality of schemes, the invention considers the correlation among technical economy such as boundary grade, industrial grade, geological grade, extraction grade, selection grade, ore dressing recovery rate and the like, and achieves the linkage of other technical indexes in the range of changing the boundary grade and the industrial grade.
(2) The invention has similar place with the existing 'grade index integral dynamic optimization method', but the superior place is reflected on the model of ore body resource quantity and geological grade, the invention automatically delineating ore body through RBF, when the boundary grade and industrial grade change, the resource quantity and geological grade of the ore body are obtained, the invention can also provide three-dimensional ore body models under different boundary grades and industrial grades.
(3) Compared with the existing mixed-circle method and the comprehensive grade method, the low-grade multi-metal ore circle ore body fixing method can utilize more mineral resources, and the loss of the mineral resources is reduced to the greatest extent.
(4) When the grade conversion coefficient is calculated, the internal relation between the input grade and the beneficiation recovery rate is considered, and a model for solving the beneficiation recovery rate according to the input grade is established. When the boundary grade and the industrial grade change to cause the change of the geological grade, the extraction grade and the selection grade, the comprehensive grade can dynamically change.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for determining a grade index of a low-grade multi-metal ore according to an embodiment of the present invention;
FIG. 2 is a flow chart for estimating the resource quantity and the geological grade based on a mixed circle method and a dynamic comprehensive grade method provided by the embodiment of the invention;
FIG. 3 is a block diagram of a system for determining a grade index of a low-grade multi-metal ore according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the technical problems, technical solutions and advantages to be solved more apparent, the following detailed description will be given with reference to the accompanying drawings and specific embodiments.
As shown in fig. 1, the embodiment of the invention provides a method for determining a grade index of low-grade multi-metal ore, which comprises the following steps:
s1, acquiring test data of a mineral exploration engineering sample of low-grade multi-metal ores, wherein the test data records the test grade of each component;
s2, establishing a three-dimensional ore body model by using a radial basis function based on a mixed circle method and a dynamic comprehensive grade method according to the tested grades of the components, calculating to obtain a plurality of groups of resource quantity and geological grade of the ore body under different boundary grades and industrial grades, and calling each calculation result as a scheme;
s3, establishing a comprehensive technical and economic analysis model for solving the total concentrate amount and the total profit of each scheme ore body, and calculating the total concentrate amount and the total profit of each scheme ore body;
s4, taking the total concentrate amount and the total profit of the ore body of each scheme as decision targets, and using a fuzzy comprehensive judgment method to preferably select the optimal scheme of the boundary grade and the industrial grade of the ore body.
The method for determining the grade index of the low-grade multi-metal ore provided by the embodiment of the invention is described in detail below with reference to fig. 2 (solid line in the figure is based on a mixed circle method, and dotted line is based on a dynamic comprehensive grade method), and comprises the following steps:
s1, acquiring test data of a mineral exploration engineering sample of low-grade multi-metal ores, wherein the test data records the test grade of each component;
specifically, according to the data obtained by mineral exploration engineering such as drilling, pit exploration and the like, a geological database is built in three-dimensional mining software, and the database comprises open-pore coordinates, inclinometry data, lithology data and test data. In the assay data, a row is reserved for recording the comprehensive grade except for recording the assay grade of each component.
S2, establishing a three-dimensional ore body model by using a radial basis function based on a mixed circle method and a dynamic comprehensive grade method according to the tested grades of the components, calculating to obtain a plurality of groups of resource quantity and geological grade of the ore body under different boundary grades and industrial grades, and calling each calculation result as a scheme;
optionally, the mixed circle method in S2 refers to that the grade index of any one of the useful components reaches the index requirement, and is determined to meet the requirement of delineating the grade index of the ore body;
The dynamic comprehensive grade method in the S2 means that the grade index of the main component can not meet the requirement, but the grade of the co-associated component can meet the grade index requirement after being converted into the main component, and the grade index meets the requirement of the grade index of the delineating ore body; the converted grade is called as comprehensive grade; the dynamic comprehensive grade method means that the comprehensive grade changes along with the change of boundary grade and industrial grade;
the mixed circle method and the dynamic comprehensive grade method in the S2 refer to the condition that the grade index of any one useful component meets the index requirement and is determined to meet the requirement of the grade index of the defined ore body; or the grade index of the main component can not meet the requirement, but the grade of the co-associated component can meet the requirement of the grade index after being converted into the main component, and the grade index meets the requirement of the grade index of the delineating ore body.
If only the main component ring is adopted to fix ore body, the resources which can meet the requirements of the grade index of the main component are lost, and the grade index of the co-associated component is converted into the main component although the grade index of the main component can not meet the requirements. The comprehensive grade is used for delineating the ore body according to the grade index of the main component, and obviously, the resource quantity of the ore body delineating the comprehensive grade is larger than that of the ore body delineating the main component.
However, the existing comprehensive grade method calculates the comprehensive gradeIn particular in calculating the grade conversion coefficient of the co-associated component i>The beneficiation recovery of co-associated component i was considered to be constant.
(equation 5)
(equation 6)
Wherein:-comprehensive grade;
-grade of the main component;
n-the amount of co-associated components available;
-the grade conversion factor of the co-associated component i;
-grade of co-associated component i;
mineral recovery of co-associated component i;
-the ton concentrate price of the co-associated component i;
-concentrate grade of the main component;
mineral recovery of the main component;
-the ton concentrate price of the main component;
the concentrate grade of component i is co-associated.
When the prior comprehensive grade method does not consider the change of the selected grade, the ore dressing recovery rate can change correspondingly, so the following dynamic relationship cannot be reflected: when the boundary grade and the industrial grade change, the geological grade and the extraction grade change, the selected grade changes, and the selected recovery rate changes. That is, the comprehensive grade should be able to change along with the change of the boundary grade and the industrial grade, and should not be static, so the embodiment of the invention adopts a dynamic comprehensive grade method, and the calculation formula of the dynamic comprehensive grade method will be described below.
In addition, the problems exist in the single comprehensive grade method or the mixed circle method: ore body defined by comprehensive grade methodThe ore body defined by the mixed circle method is not necessarily included, and the ore body defined by the comprehensive grade is not necessarily included, so that the resource amount can be lost when the ore body is defined by either method, and the method is described below by way of example. Table 1 shows a certain WO 3 Industrial indexes of the multi-metal ore ring ore body of Mo and Bi.
TABLE 1 Industrial index of round ore body
WO is incorporated into 3 As the main component, mo and Bi are used as co-associated components, after grade conversion through a formula 5 and a formula 6, grade conversion coefficients of Mo and Bi are respectively 1.058 and 0.1658, and the resource amount defined by the existing mixed-loop method, comprehensive grade method and mixed-loop method and dynamic comprehensive grade method provided by the embodiment of the invention is shown in a table 2. As can be seen from the table, the resource quantity circled by the ore mixing circle method is larger than the resource quantity circled by the comprehensive grade method, but smaller than the resource quantity circled by the mixing circle method and the dynamic comprehensive grade method. After the ore body is defined by adopting the mixed circle method and the dynamic comprehensive grade method, the volume of the ore body is increased, and the metal quantity is increased, so that the ore body defined by the mixed circle method does not completely cover the ore body defined by the comprehensive grade method, and the mixed circle method possibly loses the part of resources of which the single metal grade index does not meet the grade index, but the comprehensive grade can meet the grade index requirement. For low-grade multi-metal ores, the amount of resources defined by the comprehensive grade method may be more than that of the mixed-grade method, but may be less than that of the mixed-grade method, depending on the grade distribution of multi-metal elements in the sample. The embodiment of the invention adopts a dynamic comprehensive grade method, which can possibly lead the WO 3 Samples with a grade lower than 0.1% are converted to a comprehensive grade higher than 0.1% by Mo and Bi, such as sample No. 1 in Table 3, but may not meet the grade requirement in the mixed-loop method, such as sample No. 2,3 and 4 in Table 3, and the mixed-loop method marks the 3 samples asore, however, the integrated grade method is marked as waste because it is lower than the requirement of 0.1% of the main component, so that ore bodies defined by the integrated grade method may be smaller than those defined by the mixed-grade method.
TABLE 2 Mixed circle method+comprehensive grade method circle defined resource quantity
TABLE 3 conversion by dynamic Integrated grade method
The embodiment of the invention provides a method for determining and delineating grade indexes of low-grade multi-metal ores by using a Radial Basis Function (RBF) to establish a three-dimensional ore body model based on a mixed circle method and a dynamic comprehensive grade method and taking comprehensive benefits (including resource benefits and economic benefits) as decision targets.
Optionally, the S2 specifically includes:
s21, according to the assay grades of the components, counting and analyzing grade distribution rules of the components of the sample;
specifically, the embodiment of the invention draws a grade distribution histogram of each component and analyzes the grade distribution rule of each component. If the grade distribution of each component does not follow normal distribution, taking the logarithm to basically follow the lognormal distribution.
Analyzing whether an abnormal value exists in the sample and processing the abnormal value. Calculating the average value of the grades of all the components, and setting the average value to be more than 6 times as suspicious abnormal value. Looking at the spatial distribution of the suspected outliers, if samples of the suspected outliers are scattered at various locations in the sample, indicating that there are outliers in the sample, replacing the outliers with an average value. If samples of suspected outliers are clustered together, it is indicated that these samples are from rich sites and cannot be rejected as outliers.
S22, setting a first group of boundary grades and industrial grades of each component on the basis of the grade distribution rule of each component of the sample, and carrying out sample length combination according to the first group of boundary grades, the minimum recoverable thickness and the stone-clamping rejecting thickness of each component;
according to the embodiment of the invention, on the basis of the grade distribution rule of each component of the sample, a first group of boundary grade and industrial grade of each component are set, sample length combination is carried out according to the requirements of industrial indexes such as the first group of boundary grade, the minimum recoverable thickness, the stone-clamping rejection thickness and the like of each component, an ore sample section is marked as ore, and a waste stone section is marked as waste.
S23, performing spatial interpolation on each component through RBF, and establishing a solid model of each component;
Specifically, the step S23 specifically includes:
s231, setting a range for delineating the ore body;
and setting the range of the delineating the ore body by setting the range of x, y and z coordinates of the delineating the ore body according to the range of the exploration engineering control.
S232, modeling through a variation function in three-dimensional mining software surfac, and determining the lengths of a principal axis, a semi-principal axis and a secondary principal axis of a search ellipsoid of each component by referring to the shape of a ore body; simulating the value of nugget, sill, range parameters of the spherical model;
s233, utilizing the RBF interptant function of the three-dimensional modeling software to carry out spatial interpolation on each component through RBF.
And (3) during interpolation, using a combined sample after sample length combination, and performing spatial interpolation in a set range by utilizing a mineral sample section marked as ore according to an ellipsoid search parameter and an ellipsoid model to form a solid model of each component.
S24, dividing the solid model of each component into a plurality of small blocks, carrying out or relation graphic constraint on the block model by using the solid model of each component, and estimating the grade value of each component of each small mineral block based on a Ke-Lige method;
specifically, the embodiment of the invention divides the whole solid model of each component into a plurality of small blocks, the size of the blocks in the horizontal direction is generally 1/3-1/5 of the exploration interval, and the size of the blocks in the vertical direction is 2-3 times of the combined sample length. The properties of each component are added to the block model.
And (3) carrying out or relation graphic constraint on the block model by using the solid model of each component, estimating the grade value of each component of each small mineral block based on a Kriging method, and storing the grade value into the corresponding attribute. The values of the ellipsoid parameter and the spherical model parameter are the same as the requirements in S232.
S25, selecting the jth mineral blocks with the grade value of any component i larger than the first group of industrial grade set by the component from the block model, calculating the resource quantity of the mineral blocks according to a formula 1, summing to obtain the resource quantity of the whole mineral body, and then calculating the geological grade of each component by using a formula 2;
(equation 1)
Wherein:-the volume of the jth ore block;
d—average density of ore;
q is the resource quantity of the whole ore body;
t is the total number of ore blocks;
(equation 2)
Wherein:-geological grade of component i;
-grade of j-th ore block of component i;
-the resource amount of the jth ore block;
q is the resource quantity of the whole ore body;
i, j-denote the constituent and the lump ore, respectively.
Optionally, the S2 further includes:
s26, dynamically calculating the comprehensive grade of the ore body, and carrying out sample length combination of the comprehensive components consisting of the main component and the co-associated components;
s27, performing spatial interpolation on the comprehensive components through RBF, and establishing a solid model of the comprehensive components;
S28, estimating the grade value of each component of each small mineral block based on a Ke-Lig method by using the or relation operation of the solid model of each component and the solid model of the comprehensive component as the constraint condition of the block model;
s29, screening j-th mineral blocks with the grade value of any component i larger than the first group of industrial grade given by the component from the block model, calculating the resource quantity of the mineral blocks according to a formula 1, summing to obtain the resource quantity of the whole mineral body, and then calculating the geological grade of each component by using a formula 2;
setting different boundary grades and industrial grades of each component, repeating the steps S22 to S29, calculating to obtain the resource quantity and geological grade of each component under each group of boundary grades and industrial grades, and calling each calculation result as a scheme.
Optionally, the step S26 of dynamically calculating the comprehensive grade of the ore body specifically includes:
s261, establishing a model for solving the mineral separation recovery rate by taking the input grade of each component as an independent variable according to mineral separation production daily report or mineral separation experimental data, wherein the model is shown in a formula 3:
(equation 3)
Wherein:beneficiation recovery of component i;
-a model for mineral separation recovery according to the grade of the beneficiation;
-grade of choice of component i;
S262, calculating the extraction grade by using a formula 4 according to the geological grade of each component under each group of boundary grade and industrial grade, and assuming that the ore of the concentrating mill is only from the mine with the boundary grade and industrial grade to be determined, considering the amount of the selected ore as the extraction ore amount and the extraction grade as the extraction grade;
(equation 4)
Wherein:-the produced grade of component i;
-geological grade of component i, obtained from equation 2;
h-rock mixing rate;
w is the grade of surrounding rock;
will beThe visual identity is +.>Substituting formula 3 to obtain->Then, the comprehensive grade is obtained by using the formula 5 and the formula 6>；
(equation 5)
(equation 6)
Wherein:-comprehensive grade;
-grade of the main component;
n-the amount of co-associated components available;
-the grade conversion factor of the co-associated component i;
-grade of co-associated component i;
mineral recovery of co-associated component i;
-the ton concentrate price of the co-associated component i;
-concentrate grade of the main component;
mineral recovery of the main component; />
-the ton concentrate price of the main component;
-co-associated componentsi concentrate grade.
According to the embodiment of the invention, the resource quantity and the geological grade of the ore body can be obtained under a plurality of groups of different boundary grades and industrial grades, the data in the form of table 4 are obtained (the total concentrate quantity and the total profit are required to be calculated in the follow-up process), and each calculation result is called as a scheme.
TABLE 4 resource quantity and geological grade under different scenarios
S3, establishing a comprehensive technical and economic analysis model for solving the total concentrate amount and the total profit of each scheme ore body, and calculating the total concentrate amount and the total profit of each scheme ore body;
the invention uses Excel or some programming language such as Python, uses the data in table 4 to calculate the ore extraction amount, the extraction grade of each component, the concentrate amount of each component, the cost and income of each component, etc. under different schemes in table 4 through the following formulas and models, and finally calculates the total concentrate amount and the total profit.
Optionally, the step S3 specifically includes:
s31, calculating the amount of the extracted ore by using a formula 7;
(equation 7)
Wherein:-the k-th option produces ore quantities;
-a kth scheme resource amount;
l-loss rate;
h-rock mixing rate;
s32, calculating the concentrate quantity of the j-th scheme of the i component by using a formula 8, and calculating the concentrate total quantity of each scheme;
(equation 8)
Wherein:-i component concentrate quantity of the kth variant;
mineral recovery of the kth variant of component i,/->The method is characterized in that i component is a model for solving the ore dressing recovery rate according to the input grade;
-the amount of ore selected in the kth regime, as the amount of ore mined; / >
-i component k scheme grade, as produced grade;
summing the concentrate amounts of all components of the same scheme to obtain the concentrate total amount of each scheme；
S33, calculating the income of each component of each scheme according to the unit price and concentrate quantity of each component, and the formula 9;
(equation 9)
Wherein:revenue for the kth scenario of the i component;
-unit price of component i;
-i component concentrate quantity of the kth variant;
summing the revenues of all the components of the same scheme to obtain the total revenues of each scheme；
S34, calculating the total cost of each scheme;
according to the manufacturing cost method, the operation cost of each component is calculated firstly, and then other costs are added to obtain the total cost;
the common process, such as crushing process in ore dressing, exists in the multi-metal ore treatment process, and the common process cost is shared. In view of the principle that the common flow is required as it is when only the main component is used even if the common accompanying component is not used, the cost of the common flow in calculation is completely borne by the main component.
S35, calculating the total profit of each scheme
The formula for calculating the total profit: total profit = total revenue-total cost-tax sales and additional.
S4, taking the total concentrate amount and the total profit of the ore body of each scheme as decision targets, and using a fuzzy comprehensive judgment method to preferably select the optimal scheme of the boundary grade and the industrial grade of the ore body.
By means of the first and second sections, a number of different boundary and industrial grade solutions can be created, among which the best solution can be preferred by means of comprehensive evaluation, since the total profit maximum and the concentrate total maximum do not necessarily occur in the same solution. The embodiment of the invention realizes the scheme optimization by using a fuzzy comprehensive judgment method.
Optionally, the S4 specifically includes:
s41, defining a membership function shown in a formula 10:
(equation 10)
Wherein:-degree of membership of the kth scheme to the p-th target, fuzzy degree of membership;
X kp -a calculation of a p-th decision target under a kth scheme;
X pmin -the minimum of the p-th decision target in all schemes;
X pmax -maximum value of p-th decision target in all schemes;
the decision targets comprise the total concentrate amount and the total profit of ore bodies;
s42, calculating the comprehensive membership degree;
after the membership degree of each decision target value of different schemes is calculated by using the membership function, the weighted average value of each decision target value is calculated, and for the kth specific scheme, the comprehensive membership degree is calculated by using a formula 11:
(equation 11)
Wherein:-the comprehensive membership of the kth scheme;
-the weight of the p-th object;
q-number of decision targets;
the scheme with the largest comprehensive membership in all schemes is the optimal scheme of the optimized boundary grade and industrial grade.
As shown in fig. 3, the embodiment of the invention further provides a system for determining a grade index of low-grade multi-metal ore, which comprises:
an acquisition module 310, configured to acquire assay data of a prospecting engineering sample of the low-grade multi-metal ore, where the assay data records an assay grade of each component;
the first calculation module 320 is configured to establish a three-dimensional ore body model by using a radial basis function based on a mixed circle method and a dynamic comprehensive grade method according to the assay grade of each component, calculate to obtain a plurality of groups of resource amounts and geological grades of the ore body under different boundary grades and industrial grades, and refer to each calculation result as a scheme;
a second calculation module 330 for establishing a comprehensive technical economic analysis model for obtaining the total concentrate amount and the total profit of each scheme ore body, and calculating the total concentrate amount and the total profit of each scheme ore body;
a preference module 340 is used for taking the total concentrate amount and the total profit of the ore body of each scheme as decision targets, and using a fuzzy comprehensive judgment method, the best scheme of the boundary grade and the industrial grade of the ore body is preferred.
The functional structure of the system for determining the grade index of the low-grade multi-metal ore provided by the embodiment of the invention corresponds to the method for determining the grade index of the low-grade multi-metal ore provided by the embodiment of the invention, and is not repeated here.
Fig. 4 is a schematic structural diagram of an electronic device 400 according to an embodiment of the present invention, where the electronic device 400 may have a relatively large difference due to different configurations or performances, and may include one or more processors (central processing units, CPU) 401 and one or more memories 402, where the memories 402 store instructions, and the instructions are loaded and executed by the processors 401 to implement the steps of the method for determining a grade index of a low-grade multi-metal mine.
In an exemplary embodiment, a computer readable storage medium, such as a memory comprising instructions executable by a processor in a terminal to perform the above-described method of determining a grade indicator of a low-grade multi-metal mine, is also provided. For example, the computer readable storage medium may be ROM, random Access Memory (RAM), CD-ROM, magnetic tape, floppy disk, optical data storage device, etc.
It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program for instructing relevant hardware, where the program may be stored in a computer readable storage medium, and the storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.
The foregoing description of the preferred embodiments of the invention is not intended to limit the invention to the precise form disclosed, and any such modifications, equivalents, and alternatives falling within the spirit and scope of the invention are intended to be included within the scope of the invention.
Claims (6)
1. A method of determining a grade indicator for a low grade multi-metal ore, the method comprising:
s1, acquiring test data of a mineral exploration engineering sample of low-grade multi-metal ores, wherein the test data records the test grade of each component;
s2, establishing a three-dimensional ore body model by using a radial basis function based on a mixed circle method and a dynamic comprehensive grade method according to the tested grades of the components, calculating to obtain a plurality of groups of resource quantity and geological grade of the ore body under different boundary grades and industrial grades, and calling each calculation result as a scheme;
S3, establishing a comprehensive technical and economic analysis model for solving the total concentrate amount and the total profit of each scheme ore body, and calculating the total concentrate amount and the total profit of each scheme ore body;
s4, taking the total concentrate amount and the total profit of the ore body of each scheme as decision targets, and using a fuzzy comprehensive judgment method to preferably select the optimal scheme of the boundary grade and the industrial grade of the ore body;
the mixed circle method in S2 refers to that the grade index of any one useful component meets the index requirement, and is determined to meet the requirement of the grade index of the defined ore body;
the dynamic comprehensive grade method in the S2 means that the grade index of the main component can not meet the requirement, but the grade of the co-associated component can meet the grade index requirement after being converted into the main component, and the grade index meets the requirement of the grade index of the delineating ore body; the converted grade is called as comprehensive grade; the dynamic comprehensive grade method means that the comprehensive grade changes along with the change of boundary grade and industrial grade;
the mixed circle method and the dynamic comprehensive grade method in the S2 refer to the condition that the grade index of any one useful component meets the index requirement and is determined to meet the requirement of the grade index of the defined ore body; or the grade index of the main component can not meet the requirement, but the grade of the co-associated component can meet the grade index requirement after being converted into the main component, and the grade index meets the requirement of the grade index of the delineating ore body;
The step S2 specifically comprises the following steps:
s21, according to the assay grades of the components, counting and analyzing grade distribution rules of the components of the sample;
s22, setting a first group of boundary grades and industrial grades of each component on the basis of the grade distribution rule of each component of the sample, and carrying out sample length combination according to the first group of boundary grades, the minimum recoverable thickness and the stone-clamping rejecting thickness of each component;
s23, performing spatial interpolation on each component through RBF, and establishing a solid model of each component;
s24, dividing the solid model of each component into a plurality of small blocks, carrying out or relation graphic constraint on the block model by using the solid model of each component, and estimating the grade value of each component of each small mineral block based on a Ke-Lige method;
s25, selecting the jth mineral blocks with the grade value of any component i larger than the first group of industrial grade set by the component from the block model, calculating the resource quantity of the mineral blocks according to a formula 1, summing to obtain the resource quantity of the whole mineral body, and then calculating the geological grade of each component by using a formula 2;
(equation 1)
Wherein:-j-th ore blockIs defined by the volume of (2);
d—average density of ore;
q is the resource quantity of the whole ore body;
t is the total number of ore blocks;
(equation 2)
Wherein:-geological grade of component i;
-grade of j-th ore block of component i;
-the resource amount of the jth ore block;
q is the resource quantity of the whole ore body;
i, j-represent the constituent and the ore block, respectively;
the S2 further includes:
s26, dynamically calculating the comprehensive grade of the ore body, and carrying out sample length combination of the comprehensive components consisting of the main component and the co-associated components;
s27, performing spatial interpolation on the comprehensive components through RBF, and establishing a solid model of the comprehensive components;
s28, estimating the grade value of each component of each small mineral block based on a Ke-Lig method by using the or relation operation of the solid model of each component and the solid model of the comprehensive component as the constraint condition of the block model;
s29, screening j-th mineral blocks with the grade value of any component i larger than the first group of industrial grade given by the component from the block model, calculating the resource quantity of the mineral blocks according to a formula 1, summing to obtain the resource quantity of the whole mineral body, and then calculating the geological grade of each component by using a formula 2;
setting different boundary grades and industrial grades of each component, repeating the steps S22 to S29, calculating to obtain the resource quantity and geological grade of each component under each group of boundary grades and industrial grades, and calling each calculation result as a scheme;
And in the step S26, dynamically calculating the comprehensive grade of the ore body, wherein the method specifically comprises the following steps:
s261, establishing a model for solving the mineral separation recovery rate by taking the input grade of each component as an independent variable according to mineral separation production daily report or mineral separation experimental data, wherein the model is shown in a formula 3:
(equation 3)
Wherein:beneficiation recovery of component i;
-a model for mineral separation recovery according to the grade of the beneficiation;
-grade of choice of component i;
s262, calculating the extraction grade by using a formula 4 according to the geological grade of each component under each group of boundary grade and industrial grade, and assuming that the ore of the concentrating mill is only from the mine with the boundary grade and industrial grade to be determined, considering the amount of the selected ore as the extraction ore amount and the extraction grade as the extraction grade;
(equation 4)
Wherein:-the produced grade of component i;
-geological grade of component i, obtained from equation 2;
h-rock mixing rate;
w is the grade of surrounding rock;
will beThe visual identity is +.>Substituting formula 3 to obtain->Then, the comprehensive grade is obtained by using the formula 5 and the formula 6>；
(equation 5)
(equation 6)
Wherein:-comprehensive grade;
-grade of the main component;
n-the amount of co-associated components available;
-the grade conversion factor of the co-associated component i;
-grade of co-associated component i;
mineral recovery of co-associated component i;
-the ton concentrate price of the co-associated component i;
-concentrate grade of the main component;
mineral recovery of the main component;
-the ton concentrate price of the main component;
the concentrate grade of component i is co-associated.
2. The method according to claim 1, wherein S3 specifically comprises:
s31, calculating the amount of the extracted ore by using a formula 7;
(equation 7)
Wherein:-the k-th option produces ore quantities;
-a kth scheme resource amount;
l-loss rate;
h-rock mixing rate;
s32, calculating the concentrate quantity of the j-th scheme of the i component by using a formula 8, and calculating the concentrate total quantity of each scheme;
(equation 8)
Wherein:-i component concentrate quantity of the kth variant;
mineral recovery of the kth variant of component i,/->The method is characterized in that i component is a model for solving the ore dressing recovery rate according to the input grade;
-the amount of ore selected in the kth regime, as the amount of ore mined;
-i component k scheme grade, as produced grade;
summing the concentrate amounts of all components of the same scheme to obtain the concentrate total amount of each scheme；
S33, calculating the income of each component of each scheme according to the unit price and concentrate quantity of each component, and the formula 9;
(equation 9)
Wherein:revenue for the kth scenario of the i component;
-unit price of component i;
-i component concentrate quantity of the kth variant;
summing the revenues of all the components of the same scheme to obtain the total revenues of each scheme；
S34, calculating the total cost of each scheme;
according to the manufacturing cost method, the operation cost of each component is calculated firstly, and then other costs are added to obtain the total cost;
s35, calculating the total profit of each scheme
The formula for calculating the total profit: total profit = total revenue-total cost-tax sales and additional.
3. The method according to claim 1, wherein S4 specifically comprises:
s41, defining a membership function shown in a formula 10:
(equation 10)
Wherein:-degree of membership of the kth scheme to the p-th target, fuzzy degree of membership;
X kp -a calculation of a p-th decision target under a kth scheme;
X pmin -the minimum of the p-th decision target in all schemes;
X pmax -maximum value of p-th decision target in all schemes;
the decision targets comprise the total concentrate amount and the total profit of ore bodies;
s42, calculating the comprehensive membership degree;
after the membership degree of each decision target value of different schemes is calculated by using the membership function, the weighted average value of each decision target value is calculated, and for the kth specific scheme, the comprehensive membership degree is calculated by using a formula 11:
（11）
Wherein:-the comprehensive membership of the kth scheme;
-the weight of the p-th object;
q-number of decision targets;
the scheme with the largest comprehensive membership in all schemes is the optimal scheme of the optimized boundary grade and industrial grade.
4. A system for determining a grade indicator for a low grade multi-metal ore, the system comprising:
the acquisition module is used for acquiring the test data of the exploring engineering sample of the low-grade multi-metal ore, and the test data records the test grade of each component;
the first calculation module is used for establishing a three-dimensional ore body model by using a radial basis function based on a mixed circle method and a dynamic comprehensive grade method according to the assay grade of each component, calculating to obtain the resource quantity and the geological grade of the ore body under a plurality of groups of different boundary grades and industrial grades, and referring each calculation result as a scheme;
the second calculation module is used for establishing a comprehensive technical and economic analysis model for solving the total concentrate amount and the total profit of each scheme ore body and calculating the total concentrate amount and the total profit of each scheme ore body;
the optimization module is used for taking the total concentrate amount and the total profit of the ore body of each scheme as decision targets, and optimizing an optimal scheme of the boundary grade and the industrial grade of the ore body by using a fuzzy comprehensive judgment method;
The circle mixing method refers to that the grade index of any one useful component meets the index requirement, and is determined to meet the requirement of the grade index of the defined ore body;
the dynamic comprehensive grade method is characterized in that although the grade index of the main component can not meet the requirement, the grade of the co-associated component can meet the grade index requirement after being converted into the main component, and the grade index meets the requirement of the delineating ore body grade index; the converted grade is called as comprehensive grade; the dynamic comprehensive grade method means that the comprehensive grade changes along with the change of boundary grade and industrial grade;
the mixed circle method and the dynamic comprehensive grade method refer to the condition that the grade index of any one useful component meets the index requirement and is determined to meet the requirement of the grade index of the defined ore body; or the grade index of the main component can not meet the requirement, but the grade of the co-associated component can meet the grade index requirement after being converted into the main component, and the grade index meets the requirement of the grade index of the delineating ore body;
the first computing module is specifically configured to:
according to the assay grade of each component, counting and analyzing the grade distribution rule of each component of the sample;
setting a first group of boundary grades and industrial grades of each component on the basis of the grade distribution rule of each component of the sample, and carrying out sample length combination according to the first group of boundary grades, the minimum recoverable thickness and the stone-clamping removing thickness of each component;
Performing spatial interpolation on each component through RBF, and establishing a solid model of each component;
dividing the solid model of each component into a plurality of small blocks, carrying out or relation graphic constraint on the block model by using the solid model of each component, and estimating the grade value of each component of each small mineral block based on a Kriging method;
selecting the jth mineral blocks with the grade value of any component i larger than the first group of industrial grade set by the component from the block model, calculating the resource quantity of the mineral blocks according to a formula 1, summing to obtain the resource quantity of the whole mineral body, and then calculating the geological grade of each component by using a formula 2;
(equation 1)
Wherein:-the volume of the jth ore block;
d—average density of ore;
q is the resource quantity of the whole ore body;
t is the total number of ore blocks;
(equation 2)
Wherein:-geological grade of component i;
-grade of j-th ore block of component i;
-the resource amount of the jth ore block;
q is the resource quantity of the whole ore body;
i, j-represent the constituent and the ore block, respectively;
dynamically calculating the comprehensive grade of the ore body, and carrying out sample length combination of the comprehensive components consisting of the main component and the co-associated components;
performing spatial interpolation on the comprehensive components through RBF, and establishing a solid model of the comprehensive components;
The method comprises the steps of estimating the grade value of each component of each small mineral block based on a Kriging method by using the or relation operation of the solid model of each component and the solid model of the comprehensive component as the constraint condition of a block model;
selecting the jth mineral blocks with the grade value of any component i larger than the first group of industrial grade given by the component from the block model, calculating the resource quantity of the mineral blocks according to a formula 1, summing to obtain the resource quantity of the whole mineral body, and then calculating the geological grade of each component by using a formula 2;
setting different boundary grades and industrial grades of each component, repeating the steps S22 to S29, calculating to obtain the resource quantity and geological grade of each component under each group of boundary grades and industrial grades, and calling each calculation result as a scheme;
the dynamic calculation of the comprehensive grade of the ore body specifically comprises the following steps:
according to mineral separation production daily report or mineral separation experimental data, establishing a model for solving mineral separation recovery rate by taking the input grade of each component as an independent variable, wherein the model is shown in a formula 3:
(equation 3)
Wherein:beneficiation recovery of component i;
-a model for mineral separation recovery according to the grade of the beneficiation;
-grade of choice of component i;
according to the geological grades of the components obtained by calculation under each group of boundary grades and industrial grades, calculating the extracted grades by using a formula 4, and assuming that the ore of the concentrating mill is only sourced from the mine with the boundary grades and industrial grades to be determined, considering the amount of the ore to be selected as the extracted ore amount and the amount of the ore to be selected as the extracted grade;
(equation 4)
Wherein:-the produced grade of component i;
-geological grade of component i, obtained from equation 2;
h-rock mixing rate;
w is the grade of surrounding rock;
will beThe visual identity is +.>Substituting formula 3 to obtain->Then, the comprehensive grade is obtained by using the formula 5 and the formula 6>；
(equation 5)
(equation 6)
Wherein:-comprehensive grade;
-grade of the main component;
n-the amount of co-associated components available;
-the grade conversion factor of the co-associated component i;
-grade of co-associated component i;
mineral recovery of co-associated component i;
-the ton concentrate price of the co-associated component i;
-concentrate grade of the main component;
mineral recovery of the main component;
-the ton concentrate price of the main component;
the concentrate grade of component i is co-associated.
5. An electronic device comprising a processor and a memory having instructions stored therein, wherein the instructions are loaded and executed by the processor to implement the method of determining a low grade multi-metal ore grade index of any one of claims 1-3.
6. A computer readable storage medium having instructions stored therein, wherein the instructions are loaded and executed by a processor to implement the method of determining a low grade multi-metal ore grade index of any one of claims 1-3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311324867.8A CN117057633B (en) | 2023-10-13 | 2023-10-13 | Method and system for determining grade index of low-grade multi-metal ore |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311324867.8A CN117057633B (en) | 2023-10-13 | 2023-10-13 | Method and system for determining grade index of low-grade multi-metal ore |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117057633A CN117057633A (en) | 2023-11-14 |
| CN117057633B true CN117057633B (en) | 2024-01-02 |
Family
ID=88666754
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311324867.8A Active CN117057633B (en) | 2023-10-13 | 2023-10-13 | Method and system for determining grade index of low-grade multi-metal ore |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117057633B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632794A (en) * | 2018-12-18 | 2019-04-16 | 核工业北京地质研究院 | A kind of comprehensive identification method suitable for gold occurrence state in low-grade gold mineralising area |
| CN113338935A (en) * | 2021-04-21 | 2021-09-03 | 铜陵有色金属集团股份有限公司 | Method for delineating ore body |
| CN115660883A (en) * | 2022-07-20 | 2023-01-31 | 中国矿业大学 | Scientific mining planning method for coal and co-associated mineral products in whole life cycle |
| CN116342159A (en) * | 2023-02-24 | 2023-06-27 | 昆明理工大学 | Construction method of multi-metal ore recoverable economic grade model |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114266179B (en) * | 2021-12-23 | 2023-08-04 | 成都理工大学 | Finite element-based deposit drilling information processing and analyzing method and finite element-based deposit drilling information processing and analyzing device |
-
2023
- 2023-10-13 CN CN202311324867.8A patent/CN117057633B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632794A (en) * | 2018-12-18 | 2019-04-16 | 核工业北京地质研究院 | A kind of comprehensive identification method suitable for gold occurrence state in low-grade gold mineralising area |
| CN113338935A (en) * | 2021-04-21 | 2021-09-03 | 铜陵有色金属集团股份有限公司 | Method for delineating ore body |
| CN115660883A (en) * | 2022-07-20 | 2023-01-31 | 中国矿业大学 | Scientific mining planning method for coal and co-associated mineral products in whole life cycle |
| CN116342159A (en) * | 2023-02-24 | 2023-06-27 | 昆明理工大学 | Construction method of multi-metal ore recoverable economic grade model |
Non-Patent Citations (2)
| Title |
|---|
| 提高矿产资源利用率的探讨;朱天云;;矿产保护与利用(02);10-11 * |
| 矿山地球化学勘探与矿体主共伴生组分综合评价;方维萱;;金属矿山(08);110-116 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117057633A (en) | 2023-11-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Knox-Robinson | Vectorial fuzzy logic: a novel technique for enhanced mineral prospectivity mapping, with reference to the orogenic gold mineralisation potential of the Kalgoorlie Terrane, Western Australia | |
| Wu et al. | A data mining approach for spatial modeling in small area load forecast | |
| CN103049733B (en) | Method for detecting human face and human-face detection equipment | |
| CN107194803A (en) | A kind of P2P nets borrow the device of borrower's assessing credit risks | |
| CN112541532B (en) | Target detection method based on dense connection structure | |
| CN107122860B (en) | Rock burst danger level prediction method based on grid search and extreme learning machine | |
| Ghezelbash et al. | An improved data-driven multiple criteria decision-making procedure for spatial modeling of mineral prospectivity: adaption of prediction–area plot and logistic functions | |
| Bashari et al. | Estimation of deformation modulus of rock masses by using fuzzy clustering-based modeling | |
| CN116665067B (en) | Ore finding target area optimization system and method based on graph neural network | |
| CN110634060A (en) | User credit risk assessment method, system, device and storage medium | |
| CN108334988A (en) | A kind of short-term Load Forecasting based on SVM | |
| CN116522250A (en) | Mineral deposit type prediction method and system based on various algorithms | |
| CN117057633B (en) | Method and system for determining grade index of low-grade multi-metal ore | |
| CN110196797B (en) | Automatic optimization method and system suitable for credit scoring card system | |
| CN108876169A (en) | Uranium deposit Technological Economy dynamic evaluation method | |
| CN107908834B (en) | Three-dimensional positioning mineralization prediction method and system for blind ore body | |
| Mikaeil et al. | Ranking and evaluating the coal seam mechanization based on geological conditions | |
| CN112288571B (en) | Personal credit risk assessment method based on rapid construction of neighborhood coverage | |
| Samanta et al. | Sparse data division using data segmentation and Kohonen network for neural network and geostatistical ore grade modeling in Nome offshore placer deposit | |
| Esmaeiloghli et al. | Optimizing the grade classification model of mineralized zones using a learning method based on harmony search algorithm | |
| CN103886007A (en) | Mutual constraint based fuzzy data classification method | |
| CN116244426A (en) | Geographic function area identification method, device, equipment and storage medium | |
| Afradi et al. | Selecting the best mining method using analytical and numerical methods | |
| CN110929365A (en) | Oil sand resource amount calculation system | |
| CN108961056A (en) | Stock market's investor sentiment based on entropy dimension and maximum entropy ratio acts on behalf of index selection method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |