CN110130883A - 岩层参数的确定方法及装置 - Google Patents
岩层参数的确定方法及装置 Download PDFInfo
- Publication number
- CN110130883A CN110130883A CN201910255972.8A CN201910255972A CN110130883A CN 110130883 A CN110130883 A CN 110130883A CN 201910255972 A CN201910255972 A CN 201910255972A CN 110130883 A CN110130883 A CN 110130883A
- Authority
- CN
- China
- Prior art keywords
- rock
- drilling
- drilling rod
- value
- data
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 28
- 238000005553 drilling Methods 0.000 claims abstract description 220
- 239000011435 rock Substances 0.000 claims abstract description 199
- 238000012360 testing method Methods 0.000 claims abstract description 59
- 239000000523 sample Substances 0.000 claims abstract description 9
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 230000035515 penetration Effects 0.000 claims description 40
- 230000006835 compression Effects 0.000 claims description 12
- 238000007906 compression Methods 0.000 claims description 12
- 238000012417 linear regression Methods 0.000 claims description 9
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 241001074085 Scophthalmus aquosus Species 0.000 claims 2
- 238000006467 substitution reaction Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000003245 coal Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000004570 mortar (masonry) Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000001755 vocal effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0218—Calculations based on experimental data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/025—Geometry of the test
- G01N2203/0252—Monoaxial, i.e. the forces being applied along a single axis of the specimen
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Earth Drilling (AREA)
Abstract
本发明提供的一种岩层参数的确定方法及装置，属于地质勘探工程技术领域。通过将获取的钻进岩层过程中在不同的数据采集时间采集到的多个试验数据组中的钻进单位体积岩石时的钻杆轴力做功值、钻进单位体积岩石时的钻杆扭矩做功值、钻进单位体积岩石时的钻头与孔底摩擦做功值代入岩石单轴抗压强度值Rc的估算公式，来得到多个岩石单轴抗压强度值，并通过K‑means聚类分析软件来确定不同岩层分类中的岩石单轴抗压强度值组成的子序列，通过求子序列中岩石单轴抗压强度值的均值，即得到该岩层的单轴抗压强度均值，并可进一步确定各岩层的厚度。从而极大的简化了地质勘探程，且对于井巷工程实现了岩性连续探测，节约了人力和时间成本。
Description
技术领域
本发明涉及一种岩层参数的确定方法及装置，属于地质勘探工程技术领域。
背景技术
在煤矿、隧道等地下工程中，岩层情况对工程的安全有较大影响。目前，主要通过岩石取芯、钻孔窥视等了解岩层信息，但上述方法影响工程进展速度，且费时费力，不能够实时了解工程所处岩层环境。
目前，石油行业的随钻测井技术相对成熟，但局限性较大，在煤矿等地下工程中还没广泛应用。也有学者提出了基于钻孔岩屑的地层探测法，但岩屑收集困难，且对于岩性易溶于水的岩层，不能够收集到该岩层的岩屑，导致对地层结构的判断不连续，且该方法不能对岩石的单轴抗压强度进行预测。
因此，许多学者提出了基于锚杆钻机的随钻参数实时识别岩层的方法，但目前国内还未找到有效方法使该技术得到广泛应用。
发明内容
为了解决现有技术中存在的上述技术问题，本发明提供了一种岩层参数的确定方法及装置。
本发明提供的岩层参数的确定方法，包括获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，所述试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf；将所述g个试验数据组中的试验数据分别代入岩石单轴抗压强度值Rc的估算公式，以确定在不同的数据采集时间t采集到的所述g个试验数据组分别对应的g个岩石单轴抗压强度值Rci，其中Rci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；将确定的g个岩石单轴抗压强度值Rci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(Rc1，Rc2，Rc3，……，Rci，……，Rcg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；所述最终分类结果包括最佳分类数ds和每一分类的子序列；其中最佳分类数ds即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为Ae＝(Rca，Rca+1，……，Rcb)时，则第e层岩层的单轴抗压强度均值R为：
可选地，本发明提供的岩层参数的确定方法还包括获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，所述钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；所述钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；确定第e个子序列Ae＝(Rca，Rca+1，……，Rcb)中岩石单轴抗压强度值所对应的数据采集时间t的区间Te＝[ta，tb]；则第e层岩层的厚度为：
可选地，本发明提供的岩层参数的确定方法还包括：根据公式(3)计算所述钻进单位体积岩石时的钻杆轴力做功值wF：
根据公式(4)计算所述钻进单位体积岩石时的钻杆扭矩做功值wM：
根据公式(5)计算所述钻进单位体积岩石时的钻头与孔底摩擦做功值wf：
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
可选地，本发明提供的岩层参数的确定方法中所述岩石单轴抗压强度值Rc的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf为自变量，以岩石单轴抗压强度值Rc为因变量，通过线性回归方法得到。
此外，本发明还提供了岩层参数的确定装置，包括获取模块，所述获取模块用于获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，所述试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf；还包括确定模块，所述确定模块用于将所述g个试验数据组中的试验数据分别代入岩石单轴抗压强度值Rc的估算公式，以确定在不同的数据采集时间t采集到的所述g个试验数据组分别对应的g个岩石单轴抗压强度值Rci，其中Rci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；所述确定模块还用于将确定的g个岩石单轴抗压强度值Rci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(Rc1，Rc2，Rc3，……，Rci，……，Rcg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；所述最终分类结果包括最佳分类数ds和每一分类的子序列；其中最佳分类数ds即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为Ae＝(Rca，Rca+1，……，Rcb)时，则第e层岩层的单轴抗压强度均值R为：
可选地，本发明提供的岩层参数的确定装置中所述获取模块还用于获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，所述钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；所述钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；确定第e个子序列Ae＝(Rca，Rca+1，……，Rcb)中岩石单轴抗压强度值所对应的数据采集时间t的区间Te＝[ta，tb]；则第e层岩层的厚度为：
可选地，本发明提供的岩层参数的确定装置中，所述获取模块还用于根据公式(3)获取所述钻进单位体积岩石时的钻杆轴力做功值wF：
所述获取模块还用于根据公式(4)获取所述钻进单位体积岩石时的钻杆扭矩做功值wM：
所述获取模块还用于根据公式(5)获取所述钻进单位体积岩石时的钻头与孔底摩擦做功值wf：
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
可选地，本发明提供的岩层参数的确定装置中所述确定模块还用于确定所述岩石单轴抗压强度值Rc的估算公式；所述岩石单轴抗压强度值Rc的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf为自变量，以岩石单轴抗压强度值Rc为因变量，通过线性回归方法得到。
可选地，本发明提供的岩层参数的确定装置还包括数据采集器，所述数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；所述数据采集器包括钻速传感器、转速传感器、压力传感器和扭矩传感器；所述钻速传感器用于采集钻进速度V；所述转速传感器用于采集钻杆的转速n；所述压力传感器用于采集钻杆的推力F；所述扭矩传感器用于采集钻杆的扭矩M；还包括显示器，所述显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块获取的数据和确定模块确定的数据；所述获取模块用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值Rc；所述确定模块对获取模块获取的数据信息进行处理。
本发明提供的一种岩层参数的确定方法及装置，通过将获取的钻进岩层过程中在不同的数据采集时间采集到的多个试验数据组中的钻进单位体积岩石时的钻杆轴力做功值、钻进单位体积岩石时的钻杆扭矩做功值、钻进单位体积岩石时的钻头与孔底摩擦做功值代入岩石单轴抗压强度值Rc的估算公式，来得到多个岩石单轴抗压强度值，并通过K-means聚类分析软件来确定不同岩层分类中的岩石单轴抗压强度值组成的子序列，通过求子序列中岩石单轴抗压强度值的均值，即得到该岩层的单轴抗压强度均值。从而极大的简化了地质勘探程，且对于井巷工程实现了岩性连续探测，节约了人力和时间成本。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。
图1为本实施例的岩层参数的确定方法的流程示意图；
图2为本实施例的岩层参数的确定装置的连接示意图。
附图标记说明：1-获取模块；2-确定模块；3-显示器；
11-钻速传感器；12-转速传感器；13-压力传感器；14-扭矩传感器。
通过上述附图，已示出本发明明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本发明构思的范围，而是通过参考特定实施例为本领域技术人员说明本发明的概念。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。
基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本实施例保护的范围。在不冲突的情况下，下述的实施例及实施例中的特征可以相互组合。
如图1所示，本实施例提供的岩层参数的确定方法，包括：
S1、获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf；
示例性的，在钻机对岩层进行钻探的过程中，通过实验装置可以直接或者间接的获取单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf这类的实验数据；当实验装置在获取这些数据时，可以同时记录下获取这些数据时的时间，这些数据获取的时间可以是同步的。
其中，g个试验数据组表示多个试验数据组，并不是对数据组的个数作具体的限定，只是为了后续的表述方便，对获取的试验数据组的个数所做的表达。
S2、将g个试验数据组中的试验数据分别代入岩石单轴抗压强度值Rc的估算公式，以确定在不同的数据采集时间t采集到的g个试验数据组分别对应的g个岩石单轴抗压强度值Rci，其中Rci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；
示例性的，岩石单轴抗压强度值Rc的估算公式是预先已经确定的，其中岩石单轴抗压强度值Rc的估算公式可以是直接用其他文献中已经确定的参考公式，也可以通过对另外的实验数据进行线性回归来得到。
S3、将确定的g个岩石单轴抗压强度值Rci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(Rc1，Rc2，Rc3，……，Rci，……，Rcg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；最终分类结果包括最佳分类数ds和每一分类的子序列；其中最佳分类数ds即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为Ae＝(Rca，Rca+1，……，Rcb)时，则第e层岩层的单轴抗压强度均值R为：
示例性的，由于在钻探过程中，随着钻探深度的增加，岩层的岩性发生变化，在这个过程中所获取的g个岩石单轴抗压强度值Rci自身具有时序性，为了得到连续的岩层的岩性，当对这些数据进行分析的时候，也有时序性的条件限制。子序列为数据序列A中的某一段序列，且多个子序列之间不出现重叠的情况。在最终分类结果中，数据序列A分为ds个子序列，且将多各子序列按照时间顺序排列后，其数据的连续性和数据序列A可以完全相同。
K-means算法是硬聚类算法，是典型的基于原型的目标函数聚类方法的代表，它是数据点到原型的某种距离作为优化的目标函数，利用函数求极值的方法得到迭代运算的调整规则。算法采用误差平方和准则函数作为聚类准则函数。具有K-means算法功能的软件常用的有R软件，Matlab。聚类分析的具体过程是依靠现有的软件完成的，本实施例不对聚类分析的过程和原理进行详细的描述，仅对所需输入的原始数据和所要得到的数据类型作示例性说明。
可选地，在上述实施例的基础上，本实施例提供的岩层参数的确定方法还包括以下步骤：
S4、获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；确定第e个子序列Ae＝(Rca，Rca+1，……，Rcb)中岩石单轴抗压强度值所对应的数据采集时间t的区间Te＝[ta，tb]；则第e层岩层的厚度为：
示例性的，受岩石单轴抗压强度的影响，钻杆在岩层中的钻速会发生变化，钻速传感器可以每隔一定时间，比如几毫秒，记录一次钻杆的钻速值，根据实验中采集得到的某个时间段中不同时间对应的多个钻速值，通过非线性拟合，可以得到随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)。
也可以在得到数据采集时间t的区间Te＝[ta，tb]的前提下，结合在此时间段钻探的总深度，将钻速的平均值作为钻杆的钻进速度V。
可选地，在上述实施例的基础上，步骤S1还可以包括：根据公式(3)计算钻进单位体积岩石时的钻杆轴力做功值wF：
根据公式(4)计算钻进单位体积岩石时的钻杆扭矩做功值wM：
根据公式(5)计算钻进单位体积岩石时的钻头与孔底摩擦做功值wf：
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
示例性的，有些测试装置可以直接得到钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf；而在实际中，也可以通过间接测量的方式来得到。
可选地，在上述实施例的基础上，步骤S2还可以包括：本实施例提供的岩层参数的确定方法中岩石单轴抗压强度值Rc的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf为自变量，以岩石单轴抗压强度值Rc为因变量，通过线性回归方法得到。
示例性的，线性回归方法时科学研究中常用的一种数据处理方法，即通过预先获取可靠的实验数据来得到经验公式，然后利用该经验公式对结果进行估算的一种方法。多元线性回归的软件可以使用SPSS(Statistical Product and Service Solutions，统计产品与服务解决方案)、Microsoft Excel等可以进行数据处理的软件。
此外，如图2所示，本实施例还提供了岩层参数的确定装置，包括获取模块1，获取模块1用于获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf；还包括确定模块2，确定模块2用于将g个试验数据组中的试验数据分别代入岩石单轴抗压强度值Rc的估算公式，以确定在不同的数据采集时间t采集到的g个试验数据组分别对应的g个岩石单轴抗压强度值Rci，其中Rci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；确定模块2还用于将确定的g个岩石单轴抗压强度值Rci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(Rc1，Rc2，Rc3，……，Rci，……，Rcg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；最终分类结果包括最佳分类数ds和每一分类的子序列；其中最佳分类数ds即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为Ae＝(Rca，Rca+1，……，Rcb)时，则第e层岩层的单轴抗压强度均值R为：
可选地，本实施例提供的岩层参数的确定装置中获取模块1还用于获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；确定第e个子序列Ae＝(Rca，Rca+1，……，Rcb)中岩石单轴抗压强度值所对应的数据采集时间t的区间Te＝[ta，tb]；则第e层岩层的厚度为：
可选地，本实施例提供的岩层参数的确定装置中，获取模块1还用于根据公式(3)获取钻进单位体积岩石时的钻杆轴力做功值wF：
获取模块1还用于根据公式(4)获取钻进单位体积岩石时的钻杆扭矩做功值wM：
获取模块1还用于根据公式(5)获取钻进单位体积岩石时的钻头与孔底摩擦做功值wf：
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
示例性的，钻杆的推力F可以通过钻机与地面接触的压力减去钻机自身的重力得到；扩孔系数λ为钻杆成孔后钻孔横截面与钻杆横截面的面积的比值；钻杆的扭矩M可以通过钻机上连接的扭矩传感器14测量得到；钻杆的转速n可以通过钻机上连接的转速传感器12测量得到；钻杆的钻进速度V可以通过钻机上连接的钻速传感器11测量得到。
可选地，本实施例提供的岩层参数的确定装置中确定模块2还用于确定岩石单轴抗压强度值Rc的估算公式；岩石单轴抗压强度值Rc的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf为自变量，以岩石单轴抗压强度值Rc为因变量，通过线性回归方法得到。
可选地，本实施例提供的岩层参数的确定装置还包括数据采集器，数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；数据采集器包括钻速传感器11、转速传感器12、压力传感器13和扭矩传感器14；钻速传感器11用于采集钻进速度V；转速传感器12用于采集钻杆的转速n；压力传感器13用于采集钻杆的推力F；扭矩传感器14用于采集钻杆的扭矩M；还包括显示器，显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块1获取的数据和确定模块2确定的数据；获取模块1用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值Rc；确定模块2对获取模块1获取的数据信息进行处理。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质。该程序在执行时，执行包括上述个方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
有学者在实验室进行了相关研究，采用直径为60mm的PDC钻头进行钻进实验。钻进试件分为不同强度的砂浆试件28组，编号J1-J28，砂岩试件8组，编号为S1-S8。如表1所示为其中的30组实验数据。
以表1中数据为基础，通过岩石单轴抗压强度值Rc的估算模型预设为公式(1)进行多元线性回归，得到岩石单轴抗压强度值Rc的估算公式：
表1钻进实验数据
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (9)
1.一种岩层参数的确定方法，其特征在于，包括：
获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，所述试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf；
将所述g个试验数据组中的试验数据分别代入岩石单轴抗压强度值Rc的估算公式，以确定在不同的数据采集时间t采集到的所述g个试验数据组分别对应的g个岩石单轴抗压强度值Rci，其中Rci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；
将确定的g个岩石单轴抗压强度值Rci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(Rc1，Rc2，Rc3，……，Rci，……，Rcg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；所述最终分类结果包括最佳分类数ds和每一分类的子序列；其中最佳分类数ds即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为Ae＝(Rca，Rca+1，……，Rcb)时，则第e层岩层的单轴抗压强度均值R为：
2.根据权利要求1所述的确定方法，其特征在于，还包括：
获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，所述钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；所述钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；
确定第e个子序列Ae＝(Rca，Rca+1，……，Rcb)中岩石单轴抗压强度值所对应的数据采集时间t的区间Te＝[ta，tb]；
则第e层岩层的厚度为：
3.根据权利要求1所述的确定方法，其特征在于：
根据公式(3)计算所述钻进单位体积岩石时的钻杆轴力做功值wF：
根据公式(4)计算所述钻进单位体积岩石时的钻杆扭矩做功值wM：
根据公式(5)计算所述钻进单位体积岩石时的钻头与孔底摩擦做功值wf：
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
4.根据权利要求1所述的确定方法，其特征在于，所述岩石单轴抗压强度值Rc的估算公式的确定方法包括：
以钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf为自变量，以岩石单轴抗压强度值Rc为因变量，通过线性回归方法得到。
5.一种岩层参数的确定装置，其特征在于，包括：
获取模块，所述获取模块用于获取钻进岩层过程中在不同的数据采集时间t采集到的g个试验数据组，所述试验数据组中的试验数据包括钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf；
确定模块，所述确定模块用于将所述g个试验数据组中的试验数据分别代入岩石单轴抗压强度值Rc的估算公式，以确定在不同的数据采集时间t采集到的所述g个试验数据组分别对应的g个岩石单轴抗压强度值Rci，其中Rci为按照采集时间顺序获取的第i个试验数据组所对应的岩石单轴抗压强度值；
所述确定模块还用于将确定的g个岩石单轴抗压强度值Rci按照数据采集的时间顺序形成一岩石单轴抗压强度值的数据序列A＝(Rc1，Rc2，Rc3，……，Rci，……，Rcg)，将该数据序列A通过K-means聚类分析软件来确定最终分类结果；所述最终分类结果包括最佳分类数ds和每一分类的子序列；其中最佳分类数ds即为岩层按照岩石单轴抗压强度值的区别划分的总层数；当第e个子序列为Ae＝(Rca，Rca+1，……，Rcb)时，则第e层岩层的单轴抗压强度均值R为：
6.根据权利要求5所述的确定装置，其特征在于：
所述获取模块还用于获取随数据采集时间t变化的钻杆的钻进速度V的表达式V(t)，所述钻杆的钻进速度为钻杆每秒进入岩层的深度，单位为m/s；所述钻杆的钻进速度V的表达式V(t)的获取方法包括：将不同的钻杆的钻进速度V与钻进速度V的采集时间t通过数据拟合得到表达式V(t)；
确定第e个子序列Ae＝(Rca，Rca+1，……，Rcb)中岩石单轴抗压强度值所对应的数据采集时间t的区间Te＝[ta，tb]；
则第e层岩层的厚度为：
7.根据权利要求5所述的确定装置，其特征在于：
所述获取模块还用于根据公式(3)获取所述钻进单位体积岩石时的钻杆轴力做功值wF：
所述获取模块还用于根据公式(4)获取所述钻进单位体积岩石时的钻杆扭矩做功值wM：
所述获取模块还用于根据公式(5)获取所述钻进单位体积岩石时的钻头与孔底摩擦做功值wf：
式中，F为钻杆的推力，单位为N；λ为扩孔系数，即钻孔的横截面的面积与钻杆的横截面的面积之比；r为钻杆的钻头直径，单位为m；M为钻杆的扭矩，单位为N·m；n为钻杆的转速，单位为r/min；V为钻杆的钻进速度，即为钻杆每秒进入岩层的深度，单位为m/s；μ为钻头与孔底的摩擦系数，μ取值为0.21。
8.根据权利要求5所述的确定装置，其特征在于：
所述确定模块还用于确定所述岩石单轴抗压强度值Rc的估算公式；
所述岩石单轴抗压强度值Rc的估算公式的确定方法包括：以钻进单位体积岩石时的钻杆轴力做功值wF、钻进单位体积岩石时的钻杆扭矩做功值wM、钻进单位体积岩石时的钻头与孔底摩擦做功值wf为自变量，以岩石单轴抗压强度值Rc为因变量，通过线性回归方法得到。
9.根据权利要求5所述的确定装置，其特征在于，还包括：
数据采集器，所述数据采集器用于采集钻杆钻进单位体积岩石时的试验数据；所述数据采集器包括钻速传感器、转速传感器、压力传感器和扭矩传感器；所述钻速传感器用于采集钻进速度V；所述转速传感器用于采集钻杆的转速n；所述压力传感器用于采集钻杆的推力F；所述扭矩传感器用于采集钻杆的扭矩M；
显示器，所述显示器用于显示数据采集器采集的钻杆钻进单位体积岩石时的试验数据、获取模块获取的数据和确定模块确定的数据；
所述获取模块用于接收数据采集器采集的钻杆钻进单位体积岩石时的试验数据和对应的单轴抗压强度值Rc；所述确定模块对获取模块获取的数据信息进行处理。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910255972.8A CN110130883A (zh) | 2019-04-01 | 2019-04-01 | 岩层参数的确定方法及装置 |
| AU2019439997A AU2019439997A1 (en) | 2019-04-01 | 2019-09-03 | Method and device for determining rock stratum parameters |
| RU2020141122A RU2762675C9 (ru) | 2019-04-01 | 2019-09-03 | Способ и устройство для определения параметров пласта горной породы |
| PCT/CN2019/104107 WO2020199495A1 (zh) | 2019-04-01 | 2019-09-03 | 岩层参数的确定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910255972.8A CN110130883A (zh) | 2019-04-01 | 2019-04-01 | 岩层参数的确定方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110130883A true CN110130883A (zh) | 2019-08-16 |
Family
ID=67569123
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910255972.8A Pending CN110130883A (zh) | 2019-04-01 | 2019-04-01 | 岩层参数的确定方法及装置 |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN110130883A (zh) |
| AU (1) | AU2019439997A1 (zh) |
| WO (1) | WO2020199495A1 (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110486007A (zh) * | 2019-08-29 | 2019-11-22 | 武汉长盛煤安科技有限公司 | 煤矿随钻围岩力学参数原位测试装置及方法 |
| CN111075424A (zh) * | 2019-12-25 | 2020-04-28 | 中国石油大学(华东) | 一种修正随钻测量参数测量结果的方法 |
| CN111706322A (zh) * | 2020-07-17 | 2020-09-25 | 中国铁建重工集团股份有限公司 | 一种岩石钻进响应的预测方法和预测系统 |
| WO2020199495A1 (zh) * | 2019-04-01 | 2020-10-08 | 中国矿业大学 | 岩层参数的确定方法及装置 |
| CN111855961A (zh) * | 2020-07-24 | 2020-10-30 | 中南大学 | 一种岩体钻进质量检测方法、钻机、服务器、存储介质 |
| RU2762675C9 (ru) * | 2019-04-01 | 2022-08-04 | Чайна Юниверсити Оф Майнинг Энд Текнолоджи | Способ и устройство для определения параметров пласта горной породы |
| CN117763466A (zh) * | 2024-02-22 | 2024-03-26 | 中石化经纬有限公司 | 一种基于聚类算法的地层可钻性评价方法及系统 |
| CN118392662A (zh) * | 2024-06-26 | 2024-07-26 | 中国科学院武汉岩土力学研究所 | 基于岩芯回弹和随钻参数的岩体抗压强度测定方法及系统 |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112412394B (zh) * | 2020-11-11 | 2022-12-06 | 安徽理工大学 | 一种钻孔分层充填方法 |
| CN113482599A (zh) * | 2021-07-08 | 2021-10-08 | 中国科学院武汉岩土力学研究所 | 便携式原位测试与识别数据采集系统及设备 |
| CN116485225B (zh) * | 2023-03-15 | 2023-11-10 | 西南交通大学 | 基于钻进参数的施工阶段围岩bq值自动获取方法及系统 |
| CN116499827A (zh) * | 2023-04-18 | 2023-07-28 | 中国石油大学(北京) | 多岩性互层人造岩样制备装置及方法 |
| CN117288587B (zh) * | 2023-11-24 | 2024-02-20 | 中国矿业大学(北京) | 岩体抗拉强度随钻测试方法与系统 |
Citations (5)
| 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施工围岩分级方法 |
| CN108519284A (zh) * | 2018-03-30 | 2018-09-11 | 山东大学 | 地下工程围岩数字钻探实时分级方法 |
| CN108663269A (zh) * | 2018-03-30 | 2018-10-16 | 山东大学 | 基于等效岩体强度的地下工程围岩数字钻探分区方法 |
| CN108760487A (zh) * | 2018-04-16 | 2018-11-06 | 西安科技大学 | 基于锚杆受力分析的深埋洞室围岩分区破裂演化分析方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6672386B2 (en) * | 2002-06-06 | 2004-01-06 | Baker Hughes Incorporated | Method for in-situ analysis of formation parameters |
| CN101095143B (zh) * | 2004-11-04 | 2010-06-16 | 贝克休斯公司 | 多精度多维测井数据反演以及深岩层成像方法 |
| CN102587869A (zh) * | 2005-05-25 | 2012-07-18 | 地质力学国际公司 | 分析和控制在钻孔中水锤产生的波传播的方法和装置 |
| US8274289B2 (en) * | 2006-12-15 | 2012-09-25 | Halliburton Energy Services, Inc. | Antenna coupling component measurement tool having rotating antenna configuration |
| US8768627B2 (en) * | 2011-03-11 | 2014-07-01 | Landmark Graphics Corporation | Methods and systems of estimating formation parameters |
| CN110130883A (zh) * | 2019-04-01 | 2019-08-16 | 中国矿业大学 | 岩层参数的确定方法及装置 |
-
2019
- 2019-04-01 CN CN201910255972.8A patent/CN110130883A/zh active Pending
- 2019-09-03 AU AU2019439997A patent/AU2019439997A1/en not_active Abandoned
- 2019-09-03 WO PCT/CN2019/104107 patent/WO2020199495A1/zh active Application Filing
Patent Citations (5)
| 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施工围岩分级方法 |
| CN108519284A (zh) * | 2018-03-30 | 2018-09-11 | 山东大学 | 地下工程围岩数字钻探实时分级方法 |
| CN108663269A (zh) * | 2018-03-30 | 2018-10-16 | 山东大学 | 基于等效岩体强度的地下工程围岩数字钻探分区方法 |
| CN108760487A (zh) * | 2018-04-16 | 2018-11-06 | 西安科技大学 | 基于锚杆受力分析的深埋洞室围岩分区破裂演化分析方法 |
Non-Patent Citations (1)
| Title |
|---|
| 张泽荫: "松散矿石流动性分析及其系统聚类法分类", 《金属矿山》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020199495A1 (zh) * | 2019-04-01 | 2020-10-08 | 中国矿业大学 | 岩层参数的确定方法及装置 |
| RU2762675C1 (ru) * | 2019-04-01 | 2021-12-21 | Чайна Юниверсити Оф Майнинг Энд Текнолоджи | Способ и устройство для определения параметров пласта горной породы |
| RU2762675C9 (ru) * | 2019-04-01 | 2022-08-04 | Чайна Юниверсити Оф Майнинг Энд Текнолоджи | Способ и устройство для определения параметров пласта горной породы |
| CN110486007A (zh) * | 2019-08-29 | 2019-11-22 | 武汉长盛煤安科技有限公司 | 煤矿随钻围岩力学参数原位测试装置及方法 |
| CN110486007B (zh) * | 2019-08-29 | 2023-02-03 | 武汉长盛煤安科技有限公司 | 煤矿随钻围岩力学参数原位测试装置及方法 |
| CN111075424A (zh) * | 2019-12-25 | 2020-04-28 | 中国石油大学(华东) | 一种修正随钻测量参数测量结果的方法 |
| CN111706322A (zh) * | 2020-07-17 | 2020-09-25 | 中国铁建重工集团股份有限公司 | 一种岩石钻进响应的预测方法和预测系统 |
| CN111855961A (zh) * | 2020-07-24 | 2020-10-30 | 中南大学 | 一种岩体钻进质量检测方法、钻机、服务器、存储介质 |
| CN117763466A (zh) * | 2024-02-22 | 2024-03-26 | 中石化经纬有限公司 | 一种基于聚类算法的地层可钻性评价方法及系统 |
| CN118392662A (zh) * | 2024-06-26 | 2024-07-26 | 中国科学院武汉岩土力学研究所 | 基于岩芯回弹和随钻参数的岩体抗压强度测定方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| AU2019439997A1 (en) | 2021-01-07 |
| RU2762675C1 (ru) | 2021-12-21 |
| WO2020199495A1 (zh) | 2020-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110130883A (zh) | 岩层参数的确定方法及装置 | |
| CN104329079B (zh) | 一种识别气测录井油气层的方法以及系统 | |
| CN104514552B (zh) | 一种煤层气储层识别与丰度预测的方法 | |
| CN110644980A (zh) | 一种超低渗透油藏储层综合分类评价方法 | |
| WO2009129083A1 (en) | Continuous measurement of heterogeneity of geomaterials | |
| CN104747183A (zh) | 一种碳酸盐岩储层综合分类方法 | |
| CN105697002A (zh) | 一种用于识别煤系地层岩性的方法 | |
| CN108561126B (zh) | 一种确定页岩气储层有机孔隙度的简易方法 | |
| CN101750633A (zh) | 利用核磁共振录井解释参数进行储层评价的方法 | |
| CN107505344A (zh) | 利用“最小二乘积”法的岩性解释方法 | |
| CA2717596C (en) | Pay zone prediction | |
| CN104153768A (zh) | 一种评价花岗岩储层储集性能的方法 | |
| CN104428492A (zh) | 使用矿物学确定压裂刺激点的布置的系统、方法、和计算机程序产品 | |
| CN104374827A (zh) | 横观各向同性岩体原位动弹性模量的各向异性系数的测量方法 | |
| CN115030707A (zh) | 一种油页岩“甜点”快速评价方法 | |
| CN105064987B (zh) | 利用随钻录井q参数进行水层识别的解释评价方法 | |
| CN103306672B (zh) | 一种预测页岩地层不同钻进方向的研磨性的方法 | |
| CN104463686A (zh) | 利用判别分析法随钻识别页岩气储层的方法 | |
| CN113323647B (zh) | 用于监测岩石抗压强度和可钻性的方法及装置 | |
| CN105182421A (zh) | 一种定量评价地层脆性特征的方法 | |
| Sarkheil et al. | The fracture network modeling in naturally fractured reservoirs using artificial neural network based on image loges and core measurements | |
| CN106568918A (zh) | 页岩有机碳含量toc预测方法 | |
| CN112184033B (zh) | 一种碳酸盐岩岩层断层封闭性评价方法 | |
| CN108661629A (zh) | 一种针对页岩地层的工程甜点定量评价方法 | |
| US20210080605A1 (en) | Methods for the Regional-Scale Simulation of Hydraulic Fracturing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190816 |