Mine digital model construction method based on three-dimensional point cloud database

Abstract

The present invention belongs to the field of mine three-dimensional digital modeling, and disclosed is a mine digital model construction method based on a three-dimensional point cloud database. Mine

Description

Description MINE DIGITAL MODEL CONSTRUCTION METHOD BASED ON THREE-DIMENSIONAL POINT CLOUD DATABASE Technical field The present disclosure relates to a digital mine modeling method based on a three-dimensional point cloud database, and belongs to the technical field of three-dimensional digital mine modeling. Background Art With the continuous exploitation of underground mineral resources, mining depth and intensity are gradually increasing. Once the mine dynamic disasters occur, they will have serious impacts on personnel safety and mine production. Ground pressure observation is the main decision-making basis for solving mining problems at present, and has influence on the decision-making and deployment for major technical problems such as pressure control of roadway, design of mining program, reasonable selection and maintenance of roadway location, prediction, prevention and control of rock burst and outburst, and safe coal mining above confined aquifer. Ground pressure observation is limited by test location and high cost. Therefore, it is difficult to carry out large-scale monitoring in mining areas with ground pressure observation. The development of simulation technology has solved the problem of limited actual measurement and can provide comprehensive and scientific guidance for field practice. One of the keys for simulation is modeling, and the accuracy of the model is an important basis for the reliability of the simulation result. The establishment of large-scale complex models has always been the most important concern of scholars. The models established previously have three shortcomings: 1) the models are small, and mostly models of local mining areas; 2) the models are simple, in which the actual geological structures such as coal seam folds and faults, etc. are ignored; 3) professional modeling software is rarely used in the mine models established previously, and the built-in modeling modules of most simulation software are relatively simple, and cannot meet the modeling requirements under complex geological conditions. The present disclosure employs professional 3D modeling software Rhino, which can reproduce the actual appearance of the coal seam, improve the fitness with the actual field, and meet the requirements of complex large-scale mine modeling. A point cloud model is often obtained directly from measurement, and each point corresponds to a measurement point without other processing means. Therefore, the point cloud model contains a maximum amount of information. However, the information concealed in the point cloud needs to be extracted by other extraction methods, and the process of extracting information in the point cloud is three-dimensional image processing. The idea of three-dimensional point cloud modeling has been applied in many fields. The surface attribute point information of the target object is usually acquired in a non-contact method to obtain three-dimensional point cloud data, but the internal structure of the target object is not finely modeled. At present, in the mining field, three-dimensional laser scanning technology is applied to the above-ground model of opencast mines and side slopes, etc., but there is no detailed research on underground complex modeling. Therefore, the present disclosure introduces in detail a digital modeling method for complex mining inversion and focuses on the application of an idea of three-dimensional point cloud modeling in the establishment of an underground model. Summary of the Invention In view of the shortcomings in the prior art, the present disclosure provides a digital mine modeling method based on a three-dimensional point cloud database. By establishing a three-dimensional point cloud database for a mine and using a three-dimensional modeling software, a digital model for mining inversion can be established to reflect the rock relief and extension, mining areas and roadway locations, and geological structures. To attain the objective described above, the present disclosure provides a digital mine modeling method based on a three-dimensional point cloud database, including: collecting mine data, extracting effective coordinate point position information of rock interfaces according to criteria for extracting point cloud data, establishing a three-dimensional point cloud database of the mine; establishing a model of rock interfaces, mining areas and roadways by using the three-dimensional point cloud data of the mine, so as to establish a three-dimensional geological model of the ground surface and establish a three-dimensional geological model of the mine; performing grid division on the three-dimensional geological model of the mine, to obtain a digital model for mining inversion that reproduces each identified rock strata relief, mining areas and roadway locations, fault morphology and surface morphology. Thus, the fitness and accuracy of the model are improved. The method includes the following steps: S1. Determining dimensions of the model by collecting geological data of a target mining area, and analyzing and obtaining spatial coordinate sets of upper and lower interfaces of each rock stratum, a spatial coordinate set of inflection points of mining sections, a spatial coordinate set of inflection points of key roadways, a spatial coordinate set of contour lines of the ground surface, fault locations, fault displacements and fault angles; S2. Processing the spatial coordinate sets of the upper and lower interfaces of each rock stratum, extracting effective coordinate point position information, and establishing a three-dimensional point cloud database of the mine; $2.1. Determining the model accuracy € (0 < € < 1) according to simulation requirements; $2.2. Determining a criteria for extracting the three-dimensional point cloud data of the mine and a minimum significant coefficient for the extracted spatial coordinate sets of the rock interfaces; supposing there are m spatial coordinate sets of rock interfaces and » criteria for extracting the three-dimensional point cloud data of the mine, calculating a significant coefficient Px for the &” spatial coordinate set of rock interfaces with the following formula: Pe=YwX+l (sk<m 1€ign) Jul Sw=l, Xe(010), 1e{0.1) F=} where X; is the significant coefficient of a criterion 7 in the » criteria for extracting the three-dimensional point cloud data of the mine, w; is a weight of the criterion 7, and / is a correction coefficient; determining a minimum significant coefficient Pi» for the extracted spatial coordinate sets of the Poi : min rock interfaces, which meets # sx , according to the simulation requirements, where Pax is a maximum significant coefficient; $2.3. Directly extracting the spatial coordinate sets of the rock interfaces that meet Px > Pin into the three-dimensional point cloud database of the mine and designating names to them respectively; S2.4. Categorizing spatial coordinate sets for remaining rock interfaces into two categories, i.e., AN Pr ee IE HE en CL, HE coordinate sets that meet & wes or coordinate sets that meet £ wx , where a is an Pr mmnneeneonn. Sn, fl ES adjustable coefficient and a € (0, 1); deleting rock strata that meet £ wax Lie, combining the Ps: a — EF rock strata with neighboring rock strata that meet £ == to form new coordinate sets for rock interfaces, and extracting the coordinate sets into the three-dimensional point cloud database of the mine and designating names to them respectively; $3. Establishing an underground part of a three-dimensional geological model that reproduces fault morphology, mining areas and roadways; S3.1. Importing the spatial coordinate sets of upper and low interfaces of the identified rock strata into Rhino software, and fitting the rock interfaces hierarchically from bottom to top with the Rhino software; establishing a cuboid model according to the dimensions of the model, cutting the rock strata according to the rock interfaces obtained through the fitting and designating names to the rock strata respectively, so as to obtain a three-dimensional underground geological model in a cuboid structure; $3.2. Establishing fault boundaries according to the fault locations, cutting the coal seam according to the fault angles; and adjusting the cut coal seam according to the fault displacements to reproduce fault morphology; S3.3. Inputting the spatial coordinate set of the inflection points of mining sections and the spatial coordinate set of the inflection points of key roadways into the Rhino software, and fitting mining sections and key roadways at the working faces and designating names to them respectively; S3.4. Performing Boolean operations on the spatial model of the rock strata, the mining sections and the key roadways, removing a redundant part in the spatial model of the rock strata, so as to form the underground part of the three-dimensional geological model of the entire mining area; S4. Establishing a three-dimensional geological model of the ground surface according to the spatial coordinate set of the contour lines of the ground surface, designating a name to it, and combining the two models together, so as to establish the three-dimensional geological model of the entire mining area; S35. Performing grid division on the three-dimensional geological model of the entire mining area by using the grid division tool that comes with the Rhino software, setting a grid dimension according to actual requirements, performing grid division on the model again by using a plug-in Griddle, and modifying grids by limiting a grid shape and maximum and minimum edge dimensions, so as to make the grid change uniform; exporting the model as a grid file that can be recognized by FLAC 3D software, and importing the grid file into the simulation software, so as to obtain a digital model for mining inversion of the entire mine. The geological data collected in the step S1 includes a mining engineering plan, an engineering geological profile, a geological borehole histogram and a contour map of the coal seam floor of the target mining area. The spatial coordinate sets of the upper and lower interfaces of the rock strata, the spatial coordinate set of the inflection points of mining sections, the spatial coordinate sets of the inflection points of key roadways, the spatial coordinate set of the contour lines of the ground surface, the dimensions of the model, the fault locations, the fault angles and the fault displacements are obtained by analyzing the geological data of the target mining area. The key roadways mainly refer to the return-air roadways and haulage roadways of the sections. The spatial coordinate sets of the upper and lower interfaces of the rock strata in the step S1 are obtained by analyzing the geological borehole histogram. Because of fixed latitude distance and longitude distance of each borehole, the horizontal coordinate X and the vertical coordinate Y in the coordinate set of the boreholes can be directly obtained. According to the number p of rock strata and a buried depth 4, the coordinate set of the boreholes can be obtained as (X, Y, hu), (X,Y, A2), ..., X,Y, hy), ..., (X,Y, hp-1) from top to bottom, where / <hp.,. All obtained coordinates are classified according to the geological unit "section" first, and further classified according to different rock strata in each "section", so that spatial coordinate sets of the upper and lower interfaces of the rock strata are obtained. The model accuracy € in the step S2 is the precision of the simulation model reproduction determined according to the simulation requirements, which can be the precision of some local continuous strata reproduction or the precision of all the strata reproduction. The richer the geological data and the denser the geological boreholes are, the higher the maximum reproduction accuracy of the model is. Because it is impossible to achieve 100% reproduction of the underground appearance of the mine, € is smaller than 1. The criteria for extracting three-dimensional point cloud data of the mine in the step S2 are detailed requirements determined according to the simulation requirements, the underground rock morphology of the mine, the geological borehole histogram, and the spatial coordinate sets of rock interfaces. The precision of the detailed requirements directly influences the establishment of the three-dimensional point cloud database of the mine. The identified rock strata in the step S3 refer to the rock strata identified in the three-dimensional point cloud database of the mine. The higher the model accuracy is and the more detailed the criteria for extracting the three-dimensional point cloud data of the mine is, the more the number of identified rock strata is. The identified rock strata gradually formed chronologically in the long geological period can be used to distinguish stratigraphic information and geological age in the process of geological investigation, borehole exploration and mining. In the step S3, since the mining sections and the key roadways spatially overlap with the spatial model of the rock strata respectively, the difference operation in the Boolean operation function of the software is used to perform difference operation on the spatial model of the rock strata, the mining sections and the key roadways, so as to cut out the mining sections and the key roadways in the spatial model of the rock strata and avoid overlapping of spatial locations . Thus, the spatial model of the rock strata, the mining sections and the key roadways constitute the underground part of the three-dimensional geological model of the entire mining area. Beneficial effects: 1) in the method, an idea of three-dimensional point cloud modeling is applied to strata, and a three-dimensional point cloud database of the mine can be obtained according to the method; 2) The method can not only reproduce the actual geological structure of the coal seam, the relief and extension of the rock strata and the surface morphology to a maximum extent, but also establish digital mining inversion models with different precision according to the simulation requirements. Using this method to model and perform the ground stress field inversion can provide more comprehensive and scientific guidance for field practice such as mine development, roadway layout and support, etc. Brief Description of Drawings Fig. 11s a flow chart of the digital mine modeling method based on a three-dimensional point cloud database in the present invention; Fig. 2 is a flow chart of establishing a point cloud database in an embodiment of the present invention; Fig. 3 is a schematic diagram of a equal-proportion three-dimensional geological model in an embodiment of the present invention. Embodiments Some embodiments of the present invention will be further described below with reference to the accompanying drawings. Embodiment 1 As shown in Fig. 1, a digital mine modeling method based on a three-dimensional point cloud database, including the following steps: S1. Collecting data; Collecting a mining engineering plan, an engineering geological profile map, a geological borehole histogram, and a contour map of the coal seam floor of the target mining area, determining the dimensions of the model, and analyzing and obtaining spatial coordinate sets of the upper and lower rock interfaces, a spatial coordinate set of the inflection points of mining sections, a spatial coordinate set of the inflection points of key roadways, fault locations, fault displacements and fault angles. The spatial coordinate sets of the upper and lower interfaces of the rock strata are obtained by analyzing the geological borehole histogram. Because of fixed latitude distance and longitude distance of each borehole, the horizontal coordinate X and the vertical coordinate Y in the coordinate set of the boreholes can be directly obtained. According to the number p of rock strata and a buried depth 4, the coordinate set of the boreholes can be obtained as (X, Y, 41), (X, Y, 42), ..., X,Y, hy), ..., (X,Y, hp-1) from top to bottom, where A;<hp+1. All obtained coordinates are classified according to the geological unit "section" first, and further classified according to different rock stratum in each "section", so that spatial coordinate sets of the upper and lower interfaces of the rock strata are obtained. S2. Establishing a three-dimensional point could database of the mine; 1) determining the model accuracy ¢ (¢ < 1) according to the simulation requirements; 2) determining a criteria for extracting the three-dimensional point cloud data of the mine and a minimum significant coefficient for the extracted spatial coordinate sets of the rock interfaces; Supposing there are m spatial coordinate sets of rock interfaces and n criteria for extracting the three-dimensional point cloud data of the mine, Px is a significant coefficient of the A” spatial coordinate set of the rock interfaces, JX; is a significant coefficient of a criterion 7 in the » criteria for extracting the three-dimensional point cloud data of the mine, w; is the weight of the criterion 7, and 1 is a correction coefficient, then: Pe=3"wXi+l (ISk<m 18i<n) Sw=l, X=(010), 7e(0.]) where determining a minimum significant coefficient Pui» for the extracted spatial coordinate sets of the P wie | Rr rock interfaces, which meets # 2x, according to the simulation requirements. 3) Directly extracting the spatial coordinate sets of the rock interfaces that meet Py > Pin into the three-dimensional point cloud database of the mine and designate names to them respectively; 4) Categorize the spatial coordinate sets for remaining rock interfaces into two categories, i.e., Pe Ps J. SY i MBq ren 533 31 ee SRE coordinate sets that meet # wx or coordinate sets that meet F max . where a is an Ps — EE adjustable coefficient and a € (0, 1), deleting rock strata that meet sex ,1.e., combining the Pi LR ny e—— 0 FLY rock strata with neighboring rock strata that meet Fw==x to obtain new coordinate sets for rock interfaces, and extracting the coordinate sets into the three-dimensional point cloud database of the mine and designate names to them respectively. $3. Constructing an underground part of a three-dimensional geological model that reproduces fault morphology, mining areas and roadways; 1) Importing the spatial coordinate sets of upper and low interfaces of the identified rock strata into Rhino software, and fit the rock interfaces hierarchically from bottom to top with the Rhino software; establish a cuboid model according to the dimensions of the model, cut the rock strata according to the rock interfaces obtained through the fitting and designate names to the rock strata respectively, so as to obtain an three-dimensional underground geological model; 2) Establishing fault boundaries according to the fault locations, and cutting the coal seam according to the fault angles; and adjusting the cut coal seam according to the fault displacements to reproduce fault morphology; 3) Inputting the spatial coordinate set of the inflection points of mining sections and the spatial coordinate set of the inflection points of key roadways into the Rhino software, and fit the mining sections and key roadways at working faces and designate names to them respectively; 4) Performing Boolean operations on the spatial model of the rock strata, the mining sections and the key roadways, removing a redundant part in the spatial model of the rock strata, so as to form the underground part of the three-dimensional geological model of the entire mining area; S4. Establishing a three-dimensional geological model of the mine; establishing a three-dimensional geological model of the ground surface according to the spatial coordinate set of the contour lines of the ground surface, designating a name to it, and combining the two models together, so as to establish the three-dimensional geological model of the entire mining area; SS. Performing grid division on the model to obtain an digital model for mining inversion; Performing grid division on the three-dimensional geological model of the entire mining area by using the grid division tool that comes with the Rhino software, where the grid dimension is set according to actual requirements, performing grid division on the model again by using a plug-in Griddle, and modifying grids by limiting a grid shape and maximum and minimum edge dimensions, so as to make the grid change uniform; exporting the model as a grid file that can be recognized by a simulation software, and importing the grid file into the simulation software, so as to obtain a digital model for mining inversion of the entire mine. The present invention provides a digital mine modeling method based on a three-dimensional point cloud database. The method provided by the present disclosure will be explained in detail by taking a mine as an example. S1. Collecting data; collecting a mining engineering plan, an engineering geological profile map, a geological borehole histogram, and a contour map of the coal seam floor of the target mining area, determining the dimensions of the model, and analyzing and obtaining spatial coordinate sets of the upper and lower rock interfaces, a spatial coordinate set of the inflection points of mining sections, a spatial coordinate sets of the inflection points of key roadways, fault locations, fault displacements and fault angles. The spatial coordinate sets of the upper and lower interfaces of the rock strata are obtained by analyzing the geological borehole histogram. Because of fixed latitude distance and longitude distance of each borehole, the horizontal coordinate X and the vertical coordinate Y in the coordinate set of the boreholes can be directly obtained. According to the number p of rock strata and a buried depth / ;, the coordinate set of the boreholes can be obtained as (X, Y, Au), (X, Y, ha), XY, hy), ..., (X,Y, hp-1) from top to bottom, where #<hp-;. All obtained coordinates are classified according to the geological unit "section" first, and further classified according to different rock stratum in each "section", so as to obtain m (m>p+1) spatial coordinate sets of the upper and lower interfaces of the rock strata. S2. Establishing a three-dimensional point cloud database of the mine according to the flow in Fig.2 which is a flow chart for establishing a three-dimensional point cloud database according to an embodiment; 1) Determining the model accuracy € as 0.4, an adjustable coefficient a as 0.8 according to the simulation requirements; 2) Determining the criteria for extracting the three-dimensional point cloud data of the mine and a minimum significant coefficient for the extracted spatial coordinate sets of the rock interfaces; Suppose there are 4 criteria for extracting the three-dimensional point cloud data of the mine, Pris a significant coefficient of the 4” spatial coordinate set of the rock interfaces, and JX; is a significant coefficient of a criterion 7 in the 4 criteria for extracting the three-dimensional point cloud data of the mine, wherein the criterion 1 is a thickness of rock stratum, and w;=0.4; the criterion 2 is stratigraphic stability, and w>=0.3; the criterion 3 is distribution range, and ws=0.2; the criterion 4 is borehole density, and w= 0.1; the specific values of X}, X;, Xz and Xs are shown in Table 1, Table 2, Table 3 and Table 4. The correction coefficient / is 0.5. 4 Pe=3 willis! =04X (+ 03X24 02X:+ 015 « +03 (Isk<s4) ani A minimum significant coefficient Pm» for the extracted three-dimensional coordinate sets of the rock interfaces is determined as 7 according to the simulation requirements, and Prax is 10.5, which Pin os meets £ max Pe Ps ns 3A noni ie, CCE FP owas Jie, P>3.360 Fax Lie, Pe<<3.36. Table 1. Thickness of Rock Stratum and Value of Xi Thickness of <0.5m 0.5m to Sm 5m to 10m 10m to 20m > 20m rock stratum xX 2 5 7 10 Table 2. Stratigraphic Stability and Value of X, Stratigraphic | Very unstable | ratigraphic | Very unstable | Less Unstable | Normal More Stable Very stable stability Xs 2 5 7 10 Table 3. Distribution Range and Value of X3 Distribution <30% 30% -50% 50% - 70% 70% - 90% 290% range X3 3 5 7 10 Table 4. Borehole Density and Value of X4 Borehole Very Disperse | Disperse Normal More Dense Very dense density Xa 3 3 7 10 3) Directly extracting the spatial coordinate sets of the rock interfaces that meet Pe>7 into the three-dimensional point cloud database of the mine and designate names to them respectively; 4) Categorizing the spatial coordinate sets for remaining rock interfaces into two categories, i.e., coordinate sets that meet P«>3.36 or coordinate sets that meet Pr<<3.36, deleting rock strata that meet Px<<3.36, i.e, combining the rock strata with neighboring rock strata that meet Px>3.36 to obtain new coordinate sets for rock interfaces, and extracting the coordinate sets into the three-dimensional point cloud database of the mine and designate names to them respectively. S3. Establishing an underground part of a three-dimensional geological model that reproduces fault morphology, mining areas and roadways; 1) Importing the spatial coordinate sets of upper and low interfaces of the identified rock strata into Rhino software, and fitting the rock interfaces hierarchically from bottom to top with the Rhino software; establishing a cuboid model according to the dimensions of the model, cutting the rock strata according to the rock interfaces obtained through the fitting the rock strata and designate names to the rock strata respectively, so as to obtain a three-dimensional underground geological model; 2) Establishing fault boundaries according to the fault locations and cutting the coal seam according to the fault angles; and adjusting the cut coal seam according to the fault displacements to reproduce fault morphology; 3) Inputting the spatial coordinate set of the inflection points of mining sections and the spatial coordinate set of the inflection points of key roadways into the Rhino software, and fitting the mining sections and the key roadways at working faces and designating names to them respectively; 4) Performing Boolean operations on the spatial model of the rock strata, the mining sections and the key roadways, remove a redundant part in the spatial model of the rock strata, so as to form the underground part of the three-dimensional geological model of the entire mining area; S4. Establishing a three-dimensional geological model of the mine; Establishing a three-dimensional geological model of the ground surface according to the spatial coordinate set of the contour lines of the ground surface, designating a name to it, and combining the two models together, so as to establish the three-dimensional geological model of the entire mining area; S35. Perform grid division on the model to obtain an inversed digital model for mining in the mine; Preliminarily grid the three-dimensional geological model of the entire mining area by using the grid division tool that comes with the Rhino software, where the grid dimension is set according to the actual requirements, performing grid division on the model again by using a plug-in Griddle, and modifying grids by limiting a grid shape and maximum and minimum edge dimensions, so as to make grid change uniform; exporting the model as a grid file that can be recognized by a simulation software, and import the grid file into the simulation software, so as to obtain a digital model for mining inversion of the entire mine as shown in Fig. 3. The above embodiments are used only to illustrate the technical solution of the present invention, and not to limit the present invention. While the present invention is described in detail according to the above embodiments, those skilled in the art should appreciate that various modifications may be made to the technical solution described in the above embodiments or equivalent replacements may be made to some or all technical features in the technical solution; however, all such modifications or replacements me née of the technical sc depart from the scope defined by the claims.

Claims

Claims A digital mine modeling method based on a three-dimensional point cloud database, comprising: collecting mine data, extracting effective coordinate point position information of rock interfaces according to criteria for extracting point cloud data, establishing a three-dimensional point cloud database of a mine; establishing a model of rock interfaces, mining areas and roadways by using the three-dimensional point cloud data of the mine, establishing a three-dimensional geological model of ground surface, so as to establish a three-dimensional geological model of the mine, perform grid division on the three-dimensional geological model of the mine, so as to obtain a digital model for mining inversion that reproduces each identified rock strata relief, mining areas and roadway locations, fault morphology and surface morphology, and improve fitness and accuracy of the model; the method comprising: S1. determining dimensions of the model by collecting geological data of a target mining area, and analyzing and obtaining spatial coordinate sets of upper and lower interfaces of each rock stratum, a spatial coordinate set of inflection points of mining sections, a spatial coordinate set of inflection points of key roadways, a spatial coordinate set of contour lines of the ground surface, fault locations, fault displacements and fault angles; S2. processing the spatial coordinate sets of the upper and lower interfaces of each rock stratum, extracting effective coordinate point position information, and establishing a three-dimensional point cloud database of the mine; S2.

1. determining the accuracy (0 < € < 1) of the model according to simulation requirements; $2.

2. determining a criteria for extracting the three-dimensional point cloud data of the mine and a minimum significant coefficient for the extracted spatial coordinate sets of the rock interfaces; supposing there are m spatial coordinate sets of rock interfaces and » criteria for extracting the three-dimensional point cloud data of the mine, calculating a significant coefficient Py for a k” spatial coordinate set of rock interfaces with a following formula: Pe=Y w+! (Isksm. 1gism Swe=l, Xee(010), Je(01); imi wherein X; is a significant coefficient of a criterion 7 in the » criteria for extracting the three-dimensional point cloud data of the mine, w; is a weight of the criterion 7, and / is a correction coefficient; determining a minimum significant coefficient Pi, for the extracted spatial coordinate sets of Powis A ng the rock interfaces, which meet # max , according to the simulation requirements, wherein Prax is a maximum significant coefficient; $2.

3. directly extracting the spatial coordinate sets of the rock interfaces that meet Pr 2 Pin into the three-dimensional point cloud database of the mine and designating names to them respectively; S2.

4. categorizing spatial coordinate sets for remaining rock interfaces into two categories, i.e., Ps Pr ee EE XE eee AE coordinate sets that meet & ws or coordinate sets that meet # wax , wherein a is Px ee SU XE an adjustable coefficient and a €(0, 1), deleting the rock strata that meet # m= , le, Ps e335 GE combining the rock strata with neighboring rock strata that meet 2 ws to form new coordinate sets of rock interfaces, and extracting the coordinate sets into the three-dimensional point cloud database of the mine and designating names to them respectively; S3. establishing an underground part of a three-dimensional geological model that reproduces fault morphology, mining areas and roadways; $3.

1. importing spatial coordinate sets of upper and low interfaces of identified rock strata into Rhino software, and fitting the rock interfaces hierarchically from bottom to top with the Rhino software; establishing a cuboid model according to the dimensions of the model, cutting the rock strata according to the rock interfaces obtained through the fitting and designating names to the rock strata respectively, so as to obtain a three-dimensional underground geological model in a cuboid structure; S3.

2. establishing fault boundaries according to the fault locations, cutting coal seam according to the fault angles; and adjusting the cut coal seam according to the fault displacements to reproduce fault morphology; $3.

3. inputting the spatial coordinate set of the inflection points of the mining sections and the spatial coordinate set of the inflection points of the key roadways into the Rhino software, and fitting the mining sections and the key roadways at working faces and designating names to them respectively; S3.

4. performing Boolean operations on the spatial model of the rock strata, the mining sections and the key roadways, removing a redundant part in the spatial model of the rock strata, so as to form an underground part of a three-dimensional geological model of an entire mining area; S4. establishing a three-dimensional geological model of the ground surface according to the spatial coordinate set of the contour lines of the ground surface, designating a name to it, and combining the two models together, so as to establish the three-dimensional geological model of the entire mining area; S35. performing grid division on the three-dimensional geological model of the entire mining area by using a grid division tool that comes with the Rhino software, setting a grid dimension according to actual requirements, performing grid division on the model again by using a plug-in Griddle, and modifying grids by limiting a grid shape and maximum and minimum edge dimensions, so as to make grid change uniform; exporting the model as a grid file that can be recognized by FLAC 3D software, and importing the grid file into the simulation software, so as to obtain a digital model for mining inversion of the entire mine.

2. The digital mine modeling method based on a three-dimensional point cloud database according to claim 1, wherein the geological data collected in the step S1 comprises a mining engineering plan, an engineering geological profile , a geological borehole histogram and a contour map of coal seam floor of the target mining area; the spatial coordinate sets of the upper and lower interfaces of the rock strata, the spatial coordinate set of the inflection points of the mining sections, the spatial coordinate set of the inflection points of key roadways, the spatial coordinate set of the contour lines of the ground surface, the dimensions of the model, the fault locations, the fault angles and the fault displacements are obtained by analyzing the geological data of the target mining area, and the key roadways refer to return-air roadways and haulage roadways of sections.

3. The digital mine modeling method based on a three-dimensional point cloud database according to claim 1, wherein the spatial coordinate sets of the upper and lower interfaces of the rock strata in the step S1 are obtained by analyzing the geological borehole histogram; with fixed latitude distance and longitude distance of each borehole, horizontal coordinate X and vertical coordinate Y in a coordinate set of boreholes can be directly obtained; the coordinate set of the boreholes can be obtained as (X, Y, 41), (X,Y, 42), ..., (X,Y, Ay), ..., (X,Y, hp+1) from top to bottom according to a number p of rock strata and a buried depth A;,, wherein hA<h,-; all obtained coordinates are classified according to a geological unit "section", and classified again according to different rock strata in each "section", so that the spatial coordinate sets of the upper and lower interfaces of the rock strata are obtained.

4. The digital mine modeling method based on a three-dimensional point cloud database according to claim 1, wherein the accuracy € of the model in the step S2 is precision of simulation model reproduction determined according to the simulation requirements, which can be precision of some local continuous strata reproduction or precision of all the strata reproduction; the richer the geological data is and the denser the geological boreholes are, the higher the maximum accuracy of the model reproduction is; because it is impossible to achieve 100% reproduction of underground appearance of the mine, € is smaller than 1.

5. The digital mine modeling method based on a three-dimensional point cloud database according to claim 1, wherein the criteria for extracting three-dimensional point cloud data of the mine in the step S2 are detailed requirements determined according to the simulation requirements, the underground rock morphology of the mine, the geological borehole histogram, and the spatial coordinate sets of rock interfaces, and the precision of the detailed requirements directly influences the establishment of the three-dimensional point cloud database of the mine.

6. The digital mine modeling method based on a three-dimensional point cloud database according to claim 1, wherein the identified rock strata in the step S3 refer to rock strata identified in the three-dimensional point cloud database of the mine; the higher the accuracy of the model is and the more detailed the criteria for extracting the three-dimensional point cloud data of the mine is, the more a number of identified rock strata is; identified strata gradually formed chronologically in long geological period can be used to distinguish stratigraphic information and geological age in process of geological investigation, borehole exploration and mining.

7. The digital mine modeling method based on a three-dimensional point cloud database according to claim 1, wherein in the step S3, since the mining sections and the key roadways spatially overlap with the spatial model of the rock strata respectively, a difference operation in the Boolean operation function of software is used to perform the difference operation on the spatial model of the rock strata, the mining sections and the key roadways, so as to cut out the mining sections and the key roadways in the spatial model of the rock strata and avoid overlapping of spatial locations, so that the spatial model of the rock strata, the mining sections and the key roadways constitute the underground part of the three-dimensional geological model of the entire mining area.

Images