南四湖流域水资源分配与优化调控：构建高效权重评估模型

南四湖流域水资源优化布局模型构建

2.3.1.1流域水资源体系简化

鉴于水资源体系的繁杂性，水资源布局需依据其目标与需求，对水资源体系的特性和演变规律进行恰当的简化。简化即指将现实的水资源体系转化为计算机能够识别的网络体系。简化的原则是需能较为真实地展现水资源计算分区的水量传递，并有利于揭示供需矛盾；同时，还需适应基础资料与数据源的精确程度，并便于分析计算。

根据南四湖流域实际情况及水资源计算分区，将南四湖流域水资源体系简化为由节点和有向线段组合的网络，形成了水资源体系网络图（图2.4）。图中节点包括关键计算单元、河渠等，各种水源的供水均基于各计算单元之上进行，有向线段代表天然河道或人工输水渠道，它们反映了节点之间的水流传输关系。

图2.4南四湖流域水资源体系网络示意图

2.3.1.2水资源优化布局模型描述

（1）布局方法

水资源优化布局需运用水资源体系分析方法来解决问题。在水资源体系分析中，数学模型扮演着至关重要的角色，水资源优化布局问题可通过构建数学模型来得到解决。水资源体系的数学模型一般包括目标函数和约束条件两部分。对于不同的体系和不同的水资源问题，数学模型是有所区别的。数学模型通常是根据体系的实际需求来设计目标函数，使目标函数值达到最大或最小，即体系达到最佳状态时获得的水资源优化布局方案。

水资源优化布局的目标是为了支撑全流域社会、经济、环境的全面协调和持续发展。水资源利用是多目标的，水资源优化布局就是多目标优化问题，其目标不是追求某一方面或对象的效益最佳，而应追求整体效益最佳。因此，水资源优化布局问题实际上是一种复杂的多目标决策问题。

根据南四湖流域自然地理与人文地理特点，采用多目标规划模型对流域的水资源进行合理布局。水资源优化布局多目标问题一般表达式如下：

变环境条件下的水资源保护与可持续利用研究

式中：x为决策变量；fp(x)为p个独立的目标函数向量；gi(x)为约束条件组；bi为右端常数项向量。

（2）目标函数

目标函数体现模型系统的目标要求。针对研究的问题不同，要求目标函数实现最大化或最小化。根据水资源优化布局的科学内涵，水资源优化布局是通过科学合理分配有限的水资源，以解决水资源的短缺和用水竞争问题，更好地满足生活、工农业生产及生态环境的需求。

对于南四湖流域来说，水资源布局的主要目的是着重研究如何联合运用多种水源（包括当地水源和外调水源），以缓解水资源短缺而引起的争水问题，结合本流域“十一五”规划要求，最终设定了以下2个目标函数。

1）经济效益目标。以区域供水带来的直接经济效益最大来表示。函数表达形式为

变环境条件下的水资源保护与可持续利用研究

式中：

、

为独立水源i、公共水源c向k子区j用户的供水量，万m3；

、

为独立水源i、公共水源c向k子区j用户供水的效益系数，元/m3；

、

为独立水源i、公共水源c向k子区j用户供水的费用系数，元/m3；

、

为k子区独立水源i、公共水源c的供水次序系数；

为第k子区的第j用户的用水公平系数；βk为第k子区的权重系数。

2）社会效益目标。由于社会效益难以衡量，而区域缺水量的大小或缺水程度会直接影响到社会的发展和稳定，故采用区域供水系统总缺水量最小来间接反映社会效益。函数表达形式为

变环境条件下的水资源保护与可持续利用研究

式中：

为k子区j用户的需水量，万m3，

、

分别为独立水源i、公共水源c向k子区j用户的供水量，万m3。

（3）约束条件

约束条件体现目标函数的限制条件。推求目标函数达到最优时的决策变量，应是在约束条件下求得的。在水资源优化布局中，产水量、供水量、输水建筑物的过水能力等都可能成为约束条件。

1）供水量约束。根据资源节约和有效利用的原则，不同水源供给计算分区各用水户的总水量不应多于其可供水量，如下式

变环境条件下的水资源保护与可持续利用研究

式中：

、

分别为规划水平年内独立水源i、公共水源c对k子区第j用户的供水量；

、Wc分别为k子区独立水源i及公共水源c的可供水量。

2）供水能力约束。各分区的输水河道及泵站都有各自的输水能力。因此，在水资源布局计算时，供水水源对计算分区各用水户的供水量不应大于其最大输水能力，如下式:

变环境条件下的水资源保护与可持续利用研究

式中：Qmax，ik、Qmax，ck分别为规划水平年内独立水源i及公共水源c对第k水资源分区输水能力。

3）部门用水量约束。各水源提供给各分区各用水户的水量不低于该部门的最低用水量，如下式:

变环境条件下的水资源保护与可持续利用研究

式中：

、

分别为规划水平年内k子区j用户的最小需水量和最大需水量。

4）变量非负约束。各个分区的任何用水户的用水量都不为负，所能提供的水量能满足每个用水户的需要，如下式:

变环境条件下的水资源保护与可持续利用研究

式中：

，

分别表示独立水源i、公共水源c向k子区第j用户供水量。

2.3.1.3模型的建立

将上述目标函数及各种约束条件组合在一起，即构成南四湖流域水资源优化布局的总体模型

变环境条件下的水资源保护与可持续利用研究

变环境条件下的水资源保护与可持续利用研究

式中：各符号的意义同式(2.10)至式(2.17)。

该水资源优化分配模型充分考虑了流域水资源的可持续利用和经济社会的可持续发展，以满足流域供水经济效益最大及水资源系统总缺水量最小为目标，并在流域划分水资源计算分区中，考虑各用水户及各供水水源之间的相互协调作用，使得相应的水资源分配最优。

该水资源合理配置模型全面考虑了流域水资源的持续使用和经济社会的持续发展，旨在实现流域供水经济效益最大化及水资源系统总缺水量最小化。在流域划分水资源计算分区时，考虑到各用水户及各供水水源之间的相互协调作用，以确保水资源分配达到最佳状态。

2.3.1.4模型参数确定

南四湖流域供水水源包括4个，分别是本地地表水、地下水以及跨流域调水。跨流域调水包括引黄河水和南水北调工程引长江水。根据流域水资源计算分区划分（图2.3），计算分区共4个。其中，引长江水供应所有分区，作为全区公共水源；引黄河水分别供应济宁及湖区和湖西菏泽区，作为两分区公共水源；各分区的当地地表水和地下水为单分区水源。

根据流域内各计算分区的实际情况，将用水部门具体划分为生活用水、工业用水、农业用水和生态环境用水等4个用水部门。

在水资源优化配置模型中，以k个计算分区内i水源提供给j用水户的供水量作为决策变量。根据上述分析，南四湖流域计算分区数k=4；其中，各分区独立水源数i=2，公共水源数c=1或2，用水部门数j=4，因此南四湖流域水资源优化配置模型拥有56个决策变量，44个约束条件。

(1)用水部门公平系数

确定用水优先权是模型分析计算的前提。根据流域用水部门的性质和重要程度，按照“先生活，后生产”的原则，在同一计算分区中把用水部门划分为不同的级别。

据《山东省水利发展和改革“十一五”规划》(2006年)、《江苏省水利发展“十一五”规划》(〔2006〕147号文)和《山东省经济发展“十一五”规划》(2006年)，流域内经济发展以工业为主，兼顾农业发展，优先保障生活用水，然后保障工业用水，最后安排农业用水，使其供水保证率分别达到98%以上、75%～90%、50%～75%。目前南四湖流域水污染较为严重[5，27]，生态系统遭到破坏，在保障生活用水的同时，还要充分考虑生态环境用水。

由此，将南四湖流域各用水部门划分为如下四个等级：第一级为生活用水，第二级为生态环境用水，第三级为工业用水，第四级为农业用水。当发生缺水时，级别低的用水部门先缩减供水，以确保级别高的用水户正常供水。

表示为第k子区、第j用水部门相对其他用水部门优先得到水资源供给的重要程度，它与优先得到供水的次序有关，可将各用水部门用水的优先程度转化为[0，1]区间上的系数。根据用水部门的性质和重要程度，可以由决策者参照下式并根据流域水资源实际情况确定。

变环境条件下的水资源保护与持续利用研究

式中:

为第k水资源区第j用水部门的用水次序序号;

为第k水资源区用水次序序号最大值，对南四湖流域而言，

。

经计算，各用水户的用水公平系数

分别为生活用水0.4，生态用水0.3，工业用水0.2，农业用水0.1。

(2)供水次序系数

供水次序系数

可反映k子区i水源相对于其他水源供水的优先程度。对于当地水，根据各水源调节能力的不同，将当地水资源的供水次序划分成不同的优先级，具有较低调节能力的水源先供水，具有较高调节能力的水源后供水。供水次序为河道提引河水和小型塘坝供水、山区水库供水，地下水供水；对于外调水，由于其成本相对较高，因而通常作为当地水资源的补充水源。本节根据外调水的成本和实际情况，在当地水源供水的基础上，确定外调水的供水次序为先黄河水供水，再长江水供水。

根据以上原则，确定南四湖流域各水源的供水次序为：①地表水；②地下水；③引黄河水；④引长江水。供水次序系数

可参照

的计算公式确定。

对于济宁及湖区、湖西菏泽区两水资源计算分区，

，经计算，地表水、地下水、引黄河水和引长江水的供水次序系数

分别为0.4，0.3，0.2，0.1。

对于湖东枣庄区和湖西徐州区两水资源计算分区，

，经计算，地表水、地下水和引长江水的供水次序系数

分别为0.5，0.33，0.17。

(3)供水效益系数及费用系数

1)效益系数。水资源优化配置数学模型，涉及各类水源的供水经济效益，是分析水资源优化配置的主要依据条件。水资源优化配置过程中，在满足生活用水、生态环境用水及各类生产部门用水最小需水量的前提下，将水资源量尽可能分配到经济效益较大的用水部门中去，最大限度的发挥水资源的经济效益；同时，要使整个研究区的缺水量最小。

南四湖流域农业生产用水的经济效益近年来有所增长，据统计，随着节水灌溉方式的普及，农业生产用水的经济效益显著提高。目前，南四湖流域农业生产用水的效益系数在15～20元/m3之间，至规划水平年(2015年)，农业生产用水的效益将成倍增长。工业用水经济效益较大，根据经济发展水平，工业用水效益多在150～250元/m3之间。

在计算过程中，农业用水、工业用水的经济效益系数采用山东省和江苏省统计部门提供的数据进行计算；生活、环境的效益是间接而复杂的，不仅有经济方面的因素，还有社会效益存在，其效益系数较难确定。根据生活、生态环境用水优先满足的配置原则，在计算中赋以较大的权值，用以表示其效益系数。由此，得出南四湖流域规划水平年(2015年)各用水部门的用水效益系数，见表2.17。

表2.17南四湖流域2015年各用水部门效益系数单位:元/m3

2)费用系数[57，62，63]。不同水源供水给各用水部门费用系数，参考水费征收标准确定。对有资料水源工程，根据资料计算确定；缺乏资料时，参考邻近地区同类水源工程选取。

2) 费用比率[57，62，63]。不同水源向各用水单位分配的费用比率，依据水费征收规范进行设定。对于拥有资料的水源工程，依据资料进行计算确定；资料不足时，参照邻近区域同类水源工程进行选择。

a. 从水厂获取水源的用户以水价为费用比率。

b. 从自备井获取水源的用户以水资源费、污水处理费与抽水成本之和作为费用比率。

c. 从水利工程获取水源的用户以水资源费、污水费与输送成本之和作为费用比率。

d. 农业用户的费用比率参照水费征收规范确定。

根据上述原则，并配合南四湖流域南水北调东线工程调水费用分析，得出流域规划基准年（2015年）供水费用比率，详见表2.18。

表2.18南四湖流域2015年供水费用比率单位:元/m3

(4) 需水量范围

设k子区j用户的需水量上限和下限分别为

、

。

其确定方式如下：

1) 生活需水量范围。根据生活用水特性，其上限和下限均取为生活需水量，即

变环境条件下的水资源保护与可持续利用研究

式中:

、

分别为生活用水的上限和下限;

为规划基准年的生活需水量。

2) 环境需水量范围。考虑到人们对环境用水的重视，环境用水的上限和下限也均取为环境需水量，即

变环境条件下的水资源保护与可持续利用研究

式中:

、

分别为环境用水的上限和下限;

为规划基准年的环境需水量。

3) 工业需水量范围。考虑工业用水的特性，工业需水量的上限和下限按以下公式确定

变环境条件下的水资源保护与可持续利用研究

式中:

、

分别为工业用水上限和下限;

为工业需水量。

4) 农业需水量范围。农业灌溉需水量的上限和下限需根据有效灌溉面积、保证灌溉面积和综合灌溉定额来确定，即

变环境条件下的水资源保护与可持续利用研究

式中:

、

分别为k子区农业需水量的上限和下限;

、

分别为k子区有效灌溉面积和保证灌溉面积;Gk为k子区的综合灌溉定额。根据式(2.19)～式(2.22)，可计算出南四湖流域各计算分区不同用水部门需水量的上限和下限，详见表2.19。

表2.19南四湖流域各计算分区不同部门需水量范围单位:万m3

(5) 权重比率

目标权重比率λp表示p个目标对其他目标的重要性程度；子区权重比率βk表示k子区对整个区域的重要性程度。本节采用层次分析法确定权重比率βk和λp。

1) 层次分析法求解的基本思路。层次分析法是美国运筹学家T.L.Saaty在20世纪70年代提出的一种多目标决策分析方法，属于定性与定量分析相结合的方法，是一种将决策者的复杂系统决策思维过程模型化、定量化的过程。应用层次分析法，决策者可以将复杂问题分解为若干层次，每个层次包含若干因素；在各层次、因素间进行比较和计算，可以得到表示方案重要性程度的权重，为最优方案的选择提供依据。层次分析法适用于多目标、多层次的非结构化、半结构化决策问题，在系统评价、方案比较等方面得到了广泛的应用。在流域水资源规划方案、工程设计方案、工程施工方案等的比较与优选中，均可以考虑采用层次分析法[56，64]。

本节将应用层次分析法确定经济效益和社会效益两个子目标的权重比率λp及四个水资源计算分区在整个流域中的权重比率βk。确定因子权重的具体步骤如下[66～68]:

a. 建立层次结构模型，如图2.5所示。

图2.5南四湖流域多目标优化递阶层次结构图

b. 构建判断矩阵。对于建立的层次结构模型，需要逐层计算相关因素间的重要性，并予以量化，构成判断矩阵，作为进一步分析的基础。对各因素因子间两两进行比较，用bij表示针对上一层次的某因素而言，本层次与之有关因素之间的相对重要性，采用Saaty提出的9级标度法进行量化，见表2.20。

表2.20 Saaty标度法及其含义

c. 层次单排序及一致性检验。

(a) 计算各指标权重值。

第一步，计算判断矩阵中每行元素的几何平均值

变环境条件下的水资源保护与可持续利用研究

第二步，将

归一化，即

变环境条件下的水资源保护与可持续利用研究

可得到近似特征向量ω=[ω1，ω2，…，ωn]T

第三步，计算判断矩阵的最大特征值λmax

变环境条件下的水资源保护与可持续利用研究

式中:(AW)i为向量BW的第i个元素

(b) 判断矩阵偏差一致性检验。由于判断矩阵的构造是由决策人员的定性分析转入定量描述的，因而无法保证完全一致性，需进行检验，目的是使差异不致过大。

由判断矩阵的偏差一致性指标CI的表达式

CI=(λmax-n)/(n-1)(2.26)

引入判断矩阵的随机一致性比率CR=CI/RI，判断矩阵是否具有满意的一致性，其中RI为平均随机一致性指标，其值可从表2.21查得。若CR＜0.1，则认为判断矩阵具有满意的一致性；否则，需要对判断矩阵进行适当调整直到具有满意的一致性为止。

表2.21平均随机一致性指标RI

d. 层次总排序及一致性检验。从层次结构模型的第二层开始，逐层计算各层相对于最高层相对重要性的排序权值，称为层次总排序。

由上述步骤得到每一个要素相对于上一层次对应要素的权重值后，通过层次总排序计算出每一个评价指标相对于总目标整个研究区水资源开发利用评价的权重值。最后，计算各层次所有元素对总目标相对重要性的排序权值。

层次总排序后同样要进行一致性检验，假设第K层层次总排序权值为αi(i=1，2，…， n)，一致性指标为CIi，相应的平均随机一致性指标是RIi，则总排序的一致性指标

变环境条件下的水资源保护与可持续利用研究

总排序的平均随机一致性指标

变环境条件下的水资源保护与可持续利用研究

变环境条件下的水资原维护与持续开发研究

1) 层次总排序一致性评价。当CR=CI/RI<0.1时，视为层次总排序的一致性令人满意，否则，需重新调整判断矩阵，直至满意。

2) 权重求解计算。依据上述层次分析法的求解步骤，对南四湖流域优化布局中的各个权重值进行求解。

a. 目标权重系数λp确定。

(a) 如图2.5所示的南四湖流域水资原系统规划的层次结构模型，根据各因素的类别及支配关系，分为目标层、准则层、措施层。目标层为流域水资原规划的总体目标，即南四湖流域水资原的优化布局；准则层是为衡量总体目标能否实现的标准，模型确定了经济效益和社会效益两准则；措施层是根据流域具体情况及发展规划等所设置的若干个技术经济可行的布局方案，模型为经济效益和社会效益两准则分别设立了9个措施。

(b) 对经济效益子目标(B1)和社会效益子目标(B2)在影响水资原优化布局综合评价结果的关键程度方面进行两两比较，其结果见表2.22。

其中，λmax=2，CI=0，判断矩阵具有完全一致性。对于二阶矩阵而言，总是一致的，不必检验。

表2.22水资原优化布局综合评价下判断矩阵A-B

对措施层中的9个措施(C1～C4)、(C5～C9)在影响水资原优化布局综合评价结果的关键程度方面进行两两比较，其结果见表2.23、表2.24。

表2.23水资原优化布局综合评价下判断矩阵B1-C

其中，λmax=4.1755，CI=0.0585，CR=0.0650<0.1，具有满意一致性。

表2.24水资原优化布局综合评价下判断矩阵B2-C

其中，λmax=5.3522，CI=0.0881，CR=0.0786<0.1，具有满意一致性。

将所有判断矩阵进行一致性检验，由上述分析可知，判断矩阵A-B、判断矩阵B1-C和判断矩阵B2-C这三个判断矩阵均具有满意一致性。

(c) 在得到每一个要素相对于上一层次对应要素的权重后，通过层次总排序，计算出每一个措施相对于总目标南四湖流域水资原优化布局综合评价权重值。计算各层次所有元素对总目标相对重要性的排序权重，总排序的结果见表2.25。

表2.25水资原优化布局层次总排序分析结果表

续表

在层次总排序之后，利用CR=CI/RI进行整个层次的一致性检验，总排序随机一致性比率为，所以层次总排序具有满意的一致性。该流域应用层次分析法求解有效。

故在水资源优化布局总目标下，经济效益和社会效益两个子目标的排序权重系数λp分别为0.1667、0.8333。

在决策过程中，权重还可以与决策者交互调整，不同的权重值可得出原多目标规划问题的一个非劣解，为决策者提供更多的有关目标权衡比较的信息，以便选择最佳权衡解[56]。

b. 子区权重系数βk确定。针对南四湖流域四个水资原计算分区(济宁及湖区、湖东枣庄区、湖西菏泽区及湖西徐州区)，同样采用层次分析法确定其权重系数βk。针对该流域的具体情况，经过层次分析法分析计算后，分别拟定为β1=0.1299，β2=0.5567，β3=0.2556，β4=0.0577。

以2015年为规划水平年，将前述分析计算的各相关参数代入模型中，进行水资原优化布局多目标规划模型求解。

权重差值是什么意思

在统计学中，我们常常会使用权重来表示数据的重要性。权重越大的数据，对整个数据集的贡献也就越大。而权重差值则是比较两个数据点的权重大小差异。例如，在进行市场调研时，我们可能会为不同年龄段的人群赋予不同的权重，这样就可以更准确地反映出不同年龄段人群对产品的看法和需求。如果我们比较两个不同年龄段的权重，可以计算出它们之间的权重差值，从而判断哪个年龄段对该产品更为重要。

权重差值还可以用于分析机器学习模型的效果。在模型训练过程中，我们会将不同属性赋予不同的权重，以便模型能够更好地预测结果。而权重差值可以帮助我们判断哪些属性对模型预测结果的贡献更大。例如，我们训练了一个人脸识别模型，可以计算出每个像素点的权重，从而计算出不同区域的权重差值。如果我们发现某个区域的权重差值很大，说明该区域对人脸识别的贡献很大，我们就可以在优化模型时更加注重这个区域的属性。

权重差值还可以用于评估投资组合的表现。在投资时，我们会对不同的资产进行权重分配，以便达到最佳的风险收益比。例如，我们可能会将资产组合分配为60%的股票和40%的债券。如果股票市场表现不佳，可能会使得整个资产组合的表现下降，我们需要通过计算权重差值来评估股票市场对组合的影响。如果股票市场权重差值很大，说明股票对组合的收益贡献很大，我们就需要重新考虑资产组合的分配。