A rapid implicit geological modeling method for ore bodies of weathering crust elution deposits

Abstract

The present invention discloses a method for rapid implicit geological modeling of ore bodies in weathering crust eluvial deposits, including steps such as generating a triangulated mesh of the ore body's roof and floor surfaces within the modeling range, importing a surface DTM model within the modeling range, constructing a surface DTM surface node search point set for the ore body's roof and floor DTM nodes, and performing estimated elevation adjustments on the ore body's roof and floor DTM surface nodes. This method omits a series of complex manual operations, such as ore body interpretation on profiles and connecting ore bodies between profiles, in traditional geological modeling methods. The geological modeling process is completed automatically, eliminating the need for manual drawing operations. Furthermore, the method embodies the mineralization principles of eluvial deposits, features a simple and clear algorithm, and offers rapid calculations. Ore body reconstruction is simple and rapid after exploration data updates and changes to ore body interpretation standards.

Description

Technical Field

[0001] The invention belongs to the technical field of three-dimensional digital geological modeling and relates to a method for rapid implicit geological modeling of an ore body of a weathering crust elution type deposit. Background Art

[0002] Three-dimensional geological models are the foundation of modern digital mining. Three-dimensional geological modeling utilizes computer graphics technology to generate a three-dimensional model of a geological object by integrating various geological information and geological interpretation results. Implicit modeling is a parametric spatial interpolation technique. One of its key features is the automated geological modeling achieved through the use of spatial interpolation algorithms. The reliability and practicality of implicit modeling have been proven in geological modeling of various types of mineral deposits, making it a key development and application direction in the field of geological modeling, both now and in the future.

[0003] Weathering crust eluvial deposits are a type of mineral resource with very distinct characteristics. Influenced by weathering and eluvial mineralization, this type of deposit generally has the following occurrence characteristics: 1. The ore body is located in the surface weathering crust. The geological profile can generally be divided into the topsoil layer, the fully weathered layer, the semi-weathered layer, and the bedrock from top to bottom; 2. The ore body is often layered and tilts from the ridge to the foot of the mountain. The ore body occurrence elevation is generally higher than the local lowest erosion base level; 3. The occurrence of the ore body is jointly influenced and controlled by the strata and surface morphology. Generally speaking, weathering crust eluvial deposits have layered ore bodies and simple morphology, which are controlled by the weathering crust strata and surface topography. Current implicit geological modeling techniques primarily focus on exploration data obtained through drilling and other exploration projects. While developing complex algorithms and computational processes, they completely ignore the combined influence and control of surface topography and strata on ore body morphology during the formation of weathering crust-eluvial deposits. This results in significant discrepancies between implicit geological modeling results and the actual occurrence of the deposits. Consequently, current software supporting implicit geological modeling is generally insufficiently practical for modeling weathering crust-eluvial deposits. Summary of the Invention

[0004] The main purpose of the present invention is to provide a rapid implicit geological modeling method for the ore body of a weathering crust elution type deposit. This method automatically and quickly calculates the elevations of all nodes on the top and bottom surfaces of the ore body within the modeling range based on the borehole grade information revealed by geological exploration drilling and combined with the industrial grade. At the same time, it takes into account the influence of the surface topography on the weathering elution mineralization and corrects the elevations of each node on the top and bottom of the ore body. It fully considers the mineralization principle and occurrence characteristics of the weathering crust elution type deposit, greatly improving the accuracy and reliability of engineering geological three-dimensional modeling.

[0005] To this end, the present invention provides a method for rapid implicit geological modeling of ore bodies of weathering crust elution-type deposits, comprising:

[0006] Generate the DTM triangulated network of the roof and floor of the ore body according to the modeling scope and geological modeling accuracy requirements; import the surface DTM model within the modeling scope;

[0007] Taking each node k on the DTM surface of the ore body roof and floor as the center, increase the search radius from small to large, and search the node closest to the node k on the surface DTM surface in the horizontal direction to form the search point set of node k;

[0008] Using the search point set of node k and the inverse distance power method, the elevation of each point in the search point set is regarded as a known value. The elevation of node k is adjusted for the first time, and the elevation of each node on the ore body roof and floor DTM is adjusted to fit with the surface DTM surface. If the fitting accuracy does not meet the design requirements, the triangulation density is adjusted and the calculation is repeated again after returning to the step of generating the triangulation of the ore body roof and floor DTM surface.

[0009] Using drilling exploration data and industrial grade values, determine the top and bottom depth points of the ore body in each drill hole in three-dimensional space;

[0010] Using the top and bottom depth points of the ore body on the drill hole, calculate the sinking depth values ​​of the top and bottom of the ore body relative to the drill hole opening;

[0011] Taking each node k on the DTM surface of the ore body roof and floor as the center, increase the search radius from small to large, and search for the drill hole closest to the node in the horizontal direction. Stop searching when the number of drill holes meets the design requirements, and use the searched drill holes to construct the search drill hole set of the node;

[0012] Using the search drill hole set, a similar method to searching for surface DTM nodes is used to perform a second estimation adjustment on the elevation of each node on the ore body roof and floor DTM after the first estimation adjustment. The elevation of each node on the ore body roof and floor DTM is adjusted to match the buried depth of the ore body roof and floor. If the matching accuracy does not meet the engineering requirements, the drill hole search parameters are adjusted and the calculation is returned to the search drill hole set construction step for recalculation.

[0013] Using the mining area boundary, the DTM surface of the ore body top and bottom plates after the second valuation adjustment is cut to obtain the final DTM model of the ore body top and bottom plates, completing the rapid implicit geological modeling of the ore body model in the mining area.

[0014] Specifically, for the node kr on the DTM surface of the ore body roof, based on the obtained search drill hole set, the starting depth of the ore body on each drill hole is taken as a known value, and the inverse power method of distance is used to estimate the elevation drop value of the node kr, and the node kr elevation is adjusted accordingly;

[0015] For the node kf on the DTM surface of the ore body bottom plate, according to the obtained search drill hole set, the starting depth of the ore body on each drill hole + the ore body thickness is taken as a known value, and the inverse power method of distance is used to estimate the elevation drop value of the node kf, and the node kf elevation is adjusted accordingly.

[0016] Specifically, for the node kr on the DTM surface of the ore body roof, if there is only one drill hole in the search drill hole set, the ore body starting depth value of the drill hole is directly used to perform the subsidence adjustment calculation for the elevation of the current node;

[0017] For the node kf on the DTM surface of the ore body bottom plate, if there is only one drill hole in the search drill hole set, the ore body starting depth + ore body thickness value of the drill hole is directly used to perform the subsidence adjustment calculation for the elevation of the current node.

[0018] Specifically, according to the scope of the mining area and the mining conditions, a rectangular geological modeling range that completely covers the mining area is determined. The rectangular geological modeling range is east-west and north-south.

[0019] Specifically, according to the geological modeling accuracy requirements, an east-west and north-south triangulation value is selected for the DTM surface of the ore body top and bottom plates, and a triangulation network of the DTM surface of the ore body top and bottom plates within the modeling range is generated.

[0020] Specifically, if there is only one node in the search point set, the elevation of the point in the search point set is directly used to assign elevation values ​​to the corresponding nodes of the ore body roof and floor DTM.

[0021] Specifically, the surface DTM model of the mining area is constructed using surface data such as contour lines and elevation points using three-dimensional mining software.

[0022] Specifically, the distance power inverse method is used to calculate the elevation value of each node k on the DTM surface of the ore body roof and floor for the first time. The calculation formula is as follows:

[0023]

[0024] Where: Z is the elevation value of point k adjusted for the first time

[0025] z1, z2, z3…z m The elevation value of the point in the search point set

[0026] d1, d2, d3…d m is the horizontal distance between the point in the search point set and point k

[0027] Q is the power, which is 2;

[0028] If d1, d2, d3…d m Existence d h =0, then Z is directly assigned to z h; If there is only one point in the search point set and its elevation is z, then z is also directly assigned to Z.

[0029] Specifically, the distance power inverse method is used to calculate the k of each node on the DTM surface of the ore body roof. r and each node k on the DTM surface of the ore body bottom f The calculation formulas for elevation correction values ​​are as follows:

[0030]

[0031] Where: R Δz is the node k on the DTM surface of the ore body roof r The elevation reduction value of the second elevation adjustment, i.e., the subsidence value;

[0032] rd1, rd2, rd3…rd n The depth of the ore body roof of the drill hole in the search drill hole set, that is, the distance between the ore body roof and the hole mouth in the drill hole;

[0033] F Δz is the node k on the DTM surface of the ore body bottom f The elevation reduction value of the second elevation adjustment, i.e., the subsidence value;

[0034] fd1, fd2, fd3…fd n The depth of the ore body floor of the drill hole in the search drill hole set, that is, the distance between the ore body floor and the hole mouth in the drill hole;

[0035] d1, d2, d3…d n is the horizontal distance between the drilling hole opening in the search drilling set and point k

[0036] p is the power, and its value is 3;

[0037] If d1, d2, d3…d n Existence d t =0, node k r Elevation directly reduces rd t , and node k f Elevation directly reduces fd t If there is only one drill hole in the search drill hole set, and the top plate depth of the ore body in the drill hole is rd and the bottom plate depth is fd, then the elevation of node kr directly drops to rd, and node k f The altitude drops directly to fd.

[0038] Specifically, when searching for the nearest node K or the nearest drill hole D on the surface DTM for the ore body roof and floor DTM node k, the search and judgment method is as follows:

[0039] For an ellipse centered at node k, the inequality for determining whether node K or borehole D falls within the ellipse is:

[0040]

[0041] Where: x0, y0 are the horizontal coordinates of node k;

[0042] x, y are the horizontal coordinates of node K or the nearest drill hole D;

[0043] r is the length of the north-south axis of the ellipse, i.e. the search radius;

[0044] f is the ratio of the east-west axis to the north-south axis; if f = 1, the ellipse is a circle

[0045] r、f×r≤R max , R max is the maximum search radius;

[0046] K or D that satisfies the above inequality, add k search points / drilling sets, the search radius starts from 0 and gradually increases to the maximum search radius R max ;

[0047] When the search radius is 0, the search judgment inequality is:

[0048] |x-x0|+|y-y0|<∈

[0049] Where: ∈ is a very small positive number;

[0050] When the search radius is 0, the search hits point K / borehole D, then it is considered that point k overlaps with point K or k is on borehole D. At this time, only K / D is retained in the search set.

[0051] Compared with previous traditional geological modeling and implicit geological modeling technologies / methods, this application automatically and quickly calculates the elevations of all nodes on the top and bottom surfaces of the ore body within the modeling range based on the drill hole grade information revealed by geological exploration drilling and combined with the industrial grade. At the same time, it takes into account the influence of surface topography on weathering and elution mineralization, and corrects the elevations of each node on the top and bottom of the ore body. It fully considers the mineralization principle and occurrence characteristics of weathering crust elution type deposits, greatly improving the accuracy and reliability of engineering geological three-dimensional modeling.

[0052] The present invention has the following technical features: 1. The present invention is an implicit geological modeling technology / method, and the geological modeling process discards many tedious drawing operations in the traditional geological modeling process; 2. The present invention greatly simplifies the implicit modeling process based on the mineralization principle and occurrence characteristics of weathering crust elution type deposits, and the algorithm is simple to implement; 3. The geological modeling calculation process of the present invention fully considers the influence and control of the mining area surface on the ore body mineralization process and the valuable element enrichment process, and solves the basic problem that traditional geological modeling and general implicit geological modeling do not conform to the characteristics of weathering crust elution type deposits. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] In order to more clearly illustrate the technology / method of the present invention, a brief introduction to the accompanying drawings is given below.

[0054] Figure 1 Schematic diagram of the operation flow of the technology / method of the present invention;

[0055] Figure 2 Schematic diagram of the initial DTM nodes of the roof and floor of the ore body within the mining area of ​​the present invention;

[0056] Figure 3 The initial DTM surface of the roof and floor of the ore body within the mining area of ​​the present invention;

[0057] Figure 4 A schematic diagram of the process for generating the initial DTM surface of the roof and floor of the ore body within the mining area of ​​the present invention;

[0058] Figure 5 This is a schematic diagram of the surface DTM surface within the mining area of ​​the present invention;

[0059] Figure 6 This is an example of production prospecting drilling data within the mining area of ​​the present invention;

[0060] Figure 7 Schematic diagram of the spatial position relationship between the production prospecting drill holes and the surface DTM within the mining area of ​​the present invention;

[0061] Figure 8 This is a flow chart for calculating the DTM surface node elevations of the ore body roof and floor within the mining area of ​​the present invention;

[0062] Figure 9 Schematic diagram of the calculation results of the DTM surface of the ore body roof and floor within the mining area of ​​the present invention;

[0063] Figure 10 It is a cross-sectional schematic diagram of the coincidence relationship between the calculated results of the DTM surface of the ore body roof and floor within the mining area of ​​the present invention and the surface DTM surface and drilling data;

[0064] Figure 11 This is a solid model diagram of the actual ore-bearing body within the mining area of ​​the present invention. DETAILED DESCRIPTION

[0065] The following is a clear and complete description of the technical solutions in the embodiments of the present invention, in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts are within the scope of protection of the present invention.

[0066] See also Figure 1The present invention provides a method for rapid implicit geological modeling of an ore body of a weathering crust elution type deposit, comprising the following steps:

[0067] Step 1: Determine the scope of ore body modeling in the mining area

[0068] According to the mining area scope and mining conditions, a rectangular geological modeling scope that completely covers the mining area is determined. The rectangular scope is east-west and north-south.

[0069] Step 2: Generate the DTM (Digital Terrain Model) triangulation network of the ore body roof and floor within the modeling range

[0070] Based on the modeling accuracy requirements, select an east-west and north-south triangulation value for the ore body roof and floor DTM surface, that is, the minimum side length value in both directions of the initial DTM surface triangulation. Generate the initial triangulation of the ore body roof and floor DTM surface within the modeling range.

[0071] Step 3: Import the surface DTM model within the modeling range

[0072] Given the availability of surface contours and discrete elevation points, it's very convenient to generate a surface DTM within the modeling range using 3D mining software such as 3DMine, which can be directly imported. Alternatively, a Delaunay triangulation algorithm can be used to complete the modeling of the surface DTM within the modeling range. This method ultimately utilizes all nodes of the surface DTM model, so nodes on the surface contours and discrete elevation points can also be directly imported.

[0073] Step 4: Search the surface DTM node search point set corresponding to the ore body roof and floor DTM surface nodes

[0074] On the horizontal plane, with each node of the DTM triangulation network of the top and bottom plates of the ore body as the center, search for the nearest node on the surface DTM surface in a circular or elliptical manner (if anisotropy in the horizontal direction needs to be considered), and gradually increase the search radius from small to large. When enough nodes on the surface DTM surface are searched within the maximum search radius and the number of points (minimum number of points) required by the design is met, the search is terminated, and the search point set is constructed using the searched surface DTM surface nodes.

[0075] When searching for the search point set of the top and bottom DTM surface nodes of the ore body, if there is a node on the surface DTM surface whose horizontal distance to the current top and bottom DTM surface nodes is so small that it can be regarded as two points overlapping, only this closest point will be included in the search point set to ensure that if there is a surface DTM surface node with a horizontal distance close to zero, there is only this one point in the search point set.

[0076] Step 5: Using the search point set found in step 4, make the first estimation adjustment to the elevation of the DTM surface nodes of the top and bottom plates of the ore body.

[0077] Using the elevation values ​​of each node in the search point set as known values, and the inverse of the power of the distance between each node in the search point set and the current node in the orebody roof and floor DTM as weights, the inverse distance power method is used to perform a first elevation calculation for the corresponding node in the orebody roof and floor DTM, modifying the elevations of all orebody roof and floor DTM nodes. After this calculation, the orebody roof and floor DTM surfaces will fully align with the surface DTM. If there is only one node in the search point set, the inverse distance power method is not used for the estimation calculation, and the elevation of the corresponding node in the orebody roof and floor DTM is directly used to assign the elevation value. This means that when the distance is 0, the weight is infinite.

[0078] After this step, the elevations of each node on the orebody roof and floor DTM are adjusted to align with the surface DTM surface. The accuracy of this alignment is related to the mesh size of the orebody roof and floor DTM triangulation determined in step 2. If the alignment accuracy does not meet the requirements, return to step 2 and reduce the triangle side lengths of the initial DTM surface before continuing.

[0079] Step 6: Normalize and import drilling exploration data in the mining area

[0080] Using drilling exploration data and industrial grade values, the top and bottom depth points of the ore body are determined in three-dimensional space for each borehole. Using these top and bottom depth points, the sinking depth values ​​of the top and bottom of the ore body relative to the borehole mouth are calculated.

[0081] The exploration data for all exploratory boreholes within the mining area was standardized and imported into Excel. This information included: borehole number, XYZ coordinates of the borehole mouth, maximum hole depth, starting depth of the ore body, ore body thickness, reason for hole termination, sampling depth and grade value for each sample, and borehole inclination data (inclination depth, azimuth, and inclination). Since production exploration boreholes are drilled vertically downward in actual mining operations with little deviation, inclination data is generally not required. In other words, each borehole is directly identified as having only inclination data at the hole mouth: inclination depth 0 meters (i.e., at the hole mouth), azimuth 0 degrees, and inclination -90 degrees (vertically downward).

[0082] The initial and final burial depths of the ore body for each drill hole are determined based on the industrial grade and the initial depth and thickness of the ore body interpreted in the original drill hole. The initial burial depth of the ore body for a particular drill hole is the depth that the ore body roof has sunk relative to the ground surface at that level, while the final burial depth is the depth that the ore body floor has sunk relative to the ground surface at that level.

[0083] Step 7: Search for the affected drill hole sets of the ore body roof and floor DTM surface nodes

[0084] In the horizontal direction, with each node of the DTM triangulation network of the top and bottom plates of the ore body as the center, search for the nearest drill hole in the form of a circle or ellipse (if anisotropy in the horizontal direction needs to be considered), and gradually increase the search radius from small to large. When enough drill holes are found within the maximum search radius and the number of drill holes required by the design (minimum number of drill holes) is met, the search is stopped and the searched drill holes are used to construct the search drill hole set of the node.

[0085] When searching for the search drill hole set of the top and bottom DTM surface nodes of the ore body, if there is a drill hole whose horizontal distance from the current top and bottom DTM surface nodes is small enough to be considered as the node being on a certain drill hole, that is, a certain node just overlaps with a certain drill hole position in the horizontal direction, that is, when the search radius is 0, the search hits a drill hole, then the search is stopped and only this drill hole is added to the search drill hole set to ensure that if the node is on a drill hole, there is only this drill hole in the search drill hole set.

[0086] Step 8: Using the search drill hole set found in step 7, perform a second estimation adjustment on the elevation of the DTM surface nodes of the ore body top and bottom plates.

[0087] For each node in the ore body roof DTM, the search drillholes found in the search set are used. The starting ore body depth of each drillhole is used as a known value. The inverse power of the horizontal distance between each drillhole in the search drillhole set and the current node in the ore body roof is used as the weight. The inverse distance power method is used to estimate the elevation subsidence value of the ore body roof DTM node. If there is only one drillhole in the drillhole set, the starting ore body depth of that drillhole is directly used to correct the elevation subsidence of the estimated point, which means that the distance is 0 and the weight is infinite.

[0088] For each node in the orebody floor DTM, the final orebody depth of each drillhole is used as a known value. The inverse distance power method is used to estimate the elevation subsidence value of the orebody floor DTM node. If there is only one drillhole in the drillhole set, the final orebody depth of that drillhole is directly used to correct the elevation subsidence of the estimated point.

[0089] After Step 8, the elevation of each node on the orebody roof and floor DTM is adjusted to the orebody roof and floor depth. The degree of agreement between the orebody roof and floor DTM depth at each location and the orebody roof and floor depth in each drillhole is related to parameters such as the drillhole search parameters and the power of the inverse distance power method. If the orebody roof and floor DTM surface accuracy after elevation correction does not meet the requirements, return to Step 7 and adjust the search and calculation parameters before continuing.

[0090] Step 9: Use the mining area boundary to cut the ore body top and bottom DTM surface to obtain the final ore body top and bottom DTM model.

[0091] Using the actual boundary of the mining area, a closed polygon of the mining area range is constructed. Using this closed polygon, the DTM surfaces of the ore body roof and floor within the mining area are cut, and finally the actually available DTM surfaces of the ore body roof and floor in this mining area are obtained, completing the implicit geological modeling of the ore body in the weathered crust eluvial deposit within this mining area.

[0092] Specifically, in step 2, the method for generating the initial triangular network of the DTM surfaces of the ore body roof and floor is as follows:

[0093] According to the modeling range determined in step 1 and the east-west and north-south triangular network densities determined in step 2, nodes arranged in a rectangular array of a rows and b columns are generated.

[0094]

[0095] Where: a is the number of node rows

[0096] b is the number of node columns

[0097] trunc() represents truncating and taking the integer

[0098] LenN is the north-south length of the modeling range

[0099] EdgeN is the side length of the north-south triangle <00002​​​​​​​​​​​​​​​​​​​​​​​​In step 3, since the arrangement of the surface contour nodes and discrete elevation points on the horizontal plane does not have a regular pattern of aligned rows and columns, the method in step 2 cannot be used to generate the surface DTM model. Instead, the Delaunay triangulation algorithm can be used. Of course, for simplicity, the surface DTM model can be directly generated using 3DMine or other software and then imported for use. The Delaunay triangulation algorithm can also be used to generate the triangulated network of the top and bottom plates of the ore body in step 2, but since the nodes are neatly arranged in rows and columns on the horizontal plane, the simpler and more effective algorithm mentioned above can be used to complete the triangulation.

[0109] In steps 4 and 7, when searching for the nearest node K or the nearest drill hole D (horizontal coordinates (x, y)) on the surface DTM for the ore body roof and floor DTM node k (horizontal coordinates (x0, y0)), the search and judgment method is as follows:

[0110] For an ellipse with its center at point k(x0,y0), the inequality for determining whether node K or borehole D falls within the ellipse is:

[0111]

[0112] Where: r is the length of the north-south axis of the ellipse, that is, the search radius

[0113] f is the ratio of the east-west axis length to the north-south axis length. If f = 1, the ellipse is a circle

[0114] r、f×r≤R max , R max is the maximum search radius

[0115] The search radius starts from 0 and increases gradually to the maximum search radius R. max ,When the number of search points / drilling set elements reaches the designed number, the search is terminated early.

[0116] In particular, when the search radius is 0, the search judgment inequality is:

[0117] |x-x0|+|y-y0|<∈

[0118] Where: ∈ is a very small positive number, which can be specifically set to 0.1;

[0119] When the search radius is 0, the search hits point K / hole D, then let k overlap with point K or k is on hole D. At this time, only this K / D is retained in the search set.

[0120] In step 5, the distance power inverse method is used to calculate the elevation value of each node k on the top and bottom plate DTM surface of the ore body for the first adjustment. The calculation formula is as follows:

[0121]

[0122] Where: Z is the elevation value of point k adjusted for the first time

[0123] z1, z2, z3…z m The elevation value of the point in the search point set

[0124] d1, d2, d3…d m is the horizontal distance between point 1, point 2, point 3…point m and point k in the search point set

[0125] Q is the power, which takes the value of 2; m represents the number of points in the search point set;

[0126] If d1, d2, d3…d m Existence d h =0, then Z is directly assigned to z h ;d h is the horizontal distance between points h and k in the search point set, h≤m. If there is only one point in the search point set and its elevation is z, then z is also directly assigned to Z.

[0127] In step 8, the distance power inverse method is used to calculate the k of each node on the DTM surface of the ore body roof. r and each node k on the DTM surface of the ore body bottom f The calculation formulas for elevation correction value (sinking value) are as follows:

[0128]

[0129] Where: R Δz is the node k on the DTM surface of the ore body roof r The elevation reduction value of the second elevation adjustment, that is, the sinking value

[0130] rd1, rd2, rd3…rd n To search for the ore body roof depth value of borehole 1, borehole 2, borehole 3... borehole n in the borehole set, that is, the distance between the ore body roof and the hole mouth in the borehole

[0131] F Δz is the node k on the DTM surface of the ore body bottom f The elevation reduction value of the second elevation adjustment, that is, the sinking value

[0132] fd1, fd2, fd3…fd n The depth of the ore body bottom plate of the drilling hole 1, drilling hole 2, drilling hole 3... drilling hole n in the drilling set is searched, that is, the distance between the ore body bottom plate and the hole mouth in the drilling hole

[0133] d1, d2, d3…d nTo search for the horizontal distance between the opening of borehole 1, borehole 2, borehole 3, etc. and point k in the borehole set

[0134] p is the power, which can be 3 here.

[0135] d1, d2, d3…d n is the horizontal distance between the drilling hole opening in the search drilling set and point k

[0136] p is the power, which can be 3 here.

[0137] If d1, d2, d3…d n Existence d t =0,d t The horizontal distance between the hole opening of the searched hole set and point k, node k r Elevation directly reduces rd t , and node k f Elevation directly reduces fd t If there is only one drill hole in the search drill hole set, and the top plate depth of the ore body in the drill hole is rd and the bottom plate depth is fd, then the elevation of node kr directly drops to rd, and node k f The altitude drops directly to fd.

[0138] Specifically, in step 9, weathering crust eluvial deposits are often continuously mineralized in areas with mineralization conditions, and when the mine actually organizes production, the mining area will be divided according to the original mountain site conditions. The top and bottom DTM model of the ore body in the mining range of a certain mining area uses the mining range demarcation line of the mining area to cut the top and bottom DTM of the ore body within the rectangular modeling range obtained in the previous step, and finally obtain the top and bottom DTM model of the ore body within the mining range of the mining area. The function of cutting the DTM surface by the mining range demarcation line is provided by many software such as 3DMine and Surpac, and the present invention will no longer explain the function implementation process.

[0139] Compared with previous traditional geological modeling and implicit geological modeling technologies / methods, the present invention fully considers the mineralization principle and occurrence characteristics of weathering crust eluvial deposits, and combines the actual conditions and needs of on-site production exploration and mining work in such resource mines. The present invention has the following technical features: 1. The present invention is an implicit geological modeling technology / method, and the geological modeling process discards many tedious drawing operations in the traditional geological modeling process; 2. The present invention greatly simplifies the implicit modeling process based on the mineralization principle and occurrence characteristics of weathering crust eluvial deposits, and the algorithm is simple to implement; 3. The geological modeling calculation process of the present invention fully considers the influence and control of the mining area surface on the ore body mineralization process and the valuable element enrichment process, solving the basic problem that traditional geological modeling and general implicit geological modeling do not conform to the characteristics of weathering crust eluvial deposits.

[0140] In summary, the present invention realizes a simplified and rapid implicit geological model technology / method based on the unique mineralization principles and occurrence characteristics of weathering crust eluvial deposits, which not only solves the problem of simplifying the operation process and algorithm, but also conforms to the actual geological characteristics and production needs.

[0141] The present invention will be further described below with reference to specific examples:

[0142] The method for rapid implicit geological modeling of an ore body of a weathering crust elution-type deposit provided in an embodiment of the present invention specifically includes the following steps:

[0143] Step S101: According to the mining area division of the weathering crust eluvial deposit, a rectangular area with an east-west and north-south direction that can completely cover the mining area is selected and determined as the modeling range: the southwest corner coordinates are (19,14,0), the northeast corner coordinates are (269,214,0), the east-west length is 250 meters, and the north-south length is 200 meters;

[0144] Step S201: Based on the rectangular modeling range, generate a node matrix of 21 rows and 26 columns with a total of 546 nodes, with a node step distance of 10 meters in the east-west and north-south directions. Figure 2 ;

[0145] Step S202: Generate the initial triangulation of the DTM surface of the top and bottom plates of the mining area using the node matrix, see the attached Figure 3 The generation process of the initial triangulated network is shown in the attached Figure 4 ;

[0146] Step S301: In 3DMine software, use the contour lines within the mining area to generate the surface DTM surface within the modeling range. Figure 5 ;

[0147] Step S401: The DTM surface of the roof and floor of the mining area ore body generated in step S202 is generated based on a node matrix of 21 rows × 26 columns with a total of 546 nodes. With each node as the center, a node on the surface DTM is searched for that node in the horizontal direction to construct a search point set for each node of the roof and floor DTM surface of the mining area ore body. The search radius is gradually increased until the search point set meets the minimum number of points.

[0148] Step S501: Using the surface DTM surface node exploration point set of each node of the DTM surface of the ore body roof and floor in the mining area, a first estimation adjustment is performed on the elevation of the node. The estimation method uses the inverse distance power method. After the elevation estimation adjustment, the elevation of the DTM surface of the ore body roof and floor in the mining area is basically consistent with the elevation of the surface DTM surface in the mining area.

[0149] Step S601: Standardize and catalog the geological drilling data within the mining area in an EXCEL file. Figure 6 ;

[0150] Step S602: Read the data in the normalized EXCEL file of geological drilling data and display the ore body line segment of each drill hole in the 3D graphics area - the sample segment in the drill hole where the ore body grade reaches the cut-off grade (or industrial grade). At the same time, the surface DTM surface is displayed in the 3D graphics area. The spatial position relationship between the geological drilling hole and the surface DTM is shown in the attached figure. Figure 7 ;

[0151] Step S701: In the horizontal direction, search for influencing drill hole datasets for each node of the roof and floor DTM surface of the mining area. When searching for drill holes, determine the appropriate maximum search radius based on the geological drill hole grid. The search radius increases from 0 to 1. When the number of drill holes found meets the minimum number of drill holes, stop searching.

[0152] Step S801: Using the distance power inverse method, the starting depth and thickness of the ore body in each drill hole of the ore body top and bottom plate DTM surface node of each mining area are concentrated to perform a second elevation correction calculation for the ore body top and bottom plate DTM nodes. The top plate DTM node uses the starting depth value, and the bottom plate DTM node uses the starting depth and ore body thickness value to participate in the distance power inverse method calculation, that is, the top and bottom plate nodes are respectively sunken based on the original elevation. The calculation process of this step is shown in the attached figure. Figure 8 After the second elevation correction calculation, the top and bottom DTM of the ore body are the DTM models of the top and bottom of the ore body within the modeling range of the mining area. Figure 9 ;

[0153] Step S901: Based on the specific on-site production conditions such as the mining area topography, use a multi-segment line to draw the boundary of the ore body that can be mined in the mining area, and use the closed polyline to cut the DTM of the top and bottom plates of the ore body to finally obtain the DTM surface of the top and bottom plates of the ore body that can be mined in the mining area, thus completing the implicit geological modeling of the ore body that can be mined in the mining area. Figure 9 ;

[0154] Step S902: Cut the longitudinal section along a certain direction and verify on the longitudinal section that the implicit modeling result of the ore body is very consistent with the surface DTM and geological drilling data. Figure 10 .

[0155] The above description is only a specific embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the invention should be included in the scope of protection of the present invention.

Claims

1. A method for rapid implicit geological modeling of ore bodies of weathering crust elution type deposits, characterized in that: include: Generate the DTM triangulated network of the roof and floor of the ore body according to the modeling scope and geological modeling accuracy requirements; import the surface DTM model within the modeling scope; Taking each node k on the DTM surface of the ore body roof and floor as the center, increase the search radius from small to large, and search the node closest to the node k on the surface DTM surface in the horizontal direction to form the search point set of node k; Using the search point set of node k and the inverse distance power method, the elevation of each point in the search point set is regarded as a known value. The elevation of node k is adjusted for the first time, and the elevation of each node on the ore body roof and floor DTM is adjusted to fit with the surface DTM surface. If the fitting accuracy does not meet the design requirements, the triangulation density is adjusted and the calculation is repeated again after returning to the step of generating the triangulation of the ore body roof and floor DTM surface. Using drilling exploration data and industrial grade values, determine the top and bottom depth points of the ore body in each drill hole in three-dimensional space; Using the top and bottom depth points of the ore body on the drill hole, calculate the sinking depth values ​​of the top and bottom of the ore body relative to the drill hole opening; Taking each node k on the DTM surface of the ore body roof and floor as the center, increase the search radius from small to large, and search for the drill hole closest to the node in the horizontal direction. Stop searching when the number of drill holes required by the design is reached, and use the searched drill holes to construct the search drill hole set of the node; Using the search drill hole set, a similar method to searching for surface DTM nodes is used to perform a second estimation adjustment on the elevation of each node on the ore body roof and floor DTM after the first estimation adjustment. The elevation of each node on the ore body roof and floor DTM is adjusted to match the buried depth of the ore body roof and floor. If the matching accuracy does not meet the project requirements, adjust the search parameters and return to the search drill hole set construction step to recalculate. Using the mining area boundary, the DTM surface of the ore body top and bottom plates after the second valuation adjustment is cut to obtain the final DTM model of the ore body top and bottom plates, completing the rapid implicit geological modeling of the ore body model in the mining area.

2. The method for rapid implicit geological modeling of an ore body according to claim 1, characterized in that: For the node kr on the DTM surface of the ore body roof, based on the obtained search drill hole set, the starting depth of the ore body on each drill hole is taken as a known value, and the inverse power method of distance is used to estimate the elevation drop value of the node kr, and the node kr elevation is adjusted accordingly; For the node kf on the DTM surface of the ore body bottom plate, according to the obtained search drill hole set, the starting depth of the ore body on each drill hole + the ore body thickness is taken as a known value, and the inverse power method of distance is used to estimate the elevation drop value of the node kf, and the node kf elevation is adjusted accordingly.

3. The method for rapid implicit geological modeling of an ore body according to claim 2, characterized in that: For the node kr on the ore body roof DTM surface, if there is only one drill hole in the search drill hole set, the ore body starting depth value of the drill hole is directly used to perform subsidence adjustment calculation on the elevation of the current node; For the node kf on the DTM surface of the ore body bottom plate, if there is only one drill hole in the search drill hole set, the ore body starting depth + ore body thickness value of the drill hole is directly used to perform the subsidence adjustment calculation for the elevation of the current node.

4. The method for rapid implicit geological modeling of an ore body according to claim 1, wherein: According to the mining area scope and mining conditions, a rectangular geological modeling range that completely covers the mining area is determined. The rectangular geological modeling range is east-west and north-south.

5. The method for rapid implicit geological modeling of an ore body according to claim 4, characterized in that: According to the requirements of geological modeling accuracy, an east-west and north-south triangulation value is selected for the DTM surface of the ore body roof and bottom plates to generate the triangulation of the DTM surface of the ore body roof and bottom plates within the modeling range.

6. The method for rapid implicit geological modeling of an ore body according to claim 1, characterized in that: If there is only one node in the search point set, the elevation of the point in the search point set is directly used to assign elevation values ​​to the corresponding nodes of the ore body roof and floor DTM.

7. The method for rapid implicit geological modeling of an ore body according to claim 1, characterized in that: Using surface data such as contour lines and elevation points, a DTM model of the mining area surface is constructed using 3D mining software.

8. The method for rapid implicit geological modeling of an ore body according to any one of claims 1 to 7, characterized in that: The distance power inverse method is used to calculate the elevation value of each node k on the DTM surface of the ore body roof and floor for the first time. The calculation formula is as follows: Where: Z is the elevation value of point k adjusted for the first time z1, z2, z3…z m The elevation value of the point in the search point set d1, d2, d3…d m is the horizontal distance between the point in the search point set and point k Q is the power, which is 2; If d1, d2, d3…d m Existence d h =0, then Z is directly assigned to z h ; If there is only one point in the search point set and its elevation is z, then z is also directly assigned to Z.

9. The method for rapid implicit geological modeling of an ore body according to claim 8, characterized in that: The distance power inverse method is used to calculate the elevation correction values ​​of each node kr on the ore body roof DTM surface and each node kf on the ore body bottom DTM surface. The calculation formulas are as follows: Where: R Δz The elevation reduction value of the second elevation adjustment of node kr on the DTM surface of the ore body roof, that is, the subsidence value; rd1, rd2, rd3…rd n The depth of the ore body roof of the drill hole in the search drill hole set, that is, the distance between the ore body roof and the hole mouth in the drill hole; F Δz The elevation reduction value of the node kf on the DTM surface of the ore body bottom plate for the second elevation adjustment, that is, the subsidence value; fd1, fd2, fd3…fd n The depth of the ore body floor of the drill hole in the search drill hole set, that is, the distance between the ore body floor and the hole mouth in the drill hole; d1, d2, d3…d n is the horizontal distance between the drilling hole opening in the search drilling set and point k p is the power, and its value is 3; If d1, d2, d3…d n Existence d t =0, the node kr elevation is directly reduced to rd t , and the node kf elevation directly reduces fd t If there is only one borehole in the search borehole set, and the top plate depth of the ore body on this borehole is rd and the bottom plate depth is fd, then the elevation of node kr directly drops by rd, and the elevation of node kf directly drops by fd.

10. The method for rapid implicit geological modeling of an ore body according to any one of claims 1 to 7, characterized in that: When searching for the nearest node K or the nearest drill hole D on the surface DTM for the ore body roof and floor DTM node k, the search and judgment method is as follows: For an ellipse centered at node k, the inequality for determining whether node K or borehole D falls within the ellipse is: Where: x0, y0 are the horizontal coordinates of node k; x, y are the horizontal coordinates of node K or the nearest drill hole D; r is the length of the north-south axis of the ellipse, i.e. the search radius; f is the ratio of the east-west axis to the north-south axis; if f = 1, the ellipse is a circle K or D that satisfies the above inequality, add k search points / drilling sets, the search radius starts from 0 and gradually increases to the maximum search radius R max .

Images