1.一种在不确定环境下考虑综合需求响应的园区综合能源系统优化调度方法，其特征是，它包括以下步骤：

1)构建园区综合能源系统物理模型；参见图2，风机、光伏、储能装置(energy storage device,ESD)、外部电网和燃气轮机(micro‑gas turbine,MT)联合为用户提供电需求，电锅炉(electric boiler,EB)、MT和储热装置(heat storage device,HSD)联合为用户提供相应的热需求，外部气网和电转气装置(power to gas,P2G)联合为用户提供相应的天然气需求。

2)建立考虑电‑气‑热柔性负荷和电动汽车的综合需求响应机制；

①电力需求响应

电力负荷由固定负荷和柔性负荷组成，根据需求响应的特点，电力柔性负荷分为可转移负荷和可中断负荷；

a)可时移的电力负荷

可时移负荷的特征是总耗电量是恒定的，并且消耗时间可灵活改变，用式(1)、式(2)描述:

TSE

式中，Pt 是时移的电力负荷，αTSE表示可时移负荷所占比例， 和 是t时段可时移负荷的上下限；

b)可中断电力负荷

在电力供应不足或电价高的时期，用户可中断部分负荷以缓解电力供应压力，可中断负荷的相关约束用式(3)描述:EIL

式中，Pt 和 分别是t时段中断的电力负荷及其上限值，αIE表示可中断负荷所占比例；

②热需求响应

将建筑热需求视为热负荷，利用暂态热平衡方程将建筑温度与热需求联系起来，采用预测平均投票指标(predicted mean vote,PMV)来描述用户对室内温度变化的舒适体验，供热可中断负荷的相关约束为式(4):式中， 是中断的热负荷， 是在t时段的中断热负荷最大值。通过引入PMV来描述用户可接受的热舒适范围；

式中，M为人体能量代谢率；Icl是服装的热阻；Ts是处于舒适状态的人体皮肤的平均温度；Tin,t是室内温度；

室内温度变化范围是

建筑物的热需求可由下式计算得到：式中，Tout,t为t时段室外温度；K为建筑物综合传热系数；F为建筑物表面积；V为建筑物体积；cair为室内空气的比热容；ρair为室内空气的密度；

③天然气需求响应

天然气需求响应与电力需求响应相似，天然气需求响应负荷可由下式表示：式中， 和 分别表示t时段可时移气负荷及其上下限， 和表示t时段可中断气负荷及其上限值，αTSQ和αIQ分别为可时移气负荷和可中断气负荷所占比例；

④电动汽车模型

电动汽车到达充电站时间的概率密度函数为：式中，μs和σs分别为电动汽车到达充电站时间的平均值和标准差；

电动汽车的日行驶里程服从正态分布，其概率密度函数为：式中，x表示电动汽车的日行驶里程，μd和σd分别为日行驶里程的平均值和标准差；

根据电动汽车的行驶里程，电动汽车的初始充电状态为：式中， 为电动汽车n的初始电池容量，W100表示电动汽车100公里耗能， 为电动汽车额定容量；电动汽车的充电时间计算为：式中，Tn为电动汽车的充电时间， 和 为电动汽车的额定充电功率和效率；

为了衡量IDR对用户体验的影响，设计了用户的综合满意度指标：式中，ms,t为用户的综合满意度， 和Pload,t为t时段初始电力负荷和实际电力负荷，Hload,t和 为t时段初始热负荷和实际热负荷，Qload,t和 为t时段初始天然气负荷和实际天然气负荷；

2)建立基于机会约束规划的园区综合能源系统优化调度模型；

①选取优化目标，模型以最小化CIES总运行成本作为目标函数，目标函数的表达式为：式中，C1为能源交易成本，C2为旋转备用成本，C3为运行维护成本，C4为环境成本，C5为IDR补偿成本， 和 为分时电价和天然气价格， 为电力负荷消耗电网功率， 为气ESSD

负荷消耗气网天然气体积， 和δ 分别表示为旋转备用成本系数， 和 表示电网和EESD提供的旋转备用容量，表示设备i的运维成本，Pi,t表示设备i的功率， 和 表示所购电力和天然气中污染物j的排放系数， 表示P2G过程中吸收CO2的系数， 是污染物j的价格，γIE、γTSE、γCH、γIQ和γTSQ分别为可中断电负荷、可时移电负荷、可中断热负荷、可中断气负荷和可时移气负荷的单位补偿成本。

②确定约束条件，调度模型的约束条件包括额能量平衡约束、储能设备运行约束、电锅炉约束、P2G运行约束、燃气轮机运行约束、旋转备用约束，具体如下：能量平衡约束：包括电力供需平衡、热功率约束和天然气平衡约束EV s

式中， 为t时段电力负荷消耗RG出力，Pt 为t时段电动汽车充电功率，Pt为弃风弃光功率；

储能设备约束：包括储能装置功率约束和容量约束，式中，Ct+1和Ct分别为t+1和t时段储能容量，PCH,max和PDC,max表示储能装置的最大充放功率，Cmin和Cmax分别为最大和最小储能容量，C0和 分别为一个调度周期内初始和结束时的储能容量；

电锅炉运行约束：

EB EB

式中，Pt 和 表示耗电功率和供热功率，η 为电锅炉效率， 和 为电锅炉消耗电网和可再生能源功率， 为额定供热功率；

P2G运行约束：

P2G P2G

式中，Pt 为t时段耗电功率，η 为P2G效率系数， 和 分别为P2G消耗电网功率和可再生能源功率， 为t时段P2G产气体积，HHV为天然气热值， 和 为P2G输入功率上下限，θt为P2G运行状态变量， 和 为爬坡功率上下限；

燃气轮机运行约束：

式中， 为MT输入天然气的体积， 和 分别为MT供电和供热效率系数， 为热损MT

失系数，Pt 和 分别为MT供电和供热功率， 和 为MT消耗气网和P2G气体体积，和 为MT输入气体体积上下限，ψt为MT运行状态变量， 和 为MT爬坡能力上下限；

旋转备用约束：

DG

式中， 为电网提供的最大旋转备用，E(Pt )为可再生能源联合处理期望值，α为置信WT PV

水平；Pt 和Pt 分别为风机和光伏出力；

4)利用序列运算理论将光伏、风机输出功率的概率分布进行离散化处理，得到其对应的概率性序列分别为a(iat)和b(ibt)。

5)通过概率性序列，获得各时段风光联合出力的期望值；t时段预测的间歇性风光共同DG

出力的期望值E(Pt )计算式为：式中，Nat为光伏出力概率序列长度，Nbt为风机出力概率序列长度，q为离散化步长，uatq为光伏第ua种状态的出力值，ubtq为风机第ub种状态的出力值；

6)对旋转备用的机会约束形式进行转化，得到确定性约束形式；t时段风光共同出力所对应的概率性序列c(ict)可以利用概率性序列a(iat)和b(ibt)的卷和获得，根据卷和的定义有：

为了方便处理旋转备用约束，定义一类新的0‑1变量 它满足以下关系：在t时段，当系统旋转备用容量大于风光出力期望值与风光第uct种出力的差值时取1，否则为0；

因此旋转备用的机会约束形式可简化为：通过上述方法，我们将基于机会约束规划的模型转换为确定性模型；

7)利用线性化方法对调度模型中的分段函数和最小化因子进行处理，将模型转换为混合整数线性规划形式；

①由于式(38)中 的表达式不兼容混合整数规划的求解形式，需要对其进行转化：式中χ是很大的正数，当 时，式(40)等价为 χ是一个非常小的正数，由于 是一个0‑1变量，所以只能等于1，否则为0。

②消除最小化因子

TSE

为了处理目标函数中的最小化因子，我们引入了辅助变量，以min{Pt ,0}为例，定义TSE

所以g(Pt )可以表示为：通过引入三个连续变量w1,w2,w3和三个0‑1变量z1,z2,z3，式(41)可转换为线性形式：利用相同方法对 进行处理，确定性模型将转换为混合整数线性规划模型；

8)输入初始参数包括：风机参数，光伏组件参数，微燃机参数，电转气设备参数，建筑物参数，电动汽车参数，储能设备参数，电锅炉参数，调度时段数，各负荷预测值以及各变量的上下限值；

9)设置备用容量的置信水平；

10)采用CPLEX求解器对园区综合能源系统调度模型进行求解；

11)检查解决方案是否存在：若存在，则终止流程；否则更新置信度，转到步骤10)再次进行求解；

12)输出园区综合能源系统优化调度方案。