CN112564127B - 一种光伏并网电力系统频率暂态稳定性分析方法 - Google Patents
一种光伏并网电力系统频率暂态稳定性分析方法 Download PDFInfo
- Publication number
- CN112564127B CN112564127B CN202011188342.2A CN202011188342A CN112564127B CN 112564127 B CN112564127 B CN 112564127B CN 202011188342 A CN202011188342 A CN 202011188342A CN 112564127 B CN112564127 B CN 112564127B
- Authority
- CN
- China
- Prior art keywords
- frequency
- photovoltaic
- disturbance
- node
- 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.)
- Active
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 26
- 238000004458 analytical method Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 49
- 230000035699 permeability Effects 0.000 claims abstract description 31
- 238000011156 evaluation Methods 0.000 claims abstract description 22
- 238000005457 optimization Methods 0.000 claims abstract description 14
- 230000006872 improvement Effects 0.000 claims abstract description 13
- 238000010248 power generation Methods 0.000 claims description 26
- 238000004088 simulation Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 8
- 238000012795 verification Methods 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 6
- 230000035515 penetration Effects 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000001276 controlling effect Effects 0.000 claims description 4
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000005094 computer simulation Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000005034 decoration Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000010534 mechanism of action Effects 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
Images
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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- 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
- 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- 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
- H02J3/381—Dispersed generators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
-
- 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)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Marketing (AREA)
- Power Engineering (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Computation (AREA)
- Entrepreneurship & Innovation (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种光伏并网电力系统频率暂态稳定性分析方法。该模型首先对规模化光伏并网系统进行建模和分析，构建了一种规模化光伏并网系统建模方法及光伏并网系统多重故障扰动模型；在此基础上，构建了一种不同光伏渗透率下大扰动对系统频率波动的影响分析方法，研究了光伏高渗透率电网大扰动对系统频率稳定性的影响，并以最大频率偏差为指标确定出了系统可接纳光伏最大容量的区间；其次，提出了基于曲面拟合的光伏接入节点频率动态特性的评价分析与预测方法，并给出了不同并网点高渗透率光伏接入的临界值；最后，构建了基于大容量机组调频参数优化的改善措施，提高了规模化光伏接入后的系统暂态频率稳定性。
Description
技术领域
本发明涉及电力系统中频率稳定性分析领域，具体地涉及一种光伏并网电力系统频率暂态稳定性分析方法。
背景技术
太阳能光伏发电作为最具利用潜能的、普遍的、无尽的清洁能源发电方式，备受电力从业者和研究人员瞩目。然而，光伏出力具有随机性、间歇性、无机械转动惯量、不参与系统调频等特点，光伏发电的大规模接入影响着电网的安全稳定运行，给电网控制带来新挑战，使得现有控制模式难以保证电力系统安全。因此，在机理分析层面，深入研究规模化光伏接入对电力系统暂态频率稳定性的作用机理，并在光伏发电单元建模的基础上进行仿真验证，对进一步提出如何改善并网后系统暂态稳定性具有基础的理论指导作用。在稳定性评价层面，研究大规模光伏并网后系统节点频率特性，提出系统频率支撑能力评价方法对工程实践具有参考价值。在稳定性改善措施层面，提出大规模光伏并网后调频切实可行的措施，对电力系统安全稳定运行具有实际意义。
发明内容
本发明的目的，在于提供一种光伏并网电力系统频率暂态稳定性分析方法，本发明能够实现对光伏并网条件下系统受扰动后的频率动态特性的有效分析，实现光伏接入节点频率动态特性的评价与预测，并能增大系统临界扰动量，明显改善系统的频率稳定性。
为实现以上发明目的，本发明所采用的技术方案是：一种光伏并网电力系统频率暂态稳定性分析方法，本发明所述的方法具体地包括以下几个部分：
Ⅰ.一种规模化光伏并网系统建模方法及光伏并网系统多重故障扰动模型；
Ⅱ.一种不同光伏渗透率下大扰动对系统频率波动的影响分析方法；
Ⅲ.一种光伏接入节点的频率动态特性评价与预测方法；
Ⅳ.一种基于大容量机组调频参数优化的频率稳定性改善措施。
进一步地，本发明所述的一种规模化光伏并网系统模型的构建方法为：基于 PSD-BPA软件搭建规模化光伏并网系统模型，模型由多个光伏发电单元和并网系统模型构成。所述的光伏发电单元由三部分模块组成：光伏电池阵列、低压箱式变压器和电压型逆变器及其控制系统；所述光伏并网系统模型结构为：多个光伏发电单元并联接入，首先经过升压变电器升压，然后通过通道110kV汇集线路将光伏发电系统的有功输送到汇集站，规模化光伏并网系统包括了若干条光伏发电汇集支路的并联传输。
进一步地，所述的一种光伏并网系统多重故障扰动模型的构建方法为：基于n台常规发电机组系统在节点k处发生多重故障的扰动模型，将其中m台常规机组替换为光伏发电机组，得到光伏并网系统多重故障扰动模型：
式中：ΔPei和ΔP分别表示第i台发电机组电磁功率的变化量和扰动功率；Psik表示节点i、k间的整步功率系数；TJi和Ri分别表示第i台发电机的惯性时间常数和调差系数。n表示剩余常规机组的台数。
进一步地，所述的一种不同光伏渗透率下大扰动对系统频率波动的影响分析方法，包括以下步骤：
步骤1：设置电力系统内某节点为光伏接入节点；
步骤2：确定系统常规机组装机容量最大的节点，令光伏渗透率分别为该节点常规机组装机容量的0～100％；
步骤3：设置大扰动故障，借助PSD-BPA软件仿真得到不同渗透率下发生大扰动时系统的频率动态特性曲线。
进一步地，所述的频率动态特性包括：系统的节点最低最高频率动态过程、公共接入节点的节点频率动态过程、负荷最大且离并网较近节点的节点频率动态过程。
进一步地，所述的一种光伏接入节点的频率动态特性评价与预测方法，包括以下步骤：
步骤1：确定节点频率动态特性的评价指标：扰动功率、最大节点频率偏差绝对值、电气距离；
步骤2：通过在PSD-BPA中进行遍历仿真的方法获取光伏接入节点的频率动态特性的评价指标数据；
步骤3：通过改进B样条法对由BPA仿真获得的数据进行插值后拟合，生成节点频率特性曲面并进行精度校验；
步骤4：保持最大节点频率偏差绝对值不变，选取适当的电气距离，并在拟合的曲面上取点得到对应的扰动功率，将其近似作为光伏渗透容量的临界值；
进一步地，所述扰动功率ΔP定义为系统受扰动后有功功率的波动量；
所述最大节点频率偏差绝对值Δf定义为：
Δf＝|Δfmax|＝max(|fmax-fN|,|fmin-fN|) (9)
式中，fN为额定频率，fmax为系统受到功率扰动ΔP后的最大频率，fmin为系统受到功率扰动后的最大频率；
所述电气距离S定义为两点之间的转移电抗的大小。当大扰动产生有功缺额后，在扰动初瞬时不平衡功率与电气距离成反比在各机组间进行分配。
进一步地，所述的一种基于大容量机组调频参数优化的频率稳定性改善措施，包括以下步骤：
步骤1：构建节点频率支撑能力评价指标：频率裕度、临界扰动量；
步骤2：分析常规机组的一次调频参数，机组的动态特性主要影响因素为：调差系数、调速器响应时间、引导阀门时间常数和调速器死区；
步骤3：从优化调差系数和调速器死区两个方面，分析基于大容量机组调频参数优化对频率稳定性的改善情况。
进一步地，所述频率裕度指标定义为：
式中，fN、fcr和f’cr分别表示额定频率值、发电机高低频保护值和临界安全阈值；η＞0则节点暂态频率安全，反之不安全，其中η＝0为临界安全。
所述临界扰动量定义为：
该值可通过插值计算得到。
本发明的有益效果是：
本发明是一种光伏并网电力系统频率暂态稳定性分析方法。首先构建了规模化光伏并网系统模型及光伏并网系统多重故障扰动模型；在此基础上，构建了一种不同光伏渗透率下大扰动对系统频率波动的影响分析方法，研究了光伏高渗透率电网大扰动对系统频率稳定性的影响，并以最大频率偏差为指标确定出了系统可接纳光伏最大容量的区间；其次，提出了基于曲面拟合的光伏接入节点频率动态特性的评价分析与预测方法，并给出了不同并网点高渗透率光伏接入的临界值；最后，构建了基于大容量机组调频参数优化的改善措施，提高了规模化光伏接入后的系统暂态频率稳定性；通过本发明所述的一种光伏电力系统频率稳定性分析方法，可以有效分析光伏发电系统接入条件下系统受扰动后的频率动态特性，实现光伏接入节点频率动态特性的评价与预测，明显改善系统的稳定性。
附图说明
图1是光伏发电单元的拓扑结构示意图。
图2是规模化光伏并网系统模型示意图。
图3是IEEE39节点标准算例地理接线图。
图4是不同渗透率发生故障时系统最低最高频率动态特性。
图5a是不同渗透率发生故障时节点1的频率动态特性。
图5b是不同渗透率发生故障时节点39的频率动态特性。
图6是受扰情况下节点频率特性曲面。
图7a是节点暂态频率偏移安全性示意图一。
图7b是节点暂态频率偏移安全性示意图二。
图8是水轮机调速器和原动机模型的典型结构。
图9是调差系数优化设置后临界扰动量变化曲线。
图10是调速器死区优化设置后临界扰动量变化曲线。
具体实施方式
下面给出本发明的具体实施例，对本发明的技术方案作进一步具体的说明。
本具体实施方式公开了一种光伏并网电力系统频率暂态稳定性分析方法，所述方法基于一个光伏并网系统多重故障扰动模型,该扰动模型基于n台常规发电机组系统在节点k处发生多重故障的扰动模型，将其中m台常规机组替换为光伏发电机组，光伏发电机组为规模化光伏并网系统模型，模型由多个光伏发电单元和并网系统模型构成。所述的光伏发电单元由三部分模块组成：光伏电池阵列、低压箱式变压器和电压型逆变器及其控制系统；所述光伏并网系统模型结构为：多个光伏发电单元并联接入，首先经过升压变电器升压，然后通过通道110kV汇集线路将光伏发电系统的有功输送到汇集站，规模化光伏并网系统包括了若干条光伏发电汇集支路的并联传输，光伏并网系统多重故障扰动模型：
式中：ΔPei和ΔP分别表示第i台发电机组电磁功率的变化量和扰动功率；Psik表示节点i、k间的整步功率系数；TJi和Ri分别表示第i台发电机的惯性时间常数和调差系数。n表示剩余常规机组的台数，分析方法具体包括：
不同光伏渗透率下大扰动对系统频率波动的影响分析步骤；
步骤1：设置电力系统内某节点为光伏接入节点；
步骤2：确定系统常规机组装机容量最大的节点，令光伏渗透率分别为该节点常规机组装机容量的0～100％；
步骤3：设置大扰动故障，借助PSD-BPA软件仿真得到不同渗透率下发生大扰动时系统的频率动态特性曲线。
光伏接入节点的频率动态特性评价与预测步骤：
步骤1：通过在PSD-BPA中进行遍历仿真的方法获取光伏接入节点的频率动态特性的评价指标数据，评价指标包括：扰动功率、最大节点频率偏差绝对值、电气距离；
步骤2：通过改进B样条法对由BPA仿真获得的数据进行插值后拟合，生成节点频率特性曲面并进行精度校验；
步骤3：保持最大节点频率偏差绝对值不变，根据光伏并网点的位置确定电气距离，并在拟合的曲面上取点得到对应的扰动功率，将其近似作为光伏渗透容量的临界值；
基于大容量机组调频参数优化的频率稳定性改善步骤：
步骤1：构建节点频率支撑能力评价指标：频率裕度、临界扰动量；
步骤2：分析常规机组的一次调频参数，机组的动态特性主要影响因素为：调差系数、调速器响应时间、引导阀门时间常数和调速器死区；
步骤3：从优化调差系数和调速器死区两个方面，分析基于大容量机组调频参数优化对频率稳定性的改善情况。
所述三个步骤是基于光伏扰动模型的独立步骤。通过对不同光伏渗透率下大扰动对系统频率波动的影响进行分析，可以确定系统频率动态特性与光伏渗透率的关系；通过光伏接入节点的频率动态特性评价与预测，可以确定系统内某点光伏装机容量的临界值；通过基于大容量机组调频参数优化的频率稳定性改善措施，可以增大并网节点的临界光伏容量，提高光伏并网后系统的频率稳定性。
下面将详细说明本发明的实施方式。
1，在本实施方式中，构建了一种规模化光伏并网系统建模方法及光伏并网系统多重故障扰动模型。下面给出本发明第一部分的具体实施例：
(1)规模化光伏并网系统建模。该模型由多个光伏发电单元和并网系统模型构成，所述的光伏发电单元由三部分模块组成：光伏电池阵列、低压箱式变压器和电压型逆变器及其控制系统，如图1所示；所述光伏并网系统模型结构为：多个光伏发电单元并联接入，首先经过升压变电器升压，然后通过通道110kV汇集线路将光伏发电系统的有功输送到汇集站，规模化光伏并网系统包括了若干条光伏发电汇集支路的并联传输，如图2 所示。
(2)光伏并网系统多重故障扰动建模。所述模型的构建方法为：基于n台常规发电机组系统在节点k处发生多重故障的扰动模型，将其中m台常规机组替换为光伏发电机组，得到光伏并网系统多重故障扰动模型：
式中：ΔPei和ΔP分别表示第i台发电机组电磁功率的变化量和扰动功率；Psik表示节点i、k间的整步功率系数；TJi和Ri分别表示第i台发电机的惯性时间常数和调差系数。n表示剩余常规机组的台数。
所述模型用以分析光伏高渗透率电网在大扰动下频率稳定性的影响机理。根据模型，受大扰动情况下，系统随着光伏渗透率增大，剩余常规机组将承担更多功率缺额，其电磁功率的瞬时变化和稳态偏差也将更为明显，即系统发生暂态频率稳定问题愈加容易。
2，在本实施方式中，构建了一种不同光伏渗透率下大扰动对系统频率波动的影响分析方法。下面给出本发明第二部分的具体实施例：
(1)基于IEEE39节点标准算例，设置bus2节点为光伏电站并网节点，IEEE39节点标准算例地理接线图如图3所示；
(2)确定系统常规机组装机容量最大的节点为bus39，令光伏渗透率分别为该节点常规机组装机容量的0～100％；
(3)设置如下两个大扰动故障：
①支路bus6-bus11在第5个周波时发生三相短路故障，在第10个周波时清除故障；
②增加光伏有功出力时，等比切除常规机组有功出力。
借助PSD-BPA软件仿真得到不同渗透率下发生大扰动时该系统的频率动态特性曲线，包括系统的节点最低最高频率动态过程和公共接入点bus1及负荷最大且离并网点较近的节点bus39的节点频率动态过程，如表1和图4、5所示。
表1不同光伏渗透率发生故障时系统的频率动态响应
根据仿真结果，系统的最低最高频率及母线1和母线39节点动态频率随着光伏渗透率的不同而发生变化。系统的最低最高动态频率变化趋势相同，且在受扰动情况下，光伏高渗透率系统既有发生频率上限越限的可能，也有发生频率下限越限的可能，即会出现高频和低频问题。其原因主要是随着光伏渗透率的增大，系统惯性减小，功角越容易失稳。
3，在本实施方式中，构建了一种光伏接入节点的频率动态特性评价与预测方法。下面给出本发明第三部分的具体实施例：
(1)确定节点频率动态特性的评价指标：扰动功率、最大节点频率偏差绝对值、电气距离。
①扰动功率ΔP：系统受扰动后有功功率的波动量；
②最大节点频率偏差绝对值Δf：
Δf＝|Δfmax|＝max(|fmax-fN|,|fmin-fN|) (14)
式中，fN为额定频率，fmax为系统受到功率扰动ΔP后的最大频率，fmin为系统受到功率扰动后的最大频率；
③电气距离S：两点之间的转移电抗大小。当大扰动产生有功缺额后，在扰动初瞬时不平衡功率与电气距离成反比在各机组间进行分配。
(2)基于IEEE39节点标准算例，选取光伏接入节点作为观测点，依次设定切机点和扰动量，借助PSD-BPA软件，通过遍历仿真的方法多次改变潮流文件和稳态文件中的扰动点位置和扰动量大小，得到(ΔP,S,Δf)的一系列实验数据。
(3)对得到的实验数据利用三次样条插值方法，通过MATLAB的程序编写实现“频率偏差-电气距离-光伏渗透率”的曲面拟合，并利用改进的B样条法对由BPA仿真获得的有限个数据点所拟合出的不光滑曲面进行修正，得到受扰情况下的节点特性曲面如图 6所示。最后按照下式对拟合曲面进行精确度校验：
式中，Δf和Δf’分别为频率指标的拟合值和仿真值。
按上式检验拟合精度，具体方法为：首先确定一个电气距离，用它在不同的扰动量下对应的频率偏差绝对值进行精度校验；然后以同样的方法，再确定一个扰动量，用它在不同的电气距离下对应的频率偏差绝对值进行精度校验。预测值可以根据MATLAB 中的句柄函数从绘制的曲面中读取，具体精度校验结果如表2所示：
表2分别固定扰动量和电气距离时的频率指标值校验
经过精确度校验可知，本发明所提的一种光伏接入节点的频率动态特性评价与预测方法满足足够的精确度，误差在2％以内。
(4)保持最大节点频率偏差绝对值不变，选取适当的电气距离，通过在拟合的曲面上取点得到对应的扰动功率，即可将其近似作为光伏渗透容量的临界值。如，当S取0.0012时，本具体实施例中频率指标所对应的临界光伏渗透容量为272.4MW。
4，在本实施方式中，构建了一种基于大容量机组调频参数优化的频率稳定性改善措施。下面给出本发明第四部分的具体实施例：
(1)构建节点频率支撑能力评价指标：频率裕度、临界扰动量。节点暂态频率偏移安全性如图7所示。
①频率裕度指标定义为：
式中，fN、fcr和f’cr分别表示额定频率值、发电机高低频保护值和临界安全阈值；η＞0则节点暂态频率安全，反之不安全，其中η＝0为临界安全。
②临界扰动量定义为：
该值可通过插值计算得到。
(2)以常规水轮机组调速系统为例，分析一次调频参数优化对改善系统频率稳定性的影响。水轮机调速器和原动机模型如图8所示。
调速器的闭环传递函数为：
根据发电机动态稳定特征方程的判据推导可得：
式中，软反馈系数与水锤效应时间常数有关，因此机组的动态特性主要影响因素为调差系数、调速器响应时间、引导阀门时间常数和调速器死区。
(3)从优化调差系数和调速器死区两个方面，分析基于大容量机组调频参数优化对频率稳定性的改善情况。基于IEEE39节点标准算例，选取bus32进行仿真，设置如下步骤：
①控制调速器死区不变，分别设置调差系数，计算该系统的临界扰动量；
②控制调差系数不变，分别设置调速器死区，计算该系统的临界扰动量。
优化设置后的结果分别如图9和图10所示。根据结果，满足频率暂态稳定裕度指标的临界扰动量随着调差系数的增大，呈现出先增大后减小的趋势；随着调速器死区的增大，呈现出单调减小的趋势。本具体实施例中，通过优化大容量机组的调差系数可增大基于节点频率特性曲面的临界扰动量，使频率稳定性得到明显改善。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在本发明所述技术方案的基础上，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种光伏并网电力系统频率暂态稳定性分析方法，其特征在于，基于一个光伏并网系统多重故障扰动模型,该扰动模型基于n台常规发电机组系统在节点k处发生多重故障的扰动模型，将其中m台常规机组替换为光伏发电机组，光伏发电机组为规模化光伏并网系统模型，模型由多个光伏发电单元和并网系统模型构成；所述的光伏发电单元由三部分模块组成：光伏电池阵列、低压箱式变压器和电压型逆变器及其控制系统；所述光伏并网系统模型结构为：多个光伏发电单元并联接入，首先经过升压变电器升压，然后通过通道110kV汇集线路将光伏发电系统的有功输送到汇集站，规模化光伏并网系统包括了若干条光伏发电汇集支路的并联传输，光伏并网系统多重故障扰动模型：
式中：ΔPei和ΔP分别表示第i台发电机组电磁功率的变化量和扰动功率；Psik表示节点i、k间的整步功率系数；TJi和Ri分别表示第i台发电机的惯性时间常数和调差系数；n'表示剩余常规机组的台数，分析方法具体包括：
不同光伏渗透率下大扰动对系统频率波动的影响分析步骤；
步骤1：设置电力系统内某节点为光伏接入节点；
步骤2：确定系统常规机组装机容量最大的节点，令光伏渗透率分别为该节点常规机组装机容量的0～100％；
步骤3：设置大扰动故障，借助PSD-BPA软件仿真得到不同渗透率下发生大扰动时系统的频率动态特性曲线；
光伏接入节点的频率动态特性评价与预测步骤：
步骤1：通过在PSD-BPA中进行遍历仿真的方法获取光伏接入节点的频率动态特性的评价指标数据，评价指标包括：扰动功率、最大节点频率偏差绝对值、电气距离；
步骤2：通过改进B样条法对由BPA仿真获得的数据进行插值后拟合，生成节点频率特性曲面并进行精度校验；
步骤3：保持最大节点频率偏差绝对值不变，根据光伏并网节点的位置确定电气距离，并在拟合的曲面上取点得到对应的扰动功率，将其作为光伏渗透容量的临界值；
基于大容量机组调频参数优化的频率稳定性改善步骤：
步骤1：构建节点频率支撑能力评价指标：频率裕度、临界扰动量；
步骤2：分析常规机组的一次调频参数，机组的动态特性主要影响因素为：调差系数、调速器响应时间、引导阀门时间常数和调速器死区；
步骤3：从优化调差系数和调速器死区两个方面，分析基于大容量机组调频参数优化对频率稳定性的改善情况；
所述三个步骤是基于光伏扰动模型的独立步骤；通过对不同光伏渗透率下大扰动对系统频率波动的影响进行分析，可以确定系统频率动态特性与光伏渗透率的关系；通过光伏接入节点的频率动态特性评价与预测，可以确定系统内某点光伏装机容量的临界值；通过基于大容量机组调频参数优化的频率稳定性改善措施，可以增大并网节点的临界光伏容量，提高光伏并网后系统的频率稳定性。
2.如权利要求1所述的一种光伏并网电力系统频率暂态稳定性分析方法，其特征在于：所述影响分析步骤的步骤3中，频率动态特性包括：系统的节点最低最高频率动态过程、公共接入节点的节点频率动态过程、负荷最大且离并网较近节点的节点频率动态过程。
3.如权利要求2所述的一种光伏并网电力系统频率暂态稳定性分析方法，其特征在于：光伏接入节点的频率动态特性评价与预测步骤包括：
步骤1：通过在PSD-BPA中进行遍历仿真的方法获取光伏接入节点的频率动态特性的评价指标数据，评价指标包括：扰动功率、最大节点频率偏差绝对值、电气距离；
步骤2：通过改进B样条法对由BPA仿真获得的数据进行插值后拟合，生成节点频率特性曲面并进行精度校验；
所述步骤的具体方式是对得到的(ΔP,S,Δf)数据利用三次样条插值方法，通过MATLAB程序编写实现“频率偏差-电气距离-新能源渗透率”的曲面拟合；并利用改进的B样条法对由BPA仿真获得的有限个数据点所拟合出的不光滑曲面进行修正；
所述改进B样条法的具体思路是将平面分布构造成可柔性光滑变化的B样条曲面，当实际数据位于平面区域外时，柔性自适应地调整计算面的形状，使该数据点被包含到对应的曲面区域内，调整后的曲面仍然清晰，但其形状不再是平面，而是随数据样本形状而改变，对应的曲面区域包含了所有的数据样本，从而提高了曲面构造的拟合精度；
精确度校验公式是：
式中，Δf和Δf'分别为频率指标的拟合值和仿真值；
步骤3：保持最大节点频率偏差绝对值不变，根据光伏并网节点的位置确定合适的电气距离S，并在拟合的曲面上取点得到对应的扰动功率ΔP，将其作为光伏渗透容量的临界值。
4.如权利要求3所述的一种光伏并网电力系统频率暂态稳定性分析方法，其特征在于：光伏接入节点的频率动态特性评价与预测步骤中，
扰动功率ΔP定义为系统受扰动后有功功率的波动量；
最大节点频率偏差绝对值Δf定义为：
Δf＝|Δfmax|＝max(|fmax-fN|,|fmin-fN|) (3)
式中，fN为额定频率，fmax为系统受到功率扰动ΔP后的最大频率，fmin为系统受到功率扰动后的最大频率；
电气距离S定义为两点之间的转移电抗的大小；当大扰动产生有功缺额后，在扰动初瞬时不平衡功率与电气距离成反比在各机组间进行分配。
5.如权利要求4所述的一种光伏并网电力系统频率暂态稳定性分析方法，其特征在于：
所述频率裕度指标和临界扰动量两项指标用于评估系统频率稳定性，因为在暂态频率失稳事故中系统频率必然先达到触发发电机组的跳闸条件；
频率裕度指标定义为：
式中，fN、fcr和f′cr分别表示额定频率值、发电机高低频保护值和临界安全阈值；η＞0则节点暂态频率安全，反之不安全，其中η＝0为临界安全；
所述临界扰动量定义为：
该值可通过插值计算得到；
所述常规机组的一次调频参数分析方法为，基于水轮机调速器和原动机模型，得出调速器的闭环传递函数为：
根据发电机动态稳定特征方程的判据推导可得：
Td 2(R+Dd)+Td(TG+DdTP)+TGTP＞0 (7)
所述常规机组的动态特性主要影响因素为调差系数、调速器响应时间、引导阀门时间常数和调速器死区；所述基于调差系数和调速器死区两个方面进行调频参数优化的光伏并网电力系统频率稳定性改善方法的具体思路为：
①控制调速器死区不变，分别设置调差系数，计算该系统的临界扰动量；
②控制调差系数不变，分别设置调速器死区，计算该系统的临界扰动量；
③通过仿真得到调差系数、调速器死区与临界扰动量的关系曲线，根据关系曲线设定适当的调差系数和调速器死区数值，提高并网点的临界光伏容量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011188342.2A CN112564127B (zh) | 2020-10-30 | 2020-10-30 | 一种光伏并网电力系统频率暂态稳定性分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011188342.2A CN112564127B (zh) | 2020-10-30 | 2020-10-30 | 一种光伏并网电力系统频率暂态稳定性分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112564127A CN112564127A (zh) | 2021-03-26 |
| CN112564127B true CN112564127B (zh) | 2022-07-05 |
Family
ID=75042712
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011188342.2A Active CN112564127B (zh) | 2020-10-30 | 2020-10-30 | 一种光伏并网电力系统频率暂态稳定性分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112564127B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113541135B (zh) * | 2021-07-23 | 2022-09-02 | 哈尔滨工业大学 | 以暂态频率跌落基准轨迹为统一尺度的广义惯性测度方法 |
| CN114123233B (zh) * | 2021-08-20 | 2023-10-03 | 华北电力大学（保定） | 含不同光伏占比系统的暂态频率稳定性量化评估方法 |
| CN114723273B (zh) * | 2022-03-31 | 2024-06-18 | 华中科技大学 | 并网水电站系统不确定非线性模型构建、分析方法和系统 |
| CN114696323B (zh) * | 2022-04-12 | 2023-05-26 | 国网湖南省电力有限公司 | 高比例新能源电力系统频率安全判定方法 |
| CN114725939B (zh) * | 2022-05-31 | 2022-08-16 | 中国电力科学研究院有限公司 | 一种确定电网系统可承受的最大扰动功率的方法及系统 |
| CN115693708A (zh) * | 2022-11-01 | 2023-02-03 | 南通理工学院 | 一种估算电力系统调频暂态频率最低点的方法 |
| CN117439518B (zh) * | 2023-12-19 | 2024-02-20 | 成都纺织高等专科学校 | 一种光伏并网系统及方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104361260A (zh) * | 2014-12-05 | 2015-02-18 | 武汉大学 | 基于曲面拟合的电力系统节点动态频率响应特性分析方法 |
| KR20150061133A (ko) * | 2013-11-26 | 2015-06-04 | 한국전기연구원 | 신재생 발전 단지 출력 제어 방법 및 장치 |
| CN109038664A (zh) * | 2018-08-02 | 2018-12-18 | 山东大学 | 计及风机涉网保护高频分轮次整定方法 |
| CN109936152A (zh) * | 2018-11-09 | 2019-06-25 | 西南交通大学 | 高渗透率风电并网后的电网调频控制方法、风力发电机 |
| CN111224434A (zh) * | 2020-03-12 | 2020-06-02 | 安徽工程大学 | 光火储混合发电系统的负荷频率协调优化控制方法 |
-
2020
- 2020-10-30 CN CN202011188342.2A patent/CN112564127B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20150061133A (ko) * | 2013-11-26 | 2015-06-04 | 한국전기연구원 | 신재생 발전 단지 출력 제어 방법 및 장치 |
| CN104361260A (zh) * | 2014-12-05 | 2015-02-18 | 武汉大学 | 基于曲面拟合的电力系统节点动态频率响应特性分析方法 |
| CN109038664A (zh) * | 2018-08-02 | 2018-12-18 | 山东大学 | 计及风机涉网保护高频分轮次整定方法 |
| CN109936152A (zh) * | 2018-11-09 | 2019-06-25 | 西南交通大学 | 高渗透率风电并网后的电网调频控制方法、风力发电机 |
| CN111224434A (zh) * | 2020-03-12 | 2020-06-02 | 安徽工程大学 | 光火储混合发电系统的负荷频率协调优化控制方法 |
Non-Patent Citations (2)
| Title |
|---|
| 光伏渗透率对电力系统频率的影响及控制策略;蒋小亮 等;《武汉大学学报(工学版)》;20151031;第686-696页 * |
| 调频关键参数对电网频率特性的影响及其灵敏度分析;王凡 等;《电力系统保护与控制》;20201016;第1-8页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112564127A (zh) | 2021-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112564127B (zh) | 一种光伏并网电力系统频率暂态稳定性分析方法 | |
| CN101931241B (zh) | 风电场并网协调控制方法 | |
| CN111799813B (zh) | 计及风电机组无功调节的海上风电场无功优化配置方法 | |
| CN102299527B (zh) | 一种风电场无功功率控制方法和系统 | |
| CN108695857B (zh) | 风电场自动电压控制方法、装置及系统 | |
| CN109361242B (zh) | 一种光伏发电自动电压控制方法 | |
| CN102611118A (zh) | 一种引入预测信息的风电场综合无功电压控制方法 | |
| CN103346577A (zh) | 降低风电场功率损耗的风电场avc无功控制系统及方法 | |
| CN108711868A (zh) | 一种计及孤岛运行电压安全的配电网无功优化规划方法 | |
| CN107065556A (zh) | 一种堆芯机组变功率运行策略优化方案的自动搜索方法 | |
| KR102133897B1 (ko) | 실시간 계측 기반 전압 및 무효전력 제어 방법 및 그 프로그램 | |
| CN103490426A (zh) | 大规模风电接入的弱送端电力系统电压稳定综合控制方法 | |
| CN104104102A (zh) | 电压源换流器型多端直流输电系统稳态工作点优化方法 | |
| CN112491043A (zh) | 一种新能源富集电网电源规划方法及系统 | |
| CN113937825A (zh) | 一种基于E-C-Kmeans聚类和SOP优化的DG双层优化配置方法 | |
| CN115017854A (zh) | 一种基于多维评估指标体系的配电网dg最大准入容量计算方法 | |
| CN109546677A (zh) | 一种规模化海上风电柔性送出系统安全控制策略求解方法 | |
| CN115842376A (zh) | 电力系统等效惯量趋势与安全状态评估方法、设备和介质 | |
| CN109980660B (zh) | 基于发电机改进频率响应特性的鲁棒区间风电调度方法 | |
| CN111030187A (zh) | 一种多能源分频互补新能源的补偿容量优化计算方法 | |
| CN105262112A (zh) | 风电场集群式静止型无功补偿装置控制方法 | |
| US11881713B1 (en) | Digital twin advanced distribution management systems (ADMS) and methods | |
| Hong et al. | A study of evaluation system based on large scale photovoltaic power generation | |
| CN112736913A (zh) | 含分布式电源的配电网功率优化模式影响因素分析方法 | |
| Elahidoost et al. | Reducing HVDC network oscillations considering wind intermittency through optimized grid expansion decision |
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 |