WO2020199495A1 - 岩层参数的确定方法及装置 - Google Patents

岩层参数的确定方法及装置 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
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.)
Ceased
Application number
PCT/CN2019/104107
Other languages
English (en)
French (fr)
Inventor
刘灿灿
郑西贵
艾德春
刘洪洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology CUMT
China University of Mining and Technology Beijing CUMTB
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 China University of Mining and Technology CUMT, China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology CUMT
Priority to RU2020141122A priority Critical patent/RU2762675C9/ru
Priority to AU2019439997A priority patent/AU2019439997A1/en
Publication of WO2020199495A1 publication Critical patent/WO2020199495A1/zh
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

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

一种岩层参数的确定方法及装置被公开。该方法通过将获取的钻进岩层过程中在不同的数据采集时间采集到的多个试验数据组中的钻进单位体积岩石时的钻杆轴力做功值、钻进单位体积岩石时的钻杆扭矩做功值、钻进单位体积岩石时的钻头与孔底摩擦做功值代入岩石单轴抗压强度值R c的估算公式，来得到多个岩石单轴抗压强度值，并通过K-means聚类分析软件来确定不同岩层分类中的岩石单轴抗压强度值组成的子序列，通过求子序列中岩石单轴抗压强度值的均值，即得到该岩层的单轴抗压强度均值，并可进一步确定各岩层的厚度。该方法对于井巷工程实现了岩性连续探测，节约了人力和时间成本。

Description

岩层参数的确定方法及装置 技术领域
本发明涉及一种岩层参数的确定方法及装置，属于地质勘探工程技术领域。
背景技术
在煤矿、隧道等地下工程中，岩层情况对工程的安全有较大影响。目前，主要通过岩石取芯、钻孔窥视等了解岩层信息，但上述方法影响工程进展速度，且费时费力，不能够实时了解工程所处岩层环境。
目前，石油行业的随钻测井技术相对成熟，但局限性较大，在煤矿等地下工程中还没广泛应用。也有学者提出了基于钻孔岩屑的地层探测法，但岩屑收集困难，且对于岩性易溶于水的岩层，不能够收集到该岩层的岩屑，导致对地层结构的判断不连续，且该方法不能对岩石的单轴抗压强度进行预测。
因此，许多学者提出了基于锚杆钻机的随钻参数实时识别岩层的方法，但目前国内还未找到有效方法使该技术得到广泛应用。
发明内容
为了解决现有技术中存在的上述技术问题，本发明提供了一种岩层参数的确定方法及装置。
本发明提供的岩层参数的确定方法，包括获取钻进岩层过程中在不同的数据采集时间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为：
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层岩层的厚度为：
Figure PCTCN2019104107-appb-000002
可选地，本发明提供的岩层参数的确定方法还包括：根据公式(3)计算所述钻进单位体积岩石时的 钻杆轴力做功值w F：
Figure PCTCN2019104107-appb-000003
根据公式(4)计算所述钻进单位体积岩石时的钻杆扭矩做功值w M：
Figure PCTCN2019104107-appb-000004
根据公式(5)计算所述钻进单位体积岩石时的钻头与孔底摩擦做功值w f：
Figure PCTCN2019104107-appb-000005
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
可选地，本发明提供的岩层参数的确定方法中所述岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。
此外，本发明还提供了岩层参数的确定装置，包括获取模块，所述获取模块用于获取钻进岩层过程中在不同的数据采集时间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为：
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层岩层的厚度为：
Figure PCTCN2019104107-appb-000007
可选地，本发明提供的岩层参数的确定装置中，所述获取模块还用于根据公式(3)获取所述钻进单位体积岩石时的钻杆轴力做功值w F：
Figure PCTCN2019104107-appb-000008
所述获取模块还用于根据公式(4)获取所述钻进单位体积岩石时的钻杆扭矩做功值w M：
Figure PCTCN2019104107-appb-000009
所述获取模块还用于根据公式(5)获取所述钻进单位体积岩石时的钻头与孔底摩擦做功值w f：
Figure PCTCN2019104107-appb-000010
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
可选地，本发明提供的岩层参数的确定装置中所述确定模块还用于确定所述岩石单轴抗压强度值R c的估算公式；所述岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。
可选地，本发明提供的岩层参数的确定装置还包括数据采集器，所述数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；所述数据采集器包括钻速传感器、转速传感器、压力传感器和扭矩传感器；所述钻速传感器用于采集钻进速度V；所述转速传感器用于采集钻杆的转速n；所述压力传感器用于采集钻杆的推力F；所述扭矩传感器用于采集钻杆的扭矩M；还包括显示器，所述显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块获取的数据和确定模块确定的数据；所述获取模块用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值R c；所述确定模块对获取模块获取的数据信息进行处理。
本发明提供的一种岩层参数的确定方法及装置，通过将获取的钻进岩层过程中在不同的数据采集时间采集到的多个试验数据组中的钻进单位体积岩石时的钻杆轴力做功值、钻进单位体积岩石时的钻杆扭矩做功值、钻进单位体积岩石时的钻头与孔底摩擦做功值代入岩石单轴抗压强度值R c的估算公式，来得到多个岩石单轴抗压强度值，并通过K-means聚类分析软件来确定不同岩层分类中的岩石单轴抗压强度值组成的 子序列，通过求子序列中岩石单轴抗压强度值的均值，即得到该岩层的单轴抗压强度均值。从而极大的简化了地质勘探程，且对于井巷工程实现了岩性连续探测，节约了人力和时间成本。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。
图1为本实施例的岩层参数的确定方法的流程示意图；
图2为本实施例的岩层参数的确定装置的连接示意图。
通过上述附图，已示出本发明明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本发明构思的范围，而是通过参考特定实施例为本领域技术人员说明本发明的概念。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。
基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本实施例保护的范围。在不冲突的情况下，下述的实施例及实施例中的特征可以相互组合。
如图1所示，本实施例提供的岩层参数的确定方法，包括：
S1、获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f；
示例性的，在钻机对岩层进行钻探的过程中，通过实验装置可以直接或者间接的获取单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f这类的实验数据；当实验装置在获取这些数据时，可以同时记录下获取这些数据时的时间，这些数据获取的时间可以是同步的。
其中，g个试验数据组表示多个试验数据组，并不是对数据组的个数作具体的限定，只是为了后续的表述方便，对获取的试验数据组的个数所做的表达。
S2、将g个试验数据组中的试验数据分别代入岩石单轴抗压强度值R c的估算公式，以确定在不同的数据采集时间t采集到的g个试验数据组分别对应的g个岩石单轴抗压强度值R ci，其中R ci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；
示例性的，岩石单轴抗压强度值R c的估算公式是预先已经确定的，其中岩石单轴抗压强度值R c的估算公式可以是直接用其他文献中已经确定的参考公式，也可以通过对另外的实验数据进行线性回归来得到。
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为：
Figure PCTCN2019104107-appb-000011
示例性的，由于在钻探过程中，随着钻探深度的增加，岩层的岩性发生变化，在这个过程中所获取的g个岩石单轴抗压强度值R ci自身具有时序性，为了得到连续的岩层的岩性，当对这些数据进行分析的时候，也有时序性的条件限制。子序列为数据序列A中的某一段序列，且多个子序列之间不出现重叠的情况。在最终分类结果中，数据序列A分为d s个子序列，且将多各子序列按照时间顺序排列后，其数据的连续性和数据序列A可以完全相同。
K-means算法是硬聚类算法，是典型的基于原型的目标函数聚类方法的代表，它是数据点到原型的某种距离作为优化的目标函数，利用函数求极值的方法得到迭代运算的调整规则。算法采用误差平方和准则函数作为聚类准则函数。具有K-means算法功能的软件常用的有R软件，Matlab。聚类分析的具体过程是依靠现有的软件完成的，本实施例不对聚类分析的过程和原理进行详细的描述，仅对所需输入的原始数据和所要得到的数据类型作示例性说明。
可选地，在上述实施例的基础上，本实施例提供的岩层参数的确定方法还包括以下步骤：
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层岩层的厚度为：
Figure PCTCN2019104107-appb-000012
示例性的，受岩石单轴抗压强度的影响，钻杆在岩层中的钻速会发生变化，钻速传感器可以每隔一定时间，比如几毫秒，记录一次钻杆的钻速值，根据实验中采集得到的某个时间段中不同时间对应的多个钻速值，通过非线性拟合，可以得到随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)。
也可以在得到数据采集时间t的区间T e＝[t a，t b]的前提下，结合在此时间段钻探的总深度，将钻速的平均值作为钻杆的钻进速度V。
可选地，在上述实施例的基础上，步骤S1还可以包括：根据公式(3)计算钻进单位体积岩石时的钻杆轴力做功值w F：
Figure PCTCN2019104107-appb-000013
根据公式(4)计算钻进单位体积岩石时的钻杆扭矩做功值w M：
Figure PCTCN2019104107-appb-000014
根据公式(5)计算钻进单位体积岩石时的钻头与孔底摩擦做功值w f：
Figure PCTCN2019104107-appb-000015
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
示例性的，有些测试装置可以直接得到钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f；而在实际中，也可以通过间接测量的方式来得到。
可选地，在上述实施例的基础上，步骤S2还可以包括：本实施例提供的岩层参数的确定方法中岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。
示例性的，线性回归方法时科学研究中常用的一种数据处理方法，即通过预先获取可靠的实验数据来得到经验公式，然后利用该经验公式对结果进行估算的一种方法。多元线性回归的软件可以使用SPSS(Statistical Product and Service Solutions，统计产品与服务解决方案)、Microsoft Excel等可以进行数据处理的软件。
此外，如图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为：
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层岩层的厚度为：
Figure PCTCN2019104107-appb-000017
可选地，本实施例提供的岩层参数的确定装置中，获取模块还用于根据公式(3)获取钻进单位体积岩石时的钻杆轴力做功值w F：
Figure PCTCN2019104107-appb-000018
获取模块还用于根据公式(4)获取钻进单位体积岩石时的钻杆扭矩做功值w M：
Figure PCTCN2019104107-appb-000019
获取模块还用于根据公式(5)获取钻进单位体积岩石时的钻头与孔底摩擦做功值w f：
Figure PCTCN2019104107-appb-000020
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
示例性的，钻杆的推力F可以通过钻机与地面接触的压力减去钻机自身的重力得到；扩孔系数λ为钻杆成孔后钻孔横截面与钻杆横截面的面积的比值；钻杆的扭矩M可以通过钻机上连接的扭矩传感器测量得到；钻杆的转速n可以通过钻机上连接的转速传感器测量得到；钻杆的钻进速度V可以通过钻机上连接的钻速传感器测量得到。
可选地，本实施例提供的岩层参数的确定装置中确定模块还用于确定岩石单轴抗压强度值R c的估算公式；岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。
可选地，本实施例提供的岩层参数的确定装置还包括数据采集器，数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；数据采集器包括钻速传感器、转速传感器、压力传感器和扭矩传感器；钻速传感器用于采集钻进速度V；转速传感器用于采集钻杆的转速n；压力传感器用于采集钻杆的推力F；扭矩传感器用于采集钻杆的扭矩M；还包括显示器，显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块获取的数据和确定模块确定的数据；获取模块用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值R c；确定模块对获取模块获取的数据信息进行处理。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件 来完成。前述的程序可以存储于一计算机可读取存储介质。该程序在执行时，执行包括上述个方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
有学者在实验室进行了相关研究，采用直径为60mm的PDC钻头进行钻进实验。钻进试件分为不同强度的砂浆试件28组，编号J1-J28，砂岩试件8组，编号为S1-S8。如表1所示为其中的30组实验数据。
以表1中数据为基础，通过岩石单轴抗压强度值R c的估算模型预设为公式(1)进行多元线性回归，得到岩石单轴抗压强度值R c的估算公式：
Figure PCTCN2019104107-appb-000021
表1 钻进实验数据
Figure PCTCN2019104107-appb-000022

Claims (9)

  1. 一种岩层参数的确定方法，其特征在于，包括：
    获取钻进岩层过程中在不同的数据采集时间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为：
    Figure PCTCN2019104107-appb-100001
  2. 根据权利要求1所述的确定方法，其特征在于，还包括：
    获取随数据采集时间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层岩层的厚度为：
    Figure PCTCN2019104107-appb-100002
  3. 根据权利要求1所述的确定方法，其特征在于：
    根据公式(3)计算所述钻进单位体积岩石时的钻杆轴力做功值w F：
    Figure PCTCN2019104107-appb-100003
    根据公式(4)计算所述钻进单位体积岩石时的钻杆扭矩做功值w M：
    Figure PCTCN2019104107-appb-100004
    根据公式(5)计算所述钻进单位体积岩石时的钻头与孔底摩擦做功值w f：
    Figure PCTCN2019104107-appb-100005
    式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
  4. 根据权利要求1所述的确定方法，其特征在于，所述岩石单轴抗压强度值R c的估算公式的确定方法包括：
    以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。
  5. 一种岩层参数的确定装置，其特征在于，包括：
    获取模块，所述获取模块用于获取钻进岩层过程中在不同的数据采集时间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为：
    Figure PCTCN2019104107-appb-100006
  6. 根据权利要求5所述的确定装置，其特征在于：
    所述获取模块还用于获取随数据采集时间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层岩层的厚度为：
    Figure PCTCN2019104107-appb-100007
  7. 根据权利要求5所述的确定装置，其特征在于：
    所述获取模块还用于根据公式(3)获取所述钻进单位体积岩石时的钻杆轴力做功值w F：
    Figure PCTCN2019104107-appb-100008
    所述获取模块还用于根据公式(4)获取所述钻进单位体积岩石时的钻杆扭矩做功值w M：
    Figure PCTCN2019104107-appb-100009
    所述获取模块还用于根据公式(5)获取所述钻进单位体积岩石时的钻头与孔底摩擦做功值w f：
    Figure PCTCN2019104107-appb-100010
    式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
  8. 根据权利要求5所述的确定装置，其特征在于：
    所述确定模块还用于确定所述岩石单轴抗压强度值R c的估算公式；
    所述岩石单轴抗压强度值R c的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值w F、钻进单位体积岩石时的钻杆扭矩做功值w M、钻进单位体积岩石时的钻头与孔底摩擦做功值w f为自变量，以岩石单轴抗压强度值R c为因变量，通过线性回归方法得到。
  9. 根据权利要求5所述的确定装置，其特征在于，还包括：
    数据采集器，所述数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；所述数据采集器包括钻速传感器、转速传感器、压力传感器和扭矩传感器；所述钻速传感器用于采集钻进速度V；所述转速传感器用于采集钻杆的转速n；所述压力传感器用于采集钻杆的推力F；所述扭矩传感器用于采集钻杆的扭矩M；
    显示器，所述显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块获取的数据和确定模块确定的数据；
    所述获取模块用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值R c；所述确定模块对获取模块获取的数据信息进行处理。
PCT/CN2019/104107 2019-04-01 2019-09-03 岩层参数的确定方法及装置 Ceased WO2020199495A1 (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
RU2020141122A RU2762675C9 (ru) 2019-04-01 2019-09-03 Способ и устройство для определения параметров пласта горной породы
AU2019439997A AU2019439997A1 (en) 2019-04-01 2019-09-03 Method and device for determining rock stratum parameters

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201910255972.8A CN110130883A (zh) 2019-04-01 2019-04-01 岩层参数的确定方法及装置
CN201910255972.8 2019-04-01

Publications (1)

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

Family

ID=67569123

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/104107 Ceased WO2020199495A1 (zh) 2019-04-01 2019-09-03 岩层参数的确定方法及装置

Country Status (3)

Country Link
CN (1) CN110130883A (zh)
AU (1) AU2019439997A1 (zh)
WO (1) WO2020199495A1 (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112412394A (zh) * 2020-11-11 2021-02-26 安徽理工大学 一种钻孔分层充填方法
CN113482599A (zh) * 2021-07-08 2021-10-08 中国科学院武汉岩土力学研究所 便携式原位测试与识别数据采集系统及设备
CN115266339A (zh) * 2022-07-18 2022-11-01 山东黄金矿业科技有限公司深井开采实验室分公司 花岗岩单轴抗压强度的计算方法
CN115824853A (zh) * 2022-11-18 2023-03-21 国能神东煤炭集团有限责任公司 缩短地下水库与煤炭开采空间安全距离的方法
CN115982640A (zh) * 2022-12-20 2023-04-18 鞍钢集团矿业有限公司 一种钻孔过程实时岩石种类判断方法
CN116165268A (zh) * 2022-11-23 2023-05-26 中国电建集团成都勘测设计研究院有限公司 岩石强度检测定距装置及岩石强度检测仪、使用方法
CN116485225A (zh) * 2023-03-15 2023-07-25 西南交通大学 基于钻进参数的施工阶段围岩bq值自动获取方法及系统
CN116499827A (zh) * 2023-04-18 2023-07-28 中国石油大学(北京) 多岩性互层人造岩样制备装置及方法
CN116609184A (zh) * 2023-05-22 2023-08-18 中国石油大学(华东) 一种由上返块体预测原位条件下破碎性地层抗压强度的方法
CN117288587A (zh) * 2023-11-24 2023-12-26 中国矿业大学(北京) 岩体抗拉强度随钻测试方法与系统
CN119712228A (zh) * 2024-12-18 2025-03-28 中国矿业大学 一种基于随钻参数响应特征的巷道顶板稳定性评价方法

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110130883A (zh) * 2019-04-01 2019-08-16 中国矿业大学 岩层参数的确定方法及装置
CN110486007B (zh) * 2019-08-29 2023-02-03 武汉长盛煤安科技有限公司 煤矿随钻围岩力学参数原位测试装置及方法
CN111075424B (zh) * 2019-12-25 2022-11-18 中国石油大学(华东) 一种修正随钻测量参数测量结果的方法
CN111706322B (zh) * 2020-07-17 2023-03-21 中国铁建重工集团股份有限公司 一种岩石钻进响应的预测方法和预测系统
CN111855961B (zh) * 2020-07-24 2021-10-26 中南大学 一种岩体钻进质量检测方法、钻机、服务器、存储介质
CN116879527B (zh) * 2023-06-21 2026-02-24 中铁第四勘察设计院集团有限公司 基于随钻参数实时确定岩体基本质量分级的方法及系统
CN117763466B (zh) * 2024-02-22 2024-07-09 中石化经纬有限公司 一种基于聚类算法的地层可钻性评价方法及系统
CN118392662B (zh) * 2024-06-26 2024-11-05 中国科学院武汉岩土力学研究所 基于岩芯回弹和随钻参数的岩体抗压强度测定方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1671946A (zh) * 2002-06-06 2005-09-21 贝克休斯公司 就地分析岩层参数的方法
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 (zh) * 2006-12-15 2009-06-17 哈里伯顿能源服务公司 具有旋转天线结构的天线耦合元件测量工具
CN103443657A (zh) * 2011-03-11 2013-12-11 兰德马克绘图国际公司 估计岩层参数的方法和系统
CN110130883A (zh) * 2019-04-01 2019-08-16 中国矿业大学 岩层参数的确定方法及装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007327880A (ja) * 2006-06-08 2007-12-20 Tokyo Electric Power Co Inc:The コンクリート圧縮強度測定装置及び測定方法
CN104217124A (zh) * 2014-09-15 2014-12-17 天津大学 以工程样本数据为依托的tbm施工围岩分级方法
CN108519284B (zh) * 2018-03-30 2019-09-27 山东大学 地下工程围岩数字钻探实时分级方法
CN108663269B (zh) * 2018-03-30 2019-12-24 山东大学 基于等效岩体强度的地下工程围岩数字钻探分区方法
CN108760487A (zh) * 2018-04-16 2018-11-06 西安科技大学 基于锚杆受力分析的深埋洞室围岩分区破裂演化分析方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1671946A (zh) * 2002-06-06 2005-09-21 贝克休斯公司 就地分析岩层参数的方法
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 (zh) * 2006-12-15 2009-06-17 哈里伯顿能源服务公司 具有旋转天线结构的天线耦合元件测量工具
CN103443657A (zh) * 2011-03-11 2013-12-11 兰德马克绘图国际公司 估计岩层参数的方法和系统
CN110130883A (zh) * 2019-04-01 2019-08-16 中国矿业大学 岩层参数的确定方法及装置

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112412394A (zh) * 2020-11-11 2021-02-26 安徽理工大学 一种钻孔分层充填方法
CN113482599A (zh) * 2021-07-08 2021-10-08 中国科学院武汉岩土力学研究所 便携式原位测试与识别数据采集系统及设备
CN115266339A (zh) * 2022-07-18 2022-11-01 山东黄金矿业科技有限公司深井开采实验室分公司 花岗岩单轴抗压强度的计算方法
CN115824853A (zh) * 2022-11-18 2023-03-21 国能神东煤炭集团有限责任公司 缩短地下水库与煤炭开采空间安全距离的方法
CN116165268A (zh) * 2022-11-23 2023-05-26 中国电建集团成都勘测设计研究院有限公司 岩石强度检测定距装置及岩石强度检测仪、使用方法
CN115982640A (zh) * 2022-12-20 2023-04-18 鞍钢集团矿业有限公司 一种钻孔过程实时岩石种类判断方法
CN116485225A (zh) * 2023-03-15 2023-07-25 西南交通大学 基于钻进参数的施工阶段围岩bq值自动获取方法及系统
CN116485225B (zh) * 2023-03-15 2023-11-10 西南交通大学 基于钻进参数的施工阶段围岩bq值自动获取方法及系统
CN116499827A (zh) * 2023-04-18 2023-07-28 中国石油大学(北京) 多岩性互层人造岩样制备装置及方法
CN116609184A (zh) * 2023-05-22 2023-08-18 中国石油大学(华东) 一种由上返块体预测原位条件下破碎性地层抗压强度的方法
CN117288587A (zh) * 2023-11-24 2023-12-26 中国矿业大学(北京) 岩体抗拉强度随钻测试方法与系统
CN117288587B (zh) * 2023-11-24 2024-02-20 中国矿业大学(北京) 岩体抗拉强度随钻测试方法与系统
CN119712228A (zh) * 2024-12-18 2025-03-28 中国矿业大学 一种基于随钻参数响应特征的巷道顶板稳定性评价方法

Also Published As

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

Similar Documents

Publication Publication Date Title
CN110130883A (zh) 岩层参数的确定方法及装置
CN110058323A (zh) 一种致密砂岩地层脆性指数计算方法
CN102220865B (zh) 一种检测灰岩地层孔隙压力的方法
US9176255B2 (en) Permeability prediction systems and methods using quadratic discriminant analysis
CN103590828B (zh) 一种录井dck指数法评价地层压力的方法
CN115032062B (zh) 页岩油储层的可压性确定方法、装置和可压性确定系统
Bai et al. Intelligent real-time predicting method for rock characterization based on multi-source information integration while drilling
WO2014099832A1 (en) Methods and systems for analyzing the quality of a wellbore
US20250075616A1 (en) Noise characterization in formation testing
CN116956632A (zh) 一种基于改进Eaton法的地层孔隙压力计算方法
NO344538B1 (no) Fremgangsmåte og datamskinlesbart medium for prediksjon av produksjonssone
CN111625916B (zh) 井壁稳定性值计算方法及系统
CN118916962A (zh) 一种基于多源勘察信息的深部岩体力学参数反演方法
CN107091084A (zh) 建立虚拟数据井提取页岩气水平井储层物性参数的方法
CN104459790B (zh) 含油气性盆地有效储层的分析方法和装置
CN110781571A (zh) 一种钻头工作效率评估系统
CN116624144A (zh) 一种利用工程参数计算泥岩地层纵波速度的方法
CN116071545A (zh) 一种巷道顶板稳定性随钻感知方法
WO2024096889A1 (en) Automatic landing of formation testing tools
CN115434696A (zh) 一种岩性识别数据库构建方法及岩性识别方法
CN116104470A (zh) 地质导向地层识别预测方法、系统、设备和存储介质
Chao et al. A prediction of borehole stability while drilling preliminary prospecting wells based on seismic impedance
CN115408361A (zh) 岩层危险性评估数据库构建方法及岩层的危险性评估方法
CN119986795A (zh) 一种基于孔压静力触探试验的剪切波速估算方法
CN104047593A (zh) 一种利用速度定位钻井微芯片示踪器深度的方法

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