CN114519257A - 一种地表沉陷类别预计方法 - Google Patents

一种地表沉陷类别预计方法 Download PDF

Info

Publication number
CN114519257A
CN114519257A CN202210073917.9A CN202210073917A CN114519257A CN 114519257 A CN114519257 A CN 114519257A CN 202210073917 A CN202210073917 A CN 202210073917A CN 114519257 A CN114519257 A CN 114519257A
Authority
CN
China
Prior art keywords
layer
key layer
arch
loose
key
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.)
Withdrawn
Application number
CN202210073917.9A
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.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
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 Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN202210073917.9A priority Critical patent/CN114519257A/zh
Publication of CN114519257A publication Critical patent/CN114519257A/zh
Priority to LU503279A priority patent/LU503279B1/en
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Primary Health Care (AREA)
  • Game Theory and Decision Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • Operations Research (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Agronomy & Crop Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Mining & Mineral Resources (AREA)
  • Animal Husbandry (AREA)
  • Development Economics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明提供一种地表沉陷类别预计方法，包括：收集钻孔周围松散层的力学参数和各个岩层结构的力学参数；根据松散层拱理论，判定松散层拱是否可以形成以及形成后的稳定性；通过对上覆岩层发挥的控制作用判别关键层的位置；当松散层中无法形成松散层拱结构时，判断关键层为远场关键层还是近场关键层；当关键层破断形成砌体梁时且Δ＞ΔT时，则砌体梁会发生滑落失稳；再根据近场关键层和远场关键层以及砌体梁稳定性和松散层拱稳定性判定地表均匀下沉盆地、地表倾向台阶裂缝、地表走向斑裂、地表塌陷坑。本发明根据得出的地表沉陷类型针对不同情况提出防治措施，对矿山开采中实现经济发展、财产安全、保护环境等方面的协调发展有着一定的指导意义。

Description

一种地表沉陷类别预计方法
技术领域
本发明涉及矿山开采地表沉陷技术领域，尤其涉及一种地表沉陷 类别预计方法。
背景技术
煤矿开采过程中，上覆岩层内部的移动变形是引起地表沉陷的主要因素， 地表沉陷是开采过程中覆岩运动由下往上传递到地表的最终反映。因煤炭开采 产生的地表沉陷，使地层稳定受损，威胁到地表的安全，并破坏了生态环境， 限制了经济发展，危害了社会和谐，让城市近郊沉陷地成为社会经济发展、城 市扩展以及基础建设的用地瓶颈。
已有研究表明，煤炭开采后采场上覆岩层运动受控于采动覆岩承 载结构，松散层中的承载结构为松散层拱，基岩中的承载结构为关键 层。当关键层初次破断时，地表的最大下沉值和最大下沉速度尚未达 到最大值，随着松散层拱周期性失稳逐渐增大；当松散层拱完全失稳 破坏后，地表最大下沉速度逐渐减小。然而我国地质构造条件较为复 杂，地表沉陷种类繁多，如：地表均匀下沉盆地、地表倾向台阶裂缝、 地表走向斑裂和地表塌陷坑等。现有研究地表沉陷的方法有沉陷预计 和地表沉陷观测等方法，沉陷预计运用概率积分法，就地表将出现的 移动和变形规律进行预计，包括地表下沉值，地表下沉范围及地表下沉占总面积的百分比，但无法对地表沉陷特征进行超前判别；而地表 沉陷观测需在开采区域上方地表设置地表岩移观测站定量监测地面 的移动变形，获得地表移动变形参数，最大下沉点以及地面的累计下 沉量，对地表沉陷进行实时监控，但无法对地表沉陷进行超前预计。
发明内容
本发明的目的在于解决上述现有技术存在的缺陷，提供一种地表 沉陷类别预计方法。
一种地表沉陷类别预计方法，包括以下步骤：
步骤一，收集钻孔周围松散层的力学参数和各个岩层结构的力学 参数；
步骤二，根据松散层拱理论，判定松散层拱是否可以形成以及形 成后的稳定性；
步骤三，根据关键层理论，通过对上覆岩层发挥的控制作用判别 关键层的位置；
步骤四，当松散层中无法形成松散层拱结构时，根据关键层理论， 判断关键层为远场关键层还是近场关键层；当关键层破断形成砌体梁 时，进入步骤五；
步骤五，当关键层破断形成砌体梁时，且Δ＞ΔT时，则砌体梁会 发生滑落失稳；其中，Δ为工作面回采后关键层与垮落的直接顶的间 距，ΔT为砌体梁保持自身稳定的临界回转量；
步骤六，当上覆岩层中关键层和松散层拱同时存在，且上覆岩层 中关键层和松散层拱保持稳定时，则地表会产生均匀下沉盆地；
当上覆岩层中只有近场关键层而无松散层拱，且上覆岩层中近场 关键层破断后，则地表会产生倾向台阶裂缝；
当上覆岩层中只存在远场关键层而无松散层拱，且上覆岩层中远 场关键层破断后，则地表会产生走向斑裂；
当关键层破断后的砌体梁和松散层拱同时失稳时，则地表会产生 塌陷坑。
进一步地，如上所述的地表沉陷类别预计方法，其特征在于，步 骤一中，所述松散层的力学参数包括：松散层的容重、厚度、矢高、 内聚力、内摩擦角、侧压系数、基岩破断角、煤层顶界面与关键层底 界面间的距离；
所述各个岩层的力学参数包括：各个岩层的容重、关键层的厚度、 破断距、关键层的抗拉强度、所受载荷、直接顶的碎胀系数、直接顶 的垮落高度，煤层的厚度和埋深、基岩中各个岩层厚度及岩性。
进一步地，如上所述的地表沉陷类别预计方法，在煤炭开采区域 地表施工地面钻孔取芯保存，并绘制钻孔综合柱状图；根据钻孔综合 柱状图，确定煤炭开采区域开采煤层的厚度和埋深、松散层厚度和基 岩中各岩层厚度及岩性；
通过对松散层进行取芯，获得的试样进行物理力学参数测试，从 而得到松散层容重、内摩擦角和内聚力；
通过对基岩中各个岩层进行取芯，获得的试样进行物理力学参数 测试，从而得到基岩中各个岩层的容重、弹性模量和各个岩层的抗拉 强度；直接顶的碎胀系数由基岩中岩石力学性质和岩层的破坏形式确 定；
在煤层开采区域，通过井下水压致裂地应力测量装置，确定松散 层的侧压系数；
利用关键层判别软件KSPB，基于松散层和基岩力学的特性，在 钻孔综合柱状图进行关键层的判别，从而得到关键层的厚度、位置、 破断距、所受载荷、关键层的抗拉强度、关键层底界面与煤层顶界面 间的距离；
利用钻孔测井分析仪，在工作面巷道顶板上方打钻孔，通过窥视 结果中钻孔的深度，绘制出窥视钻孔破断演化图，得出基岩破断角。
进一步地，如上所述的地表沉陷类别预计方法，所述步骤二中， 以松散层拱左侧拱基为原点，拱基所在直线为X轴，X轴向右为正 方向，垂直于拱基所在直线的直线为Y轴，Y轴向上为正方向建立 松散层拱二维力学模型；
根据该模型计算使松散层拱形成及稳定的松散层临界厚度：
Figure BDA0003483244450000051
式中：Larch、Harch、δarch为松散层拱的跨度、矢高和厚度；γ、α 为基岩容重和破断角；C、
Figure BDA0003483244450000052
为松散层的内聚力和内摩擦角；λ为侧压 系数；Σh为主关键层底界面与煤层顶界面间的距离；h0为松散层拱 上覆松散层厚度；HC为松散层临界厚度；
根据松散层拱的形成及稳定条件：当Σh＞HC时，松散层中能形 成松散层拱；当松散层拱稳定时，据此推出煤矿开采参数中的工作面 临界宽度LM。
进一步地，如上所述的地表沉陷类别预计方法，所述步骤三中， 由下式判别关键层的强度准则和刚度准则：
Figure BDA0003483244450000053
式中：q1|n+1、q1|n为第n+1和n层岩层对第1层岩层的载荷；ln+1、 ln为第n+1和n层岩层的破断距；
借助关键层判别位置软件KSPB对工作面的关键层进行判别，推 断出关键层在柱状图中的位置，从而判定关键层在工作面上覆岩层中 的位置。
进一步地，如上所述的地表沉陷类别预计方法，所述步骤四中， 当开采区域地质条件确定存在关键层，无法形成松散层拱时，根据远 近场关键层划分依据，将工作面上覆关键层划分为远场关键层和近场 关键层；
由下式确定a和b：
Figure BDA0003483244450000061
式中，a为工作面长度；b为关键层倾向悬露长度；lm关键层极 限跨距；
根据关键层破断规律：3lm>b>1.414lm，来压步距lm<a<1.414lm， 此时关键层“竖O-X”破断，关键层为近场关键层；lm<b<1.414lm，来 压步距a>1.414lm>b，此时关键层“横O-X”破断，关键层为远场关键 层。据此推出工作面上覆关键层为近场关键层或远场关键层。
进一步地，如上所述的地表沉陷类别预计方法，所述步骤五中， 由下式确定Δ和ΔT：
Figure BDA0003483244450000071
式中：MB为工作面临界采高；KP为直接顶的碎胀系数；Σhi为直 接顶的垮落高度；σC为关键层的抗压强度；hKS、qKS、lKS为关键层的 厚度、载荷和破断距。
有益效果：
本发明为了预计地表沉陷类型，确定了松散层拱的形成条件并进 行稳定性判别，确定了关键层稳定性判别及临界采高，确定了关键层 为远场关键层或近场关键层。在不同地质条件下，采动覆岩承载结构 的组成存在多样性，故需要一种通过采动覆岩承载结构特征预计地表 沉陷类型，并且根据得出的地表沉陷类型针对不同情况提出防治措施， 预测地表沉陷各项参数，减少或杜绝开采损害带来的损失，对矿山开 采中实现经济发展、财产安全、保护环境等方面的协调发展有着一定 的指导意义。
附图说明
图1为本发明地表沉陷类别预计方法流程图；
图2为本发明地表沉陷类别预计方法原理图；
图3为松散层拱结构二维力学模型图；
图4为地表下沉曲线图；
图5为地表均匀下沉盆地现场图；
图6为钻孔柱状图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面本发明 中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实 施例，本领域普通技术人员在没有作出创造性劳动前提下所获得 的所有其他实施例，都属于本发明保护的范围。
如图1、图2所示，本发明提供的方法，包括以下步骤：
步骤一，收集钻孔周围松散层的力学参数和各个岩层结构的力学 参数；
步骤二，根据松散层拱理论，判定松散层拱是否可以形成以及形 成后的稳定性；
步骤三，根据关键层理论，通过对上覆岩层发挥的控制作用判别 关键层的位置；
步骤四，当松散层中无法形成松散层拱结构时，根据关键层理论， 判断关键层为远场关键层还是近场关键层；当关键层破断形成砌体梁 时，进入步骤五；
步骤五，当关键层破断形成砌体梁时，且Δ＞ΔT时，则砌体梁会 发生滑落失稳；其中，Δ为工作面回采后关键层与垮落的直接顶的间 距，ΔT为砌体梁保持自身稳定的临界回转量；
步骤六，当上覆岩层中关键层和松散层拱同时存在，且上覆岩层 中关键层和松散层拱保持稳定时，则地表会产生均匀下沉盆地；
当上覆岩层中只有近场关键层而无松散层拱，且上覆岩层中近场 关键层破断后，则地表会产生倾向台阶裂缝；
当上覆岩层中只存在远场关键层而无松散层拱，且上覆岩层中远 场关键层破断后，则地表会产生走向斑裂；
当关键层和松散层拱同时破断消失时，若砌体梁发生滑落失稳， 则地表会产生塌陷坑。
进一步，所述步骤一中，所述力学参数为松散层的容重、厚度、 矢高、内聚力、内摩擦角、侧压系数、基岩破断角、煤层顶界面与关 键层底界面间的距离；所述各个岩层的力学参数具体为：各个岩层的 容重、关键层的厚度、破断距、抗拉强度、所受载荷和直接顶的碎胀 系数和垮落高度，煤层的厚度、埋深和基岩中各个岩层厚度及岩性。 直接顶的垮落高度一般为1.67～2.17倍的采高；
进一步，所述步骤一中，具体步骤如下：
在煤炭开采区域地表施工地面钻孔取芯保存，并绘制钻孔综合柱 状图；根据钻孔综合柱状图，确定煤炭开采区域开采煤层的厚度、埋 深、松散层厚度和基岩中各岩层厚度及岩性；
通过对松散层进行取芯，获得的试样进行物理力学参数测试，从 而得到松散层容重、内摩擦角和内聚力；
通过对基岩中各个岩层进行取芯，获得的试样进行物理力学参数 测试，从而得到基岩中各个岩层的容重、弹性模量和各个岩层的抗拉 强度；直接顶的碎胀系数由基岩中岩石力学性质和岩层的破坏形式确 定。
在煤层开采区域，通过井下水压致裂地应力测量装置，确定松散 层的侧压系数；
利用关键层判别软件KSPB，基于松散层和基岩力学的特性，在 钻孔综合柱状图进行关键层的判别，从而得到关键层的厚度、位置、 破断距、所受载荷、关键层的抗拉强度、煤层顶界面与关键层底界面 间的距离；
利用钻孔测井分析仪，在工作面巷道顶板上方打钻孔，通过窥视 结果中钻孔的深度，绘制出窥视钻孔破断演化图，得出基岩破断角。
进一步，所述步骤二中，当开采区域地质条件确定后，松散层拱 的形成主要由松散层厚度决定，故松散层拱形成时计算得到的松散层 厚度即为松散层临界厚度HC，松散层拱稳定时计算得到的工作面采 宽为临界采宽LM。
进一步，所述步骤二中，具体过程如下：
图3为松散层拱二维力学模型图，该松散层拱二维力学模型的 建立方式方式为：以松散层拱左侧拱基为原点，拱基所在直线为X 轴，X轴向右为正方向，垂直于拱基所在直线的直线为Y轴，Y轴 向上为正方向建立松散层拱二维力学模型。
计算使松散层拱形成及稳定的松散层临界厚度：
Figure BDA0003483244450000111
式中：Larch、Harch、δarch为松散层拱的跨度、矢高和厚度；γ、α 为基岩容重和破断角；C、
Figure BDA0003483244450000112
为松散层的内聚力和内摩擦角；λ为侧压 系数；Σh为主关键层底界面与煤层顶界面间的距离；h0为松散层拱 上覆松散层厚度；HC为松散层临界厚度。
进一步，所述步骤二中，根据松散层拱的形成及稳定条件：当钻 孔收集的工作面上覆松散层厚度参数ΣH大于松散层拱形成时的临界 厚度HC时，松散层中能形成松散层拱；当松散层拱稳定时，据此推 出煤矿开采参数中的工作面临界宽度LM。
进一步，所述步骤三中，当确定了工作面钻孔柱状及柱状中不同 岩层的强度参数后，根据变形特征中的刚度条件和强度条件式，借助 关键层判别位置软件KSPB对工作面的关键层位置进行判别。
进一步，所述步骤三中，具体步骤如下：
由下式判别关键层的强度准则和刚度准则：
Figure BDA0003483244450000113
式中：q1|n+1、q1|n为第n+1和n层岩层对第1层岩层的载荷；ln+1、 ln为第n+1和n层岩层的破断距。
借助关键层判别位置软件KSPB对工作面的关键层进行判别，推 断出关键层在柱状图中的位置，从而判定关键层在工作面上覆岩层中 的位置。
进一步，所述步骤四中，当开采区域地质条件确定存在关键层， 无法形成松散层拱时，根据远近场关键层划分依据，将工作面上覆关 键层划分为远场关键层和近场关键层。
进一步，所述步骤四中，具体过程如下：
由下式确定a和b：
Figure BDA0003483244450000121
式中，a为工作面长度；b为关键层倾向悬露长度；lm关键层极 限跨距。
根据关键层破断规律：3lm>b>1.414lm，来压步距lm<a<1.414lm， 此时关键层“竖O-X”破断，关键层为近场关键层；lm<b<1.414lm，来 压步距a>1.414lm>b，此时关键层“横O-X”破断，关键层为远场关键 层。据此推出工作面上覆关键层为近场关键层或远场关键层。
进一步，所述步骤五中，当开采区域地质条件确定后，工作面采 高决定了砌体梁的稳定，故关键层中砌体梁稳定时计算得到的工作面 采高即为工作面临界采高MB，从而判断砌体梁的稳定性。
进一步，所述步骤五中，具体步骤如下：
由下式确定Δ和ΔT：
Figure BDA0003483244450000131
式中：MB为工作面临界采高；KP为直接顶的碎胀系数；Σhi为直 接顶的垮落高度；σC为关键层的抗压强度；hKS、qKS、lKS为关键层的 厚度、载荷和破断距；
根据关键层中砌体梁稳定性条件Δ＜ΔT，推出当工作面的实际采 高M小于MB时，砌体梁稳定。
将工作面临界采高MB与工作面实际采高M相比，即可判定工作 面上覆砌体梁的稳定性。
下面结合具体实施例对本发明的技术方案作进一步的描述。
山东某矿的松散层厚度为165.4～282.3m，平均为223.9m，基岩 厚度为35.45～85.15m，平均为60.3m。现以该矿3#煤层的6304首采 工作面钻孔提供的地质资料为基础，判定6304工作面在开采过程中 地表沉陷形式，具体步骤如下：
a.在3#煤层的首采工作面所在区域地表施工地面上钻孔取芯保 存，并绘制钻孔综合柱状图，通过钻孔综合柱状图(具体见图6)， 确定首采工作面所开采的3#煤层埋深319.14m，厚度平均3.40m，松 散层厚度223.9m，基岩中个岩层厚度以及岩性；
b.对松散层取芯试样物理力学参数测试，获得松散层容重为 20.25kN/m3，内摩擦角为7.5°，内聚力为1MPa；
c.对基岩中各个岩层进行取芯，对其试样进行物理力学参数测试， 得到基岩中各个岩层的抗压强度和容重，基于钻孔柱状图可知，基岩 中的岩性主要为泥岩、细砂岩和粉砂岩，强度相对较低，据此直接顶 的碎胀系数取1.25；
d.根据基岩中各个岩层容重和松散层容重计算覆岩的平均容重 为20.45kN/m3；
e.在煤矿开采区通过煤矿井下专用的小孔径水压致裂地应力测 量装置进行地应力测试，确定侧压系数为0.7；
f.根据关键层判别软件KSPB(国家计算机软件著作权编号为 2008SR34419)，基于松散层和基岩力学特性，根据钻孔综合柱状图 进行关键层判别，确定主关键层的埋深为282.30m，厚度为10m，抗 压强度为37MPa，所受荷载为7MPa，破断距为24.96m，关键层底界面与煤层顶界面间的距离为33.82m；
g.钻孔窥视采用钻孔测井分析仪，在轨道顺槽顶板上方打钻孔， 通过窥视结果钻孔的深度，绘制窥视钻孔破断位置演化图，获得基岩 破断角为75°。
将松散层拱参数(具体见表1)代入公式(1)求得工作面临界 宽度为97m，将关键层中砌体梁参数(具体见表2)代入公式(3) 得到工作面临界采高为3.49m。由于6304工作面上方覆岩既有松散 层又有关键层，因此不需要对关键层远近场进行判别。
根据本发明的发明内容，6304工作面上覆岩层既有松散层拱又 有关键层且松散层拱和关键层均能保持稳定，因此预计的沉陷类型为 地表均匀下沉盆地。为了现场监测地表沉陷规律，在6304工作面上 方布置两条观测线A和观测线B。经现场观测，观测线A累计最大 下沉2579mm，观测线B累计最大下沉2597mm(图4)，地表沉陷 表现为均匀下沉盆地(图5)。综上，当受采动影响时，地表移动只 会产生均匀下沉盆地，发明内容比较符合工程实践结果。
表1承载结构参数
Figure BDA0003483244450000151
表2关键层中砌体梁参数
Figure BDA0003483244450000161
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而 非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领 域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技 术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修 改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方 案的精神和范围。

Claims (7)

1.一种地表沉陷类别预计方法，其特征在于，包括以下步骤：
步骤一，收集钻孔周围松散层的力学参数和各个岩层结构的力学参数；
步骤二，根据松散层拱理论，判定松散层拱是否可以形成以及形成后的稳定性；
步骤三，根据关键层理论，通过对上覆岩层发挥的控制作用判别关键层的位置；
步骤四，当松散层中无法形成松散层拱结构时，根据关键层理论，判断关键层为远场关键层还是近场关键层；当关键层破断形成砌体梁时，进入步骤五；
步骤五，当关键层破断形成砌体梁且Δ＞ΔT时，则砌体梁会发生滑落失稳；其中，Δ为工作面回采后关键层与垮落的直接顶的间距，ΔT为砌体梁保持自身稳定的临界回转量；
步骤六，当上覆岩层中关键层和松散层拱同时存在，且上覆岩层中关键层和松散层拱保持稳定时，则地表会产生均匀下沉盆地；
当上覆岩层中只有近场关键层而无松散层拱，且上覆岩层中近场关键层破断后，则地表会产生倾向台阶裂缝；
当上覆岩层中只存在远场关键层而无松散层拱，且上覆岩层中远场关键层破断后，则地表会产生走向斑裂；
当关键层和松散层拱同时破断消失时，若砌体梁发生滑落失稳，则地表会产生塌陷坑。
2.根据权利要求1所述的地表沉陷类别预计方法，其特征在于，步骤一中，所述松散层的力学参数包括：松散层的容重、厚度、矢高、内聚力、内摩擦角、侧压系数、基岩破断角、煤层顶界面与关键层底界面间的距离；
所述各个岩层的力学参数包括：各个岩层的容重、关键层的厚度、破断距、关键层的抗拉强度、所受载荷、直接顶的碎胀系数、直接顶的垮落高度，煤层的厚度和埋深、基岩中各个岩层厚度及岩性。
3.根据权利要求2所述的地表沉陷类别预计方法，其特征在于，在煤炭开采区域地表施工地面钻孔取芯保存，并绘制钻孔综合柱状图；根据钻孔综合柱状图，确定煤炭开采区域开采煤层的厚度和埋深、松散层厚度和基岩中各岩层厚度及岩性；
通过对松散层进行取芯，获得的试样进行物理力学参数测试，从而得到松散层容重、内摩擦角和内聚力；
通过对基岩中各个岩层进行取芯，获得的试样进行物理力学参数测试，从而得到基岩中各个岩层的容重、弹性模量和各个岩层的抗拉强度；直接顶的碎胀系数由基岩中岩石力学性质和岩层的破坏形式确定；
在煤层开采区域，通过井下水压致裂地应力测量装置，确定松散层的侧压系数；
利用关键层判别软件KSPB，基于松散层和基岩力学的特性，在钻孔综合柱状图进行关键层的判别，从而得到关键层的厚度、位置、破断距、所受载荷、关键层的抗拉强度、关键层底界面与煤层顶界面间的距离；
利用钻孔测井分析仪，在工作面巷道顶板上方打钻孔，通过窥视结果中钻孔的深度，绘制出窥视钻孔破断演化图，得出基岩破断角。
4.根据权利要求1所述的地表沉陷类别预计方法，其特征在于，所述步骤二中，以松散层拱左侧拱基为原点，拱基所在直线为X轴，X轴向右为正方向，垂直于拱基所在直线的直线为Y轴，Y轴向上为正方向建立松散层拱二维力学模型；
根据该模型计算使松散层拱形成及稳定的松散层临界厚度：
Figure FDA0003483244440000031
式中：Larch、Harch、δarch为松散层拱的跨度、矢高和厚度；γ、α为基岩容重和破断角；C、
Figure FDA0003483244440000032
为松散层的内聚力和内摩擦角；λ为侧压系数；Σh为主关键层底界面与煤层顶界面间的距离；h0为松散层拱上覆松散层厚度；HC为松散层临界厚度；
根据松散层拱的形成及稳定条件：当Σh＞HC时，松散层中能形成松散层拱；当松散层拱稳定时，据此推出煤矿开采参数中的工作面临界宽度LM。
5.根据权利要求1所述的地表沉陷类别预计方法，其特征在于，所述步骤三中，由下式判别关键层的强度准则和刚度准则：
Figure FDA0003483244440000041
式中：q1|n+1、q1|n为第n+1和n层岩层对第1层岩层的载荷；ln+1、ln为第n+1和n层岩层的破断距；
借助关键层判别位置软件KSPB对工作面的关键层进行判别，推断出关键层在柱状图中的位置，从而判定关键层在工作面上覆岩层中的位置。
6.根据权利要求1所述的地表沉陷类别预计方法，其特征在于，所述步骤四中，当开采区域地质条件确定存在关键层，无法形成松散层拱时，根据远近场关键层划分依据，将工作面上覆关键层划分为远场关键层和近场关键层；
由下式确定a和b：
Figure FDA0003483244440000051
式中，a为工作面长度；b为关键层倾向悬露长度；lm关键层极限跨距；
根据关键层破断规律：3lm>b>1.414lm，来压步距lm<a<1.414lm，此时关键层“竖O-X”破断，关键层为近场关键层；lm<b<1.414lm，来压步距a>1.414lm>b，此时关键层“横O-X”破断，关键层为远场关键层，据此推出工作面上覆关键层为近场关键层或远场关键层。
7.根据权利要求1所述的地表沉陷类别预计方法，其特征在于，所述步骤五中，由下式确定Δ和ΔT：
Figure FDA0003483244440000052
式中：MB为工作面临界采高；KP为直接顶的碎胀系数；Σhi为直接顶的垮落高度；σC为关键层的抗压强度；hKS、qKS、lKS为关键层的厚度、载荷和破断距。
CN202210073917.9A 2022-01-21 2022-01-21 一种地表沉陷类别预计方法 Withdrawn CN114519257A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210073917.9A CN114519257A (zh) 2022-01-21 2022-01-21 一种地表沉陷类别预计方法
LU503279A LU503279B1 (en) 2022-01-21 2022-12-29 Surface subsidence category prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210073917.9A CN114519257A (zh) 2022-01-21 2022-01-21 一种地表沉陷类别预计方法

Publications (1)

Publication Number Publication Date
CN114519257A true CN114519257A (zh) 2022-05-20

Family

ID=81596987

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210073917.9A Withdrawn CN114519257A (zh) 2022-01-21 2022-01-21 一种地表沉陷类别预计方法

Country Status (2)

Country Link
CN (1) CN114519257A (zh)
LU (1) LU503279B1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115495879A (zh) * 2022-08-18 2022-12-20 山东科技大学 一种工作面近正断层开采保护煤柱宽度设计方法
CN117436291A (zh) * 2023-12-21 2024-01-23 山东科技大学 一种基于覆岩运动学理论的采场三带辨别方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115495879A (zh) * 2022-08-18 2022-12-20 山东科技大学 一种工作面近正断层开采保护煤柱宽度设计方法
CN115495879B (zh) * 2022-08-18 2024-04-30 山东科技大学 一种工作面近正断层开采保护煤柱宽度设计方法
CN117436291A (zh) * 2023-12-21 2024-01-23 山东科技大学 一种基于覆岩运动学理论的采场三带辨别方法

Also Published As

Publication number Publication date
LU503279B1 (en) 2023-07-04

Similar Documents

Publication Publication Date Title
Yang et al. Study on presplitting blasting the roof strata of adjacent roadway to control roadway deformation
CN112446143B (zh) 一种厚松散层薄基岩下条带充填开采参数设计方法
CN109798149B (zh) 一种坚硬顶板厚煤层采场来压强度分级预测方法
CN106919770A (zh) 一种基于数值模拟的损伤变量确定方法
CN114519257A (zh) 一种地表沉陷类别预计方法
Winn et al. Multi-approach stability analyses of large caverns excavated in low-angled bedded sedimentary rock masses in Singapore
CN112901169B (zh) 一种不对称孤岛工作面的冲击失稳判别与防冲开采方法
Liang et al. Study on Roof‐Coal Caving Characteristics with Complicated Structure by Fully Mechanized Caving Mining
Zhen et al. Comparative study on two types of nonpillar mining techniques by roof cutting and by filling artificial materials
Wu et al. Failure mechanism and stability control of surrounding rock in mining roadway with gentle slope and close distance
Liu et al. Stability analysis and monitoring method for the key block structure of the basic roof of noncoal pillar mining with automatically formed gob‐side entry
Xie et al. Formation Mechanism and the Height of the Water‐Conducting Fractured Zone Induced by Middle Deep Coal Seam Mining in a Sandy Region: A Case Study from the Xiaobaodang Coal Mine
Zhu et al. Study on the movement characteristics of the overlying stratum and surrounding rock control in ultraclose coal seams: a case study
Shi et al. Research on key technologies of floor heave control in soft rock roadway
Ozturk Support design of underground openings in an asphaltite mine
Hu et al. Presplitting blasting the roof strata to control large deformation in the deep mine roadway
Yang et al. Movement laws of overlying strata above a fully mechanized coal mining face backfilled with gangue: a case study in jiulishan coal mine in henan Province, China
Guo et al. Research on the settlement regulation and stability of large building foundation over gobs: a case study in the xiangcheng coal mine, China
Yu et al. Innovative Control Technique for the Floor Heave in Goaf‐Side Entry Retaining Based on Pressure Relief by Roof Cutting
Zhang et al. Tunnel collapse mechanism and its control strategy in fault fracture zone
Ju et al. Study on Fracture and Stress Evolution Characteristics of Ultra‐Thick Hard Sandstone Roof in the Fully Mechanized Mining Face with Large Mining Height: A Case Study of Xiaojihan Coal Mine in Western China
Xue et al. Numerical investigation on overburden migration behaviors in stope under thick magmatic rocks
Wang et al. Numerical investigation into evolution of crack and stress in residual coal pillars under the influence of longwall mining of the adjacent underlying coal seam
Hu et al. Study on evolution characteristics of regenerated roof structure in downward mining of bifurcated coal seam
Liu et al. Study on stability of underlying room and pillar old goaf in close coal seam and mining of the upper coal seam

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
WW01 Invention patent application withdrawn after publication

Application publication date: 20220520