CN113153418A - 一种基于综合指数法的三维空间充填效果评价方法 - Google Patents

一种基于综合指数法的三维空间充填效果评价方法 Download PDF

Info

Publication number
CN113153418A
CN113153418A CN202110374700.7A CN202110374700A CN113153418A CN 113153418 A CN113153418 A CN 113153418A CN 202110374700 A CN202110374700 A CN 202110374700A CN 113153418 A CN113153418 A CN 113153418A
Authority
CN
China
Prior art keywords
filling
dimensional space
degree
weight
membership
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
Application number
CN202110374700.7A
Other languages
English (en)
Other versions
CN113153418B (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.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
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 CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202110374700.7A priority Critical patent/CN113153418B/zh
Publication of CN113153418A publication Critical patent/CN113153418A/zh
Application granted granted Critical
Publication of CN113153418B publication Critical patent/CN113153418B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH 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
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH 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 DRILLING; MINING
    • E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
    • E21F17/18—Special adaptations of signalling or alarm devices

Abstract

本发明公开了一种基于综合指数法的三维空间充填效果评价方法，具体涉及矿业工程领域。本发明通过选取充填体充填质量、覆岩运动程度和地表变形特征作为三维空间充填效果评价参数，根据三维空间充填稳定后的现场实际测量结果、钻孔电视资料和地表监测结果，获取充填体充填质量、覆岩运动程度和地表变形特征，分析确定三维空间充填效果评价参数的影响因素，基于综合指数熵权法计算各影响因素的权重，构建三维空间充填开采评价函数，计算三维空间的充填开采评价结果，并根据充填开采评价结果评价三维空间的充填效果。本发明综合各影响因素对三维空间充填效果的影响，实现了对三维空间充填效果的准确评价，为充填开采的安全生产提供了依据。

Description

一种基于综合指数法的三维空间充填效果评价方法
技术领域
本发明涉及矿业工程领域，具体涉及一种基于综合指数法的三维空间充填效果评价方法。
背景技术
随着我国东部地区煤炭资源的大量开采，剩余煤炭资源多为“三下压煤”，大量优质煤炭资源急需充填开采，为了在提高采出率的同时延长矿井服务年限，充分贯彻绿色发展的战略目标，因此亟需从充填体、工作面围岩及地面三个方面对采煤作业面地下三维空间的充填效果进行综合评价，综合考虑各种影响因素准确分析三维空间的充填效果，为充填开采的安全生产提供更加科学有效的判断依据。
发明内容
本发明旨在解决上述问题，提出了一种基于综合指数法的三维空间充填效果评价方法。该方法从充填体、工作面围岩以及地面三个方面对三维空间的充填效果进行综合评价，充分考虑各种因素对三维空间充填效果的影响，有利于三维空间充填效果的准确评价，为充填开采的安全生产提供了有效的判断依据。
本发明采用以下的技术方案：
一种基于综合指数法的三维空间充填效果评价方法，具体包括以下步骤：
步骤1，选取充填体充填质量、覆岩运动程度和地表变形特征作为三维空间充填效果评价参数，根据采煤工作面三维空间充填稳定后的现场实际测量结果、钻孔电视资料和地表监测结果，确定充填体充填质量、覆岩运动程度和地表变形特征；
步骤2，根据获得的充填体充填质量、覆岩运动程度和地表变形特征，分别确定影响充填体充填质量、覆岩运动程度和地表变形特征的影响因素，得到三维空间充填效果评价参数的影响因素；
所述充填体充填质量受等价采厚、充填结石体的物理力学特征和充填程度的影响，其中，等价采厚的影响因素包括顶板超前下沉量、欠接顶量和充填体压实量，充填结石体的物理力学特征的影响因素包括抗压强度、弹性模量和孔隙度，充填程度的影响因素包括充填体应力、顶底板移近量和巷道两帮移近量；
所述覆岩运动程度受顶板离层裂隙、工作面侧向支承压力和工作面超前支承压力的影响，其中，顶板离层裂隙的影响因素包括空间裂隙相对面积比、顶板岩层裂隙漏失量和顶板岩层相对位移量，工作面侧向支承压力的影响因素包括孔深5m处的侧向支承压力、孔深10m处的侧向支承压力和孔深15m处的侧向支承压力，工作面超前支承压力的影响因素包括距离工作面20m处的工作面超前支承压力、距离工作面40m处的工作面超前支承压力和距离工作面60m处的工作面超前支承压力；
所述地表变形特征的影响因素包括地表下沉、地表倾斜变形、地表水平移动、曲率和地表水平变形；
步骤3，基于综合指数熵权法，结合三维空间充填稳定后的现场实际测量结果，分别计算三维空间充填效果评价参数中各影响因素的权重；
步骤4，根据三维空间充填效果评价参数中各影响因素的权重，构建三维空间充填开采评价函数；
步骤5，根据现场实际测量结果、钻孔电视资料和地表监测结果，利用三维空间充填开采评价函数计算三维空间的充填开采评价结果θ，并根据充填开采评价结果对三维空间的充填效果进行评价。
优选地，所述步骤3中，三维空间充填效果评价参数中各影响因素的权重计算公式，如式(1)所示：
Figure BDA0003010718500000021
其中，
Figure BDA0003010718500000022
Figure BDA0003010718500000023
式中，yij为三维空间充填效果评价参数中第i个参数中第j个影响因素所占的比重；m为影响因素总数；ej为第i个参数中第j个影响因素的信息熵；K为常数，
Figure BDA0003010718500000024
wj表示第i个参数中第j个影响因素的权重。
优选地，所述步骤4中，具体包括以下步骤：
步骤4.1，构建充填体充填质量评价函数α，如式(5)所示：
α＝α等价+α物理+α充填 (5)
其中，等价采厚的隶属度α等价如式(6)所示：
Figure BDA0003010718500000031
式中，α超前为顶板超前下沉量的隶属度，α欠接为欠接顶量的隶属度，α压实为充填体压实量的隶属度；a1为顶板超前下沉量的权重，a2为欠接顶量的权重，a3为充填体压实量的权重；H为开采厚度，单位为m；Δh超前为顶板超前下沉量，单位为m；Δh欠接为欠接顶量，单位为m；Δh压实为充填体压实量，单位为m；
充填结石体物理力学特征的隶属度α物理如式(7)所示：
Figure BDA0003010718500000032
式中，αR为抗压强度的隶属度，αE为弹性模量的隶属度，αn为孔隙率的隶属度；b1为抗压强度的权重，b2为弹性模量的权重，b3为孔隙率的权重；R结石为充填结石体的抗压强度，单位为MPa；R煤为煤的抗压强度，单位为MPa；E结石为充填结石体的弹性模量，单位为MPa；n结石为充填结石体的孔隙率；n煤为煤的孔隙率；
充填程度的隶属度α充填如式(8)所示：
Figure BDA0003010718500000033
式中，
Figure BDA0003010718500000034
为充填体应力的隶属度，
Figure BDA0003010718500000035
为顶底板移近量的隶属度，
Figure BDA0003010718500000036
为巷道两帮移近量的隶属度；c1为充填体应力的权重，c2为顶底板移近量的权重，c3为巷道两帮移近量的权重；
Figure BDA0003010718500000037
为监测到的最大应力，单位为MPa；
Figure BDA0003010718500000038
为工作面处的垂直地应力，单位为MPa；H为开采厚度，单位为m；ΔH顶底为顶底板移近量，单位为m；W为巷道宽度，单位为m；ΔW两帮为巷道两帮移近量，单位为m；
步骤4.2，构建覆岩运动程度评价函数β，如式(9)所示：
β＝β离层+β侧向+β超前 (9)
其中，顶板离层裂隙β离层的隶属度如式(10)所示：
Figure BDA0003010718500000041
式中，βS为空间裂隙相对面积比的隶属度，βV为顶板岩层裂隙漏失量的隶属度，βL为顶板岩层相对位移量的隶属度；d1为空间裂隙相对面积比的权重，d2为顶板岩层裂隙漏失量的权重，d3为顶板岩层相对位移量的权重；S钻孔为两带内的钻孔表面积，单位为m2；S裂隙为两带内的裂隙表面积，单位为m2；Vmax为两带内的阶段最大漏失量，单位为m3；V平均为两带内的阶段平均漏失量，单位为m3；L为两带高度，单位为m；L相对为两带内相对位移量之和，单位为m；
工作面侧向支承压力β侧向的隶属度如式(11)所示：
Figure BDA0003010718500000042
式中，β5为孔深5m处侧向支承压力的隶属度，β10为孔深10m处侧向支承压力的隶属度，β15为孔深15m处侧向支承压力的隶属度；e1为孔深5m处侧向支承压力的权重，e2为孔深10m处侧向支承压力的权重，e3为孔深15m处侧向支承压力的权重；Pmax为最大侧向支承压力，单位为MPa；P5为孔深5m处的侧向支承压力，单位为MPa；P10为孔深10m处的侧向支承压力，单位为MPa；P15为孔深15m处的侧向支承压力，单位为MPa；
工作面超前支承压力β超前的隶属度如式(12)所示：
Figure BDA0003010718500000043
式中，βC20为距离采煤工作面20m处超前支承压力的隶属度，βC40为距离采煤工作面40m处超前支承压力的隶属度，βC60为距离采煤工作面60m处超前支承压力的隶属度；f1为距离采煤工作面20m处超前支承压力的权重，f2为距离采煤工作面40m处超前支承压力的权重，f3为距离采煤工作面60m处超前支承压力的权重；Nmax为最大超前支承压力，单位为MPa；N20为距离采煤工作面20m处的超前支承压力，单位为MPa；N40为距离采煤工作面40m处的超前支承压力，单位为MPa；N60为距离采煤工作面60m处的超前支承压力，单位为MPa；
步骤4.3，构建地表变形特征评价函数γ，如式(13)所示：
Figure BDA0003010718500000051
式中，γ下沉为地表下沉的隶属度，γ倾斜为地表倾斜的隶属度，γ水平为地表水平移动的隶属度，γ曲率为曲率的隶属度，γ变形为地表水平变形的隶属度；g为地表下沉的权重，h为地表倾斜的权重，I为地表水平移动的权重，J为曲率的权重，k为地表水平变形的权重；Hmax为地表最大下沉量，单位为m；H为开采厚度，单位为m；εmax为地表最大水平变形量，单位为m；Fmax为地表最大倾斜值；
步骤4.4，基于充填体充填质量评价函数α、覆岩运动程度评价函数β和地表变形特征评价函数γ，构建三维空间充填开采评价函数，如式(14)所示：
θ＝α+β+γ (14)
式中，θ为三维空间的充填开采评价结果。
优选地，所述步骤4.2中，两带为三维空间开采后因覆岩破坏和移动而形成的冒落带和裂隙带。
优选地，所述步骤5中，当0.9≤θ＜1时，三维空间的充填效果评价为优秀；当0.8≤θ＜0.9时，三维空间的充填效果为良好；当0.7≤θ＜0.8时，三维空间的充填效果为一般；当0.6≤θ＜0.7时，三维空间的充填效果为差；当θ＜0.6时，三维空间充填不合格。
优选地，所述步骤5中，当三维空间的充填效果评价为优秀或良好时，可以对三维空间进行正常开采；当三维空间的充填效果评价为一般时，需要对三维空间中的部分充填区域改进，提高充填质量；当三维空间的充填效果评价为差时，需要对三维空间中的部分生产区域停工改进，在保证工作面能够安全生产的前提下，三维空间才能复工复产；当三维空间的充填效果评价为不合格时，需要对三维空间中的所有工作区域停工停产，并对三维空间内的所有工作区域进行全面检查，对三维空间进行二次充填待生产区域围岩稳定后，经相关部门批准才能复工复产。
本发明具有如下有益效果：
1、本发明方法基于现场实际测量结果、钻孔电视资料和地表监测结果，综合考虑影响三维空间充填效果的影响因素，建立三维空间充填开采评价函数，通过计算各影响因素的权重及其隶属度，叠加后得到三维空间的充填开采评价结果，该充填开采评价结果更加符合三维空间充填后的实际情况，参照性强、可信度高，具有较强的实际应用价值。
2、本发明方法从充填体、工作面围岩以及地面三个方面对三维空间的充填效果进行综合评价，充分考虑了多种因素对三维空间充填效果的影响，实现了对三维空间充填效果的准确评价，为充填开采的安全生产提供了有效的判断依据，有利于复杂地质条件下矿产的高效开发和安全开采。
附图说明
图1为三维空间充填效果评价参数影响因素的示意图。
具体实施方式
下面结合附图和某碳酸岩区块例井为例，对本发明的具体实施方式做进一步说明：
以采煤作业后的充填采煤工作面为例，采用本发明提出的一种基于综合指数法的三维空间充填效果评价方法，如图1所示，具体包括以下步骤：
步骤1，选取充填体充填质量、覆岩运动程度和地表变形特征作为三维空间充填效果评价参数，根据采煤工作面三维空间充填稳定后的现场实际测量结果、钻孔电视资料和地表监测结果，确定充填体充填质量、覆岩运动程度和地表变形特征。
步骤2，根据获得的充填体充填质量、覆岩运动程度和地表变形特征，分别确定影响充填体充填质量、覆岩运动程度和地表变形特征的影响因素。
结合现场实际测量结果、钻孔电视资料和地表监测结果，分别对充填体充填质量、覆岩运动程度和地表变形特征的影响因素进行分析，如图1所示，得到充填体充填质量受等价采厚、充填结石体的物理力学特征和充填程度的影响，其中，等价采厚受顶板超前下沉量、欠接顶量和充填体压实量的影响，充填结石体的物理力学特征受抗压强度、弹性模量和孔隙度的影响，充填程度受充填体应力、顶底板移近量和巷道两帮移近量的影响；覆岩运动程度受顶板离层裂隙、工作面侧向支承压力和工作面超前支承压力的影响，其中，顶板离层裂隙受空间裂隙相对面积比、顶板岩层裂隙漏失量和顶板岩层相对位移量的影响，工作面侧向支承压力受孔深5m处的侧向支承压力、孔深10m处的侧向支承压力和孔深15m处的侧向支承压力的影响，工作面超前支承压力受距离工作面20m处的工作面超前支承压力、距离工作面40m处的工作面超前支承压力和距离工作面60m处的工作面超前支承压力的影响；地表变形特征的影响因素受地表下沉、地表倾斜变形、地表水平移动、曲率和地表水平变形的影响。
由此可得，三维空间充填效果评价参数的影响因素包括顶板超前下沉量、欠接顶量、充填体压实量、抗压强度、弹性模量、孔隙度、充填体应力、顶底板移近量、巷道两帮移近量、空间裂隙相对面积比、顶板岩层裂隙漏失量、顶板岩层相对位移量、孔深5m处的侧向支承压力、孔深10m处的侧向支承压力、孔深15m处的侧向支承压力、距离工作面20m处的工作面超前支承压力、距离工作面40m处的工作面超前支承压力、距离工作面60m处的工作面超前支承压力、地表下沉、地表倾斜变形、地表水平移动、曲率和地表水平变形。
步骤3，基于综合指数熵权法，结合三维空间充填稳定后的现场实际测量结果，利用公式(1)至公式(3)，分别计算三维空间充填效果评价参数中各影响因素的权重。
计算得到顶板超前下沉量的权重a1＝3％，欠接顶量的权重a2＝3％，充填体压实量的权重a3＝4％；抗压强度的权重b1＝3％，弹性模量的权重b2＝1％，孔隙度的权重b3＝1％；充填体应力的权重c1＝3％、顶底板移近量的权重c2＝4％、巷道两帮移近量的权重c3＝3％；空间裂隙相对面积比的权重d1＝7％、顶板岩层裂隙漏失量的权重d2＝7％、顶板岩层相对位移量的权重d3＝5％；孔深5m处侧向支承压力的权重e1＝1％、孔深10m处侧向支承压力的权重e2＝1％、孔深15m处侧向支承压力的权重e3＝1％；距离工作面20m处工作面超前支承压力的权重f1＝1％、距离工作面40m处工作面超前支承压力的权重f2＝1％、距离工作面60m处工作面超前支承压力的权重f3＝1％；地表下沉的权重g＝15％、地表倾斜变形的权重h＝10％、地表水平移动的权重I＝5％、曲率的权重J＝10％和地表水平变形的权重k＝10％。
由此得到，充填体充填质量的权重为25％，其中，等价采厚的权重为10％，充填结石体的物理力学特征的权重为5％，充填程度的权重为10％；覆岩运动程度的权重为25％，其中，顶板离层裂隙的权重为19％，工作面侧向支承压力的权重为3％，工作面超前支承压力的权重为3％；地表变形特征的权重为50％，其中，地表下沉的权重为15％、地表倾斜变形的权重为10％，地表水平移动的权重为5％，曲率的权重为10％，地表水平变形的权重为10％。
步骤4，根据三维空间充填效果评价参数中各影响因素的权重，构建三维空间充填开采评价函数，具体包括以下步骤：
步骤4.1，构建充填体充填质量评价函数α，如式(5)所示：
α＝α等价+α物理+α充填 (5)
其中，
Figure BDA0003010718500000081
α物理＝(0.03R结石/R煤+0.01E结石/E煤+0.01n结石/n煤)×100％ (16)
Figure BDA0003010718500000082
式中，α等价为等价采厚的隶属度，α物理为充填结石体物理力学特征的隶属度，α充填为充填程度的隶属度；H为开采厚度，单位为m；Δh超前为顶板超前下沉量，单位为m；Δh欠接为欠接顶量，单位为m；Δh压实为充填体压实量，单位为m；R结石为充填结石体的抗压强度，单位为MPa；R煤为煤的抗压强度，单位为MPa；E结石为充填结石体的弹性模量，单位为MPa；n结石为充填结石体的孔隙率；n煤为煤的孔隙率；
Figure BDA0003010718500000083
为监测到的最大应力，单位为MPa；
Figure BDA0003010718500000084
为工作面处的垂直地应力，单位为MPa；ΔH顶底为顶底板移近量，单位为m；W为巷道宽度，单位为m；ΔW两帮为巷道两帮移近量，单位为m。
步骤4.2，构建覆岩运动程度评价函数β，如式(9)所示：
β＝β离层+β侧向+β超前 (9)
其中，
Figure BDA0003010718500000085
Figure BDA0003010718500000086
Figure BDA0003010718500000087
式中，β离层为顶板离层裂隙的隶属度，βS为空间裂隙相对面积比的隶属度，β侧向为工作面侧向支承压力的隶属度；S钻孔为冒落带和裂隙带内的钻孔表面积，单位为m2；S裂隙为冒落带和裂隙带内的裂隙表面积，单位为m2；Vmax为冒落带和裂隙带内的阶段最大漏失量，单位为m3；V平均为冒落带和裂隙带内的阶段平均漏失量，单位为m3；L为冒落带和裂隙带的高度，单位为m；L相对为冒落带和裂隙带内相对位移量之和，单位为m；Pmax为最大侧向支承压力，单位为MPa；P5为孔深5m处的侧向支承压力，单位为MPa；P10为孔深10m处的侧向支承压力，单位为MPa；P15为孔深15m处的侧向支承压力，单位为MPa；Nmax为最大超前支承压力，单位为MPa；N20为距离采煤工作面20m处的超前支承压力，单位为MPa；N40为距离采煤工作面40m处的超前支承压力，单位为MPa；N60为距离采煤工作面60m处的超前支承压力，单位为MPa。
步骤4.3，构建地表变形特征评价函数γ，如式(13)所示：
Figure BDA0003010718500000091
式中，Hmax为地表最大下沉量，单位为m；H为开采厚度，单位为m；εmax为地表最大水平变形量，单位为m；Fmax为地表最大倾斜值。
步骤4.4，基于充填体充填质量评价函数α、覆岩运动程度评价函数β和地表变形特征评价函数γ，构建三维空间充填开采评价函数，如式(14)所示：
θ＝α+β+γ (14)
式中，θ表示充填开采评价结果。
步骤5，根据三维空间充填开采评价函数，计算得到三维空间的充填开采评价结果θ，并根据充填开采评价结果对三维空间的充填效果进行评价。
当0.9≤θ＜1时，三维空间的充填效果评价为优秀；当0.8≤θ＜0.9时，三维空间的充填效果为良好；当0.7≤θ＜0.8时，三维空间的充填效果为一般；当0.6≤θ＜0.7时，三维空间的充填效果为差；当θ＜0.6时，三维空间充填不合格。
当三维空间的充填效果评价为优秀或良好时，三维空间可以正常开采；当三维空间的充填效果评价为一般时，需要对三维空间进行改进，提高部分充填区域的充填质量；当三维空间的充填效果评价为差时，需要对三维空间中的部分生产区域停工改进，在保证工作面能够安全生产的前提下，才能对三维空间复工复产；当三维空间的充填效果评价为不合格时，三维空间中的所有工作区域需要停工停产，通过对三维空间内的所有工作区域进行全面检查，对三维空间进行二次充填，待生产区域围岩稳定后，经相关部门批准后三维空间才能复工复产。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。

Claims (6)

1.一种基于综合指数法的三维空间充填效果评价方法，其特征在于，具体包括以下步骤：
步骤1，选取充填体充填质量、覆岩运动程度和地表变形特征作为三维空间充填效果评价参数，根据采煤工作面三维空间充填稳定后的现场实际测量结果、钻孔电视资料和地表监测结果，确定充填体充填质量、覆岩运动程度和地表变形特征；
步骤2，根据获得的充填体充填质量、覆岩运动程度和地表变形特征，分别确定影响充填体充填质量、覆岩运动程度和地表变形特征的影响因素，得到三维空间充填效果评价参数的影响因素；
所述充填体充填质量受等价采厚、充填结石体的物理力学特征和充填程度的影响，其中，等价采厚的影响因素包括顶板超前下沉量、欠接顶量和充填体压实量，充填结石体的物理力学特征的影响因素包括抗压强度、弹性模量和孔隙度，充填程度的影响因素包括充填体应力、顶底板移近量和巷道两帮移近量；
所述覆岩运动程度受顶板离层裂隙、工作面侧向支承压力和工作面超前支承压力的影响，其中，顶板离层裂隙的影响因素包括空间裂隙相对面积比、顶板岩层裂隙漏失量和顶板岩层相对位移量，工作面侧向支承压力的影响因素包括孔深5m处的侧向支承压力、孔深10m处的侧向支承压力和孔深15m处的侧向支承压力，工作面超前支承压力的影响因素包括距离工作面20m处的工作面超前支承压力、距离工作面40m处的工作面超前支承压力和距离工作面60m处的工作面超前支承压力；
所述地表变形特征的影响因素包括地表下沉、地表倾斜变形、地表水平移动、曲率和地表水平变形；
步骤3，基于综合指数熵权法，结合三维空间充填稳定后的现场实际测量结果，分别计算三维空间充填效果评价参数中各影响因素的权重；
步骤4，根据三维空间充填效果评价参数中各影响因素的权重，构建三维空间充填开采评价函数；
步骤5，根据现场实际测量结果、钻孔电视资料和地表监测结果，利用三维空间充填开采评价函数计算三维空间的充填开采评价结果θ，并根据充填开采评价结果对三维空间的充填效果进行评价。
2.根据权利要求1所述的一种基于综合指数法的三维空间充填效果评价方法，其特征在于，所述步骤3中，三维空间充填效果评价参数中各影响因素的权重计算公式，如式(1)所示：
Figure FDA0003010718490000011
其中，
Figure FDA0003010718490000021
Figure FDA0003010718490000022
式中，yij为三维空间充填效果评价参数中第i个参数中第j个影响因素所占的比重；m为影响因素总数；ej为第i个参数中第j个影响因素的信息熵；K为常数，
Figure FDA0003010718490000023
wj表示第i个参数中第j个影响因素的权重。
3.根据权利要求1所述的一种基于综合指数法的三维空间充填效果评价方法，其特征在于，所述步骤4中，具体包括以下步骤：
步骤4.1，构建充填体充填质量评价函数α，如式(5)所示：
α＝α等价+α物理+α充填 (5)
其中，等价采厚的隶属度α等价如式(6)所示：
Figure FDA0003010718490000024
式中，α超前为顶板超前下沉量的隶属度，α欠接为欠接顶量的隶属度，α压实为充填体压实量的隶属度；a1为顶板超前下沉量的权重，a2为欠接顶量的权重，a3为充填体压实量的权重；H为开采厚度，单位为m；Δh超前为顶板超前下沉量，单位为m；Δh欠接为欠接顶量，单位为m；Δh压实为充填体压实量，单位为m；
充填结石体物理力学特征的隶属度α物理如式(7)所示：
Figure FDA0003010718490000025
式中，αR为抗压强度的隶属度，αE为弹性模量的隶属度，αn为孔隙率的隶属度；b1为抗压强度的权重，b2为弹性模量的权重，b3为孔隙率的权重；R结石为充填结石体的抗压强度，单位为MPa；R煤为煤的抗压强度，单位为MPa；E结石为充填结石体的弹性模量，单位为MPa；n结石为充填结石体的孔隙率；n煤为煤的孔隙率；
充填程度的隶属度α充填如式(8)所示：
Figure FDA0003010718490000031
式中，
Figure FDA0003010718490000032
为充填体应力的隶属度，
Figure FDA0003010718490000033
为顶底板移近量的隶属度，
Figure FDA0003010718490000034
为巷道两帮移近量的隶属度；c1为充填体应力的权重，c2为顶底板移近量的权重，c3为巷道两帮移近量的权重；
Figure FDA0003010718490000035
为监测到的最大应力，单位为MPa；
Figure FDA0003010718490000036
为工作面处的垂直地应力，单位为MPa；H为开采厚度，单位为m；ΔH顶底为顶底板移近量，单位为m；W为巷道宽度，单位为m；ΔW两帮为巷道两帮移近量，单位为m；
步骤4.2，构建覆岩运动程度评价函数β，如式(9)所示：
β＝β离层+β侧向+β超前 (9)
其中，顶板离层裂隙β离层的隶属度如式(10)所示：
Figure FDA0003010718490000037
式中，βS为空间裂隙相对面积比的隶属度，βV为顶板岩层裂隙漏失量的隶属度，βL为顶板岩层相对位移量的隶属度；d1为空间裂隙相对面积比的权重，d2为顶板岩层裂隙漏失量的权重，d3为顶板岩层相对位移量的权重；S钻孔为两带内的钻孔表面积，单位为m2；S裂隙为两带内的裂隙表面积，单位为m2；Vmax为两带内的阶段最大漏失量，单位为m3；V平均为两带内的阶段平均漏失量，单位为m3；L为两带高度，单位为m；L相对为两带内相对位移量之和，单位为m；
工作面侧向支承压力β侧向的隶属度如式(11)所示：
Figure FDA0003010718490000041
式中，β5为孔深5m处侧向支承压力的隶属度，β10为孔深10m处侧向支承压力的隶属度，β15为孔深15m处侧向支承压力的隶属度；e1为孔深5m处侧向支承压力的权重，e2为孔深10m处侧向支承压力的权重，e3为孔深15m处侧向支承压力的权重；Pmax为最大侧向支承压力，单位为MPa；P5为孔深5m处的侧向支承压力，单位为MPa；P10为孔深10m处的侧向支承压力，单位为MPa；P15为孔深15m处的侧向支承压力，单位为MPa；
工作面超前支承压力β超前的隶属度如式(12)所示：
Figure FDA0003010718490000042
式中，βC20为距离采煤工作面20m处超前支承压力的隶属度，βC40为距离采煤工作面40m处超前支承压力的隶属度，βC60为距离采煤工作面60m处超前支承压力的隶属度；f1为距离采煤工作面20m处超前支承压力的权重，f2为距离采煤工作面40m处超前支承压力的权重，f3为距离采煤工作面60m处超前支承压力的权重；Nmax为最大超前支承压力，单位为MPa；N20为距离采煤工作面20m处的超前支承压力，单位为MPa；N40为距离采煤工作面40m处的超前支承压力，单位为MPa；N60为距离采煤工作面60m处的超前支承压力，单位为MPa；
步骤4.3，构建地表变形特征评价函数γ，如式(13)所示：
Figure FDA0003010718490000043
式中，γ下沉为地表下沉的隶属度，γ倾斜为地表倾斜的隶属度，γ水平为地表水平移动的隶属度，γ曲率为曲率的隶属度，γ变形为地表水平变形的隶属度；g为地表下沉的权重，h为地表倾斜的权重，I为地表水平移动的权重，J为曲率的权重，k为地表水平变形的权重；Hmax为地表最大下沉量，单位为m；H为开采厚度，单位为m；εmax为地表最大水平变形量，单位为m；Fmax为地表最大倾斜值；
步骤4.4，基于充填体充填质量评价函数α、覆岩运动程度评价函数β和地表变形特征评价函数γ，构建三维空间充填开采评价函数，如式(14)所示：
θ＝α+β+γ (14)
式中，θ为三维空间的充填开采评价结果。
4.根据权利要求1所述的一种基于综合指数法的三维空间充填效果评价方法，其特征在于，所述步骤4.2中，两带为三维空间开采后因覆岩破坏和移动而形成的冒落带和裂隙带。
5.根据权利要求1所述的一种基于综合指数法的三维空间充填效果评价方法，其特征在于，所述步骤5中，当0.9≤θ＜1时，三维空间的充填效果评价为优秀；当0.8≤θ＜0.9时，三维空间的充填效果为良好；当0.7≤θ＜0.8时，三维空间的充填效果为一般；当0.6≤θ＜0.7时，三维空间的充填效果为差；当θ＜0.6时，三维空间充填不合格。
6.根据权利要求1所述的一种基于综合指数法的三维空间充填效果评价方法，其特征在于，所述步骤5中，当三维空间的充填效果评价为优秀或良好时，可以对三维空间进行正常开采；当三维空间的充填效果评价为一般时，需要对三维空间中的部分充填区域改进，提高充填质量；当三维空间的充填效果评价为差时，需要对三维空间中的部分生产区域停工改进，在保证工作面能够安全生产的前提下，三维空间才能复工复产；当三维空间的充填效果评价为不合格时，需要对三维空间中的所有工作区域停工停产，并对三维空间内的所有工作区域进行全面检查，对三维空间进行二次充填待生产区域围岩稳定后，经相关部门批准才能复工复产。
CN202110374700.7A 2021-04-08 2021-04-08 一种基于综合指数法的三维空间充填效果评价方法 Active CN113153418B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110374700.7A CN113153418B (zh) 2021-04-08 2021-04-08 一种基于综合指数法的三维空间充填效果评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110374700.7A CN113153418B (zh) 2021-04-08 2021-04-08 一种基于综合指数法的三维空间充填效果评价方法

Publications (2)

Publication Number Publication Date
CN113153418A true CN113153418A (zh) 2021-07-23
CN113153418B CN113153418B (zh) 2022-07-19

Family

ID=76889264

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110374700.7A Active CN113153418B (zh) 2021-04-08 2021-04-08 一种基于综合指数法的三维空间充填效果评价方法

Country Status (1)

Country Link
CN (1) CN113153418B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116241326A (zh) * 2022-11-09 2023-06-09 华能煤炭技术研究有限公司 保护层充填开采关键参数设计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102011611A (zh) * 2010-11-30 2011-04-13 淄博市王庄煤矿 控制上覆岩层移动变形的高水膨胀材料条带式充填方法
CN103498696A (zh) * 2013-10-02 2014-01-08 中国矿业大学 一种从顶板巷道注浆充填压实采空区充填体的方法
CN103902780A (zh) * 2014-04-08 2014-07-02 中国矿业大学 固体充填采煤地表变形预计方法
CN105699432A (zh) * 2016-01-13 2016-06-22 太原理工大学 一种膏体充填效果的评价方法
CN111042821A (zh) * 2019-11-12 2020-04-21 山东科技大学 一种定量判断特厚煤层沿空巷道开掘时间的双指标方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102011611A (zh) * 2010-11-30 2011-04-13 淄博市王庄煤矿 控制上覆岩层移动变形的高水膨胀材料条带式充填方法
CN103498696A (zh) * 2013-10-02 2014-01-08 中国矿业大学 一种从顶板巷道注浆充填压实采空区充填体的方法
CN103902780A (zh) * 2014-04-08 2014-07-02 中国矿业大学 固体充填采煤地表变形预计方法
CN105699432A (zh) * 2016-01-13 2016-06-22 太原理工大学 一种膏体充填效果的评价方法
CN111042821A (zh) * 2019-11-12 2020-04-21 山东科技大学 一种定量判断特厚煤层沿空巷道开掘时间的双指标方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116241326A (zh) * 2022-11-09 2023-06-09 华能煤炭技术研究有限公司 保护层充填开采关键参数设计方法
CN116241326B (zh) * 2022-11-09 2024-04-26 华能煤炭技术研究有限公司 保护层充填开采关键参数设计方法

Also Published As

Publication number Publication date
CN113153418B (zh) 2022-07-19

Similar Documents

Publication Publication Date Title
CA3060277C (en) Mine exploitation method based on stopping, seperation and filling control
CN105626071B (zh) 一种缓倾斜薄矿体采矿方法
CN112446143B (zh) 一种厚松散层薄基岩下条带充填开采参数设计方法
CN104712358A (zh) 基于首采全岩卸压工作面沿空留巷的高瓦斯煤层群卸压共采方法
CN110295909A (zh) 只留齿形矿柱的充填体底下采矿方法
CN110030013B (zh) 一种过渡支架区三缝周期切顶自成巷帮的沿空留巷方法
CN108049870B (zh) 上盘含不稳岩层的急倾斜中厚矿体的诱导冒落采矿方法
CN113153418B (zh) 一种基于综合指数法的三维空间充填效果评价方法
CN106677781B (zh) 一种极薄煤层钻采卸压增透方法
CN114251103A (zh) 一种定向切缝压裂顶板大巷防冲护巷方法及安全采矿方法
CN108150171A (zh) 一种大埋深薄基岩区厚煤层综放开采压架突水防治方法
CN114329922B (zh) 基于结构性覆岩的导水裂隙带高度确定方法
CN111042821B (zh) 一种定量判断特厚煤层沿空巷道开掘时间的双指标方法
CN114809992A (zh) 一种低渗储层煤系气全生命周期高效抽采方法
CN110863828B (zh) 一种特厚煤层沿空掘巷紧张接替的定量解决方法
CN111008758B (zh) 基于双指标权重法的特厚煤层沿空煤巷开掘时间设计方法
CN114444266A (zh) 一种近松散层下煤层提高回采上限安全开采的可行性评价方法
CN113236248B (zh) 一种考虑悬臂结构破断过程的切顶角度计算方法
CN113982581B (zh) 基于低碳开采的覆岩渗流隔离带稳定性控制方法
CN115199273A (zh) 基于地表沉陷控制要求的采煤-充填工作面布局方法
Pu et al. Determination of reasonable width of filling body for gob-side entry retaining in mining face with large cutting height
Li et al. Isolated overburden grout injection technology mining and grouting parameters discussion and optimization
Chapula et al. The evaluation of large deformation in deep underground mine excavations using RQD
CN117823044A (zh) 一种千米水平钻压裂提高本工作面瓦斯抽采效果与邻近巷道围岩稳定性的方法
CN118014414A (zh) 一种考虑岩层强度的采深采厚比计算与评价方法

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