CN105914758B - 考虑分布式电源无功补偿成本的电网无功优化方法 - Google Patents

考虑分布式电源无功补偿成本的电网无功优化方法 Download PDF

Info

Publication number
CN105914758B
CN105914758B CN201610326301.2A CN201610326301A CN105914758B CN 105914758 B CN105914758 B CN 105914758B CN 201610326301 A CN201610326301 A CN 201610326301A CN 105914758 B CN105914758 B CN 105914758B
Authority
CN
China
Prior art keywords
mrow
msub
msup
mtd
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.)
Expired - Fee Related
Application number
CN201610326301.2A
Other languages
English (en)
Other versions
CN105914758A (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
North China Electric Power University
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power University
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, North China Electric Power University filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201610326301.2A priority Critical patent/CN105914758B/zh
Publication of CN105914758A publication Critical patent/CN105914758A/zh
Application granted granted Critical
Publication of CN105914758B publication Critical patent/CN105914758B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/1871—Methods for planning installation of shunt reactive power compensators
    • 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]
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明属于电网运行配电技术领域，尤其涉及一种考虑分布式电源无功补偿成本的电网无功优化方法，包括：分析分布式电源的无功服务成本，收集分布式电源技术参数，建立分布式电源无功服务成本模型；计及节点电压约束条件、可调变压器变比约束条件，兼顾经济性与安全性建立了以运行成本最小和电压偏差最小为目标函数的电网多目标无功优化模型；运用非支配排序遗传算法对无功优化模型进行求解，采用优劣解距离法在最优解集中求取分布式电源的最优无功出力。本发明对分布式电源的无功服务进行了经济补偿，在增加运营商自身收益的同时，使电网的无功优化调度方式更加灵活，大大减小电网的有功网损和电压偏差，有助于提高电网运行的经济性和安全性。

Description

考虑分布式电源无功补偿成本的电网无功优化方法
技术领域
本发明属于电网运行配电技术领域，尤其涉及一种考虑分布式电源无功补偿成本的电网无功优化方法。
背景技术
随着分布式电源的快速发展，分布式电源已成为电网调度环节中不可或缺的角色。在传统的电力工业结构模式下，发电商和电网公司有着共同的经济利益目标，因此传统无功优化模型中网损最小化目标能够满足他们的共同利益要求，从而使传统无功优化模型具有现实的可行性。但是在智能电网架构下，电力市场机制形成以后，能够提供无功辅助服务的分布式电源运营商和电网公司属于不同的经济利益主体，具有不同的利益需求，无功功率不再是有功功率的附属品，而是和有功功率一样成为可消费的商品，需要电网公司向分布式电源运营商支付无功服务的费用。
在不对分布式电源发出的无功进行经济补偿时，分布式电源对电网提供无功支撑不仅没有利益，而且有时增发无功还会降低有功输出从而影响有功收益，因此这种情况下分布式电源运营商并不愿意参与到电网的无功优化中来。为了充分调动分布式电源运营商参与电网无功优化的积极性，有必要从市场的角度出发考虑对分布式电源运营商提供无功服务的经济补偿问题。由于分布式电源提供的无功功率不再是无偿的，因此虽然分布式电源多发无功可以减小网络损耗，但是电网公司需要为多发的无功付费，而且电网的电压偏移也跟无功分布情况息息相关。因此有必要研究电网公司在竞争电力市场环境下，兼顾经济性与安全性，购买多种无功辅助服务的无功优化调度模型。
目前国内外学者针对电力市场环境下的无功优化研究基本分为两类：一类是无功辅助服务的定价研究，一类是考虑无功成本补偿的无功优化调度研究。关于第二类研究的一些文献侧重于不同类型分布式发电机组的成本分析，然后在成本分析的基础上建立电网的无功优化调度经济模型，而忽视了优化结果对电网安全性的影响。有些研究兼顾了无功优化的经济性和安全性，但其无功优化模型并没有考虑多种分布式电源的参与。
发明内容
为了解决目前电网无功优化方法无法同时兼顾经济性和安全性的问题，本发明提出了一种考虑分布式电源无功补偿成本的电网无功优化方法，包括如下步骤：
步骤1：分析分布式电源的无功服务成本，收集分布式电源技术参数，建立分布式电源无功服务成本模型；
步骤2：计及节点电压约束条件、可调变压器变比约束条件，兼顾经济性与安全性建立了以运行成本最小和电压偏差最小为目标函数的电网多目标无功优化模型；
步骤3：运用非支配排序遗传算法(NSGA-II)对步骤2的无功优化模型进行求解，采用优劣解距离法(TOPSIS法)在最优解集中求取分布式电源的最优无功出力。
所述步骤1建立分布式电源无功服务成本模型包括：光伏无功服务补偿成本和双馈风机无功服务补偿成本，具体计算方法包括：
1)光伏无功服务补偿成本的计算
光伏电站无功运行费用以逆变器的有功功率损耗费用为主，逆变器损耗可近似由二阶多项式函数表示：
PL,PV(SPV)＝l0+lV×SPV+lR×S2 PV (1)
式中：SPV是逆变器输出的视在功率；l0、lV、lR分别表示待机损耗、电压相关损耗(与电流I成比例)和电流相关损耗的损耗系数(与电流的平方I2成比例),这些损耗系数可由逆变器厂商提供的逆变器效率曲线估算得到；
光伏电站日间既可发出有功又可发出无功，发出无功而增加的损耗可以由逆变器发出和不发出无功时的损耗差获得，光伏电站夜间不发出有功但会由从电网吸收无功而产生损耗：
PPV是逆变器输出的有功功率，QPV是逆变器输出的无功功率；
因此对光伏电站无功服务的支付公式如下：
式中为光伏电站无功服务费用，FITPV是光伏的上网电价，指每kWh可再生能源发电获得的收益；
2)双馈风机无功服务补偿成本的计算
双馈风机的无功服务运行成本以机械和电力电子器件产生的额外有功损耗为主，假设网侧变换器始终工作在单位功率因数下，因此不考虑网侧变换器的功率损耗；机侧的变换器的功率损耗函数为：
PL,Inv,rotor＝a1,r+a2,r·|Sr|+a3,r|Sr|2 (4)
式中：PL,Inv,rotor为转子侧逆变器损耗；a1,r、a2,r、a3,r分别为逆变器的损耗系数；Sr为转子侧的视在功率；和分别为转子侧的电压和电流；
感应电动机的损耗包括绕组的铜损和磁化损耗；
Pr＝sPs (6)
Pelec＝Ps-Pr (7)
Qs＝Qelec (9)
式中：Ps和Qs分别为定子侧有功和无功功率；s为转差率；Pr为转子侧有功功率；Pelec和Qelec分别为风机的有功和无功输出；
风机输出有功功率与转差率之间的关系：
pu为标幺值；
PL,Cu＝3·Rs·|Is|2+3·Rr′·|Ir′|2 (14)
PL,Fe＝3·Rm·|Is+Ir′|2 (15)
式中：转换到定子侧的转子等效阻抗Zr′＝Rr′+jsXr′，Rr′为转换到定子侧的转子等效电阻，Xr′为转换到定子侧的转子等效电抗；转换到定子侧的转子励磁阻抗Zm′＝Rm+jsXm，Rm为励磁电阻，Xm为励磁电抗；定子侧励磁阻抗Zm＝Rm+jXm；定子侧阻抗Zs＝Rs+jXs，Rs为定子侧电阻，Xs为定子侧电抗；Ur′和Ir′分别为转换到定子侧的转子电压和电流；Us和Is分别为定子侧电压和电流；PL,Cu为风机铜损；PL,Fe为风机铁损；S* s为定子侧复功率，Ps为定子侧有功功率，Qs为定子侧无功功率；
因此整个双馈机组由多发无功而产生的损耗为以上三部分的和：
PL,WT(PWT,QWT)＝PL,Inv,rotor+PL,Cu+PL,Fe (16)
PWT是双馈机组输出的有功功率，QWT是双馈机组输出的无功功率；
配网需支付给双馈风机的无功服务费用包括风机提供无功的运行成本，计算公式如下式所示：
式中：为双馈风机的无功服务费用，FITWT是配网中的风电上网电价。
所述步骤2中的电网多目标无功优化模型包括：
目标函数1：电网运行成本最小；
运行成本包括为购买无功电量所支付的费用以及系统的网损成本，从而获得最大的综合经济效益:
式中ρp为有功清算价格，Ploss为配网有功损耗，为光伏电站无功服务费用，为双馈风机的无功服务费用；
目标函数2：电压偏移最小
电压偏移最小的目标函数就是将各节点的电压与理想电压值范围的偏移量总和最小化，即提高负荷节点的电压水平：
式中，Ui为i节点的电压，为节点i的理想电压值，ΔUi为节点i可以接受的最大电压偏离范围，ΔU为各节点的电压与理想电压值范围的偏移量总和，中间变量N为节点总数；
约束条件：
1)潮流等式约束：
式中，PGi为节点i的有功出力，PDi为节点i的有功负荷，QGi为节点i的无功出力，QDi为节点i的无功负荷，QCi为节点i的无功功率补偿，Ui、Uj分别为节点i、j的电压，Gij、Bij、δij为节点i、j之间的电导、电纳和相角差；
2)不等式约束：
式中，Uimax、Uimin为节点i电压上下限，QCmin、QCmax为并联电容器组的无功出力上下限，QC为并联电容器组的无功出力；
3)光伏电池的无功出力约束：
光伏电池的无功出力极限由并网逆变器的容量决定：
式中，QPV为光伏电池在有功出力为P的情况下能够吸收和发出的最大无功功率，S为光伏电池的视在功率；
4)双馈风机的无功出力约束：
双馈风机能够吸收的最大无功功率值取决于双馈风机的定子电流极限：
式中：U为双馈风机节点出口处电压，为定子电流极限值，PWT为风机输出的有功功率，即为Pelec，s为转速差；
双馈风机能够发出的最大无功出力取决于转子电流极限：
式中，Irmax为转子电流极限值，γ为功率因数角，Zs为定子侧阻抗，Zm为定子侧励磁阻抗，Xm为励磁电抗，Xs为定子侧电抗；
Rs为定子侧电阻，Rm为励磁电阻。
所述步骤3中运用非支配排序遗传算法对步骤2的无功优化模型进行求解，具体步骤如下：
步骤301、输入网络参数、负荷预测和分布式电源出力预测参数，进行参数初始化及种群初始化；
步骤302、进行选择、交叉、变异操作，以运行成本最小和电压偏差最小为目标函数，计及节点电压约束条件、可调变压器变比约束条件，求得帕累托最优解集。
步骤303、采用TOPSIS法求得最优折衷解，选出最合适的分布式电源无功出力方案。
本发明的有益效果是本发明提出的考虑分布式电源无功补偿成本的电网无功优化方法，是在智能电网电力市场环境下，兼顾电网安全性和经济性，在分布式电源提供无功辅助服务的前提下建立无功优化模型，提高分布式电源运营商参与电网无功优化的积极性，并改善电网经济效益。
附图说明
图1为求解无功优化模型算法流程图。
图2为并入2个分布式电源的33节点电网系统结构图。
图3为非支配排序遗传算法得到帕累托解集示意图。
图4为两种场景下双馈风机无功出力结果对比图。
图5为两种场景下光伏无功出力结果对比图。
图6为两种场景下电压水平对比图。
具体实施方式
下面结合附图，详细说明实施方案。
本发明提供一种考虑分布式电源无功补偿成本的电网无功优化方法，包括如下步骤：
步骤1：分析分布式电源的无功服务成本，收集分布式电源技术参数，建立分布式电源无功服务成本模型；
步骤2：计及节点电压约束、可调变压器变比约束等约束条件，兼顾经济性与安全性建立了以运行成本最小和电压偏差最小为目标函数的电网多目标无功优化模型；
步骤3：运用非支配排序遗传算法(NSGA-II)对步骤2的无功优化模型进行求解，采用优劣解距离法(TOPSIS法)在最优解集中求取分布式电源的最优无功出力。
下面结合附图，对本发明作详细说明。
所述步骤1建立分布式电源无功服务成本模型，具体包括：
(1)光伏无功服务补偿成本
光伏电站无功运行费用主要是逆变器的有功功率损耗费用，逆变器损耗可近似由二阶多项式函数表示：
PL,PV(SPV)＝l0+lV×SPV+lR×S2 PV (1)
式中：SPV是逆变器输出的视在功率；l0、lV、lR分别表示待机损耗、电压相关损耗(与电流I成比例)和电流相关损耗的损耗系数(与电流的平方I2成比例)。这些参数可由逆变器厂商提供的逆变器效率曲线估算得到。
由于发出无功而增加的损耗可以由逆变器发出和不发出无功时的损耗差获得，光伏电站在夜晚有功输出为0，因此可以将光伏电站的运行模式分为两种：日间模式和夜间模式。光伏电站处于日间模式时既可发出有功又可发出无功，增加的损耗为式(2)中PPV≠0所对应的损耗。光伏电站处于夜间模式时不发出有功，但是会由于从电网吸收无功功率而产生损耗，如式(2)所示。
因此，对光伏电站无功服务的支付公式如下：
式中FITPV是光伏的上网电价，指每kWh可再生能源发电获得的收益。可再生能源的上网电价取决于很多因素，包括电厂类型和装机容量、合同时间等。在本发明模型中采用的是上网电价，功率损耗按照电能购买价格报销，从而反映从电网吸收功率的成本。式(3)中还包含了光伏电站提供无功的运行成本，当光伏电站停产时或者不允许向电网注入无功时，对光伏电站的无功服务支付费用则为0。
(2)双馈风机无功服务补偿成本
双馈风机的有功和无功输出并不存在耦合关系，可各自独立控制。双馈风机的无功服务运行成本主要来源于机械和电力电子器件产生的额外有功损耗。假设网侧变换器始终工作在单位功率因数下，因此不考虑网侧变换器的功率损耗。机侧的变换器的功率损耗函数为：
PL,Inv,rotor＝a1,r+a2,r·|Sr|+a3,r|Sr|2 (4)
式中：PL,Inv,rotor为转子侧逆变器损耗；a1,r、a2,r、a3,r分别为逆变器的损耗系数；Sr为转子侧的视在功率；和分别为转子侧的电压和电流。
感应电动机的损耗包括绕组的铜损和磁化损耗。
Pr＝sPs (6)
Pelec＝Ps-Pr (7)
Qs＝Qelec (9)
式中：Ps和Qs分别为定子侧有功和无功功率；s为转差率；Pr为转子侧有功功率；Pelec和Qelec分别为风机的有功和无功输出。
输出功率与转差率之间的关系：
pu为标幺值；
PL,Cu＝3·Rs·|Is|2+3·Rr′·|Ir′|2 (14)
PL,Fe＝3·Rm·|Is+Ir′|2 (15)
式中：Zr′＝Rr′+jsXr′，Rr′为转换到定子侧的转子等效电阻，Xr′为转换到定子侧的转子等效电抗；Zm′＝Rm+jsXm，Rm为励磁电阻，Xm为励磁电抗；Zm＝Rm+jXm；Zs＝Rs+jXs，Rs为定子侧电阻，Xs为定子侧电抗；Ur′和Ir′分别为转换到定子侧的转子电压和电流；Us和Is分别为定子侧电压和电流；PL,Cu为风机铜损；PL,Fe为风机铁损。
因此整个双馈机组由多发无功而产生的损耗为以上三部分的和：
PL,WT(PWT,QWT)＝PL,Inv,rotor+PL,Cu+PL,Fe (16)
配网需支付给双馈风机的无功服务费用包括风机提供无功的运行成本，计算公式如下式所示：
式中：FITWT是配网中的风电上网电价。
所述步骤2中考虑分布式电源无功补偿成本的电网无功优化模型的建立所涉及运行成本最小、电压偏移最小的目标函数以及一系列约束条件，以兼顾电网运行的安全性和经济性：
目标函数1：电网运行成本最小
在电力市场化的环境下所进行的无功优化必须考虑无功补偿所带来的一系列费用，电网公司由原来所关注的系统网损最小转变为综合成本的最小化，所以本发明无功优化的第一个目标是使电网公司运行成本最小，运行成本包括为购买无功电量所支付的费用以及系统的网损成本，从而获得最大的综合经济效益:
式中ρp为有功清算价格，Ploss为配网有功损耗。
目标函数2：电压偏移最小
电压偏移最小的目标函数就是将各节点的电压与理想电压值范围的偏移量总和最小化，即提高负荷节点的电压水平。
式中，Ui为i节点的电压，为节点i的理想电压值，ΔUi为节点i可以接受的最大电压偏离范围。
约束条件：
1)潮流等式约束：
式中，PGi为节点i的有功出力，PDi为节点i的有功负荷，QGi为节点i的无功出力，QDi为节点i的无功负荷，QCi为节点i的无功功率补偿，Ui、Uj分别为节点i、j的电压，Gij、Bij、δij为节点i、j之间的电导、电纳和相角差；
2)不等式约束：
式中，Uimax、Uimin为节点电压上下限，QCmin、QCmax为并联电容器组的无功出力上下限。
除此之外，还有光伏和双馈风电机组的无功出力范围约束。
3)光伏电池的无功出力约束：
现有光伏电站一般按照单位功率因数运行，忽略了光伏逆变器的无功出力能力。光伏电池的无功出力极限由并网逆变器的容量决定：
式中，QPV为光伏电池在有功出力为P的情况下能够吸收和发出的最大无功功率。
4)双馈风机的无功出力约束：
双馈风机能够吸收的最大无功功率值取决于双馈风机的定子电流极限：
式中：U为双馈风机节点出口处电压，为定子电流极限值，PWT为风机输出的有功功率，即为Pelec。
双馈风机能够发出的最大无功出力取决于转子电流极限：
式中，Irmax为转子电流极限值，γ为功率因数角。
所述步骤3中运用非支配排序遗传算法(NSGA-II)求解无功优化模型的优化流程，如图1所示，具体步骤如下：
(1)输入网络参数、负荷预测和分布式电源出力预测参数，参数初始化及种群初始化。
(2)进行选择、交叉、变异等操作，以运行成本最小和电压偏差最小为目标函数，计及节点电压约束、可调变压器变比约束等约束条件，由智能算法求得帕累托最优解集。
(3)采用TOPSIS法求得最优折衷解，选出最合适的分布式电源无功出力方案。
具体算例如下：
采用IEEE33节点配电系统作为算例，该配电系统共有33个节点，如图2所示。参与无功辅助服务市场的无功源包括风机、光伏电池、微型燃气轮机和并联电容器组，风机安装位置为13节点，额定容量为2MW，光伏电站的安装位置为31节点，安装容量为1MW，并联电容器组安装位置为3节点和26节点，每组额定容量为10kvar，每个安装位置各安装30组。假设分布式电源按照发多少用多少的原则进行发电，本发明所做的无功优化调度是在有功调度结果的基础上进行的。根据现有国家规定，光伏上网价格0.92元/kW·h，风电上网价格0.52元/kW·h，有功清算的价格为0.44元/kW·h。
光伏的有功损耗参数为l0＝4.5kW，lV＝0.00352kW/kVA，lR＝0.0000125kW/kVA2。双馈风机的各个参数见下表：
表1风机参数
对于上述优化模型，本发明使用MATLAB编程求解。NSGA-II算法中参数设置为：种群数量50，最大迭代次数200,交叉率0.8，变异率0.3。
表2给出了在考虑向分布式电源运营商提供无功经济补偿后各个时段的无功优化调度结果。根据结果可以看出，双馈风机以及光伏电池都向电网提供了较多的无功功率，为电网电压提供了支撑。
表2含分布式电源的无功优化调度结果
图3为使用NSGA-II在第21时段对本发明模型进行优化计算所得到的Pareto前沿.通过NSGA-II算法所得到的Pareto前沿分布均匀且稳定，说明了该算法对于求解本发明多目标优化模型的可行性。Pareto前沿分布图横轴为计及分布式电源无功服务的电网的运行成本，纵轴为电网电压偏差，随着电网运行成本的增加，电网的电压偏差会逐渐减小，前沿分布图则为这两个冲突的目标函数的折衷提供了丰富的决策信息。同时也证明了本发明所提出的成本计算方法能够很好的描述无功辅助服务成本信息。
为了研究不同分布式无功补偿方式对结果影响，假设两种场景对本发明所提模型进行分析：
场景1：分布式电源按照恒定功率因数运行(其中双馈风机功率因数范围为±0.95，光伏电池功率因数为1)；
场景2：依照本发明模型，双馈风机以及光伏电池的无功出力在其无功出力极限范围内波动。
同样以有功损耗最小和电压偏差最小为目标函数对上述两种场景进行优化计算。
图4和图5分别给出了两种场景下双馈机组和光伏电池的无功出力结果对比图。根据图中结果可以明显看出将分布式电源按照无功出力极限参与无功调度时，分布式电源可以发出更多的无功，更能有效的支撑电网的运行。而且在本发明模型中还对分布式电源的无功服务进行了经济补偿，在这种情况下分布式电源运营商更愿意参与到电网的无功优化中来，从而充分利用了分布式电源的无功调节能力，在增加分布式电源运营商自身收益的同时，使电网的无功优化调度方式更加灵活。
优化后两个场景下各个目标函数值的大小见表3：
表3两种场景下优化目标对比
从表3可以看出，场景2的有功网损大小和电压水平都要优于场景1的优化结果，说明调动分布式电源运营商积极参与电网无功优化后，可以大大减小电网的有功网损，同时减小了电压偏差，有助于提高电网运行的经济性和安全性。
两种场景下电网各个节点的电压水平对比见图6。从图6中可以看出，场景2下的电压曲线明显比场景1下的电压曲线更平缓，各节点电压基本在理想电压附近波动，分布更为均匀。而场景1下的各节点电压波动幅度较大，不利于电网的安全稳定运行。
对于表3中的第三个运行成本目标，虽然在场景2下电网需要向分布式电源运营商多支付一部分的无功补偿费用，但与网损费用中和后总费用增加的并不多，因此相较于第三个运行成本的增加，场景2下在网损和电压偏差方面的益处要更加突出，因此向分布式电源运营商提供无功辅助服务经济补偿，调动其参与电网无功优化是有必要且有益的。
此实施例仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

Claims (2)

1.一种考虑分布式电源无功补偿成本的电网无功优化方法，其特征在于，包括如下步骤：
步骤1：分析分布式电源的无功服务成本，收集分布式电源技术参数，建立分布式电源无功服务成本模型；
步骤2：计及节点电压约束条件、可调变压器变比约束条件，兼顾经济性与安全性建立了以运行成本最小和电压偏差最小为目标函数的电网多目标无功优化模型；
步骤3：运用非支配排序遗传算法对步骤2的无功优化模型进行求解，采用优劣解距离法在最优解集中求取分布式电源的最优无功出力；
所述步骤1建立分布式电源无功服务成本模型包括：光伏无功服务补偿成本和双馈风机无功服务补偿成本，具体计算方法包括：
1)光伏无功服务补偿成本的计算
光伏电站无功运行费用以逆变器的有功功率损耗费用为主，逆变器损耗可近似由二阶多项式函数表示：
PL,PV(SPV)＝l0+lV×SPV+lR×S2 PV (1)
式中：SPV是逆变器输出的视在功率；l0、lV、lR分别表示待机损耗、电压相关损耗和电流相关损耗的损耗系数,这些损耗系数可由逆变器厂商提供的逆变器效率曲线估算得到；
光伏电站日间既可发出有功又可发出无功，发出无功而增加的损耗可以由逆变器发出和不发出无功时的损耗差获得，光伏电站夜间不发出有功但会由从电网吸收无功而产生损耗：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;P</mi> <mrow> <mi>L</mi> <mo>.</mo> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>Q</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>Q</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>Q</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>=</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>&amp;NotEqual;</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>=</mo> <mn>0</mn> <mo>,</mo> <msub> <mi>Q</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
PPV是逆变器输出的有功功率，QPV是逆变器输出的无功功率；
因此对光伏电站无功服务的支付公式如下：
<mrow> <msubsup> <mi>C</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> <mi>Q</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>Q</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>FIT</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>&amp;times;</mo> <msub> <mi>&amp;Delta;P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>Q</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中为光伏电站无功服务费用，FITPV是光伏的上网电价，指每kWh可再生能源发电获得的收益；
2)双馈风机无功服务补偿成本的计算
双馈风机的无功服务运行成本以机械和电力电子器件产生的额外有功损耗为主，假设网侧变换器始终工作在单位功率因数下，因此不考虑网侧变换器的功率损耗；机侧的变换器的功率损耗函数为：
<mrow> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>I</mi> <mi>n</mi> <mi>v</mi> <mo>,</mo> <mi>r</mi> <mi>o</mi> <mi>t</mi> <mi>o</mi> <mi>r</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>a</mi> <mrow> <mn>1</mn> <mo>,</mo> <mi>r</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>a</mi> <mrow> <mn>2</mn> <mo>,</mo> <mi>r</mi> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mo>|</mo> <msub> <mi>S</mi> <mi>r</mi> </msub> <mo>|</mo> <mo>+</mo> <msub> <mi>a</mi> <mrow> <mn>3</mn> <mo>,</mo> <mi>r</mi> </mrow> </msub> <mo>|</mo> <msub> <mi>S</mi> <mi>r</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>S</mi> <mi>r</mi> </msub> <mo>=</mo> <mn>3</mn> <mo>&amp;times;</mo> <msub> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>r</mi> </msub> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>r</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
式中：PL,Inv,rotor为转子侧逆变器损耗；a1,r、a2,r、a3,r分别为逆变器的损耗系数；Sr为转子侧的视在功率；和分别为转子侧的电压和电流；
感应电动机的损耗包括绕组的铜损和磁化损耗；
Pr＝sPs (6)
Pelec＝Ps-Pr (7)
<mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>1</mn> <mo>-</mo> <mi>s</mi> </mrow> </mfrac> <msub> <mi>P</mi> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
Qs＝Qelec (9)
式中：Ps和Qs分别为定子侧有功和无功功率；s为转差率；Pr为转子侧有功功率；Pelec和Qelec分别为风机的有功和无功输出；
风机输出有功功率与转差率之间的关系：
<mrow> <mi>s</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>0.3</mn> </mtd> <mtd> <mrow> <mn>0</mn> <mi>p</mi> <mi>u</mi> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> <mi>c</mi> </mrow> </msub> <mo>&amp;le;</mo> <mn>0.1</mn> <mi>p</mi> <mi>u</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mn>13</mn> <mn>30</mn> </mfrac> <mo>-</mo> <mfrac> <mn>4</mn> <mn>3</mn> </mfrac> <msub> <mi>P</mi> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> <mi>c</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mn>0.1</mn> <mi>p</mi> <mi>u</mi> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> <mi>c</mi> </mrow> </msub> <mo>&amp;le;</mo> <mn>0.4</mn> <mi>p</mi> <mi>u</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>0.1</mn> </mrow> </mtd> <mtd> <mrow> <mn>0.4</mn> <mi>p</mi> <mi>u</mi> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> <mi>c</mi> </mrow> </msub> <mo>&amp;le;</mo> <mn>1</mn> <mi>p</mi> <mi>u</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
pu为标幺值；
<mrow> <mtable> <mtr> <mtd> <mrow> <msup> <msub> <mi>U</mi> <mi>r</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <msup> <msub> <mi>Z</mi> <mi>r</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>+</mo> <msup> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <msub> <mi>Z</mi> <mi>m</mi> </msub> </mfrac> <msub> <mi>U</mi> <mi>s</mi> </msub> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;lsqb;</mo> <msup> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <mfrac> <mrow> <mo>(</mo> <msup> <msub> <mi>Z</mi> <mi>r</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>+</mo> <msup> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>)</mo> <mo>(</mo> <msub> <mi>Z</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>Z</mi> <mi>m</mi> </msub> </mfrac> <mo>&amp;rsqb;</mo> <msub> <mi>I</mi> <mi>s</mi> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msup> <msub> <mi>I</mi> <mi>r</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>Z</mi> <mi>m</mi> </msub> </mfrac> <msub> <mi>U</mi> <mi>s</mi> </msub> <mo>-</mo> <mfrac> <mrow> <msub> <mi>Z</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> </mrow> <msub> <mi>Z</mi> <mi>m</mi> </msub> </mfrac> <msub> <mi>I</mi> <mi>s</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>I</mi> <mi>s</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <msubsup> <mi>S</mi> <mi>s</mi> <mo>*</mo> </msubsup> <mrow> <mn>3</mn> <mo>&amp;CenterDot;</mo> <msub> <mi>U</mi> <mi>s</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>jQ</mi> <mi>s</mi> </msub> </mrow> <mrow> <mn>3</mn> <mo>&amp;CenterDot;</mo> <msub> <mi>U</mi> <mi>s</mi> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mfrac> <mn>1</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>s</mi> <mo>)</mo> </mrow> </mfrac> <msub> <mi>P</mi> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> <mi>c</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>jQ</mi> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> <mi>c</mi> </mrow> </msub> </mrow> <mrow> <mn>3</mn> <mo>&amp;CenterDot;</mo> <msub> <mi>U</mi> <mi>s</mi> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>C</mi> <mi>u</mi> </mrow> </msub> <mo>=</mo> <mn>3</mn> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <mo>&amp;CenterDot;</mo> <mo>|</mo> <msub> <mi>I</mi> <mi>s</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>+</mo> <mn>3</mn> <mo>&amp;CenterDot;</mo> <msup> <msub> <mi>R</mi> <mi>r</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;CenterDot;</mo> <mo>|</mo> <msup> <msub> <mi>I</mi> <mi>r</mi> </msub> <mo>&amp;prime;</mo> </msup> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>F</mi> <mi>e</mi> </mrow> </msub> <mo>=</mo> <mn>3</mn> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>m</mi> </msub> <mo>&amp;CenterDot;</mo> <mo>|</mo> <msub> <mi>I</mi> <mi>s</mi> </msub> <mo>+</mo> <msup> <msub> <mi>I</mi> <mi>r</mi> </msub> <mo>&amp;prime;</mo> </msup> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
式中：转换到定子侧的转子等效阻抗Zr′＝Rr′+jsXr′，Rr′为转换到定子侧的转子等效电阻，Xr′为转换到定子侧的转子等效电抗；转换到定子侧的转子励磁阻抗Zm′＝Rm+jsXm，Rm为励磁电阻，Xm为励磁电抗；定子侧励磁阻抗Zm＝Rm+jXm；定子侧阻抗Zs＝Rs+jXs，Rs为定子侧电阻，Xs为定子侧电抗；Ur′和Ir′分别为转换到定子侧的转子电压和电流；Us和Is分别为定子侧电压和电流；PL,Cu为风机铜损；PL,Fe为风机铁损；S* s为定子侧复功率，Ps为定子侧有功功率，Qs为定子侧无功功率；
因此整个双馈机组由多发无功而产生的损耗为以上三部分的和：
PL,WT(PWT,QWT)＝PL,Inv,rotor+PL,Cu+PL,Fe (16)
ΔPL，WT(PWT,QWT)＝
PL,WT(PWT,QWT)-PL,WT(PWT,QWT＝0) (17)
PWT是双馈机组输出的有功功率，QWT是双馈机组输出的无功功率；
配网需支付给双馈风机的无功服务费用包括风机提供无功的运行成本，计算公式如下式所示：
<mrow> <msubsup> <mi>C</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> <mi>Q</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>Q</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>FIT</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msub> <mo>&amp;times;</mo> <msub> <mi>&amp;Delta;P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>W</mi> <mi>T</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>Q</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
式中：为双馈风机的无功服务费用，FITWT是配网中的风电上网电价；
所述步骤2中的电网多目标无功优化模型包括：
目标函数1：电网运行成本最小；
运行成本包括为购买无功电量所支付的费用以及系统的网损成本，从而获得最大的综合经济效益:
<mrow> <msub> <mi>minF</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>&amp;rho;</mi> <mi>p</mi> </msub> <msub> <mi>P</mi> <mrow> <mi>l</mi> <mi>o</mi> <mi>s</mi> <mi>s</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>C</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> <mi>Q</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>Q</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>C</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> <mi>Q</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>Q</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
式中ρp为有功清算价格，Ploss为配网有功损耗，为光伏电站无功服务费用，为双馈风机的无功服务费用；
目标函数2：电压偏移最小
电压偏移最小的目标函数就是将各节点的电压与理想电压值范围的偏移量总和最小化，即提高负荷节点的电压水平：
<mrow> <msub> <mi>minF</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&amp;Delta;</mi> <mi>U</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mfrac> <mrow> <mi>&amp;psi;</mi> <mrow> <mo>(</mo> <mo>|</mo> <msub> <mi>U</mi> <mi>i</mi> </msub> <mo>-</mo> <msubsup> <mi>U</mi> <mi>i</mi> <mrow> <mi>i</mi> <mi>d</mi> <mi>e</mi> <mi>a</mi> <mi>l</mi> </mrow> </msubsup> <mo>|</mo> <mo>-</mo> <msub> <mi>&amp;Delta;U</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> <msub> <mi>U</mi> <mi>i</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>&amp;psi;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>x</mi> <mo>&lt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>x</mi> </mtd> <mtd> <mrow> <mi>x</mi> <mo>&gt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
式中，Ui为i节点的电压，为节点i的理想电压值，ΔUi为节点i可以接受的最大电压偏离范围，ΔU为各节点的电压与理想电压值范围的偏移量总和，中间变量N为节点总数；
约束条件：
1)潮流等式约束：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>D</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>U</mi> <mi>i</mi> </msub> <munderover> <mo>&amp;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&amp;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&amp;delta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>Q</mi> <mrow> <mi>D</mi> <mi>i</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>Q</mi> <mrow> <mi>C</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>U</mi> <mi>i</mi> </msub> <munderover> <mo>&amp;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&amp;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&amp;delta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>
式中，PGi为节点i的有功出力，PDi为节点i的有功负荷，QGi为节点i的无功出力，QDi为节点i的无功负荷，QCi为节点i的无功功率补偿，Ui、Uj分别为节点i、j的电压，Gij、Bij、δij为节点i、j之间的电导、电纳和相角差；
2)不等式约束：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mrow> <mi>C</mi> <mi>min</mi> </mrow> </msub> <mo>&amp;le;</mo> <msub> <mi>Q</mi> <mi>C</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>Q</mi> <mrow> <mi>C</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>min</mi> </mrow> </msub> <mo>&amp;le;</mo> <msub> <mi>U</mi> <mi>i</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
式中，Uimax、Uimin为节点i电压上下限，QCmin、QCmax为并联电容器组的无功出力上下限，QC为并联电容器组的无功出力；
3)光伏电池的无功出力约束：
光伏电池的无功出力极限由并网逆变器的容量决定：
<mrow> <msup> <mi>Q</mi> <mrow> <mi>P</mi> <mi>V</mi> </mrow> </msup> <mo>=</mo> <mo>&amp;PlusMinus;</mo> <msqrt> <mrow> <msup> <mi>S</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mi>P</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>
式中，QPV为光伏电池在有功出力为P的情况下能够吸收和发出的最大无功功率，S为光伏电池的视在功率；
4)双馈风机的无功出力约束：
双馈风机能够吸收的最大无功功率值取决于双馈风机的定子电流极限：
<mrow> <msubsup> <mi>Q</mi> <mrow> <mi>a</mi> <mi>b</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msubsup> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mo>|</mo> <mi>U</mi> <mo>|</mo> <msubsup> <mi>I</mi> <mrow> <mi>s</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>P</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msub> <mrow> <mn>1</mn> <mo>-</mo> <mi>s</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow>
式中：U为双馈风机节点出口处电压，为定子电流极限值，PWT为风机输出的有功功率，即为Pelec，s为转速差；
双馈风机能够发出的最大无功出力取决于转子电流极限：
<mrow> <msubsup> <mi>Q</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msubsup> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>|</mo> <mi>U</mi> <msup> <mo>|</mo> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>X</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>|</mo> <msub> <mi>Z</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>|</mo> <mi>U</mi> <mo>|</mo> <mo>|</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>|</mo> <msub> <mi>I</mi> <mrow> <mi>r</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;gamma;</mi> </mrow> <mrow> <mo>|</mo> <msub> <mi>Z</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>26</mn> <mo>)</mo> </mrow> </mrow>
式中，Irmax为转子电流极限值，γ为功率因数角，Zs为定子侧阻抗，Zm为定子侧励磁阻抗，Xm为励磁电抗，Xs为定子侧电抗；
<mrow> <mi>&amp;gamma;</mi> <mo>=</mo> <msup> <mi>cos</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mrow> <mi>W</mi> <mi>T</mi> </mrow> </msub> <mo>|</mo> <msub> <mi>Z</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>+</mo> <mo>|</mo> <mi>U</mi> <msup> <mo>|</mo> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>s</mi> <mo>)</mo> <mo>|</mo> <mi>U</mi> <mo>|</mo> <mo>|</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>|</mo> <mo>|</mo> <msub> <mi>Z</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mi>m</mi> </msub> <mo>|</mo> <msub> <mi>I</mi> <mrow> <mi>r</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>
Rs为定子侧电阻，Rm为励磁电阻。
2.根据权利要求1所述方法，其特征在于，所述步骤3中运用非支配排序遗传算法对步骤2的无功优化模型进行求解，具体步骤如下：
步骤301、输入网络参数、负荷预测和分布式电源出力预测参数，进行参数初始化及种群初始化；
步骤302、进行选择、交叉、变异操作，以运行成本最小和电压偏差最小为目标函数，计及节点电压约束条件、可调变压器变比约束条件，求得帕累托最优解集；
步骤303、采用TOPSIS法求得最优折衷解，选出最合适的分布式电源无功出力方案。
CN201610326301.2A 2016-05-17 2016-05-17 考虑分布式电源无功补偿成本的电网无功优化方法 Expired - Fee Related CN105914758B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610326301.2A CN105914758B (zh) 2016-05-17 2016-05-17 考虑分布式电源无功补偿成本的电网无功优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610326301.2A CN105914758B (zh) 2016-05-17 2016-05-17 考虑分布式电源无功补偿成本的电网无功优化方法

Publications (2)

Publication Number Publication Date
CN105914758A CN105914758A (zh) 2016-08-31
CN105914758B true CN105914758B (zh) 2018-04-10

Family

ID=56749143

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610326301.2A Expired - Fee Related CN105914758B (zh) 2016-05-17 2016-05-17 考虑分布式电源无功补偿成本的电网无功优化方法

Country Status (1)

Country Link
CN (1) CN105914758B (zh)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107196315A (zh) * 2017-06-09 2017-09-22 华南理工大学 含光储系统配电网的动态拓展无功优化控制方法
CN107611993B (zh) * 2017-08-28 2019-09-10 国网山东省电力公司电力科学研究院 一种适用于特高压半波长输电系统的无功优化方法
CN107480837B (zh) * 2017-09-29 2020-10-09 山东理工大学 一种基于gps同步定频的孤岛运行微电网协调控制方法
US10570882B2 (en) 2017-11-13 2020-02-25 General Electric Company Dynamic active and reactive power capability for wind farms
CN108646552B (zh) * 2018-04-16 2020-12-11 杭州电子科技大学信息工程学院 基于遗传算法的天然气分布式能源机组参数的多目标优化方法
CN110738344B (zh) * 2018-07-20 2022-07-05 中国农业大学 电力系统负荷预测的分布式无功优化方法及装置
CN109149564A (zh) * 2018-08-31 2019-01-04 国网浙江省电力有限公司经济技术研究院 一种交直流混合配电网分布式电源优化配置方法
CN110277799A (zh) * 2019-06-17 2019-09-24 国网江苏省电力有限公司淮安供电分公司 计及控制设备经济性及动作次数限制的电网无功优化方法
CN110245799B (zh) * 2019-06-18 2023-04-07 国网江西省电力有限公司经济技术研究院 考虑负荷柔性需求的配电网架结构过渡的多目标规划方法
CN110323785B (zh) * 2019-07-26 2022-02-01 东北电力大学 源网荷储互动的多电压等级直流配电网优化调度方法
CN111682570B (zh) * 2020-06-30 2022-02-15 广东电网有限责任公司电力科学研究院 一种电力系统变压器运行效率寻优方法、系统及设备
CN111987747B (zh) * 2020-07-14 2022-02-15 湖南大学 一种大型双馈风电场的无功优化控制方法
CN112003270B (zh) * 2020-08-06 2022-03-18 国网浙江省电力有限公司衢州供电公司 一种应用于电力市场环境下的无功辅助服务优化调度方法
CN112366725A (zh) * 2020-12-03 2021-02-12 国网青海省电力公司果洛供电公司 配电网用静止无功发生器容量及安装位置选择方法及系统
CN112865080B (zh) * 2021-01-15 2022-12-16 河海大学 一种电池储能参与电网调压辅助服务的补偿方法
CN112928760A (zh) * 2021-02-04 2021-06-08 中国电建集团青海省电力设计院有限公司 计及电网负荷波动的风电场主动无功补偿控制方法
CN114552653B (zh) * 2022-03-21 2024-03-26 合肥工业大学 基于智能算法的并网逆变器系统控制参数多目标优化方法
CN116683471B (zh) * 2023-04-28 2024-06-11 国网河北省电力有限公司电力科学研究院 无功功率补偿资源的配置方法、装置以及设备
CN116362141B (zh) * 2023-06-01 2023-08-01 中诚华隆计算机技术有限公司 基于量子遗传算法的配电网无功优化方法及装置
CN117424294B (zh) * 2023-12-18 2024-03-01 国网辽宁省电力有限公司经济技术研究院 一种配电网高效无功规划方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3039771B1 (en) * 2013-08-28 2018-05-09 Robert Bosch GmbH System and method for energy asset sizing and optimal dispatch
CN104135025B (zh) * 2014-05-30 2017-01-18 国家电网公司 基于模糊粒子群算法的微网并网经济优化方法

Also Published As

Publication number Publication date
CN105914758A (zh) 2016-08-31

Similar Documents

Publication Publication Date Title
CN105914758B (zh) 考虑分布式电源无功补偿成本的电网无功优化方法
Pan et al. Optimal design and operation of multi-energy system with load aggregator considering nodal energy prices
Tavakoli et al. Impacts of grid integration of solar PV and electric vehicle on grid stability, power quality and energy economics: A review
CN110244566B (zh) 计及柔性负荷的冷热电联供系统容量优化配置方法
CN107546743B (zh) 一种辐射状配电网分布式优化潮流方法
Caldon et al. Optimal control of a distribution system with a virtual power plant
CN102541621B (zh) 风光储联合发电系统的仿真等值方法
CN110601198B (zh) 计及谐波和电压不平衡约束的混合微电网优化运行方法
CN109740786A (zh) 一种多端柔性互联配电网规划方法及系统
CN108304972B (zh) 一种基于供需互动和dg运行特性的主动配电网网架规划方法
CN111769543B (zh) 一种含多微网的区域配电网自律协同运行优化方法
CN112163968B (zh) 一种综合能源储能系统优化调度方法
Liu et al. Energy storage in China: Development progress and business model
CN112149914A (zh) 一种多约束条件下电力市场资源优化配置的方法
CN111952958A (zh) 一种考虑控制模式转换的配电网柔性软开关优化配置方法
CN111445092A (zh) 基于改进jaya算法的多微电网优化方法
CN112311017A (zh) 一种虚拟电厂与主网最优协同调度方法
CN116345466A (zh) 一种考虑配网重构的主动配电网两阶段潮流优化方法
CN114884110A (zh) 一种源网荷多约束条件下的电力系统储能优化运行方法
CN115149540A (zh) 一种基于分区自治的多区域虚拟电厂协同优化运行方法
Jayawardena et al. Microgrid capability diagram: A tool for optimal grid-tied operation
CN116707023A (zh) 基于源荷相关性聚类的主动配电网分层分区综合优化方法
CN112186801A (zh) 一种采用p2p能源交易机制提高农村配网分布式光伏并网容量的方法
Kermani et al. Optimal operation of a real power hub based on PV/FC/GenSet/BESS and demand response under uncertainty
CN110783928B (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180410

Termination date: 20190517

CF01 Termination of patent right due to non-payment of annual fee