CN114817826A - 一种基于等效下沉的采动覆岩离层计算方法 - Google Patents

一种基于等效下沉的采动覆岩离层计算方法 Download PDF

Info

Publication number
CN114817826A
CN114817826A CN202210443353.3A CN202210443353A CN114817826A CN 114817826 A CN114817826 A CN 114817826A CN 202210443353 A CN202210443353 A CN 202210443353A CN 114817826 A CN114817826 A CN 114817826A
Authority
CN
China
Prior art keywords
rock stratum
mining
layer
separation
stratum
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.)
Pending
Application number
CN202210443353.3A
Other languages
English (en)
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.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
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 China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN202210443353.3A priority Critical patent/CN114817826A/zh
Publication of CN114817826A publication Critical patent/CN114817826A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • 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
    • E21F15/00—Methods or devices for placing filling-up materials in underground workings
    • 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
    • E21F15/00—Methods or devices for placing filling-up materials in underground workings
    • E21F15/005—Methods or devices for placing filling-up materials in underground workings characterised by the kind or composition of the backfilling material

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geology (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明提出一种基于等效下沉的采动覆岩离层计算方法，属于矿业工程领域中矿山开采沉陷预计、覆岩离层注浆减沉与岩层控制领域。该方法首先计算开采工作面每个岩层的采动下沉值；从表土层下第一个岩层开始，将每个岩层的采动下沉值与该岩层上位岩层的采动下沉值作差，得到每个岩层的覆岩离层量；按照所述离层量从大到小的顺序对各岩层进行排序，将序列中所述离层量的最大值对应的岩层作为当前岩层；对所述当前岩层进行判定，以确定覆岩离层注浆的岩层位置。本发明可方便地计算出上覆岩层的离层量和离层层位，为覆岩离层注浆减沉设计和实施提供简洁的判断依据。

Description

一种基于等效下沉的采动覆岩离层计算方法
技术领域
本发明属于矿业工程领域中矿山开采沉陷预计、覆岩离层注浆减沉与岩层控制领域，具体地涉及一种基于等效下沉的采动覆岩离层计算方法。
背景技术
充填开采是控制覆岩与地表沉陷、减小土地及生态环境破坏的有效手段，是实现矿区绿色可持续发展的重要途径。充填开采包括采空区充填和覆岩离层注浆充填两种类型。采空区充填又可细分为水砂充填、矸石自溜充填、高水材料充填、固体密实充填等，尽管采空区充填可以达到很好的控制覆岩与地表移动效果，但也存在充填系统与生产系统并行、工艺复杂、成本高等缺点。而覆岩离层注浆则是通过地面钻孔作为充填通道，以电厂粉煤灰、矿渣等为充填材料，该充填系统与生产系统分离，互不干扰，当注浆层位选择准确、注浆时机把握及时，也会取得良好的覆岩移动控制效果。覆岩离层注浆的难点在于离层层位和离层量的确定。目前的确定方法大多利用岩层控制的关键层理论，通过复杂的力学计算，确定离层层位，但还难以定量确定离层量；此外也有通过有限元、边界元、离散元、快速拉格朗日分析等数值模拟方法进行计算判定，由于边界条件的限制、物理力学参数选择的不确定性，往往需要有一定的工程经验、反复调试，模拟计算结果不尽人意。
发明内容
本发明的目的是为克服已有技术的不足之处，提出一种基于等效下沉的采动覆岩离层计算方法。本发明可方便地计算出上覆岩层的离层量和离层层位，为覆岩离层注浆减沉设计和实施提供简洁的判断依据。
本发明是实力提出一种基于等效下沉的采动覆岩离层计算方法，包括：
计算开采工作面每个岩层的采动下沉值；
从表土层下第一个岩层开始，将每个岩层的采动下沉值与该岩层上位岩层的采动下沉值作差，得到每个岩层的覆岩离层量；
按照所述离层量从大到小的顺序对各岩层进行排序，将序列中所述离层量的最大值对应的岩层作为当前岩层；
对所述当前岩层进行判定，以确定覆岩离层注浆的岩层位置；
其中，对于从地表向下的第j层岩层，所述采动下沉值计算表达式如下：
Figure BDA0003615492180000021
式中，xo为开采工作面地表设定的参照点po的坐标，点po的下沉值为wo(xo)；点po在向采空区移动方向上对应的第j层岩层的覆岩点
Figure BDA0003615492180000022
坐标为
Figure BDA0003615492180000023
覆岩点
Figure BDA0003615492180000024
的下沉值为
Figure BDA0003615492180000025
b为开采工作面宽度；δ0为覆岩与地表移动的综合边界角；hi为从地表开始向下的第i层的岩层厚度，i＝1,2,…,n，n为岩层总数；当i＝1时，h1代表表土层厚度；
Figure BDA0003615492180000026
为开采深度；
Figure BDA0003615492180000027
表示从地表向下至第j层岩层底界面的累积深度，j＝0,1,2,3,…,n；当j＝0时，h0＝0。
在本发明的一个具体实施例中，对所述当前岩层进行判定，以确定覆岩离层注浆的岩层位置，包括：
若所述当前岩层的底界面距离开采层顶板底部的高度小于导水裂缝带的高度值，则所述当前岩层的离层为无效离层，从所述序列中剔除所述当前岩层；从更新后的序列中选出离层量最大值对应的岩层作为新的当前岩层，然后重新对所述当前岩层进行判定；
若所述当前岩层的底界面距离开采层顶板底部的高度大于等于所述导水裂缝带的高度值，则所述当前岩层的离层为有效离层，将所述离层作为覆岩离层注浆层位；
若所述序列遍历完毕后，所有岩层的离层均为无效离层，则所述开采工作面无法进行覆岩离层注浆。
在本发明的一个具体实施例中，所述导水裂缝带的高度值计算表达式如下：
Figure BDA0003615492180000028
式中，H导为导水裂缝带高度，∑M为累积开采厚度；
将H导的上限值作为最终的导水裂缝带的高度值。
本发明的特点及有益效果在于：
本发明提出的一种基于等效下沉的采动覆岩离层计算方法，该方法计算涉及的参数少，仅与工作面宽度、开采深度、岩层厚度、岩层埋藏深度、覆岩与地表移动的综合边界角有关，不涉及各岩层的物理力学特性，且能反映采动程度的影响；计算过程简便。可以方便地计算出各岩层的下沉，利用相邻岩层的非一致性下沉，可以得到离层量；通过剔除无效离层，可以得到有效的封闭离层和层位，为离层注浆减沉提供技术判据。
附图说明
图1为本发明实施例中一种基于等效下沉的采动覆岩离层计算方法的整体流程图。
具体实施方式
本发明提出一种基于等效下沉的采动覆岩离层计算方法，下面结合具体实施例对本发明做进一步详细说明。
本发明提出一种基于等效下沉的采动覆岩离层计算方法，整体流程如图1所示，包括以下步骤：
1)根据开采工作面的地质采矿条件，计算各岩层的采动下沉值；
其中，对于从地表向下的第j层岩层，其下沉值计算表达式如下：
Figure BDA0003615492180000031
式中，以地表下沉盆地中心为坐标原点，开采工作面地表任一点po的坐标为xo，将点po作为参照点，点po下沉值记为wo(xo)；地表点po在向采空区移动方向上对应的第j层岩层的覆岩点
Figure BDA0003615492180000032
的坐标为
Figure BDA0003615492180000033
覆岩点
Figure BDA0003615492180000034
的下沉值记为
Figure BDA0003615492180000035
需要说明的是，本发明计算时每个岩层都是参照某个地表点po进行计算，但不同埋深岩层上对应的覆岩点坐标不同，即同一个地表点po，对于n个岩层，则有n个ps点，每个岩层的ps点坐标不同。
b为开采工作面宽度，单位为米；δ0为覆岩与地表移动的综合边界角，单位为度；hi为从地表开始向下的第i层的岩层厚度，i＝1,2,…,n，n为岩层总数；当i＝1时，h1代表表土层厚度，单位为米；
Figure BDA0003615492180000036
为开采深度，单位为米,为一定值；
Figure BDA0003615492180000037
表示从地表向下至第j层岩层底界面的累积深度(该深度也是从地表向下至第j+1层岩层顶界面的累积深度)，j＝0,1,2,3,…,n；当j＝0时，h0＝0，代表的计算深度是地表；当j＝1时，
Figure BDA0003615492180000038
Figure BDA0003615492180000041
代表的计算深度是第1层岩层即表土层的底界面(也是第2层岩层的顶界面)累积深度，单位为米；当j＝n时，
Figure BDA0003615492180000042
等于开采深度，单位为米，代表从地表向下至第n层岩层底部界面(即开采层的顶板底部，本实施例为开采煤层)的累积深度。
2)利用步骤1)的结果，计算各岩层的覆岩离层量。
由于同一岩层具有整体性，下沉相同；而不同层位岩层的下沉则具有不一致性。从表土层下第一个岩层开始，将该岩层的下沉值减去该岩层的上位岩层的下沉值，即为该岩层的覆岩离层量，其中上位岩层为位于该岩层上方的相邻岩层。
3)按照离层量从大到小的顺序对各岩层进行排序，将序列中离层量最大值对应的岩层作为当前岩层；
4)通过对当前岩层进行判定，确定覆岩离层注浆的岩层位置。
具体方法为：
4-1)按《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规范》中的经验公式计算导水裂缝带高度，对于坚硬覆岩表达式如下：
Figure BDA0003615492180000043
式中，H导为导水裂缝带高度，单位为米；∑M为累积开采厚度，单位为米。
4-2)判定：
若当前岩层的底界面距离开采层顶板底部的高度小于导水裂缝带高度值(本实施例采用利用式(2)计算得到的上限值作为导水裂缝带高度值)，则当前岩层的离层为无效离层，从步骤3)得到的序列中剔除该岩层，从更新后的序列中选出离层量最大值对应的岩层作为新的当前岩层，然后重新返回步骤3-4)；
若当前岩层的底界面距离开采层顶板底部的高度大于等于导水裂缝带高度，则当前岩层的离层为有效离层，将该离层作为覆岩离层注浆层位。
若序列遍历完毕后，所有当前岩层的离层均为无效离层，则该工作面无法进行覆岩离层注浆。
下面结合一个具体实施例对本发明方法进一步详细说明如下。
本发明一个具体实施例中，开采工作面宽80米，走向达到充分采动，开采深度200米，宽深比为0.4；开采厚度2米，覆岩与地表移动的综合边界角δ0＝65°。上覆岩层的岩性自地表到顶板如表1第2列所示，以砂岩、粉砂岩为主，单向抗拉强度为55MPa，属坚硬岩层；各岩层厚度如表1第3列所示，累积深度如表1第4列所示。
本实施例提出一种基于等效下沉的采动覆岩离层计算方法，包括以下步骤：
1)根据开采工作面的地质采矿条件，计算各岩层的采动下沉值；
本实施例中，煤层开采后，采动下沉自顶板向上覆岩层传递，顶板以垮落形式冒落堆积于采空区，顶板下沉为2000毫米，等于煤层开采厚度。
利用式(1)计算各岩层的下沉值，结果见表1中的第5列。
2)利用步骤1)的结果，计算各岩层的覆岩离层量。
由于不同层位层的下沉不一致性，从表土层下第一个岩层开始，下位岩层的下沉值减去其上位岩层的下沉值，即为下位岩层的覆岩离层量，本实施例中各岩层的覆岩离层量结果见表1中的第6列。
3)按照离层量从大到小的顺序对各岩层进行排序，将序列中离层量最大值对应的下位岩层作为当前岩层；
由表1第6列可知，经过排序后，离层量最大值对应第39层页岩层，其上位岩层为第38层粉砂岩层，第38层粉砂岩层与第39层页岩层间将产生138毫米的离层；离层量第二最大值对应第26层泥岩层，其上位岩层为第25层砂岩层，第25层砂岩层与第26层泥岩层间将产生125毫米的离层。
4)通过对当前岩层进行判定，确定覆岩离层注浆的岩层位置。
具体方法为：
4-1)按照式(2)计算导水裂缝带高度；本实施例中，计算可得导水裂缝带高度为45.45±8.9米，则导水裂缝带高度的上限值为54.35米。
4-2)判定：
根据表1，将第39层页岩层作为当前岩层时，第39层页岩层的底界面距离煤层顶板底部高度为200-182＝18米，小于导水裂缝带高度的上限值，表明第39层页岩层位于该裂缝带内，第38层粉砂岩层与第39层页岩层间无法形成封闭的充填空间，故不予考虑、应剔除。
将第26层泥岩层作为当前岩层时，第26层泥岩层的底界面距离煤层顶板底部高度为200-134＝66米，大于导水裂缝带高度的上限值，可选择第25层砂岩层和第26层泥岩层间的125毫米离层作为覆岩离层注浆层位。
表1本发明一个具体实施例中采动覆岩离层的参数表
Figure BDA0003615492180000051
Figure BDA0003615492180000061

Claims (3)

1.一种基于等效下沉的采动覆岩离层计算方法，其特征在于，包括：
计算开采工作面每个岩层的采动下沉值；
从表土层下第一个岩层开始，将每个岩层的采动下沉值与该岩层上位岩层的采动下沉值作差，得到每个岩层的覆岩离层量；
按照所述离层量从大到小的顺序对各岩层进行排序，将序列中所述离层量的最大值对应的岩层作为当前岩层；
对所述当前岩层进行判定，以确定覆岩离层注浆的岩层位置；
其中，对于从地表向下的第j层岩层，所述采动下沉值计算表达式如下：
Figure FDA0003615492170000011
式中，xo为开采工作面地表设定的参照点po的坐标，点po的下沉值为wo(xo)；点po在向采空区移动方向上对应的第j层岩层的覆岩点
Figure FDA0003615492170000012
坐标为
Figure FDA0003615492170000013
覆岩点
Figure FDA0003615492170000014
的下沉值为
Figure FDA0003615492170000015
b为开采工作面宽度；δ0为覆岩与地表移动的综合边界角；hi为从地表开始向下的第i层的岩层厚度，i＝1,2,…,n，n为岩层总数；当i＝1时，h1代表表土层厚度；
Figure FDA0003615492170000016
为开采深度；
Figure FDA0003615492170000017
表示从地表向下至第j层岩层底界面的累积深度，j＝0,1,2,3,…,n；当j＝0时，h0＝0。
2.根据权利要求1所述的方法，其特征在于，对所述当前岩层进行判定，以确定覆岩离层注浆的岩层位置，包括：
若所述当前岩层的底界面距离开采层顶板底部的高度小于导水裂缝带的高度值，则所述当前岩层的离层为无效离层，从所述序列中剔除所述当前岩层；从更新后的序列中选出离层量最大值对应的岩层作为新的当前岩层，然后重新对所述当前岩层进行判定；
若所述当前岩层的底界面距离开采层顶板底部的高度大于等于所述导水裂缝带的高度值，则所述当前岩层的离层为有效离层，将所述离层作为覆岩离层注浆层位；
若所述序列遍历完毕后，所有岩层的离层均为无效离层，则所述开采工作面无法进行覆岩离层注浆。
3.根据权利要求2所述的方法，其特征在于，所述导水裂缝带的高度值计算表达式如下：
Figure FDA0003615492170000021
式中，H导为导水裂缝带高度，∑M为累积开采厚度；
将H导的上限值作为最终的导水裂缝带的高度值。
CN202210443353.3A 2022-04-26 2022-04-26 一种基于等效下沉的采动覆岩离层计算方法 Pending CN114817826A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210443353.3A CN114817826A (zh) 2022-04-26 2022-04-26 一种基于等效下沉的采动覆岩离层计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210443353.3A CN114817826A (zh) 2022-04-26 2022-04-26 一种基于等效下沉的采动覆岩离层计算方法

Publications (1)

Publication Number Publication Date
CN114817826A true CN114817826A (zh) 2022-07-29

Family

ID=82507449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210443353.3A Pending CN114817826A (zh) 2022-04-26 2022-04-26 一种基于等效下沉的采动覆岩离层计算方法

Country Status (1)

Country Link
CN (1) CN114817826A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116644614A (zh) * 2023-06-30 2023-08-25 中国矿业大学 一种采动覆岩突水离层带层位识别方法
CN117370727A (zh) * 2023-10-24 2024-01-09 山东能源集团有限公司 覆岩离层注浆减沉降震评价方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116644614A (zh) * 2023-06-30 2023-08-25 中国矿业大学 一种采动覆岩突水离层带层位识别方法
CN116644614B (zh) * 2023-06-30 2023-11-14 中国矿业大学 一种采动覆岩突水离层带层位识别方法
CN117370727A (zh) * 2023-10-24 2024-01-09 山东能源集团有限公司 覆岩离层注浆减沉降震评价方法

Similar Documents

Publication Publication Date Title
CN103902780B (zh) 固体充填采煤地表变形预计方法
CN100467829C (zh) 利用结构关键层作为隔水层的保水采煤方法
Hoek Putting numbers to geology—an engineer's viewpoint
CN114817826A (zh) 一种基于等效下沉的采动覆岩离层计算方法
WO2016041389A1 (zh) 一种固体充填采煤工作面充填质量评价方法
CN104408277B (zh) 矿区新建建筑物引起的地表残余移动变形预测及防治方法
CN106150546B (zh) 一种根据地表沉陷特征进行采空区部分充填的方法
CN108663724B (zh) 一种基于水资源转移路径的煤矿地下水库位置选择方法
Zhang et al. A case study on integrated modeling of spatial information of a complex geological body
Ulusay et al. Improvement of slope stability based on integrated geotechnical evaluations and hydrogeological conceptualisation at a lignite open pit
Sun et al. A study of surface subsidence and coal pillar safety for strip mining in a deep mine
CN110851991B (zh) 一种地下水流数值模拟方法
CN117167080A (zh) 一种三软厚煤层开采覆岩多层位离层注浆减沉方法
CN115199273A (zh) 基于地表沉陷控制要求的采煤-充填工作面布局方法
CN111159643B (zh) 一种基于实测数据的Weibull时间函数对地面沉陷影响的方法
Guang-Li et al. Similar material and numerical simulation of strata movement laws with long wall fully mechanized gangue backfilling
CN114370250B (zh) 基于地层勘探的采空区灾害治理方法
Kortas et al. Using the GRAM model to reconstruct the important factors in historic groundwater rebound in part of the Durham Coalfield, UK
Zhang et al. Study on monitoring and variation law of strata movement induced by caving mining of slowly inclined large and thick orebody
Lu et al. Monitoring and analysis of ground subsidence and backfill stress distribution in Jinchuan Mine, China
BELL The history and techniques of coal mining and the associated effects and influence on construction
Romero Benitez et al. Prediction of ground movements under dipping seam conditions
Fan et al. The characteristics of grout diffusion in reduce subsidence mining with abscission layer grouting
Sheets et al. Remediation of large-scale slope
Li et al. Rock strata failure and subsidence characteristics under the mining of short distance thick coal seams: a case in west China

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