CN105808935B - 多煤层开采下地面积水时间的确定方法 - Google Patents

多煤层开采下地面积水时间的确定方法 Download PDF

Info

Publication number
CN105808935B
CN105808935B CN201610121402.6A CN201610121402A CN105808935B CN 105808935 B CN105808935 B CN 105808935B CN 201610121402 A CN201610121402 A CN 201610121402A CN 105808935 B CN105808935 B CN 105808935B
Authority
CN
China
Prior art keywords
mining
seam
coal seam
coal
mutil
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
Application number
CN201610121402.6A
Other languages
English (en)
Other versions
CN105808935A (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.)
Pingan Coal Mining Exploitation Engineering Technology Institute Co ltd
China University of Mining and Technology Beijing CUMTB
Original Assignee
Pingan Coal Mining Exploitation Engineering Technology Institute Co ltd
China University of Mining and Technology Beijing CUMTB
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 Pingan Coal Mining Exploitation Engineering Technology Institute Co ltd, China University of Mining and Technology Beijing CUMTB filed Critical Pingan Coal Mining Exploitation Engineering Technology Institute Co ltd
Priority to CN201610121402.6A priority Critical patent/CN105808935B/zh
Publication of CN105808935A publication Critical patent/CN105808935A/zh
Application granted granted Critical
Publication of CN105808935B publication Critical patent/CN105808935B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

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。
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 (1)

1.一种多煤层开采下地面积水时间的确定方法，其特征在于，该方法包括：根据地下煤炭赋存条件，利用概率积分法分别预计各煤层开采引起的地面下沉，依据煤层开采计划，通过叠加确定多煤层重复采动下地表累积下沉量，然后将其与地表原始高程进行空间叠加，模拟多煤层开采后地面沉陷情况，最后通过比较不同煤层开采结束后地面高程和当地潜水位埋深的关系，确定多煤层开采地面积水时间；
该方法具体包括以下步骤：
1)根据地下煤炭赋存条件，利用概率积分法分别预计各煤层开采引起的地面下沉：设开采区域共n个煤层，且开采顺序为煤层1、煤层2、…、煤层i、…、煤层n，地面坐标系为x、o、y，煤层i投影为水平煤层后的坐标系为εi、o、ηi，且坐标系的水平投影重合，地面任意点P的坐标为(x,y)，则煤层i开采区域Si，引起点P的下沉值Wi(x,y)如式(1)所示：
式(1)中：Wi(x,y)为煤层i开采引起地面任意点P的下沉值，单位mm；
qi为煤层i开采的下沉系数；mi为煤层i的煤厚，单位mm；
αi为煤层i的煤层倾角，单位°；
ri为煤层i的主要影响半径，单位m；i＝1,2,…,n；
2)根据煤层开采计划，确定多煤层重复采动下地表累积下沉量：m个煤层重复采动引起的地表累积下沉如式(2)所示：
式(2)中：为m个煤层重复采动引起地面任意点P的累积下沉值，单位m；
Wj(x,y)为煤层j开采引起地面任意点P的下沉值，单位mm；
m＝1,2,…,n；j＝1,2,…,m；
3)结合地表原始高程，模拟多煤层开采地面沉陷：将多煤层开采地表累积下沉与煤炭开采前地表的原始高程进行空间叠加分析，模拟多煤层开采后地面沉陷情况，m个煤层开采结束后，地面任意点P的高程如式(3)所示：
式(3)中：Hm(x,y)为m个煤层开采结束后地面任意点P的高程，单位m；
H0(x,y)为煤炭开采前地面任意点P的原始高程，单位m；
m＝1,2,…,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。
CN201610121402.6A 2016-03-03 2016-03-03 多煤层开采下地面积水时间的确定方法 Active CN105808935B (zh)

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 CN105808935A (zh) 2016-07-27
CN105808935B true 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)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257870B (zh) * 2020-02-26 2022-12-06 安徽大学 一种利用InSAR监测数据的采煤沉陷积水区水下地形反演方法
CN112749477B (zh) * 2020-11-30 2021-09-24 安徽理工大学 厚表土薄基岩煤矿立井井筒偏斜量的预计方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101422094B (zh) * 2007-11-02 2012-07-04 煤炭科学研究总院唐山研究院 平原矿区厚煤层开采塌陷地动态预复垦方法
CN101709645A (zh) * 2009-12-17 2010-05-19 煤炭科学研究总院唐山研究院 一种矿区厚煤层开采沉陷区就地取土复垦方法
CN104564069B (zh) * 2015-01-09 2016-08-17 中国矿业大学(北京) 一种基于方格网法的地面动态沉陷预测与复垦方法

Also Published As

Publication number Publication date
CN105808935A (zh) 2016-07-27

Similar Documents

Publication Publication Date Title
Hu et al. Optimization of concurrent mining and reclamation plans for single coal seam: a case study in northern Anhui, China
Chi et al. Simulation analysis of water resource damage feature and development degree of mining-induced fracture at ecologically fragile mining area
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
CN111750822A (zh) 一种采煤诱发的覆岩与地表沉陷协同动态预测方法
CN104453903A (zh) 一种近距煤层群保水开采方法
Shi Mine green mining
CN103541376B (zh) 采煤沉陷区地基在重复开采条件下的基础变形预测方法
CN103235853B (zh) 用于采煤沉陷地边采边复的基于土方平衡的基塘布局方法
CN107101617A (zh) 一种基于遥感监测的冲沟地形下采场覆岩活动控制方法
Yuan et al. Application and prospect of curtain grouting technology in mine water safety management in China: A review
CN105808935B (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
CN105808934B (zh) 多煤层开采下地表受损范围及扰动次数确定方法
Li et al. A new approach to increased land reclamation rate in a coal mining subsidence area: A case‐study of Guqiao Coal Mine, China
CN104564069A (zh) 一种基于方格网法的地面动态沉陷预测与复垦方法
Zhao et al. Discrete element analysis of deformation features of slope controlled by karst fissures under the mining effect: a case study of Pusa landslide, China
CN111400659A (zh) 采-选-充一体化开采技术耕地保护率计算方法
CN110675272A (zh) 一种煤炭超大工作面地下开采减损关键参数的确定方法
Xie et al. Stability analysis of seabed strata and casing structure during the natural gas hydrates exploitation by depressurization in horizontal wells in South China Sea
Tianwei et al. Numerical study on 4-1 coal seam of xiaoming mine in ascending mining
Chunhu et al. Water inflow characteristics of coal seam mining aquifer in Yushen mining area, China
Zeng Research on risk evaluation methods of groundwater bursting from aquifers underlying coal seams and applications to coalfields of North China
Robins et al. Groundwater flow in the South Wales coalfield: historical data informing 3D modelling
Yao et al. Sliding surface searching method for slopes containing a potential weak structural surface
CN105045962B (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