CN117728407B - 新能源逆变器异构并网暂态稳定问题二维相图分析方法 - Google Patents
新能源逆变器异构并网暂态稳定问题二维相图分析方法 Download PDFInfo
- Publication number
- CN117728407B CN117728407B CN202410179651.5A CN202410179651A CN117728407B CN 117728407 B CN117728407 B CN 117728407B CN 202410179651 A CN202410179651 A CN 202410179651A CN 117728407 B CN117728407 B CN 117728407B
- Authority
- CN
- China
- Prior art keywords
- grid
- inverter
- power
- heterogeneous
- delta
- 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.)
- Active
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 51
- 238000010587 phase diagram Methods 0.000 title claims abstract description 17
- 238000004458 analytical method Methods 0.000 title claims abstract description 13
- 238000010586 diagram Methods 0.000 claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims abstract description 19
- 230000007246 mechanism Effects 0.000 claims abstract description 11
- 238000013507 mapping Methods 0.000 claims abstract description 9
- 102100020865 EKC/KEOPS complex subunit LAGE3 Human genes 0.000 claims description 25
- 101001137983 Homo sapiens EKC/KEOPS complex subunit LAGE3 Proteins 0.000 claims description 25
- MOJZMWJRUKIQGL-XILRTYJMSA-N procyanidin C1 Chemical compound C1([C@@H]2[C@H](O)[C@H](C3=C(O)C=C(O)C=C3O2)C2=C3O[C@@H]([C@H](O)[C@H](C3=C(O)C=C2O)C=2C(O)=CC(O)=C3C[C@H]([C@H](OC3=2)C=2C=C(O)C(O)=CC=2)O)C=2C=C(O)C(O)=CC=2)=CC=C(O)C(O)=C1 MOJZMWJRUKIQGL-XILRTYJMSA-N 0.000 claims description 25
- 101100122529 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) GON7 gene Proteins 0.000 claims description 12
- 230000001360 synchronised effect Effects 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 9
- 238000002347 injection Methods 0.000 claims description 8
- 239000007924 injection Substances 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000013178 mathematical model Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 239000000243 solution Substances 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001687 destabilization Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Landscapes
- Inverter Devices (AREA)
Abstract
本发明涉及新能源电力系统技术领域，公开了一种新能源逆变器异构并网暂态稳定问题二维相图分析方法，根据异构并网系统的结构和参数，分析无穷大电网故障的主要特征，得到暂态故障工况；基于所述暂态故障工况对所述异构并网系统进行准静态建模；基于保结构降维映射，利用等高线将三维空间中的构网功率曲面与跟网q轴电压曲面降至二维平面，得到异构并网系统的关键区域划分图；在关键区域划分图中，绘制异构并网系统的二维相轨迹，根据各区域的功角物理特性，洞察系统在所述暂态故障工况下的物理动态过程，得到异构并网系统的暂态失稳机理。为新能源逆变器异构并网系统暂态稳定相关的优化与控制提供了参考。
Description
技术领域
本发明涉及新能源电力系统技术领域，具体为一种新能源逆变器异构并网暂态稳定问题二维相图分析方法。
背景技术
电力电子逆变器是风、光等新能源接入电网的关键接口。目前新能源并网逆变器主要分为跟网型和构网型两种类型，随着新型电力系统建设速度的加快，在未来的新能源场站中，跟网型和构网型逆变器混合并联于主网形成异构系统将成为一种常态。当电网电压发生严重暂降时，异构系统中的两类逆变器之间会发生非常复杂的相互作用，其动态特性可用一组高阶非线性微分方程组来描述，是典型的复杂高阶非线性问题。由于该系统在此过程中的强非线性与强交互耦合特性，现有的传统相图分析法、Lyapunov直接法等分析方法都无法清晰地洞察暂态过程中整个系统的物理动态特性，无法直观地解释其暂态失稳机理。因此，需要明确电网电压暂降后系统中跟网、构网型逆变器功角的物理动态过程，更加清晰地阐述异构并网系统的暂态失稳机理，为实现系统运行的优化设计提供理论支撑。
发明内容
针对背景技术中所提出的问题，本发明的目的在于提供一种新能源逆变器异构并网暂态稳定问题二维相图分析方法，可在电网电压严重暂降时，清晰地洞察跟网型、构网型逆变器功角的物理动态过程，更直观地揭示整个系统的暂态失稳机理。为新能源基地并网运行时，各逆变器的协调运行与致稳控制提供参考依据，为新能源基地调度运行和稳定优化提供理论支撑，提升系统稳定性。技术方案如下：
新能源逆变器异构并网暂态稳定问题二维相图分析方法，包括以下步骤：
步骤1：将构网型逆变器与跟网型逆变器混合并联于无穷大电网中形成异构并网系统，并定义异构并网系统的结构和参数；
步骤2：分析无穷大电网故障的特征，得到暂态故障工况的典型特性，并基于所述暂态故障工况对所述异构并网系统进行准静态建模；
步骤3：基于结构降维映射的思想，利用等高线将三维空间中的构网功率曲面与跟网q轴电压曲面降至二维平面，得到异构并网系统的关键区域划分图；
步骤4：在关键区域划分图中，绘制异构并网系统的二维相轨迹，洞察系统在所述暂态故障工况下的物理动态过程，得到异构并网系统的暂态失稳机理。
进一步的，所述步骤1中定义异构并网系统的结构和参数具体包括：
定义跟网型逆变器的并网点为PCC1，并网点PCC1经等效电阻为R1和等效电感为L1的集电线路连接至汇流母线Bus-C，且该集电线路的等效导纳为；
定义构网型逆变器的并网点为PCC2，并网点PCC2经等效电阻为R2和等效电感为L2的集电线路连接至汇流母线Bus-C，且该集电线路的等效导纳为所述汇流母线Bus-C经等效电阻为R3和等效电感为L3的输电线路与无穷大电网相连，且输电线路对应的等效导纳为/>
定义无穷大电网的电压为E∠0°，注入汇流母线Bus-C的等效电流为汇流母线Bus-C的电压为Vc∠θc；/>为电网注入汇流母线电流的相角；θc为汇流母线Bus-C的电压相角；
跟网型逆变器并网点PCC1的电压为VPCC1∠δPCC1，输出的电流为输出有功功率和无功功率分别为P1和Q1；其中，δPCC1为跟网型逆变器并网点PCC1的实际电压相角；δ1为跟网型逆变器的功角；/>为跟网型逆变器注入电流的期望功率因数角；
构网型逆变器并网点PCC2的电压为Vref∠δ2，输出的电流为输出有功功率和无功功率分别为P2和Q2，其中，δ2为构网型逆变器的功角，/>为构网型逆变器注入汇流母线电流的相角。
更进一步的，所述步骤2中暂态故障工况为：三相对称的电网电压暂降，且故障只影响电网电压的幅值，电网电压不存在负序分量；
基于所述暂态故障工况对所述异构并网系统进行准静态建模具体为：
综合基尔霍夫电压、电流定律，对汇流母线节点列写基本电路方程得到：
公式中，θ2为构网型逆变器与汇流母线之间的线路导纳角；θ3为电网与汇流母线之间的线路导纳角；Y2和Y3分别为等效导纳和/>的模值，即/>
进而得跟网型逆变器并网点PCC1的电压为：
构网型逆变器注入汇流母线的电流为：
根据等式(2)计算得到锁相环同步参考坐标系下的跟网型逆变器并网点PCC1电压的q轴分量为：
其中，Y1为等效导纳的模值，θ1为跟网型逆变器与汇流母线之间的线路导纳角；Y10、Y11和Y12分别为复导纳/>和/>的模值；θ10、θ11和θ12分别为复导纳/>和/>的导纳角；
根据等式(3)计算得到构网型逆变器注入汇流母线的有功功率为：
公式中，Re表示取实部；V2为构网型逆变器并网点PCC2的电压相量；为构网型逆变器注入汇流母线电流相量的共轭；Yt为复导纳/>的模值；θt为复导纳/>的导纳角；
结合跟网型逆变器锁相环模型与构网型逆变器功率同步环模型，得到所述异构并网系统的准静态数学模型的微分代数方程组：
公式中，Ki-PLL为锁相环中的积分增益系数；Kp-PLL为锁相环中的比例增益系数；KP-PSC为构网型逆变器功率同步控制环的同步系数；Pref为构网型逆变器的有功功率参考值；
复导纳和/>分别表示为：
更进一步的，所述步骤3具体为：
步骤3.1：根据所述锁相环同步参考坐标系下的跟网型逆变器并网点PCC1电压的q轴分量的表达式，以跟网型逆变器功角δ1为x轴，构网型逆变器功角δ2为y轴，VPCC1q为z轴，绘制跟网型逆变器的q轴电压曲面与VPCC1q＝0平面，并做出其交线，当VPCC1q>0时，δ1的二阶导数大于0，跟网型逆变器功角δ1做加速运动，反之当VPCC1q<0时，跟网型逆变器功角δ1做减速运动；
步骤3.2：根据构网型逆变器注入汇流母线的有功功率表达式，以跟网型逆变器功角δ1为x轴，构网型逆变器功角δ2为y轴，P2为z轴，绘制构网型逆变器的有功功率曲面与P2＝Pref平面，并做出其交线，当P2<Pref时，δ2的一阶导数大于0，构网型逆变器功角δ2做增值运动，反之当P2>Pref时，构网型逆变器功角δ2做减值运动；
步骤3.3：利用降维映射思想，以跟网型逆变器和构网型逆变器的平衡交线为关键区域边界，将三维曲面降维至二维平面，同时保留各区域内功角δ1、δ2的关键导数信息，进而得到跟网型逆变器功角δ1的加速/减速区域与构网型逆变器功角δ2的增值/减值区域，得到异构并网系统在δ1-δ2平面的关键区域划分图。
更进一步的，所述步骤4具体包括：
求解所述异构并网系统的准静态数学模型的微分代数方程组，得到跟网型逆变器功角δ1和构网型逆变器功角δ2的一组数值解，基于该数值解在步骤3.3所得到的δ1-δ2平面的关键区域划分图上绘制关于δ1-δ2的二维相轨迹；
根据关键区域划分图中功角δ1的加速/减速区域与δ2的增值/减值区域，分析所绘制的二维相轨迹的运动方向与变化趋势，并以此为基础，分析导致系统失稳的可能因素，明确并网异构系统的暂态失稳机理。
本发明与现有技术相比，具有如下的优点和有益效果：
1、本发明将跟网型逆变器与构网型逆变器混合并联于无穷大电网形成异构并网系统，对电网故障的主要特征进行分析得到暂态故障工况，并基于暂态故障工况对异构并网系统进行准静态建模，可以得异构并网系统的关键平衡方程。
2、本发明利用降维映射将跟网q轴电压曲面与构网有功功率曲面由三维空间降至二维平面，结合关键平衡方程得到关键区域划分图，通过在区域划分图中绘制二维相轨迹，明确了异构并网系统在暂态故障工况下动态特性，实现对系统物理动态过程的清晰描述，得到了异构并网系统的暂态失稳机理，克服了现有技术的缺陷。
3、本发明在电网电压严重暂降时，可清晰地描述跟网型、构网型逆变器功角的物理动态过程，具有很强的物理洞察力，更直观地揭示整个系统的暂态失稳机理。为新能源基地并网运行时，各逆变器的协调运行与致稳控制提供参考依据，为新能源基地调度运行和稳定优化提供理论支撑，提升系统稳定性。
附图说明
图1为本发明实施例提供的跟网型逆变器与构网型逆变器异构并网系统示意图。
图2(a)为跟网型逆变器三维空间中的电压平衡关系示意图。
图2(b)为构网型逆变器三维空间中的功率平衡关系示意图。
图3为本发明实施例提供的异构并网系统的关键区域划分图。
图4为Pref＝11000W时异构系统的二维相图。
图5(a)为Pref＝11000W时跟网型逆变器功角曲线图。
图5(b)为Pref＝11000W时构网型逆变器功角曲线图。
图6为Pref＝14600W时异构系统的二维相图。
图7(a)为Pref＝14600W时跟网型逆变器功角曲线。
图7(b)为Pref＝14600W时构网型逆变器功角曲线。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
本发明提供了一种新能源逆变器异构并网暂态稳定问题二维相图分析方法，包括以下步骤：
步骤S1、将构网型逆变器与跟网型逆变器混合并联于无穷大电网形成异构并网系统，并定义所述异构并网系统的结构和参数。
图1给出了异构并网系统的拓扑图，其中参数分别为：L1＝1mH，L2＝10mH，L3＝1mH，正常运行时的电网电压E的幅值为311V，跟网逆变器采用纯有功注入，即I1＝100A，在一种可选的实施例中，定义所述混合并联系统的结构和参数包括：
跟网型逆变器的并网点为PCC1，该并网点经等效电感为L1的集电线路连接至汇流母线Bus-C，且该集电线路的等效导纳为
构网型逆变器的并网点为PCC2，该并网点经等效电感为L2的集电线路连接至汇流母线Bus-C，且该集电线路的等效导纳为
汇流母线Bus-C经等效电感为L3的输电线路与无穷大电网相连，且输电线路对应的等效导纳为
定义无穷大电网的电压为E∠0°，注入汇流母线的等效电流为汇流母线的电压为Vc∠θc；/>为电网注入汇流母线电流的相角；θc为汇流母线Bus-C的电压相角。
跟网型逆变器并网点PCC1的电压为VPCC1∠δPCC1，输出的电流为输出有功功率和无功功率分别为P1和Q1；δPCC1为跟网型逆变器并网点PCC1的实际电压相角；其中，δ1定义为跟网型逆变器的功角；/>为跟网型逆变器注入电流的期望功率因数角；
构网型逆变器的并网点PCC2的电压为Vref∠δ2，输出的电流为输出有功功率和无功功率分别为P2和Q2，其中，δ2定义为构网型逆变器的功角；/>为构网型逆变器注入汇流母线电流的相角。
需要说明的是，在本实施例中，跟网逆变器和构网逆变器的运行方式如下：
跟网逆变器通过锁相环获得其并网点PCC1的电压相位，并以此相位为参考，向汇流母线注入相对相角为幅值为I1的电流，其中/>又称为跟网电流和电压的期望相角差，锁相环中的比例增益系数为Kp-PLL，锁相环中的积分增益系数为Ki-PLL。
构网逆变器通过采集其并网点PCC2的有功功率P2和无功功率Q2，通过功率同步控制发出角频率为ω，电压幅值为常数Vref的电压Vref∠δ2，功率同步控制有功环中的同步系数为KP-PSC。
步骤S2、分析无穷大电网故障的主要特征，得到暂态故障工况，基于所述暂态故障工况对所述异构并网系统进行准静态建模；
在一种可选的实施例中，对主网故障的主要特征进行分析得到暂态故障工况包括：
故障类型为电网电压三相对称的电压暂降，且故障仅仅影响电压E的幅值；
故障电网电压不存在负序分量。
在一种可选的实施例中，基于所述暂态故障工况对所述异构并网系统进行准静态建模包括：
综合基尔霍夫电压、电流定律，对汇流母线节点列写基本电路方程得到：
综合等式(1)和等式(2)得到：
公式中，θ2为构网型逆变器与汇流母线之间的线路导纳角；θ3为电网与汇流母线之间的线路导纳角；Y2和Y3分别为等效导纳和/>的模值，即/>
进而可得跟网型逆变器并网点PCC1的电压为：
将等式(3)带入等式(1)可得构网型逆变器注入汇流母线的电流为：
根据等式(4)可计算得到锁相环同步参考坐标系下的跟网型逆变器并网点PCC1电压的q轴分量为:
其中，Y1为等效导纳的模值；θ1为跟网型逆变器与汇流母线之间的线路导纳角；Y10、Y11和Y12分别为复导纳/>和/>的模值；θ10、θ11、θ12分别为复导纳/>和/>的导纳角。
根据等式(5)可计算得到构网型逆变器注入汇流母线的有功功率为：
公式中，Re表示取实部；V2为构网型逆变器并网点PCC2的电压相量；为构网型逆变器注入汇流母线电流相量的共轭；Yt为复导纳/>的模值，θt为复导纳/>的导纳角。
结合跟网型逆变器锁相环模型与构网型逆变器功率同步环模型，可得到该异构并网系统的准静态数学模型为：
公式中，Ki-PLL为锁相环中的积分增益系数；Kp-PLL为锁相环中的比例增益系数；KP-PSC为构网型逆变器功率同步控制环的同步系数；Pref为构网型逆变器的有功功率参考值。
在本实施例中，等式(6)、等式(7)和等式(8)中的参数可以表示为：
上述准静态模型可以描绘跟网型、构网型逆变器并网异构系统在电网电压暂降时的暂态过程，其中，等式(8)中的代数方程描绘了两类逆变器平衡关系的静态特性，揭示了系统平衡点的存在性；等式(8)中的微分方程描述了各类逆变器的功角动态特性。
步骤S3、基于结构降维映射的思想，利用等高线将三维空间中的构网功率曲面与跟网q轴电压曲面降至二维平面，得到异构并网系统的关键区域划分图。
在一种可选的实施例中，采用降维映射做出异构并网系统的关键区域划分图，如图2(a)所示，根据等式(6)，以δ1为x轴，δ2为y轴，VPCC1q为z轴，绘制跟网型逆变器的q轴电压曲面与VPCC1q＝0平面，并做出其交线，当VPCC1q>0时δ1的二阶导数大于0，跟网型逆变器功角δ1做加速运动，反之当VPCC1q<0时，跟网型逆变器功角δ1做减速运动。
如图2(b)所示，根据等式(7)，以δ1为x轴，δ2为y轴，P2为z轴，绘制构网型逆变器的有功功率曲面与P2＝Pref平面，并做出其交线，当P2<Pref时，δ2的一阶导数大于0，构网型逆变器功角δ2做增值运动，反之当P2>Pref时，构网型逆变器功角δ2做减值运动。
利用降维映射思想，以跟网型、构网型逆变器的平衡交线为关键区域边界，将三维曲面降维至二维平面，可同时保留功角δ1、δ2的关键导数信息，进而可得跟网型逆变器的功角加速/减速区域与构网型逆变器的功角增值/减值区域，得到异构并网系统的关键区域划分图。图3示出了系统的关键区域划分图，根据关键区域划分图中所示的区域导数值，可以很直观地预测δ1-δ2平面内任一点的功角运动趋势。
步骤S4、在关键区域划分图中，绘制异构并网系统的二维相轨迹，洞察系统在所述暂态故障工况下的物理动态过程，得到异构并网系统的暂态失稳机理。
在一种可选的实施例中，利用二维相图分析异构并网系统的暂态过程包括：求解等式(8)所示的微分代数方程组，并在异构系统关键区域划分图中做出系统的二维相轨迹，得到系统暂态过程中的功角动态特性，根据相轨迹所处的区域对相轨迹地运动方向进行预测，进而清晰地洞察整个系统的物理动态过程。
本实施例还提供了一种适用于新能源逆变器异构并网系统暂态稳定问题的二维相图分析方法的实施实例。
在实例1中，故障时电网电压暂降为120V，跟网逆变器采用纯有功电流注入，即I1＝100A，，Pref＝11000W。根据给定参数，利用跟网型、构网型逆变器异构并网系统的准静态数学模型和二维相图分析方法，可得到异构系统在该工况下的二维相图(如图4所示)，二维相轨迹显示系统将从故障前的稳态工作点a过渡到稳态工作点b，且其过渡过程中的动态物理过程都可以根据功角区域划分图中各关键区域的特性被清晰地洞察，相轨迹反应系统最终可以保持暂态稳定。
图4为本发明实施例提供的电网电压E由311V暂降至120V，跟网型逆变器纯有功电流注入Id＝100A，Pref＝11000W时系统的二维相轨迹，此时有功参考的大小合适，构网平衡曲线与跟网平衡曲线存在交点，系统存在稳定的平衡点b，因此，异构系统的二维相轨迹将从故障前的工作点a过渡到故障后的稳态工作点b，且根据其相轨迹在关键区域划分图中的位置，可以清晰直观地分析相轨迹的走向。相轨迹反应系统最终可以保持暂态稳定。
图5(a)和图5(b)为本发明实施例提供的电网电压E由311V暂降至120V，跟网型逆变器纯有功电流注入Id＝100A，Pref＝11000W时，系统仿真得到此工况下异构系统暂态过程的功角曲线仿真结果，仿真结果表明，图4所示的二维相轨迹完全正确地预测了跟网型、构网型逆变器的暂态功角特性，清晰地洞察了异构并网系统的物理动态过程。
在实施例2中，故障时电网电压暂降为120V，跟网逆变器采用纯有功电流注入，即I1＝100A，Pref＝14600W。根据给定参数，利用跟网型、构网型逆变器异构并网系统的准静态数学模型和二维相图分析方法，可得到系统在该工况下的二维相图。
图6为本发明实施例提供的电网电压E由311V暂降至120V，跟网型逆变器纯有功电流注入Id＝100A，Pref＝14600W时系统的二维相轨迹，此时有功参考过大，导致构网平衡曲线与跟网平衡曲线在相轨迹可能经过的区间内丧失了交点，即系统在此工况下不存在稳定的平衡点，且根据其相轨迹在关键区域划分图中的位置，可以清晰直观地分析相轨迹的走向：相轨迹在几个震荡周期后被拉向跟网平衡曲线VPCC1q＝0，由于系统不存在稳定的平衡点，所以相轨迹将延跟网平衡曲线向δ2趋于无穷的方向运动，最终导致构网型逆变器的功角发散至无穷大，而跟网型逆变器的功角将在0.25-0.95rad之间做周期性震荡。
图7(a)和图7(b)为本发明实施例提供的电网电压E由311V暂降至120V，跟网型逆变器纯有功电流注入Id＝100A，Pref＝14600W时，系统仿真得到的此工况下异构系统暂态过程的功角曲线仿真结果，仿真结果表明，图6所示的二维相轨迹完全正确地预测了跟网型、构网型逆变器的暂态功角特性，清晰地洞察了异构并网系统的物理动态过程。
Claims (1)
1.新能源逆变器异构并网暂态稳定问题二维相图分析方法，其特征在于，包括以下步骤：
步骤1：将构网型逆变器与跟网型逆变器混合并联于无穷大电网中形成异构并网系统，并定义异构并网系统的结构和参数；
步骤2：分析无穷大电网故障的特征，得到暂态故障工况的典型特性，并基于所述暂态故障工况对所述异构并网系统进行准静态建模；
步骤3：基于结构降维映射的思想，利用等高线将三维空间中的构网功率曲面与跟网q轴电压曲面降至二维平面，得到异构并网系统的关键区域划分图；
步骤4：在关键区域划分图中，绘制异构并网系统的二维相轨迹，洞察系统在所述暂态故障工况下的物理动态过程，得到异构并网系统的暂态失稳机理；
所述步骤1中，定义异构并网系统的结构和参数具体包括：
定义跟网型逆变器的并网点为PCC1，并网点PCC1经等效电阻为R1和等效电感为L1的集电线路连接至汇流母线Bus-C，且该集电线路的等效导纳为定义构网型逆变器的并网点为PCC2，并网点PCC2经等效电阻为R2和等效电感为L2的集电线路连接至汇流母线Bus-C，且该集电线路的等效导纳为/>所述汇流母线Bus-C经等效电阻为R3和等效电感为L3的输电线路与无穷大电网相连，且输电线路对应的等效导纳为/>
定义无穷大电网的电压为E∠0°，注入汇流母线Bus-C的等效电流为汇流母线Bus-C的电压为Vc∠θc；/>为电网注入汇流母线电流的相角，θc为汇流母线Bus-C的电压相角；跟网型逆变器并网点PCC1的电压为VPCC1∠δPCC1，输出的电流为/>输出有功功率和无功功率分别为P1和Q1；其中，δPCC1为跟网型逆变器并网点PCC1的实际电压相角；δ1为跟网型逆变器的功角，/>为跟网型逆变器注入电流的期望功率因数角；
构网型逆变器并网点PCC2的电压为Vref∠δ2，输出的电流为输出有功功率和无功功率分别为P2和Q2，其中，δ2为构网型逆变器的功角，/>为构网型逆变器注入汇流母线电流的相角；
所述步骤2中，暂态故障工况为：三相对称的电网电压暂降，且故障只影响电网电压的幅值，电网电压不存在负序分量；
基于所述暂态故障工况对所述异构并网系统进行准静态建模具体为：
综合基尔霍夫电压、电流定律，对汇流母线节点列写基本电路方程得到：
公式中，θ2为构网型逆变器与汇流母线之间的线路导纳角，θ3为电网与汇流母线之间的线路导纳角；Y2和Y3分别为等效导纳和/>的模值，即/>
进而得跟网型逆变器并网点PCC1的电压为：
构网型逆变器注入汇流母线的电流为：
根据等式(2)计算得到锁相环同步参考坐标系下的跟网型逆变器并网点PCC1电压的q轴分量为：
其中，Y1为等效导纳的模值，θ1为跟网型逆变器与汇流母线之间的线路导纳角；Y10、Y11和Y12分别为复导纳/>和/>的模值；θ10、θ11和θ12分别为复导纳/>和/>的导纳角；根据等式(3)计算得到构网型逆变器注入汇流母线的有功功率为：
公式中，Re表示取实部；V2为构网型逆变器并网点PCC2的电压相量；为构网型逆变器注入汇流母线电流相量的共轭；Y为复导纳/>的模值；θt为复导纳/>的导纳角；
结合跟网型逆变器锁相环模型与构网型逆变器功率同步环模型，得到所述异构并网系统的准静态数学模型的微分代数方程组：
公式中，Ki-PLL为锁相环中的积分增益系数；Kp-PLL为锁相环中的比例增益系数；KP-PSC为构网型逆变器功率同步控制环的同步系数；Pref为构网型逆变器的有功功率参考值；
复导纳和/>分别表示为：
所述步骤3具体为：
步骤3.1：根据所述锁相环同步参考坐标系下的跟网型逆变器并网点PCC1电压的q轴分量的表达式，以跟网型逆变器功角δ1为x轴，构网型逆变器功角δ2为y轴，VPCC1q为z轴，绘制跟网型逆变器的q轴电压曲面与VPCC1q＝0平面，并做出其交线，当VPCC1q>0时，δ1的二阶导数大于0，跟网型逆变器功角δ1做加速运动，反之当VPCC1q<0时，跟网型逆变器功角δ1做减速运动；
步骤3.2：根据构网型逆变器注入汇流母线的有功功率表达式，以跟网型逆变器功角δ1为x轴，构网型逆变器功角δ2为y轴，P2为z轴，绘制构网型逆变器的有功功率曲面与P2＝Pref平面，并做出其交线，当P2<Pref时，δ2的一阶导数大于0，构网型逆变器功角δ2做增值运动，反之当P2>Pref时，构网型逆变器功角δ2做减值运动；
步骤3.3：利用降维映射思想，以跟网型逆变器和构网型逆变器的平衡交线为关键区域边界，将三维曲面降维至二维平面，同时保留各区域内功角δ1、δ2的关键导数信息，进而得到跟网型逆变器功角δ1的加速或减速区域与构网型逆变器功角δ2的增值或减值区域，得到异构并网系统在δ1-δ2平面的关键区域划分图；
所述步骤4具体包括：
求解所述异构并网系统的准静态数学模型的微分代数方程组，得到跟网型逆变器功角δ1和构网型逆变器功角δ2的一组数值解，基于该数值解在步骤3.3所得到的δ1-δ2平面的关键区域划分图上绘制关于δ1-δ2的二维相轨迹；
根据关键区域划分图中功角δ1的加速或减速区域与δ2的增值或减值区域，分析所绘制的二维相轨迹的运动方向与变化趋势，并以此为基础，分析导致系统失稳的可能因素，明确并网异构系统的暂态失稳机理。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410179651.5A CN117728407B (zh) | 2024-02-18 | 2024-02-18 | 新能源逆变器异构并网暂态稳定问题二维相图分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410179651.5A CN117728407B (zh) | 2024-02-18 | 2024-02-18 | 新能源逆变器异构并网暂态稳定问题二维相图分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117728407A CN117728407A (zh) | 2024-03-19 |
| CN117728407B true CN117728407B (zh) | 2024-05-03 |
Family
ID=90209284
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410179651.5A Active CN117728407B (zh) | 2024-02-18 | 2024-02-18 | 新能源逆变器异构并网暂态稳定问题二维相图分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117728407B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116094025A (zh) * | 2022-09-07 | 2023-05-09 | 南瑞集团有限公司 | 一种含跟网型和构网型变流器功率系统暂态稳定分析方法 |
| WO2023108925A1 (zh) * | 2021-12-15 | 2023-06-22 | 新疆金风科技股份有限公司 | 风力发电机组、风电变流器的控制方法和控制装置 |
| CN116961116A (zh) * | 2023-09-20 | 2023-10-27 | 四川大学 | 基于自适应q轴电压反馈的构网逆变器暂态稳定提升方法 |
| CN117498440A (zh) * | 2023-11-17 | 2024-02-02 | 国网四川省电力公司电力科学研究院 | 并网逆变器混合并联系统暂态电流边界的分析方法及装置 |
| CN117498433A (zh) * | 2023-11-03 | 2024-02-02 | 国网四川省电力公司电力科学研究院 | 一种混合并联系统暂态稳定功率极限量化方法 |
-
2024
- 2024-02-18 CN CN202410179651.5A patent/CN117728407B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023108925A1 (zh) * | 2021-12-15 | 2023-06-22 | 新疆金风科技股份有限公司 | 风力发电机组、风电变流器的控制方法和控制装置 |
| CN116094025A (zh) * | 2022-09-07 | 2023-05-09 | 南瑞集团有限公司 | 一种含跟网型和构网型变流器功率系统暂态稳定分析方法 |
| CN116961116A (zh) * | 2023-09-20 | 2023-10-27 | 四川大学 | 基于自适应q轴电压反馈的构网逆变器暂态稳定提升方法 |
| CN117498433A (zh) * | 2023-11-03 | 2024-02-02 | 国网四川省电力公司电力科学研究院 | 一种混合并联系统暂态稳定功率极限量化方法 |
| CN117498440A (zh) * | 2023-11-17 | 2024-02-02 | 国网四川省电力公司电力科学研究院 | 并网逆变器混合并联系统暂态电流边界的分析方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| Black-Start and Service Restoration in Resilient Distribution Systems With Dynamic Microgrids;Yuhua Du et al;《IEEE JOURNAL OF EMERGING AND SELECTED TOPICS IN POWER ELECTRONICS》;20220831;第10卷(第4期);3975-3986 * |
| 基于等效同步功率的孤岛并联构网变流器系统暂态稳定性分析;邹赵悦;《电力系统自动化》;20240125;第48卷(第2期);140-150 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117728407A (zh) | 2024-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huang et al. | A virtual synchronous control for voltage-source converters utilizing dynamics of DC-link capacitor to realize self-synchronization | |
| Yang et al. | Nonlinear modeling and analysis of grid-connected voltage-source converters under voltage dips | |
| CN109149605B (zh) | 一种基于vsg的微电网暂态自适应参数控制策略 | |
| CN107465189A (zh) | 基于自适应旋转惯量的虚拟同步发电机控制方法 | |
| CN110601268B (zh) | 一种双馈风机并网端口输出阻抗建模及稳定性分析方法 | |
| CN109638881B (zh) | 电网强度自适应优化的储能逆变器虚拟同步方法及系统 | |
| CN114640141B (zh) | 海上风电二极管整流单元送出系统的构网型风机控制方法 | |
| CN106300444A (zh) | 一种三相pwm变换器功率预测分层控制方法 | |
| CN117728407B (zh) | 新能源逆变器异构并网暂态稳定问题二维相图分析方法 | |
| CN113193604A (zh) | 一种引入无功功率前馈解耦的新能源并网控制方法和装置 | |
| Wu et al. | Subsynchronous oscillation analysis of grid-connected converter based on MIMO transfer functions | |
| CN109698512B (zh) | 一种光伏集中并网系统失步振荡中心迁移特性分析方法 | |
| CN109687482B (zh) | 面向多虚拟同步机微电网有功振荡阻尼控制方法 | |
| CN114069705B (zh) | 一种判断并网电压源型换流器大扰动稳定性的方法及系统 | |
| CN112421976B (zh) | 基于混杂系统理论的三电平逆变器电源降阶建模方法 | |
| Liu et al. | Control loop stability criterion and interaction law analysis for grid-connected inverter in weak grid | |
| CN113162117A (zh) | 一种弱电网下并网逆变器控制器带宽设计方法 | |
| Zhu et al. | Research on improved virtual synchronous generator based on differential compensation link | |
| Zhao et al. | Transient Stability Analysis of Virtual Synchronous Generator Integrated Power Systems | |
| Fang et al. | Stability analysis of photovoltaic and battery energy storage systems integrated to weak grid | |
| Pan et al. | Linear hall position calculation method of PMLSM based on Type-III PLL | |
| Ren et al. | Influence of PLL on broadband oscillation under weak grid | |
| CN113972676B (zh) | 一种改善电力系统暂态稳定性的分布式bess协调控制方法 | |
| CN113659622B (zh) | 一种实现两台并联运行逆变器主从机无缝切换控制方法 | |
| Xu et al. | A virtual inertia control strategy based on multi-order lead-lag compensation |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Shunliang Inventor after: Jiang Aobo Inventor after: Ma Junpeng Inventor after: Jiao Ning Inventor after: Tu Hao Inventor after: Liu Tianqi Inventor before: Ma Junpeng Inventor before: Jiang Aobo Inventor before: Wang Shunliang Inventor before: Jiao Ning Inventor before: Tu Hao Inventor before: Liu Tianqi |