CN105932667B - 一种考虑网损分配的动态潮流仿真方法 - Google Patents
一种考虑网损分配的动态潮流仿真方法 Download PDFInfo
- Publication number
- CN105932667B CN105932667B CN201610331139.3A CN201610331139A CN105932667B CN 105932667 B CN105932667 B CN 105932667B CN 201610331139 A CN201610331139 A CN 201610331139A CN 105932667 B CN105932667 B CN 105932667B
- Authority
- CN
- China
- Prior art keywords
- unit
- power
- load
- network loss
- current
- 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
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000005283 ground state Effects 0.000 claims abstract description 10
- 230000033228 biological regulation Effects 0.000 claims description 23
- 238000000034 method Methods 0.000 claims description 16
- 238000004088 simulation Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 4
- 230000035772 mutation Effects 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000005611 electricity Effects 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- CVOFKRWYWCSDMA-UHFFFAOYSA-N 2-chloro-n-(2,6-diethylphenyl)-n-(methoxymethyl)acetamide;2,6-dinitro-n,n-dipropyl-4-(trifluoromethyl)aniline Chemical compound CCC1=CC=CC(CC)=C1N(COC)C(=O)CCl.CCCN(CCC)C1=C([N+]([O-])=O)C=C(C(F)(F)F)C=C1[N+]([O-])=O CVOFKRWYWCSDMA-UHFFFAOYSA-N 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑网损分配的动态潮流仿真方法,包括以下步骤:步骤1:在电力系统基态潮流基础上,计算系统当前的网损总量;步骤2:添加系统扰动,利用步骤1获得的网损作为当前系统网损近似值,进而求解当前系统不平衡功率近似值;步骤3:依据发电机组的机电特性求得系统当前近似频差,进而修正节点信息;步骤4:利用步骤3得到的节点信息进行潮流计算,得到新的电力系统网损总量,步骤5:如果网损变量较大,修正不平衡功率变化量,返回步骤3进行迭代计算;如果网损变量足够小,得到的潮流计算结果消除了平衡节点引入的误差,达到计算精度,迭代结果足够接近真实值,迭代计算停止。
Description
技术领域
本发明属于电力系统仿真领域,尤其涉及一种考虑网损分配的动态潮流仿真方法。
背景技术
目前在互联电网有功平衡层面,所涉及的仿真主要有两大类,一类是潮流计算,另一类是有功功率平衡仿真。然而,常规潮流计算只能单纯计算电力系统稳态运行情况下的潮流分布,不能反映系统的动态变化,且缺乏系统频率仿真。动态潮流这一概念的提出,刚好结合了这两大算法的优点。动态潮流的核心是潮流计算和频率计算。在动态潮流计算中,系统中由于功率扰动(切负荷、发电机增减出力)而产生的不平衡功率,按照各发电机和负荷的功频静特性系数在多台发电机及负荷之间进行分配,得到调整后的发电机出力和负荷的大小以及系统的频率连续变化的情况。相对于常规潮流算法来说,动态潮流算法能够在一段时间范围内对系统运行的实际情况进行更有效的连续潮流模拟,同时还能给出系统的频率波动情况。
现有的关于动态潮流的研究,主要是在常规潮流计算的基础上考虑了负荷和发电机的频率动态特性。然而,电力系统网损作为一种系统网络状态的依从变量,在潮流计算之前是未知的。因此,在常规潮流计算过程中,由网损变化导致的系统不平衡功率将全部由人工选定的平衡节点承担,使得计算结果依赖于平衡节点的选择。随着系统扰动的持续增加,系统中的网损也会非线性增长,而增长的网损造成的功率不平衡量仍由人工选定的平衡节点全部承担,会影响系统潮流计算以及系统频率仿真的精度,不甚合理。
发明内容
针对常规动态潮流在处理网损变量上的不足,本发明利用迭代分配方法合理地分配了网损变量。网损变量是在潮流计算的基础上得到的,是系统状态的依从变量。系统的基态潮流(即扰动发生前系统供需平衡且频率稳定的状态下的系统潮流)及添加的扰动(即节点信息变化主要指电力系统中节点i的发电机有功出力变化或有功负荷变化)是本发明中的输入量,也是已知量。具体的技术方案是:一种考虑网损分配的动态潮流仿真方法,包括以下步骤:
S1:在n节点电力系统基态潮流基础上计算系统当前的网损总量节点i的机组有功出力为PGi 0,负荷有功量为PLi 0;
S2:添加系统扰动:即节点i的发电机有功出力或负荷有功量突变,其中,发电机有功突变量为ΔPGr-i,负荷有功突变量为ΔPLr-i,此时系统中网损总量未知,利用步骤1获得的网损作为当前系统网损近似值,进而求解当前系统不平衡功率近似值,则当前系统不平衡功率Pacc为:
S3:依据发电机组的机电特性可知:
则系统当前近似频差Δω:
根据发电机组机电特性、一次调节特性以及负荷静频特性,经过时间Δt,电力系统中发电机组惯性出力变化量ΔPGTi,一次调节出力变化量ΔPGi,以及负荷在静频特性作用下的变化量ΔPLi:
对应当前电力系统节点i的机组有功出力为PGi 1、负荷有功量为PLi 1如下:
PGi 1=PGi 0+ΔPGr-i+ΔPGi+ΔPGTi
PLi 1=PLi 0+ΔPLr-i+ΔPLi
其中:Δt为仿真时间步长,单位:s,TJi为系统发电机组i的惯性时间常数,单位:s,TJ∑为系统所有发电机组的惯性时间常数之和,单位:s,KGi为发电机组i的单位调节功率,单位:MW/Hz,KGΣ为系统所有发电机组的单位调节功率之和,单位:MW/Hz,KLi为负荷的单位调节功率,单位:MW/Hz,KLi为系统所有负荷的单位调节功率之和,单位:MW/Hz;
S4:利用S3中得到的各机组有功出力以及负荷有功量进行潮流计算,得到新的电力系统网损总量则相对于Δt时间前的基态潮流,当前电力系统网损的总变化量为:
其中:k为扰动发生后潮流计算迭代次数,基态潮流对应迭代次数k=0;
S5:设置迭代精度值ε,精度值ε越小,计算结果精确度越高;如果说明经过时间Δt后,系统网损发生变化,得到的电力系统不平衡功率近似值与实际相比误差大时,该误差近似与电力系统网损的总变化量相等,则令:
ΔPGr-i=0,
ΔPLr-i=0
潮流计算迭代次数k加1,返回步骤3进行迭代计算;
如果迭代计算停止。
网损分配系数Kloss_i取值为各个机组的单位调节功率KGi,单位:MW/Hz,占系统所有机组单位调节功率之和ΣKGi的百分比,即:
当机组i参与频率调节的调节裕量为0时,令KGi=0,即:机组i对应当前网损分配系数Kloss_i=0。
节点信息变化主要指电力系统中节点i的发电机有功出力变化或有功负荷变化,且满足:
该方法中的电力系统基态潮流是指该电力系统在扰动发生前,整个系统供需平衡、频率稳定状态下的系统潮流。
由于采用了上述技术方案,本发明提供的一种考虑网损分配的动态潮流仿真方法,将动态潮流计算过程中由于网损变化导致的部分不平衡功率通过迭代分配给系统相应机组,而不是完全由人工选定的平衡节点承担,进而提高了潮流计算以及系统频率仿真的精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本方法的流程图;
图2为本发明方法的逻辑结构示意图。
具体实施方式
为使本发明的技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述:
如图1和图2所示的一种考虑网损分配的动态潮流仿真方法,具体包括以下步骤:
S1:在n节点电力系统基态潮流基础上计算系统当前的网损总量节点i的机组有功出力为PGi 0,负荷有功量为PLi 0;
S2:添加系统扰动:即节点i的发电机有功出力或负荷有功量突变,其中,发电机有功突变量为ΔPGr-i,负荷有功突变量为ΔPLr-i,此时系统中网损总量未知,利用步骤1获得的网损作为当前系统网损近似值,进而求解当前系统不平衡功率近似值,则当前系统不平衡功率Pacc为:
S3:依据发电机组的机电特性可知:
则系统当前近似频差Δω:
根据发电机组机电特性、一次调节特性以及负荷静频特性,经过时间Δt,电力系统中发电机组惯性出力变化量ΔPGTi,一次调节出力变化量ΔPGi,以及负荷在静频特性作用下的变化量ΔPLi:
对应当前电力系统节点i的机组有功出力为PGi 1、负荷有功量为PLi 1如下:
PGi 1=PGi 0+ΔPGr-i+ΔPGi+ΔPGTi
PLi 1=PLi 0+ΔPLr-i+ΔPLi
其中:Δt为仿真时间步长,单位:s,TJi为系统发电机组i的惯性时间常数,单位:s,TJΣ为系统所有发电机组的惯性时间常数之和,单位:s,KGi为发电机组i的单位调节功率,单位:MW/Hz,KGΣ为系统所有发电机组的单位调节功率之和,单位:MW/Hz,KLi为负荷的单位调节功率,单位:MW/Hz,KLi为系统所有负荷的单位调节功率之和,单位:MW/Hz;
S4:利用S3中得到的各机组有功出力以及负荷有功量进行潮流计算,得到新的电力系统网损总量则相对于Δt时间前的基态潮流,当前电力系统网损的总变化量为:
其中:k为扰动发生后潮流计算迭代次数,基态潮流对应迭代次数k=0;
S5:设置迭代精度值ε,精度值ε越小,计算结果精确度越高;如果说明经过时间Δt后,系统网损发生变化,得到的电力系统不平衡功率近似值与实际相比误差大时,该误差近似与电力系统网损的总变化量相等,则令:
ΔPGr-i=0,
ΔPLr-i=0
潮流计算迭代次数k加1,返回步骤3进行迭代计算;
如果迭代计算停止。
网损分配系数Kloss_i取值为各个机组的单位调节功率KGi,单位:MW/Hz,占系统所有机组单位调节功率之和ΣKGi的百分比,即:
当机组i参与频率调节的调节裕量为0时,令KGi=0,即:机组i对应当前网损分配系数Kloss_i=0。
该方法所采用的动态潮流算法应用于节点信息发生变化的电力系统当中,且对系统的实时变化进行仿真计算。而其中涉及的节点信息变化主要指电力系统中节点i的发电机有功出力变化或有功负荷变化,且满足:
该方法中的电力系统基态潮流是指该电力系统在扰动发生前,整个系统供需平衡、频率稳定状态下的系统潮流;
该方法涉及的仿真步长Δt应足够短,使得在实时仿真计算过程中,只能考虑一次调节作用,而二次调节并未开始作用,且认为一次调节在扰动发生时立即作用,即忽略一次调节的延迟时间。
实施例:
现以10机39节点系统为例,现系统负荷突减(即系统扰动)300MW。系统负荷突减时刻为0时刻,突减前系统稳定运行情况下网损为42.34MW;系统负荷突减后,不对网损进行处理以及采用本方法对网损进行处理,结果如下:
本发明通过对系统网损总变化量的迭代分配,消除了由人工选定的平衡节点承担系统网损变化量导致的计算误差,同时通过迭代算法,使得网损变化量的分配结果达到无限趋近于系统实际运行状态的网损分配结果,提高了动态潮流在仿真过程中的精度。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (1)
1.一种考虑网损分配的动态潮流仿真方法,其特征在于包括以下步骤:
S1:在n节点电力系统基态潮流基础上计算系统当前的网损总量节点i的机组有功出力为PGi 0,负荷有功量为PLi 0;
S2:添加系统扰动:即节点i的发电机有功出力或负荷有功量突变,其中,发电机有功突变量为ΔPGr-i,负荷有功突变量为ΔPLr-i,此时系统中网损总量未知,利用步骤1获得的网损作为当前系统网损近似值,进而求解当前系统不平衡功率近似值,则当前系统不平衡功率Pacc为:
S3:依据发电机组的机电特性可知:
则系统当前近似频差Δω:
根据发电机组机电特性、一次调节特性以及负荷静频特性,经过时间Δt,电力系统中发电机组惯性出力变化量ΔPGTi,一次调节出力变化量ΔPGi,以及负荷在静频特性作用下的变化量ΔPLi:
对应当前电力系统节点i的机组有功出力为PGi 1、负荷有功量为PLi 1如下:
PGi 1=PGi 0+ΔPGr-i+ΔPGi+ΔPGTi
PLi 1=PLi 0+ΔPLr-i+ΔPLi
其中:n为电力系统中的节点个数;N为电力系统中的发电机组个数,Δt为仿真时间步长,单位:s,TJi为系统发电机组i的惯性时间常数,单位:s,TJ∑为系统所有发电机组的惯性时间常数之和,单位:s,KGi为发电机组i的单位调节功率,单位:MW/Hz,KGΣ为系统所有发电机组的单位调节功率之和,单位:MW/Hz,KLi为负荷的单位调节功率,单位:MW/Hz,KLΣ为系统所有负荷的单位调节功率之和,单位:MW/Hz;
S4:利用S3中得到的各机组有功出力以及负荷有功量进行潮流计算,得到新的电力系统网损总量则相对于Δt时间前的基态潮流,当前电力系统网损的总变化量为:
其中:k为扰动发生后潮流计算迭代次数,基态潮流对应迭代次数k=0;
S5:设置迭代精度值ε,精度值ε越小,计算结果精确度越高;如果说明经过时间Δt后,系统网损发生变化,得到的电力系统不平衡功率近似值与实际相比误差大时,该误差近似与电力系统网损的总变化量相等,则令:
潮流计算迭代次数k加1,返回步骤3进行迭代计算;
如果迭代计算停止;
网损分配系数Kloss_i取值为各个机组的单位调节功率KGi,单位:MW/Hz,占系统所有机组单位调节功率之和ΣKGi的百分比,即:
当机组i参与频率调节的调节裕量为0时,令KGi=0,即:机组i对应当前网损分配系数Kloss_i=0。
节点信息变化主要指电力系统中节点i的发电机有功出力变化或有功负荷变化,且满足:
或
该方法中的电力系统基态潮流是指该电力系统在扰动发生前,整个系统供需平衡、频率稳定状态下的系统潮流。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610331139.3A CN105932667B (zh) | 2016-05-18 | 2016-05-18 | 一种考虑网损分配的动态潮流仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610331139.3A CN105932667B (zh) | 2016-05-18 | 2016-05-18 | 一种考虑网损分配的动态潮流仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105932667A CN105932667A (zh) | 2016-09-07 |
CN105932667B true CN105932667B (zh) | 2019-02-12 |
Family
ID=56840640
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610331139.3A Expired - Fee Related CN105932667B (zh) | 2016-05-18 | 2016-05-18 | 一种考虑网损分配的动态潮流仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105932667B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106451431B (zh) * | 2016-11-01 | 2022-07-08 | 中国电力科学研究院 | 一种混合交易模式下的网损分摊方法 |
CN108075475A (zh) * | 2016-11-14 | 2018-05-25 | 中国电力科学研究院 | 电力系统直流故障后受端电网动态潮流分析方法及系统 |
CN111049124B (zh) * | 2018-10-12 | 2022-09-23 | 中国电力科学研究院有限公司 | 一种电网直流闭锁后系统功率调整方法及系统 |
CN109586346B (zh) * | 2018-12-04 | 2022-06-07 | 深圳供电局有限公司 | 一种交直流大电网动态等值的平衡机出力调整方法和系统 |
CN110474369A (zh) * | 2019-08-16 | 2019-11-19 | 南京理工大学 | 考虑风电参与一次调频的动态潮流分析方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100461578C (zh) * | 2007-04-06 | 2009-02-11 | 清华大学 | 基于动态连续潮流的电力系统静态电压稳定的评估方法 |
US9281756B2 (en) * | 2011-11-11 | 2016-03-08 | Varentec, Inc. | Power flow controller with a fractionally rated back-to-back converter |
-
2016
- 2016-05-18 CN CN201610331139.3A patent/CN105932667B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105932667A (zh) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105932667B (zh) | 一种考虑网损分配的动态潮流仿真方法 | |
CN109256801B (zh) | 虚拟同步发电机虚拟惯量和虚拟阻尼系数自适应控制方法 | |
US20150149143A1 (en) | Apparatus, methods and systems for parallel power flow calculation and power system simulation | |
CN107240918B (zh) | 一种风电接入的电力系统等值简化方法 | |
CN108923460A (zh) | 微网虚拟同步机多机并联动态一致响应的参数配置方法 | |
CN109149620B (zh) | 一种自储能多端柔直系统控制方法及系统 | |
CN109962495B (zh) | 一种超低频振荡扰动源定位及抑制方法 | |
Apostolopoulou et al. | Balancing authority area model and its application to the design of adaptive AGC systems | |
CN107579543A (zh) | 一种基于分层控制策略的孤岛微电网分布式协调控制方法 | |
Ademola-Idowu et al. | Optimal design of virtual inertia and damping coefficients for virtual synchronous machines | |
CN104485670B (zh) | 孤立电网中电压敏感性工业负荷时变阻尼特性的控制方法 | |
Hansen et al. | A dynamic market mechanism for markets with shiftable demand response | |
CN104036320A (zh) | 一种基于改进粒子群算法的微网系统动态经济调度方法 | |
CN109103927A (zh) | 提高一次调频动态响应特性的调速系统pid控制器参数整定方法 | |
CN104979833A (zh) | 考虑地区电网无功支撑能力的静态电压稳定性判断方法 | |
Wang et al. | Phase balancing in power distribution network with data center | |
CN109301814A (zh) | 一种接入电网风电容量分析方法和系统 | |
Dou et al. | An improved CPF for static stability analysis of distribution systems with high DG penetration | |
CN108964061A (zh) | 一种考虑负荷频率电压静特性的含风电交直流电力系统概率动态连续潮流新方法 | |
CN103618339B (zh) | 自动发电控制中总调节功率分配方法 | |
Duan et al. | Hierarchical power flow control in smart grids: Enhancing rotor angle and frequency stability with demand-side flexibility | |
Panda et al. | Real-coded genetic algorithm for robust power system stabilizer design | |
CN108206544A (zh) | 基于一致性协议的多微网协调控制方法 | |
CN107402515A (zh) | 一种优化同步发电机励磁系统pid参数的方法 | |
CN112332409A (zh) | 一种电力系统输电断面潮流调整方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20190212 Termination date: 20210518 |
|
CF01 | Termination of patent right due to non-payment of annual fee |