[发明专利]一种松散承压含水层区域富水性动态确定方法

权利要求书

1.一种松散承压含水层区域富水性动态确定方法，其特征在于，所述方法包括：

步骤A：根据《煤矿防治水细则》附录一，确定松散承压含水层富水性的分区类别；

步骤B：归纳影响松散承压含水层富水性所属分区类别的主控因素，并建立层次结构模型；

步骤C：对松散承压含水层抽放水试验孔的主控因素数据进行归一化处理，并计算主观权重向量W₁；

步骤D：计算两两主控因素间的相关系数、各主控因素标准差，并计算客观权重向量W₂；

步骤E：建立针对所述主观权重向量以及所述客观权重向量的最优化一阶导数矩阵，确定优化系数，得到主客观耦合权重向量W；

步骤F：归一化后的松散承压含水层抽放水试验孔数据与主客观耦合权重向量相乘得到富水性值V_k，对富水性值V_k聚类分析后确定富水性分区类别的中断值；

步骤G：根据待评价含水层的地质勘探钻孔主控因素数据计算富水性值V_k，进而根据中断值评价该含水层的富水性。

2.根据权利要求1所述的一种松散承压含水层区域富水性动态确定方法，其特征在于，所述步骤A：包括：

根据《煤矿防治水细则》中附录一，将含水层富水性划分为弱富水性、中等富水性、强富水性以及极强富水性。

3.根据权利要求1所述的一种松散承压含水层区域富水性动态确定方法，其特征在于，所述主控因素包括：水头系数、级配系数、沉积微相系数、含水层厚度、最厚砂层厚度以及隔水系数。

4.根据权利要求3所述的一种松散承压含水层区域富水性动态确定方法，其特征在于，所述步骤C：包括：

C1：采用极大值法对最厚砂层厚度、含水层厚度、水头系数、级配系数进行归一化处理，采用极小值法对沉积微相系数、隔水系数进行归一化处理；

C2：运用Excel软件中CORREL函数对各主控因素和评价标准值进行回归分析，计算两者的CORREL相关系数C_i；根据CORREL相关系数C_i的大小，按预设的期望标度法对各主控因素进行重要程度标度；

C3：计算第一判断矩阵未知元素的值，构建出第一判断矩阵B₁；

C4：计算各主控因素的样本标准差S(i)，并依据所述样本标准差计算相对重要程度参数值；

C5：根据所述相对重要程度参数值构建第二判断矩阵B₂；

C6：分别计算第一判断矩阵对应的第一权重向量W₀₁以及第二判断矩阵对应的第二权重向量W₀₂，根据第一权重向量与第二权重向量的平均值获取主观权重向量W₁。

5.根据权利要求4所述的一种松散承压含水层区域富水性动态确定方法，其特征在于，所述C4中的依据所述样本标准差计算相对重要程度参数值步骤，包括：

利用公式，计算相对重要程度参数值，其中，b_ij为主控因素i相对于主控因素j的相对重要程度参数；s(i)为主控因素i的样本标准差；s(j)为主控因素j的样本标准差；b_m＝min{9,int[S_max/S_min+0.5]}，且int为取整函数；s_max为主控因素i的样本标准差最大值；s_min为主控因素i的样本标准差最小值。