煤层底板突水评价的新型实用方法-脆弱性指数法 所属技术领域
本发明涉及一种煤矿山煤层底板突水的评价预测方法,尤其是一种煤层底板突水评价预测的“脆弱性指数”新型实用方法。
背景技术
煤层底板突水是人为采掘工程活动导致煤层围岩体应力场能量释放、煤层底板隔水岩体结构破坏、矿井局部充水水文地质条件突变的一种地下岩体失稳现象,是一种受控于多因素影响且具有非常复杂形成机理的非线性动力现象。由于现代数学还未能用一个确定性数学方程详实地描述这种机理复杂且受控因素颇多的非线性动力现象,因而经验、半经验或统计等模型是刻画这种非线性动力现象的主要方法,如1964年焦作矿区水文地质大会借鉴匈牙利韦格弗伦斯相对系数概念(隔水层厚度与水压力之比)提出的突水系数评价法,由于该方法物理概念简单,计算方便,现场易于操作,因而一直是我国煤层底板突水评价的主要方法。但是,该方法仅考虑充水含水层水压和煤层底板隔水岩段厚度这2个控制底板突水的影响因素,虽经几次修改完善,但仍未能摆脱仅考虑这2个控制因素且无影响“权重”概念等重大缺陷。显然,该方法所能考虑的突水控制因素极为有限,未能描述煤层底板突水这种受控于多因素且具有非常复杂机理的非线性动力现象,歪曲了煤层底板十分复杂的突水机理,已不能适应新的采矿方法和新的地质环境条件下的煤层底板突水评价。
为了解决煤层底板突水预测预报评价难题,本发明人于2007年曾系统地建立了能够全面真实反映煤层底板复杂突水机理和演变过程的底板突水的主控指标体系,详细论述了体系中各突水主控因素在矿井突水过程中的作用方式与特征。但是,这种详实系统地分析确定控制煤层底板突水的主控因素,正确建立煤层底板突水的物理概念模型,仅仅是完成了解决底板突水预测预报评价难题的第一步。
长期以来,煤炭工业迫切需要系统地解决煤层底板突水预测预报评价难题,对煤层底板突水做出科学的区划和预测预报评价。
经过长期的研究和实践,本发明满足了上述煤炭工业的需求。
【发明内容】
本发明目的之一在于,满足煤炭工业的需求,提供一种能够真实描述受控于多因素影响且具有非常复杂形成机理的非线性动力现象的煤层底板突水的数学模型和评价方法,完成系统地解决煤层底板突水预测预报评价难题的第二步。
本发明进一步的目的在于,满足煤炭工业的需求,提供一种在现场工程实际应用中能够概念清晰明确,计算程序简单实用,易于现场工程技术人员掌握,操作便利的现代先进的数学模型和评价方法,完成系统地解决煤层底板突水预测预报评价难题的关键的第三步。
为了实现上述目的,本发明采用了以下技术方案:煤层底板突水评价的新型实用方法-脆弱性指数法,是可以考虑煤层底板突水众多主控因素相互复杂作用关系和相对“权重”比例并具有多级分区特点的一种评价方法,包括以下步骤:
1.基于多源信息融合理论,以GIS为操作平台,确定煤层底板突水主控因素,经过数据采集、分析和处理,建立各主控因素子专题层图;
2.应用多源地学数据复合叠加原理,采用现代线性或非线性数学方法,通过模型的反演识别或学习训练,确定出煤层底板突水的各主控因素对复杂突水过程的“贡献”或“权重”,建立煤层底板突水预测预报评价模型;
3.根据研究区各单元计算的突水脆弱性指数频率直方图分析,合理确定突水脆弱性分区阈值,最终对煤层底板突水脆弱性做出科学的区划和预测预报评价。
由于采用了上述的“脆弱性指数法”技术方案,本发明具有的有益效果在于:
1.由于现代科学技术尚未提出一个能够详实描述煤层底板这种复杂突水机理和受控因素颇多的非线性动力现象的确定性数学方程与模型,因而应用现代多源信息融合理论,采用具有强大空间数据统计分析处理功能的地理信息系统(GIS)与线性或非线性数学的耦合方法,提出的煤层底板突水评价的脆弱性指数法,较好地解决了煤层底板突水预测预报评价的难题。
2.突破了煤层底板突水评价传统的突水系数法仅能考虑两个控制因素且无影响“权重”概念等重大缺陷的束缚,采用GIS与现代非线性数学的耦合方法,真实地刻画了煤层底板突水这种受控于多因素且具有非常复杂形成机理的非线性动力过程。
【附图说明】
附图1,本发明煤层底板突水评价地新型实用方法-脆弱性指数法工作流程图。
具体实施例
下面将结合附图对本发明作详细描述。
本发明一种煤层底板突水评价的新型实用方法-脆弱性指数法的具体实施技术方案包括以下步骤:
1.建立各主控因素子专题层图:
基于多源信息集成理论,以GIS为操作平台,确定煤层底板突水主控因素,经过数据采集、分析和处理,建立能够全面真实反映煤层底板复杂突水机理和演变过程的底板突水的主控指标体系,详细论述了体系中各突水主控因素在矿井突水过程中的作用方式与特征。
2.建立煤层底板突水预测预报评价模型:
应用多源地学数据复合叠加原理,采用地理信息系统(GIS)与线性或非线性数学方法的耦合技术,通过模型的反演识别或学习训练,确定出煤层底板突水的各主控因素对复杂突水过程的“贡献”或“权重”,建立煤层底板突水预测预报评价模型。
所述的采用地理信息系统(GIS)与线性或非线性数学方法的耦合,可以是地理信息系统(GIS)与人工神经网络(Artificial Neural Network,ANN)的耦合。煤层底板突水是一个不断与外部环境进行物质、能量、信息交换的开放系统,具有协同性、自组织性、主控因素多、条件不准确的特点,突水的基础空间信息量大。人工神经网络(ANN)具有大规模并行、分布式存储与处理、自组织和自适应能力,能够同时处理复杂的非线性系统中多种主控因素在控制机理不清楚、影响参数模糊状况下彼此之间的相互作用关系和相对权重,但无法解决煤层底板突水的空间信息分析问题;而地理信息系统(GIS)则具有强大的空间信息处理分析和统计功能。将GIS与ANN耦合,就是整合他们各自的优点,以GIS的空间数据统计分析功能为操作平台,利用它强大的空间信息处理分析能力,对各种图形信息进行量化,建立煤层底板突水各主控因素的子专题层图,并对它们实施复合叠加处理,确定研究区计算单元;通过ANN对采集到的有关煤层底板突水样本数据的反复训练学习,确定各主控因素在复杂的底板突水非线性系统中的“贡献”或“权重”。GIS与ANN耦合技术的应用,建立了煤层底板突水的评价模型,经过模型验证,即可确定煤层底板突水的最终预测模型,计算煤层底板突水的脆弱性指数。
所述的采用地理信息系统(GIS)与线性或非线性数学方法的耦合,也可以是GIS与证据权法耦合。证据权法实现方法:①将研究区内已知突水点根据证据权方法按照一定面积将研究区进行网格划分,生成具有面积属性的训练图层或底板突水网格图层;②分析煤层底板突水的主控因素,将这些主控因素作为证据层,并对这些证据层进行预处理,将其转化为分类数据,形成证据专题图层;③将各证据因子的专题图层分别与训练图层进行叠加,计算每个证据图层的先验概率及权重或模糊权重;④对证据图层关于煤层底板突水条件的独立性进行检验,并根据前验概率及权重,筛选出最合理的证据因子专题图层,进行后验概率计算;⑤根据后验概率计算的结果确定煤层底板危险突水区。
基于GIS的煤层底板突水预测方法分为数据驱动和知识驱动。证据权重法属于数据驱动方法,数据驱动避免了权值选择的主观性,易于编程,易于组合多图件模式,但是用于分析的证据图层必须是条件独立的。一方面GIS具有强大空间信息处理与分析和各种图形信息量化功能,可以解决证据权重法实现过程中训练图层和证据图层的生成问题,生成各种煤层底板突水预测所需的信息图层,建立煤层底板突水各主要控制因素的子专题层图,而且可将预测评价的突水脆弱性指数的后验概率图层按照多种方法进行脆弱性分区;另一方面,证据权重方法的多源信息集成和综合能力,极大提高了GIS的空间信息分析和处理功能。因此,将证据权重法与GIS耦合应用,提高了受控多因素影响的煤层底板突水预测评价的可靠度。
所述的采用地理信息系统(GIS)与线性或非线性数学方法的耦合,还可以是GIS与Logistic回归法耦合。GIS中的Logistic回归法首先将各个证据层空间叠加处理,确定出互不相交且互不相同的突水预测区域。假设将m个证据层进行空间叠加,共产生n个互不相交且不同的子区域,有n1区域突水点和n2区域非突水点,n=n1+n2,设仅前n1个区域有突水点,计算结果见表1。
表1多证据层空间叠加结果
区域 序号 常 数 项 证据层X1 在每个子 区域内的 取值 证据层X2 在每个子 区域内的 取值... 证据层Xm 在每个子 区域内的 取值 每子区 域包含 的水点 每子区域 包含的面 积单元 子区域内是 否包含突水 点标示 1 1 x1(1) x2(1)... xm(1) m1 A1 1
区域 序号 常 数 项 证据层X1 在每个子 区域内的 取值 证据层X2 在每个子 区域内的 取值... 证据层Xm 在每个子 区域内的 取值 每子区 域包含 的水点 每子区域 包含的面 积单元 子区域内是 否包含突水 点标示 2 1 x1(2) x2(2)... xm(2) m2 A2 1 ... ... ... ...... ... ... ... ... n1 1 x1(n1) x2(n1) xm(n1) mn1 An1 1 n1+1 1 x1(n1+1) x2(n1+1) xm(n1+1) 0 An1+1 0 n1+2 1 x1(n1+2) x2(n1+2) xm(n1+m) 0 An1+2 0 ... ... ... ...... ... ... ... ... n 1 x1(n) x2(n) xm(n) 0 An 0
由表1的值计算在各个子区域的ln[P/(1-P)]的估计值。由于在第i个子区域上的概率Pi可用估计,ln[P/(1-P)]的估计:当mi=0或Ai=mi时,上式无法计算。则将上式修正为:
令则z=a0+a1x1+a2x2+…+akxk为线性回归模型,利用表1和的估计值,用多元线性回归的方法获得a1,a2,…,ak的估计值
所述的GIS与Logistic回归法耦合,可以是Logistic回归与GIS分开使用,单独发挥自己功能,首先在GIS中对证据层数据进行分类,形成多分类数据图层,在此基础上进行多证据层的空间叠加;然后根据叠加结果,运用统计分析方法,估计Logistic回归的回归系数并进行假设检验,计算在研究区内各点预测对象发生的预测概率;最后,将预测概率结果在GIS中显示、成图,并进行突水脆弱性多级分区预测。
所述的GIS与Logistic回归法耦合,还可以是将Logistic回归嵌入于GIS中,在GIS平台上进行Logistic回归的所有工作。
所述的采用地理信息系统(GIS)与线性或非线性数学方法的耦合,还可以是GIS和层次分析法耦合。即以GIS的空间数据统计分析功能为操作平台,利用GIS强大的空间信息分析处理功能,对主控因素信息进行空间复合叠加处理。应用层次分析法确定煤层底板突水各主控因素对突水的“贡献”或权重系数,建立目标评价模型.根据现场采集的大量煤层底板突水样本,经模型反演识别,建立煤层底板突水脆弱性评价模型,计算煤层底板突水脆弱性指数。依据研究区各单元计算的脆弱性指数值,通过累计频率统计直方图分析,确定煤层底板突水脆弱性分区阈值。最后运用GIS功能,提出研究区煤层底板突水脆弱性的分区方案。
3.根据研究区各单元计算的突水脆弱性指数频率直方图分析,合理确定突水脆弱性分区阈值,最终对煤层底板突水脆弱性做出科学的区划和预测预报评价。