[发明专利]一种基于实时碳排放控制的电气互联系统低碳调度方法

权利要求书

1.一种基于实时碳排放控制的电气互联系统低碳调度方法，其特征在于，包括如下步骤：

步骤1，建立零和博弈的低碳激励机制；

步骤2，构建碳排放P2P交易框架；

步骤3，构建电气互联系统双层组合优化模型，上层模型考虑电气互联系统的实时碳排放总量控制，以电气互联系统低碳经济运行成本最小为目标建立目标函数，包括经济成本最优和碳排放控制效果最优，约束条件包括电力网络约束和天然气网络约束，下层模型以低碳激励机制下多能园区的运行成本和环境成本最低为目标，约束条件包括电、热、气功率平衡约束、电储能约束、储热约束、能源转换约束及多能园区的碳排放约束；

步骤4，对双层模型进行迭代运算，获取考虑P2P交易的低碳经济调度结果。

2.根据权利要求1所述的一种基于实时碳排放控制的电气互联系统低碳调度方法，其特征在于，步骤1建立零和博弈的低碳激励机制的具体方法如下：

在低碳激励机制下，设定所有多能园区碳配额的总和为定值，如下式

式中，Ω_EH(j)为多能园区的集合；E_h，t为多能园区h在时刻t的碳配额；为所有多能园区的碳配额总和；为碳管理商所获碳收益；为碳市场碳交易价格；S_h，t为多能园区h在时刻t的实际碳排放；β_t和Υ_t分别为碳排放奖励因子以及惩罚因子。

3.根据权利要求1所述的一种基于实时碳排放控制的电气互联系统低碳调度方法，其特征在于，步骤2构建碳排放P2P交易框架的步骤为：

步骤2.1，数据初始化，电气互联系统给各多能园区下发初始电价、气价、碳配额信息；

步骤2.2，内部工序平衡，具体多能园区h内部供需平衡目标函数为：

式中，为多能园区h的成本函数；P_h，t表示时刻t配网给多能园区h提供的有功功率；为市场上的电价；G_h，t表示多能园区h的耗气量；为气价；和为多能园区内部自平衡优化形成的对外购售电功率、热功率、碳配额的需求；和分别为多能园区从其他多能园区购售电功率、热功率、碳配额的价格；

步骤2.3，各多能园区以最大化预期经济效益、最小化预期经济成本为目标确定报价/投标价格，其最优竞价目标函数为：

式中，Ψ_e，t为多能园区作为卖方的收益函数，和为待优化的电、热、碳报价价格；和ρ、和h、和s分别为电、热、碳报价价格的上下限；和为卖方多能园区选择买方多能园区的预期交易电功率、热功率、碳配额；为市场上的电价、热价；

在竞价过程中需要满足的约束条件为：

步骤2.4，卖方和买方最佳匹配的社会福利ψ_e，f，t可表示为：

式中，和为买方实际成交的电功率、热功率、碳配额；为卖方实际成交的电功率、热功率、碳配额；

电功率、热功率与碳配额交易需满足的约束类似，电功率交易约束如下：

步骤2.5，交割，作为卖方的多能园区需根据交易结果发电，配网负责将该电量配送到对应的买方处，当买方和卖方出现交易偏差时，由配网承担兜底责任。

4.根据权利要求1所述的一种基于实时碳排放控制的电气互联系统低碳调度方法，其特征在于，步骤3构建电气互联系统双层组合优化模型的具体步骤为：

步骤3.1，上层优化模型考虑电气互联系统低碳经济运行成本最小，包括经济成本最优和碳排放控制效果最优，关于经济成本的目标函数可表示为：

式中：右边第一项为电气互联系统在电力市场上的购电成本，Ω_UP(j)表示与节点j相连的上级电网节点的集合，P_x，t表示时刻t上级电网节点向节点j注入的有功；右边第二项为电气互联系统在天然气市场上的购气成本，Ω_NW(m)表示与节点m相连的天然气井的集合，G_w，t表示天然气井w提供的天然气流量；右边第三项为配网向能源集成商出售电能、天然气的获取的收益，和分别为配网向能源集成商售电、售气的价格；

系统碳排放控制效果可利用系统实时碳排放目标与实际碳排放之间的关系进行建模，可表示为：

式中：为电气互联系统的实时低碳排放目标值；为单位电功率的碳排放强度；为单位天然气的碳排放强度；

步骤3.2，根据电气互联系统由配电网和配气网构成，建立系统约束，系统约束包括电力网络约束和天然气网络约束，具体如下：

电力网络约束：

式中：Ω_PL表示电力网络中所有支路的集合；Ω_OL表示含有载调压变压器支路的集合；P_j，t和P_ij，t分别表示时刻t由节点j→k和节点i→j的有功；r_ij表示线路ij的电阻；I_ij，t表示时刻t流经线路ij电流I_ij，t的平方；P_j，t表示时刻t节点j处的有功负荷；Ω_wid(j)、Ω_PV(j)、Ω_GU(j)、Ω_PG(j)、Ω_EH(j)分别表示与节点j相连的上级电网节点、风机、光伏、燃气机组、电转气设备、EH的集合；P_a，t、P_b，t和P_u，t分别表示时刻t风机a、光伏b和燃气机组u向节点j注入的有功；P_g，t表示时刻t电转气设备g消耗的有功；Ω_B表示配电网中所有节点的集合；

式中：Q_j，t和Q_ij，t分别表示时刻t从节点j→k和节点i→j的无功；x_ij表示线路ij的电抗；Q_j，t表示时刻t节点j处的无功负荷；Q_a，t表示时刻t风机a吸收的无功；Q_c，t表示时刻t静止无功补偿装置c提供的无功；Ω_SV(j)表示与节点j相连SVC装置的集合。

式中：v_i和v_j表示节点i和节点j的电压的平方；Z_ij和S_ij表示支路ij的阻抗和视在功率；表示SVC装置的可调容量上限；

天然气网络约束：

天然气节点m的能量平衡约束可表示为：

式中：Ω_PG表示与节点m相连的P2G设备的集合；G_g′，t表示P2G设备g′提供的天然气流量；Ω_EH(m)表示与节点m相连的多能园区的集合；G_m，t表示节点m的天然气负荷；Ω_GU(m)表示与节点m相连的燃气机组的集合；G_u′，t表示燃气机组u′的耗气量；ψ_b(m)、ψ_f(m)分别表示与节点m相连的上游节点、下游节点的集合；G_mn′，t和G_mn，t表示由节点n′→m和节点m→n的天然气流量；G_mn，t表示管道mn的天然气流量；Ω_G(mn)表示管道mn中的压缩机的集合；α_c′表示压缩机c′的燃料消耗系数；Ω_GB表示所有天然气节点的集合；

天然气管道潮流约束可表示为：

式中：表示管道mn的流量上限。

假定天然气管道中潮流的流向是由节点m流向节点n，则有：

G_u，t＝P_u，t/(η_u·I_gas)

G_g，t＝P_g，t·η_g/I_gas

式中：γ_c′为天然气压缩机c′的压缩因子；η_u为GU的气转电效率；η_g为P2G的电转气效率；I_gas为天然气的燃烧热值；

步骤3.3，建立下层优化模型目标函数

下层优化以能源集成商运行成本最小为目标，可表示为：

步骤3.4，设置下层优化模型约束条件；

系统的电、热、气功率平衡约束可表示为：

式中：为风电机组出力；为光伏出力；表示EH中CHP机组提供的有功功率；和分别表示电储能装置的充、放电功率；表示EH中电锅炉消耗的电功率；表示EH的电负荷需求；

式中：和分别表示CHP机组、燃气锅炉、电锅炉提供的热功率；和表示储热装置的储热功率、放热功率；表示EH的热负荷需求；

式中：和表示CHP机组和燃气锅炉所需的天然气量；

电储能装置的模型及约束条件：

式中：η_EC、η_E表示电储能装置的充放电效率；表示电储能装置蓄电量的始末值；

储热装置的模型及约束条件：

式中，η_HC、η_H表示储热装置的储热、放热效率；表示储热装置蓄热量的始末值；

能源转换装置的模型及约束条件：

式中：η_CE、η_EH和η_G分别表示CHP机组的电效率、电锅炉的效率和燃气锅炉的效率；

多能园区的碳排放模型及约束条件：