Ore body rapid implicit geological modeling method for weathering crust elution deposit

Abstract

The invention discloses an ore body rapid implicit geological modeling method for a weathering crust eluviation type ore deposit. The method comprises the steps of generating an ore body top and bottom plate surface triangulation network in a modeling range, importing an earth surface DTM model in the modeling range, constructing an earth surface DTM surface node search point set of an ore body top and bottom plate DTM node, performing estimation adjustment on the elevation of the ore body top and bottom plate DTM surface node, and the like. According to the method, a series of complex manual operations such as interpretation of the ore body on the section and connection of the ore body between the sections in a traditional geological modeling method are omitted, the geological modeling process is automatically completed, manual drawing operation is not needed, the metallogenic principle of the illuvial deposit is reflected, the algorithm is simple and clear, calculation is rapid, exploration data are updated, and after the ore body interpretation standard is changed, the geological modeling efficiency is greatly improved. And ore body reconstruction is simple and rapid.

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 ore deposit. Background Art

[0002] Three-dimensional geological model is the foundation of modern digital mine. Three-dimensional geological modeling is to generate three-dimensional models of geological objects by using computer graphics technology after summarizing various geological information and geological interpretation results. Implicit modeling is a parametric spatial interpolation technology. One of the keys of this technology is to realize automatic geological modeling by using spatial interpolation algorithm. At present, the reliability and practicality of implicit modeling technology have been verified in geological modeling of various types of mineral deposits. It is also an important development and promotion and application direction in the field of geological modeling now and in the future.

[0003] Weathering crust eluvial deposits are a type of mineral resource with very obvious characteristics. Affected by weathering 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 produced in layers, tilted from the ridge to the foot of the mountain, and the ore body occurrence elevation is generally higher than the local lowest erosion datum level; 3. The occurrence of the ore body is jointly affected and controlled by the strata and the surface topography. In general, weathering crust eluvial deposits have layered ore bodies and simple morphology, which are controlled by the weathering crust strata and the surface topography. At present, implicit geological modeling technology mainly focuses on exploration data obtained from drilling and other exploration projects. While studying complex algorithms and calculation processes, it completely ignores the joint influence and control of surface topography and strata on the morphology of ore bodies during the mineralization of weathering crust elution-type deposits, resulting in serious deviations between the results of implicit geological modeling and the actual occurrence of ore deposits. Therefore, the current software that supports implicit geological modeling is generally insufficient in practicality for geological modeling of weathering crust elution-type deposits. Summary of the invention

[0004] The main purpose of the present invention is to provide a method for rapid implicit geological modeling of ore bodies of weathering crust elution type deposits. This method automatically and rapidly 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, the elevations of each node on the top and bottom surfaces of the ore body are corrected taking into account the influence of the surface topography on the weathering elution mineralization, fully considering the mineralization principle and occurrence characteristics of the weathering crust elution type deposits, and greatly improving the accuracy and reliability of engineering geological three-dimensional modeling.

[0005] To this end, the invention provides a method for rapid implicit geological modeling of an ore body of a weathering crust elution type ore deposit, comprising:

[0006] Generate the DTM triangulated network of the top and bottom plates 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 top and bottom DTM surface of the ore body as the center, increasing the search radius from small to large, searching for the node closest to node k on the surface DTM surface in the horizontal direction, forming a search point set for 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, and the elevation of node k is adjusted for the first time. The elevation of each node on the top and bottom plates of the ore body 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 triangulation is returned to the step of generating the triangulation of the top and bottom plates of the ore body DTM surface for recalculation.

[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 borehole, calculate the sinking depth values ​​of the top and bottom plates of the ore body on the borehole relative to the borehole opening;

[0011] Taking each node k on the DTM surface of the top and bottom plates of the ore body as the center, increase the search radius from small to large, 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] By searching the borehole set, a similar method to searching the surface DTM nodes is adopted to make a second estimation adjustment on the elevation of each node on the top and bottom plates of the ore body after the first estimation adjustment, and the elevation of each node on the top and bottom plates of the ore body is adjusted to match the buried depth of the top and bottom plates of the ore body. If the matching accuracy does not meet the engineering requirements, the borehole search parameters are adjusted and the calculation is returned to the step of building the search borehole set 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, according to 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 thickness of the ore body 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 borehole in the search borehole set, the ore body starting depth value of the borehole is directly used to perform a sinking adjustment calculation on 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 borehole in the search borehole set, the ore body starting depth + ore body thickness value of the borehole 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, and the rectangular geological modeling range is east-west and north-south.

[0019] Specifically, according to the requirements of geological modeling accuracy, an east-west and north-south triangulation grid value is selected for the DTM surface of the top and bottom plates of the ore body, and a triangulation network of the DTM surface of the top and bottom plates of the ore body 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 elevations to corresponding nodes of the ore body top and bottom plate 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 adjusted for the first time on the DTM surface of the top and bottom plates of the ore body. The calculation formula is as follows:

[0023]

[0024] Where: Z is the elevation value adjusted to point k 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, and its value is 2;

[0028] If d1, d2, d3…d m Existence 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 sinking value;

[0032] rd1, rd2, rd3…rd n The depth value 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 sinking value;

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

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

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

[0037] If d1, d2, d3…d n Existence t =0, node k r Elevation directly reduces rd t , and node k f 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 in the borehole 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 directly drops fd.

[0038] Specifically, when searching for the nearest node K or the nearest borehole D on the surface DTM surface 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 borehole 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 length to the north-south axis length; 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 according to the step size. max ;

[0047] When the search radius is 0, the search decision 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 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 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 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 drawings is given below.

[0054] Figure 1 It is a schematic diagram of the operation flow of the technology / method of the present invention;

[0055] Figure 2 It is a 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 of 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 It 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 It is a schematic diagram of the spatial position relationship between the production prospecting drill hole and the surface DTM within the mining area of ​​the present invention;

[0061] Figure 8 It is a flow chart of calculating the elevation of the DTM surface nodes of the top and bottom plates of the ore body within the mining area of ​​the present invention;

[0062] Fig. 9 It is a 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] Fig.10 It is a cross-sectional schematic diagram of the matching relationship between the calculation results of the DTM surface of the roof and floor of the ore body within the mining area of ​​the present invention and the surface DTM surface and drilling data;

[0064] Fig.11 It is a solid model diagram of the actual ore-carrying 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 drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

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

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

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

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

[0070] According to the modeling accuracy requirements, select an east-west and north-south triangulation value for the DTM surface of the ore body top and bottom plates, that is, the minimum side length value in two directions of the initial DTM surface triangulation. Generate the initial triangulation of the DTM surface of the ore body top and bottom plates within the modeling range.

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

[0072] Under the condition of surface contour lines and discrete elevation points, it is very convenient to generate the surface DTM surface within the modeling range with the help of mining 3D software such as 3DMine, which can be directly imported. Of course, the Delaunay triangulation algorithm can also be used to complete the surface DTM surface modeling within the modeling range. The present invention ultimately utilizes all the nodes of the surface DTM model, so the nodes on the surface contour lines and discrete elevation points can also be directly imported.

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

[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) that meet the design requirements are found, 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 of the ore body is small enough to be regarded as two points overlapping, only this nearest point will be added to 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] The elevation values ​​of each node in the search point set of each node of the top and bottom plate DTM of the ore body are used as known values, and the inverse of the power of the distance between each node in the search point set and the current node of the top and bottom plate DTM of the ore body is used as the weight. The inverse power method of distance is used to calculate the first elevation value of the elevation of the corresponding node of the top and bottom plate DTM of the ore body, and the elevation of all the nodes of the top and bottom plate DTM of the ore body is modified. After the calculation is completed, the surface of the top and bottom plate DTM of the ore body will fully fit with the surface DTM. If there is only one node in the search point set, the inverse power method of distance is not used for valuation calculation, and the elevation of the point in the search point set is directly used to assign elevation to the corresponding node of the top and bottom plate DTM of the ore body, that is, when the distance is 0, the weight is infinite.

[0078] After this step, the elevation of each node on the top and bottom plates of the ore body is adjusted to fit with the surface DTM surface, and the fitting accuracy is related to the size of the triangulated network of the top and bottom plates of the ore body determined in step 2. If the fitting accuracy does not meet the requirements, return to step 2, reduce the triangle side length of the initial DTM surface and continue.

[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 the top and bottom depth points of the ore body, the sinking depth values ​​of the top and bottom of the ore body relative to the borehole mouth are calculated.

[0081] The exploration data of all exploration boreholes in the mining area are standardized and sorted in Excel and then imported. The imported information includes: borehole number, XYZ coordinates of borehole mouth, maximum hole depth, starting depth of ore body, ore body thickness, reason for ending hole, sampling depth and grade value of each sample, and borehole inclination data (inclination depth, azimuth, inclination). Since in actual mining work, production exploration boreholes are drilled vertically downward and basically without deviation, inclination measurement is generally not required, that is, each borehole is directly identified as having only inclination measurement data at the hole mouth: inclination depth 0 meters (i.e., hole mouth), azimuth 0 degrees, inclination -90 degrees (vertically downward).

[0082] According to the industrial grade value, the initial depth and thickness of the ore body interpreted on the original borehole are comprehensively referenced to determine the initial and final burial depths of the ore body on each borehole. It can be considered that the initial burial depth of the ore body on a certain borehole is the sinking depth of the ore body top plate at that horizontal position relative to the ground surface; the final burial depth of the ore body is the sinking depth of the ore body bottom plate at that horizontal position relative to the ground surface.

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

[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 the 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 to meet the design requirements (minimum number of drill holes), the search is stopped and the searched drill holes are used to build the search drill hole set of the node.

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

[0086] Step 8: Use the searched drill hole set found in step 7 to make a second estimation adjustment on the elevation of the DTM surface nodes of the top and bottom plates of the ore body.

[0087] In the search borehole set for each node in the ore body roof DTM surface, the starting burial depth of each borehole is taken as a known value, and the inverse of the power of the horizontal distance between each borehole in the search borehole set and the current node of the ore body roof is taken as the weight. The inverse distance power method is used to estimate the elevation subsidence value of the ore body roof DTM surface node. If there is only one borehole in the borehole set, the starting burial depth of the ore body of the borehole is directly used to correct the elevation subsidence of the estimated point, that is, the distance is 0 and the weight is infinite.

[0088] The final buried depth of the ore body of each drill hole in the searched drill hole set is taken as a known value, and the elevation subsidence value of the node of the ore body bottom plate DTM surface is estimated by using the inverse power method of distance. If there is only one drill hole in the drill hole set, the final buried depth of the ore body of the drill hole is directly used to correct the elevation subsidence of the estimated point.

[0089] After step 8, the elevation of each node on the ore body roof and floor DTM is adjusted to the ore body roof and floor depth. The depth of the ore body roof and floor DTM at each location is consistent with the depth of the ore body roof and floor in each borehole, which is related to the borehole search parameters, the power of the inverse distance power method and other parameters. If the surface accuracy of the ore body roof and floor DTM after the 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 DTM surface of the ore body top and bottom plates to obtain the final DTM model of the ore body top and bottom surfaces.

[0091] Using the actual boundary of the mining area, construct a closed polygon for the mining area. Using this closed polygon, cut the DTM surfaces of the ore body roof and floor within the mining area, and finally obtain the actually available DTM surfaces of the ore body roof and floor in this mining area, completing the implicit geological modeling of the ore body in this weathering crust eluvial deposit within the 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, generate nodes arranged in a rectangular pattern of a rows and b columns.

[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

[0100] LenE is the east-west length of the modeling range

[0101] EdgeE is the side length of the east-west triangle

[0102] ε is a relatively small positive number

[0103] For the a×b nodes arranged in a rectangular pattern, construct a triangular network in the way that for every four closely adjacent nodes on the horizontal plane (the 4 vertices of a small rectangle on the plane), two right triangles are constructed, completing the triangulation of the a×b nodes. For example, for four nodes on the horizontal plane:

[0104] P1: The j-th node in the i-th row

[0105] P2: The (j + 1)-th node in the i-th row

[0106] P3: The j-th node in the (i + 1)-th row

[0107] P4: The (j + 1)-th node in the (i + 1)-th row, where 1 ≤ i < a and 1 ≤ j < b. Two triangles are formed, sorted counterclockwise. The nodes of these two triangles are respectively: (P1 P2 P3), (P2 P4 P3).

[0108] In step 3, since the arrangement of the surface contour nodes and discrete elevation points on the horizontal plane does not have a regular arrangement of rows and columns, the method in step 2 cannot be used to generate the surface DTM model. The Delaunay triangulation algorithm can be used to complete it. Of course, for simplicity, you can directly use 3DMine or other software to generate the surface DTM model and then import it 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 aforementioned simpler and more effective algorithm can be used to complete the triangulation.

[0109] In step 4 and step 7, when searching for the nearest node K or the nearest borehole 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, i.e. 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] K or D that satisfies the above inequality, add k's search point / drilling set. The search radius starts from 0 and gradually increases to the maximum search radius R in steps. max , when the number of elements in the search point / drilling set 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 taken as 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 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 adjusted for the first time on the DTM surface of the top and bottom plates of the ore body. The calculation formula is as follows:

[0121]

[0122] Where: Z is the elevation value adjusted to point k 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 h =0, then Z is directly assigned to z h ;d h is the horizontal distance between point h and point 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, i.e. 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, i.e. the sinking value

[0132] fd1, fd2, fd3…fd n To search for the ore body bottom plate depth value of borehole 1, borehole 2, borehole 3... borehole n in the borehole set, that is, the distance between the ore body bottom plate and the hole mouth in the borehole

[0133] d1, d2, d3…d nTo search for the horizontal distance between the opening of borehole 1, the opening of borehole 2, the opening of borehole 3, etc. and the 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 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 t =0,d t To search for the horizontal distance between the opening of borehole t and point k in the borehole set, node k r Elevation directly reduces rd t , and node k f 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 in the borehole 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 directly drops 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 delineation 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 delineation line is provided by many software such as 3DMine and Surpac, and the present invention will no longer explain the function realization process.

[0139] Compared with the previous traditional geological modeling and implicit geological modeling technology / method, the present invention fully considers the mineralization principle and occurrence characteristics of weathering crust eluvial deposits, and combines the actual situation and needs of on-site production exploration and recovery work of such resource mines. The present invention has the following technical characteristics: 1. The present invention is an implicit geological modeling technology / method, and the geological modeling process discards many cumbersome drawing operations in the traditional geological modeling process; 2. The present invention greatly simplifies the implicit modeling process according to 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, and solves 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 fast 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 in conjunction with specific examples:

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

[0143] Step S101: According to the mining area division of the weathering crust eluvial deposit, a rectangular area in the east-west-north-south direction that can completely cover the mining area is selected 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, a node matrix of 21 rows × 26 columns is generated with a node step distance of 10 meters in the east-west and north-south directions, with a total of 546 nodes. Figure 2 ;

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

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

[0147] Step S401: The DTM surface of the top and bottom plates of the ore body in the mining area 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, the nodes on the surface DTM are searched for the node in the horizontal direction to construct the search point set of the nodes of the top and bottom plates of the ore body in each mining area. The search radius increases gradually 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 top and bottom plates of the ore body in the mining area, the elevation of the node is first estimated and adjusted. The estimation method uses the inverse distance power method. After the elevation estimation adjustment, the elevation of the DTM surface of the top and bottom plates of the ore body in the mining area is basically consistent with the elevation of the surface DTM of 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 borehole in the 3D graphics area - the sample segment in the borehole where the ore body grade reaches the boundary grade (or industrial grade), and display the surface DTM surface 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 the influencing drill hole data set for each node of the DTM surface of the top and bottom plates of the ore body in the mining area. When searching for drill holes, determine the appropriate maximum search radius according to the geological drill hole grid. The search radius increases from 0 to , and stop searching when the number of drill holes searched meets the minimum number of drill holes;

[0152] Step S801: Using the distance power inverse method, the starting depth and thickness of the ore body in each borehole in the drilling search set of the top and bottom plate DTM surface nodes of each mining area are used to perform the second elevation correction calculation for the top and bottom plate DTM nodes of the ore body. 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 sink the calculated values ​​of this step respectively on the basis of 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 is the top and bottom DTM model of the ore body within the modeling range of the mining area. Fig. 9 ;

[0153] Step S901: According to 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, and finally obtain the DTM surface of the top and bottom plates of the mineable ore body in the mining area, completing the implicit geological modeling of the mineable ore body in the mining area. Fig. 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. Fig.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 principle of the invention should be included in the protection scope of the present invention.

Claims

1. A method for rapid implicit geological modeling of ore bodies of weathering crust eluvial deposits, characterized in that: include: Generate the DTM triangulated network of the top and bottom plates 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 top and bottom DTM surface of the ore body as the center, increasing the search radius from small to large, searching for the node closest to node k on the surface DTM surface in the horizontal direction, forming a search point set for 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, and the elevation of node k is adjusted for the first time. The elevation of each node on the top and bottom plates of the ore body 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 triangulation is returned to the step of generating the triangulation of the top and bottom plates of the ore body DTM surface for recalculation. 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 borehole, calculate the sinking depth values ​​of the top and bottom plates of the ore body on the borehole relative to the borehole opening; Taking each node k on the DTM surface of the top and bottom plates of the ore body as the center, increase the search radius from small to large, 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; By using the search borehole set, a similar method to the search for surface DTM nodes is used to perform a second estimation adjustment on the elevation of each node on the top and bottom plates of the ore body after the first estimation adjustment, and the elevation of each node on the top and bottom plates of the ore body is adjusted to match the buried depth of the top and bottom plates of the ore body. If the matching accuracy does not meet the engineering requirements, the search parameters are adjusted and the calculation is returned to the search borehole set construction step for recalculation. 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, according to 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 thickness of the ore body 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 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; For the node kf on the DTM surface of the ore body bottom plate, if there is only one borehole in the search borehole set, the ore body starting depth + ore body thickness value of the borehole 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, characterized in that: According to the scope of the mining area 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 top and bottom plates to generate the triangulation of the DTM surface of the ore body top 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 elevations 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 3D mining software is used to construct a DTM model of the mining area surface.

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 adjusted for the first time on the DTM surface of the top and bottom plates of the ore body. The calculation formula is as follows: Where: Z is the elevation value adjusted to point k 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, and its value is 2; If d1, d2, d3…d m Existence 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 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: 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 sinking value; rd1, rd2, rd3…rd n The depth value 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 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 sinking value; fd1, fd2, fd3…fd n The depth value of the ore body bottom plate of the drill hole in the search drill hole set, that is, the distance between the ore body bottom plate and the hole mouth in the drill hole; d1, d2, d3…d n is the horizontal distance between the hole opening in the search drilling set and point k p is the power, the value is 3; If d1, d2, d3…d n Existence t =0, node k r Elevation directly reduces rd t , and node k f 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 in the borehole 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 directly drops 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 borehole D on the surface DTM surface 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 borehole 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 length to the north-south axis length; 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 according to the step size. max .

Images