基于连续演化半偏减法集对势的区域旱灾风险评估方法.pdf

收稿日期：基金项目：国家自然科学基金资助项目（，）；安徽省自然科学基金资助项目（，）；安徽省高等科研计划重点项目（）作者简介：金菊良（），男，江苏苏州人，教授，主要从事结构水资源学研究工作 通信作者：崔毅（），男，安徽铜陵人，讲师，博士，研究方向为水资源系统工程：【“黄河流域生态保护和高质量发展”专题】基于连续演化半偏减法集对势的区域旱灾风险评估方法金菊良，周 乐，崔 毅，汪 洁，吴啸骎，张诗琪，赵齐雅（合肥工业大学 土木与水利工程学院，安徽 合肥；合肥工业大学 水资源与环境系统工程研究所，安徽 合肥）摘 要：为增强旱灾风险评估的可靠性，从微观层次考虑联系数分量的迁移转化过程，更细致地刻画连续半偏减法集对势中确定性与不确定性之间的关系，修正偏同势和偏反势转移信息，提出一种演化半偏联系数，构建了基于连续演化半偏减法集对势的区域旱灾风险评估模型，并应用于济南市旱灾风险评估。结果表明：修正后的连续半偏减法集对势与五元减法集对势对济南市旱灾风险和各子系统评估结果的变化趋势均一致；与五元减法集对势相比，连续演化半偏减法集对势变化较为缓和，计算出的各子系统指标值集对势变化趋势差异更为显著，便于识别影响子系统风险发展趋势的关键因子。关键词：旱灾风险评估；联系数；半偏减法集对势；因子诊断；济南市中图分类号：文献标志码：引用格式：金菊良，周乐，崔毅，等基于连续演化半偏减法集对势的区域旱灾风险评估方法人民黄河，（）：，（，；，）：，：；我国地处东亚季风区，多年来旱灾频发，其中严重旱灾 发生 次。随着经济社会发展水平不断提高，我国水资源紧缺形势日趋严重。旱灾风险反映了不同程度的干旱事件在未来可能出现的概率与造成的承灾体损失之间的关系。金菊良等提出干旱灾害损失风险的基本过程是“危险性脆弱性损失风险”的三元链式结构。旱灾风险评估可为防灾减灾工作提供科学依据，对有效指导旱灾防御工作具有重要意义。目前，旱灾风险评估方法主要有 大类：）根据旱灾风险影响因子构成，确定旱灾风险评价指标体系，建立旱灾风险指数模型。如集对分析，在此基础上，王雅燕等建立了安徽省的农业旱灾风险评价模型，其评价结果符合现实状况。）分析历史旱灾损失频率，建立旱灾风险曲线模型。如信息扩散理论，刘亚彬等利用信息扩散理论和水旱成灾面积指数进行风险概率计算，其优势在于能够弥补小样本信息的缺失，挖掘潜在信息，提高计算准确度。）模拟旱灾损失发生过程的旱灾损失风险曲线方法。如第 卷第 期 人 民 黄 河 ，年 月 ，脆弱性曲线法，徐昆等采用 曲线拟合不同情景下的干旱灾损敏感性曲线，测算旱情对玉米产量的影响程度。其中，针对不确定性和复杂性的典型旱灾风险评估问题，集对分析法以联系数从同异反 个方面，深入探讨评价样本与等级标准间的不确定性关系。联系数方法客观、计算方便，可对问题进行全面分析，具有很好的适用性。当前研究已从多个角度对旱灾风险进行评估，但对旱灾风险子系统的发展趋势和影响因素的诊断方法仍显不足，因此金菊良等提出了减法集对势、半偏减法集对势方法。为进一步验证和提高连续半偏减法集对势评估区域旱灾风险的准确性，本文从微观层面完善了半偏联系数公式，提出演化半偏联系数和改进连续半偏减法集对势，以济南市为研究区域，对比五元减法集对势，应用连续演化半偏减法集对势进行区域旱灾风险评估，识别旱灾风险发展趋势和诊断脆弱性指标，以期为旱灾风险评估提供借鉴。基于连续演化半偏减法集对势的区域旱灾风险评估模型构建基于连续演化半偏减法集对势的区域旱灾风险评估模型需通过以下 个步骤来构建。）建立风险评估指标体系。本研究将旱灾风险评估系统分为敏感性、适应性、危险性 个子系统，基于区域干旱致灾因子，建立区域干旱灾害风险评估指标体系，；，其中为评估指标体系中第 个样本值的第 个评价指标、为样本数目、为评估指标数目。基于科学性、相对性，结合研究区域实际情况，建立指标评价等级标准，；，其中 为评估标准等级数，。本文把区域旱灾风险分为轻度、较轻、中度、较严重、重度个级别。）计算样本值联系数。为刻画样本值与评价等级标准之间关系程度的可变模糊性，参考文献改进的分段函数，首先计算单指标样本 与 个等级的联系数，然后计算样本 的指标值联系数。（正向指标 或反向指标）（）（）（正向指标 或反向指标）（正向指标 或反向指标）（）（）（）（正向指标 或反向指标）（正向指标 或反向指标）（）（）（正向指标 或反向指标）（正向指标 或反向指标）（）（）（）（正向指标 或反向指标）（正向指标 或反向指标）（）（）（正向指标 或反向指标）（）（正向指标 或反向指标）（）（）（正向指标 或反向指标）（正向指标 或反向指标）（）（）（正向指标 或反向指标）（）（正向指标 或反向指标）（）（）（正向指标 或反向指标）（正向指标 或反向指标）（）式中：、分别为旱灾风险评价等级标准级的阈值。根据式（）式（）得到的联系数，计算评估样本值 隶属于评估标准等级 的相对隶属度：（，）（）将相对隶属度归一化得单指标值联系数分量：（）由式（）得到五元联系数分量，分别作为同一度，差异度 个子项、和对立度，于是可得评估指标值联系数，：（）人 民 黄 河 年第 期式中：、分别为 个差异度子项的系数；为对立度系数，一般取。同理，由各评估指标权重 可得样本 的指标值联系数：（）用级别特征值法，估计样本 指标值联系数所对应的风险等级：（）（）构建连续演化半偏减法集对势，并计算旱灾评估样本 各联系数集对势值和诊断识别评估指标。文献通过将三元联系数中的不确定项进一步细化，定义了新的五元减法集对势计算公式：（）（）（）（）（）（）式中：、为联系数分量；、分别为差异度、的系数，按经验分别取 、，可得改进的五元减法集对势公式：（）（）（）（）（）（）（）的差异度子项、分别对联系数确定性整体取值有偏同和偏反的影响，的影响是相对最不确定的。因此，在式（）中系数的取值上，考虑 取、和 取各自区间的平均值，但这是设定的一种理想取值，个不确定项具体偏同和偏反的趋势是不确定的。文献在五元半偏减法集对势的基础上，以 （）、（）的比例趋势，利用半偏联系数将 分配到偏同势 和偏反势 上，得到新的偏同势和偏反势为（）（）（）文献运用半偏联系数思想将 和 进行第二次半偏，考虑将偏同项进一步偏向同一项、偏反项进一步偏向对立度 的可能，即 和 分别以 （）和 （）的比例分到 和 上，得到连续半偏减法集对势公式为（）（）（）（）（）（）（）（）式（）利用半偏联系数的思想实现了 向 和转化，但作为新的偏同势和偏反势，需考虑不确定项对目标项分量转移的所有可能性，即 不只从 本身和 转化而来，也可由 迁移而来，同理 还应考虑偏向 的可能。同时，当不相邻的两项彼此进行半偏转移时，往往是以相邻项为媒介逐步完成的，在这一过程中，每一次半偏都会对分量产生影响，如计算 偏向 的分量时，会由 向 转化得到 （），而这一部分分量应与 协同进行第二次半偏迁移到，（）可看作是由上一步半偏演化得到的新，故在下一步向 半偏时，作为演化来源的分母应该是 本身及上一步经半偏演化组成的新分量，由此可将 （）看作一个新的联系数伴随函数，称为演化半偏联系数，将其应用于非相邻联系项相互迁移的过程中。综上可将 和 修正为（）（）（）式（）中：（）反映在 经再向 转化的过程中 协同前一次转化的信息向 转化的可能性，（）反映在 经再向 转化的过程中 协同前一次转化的信息向 转化的可能性，它们均反映非相邻联系项迁移演化过程的半偏联系数。此外，式（）反映了 偏向 和 偏向 的情况，而作为不确定项，偏同项 不仅会偏向，而且可偏向；同理，偏反项 不仅会偏向，而且可偏向。因此，结合联系数同异反关系相互联系转换的不确定性特点和演化半偏联系数的思想，可完善式（），得到连续演化半偏减法集对势：（）（）（）（）考虑到“”“”在式（）中作为分母，且、为非负数，故当二者为 时，有 或，即在进行第二次半偏过程中，不用考虑偏同项或偏反项的转化迁移信息，因此式（）可改写为人 民 黄 河 年第 期（）（）（）（且）（）（但）（）（但）（且）（）由式（）可知：当同一度项，不确定项和对立度项为 时，集对势取最大值（）；当对立度，其他项为 时，集对势取最小值（），故式（）连续演化半偏减法集对势（）的取值范围为，。为对应本文评价标准等级数目，可定义（）分为 个集对势的势级：反势（），），偏反势（），均势（）（，偏同势（）（，同势（）（，。当指标值联系数的集对势处于偏反势或反势状态时（），可将该指标诊断为旱灾风险评估等级高的主要原因，是加强区域抗旱减灾能力的重要指标和政府开展旱灾防控工作的着力点。为进一步检验连续演化半偏减法集对势是否合理，引入平均绝对误差公式。与五元减法集对势方法所得结果对比，式（）计算的连续演化半偏减法集对势的平均绝对误差为 （）（）（）式中：为随机模拟五元联系数的次数，（）、（）分别为第 次模拟的五元减法集对势、连续演化半偏减法集对势。使用程序随机模拟得到、和个联系数两种集对势（）和（）间的平均绝对误差 均为 左右，说明两种集对势方法计算结果比较接近，连续演化半偏减法集对势可用于反映联系数确定性整体发展趋势，同时也验证了式（）中差异度系数、经验取值具有合理性，而（）的计算过程中没有人为设定差异度系数，考虑了半偏过程中分量的演化特性，可见（）计算结果具有客观性和可解释性。实例分析济南市地处山东省中部，属暖温带季风气候，客水水源主要为黄河、长江，拥有丰富的地下水资源和独特的泉水景观，但人均水资源量仅为 人，是典型的缺水型城市，同时经济社会的快速发展导致水资源需求增加，进一步加剧了济南市供需不平衡的情况。综合考虑济南市地理环境、经济社会等因素对旱灾风险等级的影响，参考文献，济南市旱灾风险评估指标体系包括敏感性、适应性和危险性 个子系统，各子系统指标值选取、等级标准及权重见表。表 济南市 年旱灾风险评估指标系统、等级标准及权重子系统指标权重风险等级标准轻度（级）较轻（级）中度（级）较严重（级）重度（级）敏感性人均水资源量 （人）地下水水位 森林覆盖率 第三产业产值比率 城市日用水量 （）适应性水库调蓄能力 （）城市水资源协调能力 万元 用水量 （万元）污水达标排放率 工业用水循环利用率 危险性降雨量 致灾强度 出现频率 重现期 人 民 黄 河 年第 期 根据表 中济南市 年旱灾风险评估指标体系中的子系统划分和指标值，将指标值代入式（）式（），得到指标值与各等级的联系数，进一步计算相对隶属度并归一化可得各联系数分量，代入式（）得评估指标值联系数。同理，将表 中各指标权重代入式（）得样本 的指标值联系数，用级别特征值法得到样本所属风险等级；再用五元减法集对势、连续半偏减法集对势、连续演化半偏减法集对势方法计算这些联系数的集对势并判断态势。样本联系数及 种方法集对势计算结果见表，多年趋势变化对比如图（）所示。表 济南市 年旱灾风险评估样本联系数及 种方法集对势值年份样本联系数五元减法集对势态势连续半偏减法集对势态势连续演化半偏减法集对势态势综合评价等级均势均势均势偏同势偏同势偏同势偏同势偏同势偏同势均势均势均势同势偏同势同势同势偏同势同势同势偏同势同势 然后，将文献，中的各子系统指标值和权重代入式（）式（），计算出济南市 年旱灾敏感性、适应性和危险性 个子系统逐年连续演化半偏减法集对势和五元减法集对势，由折线图对比分析两种方法所得结果，如图 （）图（）所示。由表、图 可知：（）种方法评价结果对比（）敏感性（）适应性（）危险性图 济南市 年旱灾风险评估结果对比 ）经连续演化半偏减法集对势方法计算所得结果与其他两种方法基本一致，敏感性、适应性和危险性 个子系统在连续半偏减法集对势、五元减法集对势、连续演化半偏减法集对势方法中趋势相同、数值接近。说明本文方法体现了五元减法集对势的计算原理，计算结果合理有效。）从表 中 种方法计算的集对势来看，三者计算结果差距较小，所得态势也基本一致，但总体上连续演化半偏减法集对势比连续半偏减法集对势偏高，原因是前者多考虑计算了偏反项偏同一和偏同项偏对立的迁移信息。由评估结果可知，济南市 年的旱灾风险综合评价等级为 级，偏同的趋势要多于偏反的趋势，计算时补充项加的部分多于减的部分，反映在本文所用的数据资料下改进后的集对势略大于原来的连续半偏减法集对势，且更接近于五元减法集对势。）由图 可知，集对势越偏向两端，五元减法集对势和连续演化半偏减法集对势结果相差越大，后者结果更靠近于，因此反映在折线图上连续演化半偏减法集对势变化趋势偏缓一些，主要原因是连续演化半偏减法集对势将不确定项划分得更详细，从微观层面运用联系数分量的潜在转化机制将不确定项分量向确定项转移，缓和了集对势的变化趋势，相比主观设定差异度系数，本文计算方法所得结果更显合理。）济南市 年旱灾风险总体态势从均势上逐渐转好，但在 年发生恶化，同时图（）和图（）集对势也有下降趋势，尤其是危险性子系统下降幅度较大，可判断出敏感性子系统和危险性子系统对整个系统的集对势影响偏大。从单个子系统来看，敏感性子系统比较稳定，从偏反势逐渐好转发展为均势。适应性子系统变化大但后期发展稳定，由偏反势发展为同势。危险性子系统波动较大，但一直处于偏同势和同势，最后稳定在同势。评估整个旱灾风险和子系统风险后，诊断识别影响旱灾风险评估的重要指标，进一步分析计算各指标值联系数的集对势变化趋势。将数据代入式（）式（）计算各指标值联系数，再计算整理各联系数的连续演化半偏减法集对势值，绘制集对势逐年变化折线图，分析济南市 年各子系统指标项发展趋势和波动状况（见图）。人 民 黄 河 年第 期图 济南市 年各子系统指标连续演化半偏减法集对势变化趋势 由图（）可知：在敏感性子系统中，人均水资源量稳定地处在偏反势，说明济南市未有效改善水资源紧缺问题，是增高敏感性子系统风险等级的重要指标。地下水水位变化趋势与敏感性子系统相似且总体上变好，说明地下水的变化对敏感性子系统集对势有较大影响，合理开采地下水可有效降低旱灾风险敏感性；另外，森林覆盖率从偏反势逐步向均势发展，说明济南市逐渐重视城市绿化发展建设。第三产业产值比率由均势慢慢变为偏同势，也在逐渐转好。城市日用水量一直处于偏反势状态，且逐渐向反势靠近，对敏感性子系统有明显的消极作用。由图（）可知：在适应性子系统中，个指标整体呈现逐渐偏好的发展趋势。一方面，水库调蓄能力从反势逐渐转好，最终稳定在同势，同时其发展趋势与适应性子系统集对势折线图相似，说明该指标是影响济南市适应性水平的重要指标。城市水资源协调能力从偏反势逐步发展为同势，说明城市水资源协调能力对降低旱灾风险有重要作用。另一方面，万元 用水量一直保持同势状态，反映了济南市在水资源利用方面的高效和可持续性，在经济发展的同时能够有效地管理和利用有限的水资源。污水排放达标率和工业用水循环利用率从最初的均势或偏同势向同势变化，表明济南市工业部门在污水处理和节约用水方面的管理能力逐年提高，对降低适应性风险等级有积极作用。由图（）可知：在危险性子系统中，个指标集对势在 年以不同的趋势下降到均势范围，在 年达到最差状态，而后上升并稳定在同势状态，说明济南市经历了一段干旱强度变高的时期。其中降雨量波动幅度最小，总体上态势逐渐向同势发展，致灾强度和出现频率变化趋势相似，每年的态势相同，重现期一直处于偏同势或同势状态，说明济南市 年期间未遇重大旱灾状况。结论）本文结合半偏联系数和联系数分量迁移转化的思想，遵循微观运动跟随宏观变化的规律，充分考虑不确定项转向确定项和不确定项之间互相转化的潜在可能。在连续半偏减法集对势的基础上，对偏同势 和偏反势 进行分量补充，除确定项外，考虑将 分配到对立度、分配到同一度 的分量补充到连续半偏减法集对势中，并在半偏联系数的基础上，针对非相邻联系分量的转化情况，修正第二次半偏公式，提出一种演化半偏联系数，细致刻画半偏过程中分量的变化，使结果符合客观实际。经误差计算可知，本文改进方法与五元减法集对势计算结果十分相近，且对研究区域进行旱灾风险评估得到的风险系统发展趋势也是相同的，验证了本文方法的合理性。此外，本文方法计算得到的集对势与五元减法集对势相比，范围更小，发展趋势较为缓和，根据联系数分量之间的关系隐式确定差异度系数，比经验主观确定差异度系数更显合理、便捷。）由上述评价等级结果可知：济南市 年旱灾风险等级维持在较轻和中度附近（从 年的 级逐渐降低到 年的 级），济南市政府逐渐重视城市建设，从各方面采取措施降低旱灾风险等级，优化环境，取得了一定效果。通过样本与等级的联系数和 种方法计算得到的集对势值可看出，只有 年集对势是负值、但态势仍属于均势，年转变为偏同势，年略有恶化，从 年变回均势后又开始好转，态势一直稳定为同势。）由济南市旱灾风险敏感性、适应性和危险性 个子系统的发展趋势可知：在旱灾风险评估中，敏感性子系统变化较微弱，由偏反势逐渐转为均势状态；适应性子系统从偏反势状态转变为同势状态，并且在 年发生了跃迁，此后维持在同势状态；危险性子系统集对势波动较大，但状态一直在均势和同势之间变化，说明济南市采取的抗旱减灾措施取得了一定效果。经连续演化半偏减法集对势法诊断可知，影响济南市旱灾风险的主要指标有敏感性子系统中的人均水资源、森林覆盖率和城市日用水量，适应性子系人 民 黄 河 年第 期统中的水库调蓄能力及城市水资源协调能力。参考文献：王利娜，苏静，郑晓东，等我国农业干旱风险研究进展简述水科学与工程技术，（）：赵丽，冯宝平，张书花国内外干旱及干旱指标研究进展江苏农业科学，（）：金菊良，宋占智，崔毅，等旱灾风险评估与调控关键技术研究进展水利学报，（）：张继权，刘兴明，严登华综合灾害风险管理导论北京：北京大学出版社，：金菊良，马强，崔毅，等基于三元链式传递结构的区域旱灾实际风险综合防范机制分析灾害学，（）：章国材气象灾害风险评估与区划方法北京：气象出版社，：倪长健论自然灾害风险评估的途径灾害学，（）：金菊良，郦建强，周玉良，等旱灾风险评估的初步理论框架灾害学，（）：赵克勤集对分析及其初步应用大自然探索，（）：王雅燕，周玉良，周平，等基于集对分析的安徽省农业旱灾风险评估水电能源科学，（）：黄崇福，刘新立，周国贤，等以历史灾情资料为依据的农业自然灾害风险评估方法自然灾害学报，（）：刘亚彬，刘黎明，许迪，等基于信息扩散理论的中国粮食主产区水旱灾害风险评估农业工程学报，（）：，（）：，：，（）：徐昆，朱秀芳，刘莹，等采用 作物生长模型研究中国玉米干旱脆弱性农业工程学报，（）：金菊良，孔令茹，崔毅，等基于五元半偏减法集对势的区域农业旱灾脆弱性评价农业机械学报，（）：金菊良，沈时兴，郦建强，等基于联系数的区域水资源承载力评价与诊断分析方法华北水利水电大学学报（自然科学版），（）：董涛，王振龙，金菊良，等基于风险矩阵和五元减法集对势的区域旱灾风险链式传递诊断评估方法灾害学，（）：汪洁，金菊良，周戎星，等基于连续半偏减法集对势的区域旱灾风险评估方法灾害学，（）：陈守煜水资源与防洪系统可变模糊集理论与方法大连：大连理工大学出版社，：金菊良，吴开亚，魏一鸣基于联系数的流域水安全评价模型水利学报，（）：赵克勤集对分析及其初步应用杭州：浙江科技出版社，：金菊良，陈鹏飞，陈梦璐，等水文水资源学家陈守煜先生学术研究的知识图谱分析水利学报，（）：金菊良，陈鹏飞，张浩宇，等五元减法集对势及其在水资源承载力趋势分析中的应用华北水利水电大学学报（自然科学版），（）：，金菊良，赵行玥，崔毅，等半偏减法集对势在区域旱灾风险评估中的应用水利水运工程学报，（）：李健，夏文菊，赵剑锋济南市水资源开发利用存在问题与建议山东水利，（）：，桑国庆区域干旱风险管理研究济南：山东大学，：【责任编辑 简 群】（上接第 页），（）：，：马雪燕植被恢复对黄土高塬沟壑区产流模式的驱动效应杨凌：西北农林科技大学，：刘卓昕植被恢复对黄土丘陵沟壑区土壤水文性质及产流模式的影响研究 杨凌：西北农林科技大学，：，：刘姗姗，李小兵，刘思君，等黄土丘陵沟壑区典型小流域不同空间尺度水沙变化研究泥沙研究，（）：李蔚，冉启华沟谷横断面类型对产流过程的影响水利学报，（增刊）：芮孝芳，宫兴龙，张超，等流域产流分析及计算水力发电学报，（）：【责任编辑 张 帅】人 民 黄 河 年第 期