Method and device for predicting mineral reserves

Abstract

The application provides a mineral reserves prediction method and device, which can be used in the field of minerals. The mineral reserves prediction method and device provided by the embodiment of the application can perform reserves calculation according to a pre-established three-dimensional reserves model, obtain initial mineral reserves of each grade interval in a pre-divided elevation interval of a mining area, determine a boundary grade according to a current mineral market price and a pre-established price-boundary grade calculation model, remove a part of the initial mineral reserves that does not reach the boundary grade, and obtain predicted mineral reserves, which can quickly respond to market price fluctuations and make accurate economic evaluation and decision.

Description

Technical Field

[0001] This application relates to the field of mineral resources, and more particularly to a method and apparatus for predicting mineral reserves. Background Technology

[0002] Traditional methods for calculating reserves often treat the cutoff grade as a relatively fixed value, neglecting the impact of constantly changing market prices on it. In mining operations, the cutoff grade is one of the most fundamental and crucial parameters. Its selection directly affects the morphology of the delineated ore body and the ore reserves, thus influencing mine construction, production scale, profitability, equipment selection, and operational lifespan, while also impacting the utilization rate of non-renewable resources. The cutoff grade is a techno-economic parameter with a significant impact on the economic and social benefits of a mine. Therefore, traditional reserve estimation methods are often relatively lagging, unable to respond quickly to market price fluctuations, leading to difficulties for mining companies in making accurate economic assessments and decisions. Summary of the Invention

[0003] In view of the problems in the prior art, this application provides a mineral reserve prediction method and apparatus, which can at least partially solve the problems existing in the prior art.

[0004] To solve the above-mentioned technical problems, this application provides the following technical solution:

[0005] Firstly, this application provides a method for predicting mineral reserves, including:

[0006] Based on the pre-established three-dimensional reserve model, the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area are calculated.

[0007] The boundary grade is determined based on the current mineral market price and a pre-established price-boundary grade calculation model.

[0008] The predicted mineral reserves are obtained by removing the portion of the initial mineral reserves that does not reach the specified grade.

[0009] Secondly, this application provides a mineral reserves prediction device, comprising:

[0010] Thirdly, this application provides an electronic device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the mineral reserve prediction method described in any of the above embodiments.

[0011] Fourthly, this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the mineral reserve prediction method described in any of the above embodiments.

[0012] Fifthly, this application provides a computer program product, which includes a computer program that, when executed by a processor, implements the mineral reserve prediction method described in any of the above embodiments.

[0013] The mineral reserve prediction method and apparatus provided in this application can calculate reserves based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within a pre-divided elevation interval of the mining area; determine the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model; remove the portion of the initial mineral reserves that does not reach the boundary grade to obtain the predicted mineral reserves, and quickly respond to market price fluctuations to make accurate economic evaluations and decisions. Attached Figure Description

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

[0015] Figure 1 A flowchart of the mineral reserves prediction method provided in the embodiments of this application;

[0016] Figure 2 A flowchart for establishing a three-dimensional reserve model provided in this application embodiment;

[0017] Figure 3 A flowchart for establishing a digital terrain model is provided for embodiments of this application;

[0018] Figure 4 A flowchart for processing abnormal points or abnormal lines provided in the embodiments of this application;

[0019] Figure 5 A flowchart for drawing cross-sectional lines provided in this application embodiment;

[0020] Figure 6 A flowchart for establishing a three-dimensional ore body distribution model is provided for embodiments of this application;

[0021] Figure 7 A flowchart for processing openings or intersecting triangles provided in an embodiment of this application;

[0022] Figure 8 A flowchart for establishing a three-dimensional reserve model provided in this application embodiment;

[0023] Figure 9 A flowchart for obtaining initial mineral reserves provided in this application embodiment;

[0024] Figure 10 A flowchart for determining the boundary grade is provided for embodiments of this application;

[0025] Figure 11 A flowchart of the mineral reserves prediction method provided in the embodiments of this application;

[0026] Figure 12 A schematic diagram of the device structure of the mineral reserves prediction device provided in the embodiments of this application;

[0027] Figure 13 A schematic diagram of the device structure of the mineral reserves prediction device provided in the embodiments of this application;

[0028] Figure 14 A schematic diagram of the device structure of the mineral reserves prediction device provided in the embodiments of this application;

[0029] Figure 15 A schematic diagram of the device structure of the mineral reserves prediction device provided in the embodiments of this application;

[0030] Figure 16 A schematic diagram of the device structure of the mineral reserves prediction device provided in the embodiments of this application; Detailed Implementation

[0031] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the embodiments of this application will be further described in detail below with reference to the accompanying drawings. Here, the illustrative embodiments and their descriptions are used to explain this application, but are not intended to limit this application. It should be noted that, unless otherwise specified, the embodiments and features in the embodiments of this application can be arbitrarily combined with each other.

[0032] Figure 1 This is a flowchart of the mineral reserve prediction method provided in the embodiments of this application. Figure 1 As shown in the embodiments of this application, the mineral reserve prediction method includes:

[0033] S101: Calculate reserves based on the pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area.

[0034] S102: Determine the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model;

[0035] S103: Remove the portion of the initial mineral reserves that does not reach the cutoff grade to obtain the predicted mineral reserves.

[0036] Specifically, the grade range has an upper limit grade and a lower limit grade. For grade ranges where the upper limit grade does not reach the boundary grade, the reserves within that grade range are directly removed. For grade ranges where the lower limit grade does not reach the boundary grade but the upper limit grade does reach the boundary grade, the reserves within that grade range are removed, and the boundary grade is used as the new lower limit grade, while the upper limit grade remains unchanged, resulting in a new grade range. The reserves are then recalculated based on the three-dimensional reserve model to obtain the mineral reserves within the pre-divided elevation range of the mining area for that grade range.

[0037] The mineral reserve prediction method and apparatus provided in this application can calculate reserves based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within a pre-divided elevation interval of the mining area; determine the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model; remove the portion of the initial mineral reserves that does not reach the boundary grade to obtain the predicted mineral reserves, and quickly respond to market price fluctuations to make accurate economic evaluations and decisions.

[0038] like Figure 2 As shown, based on the above embodiments, the further step of establishing the three-dimensional reserve model includes:

[0039] S201: Establish a three-dimensional ore body distribution model based on the pre-established digital ground model and the pre-drawn ore body profile lines;

[0040] S202: Estimate the grade distribution within the deposit in the three-dimensional orebody distribution model and calculate the reserves based on the grade values ​​of the mineral assemblage sample and some basic blocks of the three-dimensional orebody distribution model, thus obtaining the three-dimensional reserve model.

[0041] like Figure 3 As shown, based on the above embodiments, the steps for establishing the digital ground model further include:

[0042] S301: Obtain the mine area plan and convert the curve segments in the mine area plan into polylines;

[0043] Specifically, the mining area plan is obtained, and the curve segments in the mining area plan are converted into multi-line segments. Since many mining-specific software do not provide support for curve segments, it is necessary to convert the curve segments into multi-line segments. The mining area plan is drawn based on historical exploration data. The mining area plan has contour lines and measuring points, and the contour lines and measuring points have corresponding heights, which can be used for subsequent operations.

[0044] S302: Assign elevation and color configuration to multiple line segments to generate a digital ground model.

[0045] Specifically, for multiple line segments corresponding to contour lines, elevations are assigned based on the marked heights. For multiple line segments corresponding to line segments without marked heights on the mining area plan, the heights can be calculated proportionally based on the measuring points and contour lines within a certain range around the unmarked line segments, and the calculated elevations are then assigned. For multiple line segments with assigned elevations, the heights can be manually edited. After assigning elevations, color ranges can be configured so that the generated digital ground model displays different heights through gradient colors.

[0046] like Figure 4 As shown, based on the above embodiments, before performing elevation and color configuration operations on the multi-line segments to generate the digital ground model, the following steps are included:

[0047] S401: Detects whether there are abnormal points or abnormal lines in a multi-line segment;

[0048] Specifically, it detects whether there are abnormal points or abnormal lines in the multi-line segment. Abnormal points are points that will cause large protrusions in the generated digital ground model, and abnormal lines are broken lines with excessively small included angles.

[0049] S402: Correct abnormal points or lines;

[0050] Specifically, abnormal points or lines can be corrected by deleting or trimming them.

[0051] The step of assigning elevation and color configuration to the multi-line segments to generate the digital ground model includes:

[0052] S403: Assign elevation and color configuration to the corrected multi-line segments to generate a digital ground model.

[0053] like Figure 5 As shown, based on the above embodiments, the further step of drawing the ore body profile line includes:

[0054] S501: Obtain the ore body profile and rotate and translate it so that each point in the ore body profile coincides with the corresponding point in the digital ground model;

[0055] Specifically, a ore body profile is obtained and rotated and translated so that each point in the ore body profile coincides with the corresponding point in the digital ground model. The ore body profile is drawn based on historical exploration data. The rotation and translation method can be selected as needed, and this application does not impose any restrictions.

[0056] For example, the angle of deviation between the positive direction of the ore body profile and the corresponding positive direction of the digital ground model can be measured. The ore body profile can be rotated by setting the rotation base point coordinates and rotation azimuth angle, and the translation base point coordinates and target point coordinates can be selected. By moving the translation base point to the position of the target point, the ore body profile can be translated so that each point in the ore body profile coincides with the corresponding point in the digital ground model.

[0057] S502: Convert the curve segments in the ore body profile into multiple line segments, which are the ore body profile lines.

[0058] Based on the above embodiments, the mineral reserve prediction method provided in this application further includes:

[0059] Adjust the range, distribution density, and density of points on the ore body profile lines.

[0060] Specifically, the adjustments include at least one or a combination of multi-segment thinning, multi-segment densification, multi-segment extension, and adding a spatial fitting line. A multi-segment is composed of lines connecting points, and the ore body profile line is a multi-segment. Multi-segment thinning involves deleting some points on the selected multi-segment without changing its basic shape, thus reducing the point density. Multi-segment densification involves inserting a certain number of points into the selected multi-segment to meet operational requirements. Multi-segment extension extends one segment or increases the number of segments when the selected multi-segment is insufficient. Adding a spatial fitting line adds a new multi-segment between two existing multi-segments according to the spacing ratio.

[0061] like Figure 6 As shown, based on the above embodiments, the step of establishing a three-dimensional ore body distribution model according to a pre-established digital ground model and pre-drawn ore body profile lines further includes:

[0062] S601: Couple the digital ground model with the ore body profile line;

[0063] Specifically, the digital ground model is coupled with the ore body profile line. In the above operation, the points in the ore body profile map have been made to coincide with the corresponding points in the digital ground model through rotation and translation. Therefore, the ore body profile line obtained from the ore body profile map is also located at the corresponding position in the digital ground model and can be directly coupled with the ore body profile model.

[0064] S602: Connect adjacent ore body profiles in the coupling results to form an ore body, thus obtaining a three-dimensional ore body distribution model.

[0065] Specifically, adjacent orebody profile lines in the coupling results are connected to form an orebody, resulting in a three-dimensional orebody distribution model. Since the orebody profile lines obtained from each orebody profile are located at their corresponding positions in the digital terrain model, the profile lines from different profiles of the same orebody are adjacent in the digital terrain model. By connecting adjacent profile lines, a three-dimensional model of the orebody can be obtained. Because the orebody profile lines converted from the orebody profiles of multiple orebody profiles are all coupled with the digital terrain model, there is more than one three-dimensional model of the orebody formed by connecting them, thus obtaining a three-dimensional orebody distribution model. The three-dimensional model of the orebody is hereinafter referred to as the orebody.

[0066] like Figure 7 As shown, based on the above embodiments, after connecting adjacent ore body profiles in the coupling result to form an ore body and obtaining the three-dimensional ore body distribution model, the method further includes:

[0067] S701: Detect whether there are openings or intersecting triangles in the ore body distribution model in the three-dimensional ore body distribution model;

[0068] Specifically, it detects whether there are openings or intersecting triangles in the ore body distribution model. If there are openings or intersecting triangles, it will lead to the inability to determine the minimum outer perimeter of the ore body or the repeated calculation of some ore bodies.

[0069] S702: Adjust the local triangular mesh of the ore body.

[0070] Specifically, if open or intersecting triangles exist, adjust the local triangulation of the ore body to eliminate them. When connecting adjacent profile lines, some profile lines may not be connected correctly, resulting in open areas. This can be addressed by merging entities to adjust the triangulation at the open areas. Merging two or more entities belonging to the same ore body that are not correctly connected due to partial profile line misconnection will yield the correct ore body. When intersecting triangles appear, delete the triangulation at the intersection and stitch in a new triangulation to eliminate them.

[0071] like Figure 8 As shown, based on the above embodiments, further, the step of estimating the grade distribution within the deposit in the three-dimensional orebody model and calculating reserves based on the grade values ​​of the mineral assemblage sample and the grade values ​​of some basic blocks of the three-dimensional orebody distribution model to obtain a three-dimensional reserve model includes:

[0072] S801: Determine the sample length of mineral assemblage samples by statistically analyzing the sample length distribution based on the sample file data in the pre-established geological engineering database.

[0073] Specifically, the sample length distribution is statistically analyzed based on the sample file data in the pre-established geological engineering database to obtain the length of the concentrated distribution of samples. This length is then used as the sample length of the mineral assemblage sample. The method of statistically analyzing the sample length distribution can be selected as needed, and this application does not impose any restrictions.

[0074] For example, a bar chart can be drawn to visually show the length of the concentrated distribution of samples.

[0075] S802: Combine grades based on sample length and calculate the grade value of the mineral combination sample;

[0076] Specifically, since the grade data of the samples obtained from the borehole are all continuous sample grade data, the samples can be combined in terms of grade. The borehole depth is divided according to the sample length, and the weighted average of the grade of each sample within the sample length is calculated to obtain the grade value of the mineral combination sample, so as to reduce the impact of extreme data on the estimated grade value.

[0077] S803: Establish the minimum envelope of ore bodies in the three-dimensional ore body distribution model;

[0078] S804: Add constraints to the minimum outer package to obtain the block model;

[0079] Specifically, by adding constraints, the minimum outer surface is graded to make it fit the ore body better, resulting in a minimum outer surface with added constraints, i.e., a block model.

[0080] S805: Divide the block model into multiple basic blocks. The length of the basic block is equal to the sample length of the mineral assemblage sample. The basic block includes internal blocks and boundary blocks. The internal blocks have the same size, the boundary blocks have the same size, and the size of the boundary blocks is smaller than that of the internal blocks.

[0081] Specifically, the block model is divided into multiple basic blocks. The length of a basic block is equal to the sample length of a mineral assemblage sample. The basic blocks include internal blocks and boundary blocks. The internal blocks have the same size, and the boundary blocks have the same size. The size of the boundary blocks is smaller than that of the internal blocks, making the boundaries more precise.

[0082] S806: Determine whether the grade and reserve values ​​of all the aforementioned basic blocks have been obtained;

[0083] Specifically, since there may be no valued base blocks and / or mineral combination samples within the search ellipsoid centered on the current base block, not every base block can be valued after one round of valuation, requiring multiple rounds of iterative valuation. If the grade and reserve values ​​of all the base blocks are not obtained, proceed to S807; if the grade and reserve values ​​of all the base blocks have been obtained, proceed to S816.

[0084] S807: Creates a search ellipsoid based on a preset search radius;

[0085] Specifically, a search ellipsoid is created based on a preset search radius, which includes the two equatorial radii and one polar radius of the ellipsoid.

[0086] S808: Move the search ellipsoid so that the current base block to be valued is located at the center of the search ellipsoid;

[0087] S809: Search for whether there are pre-valued basic blocks and / or the mineral combination samples within the space defined by the search ellipsoid;

[0088] Specifically, if the space defined by the search ellipsoid contains both valued basic blocks and mineral combination samples, proceed to S810; if the space defined by the search ellipsoid contains only valued basic blocks, proceed to S811; if the space defined by the search ellipsoid contains only the mineral combination samples, proceed to S812.

[0089] S810: Estimate the grade value of the current base block to be valued based on the grade value of the mineral assemblage sample and the grade value of the already valued base block;

[0090] S811: Estimate the grade of the current base block to be valued based on the grade of the base block that has already been valued;

[0091] S812: Estimate the grade of the current base block to be valued based on the grade values ​​of the mineral assemblage sample;

[0092] S813: Calculate the reserve value of the base block based on the volume and grade value of the base block to be valued;

[0093] Specifically, the reserves of a base block can be calculated using its volume, density, and grade, where density is a constant.

[0094] S814: Determine if there are any remaining base blocks that have not been searched as the center of the search ellipsoid in this round and are yet to be valued;

[0095] Specifically, if there are still basic blocks to be valued that have not been searched as the center of the search ellipsoid, the process will cycle through to S808; if there are no basic blocks to be valued that have not been searched as the center of the search ellipsoid, the process will cycle through to S806.

[0096] S815: Save the base block grade and reserve values ​​obtained from this estimation, and expand the search radius of the search ellipsoid.

[0097] Specifically, the base block grade and reserve values ​​obtained from this estimation are saved, the search radius of the search ellipsoid is expanded, and the process is repeated in S806.

[0098] S816: Obtain a three-dimensional reserve model.

[0099] like Figure 9 As shown, based on the above embodiments, the step of calculating reserves according to a pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area includes:

[0100] S901: Based on the mapping relationship between the spatial location and elevation interval of all basic blocks, and the mapping relationship between the grade value and grade interval of all basic blocks, the sample blocks contained in each grade interval within the elevation interval are obtained.

[0101] S902: Sum the reserve values ​​of the sample blocks contained in each grade interval within the elevation interval to obtain the mineral reserves of each grade interval within the elevation interval.

[0102] like Figure 10 As shown, based on the above embodiments, the step of determining the boundary grade according to the current mineral market price and the pre-established price-boundary grade calculation model further includes:

[0103] S1001: Based on historical data of the mining area, obtain the mineral processing cost, mining cost, mineral processing recovery rate, and mining recovery rate of the mining area;

[0104] Specifically, based on historical data of the mining area, the beneficiation cost, mining cost, beneficiation recovery rate, and mining recovery rate of the mining area are obtained. The beneficiation cost, mining cost, beneficiation recovery rate, and mining recovery rate of the mining area are all weighted averages of historical data.

[0105] S1002: When the profit break-even point is reached, establish a balance between cost and profit, and further obtain the balance between the beneficiation cost of the mining area, the mining cost of the mining area, the beneficiation recovery rate of the mining area, the mining recovery rate of the mining area, the boundary grade and the market price.

[0106] S1003: Based on the balance between the mineral processing cost, mining cost, mineral processing recovery rate, mining recovery rate, cutoff grade, and market price of the mining area, a price-cutoff grade calculation model is obtained.

[0107] Specifically, by organizing the balance between the beneficiation cost, mining cost, beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area, a price-cutoff grade calculation model is obtained.

[0108] S1004: Input the current mineral market price into the price-boundary grade calculation model to obtain the boundary grade.

[0109] Based on the above embodiments, the step of establishing a three-dimensional ore body distribution model based on the digital ground model and profile lines further includes:

[0110] Drilling data is obtained through exploration drilling to establish a geological engineering database. The drilling data includes borehole file data and sample file data.

[0111] Specifically, the orifice file data includes the orifice number, orifice coordinates, and maximum orifice depth, while the sample file data includes the orifice number from which the sample was obtained, the sample start depth, the sample end depth, the sample number, the sample length, and the sample grade.

[0112] The following uses the prediction of graphite mineral reserves in Yunshan as an example to illustrate the specific implementation process of the mineral reserve prediction method provided in this application embodiment. The specific implementation process is as follows: Figure 11 As shown.

[0113] S1101: Obtain the mine area plan and convert the curve segments in the mine area plan into polylines;

[0114] S1102: Assign elevation and color configuration to multiple line segments to generate a digital ground model;

[0115] S1103: Obtain the ore body profile and rotate and translate it so that each point in the ore body profile coincides with the corresponding point in the digital ground model;

[0116] S1104: Convert the curve segments in the ore body profile diagram into multi-line segments, which are the ore body profile lines;

[0117] S1105: Couple the digital ground model with the ore body profile line;

[0118] S1106: Connect adjacent ore body profiles in the coupling results to form an ore body, and obtain a three-dimensional ore body distribution model;

[0119] S1107: Based on the statistical distribution of sample lengths in the sample files of the pre-established geological engineering database, the sample length of the mineral assemblage sample is determined to be 2.0m;

[0120] Specifically, by drawing a statistical chart, it was found that the sample lengths were concentrated at 2.0m, so the sample length of the mineral assemblage sample was determined to be 2.0m.

[0121] S1108: Combine grades based on sample length and calculate the grade value of the mineral combination sample;

[0122] S1109: Establish the minimum bounding box of ore bodies in the three-dimensional ore body distribution model;

[0123] S1110: Add constraints to the minimum outer package to obtain the block model;

[0124] S1111: Divide the block model into multiple basic blocks. The length of the basic block is equal to the sample length of the mineral assemblage sample. The basic block includes internal blocks and boundary blocks. The internal blocks have the same size, the boundary blocks have the same size, and the size of the boundary blocks is smaller than that of the internal blocks.

[0125] Specifically, the dimensions of the internal blocks are 2.0m × 4.0m × 8.0m, and the dimensions of the boundary blocks are 2.0m × 2.0m × 2.0m.

[0126] S1112: Determine whether the grade and reserve values ​​of all the aforementioned basic blocks have been obtained;

[0127] Specifically, if the grade and reserve values ​​of all the said basic blocks are not obtained, proceed to S1113; if the grade and reserve values ​​of all the said basic blocks have been obtained, proceed to S1122.

[0128] S1113: Create a search ellipsoid based on the preset search radius;

[0129] Specifically, the initial search radius is set to 16m, and a new search radius is set for each subsequent search round.

[0130] S1114: Move the search ellipsoid so that the current base block to be valued is located at the center of the search ellipsoid;

[0131] S1115: Search for whether there are pre-valued basic blocks and / or the mineral combination samples within the space defined by the search ellipsoid;

[0132] Specifically, if the space defined by the search ellipsoid contains both the valued basic block and the mineral combination sample, proceed to S1116; if the space defined by the search ellipsoid contains only the valued basic block, proceed to S1117; if the space defined by the search ellipsoid contains only the mineral combination sample, proceed to S1118.

[0133] S1116: Estimate the grade value of the current base block to be valued based on the grade value of the mineral assemblage sample and the grade value of the already valued base block;

[0134] Specifically, assuming a mineral composite sample with a grade of 0.5 is found and its distance from the current base block to be valued is 20, and two base blocks have already been valued, one with a grade of 0.46 and a distance of 4 from the current base block to be valued, and the other with a grade of 0.48 and a distance of 5 from the current base block to be valued, the grade of the current base block to be valued can be estimated as 0.47 by using the inverse power of distance method to estimate the value of the current base block to be valued.

[0135] S1117: Estimate the grade of the current base block to be valued based on the grade of the base block that has already been valued;

[0136] Specifically, assuming two valued basic blocks are found, one with a quality value of 0.24 and a distance of 4 from the current basic block to be valued, and the other with a quality value of 0.25 and a distance of 5 from the current basic block to be valued, the quality value of the current basic block to be valued can be estimated by using the inverse power of distance method.

[0137] S1118: Estimate the grade of the current base block to be valued based on the grade values ​​of the mineral assemblage sample;

[0138] Specifically, assuming a mineral combination sample is found with a grade of 0.36 and a distance of 4 from the current base block to be valued, the grade of the current base block to be valued can be estimated by using the inverse power of distance method.

[0139] S1119: Calculate the reserve value of the base block based on the volume and grade value of the base block to be valued;

[0140] Specifically, assuming the current base block is an internal block, the estimated grade is 0.24, and the density of the graphite ore is 0.22 g / cm³. 3 Therefore, the current reserve value of the base block to be valued is 3.38 tons.

[0141] S1120: Determine if there are any remaining base blocks that have not been searched as the center of the search ellipsoid in this round and are yet to be valued;

[0142] Specifically, if there are still basic blocks to be valued that have not been searched as the center of the search ellipsoid, the process loops to S1114; if there are no basic blocks to be valued that have not been searched as the center of the search ellipsoid, the process loops to S1121.

[0143] S1121: Save the base block grade and reserve values ​​obtained in this estimation, and expand the search radius of the search ellipsoid.

[0144] Specifically, the base block grade and reserve values ​​obtained from this estimation are saved, the search radius of the search ellipsoid is expanded, and the process is repeated in S1112.

[0145] S1122: Obtain the three-dimensional reserve model;

[0146] S1123: Based on the mapping relationship between the spatial location and elevation interval of all basic blocks, and the mapping relationship between the grade value and grade interval of all basic blocks, the sample blocks contained in each grade interval within the elevation interval are obtained.

[0147] S1124: Sum the reserve values ​​of the sample blocks contained in each grade interval within the elevation interval to obtain the mineral reserves of each grade interval within the elevation interval.

[0148] S1125: Based on historical data of the mining area, obtain the mineral processing cost, mining cost, mineral processing recovery rate, and mining recovery rate of the mining area;

[0149] Specifically, the beneficiation cost is 2,372 yuan / ton concentrate, the mining cost is 32.06 yuan / ton, the beneficiation recovery rate is 0.85, and the mining recovery rate is 0.91.

[0150] S1126: When the profit break-even point is reached, establish a balance between cost and profit. Based on the relationship between cost and the beneficiation cost, mining cost, beneficiation recovery rate, mining recovery rate, and cutoff grade of the mining area, and the relationship between profit and the beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area, further obtain the balance between the beneficiation cost, mining cost, beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area.

[0151] S1127: Based on the balance between the mineral processing cost, mining cost, mineral processing recovery rate, mining recovery rate, cutoff grade, and market price of the mining area, a price-cutoff grade calculation model is obtained.

[0152] S1128: Input the current mineral market price into the price-cutoff grade calculation model to obtain the cutoff grade;

[0153] S1129: Remove the portion of the initial mineral reserves that does not reach the cutoff grade to obtain the predicted mineral reserves.

[0154] The mineral reserve prediction method provided in this application can calculate reserves based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves for each grade interval within a pre-divided elevation interval of the mining area; based on the current mineral market price of 3,500 yuan / ton and the pre-established price-boundary grade calculation model, the boundary grade is determined to be 0.04; after removing the portion of the initial mineral reserves of 227 million tons that does not reach the boundary grade, the predicted mineral reserves are obtained as 216 million tons, which can quickly react to market price fluctuations and make accurate economic evaluations and decisions.

[0155] Figure 12 This is a schematic diagram of the mineral reserve prediction device provided in the embodiments of this application, as shown below. Figure 12 As shown, the mineral reserve prediction device provided in this application embodiment includes a boundary grade calculation unit 1210 and a reserve calculation unit 1220, wherein:

[0156] The boundary grade calculation unit 1210 is used to determine the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model.

[0157] The reserve calculation unit 1220 is used to perform reserve calculation based on a pre-established three-dimensional reserve model, to obtain the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area, and to remove the portion of the initial mineral reserves that does not reach the boundary grade, so as to obtain the predicted mineral reserves.

[0158] The mineral reserve prediction device provided in this application embodiment can calculate reserves based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area; determine the boundary grade of 0.04 based on the current mineral market price of 3,500 yuan and the pre-established price-boundary grade calculation model; remove the portion of the initial mineral reserves of 227 million tons that does not reach the boundary grade to obtain the predicted mineral reserves of 216 million tons, which can quickly react to market price fluctuations and make accurate economic evaluations and decisions.

[0159] Figure 13 This is a schematic diagram of the mineral reserve prediction device provided in the embodiments of this application, as shown below. Figure 13 As shown, based on the above embodiments, the mineral reserve prediction device provided in this embodiment further includes a model building unit 1230, which includes a three-dimensional modeling module 1231 for:

[0160] A three-dimensional orebody distribution model is established based on a pre-established digital ground model and pre-drawn orebody profiles. The grade distribution within the ore deposit in the three-dimensional orebody distribution model is estimated and the reserves are calculated based on the grade values ​​of mineral assemblage samples and the grade values ​​of some basic blocks of the three-dimensional orebody distribution model, thus obtaining a three-dimensional reserve model.

[0161] Based on the above embodiments, the three-dimensional modeling module 1231 in the model building unit 1230 is further used for:

[0162] Obtain the mining area plan, convert the curve segments in the mining area plan into multi-line segments, assign elevation and color configuration operations to the multi-line segments, and generate a digital ground model.

[0163] Figure 14 This is a schematic diagram of the mineral reserve prediction device provided in the embodiments of this application, as shown below. Figure 14 As shown, based on the above embodiments, the mineral reserve prediction device provided in this embodiment of the invention further includes a model correction unit 1240, used for:

[0164] Detect whether there are abnormal points or lines in the multi-line segment; if so, correct the abnormal points or lines.

[0165] Based on the above embodiments, the three-dimensional modeling module 1231 in the model building unit 1230 is further used for:

[0166] Obtain the ore body profile and rotate and translate it so that each point in the ore body profile coincides with the corresponding point in the digital ground model; convert the curve segments in the ore body profile into multi-line segments, which are the ore body profile lines.

[0167] Based on the above embodiments, the model correction unit 1240 is further configured to:

[0168] Adjustments are made to the extent, distribution density, and density of points on the ore body profile lines; the adjustments include at least one or a combination of multiple line segment thinning, multiple line segment densification, multiple line segment extension, and adding spatial fitting lines.

[0169] Based on the above embodiments, the three-dimensional modeling module 1231 in the module establishment unit 1230 is further used for:

[0170] By coupling the digital ground model with the ore body profile lines, and connecting adjacent ore body profile lines in the coupling result to form an ore body, a three-dimensional ore body distribution model is obtained.

[0171] Based on the above embodiments, the model correction unit 1240 is further configured to:

[0172] The system detects whether there are open or intersecting triangles in the ore body distribution model. If open or intersecting triangles exist, the local triangulation of the ore body is adjusted to eliminate them.

[0173] Based on the above embodiments, further, the step of estimating the grade distribution within the deposit in the three-dimensional orebody model and calculating reserves based on the grade values ​​of the mineral assemblage sample and the grade values ​​of some basic blocks of the three-dimensional orebody distribution model to obtain the three-dimensional reserve model includes:

[0174] Based on the statistical distribution of sample lengths in the sample files of the pre-established geological engineering database, the sample length of the mineral assemblage sample is determined.

[0175] Grade combinations are performed based on the sample length, and the grade value of the mineral combination sample is calculated.

[0176] Establish the minimum bounding box of the ore body in the three-dimensional ore body distribution model;

[0177] Add constraints to the minimum outer package to obtain the block model;

[0178] The block model is divided into multiple basic blocks. The length of each basic block is equal to the sample length of the mineral assemblage sample. Each basic block includes internal blocks and boundary blocks. The internal blocks have the same size, the boundary blocks have the same size, and the size of the boundary blocks is smaller than that of the internal blocks.

[0179] Perform the following iterative operations until the grade and reserve values ​​of all the aforementioned base blocks are obtained, thus obtaining a three-dimensional reserve model:

[0180] Create a search ellipsoid based on the preset search radius;

[0181] Iterate through all the underlying blocks to be valued and perform the following operations:

[0182] Move the search ellipsoid so that the current base block to be valued is located at the center of the search ellipsoid;

[0183] Search for whether there are valued basic blocks and / or mineral combination samples within the space defined by the search ellipsoid;

[0184] If there are already valued basic blocks and the mineral combination sample within the space defined by the search ellipsoid, estimate the grade value of the current basic block to be valued based on the grade value of the mineral combination sample and the grade value of the already valued basic blocks.

[0185] The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued.

[0186] If only pre-valued basic blocks exist within the space defined by the search ellipsoid, estimate the quality value of the current basic block to be valued based on the quality value of the pre-valued basic blocks.

[0187] The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued.

[0188] If only the mineral assemblage sample exists within the space defined by the search ellipsoid, the grade value of the current base block to be valued is estimated based on the grade value of the mineral assemblage sample.

[0189] The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued.

[0190] After the traversal is completed, the base block grade and reserve values ​​obtained in this estimation are saved, and the search radius of the search ellipsoid is expanded.

[0191] Based on the above embodiments, the reserve calculation unit 1220 is further specifically used for:

[0192] Based on the mapping relationship between the spatial location and elevation interval of all basic blocks, and the mapping relationship between the grade value and grade interval of all basic blocks, the sample blocks contained in each grade interval within the elevation interval are obtained; the reserve values ​​of the sample blocks contained in each grade interval within the elevation interval are summed to obtain the mineral reserves of each grade interval within the elevation interval.

[0193] Figure 15 This is a schematic diagram of the mineral reserve prediction device provided in the embodiments of this application, as shown below. Figure 15 As shown, based on the above embodiments, the model building unit 1230 of the mineral reserve prediction device provided in this embodiment further includes a calculation model modeling module 1232, used for:

[0194] Based on historical data from the mining area, the beneficiation cost, mining cost, beneficiation recovery rate, and mining recovery rate of the mining area are obtained. When the break-even point is reached, a balance relationship between cost and profit is established. Based on the relationships between cost and the beneficiation cost, mining cost, beneficiation recovery rate, mining recovery rate, and cutoff grade of the mining area, and the relationships between profit and the beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area, a balance relationship between the beneficiation cost, mining cost, beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area is further obtained. Based on the balance relationship between the beneficiation cost, mining cost, beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area, the price-cutoff grade calculation model is obtained.

[0195] Figure 16 This is a schematic diagram of the mineral reserve prediction device provided in the embodiments of this application, as shown below. Figure 16 As shown, based on the above embodiments, the mineral reserve prediction device provided in this embodiment of the invention further includes a data storage unit 1250, used for:

[0196] The system stores the calculated reserve data and the borehole data obtained through exploration drilling to establish a geological engineering database. The borehole data includes borehole file data and sample file data.

[0197] This embodiment provides a computer electronic device, which may include a processor, a communication interface, memory, and a communication bus. The processor, communication interface, and memory communicate with each other via the communication bus. The processor can invoke logical instructions from the memory to execute the following methods:

[0198] Based on the pre-established three-dimensional reserve model, the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area are calculated.

[0199] The boundary grade is determined based on the current mineral market price and a pre-established price-boundary grade calculation model.

[0200] The predicted mineral reserves are obtained by removing the portion of the initial mineral reserves that does not reach the specified grade.

[0201] The computer electronic device provided in this application embodiment can perform reserve calculations based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within a pre-divided mining area elevation interval; determine the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model; remove the portion of the initial mineral reserves that does not reach the boundary grade to obtain the predicted mineral reserves, and quickly react to market price fluctuations to make accurate economic evaluations and decisions.

[0202] Furthermore, the logical instructions in the aforementioned memory can be implemented as software functional units and sold or used as independent products, and can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0203] This embodiment discloses a computer program product, which includes a computer program stored on a computer-readable storage medium. The computer program includes program instructions, and when the program instructions are executed by a computer, the computer can perform the methods provided in the above-described method embodiments, such as:

[0204] Based on the pre-established three-dimensional reserve model, the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area are calculated.

[0205] The boundary grade is determined based on the current mineral market price and a pre-established price-boundary grade calculation model.

[0206] The predicted mineral reserves are obtained by removing the portion of the initial mineral reserves that does not reach the specified grade.

[0207] The computer program product provided in this application embodiment can perform reserve calculations based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within a pre-divided mining area elevation interval; determine the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model; remove the portion of the initial mineral reserves that does not reach the boundary grade to obtain the predicted mineral reserves, and quickly react to market price fluctuations to make accurate economic evaluations and decisions.

[0208] This embodiment provides a computer-readable storage medium storing a computer program that causes the computer to execute the methods provided in the above-described method embodiments, including, for example:

[0209] Based on the pre-established three-dimensional reserve model, the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area are calculated.

[0210] The boundary grade is determined based on the current mineral market price and a pre-established price-boundary grade calculation model.

[0211] The predicted mineral reserves are obtained by removing the portion of the initial mineral reserves that does not reach the specified grade.

[0212] The readable storage medium provided in this application embodiment can perform reserve calculations based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within a pre-divided mining area elevation interval; determine the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model; remove the portion of the initial mineral reserves that does not reach the boundary grade to obtain the predicted mineral reserves, and quickly respond to market price fluctuations to make accurate economic evaluations and decisions.

[0213] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0214] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0215] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0216] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0217] In the description of this specification, the references to terms such as "an embodiment," "a specific embodiment," "some embodiments," "for example," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples.

[0218] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of this application. It should be understood that the above descriptions are merely specific embodiments of this application and are not intended to limit the scope of protection of this application. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of protection of this application.

Claims

1. A method for predicting mineral reserves, characterized in that, include: Based on the pre-established three-dimensional reserve model, the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area are calculated. The following iterative operations are performed until the grade and reserve values ​​of all basic blocks are obtained, resulting in a three-dimensional reserve model: Create a search ellipsoid based on the preset search radius; Iterate through all the underlying blocks to be valued and perform the following operations: Move the search ellipsoid so that the current base block to be valued is located at the center of the search ellipsoid; Search for whether there are pre-valued basic blocks and / or mineral combination samples within the space defined by the search ellipsoid; If there are already valued basic blocks and the mineral combination sample within the space defined by the search ellipsoid, estimate the grade value of the current basic block to be valued based on the grade value of the mineral combination sample and the grade value of the already valued basic blocks. The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued. If only pre-valued basic blocks exist within the space defined by the search ellipsoid, estimate the quality value of the current basic block to be valued based on the quality value of the pre-valued basic blocks. The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued. If only the mineral assemblage sample exists within the space defined by the search ellipsoid, the grade value of the current base block to be valued is estimated based on the grade value of the mineral assemblage sample. The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued. After the traversal is completed, save the basic block grade value and reserve value obtained in this estimation, and expand the search radius of the search ellipsoid; The boundary grade is determined based on the current mineral market price and a pre-established price-boundary grade calculation model. The portion of the initial mineral reserves that does not reach the cutoff grade is removed to obtain the predicted mineral reserves. The step of calculating reserves based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves for each grade interval within the pre-divided elevation interval of the mining area includes: Based on the mapping relationship between the spatial location and elevation interval of all basic blocks, and the mapping relationship between the grade value and grade interval of all basic blocks, the sample blocks contained in each grade interval within the elevation interval are obtained. The mineral reserves of each grade interval within the elevation interval are obtained by summing the reserve values ​​of the sample blocks contained in each grade interval within the elevation interval. Specifically, the grade range has an upper limit grade and a lower limit grade. For grade ranges where the upper limit grade does not reach the boundary grade, the reserves within that grade range are directly removed. For grade ranges where the lower limit grade does not reach the boundary grade but the upper limit grade does reach the boundary grade, the reserves within that grade range are removed, and the boundary grade is used as the new lower limit grade, while the upper limit grade remains unchanged, resulting in a new grade range. The reserves are then recalculated based on the three-dimensional reserve model to obtain the mineral reserves within the pre-divided elevation range of the mining area for that grade range.

2. The mineral reserve prediction method according to claim 1, characterized in that, The steps for establishing the three-dimensional reserve model include: A three-dimensional ore body distribution model is established based on a pre-established digital ground model and pre-drawn ore body profile lines; Based on the grade values ​​of the mineral assemblage sample and the grade values ​​of the estimated base blocks of the three-dimensional orebody distribution model, the grade distribution within the ore deposit in the three-dimensional orebody distribution model is estimated and the reserves are calculated to obtain a three-dimensional reserve model.

3. The mineral reserve prediction method according to claim 2, characterized in that, The steps for establishing the digital terrain model include: Obtain a mining area plan and convert the curve segments in the mining area plan into polylines; The digital ground model is generated by assigning elevations and colors to the multi-line segments.

4. The mineral reserve prediction method according to claim 3, characterized in that, Before generating the digital terrain model by assigning elevation and color configurations to the multi-line segments, the process includes: Detect whether there are any abnormal points or abnormal lines in the multi-line segment; If they exist, correct the abnormal points or abnormal lines; The step of assigning elevation and color configuration to the multi-line segments to generate the digital ground model includes: The corrected multi-line segments are then assigned elevations and colors to generate the digital ground model.

5. The mineral reserve prediction method according to claim 3 or 4, characterized in that, The steps for drawing the profile lines of the ore body include: Obtain a cross-sectional view of the ore body and rotate and translate it so that each point in the cross-sectional view of the ore body coincides with the corresponding point in the digital ground model; The curve segments in the ore body profile are converted into multiple line segments, which are the ore body profile lines.

6. The mineral reserve prediction method according to claim 5, characterized in that, Also includes: The range, distribution density, and density of points on the ore body profile line are adjusted; the adjustment includes at least one or a combination of multiple line segment thinning, multiple line segment densification, multiple line segment extension, and adding spatial fitting lines.

7. The mineral reserve prediction method according to claim 2, characterized in that, The process of establishing a three-dimensional ore body distribution model based on a pre-established digital ground model and pre-drawn ore body profile lines includes: Couple the digital ground model with the ore body profile line; Connecting adjacent ore body profiles in the coupling results forms an ore body, thus obtaining the three-dimensional ore body distribution model.

8. The mineral reserve prediction method according to claim 7, characterized in that, After connecting adjacent orebody profiles in the coupling results to form an orebody, thus obtaining the three-dimensional orebody distribution model, the process further includes: Detect whether there are openings or intersecting triangles in the ore body distribution model; If open or intersecting triangles exist, adjust the local triangulation of the ore body to eliminate the open or intersecting triangles.

9. The mineral reserve prediction method according to claim 7, characterized in that, The process of estimating the grade distribution within the ore deposit in the 3D ore body model and calculating reserves based on the grade values ​​of the mineral assemblage samples and the grade values ​​of some basic blocks of the 3D ore body distribution model, to obtain a 3D reserve model, includes: Based on the statistical distribution of sample lengths in the sample files of the pre-established geological engineering database, the sample length of the mineral assemblage sample is determined. Grade combinations are performed based on the sample length, and the grade value of the mineral combination sample is calculated. Establish the minimum bounding box of the ore body in the aforementioned three-dimensional ore body distribution model; Add constraints to the minimum outer package to obtain the block model; The block model is divided into multiple basic blocks. The length of each basic block is equal to the sample length of the mineral assemblage sample. Each basic block includes internal blocks and boundary blocks. The internal blocks have the same size, the boundary blocks have the same size, and the size of the boundary blocks is smaller than that of the internal blocks.

10. The mineral reserves prediction method according to claim 1, characterized in that, The determination of the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model includes: Based on historical data of the mining area, the mineral processing cost, mining cost, mineral processing recovery rate, and mining recovery rate of the mining area are obtained. When the break-even point is reached, a balance relationship between cost and profit is established. Based on the relationship between the cost and the beneficiation cost, mining cost, beneficiation recovery rate, mining recovery rate, and cutoff grade of the mining area, and the relationship between the profit and the beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area, a balance relationship between the beneficiation cost, mining cost, beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area is further obtained. Based on the balance between the ore beneficiation cost, mining cost, ore beneficiation recovery rate, mining recovery rate, cutoff grade, and market price of the mining area, the price-cutoff grade calculation model is obtained. The current mineral market price is input into the price-boundary grade calculation model to obtain the boundary grade.

11. The mineral reserves prediction method according to claim 2, characterized in that, Before establishing a three-dimensional ore body distribution model based on the digital terrain model and profile lines, the following steps are also included: Drilling data is obtained through exploration drilling to establish a geological engineering database. The drilling data includes borehole file data and sample file data.

12. A mineral reserve prediction device, characterized in that, include: The reserve calculation unit is used to calculate reserves based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves of each grade interval within the pre-divided elevation interval of the mining area. The following iterative operations are performed until the grade and reserve values ​​of all basic blocks are obtained, resulting in a three-dimensional reserve model: Create a search ellipsoid based on the preset search radius; Iterate through all the underlying blocks to be valued and perform the following operations: Move the search ellipsoid so that the current base block to be valued is located at the center of the search ellipsoid; Search for whether there are pre-valued basic blocks and / or mineral combination samples within the space defined by the search ellipsoid; If there are already valued basic blocks and the mineral combination sample within the space defined by the search ellipsoid, estimate the grade value of the current basic block to be valued based on the grade value of the mineral combination sample and the grade value of the already valued basic blocks. The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued. If only pre-valued basic blocks exist within the space defined by the search ellipsoid, estimate the quality value of the current basic block to be valued based on the quality value of the pre-valued basic blocks. The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued. If only the mineral assemblage sample exists within the space defined by the search ellipsoid, the grade value of the current base block to be valued is estimated based on the grade value of the mineral assemblage sample. The reserve value of the base block is calculated based on the volume and grade value of the base block to be valued. After the traversal is completed, save the basic block grade value and reserve value obtained in this estimation, and expand the search radius of the search ellipsoid; The grade calculation unit is used to determine the boundary grade based on the current mineral market price and a pre-established price-boundary grade calculation model. The reserve calculation unit is also used to remove the portion of the initial mineral reserves that does not reach the boundary grade, so as to obtain the predicted mineral reserves. The step of calculating reserves based on a pre-established three-dimensional reserve model to obtain the initial mineral reserves for each grade interval within the pre-divided elevation interval of the mining area includes: Based on the mapping relationship between the spatial location and elevation interval of all basic blocks, and the mapping relationship between the grade value and grade interval of all basic blocks, the sample blocks contained in each grade interval within the elevation interval are obtained. The mineral reserves of each grade interval within the elevation interval are obtained by summing the reserve values ​​of the sample blocks contained in each grade interval within the elevation interval. Specifically, the grade range has an upper limit grade and a lower limit grade. For grade ranges where the upper limit grade does not reach the boundary grade, the reserves within that grade range are directly removed. For grade ranges where the lower limit grade does not reach the boundary grade but the upper limit grade does reach the boundary grade, the reserves within that grade range are removed, and the boundary grade is used as the new lower limit grade, while the upper limit grade remains unchanged, resulting in a new grade range. The reserves are then recalculated based on the three-dimensional reserve model to obtain the mineral reserves within the pre-divided elevation range of the mining area for that grade range.

13. A computer electronic device, comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the mineral reserves prediction method according to any one of claims 1 to 11.

14. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the mineral reserve prediction method according to any one of claims 1 to 11.

15. A computer program product comprising a computer program / instructions, characterized in that, When the computer program / instructions are executed by the processor, they implement the steps of the mineral reserve prediction method according to any one of claims 1 to 11.

Images