CN109149645B - 一种含有双馈感应式风电机组电网的暂态稳定计算方法 - Google Patents

一种含有双馈感应式风电机组电网的暂态稳定计算方法 Download PDF

Info

Publication number
CN109149645B
CN109149645B CN201811150207.1A CN201811150207A CN109149645B CN 109149645 B CN109149645 B CN 109149645B CN 201811150207 A CN201811150207 A CN 201811150207A CN 109149645 B CN109149645 B CN 109149645B
Authority
CN
China
Prior art keywords
power
network
generator
equation
dfig
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
Application number
CN201811150207.1A
Other languages
English (en)
Other versions
CN109149645A (zh
Inventor
孙淑琴
苟鑫
李永彬
梁伟强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN201811150207.1A priority Critical patent/CN109149645B/zh
Publication of CN109149645A publication Critical patent/CN109149645A/zh
Application granted granted Critical
Publication of CN109149645B publication Critical patent/CN109149645B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H02J3/386—
    • 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/24—Arrangements for preventing or reducing oscillations of power in networks
    • 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/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及风力发电领域，具体地来讲为一种含有双馈感应式风电机组电网的暂态稳定计算方法，通过在计算时间段的每一迭代步长内交替切换求解描述同步电机转子运动状态的微分方程组和描述网络中各节点电压电流约束关系的代数方程组，计算结束时即可得到网络中各同步发电机功角随时间的变化曲线，进而分析系统在受到各种干扰后的暂态稳定性。将暂态稳定分析分解成代数方程组求解和微分方程组求解，降低了求解不收敛的风险；同时对于高阶代数方程组的求解，采用LU分解法，避免矩阵求逆运算，提高计算速度并减小占用计算机内存空间；对于微分方程的求解，采用改进欧拉法构建一个预估‑校正系统，提高计算精度。

Description

一种含有双馈感应式风电机组电网的暂态稳定计算方法
技术领域
本发明涉及风力发电领域，具体地来讲为一种含有双馈感应式风电机组电网的暂态稳定计算方法。
背景技术
科学技术作为新型产业形成与应用的侦察兵，研究具备一定应用规模的风电容量接入传统电网的并网条件、运行与控制理论、紧急故障处理方式等课题，对风电行业的健康快速发展以及已并网风电在传统电网中安全可靠经济运行具有重大意义。
以风力机为原动机的发电方式是一种不同于利用汽轮机或水轮机作为原动机的新型发电方式，双馈感应式风力发电机(Doubly Fed Induction Generator，DFIG)是现阶段实际使用最广泛的风电机组。随着容量不断增加，风电并入传统电网必将呈现新的问题，工程研究人员必须深入研究并采取有效应对措施。
总的来看，大规模风电并入网络对系统产生各种影响的原因可以总结为以下两个方面。一方面，单台风力发电机组的输出功率较小(目前一般是兆瓦级)，为了达到像汽轮机那样数百兆瓦的输出功率，往往将数百台风力发电机组集中互联，形成具有较大地理分布面积的风力发电场。然而，自然界的风能本身具有波动性和突变性，而且，受风电场区域内部地形环境、风电机组排列位置等因素的影响，在某一时刻到达风电场内每台风力发电机迎风面的风速都可能有所差异，使得风电场的输出功率不能一直保持恒定。因此，并网区域电网的电压可能产生波动并带来谐波污染，也正是基于这些问题，风力发电机组不能像汽轮发电机组和水轮发电机组那样承担整个电力系统的调频和调压任务。另一方面，风电机组一般具有特殊的运行控制策略，使得它的功率输出外特性会随着自然环境的改变发生变化。风电机组并入电网后改变了原电网的功率分布格局，原电网中同步发电机的输出功率发生变化，重新构成的电网应对各种干扰或紧急故障的能力也随之发生改变。
在由双馈感应式风电机组构成的大规模风电场集中接入对电力系统功角稳定性的影响研究方面，国内外也取得了初步结果。风电场接入点附近发生故障会导致本地或其他地区的风电机组保护动作，系统的暂态稳定性相比于风电机组接入前降低；而常规同步发电机附近发生故障，风电机组的接入对提高系统的稳定性有所帮助。还有研究采用相同容量的双馈风电机组风电场替换常规同步发电机，发现系统的暂态稳定性趋于变好。不少研究者在研究过程中还发现，双馈风电机组接入对电力系统暂态稳定性的影响不仅与风电机组本身的参量有关，还与所连接电网的拓扑结构、干扰或故障发生位置、干扰或故障类型等因素有关。
发明内容
本发明所要解决的技术问题在于提供一种含有双馈感应式风电机组电网的暂态稳定分析交替求解方法，提高计算精度，求解结果可根据分析系统在受到各种干扰后的暂态稳定性。
本发明是这样实现的，
一种含有双馈感应式风电机组电网的暂态稳定计算方法，将双馈感应式风电机组处理为负导纳模型，包括如下的步骤：
1)导入双馈感应式风电机组接入电网的拓扑结构以及潮流计算数据，根据网络稳态潮流计算结果和暂态导纳矩阵，计算发电机初始值和负荷等效导纳；
2)判断系统的运行情况，根据系统的实际运行情况计算对应的网络节点导纳矩阵，在此基础之上加入负荷等效导纳以修改网络节点导纳矩阵；
3)网络中接入的DFIG采用负导纳简化模型，根据系统的实际运行情况计算DFIG端口母线在相应运行情况下注入电网的有功功率和无功功率，并据此等效为负导纳模型；
4)将步骤3中DFIG等效的负导纳模型加入到所述步骤2)的网络节点导纳矩阵中，得到除开同步发电机全网络的节点导纳矩阵；
5)根据xy0坐标系和dq0坐标系之间的变换关系，将同步发电机在dq0坐标系中的电压平衡方程变换到xy0坐标系下；将基于步骤4中节点导纳矩阵的网络方程变换到xy0坐标系下；联立xy0坐标系下的电压平衡方程和节点导纳矩阵的网络方程到暂态稳定计算用的网络代数方程；
6)求解步骤5中得到的暂态稳定计算用的网络代数方程，得到网络中各节点的电压，进而计算各同步发电机的注入电网的电流和电磁功率；
7)用改进Euler法求解描述同步电机转子运动状态的微分方程，通过构建一个预估-校正系统计算每一迭代步长结束时刻的同步电机转子功角和角频率；
8)在计算时间段的每一迭代步长内交替切换求解描述同步电机转子运动状态的微分方程组和描述网络中各节点电压电流约束关系的代数方程组，将微分方程组的计算结果作为解代数方程组的已知值，将代数方程组的计算结果作为求解微分方程组的已知值，计算结束时得到网络中各同步发电机功角随时间的变化曲线，进而分析系统的暂态稳定性。
进一步地，在求解整个网络的代数方程时，包括如下的步骤：
21)同步发电机在dq0坐标系中的电压平衡方程变换到xy0坐标系后的方程如下：
Figure GDA0001849343510000041
利用上一迭代步求解微分方程的数值计算结果，计算方程中各参数
Figure GDA0001849343510000042
式中
Rai：第i台同步发电机定子绕组电阻，
X'di：第i台同步发电机d轴暂态电抗，
Xqi：第i台同步发电机q轴同步电抗，
如果是首次求解网络代数方程，δ值按式
Figure GDA0001849343510000043
进行计算；
22)根据电网的实际运行状况按式
Figure GDA0001849343510000044
列出除同步发电机以外网络的代数方程；
23)联立步骤步骤21)和步骤22)中的方程组，得到暂态稳定计算用的网络方程
Figure GDA0001849343510000045
网络方程中的节点导纳矩阵包含负荷、双馈感应式风电机组以及同步发电机的导纳。
24)利用步骤21)中计算出的系数代入式
Figure GDA0001849343510000046
计算各同步发电机机端节点注入的虚拟电流；
25)将24)中求得的发电机节点虚拟电流代入式
Figure GDA0001849343510000051
利用LU三角分解法，求解此高阶线性方程组，得到系统中所有节点电压的x分量和y分量；
26)将25)中计算所得的节点电压代入式
Figure GDA0001849343510000052
结合21)中计算出的参数，求得系统中所有节点注入电流的x分量和y分量；
27)根据注入电压和注入电流，按式
Figure GDA0001849343510000053
计算发电机输出的电磁功率。
进一步地，所述的LU三角分解法具体为：对系数矩阵A为n阶非奇异矩阵的方程Ax＝b，将系数矩阵A分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，使n阶线性方程组求解中不存在矩阵求逆运算，只有乘积与求和运算，最终得到高阶线性方程组Ax＝b的解，利用LU分解法实现占据内存较少并节约计算时间。
进一步地，当网络代数方程求解结束以后，步骤7)中用改进Euler法求解描述同步电机转子运动状态的微分方程式，包括如下步骤：
41)根据式
Figure GDA0001849343510000054
得到在微小时间间隔tk～tk+1初始时刻tk处状态变量δ和ω的变化率为
Figure GDA0001849343510000055
其中，发电机的电磁功率在本迭代步长内求解网络代数方程时得到；
42)按式
Figure GDA0001849343510000056
Δt＝tk+1-tk以得到在时间间隔tk～tk+1终止时刻tk+1处状态变量δ和ω的预估值为
Figure GDA0001849343510000061
43)将δ预估值
Figure GDA0001849343510000062
作为已知量，重新求解整个网络的代数方程，最后得到按预估值
Figure GDA0001849343510000063
计算所得的发电机端口电压和电流的预估值
Figure GDA0001849343510000064
Figure GDA0001849343510000065
和
Figure GDA0001849343510000066
进而，按式
Figure GDA0001849343510000067
计算发电机电磁功率的预估值
Figure GDA0001849343510000068
44)类比步骤41)，根据式
Figure GDA0001849343510000069
求得在时间间隔tk～tk+1终止时刻tk+1处状态变量δ和ω变化率的预估值为
Figure GDA00018493435100000610
45)根据式
Figure GDA00018493435100000611
求得在时间间隔tk～tk+1终止时刻tk+1处状态变量δ和ω变化率的校正值为
Figure GDA00018493435100000612
在时间间隔tk～tk+1内微分方程的求解完成后，所求状态变量δ的值作为下一时间间隔求解网络代数方程的已知条件，如此循环，完成整个仿真时段的计算。
进一步地，步骤2)中判断系统的运行情况包括：网络正常运行、发生某种故障以及故障清除。
进一步地，所述潮流计算数据的计算过程包括：将DFIG接入电网的母线处理为“类PQ”母线，具体为如下的步骤：
S1：在传统电力网络基础上，根据DFIG接入电网的实际情况更改整个电力网络的结构；
S2：置循环次数k＝0，输入DFIG的工作风速V，设置传统电力网络中各母线电压的初值以及风电场端口母线电压的初值
Figure GDA0001849343510000071
S3：根据DFIG的电磁功率-风速曲线求得风电机组输出的电磁功率Pe；
S4：根据DFIG的电磁功率-旋转角速度曲线求得风力发电机的转子旋转角速度ω；
S5：将DFIG等效为绕线式异步发电机模型，根据风力发电机的转子旋转角速度ω计算等效的绕线异步发电机的转差率s；
S6：根据S2、S3以及S4的DFIG出口母线电压、转差率s以及风电机组输出的电磁功率Pe，计算等效的绕线异步发电机的定子绕组向电网馈送的有功功率Ps和无功功率Qs；
S7：将DFIG的出口母线等效为稳态“类PQ”母线，将上述计算所得的有功功率Ps和无功功率Qs以PQ节点的形式作为DFIG出口母线注入电网的功率，采用牛顿-拉夫逊迭代算法计算整个网络中各母线电压，风电场端口母线电压
Figure GDA0001849343510000072
S8：对计算得到的风电场端口母线电压与循环次数k的风电场端口母线电压的初值进行比较判断，若满足
Figure GDA0001849343510000073
ε指双馈风电机组所接母线电压的迭代计算允许误差，则得到最终的风电场端口母线电压Us，并计算有功功率Ps和无功功率Qs；
S9：根据最终的有功功率Ps和无功功率Qs，计算得到网络中各发电机母线的注入功率，各支路传输功率及损耗功率。
进一步地，有功功率Ps通过解一元二次方程：aPs 2+bPs+c＝0得到，其中系数的计算公式如下：
Figure GDA0001849343510000081
其中
Xss＝Xs+Xm：定子绕组的漏抗与激磁漏抗之和，
Figure GDA0001849343510000082
DFIG的功率因素，
s：等效的绕线式异步发电机转差率，
Us：DFIG端口母线电压。
进一步地，无功功率Qs通过下列公式计算：
Figure GDA0001849343510000083
进一步地，将DFIG端口母线处理为“类PQ”母线包括：“类PQ”母线的注入功率在每一迭代步中是变化的，会随着该母线电压上一迭代步的迭代结果而变化。
进一步地，步骤S7中设迭代计算时DFIG出口母线电压第k次迭代结果为
Figure GDA0001849343510000084
第(k+1)次迭代结果为
Figure GDA0001849343510000085
若不满足
Figure GDA0001849343510000086
则，令k＝k+1，再次计算等效的绕线异步发电机定子绕组向电网馈送的有功功率Ps和无功功率Qs，继续迭代直到满足误差要求。
本发明在含有双馈感应式风电机组电网的稳态潮流计算基础之上，建立了DFIG暂态稳定求解简化模型，通过在计算时间段的每一迭代步长内交替切换求解描述同步电机转子运动状态的微分方程组和描述网络中各节点电压电流约束关系的代数方程组，计算结束时即可得到网络中各同步发电机功角随时间的变化曲线，进而分析系统在受到各种干扰后的暂态稳定性。
本发明与现有技术相比，有益效果在于：
1)将暂态稳定分析分解成代数方程组求解和微分方程组求解，降低了求解不收敛的风险；同时对于高阶代数方程组的求解，采用LU分解法，避免矩阵求逆运算，提高计算速度并减小占用计算机内存空间；对于微分方程的求解，采用改进欧拉法构建一个预估-校正系统，提高计算精度。
2)计算结果输出网络中各台同步发电机的功角摇摆曲线，通过对曲线进行分析，能够直观地判断系统在收到各种干扰之后是否保持同步稳定。
附图说明
图1为本发明实施例方法的流程图；
图2为本发明实施例中提供的xy0坐标系中节点电压和电流转换到dq0坐标系中，两坐标系的关系图；
图3为传统电网结构图；
图4为本发明一实施例中设定的DFIG的故障穿越特性图；
图5为本发明一实施例中经过交替求解输出4台同步发电机相对功角的摇摆曲线图；
图6为本发明一实施例中经过交替求解输出4台同步发电机转子角速度变化曲线图；
图7为本发明一实施例中新的同步发电机相对功角的摇摆曲线图
图8为本发明一实施例中新的转子角速度变化曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
参见图1，将双馈感应式风电机组处理为负导纳模型，包括如下的步骤：
1)导入双馈感应式风电机组接入电网的拓扑结构以及潮流计算数据，根据网络稳态潮流计算结果和暂态导纳矩阵，计算发电机初始值和负荷等效导纳；
2)判断系统的运行情况，根据系统的实际运行情况(网络正常运行、发生某种故障、故障清除)计算对应的网络节点导纳矩阵，在此基础之上加入负荷等效导纳以修改网络节点导纳矩阵；
3)网络中接入的DFIG采用负导纳简化模型，根据系统的实际运行情况(网络正常运行、发生某种故障、故障清除)计算DFIG端口母线在相应运行情况下注入电网的有功功率和无功功率，并据此等效为负导纳模型；步骤3是对步骤2的进一步补充，由于步骤2中只是对电网负荷和拓扑结构变化做出处理，并将这些改变加入到网络导纳矩阵中，没有对电网中的风电机组进行处理，所以在步骤3中，将DFIG利用负导纳简化模型加入到网络节点导纳矩阵中。
4)将步骤3中DFIG等效的负导纳模型加入到所述步骤2)的网络节点导纳矩阵中，得到除开同步发电机全网络的节点导纳矩阵；
5)根据xy0坐标系和dq0坐标系之间的变换关系，将同步发电机在dq0坐标系中的电压平衡方程变换到xy0坐标系下；将基于步骤4中节点导纳矩阵的网络方程变换到xy0坐标系下；联立xy0坐标系下的电压平衡方程和节点导纳矩阵的网络方程到暂态稳定计算用的网络代数方程；
6)求解步骤5中得到的暂态稳定计算用的网络代数方程，得到网络中各节点(包含同步发电机和双馈感应式风力发电机出口节点)的电压，进而计算各同步发电机的注入电网的电流和电磁功率；
7)用改进Euler法求解描述同步电机转子运动状态的微分方程，通过构建一个预估-校正系统计算每一迭代步长结束时刻的同步电机转子功角和角频率；
8)在计算时间段的每一迭代步长内交替切换求解描述同步电机转子运动状态的微分方程组和描述网络中各节点电压电流约束关系的代数方程组，将微分方程组的计算结果作为解代数方程组的已知值，将代数方程组的计算结果作为求解微分方程组的已知值，计算结束时得到网络中各同步发电机功角随时间的变化曲线，进而分析系统的暂态稳定性。
对网络中的双馈风电机组，等效处理为负导纳模型，模型具体建立过程如下：
通过采集国产某型号输出功率为1.5MW的DFIG在不同母线电压时的功率输出数据，归纳DFIG在故障期间以及故障清除后，其输出功率随时间的变化规律：设DFIG稳态正常运行时输出的有功为P0，无功为Q0，端口母线电压为U0，发生故障期间各量分别为P1、Q1和U1，故障清除后各量分别为P2、Q2和U2，则故障期间DFIG的功率变化规律可表示如下
Figure GDA0001849343510000111
Q1＝Q0·(1+γ) (2)
其中η为有功修正系数，可取范围为0.85-0.95，本发明中取值为0.90；γ为无功修正系数，可取范围为0-2.0，本发明中取值为1.0。
故障清除后DFIG的功率变化规律为
P2＝P1+p·t (3)
Q2＝Q0 (4)
其中p为有功爬升速度，一般限制在0.3MW/s以下，本发明中设置的爬升速度p＝0.2MW/s。
得到每个时刻DFIG注入电网的有功功率和无功功率后，将其等效为负导纳
Figure GDA0001849343510000121
以便将DFIG暂态等效模型连接到整个网络的节点导纳矩阵中。
步骤8)中，采用在计算时间段的每一迭代步长内交替切换求解描述同步电机转子运动状态的微分方程组和描述网络中各节点电压电流约束关系的代数方程组，将微分方程组的计算结果作为解代数方程组的已知值，将代数方程组的计算结果作为求解微分方程组的已知值，计算结束时得到网络中各同步发电机功角随时间的变化曲线，进而分析系统的暂态稳定性。
所谓交替求解法就是在每一个微小时间间隔tk～tk+1中，交替切换求解描述各同步发电机转子运动状态的微分方程组与描述电力网络电压电流约束关系的代数方程组，微分方程组的计算结果作为解代数方程组的已知值，同样，代数方程组的计算结果作为求解微分方程组的已知值。
首先：初始条件计算
在进行第一个迭代步长t0～t1的交替计算之前，必须根据稳态潮流计算结果计算t0时刻各发电机的初始条件。发电机所在节点t0时刻的电压(即端电压)为
Figure GDA0001849343510000122
和注入功率
Figure GDA0001849343510000123
则该节点t0时刻的注入电流为
Figure GDA0001849343510000124
为方便计算功角构建的虚拟电势
Figure GDA0001849343510000125
功角初始条件
Figure GDA0001849343510000126
旋转角速度初始条件
ω(0)＝1 (8)
为方便求得同步发电机q轴的暂态电动势，需要将上述xy0坐标系中节点电压和电流转换到dq0坐标系中，两坐标系的关系如图2所示，发电机端口电压和电流从xy0坐标系变换到dq0坐标系的变换公式为
Figure GDA0001849343510000131
Figure GDA0001849343510000132
现给出发电机q轴暂态电动势初始值的计算表达式(11)，在整个计算过程中可以认为发电机励磁系统的励磁很强，使q轴暂态电势一直恒定，为初始值。
E'q(0)＝Uq(0)+RaIq(0)+X'dId(0) (11)
输入到发电机的机械功率(即原动机的输出功率)按式(12)计算，在整个计算过程中，由于暂态稳定计算时间一般不是太长，可认为该机械功率一直保持初始值不变。
Figure GDA0001849343510000133
其次，求解网络代数方程组
设某台发电机所在节点的编号为i，若采用E'q＝C模型来描述同步发电机，则定子电压平衡方程式为
Figure GDA0001849343510000134
结合变换公式(9)和(10)，将式(9)转化到xy0坐标下可得
Figure GDA0001849343510000135
式(14)中各参数的具体表达式为
Figure GDA0001849343510000141
式中
Rai：第i台同步发电机定子绕组电阻
X'di：第i台同步发电机d轴暂态电抗
Xqi：第i台同步发电机q轴同步电抗
如果是首次求解网络代数方程，δ值按式
Figure GDA0001849343510000142
进行计算。
不考虑同步发电机，描述电力网络运行状态的代数方程如式(16)所示，等号左边为网络中节点的注入电流。
Figure GDA0001849343510000143
联合式(14)与(16)，消去同步发电机所连接节点的注入电流可得暂态稳定计算用的网络方程
Figure GDA0001849343510000144
方程式等号左端是方便计算用的假想电流，它的计算公式为
Figure GDA0001849343510000145
在求解整个网络的代数方程时，求解式(14)到(18)必须按照一定的顺序：
①利用微分方程的数值计算结果，按式(15)计算将会用到的参数，如果是首次求解网络代数方程，δ值按式
Figure GDA0001849343510000146
计算；
②将①中计算出的系数代入式(18)计算各发电机节点注入的虚拟电流；
③根据电网的实际运行状况按式(16)列出除发电机以外网络的代数方程(如果发生故障或为清除故障的断线操作，式中节点导纳矩阵应该是短路或断线时电力网络的节点导纳矩阵)；
④利用①中计算出的系数，按式(17)给出的规则修改③中求得的网络导纳矩阵，修改后的节点导纳矩阵实际上已经包含了负荷、DFIG以及同步发电机的等效导纳；
⑤将②中求得的发电机节点假想电流代入式(17)左端，可利用LU三角分解法，求解此高阶线性方程组，得到系统中所有节点电压的x分量和y分量；
⑥将⑤中计算所得的节点电压代入式(16)，结合①中计算出的参数，求得系统中所有节点注入电流的x分量和y分量；
⑦按式(19)计算发电机的输出功率
Figure GDA0001849343510000151
网络代数方程求解结束以后，便可以用改进Euler法求解描述同步电机转子运动状态的微分方程式，下面推导其求解的具体过程，
41)根据式
Figure GDA0001849343510000156
可以得到在微小时间间隔tk～tk+1初始时刻tk处状态变量δ和ω的变化率为
Figure GDA0001849343510000152
其中，同步发电机的电磁功率在求解网络代数方程时已按式(15)计算给出。
42)按式
Figure GDA0001849343510000153
Δt＝tk+1-tk可以得到在时间间隔tk～tk+1的终止时刻tk+1处状态变量δ和ω的预估值为
Figure GDA0001849343510000154
43)将按式(21)计算得到的δ预估值
Figure GDA0001849343510000155
作为已知量，重新求解整个网络的代数方程，最后得到按预估值
Figure GDA0001849343510000161
计算所得的发电机端口电压和电流的预估值
Figure GDA0001849343510000162
和
Figure GDA0001849343510000163
进而，按式(19)计算发电机电磁功率的预估值
Figure GDA0001849343510000164
44)类比步骤41)，根据式
Figure GDA0001849343510000165
可以求得在时间间隔tk～tk+1终止时刻tk+1处状态变量δ和ω变化率的预估值为
Figure GDA0001849343510000166
45)根据式
Figure GDA0001849343510000167
可以求得在时间间隔tk～tk+1终止时刻tk+1处状态变量δ和ω变化率的校正值为
Figure GDA0001849343510000168
至此，在时间间隔tk～tk+1内微分方程的求解已经完成，所求状态变量δ的值可作为下一时间间隔求解代数方程的已知条件，如此循环，便可完成整个仿真时段的计算。
本实施例中，潮流计算数据的计算方法为：将DFIG接入电网的母线处理为“类PQ”母线，包括如下的步骤：
S1：在传统电力网络基础上，根据DFIG接入电网的实际情况更改整个电力网络的结构；
S2：置循环次数k＝0，输入DFIG的工作风速V，设置传统电力网络中各母线电压的初值以及风电场端口母线电压的初值
Figure GDA0001849343510000169
S3：根据DFIG的电磁功率-风速曲线求得风电机组输出的电磁功率Pe；
S4：根据DFIG的电磁功率-旋转角速度曲线求得风力发电机的转子旋转角速度ω；
S5：将DFIG等效为绕线式异步发电机模型，根据风力发电机的转子旋转角速度ω计算等效的绕线异步发电机的转差率s；
S6：根据S2、S3以及S4的DFIG出口母线电压、转差率s以及风电机组输出的电磁功率Pe，计算等效的绕线异步发电机的定子绕组向电网馈送的有功功率Ps和无功功率Qs；
S7：将DFIG的出口母线等效为稳态“类PQ”母线，将上述计算所得的有功功率Ps和无功功率Qs以PQ节点的形式作为DFIG出口母线注入电网的功率，采用牛顿-拉夫逊迭代算法计算整个网络中各母线电压，风电场端口母线电压
Figure GDA0001849343510000171
S8：对计算得到的风电场端口母线电压与循环次数k的风电场端口母线电压的初值进行比较判断，若满足
Figure GDA0001849343510000172
ε指双馈风电机组所接母线电压的迭代计算允许误差，则得到最终的风电场端口母线电压Us，并计算有功功率Ps和无功功率Qs；
S9：根据最终的有功功率Ps和无功功率Qs，计算得到网络中各发电机母线的注入功率，各支路传输功率及损耗功率。
有功功率Ps通过解一元二次方程：aPs 2+bPs+c＝0得到，其中系数的计算公式如下：
Figure GDA0001849343510000173
其中
Xss＝Xs+Xm：定子绕组的漏抗与激磁漏抗之和，
Figure GDA0001849343510000181
DFIG的功率因素，
s：等效的绕线式异步发电机转差率，
Us：DFIG端口母线电压。
无功功率Qs通过下列公式计算：
Figure GDA0001849343510000182
将DFIG端口母线处理为“类PQ”母线包括：“类PQ”母线的注入功率在每一迭代步中是变化的，会随着该母线电压上一迭代步的迭代结果而变化。
步骤S7中设迭代计算时DFIG出口母线电压第k次迭代结果为
Figure GDA0001849343510000183
第(k+1)次迭代结果为
Figure GDA0001849343510000184
若不满足
Figure GDA0001849343510000185
则，令k＝k+1，再次计算等效的绕线异步发电机定子绕组向电网馈送的有功功率Ps和无功功率Qs，继续迭代直到满足误差要求。
应用实施例
本实施例将30台1.5MW的DFIG功率输出之和接入传统电网，传统电网结构见图3。每台DFIG的参数见表1，匹配变压器采用0.69kV/35kV电压等级，为简化计算，设匹配变压器为理想变压器并忽略损耗，风电场出口处的主变压器采用35kV/110kV电压等级，所用主变压器的具体参数如表2所示。从主变压器高压侧架设一条40km的LGJ-240型电力传输线，具体参数见表3，接入常规电网中编号为13的母线上。
DFIG风电场接入区域电网后，将风电场出口母线编号为14，在13和14号母线之间将增加一条变压器支路和一条传输线支路，将两者标幺处理后在本质上没有区别，所以将两者视为一条支路，主变压器高压侧母线不再单独编号。
设定DFIG的故障穿越特性如图4所示。现假设该案例发生故障，具体故障情况是：在仿真开始后的0.5秒，10号母线到11号母线这条支路靠近10号母线处发生永久性三相金属性接地故障；在仿真开始后的0.70秒，故障支路上的继电保护装置动作，切除10号母线到11号母线这条支路。经过交替求解方法输出4台同步发电机相对功角的摇摆曲线和转子角速度变化曲线如图5和图6所示。
改变故障支路上继电保护装置动作的时间，使其在仿真开始后0.83秒动作，其余故障参数保持不变，重新输入程序计算，得到新的同步发电机相对功角的摇摆曲线和转子角速度变化曲线分别如图7和图8所示。
暂态稳定分析后得出以下结论：
⑴随着故障切除时间的延长，相对功角摇摆曲线中的部分曲线幅值逐渐增大，转子角速度的波动情况也愈加剧烈，系统将逐渐失去稳定；
⑵通过比较DFIG接入电网前后，系统在相同故障干扰下失去稳定时间，可以判定双馈感应式风电机组接入后对系统暂态稳定性的影响。
表1 DFIG基本参数
Figure GDA0001849343510000191
表2主变压器基本参数
Figure GDA0001849343510000201
表3 LGJ-240型电力传输线参数
Figure GDA0001849343510000202
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种含有双馈感应式风电机组电网的暂态稳定计算方法，其特征在于，将双馈感应式风电机组处理为负导纳模型，包括如下的步骤：
1)导入双馈感应式风电机组接入电网的拓扑结构以及潮流计算数据，根据网络稳态潮流计算结果和暂态导纳矩阵，计算发电机初始值和负荷等效导纳；
2)判断系统的运行情况，根据系统的实际运行情况计算对应的网络节点导纳矩阵，在此基础之上加入负荷等效导纳以修改网络节点导纳矩阵；
3)网络中接入的DFIG采用负导纳简化模型，根据系统的实际运行情况计算DFIG端口母线在相应运行情况下注入电网的有功功率和无功功率，并据此等效为负导纳模型；
4)将步骤3中DFIG等效的负导纳模型加入到所述步骤2)的网络节点导纳矩阵中，得到除开同步发电机全网络的节点导纳矩阵；
5)根据xy0坐标系和dq0坐标系之间的变换关系，将同步发电机在dq0坐标系中的电压平衡方程变换到xy0坐标系下；将基于步骤4中节点导纳矩阵的网络方程变换到xy0坐标系下；联立xy0坐标系下的电压平衡方程和节点导纳矩阵的网络方程到暂态稳定计算用的网络代数方程；
6)求解步骤5)中得到的暂态稳定计算用的网络代数方程，得到网络中各节点的电压，进而计算各同步发电机的注入电网的电流和电磁功率；
7)用改进Euler法求解描述同步电机转子运动状态的微分方程，通过构建一个预估-校正系统计算每一迭代步长结束时刻的同步电机转子功角和角频率；
8)在计算时间段的每一迭代步长内交替切换求解描述同步电机转子运动状态的微分方程组和描述网络中各节点电压电流约束关系的代数方程组，将微分方程组的计算结果作为解代数方程组的已知值，将代数方程组的计算结果作为求解微分方程组的已知值，计算结束时得到网络中各同步发电机功角随时间的变化曲线，进而分析系统的暂态稳定性；
其中，在求解整个网络的代数方程时，包括如下的步骤：
21)同步发电机在dq0坐标系中的电压平衡方程变换到xy0坐标系后的方程如下：
Figure FDA0002399275410000021
利用上一迭代步求解微分方程的数值计算结果，计算方程中各参数
Figure FDA0002399275410000022
式中
Rai：第i台同步发电机定子绕组电阻，
X'di：第i台同步发电机d轴暂态电抗，
Xqi：第i台同步发电机q轴同步电抗，
如果是首次求解网络代数方程，δ值按式
Figure FDA0002399275410000023
进行计算；
22)根据电网的实际运行状况按式
Figure FDA0002399275410000024
列出除同步发电机以外网络的代数方程；
23)联立步骤步骤21)和步骤22)中的方程组，得到暂态稳定计算用的网络方程
Figure FDA0002399275410000025
网络方程中的节点导纳矩阵包含负荷、双馈感应式风电机组以及同步发电机的导纳；
24)利用步骤21)中计算出的系数代入式
Figure FDA0002399275410000031
计算各同步发电机机端节点注入的虚拟电流；
25)将24)中求得的发电机节点虚拟电流代入式
Figure FDA0002399275410000032
利用LU三角分解法，求解此高阶线性方程组，得到系统中所有节点电压的x分量和y分量；
26)将25)中计算所得的节点电压代入式
Figure FDA0002399275410000033
结合21)中计算出的参数，求得系统中所有节点注入电流的x分量和y分量；
27)根据注入电压和注入电流，按式
Figure FDA0002399275410000034
计算发电机输出的电磁功率；
当网络代数方程求解结束以后，步骤7中用改进Euler法求解描述同步电机转子运动状态的微分方程式，包括如下步骤：
41)根据式
Figure FDA0002399275410000035
得到在微小时间间隔tk～tk+1初始时刻tk处状态变量δ和ω的变化率为
Figure FDA0002399275410000036
其中，发电机的电磁功率在本迭代步长内求解网络代数方程时得到；
42)按式
Figure FDA0002399275410000037
Δt＝tk+1-tk以得到在时间间隔tk～tk+1终止时刻tk+1处状态变量δ和ω的预估值为
Figure FDA0002399275410000041
43)将δ预估值
Figure FDA0002399275410000042
作为已知量，重新求解整个网络的代数方程，最后得到按预估值
Figure FDA0002399275410000043
计算所得的发电机端口电压和电流的预估值
Figure FDA0002399275410000044
Figure FDA0002399275410000045
和
Figure FDA0002399275410000046
进而，按式
Figure FDA0002399275410000047
计算发电机电磁功率的预估值
Figure FDA0002399275410000048
44)类比步骤41)，根据式
Figure FDA0002399275410000049
求得在时间间隔tk～tk+1终止时刻tk+1处状态变量δ和ω变化率的预估值为
Figure FDA00023992754100000410
45)根据式
Figure FDA00023992754100000411
求得在时间间隔tk～tk+1终止时刻tk+1处状态变量δ和ω变化率的校正值为
Figure FDA00023992754100000412
在时间间隔tk～tk+1内微分方程的求解完成后，所求状态变量δ的值作为下一时间间隔求解网络代数方程的已知条件，如此循环，完成整个仿真时段的计算。
2.按照权利要求1所述的方法，其特征在于，所述的LU三角分解法具体为：对系数矩阵A为n阶非奇异矩阵的方程Ax＝b，将系数矩阵A分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，使n阶线性方程组求解中不存在矩阵求逆运算，只有乘积与求和运算，最终得到高阶线性方程组Ax＝b的解，利用LU分解法实现占据内存较少并节约计算时间。
3.按照权利要求1所述的方法，其特征在于，步骤2)中判断系统的运行情况包括：网络正常运行、发生某种故障以及故障清除。
4.按照权利要求1所述的方法，其特征在于，所述潮流计算数据的计算过程包括：将DFIG接入电网的母线处理为“类PQ”母线，具体为如下的步骤：
S1：在传统电力网络基础上，根据DFIG接入电网的实际情况更改整个电力网络的结构；
S2：置循环次数k＝0，输入DFIG的工作风速V，设置传统电力网络中各母线电压的初值以及风电场端口母线电压的初值
Figure FDA0002399275410000051
S3：根据DFIG的电磁功率-风速曲线求得风电机组输出的电磁功率Pe；
S4：根据DFIG的电磁功率-旋转角速度曲线求得风力发电机的转子旋转角速度ω；
S5：将DFIG等效为绕线式异步发电机模型，根据风力发电机的转子旋转角速度ω计算等效的绕线异步发电机的转差率s；
S6：根据S2、S3以及S4的DFIG出口母线电压、转差率s以及风电机组输出的电磁功率Pe，计算等效的绕线异步发电机的定子绕组向电网馈送的有功功率Ps和无功功率Qs；
S7：将DFIG的出口母线等效为稳态“类PQ”母线，将上述计算所得的有功功率Ps和无功功率Qs以PQ节点的形式作为DFIG出口母线注入电网的功率，采用牛顿-拉夫逊迭代算法计算整个网络中各母线电压，风电场端口母线电压
Figure FDA0002399275410000052
S8：对计算得到的风电场端口母线电压与循环次数k的风电场端口母线电压的初值进行比较判断，若满足
Figure FDA0002399275410000061
ε指双馈风电机组所接母线电压的迭代计算允许误差，则得到最终的风电场端口母线电压Us，并计算有功功率Ps和无功功率Qs；
S9：根据最终的有功功率Ps和无功功率Qs，计算得到网络中各发电机母线的注入功率，各支路传输功率及损耗功率。
5.按照权利要求4所述的方法，其特征在于，有功功率Ps通过解一元二次方程：
Figure FDA0002399275410000062
得到，其中系数的计算公式如下：
Figure FDA0002399275410000063
其中
Xss＝Xs+Xm：定子绕组的漏抗与激磁漏抗之和，
Figure FDA0002399275410000064
DFIG的功率因素，
s：等效的绕线式异步发电机转差率，
Us：DFIG端口母线电压。
6.按照权利要求4所述的方法，其特征在于，无功功率Qs通过下列公式计算：
Figure FDA0002399275410000065
7.按照权利要求4所述的方法，其特征在于，将DFIG端口母线处理为“类PQ”母线包括：“类PQ”母线的注入功率在每一迭代步中是变化的，会随着该母线电压上一迭代步的迭代结果而变化。
8.按照权利要求4所述的方法，其特征在于，步骤S7中设迭代计算时DFIG出口母线电压第k次迭代结果为
Figure FDA0002399275410000066
第(k+1)次迭代结果为
Figure FDA0002399275410000067
若不满足
Figure FDA0002399275410000071
则，令k＝k+1，再次计算等效的绕线异步发电机定子绕组向电网馈送的有功功率Ps和无功功率Qs，继续迭代直到满足误差要求。
CN201811150207.1A 2018-09-29 2018-09-29 一种含有双馈感应式风电机组电网的暂态稳定计算方法 Expired - Fee Related CN109149645B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811150207.1A CN109149645B (zh) 2018-09-29 2018-09-29 一种含有双馈感应式风电机组电网的暂态稳定计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811150207.1A CN109149645B (zh) 2018-09-29 2018-09-29 一种含有双馈感应式风电机组电网的暂态稳定计算方法

Publications (2)

Publication Number Publication Date
CN109149645A CN109149645A (zh) 2019-01-04
CN109149645B true CN109149645B (zh) 2020-05-05

Family

ID=64813803

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811150207.1A Expired - Fee Related CN109149645B (zh) 2018-09-29 2018-09-29 一种含有双馈感应式风电机组电网的暂态稳定计算方法

Country Status (1)

Country Link
CN (1) CN109149645B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110119523A (zh) * 2019-01-29 2019-08-13 中国电力科学研究院有限公司 一种基于功角等效法的网络代数方程求解预处理方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102609576B (zh) * 2012-01-19 2013-12-04 浙江大学 预估-校正数值积分的电力系统暂态稳定仿真方法
CN102611132B (zh) * 2012-02-27 2014-08-13 山东大学 双馈变速风电机组的附加频率控制器参数整定方法
CN102957166B (zh) * 2012-09-13 2015-04-29 中国电力科学研究院 一种基于轨迹灵敏度的风电配置比例快速计算方法
JP2017200275A (ja) * 2016-04-26 2017-11-02 株式会社日立製作所 電力系統解析システム及び方法
CN107317354B (zh) * 2017-06-30 2019-12-17 天津大学 一种含风电场的多机系统暂态功角稳定性分析方法

Also Published As

Publication number Publication date
CN109149645A (zh) 2019-01-04

Similar Documents

Publication Publication Date Title
Kalcon et al. Small-signal stability analysis of multi-terminal VSC-based DC transmission systems
CN107017646B (zh) 基于虚拟阻抗控制的双馈风机次同步振荡抑制方法
CN112671006A (zh) 一种评估海上风电场柔性直流输电系统谐振稳定性的方法
CN106026113A (zh) 一种具有无功自动补偿的微电网系统的监控方法
CN102856917A (zh) 一种配电网无功优化方法
CN109586337A (zh) 基于频域建模的vsc并网系统次同步振荡风险评估方法
CN105762816B (zh) 一种抑制次同步相互作用的方法
CN116706927B (zh) 一种四桥臂逆变器最优电压支撑参考电流计算方法
CN111654052A (zh) 基于动态相量法的柔直换流器建模装置及方法
Gao et al. A fast high-precision model of the doubly-fed pumped storage unit
CN111130135A (zh) 一种适用于高比例新能源接入下的电力系统惯量计算方法
Xu et al. Sub-synchronous frequency domain-equivalent modeling for wind farms based on rotor equivalent resistance characteristics
CN111027179A (zh) 一种计及辅助调频服务的双馈风电场的等值建模方法
CN114491886A (zh) 含多类型分布式新能源的主动配电网通用建模方法和装置
CN109149645B (zh) 一种含有双馈感应式风电机组电网的暂态稳定计算方法
CN109830987A (zh) 计及分布式光伏随机性的主动配电网概率稳定性分析方法
CN105958530A (zh) 一种具有无功自动补偿的微电网系统
Li et al. Structure preserving aggregation method for doubly-fed induction generators in wind power conversion
CN110336299B (zh) 一种考虑综合能源系统小干扰稳定的配网重构方法
CN117375068A (zh) 一种基于非线性仿射的风电场电压状态计算和优化方法及计算机可读介质
CN109740906B (zh) 一种含双馈风力发电机的配电网短路电流计算方法
CN113872189B (zh) Vsc接入薄弱电网时低频振荡特性的等效pll分析方法
CN112018787B (zh) 分布式电源的三相不平衡时序仿真方法
CN109815638B (zh) 一种结合模型切换和变步长的双馈风电模型仿真方法
CN112952863A (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200505

Termination date: 20210929