Tunnel face rock mass structural surface identification method and device and storage medium

Abstract

The invention discloses a tunnel face rock mass structural surface recognition method and device and a storage medium, and relates to the technical field of rock mass recognition. The method specifically comprises the steps that the distance between a scanning center point and each triangular patch on a tunnel face is obtained; a distance condition is set, the distance between the scanning center point and each triangular patch on the tunnel face is screened, and the triangular patches meeting the distance condition are obtained; the distance between the scanning center point and each triangular patch on the tunnel face is screened, and the triangular patches meeting the distance condition are obtained; filling the triangular patch of the first structural surface; and automatically combining the second structural surfaces with the regional relationship, and identifying a third structural surface with a complete region. Automatic recognition and geometric information extraction of the structural plane are completed through cooperation of a computer algorithm and manual operation, subjective influences of people on occurrence measurement of the structural plane are eliminated, the data size of redundant structural planes is reduced, the processing efficiency is improved, and meanwhile the measurement accuracy of the structural plane is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of rock mass identification, and particularly to a method, device and storage medium for identifying the structural planes of the rock mass at the tunnel face. Background Art

[0002] As the main body of the rock mass structure, the structural plane is an important part of the rock mass and also the main factor controlling the strength and stability of the rock mass. Its information is the basis for the division of the rock mass structure and the analysis of the rock mass stability, and largely determines the stability of the tunnel surrounding rock. Therefore, the statistics and evaluation of the structural planes of the rock mass in the construction area are extremely important in geological engineering.

[0003] The structural plane is a general term for various geological interfaces that cut the rock mass. The structural plane is formed during a long geological history. After the rock mass is subjected to geological stress, the structural plane also has certain morphological characteristics. Traditional manual measurement and analysis are also based on some characteristics of the structural plane, and the characteristics of the structural plane are important factors affecting the strength and other properties of the structural plane.

[0004] In the prior art, for the identification of the structural plane, tools such as a geological compass and a measuring rope are used for on-site measurement, and parameters such as the attitude, trace length, spacing and aperture of the rock mass structural plane are measured, and the data obtained from the on-site measurement are recorded. These data include the azimuth, dip angle, extension length, etc. of the structural plane. However, in the above solution, to ensure a detailed investigation of the entire rock mass surface to ensure that no structural plane that may affect the rock mass stability is missed, it results in including the measurement of redundant structural planes, leading to low overall measurement efficiency. At the same time, it is difficult to achieve automated analysis under complex geological conditions, and the results are subjective. Summary of the Invention

[0005] The main object of the present invention is to provide a method, device and storage medium for identifying the structural planes of the rock mass at the tunnel face, aiming to reduce the measurement of redundant structural planes and improve the accuracy of the structural plane measurement at the same time.

[0006] To achieve the above object, the present invention proposes a method for identifying the structural planes of the rock mass at the tunnel face, and the specific steps include:

[0007] Step S1: Divide the rock mass at the tunnel face into several triangular patches, and obtain the distances between the scanning center point and each triangular patch on the tunnel face;

[0008] Step S2: Set distance conditions, and screen the distances between the scanning center point and each triangular patch on the tunnel face to obtain the triangular patches that meet the distance conditions;

[0009] Step S3: Conduct connectivity analysis and size screening on the triangular facets that meet the distance condition in sequence to form the first structural plane, and mark the triangular facets that meet the distance condition but fail the connectivity analysis as background triangular facets;

[0010] Step S4: Fill in the triangular facets of the first structural plane to form a second structural plane with regional connectivity;

[0011] Step S5: Automatically merge the second structural planes with regional relationships to form a third structural plane with complete regions.

[0012] Furthermore, the said Step S1 includes:

[0013] Step S11: Obtain the three vertices (x k , y k , z k ) of the triangular facet, where k = 1, 2, 3;

[0014] Step S12: Obtain the plane equation of the triangular facet through the three vertices of the triangular facet;

[0015] Step S13: Obtain the distance from the scanning center point (x0, y0, z0) to the triangular facet through the plane equation of the triangular facet.

[0016] Furthermore, in the said Step S3, the connectivity analysis includes: traversing the triangular facets that meet the distance condition. When the included angle between two adjacent triangular facets that meet the distance condition is less than the angle threshold, it indicates that the two adjacent triangular facets that meet the distance condition are connected to each other, and mark the connected and distance - condition - meeting triangular facets as an intermediate structural plane.

[0017] Furthermore, in Step S3, the size screening includes setting a quantity threshold, traversing the triangular facets of the intermediate structural plane. If the number of triangular facets is greater than or equal to the quantity threshold, mark that an intermediate structural plane is found; if the number of triangular facets of the intermediate structural plane is less than the quantity threshold, delete the current intermediate structural plane.

[0018] Furthermore, in the said Step S4, filling in the triangular facets of the first structural plane includes: selecting the hollowed - out triangular facets, and judging the three adjacent faces of the hollowed - out triangular facets. If at least two adjacent faces belong to the same structural plane, fill the hollowed - out triangular facet into the first structural plane to form the second structural plane, and loop to judge the three adjacent faces of the hollowed - out triangular facet until the hollowed - out triangular facet cannot be filled into the first structural plane, where the hollowed - out triangular facet refers to all triangular facets outside the first structural plane.

[0019] Furthermore, after Step S4, it also includes aggregating the boundaries of the second structural plane, specifically including:

[0020] Step S41: Sequentially obtain background triangular patches. If the adjacent triangular patches at two boundaries of a background triangular patch belong to a certain second structural plane, fill the current background triangular patch into this second structural plane.

[0021] Step S42: Analyze the background triangular patches in a loop according to Step S41 until there are no background triangular patches that meet the conditions.

[0022] Furthermore, it further includes Step S5. The automatic merging of multiple said second structural planes specifically includes:

[0023] Step S51: Set the small structural planes in each second structural plane as Si, and sequentially perform plane fitting on Si to obtain the fitting plane Pi.

[0024] Step S52: Sequentially analyze the spatial position relationship between the current structural plane and other structural planes to determine whether to merge.

[0025] Furthermore, in Step S52, the parameters for determining whether to merge include at least one of the included angle between the normal vectors of the structural planes, the minimum net distance between the fitting planes of the structural planes, the planar distance between the structural planes, and the optimal structural plane merging.

[0026] This application also discloses a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, it implements the above-mentioned tunnel face rock mass structural plane recognition method.

[0027] This application also discloses a storage medium, on which a rock mass structural plane recognition program is stored. When this program is executed by a processor, it implements the above-mentioned tunnel face rock mass structural plane recognition method.

[0028] Adopting the above technical solutions has the following advantages:

[0029] The present invention calculates the distance from the scanning center point to each triangular patch, sets distance conditions, obtains the triangular patches that meet the distance conditions, then sequentially performs connectivity analysis and scale screening on the triangular patches that meet the distance conditions, screens the first structural plane, and forms a complete third structural plane in the region through filling and automatic merging of the first structural plane. Through the cooperation of computer algorithms and manual operations, the automatic recognition of the structural plane and the extraction of geometric information are completed, reducing the redundant structural plane data volume, improving the processing efficiency, and at the same time forming a complete and continuous structural plane network system, eliminating the pseudo-structural planes caused by discrete measurement points, and improving the accuracy of structural plane measurement.

[0030] The present invention proposes to perform connectivity analysis and scale screening on triangular patches that meet the distance condition in sequence, and mark the combination of connected triangular patches that meet the distance condition as intermediate structural planes. Traverse the triangular patches of the intermediate structural planes, and mark an intermediate structural plane that meets the quantity threshold, and record this intermediate structural plane as the first structural plane. Connectivity analysis can merge adjacent triangular patches with similar characteristics into a complete structural plane, thereby improving the integrity and accuracy of the recognition result, and retaining large-scale structural planes that have an important impact on tunnel stability. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] The following will describe the present invention in detail with reference to specific embodiments and the accompanying drawings, where:

[0032] Figure 1 is a flowchart of the scale screening and connectivity analysis of the present invention;

[0033] Figure 2 is a flowchart of the structural plane filling algorithm of the present invention;

[0034] Figure 3 is a point-to-plane distance diagram of the present invention;

[0035] Figure 4 is a triangular patch connectivity diagram of the present invention;

[0036] Figure 5 is a diagram of the hollowing situation of the structural plane of the present invention;

[0037] Figure 6 is a diagram of the recognition and filling result of the structural plane of the present invention;

[0038] Figure 7 is an effect diagram of the aggregation of the boundaries of the structural plane of the present invention;

[0039] Figure 8 is a diagram of the convex and concave structural planes on the excavation surface of the present invention;

[0040] Figure 9 is a diagram of the structural plane with weak integrity of the present invention;

[0041] Figure 10 is a diagram of the included angle between the normal vectors of the structural planes on the heading face of the present invention;

[0042] Figure 11 is a schematic diagram of the distance between adjacent structural planes of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0043] In order to make the objectives, technical solutions and advantages of the present invention clearer, the following will describe the present invention in detail with reference to the accompanying drawings and embodiments. It should be understood that the following specific embodiments are only used to explain the present invention and do not constitute a limitation to the present invention.

[0044] A method for identifying rock mass structural planes at the tunnel face, the specific steps including:

[0045] Step S1: Divide the rock mass at the tunnel face into several triangular patches, that is, construct the spatial relationship of the structural planes using topological relationships, and obtain the distances from the scanning center point to each triangular patch on the tunnel face.

[0046] As Figure 3 It can be seen from the simulation that points E and F are two points at different positions in the plane ABCD, and the distances from the point G outside the plane to the plane where E and F are located are both GH. Taking the scanning center of the 3D laser scanner as the center point, it can be inferred that the distance values from the scanning center point to each triangular patch in a certain structural plane should also be approximately equal. Then, analyze the connectivity between the triangular patches with approximately equal distances to identify and find each structural plane.

[0047] Step S11: Obtain the three vertices (x k , y k , z k ) of the triangular patch, k = 1, 2, 3; calculate the plane formula of the spatial triangular patch. Each triangular patch corresponds to three vertices. According to the axiom 3 of solid geometry: Through three points not on a straight line, there is and only one plane. That is, three non-collinear points determine a unique plane, obtaining the basis for determining a plane.

[0048] Step S12: Obtain the plane equation of the triangular patch through the three vertices of the triangular patch; the plane equation of this application is Ax + By + Cz + D = 0, where (x, y, z) is any point in the plane, A, B, and C are the normal vectors of the plane respectively, and D is the constant term, which determines the position of the plane in space. Among them, the coefficients of x, y, and z are the coordinates of a normal vector of the plane, that is Substitute (x k , y k , z k ), k = 1, 2, 3 into the general form Ax + By + Cz + D = 0 to obtain three plane equations. After solving, the values of A, B, C, and D can be obtained. Three consecutive vertices of a counterclockwise convex polygon can be selected and the following system of simultaneous equations can be solved to find A / D, B / D, and C / D.

[0049]

[0050] It can be solved using Cramer's rule,

[0051]

[0052] Step S13: Through the scanning center point (x0, y0, z0) and the plane equation of the triangular patch, obtain the distance of the triangular patch and input it into the Delaunay triangular mesh;

[0053] Expanding the determinant, the expression for calculating the plane coefficients can be obtained:

[0054]

[0055] Immediately afterwards, calculate the distance d from the scanning center point (x0, y0, z0) to each triangular plane. The formula is:

[0056]

[0057] In addition, the distance formula from a point to a plane in vector form can also be used, which will not be elaborated here; each distance d corresponds to the order of the triangular patch in the triangular mesh.

[0058] Step S2: Set the distance condition, screen the distances from the scanning center point to each triangular patch on the heading face, and obtain the triangular patches that meet the distance condition; among them, use the distance from the scanning center point to the triangular patch to identify the structural plane similar plane features. If the triangular patches in a certain area form a structural plane, then the distance values from the scanning center point to each triangular patch in this area should be close. At the same time, due to local concavity and convexity, it is speculated that this distance should be an interval. Therefore, the distance condition is an interval threshold. If the distance of the triangular patch is within the interval of the distance condition, it is a triangular patch that meets the distance condition; for example, the distance condition can be set between 0 and 0.7d. During actual operation, the distance threshold can be adjusted according to the effect.

[0059] As Figure 1 shown, Step S3: Perform scale screening and connectivity analysis on the triangular patches that meet the distance condition to form the first structural plane, and mark the triangular patches that meet the distance condition but fail the connectivity analysis as background triangular patches. Among them, by setting different distance conditions, find the triangular patches that meet the distance condition corresponding to different structural planes, and perform corresponding screening on the scale of the triangular patches that meet the distance condition. Secondly, perform connectivity analysis. Then, through multiple datasets of triangular patches that have undergone scale screening and connectivity analysis, in the corresponding distance conditions, the first structural plane composed of the triangular patch datasets, and the background triangular patches are the triangular patches that cannot be recognized and form a structural plane or are discarded triangular patches.

[0060] If two triangular patches are connected, but the difference in their distance values from the scanning center point is large, then they basically cannot form a structural plane; if the distance values of two triangular patches from the scanning center point are both close, but they are not connected, or there are no triangular patches that meet the connection conditions between them, they cannot form a plane together. As Figure 4As shown, if all the black triangular patches meet the threshold conditions and the triangular patches in the upper part are connected, it can be seen that a structural plane can basically be formed. However, the three triangular patches in the lower part are not connected to the structural plane in the upper part and cannot jointly form a structural plane. At the same time, although the three triangular patches in the lower part are connected and seem to form a structural plane, based on the characteristic that a structural plane has a certain scale, the number of them is too small and will ultimately be discarded. For a structural plane, the specific scheme is as follows:

[0061] The connectivity analysis includes: if the included angle between two adjacent triangular patches that meet the distance condition is less than a given angle threshold, it indicates that the two adjacent triangular patches that meet the distance condition are connected to each other. The given angle threshold can be set according to different scenarios. The triangular patches that are connected and meet the distance condition are grouped and marked as the intermediate structural plane.

[0062] The scale screening includes setting a quantity threshold and traversing the triangular patches of the intermediate structural plane. If the number of triangular patches of the intermediate structural plane is greater than or equal to the quantity threshold, it is marked that an intermediate structural plane that meets the quantity threshold is found. This intermediate structural plane is recorded as the first structural plane, and the number of the current first structural plane is marked and the number of triangular patches in the first structural plane is obtained. If the number of triangular patches of the intermediate structural plane is less than the quantity threshold, the current intermediate structural plane is deleted, and the connected triangular patches are not regarded as a structural plane. The number of triangular patches in each first structural plane starts counting from 1. The results will be visually displayed, and it can be seen from the naked eye whether a plane is formed. Through multiple visual tests, the angle threshold and the quantity threshold are obtained. The recommended values of the angle threshold and the quantity threshold are 40° and 20 respectively, and they can also be adjusted according to the data situation.

[0063] Structural plane filling:

[0064] As an implementation manner of the present application, due to the influence of the excavation machinery scraping on the surface of the tunnel rock mass, there may be local unevenness or even large undulations on the structural plane, which can be said to be pitted. Due to the limitation of the quantity threshold, this will cause single or small-area hollowing in the identified structural plane, and filling is required to make the identified structural plane regionally connected.

[0065] As Figure 2 shown, obtain the first structural plane in the Delaunay triangular mesh and determine the three adjacent faces of the hollow triangular patch. If at least two adjacent faces belong to the same first structural plane, then fill the hollow triangular patch into the first structural plane; otherwise, discard the hollow triangular patch, and loop to determine the three adjacent faces of the hollow triangular patch until the hollow triangular patch cannot be filled into the first structural plane. The first structural plane with the hollow triangular patch filled is recorded as the second structural plane, where the hollow triangular patch refers to all triangular patches other than the first structural plane. Select the hollow triangular patch. Specifically, the situations of the hollow triangular patch are mainly divided into the following several types:

[0066] (1) There is exactly one hollow triangular patch in the first structural plane. This situation is the simplest. The three adjacent faces of the hollow triangular patch can be analyzed. If the three adjacent faces form the same first structural plane, then directly add this hollow triangular patch to the first structural plane of the adjacent faces.

[0067] (2) The first structural plane has local small - area hollows. This situation is similar to the first one. Analyze each hollow triangular patch one by one. If two of its three adjacent faces form the same first structural plane, then this hollow triangular patch is added to the first structural plane formed by its adjacent faces. Analyze the hollow triangles in a loop until no more can be filled into the first structural plane. As shown in the left - hand side figure of the appendix, the hollow part can be filled in the order of 1, 5, 6, 2, 4, 3. Figure 5 The left - hand side figure of the appendix can fill the hollow part in the order of 1, 5, 6, 2, 4, 3.

[0068] Analyze and fill the adjacent structural planes according to the above two situations, and finally form the second structural plane that fills the hollow triangular patches. Finally, update and save the second - structural - plane array, and update the number of triangular patches in each second structural plane.

[0069] After step S4 of this application, it also includes aggregating the boundaries of the second structural plane, specifically including:

[0070] As Figure 6 and Figure 7 shown, first find all the triangular patches in the triangular mesh that do not form the second structural plane, that is, the background triangular patches of this application. The boundaries of the second structural plane formed after non - aggregation are likely to form jagged and uneven edges, which bite together with the background triangular patches, resulting in a poor display effect of the structural plane, being not conducive to further analysis, nor to 3D display and image output. It is necessary to process the non - smooth boundaries through a boundary aggregation algorithm. Specifically as follows:

[0071] Step S41: Sequentially obtain the background triangular patches. If the adjacent triangular patches at two boundaries of the background triangular patch belong to a certain second structural plane, fill the current background triangular patch into this second structural plane;

[0072] Step S42: Analyze the background triangular patches in a loop according to step S41 until there are no background triangular patches that meet the conditions.

[0073] Through this method, process the Figure 6 right - hand side surface model, and the generated effect of the triangular - patch structural plane is as Figure 7 shown. It can be seen from the figure that the marked boundary part of the structural plane becomes smoother.

[0074] As Figure 8 and Figure 9As shown in the figure, the present application further includes step S5, which automatically merges a plurality of the second structural planes to identify a third structural plane with a complete area. For the second structural planes formed on the tunnel construction excavation surface, affected by surrounding rock characteristics, rock mass geological structure characteristics, blasting drilling, blasting construction, mechanical excavation, etc., they generally exhibit the following characteristics:

[0075] (1) The area of a single second structural plane is relatively small;

[0076] (2) The second structural plane protrudes from or is recessed into the excavation surface, as shown in Figure 8 the figure.

[0077] (3) The integrity of the second structural plane is not strong. On the entire construction excavation surface, the exposed areas of the same second structural plane are not directly suitable for forming a complete structural plane visualization area, as shown in Figure 9 the figure.

[0078] (4) The second structural plane is approximately on a plane within the tunnel construction excavation area. The area of the tunnel construction excavation surface is generally within 100 m 2 , which is relatively small compared to large-scale projects such as slopes and foundation pits. Within this range, the second structural plane is usually an approximate plane.

[0079] According to the recognition results of the second structural planes, the rock strata attitudes of the second structural planes with the same color are basically similar, but they may not be the small structural planes on the same second structural plane. It is necessary to analyze the mutual relationship of each second structural plane area of the rock mass to achieve the automatic merger of the second structural planes.

[0080] Step S51: Set the small structural planes in each second structural plane as Si, and perform plane fitting on Si in sequence to obtain the fitting plane Pi; the CGAL algorithm library is used to establish a data structure for plane fitting, and the least squares method is used for the calculation of the plane fitting equation; among them, the small structural plane is a part of the second structural plane.

[0081] Let the equation of the fitting plane Pi be: Ax + By + Cz + D = 0 (C ≠ 0);

[0082] Then

[0083] Let:

[0084] Then: z = a0x + a1y + a2;

[0085] The coordinates of n points (n ≥ 3) on the known structural plane: (xi, yi, zi), i = 0, 1,..., n - 1;

[0086] Let Let k = 0, 1, 2;

[0087] That is:

[0088] Then:

[0089] That is:

[0090] Solve the system of equations to obtain a0, a1, and a2. At this time, the plane direction z = a0x + a1y + a2 is obtained.

[0091] Step S52: Analyze the spatial position relationship between the current structural plane and other structural planes in sequence to determine whether to merge. The selected parameters are as follows:

[0092] 1) Angle between the normal vectors of the structural planes

[0093] If each exposed small structural plane S i is on the same second structural plane, the shape of the second structural plane should be basically the same. Let the normal vectors of the small structural plane S i and the small structural plane S j be N i and N j respectively, where j = 0, 1,..., n - 1. Calculate the angle α i between N j and N ij . If α ij is less than the given threshold α T , then continue with the next analysis, as shown in Appendix Figure 10 .

[0094] The calculation formula for the angle α ij between the normal vectors is as follows:

[0095]

[0096] Set the threshold α T . If α ij < α T , then the two structural planes S i , S j may correspond, and continue with the next analysis.

[0097] 2) Minimum net distance between the fitting planes of the structural planes

[0098] Calculate the distances from each boundary point on the structural plane S j to the fitting plane P i of the structural plane S i respectively, and calculate the average value of the distance values If is less than the given threshold , then continue with the next analysis, as shown in Appendix Figure 11 .

[0099] 3) Structural surface plane distance

[0100] The structural surface S i , S j The points on the outer boundary are projected onto the fitting plane P i , let the projection boundary polygons be L i , L j , if L i , L j If they do not intersect, then calculate L i , L j The minimum distance D ij , if it is less than a given threshold D T , then proceed to the next step of analysis.

[0101] 4) Optimal structural surface merging

[0102] make Calculate the structural surface S i And the R value between all structural faces that may be merged with it, find the minimum value R min The structural surface corresponding to this value is the one that can be compared with S i Merged structural surfaces.

[0103] Step S6: Use the measuring tool to measure the third structural surface and obtain the characteristic information of the third structural surface. Through the above steps, the rock structural surface of the tunnel face can be obtained, which is convenient for measuring and identifying the structural surface using different tools. Specifically, the parameters such as the occurrence, trace length, spacing and opening of different structural surfaces can be accurately obtained, which is helpful for the data statistics of the later tunnel and the exploration before construction.

[0104] A computer device comprises a memory, a processor and a computer program stored in the memory and executable on the processor. When the processor executes the program, the above-mentioned method for identifying the rock structure surface of a tunnel face is implemented.

[0105] A storage medium stores a rock mass structural surface recognition program, which implements the above-mentioned tunnel face rock mass structural surface recognition method when executed by a processor.

[0106] Through computer-readable storage media and computer equipment, the problem of further promotion and application of the tunnel face rock structure surface identification method in engineering is solved, so as to reduce the measurement of redundant structure surfaces and improve the accuracy of structure surface measurement.

[0107] Those skilled in the art should understand that the embodiments of the present invention may provide a method, a system or a computer program product. Therefore, the present invention may take the form of a completely hardware embodiment, a completely software embodiment or an embodiment combining software and hardware aspects. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0108] The above are only the preferred embodiments of the present invention, and do not limit the patent scope of the present invention accordingly. Any equivalent structural transformation made under the inventive concept of the present invention by using the content of the specification and drawings of the present invention, or any direct / indirect application in other related technical fields is included in the patent protection scope of the present invention.

Claims

1. A method for identifying a rock mass structural surface at a tunnel face, characterized in that: The specific steps include: Step S1: Divide the tunnel face rock mass into a number of triangular facets, and obtain the distance between the scanning center point and each triangular facet on the tunnel face; Step S2: setting a distance condition, and screening the distance between the scanning center point and each triangular facet on the tunnel face to obtain triangular facets that meet the distance condition; Step S3: Connectivity analysis and scale screening are performed on the triangular facets that meet the distance condition in turn to form a first structural surface, and the remaining triangular facets that meet the distance condition but fail the connectivity analysis are recorded as background triangular facets; Step S4: filling the triangular facets of the first structural surface to form a second structural surface with connected regions; Step S5: automatically merging the second structural surfaces having regional relationships to identify a third structural surface with a complete region.

2. The method for identifying rock mass structural surface at a tunnel face according to claim 1, characterized in that: The step S1 comprises: Step S11: Get the three vertices (x k ,y k , z k ), k = 1, 2, 3; Step S12: Obtain the plane equation of the triangular face through the three vertices of the triangular face; Step S13: Obtain the distance from the scanning center point to the triangular face by scanning the center point (x0, y0, z0) and the plane equation of the triangular face.

3. The method for identifying rock mass structural surface at a tunnel face according to claim 1, characterized in that: In step S3, the connectivity analysis includes: traversing the triangular facets that meet the distance condition, when the angle between two adjacent triangular facets that meet the distance condition is less than the angle threshold, it indicates that the two adjacent triangular facets that meet the distance condition are connected to each other, and marking the combination of triangular facets that are connected and meet the distance condition as an intermediate structural surface.

4. The method for identifying rock mass structural surface at a tunnel face according to claim 3, characterized in that: In step S3, the scale screening includes setting a quantity threshold, traversing the triangles of the intermediate structure surface, and if the number of triangles is greater than or equal to the quantity threshold, marking that an intermediate structure surface is found; if the number of triangles of the intermediate structure surface is less than the quantity threshold, deleting the current intermediate structure surface.

5. The method for identifying rock mass structural surface at a tunnel face according to claim 1, characterized in that: In the step S4, filling the triangular facets of the first structural surface includes: selecting a hollow triangular facet, and judging the three adjacent faces of the hollow triangular facet; if at least two adjacent faces belong to the same structural surface, the hollow triangular facet is filled into the first structural surface to form a second structural surface, and the three adjacent faces of the hollow triangular facet are judged repeatedly until the hollow triangular facet can no longer be filled into the first structural surface; wherein the hollow triangular facets refer to all triangular facets other than the first structural surface.

6. The method for identifying rock mass structural surface at a tunnel face according to claim 1, characterized in that: After step S4, the process further includes the following: Step S41: obtaining background triangles in sequence, and if the adjacent triangles at two boundaries of the background triangles belong to a second structural surface, filling the current background triangles into the second structural surface; Step S42: Analyze the background triangles in a loop according to step S41 until there are no background triangles that meet the conditions.

7. The method for identifying rock mass structural surface at a tunnel face according to claim 1, characterized in that: The step S5 is also included, wherein the automatic merging of the plurality of the second structural surfaces specifically comprises: Step S51: setting the small structural planes in each second structural plane as Si, and performing plane fitting on Si in turn to obtain fitting planes Pi; Step S52: Analyze the spatial position relationship between the current structural plane and other structural planes in turn to determine whether to merge them.

8. The method for identifying rock mass structural surface at a tunnel face according to claim 7, characterized in that: In step S52, the parameters for determining whether to merge include at least one of the following: the angle of the normal vector of the structural surface, the minimum clear distance between the structural surface fitting planes, the plane distance between the structural surfaces, and the optimal structural surface merging.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the method for identifying the rock structure surface of a tunnel face as described in any one of claims 1 to 8 is implemented.

10. A storage medium, characterized in that: A rock mass structure surface identification program is stored thereon, and when the program is executed by the processor, the method for identifying the rock mass structure surface of a tunnel face as described in any one of claims 1 to 8 is implemented.

Images