1.一种多模态多目标分层分布式综合能源系统经济调度方法，其特征在于，将解析式多模态多目标方法和分层分布式一致性方法进行结合，用于综合能源系统的经济调度，增加经济调度的多样性和稳定性，提高计算速度和精度，为决策者提供丰富的能进行替换的决策方案；在使用过程中的步骤为：步骤(1)：综合能源系统的经济调度模型的构建包含火能、水能、风能、光能和地热能，并以成本花费和碳排放量为多目标，遵循经济调度的等式约束和不等式约束；

总发电成本的目标函数为：

其中， 为总的发电成本； 表示在t时刻第i台发电机组的成本； 表示在t时刻第j台火电发电机组的成本； 表示在t时刻第k台水电发电机组的成本；

表示在t时刻第l台风力发电机组的成本； 表示在t时刻第r台光伏发电机组的成本； 表示在t时刻第s台地热能发电机组的成本； 表示在t时刻第u台清洁能源发电机组的成本；N为总发电机组的数量；NTG为火电发电机机组的数量；NHG为水电发电机机组的数量；NWG为风力发电机机组的数量；NPG为光伏发电机机组的数量；NOG为地热能发电机机组的数量；NCleanG为清洁能源发电机机组的数量；T为目标函数的统计时间；

为t时刻第i台发电机组的发电量； 为t时刻第j台火电发电机组的发电量； 为t时刻第k台水电发电机组的发电量； 为t时刻第l台风力发电机组的发电量； 为t时刻第r台光伏发电机组的发电量； 为t时刻第s台地热能发电机组的发电量； 为t时刻第u台清洁能源发电机组的发电量；aj、bj和cj为第j台火电发电机组成本系数；dk为第k台水电发电机组成本系数；el为第l台风力发电机组成本系数；gr为第r台光伏发电机组成本系数；

ms为第s台地热能发电机组成本系数；ku为第u台清洁能源发电机组成本系数；用清洁能源机组成本 替换水电成本 风电成本 光能成本 和地热能成本

总碳排放量的目标函数为：

其中， 为总的碳排放量； 表示在t时刻第j台火电发电机组的碳排放量；αj、βj和γj为第j台火电发电机组的碳排放系数；

功率平衡的等式约束为：

其中， 为t时刻总的负荷需求值；

功率的上下限约束为：

其中， 为第j台火电发电机组的发电量的下限； 为第j台火电发电机组的发电量的上限； 为第k台水电发电机组的发电量的下限； 为第k台水电发电机组的发电量的上限； 为第l台风力发电机组的发电量的下限； 为第l台风力发电机组的发电量的上限； 为第r台光伏发电机组的发电量的下限； 为第r台光伏发电机组的发电量的上限； 为第s台地热能发电机组的发电量的下限； 为第s台地热能发电机组的发电量的上限；

火电机组爬坡速率约束为：

其中， 为在t‑1时刻第j台火电发电机组的发电量； 为第j台火电机组的向下爬坡速率； 为第j台火电机组的向上爬坡速率；T60为60分钟；

光电机组的时间约束为：

步骤(2)：按比例R对权重因子η进行初始化，令

其中，R为对权重因子取值的比例；η为权重因子，其取值范围为0≤η≤1；nη表示对权重因子取值范围分段的段数；

步骤(3)：通过线性加权的方法将多目标问题转变为单目标问题；

其中，fTotal为线性加权后的总目标值； 为t时刻第j个火电机组线加权后的目标值；

步骤(4)：设置初始迭代时刻t＝1；

步骤(5)：输入预测负荷值，并令初始迭代步数k＝0；

步骤(6)：将综合能源系统分为四层：第一层分为区域A、区域B和区域C；第二层分为区域Aa、区域Ab、区域Ac、区域Ba、区域Bb、区域Bc、区域Ca和区域Cb；第三层分为区域Aa1、区域Aa2、区域Ab1、区域Ab2、区域Ac1、区域Ac2、区域Ba1、区域Ba2、区域Ba3、区域Bb1、区域Bb2、区域Bc1、区域Bc2、区域Ca1、区域Ca2、区域Cb1、区域Cb2和区域Cb3；第四层为第三层各区域内的发电机组；步骤(7)：选出各层的领导者和跟随者智能体，第一层的区域A为领导者，区域B和区域C为跟随者；第二层的区域Aa、区域Ba和区域Ca为领导者，其他的二层区域为跟随者；第三层区域Aa1、区域Ab1、区域Ac1、区域Ba1、区域Bb1、区域Bc1、区域Ca1和区域Cb1为领导者，其他第三层的区域为跟随者；步骤(8)：结合智能体的拓扑结构图以及公式(15)和公式(16)对第一层智能体的一致性变量进行更新；

分布式网络拓扑结构的有向图G＝(V,E)和相应智能体之间权重关系的邻接矩阵A＝(aij)构成多智能体的网络拓扑结构；有向图G的Laplace矩阵L＝[lij]为：其中，lii和lij是Laplace矩阵L中的元素；aij是邻接矩阵A中的元素；

智能体的行随机矩阵Hij为：

一致性变量λ表示为：

λ＝[λT，λClean]   (12)

其中，λ表示综合能源系统的一致性变量；λT表示火力发电机组的一致性变量；λClean表示清洁能源发电机组的一致性变量；

跟随者和领导者的一致性变量更新公式分别为：

其中，NIB表示智能体的个数；k表示第k次迭代； 表示t时刻第n个智能体在第k次迭代的一致性变量； 表示t时刻第m个智能体在第k+1次迭代的一致性变量；τ表示分t布式一致性方法的功率平衡调节因子；ΔP (k)表示t时刻第k次迭代的功率偏差；Hmn(k)表示智能体的行随机矩阵；

步骤(9)：根据公式(17)求出第一层各智能体的有功出力；

各机组的有功功率公式为：

其中， 表示t时刻第i个智能体在k+1次迭代的有功出力； 表示t时刻第j个火电发电机组在k+1次迭代的有功出力； 表示t时刻第u个清洁能源发电机组在k+1次迭代的有功出力； 表示为第ω层第μ区域的有功上限，为该区域所分的下层子区域或机组有功上限之和； 表示t时刻第i个智能体在k+1次迭代的一致性变量； 表示t时刻第j个火电发电机组在k+1次迭代的一致性变量； 表示t时刻第u个清洁能源机组在k+1次迭代的一致性变量； 表示第ω层第μ区域的火电机组的数量； 表示第ω层第μ区域的清洁能源发电机组的数量；U表示该层的区域总数；W表示层总数；

步骤(10)：根据公式(18)对步骤(9)求出的第一层有功出力进行修正；

有功功率的修正公式为：

其中， 表示第i个发电智能体的有功出力下限； 表示第i个发电智能体的有功出力上限；

步骤(11)：根据公式(19)求出第一层的有功功率偏差值；

t

有功功率偏差ΔP(k+1)为：

步骤(12)：判断步骤(11)求出的第一层的多模态有功功率偏差值是否满足精度要求其中，δ表示功率偏差允许最大值；表示多模态功率偏差允许值；

如果 则令迭代步数k＝k+1，并返回步骤(8)，继续进行第一层的迭代计算；

如果 将得到的一组智能体的有功出力存储起来，然后另迭代步数k＝k+1，并返回步骤(8)，继续进行第一层的迭代计算；

如果 则第一层的多模态有功偏差值满足要求，迭代结束，并获得多组智能体的有功出力；

步骤(13)：对步骤(12)得到的多组智能体发电机组的有功出力进行多模态判断和决策拥挤距离筛选，获得满足决策者要求的多模态有功出力；

多模态判断为：

|f(x1)‑f(x2)|<ε   (21)其中，x1和x2就是多模态的两组解决方案；ε表示多模态的临界阀值；

决策空间拥挤距离筛选为：

其中， 表示第i组解的第j个决策变量的拥挤距离；xi+1,j表示第i+1组解的第j个决策变量；xi,j表示第i组解的第j个决策变量； 表示最大目标值的第j个决策变量； 表示最小目标值的第j个决策变量；xi‑1,j表示第i‑1组解的第j个决策变量； 表示第i组解的拥挤距离；D表示决策变量的个数；ξ表示多模态允许的最小决策拥挤距离；

步骤(14)：将步骤(13)得到的第一层的多组多模态发电量作为负荷值输入到得二层，并令初始迭代步数k＝0；

步骤(15)：结合智能体的拓扑结构图以及公式(15)和公式(16)对第二层智能体的一致性变量进行更新；

步骤(16)：根据公式(17)求出第二层各智能体机组的有功出力；

步骤(17)：根据公式(18)对步骤(16)求出的有功出力进行修正；

步骤(18)：根据公式(19)求出第二层有功功率偏差值；

步骤(19)：判断步骤(18)求出的多模态有功功率偏差值是否满足精度要求；

如果 则令迭代步数k＝k+1，并返回步骤(15)，继续进行第二层的迭代计算；

如果 将得到的一组智能体的有功出力存储起来，然后另迭代步数k＝k+1，并返回步骤(15)，继续进行第二层的迭代计算；

如果 则第二层的多模态有功偏差值满足要求，迭代结束，并获得多组智能体发电机组的有功出力；

步骤(20)：对步骤(19)得到的第二层的多组智能体的有功出力进行多模态判断和决策拥挤距离筛选，获得满足决策者要求的多模态有功出力；

步骤(21)：将步骤(20)得到的第二层的多组多模态发电量作为负荷值输入到得三层，并令初始迭代步数k＝0；

步骤(22)：结合智能体的拓扑结构图以及公式(15)和公式(16)对第三层智能体的一致性变量进行更新；

步骤(23)：根据公式(17)求出第三层各智能体机组的有功出力；

步骤(24)：根据公式(18)对步骤(23)求出的有功出力进行修正；

步骤(25)：根据公式(19)求出第三层有功功率偏差值；

步骤(26)：判断步骤(25)求出的多模态有功功率偏差值是否满足精度要求；

如果 则令迭代步数k＝k+1，并返回步骤(22)，继续进行第三层的迭代计算；

如果 将得到的一组智能体的有功出力存储起来，然后另迭代步数k＝k+1，并返回步骤(22)，继续进行第三层的迭代计算；

如果 则第三层的多模态有功偏差值满足要求，迭代结束，并获得多组智能体发电机组的有功出力；

步骤(27)：对步骤(26)得到的第三层的多组智能体的有功出力进行多模态判断和决策拥挤距离筛选，获得满足决策者要求的多模态有功出力；

步骤(28)：将步骤(27)得到的第三层的多组多模态发电量作为负荷值输入到得四层，并令初始迭代步数k＝0；

步骤(29)：结合智能体的拓扑结构图以及公式(15)和公式(16)对第四层各发电机组的一致性变量进行更新；步骤(30)：根据公式(17)求出第四层各发电机组的有功出力；

步骤(31)：根据公式(18)对步骤(30)求出的有功出力进行修正；

步骤(32)：根据公式(19)求出第四层有功功率偏差值；

步骤(33)：判断步骤(32)求出的多模态有功功率偏差值是否满足精度要求；

如果 则令迭代步数k＝k+1，并返回步骤(29)，继续进行第四层的迭代计算；

如果 将得到的一组发电机组的有功出力存储起来，然后另迭代步数k＝k+1，并返回步骤(29)，继续进行第四层的迭代计算；

如果 则第四层的多模态有功偏差值满足要求，迭代结束，并获得多组发电机组的有功出力；

步骤(34)：对步骤(33)得到的第四层的多组发电机组的有功出力进行多模态判断和决策拥挤距离筛选，获得满足决策者要求的多模态有功出力决策方案；

步骤(35)：判断t

步骤(36)：判断η<1是否满足，若满足，则将取值为η时得到的多模态解储存起来，并令η＝η+R，并转入步骤(4)；若不满足，执行下一步；

步骤(37)：将取不同η值得到的所有多模态解反求得到多模态多目标值，所有η值下得到的多模态多目标值在目标空间组成帕累托前沿。