CN107516903A - 一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法 - Google Patents
一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法 Download PDFInfo
- Publication number
- CN107516903A CN107516903A CN201710774974.9A CN201710774974A CN107516903A CN 107516903 A CN107516903 A CN 107516903A CN 201710774974 A CN201710774974 A CN 201710774974A CN 107516903 A CN107516903 A CN 107516903A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- frequency
- voltage
- broad sense
- 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
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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
-
- 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/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
- Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
-
- 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
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
-
- 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
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
- Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法，属于电力系统及其自动化技术领域。本发明针对故障发展过程中系统不同时间尺度面临的安全稳定问题，计及切负荷措施代价和电力安全事故责任代价，建立了暂态广义经济损失、准稳态广义经济损失及稳态广义经济损失指标，实现了统筹毫秒级、秒级及分钟级的精准负荷控制。本发明有利于电力系统调度运行人员把握复杂故障后系统不同时间尺度安全运行规律，提升现有三道防线负荷控制手段的适应性，从而提高大电网安全稳定控制的有效性和精益化水平。
Description
技术领域
本发明属于电力系统及其自动化技术领域，更准确地说本发明涉及一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法。
背景技术
在特高压电网建设的过渡期，电网“强直弱交”矛盾突出，特高压直流输送容量在受端系统规模中的比例不断增大，直流闭锁故障造成大功率冲击及潮流窜动，全网性电压问题及频率问题突显。针对受端电网由于大扰动产生的功率缺额，为防止发生稳定破坏，将按预定的切负荷策略，采用安全稳定控制装置实施切负荷控制措施；针对电力系统发生频率、电压异常等事故，采用低频低压减负荷措施，防止系统频率、电压崩溃；为了解决断面越限、联络线超用、旋转备用不足等，采用预设的过载值根据事先制定的限电序位表进行拉路控制等。由于国务院第599号令《电力安全事故应急处置和调查处理条例》(以下简称“条例”)明确了稳控系统切负荷等同于故障损失负荷，区域切负荷比例过高或负荷量分配的不合理均会导致更为严重的事故等级评级和事后追责，因此，在执行紧急负荷控制时一般又会考虑经济代价。
现有负荷控制一般是以最小切负荷量或最小经济损失为目标，单一解决故障后某一时间尺度稳定问题，如解决暂态电压稳定问题、频率稳定或稳态低电压或过载等问题。其本质上大多按照事件、轨迹驱动，以切除主变和高压负荷线为主，没有从不同时间尺度控制角度对负荷控制进行协调优化，对不同切负荷对象也不加选择，当有大规模切负荷控制要求时，其经济性和社会接受度均较低。因此，有必要研究多时间尺度精准负荷控制技术，根据系统暂态及稳态不同时间尺度综合安全稳定需求及经济性代价，对大范围大规模的负荷控制进行统筹管理，提出负荷控制经济性与电网安全性相协调的优化控制方案，从而提升现有三道防线负荷控制手段的适应性，提高过渡期特高压交直流混联电网安全稳定和电网控制管理精益化水平。
发明内容
本发明的目的是：针对现有技术的不足，提高负荷控制方法对于暂态及稳态综合安全稳定、经济性的适应性，避免解决故障后单一稳定问题而导致的后续不同时间尺度稳定问题、切负荷效率低及经济代价高等问题，给出一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法。
具体地说，本发明采用以下的技术方案来实现的，包括下列步骤：
1)实时测量当地电网电气量，在故障发生到实测结束时刻的给定时间窗口内，获取系统频率及母线电压的动态响应曲线，确定切负荷措施动作后导致的切负荷后常规经济损失与电力安全事故责任代价；
2)将故障后不同时间尺度的安全稳定问题折算成广义经济损失，确定广义经济损失代价，所述不同时间尺度包括毫秒级、秒级、分钟级；
3)综合考虑切负荷后常规经济损失、电力安全事故责任代价及广义经济损失代价，得到紧急负荷控制数学模型，并采用原对偶内点法进行不同时间尺度最优切负荷求解，求解次数由不同时间尺度安全约束条件是否满足决定，且下一次模型求解的初值为上一次模型最优解；
4)将不同时间尺度最优切负荷解提交负荷控制中心站，当负荷控制中心站判定不同时间尺度安全稳定约束低至动作门槛时，实时向子站下发决策结果，子站结合就地频率，对负控终端实时采集的用户所有380V分支回路功率进行负荷精准控制。
上述技术方案的进一步特征在于，所述步骤1)中，切负荷措施动作后导致的切负荷后常规经济损失与电力安全事故责任代价，按下式表示：
FJ＝FJ1+FJ2
式中：FJ为切负荷后常规经济损失与电力安全事故责任代价之和，FJ1为稳控系统切负荷措施动作导致的切负荷后常规经济损失，FJ2为根据不同等级电网导致相应电力安全事故的责任代价。
上述技术方案的进一步特征在于，所述步骤2)中，具体包括以下步骤：
2-1)根据暂态功角稳定性、暂态频率安全性和暂态电压安全性将毫秒级暂态功角、电压、频率稳定情况折算成暂态广义经济损失，如下式所示：
Fms＝αFms,δ+βFms,V+γFms,f
式中：Fms为将毫秒级暂态功角、电压、频率稳定情况折算成的暂态广义经济损失，Fms,δ、Fms,V、Fms,f分别为通过时域仿真得到故障仿真轨迹并从中提取出暂态功角、暂态电压、暂态频率的量化信息的暂态安全稳定性指标，α、β、γ分别为将故障后系统暂态功角、暂态电压、暂态频率指标换算成相应广义经济损失的折算系数；Adec为故障后一定仿真时间内减速面积，Ainc为加速面积；V、f分别为实际电压、频率仿真轨迹，Vcr和tcr为预设的电压偏移二元表(Vcr,tcr)的元素，fcr和tcr'为预设的频率偏移二元表(fcr,tcr')的元素，其中Vcr、fcr分别为电压偏移、频率偏移门槛值，t、t'为积分起始时间，分别取为暂态过程中实际电压曲线与Vcr首次交叉的时间、频率仿真曲线与fcr首次交叉的时间，tcr、tcr'分别为电压、频率偏移最大容许时间，VN为系统的额定电压，fN为系统的额定频率；
2-2)根据事故后系统频率恢复及电压恢复情况将秒级电压、频率稳定情况折算成准稳态广义经济损失，如下式所示：
Fs＝λFs,V+μFs,f
式中：Fs为将秒级电压、频率稳定情况折算成的准稳态广义经济损失，Fs,V、Fs,f分别为满足故障后系统暂态稳定并从准稳态仿真轨迹提取系统电压、频率的量化信息的准稳态安全稳定性指标，λ、μ为将故障后系统恢复电压、恢复频率指标换算成相应广义经济损失的折算系数；k为不满足准稳态恢复电压的薄弱节点总数，Vi为第i个节点的恢复电压，为恢复电压下限；fs为系统恢复频率，为恢复频率下限；
2-3)根据事故后系统安全备用、联络断面稳定情况将分钟级安全备用不足、大区联络线超用稳定情况折算成稳态广义经济损失，如下式所示：
式中：Fmin为将分钟级安全备用不足、大区联络线超用稳定情况折算成的稳态广义经济损失，Fmin,s、Fmin,p分别为满足故障后系统暂态稳定及准稳态稳定并从系统稳态恢复中提取系统安全备用、联络断面稳定情况的量化信息的稳态安全指标，χ分别为将系统稳态安全备用、联络断面指标换算成相应广义经济损失的折算系数；S为事故后系统恢复安全备用，SN为满足系统安全要求的备用水平；m为联络断面过载的线路总数，Pj为第j个过载线路的稳态恢复功率，PjN为第j个过载线路的额定功率。
上述技术方案的进一步特征在于，所述步骤3)中，紧急负荷控制数学模型表示如下：
minF＝FJ+Fms+Fs+Fmin
式中：F为将切负荷后常规经济损失与电力安全事故责任代价、暂态广义经济损失、准稳态广义经济损失、稳态广义经济损失进行加权的和。
紧急负荷控制数学模型的约束条件包括：暂态安全稳定性临界值、准稳态电压及频率下限、系统安全备用及联络断面约束以及各切负荷点的减负荷上限约束。
本发明的有益效果如下：本发明针对故障发展过程中系统不同时间尺度面临的安全稳定问题，计及切负荷措施代价和电力安全事故责任代价，建立了暂态广义经济损失、准稳态广义经济损失及稳态广义经济损失指标，实现了统筹毫秒级、秒级及分钟级的精准负荷控制。本发明有利于电力系统调度运行人员把握复杂故障后系统不同时间尺度安全运行规律，提升现有三道防线负荷控制手段的适应性，从而提高大电网安全稳定控制的有效性和精益化水平。
附图说明
图1为本发明方法的流程图。
具体实施方式
下面参照附图并结合实例对本发明作进一步详细描述。
实施例1：
本实施例为一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法，具体步骤如图1所示。
图1中步骤1描述的是，从在线决策系统平台实时测量当地电网电气量，在故障发生到实测结束时刻的给定时间窗口内，获取系统频率、母线电压等动态响应曲线，确定切负荷措施动作后导致的常规经济损失和电力安全事故责任代价。其中，常规经济损失和电力安全事故责任代价的表达式为：
FJ＝FJ1+FJ2
式中：FJ为切负荷后常规经济损失与电力安全事故责任代价之和，FJ1为稳控系统切负荷措施动作导致的常规经济损失；FJ2为根据不同等级电网导致相应电力安全事故的责任代价函数。
FJ越大，说明对应的切负荷措施所造成的损失越大，当FJ过大时，应考虑对运行方式和稳控措施的配置做出适当调整。
图1中步骤2描述的是，将故障后不同时间尺度的安全稳定问题折算成广义经济损失，确定广义经济损失代价，所述不同时间尺度包括毫秒级、秒级、分钟级。具体包括以下步骤：
图1中步骤2-1描述的是，将毫秒级暂态功角、电压、频率稳定情况折算成暂态广义经济损失，此过程需要考察三个方面，即暂态功角稳定性、暂态频率安全性和暂态电压安全性，毫秒级暂态稳定性越差，其折算成的暂态广义经济损失越大，表达式为：
Fms＝αFms,δ+βFms,V+γFms,f
式中：Fms为将毫秒级暂态功角、电压、频率稳定情况折算成的暂态广义经济损失，Fms,δ、Fms,V、Fms,f分别为通过时域仿真得到故障仿真轨迹并从中提取出暂态功角、暂态电压、暂态频率的量化信息的暂态安全稳定性指标。α、β、γ为将故障后系统暂态功角、暂态电压、暂态频率指标换算成相应广义经济损失的折算系数，其取值取决于用户对三种暂态稳定问题的容许程度；Adec为故障后一定仿真时间内减速面积，Ainc为加速面积；V、f分别为实际电压、频率仿真轨迹，Vcr和tcr为预设的电压偏移二元表(Vcr,tcr)的元素，fcr和tcr'为预设的频率偏移二元表(fcr,tcr')的元素，其中Vcr、fcr分别为电压偏移、频率偏移门槛值，t、t'为积分起始时间，分别取为暂态过程中实际电压曲线与Vcr首次交叉的时间、频率仿真曲线与fcr首次交叉的时间，tcr、tcr'分别为电压、频率偏移最大容许时间，VN为系统的额定电压，fN为系统的额定频率。
由于系统不同母线节点电压具有差异性，而不同中枢节点频率差异无几，故上述V可通过戴维南等值阻抗法判断系统电压最薄弱的母线节点，再读取该节点的实际电压仿真轨迹得到，而f可直接读取为高电压等级节点的实际频率仿真轨迹得到。
毫秒级负荷控制措施实施对象一般由客户自主选择一部分非核心可短时中断的用电负荷，如启停方便的生产线和空调用电部分照明用电等可中断负荷，当施加毫秒级负荷控制措施后，Fms,δ、Fms,V、Fms,f均达到暂态安全稳定性临界值(一般取为0)后，进入秒级负荷控制实施范畴。
图1中步骤2-2描述的是，将秒级电压、频率稳定情况折算成准稳态广义经济损失，此过程需要考察两个方面，即事故后系统频率恢复及电压恢复。准稳态稳定问题恢复越差，其折算成的准稳态广义经济损失越大，表达式为：
Fs＝λFs,V+μFs,f
式中：Fs为将秒级电压、频率稳定情况折算成的准稳态广义经济损失，Fs,V、Fs,f分别为满足故障后系统暂态稳定并从准稳态仿真轨迹提取系统电压、频率的量化信息的准稳态安全稳定性指标。λ、μ为将故障后系统恢复电压、恢复频率指标换算成相应广义经济损失的折算系数，其取值取决于用户对两种准稳态稳定问题的容许程度；k为不满足准稳态恢复电压的薄弱节点总数，Vi为第i个节点的恢复电压。一般而言，准稳态恢复电压应达到{0.9p.u.,1.1p.u.}，即应恢复至额定电压的90％至110％之间，为恢复电压下限；fs为系统恢复频率，为恢复频率下限。
通过秒级负荷控制措施的精准优选，避免原有低频低压减载装置集中切除负荷的弊端，当施加秒级负荷控制措施后，系统恢复电压、频率均大于电压、频率下限，进入分钟级负荷控制实施范畴。
图1中步骤2-3描述的是，将分钟级安全备用不足、大区联络线超用稳定情况折算成稳态广义经济损失，此过程需要考察两个方面，即事故后系统安全备用、联络断面稳定问题恢复越差，其折算成的稳态广义经济损失越大，表达式为：
式中：Fmin为将分钟级安全备用不足、大区联络线超用稳定情况折算成的稳态广义经济损失，Fmin,s、Fmin,p分别为满足故障后系统暂态稳定及准稳态稳定并从系统稳态恢复中提取系统安全备用、联络断面稳定情况的量化信息的稳态安全指标。χ为将系统稳态安全备用、联络断面指标换算成相应广义经济损失的折算系数，其取值取决于用户对两种稳态稳定问题的容许程度；S为事故后系统恢复安全备用，SN为满足系统安全要求的备用水平；m为联络断面过载的线路总数，Pj为第j个过载线路的稳态恢复功率，PjN为第j个过载线路的额定功率。
当施加分钟级负荷控制措施后，系统安全备用大于安全要求备用水平、联络断面功率小于线路额定功率，即为满足多时间尺度安全稳定性要求的精准负荷控制。
图1中步骤3描述的是，综合考虑切负荷后常规经济损失、电力安全事故责任代价及广义经济损失代价，得到紧急负荷控制数学模型，并采用原对偶内点法进行不同时间尺度最优切负荷求解。具体包括以下步骤：
图1中步骤3-1描述的是，由于精准负荷控制的目标是在控制代价和控制后系统多时间安全稳定性之间求得平衡，在不过多增加控制代价的前提下，尽量提高事故后不同时间尺度内系统的稳定性，故将紧急负荷控制数学模型表示如下：
minF＝FJ+Fms+Fs+Fmin
式中：F为将切负荷后常规经济损失与电力安全事故责任代价、暂态广义经济损失、准稳态广义经济损失、稳态广义经济损失进行加权的和。
紧急负荷控制数学模型的约束条件包括：暂态安全稳定性临界值、准稳态电压及频率下限、系统安全备用及联络断面约束以及各切负荷点的减负荷上限约束。
图1中步骤3-2描述的是，将不同时间尺度的负荷控制分阶段优化，并分别代入计算切负荷措施代价和电力安全事故责任代价，采用常规的原对偶内点法求解最优减负荷率。由于紧急负荷控制数学模型涉及了毫秒级、秒级及分钟级三种不同时间尺度，需要分别进行至少一次、最多三次的模型求解，具体求解次数由不同时间尺度安全约束条件是否满足决定，且下一次模型求解的初值为上一次模型最优解。
图1中步骤4描述的是，将不同时间尺度最优切负荷解提交负荷控制中心站，当负荷控制中心站判定不同时间尺度安全稳定约束低至动作门槛时，实时向子站下发决策结果，子站结合就地频率，对负控终端实时采集的用户所有380V分支回路功率进行负荷精准控制。
虽然本发明已以较佳实施例公开如上，但实施例并不是用来限定本发明的。在不脱离本发明之精神和范围内，所做的任何等效变化或润饰，同样属于本发明之保护范围。因此本发明的保护范围应当以本申请的权利要求所界定的内容为标准。
Claims (4)
1.一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法，其特征在于，包括如下步骤：
1)实时测量当地电网电气量，在故障发生到实测结束时刻的给定时间窗口内，获取系统频率及母线电压的动态响应曲线，确定切负荷措施动作后导致的切负荷后常规经济损失与电力安全事故责任代价；
2)将故障后不同时间尺度的安全稳定问题折算成广义经济损失，确定广义经济损失代价，所述不同时间尺度包括毫秒级、秒级、分钟级；
3)综合考虑切负荷后常规经济损失、电力安全事故责任代价及广义经济损失代价，得到紧急负荷控制数学模型，并采用原对偶内点法进行不同时间尺度最优切负荷求解，求解次数由不同时间尺度安全约束条件是否满足决定，且下一次模型求解的初值为上一次模型最优解；
4)将不同时间尺度最优切负荷解提交负荷控制中心站，当负荷控制中心站判定不同时间尺度安全稳定约束低至动作门槛时，实时向子站下发决策结果，子站结合就地频率，对负控终端实时采集的用户所有380V分支回路功率进行负荷精准控制。
2.根据权利要求1所述的考虑经济性及多时间尺度安全稳定性的精准负荷控制方法，其特征在于，所述步骤1)中，切负荷措施动作后导致的切负荷后常规经济损失与电力安全事故责任代价，按下式表示：
FJ＝FJ1+FJ2
式中：FJ为切负荷后常规经济损失与电力安全事故责任代价之和，FJ1为稳控系统切负荷措施动作导致的切负荷后常规经济损失，FJ2为根据不同等级电网导致相应电力安全事故的责任代价。
3.根据权利要求2所述的考虑经济性及多时间尺度安全稳定性的精准负荷控制方法，其特征在于，所述步骤2)中，具体包括以下步骤：
2-1)根据暂态功角稳定性、暂态频率安全性和暂态电压安全性将毫秒级暂态功角、电压、频率稳定情况折算成暂态广义经济损失，如下式所示：
Fms＝αFms,δ+βFms,V+γFms,f
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
<mo>,</mo>
<mi>&delta;</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
<mo>,</mo>
<mi>V</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mi>t</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>m</mi>
<mi>s</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msup>
<mi>t</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<msup>
<mi>t</mi>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>t</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
<msup>
<msub>
<mi>t</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
</mrow>
式中：Fms为将毫秒级暂态功角、电压、频率稳定情况折算成的暂态广义经济损失，Fms,δ、Fms,V、Fms,f分别为通过时域仿真得到故障仿真轨迹并从中提取出暂态功角、暂态电压、暂态频率的量化信息的暂态安全稳定性指标，α、β、γ分别为将故障后系统暂态功角、暂态电压、暂态频率指标换算成相应广义经济损失的折算系数；Adec为故障后一定仿真时间内减速面积，Ainc为加速面积；V、f分别为实际电压、频率仿真轨迹，Vcr和tcr为预设的电压偏移二元表(Vcr,tcr)的元素，fcr和tcr'为预设的频率偏移二元表(fcr,tcr')的元素，其中Vcr、fcr分别为电压偏移、频率偏移门槛值，t、t'为积分起始时间，分别取为暂态过程中实际电压曲线与Vcr首次交叉的时间、频率仿真曲线与fcr首次交叉的时间，tcr、tcr'分别为电压、频率偏移最大容许时间，VN为系统的额定电压，fN为系统的额定频率；
2-2)根据事故后系统频率恢复及电压恢复情况将秒级电压、频率稳定情况折算成准稳态广义经济损失，如下式所示：
Fs＝λFs,V+μFs,f
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>V</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>V</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<mover>
<mi>f</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
式中：Fs为将秒级电压、频率稳定情况折算成的准稳态广义经济损失，Fs,V、Fs,f分别为满足故障后系统暂态稳定并从准稳态仿真轨迹提取系统电压、频率的量化信息的准稳态安全稳定性指标，λ、μ为将故障后系统恢复电压、恢复频率指标换算成相应广义经济损失的折算系数；k为不满足准稳态恢复电压的薄弱节点总数，Vi为第i个节点的恢复电压，为恢复电压下限；fs为系统恢复频率，为恢复频率下限；
2-3)根据事故后系统安全备用、联络断面稳定情况将分钟级安全备用不足、大区联络线超用稳定情况折算成稳态广义经济损失，如下式所示：
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>,</mo>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>S</mi>
<msub>
<mi>S</mi>
<mi>N</mi>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中：Fmin为将分钟级安全备用不足、大区联络线超用稳定情况折算成的稳态广义经济损失，Fmin,s、Fmin,p分别为满足故障后系统暂态稳定及准稳态稳定并从系统稳态恢复中提取系统安全备用、联络断面稳定情况的量化信息的稳态安全指标，χ分别为将系统稳态安全备用、联络断面指标换算成相应广义经济损失的折算系数；S为事故后系统恢复安全备用，SN为满足系统安全要求的备用水平；m为联络断面过载的线路总数，Pj为第j个过载线路的稳态恢复功率，PjN为第j个过载线路的额定功率。
4.根据权利要求2所述的考虑经济性及多时间尺度安全稳定性的精准负荷控制方法，其特征在于，所述步骤3)中，紧急负荷控制数学模型表示如下：
minF＝FJ+Fms+Fs+Fmin
式中：F为将切负荷后常规经济损失与电力安全事故责任代价、暂态广义经济损失、准稳态广义经济损失、稳态广义经济损失进行加权的和；
紧急负荷控制数学模型的约束条件包括：暂态安全稳定性临界值、准稳态电压及频率下限、系统安全备用及联络断面约束以及各切负荷点的减负荷上限约束。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710774974.9A CN107516903B (zh) | 2017-08-31 | 2017-08-31 | 一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710774974.9A CN107516903B (zh) | 2017-08-31 | 2017-08-31 | 一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107516903A true CN107516903A (zh) | 2017-12-26 |
| CN107516903B CN107516903B (zh) | 2020-08-14 |
Family
ID=60725025
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710774974.9A Active CN107516903B (zh) | 2017-08-31 | 2017-08-31 | 一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107516903B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108551177A (zh) * | 2018-04-02 | 2018-09-18 | 河海大学 | 基于灵敏度分析的直流受端系统暂态切负荷控制优化方法 |
| CN108683196A (zh) * | 2018-06-12 | 2018-10-19 | 烟台职业学院 | 一种电力系统负荷调节控制方法 |
| CN109038594A (zh) * | 2018-08-13 | 2018-12-18 | 深圳供电局有限公司 | 一种稳控装置切负荷定值的分析方法及系统 |
| CN109787244A (zh) * | 2019-03-18 | 2019-05-21 | 重庆大学 | 一种确定配电网低压加速减载总量的方法 |
| CN113346476A (zh) * | 2021-04-23 | 2021-09-03 | 浙江华云电力工程设计咨询有限公司 | 一种多时间尺度控制融合的负荷终端装置及其工作方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102664416A (zh) * | 2012-05-16 | 2012-09-12 | 中国电力科学研究院 | 一种切负荷导致电力安全事故风险的防控方法 |
| CN103440401A (zh) * | 2013-07-31 | 2013-12-11 | 南京南瑞集团公司 | 提高电网输电能力的紧急控制措施风险量化评估方法 |
| CN103532142A (zh) * | 2013-10-24 | 2014-01-22 | 广州电力设计院 | 电力系统精细化稳定控制方法和系统 |
| US20160099567A1 (en) * | 2014-10-02 | 2016-04-07 | Mitsubishi Electric Research Laboratories, Inc. | Dynamic and Adaptive Configurable Power Distribution System |
-
2017
- 2017-08-31 CN CN201710774974.9A patent/CN107516903B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102664416A (zh) * | 2012-05-16 | 2012-09-12 | 中国电力科学研究院 | 一种切负荷导致电力安全事故风险的防控方法 |
| CN103440401A (zh) * | 2013-07-31 | 2013-12-11 | 南京南瑞集团公司 | 提高电网输电能力的紧急控制措施风险量化评估方法 |
| CN103532142A (zh) * | 2013-10-24 | 2014-01-22 | 广州电力设计院 | 电力系统精细化稳定控制方法和系统 |
| US20160099567A1 (en) * | 2014-10-02 | 2016-04-07 | Mitsubishi Electric Research Laboratories, Inc. | Dynamic and Adaptive Configurable Power Distribution System |
Non-Patent Citations (2)
| Title |
|---|
| 刘福锁 等: "计及电力安全事故责任的稳定控制系统风险管理", 《电力系统自动化》 * |
| 李兆伟 等: "计及电力安全事故风险的电网动态分区与切负荷协调控制研究", 《电力系统保护与测控》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108551177A (zh) * | 2018-04-02 | 2018-09-18 | 河海大学 | 基于灵敏度分析的直流受端系统暂态切负荷控制优化方法 |
| CN108683196A (zh) * | 2018-06-12 | 2018-10-19 | 烟台职业学院 | 一种电力系统负荷调节控制方法 |
| CN109038594A (zh) * | 2018-08-13 | 2018-12-18 | 深圳供电局有限公司 | 一种稳控装置切负荷定值的分析方法及系统 |
| CN109787244A (zh) * | 2019-03-18 | 2019-05-21 | 重庆大学 | 一种确定配电网低压加速减载总量的方法 |
| CN109787244B (zh) * | 2019-03-18 | 2021-07-30 | 重庆大学 | 一种确定配电网低压加速减载总量的方法 |
| CN113346476A (zh) * | 2021-04-23 | 2021-09-03 | 浙江华云电力工程设计咨询有限公司 | 一种多时间尺度控制融合的负荷终端装置及其工作方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107516903B (zh) | 2020-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107516903A (zh) | 一种考虑经济性及多时间尺度安全稳定性的精准负荷控制方法 | |
| CN104882905B (zh) | 一种考虑暂态安全约束的新能源接纳能力评估方法 | |
| CN103384068B (zh) | 电力系统暂态安全稳定紧急控制在线策略优化计算方法 | |
| CN107317345A (zh) | 一种电解类负荷参与孤立电网频率控制的方法 | |
| CN108122068A (zh) | 一种配电网风险规划方法及系统 | |
| US11201473B1 (en) | Coordinated control system and method of wind turbine and STATCOM for suppressing unbalanced voltage in dispersed wind farm | |
| CN106230024B (zh) | 含双馈风电场的电力系统停电风险计算方法 | |
| CN102868164B (zh) | 一种低频电压减载联动协调控制方法 | |
| CN102882229A (zh) | 风电场动态电压自动控制系统 | |
| CN104778550A (zh) | 一种基于实时运行数据的电网质量分析方法 | |
| CN101431238A (zh) | 一种计及“n-1”准则的配电系统供电能力评估方法 | |
| CN103337904A (zh) | 基于多源数据的电网电压越限和设备过载紧急控制方法 | |
| CN107194574A (zh) | 一种基于负荷损失的电网安全风险评估方法 | |
| CN107181253A (zh) | 一种基于电网动态可靠性概率指标的电网规划方法 | |
| CN104505855B (zh) | 一种风电场并网特性在线评价系统 | |
| CN106875287A (zh) | 一种电力系统分区供电能力确定方法 | |
| CN112865083B (zh) | 一种含可再生能源电力系统的实时运行状态评估方法 | |
| CN106026108A (zh) | 一种考虑电网安全事故风险的稳控系统切负荷方法 | |
| CN104283222A (zh) | 一种地区电网无功电压控制系统 | |
| CN103107544A (zh) | 一种海上油田群电网紧急控制的在线预决策方法 | |
| CN104993466B (zh) | 一种适用于交直流电网的连锁故障快动态仿真方法 | |
| Ballal et al. | Improvements in Existing System Integrity Protection Schemes Under Stressed Conditions by Synchrophasor Technology—Case Studies | |
| CN108092272A (zh) | 一种基于渐消卡尔曼滤波的电压稳定在线监测方法 | |
| CN104732457A (zh) | 一种基于动态ace的区域电网一次调频能力评价方法 | |
| Novosel et al. | Practical protection and control strategies during large power-system disturbances |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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 |