CN113338929B - 一种单一窄工作面极不充分开采地表残余变形预测方法 - Google Patents
一种单一窄工作面极不充分开采地表残余变形预测方法 Download PDFInfo
- Publication number
- CN113338929B CN113338929B CN202110765425.1A CN202110765425A CN113338929B CN 113338929 B CN113338929 B CN 113338929B CN 202110765425 A CN202110765425 A CN 202110765425A CN 113338929 B CN113338929 B CN 113338929B
- Authority
- CN
- China
- Prior art keywords
- working face
- single narrow
- narrow working
- deformation
- coal
- 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
- 238000005065 mining Methods 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000003245 coal Substances 0.000 claims abstract description 96
- 230000010354 integration Effects 0.000 claims description 15
- 238000009825 accumulation Methods 0.000 claims description 4
- 208000029154 Narrow face Diseases 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims 1
- 239000011435 rock Substances 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000006578 abscission Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21C—MINING OR QUARRYING
- E21C39/00—Devices for testing in situ the hardness or other properties of minerals, e.g. for giving information as to the selection of suitable mining tools
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21C—MINING OR QUARRYING
- E21C41/00—Methods of underground or surface mining; Layouts therefor
- E21C41/26—Methods of surface mining; Layouts therefor
- E21C41/28—Methods of surface mining; Layouts therefor for brown or hard coal
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种单一窄工作面极不充分开采地表残余变形预测方法,适用于煤矿采空区场地稳定性评价使用。首先计算出单一窄工作面采空区两侧煤壁最大屈服和剥离宽度的总和;然后预测两侧煤壁屈服剥离前单一窄工作面采空区地表的移动数与变形值;接着计算单一窄工作面采空区宽度加上采空区两侧煤壁的最大屈服和剥离宽度之和,再预测两侧煤壁屈服剥离后的该采空区地表的移动数与变形值,最后通过对该单一工作面采空区两侧煤壁屈服数值和剥离前后预测的地表的移动与变形值作差,得到该单一窄工作面极不充分开采地表残余变形量,实现单一窄工作面极不充分开采残余变形的预测。其步骤简单,预测精度高,为煤矿采空区场地稳定性评价提供了理论依据。
Description
技术领域
本发明涉及一种地表残余变形预测方法,尤其适用于煤矿采空区场地稳定性评价使用的一种单一窄工作面极不充分开采地表残余变形预测方法。
背景技术
我国长期大规模煤炭开采在保证国家能源安全的同时,形成了大面积采空塌陷区。虽然经过长时间的自然压实,煤矿采空区及其上方破裂岩体处于一个相对平衡状态,但采空区覆岩内的裂隙、离层和空洞依然存在,在各种内外荷载(建筑荷载、振动荷载、地震力、地下水等)的作用下,上覆岩层的这种相对平衡状态将再次被打破,从而形成地表的二次沉降,这一过程称为老采空区的“活化”。矿区大面积采空塌陷区不良地基场地,对矿业城市转型和现代化建设迫切需要的建设用地构成了严重制约和安全隐患。
采空区场地稳定性评价是进行场地建设利用的依据,其中采空区场地地表残余变形预测是进行采空区场地稳定性评价的关键。目前国内外学者针对长壁式采空区场地地表残余变形预测方法做了大量的研究工作,建立了面向不同地质采矿条件的长壁式采空区场地地表残余变形预测方法,且已编入相应的国家规范中,已基本解决了长壁式采空区场地地表残余变形预测问题。但单一窄工作面采空区覆岩中能够形成相应比较稳定的自然平衡拱式次生结构体系,与长壁式采空区覆岩结构有显著差异,进而导致单一窄工作采空区和长壁式采空区场地地表残余变形量有较大差异。但目前尚未有计算单一窄工作面极不充分开采地表残余变形预测方法,导致单一窄工作面采空区场地稳定性评价缺乏科学依据。
发明内容
针对现有技术的不足之处,提供一种步骤简单,精准度高,兼顾了采空区两侧煤壁屈服和剥离影响的单一窄工作面极不充分开采地表残余变形预测方法。
为实现上述技术目的,本发明的一种单一窄工作面极不充分开采地表残余变形预测方法,首先计算出单一窄工作面采空区两侧煤壁最大屈服和剥离宽度的总和;然后利用概率积分法预测两侧煤壁屈服剥离前单一窄工作面采空区地表的移动数与变形值;接着计算单一窄工作面采空区宽度加上采空区两侧煤壁的最大屈服和剥离宽度之和,再通过概率积分法预测两侧煤壁屈服剥离后的该采空区地表的移动数与变形值,最后通过对该单一工作面采空区两侧煤壁屈服数值和剥离前后预测的地表的移动与变形值作差,得到该单一窄工作面极不充分开采地表残余变形量,实现单一窄工作面极不充分开采残余变形预测。
具体步骤如下:
利用矿区基础资料确定研究区域的地层情况,包括研究区域煤层的采深、厚度等信息;
根据煤层的采深、采厚数值计算出单一窄工作面采空区两侧煤壁最大屈服宽度;
根据煤层的采厚、碎胀系数、堆积体休止角计算出单一窄工作面采空区两侧煤壁最大剥离宽度;
之后根据单一窄工作面采空区宽度值,基于概率积分法预测单一窄工作面采空区的移动与变形值;
接着将单一窄工作面采空区宽度加上采空区两侧煤壁的最大屈服和剥离宽度之和,再基于概率积分法预测采空区和两侧煤壁上方地表移动与变形值;
最后通过对单一工作面采空区两侧煤壁屈服和剥离前后预测的地表移动与变形值作差,从而预测出单一窄工作面极不充分开采地表残余变形量,实现单一窄工作面极不充分开采残余变形预测。
进一步,根据所述地层情况包括研究区域地质采矿条件及煤层赋存条件,确定煤层厚度m和煤层深度H,之后利用公式Y=0.0049mH计算出该区域单一窄工作面采空区两侧煤壁最大屈服宽度Y,式中m为煤层厚度,H为煤层深度。
进一步,单一窄工作面采空区两侧煤壁剥离的最大宽度的计算步骤为:
进一步,利用公式:bm=Y+dm计算单一窄工作面采空区两侧煤壁最大屈服和剥离宽度和,其中Y为单一窄工作面采空区两侧煤壁最大屈服宽度,dm为最大剥离宽度。
计算任意点A(x,y)的下沉值W(x,y)具体为
W(x,y)=WcmC'xC'y (1)
其中:
式中:Wcm为充分采动条件下地表最大下沉值,Wcm=mqcosα;m为采出煤层厚度;α为煤层倾角;q为地表下沉系数;Cx′,Cy′为待求点在走向和倾向主断面上投影点处的下沉分布系数;l,b为采区拐点平移后走向长度及倾斜方向在地表的计算开采宽度;r,r1,r2分别为走向、下山、上山方向的主要影响半经;x、y为待求点坐标;
式中:
K(x,y)max,K(x,y)min分别为待求点最大、最小曲率变形值;Kx,Ky分别为待求点沿走向及倾向在主断面投影处迭加后的曲率值;
式中:ε(x,y)max,ε(x,y)min为待求点最大、最小水平变形值。
进一步,利用公式:b′=b+bm计算采空区两侧煤壁屈服和剥离后的采空区实际宽度,式中b为单一窄工作面采空区宽度,bm为两侧煤壁最大屈服和剥离宽度和,同样基于概率积分法可以预测得到工作面采空区两侧煤壁屈服和剥离后的地表移动与变形值:下沉值W(x,y)′、倾斜变形值曲率变形水平移动值水平变形值
有益效果:本申请能兼顾单一窄工作面采空区两侧煤壁的屈服和剥离特征,推导出单一窄工作面的实际采空区宽度,从而计算出单一窄工作面采空区场地地表残余量,实现了单一窄工作面采空区场地地表残余变形预测,为单一窄工作面采空区场地稳定性评价提供了科学依据。
附图说明
图1为本发明单一窄工作面极不充分开采残余变形预测方法流程图
图2为本发明单一窄工作面两侧煤壁应力分区示意图
图3为本发明单一窄工作面采空区两侧煤壁剥离模型示意图
图4为本发明单一窄工作面采空区场地地表残余下沉分布等值线图
具体实施方式
下面将结合图和具体实施过程对本发明做进一步详细说明:
如图1所示,本发明的一种单一窄工作面极不充分开采地表残余变形预测方法,首先计算出单一窄工作面采空区两侧煤壁最大屈服和剥离宽度的总和;然后利用概率积分法预测两侧煤壁屈服剥离前单一窄工作面采空区地表的移动数与变形值;接着计算单一窄工作面采空区宽度加上采空区两侧煤壁的最大屈服和剥离宽度之和,再通过概率积分法预测两侧煤壁屈服剥离后的该采空区地表的移动数与变形值,最后通过对该单一工作面采空区两侧煤壁屈服数值和剥离前后预测的地表的移动与变形值作差,得到该单一窄工作面极不充分开采地表残余变形量,实现单一窄工作面极不充分开采残余变形预测。
具体步骤如下:
利用矿区基础资料确定研究区域的地层情况,包括研究区域煤层的采深、厚度等信息;
根据煤层的采深、采厚数值计算出单一窄工作面采空区两侧煤壁最大屈服宽度;
如图2、图3所示,首先利用矿区的地质采矿资料,获得煤层的采深H、厚度m,利用公式Y=0.0049mH,则可以计算该区域单一窄工作面采空区两侧煤壁最大屈服宽度Y;
(1)任意点A(x,y)的下沉值W(x,y)
W(x,y)=WcmC'xC'y (15)
其中:
式中:
Wcm—充分采动条件下地表最大下沉值,Wcm=mqcosα;
m—采出煤层厚度;α—煤层倾角;
q—地表下沉系数;
Cx′,Cy′—为待求点在走向和倾向主断面上投影点处的下沉分布系数;
l,b—为采区拐点平移后走向长度及倾斜方向在地表的计算开采宽度;
r,r1,r2——分别为走向、下山、上山方向的主要影响半经;
x、y——待求点坐标。
式中:
T(x,y)m——待求点的最大倾斜值,mm/m;
Tx,Ty——分别为待求点沿走向和倾向主断面上投影点处迭加后的倾斜变形值,mm/m。
式中:
K(x,y)max,K(x,y)min——分别为待求点最大、最小曲率变形值;
Kx,Ky——分别为待求点沿走向及倾向在主断面投影处迭加后的曲率值。
Ux,Uy——分别为待求点沿走向和倾向在主断面投影点处的水平移动值,mm。对于倾斜方向需加Cy′·Wcm·ctgθ。
ε(x,y)max,ε(x,y)min——为待求点最大、最小水平变形值;
然后单一窄工作面采空区宽度b和两侧煤壁最大屈服和剥离宽度和bm,可以计算得到采空区两侧煤壁屈服和剥离后的采空区实际宽度b′=b+bm,同样基于概率积分法可以预测得到工作面采空区两侧煤壁屈服和剥离后的地表移动与变形值
最后通过对单一工作面采空区两侧煤壁屈服和剥离前后预测的地表移动与变形值作差,可以得到单一窄工作面极不充分开采地表残余变形量,实现单一窄工作面极不充分开采残余变形预测。
实施例一、
步骤1:收集某研究区域地质采矿条件资料,确定煤层采厚为2.47m、采深为466.6m。
步骤2:利用公式Y=0.0049mH,则可以计算该区域单一窄工作面采空区两侧煤壁最大屈服宽度Y=5.647m。
步骤5:根据单一窄工作面采空区宽度50m和两侧煤壁最大屈服和剥离宽度和6.792m,可以计算得到采空区两侧煤壁屈服和剥离后的采空区实际宽度56.792m。
步骤7:通过对单一工作面采空区两侧煤壁屈服和剥离前后预测的地表移动与变形值作差,可以得到单一窄工作面极不充分开采地表残余变形量。生成如图4所示该实例单一窄工作面采空区场地地表残余下沉分布等值线图。
Claims (7)
1.一种单一窄工作面极不充分开采地表残余变形预测方法,其特征在于:首先计算出单一窄工作面采空区两侧煤壁最大屈服和剥离宽度的总和;然后利用概率积分法预测两侧煤壁屈服剥离前单一窄工作面采空区地表的移动数与变形值;接着计算单一窄工作面采空区宽度加上采空区两侧煤壁的最大屈服和剥离宽度之和,再通过概率积分法预测两侧煤壁屈服剥离后的该采空区地表的移动数与变形值,最后通过对单一窄工作面采空区两侧煤壁屈服数值和剥离前后预测的地表的移动与变形值作差,得到该单一窄工作面极不充分开采地表残余变形量,实现单一窄工作面极不充分开采残余变形预测。
2.根据权利要求1所述的单一窄工作面极不充分开采地表残余变形预测方法,其特征在于具体步骤如下:
利用矿区基础资料确定研究区域的地层情况,包括研究区域煤层的采深、厚度信息;
根据煤层的采深、采厚数值计算出单一窄工作面采空区两侧煤壁最大屈服宽度;
根据煤层的采厚、碎胀系数、堆积体休止角计算出单一窄工作面采空区两侧煤壁最大剥离宽度;
之后根据单一窄工作面采空区宽度值,基于概率积分法预测单一窄工作面采空区的移动与变形值;
接着将单一窄工作面采空区宽度加上采空区两侧煤壁的最大屈服和剥离宽度之和,再基于概率积分法预测采空区和两侧煤壁上方地表移动与变形值;
最后通过对单一窄工作面采空区两侧煤壁屈服和剥离前后预测的地表移动与变形值作差,从而预测出单一窄工作面极不充分开采地表残余变形量,实现单一窄工作面极不充分开采残余变形预测。
3.根据权利要求2所述的单一窄工作面极不充分开采地表残余变形预测方法,其特征在于:根据所述地层情况包括研究区域地质采矿条件及煤层赋存条件,确定煤层厚度m和煤层深度H,之后利用公式Y=0.0049mH计算出该区域单一窄工作面采空区两侧煤壁最大屈服宽度Y,式中m为煤层厚度,H为煤层深度。
5.根据权利要求4所述的单一窄工作面极不充分开采地表残余变形预测方法,其特征在于:利用公式:bm=Y+dm计算单一窄工作面采空区两侧煤壁最大屈服和剥离宽度和,其中Y为单一窄工作面采空区两侧煤壁最大屈服宽度,dm为最大剥离宽度。
6.根据权利要求4所述的单一窄工作面极不充分开采地表残余变形预测方法,其特征在于根据单一窄工作面采空区宽度b,基于概率积分法预测地表移动与变形值:下沉值W(x,y)、倾斜变形值曲率变形水平移动值水平变形值具体预测模型如下:
计算任意点A(x,y)的下沉值W(x,y)具体为
W(x,y)=WcmC'xC'y (1)
其中:
式中:Wcm为充分采动条件下地表最大下沉值,Wcm=mqcosα;m为采出煤层厚度;α为煤层倾角;q为地表下沉系数;Cx′,Cy′为待求点在走向和倾向主断面上投影点处的下沉分布系数;l,b为采区拐点平移后走向长度及倾斜方向在地表的计算开采宽度;r,r1,r2分别为走向、下山、上山方向的主要影响半经;x、y为待求点坐标;
式中:
K(x,y)max,K(x,y)min分别为待求点最大、最小曲率变形值;Kx,Ky分别为待求点沿走向及倾向在主断面投影处迭加后的曲率值;
式中:ε(x,y)max,ε(x,y)min为待求点最大、最小水平变形值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110765425.1A CN113338929B (zh) | 2021-07-07 | 2021-07-07 | 一种单一窄工作面极不充分开采地表残余变形预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110765425.1A CN113338929B (zh) | 2021-07-07 | 2021-07-07 | 一种单一窄工作面极不充分开采地表残余变形预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113338929A CN113338929A (zh) | 2021-09-03 |
CN113338929B true CN113338929B (zh) | 2022-08-30 |
Family
ID=77482955
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110765425.1A Active CN113338929B (zh) | 2021-07-07 | 2021-07-07 | 一种单一窄工作面极不充分开采地表残余变形预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113338929B (zh) |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2162149C1 (ru) * | 1999-06-22 | 2001-01-20 | Институт горного дела СО РАН | Способ определения устойчивости пород при креплении горных выработок |
CN105926569B (zh) * | 2016-04-28 | 2018-02-16 | 河北地质大学 | 一种基于沉降监测数据的煤矿老采空区场地稳定性定量评价方法 |
CN107092740B (zh) * | 2017-04-12 | 2020-08-04 | 乌鲁木齐城市轨道集团有限公司 | 一种急倾斜特厚煤层老采空区残余变形预计方法 |
CN108021763A (zh) * | 2017-12-15 | 2018-05-11 | 济南大学 | 一种老采空区残余沉降预测方法 |
CN107989612B (zh) * | 2018-01-08 | 2019-11-15 | 西安科技大学 | 一种高潜水位煤矿区动态预复垦新方法 |
CN110750866B (zh) * | 2019-09-17 | 2022-04-26 | 中国矿业大学 | 一种利用无人机技术快速获得矿区开采沉陷预计参数的方法 |
CN111400659A (zh) * | 2019-12-03 | 2020-07-10 | 中国矿业大学 | 采-选-充一体化开采技术耕地保护率计算方法 |
AU2020102336A4 (en) * | 2020-09-18 | 2020-10-29 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method, system and device for identifying water-filled goaf area in coal mine |
CN113077088A (zh) * | 2021-04-06 | 2021-07-06 | 安徽大学 | 城中采煤沉陷区水土资源空间格局重构方法 |
-
2021
- 2021-07-07 CN CN202110765425.1A patent/CN113338929B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113338929A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111750822B (zh) | 一种采煤诱发的覆岩与地表沉陷协同动态预测方法 | |
Zhang et al. | Study on overlying strata movement patterns and mechanisms in super-large mining height stopes | |
CN107330194B (zh) | 盾构隧道开挖对邻近群桩基础影响的数值模拟分析方法 | |
CN112446143B (zh) | 一种厚松散层薄基岩下条带充填开采参数设计方法 | |
CN106150546B (zh) | 一种根据地表沉陷特征进行采空区部分充填的方法 | |
WO2023000837A1 (zh) | 深埋采场覆岩类型评价标准及薄基岩加厚改造设计方法 | |
CN115030722A (zh) | 一种采空区滞后充填高效保水采煤方法 | |
CN105971630A (zh) | 一种近距离煤层巷道顶板冒落防治方法 | |
CN110046470A (zh) | 一种盾构掘进引起工后地表沉降量的确定方法 | |
Zhou et al. | Stress-release law and deformation characteristics of large-span tunnel excavated with semi central diaphragm method | |
CN106988744A (zh) | 一种巷旁“支护体‑充填体”协调支护沿空留巷方法 | |
CN117167080A (zh) | 一种三软厚煤层开采覆岩多层位离层注浆减沉方法 | |
Lei et al. | Study on the geological adaptability of the arch cover method for shallow-buried large-span metro stations | |
Wang et al. | Improving stress environment in development entries through an alternate longwall mining layout | |
Li et al. | Mechanism of overlying strata migration and failure during underground mining in the mountainous carbonate areas in southwestern China | |
CN113944510B (zh) | 一种岩层移动共轭内、外“类双曲线”模型判定方法 | |
CN114017103A (zh) | 一种长壁逐巷胶结充填采煤工作面参数设计方法 | |
CN113338929B (zh) | 一种单一窄工作面极不充分开采地表残余变形预测方法 | |
Xue et al. | Numerical investigation on overburden migration behaviors in stope under thick magmatic rocks | |
CN114075983A (zh) | 一种厚松散层下采动松散层结构的分类方法 | |
CN115344928B (zh) | 一种耕地下采煤区域变形控制的地表沉陷叠加计算方法 | |
CN110630316A (zh) | 一种废弃工作面导气裂隙带发育特征判别方法 | |
CN115829133A (zh) | 一种充填工艺回收条带煤柱的地表沉陷预测方法 | |
CN116502390A (zh) | 一种基于flac3d计算矿区应力分布和位移变化的方法 | |
Pan et al. | Construction mechanics behavior of double-layer primary support arch cover method |
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 |