A method for calculating deviation of overbreakage contour of slope blasting

Abstract

The application discloses a kind of slope blasting overbreak and underbreak profile deviation calculation method, specifically related to the field of slope blasting measurement, for solving the problem that local anomaly of slope toe is easy to interfere with main deviation judgment;By obtaining the slope design profile model and the post-blast measured profile point set, the post-blast measured profile point set is registered to the slope design profile model to generate the same coordinate profile data;Set profile sampling line along the design slope toe line, organize height sequence according to the design slope direction, and calculate signed deviation along the normal deviation calculation direction of the design slope surface to form the slope surface height direction deviation sequence;Determine the slope toe hole bottom influence zone according to the deviation change pattern that gradually weakens from the slope toe to the upper part of the slope surface, and determine the data that does not fall into the slope toe hole bottom influence zone as the slope surface main data;Calculate the slope surface main overbreak and underbreak deviation and the slope toe local overbreak and underbreak deviation respectively, and check the slope toe hole bottom influence zone boundary attribution, output the slope blasting overbreak and underbreak profile deviation result.

Description

Technical Field

[0001] This invention relates to the field of slope blasting measurement, and more specifically, to a method for calculating the over- and under-excavation profile deviation in slope blasting. Background Technology

[0002] In pre-splitting blasting or smooth blasting of roadbed slopes, the post-blast slope profile usually needs to be compared with the design profile to evaluate whether there is over-excavation or under-excavation. With the application of 3D laser scanning, total station surveying, and photogrammetry in construction inspection, the post-blast slope can form a set of measured post-blast profile points containing a large number of discrete spatial points. Design drawings, design sections, or construction layout data can then form a slope design profile model. The existing processing method generally involves comparing the distance between the post-blast measured profile point set and the slope design profile model after coordinate registration, or calculating the deviation of the measured profile relative to the design slope along several sections, thereby obtaining the over-excavation / under-excavation evaluation results.

[0003] However, in slope blasting operations, post-blast profile deviations are not always uniformly distributed across the entire slope. Because the slope toe is close to the borehole bottom control zone, factors such as drilling depth, borehole bottom alignment, over-depth control, charge structure, filler length, and borehole spacing adjustments can easily create localized outward protrusions, residual rock, rock ledges, or irregular contour shapes near the slope toe. These deviations typically begin to become apparent near the design slope toe line and gradually weaken along the upper slope. Their causes and distribution characteristics differ from the main shaping deviations in the middle and upper parts of the slope.

[0004] If existing calculations of over- and under-excavation profile deviations directly treat the post-blast slope as a whole, local anomalies near the slope toe caused by borehole bottom formation will be mixed with the main slope data in the calculation. When using overall fitting, overall averaging, or unified deviation statistics for the entire slope, local outward protrusions or residual points at the slope toe may change the calculation benchmark of the main slope, amplifying or masking the true deviation of the main slope. When using single-section comparison, local anomalies near the slope toe may also be misjudged as continuous deviations of the entire sampling section, making it impossible to accurately distinguish between the main slope formation problem and the local borehole bottom formation problem at the slope toe.

[0005] Furthermore, when the local anomalies at the slope toe and the main deviation of the slope surface are output together, the results only show the overall over-excavation or under-excavation values, lacking a zonal representation of the deviation source area. For the same construction section, the main slope surface may basically conform to the design slope direction, while the area near the slope toe may have local deviations due to borehole bottom control; or the deviation near the slope toe may be small, while the main slope surface may have an overall offset. If a uniform result expression is still used, it is easy to cause unclear evaluation objects, unclear sources of deviation, and unstable boundary attribution. Therefore, a method for calculating slope blasting over-excavation and under-excavation profile deviation is needed that can identify the influence zone of the borehole bottom at the slope toe under the same coordinate conditions and calculate and output the main slope surface deviation and the local deviation at the slope toe separately.

[0006] To address the aforementioned problems, a technical solution is provided. Summary of the Invention

[0007] To overcome the aforementioned deficiencies of the prior art, embodiments of the present invention provide a method for calculating the over- and under-excavation profile deviation of slope blasting. This method involves acquiring a slope design profile model and a set of post-blast measured profile points, registering the post-blast measured profile points to the slope design profile model, and generating profile data with the same coordinates. Profile sampling lines are set along the design slope toe line, and the height sequence is organized according to the design slope direction. Signed deviations are calculated along the normal deviation calculation direction of the design slope, forming a slope height direction deviation sequence. The influence zone at the bottom of the borehole at the slope toe is determined based on the deviation change pattern that starts from the slope toe and gradually weakens towards the upper part of the slope. Data that does not fall into the influence zone at the bottom of the borehole at the slope toe is identified as the main slope data. The over- and under-excavation deviations of the main slope and the local over- and under-excavation deviations at the slope toe are calculated separately, and the boundary assignment of the influence zone at the bottom of the borehole at the slope toe is verified. The results of the over- and under-excavation profile deviation of slope blasting are then output to solve the problems mentioned in the background art.

[0008] To achieve the above objectives, the present invention provides the following technical solution: S1. Obtain the slope design profile model and the post-blast measured profile point set, register the post-blast measured profile point set to the slope design profile model, and generate profile data with the same coordinates based on the design slope toe line and design slope direction in the slope design profile model. S2, set a contour sampling line along the design slope toe line, so that the contour sampling line extends from the slope toe to the upper part of the slope along the design slope direction. Project the registered measured points on the contour sampling line onto the design slope along the deviation calculation direction to obtain the corresponding design slope reference points. Calculate the signed deviation based on the projection distance between the registered measured points and the corresponding design slope reference points along the deviation calculation direction to generate a slope height direction deviation sequence; where the deviation calculation direction is the normal direction of the design slope. S3. Based on the slope height direction deviation sequence, the local deviation segment that starts from the slope toe and gradually weakens towards the upper part of the slope is identified as the slope toe hole bottom influence zone. The data in the same coordinate contour data that do not fall into the slope toe hole bottom influence zone are identified as the main slope data. S4. Under the condition of keeping the design slope direction unchanged, generate the main deviation benchmark based on the main slope data, and calculate the over-excavation and under-excavation deviation of the main slope. Then, determine the local reference point of the slope toe based on the design slope toe line, the sampling benchmark point of the corresponding contour sampling line, the distance of the corresponding height node, and the design slope direction. Calculate the local over-excavation and under-excavation deviation of the slope toe based on the projection distance between the registered measured point in the influence zone of the slope toe hole bottom and the local reference point of the slope toe along the deviation calculation direction. S5 verifies the boundary assignment based on the deviation change trend on both sides of the upper boundary of the slope toe hole bottom influence zone, and outputs the slope blasting over- and under-excavation profile deviation results, including the over- and under-excavation deviation of the main slope surface and the local over- and under-excavation deviation of the slope toe.

[0009] Furthermore, the post-blast measured profile point set is registered to the slope design profile model, including: The design control point sequence is obtained from the slope design profile model, and the measured control point sequence representing the same field control position is obtained from the post-blast measurement results. Based on the coordinate centers of the design control point sequence and the measured control point sequence, a covariance matrix is ​​constructed. Singular value decomposition is performed on the covariance matrix to obtain the rotation matrix and translation vector. The coordinates of the measured points in the post-explosion measured contour point set are transformed using a rotation matrix and a translation vector to obtain the registered measured contour point set.

[0010] Furthermore, the influence zone at the bottom of the borehole at the slope toe is determined based on the slope height direction deviation sequence, including: The signed deviations on the same contour sampling line are organized into an initial deviation discrimination sequence according to the height nodes, and the maximum absolute value of the signed deviation within the same height node is taken as the node deviation. The initial deviation boundary is determined based on the upper quartile of the deviation of the nodes in the upper part of the slope. The scanning starts from the first height node on the slope toe side. When the node deviation of a height node reaches the initial deviation boundary, and there are no two consecutive adjacent height nodes where the node deviation of the latter height node is greater than that of the former height node, the corresponding height node is assigned to the candidate local deviation segment. When there are two consecutive height nodes where the node deviation is lower than the initial deviation boundary, the scanning stops, and the height node previously assigned to the candidate local deviation segment is used as the upper boundary node of the candidate local deviation segment.

[0011] Furthermore, a main body deviation benchmark is generated based on the main slope data, and the over-excavation and under-excavation deviations of the main slope are calculated, including: Read the signed deviations from the main slope data and organize them into a main deviation input sequence. Determine the main reference offset based on the median offset of the main deviation input sequence. The design slope surface in the slope design profile model is translated along the deviation calculation direction by the main reference offset to generate the main deviation reference. After registering the measured points in the main slope data, project them onto the main deviation benchmark to obtain the corresponding reference points. Based on the projection distance between the registered measured points and the corresponding reference points along the deviation calculation direction, the over- or under-excavation deviation of the main slope is formed.

[0012] Furthermore, contour sampling lines are set along the design slope toe line, including: If the design slope in the slope design profile model is a single slope, then each sampling reference point is expanded along the same design slope direction to form a profile sampling line. If the design slope in the slope design profile model consists of multiple slope segments, the corresponding slope segment is determined according to the design slope toe line segment where the sampling reference point is located, and the profile sampling line is formed by unfolding along the design slope direction of this corresponding slope segment.

[0013] Furthermore, contour sampling lines are set along the design slope toe line, including: The sampling reference point interval is determined based on the projection scalar sequence of the registered measured contour point set in the direction of the design slope toe line; If the projected scalar sequence can form an adjacent difference sequence, then the median of the adjacent difference sequence is taken as the sampling reference point interval. If the projected scalar sequence cannot form an adjacent difference sequence, the design section spacing in the slope design profile model shall be used as the sampling reference point interval. If the slope design profile model does not provide the design section spacing, the hole spacing in the blasting construction parameters shall be used as the sampling reference point interval.

[0014] Furthermore, the calculation methods also include: When the vertical distances from the same registered measured point to two adjacent contour sampling lines are equal, the corresponding registered measured point will be assigned to the contour sampling line with the smaller sampling reference point number. When there are fewer than 2 extractable height nodes on either side of the upper boundary of the slope toe hole bottom influence zone, the corresponding contour sampling line is marked as the boundary line to be maintained, and the boundary assignment determined in step S3 is maintained. When the boundary assignment changes, the main deviation benchmark is regenerated based on the verified main slope data; when the boundary assignment does not change, the main deviation benchmark formed in step S4 is retained.

[0015] Furthermore, the results of the slope blasting over- and under-excavation profile deviation include the profile sampling line assignment, height node order, deviation sign, and upper boundary node after verification corresponding to each registered measured point; Among them, the contour sampling line belongs to the contour sampling line to which the corresponding registered measured point belongs, the height node order indicates the arrangement position of the corresponding registered measured point along the design slope direction from the slope foot to the upper part of the slope, and the deviation sign indicates the over-excavation and under-excavation direction of the corresponding registered measured point relative to the reference object.

[0016] Furthermore, the results of slope blasting over- and under-excavation profile deviations include zoning deviation results; In the zoning deviation results, the over- or under-excavation deviation of the main slope surface corresponds to the main deviation benchmark, and the local over- or under-excavation deviation of the slope toe corresponds to the local reference point of the slope toe determined by the design slope toe line, the corresponding contour sampling line and the design slope direction. For the deviation results of the upper boundary node of the influence zone near the bottom of the hole at the slope toe, the partitioned deviation results retain the height node order and adjacency relationship.

[0017] Furthermore, the boundary assignment is verified based on the deviation variation trend on both sides of the upper boundary of the slope toe hole bottom influence zone, including: Based on the sign of the deviation, the degree of similarity of the deviation amplitude, and the connection relationship of the direction of change on both sides of the upper boundary of the influence zone at the bottom of the slope hole, one of the boundary stability mark and the boundary continuation mark is generated; When generating boundary stability markers, the upper boundary nodes of the slope toe hole bottom influence zone remain unchanged; When generating boundary continuation markers, the continuous height nodes that still appear as local deviations at the toe of the slope will be included in the influence zone at the bottom of the hole after verification, and the updated partition deviation results will be output.

[0018] The technical effects and advantages of the present invention regarding the method for calculating the over-excavation and under-excavation profile deviation in slope blasting are as follows: Before calculating the main slope deviation, this invention first identifies the influence zone of the borehole bottom at the slope toe from the slope height direction deviation sequence. This prevents local outward protrusions, remnants, or irregular lines formed near the slope toe by the borehole bottom from directly entering the main slope data. Compared with the unified comparison method for the entire slope, this processing can reduce the influence of local anomalies at the slope toe on the main deviation benchmark, making the over-excavation and under-excavation deviations of the main slope more stably reflect the main slope forming state.

[0019] This invention lays out contour sampling lines along the designed slope toe line, organizes height nodes in sequence according to the designed slope direction, and organizes signed deviations along the normal deviation calculation direction of the designed slope, ensuring that each contour sampling line retains the deviation arrangement relationship from the slope toe to the upper part of the slope. This arrangement relationship identifies the local deviation segments that begin at the slope toe and gradually decrease, reducing inconsistencies caused by relying solely on manual visual inspection or single-section judgment, and providing a clear data source and processing path for the delineation of the slope toe and borehole bottom influence zone.

[0020] This invention outputs the over- and under-excavation deviations of the main slope surface and the local over- and under-excavation deviations at the slope toe separately, and verifies the boundary assignment of the influence zone at the bottom of the borehole at the slope toe before outputting the results. This result can simultaneously express the main slope surface deviation and the local slope toe deviation, so that the results of slope blasting over- and under-excavation profile deviations are no longer mixed, which helps to reduce the mutual interference of different deviation sources in the same evaluation result. Attached Figure Description

[0021] Figure 1 This is a flowchart illustrating a method for calculating the profile deviation of over-excavation and under-excavation in slope blasting according to the present invention. Figure 2 A schematic diagram of the spatial arrangement of contour sampling lines along the design slope line of the present invention; Figure 3 This is a schematic diagram illustrating the identification of the influence zone at the bottom of the slope toe hole according to the present invention. Detailed Implementation

[0022] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0023] Please see Figure 1 - Figure 3 This invention provides a method for calculating the over- and under-excavation profile deviation in slope blasting, comprising: This invention uses a slope design profile model and a set of post-blast measured profile points as input. First, coordinate registration is used to bring the set of post-blast measured profile points into the design coordinate reference of the slope design profile model, and the design toe line and design slope direction are fixed in the same coordinate profile data. Then, a profile sampling line is set from the design toe line, so that the registered set of measured profile points forms a slope height direction deviation sequence according to the profile sampling line and height nodes. In this sequence, a local deviation segment starting from the toe side and gradually weakening along the upper part of the slope is identified, forming the toe hole bottom influence zone, and the data that does not fall into this zone is used as the main slope data. Then, the main slope deviation reference is generated using the main slope data and the over-excavation and under-excavation deviation of the main slope is calculated. The local over-excavation and under-excavation deviation of the toe hole bottom influence zone is calculated with the design toe line as a reference. Finally, the boundary assignment is checked by the deviation change trend on both sides of the upper boundary of the toe hole bottom influence zone, forming a slope blasting over-excavation and under-excavation profile deviation result containing two types of partition deviations.

[0024] This embodiment provides a method for calculating the over- and under-excavation profile deviation in slope blasting. All the following steps are executed by a data processing server. The data processing server first acquires the slope design profile model and the post-blast measured profile point set. The slope design profile model is a three-dimensional geometric model formed based on slope design drawings, design sections, or construction layout data, used to represent the spatial relationship between the design slope surface, the design slope toe line, and the design slope direction. The post-blast measured profile point set is a set of discrete spatial points on the slope surface collected after slope blasting using three-dimensional laser scanning, total station measurement, or photogrammetry, used to represent the actual formed profile near the slope toe, in the middle of the slope, and at the top of the slope. The data processing server registers the post-blast measured contour point set to the slope design contour model, forming contour data with the same coordinates in the same spatial reference. Then, contour sampling lines are set along the design slope toe line, and the signed deviation of the measured contour relative to the design slope surface is calculated along the design slope surface direction. Then, the influence zone of the slope toe hole bottom is identified based on the change pattern of the deviation in the slope height direction, separating the local anomaly at the slope toe from the main slope surface data. After that, the over-excavation and under-excavation deviations of the main slope surface and the local over-excavation and under-excavation deviations at the slope toe are calculated separately, and the boundary assignment of the influence zone of the slope toe hole bottom is checked before output.

[0025] S1. Obtain the slope design profile model and the post-blast measured profile point set. Register the post-blast measured profile point set to the slope design profile model, and generate contour data with the same coordinates based on the design slope toe line and design slope direction in the slope design profile model; specifically: The data processing server acquires the slope design profile model and the post-blast measured profile point set. The slope design profile model is a geometric model representing the spatial location of the designed slope surface and the location of the designed slope toe. The post-blast measured profile point set is a set of discrete spatial points measured along the slope surface after blasting. The data processing server reads the designed slope toe line and the geometric basis of the designed slope surface adjacent to the designed slope toe line from the slope design profile model. The designed slope toe line serves as the spatial basis for subsequently locating the influence position of the borehole bottom near the slope toe, and the geometric basis of the designed slope surface is used to determine the design slope direction. The design slope direction is the direction along the designed slope surface from the slope toe to the upper part of the slope surface, and is used to subsequently set the profile sampling line and organize the deviation order in the slope height direction.

[0026] The data processing server registers the post-blast measured contour point set to the slope design contour model. During registration, a sequence of design control points that can express the overall position of the slope is first selected in the slope design contour model, and a sequence of measured control points corresponding to the design control point sequence is determined in the post-blast measured contour point set. The design control point sequence and the measured control point sequence are arranged in the same spatial order, and design control vector sequences and measured control vector sequences are constructed respectively. The design control point sequence is obtained from the design slope toe endpoint, slope top control point, construction control pile corresponding point, or pre-deployed measurement target corresponding point in the slope design contour model, which have fixed spatial positions. The measured control point sequence consists of measured points in the post-blast measurement results that represent the same field control position as the aforementioned design control points. The data processing server establishes a sequence relationship between design control points and measured control points according to the one-to-one correspondence of field control positions, so that the... The first design control point and the first Each measured control point corresponds to the same slope toe endpoint, slope crest control point, construction control stake, or measurement target location. If multiple measured points exist at the same field control location, the data processing server takes the arithmetic mean of the coordinates of the measured points within the neighborhood of that field control location as the coordinates of the corresponding measured control point. The data processing server uses a rigid registration algorithm, taking the spatial difference between the design control vector and the measured control vector under the same sequence number as the registration residual, and calculates the rotation matrix and translation vector that reduce the overall registration residual. The rotation matrix is ​​used to correct the direction difference between the measured coordinates and the design coordinates, and the translation vector is used to correct the origin difference between the measured coordinates and the design coordinates.

[0027] When the data processing server performs rigid registration, both the design control point sequence and the measured control point sequence must contain at least three non-collinear control points, and design control points and measured control points with the same sequence number represent the same field control location. Let the measured control point sequence be... The design control point sequence is as follows ,in For the first Spatial coordinates of the measured control points To and The corresponding number The spatial coordinates of the design control points. Assume the number of control points is... The data processing server calculates the coordinate center according to the following relationship: ; ; in, The coordinate center of the measured control point sequence, To design the coordinate center of the control point sequence; Number of control points; For the first Spatial coordinates of one measured control point; For the first Spatial coordinates of the design control points.

[0028] The data processing server constructs the covariance matrix: ; in, The covariance matrix used to calculate the rotation matrix is ​​derived from the centered coordinates of the measured control point sequence and the design control point sequence; For the first The vector of each measured control point relative to the coordinate center of the measured control point sequence; For the first A vector relative to the coordinate center of the design control point sequence; symbol This indicates the matrix transpose.

[0029] Data processing server for covariance matrix Performing singular value decomposition, we obtain: ; in, and The orthogonal matrix obtained from singular value decomposition. This is a singular value matrix. Rotation matrix. Translation vector Determined according to the following relationship: ; If the calculation yields satisfy ,in Representing the rotation matrix If the determinant is such that the data processing server will... Invert the column vector corresponding to the minimum singular value and then recalculate. ,make It is a positive value. This processing is used to ensure... This indicates rotation in the coordinate direction, not a mirror reflection.

[0030] ; in, The rotation matrix is ​​used to transform the post-explosion measured contour point set from the measurement coordinate direction to the design coordinate direction; This is a translation vector used to transform the post-explosion measured profile point set from the measurement coordinate origin to the design coordinate origin. Then, the data processing server processes each measured point in the post-explosion measured profile point set according to the following relationship. Perform coordinate transformation: ; in, For actual measurement points The measured points obtained after rigid registration.

[0031] After obtaining the rotation matrix and translation vector, the data processing server applies the same transformation relationship to every measured point in the post-blast measured contour point set, generating a registered measured contour point set. The point order in the registered measured contour point set is consistent with that in the post-blast measured contour point set, and the coordinate reference is converted to the design coordinate reference of the slope design contour model. During on-site implementation, after blasting, the slope surface can be measured using 3D laser scanning or a total station to obtain multiple measured points near the slope toe, the middle of the slope, and the upper part of the slope. The data processing server calculates the rotation matrix and translation vector based on the correspondence between the slope top control point, slope toe control point, or construction control stake, and then moves the entire post-blast point cloud into the design coordinate system, so that the measured points near the slope toe and the design slope toe line are in the same comparable space.

[0032] Within the same design coordinate datum, the data processing server reads the geometric basis of the design slope surface adjacent to the design slope toe line in the slope design profile model, establishes a direction pointing from the design slope toe line to the upper part of the design slope surface, and determines the design slope direction. If the design slope surface is a plane, the design slope direction is determined by the shortest slope direction from the design slope toe line to the upper boundary of the design slope surface; if the design slope surface is a design surface controlled by a polyline, the data processing server determines the design slope direction along the slope direction corresponding to each segment of the design slope surface, and ensures that the direction of each segment of the design slope surface points from the slope toe to the upper part of the slope surface.

[0033] The coordinated contour data is formed by the slope design contour model, the design toe line, the design slope direction, and the registered set of measured contour points. When generating the coordinated contour data, the data processing server does not change the position of the design toe line in the slope design contour model, does not change the design slope direction, and does not perform over- or under-excavation judgment in this step. It only saves the registered set of measured contour points according to the design coordinate reference and places each measured point in the same spatial reference as the slope design contour model. Step S1 outputs the coordinated contour data, which can be directly called in step S2.

[0034] S2, set a contour sampling line along the design slope toe line, extending the contour sampling line from the slope toe towards the upper part of the slope surface along the design slope direction. Project the registered measured points on the contour sampling line onto the design slope surface along the deviation calculation direction to obtain the corresponding design slope surface reference points. Calculate the signed deviation based on the projection distance between the registered measured points and the corresponding design slope surface reference points along the deviation calculation direction, generating a slope height direction deviation sequence; where the deviation calculation direction is the normal direction of the design slope surface; specifically: The data processing server uses the same-coordinate contour data output in step S1 as the starting point for processing. It reads the design slope toe line, design slope direction, and the registered measured contour point set from the same-coordinate contour data. The data processing server organizes the design slope toe line into ordered line segments according to its extension direction and sets sampling reference points on the design slope toe line. The sampling reference points are the starting points located on the design slope toe line and used to generate the contour sampling lines, arranged sequentially along the extension direction of the design slope toe line.

[0035] When the data processing server determines the sampling reference point interval, it first selects an endpoint on the design slope toe line as the starting point of the slope toe line. And determine the unit vector of the slope toe line direction based on the extension direction of the designed slope toe line. For any registered measured point in the set of registered measured contour points... The data processing server calculates its projected scalar in the direction of the design slope toe line according to the following relationship: ; in, For the first The projection scalar of a registered measured point on the design slope toe line direction is used to represent the position of the registered measured point along the extension direction of the design slope toe line. For the first Spatial coordinates of the registered points; The starting point for designing the slope toe line; The unit vector for the design of the slope toe line direction; symbol This represents the vector dot product. The data processing server will convert each projected scalar... Arrange the values ​​in ascending order, and calculate the difference between adjacent projected scalars. The adjacent difference sequence consists of differences between adjacent projected scalars that are greater than 0, and the median of this adjacent difference sequence is taken as the sampling reference point interval. When the number of projected scalars is insufficient to form an adjacent difference sequence, or when there are no adjacent projected scalar differences greater than 0, the data processing server uses the existing design section spacing in the slope design profile model as the sampling reference point interval; if the slope design profile model does not provide design section spacing, the hole spacing in the blasting construction parameters is used as the sampling reference point interval. The sampling reference points start from one end of the design slope toe line and are set sequentially along the design slope toe line according to the sampling reference point interval until the construction section corresponding to the design slope toe line is covered.

[0036] The data processing server uses each sampling reference point as the starting point of a contour sampling line, extending the contour sampling line along the design slope direction from the toe of the slope to the upper part of the slope. The contour sampling line is a spatial line extending along the design slope direction from the sampling reference point on the design slope toe line, used to represent the contour changes from the toe of the slope to the upper part of the slope at that location. If the... The sampling reference point corresponding to each contour sampling line is The design slope direction unit vector is Then the contour sampling line is represented as: ; in, For the first Points on the contour sampling line; For the first One sampling reference point; For the first The unit vector of the design slope direction corresponding to each contour sampling line; This is the distance parameter extending from the toe of the slope to the upper part of the slope along the design slope direction.

[0037] If the design slope in the slope design profile model is a single slope, all sampling reference points are expanded along the same design slope direction; if the design slope is composed of multiple slope segments, the data processing server determines the corresponding slope segment based on the design slope toe line segment where the sampling reference point is located, and expands the profile sampling line along the design slope direction of that slope segment.

[0038] After the contour sampling lines are generated, the data processing server sends the registered measured contour point set to the neighborhood of the contour sampling lines to determine the measured contour points corresponding to each contour sampling line. For any registered measured point, the data processing server calculates the positional relationship of the registered measured point in the extension direction of the design slope toe line and assigns it to the nearest contour sampling line; when a registered measured point is located between two adjacent contour sampling lines, its assignment is determined according to its vertical distance to the two adjacent contour sampling lines, and the contour sampling line with the smaller distance receives the registered measured point.

[0039] For the registered measured points and the Contour sampling line The data processing server calculates the vertical distance from the point to the contour sampling line. If the contour sampling line is represented as ,but: ; in, For the actual measurement points after registration To the The vertical distance of the contour sampling lines; For the first The sampling reference point for the contour sampling line; For the first The unit vector of the design slope direction corresponding to each contour sampling line; Indicates the vector length. The data processing server will... Into The smallest contour sampling line. When the vertical distance from the same registered measured point to two adjacent contour sampling lines is equal, the data processing server assigns the registered measured point to the contour sampling line with the smaller sequence number according to the arrangement order of the sampling reference points on the design slope toe line. The above assignment process is only used to determine the subset of measured contour points to which the registered measured point belongs, and does not change the spatial coordinates of the registered measured point in the contour data of the same coordinate system. After the assignment is completed, each contour sampling line obtains a corresponding subset of measured contour points, and the registered measured points in the subset of measured contour points are arranged according to the order from the slope toe to the upper part of the slope along the design slope direction.

[0040] The data processing server establishes a design slope reference point sequence on each contour sampling line. This sequence consists of points located on the design slope surface of the slope design contour model, arranged from the slope toe to the top of the slope in the design slope direction. These points are used to establish a positional correspondence with the subset of measured contour points. When establishing the design slope reference point sequence, the data processing server projects each registered measured point from the subset of measured contour points onto the design slope surface of the slope design contour model along the deviation calculation direction. The projected point serves as the corresponding design slope reference point.

[0041] The deviation calculation direction is the normal direction of the designed slope. In this embodiment, the deviation calculation direction... The positive direction is defined as the direction from the design slope surface towards the excavated space. According to this notation, a positive sign indicates that the measured profile deviates from the design slope surface towards the excavated space, corresponding to under-excavation; a negative sign indicates that the measured profile deviates from the design slope surface towards the interior of the slope, corresponding to over-excavation. This notation is also used for subsequent over-excavation / under-excavation deviations of the main slope surface and local over-excavation / under-excavation deviations at the slope toe.

[0042] When the designed slope is a plane, the data processing server selects any point on the designed slope. As the reference point for the design slope, the first... Registered measured points Corresponding design slope reference point : ; in, for The design slope reference point is obtained by projecting it onto the design slope along the deviation calculation direction; For the first Spatial coordinates of the registered points; For any point on the slope, it is derived from the slope design outline model; This is the unit vector for calculating the deviation. express The projected distance relative to the design slope along the deviation calculation direction. When the design slope consists of multiple slope segments, the data processing server first... The corresponding slope segment is determined by the contour sampling line, and then the above relationship is calculated within that slope segment. .

[0043] For any contour sampling line, the first For each registered measured point, the data processing server calculates the signed bias according to the following relationship: ; in, For the first The signed deviation corresponding to each registered measured point comes from a subset of measured contour points on the current contour sampling line, and is used to represent the over- or under-excavation direction and deviation of the measured contour relative to the design slope. For the first The spatial coordinates of the registered measured points are derived from the registered measured contour point set formed in step S1. for The design slope reference point, obtained by projecting it onto the design slope along the deviation calculation direction, is used to represent the design position of the same slope height. The unit vector representing the direction of deviation calculation originates from the normal direction of the designed slope surface in the slope design profile model; symbol This represents the dot product of vectors, used to obtain... relatively The projection along the deviation calculation direction.

[0044] The data processing server uses the sampling reference point as the zero starting point to calculate the sequential position of each design slope reference point relative to the sampling reference point along the contour sampling line. For the first... The first on the contour sampling line Design slope reference points Their order is determined according to the following relationship: ; in, For the first Design slope reference points Along the first The contour sampling line is relative to the sampling reference point The order of positions; For the first The unit vector of the design slope direction corresponding to each contour sampling line. The data processing server follows... Arrange the corresponding signed deviations from smallest to largest to generate a sequence of slope height direction deviations.

[0045] The slope height direction deviation sequence is a signed deviation sequence arranged along the design slope direction from the slope toe to the upper part of the slope. It is used to represent the shape of the deviation of the measured profile relative to the design slope as a function of slope height. When generating the slope height direction deviation sequence, the data processing server first reads the sequential positions of each design slope reference point within the same profile sampling line. and in accordance with Arrange the samples in ascending order; then calculate the positive difference between adjacent sequential positions, and take the median of this difference sequence as the height node interval of the contour sampling line. The height node interval is positive. When the number of sequential positions within the same contour sampling line is insufficient to form an adjacent difference sequence, or when the height node interval obtained from the sequential position difference sequence is not greater than 0, the data processing server uses the existing design sampling interval within the slope segment corresponding to the contour sampling line as the height node interval. If the design sampling interval is not given for the corresponding slope segment, the scalar length of the sampling reference point interval shall be used as the height node interval. For the first For each design slope reference point, the data processing server determines its corresponding height node number according to the following relationship: ; in, For the first The height node number corresponding to each design slope reference point; For the first The sequential positions of the design slope reference points relative to the sampling benchmark points along the contour sampling line; For the first The height node spacing of the contour sampling line; This indicates rounding down. The data processing server groups signed deviations with the same height node number into the same height node and retains a set of signed deviations within the same height node. Step S2 outputs the slope height direction deviation sequence corresponding to each contour sampling line, for step S3 to identify the influence zone at the bottom of the hole at the slope toe.

[0046] S3, based on the slope height direction deviation sequence, the local deviation segment starting from the slope toe and gradually weakening towards the upper part of the slope is identified as the slope toe borehole bottom influence zone, and the data in the same coordinate contour data that do not fall into the slope toe borehole bottom influence zone are identified as the main slope data; specifically: The data processing server uses the slope height direction deviation sequence output in step S2 as the starting point for processing and continues to reference the same coordinate contour data output in step S1. The data processing server reads the slope height direction deviation sequence line by line according to the association relationship of the contour sampling lines, so that the deviation values ​​on each contour sampling line maintain the arrangement order starting from the design slope toe line and pointing to the upper part of the slope along the design slope direction.

[0047] The data processing server uses height nodes in the slope height direction deviation sequence as analysis units. Height nodes are position nodes arranged according to the design slope direction in the slope height direction deviation sequence, representing the order of slope height relative to the design slope toe line on the same profile sampling line. Each height node corresponds to one or more signed deviations. When multiple signed deviations exist for the same height node, the data processing server first preserves their internal and external relationships according to the deviation signs, and then uses the maximum absolute value of the deviations within that height node as the node deviation. The node deviation is the maximum degree of deviation of the measured profile relative to the design slope at the same height node, used in the continuous discrimination of the initial local deviation segment at the slope toe. The node deviation of each height node is determined according to the following relationship: ; in, For the first The node deviation of each height node; For the first The set of signed biases within each height node; Signed bias The absolute value of.

[0048] When multiple signed deviations exist within the same height node, the sign corresponding to the signed deviation with the largest absolute value is taken as the node sign of that height node; if the absolute values ​​of positive and negative deviations are the same and both are the largest absolute values ​​within that height node, then the positive and negative deviations of that height node are retained in the initial deviation discrimination sequence of the corresponding signs, and participate in subsequent continuous discrimination and strip region connection respectively.

[0049] The data processing server organizes the height nodes on each contour sampling line into an initial deviation discrimination sequence, arranged from the toe of the slope to the upper part of the slope. This initial deviation discrimination sequence consists of the height node order, the node deviation amount, and the corresponding signed deviation sign, used to determine whether the deviation originates and concentrates near the design toe of the slope. Subsequently, the data processing server uses a continuous initial segment scanning method to identify local deviation segments. This method starts from the height nodes on the toe side, compares the continuity and reduction relationship of node deviations point by point along the slope height direction, and outputs a sequence identification method for initial continuous node segments that satisfy the borehole bottom influence morphology.

[0050] The gradual weakening towards the upper part of the slope means that the node deviation of adjacent height nodes along the design slope direction satisfies that the node deviation of the next height node is not greater than the node deviation of the previous height node, or only local undulations that do not constitute a continuous increase appear; where continuous increase means that two consecutive adjacent height nodes satisfy the condition that the node deviation of the next height node is greater than the node deviation of the previous height node.

[0051] Before executing the continuous initial segment scanning method, the data processing server determines the initial deviation boundary for each contour sampling line. The initial deviation boundary is the lower limit of the node deviation used to determine whether a height node has entered a local deviation segment, and is determined by the upper quartile of the node deviations located in the upper-middle part of the slope within the same contour sampling line, arranged in ascending order. For lines containing... The initial deviation discrimination sequence of each height node, the data processing server will set the sequence number to be no less than... The height node is determined to be a node in the upper middle part of the slope, where... The total number of height nodes on the same contour sampling line. Indicates to The node number is rounded up. The data processing server reads the node deviations corresponding to the upper nodes in the slope, sorts them in ascending order, and takes the node deviation at the 75th percentile as the starting deviation boundary. If the number of nodes in the upper part of the slope is... The data processing server sorts the deviations of the upper nodes on the slope from smallest to largest, and then... Determine the 75th percentile position, where, The sorted value sequence number. This represents the number of nodes in the upper part of the slope. This indicates rounding up. The data processing server retrieves the integer part after sorting. The deviation of each node serves as the initial deviation boundary. This boundary is only used in step S3 to identify local deviation segments and does not change the value of the sign deviation in step S2.

[0052] When the data processing server executes the continuous starting segment scanning method, the deviation of nodes on the same contour sampling line is recorded as... ,in For the first The node deviation of each height node. The order increases from the toe of the slope to the upper part of the slope. If the... The height node and the first Each height node satisfies If the deviation decreases along the upper part of the slope, then it is determined that the deviation decreases; if the condition is met... If no two consecutive adjacent height nodes satisfy the condition that the deviation of the latter height node is greater than that of the former height node, it is considered a local fluctuation. If two consecutive adjacent height nodes satisfy the condition that the deviation of the latter height node is greater than that of the former height node, it is considered a continuous increase. The continuous starting segment scanning method starts from the first height node. When the node deviation of a continuous height node reaches the starting deviation boundary and is not considered a continuous increase, it is classified into a candidate local deviation segment. When the node deviation of two consecutive height nodes is lower than the starting deviation boundary, the scanning stops, and the height node previously classified into the candidate local deviation segment is used as the upper boundary node of the candidate local deviation segment.

[0053] When the initial height node does not reach the initial deviation boundary, the data processing server records the local deviation segment of the contour sampling line as an empty segment. When the initial height node reaches the initial deviation boundary but the deviation of subsequent nodes continues to increase along the upper slope, the data processing server treats this deviation pattern as a manifestation of deviation in the main direction of the slope, and does not treat it as a local deviation segment that starts from the slope toe and gradually weakens towards the upper slope. During field implementation, if residual rock ridges are formed near the slope toe due to insufficient control of hole bottom leveling, the node deviation of the initial height node of the contour sampling line is relatively large. As the height node moves towards the upper slope, the node deviation gradually decreases, and the data processing server marks this initial continuous node segment as a candidate local deviation segment.

[0054] The data processing server organizes candidate local deviation segments of each contour sampling line along the design slope toe line direction, arranging the candidate local deviation segments on adjacent contour sampling lines according to the spatial order of the design slope toe line, and determines the connection relationship of adjacent candidate local deviation segments in the slope height direction. The same slope toe starting characteristic is determined by the starting height node and starting deviation sign of the candidate local deviation segment. If the candidate local deviation segments on adjacent contour sampling lines all start from the first height node on the slope toe side, or all start from the slope toe side height node adjacent to the first height node on the slope toe side, and the signed deviation signs corresponding to the starting height nodes of each candidate local deviation segment are consistent, then the data processing server determines that the adjacent candidate local deviation segments have the same slope toe starting characteristic. If the starting deviation signs of adjacent candidate local deviation segments are different, even if their upper boundary nodes are continuous in the slope height direction, the data processing server segments them to form different slope toe hole bottom influence zones. When determining whether the candidate local deviation segments of adjacent contour sampling lines are continuous, the data processing server reads the upper boundary node sequence number of the candidate local deviation segments on the two adjacent contour sampling lines. If the difference between the indices of two upper boundary nodes is no greater than one height node, they are determined to be continuously distributed in the slope height direction and are included in the same slope toe hole bottom influence zone. If the difference between the indices of two upper boundary nodes is greater than one height node, they are determined to be discontinuous in the slope height direction and are segmented according to the design slope toe line to form different slope toe hole bottom influence zones. For cases where the height node spacing is inconsistent, the data processing server uses the median distance between adjacent height nodes along the design slope direction as the height node spacing and compares the distance difference between two upper boundary nodes along the design slope direction with this height node spacing. When the distance difference is no greater than the height node spacing, it is determined to be a continuous distribution.

[0055] The data processing server connects candidate local deviation segments with the same slope toe starting characteristics and whose upper boundary nodes are continuously distributed between adjacent contour sampling lines into a strip-shaped region, forming the slope toe hole bottom influence zone. The slope toe hole bottom influence zone is a coordinate-based data area that starts near the design slope toe line, extends upwards along the design slope direction, and is jointly defined by the local deviation segments on each contour sampling line. The lower boundary of this area is based on the design slope toe line, the upper boundary is formed by connecting the upper boundary nodes of the candidate local deviation segments on each contour sampling line, and the lateral boundary is determined by the continuous range of the contour sampling line with candidate local deviation segments in the direction of the design slope toe line.

[0056] The data processing server performs a consistent partitioning of the registered measured contour point set, the signed deviation point sequence on the contour sampling line, and the slope height direction deviation sequence based on the spatial range of the slope toe hole bottom influence zone in the same coordinate contour data. For each registered measured point, the data processing server determines whether the registered measured point is located within the lower, upper, and lateral boundaries of the slope toe hole bottom influence zone based on its corresponding contour sampling line and height node. If it is within these boundaries, the registered measured point is included in the slope toe hole bottom influence zone; if it is not within these boundaries, it is included in the main slope data. The main slope data consists of data in the same coordinate contour data that does not fall within the slope toe hole bottom influence zone, and is used in step S4 to generate the main deviation benchmark while keeping the design slope direction unchanged. Step S3 outputs the slope toe hole bottom influence zone, the main slope data, and the upper boundary node of the slope toe hole bottom influence zone.

[0057] S4, while keeping the design slope direction unchanged, generate the main deviation benchmark based on the main slope data, and calculate the over-excavation and under-excavation deviation of the main slope; and determine the local reference point at the slope toe based on the design slope toe line, the sampling benchmark point of the corresponding contour sampling line, the distance of the corresponding height node, and the design slope direction. Calculate the local over-excavation and under-excavation deviation at the slope toe based on the projection distance between the registered measured point within the influence zone at the bottom of the slope toe hole and the local reference point at the slope toe along the deviation calculation direction; specifically: The data processing server uses the slope main data output in step S3 as the starting point for processing, and continues to reference the slope design profile model, design slope direction, and coordinate profile data output in step S1. The data processing server reads the slope main data according to the profile sampling line attribution relationship, so that the slope main data maintains the correspondence between its source profile sampling line, corresponding height node, and signed deviation.

[0058] The main deviation benchmark is a calculation reference surface generated based on the offset state of the main slope data relative to the design slope, while keeping the design slope direction unchanged. It is used to calculate the over-excavation and under-excavation deviations of the main slope data. When generating the main deviation benchmark, the data processing server does not change the slope aspect of the design slope in the slope design profile model, but instead counts the overall offset of the main slope data relative to the design slope along the deviation calculation direction determined in step S2. The data processing server reads the signed deviations corresponding to each registered measured point from the main slope data and organizes these signed deviations into a main deviation input sequence according to the contour sampling line and height node order. The main deviation input sequence is a sequence composed of signed deviations in the main slope data, excluding deviations within the influence zone of the slope toe hole bottom, and is used to calculate the translation amount of the main deviation benchmark.

[0059] The data processing server uses the median offset estimation algorithm to calculate the main reference offset. The median offset estimation algorithm sorts the main deviation input sequence by value and takes the middle position deviation value. When the main deviation input sequence contains an even number of signed deviations, the average of the two middle adjacent deviation values ​​is taken as the main reference offset. When the slope main data is empty or there are no signed deviations in the main deviation input sequence, the data processing server initializes the main reference offset to 0 and uses the design slope in the slope design profile model as the main deviation reference. In this case, the slope main over-excavation / under-excavation deviation output is empty, but the local over-excavation / under-excavation deviation of the slope toe hole bottom influence zone is still calculated according to the design slope toe line and the local slope toe reference point. The main reference offset is a representative offset of the slope main data relative to the design slope, used to translate the design slope along the deviation calculation direction. The main deviation reference is generated according to the following relationship: ; in, As the main deviation benchmark, it represents the design slope. Every point on Along the deviation calculation direction Translation of the main body reference offset The resulting geometric surface; For designing slope any point on; The design slope surface in the slope design profile model is derived from the same coordinate profile data formed in step S1. The main reference offset is derived from the main deviation input sequence in the main slope data and is used to represent the representative offset of the main slope relative to the design slope. The unit vector for the deviation calculation direction follows the normal direction of the designed slope in step S2. This formula indicates that the main deviation reference is only translated normally relative to the designed slope, while the direction of the designed slope remains unchanged.

[0060] Subsequently, the data processing server calculates the over- and under-excavation deviations of the slope main body based on the main body deviation benchmark. For each registered measured point in the slope main body data, the data processing server first determines the corresponding reference point on the main body deviation benchmark for that registered measured point. Registered measured points The data processing server first obtains the corresponding design slope reference point according to the projection method of the design slope reference point in step S2. Then, determine the corresponding reference point on the main deviation benchmark according to the following relationship: ; in, for The corresponding reference point on the main deviation benchmark; for The design slope reference point is obtained by projecting it onto the design slope along the deviation calculation direction; This is the main reference offset; This is the unit vector for calculating the deviation direction.

[0061] For the first in the main data of the slope The deviation of the main body point after registration and actual measurement is calculated according to the following relationship: ; in, For the first The deviation of the main point corresponding to each registered measured point comes from the main slope data and is used to constitute the over- or under-excavation deviation of the main slope. For the first The spatial coordinates of each registered measured point are derived from the set of registered measured contour points in the same coordinate contour data. for The corresponding reference point on the main body deviation datum is used to indicate the main body datum position at that location; The unit vector for calculating the deviation direction continues to use the direction definition from step S2; symbol This represents the vector dot product, used to obtain the signed projection along the deviation calculation direction. After the data processing server completes the deviation calculation for all main points, it forms the over- or under-excavation deviation of the slope main body according to the deviation sign, and retains the outline sampling line assignment and height node order.

[0062] The data processing server further calculates the local over- or under-excavation deviation at the toe of the slope based on the slope toe and borehole bottom influence zone output in step S3, and by referencing the design slope toe line, design slope direction, and coordinate contour data output in step S1. The data processing server uses the registered measured points within the slope toe and borehole bottom influence zone as the calculation object, reading their spatial positions according to the contour sampling line and height node to which these registered measured points belong, ensuring that the local deviation calculation of the slope toe and borehole bottom influence zone is consistent with the sampling order formed in step S2.

[0063] Before calculating the local over- or under-excavation deviation at the toe of the slope, the data processing server determines a local reference point at the toe for each registered measured point within the influence zone of the borehole bottom at the toe. The local reference point is a point located on the design toe line or at the design toe reference position formed by the design toe line along the design slope direction, used to indicate the design starting position corresponding to the measured point within the influence zone of the borehole bottom. For the first... After registration, the data processing server first determines the contour sampling line to which each measured point belongs, and records the sampling reference point of that contour sampling line on the design slope toe line as... Then read the height node corresponding to the registered measured point along the design slope direction relative to the actual height node. The distance is denoted as The unit vector of the design slope direction corresponding to the contour sampling line is denoted as... The local reference point at the toe of the slope is determined according to the following relationship: ; in, For the first Local reference points at the toe of the slope corresponding to the registered measured points; This is the sampling reference point on the design slope toe line of the contour sampling line to which the registered measured point belongs; This is the distance of the registered measured point relative to the sampling reference point along the design slope direction; This is the unit vector of the design slope direction corresponding to the sampling line of the profile. If the slope design profile model is a segmented slope, this is the unit vector of the design slope direction. The slope is divided into sections according to the design slope direction within the corresponding slope segment.

[0064] For the first in the influence zone of the slope toe hole bottom After registration, the local point deviation of the measured points is calculated according to the following relationship: ; in, For the first The local point deviation corresponding to each registered measured point originates from the influence zone of the bottom of the hole at the toe of the slope and is used to constitute the local over- or under-excavation deviation at the toe of the slope. For the first The spatial coordinates of each registered measured point are derived from the set of registered measured contour points in the same coordinate contour data. for The corresponding local reference point at the toe of the slope is derived from the design reference position jointly determined by the design toe line, the outline sampling line, and the design slope direction. The unit vector for the deviation calculation direction continues to use the direction definition from step S2. After the data processing server completes the local point deviation calculation for all registered measured points within the influence zone of the slope toe hole bottom, it forms the local over-excavation / under-excavation deviation at the slope toe and retains the correspondence between it and the boundary nodes, contour sampling lines, and design slope toe line direction ranges of the influence zone at the slope toe hole bottom.

[0065] The data processing server organizes the over- and under-excavation deviations of the main slope surface and the local over- and under-excavation deviations at the slope toe into zoned deviation results. In the zoned deviation results, the over- and under-excavation deviations of the main slope surface correspond to the main deviation benchmark, while the local over- and under-excavation deviations at the slope toe correspond to the local reference points at the slope toe determined by the design slope toe line, the corresponding contour sampling line, and the design slope direction. The two types of results are distinguished in terms of reference objects, but remain consistent in terms of deviation calculation direction and symbol conventions. For deviation results near the upper boundary node of the influence zone at the bottom of the slope toe hole, the data processing server retains its height node order and adjacency relationship for step S5 to read the deviation change trend on both sides of the upper boundary.

[0066] S5, based on the deviation change trend on both sides of the upper boundary of the slope toe hole bottom influence zone, verifies the boundary assignment, and outputs the slope blasting over- and under-excavation profile deviation results, including the main slope over- and under-excavation deviation and the local over- and under-excavation deviation at the slope toe; specifically: The data processing server takes the partition deviation result output in step S4 as the starting point for processing, and continues to reference the upper boundary nodes of the slope toe hole bottom influence zone output in step S3, the contour sampling lines formed in step S2, and the height node sequence. The data processing server reads the upper boundary nodes of the slope toe hole bottom influence zone one by one according to the contour sampling lines, so that each upper boundary node maintains its position in the slope height direction on the corresponding contour sampling line.

[0067] The data processing server extracts deviation values ​​from the slope toe and upper slope sides along the corresponding contour sampling lines, centered on the upper boundary node of the slope toe borehole bottom influence zone. Deviation values ​​on the slope toe side originate from local over- or under-excavation deviations at the slope toe, forming an inner-zone deviation sequence; deviation values ​​on the upper slope side originate from the main slope over- or under-excavation deviations, forming a main-side deviation sequence. The inner-zone deviation sequence is a local point deviation sequence near the upper boundary of the slope toe borehole bottom influence zone and still belonging to the slope toe borehole bottom influence zone, used to represent the change pattern of local deviations at the slope toe when approaching the upper boundary; the main-side deviation sequence is a main point deviation sequence near the upper boundary of the slope toe borehole bottom influence zone and belonging to the main slope data, used to represent the change pattern of the main slope deviations outside the upper boundary.

[0068] The data processing server extracts local point deviations of at least two consecutive height nodes towards the slope foot as an inner deviation sequence, centered on the upper boundary node of each contour sampling line, and extracts main point deviations of at least two consecutive height nodes towards the upper slope as a main deviation sequence. When there are fewer than two extractable height nodes on either side, the contour sampling line is marked as a boundary line to be maintained, and the boundary belongs to the result determined in step S3.

[0069] The data processing server calculates the changes in adjacent deviations along the design slope direction for both the inner-band deviation sequence and the main-side deviation sequence, resulting in the inner-band change sequence and the main-side change sequence. The inner-band change sequence indicates whether the local over- or under-excavation deviation at the slope toe still weakens towards the upper slope when approaching the upper boundary, while the main-side change sequence indicates whether the main over- or under-excavation deviation of the slope has entered a relatively smooth main change pattern outside the upper boundary. The data processing server calculates the change sequence by subtracting the absolute values ​​of adjacent deviations along the design slope direction from the slope toe towards the upper slope. For two adjacent deviation values... and Its amplitude change is: ; in, This represents the magnitude change in the deviation between adjacent height nodes; For the first The deviation value of each height node; For the first The deviation value of each height node. If This indicates that the deviation amplitude decreases or remains constant along the upper part of the slope; if This indicates that the deviation amplitude increases along the upper part of the slope.

[0070] The data processing server reads the deviation values ​​closest to the upper boundary from the inner-band deviation sequence and the deviation values ​​closest to the upper boundary from the main-side deviation sequence, and calculates their sign consistency, deviation amplitude similarity, and direction of change correlation. Sign consistency is used to determine whether the two sides of the upper boundary are still in the same over- or under-excavation direction; deviation amplitude similarity is used to determine whether the deviation on the outer side of the upper boundary has significantly deviated from the local deviation at the slope toe; and direction of change correlation is used to determine whether the main-side deviation still continues along the weakening direction of the inner-band deviation. The boundary verification discriminant is the judgment result formed by the above three relationships, and the output form is a boundary stability marker or a boundary continuation marker.

[0071] The data processing server records the deviation value closest to the upper boundary in the inner deviation sequence as... The deviation value closest to the upper boundary in the main body side deviation sequence is denoted as... ,in, This refers to the local point deviation of the height node inside the upper boundary. This represents the subject point deviation between adjacent height nodes outside the upper boundary. The data processing server records the magnitude difference of the subject point deviations in the subject-side deviation sequence as: ; in, The magnitude difference between adjacent height nodes in the main body side deviation sequence; The first in the main side deviation sequence The deviation of the main point of each height node; The first in the main side deviation sequence The main point deviation of each height node. The data processing server will Sort by numerical value in ascending order, and take the value at the 75th percentile as the boundary amplitude tolerance. It is used to determine whether the deviation amplitudes on both sides of the upper boundary belong to the same continuous range of change.

[0072] when and When the signs are the same and the following relationship is satisfied, the data processing server determines that the deviation magnitudes on both sides of the upper boundary are close: ; in, This represents the absolute value of the deviation between adjacent main points on the outer side of the upper boundary. This represents the absolute value of the deviation at a local point inside the upper boundary. This represents the absolute value of the difference in deviation amplitude between the two sides of the upper boundary. This is the boundary magnitude tolerance. If it simultaneously satisfies the self-... The absolute value of the deviation of the next height node after the initial consecutive height nodes is not greater than the absolute value of the deviation of the previous height node. If the sum of the above conditions is met, the data processing server generates a boundary continuation marker; if the above conditions are not met, the data processing server generates a boundary stability marker.

[0073] The data processing server processes the boundary verification discriminant along each contour sampling line, performing a classification verification on the upper boundary nodes of the slope toe and hole bottom influence zone on each contour sampling line. When the boundary verification discriminant is a boundary stability marker, the data processing server keeps the upper boundary nodes of the slope toe and hole bottom influence zone on that contour sampling line unchanged; the data on the original upper boundary slope toe side continues to be classified as belonging to the slope toe and hole bottom influence zone, and the data on the original upper boundary slope upper side continues to be classified as belonging to the slope main data. When the boundary verification discriminant is a boundary continuation marker, the data processing server starts judging point by point from the adjacent height node outside the upper boundary. For the first deviation in the main body side sequence... A height node, if the corresponding main point of the height node deviates Deviation from local point inside the upper boundary The symbols are the same, and they satisfy: ; Meanwhile, if the deviation amplitude of the subsequent height node does not show a continuous increase, then this height node is included in the influence zone of the slope toe hole bottom after verification; when different signs appear, or If, or if the deviation amplitude of two consecutive main body side height nodes increases, the boundary extension is stopped, and the height node of the previously included slope toe hole bottom influence zone after verification is used as the upper boundary node after verification. The first in the main side deviation sequence The main point deviation of each height node, This refers to the local point deviation of the height node inside the upper boundary. This is the allowable range of boundary amplitude.

[0074] After completing the boundary assignment verification of a single contour sampling line, the data processing server compares the verified upper boundary nodes of adjacent contour sampling lines along the design slope toe line direction. If the verified upper boundary nodes of adjacent contour sampling lines are continuously distributed in the slope height direction, they are connected to form the upper boundary of the verified slope toe hole bottom influence zone; if a local contour sampling line is marked as a boundary line to be maintained due to missing data, the upper boundary node determined in step S3 is used at that location, and the original lateral range is maintained with the adjacent valid contour sampling lines. The verified slope toe hole bottom influence zone is the slope toe hole bottom influence zone after the upper boundary assignment verification, and the verified slope main data is the data in the same coordinate contour data that does not fall into the verified slope toe hole bottom influence zone.

[0075] If the boundary assignment changes, the data processing server first reorganizes the main deviation input sequence based on the verified slope main data, and redetermines the verified main reference offset according to the calculation method of the main reference offset in step S4. Then, it regenerates the main deviation benchmark based on the verified main reference offset, and calculates the main point deviation in the verified slope main data according to the regenerated main deviation benchmark. Height nodes assigned to the influence zone of the hole bottom at the slope toe after verification obtain local point deviations according to the calculation method of the local over-excavation / under-excavation deviation at the slope toe in step S4, and are removed from the main over-excavation / under-excavation deviation of the slope. Height nodes still assigned to the main data of the verified slope obtain main point deviations according to the regenerated main deviation benchmark. If the boundary assignment remains unchanged, the data processing server retains the main deviation benchmark, the main over-excavation / under-excavation deviation of the slope, and the local over-excavation / under-excavation deviation at the slope toe formed in step S4. This results in an updated zoning deviation result.

[0076] The data processing server outputs the over- and under-excavation profile deviation results for slope blasting, based on the attribution relationship between the influence zone at the bottom of the borehole at the slope toe and the main data of the slope surface after verification. These results include over- and under-excavation deviations for the main slope surface and local over- and under-excavation deviations at the slope toe. The main slope surface over- and under-excavation deviation corresponds to the main data of the slope surface after verification, using the main deviation benchmark as a reference. The local over- and under-excavation deviation corresponds to the influence zone at the bottom of the borehole at the slope toe after verification, using the local reference point at the slope toe determined by the design slope toe line, the corresponding profile sampling line, and the design slope direction as a reference. The data processing server retains the attribution of the profile sampling line, the order of height nodes, the deviation sign, and the upper boundary node after verification for each deviation result, enabling the main slope surface issues and local borehole bottom issues within the same construction section to be read separately in the same coordinate benchmark.

[0077] The above embodiments are merely preferred embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for calculating the over- and under-excavation profile deviation in slope blasting, characterized in that, Performed by the data processing server, including: S1. Obtain the slope design profile model and the post-blast measured profile point set, register the post-blast measured profile point set to the slope design profile model, and generate profile data with the same coordinates based on the design slope toe line and design slope direction in the slope design profile model. S2, set a contour sampling line along the design slope toe line, so that the contour sampling line extends from the slope toe to the upper part of the slope along the design slope direction. Project the registered measured points on the contour sampling line onto the design slope along the deviation calculation direction to obtain the corresponding design slope reference points. Calculate the signed deviation based on the projection distance between the registered measured points and the corresponding design slope reference points along the deviation calculation direction to generate a slope height direction deviation sequence; where the deviation calculation direction is the normal direction of the design slope. S3. Based on the slope height direction deviation sequence, the local deviation segment that starts from the slope toe and gradually weakens towards the upper part of the slope is identified as the slope toe hole bottom influence zone. The data in the same coordinate contour data that do not fall into the slope toe hole bottom influence zone is identified as the main slope data. S4. Under the condition of keeping the design slope direction unchanged, generate the main deviation benchmark based on the main slope data, and calculate the over-excavation and under-excavation deviation of the main slope. Then, determine the local reference point of the slope toe based on the design slope toe line, the sampling benchmark point of the corresponding contour sampling line, the distance of the corresponding height node, and the design slope direction. Calculate the local over-excavation and under-excavation deviation of the slope toe based on the projection distance between the registered measured point in the influence zone of the slope toe hole bottom and the local reference point of the slope toe along the deviation calculation direction. S5 verifies the boundary assignment based on the deviation change trend on both sides of the upper boundary of the slope toe hole bottom influence zone, and outputs the slope blasting over- and under-excavation profile deviation results, including the over- and under-excavation deviation of the main slope surface and the local over- and under-excavation deviation of the slope toe.

2. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, Register the post-blast measured profile point set to the slope design profile model, including: The design control point sequence is obtained from the slope design profile model, and the measured control point sequence representing the same field control position is obtained from the post-blast measurement results. Based on the coordinate centers of the design control point sequence and the measured control point sequence, a covariance matrix is ​​constructed. Singular value decomposition is performed on the covariance matrix to obtain the rotation matrix and translation vector. The coordinates of the measured points in the post-explosion measured contour point set are transformed using a rotation matrix and a translation vector to obtain the registered measured contour point set.

3. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, The influence zone at the bottom of the borehole at the toe of the slope is determined based on the slope height direction deviation sequence, including: The signed deviations on the same contour sampling line are organized into an initial deviation discrimination sequence according to the height nodes, and the maximum absolute value of the signed deviation within the same height node is taken as the node deviation. The initial deviation boundary is determined based on the upper quartile of the deviation of the nodes in the upper part of the slope. The scanning starts from the first height node on the slope toe side. When the node deviation of a height node reaches the initial deviation boundary, and there are no two consecutive adjacent height nodes where the node deviation of the latter height node is greater than that of the former height node, the corresponding height node is assigned to the candidate local deviation segment. When there are two consecutive height nodes where the node deviation is lower than the initial deviation boundary, the scanning stops, and the height node previously assigned to the candidate local deviation segment is used as the upper boundary node of the candidate local deviation segment.

4. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, A baseline for main slope deviation is generated based on the main slope data, and the over- or under-excavation deviation of the main slope is calculated, including: Read the signed deviations from the main slope data and organize them into a main deviation input sequence. Determine the main reference offset based on the median offset of the main deviation input sequence. The design slope surface in the slope design profile model is translated along the deviation calculation direction by the main reference offset to generate the main deviation reference. After registering the measured points in the main slope data, project them onto the main deviation benchmark to obtain the corresponding reference points. Based on the projection distance between the registered measured points and the corresponding reference points along the deviation calculation direction, the over- or under-excavation deviation of the main slope is formed.

5. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, Set contour sampling lines along the design slope toe line, including: If the design slope in the slope design profile model is a single slope, then each sampling reference point is expanded along the same design slope direction to form a profile sampling line. If the design slope in the slope design profile model consists of multiple slope segments, the corresponding slope segment is determined according to the design slope toe line segment where the sampling reference point is located, and the profile sampling line is formed by unfolding along the design slope direction of this corresponding slope segment.

6. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, Set contour sampling lines along the design slope toe line, including: The sampling reference point interval is determined based on the projection scalar sequence of the registered measured contour point set in the direction of the design slope toe line; If the projected scalar sequence can form an adjacent difference sequence, then the median of the adjacent difference sequence is taken as the sampling reference point interval. If the projected scalar sequence cannot form an adjacent difference sequence, the design section spacing in the slope design profile model shall be used as the sampling reference point interval. If the slope design profile model does not provide the design section spacing, the hole spacing in the blasting construction parameters shall be used as the sampling reference point interval.

7. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, The calculation methods also include: When the vertical distances from the same registered measured point to two adjacent contour sampling lines are equal, the corresponding registered measured point will be assigned to the contour sampling line with the smaller sampling reference point number. When there are fewer than 2 extractable height nodes on either side of the upper boundary of the slope toe hole bottom influence zone, the corresponding contour sampling line is marked as the boundary line to be maintained, and the boundary assignment determined in step S3 is maintained. When the boundary assignment changes, the main deviation benchmark is regenerated based on the verified main slope data; when the boundary assignment does not change, the main deviation benchmark formed in step S4 is retained.

8. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, The results of the slope blasting over- and under-excavation profile deviation include the profile sampling line assignment, height node order, deviation sign, and upper boundary node after verification, corresponding to each registered measured point; Among them, the contour sampling line belongs to the contour sampling line to which the corresponding registered measured point belongs, the height node order indicates the arrangement position of the corresponding registered measured point along the design slope direction from the slope foot to the upper part of the slope, and the deviation sign indicates the over-excavation and under-excavation direction of the corresponding registered measured point relative to the reference object.

9. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, The results of over- and under-excavation profile deviation in slope blasting include zoning deviation results; In the zoning deviation results, the over- or under-excavation deviation of the main slope surface corresponds to the main deviation benchmark, and the local over- or under-excavation deviation of the slope toe corresponds to the local reference point of the slope toe determined by the design slope toe line, the corresponding contour sampling line and the design slope direction. For the deviation results of the upper boundary node of the influence zone near the bottom of the hole at the slope toe, the partitioned deviation results retain the height node order and adjacency relationship.

10. The method for calculating the over- and under-excavation profile deviation in slope blasting according to claim 1, characterized in that, The boundary assignment is verified based on the deviation trend on both sides of the upper boundary of the influence zone at the toe of the slope hole, including: Based on the deviation sign, deviation amplitude proximity, and change direction connection relationship on both sides of the upper boundary of the slope toe hole bottom influence zone, generate one of the boundary stability mark and boundary continuation mark; When generating boundary stability markers, the upper boundary nodes of the slope toe hole bottom influence zone remain unchanged; When generating boundary continuation markers, the continuous height nodes that still appear as local deviations at the toe of the slope will be included in the influence zone at the bottom of the hole after verification, and the updated partition deviation results will be output.

Images