AU2019439997A1 - Method and device for determining rock stratum parameters - Google Patents
Method and device for determining rock stratum parameters Download PDFInfo
- Publication number
- AU2019439997A1 AU2019439997A1 AU2019439997A AU2019439997A AU2019439997A1 AU 2019439997 A1 AU2019439997 A1 AU 2019439997A1 AU 2019439997 A AU2019439997 A AU 2019439997A AU 2019439997 A AU2019439997 A AU 2019439997A AU 2019439997 A1 AU2019439997 A1 AU 2019439997A1
- Authority
- AU
- Australia
- Prior art keywords
- rock
- drill rod
- drilling
- compressive strength
- unit volume
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 239000011435 rock Substances 0.000 title claims abstract description 214
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000005553 drilling Methods 0.000 claims description 142
- 238000012360 testing method Methods 0.000 claims description 57
- 238000012417 linear regression Methods 0.000 claims description 10
- 238000003064 k means clustering Methods 0.000 claims description 7
- 230000001419 dependent effect Effects 0.000 claims description 6
- 230000000149 penetrating effect Effects 0.000 claims description 6
- 230000035515 penetration Effects 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005520 cutting process Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000003245 coal Substances 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000004570 mortar (masonry) Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 208000008918 voyeurism Diseases 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0218—Calculations based on experimental data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/025—Geometry of the test
- G01N2203/0252—Monoaxial, i.e. the forces being applied along a single axis of the specimen
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Earth Drilling (AREA)
Abstract
Disclosed are a method and device for determining rock stratum parameters. The method comprises: substituting into an estimation formula of rock uniaxial compressive strength values R
Description
Method and Device for Determining Rock Stratum Parameters
Technical Field
[0001] The present invention relates to a method and a device for determining rock stratum parameters, and belongs to the technical field of geological exploration engineering.
Background Art
[0002] In underground projects such as coal mines and tunnels, etc., the condition of the rock strata has a great influence on the safety of the project. At present, the information on the rock strata is mainly learned through rock coring, drilling and peeping, etc., but those methods affect the progress of the project and are time and labor consuming, and can't be used to learn about the rock stratum environment of the project in real time.
[0003] Presently, the logging while drilling technology is relatively matured in the oil industry, but hasn't been widely applied in underground projects such as coal mines, owing to its severe limitations. Although some researchers have put forward a formation detection method based on drilling cuttings, but it is difficult to collect rock cuttings; especially for the rock strata with water-dissolvable lithology, the rock cuttings thereof can't be collected. Consequently, continuous judgment can't be made for the formation structure, and the method can't be used to predict the uniaxial compressive strength of the rock.
[0004] Therefore, many researchers have put forward a method for real-time identification of rock strata based on the parameters of a jumbolter while drilling, but there is no effective way to apply the method widely in China at present.
Contents of the Invention
[0005] To solve the above technical problems in the prior art, the present invention provides a method and a device for determining rock stratum parameters.
[0006] The method for determining rock stratum parameters provided by the present invention includes: obtaining g test data groups acquired at different data acquisition times t during rock stratum drilling, wherein the test data in the test data groups includes work values WF of the axial force of the drill rod when drilling through unit volume of rock, work values WM of the torque of the drill rod when drilling through unit volume of rock, and work values Wf of the friction between the drill bit and borehole bottom when drilling through unit volume of rock; substituting the test data in the g test data groups into an estimation formula of uniaxial compressive strength R, of the rock to determine g uniaxial compressive strength values Rei of the rock that respectively correspond to the g test data groups acquired at different data acquisition times t, wherein Rei is the uniaxial compressive strength value of the rock that corresponds to the ithtest data group acquired in the sequence of acquisition times; forming a data sequence A= (R 1 , R,2 ,
R 3 ,..., Rei,..., R,) from the determined g uniaxial compressive strength values Rei of the rock in accordance with the sequence of data acquisition times, and analyzing the data sequence A with K-means clustering analysis software to determine a final categorization result, wherein the final categorization result includes a best number of categories d, and subsequences in each category; wherein the best number of categorizes d, is the total number of rock strata divided according to the difference in the uniaxial compressive strength of the rock; if the eth subsequence is Ae= (R,, Rca ,..., Re), the average value R of the uniaxial compressive strength of the eth layer of rock strata is:
1 R= 1 R (1) b-a+1 i,
[0007] Optionally, the method for determining rock stratum parameters provided by the present invention further includes: obtaining an expression V(t) of the drilling rate V of the drill rod that changes with the data acquisition time t, wherein the drilling rate of the drill rod is the depth of the drill rod penetrating into the rock stratum per second, in unit of m/s; the method for obtaining the expression V(t) of the drilling rate V of the drill rod includes: fitting the drilling rates V of different drill rods and the acquisition times t of the drilling rates V to obtain the expression V(t); determining an interval T= [t,, tb ] of the data acquisition time t corresponding to the uniaxial compressive strength value of the rock in the eth subsequence A,= (Ra, Rca,..., R); then the thickness of the eth layer of rock strata is:
L= fhV(t) dt (2)
[0008] Optionally, the method for determining rock stratum parameters provided by the present invention further includes: calculating the work value WF of the axial force of the drill rod when drilling through unit volume of rock according to a formula (3):
F WF = A;Tr2 (3)
calculating the work value WM of the torque of the drill rod when drilling through unit volume of rock according to a formula (4):
Mn WM M (4) 3u 30r2V
calculating the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock according to a formula (5):
p7,Fn W = O r (5) 60AYrV
where, F is the thrust force of the drill rod, in unit of N; A is a borehole enlargement coefficient, i.e., a ratio of the cross-sectional area of the borehole to the cross-sectional area of the drill rod; r is the bit diameter of the drill rod, in unit of m; M is the torque of the drill rod, in unit of N-m; n is the rotation speed of the drill rod, in unit of r/min.; V is the drilling rate of the drill rod, i.e., the penetration depth of the drill rod into the rock stratum per second, in unit of m/s; p is a coefficient of friction between the drill bit and the borehole bottom, and is determined to be 0.21.
[0009] Optionally, the method for determining the estimation formula of uniaxial compressive strength R, of the rock in the method for determining rock stratum parameters provided by the present invention includes: taking the work value WF of the axial force of the drill rod when drilling through unit volume of rock, the work value WM of the torque of the drill rod when drilling through unit volume of rock and the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock as independent variables, taking the uniaxial compressive strength R, of the rock as a dependent variable, and obtaining the estimation formula with a linear regression method.
[0010] In addition, the present invention further provides a device for determining rock stratum parameters, which includes: an acquisition module configured to obtain g test data groups acquired at different data acquisition times t during rock drilling, wherein the test data in the test data groups includes work values WF of the axial force of the drill rod when drilling through unit volume of rock, work values WM of the torque of the drill rod when drilling through unit volume of rock, and work values Wf of the friction between the drill bit and borehole bottom when drilling through unit volume of rock; the device further includes a determination module configured to substitute the test data in the g test data groups into an estimation formula of uniaxial compressive strength R, of the rock to determine g uniaxial compressive strength values Rej of the rock that respectively correspond to the g test data groups acquired at different data acquisition times t, wherein Rej is the uniaxial compressive strength value of the rock that corresponds to the ith test data group acquired in accordance with the sequence of acquisition times; the determination module is further configured to form a data sequence A= (R 1 , R, 2 , R,..., Rei,..., R,) from the determined g uniaxial compressive strength values Rej of the rock in the sequence of data acquisition times, and analyze the data sequence A with K-means clustering analysis software to determine a final categorization result, wherein the final categorization result includes a best number of categories ds and subsequences in each category; wherein the best number of categorizes d, is the total number of rock strata divided according to the difference in the uniaxial compressive strength values of the rock; if the eth subsequence is A= (R, R,,,..., Re), the average value R of the uniaxial compressive strength of the eth layer of rock strata is:
1 b RE=Ra (1) b- a+1 =
[0011] Optionally, the acquisition module in the device for determining rock stratum parameters provided by the present invention is further configured to obtain an expression V(t) of the drilling rate V of the drill rod that changes with the data acquisition time t, wherein the drilling rate of the drill rod is the depth of the drill rod penetrating into the rock stratum per second, in unit of m/s; the method for obtaining the expression V(t) of the drilling rate V of the drill rod includes: fitting the drilling rates V of different drill rods and the acquisition times t of the drilling rates V to obtain the expression V(t); determining an interval T= [ t, , tb ] of the data acquisition time t corresponding to the uniaxial compressive strength value of the rock in the eth subsequence A,= (Rea, Reg,. Re, ); then the thickness of the eth layer of rock strata is:
L =fthV(t) ta dt (2)
[0012] Optionally, in the device for determining rock stratum parameters provided by the present invention, the acquisition module is further configured to obtain the work value WF of the axial force of the drill rod when drilling through unit volume of rock according to a formula (3):
F WF= 2 (3)
the acquisition module is further configured to obtain the work value Wm of the torque of the drill rod when drilling through unit volume of rock according to a formula (4):
WM = M (4) 330Ar2V 'u
the acquisition module is further configured to obtain the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock according to a formula (5):
p Fn W I ,f= V (5) 60AX rV
where, F is the thrust force of the drill rod, in unit of N; A is a borehole enlargement coefficient, i.e., a ratio of the cross-sectional area of the borehole to the cross-sectional area of the drill rod; r is the bit diameter of the drill rod, in unit of m; M is the torque of the drill rod, in unit of N-m; n is the rotation speed of the drill rod, in unit of r/min.; V is the drilling rate of the drill rod, i.e., the penetration depth of the drill rod into the rock stratum per second, in unit of m/s; U is a coefficient of friction between the drill bit and the borehole bottom, and is determined to be 0.21.
[0013] Optionally, in the device for determining rock stratum parameters provided by the present invention, the determination module is further configured to determine the estimation formula of uniaxial compressive strength R, of the rock; the method for determining the estimation formula of the uniaxial compressive strength R, of the rock includes: taking the work value WF of the axial force of the drill rod when drilling through unit volume of rock, the work value W of the torque of the drill rod when drilling through unit volume of rock and the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock as independent variables, taking the uniaxial compressive strength R, of the rock as a dependent variable, and obtaining the estimation formula with a linear regression method.
[0014] Optionally, the device for determining the rock stratum parameters provided by the present invention further includes a data acquisitor configured to acquire test data when the drill rod drills through unit volume of rock; the data acquisitor includes a drilling rate sensor, a rotation speed sensor, a pressure sensor and a torque sensor; the drilling rate sensor is configured to acquire the drilling rate V; the rotation speed sensor is configured to acquire the rotation speed n of the drill rod; the pressure sensor is configured to acquire the thrust force F of the drill rod; the torque sensor is configured to acquire the torque M of the drill rod; the device further includes a display unit configured to display the test data acquired by the data acquisitor when the drill rod drills through unit volume of rock, the data acquired by the acquisition module, and the data determined by the determination module; the acquisition module is configured to receive the test data acquired by the data acquisitor when the drill rod drills through unit volume of rock and the corresponding uniaxial compressive strength value R,; the determination module processes the data information acquired by the acquisition module.
[0015] With the method and device for determining rock stratum parameters provided by the present invention, the work values of the axial force of the drill rod when drilling through unit volume of rock, the work values of the torque of the drill rod when drilling through unit volume of rock, and the work values of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock in a plurality of test data groups acquired at different data acquisition times in a rock drilling process are substituted into an estimation formula of uniaxial compressive strength value Rc of the rock to obtain a plurality of uniaxial compressive strength values of the rock, the plurality of uniaxial compressive strength values of the rock are analyzed with K-means clustering analysis software to determine subsequences composed of uniaxial compressive strength values of rocks in different rock stratum categories, the average value of uniaxial compressive strength value of a rock stratum is obtained by calculating the average value of uniaxial compressive strength value of the rock in a subsequence. Thus, the geological exploration process is greatly simplified, continuous lithology detection is realized for mine engineering, and manpower and time costs are saved.
Description of Drawings
[0016] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments according to the present invention, and are intended to be used together with the specification to explain the principle of the present invention.
[0017] FIG. 1 is a flow chart of the method for determining rock stratum parameters in an embodiment;
[0018] FIG. 2 is a schematic diagram of the connections of the device for determining rock stratum parameters in an embodiment.
[0019] Specific embodiments of the present invention are shown in the above-mentioned accompanying drawings, and will be described in more detail below. Those drawings and description are not intended to limit the scope of the inventive concept of the present invention in any way, but are used to explain the inventive concept of the present invention to those skilled in the art with reference to specific embodiments. Embodiments
[0020] To make the object, technical solution, and advantages of the embodiments of the present invention understood more clearly, hereunder the technical solution in the embodiments of the present invention will be detailed clearly and completely, with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the embodiments described here are only some embodiments of the present invention rather than all possible embodiments of the present invention.
[0021] On the basis of the embodiments provided in the present invention, all other embodiments obtained by those having ordinary skills in the art without expending any creative labor shall be deemed as falling in the protection scope of the present invention. The following embodiments and the features in the embodiments can be combined with each other, provided there is no conflict.
[0022] As shown in FIG. 1, the method for determining rock stratum parameters provided in this embodiment includes:
[0023] SI: obtaining g test data groups acquired at different data acquisition times t during rock drilling, wherein the test data in the test data groups includes work values WF of the axial force of the drill rod when drilling through unit volume of rock, work values WM of the torque of the drill rod when drilling through unit volume of rock, and work values Wf of the friction between the drill bit and borehole bottom when drilling through unit volume of rock;
exemplarily, in the process of drilling the rock formation with a drilling rig, test data, such as the work value WF of the axial force of the drill rod when drilling through unit volume of rock, the work value WM of the torque of the drill rod when drilling through unit volume of rock, and the work value Wf of the friction between the drill bit and borehole bottom when drilling through unit volume of rock, can be obtained directly or indirectly with the experimental device; when the experimental device acquires such data, it can log the data acquisition time at the same time, and the data acquisition time may be synchronous. Wherein, the "g test data groups" represents multiple test data groups rather than specifically define the number of the data groups. It is only an expression of the number of the acquired test data groups for the convenience of description later.
[0024] S2: substituting the test data in the g test data groups into an estimation formula of uniaxial compressive strength value R, of the rock to determine g uniaxial compressive strength values Ra of the rock that respectively correspond to the g test data groups acquired at different data acquisition times t, wherein Ri is the uniaxial compressive strength value of the rock that corresponds to the ith test data group acquired in accordance with the sequence of acquisition times;
exemplarily, the estimation formula of uniaxial compressive strength value R, of the rock is determined in advance, and the estimation formula of uniaxial compressive strength value R, of the rock may be obtained by directly using a reference formula determined in other literatures, or may be obtained through linear regression of other test data.
[0025] S3: forming a data sequence A= (R,, R,2 , R,..., Ri, ... , R,) from the determined g uniaxial compressive strength values Rei of the rock in accordance with the sequence of data acquisition times, and analyzing the data sequence A with K-means clustering analysis software to determine a final categorization result, wherein the final categorization result includes a best number of categories d, and subsequences in each category; wherein the best number of categorizes ds is the total number of rock strata divided according to the difference in the uniaxial compressive strength value of the rock; if the eth subsequence is Ae= (R, R R,,), the average value R of the uniaxial compressive strength of the eth layer of rock strata is:
1 b RY= Re (1) b- a+1 =
exemplarily, in the drilling process, the lithology of the rock strata changes as the drilling depth is increased, and the g uniaxial compressive strength values Rei of the rock obtained in the process have their own time sequence. In order to obtain continuous lithology of the rock strata, a time sequence constraint is also applied when the data is analyzed. A subsequence is a segment of the data sequence A, and there is no overlap among multiple subsequences. In the final categorization result, the data sequence A is divided into d, subsequences, and the continuity of the data may be exactly the same as that of the data sequence A after multiple subsequences are arranged in a chronological order.
The K-means algorithm is a hard clustering algorithm and is a representative of typical prototype-based objective function clustering methods. It takes a distance from the data point to the prototype as an objective function for optimization, and obtains adjustment rules for iterative operations by using a method of finding an extremum of the function. The algorithm uses a Sum of Squared Errors (SSE) criterion function as the clustering criterion function. R software or Matlab is commonly used software with a K-means algorithm function. The specific process of clustering analysis is completed by means of existing software. In this embodiment, the process and principle of clustering analysis are not described in detail; instead, only the raw data to be inputted and the data types to be obtained are explained exemplarily.
Optionally, on the basis of the above embodiment, the method for determining rock stratum parameters provided in this embodiment further includes the following steps:
[0026] S4: obtaining an expression V(t) of the drilling rate V of the drill rod that changes with the data acquisition time t, wherein the drilling rate of the drill rod is the depth of the drill rod penetrating into the rock stratum per second, in unit of m/s; the method for obtaining the expression V(t) of the drilling rate V of the drill rod includes: fitting the drilling rates V of different drill rods and the acquisition times t of the drilling rates V to obtain the expression V(t); determining an interval T = [t.,, t6 ] of the data acquisition time t corresponding to the uniaxial compressive strength value of the rock in the eth subsequence Ae= (Rca, Rcai,..., R,,); then the thickness of the eth layer of rock strata is: ta L =ftbV(t) dt (2) exemplarily, affected by the uniaxial compressive strength of the rock, the drilling rate of the drill rod in the rock strata may change. The drilling rate sensor may record the drilling rate of the drill rod periodically, e.g., once several milliseconds. According to a plurality of acquired drilling rate values corresponding to different times in a certain time period in the experiment, an expression V(t) of the drilling rate V of the drill rod changing with data acquisition time t can be obtained through nonlinear fitting.
[0027] Alternatively, the average drilling rate may be taken as the drilling rate V of the drill rod, on the premise of obtaining an interval T=[t,, tb ] of data acquisition time t, in conjunction with
the total drilling depth in that time period.
[0028] Optionally, on the basis of the above embodiment, the step Sl may further include: calculating the work value WF of the axial force of the drill rod when drilling through unit volume of rock according to a formula (3):
F WF=pAjrr2 (3)
calculating the work value WM of the torque of the drill rod when drilling through unit volume of rock according to a formula (4):
Mn WM M (4) 'u30Ar2V
calculating the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock according to a formula (5):
p7,Fn Wf = P (5) 60AYLrV
where, F is the thrust force of the drill rod, in unit of N; A is a borehole enlargement coefficient, i.e., a ratio of the cross-sectional area of the borehole to the cross-sectional area of the drill rod; r is the bit diameter of the drill rod, in unit of m; M is the torque of the drill rod, in unit of N-m; n is the rotation speed of the drill rod, in unit of r/min.; V is the drilling rate of the drill rod, i.e., the penetration depth of the drill rod into the rock stratum per second, in unit of m/s; P is a coefficient of friction between the drill bit and the borehole bottom, and is determined to be 0.21.
[0029] Exemplarily, some testing devices can directly obtain the work value WF of the axial force of the drill rod when drilling through unit volume of rock, the work value W of the torque of the drill rod when drilling through unit volume of rock, and the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock; in practice, those values may also be obtained by indirect measurement.
[0030] Optionally, on the basis of the above embodiment, the step S2 may further include: the method for determining the estimation formula of uniaxial compressive strength value R, of the rock in the method for determining rock stratum parameters provided in this embodiment includes: taking the work value WF of the axial force of the drill rod when drilling through unit volume of rock, the work value WM of the torque of the drill rod when drilling through unit volume of rock and the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock as independent variables, taking the uniaxial compressive strength value R, of the rock as a dependent variable, and obtaining the estimation formula with a linear regression method.
[0031] Exemplarily, a linear regression method is a data processing method commonly used in scientific researches, and is a method of working out an empirical formula by obtaining reliable experimental data in advance and then using the empirical formula to estimate the result. For multi-element linear regression, Statistical Product and Service Solutions (SPSS), Microsoft Excel, or other software capable of data processing may be used.
[0032] In addition, as shown in FIG. 2, this embodiment further provides a device for determining rock stratum parameters, which includes: an acquisition module configured to obtain g test data groups acquired at different data acquisition times t during rock drilling, wherein the test data in the test data groups includes work values WF of the axial force of the drill rod when drilling through unit volume of rock, work values WM of the torque of the drill rod when drilling through unit volume of rock, and work values Wf of the friction between the drill bit and borehole bottom when drilling through unit volume of rock; the device further includes a determination module configured to substitute the test data in the g test data groups into an estimation formula of uniaxial compressive strength value R, of the rock to determine g uniaxial compressive strength values Rei of the rock that respectively correspond to the g test data groups acquired at different data acquisition times t, wherein Rej is the uniaxial compressive strength value of the rock that corresponds to the ith test data group acquired in accordance with the sequence of acquisition times; the determination module is further configured to form a data sequence A= (R, R, 2 , R,3 ,..., Re1 ,
. . , R,,) from the determined g uniaxial compressive strength values Rej of the rock in accordance with the sequence of data acquisition times, and analyze the data sequence A with K-means clustering analysis software to determine a final categorization result, wherein the final categorization result includes a best number of categories ds and subsequences in each category; wherein the best number of categorizes ds is the total number of rock strata divided according to the difference in the uniaxial compressive strength value of the rock; if the eth subsequence is A,= (
Rca, R ,cai..., Rc), the average value R of the uniaxial compressive strength of the eth layer of rock strata is:
R= IR, (1) b-a+1 i,
[0033] Optionally, the acquisition module in the device for determining rock stratum parameters provided in this embodiment is further configured to obtain an expression V(t) of the drilling rate V of the drill rod that changes with the data acquisition time t, wherein the drilling rate of the drill rod is the depth of the drill rod penetrating into the rock stratum per second, in unit of m/s; the method for obtaining the expression V(t) of the drilling rate V of the drill rod includes: fitting the drilling rates V of different drill rods and the acquisition times t of the drilling rates V to obtain the expression V(t); determining an interval T,= [ta, tb ] of the data acquisition time t corresponding to the uniaxial compressive strength value of the rock in the eth subsequence A.= (R,, Ra..I. Re, ); then the thickness of the eth layer of rock strata is:
L= V(t) dt (2)
[0034] Optionally, in the device for determining rock stratum parameters provided in this embodiment, the acquisition module is further configured to obtain the work value WF of the axial force of the drill rod when drilling through unit volume of rock according to a formula (3):
F WF=pAjrr2 (3)
the acquisition module is further configured to obtain the work value Wm of the torque of the drill rod when drilling through unit volume of rock according to a formula (4):
Mn WM M (4) "'=30Ar 2V
the acquisition module is further configured to obtain the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock according to a formula (5):
p7 uFn W = O r (5)
where, F is the thrust force of the drill rod, in unit of N; A is a borehole enlargement coefficient, i.e., a ratio of the cross-sectional area of the borehole to the cross-sectional area of the drill rod; r is the bit diameter of the drill rod, in unit of m; M is the torque of the drill rod, in unit of N-m; n is the rotation speed of the drill rod, in unit of r/min.; V is the drilling rate of the drill rod, i.e., the penetration depth of the drill rod into the rock stratum per second, in unit of m/s; p is a coefficient of friction between the drill bit and the borehole bottom, and is determined to be 0.21.
[0035] Exemplarily, the thrust force F of the drill rod may be obtained by subtracting the gravity of the drilling rig from the pressure of the drilling rig contacting with the ground; the borehole enlargement coefficient X is the ratio of the cross-sectional area of the borehole to the cross-sectional area of the drill rod after the borehole is drilled; the torque M of the drill rod may be measured with a torque sensor connected to the drilling rig; the rotation speed n of the drill rod may be measured with a rotation speed sensor connected to the drilling rig; the drilling rate V of the drill rod may be measured with a drilling rate sensor connected to the drilling rig.
[0036] Optionally, in the device for determining rock stratum parameters provided in this embodiment, the determination module is further configured to determine the estimation formula of uniaxial compressive strength value R, of the rock; the method for determining the estimation formula of the uniaxial compressive strength value R, of the rock includes: taking the work value WF of the axial force of the drill rod when drilling through unit volume of rock, the work value Wm of the torque of the drill rod when drilling through unit volume of rock and the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock as independent variables, taking the uniaxial compressive strength value R, of the rock as a dependent variable, and obtaining the estimation formula with a linear regression method.
[0037] Optionally, the device for determining the rock stratum parameters provided in this embodiment further includes a data acquisitor configured to acquire test data when the drill rod drills through unit volume of rock; the data acquisitor includes a drilling rate sensor, a rotation speed sensor, a pressure sensor and a torque sensor; the drilling rate sensor is configured to acquire the drilling rate V; the rotation speed sensor is configured to acquire the rotation speed n of the drill rod; the pressure sensor is configured to acquire the thrust force F of the drill rod; the torque sensor is configured to acquire the torque M of the drill rod; the device further includes a display unit configured to display the test data acquired by the data acquisitor when the drill rod drills through unit volume of rock, the data acquired by the acquisition module, and the data determined by the determination module; the acquisition module is configured to receive the test data acquired by the data acquisitor when the drill rod drills through unit volume of rock and the corresponding uniaxial compressive strength value R,; the determination module processes the data information acquired by the acquisition module.
[0038] Those having ordinary skills in the art can easily understand that all or a part of the steps that implement the above embodiment of method can be accomplished with hardware related to programs and instructions. The program may be stored in computer readable storage media. The program executes the steps of the above embodiment of method when it is executed; the aforementioned storage media include ROM, RAM, magnetic disk or optical disk, etc., which can store program codes.
[0039] Some researchers have made relevant researches in the laboratory and used PDC drill bits with 60mm diameter to carry out drilling experiments. The drilling specimens are divided into 28 groups of mortar specimens different in strength and numbered as J1-J28 and 8 groups of sandstone specimens numbered as S1-S8. Table 1 shows 30 groups of experimental data.
[0040] Based on the data in the Table 1, an estimation model of uniaxial compressive strength value R, of rock is preset as formula (1) for multi-element linear regression, and an estimation formula of uniaxial compressive strength value R, of rock is obtained:
R = (5.82067e-6)* F +(5.4583e - 7)* +(7.78417e -7)* pFN - 4.72111. ArLcr 30Ar 2V 60I rV
Table 1 - Drilling Experiment Data
Type of Drilling rate, Rotation speed, Thrust specimen mm/min r/min. force, kN Ji 148.46 50 18.92 0.02 1.90 J2 177.89 100 14.43 0.01 2.00 M5 J3 185.07 100 12.98 0.01 1.94 J4 174.69 100 13.99 0.01 1.99 J5 124.8 50 17.01 0.03 2.37 M7.5 J6 87.94 100 7.30 0.02 2.58 J7 105.70 100 6.19 0.03 3.29 J9 81.52 50 17.05 0.03 6.70 103.09 100 12.22 0.02 7.20 M10 J1O J11 112.75 100 10.13 0.03 6.24 J12 130.77 100 14.67 0.03 6.99 J13 83.85 50 28.77 2.79 10.23 83.95 100 16.42 2.66 10.05 M15 J14 J15 132.05 100 22.34 2.15 10.54 J16 118.08 100 21.44 2.07 10.60 J17 83.91 50 45.43 2.23 23.54 M20 J18 83.65 100 29.89 2.33 30.81 J19 111.42 100 35.88 3.01 22.43 J21 83.46 50 43.66 3.44 21.66 M25 J22 84.23 100 26.15 1.20 27.8 J23 137.82 100 38.30 3.16 22.22 J25 84.64 50 41.49 0.85 27.71 M30 J26 82.88 100 25.99 0.51 35.21 J27 137.84 100 31.36 1.49 22.73 Sl 79.97 50 103.08 5.59 58.09 S3 84.07 150 40.69 2.65 59.95 Sandstone S5 84.90 250 28.00 2.15 61.91 S6 84.41 300 23.78 2.05 60.88 S7 110.72 100 65.20 5.48 51.41 S8 136.47 100 73.96 6.10 49.80
Claims (9)
- Claims 1. A method for determining rock stratum parameters, comprising:obtaining g test data groups acquired at different data acquisition times t during rock drilling, wherein the test data in the test data groups comprises work values WF of the axial force of the drill rod when drilling through unit volume of rock, work values WM of the torque of the drill rod when drilling through unit volume of rock, and work values Wf of the friction between the drill bit and borehole bottom when drilling through unit volume of rock;substituting the test data in the g test data groups into an estimation formula of uniaxial compressive strength value R, of the rock to determine g uniaxial compressive strength values Rei of the rock that respectively correspond to the g test data groups acquired at different data acquisition times t, wherein Rej is the uniaxial compressive strength value of the rock that corresponds to the ith test data group acquired in accordance with the sequence of acquisition times;forming a data sequence A= (R 1 , R,2 , Rc3 ,..., Rei,..., R,) from the determined g uniaxial compressive strength values Rci of the rock in accordance with the sequence of data acquisition times, and analyzing the data sequence A with K-means clustering analysis software to determine a final categorization result, wherein the final categorization result comprises a best number of categories ds and subsequences in each category; wherein the best number of categorizes d, is the total number of rock strata divided according to the difference in the uniaxial compressive strength value of the rock; if the eth subsequence is Ae = (Rca, Rcai,..., Rch,), the average value R of the uniaxial compressive strength of the ethlayer of rock strata is:1 b R =R, (1. b-a+1 =
- 2. The determination method according to claim 1, further comprising: obtaining an expression V(t) of the drilling rate V of the drill rod that changes with the data acquisition time t, wherein the drilling rate of the drill rod is the depth of the drill rod penetrating into the rock stratum per second, in unit of m/s; the method for obtaining the expression V(t) of the drilling rate V of the drill rod comprises: fitting the drilling rates V of different drill rods and the acquisition times t of the drilling rates V to obtain the expression V(t);determining an intervallT= t,., t6 ] of the data acquisition time t corresponding to the uniaxial compressive strength value of the rock in the eth subsequence Ae= (R,, Rca+,..., Rc );then the thickness of the eth layer of rock strata is:L= fbV(t) dt (2).
- 3. The determination method according to claim 1, whereinthe work value WF of the axial force of the drill rod when drilling through unit volume of rock is calculated according to a formula (3):F WF= 2 (3)the work value WM of the torque of the drill rod when drilling through unit volume of rock is calculated according to a formula (4):WM =M (4) 30A~r2Vthe work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock is calculated according to a formula (5):p Fn Wf =,=un (5) 60Y ArVwhere, F is the thrust force of the drill rod, in unit of N; A is a borehole enlargement coefficient, i.e., a ratio of the cross-sectional area of the borehole to the cross-sectional area of the drill rod; r is the bit diameter of the drill rod, in unit of m; M is the torque of the drill rod, in unit of N-m; n is the rotation speed of the drill rod, in unit of r/min.; V is the drilling rate of the drill rod, i.e., the penetration depth of the drill rod into the rock stratum per second, in unit of m/s; P is a coefficient of friction between the drill bit and the borehole bottom, and is determined to be 0.21.
- 4. The determination method according to claim 1, wherein the method for determining the estimation formula of the uniaxial compressive strength value R, of the rock comprises:taking the work value WF of the axial force of the drill rod when drilling through unit volume of rock, the work value WM of the torque of the drill rod when drilling through unit volume of rock and the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock as independent variables, taking the uniaxial compressive strength value R, of the rock as a dependent variable, and obtaining the estimation formula with a linear regression method.
- 5. A device for determining rock stratum parameters, comprising:an acquisition module configured to obtain g test data groups acquired at different data acquisition times t during rock drilling, wherein the test data in the test data groups comprises work values WF of the axial force of the drill rod when drilling through unit volume of rock, work values WM of the torque of the drill rod when drilling through unit volume of rock, and work values Wf of the friction between the drill bit and borehole bottom when drilling through unit volume of rock; a determination module configured to substitute the test data in the g test data groups into an estimation formula of uniaxial compressive strength value R, of the rock to determine g uniaxial compressive strength values Rj of the rock that respectively correspond to the g test data groups acquired at different data acquisition times t, wherein Rj is the uniaxial compressive strength value of the rock that corresponds to the ith test data group acquired in accordance with the sequence of acquisition times; the determination module is further configured to form a data sequence A= (R 1 , Rc 2 , R 3 ,..., Ri,..., Rg) from the determined g uniaxial compressive strength values Rj of the rock in accordance with the sequence of data acquisition times, and analyze the data sequence A with K-means clustering analysis software to determine a final categorization result, wherein the final categorization result comprises a best number of categories d, and subsequences in each category; wherein the best number of categorizes d, is the total number of rock strata divided according to the difference in the uniaxial compressive strength value of the rock; if the eth subsequence is A,= (R, Ra,+,..., Re), the average value R of the uniaxial compressive strength of the eth layer of rock strata is:I b R= R (1). b-a+1 j,
- 6. The determination device according to claim 5, whereinthe acquisition module is further configured to obtain an expression V(t) of the drilling rate V of the drill rod that changes with the data acquisition time t, wherein the drilling rate of the drill rod is the depth of the drill rod penetrating into the rock stratum per second, in unit of m/s; the method for obtaining the expression V(t) of the drilling rate V of the drill rod comprises: fitting the drilling rates V of different drill rods and the acquisition times t of the drilling rates V to obtain the expression V(t);determining an interval T,= [t, , tb ] of the data acquisition time t corresponding to the uniaxial compressive strength value of the rock in the eth subsequence A,= (Rca, R.,..., Ry )then the thickness of the eth layer of rock strata is:L= bV(t) dt (2).
- 7. The determination device according to claim 5, wherein the acquisition module is further configured to obtain the work value WF of the axial force of the drill rod when drilling through unit volume of rock according to a formula (3):F WF 2 (3)the acquisition module is further configured to obtain the work value WM of the torque of the drill rod when drilling through unit volume of rock according to a formula (4):Mn WM 3 O>{r2 (4) 30Ar2Vthe acquisition module is further configured to obtain the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock according to a formula (5):W= Fn (5) 6OA2rVwhere, F is the thrust force of the drill rod, in unit of N; is a borehole enlargement coefficient, i.e., a ratio of the cross-sectional area of the borehole to the cross-sectional area of the drill rod; r is the bit diameter of the drill rod, in unit of m; M is the torque of the drill rod, in unit of N-m; n is the rotation speed of the drill rod, in unit of r/min.; V is the drilling rate of the drill rod, i.e., the penetration depth of the drill rod into the rock stratum per second, in unit of m/s; p is a coefficient of friction between the drill bit and the borehole bottom, and is determined to be 0.21.
- 8. The determination device according to claim 5, wherein the determination module is further configured to determine the estimation formula of the uniaxial compressive strength value R, of the rock;the method for determining the estimation formula of the uniaxial compressive strength value R, of the rock comprises: taking the work value WF of the axial force of the drill rod when drilling through unit volume of rock, the work value WM of the torque of the drill rod when drilling through unit volume of rock and the work value Wf of the friction between the drill bit and the borehole bottom when drilling through unit volume of rock as independent variables, taking the uniaxial compressive strength value R, of the rock as a dependent variable, and obtaining the estimation formula with a linear regression method.
- 9. The determination device according to claim 5, further comprising: a data acquisitor configured to acquire the test data when the drill rod drills through unit volume of rock; the data acquisitor comprises a drilling rate sensor, a rotation speed sensor, a pressure sensor and a torque sensor; the drilling rate sensor is configured to acquire the drilling rate V; the rotation speed sensor is configured to acquire the rotation speed n of the drill rod; the pressure sensor is configured to acquire the thrust force F of the drill rod; the torque sensor is configured to acquire the torque M of the drill rod; a display unit configured to display the test data acquired by the data acquisitor when the drill rod drills through unit volume of rock, the data acquired by the acquisition module, and the data determined by the determination module; the acquisition module is configured to receive the test data acquired by the data acquisitor when the drill rod drills through unit volume of rock and the corresponding uniaxial compressive strength value R,; the determination module processes the data information acquired by the acquisition module.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910255972.8 | 2019-04-01 | ||
| CN201910255972.8A CN110130883A (en) | 2019-04-01 | 2019-04-01 | The determination method and device of formation parameters |
| PCT/CN2019/104107 WO2020199495A1 (en) | 2019-04-01 | 2019-09-03 | Method and device for determining rock stratum parameters |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| AU2019439997A1 true AU2019439997A1 (en) | 2021-01-07 |
Family
ID=67569123
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| AU2019439997A Abandoned AU2019439997A1 (en) | 2019-04-01 | 2019-09-03 | Method and device for determining rock stratum parameters |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN110130883A (en) |
| AU (1) | AU2019439997A1 (en) |
| WO (1) | WO2020199495A1 (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110130883A (en) * | 2019-04-01 | 2019-08-16 | 中国矿业大学 | The determination method and device of formation parameters |
| CN110486007B (en) * | 2019-08-29 | 2023-02-03 | 武汉长盛煤安科技有限公司 | In-situ testing device and method for mechanical parameters of coal mine surrounding rock while drilling |
| CN111075424B (en) * | 2019-12-25 | 2022-11-18 | 中国石油大学(华东) | Method for correcting measurement results of measurement-while-drilling parameters |
| CN111706322B (en) * | 2020-07-17 | 2023-03-21 | 中国铁建重工集团股份有限公司 | Rock drilling response prediction method and prediction system |
| CN111855961B (en) * | 2020-07-24 | 2021-10-26 | 中南大学 | Rock mass drilling quality detection method, drilling machine, server and storage medium |
| CN112412394B (en) * | 2020-11-11 | 2022-12-06 | 安徽理工大学 | Drilling layered filling method |
| CN113482599A (en) * | 2021-07-08 | 2021-10-08 | 中国科学院武汉岩土力学研究所 | Portable in-situ test and identification data acquisition system and equipment |
| CN116485225B (en) * | 2023-03-15 | 2023-11-10 | 西南交通大学 | Automatic acquisition method and system for BQ value of surrounding rock in construction stage based on drilling parameters |
| CN116499827A (en) * | 2023-04-18 | 2023-07-28 | 中国石油大学(北京) | Multi-lithology interbed artificial rock sample preparation device and method |
| CN117288587B (en) * | 2023-11-24 | 2024-02-20 | 中国矿业大学(北京) | Rock tensile strength while drilling test method and system |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6672386B2 (en) * | 2002-06-06 | 2004-01-06 | Baker Hughes Incorporated | Method for in-situ analysis of formation parameters |
| EP1836644B1 (en) * | 2004-11-04 | 2013-10-23 | Baker Hughes Incorporated | Multiscale multidimensional well log data inversion and deep formation imaging method |
| CN102587869A (en) * | 2005-05-25 | 2012-07-18 | 地质力学国际公司 | Methods and devices for analyzing and controlling the propagation of waves in a borehole generated by water hammer |
| JP2007327880A (en) * | 2006-06-08 | 2007-12-20 | Tokyo Electric Power Co Inc:The | Instrument for measuring compression strength of concrete, and method of measuring compression strength of concrete |
| WO2008076130A1 (en) * | 2006-12-15 | 2008-06-26 | Halliburton Energy Services, Inc. | Antenna coupling component measurement tool having rotating antenna configuration |
| WO2012125139A1 (en) * | 2011-03-11 | 2012-09-20 | Landmark Graphics Corporation | Methods and systems of estimating formation parameters |
| CN104217124A (en) * | 2014-09-15 | 2014-12-17 | 天津大学 | TBM (Tunnel Boring Machine) construction surrounding rock classification method depending on engineering sample data |
| CN108663269B (en) * | 2018-03-30 | 2019-12-24 | 山东大学 | Underground engineering surrounding rock digital drilling zoning method based on equivalent rock mass strength |
| CN108519284B (en) * | 2018-03-30 | 2019-09-27 | 山东大学 | Underground engineering wall rock digital rig real-time graded method |
| CN108760487A (en) * | 2018-04-16 | 2018-11-06 | 西安科技大学 | Deep cavern country rock subregion failure evolvement analysis method based on rockbolt stress analysis |
| CN110130883A (en) * | 2019-04-01 | 2019-08-16 | 中国矿业大学 | The determination method and device of formation parameters |
-
2019
- 2019-04-01 CN CN201910255972.8A patent/CN110130883A/en active Pending
- 2019-09-03 WO PCT/CN2019/104107 patent/WO2020199495A1/en active Application Filing
- 2019-09-03 AU AU2019439997A patent/AU2019439997A1/en not_active Abandoned
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020199495A1 (en) | 2020-10-08 |
| RU2762675C1 (en) | 2021-12-21 |
| CN110130883A (en) | 2019-08-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2019439997A1 (en) | Method and device for determining rock stratum parameters | |
| CN108663269B (en) | Underground engineering surrounding rock digital drilling zoning method based on equivalent rock mass strength | |
| RU2723805C9 (en) | Method and computer system for control of drilling of the wells | |
| NO333278B1 (en) | Method and data processing system for establishing a model of underground formations during drilling | |
| CN103244101A (en) | Method for evaluating quality of rock stratum while drilling | |
| CN103308946A (en) | Tunnel advance geology forecast method based on blast hole drilling information | |
| CN116522692B (en) | Underground engineering surrounding rock structural feature in-situ detection and classification method | |
| CN104374827A (en) | Measuring method of anisotropy coefficient of transverse isotropic rock in-situ dynamic elasticity modulus | |
| NO344538B1 (en) | Procedure and computer-readable medium for prediction of production zone | |
| CN112664174A (en) | Tunnel surrounding rock grade determination method and system based on multiple drill holes | |
| CN113074695A (en) | Tunnel stability interpretation method based on tunnel face three-dimensional deformation monitoring | |
| CN104153768A (en) | Granite reservoir stratum reservoir performance evaluation method | |
| CN105093349A (en) | Method for actually measuring growth and development rule of crack in tunnel roof | |
| CN103590828B (en) | A kind of well logging dck index method evaluates the method for strata pressure | |
| Ibrahim et al. | Integration of pressure-transient and fracture area for detecting unconventional wells interference | |
| CN112948924B (en) | Method for determining height of near-unconsolidated formation mining water-guiding fracture zone based on overburden structure | |
| CN103195425A (en) | System for rapidly measuring in-situ wall rock strength of coal mine tunnel | |
| CN116451013B (en) | Deep stratum rock in-situ drillability grade value prediction method | |
| CN105064987A (en) | Interpretation and evaluation method for water layer recognition by logging while drilling Q parameter | |
| CN109736771B (en) | Surrounding rock analysis method and system based on drilling parameters of drill jumbo | |
| CN110107274B (en) | Real-time online monitoring system and monitoring method for TBM (Tunnel boring machine) jumbolter based on hydraulic system | |
| CN104727815A (en) | Real-time well drilling formation correction method and device | |
| CN116432855A (en) | Tunnel collapse condition prediction method based on tunneling data | |
| US10677959B2 (en) | Method of and apparatus for calculating UCS and CCS | |
| CN104314560A (en) | Comparing method and comparing device for non-proportional stratum before window entering of horizontal well |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| DA3 | Amendments made section 104 |
Free format text: THE NATURE OF THE AMENDMENT IS: AMEND THE NAME OF THE INVENTOR TO READ AI, DECHUN; ZHENG, XIGUI; LIU, CANCAN; YANG, JUNWEI; WEI, ZHONGJU AND LI, PENG |
|
| HB | Alteration of name in register |
Owner name: CHINA UNIVERSITY OF MINING AND TECHNOLOGY Free format text: FORMER NAME(S): CHINA UNIVERSITY OF MINING AND TECHNOLOGY Owner name: LIUPANSHUI NORMAL UNIVERSITY Free format text: FORMER NAME(S): CHINA UNIVERSITY OF MINING AND TECHNOLOGY |
|
| MK5 | Application lapsed section 142(2)(e) - patent request and compl. specification not accepted |