一种巷道过陷落柱围岩支护参数智能优化方法及系统

本发明涉及煤矿巷道支护，尤其涉及一种巷道过陷落柱围岩支护参数智能优化方法及系统。

背景技术：

1、陷落柱主要分布在我国华北煤田，属于煤矿地质构造。在陷落柱比较发育的地区，煤系地层中的煤层及其周围的岩石常遭到严重的破坏，给煤矿安全生产造成了很大的困难。在巷道掘进前方，受采动影响，陷落柱破碎区持续发生破坏，当掘进工作面通过陷落柱破碎区时，容易发生安全事故。在巷道过陷落柱区域采取合理的支护设计，能够解决该难题。

2、针对巷道过陷落柱围岩支护，目前主要采用工程类比方法。然而，工程类比方法并未同时考虑到围岩变形、锚杆变形、支护成本和掘进速度等指标。采取传统的工程类别方法确定支护方案，无法使技术参数和设计方案达到最优。为此，需要一种巷道过陷落柱围岩支护参数智能优化方法，来实现巷道过陷落柱围岩支护参数的智能优化，确保巷道顶板安全。

技术实现思路

1、本发明要解决的技术问题是针对上述现有技术的不足，提供一种巷道过陷落柱围岩支护参数智能优化方法，实现巷道过陷落柱围岩支护参数的智能优化。

2、为解决上述技术问题，本发明所采取的技术方案是：一方面，本发明提供一种巷道过陷落柱围岩支护参数智能优化方法，基于flac3d数值模拟软件与模糊综合评价方法实现，包括以下步骤：

3、步骤1：基于加固拱、自然平衡拱理论，初步选择巷道过陷落柱围岩的支护参数，设计三种不同支护方案；

4、步骤2：采用数值模拟软件，对三种支护方案进行数值模拟，分析巷道过陷落柱围岩顶板变形、两帮变形和锚杆变形数据，计算不同方案的支护成本，并考虑不同支护方案的影响确定掘进速度；

5、步骤3：构建巷道过陷落柱围岩支护效果评价指标体系；

6、步骤3.1：以支护方案的支护技术、经济为评价指标，综合评价巷道过陷落柱围岩控制效果；定义巷道顶板下沉率c1、帮部变形率c2、底板底鼓率c3、顶板锚杆变形率c4、帮部锚杆变形率c5、掘进速度b3和支护成本b4七种评价指标；

7、其中，巷道顶板下沉率c1为顶板变形量与巷道断面高度的百分比；帮部变形率c2为两帮变形量与巷道断面宽度的百分比；底板底鼓率c3：底板底鼓量与巷道断面高度的百分比顶板锚杆变形率c4和帮部锚杆变形率c5均为锚杆变形与锚杆长度的百分比；

8、步骤3.2：建立巷道过陷落柱围岩支护效果三层模糊因素综合评价体系；其中，第一层模糊因素集a＝{a1,a2}＝{支护技术指标，支护经济指标}；第二层模糊因素集b＝{b1,b2,b3,b4}＝{围岩变形，锚杆变形，掘进速度，支护成本}；第三层模糊因素集c＝{c1,c2,c3,c4,c5}＝{巷道顶板下沉率，帮部变形率，底板底鼓率，顶板锚杆变形率，帮部锚杆变形率}；

9、步骤4：建立巷道围岩支护效果综合评价指标分级标准并进行评价指标模糊量化，求解各指标各级的隶属度函数；

10、步骤4.1：建立评价指标分级标准，对围岩支护效果进行等级划分；

11、采用10分制度对巷道过陷落柱围岩支护效果评价进行模糊综合评价打分，构建模糊综合评价打分表；

12、步骤4.2：采用半梯形隶属度函数来描述各项指标，对评价指标进行模糊量化处理；

13、步骤4.3：将模糊综合评价打分表中各评价指标分级的数值带入半梯形隶属度函数，从而求解各评价指标的5级隶属度函数；

14、步骤5：求解巷道过陷落柱围岩支护效果评价模糊矩阵；

15、步骤5.1：整理步骤2中三种支护方案的7种支护评价指标数值；

16、步骤5.2：将各支护评价指标数值依次带入对应的隶属度函数进而求解出围岩变形、锚杆变形和支护经济评价模糊矩阵；

17、步骤6：求解巷道过陷落柱围岩支护效果各评价指标的ahp法权重和熵权法权重，提高权重的客观准确性，并将ahp-熵权进行组合赋权；

18、步骤6.1：求解巷道过陷落柱围岩支护效果评价指标ahp法权重；

19、步骤6.1.1：确定各级评价指标的重要程度；

20、确定一级评价指标重要程度为a1＞a2；二级评价指标重要程度依次为b1＞b2＞b3＞b4；三级评价指标中围岩变形部分重要程度依次为c1＞c2＞c3；锚杆变形部分重要程度依次为c4＞c5；

21、步骤6.1.2：对各评价指标进行重要程度比较，引入1～9级标度，确定围岩支护效果评价指标的一级～三级判断矩阵；

22、步骤6.1.3：将各级判断矩阵进行归一化、最大特征值和一致性检验，求解得到各支护效果评价指标的ahp法权重；

23、步骤6.2：求解巷道过陷落柱围岩支护效果评价指标熵权法权重；

24、步骤6.2.1：将步骤6.1.2的各级判断矩阵进行归一化处理；

25、步骤6.2.2：求解一级～三级评价指标的信息熵；

26、步骤6.2.3：确定一级～三级评价指标的熵权法权重；

27、步骤6.3：采用ahp-熵权组合进行权重求解；利用平均值法将ahp法与熵权法得到的权重进行耦合，从而得到综合权重：

28、

29、式中：wj为综合权重；为ahp法权重；为熵权法权重；α，β均为权重系数；

30、步骤7：求解一级、二级评价指标的隶属度函数，进而求解目标层隶属度函数；

31、步骤7.1：求解二级指标隶属度函数；将三级指标c1～c3的综合权重与围岩变形模糊矩阵相乘得到二级指标b1的隶属度函数；将三级指标c4～c5的综合权重与锚杆变形模糊矩阵相乘得到二级指标b2的隶属度函数；

32、步骤7.2：求解一级指标隶属度函数；将一级指标a1的综合权重与二级指标b1、b2的隶属度函数相乘得到一级指标a1的隶属度函数；将二级指标b3、b4的综合权重与支护经济模糊矩阵相乘得到一级指标a2隶属度函数；

33、步骤7.3：求解目标层隶属度函数；将一级指标a1和a2的综合权重与一级指标隶属度函数相乘得到目标层隶属度函数；

34、步骤8：进行支护方案模糊综合打分；将目标层隶属度函数与1-5级别的模糊综合评价打分值进行相乘，得到方案ⅰ、ⅱ和ⅲ的模糊综合打分；进而根据各方案模糊综合打分分值大小确定巷道过陷落柱围岩支护最优参数。

35、另一方面，本发明还提供一种巷道过陷落柱围岩支护参数智能优化系统，该系统界面基于node-red语言开发，flac3d软件为该系统的底层服务，界面模块化输入支护参数并储存在云数据库mysql中；界面化输入支护设计参数，自动生成dat文件，采用python语言调用dat文件，dos模式下运行flac3d软件，界面模块化展示数值计算结果；采用python语言将巷道过陷落柱围岩支护效果评价的指标体系和各评价指标的ahp-熵权组合权重嵌入巷道过陷落柱围岩支护参数优化系统；最终，基于综合评价打分，确定最优支护参数。

36、采用上述技术方案所产生的有益效果在于：本发明提供的一种巷道过陷落柱围岩支护参数智能优化方法及系统，能够选取巷道过陷落柱围岩支护最优参数，同时满足支护技术指标和支护成本最优，保证巷道过陷落柱围岩安全、高效作业。