一种综合能源系统RLC暂态模型的能流仿真方法和装置
技术领域
[0001] 本
发明涉及综合能源系统技术领域，具体涉及一种综合能源系统RLC暂态模型的能流仿真方法和装置。
背景技术
[0002] 综合能源系统能够完成电、气、热\冷等各类异质
能量的耦合、转化和分配，打破传统能源系统各自独立运行的状态，实现能量的
梯级利用，有效提升能源供应的灵活度，提高能源综合利用率和经济收益。
[0003] 能流仿真是开展综合能源系统运行状态评估、多异质能流协同优化的重要
支撑。当前，综合能源系统的能流仿真方法主要沿用传统电、气、热\冷等各领域典型的能量潮流计算方法，如电
力系统多采用PQ分解法，
流体系统多采用达西公式和Weymouth公式结合
牛顿网孔-
节点法，热力系统则多采用苏霍夫公式结合工质流法。随后通过能量耦合元件将信息传递给相邻的系统以完成综合能源系统能流计算。由于不同系统之间的耦合性不强，现有综合能源能系统流仿真技术存在以下问题：
[0004] 1)异质能流系统的建模过程复杂，仿真难度大，所建模型的可信度低，尤其是综合能源系统规模较大时，求解难度急剧增大；
[0005] 2)不同异质能流系统在时间尺度存在较大差异，开展综合能源系统动态过程数值仿真时难以协调
算法步长，计算
精度低，收敛速率慢。
发明内容
[0006] 为了克服上述
现有技术中建模仿真难度大、可信度低、计算精度低且收敛速率慢的不足，本发明提供一种综合能源系统RLC暂态模型的能流仿真方法，根据预先构建的综合能源系统RLC暂态模型，确定动态仿真初值；根据动态仿真初值，对综合能源系统RLC暂态模型进行仿真得到综合能源系统中各支路首末端参数；通过支路的首末端参数，确定支路的能流
波动；综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建，降低了建模仿真难度，提高了可信度和计算精度，且收敛速率快。
[0007] 为了实现上述发明目的，本发明采取如下技术方案：
[0008] 一方面，本发明提供一种综合能源系统RLC暂态模型的能流仿真方法，包括：
[0009] 根据预先构建的综合能源系统RLC暂态模型，确定动态仿真初值；
[0010] 根据所述动态仿真初值，对综合能源系统RLC暂态模型进行仿真得到综合能源系统中各支路首末端参数；
[0011] 通过所述支路的首末端参数，确定支路的能流波动；
[0012] 所述综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建。
[0013] 所述综合能源系统包括流体系统、热力系统和电力系统中的任意两种至全部，所述RLC元件模型包括流体系统的RLC元件模型、热力系统的RLC元件模型和电力系统的RLC元件模型中的任意两种至全部；
[0014] 所述综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建，包括：
[0015] 根据RLC元件模型和能量转化元件模型，描述综合能源系统的能量耦合设备；
[0016] 根据综合能源系统的能流路线连接各能量耦合设备，构成综合能源系统的网络拓扑；
[0017] 根据网络拓扑中随机选择的参考节点，采用广义基尔霍夫定律建立综合能源系统RLC暂态模型。
[0018] 所述根据预先构建的RLC元件模型和能量转化元件模型，描述综合能源系统的能量耦合设备，包括：
[0019] 基于能量转化元件模型，描述不同能源之间的转化过程；
[0020] 基于流体系统的RLC元件模型、热力系统的RLC元件模型和电力系统的RLC元件模型描述能量耦合设备的耗散过程、惯性过程和弹性过程。
[0021] 所述流体系统的RLC元件模型包括不可压缩流体的RLC元件模型和可压缩流体的RLC元件模型；
[0022] 所述不可压缩流体的RLC元件模型按下式构建：
[0023]
[0024]
[0025]
[0026] 式中，ρV为不可压缩流体的
密度，lV为不可压缩流体所经过的管道长度，μV为不可压缩流体运动粘性系数，AV为不可压缩流体所经过的管道横截面积，HV为不可压缩流体的广延量，XV为不可压缩流体的强度量，RV为不可压缩流体流动过程的阻性元件，LV为不可压缩流体流动过程的感性元件；
[0027] 所述可压缩流体的RLC元件模型按下式构建：
[0028]
[0029]
[0030]
[0031]
[0032] 式中，lM为可压缩流体所经过的管道长度，kM为可压缩流体所经过的管道沿程阻力系数,aM为扰动波传递的速度，AM是可压缩流体所经过的管道横截面积，HM为可压缩流体的广延量，XM为可压缩流体的强度量，RM为可压缩流体流动过程的阻性元件，LM为可压缩流体流动过程的感性元件，CM为可压缩流体流动过程的容性元件。
[0033] 所述热力系统的RLC元件模型包括导热
传热过程RLC元件模型和
对流传热过程RLC元件模型；
[0034] 所述导热传热过程RLC元件模型按下式构建：
[0035]
[0036]
[0037] λe＝KeXe
[0038] 式中，Ae为导热传热过程中流体流经管道的横截面积，Ke为导热传热过程中流体流经管道的广延量的传递系数，λe为导热传热过程中等效能量传递系数；r1为导热传热过程中流体流经管道内径，r2为导热传热过程中流体流经管道外径，Xe为导热传热过程中流体的强度量， 为导热传热过程中流体流经管道内径处的强度量， 为导热传热过程中流体流经管道外径处的强度量，He为导热传热过程中流体的广延量，Re为导热传热过程的阻性元件，dre为导热传热过程中流体流经管道管径积分变量；
[0039] 所述对流传热过程RLC元件模型按下式构建：
[0040]
[0041] CT＝ρTrTcT
[0042] 式中，cT为对流传热过程中流体
比热容，HT为对流传热过程中流体的广延量，XT,i为对流传热流入控制体的强度量，XT,o为对流传热过程流出控制体的强度量，XT,a为环境的强度量，rT为对流传热过程中流体流经管道半径，ρT为对流传热过程中流体密度，CT为对流传热过程的容性元件，RT为对流传热过程的阻性元件。
[0043] 所述能量转化元件模型包括广义
变压器模型和广义回转器模型；
[0044] 所述广义变压器模型按下式构建：
[0045]
[0046] 式中，ΔX2为广义变压器出口端的强度量差，ΔX1为广义变压器入口端的强度量差，H2为广义变压器出口端的广延量，H1为广义变压器入口端的广延量，r为广义变压器的变换比；
[0047] 所述广义回转器模型按下式构建：
[0048]
[0049] 式中，ΔX4为广义回转器出口端的强度量差，ΔX3为广义回转器入口端的强度量差，H4为广义回转器出口端的广延量，H3为广义回转器入口端的广延量，k为广义回转器的变换比。
[0050] 所述支路的能流波动按下式确定：
[0051] ΔE＝XEHxE-XAHxA
[0052] 式中，ΔE为支路的能流波动，XE为管道末端截面的强度量，XA为管道首端截面的强度量，HxE为流经管道末端截面广延量的流率，HxA为流经管道首端截面广延量的流率。
[0053] 另一方面，本发明还提供一种综合能源系统RLC暂态模型的能流仿真装置，包括：
[0054] 动态仿真初值确定模
块，用于根据预先构建的综合能源系统RLC暂态模型，确定动态仿真初值；
[0055] 仿真模块，用于根据所述动态仿真初值，对综合能源系统RLC暂态模型进行仿真得到综合能源系统中各支路首末端参数；
[0056] 能流波动确定模块，用于通过所述支路的首末端参数，确定支路的能流波动；
[0057] 所述综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建。
[0058] 所述综合能源系统包括流体系统、热力系统和电力系统中的任意两种至全部，所述RLC元件模型包括流体系统的RLC元件模型、热力系统的RLC元件模型和电力系统的RLC元件模型中的任意两种至全部；
[0059] 所述装置还包括建模模块，所述建模模块包括：
[0060] 描述单元，用于根据RLC元件模型和能量转化元件模型，描述综合能源系统的能量耦合设备；
[0061] 网络拓扑确定单元，用于根据综合能源系统的能流路线连接各能量耦合设备，构成综合能源系统的网络拓扑；
[0062] 建模单元，用于根据网络拓扑中随机选择的参考节点，采用广义基尔霍夫定律建立综合能源系统RLC暂态模型。
[0063] 所述描述单元具体用于：
[0064] 基于能量转化元件模型，描述不同能源之间的转化过程；
[0065] 基于流体系统的RLC元件模型、热力系统的RLC元件模型和电力系统的RLC元件模型描述能量耦合设备的耗散过程、惯性过程和弹性过程。
[0066] 所述流体系统的RLC元件模型包括不可压缩流体的RLC元件模型和可压缩流体的RLC元件模型；
[0067] 所述建模单元包括第一建模单元和第二建模单元；
[0068] 所述第一建模单元按下式构建不可压缩流体的RLC元件模型：
[0069]
[0070]
[0071]
[0072] 式中，ρV为不可压缩流体的密度，lV为不可压缩流体所经过的管道长度，μV为不可压缩流体运动粘性系数，AV为不可压缩流体所经过的管道横截面积，HV为不可压缩流体的广延量，XV为不可压缩流体的强度量，RV为不可压缩流体流动过程的阻性元件，LV为不可压缩流体流动过程的感性元件；
[0073] 所述第二建模单元按下式构建可压缩流体的RLC元件模型：
[0074]
[0075]
[0076]
[0077]
[0078] 式中，lM为可压缩流体所经过的管道长度，kM为可压缩流体所经过的管道沿程阻力系数,aM为扰动波传递的速度，AM是可压缩流体所经过的管道横截面积，HM为可压缩流体的广延量，XM为可压缩流体的强度量，RM为可压缩流体流动过程的阻性元件，LM为可压缩流体流动过程的感性元件，CM为可压缩流体流动过程的容性元件。
[0079] 所述建模单元包括第三建模单元和第四建模单元；
[0080] 所述第三建模单元按下式构建导热传热过程RLC元件模型：
[0081]
[0082]
[0083] λe＝KeXe
[0084] 式中，Ae为导热传热过程中流体流经管道的横截面积，Ke为导热传热过程中流体流经管道的广延量的传递系数，λe为导热传热过程中等效能量传递系数；r1为导热传热过程中流体流经管道内径，r2为导热传热过程中流体流经管道外径，Xe为导热传热过程中流体的强度量， 为导热传热过程中流体流经管道内径处的强度量， 为导热传热过程中流体流经管道外径处的强度量，He为导热传热过程中流体的广延量，Re为导热传热过程的阻性元件，dre为导热传热过程中流体流经管道管径积分变量；
[0085] 所述第四建模单元按下式构建对流传热过程RLC元件模型：
[0086]
[0087] CT＝ρTrTcT
[0088] 式中，cT为对流传热过程中流体比
热容，HT为对流传热过程中流体的广延量，XT,i为对流传热流入控制体的强度量，XT,o为对流传热过程流出控制体的强度量，XT,a为环境的强度量，rT为对流传热过程中流体流经管道半径，ρT为对流传热过程中流体密度，CT为对流传热过程的容性元件，RT为对流传热过程的阻性元件。
[0089] 所述建模单元包括第五建模单元和第六建模单元；
[0090] 所述第五建模单元按下式构建广义变压器模型：
[0091]
[0092] 式中，ΔX2为广义变压器出口端的强度量差，ΔX1为广义变压器入口端的强度量差，H2为广义变压器出口端的广延量，H1为广义变压器入口端的广延量，r为广义变压器的变换比；
[0093] 所述第六建模单元按下式构建广义回转器模型：
[0094]
[0095] 式中，ΔX4为广义回转器出口端的强度量差，ΔX3为广义回转器入口端的强度量差，H4为广义回转器出口端的广延量，H3为广义回转器入口端的广延量，k为广义回转器的变换比。
[0096] 所述能流波动确定模块按下式确定支路的能流波动：
[0097] ΔE＝XEHxE-XAHxA
[0098] 式中，ΔE为支路的能流波动，XE为管道末端截面的强度量，XA为管道首端截面的强度量，HxE为流经管道末端截面广延量的流率，HxA为流经管道首端截面广延量的流率。
[0099] 与最接近的现有技术相比，本发明提供的技术方案具有以下有益效果：
[0100] 本发明提供的综合能源系统RLC暂态模型的能流仿真方法中，根据预先构建的综合能源系统RLC暂态模型，确定动态仿真初值；根据动态仿真初值，对综合能源系统RLC暂态模型进行仿真得到综合能源系统中各支路首末端参数；通过支路的首末端参数，确定支路的能流波动；综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建，降低了建模仿真难度，提高了可信度和计算精度，且收敛速率快；
[0101] 本发明提供的综合能源系统RLC暂态模型的能流仿真装置中，动态仿真初值确定模块，用于根据预先构建的综合能源系统RLC暂态模型，确定动态仿真初值；仿真模块，用于根据动态仿真初值，对综合能源系统RLC暂态模型进行仿真得到综合能源系统中各支路首末端参数；能流波动确定模块，用于通过支路的首末端参数，确定支路的能流波动；综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建，降低了建模仿真难度，提高了可信度和计算精度，且收敛速率快；
[0102] 本发明提供的综合能源系统RLC暂态模型完成了异质能流之间的统一建模，解决了综合能源系统动态仿真的多时间尺度问题；
[0103] 本发明可以采用成熟的龙格-库塔算法求解
软件求解综合能源系统RLC暂态模型，降低综合能源系统的仿真计算难度。
附图说明
[0104] 图1是本发明
实施例1中综合能源系统RLC暂态模型的能流仿真方法
流程图；
[0105] 图2是本发明实施例1中感应
电动机的能量网络等值模型示意图；
[0106] 图3是本发明实施例1中
离心泵的能量网络模型示意图；
[0107] 图4是本发明实施例1中电力-流体耦合综合能源系统示意图；
[0108] 图5是本发明实施例1中综合能源能量系统的网络拓扑结构图；
[0109] 图6是本发明实施例1中
离心泵出口处体积流率动态仿真结果示意图。
具体实施方式
[0110] 下面结合附图对本发明作进一步详细说明。
[0111] 实施例1
[0112] 本发明实施例1提供了一种综合能源系统RLC暂态模型的能流仿真方法，具体流程图如图1所示，具体过程如下：
[0113] S101：根据预先构建的综合能源系统RLC暂态模型，确定动态仿真初值；
[0114] S102：根据动态仿真初值，采用龙格-库塔算法对综合能源系统RLC暂态模型进行仿真得到综合能源系统中各支路首末端参数；
[0115] S103：通过所述支路的首末端参数，确定支路的能流波动；
[0116] 上述综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建。
[0117] 流体系统和热力系统中，支路为管道，电力系统中，支路为线路。
[0118] 综合能源系统包括流体系统、热力系统和电力系统中的任意两种至全部，RLC元件模型包括流体系统的RLC元件模型、热力系统的RLC元件模型和电力系统的RLC元件模型中的任意两种至全部；
[0119] 上述的电力系统的RLC元件模型基于RLC暂态
电路构建。
[0120] 综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建，包括：
[0121] 根据RLC元件模型和能量转化元件模型，描述综合能源系统的能量耦合设备；
[0122] 根据综合能源系统的能流路线连接各能量耦合设备，构成综合能源系统的网络拓扑；
[0123] 根据网络拓扑中随机选择的参考节点，采用广义基尔霍夫定律建立综合能源系统RLC暂态模型。
[0124] 根据预先构建的RLC元件模型和能量转化元件模型，描述综合能源系统的能量耦合设备，包括：
[0125] 基于能量转化元件模型，描述不同能源之间的转化过程；
[0126] 基于流体系统的RLC元件模型、热力系统的RLC元件模型和电力系统的RLC元件模型描述能量耦合设备的耗散过程、惯性过程和弹性过程。
[0127] 上述流体系统的RLC元件模型包括不可压缩流体的RLC元件模型和可压缩流体的RLC元件模型；
[0128] 由于阻性元件描述耗散过程，感性元件描述惯性过程，容性元件描述弹性过程，所以本发明实施例1中的不可压缩流体的RLC元件模型用于描述能量耦合设备的耗散过程和惯性过程，可压缩流体的RLC元件模型用于描述能量耦合设备的耗散过程、惯性过程和弹性过程。
[0129] 针对不可压缩流体，以体积流率为广延量，压强为强度量，根据Navier-Stokes方程和能量网络理论，得到不可压缩流体的RLC元件模型，不可压缩流体的RLC元件模型按下式构建：
[0130]
[0131]
[0132]
[0133] 式中，ρV为不可压缩流体的密度，lV为不可压缩流体所经过的管道长度，μV为不可压缩流体运动粘性系数，AV为不可压缩流体所经过的管道横截面积，HV为不可压缩流体的广延量，XV为不可压缩流体的强度量，RV为不可压缩流体流动过程的阻性元件，LV为不可压缩流体流动过程的感性元件；
[0134] 针对可压缩流体，以
质量流率为广延量，压强为强度量，根据连续性公理，得可压缩流体的RLC元件模型，可压缩流体的RLC元件模型按下式构建：
[0135]
[0136]
[0137]
[0138]
[0139] 式中，lM为可压缩流体所经过的管道长度，kM为可压缩流体所经过的管道沿程阻力系数,aM为扰动波传递的速度，AM是可压缩流体所经过的管道横截面积，HM为可压缩流体的广延量，XM为可压缩流体的强度量，RM为可压缩流体流动过程的阻性元件，LM为可压缩流体流动过程的感性元件，CM为可压缩流体流动过程的容性元件。
[0140] 综合能源系统考虑传热时，热力系统只会单向地趋于平衡态，不会发生振荡运动，故传热学中没有热感这一概念，因此热力系统建模不考虑感性元件模型。上述热力系统的RLC元件模型包括导热传热过程RLC元件模型和对流传热过程RLC元件模型；由于阻性元件描述耗散过程，容性元件描述弹性过程，所以本发明实施例1中的导热传热过程RLC元件模型用于描述能量耦合设备的耗散过程，对流传热过程RLC元件模型用于描述能量耦合设备的耗散过程和弹性过程。
[0141] 导热传热过程，以
温度为强度量，以热流率为广延量，根据传递公理，得导热传热过程RLC元件模型，导热传热过程RLC元件模型按下式构建：
[0142]
[0143]
[0144] λe＝KeXe
[0145] 式中，Ae为导热传热过程中流体流经管道的横截面积，Ke为导热传热过程中流体流经管道的广延量的传递系数，λe为导热传热过程中等效能量传递系数；r1为导热传热过程中流体流经管道内径，r2为导热传热过程中流体流经管道外径，Xe为导热传热过程中流体的强度量， 为导热传热过程中流体流经管道内径处的强度量， 为导热传热过程中流体流经管道外径处的强度量，He为导热传热过程中流体的广延量，Re为导热传热过程的阻性元件，dre为导热传热过程中流体流经管道管径积分变量；
[0146] 对流传热过程，以温度为强度量，以质量流率为广延量，得对流传热过程RLC元件模型，对流传热过程RLC元件模型按下式构建：
[0147]
[0148] CT＝ρTrTcT
[0149] 式中，cT为对流传热过程中流体比热容，HT为对流传热过程中流体的广延量，XT,i为对流传热流入控制体的强度量，XT,o为对流传热过程流出控制体的强度量，XT,a为环境的强度量，rT为对流传热过程中流体流经管道半径，ρT为对流传热过程中流体密度，CT为对流传热过程的容性元件，RT为对流传热过程的阻性元件。
[0150] 上述能量转化元件模型包括广义变压器模型和广义回转器模型；
[0151] 本领域技术人员公知，变压器特指电力系统的变压器，本发明实施例1中的广义变压器为适用于电力系统、热力系统和流体系统的变压器，同样，本发明实施例1中的广义回转器也为适用于电力系统、热力系统和流体系统的回转器。
[0152] 广义变压器是指：能量流过时，功率不变，广义变压器入口端的广延量和强度量差与出口端的广延量和强度量差均以预设的变换比变化；上述广义变压器模型按下式构建：
[0153]
[0154] 式中，ΔX2为广义变压器出口端的强度量差，ΔX1为广义变压器入口端的强度量差，H2为广义变压器出口端的广延量，H1为广义变压器入口端的广延量，r为广义变压器的变换比；
[0155] 广义回转器包括：能量流过时，广义回转器入口端的广延量和出口端的的强度量差以预设的变换比变化，且广义回转器出口端的广延量和入口端的强度量差以预设的变换比变化；上述广义回转器模型按下式构建：
[0156]
[0157] 式中，ΔX4为广义回转器出口端的强度量差，ΔX3为广义回转器入口端的强度量差，H4为广义回转器出口端的广延量，H3为广义回转器入口端的广延量，k为广义回转器的变换比。
[0158] 可以利用上述过程建立包括感应电动机的能量网络等值模型和离心泵的能量网络模型的能量耦合设备模型。以阻性元件代表电动机转速差、
轴承摩擦及
转子摩擦特性，以感性元件代表转子及轴的
转动惯量特性，以容性元件代表机轴的柔度特性，根据能量传输的顺序连接基本元件，构建感应电动机的能量网络等值模型，如图2所示，其中Rmw是阻性元件，代表电动机转速差；Rmf是阻性元件，代表电动机轴承摩擦及转子摩擦特性；Lmx是能量惯性的感性元件，代表转子及轴的转动惯量特性；Cmk是储存能量的容性元件，代表机轴的柔度等特性；GY为广义回转器元件，表示机电转换关系；km为感应
电机内广义回转器的转矩转换系数；Hi为感应电机
定子电流；Hω为感应电机的转子
角速度；ΔXE为感应电机定子两端的
电压差；ΔXT为感应电机转子转矩。
[0159] 以阻性元件代表应离心泵转子摩擦特性、泄露特性与
阀门出口压力损失特性，以感性元件对应离心泵转子转动惯量特性，根据能量传输的顺序连接基本元件，构建离心泵的能量网络模型，如图3所示，图3中，Rpf为阻性元件，代表离心泵转子摩擦特性；Rpξ为阻性元件，代表离心泵泄露特性；Rpr为阻性元件，代表离心泵阀门出口压力损失特性；Lpx为储存能量的感性元件，代表离心泵转子转动惯量特性；ΔXT为离心泵的输入转矩，Hω为离心泵的输入角速度，XP为离心泵的输出压力；HqV为离心泵的输出体积流量，rp为
叶轮回旋器模数，即离心泵内广义回转器GY的回转系数。
[0160] 某电力-流体耦合综合能源系统如图4所示，实线代表电力网络，传输
电能；虚线代表流体网络，传输
热能；电力网络与流体网络依靠离心泵TM(由感应电动机提供动力)和热负荷(PT＝PTh+PTe)耦合在一起。PTh代表用户侧换热器提供的热量；PTe代表电热器提供的热量；Se为电源；ST为热源，输入流体的温度保持不变；n1、n2为电力网络各节点编号；n3～n6为流体网络各节点编号；b1、b2为电力网络各支路编号；b3～b7为流体网络各支路编号；①、②为电力网络各回路编号；③为流体网络的回路编号。
[0161] 基于图4，确定如图5所示的综合能源能量系统的网络拓扑，图4中，n1、n2为电力网络各节点编号；n3～n6为流体网络各节点编号；n7～n13为由离心泵与感应电机组成的能量耦合设备内的网络节点编号。b1、b2为电力网络各支路编号；b3～b7为流体网络各支路编号；b8～b16为由离心泵与感应电机组成的能量耦合设备的各元件支路编号。①、②为电力网络各回路编号；③为流体网络的回路编号；④、⑤为能量耦合设备内感应电机的网络回路编号；⑥、⑦则为能量耦合设备内离心泵的网络回路编号。
[0162] 上述S103中，支路的能流波动按下式确定：
[0163] ΔE＝XEHxE-XAHxA
[0164] 式中，ΔE为支路的能流波动，XE为管道末端截面的强度量，XA为管道首端截面的强度量，HxE为流经管道末端截面广延量的流率，HxA为流经管道首端截面广延量的流率。
[0165] 为了验证本发明实施例1的有效性，基于图4，对三种不同扰动发生后的能量网络响应进行时域仿真，该三种扰动分别为：(1)
电网侧电压源Se突变，其电压有效值由220V下降至218V；(2)增大热网侧离心泵的出口阀门开度，阀门阻力系数突变为Rpr＝10；(3)减小热网侧离心泵的出口阀门开度，阀门阻力系数突变为Rpr＝20。
[0166] 以流体网络出口处的流量作为仿真结果，得到如图6的离心泵出口处体积流率动态仿真结果，从图6可以看出，当扰动1发生时，电压突变首先冲击感应电动机，由于感应电机自身的惯性作用，振荡将持续十多秒钟，电动机转速与转矩的大幅度下降导致离心泵出口的体积流量H4随着震荡逐步减小；当扰动2发生时，离心泵出口阀门开度突然增大，阀门对离心泵产生的阻力减小，离心泵体积流量H4增大，但由于感应电机以及离心泵的惯性作用，存在震荡现象；当扰动3发生时，离心泵出口阀门开度突然减小，与扰动2相反，阀门对离心泵产生的阻力增加，离心泵体积流量H4减小，并由于惯性作用发生十数秒的震荡。
[0167] 从仿真结果来看，采用综合能源系统暂态能流仿真方法所实现的综合能源系统动态仿真与实际工况十分接近，证明了能量网络理论体系的简易性以及在一定精度范围内的可行性。
[0168] 实施例2
[0169] 基于同一发明构思，本发明实施例2还提供一种综合能源系统RLC暂态模型的能流仿真装置，包括动态仿真初值确定模块、仿真模块和能流波动确定模块，下面对上述几个模块的功能进行详细说明：
[0170] 动态仿真初值确定模块，用于根据预先构建的综合能源系统RLC暂态模型，确定动态仿真初值；
[0171] 仿真模块，用于根据动态仿真初值，对综合能源系统RLC暂态模型进行仿真得到综合能源系统中各支路首末端参数；
[0172] 能流波动确定模块，用于通过所述支路的首末端参数，确定支路的能流波动；
[0173] 上述综合能源系统RLC暂态模型基于RLC元件模型和能量转化元件模型构建。
[0174] 综合能源系统包括流体系统、热力系统和电力系统中的任意两种至全部，所述RLC元件模型包括流体系统的RLC元件模型、热力系统的RLC元件模型和电力系统的RLC元件模型中的任意两种至全部；
[0175] 本发明实施例2提供的装置还包括建模模块，建模模块包括：
[0176] 描述单元，用于根据RLC元件模型和能量转化元件模型，描述综合能源系统的能量耦合设备；
[0177] 网络拓扑确定单元，用于根据综合能源系统的能流路线连接各能量耦合设备，构成综合能源系统的网络拓扑；
[0178] 建模单元，用于根据网络拓扑中随机选择的参考节点，采用广义基尔霍夫定律建立综合能源系统RLC暂态模型。
[0179] 描述单元具体用于：
[0180] 基于能量转化元件模型，描述不同能源之间的转化过程；
[0181] 基于流体系统的RLC元件模型、热力系统的RLC元件模型和电力系统的RLC元件模型描述能量耦合设备的耗散过程、惯性过程和弹性过程。
[0182] 上述流体系统的RLC元件模型包括不可压缩流体的RLC元件模型和可压缩流体的RLC元件模型；
[0183] 建模单元包括第一建模单元和第二建模单元；
[0184] 上述第一建模单元按下式构建不可压缩流体的RLC元件模型：
[0185]
[0186]
[0187]
[0188] 式中，ρV为不可压缩流体的密度，lV为不可压缩流体所经过的管道长度，μV为不可压缩流体运动粘性系数，AV为不可压缩流体所经过的管道横截面积，HV为不可压缩流体的广延量，XV为不可压缩流体的强度量，RV为不可压缩流体流动过程的阻性元件，LV为不可压缩流体流动过程的感性元件；
[0189] 上述第二建模单元按下式构建可压缩流体的RLC元件模型：
[0190]
[0191]
[0192]
[0193]
[0194] 式中，lM为可压缩流体所经过的管道长度，kM为可压缩流体所经过的管道沿程阻力系数,aM为扰动波传递的速度，AM是可压缩流体所经过的管道横截面积，HM为可压缩流体的广延量，XM为可压缩流体的强度量，RM为可压缩流体流动过程的阻性元件，LM为可压缩流体流动过程的感性元件，CM为可压缩流体流动过程的容性元件。
[0195] 建模单元包括第三建模单元和第四建模单元；
[0196] 上述第三建模单元按下式构建导热传热过程RLC元件模型：
[0197]
[0198]
[0199] λe＝KeXe
[0200] 式中，Ae为导热传热过程中流体流经管道的横截面积，Ke为导热传热过程中流体流经管道的广延量的传递系数，λe为导热传热过程中等效能量传递系数；r1为导热传热过程中流体流经管道内径，r2为导热传热过程中流体流经管道外径，Xe为导热传热过程中流体的强度量， 为导热传热过程中流体流经管道内径处的强度量， 为导热传热过程中流体流经管道外径处的强度量，He为导热传热过程中流体的广延量，Re为导热传热过程的阻性元件，dre为导热传热过程中流体流经管道管径积分变量；
[0201] 上述第四建模单元按下式构建对流传热过程RLC元件模型：
[0202]
[0203] CT＝ρTrTcT
[0204] 式中，cT为对流传热过程中流体比热容，HT为对流传热过程中流体的广延量，XT,i为对流传热流入控制体的强度量，XT,o为对流传热过程流出控制体的强度量，XT,a为环境的强度量，rT为对流传热过程中流体流经管道半径，ρT为对流传热过程中流体密度，CT为对流传热过程的容性元件，RT为对流传热过程的阻性元件。
[0205] 建模单元包括第五建模单元和第六建模单元；
[0206] 上述第五建模单元按下式构建广义变压器模型：
[0207]
[0208] 式中，ΔX2为广义变压器出口端的强度量差，ΔX1为广义变压器入口端的强度量差，H2为广义变压器出口端的广延量，H1为广义变压器入口端的广延量，r为广义变压器的变换比；
[0209] 上述第六建模单元按下式构建广义回转器模型：
[0210]
[0211] 式中，ΔX4为广义回转器出口端的强度量差，ΔX3为广义回转器入口端的强度量差，H4为广义回转器出口端的广延量，H3为广义回转器入口端的广延量，k为广义回转器的变换比。
[0212] 能流波动确定模块按下式确定支路的能流波动：
[0213] ΔE＝XEHxE-XAHxA
[0214] 式中，ΔE为支路的能流波动，XE为管道末端截面的强度量，XA为管道首端截面的强度量，HxE为流经管道末端截面广延量的流率，HxA为流经管道首端截面广延量的流率。
[0215] 为了描述的方便，以上所述装置的各部分以功能分为各种模块或单元分别描述。当然，在实施本
申请时可以把各模块或单元的功能在同一个或多个软件或
硬件中实现。
[0216] 本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或
计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘
存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
[0217] 本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方
框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程
数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中
指定的功能的装置。
[0218] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0219] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0220] 最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行
修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的
权利要求保护范围之内。
