和负荷平衡Agent。
[0075] 能量路由器的经济调度模型如式（2)所示：
[0076]
[0077] 其中，Totalcost为总费用，ae(t)为实时的电费用，Pe(t)为t时刻公共电网输 入电功率，ag为天然气费用，Pg(t)为t时刻天然气输入功率，a。为石油费用，Pjt)为t 时刻石油输入功率，if(〇为t时刻电能存储装置的充电功率，if5(〇为t时刻电能存储装 置的放电功率，af为电能存储装置操作费用，<为热能存储装置操作费用，EENS°为能量 路由器的电能负荷损失能量，惩罚成本系数，取30美分/KWh，aDR为用电补偿费用， ^(0为t时刻内增加的电能负荷功率，/T⑴为t时刻内中断的电能负荷功率，T= 24h为 总时间。
[0078] 能量路由器模型的多智能体系统的约束条件分别为：
[0079] 电能Agent的约束条件为：当前时刻公共电网输入电功率在允许公共电网输入电 功率的最小值和最大值之间；当前时刻公共电网输出的电能负荷功率为输电网稳定性概 率、变压器转换效率和当前时刻公共电网输入电功率的乘积。
[0080] 本实施方式中，电能Agent的约束条件如式（3)所示：
[0081]
[0082] 其中，ZT⑴为t时刻公共电网输出的电能负荷功率，Sn6t为输电网络稳定性概率， <为变压器转换效率，Pe(t)为t时刻公共电网输入电功率，Pemin(t)为t时刻允许公共电 网输入电功率的最小值，P广x(t)为t时刻允许公共电网输入电功率的最大值。
[0083]本实施方式中，Snet= 0? 98，<=0.95,Pe_(t) = 1500kw，Pe_(t) = -200kw。
[0084] 风力发电Agent的约束条件为：当前时刻风力发电输入电功率在允许风力发电输 入电功率的最小值和最大值之间；当前时刻风力发电输出的电能负荷功率为风力发电机组 稳定性概率、AC/AC转换器的转换效率和当前时刻风力发电输入电功率的乘积。
[0085] 本实施方式中，风力发电Agent的约束条件如式（4)所示：
[0086] (43：
[0087]其中，为t时刻风力发电输出的电能负荷功率，风力发电机组稳定性 概率，C"为AC/AC转换器转换效率，Pw⑴为t时刻风力发电输入电功率，Pwmin⑴为t时刻 允许风力发电输入电功率的最小值，Pwmax(t)为t时刻允许风力发电输入电功率的最大值。
[0088]本实施方式中，Swind= 0? 95, _=〇爲:Pw_(t) = 500kw，Pw_(t) = 0。
[0089] 光伏发电Agent的约束条件为：当前时刻光伏发电输入电功率在允许光伏发电输 入电功率的最小值和最大值之间；当前时刻光伏发电输出的电能负荷功率为光伏列阵发电 稳定性概率、DC/AC转换器转换效率和当前时刻光伏发电输入电功率的乘积。
[0090] 本实施方式中，光伏发电Agent的约束条件如式（5)所示：
[0091]
(3)
[0092] 其中，1广(〇为t时刻光伏发电输出的电能负荷功率，Ss^为光伏阵列发电稳定 性概率，为DC/AC转换器转换效率，Ps(t)为光伏发电输入电功率，Psmin(t)为允许光伏 发电输入电功率的最小值，Psmax(t)为允许光伏发电输入电功率的最大值。
[0093]本实施方式中，Ss°lar=0?95，C"=0.85,Ps_(t) =0,Ps_(t) = 450kw。
[0094] 石油Agent的约束条件为：当前时刻石油输入功率在允许石油输入功率的最小值 和最大值之间；当前时刻石油输出的负荷功率为输油管道稳定性概率、当前时刻石油用于 石油用户负载的调度参数和当前时刻石油输入功率的乘积；当前时刻石油用于石油用户负 载的调度参数和当前时刻石油转换为热能的调度参数之和为1 ;当前时刻石油用于石油用 户负载的调度参数大于等于〇且小于等于1。
[0095] 本实施方式中，不考虑功率损耗，石油Agent的约束条件如式（6)所示：
[0096]
[0097] 其中，Lftans⑴为t时刻石油输出的负荷功率，ST?S为输油管道稳定性概率，v。3⑴ 为t时刻石油用于石油用户负载的调度参数，VciF(t)为t时刻石油转换为热能的调度参数， Pjt)为t时刻石油输入功率，P^t)为t时刻允许石油输入功率的最小值，P^t)为t 时刻允许石油输入功率的最大值。
[0098]本实施方式中，STrans= 0? 98,PDmin(t) = 0,PDmax(t) = 800kw〇
[0099] 热电联产设备Agent的约束条件为：当前时刻热电联产设备电输出功率在允许热 电联产设备电输出功率的最大值之下；当前时刻天然气输入功率在允许天然气输入功率的 最小值和最大值之间；当前时刻热电联产设备输出的电能负荷功率为热电联产设备运行稳 定性概率、热电联产设备的天然气转换为电能的转换效率、当前时刻天然气转换为电能的 调度参数和当前时刻天然气输入功率的乘积；当前时刻热电联产设备输出热能负荷功率为 热电联产设备运行稳定性概率、热电联产设备的天然气转换为热能的转换效率、当前时刻 天然气转换为电能的调度参数和当前时刻天然气输入功率的乘积；当前时刻天然气转换为 电能的调度参数大于等于〇且小于等于1。
[0100] 本实施方式中，热电联产设备Agent的约束条件如式（7)所示：
[0101]
[0102] 其中，If(0为t时刻热电联产设备输出的电能负荷功率，SeHP为热电联产设备运 行稳定性概率，为热电联产设备的天然气转换为电能的转换效率，Vg~(t)为t时刻天 然气转换为电能的调度参数，pg(t)为t时刻天然气输入功率，If(0为t时刻热电联产设 备输出的热能负荷功率，为热电联产设备的天然气转换为热能的转换效率，为热电 联产设备允许最大电输出功率，为t时刻允许天然气输入功率的最小值，pgmax(t)为t时刻允许天然气输入功率的最小值。
[0103] 本实施方式中，SCHP=0.9,f=〇35, %; :=〇45,Pchp= 500，Pgmin(t) = 150kw， Pgmax(t) = 1800kw〇
[0104] 加热装置Agent的约束条件为：当前时刻加热装置将天然气转换为热能的热能负 荷功率为加热装置设备运行稳定性概率、加热装置的天然气转换为热能的转换效率、当前 时刻天然气转换为热能的调度参数和当前时刻天然气输入功率的乘积；当前时刻加热装置 将石油转换为热能的热能负荷功率为加热装置设备运行稳定性概率、加热装置的石油转换 为热能的转换效率、当前时刻石油转换为热能的调度参数和当前时刻石油流向加热装置的 输入功率的乘积；当前时刻天然气转换为电能的调度参数大于等于〇且小于等于1 ;当前时 刻石油转换为热能的调度参数大于等于〇且小于等于1，当前时刻天然气转换为电能的调 度参数和当前时刻天然气转换为热能的调度参数之和为1。
[0105] 本实施方式中，加热装置Agent的约束条件如式（8)所示：
[0106]
[0107]其中，Zf⑴为t时刻加热装置将天然气转换为热能的热能负荷功率，#为加热装 置设备运行稳定性概率，为加热装置的天然气转换为热能的转换效率，vgF(t)为t时刻 天然气转换为热能的调度参数，If的为t时刻石油流向加热装置的输入功率为t时 刻加热装置将石油转换为热能的热能负荷功率，'为加热装置的石油转换为热能的转换 效率，VtfU)为t时刻石油转换为热能的调度参数。
[0108]本实施方式中，SF= 0.95，<=〇.85, d=〇95i:
[0109] 电能存储装置Agent的约束条件为：当前时刻电能存储装置的充放电功率平衡； 当前时刻电能存储装置存储功率在电能存储装置存储功率最小值和最大值之间；当前时刻 电能存储装置的充电功率在电能存储装置充电功率最小值和最大值之间；当前时刻电能存 储装置的放电功率在电能存储装置放电功率的最小值和最大值之间；当前时刻电能存储装 置的充电状态变量和放电状态变量之和大于等于〇且小于等于1。
[0110] 本实施方式中，电能存储装置Agent的约束条件如式（9)所示：
[0111]
[0112] 其中，P/(〇为t时刻电能存储装置电能存储功率，/f(?-1)为t-1时刻电能存储 装置电能存储功率，/f(f)为t时刻电能存储装置充电功率，P/IO为t时刻电能存储装 置放电功率，为电能存储装置最大存储功率，为电能存储装置最小存储率，//f*为 电能存储装置最大存储率，为电能存储装置充电效率，为电能存储装置放电效率， Sf切为t时刻电能存储装置的充电状态变量，为0或1的整型变量，(〇为t时刻 电能存储装置的放电状态变量，5："'为0或1的整型变量。
[0113]本实施方式中，#=lft0kw，//fmin=0J7,爲_=0斑,：#=0,終
[0114] 热能存储装置Agent的约束条件为：当前时刻热能存储装置的充放热功率平衡； 当前时刻热能存储装置存储功率在热能存储装置存储功率最小值和最大值之间；当前时刻 热能存储装置的充热功率在电能存储装置充热功率最小值和最大值之间；当前时刻热能存 储装置的放热功率在热能存储装置放热