CN111101941B - 一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法 - Google Patents
一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法 Download PDFInfo
- Publication number
- CN111101941B CN111101941B CN202010012470.5A CN202010012470A CN111101941B CN 111101941 B CN111101941 B CN 111101941B CN 202010012470 A CN202010012470 A CN 202010012470A CN 111101941 B CN111101941 B CN 111101941B
- Authority
- CN
- China
- Prior art keywords
- subsidence
- shallow
- pad
- area
- elevation
- 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
- 238000009412 basement excavation Methods 0.000 title claims abstract description 138
- 238000005065 mining Methods 0.000 title claims abstract description 136
- 239000003245 coal Substances 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000005457 optimization Methods 0.000 title claims abstract description 12
- 238000010276 construction Methods 0.000 claims abstract description 42
- 239000002689 soil Substances 0.000 claims abstract description 31
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 18
- 230000010354 integration Effects 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 230000000203 accumulation affect Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000009189 diving Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000035558 fertility Effects 0.000 abstract description 3
- 230000002093 peripheral effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 239000005442 atmospheric precipitation Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000994 depressogenic effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21C—MINING OR QUARRYING
- E21C41/00—Methods of underground or surface mining; Layouts therefor
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法,适用于高潜水位矿区沉陷积水区地质环境治理。首先收集矿区地质参数,设计地下采煤工作面走向主断面正上方各地表点挖深垫浅后高程;根据地下采煤工作面走向主断面正上方各地表点的下沉值与深垫浅后高程值,求出走向主断面上挖深垫浅的边界点及沉陷盆地的挖深垫浅边界线;设计沉陷盆地内挖深垫浅区域与垫浅区域的挖填土高程与治理后的最终高程,并完成沉陷区的填土施工。其工艺简单,成本低廉,有效缩小了沉陷面积,无需从周边土地二次取土,既能保障覆土的肥力,又降低了二次污染的风险。
Description
技术领域
本发明主要涉及一种超前治理方法,尤其适用于一种高潜水位矿区使用的基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法。
背景技术
采煤沉陷积水区是东部高潜水位矿区开采沉陷造成的一种典型的地质灾害。东部矿区普遍具有储量大、多煤层、高潜水位、地表叠加沉降等特点,开采造成地表大量土地资源丧失,原有耕地破坏,地表形成大面积塌陷坑,同时受到地下潜水外渗、大气降水的综合影响,塌陷坑内常年积水,历经几十年重复开采,大片开采沉陷区形成地表塌陷塘,地表原有耕地、道路、村庄等被淹没,塌陷区原有生态系统发生不可逆的破坏。为保证国民经济建设,大量煤炭资源被采出的同时,采煤沉陷造成大量耕地、村庄丧失,以两淮矿区为例,截至2015年底,已形成的塌陷区面积达508km2,且以20km2/年的速度在增加,同时受到高潜水位和地表降水的综合影响,塌陷区形成大量水域,原有陆生生态系统发生颠覆性变化,转变为水生生态系统和水陆复合生态系统,动态沉陷形成的生态系统服务功能低下,导致生态环境恶化,大量沉陷区土地撂荒,社会问题突出,严重影响了沉陷区可持续发展。
常规沉陷区挖深垫浅一般是等采煤沉陷区完全稳定后,再进行改造,且仅凭工程需要对挖深区域与填土区域进行区分,该方法的缺点在于采煤沉陷区稳沉时间较长,待环境治理完全之后,大量沉陷区土地长期撂荒,未得到有效利用;粗暴地挖填未设计合理边界用于水土保持,且挖深区域与垫浅区域划分不明确,需要从外界取土,对环境造成二次扰动。
发明内容
本发明的目的是提供一种工艺简单、成本低廉、在沉陷区非稳沉条件下进行治理的基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法。
为实现上述技术目的,本发明的基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法,首先根据矿区地质采矿条件参数设计高潜水位采沉陷区挖深垫浅后地表高程;计算采煤工作面走向主断面上挖深垫浅的边界点及采空后沉陷盆地的挖深垫浅边界线;预计采煤工作面采空后形成沉陷区内进行治理需要的填土高程;利用采煤沉陷区内的填土高程,在地下矿区工作面开采前对地表对应区域进行超前治理;当地下开采完成后地表自动形成预计的治理效果,平整地面,从而完成沉陷区挖深垫浅的治理方法;
具体步骤如下
a..根据开采工作面所在矿区地质采矿条件参数,包括地表地形信息和地下开采计划,预测地表沉陷信息,利用地表沉陷信息结合原始地表地形信息,获取地表沉陷等值线图,设计出沉陷区挖深垫浅后地表高程值计算公式:
利用概率积分法地表沉陷预计,计算开采工作面走向主断面上方地表点下沉值W(x);选取开采工作面走向主断面进行挖深垫浅边界设计的参照面,选取走向主断面地表最大下沉值Wmax为最小设计高程,以走向主断面土方量平衡为条件,计算出挖深垫浅后地表高程值中变量b与k的关系式,根据实际工程施工、采煤沉陷地表变形规律、变量b与k的关系式分析设计挖深垫浅后地表高程值公式中的变量k,进而得出挖深垫浅后地表高程值Z(x)=kx+b;
下式为采煤工作面走向主断面正上方地表各点下沉值计算公式W(x),工作面走向主断面正上方地表最大下沉值Wmax:
式中:x为采煤工作面走向主断面正上方地表点相对于工作面的坐标;W0为概率积分法采煤沉陷预计地表最大下沉值;l为地下采煤工作面走向计算长度;L为地下采煤工作面倾向计算长度;D1为地下采煤工作面走向实际长度;Cym为倾向采动程度系数;Cy为倾向采动系数;s为工作面拐点偏移距;r为工作面主要影响半径;erf()为高斯误差函数;
b.为了减小积水区域,最大程度的改造土地,需要计算走向主断面上挖深垫浅的各边界点,获得沉陷盆地的挖深垫浅的地表内边界信息、挖深垫浅边界信息、挖深盆地边界信息:
首先使用二分法找出开采工作面走向主断面正上方地表下沉值W(x)与挖深垫浅后地表高程值Z(x)相等的点i坐标即挖深垫浅的边界点;挖深垫浅后地表高程值Z(x)与地面原始平均高程H0相等的地表内边界点p即为地表内边界点,利用挖深垫浅后地表高程值公式计算出地表内边界点p横坐标为挖深垫浅后地表高程值Z(x)与地表最大下沉值Wmax相等的点f即为挖深盆地边界点,根据挖深垫浅后地表高程值公式计算出挖深盆地边界点f横坐标为点i′为点i在地表的投影,点f′为点f在地表的投影,则地表处点p,点i′,点f′由采煤工作面走向主断面正上方地表下沉计算公式W(x)可得出上述三点处下沉值W(p),W(i′),W(f′),地表点p,点i′,点f′处下沉等值线即为沉陷盆地的挖深垫浅的地表内边界、挖深垫浅边界、挖深盆地边界;
c.在开采工作面开采前或沉陷区稳沉前,根据填土高程设计进行施工,填土高程施工设计分为填土设计高程H1与治理后的最终高程H2,在工作面走向主断面与地表相交处的每个地表点上的填土设计高程H1就等于地表原始高程H0加上挖深垫浅边界高程值Z(x)减去该点处地表移动下沉值W(x),值为正时为填方高程,值为负时为挖方高程;在整个沉陷盆地内,根据沉陷预计等值线图与走向主断面的挖填方高程的设计,下沉值相同处的高程设计相同;治理后的最终高程H2等于填土设计高程H1减去地表移动下沉值W(x),完成沉陷区的挖深垫浅施工设计;
d.利用沉陷区的挖深垫浅施工设计,对整个地表沉陷盆地的挖深区域和垫浅区域进行改造,根据填土高程设计施工,地表点i′处的沉陷等值线即为挖深垫浅边界,边界内部沉陷盆地为挖深区域,边界外侧沉陷边缘为垫浅区域;在挖深区域,地表点f′处的沉陷等值线即为挖深盆地边界,挖深盆地边界内部开挖区填土设计高程H1=H0+W(x)-Wmax,平整开发区底部地面;在挖深盆地边界与挖深垫浅边界之间的环形区域,开挖区填土设计高程H1=H0+W(x)-Z(x),将挖掘出的土方填入垫浅区;在垫浅区域,地表点p处的沉陷等值线为地表内边界,在地表内边界与挖深垫浅边界之间的环形区域,垫浅区填土设计高程H1=H0+W(x)-Z(x),在沉陷边界与地表内边界之间的环形区域,垫浅区填土设计高程H1=H0+W(x),平整垫浅区域地面,从而完成对治理矿区采煤工作面上方沉陷区的地质区域改造。
挖深垫浅后地表高程值Z(x)计算公式中变量b与k关系式的求取:
根据土方量平衡条件:
式中:为挖深垫浅后地表高程值Z(x)与地面平均高程相等的挖深盆地边界点p处的横坐标,点x0地下采煤工作面走向主断面正上方地表下沉值W(x)与挖深垫浅后地表高程值Z(x)相等的点i处的横坐标,为挖深垫浅后地表高程值Z(x)与地表最大下沉值Wmax相等的点f横坐标。
挖深垫浅后地表高程值计算公式中b与k的关系式:
挖深垫浅后地表高程值Z(x)计算公式中变量k的选取:
k为挖深垫浅后地表高程值计算公式中的变量未知,结合实际工程施工、采煤沉陷地表变形规律、变量b与k的关系式可知:在走向主断面正上方的地表点,地表下沉值与挖深垫浅后地表高程值相等时的点i有且仅有一个点i;在点i左侧W(x)>Z(x),在点i右侧W(x)<Z(x),结合概率积分法主断面地表移动与变形规律分析,k必须大于地表下沉曲线W(x)的最大斜率,即大于x方向的地表最大倾斜值,即根据挖深垫浅后地表高程值计算公式中变量b与k的关系式分析,截距b随挖深垫浅后地表高程曲线的斜率k增大而减小,当k偏大时,此时两曲线交点偏上,根据实际验证,明显不符合挖填方平衡原则。综上分析,挖深垫浅后地表高程曲线斜率k的取值需略大于地表下沉曲线最大斜率;
由此得出k的取值,进而由b与k的关系式得出b的值,最后求出挖深垫浅挖深垫浅后地表高程曲线Z(x)=kx+b;同时挖深垫浅后地表高程值Z(x)计算公式中变量k的选取,使得采煤沉陷区内积水区与地面之间存在适宜的坡度,更有利于沉陷区水土保持与农业利用。
在工作面开采前或沉陷区稳沉前进行沉陷区的挖深垫浅填土设计高程与施工;通常为了使土体不被闲置,沉陷区稳沉后进行施工挖深垫浅,但当高潜水位矿区在地表沉陷时沉陷区产生积水时,需要在工作面开采前或沉陷区稳沉前进行沉陷区的挖深垫浅填土设计高程与施工,避免了沉陷期土地闲置及积水影响施工的问题;常规的边采边复高程模型中高程受开采进度影响,在工作面开采前或沉陷区稳沉前根据开采沉陷设计填土高程,治理后的地表高程不受开采工作的影响。
有益效果:本发明提供了基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法,适用于我国东部高潜水位矿区沉陷积水区地质环境治理。有效缩小了沉陷面积,无需从周边土地二次取土,既能保障覆土的肥力,又降低了二次污染的风险;在沉陷区稳沉前治理,避免了沉陷期土地闲置的问题,同时挖深垫浅后高程的设计,明确了挖深区域与垫浅区域的边界,并使得积水区与地面之间存在适宜的坡度,更有利于沉陷区水土保持与农业利用,大大提高了沉陷区土地的利用率。
在沉陷区非稳沉状态下进行超前治理的方法,解决了非稳沉条件下沉陷区域挖深垫浅土地复垦的问题,最大程度地改造矿区生态环境,实现了沉陷区水土综合利用。
本发明工艺简单,成本低廉,有效缩小了沉陷面积,无需从周边土地二次取土,既能保障覆土的肥力,又降低了二次污染的风险;在沉陷区稳沉前治理,避免了沉陷期土地闲置的问题,同时挖深垫浅后高程的设计,明确了挖深区域与垫浅区域的边界,并使得积水区与地面之间存在适宜的坡度,更有利于沉陷区水土保持与农业利用,大大提高了沉陷区土地的利用率。
说明书附图
图1为本发明的预计煤矿开采后地表沉陷示意图;
图2为本发明的沉陷积水区俯视图;
图3为本发明的应用高潜水位采煤沉陷区超前治理方法边界设计示意图;
图4为本发明的应用高潜水位采煤沉陷区超前治理方法填土高程设计示意图;
图中:1-开采工作面,2-煤层,3-表土层,4-沉降后地表,5-未沉降前地表,6-工作面开采区域,7-挖深区域,8-垫浅区域,9-挖深垫浅边界,10-挖深垫浅边界线,11-治理后地表,12-填土设计高程,i-挖深垫浅边界点,p-地表内边界点,f为挖深盆地边界点,i′-i点在走向主断面正上方地表的投影,f′-f点在走向主断面正上方地表的投影,Wmax-走向主断面实际地表下沉最大值。
具体实施方式
下面结合附图对本申请的实施例做进一步说明:
本发明的一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法,首先根据矿区地质采矿条件参数设计高潜水位采沉陷区挖深垫浅后地表高程;计算采煤工作面走向主断面上挖深垫浅的边界点及采空后沉陷盆地的挖深垫浅边界线10;预计采煤工作面采空后形成沉陷区内进行治理需要的填土高程;利用采煤沉陷区内的填土高程,在地下矿区工作面开采前对地表对应区域进行超前治理;当地下开采完成后地表自动形成预计的治理效果,平整地面,从而完成沉陷区挖深垫浅的治理方法;
具体步骤如下
a.如图1和图2所示,根据开采工作面1所在矿区地质采矿条件参数,包括地表地形信息和地下开采计划,预测地表沉陷信息,利用地表沉陷信息结合原始地表地形信息,获取地表沉陷等值线图,设计出沉陷区挖深垫浅后地表高程值计算公式:
利用概率积分法地表沉陷预计,计算开采工作面1走向主断面上方地表点下沉值W(x);选取开采工作面1走向主断面进行挖深垫浅边界设计的参照面,选取走向主断面地表最大下沉值Wmax为最小设计高程,以走向主断面土方量平衡为条件,计算出挖深垫浅后地表高程值中变量b与k的关系式,根据实际工程施工、采煤沉陷地表变形规律、变量b与k的关系式分析设计挖深垫浅后地表高程值公式中的变量k,进而得出挖深垫浅后地表高程值Z(x)=kx+b;
下式为采煤工作面走向主断面正上方地表各点下沉值计算公式W(x),工作面走向主断面正上方地表最大下沉值Wmax:
式中:x为采煤工作面走向主断面正上方地表点相对于工作面的坐标;W0为概率积分法采煤沉陷预计地表最大下沉值;l为地下采煤工作面走向计算长度;L为地下采煤工作面倾向计算长度;D1为地下采煤工作面走向实际长度;Cym为倾向采动程度系数;Cy为倾向采动系数;s为工作面拐点偏移距;r为工作面主要影响半径;erf()为高斯误差函数;
于挖深垫浅后地表高程值Z(x)计算公式中变量b与k关系式的求取:
根据土方量平衡条件:
式中:为挖深垫浅后地表高程值Z(x)与地面平均高程相等的挖深盆地边界点p处的横坐标,点x0地下采煤工作面走向主断面正上方地表下沉值W(x)与挖深垫浅后地表高程值Z(x)相等的点i处的横坐标,为挖深垫浅后地表高程值Z(x)与地表最大下沉值Wmax相等的点f横坐标。
挖深垫浅后地表高程值计算公式中b与k的关系式:
b.为了减小积水区域,最大程度的改造土地,需要计算走向主断面上挖深垫浅的各边界点,获得沉陷盆地的挖深垫浅的地表内边界信息、挖深垫浅边界信息、挖深盆地边界信息:
首先使用二分法找出开采工作面1走向主断面正上方地表下沉值W(x)与挖深垫浅后地表高程值Z(x)相等的点i坐标即挖深垫浅的边界点;挖深垫浅后地表高程值Z(x)与地面原始平均高程H0相等的地表内边界点p即为地表内边界点,利用挖深垫浅后地表高程值公式计算出地表内边界点p横坐标为挖深垫浅后地表高程值Z(x)与地表最大下沉值Wmax相等的点f即为挖深盆地边界点,根据挖深垫浅后地表高程值公式计算出挖深盆地边界点f横坐标为点i′为点i在地表的投影,点f′为点f在地表的投影,则地表处点p,点i′,点f′由采煤工作面走向主断面正上方地表下沉计算公式W(x)可得出上述三点处下沉值W(p),W(i′),W(f′),地表点p,点i′,点f′处下沉等值线即为沉陷盆地的挖深垫浅的地表内边界、挖深垫浅边界、挖深盆地边界;
挖深垫浅后地表高程值Z(x)计算公式中变量k的选取:
k为挖深垫浅后地表高程值计算公式中的变量未知,结合实际工程施工、采煤沉陷地表变形规律、变量b与k的关系式可知:在走向主断面正上方的地表点,地表下沉值与挖深垫浅后地表高程值相等时的点i有且仅有一个点i;在点i左侧W(x)>Z(x),在点i右侧W(x)<Z(x),结合概率积分法主断面地表移动与变形规律分析,k必须大于地表下沉曲线W(x)的最大斜率,即大于x方向的地表最大倾斜值,即根据挖深垫浅后地表高程值计算公式中变量b与k的关系式分析,截距b随挖深垫浅后地表高程曲线的斜率k增大而减小,当k偏大时,此时两曲线交点偏上,根据实际验证,明显不符合挖填方平衡原则。综上分析,挖深垫浅后地表高程曲线斜率k的取值需略大于地表下沉曲线最大斜率;
如图3和图4所示,由此得出k的取值,进而由b与k的关系式得出b的值,最后求出挖深垫浅挖深垫浅后地表高程曲线Z(x)=kx+b;同时挖深垫浅后地表高程值Z(x)计算公式中变量k的选取,使得采煤沉陷区内积水区与地面之间存在适宜的坡度,更有利于沉陷区水土保持与农业利用。
c.在开采工作面1开采前或沉陷区稳沉前,根据填土高程设计进行施工,填土高程施工设计分为填土设计高程H1与治理后的最终高程H2,在工作面走向主断面与地表相交处的每个地表点上的填土设计高程H1就等于地表原始高程H0加上挖深垫浅边界9高程值Z(x)减去该点处地表移动下沉值W(x),值为正时为填方高程,值为负时为挖方高程;在整个沉陷盆地内,根据沉陷预计等值线图与走向主断面的挖填方高程的设计,下沉值相同处的高程设计相同;治理后的最终高程H2等于填土设计高程H1减去地表移动下沉值W(x),完成沉陷区的挖深垫浅施工设计;具体的填土设计高程H1与治理后的最终高程H2如下表
表1填土设计高程H1与治理后的最终高程H2计算公式
d.利用沉陷区的挖深垫浅施工设计,对整个地表沉陷盆地的挖深区域7和垫浅区域8进行改造,根据填土高程设计施工,地表点i′处的沉陷等值线即为挖深垫浅边界,边界内部沉陷盆地为挖深区域7,边界外侧沉陷边缘为垫浅区域8;在挖深区域7,地表点f′处的沉陷等值线即为挖深盆地边界,挖深盆地边界内部开挖区填土设计高程H1=H0+W(x)-Wmax,平整开发区底部地面;在挖深盆地边界与挖深垫浅边界之间的环形区域,开挖区填土设计高程H1=H0+W(x)-Z(x),将挖掘出的土方填入垫浅区;在垫浅区域8,地表点p处的沉陷等值线为地表内边界,在地表内边界与挖深垫浅边界之间的环形区域,垫浅区填土设计高程H1=H0+W(x)-Z(x),在沉陷边界与地表内边界之间的环形区域,垫浅区填土设计高程H1=H0+W(x),平整垫浅区域8地面,从而完成对治理矿区采煤工作面上方沉陷区的地质区域改造。
在工作面开采前或沉陷区稳沉前进行沉陷区的挖深垫浅填土设计高程与施工;通常为了使土体不被闲置,沉陷区稳沉后进行施工挖深垫浅,但当高潜水位矿区在地表沉陷时沉陷区产生积水时,需要在工作面开采前或沉陷区稳沉前进行沉陷区的挖深垫浅填土设计高程与施工,避免了沉陷期土地闲置及积水影响施工的问题;常规的边采边复高程模型中高程受开采进度影响,在工作面开采前或沉陷区稳沉前根据开采沉陷设计填土高程,治理后的地表高程不受开采工作的影响。
实施例:
a.设计高潜水位采沉陷区挖深垫浅后地表高程:根据矿区地质采矿条件参数,包括开采工作面1、煤层2、表土层3和未沉降前地表5的信息利用概率积分法地表沉陷预计沉降后地表4,计算地下工作面开采区域6走向主断面上方地表点下沉值W(x);选取采煤工作面走向主断面进行挖深垫浅边界设计的参照面,选取走向主断面地表最大下沉值Wmax为最小设计高程,地下采煤工作面走向主断面上方地表各点挖深垫浅后地表高程值Z(x),计算公式为Z(x)=kx+b;以走向主断面土方量平衡为条件,计算出挖深垫浅后地表高程值中变量b与k的关系式,根据实际工程施工、采煤沉陷地表变形规律、变量b与k的关系式分析设计挖深垫浅后地表高程值公式中的变量k,进而得出挖深垫浅后地表高程值计算公式Z(x)=kx+b。
煤层2平均采厚m为3.1m,平均倾角α为10°,平均采深H为357m,其走向长为D1为575m,倾向斜宽为D2为150m,下沉系数q为0.98,主要影响角正切tanβ为2.16,最大下沉角θ0为88°,主要影响传播角θ为86.3°,拐点偏移距经验公式s/H等于0.07。
得出工作面走向主断面正上方地表下沉的计算公式W(x),工作面走向主断面最大下沉值Wmax,其中Cy为倾向采动系数:
Cy=0.8686
设计挖深垫浅后地表高程值:
Z(x)=kx+b (3)
土方量平衡条件:
式中:为挖深垫浅后地表高程值Z(x)与地面平均高程相等的p处的横坐标,点x0地下采煤工作面走向主断面正上方地表下沉值W(x)与挖深垫浅后地表高程值Z(x)相等的点i处的横坐标,为挖深垫浅后地表高程值Z(x)与地表最大下沉值Wmax相等的点f横坐标。
挖深垫浅后地表高程值计算公式中b与k的关系式:
k为挖深垫浅后地表高程值计算公式中的变量未知,结合实际工程施工、采煤沉陷地表变形规律、变量b与k的关系式可知:1)在走向主断面正上方的地表点,地表下沉值与挖深垫浅后地表高程值相等时的点i有且仅有一个点i;2)在点i左侧W(x)>Z(x),在点i右侧W(x)<Z(x),结合概率积分法主断面地表移动与变形规律分析,k必须大于地表下沉曲线W(x)的最大斜率,即大于x方向的地表最大倾斜值,即3)根据挖深垫浅后地表高程值计算公式中变量b与k的关系式分析,截距b随挖深垫浅后地表高程曲线的斜率k增大而减小,当k偏大时,此时两曲线交点偏上,根据实际验证,明显不符合挖填方平衡原则。综上分析,挖深垫浅后地表高程曲线斜率k的取值需略大于地表下沉曲线最大斜率。
b.计算走向主断面上挖深垫浅的边界点及沉陷盆地的挖深垫浅边界线10:使用二分法找出地下采煤工作面走向主断面正上方地表下沉值W(x)与挖深垫浅后地表高程值Z(x)相等的点i坐标等于(25,1.0),即挖深垫浅的边界点;挖深垫浅后地表高程值Z(x)与地面原始平均高程H0相等的点p即为地表内边界点,根据挖深垫浅高程公式计算出点p横坐标为挖深垫浅后地表高程值Z(x)与地表最大下沉值Wmax相等的点f即为挖深盆地边界点,根据挖深垫浅高程公式计算出点f横坐标为将上述三个交点pif处横坐标带入采煤工作面走向主断面正上方地表下沉计算公式W(x)可得出三点处下沉值W(p),W(i),W(f)分别等于0.4m,1.0m,1.6m;p、i、f处下沉等值线即为沉陷盆地的挖深垫浅的地表内边界、挖深垫浅边界、挖深盆地边界。
c.采煤沉陷区内的填土高程施工设计:在地下采煤工作面开采前或沉陷区稳沉前,根据填土高程设计进行施工,填土高程施工设计分为填土设计高程H112与治理后的最终高程H2,在走向主断面正上方地表每个点上的填土设计高程H112就等于地表原始平均高程H0加上挖深垫浅边界高程值Z(x)减去该点处地表移动下沉值W(x),值为正时为填方高程,值为负时为挖方高程;在整个沉陷盆地内,根据沉陷预计等值线图与走向主断面的挖填方高程的设计,下沉值相同处的高程设计相同;治理后的最终高程H2等于填土设计高程H112减去该点处地表移动下沉值W(x);
d.利用沉陷区的挖深垫浅施工设计,对整个地表沉陷盆地的挖深区域7和垫浅区域8进行改造,根据填土高程设计施工,地表点i′处的沉陷等值线即为挖深垫浅边界9,边界内部沉陷盆地为挖深区域7,边界外侧沉陷边缘为垫浅区域8;在挖深区域7,地表点f′处的沉陷等值线即为挖深盆地边界,挖深盆地边界内部开挖区填土设计高程H112等于H0+W(x)-Wmax,平整开发区底部地面;在挖深盆地边界与挖深垫浅边界9之间的环形区域,开挖区填土设计高程H112等于H0+W(x)-Z(x),将挖掘出的土方填入垫浅区;在垫浅区域8,地表点p处的沉陷等值线为地表内边界,在地表内边界与挖深垫浅边界9之间的环形区域,垫浅区填土设计高程H112等于H0+W(x)-Z(x),在沉陷边界与地表内边界之间的环形区域,垫浅区填土设计高程H112等于H0+W(x),平整垫浅区域8地面,从而完成对治理矿区采煤工作面上方沉陷区的地质区域治理后地表的11改造。
Claims (4)
1.一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法,其特征在于:首先根据矿区地质采矿条件参数设计高潜水位采煤沉陷区挖深垫浅后地表高程;计算采煤工作面走向主断面上挖深垫浅的边界点及采空后沉陷盆地的挖深垫浅边界线;预计采煤工作面采空后形成沉陷区内进行治理需要的填土高程;利用采煤沉陷区内的填土高程,在地下矿区工作面开采前对地表对应区域进行超前治理;当地下开采完成后地表自动形成预计的治理效果,平整地面,从而完成沉陷区挖深垫浅的治理方法;
具体步骤如下
a.根据开采工作面1所在矿区地质采矿条件参数,包括地表地形信息和地下开采计划,预测地表沉陷信息,利用地表沉陷信息结合原始地表地形信息,获取地表沉陷等值线图,设计出沉陷区挖深垫浅后地表高程值计算公式:
利用概率积分法地表沉陷预计,计算开采工作面1走向主断面上方地表点下沉值W(x);选取开采工作面1走向主断面进行挖深垫浅边界设计的参照面,选取走向主断面地表最大下沉值Wmax为最小设计高程,以走向主断面土方量平衡为条件,计算出挖深垫浅后地表高程值中变量b与k的关系式,根据实际工程施工、采煤沉陷地表变形规律、变量b与k的关系式分析设计挖深垫浅后地表高程值公式中的变量k,进而得出挖深垫浅后地表高程值Z(x)=kx+b;
下式为采煤工作面走向主断面正上方地表各点下沉值计算公式W(x),工作面走向主断面正上方地表最大下沉值Wmax:
式中:x为采煤工作面走向主断面正上方地表点相对于工作面的坐标;W0为概率积分法采煤沉陷预计地表最大下沉值;l为地下采煤工作面走向计算长度;L为地下采煤工作面倾向计算长度;D1为地下采煤工作面走向实际长度;Cym为倾向采动程度系数;Cy为倾向采动系数;s为工作面拐点偏移距;r为工作面主要影响半径;erf()为高斯误差函数;
b.为了减小积水区域,最大程度的改造土地,需要计算走向主断面上挖深垫浅的各边界点,获得沉陷盆地的挖深垫浅的地表内边界信息、挖深垫浅边界信息、挖深盆地边界信息:
首先使用二分法找出开采工作面1走向主断面正上方地表下沉值W(x)与挖深垫浅后地表高程值Z(x)相等的点i坐标即挖深垫浅的边界点;挖深垫浅后地表高程值Z(x)与地面原始平均高程H0相等的地表内边界点p即为地表内边界点,利用挖深垫浅后地表高程值公式计算出地表内边界点p横坐标为挖深垫浅后地表高程值Z(x)与地表最大下沉值Wmax相等的点f即为挖深盆地边界点,根据挖深垫浅后地表高程值公式计算出挖深盆地边界点f横坐标为点i′为点i在地表的投影,点f′为点f在地表的投影,则地表处点p,点i′,点f′由采煤工作面走向主断面正上方地表下沉计算公式W(x)可得出上述三点处下沉值W(p),W(i′),W(f′),地表点p,点i′,点f′处下沉等值线即为沉陷盆地的挖深垫浅的地表内边界、挖深垫浅边界、挖深盆地边界;
c.在开采工作面1开采前或沉陷区稳沉前,根据填土高程设计进行施工,填土高程施工设计分为填土设计高程H1与治理后的最终高程H2,在工作面走向主断面与地表相交处的每个地表点上的填土设计高程H1就等于地表原始高程H0加上挖深垫浅边界高程值Z(x)减去该点处地表移动下沉值W(x),值为正时为填方高程,值为负时为挖方高程;在整个沉陷盆地内,根据沉陷预计等值线图与走向主断面的挖填方高程的设计,下沉值相同处的高程设计相同;治理后的最终高程H2等于填土设计高程H1减去地表移动下沉值W(x),完成沉陷区的挖深垫浅施工设计;
d.利用沉陷区的挖深垫浅施工设计,对整个地表沉陷盆地的挖深区域和垫浅区域进行改造,根据填土高程设计施工,地表点i′处的沉陷等值线即为挖深垫浅边界,边界内部沉陷盆地为挖深区域,边界外侧沉陷边缘为垫浅区域;在挖深区域,地表点f′处的沉陷等值线即为挖深盆地边界,挖深盆地边界内部开挖区填土设计高程H1=H0+W(x)-Wmax,平整开发区底部地面;在挖深盆地边界与挖深垫浅边界之间的环形区域,开挖区填土设计高程H1=H0+W(x)-Z(x),将挖掘出的土方填入垫浅区;在垫浅区域,地表点p处的沉陷等值线为地表内边界,在地表内边界与挖深垫浅边界之间的环形区域,垫浅区填土设计高程H1=H0+W(x)-Z(x),在沉陷边界与地表内边界之间的环形区域,垫浅区填土设计高程H1=H0+W(x),平整垫浅区域地面,从而完成对治理矿区采煤工作面上方沉陷区的地质区域改造。
3.根据权利要求1或2所述的基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法,其特征在于挖深垫浅后地表高程值Z(x)计算公式中变量k的选取:
k为挖深垫浅后地表高程值计算公式中的变量未知,结合实际工程施工、采煤沉陷地表变形规律、变量b与k的关系式可知:在走向主断面正上方的地表点,地表下沉值与挖深垫浅后地表高程值相等时的点i有且仅有一个点i;在点i左侧W(x)>Z(x),在点i右侧W(x)<Z(x),结合概率积分法主断面地表移动与变形规律分析,k必须大于地表下沉曲线W(x)的最大斜率,即大于x方向的地表最大倾斜值,即根据挖深垫浅后地表高程值计算公式中变量b与k的关系式分析,截距b随挖深垫浅后地表高程曲线的斜率k增大而减小,当k偏大时,此时两曲线交点偏上,根据实际验证,明显不符合挖填方平衡原则;综上分析,挖深垫浅后地表高程曲线斜率k的取值需略大于地表下沉曲线最大斜率;
由此得出k的取值,进而由b与k的关系式得出b的值,最后求出挖深垫浅后地表高程曲线Z(x)=kx+b;同时挖深垫浅后地表高程值Z(x)计算公式中变量k的选取,使得采煤沉陷区内积水区与地面之间存在适宜的坡度,更有利于沉陷区水土保持与农业利用。
4.根据权利要求1所述的基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法,其特征在于:在工作面开采前或沉陷区稳沉前进行沉陷区的挖深垫浅填土设计高程与施工;通常为了使土体不被闲置,沉陷区稳沉后进行施工挖深垫浅,但当高潜水位矿区在地表沉陷时沉陷区产生积水时,需要在工作面开采前或沉陷区稳沉前进行沉陷区的挖深垫浅填土设计高程与施工,避免了沉陷期土地闲置及积水影响施工的问题;常规的边采边复高程模型中高程受开采进度影响,在工作面开采前或沉陷区稳沉前根据开采沉陷设计填土高程,治理后的地表高程不受开采工作的影响。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010012470.5A CN111101941B (zh) | 2020-01-07 | 2020-01-07 | 一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010012470.5A CN111101941B (zh) | 2020-01-07 | 2020-01-07 | 一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111101941A CN111101941A (zh) | 2020-05-05 |
CN111101941B true CN111101941B (zh) | 2021-04-30 |
Family
ID=70425729
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010012470.5A Active CN111101941B (zh) | 2020-01-07 | 2020-01-07 | 一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111101941B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112282757B (zh) * | 2020-10-26 | 2021-09-21 | 中国矿业大学(北京) | 一种多煤层时空协调开采布局方法 |
CN113530597B (zh) * | 2021-07-22 | 2023-05-16 | 淮北矿业股份有限公司 | 一种应力阻隔式采矿沉陷区地表构筑物防护方法 |
CN113836740B (zh) * | 2021-09-30 | 2024-03-19 | 安徽大学 | 一种高潜水位采煤沉陷积水区历史空间信息计算方法 |
CN115829133B (zh) * | 2022-12-12 | 2023-07-18 | 通用技术集团工程设计有限公司 | 一种充填工艺回收条带煤柱的地表沉陷预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2387836C1 (ru) * | 2008-12-29 | 2010-04-27 | Институт угля и углехимии Сибирского отделения Российской Академии Наук (ИУУ СО РАН) | Способ открытой добычи угля с использованием инфраструктуры смежной шахты |
CN102726187A (zh) * | 2012-06-25 | 2012-10-17 | 西北综合勘察设计研究院 | 一种煤矿开采沉陷区地表生态治理方法 |
CN103047970A (zh) * | 2012-12-18 | 2013-04-17 | 中国矿业大学(北京) | 一种确定塌陷耕地损害边界的方法 |
CN107023297A (zh) * | 2016-02-02 | 2017-08-08 | 中国矿业大学(北京) | 煤矿巨型塌陷坑的修复方法 |
CN107882561A (zh) * | 2017-11-10 | 2018-04-06 | 安徽省交通航务工程有限公司 | 一种高潜水位采煤沉陷区超前治理方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10960444B2 (en) * | 2018-04-06 | 2021-03-30 | Karl William Yost | Closure methods for mines |
-
2020
- 2020-01-07 CN CN202010012470.5A patent/CN111101941B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2387836C1 (ru) * | 2008-12-29 | 2010-04-27 | Институт угля и углехимии Сибирского отделения Российской Академии Наук (ИУУ СО РАН) | Способ открытой добычи угля с использованием инфраструктуры смежной шахты |
CN102726187A (zh) * | 2012-06-25 | 2012-10-17 | 西北综合勘察设计研究院 | 一种煤矿开采沉陷区地表生态治理方法 |
CN103047970A (zh) * | 2012-12-18 | 2013-04-17 | 中国矿业大学(北京) | 一种确定塌陷耕地损害边界的方法 |
CN107023297A (zh) * | 2016-02-02 | 2017-08-08 | 中国矿业大学(北京) | 煤矿巨型塌陷坑的修复方法 |
CN107882561A (zh) * | 2017-11-10 | 2018-04-06 | 安徽省交通航务工程有限公司 | 一种高潜水位采煤沉陷区超前治理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111101941A (zh) | 2020-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111101941B (zh) | 一种基于边界优化的高潜水位采煤沉陷区挖深垫浅超前治理方法 | |
CN102609782B (zh) | 基于采煤沉陷模拟和复耕率的边采边复时机的优选方法 | |
CN107882561A (zh) | 一种高潜水位采煤沉陷区超前治理方法 | |
CN109372511B (zh) | 一种利用黄河泥沙动态充填采煤塌陷地的方法 | |
CN109555556B (zh) | 一种面向耕地保护的固体充填采煤充实率设计方法 | |
CN111750822A (zh) | 一种采煤诱发的覆岩与地表沉陷协同动态预测方法 | |
CN110174506A (zh) | 一种喀斯特地区土壤有机碳估算方法 | |
CN105093307B (zh) | 下古生界斜地层真厚度求取方法 | |
CN111798052B (zh) | 一种高潜水位采煤沉陷积水区三维空间信息动态预测方法 | |
CN105551083B (zh) | 一种刻画水平井井眼轨迹与地层关系的方法 | |
CN102968806A (zh) | 一种不规则边界区域内自适应网格的生成方法 | |
CN109063285B (zh) | 一种土坡钻孔布置方案设计方法 | |
CN113077088A (zh) | 城中采煤沉陷区水土资源空间格局重构方法 | |
CN111551932A (zh) | 准确获取开采影响边界及确定建筑损害等级的方法 | |
CN103541376B (zh) | 采煤沉陷区地基在重复开采条件下的基础变形预测方法 | |
CN105975782A (zh) | 一种复合软弱起伏基底排土场边坡形态优化方法 | |
CN111257870A (zh) | 一种利用InSAR监测数据的采煤沉陷积水区水下地形反演方法 | |
CN103235853A (zh) | 用于采煤沉陷地边采边复的基于土方平衡的基塘布局方法 | |
CN111400659A (zh) | 采-选-充一体化开采技术耕地保护率计算方法 | |
CN107131867A (zh) | 厚松散层条件煤层开采地表下沉率交迭双波兹曼拟合方法 | |
CN104541643A (zh) | 一种平原区的采煤沉陷区土地整治方法 | |
CN117706655A (zh) | 基于气藏气势梯度判断输导模型的方法 | |
CN115344928B (zh) | 一种耕地下采煤区域变形控制的地表沉陷叠加计算方法 | |
RU2360119C1 (ru) | Способ естественного возобновления болото- и торфообразовательного процесса при разработке торфяных месторождений | |
Shi et al. | Study on numerical models in predicting surface deformation caused by underground coal mining |
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 | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 230000 No.20 Tongling Road, Yaohai District, Hefei City, Anhui Province Patentee after: Anhui Construction Engineering Transportation and Navigation Group Co.,Ltd. Patentee after: ANHUI University Address before: No. 20 Tongling Road, Yaohai District, Hefei City, Anhui Province, 230601 Patentee before: ANHUI PROVINCE TRANSPORTATION HARBOR ENGINEERING CO.,LTD. Patentee before: ANHUI University |