CN112906148B - 一种煤巷穿层瓦斯防治钻孔设计计算方法 - Google Patents
一种煤巷穿层瓦斯防治钻孔设计计算方法 Download PDFInfo
- Publication number
- CN112906148B CN112906148B CN202110068860.9A CN202110068860A CN112906148B CN 112906148 B CN112906148 B CN 112906148B CN 202110068860 A CN202110068860 A CN 202110068860A CN 112906148 B CN112906148 B CN 112906148B
- Authority
- CN
- China
- Prior art keywords
- drilling
- cosβsinα
- sinβ
- cosβcosα
- cosγsinα
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
本发明提供了一种煤巷穿层瓦斯防治钻孔设计计算方法，所述计算方法包括：S1、预抽煤巷穿层钻孔的设计模型确定以及设计；S11、穿层钻孔设计：依据数据库、规则库、输出参数进行设计；S12、各钻场终孔点坐标；S13、各钻场开孔点坐标；S2、穿层钻孔设计模型确定以及设计；S21：穿层钻孔设计规则；S22：穿层钻孔计算。本发明基于动态三维瓦斯地质模型，建立瓦斯防治钻孔信息数据库、提出瓦斯防治无人化智能钻孔设计原则与方法，实现了三维瓦斯地质模型的自动建立与精准纠正、抽采单元的自动划分与更新、钻孔智能设计与动态调整、钻孔合理性评判以及智能补救。
Description
技术领域
本发明涉及煤矿瓦斯抽采钻孔技术领域，特别涉及一种煤巷穿层瓦斯防治钻孔设计计算方法。
背景技术
瓦斯防治钻孔的设计考虑多个方面的因素，一是煤岩层结构、倾角、构造等地质特征；二是煤层透气性、瓦斯含量等瓦斯赋存特征；三是煤巷、底板巷或顶板巷道布置方式；四是回采设计、采掘接替时空要求等相关的采掘计划；五是瓦斯治理过程中发现的局部不达标或钻孔不到位问题。
目前瓦斯防治钻孔的设计主要是技术人员依据以上因素进行人工设计，利用CAD等软件绘制钻孔布置图，其效率低、工作量大，对于复杂的抽采设计，则不能很好反映设计思想和清楚表达钻孔工程与岩层之间的相互关系。在进行抽采设计时，设计人员首先选用设计模版或者惯用的设计参数，而这种选择不一定是当前条件下的最优设计参数，甚至不符合煤矿井下复杂变化的赋存条件，导致产生空白带、不符合规程等不良后果。另外，人工设计必然带来一些问题，首先，钻孔设计时考虑的因素较多，人工设计不可避免会出现遗漏某种因素的问题，特别是面临构造等复杂地质条件时，人工设计难以很好的把握构造影响范围，尤其是三维尺度的影响范围；其次，实际煤层的赋存状态不是均一的，三维尺度上煤层的赋存状态更为复杂，人工设计时普遍忽略煤层三维尺度的变化，这就导致钻孔施工的轨迹出现与设计轨迹不一致的现象；再次，钻孔施工过程受到煤岩层物理性质差异的影响，也会受到钻机钻杆本身性能方面的影响，钻孔实际轨迹会出现偏移现象，而目前人工钻孔施工还无法实现钻孔轨迹的自动修正，也无法回馈钻孔轨迹参数。这就需要根据前期施工钻孔的实际抽采特征进行智能分析或诊断，及时修正钻孔设计或进行补救钻孔设计等工作。
这些问题的出现，不是技术人员的专业技术能力的问题，更多是现有技术水平的限制而导致的。特别是考虑众多因素约束条件时的钻孔合理设计、前期施工的钻孔群瓦斯抽采数据的智能分析、钻孔实际轨迹的三维成像及合理性评判等，均需要极大的工作量，而人工智能则成为最可能的解决办法。但是，目前还没有智能钻孔设计的先例，特别是针对含有瓦斯气体的复杂煤层环境。
发明内容
为解决上述问题，本发明旨在提出一种煤巷穿层瓦斯防治钻孔设计计算方法，通过数据库、规则库以及输出参数确定设计模型，实现三维瓦斯地质模型的自动建立与精准纠正、抽采单元的自动划分与更新、钻孔智能设计与动态调整、钻孔合理性评判与智能补救等。
为达到上述目的，本发明的技术方案是这样实现的：
一种煤巷穿层瓦斯防治钻孔设计计算方法，所述计算方法包括：
S1、预抽煤巷穿层钻孔的设计模型确定以及设计；
S11、穿层钻孔设计：
依据数据库、规则库、输出参数进行设计；
具体计算方法见实施例；
S12、各钻场终孔点坐标，具体计算方法见实施例：
S13、各钻场开孔点坐标，具体计算方法见实施例。
进一步的，所述计算方法还包括：
S2、穿层钻孔设计模型确定以及设计：
S21：穿层钻孔设计规则
(1)所有设计的钻孔孔长需要≤钻机的实际施工能力，以钻机最大的施工能力作为单条底板巷的控制范围；
(2)三条底板巷设计的钻孔，均按照最大的控制范围设计，当存在部分钻孔设计重叠，此时需要钻孔优化，相邻底板巷重叠的钻孔，按照钻孔长度进行比较，优选重叠部分钻孔长度较短的进行施工；
S22：穿层钻孔计算，具体计算方法见实施例。
进一步的，所述S11中穿层钻孔设计中的数据库包括空间坐标确定、抽采参数、煤层厚度、煤层倾角、瓦斯含量、瓦斯压力、施工钻场；规则库包括钻孔倾角、钻孔方位角、巷道方位角、钻孔长度计算、钻孔开孔位置、钻孔终孔位置、钻孔数量确定；输出参数包括钻孔倾角、钻孔方位角、巷道方位角、钻孔长度、钻孔开孔点坐标、钻孔终孔点坐标、钻孔数量。
有益效果：本发明基于动态三维瓦斯地质模型，建立瓦斯防治钻孔信息数据库、提出瓦斯防治无人化智能钻孔设计原则与方法，实现了三维瓦斯地质模型的自动建立与精准纠正、抽采单元的自动划分与更新、钻孔智能设计与动态调整、钻孔合理性评判以及智能补救。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明实施例所述的煤巷穿层瓦斯防治钻孔设计计算方法的流程图；
图2为本发明实施例所述的预抽煤巷钻孔控制范围示意图；
图3为本发明实施例所述的各钻场终孔点位置平面示意图；
图4为本发明实施例所述的穿层钻孔钻场开孔点位置平面示意图；
图5为本发明实施例所述的钻场控制点几何坐标位置示意图；
图6为本发明实施例所述的钻场内开孔点位置平面示意图；
图7为本发明实施例所述的目标煤层理想状态下的空间几何坐标示意图；
图8为本发明实施例所述的底板巷设计平面图；
图9为本发明实施例所述的底板巷钻孔设计剖面图；
图10为本发明实施例所述的第一条巷道钻孔布置剖面图；
图11为本发明实施例所述的第二条巷道钻孔布置剖面图；
图12为本发明实施例所述的第三条巷道钻孔布置剖面图。
具体实施方式
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面将参考附图并结合实施例来详细说明本发明。
实施例1
参见图1-12：一种煤巷穿层瓦斯防治钻孔设计计算方法，所述计算方法包括：
S1、预抽煤巷穿层钻孔的设计模型确定以及设计；
S11、穿层钻孔设计：
依据数据库、规则库、输出参数进行设计；所述S11中穿层钻孔设计中的数据库包括空间坐标确定、抽采参数、煤层厚度、煤层倾角、瓦斯含量、瓦斯压力、施工钻场；规则库包括钻孔倾角、钻孔方位角、巷道方位角、钻孔长度计算、钻孔开孔位置、钻孔终孔位置、钻孔数量确定；输出参数包括钻孔倾角、钻孔方位角、巷道方位角、钻孔长度、钻孔开孔点坐标、钻孔终孔点坐标、钻孔数量。
假定在理想条件下，施工穿层钻孔的钻场所在的巷道，和目标煤层的倾角，煤层沿走向夹角，以及方位角的都相同；
预抽煤巷条带沿工作面走向长：Lz＝LAB，m
预抽煤巷条带宽度：L宽＝LAD，m
每个钻场钻孔个数，个；
n＝i*j
式中，i为每个钻场沿走向钻孔个数，个；
j为每个钻场沿倾向钻孔个数，个；
巷道方位角α：
煤层走向与水平面夹角β：
煤层倾角γ：
预抽区段煤层瓦斯区域防突措施的钻孔应当控制区段内的整个开采块段、两侧回采巷道及其外侧一定范围内的煤层；要求钻孔控制回采巷道外侧的范围是：倾斜、急倾斜煤层巷道上帮轮廓线外至少20m，下帮至少10m；其他为巷道两侧轮廓线外至少各15m；考虑理想条件下，在预抽煤层的整个煤层工作面外围增加Lp的抽采范围；
参见图2为预抽煤巷钻孔控制范围示意图；那么控制点坐标可重新计算确定：
A1：(x1′，y1′，z1′)
x1′＝x1-Lp cosβcosα+Lp cosγsinα
y1′＝y1-Lp cosβsinα-Lp cosγsinα
z1′＝z1-Lp sinβ-Lp sinγ
B1：(x2′，y2′，z2′)
x2′＝x2+Lp cosβcosα+Lp cosγsinα
y2′＝y2+Lp cosβsinα-Lp cosγsinα
z2′＝z2+Lp sinβ-Lp sinγ
C1：(x3′，y3′，z3′)
x3′＝x3+Lp cosβcosα-Lp cosγsinα
y3′＝y3+Lp cosβsinα+Lp cosγsinα
z3′＝z3+Lp sinβ+Lp sinγ
D1：(x4′，y4′，z4′)
x4′＝x4-Lp cosβcosα-Lp cosγsinα
y4′＝y4-Lp cosβsinα+Lp cosγsinα
z4′＝z4-Lp sinβ+Lp sinγ
加上预抽控制范围之后确定的走向长度：
加上预抽控制范围之后确定的预抽煤巷条带宽度：
钻场个数：
S12、各钻场终孔点坐标：
参见图3为各钻场终孔点位置平面示意图；第1个钻场(j，i)#钻孔终孔点坐标：
xj，i＝x1′-(j-1)d2 cosγsinα+(i-1)d1 cosβcosα
yj，i＝y1′+(j-1)d2 cosγcosα+(i-1)d1 cosβsinα
zj，i＝z1′+(j-1)d2 sinγ+(i-1)d1 sinβ
第2个钻场(1，i)#钻孔终孔点坐标：
x1，i＝x1′+(2i-1)d1 cosβcosα
y1，i＝y1′+(2i-1)d1 cosβsinα
z1，i＝z1′+(2i-1)d1 sinβ
第2个钻场(j，i)#钻孔终孔点坐标：
xj，i＝x1′+(2i-1)d1 cosβcosα-(j-1)d2 cosγsinα
yj，i＝y1′+(2i-1)d1 cosβsinα+(j-1)d2 cosγcosα
zj，i＝z1′+(2i-1)d1 sinβ+(j-1)d2 sinγ
第m个钻场(j，i)#钻孔终孔点坐标：
xj，i＝x1′+(mi-1)d1 cosβcosα-(j-1)d2 cosγsinα
yj，i＝y1′+(mi-1)d1 cosβsinα+(j-1)d2 cosγcosα
zj，i＝z1′+(mi-1)d1 sinβ+(j-1)d2 sinγ
S13、各钻场开孔点坐标：
参见图4为穿层钻孔钻场开孔点位置平面示意图，图5为钻场控制点几何坐标位置示意图；
图中：A2＝E1
假定第1个钻场的控制点坐标E1：(x0，y0，z0)
那么每个钻场的坐标可表示为：
第2个钻场控制点位置坐标：E2：(xe2，ye2，ze2)
xe2＝x0+id1 cosβcosα
ye2＝yn+id1 cosβsinα
ze2＝z0+id1 sinβ
第m个钻场控制点位置坐标Em：(xem，yem，zem)
xem＝x0+(m-1)id1 cosβcosα
yem＝y0+(m-1)id1 cosβsinα
zem＝z0+(m-1)id1 sinβ
参见图6为钻场内开孔点位置平面示意图；第1个钻场(j，i)#钻孔开孔点坐标：
xj，i′＝x0-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
yj，i′＝y0+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
zj，i′＝z0+(j-1)d4 sinγ+(i-1)d3 sinβ
第2个钻场(j，1)#钻孔开孔点坐标：
xj，1′＝x0+id1 cosβcosα-(j-1)d4 cosγsinα
yj，1′＝y0+id1 cosβsinα+(j-1)d4 cosγcosα
zj，1′＝z0+id1 sinβ+(j-1)d4 sinγ
第2个钻场(j，i)#钻孔开孔点坐标：
xj，i′＝xe2-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
＝x0+id1 cosβcosα-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
＝x0+[id1+(i-1)d3]cosβcosα-(j-1)d4 cosγsinα
yj，i′＝ye2+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
＝y0+id1 cosβsinα+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
＝y0+[id1+(i-1)d3]cosβsinα+(j-1)d4 cosγcosα#
zj，i′＝ze2+(j-1)d4 sinγ+(i-1)d3 sinβ
＝z0+id1 sinβ+(j-1)d4 sinγ+(i-1)d3 sinβ
＝z0+[id1+(i-1)d3]sinβ+(j-1)d4 sinγ#
第m钻场(j，i)#钻孔开孔点坐标：
xj，i′＝xem-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
＝x0+(m-1)id1 cosβcosα-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
＝x0+[(m-1)id1+(i-1)d3]cosβcosα-(j-1)d4 cosγsinα#
yj，i′＝yem+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
＝y0+(m-1)id1 cosβsinα+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
＝y0+[(m-1)id1+(i-1)d3]cosβsinα+(j-1)d4 cosγcosα#
zj，i′＝zem+(j-1)d4 sinγ+(i-1)d3 sinβ
＝z0+(m-1)id1 sinβ+(j-1)d4 sinγ+(i-1)d3 sinβ
＝z0+[(m-1)id1+(i-1)d3]sinβ+(j-1)d4 sinγ#
各钻孔的钻孔长度：L孔长
各钻孔的钻孔方位角：α1
本实施例通过依据瓦斯防治钻孔信息数据库、规则库以及输出参数，建立预抽煤巷穿层钻孔的设计模型，从而构建动态的三维瓦斯地质模型。
进一步的，所述计算方法还包括：
S2、穿层钻孔设计模型确定以及设计：
S21：穿层钻孔设计规则
(3)所有设计的钻孔孔长需要≤钻机的实际施工能力，以钻机最大的施工能力作为单条底板巷的控制范围；
(4)三条底板巷设计的钻孔，均按照最大的控制范围设计，当存在部分钻孔设计重叠，此时需要钻孔优化，相邻底板巷重叠的钻孔，按照钻孔长度进行比较，优选重叠部分钻孔长度较短的进行施工；
S22：穿层钻孔计算
A：(x1，y1，z1)，控制点坐标A；
B：(x2，y2，z2)，控制点坐标B；
C：(x3，y3，z3)，控制点坐标C；
D：(x4，y4，z4)，控制点坐标D；
d：钻孔间距
Lz：预抽煤层工作面沿走向长：Lz＝LAB，m
L面长：预抽煤层工作面长度：L面长＝LAD，m
L孔长：钻孔长度，m；
α：巷道方位角(单位°)；
β：煤层走向与水平面夹角(单位°)；
γ：：煤层倾角(单位°)；
假定在理想条件下，底板巷和目标煤层的倾角，煤层沿走向夹角，以及方位角的都相同；
巷道方位角α：
煤层走向与水平面夹角β：
煤层倾角γ：
预抽区段煤层瓦斯区域防突措施的钻孔应当控制区段内的整个开采块段外侧一定范围内的煤层；
那么控制点坐标可重新计算确定：
A1：(x1′，y1′，z1′)
x1′＝x1-Lp cosβcosα+Lp cosγsinα
y1′＝y1-Lp cosβsinα-Lp cosγsinα
z1′＝z1-Lpsinβ-Lp sinγ
B1：(x2′，y2′，z2′)
x2′＝x2+Lp cosβcosα+Lp cosγsinα
y2′＝y2+Lp cosβsinα-Lp cosγsinα
z2′＝z2+Lp sinβ-Lp sinγ
C1：(x3′，y3′，z3′)
x3′＝x3+Lp cosβcosα-Lp cosγsinα
y3′＝y3+Lp cosβsinα+Lp cosγsinα
z3′＝z3+Lp sinβ+Lp sinγ
D1：(x4′，y4′，z4′)
x4′＝x4-Lp cosβcosα-Lp cosγsinα
y4′＝y4-Lp cosβsinα+Lp cosγsinα
z4′＝z4-Lp sinβ+Lp sinγ
加上预抽控制范围之后确定的走向长度：
加上预抽控制范围之后确定的煤层工作面长度：
参见图7为目标煤层理想状态下的空间几何坐标示意图；图中各个点的三维坐标为：
A(x1，y1，z1)，B(x2，y2，z2)，C(x3，y3，z3)，D(x4，y4，z4)
A1(x1′，y1′，z1′)，B1(x2′，y2′，z2′)，C1(x3′，y3′，z3′)，D1(x4′，y4′，z4′)
参见图8为底板巷设计平面图；图9为底板巷钻孔设计剖面图；假设三条底板巷的控制点坐标分别为(x0，y0，20)，(x′0，y′0，z′0)，(x″0，y″0，z″0)
每条底板巷内设计施工m组钻孔
每组钻孔设计2n-1个钻孔，且：
且
式中：D为钻机打钻能力，m；
h为底板中心线开孔点到煤层的垂直距离，即开孔点对应的中心点终孔坐标的距离；
设计的钻孔长度最长不超过钻机可以施工的最大长度，假设钻机可以施工的最大长度为80m，则设计中的所有孔长均小于80m；
第一条底板巷内：
参见图10为第一条巷道钻孔布置剖面图。
第m个开孔点坐标：
x0+(m-1)d cosβcosα
y0+(m-1)d cosβsinα
z0+(m-1)d sinβ
那么第m个开孔点对应的第m组终孔点坐标为：
中心点坐标
x0+(m-1)d cosβcosα
y0+(m-1)d cosβsinα
z0+(m-1)d sinβ+h
中心点向下第2个终孔点坐标：
x0+(m-1)d cosβcosα+d cosγcosα
y0+(m-1)d cosβsinα-d cosγsinα
z0+(m-1)d sinβ+h-d sinγ
中心点向下第n个终孔点坐标：
x0+(m-1)d cosβcosα+(n-1)d cosγcosα
y0+(m-1)d cosβsinα-(n-1)d cosγsinα
z0+(m-1)d sinβ+h-(n-1)d sinγ
中心点向上第2个终孔点坐标：
x0+(m-1)d cosβcosα-d cosγcosα
y0+(m-1)d cosβsinα+d cosγsinα
z0+(m-1)d sinβ+h+d sinγ
中心点向上第n个终孔点坐标：
x0+(m-1)d cosβcosα-(n-1)d cosγcosα
y0+(m-1)d cosβsinα+(n-1)d cosγsinα
z0+(m-1)d sinβ+h+(n-1)d sinγ
钻孔长度计算：
且L孔长小于钻机的打钻能力，一般为80m；
第二条底板巷内：
参见图11为第二条巷道钻孔布置剖面图。
第m个开孔点坐标：
x′0+(m-1)d cosβcosα
y′0+(m-1)d cosβsinα
z′0+(m-1)d sinβ
那么第m个开孔点对应的第m组终孔点坐标为：
中心点坐标
x′0+(m-1)d cosβcosα
y′0+(m-1)d cosβsinα
z′0+(m-1)d sinβ+h
中心点向下第2个终孔点坐标：
x′0+(m-1)d cosβcosα+d cosγcosα
y′0+(m-1)d cosβsinα-d cosγsinα
z′0+(m-1)d sinβ+h-d sinγ
中心点向下第n个终孔点坐标：
x′0+(m-1)d cosβcosα+(n-1)d cosγcosα
y′0+(m-1)d cosβsinα-(n-1)d cosγsinα
z′0+(m-1)d sinβ+h-(n-1)d sinγ
中心点向上第2个终孔点坐标：
x′0+(m-1)d cosβcosα-2 cosγcosα
y′0+(m-1)d cosβsinα+d cosγsinα
z′0+(m-1)d sinβ+h+d sinγ
中心点向上第n个终孔点坐标：
x′0+(m-1)d cosβcosα-(n-1)d cosγcosα
y′0+(m-1)d cosβsinα+(n-1)d cosγsinα
z′0+(m-1)d sinβ+h+(n-1)d sinγ
第三条底板巷内：
参见图12为第三条巷道钻孔布置剖面图。
第m个开孔点坐标：
x″0+(m-1)d cosβcosα
y″0+(m-1)d cosβsinα
z″0+(m-1)d sinβ
那么第m个开孔点对应的第m组终孔点坐标为：
中心点坐标
x″0+(m-1)d cosβcosα
y″0+(m-1)d cosβsinα
z′0+(m-1)d sinβ+h
中心点向下第2个终孔点坐标：
x″0+(m-1)d cosβcosα+d cosγcosα
y″0+(m-1)d cosβsinα-d cosγsinα
z″0+(m-1)d sinβ+h-d sinγ
中心点向下第n个终孔点坐标：
x″0+(m-1)d cosβcosα+(n-1)d cosγcosα
y″0+(m-1)d cosβsinα-(n-1)d cosγsinα
z″0+(m-1)d sinβ+h-(n-1)d sinγ
中心点向上第2个终孔点坐标：
x″0+(m-1)d cosβcosα-d cosγcosα
y″0+(m-1)d cosβsinα+d cosγsinα
z″0+(m-1)d sinβ+h+d sinγ
中心点向上第n个终孔点坐标：
x″0+(m-1)d cosβcosα-(n-1)d cosγcosα
y″0+(m-1)d cosβsinα+(n-1)d cosγsinα
z″0+(m-1)d sinβ+h+(n-1)d sinγ
钻孔长度计算：
且L孔长小于钻机的打钻能力，一般为80m；
各钻孔的钻孔方位角α1：
本实施例通过建立瓦斯防治无人化智能钻孔设计规则，从而确定穿层钻孔设计模型，实现了三维瓦斯地质模型的自动建立与精准纠正、抽采单元的自动划分与更新、钻孔智能设计与动态调整、钻孔合理性评判以及智能补救。
需要说明的是，我国煤矿穿层瓦斯钻孔施工普遍采用等间距的钻孔设计方法，其精度差、设计工作量大，不利于矿井的高效、安全生产。因此，研究煤矿瓦斯抽采钻孔的智能设计方法，具有十分重要的意义。本实施例具有以下有益效果：1)技术层面的预期效益：围绕瓦斯防治钻孔需求，进行瓦斯地质单元划分，提出钻孔优选判别原则，构建典型瓦斯抽采单元钻孔设计数据库，基于煤层瓦斯抽采单元划分结论及采掘时空布置，提出区域瓦斯防治钻孔无人化智能设计方法，设计面向钻孔机器人的瓦斯防治钻孔信息数据库；编制钻孔无人化智能设计与分析系统与钻孔机器人数据交互通信协议；开发钻孔无人化智能设计与分析系统。提高瓦斯防治钻孔设计的效率和科学性，提升企业经济效益，促进行业安全技术装备的升级换代，具有巨大的社会和经济效益。2)产业层面的预期指标及效益：在深刻认识主动扰动条件下瓦斯防治钻孔群周围多场耦合规律的基础上，形成瓦斯防治钻孔的智能设计方法，建成钻孔无人化智能设计系统软件平台。通过煤矿现场应用，实现对新技术、新装备的验证、改进，从而引领行业整体创新，推动行业技术全面进步，产生巨大的经济、社会效益。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (3)
1.一种煤巷穿层瓦斯防治钻孔设计计算方法，其特征在于，所述计算方法包括：
S1、预抽煤巷穿层钻孔的设计模型确定以及设计；
S11、穿层钻孔设计：
依据数据库、规则库、输出参数进行设计；
假定在理想条件下，施工穿层钻孔的钻场所在的巷道，和目标煤层的倾角，煤层沿走向夹角，以及方位角的都相同；
预抽煤巷条带沿工作面走向长：Lz＝LAB，单位为m；
预抽煤巷条带宽度：
L宽＝LAD，单位为m；
每个钻场钻孔个数，单位为个；
n＝i*j
式中，i为每个钻场沿走向钻孔个数，单位为个；
j为每个钻场沿倾向钻孔个数，单位为个；
巷道方位角α：
煤层走向与水平面夹角β：
煤层倾角γ：
预抽区段煤层瓦斯区域防突措施的钻孔应当控制区段内的整个开采块段、两侧回采巷道及其外侧一定范围内的煤层；要求钻孔控制回采巷道外侧的范围是：倾斜、急倾斜煤层巷道上帮轮廓线外至少20m，下帮至少10m；其他为巷道两侧轮廓线外至少各15m；考虑理想条件下，在预抽煤层的整个煤层工作面外围增加Lp的抽采范围；
那么控制点坐标可重新计算确定：
A1：(x1′，y1′，z1′)
x1′＝x1-Lp cosβcosα+Lp cosγsinα
y1′＝y1-Lp cosβsinα-Lp cosγsinα
z1′＝z1-Lp sinβ-Lp sinγ
B1：(x2′，y2′，z2′)
x2′＝x2+Lp cosβcosα+Lp cosγsinα
y2′＝y2+Lp cosβsinα-Lp cosγsinα
z2′＝z2+Lp sinβ-Lp sinγ
C1：(x3′，y3′，z3′)
x3′＝x3+Lp cosβcosα-Lp cosγsinα
y3′＝y3+Lp cosβsinα+Lp cosγsinα
z3′＝z3+Lp sinβ+Lp sinγ
D1：(x4′，y4′，z4′)
x4′＝x4-Lp cosβcosα-Lp cosγsinα
y4′＝y4-Lp cosβsinα+Lp cosγsinα
z4′＝z4-Lp sinβ+Lp sinγ
加上预抽控制范围之后确定的走向长度：
加上预抽控制范围之后确定的预抽煤巷条带宽度：
钻场个数：
其中，字符″x1，x2，x3，x4，y1，y2，y3，y4，z1，z2，z3，z4″指的是预抽煤巷条带控制范围的四个顶端控制点A、B、C、D的三维坐标A(x1，y1，z1，)、B(x2，y2，z2，)、C(x3，y3，z3)、D(x4，y4，z4)；
S12、各钻场终孔点坐标：
第1个钻场(j，i)#钻孔终孔点坐标：
xj，i＝x1′-(j-1)d2 cosγsinα+(i-1)d1 cosβcosα
yj，i＝y1′+(j-1)d2 cosγcosα+(i-1)d1 cosβsinα
zj，i＝z1′+(j-1)d2 sinγ+(i-1)d1 sinβ
第2个钻场(1，i)#钻孔终孔点坐标：
x1，i＝x1′+(2i-1)d1 cosβcosα
y1，i＝y1′+(2i-1)d1 cosβsinα
z1，i＝z1′+(2i-1)d1 sinβ
第2个钻场(j，i)#钻孔终孔点坐标：
xj，i＝x1′+(2i-1)d1 cosβcosα-(j-1)d2 cosγsinα
yj，i＝y1′+(2i-1)d1 cosβsinα+(j-1)d2 cosγcosα
zj，i＝z1′+(2i-1)d1 sinβ+(j-1)d2 sinγ
第m个钻场(j，i)#钻孔终孔点坐标：
xj，i＝x1′+(mi-1)d1 cosβcosα-(j-1)d2 cosγsinα
yj，i＝y1′+(mi-1)d1 cosβsinα+(j-1)d2 cosγcosα
zj，i＝z1′+(mi-1)d1 sinβ+(j-1)d2 sinγ
S13、各钻场开孔点坐标：
A2＝E1
假定第1个钻场的控制点坐标E1：(x0，y0，z0)
那么每个钻场的坐标可表示为：
第2个钻场控制点位置坐标：E2：(xe2，ye2，ze2)
xe2＝x0+id1 cosβcosα
ye2＝y0+id1 cosβsinα
ze2＝z0+id1 sinβ
第m个钻场控制点位置坐标Em：(xem，yem，zem)
xem＝x0+(m-1)id1 cosβcosα
yem＝y0+(m-1)id1 cosβsinα
zem＝z0+(m-1)id1 sinβ
第1个钻场(j，i)#钻孔开孔点坐标：
xj，i′＝x0-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
yj，i′＝y0+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
zj，i′＝z0+(j-1)d4 sinγ+(i-1)d3 sinβ
第2个钻场(j，1)#钻孔开孔点坐标：
xj，1′＝x0+id1 cosβcosα-(j-1)d4 cosγsinα
yj，1′＝y0+id1 cosβsinα+(j-1)d4 cosγcosα
zj，1′＝z0+id1 sinβ+(j-1)d4 sinγ
第2个钻场(j，i)#钻孔开孔点坐标：
xj，i′＝xe2-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
＝x0+id1 cosβcosα-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
＝x0+[id1+(i-1)d3]cosβcosα-(j-1)d4 cosγsinα
yj，i′＝ye2+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
＝y0+id1 cosβsinα+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
＝y0+[id1+(i-1)d3]cosβsinα+(j-1)d4 cosγcosα
zj，i′＝ze2+(j-1)d4 sinγ+(i-1)d3 sinβ
＝z0+id1 sinβ+(j-1)d4 sinγ+(i-1)d3 sinβ
＝z0+[id1+(i-1)d3]sinβ+(j-1)d4 sinγ
第m钻场(j，i)#钻孔开孔点坐标：
xj，i′＝xem-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
＝x0+(m-1)id1cosβcosα-(j-1)d4 cosγsinα+(i-1)d3 cosβcosα
＝x0+[(m-1)id1+(i-1)d3]cosβcosα-(j-1)d4 cosγsinα
yj，i′＝yem+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
＝y0+(m-1)id1 cosβsinα+(j-1)d4 cosγcosα+(i-1)d3 cosβsinα
＝y0+[(m-1)id1+(i-1)d3]cosβsinα+(j-1)d4 cosγcosα
zj，i′＝zem+(j-1)d4 sinγ+(i-1)d3 sinβ
＝z0+(m-1)id1 sinβ+(j-1)d4 sinγ+(i-1)d3 sinβ
＝z0+[(m-1)id1+(i-1)d3]sinβ+(j-1)d4 sinγ
各钻孔的钻孔长度：L孔长
各钻孔的钻孔方位角：α1
2.根据权利要求1所述的煤巷穿层瓦斯防治钻孔设计计算方法，其特征在于，所述计算方法还包括：
S2、穿层钻孔设计模型确定以及设计：
S21：穿层钻孔设计规则
(1)所有设计的钻孔孔长需要≤钻机的实际施工能力，以钻机最大的施工能力作为单条底板巷的控制范围；
(2)三条底板巷设计的钻孔，均按照最大的控制范围设计，当存在部分钻孔设计重叠，此时需要钻孔优化，相邻底板巷重叠的钻孔，按照钻孔长度进行比较，选重叠部分钻孔长度较短的进行施工；
S22：穿层钻孔计算
A：(x1，y1，z1)，控制点坐标A；
B：(x2，y2，z2)，控制点坐标B；
C：(x3，y3，z3)，控制点坐标C；
D：(x4，y4，z4)，控制点坐标D；
d：钻孔的开孔间距
Lz：预抽煤层工作面沿走向长：Lz＝LAB，单位为m
L面长：预抽煤层工作面长度：L面长＝LAD，单位为m
L孔长：钻孔长度，单位为m；
α：巷道方位角，单位为°；
β：煤层走向与水平面夹角，单位为°；
γ：煤层倾角，单位为°；
假定在理想条件下，底板巷和目标煤层的倾角，煤层沿走向夹角，以及方位角的都相同；
巷道方位角α：
煤层走向与水平面夹角β：
煤层倾角γ：
预抽区段煤层瓦斯区域防突措施的钻孔应当控制区段内的整个开采块段外侧一定范围内的煤层；
那么控制点坐标可重新计算确定：
A1：(x1′，y1′，z1′)
x1′＝x1-Lp cosβcosα+Lp cosγsinα
y1′＝y1-Lp cosβsinα-Lp cosγsinα
z1′＝z1-Lp sinβ-Lp sinγ
B1：(x2′，y2′，z2′)
x2′＝x2+Lp cosβcosα+Lp cosγsinα
y2′＝y2+Lp cosβsinα-Lp cosγsinα
z2′＝z2+Lp sinβ-Lp sinγ
C1：(x3′，y3′，z3′)
x3′＝x3+Lp cosβcosα-Lp cosγsinα
y3′＝y3+Lp cosβsinα+Lp cosγsinα
z3′＝z3+Lp sinβ+Lp sinγ
D1：(x4′，y4′，z4′)
x4′＝x4-Lp cosβcosα-Lp cosγsinα
y4′＝y4-Lp cosβsinα+Lp cosγsinα
z4′＝z4-Lp sinβ+Lp sinγ
加上预抽控制范围之后确定的走向长度：
加上预抽控制范围之后确定的煤层工作面长度：
各个点的三维坐标为：
A(x1，y1，z1)，B(x2，y2，z2)，C(x3，y3，z3)，D(x4，y4，z4)
A1(x1′，y1′，z1′)，B1(x2′，y2′，z2′)，C1(x3′，y3′，z3′)，D1(x4′，y4′，z4′)
假设三条底板巷的控制点坐标分别为(x0，y0，z0)，(x′0，y′0，z′0)，(x″0，y″0，z″0)
每条底板巷内设计施工m组钻孔
每组钻孔设计2n-1个钻孔，且：
且
式中：D为钻机打钻能力，m；
h为底板中心线开孔点到煤层的垂直距离，即开孔点对应的中心点终孔坐标的距离；
设计的钻孔长度最长不超过钻机可以施工的最大长度，假设钻机可以施工的最大长度为80m，则设计中的所有孔长均小于80m；
第一条底板巷内：
第m个开孔点坐标：
x0+(m-1)d cosβcosα
y0+(m-1)d cosβsinα
z0+(m-1)d sinβ
那么第m个开孔点对应的第m组终孔点坐标为：
中心点坐标
x0+(m-1)d cosβcosα
y0+(m-1)d cosβsinα
z0+(m-1)d sinβ+h
中心点向下第2个终孔点坐标：
x0+(m-1)d cosβcosα+d cosγcosα
y0+(m-1)d cosβsinα-d cosγsinα
z0+(m-1)d sinβ+h-d sinγ
中心点向下第n个终孔点坐标：
x0+(m-1)d cosβcosα+(n-1)d cosγcosα
y0+(m-1)d cosβsinα-(n-1)d cosγsinα
z0+(m-1)d sinβ+h-(n-1)d sinγ
中心点向上第2个终孔点坐标：
x0+(m-1)d cosβcosα-d cosγcosα
y0+(m-1)d cosβsinα+d cosγsinα
z0+(m-1)d sinβ+h+d sinγ
中心点向上第n个终孔点坐标：
x0+(m-1)dcosβcosα-(n-1)d cosγcosα
y0+(m-1)d cosβsinα+(n-1)d cosγsinα
z0+(m-1)d sinβ+h+(n-1)d sinγ
钻孔长度计算：
且L孔长小于钻机的打钻能力，取值为80m；
第二条底板巷内：
第m个开孔点坐标：
x′0+(m-1)d cosβcosα
y′0+(m-1)d cosβsinα
z′0+(m-1)d sinβ
那么第m个开孔点对应的第m组终孔点坐标为：
中心点坐标
x′0+(m-1)d cosβcosα
y′0+(m-1)d cosβsinα
z′0+(m-1)d sinβ+h
中心点向下第2个终孔点坐标：
x′0+(m-1)d cosβcosα+d cosγcosα
y′0+(m-1)d cosβsinα-d cosγsinα
z′0+(m-1)d sinβ+h-d sinγ
中心点向下第n个终孔点坐标：
x′0+(m-1)d cosβcosα+(n-1)d cosγcosα
y′0+(m-1)d cosβsinα-(n-1)d cosγsinα
z′0+(m-1)d sinβ+h-(n-1)d sinγ
中心点向上第2个终孔点坐标：
x′0+(m-1)d cosβcosα-2cosγcosα
y′0+(m-1)d cosβsinα+d cosγsinα
z′0+(m-1)d sinβ+h+d sinγ
中心点向上第n个终孔点坐标：
x′0+(m-1)d cosβcosα-(n-1)d cosγcosα
y′0+(m-1)d cosβsinα+(n-1)d cosγsinα
z′0+(m-1)d sinβ+h+(n-1)d sinγ
第三条底板巷内：
第m个开孔点坐标：
x″0+(m-1)d cosβcosα
y″0+(m-1)d cosβsinα
z″0+(m-1)d sinβ
那么第m个开孔点对应的第m组终孔点坐标为：
中心点坐标
x″0+(m-1)d cosβcosα
y″0+(m-1)d cosβsinα
z″0+(m-1)d sinβ+h
中心点向下第2个终孔点坐标：
x″0+(m-1)d cosβcosα+d cosγcosα
y″0+(m-1)d cosβsinα-d cosγsinα
z″0+(m-1)d sinβ+h-d sinγ
中心点向下第n个终孔点坐标：
x″0+(m-1)d cosβcosα+(n-1)d cosγcosα
y″0+(m-1)d cosβsinα-(n-1)d cosγsinα
z″0+(m-1)d sinβ+h-(n-1)d sinγ
中心点向上第2个终孔点坐标：
x″0+(m-1)d cosβcosα-d cosγcosα
y″0+(m-1)d cosβsinα+d cosγsinα
z″0+(m-1)d sinβ+h+d sinγ
中心点向上第n个终孔点坐标：
x″0+(m-1)d cosβcosα-(n-1)d cosγcosα
y″0+(m-1)d cosβsinα+(n-1)d cosγsinα
z″0+(m-1)d sinβ+h+(n-1)d sinγ
钻孔长度计算：
且L孔长小于钻机的打钻能力，取值为80m；
各钻孔的钻孔方位角α1：
3.根据权利要求1所述的煤巷穿层瓦斯防治钻孔设计计算方法，其特征在于，所述S11中穿层钻孔设计中的数据库包括空间坐标确定、抽采参数、煤层厚度、煤层倾角、瓦斯含量、瓦斯压力、施工钻场；规则库包括钻孔倾角、钻孔方位角、巷道方位角、钻孔长度计算、钻孔开孔位置、钻孔终孔位置、钻孔数量确定；输出参数包括钻孔倾角、钻孔方位角、巷道方位角、钻孔长度、钻孔开孔点坐标、钻孔终孔点坐标、钻孔数量。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110068860.9A CN112906148B (zh) | 2021-01-19 | 2021-01-19 | 一种煤巷穿层瓦斯防治钻孔设计计算方法 |
| ZA2021/02316A ZA202102316B (en) | 2021-01-19 | 2021-04-08 | Calculation method for design of gas-control crossing holes in coal roadway |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110068860.9A CN112906148B (zh) | 2021-01-19 | 2021-01-19 | 一种煤巷穿层瓦斯防治钻孔设计计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112906148A CN112906148A (zh) | 2021-06-04 |
| CN112906148B true CN112906148B (zh) | 2022-09-13 |
Family
ID=76033891
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110068860.9A Active CN112906148B (zh) | 2021-01-19 | 2021-01-19 | 一种煤巷穿层瓦斯防治钻孔设计计算方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112906148B (zh) |
| ZA (1) | ZA202102316B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107100663A (zh) * | 2017-05-02 | 2017-08-29 | 中国矿业大学 | 一种煤矿瓦斯的精准抽采方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE112006002137T5 (de) * | 2005-08-08 | 2008-06-26 | Halliburton Energy Services, Inc., Houston | Verfahren und Systeme zum Konstruieren und/oder Auswählen von Bohrausrüstung mit einer gewünschten Bohrmeißellenkbarkeit |
| CN102953679B (zh) * | 2011-08-22 | 2015-09-09 | 淮南矿业(集团)有限责任公司 | 钻孔设计方法 |
| CN102352769A (zh) * | 2011-10-21 | 2012-02-15 | 河南煤业化工集团研究院有限责任公司 | 高突矿井煤与瓦斯共采一体化开采方法 |
| CN103184887B (zh) * | 2013-03-08 | 2015-03-11 | 淮南矿业(集团)有限责任公司 | 井下巷道钻孔施工及地质信息反演方法 |
| CN107100666B (zh) * | 2017-06-15 | 2020-02-07 | 安徽理工大学 | 一种瓦斯突出煤层煤巷掘进的下爆上挤消突与降尘方法 |
| US11078786B2 (en) * | 2018-05-31 | 2021-08-03 | Saudi Arabian Oil Company | Salt mobility assessment and review technique (smart) for exploratory wells |
| CN110245427A (zh) * | 2019-06-17 | 2019-09-17 | 华夏天信（北京）智能低碳技术研究院有限公司 | 基于bim的煤矿工作面煤层gim数字模型创建方法 |
| CN110259378B (zh) * | 2019-07-03 | 2020-07-28 | 中煤科工集团重庆研究院有限公司 | 一种消除底抽巷穿层钻孔预抽煤层空白带的钻孔设计方法 |
-
2021
- 2021-01-19 CN CN202110068860.9A patent/CN112906148B/zh active Active
- 2021-04-08 ZA ZA2021/02316A patent/ZA202102316B/en unknown
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107100663A (zh) * | 2017-05-02 | 2017-08-29 | 中国矿业大学 | 一种煤矿瓦斯的精准抽采方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| ZA202102316B (en) | 2021-05-26 |
| CN112906148A (zh) | 2021-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114708393B (zh) | 全时空采掘过程的矿山应力场孪生建模同化系统及方法 | |
| CN101617101B (zh) | 井和泄流位置的自动化油田开发方案 | |
| CN106245852B (zh) | 异形混凝土拱柱钢筋定位方法 | |
| CN105626054A (zh) | 厚层断块油藏断棱精细刻画的方法 | |
| CN108442882A (zh) | 基于二维地震资料的页岩气大位移水平井随钻地质导向方法 | |
| CN113294137A (zh) | 水热型地热井工厂的建立方法及水热型地热井工厂 | |
| CN111411892B (zh) | 一种大井组丛式井轨迹防碰优化设计的方法 | |
| CN109635367A (zh) | 一种掘进机三维模拟方法、装置及系统 | |
| CN112906148B (zh) | 一种煤巷穿层瓦斯防治钻孔设计计算方法 | |
| CN106320978A (zh) | 基于井下定向孔的煤层等高线绘制方法 | |
| CN105551074B (zh) | 一种用于无人化采掘工作面的三维模型自动更新方法 | |
| CN110566177A (zh) | 一种工厂化钻井的布井方法 | |
| CN111090918B (zh) | 一种水平井井眼轨迹的设计方法及系统 | |
| CN108829994B (zh) | 多层页岩拟三维压裂模型缝高求解方法 | |
| CN210068123U (zh) | 页岩气藏多分支井开发井网结构 | |
| CN112963139B (zh) | 基于单筒多井的井口确定方法、装置、设备及存储介质 | |
| RU2696690C1 (ru) | Способ разработки многообъектного нефтяного месторождения (варианты) | |
| CN111612900B (zh) | 一种回采工作面煤层透明化三维地质模型构建方法及装置 | |
| RU2285795C1 (ru) | Способ разработки залежи нефти | |
| CN111485825B (zh) | 采煤工作面煤岩界面探测定向孔设计施工与数据处理方法 | |
| CN112819198B (zh) | 一种基于成本分析模型的工厂化钻井优化配置方法和系统 | |
| CN107859482B (zh) | 一种戟式煤层气水平井钻进方法 | |
| CN216741476U (zh) | 一种多层系低渗透油藏水平井、直井开发井网结构 | |
| CN116595700A (zh) | 一种基于三维地质模型的大平台水平井设计方法 | |
| CN113901632A (zh) | 一种提高页岩气藏三维地质模型构造精度的多级质控方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |