CN108365629A - 一种基于直算法的电力系统机电暂态仿真方法 - Google Patents

一种基于直算法的电力系统机电暂态仿真方法 Download PDF

Info

Publication number
CN108365629A
CN108365629A CN201810219284.1A CN201810219284A CN108365629A CN 108365629 A CN108365629 A CN 108365629A CN 201810219284 A CN201810219284 A CN 201810219284A CN 108365629 A CN108365629 A CN 108365629A
Authority
CN
China
Prior art keywords
generator
sqi
sdi
frequency
power
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
Application number
CN201810219284.1A
Other languages
English (en)
Other versions
CN108365629B (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201810219284.1A priority Critical patent/CN108365629B/zh
Publication of CN108365629A publication Critical patent/CN108365629A/zh
Priority to US16/969,583 priority patent/US20210057912A1/en
Priority to PCT/CN2019/077729 priority patent/WO2019174555A1/zh
Application granted granted Critical
Publication of CN108365629B publication Critical patent/CN108365629B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
    • H02J3/00125—Transmission line or load transient problems, e.g. overvoltage, resonance or self-excitation of inductive loads
    • 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • 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
    • 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
    • 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/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • 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
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
    • Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开了一种基于直算法的电力系统机电暂态仿真方法，包括以下步骤：A初始化电力系统参数；B采用直算法计算电力系统潮流分布；C计算每台发电机下一帧角加速度、角速度、转子角度、频率；D进行频率调节；E进行励磁调节；F根据发电机下一帧转子角度和频率计算发电机旋转电动势；G根据每台发电机的频率fi,n+1，计算全网的平均频率fW,n+1；H将平均频率fW,n+1带入每一个节点并计算其相应的电抗和电纳，并将旋转电动势代入到相应的发电机中，最后计算得出所有节点的新矩阵；I重新进入步骤B。本发明克服了传统仿真法的缺陷，无迭代、计算速度快、精度高且误差小，真实地反映了电网的变化特性。

Description

一种基于直算法的电力系统机电暂态仿真方法
技术领域
本发明属于电力系统仿真技术领域，具体涉及基于直算法的电力系统机电暂态仿真方法。
背景技术
电力系统机电暂态仿真是电力系统非常重要的分析方法。对运行中的电力系统，通过仿真可以预知在大扰动下(如短路故障、切除线路、发电机、负载、发电机调节励磁、冲击性负载以及变压器调档位等)会不会危及电力系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(如线路、变压器等)是否会出现过负载，以及可能出现过负载时应事先采取哪些预防措施。对规划中的电力系统，也可通过仿真来检验所提出的电力系统规划方案能否满足各种运行方式的要求。
传统电力系统机电暂态仿真的算法是联立求解电力系统微分方程组和代数方程组，以获得物理量的时域解。微分方程组的求解方法主要有隐式梯形积分法、改进尤拉法和龙格-库塔法等。代数方程组的求解法主要采用适用于求解非线性代数方程组的牛顿法，即迭代法。
传统机电暂态仿真方法的缺陷是计算结果误差大、精度低以及速度慢、还时常不收敛；且其建立在所有发电机的频率是同步一致的基础之上，这与电网的实际情况不符，难以反映电网的真实变化。
发明内容
为了解决现有技术存在的上述问题，本发明目的在于提供一种基于直算法的电力系统机电暂态仿真方法，其克服了传统仿真法的缺陷，无迭代、计算速度快、精度高且误差小，真实地反映了电网的变化特性，如各线路、负载和变压器的阻抗随着电网的频率变化而变化，电网中各发电机的频率也能按各自的规律动态地变化。
本发明所采用的技术方案为：一种基于直算法的电力系统机电暂态仿真方法，包括以下步骤：
A、初始化电力系统参数，计算电力系统中各节点初始矩阵；
B、采用直算法计算电力系统潮流分布，其中第i节点发电机第n帧的输出有功功率为Pi,n、输出无功功率为Qi,n、发电机端电压为Ui,n和功率因数为COSi,n；
C、进行频率调节，调节第i节点发电机第n+1帧动力矩Ti,n+1；
D、进行励磁调节，调节第i节点发电机第n+1帧励磁电流ILi,n+1；
E、计算第i节点发电机的第n+1帧角加速度角速度ωi,n+1、转子角度θi,n+1和频率fi,n+1：
其中，Ji是第i节点发电机的转动惯量，ΔT是帧计算时间间隔，n是帧序号，n的初始值为0；
F、根据发电机第n+1帧的转子角度θi,n+1、频率fi,n+1以及励磁电流ILi,n+1，计算发电机的第n+1帧旋转电动势：
发电机旋转电动势的绝对值|Ei,n+1|＝KLi×ILi,n+1×fi,n+1，则旋转电动势向量值其中KLi是i节点发电机的电动势系数；
G、用电网上所有发电机第n+1帧的平均频率代替电网的第n+1帧频率，其中m为电力系统中发电机的总数量；
H、根据fW,n+1计算各节点的电抗和电纳，将以及计算得到的各节点的电抗和电纳带入到各节点的初始矩阵中，最后计算得出各节点的新矩阵；
I、重新进入步骤B。
进一步地，步骤A中初始化电力系统参数具体包括如下过程：
A1、设置帧计算时间间隔ΔT；
A2、设置每台发电机的初始频率fi,0＝50Hz，初始角加速度初始角度θi,0＝0，初始角速度ωi,0＝2×π×fi,0；
A3、设置每台发电机的励磁电流ILi,0为额定值以及励磁系数KLi，则其旋转电动势初始值
A4、设网上系统频率初始值为fW,0＝50Hz,根据fW,0确定每一个节点的电抗和电纳，最后计算出所有节点的初始矩阵；
A5、设置并网发电机的动力矩，或按并网发电机的容量占比分摊全网负载，确定每台发电机的第1帧和起始帧动力矩
A6、设置调频发电机的调频系数KTi和死区频率Δfsqi；
A7、设置电压调节发电机的电压调节系数Kui、设置电压Usdi和死区电压ΔUsqi；
A8、设置无功调节发电机的无功调节系数KQi、设置无功功率Qsdi和死区无功功率ΔQsqi；
A9、设置功率因数调节发电机的功率因数调节系数Kcosi、设置功率因数COSsdi和死区功率因数ΔCOSsqi。
进一步地，步骤A4具体包括以下过程：设电力系统中，在初始频率fW,0＝50Hz时，负载的电阻为Ri,0和电抗为Xi,0；线路每公里电阻为ri,0、每公里电抗为xi,0、每公里电导为gi,0、每公里电纳为bi,0和线路长度为li；变压器的电导为Gti,0、电纳为Bti,0、电阻为Rti,0、电抗为Xti,0，原边匝数为ni,1和副边匝数为ni,2；发电机的内阻为r′i，0和电抗为x′i，0；则各节点初始矩阵如下：
负载初始矩阵为：
线路初始矩阵为：
其中，zi,0＝ri,0+jxi,0，yi,0＝gi,0+jbi,0，
变压器初始矩阵为：
发电机初始矩阵为：
进一步地，步骤H中，根据fW,n+1计算各节点的电抗和电纳的具体过程如下：
负载在频率fW,n+1下的电抗为
线路在频率fW,n+1下的每公里电抗为每公里电纳为
变压器在频率fW,n+1下的电抗为电纳为
发电机在频率fW,n+1下的电抗为
进一步地，步骤H中，各节点在频率fW,n+1下的新矩阵如下：
负载新矩阵为：
线路新矩阵为：
其中，zi,n+1＝ri,0+jxi,n+1，yi,n+1＝gi,0+jbi,n+1，
变压器新矩阵为：
发电机新矩阵为：
进一步地，所述步骤C中进行频率调节时，若发电机是非调频发电机，则有：Ti,n+1＝Ti,n；若发电机是调频发电机，则有：
当fi,n+1>50+Δfsqi时，Ti,n+1＝Ti,n-KTi×[fi,n+1-(50+Δfsqi)]；
当fi,n+1＜50-Δfsqi时，Ti,n+1＝Ti,n+KTi×[(50-Δfsqi)-fi,n+1]；
当fi,n+1≥50-Δfsqi且fi,n+1≤50+Δfsqi时，Ti,n+1＝Ti,n；
其中，KTi是调频系数，Δfsqi是死区频率。
进一步地，所述步骤D中进行励磁调节时，发电机的励磁调节只能在不调节、电压调节、无功调节和功率因数调节这四种之中选择一种：
若发电机不参与励磁调节，则ILi,n+1＝ILi,n；
若是电压调节发电机，则有：
当Ui,n>Usdi+ΔUsqi时，ILi,n+1＝ILi,n-Kui×[Ui,n-(Usdi+ΔUsqi)]；
当Ui,n＜Usdi-ΔUsqi时，ILi,n+1＝ILi,n+Kui×[(Usdi-ΔUsqi)-Ui,n]；
当Ui,n≥Usdi-ΔUsqi且Ui,n≤Usdi+ΔUsqi时，ILi,n+1＝ILi,n；
若是无功调节发电机，则有：
当Qi,n>Qsdi+ΔQsqi时，ILi,n+1＝ILi,n-KQi×[Qi,n-(Qsdi+ΔQsqi)]；
当Qi,n＜Qsdi-ΔQsqi时，ILi,n+1＝ILi,n+KQi×[(Qsdi-ΔQsqi)-Qi,n]；
当Qi,n≥Qsdi-ΔQsqi且Qi,n≤Qsdi+ΔQsqi时，ILi,n+1＝ILi,n；
若是功率因数调节发电机，则有：
当COSi,n>COSsdi+ΔCOSsqi时，ILi,n+1＝ILi,n+KCOSi×[COSi,n-(COSsdi+ΔCOSsqi)]；
当COSi,n＜COSsdi-ΔCOSsqi时，ILi,n+1＝ILi,n-KCOSi×[(COSsdi-ΔCOSsqi)-COSi,n]；
当COSi,n≥COSsdi-ΔCOSsqi且COSi,n≤COSsdi+ΔCOSsqi时，ILi,n+1＝ILi,n；
其中，Kui是发电机的电压调节系数，Usdi是设定电压，ΔUsqi是死区电压，Ui,n是第i节点发电机第n帧端口电压，KQi是发电机的无功调节系数，Qsdi是设定无功功率，ΔQsqi是死区无功功率，Qi,n是第i节点发电机第n帧输出无功功率，KCOSi是发电机的功率因数调节系数，COSsdi是设定功率因数，ΔCOSsqi是死区功率因数，COSi,n是第i节点发电机第n帧功率因数。
本发明的有益效果为：
(1)本仿真方法以帧为计算单位，一帧一帧不断重复地计算，此计算方法没有终结，每一帧都有输出结果。在数学上当所有发电机的角加速度为零时，即认为潮流进入稳态；在工程应用方面，当前后两帧计算的结果相差较小时(电压和电流的辐角除外)，可认为潮流进入稳态。
(2)每帧计算的值即为机电暂态过程的中间值。
(3)采用这种计算方式，可在计算时修改参数，模仿电网中的各种扰动，使潮流发生变化。如修改负载的Y值，可仿真负载变化对电网的影响；修改变压器的ni,1、ni,2,可仿真变压器的有载调压；修改发电机的励磁电流，可仿真发电机励磁改变后的电网变化过程；修改发电机的动力矩，可仿真发电机出力变化对电网的影响；若将线路分成两段，在中间插入负载，令负载的Y＝0，则是正常线路，若令负载Y＝10000，则可以仿真线路短路。
(4)本仿真法克服了传统仿真法的缺陷，无迭代、计算速度快、精度高且误差小，真实地反映了电网的变化特性，如各线路、负载和变压器的阻抗随着电网的频率变化而变化，电网中各发电机的频率也能按各自的规律动态地变化。
附图说明
图1是基于直算法的电力系统机电暂态仿真方法的流程图。
具体实施方式
下面结合附图及具体实施例对本发明做进一步阐释。
实施例
如图1所示，本发明所采用的技术方案为：一种基于直算法的电力系统机电暂态仿真方法，包括以下步骤：
A、初始化电力系统参数，计算电力系统中各节点初始矩阵；
B、采用直算法计算电力系统潮流分布，其中第i节点发电机第n帧的输出有功功率为Pi,n、输出无功功率为Qi,n、发电机端电压为Ui,n和功率因数为COSi,n；直算法计算电力系统潮流分布为现有技术，请参考申请人的专利CN201410142938.7以及专利申请CN201610783305.3；
C、进行频率调节，调节第i节点发电机第n+1帧动力矩Ti,n+1；
D、进行励磁调节，调节第i节点发电机第n+1帧励磁电流ILi,n+1；
E、计算第i节点发电机的第n+1帧角加速度aωi,n+1、角速度ωi,n+1、转子角度θi,n+1和频率fi,n+1：
其中，Ji是第i节点发电机的转动惯量，ΔT是帧计算时间间隔，n是帧序号，n的初始值为0；
F、根据发电机第n+1帧的转子角度θi,n+1、频率fi,n+1以及励磁电流ILi,n+1，计算发电机的第n+1帧旋转电动势：
发电机旋转电动势的绝对值|Ei，n+1|＝KLi×ILi，n+1×fi，n+1，则旋转电动势向量值其中KLi是i节点发电机的电动势系数；
G、用电网上所有发电机第n+1帧的平均频率(其中m为电力系统中发电机的总数量)代替电网的第n+1帧频率；
H、根据fW,n+1计算各节点的电抗和电纳，将以及计算得到的各节点的电抗和电纳带入到各节点的初始矩阵中，最后计算得出各节点的新矩阵；
I、重新进入步骤B。
在另一实施例中，步骤A中初始化电力系统参数具体包括如下过程：
A1、设置帧计算时间间隔ΔT；
A2、设置每台发电机的初始频率fi,0＝50Hz，初始角加速度aωi,0＝0，初始角度θi,0＝0，初始角速度ωi,0＝2×π×fi,0；
A3、设置每台发电机的励磁电流ILi,0为额定值以及励磁系数KLi，则其旋转电动势
A4、设网上系统频率初始值为fW,0＝50Hz,根据fW,0确定每一个节点的电抗和电纳，最后计算得出所有节点的初始矩阵；
A5、设置并网发电机的动力矩，或按并网发电机的容量占比分摊全网负载，确定每台发电机的第1帧和起始帧动力矩
A6、设置调频发电机的调频系数KTi和死区频率Δfsqi；
A7、设置电压调节发电机的电压调节系数Kui、设置电压Usdi和死区电压ΔUsqi；
A8、设置无功调节发电机的无功调节系数KQi、设置无功功率Qsdi和死区无功功率ΔQsqi；
A9、设置功率因数调节发电机的功率因数调节系数Kcosi、设置功率因数COSsdi和死区功率因数ΔCOSsqi。
在另一实施例中，步骤A4具体包括以下过程：设电力系统中，在初始频率fW,0＝50Hz时，负载的电阻为Ri,0和电抗为Xi,0；线路每公里电阻为ri,0、每公里电抗为xi,0、每公里电导为gi,0、每公里电纳为bi,0和线路长度为li；变压器的电导为Gti,0、电纳为Bti,0、电阻为Rti,0、电抗为Xti,0、原边匝数为ni,1和副边匝数为ni,2；发电机的内阻为r′i,0和电抗为x′i,0；则各节点初始矩阵如下：
负载初始矩阵为：
线路初始矩阵为：
其中，zi,0＝ri,0+jxi,0，yi,0＝gi,0+jbi,0，
变压器初始矩阵为：
发电机初始矩阵为：
在初始频率fW,0＝50Hz时，Ri,0、Xi,0、ri,0、xi,0、gi,0、bi,0、Gti,0、Bti,0、Rti,0、Xti,0、r′i,0和x′i,0的值可确定，为已知条件，将具体值带入相应初始矩阵。
在另一实施例中，步骤H中，根据fW,n+1计算各节点的电抗和电纳的具体过程如下：
负载在频率fW,n+1下的电抗为
线路在频率fW,n+1下的每公里电抗为每公里电纳为
变压器在频率fW,n+1下的电抗为电纳为发电机在频率fW,n+1下的电抗为
在另一实施例中，步骤H中，各节点在频率fW,n+1下的新矩阵如下：
负载新矩阵为：
线路新矩阵为：其中：
变压器新矩阵为：
发电机新矩阵为：
在另一实施例中，所述步骤C中进行频率调节时，
若发电机是非调频发电机，则Ti,n+1＝Ti,n；
若发电机是调频发电机，则有：
当fi,n+1>50+Δfsqi时，Ti,n+1＝Ti,n-KTi×[fi,n+1-(50+Δfsqi)]；
当fi,n+1＜50-Δfsqi时，Ti,n+1＝Ti,n+KTi×[(50-Δfsqi)-fi,n+1]；
当fi,n+1≥50-Δfsqi且fi,n+1≤50+Δfsqi时，Ti,n+1＝Ti,n；
其中，KTi是调频系数，Δfsqi是死区频率。
在另一实施例中，所述步骤D中进行励磁调节时，发电机的励磁调节只能在不调节、电压调节、无功调节和功率因数调节这四种之中选择一种：
若发电机不参与励磁调节，则ILi,n+1＝ILi,n；
若是电压调节发电机，则有：
当Ui,n>Usdi+ΔUsqi时，ILi,n+1＝ILi,n-Kui×[Ui,n-(Usdi+ΔUsqi)]；
当Ui,n＜Usdi-ΔUsqi时，ILi,n+1＝ILi,n+Kui×[(Usdi-ΔUsqi)-Ui,n]；
当Ui,n≥Usdi-ΔUsqi且Ui,n≤Usdi+ΔUsqi时，ILi,n+1＝ILi,n；
若是无功调节发电机，则有：
当Qi,n>Qsdi+ΔQsqi时，ILi,n+1＝ILi,n-KQi×[Qi,n-(Qsdi+ΔQsqi)]；
当Qi,n＜Qsdi-ΔQsqi时，ILi,n+1＝ILi,n+KQi×[(Qsdi-ΔQsqi)-Qi,n]；
当Qi,n≥Qsdi-ΔQsqi且Qi,n≤Qsdi+ΔQsqi时，ILi,n+1＝ILi,n；
若是功率因数调节发电机，则有：
当COSi,n>COSsdi+ΔCOSsqi时，ILi,n+1＝ILi,n+KCOSi×[COSi,n-(COSsdi+ΔCOSsqi)]；
当COSi,n＜COSsdi-ΔCOSsqi时，ILi,n+1＝ILi,n-KCOSi×[(COSsdi-ΔCOSsqi)-COSi,n]；
当COSi,n≥COSsdi-ΔCOSsqi且COSi,n≤COSsdi+ΔCOSsqi时，ILi,n+1＝ILi,n；
其中，Kui是发电机的电压调节系数，Usdi是设定电压，ΔUsqi是死区电压，Ui,n是第i节点发电机第n帧端口电压，KQi是发电机的无功调节系数，Qsdi是设定无功功率，ΔQsqi是死区无功功率，Qi,n是第i节点发电机第n帧输出无功功率，KCOSi是发电机的功率因数调节系数，COSsdi是设定功率因数，ΔCOSsqi是死区功率因数，COSi,n是第i节点发电机第n帧功率因数。
本发明不局限于上述可选的实施方式，任何人在本发明的启示下都可得出其他各种形式的产品。上述具体实施方式不应理解成对本发明的保护范围的限制，本发明的保护范围应当以权利要求书中界定的为准，并且说明书可以用于解释权利要求书。

Claims (7)

1.一种基于直算法的电力系统机电暂态仿真方法，其特征在于，包括以下步骤：
A、初始化电力系统参数，计算电力系统中各节点初始矩阵；
B、采用直算法计算电力系统潮流分布，其中第i节点发电机第n帧的输出有功功率为Pi,n、输出无功功率为Qi,n、发电机端电压为Ui,n和功率因数为COSi,n；
C、进行频率调节，调节第i节点发电机第n+1帧动力矩Ti,n+1；
D、进行励磁调节，调节第i节点发电机第n+1帧励磁电流ILi,n+1；
E、计算第i节点发电机的第n+1帧角加速度角速度ωi,n+1、转子角度θi,n+1和频率fi,n+1：
其中，Ji是第i节点发电机的转动惯量，ΔT是帧计算时间间隔，n是帧序号，n的初始值为0；
F、根据发电机第n+1帧的转子角度θi,n+1、频率fi,n+1以及励磁电流ILi,n+1，计算发电机的第n+1帧旋转电动势：
发电机旋转电动势的绝对值|Ei,n+1|＝KLi×ILi,n+1×fi,n+1，则旋转电动势向量值其中KLi是i节点发电机的电动势系数；
G、用电网上所有发电机第n+1帧的平均频率代替电网的第n+1帧频率，其中m为电力系统中发电机的总数量；
H、根据fW,n+1计算各节点的电抗和电纳，将以及计算得到的各节点的电抗和电纳带入到各节点的初始矩阵中，最后计算得出各节点的新矩阵；
I、重新进入步骤B。
2.根据权利要求1所述的基于直算法的电力系统机电暂态仿真方法，其特征在于，步骤A中初始化电力系统参数具体包括如下过程：
A1、设置帧计算时间间隔ΔT；
A2、设置每台发电机的初始频率fi,0＝50Hz，初始角加速度初始角度θi,0＝0，初始角速度ωi,0＝2×π×fi,0；
A3、设置每台发电机的励磁电流ILi,0为额定值以及励磁系数KLi，则其旋转电动势初始值
A4、设网上系统频率初始值为fW,0＝50Hz,根据fW,0确定每一个节点的电抗和电纳，最后计算出所有节点的初始矩阵；
A5、设置并网发电机的动力矩，或按并网发电机的容量占比分摊全网负载，确定每台发电机的第1帧和起始帧动力矩
A6、设置调频发电机的调频系数KTi和死区频率Δfsqi；
A7、设置电压调节发电机的电压调节系数Kui、设置电压Usdi和死区电压ΔUsqi；
A8、设置无功调节发电机的无功调节系数KQi、设置无功功率Qsdi和死区无功功率ΔQsqi；
A9、设置功率因数调节发电机的功率因数调节系数Kcosi、设置功率因数COSsdi和死区功率因数ΔCOSsqi。
3.根据权利要求2所述的基于直算法的电力系统机电暂态仿真方法，其特征在于，步骤A4具体包括以下过程：设电力系统中，在初始频率fW,0＝50Hz时，负载的电阻为Ri,0和电抗为Xi,0；线路每公里电阻为ri,0、每公里电抗为xi,0、每公里电导为gi,0、每公里电纳为bi,0和线路长度为li；变压器的电导为Gti,0、电纳为Bti,0、电阻为Rti,0、电抗为Xti,0，原边匝数为ni,1和副边匝数为ni,2；发电机的内阻为ri',0和电抗为xi',0；则各节点初始矩阵如下：
负载初始矩阵为：
线路初始矩阵为：
其中，zi,0＝ri,0+jxi,0，yi,0＝gi,0+jbi,0，
变压器初始矩阵为：
发电机初始矩阵为：
4.根据权利要求3所述的基于直算法的电力系统机电暂态仿真方法，其特征在于，步骤H中，根据fW,n+1计算各节点的电抗和电纳的具体过程如下：
负载在频率fW,n+1下的电抗为
线路在频率fW,n+1下的每公里电抗为每公里电纳为
变压器在频率fW,n+1下的电抗为电纳为
发电机在频率fW,n+1下的电抗为
5.根据权利要求4所述的基于直算法的电力系统机电暂态仿真方法，其特征在于，步骤H中，各节点在频率fW,n+1下的新矩阵如下：
负载新矩阵为：
线路新矩阵为：
其中，zi,n+1＝ri,0+jxi,n+1，yi,n+1＝gi,0+jbi,n+1，
变压器新矩阵为：
发电机新矩阵为：
6.根据权利要求2所述的基于直算法的电力系统机电暂态仿真方法，其特征在于，所述步骤C中进行频率调节时，
若发电机是非调频发电机，则有：Ti,n+1＝Ti,n；
若发电机是调频发电机，则有：
当fi,n+1>50+Δfsqi时，Ti,n+1＝Ti,n-KTi×[fi,n+1-(50+Δfsqi)]；
当fi,n+1＜50-Δfsqi时，Ti,n+1＝Ti,n+KTi×[(50-Δfsqi)-fi,n+1]；
当fi,n+1≥50-Δfsqi且fi,n+1≤50+Δfsqi时，Ti,n+1＝Ti,n；
其中，KTi是调频系数，Δfsqi是死区频率。
7.根据权利要求2所述的基于直算法的电力系统机电暂态仿真方法，其特征在于，所述步骤D中进行励磁调节时，发电机的励磁调节只能在不调节、电压调节、无功调节和功率因数调节这四种之中选择一种：
若发电机不参与励磁调节，则ILi,n+1＝ILi,n；
若是电压调节发电机，则有：
当Ui,n>Usdi+ΔUsqi时，ILi,n+1＝ILi,n-Kui×[Ui,n-(Usdi+ΔUsqi)]；
当Ui,n＜Usdi-ΔUsqi时，ILi,n+1＝ILi,n+Kui×[(Usdi-ΔUsqi)-Ui,n]；
当Ui,n≥Usdi-ΔUsqi且Ui,n≤Usdi+ΔUsqi时，ILi,n+1＝ILi,n；
若是无功调节发电机，则有：
当Qi,n>Qsdi+ΔQsqi时，ILi,n+1＝ILi,n-KQi×[Qi,n-(Qsdi+ΔQsqi)]；
当Qi,n＜Qsdi-ΔQsqi时，ILi,n+1＝ILi,n+KQi×[(Qsdi-ΔQsqi)-Qi,n]；
当Qi,n≥Qsdi-ΔQsqi且Qi,n≤Qsdi+ΔQsqi时，ILi,n+1＝ILi,n；
若是功率因数调节发电机，则有：
当COSi,n>COSsdi+ΔCOSsqi时，ILi,n+1＝ILi,n+KCOSi×[COSi,n-(COSsdi+ΔCOSsqi)]；
当COSi,n＜COSsdi-ΔCOSsqi时，ILi,n+1＝ILi,n-KCOSi×[(COSsdi-ΔCOSsqi)-COSi,n]；
当COSi,n≥COSsdi-ΔCOSsqi且COSi,n≤COSsdi+ΔCOSsqi时，ILi,n+1＝ILi,n；
其中，Kui是发电机的电压调节系数，Usdi是设定电压，ΔUsqi是死区电压，Ui,n是第i节点发电机第n帧端口电压，KQi是发电机的无功调节系数，Qsdi是设定无功功率，ΔQsqi是死区无功功率，Qi,n是第i节点发电机第n帧输出无功功率，KCOSi是发电机的功率因数调节系数，COSsdi是设定功率因数，ΔCOSsqi是死区功率因数，COSi,n是第i节点发电机第n帧功率因数。
CN201810219284.1A 2018-03-16 2018-03-16 一种基于直算法的电力系统机电暂态仿真方法 Active CN108365629B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201810219284.1A CN108365629B (zh) 2018-03-16 2018-03-16 一种基于直算法的电力系统机电暂态仿真方法
US16/969,583 US20210057912A1 (en) 2018-03-16 2019-03-12 Electromechanical transient simulation method for power system based on direct linear algorithm
PCT/CN2019/077729 WO2019174555A1 (zh) 2018-03-16 2019-03-12 一种基于直算法的电力系统机电暂态仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810219284.1A CN108365629B (zh) 2018-03-16 2018-03-16 一种基于直算法的电力系统机电暂态仿真方法

Publications (2)

Publication Number Publication Date
CN108365629A true CN108365629A (zh) 2018-08-03
CN108365629B CN108365629B (zh) 2020-03-24

Family

ID=63000687

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810219284.1A Active CN108365629B (zh) 2018-03-16 2018-03-16 一种基于直算法的电力系统机电暂态仿真方法

Country Status (3)

Country Link
US (1) US20210057912A1 (zh)
CN (1) CN108365629B (zh)
WO (1) WO2019174555A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109063408A (zh) * 2018-10-31 2018-12-21 邓宏伟 一种基于直算法的电磁机电暂态仿真算法
WO2019174555A1 (zh) * 2018-03-16 2019-09-19 邓宏伟 一种基于直算法的电力系统机电暂态仿真方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113346482B (zh) * 2021-05-19 2022-10-18 电子科技大学 基于sfr模型预测广域电力系统频率时空分布的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102593820A (zh) * 2011-12-22 2012-07-18 河海大学 考虑发电机励磁电流约束和电枢电流约束的连续潮流算法
CN103956741A (zh) * 2014-04-10 2014-07-30 邓宏伟 一字链及支链式的三相对称多电源非环网电力系统直算法
CN106026083A (zh) * 2016-06-27 2016-10-12 吉林大学 一种基于matlab的电网暂态稳定分析方法
US20160329714A1 (en) * 2015-05-08 2016-11-10 The Board Of Trustees Of The University Of Alabama Systems and methods for providing vector control of a grid connected converter with a resonant circuit grid filter
CN106451456A (zh) * 2016-08-30 2017-02-22 邓宏伟 一种基于环网电力系统的直算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108365629B (zh) * 2018-03-16 2020-03-24 邓宏伟 一种基于直算法的电力系统机电暂态仿真方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102593820A (zh) * 2011-12-22 2012-07-18 河海大学 考虑发电机励磁电流约束和电枢电流约束的连续潮流算法
CN103956741A (zh) * 2014-04-10 2014-07-30 邓宏伟 一字链及支链式的三相对称多电源非环网电力系统直算法
US20160329714A1 (en) * 2015-05-08 2016-11-10 The Board Of Trustees Of The University Of Alabama Systems and methods for providing vector control of a grid connected converter with a resonant circuit grid filter
CN106026083A (zh) * 2016-06-27 2016-10-12 吉林大学 一种基于matlab的电网暂态稳定分析方法
CN106451456A (zh) * 2016-08-30 2017-02-22 邓宏伟 一种基于环网电力系统的直算方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019174555A1 (zh) * 2018-03-16 2019-09-19 邓宏伟 一种基于直算法的电力系统机电暂态仿真方法
CN109063408A (zh) * 2018-10-31 2018-12-21 邓宏伟 一种基于直算法的电磁机电暂态仿真算法
WO2020088245A1 (zh) * 2018-10-31 2020-05-07 邓宏伟 一种基于直算法的电磁机电暂态仿真算法

Also Published As

Publication number Publication date
US20210057912A1 (en) 2021-02-25
CN108365629B (zh) 2020-03-24
WO2019174555A1 (zh) 2019-09-19

Similar Documents

Publication Publication Date Title
CN107086589B (zh) 一种基于临界参数抑制水电机组超低频振荡的调速器参数优化方法
CN104866665B (zh) 基于接口等值与交互的含电力电子设备的混合仿真方法
EP2529462B1 (en) Method for emulation of synchronous machine
CN108365629A (zh) 一种基于直算法的电力系统机电暂态仿真方法
CN102682358B (zh) 一种评估风电并网规模与电网网架适应性的规划仿真方法
Zhang et al. Optimized virtual inertia of wind turbine for rotor angle stability in interconnected power systems
CN107453375A (zh) 巨型水电机组功率模式下的一次调频电量补偿精细化控制方法及装置
CN103545810A (zh) 基于小信号稳定分析的微电网逆变器下垂自动控制方法
CN106549402B (zh) Tcsc抑制交直流并联运行系统次同步振荡的方法
Mahvash et al. A new approach for power quality improvement of DFIG based wind farms connected to weak utility grid
CN102904246B (zh) 一种点对网输电系统暂态稳定输电极限的快速计算方法
Tabatabaei et al. Designing power system stabilizer with PID controller
CN104808147B (zh) 一种三相短路下双馈发电机仿真方法及其系统
CN108258707A (zh) 一种具有调峰调频功能的中压直挂式储能换流器系统及其控制方法
CN108493966A (zh) 一种基于虚拟同步技术的微电网不平衡负荷控制方法和装置
CN109494808A (zh) 一种虚拟同步发电机的离并网平滑切换控制方法
CN105281329B (zh) 一种提高2机电力系统暂态频率稳定性的单相重合时序整定方法
CN106532743A (zh) 一种基于振荡能量的互联系统低频振荡安控量获取方法
Zacharia et al. Wide area control of governors and power system stabilizers with an adaptive tuning of coordination signals
CN105224812A (zh) 一种负荷模型中的静态负荷频率因子聚合方法
Bignucolo et al. Advantages of the virtual synchronous machine regulation for integrating low-inertia variable renewable generation in transmission systems
CN106505582A (zh) 一种基于网络预测方式的动态无功功率‑电压协同控制方法
CN116402017A (zh) 一种用于电力系统频率动态分析的双馈风机简化建模方法
CN106300353B (zh) 由转差引起的双馈风机定子间谐波电流解析模型建立方法
CN102957166B (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