Method for dynamically evaluating resource consumption of open-pit mine by high-resolution satellite remote sensing

Abstract

This application discloses a method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing, relating to the field of mineral resource assessment and management. The method includes: acquiring stereo image pairs of the open-pit mine using high-resolution satellite remote sensing technology; the stereo image pairs include a front view image and a rear view image; geometrically locating the stereo image pairs; completing stereo image pair modeling based on the geometric location results to obtain an oblique photogrammetry 3D modeling model; performing 3D modeling of the ore body in the open-pit mine to obtain an ore body 3D modeling model; registering and aligning the ore body 3D modeling model and the oblique photogrammetry 3D modeling model; determining the quota earthwork volume based on the registration and aligning results; and assessing the resource utilization of the open-pit mine based on the quota earthwork volume. The resource utilization includes the amount of ore and the amount of metal utilized. This application, based on high-resolution satellite remote sensing technology, enables real-time and effective measurement of the annual resource utilization of open-pit mines.

Description

Technical Field

[0001] This application relates to the field of mineral resource assessment and management, and in particular to a method for dynamic assessment of the amount of resources used in open-pit mines using high-resolution satellite remote sensing. Background Technology

[0002] The annual resource utilization of open-pit mines is mainly assessed based on data reported by the mines. However, the accuracy of this method depends on the reported data. The three-dimensional data obtained by UAV oblique photography cannot be used to measure the annual resource utilization of open-pit mines in real time. Therefore, it is necessary to use high-resolution satellites to assess the resource utilization of open-pit mines. Summary of the Invention

[0003] The purpose of this application is to provide a method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing, which can effectively measure the annual resource utilization of open-pit mines in real time based on high-resolution satellite remote sensing technology.

[0004] To achieve the above objectives, this application provides the following solution: This application provides a method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing, including: Stereo image pairs of open-pit mines are acquired using high-resolution satellite remote sensing technology; the stereo image pairs include a front view image and a rear view image. Geometric positioning of the stereo image pair; Based on the geometric positioning results, stereo image pair modeling is completed to obtain the oblique photogrammetry 3D modeling model; A three-dimensional model of the ore body in an open-pit mine is obtained. The 3D modeling model of the ore body and the 3D modeling model of oblique photography are registered and docked. Based on the registration and docking processing results, the quota earthwork volume is determined; Based on the quota of earthwork volume, assess the amount of resources utilized in the open-pit mine; the amount of resources utilized includes the amount of ore utilized and the amount of metal utilized.

[0005] Optionally, the geometric positioning method includes: a geometric positioning method based on a physical model and a geometric positioning method based on a rational function model.

[0006] Optionally, a three-dimensional model of the ore body in the open-pit mine is performed to obtain a three-dimensional model of the ore body, including: The pinch-out point of the ore body is determined using the first method; the first method is interpolation, finite inference, or infinite inference. Based on the pinch-out point of the ore body, a second method is used to perform three-dimensional modeling interpolation of the ore body to obtain a three-dimensional modeling model of the ore body; the second method is the natural critical point method, the inverse power distance method, or the Kriging method.

[0007] Optionally, the registration and docking process includes registration processing and docking processing; The registration process is as follows: ; in, ; ; ; ; ; In the formula, where, The coordinate system for the 3D model of oblique photogrammetry; These are the X-axis, Y-axis, and Z-axis coordinates in the coordinate system of the oblique photogrammetry 3D modeling model; These are translation parameters; These are the translation parameters for the X-axis, Y-axis, and Z-axis, respectively. A coordinate system for the 3D modeling of the ore body; The X, Y, and Z coordinates of the ore body in the 3D model coordinate system; Scale variation parameters; Represents the transformation matrix; It represents the Euler angle.

[0008] Optionally, the docking process includes: The first intersection line is calculated using the triangular mesh of the 3D modeling model based on oblique photogrammetry and the triangular mesh of the ore body 3D modeling model. The region enclosed by the intersection of the first line of intersection and the triangular mesh of the oblique photogrammetry 3D model is defined as the first region; The area enclosed by the intersection of the first intersection line and the triangular mesh of the 3D model of the ore body is the second region; When there is an overlapping area between the first and second regions, it is determined that the ore body of the open-pit mine intersects with the oblique photogrammetry 3D model; When there is no overlap between the first and second regions, it is determined that the ore body of the open-pit mine does not intersect with the oblique photogrammetry 3D model.

[0009] Optionally, after determining that the ore body of the open-pit mine intersects with the oblique photogrammetry 3D model, the method further includes: cutting out the intersecting ore body portion in the ore body 3D model; The intersecting ore body portions in the 3D model of the ore body are cut out, specifically including: The intersection line between the ore body of the open-pit mine and the oblique photogrammetry 3D model is obtained as the second intersection line; The surface within the second intersection line in the oblique photogrammetry 3D model is determined as the top surface of the ore body, and modeling is performed accordingly; Obtain the closed polygon formed by the second intersection line and the triangular mesh of the oblique photogrammetry 3D modeling model; Construct a triangular mesh for the vertices of a closed polygon to obtain the polygon vertex set; Take three points from the set of polygon vertices in turn to determine if they are collinear. If they are collinear, delete the middle point. Traverse the set of polygon vertices to obtain the optimized set of polygon vertices; The concavity and convexity of each vertex in the polygon vertex set are calculated using normal vectors; Based on the concavity and convexity of each vertex in the polygon vertex set, a TIN-constructed surface is used as the top surface of the 3D model of the ore body after docking.

[0010] Optionally, based on the registration and docking processing results, the quota earthwork volume is determined, including: The top surface of the 3D model of the ore body obtained after the first docking moment is taken as the first top surface; The top surface of the ore body 3D model obtained after the second docking moment is the second top surface; the second moment is later than the first moment. A three-dimensional model is constructed using the first top surface as the top surface, the second top surface as the bottom surface, and the side surface of the open-pit mine ore body as the side surface. Divide the side of the 3D model into multiple movable blocks along the direction of the 3D side surface; Determine any mobilization block as the current mobilization block; Determine the upper and lower cross-sectional areas of the currently used block based on the cross-sectional area formula; Based on the current block's guide, upper cross-sectional area, and lower cross-sectional area, determine the volume of the current block. The sum of the volumes of all the blocks used is determined as the quota earthwork volume from the first time point to the second time point.

[0011] Optionally, the quota earthwork volume is: ; ; ; in, For quota earthwork volume; This represents the i-th block to be used; n is the number of blocks to be used. The distance between two sections in the i-th active block; m is the area of ​​the upper cross section; m is the number of triangular meshes in the cross section; Let be the area of ​​the i-th triangular mesh in the upper section; These are the distances to the three sides of the triangular mesh on the upper section; The cross-sectional area is the lower section area. Let be the area of ​​the i-th triangular mesh in the lower section; These are the distances to the three sides of the triangular mesh in the lower section.

[0012] Optionally, the amount of ore used is: ; in, To determine the amount of ore to be used; For the first Density of individual active blocks; The amount of metal used is: ; in, To utilize the amount of metal; To utilize the grade of the block segment.

[0013] Optionally, after assessing the amount of resources utilized by the open-pit mine based on the quota earthwork volume, the method further includes: Return to the step "Use high-resolution satellite remote sensing technology to acquire stereo image pairs of open-pit mines" until the absolute value of the difference between the amount of ore used and the quota of earthwork is less than the preset threshold. The amount of resources utilized in the open-pit mine is used as a vector input, which is then overlaid with the triangular mesh of the 3D model of the ore body to output a GIS vector-raster overlay image.

[0014] According to the specific embodiments provided in this application, the following technical effects are disclosed: This application provides a method for dynamically assessing the resource utilization of open-pit mines using high-resolution satellite remote sensing. It utilizes high-resolution satellite remote sensing to investigate the annual resource utilization of open-pit mines. For cases where annual UAV oblique photography is unavailable, the method acquires the three-dimensional spatial information of the open-pit mine annually through high-resolution satellite remote sensing. Then, by comparing the annual earthwork volume of the open-pit mine, the resource utilization is calculated using a block method. The core technology lies in the integration of the three-dimensional data from high-resolution satellite remote sensing with the utilized ore body. This application utilizes high-resolution satellite remote sensing technology to investigate the annual resource utilization of open-pit mines, providing a new method for investigating annual resource utilization in open-pit mines. It primarily addresses the challenge of real-time investigation of annual resource utilization in open-pit mines, offering a rapid new method for assessing annual mineral resource consumption. This method is of great significance for the investigation and evaluation of mineral resource utilization. Attached Figure Description

[0015] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments 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.

[0016] Figure 1 This is a schematic diagram illustrating the principle of a method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing, as described in one embodiment of this application. Figure 2 This is a flowchart of a method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing, as described in one embodiment of this application. Figure 3 This is a schematic diagram of satellite stereo image pair acquisition in one embodiment of this application; Figure 4 This is a schematic diagram of calibration points and ranging control points in one embodiment of this application; Figure 5 This is an intersection and inclusion model in one embodiment of this application; Figure 6 This is a schematic diagram of the original stereo image pair in one embodiment of this application; Figure 7 This is a stereoscopic image of a mining area overlay in one embodiment of this application. Figure 8 This is a schematic diagram of the GCP points for geometric correction in one embodiment of this application; Figure 9 This is a schematic diagram of a geometrically corrected stereo image pair in one embodiment of this application; Figure 10 This is a schematic diagram of a three-dimensional relative model in one embodiment of this application; Figure 11 This is a schematic diagram of a three-dimensional model of the ore body in one embodiment of this application; Figure 12 This is a schematic diagram of the three-dimensional docking of the ore body and the stereo image pair in one embodiment of this application; Figure 13 This is a schematic diagram illustrating the earthwork calculation principle in one embodiment of this application; Figure 14 This is the final drawing in one embodiment of this application. Detailed Implementation

[0017] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0018] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, the application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0019] In one exemplary embodiment, such as Figure 1As shown, a method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing is provided, including steps such as stereo image pair data acquisition (step 201), geometric positioning (step 202), stereo image pair modeling (step 203), 3D docking of ore body (steps 205 and 206), earthwork quota calculation (step 206), resource utilization estimation (step 207), and final mapping (step 208).

[0020] like Figure 2 As shown in this embodiment, the method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing specifically includes: Step 201: Acquire stereo image pairs of the open-pit mine using high-resolution satellite remote sensing technology. The stereo image pairs include a front view image and a rear view image.

[0021] The imaging principle of satellite stereo image pairs is that when a satellite is photographing a ground target, it simultaneously acquires forward-looking and rear-looking images. Stereo image pair products include one with a low shooting angle and one with a high shooting angle. For example... Figure 3 Stereo image pairs have a certain degree of overlap, satellite data has a certain degree of instability, the viewing angle of satellite images is limited by the orbit and has dialysis errors, satellite orbit time exposure is inconsistent and there are potential photometric differences, and the radiometric quality is poor, especially in shadowed areas. Therefore, stereo image pairs with good lighting, no clouds and few shadows are usually required.

[0022] Step 202: Perform geometric localization on the stereo image. Geometric localization methods include: physical model-based geometric localization methods and rational function model-based geometric localization methods. One method can be chosen, or both can be used. When using both, the physical model takes precedence, followed by the rational function model.

[0023] There are two methods for geometric positioning: one is based on physical model correction, and the other is based on rational function model. This method is mainly for non-public stereo image pairs such as satellite orbit parameters and sensor parameters at imaging time.

[0024] (1) Physical model.

[0025] Single-line array CCD pushbroom satellite remote sensing acquires two-dimensional images sequentially, row by row, along the vertical direction of the orbit. It then advances forward along the flight direction, scanning line by line to ultimately obtain an image. Each row of pixels in the image is projected onto the row center simultaneously. The image is tilted forward and backward at certain angles along the flight direction, forming back-view and forward-view images. When performing stereo imaging with single-line array CCD pushbroom remote sensing, the corresponding forward-view and back-view image plane coordinate systems can be derived from the image plane coordinates Ox of the downward-view image. n y n z nObtained through matrix translation and rotation. The flight direction tilt angle θ, and the coordinates of ground point p in the fore-and-aft images of a single-line CCD pushbroom remote sensing stereo image pair are (x...). F ,y F ,z F ) and (x B ,y B ,z B At the same time, the coordinates in the downward-looking image coordinates are (x N ,y N ,z N Then, we have: .

[0026] .

[0027] In three-line array CCD pushbroom remote sensing stereo imaging, when performing forward and backward viewing imaging, the corresponding forward and backward image coordinate systems can be considered as being derived from the image plane coordinates Ox of the downward viewing imaging. n y n z n The coordinates of p in the downward, forward, and backward viewing imaging coordinate systems are: ,and The transformation relationship is as follows: 。

[0028] .

[0029] The coordinates (X, Y, Z) of point P in the GIS coordinate system are different from the coordinates of the image point (x, y, z). p ,y p The coordinates (X, Y) of the CCD platform in the GIS coordinate system at the instant of imaging. s ,Y s Z s The camera's side tilt angle is Ψ. x ,Ψ y The camera's yaw, pitch, and roll angles are ω. Given that κ is the imaging scale factor M, the sensor model is: .

[0030] .

[0031] .

[0032] .

[0033] .

[0034] .

[0035] .

[0036] .

[0037] .

[0038] A rigorous physical sensor imaging model requires satellite orbit models, attitude parameters, imaging geometry models, etc., to describe the pushbroom remote sensing stereo imaging process.

[0039] (2) Rational geometric function model.

[0040] Because the satellite's orbital parameters, imaging timing, sensor parameters, and other related information need not be disclosed, a rigorous physical model cannot be used. Therefore, a rational function model is required, which is a geometric model of imaging. This model assumes that the pixel coordinates (r, c) represent the corresponding ground point, and the spatial coordinates (X, Y, Z) of the point p are polynomials, with the polynomial ratio as the independent variable. .

[0041] n is any pixel. It is a regular polynomial with a maximum degree of 3.

[0042] .

[0043] Where a0,…,a 19 The coefficients of a rational function, (r) n ,c n ) (X n ,Y n Z n () are the pixel coordinates and ground coordinates after translation and scaling, respectively.

[0044] .

[0045] X0,…,c0 are translation parameters, X s ,…,c s This is a standardized proportional parameter.

[0046] Ground control points utilize both calibration points and distance measurement control points to control accuracy. Calibration points and distance measurement control points are as follows: Figure 4 .

[0047] Step 203: Based on the geometric positioning results, complete the stereo image pair modeling to obtain the oblique photogrammetry 3D modeling model.

[0048] Calculate the GIS coordinates of the image points based on the coordinates of the two images: .

[0049] make .

[0050] .

[0051] 。

[0052] Substituting into Taylor's formula and expanding, we get: .

[0053] Error equation: .

[0054] The rational function models RFM1 and RFM2 for stereo image pairs have corresponding model parameters RPC1 and RPC2. The row and column coordinates of the feature points in the stereo image pair pixels are (ROW1, COL1) and (ROW2, COL2), respectively. The standardized coordinates are (r1, c1) and (r2, c2). The error equation is as follows: .

[0055] Where v = Ad - l, then the ground point coordinate correction is: .

[0056] Ground control points are set to correct the budget value, and a maximum threshold and a maximum number of iterations are set. The calculation stops when the correction value is less than the threshold or the maximum number of iterations is reached, and an approximate numerical solution of the ground point coordinates is obtained. The smaller the threshold, the larger the number of iterations and the longer the calculation time.

[0057] Step 204: Perform 3D modeling of the ore body in the open-pit mine to obtain a 3D model of the ore body.

[0058] A 3D model of the ore body in an open-pit mine is generated, including: determining the pinch-out point of the ore body using a first method. The first method can be interpolation, finite inference, or infinite inference. Based on the pinch-out point, a second method is used for 3D modeling interpolation of the ore body to obtain the 3D model. The second method can be the natural critical point method, the inverse power distance method, or the Kriging method.

[0059] Step 205: Register and connect the 3D model of the ore body and the 3D model of the oblique photogrammetry.

[0060] Registration and docking processing includes registration processing and docking processing.

[0061] The registration process is as follows: .

[0062] in, .

[0063] .

[0064] .

[0065] .

[0066] .

[0067] In the formula, where, The coordinate system for the 3D model of oblique photogrammetry. These are the X-axis, Y-axis, and Z-axis coordinates in the coordinate system of the oblique photogrammetry 3D model. These are the translation parameters. These are the translation parameters for the X-axis, Y-axis, and Z-axis, respectively. The coordinate system for the 3D modeling of the ore body. The coordinates of the X, Y, and Z axes in the coordinate system of the 3D model of the ore body. Scale variation parameters. This represents the transformation matrix. It represents the Euler angle.

[0068] The docking process includes: calculating the first intersection line using the triangular meshes of the oblique photogrammetry 3D model and the ore body 3D model; defining the area enclosed by the intersection points of the first intersection line and the triangular mesh of the oblique photogrammetry 3D model as the first region; and defining the area enclosed by the intersection points of the first intersection line and the triangular mesh of the ore body 3D model as the second region. If the first and second regions overlap, the ore body of the open-pit mine is determined to intersect with the oblique photogrammetry 3D model. If the first and second regions do not overlap, the ore body of the open-pit mine is determined not to intersect with the oblique photogrammetry 3D model.

[0069] After determining that the ore body of the open-pit mine intersects with the oblique photogrammetry 3D model, the process also includes: cutting out the intersecting ore body portion in the ore body 3D model.

[0070] The process involves cutting out the intersecting ore body portions in the 3D model of the ore body, specifically including: obtaining the intersection line between the ore body of the open-pit mine and the oblique photogrammetry 3D model as the second intersection line; determining the surface within the range of the second intersection line in the oblique photogrammetry 3D model as the top surface of the ore body and modeling it; obtaining the closed polygon formed by the second intersection line and the triangular mesh of the oblique photogrammetry 3D model; constructing triangular meshes for the vertices of the closed polygon to obtain the polygon vertex set; sequentially selecting three points from the polygon vertex set to determine if they are collinear, deleting the middle point if they are; traversing the polygon vertex set to obtain the optimized polygon vertex set; calculating the concavity / convexity of each vertex in the polygon vertex set using normal vectors; and constructing a surface using TIN as the top surface of the ore body 3D model after docking, based on the concavity / convexity of each vertex in the polygon vertex set.

[0071] 3D docking of ore bodies: To estimate the amount of ore resources to be utilized, it is first necessary to dock the ore body with the oblique photogrammetry 3D model. After aerial triangulation modeling, the obtained result is a 3D spatial coordinate in a coordinate system. Most mines use cross-sectional descriptions for ore bodies, which require the creation of a 3D model. The oblique photogrammetry real-scene 3D model and the ore body 3D model are visually surface features. After docking, there are various relationships such as intersection, overlap, and inclusion. These relationships need to be interpreted before remodeling.

[0072] ① Ore body modeling: First, determine the pinch-out point of the ore body, and then use interpolation, finite inference and infinite inference methods.

[0073] The interpolation method mainly obtains the boundary reference points through interpolation. The specific formula is as follows: .

[0074] In the formula, X is the distance from the borehole to the ore body boundary, and C min C represents the grade at the boundary of the ore body. A C B Let L represent the ore grade of industrially mined borehole A and non-industrially mined borehole B, and L be the distance between the two boreholes. If the ore body thickness in the non-mineralized borehole is less than the minimum minable thickness but the grade meets the requirements, then: .

[0075] In the formula, M A M B M represents the ore body thickness of the ore-bearing borehole A and the non-industrial ore-bearing borehole B. min This is the minimum mineable thickness.

[0076] ② Common interpolation methods for 3D modeling of ore bodies include the natural critical point method, the inverse power distance method, and the Kriging method.

[0077] The natural proximity electric field method uses the dual graph of a triangulation network to perform arbitrary point interpolation, as shown in the following formula: 。

[0078] In the formula, Z(Y) represents the information value of a known point; D is the distance between the point to be interpolated and its naturally neighboring known points; D0 is the critical influence distance; and n is the number of relevant known points. k (A) is the value of the interpolation point obtained separately for the k-th known point.

[0079] The inverse power method of distance assumes that the closer the known point and the interpolation point are, the greater the correlation. The main calculation formula is: .

[0080] D k denoted as , where is the distance from the interpolation point to the known point, and μ is the exponent.

[0081] Kriging minimizes the interpolation points using linear regression, based on covariance.

[0082] 。

[0083] Z k Given the value of ω at a point. k These are weights, which can be solved using a combination of functions. It represents the value of the distance model variable graph between points i and j.

[0084] .

[0085] ③ Registration of the ore body with the 3D image model.

[0086] First, the ore body and the image need to be transformed into the same coordinate system. The two projected coordinate systems (X1, Y1, Z1) are transformed into (X2, Y2, Z2). The origins O1 and O2 of the two systems are different, and the coordinate axes are not parallel. During the transformation, in addition to the translation parameters (ΔX, ΔY, ΔZ), the rotation parameters (α, β, γ) corresponding to the Euler angles (εX, εY, εZ) are also required. In addition, the scale transformation parameter κ needs to be set, for a total of 7 transformation parameters. The transformed coordinates are consistent with the coordinates of the 3D model.

[0087] ④ Processing for docking orebody and image 3D model: The docking of orebody 3D model with open-pit mine image 3D model mainly includes intersection, inclusion, and non-correlation. Open-pit mine activity monitoring typically involves intersection or non-correlation (e.g., ...). Figure 5 If the ore body is not fully mined, it is considered intersecting; if the ore body is fully mined, it is considered unrelated.

[0088] To determine if they intersect, the intersection line L is calculated using the triangular mesh T1 from the image model and the triangular mesh T2 from the ore body model. The intersection line L intersects with the triangular meshes T1 and T2 at points A, B, C, D, forming intervals [A, B] and [C, D]. If they overlap, they intersect; otherwise, they do not intersect.

[0089] For a ore body A that is not related to the 3D model of the open-pit mine image, there are n blocks. Each block has a volume Vi, density pi, grade wi, ore quantity dA, and metal quantity mA. The total resource quantity of the entire ore body is directly included in the utilization quantity.

[0090] .

[0091] For intersecting ore bodies, it is necessary to cut out the intersecting ore body portion and then perform calculations.

[0092] The intersection involves the line of intersection between the ore body and the ground of the 3D model image. The surface within the intersection line is used as the top surface of the ore body for modeling. The intersection line and the triangular mesh form a closed polygon, and triangular meshes are constructed for the vertices of the polygon. For the vertices of the formed polygon [V1,V2,V3,…,Vn], three points are taken sequentially to determine if they are collinear. If they are collinear, the intermediate points are deleted, and the above process is repeated until the end. The concavity and convexity of the vertices are determined using normal vectors. These triangular meshes are used to construct surfaces through TIN, which serve as the top surface of the docking 3D ore body model.

[0093] Step 206: Determine the quota earthwork volume based on the registration and docking processing results.

[0094] Step 206 includes: obtaining the top surface of the 3D model of the ore body after docking at the first moment as the first top surface; obtaining the top surface of the 3D model of the ore body after docking at the second moment as the second top surface. The second moment is later than the first moment. Constructing a 3D model using the first top surface as the top surface, the second top surface as the bottom surface, and the side surface of the open-pit mine ore body as the side surface. Dividing the side surface of the 3D model into multiple movable blocks along the direction of the 3D side surface. Determining any movable block as the current movable block. Determining the upper and lower cross-sectional areas of the current movable block according to the cross-sectional area formula. Determining the volume of the current movable block based on the current movable block's guide, upper and lower cross-sectional areas. Determining the sum of the volumes of all movable blocks as the quota earthwork volume from the first moment to the second moment.

[0095] The quota for earthwork volume is: ; ; ; in, For quota earthwork volume; This represents the i-th block to be used; n is the number of blocks to be used. The distance between two sections in the i-th active block; m is the area of ​​the upper cross section; m is the number of triangular meshes in the cross section; Let be the area of ​​the i-th triangular mesh in the upper section; These are the distances to the three sides of the triangular mesh on the upper section; The cross-sectional area is the lower section area. Let be the area of ​​the i-th triangular mesh in the lower section; These are the distances to the three sides of the triangular mesh in the lower section.

[0096] Open-pit mining involves not only the extraction of the ore body but also the surrounding rock. Therefore, the earthwork quota refers to the earthwork volume of the ore body, and the calculation mainly focuses on the earthwork volume T of the mined portion of the ore body. Assume the top surface constructed in the first year is S1, and the top surface constructed in the second year is S2. Using the second-year top surface as the bottom surface, the first-year top surface as the top surface, and the ore body side surface as the side surface, a volume is formed. This volume is the earthwork quota calculation quantity. During the calculation, the earthwork volume formed by S1, S2, and the ore body side surface is divided into n segments along a certain direction. Each segment has two cross-sections, D1 and D2, whose cross-sectional areas are calculated using a DTM network and then divided into m small triangular mesh areas.

[0097] Step 207: Assess the amount of resources to be utilized in the open-pit mine based on the quota earthwork volume. The amount of resources to be utilized includes the amount of ore and the amount of metal to be utilized.

[0098] The amount of ore used is: .

[0099] in, The amount of ore to be used. For the first The density of the active block segment is the volume of the ore body composed of S1 and S2.

[0100] The amount of metal used is: .

[0101] in, To utilize the amount of metal. To utilize the grade of the block segment.

[0102] ① Resource estimation should be conducted, and the estimation parameters, methods, and block divisions for all types of resources should be consistent with the exploration report. The industrial indicators for the deposit should be consistent with current industrial indicators; any inconsistencies between the industrial indicators in the report and current indicators should be corrected to reflect current indicators.

[0103] Basic parameters for resource estimation include the area, average thickness, dip angle, grade, and average weight of the ore body within the utilization range. Sometimes, ore moisture content and mineralization coefficient are also included. The most recent exploration report data should be used. Based on factors such as geological body occurrence, faults, and exploration projects, profile locations are selected. The vertical profile method is used to create the boundary lines of each profile, which are then projected onto the reserve estimation map. Connecting the boundary points of each profile line delineates the utilization range.

[0104] Depending on the scope of the impact, geometric methods (arithmetic mean method, geological block method, mining block method, cross-section method, contour line method, linear reserve method, trigonometric method, nearest area method, polygonal method), statistical analysis methods (distance-weighted method, kriging method), and SD method are used to calculate the amount of land covered. The geological block method and cross-section method are commonly used.

[0105] The basic formula for calculating ore quantity: The basic formula for calculating metal content: Quantity of ore used The amount of metal used is M, the total number of blocks used is n, and the first... Each moving block volume , No. Individual segment density The i-th segment quality w i .

[0106] ② Verification of resource utilization: Verify whether the estimated volume of ore to be utilized is consistent with the calculated volume of earthwork quota. If the difference between the two is less than the given limit ε, the calculation is considered valid; otherwise, start from scratch.

[0107] .

[0108] Step 208: Return to step 201 until the absolute value of the difference between the amount of ore used and the quota of earthwork is less than the preset threshold. Use the amount of resources used in the open-pit mine as a vector input, overlay it with the triangular mesh of the 3D model of the ore body, and output a GIS vector-raster overlay image.

[0109] 1. Vector and raster overlay device.

[0110] The optimized annual mobilization range vector is input, along with the bands R (red): 669.43 nm, G (green): 538.96 nm, and B (blue): 479.25 nm from the oblique photography. These are used as inputs for the RGB (red, green, blue) color composite image. The composite color image is generated using points, lines, and surfaces with the same projection. The raster and vector are then overlaid using coordinate layering.

[0111] vector x and y are the corresponding coordinates, and z is the eigenvalue. Vector values, raster , and For the corresponding coordinates, This represents the raster grayscale value.

[0112] make, Thus achieving a grid Grayscale values ​​and vectors The superposition of.

[0113] By using corresponding coordinates to overlay raster and vector data, the output is a GIS vector and raster overlay image, thus creating an image that is suitable for human vision.

[0114] 2. Final Image: Input the UAV oblique photogrammetry image and the annual usage range, and output a superimposed image (GIS) suitable for human visual observation. The final image can be output as a JPG or TIF format through software.

[0115] Below, with Figure 6 Using the satellite stereo image pair shown as an example, this invention provides a detailed explanation of the method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing.

[0116] Step 1: Obtain stereo image pair data. Figure 6 This is a 15-meter resolution satellite stereo image pair acquired on October 8, 2025. The resulting overlay image of the mining area's stereo image pair is as follows. Figure 7 .

[0117] Step 2: Geometric Positioning: Geometric positioning is performed using the GCP rational geometric function method; the geometrically corrected GCP points are as follows: Figure 8 Geometrically corrected stereo image pairs, such as Figure 9 As shown.

[0118] Step 3: Establishing the 3D model of the stereo pair. The stereo pair 3D model is as follows: Figure 10 As shown.

[0119] Step 4: Construct a 3D model of the ore body; the 3D model of the ore body is as follows: Figure 11 As shown.

[0120] Step 5: Perform 3D docking between the ore body and the stereo image pair. The docking result is as follows: Figure 12 As shown.

[0121] Step 6: Calculate earthwork and stonework. The calculation principle is as follows: Figure 13 As shown.

[0122] .

[0123] Step 7: Calculation of mine mobilization: .

[0124] Step 8: Final image, as shown Figure 14 As shown.

[0125] In one exemplary embodiment, a computer device is provided, which may be a server or a terminal. The computer device includes a processor, memory, input / output interfaces (I / O), and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is connected to the system bus via the I / O interfaces. The processor of the computer device provides computing and control capabilities. The memory of the computer device includes non-volatile storage media and internal memory. The non-volatile storage media stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The I / O interfaces of the computer device are used for exchanging information between the processor and external devices. The communication interface of the computer device is used for communicating with external terminals via a network connection.

[0126] In one exemplary embodiment, a computer device is also provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps in the above-described method embodiments.

[0127] In one exemplary embodiment, a computer-readable storage medium is provided storing a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0128] In one exemplary embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0129] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0130] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments described above. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0131] The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0132] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0133] This document uses specific examples to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. Furthermore, those skilled in the art will recognize that, based on the ideas of this application, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A high-resolution satellite remote sensing method for dynamically evaluating the amount of resources used in an open-pit mine, characterized by, include: Use high-resolution satellite remote sensing technology to acquire stereo image pairs of open-pit mines; The stereo image pair includes a front view image and a rear view image; Geometric positioning of the stereo image pair; Based on the geometric positioning results, stereo image pair modeling is completed to obtain the oblique photogrammetry 3D modeling model; A three-dimensional model of the ore body in an open-pit mine is obtained. The 3D modeling model of the ore body and the 3D modeling model of oblique photography are registered and docked. Based on the registration and docking processing results, the quota earthwork volume is determined; Based on the quota of earthwork volume, assess the amount of resources utilized in the open-pit mine; the amount of resources utilized includes the amount of ore utilized and the amount of metal utilized.

2. The method according to claim 1, characterized in that, The geometric positioning methods include: a geometric positioning method based on a physical model and a geometric positioning method based on a rational function model.

3. The method according to claim 1, characterized in that, A 3D model of the ore body in an open-pit mine is created, resulting in a 3D model of the ore body, including: The pinch-out point of the ore body is determined using the first method; the first method is interpolation, finite inference, or infinite inference. Based on the pinch-out point of the ore body, a second method is used to perform three-dimensional modeling interpolation of the ore body to obtain a three-dimensional modeling model of the ore body; the second method is the natural critical point method, the inverse power distance method, or the Kriging method.

4. The method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing according to claim 1, characterized in that, The registration and docking process includes registration processing and docking processing; The registration process is as follows: ； in, ； ； ； ； ； In the formula, where, The coordinate system for the 3D model of oblique photogrammetry; These are the X-axis, Y-axis, and Z-axis coordinates in the coordinate system of the oblique photogrammetry 3D modeling model; These are translation parameters; These are the translation parameters for the X-axis, Y-axis, and Z-axis, respectively. A coordinate system for the 3D modeling of the ore body; The X, Y, and Z coordinates of the ore body in the 3D model coordinate system; Scale variation parameters; Represents the transformation matrix; It represents the Euler angle.

5. The method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing according to claim 4, characterized in that, The docking process includes: The first intersection line is calculated using the triangular mesh of the 3D modeling model based on oblique photogrammetry and the triangular mesh of the ore body 3D modeling model. The region enclosed by the intersection of the first line of intersection and the triangular mesh of the oblique photogrammetry 3D model is defined as the first region; The area enclosed by the intersection of the first intersection line and the triangular mesh of the 3D model of the ore body is the second region; When there is an overlapping area between the first and second regions, it is determined that the ore body of the open-pit mine intersects with the oblique photogrammetry 3D model; When there is no overlap between the first and second regions, it is determined that the ore body of the open-pit mine does not intersect with the oblique photogrammetry 3D model.

6. The method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing according to claim 4, characterized in that, After determining that the ore body of the open-pit mine intersects with the oblique photogrammetry 3D model, the process also includes: cutting out the intersecting ore body portion in the ore body 3D model; The intersecting ore body portions in the 3D model of the ore body are cut out, specifically including: The intersection line between the ore body of the open-pit mine and the oblique photogrammetry 3D model is obtained as the second intersection line; The surface within the second intersection line in the oblique photogrammetry 3D model is determined as the top surface of the ore body, and modeling is performed accordingly; Obtain the closed polygon formed by the second intersection line and the triangular mesh of the oblique photogrammetry 3D modeling model; Construct a triangular mesh for the vertices of a closed polygon to obtain the polygon vertex set; Take three points from the set of polygon vertices in turn to determine if they are collinear. If they are collinear, delete the middle point. Traverse the set of polygon vertices to obtain the optimized set of polygon vertices; The concavity and convexity of each vertex in the polygon vertex set are calculated using normal vectors; Based on the concavity and convexity of each vertex in the polygon vertex set, a TIN-constructed surface is used as the top surface of the 3D model of the ore body after docking.

7. The method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing according to claim 1, characterized in that, Based on the registration and docking processing results, the quota earthwork volume is determined, including: The top surface of the 3D model of the ore body obtained after the first docking moment is taken as the first top surface; The top surface of the ore body 3D model obtained after the second docking moment is the second top surface; the second moment is later than the first moment. A three-dimensional model is constructed using the first top surface as the top surface, the second top surface as the bottom surface, and the side surface of the open-pit mine ore body as the side surface. Divide the side of the 3D model into multiple movable blocks along the direction of the 3D side surface; Determine any mobilization block as the current mobilization block; Determine the upper and lower cross-sectional areas of the currently used block based on the cross-sectional area formula; Based on the current block's guide, upper cross-sectional area, and lower cross-sectional area, determine the volume of the current block. The sum of the volumes of all the blocks used is determined as the quota earthwork volume from the first time point to the second time point.

8. The method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing according to claim 7, characterized in that, The quota for earthwork volume is: ； ； ； in, For quota earthwork volume; This represents the i-th block to be used; n is the number of blocks to be used. The distance between two sections in the i-th active block; m is the area of ​​the upper cross section; m is the number of triangular meshes in the cross section; Let be the area of ​​the i-th triangular mesh in the upper section; These are the distances to the three sides of the triangular mesh on the upper section; The cross-sectional area is the lower section area. Let be the area of ​​the i-th triangular mesh in the lower section; These are the distances to the three sides of the triangular mesh in the lower section.

9. The method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing according to claim 7, characterized in that, The amount of ore used is: ; in, To determine the amount of ore to be used; For the first Density of individual active blocks; The amount of metal used is: ; in, To utilize the amount of metal; To utilize the grade of the block segment.

10. The method for dynamic assessment of resource utilization in open-pit mines using high-resolution satellite remote sensing according to claim 1, characterized in that, After assessing the amount of resources utilized by the open-pit mine based on the quota earthwork volume, the process also includes: Return to the step "Use high-resolution satellite remote sensing technology to obtain stereo image pairs of open-pit mines" until the absolute value of the difference between the amount of ore used and the quota of earthwork is less than the preset threshold. The amount of resources utilized in the open-pit mine is used as a vector input, which is then overlaid with the triangular mesh of the 3D model of the ore body to output a GIS vector-raster overlay image.

Images