CN105262101B - 基于极限运行方式的无功电压校核方法 - Google Patents
基于极限运行方式的无功电压校核方法 Download PDFInfo
- Publication number
- CN105262101B CN105262101B CN201510781627.XA CN201510781627A CN105262101B CN 105262101 B CN105262101 B CN 105262101B CN 201510781627 A CN201510781627 A CN 201510781627A CN 105262101 B CN105262101 B CN 105262101B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- voltage
- delta
- mtd
- 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
- 238000000034 method Methods 0.000 title claims abstract description 84
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims description 22
- 238000012937 correction Methods 0.000 claims description 19
- 238000002347 injection Methods 0.000 claims description 16
- 239000007924 injection Substances 0.000 claims description 16
- 230000001939 inductive effect Effects 0.000 claims description 6
- 238000004804 winding Methods 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000011161 development Methods 0.000 description 9
- 239000011800 void material Substances 0.000 description 3
- 238000000205 computational method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Classifications
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Abstract
本发明属于电网无功电压控制技术领域，尤其涉及一种基于极限运行方式的无功电压校核方法。通过读入电网数据信息，形成四种运行方式，比较电压计算结果选择极限运行方式；选择一个电压越限最严重的厂站；检查该厂站变压器两侧或三侧电压，是否同时越电压上限Umax，或电压下限Umin；设定补偿因子ε，计算该厂站无功应补偿量QC或QL；判断无功应补偿量QC/QL与无功补偿实际已投入量QCr/QLr之和，是否超出无功补偿配置容量QCmax/QLmax；经上述调整后，重新计算电压。该流程循环进行，直至未有电压越限。本发明通过计算极限运行方式，可有效校核无功补偿设备配置现状是否满足相关规定及电压调整需要，指导电网运行决策；填补了无功电压校核方法的空白，贴近电网运行实际，适合在线应用。
Description
技术领域
本发明属于电网无功电压控制技术领域，尤其涉及一种基于极限运行方式的无功电压校核方法。
背景技术
在目前的电力系统中，随着负荷不断增长，特高压等各项新增工程不断开展，无功电压控制面临各种新的挑战。为保障电网电压运行质量，有必要对电网无功电压情况进行校核，掌握电网无功补偿设备配置情况，校核是否满足相关规定和电压调整需要，对无功优化、节能降损和AVC系统建设奠定基础。无功电压校核是目前建设智能电网的一项不可或缺的工作，对电网的安全、稳定、经济运行具有重要意义。
目前，对无功电压控制的研究大多集中在无功优化方面,而对无功电压校核方法的研究较为缺乏，并未提出一种切实可行的电网无功电压校核方法。虽然各地方电网可能开展过无功电压校核工作，但并未对无功电压校核方法进行提炼总结，该项工作纯靠技术人员经验，目前并无方法可循。
因此，发明一种切实可行的、易于理解和实现、并贴近电网运行实际的无功电压校核方法具有重要意义。
发明内容
为克服现有技术中存在的缺陷，本发明提供一种基于极限运行方式的无功电压校核方法，其目的是为了提供一种可以填补无功电压校核方法的空白，为电网运行决策提供指导的无功电压校核方法。
为达到上述目的，本发明所采用的技术方案是：
基于极限运行方式的无功电压校核方法，包括以下步骤：
步骤1：读入所研究电网数据信息。
所述数据信息包括以下方面：所研究电网范围内的厂站夏季最大、最小负荷、各机组夏季最大、最小出力，以及厂站冬季最大、最小负荷、机组冬季最大、最小出力；所研究电网各厂站无功装置配置情况：无功补偿装置配置情况QCmax、QLmax，变压器档位分接头设置情况；实际无功运行情况：目前无功补偿装置实际投退情况QCr、QLr，目前电压器分接头所处档位信息；无功控制限制：各电压等级上下限要求Umax、Umin，各发电机功率因数限制；其他用于潮流计算的电网数据。
步骤2：从读入的数据信息中形成四种运行方式，分别计算电压，并选择极限运行方式。
所述运行方式是指电力系统运行方式，根据本系统实际情况，合理使用资源，使整个系统在安全、优质、经济运行情况下的决策。运行方式表征电网的一种运行状态。本专利中的运行方式侧重于所研究电网的发电机组出力、运行负荷、无功补偿装置配置等情况。
所述形成的四种运行方式包括：
①夏季最大负荷，夏季最小出力，实际无功运行情况
②冬季最大负荷，冬季最小出力，实际无功运行情况
③夏季最小负荷，夏季最大出力，实际无功运行情况
④冬季最小负荷，冬季最大出力，实际无功运行情况
所述计算电压，是指用高斯-赛德尔法对四种运行方式进行潮流计算，得到各母线电压，采用公式其中，Yii、Yij分别为节点导纳矩阵的导纳值，Yii(i＝1,2,…,n)为节点导纳矩阵的对角元，称为自导纳；Yij(i＝1,2,…,n；j＝1,2,…,n；i≠j)为节点导纳矩阵的非对角元，称为互导纳；分别为节点i、j的电压，初值均设为1∠0°；Pi、Qi分别为i节点注入的有功功率、无功功率幅值；jQi为i节点注入功率的虚部，表征注入无功功率幅值和相位；-jQi是由于公式推导时有反相，所以前为负号；n为节点导纳矩阵的阶数。
所述选择极限运行方式是指，在方式①和②中，选择电压水平最低的情况，作为低电压极限运行方式；在方式③和④中，选择电压水平最高的情况，作为高电压极限运行方式。
步骤3：分别对两种极限运行方式的电压计算结果进行排序，并与电压运行规定上下限进行比较，选择一个电压越限最严重的厂站。
步骤4：检查该厂站变压器两侧或三侧电压，是否同时越电压上限Umax，或者同时越电压下限Umin。若电压同时越上限，或同时越下限，则进行下一步操作；否则，在可调整范围内，调整变压器分接头，之后重新计算该方式电压，回至步骤3。
所述检查该厂站变压器两侧或三侧电压是指，若该厂站为两绕组变压器，则检查变压器两侧电压；若该厂站为三绕组变压器，则检查变压器三侧电压。
步骤5：设定补偿因子ε，计算该厂站无功应补偿量QC或QL。
所述计算该厂站无功应补偿量采用公式QC＝(Umin-U)ε或QL＝(U-Umax)ε，其中，QC为容性无功应补偿量；QL为感性无功应补偿量；U为该厂站电压值；Umin为电压运行下限；Umax为电压运行下限；ε为设定的补偿因子，通常取1-10，可通过实际电网情况得到。若该厂站电压U低于运行电压下限Umin，则利用公式QC＝(Umin-U)ε计算容性无功应补偿量QC；若该厂站电压U高于运行电压上限Umax，则利用公式QL＝(U-Umax)ε计算感性无功应补偿量QL。
步骤6：判断无功应补偿量QC/QL与无功补偿实际已投入量QCr/QLr之和，是否超出无功补偿配置容量QCmax/QLmax。
1)若QC+QCr＞QCmax或QL+QLr＞QLmax成立，判断该变电站附近是否有发电机可在功率因数限制内调整无功出力，帮助调节电压。
若这样的发电机存在，则在功率因数限制内调节发电机无功出力，并投入现有无功补偿配置容量QCmax/QLmax，之后重新计算该方式电压，回至步骤3；
若不存在这样的发电机，投入容量为无功应补偿量QC/QL的无功补偿，并输出记录该变电站名称，以及该变电站无功应补偿量QC/QL、无功补偿已投入量QCr/QLr、无功补偿配置容量QCmax/QLmax。之后进行下一步操作。
2)否则，投入最接近无功应补偿量QC/QL的无功补偿容量，进行下一步操作。
步骤7：经上述调整后，重新计算电压，之后回至步骤3；
上述步骤循环进行，一旦步骤3判断未有电压越限时，则退出程序。
所述步骤2中：计算电压也可用牛顿-拉夫逊法等各种计算方法对四种运行方式进行潮流计算，得到各母线电压，不局限于一种潮流计算方法；
牛顿拉夫逊极坐标法中，潮流方程的修正方程为式(1)，由式(1)可推导出电压计算公式(2)，式(2)各项的具体计算公式为(3)-(10)；
其中，(2)式中，ΔP为注入有功功率的修正值；ΔQ为注入无功功率的修正值；H、N、J、L为雅克比矩阵的元素；Δδ为电压相角修正值，Δδi表征节点i的电压相角修正值，计算公式为(9)，为i节点电压相角k+1、k次迭代值，一般初值为0；ΔU为电压幅值修正值，ΔUi表征节点i的电压幅值修正值，计算公式为(10)，为i节点电压幅值k+1、k次迭代值，一般初值为1；
(3)、(4)式中ΔPi、ΔQi分别为节点i注入有功功率、无功功率的修正值；Pi、Qi为节点i注入有功功率、无功功率实际值，为常数；Pis、Qis为节点i注入有功功率、无功功率的计算值；Ui、Uj为节点i、j的电压幅值；δij为节点i、j的电压相角差，δij＝δi-δj，δi、δj为节点i、j的电压相角；Yij＝Gij+jBij，Yij为节点导纳矩阵第i行、第j列的元素，Gij、Bij为Yij的实部、虚部；
(5)-(8)式中，Hij、Nij、Jij、Lij分别雅克比矩阵H、N、J、L部分第i行、第j列的元素。
本发明的优点及有益效果是：
本发明提供的方法，通过对电压极限运行方式进行计算，有效地校核了电网无功补偿设备配置情况是否满足相关规定和电压调整需要，从而对电网运行决策提供指导；另外，本方法填补了无功电压校核方法的空白，易于理解和实现、并贴近电网运行实际，适合在线应用。
附图说明
图1是本发明基于极限运行方式的无功电压校核方法流程图；
图2是本发明实施例提供的电网接线图。
具体实施方式
下面结合附图，对优选实施例作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。
实施例1：
图1是基于极限运行方式的无功电压校核方法流程图。如图1所示，本发明提出的基于极限运行方式的无功电压校核方法包括：
步骤1：读入所研究电网数据信息。
所述数据信息包括所研究电网范围内的厂站夏季最大、最小负荷、各机组夏季最大、最小出力，以及厂站冬季最大、最小负荷、机组冬季最大、最小出力；所研究电网各厂站无功装置配置情况：无功补偿装置配置情况QCmax、QLmax，变压器档位分接头设置情况；实际无功运行情况：目前无功补偿装置实际投退情况QCr、QLr，目前电压器分接头所处档位信息；无功控制限制：各电压等级上下限要求Umax、Umin，各发电机功率因数限制；用于潮流计算的其他电网数据。
可将上述电网数据信息输入电力系统分析综合程序(PSASP)，形成仿真数据库，用于开展以下无功电压校核工作。
步骤2：从读入的数据信息中形成四种运行方式，分别计算电压，并选择极限运行方式。
所述运行方式是指电力系统运行方式，根据本系统实际情况，合理使用资源，使整个系统在安全、优质、经济运行情况下的决策。运行方式表征电网的一种运行状态。本专利中的运行方式侧重于所研究电网的发电机组出力、运行负荷、无功补偿装置配置等情况。
所述形成的四种运行方式包括：
①夏季最大负荷，夏季最小出力，实际无功运行情况
②冬季最大负荷，冬季最小出力，实际无功运行情况
③夏季最小负荷，夏季最大出力，实际无功运行情况
④冬季最小负荷，冬季最大出力，实际无功运行情况
所述计算电压，是指用高斯-赛德尔法对四种运行方式进行潮流计算，得到各母线电压。采用高斯-赛德尔法，电压计算公式为：
其中，Yii、Yij分别为节点导纳矩阵的导纳值，Yii(i＝1,2,…,n)为节点导纳矩阵的对角元，称为自导纳；Yij(i＝1,2,…,n；j＝1,2,…,n；i≠j)为节点导纳矩阵的非对角元，称为互导纳；节点导纳矩阵可通过电网数据信息得到；分别为节点i、j的电压，初值均设为1∠0°；Pi、Qi分别为i节点注入的有功功率、无功功率幅值；jQi为i节点注入功率的虚部，表征注入无功功率幅值和相位；-jQi是由于公式推导时有反相，所以前为负号；n为节点导纳矩阵的阶数，也等于网络中除参考节点外的节点数。
采用电力系统分析综合程序(PSASP)搭建上述四种运行方式，进行潮流计算，可方便得到并输出电压计算结果。
在方式①和②中，选择电压水平最低的情况，作为低电压极限运行方式；
在方式③和④中，选择电压水平最高的情况，作为高电压极限运行方式。
步骤3：分别对两种极限运行方式的电压计算结果进行排序，并与电压运行规定上下限进行比较，选择一个电压越限最严重的厂站。
PSASP中电压计算结果可用Excel输出，并且越限电压会用星号标示。利用Excel排序功能，可方便排序，表格最上方或最大方即为越限最严重的厂站。
步骤4：检查该厂站变压器两侧或三侧电压，是否同时越电压上限Umax，或者同时越电压下限Umin。
若该厂站为两绕组变压器，则检查变压器两侧电压；若该厂站为三绕组变压器，则检查变压器三侧电压。
若电压同时越上限，或同时越下限，则进行下一步操作；否则，在可调整范围内，调整变压器分接头，之后重新计算该方式电压，回至步骤3。
PSASP以Excel形式输出各电压等级电压计算结果，可方便筛选同一厂站不同电压等级侧电压，可用于检查是否同时越上限，或越下限。
采用PSASP可方便查看并调整变压器分接头，重新进行潮流计算即可得分接头调整后的电压计算结果。
步骤5：设定补偿因子ε，计算该厂站无功应补偿量QC或QL。
计算该厂站无功应补偿量采用公式
QC＝(Umin-U)ε (2)
QL＝(U-Umax)ε (3)
其中，QC为容性无功应补偿量；QL为感性无功应补偿量；U为该厂站电压值；Umin为电压运行下限；Umax为电压运行下限；ε为设定的补偿因子，通常取1-10，可通过实际电网情况得到。若该厂站电压U低于运行电压下限Umin，则利用公式(2)计算容性无功应补偿量QC；若该厂站电压U高于运行电压上限Umax，则利用公式(3)计算感性无功应补偿量QL。
步骤6：判断无功应补偿量QC/QL与无功补偿实际已投入量QCr/QLr之和，是否超出无功补偿配置容量QCmax/QLmax。
1)若QC+QCr＞QCmax或QL+QLr＞QLmax成立，判断该变电站附近是否有发电机可在功率因数限制内调整无功出力，帮助调节电压。
若这样的发电机存在，则在功率因数限制内调节发电机无功出力，并投入现有无功补偿配置容量QCmax/QLmax，之后重新计算该方式电压，回至步骤3；
若不存在这样的发电机，投入容量为无功应补偿量QC/QL的无功补偿，并输出记录该变电站名称，以及该变电站无功应补偿量QC/QL、无功补偿已投入量QCr/QLr、无功补偿配置容量QCmax/QLmax。之后进行下一步操作。
2)否则，投入最接近无功应补偿量QC/QL的无功补偿容量，进行下一步操作。
采用PSASP可方便对厂站无功补偿装置进行投退，并可以按照需要修改无功补偿容量。
步骤7：经上述调整后，重新计算电压，之后回至步骤3。
上述步骤循环进行，一旦步骤3判断未有电压越限时，则退出程序。输出记录的变电站，即为校核得到的无功补偿容量不足的变电站。输出的无功应补偿量QC/QL、无功补偿已投入量QCr/QLr、无功补偿配置容量QCmax/QLmax，可计算该变电站应配置的无功总容量(QC+QCr、QL+QLr)，以及需增设的无功容量(QC+QCr-QCmax、QL+QLr-QLmax)。
实施例2：
下面以2015年内蒙古通辽地区电网作为本发明的一个实施例，对本发明的发明内容做进一步说明。
内蒙古通辽地区电网接线图如图2所示。其中通辽电网电压控制要求：500kV母线电压控制在520-540kV之间，220kV变电站母线电压控制在220-237之间；风电场升压变电所有220kV母线电压需控制在225-235kV之间。
通过对四种运行方式进行计算，选出通辽电网电压极限运行方式分别为：
夏季最大负荷，夏季机组最小出力，作为低电压极限运行方式；
夏季最小负荷，夏季机组最大出力，作为高电压极限运行方式。
对两种方式电压计算结果进行排序，并与电压运行规定上下限进行比较，可得夏小负荷的高电压极限运行方式中，各母线电压均未越限；夏大负荷的低电压极限运行方式中，电压存在越下限情况。
在夏大负荷方式下，220kV开鲁变、开鲁开发区变、奈曼变、红星变和和谊变电压均低于220kV，其中奈曼变电压216kV，为通辽全网最低电压。奈曼接入风电母线电压均低于220kV，难以满足正常运行。
检查发现奈曼变变压器两侧电压同时越下限，无需调整变压器分接头。设定补偿因子ε＝5，Umin此时为220kV，计算得奈曼变容性无功应补偿量QC＝20Mvar。目前运行方式中，奈曼变配置的21.3Mvar电容已全部投入，并且奈曼变附近的吉利、萨日、永兴风场功率因数已达最低限值0.98。因此，奈曼变再增投20Mvar电容，重新计算该方式电压。
奈曼变增投20Mvar电容后，奈曼变、红星变、和谊变电压均上升至220kV以上。开鲁开发区变、开鲁变电压仍低于220kV，其中开鲁开发区变电压最低，为218kV。判断无需调整变压器分接头。计算得开鲁开发区变容性无功应补偿量QC＝10Mvar。目前开鲁开发区未配置无功补偿装置。开鲁开发区变附近的北清河、木仁、义和、塔拉风场功率因数已达最低限值0.98。因此，开鲁开发区变增投10Mvar电容，重新计算该方式电压。
经上述电压调整后，开鲁、奈曼地区电压均提高至220kV以上。奈曼变和开鲁变接入风电的母线电压，以及红星变、开鲁变、和谊变66kV母线电压均提高至正常运行范围内。开鲁、奈曼地区电压偏低问题得以解决。输出的电压信息如表1：输出的无功补偿信息表。
通过本方法校核可知，开鲁、奈曼地区存在电压偏低问题，奈曼变、开鲁开发区变无功补偿容量不足。需在奈曼变新增20Mvar的容性无功补偿，开鲁开发区变新增10Mvar的容性无功补偿。
实施例3：
本发明实施例1中的步骤2中，计算电压也可用牛顿-拉夫逊法等各种计算方法对四种运行方式进行潮流计算，得到各母线电压，不局限于一种潮流计算方法；
牛顿拉夫逊极坐标法中，潮流方程的修正方程为式(1)，由式(1)可推导出电压计算公式(2)，式(2)各项的具体计算公式为(3)-(10)；
其中，(2)式中，ΔP为注入有功功率的修正值；ΔQ为注入无功功率的修正值；H、N、J、L为雅克比矩阵的元素；Δδ为电压相角修正值，Δδi表征节点i的电压相角修正值，计算公式为(9)，为i节点电压相角k+1、k次迭代值，一般初值为0；ΔU为电压幅值修正值，ΔUi表征节点i的电压幅值修正值，计算公式为(10)，为i节点电压幅值k+1、k次迭代值，一般初值为1；
(3)、(4)式中ΔPi、ΔQi分别为节点i注入有功功率、无功功率的修正值；Pi、Qi为节点i注入有功功率、无功功率实际值，为常数；Pis、Qis为节点i注入有功功率、无功功率的计算值；Ui、Uj为节点i、j的电压幅值；δij为节点i、j的电压相角差，δij＝δi-δj，δi、δj为节点i、j的电压相角；Yij＝Gij+jBij，Yij为节点导纳矩阵第i行、第j列的元素，Gij、Bij为Yij的实部、虚部；
(5)-(8)式中，Hij、Nij、Jij、Lij分别雅克比矩阵H、N、J、L部分第i行、第j列的元素。
其它同实施例1。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
表1：输出的无功补偿信息表。
Claims (2)
1.基于极限运行方式的无功电压校核方法，其特征在于包括以下步骤：
步骤1：读入所研究电网数据信息；
所述数据信息包括以下方面：所研究电网范围内的厂站夏季最大、最小负荷、各机组夏季最大、最小出力，以及厂站冬季最大、最小负荷、机组冬季最大、最小出力；所研究电网各厂站无功装置配置情况：无功补偿装置配置情况QC max、QL max，变压器档位分接头设置情况；实际无功运行情况：目前无功补偿装置实际投退情况QCr、QLr，目前变压器分接头所处档位信息；无功控制限制：各电压等级上下限要求Umax、Umin，各发电机功率因数限制；其他用于潮流计算的电网数据；
步骤2：从读入的数据信息中形成四种运行方式，分别计算电压，并选择极限运行方式；
所述运行方式是指电力系统运行方式，根据本系统实际情况，合理使用资源，使整个系统在安全、优质、经济运行情况下的决策；运行方式表征电网的一种运行状态；这里的运行方式侧重于所研究电网的发电机组出力、运行负荷、无功补偿装置配置情况；
所述形成的四种运行方式包括：
①夏季最大负荷，夏季最小出力，实际无功运行情况；
②冬季最大负荷，冬季最小出力，实际无功运行情况；
③夏季最小负荷，夏季最大出力，实际无功运行情况；
④冬季最小负荷，冬季最大出力，实际无功运行情况；
所述计算电压，用高斯-赛德尔法对四种运行方式进行潮流计算，得到各母线电压，采用公式其中，Yii、Yij分别为节点导纳矩阵的导纳值，Yii(i＝1,2,…,n)为节点导纳矩阵的对角元，称为自导纳；Yij(i＝1,2,…,n；j＝1,2,…,n；i≠j)为节点导纳矩阵的非对角元，称为互导纳；分别为节点i、j的电压，初值均设为1∠0°；Pi、Qi分别为i节点注入的有功功率、无功功率幅值；jQi为i节点注入功率的虚部，表征注入无功功率幅值和相位；-jQi是由于公式推导时有反相，所以前为负号；n为节点导纳矩阵的阶数；
所述选择极限运行方式是指，在方式①和②中，选择电压水平最低的情况，作为低电压极限运行方式；在方式③和④中，选择电压水平最高的情况，作为高电压极限运行方式；
步骤3：分别对两种极限运行方式的电压计算结果进行排序，并与电压运行规定上下限进行比较，选择一个电压越限最严重的厂站；
步骤4：检查该厂站变压器两侧或三侧电压，是否同时越电压上限Umax，或者同时越电压下限Umin；若电压同时越上限，或同时越下限，则进行下一步操作；否则，在可调整范围内，调整变压器分接头，之后重新计算该方式电压，回至步骤3；
所述检查该厂站变压器两侧或三侧电压是指，若该厂站为两绕组变压器，则检查变压器两侧电压；若该厂站为三绕组变压器，则检查变压器三侧电压；
步骤5：设定补偿因子ε，计算该厂站无功应补偿量QC或QL；
所述计算该厂站无功应补偿量采用公式QC＝(Umin-U)ε或QL＝(U-Umax)ε，其中，QC为容性无功应补偿量；QL为感性无功应补偿量；U为该厂站电压值；Umin为电压运行下限；Umax为电压运行上限；ε为设定的补偿因子，取1-10，可通过实际电网情况得到；若该厂站电压U低于运行电压下限Umin，则利用公式QC＝(Umin-U)ε计算容性无功应补偿量QC；若该厂站电压U高于运行电压上限Umax，则利用公式QL＝(U-Umax)ε计算感性无功应补偿量QL；
步骤6：判断无功应补偿量QC/QL与无功补偿实际已投入量QCr/QLr之和，是否超出无功补偿配置容量QC max/QL max；
1)若QC+QCr＞QC max或QL+QLr＞QL max成立，判断变电站附近是否有发电机可在功率因数限制内调整无功出力，帮助调节电压；
若这样的发电机存在，则在功率因数限制内调节发电机无功出力，并投入现有无功补偿配置容量QC max/QL max，之后重新计算该方式电压，回至步骤3；
若不存在这样的发电机，投入容量为无功应补偿量QC/QL的无功补偿，并输出记录变电站名称，以及变电站无功应补偿量QC/QL、无功补偿已投入量QCr/QLr、无功补偿配置容量QC max/QL max；之后进行下一步操作；
2)否则，投入最接近无功应补偿量QC/QL的无功补偿容量，进行下一步操作；
步骤7：经上述调整后，重新计算电压，之后回至步骤3；
上述步骤循环进行，一旦步骤3判断未有电压越限时，则退出程序。
2.根据权利要求1所述的基于极限运行方式的无功电压校核方法，其特征在于：所述步骤2中：计算电压用牛顿-拉夫逊法对四种运行方式进行潮流计算，得到各母线电压；”
牛顿拉夫逊极坐标法中，潮流方程的修正方程为式(1)，由式(1)可推导出电压计算公式(2)，式(2)各项的具体计算公式为(3)-(10)；
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>P</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>Q</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>H</mi>
</mtd>
<mtd>
<mi>N</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>J</mi>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>&delta;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>U</mi>
<mo>/</mo>
<mi>U</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>&delta;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>U</mi>
<mo>/</mo>
<mi>U</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>H</mi>
</mtd>
<mtd>
<mi>N</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>J</mi>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>P</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>Q</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Delta;Q</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&delta;</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>J</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;Q</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&delta;</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<mi>B</mi>
<mi> </mi>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;Q</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Delta;&delta;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>&delta;</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&delta;</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Delta;U</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，(2)式中，ΔP为注入有功功率的修正值；ΔQ为注入无功功率的修正值；H、N、J、L为雅克比矩阵的元素；Δδ为电压相角修正值，Δδi表征节点i的电压相角修正值，计算公式为(9)，为i节点电压相角k+1、k次迭代值，初值为0；ΔU为电压幅值修正值，ΔUi表征节点i的电压幅值修正值，计算公式为(10)，为i节点电压幅值k+1、k次迭代值，初值为1；
(3)、(4)式中ΔPi、ΔQi分别为节点i注入有功功率、无功功率的修正值；Pi、Qi为节点i注入有功功率、无功功率实际值，为常数；Pis、Qis为节点i注入有功功率、无功功率的计算值；Ui、Uj为节点i、j的电压幅值；δij为节点i、j的电压相角差，δij＝δi-δj，δi、δj为节点i、j的电压相角；Yij＝Gij+jBij，Yij为节点导纳矩阵第i行、第j列的元素，Gij、Bij为Yij的实部、虚部；
(5)-(8)式中，Hij、Nij、Jij、Lij分别雅克比矩阵H、N、J、L部分第i行、第j列的元素。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510781627.XA CN105262101B (zh) | 2015-11-12 | 2015-11-12 | 基于极限运行方式的无功电压校核方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510781627.XA CN105262101B (zh) | 2015-11-12 | 2015-11-12 | 基于极限运行方式的无功电压校核方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105262101A CN105262101A (zh) | 2016-01-20 |
| CN105262101B true CN105262101B (zh) | 2018-02-02 |
Family
ID=55101657
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510781627.XA Active CN105262101B (zh) | 2015-11-12 | 2015-11-12 | 基于极限运行方式的无功电压校核方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105262101B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102522756A (zh) * | 2011-12-14 | 2012-06-27 | 华南理工大学 | 一种避免电压越限风险的电网感性无功补偿方法 |
| CN102545230A (zh) * | 2012-01-10 | 2012-07-04 | 四川省电力公司 | 一种电压越限调整方法 |
| CN102709923A (zh) * | 2012-06-14 | 2012-10-03 | 济南大学 | 基于负荷波动规律的时变限值区域电网无功优化控制方法 |
| CN104362650A (zh) * | 2014-11-14 | 2015-02-18 | 国家电网公司 | 一种考虑成本因素的电力系统无功优化方法 |
-
2015
- 2015-11-12 CN CN201510781627.XA patent/CN105262101B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102522756A (zh) * | 2011-12-14 | 2012-06-27 | 华南理工大学 | 一种避免电压越限风险的电网感性无功补偿方法 |
| CN102545230A (zh) * | 2012-01-10 | 2012-07-04 | 四川省电力公司 | 一种电压越限调整方法 |
| CN102709923A (zh) * | 2012-06-14 | 2012-10-03 | 济南大学 | 基于负荷波动规律的时变限值区域电网无功优化控制方法 |
| CN104362650A (zh) * | 2014-11-14 | 2015-02-18 | 国家电网公司 | 一种考虑成本因素的电力系统无功优化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105262101A (zh) | 2016-01-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103199542B (zh) | 一种风电场无功电压优化控制方法 | |
| CN100557918C (zh) | 基于协调约束实时更新的电厂与变电站协调电压控制方法 | |
| CN102832627B (zh) | 一种基于电力网络模型的风电场自动电压控制方法 | |
| CN107134799A (zh) | 一种基于短路比的特高压电网换流器无功配置方法 | |
| CN103515967B (zh) | 风电场无功补偿的配置方法 | |
| CN108418255A (zh) | 一种适用于含高渗透率新能源的特高压直流外送电网规划方法及系统 | |
| CN102570477B (zh) | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 | |
| CN103150606A (zh) | 一种分布式电源最优潮流优化方法 | |
| CN105490282B (zh) | 考虑微电源无功出力均衡度的微电网实时电压控制方法 | |
| CN105226665A (zh) | 通过自动厂站协调控制降低母线电压日波动率的方法 | |
| CN104158199B (zh) | 对电力系统实时状态进行无功电压优化控制的系统和方法 | |
| CN106026113A (zh) | 一种具有无功自动补偿的微电网系统的监控方法 | |
| CN103793755B (zh) | 基于原对偶内点法的混合直流输电系统最优潮流方法 | |
| CN103701132B (zh) | 一种基于有功-无功协调的控制主动配电网过电压的方法 | |
| CN101917002B (zh) | 基于特高压电网联络线运行状态的电压控制方法 | |
| CN105870968A (zh) | 一种计及系统负序电压的三相不平衡无功电压控制方法 | |
| CN103400011A (zh) | 一种新型磁控式并联电抗器的仿真建模方法 | |
| CN105375513A (zh) | 一种基于实时在线等值的110千伏风电场自动电压控制方法 | |
| CN106374520A (zh) | 一种安全互动调度的分布式电源可接纳功率评估方法 | |
| CN103346573B (zh) | 基于黄金分割云粒子群优化算法的风电系统无功规划方法 | |
| CN104392285B (zh) | 一种含混合直流输电的电力系统最优潮流获取方法 | |
| Yang et al. | Optimal coordinated voltage control of AC/DC power systems for voltage stability enhancement | |
| CN106295952A (zh) | 一种高阻抗变压器限制短路电流效果的模糊综合评价方法 | |
| CN107196343A (zh) | 一种多端柔性直流孤岛电网送端的日前电压计划生成方法 | |
| CN102904265B (zh) | 一种基于区间潮流的变电站电压无功控制方法和系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |