A mineral resource overburden query method and system based on a three-dimensional model of an ore body

Abstract

The present application belongs to the technical field of mineral resources information processing, and discloses a mineral resources overburden query method and system based on a three-dimensional model of ore bodies. The method comprises the following steps: collecting all mine area geology and mineral resources data in the query area; performing three-dimensional modeling on all ore bodies in the query area; collecting DEM data of the query area to establish a surface model of the query area; obtaining a three-dimensional model of the maximum protection range of the ore body according to the occurrence of the three-dimensional model of the ore body; extending the three-dimensional model of the maximum protection range of the ore body to the surface model surface; determining the containment zone range of the queried construction project and falling onto the same surface model surface; and obtaining the state of the overburden mineral resources in the query area caused by the queried construction project according to the spatial position relationship between the three-dimensional model of the maximum protection range of the ore body and the containment zone range of the queried construction project. The present application can quickly and intuitively determine the state of the overburden mineral resources of the construction project land.

Description

Technical Field

[0001] This invention belongs to the field of mineral resource information processing technology, and in particular relates to a method and system for querying mineral resource overburden based on a three-dimensional model of an ore body. Background Technology

[0002] In the scientific guidance of mineral resource management, the inquiry of mineral resources covered by construction projects, the approval of construction projects, and the improvement of the protection and management of proven mineral resources reserves, the inquiry area is used as the whole inquiry unit to inquire about the status of mineral resources covered by surface and underground construction project land in three-dimensional space, and to qualitatively determine the coverage status of the mineral resources in the mining area (divided into full coverage, partial coverage, and no coverage), so as to provide basic information for mineral resource management and construction project approval.

[0003] Currently, the commonly used method for investigating mineral resources overlaid by construction projects is the traditional geological survey and assessment method. This method is technically cumbersome and complex, the results are not intuitive, and the calculation process is prone to errors that are not easily detected, leading to deviations in the survey results. Traditional methods also involve long survey areas, which is not conducive to the efficient work of natural resource management departments under the new circumstances. Summary of the Invention

[0004] To overcome the problems existing in related technologies, the present invention discloses a method and system for querying mineral resource overburden based on a three-dimensional model of the ore body.

[0005] The technical solution is as follows: a method for querying mineral resource overburden based on a three-dimensional ore body model, comprising the following steps:

[0006] S1 collects and queries geological and mineral data for all mining areas within the region;

[0007] S2, performs 3D modeling of all ore bodies in the query area;

[0008] S3, collect DEM data of the query area and build a surface model of the query area;

[0009] S4. Obtain the three-dimensional model of the maximum protection range of the ore body based on the occurrence of the three-dimensional model of the ore body.

[0010] S5 extends the three-dimensional model of the maximum protection range of the ore body to the surface model;

[0011] S6, determine the scope of the retaining wall of the queried construction project and place it on the same surface model. If the construction project is an underground project, construct a three-dimensional model of the retaining wall of the construction project and import it into the same surface model space.

[0012] S7. Establish a query system to obtain the status of mineral resources covered by the project within the query area based on the spatial relationship between the three-dimensional model of the maximum protection range of the ore body and the protection zone of the project being queried.

[0013] In step S2, three-dimensional modeling is performed on all ore bodies in the query area, including: three-dimensional modeling of different ore bodies in the same mining area using one of the following methods: profile method, surface outcrop extrapolation method, top / bottom plate contour line method, and borehole data modeling method.

[0014] In step S3, a surface model of the query area is established using ArcGIS or 3Dmine software via the DEM method. This specifically includes:

[0015] S3.1, Collect DEM digital elevation information data within the query area;

[0016] S3.2, Extract DEM elevation point cloud data within the query area;

[0017] S3.3, use 3DMine or ArcGIS software to import point clouds and build a surface model.

[0018] In step S4, the maximum protection range three-dimensional model of the ore body is obtained based on the occurrence of the ore body three-dimensional model, including: loading the three-dimensional models of all ore bodies in a mining area into the three-dimensional mining modeling software, constructing the maximum protection range three-dimensional model for each ore body, and using the software entity module-Boolean calculation to combine the maximum protection range three-dimensional models of all ore bodies to obtain the maximum protection range three-dimensional model of the ore body in the mining area.

[0019] Furthermore, using the software entity module - Boolean calculation, the three-dimensional models of the maximum protection range of all ore bodies are combined to obtain the three-dimensional model of the maximum protection range of the ore bodies in the mining area, including:

[0020] 1) Enable editing mode for a single ore body model in the 3D mining modeling software;

[0021] 2) Determine the overall strike and dip direction of the ore body;

[0022] 3) By using geological data of the mining area, obtain characteristic parameters of lithology, bedding, joints, fracture coefficient, and metamorphic degree of the ore body and surrounding rocks, and determine the strata movement angle of the mining area;

[0023] 4) Around the ore body, draw line segments in four directions according to the dip direction and the angle of the ore body's uphill movement. The length of the line segments should be from the lowest point to the lowest elevation of the ore body and from the highest point to above the ground surface.

[0024] 5) Copy the four line segments to the four different locations of the ore body, and paste them one by one to the tangent points of all the turning points and protrusions of the ore body in the corresponding locations;

[0025] 6) Check the length of the tangent segments around the ore body to see if the lowest point is the lowest elevation of the ore body and the highest point is above the ground surface;

[0026] 7) In the 3D mining software, select the snap point to draw a polyline, and draw two closed lines at the lowest and highest points of all tangent segments around the ore body;

[0027] 8) Close the ore body model; only the tangent segments around the ore body and the closed lines at the lowest and highest points are displayed in the software interface.

[0028] 9) Generate triangular meshes within and between the two closed lines, merge all triangular meshes, generate a solid model, and verify the solid model. This solid model is the three-dimensional model of the maximum protection range of the single ore body.

[0029] In step S5, the three-dimensional model of the maximum protection range of the ore body is extended to the surface model surface, including: loading the three-dimensional model of the maximum protection range of the ore body of all mining areas and the surface model of the query area into the same software interface, using the surface clipping function to clip and delete the three-dimensional model of the maximum protection range of the ore body that is higher than the surface range, and obtaining the three-dimensional model of the maximum protection range of the ore body below the surface.

[0030] In step S7, the status of mineral resources covered by the queried construction project within the query area is obtained, including: overlaying the physical model of the delineated surface or underground retaining zone of the construction project with the three-dimensional model of the maximum protection range of the ore body in the mining area in three-dimensional space, and analyzing and determining the status of ore body coverage in the mining area; the specific determination steps are as follows:

[0031] S7.1, Covering a single mining area: If the maximum protection range 3D model of a single mining area within the query area intersects with the scope of the project's retaining wall or the entity model of the retaining wall, then the project is determined to cover a single mining area within the query area.

[0032] S7.2, Covering Multiple Mining Areas: If the maximum protection range 3D model of two or more mining areas within the query area intersects with the scope or entity model of the project's retaining wall, then the project is determined to cover multiple mining areas within the query area.

[0033] S7.3, No Overburden: If the maximum protection range 3D model of the ore body in the query area does not intersect with the scope of the project's retaining belt or the entity model of the retaining belt, then the project is determined to not overburden the ore body in the query area.

[0034] Step S7.2 specifically includes:

[0035] Step 1: Decompose the non-intersecting areas and intersecting areas of the three-dimensional models of the maximum protection range of the ore bodies of two or more mining areas within the query area with the scope of the project's retaining belt or the entity model of the retaining belt, and obtain different regional components;

[0036] Step 2: Use image grayscale difference spectrum analysis to remove outliers from each intrinsic region component, and calculate the multi-size arrangement weights of each region component to construct a feature vector;

[0037] Step 3: Use intersection region analysis to reduce the order of feature vectors, and use a support vector machine recognizer to identify all or part of the three-dimensional models of the maximum protection range of ore bodies in two or more mining areas with the range of the project's retaining belt or a single sample of the retaining belt entity model.

[0038] New feature vectors are formed by linear projection of the feature vectors, and the intersection region of the features is calculated. The expression is:

[0039]

[0040] In the formula, Φ represents the feature of the intersection region, and v k Let covariance matrix be the variance matrix. Features that require downgrading in order to represent all or part of the 3D model of the maximum protection range of ore bodies in two or more mining areas compared to the scope of the retaining belt or the solid model of the retaining belt of the construction project. The feature mean of the training samples

[0041] The value of the feature Φ of the intersection area is determined according to the principle of choosing the higher protection level. If it is a low protection level, it is a mining area covered by mining; if it is a high protection level, it is a mining area not covered by mining.

[0042] The covariance matrix v k The calculation formula is as follows:

[0043]

[0044] In the formula, X is the total number of three-dimensional models of the maximum protection range of ore bodies in two or more mining areas that are all or part of the project's retaining zone or retaining zone entity model, x is any single sample among the three-dimensional models of the maximum protection range of ore bodies in two or more mining areas that are all or part of the project's retaining zone or retaining zone entity model, and K is the exponent.

[0045] Furthermore, for two or more mining areas, the three-dimensional models of the maximum protection range of the ore bodies are all or partly intersected with the project's retaining wall area or retaining wall entity model X(t) by performing a non-overlapping region decomposition, expressed as:

[0046]

[0047] In the formula, y x(t) represents the x-th regional component, where X is the total number of three-dimensional models of the maximum protection range of ore bodies in two or more mining areas, including all or part of the model and the project's retaining wall range or retaining wall entity model, i.e., the total number of regions to be decomposed, and f(t) is the remainder; in the non-overlapping region decomposition method, y is obtained through decomposition. x The (t) component and the residual component f(t) can completely recover the original signal X(t).

[0048] In step two, the multi-size arrangement weights of each regional component are calculated to construct a feature vector, including:

[0049] Inherent region components A one-dimensional size sequence of length M The sequence was obtained by coarsening:

[0050]

[0051] In the formula, v is the scaling factor, Q = 1, 2, ... [M / χ], and [M / χ] represents the integer division of M / χ; when χ = 1, the coarse-grained sequence is the original sequence; M is the intrinsic region component. The length of is given by E, where E is a one-dimensional size sequence and i is a certain length of an intrinsic region component. For coarse-grained sequences, Qχ is the integer value of a certain scale factor, C is the length value of a certain scale factor after linear calculation, i is the linear calculation coefficient of the scale factor, and Q is the integer value of M / χ.

[0052] For coarse-grained sequences After reconstructing the dimensions, we get:

[0053]

[0054] In the formula, For the sequence reconstructed from coarse-grained sequence size, This is the initial sequence for reconstructing the coarse-grained sequence size. For reconstructing a sequence with a delayed size after initialization from a coarse-grained sequence size, For the reconstructed coarse-grained sequence size, the sequence is delayed by n-1 embedding dimensions after initialization, where ε is the delay scale factor bias. The deviation value in the lag scale factor deviation is the deviation value after lag of m-1 dimensions, where m is the dimension of the size reconstruction; r is the r-th reconstruction component, r = 1, 2…M-(n-1), where n is the embedding dimension, and Γ is the lag size; r1, r2…r n Represents reconstructed components The index of the column containing each element will Arranged in ascending order as follows:

[0055]

[0056] In the formula, To be Reconstructed components arranged in ascending order The first delayed index sequence of each element in the column. To be Reconstructed components arranged in ascending order The delayed index sequence of each element in the column. To be Reconstructed components arranged in ascending order The nth delayed index sequence of each element in the column;

[0057] If equal values ​​exist in the reconstructed components, they are arranged in chronological order; for any coarse-grained sequence We obtain a sequence of symbols B(l) = [c1, c2, ..., c1] n ], where l = 1, 2, ..., N, and N ≤ n! ; the reconstructed sequence with an embedding dimension of n has n! permutations, and the symbol sequence B(l) is one of these permutations. Calculate the probability O of each symbol sequence. l (l = 1, 2, ..., N); B(l) represents the value for any coarse-grained sequence. A sequence of symbols, c, is obtained. n Let l be the sequence of symbols, N be the number of symbols, n be the embedding dimension, and O be the number of symbols. l The probability of each symbol sequence occurring;

[0058] Define the permutation entropy R of different symbol sequences in the form of information entropy. C (n) is:

[0059]

[0060] When O r = 1 / n!, R C (n) Take the maximum value ln(n!), and let R C (n) After normalization, we get:

[0061] R C =R C (n) / ln(n!)

[0062] Among them, R C The normalized permutation weights, 0 <R C <1,R C The value can reflect and amplify minute changes in a size sequence;

[0063] The feature vector is constructed by calculating the multi-size permutation weights of each intrinsic region component:

[0064]

[0065] In the formula, ∧ represents the feature vector value constructed by the multi-size arrangement weights of each intrinsic region component. The permutation weights after multi-size normalization of the first intrinsic region component. The permutation weights of the Xth intrinsic region component after multi-size normalization are given. Reconstruct the permutation weights for the multi-size normalization process of the first intrinsic region component. The permutation weights are reconstructed using multi-size normalization processing for the Xth intrinsic region component.

[0066] Another objective of this invention is to provide a mineral resource overburden query system based on a three-dimensional ore body model. This system implements the aforementioned mineral resource overburden query method based on a three-dimensional ore body model. The system includes:

[0067] The data collection module is used to collect and query geological and mineral data for all mining areas within the region.

[0068] The 3D modeling module is used to create 3D models of all ore bodies within the query area.

[0069] The query area surface model building module is used to collect DEM data of the query area and build a surface model of the query area.

[0070] The module for obtaining the 3D model of the maximum protection range of the ore body is used to obtain the 3D model of the maximum protection range of the ore body based on the occurrence of the 3D model of the ore body.

[0071] The 3D model extension module is used to extend the 3D model of the maximum protection range of the ore body to the surface model.

[0072] The retaining zone range determination module is used to determine the retaining zone range of the queried construction project and place it on the same surface model.

[0073] The module for obtaining the status of mineral resources covered by the ore body is used to establish a query system. Based on the spatial relationship between the three-dimensional model of the maximum protection range of the ore body and the range of the retaining belt of the project being queried, the status of mineral resources covered by the ore body within the query area is obtained.

[0074] Combining all the above technical solutions, the advantages and positive effects of this invention are as follows: To improve the efficiency of mineral resource management and achieve rapid and high-quality development under the new circumstances, this invention innovatively proposes a method for querying mineral resource overburden. It eliminates the need to estimate the quantity of overburdened mineral resources or compile reports on overburdened mineral resources, and uses three-dimensional technology to quickly and intuitively determine the status of mineral resources overburdened by construction project land.

[0075] This invention innovatively integrates 3D orebody modeling, a 3D comprehensive display platform, and DEM (Digital Elevation Model) technology to query mineral resources covering areas. It integrates 3D models of all proven orebody reserves within the query area, constructing a unified and rapid query method. This invention collects geological and mineral-related data from all mining areas within the query area and performs reliability analysis on the collected data. Data from questionable mining areas is scientifically and systematically verified, and supplementary on-site investigations are conducted to obtain accurate results, ensuring the authenticity and reliability of the selected geological and mineral data and meeting the requirements for 3D orebody modeling. When multiple sets of data from different periods exist for the same mining area, the data that is reasonable, authentic, valid, comprehensive, and closest to the current reality is selected as the basis for 3D orebody modeling. Attached Figure Description

[0076] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this disclosure and, together with the description, serve to explain the principles of this disclosure;

[0077] Figure 1 This is a flowchart of a mineral resource overlay query method based on a three-dimensional ore body model provided by an embodiment of the present invention;

[0078] Figure 2 This is a schematic diagram of the uphill (downhill) movement angle provided in an embodiment of the present invention;

[0079] Figure 3 This is a three-dimensional spatial cross-sectional view of an iron ore mine provided in an embodiment of the present invention;

[0080] Figure 4 This is a three-dimensional model of an iron ore body provided in an embodiment of the present invention;

[0081] Figure 5 This is a three-dimensional model of the query area provided in this embodiment of the invention;

[0082] Figure 6 This is a single ore body strata movement angle diagram provided in an embodiment of the present invention;

[0083] Figure 7 This invention provides a method for checking the length of tangent segments around a ore body to determine whether the lowest point is the lowest elevation of the ore body and the highest point is above the surface.

[0084] Figure 8 This is a three-dimensional model diagram of the maximum protection range of a single ore body provided in an embodiment of the present invention;

[0085] Figure 9 This is a three-dimensional model diagram of the maximum protection range of the uncut ore body provided in the embodiments of the present invention;

[0086] Figure 10 This is a three-dimensional model diagram of the maximum protection range of the ore body provided in the embodiments of the present invention;

[0087] Figure 11 This is a three-dimensional model diagram of the ore body and its maximum protection range provided in the embodiments of the present invention;

[0088] Figure 12 This is a query diagram of the enclosure zone of a construction project provided by an embodiment of the present invention. Detailed Implementation

[0089] To make the above-mentioned objects, features, and advantages of the present invention more apparent and understandable, specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. Many specific details are set forth in the following description to provide a thorough understanding of the present invention. However, the present invention can be practiced in many other ways different from those described herein, and those skilled in the art can make similar modifications without departing from the spirit of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below.

[0090] The innovative aspects of the mineral resource overlay query method and system based on a 3D ore body model provided by this invention are as follows: Collecting geological and mineral data of all mining areas within the query area; performing 3D modeling of all ore bodies within the query area; collecting DEM data of the query area and establishing a surface model of the query area; obtaining a 3D model of the maximum protection range of the ore body based on the occurrence of the ore body 3D model; extending the 3D model of the maximum protection range of the ore body to the surface model; determining the scope of the retaining wall of the queried construction project and placing it on the same surface model (if the construction project is an underground project, constructing a 3D model of the project's retaining wall and importing it into the same surface model space); and obtaining the status of the mineral resources overlaid by the queried construction project within the query area based on the spatial relationship between the 3D model of the maximum protection range of the ore body and the scope of the retaining wall of the queried construction project.

[0091] Example 1, as Figure 1 As shown in the figure, the mineral resource overlay query method based on a three-dimensional ore body model provided by this invention includes:

[0092] S1 collects and queries geological and mineral data for all mining areas within the region;

[0093] S2, performs 3D modeling of all ore bodies in the query area;

[0094] S3, collect DEM data of the query area and build a surface model of the query area;

[0095] S4. Obtain the three-dimensional model of the maximum protection range of the ore body based on the occurrence of the three-dimensional model of the ore body.

[0096] S5 extends the three-dimensional model of the maximum protection range of the ore body to the surface model;

[0097] S6, determine the scope of the retaining wall of the queried construction project and place it on the same surface model. If the construction project is an underground project, construct a three-dimensional model of the retaining wall of the construction project and import it into the same surface model space.

[0098] S7. Establish a query system to obtain the status of mineral resources covered by the project within the query area based on the spatial relationship between the three-dimensional model of the maximum protection range of the ore body and the protection zone of the project being queried.

[0099] This invention also provides a mineral resource overlay query system based on a three-dimensional ore body model, comprising:

[0100] The data collection module is used to collect and query geological and mineral data for all mining areas within the region.

[0101] The 3D modeling module is used to create 3D models of all ore bodies within the query area.

[0102] The query area surface model building module is used to collect DEM data of the query area and build a surface model of the query area.

[0103] The module for obtaining the 3D model of the maximum protection range of the ore body is used to obtain the 3D model of the maximum protection range of the ore body based on the occurrence of the 3D model of the ore body.

[0104] The 3D model extension module is used to extend the 3D model of the maximum protection range of the ore body to the surface model.

[0105] The retaining zone range determination module is used to determine the retaining zone range of the queried construction project and place it on the same surface model.

[0106] The module for obtaining the status of mineral resources covered by the ore body is used to establish a query system. Based on the spatial relationship between the three-dimensional model of the maximum protection range of the ore body and the range of the retaining belt of the project being queried, the status of mineral resources covered by the ore body within the query area is obtained.

[0107] Example 2, as another embodiment of the present invention, provides a method for querying mineral resource overburden based on a three-dimensional ore body model, comprising:

[0108] Step 1, Data Preparation: Comprehensively collect and query all geological and mineral-related data and information from all mining areas within the region, and first conduct a reliability analysis on the collected data. Then, organize and analyze the data collected for each mining area to determine the basic data for 3D modeling of the ore body. Vectorize the selected data for each mining area using different mining 3D modeling software (such as 3D Mine, Micromine) to provide basic data for 3D modeling of the ore body.

[0109] Step 2, 3D modeling of the ore body:

[0110] All ore bodies within the query area will be modeled in 3D based on reliable data, and the appropriate 3D modeling method will be selected according to the data type.

[0111] Step 2.1, preparation work, including: first, preparing the tools required for 3D modeling, mainly including computers, MapGIS software, AutoCAD software, 3Dmine mining engineering software, ArcGIS geographic information system software or Micromine software, etc.

[0112] Secondly, prepare the basic data required for 3D modeling, mainly including the planar geological map, exploration line profile, reserve calculation map (table), and reserve calculation method in the exploration report.

[0113] Third, coordinate and elevation conversion should be carried out to unify the data into the 2000 National Geodetic Coordinate System and the 1985 National Elevation Datum.

[0114] Step 2.2, Modeling methods: There are four main methods for establishing orebody models: profile method, surface outcrop extrapolation method, top / bottom plate contour line method, and borehole data modeling method. The method to be used depends on the orebody delineation and reserve estimation method in the mining area exploration report. Different orebody models in the same mining area can adopt different modeling methods according to the report.

[0115] The process for establishing a 3D mining area resource reserve ore body model is as follows (using 3Dmine and Micromine software as examples):

[0116] a. Sectioning method:

[0117] 1) Sectional view processing: ① Import the original sectional view in JPG format into 3Dmine software and adjust the scale to 1:1000;

[0118] ② Use tools such as coordinate transformation and creating projected sections to transform the section view to the actual coordinate position.

[0119] 2) Orebody Boundary Drawing: Define a projection surface, and then draw the boundary outline of the orebody on the projection surface. The orebody drawn must be one that has participated in the resource reserve estimation; orebodies whose reserves have not been calculated are not drawn.

[0120] 3) Create the ore body model: ① Use the 3Dmine software entity → connect triangular mesh function to connect the ore bodies between the profiles one by one. When connecting, pay attention to the correspondence of the ore bodies and avoid connecting different ore bodies together.

[0121] ② Use the 3Dmine software entity → connect triangulation function to create the extrapolated portions at both ends of the ore body. The ore body inference should conform to the ore body inference principles and reserve estimation methods in the exploration report.

[0122] ③ Use the 3Dmine software's merging triangular mesh function to merge the ore body into a solid and perform solid verification. Once the verification is successful, the 3D modeling of the ore body is complete.

[0123] b. Outcropping method for surface outcrops:

[0124] 1) Determine the outcrop range of the ore body on the topographic and geological map, and extract the outline of the ore body;

[0125] 2) Use the 3Dmine software's Surface → Points and Lines onto DTM Surface function to project the ore body outline onto the surface model;

[0126] 3) Determine the thickness, length, and direction of extension of the ore body;

[0127] 4) Use the 3Dmine software's Entity → Connect Triangulation Network function to connect the ore body outcrops;

[0128] 5) Use the 3Dmine software entity → extension extrapolation function to extend the ore body;

[0129] 6) Use the 3Dmine software's merging triangular mesh function to merge the ore body into a solid and perform solid verification. Once the verification is successful, the 3D modeling of the ore body is complete.

[0130] c. Top / bottom plate contour line method:

[0131] 1) Extract contour lines and their corresponding elevation values ​​from the top / bottom plate reserve calculation map;

[0132] 2) The extracted contour lines are processed using 3DMine to assign elevations;

[0133] 3) Use surface modeling in 3DMine to create the top / bottom plate model;

[0134] 4) Determine the thickness of the ore body;

[0135] 5) Use the 3Dmine software entity → extension extrapolation function to extend the ore body;

[0136] 6) Use the 3Dmine software's merging triangular mesh function to merge the ore body into a solid and perform solid verification. Once the verification is successful, the 3D modeling of the ore body is complete.

[0137] d-Drilling data modeling method (Micromine software):

[0138] 1) Create data files in an Excel spreadsheet. The initial data files required include (drilling) wellhead files, lithology files, inclination survey files, sample analysis files, and trench (tunnel) sampling point files, etc.

[0139] 2) Import data from various Excel spreadsheets into Micromine to create a drilling (and other engineering) database;

[0140] 3) Display engineering data in Micromine's 3D space and cut exploration line profiles;

[0141] 4) On the exploration line profile, the ore bodies are connected in circles, and the ore bodies are connected from the exploration lines to the outer end to form a single ore body wireframe, i.e., a three-dimensional model of the ore body.

[0142] 5) Create wireframe libraries for multiple ore bodies (layers) in the mining area to form a collection of three-dimensional models of the ore bodies in the mining area.

[0143] Step 3, collect DEM data of the query area and build a surface model of the query area: The process of building a three-dimensional surface model using ArcGIS or 3Dmine software using the DEM method is as follows:

[0144] 1) Collect DEM digital elevation information data within the query area;

[0145] 2) Extract DEM elevation point cloud data within the query area;

[0146] 3) Import point clouds using 3DMine or ArcGIS software to create a surface model.

[0147] Step 4: Obtain the 3D model of the maximum protection range of the ore body: Load the 3D models of all ore bodies in a mining area into the 3D mining modeling software, construct the 3D model of the maximum protection range for each ore body, and finally use the Boolean calculation function of the software's entity module to combine the 3D models of the maximum protection range of all ore bodies to obtain the 3D model of the maximum protection range of the ore bodies in the mining area.

[0148] Step 4.1, Method for constructing a 3D model of the maximum protection range for a single ore body:

[0149] 1) In the 3D mining modeling software, open a single ore body model, enable editing mode, and keep other ore body models closed.

[0150] 2) Determine the overall strike and dip direction of the ore body. When the ore body has a complex and varied shape, it can be segmented to determine the occurrence data.

[0151] 3) By using geological data of the mining area, obtain characteristic parameters of ore body and surrounding rock lithology, bedding, joints, fracture coefficient and metamorphic degree, and determine the strata movement angle of the mining area.

[0152] 4) Around the ore body, the dip direction is determined by the angle of inclination of the ore body (the angle of inclination refers to the elevation angle of the area affected by the upward movement of the ore body along its dip direction, such as...). Figure 2Draw a line segment representing the angle and azimuth of the ore body's uphill (downhill) movement angle. Draw another line segment representing the angle and azimuth of the ore body's downhill movement angle in the opposite direction. Draw two line segments for each of the two azimuths corresponding to the strike movement angle (note that the same movement angle parameters are used for the same mining area). The length of the line segments for all four azimuths should be from the lowest point to the lowest elevation of the ore body and from the highest point to above the surface. Figure 2 In the figure, β is the angle of the ore body moving uphill, and γ is the angle of the ore body moving downhill.

[0153] 5) Copy the four line segments to the four different locations of the ore body, and paste them one by one to the tangent points of all the turning points and protrusions of the ore body in the corresponding locations, making sure not to miss any points.

[0154] 6) Check the length of the tangent segments around the ore body to see if the lowest point is the lowest elevation of the ore body and the highest point is above the ground surface. If not, it needs to be modified.

[0155] 7) In the 3D mining software, select the snap point to draw a polyline, and draw two closed lines at the lowest and highest points of all tangent segments around the ore body.

[0156] 8) Close the ore body model. Only the tangent segments around the ore body and the closed lines at the lowest and highest points are retained in the software interface.

[0157] 9) Generate triangular meshes within and between the two closed lines generated in 7), merge all triangular meshes, generate a solid model, and verify the solid model. This solid model is the three-dimensional model of the maximum protection range of the single ore body.

[0158] Step 4.2, construct a three-dimensional model of the maximum protection range of the ore body in the mining area: use the entity module-Boolean calculation in 3Dmine software to combine the three-dimensional models of the maximum protection range of all ore bodies to obtain the three-dimensional model of the maximum protection range of the ore body in the mining area without trimming.

[0159] Step 5: Extend the 3D model of the maximum protection range of the ore body to the surface model: Load the 3D models of the maximum protection range of all mining areas and the surface model of the query area into the same software interface. Use the surface clipping function to clip and delete the 3D model of the maximum protection range of the ore body that extends above the surface. Obtain the 3D model of the maximum protection range of the ore body below the surface.

[0160] Step 6: Determine the maximum boundary line of the retaining wall for the queried construction project: Determine the width of the retaining wall based on the nature and level of the queried construction project, delineate the scope of the retaining wall, and place the retaining wall scope onto the surface model of the query area. When the construction project is an underground construction project, construct an underground three-dimensional solid model of the project's retaining wall based on the project's lowest and highest elevations.

[0161] Retaining strip width: Refer to the retaining strip width for major buildings (structures) in the "Specifications for Coal Pillar Retention and Coal Mining under Pressure in Buildings, Water Bodies, Railways and Main Shafts". Specific situations should be analyzed during the investigation to select a reasonable (highest level of protection) protection grade.

[0162] Step 7, determine the status of mineral resources covered by the construction project land: Overlay the above-delineated surface or underground construction project retaining belt physical model with the three-dimensional model of the maximum protection range of the ore body in the mining area in three-dimensional space, and analyze and determine the status of ore body coverage in the mining area.

[0163] Based on the positional relationship between the physical model of the surface retaining zone or underground retaining zone of the construction project and the three-dimensional model of the maximum protection range of the ore body in the mining area in three-dimensional space, the overlay status of the ore body in the mining area is determined. The specific determination process is as follows:

[0164] Step 7.1, Overlapping a Single Mining Area: If the maximum protection range 3D model of a single mining area within the query area intersects entirely or partially with the scope of the project's retaining wall or the entity model of the retaining wall, then it is determined that the project overlaps a single mining area within the query area.

[0165] Step 7.2, Overlapping Multiple Mining Areas: If the maximum protection range 3D model of two or more mining areas within the query area intersects entirely or partially with the scope or entity model of the project's retaining wall, then the project is determined to have overlapped multiple mining areas within the query area.

[0166] Step 7.3, No Overburden: If the maximum protection range 3D model of the ore body in the query area does not intersect with the scope of the project's retaining belt or the entity model of the retaining belt, then it is determined that the project does not overburden the ore body in the query area.

[0167] In step S7.2 of this embodiment of the invention, if the three-dimensional models of the maximum protection range of two or more mining areas within the query area have all or part of an intersection with the scope or entity model of the retaining wall of the construction project, then it is determined that the construction project covers multiple mining areas within the query area, including:

[0168] Step 1: Decompose the non-intersecting areas and intersecting areas of the three-dimensional models of the maximum protection range of the ore bodies of two or more mining areas within the query area with the scope of the project's retaining belt or the entity model of the retaining belt, and obtain different regional components;

[0169] Step 2: Use image grayscale difference spectrum analysis to remove outliers from each intrinsic region component, and calculate the multi-size arrangement weights of each region component to construct a feature vector;

[0170] Step 3: Use intersection region analysis to reduce the order of feature vectors, and use a support vector machine recognizer to identify all or part of the three-dimensional models of the maximum protection range of ore bodies in two or more mining areas with the range of the project's retaining belt or a single sample of the retaining belt entity model.

[0171]

[0172] In the formula, Φ represents the feature of the intersection region, and v k Let covariance matrix be the variance matrix. Features that require downgrading in order to represent all or part of the 3D model of the maximum protection range of ore bodies in two or more mining areas compared to the scope of the retaining belt or the solid model of the retaining belt of the construction project. The feature mean of the training samples;

[0173] The value of the feature Φ of the intersection area is determined according to the principle of choosing the higher protection level. If it is a low protection level, it is a mining area covered by mining; if it is a high protection level, it is a mining area not covered by mining.

[0174] The covariance matrix v k The calculation formula is as follows:

[0175]

[0176] In the formula, X is the total number of three-dimensional models of the maximum protection range of ore bodies in two or more mining areas that are all or part of the project's retaining zone or retaining zone entity model, x is any single sample among the three-dimensional models of the maximum protection range of ore bodies in two or more mining areas that are all or part of the project's retaining zone or retaining zone entity model, and K is the exponent.

[0177] For two or more mining areas, the non-overlapping region decomposition is performed on all or part of the 3D model of the maximum protection range of the ore body and the scope of the project's retaining wall or the solid model of the retaining wall X(t). The expression is as follows:

[0178]

[0179] In the formula, y x (t) represents the x-th regional component, where X is the total number of three-dimensional models of the maximum protection range of ore bodies in two or more mining areas, including all or part of the model and the project's retaining wall range or retaining wall entity model, i.e., the total number of regions to be decomposed, and f(x) is the remainder; in the non-intersecting region decomposition method, y is obtained through decomposition. x The (t) component and the residual component f(t) can completely recover the original signal X(t).

[0180] The specific algorithms for decomposing disjoint regions include:

[0181] (1) Find all local maxima and minima on the original signal, and use cubic spline interpolation function to interpolate the corresponding upper and lower envelopes respectively;

[0182] (2) Calculate the mean of the two envelopes n1(t) = [p min (t)+p max [(tr)] / 2; where p min (t) represents the minimum interpolation of the upper and lower envelopes at time t, p max (t) represents the maximum interpolation of the upper and lower envelopes at time t;

[0183] (3) Calculate the difference between the original signal and the mean of the upper and lower envelopes, o1(t) = X(t) - n1(t); where X(t) represents the original regional component value of the three-dimensional model of the maximum protection range of the ore body of two or more mining areas at time t, which is the difference between the original regional component value of the project enclosure range or the enclosure entity model.

[0184] (4) Determine whether o1(t) is an intrinsic region component; if not, replace X(t) with o1(t) as the original signal and repeat the loop a times to obtain o. 1a =o 1(a-1) -n 1a This continues until the intrinsic region component condition is met, at which point it becomes the first intrinsic region component, denoted as y1 = o. 1a In the formula, o 1a This represents the difference between the mean values ​​of the upper and lower envelopes after repeating the cycle a times. 1(a-1) n represents the difference between the mean values ​​of the upper and lower envelopes after repeating the cycle a-1 times. 1a This represents the mean of the two envelopes after repeating the cycle a times.

[0185] (5) Separate y1 from the original signal to obtain the first-order residual signal f1(t) = X(t) - y1(t); where y1(t) represents the difference between the mean values ​​of the upper and lower envelopes that are repeated a times at time t;

[0186] (6) Treating f1(t) as a new original signal and repeating it, we obtain the l-th intrinsic region component y. a (t) and the a-order residual signal f a (t), until f a (t) is a monotonic function, and the decomposition ends.

[0187] The intersection region decomposition is performed on all or part of the 3D model of the maximum protection range of ore bodies in two or more mining areas and the scope of the retaining belt or the solid model of the retaining belt of the construction project, as expressed below:

[0188]

[0189] Among them, y x (t) represents the x-th regional component, where X is the total number of three-dimensional models of the maximum protection range of ore bodies in two or more mining areas, including all or part of the model and the project's retaining wall range or retaining wall entity model. Here, it can be the total number of decomposed regions, and g(t) is the residual term. In the intersection region decomposition method, the intrinsic region function y is defined. x (t) is an amplitude-frequency modulation (AM-FM) signal, expressed as:

[0190] y a (t)=B a (t)cos(τ a (t))

[0191] Among them, B a (t) represents the instantaneous amplitude of . For b a The instantaneous frequency of (t), and For the regional function b a (t), using Hilbert transform to obtain the analytic signals of each region, and finding H region functions that minimize the sum of the estimated enclosure ranges for each region; for b a Transform (t) to obtain its analytic signal and one-sided spectrum:

[0192]

[0193] Among them, b a (t) represents the intrinsic region function, and η(t) represents the impulse function;

[0194] By multiplying by an exponential function Adjust the center frequency of each intrinsic region function estimate to modulate the spectrum of each intrinsic region function to the corresponding fundamental frequency band:

[0195]

[0196] In the formula, c represents the amplitude of a certain region, and v i Indicates the amplitude frequency of a certain region;

[0197] Calculate the squared D² norm of the modulated signal gradient to estimate the range of the guard zone for each eigenfunction region:

[0198]

[0199] Where, Θ t [·] represents gradient operation. Denotes the squared D² norm; Represents an exponential function;

[0200] In the intersection region decomposition, to minimize the sum of bandwidths of each region, the following constrained variational model is established:

[0201]

[0202]

[0203] Among them, {b a}={b1…b X} represents a set of X regional components; {v a}={v1…v X} represents the set of center frequencies of each regional component. t represents time, and h represents the set of values ​​for the X regional components.

[0204] In step two, the method for calculating the multi-size arrangement weights of different regional components includes:

[0205] Inherent region components A one-dimensional size sequence of length M The sequence was obtained by coarsening:

[0206]

[0207] In the formula, χ is the scaling factor, Q = 1, 2, ... [M / χ], and [M / χ] represents the integer division of M / χ; when χ = 1, the coarse-grained sequence is the original sequence; M is the intrinsic region component. The length of is given by E, where E is a one-dimensional size sequence and i is a certain length of an intrinsic region component. For coarse-grained sequences, Qχ is the integer value of a certain scale factor, C is the length value of a certain scale factor after linear calculation, Y is the linear calculation coefficient of the scale factor, and Q is the integer value of M / χ.

[0208] For coarse-grained sequences After reconstructing the dimensions, we get:

[0209]

[0210] In the formula, For the sequence reconstructed from coarse-grained sequence size, This is the initial sequence for reconstructing the coarse-grained sequence size. For reconstructing the sequence size from the initial post-delay size using coarse-grained sequence size For the reconstructed coarse-grained sequence size, the sequence is delayed by n-1 embedding dimensions after initialization, where ε is the delay scale factor bias. The deviation value in the lag scale factor deviation is the deviation value after lag of m-1 dimensions, where m is the dimension of the size reconstruction; r is the r-th reconstruction component, r = 1, 2…M-(n-1), where n is the embedding dimension, and Γ is the lag size; r1, r2…r n Represents reconstructed components The index of the column containing each element will Arranged in ascending order as follows:

[0211]

[0212] In the formula, To be Reconstructed components arranged in ascending order The first delayed index sequence of each element in the column. To be Reconstructed components arranged in ascending order The delayed index sequence of each element in the column. To be Reconstructed components arranged in ascending order The nth delayed index sequence of each element in the column;

[0213] If equal values ​​exist in the reconstructed components, they are arranged in chronological order; for any coarse-grained sequence We obtain a sequence of symbols B(l) = [c1, c2, ..., c1] n ], where l = 1, 2, ..., N, and N ≤ n! ; the reconstructed sequence with an embedding dimension of n has n! permutations, and the symbol sequence B(l) is one of these permutations. Calculate the probability O of each symbol sequence. l (l = 1, 2, ..., N); B(l) represents the value for any coarse-grained sequence. A sequence of symbols, c, is obtained. n Let l be the sequence of symbols, N be the number of symbols, n be the embedding dimension, and O be the number of symbols. l The probability of each symbol sequence occurring;

[0214] Define the permutation entropy R of different symbol sequences in the form of information entropy. C (n) is:

[0215]

[0216] When O r = 1 / n!, R C (n) Take the maximum value ln(n!), and let R C (n) After normalization, we get:

[0217] R C =R C (n) / ln(n!)

[0218] Among them, R C The normalized permutation weights, 0 <R C <1,R C The value can reflect and amplify minute changes in a size sequence;

[0219] The feature vector is constructed by calculating the multi-size permutation weights of each intrinsic region component:

[0220]

[0221] In the formula, ∧ represents the feature vector value constructed by the multi-size arrangement weights of each intrinsic region component. The permutation weights after multi-size normalization of the first intrinsic region component. The permutation weights of the Xth intrinsic region component after multi-size normalization are given. Reconstruct the permutation weights for the multi-size normalization process of the first intrinsic region component. The permutation weights are reconstructed using multi-size normalization processing for the Xth intrinsic region component.

[0222] Example 2, taking an iron mine as an example, provides a detailed introduction to the workflow and technical principles of a mineral resource overburden query system based on 3D ore body modeling.

[0223] The first step is preparation.

[0224] 1) Comprehensively collect geological and mineral-related data and information for a certain iron ore area, and first conduct a reliability analysis on the collected data. Then, organize and analyze the collected mining area data to determine the basic data for 3D modeling of the ore body.

[0225] The iron ore mine is located in a low mountain and hilly area with stable natural slopes of 20°–40°. The Quaternary strata are 13–26 m thick and structurally well-developed. The surrounding rock dips at 290–310° NW with an angle of 50–70°, while the ore body dips at 290°–300° NW with an angle of 60–70°.

[0226] This iron ore mine, along with other nearby iron ore mines, belongs to the Anshan-type sedimentary metamorphic magnetite-quartzite deposit. The ore type, ore structure, texture, mineral composition, chemical composition, magnetite grain size, and occurrence state are basically similar. The final slope angles of the open-pit mines in nearby iron ore mines are: south slope 39.7°, north slope 43.1°, east slope 41.6°, and west slope 44.2°. Considering the actual conditions of this iron ore mine and comparing it with the final slope angle data of nearby iron ore mines, the strata movement angle of this iron ore mine area is determined to be 45°.

[0227] 2) Prepare the tools required for 3D modeling, mainly including computers, MapGIS software, AutoCAD software, 3Dmine mining engineering software, Micromine software, and ArcGIS geographic information system software.

[0228] 3) Prepare the basic data required for 3D modeling, mainly including the geological plan map, exploration line profile, reserve calculation map (table), and reserve calculation method from the exploration report. Vectorize the selected modeling data and maps to provide basic data for 3D modeling of the ore body or verification of the 3D model of the ore body.

[0229] 4) Perform coordinate and elevation conversion, and convert the data to the 2000 National Geodetic Coordinate System and the 1985 National Elevation Datum.

[0230] The second step involves modeling methods. There are four main methods for establishing orebody models: the profile method, the surface outcrop extrapolation method, the top / bottom contour line method, and the borehole data modeling method. The method chosen depends on the orebody delineation and reserve estimation methods in the mining area exploration report. Different orebodies within the same mining area can employ different modeling methods based on the report. The profile method was used for the 3D modeling of the orebody in a certain iron ore area.

[0231] The modeling methods and processes are as follows (using 3Dmine software as an example):

[0232] 1) Sectional view processing: ① Import the original sectional view in JPG format into 3Dmine software and adjust the scale to 1:1000;

[0233] ② Use tools such as coordinate transformation and creating projected sections to transform the section view to its actual coordinate position (see...). Figure 3 (Three-dimensional spatial cross-section of an iron mine).

[0234] 2) Drawing the boundary of the ore body:

[0235] Define a projection plane, and then draw the boundary outline of the ore body on the projection plane. The ore body drawn must be one that has participated in the resource reserve estimation; ore bodies whose reserves have not been calculated should not be drawn.

[0236] 3) Create an ore body model:

[0237] ① Use the 3Dmine software's Entity → Connect Triangular Mesh function to connect the ore bodies between profiles one by one. When connecting, pay attention to the correspondence of the ore bodies and avoid connecting different ore bodies together.

[0238] ② Use the 3Dmine software entity → connect triangulation function to create the extrapolated portions at both ends of the ore body. The ore body inference should conform to the ore body inference principles and reserve estimation methods in the exploration report.

[0239] ③ Use the merging triangulation function of 3Dmine software to merge the ore body into a solid and perform solid verification. Once the verification is successful, the 3D modeling of the ore body is complete (see...). Figure 4 ).

[0240] The third step is to collect DEM data for the query area and build a surface model of the query area:

[0241] The process of creating a 3D surface model using ArcGIS or 3Dmine software via the DEM method is as follows:

[0242] 1) Collect DEM digital elevation information data within the query range;

[0243] 2) Extract DEM elevation point cloud data within the query area;

[0244] 3) Import point clouds using 3DMine or ArcGIS software to create a surface model (see...). Figure 5 ).

[0245] Step 4: Obtain a 3D model of the maximum protection range of the ore body:

[0246] The three-dimensional models of all ore bodies in a certain iron ore area are loaded into a three-dimensional mining modeling software. A three-dimensional model of the maximum protection range is constructed for each ore body. The software’s [Entity Module] - Boolean calculation function is used to combine the three-dimensional models of the maximum protection range of all ore bodies to obtain the three-dimensional model of the maximum protection range of the ore bodies in the certain iron ore area.

[0247] 1) Method for constructing a 3D model of the maximum protection range of a single ore body:

[0248] ① In the 3D mining modeling software, open a single ore body model, enable editing mode, and keep other ore body models closed.

[0249] ② Determine the overall strike and dip direction of the ore body. When the ore body has a complex and varied shape, it can be segmented to determine the occurrence data.

[0250] ③ By obtaining the geological data of the mining area, the characteristic parameters of lithology, bedding, joints, fracture coefficient and metamorphism of the ore body and surrounding rocks are obtained, and the strata movement angle of the mining area (note that the same movement angle parameters are used in the same mining area, including the angle of the downhill angle of the ore body, the angle of the uphill angle of the ore body, and the angle of the strike movement angle) is determined to be 45°.

[0251] ④ Around the ore body, draw a line segment along the dip direction using the angle and direction of the ore body's downhill corner, and another line segment along the opposite dip direction using the angle and direction of the ore body's uphill corner. At each end of the strike direction, draw two line segments along the angle and direction of the corresponding strike movement angle. The length of the line segments in all four directions should extend from the lowest point to the lowest elevation of the ore body and from the highest point to above the surface.

[0252] ⑤ Copy the four line segments to four different locations on the ore body, and paste them one by one to the tangent points of all the turning points and protrusions of the ore body at the corresponding locations, making sure not to miss any points (see...). Figure 6 ).

[0253] ⑥ Check the length of the tangent segments around the ore body to see if the lowest point is the lowest elevation of the ore body and the highest point is above the ground surface. If not, it needs to be modified (see...). Figure 7 ).

[0254] ⑦ In the 3DMine 3D mining software, select the snap point to draw a polyline, and draw two closed lines at the lowest and highest points of all tangent segments around the ore body.

[0255] ⑧ Close the ore body model. Only retain the tangent segments around the ore body and the closed lines at the lowest and highest points in the software interface.

[0256] ⑨ Generate triangular meshes within and between the two closed lines generated in step ⑦, merge all triangular meshes, generate a solid model, and verify the solid model. This solid model is the three-dimensional model of the maximum protection range of this single ore body (see...). Figure 8 ).

[0257] 2) Construct a 3D model of the maximum protection range of an iron ore body: Using the Boolean calculation function in the [Solid Module] of 3Dmine software, combine all the 3D models of the maximum protection range of the ore bodies to obtain the 3D model of the maximum protection range of the ore bodies of the iron ore body without clipping (see...). Figure 9 ).

[0258] The fifth step is to extend the three-dimensional model of the maximum protection range of the ore body onto the surface model.

[0259] Load the 3D model of the maximum protection range of the ore body in the mining area and the surface model of the query area into the same interface of the 3Dmine software. Use the "Surface Clipping" function to clip and delete the 3D model of the maximum protection range of the ore body that extends above the surface. Obtain the 3D model of the maximum protection range of the ore body below the surface (see...). Figure 10 , Figure 11 ).

[0260] Step 6: Determine the maximum enclosure zone of the construction project being queried.

[0261] The width of the retaining wall is determined based on the nature and level of the queried construction project, thus defining the scope of the retaining wall. The scope of the retaining wall is then mapped onto the surface model of the queried area using coordinate data (see [link]). Figure 12 (In the diagram, the double lines in the diagonal area represent the queried construction project, the area outside the double lines in the diagonal area represents the scope of the construction project's retaining wall, and the polygonal area represents the maximum surface protection area of ​​a certain iron ore body.) When the construction project is an underground construction project, a three-dimensional underground solid model of the construction project's retaining wall is constructed based on the project's lowest and highest elevations.

[0262] Retaining strip width: Refer to the retaining strip width for major buildings (structures) in the "Specifications for Coal Pillar Retention and Coal Mining under Pressure in Buildings, Water Bodies, Railways and Main Shafts". Specific situations should be analyzed during the investigation to select a reasonable (highest level of protection) protection grade.

[0263] Step 7: Determine the status of mineral resources covered by the construction project. Overlay the physical model of the surface or underground retaining zone of the construction project with the three-dimensional model of the maximum protection range of the ore body in a certain mining area in three-dimensional space, and analyze and determine the status of ore body coverage in the mining area.

[0264] Based on the positional relationship in three-dimensional space between the physical model of the surface or underground retaining zone of the construction project and the three-dimensional model of the maximum protection range of the ore body in the mining area, the overlay status of the ore body in the query area by the construction project is determined. The specific determination process is as follows:

[0265] a. Overlapping a single mining area:

[0266] If the maximum protection range 3D model of a single mining area within the query area intersects entirely or partially with the scope of the project's retaining wall or the entity model of the retaining wall, then the project is determined to cover a single mining area within the query area.

[0267] b. Overlapping multiple mining areas:

[0268] If the maximum protection range 3D model of two or more mining areas within the query area intersects entirely or partially with the scope or physical model of the project's retaining wall, then the project is determined to cover multiple mining areas within the query area.

[0269] c. No overlay is present:

[0270] If the maximum protection range 3D model of the ore body in the query area does not intersect with the scope of the project's retaining belt or the entity model of the retaining belt, then it is determined that the project does not cover any ore body in the query area.

[0271] Based on the above query process, and considering the intersection between the spatial relationship between the scope of the project's retaining wall and the three-dimensional model of the maximum protection range of the ore body in a certain iron ore area, it is determined that the project in question covers a certain iron ore within the query area.

[0272] In the above embodiments, the descriptions of each embodiment have different focuses. For parts that are not described in detail or recorded in a certain embodiment, please refer to the relevant descriptions of other embodiments.

[0273] The information interaction and execution process between the above-mentioned devices / units are based on the same concept as the method embodiments of the present invention. For details on their specific functions and technical effects, please refer to the method embodiments section, and they will not be repeated here.

[0274] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the above-described division of functional units and modules is merely an example. In practical applications, the above functions can be assigned to different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiments can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit. Furthermore, the specific names of the functional units and modules are only for easy differentiation and are not intended to limit the scope of protection of this invention. The specific working process of the units and modules in the above system can be referred to the corresponding process in the foregoing method embodiments.

[0275] This invention also provides a computer device comprising: at least one processor, a memory, and a computer program stored in the memory and executable on the at least one processor, wherein the processor executes the computer program to implement the steps in any of the above method embodiments.

[0276] This invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, can implement the steps described in the various method embodiments above.

[0277] This invention also provides an information data processing terminal, which, when executed on an electronic device, provides a user input interface to implement the steps described in the above method embodiments. The information data processing terminal is not limited to mobile phones, computers, or switches.

[0278] This invention also provides a server that, when executed on an electronic device, provides a user input interface to implement the steps described in the above method embodiments.

[0279] This invention provides a computer program product that, when run on an electronic device, enables the electronic device to implement the steps described in the various method embodiments above.

[0280] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the methods of the above embodiments of this application can be implemented by a computer program instructing related hardware. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms. The computer-readable medium can include at least: any entity or device capable of carrying the computer program code to a photographing device / terminal device, a recording medium, a computer memory, a read-only memory (ROM), a random access memory (RAM), an electrical carrier signal, a telecommunication signal, and a software distribution medium. Examples include USB flash drives, portable hard drives, magnetic disks, or optical disks.

[0281] The above description is merely a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any modifications, equivalent substitutions, and improvements made by those skilled in the art within the scope of the technology disclosed in the present invention, and within the spirit and principles of the present invention, should be covered within the scope of protection of the present invention.

Claims

1. A method for querying mineral resource overlay based on a three-dimensional ore body model, characterized in that, The method includes the following steps: S1 collects and queries geological and mineral data for all mining areas within the region; S2, performs 3D modeling of all ore bodies in the query area; S3, collect DEM data of the query area and build a surface model of the query area; S4. Obtain the three-dimensional model of the maximum protection range of the ore body based on the occurrence of the three-dimensional model of the ore body. S5 extends the three-dimensional model of the maximum protection range of the ore body to the surface model; S6, determine the scope of the retaining wall of the queried construction project and place it on the same surface model. If the construction project is an underground project, construct a three-dimensional model of the retaining wall of the construction project and import it into the same surface model space. S7. Establish a query system to obtain the status of mineral resources covered by the project within the query area based on the spatial relationship between the three-dimensional model of the maximum protection range of the ore body and the protection zone of the project being queried. In step S4, the maximum protection range three-dimensional model of the ore body is obtained based on the occurrence of the ore body three-dimensional model, including: loading all the ore body three-dimensional models of a mining area into the three-dimensional mining modeling software, constructing the maximum protection range three-dimensional model for each ore body, and using the software entity module-Boolean calculation to combine all the maximum protection range three-dimensional models of the ore body to obtain the maximum protection range three-dimensional model of the ore body in the mining area. Using the software entity module - Boolean calculation, the three-dimensional models of the maximum protection range of all ore bodies are combined to obtain the three-dimensional model of the maximum protection range of the ore bodies in this mining area, including: 1) Enable editing mode for a single ore body model in the 3D mining modeling software; 2) Determine the overall strike and dip direction of the ore body; 3) By using geological data of the mining area, obtain characteristic parameters of lithology, bedding, joints, fracture coefficient, and metamorphic degree of the ore body and surrounding rocks, and determine the strata movement angle of the mining area; 4) Around the ore body, draw line segments in four directions according to the dip direction and the angle of the ore body's uphill movement. The length of the line segments should be from the lowest point to the lowest elevation of the ore body and from the highest point to above the ground surface. 5) Copy the four line segments to the four different locations of the ore body, and paste them one by one to the tangent points of all the turning points and protrusions of the ore body in the corresponding locations; 6) Check the length of the tangent segments around the ore body to see if the lowest point is the lowest elevation of the ore body and the highest point is above the ground surface; 7) In the 3D mining software, select the snap point to draw a polyline, and draw two closed lines at the lowest and highest points of all tangent segments around the ore body; 8) Close the ore body model; only the tangent segments around the ore body and the closed lines at the lowest and highest points are displayed in the software interface. 9) Generate triangular meshes within and between the two closed lines, merge all triangular meshes, generate a solid model, and verify the solid model. This solid model is the three-dimensional model of the maximum protection range of the single ore body. In step S5, the three-dimensional model of the maximum protection range of the ore body is extended to the surface model surface, including: loading the three-dimensional model of the maximum protection range of the ore body of all mining areas and the surface model of the query area into the same software interface, using the surface clipping function to clip and delete the three-dimensional model of the maximum protection range of the ore body that is higher than the surface range, and obtaining the three-dimensional model of the maximum protection range of the ore body below the surface. In step S7, the status of mineral resources covered by the queried construction project within the query area is obtained, including: overlaying the physical model of the delineated surface or underground retaining zone of the construction project with the three-dimensional model of the maximum protection range of the ore body in the mining area in three-dimensional space, and analyzing and determining the status of ore body coverage in the mining area; the specific determination steps are as follows: S7.1, Covering a single mining area: If the maximum protection range 3D model of a single mining area within the query area intersects with the scope of the project's retaining wall or the entity model of the retaining wall, then the project is determined to cover a single mining area within the query area. S7.2, Covering Multiple Mining Areas: If the maximum protection range 3D model of two or more mining areas within the query area intersects with the scope or entity model of the project's retaining wall, then the project is determined to cover multiple mining areas within the query area. S7.3, No Overburden: If the maximum protection range 3D model of the ore body in the query area does not intersect with the scope of the project's retaining belt or the entity model of the retaining belt, then the project is determined to not overburden the ore body in the query area.

2. The mineral resource overlay query method based on a three-dimensional ore body model according to claim 1, characterized in that, In step S2, three-dimensional modeling is performed on all ore bodies in the query area, including: three-dimensional modeling of different ore bodies in the same mining area using one of the following methods: profile method, surface outcrop extrapolation method, top / bottom plate contour line method, and borehole data modeling method.

3. The mineral resource overlay query method based on a three-dimensional ore body model according to claim 1, characterized in that, In step S3, a surface model of the query area is established using ArcGIS or 3Dmine software via the DEM method. This specifically includes: S3.1, Collect DEM digital elevation information data within the query area; S3.2, Extract DEM elevation point cloud data within the query area; S3.3, use 3DMine or ArcGIS software to import point clouds and build a surface model.

4. The mineral resource overlay query method based on a three-dimensional ore body model according to claim 1, characterized in that, Step S7.2 specifically includes: Step 1: Decompose the non-intersecting areas and intersecting areas of the three-dimensional models of the maximum protection range of the ore bodies of two or more mining areas within the query area with the scope of the project's retaining belt or the entity model of the retaining belt, and obtain different regional components; Step 2: Use image grayscale difference spectrum analysis to remove outliers from each intrinsic region component, and calculate the multi-size arrangement weights of each region component to construct a feature vector; Step 3: Use intersection region analysis to reduce the order of feature vectors, and use a support vector machine recognizer to identify all or part of the three-dimensional models of the maximum protection range of ore bodies in two or more mining areas with the range of the project's retaining belt or a single sample of the retaining belt entity model. New feature vectors are formed by linear projection of the feature vectors, and the intersection region of the features is calculated. The expression is: ； In the formula, Features of the intersection region Let covariance matrix be the variance matrix. Features that require downgrading in order to represent all or part of the 3D model of the maximum protection range of ore bodies in two or more mining areas compared to the scope of the retaining belt or the solid model of the retaining belt of the construction project. The feature mean of the training samples Intersection region characteristics The value is determined according to the highest protection level. If it is a low protection level, it is a mining area covered by mining; if it is a high protection level, it is a mining area not covered by mining. Covariance matrix The calculation formula is as follows: ； In the formula, The total number of 3D models representing all or part of the maximum protection area of ​​ore bodies in two or more mining areas, and the total number of project enclosure belt ranges or enclosure belt entity models. The maximum protection range of two or more mining areas is represented by all or part of the three-dimensional model and any single sample from the project's retaining wall range or retaining wall entity model.

5. The mineral resource overlay query method based on a three-dimensional ore body model according to claim 4, characterized in that, For two or more mining areas, the maximum protection range of the 3D model of the ore body is wholly or partially related to the scope of the project's retaining wall or the solid model of the retaining wall. The expression for decomposing the non-overlapping regions is: ； In the formula, For the first Each region component, The total number of 3D models representing the maximum protection area of ​​ore bodies in two or more mining areas, including all or part of the model, and the total number of project enclosure belts or enclosure belt entity models, i.e., the total number of areas to be decomposed. The remainder; in the non-overlapping region decomposition method, the terms obtained through decomposition are... Components and residual components Able to convert the original signal Fully recovered.

6. The mineral resource overlay query method based on a three-dimensional ore body model according to claim 4, characterized in that, In step two, the multi-size arrangement weights of each regional component are calculated to construct a feature vector, including: Inherent region components Length is One-dimensional size sequence The coarse-graining process yields the following sequence: ； In the formula, As a scale factor, , Indicates to Rounding; when At that time, the coarse-grained sequence is the original sequence; For inherent region components, For intrinsic region components Length, It is a one-dimensional size sequence. For a certain length of the intrinsic region component, It is a coarse-grained sequence. The value after rounding down to a certain scale factor. The length value after linear calculation for a certain scale factor. The coefficients are calculated linearly for the scaling factor. To The integer value; For coarse-grained sequences After reconstructing the dimensions, we get: ； In the formula, For the sequence reconstructed from coarse-grained sequence size, This is the initial sequence for reconstructing the coarse-grained sequence size. For reconstructing a sequence with a delayed size after initialization from a coarse-grained sequence size, For the initial post-delay of coarse-grained sequence size reconstruction A sequence of embedding dimensions, For the delay scale factor bias, Delay in the lag scaling factor bias The deviation value after dimension, The dimension for size reconstruction; For the first One reconstructed component, , For the embedding dimension, For delay dimensions; using Represents reconstructed components The index of the column containing each element will Arranged in ascending order as follows: ； In the formula, To be Reconstructed components arranged in ascending order The first delayed index sequence of each element in the column. To be Reconstructed components arranged in ascending order The nth delayed index sequence of each element in the column; If equal values ​​exist in the reconstructed components, they are arranged in chronological order; for any coarse-grained sequence A set of symbol sequences is obtained ,in, ,and Embedding dimension is The size reconstruction sequence has a total of Permutations, symbol sequences One of the permutations involves calculating the probability of each symbol sequence occurring. ; For any coarse-grained sequence A sequence of symbols is obtained. For symbol sequences, For signed numbers, For a sign positive integer, For the embedding dimension, The probability of each symbol sequence occurring; Define the permutation entropy of different symbol sequences in the form of information entropy. for: ； when , Take the maximum value ,Will After normalization, we get: ； in, These are the normalized permutation weights. , The value can reflect and amplify minute changes in a size sequence; The feature vector is constructed by calculating the multi-size permutation weights of each intrinsic region component: ； In the formula, Feature vector values ​​are constructed by arranging the weights of each intrinsic region component using multiple dimensions. The permutation weights after multi-size normalization of the first intrinsic region component. The permutation weights of the Xth intrinsic region component after multi-size normalization are given. Reconstruct the permutation weights for the multi-size normalization process of the first intrinsic region component. The permutation weights are reconstructed using multi-size normalization processing for the Xth intrinsic region component.

7. A mineral resource overlay query system based on a three-dimensional ore body model, characterized in that, This system implements the mineral resource overburden query method based on a three-dimensional ore body model as described in any one of claims 1-6. The system includes: The data collection module is used to collect and query geological and mineral data for all mining areas within the region. The 3D modeling module is used to create 3D models of all ore bodies within the query area. The query area surface model building module is used to collect DEM data of the query area and build a surface model of the query area. The module for obtaining the 3D model of the maximum protection range of the ore body is used to obtain the 3D model of the maximum protection range of the ore body based on the occurrence of the 3D model of the ore body. The 3D model extension module is used to extend the 3D model of the maximum protection range of the ore body to the surface model. The retaining zone range determination module is used to determine the retaining zone range of the queried construction project and place it on the same surface model. The module for obtaining the status of mineral resources covered by the ore body is used to establish a query system. Based on the spatial relationship between the three-dimensional model of the maximum protection range of the ore body and the range of the retaining belt of the project being queried, the status of mineral resources covered by the ore body within the query area is obtained.

Images