[发明专利]一种面向电-热综合能源系统的量子化事件驱动仿真方法

权利要求书

1.一种面向电-热综合能源系统的量子化事件驱动仿真方法，其特征在于，包括如下步骤：

1)针对待仿真的电-热综合能源系统，输入系统的拓扑连接关系、元件参数、控制参数、仿真计算参数，进行仿真初始化；其中控制参数包括热源出口温度控制区间、用户室内温度控制区间，仿真计算参数包括仿真终止时间、量子化积分阈值、积分误差容限、电压变化阈值、各仿真变量及仿真变量历史量的初值、热力状态变量的一阶、二阶、三阶导数初值、量子化热力状态变量的初值和一阶、二阶导数初值；将所有变量的上次更新时刻和下次更新时刻置为0，置当前仿真时刻t＝0；

2)建立待仿真的电-热综合能源系统量子化耦合模型，包括热力环节微分方程组、热力控制环节离散代数方程组、电力环节微分代数方程组、耦合接口代数方程组；

3)生成仿真事件时刻表T_list＝[t^h,s,t^h,d,t^e,t^he,t^eh]，其中t^h,s为热力状态变量的下次更新时刻，t^h,d为热力离散代数变量的下次更新时刻，t^e为电力状态变量x^e和电力代数变量y^e的下次更新时刻，t^he为热-电耦合接口变量的下次更新时刻，t^eh为电-热耦合接口变量的下次更新时刻；取下一步最先发生的仿真事件时刻t＝min{T_list}；

4)若当前仿真时刻t＝min{t^h,s}，则采用三阶量子化状态系统积分方法将热力环节微分方程组积分至当前仿真时刻t，并更新热力状态变量x^h的下次更新时刻t^h,s、热力离散代数变量u^h的下次更新时刻t^h,d，执行步骤5)；否则直接执行步骤5)；

5)若当前仿真时刻t＝t^h,d，则根据热力控制环节离散代数方程组计算当前仿真时刻t的热力离散代数变量u^h(t)，并更新热力离散代数变量u^h的下次更新时刻t^h,d、热-电耦合接口变量u^he的下次更新时刻t^he，执行步骤6)；否则直接执行步骤6)；

6)若当前仿真时刻t＝t^e，则采用变步长积分算法将电力环节微分代数方程组积分至当前仿真时刻t，并更新电力状态变量x^e和电力代数变量y^e的下次更新时刻t^e、电-热耦合接口变量u^eh的下次更新时刻t^eh，执行步骤7)；否则直接执行步骤7)；

7)若当前仿真时刻t＝min{t^he}，则根据耦合接口代数方程组计算当前仿真时刻t的热-电耦合接口变量u^he(t)，并更新电-热耦合接口变量u^eh下次更新时刻t^he，执行步骤8)；否则直接执行步骤8)；

8)若当前仿真时刻t＝min{t^eh}，则根据耦合接口代数方程组计算当前仿真时刻t电-热耦合接口变量u^eh(t)，并更新u^eh(t)下次更新时刻t^eh，执行步骤9)；否则直接执行步骤9)；

9)判断当前仿真时刻t是否达到输入的仿真终止时间T，若t≥T，则仿真结束，输出仿真结果；否则返回步骤3)。