Mine drilling group track intelligent comparative analysis method based on multi-source data fusion

Abstract

The invention discloses a mine drilling group track intelligent comparative analysis method based on multi-source data fusion, and belongs to the technical field of mine engineering digitization. The method comprises the following steps: acquiring a drilling group track and performing track data analysis to obtain a track inclination angle, and numbering the track based on the size of the inclination angle to obtain a design track; actual drilling track data are obtained and preprocessed, the inclination angle average value of local and global actual drilling tracks is obtained, and an inclination angle matching algorithm is used for matching with a design track; and a dynamic coordinate translation mode is used for achieving alignment of the actual drilling track and the design track, and a final result is obtained. Through fusion processing, the method is suitable for analyzing drilling trajectories of different data formats, and the application range of the system is widened; automatic matching of a measurement track and a design track is achieved through the inclination angle, manual alignment errors are avoided, and the comparison precision is improved. Manual operation is reduced, and the automation level of engineering data management is improved.

Description

Technical Field

[0001] The present invention belongs to the technical field of mine engineering digitization, and in particular relates to an intelligent comparative analysis method for mine drilling group trajectories based on multi-source data fusion. Background Art

[0002] Drilling trajectory accuracy is a crucial factor in mine safety and mineral resource development. However, due to factors such as geological conditions and drilling techniques, the actual drilling trajectory often deviates from the designed trajectory, particularly in coal mines. Therefore, analyzing the designed and measured trajectory is crucial for ensuring construction quality, improving efficiency, and ensuring safe production during drilling operations. This provides a strong foundation for ensuring efficient and safe mine production.

[0003] However, there are already some methods for comparing and analyzing design and measured trajectories. However, design trajectories are usually presented in CAD drawings, while measured trajectory data are stored in formats such as Excel and CVS. The data formats differ greatly, making direct comparison and analysis difficult. In addition, the design drawings contain not only drilling trajectories, but also other line segment annotations, making it difficult to directly extract drilling trajectory line segments from the design drawings for comparison. At the same time, the coordinate systems of the design trajectory and the measured trajectory are inconsistent, requiring complex geometric calculations to achieve spatial alignment. Manual operations are prone to errors and inefficient. In addition, the comparative analysis results are often presented in the form of two-dimensional charts, lacking intuitive three-dimensional visualization effects, and the measured trajectory data cannot be directly written into the original design drawings. There is a lack of batch automated matching and error analysis functions, making it difficult to meet the processing needs of large-scale drilling trajectory data. Summary of the Invention

[0004] In order to solve the above technical problems, the present invention provides an intelligent comparative analysis method for mine drilling group trajectories based on multi-source data fusion.

[0005] To implement the above technology, the following steps are included: S1. Obtain the borehole group trajectory based on the design drawing data, perform trajectory data analysis on the borehole group trajectory, and obtain the spatial coordinates of the trajectory; obtain the trajectory inclination angle based on the spatial coordinates of the trajectory, number the trajectory based on the magnitude of the inclination angle, construct a structured format based on the numbering, and perform storage operations to obtain the designed trajectory; the steps are as follows: S1.1. Based on the CAD design drawing in DXF format, obtain the drilling group trajectory segments and trench line segments; Parsing CAD drawings in DXF format to analyze the cluttered lines of the design. The analysis method is as follows: based on the fact that all drilling trajectories intersect with the trenches, and their extension lines intersect at one point, the drilling trajectory segments and trench segments can be determined based on their topological relationship. S1.2. Define the intersection of the drilling trajectory segment and the trench segment as the starting point, and the end point of the drilling trajectory segment as the final point; and convert the two-dimensional design drawing into a three-dimensional drawing to facilitate subsequent comparison and analysis; The way to convert a two-dimensional design drawing into a three-dimensional drawing is: convert the two-dimensional design drawing into a three-dimensional drawing with the X coordinate of 0, that is, the starting point of the trajectory is represented as , the end point coordinates are expressed as ; S1.3. Calculate the inclination angle of the trajectory based on the three-dimensional graph; The calculation method is: based on the straight line segment obtained from the drilling group trajectory, the line segment is projected onto the YZ plane in the three-dimensional graph, and the inclination angle on the YZ plane is calculated. , the expression is as follows: ; S1.4. Calculate the length of the trajectory based on the obtained starting and final points , the expression is as follows: ; S1.5, numbering the drilling trajectories in descending order based on the obtained inclination angles, constructing a structured format after numbering, and performing a storage operation; The structured format includes: the number of the drilling trajectory, the coordinates of the starting point of the drilling trajectory, the coordinates of the final point of the drilling trajectory, the inclination angle of the drilling trajectory and the hole depth.

[0006] S2. Acquire actual drilling trajectory data based on open source data, perform preprocessing on the actual drilling trajectory data, and calculate the average inclination angle of the local and global actual drilling trajectory based on the preprocessed data; Here are the steps: S2.1. Obtain measured drilling trajectory data in multiple sheets in Excel format from open source data; S2.2. Preprocessing the actual drilling trajectory data, including: S2.2.1. Calculate the rate of change of the plane slope between adjacent coordinate points. If it exceeds a set threshold, it is determined to be an outlier and removed. The set threshold is 30%. S2.2.2. For the missing coordinate points after elimination, use the cubic spline interpolation method to fill them in and generate a continuous measurement trajectory curve; S2.3. Calculate the average inclination angle of the local and global actual drilling trajectories for the pre-processed data; wherein the local actual drilling trajectory is defined as: ; The expression for the average local tilt angle is as follows: Where, represents the average local tilt angle; Indicates the number of continuous line segments selected when calculating the local inclination (e.g. k=10 means taking line segments between 10 adjacent points); Indicates the starting point index of the local calculation (e.g. m=10 means calculation starts from the 10th point); Indicates the current calculation Segment index; and Respectively represent the Z coordinate and Y coordinate of the midpoint of the trajectory; Indicates the indexes; The expression of the tilt average of the global measurement trajectory is as follows: Where, represents the global average tilt angle; represents the number of bus segments of the trajectory, ∈ ; S2.4. Calculate the hole depth based on the pre-processed data. The expression is as follows: Where, The X, Y, and Z coordinates of the final point; the opening point of each trajectory is used as the starting point coordinate , both are (0,0,0).

[0007] S3, based on the obtained average inclination angles of the local and global actual drilling trajectories, matching them with the designed trajectory using an inclination matching algorithm; The inclination matching algorithm is: calculate the inclination of the designed trajectory and the local average inclination The difference between the two is used to find the design trajectory with the smallest inclination difference for matching; if an actual drilling trajectory matches multiple design trajectories, the inclination angle is calculated and the average inclination angle of the entire measurement trajectory is calculated. The difference is used to determine the optimal match.

[0008] S4. Based on the matching results, the actual drilling trajectory is aligned with the designed trajectory using dynamic coordinate translation to obtain the final result and complete the mine drilling comparison analysis; The alignment formula is as follows: Where, Represents all coordinate point data of the actual trajectory; Indicates the coordinate point after alignment.

[0009] S5. Use the error analysis method to verify the inclination angle and hole depth errors between the designed trajectory and the actual trajectory.

[0010] Beneficial effects of the present invention: The present invention integrates CAD design data with Excel measured data, which is applicable to drilling trajectory analysis requirements in different data formats and improves the applicability of the system. It uses the inclination angle to achieve automatic matching of the measured trajectory and the designed trajectory, avoiding manual alignment errors and improving comparison accuracy.

[0011] The present invention adopts dynamic coordinates to automatically and accurately align the measured trajectory with the designed trajectory in batches.

[0012] The present invention is expected to display trajectory comparison diagrams and error information, which improves the intuitiveness of analysis and decision-making efficiency; and the present invention supports directly writing the measurement trajectory after translation and alignment back into the CAD design drawing, reducing manual operations and improving the automation level of engineering data management. BRIEF DESCRIPTION OF THE DRAWINGS

[0013] Figure 1 It is a schematic diagram of the process of the present invention; Figure 2 is the error analysis result of the embodiment of the present invention; Figure 3 This is a 3D trajectory comparison result diagram of an embodiment of the present invention; Figure 4 2 is a diagram showing the comparison results of the YZ profile trajectory of an embodiment of the present invention. DETAILED DESCRIPTION

[0014] The present invention is further described in detail below with reference to specific embodiments.

[0015] like Figure 1 As shown, a method for intelligent comparative analysis of mine drilling group trajectories based on multi-source data fusion includes the following steps: S1. Obtain the borehole group trajectory based on the design drawing data, perform trajectory data analysis on the borehole group trajectory, and obtain the spatial coordinates of the trajectory; obtain the trajectory inclination angle based on the spatial coordinates of the trajectory, number the trajectory based on the magnitude of the inclination angle, construct a structured format based on the numbering, and perform storage operations to obtain the designed trajectory; the steps are as follows: S1.1. Based on the CAD design drawing in DXF format, obtain the drilling group trajectory segments and trench line segments; Parsing CAD drawings in DXF format to analyze the cluttered lines of the design. The analysis method is as follows: based on the fact that all drilling trajectories intersect with the trenches, and their extension lines intersect at one point, the drilling trajectory segments and trench segments can be determined based on their topological relationship. S1.2. Define the intersection of the drilling trajectory segment and the trench segment as the starting point, and the end point of the drilling trajectory segment as the final point; and convert the two-dimensional design drawing into a three-dimensional drawing to facilitate subsequent comparison and analysis; The coordinates of the starting point of the drilling design trajectory are not the starting point of the trajectory segment, but the intersection of the segment and the trench segment. The intersection is the drilling point. The way to convert a two-dimensional design drawing into a three-dimensional drawing is: convert the two-dimensional design drawing into a three-dimensional drawing with the X coordinate of 0, that is, the starting point of the trajectory is represented as , the end point coordinates are expressed as ; S1.3. Calculate the inclination angle of the trajectory based on the three-dimensional graph; The calculation method is: based on the straight line segment obtained from the drilling group trajectory, the line segment is projected onto the YZ plane in the three-dimensional graph, and the inclination angle on the YZ plane is calculated. , the expression is as follows: ; S1.4. Calculate the length of the trajectory based on the obtained starting and final points , the expression is as follows: ; S1.5, numbering the drilling trajectories in descending order based on the obtained inclination angles, constructing a structured format after numbering, and performing a storage operation; The structured format includes: the number of the drilling trajectory, the coordinates of the starting point of the drilling trajectory, the coordinates of the final point of the drilling trajectory, the inclination angle of the drilling trajectory and the hole depth.

[0016] S2. Acquire actual drilling trajectory data based on open source data, perform preprocessing on the actual drilling trajectory data, and calculate the average inclination angle of the local and global actual drilling trajectory based on the preprocessed data; Here are the steps: S2.1. Obtain measured drilling trajectory data in multiple sheets in Excel format from open source data; S2.2. Preprocessing the actual drilling trajectory data, including: S2.2.1. Calculate the rate of change of the plane slope between adjacent coordinate points. If it exceeds a set threshold, it is determined to be an outlier and removed. The set threshold is 30%. S2.2.2. For the missing coordinate points after elimination, use the cubic spline interpolation method to fill them in and generate a continuous measurement trajectory curve; S2.3. Calculate the average inclination angle of the local and global actual drilling trajectories for the pre-processed data; wherein the local actual drilling trajectory is defined as: ; The expression for the average local tilt angle is as follows: Where, represents the average local tilt angle; Indicates the number of continuous line segments selected when calculating the local inclination (e.g. k=10 means taking line segments between 10 adjacent points); Indicates the starting point index of the local calculation (e.g. m=10 means calculation starts from the 10th point); Indicates the current calculation Segment index; and Respectively represent the Z coordinate and Y coordinate of the midpoint of the trajectory; Indicates the indexes; The expression of the tilt average of the global measurement trajectory is as follows: Where, represents the global average tilt angle; represents the number of bus segments of the trajectory, ∈ ; S2.4. Calculate the hole depth based on the pre-processed data. The expression is as follows: Where, The X, Y, and Z coordinates of the final point; the opening point of each trajectory is used as the starting point coordinate , both are (0,0,0).

[0017] S3, based on the obtained average inclination angles of the local and global actual drilling trajectories, matching them with the designed trajectory using an inclination matching algorithm; The inclination matching algorithm matches the actual drilling trajectory with the designed trajectory by inclination similarity, giving priority to local matching and then global matching; The inclination matching algorithm is: calculate the inclination of the designed trajectory and the local average inclination The difference between the two is used to find the design trajectory with the smallest inclination difference for matching; if an actual drilling trajectory matches multiple design trajectories, the inclination angle is calculated and the average inclination angle of the entire measurement trajectory is calculated. The difference is used to determine the optimal match.

[0018] S4. Based on the matching results, the actual drilling trajectory is aligned with the designed trajectory using dynamic coordinate translation to obtain the final result and complete the mine drilling comparison analysis; In trajectory alignment, dynamic coordinate translation is used to achieve precise alignment between the measured trajectory and the designed trajectory. Specifically, the starting coordinates of the designed trajectory of the drilling group trajectory are different and irregular for each trajectory, while the starting coordinates of the measured trajectory are all (0,0,0), making it difficult to complete the translation alignment of the drilling group trajectory using a single reference point. Therefore, it is proposed to achieve precise alignment between the measured trajectory and the designed trajectory through dynamic coordinate translation. The dynamic coordinates are the starting coordinates of each designed trajectory. According to the matching results, all coordinate point data of each measured trajectory are added to the matching starting coordinates of the designed trajectory to perform trajectory translation, achieving precise alignment. The alignment formula is as follows: Where, Represents all coordinate point data of the actual trajectory; Indicates the coordinate point after alignment.

[0019] S5. Use error analysis method to verify the inclination angle and hole depth errors between the designed trajectory and the actual trajectory; like Figure 2 As shown, in this embodiment, 14 error analyses were performed. From the inclination error, it can be seen that the inclination error range is: -2.83°~5.13°; 10 / 14 of the inclination errors are less than 2°, indicating that the drilling direction control is good and meets the design requirements. In addition, except for design trajectory 1, the inclination errors of 13 trajectories are less than 3°, indicating that the designed drilling direction distribution is relatively uniform; Among them, the reason for the large inclination deviation of design trajectory 1 is due to complex geological conditions or equipment offset; It can be seen from the hole depth error that the hole depth error distribution is relatively uniform, the number of positive and negative errors is close (7 positive vs. 7 negative), and the hole depth error range is: -1.22 meters to 1.92 meters, indicating that the present invention has good control over the hole depth error and meets the design requirements.

[0020] like Figure 3 and Figure 4 As shown, this embodiment shows the comparison diagram of the 3D and YZ cross-sections of the designed trajectory and the measured trajectory and the error information; the present invention can also support writing the measured trajectory after translation alignment into the original CAD design drawing with one click.

[0021] It should be noted that the above are only preferred embodiments of the present application and do not limit the scope of patent protection of the present application. Any equivalent structure or equivalent process transformation made using the contents of the description and drawings of this application, or directly or indirectly applied in other related technical fields, are also included in the scope of patent protection of the present application.

Claims

1. An intelligent comparative analysis method for mine drilling group trajectories based on multi-source data fusion, characterized in that: The following steps are involved: S1. Obtain the trajectory of the drilling group according to the design data, perform trajectory data analysis on the trajectory of the drilling group, and obtain the spatial coordinates of the trajectory; The inclination angle of the trajectory is obtained through the spatial coordinates of the trajectory, the trajectory is numbered based on the size of the inclination angle, a structured format is constructed based on the numbering and a storage operation is performed to obtain the designed trajectory; S2. Acquire actual drilling trajectory data based on open source data, perform preprocessing on the actual drilling trajectory data, and calculate the average inclination angle of the local and global actual drilling trajectory based on the preprocessed data; S3, based on the obtained average inclination angles of the local and global actual drilling trajectories, matching them with the designed trajectory using an inclination matching algorithm; S4. Based on the matching results, the actual drilling trajectory is aligned with the designed trajectory using dynamic coordinate translation to obtain the final result and complete the mine drilling comparison analysis; S5. Use the error analysis method to verify the inclination angle and hole depth errors between the designed trajectory and the actual trajectory.

2. The method for intelligent comparative analysis of mine drilling group trajectories based on multi-source data fusion according to claim 1, characterized in that: The steps of obtaining a borehole group trajectory according to the design drawing data, parsing the trajectory data through the borehole group trajectory to obtain the spatial coordinates of the trajectory; obtaining the inclination angle of the trajectory through the spatial coordinates of the trajectory, numbering the trajectory based on the magnitude of the inclination angle, constructing a structured format based on the numbering, and performing a storage operation, to obtain the design trajectory include: S1.

1. Based on the CAD design drawing in DXF format, obtain the drilling group trajectory segments and trench line segments; S1.

2. Define the intersection of the drilling trajectory segment and the trench segment as the starting point, and the end point of the drilling trajectory segment as the final point; Converting a two-dimensional design drawing into a three-dimensional drawing; converting a two-dimensional design drawing into a three-dimensional drawing by: converting the two-dimensional design drawing into a three-dimensional drawing with all X coordinates being 0; S1.

3. Calculate the inclination angle of the trajectory based on the three-dimensional graph; The calculation method is as follows: based on the straight line segments obtained from the drilling group trajectory, the line segments are projected onto the YZ plane in the three-dimensional graph, and the inclination angle on the YZ plane is calculated; S1.

4. Calculate the length of the trajectory based on the obtained starting and final points; S1.5, numbering the drilling trajectories in descending order based on the obtained inclination angles, constructing a structured format after numbering, and performing a storage operation; The structured format includes: the number of the drilling trajectory, the coordinates of the starting point of the drilling trajectory, the coordinates of the final point of the drilling trajectory, the inclination angle of the drilling trajectory and the hole depth.

3. The method for intelligent comparative analysis of mine drilling group trajectories based on multi-source data fusion according to claim 1, characterized in that: The steps of obtaining actual drilling trajectory data according to open source data, performing preprocessing on the actual drilling trajectory data, and calculating the average values ​​of the inclination angles of the local and global actual drilling trajectories based on the preprocessed data are as follows: S2.

1. Obtain measured drilling trajectory data from open source data; S2.

2. Preprocessing the actual drilling trajectory data, including: S2.2.

1. Calculate the rate of change of the plane slope between adjacent coordinate points. If it exceeds the set threshold, it is determined to be an outlier and removed. S2.2.

2. For the missing coordinate points after elimination, use the cubic spline interpolation method to fill them in and generate a continuous measurement trajectory curve; S2.3, calculating the average value of the inclination angles of the local and global actual drilling trajectories for the pre-processed data; S2.

4. Calculate the hole depth of the pre-processed data.

4. The method for intelligent comparative analysis of mine drilling group trajectories based on multi-source data fusion according to claim 1, characterized in that: The obtained average values ​​of the inclination angles of the local and global actual drilling trajectories are matched with the designed trajectory using an inclination matching algorithm; The inclination matching algorithm is as follows: calculate the difference between the inclination of the designed trajectory and the average local inclination angle, and find the designed trajectory with the smallest inclination difference for matching; If an actual drilling trajectory matches multiple designed trajectories, the difference between the inclination angle and the average inclination angle of the entire measured trajectory is calculated to determine the optimal match.

5. The method for intelligent comparative analysis of mine drilling group trajectories based on multi-source data fusion according to claim 1, characterized in that: Based on the matching results, the dynamic coordinate translation method is used to achieve the alignment of the actual drilling trajectory and the designed trajectory to obtain the final result and complete the mine drilling comparison analysis. The dynamic coordinate translation method is: all the coordinate point data of the actual trajectory are added to the matching starting point coordinates of the designed trajectory to translate the trajectory to achieve alignment.

Images