CN112329217B - 一种电-水综合能源系统的准稳态能流分析方法 - Google Patents
一种电-水综合能源系统的准稳态能流分析方法 Download PDFInfo
- Publication number
- CN112329217B CN112329217B CN202011148025.8A CN202011148025A CN112329217B CN 112329217 B CN112329217 B CN 112329217B CN 202011148025 A CN202011148025 A CN 202011148025A CN 112329217 B CN112329217 B CN 112329217B
- Authority
- CN
- China
- Prior art keywords
- water
- pipeline
- node
- flow
- water pump
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 288
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000005206 flow analysis Methods 0.000 title claims abstract description 14
- 238000009826 distribution Methods 0.000 claims abstract description 102
- 230000001052 transient effect Effects 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 claims abstract description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 44
- 230000011218 segmentation Effects 0.000 claims description 28
- 238000004088 simulation Methods 0.000 claims description 17
- 230000001133 acceleration Effects 0.000 claims description 7
- -1 hydrogen Chemical class 0.000 claims description 5
- 239000001257 hydrogen Substances 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 230000018109 developmental process Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003809 water extraction Methods 0.000 description 1
- 238000013014 water-saving technology Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Pipeline Systems (AREA)
Abstract
本发明公开一种电‑水综合能源系统的准稳态能流分析方法，步骤包括：1)获取电‑水综合能源系统的网络结构参数；2)获取电‑水综合能源系统的初始稳态能流，即时间t＝0时配水网的能流分布和配电网的潮流分布。3)利用特征线法计算配水网在下一时段的暂态能流和水泵功率；4)将下一时间步长的水泵功率代入配电网中，作为电负荷，并计算下一时间步长的配电网能流；5)保存能流计算结果，若当前时间t小于总时长T，则返回步骤3)，否则，输出配电网能流。本发明在分析电‑水综合能源系统能流过程中综合考虑了配水网的暂态特性，能够分析时间尺度为分钟级或秒级的水力系统暂态过程对电‑水综合能源系统的影响。
Description
技术领域
本发明涉及综合能源技术领域，具体是一种电-水综合能源系统的准稳态能流分析方法。
背景技术
世界范围内日益严峻的能源与环境危机催生了能源互联网、综合能源系统(Integrated Energy System，IES)等新型能源系统的快速发展。综合能源系统作为能源互联网的物理载体，以打破多种异质能源系统的行业壁垒，实现多种能源互联互济、综合利用为基本特征，成为未来能源系统的重要发展方向。
虽然电力系统和输配水系统之间长期处于独立规划、分立运行的模式，但是这两个系统密切关联、相互依存，主要表现在电力生产过程需要大量取水，而水的提取、处理和输送又需要消耗大量电力。此前国内外已有大量研究关注水与能源之间的依存关系(一般称为water-energy nexus)和协调发展问题，这些研究大多从能源与环境的宏观视角研究某个国家/地区水资源和能源的供需关系、能源开发利用的水资源约束、能源生产的水足迹及节水技术，以及社会水循环过程中的能耗和节能问题。近2年，随着对综合能源系统研究的不断发展和深化，国内外学者开始立足于电力系统与输配水系统的相互影响来研究水-能关系，提出将输配水系统纳入综合能源系统的框架，用集成、融合的新思路来构建和发展电-水联合系统 (Integrated Electricity-Water System，电-水综合能源系统)。
然而，现有面向电-水综合能源系统的研究都是基于对电力系统和输配水系统的稳态分析，认为管道中的水流变化十分缓慢，因此一个点的水力状态不会随时间变化，研究的时间尺度多为小时级。上述模型难以描述和分析时间尺度为分钟级或秒级的水力系统暂态过程对电-水综合能源系统的影响。
发明内容
本发明的目的是提供一种电-水综合能源系统的准稳态能流分析方法，包括以下步骤：
1)获取电-水综合能源系统的网络结构参数，并设置管道支路分段点数为I，管道分段步长Δx，时间步长为Δt，仿真时长为T。所述分段节点包括上、下游端节点。
管道分段步长Δx和时间步长为Δt满足如下关系：
式中，c为水锤波速。
2)获取电-水综合能源系统的初始稳态能流，即时间t＝0时配水网的能流分布和配电网的潮流分布。
3)利用特征线法计算配水网在下一时段的暂态能流和水泵功率，步骤包括：
3.1)基于已知的配水网的能流分布计算管道内下一时刻各分段点位置的管道流量q(t+1)和管道水头H(t+1)，步骤包括：
3.1.1)利用特征线法计算配水网暂态能流的主方程，即：
式中，q表示管道流量。H表示管道水头。g为重力加速度。c 为水锤波速。D为管道直径。A为管道断面积。f为水力阻力系数。x为空间。t为时间。C+和C-分别表示沿管道两个不同方向的正特性方程和负特性方程。
3.1.2)将公式(2)和公式(3)改写如下：
3.1.3)联立公式(4)和公式(5)，计算管道在下一时刻各分段点位置的管道流量和管道水头，得到：
式中，(t+1)表示下一时间步长。下标i表示支路的分段节点。i＝ 1、2、3，...，I。(i-1)和(i+1)表示管道分段节点i的相邻节点。qi(t+1)表示下一时间步长节点i的管道流量；qi-1(t)、qi+1(t)表示时间t 节点i-1和节点i+1的管道流量；Hi-1(t)、Hi+1(t)表示时间t节点i-1 和节点i+1的管道水头；Hi(t+1)表示下一时间步长节点i的管道水头。
3.2)基于与支路两端分段点相连的元件确定相应端节点的连续性条件。
节点j的连续性条件分别如公式(8)和公式(9)所示，即：
式中，下标ta、ws、d、u和p分别表示蓄水池、水源、水负荷、水泵和管道。Nta(j)、Nws(j)、Nd(j)、Nu(j)和Np(j)分别表示蓄水池、水源、水负荷、水泵和管道的节点集合。
分别表示t时刻蓄水池节点、水源节点、水负荷节点、水泵上游节点、水泵下游节点、管道上游端节点、管道下游端节点的流量。
分别表示t时刻蓄水池节点、水源节点、水负荷节点、管道上游端节点、管道下游端节点的水头。
分别表示t时刻水泵出水和水泵进水。1表示上游端节点，I表示下游端节点。
3.3)根据管道两端连接的不同设备确定边界条件，步骤包括：
3.3.1)确定管道上游端和下游端的边界条件。
其中，管道上游端边界条件如下所示：
式中，qjk,1(t+1)和Hjk,1(t+1)分别表示支路j-k下一时间步长上游端节点的未知管道流量和管道水头。
其中，变量Cajk和变量Cnjk分别如下所示：
式中，qjk,2和Hjk,2分别表示支路j-k上游端节点的相邻分段节点的已知管道流量和管道水头。g为重力加速度。c为水锤波速。D为管道直径。A为管道断面积。
管道下游端边界条件如下所示：
式中，qjk,I(t+1)和Hjk,I(t+1)分别表示支路j-k下一时间步长下游端节点的未知管道流量和管道水头。
其中，变量Cpjk如下所示：
式中，qjk,I-1和Hjk,I-1分别表示支路j-k下游端节点的相邻分段节点的已知管道流量和管道水头。
3.3.2)确定水源边界条件，即：
式中，Hj(t)和qj(t)分别表示t时刻水源节点j的水头和流量。Hws(t)、 qws(t)分别表示t时刻水源的管道水头和管道流量。
3.3.3)确定蓄水池的边界条件，步骤如下：
3.3.3.1)获取下一个时间步长与蓄水池相连支路的相应端的分段点的水头。
3.3.3.2)建立蓄水池的边界条件，即：
式中，Hj(t+1)为下一个时间步长与蓄水池节点j相连支路的相应端的分段点的水头。An为蓄水池相连支路的连接处的断面积。Co为蓄水池的孔口流量系数。qta(t)为t时刻蓄水池的管道流量；
基于蓄水池的边界条件，计算下一时间步长的蓄水池水头，即：
式中，Ata为蓄水池的底面积。
3.3.4)确定变速离心水泵边界条件，即：
式中，ΔHu,jk(t)为t时刻水泵j-k扬程，其值等于水泵的下游端节点与上游端节点的水头差。qu,jk(t)为t时刻水泵流量，其值等于水泵的上游端节点和下游端节点的流量。H0jk、cjk和mjk分别为水泵j-k 静扬程、内阻系数和扬程指数。ωjk(t)为t时刻水泵j-k的相对转速，即工作转速与额定转速之比。
3.4)将节点连续性方程、各管道中已知分段点的流量q(t)和水头H(t)代入边界条件中，求解得到下一时间步长各端节点的流量和水头。
3.5)根据步骤3.1)至步骤3.4)计算出的配水网能流分布，利用水泵功率计算公式确定下一时间步长的水泵功率。
t时刻的水泵j-k功率Pu,jk如下所示：
式中，ρ为水的密度。ηu,jk为水泵j-k的效率。
4)将下一时间步长的水泵功率代入配电网中，作为电负荷，并计算下一时间步长的配电网能流。
t时刻的配电网能流如下所示：
式中，
ψ∈{a,b,c}。
和
分别为t时刻节点j的
相注入有功和无功功率。
和
分别为t时刻节点j的
相电压模值和相角。
和
分别为配电网节点导纳矩阵对应节点j、节点k、
相、ψ相的电导和电纳。Ne为配电网节点数。
分别表示t时刻节点j的
相角和节点k的ψ相角。
和
分别为t时刻节点k的ψ相电压模值和相角；
5)保存能流计算结果，若当前时间t小于总时长T，则返回步骤3)，否则，输出配电网能流。
本发明的技术效果是毋庸置疑的，本发明在分析电-水综合能源系统能流过程中综合考虑了配水网的暂态特性，与现有电-水稳态能流模型相比能够分析时间尺度为分钟级或秒级的水力系统暂态过程对电-水综合能源系统的影响。
附图说明
图1为本发明方法的流程框图；
图2为x-t平面的特征线示意图；
图3为本发明方法所采用的电-水综合能源系统结构示意图；
图4为节点i连续性条件的示意图；
图5为管道j-k上游端边界条件的示意图；
图6为管道j-k下游端边界条件的示意图；
图7为水源节点j边界条件的示意图；
图8为蓄水池节点j边界条件的示意图；
图9为水泵j-k两节点边界条件的示意图；
图10为场景1和场景2中水泵N1-N2的流量，扬程和功率的仿真结果图；
图11为场景1和场景2中母线B8的A相电压模值的仿真结果图。
具体实施方式
下面结合实施例对本发明作进一步说明，但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下，根据本领域普通技术知识和惯用手段，做出各种替换和变更，均应包括在本发明的保护范围内。
实施例1：
参见图1至图11，一种电-水综合能源系统的准稳态能流分析方法，包括以下步骤：
1)获取电-水综合能源系统(IEWS)的网络结构参数，并设置管道支路分段点数为I，管道分段步长Δx，时间步长为Δt，仿真时长为T。所述分段节点包括上、下游端节点。
管道分段步长Δx和时间步长为Δt满足如下关系：
式中，c为水锤波速。
2)获取电-水综合能源系统的初始稳态能流，即时间t＝0时配水网的能流分布和配电网的潮流分布。
3)利用特征线法计算配水网在下一时段的暂态能流和水泵功率，步骤包括：
3.1)基于已知的配水网的能流分布计算管道内下一时刻各分段点位置的管道流量q(t+1)和管道水头H(t+1)。
3.1.1)利用特征线法计算配水网暂态能流的主方程，即：
式中，q表示管道流量。H表示管道水头。g为重力加速度。c 为水锤波速。D为管道直径。A为管道断面积。f为水力阻力系数。 x为空间。t为时间。C+和C-分别表示沿管道两个不同方向的正特性方程和负特性方程。
3.1.2)将公式(2)和公式(3)改写如下：
式中，q(x,t+Δt)表示在空间x、t至t+Δt时间段的管道流量； q(x-Δx,t)表示在空间x-Δx、时间t的管道流量；q(x+Δx,t)表示在空间 x+Δx、时间t的管道流量；H(x,t+Δt)表示在空间x、t至t+Δt时间段的管道水头；H(x-Δx,t)表示在空间x-Δx、时间t的管道水头；H(x+Δx,t) 表示在空间x+Δx、时间t的管道水头。
3.1.3)联立公式(4)和公式(5)，计算管道在下一时刻各分段点位置的管道流量和管道水头，得到：
式中，(t+1)表示下一时间步长。下标i表示支路的分段节点。i＝ 1、2、3，...，I。(i-1)和(i+1)表示管道分段节点i的相邻节点。 qi(t+1)表示下一时间步长节点i的管道流量；qi-1(t)、qi+1(t)表示时间t 节点i-1和节点i+1的管道流量；Hi-1(t)、Hi+1(t)表示时间t节点i-1 和节点i+1的管道水头；Hi(t+1)表示下一时间步长节点i的管道水头。
3.2)基于与支路两端分段点相连的元件确定相应节点的连续性条件。
节点j的连续性条件分别如公式(8)和公式(9)所示，即：
式中，下标ta、ws、d、u和p分别表示蓄水池、水源、水负荷、水泵和管道。Nta(j)、Nws(j)、Nd(j)、Nu(j)和Np(j)分别表示蓄水池、水源、水负荷、水泵和管道的节点集合。
分别表示t时刻蓄水池节点、水源节点、水负荷节点、水泵上游节点、水泵下游节点、管道上游端节点、管道下游端节点的流量。
分别表示t时刻蓄水池节点、水源节点、水负荷节点、管道上游端节点、管道下游端节点的水头。
分别表示t时刻水泵出水和水泵进水。1表示上游端节点，I表示下游端节点。
3.3)根据管道两端连接的不同设备确定边界条件，步骤包括：
3.3.1)确定管道上游端和下游端的边界条件。
其中，管道上游端边界条件如下所示：
式中，qjk,1(t+1)和Hjk,1(t+1)分别表示支路j-k下一时间步长上游端节点的未知管道流量和管道水头。
其中，变量Cajk和变量Cnjk分别如下所示：
式中，qjk,2和Hjk,2分别表示支路j-k上游端节点的相邻分段节点的已知管道流量和管道水头。g为重力加速度。c为水锤波速。D为管道直径。A为管道断面积。
管道下游端边界条件如下所示：
式中，qjk,I(t+1)和Hjk,I(t+1)分别表示支路j-k下一时间步长下游端节点的未知管道流量和管道水头。
其中，变量Cpjk如下所示：
式中，qjk,I-1和Hjk,I-1分别表示支路j-k下游端节点的相邻分段节点的已知管道流量和管道水头。
3.3.2)确定水源边界条件，即：
式中，Hj(t)和qj(t)分别表示t时刻水源节点j的水头和流量。Hws(t)、 qws(t)分别表示t时刻水源的管道水头和管道流量。
3.3.3)确定蓄水池的边界条件，步骤如下：
3.3.3.1)获取下一个时间步长与蓄水池相连支路的相应端的分段点的水头。
3.3.3.2)建立蓄水池的边界条件，即：
式中，Hj(t+1)为下一个时间步长与蓄水池节点j相连支路的相应端的分段点的水头。An为蓄水池相连支路的连接处的断面积。Co为蓄水池的孔口流量系数。qta(t)为t时刻蓄水池的管道流量；
基于蓄水池的边界条件，计算下一时间步长的蓄水池水头，即：
式中，Ata为蓄水池的底面积。
3.3.4)确定变速离心水泵边界条件，即：
式中，ΔHu,jk(t)为t时刻水泵j-k扬程，其值等于水泵的下游端节点与上游端节点的水头差。qu,jk(t)为t时刻水泵流量，其值等于水泵的上游端节点和下游端节点的流量。H0jk、cjk和mjk分别为水泵j-k 静扬程、内阻系数和扬程指数。ωjk(t)为t时刻水泵j-k的相对转速，即工作转速与额定转速之比。
3.4)将节点连续性方程、各管道中已知分段点的流量q(t)和水头H(t)代入边界条件中，求解得到下一时间步长各端节点的流量和水头。
3.5)根据步骤3.1)至步骤3.4)计算出的配水网能流分布，利用水泵功率计算公式确定下一时间步长的水泵功率。
t时刻的水泵j-k功率Pu,jk如下所示：
式中，ρ为水的密度。ηu,jk为水泵j-k的效率。
4)将下一时间步长的水泵功率代入配电网中，作为电负荷，并计算下一时间步长的配电网能流。
t时刻的配电网能流如下所示：
式中，
ψ∈{a,b,c}。
和
分别为t时刻节点j的
相注入有功和无功功率。
和
分别为t时刻节点j的
相电压模值和相角。
和
分别为配电网节点导纳矩阵对应节点j、节点k、
相、ψ相的电导和电纳。Ne为配电网节点数。
分别表示t时刻节点j的
相角和节点k的ψ相角。
和
分别为t时刻节点k的ψ相电压模值和相角；
5)保存能流计算结果，若当前时间t小于总时长T，则返回步骤3)，否则，输出配电网能流。
实施例2：
一种考虑输配水系统暂态特性的IEWS准稳态能流分析方法，步骤如下：
1)输入IEWS的网络结构参数，设置管道支路分段点数为3(包括上、下游端节点)，取水锤波速为1000m/s，管道分段步长为1000m，仿真时间步长为1s，仿真时长为120s。
为了保证解的精度，必须使Δt和Δx满足如下Courant条件的时间-空间步长关系：
2)输入t＝0时配水网的能流分布和配电网的潮流分布。
3)基于已知的上一时段的配水网能流分布，采用特征线法计算配水网在下一仿真步长的暂态能流。
3.1)基于主方程计算管道中在任意时刻各分段点位置的流量和水头。采用特征线法计算配水网暂态能流的主方程是分别沿管道两个方向的两个一阶线性微分方程，即：
式中，q表示管道流量。H表示管道水头。g为重力加速度。c 为水锤波速；D为管道直径。A为管道断面积。f为水力阻力系数，与水流雷诺数、管材粗糙度有关，一般由经验公式计算得到或根据管道材质直接查表获取。x为空间。t为时间。C+和C-是特征线方程的名称，表示方程式(2)和(3)分别为沿管道两个不同方向的正特性方程和负特性方程。若满足dx/dt＝c条件，等式(2)成立。反之则等式(3)成立。式(2)和(3)还可以写为：
如图2所示，线C+和C-分别表示式(4)和式(5)。对于时间 (t+Δt)，位于空间x处的坐标点P的流量q和水头H可以分别通过时间t，位于空间(x-Δx)和(x+Δx)上的坐标点A和B获得。
联立公式4和5，获得管道中在任意时刻各分段点位置的流量和水头的求解式如下所示：
式中，(t+1)表示下一时间步长。下标i(i＝1、2、3，...，I)表示支路的分段节点(包括上、下游端节点)。例如，下标“1”代表上游端节点。下标“I”代表下游端节点。(i-1)和(i+1)分别为管道中的上一个和下一个分段节点。
3.2)基于与支路两端分段点相连的元件确定不同节点的连续性条件，如图4所示。对于在任何时间的任意支路一端分段节点j，节点连续性条件可以表示如下：
式中，下标ta，ws，d，u和p分别表示蓄水池，水源，水负荷，水泵和管道，Nta(j)，Nws(j)，Nd(j)，Nu(j)和Np(j)分别为相应类型的节点集合,j的不同下标1和I分别表示j为上游端节点和下游端节点。
3.3)根据管道两端连接和设备的不同确定边界条件。为了简化推导,对于支路j-k，定义变量Cajk,Cnjk和Cpjk的等价式如下：
3.3.1)确定管道上游端边界条件，如图5所示为管道j-k上游端边界条件示意图，水流在管道中从节点j流向k，其边界条件方程为：
3.3.2)确定管道下游端边界条件，如图6所示为管道j-k下游端边界条件示意图，水流在管道中从节点j流向k，其边界条件方程为：
3.3.3)确定水源边界条件，如图7所示，水源节点j的水头和流量符合以下关系式
3.3.4)如图8所示，Hta代表蓄水池水头。为了确定蓄水池的边界条件，需要首先获取下一个时间步长的与蓄水池相连支路的相应端的分段点的水头Hj(t+1)(Hta(t)和Hj(t+1)不相等)。当[Hta(t)-Hj(t+1)] 为正,蓄水池的边界条件如下:
式中，An为蓄水池相连支路的连接处的断面积；Co为蓄水池的孔口流量系数。
根据蓄水池的边界条件，求出下一时间步长的蓄水池水头，即：
式中，Ata为蓄水池的底面积。
当[Hta(t)-Hj(t+1)]为负,则蓄水池的边界条件如下:
根据蓄水池的边界条件，求出下一时间步长的蓄水池水头，即：
3.3.5)如图9所示，变速离心水泵j-k形成了两节点的边界条件。根据变速离心水泵j-k的水泵扬程计算公式计算水泵扬程，即：
式中，ΔHu,jk(t)为t时刻水泵扬程，其值等于水泵的下游端节点与上游端节点的水头差。qu,jk(t)为t时刻水泵流量，其值等于水泵的上游端节点和下游端节点的流量。H0jk、cjk和mjk分别为水泵j-k静扬程、内阻系数和扬程指数。ωjk(t)为t时刻水泵j-k的相对转速，即工作转速与额定转速之比。
3.4)根据节点连续性方程和步骤3.1)中求得的下一时间步长各管道中分段点的流量和水头，代入相应边界条件求解各节点流量和水头。
3.5)根据步骤3.1)3.4)计算出的配水网能流分布，采用水泵功率计算公式确定下一时间步长水泵j-k的功率，即：
式中，ρ为水的密度。ηu,jk为水泵j-k的效率。
4)将计算得到的下一时间步长的水泵功率代入配电网中，作为电负荷。以三相三线接地配电网为例，假设分布式电源三相独立可控、负荷均采用中性点接地的Y形连接方式，则可采用如下三相潮流模型求解下一时间步长的配电网潮流：
式中，
，ψ∈{a,b,c}。
和
分别为t时刻节点j的
相注入有功和无功功率。
和
分别为t时刻节点j的
相电压模值和相角。
分别为t时刻配电网节点导纳矩阵对应节点j、 k及
、ψ相的电导和电纳。Ne为配电网节点数。
5)保存步骤3)，4)计算得到的能流数据。若仿真时间t小于仿真时长120s，返回步骤3)，若仿真时间t大于等于120s，仿真结束并输出数据。
实施例3：
参见附图3，使用一个由IEEE-13配电网和改进的5节点配水网组成的测试IEWS来说明实施例1所提出方法的有效性，步骤如下：
1)获取配电网数据：以字母B和N分别表示配电网母线和配水网节点，配水网的水泵接入配电网的母线B8，配水网的详细数据如表1，表2所示：
表1配水网节点数据
| 节点ID | 节点类型 | 标高(m) | 水负荷(m<sup>3</sup>/s) |
| N1(水源) | 水源 | 0 | - |
| N2 | 端节点 | 0 | 0 |
| N3 | 端节点 | 5 | 0.22 |
| N4 | 端节点 | 5 | 0.22 |
| N5 | 端节点 | 5 | 0.22 |
表2配水网支路数据
2)实验：
以附图3所示的由IEEE-13配电网和改进的5节点配水网组成的测试IEWS为仿真对象，以水泵支路N1-N2的水泵流量，水泵扬程，水泵功率和配电网中的水泵接入母线B8的A相电压模值为输出变量，并设置了如下两种场景：
场景1：基于稳态能流分析方法的电-水综合能源系统仿真结果；
场景2：基于电-水综合能源系统的准稳态能流分析方法的仿真结果。
如附录图10所示，在t＝10s时，配水网各节点水负荷减小至原来的80％。水负荷改变后，场景1(图中各子图的黑色虚线)中的水泵流量从0.67m3/s降低至0.53m3/s；水泵扬程从51.2m提升至 61.57m；水泵功率从371kW减小至357.6kW。由于稳态计算无法考虑两种稳态状态之间的暂态过渡过程，上述参数变化均发生在 10s-11s，且在后续时间均保持恒定。在场景2中，由于水负荷扰动而产生的水锤波在管道中以1000m/s的速度传播。距水泵最近的水负荷位于N3节点，两者距离2,000m管长(N2-N3)。因此，水泵流量，水泵扬程和水泵功率的变化(图中各子图的黑色实线)相比稳态计算结果晚2s。水锤波产生后，水泵功率的峰值和谷值分别为 372.5kW和335.6kW，与稳态计算结果的相对偏差分别为4.2％和 -6.2％。然后，水锤波在配水网中作往复运动，直到由于沿程水头损失而消散为止。因此，水泵流量，水泵扬程和水泵功率呈现幅值衰减的波动，并在稳态计算结果附近上下摆动。在t＝87s时，水泵流量，水泵扬程和水泵功率分别为61.57m，0.53m3/s和357.6kW，与稳态计算结果一致。
配水网水负荷的扰动通过水泵耦合环节的电功率变化传至配电网，如附录图11所示。由于场景1中配水网采用的是稳态计算模型，对应的配电网母线B8的A相电压模值反映出的变化同样发生在 10-11s，仅能反映变化前后的两种稳态潮流结果，即电压模值从0.972 升高至0.973。相比之下，场景2中，采用了本发明后能够进行考虑配水网暂态特性的IEWS准稳态能流分析。从图中黑色实线部分可知，受到配水网的暂态过渡过程影响，配电网的潮流会产生与水泵功率周期性变化相对应的波动，波动的峰值0.974和谷值0.972分别发生在19s和13s，当t＝87s时，电压模值趋于稳定值0.973，与场景1结果一致。上述仿真结果验证了所提方法的正确性和有效性。
Claims (6)
1.一种电-水综合能源系统的准稳态能流分析方法，其特征在于，包括以下步骤：
1)获取所述电-水综合能源系统的网络结构参数，并设置管道支路分段点数为I，管道分段步长Δx，时间步长为Δt，仿真时长为T；所述分段点包括上、下游端节点；
2)获取电-水综合能源系统的初始稳态能流，即时间t＝0时配水网的能流分布和配电网的潮流分布；
3)利用特征线法计算配水网在下一时段的暂态能流和水泵功率；
利用特征线法计算配水网在下一仿真时段的暂态能流的步骤包括：
3.1)基于已知的配水网的能流分布计算管道内下一时刻各分段点位置的管道流量q(t+1)和管道水头H(t+1)；
计算管道在任意时刻各分段点位置的管道流量q和管道水头H，步骤包括：
3.1.1)利用特征线法计算配水网暂态能流的主方程，即：
式中，q表示管道流量；H表示管道水头；g为重力加速度；c为水锤波速；D为管道直径；A为管道断面积；f为水力阻力系数；x为空间；t为时间；C+和C-分别表示沿管道两个不同方向的正特性方程和负特性方程；
3.1.2)将公式(1)和公式(2)改写如下：
3.1.3)联立公式(3)和公式(4)，计算管道在下一时刻各分段点位置的管道流量和管道水头，得到：
式中，(t+1)表示下一时间步长；下标i表示支路的分段节点；i＝1、2、3，...，I；(i-1)和(i+1)表示管道分段节点i的相邻节点；qi(t+1)表示下一时间步长节点i的管道流量；qi-1(t)、qi+1(t)表示时间t节点i-1和节点i+1的管道流量；Hi-1(t)、Hi+1(t)表示时间t节点i-1和节点i+1的管道水头；Hi(t+1)表示下一时间步长节点i的管道水头；
3.2)基于与支路两端分段点相连的元件确定相应端节点的连续性条件；
3.3)根据管道两端连接的不同设备确定边界条件；
3.4)将节点连续性方程、各管道中已知分段点的流量q(t)和水头H(t)代入边界条件中，求解得到下一时间步长各端节点的流量和水头；
3.5)根据步骤3.1)至步骤3.4)计算出的配水网能流分布，利用水泵功率计算公式确定下一时间步长的水泵功率；
4)将下一时间步长的水泵功率代入配电网中，作为电负荷，并计算下一时间步长的配电网能流；
5)保存能流计算结果，若当前时间t小于总时长T，则返回步骤3)，否则，输出配电网能流。
2.根据权利要求1所述的一种电-水综合能源系统的准稳态能流分析方法，其特征在于，管道分段步长Δx和时间步长为Δt满足如下关系：
式中，c为水锤波速。
3.根据权利要求1所述的一种电-水综合能源系统的准稳态能流分析方法，其特征在于，节点j的连续性条件分别如公式(8)和公式(9)所示，即：
式中，下标ta、ws、d、u和p分别表示蓄水池、水源、水负荷、水泵和管道；Nta(j)、Nws(j)、Nd(j)、Nu(j)和Np(j)分别表示蓄水池、水源、水负荷、水泵和管道的节点集合；
分别表示t时刻蓄水池节点、水源节点、水负荷节点、水泵上游节点、水泵下游节点、管道上游端节点、管道下游端节点的流量；
分别表示t时刻蓄水池节点、水源节点、水负荷节点、管道上游端节点、管道下游端节点的水头；
分别表示t时刻水泵出水和水泵进水；1表示上游端节点，I表示下游端节点。
4.根据权利要求1所述的一种电-水综合能源系统的准稳态能流分析方法，其特征在于，确定边界条件的步骤包括：
1)确定管道上游端和下游端的边界条件；
其中，管道上游端边界条件如下所示：
式中，qjk,1(t+1)和Hjk,1(t+1)分别表示支路j-k下一时间步长上游端节点的未知管道流量和管道水头；
其中，变量Cajk和变量Cnjk分别如下所示：
式中，qjk,2和Hjk,2分别表示支路j-k上游端节点的相邻分段节点的已知管道流量和管道水头；g为重力加速度；c为水锤波速；D为管道直径；A为管道断面积；f为水力阻力系数；
管道下游端边界条件如下所示：
式中，qjk,I(t+1)和Hjk,I(t+1)分别表示支路j-k下一时间步长下游端节点的未知管道流量和管道水头；
其中，变量Cpjk如下所示：
式中，qjk,I-1和Hjk,I-1分别表示支路j-k下游端节点的相邻分段节点的已知管道流量和管道水头；
2)确定水源边界条件，即：
式中，Hj(t)和qj(t)分别表示t时刻水源节点j的水头和流量；Hws(t)、qws(t)分别表示t时刻水源的管道水头和管道流量；
3)确定蓄水池的边界条件，步骤如下：
3.1)获取下一个时间步长与蓄水池相连支路的相应端的分段点的水头；
3.2)建立蓄水池的边界条件，即：
式中，Hj(t+1)为下一个时间步长与蓄水池节点j相连支路的相应端的分段点的水头；An为蓄水池相连支路的连接处的断面积；Co为蓄水池的孔口流量系数；qta(t)为t时刻蓄水池的管道流量；
基于蓄水池的边界条件，计算下一时间步长的蓄水池水头，即：
式中，Ata为蓄水池的底面积；
4)确定变速离心水泵边界条件，即：
式中，ΔHu,jk(t)为t时刻水泵j-k扬程，其值等于水泵的下游端节点与上游端节点的水头差；qu,jk(t)为t时刻水泵流量，其值等于水泵的上游端节点和下游端节点的流量；H0jk、cjk和mjk分别为水泵j-k静扬程、内阻系数和扬程指数；ωjk(t)为t时刻水泵j-k的相对转速，即工作转速与额定转速之比。
5.根据权利要求1所述的一种电-水综合能源系统的准稳态能流分析方法，其特征在于，t时刻的水泵j-k功率Pu,jk如下所示：
式中，ρ为水的密度；ηu,jk为水泵j-k的效率；qu,jk(t)为t时刻水泵流量，其值等于水泵的上游端节点和下游端节点的流量；ΔHu,jk(t)为t时刻水泵j-k扬程，其值等于水泵的下游端节点与上游端节点的水头差。
6.根据权利要求1所述的一种电-水综合能源系统的准稳态能流分析方法，其特征在于，t时刻的配电网能流如下所示：
式中，
ψ∈{a,b,c}；
和
分别为t时刻节点j的
相注入有功和无功功率；
和
分别为t时刻节点j的
相电压模值和相角；
为t时刻节点k的ψ相电压模值；
和
分别为配电网节点导纳矩阵对应节点j、节点k、
相、ψ相的电导和电纳；Ne为配电网节点数；
分别表示t时刻节点j的
相角和节点k的ψ相角。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011148025.8A CN112329217B (zh) | 2020-10-23 | 2020-10-23 | 一种电-水综合能源系统的准稳态能流分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011148025.8A CN112329217B (zh) | 2020-10-23 | 2020-10-23 | 一种电-水综合能源系统的准稳态能流分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112329217A CN112329217A (zh) | 2021-02-05 |
| CN112329217B true CN112329217B (zh) | 2022-07-15 |
Family
ID=74312242
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011148025.8A Expired - Fee Related CN112329217B (zh) | 2020-10-23 | 2020-10-23 | 一种电-水综合能源系统的准稳态能流分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112329217B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114154276B (zh) * | 2021-10-29 | 2024-09-24 | 天津大学 | 考虑蒸汽空腔断流弥合的长输供热系统停泵水锤数值仿真算法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103134672A (zh) * | 2013-02-04 | 2013-06-05 | 哈尔滨工业大学 | 一种管道排气节能效果的测试评价方法 |
| WO2015155555A1 (en) * | 2014-04-07 | 2015-10-15 | Horvath István | Water flow energy converter and blade wheel |
| CN106499001A (zh) * | 2016-12-23 | 2017-03-15 | 山东科源供排水设备工程有限公司 | 流量分区智能供水系统及其使用方法 |
| CN106874595A (zh) * | 2017-02-13 | 2017-06-20 | 中国水利水电科学研究院 | 基于节点参数化技术的输水管网计算方法 |
| CN107228042A (zh) * | 2017-07-31 | 2017-10-03 | 中国水利水电科学研究院 | 一种基于水锤泵的近岸波浪能发电系统和方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104659780B (zh) * | 2015-03-10 | 2018-01-19 | 国家电网公司 | 适用于消纳大规模分布式电源的能源互联网的实现方法 |
| US10276321B2 (en) * | 2015-04-06 | 2019-04-30 | Atom Power, Inc. | Dynamic coordination of protection devices in electrical distribution systems |
| CN107730129B (zh) * | 2017-10-24 | 2021-06-22 | 重庆大学 | 考虑光热热电联产与电锅炉的电-气-热互联系统风险评估方法 |
| CN108875269B (zh) * | 2018-07-09 | 2022-06-10 | 重庆大学 | 一种考虑电力系统多平衡机与天然气系统慢动态特性的电-气互联系统多时段能流计算方法 |
| CN109162810A (zh) * | 2018-08-23 | 2019-01-08 | 中国电力工程顾问集团西南电力设计院有限公司 | 一种用于分布式能源站的多能流测点系统 |
| CN110752605B (zh) * | 2019-10-29 | 2021-05-14 | 国家电网有限公司 | 一种电热耦合综合能源系统最优潮流计算方法 |
-
2020
- 2020-10-23 CN CN202011148025.8A patent/CN112329217B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103134672A (zh) * | 2013-02-04 | 2013-06-05 | 哈尔滨工业大学 | 一种管道排气节能效果的测试评价方法 |
| WO2015155555A1 (en) * | 2014-04-07 | 2015-10-15 | Horvath István | Water flow energy converter and blade wheel |
| CN106499001A (zh) * | 2016-12-23 | 2017-03-15 | 山东科源供排水设备工程有限公司 | 流量分区智能供水系统及其使用方法 |
| CN106874595A (zh) * | 2017-02-13 | 2017-06-20 | 中国水利水电科学研究院 | 基于节点参数化技术的输水管网计算方法 |
| CN107228042A (zh) * | 2017-07-31 | 2017-10-03 | 中国水利水电科学研究院 | 一种基于水锤泵的近岸波浪能发电系统和方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112329217A (zh) | 2021-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jiyun et al. | Development of an inline vertical cross-flow turbine for hydropower harvesting in urban water supply pipes | |
| Wang et al. | Water hammer simulation using explicit–implicit coupling methods | |
| Jain et al. | Investigations on pump running in turbine mode: A review of the state-of-the-art | |
| Satrio et al. | The influence of time step setting on the CFD simulation result of vertical axis tidal current turbine | |
| CN110783967A (zh) | 基于约束聚合的虚拟电厂出力可行域辨识方法及装置 | |
| Liu et al. | Evaluating the pressure fluctuations during load rejection of two pump-turbines in a prototype pumped-storage system by using 1D-3D coupled simulation | |
| CN108062599B (zh) | 一种电力系统和天然气系统协同运行的建模方法 | |
| Shahverdi et al. | CFD numerical simulation of Archimedes screw turbine with power output analysis | |
| CN106874595A (zh) | 基于节点参数化技术的输水管网计算方法 | |
| CN105069236A (zh) | 考虑风电场节点空间相关性的广义负荷联合概率建模方法 | |
| CN112329217B (zh) | 一种电-水综合能源系统的准稳态能流分析方法 | |
| Ye et al. | High-fidelity CFD simulations for the wake characteristics of the NTNU BT1 wind turbine | |
| Yu et al. | Hydraulic transients in the long diversion‐type hydropower station with a complex differential surge tank | |
| Naung et al. | High-fidelity CFD simulations of two wind turbines in arrays using nonlinear frequency domain solution method | |
| Fu et al. | Multi-objective optimization of guide vane closure scheme in clean pumped-storage power plant with emphasis on pressure fluctuations | |
| Yan et al. | Global matrix method for frequency-domain stability analysis of hydropower generating system | |
| CN110188996B (zh) | 水驱油藏能耗-产量-效益一体化表征方法 | |
| Guan et al. | Towards the integration of distributed renewables: Operation analysis of pumped storage system under off-design condition based on CFD | |
| Shojaeefard et al. | Analyzing the impact of blade geometrical parameters on energy recovery and efficiency of centrifugal pump as turbine installed in the pressure-reducing station | |
| Huang et al. | Transient simulation of hydropower station with consideration of three-dimensional unsteady flow in turbine | |
| Zhang et al. | Data-model-interactive enhancement-based Francis turbine unit health condition assessment using graph driven health benchmark model | |
| Liu et al. | Potential risks in balancing flexibility and investment of pumped storage plants: Hydraulic disturbances during transient processes in parallel operation of fixed-speed and variable-speed units sharing a diversion tunnel | |
| Schleicher et al. | Numerical optimization of a portable hydrokinetic turbine | |
| Chen et al. | Design and development of an eco-gate pump installation based on computational fluid dynamics | |
| Zeng et al. | Numerical prediction of hydraulic performance in model and homologous prototype Pelton turbine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220715 |