CN105808935A - 多煤层开采下地面积水时间的确定方法 - Google Patents
多煤层开采下地面积水时间的确定方法 Download PDFInfo
- Publication number
- CN105808935A CN105808935A CN201610121402.6A CN201610121402A CN105808935A CN 105808935 A CN105808935 A CN 105808935A CN 201610121402 A CN201610121402 A CN 201610121402A CN 105808935 A CN105808935 A CN 105808935A
- Authority
- CN
- China
- Prior art keywords
- mining
- coal seam
- coal
- seam
- ground
- 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
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Revetment (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及多煤层开采下地面积水时间的确定方法,属于采矿技术、土地利用和土地复垦领域。该方法包括:根据地下煤炭赋存条件,利用概率积分法分别预计各煤层开采引起的地面下沉,然后依据煤层开采计划,确定多煤层重复采动下地表累积下沉量,同时结合地表原始高程,模拟多煤层开采地面沉陷情况,最后根据当地潜水位埋深,确定多煤层开采地面积水时间。本发明综合考虑了多煤层重复采动下各煤层开采引起的地表下沉、地面原始高程、地下潜水位标高,从而确定了多煤层开采下地面积水的时间,简单易懂,便于实际应用,为后续多煤层开采沉陷地土地复垦工作提供了基础信息。
Description
技术领域
本发明属于采矿技术、土地利用和土地复垦技术领域,特别涉及中高潜水位平原矿区多煤层重复采动下地面积水时间的确定方法。
背景技术
我国煤炭资源储量丰富,已探明的煤炭资源储量高达1.48万亿吨。其中,部分地区煤层层数多、累积可采厚度大,例如安徽省淮南矿区可采煤层9-18层,累积可采厚度达25-34m;山东省兖州、济宁等矿区可采煤层大部分在7-10层,可采煤层总厚度约8-12m;江苏省徐州矿区主要可采煤层9层,累积采厚达15m。众所周知,我国是煤炭生产大国,据《2014年国民经济和社会发展统计公报》数据显示,2014年全年我国原煤产量为38.7亿吨,其中90%以上的煤炭产量来自于井工煤矿,而且大多采用走向长壁式全部跨落法管理顶板,从而形成了众多的采煤沉陷土地,据相关资料,到2012年底,我国因采煤而沉陷的土地面积约为156万公顷,采煤沉陷地复垦工作迫在眉睫。
随着煤炭的不断开采和开采深度的增加,多煤层开采更加普遍。多煤层开采下,由于地表受到多次重复扰动,地面损毁情况复杂多变,不利于后续土地复垦工作的进行。尤其是在中高潜水位平原矿区,多煤层开采后沉陷深度大,加之当地较高的地下潜水位,采煤沉陷后,地面极易形成积水,从而丧失大面积耕地而无法挽救。而且研究表明,土壤浸水后,其中部分营养物质将会有所流失,不利于复垦后耕地生产力的迅速恢复。因此,为了保护珍贵的表土资源,首先需要确定多煤层开采下地面出现积水的时间,才能提前安排对表土的保护和剥离措施,同时还可以为后续复垦提供优质的表土资源。
目前,对于单一工作面或单一煤层开采的地面损毁及积水等有了大量深入的探索和研究,并取得了丰硕的成果。但对于多煤层开采地面积水情况的研究较少,在多煤层重复采动下,如何确定地面积水时间,从而为表土剥离工作提供基础信息,对后续多煤层开采沉陷地复垦至关重要。
发明内容
本发明的目的是为解决上述难题,提出一种多煤层开采下地面积水时间的确定方法,本发明通过模拟分析多煤层重复采动影响下地面沉陷情况,同时结合当地潜水位埋深,确定地表积水时间,从而为后续多煤层开采土地复垦工作提供基础信息。
本发明提出的多煤层开采下地面积水时间的确定方法,该方法主要适用于中高潜水位多煤层开采矿区,其特征在于,该方法包括:根据地下煤炭赋存条件,利用概率积分法分别预计各煤层开采引起的地面下沉,然后依据煤层开采计划,确定多煤层重复采动下地表累积下沉量,同时结合地表原始高程,模拟多煤层开采地面沉陷情况,最后根据当地潜水位埋深,确定多煤层开采地面积水时间。
该方法具体包括以下步骤:
1)根据地下煤炭赋存条件,利用概率积分法分别预计各煤层开采引起的地面下沉:设开采区域共n个煤层,且开采顺序为煤层1、煤层2、…、煤层i、…、煤层n,地面任意点P的坐标为(x,y),则煤层i开采引起点P的下沉值Wi(x,y)如表达式(1)所示:
式(1)中:Wi(x,y)为煤层i开采引起地面任意点P的下沉值,单位mm;
qi为煤层i开采的下沉系数;mi为煤层i的煤厚,单位mm;
αi为煤层i的煤层倾角,单位°;Si为煤层i的开采区域;
ri为煤层i的主要影响半径,单位m;i=1,2,3,…,i,…n;
2)根据煤层开采计划,确定多煤层重复采动下地表累积下沉量:m个煤层重复采动引起的地表累积下沉如表达式(2)所示:
式(2)中:为m个煤层重复采动引起地面任意点P的累积下沉值,单位m;
Wj(x,y)为煤层j开采引起地面任意点P的下沉值,单位mm;
m=1,2,3,…,m,…n;j=1,2,3,…,j,…m;
3)结合地表原始高程,模拟多煤层开采地面沉陷:将多煤层开采地表累积下沉与煤炭开采前地表的原始高程进行叠加分析,模拟多煤层开采地面沉陷情况,m个煤层开采结束后,地面任意点P的高程如表达式(3)所示:
式(3)中:Hm(x,y)为m个煤层开采结束后地面任意点P的高程,单位m;
H0(x,y)为煤炭开采前地面任意点P的原始高程,单位m;
m=1,2,3,…,m,…n;
4)根据当地潜水位埋深,确定多煤层开采地面积水时间:由于采煤沉陷影响,如果t-1个煤层开采结束后,地面任意点P的高程高于地下潜水位高程,且t个煤层开采结束后,点P的高程等于或低于地下潜水位时,地面任意点P将在t个煤层开采结束后出现积水,则
式(4)中:Ht(x,y)为t个煤层开采结束后地面任意点P的高程,单位m;
Ht-1(x,y)为t-1个煤层开采结束后地面任意点P的高程,单位m;
HQ为地下潜水位高程,单位m。
本发明主要具有以下技术优点:
本发明通过利用概率积分法预计了多煤层开采下,各煤层开采引起的地表下沉,同时结合煤炭开采前地面的原始高程,准确的模拟分析了地面采煤沉陷情况,在此基础上,结合当地潜水位埋深,综合确定多煤层开采地面积水时间。该方法简单易懂,便于实际应用,为后续多煤层开采土地复垦工作提供了基础信息。
附图说明
图1为多煤层开采下地面积水时间确定方法流程图。
图2为实施例四煤层开采地面点P1、P2下沉示意图。
具体实施方式
本发明提出的多煤层开采下地面积水时间的确定方法,结合附图及实施例详细说明如下:
本发明的方法主要适用于中高潜水位多煤层开采矿区,包括根据地下煤炭赋存条件,利用概率积分法别预计各煤层开采引起的地面下沉,然后依据煤层开采计划,确定多煤层重复采动下地表累积下沉量,同时结合地表原始高程,模拟多煤层开采地面沉陷情况,最后根据当地潜水位埋深,确定多煤层开采地面积水时间
该方法具体包括以下步骤:
1)根据地下煤炭赋存条件,利用概率积分法分别预计各煤层开采引起的地面下沉:设开采区域共n个煤层,且开采顺序为煤层1、煤层2、…、煤层i、…、煤层n,地面任意点P的坐标为(x,y),则煤层i开采引起点P的下沉值Wi(x,y)如表达式(1)所示:
式(1)中:Wi(x,y)为煤层i开采引起地面任意点P的下沉值,单位mm;
qi为煤层i开采的下沉系数;mi为煤层i的煤厚,单位mm;
αi为煤层i的煤层倾角,单位°;Si为煤层i的开采区域;
ri为煤层i的主要影响半径,单位m;i=1,2,3,…,i,…n;
2)根据煤层开采计划,确定多煤层重复采动下地表累积下沉量:m个煤层重复采动引起的地表累积下沉如表达式(2)所示:
式(2)中:为m个煤层重复采动引起地面任意点P的累积下沉值,单位m;
Wj(x,y)为煤层j开采引起地面任意点P的下沉值,单位mm;
m=1,2,3,…,m,…n;j=1,2,3,…,j,…m;
3)结合地表原始高程,模拟多煤层开采地面沉陷:将多煤层开采地表累积下沉与煤炭开采前地表的原始高程进行叠加分析,模拟多煤层开采地面沉陷情况,m个煤层开采结束后,地面任意点P的高程如表达式(3)所示:
式(3)中:Hm(x,y)为m个煤层开采结束后地面任意点P的高程,单位m;
H0(x,y)为煤炭开采前地面任意点P的原始高程,单位m;
m=1,2,3,…,m,…n;
4)根据当地潜水位埋深,确定多煤层开采地面积水时间:由于采煤沉陷影响,如果t-1个煤层开采结束后,地面任意点P的高程高于地下潜水位高程,且t个煤层开采结束后,点P的高程等于或低于地下潜水位时,地面任意点P将在t个煤层开采结束后出现积水,则
式(4)中:Ht(x,y)为t个煤层开采结束后地面任意点P的高程,单位m;
Ht-1(x,y)为t-1个煤层开采结束后地面任意点P的高程,单位m;
HQ为地下潜水位高程,单位m。
实施例:
本实施例为东部某一高潜水位平原矿区,地面平坦开阔,地表自然标高在+22.1m~+23.6m之间,平均+22.7m,地下潜水位埋深为2.0m左右,区内主要有四个开采煤层,分别为煤层1、煤层2、煤层3和煤层4,这四个煤层均为全区可采,且结构简单,其中煤层1平均厚度为2.0m,煤层2为2.5m,煤层3为2.0m,煤层4为3.0m。
本实施例在计算机软件环境下模拟确定多煤层开采下地面积水时间,其方法流程如图1所示,包括以下步骤:
1)根据地下煤炭赋存条件,利用概率积分法分别预计各煤层开采引起的地面下沉:本实施例中,开采区域共4个煤层,且开采顺序为煤层1、煤层2、煤层3、煤层4,地面点P1的坐标为(362,211),P2的坐标为(538,211),煤层1开采引起点P1、P2的下沉值分别为W1(362,211)=700mm、W1(538,211)=500mm,煤层2开采引起点P1、P2的下沉值分别为W2(362,211)=820mm、W2(538,211)=510mm,煤层3开采引起点P1、P2的下沉值分别为W3(362,211)=710mm、W3(538,211)=500mm,煤层4开采引起点P1、P2的下沉值分别为W4(362,211)=9200mm、W4(538,211)=530mm。
2)根据煤层开采计划,确定多煤层重复采动下地表累积下沉量:煤层1开采引起点P1、P2的下沉分别为 煤层2开采结束后,引起的地表累积下沉 煤层3开采结束后,引起的地表累积下沉 煤层4开采结束后,引起的地表累积下沉
3)结合地表原始高程,模拟多煤层开采地面沉陷:将多煤层开采地表累积下沉与煤炭开采前地表的原始高程进行叠加分析,模拟多煤层开采地面采煤沉陷情况,本实施例中,煤炭开采前点P1、P2的原始高程分别为H0(362,211)=22.4m、H0(538,211)=22.7m,煤层1开采结束后,点P1、P2的高程分别为H1(362,211)=(22.4-0.71)m=21.69m、H1(538,211)=(22.6-0.49)m=22.11m;煤层2开采结束后,点P1、P2的高程分别为H2(362,211)=(22.4-1.53)m=20.87m、H2(538,211)=(22.6-1.0)m=21.6m;煤层3开采结束后,点P1、P2的高程分别为H3(362,211)=(22.4-2.25)m=20.15m、H3(538,211)=(22.6-1.5)m=21.1m;煤层4开采结束后,点P1、P2的高程分别为H4(362,211)=(22.4-3.15)m=19.25m、H2(538,211)=(22.6-2.07)m=20.53m。
4)根据当地潜水位埋深,确定多煤层开采地面积水时间:由于采煤沉陷影响,如果t-1个煤层开采结束后,地面任意点P的高程高于地下潜水位高程,且t个煤层开采结束后,点P的高程等于或低于地下潜水位时,地面任意点P将在煤层t开采结束后出现积水,本实施例中,地下潜水位高程HQ=22.7-2.0=20.7m,由于H3(362,211)≤HQ<H2(362,211),H4(538,211)≤HQ<H3(538,211),所以当煤层3开采结束后,地面点P1将出现积水,当煤层4开采结束后,地面点P2将出现积水。多煤层开采地面点P下沉剖面如图2所示,图中A-B表示原始地面,C表示当地潜水位,S1、S2、S3、S4分别表示煤层1、煤层2、煤层3、煤层4,曲线L1、L2、L3、L4分别表示煤层1、煤层2、煤层3、煤层4开采结束后下沉盆地剖面,区域R1、R2、R3、R4分别表示煤层1、煤层2、煤层3、煤层4开采结束后地表下沉盆地范围。
Claims (2)
1.一种多煤层开采下地面积水时间的确定方法,其特征在于,该方法包括:根据地下煤炭赋存条件,利用概率积分法分别预计各煤层开采引起的地面下沉,然后依据煤层开采计划,确定多煤层重复采动下地表累积下沉量,同时结合地表原始高程,模拟多煤层开采地面沉陷情况,最后根据当地潜水位埋深,确定多煤层开采地面积水时间。
2.如权利要求1所述的方法,其特征在于,该方法具体包括以下步骤:
1)根据地下煤炭赋存条件,利用概率积分法分别预计各煤层开采引起的地面下沉:设开采区域共n个煤层,且开采顺序为煤层1、煤层2、…、煤层i、…、煤层n,地面任意点P的坐标为(x,y),则煤层i开采引起点P的下沉值Wi(x,y)如表达式(1)所示:
式(1)中:Wi(x,y)为煤层i开采引起地面任意点P的下沉值,单位mm;
qi为煤层i开采的下沉系数;mi为煤层i的煤厚,单位mm;
αi为煤层i的煤层倾角,单位°;Si为煤层i的开采区域;
ri为煤层i的主要影响半径,单位m;i=1,2,3,…,i,…n;
2)根据煤层开采计划,确定多煤层重复采动下地表累积下沉量:m个煤层重复采动引起的地表累积下沉如表达式(2)所示:
式(2)中:为m个煤层重复采动引起地面任意点P的累积下沉值,单位m;
Wj(x,y)为煤层j开采引起地面任意点P的下沉值,单位mm;
m=1,2,3,…,m,…n;j=1,2,3,…,j,…m;
3)结合地表原始高程,模拟多煤层开采地面沉陷:将多煤层开采地表累积下沉与煤炭开采前地表的原始高程进行叠加分析,模拟多煤层开采地面沉陷情况,m个煤层开采结束后,地面任意点P的高程如表达式(3)所示:
式(3)中:Hm(x,y)为m个煤层开采结束后地面任意点P的高程,单位m;
H0(x,y)为煤炭开采前地面任意点P的原始高程,单位m;
m=1,2,3,…,m,…n;
4)根据当地潜水位埋深,确定多煤层开采地面积水时间:由于采煤沉陷影响,如果t-1个煤层开采结束后,地面任意点P的高程高于地下潜水位高程,且t个煤层开采结束后,点P的高程等于或低于地下潜水位时,地面任意点P将在t个煤层开采结束后出现积水,则
式(4)中:Ht(x,y)为t个煤层开采结束后地面任意点P的高程,单位m;
Ht-1(x,y)为t-1个煤层开采结束后地面任意点P的高程,单位m;
HQ为地下潜水位高程,单位m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610121402.6A CN105808935B (zh) | 2016-03-03 | 2016-03-03 | 多煤层开采下地面积水时间的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610121402.6A CN105808935B (zh) | 2016-03-03 | 2016-03-03 | 多煤层开采下地面积水时间的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105808935A true CN105808935A (zh) | 2016-07-27 |
CN105808935B CN105808935B (zh) | 2018-08-10 |
Family
ID=56466129
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610121402.6A Active CN105808935B (zh) | 2016-03-03 | 2016-03-03 | 多煤层开采下地面积水时间的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105808935B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257870A (zh) * | 2020-02-26 | 2020-06-09 | 安徽大学 | 一种利用InSAR监测数据的采煤沉陷积水区水下地形反演方法 |
CN112749477A (zh) * | 2020-11-30 | 2021-05-04 | 安徽理工大学 | 厚表土薄基岩煤矿立井井筒偏斜量的预计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101422094A (zh) * | 2007-11-02 | 2009-05-06 | 煤炭科学研究总院唐山研究院 | 平原矿区厚煤层开采塌陷地动态预复垦方法 |
CN101709645A (zh) * | 2009-12-17 | 2010-05-19 | 煤炭科学研究总院唐山研究院 | 一种矿区厚煤层开采沉陷区就地取土复垦方法 |
CN104564069A (zh) * | 2015-01-09 | 2015-04-29 | 中国矿业大学(北京) | 一种基于方格网法的地面动态沉陷预测与复垦方法 |
-
2016
- 2016-03-03 CN CN201610121402.6A patent/CN105808935B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101422094A (zh) * | 2007-11-02 | 2009-05-06 | 煤炭科学研究总院唐山研究院 | 平原矿区厚煤层开采塌陷地动态预复垦方法 |
CN101709645A (zh) * | 2009-12-17 | 2010-05-19 | 煤炭科学研究总院唐山研究院 | 一种矿区厚煤层开采沉陷区就地取土复垦方法 |
CN104564069A (zh) * | 2015-01-09 | 2015-04-29 | 中国矿业大学(北京) | 一种基于方格网法的地面动态沉陷预测与复垦方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257870A (zh) * | 2020-02-26 | 2020-06-09 | 安徽大学 | 一种利用InSAR监测数据的采煤沉陷积水区水下地形反演方法 |
CN112749477A (zh) * | 2020-11-30 | 2021-05-04 | 安徽理工大学 | 厚表土薄基岩煤矿立井井筒偏斜量的预计方法 |
CN112749477B (zh) * | 2020-11-30 | 2021-09-24 | 安徽理工大学 | 厚表土薄基岩煤矿立井井筒偏斜量的预计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105808935B (zh) | 2018-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
He et al. | The synthetic geo-ecological environmental evaluation of a coastal coal-mining city using spatiotemporal big data: a case study in Longkou, China | |
Dawei et al. | Formation and development mechanism of ground crack caused by coal mining: effects of overlying key strata | |
Hu et al. | Optimization of concurrent mining and reclamation plans for single coal seam: a case study in northern Anhui, China | |
Keim et al. | A numerical study on optimization of multilateral horizontal wellbore patterns for coalbed methane production in Southern Shanxi Province, China | |
Xiao et al. | Dynamic subsidence simulation and topsoil removal strategy in high groundwater table and underground coal mining area: a case study in Shandong Province | |
Jiang et al. | Landslide risk prediction by using GBRT algorithm: Application of artificial intelligence in disaster prevention of energy mining | |
Li et al. | Innovation for sustainable mining: Integrated planning of underground coal mining and mine reclamation | |
CN105114068A (zh) | 一种利用测井资料预测煤层气区块高产水区的方法 | |
Chi et al. | Simulation analysis of water resource damage feature and development degree of mining-induced fracture at ecologically fragile mining area | |
Liu et al. | A boundary model of terrain reconstruction in a coal-mining subsidence waterlogged area | |
CN111798052B (zh) | 一种高潜水位采煤沉陷积水区三维空间信息动态预测方法 | |
Li et al. | The optimal framework and model to balance underground coal mining and cropland protection in Jining, eastern China | |
CN103541376B (zh) | 采煤沉陷区地基在重复开采条件下的基础变形预测方法 | |
Feng et al. | Improving mine reclamation efficiency for farmland sustainable use: Insights from optimizing mining scheme | |
CN102877846B (zh) | 基于单元法的边采边复表土剥离时机确定方法 | |
Qiao et al. | Prediction of floor water disasters based on fractal analysis of geologic structure and vulnerability index method for deep coal mining in the Yanzhou mining area | |
Dai et al. | A sensitivity analysis of factors affecting in geologic CO2 storage in the Ordos Basin and its contribution to carbon neutrality | |
Xue et al. | A prediction model for overlying rock thickness of subsea tunnel: A hybrid intelligent system | |
CN105808935A (zh) | 多煤层开采下地面积水时间的确定方法 | |
CN105808934B (zh) | 多煤层开采下地表受损范围及扰动次数确定方法 | |
Popovs et al. | Bedrock surface topography of Latvia | |
CN103742151B (zh) | 应用于边采边复技术的土壤重构方法 | |
Li et al. | Enhancing stability analysis of open-pit slopes via integrated 3D numerical modeling and data monitoring | |
He et al. | Water Inrush Mechanism and Treatment Measures in Huali Highway Banyanzi Tunnel—A Case Study | |
CN110675272A (zh) | 一种煤炭超大工作面地下开采减损关键参数的确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |