[发明专利]热电耦合系统连锁故障仿真方法和系统

权利要求书

1.一种热电耦合系统连锁故障仿真方法，其特征在于，包括：

S1、获取t时刻的自然灾害气象数据和系统负荷数据；

S2、基于所述自然灾害气象数据计算出线路的故障概率，并更新热电设备的性能状态；

S3、基于所述线路的故障概率随机断开线路，更新热电耦合系统网络的拓扑结构，并确定当前的运行方式；

S4.1、基于所述当前的运行方式，进行热电耦合潮流计算；

S4.2、分析热电耦合潮流计算结果，统计出功率越限的线路；若无功率越限线路，执行S5；

S4.3、对功率越限的线路以预设概率随机断开，更新系统的网络拓扑，确定当前的运行方式之后执行所述S4.1；

S5、统计t时刻内总的线路开断情况，并基于所述系统负荷数据，确定热电耦合系统中的负荷损失。

2.根据权利要求1所述的热电耦合系统连锁故障仿真方法，其特征在于，所述S5之后，还包括；

S6、判断t是否等于t_d，若否，令t＝t+1，并执行所述S1，

其中，t_d为预设的自然灾害发生的最后时刻，1表示自然灾害发生的过程中的预设步长。

3.根据权利要求1所述的热电耦合系统连锁故障仿真方法，其特征在于，所述热电设备包括：电源设备、冷热源设备、储能设备、配电系统传输设备、供热系统传输设备、热电负荷侧的设备中的至少一种。

4.根据权利要求2所述的热电耦合系统连锁故障仿真方法，其特征在于，所述S4.1至S4.3表示的短时间尺度过程的马尔可夫状态空间为：

其中，t_m表示短时间尺度的时刻，L、G、M分别表示短时间尺度过程的线路、机组和储能状态。

5.根据权利要求4所述的热电耦合系统连锁故障仿真方法，其特征在于，所述S6表示的长时间尺度过程的马尔可夫状态空间为：

其中，t_h表示长时间尺度的时刻，W、D、γ分别表示长时间尺度过程的线路、机组和储能状态。

6.根据权利要求5所述的热电耦合系统连锁故障仿真方法，其特征在于，所述长时间尺度过程的马尔可夫状态空间的状态转移概率可由下式给定：

其中n为短时间尺度中的状态转移次数，表示t_m,1时刻的线路状态转移概率。

7.根据权利要求1所述的热电耦合系统连锁故障仿真方法，其特征在于，所述当前的运行方式包括热电是否解耦、热电耦合系统是否分岛。

8.一种热电耦合系统连锁故障仿真系统，其特征在于，包括：

数据获取模块，用于获取t时刻的自然灾害气象数据和系统负荷数据；

状态更新模块，用于基于所述自然灾害气象数据计算出线路的故障概率，并更新热电设备的性能状态；

运行方式确定模块，用于基于所述线路的故障概率随机断开线路，更新热电耦合系统网络的拓扑结构，并确定当前的运行方式；

短时间尺度过程模块，所述短时间尺度过程模块包括：

耦合潮流计算单元，用于基于所述当前的运行方式，进行热电耦合潮流计算；

线路统计单元，用于分析热电耦合潮流计算结果，统计出功率越限的线路；若无功率越限线路，跳转至负荷损失确定模块；

网络更新单元，用于对功率越限的线路以预设概率随机断开，更新系统的网络拓扑，确定当前的运行方式之后跳转至所述耦合潮流计算单元；

还包括负荷损失确定模块，用于统计t时刻内总的线路开断情况，并基于所述系统负荷数据，确定热电耦合系统中的负荷损失。

9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至7任一项所述的热电耦合系统连锁故障仿真方法的步骤。

10.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1至7任一项所述的热电耦合系统连锁故障仿真方法的步骤。