顾及指标权重变化的露天矿边坡稳定性评价系统的设计与实现

王丽丽， 虞列沛

（广东省有色地质测绘院，广州 510080）

目前，我国矿山开采主要有地下开采和露天开采两种形式，据不完全统计，我国露天开采占矿山开采总量的87.1%，而中小型露天矿山约占94.9%[1]。 但由于露天矿边坡处于开放的外界环境中， 受自然条件、地质条件、外部条件（降雨、扰动等）以及管理水平和资金投入等各种内外部环境的影响，导致边坡发生滑坡等地质灾害，给周边人员及设备造成损失[2-4]。 因此，开展露天矿边坡稳定性分析与评价已成为国内外相关领域学者和企业管理人员研究的热点。

露天矿边坡的稳定性是多种因素耦合的结果，国内外学者及企业针对这一问题做了大量且富有成效的研究。 目前，对其稳定性评价和分析的常用方法有工程类比法[5]、极限平衡法[6-7]、数值模拟法[8-9]、相似材料模拟法[10-11]、非确定性分析法[12-13]等。 考虑露天矿边坡稳定性影响因素和评价方法的繁琐性和复杂性，计算效率较低，分析结果具有一定的滞后性。 为保证评价结果的快速、及时、准确，邬凯等以四川红层丘陵区公路边坡为例，考虑其坡体结构及失稳特征，开发了边坡稳定性在线量化评价系统[14]；黄健等在Android二次开发的基础上，采用灰色关联分析及综合指数相结合的方法，开发了一套边坡稳定性快速评价系统，并应用于安徽汤屯高速公路边坡，取得了良好的效果[15]；戴哲等以乌尉高速公路边坡为例，采用模糊综合评价与层次分析相结合的边坡稳定性评价方法，开发了基于ArcGISEngine 的快速评价系统，评价结果与实际基本吻合[16]；田维刚等开发了基于模糊综合评判的边坡稳定性评价系统，实现了对王家坪滑坡灾害危险性等级的评价[17]；陈强等通过对松散体边坡体进行调研分析，确定影响其稳定性的评价指标，开发了基于模糊分析的快速评价系统，研究表明，系统评价结果具有较高的可靠性[18]。 在上述研究结果中，尽管实现了对边坡稳定性程度和等级的快速、及时评价，但由于影响边坡稳定性状态的评价指标取值经常会发生变化，导致各指标权重也在动态变化，其评价结果的准确性存在一定的局限性。 LIN 等针对公路边坡11 个参数的评价指标体系， 提出了权重系数分布的指标重要性排序，研究表明，排序方法大大减少了评级过程中的主观性[19]。 XIA 等利用群决策理论和层次分析法确定各评价要素的指标权重分配，有效避免了单个专家主观偏见造成的评价误差[20]。QIAO 等利用可拓理论和改进的熵权法， 将评价指标从单一的确定值变为区间值，更全面地分析评价指标的影响程度[21]。

根据影响边坡稳定性评价指标的不确定性以及指标权重动态变化的特性，为解决常权重系数的确定过度依赖领域专家的主观意见和经验的问题，以越堡露天矿HP1 边坡为例， 利用在解决不确定信息或数据方面有着特有优势的未确知测度评价模型，以及在确定指标权重具有更为客观优势的信息熵法， 设计、开发露天矿边坡稳定性动态评价系统，实现边坡评价指标权重的动态计算、 边坡稳定性智能分析与评价等，以提高稳定性评价的准确性，为提高边坡灾害预警能力提供技术指导和科学依据。

1 顾及指标权重变化露天矿边坡稳定性评价方法

以越堡露天矿边坡为应用案例，对其进行稳定性程度的分析和评价。根据实际调查，2016 年10 月，该矿西部HP1 边坡在受到连续降雨的影响下， 其稳定性变差，最终导致了滑坡地质灾害的发生。

1.1 影响矿山边坡稳定性的因素分析

影响矿山边坡稳定性的因素一般有： 地质因素、水文因素、外部环境及人为因素的影响。 鉴于此，本文结合越堡露天矿边坡上述矿区基本情况及研究区调查结果，分析得到影响研究区山体稳定性的因素主要有：

1）地形地貌：项目研究区地势较为平坦开阔，但地貌较为复杂，经多年的采矿削坡之后，形成了上陡中缓下陡的形态。

2）岩土类型：根据研究区现场地质勘查和调查结果，该区域主要包含了石炭系和第四系覆盖层。

3）地质构造：研究区附近地质构造主要是北东向褶皱，其次是北西向褶皱和断裂，受此构造因素的影响，研究区还生成一些次生断裂。

4）水文因素：项目研究区地处于丘陵与河流交接处，地表面鱼塘较多，与地下水形成连通，导致经常出现地表水向矿坑底部下渗。此外，研究区雨季长、降雨量大，在雨水渗透作用下，对边坡的稳定性影响较大。

5）外部环境及人为因素的影响：研究区属于抗震烈度6 度区，外部营力及人为因素对边坡稳定性影响较小，对矿区安全生产影响不大。

1.2 露天矿边坡稳定性评价指标的提取

在分析总结国内外露天矿边坡稳定性影响因素研究成果的基础上[22-23]，结合越堡露天矿HP1 边坡的特点，依据现场地质勘查和调查结果，将其危险性评定为I 级（C1）、Ⅱ级（C2）、Ⅲ级（C3）、Ⅳ级（C4） 4 个等级，分别代表发生滑坡的危险性极高、较高、一般和较低， 并建立了越堡露天矿边坡稳定性评价指标体系，提取了3 个一级指标：环境条件、地质条件和工程条件；10 个二级指标： 月最大降雨量 （X1)、 水文条件（X2）、地下水体（X3）、岩体结构（X4）、地质构造（X5）、地层岩性（X6）、黏聚力（X7）、内摩擦角（X8）、边坡高度（X9）和边坡坡度（X10），其中X1、X7、X8、X9、X105 个指标为定量指标，其他5 个为定性指标，并分别对其进行分级、赋值，如表1 所列。

表1 HP1 边坡稳定性定性评价指标及其赋值Table 1 Qualitative evaluation index of HP1 slope stability and its assignment

1.3 评价指标动态变权重计算方法

权重是衡量影响因素（评价指标）对边坡稳定性程度或大小的一个重要指标。在不同影响因素的作用下，由于各评价指标权重对露天矿边坡稳定性影响程度的不同，边坡呈现出不同的危险性程度。 考虑到边坡是处于一个开放的现实空间，其影响因素赋值将发生动态变化，其指标权重不再是固定不变的。 为克服定权对评价结果产生的偏差，合理、准确地评价边坡稳定性程度，计算出评价指标的动态权重是关键。 由此，本文在依据信息熵理论计算出各评价指标的常权向量的基础上， 通过构造出变权向量的均衡函数，计算出其动态变权向量，具体步骤及方法如下：

1）在对评价指标进行归一化处理的基础上，计算出第j个评价指标Xj（j=1,2, …,m）的常权重为wj。 由此，根据信息熵基本原理：

式（1）—式（3）中：ej代表的是评价指标熵值；k表示的是标准化系数；pj表示的是传递系数；aj表示的是评价指标值；bj表示的是评价指标归一化值。但是，考虑到在计算时，可能会出现pj=0 的情况，那么lgpj则不存在，无法计算。 为保证计算的可用性，假设pj=0时，lgpj=0。 与此同时，考虑到ej→1 时，即使熵值变化一个非常小的量， 其权重也会出现巨大的不同。为此，构建了改进的指标常权重计算公式，如式（4）所示：

式（4）中：ej、el、et表示的是不同影响因子j、l和t的熵值。

2）变权向量均衡函数的建立

假设有：

S（X）=（S1（xj），S2（xj），…，Sm（xj）），j=1，2，…，m；v≥0；0＜u≤1

式（5）中：u表示否定因子，当xj≤u时，其权重值可以依据变权的方法予以增加。v表示惩罚因子，当v值增大时，其惩罚程度随之增加，其中，u和v的取值依据具体边坡情况选定。

3）变权重模型的构建

利用计算得到的评价指标常权向量w，以及建立的变权向量均衡函数S（X），再依据式（6）计算得到变权向量矩阵W。

1.4 顾及指标权重变化露天矿边坡稳定性评价方法的建立

依据文献[1]中关于未确知测度理论建模的基本知识，建立顾及指标权重变化的边坡稳定性评价模型，将其危险性评定为I 级（C1）、Ⅱ级（C2）、Ⅲ级（C3）、Ⅳ级（C4） 4 个等级，分别代表发生滑坡的危险性极高、较高、一般和较低。其建模基本流程如图1所示。

图1 顾及指标权重变化露天矿边坡稳定性评价建模流程Fig. 1 Modeling process of open pit slope stability evaluation considering the change of index weight

1）单指标测度函数及单指标评价矩阵的构建

假如μjk=μ（xj∈Ck）表示评价指标在其相应等级Ck中的值xj所属程度，而且要使得μ 满足：

如果μ 能够同时满足式（7）和式（8）条件的要求，那么μ 则称为未确知测度，式（10）中的（μjk）m×p则称为单指标测度评价矩阵。

2）评价指标动态变权重的确定

根据1.3 节中动态变权重的计算思路，并考虑边坡实际，项目在对评价指标进行归一化处理和分析的基础上，将v设定为0.5，u设定为0.3，并针对不同月最大降雨量的变化， 分别得到HP1 边坡稳定性评价指标的常权向量w和变权向量W。

3）多指标未确知综合测度μ 的计算

结合上述方法计算出的W及相应的（μjk）m×p，可以得到露天矿边坡HP1 的稳定性多指标未确知综合测度如式（11）：

那么，{μ1，μ2，…，μp}则表示为HP1 边坡稳定性的多指标测度评价向量。

4）置信度识别准则的确定

根据越堡露天矿HP1 边坡的特点， 设置信度0.5≤λ＜1，并使得：

则认为评价对象属于第k0个评价类Ck0。

5）边坡危险性等级评价

考虑到边坡稳定性评价等级为一个宏观的概念，即使评价指标取值不同， 其评价等级也可能一样，导致无法准确判断边坡的危险性状态。 因此，为定量分析不同评价指标取值情况下边坡的危险性程度，提出危险性相对重要度指标q，如式（13）。

其中，I的取值依据边坡稳定性评价等级， 并按照C1＞C2＞C3＞C4的排序， 分别赋值为I1=4，I2=3，I3=2，I4=1。

根据越堡露天矿HP1 边坡实际， 并结合上述计算流程和方法，计算得到的HP1 边坡常权、变权、多指标测度向量及评价结果在系统运行结果中展示，具体结果后文将描述。

2 顾及指标权重变化的露天矿边坡稳定性评价系统的设计与实现

2.1 设计平台及系统架构

系统开发工具： 本系统主要采用微软Visual Studio 2012 编译工具，采用Windows 窗体应用程序开发模式， 编程语言为C#4.0， 环境为.NET Framework 4.5，编程方式较为简单，灵活性、兼容性较好。编译生成的可执行文件，具有良好的可移植性，同时windows 窗体界面通俗易懂。

露天矿边坡稳定性动态评价系统主要由数据管理、数据处理和评价分析预警3 部分组成，其中，系统中采用最常见的txt 文本文件进行数据的读取和输出，能极大减小系统资源的占用空间，同时也方便对数据进行修改和管理。系统在对所采集的数据进行转换、分析、处理的基础上，采用信息熵和未确知理论，实现评价指标动态变权重的计算、边坡稳定性评价以及边坡变形预警等，其系统架构如图2 所示。

图2 顾及指标权重变化的露天矿边坡稳定性评价系统架构Fig. 2 Structure of open pit slope stability evaluation system considering the change of index weight

2.2 系统功能的设计及实现

根据越堡露天矿企业的要求，课题组开发出的顾及指标权重变化的露天矿边坡稳定性智能评价系统，尽量做到简洁、明了、适用，不仅方便企业的使用，还能指导矿区的施工和管理，实现企业的高效、安全生产。 为突出本系统特点，本文将系统部分主要功能的实现及成果分析界面展示如图3 所示。图3（a）为系统主界面，用户通过输入用户名、密码，登入系统。 图3（b）是依据未确知测度理论方法构建出的露天矿边坡各影响因素单指标测度函数， 避免了人为计算的繁杂，提高了工作效率。在此基础上，计算出HP1 边坡的多指标测度向量，如图3（c）所示，为边坡稳定性等级及其重要度的计算提供了基础。 考虑到越堡露天矿边坡受降雨等因素的影响， 其稳定性评价指标权重随之发生变化，系统依据1.2 节中评价指标动态权重确定方法， 计算出各评价指标的常权重和变权重向量，如图3（d）所示，在此基础上，再依据1.3节中露天矿边坡稳定性评价模型的建立方法，计算出边坡稳定性评价等级， 其研究结果与边坡实际一致。

图3 系统功能的实现与应用Fig. 3 Implementation and application of system functions

根据该系统的评价和分析，不仅可以提前预测边坡在不同状态下的稳定性状态，而且在一定程度上能够指导矿山企业的设计、生产和管理，保障施工设备及人员的安全。

3 结束语

1） 针对传统边坡稳定性评价指标权重受主观经验影响，通常确定其为常权向量的问题，在充分考虑评价指标权重的动态变化特性的基础上，基于信息熵方法，构造评价指标均衡函数，实现评价指标动态变权重的合理性分析。

2） 针对露天矿边坡稳定性影响因素复杂、 不确定，导致其评价准确性不足的问题，结合评价指标动态变权重模型，利用在解决不确定信息或数据方面具有特有优势的未确知测度方法，构建边坡稳定性动态评价模型，提高了评价的准确性，为其他相关领域稳定性评价提供了重要的参考。

3） 针对边坡稳定性评价模型计算复杂、 效率低的问题，开发评价智能系统，使边坡稳定性评价流程更为简单、更易操作，降低了计算工作量，提高了工作效率。通过在越堡露天矿边坡中的应用，结果表明：其评价结果与边坡实际基本一致， 具有较强的实用性、可靠性和可行性。

4） 本文所建立的系统评价模型主要是针对越堡露天矿边坡稳定性主要影响因素而言，其普适性有所欠缺。 为提高系统评价的普适性，系统应进一步考虑不同边坡的特定性，在数据更新等方面也有待进一步完善。