CN110502859A - 用于电气耦合园区综合能源系统的多速率动态仿真方法 - Google Patents
用于电气耦合园区综合能源系统的多速率动态仿真方法 Download PDFInfo
- Publication number
- CN110502859A CN110502859A CN201910803105.3A CN201910803105A CN110502859A CN 110502859 A CN110502859 A CN 110502859A CN 201910803105 A CN201910803105 A CN 201910803105A CN 110502859 A CN110502859 A CN 110502859A
- Authority
- CN
- China
- Prior art keywords
- model
- electrical couplings
- energy system
- gas
- distribution network
- 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.)
- Granted
Links
- 238000010168 coupling process Methods 0.000 title claims abstract description 91
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 91
- 230000008878 coupling Effects 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005094 computer simulation Methods 0.000 title abstract description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 84
- 239000003345 natural gas Substances 0.000 claims abstract description 44
- 238000009826 distribution Methods 0.000 claims abstract description 38
- 238000013178 mathematical model Methods 0.000 claims abstract description 34
- 239000007789 gas Substances 0.000 claims description 78
- 230000005611 electricity Effects 0.000 claims description 30
- 238000004088 simulation Methods 0.000 claims description 30
- 239000000243 solution Substances 0.000 claims description 29
- 238000002347 injection Methods 0.000 claims description 8
- 239000007924 injection Substances 0.000 claims description 8
- 239000000446 fuel Substances 0.000 claims description 5
- 238000002485 combustion reaction Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract 2
- 230000001360 synchronised effect Effects 0.000 description 15
- 239000011159 matrix material Substances 0.000 description 8
- 238000013016 damping Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 239000000567 combustion gas Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及电气耦合技术领域,具体公开一种用于电气耦合园区综合能源系统的多速率动态仿真方法,所述方法包括:包括:建立电气耦合园区综合能源系统及其对应的电气耦合数学模型;采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解。本发明提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法,既具有较高的计算精度又有较高计算效率。
Description
技术领域
本发明涉及电气耦合技术领域,尤其涉及一种用于电气耦合园区综合能源系统的多速率动态仿真方法。
背景技术
随着能源与环境问题的日益突出,能够实现多种能源形式融合并高效利用的电气耦合综合能源系统已经成为当今能源技术领域的研究热点。其中,电力能源清洁高效、易于转化利用、可控性强,是综合能源系统的核心能源形式;而天然气因其经济环保、储量丰富等特征,同样被认为是未来综合能源体系中的重要组成部分。由于电力和燃气的动态特征和行为模式差异巨大,在综合能源系统中,对这两种能源形式进行准确、可靠的协同建模与仿真分析成为成为研究重点。
目前,电气系统由配电网、微型燃气轮机机组及P2G装置(P2G,即Power to Gas,概念源于德国可再生能源导入量的扩大,利用风力发电、太阳能发电等的剩余电力电解水生成氢,然后提供给现有的燃气管道网络;或者利用电力、水及大气中的CO2,通过甲烷化反应制造甲烷,提供燃气)组成,配电网模型多种多样,常用的一种是将配电网中各元件及负荷等效为恒定的阻抗,利用节点分析法求解整个系统如下所示:
I=Y
式中:I表示阻抗网络中各节点注入电流矩阵,Y表示网络的节点导纳矩阵,U表示网络各点的电压矩阵。
天然气管网模型主要以下式表示:
式中:A为管道截面积,p为管道节点压力,M为管道的流量,D为管道的直径,c为管道内的声速,λ为摩擦系数。
配电网和天然气系统采用微型燃气轮机机组和P2G装置进行双向耦合,微型燃气轮机机组由原动机和同步发电机组成,具体模型如下:
同步发电机模型:
式中:δ为同步电机的功角,ω*为电机的转速标幺值,ωN为电机的额定转速,Pm为输入的机械功率的标幺值,Pe为电磁功率的标幺值,D为同步电机阻尼系数,H为惯性时间常数的一半。
原动机模型:
Pm=f(p,M)
式中:Pm表示燃气轮机的输出功率,p表示电磁功率,M为管道流量信号。
将P2G装置看作配电网中的用电器,P2G模型如下所示:
式中:M表示P2G装置产生的气体流量,P为P2G装置需要的电功率,μ为转换效率,H为天然气燃烧的热值。
对电气耦合园区综合能源系统进行仿真,就是求解电气耦合园区综合能源系统各个节点处的电流或者电压等电力状态参数以及流量或者压力等气体状态参数的过程。
通过仿真,可以预知整个电气耦合园区综合能源系统各个节点处的状态参数,从而指导实际生产和运行。
现有的仿真方法大多是采用较小的步长,对天然气系统和电气系统在相同的时间尺度上依次进行求解,虽然求解的精度较高,但是计算所耗费的时间长,占用的空间和资源较多。
因此,急需针对研究出一种既能满足计算精度又有较高计算效率的仿真方法。
发明内容
本发明的一个目的在于,提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法,既具有较高的计算精度又有较高计算效率。
为达以上目的,本发明提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法,包括:
建立电气耦合园区综合能源系统及其对应的电气耦合数学模型;
采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解。
优选的,所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型包括:
确定配电网中线路的阻抗及变压器参数,然后建立电力能源系统拓扑连接模型;
确定天然气管网中各支路的长度、直径及摩擦系数,然后建立天然气能源系统拓扑连接模型;
分别确定微型燃气轮机机组和P2G装置二者的位置、容量及动态参数,通过所述微型燃气轮机机组和P2G装置耦合所述天然气能源系统拓扑连接模型和电力能源系统拓扑连接模型,得到所述电气耦合园区综合能源系统对应的电气耦合拓扑连接模型;
针对所述电气耦合拓扑连接模型建立电气耦合数学模型。
优选的,所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型还包括:
指定所述电气耦合数学模型的仿真参数;所述仿真参数包括电气耦合数学模型的仿真终止时间T、配电网模型的仿真步长Δt1、天然气管网模型的仿真步长Δt2、天然气管网模型的空间步长Δx,且满足Δt2=NΔt1,N为正整数,并设置仿真的扰动事件信息。
优选的,电气耦合数学模型包括配电网模型、天然气管网模型、同步发电机模型、原动机模型和P2G模型,所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型之后,还包括:
将所述天然气管网模型和同步发电机模型转化为差分化形式。
优选的,所述采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解包括:
设定外部循环仿真时刻初值tn=0以及设定电气耦合数学模型中若干状态变量的初值;
设置配电网模型的求解次数s=1且s为正整数;
设置当前时刻t=tn+sΔt1,以tn时刻天然气管网末端出口流量Mk(tn)作为微型燃气轮机机组燃料系统的注入流量,以仿真步长Δt1求解配电网模型,得到微型燃气轮机机组的转速、原动机机械功率、发电机电磁功率、功角、天然气管网与微型燃气轮机连接点的压力p(t)、电网中各点的电压及P2G流向天然气管网中的流量;使s=s+1;
判断t时刻配电网模型是否有扰动发生,若有则修改电力能源系统拓扑连接模型的拓扑结构和动态参数,否则直接进入下一步;
判断s是否大于或等于N+1,若是,进入下一步;
令外部循环仿真时刻tn=tn+Δt2,将求得的天然气管网与微型燃气轮机连接点的压力p(tn)作为天然气管网的末端压力,P2G注入天然气管网中的流量作为边界条件,以仿真步长Δt2求解天然气管网模型,得到天然气管网系统各节点压力及各支路流量,并更新天然气管网流向微型燃气轮机的流量Mk(tn);
判断在tn+Δt2时刻天然气管网模型是否发生扰动,若有扰动,则根据扰动信息修改天然气能源系统拓扑连接模型的拓扑结构和动态参数,否则直接进行下一步骤;
判断仿真时间tn是否大于仿真终止时刻T,若是则进入下一步;
输出仿真结果,仿真任务结束。
本发明的有益效果在于:提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法,可以用于电气耦合技术领域中的电气耦合数学模型求解情形,通过适当减少天然气管网模型的求解次数,可以以较小的计算精度损失换取较大的计算速度提升,能有效提高仿真效率。
附图说明
为了更清楚地说明本实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1是本发明实施例提供的用于电气耦合园区综合能源系统的多速率动态仿真方法的流程图;
图2是本发明实施例提供的电气耦合园区综合能源系统的拓扑结构图;
图3~图10是采用本实施例所提供的多速率动态仿真方法和采用传统的单速率求解方法求解所得的各状态参数的对比示意图。
具体实施方式
为使得本发明的目的、特征、优点能够更加的明显和易懂,下面将结合本实施例中的附图,对本实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本实施例提供一种用于电气耦合园区综合能源系统的多速率动态仿真方法,适用于电气耦合领域中的应用场景,可以提高电气耦合园区综合能源系统模型的求解效率,所述用于电气耦合园区综合能源系统的多速率动态仿真方法由一种监控装置来执行,通过软件和/或硬件实现。
图1是本实施例一提供的用于电气耦合园区综合能源系统的多速率动态仿真方法的流程图。
参见图1,所述用于电气耦合园区综合能源系统的多速率动态仿真方法包括如下步骤:
S10:参见图1,建立电气耦合园区综合能源系统及其对应的电气耦合数学模型;
具体地,S10包括:
S101:确定配电网中线路的阻抗及变压器参数,然后建立电力能源系统拓扑连接模型1;
S102:确定天然气管网中各支路的长度、直径及摩擦系数,然后建立天然气能源系统拓扑连接模型2;
本实施例中,天然气管网的各个参数可参照下表1进行确定。
表1传输管道参数表
S103:分别确定微型燃气轮机机组3和P2G装置4二者的位置、容量及动态参数,通过所述微型燃气轮机机组和P2G装置耦合所述天然气能源系统拓扑连接模型2和电力能源系统拓扑连接模型1,得到所述电气耦合园区综合能源系统对应的电气耦合拓扑连接模型;
本实施例中,同步电机的各个参数可参照下表2进行确定。
表2同步电机参数表
本实施例中,微型燃气轮机的各个参数可参照下表3进行确定。
表3微型燃气轮机参数表
参数名称 | 参数值 |
微型燃气轮机基准功率 | 250kW |
有功功率参考值 | 0.85p.u. |
阻尼系数 | 0.03 |
燃料系统延迟时间常数 | 0.1s |
燃料系统延迟时间常数 | 1.0s |
负荷限定时间常数 | 3.0s |
负荷限制 | 1.2 |
温度循环控制增益 | 1.0 |
功率控制比例增益 | 0.01 |
功率控制积分增益 | 5 |
本实施例中,P2G装置的各个参数可参照下表4进行确定。
表4 P2G装置参数表
参数名称 | 参数值 |
天然气热值 | 35588kJ/m<sup>3</sup> |
P2G转化效率 | 0.62 |
天然气密度 | 0.7174kg/m<sup>3</sup> |
S104:针对所述电气耦合拓扑连接模型建立电气耦合数学模型;
一般地,电气耦合数学模型包括配电网模型、天然气管网模型、同步发电机模型、原动机模型和P2G模型。可以理解的是,不同的电气耦合拓扑连接模型可以对应不同的电气耦合数学模型,为了方便说明,本实施例以背景技术中所提及的常用的电气耦合数学模型进行说明。即本实施例中:
配电网模型为:
I=YU
式中:I表示阻抗网络中各节点注入电流矩阵,Y表示网络的节点导纳矩阵,U表示网络各点的电压矩阵。
天然气管网模型为:
式中:A为管道截面积,p为管道节点压力,M为管道的流量,D为管道的直径,c为管道内的声速,λ为摩擦系数。
同步发电机模型:
式中:δ为同步电机的功角,ω*为电机的转速标幺值,ωN为电机的额定转速,Pm为输入的机械功率的标幺值,Pe为电磁功率的标幺值,D为同步电机阻尼系数,H为惯性时间常数的一半。
原动机模型为:
Pm=f(p,M)
式中:Pm表示燃气轮机的输出功率,p表示电磁功率,M为管道流量信号。
P2G模型为:
式中:M表示P2G装置产生的气体流量,P为P2G装置需要的电功率,μ为转换效率,H为天然气燃烧的热值。
S105:指定所述电气耦合数学模型的仿真参数;所述仿真参数包括电气耦合数学模型的仿真终止时间T、配电网模型的仿真步长Δt1、天然气管网模型的仿真步长Δt2、天然气管网模型的空间步长Δx,且满足Δt2=NΔt1,N为正整数,并设置仿真的扰动事件信息。
S20:电气耦合数学模型包括配电网模型、天然气管网模型、同步发电机模型、原动机模型和P2G模型,将所述天然气管网模型和同步发电机模型转化为差分化形式;
具体地,
1)同步电机差分化模型为:
式中,δ为同步电机的功角,ω*为电机的转速标幺值,ωN为电机的额定转速,Pm为输入的机械功率的标幺值,Pe为电磁功率的标幺值,D为同步电机阻尼系数,H为惯性时间常数的一半。
2)天然气管网差分化模型为:
式中,pm表示每一微元管道的末端压力,pk表示每一微元管道的始端压力,Mk表示每一微元管道的流入流量,Mm表示每一微元管道的流出流量,A为管道截面积,D为管道的直径,c为管道内的声速,λ为摩擦系数。
S30:采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解。
具体地,S30包括:
S301:设定外部循环仿真时刻初值tn=0以及设定电气耦合数学模型中若干状态变量的初值。
具体地,若干状态变量可以包括天然气管道的支路流量和节点压力以及微型燃气轮机机组的转速、原动机机械功率、发电机电磁功率和功角等。
S302:设置配电网模型的求解次数s=1且s为正整数;
S303:设置当前时刻t=tn+sΔt1,以tn时刻天然气管网末端出口流量Mk(tn)作为微型燃气轮机机组燃料系统的注入流量,以仿真步长Δt1求解配电网模型,得到微型燃气轮机机组的转速、原动机机械功率、发电机电磁功率、功角、天然气管网与微型燃气轮机连接点的压力p(t)、电网中各点的电压及P2G流向天然气管网中的流量;使s=s+1;
具体地,S303包括:
S3031:求解同步发电机模型,根据t时刻微型燃气轮机输出的机械功率Pm(t)和电网电磁功率Pe(t),求解同步电机差分化模型,得到t+Δt1时刻电机的功角δ(t+Δt1)以及微型燃气轮机机组的转速ωr*(t+Δt1);
S3032:求解注入电网侧的电流,根据同步电机功角δ(t+Δt1)和端口电压U(t),求得端口有功功率Pe(t+Δt1)、无功功率Qe(t+Δt1)及电流I(t+Δt1),如下式所示:
式中,Pe和Qe为微型燃气轮机机组注入电网的有功和无功功率,U为注入点的电压,δ为功角,为注入电流的共轭,Se为注入电网的视在功率,Xd为同步电机的d轴同步电抗,Eq为发电机空载电动势;
S3033:求解配电网模型,根据同步电机端口电流I(t+Δt1),更新电网节点注入电流矩阵,求解配电网模型得到电网中各节点的电压;
S3034:求解原动机模型,根据天然气管网末端出口流向微型燃气轮机的流量Mk(tn),求解差分化的微型燃气轮机模型,得到微型燃气轮机输出的机械功率Pm(t+Δt1)以及天然气管网与微型燃气轮机连接点的压力p(t+Δt1);
S3035:求解P2G模型,根据P2G装置在电力能源系统拓扑连接模型接入点的电压,求解P2G模型得到其注入天然气管网中的流量。
S304:判断t时刻配电网模型是否有扰动发生,若有则修改电力能源系统拓扑连接模型的拓扑结构和动态参数,否则直接进入下一步;
S305:判断s是否大于或等于N+1,若是,进入下一步,否则跳转至步骤S303;
S306:令外部循环仿真时刻tn=tn+Δt2,将求得的天然气管网与微型燃气轮机连接点的压力p(tn)作为天然气管网的末端压力,P2G注入天然气管网中的流量作为边界条件,以仿真步长Δt2求解天然气管网模型,得到天然气管网系统各节点压力及各支路流量,并更新天然气管网流向微型燃气轮机的流量Mk(tn);
具体地,以仿真步长Δt_2求解所得的天然气管网模型为联立天然气管网系统的元件约束和拓扑约束后的系统方程:
式中,A为系数矩阵,b为常数向量,p表示各点的压力向量,M表示各支路的流量向量,其中,M又可以写为M=[Mk;Mm]T
式中,Mk表示管道的流入流量向量,Mm表示管道的流出流量向量。
S307:判断在tn+Δt2时刻天然气管网模型是否发生扰动,若有扰动,则根据扰动信息修改天然气能源系统拓扑连接模型的拓扑结构和动态参数,否则直接进行下一步骤;
S308:判断仿真时间tn是否大于仿真终止时刻T,若是则进入下一步,否则跳转至步骤S302;
S309:输出仿真结果,仿真任务结束。
可以理解的是,所述仿真结果包括电气耦合拓扑连接模型以及求解所得的各状态参数。
结合上述步骤描述,下面给出具体算例:
本例以Matlab编程语言实现了用于电气耦合园区综合能源系统的多速率动态仿真方法,对本实施例提供的多速率动态仿真方法进行验证与分析,并与传统的单速率求解仿真方法仿真结果进行比较。仿真测试的硬件平台为Intel Core(TM)i5-6300HQ CPU@2.30GHz,8GB RAM的4核PC机;软件环境为64位Windows 10操作系统。
本算例中天然气经过50m管道后一部分再经100m管道流向一个恒定的天然气负荷,另一部分经过150m管道流向微型燃气轮机机组发电,配电网模型采用IEEE 33节点算例。
仿真场景设置如下:
设置仿真时间T=1000s,仿真步长Δt1=0.001s,Δt2=0.05s,Δx=25m,天然气流量的基值为50kg/s,600s时分别发生天然气负荷流量由50kg/s减少到45kg/s的扰动和电网25号节点阻抗标幺值增加0.1+j0.1的扰动,并将结果与单速率求解的结果进行对比,对比结果参见图3~图10。
需要说明的是,在传统的单速率求解步骤中,总是对配电网模型和天然气管网模型进行同步求解,即每次对配电网模型进行求解后,立刻使用求解所得的配电网模型对天然气管网模型进行求解,这样的做法可以提高求解精度,但是会导致运算次数较多,极大地降低了运算效率。本实施例采用的求解方法是,每对配电网模型进行N次求解才对天然气管网模型进行一次求解,这样虽然会损失一定的计算精度,但是却可以极大的减少运算次数,从而提高运算效率。针对上述参数进行对比,采用本实施例提供的多速率动态仿真方法所使用的运算时间为556.58s,用传统的单速率求解方法所使用的运算时间为827.22s,运算时间减少了大概30%。图3~图10为采用本实施例所提供的多速率动态仿真方法和采用传统的单速率求解方法求解所得的各状态参数的对比示意图,由图3~图10可知,采用本实施例所提供的多速率动态仿真方法后,各状态参数与采用传统的单速率求解方法求解所得的数值虽有差异,但是差异极小,均未达1%,但是运算时间却可以减少超过30%,其实际应用价值较高。
本实施例提供的用于电气耦合园区综合能源系统的多速率动态仿真方法,可以用于电气耦合技术领域中的电气耦合数学模型求解情形,减少天然气管网模型的求解次数,可以以较小的计算精度损失换取较大的计算速度提升,能有效提高仿真效率。
在本申请所提供的实施例中,应该理解到,所揭露的系统,单元,装置和方法,可以通过其它的方式实现。例如,以上所描述的所有实施例仅仅是示意性的,例如,上述单元或者模块等的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元、模块以及组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个计算机可读存储介质中,包括若干指令用以使得一台终端设备(可以是手机,笔记本,或者其他电子设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种用于电气耦合园区综合能源系统的多速率动态仿真方法,其特征在于,包括:
建立电气耦合园区综合能源系统及其对应的电气耦合数学模型;
采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解。
2.根据权利要求1所述的用于电气耦合园区综合能源系统的多速率动态仿真方法,其特征在于,所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型包括:
确定配电网中线路的阻抗及变压器参数,然后建立电力能源系统拓扑连接模型;
确定天然气管网中各支路的长度、直径及摩擦系数,然后建立天然气能源系统拓扑连接模型;
分别确定微型燃气轮机机组和P2G装置二者的位置、容量及动态参数,通过所述微型燃气轮机机组和P2G装置耦合所述天然气能源系统拓扑连接模型和电力能源系统拓扑连接模型,得到所述电气耦合园区综合能源系统对应的电气耦合拓扑连接模型;
针对所述电气耦合拓扑连接模型建立电气耦合数学模型。
3.根据权利要求2所述的用于电气耦合园区综合能源系统的多速率动态仿真方法,其特征在于,所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型还包括:
指定所述电气耦合数学模型的仿真参数;所述仿真参数包括电气耦合数学模型的仿真终止时间T、配电网模型的仿真步长Δt1、天然气管网模型的仿真步长Δt2、天然气管网模型的空间步长Δx,且满足Δt2=NΔt1,N为正整数,并设置仿真的扰动事件信息。
4.根据权利要求3所述的用于电气耦合园区综合能源系统的多速率动态仿真方法,其特征在于,电气耦合数学模型包括配电网模型、天然气管网模型、同步发电机模型、原动机模型和P2G模型,所述建立电气耦合园区综合能源系统及其对应的电气耦合数学模型之后,还包括:
将所述天然气管网模型和同步发电机模型转化为差分化形式。
5.根据权利要求4所述的用于电气耦合园区综合能源系统的多速率动态仿真方法,其特征在于,所述采用多次求解配电网模型后才进行一次天然气管网模型求解的方式对所述电气耦合数学模型进行参数求解包括:
设定外部循环仿真时刻初值tn=0以及设定电气耦合数学模型中若干状态变量的初值;
设置配电网模型的求解次数s=1且s为正整数;
设置当前时刻t=tn+sΔt1,以tn时刻天然气管网末端出口流量Mk(tn)作为微型燃气轮机机组燃料系统的注入流量,以仿真步长Δt1求解配电网模型,得到微型燃气轮机机组的转速、原动机机械功率、发电机电磁功率、功角、天然气管网与微型燃气轮机连接点的压力p(t)、电网中各点的电压及P2G流向天然气管网中的流量;使s=s+1;
判断t时刻配电网模型是否有扰动发生,若有则修改电力能源系统拓扑连接模型的拓扑结构和动态参数,否则直接进入下一步;
判断s是否大于或等于N+1,若是,进入下一步;
令外部循环仿真时刻tn=tn+Δt2,将求得的天然气管网与微型燃气轮机连接点的压力p(tn)作为天然气管网的末端压力,P2G注入天然气管网中的流量作为边界条件,以仿真步长Δt2求解天然气管网模型,得到天然气管网系统各节点压力及各支路流量,并更新天然气管网流向微型燃气轮机的流量Mk(tn);
判断在tn+Δt2时刻天然气管网模型是否发生扰动,若有扰动,则根据扰动信息修改天然气能源系统拓扑连接模型的拓扑结构和动态参数,否则直接进行下一步骤;
判断仿真时间tn是否大于仿真终止时刻T,若是则进入下一步;
输出仿真结果,仿真任务结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910803105.3A CN110502859B (zh) | 2019-08-28 | 2019-08-28 | 用于电气耦合园区综合能源系统的多速率动态仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910803105.3A CN110502859B (zh) | 2019-08-28 | 2019-08-28 | 用于电气耦合园区综合能源系统的多速率动态仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110502859A true CN110502859A (zh) | 2019-11-26 |
CN110502859B CN110502859B (zh) | 2023-04-07 |
Family
ID=68588684
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910803105.3A Active CN110502859B (zh) | 2019-08-28 | 2019-08-28 | 用于电气耦合园区综合能源系统的多速率动态仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110502859B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111416349A (zh) * | 2020-04-22 | 2020-07-14 | 南方电网科学研究院有限责任公司 | 一种电气互联综合能源系统的协同规划方法 |
CN112257244A (zh) * | 2020-10-15 | 2021-01-22 | 国网江苏省电力有限公司扬州供电分公司 | 一种微型燃气轮机发电系统的建模方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107291990A (zh) * | 2017-05-24 | 2017-10-24 | 河海大学 | 基于电‑气互联综合能源系统暂态模型的能量流仿真方法 |
CN109242217A (zh) * | 2018-10-31 | 2019-01-18 | 广东电网有限责任公司 | 一种电-气能源多目标分布式交易协同优化方法和装置 |
CN109242366A (zh) * | 2018-11-06 | 2019-01-18 | 国电南瑞科技股份有限公司 | 一种电-气互联综合能源系统的多时段潮流优化方法 |
US20190067939A1 (en) * | 2016-05-18 | 2019-02-28 | China Electric Power Research Institute Company Limited | Multi-time-scale digital/analog hybrid simulation system and method for power distribution network and storage medium |
CN109886523A (zh) * | 2018-12-25 | 2019-06-14 | 清华大学 | 一种综合能源网动态模型多速率计算方法 |
CN109919493A (zh) * | 2019-03-11 | 2019-06-21 | 长沙理工大学 | 一种基于多能流故障耦合传播的综合能源系统风险识别方法 |
-
2019
- 2019-08-28 CN CN201910803105.3A patent/CN110502859B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190067939A1 (en) * | 2016-05-18 | 2019-02-28 | China Electric Power Research Institute Company Limited | Multi-time-scale digital/analog hybrid simulation system and method for power distribution network and storage medium |
CN107291990A (zh) * | 2017-05-24 | 2017-10-24 | 河海大学 | 基于电‑气互联综合能源系统暂态模型的能量流仿真方法 |
CN109242217A (zh) * | 2018-10-31 | 2019-01-18 | 广东电网有限责任公司 | 一种电-气能源多目标分布式交易协同优化方法和装置 |
CN109242366A (zh) * | 2018-11-06 | 2019-01-18 | 国电南瑞科技股份有限公司 | 一种电-气互联综合能源系统的多时段潮流优化方法 |
CN109886523A (zh) * | 2018-12-25 | 2019-06-14 | 清华大学 | 一种综合能源网动态模型多速率计算方法 |
CN109919493A (zh) * | 2019-03-11 | 2019-06-21 | 长沙理工大学 | 一种基于多能流故障耦合传播的综合能源系统风险识别方法 |
Non-Patent Citations (1)
Title |
---|
兰晓明等: "基于动态降阶模型的电力系统非线性电压预测控制", 《电力建设》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111416349A (zh) * | 2020-04-22 | 2020-07-14 | 南方电网科学研究院有限责任公司 | 一种电气互联综合能源系统的协同规划方法 |
CN111416349B (zh) * | 2020-04-22 | 2021-12-03 | 南方电网科学研究院有限责任公司 | 一种电气互联综合能源系统的协同规划方法 |
CN112257244A (zh) * | 2020-10-15 | 2021-01-22 | 国网江苏省电力有限公司扬州供电分公司 | 一种微型燃气轮机发电系统的建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110502859B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111008468B (zh) | 综合能源能量管理系统的测试方法及其测试系统 | |
CN103700036B (zh) | 一种适于电力系统多时间尺度的暂态稳定性投影积分方法 | |
CN102012956B (zh) | 一种考虑风电场输入风速和风向随机波动的风电场等值方法 | |
CN103810646B (zh) | 一种基于改进投影积分算法的有源配电系统动态仿真方法 | |
CN110851960B (zh) | 一种计及天然气网络影响的配电网运行可靠性评估方法 | |
CN107516895A (zh) | 配电网快速仿真方法、装置、存储介质及其计算机设备 | |
CN107390547B (zh) | 一种包含微网群的主动配电网性能测试方法 | |
CN103678900B (zh) | 一种用于区域配电网实时仿真的网络解耦计算方法 | |
CN106532709A (zh) | 一种含配电网重构的区域综合能源系统最优潮流计算方法 | |
CN110502859A (zh) | 用于电气耦合园区综合能源系统的多速率动态仿真方法 | |
Li et al. | Gradient descent iterative method for energy flow of integrated energy system considering multiple modes of compressors | |
CN204314716U (zh) | 一种用于风电场控制系统硬件在环测试平台 | |
CN115392035A (zh) | 电-气互联综合能源系统优化调度的改进方法 | |
Huang et al. | A multi-rate dynamic energy flow analysis method for integrated electricity-gas-heat system with different time-scale | |
CN104218681B (zh) | 一种用于降低孤岛微电网切负荷成本的控制方法 | |
CN107887912A (zh) | 一种静态潮流分析方法及静态潮流分析系统 | |
CN110705066A (zh) | 基于投影积分的气电耦合园区综合能源系统动态仿真方法 | |
CN112182905B (zh) | 一种用于综合能源系统的供热管网仿真方法和装置 | |
CN105787604A (zh) | 一种输配协同的配网态势快速感知方法 | |
Zhai et al. | Reliability analysis of power-gas integrated energy system based on dynamic simulation | |
CN113536565B (zh) | 一种综合能源多时空尺度协同仿真方法与系统 | |
Gusain et al. | Energy flexibility analysis using fmuworld | |
CN109149630A (zh) | 一种新能源电量消纳能力分析规划方法 | |
CN105226635B (zh) | 双馈式风力发电仿真系统 | |
CN105470946A (zh) | 一种含分布式电源的配电网机电暂态仿真算法 |
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 |