一种多智能体配置方法与流程

本发明涉及水资源配置领域，更具体地，涉及一种多智能体配置方法。

背景技术：

水资源优化配置是人类可持续开发、利用水资源的有效调控措施，是解决水资源供需矛盾、增强水资源系统抗风险性的根本途径。由于水资源配置涉及多个决策层次，地区与部门多个决策主体，近期与远期多个决策时段，社会、经济、生态环境多个决策目标，以及水文、生态、工程、市场等多领域，因此是一个复杂性的决策系统。近年来，尽管水资源配置系统的研究成果卓越，但在将系统中各种各样的个体行为与系统的整体演化有机的联系方面还难以得到有效的描述。

自20世纪70年代末，出现了多智能体系统(mas:multi-agentsystem)的理论与方法，并且由于mas能体现人类的社会智能，具有更大的灵活性和适应性，因此已成为一种复杂系统分析和模拟的思想方法和工具。国外已运用mas技术在需水预测、水资源管理等方面，但还只是处于初始阶段，真正运用于水资源配置系统分析和优化配置建模上还非常有限。如becu等建立了基于agent的catchscape模型，用于考察流域内的水力学、农户行为和水管理措施，并在泰国北部得到应用，feuillette等针对突尼斯的凯万市地下水过量开采问题，在水资源需求与地下水含水层的来水之间建立了多智能体的统一管理的协商模型；lebars等基于主体仿真模型研究了集体农业用水户行为对水资源配置的影响，lucillagiannetti等建立了基于智能体的城市供水系统的优化模型。然而，目前尚缺乏一个比较完整的多智能体模型来分析水资源配置系统中各层次的相互关系，描述其中各主体行为与整体演化的有机联系，从而模拟水资源分配过程，进行水资源的优化配置。

技术实现要素：

本发明为克服上述现有技术所述的至少一种缺陷，提供一种多智能体配置方法，将智能体系统的理论与方法应用于水资源配置领域，确立了一套比较完整的多智能体模型，有助于分析水资源配置系统中各层次的相互关系，描述其中各主体行为与整体演化的有机联系，进而模拟水资源分配过程，进行水资源的优化配置。

为解决上述技术问题，本发明采用的技术方案是：一种多智能体配置方法，包括以下步骤：

s1.根据研究区人类用水的特点，确立水资源配置系统的层次结构；

s2.建立基于多智能体系统mas的水资源优化配置模型框架：首先确立水资源配置系统中智能体的组成与系统结构，并在智能体组成与层次结构的基础上构建水资源优化配置模型框架结构，其次确定智能体之间的协作关系、协商机制；

s3.分析水资源优化配置模型中各智能体agent的模型行为，确定各智能体的约束和边界条件；

s4.利用优化算法对模型进行求解。

进一步地，所述的s1步骤中系统的层次结构为由流域级、部门或地区级、用户级所组成的结构，分别对应于mas的系统层、群体层和agent层，能够反应部门或用户用水与流域水量之间的循环反馈关系，进而研究局部的用水变化对流域全局水资源配置的影响。

进一步地，所述的s2步骤中多智能体的组成包括：水源agent、需水agent、供水agent、水资源调配agent；所述的水源agent用于分析水文气象变化特点，给出不同时空水资源量的大小、水质情况，与供水agent协商，实现任务委派，监测水质变化，执行调配计划；所述的需水agent用于分析社会经济发展状况，预测水资源需求量，将水资源需求的统计信息传递给水资源配置agent；所述的供水agent用于考虑不同区域工程特点和规模，形成供水任务分配联盟，与水源agent协商，形成任务委派，向水资源调配agent传递任务执行情况，与水资源调配agent协商，完成供水agent系统优化，确定水资源配置目标；所述的水资源调配agent用于接受需水agent、水源agent、供水agent的约束，指导整个水资源系统的优化配置，将分配结果传递给供水agent和用水agent，指导水资源的应急调度。

进一步地，所述的流域级与区域级之间包括对话层、问题求解层和控制层，其中，对话层由通信模块构成，主要完成与其他agent或外在环境的信息交互；问题求解层由学习机、推理机、规则库和知识库四个部分组成，主要完成和其他agent进行协调并生成最终决策的功能；控制层由控制模块构成，主要完成指导控制任务，并将控制任务的信息通过通信层传递给其他agent。

进一步地，所述的s2步骤中，智能体之间的协作关系主要基于每个智能体之间本身的水利联系，以及蓄水工程的影响，协商机制采用基于多agent的合同网协议。合同网是agent间建立合作机制的主要方法，其主要原理是采用市场机制进行任务通告、投标，最后签订合同来实现任务分配(资源分配)。合同网系统由多个节点组成，两个节点就任务的委托和承揽构成合同关系，一组这样的节点构成了合同网。在该系统中，每个节点代表一个agent，agent间通过招标一投标一中标过程进行任务分配和解决资源、知识冲突。在基于mas水资源配置系统中，控制agent是招标方，任务agent则作为投标方，控制agent将待分配的水源多少、水质状态、供水工程大小、各区域agent需水量多少在系统内公布，各个任务agent根据其内部资源(即需水量)采取相应的竞标策略，主动提出各自的投标方案，经过相互磋商，控制agent将各个时段的不同水质的水资源分配给合适的任务agent(即分配到不同的区域和不同的行业)，使得agent之间以一个较高的概率达成一个足够好的协定。

在本发明中，智能体的组成分为：水源agent(包括水库调蓄agent)、需水agent、供水agent、水资源调配agent四个部分。水资源优化配置模型框架结构在流域级与区域级之间包括对话层、问题求解层和控制层。而在每个区域agent里又包括主要由供水agent、水源agent、需水agent、水资源调配agent等构成的群体行为层。智能体之间的协作关系主要基于每个智能体之间本身的水力联系，以及蓄水工程的影响。而协商机制采用基于多agent的合同网协议。

进一步地，所述的s3步骤中，所述的智能体的约束和边界条件，对于水源agent主要考虑区域水资源量约束、河渠节点水量平衡约束、水库枢纽水量平衡约束、地下水库水量平衡约束及水质要求；对于需水agent主要考虑农业需水要求、工业需水要求、生活需水要求、最小生态需水要求、各行业用水不能超过需水量约束、水质要求；对于供水agent主要考虑水库蓄水库容、引提水工程能力、调水工程能力的约束；对于水资源调配agent主要考虑使流域整体经济效益最大、造成的污染最小、区域总缺水量最小三个目标，以及区域供需水量平衡约束、水费约束及非负约束。

进一步地，所述的优化算法为人工免疫算法。本专利在人工免疫优化求解中各个抗体、抗原(可行解)均采用实数编码，即每个分量用0.ai1ai2…ail，而其中的每个基因ai的取值为0到9之间的整数，而各个抗体、抗原与各区各行业的需水之积为各个智能体的决策变量。抗体群、抗原群初始规模分别为100、20；记忆细胞的规模设为20；调节因子β为2；终止迭代次数设为200。

与现有技术相比，有益效果是：

1.本发明将智能体系统的理论与方法应用于水资源配置领域，确立了一套比较完整的多智能体模型，有助于分析水资源配置系统中各层次的相互关系，描述其中各主体行为与整体演化的有机联系，进而模拟水资源分配过程，进行水资源的优化配置。

2.本发明提出的基于多智能体系统的水资源优化配置模型在分析水资源配置系统层次的基础上，构建了水资源优化配置模型框架，描述了同一层次智能体之间的行为和相互关系以及不同层次间的信息传递和作用机制，既考虑到水资源配置系统各智能体微观的决策行为，又具有强大的整体自组织能力，能够更加科学地反映系统的结构和内在运行机制，可以被广泛应用于水资源配置领域。

附图说明

图1是本发明方法流程图。

图2是本发明基于mas的水资源优化配置模型中流域与区域级层次间的系统结构。

图3是本发明步骤二中区域agent的系统结构图。

图4是本发明步骤四中人工免疫优化算法的计算流程图。

具体实施方式

附图仅用于示例性说明，不能理解为对本发明的限制；为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。附图中描述位置关系仅用于示例性说明，不能理解为对本发明的限制。

实施例1：

如图1所示，本发明提供的一种多智能体配置方法，主要包括以下步骤：

步骤一.确立水资源配置系统的层次结构。具体来说就是分析根据人类用水过程的特点，将水资源配置系统分为流域级、部门(或地区)级、用户级三个层次，分别对应于mas的系统层、群体层和agent层。其中，agent为各智能体行为层的各个节点，即各个用水户；群体代表智能体的聚合体，即不同地区或不同行业(或产业)部门；在智能体和群体行为基础上形成流域水资源配置的系统演变行为，即系统层。

步骤二.建立基于多智能体系统(mas)的水资源优化配置模型框架。

2.1确立水资源配置系统中智能体的组成与系统结构

根据水资源配置系统的组成可以将配置模型中的智能体主要简化为水源agent(包括水库调蓄agent)、需水agent、供水agent、水资源调配agent等几个部分，见表1所示。

表1水资源配置系统中智能体的组成与功能

以上每个agent又可主要分为三个层次：通信层、协作层和控制层。通信层由通信模块构成，主要完成与其他agent或外在环境的信息交互；协作层由学习机、推理机、规则库和知识库四个部分组成，主要完成和其他agent进行协调并生成最终决策的功能；控制层由控制模块构成，主要完成指导控制任务，并将控制任务的信息通过通信层传递给其他agent。此外，所有的agent的共同目标是使整个流域的水资源配置综合效益达到最大；而每一个行业部门、每一个用户agent都有自己的局部目标－各自部门效益；每个区域agent的局部目标就是自己区域内的综合效益达到最大。不用的用水户agent之间、用水行业agent之间、区域agent之间是相互影响、相互作用的

2.2构建基于mas的水资源优化配置模型框架结构

流域与区域级之间的求解系统结构可分为对话层、控制层和问题求解层，如图2所示：对话层主要接收需要优化配置的不同时段各区域(计算单元)的正效益和负效益，并将每个区域对水资源的需求、水资源可利用量等信息传递给控制层；控制层可以看作是水资源量、质的分配者，主要根据任务agent与其协商的结果给不同区域agent分配水资源；同时控制agent对各区域的配水情况进行监督和统计；问题求解层是由多个任务agent组成，每个任务agent对应于一个区域的需水问题，该层主要是针对各个区域的水资源优化配置总体目标，由多个任务agent协商合作进行问题求解。而每个区域agent主要由供水agent、水源agent、需水agent、水资源调配agent等，其系统结构图如图3所示。因此，图2和图3构成了基于mas的水资源优化配置模型的基本框架。

2.3确定智能体之间的协作关系、协商机制

同一层次内各个智能体之间、不同层次的智能体之间的关系通过物质流、能量流和信息流等多种方式进行交互作用产生的，在水资源配置系统中物质流主要包括水量、污染物等；能量流主要是指资金流；信息流主要是指水资源配置信息和政策法规信息，其中物质流中的水量和污染物是各个agent在配置中进行协作和协商的基础。在提取各个agent之间的水力关系时，要以系统中所存在的所有可能水量传递、转化和影响关系为主线，分类提取系统中各类agent并建立联系，如河道、湖泊(湿地、沼泽等)、水汇(海洋等)、地下水含水层、地下水侧渗传输渠道；还包括各类人工修建的水利工程，如蓄引提水工程、调水工程、海水直接利用(淡化)工程、集雨工程等；此外还包括各类用水户等。此外，蓄水工程对水资源时空分配有着重要调节作用，蓄水工程的运行过程就是水资源优化配置中控制设备的动作策略，也为各个agent在优化配置中进行协商奠定了基础。

水资源配置系统中的协商就是指供水agent、水源agent、需水agent、水资源配置agent等之间通过交换相关的各自结构化信息，形成一致的观点和规划，达成意图上的一致，使水资源配置目标在时空上达到最优。本发明采用smith和davis提出的合同网协议来建立各智能体之间的协商机制。

步骤三、分析水资源优化配置模型中各智能体(agent)的模型行为。

具体来说就是利用智能体的自利和永不满足的特性，实现局部的动态优化，即在一定约束条件和边界条件下，智能体具有自己的特定的目标，根据外部(包括环境和其他智能体)输入的不同—边界条件的变化，选择不同的参数，达到自身的目标和最佳状态。各agent的约束条件分别确立如下：

3.1水源agent必须满足以下约束条件：

(1)区域水资源量约束：

其中，awi(m,z)表示第m时段z区农业、工业、生活、生态分配的水资源量；wq1、wq2、wq3、wq4分别表示本地水资源量、调入水资源量、上游区域的回归到z区的水资源量、调出水资源量；wq表示水资源可利用量；λ表示回归系数。

(2)河渠节点水量平衡约束：

wh(m,z)表示第m时段第z节点的过水量；wqh表示区间来水量；wqrx表示第i个水库的下泄水量；wqrc表示给该区域供水的第i个水库的蓄水变化量；ql表示蒸发渗漏损失量。

(3)水库枢纽水量平衡约束：

vr(m+1,i)＝vr(m,i)+wqrc(m,i)-wqrx(m,i)-wqvl(m,i)

vr(m+1,i)表示第m+1时段第i个水库枢纽的库容；wqvl表示水库的水量损失。

(4)地下水库水量平衡约束：

wg(m,z)表示第m时段第z区地下水开采量；wgmax(z)表示z区地下水年最大开采量。

(5)水质变化模拟行为表现：

c(m,z)表示第m时段第z区的水质浓度；ww(m,z-1)表示第m时段第z-1区的污染物排放量；k表示染物的衰减系数；x表示第z-1区与第z区之间的概化距离；u表示流速；ξ表示污染物排放系数。

3.2需水agent必须满足以下约束条件：

(6)农业的需水要求：

aw1(m,z)≤awu1(z)wq(m,z)

aw1(m,z)≥awl1(z)wq(m,z)

awu1、awl1分布表示农业分配水上限、下限。

(7)工业的需水要求：

aw2(m,z)≤awu2(z)wq(m,z)

aw2(m,z)≥awl2(z)wq(m,z)

awu2、awl2分布表示工业分配水上限、下限。

(8)生活的需水要求：

aw3(m,z)≥awl3(z)wq(m,z)

awl3分布表示生活配水下限。

(9)最小生态需水要求：

aw4(m,z)≥daw4(m,z)min

daw4(m,z)表示第m时段第z区生态需水量的下限。

(10)各行业用水不能超过需水量：

awwi(m,z)≤dawi(m,z)

awwi(m,z)表示第m时段第z区第i行业的实际净供水量；dawi(m,z)第m时段第z区第i个行业需水量。

(11)水质要求：

其中n表示水质级别，1，2，3，4，5对应水质ⅰ，ⅱ，ⅲ，ⅳ，ⅴ类，6代表劣ⅴ类。

3.3供水agent模型应满足以下约束条件：

(12)水库蓄水库容：

vrmin(i)≤vr(m,i)≤vrmax(i)

vrmin(i)≤vr(m,i)≤vr′max(i)

vrmin(i)表示第i个水库的死库容、兴利库容；vr′min表示第i个水库的汛限库容。

(13)引提水工程能力：

wqp(m,z)≤wqpmax(z)

wqp(m,z)、wqpmax(m,z)分别表示第m时段第z区的引提水量、最大引提水能力。

(14)调水工程能力：

wd(m,z)≤wdmax(z)

另外，在保障最小河道内生态需水的前提下，不同行业供水的供水优先顺序依次为：生活、工业、农业和河道外生态。

3.4水资源调配agent模型：

(15)经济目标行为表现。要求流域整体的经济效益达到最大，用f1表示：

ω(z)表示第z区的权重，意为在流域中重要程度；b(z)表示水资源综合利用效益；c(z)表示水资源利用损耗及费用；bi(m,z)表示第m时段第z区第i行业的净效益。

(16)环境目标行为表现。选用水污染造成水环境质量损失最小作为水资源优化配置的环境目标，即可用重要污染物的排放量最小体现，用f2表示：

di(z)表示第z区第i行业的单位废水排放量中重要污染因子的含量(mg/l)；pj(z)表示第z区第i行业的污水排放系数。

(17)社会目标行为表现。不同区域的缺水量大小或缺水程度影响到社会的发展和稳定，故采用区域总缺水量来间接反映社会效益，用f3表示：

lawwi(m,z)表示第m时段第z区第i个行业缺水量；表示第z区第i行业的重要程度，可根据供水优先顺序来确定。

调配约束条件有：

(18)区域供需水量平衡：

lawwi(m,z)＝dawi(m,z)-awwi(m,z)

(19)水费约束：

aw3(m,z)ewpurb(z)-tp(z)aw3(m,z)(ewpurb(z)-bwpurb(z))≤r(z)income(z)

ewpurb(z)表示第z区城市超额水价；tp(z)表示城市基本生活用水量占生活用水中的比例；bwpurb(z)表示城市基础水价；r(z)表示第z区的生活用水水费支出与人均可支配收入的最高比例；income(z)表示第z区的人均可支配支出。

(20)非负约束，即所有的变量均为非负数。

步骤四、利用优化算法对模型进行求解。本发明提出的基于多智能体系统的水资源优化配置模型求解方法采用人工免疫优化算法，其计算流程图见图4。

本发明多智能体配置方法的具体实施例如下：

根据本发明所提供的步骤，将基于多智能体系统的水资源优化配置模型应用于泛东江流域。首先建立基于多智能体系统(mas)的水资源优化配置模型框架。模型的节点概化首先通过确定研究区域范围确定水资源系统模拟中需要考虑的实体类别，选择的依据为是否与水量水质转换有直接或间接联系，一般来说，包括河道、湖泊、地下水含水层等天然水源承载体；基本的计算分区和城市等构成的用水户；河流分水工程、蓄引提调水工程、供用耗排渠道等人工修建的水利工程；再次根据选出的实体及其类别用简练的参数抽象出各实体对水资源运动过程其作用的性能。最后，在对概化的各个实体之间水力联系进行描述的基础上，做出反映实际水资源系统及其内部关系整体框架的节点图。对泛东江流域模型的节点概化结果如图4所示。

其次，要确定研究区各智能体的约束和边界条件。水资源的需求量是以2005年为基准年，应用2010年泛东江流域各个计算单元各用水部门的需水预测成果，并结合各单元1956～2005年各年的来水频率和不同频率下水资源需求量月分配系数得到2010水平年需水长系列资料。各水利工程数据以现状工程数据为基础，模型以月为配水时段和年为计算周期进行1956～2005年间的长系列调算。求解过程中，大中型水库在流域内根据水力联系进行调算，各计算单元的小型水库则打捆计算。水质变化模拟中，主要污染物选取cod，应用公式(5)时采用月平均流量，并根据各主要断面的水位～面积曲线、水位～流量曲线组合得到流量～流速曲线，cod综合衰减系数取0.15/d，入境水、区间来水、水库集雨面积上的产水水质设为ⅰ类(或ⅱ类)水标准。

求解模型时，在人工免疫优化求解中各个抗体、抗原(可行解)均采用实数编码，即每个分量用0.ai1ai2…ail，而其中的每个基因ai的取值于0到9之间的整数，而各个抗体、抗原与各区各行业的需水之积则为第五章中各个智能体的决策变量。抗体群、抗原群初始规模分别为100、20；记忆细胞的规模设为20；调节因子β为2；终止迭代次数设为200。

综合以上步骤，以2005年的水利工程现状为基础，针对2010年社会经济发展规划，进行水资源的优化配置，其配置的主要指标结果如表2所示。

表22010年泛东江流域水资源优化配置结果(单位：亿m³)

综上，本发明多智能体配置方法是一种较为完整的基于多智能体系统的水资源优化配置方法，既考虑到水资源配置系统各智能体微观的决策行为，又具有强大的整体自组织能力,因此在建立水资源优化配置模型时有着巨大的优势。

显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。