1.一种能源互联网协同优化运行方法,其特征在于,包括以下步骤:1)建立精细化热网模型对于管道k而言,定义t时刻的流入热流量 与流出热流量 如下:考虑t时刻管道k的输送延时τ k,t ,t时刻的流入热流量 与t+τ k,t 时刻的流出热流量 二者满足:τ k,t =F k /G k,t式中:G k,t 为t时刻管道k的流量,C为工质体积比热, 分别为t时刻管道k的入口温度与出口温度,ε k 为管道k的损失系数,F k 为管道k的特征量,由管道长度、截面积等参数决定;热力网络平衡方程为:式中:A 1 ,A 2 分别为热力网络的起点关联矩阵与终点关联矩阵, 分别为各管道在时段t的流入热流量构成的向量与流出热流量构成的向量,Q t 为所有节点在时段t的热流量流入量构成的向量,其中热源节点流入量为正,热负荷节点流入量为负;2)电-热互联综合能源系统中,电网由微型燃气轮机和风力发电设备供电,电能出现不足或富余时,与主网进行购售电交易;热网由微型燃气轮机和燃气锅炉供热;将微型燃气轮机作为热网与电网的耦合载体,建立电-热互联综合能源系统分布式能量优化模型;模型以整个系统的运行成本最小为优化目标,可以得到目标函数为:式中:f为系统运行总成本,由电网运行成本和f e 热网运行成本f h 构成, 分别为第t个调度周期时刻对应的购售电价, 分别为第t个调度周期内电网购电功率与售电功率,Δt为调度周期, 为第t个调度周期时刻对应的天然气价格, 为第t个调度周期内热网购气量,分别用于微型燃气轮机和燃气锅炉;电网约束条件考虑电力平衡与电力网络约束,其中电力平衡约束如下:式中:η T 为变压器效率, 分别为第t个调度周期内微型燃气轮机发电量和风力发电量, 为第t个调度周期内电负荷用电量;热网约束条件除了考虑上述热网精细化模型中建立的热力支路特性方程,还需考虑热力平衡约束、热力网络约束以及微型燃气轮机、燃气锅炉的设备运行约束;其中热力平衡约束如下:式中: 分别为第t个调度周期内微型燃气轮机、燃气锅炉产热量, 为第t个调度周期内热负荷耗热量;微型燃气轮机作为电-热耦合元件,具有如下耦合方程:式中,υ为分配系数,即用于微型燃气轮机的天然气比例,η E,MT 为微型燃气轮机的电转化效率;利用上述耦合载体对热网模型和电网模型进行耦合,形成电-热互联综合能源系统分布式能量优化模型;min f e (x)+f h (z)式中:区域e代表电网,区域h代表热网,x、z分别为电网、热网内部优化变量,f e (x)、f h (z)分别为电网、热网进行独立优化的目标函数,h e (x)=0、h h (z)=0分别为电网与热网内部等式约束,g e (x)≤0、g h (z)≤0分别为电网与热网内部不等式约束,x B 、z B 分别为电网、热网的边界变量,由耦合方程可知, 3)采用同步型ADMM算法求解电-热互联综合能源系统分布式能量优化模型;初始化电网内优化变量x 0 ,热网内优化变量z 0 ,并分别初始化电网和热网内广义拉格朗日乘子u x 0 和u z 0 ,置权重值λ=0.5,惩罚因子ρ=1.5,初始化固定参考值初值 和 根据下式进行电网与热网区域的优化迭代计算,得到优化变量的更新值x t+1 ,z t+1 ;将电网边界变量更新值 传输给热网,并将热网边界变量更新值 传输给电网;分别计算电网和热网下次迭代需要的参考固定值 和本区的广义拉格朗日乘子的更新值 通过收敛判据公式判断当前迭代优化得到的解x与z是否满足精度要求,如满足则迭代结束,输出优化结果,如不满足预定的收敛精度,那么令t=t+1,并且继续迭代计算。
2.根据权利要求1所述能源互联网协同优化运行方法,其特征在于:所述固定参考值的设置依据为
3.根据权利要求2所述能源互联网协同优化运行方法,其特征在于:广义拉格朗日乘子的设置依据为
4.根据权利要求3所述能源互联网协同优化运行方法,其特征在于:所述收敛判据公式为