CN112906148A - 一种煤巷穿层瓦斯防治钻孔设计计算方法 - Google Patents

一种煤巷穿层瓦斯防治钻孔设计计算方法 Download PDF

Info

Publication number
CN112906148A
CN112906148A CN202110068860.9A CN202110068860A CN112906148A CN 112906148 A CN112906148 A CN 112906148A CN 202110068860 A CN202110068860 A CN 202110068860A CN 112906148 A CN112906148 A CN 112906148A
Authority
CN
China
Prior art keywords
drilling
cosβsinα
cosβcosα
sinβ
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.)
Granted
Application number
CN202110068860.9A
Other languages
English (en)
Other versions
CN112906148B (zh
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.)
Anhui University of Science and Technology
Original Assignee
Anhui University of Science and Technology
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 Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN202110068860.9A priority Critical patent/CN112906148B/zh
Priority to ZA2021/02316A priority patent/ZA202102316B/en
Publication of CN112906148A publication Critical patent/CN112906148A/zh
Application granted granted Critical
Publication of CN112906148B publication Critical patent/CN112906148B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Earth Drilling (AREA)
  • Road Paving Structures (AREA)

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
Figure BDA0002905253560000061
预抽煤巷条带宽度:L=LAD,m
Figure BDA0002905253560000062
每个钻场钻孔个数,个;
n=i*j
式中,i为每个钻场沿走向钻孔个数,个;
j为每个钻场沿倾向钻孔个数,个;
巷道方位角α:
Figure BDA0002905253560000063
煤层走向与水平面夹角β:
Figure BDA0002905253560000071
煤层倾角γ:
Figure BDA0002905253560000072
预抽区段煤层瓦斯区域防突措施的钻孔应当控制区段内的整个开采块段、两侧回采巷道及其外侧一定范围内的煤层;要求钻孔控制回采巷道外侧的范围是:倾斜、急倾斜煤层巷道上帮轮廓线外至少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γ
加上预抽控制范围之后确定的走向长度:
Figure BDA0002905253560000081
加上预抽控制范围之后确定的预抽煤巷条带宽度:
Figure BDA0002905253560000091
钻场个数:
Figure BDA0002905253560000092
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孔长
Figure BDA0002905253560000131
各钻孔的钻孔方位角:α1
Figure BDA0002905253560000132
Figure BDA0002905253560000133
Figure BDA0002905253560000134
本实施例通过依据瓦斯防治钻孔信息数据库、规则库以及输出参数,建立预抽煤巷穿层钻孔的设计模型,从而构建动态的三维瓦斯地质模型。
进一步的,所述计算方法还包括:
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
Figure BDA0002905253560000141
L面长:预抽煤层工作面长度:L面长=LAD,m
Figure BDA0002905253560000142
L孔长:钻孔长度,m;
α:巷道方位角(单位°);
β:煤层走向与水平面夹角(单位°);
γ::煤层倾角(单位°);
假定在理想条件下,底板巷和目标煤层的倾角,煤层沿走向夹角,以及方位角的都相同;
巷道方位角α:
Figure BDA0002905253560000151
煤层走向与水平面夹角β:
Figure BDA0002905253560000152
煤层倾角γ:
Figure BDA0002905253560000153
预抽区段煤层瓦斯区域防突措施的钻孔应当控制区段内的整个开采块段外侧一定范围内的煤层;
那么控制点坐标可重新计算确定:
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γ
加上预抽控制范围之后确定的走向长度:
Figure BDA0002905253560000161
加上预抽控制范围之后确定的煤层工作面长度:
Figure BDA0002905253560000162
参见图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组钻孔
Figure BDA0002905253560000171
每组钻孔设计2n-1个钻孔,且:
Figure BDA0002905253560000172
Figure BDA0002905253560000173
式中: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γ
钻孔长度计算:
Figure BDA0002905253560000181
且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γ
Figure BDA0002905253560000191
中心点向下第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γ
Figure BDA0002905253560000201
中心点向上第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γ
Figure BDA0002905253560000211
中心点向下第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γ
Figure BDA0002905253560000212
中心点向上第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γ
钻孔长度计算:
Figure BDA0002905253560000221
且L孔长小于钻机的打钻能力,一般为80m;
各钻孔的钻孔方位角α1
Figure BDA0002905253560000222
Figure BDA0002905253560000223
本实施例通过建立瓦斯防治无人化智能钻孔设计规则,从而确定穿层钻孔设计模型,实现了三维瓦斯地质模型的自动建立与精准纠正、抽采单元的自动划分与更新、钻孔智能设计与动态调整、钻孔合理性评判以及智能补救。
需要说明的是,我国煤矿穿层瓦斯钻孔施工普遍采用等间距的钻孔设计方法,其精度差、设计工作量大,不利于矿井的高效、安全生产。因此,研究煤矿瓦斯抽采钻孔的智能设计方法,具有十分重要的意义。本实施例具有以下有益效果:1)技术层面的预期效益:围绕瓦斯防治钻孔需求,进行瓦斯地质单元划分,提出钻孔优选判别原则,构建典型瓦斯抽采单元钻孔设计数据库,基于煤层瓦斯抽采单元划分结论及采掘时空布置,提出区域瓦斯防治钻孔无人化智能设计方法,设计面向钻孔机器人的瓦斯防治钻孔信息数据库;编制钻孔无人化智能设计与分析系统与钻孔机器人数据交互通信协议;开发钻孔无人化智能设计与分析系统。提高瓦斯防治钻孔设计的效率和科学性,提升企业经济效益,促进行业安全技术装备的升级换代,具有巨大的社会和经济效益。2)产业层面的预期指标及效益:在深刻认识主动扰动条件下瓦斯防治钻孔群周围多场耦合规律的基础上,形成瓦斯防治钻孔的智能设计方法,建成钻孔无人化智能设计系统软件平台。通过煤矿现场应用,实现对新技术、新装备的验证、改进,从而引领行业整体创新,推动行业技术全面进步,产生巨大的经济、社会效益。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种煤巷穿层瓦斯防治钻孔设计计算方法,其特征在于,所述计算方法包括:
S1、预抽煤巷穿层钻孔的设计模型确定以及设计;
S11、穿层钻孔设计:
依据数据库、规则库、输出参数进行设计;
假定在理想条件下,施工穿层钻孔的钻场所在的巷道,和目标煤层的倾角,煤层沿走向夹角,以及方位角的都相同;
预抽煤巷条带沿工作面走向长:Lz=LAB,m
Figure FDA0002905253550000011
预抽煤巷条带宽度:L=LAD,m
Figure FDA0002905253550000012
每个钻场钻孔个数,个;
n=i*j
式中,i为每个钻场沿走向钻孔个数,个;
j为每个钻场沿倾向钻孔个数,个;
巷道方位角α:
Figure FDA0002905253550000021
煤层走向与水平面夹角β:
Figure FDA0002905253550000022
煤层倾角γ:
Figure FDA0002905253550000023
预抽区段煤层瓦斯区域防突措施的钻孔应当控制区段内的整个开采块段、两侧回采巷道及其外侧一定范围内的煤层;要求钻孔控制回采巷道外侧的范围是:倾斜、急倾斜煤层巷道上帮轮廓线外至少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γ
加上预抽控制范围之后确定的走向长度:
Figure FDA0002905253550000041
加上预抽控制范围之后确定的预抽煤巷条带宽度:
Figure FDA0002905253550000042
钻场个数:
Figure FDA0002905253550000043
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)d1cosβ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)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)d3sinβ
=z0+(m-1)id1sinβ+(j-1)d4 sinγ+(i-1)d3 sinβ
=z0+[(m-1)id1+(i-1)d3]sinβ+(j-1)d4 sinγ#
各钻孔的钻孔长度:L孔长
Figure FDA0002905253550000081
各钻孔的钻孔方位角:α1
Figure FDA0002905253550000082
Figure FDA0002905253550000083
Figure FDA0002905253550000084
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
Figure FDA0002905253550000091
L面长:预抽煤层工作面长度:L面长=LAD,m
Figure FDA0002905253550000092
L孔长:钻孔长度,m;
α:巷道方位角(单位°);
β:煤层走向与水平面夹角(单位°);
γ::煤层倾角(单位°);
假定在理想条件下,底板巷和目标煤层的倾角,煤层沿走向夹角,以及方位角的都相同;
巷道方位角α:
Figure FDA0002905253550000101
煤层走向与水平面夹角β:
Figure FDA0002905253550000102
煤层倾角γ:
Figure FDA0002905253550000103
预抽区段煤层瓦斯区域防突措施的钻孔应当控制区段内的整个开采块段外侧一定范围内的煤层;
那么控制点坐标可重新计算确定:
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γ
加上预抽控制范围之后确定的走向长度:
Figure FDA0002905253550000111
加上预抽控制范围之后确定的煤层工作面长度:
Figure FDA0002905253550000112
各个点的三维坐标为:
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组钻孔
Figure FDA0002905253550000113
每组钻孔设计2n-1个钻孔,且:
Figure FDA0002905253550000121
Figure FDA0002905253550000122
式中: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+dsinγ
中心点向上第n个终孔点坐标:
x0+(m-1)d cosβcosα-(n-1)d cosγcosα
y0+(m-1)d cosβsinα+(m-1)d cosγsinα
z0+(m-1)d sinβ+h+(n-1)d sinγ
钻孔长度计算:
Figure FDA0002905253550000131
且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组终孔点坐标为:
Figure FDA0002905253550000141
Figure FDA0002905253550000151
第三条底板巷内:
第m个开孔点坐标:
x″0+(m-1)d cosβcosα
y″0+(m-1)d cosβsinα
z″0+(m-1)d sinβ
那么第m个开孔点对应的第m组终孔点坐标为:
Figure FDA0002905253550000152
Figure FDA0002905253550000161
钻孔长度计算:
Figure FDA0002905253550000162
且L孔长小于钻机的打钻能力,一般为80m;
各钻孔的钻孔方位角α1
Figure FDA0002905253550000163
Figure FDA0002905253550000164
3.根据权利要求1所述的煤巷穿层瓦斯防治钻孔设计计算方法,其特征在于,所述S11中穿层钻孔设计中的数据库包括空间坐标确定、抽采参数、煤层厚度、煤层倾角、瓦斯含量、瓦斯压力、施工钻场;规则库包括钻孔倾角、钻孔方位角、巷道方位角、钻孔长度计算、钻孔开孔位置、钻孔终孔位置、钻孔数量确定;输出参数包括钻孔倾角、钻孔方位角、巷道方位角、钻孔长度、钻孔开孔点坐标、钻孔终孔点坐标、钻孔数量。
CN202110068860.9A 2021-01-19 2021-01-19 一种煤巷穿层瓦斯防治钻孔设计计算方法 Active CN112906148B (zh)

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 true CN112906148A (zh) 2021-06-04
CN112906148B 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 (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070029111A1 (en) * 2005-08-08 2007-02-08 Shilin Chen Methods and systems for designing and/or selecting drilling equipment using predictions of rotary drill bit walk
CN102352769A (zh) * 2011-10-21 2012-02-15 河南煤业化工集团研究院有限责任公司 高突矿井煤与瓦斯共采一体化开采方法
CN102953679A (zh) * 2011-08-22 2013-03-06 淮南矿业(集团)有限责任公司 钻孔设计方法
CN103184887A (zh) * 2013-03-08 2013-07-03 淮南矿业(集团)有限责任公司 井下巷道钻孔施工及地质信息反演方法
CN107100663A (zh) * 2017-05-02 2017-08-29 中国矿业大学 一种煤矿瓦斯的精准抽采方法
CN107100666A (zh) * 2017-06-15 2017-08-29 安徽理工大学 一种瓦斯突出煤层煤巷掘进的下爆上挤消突与降尘方法
CN110245427A (zh) * 2019-06-17 2019-09-17 华夏天信(北京)智能低碳技术研究院有限公司 基于bim的煤矿工作面煤层gim数字模型创建方法
CN110259378A (zh) * 2019-07-03 2019-09-20 中煤科工集团重庆研究院有限公司 一种消除底抽巷穿层钻孔预抽煤层空白带的钻孔设计方法
US20190368348A1 (en) * 2018-05-31 2019-12-05 Saudi Arabian Oil Company Salt Mobility Assessment and Review Technique (SMART) for Exploratory Wells

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070029111A1 (en) * 2005-08-08 2007-02-08 Shilin Chen Methods and systems for designing and/or selecting drilling equipment using predictions of rotary drill bit walk
CN102953679A (zh) * 2011-08-22 2013-03-06 淮南矿业(集团)有限责任公司 钻孔设计方法
CN102352769A (zh) * 2011-10-21 2012-02-15 河南煤业化工集团研究院有限责任公司 高突矿井煤与瓦斯共采一体化开采方法
CN103184887A (zh) * 2013-03-08 2013-07-03 淮南矿业(集团)有限责任公司 井下巷道钻孔施工及地质信息反演方法
CN107100663A (zh) * 2017-05-02 2017-08-29 中国矿业大学 一种煤矿瓦斯的精准抽采方法
CN107100666A (zh) * 2017-06-15 2017-08-29 安徽理工大学 一种瓦斯突出煤层煤巷掘进的下爆上挤消突与降尘方法
US20190368348A1 (en) * 2018-05-31 2019-12-05 Saudi Arabian Oil Company Salt Mobility Assessment and Review Technique (SMART) for Exploratory Wells
CN110245427A (zh) * 2019-06-17 2019-09-17 华夏天信(北京)智能低碳技术研究院有限公司 基于bim的煤矿工作面煤层gim数字模型创建方法
CN110259378A (zh) * 2019-07-03 2019-09-20 中煤科工集团重庆研究院有限公司 一种消除底抽巷穿层钻孔预抽煤层空白带的钻孔设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WEIDONG ZHANG 等: "Research on Design and Construction of Directional Drilling of Shen-Wei Coal Pipeline", 《2016 INTERNATIONAL CONFERENCE ON SMART CITY AND SYSTEMS ENGINEERING (ICSCSE)》 *
罗亚涛: "顺层瓦斯抽采钻孔设计研究", 《内蒙古煤炭经济》 *

Also Published As

Publication number Publication date
ZA202102316B (en) 2021-05-26
CN112906148B (zh) 2022-09-13

Similar Documents

Publication Publication Date Title
CN114708393B (zh) 全时空采掘过程的矿山应力场孪生建模同化系统及方法
CN101617101B (zh) 井和泄流位置的自动化油田开发方案
US20120232873A1 (en) Systems and Methods for Improved Positioning of Pads
CN108442882A (zh) 基于二维地震资料的页岩气大位移水平井随钻地质导向方法
CN103362519B (zh) 煤巷安全快速掘进方法
CN103790577B (zh) 基于水平井水平段虚拟直井化的深度域约束反演方法
CN105626054A (zh) 厚层断块油藏断棱精细刻画的方法
RU2447270C1 (ru) Способ разработки нефтяной малоразведанной залежи
CN109630088A (zh) 一种高位钻孔布置位置的确定方法
CN113294137A (zh) 水热型地热井工厂的建立方法及水热型地热井工厂
CN111411892A (zh) 一种大井组丛式井轨迹防碰优化设计的方法
CN112906148B (zh) 一种煤巷穿层瓦斯防治钻孔设计计算方法
CN105551074B (zh) 一种用于无人化采掘工作面的三维模型自动更新方法
CN111090918B (zh) 一种水平井井眼轨迹的设计方法及系统
CN106320978A (zh) 基于井下定向孔的煤层等高线绘制方法
CN110566177A (zh) 一种工厂化钻井的布井方法
CN115077315B (zh) 一种基于凿岩台架结构的隧道炮孔分布方案获取方法
RU2696690C1 (ru) Способ разработки многообъектного нефтяного месторождения (варианты)
CN210068123U (zh) 页岩气藏多分支井开发井网结构
CN107060763A (zh) 一种实现极近距煤层经济、安全与高效开采的方法
RU2285795C1 (ru) Способ разработки залежи нефти
CN112963139A (zh) 基于单筒多井的井口确定方法、装置、设备及存储介质
CN112343656A (zh) 一种煤矿大数据下预测煤层底板标高应用方法与系统
CN111485825B (zh) 采煤工作面煤岩界面探测定向孔设计施工与数据处理方法
CN221856626U (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