Method and system for generating rock three-dimensional space explosibility characteristic cloud picture

Abstract

The invention provides a method and system for generating a rock three-dimensional space explosibility characteristic cloud picture, and belongs to the technical field of blasting engineering.The method comprises the steps that a drill jumbo is used for drilling operation, and a depth sensor obtains measurement-while-drilling parameters as vectors to be stored in a database every time d is drilled; original data are calculated, and hole data are converted into three-dimensional space coordinates to be stored in a database; according to the relation between the drilling specific energy and the rock strength, the rock compressive strength of each measuring point is calculated; the corresponding rock tensile strength is further obtained and stored in a database; carrying out cross section acquisition and dichotomy interpolation cycle calculation on the new vector to obtain n complete two-dimensional cross section vector data; drawing a scatter diagram to obtain an excavation shape; and a two-dimensional partition map is obtained through filling by a Kriging interpolation method, and the three-dimensional vector data is cut into n two-dimensional explosive thermodynamic diagrams after n times of circulation. According to the invention, the collected data is automatically processed to obtain the explosibility result, and the engineering efficiency is improved.

Description

Technical Field

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

[0002] Tunnel blasting is a type of blasting engineering technology. It is an efficient and adaptable underground engineering construction method that is widely used in projects that cross complex geological conditions. Its main features include rapid crushing of large volumes of rock, significant improvement in construction efficiency, reduced labor intensity, and good economic efficiency. Modern blasting technology can effectively reduce the impact on the surrounding environment and ensure construction safety through precise design and control.

[0003] With the continuous development of blasting engineering technology, in drilling technology, traditional drilling technology needs to be measured and collected, which not only takes a lot of time, but also has a lot of errors. At present, in drilling technology, sensors and intelligent collection systems have gradually appeared, which can collect and transmit back the parameters of the drilling trolley: such as drilling propulsion speed, drilling propulsion pressure, average rotary pressure, water pressure, water flow, depth and other parameters in real time, and in tunnel drilling, the relative position of the drilled hole can be set, and the progress of the drilled tunnel can be monitored in real time. This technology has greatly improved the efficiency of drilling work and reduced the errors caused by data collection and analysis. However, there are still some problems in data analysis. In each tunnel and even in each project, the amount of data collected and transmitted back by the drilling trolley is huge. In the analysis of data, a lot of time will be wasted, so that the time of the next round of drilling or blasting will be delayed, causing huge time, personnel and even economic losses. Summary of the invention

[0004] The present invention provides a method and system for generating a three-dimensional rock explosiveness characteristic cloud map, which is used to solve the defect in the prior art that it is difficult to efficiently and effectively process the data collected by a rock drilling rig.

[0005] In a first aspect, the present invention provides a method for generating a three-dimensional rock explosiveness characteristic cloud map, comprising: using a rock drilling rig to perform a drilling operation, obtaining a drilling parameter once every time a drilling depth d is reached and storing it in a database as a vector; the drilling parameter comprises: a measuring point number, hole data and drilling measurement data, the hole data comprises a three-dimensional spatial coordinate of a borehole opening corresponding to the measuring point number and a three-dimensional azimuth of a drilling direction, and the drilling measurement data comprises a plurality of measurement parameters corresponding to the measuring point number; using a spatial coordinate conversion formula, converting the hole data in the database to obtain a three-dimensional coordinate of all boreholes in an operating direction perpendicular to an excavation surface. space vector; according to the relationship between drilling specific energy and rock strength, the rock tensile strength of each three-dimensional space coordinate equidistantly arranged on the three-dimensional space vector along the borehole is obtained to replace the measurement while drilling 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 according to the three-dimensional space coordinates of the intersection of each borehole at the two-dimensional section and the rock tensile strength value; the rock tensile strength value of each scattered point in the scatter plot is filled by the Kriging interpolation method to obtain a two-dimensional partition map; all the two-dimensional partition maps are used to draw a three-dimensional space explosive characteristic cloud map of the rock.

[0006] According to the method for generating a three-dimensional rock blastability characteristic cloud map provided by the present invention, the spatial coordinate conversion formula includes:

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

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

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

[0010] Among them, X0, Y0, Z0 represent the three-dimensional spatial coordinates of the borehole mouth, 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 γ represent the Euler angles of the borehole offset in the X direction, Y direction and Z direction respectively.

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

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

[0013]

[0014]

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

[0016] According to the method for generating a three-dimensional rock blastability characteristic cloud map 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 is missing data in the two-dimensional section, the missing data is supplemented to obtain multiple complete two-dimensional sections.

[0017] According to the method for generating a three-dimensional rock explosiveness characteristic cloud map provided by the present invention, a scatter plot is drawn according to the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional section and the rock tensile strength value, including: determining the position of the scatter point according to the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional section; determining the color value of the scatter point according to the rock tensile strength value; wherein the color value is linearly related to the rock tensile strength value.

[0018] According to the method for generating a three-dimensional rock blastability characteristic cloud map provided by the present invention, the color value is a grayscale value.

[0019] In a second aspect, the present invention further provides a system for generating a three-dimensional rock blastability characteristic cloud map, comprising:

[0020] The first processing module is used to perform drilling operations using a rock drilling rig, and obtain a while-drilling parameter once every drilling depth d and store it in a database as a vector; the while-drilling parameter includes: a measuring point number, hole data, and while-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 of the drilling direction, and the while-drilling measurement data includes a variety of measurement parameters corresponding to the measuring point number;

[0021] The second processing module is used to transform the hole data in the database using a spatial coordinate transformation formula to obtain the three-dimensional space vectors of all the holes in the working direction perpendicular to the excavation surface;

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

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

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

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

[0026] In a third aspect, the present invention provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of a method for generating a three-dimensional rock explosiveness characteristic cloud map as described in any one of the above-mentioned methods are implemented.

[0027] In a fourth aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of a method for generating a three-dimensional rock explosiveness characteristic cloud map as described in any one of the above.

[0028] The method and system for generating a three-dimensional rock explosiveness characteristic cloud map provided by the present invention can draw a three-dimensional rock explosiveness characteristic cloud map based on the drilling parameters of the rock drilling rig and the rock tensile strength. During the drilling process, the rock's solidity and strength can be effectively predicted and the rock explosiveness can be visualized. This has a significant practical impact on improving construction efficiency.

[0029] The present invention can obtain the parameters and cloud maps after coordinate conversion in real time. Engineering personnel can monitor the drilling rig, view and analyze the data while working, and can timely discover and solve problems, as well as improve the efficiency of charge calculation in the charging process in the blasting field; and reduce the personnel requirements, time and economic requirements in the drilling and blasting process, thereby improving engineering benefits. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0031] Figure 1 It is a schematic flow chart of a method for generating a three-dimensional rock blastability characteristic cloud map provided by the present invention;

[0032] Figure 2 It is a number missing two-dimensional scatter plot provided by the present invention;

[0033] Figure 3 It is a two-dimensional scatter plot after supplementation provided by the present invention;

[0034] Figure 4 It is the explosive area division diagram provided by the present invention;

[0035] Figure 5 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0036] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0037] It should be noted that, in the description of the embodiments of the present invention, the terms "include", "comprise" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also includes other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0038] The terms "first", "second", etc. in this application are used to distinguish similar objects, and are not used to describe a specific order or sequence. It should be understood that the terms used in this way can be interchangeable under appropriate circumstances, so that the embodiments of the present application can be implemented in an order other than those illustrated or described here, and the objects distinguished by "first", "second", etc. are generally of the same type, and the number of objects is not limited. For example, the first object can be one or more.

[0039] Combine the following Figure 1-Figure 5 The invention describes a method and system for generating a three-dimensional rock explosiveness characteristic cloud map provided by an embodiment of the invention.

[0040] Figure 1 It is a schematic diagram of the process of the method for generating the three-dimensional rock blastability characteristic cloud map provided by the present invention, such as Figure 1 As shown, including but not limited to the following steps:

[0041] S1. Use a drilling rig to perform rock drilling operations, and use the Bever Team Online software to analyze and summarize the drilling parameters of the rock drilling rig. The drilling rig obtains the drilling parameters once through the sensor every time it drills a depth of d. The drilling parameters include: measuring point number, hole data, and drilling measurement data. The drilling measurement data includes the depth of the borehole corresponding to the measuring point number (m), average advancement speed (m / min), average impact pressure (bar), average advancement pressure (bar), average rotation pressure (bar), average water pressure (bar), average water flow (L / min) and other measurement parameters. 1bar = 100kPa. The hole data includes the three-dimensional spatial coordinates (X0, Y0, Z0) of the borehole mouth corresponding to the measuring point number, and the three-dimensional azimuth (α, β, γ) of the drilling direction.

[0042] S2. The measuring point hole data k corresponding to each measuring point serial number Q and the downhole measurement data Data (depth L, thrust speed V, impact pressure Pi, thrust 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 conversion formula, the three-dimensional spatial coordinates, three-dimensional azimuth and depth of each measuring point are converted according to formulas (1) to (3) to obtain the three-dimensional spatial coordinates (Xi, Yi, Zi) of the i-th drilling depth d of the borehole at each measuring point, and the three-dimensional spatial coordinates (X0, Y0, Z0) of the vector A are updated to i sets of three-dimensional spatial coordinates as the vector A of the i-th drilling depth d of the borehole at each measuring point. di Save into the database.

[0044] The hole data directly extracted from the drilling rig, including the depth of the borehole, the three-dimensional spatial coordinates of the borehole, and the three-dimensional azimuth, need to be converted using a special algorithm into the format used in the method of drawing the three-dimensional space explosive characteristic cloud map. 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(β) Formula (3)

[0048] Where X0, Y0, Z0 represent the three-dimensional spatial coordinates of the borehole mouth, 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 (m), α, β and γ represent the Euler angles of the borehole offset in the X direction, Y direction and Z direction respectively. The hole is not offset and rotated in the Z direction when drilling, so the γ Euler angle is not used in this space transformation method.

[0049] α: The angle of rotation around the Z axis, used to represent the change in the X and Y directions of the drilling.

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

[0051] γ: The rotation angle around the X axis.

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

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

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

[0055] S4. Process the measurement while drilling data, and calculate the rock compressive strength at the depth d of the first to i drillings of the borehole at each measuring point according to the relationship between the drilling specific energy and the rock strength.

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

[0057] The formula is as follows:

[0058]

[0059]

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

[0061] The drilling parameters of the drilling rig and the parameters of the drilling rig are collected to calculate the rock breaking energy of each part, and finally the drilling specific work and the percentage of each part of energy can be calculated. The impact of the rock drilling rig can produce stress concentration, and the average proportion of this energy in the rock breaking process is 88%, that is, the impact energy of the drilling specific energy e i Play a leading role.

[0062] Further, according to the correlation between the impact energy ei, which plays a dominant role in the drilling specific energy, and the crushing specific work a; a = 0.86e i , the specific work of rock crushing can be calculated.

[0063] Where: e i —Impact energy (MJ / m 3 ); a—crushing specific work (MJ / m 3 ).

[0064] Collect a large amount of information about rock crushing specific work a and rock saturated uniaxial compressive strength δ c The research data of the rock crushing specific work and the rock saturated uniaxial compressive strength formula are obtained: δ c =0.205α;

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

[0066] According to the relationship between the specific work of crushing and the saturated uniaxial compressive strength of rock, the rock compressive strength δ at the 1st to ith drilling depth d of the borehole at each measuring point is calculated. c .

[0067] S5. According to the ratio of the rock compressive strength to the tensile strength, calculate the rock tensile strength σ at the 1st to ith drilling depth d of each measuring point, and replace the vector A at the 1st to ith drilling depth d of the measuring point di The measurement data while drilling Data in the same file are saved together as vector A. di σ{Q,k,σ}, as an index to characterize the explosiveness of rock, is stored in the database;

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

[0069] Add the serial numbers Q of all measurement points to the array N, and traverse all A according to Q di σ, find the vector A di The maximum value Max and the minimum value min of the X coordinate value of σ are recorded every time the drilling rig drills 2 cm along the straight line vertical to the excavation surface. Therefore, the vector A is converted into 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 the array N j Xi, Yi, Zi, σ i If the hole data and measurement while drilling data of the borehole corresponding to the jth measuring point number are missed during the i-th collection due to communication failure, the two-dimensional vector C i The jth measuring point number q j If missing, the following algorithm is used for interpolation:

[0070] For each two-dimensional vector C i Loop through and check each two-dimensional vector C i All q j Compare with Q. If there is a missing measurement point number, add the missing measurement point number to array k in the loop. i In the table, record the X of the boreholes that do not have missing measuring point numbers. i The value is also equal to the X value of the borehole with the missing point number. i value; according to k i The point number in the three-dimensional vector A di Find k in i The drilling measurement data and spatial 3D coordinates of all the measuring point numbers in the same place are found in the drilling measurement data and spatial 3D coordinates with the recorded X coordinates. i The measurement while drilling data and spatial three-dimensional coordinates whose absolute value of the difference is less than Δd (Δd = 2 cm, the measurement while drilling data includes the advancement speed V, impact pressure Pi, advancement pressure Pt, rotation pressure Pn, water pressure, and water flow of the borehole of the measuring point sequence 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} and {X i , Y i , Z i ,Data i} are the same), and the two sets of data are used to calculate the data corresponding to the missing measurement point number {X i , Y i , Z i ,Data i}, then according to Data i Calculate the rock tensile strength σ according to steps S4 to S5 i After completing the two-dimensional vector C i The data corresponding to each missing measurement point number in {X i , Yi , Z i , σ i}.

[0071] If there is a drill hole with a missing measuring point number, i If the value is greater than the maximum value of the X-coordinate depth of the borehole or less than the minimum value of the X-coordinate depth of the borehole, the measurement while drilling data and the spatial three-dimensional coordinates corresponding to the maximum or minimum value of the X-coordinate depth of the borehole are used as the three-dimensional vector A di The drilling measurement data and spatial three-dimensional coordinates of the borehole with the missing measurement point number are then i Calculate the rock tensile strength σ according to steps S4 to S5 i , after completing the two-dimensional vector C i The data corresponding to each missing measurement point number in {X i , Y i , Z i , σ i}, because the drill did not pass through this point during drilling, but in order to ensure the integrity of the two-dimensional map, 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] The partial data {X i ,{Yi,Zi},σ i}Create arrays respectively, and use the scatter method in the matplotlib library in Python (the scatter method is a method for drawing scatter plots in the matplotlib library, which can draw points according to the x, y, and z coordinates, and determine the color of the points according to the value of the rock tensile strength) to draw the scatter plot (the horizontal coordinate of the scatter plot is the y value of the two-dimensional vector, and the vertical coordinate is the z value).

[0075] S8. Use the Kriging interpolation method to fill in the rock tensile strength of each point in the scatter plot to obtain a two-dimensional partition map. (The Kriging interpolation method is also called the spatial local interpolation method. In this step, the rock tensile strength of each point is used to interpolate the value of the space near the point to ensure the integrity of the two-dimensional area value and obtain a complete partition map).

[0076] S9. Repeat steps S7-S8 to obtain the partition diagrams of all two-dimensional cross sections. A three-dimensional space blastability characteristic cloud diagram of the excavation area can be drawn based on the partition diagrams of all two-dimensional cross sections.

[0077] By using the method for drawing a three-dimensional rock explosiveness characteristic cloud map provided by the present invention, a three-dimensional rock explosiveness characteristic cloud map can be drawn according to the drilling parameters of the rock drilling rig and the rock tensile strength. During the drilling process, the rock solidity and strength can be effectively predicted and the rock explosiveness can be visualized. This has a significant practical impact on improving construction efficiency.

[0078] In order to more clearly illustrate the technical solution of the present invention, it is described below in conjunction with a specific engineering scenario.

[0079] Overview of a certain project: the top of the straight wall of the tunnel is a semicircular arch, and the tunnel strata are mainly granite. The construction method of full-section one-time excavation is adopted. The tunnel surrounding rock is good. According to the design requirements of the railway department, excavation construction is carried out in the working direction of the tunnel face. The working direction is perpendicular to the tunnel face and points to the rock.

[0080] The surrounding rock conditions of the tunnel are relatively good. Drilling rigs are used to enter the tunnel from the entrance and exit, and a full-section one-time drilling and blasting plan is adopted, that is, a middle slot hole, surrounding auxiliary holes, and peripheral holes are excavated according to the designed contour and smooth blasting is used.

[0081] The characteristic cloud diagram is drawn according to the method of the present invention:

[0082] Step 1: Use a drilling rig to drill rocks in the tunnel, and use the Bever Team Online software to analyze and summarize the drilling parameters of the rock drilling rig. The drilling rig obtains the drilling measurement parameters once through the sensor every time it drills a depth of d. The drilling measurement parameters include: 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 advancement speed (m / min), average impact pressure (bar), average advancement pressure (bar), average rotation pressure (bar), average water pressure (bar), average water flow (L / min) and other drilling parameters. 1bar = 100kPa. As shown in Table 1:

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

[0084]

[0085]

[0086] Step 2: The measuring point hole data k corresponding to each measuring point serial number Q and the downhole measurement data Data (depth d, thrust speed V, impact pressure Pi, thrust 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}.

[0087] Step 3: Use the spatial coordinate conversion formula to convert the three-dimensional spatial coordinates, three-dimensional azimuth 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 of the borehole 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 of the borehole at each measuring point. di Save it into the database, as shown in Table 2 below:

[0088] Table 2 Parameters after three-dimensional coordinate transformation of original data

[0089] Hole No. Depth m Click 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 and calculate the rock compressive strength at the depth d of the first to i drillings of each measuring point according to the relationship between drilling specific energy and rock strength. See Table 3.

[0091] Rock drilling rig impact energy:

[0092] Chiseling specific work:

[0093] Uniaxial compressive strength δ c =0.205α

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

[0095]

[0096] Step 5: According to the ratio of the rock compressive strength to the tensile strength, calculate the rock tensile strength σ at the 1st to ith drilling depth d of each measuring point, and replace the vector A of the 1st to ith drilling depth d of the measuring point di The measurement data while drilling Data in the same file are saved together as vector A. di σ{Q,k,σ}, as an index to characterize the explosiveness of rock, is stored in the database;

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

[0098]

[0099]

[0100] Step 6: For the new vector, perform cross-section acquisition and bisection interpolation loop calculation to obtain n complete two-dimensional cross-section vector data. Figure 2 As shown, the missing values ​​were filled using the before and after group data using the dichotomy method.

[0101] For example: If the point numbered 135 is missing in the obtained two-dimensional vector group Xi = -1.43, then use the X in the above table. i-1 = -1.42 and X i+1 =1.44 are evaluated using the binary method.

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

[0103] Step 8. Use the rock tensile strength of each point in the scatter plot to fill in a two-dimensional partition map using the Kriging interpolation method. (The Kriging interpolation method is also called the spatial local interpolation method. In this step, the rock tensile strength of each point is used to interpolate the value of the space near the point to ensure the integrity of the two-dimensional area value and obtain a complete partition map.)

[0104] Step 9: Figure 4 As shown: Repeat steps 7 and 8 to obtain the partition diagrams of all two-dimensional sections. According to the partition diagrams of all two-dimensional sections, a cloud diagram of the three-dimensional blastability characteristics of the tunnel rock can be drawn.

[0105] On the other hand, the present invention also provides a system for generating a three-dimensional rock blastability characteristic cloud map, the system comprising:

[0106] The first processing module is used to perform drilling operations using a rock drilling rig, and obtain a while-drilling parameter once every drilling depth d and store it in a database as a vector; the while-drilling parameter includes: a measuring point number, hole data, and while-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 of the drilling direction, and the while-drilling measurement data includes a variety of measurement parameters corresponding to the measuring point number;

[0107] The second processing module is used to transform the hole data in the database using a spatial coordinate transformation formula to obtain the three-dimensional space vectors of all the holes in the working direction perpendicular to the excavation surface;

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

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

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

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

[0112] It should be noted that the system for generating a three-dimensional rock explosiveness characteristic cloud map provided in an embodiment of the present invention can, when in operation, execute the method for generating a three-dimensional rock explosiveness characteristic cloud map described in any of the above embodiments, which will not be elaborated in this embodiment.

[0113] Figure 5 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, wherein the processor 510, the communication interface 520 and the memory 530 communicate with each other through the communication bus 540. The processor 510 may call the logic instructions in the memory 530 to execute the method for generating the three-dimensional rock blastability characteristic cloud map.

[0114] In addition, the logic instructions in the above-mentioned memory 530 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0115] On the other hand, the present invention also provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the method for generating a three-dimensional rock explosiveness characteristic cloud map provided in the above-mentioned embodiments.

[0116] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to execute the method for generating a three-dimensional rock explosiveness characteristic cloud map provided in the above-mentioned embodiments.

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

[0118] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions 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 rock blastability characteristic cloud map, characterized in that: include: The drilling operation is performed using a drilling rig, and the drilling parameters are obtained once every time the drilling depth reaches d, and stored in the database as a vector; 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 and the three-dimensional azimuth of the drilling direction corresponding to the measuring point number. The drilling measurement data includes multiple measurement parameters corresponding to the measuring point number. Using the spatial coordinate conversion formula, the hole data in the database is converted to obtain the three-dimensional space vectors of all the holes in the working direction perpendicular to the excavation surface; According to the relationship between drilling specific energy and rock strength, the rock tensile strength of each three-dimensional space coordinate arranged equidistantly on the three-dimensional space vector along the borehole is obtained to replace the while-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 at equal intervals along the working direction, and a scatter plot is drawn according to 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 scattered point in the scatter plot is filled by Kriging interpolation method to obtain a two-dimensional partition map; All two-dimensional zoning maps are used to draw a three-dimensional rock blastability characteristic cloud map.

2. The method for generating a three-dimensional rock explosiveness characteristic cloud map according to claim 1, characterized in that: The spatial coordinate conversion formula includes: Xi=X0+L×cos(α)×cos(β); Yi=Y0+L×cos(β)×sin(α); Zi = Z0 - L × sin (β); Among them, X0, Y0, Z0 represent the three-dimensional spatial coordinates of the borehole mouth, 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 γ represent the Euler angles of the borehole offset in the X direction, Y direction and Z direction respectively.

3. The method for generating a three-dimensional rock blastability characteristic cloud map according to claim 1, characterized in that: According to the relationship between drilling specific energy and rock strength, the rock tensile strength of each three-dimensional space coordinate equidistantly arranged on the three-dimensional space vector along the borehole is obtained, including: Through the analysis of rock breaking mechanics of various drill bits, the drilling specific energy and impact energy required by the drilling rig to drill a unit volume of rock are calculated; According to the conversion relationship between impact energy and specific work of chipping, the specific work of chipping of rock is calculated; According to the relationship between the specific work of crushing and the saturated uniaxial compressive strength of rock, the compressive strength of rock at the depth d of the first to i drilling of the borehole at each measuring point is calculated.

4. The method for generating a three-dimensional rock blastability characteristic cloud map according to claim 1, characterized in that: Through the analysis of various drill bit rock breaking mechanics, the drilling specific energy and impact energy required for the drilling rig to drill a unit volume of rock are calculated, including: Where, e represents the drilling specific energy, e i impact energy, Pi represents the average impact pressure, Pt represents the average propulsion pressure, Pn represents the average rotation pressure, and V represents the average propulsion speed.

5. The method for generating a three-dimensional rock blastability characteristic cloud map according to claim 1, characterized in that: The obtained new vector is divided into multiple two-dimensional sections at equal intervals along the working direction, including: When there is missing data in the two-dimensional section, the missing data is supplemented to obtain multiple complete two-dimensional sections.

6. The method for generating a three-dimensional rock blastability characteristic cloud map according to claim 1, characterized in that: Draw a scatter plot based on the three-dimensional spatial coordinates of the intersection of each borehole at the two-dimensional section and the rock tensile strength value, including: Determine the position of the scattered points according to the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional cross section; According to the rock tensile strength value, the color value of the scattered points is determined; Therein, the color value is in a linear relationship with the rock tensile strength value.

7. The method for generating a three-dimensional rock blastability characteristic cloud map according to claim 6, characterized in that: The color value is a grayscale value.

8. A system for generating a three-dimensional rock blastability characteristic cloud map, characterized in that: include: The first processing module is used to perform drilling operations using a rock drilling rig, and obtain a while-drilling parameter every time a drilling depth d is reached and store it in a database as a vector; 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 and the three-dimensional azimuth of the drilling direction corresponding to the measuring point number. The drilling measurement data includes multiple measurement parameters corresponding to the measuring point number. The second processing module is used to transform the hole data in the database using a spatial coordinate transformation formula to obtain the three-dimensional space vectors of all the drill holes in the working direction perpendicular to the excavation surface; The third processing module is used to obtain the rock tensile strength of each three-dimensional space coordinate arranged equidistantly on the three-dimensional space vector along the borehole according to the relationship between the drilling specific energy and the rock strength, so as to replace the while drilling measurement data in the database and form a new vector in the database; The fourth processing module is used to divide the obtained new vector into multiple two-dimensional sections at equal intervals along the working direction, and draw a scatter plot according to the three-dimensional spatial coordinates of the intersection points of each borehole at the two-dimensional section and the rock tensile strength value; The fifth processing module is used to fill the rock tensile strength value of each scattered point in the scatter plot by Kriging interpolation method to obtain a two-dimensional partition map; The sixth processing module is used to draw a three-dimensional rock blastability characteristic cloud map using all the two-dimensional partition maps.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the steps of the method for generating a three-dimensional rock explosiveness characteristic cloud map as described in any one of claims 1 to 7 are implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for generating a three-dimensional rock explosiveness characteristic cloud map as described in any one of claims 1 to 7 are implemented.

Images