A method and system for generating a three-dimensional spatial explosibility characteristic cloud chart of rock

Abstract

The application provides a rock three-dimensional space explosibility characteristic cloud chart generation method and system, and belongs to the technical field of blasting engineering. The method comprises the following steps: drilling operation is performed by using a rock drilling jumbo, a drilling parameter is obtained once every drilling depth, and the drilling parameter is stored in a database as a vector; original data is calculated, hole data is converted into three-dimensional space coordinates, and the three-dimensional space coordinates are stored in the database; rock compressive strength of each measuring point is calculated according to the relationship between drilling specific energy and rock strength; corresponding rock tensile strength is further obtained and stored in the database; a new vector is subjected to cyclic calculation of cross section acquisition and dichotomy interpolation, and n complete two-dimensional cross section vector data are obtained; a scatter diagram is drawn, and an excavation shape is obtained; a two-dimensional partition map is obtained by filling through Kriging interpolation, and three-dimensional vector data are cut into n two-dimensional explosibility thermal maps. The application automatically processes collected data to obtain an explosibility result, and improves engineering efficiency.

Description

Technical Field

[0001] This invention relates to the field of blasting engineering technology, and in particular to a method and system for generating a three-dimensional spatial blastability characteristic cloud map of rock. Background Technology

[0002] Tunnel blasting is a type of blasting engineering technology. It is a highly efficient and adaptable underground engineering construction method, widely used in projects traversing complex geological conditions. Its main characteristics include rapid breaking of large volumes of rock, significantly improved construction efficiency, reduced labor intensity, and good economic benefits. Modern blasting technology, through precise design and control, can effectively reduce the impact on the surrounding environment and ensure construction safety.

[0003] With the continuous development of blasting engineering technology, traditional drilling techniques, which require manual measurement and data collection, are not only time-consuming but also prone to errors. Currently, intelligent data acquisition systems using sensors are gradually emerging to collect and transmit parameters of the drilling rig in real time, such as drilling speed, drilling pressure, average rotational pressure, water pressure, water flow, and depth. Furthermore, the relative position of the drilled hole can be set in tunnel drilling, allowing for real-time monitoring of the tunnel's progress. This technology significantly improves drilling efficiency and reduces errors caused by data collection and analysis. However, some problems remain in data analysis. The amount of data collected and transmitted by the drilling rig in each tunnel and even each project is enormous. Data analysis can lead to significant time waste, potentially delaying the next drilling or blasting operation and causing substantial losses in time, personnel, and even finances. Summary of the Invention

[0004] This invention provides a method and system for generating three-dimensional spatial explosiveness feature cloud maps of rocks, which solves the problem of the difficulty in efficiently and effectively processing the data collected by rock drilling rigs in the prior art.

[0005] In a first aspect, the present invention provides a method for generating a three-dimensional spatial explosiveness characteristic cloud map of rock, comprising: using a rock drilling rig to perform drilling operations, acquiring drilling parameters as a vector and storing them in a database every time the drilling depth is d; the drilling parameters include: measuring point number, borehole data, and drilling measurement data, wherein the borehole data includes the three-dimensional spatial coordinates of the borehole opening corresponding to the measuring point number and the three-dimensional azimuth angle of the drilling direction, and the drilling measurement data includes various measurement parameters corresponding to the measuring point number; and using a spatial coordinate transformation formula, transforming the borehole data in the database to obtain the three-dimensional spatial data of all boreholes in the working direction perpendicular to the excavation face. Spatial vector; based on the relationship between drilling specific energy and rock strength, the rock tensile strength of each three-dimensional spatial coordinate equidistantly arranged along the three-dimensional spatial vector of the borehole is obtained to replace the drilling measurement data in the database and form a new vector in the database; the obtained new vector is divided into multiple two-dimensional sections equidistantly along the working direction, and a scatter plot is drawn based on the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional section and the rock tensile strength value; the rock tensile strength value of each scatter point in the scatter plot is filled by Kriging interpolation to obtain a two-dimensional partition map; all two-dimensional partition maps are used to draw a three-dimensional spatial blastability characteristic cloud map of the rock.

[0006] According to the method for generating a three-dimensional spatial explosiveness characteristic cloud map of rocks provided by the present invention, the spatial coordinate transformation formula includes:

[0007] Xi=X0+L×cos(α)×cos(β);

[0008] Yi=Y0+L×cos(β)×sin(α);

[0009] Zi = Z0 - L × sin(β);

[0010] In the formula, X0, Y0, Z0 represent the three-dimensional spatial coordinates of the borehole opening, Xi, Yi, Zi represent the three-dimensional spatial coordinates of the i-th drilling depth d, L represents the cumulative depth of the borehole after the i-th drilling depth d, and α, β and γ represent the Euler angles of the borehole offset in the X, Y and Z directions, respectively.

[0011] According to the method for generating a three-dimensional explosiveness characteristic cloud map of rock provided by the present invention, the tensile strength of rock is obtained by calculating the rock tensile strength of each three-dimensional spatial coordinate equidistantly arranged along the three-dimensional spatial vector of the borehole based on the relationship between drilling specific energy and rock strength. This includes: calculating the drilling specific energy and impact energy required for the rock drilling rig to drill a unit volume of rock through rock breaking mechanics analysis of various drill bits; calculating the rock breaking specific energy based on the relevant conversion relationship between impact energy and breaking specific energy; and calculating the rock compressive strength at the first to i-th drilling depths d of each measuring point based on the relationship between breaking specific energy and rock saturated uniaxial compressive strength.

[0012] According to the method for generating a three-dimensional explosiveness characteristic cloud map of rock provided by the present invention, the drilling specific energy and impact energy required for a rock drilling rig to drill a unit volume of rock are calculated through rock breaking mechanics analysis of various drill bits, including:

[0013]

[0014]

[0015] Where e represents the drilling specific energy, e i Impact energy, Pi represents average impact pressure, Pt represents average propulsion pressure, Pn represents average rotational pressure, and V represents average propulsion speed.

[0016] According to the method for generating a three-dimensional explosiveness feature cloud map of rocks provided by the present invention, the obtained new vector is divided into multiple two-dimensional sections at equal intervals along the working direction, including: when there are missing data in the two-dimensional sections, the missing data is filled in to obtain multiple complete two-dimensional sections.

[0017] The method for generating a three-dimensional explosiveness feature cloud map of rock according to the present invention, which plots a scatter plot based on the three-dimensional spatial coordinates of the intersection points of each borehole at a two-dimensional cross-section and the tensile strength value of the rock, includes: determining the position of the scatter points based on the three-dimensional spatial coordinates of the intersection points of each borehole at a two-dimensional cross-section; and determining the color value of the scatter points based on the tensile strength value of the rock; wherein the color value and the tensile strength value of the rock have a linear relationship.

[0018] According to the method for generating a three-dimensional spatial explosiveness feature cloud map of rocks provided by the present invention, the color value is a grayscale value.

[0019] Secondly, the present invention also provides a system for generating a three-dimensional spatial explosiveness feature cloud map of rocks, comprising:

[0020] The first processing module is used to perform drilling operations using a rock drilling rig. Every time the drilling depth is d, the drilling parameters are acquired once and stored as a vector in the database. The drilling parameters include: measuring point number, hole data and drilling measurement data. The hole data includes the three-dimensional spatial coordinates of the borehole opening corresponding to the measuring point number and the three-dimensional azimuth angle of the drilling direction. The drilling measurement data includes various measurement parameters corresponding to the measuring point number.

[0021] The second processing module is used to transform the hole data in the database using the spatial coordinate transformation formula to obtain the three-dimensional spatial vector of all boreholes in the working direction perpendicular to the excavation face.

[0022] The third processing module is used to calculate the rock tensile strength of each three-dimensional spatial coordinate equidistantly arranged along the three-dimensional spatial vector of the borehole based on the relationship between drilling specific energy and rock strength, so as to replace the drilling measurement data in the database and form a new vector in the database.

[0023] The fourth processing module is used to divide the new vector obtained into multiple two-dimensional sections at equal intervals along the working direction, and to draw a scatter plot based on the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional sections and the tensile strength value of the rock.

[0024] The fifth processing module is used to fill in the rock tensile strength value of each scatter point in the scatter plot using Kriging interpolation to obtain a two-dimensional partition map;

[0025] The sixth processing module is used to draw a three-dimensional spatial explosiveness characteristic cloud map of the rock using all the two-dimensional partition maps.

[0026] Thirdly, the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of the method for generating a three-dimensional spatial explosiveness feature cloud map of rock as described above.

[0027] Fourthly, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method for generating a three-dimensional spatial explosiveness feature cloud map of rock as described above.

[0028] The method and system for generating three-dimensional spatial explosiveness characteristic cloud maps of rocks provided by this invention can draw three-dimensional spatial explosiveness characteristic cloud maps of rocks based on drilling rig parameters and rock tensile strength. During drilling, it effectively predicts the rock's solidity and strength and visualizes its explosiveness. This has a significant practical impact on improving construction efficiency.

[0029] This invention can obtain the parameters and cloud map after coordinate transformation in real time. Engineers can monitor the rock drilling rig during operation, view and analyze the data, identify and resolve problems in a timely manner, and improve the efficiency of charge calculation in the blasting process. Furthermore, it reduces the personnel, time and economic requirements in the drilling and blasting process, thereby improving engineering efficiency. Attached Figure Description

[0030] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0031] Figure 1 This is a flowchart illustrating the method for generating a three-dimensional spatial explosiveness feature cloud map of rocks provided by the present invention.

[0032] Figure 2 This is a two-dimensional scatter plot with missing numbers provided by the present invention;

[0033] Figure 3 This is the complete two-dimensional scatter plot provided by the present invention;

[0034] Figure 4 This is the explosive zone division diagram provided by the present invention;

[0035] Figure 5 This is a schematic diagram of the structure of the electronic device provided by the present invention. Detailed Implementation

[0036] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0037] It should be noted that, in the description of the embodiments of the present invention, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Those skilled in the art can understand the specific meaning of the above terms in the present invention according to the specific circumstances.

[0038] The terms "first," "second," etc., used in this application are used to distinguish similar objects and not to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that embodiments of this application can be implemented in orders other than those illustrated or described herein, and the objects distinguished by "first," "second," etc., are generally of the same class and the number of objects is not limited; for example, a first object can be one or more.

[0039] The following is combined with Figures 1-5 This invention describes a method and system for generating three-dimensional spatial explosiveness feature cloud maps of rocks, as provided in embodiments of the present invention.

[0040] Figure 1 This is a flowchart illustrating the method for generating a three-dimensional spatial explosiveness feature cloud map of rocks provided by the present invention, as shown below. Figure 1 As shown, including but not limited to the following steps:

[0041] S1. Drilling operations are performed on the rock using a drilling rig. Bever Team Online software is used to analyze and summarize the drilling parameters of the drilling rig. The drilling rig acquires drilling parameters once per drilling depth (d), including: measuring point number, hole data, and drilling measurement data. The drilling measurement data includes the borehole depth (m), average advance speed (m / min), average impact pressure (bar), average advance pressure (bar), average rotational pressure (bar), average water pressure (bar), and average water flow rate (L / min) corresponding to the measuring point number. 1 bar = 100 kPa. The hole data includes the three-dimensional spatial coordinates (X0, Y0, Z0) of the borehole opening corresponding to the measuring point number, and the three-dimensional azimuth angles (α, β, γ) of the drilling direction.

[0042] S2. The borehole data k corresponding to each measuring point number Q and the drilling measurement data Data (depth L, propulsion speed V, impact pressure Pi, propulsion pressure Pt, rotation pressure Pn, water pressure, water flow) at the measuring point are stored in the database as a vector A{Q, k, Data}.

[0043] S3. Using the spatial coordinate transformation formula, the three-dimensional spatial coordinates, three-dimensional azimuth angle, and depth of each measuring point are transformed according to formulas (1) to (3) to obtain the three-dimensional spatial coordinates (Xi, Yi, Zi) of the i-th drilling depth d at each measuring point. The three-dimensional spatial coordinates (X0, Y0, Z0) of vector A are updated to i sets of three-dimensional spatial coordinates as vector A of the i-th drilling depth d at each measuring point. di Save it to the database.

[0044] The borehole data extracted directly from the rock drilling rig, including the borehole depth, three-dimensional spatial coordinates, and three-dimensional azimuth, needs to be converted using a special algorithm into the format used in the method for drawing three-dimensional blastability feature cloud maps. The spatial coordinate conversion method used in this method is as follows:

[0045] Xi=X0+L×cos(α)×cos(β) Formula (1)

[0046] Yi=Y0+L×cos(β)×sin(α) Formula (2)

[0047] Zi=Z0-L×sin(β) Equation (3)

[0048] In the formula, X0, Y0, Z0 represent the three-dimensional spatial coordinates of the borehole opening, Xi, Yi, Zi represent the three-dimensional spatial coordinates of the i-th drilling depth d, L represents the cumulative depth (m) of the borehole after the i-th drilling depth d, and α, β, and γ represent the Euler angles of the borehole offset in the X, Y, and Z directions, respectively. Since there was no offset or rotation in the Z direction during drilling, this spatial transformation method does not use the γ Euler angle.

[0049] α: An angle of rotation around the Z-axis, used to represent the changes in borehole diameter in the X and Y directions.

[0050] β: An angle of rotation around the Y-axis, used to represent the vertical offset of the borehole in the Z-direction.

[0051] γ: Angle of rotation around the X-axis.

[0052] The Y-axis is parallel to the horizontal plane and the excavation face;

[0053] The Z-axis is perpendicular to the horizontal plane;

[0054] The X-axis is perpendicular to the excavation face.

[0055] S4. Process the measurement-while-drilling data and calculate the rock compressive strength at the first to i-th drilling depths d for each measuring point based on the relationship between drilling specific energy and rock strength.

[0056] Through rock-breaking mechanics analysis of three types of drill bits, the energy required for a rock drilling rig to drill a unit volume of rock, i.e., the drilling specific energy e (MJ / m³), was calculated. 3 ) and impact energy e i ;

[0057] The formula is as follows:

[0058]

[0059]

[0060] Where: Pi—average impact pressure (bar), Pt—average propulsion pressure (bar), Pn—average rotation pressure, V—average propulsion speed.

[0061] By collecting drilling parameters from the rock drilling rig and the drilling machine itself, the rock-breaking energy of each component can be calculated. Ultimately, the drilling-breaking specific energy and the percentage of energy allocated to each component can be calculated. The impact of the rock drilling rig can generate stress concentration, and this energy accounts for an average of 88% during the rock-breaking process; this is the impact energy e of the drilling specific energy. i It plays a dominant role.

[0062] Furthermore, based on the correlation between the impact energy ei, which plays a dominant role in the drilling specific energy, and the breaking specific work a; a = 0.86e i It can calculate the specific work required to break rocks.

[0063] In the formula: e i —Impact energy (MJ / m 3 a—Specific work of crushing (MJ / m³) 3 ).

[0064] A large amount of data was collected on the specific work 'a' of rock crushing and the saturated uniaxial compressive strength 'δ' of rock. c Based on the research data, fitted formulas for rock crushing specific work and rock saturated uniaxial compressive strength were obtained: δ c =0.205α;

[0065] Where: δ i —Saturated uniaxial compressive strength of rock (MPa); a—Specific work of crushing (J / cm²) 3 ).

[0066] Based on the relationship between the specific work of rock breaking and the saturated uniaxial compressive strength of rock, the rock compressive strength δ at the drilling depth d of the first to i drilling operations at each measuring point is calculated. c .

[0067] S5. Based on the ratio of the rock's compressive strength to its tensile strength, calculate the rock's tensile strength σ at the drilling depth d of the first to i drilling passes at each measuring point, and replace the vector A at the drilling depth d of the first to i drilling passes at the measuring point. di The measurement-while-drilling data (Data) is collected and saved together as a vector A. di σ{Q,k,σ} is stored in the database as an index characterizing the explosiveness of rocks;

[0068] S6. For the new vector, perform iterative calculations of cross-section acquisition and bisection interpolation to obtain n complete two-dimensional cross-section vector data.

[0069] Add the index Q of all measurement points to array N, and iterate through all A points based on Q. di σ, find vector A di The maximum and minimum X-coordinate values ​​of σ are determined by the fact that the rock drilling rig records data every 2 cm of straight drilling along the vertical excavation face. Therefore, the vector A is determined based on the range of X-coordinate values. di σ is divided into (Max-min) / 2 two-dimensional sections, and the X value of each section is x. i And each two-dimensional section is a two-dimensional vector C i {X i ,q j ,{Yi,Zi},σ i Each two-dimensional vector C iThere should be a measurement point number q for each measurement point number in array N. j Xi, Yi, Zi, σ i Value. If a communication failure causes a missed acquisition of borehole data and measurement-while-drilling data corresponding to the j-th measurement point number during the i-th acquisition, then the two-dimensional vector C... i The sequence number of the j-th measuring point q j If missing, the following algorithm should be used for interpolation to complete the missing parts:

[0070] For each two-dimensional vector C i Perform a loop to examine each two-dimensional vector C. i All q j Compared with Q, if any measurement point numbers are missing, the missing measurement point numbers are added to array k in the loop. i In the middle, record the X of the boreholes that do not have missing measurement point numbers. i The value is also equivalent to the X value of the borehole with missing measurement point number. i Value; based on k i The measurement point number in the three-dimensional vector A di Search for k in i Given all the measuring point serial numbers, the drilling measurement data and three-dimensional spatial coordinates are used to find the X-axis that matches the recorded X-axis in the drilling measurement data and three-dimensional spatial coordinates, provided that the measuring point numbers are the same. i The absolute value of the difference is less than Δd. This includes drilling measurement data and three-dimensional spatial coordinates (Δd = 2cm, drilling measurement data includes borehole advance speed V, impact pressure Pi, advance pressure Pt, rotational pressure Pn, water pressure, and water flow rate at the measurement point number), for example, {X... i-1 Y i-1 Z i-1 Data i-1} and {X i+1 Y i+1 Z i+1 Data i+1 The measurement point number of} and {X i Y i Z i Data i (The two sets of data are the same), using the two sets of data, the data corresponding to the missing measurement point number {X} is obtained by using the bisection method (i.e., adding each parameter of the two sets of data and averaging). i Y i Z i Data i}, then based on Data i Calculate the tensile strength σ of the rock according to steps S4-S5. i Complete the two-dimensional vector C in sequence i The data corresponding to each missing measurement point number in {X} i Yi Z i , σ i}

[0071] If the borehole has a missing measuring point number, X i If the maximum X-coordinate depth of the borehole is greater than the maximum value or the minimum X-coordinate depth of the borehole is less than the minimum value, then the measurement-while-drilling data and spatial three-dimensional coordinates corresponding to the maximum or minimum X-coordinate depth of the borehole are used as the three-dimensional vector A. di The missing borehole measurement data and three-dimensional spatial coordinates for the corresponding measurement point number are then obtained based on the Data. i Calculate the tensile strength σ of the rock according to steps S4-S5. i Complete the two-dimensional vector C i The data corresponding to each missing measurement point number in {X} i Y i Z i , σ i Since the drilling rig did not pass through this point during drilling, but in order to ensure the integrity of the two-dimensional diagram, the above method is used to fill it in (the data at this point is not accurate, and different colors can be used to distinguish it when drawing).

[0072] Repeat the above steps to obtain (Max-min) / 2 complete two-dimensional vectors.

[0073] S7. Draw a scatter plot of the two-dimensional vectors of (Max-min) / 2 complete two-dimensional sections to obtain the excavation geometry.

[0074] Partial data {X} from a complete two-dimensional vector i ,{Yi,Zi},σ i Create arrays separately, and use the scatter method in the matplotlib library of Python (the scatter method is a way to draw scatter plots in the matplotlib library, which can draw points based on x, y, z coordinates, and the color of the points is determined by the value of the tensile strength of the rock) to draw scatter plots (the horizontal axis of the scatter plot is the y value of the two-dimensional vector, and the vertical axis is the z value).

[0075] S8. A two-dimensional partition map is obtained by using the rock tensile strength of each point in the scatter plot to fill in the gaps using Kriging interpolation. (Kriging interpolation, also known as spatial local interpolation, mainly involves interpolating the values ​​of the space near each point using the rock tensile strength to ensure the integrity of the values ​​in the two-dimensional region and obtain a complete partition map.)

[0076] S9. Repeat steps S7-S8 to obtain the partition maps of all two-dimensional sections. Based on these partition maps of all two-dimensional sections, a three-dimensional spatial blastability characteristic cloud map of the excavation area can be drawn.

[0077] The method for drawing a three-dimensional spatial explosiveness characteristic cloud map of rock provided by this invention can generate such a map based on drilling rig parameters and rock tensile strength. During drilling, it effectively predicts the rock's firmness and strength and visualizes its explosiveness. This has a significant practical impact on improving construction efficiency.

[0078] To provide a clearer explanation of the technical solution of this invention, the following description is provided in conjunction with a specific engineering scenario.

[0079] Project Overview: The tunnel has a straight wall topped with a semi-circular arch. The tunnel strata are mainly granite. The construction method involves full-section, one-time excavation. The surrounding rock is relatively good. Excavation is carried out in the working direction at the working face, according to the railway department's design requirements. The working direction is perpendicular to the working face and points into the rock.

[0080] The surrounding rock conditions of the tunnel are relatively good. The tunnel is accessed from the entrance and exit by drilling rigs. The full-section one-time drilling and blasting excavation scheme is adopted, that is, the excavation is carried out by setting the central slot hole, the surrounding auxiliary holes, and the peripheral holes according to the design outline for smooth blasting.

[0081] Feature cloud map drawing is performed according to the method of the present invention:

[0082] Step 1: Drilling is performed in the tunnel using a drilling rig. Bever Team Online software is used to analyze and summarize the drilling parameters of the drilling rig. The drilling rig acquires drilling measurement parameters once per drilling depth (d), including: measuring point number, hole data, and drilling measurement data. The drilling measurement data includes the drilling depth (m) corresponding to the measuring point number, average advance speed (m / min), average impact pressure (bar), average advance pressure (bar), average rotational pressure (bar), average water pressure (bar), and average water flow rate (L / min). 1 bar = 100 kPa. See Table 1:

[0083] Table 1. Some drilling parameters collected during the drilling process of the rock drilling rig tunnel.

[0084]

[0085]

[0086] Step 2: Store the borehole data k corresponding to each measuring point number Q and the drilling measurement data Data (depth d, propulsion speed V, impact pressure Pi, propulsion pressure Pt, rotation pressure Pn, water pressure, water flow rate) at the measuring point as a vector A{Q, k, Data} in the database.

[0087] Step 3: Using the spatial coordinate transformation formula, perform coordinate transformation on the three-dimensional spatial coordinates, three-dimensional azimuth angle, and depth of each measuring point according to formulas (1) to (3) to obtain the three-dimensional spatial coordinates (Xi, Yi, Zi) of the i-th drilling depth d at each measuring point. Update the three-dimensional spatial coordinates (X0, Y0, Z0) of vector A to i sets of three-dimensional spatial coordinates as vector A of the i-th drilling depth d at each measuring point. di Saved into the database, as shown in Table 2 below:

[0088] Table 2 shows the parameters after three-dimensional coordinate transformation of the original data.

[0089] Kong Hao Depth m Point X Point Y Point Z 135 0.334 -1.483717245 3.936524721 -0.012716331 135 0.354 -1.464117858 3.938062249 -0.016390573 135 0.374 -1.444518472 3.939599778 -0.020064815 135 0.394 -1.424919085 3.941137306 -0.023739057 135 0.414 -1.405319698 3.942674834 -0.027413299 135 0.434 -1.385720312 3.944212362 -0.031087541 135 0.454 -1.366120925 3.94574989 -0.034761782 135 0.474 -1.346521539 3.947287419 -0.038436024 135 0.494 -1.326922152 3.948824947 -0.042110266 135 0.514 -1.307322766 3.950362475 -0.045784508 135 0.534 -1.287723379 3.951900003 -0.04945875 135 0.554 -1.268123993 3.953437532 -0.053132992

[0090] Step 4: Process the measurement-while-drilling data. Based on the relationship between drilling specific energy and rock strength, calculate the rock compressive strength at the drilling depth d of the first to i drilling operations at each measuring point. As shown in Table 3.

[0091] Impact energy of rock drilling rig:

[0092] Chisel and smash to death:

[0093] Uniaxial compressive strength δ c =0.205α

[0094] Table 3 Uniaxial compressive strength at different depths of measuring point A on the rock drilling rig

[0095]

[0096] Step 5: Based on the ratio of the rock's compressive strength to its tensile strength, calculate the rock tensile strength σ at each measuring point during the first to i drilling depths d, and replace the vector A at the measuring point during the first to i drilling depths d. di The drilling measurement data (Data) is collected and saved together as a vector A. di σ{Q,k,σ} is stored in the database as an index characterizing the explosiveness of rocks;

[0097] Table 4 Tensile strength at different depths of measuring point A on the rock drilling rig

[0098]

[0099]

[0100] Step 6: For the new vector, perform iterative calculations of cross-section acquisition and bisection interpolation to obtain n complete two-dimensional cross-section vector data. (See below) Figure 2 As shown, the binary search method is used to fill in missing values ​​using the data from the previous and next groups.

[0101] For example: If the measured point with index 135 is missing in the obtained two-dimensional vector group Xi = -1.43, then use X from the table above. i-1 =-1.42 and X i+1 The two sets of data with a value of 1.44 are evaluated using the bisection method.

[0102] Step 7: Plot a scatter plot of the two-dimensional vectors of (Max-min) / 2 complete two-dimensional cross-sections to obtain the tunnel blasting shape. See below. Figure 3 As shown.

[0103] Step 8. Fill in the two-dimensional partition map by using the rock tensile strength of each point in the scatter plot using Kriging interpolation. (Kriging interpolation, also known as spatial local interpolation, mainly involves interpolating the values ​​of the space near each point using the rock tensile strength to ensure the integrity of the two-dimensional region values ​​and obtain a complete partition map.)

[0104] Step 9: As Figure 4 As shown: Repeat steps 7-8 to obtain the partition maps of all two-dimensional sections. Based on these partition maps of all two-dimensional sections, a three-dimensional spatial explosiveness characteristic cloud map of the tunnel rock can be drawn.

[0105] On the other hand, the present invention also provides a system for generating a three-dimensional spatial explosiveness feature cloud map of rocks, the system comprising:

[0106] The first processing module is used to perform drilling operations using a rock drilling rig. Every time the drilling depth is d, the drilling parameters are acquired once and stored as a vector in the database. The drilling parameters include: measuring point number, hole data and drilling measurement data. The hole data includes the three-dimensional spatial coordinates of the borehole opening corresponding to the measuring point number and the three-dimensional azimuth angle of the drilling direction. The drilling measurement data includes various measurement parameters corresponding to the measuring point number.

[0107] The second processing module is used to transform the hole data in the database using the spatial coordinate transformation formula to obtain the three-dimensional spatial vector of all boreholes in the working direction perpendicular to the excavation face.

[0108] The third processing module is used to calculate the rock tensile strength of each three-dimensional spatial coordinate equidistantly arranged along the three-dimensional spatial vector of the borehole based on the relationship between drilling specific energy and rock strength, so as to replace the drilling measurement data in the database and form a new vector in the database.

[0109] The fourth processing module is used to divide the new vector obtained into multiple two-dimensional sections at equal intervals along the working direction, and to draw a scatter plot based on the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional sections and the tensile strength value of the rock.

[0110] The fifth processing module is used to fill in the rock tensile strength value of each scatter point in the scatter plot using Kriging interpolation to obtain a two-dimensional partition map;

[0111] The sixth processing module is used to draw a three-dimensional spatial explosiveness characteristic cloud map of the rock using all the two-dimensional partition maps.

[0112] It should be noted that the rock three-dimensional spatial explosiveness feature cloud map generation system provided in this embodiment of the invention can execute the rock three-dimensional spatial explosiveness feature cloud map generation method described in any of the above embodiments during specific operation, which will not be elaborated in this embodiment.

[0113] Figure 5 This is a schematic diagram of the structure of the electronic device provided by the present invention, such as... Figure 5 As shown, the electronic device may include a processor 510, a communication interface 520, a memory 530, and a communication bus 540. The processor 510, communication interface 520, and memory 530 communicate with each other via the communication bus 540. The processor 510 can call logical instructions from the memory 530 to execute a method for generating a three-dimensional explosiveness feature cloud map of the rock.

[0114] Furthermore, the logical instructions in the aforementioned memory 530 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0115] On the other hand, the present invention also provides a computer program product, the computer program product including a computer program stored on a non-transitory computer-readable storage medium, the computer program including program instructions, and when the program instructions are executed by a computer, the computer is able to execute the method for generating a three-dimensional spatial explosiveness feature cloud map of rocks provided in the above embodiments.

[0116] In another aspect, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for generating a three-dimensional spatial explosiveness feature cloud map of rock provided in the above embodiments.

[0117] The device embodiments described above are merely illustrative, and the units described as separate components may or may not be physically separate. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0118] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for generating a three-dimensional spatial explosiveness feature cloud map of rock, characterized in that, include: Drilling operations are carried out using a rock drilling rig. Every time the drilling depth is d, the drilling parameters are acquired once and stored as a vector in the database. The parameters while drilling include: measuring point number, hole data, and measurement data while drilling. The hole data includes the three-dimensional spatial coordinates of the borehole opening corresponding to the measuring point number and the three-dimensional azimuth angle of the drilling direction. The measurement data while drilling includes various measurement parameters corresponding to the measuring point number. Using the spatial coordinate transformation formula, the hole data in the database is transformed to obtain the three-dimensional spatial vector of all boreholes in the working direction perpendicular to the excavation face; Based on the relationship between drilling specific energy and rock strength, the tensile strength of the rock is calculated for each three-dimensional spatial coordinate equidistantly arranged along the three-dimensional spatial vector of the borehole, in order to replace the drilling measurement data in the database and form a new vector in the database. The new vector obtained is divided into multiple two-dimensional sections at equal intervals along the working direction. A scatter plot is drawn based on the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional sections and the tensile strength value of the rock. A two-dimensional partition map is obtained by filling in the rock tensile strength value of each scatter point in the scatter plot using the Kriging interpolation method; Using all the two-dimensional partition maps, draw a three-dimensional spatial cloud map of the explosiveness characteristics of the rock; The step of determining the rock tensile strength at each equidistant three-dimensional spatial coordinate along the three-dimensional spatial vector of the borehole includes: calculating the drilling specific energy and impact energy required for the rock drilling rig to drill a unit volume of rock through rock breaking mechanics analysis using various drill bits; calculating the rock breaking specific energy based on the relevant conversion relationship between impact energy and breaking specific energy; calculating the rock compressive strength at the first to i-th drilling depths d at each measuring point based on the relationship between breaking specific energy and the saturated uniaxial compressive strength of the rock; and calculating the rock tensile strength at the first to i-th drilling depths d at each measuring point based on the ratio of the rock's compressive strength to its tensile strength. In cases where data is missing in a two-dimensional cross-section, the missing data is filled in to obtain multiple complete two-dimensional cross-sections.

2. The method for generating a three-dimensional spatial explosiveness characteristic cloud map of rock according to claim 1, characterized in that, The spatial coordinate transformation formula includes: Xi = X0 + L × cos(α) × cos(β); Yi = Y0 + L × cos(β) × sin(α); Zi = Z0 - L × sin(β); In the formula, X0, Y0, Z0 represent the three-dimensional spatial coordinates of the borehole opening, Xi, Yi, Zi represent the three-dimensional spatial coordinates of the i-th drilling depth d, L represents the cumulative depth of the borehole after the i-th drilling depth d, and α, β and γ represent the Euler angles of the borehole offset in the X, Y and Z directions, respectively.

3. The method for generating a three-dimensional spatial explosiveness characteristic cloud map of rock according to claim 1, characterized in that, Through rock-breaking mechanics analysis of various drill bits, the drilling specific energy and impact energy required for a rock drilling rig to drill a unit volume of rock were calculated, including: ； ； Where e represents the drilling specific energy, Impact energy, Pi represents average impact pressure, Pt represents average propulsion pressure, Pn represents average rotational pressure, and V represents average propulsion speed.

4. The method for generating a three-dimensional spatial explosiveness characteristic cloud map of rock according to claim 1, characterized in that, A scatter plot was drawn based on the three-dimensional spatial coordinates of the intersection points of each borehole at a two-dimensional cross-section and the tensile strength value of the rock, including: The location of the scattered points is determined based on the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional cross-section; The color value of the scatter plot is determined based on the tensile strength value of the rock. The color value is linearly related to the tensile strength of the rock.

5. The method for generating a three-dimensional spatial explosiveness characteristic cloud map of rock according to claim 4, characterized in that, The color value is a grayscale value.

6. A system for generating three-dimensional spatial explosiveness characteristic cloud maps of rocks, characterized in that, include: The first processing module is used to perform drilling operations using a rock drilling rig. It acquires the drilling parameters once every time the drilling depth is d and stores them as a vector in the database. The parameters while drilling include: measuring point number, hole data, and measurement data while drilling. The hole data includes the three-dimensional spatial coordinates of the borehole opening corresponding to the measuring point number and the three-dimensional azimuth angle of the drilling direction. The measurement data while drilling includes various measurement parameters corresponding to the measuring point number. The second processing module is used to transform the hole data in the database using the spatial coordinate transformation formula to obtain the three-dimensional spatial vector of all boreholes in the working direction perpendicular to the excavation face. The third processing module is used to calculate the rock tensile strength of each three-dimensional spatial coordinate equidistantly arranged along the three-dimensional spatial vector of the borehole based on the relationship between drilling specific energy and rock strength, so as to replace the drilling measurement data in the database and form a new vector in the database. The fourth processing module is used to divide the new vector obtained into multiple two-dimensional sections at equal intervals along the working direction, and to draw a scatter plot based on the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional sections and the tensile strength value of the rock. The fifth processing module is used to fill in the rock tensile strength value of each scatter point in the scatter plot using Kriging interpolation to obtain a two-dimensional partition map; The sixth processing module is used to draw a three-dimensional spatial explosiveness characteristic cloud map of the rock using all the two-dimensional partition maps; The step of determining the rock tensile strength at each equidistant three-dimensional spatial coordinate along the three-dimensional spatial vector of the borehole includes: calculating the drilling specific energy and impact energy required for the rock drilling rig to drill a unit volume of rock through rock breaking mechanics analysis using various drill bits; calculating the rock breaking specific energy based on the relevant conversion relationship between impact energy and breaking specific energy; calculating the rock compressive strength at the first to i-th drilling depths d at each measuring point based on the relationship between breaking specific energy and the saturated uniaxial compressive strength of the rock; and calculating the rock tensile strength at the first to i-th drilling depths d at each measuring point based on the ratio of the rock's compressive strength to its tensile strength. In cases where data is missing in a two-dimensional cross-section, the missing data is filled in to obtain multiple complete two-dimensional cross-sections.

7. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method for generating a three-dimensional spatial explosiveness feature cloud map of rock as described in any one of claims 1 to 5.

8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for generating a three-dimensional spatial explosiveness feature cloud map of rock as described in any one of claims 1 to 5.

Images