CN108808681A - 基于混合注入模型的光伏并网潮流计算方法 - Google Patents
基于混合注入模型的光伏并网潮流计算方法 Download PDFInfo
- Publication number
- CN108808681A CN108808681A CN201710300232.2A CN201710300232A CN108808681A CN 108808681 A CN108808681 A CN 108808681A CN 201710300232 A CN201710300232 A CN 201710300232A CN 108808681 A CN108808681 A CN 108808681A
- Authority
- CN
- China
- Prior art keywords
- node
- grid
- photovoltaic plant
- power
- mixed injection
- 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.)
- Pending
Links
- 238000002347 injection Methods 0.000 title claims abstract description 34
- 239000007924 injection Substances 0.000 title claims abstract description 34
- 238000004364 calculation method Methods 0.000 title claims description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims description 13
- 238000010586 diagram Methods 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000000243 solution Substances 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
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H02J3/383—
-
- 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于混合注入模型的光伏并网潮流确定方法。包括以下步骤:1)输入光伏电站、电网参数,形成并网系统节点导纳矩阵,并设置初始值;2)求雅克比矩阵及修正量;3)修正状态变量;4)求平衡节点注入有功功率和无功功率、PV节点无功功率、各支路潮流分布及损耗。潮流算法结果既能获得光伏电站对电网运行状态的影响,又能得到光伏电站自身的状态;将光伏电站计入网络参数,同时消去了光伏电站的一些非必要节点,简化了光伏电站模型,减少了网络节点个数;本算法既继承了常规电流注入模型迭代时雅克比矩阵PQ节点分块矩阵的非对角元素恒定不变的优点,又降低了雅克比矩阵的规模,节省存储空间。
Description
技术领域
本发明属于新能源发电领域,特别涉及一种基于混合注入模型的光伏并网潮流确定方法。
背景技术
近年来,兆瓦级大型光伏电站相继涌现,随着光伏电站容量的不断增大,大规模集中式高压并入输电网势必成为未来发展趋势。大规模光伏电站集中并网会对电网的潮流、网络损耗等带来不可忽视的影响。此外,光伏电站在储能的配合下可以有效参与电网的调峰调频工作,在未来,当光伏电站的容量足够大时,一方面为了使光伏电站产生的电能尽可能多的送入电网,另一方面又能使光伏电站能很好地实现调峰调频功能,有必要了解光伏电站内部光伏阵列潮流对并网系统潮流的影响,以给光伏电站制定调度计划提供理论依据,这样光伏电站在减少出力时,就可以视情况切掉部分光伏阵列,而不是全部停机。因此,在潮流计算过程中计入光伏电站参数,充分考虑到光伏电站潮流和电网潮流的相互影响是十分必要的。
目前对光伏电站并网潮流计算的研究基本都是针对配电网展开的,大多针对配电网中传统N-R法对网络电压初值敏感、雅克比矩阵奇异,前推回代法对PV节点失效等问题提出改进策略,而这些算法只能获得光伏电站对电网运行状态的影响,无法得到光伏电站自身的状态,因此本发明提供的基于混合注入模型的光伏并网潮流计算方法能有效地解决这一问题。
发明内容
本发明的目的在于提供一种基于混合注入模型的光伏并网潮流计算方法。
实现本发明目的的技术解决方案为:基于混合注入模型的光伏并网潮流计算方法,包括以下步骤:
步骤1、输入光伏电站、电网参数,形成并网系统节点导纳矩阵,并设置初始值;
步骤2、确定雅克比矩阵及修正量;
步骤3、对状态变量进行修正;
步骤4、确定平衡节点注入有功功率和无功功率、PV节点无功功率、各支路潮流分布及损耗。
本发明与现有技术相比,其显著优点在于:(1)潮流算法结果既能获得光伏电站对电网运行状态的影响,又能得到光伏电站自身的状态;(2)将光伏电站计入网络参数,同时消去了光伏电站的一些非必要节点,简化了光伏电站模型,减少了网络节点个数;(3)本方法既继承了常规电流注入模型迭代时雅克比矩阵PQ节点分块矩阵的非对角元素恒定不变的优点,又降低了雅克比矩阵的规模,节省存储空间;(4)针对混合注入模型的混合坐标问题,提出坐标转换策略,将潮流方程右侧均转换到同一坐标下,使潮流方程变得整齐,程序编写方便。
附图说明
图1是本发明基于混合注入模型的光伏并网潮流计算方法的过程图。
图2是本发明光伏电站并网简化接线图。
图3是本发明光伏电站并网交流潮流部分等值电路图。
图中编号所代表的含义为:1为输入光伏电站、电网参数,形成并网系统节点导纳矩阵,并设置初始值,2为求雅克比矩阵及修正量,3为修正状态变量,4为求平衡节点注入有功功率和无功功率、PV节点无功功率、各支路潮流分布及损耗。
具体实施方式
结合图1,本发明的一种基于混合注入模型的光伏并网潮流确定方法,包括以下步骤:
步骤1、输入光伏电站、电网参数,形成并网系统节点导纳矩阵,并设置初始值;所述光伏电站参数包括并网线路阻抗、并网变压器变比、光伏电站等效PV节点有功功率及节点电压幅值;电网参数包括电网平衡节点的节点电压幅值和相角、PV节点的有功功率和节点电压幅值、PQ节点的有功功率和无功功率、支路号、首末节点号及支路阻抗。
设置初始值时,将PV节点的节点电压相角设置为0,PQ节点的节点电压幅值设置为1,节点电压相角设置为0。
步骤2、确定雅克比矩阵及修正量;所述雅克比矩阵由9个分块矩阵组成,分别是注入电流的实部ΔIre、注入电流的虚部ΔIim、注入有功功率的不平衡量ΔP对三个状态变量求偏导,状态变量分别是PQ节点电压实部列向量e、PQ节点电压虚部列向量f、PV节点电压相角列向量θ。所述修正量的值由雅克比矩阵的逆矩阵与状态变量的乘积得到。
步骤3、对状态变量进行修正;对状态变量进行修正具体是确定修正量与原状态变量的和。
步骤4、确定平衡节点注入有功功率和无功功率、PV节点无功功率、各支路潮流分布及损耗。
本发明的方法将光伏电站计入网络参数,同时消去了光伏电站的一些非必要节点,简化了光伏电站模型,减少了网络节点个数。
下面进行更详细的描述。
本发明提出一种基于混合注入模型的光伏并网潮流计算方法,包括以下步骤:
步骤一、输入光伏电站、电网参数,形成并网系统节点导纳矩阵,并设置初始值;
步骤二、求雅克比矩阵及修正量;
步骤三、修正状态变量;
步骤四、求平衡节点注入有功功率和无功功率、PV节点无功功率、各支路潮流分布及损耗。
进一步,步骤一中,假设光伏电站含有d个光伏阵列,忽略光伏电站中的汇流箱、无功补偿设备等,其并网简化接线图如图2所示,忽略输电线路L和变压器等值接地支路,光伏电站并网系统交流潮流部分等值电路如图3所示。下面推导光伏电站内部化简策略,由图2不难看出,BT和BT0是联络节点,除变压器支路没有连接任何其他电源或负荷,可以消去。
(1)消去节点BT
将yL和yT0等值成一个导纳支路yLT,即可消去节点BT,yLT的表达式为
由电路关系可推得
(2)消去节点BT0
列写各变压器支路高、低侧电流方程,并运用基尔霍夫定律,消去变量可得
其中,
Y1=diag{yLT/k0 2 yT1 ··· yTc ··· yTd};
根据图2,由电路关系可推得
其中,
将式(4)代入式(3),有
其中,
式(5)即为消去了节点BT和BT0的光伏电站内部节点电压方程。令最终化简后的光伏电站内部的节点导纳矩阵为Yp,电网的节点导纳矩阵为Yg,节点电压列向量为Up,节点注入电流列向量为Ip,电网侧的节点电压列向量为Ug,节点注入电流列向量为Ig,则整个光伏电站并网系统的节点电压方程可表示为
其中,Ygp=Ypg=0。
进一步,步骤一中,将PV节点的节点电压相角设置为0,PQ节点的节点电压幅值设置为1,PQ节点的节点电压相角设置为0。
进一步,步骤二中,雅克比矩阵的数学模型可表示为:
其中,x=[eT fT θT]T,e为PQ节点电压实部列向量,f为PQ节点电压虚部列向量,θ为PV节点电压相角列向量,ΔIre和ΔIim分别为注入电流的实部和虚部不平衡量,ΔP为注入有功功率的不平衡量,各子块元素如下(i≠j):
FJ阵:
MJ阵:
NJ阵:
HJ阵:
WJ阵:
XJ阵:
RJ阵:
SJ阵:
TJ阵:
其中,
修正量的数学模型为Δx=[ΔeT ΔfT ΔθT]T=J-1f(x)。
进一步,步骤三中,修正状态变量的数学模型为xi+1=xi+Δxi,收敛要求为max|Δx|<ε,ε是给定的计算精度。
下面结合实施例对本发明做进一步详细的描述。
实施例
步骤一、以我国某光伏电站中的4列光伏阵列组成光伏电站,接入IEEE30节点经典算例的7节点,并网系统的等值电路参数如表1所示。潮流计算采用标幺值,系统的基准容量为100MW。将电网PV节点的节点电压相角设置为0,PQ节点的节点电压幅值设置为1,PQ节点的节点电压相角设置为0。
表1并网系统等值电路参数
参数 | 参数值 | 参数 | 参数值 |
yL | 10.5429-j8.9548 | U1 | 1.0943 |
yT | 0.2026-j5.8069 | U2 | 1.0315 |
yT1 | 0.0393-j0.2168 | U3 | 1.0824 |
yT2 | 0.0392-j0.3812 | U4 | 1.0356 |
yT3 | 0.0393-j0.2572 | P1 | 0.0046 |
yT4 | 0.0392-j0.3415 | P2 | 0.0036 |
k0 | 0.9167 | P3 | 0.0048 |
k1~k4 | 1.0000 | P4 | 0.0039 |
步骤二、求雅克比矩阵和修正量即
Δx=[ΔeT ΔfT ΔθT]T=J-1f(x)。
步骤三、修正状态变量xi+1=xi+Δxi。取ε=10-10,判断潮流是否满足收敛要求即max|Δx|<ε。
步骤四、求平衡节点注入有功功率和无功功率、PV节点无功功率、各支路潮流分布及损耗。
光伏电站内部潮流计算结果如表2所示。
表2光伏电站内部潮流计算结果
节点 | 电压幅值 | 电压相角 | 节点有功 | 节点无功 |
BT1 | 1.0943 | -6.3253 | 0.0046 | 0.0011 |
BT2 | 1.0315 | -6.4914 | 0.0036 | -0.0222 |
BT3 | 1.0824 | -6.3307 | 0.0048 | -0.0019 |
BT4 | 1.0356 | -6.3830 | 0.0039 | -0.0187 |
BT0 | 1.0865 | -7.2721 | 0 | 0 |
BT | 1.0023 | -7.4228 | 0 | 0 |
由上可知,本发明提出了一种基于混合注入模型的光伏并网潮流计算方法,将光伏电站计入网络参数,充分考虑光伏电站和电网两部分的相互影响。求解潮流时,基于混合注入模型列写潮流方程,采用N-R法迭代求解,既保证了潮流计算的速度,又进一步节省了内存,为开展大型光伏电站并网运行特性的实证研究提供了一个有效的技术工具。
Claims (6)
1.一种基于混合注入模型的光伏并网潮流确定方法,其特征在于,包括以下步骤:
步骤1、输入光伏电站、电网参数,形成并网系统节点导纳矩阵,并设置初始值;
步骤2、确定雅克比矩阵及修正量;
步骤3、对状态变量进行修正;
步骤4、确定平衡节点注入有功功率和无功功率、PV节点无功功率、各支路潮流分布及损耗。
2.如权利要求1所述基于混合注入模型的光伏并网潮流计算方法,其特征在于,步骤1中所述光伏电站参数包括并网线路阻抗、并网变压器变比、光伏电站等效PV节点有功功率及节点电压幅值;电网参数包括电网平衡节点的节点电压幅值和相角、PV节点的有功功率和节点电压幅值、PQ节点的有功功率和无功功率、支路号、首末节点号及支路阻抗。
3.如权利要求2所述基于混合注入模型的光伏并网潮流计算方法,其特征在于,步骤1中,设置初始值时,将PV节点的节点电压相角设置为0,PQ节点的节点电压幅值设置为1,节点电压相角设置为0。
4.如权利要求1所述基于混合注入模型的光伏并网潮流计算方法,其特征在于,步骤2中所述雅克比矩阵由9个分块矩阵组成,分别是注入电流的实部ΔIre、注入电流的虚部ΔIim、注入有功功率的不平衡量ΔP对三个状态变量求偏导,状态变量分别是PQ节点电压实部列向量e、PQ节点电压虚部列向量f、PV节点电压相角列向量θ。
5.如权利要求4所述基于混合注入模型的光伏并网潮流计算方法,其特征在于,步骤2中所述修正量的值由雅克比矩阵的逆矩阵与状态变量的乘积得到。
6.如权利要求1所述基于混合注入模型的光伏并网潮流计算方法,其特征在于,步骤3中对状态变量进行修正具体是确定修正量与原状态变量的和。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710300232.2A CN108808681A (zh) | 2017-05-02 | 2017-05-02 | 基于混合注入模型的光伏并网潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710300232.2A CN108808681A (zh) | 2017-05-02 | 2017-05-02 | 基于混合注入模型的光伏并网潮流计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108808681A true CN108808681A (zh) | 2018-11-13 |
Family
ID=64054157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710300232.2A Pending CN108808681A (zh) | 2017-05-02 | 2017-05-02 | 基于混合注入模型的光伏并网潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108808681A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109861230A (zh) * | 2018-12-28 | 2019-06-07 | 重庆大学 | 一种含光伏逆变电源的三相四线制低压配电网改进潮流计算方法 |
CN111478362A (zh) * | 2020-04-10 | 2020-07-31 | 国网河北省电力有限公司培训中心 | 一种光伏电站和电网的一体化仿真方法、系统和存储介质 |
CN112564084A (zh) * | 2020-10-22 | 2021-03-26 | 贵州电网有限责任公司 | 大规模分布式光伏接入的配电网电压稳定性快速确定方法 |
CN112803477A (zh) * | 2021-01-05 | 2021-05-14 | 国网内蒙古东部电力有限公司通辽供电公司 | 考虑光伏电站无功出力能力的光伏地区无功协调控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6154716A (en) * | 1998-07-29 | 2000-11-28 | Lucent Technologies - Inc. | System and method for simulating electronic circuits |
CN102593823A (zh) * | 2012-02-01 | 2012-07-18 | 中国电力科学研究院 | 一种基于叠加原理的配电网在线潮流计算方法 |
CN103124080A (zh) * | 2013-02-04 | 2013-05-29 | 中国电力科学研究院 | 一种光伏发电系统模型的建模方法 |
CN103199525A (zh) * | 2013-03-27 | 2013-07-10 | 河海大学 | 一种基于等效节点注入电流变换的配电网潮流计算方法 |
-
2017
- 2017-05-02 CN CN201710300232.2A patent/CN108808681A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6154716A (en) * | 1998-07-29 | 2000-11-28 | Lucent Technologies - Inc. | System and method for simulating electronic circuits |
CN102593823A (zh) * | 2012-02-01 | 2012-07-18 | 中国电力科学研究院 | 一种基于叠加原理的配电网在线潮流计算方法 |
CN103124080A (zh) * | 2013-02-04 | 2013-05-29 | 中国电力科学研究院 | 一种光伏发电系统模型的建模方法 |
CN103199525A (zh) * | 2013-03-27 | 2013-07-10 | 河海大学 | 一种基于等效节点注入电流变换的配电网潮流计算方法 |
Non-Patent Citations (1)
Title |
---|
王立国: "基于混合注入模型的光伏并网潮流计算研究", 《现代电力》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109861230A (zh) * | 2018-12-28 | 2019-06-07 | 重庆大学 | 一种含光伏逆变电源的三相四线制低压配电网改进潮流计算方法 |
CN109861230B (zh) * | 2018-12-28 | 2020-10-09 | 重庆大学 | 一种含光伏逆变电源的三相四线制低压配电网改进潮流计算方法 |
CN111478362A (zh) * | 2020-04-10 | 2020-07-31 | 国网河北省电力有限公司培训中心 | 一种光伏电站和电网的一体化仿真方法、系统和存储介质 |
CN111478362B (zh) * | 2020-04-10 | 2021-09-17 | 国网河北省电力有限公司培训中心 | 一种光伏电站和电网的一体化仿真方法、系统和存储介质 |
CN112564084A (zh) * | 2020-10-22 | 2021-03-26 | 贵州电网有限责任公司 | 大规模分布式光伏接入的配电网电压稳定性快速确定方法 |
CN112564084B (zh) * | 2020-10-22 | 2023-08-22 | 贵州电网有限责任公司 | 大规模分布式光伏接入的配电网电压稳定性快速确定方法 |
CN112803477A (zh) * | 2021-01-05 | 2021-05-14 | 国网内蒙古东部电力有限公司通辽供电公司 | 考虑光伏电站无功出力能力的光伏地区无功协调控制方法 |
CN112803477B (zh) * | 2021-01-05 | 2022-09-23 | 国网内蒙古东部电力有限公司通辽供电公司 | 考虑光伏电站无功出力能力的光伏地区无功协调控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018103317A1 (zh) | 一种含upfc电力系统的通用潮流计算方法 | |
US20060282239A1 (en) | Method of setting-up steady state model of VSC-based multi-terminal HVDC transmission system | |
CN106532711B (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
CN108808681A (zh) | 基于混合注入模型的光伏并网潮流计算方法 | |
CN103971026B (zh) | 主动配电网潮流计算通用方法 | |
CN108683191B (zh) | 一种下垂控制型孤岛微电网的三相潮流分析方法 | |
Seifi et al. | Power flow study and comparison of FACTS: Series (SSSC), Shunt (STATCOM), and Shunt-Series (UPFC) | |
CN109066694A (zh) | 含线间潮流控制器的电力系统的多目标潮流优化方法 | |
CN106159955B (zh) | 基于连续惩罚对偶分解的电力系统分布式最优潮流方法 | |
CN109301870B (zh) | 一种电力电子多馈入电力系统容量优化方法 | |
Zhao et al. | Distributed continuation power flow method for integrated transmission and active distribution network | |
CN111555359A (zh) | 一种孤岛微电网功率精准分配的二次控制方法 | |
CN111884259A (zh) | 一种考虑系统小干扰稳定特性的场站级风电机组自适应等值方法 | |
CN111799800A (zh) | 一种交直流混合配电网潮流计算方法 | |
Zhou et al. | Small-signal stability assessment of heterogeneous grid-following converter power systems based on grid strength analysis | |
CN110120673B (zh) | 基于戴维南等值参数辨识的分布式输配协同无功优化方法及系统 | |
CN112072692A (zh) | 一种新能源发电场站的阻抗等值方法及装置 | |
Miao et al. | Generalized steady-state model for energy router with applications in power flow calculation | |
CN107800137B (zh) | 一种基于mmc的upfc机电暂态仿真方法 | |
Xuefeng et al. | A brief analysis of distributed generation connected to distribution network | |
Yang et al. | Power flow calculation for unbalanced three-phase distribution network with DGs based on phase-sequence hybrid modeling | |
CN110323782B (zh) | 一种常规电源由新能源替代的最大可替换容量确定方法 | |
Obrenić et al. | A novel intervals-based algorithm for the distribution short-circuit calculation | |
CN108667008B (zh) | 一种用于电力系统次同步振荡分析的通用小信号建模方法 | |
Jianming et al. | An improved power flow algorithm using equation changing method for AC/DC power system with VSC-HVDC |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181113 |