郭 金，顾冲时，何 菁

(1.河海大学水利水电学院，江苏 南京 210098；2.河海大学水文水资源与水利工程科学国家重点实验室，江苏 南京 210098；3.南京市水利建筑工程检测中心有限公司，江苏 南京 210036)

作为防洪体系的重要组成部分，堤防工程历来受到人们的重视。近年来，长江、黄河等主要河流流域的堤防工程规模和数量不断扩大，截至2019年底，全国已建成5级及以上江河堤防32.0万km，累计达标堤防22.7万km，达标率为71.0%[1]。堤防工程在抵御洪水灾害、保护下游人民生命财产安全方面发挥着重要作用，一旦失事将造成严重后果。因此，对堤防安全进行综合评价，确定堤防风险等级，对保证堤防安全运行，发挥堤防工程效益具有重要意义。目前，一部分学者已对堤防工程安全评价开展了研究，如王恩[2]建立了基于层次分析法的模糊综合评价模型，方法简单易操作，但是主观性较强，普适性较差。王秀杰等[3]建立了具有较强鲁棒性的突变理论与BP神经网络相结合的综合评价模型，从致灾因子、孕灾环境和承灾体3个方面构建评价指标体系，但BP神经网络不易收敛且容易陷入局部最小值。蔡新等[4]在堤防安全评价中引入灰色理论，构建了安全评价数学模型，但是灰色理论存在拟合度缺陷，且对原始数据有较高要求。

二维云模型拥有出色的处理模糊性和随机性能力，目前已在多个领域得到广泛应用[5-7]，但是二维云模型在堤防风险评价中的应用甚少。在堤防的风险评价过程中，风险评价指标通常具有显著的模糊性和随机性特点，同时又不可避免地受到主观因素的影响。因此，本文将二维云模型引入到堤防风险评价中，以期能很好地解决风险因素随机性和模糊性的问题，降低主观因素的影响，为堤防的风险评价开辟新途径。

1 组合赋权二维云模型理论

1.1 组合赋权

传统的权重计算方法可分为主观赋权法和客观赋权法两大类。主观赋权方法，如层次分析法、专家咨询法、乘积标度法等主要依赖于专家的相关工程经验和主观认知，受主观因素影响较大，不能反映信息的客观真实性。客观赋权方法，如熵权法、BP神经网络法、投影寻踪法等可充分利用相关数据所蕴含的有效信息，降低主观因素的影响，但容易忽略数据间的差异性而与实际情况不符。权重计算的合理性直接影响风险评价模型的可靠性和评价结果的准确性。本文基于改进的层次分析法与熵权法进行优化组合赋权，减小了风险因素指标分析中因人员知识结构差异等造成的影响，使权重分配更加合理。

1.1.1改进层次分析法确定主观权重

层次分析法是一种将定性指标量化的多目标决策方法，具有简单、实用、灵活的特点[8]。传统的层次标度方法常采用1～9标度，专家判断的主观因素对最后权重计算结果有较大的影响，计算过程烦琐，且需要进行一致性检验。改进层次分析法采用0～2新标度法，专家更容易对指标重要性做出判断，降低了人为主观因素对计算结果的影响，且不需要进行一致性检验，使得计算更加简便。具体流程如下。

a.建立堤防工程安全风险评价指标体系。将评价指标体系分为目标层、准则层、指标层3层单向递阶层次结构。

b.对同一层次元素进行重要性比较，构造比较矩阵U=(uij)n×n，其中uij为元素i与元素j之间的相对重要性标度，n为评价指标的个数。各评价指标的重要性标度见表1。

表1 指标相对重要性标度及关系

c.计算指标重要性排序指数：

(1)

d.根据排序指数构造判断矩阵B=(bij)n×n，各元素bij为

其中rmax=max{ri}rmin=min{ri}

e.构造最优传递矩阵H，其元素hij为

(3)

f.计算判断矩阵B的拟优一致矩阵T，其元素tij为

tij=10hij

(4)

g.计算拟优一致矩阵T的最大特征值和对应的特征向量，将特征向量归一化后即为所求各因素的权重。

1.1.2熵权法确定客观权重

熵权法作为一种客观赋权方法，依据各指标所含有效信息的多少确定对应指标的权重，可尽可能降低主观因素的影响，保证指标权重确定的合理性[8]，因此得到广泛应用。熵权法确定客观权重的公式为

(5)

1.1.3组合赋权

采用差异系数法对主客观权重进行线性组合赋权，消除主客观单独赋权对权重结果的影响以达到扬长避短的效果：

W=αW′+βW″

(6)

其中

α+β=1 (α、β∈(0,1))

式中：W′、W″分别为层次分析法和熵权法确定的指标权重向量；p1、p2、…、pn为权重W′中各评价指标按升序排列的权重值。

1.2 二维云模型

二维云模型描述了更复杂的不确定信息由定性到定量的转换，由(Ex,Enx,Hex)和(Ey,Eny,Hey)两组云模型数字特征组成。设F为服从正态分布的二维随机函数，则称由满足式(7)的云滴drop(xd,yd,μd)组成的云模型为二维正态云模型[9-10]。

式中：d为云滴序号；(xd,yd)为以(Ex,Ey)为期望，(Enx,Eny)为标准差生成的随机数；(Pxd,Pyd)为以(Enx,Eny)为期望，以(Hex,Hey)为标准差生成的随机数；μd为隶属度。

2 堤防工程风险评价模型

2.1 风险评价指标体系的确定

堤防工程风险评价指标体系的建立是对堤防工程进行风险评价的基础。大多数的堤防风险评价方法中评价指标的确定依赖于主观经验判断，脱离了现行设计规范或评价准则，影响堤防评价的全面性和结果的科学准确性。评价指标体系的确定应遵循科学性原则、层次性原则、系统性原则、时效性原则、可行性与实用性原则。本文依据SL/Z 67—2015《堤防工程安全评价导则》《堤防工程设计规范》，同时参考工程实际情况及相关文献[11-14]，建立科学完整的评价指标体系。以堤防工程风险评价为总目标，一级指标包括运行管理、堤基质量、堤身质量、防洪能力、渗流安全和结构安全6个指标，一级指标再具体划分为20个二级指标，如图1所示。

图1 堤防工程风险评价指标体系

2.2 风险云

事故发生的后果和概率共同决定了事故的风险等级，故对堤防风险评价指标的发生概率和后果进行分析可确定提防的风险等级，即可明确堤防的安全状态。邀请业内专家对事故发生概率和事故后果进行打分，通过分析打分结果确定事故的等级[15]。风险云的计算公式为

(8)

式中：q为专家人数，即样本数量；xt为第t位专家的赋分值。

2.3 标准云

SL/Z 67—2015《堤防工程安全评价导则》中，堤防风险评价在结合运行管理、工程质量评价和防洪标准、渗流安全、结构安全性复核的基础上，将堤防工程的安全等级分为安全、基本安全、不安全三类。考虑到工程实际情况复杂，且存在诸多不确定性因素，故将堤防的安全等级做更细致的划分，以更准确地分析堤防的安全状况。最终拟定堤防的风险评价等级标准为高风险、较高风险、较低风险、低风险。参考我国水库大坝风险标准[16]，同时结合堤防工程实际运行条件确定堤防工程各风险等级对应风险概率。具体等级划分等见表2。

表2 风险等级评价标准

标准云的数字特征按下式计算：

(9)

式中Smin、Smax分别为风险等级标准区间的最小值和最大值。

2.4 综合风险云

各级指标权重确定后，按照从下到上的顺序计算综合风险云，即由二级指标权重和风险云矩阵计算1级风险云，由一级指标权重和风险云矩阵计算综合风险云：

(10)

3 实例验证

3.1 工程概况

某堤防工程位于黄河下游平原地区，为4级堤防工程，设计水位高于流域平原区的平均地面高程，对流域平原区有重要的防洪作用。堤防工程防洪标准为20年一遇，防洪保护面积为130 km2，其所在河段为土质岸坡，河岸遭受较严重冲刷破坏，堤身单薄。堤身横断面为梯形，堤顶宽度为4.0 m，均质土堤。堤基为多层结构，上部为淤泥质土，中下部为淤泥、含贝壳淤泥和黏土，下部为基岩。

3.2 权重计算

邀请6位业内具有丰富经验的专家对堤防工程风险评价指标进行评分，评分依据见表3；各指标的发生概率等级和事故后果等级评分结果见表4。

表3 风险等级评分依据

表4 后果等级和概率等级专家评分

根据式(1)～(4)确定主观权重，式(5)确定客观权重，式(6)进行主客观权重组合赋权，根据式(8)～(10)计算各评价指标的风险云数字特征。各指标的权重和风险云数字特征见表5。

表5 后果等级和概率等级风险云数字特征

3.3 结果分析

运用正向云发生器可得到综合风险云与标准风险云的对比，如图2所示。为明确运行管理、堤基质量、堤身质量、防洪能力、渗流安全和结构安全6个1级指标的风险等级，得出各1级指标风险云与标准云的对比，如图3所示。通过分析图2和图3，可确定该堤防工程总体风险等级为Ⅱ级，即具有较高风险；6个1级指标中，运行管理B1、堤基质量B2、堤身质量B3、防洪能力B4和渗流安全B5均处于Ⅱ级风险附近，结构安全B6风险等级更高，与Ⅰ级风险更加接近。实际情况中，由于该堤防工程建成已久，长期受水流的冲刷使堤身变得单薄，历年出现过塌岸且防洪墙有明显的变形和裂缝，已进行了除险加固设计并准备实施。评价结果与该工程实际情况基本吻合，表明该方法具有一定的可操作性和准确性。

图2 综合风险云

图3 1级风险云

4 结 语

针对传统堤防工程风险评价方法的局限与不足，本文采用以二维云模型为基础的堤防工程风险评价方法，从失事后果严重程度和发生概率2个角度对评价指标的随机性和模糊性进行综合处理，实现了定性、定量之间的转换；采用改进层次分析法与熵权法组合赋权方法，充分考虑了不同风险因素、人为主观认知对风险大小的影响程度，兼顾主观因素和客观因素，保证风险评价结果更加可靠。风险评价等级以风险云图的形式呈现，评价结果更加直观。基于此方法对黄河下游平原某堤防工程进行了风险评价，评价结果与该工程实际情况基本吻合，表明该方法具有一定的可操作性和准确性。