一种基于层次分析的砂岩赋水性评价方法与流程

本发明涉及预防煤系顶板砂岩突水
技术领域：
，具体涉及一种基于层次分析的砂岩赋水性评价方法。
背景技术：
：我国煤矿水文地质条件异常复杂，煤矿开采过程中的水害事故制约着煤炭工业的发展。砂岩裂隙水突水常常会对巷道掘进和工作面回采造成影响，是影响生产效率与生产成本的重要因素。目前对砂岩赋水性评价的方法主要有可控源音频大地电磁法、瞬变电磁法、抽、注水试验法等，这些方法大多操作较为繁琐，实施周期较长，本发明针对这一问题提供一种操作简单，且可以定性与定量评价砂岩赋水性的方法。技术实现要素：1.本发明的目的鉴于此，本发明针对煤系顶板砂岩突水技术的不足，提供一种定性与定量评价砂岩赋水性的方法。2.本发明的技术方案为实现上述目的，本发明一种基于层次分析的砂岩赋水性评价方法，该方法包括：步骤a，归纳影响砂岩赋水性的主要因素，建立层次结构模型；步骤b，依托层次结构模型，分别建立因素集与权重集；步骤c，根据开采经验结合生产实践，建立评价等级集；步骤d，确定各影响因素评价指标隶属度及模糊关系矩阵；步骤e，依据层次分析结果，结合评价等级集，综合评价砂岩赋水性。本发明基于层次分析理论，结合开采经验与生产实践，定性与定量相结合的评价砂岩赋水性，评价结果合理可靠，具有一定的新颖性，且本发明操作简单，易于实际应用，为砂岩赋水性评价提供了一种新的方法和思路。附图说明此处用来说明的附图是为了对本发明的进一步解释和说明，为本申请的一部分，但并不能限定本发明。图1为本发明基于层次分析的砂岩赋水性评价方法流程图。图2为本发明实例砂岩赋水性评价层次分析结构图。具体实施方案下面结合附图及本发明所应用的实例，对本发明进行进一步说明。图1为本发明基于层次分析的砂岩赋水性评价方法流程图。如图1所示，该方法包括：步骤a，归纳影响砂岩赋水性的主要因素，建立层次结构模型；步骤b，依托层次结构模型，分别建立因素集与权重集；步骤c，根据开采经验结合生产实践，建立评价等级集；步骤d，确定各影响因素评价指标隶属度及模糊关系矩阵；步骤e，依据层次分析结果，结合评价等级集，综合评价砂岩赋水性。以下结合具体实例，对上述每一步骤进行解释说明。步骤a，归纳影响砂岩赋水性的主要因素，建立层次结构模型。在本发明实例中，结合某矿生产实践，影响煤层顶板砂岩赋水性的主要因素归纳为3个方面：砂岩特征；赋水性探测；断裂构造特征。同时，这3个主要因素又分别受2～3个下一层次因素的影响，见图2。在图2中，最高层(目标层)为煤层顶板砂岩赋水性评价(a)；中间层(b)是3个因素，砂岩特征(b1)、赋水性探测(b2)、断裂构造特征(b3)；最底层(c)是8个因素，即砂岩厚度(c1)、砂岩百分比(c2)、砂岩裂隙度(c3)、地球物理探测(c4)、钻孔探测(c5)、断层分维值(c6)、导水裂缝带(c7)、断裂切割砂岩(c8)。步骤b，依托层次结构模型，分别建立因素集与权重集。本发明实例中，鉴于步骤a中分析，评价因素集u＝{u1，u2，u3}＝{砂岩特征、赋水性探测、断裂构造特征}。由于影响煤层顶板砂岩赋水性因素较复杂，层次多，故采用层次分析确定权值。首先建立层次结构分析模型，然后进行进行各层次各因素间比较，建立判断矩阵，完成各个因素的相对重要性排序，得到4个判断矩阵(表1、2)。表1一级判断矩阵ab1b2b3b11b12b13b2b211b23b3b31b321表2指标bi的判断矩阵bic1c2c3c11c12c13c2c211c23c3c31c321i、j的标识即对于a来说，是bi对bj的相对重要性的数值，bij＝bi/bj，通常使用1，2，3，…9及其倒数；aij＝1，3，5，7，9分别表示i因素比j因素同等重要、略为重要、比较重要、非常重要、绝对重要；aij＝2，4，6，8表示i因素比j因素重要性在1，3，5，7，9之间。对图1的最终计算结果为：b层次对a层次的权重矩阵为a，c层次对b层次的权重矩阵分别为a1、a2、a3。其数值为：a＝[b1，b2，b3]＝[0.42，0.28，0.30]a1＝[c1，c2，c3]＝[0.34，0.32，0.34]a2＝[c4，c5]＝[0.17，0.83]a3＝[c6，c7，c8]＝[0.11，0.63，0.26]上述计算出的层次单排序和总排序的一致性指标及随权一致性比率均小于0.1，所以，判断矩阵和确定的权值具有一致性。计算各层次因素的组合权重，见表3。表3某矿砂岩赋水性评价因素组合权重计算结果步骤c，根据开采经验结合生产实践，建立评价等级集。根据我国煤矿开采经验结合某矿生产实践，将煤层顶板砂岩赋水性分为n级，即v＝{v1，v2，v3，v4}＝{安全性好，安全性较好，安全性差，安全性极差}。其中元素vj(j＝1,2，…m)是若干可能作出的评判结果。评价目的就在于通过对评判对象综合考虑所影响因素，能够从评价集v中获得一个最佳的评价结果。步骤d，确定各影响因素评价指标隶属度及模糊关系矩阵。本发明实例中，对于只能定性描述的评价因素指标，采用多相模糊统计方法来构造隶属函数，对于可以定量描述的评价因素指标，则采用待定系数构造隶属函数。根据南二采区顶板砂岩各评价因素指标的具体取值，分别代入各自的隶属函数，计算出各评价因素指标的隶属度值，见表4；单因素隶属度评价值见表5。表4某矿南二采区11-2煤层顶板砂岩赋水性评价因素隶属度计算结果c层b1b2b3隶属度c10.6000.42c20.3000.26c32.44000.28c404.000.10c层b1b2b3隶属度c500.3800.69c6001.460.09c70010.00.38c8006.40.24表5单因素隶属度评价值根据上述求出的各个评价因素指标隶属度值，确定构造因素c1～c3对指标b1的评价矩阵r1，因素c4～c5对指标b2的评价矩阵r2，因素c6～c8对指标b3的评价矩阵r3，即：结合上述，建立各评价因素的权重集，用c层次中的因素c1～c10对b层次中的指标b1，b2，b3，进行综合评价，评价矩阵如下：b层次指标对目标层次指标a的评价矩阵为:步骤e，按照评价等级集，对砂岩赋水性进行评价。在本发明实例中，把上述a层次的权重矩阵a及评价矩阵r代入综合评判式b＝a·r进行计算:由综合评判结果b可知，考虑最大原则，b中最大值为0.4277，该值在评判结果中的位置对应评价等级集为“较差”，即南二采区煤层顶板砂岩赋水性较差。本发明基于层次分析理论，结合开采经验与生产实践，定性与定量相结合的评价砂岩赋水性，评价结果合理可靠，具有一定的新颖性，且本发明操作简单，易于实际应用，为砂岩赋水性评价提供了一种新的方法和思路。以上所述具体实例，对本发明的目的，过程和有益效果进行详细说明，并不用于限定本发明的限定范围，凡在本发明的精神原则之内，所做的任何修改、等同替换等，均应包含在本发明的保护范围之内。当前第1页12