WO2020199495A1 - Method and device for determining rock stratum parameters - Google Patents

Method and device for determining rock stratum parameters Download PDF

Info

Publication number
WO2020199495A1
WO2020199495A1 PCT/CN2019/104107 CN2019104107W WO2020199495A1 WO 2020199495 A1 WO2020199495 A1 WO 2020199495A1 CN 2019104107 W CN2019104107 W CN 2019104107W WO 2020199495 A1 WO2020199495 A1 WO 2020199495A1
Authority
WO
WIPO (PCT)
Prior art keywords
rock
drilling
drill pipe
compressive strength
uniaxial compressive
Prior art date
Application number
PCT/CN2019/104107
Other languages
French (fr)
Chinese (zh)
Inventor
刘灿灿
郑西贵
艾德春
刘洪洋
Original Assignee
中国矿业大学
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 中国矿业大学 filed Critical 中国矿业大学
Priority to AU2019439997A priority Critical patent/AU2019439997A1/en
Priority to RU2020141122A priority patent/RU2762675C9/en
Publication of WO2020199495A1 publication Critical patent/WO2020199495A1/en

Links

Images

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

Definitions

  • the invention relates to a method and device for determining rock parameters, belonging to the technical field of geological prospecting engineering.
  • the present invention provides a method and device for determining rock formation parameters.
  • the method for determining rock formation parameters includes obtaining g test data sets collected at different data collection times t during the process of drilling into the rock formation, and the test data in the test data set includes the data when drilling a unit volume of rock Drill pipe axial force work value w F , drill pipe torque work value w M when drilling unit volume of rock, drill bit and hole bottom friction work value w f when drilling unit volume of rock; put the g test data groups
  • the test data are respectively substituted into the estimation formula of the rock uniaxial compressive strength value R c to determine the g rock uniaxial compressive strength values R corresponding to the g test data sets collected at different data collection times t ci , where R ci is the uniaxial compressive strength value of the rock corresponding to the i-th test data set obtained in the order of collection time; the determined uniaxial compressive strength values of g rocks R ci are formed in the order of data collection
  • the data sequence A of the uniaxial compressive strength value of a rock (R )
  • the method for determining rock formation parameters provided by the present invention further includes: calculating the work value w F of the drill pipe axial force when drilling the unit volume of rock according to formula (3):
  • the method for determining the estimation formula of the rock uniaxial compressive strength value R c in the method for determining rock formation parameters provided by the present invention includes: using the drill pipe axial force value w F when drilling a unit volume of rock, The drill pipe torque work value w M when drilling a unit volume of rock, the drill bit and the hole bottom friction work value w f when drilling a unit volume of rock are the independent variables, and the rock uniaxial compressive strength value R c is the dependent variable. Obtained by linear regression method.
  • the present invention also provides a rock formation parameter determination device, including an acquisition module for acquiring g test data sets collected at different data acquisition times t during the process of drilling into the rock formation, the test data sets
  • the test data in includes the drill pipe axial work value w F when drilling a unit volume of rock, the drill pipe torque work value w M when drilling a unit volume of rock, and the friction work value between the drill bit and the bottom of the hole when drilling a unit volume of rock w f ;
  • It also includes a determination module for substituting the test data in the g test data groups into the estimation formula of the rock uniaxial compressive strength value R c to determine the different data collection time t
  • the uniaxial compressive strength values of g rocks corresponding to the g collected test data sets R ci where R ci is the uniaxial compressive strength of the rock corresponding to the i-th test data set obtained in the order of collection time
  • the determination module is also used to determine the uniaxial compressive strength values of g rocks R
  • the acquisition module in the rock formation parameter determination device is also used to acquire the expression V(t) of the drilling speed V of the drill pipe that varies with the data acquisition time t,
  • the advance rate is the depth at which the drill pipe enters the rock formation per second, in m/s;
  • the acquisition module is further configured to acquire the drill pipe axial force work value w F when drilling a unit volume of rock according to formula (3):
  • the obtaining module is also used to obtain the drill pipe torque work value w M when drilling a unit volume of rock according to formula (4):
  • the acquisition module is also used to acquire the friction work value w f between the drill bit and the hole bottom when drilling the unit volume of rock according to formula (5):
  • means for determining formation parameters in the present invention provides the determination module is further for determining the estimation formula uniaxial compressive strength value R c; and estimating the uniaxial compressive strength value R c is
  • the method for determining the formula includes: the work value w F of the drill pipe axial force when drilling a unit volume of rock, the work value w M of the drill pipe torque when drilling a unit volume of rock, and the friction between the drill bit and the bottom of the hole when drilling a unit volume of rock.
  • the work value w f is the independent variable
  • the rock uniaxial compressive strength value R c is the dependent variable, which is obtained by linear regression.
  • the rock formation parameter determination device further includes a data collector for collecting test data when the drill pipe is drilled into a unit volume of rock;
  • the data collector includes a drilling speed sensor and a rotation speed sensor , Pressure sensor and torque sensor;
  • the drilling speed sensor is used to collect the drilling speed V;
  • the rotation speed sensor is used to collect the rotation speed n of the drill rod;
  • the pressure sensor is used to collect the thrust F of the drill rod;
  • the torque sensor It is used to collect the torque M of the drill pipe; it also includes a display for displaying the test data collected by the data collector when the drill pipe is drilled into a unit volume of rock, the data obtained by the acquisition module and the data determined by the determination module;
  • the acquisition module is used to receive the test data and the corresponding uniaxial compressive strength value R c when the drill pipe is drilled into a unit volume of rock collected by the data collector;
  • the determination module processes the data information acquired by the acquisition module.
  • the method and device for determining rock formation parameters provided by the present invention are based on the drill rod axial force when drilling a unit volume of rock in a plurality of test data sets collected at different data collection times during drilling into the rock formation.
  • the work value, the drill pipe torque work value when drilling a unit volume of rock, and the drill bit and the hole bottom friction work value when drilling a unit volume of rock are substituted into the estimation formula of the rock uniaxial compressive strength value R c to obtain multiple rock units.
  • Axial compressive strength value, and K-means cluster analysis software is used to determine the sub-sequence of rock uniaxial compressive strength values in different rock formation classifications.
  • Figure 1 is a schematic flow chart of the method for determining rock formation parameters of this embodiment
  • Figure 2 is a schematic diagram of the connection of the rock formation parameter determination device of this embodiment.
  • the method for determining rock formation parameters provided in this embodiment includes:
  • test data in the test data set include the value of drill pipe axial force w F when drilling a unit volume of rock, and drilling
  • the experimental device in the process of drilling the rock formation by the rig, can directly or indirectly obtain the drill pipe axial force work value w F per unit volume of rock, and the drill pipe torque work value w when drilling unit volume of rock.
  • M Experimental data such as the friction work value w f between the drill bit and the bottom of the hole when drilling a unit volume of rock; when the experimental device is acquiring these data, the time when the data is acquired can be recorded at the same time, and the time for the data acquisition can be It is synchronized.
  • g test data sets represent multiple test data sets, which are not a specific limitation on the number of data sets, but are merely an expression of the number of acquired test data sets for the convenience of subsequent expressions.
  • the formula for estimating the uniaxial compressive strength value of rock R c has been determined in advance, and the formula for estimating the uniaxial compressive strength value of rock R c can be directly used by reference formulas already determined in other documents, or It is obtained by linear regression on other experimental data.
  • the mean value R of uniaxial compressive strength of e-layer rock is:
  • the lithology of the rock formations changes, the g rock uniaxial compressive strength values R ci obtained in this process have time series themselves, in order to obtain continuous
  • the lithology of the rock formations when analyzing these data, also has time sequence constraints.
  • the subsequence is a certain segment of the data sequence A, and there is no overlap between multiple subsequences.
  • the data sequence A is divided into d s subsequences, and after the multiple subsequences are arranged in chronological order, the continuity of the data and the data sequence A can be exactly the same.
  • K-means algorithm is a hard clustering algorithm. It is a representative of a typical prototype-based objective function clustering method. It uses a certain distance from a data point to the prototype as the optimized objective function, and iterative operation is obtained by the method of seeking extreme values of the function. Adjustment rules.
  • the algorithm uses the error square sum criterion function as the clustering criterion function.
  • Software with K-means algorithm function commonly used R software, Matlab. The specific process of cluster analysis is completed by the existing software. This embodiment does not describe the process and principle of cluster analysis in detail, and only exemplifies the raw data that needs to be input and the type of data to be obtained.
  • the method for determining rock formation parameters provided in this embodiment further includes the following steps:
  • the drilling speed of the drill pipe in the rock formation will change.
  • the drilling speed sensor can record the drilling speed value of the drill pipe every certain time, such as a few milliseconds. According to the experiment Through non-linear fitting, the multiple drilling speed values corresponding to different times in a certain time period collected in the data collection time can be obtained by the expression V(t) of the drilling speed V of the drill pipe changing with the data collection time t.
  • step S1 may further include: calculating the work value w F of the drill rod axial force when drilling a unit volume of rock according to formula (3):
  • test devices can directly obtain the work value w F of the drill pipe axial force when drilling a unit volume of rock, the work value w M of the drill pipe torque when drilling a unit volume of rock, and the bit and the value of the drill bit when drilling a unit volume of rock.
  • the friction work value w f at the bottom of the hole in practice, it can also be obtained by indirect measurement.
  • step S2 may further include: in the method for determining rock formation parameters provided in this embodiment, the method for determining the estimation formula of the rock uniaxial compressive strength value R c includes: The drill pipe axial force work value w F in volume rock, the drill pipe torque work value w M when drilling a unit volume of rock, the drill bit and the hole bottom friction work value w f when drilling a unit volume of rock are the independent variables.
  • the uniaxial compressive strength value R c is the dependent variable and is obtained by linear regression.
  • the linear regression method is a data processing method commonly used in scientific research, that is, a method of obtaining an empirical formula by obtaining reliable experimental data in advance, and then using the empirical formula to estimate the result.
  • the software for multiple linear regression can use SPSS (Statistical Product and Service Solutions), Microsoft Excel and other software that can perform data processing.
  • this embodiment also provides a rock formation parameter determination device, including an acquisition module, which is used to acquire g test data sets collected at different data acquisition times t during drilling into the rock formation.
  • the test data in the test data group includes the work value w F of the drill pipe axial force when drilling a unit volume of rock, the work value w M of the drill pipe torque when drilling a unit volume of rock, the drill bit and the bottom of the hole when drilling a unit volume of rock Friction work value w f ;
  • It also includes a determination module, which is used to substitute the test data in the g test data groups into the estimation formula of the rock uniaxial compressive strength value R c to determine the collection at different data collection times t
  • the obtained g test data sets correspond to g rock uniaxial compressive strength values R ci , where R ci is the rock uniaxial compressive strength value corresponding to the i-th test data set obtained in the order of collection time; determine;
  • the module is also used to convert the determined uniaxial compressive
  • the acquisition module in the rock formation parameter determination device is also used to acquire the expression V(t) of the drilling speed V of the drill pipe that varies with the data collection time t, and the drilling speed of the drill pipe is The depth at which the drill pipe enters the rock formation per second, in m/s;
  • the method for obtaining the expression V(t) of the drilling speed V of the drill pipe includes: the difference between the drilling speed V and the drilling speed V of different drill pipes
  • the interval of the collection time t T e [t a , t b ]; then the thickness of the e-th layer is:
  • the acquisition module is further used to acquire the value w F of the axial force of the drill pipe when drilling a unit volume of rock according to formula (3):
  • the obtaining module is also used to obtain the work value w M of drill pipe torque when drilling a unit volume of rock according to formula (4):
  • the acquisition module is also used to acquire the friction work value w f between the drill bit and the hole bottom when drilling a unit volume of rock according to formula (5):
  • the thrust F of the drill pipe can be obtained by subtracting the gravity of the drill rig from the pressure of the rig in contact with the ground; the reaming coefficient ⁇ is the ratio of the area of the cross section of the drill pipe after the drill pipe is formed;
  • the torque M of the rod can be measured by a torque sensor connected to the drill; the speed n of the drill rod can be measured by a speed sensor connected to the drill; the drilling speed V of the drill rod can be measured by a drilling speed sensor connected to the drill.
  • the independent variable, taking the rock uniaxial compressive strength value R c as the dependent variable, is obtained by linear regression.
  • the rock formation parameter determination device further includes a data collector, which is used to collect test data when the drill pipe is drilled into a unit volume of rock;
  • the data collector includes a drilling speed sensor, a rotation speed sensor, and a pressure sensor And torque sensor;
  • the drilling speed sensor is used to collect the drilling speed V;
  • the speed sensor is used to collect the rotation speed n of the drill rod;
  • the pressure sensor is used to collect the thrust F of the drill rod;
  • the torque sensor is used to collect the torque M of the drill rod;
  • the display is used to display the test data collected by the data collector when the drill pipe is drilled into a unit volume of rock, the data obtained by the acquisition module and the data determined by the determination module;
  • the acquisition module is used to receive the drill pipe drilling unit collected by the data collector
  • the determination module processes the data information acquired by the acquisition module.
  • a person of ordinary skill in the art can understand that all or part of the steps of the foregoing method embodiments can be implemented by a program instructing relevant hardware.
  • the aforementioned program can be stored in a computer readable storage medium.
  • the steps including the foregoing method embodiments are executed; and the foregoing storage medium includes: ROM, RAM, magnetic disk, or optical disk and other media that can store program codes.
  • the estimation model of rock uniaxial compressive strength value R c is preset as formula (1) to perform multiple linear regression to obtain the estimation formula of rock uniaxial compressive strength value R c :

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 Rc the axial force working values of a drill rod during drilling into rock in unit volume, torque working values of the drill rod during drilling into the rock in unit volume, and friction working values of a drill bit and hole bottom during drilling into the rock in unit volume in a plurality of obtained test data sets collected in different data collecting time in a rock stratum drilling process, so as to obtain multiple rock uniaxial compressive strength values; determining subsequences composed of the rock uniaxial compressive strength values in different rock stratum classifications by means of K-means cluster analysis software, and solving a mean value of the rock uniaxial compressive strength values in the subsequences so as to obtain a mean uniaxial compressive strength value of a current rock stratum and to further determine the thicknesses of all the rock stratums. The method is used to continuously detect the lithology of sinking and driving engineering, and thus saves the labor and time costs.

Description

岩层参数的确定方法及装置Method and device for determining rock parameters 技术领域Technical field
本发明涉及一种岩层参数的确定方法及装置，属于地质勘探工程技术领域。The invention relates to a method and device for determining rock parameters, belonging to the technical field of geological prospecting engineering.
背景技术Background technique
在煤矿、隧道等地下工程中，岩层情况对工程的安全有较大影响。目前，主要通过岩石取芯、钻孔窥视等了解岩层信息，但上述方法影响工程进展速度，且费时费力，不能够实时了解工程所处岩层环境。In underground projects such as coal mines and tunnels, the condition of the rock formation has a greater impact on the safety of the project. At present, the rock formation information is mainly understood through rock cores and drilling holes. However, the above methods affect the progress of the project and are time-consuming and laborious, and it is impossible to understand the rock formation environment of the project in real time.
目前，石油行业的随钻测井技术相对成熟，但局限性较大，在煤矿等地下工程中还没广泛应用。也有学者提出了基于钻孔岩屑的地层探测法，但岩屑收集困难，且对于岩性易溶于水的岩层，不能够收集到该岩层的岩屑，导致对地层结构的判断不连续，且该方法不能对岩石的单轴抗压强度进行预测。At present, the logging-while-drilling technology in the petroleum industry is relatively mature, but it is relatively limited and has not been widely used in underground projects such as coal mines. Some scholars have also proposed a formation detection method based on drilling cuttings, but cuttings collection is difficult, and for rock formations that are easily soluble in water, it is impossible to collect the cuttings of the rock formations, which leads to discontinuous judgments of the formation structure. Moreover, this method cannot predict the uniaxial compressive strength of rocks.
因此，许多学者提出了基于锚杆钻机的随钻参数实时识别岩层的方法，但目前国内还未找到有效方法使该技术得到广泛应用。Therefore, many scholars have proposed a method for real-time identification of rock formations based on the parameters of the bolt drill while drilling, but there is no effective method to make this technology widely used in China.
发明内容Summary of the invention
为了解决现有技术中存在的上述技术问题，本发明提供了一种岩层参数的确定方法及装置。In order to solve the above technical problems in the prior art, the present invention provides a method and device for determining rock formation parameters.
本发明提供的岩层参数的确定方法，包括获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，所述试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f；将所述g个试验数据组中的试验数据分别代入岩石单轴抗压强度值R c的估算公式，以确定在不同的数据采集时间t采集到的所述g个试验数据组分别对应的g个岩石单轴抗压强度值R ci，其中R ci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；将确定的g个岩石单轴抗压强度值R ci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(R c1，R c2，R c3，……，R ci，……，R cg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；所述最终分类结果包括最佳分类数d s和每一分类的子序列；其中最佳分类数d s即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为A e＝(R ca，R ca+1，……，R cb)时，则第e层岩层的单轴抗压强度均值R为： The method for determining rock formation parameters provided by the present invention includes obtaining g test data sets collected at different data collection times t during the process of drilling into the rock formation, and the test data in the test data set includes the data when drilling a unit volume of rock Drill pipe axial force work value w F , drill pipe torque work value w M when drilling unit volume of rock, drill bit and hole bottom friction work value w f when drilling unit volume of rock; put the g test data groups The test data are respectively substituted into the estimation formula of the rock uniaxial compressive strength value R c to determine the g rock uniaxial compressive strength values R corresponding to the g test data sets collected at different data collection times t ci , where R ci is the uniaxial compressive strength value of the rock corresponding to the i-th test data set obtained in the order of collection time; the determined uniaxial compressive strength values of g rocks R ci are formed in the order of data collection The data sequence A of the uniaxial compressive strength value of a rock = (R c1 , R c2 , R c3 ,..., R ci ,..., R cg ), use the K-means clustering analysis software to obtain the data sequence A Determine the final classification result; the final classification result includes the best classification number d s and the sub-sequence of each classification; the best classification number d s is the total number of layers divided by the rock uniaxial compressive strength value ; When the e-th sub-sequence is A e = (R ca , R ca+1 ,..., R cb ), then the uniaxial compressive strength average R of the e-th layer is:
Figure PCTCN2019104107-appb-000001
Figure PCTCN2019104107-appb-000001
可选地，本发明提供的岩层参数的确定方法还包括获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，所述钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；所述钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；确定第e个子序列A e＝(R ca，R ca+1，……，R cb)中岩石单轴抗压强度值所对应的数据采集时间t的区间T e＝[t a，t b]；则第e层岩层的厚度为： Optionally, the method for determining rock formation parameters provided by the present invention further includes obtaining the expression V(t) of the drilling speed V of the drill pipe that varies with the data collection time t, where the drilling speed of the drill pipe is The depth into the rock formation in seconds, in m/s; the method of obtaining the expression V(t) of the drilling speed V of the drill pipe includes: collecting the drilling speed V and the drilling speed V of different drill pipes Time t obtains the expression V(t) through data fitting; determine the data collection corresponding to the uniaxial compressive strength of the rock in the e-th subsequence A e = (R ca , R ca+1 ,..., R cb ) The interval of time t T e =[t a , t b ]; then the thickness of the e-th layer is:
Figure PCTCN2019104107-appb-000002
Figure PCTCN2019104107-appb-000002
可选地，本发明提供的岩层参数的确定方法还包括：根据公式(3)计算所述钻进单位体积岩石时的 钻杆轴力做功值w F： Optionally, the method for determining rock formation parameters provided by the present invention further includes: calculating the work value w F of the drill pipe axial force when drilling the unit volume of rock according to formula (3):
Figure PCTCN2019104107-appb-000003
Figure PCTCN2019104107-appb-000003
根据公式(4)计算所述钻进单位体积岩石时的钻杆扭矩做功值w M： According to formula (4), calculate the work value w M of the drill pipe torque when drilling a unit volume of rock:
Figure PCTCN2019104107-appb-000004
Figure PCTCN2019104107-appb-000004
根据公式(5)计算所述钻进单位体积岩石时的钻头与孔底摩擦做功值w f： Calculate the friction work value w f between the drill bit and the bottom of the hole when drilling the unit volume of rock according to formula (5):
Figure PCTCN2019104107-appb-000005
Figure PCTCN2019104107-appb-000005
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。In the formula, F is the thrust of the drill pipe, in N; λ is the hole expansion coefficient, that is, the ratio of the area of the cross section of the drill hole to the area of the cross section of the drill pipe; r is the drill bit diameter of the drill pipe, in m ; M is the torque of the drill pipe, in N·m; n is the rotation speed of the drill pipe, in r/min; V is the drilling speed of the drill pipe, which is the depth of the drill pipe entering the rock formation per second, in m /s; μ is the friction coefficient between the drill bit and the bottom of the hole, and the value of μ is 0.21.
可选地，本发明提供的岩层参数的确定方法中所述岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。 Optionally, the method for determining the estimation formula of the rock uniaxial compressive strength value R c in the method for determining rock formation parameters provided by the present invention includes: using the drill pipe axial force value w F when drilling a unit volume of rock, The drill pipe torque work value w M when drilling a unit volume of rock, the drill bit and the hole bottom friction work value w f when drilling a unit volume of rock are the independent variables, and the rock uniaxial compressive strength value R c is the dependent variable. Obtained by linear regression method.
此外，本发明还提供了岩层参数的确定装置，包括获取模块，所述获取模块用于获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，所述试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f；还包括确定模块，所述确定模块用于将所述g个试验数据组中的试验数据分别代入岩石单轴抗压强度值R c的估算公式，以确定在不同的数据采集时间t采集到的所述g个试验数据组分别对应的g个岩石单轴抗压强度值R ci，其中R ci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；所述确定模块还用于将确定的g个岩石单轴抗压强度值R ci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(R c1，R c2，R c3，……，R ci，……，R cg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；所述最终分类结果包括最佳分类数d s和每一分类的子序列；其中最佳分类数d s即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为A e＝(R ca，R ca+1，……，R cb)时，则第e层岩层的单轴抗压强度均值R为： In addition, the present invention also provides a rock formation parameter determination device, including an acquisition module for acquiring g test data sets collected at different data acquisition times t during the process of drilling into the rock formation, the test data sets The test data in includes the drill pipe axial work value w F when drilling a unit volume of rock, the drill pipe torque work value w M when drilling a unit volume of rock, and the friction work value between the drill bit and the bottom of the hole when drilling a unit volume of rock w f ; It also includes a determination module for substituting the test data in the g test data groups into the estimation formula of the rock uniaxial compressive strength value R c to determine the different data collection time t The uniaxial compressive strength values of g rocks corresponding to the g collected test data sets R ci , where R ci is the uniaxial compressive strength of the rock corresponding to the i-th test data set obtained in the order of collection time The determination module is also used to determine the uniaxial compressive strength values of g rocks R ci according to the time sequence of data collection to form a rock uniaxial compressive strength data sequence A=(R c1 , R c2 , R c3 , …, R ci , …, R cg ), the data sequence A is passed through K-means cluster analysis software to determine the final classification result; the final classification result includes the best classification number d s and each The sub-sequence of classification; where the best classification number d s is the total number of layers divided by the rock uniaxial compressive strength value; when the e-th sub-sequence is A e = (R ca , R ca+1 ,... …, R cb ), then the average uniaxial compressive strength R of the e-th layer is:
Figure PCTCN2019104107-appb-000006
Figure PCTCN2019104107-appb-000006
可选地，本发明提供的岩层参数的确定装置中所述获取模块还用于获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，所述钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；所述钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t 通过数据拟合得到表达式V(t)；确定第e个子序列A e＝(R ca，R ca+1，……，R cb)中岩石单轴抗压强度值所对应的数据采集时间t的区间T e＝[t a，t b]；则第e层岩层的厚度为： Optionally, the acquisition module in the rock formation parameter determination device provided by the present invention is also used to acquire the expression V(t) of the drilling speed V of the drill pipe that varies with the data acquisition time t, The advance rate is the depth at which the drill pipe enters the rock formation per second, in m/s; the method for obtaining the expression V(t) of the drill pipe’s drilling velocity V includes: comparing the drilling velocity V of different drill pipes with The acquisition time t of the drilling speed V obtains the expression V(t) through data fitting; determine the uniaxial compressive strength of the rock in the e-th subsequence A e = (R ca , R ca+1 ,..., R cb ) The interval of the data collection time t corresponding to the value T e =[t a , t b ]; then the thickness of the e-th layer is:
Figure PCTCN2019104107-appb-000007
Figure PCTCN2019104107-appb-000007
可选地，本发明提供的岩层参数的确定装置中，所述获取模块还用于根据公式(3)获取所述钻进单位体积岩石时的钻杆轴力做功值w F： Optionally, in the rock formation parameter determination device provided by the present invention, the acquisition module is further configured to acquire the drill pipe axial force work value w F when drilling a unit volume of rock according to formula (3):
Figure PCTCN2019104107-appb-000008
Figure PCTCN2019104107-appb-000008
所述获取模块还用于根据公式(4)获取所述钻进单位体积岩石时的钻杆扭矩做功值w M： The obtaining module is also used to obtain the drill pipe torque work value w M when drilling a unit volume of rock according to formula (4):
Figure PCTCN2019104107-appb-000009
Figure PCTCN2019104107-appb-000009
所述获取模块还用于根据公式(5)获取所述钻进单位体积岩石时的钻头与孔底摩擦做功值w f： The acquisition module is also used to acquire the friction work value w f between the drill bit and the hole bottom when drilling the unit volume of rock according to formula (5):
Figure PCTCN2019104107-appb-000010
Figure PCTCN2019104107-appb-000010
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。In the formula, F is the thrust of the drill pipe, in N; λ is the hole expansion coefficient, that is, the ratio of the area of the cross section of the drill hole to the area of the cross section of the drill pipe; r is the drill bit diameter of the drill pipe, in m ; M is the torque of the drill pipe, in N·m; n is the rotation speed of the drill pipe, in r/min; V is the drilling speed of the drill pipe, which is the depth of the drill pipe entering the rock formation per second, in m /s; μ is the friction coefficient between the drill bit and the bottom of the hole, and the value of μ is 0.21.
可选地，本发明提供的岩层参数的确定装置中所述确定模块还用于确定所述岩石单轴抗压强度值R c的估算公式；所述岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。 Alternatively, means for determining formation parameters in the present invention provides the determination module is further for determining the estimation formula uniaxial compressive strength value R c; and estimating the uniaxial compressive strength value R c is The method for determining the formula includes: the work value w F of the drill pipe axial force when drilling a unit volume of rock, the work value w M of the drill pipe torque when drilling a unit volume of rock, and the friction between the drill bit and the bottom of the hole when drilling a unit volume of rock The work value w f is the independent variable, and the rock uniaxial compressive strength value R c is the dependent variable, which is obtained by linear regression.
可选地，本发明提供的岩层参数的确定装置还包括数据采集器，所述数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；所述数据采集器包括钻速传感器、转速传感器、压力传感器和扭矩传感器；所述钻速传感器用于采集钻进速度V；所述转速传感器用于采集钻杆的转速n；所述压力传感器用于采集钻杆的推力F；所述扭矩传感器用于采集钻杆的扭矩M；还包括显示器，所述显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块获取的数据和确定模块确定的数据；所述获取模块用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值R c；所述确定模块对获取模块获取的数据信息进行处理。 Optionally, the rock formation parameter determination device provided by the present invention further includes a data collector for collecting test data when the drill pipe is drilled into a unit volume of rock; the data collector includes a drilling speed sensor and a rotation speed sensor , Pressure sensor and torque sensor; the drilling speed sensor is used to collect the drilling speed V; the rotation speed sensor is used to collect the rotation speed n of the drill rod; the pressure sensor is used to collect the thrust F of the drill rod; the torque sensor It is used to collect the torque M of the drill pipe; it also includes a display for displaying the test data collected by the data collector when the drill pipe is drilled into a unit volume of rock, the data obtained by the acquisition module and the data determined by the determination module; The acquisition module is used to receive the test data and the corresponding uniaxial compressive strength value R c when the drill pipe is drilled into a unit volume of rock collected by the data collector; the determination module processes the data information acquired by the acquisition module.
本发明提供的一种岩层参数的确定方法及装置，通过将获取的钻进岩层过程中在不同的数据采集时间采集到的多个试验数据组中的钻进单位体积岩石时的钻杆轴力做功值、钻进单位体积岩石时的钻杆扭矩做功值、钻进单位体积岩石时的钻头与孔底摩擦做功值代入岩石单轴抗压强度值R c的估算公式，来得到多个岩石单轴抗压强度值，并通过K-means聚类分析软件来确定不同岩层分类中的岩石单轴抗压强度值组成的 子序列，通过求子序列中岩石单轴抗压强度值的均值，即得到该岩层的单轴抗压强度均值。从而极大的简化了地质勘探程，且对于井巷工程实现了岩性连续探测，节约了人力和时间成本。 The method and device for determining rock formation parameters provided by the present invention are based on the drill rod axial force when drilling a unit volume of rock in a plurality of test data sets collected at different data collection times during drilling into the rock formation. The work value, the drill pipe torque work value when drilling a unit volume of rock, and the drill bit and the hole bottom friction work value when drilling a unit volume of rock are substituted into the estimation formula of the rock uniaxial compressive strength value R c to obtain multiple rock units. Axial compressive strength value, and K-means cluster analysis software is used to determine the sub-sequence of rock uniaxial compressive strength values in different rock formation classifications. By finding the mean value of the rock uniaxial compressive strength values in the sub-sequence, that is Obtain the average uniaxial compressive strength of the rock formation. This greatly simplifies the geological exploration process, and realizes continuous lithology detection for the tunnel engineering, saving labor and time costs.
附图说明Description of the drawings
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。The drawings herein are incorporated into the specification and constitute a part of the specification, show embodiments in accordance with the present invention, and together with the specification are used to explain the principle of the present invention.
图1为本实施例的岩层参数的确定方法的流程示意图；Figure 1 is a schematic flow chart of the method for determining rock formation parameters of this embodiment;
图2为本实施例的岩层参数的确定装置的连接示意图。Figure 2 is a schematic diagram of the connection of the rock formation parameter determination device of this embodiment.
通过上述附图，已示出本发明明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本发明构思的范围，而是通过参考特定实施例为本领域技术人员说明本发明的概念。Through the above-mentioned drawings, the specific embodiments of the present invention have been shown, which will be described in more detail below. These drawings and text descriptions are not intended to limit the scope of the inventive concept in any way, but to explain the concept of the invention to those skilled in the art by referring to specific embodiments.
具体实施方式detailed description
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。In order to make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be described clearly and completely in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments It is a part of the embodiments of the present invention, not all the embodiments.
基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本实施例保护的范围。在不冲突的情况下，下述的实施例及实施例中的特征可以相互组合。Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of this embodiment. In the case of no conflict, the following embodiments and features in the embodiments can be combined with each other.
如图1所示，本实施例提供的岩层参数的确定方法，包括：As shown in Figure 1, the method for determining rock formation parameters provided in this embodiment includes:
S1、获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f； S1. Obtain g test data sets collected at different data collection times t during the process of drilling into the rock formation. The test data in the test data set include the value of drill pipe axial force w F when drilling a unit volume of rock, and drilling The drill pipe torque work value w M per unit volume of rock, the friction work value w f between the drill bit and the bottom of the hole when drilling a unit volume of rock;
示例性的，在钻机对岩层进行钻探的过程中，通过实验装置可以直接或者间接的获取单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f这类的实验数据；当实验装置在获取这些数据时，可以同时记录下获取这些数据时的时间，这些数据获取的时间可以是同步的。 Exemplarily, in the process of drilling the rock formation by the rig, the experimental device can directly or indirectly obtain the drill pipe axial force work value w F per unit volume of rock, and the drill pipe torque work value w when drilling unit volume of rock. M. Experimental data such as the friction work value w f between the drill bit and the bottom of the hole when drilling a unit volume of rock; when the experimental device is acquiring these data, the time when the data is acquired can be recorded at the same time, and the time for the data acquisition can be It is synchronized.
其中，g个试验数据组表示多个试验数据组，并不是对数据组的个数作具体的限定，只是为了后续的表述方便，对获取的试验数据组的个数所做的表达。Among them, g test data sets represent multiple test data sets, which are not a specific limitation on the number of data sets, but are merely an expression of the number of acquired test data sets for the convenience of subsequent expressions.
S2、将g个试验数据组中的试验数据分别代入岩石单轴抗压强度值R c的估算公式，以确定在不同的数据采集时间t采集到的g个试验数据组分别对应的g个岩石单轴抗压强度值R ci，其中R ci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值； S2. Substitute the test data in the g test data groups into the estimation formula of the rock uniaxial compressive strength value R c to determine the g rocks corresponding to the g test data groups collected at different data collection times t Uniaxial compressive strength value R ci , where R ci is the uniaxial compressive strength value of the rock corresponding to the i-th test data group obtained in the order of collection time;
示例性的，岩石单轴抗压强度值R c的估算公式是预先已经确定的，其中岩石单轴抗压强度值R c的估算公式可以是直接用其他文献中已经确定的参考公式，也可以通过对另外的实验数据进行线性回归来得到。 Exemplarily, the formula for estimating the uniaxial compressive strength value of rock R c has been determined in advance, and the formula for estimating the uniaxial compressive strength value of rock R c can be directly used by reference formulas already determined in other documents, or It is obtained by linear regression on other experimental data.
S3、将确定的g个岩石单轴抗压强度值R ci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(R c1，R c2，R c3，……，R ci，……，R cg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；最终分类结果包括最佳分类数d s和每一分类的子序列；其中最佳分类数d s即为岩层 按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为A e＝(R ca，R ca+1，……，R cb)时，则第e层岩层的单轴抗压强度均值R为： S3. The determined uniaxial compressive strength values R ci of g rocks are formed according to the time sequence of data collection to form a rock uniaxial compressive strength data sequence A=(R c1 , R c2 , R c3 ,..., R ci ,..., R cg ), use K-means clustering analysis software to determine the final classification result of the data sequence A; the final classification result includes the best classification number d s and the subsequence of each classification; the best classification The number d s is the total number of layers divided according to the uniaxial compressive strength of the rock; when the e-th subsequence is A e = (R ca , R ca+1 ,..., R cb ), then the first The mean value R of uniaxial compressive strength of e-layer rock is:
Figure PCTCN2019104107-appb-000011
Figure PCTCN2019104107-appb-000011
示例性的，由于在钻探过程中，随着钻探深度的增加，岩层的岩性发生变化，在这个过程中所获取的g个岩石单轴抗压强度值R ci自身具有时序性，为了得到连续的岩层的岩性，当对这些数据进行分析的时候，也有时序性的条件限制。子序列为数据序列A中的某一段序列，且多个子序列之间不出现重叠的情况。在最终分类结果中，数据序列A分为d s个子序列，且将多各子序列按照时间顺序排列后，其数据的连续性和数据序列A可以完全相同。 Exemplarily, during the drilling process, as the drilling depth increases, the lithology of the rock formation changes, the g rock uniaxial compressive strength values R ci obtained in this process have time series themselves, in order to obtain continuous The lithology of the rock formations, when analyzing these data, also has time sequence constraints. The subsequence is a certain segment of the data sequence A, and there is no overlap between multiple subsequences. In the final classification result, the data sequence A is divided into d s subsequences, and after the multiple subsequences are arranged in chronological order, the continuity of the data and the data sequence A can be exactly the same.
K-means算法是硬聚类算法，是典型的基于原型的目标函数聚类方法的代表，它是数据点到原型的某种距离作为优化的目标函数，利用函数求极值的方法得到迭代运算的调整规则。算法采用误差平方和准则函数作为聚类准则函数。具有K-means算法功能的软件常用的有R软件，Matlab。聚类分析的具体过程是依靠现有的软件完成的，本实施例不对聚类分析的过程和原理进行详细的描述，仅对所需输入的原始数据和所要得到的数据类型作示例性说明。K-means algorithm is a hard clustering algorithm. It is a representative of a typical prototype-based objective function clustering method. It uses a certain distance from a data point to the prototype as the optimized objective function, and iterative operation is obtained by the method of seeking extreme values of the function. Adjustment rules. The algorithm uses the error square sum criterion function as the clustering criterion function. Software with K-means algorithm function commonly used R software, Matlab. The specific process of cluster analysis is completed by the existing software. This embodiment does not describe the process and principle of cluster analysis in detail, and only exemplifies the raw data that needs to be input and the type of data to be obtained.
可选地，在上述实施例的基础上，本实施例提供的岩层参数的确定方法还包括以下步骤：Optionally, on the basis of the foregoing embodiment, the method for determining rock formation parameters provided in this embodiment further includes the following steps:
S4、获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；确定第e个子序列A e＝(R ca，R ca+1，……，R cb)中岩石单轴抗压强度值所对应的数据采集时间t的区间T e＝[t a，t b]；则第e层岩层的厚度为： S4. Obtain the expression V(t) of the drilling speed V of the drill pipe that varies with the data collection time t. The drilling speed of the drill pipe is the depth of the drill pipe entering the rock formation per second, in m/s; The method for obtaining the expression V(t) of the drilling speed V includes: fitting the drilling speed V of different drill pipes and the acquisition time t of the drilling speed V to the expression V(t) by data fitting; A sub-sequence A e = (R ca , R ca+1 ,..., R cb ) in the interval of data collection time t corresponding to the uniaxial compressive strength of the rock T e = [t a , t b ]; The thickness of layer e is:
Figure PCTCN2019104107-appb-000012
Figure PCTCN2019104107-appb-000012
示例性的，受岩石单轴抗压强度的影响，钻杆在岩层中的钻速会发生变化，钻速传感器可以每隔一定时间，比如几毫秒，记录一次钻杆的钻速值，根据实验中采集得到的某个时间段中不同时间对应的多个钻速值，通过非线性拟合，可以得到随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)。Exemplarily, affected by the uniaxial compressive strength of the rock, the drilling speed of the drill pipe in the rock formation will change. The drilling speed sensor can record the drilling speed value of the drill pipe every certain time, such as a few milliseconds. According to the experiment Through non-linear fitting, the multiple drilling speed values corresponding to different times in a certain time period collected in the data collection time can be obtained by the expression V(t) of the drilling speed V of the drill pipe changing with the data collection time t.
也可以在得到数据采集时间t的区间T e＝[t a，t b]的前提下，结合在此时间段钻探的总深度，将钻速的平均值作为钻杆的钻进速度V。 It is also possible to obtain the average value of the drilling speed as the drilling speed V of the drill pipe under the premise that the interval T e =[t a , t b ] of the data collection time t is obtained, combined with the total drilling depth during this time period.
可选地，在上述实施例的基础上，步骤S1还可以包括：根据公式(3)计算钻进单位体积岩石时的钻杆轴力做功值w F： Optionally, on the basis of the foregoing embodiment, step S1 may further include: calculating the work value w F of the drill rod axial force when drilling a unit volume of rock according to formula (3):
Figure PCTCN2019104107-appb-000013
Figure PCTCN2019104107-appb-000013
根据公式(4)计算钻进单位体积岩石时的钻杆扭矩做功值w M： According to formula (4), calculate the work value w M of drill pipe torque when drilling a unit volume of rock:
Figure PCTCN2019104107-appb-000014
Figure PCTCN2019104107-appb-000014
根据公式(5)计算钻进单位体积岩石时的钻头与孔底摩擦做功值w f： Calculate the friction work value w f between the drill bit and the bottom of the hole when drilling a unit volume of rock according to formula (5):
Figure PCTCN2019104107-appb-000015
Figure PCTCN2019104107-appb-000015
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。In the formula, F is the thrust of the drill pipe, in N; λ is the hole expansion coefficient, that is, the ratio of the area of the cross section of the drill hole to the area of the cross section of the drill pipe; r is the drill bit diameter of the drill pipe, in m ; M is the torque of the drill pipe, in N·m; n is the rotation speed of the drill pipe, in r/min; V is the drilling speed of the drill pipe, which is the depth of the drill pipe entering the rock formation per second, in m /s; μ is the friction coefficient between the drill bit and the bottom of the hole, and the value of μ is 0.21.
示例性的，有些测试装置可以直接得到钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f；而在实际中，也可以通过间接测量的方式来得到。 Exemplarily, some test devices can directly obtain the work value w F of the drill pipe axial force when drilling a unit volume of rock, the work value w M of the drill pipe torque when drilling a unit volume of rock, and the bit and the value of the drill bit when drilling a unit volume of rock. The friction work value w f at the bottom of the hole; in practice, it can also be obtained by indirect measurement.
可选地，在上述实施例的基础上，步骤S2还可以包括：本实施例提供的岩层参数的确定方法中岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。 Optionally, on the basis of the foregoing embodiment, step S2 may further include: in the method for determining rock formation parameters provided in this embodiment, the method for determining the estimation formula of the rock uniaxial compressive strength value R c includes: The drill pipe axial force work value w F in volume rock, the drill pipe torque work value w M when drilling a unit volume of rock, the drill bit and the hole bottom friction work value w f when drilling a unit volume of rock are the independent variables. The uniaxial compressive strength value R c is the dependent variable and is obtained by linear regression.
示例性的，线性回归方法时科学研究中常用的一种数据处理方法，即通过预先获取可靠的实验数据来得到经验公式，然后利用该经验公式对结果进行估算的一种方法。多元线性回归的软件可以使用SPSS(Statistical Product and Service Solutions，统计产品与服务解决方案)、Microsoft Excel等可以进行数据处理的软件。Exemplarily, the linear regression method is a data processing method commonly used in scientific research, that is, a method of obtaining an empirical formula by obtaining reliable experimental data in advance, and then using the empirical formula to estimate the result. The software for multiple linear regression can use SPSS (Statistical Product and Service Solutions), Microsoft Excel and other software that can perform data processing.
此外，如图2所示，本实施例还提供了岩层参数的确定装置，包括获取模块，获取模块用于获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f；还包括确定模块，确定模块用于将g个试验数据组中的试验数据分别代入岩石单轴抗压强度值R c的估算公式，以确定在不同的数据采集时间t采集到的g个试验数据组分别对应的g个岩石单轴抗压强度值R ci，其中R ci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；确定模块还用于将确定的g个岩石单轴抗压强度值R ci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(R c1，R c2，R c3，……，R ci，……，R cg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；最终分类结果包括最佳分类数d s和每一分类的子序列；其中最佳分类数d s即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为A e＝(R ca，R ca+1，……，R cb)时，则第e层岩层的单轴抗压强度均值R为： In addition, as shown in FIG. 2, this embodiment also provides a rock formation parameter determination device, including an acquisition module, which is used to acquire g test data sets collected at different data acquisition times t during drilling into the rock formation. The test data in the test data group includes the work value w F of the drill pipe axial force when drilling a unit volume of rock, the work value w M of the drill pipe torque when drilling a unit volume of rock, the drill bit and the bottom of the hole when drilling a unit volume of rock Friction work value w f ; It also includes a determination module, which is used to substitute the test data in the g test data groups into the estimation formula of the rock uniaxial compressive strength value R c to determine the collection at different data collection times t The obtained g test data sets correspond to g rock uniaxial compressive strength values R ci , where R ci is the rock uniaxial compressive strength value corresponding to the i-th test data set obtained in the order of collection time; determine; The module is also used to convert the determined uniaxial compressive strength values R ci of g rocks into a data sequence A=(R c1 , R c2 , R c3 ,... , R ci ,..., R cg ), use the K-means clustering analysis software to determine the final classification result of the data sequence A; the final classification result includes the best classification number d s and the subsequence of each classification; the most The best classification number d s is the total number of layers divided according to the uniaxial compressive strength of the rock; when the e-th subsequence is A e = (R ca , R ca+1 ,..., R cb ), Then the average value of uniaxial compressive strength R of the e-th layer is:
Figure PCTCN2019104107-appb-000016
Figure PCTCN2019104107-appb-000016
可选地，本实施例提供的岩层参数的确定装置中获取模块还用于获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得 到表达式V(t)；确定第e个子序列A e＝(R ca，R ca+1，……，R cb)中岩石单轴抗压强度值所对应的数据采集时间t的区间T e＝[t a，t b]；则第e层岩层的厚度为： Optionally, the acquisition module in the rock formation parameter determination device provided in this embodiment is also used to acquire the expression V(t) of the drilling speed V of the drill pipe that varies with the data collection time t, and the drilling speed of the drill pipe is The depth at which the drill pipe enters the rock formation per second, in m/s; the method for obtaining the expression V(t) of the drilling speed V of the drill pipe includes: the difference between the drilling speed V and the drilling speed V of different drill pipes The acquisition time t obtains the expression V(t) through data fitting; determine the data corresponding to the uniaxial compressive strength of the rock in the e-th subsequence A e = (R ca , R ca+1 ,..., R cb ) The interval of the collection time t T e =[t a , t b ]; then the thickness of the e-th layer is:
Figure PCTCN2019104107-appb-000017
Figure PCTCN2019104107-appb-000017
可选地，本实施例提供的岩层参数的确定装置中，获取模块还用于根据公式(3)获取钻进单位体积岩石时的钻杆轴力做功值w F： Optionally, in the rock formation parameter determination device provided in this embodiment, the acquisition module is further used to acquire the value w F of the axial force of the drill pipe when drilling a unit volume of rock according to formula (3):
Figure PCTCN2019104107-appb-000018
Figure PCTCN2019104107-appb-000018
获取模块还用于根据公式(4)获取钻进单位体积岩石时的钻杆扭矩做功值w M： The obtaining module is also used to obtain the work value w M of drill pipe torque when drilling a unit volume of rock according to formula (4):
Figure PCTCN2019104107-appb-000019
Figure PCTCN2019104107-appb-000019
获取模块还用于根据公式(5)获取钻进单位体积岩石时的钻头与孔底摩擦做功值w f： The acquisition module is also used to acquire the friction work value w f between the drill bit and the hole bottom when drilling a unit volume of rock according to formula (5):
Figure PCTCN2019104107-appb-000020
Figure PCTCN2019104107-appb-000020
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。In the formula, F is the thrust of the drill pipe, in N; λ is the hole expansion coefficient, that is, the ratio of the area of the cross section of the drill hole to the area of the cross section of the drill pipe; r is the drill bit diameter of the drill pipe, in m ; M is the torque of the drill pipe, in N·m; n is the rotation speed of the drill pipe, in r/min; V is the drilling speed of the drill pipe, which is the depth of the drill pipe entering the rock formation per second, in m /s; μ is the friction coefficient between the drill bit and the bottom of the hole, and the value of μ is 0.21.
示例性的，钻杆的推力F可以通过钻机与地面接触的压力减去钻机自身的重力得到；扩孔系数λ为钻杆成孔后钻孔横截面与钻杆横截面的面积的比值；钻杆的扭矩M可以通过钻机上连接的扭矩传感器测量得到；钻杆的转速n可以通过钻机上连接的转速传感器测量得到；钻杆的钻进速度V可以通过钻机上连接的钻速传感器测量得到。Exemplarily, the thrust F of the drill pipe can be obtained by subtracting the gravity of the drill rig from the pressure of the rig in contact with the ground; the reaming coefficient λ is the ratio of the area of the cross section of the drill pipe after the drill pipe is formed; The torque M of the rod can be measured by a torque sensor connected to the drill; the speed n of the drill rod can be measured by a speed sensor connected to the drill; the drilling speed V of the drill rod can be measured by a drilling speed sensor connected to the drill.
可选地，本实施例提供的岩层参数的确定装置中确定模块还用于确定岩石单轴抗压强度值R c的估算公式；岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。 Alternatively, the formation parameter determining means provided in this embodiment in determining module further configured to determine the estimated equation uniaxial compressive strength value R c; a method of determining the uniaxial compressive strength estimation formula of the value R c comprises : Take the drill pipe axial force work value w F when drilling a unit volume of rock, the drill pipe torque work value w M when drilling a unit volume of rock, and the friction work value w f between the drill bit and the bottom of the hole when drilling a unit volume of rock. The independent variable, taking the rock uniaxial compressive strength value R c as the dependent variable, is obtained by linear regression.
可选地，本实施例提供的岩层参数的确定装置还包括数据采集器，数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；数据采集器包括钻速传感器、转速传感器、压力传感器和扭矩传感器；钻速传感器用于采集钻进速度V；转速传感器用于采集钻杆的转速n；压力传感器用于采集钻杆的推力F；扭矩传感器用于采集钻杆的扭矩M；还包括显示器，显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块获取的数据和确定模块确定的数据；获取模块用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值R c；确定模块对获取模块获取的数据信息进行处理。 Optionally, the rock formation parameter determination device provided in this embodiment further includes a data collector, which is used to collect test data when the drill pipe is drilled into a unit volume of rock; the data collector includes a drilling speed sensor, a rotation speed sensor, and a pressure sensor And torque sensor; the drilling speed sensor is used to collect the drilling speed V; the speed sensor is used to collect the rotation speed n of the drill rod; the pressure sensor is used to collect the thrust F of the drill rod; the torque sensor is used to collect the torque M of the drill rod; The display is used to display the test data collected by the data collector when the drill pipe is drilled into a unit volume of rock, the data obtained by the acquisition module and the data determined by the determination module; the acquisition module is used to receive the drill pipe drilling unit collected by the data collector The test data of the volume rock and the corresponding uniaxial compressive strength value R c ; the determination module processes the data information acquired by the acquisition module.
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件 来完成。前述的程序可以存储于一计算机可读取存储介质。该程序在执行时，执行包括上述个方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。A person of ordinary skill in the art can understand that all or part of the steps of the foregoing method embodiments can be implemented by a program instructing relevant hardware. The aforementioned program can be stored in a computer readable storage medium. When the program is executed, the steps including the foregoing method embodiments are executed; and the foregoing storage medium includes: ROM, RAM, magnetic disk, or optical disk and other media that can store program codes.
有学者在实验室进行了相关研究，采用直径为60mm的PDC钻头进行钻进实验。钻进试件分为不同强度的砂浆试件28组，编号J1-J28，砂岩试件8组，编号为S1-S8。如表1所示为其中的30组实验数据。Some scholars have conducted related research in the laboratory, using a PDC drill bit with a diameter of 60mm for drilling experiments. The drilling specimens are divided into 28 groups of mortar specimens with different strengths, numbered J1-J28, and 8 groups of sandstone specimens, numbered S1-S8. Table 1 shows the 30 sets of experimental data.
以表1中数据为基础，通过岩石单轴抗压强度值R c的估算模型预设为公式(1)进行多元线性回归，得到岩石单轴抗压强度值R c的估算公式： Based on the data in Table 1, the estimation model of rock uniaxial compressive strength value R c is preset as formula (1) to perform multiple linear regression to obtain the estimation formula of rock uniaxial compressive strength value R c :
Figure PCTCN2019104107-appb-000021
Figure PCTCN2019104107-appb-000021
表1 钻进实验数据Table 1 Drilling experiment data
Figure PCTCN2019104107-appb-000022
Figure PCTCN2019104107-appb-000022

Claims (9)

  1. 一种岩层参数的确定方法，其特征在于，包括：A method for determining rock formation parameters is characterized in that it includes:
    获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，所述试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f； Obtain g test data sets collected at different data acquisition times t during the process of drilling into the rock formation. The test data in the test data sets include the drill pipe axial work value w F when drilling a unit volume of rock, and drilling The drill pipe torque work value w M per unit volume of rock, the friction work value w f between the drill bit and the bottom of the hole when drilling a unit volume of rock;
    将所述g个试验数据组中的试验数据分别代入岩石单轴抗压强度值R c的估算公式，以确定在不同的数据采集时间t采集到的所述g个试验数据组分别对应的g个岩石单轴抗压强度值R ci，其中R ci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值； Substituting the test data in the g test data sets into the estimation formula of the rock uniaxial compressive strength value R c to determine the g corresponding to the g test data sets collected at different data collection times t The uniaxial compressive strength value R ci of the rock, where R ci is the uniaxial compressive strength value of the rock corresponding to the i-th test data group obtained in the order of collection time;
    将确定的g个岩石单轴抗压强度值R ci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(R c1，R c2，R c3，……，R ci，……，R cg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；所述最终分类结果包括最佳分类数d s和每一分类的子序列；其中最佳分类数d s即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为A e＝(R ca，R ca+1，……，R cb)时，则第e层岩层的单轴抗压强度均值R为： The determined uniaxial compressive strength values R ci of g rocks are formed according to the time sequence of data collection to form a rock uniaxial compressive strength data sequence A=(R c1 , R c2 , R c3 ,..., R ci , ……, R cg ), the final classification result of the data sequence A is determined by K-means clustering analysis software; the final classification result includes the best classification number d s and the subsequence of each classification; the best classification The number d s is the total number of layers divided according to the uniaxial compressive strength of the rock; when the e-th subsequence is A e = (R ca , R ca+1 ,..., R cb ), then the first The mean value R of uniaxial compressive strength of e-layer rock is:
    Figure PCTCN2019104107-appb-100001
    Figure PCTCN2019104107-appb-100001
  2. 根据权利要求1所述的确定方法，其特征在于，还包括：The determining method according to claim 1, characterized in that it further comprises:
    获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，所述钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；所述钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；Obtain the expression V(t) of the drilling speed V of the drill pipe that changes with the data collection time t. The drilling speed of the drill pipe is the depth of the drill pipe entering the rock formation per second, in m/s; The method for obtaining the expression V(t) of the drilling speed V of the rod includes: fitting the drilling speed V of different drill rods and the acquisition time t of the drilling speed V to the expression V(t) by data fitting;
    确定第e个子序列A e＝(R ca，R ca+1，……，R cb)中岩石单轴抗压强度值所对应的数据采集时间t的区间T e＝[t a，t b]； Determine the interval of data collection time t corresponding to the uniaxial compressive strength of the rock in the e-th subsequence A e = (R ca , R ca+1 ,..., R cb ) T e =[t a , t b ] ；
    则第e层岩层的厚度为：The thickness of the e-th rock layer is:
    Figure PCTCN2019104107-appb-100002
    Figure PCTCN2019104107-appb-100002
  3. 根据权利要求1所述的确定方法，其特征在于：The determining method according to claim 1, wherein:
    根据公式(3)计算所述钻进单位体积岩石时的钻杆轴力做功值w F： Calculate the work value w F of the drill pipe axial force when drilling the unit volume of rock according to formula (3):
    Figure PCTCN2019104107-appb-100003
    Figure PCTCN2019104107-appb-100003
    根据公式(4)计算所述钻进单位体积岩石时的钻杆扭矩做功值w M： According to formula (4), calculate the work value w M of the drill pipe torque when drilling a unit volume of rock:
    Figure PCTCN2019104107-appb-100004
    Figure PCTCN2019104107-appb-100004
    根据公式(5)计算所述钻进单位体积岩石时的钻头与孔底摩擦做功值w f： Calculate the friction work value w f between the drill bit and the bottom of the hole when drilling the unit volume of rock according to formula (5):
    Figure PCTCN2019104107-appb-100005
    Figure PCTCN2019104107-appb-100005
    式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。In the formula, F is the thrust of the drill pipe, in N; λ is the hole expansion coefficient, that is, the ratio of the area of the cross section of the drill hole to the area of the cross section of the drill pipe; r is the drill bit diameter of the drill pipe, in m ; M is the torque of the drill pipe, in N·m; n is the rotation speed of the drill pipe, in r/min; V is the drilling speed of the drill pipe, which is the depth of the drill pipe entering the rock formation per second, in m /s; μ is the friction coefficient between the drill bit and the bottom of the hole, and the value of μ is 0.21.
  4. 根据权利要求1所述的确定方法，其特征在于，所述岩石单轴抗压强度值R c的估算公式的确定方法包括： The determining method according to claim 1, wherein the method for determining the estimation formula of the uniaxial compressive strength value R c of the rock comprises:
    以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。 Take the drill pipe axial force work value w F when drilling a unit volume of rock, the drill pipe torque work value w M when drilling a unit volume of rock, and the friction work value w f between the drill bit and the bottom of the hole when drilling a unit volume of rock. The variable, taking the rock uniaxial compressive strength value R c as the dependent variable, is obtained by linear regression.
  5. 一种岩层参数的确定装置，其特征在于，包括：A device for determining rock formation parameters is characterized in that it comprises:
    获取模块，所述获取模块用于获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，所述试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f； The acquisition module is used to acquire g test data sets collected at different data acquisition times t during the process of drilling the rock formation, and the test data in the test data set includes the drill pipe when drilling a unit volume of rock Axial force work value w F , drill pipe torque work value w M when drilling a unit volume of rock, and friction work value w f between the drill bit and the hole bottom when drilling a unit volume of rock;
    确定模块，所述确定模块用于将所述g个试验数据组中的试验数据分别代入岩石单轴抗压强度值R c的估算公式，以确定在不同的数据采集时间t采集到的所述g个试验数据组分别对应的g个岩石单轴抗压强度值R ci，其中R ci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值； The determination module is used to substitute the test data in the g test data groups into the estimation formula of the uniaxial compressive strength value R c of the rock to determine the data collected at different data collection times t g rock uniaxial compressive strength values R ci corresponding to g test data sets, where R ci is the rock uniaxial compressive strength value corresponding to the i-th test data set obtained in the order of collection time;
    所述确定模块还用于将确定的g个岩石单轴抗压强度值R ci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(R c1，R c2，R c3，……，R ci，……，R cg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；所述最终分类结果包括最佳分类数d s和每一分类的子序列；其中最佳分类数d s即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为A e＝(R ca，R ca+1，……，R cb)时，则第e层岩层的单轴抗压强度均值R为： The determining module is also used to form a data sequence of the uniaxial compressive strength values of the rock A=(R c1 , R c2 , R c3 , according to the time sequence of data collection, R ci of the determined g rock uniaxial compressive strength values ,..., R ci ,..., R cg ), use the K-means clustering analysis software to determine the final classification result of the data sequence A; the final classification result includes the best classification number d s and each classification The sub-sequence; the best classification number d s is the total number of layers divided by the rock uniaxial compressive strength value; when the e-th sub-sequence is A e = (R ca , R ca+1 ,..., When R cb ), the average uniaxial compressive strength R of the e-th layer is:
    Figure PCTCN2019104107-appb-100006
    Figure PCTCN2019104107-appb-100006
  6. 根据权利要求5所述的确定装置，其特征在于：The determining device according to claim 5, characterized in that:
    所述获取模块还用于获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，所述钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；所述钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；The acquisition module is also used to acquire the expression V(t) of the drilling speed V of the drill pipe that varies with the data acquisition time t, the drilling speed of the drill pipe is the depth of the drill pipe entering the rock formation per second, and the unit is m/s; the method for obtaining the expression V(t) of the drilling speed V of the drill pipe includes: the acquisition time t of the drilling speed V and the drilling speed V of different drill pipes is expressed by data fitting Formula V(t);
    确定第e个子序列A e＝(R ca，R ca+1，……，R cb)中岩石单轴抗压强度值所对应的数据采集时间t的区间T e＝[t a，t b]； Determine the interval of data collection time t corresponding to the uniaxial compressive strength of the rock in the e-th subsequence A e = (R ca , R ca+1 ,..., R cb ) T e =[t a , t b ] ；
    则第e层岩层的厚度为：The thickness of the e-th rock layer is:
    Figure PCTCN2019104107-appb-100007
    Figure PCTCN2019104107-appb-100007
  7. 根据权利要求5所述的确定装置，其特征在于：The determining device according to claim 5, characterized in that:
    所述获取模块还用于根据公式(3)获取所述钻进单位体积岩石时的钻杆轴力做功值w F： The obtaining module is also used to obtain the value w F of the axial force of the drill pipe when drilling a unit volume of rock according to formula (3):
    Figure PCTCN2019104107-appb-100008
    Figure PCTCN2019104107-appb-100008
    所述获取模块还用于根据公式(4)获取所述钻进单位体积岩石时的钻杆扭矩做功值w M： The obtaining module is also used to obtain the drill pipe torque work value w M when drilling a unit volume of rock according to formula (4):
    Figure PCTCN2019104107-appb-100009
    Figure PCTCN2019104107-appb-100009
    所述获取模块还用于根据公式(5)获取所述钻进单位体积岩石时的钻头与孔底摩擦做功值w f： The acquisition module is also used to acquire the friction work value w f between the drill bit and the hole bottom when drilling the unit volume of rock according to formula (5):
    Figure PCTCN2019104107-appb-100010
    Figure PCTCN2019104107-appb-100010
    式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。In the formula, F is the thrust of the drill pipe, in N; λ is the hole expansion coefficient, that is, the ratio of the area of the cross section of the drill hole to the area of the cross section of the drill pipe; r is the drill bit diameter of the drill pipe, in m ; M is the torque of the drill pipe, in N·m; n is the rotation speed of the drill pipe, in r/min; V is the drilling speed of the drill pipe, which is the depth of the drill pipe entering the rock formation per second, in m /s; μ is the friction coefficient between the drill bit and the bottom of the hole, and the value of μ is 0.21.
  8. 根据权利要求5所述的确定装置，其特征在于：The determining device according to claim 5, characterized in that:
    所述确定模块还用于确定所述岩石单轴抗压强度值R c的估算公式； The determining module is also used to determine the estimation formula of the uniaxial compressive strength value R c of the rock;
    所述岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。 The method for determining the estimation formula of the rock uniaxial compressive strength value R c includes: the work value w F of the drill pipe axial force when drilling a unit volume of rock, and the work value w M of the drill pipe torque when drilling a unit volume of rock. , The friction work value w f between the drill bit and the bottom of the hole when drilling a unit volume of rock is the independent variable, and the uniaxial compressive strength value R c of the rock is the dependent variable, which is obtained by linear regression.
  9. 根据权利要求5所述的确定装置，其特征在于，还包括：The determining device according to claim 5, further comprising:
    数据采集器，所述数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；所述数据采集器包括钻速传感器、转速传感器、压力传感器和扭矩传感器；所述钻速传感器用于采集钻进速度V；所述转速传感器用于采集钻杆的转速n；所述压力传感器用于采集钻杆的推力F；所述扭矩传感器用于采集钻杆的扭矩M；A data collector, the data collector is used to collect test data when the drill pipe is drilled into a unit volume of rock; the data collector includes a drilling speed sensor, a rotation speed sensor, a pressure sensor, and a torque sensor; the drilling speed sensor is used to collect Drilling speed V; the rotational speed sensor is used to collect the rotational speed n of the drill rod; the pressure sensor is used to collect the thrust F of the drill rod; the torque sensor is used to collect the torque M of the drill rod;
    显示器，所述显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块获取的数据和确定模块确定的数据；A display, the display is used to display the test data collected by the data collector when the drill pipe is drilled into a unit volume of rock, the data obtained by the obtaining module, and the data determined by the determining module;
    所述获取模块用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值R c；所述确定模块对获取模块获取的数据信息进行处理。 The acquisition module is used to receive the test data and the corresponding uniaxial compressive strength value R c when the drill pipe is drilled into a unit volume of rock collected by the data collector; the determination module processes the data information acquired by the acquisition module.
PCT/CN2019/104107 2019-04-01 2019-09-03 Method and device for determining rock stratum parameters WO2020199495A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
AU2019439997A AU2019439997A1 (en) 2019-04-01 2019-09-03 Method and device for determining rock stratum parameters
RU2020141122A RU2762675C9 (en) 2019-04-01 2019-09-03 Method and apparatus for determining the parameters of a rock formation

Applications Claiming Priority (2)

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

Publications (1)

Publication Number Publication Date
WO2020199495A1 true WO2020199495A1 (en) 2020-10-08

Family

ID=67569123

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/104107 WO2020199495A1 (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)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112412394A (en) * 2020-11-11 2021-02-26 安徽理工大学 Drilling layered filling method
CN113482599A (en) * 2021-07-08 2021-10-08 中国科学院武汉岩土力学研究所 Portable in-situ test and identification data acquisition system and equipment
CN116485225A (en) * 2023-03-15 2023-07-25 西南交通大学 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
CN117288587A (en) * 2023-11-24 2023-12-26 中国矿业大学(北京) Rock tensile strength while drilling test method and system

Families Citing this family (5)

* Cited by examiner, † Cited by third party
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

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1671946A (en) * 2002-06-06 2005-09-21 贝克休斯公司 Method for in-situ analysis of formation parameters
WO2006127892A2 (en) * 2005-05-25 2006-11-30 Geomechanics International, Inc. Methods and devices for analyzing and controlling the propagation of waves in a borehole generated by water hammer
US7236886B2 (en) * 2004-11-04 2007-06-26 Baker Hughes Incorporated Multiscale multidimensional well log data inversion and deep formation imaging method
CN101460698A (en) * 2006-12-15 2009-06-17 哈里伯顿能源服务公司 Antenna coupling component measurement tool having rotating antenna configuration
CN103443657A (en) * 2011-03-11 2013-12-11 兰德马克绘图国际公司 Methods and systems of estimating formation parameters
CN110130883A (en) * 2019-04-01 2019-08-16 中国矿业大学 The determination method and device of formation parameters

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN104217124A (en) * 2014-09-15 2014-12-17 天津大学 TBM (Tunnel Boring Machine) construction surrounding rock classification method depending on engineering sample data
CN108519284B (en) * 2018-03-30 2019-09-27 山东大学 Underground engineering wall rock digital rig real-time graded method
CN108663269B (en) * 2018-03-30 2019-12-24 山东大学 Underground engineering surrounding rock digital drilling zoning method based on equivalent rock mass strength
CN108760487A (en) * 2018-04-16 2018-11-06 西安科技大学 Deep cavern country rock subregion failure evolvement analysis method based on rockbolt stress analysis

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1671946A (en) * 2002-06-06 2005-09-21 贝克休斯公司 Method for in-situ analysis of formation parameters
US7236886B2 (en) * 2004-11-04 2007-06-26 Baker Hughes Incorporated Multiscale multidimensional well log data inversion and deep formation imaging method
WO2006127892A2 (en) * 2005-05-25 2006-11-30 Geomechanics International, Inc. Methods and devices for analyzing and controlling the propagation of waves in a borehole generated by water hammer
CN101460698A (en) * 2006-12-15 2009-06-17 哈里伯顿能源服务公司 Antenna coupling component measurement tool having rotating antenna configuration
CN103443657A (en) * 2011-03-11 2013-12-11 兰德马克绘图国际公司 Methods and systems of estimating formation parameters
CN110130883A (en) * 2019-04-01 2019-08-16 中国矿业大学 The determination method and device of formation parameters

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112412394A (en) * 2020-11-11 2021-02-26 安徽理工大学 Drilling layered filling method
CN113482599A (en) * 2021-07-08 2021-10-08 中国科学院武汉岩土力学研究所 Portable in-situ test and identification data acquisition system and equipment
CN116485225A (en) * 2023-03-15 2023-07-25 西南交通大学 Automatic acquisition method and system for BQ value of surrounding rock in construction stage based on drilling parameters
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
CN117288587A (en) * 2023-11-24 2023-12-26 中国矿业大学(北京) Rock tensile strength while drilling test method and system
CN117288587B (en) * 2023-11-24 2024-02-20 中国矿业大学(北京) Rock tensile strength while drilling test method and system

Also Published As

Publication number Publication date
AU2019439997A1 (en) 2021-01-07
CN110130883A (en) 2019-08-16
RU2762675C1 (en) 2021-12-21

Similar Documents

Publication Publication Date Title
WO2020199495A1 (en) Method and device for determining rock stratum parameters
Li et al. Sensitivity analysis of drilling parameters in rock rotary drilling process based on orthogonal test method
CN102220865B (en) Method for detecting limestone formation pore pressure
US20160076357A1 (en) Methods for selecting and optimizing drilling systems
CN108301825B (en) High-temperature high-pressure reservoir pore pressure analysis method
BR112015009197B1 (en) METHOD AND SYSTEM FOR PERFORMING A DRILLING OPERATION
US9176255B2 (en) Permeability prediction systems and methods using quadratic discriminant analysis
WO2009129083A1 (en) Continuous measurement of heterogeneity of geomaterials
CN110058323A (en) A kind of tight sand formation brittleness index calculation method
US20140172303A1 (en) Methods and systems for analyzing the quality of a wellbore
BRPI1005003A2 (en) Method of Guiding Well Drilling Operation Using Accidental Measurements, System for Guiding Well Drilling Operation Using Rock Stress Measurements, and Computer Readable Storing Instructions for Guiding a Well Drilling Operation Using Accurate Well Measurements
NO344538B1 (en) Procedure and computer-readable medium for prediction of production zone
CN114818451A (en) Mechanical drilling rate prediction method, device, storage medium and equipment
Tavares et al. Combined effect of non-Darcy flow and formation damage on gas-well performance of dual-porosity and dual-permeability reservoirs
CN104153768A (en) Granite reservoir stratum reservoir performance evaluation method
CN115586086A (en) Borehole wall instability analysis method based on big data
CN103306672A (en) Method for predicting abrasiveness of shaly stratum in different drilling directions
US11867048B2 (en) Method and system based on quantified flowback for formation damage removal
CN111625916A (en) Method and system for calculating stability value of well wall
CN110781571A (en) Drill bit work efficiency evaluation system
CN115434696A (en) Lithology identification database construction method and lithology identification method
Chao et al. A prediction of borehole stability while drilling preliminary prospecting wells based on seismic impedance
FANG et al. Intelligent recognition of tunnel stratum based on advanced drilling tests
Bai et al. Intelligent real-time predicting method for rock characterization based on multi-source information integration while drilling
Wang et al. A case study of drilling process monitoring for geomaterial strength assessment along hydraulic rotary drillhole

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19923534

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019439997

Country of ref document: AU

Date of ref document: 20190903

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19923534

Country of ref document: EP

Kind code of ref document: A1