CN112003288A - 一种电网运行方式电压智能调整方法及装置 - Google Patents

一种电网运行方式电压智能调整方法及装置 Download PDF

Info

Publication number
CN112003288A
CN112003288A CN202010760141.9A CN202010760141A CN112003288A CN 112003288 A CN112003288 A CN 112003288A CN 202010760141 A CN202010760141 A CN 202010760141A CN 112003288 A CN112003288 A CN 112003288A
Authority
CN
China
Prior art keywords
bus
sensitivity
voltage
equal
ctr
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
CN202010760141.9A
Other languages
English (en)
Other versions
CN112003288B (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.)
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
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 China Electric Power Research Institute Co Ltd CEPRI filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN202010760141.9A priority Critical patent/CN112003288B/zh
Publication of CN112003288A publication Critical patent/CN112003288A/zh
Application granted granted Critical
Publication of CN112003288B publication Critical patent/CN112003288B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • 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/16—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
    • 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
    • H02J3/1821—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators
    • H02J3/1835—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators with stepless control
    • H02J3/1864—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators with stepless control wherein the stepless control of reactive power is obtained by at least one reactive element connected in series with a semiconductor switch
    • 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

一种电网运行方式电压智能调整方法及装置，该方法包括：首先根据容抗资源与被控母线计算容抗投切灵敏度；其次通过灵敏度组合获得被控母线与网络中电容器、电抗器无功资源的支撑关系，动态生成满足目标母线电压安全调整的预期有效灵敏度组合系数；最后通过校验获得满足目标电压调整的有效灵敏度组合系数用于电压调整控制。通过本发明实施例提供的方法及装置，解决了传统电压稳定仿真分析和无功资源配置严重依赖人工的问题，大大降低了人力成本，提高了计算效率，同时计算分析的精确度也得到了有效保证；可以为方式分析编制人员提供强有力的专业化支撑工具，提高电压控制调整的自动化智能化水平。

Description

一种电网运行方式电压智能调整方法及装置
技术领域
本发明涉及电网分析与控保领域，具体而言，涉及一种电网运行方式电压智能调整方法及装置。
背景技术
随着新能源比例提升、特高压交直流建设，电网呈现出电力电子化特征，提供电压支撑/消耗无功功率的网络元件和无功需求场景大大丰富，局部电网电压不合理现象时有发生。对于交直流混联特高压电网，存在单条特高压线路末端解环，特高压线路充电无功较大，导致线路末端电压升高而超过设备耐压水平危险。因此，在日常方式计算分析中或电网规划设计中，需要配置或确定电网低抗、低容配置及配比，以保证电网的电压运行稳定性，避免电压越限而引发系统运行安全风险。
大规模交直流电网中，常规电压调整方案除发电机出力改变、变压器有载分接开关调整外，主要依靠各类离散无功补偿装置，静态无功补偿装置有并联电容器、并联电抗器，动态无功补偿装置有静止无功补偿器(SVC)、静止同步补偿器(STATCOM)等。在获得电压稳定薄弱点及具备多类无功资源的基础上，需要进行电压安全约束快速计算以及电网无功资源控制策略的制定。
传统的电压稳定仿真分析和无功资源配置主要依靠人工，严重依赖专家经验，随着我国电网电力电子化发展和电压问题突出，针对电压调整的人工分析在诸多方面已渐显不足，主要包括大量消耗人力成本、易于引发错漏现象、难以精确调整控制误差等。
发明内容
鉴于此，本发明提出了一种电网运行方式电压智能调整方法及装置，旨在解决传统电压稳定仿真分析和无功资源配置严重依赖人工而导致人力成本高、易于引发错漏现象、难以精确调整控制误差的问题。
第一方面，本发明实施例提供了一种电网运行方式电压智能调整方法，包括：步骤101：获取基础运行方式潮流数据文件、电网容抗资源配置文件；步骤102：根据所述电网容抗资源配置文件形成可控资源矩阵MCtr，同时从所述基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定所述被控母线和所述监视母线组的电压运行限值；步骤103：基于所述基础运行方式潮流数据文件进行潮流计算并得到潮流计算结果，从所述潮流计算结果中提取所述监视母线组的当前电压值和所述被控母线的当前电压值；步骤104：对所述监视母线组的当前电压值进行限值校核，同时，基于所述被控母线的当前电压值计算所述被控母线的电压待调整值；步骤105：根据所述被控母线的电压待调整值和所述可控资源矩阵MCtr，分别进行切电容/投电抗、投电容/切电抗灵敏度计算，形成母线降压/升压灵敏度矩阵；步骤106：基于所述母线降压/ 升压灵敏度矩阵进行灵敏度组合，计算预期电压调整值，筛选生成预期有效灵敏度组合系数；步骤107：基于所述预期有效灵敏度组合系数进行容抗投退后的潮流计算以及校核，得到并输出校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制。
进一步地，所述基础运行方式潮流数据文件包括：发电机母线模型、交直流母线模型、交直流线路模型、变压器模型、电容器模型、电抗器模型、除电容器电抗器以外的其它FACTS设备模型。
进一步地，所述电网容抗资源配置文件包括：具有容抗资源的母线名称、单组电容器容量、已投电容器组数、未投电容器组数、单组电抗器容量、已投电抗器组数、未投电抗器组数。
进一步地，步骤102中根据所述电网容抗资源配置文件形成可控资源矩阵 MCtr，包括：根据所述电网容抗资源配置文件形成如下可控资源矩阵MCtr：
Figure BDA0002612837980000031
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000032
为母线i已投电容器的组数，
Figure BDA0002612837980000033
为母线 i未投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure BDA0002612837980000034
为母线i已投电抗器的组数，
Figure BDA0002612837980000035
为母线i未投电抗器的组数。
进一步地，所述步骤102中从所述基础运行方式潮流数据文件中选定母线和监视母线组，并分别设定所述被控母线和所述监视母线组的电压运行限值，包括：步骤1021：从所述基础运行方式潮流数据文件中选定1条被控母线BusC，设定所述被控母线的电压控制目标值为
Figure BDA0002612837980000036
电压控制目标值波动带宽为
Figure BDA0002612837980000037
步骤1022：从所述基础运行方式潮流数据文件中选定t条监视母线BusMi(1≤i≤t)，设定每条监视母线的电压上限值为
Figure BDA0002612837980000038
电压下限值为
Figure BDA0002612837980000039
电压监视限值调整裕度为ΔVM，同时所述t条监视母线的电压上下限值构成如下矩阵Mt：
Figure BDA00026128379800000310
进一步地，步骤104中对所述监视母线组的当前电压值进行限值校核，包括：若
Figure BDA00026128379800000311
则更新每条监视母线BusMi的电压上限值为：
Figure BDA00026128379800000312
若
Figure BDA00026128379800000313
则更新每条监视母线BusMi的电压下限值为：
Figure BDA00026128379800000314
其中，
Figure BDA00026128379800000315
为每条监视母线组BusM当前电压值，
Figure BDA00026128379800000316
为每条监视母线的电压上限值，
Figure BDA00026128379800000317
为每条监视母线的电压下限值，ΔVM为电压监视限值调整裕度。
进一步地，步骤104中基于所述被控母线的当前电压值计算所述被控母线的电压待调整值，包括：所述被控母线BusC的电压待调整值
Figure BDA0002612837980000041
其中，
Figure BDA0002612837980000042
为所述被控母线的当前电压值，
Figure BDA0002612837980000043
为所述被控母线的电压控制目标值。
进一步地，所述步骤105根据所述被控母线的电压待调整值和所述可控资源矩阵MCtr，分别进行切电容/投电抗、投电容/切电抗灵敏度计算，形成母线降压/升压灵敏度矩阵，包括：步骤1051：判断所述被控母线的电压待调整方向：步骤1052：若ΔVC＜0，根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵；步骤1053：若ΔVC＞0，根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵。
进一步地，所述根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵，包括：步骤10521：基于所述可控资源矩阵MCtr抽取信息得到母线降压资源矩阵MCtr-：
Figure BDA0002612837980000044
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000045
为母线i已投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure BDA0002612837980000046
为母线i未投电抗器的组数；
步骤10522：设置i＝1；步骤10523：根据MCtr-中可控资源，在所述基础运行方式潮流数据文件的母线BusRi退出单组电容器容量即QRiC(MVar)，执行潮流计算；步骤10524：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000047
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure BDA0002612837980000048
其中，
Figure BDA0002612837980000049
为被控母线在步骤104中的电压值，QRiC(1≤i≤m)为被控母线i单组电容器的容量；步骤10525：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000051
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure BDA0002612837980000052
步骤10526：根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电抗器容量即QR1L(MVar)，执行潮流计算；步骤10527：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000053
计算得到母线BusRi对母线 BusC的投电抗灵敏度
Figure BDA0002612837980000054
步骤10528：若潮流不收敛，则根据 MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000055
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure BDA0002612837980000056
步骤 10529：若i＜m，则i＝i+1，返回步骤10523；若i＝m，基于所述切电容灵敏度SCi-、所述投电抗灵敏度SLi+，形成母线降压灵敏度矩阵SV-：
Figure BDA0002612837980000057
进一步地，所述根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵，包括：步骤10531：基于所述可控资源矩阵MCtr抽取信息得到母线升压资源矩阵MCtr+：
Figure BDA0002612837980000058
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000061
为母线i未投电容器的组数，QRiL(1≤i≤m)为母线 i单组电抗器的容量，
Figure BDA0002612837980000062
为母线i已投电抗器的组数。
步骤10532：设置i＝1；步骤10533：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单组电抗器容量即QRiL(MVar)，执行潮流计算；步骤10534：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000063
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure BDA0002612837980000064
步骤10535：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000065
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure BDA0002612837980000066
步骤10536：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电容器容量即QR1C(MVar)，执行潮流计算；步骤10537：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000067
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure BDA0002612837980000068
步骤10538：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA00026128379800000611
，计算得到母线BusRi对母线BusC的投电容灵敏度
Figure BDA0002612837980000069
步骤10539：若i＜m，则i＝i+1，返回步骤10533；若i＝m，则基于所述切电抗灵敏度SLi-、所述投电容灵敏度SCi+，形成母线升压灵敏度矩阵SV+：
Figure BDA00026128379800000610
进一步地，所述步骤106基于所述母线降压/升压灵敏度矩阵进行灵敏度组合，计算预期电压调整值，筛选生成预期有效灵敏度组合系数，包括：步骤1061：判断所述被控母线的电压待调整方向：步骤1062：若ΔVC＜0，则采用MCtr-、 SV-对预期电压调整值
Figure BDA0002612837980000071
组合估算：
Figure BDA0002612837980000072
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure BDA0002612837980000073
共有组合数量为：
Figure BDA0002612837980000074
若满足
Figure BDA0002612837980000075
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}为预期有效灵敏度组合系数；
步骤1063：若ΔVC＞0，则采用MCtr+、SV+对预期电压调整值
Figure BDA0002612837980000076
组合估算：
Figure BDA0002612837980000077
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure BDA0002612837980000078
共有组合数量为：
Figure BDA0002612837980000079
若满足
Figure BDA00026128379800000710
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}为预期有效灵敏度组合系数。
进一步地，步骤107基于所述预期有效灵敏度组合系数进行容抗投退后的潮流计算以及校核，得到并输出校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制，包括：步骤1071：基于所述预期有效灵敏度组合系数，在所述基础运行方式潮流数据文件中投退与所述预期有效灵敏度组合系数相对应数量的电容器和电抗器；步骤1072：执行潮流计算；步骤1073：若潮流不收敛，则所述预期有效灵敏度组合系数为校验无效容抗组合系数；若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000081
提取监视母线组当前电压值
Figure BDA0002612837980000082
步骤1074：判断被控母线当前电压值
Figure BDA0002612837980000083
是否满足
Figure BDA0002612837980000084
步骤1075：若不满足，则所述预期有效灵敏度组合系数为校验无效容抗组合系数；若满足，则判断监视母线组当前电压值
Figure BDA0002612837980000085
是否满足
Figure BDA0002612837980000086
步骤1076：若不满足，则所述预期有效灵敏度组合系数为校验无效灵敏度组合系数；若满足，则所述预期有效灵敏度组合系数为校验有效灵敏度组合系数，并输出所述校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制。
第二方面，本发明实施例还提供了一种电网运行方式电压智能调整装置，包括：数据获取单元，用于获取基础运行方式潮流数据文件、电网容抗资源配置文件；矩阵形成和电压限值设定单元，用于根据所述电网容抗资源配置文件形成可控资源矩阵MCtr，同时从所述基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定所述被控母线和所述监视母线组的电压运行限值；数据提取单元，用于基于所述基础运行方式潮流数据文件进行潮流计算并得到潮流计算结果，从所述潮流计算结果中提取所述监视母线组的当前电压值和所述被控母线的当前电压值；电压校核及调整单元，用于对所述监视母线组的当前电压值进行限值校核，同时，基于所述被控母线的当前电压值计算所述被控母线的电压待调整值；灵敏度计算单元，用于根据所述被控母线的电压待调整值和所述可控资源矩阵MCtr，分别进行切电容/投电抗、投电容/切电抗灵敏度计算，形成母线降压/升压灵敏度矩阵；灵敏度组合单元，用于基于所述母线降压/升压灵敏度矩阵进行灵敏度组合，计算预期电压调整值，筛选生成预期有效灵敏度组合系数；有效灵敏度组合系数校核单元，用于基于所述预期有效灵敏度组合系数进行容抗投退后的潮流计算以及校核，得到并输出校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制。
进一步地，所述基础运行方式潮流数据文件包括：发电机母线模型、交直流母线模型、交直流线路模型、变压器模型、电容器模型、电抗器模型、除电容器电抗器以外的其它FACTS设备模型。
进一步地，所述电网容抗资源配置文件包括：具有容抗资源的母线名称、单组电容器容量、已投电容器组数、未投电容器组数、单组电抗器容量、已投电抗器组数、未投电抗器组数。
进一步地，所述根据所述电网容抗资源配置文件形成可控资源矩阵MCtr，包括：根据所述电网容抗资源配置文件形成如下可控资源矩阵MCtr：
Figure BDA0002612837980000091
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000092
为母线i已投电容器的组数，
Figure BDA0002612837980000093
为母线 i未投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure BDA0002612837980000094
为母线i已投电抗器的组数，
Figure BDA0002612837980000095
为母线i未投电抗器的组数。
进一步地，所述从基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定所述被控母线和所述监视母线组的电压运行限值，包括：步骤 1021：从所述基础运行方式潮流数据文件中选定1条被控母线BusC，设定所述被控母线的电压控制目标值为
Figure BDA0002612837980000096
电压控制目标值波动带宽为
Figure BDA0002612837980000097
步骤1022：从所述基础运行方式潮流数据文件中选定t条监视母线BusMi(1≤i≤t)，设定每条监视母线的电压上限值为
Figure BDA0002612837980000098
电压下限值为
Figure BDA0002612837980000099
电压监视限值调整裕度为ΔVM，同时所述t条监视母线的电压上下限值构成如下矩阵Mt：
Figure BDA0002612837980000101
进一步地，所述电压校核及调整单元用于对所述监视母线组的当前电压值进行限值校核，包括：若
Figure BDA0002612837980000102
则更新每条监视母线BusMi的电压上限值为：
Figure BDA0002612837980000103
若
Figure BDA0002612837980000104
则更新每条监视母线BusMi的电压下限值为：
Figure BDA0002612837980000105
其中，
Figure BDA0002612837980000106
为每条监视母线组BusM当前电压值，
Figure BDA0002612837980000107
为每条监视母线的电压上限值，
Figure BDA0002612837980000108
为每条监视母线的电压下限值，ΔVM为电压监视限值调整裕度。
进一步地，所述电压校核及调整单元用于基于所述被控母线的当前电压值计算所述被控母线的电压待调整值，包括：所述被控母线BusC的电压待调整值
Figure BDA0002612837980000109
其中，
Figure BDA00026128379800001010
为所述被控母线的当前电压值，
Figure BDA00026128379800001011
为所述被控母线的电压控制目标值。
进一步地，灵敏度计算单元，用于：步骤1051：判断所述被控母线的电压待调整方向：步骤1052：若ΔVC＜0，根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵；步骤1053：若ΔVC＞0，根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵。
进一步地，所述根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵，包括：步骤10521：基于所述可控资源矩阵MCtr抽取信息得到母线降压资源矩阵MCtr-：
Figure BDA00026128379800001012
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000111
为母线i已投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure BDA0002612837980000112
为母线i未投电抗器的组数；
步骤10522：设置i＝1；步骤10523：根据MCtr-中可控资源，在所述基础运行方式潮流数据文件的母线BusRi退出单组电容器容量即QRiC(MVar)，执行潮流计算；步骤10524：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000113
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure BDA0002612837980000114
其中，
Figure BDA0002612837980000115
为被控母线在步骤104中的电压值， QRiC(1≤i≤m)为被控母线i单组电容器的容量；步骤10525：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000116
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure BDA0002612837980000117
步骤10526：根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电抗器容量即QR1L(MVar)，执行潮流计算；步骤10527：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000118
计算得到母线BusRi对母线 BusC的投电抗灵敏度
Figure BDA0002612837980000119
步骤10528：若潮流不收敛，则根据 MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA00026128379800001110
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure BDA00026128379800001111
步骤 10529：若i＜m，则i＝i+1，返回步骤10523；若i＝m，基于所述切电容灵敏度SCi-、所述投电抗灵敏度SLi+，形成母线降压灵敏度矩阵SV-：
Figure BDA0002612837980000121
进一步地，所述根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵，包括：步骤10531：基于所述可控资源矩阵MCtr抽取信息得到母线升压资源矩阵MCtr+：
Figure BDA0002612837980000122
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000123
为母线i未投电容器的组数，QRiL(1≤i≤m)为母线 i单组电抗器的容量，
Figure BDA0002612837980000124
为母线i已投电抗器的组数。
步骤10532：设置i＝1；步骤10533：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单组电抗器容量即QRiL(MVar)，执行潮流计算；步骤10534：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000125
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure BDA0002612837980000126
步骤10535：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000127
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure BDA0002612837980000128
步骤10536：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电容器容量即QR1C(MVar)，执行潮流计算；步骤10537：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000129
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure BDA00026128379800001210
步骤10538：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000131
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure BDA0002612837980000132
步骤10539：若i＜m，则i＝i+1，返回步骤10533；若i＝m，则基于所述切电抗灵敏度SLi-、所述投电容灵敏度SCi+，形成母线升压灵敏度矩阵SV+：
Figure BDA0002612837980000133
进一步地，所述灵敏度组合单元，用于：步骤1061：判断所述被控母线的电压待调整方向：步骤1062：若ΔVC＜0，则采用MCtr-、SV-对预期电压调整值
Figure BDA0002612837980000134
组合估算：
Figure BDA0002612837980000135
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure BDA0002612837980000136
共有组合数量为：
Figure BDA0002612837980000137
若满足
Figure BDA0002612837980000138
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}为预期有效灵敏度组合系数；
步骤1063：若ΔVC＞0，则采用MCtr+、SV+对预期电压调整值
Figure BDA0002612837980000139
组合估算：
Figure BDA00026128379800001310
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure BDA0002612837980000141
共有组合数量为：
Figure BDA0002612837980000142
若满足
Figure BDA0002612837980000143
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}为预期有效灵敏度组合系数。
进一步地，所述有效灵敏度组合系数校核单元，用于：步骤1071：基于所述预期有效灵敏度组合系数，在所述基础运行方式潮流数据文件中投退与所述预期有效灵敏度组合系数相对应数量的电容器和电抗器；步骤1072：执行潮流计算；步骤1073：若潮流不收敛，则所述预期有效灵敏度组合系数为校验无效容抗组合系数；若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000144
提取监视母线组当前电压值
Figure BDA0002612837980000145
步骤1074：判断被控母线当前电压值
Figure BDA0002612837980000146
是否满足
Figure BDA0002612837980000147
步骤1075：若不满足，则所述预期有效灵敏度组合系数为校验无效容抗组合系数；若满足，则判断监视母线组当前电压值
Figure BDA0002612837980000148
是否满足
Figure BDA0002612837980000149
步骤1076：若不满足，则所述预期有效灵敏度组合系数为校验无效灵敏度组合系数；若满足，则所述预期有效灵敏度组合系数为校验有效灵敏度组合系数，并输出所述校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制。
本发明实施例提供的电压智能调整方法及装置，首先根据容抗资源与被控母线计算容抗投切灵敏度；其次通过灵敏度组合获得被控母线与网络中电容器、电抗器无功资源的支撑关系，动态生成满足目标母线电压安全调整的预期有效灵敏度组合系数；最后通过校验获得满足目标电压调整的有效灵敏度组合系数用于电压调整控制，解决了传统电压稳定仿真分析和无功资源配置严重依赖人工的问题，大大降低了人力成本，提高了计算效率，同时计算分析的精确度也得到了有效保证；可以为方式分析编制人员提供强有力的专业化支撑工具，提高电压控制调整的自动化智能化水平。
附图说明
图1示出了根据本发明实施例的电网运行方式电压智能调整方法的示例性流程图；
图2示出了根据本发明实施例的电网运行方式电压智能调整装置的结构示意图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式，然而，本发明可以用许多不同的形式来实施，并且不局限于此处描述的实施例，提供这些实施例是为了详尽地且完全地公开本发明，并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中，相同的单元/元件使用相同的附图标记。
除非另有说明，此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外，可以理解的是，以通常使用的词典限定的术语，应当被理解为与其相关领域的语境具有一致的含义，而不应该被理解为理想化的或过于正式的意义。
图1示出了根据本发明实施例的电网运行方式电压智能调整方法的示例性流程图。
如图1所示，该方法包括：
步骤101：获取基础运行方式潮流数据文件、电网容抗资源配置文件；
步骤102：根据电网容抗资源配置文件形成可控资源矩阵MCtr，同时从基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定被控母线和监视母线组的电压运行限值；
步骤103：基于基础运行方式潮流数据文件进行潮流计算并得到潮流计算结果，从潮流计算结果中提取监视母线组的当前电压值和被控母线的当前电压值；
步骤104：对监视母线组的当前电压值进行限值校核，同时，基于被控母线的当前电压值计算被控母线的电压待调整值；
步骤105：根据被控母线的电压待调整值和可控资源矩阵MCtr，分别进行切电容/投电抗、投电容/切电抗灵敏度计算，形成母线降压/升压灵敏度矩阵；
步骤106：基于母线降压/升压灵敏度矩阵进行灵敏度组合，计算预期电压调整值，筛选生成预期有效灵敏度组合系数；
步骤107：基于预期有效灵敏度组合系数进行容抗投退后的潮流计算以及校核，得到并输出校验有效灵敏度组合系数，校验有效灵敏度组合系数用于电压调整控制。
本发明实施例中，电网容抗资源配置文件可以为按照预先配置规则得到。被控母线可以为1个母线，监视母线组可以为t个母线(t为大于0的正整数)。步骤102中，根据电网容抗资源配置文件形成可控资源矩阵MCtr，可以在选定被控母线和监视母线组之前，也可以在之后，也可是同时进行，也就说，该步骤只要在步骤105之前即可；步骤102中选定被控母线和监视母线组，可以为根据预先设置的规则进行选定。校验有效灵敏度组合系数可以作为投切容抗的依据，也就是说，可以根据该系数确定需要投运的电容或电抗的数量，进而实现电压调整控制。
上述实施例中，首先根据容抗资源与被控母线计算容抗投切灵敏度；其次通过灵敏度组合获得被控母线与网络中电容器、电抗器无功资源的支撑关系，动态生成满足目标母线电压安全调整的预期有效灵敏度组合系数；最后通过校验获得满足目标电压调整的有效灵敏度组合系数用于电压调整控制，解决了传统电压稳定仿真分析和无功资源配置严重依赖人工的问题，大大降低了人力成本，提高了计算效率，同时计算分析的精确度也得到了有效保证；可以为方式分析编制人员提供强有力的专业化支撑工具，提高电压控制调整的自动化智能化水平。
进一步地，基础运行方式潮流数据文件包括：发电机母线模型、交直流母线模型、交直流线路模型、变压器模型、电容器模型、电抗器模型、除电容器电抗器以外的其它FACTS设备模型。
进一步地，电网容抗资源配置文件包括：具有容抗资源的母线名称、单组电容器容量、已投电容器组数、未投电容器组数、单组电抗器容量、已投电抗器组数、未投电抗器组数。
进一步地，步骤102中根据所述电网容抗资源配置文件形成可控资源矩阵 MCtr，包括：
根据电网容抗资源配置文件形成如下可控资源矩阵MCtr：
Figure BDA0002612837980000171
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000172
为母线i已投电容器的组数，
Figure BDA0002612837980000173
为母线 i未投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure BDA0002612837980000174
为母线i已投电抗器的组数，
Figure BDA0002612837980000175
为母线i未投电抗器的组数。
进一步地，步骤102中从基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定被控母线和监视母线组的电压运行限值，包括：
步骤1021：从基础运行方式潮流数据文件中选定1条被控母线BusC，设定被控母线的电压控制目标值为
Figure BDA0002612837980000176
电压控制目标值波动带宽为
Figure BDA0002612837980000177
步骤1022：从基础运行方式潮流数据文件中选定t条监视母线BusMi(1≤i≤t)，设定每条监视母线的电压上限值为
Figure BDA0002612837980000178
电压下限值为
Figure BDA0002612837980000179
电压监视限值调整裕度为ΔVM，同时t条监视母线的电压上下限值构成如下矩阵 Mt：
Figure BDA0002612837980000181
本发明实施例中，步骤1021和步骤1022可以如该顺序，也可以前后互换，也可以同时进行。
进一步地，步骤104中对监视母线组的当前电压值进行限值校核，包括：
若
Figure BDA0002612837980000182
则更新每条监视母线BusMi的电压上限值为：
Figure BDA0002612837980000183
若
Figure BDA0002612837980000184
则更新每条监视母线BusMi的电压下限值为：
Figure BDA0002612837980000185
其中，
Figure BDA0002612837980000186
为每条监视母线组BusM当前电压值，
Figure BDA0002612837980000187
为每条监视母线的电压上限值，
Figure BDA0002612837980000188
为每条监视母线的电压下限值，ΔVM为电压监视限值调整裕度。
进一步地，步骤104中基于被控母线的当前电压值计算被控母线的电压待调整值，包括：
被控母线BusC的电压待调整值
Figure BDA0002612837980000189
其中，
Figure BDA00026128379800001810
为被控母线的当前电压值，
Figure BDA00026128379800001811
为被控母线的电压控制目标值。
进一步地，步骤105，包括：
步骤1051：判断所述被控母线的电压待调整方向。
步骤1052：若ΔVC＜0，根据可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵。
进一步地，若ΔVC＜0，步骤1052，包括：
步骤10521：基于可控资源矩阵MCtr抽取信息得到母线降压资源矩阵MCtr-：
Figure BDA0002612837980000191
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000192
为母线i已投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure BDA0002612837980000193
为母线i未投电抗器的组数；
步骤10522：设置i＝1；
步骤10523：根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线 BusRi退出单组电容器容量即QRiC(MVar)，执行潮流计算；
步骤10524：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000194
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure BDA0002612837980000195
其中，
Figure BDA0002612837980000196
为被控母线在步骤104中的电压值，QRiC(1≤i≤m)为被控母线i单组电容器的容量；
步骤10525：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000197
计算得到母线BusRi对母线 BusC的切电容灵敏度
Figure BDA0002612837980000198
步骤10526：根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线 BusRi投入单组电抗器容量即QR1L(MVar)，执行潮流计算；
步骤10527：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000199
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure BDA00026128379800001910
步骤10528：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA00026128379800001911
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure BDA0002612837980000201
步骤10529：若i＜m，则i＝i+1，返回步骤10523；若i＝m，基于切电容灵敏度SCi-、投电抗灵敏度SLi+，形成母线降压灵敏度矩阵SV-：
Figure BDA0002612837980000202
本发明实施例中，从i＝1开始，循环步骤10523-10529，对矩阵Mctr-中每一行的元素进行计算。
步骤1053：若ΔVC＞0，根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵。
进一步地，若ΔVC＞0，步骤1053，包括：
步骤10531：基于所述可控资源矩阵MCtr抽取信息得到母线升压资源矩阵 MCtr+：
Figure BDA0002612837980000203
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000204
为母线i未投电容器的组数，QRiL(1≤i≤m)为母线 i单组电抗器的容量，
Figure BDA0002612837980000205
为母线i已投电抗器的组数。
步骤10532：设置i＝1；
步骤10533：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线 BusRi退出单组电抗器容量即QRiL(MVar)，执行潮流计算；
步骤10534：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000211
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure BDA0002612837980000212
步骤10535：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000219
计算得到母线BusRi对母线 BusC的切电抗灵敏度
Figure BDA0002612837980000213
步骤10536：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线 BusRi投入单组电容器容量即QR1C(MVar)，执行潮流计算；
步骤10537：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000214
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure BDA0002612837980000215
步骤10538：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000216
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure BDA0002612837980000217
步骤10539：若i＜m，则i＝i+1，返回步骤10533；若i＝m，则基于切电抗灵敏度SLi-、投电容灵敏度SCi+，形成母线升压灵敏度矩阵SV+：
Figure BDA0002612837980000218
本发明实施例中，从i＝1开始，循环步骤10533-10539，对矩阵Mctr+中每一行的元素进行计算。
进一步地，步骤106，包括：
步骤1061：判断被控母线的电压待调整方向：
步骤1062：若ΔVC＜0，则采用MCtr-、SV-对预期电压调整值
Figure BDA0002612837980000221
组合估算：
Figure BDA0002612837980000222
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure BDA0002612837980000223
共有组合数量为：
Figure BDA0002612837980000224
若满足
Figure BDA0002612837980000225
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}
为预期有效灵敏度组合系数；
步骤1063：若ΔVC＞0，则采用MCtr+、SV+对预期电压调整值
Figure BDA0002612837980000226
组合估算：
Figure BDA0002612837980000227
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure BDA0002612837980000228
共有组合数量为：
Figure BDA0002612837980000229
若满足
Figure BDA00026128379800002210
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}
为预期有效灵敏度组合系数。
进一步地，步骤107，包括：
步骤1071：基于预期有效灵敏度组合系数，在基础运行方式潮流数据文件中投退与预期有效灵敏度组合系数相对应数量的电容器和电抗器；
步骤1072：执行潮流计算；
步骤1073：若潮流不收敛，则预期有效灵敏度组合系数为校验无效容抗组合系数；若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000231
提取监视母线组当前电压值
Figure BDA0002612837980000232
步骤1074：判断被控母线当前电压值
Figure BDA0002612837980000233
是否满足
Figure BDA0002612837980000234
步骤1075：若不满足，则预期有效灵敏度组合系数为校验无效容抗组合系数；若满足，则判断监视母线组当前电压值
Figure BDA0002612837980000235
是否满足
Figure BDA0002612837980000236
步骤1076：若不满足，则预期有效灵敏度组合系数为校验无效灵敏度组合系数；若满足，则预期有效灵敏度组合系数为校验有效灵敏度组合系数，并输出校验有效灵敏度组合系数，校验有效灵敏度组合系数用于电压调整控制。
图2示出了根据本发明实施例的电网运行方式电压智能调整装置的结构示意图。
如图2所示，该装置包括：
数据获取单元201，用于获取基础运行方式潮流数据文件、电网容抗资源配置文件；
矩阵形成和电压限值设定单元202，用于根据电网容抗资源配置文件形成可控资源矩阵MCtr，同时从基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定被控母线和监视母线组的电压运行限值；
数据提取单元203，用于基于基础运行方式潮流数据文件进行潮流计算并得到潮流计算结果，从潮流计算结果中提取监视母线组的当前电压值和被控母线的当前电压值；
电压校核及调整单元204，用于对监视母线组的当前电压值进行限值校核，同时，基于被控母线的当前电压值计算被控母线的电压待调整值；
灵敏度计算单元205，用于根据被控母线的电压待调整值和可控资源矩阵 MCtr，分别进行切电容/投电抗、投电容/切电抗灵敏度计算，形成母线降压/升压灵敏度矩阵；
灵敏度组合单元206，用于基于母线降压/升压灵敏度矩阵进行灵敏度组合，计算预期电压调整值，筛选生成预期有效灵敏度组合系数；
有效灵敏度组合系数校核单元207，用于基于预期有效灵敏度组合系数进行容抗投退后的潮流计算以及校核，得到并输出校验有效灵敏度组合系数，校验有效灵敏度组合系数用于电压调整控制。
本发明实施例中，电网容抗资源配置文件可以为按照预先配置规则得到。被控母线可以为1个母线，监视母线组可以为t个母线(t为大于0的正整数)。矩阵形成和电压限值设定单元202，用于根据电网容抗资源配置文件形成可控资源矩阵MCtr，可以在选定被控母线和监视母线组之前，也可以在之后，也可是同时进行，也就说，该步骤只要在灵敏度计算单元205根据被控母线的电压待调整值和可控资源矩阵MCtr分别进行切电容/投电抗、投电容/切电抗灵敏度计算之前即可；矩阵形成和电压限值设定单元202，用于选定被控母线和监视母线组，可以为根据预先设置的规则进行选定。校验有效灵敏度组合系数可以作为投切容抗的依据，也就是说，可以根据该系数确定需要投运的电容或电抗的数量，进而实现电压调整控制。
上述实施例中，首先根据容抗资源与被控母线计算容抗投切灵敏度；其次通过灵敏度组合获得被控母线与网络中电容器、电抗器无功资源的支撑关系，动态生成满足目标母线电压安全调整的预期有效灵敏度组合系数；最后通过校验获得满足目标电压调整的有效灵敏度组合系数用于电压调整控制，解决了传统电压稳定仿真分析和无功资源配置严重依赖人工的问题，大大降低了人力成本，提高了计算效率，同时计算分析的精确度也得到了有效保证；可以为方式分析编制人员提供强有力的专业化支撑工具，提高电压控制调整的自动化智能化水平。
进一步地，基础运行方式潮流数据文件包括：发电机母线模型、交直流母线模型、交直流线路模型、变压器模型、电容器模型、电抗器模型、除电容器电抗器以外的其它FACTS设备模型。
进一步地，电网容抗资源配置文件包括：具有容抗资源的母线名称、单组电容器容量、已投电容器组数、未投电容器组数、单组电抗器容量、已投电抗器组数、未投电抗器组数。
进一步地，根据电网容抗资源配置文件形成可控资源矩阵MCtr，包括：
根据电网容抗资源配置文件形成如下可控资源矩阵MCtr：
Figure BDA0002612837980000251
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000252
为母线i已投电容器的组数，
Figure BDA0002612837980000253
为母线 i未投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure BDA0002612837980000254
为母线i已投电抗器的组数，
Figure BDA0002612837980000255
为母线i未投电抗器的组数。
进一步地，从基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定被控母线和监视母线组的电压运行限值，包括：
步骤1021：从基础运行方式潮流数据文件中选定1条被控母线BusC，设定被控母线的电压控制目标值为
Figure BDA0002612837980000256
电压控制目标值波动带宽为
Figure BDA0002612837980000257
步骤1022：从基础运行方式潮流数据文件中选定t条监视母线BusMi(1≤i≤t)，设定每条监视母线的电压上限值为
Figure BDA0002612837980000258
电压下限值为
Figure BDA0002612837980000259
电压监视限值调整裕度为ΔVM，同时t条监视母线的电压上下限值构成如下矩阵 Mt：
Figure BDA00026128379800002510
进一步地，电压校核及调整单元204用于对监视母线组的当前电压值进行限值校核，包括：
若
Figure 1
则更新每条监视母线BusMi的电压上限值为：
Figure BDA0002612837980000262
若
Figure BDA0002612837980000263
则更新每条监视母线BusMi的电压下限值为：
Figure BDA0002612837980000264
其中，
Figure BDA0002612837980000265
为每条监视母线组BusM当前电压值，
Figure BDA0002612837980000266
为每条监视母线的电压上限值，
Figure BDA0002612837980000267
为每条监视母线的电压下限值，ΔVM为电压监视限值调整裕度。
进一步地，电压校核及调整单元204用于基于被控母线的当前电压值计算被控母线的电压待调整值，包括：
被控母线BusC的电压待调整值
Figure BDA0002612837980000268
其中，
Figure BDA0002612837980000269
为被控母线的当前电压值，
Figure BDA00026128379800002610
为被控母线的电压控制目标值。
进一步地，灵敏度计算单元205，用于：
步骤1051：判断所述被控母线的电压待调整方向。
步骤1052：若ΔVC＜0，根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵。
进一步地，若ΔVC＜0，步骤1052，包括：
步骤10521：基于可控资源矩阵MCtr抽取信息得到母线降压资源矩阵MCtr-：
Figure BDA00026128379800002611
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA00026128379800002612
为母线i已投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure BDA00026128379800002613
为母线i未投电抗器的组数；
步骤10522：设置i＝1；
步骤10523：根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线 BusRi退出单组电容器容量即QRiC(MVar)，执行潮流计算；
步骤10524：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000271
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure BDA0002612837980000272
其中，
Figure BDA0002612837980000273
为被控母线在步骤104中的电压值，QRiC(1≤i≤m)为被控母线i单组电容器的容量；
步骤10525：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000274
计算得到母线BusRi对母线 BusC的切电容灵敏度
Figure BDA0002612837980000275
步骤10526：根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线 BusRi投入单组电抗器容量即QR1L(MVar)，执行潮流计算；
步骤10527：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000276
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure BDA0002612837980000277
步骤10528：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000278
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure BDA0002612837980000279
步骤10529：若i＜m，则i＝i+1，返回步骤10523；若i＝m，基于切电容灵敏度SCi-、投电抗灵敏度SLi+，形成母线降压灵敏度矩阵SV-：
Figure BDA00026128379800002710
本发明实施例中，从i＝1开始，循环步骤10523-10529，对矩阵Mctr-中每一行的元素进行计算。
步骤1053：若ΔVC＞0，根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵。
进一步地，若ΔVC＞0，步骤1053，包括：
步骤10531：基于所述可控资源矩阵MCtr抽取信息得到母线升压资源矩阵 MCtr+：
Figure BDA0002612837980000281
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure BDA0002612837980000282
为母线i未投电容器的组数，QRiL(1≤i≤m)为母线 i单组电抗器的容量，
Figure BDA0002612837980000283
为母线i已投电抗器的组数。
步骤10532：设置i＝1；
步骤10533：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线 BusRi退出单组电抗器容量即QRiL(MVar)，执行潮流计算；
步骤10534：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000284
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure BDA0002612837980000285
步骤10535：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000286
计算得到母线BusRi对母线 BusC的切电抗灵敏度
Figure BDA0002612837980000287
步骤10536：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电容器容量即QR1C(MVar)，执行潮流计算；
步骤10537：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000291
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure BDA0002612837980000292
步骤10538：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000293
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure BDA0002612837980000294
步骤10539：若i＜m，则i＝i+1，返回步骤10533；若i＝m，则基于切电抗灵敏度SLi-、投电容灵敏度SCi+，形成母线升压灵敏度矩阵SV+：
Figure BDA0002612837980000295
本发明实施例中，从i＝1开始，循环步骤10533-10539，对矩阵Mctr+中每一行的元素进行计算。
进一步地，灵敏度组合单元206，用于：
步骤1061：判断被控母线的电压待调整方向：
步骤1062：若ΔVC＜0，则采用MCtr-、SV-对预期电压调整值
Figure BDA0002612837980000296
组合估算：
Figure BDA0002612837980000297
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure BDA0002612837980000298
共有组合数量为：
Figure BDA0002612837980000299
若满足
Figure BDA0002612837980000301
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}
为预期有效灵敏度组合系数；
步骤1063：若ΔVC＞0，则采用MCtr+、SV+对预期电压调整值
Figure BDA0002612837980000302
组合估算：
Figure BDA0002612837980000303
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure BDA0002612837980000304
共有组合数量为：
Figure BDA0002612837980000305
若满足
Figure BDA0002612837980000306
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}
为预期有效灵敏度组合系数。
进一步地，有效灵敏度组合系数校核单元207，用于：
步骤1071：基于预期有效灵敏度组合系数，在基础运行方式潮流数据文件中投退与预期有效灵敏度组合系数相对应数量的电容器和电抗器；
步骤1072：执行潮流计算；
步骤1073：若潮流不收敛，则预期有效灵敏度组合系数为校验无效容抗组合系数；若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure BDA0002612837980000307
提取监视母线组当前电压值
Figure BDA0002612837980000308
步骤1074：判断被控母线当前电压值
Figure BDA0002612837980000309
是否满足
Figure BDA00026128379800003010
步骤1075：若不满足，则预期有效灵敏度组合系数为校验无效容抗组合系数；若满足，则判断监视母线组当前电压值
Figure BDA00026128379800003011
是否满足
Figure BDA00026128379800003012
步骤1076：若不满足，则预期有效灵敏度组合系数为校验无效灵敏度组合系数；若满足，则预期有效灵敏度组合系数为校验有效灵敏度组合系数，并输出校验有效灵敏度组合系数，校验有效灵敏度组合系数用于电压调整控制。
已经通过参考少量实施方式描述了本发明。然而，本领域技术人员所公知的，正如附带的专利权利要求所限定的，除了本发明以上公开的其他的实施例等同地落在本发明的范围内。
通常地，在权利要求中使用的所有术语都根据他们在技术领域的通常含义被解释，除非在其中被另外明确地定义。所有的参考“一个/所述/该[装置、组件等]”都被开放地解释为所述装置、组件等中的至少一个实例，除非另外明确地说明。这里公开的任何方法的步骤都没必要以公开的准确的顺序运行，除非明确地说明。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、 CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和 /或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/ 或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

Claims (24)

1.一种电网运行方式电压智能调整方法，其特征在于，所述方法包括：
步骤101：获取基础运行方式潮流数据文件、电网容抗资源配置文件；
步骤102：根据所述电网容抗资源配置文件形成可控资源矩阵MCtr，同时从所述基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定所述被控母线和所述监视母线组的电压运行限值；
步骤103：基于所述基础运行方式潮流数据文件进行潮流计算并得到潮流计算结果，从所述潮流计算结果中提取所述监视母线组的当前电压值和所述被控母线的当前电压值；
步骤104：对所述监视母线组的当前电压值进行限值校核，同时，基于所述被控母线的当前电压值计算所述被控母线的电压待调整值；
步骤105：根据所述被控母线的电压待调整值和所述可控资源矩阵MCtr，分别进行切电容/投电抗、投电容/切电抗灵敏度计算，形成母线降压/升压灵敏度矩阵；
步骤106：基于所述母线降压/升压灵敏度矩阵进行灵敏度组合，计算预期电压调整值，筛选生成预期有效灵敏度组合系数；
步骤107：基于所述预期有效灵敏度组合系数进行容抗投退后的潮流计算以及校核，得到并输出校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制。
2.根据权利要求1所述的方法，其特征在于，所述基础运行方式潮流数据文件包括：发电机母线模型、交直流母线模型、交直流线路模型、变压器模型、电容器模型、电抗器模型、除电容器电抗器以外的其它FACTS设备模型。
3.根据权利要求1所述的方法，其特征在于，所述电网容抗资源配置文件包括：具有容抗资源的母线名称、单组电容器容量、已投电容器组数、未投电容器组数、单组电抗器容量、已投电抗器组数、未投电抗器组数。
4.根据权利要求3所述的方法，其特征在于，所述步骤102中根据所述电网容抗资源配置文件形成可控资源矩阵MCtr，包括：
根据所述电网容抗资源配置文件形成如下可控资源矩阵MCtr：
Figure FDA0002612837970000021
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure FDA0002612837970000022
为母线i已投电容器的组数，
Figure FDA0002612837970000023
为母线i未投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure FDA0002612837970000024
为母线i已投电抗器的组数，
Figure FDA0002612837970000025
为母线i未投电抗器的组数。
5.根据权利要求1所述的方法，其特征在于，所述步骤102中从所述基础运行方式潮流数据文件中选定母线和监视母线组，并分别设定所述被控母线和所述监视母线组的电压运行限值，包括：
步骤1021：从所述基础运行方式潮流数据文件中选定1条被控母线BusC，设定所述被控母线的电压控制目标值为
Figure FDA0002612837970000026
电压控制目标值波动带宽为
Figure FDA0002612837970000027
步骤1022：从所述基础运行方式潮流数据文件中选定t条监视母线BusMi(1≤i≤t)，设定每条监视母线的电压上限值为
Figure FDA0002612837970000028
电压下限值为
Figure FDA0002612837970000029
电压监视限值调整裕度为ΔVM，同时所述t条监视母线的电压上下限值构成如下矩阵Mt：
Figure FDA0002612837970000031
6.根据权利要求1所述的方法，其特征在于，所述步骤104中对所述监视母线组的当前电压值进行限值校核，包括：
若
Figure FDA0002612837970000032
则更新每条监视母线BusMi的电压上限值为：
Figure FDA0002612837970000033
若
Figure FDA0002612837970000034
则更新每条监视母线BusMi的电压下限值为：
Figure FDA0002612837970000035
其中，
Figure FDA0002612837970000036
为每条监视母线组BusM当前电压值，
Figure FDA0002612837970000037
为每条监视母线的电压上限值，
Figure FDA0002612837970000038
为每条监视母线的电压下限值，ΔVM为电压监视限值调整裕度。
7.根据权利要求1所述的方法，其特征在于，所述步骤104中基于所述被控母线的当前电压值计算所述被控母线的电压待调整值，包括：
所述被控母线BusC的电压待调整值
Figure FDA0002612837970000039
其中，
Figure FDA00026128379700000310
为所述被控母线的当前电压值，
Figure FDA00026128379700000311
为所述被控母线的电压控制目标值。
8.根据权利要求4所述的方法，其特征在于，所述步骤105根据所述被控母线的电压待调整值和所述可控资源矩阵MCtr，分别进行切电容/投电抗、投电容/切电抗灵敏度计算，形成母线降压/升压灵敏度矩阵，包括：
步骤1051：判断所述被控母线的电压待调整方向：
步骤1052：若ΔVC＜0，根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵；
步骤1053：若ΔVC＞0，根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵。
9.根据权利要求8所述的方法，其特征在于，所述根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵，包括：
步骤10521：基于所述可控资源矩阵MCtr抽取信息得到母线降压资源矩阵MCtr-：
Figure FDA0002612837970000041
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure FDA0002612837970000042
为母线i已投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure FDA0002612837970000043
为母线i未投电抗器的组数；
步骤10522：设置i＝1；
步骤10523：根据MCtr-中可控资源，在所述基础运行方式潮流数据文件的母线BusRi退出单组电容器容量即QRiC(MVar)，执行潮流计算；
步骤10524：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000044
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure FDA0002612837970000045
其中，
Figure FDA0002612837970000046
为被控母线在步骤104中的电压值，QRiC(1≤i≤m)为被控母线i单组电容器的容量；
步骤10525：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000047
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure FDA0002612837970000048
步骤10526：根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电抗器容量即QR1L(MVar)，执行潮流计算；
步骤10527：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000051
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure FDA0002612837970000052
步骤10528：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000053
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure FDA0002612837970000054
步骤10529：若i＜m，则i＝i+1，返回步骤10523；若i＝m，基于所述切电容灵敏度SCi-、所述投电抗灵敏度SLi+，形成母线降压灵敏度矩阵SV-：
Figure FDA0002612837970000055
10.根据权利要求8所述的方法，其特征在于，所述根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵，包括：
步骤10531：基于所述可控资源矩阵MCtr抽取信息得到母线升压资源矩阵MCtr+：
Figure FDA0002612837970000056
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure FDA0002612837970000057
为母线i未投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure FDA0002612837970000061
为母线i已投电抗器的组数。
步骤10532：设置i＝1；
步骤10533：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单组电抗器容量即QRiL(MVar)，执行潮流计算；
步骤10534：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000062
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure FDA0002612837970000063
步骤10535：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000064
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure FDA0002612837970000065
步骤10536：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电容器容量即QR1C(MVar)，执行潮流计算；
步骤10537：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000066
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure FDA0002612837970000067
步骤10538：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000068
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure FDA0002612837970000069
步骤10539：若i＜m，则i＝i+1，返回步骤10533；若i＝m，则基于所述切电抗灵敏度SLi-、所述投电容灵敏度SCi+，形成母线升压灵敏度矩阵SV+：
Figure FDA0002612837970000071
11.根据权利要求9或10所述的方法，其特征在于，所述步骤106基于所述母线降压/升压灵敏度矩阵进行灵敏度组合，计算预期电压调整值，筛选生成预期有效灵敏度组合系数，包括：
步骤1061：判断所述被控母线的电压待调整方向：
步骤1062：若ΔVC＜0，则采用MCtr-、SV-对预期电压调整值
Figure FDA0002612837970000072
组合估算：
Figure FDA0002612837970000073
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure FDA0002612837970000074
共有组合数量为：
Figure FDA0002612837970000075
若满足
Figure FDA0002612837970000076
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}
为预期有效灵敏度组合系数；
步骤1063：若ΔVC＞0，则采用MCtr+、SV+对预期电压调整值
Figure FDA0002612837970000077
组合估算：
Figure FDA0002612837970000078
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure FDA0002612837970000079
共有组合数量为：
Figure FDA0002612837970000081
若满足
Figure FDA0002612837970000082
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}
为预期有效灵敏度组合系数。
12.根据权利要求1所述的方法，其特征在于，所述步骤107基于所述预期有效灵敏度组合系数进行容抗投退后的潮流计算以及校核，得到并输出校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制，包括：
步骤1071：基于所述预期有效灵敏度组合系数，在所述基础运行方式潮流数据文件中投退与所述预期有效灵敏度组合系数相对应数量的电容器和电抗器；
步骤1072：执行潮流计算；
步骤1073：若潮流不收敛，则所述预期有效灵敏度组合系数为校验无效容抗组合系数；若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000083
提取监视母线组当前电压值
Figure FDA0002612837970000084
步骤1074：判断被控母线当前电压值
Figure FDA0002612837970000085
是否满足
Figure FDA0002612837970000086
步骤1075：若不满足，则所述预期有效灵敏度组合系数为校验无效容抗组合系数；若满足，则判断监视母线组当前电压值
Figure FDA0002612837970000087
是否满足
Figure FDA0002612837970000088
步骤1076：若不满足，则所述预期有效灵敏度组合系数为校验无效灵敏度组合系数；若满足，则所述预期有效灵敏度组合系数为校验有效灵敏度组合系数，并输出所述校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制。
13.一种电网运行方式电压智能调整装置，其特征在于，所述装置包括：
数据获取单元，用于获取基础运行方式潮流数据文件、电网容抗资源配置文件；
矩阵形成和电压限值设定单元，用于根据所述电网容抗资源配置文件形成可控资源矩阵MCtr，同时从所述基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定所述被控母线和所述监视母线组的电压运行限值；
数据提取单元，用于基于所述基础运行方式潮流数据文件进行潮流计算并得到潮流计算结果，从所述潮流计算结果中提取所述监视母线组的当前电压值和所述被控母线的当前电压值；
电压校核及调整单元，用于对所述监视母线组的当前电压值进行限值校核，同时，基于所述被控母线的当前电压值计算所述被控母线的电压待调整值；
灵敏度计算单元，用于根据所述被控母线的电压待调整值和所述可控资源矩阵MCtr，分别进行切电容/投电抗、投电容/切电抗灵敏度计算，形成母线降压/升压灵敏度矩阵；
灵敏度组合单元，用于基于所述母线降压/升压灵敏度矩阵进行灵敏度组合，计算预期电压调整值，筛选生成预期有效灵敏度组合系数；
有效灵敏度组合系数校核单元，用于基于所述预期有效灵敏度组合系数进行容抗投退后的潮流计算以及校核，得到并输出校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制。
14.根据权利要求13所述的装置，其特征在于，所述基础运行方式潮流数据文件包括：发电机母线模型、交直流母线模型、交直流线路模型、变压器模型、电容器模型、电抗器模型、除电容器电抗器以外的其它FACTS设备模型。
15.根据权利要求13所述的装置，其特征在于，所述电网容抗资源配置文件包括：具有容抗资源的母线名称、单组电容器容量、已投电容器组数、未投电容器组数、单组电抗器容量、已投电抗器组数、未投电抗器组数。
16.根据权利要求15所述的装置，其特征在于，所述根据所述电网容抗资源配置文件形成可控资源矩阵MCtr，包括：
根据所述电网容抗资源配置文件形成如下可控资源矩阵MCtr：
Figure FDA0002612837970000101
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure FDA0002612837970000102
为母线i已投电容器的组数，
Figure FDA0002612837970000103
为母线i未投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure FDA0002612837970000104
为母线i已投电抗器的组数，
Figure FDA0002612837970000105
为母线i未投电抗器的组数。
17.根据权利要求13所述的装置，其特征在于，所述从基础运行方式潮流数据文件中选定被控母线和监视母线组，并分别设定所述被控母线和所述监视母线组的电压运行限值，包括：
步骤1021：从所述基础运行方式潮流数据文件中选定1条被控母线BusC，设定所述被控母线的电压控制目标值为
Figure FDA0002612837970000106
电压控制目标值波动带宽为
Figure FDA0002612837970000107
步骤1022：从所述基础运行方式潮流数据文件中选定t条监视母线BusMi(1≤i≤t)，设定每条监视母线的电压上限值为
Figure FDA0002612837970000108
电压下限值为
Figure FDA0002612837970000109
电压监视限值调整裕度为ΔVM，同时所述t条监视母线的电压上下限值构成如下矩阵Mt：
Figure FDA0002612837970000111
18.根据权利要求13所述的装置，其特征在于，所述电压校核及调整单元用于对所述监视母线组的当前电压值进行限值校核，包括：
若
Figure FDA0002612837970000112
则更新每条监视母线BusMi的电压上限值为：
若
Figure FDA0002612837970000114
则更新每条监视母线BusMi的电压下限值为：
Figure FDA0002612837970000115
其中，
Figure FDA0002612837970000116
为每条监视母线组BusM当前电压值，
Figure FDA0002612837970000117
为每条监视母线的电压上限值，
Figure FDA0002612837970000118
为每条监视母线的电压下限值，ΔVM为电压监视限值调整裕度。
19.根据权利要求13所述的装置，其特征在于，所述电压校核及调整单元用于基于所述被控母线的当前电压值计算所述被控母线的电压待调整值，包括：
所述被控母线BusC的电压待调整值
Figure FDA0002612837970000119
其中，VC base为所述被控母线的当前电压值，
Figure FDA00026128379700001110
为所述被控母线的电压控制目标值。
20.根据权利要求16所述的装置，其特征在于，所述灵敏度计算单元，用于：
步骤1051：判断所述被控母线的电压待调整方向：
步骤1052：若ΔVC＜0，根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵；
步骤1053：若ΔVC＞0，根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵。
21.根据权利要求20所述的装置，其特征在于，所述根据所述可控资源矩阵MCtr，进行切电容/投电抗灵敏度计算，形成母线降压灵敏度矩阵，包括：
步骤10521：基于所述可控资源矩阵MCtr抽取信息得到母线降压资源矩阵MCtr-：
Figure FDA0002612837970000121
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure FDA0002612837970000122
为母线i已投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure FDA0002612837970000123
为母线i未投电抗器的组数；
步骤10522：设置i＝1；
步骤10523：根据MCtr-中可控资源，在所述基础运行方式潮流数据文件的母线BusRi退出单组电容器容量即QRiC(MVar)，执行潮流计算；
步骤10524：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000124
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure FDA0002612837970000125
其中，
Figure FDA0002612837970000126
为被控母线在步骤104中的电压值，QRiC(1≤i≤m)为被控母线i单组电容器的容量；
步骤10525：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000127
计算得到母线BusRi对母线BusC的切电容灵敏度
Figure FDA0002612837970000128
步骤10526：根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电抗器容量即QR1L(MVar)，执行潮流计算；
步骤10527：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000131
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure FDA0002612837970000132
步骤10528：若潮流不收敛，则根据MCtr-中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000133
计算得到母线BusRi对母线BusC的投电抗灵敏度
Figure FDA0002612837970000134
步骤10529：若i＜m，则i＝i+1，返回步骤10523；若i＝m，基于所述切电容灵敏度SCi-、所述投电抗灵敏度SLi+，形成母线降压灵敏度矩阵SV-：
Figure FDA0002612837970000135
22.根据权利要求20所述的装置，其特征在于，所述根据所述可控资源矩阵MCtr，进行投电容/切电抗灵敏度计算，形成母线升压灵敏度矩阵，包括：
步骤10531：基于所述可控资源矩阵MCtr抽取信息得到母线升压资源矩阵MCtr+：
Figure FDA0002612837970000136
其中，BusRi(1≤i≤m)为具有容抗资源的母线名称，QRiC(1≤i≤m)为母线i单组电容器的容量，
Figure FDA0002612837970000137
为母线i未投电容器的组数，QRiL(1≤i≤m)为母线i单组电抗器的容量，
Figure FDA0002612837970000141
为母线i已投电抗器的组数。
步骤10532：设置i＝1；
步骤10533：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单组电抗器容量即QRiL(MVar)，执行潮流计算；
步骤10534：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000142
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure FDA0002612837970000143
步骤10535：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi退出单位电抗器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000144
计算得到母线BusRi对母线BusC的切电抗灵敏度
Figure FDA0002612837970000145
步骤10536：根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单组电容器容量即QR1C(MVar)，执行潮流计算；
步骤10537：若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000146
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure FDA0002612837970000147
步骤10538：若潮流不收敛，则根据MCtr+中可控资源，在基础运行方式潮流数据文件的母线BusRi投入单位电容器容量即1(MVar)，执行潮流计算，基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000148
计算得到母线BusRi对母线BusC的投电容灵敏度
Figure FDA0002612837970000149
步骤10539：若i＜m，则i＝i+1，返回步骤10533；若i＝m，则基于所述切电抗灵敏度SLi-、所述投电容灵敏度SCi+，形成母线升压灵敏度矩阵SV+：
Figure FDA0002612837970000151
23.根据权利要求21或22所述的装置，其特征在于，所述灵敏度组合单元，用于：
步骤1061：判断所述被控母线的电压待调整方向：
步骤1062：若ΔVC＜0，则采用MCtr-、SV-对预期电压调整值
Figure FDA0002612837970000152
组合估算：
Figure FDA0002612837970000153
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure FDA0002612837970000154
共有组合数量为：
Figure FDA0002612837970000155
若满足
Figure FDA0002612837970000156
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}
为预期有效灵敏度组合系数；
步骤1063：若ΔVC＞0，则采用MCtr+、SV+对预期电压调整值
Figure FDA0002612837970000157
组合估算：
Figure FDA0002612837970000158
其中，组合系数NCi(1≤i≤m)、NLi(1≤i≤m)均为整数，取值范围分别为
Figure FDA0002612837970000159
共有组合数量为：
Figure FDA0002612837970000161
若满足
Figure FDA0002612837970000162
则提取组合系数：
{NC1,NC2,…,NCi,…,NCm,NL1,NL2,…,NLi,…,NLm}
为预期有效灵敏度组合系数。
24.根据权利要求13所述的装置，其特征在于，所述有效灵敏度组合系数校核单元，用于：
步骤1071：基于所述预期有效灵敏度组合系数，在所述基础运行方式潮流数据文件中投退与所述预期有效灵敏度组合系数相对应数量的电容器和电抗器；
步骤1072：执行潮流计算；
步骤1073：若潮流不收敛，则所述预期有效灵敏度组合系数为校验无效容抗组合系数；若潮流收敛，则基于潮流计算结果提取被控母线BusC当前电压值
Figure FDA0002612837970000163
提取监视母线组当前电压值
Figure FDA0002612837970000164
步骤1074：判断被控母线当前电压值
Figure FDA0002612837970000165
是否满足
Figure FDA0002612837970000166
步骤1075：若不满足，则所述预期有效灵敏度组合系数为校验无效容抗组合系数；若满足，则判断监视母线组当前电压值
Figure FDA0002612837970000167
是否满足
Figure FDA0002612837970000168
步骤1076：若不满足，则所述预期有效灵敏度组合系数为校验无效灵敏度组合系数；若满足，则所述预期有效灵敏度组合系数为校验有效灵敏度组合系数，并输出所述校验有效灵敏度组合系数，所述校验有效灵敏度组合系数用于电压调整控制。
CN202010760141.9A 2020-07-31 2020-07-31 一种电网运行方式电压智能调整方法及装置 Active CN112003288B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010760141.9A CN112003288B (zh) 2020-07-31 2020-07-31 一种电网运行方式电压智能调整方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010760141.9A CN112003288B (zh) 2020-07-31 2020-07-31 一种电网运行方式电压智能调整方法及装置

Publications (2)

Publication Number Publication Date
CN112003288A true CN112003288A (zh) 2020-11-27
CN112003288B CN112003288B (zh) 2023-08-18

Family

ID=73464163

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010760141.9A Active CN112003288B (zh) 2020-07-31 2020-07-31 一种电网运行方式电压智能调整方法及装置

Country Status (1)

Country Link
CN (1) CN112003288B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113162056A (zh) * 2021-03-22 2021-07-23 华北电力科学研究院有限责任公司 特高压电网自动电压控制方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101034811A (zh) * 2007-04-27 2007-09-12 清华大学 静态混成自动电压控制最高层安全性调控方法
CN102856919A (zh) * 2012-09-21 2013-01-02 河海大学 一种混合经济压差和灵敏度分析的无功优化在线控制方法
CN103490426A (zh) * 2013-09-25 2014-01-01 国家电网公司 大规模风电接入的弱送端电力系统电压稳定综合控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101034811A (zh) * 2007-04-27 2007-09-12 清华大学 静态混成自动电压控制最高层安全性调控方法
CN102856919A (zh) * 2012-09-21 2013-01-02 河海大学 一种混合经济压差和灵敏度分析的无功优化在线控制方法
CN103490426A (zh) * 2013-09-25 2014-01-01 国家电网公司 大规模风电接入的弱送端电力系统电压稳定综合控制方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113162056A (zh) * 2021-03-22 2021-07-23 华北电力科学研究院有限责任公司 特高压电网自动电压控制方法和装置

Also Published As

Publication number Publication date
CN112003288B (zh) 2023-08-18

Similar Documents

Publication Publication Date Title
Das et al. Optimal placement of distributed energy storage systems in distribution networks using artificial bee colony algorithm
CN106451556A (zh) 一种确定配电网馈线接纳分布式电源能力的方法和装置
CN102403716A (zh) 一种多馈入交直流电网动态等值的方法
CN110323979B (zh) 一种考虑电压稳定的发电机励磁调差系数优化整定方法
CN104701858A (zh) 计及分区动态无功储备的无功电压控制方法
CN105468877B (zh) 一种光伏电站无功支撑方法
CN110970905B (zh) 一种优化风电场电压控制能力的多电源无功配合方法及系统
CN109390953A (zh) 含分布式电源和电动汽车的低压配电网无功电压协调控制方法和系统
CN103730900A (zh) 电力系统多时间尺度的省地县一体化无功优化方法
CN112003288A (zh) 一种电网运行方式电压智能调整方法及装置
CN112072692A (zh) 一种新能源发电场站的阻抗等值方法及装置
CN108376997B (zh) 一种考虑分布式电源不确定性的有源配电网孤岛划分方法
CN105279297B (zh) 一种基于工程应用的pss/e电力系统等值方法
CN110544960B (zh) 一种提升孤岛微电网无功均分能力的分布式控制方法
CN104377711B (zh) 一种动态无功补偿方法
CN110098618B (zh) 一种自动潮流调整方法及系统
CN113162056B (zh) 特高压电网自动电压控制方法和装置
CN112865117B (zh) 新能源汇集直流闭锁后无功紧急控制方法、装置和系统
CN110661268B (zh) 一种新能源汇集站点的动态无功补偿需求确定方法及系统
CN103475035B (zh) 一种发电系统无功控制方法
CN108736483B (zh) 一种电网运行方式的确定方法及确定系统
CN107069702B (zh) 基于站内拓扑用于在线安全分析的大型风电场等值方法
CN109474027A (zh) 一种下垂控制方法、装置及变流器
CN116596199B (zh) 一种配电网线损影响因子的叠加分析方法和系统
CN113946985B (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