A method for automatically updating a profile of open-pit slope stability calculation

Abstract

The application provides a method for automatically updating a profile of open pit slope stability calculation, and relates to slope stability calculation. The method firstly establishes a three-dimensional model of the open pit slope and a lithology parameter database, and assigns lithology parameters to the three-dimensional ore rock model of the open pit slope; secondly, establishes the association between the three-dimensional model of the open pit slope and the profile of slope stability calculation, and establishes the association between the three-dimensional model of the open pit slope and the profile of slope stability calculation; then, automatically updates the slope step line of the profile of open pit slope stability calculation; finally, automatically updates the ore rock polygon of the profile of open pit slope stability calculation, and the automatic update of the ore rock polygon of the profile includes the update of the physical and mechanical parameters and the update of the geometric shape. The method breaks through the defects of the existing update method, automatically updates the profile of slope stability calculation after updating the three-dimensional model of the ore rock, and improves the efficiency and reliability of the slope stability calculation and analysis and evaluation.

Description

Technical Field

[0001] This invention relates to the field of slope stability calculation, and in particular to a method for automatically updating the profile diagram for open-pit mine slope stability calculation. Background Technology

[0002] Slope stability in open-pit mines is a crucial foundation for ensuring safe production. Slope stability calculation is a commonly used technique for assessing slope stability. Regular slope stability analysis and evaluation help identify hazards, detect landslide risks, classify and manage slopes, and propose reasonable slope prevention measures to ensure efficient, safe, and sustainable production. Slope stability calculation and analysis using slope profile diagrams is a common method for slope stability evaluation. As open-pit mining progresses, slope profile diagrams constantly change and require timely updates. However, manual updates are cumbersome, inefficient, and prone to errors. How to link the three-dimensional model of the open-pit mine's ore and rock with the slope stability calculation profile diagram, so that updates to the ore and rock three-dimensional model automatically update the slope stability calculation profile diagram, improving efficiency and accuracy of slope stability calculations, is an urgent problem that needs to be researched and solved.

[0003] Currently, the commonly used method for updating slope stability calculation profiles mainly involves manually redrawing the geometry of open-pit mines, rock strata, coal seams, and geological structures, and then assigning physical and mechanical parameters to the geometry of each coal and rock strata in the profile. Because this method requires entirely manual work, the updating efficiency of the calculation profile is low, the workload is high, and the number of calculated profiles obtained is often small, significantly impacting the reliability of slope stability analysis and evaluation results. To improve the updating efficiency of slope stability calculation profiles and the reliability of slope stability analysis and evaluation results, many scholars and engineers have conducted extensive research on methods for updating and drawing slope stability calculation profiles. Currently, commonly used slope stability calculation software in China, such as ANSYS (which calculates stability coefficients using the strength reduction method), Slide (which analyzes the stability of sliding surfaces using the vertical strip limit equilibrium analysis method), and FLAC3D (based on the finite difference method), all involve manually importing the location and profile into the calculation software, assigning lithological values, and then performing the calculation. If there are updates or changes in the rock strata, it is necessary to redraw the profile and reassign lithological values ​​before performing the calculation. In 1992, Ma Wentian et al. used cubic spline curves to fit the profile lines in geological cross-sections. In 2007, Zhu Ying et al. developed a cross-section drawing system based on Arc Engine. In 2009, Wang Jiming et al. applied knowledge reasoning mechanisms to the intelligent analysis of cross-section drawing algorithms. In 2013, Zhou Liangchen et al. simplified the modeling process, making cross-section drawing more stable and efficient. In 2015, Zhang Junqiang et al. used three-dimensional geological models and mineral resource point source databases as data sources for their cross-sections, and determined the drawing style of the cross-sections through style templates, achieving rapid drawing of standardized exploration cross-sections. Due to technological limitations, none of these studies involved automatic updating of slope stability calculation cross-sections. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to address the shortcomings of the prior art by providing an automatic updating method for the slope stability calculation profile of an open-pit mine. This method enables the automatic updating of the physical and mechanical parameters and geometric shapes of the slope step lines and ore polygons in the slope stability calculation profile, thereby improving the reliability of slope stability calculation and analysis evaluation.

[0005] To solve the above-mentioned technical problems, the technical solution adopted by the present invention is: an automatic updating method for open-pit mine slope stability calculation profile, comprising the following steps:

[0006] Step 1: Establish a 3D model and lithological parameter database for the open-pit mine slope, and assign lithological parameters to the 3D rock model of the open-pit mine slope. The specific method is as follows:

[0007] Step 1.1: Establish a 3D geometric model of the open-pit mine slope. The 3D model of the open-pit mine slope consists of a closed ore rock model and an open-pit mine stope model. The specific method is as follows:

[0008] Step 1.1.1: Generate a closed slope mineral mesh model;

[0009] The Kriging method is used to interpolate the ore-rock interface to generate triangular meshes of the ore-rock interface. The meshes of each interface are then combined from top to bottom to form a closed mesh.

[0010] Step 1.1.2: Generate a 3D mesh for the open-pit mine;

[0011] The open-pit mine's bench lines, topographic contour lines, geological feature lines, and other mine feature lines are subjected to Delaunay-constrained triangulation to generate a three-dimensional mesh of the open-pit mine.

[0012] Step 1.2: Establish a rock mass physical and mechanical parameter database and store the physical and mechanical parameters of the rock mass in the database; the rock mass physical and mechanical parameter database consists of lithological parameter tables and uses a lithological data structure to store and process lithological data;

[0013] The lithological data structure includes ID, stratum name, C, Phi, Gamma, and ore name.

[0014] Step 1.3: Assigning lithological parameters to the 3D ore and rock model. Assign lithological parameters to all closed-mesh models of the ore and rock in sequence. The assigned lithological data comes from the lithological database. The lithological data of the ore and rock is added to the closed-mesh model in the form of (ID, rock layer name, C, Phi, Gamma, ore name) and used as extended data for the closed-mesh model of the ore and rock.

[0015] Step 2: Establish the association between the 3D model of the open-pit mine slope and the slope stability calculation profile. This association is established by setting the names and coordinate positions of the graphic objects. The specific method is as follows:

[0016] Step 2.1: Associating Graphical Object Names. Assign appropriate names to the graphic elements in the 3D slope model and the graphic elements on the slope stability calculation profile to achieve name association between the graphic objects. The specific method is as follows:

[0017] Step 2.1.1: In the 3D model of the slope, set the profile name to the extended data of each profile line, and set the closed mesh model primitives of each mineral rock to the layer with the corresponding mineral rock name layer name;

[0018] Step 2.1.2: In the slope stability calculation profile, set the corresponding profile name as extended data for all rock layer polygons, mining benches, horizontal coordinate lines, vertical coordinate lines, and various types of graphic elements.

[0019] Step 2.2: Spatial Position Coordinate Association. During section view updates, the objects obtained from the sectioning are automatically updated to the section view using horizontal and elevation positioning coordinates. This requires setting positioning coordinates on the section view. The specific method is as follows:

[0020] Step 2.2.1: Setting the horizontal positioning coordinates for the cross-sectional view. Let the plane... Figure 1 Section lines The two endpoints and The horizontal coordinates are respectively and ,in For the cross-section line The first point, For the cross-section line The end point; on the section line The corresponding cross-sectional view has two horizontal positioning lines, both of which are vertical lines. Let the horizontal positioning line on the left be... The horizontal positioning line on the right is , the cross-section line The first point coordinates Set as left horizontal positioning line The linked data, used as one of the horizontal positioning coordinates, will also be used for the profile lines. End point coordinates Set as the right horizontal positioning line The linked data serves as another horizontal positioning coordinate;

[0021] Step 2.2.2: Setting the elevation coordinates of the profile view. Because the elevation coordinates of the profile view are not the actual elevation values, elevation positioning lines are needed to locate the elevation in the profile view. The elevation positioning lines are multiple horizontal straight lines, and each elevation positioning line represents an actual elevation value. The actual elevation values ​​are set as the link data of each elevation positioning line, which serves as the elevation positioning coordinates.

[0022] Step 3: Automatic Update of Slope Bend Lines in Open-Pit Mine Slope Stability Calculation Profile. Based on the generated intersection lines of the profile lines and the 3D mesh of the open-pit mine, the slope bench lines are updated in the profile diagram according to the correspondence between the profile lines and the profile diagram. The specific method is as follows:

[0023] Step 3.1: Generate the intersection line of the profile line and the 3D mesh of the mining area;

[0024] Finding the intersection line between a profile line and the 3D mesh of the mining area is essentially finding the intersection line between the vertical plane containing the profile line and the triangular mesh. The vertical plane containing the profile line can be represented as two triangular faces. Therefore, finding the intersection line between the profile line and the mining area triangular mesh is transformed into finding the intersection line between the triangular faces. Finally, connecting these intersection lines in sequence gives the intersection line between the profile line and the mining area triangular mesh. Let the vertex set of the obtained intersection line be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , , where n is the total number of vertices of the intersection line.

[0025] Step 3.2: Converting the intersection line of the slope steps to the cross-sectional view. The coordinates of the intersection line vertex are converted from plane rectangular coordinates to relative coordinates of the horizontal distance to the first point of the cross-sectional line. After calculation using the horizontal and elevation positioning coordinates on the cross-sectional view, the coordinates are converted back to those coordinates on the cross-sectional view, thus converting the intersection line of the slope steps to the cross-sectional view. The specific method is as follows:

[0026] Step 3.2.1: Convert the coordinates of the intersection vertex to relative coordinates. Let the converted vertex set be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , Where n is the total number of vertices of the intersection line, and the coordinate transformation is shown in the following formula:

[0027] (1)

[0028] in, The i-th vertex of the intersection line The original coordinates, For the section line section line The first point The coordinates.

[0029] Step 3.2.2: Convert the coordinates of the intersection vertex to the coordinates of the profile view. Let the converted vertex set be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , Where n is the total number of vertices of the intersection line, and the coordinate transformation is shown in the following formula:

[0030] (2)

[0031] in, The i-th vertex of the intersection line The original coordinates, This refers to the x-coordinate of the left horizontal positioning line on the corresponding cross-sectional view. Let y be any elevation positioning line on the corresponding profile. This is the actual elevation value of the elevation positioning line.

[0032] Step 3.2.3: Update the slope step lines to the profile view. Based on the converted slope step line vertex set... The coordinates in the graph are used to draw these vertices as slope step lines, and the date of the slope step survey is set as the extended data of the slope step lines in the profile.

[0033] Step 4: Automatic Update of Ore-Rock Polygons in Open-Pit Mine Slope Stability Profile. Automatic updating of ore-rock polygons in the profile includes updating the physical and mechanical parameters and geometric shape of the ore-rock polygons. The specific method is as follows:

[0034] Step 4.1: Update the physical and mechanical parameters of the ore-rock polygon;

[0035] After updating the physical and mechanical parameters of the closed ore and rock model in the 3D model of the open-pit mine slope, all rock layer polygon objects in all cross-sections are traversed. If the rock layer name of the cross-section rock layer polygon is the same as the rock layer name in the physical and mechanical parameters of the ore and rock model, it is then determined whether the physical and mechanical parameters of the cross-section rock layer polygon are consistent with the physical and mechanical parameters of the ore and rock model. If the two physical and mechanical parameters are inconsistent, the physical and mechanical parameters of the cross-section rock layer polygon are updated using the physical and mechanical parameters of the ore and rock model.

[0036] Step 4.2: Update the geometry of the profile ore-rock polygon;

[0037] The update of the geometry of the ore-rock polygon in the profile is essentially the same as the method for updating the slope step line of the profile. The difference is that updating the slope step line of the profile involves finding the intersection line between the profile line and the 3D mesh of the stope, while updating the geometry of the ore-rock polygon in the profile involves finding the intersection line between the profile line and the closed ore-rock model. Since the closed ore-rock model is also a triangular mesh, the method is exactly the same. The method is also to generate the intersection line between the profile line and the closed ore-rock model, and then convert the intersection line to the corresponding profile. See step 3 for the specific method. Finally, the new ore-rock polygon is used to update the ore-rock polygons with the same rock layer name in the profile.

[0038] The beneficial effects of adopting the above technical solution are as follows: This invention provides an automatic updating method for open-pit mine slope stability calculation profiles, which establishes a three-dimensional model and lithological parameter database for the open-pit mine slope, assigns lithological parameter values ​​to the three-dimensional ore-rock model of the open-pit mine slope, establishes the association between the three-dimensional model of the open-pit mine slope and the slope stability calculation profile, and realizes the automatic updating of the physical and mechanical parameters and geometric shapes of the slope step lines and ore-rock polygons in the slope stability calculation profile. This overcomes the shortcomings of existing slope stability calculation profile updating methods, enabling automatic updating of the slope stability calculation profile after the ore-rock three-dimensional model is updated, thus improving the efficiency and reliability of slope stability calculation and analysis evaluation. Attached Figure Description

[0039] Figure 1 A flowchart of an automatic updating method for calculating profile diagrams of open-pit mine slope stability is provided in an embodiment of the present invention;

[0040] Figure 2 The closed ore-rock interface triangular mesh model provided in this embodiment of the invention;

[0041] Figure 3 The embodiments of this invention provide a three-dimensional mineral rock model of a slope with cross-section lines;

[0042] Figure 4 C3 is a slope stability calculation profile diagram associated with a three-dimensional mineral rock model of the slope provided in an embodiment of the present invention;

[0043] Figure 5 This is a slope stability calculation profile diagram provided by an embodiment of the present invention, showing the generated slope step lines and lithological polygons.

[0044] Figure 6 A three-dimensional mesh for an open-pit mine provided in an embodiment of the present invention;

[0045] Figure 7(a) shows the slope stability calculation profile before the update;

[0046] Figure 7(b) shows the slope stability calculation profile automatically updated using the method of the present invention;

[0047] Figure 8 The result of automatic updating of the ore-rock polygon in the profile diagram for calculating the stability of open-pit mine slopes provided in this embodiment of the invention. Detailed Implementation

[0048] The specific embodiments of the present invention will be described in further detail below with reference to the accompanying drawings and examples. The following examples are for illustrative purposes only and are not intended to limit the scope of the invention.

[0049] This embodiment is based on geological exploration borehole data and physical and mechanical parameters of soil and rock mass from an open-pit mine in Xilinhot, Inner Mongolia. It utilizes an automatic updating method for open-pit mine slope stability calculation profiles to automatically update the slope step lines and the physical and mechanical parameters and geometric shapes of the ore-rock polygons in the slope stability calculation profile. In this embodiment, an automatic updating method for open-pit mine slope stability calculation profiles is provided, such as... Figure 1 As shown, it includes the following steps:

[0050] Step 1: Establish a 3D model and lithological parameter database for the open-pit mine slope, and assign lithological parameters to the 3D rock model of the open-pit mine slope. The specific method is as follows:

[0051] Step 1.1: Establish a 3D geometric model of the open-pit mine slope. The 3D model of the open-pit mine slope consists of a closed ore rock model and an open-pit mine stope model. The specific method is as follows:

[0052] Step 1.1.1: Generate a closed slope mineral mesh model;

[0053] The Kriging method is used to interpolate the ore-rock interface to generate triangular meshes of the ore-rock interface. The meshes of each interface are then combined from top to bottom to form a closed mesh.

[0054] Step 1.1.2: Generate a 3D mesh for the open-pit mine;

[0055] The open-pit mine's bench lines, topographic contour lines, geological feature lines, and other mine feature lines are subjected to Delaunay-constrained triangulation to generate a three-dimensional mesh of the open-pit mine.

[0056] Step 1.2: Establish a rock mass physical and mechanical parameter database and store the physical and mechanical parameters of the rock mass in the database; the rock mass physical and mechanical parameter database consists of lithological parameter tables and uses a lithological data structure to store and process lithological data;

[0057] The lithological data structure includes ID, stratum name, C, Phi, Gamma, and ore name.

[0058] Step 1.3: Assigning lithological parameters to the 3D ore and rock model. Assign lithological parameters to all closed-mesh models of the ore and rock in sequence. The assigned lithological data comes from the lithological database. The lithological data of the ore and rock is added to the closed-mesh model in the form of (ID, rock layer name, C, Phi, Gamma, ore name) and used as extended data for the closed-mesh model of the ore and rock.

[0059] This implementation example uses data from 126 geological exploration boreholes in an open-pit mine in Xilinhot, Inner Mongolia. Within the modeling area, there are six strata: topsoil, Quaternary, sandstone, upper mudstone, coal, and lower mudstone. The Kriging method is used to interpolate the interfaces of these six strata, generating triangular meshes of the ore-rock interfaces. These meshes are then combined sequentially from top to bottom in the following order: topsoil, Quaternary, sandstone, upper mudstone, coal, and lower mudstone, forming a closed mesh. Figure 2 As shown in the figure; there are 364 step lines and 640 topographic contour lines within the modeling area. Delaunay-constrained triangle partitioning is performed using these step lines and topographic contour lines as constraint edges to generate a three-dimensional mesh for the open-pit mine. A rock mass physical and mechanical parameter database is established using Microsoft Access 2016. The data consists of a lithological parameter table. The lithological data structure includes ID, stratum name, C, Phi, Gamma, and ore name. The stratum parameters for this area are shown in Table 1. The lithological parameters of the waste material, Quaternary, sandstone, mudstone, coal, and mudstone in the lithological parameter table are assigned to the closed mesh models of the topsoil, Quaternary, sandstone, mudstone, coal, and mudstone, respectively, and used as extended data for the ore-rock closed mesh.

[0060] Table 1 Lithological Parameters

[0061] ID Rock strata name Cohesion (C) Angle of internal friction (Phi) Gamma (bulk weight) Mine Name 1 Discarded materials 15.13 14.00 1.8 Wumei 2 Fourth Series 0 26.9 20.3 Wumei 3 sandstone 0 14.3 20.1 Wumei 4 mudstone 50 25 19.6 Wumei 5 coal 100 30 12.2 Wumei

[0062] Step 2: Establish the association between the 3D model of the open-pit mine slope and the slope stability calculation profile. This association is established by setting the names and coordinate positions of the graphic objects. The specific method is as follows:

[0063] Step 2.1: Associating Graphical Object Names. Assign appropriate names to the graphic elements in the 3D slope model and the graphic elements on the slope stability calculation profile to achieve name association between the graphic objects. The specific method is as follows:

[0064] Step 2.1.1: In the 3D model of the slope, set the profile name to the extended data of each profile line, and set the closed mesh model primitives of each mineral rock to the layer named with the corresponding mineral rock name;

[0065] Step 2.1.2: In the slope stability calculation profile, set the corresponding profile name as extended data for all rock layer polygons, mining benches, horizontal coordinate lines, vertical coordinate lines, and various types of graphic elements.

[0066] Step 2.2: Spatial Position Coordinate Association. When updating the section view, the objects obtained from the sectioning are automatically updated to the section view using horizontal and elevation positioning coordinates. This requires setting positioning coordinates on the section view. The specific method is as follows:

[0067] Step 2.2.1: Setting the horizontal positioning coordinates for the cross-sectional view. Let the plane... Figure 1 Section lines The two endpoints and The horizontal coordinates are respectively and ,in For the cross-section line The first point, For the cross-section line The end point; on the section line The corresponding cross-sectional view has two horizontal positioning lines, both of which are vertical lines. Let the horizontal positioning line on the left be... The horizontal positioning line on the right is , the cross-section line The first point coordinates Set as left horizontal positioning line The linked data, used as one of the horizontal positioning coordinates, will also be used for the profile lines. End point coordinates Set as the right horizontal positioning line The linked data serves as another horizontal positioning coordinate;

[0068] Step 2.2.2: Setting the elevation coordinates of the profile view. Because the elevation coordinates of the profile view are not the actual elevation values, elevation positioning lines are needed to locate the elevation in the profile view. The elevation positioning lines are multiple horizontal straight lines, and each elevation positioning line represents an actual elevation value. The actual elevation values ​​are set as the link data of each elevation positioning line, which serves as the elevation positioning coordinates.

[0069] In this implementation example, the ore-rock model of an open-pit mine slope in Xilinhot, Inner Mongolia, has five profile lines named C3, C4, C5, C6, and C7. The profile names of these five lines are set as extended data for each of the five profile lines. The specific settings for the profile names and extended data are shown in Table 2. The 3D slope model contains six closed-mesh ore-rock models, from top to bottom: topsoil, Quaternary, sandstone, upper mudstone, coal, and lower mudstone. Layers with the same names as these six closed-mesh ore-rock models are created in the drawing, and these mesh models are placed into their respective layers. The 3D ore-rock model of the slope with profile lines is shown below. Figure 3 As shown.

[0070] Table 2. Settings for Profile Names and Extended Data

[0071] Section number Section name Extended data 1 C3 C3 2 C4 C4 3 C5 C5 4 C6 C6 5 C7 C7

[0072] Configure the settings for sections C3, C4, C5, C6, and C7 respectively. Taking section C3 as an example, it has 6 rock strata polygons, 1 slope step line, 2 horizontal positioning lines, and 7 elevation positioning lines. Set the corresponding section names to the extended data of these graphic objects; the starting point of section line C3 is... , the end point ,Will Set the link data for the left horizontal positioning line at the coordinates (39413336.500, 4872390.500). Set the link data for the right-side horizontal positioning line at the coordinates (39414392.726, 4872390.500); set the link data for the seven elevation values ​​from 830m to 1130m, with a spacing of 50m, as seven elevation positioning lines from bottom to top, as follows. Figure 4 As shown in the figure. The slope stability calculation profiles for slopes C3 to C7 are as follows. Figure 5 As shown.

[0073] Step 3: Automatic Update of Slope Bend Lines in Open-Pit Mine Slope Stability Calculation Profile. Based on the generated intersection lines of the profile lines and the 3D mesh of the open-pit mine, the slope bench lines are updated in the profile diagram according to the correspondence between the profile lines and the profile diagram. The specific method is as follows:

[0074] Step 3.1: Generate the intersection line of the profile line and the 3D mesh of the mining area;

[0075] Finding the intersection line between a profile line and the 3D mesh of the mining area is essentially finding the intersection line between the vertical plane containing the profile line and the triangular mesh. The vertical plane containing the profile line can be represented as two triangular faces. Therefore, finding the intersection line between the profile line and the mining area triangular mesh is transformed into finding the intersection line between the triangular faces. Finally, connecting these intersection lines in sequence gives the intersection line between the profile line and the mining area triangular mesh. Let the vertex set of the obtained intersection line be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , , where n is the total number of vertices of the intersection line.

[0076] Step 3.2: Converting the intersection line of the slope steps to the cross-sectional view. The coordinates of the intersection line vertex are converted from plane rectangular coordinates to relative coordinates of the horizontal distance to the first point of the cross-sectional line. After calculation using the horizontal and elevation positioning coordinates on the cross-sectional view, the coordinates are converted back to those coordinates on the cross-sectional view, thus converting the intersection line of the slope steps to the cross-sectional view. The specific method is as follows:

[0077] Step 3.2.1: Convert the coordinates of the intersection vertex to relative coordinates. Let the converted vertex set be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , Where n is the total number of vertices of the intersection line, and the coordinate transformation is shown in the following formula:

[0078] (1)

[0079] in, The i-th vertex of the intersection line The original coordinates, For the section line section line The first point The coordinates.

[0080] Step 3.2.2: Convert the coordinates of the intersection vertex to the coordinates of the profile view. Let the converted vertex set be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , Where n is the total number of vertices of the intersection line, and the coordinate transformation is shown in the following formula:

[0081] (2)

[0082] in, The i-th vertex of the intersection line The original coordinates, This refers to the x-coordinate of the left horizontal positioning line on the corresponding cross-sectional view. Let y be any elevation positioning line on the corresponding profile. This is the actual elevation value of the elevation positioning line.

[0083] Step 3.2.3: Update the slope step lines to the profile view. Based on the converted slope step line vertex set... The coordinates in the graph are used to draw these vertices as slope step lines, and the date of the slope step survey is set as the extended data of the slope step lines in the profile.

[0084] In this implementation example, the intersection lines of five profile lines in the 3D model of the ore and rock with the original 3D mesh of the open-pit mine are calculated. (See the 3D mesh of the open-pit mine.) Figure 6The x-coordinate of the slope step line vertex is converted to a relative coordinate of its horizontal distance to the first point of the profile line, and the y-coordinate is converted to the vertex elevation value. Using a coordinate transformation formula, these coordinates are then converted to the coordinates on the profile view, thus realizing the transformation and addition of the slope step intersection line to the profile view. Simultaneously, the slope step line uses the current date as its extended data. Following the same method, the intersection lines of the five profile lines in the 3D ore model with the 3D mesh of the new open-pit mine are calculated, and steps 3.1 and 3.2 above are repeated. Finally, the calculation results of the slope step lines from both methods can be displayed simultaneously on the profile view, achieving rapid and automatic updates of the slope step lines. When multiple slope step lines exist on the same profile view, they can be distinguished by the step date displayed in the extended data of each step line.

[0085] Taking section line C3 as an example, the first point of the section line is , the end point , The coordinates are (39413336.500, 4872390.500). The coordinates are (39414392.726,4872390.500), via the profile line Find the intersection line between the profile line and the original mining area triangular mesh, and convert the x-coordinate of the vertex of the slope step line to the starting point of the profile line. The relative coordinates of the horizontal distance point are converted from the y-coordinate to the elevation value of the step line vertex. The horizontal positioning line on the left side of the section view can be obtained from the section name "C3". With the horizontal positioning line on the right Any elevation positioning line , The horizontal positioning line on the left side of the cross-sectional view corresponding to section line C3. x-coordinate The value is 39413990.708. For any elevation positioning line on the cross-section diagram corresponding to section line C3 y-coordinate, The value is 4872649.970. This is the actual elevation of the elevation positioning line. The value is 1070m. Based on the above data, the final coordinates of the vertex of the transformed slope step line on the profile are: The vertex set data is drawn as a line, and the step date is 2022-10-03. This is set as the extended data of the slope step line to update the original step line with the step date of 2022-09-18. The result is shown in Figure 7(b).

[0086] Step 4: Automatic Update of Ore-Rock Polygons in Open-Pit Mine Slope Stability Profile. Automatic updating of ore-rock polygons in the profile includes updating the physical and mechanical parameters and geometric shape of the ore-rock polygons. The specific method is as follows:

[0087] Step 4.1: Update the physical and mechanical parameters of the ore-rock polygon;

[0088] After updating the physical and mechanical parameters of the closed ore and rock model in the 3D model of the open-pit mine slope, all rock layer polygon objects in all cross-sections are traversed. If the rock layer name of the cross-section rock layer polygon is the same as the rock layer name in the physical and mechanical parameters of the ore and rock model, it is then determined whether the physical and mechanical parameters of the cross-section rock layer polygon are consistent with the physical and mechanical parameters of the ore and rock model. If the two physical and mechanical parameters are inconsistent, the physical and mechanical parameters of the cross-section rock layer polygon are updated using the physical and mechanical parameters of the ore and rock model.

[0089] Step 4.2: Update the geometry of the profile ore-rock polygon;

[0090] The update of the geometry of the ore-rock polygon in the profile is essentially the same as the method for updating the slope step line of the profile. The difference is that updating the slope step line of the profile involves finding the intersection line between the profile line and the 3D mesh of the stope, while updating the geometry of the ore-rock polygon in the profile involves finding the intersection line between the profile line and the closed ore-rock model. Since the closed ore-rock model is also a triangular mesh, the method is exactly the same. The method is also to generate the intersection line between the profile line and the closed ore-rock model, and then convert the intersection line to the corresponding profile. See step 3 for the specific method. Finally, the new ore-rock polygon is used to update the ore-rock polygons with the same rock layer name in the profile.

[0091] In this implementation example, the physical and mechanical parameters of mudstone in the 3D model are updated from (0, 14.3, 20.1) to (50.0, 25.0, 19.6). This requires updating the physical and mechanical parameters of mudstone polygons in all cross-sections. Taking cross-section C3 as an example, by traversing the rock layer names of all rock layer polygons in this cross-section, the rock layer polygon named mudstone is obtained. The lithology parameters of this polygon are then reset to (50.0, 25.0, 19.6), thus updating the physical and mechanical parameters of the ore-rock polygon. The update of the geometry of the ore-rock polygon in the cross-section is achieved by finding the intersection line between the cross-section line and the 3D mesh of the stope, and then converting this intersection line onto the corresponding cross-section, thus updating the ore-rock polygons with the same rock layer name in the cross-section. The comparison results before and after the automatic update of the ore-rock polygon in the open-pit mine slope stability cross-section are shown below. Figure 8 As shown.

[0092] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit them. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope defined by the claims of the present invention.

Claims

1. A method for automatically updating a profile diagram for calculating the stability of an open-pit mine slope, characterized in that: Includes the following steps: Step 1: By establishing a three-dimensional model of the open-pit mine slope and a lithological parameter database, assign values ​​to the lithological parameters of the three-dimensional rock model of the open-pit mine slope; Step 2: Establish the association between the 3D model of the open-pit mine slope and the slope stability calculation profile. This is done by setting the names and coordinate positions of the graphic objects. Step 3: Automatic update of slope step lines in open-pit mine slope stability calculation profile. Based on the intersection of the profile line and the three-dimensional mesh of the open-pit mine, the slope step lines are updated in the profile according to the correspondence between the profile line and the profile. Step 4: Automatic update of the ore-rock polygon in the profile of open-pit mine slope stability. Automatic update of the ore-rock polygon in the profile includes updating the physical and mechanical parameters and the geometric shape of the ore-rock polygon. The specific method for step 2 is as follows: Step 2.1: Associating Graphical Object Names. Set corresponding names for the graphic elements in the 3D slope model and the slope stability calculation profile, thus associating the names of the graphic objects. In the 3D slope model, set the profile names as extended data for each profile line, and set the closed mesh model elements of each ore and rock to layers named with the corresponding ore and rock layer names. In the slope stability calculation profile, set the corresponding profile names as extended data for all rock layer polygons, stope benches, horizontal coordinate lines, vertical coordinate lines, and various types of labeled graphic elements. Step 2.2: Spatial location coordinate association. During the section view update, the objects obtained from the sectioning are automatically updated to the section view using horizontal and elevation positioning coordinates; assuming a section line on the plan view. The two endpoints and The horizontal coordinates are respectively and ,in For the cross-section line The first point, For the cross-section line The end point; on the section line The corresponding cross-sectional view has two horizontal positioning lines, both of which are vertical lines. Let the horizontal positioning line on the left be... The horizontal positioning line on the right is , the cross-section line The first point coordinates Set as left horizontal positioning line The linked data, used as one of the horizontal positioning coordinates, will also be used for the profile lines. End point coordinates Set as the right horizontal positioning line The linked data is used as another horizontal positioning coordinate; because the graphic elevation coordinates of the profile are not the actual elevation values, elevation positioning lines are needed to perform elevation positioning in the profile. The elevation positioning lines are multiple horizontal straight lines, and each elevation positioning line represents an actual elevation value. The actual elevation value is set as the linked data of each elevation positioning line as the elevation positioning coordinate. The specific method for step 3 is as follows: Step 3.1: Generate the intersection line between the profile line and the 3D mesh of the mining area. Finding the intersection line between the profile line and the 3D mesh of the mining area is essentially finding the intersection line between the vertical plane containing the profile line and the triangular mesh. The vertical plane containing the profile line can be represented as two triangular faces. Therefore, finding the intersection line between the profile line and the triangular mesh of the mining area is transformed into finding the intersection line between the triangular faces. Finally, connecting these intersection lines in sequence gives the intersection line between the profile line and the triangular mesh of the mining area. Let the vertex set of the obtained intersection line be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , n is the total number of vertices of the intersection line; Step 3.2: Transform the intersection line of the slope steps into a cross-sectional view. The coordinates of the intersection line vertices are converted from Cartesian coordinates to relative coordinates of the horizontal distance to the first point of the cross-sectional line. After calculation using the horizontal and elevation coordinates on the cross-sectional view, the coordinates are then converted back to those coordinates on the cross-sectional view. This completes the transformation of the slope step intersection line into a cross-sectional view. First, the coordinates of the intersection line vertices are converted to relative coordinates. Let the set of vertices after the transformation be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , Where n is the total number of vertices of the intersection line, and the coordinate transformation is shown in the following formula: （1） in, The i-th vertex of the intersection line The original coordinates, For the cross-section line The first point First, determine the coordinates of the intersection vertices; then, convert the coordinates of the intersection vertices to the coordinates of the cross section; let the converted vertex set be... ,in, Let be the i-th vertex of the intersection line, and let the coordinates of the vertex be . , Where n is the total number of vertices of the intersection line, and the coordinate transformation is shown in the following formula: （2） in, The i-th vertex of the intersection line The original coordinates, This refers to the x-coordinate of the left horizontal positioning line on the corresponding cross-sectional view. Let y be any elevation positioning line on the corresponding profile. The actual elevation value of the elevation positioning line is used; then the slope step line is updated to the profile view, according to the converted slope step line vertex set. The coordinates in the graph are used to draw these vertices as slope step lines, and the date of the slope step survey is set as the extended data of the slope step lines in the profile.

2. The method for automatically updating the profile diagram for calculating the stability of an open-pit mine slope according to claim 1, characterized in that: The specific method for step 1 is as follows: Step 1.1: Establish a 3D geometric model of the open-pit mine slope. The 3D model of the open-pit mine slope consists of a closed ore-rock model and an open-pit mine mining area model. First, the Kriging method is used to interpolate the ore-rock interface to generate triangular meshes of the ore-rock interface. The meshes of each interface are then combined from top to bottom to form a closed mesh. Second, the bench lines, topographic contour lines, geological feature point lines, and other mining area feature point lines of the open-pit mine are subjected to Delaunay constrained triangulation to generate a 3D mesh of the open-pit mine mining area. Step 1.2: Establish a rock mass physical and mechanical parameter database and store the physical and mechanical parameters of the rock mass in the database; the rock mass physical and mechanical parameter database consists of a lithological parameter table, and the storage and processing of lithological data are realized through a lithological data structure; the lithological data structure includes ID, rock layer name, C, Phi, Gamma and ore name; Step 1.3: Assign lithological parameters to the three-dimensional ore and rock model. Assign lithological parameters to all closed mesh models of ore and rock in sequence. The lithological data for assignment comes from the lithological database. The lithological data of ore and rock is added to the closed mesh model in the form of (ID, rock layer name, C, Phi, Gamma, ore name) and used as the extended data of the closed mesh of ore and rock.

3. The method for automatically updating the profile diagram for calculating the stability of an open-pit mine slope according to claim 1, characterized in that: The specific method for step 4 is as follows: Step 4.1: Update the physical and mechanical parameters of the ore-rock polygon. After updating the physical and mechanical parameters of the closed ore-rock model in the 3D model of the open-pit mine slope, traverse all rock layer polygon objects in all cross-sections. If the rock layer name of the cross-section rock layer polygon is the same as the rock layer name in the physical and mechanical parameters of the ore-rock model, then determine whether the physical and mechanical parameters of the cross-section rock layer polygon are consistent with the physical and mechanical parameters of the ore-rock model. If the two physical and mechanical parameters are inconsistent, update the physical and mechanical parameters of the cross-section rock layer polygon using the physical and mechanical parameters of the ore-rock model. Step 4.2: Update the geometry of the cross-section ore-rock polygons. The update of the geometry of the cross-section ore-rock polygons is essentially the same as the method for updating the slope step lines of the cross-section diagram. The difference is that updating the slope step lines of the cross-section diagram involves finding the intersection line between the cross-section line and the 3D mesh of the stope, while updating the geometry of the cross-section ore-rock polygons involves finding the intersection line between the cross-section line and the closed ore-rock model. Since the closed ore-rock model is also a triangular mesh, the method is exactly the same. The method is also to generate the intersection line between the cross-section line and the closed ore-rock model, then convert the intersection line to the corresponding cross-section diagram, and finally use the new ore-rock polygons to update the ore-rock polygons with the same rock layer name in the cross-section diagram.

Images