CN109995050B - 计及电压暂降经济损失不确定性的补偿设备优化配置方法 - Google Patents
计及电压暂降经济损失不确定性的补偿设备优化配置方法 Download PDFInfo
- Publication number
- CN109995050B CN109995050B CN201910403118.1A CN201910403118A CN109995050B CN 109995050 B CN109995050 B CN 109995050B CN 201910403118 A CN201910403118 A CN 201910403118A CN 109995050 B CN109995050 B CN 109995050B
- Authority
- CN
- China
- Prior art keywords
- voltage
- fault
- compensation
- node
- voltage sag
- 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 40
- 238000005457 optimization Methods 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 19
- 238000009826 distribution Methods 0.000 claims abstract description 16
- 230000003068 static effect Effects 0.000 claims abstract description 10
- 230000001360 synchronised effect Effects 0.000 claims abstract description 10
- 238000009434 installation Methods 0.000 claims abstract description 4
- 230000035772 mutation Effects 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000010367 cloning Methods 0.000 claims description 15
- 230000016507 interphase Effects 0.000 claims description 12
- 230000005764 inhibitory process Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000003213 activating effect Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000002649 immunization Methods 0.000 claims description 3
- 230000003053 immunization Effects 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 claims description 3
- 238000000819 phase cycle Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000013178 mathematical model Methods 0.000 abstract 1
- 230000002265 prevention Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/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/1842—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators with stepless control wherein at least one reactive element is actively controlled by a bridge converter, e.g. active filters
-
- 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/10—Flexible AC transmission systems [FACTS]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种计及电压暂降经济损失不确定性的补偿设备优化配置方法,计及不同故障时配电网敏感负荷电压暂降特征的不确定性,建立电压暂降随机经济损失模型;考虑静止同步补偿器和动态电压恢复器的工作特性,基于阻抗矩阵建立多类型补偿设备的数学模型,以补偿设备安装运行总成本及电压暂降经济损失最小为目标函数,应用人工免疫算法对多类型补偿设备进行优化配置,改善配电网电压质量和补偿设备配置的经济性,具有一定的工程实用价值。
Description
技术领域
本发明实施例涉及电能质量控制与分析技术领域,尤其涉及一种计及电压暂降经济损失不确定性的补偿设备优化配置方法。
背景技术
随着经济发展和高精密设备的广泛应用,用户对电压质量的要求越来越高。电压暂降会造成敏感设备的不正常运行,从而给工业生产造成巨大的经济损失。因此,不论从保证电力***的安全、经济运行角度,还是从保证设备安全的角度,迫切需要对电力***中的电压暂降加以治理,及时应对***内电压暂降的发生,分析故障发生时***内的电压暂降分布与各节点的期望电压暂降特征,进而采取必要的防治补偿措施。
为了有效应对电压暂降发生与传播对敏感电力用户影响,提高城市电网的安全可靠性,需要开展综合防治补偿设备的配置研究。考虑到电压暂降防治设备的类型和补偿能力不同与电压暂降对于不同负荷的影响特征差异,对于电压暂降补偿设备的装置的配置需对补偿设备的类型、接入位置、接入容量进行优化。此外,考虑到工程实际,优电压暂降综合防治补偿方案还需兼顾方案的投资成本与电网运行的经济效益。
发明内容
本发明提供一种计及电压暂降经济损失不确定性的补偿设备优化配置方法,考虑电压暂降发生的不确定性,建立多类型补偿设备接入后的电压暂降经济损失的随机模型。同时兼顾***内的电压暂降经济损失成本与补偿设备的配置运行成本,利用人工免疫算法对补偿设备的接入位置、接入容量和接入类型进行优化。该方法计及了电压暂降的随机性,分析了不同补偿设备对于改善电压暂降的工作特性,以一定时间周期内经济性最优为目标进行配置,更符合工程实际。
为实现上述目的,本发明提供以下的技术方案:
一种计及电压暂降经济损失不确定性的补偿设备优化配置方法,所述方法包括:
S10、根据电力***阻抗矩阵建立静止同步补偿器和动态电压恢复器的补偿模型,提出静止同步补偿器和动态电压恢复器在电力***不同位置接入时的补偿能力计算方法;
S20、结合不同故障类型下的补偿设备模型,建立电力***内任意位置发生故障时的电压暂降解析式;
S30、根据电压暂降解析式建立关于故障位置的电压暂降幅值分布的概率密度函数,以故障位置在线路上的概率分布反映电压暂降在不同暂降幅值的概率分布,计算不同电压区间的期望暂降频次及经济损失;
S40、以补偿设备的安装和运行成本及电压暂降经济损失最小为目标函数,以补偿设备容量限制为约束,建立优化配置模型,运用人工免疫算法进行优化求解;
S50、生成初始种群,进行种群免疫操作,包括选择、克隆、变异和克隆抑制,刷新种群,直至满足收敛条件,得到优化配置方案。
进一步地,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S10包括:
S101、将静止同步补偿器(Static Synchronous CoMP1908580ensator,STATCOM)视为电力***中的一个可控电流源,基于阻抗矩阵的补偿模型可表示为:
式中,N为***内的节点数,Ish为STATCOM发出的补偿电流,Ik为故障电流;
(a)三相接地短路故障
根据式(1),STATCOM并联点(设为节点t)和故障点(设为节点k)的暂降电压可表示为:
并假设:
计算矩阵M中各元素,并假设:
可得到STATCOM的补偿电路Ish的计算式如下:
因此,若Vs≥kUt fcosβ,则通过STATCOM进行电压补偿,节点t的电压暂降可补偿至稳态电压水平Ut pre,其补偿电流可按式(7)计算得到;反之,若Vs<kUt fcosβ,则经补偿后的节点t电压幅值可表示为:
计算中,若补偿电流Ish达到补偿电流阈值Irate,则令
Ish=Irate (9)
这种情况下,补偿后节点t的电压暂降幅值式(7)计算得到;
(b)单相接地短路故障
对于不对称短路故障,应用对称分量法进行分析,假设A相为故障相,则有:
式中,Ztt0、Ztt1、Ztt2和Zkk0、Zkk1、Zkk2分别为节点t和节点k的零序、正序和负序自阻抗;Ztk 0、Ztk 1、Ztk 2分别为节点t和节点k之间的零序、正序和负序互阻抗,根据单相接地短路故障的边界条件,可以得到故障相补偿电流Ish a为:
式中,Ua,Ub和Uc分别为未补偿时的各相故障电压,Ish b和Ish c分别为非故障相的补偿电流,由于A相为故障相,故BC相补偿电流Ish b和Ish c均为0,根据式(1)则可得到经STATCOM补偿后的电压暂降幅值;
(c)两相接地短路故障
与单相接地短路故障的计算方法类似,根据式(1)和该故障类型的边界条件,即可得到故障相的补偿电流,在此不再赘述;
(b)两相相间短路故障
与接地短路故障不同,相间短路故障时不存在零序故障分量,将式(1)表示为序分量形式,可表示为:
进而根据两相相间短路故障的边界条件,矩阵M变换为:
则故障相的正序补偿电流Ish 1可表示为:
类似的利用相序变换可得到各相补偿电流与补偿电压;
S102、将动态电压调节器(Dynamic Voltage Regulator,DVR)视为电力***中的一个可控电压源,基于阻抗矩阵的补偿模型可表示为:
式中,UDVR为DVR设备输出的补偿电压,经DVR参与电压暂降补偿后的节点电压可表示为:
式中,Ut为补偿后的节点电压,α为补偿后的节点电压相角,γ为故障时的故障电压相角,θ为正常运行情况下的节点电压相角,DVR参与补偿电压电压暂降所需的所提供的无功功率QDVR可表示为:
本方法仅考虑DVR运行一般情况,即仅提供无功功率进行补偿,则有:
同样的,若Ut f>Ut precosθ,则故障后电压可恢复至正常运行状态下的电压水平,此时,则Ut=Ut pre,根据式(16)可计算得到DVR输出的补偿电压幅值和相角,若DVR参与电压补偿所提供的无功功率达到其容量上限Qrate,则经补偿后的节点电压可表示为:
式中,St为节点t所接负荷的视在功率。
进一步地,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S20包括:
S201、假设故障点k位于线路m-n上的λ处,λ为故障位置在线路上的归一化距离,可表示为:
若故障点k不在线路首末节点,则将其视为虚拟节点,故障点k和其他任一节点t的互阻抗及其自阻抗可分别表示为
式中,zmn 012为线路m-n之间的零序、正序和负序线路阻抗;
S202、建立基于阻抗矩阵的不同故障类型下电压暂降方程,电压方程可表示为(假设A相为故障相)
式中,Ut f和Ut pre分别表示节点t的故障时的电压暂降幅值和稳态电压幅值;SLGF、LLF、DLGF和TPF分别表示单相接地故障、相间短路故障、相间接地短路故障和三相短路故障;
S203、根据步骤S10所述的补偿设备的补偿能力计算方法,推导节点经电压补偿后的电压暂降幅值,可表示为:
式中,UCOM为补偿设备的输出电压,根据补偿位置与补偿设备容量计算可得;Ut为故障时经补偿设备补偿后的电压水平。
进一步地,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S30包括:
S301、根据上式,不同故障类型下电压暂降幅值是关于故障位置的函数,电压暂降幅值分布与线路上不同位置的故障率有关,电压暂降的概率密度函数可表示为:
式中,g(λ)为故障位置在线路的概率分布函数,将归一化故障位置用故障电压表示,反解电压方程,即可得到故障位置关于故障电压的微分,因此,当故障发生在线路j时,故障点位于线路j上区间[λlow,λup]的概率可用节点t的电压暂降幅值在区间[Vi-1,Vi]上的概率表示,即
Pj,t(Vi-1≤|V|≤Vi)=Pj,k(λlow≤|λ|≤λup) (29)
假设Nj,t表示统计的线路j发生故障时,节点t的电压暂降频次,则处于区间[Vi-1,Vi]的暂降电压频次可表示为
S302、建立不同电压区间的电压暂降随机经济损失模型,节点t发生幅值为[Vi-1,Vi]的电压暂降时造成的经济损失可表示为
式中K1 i表示发生幅值为[Vi-1,Vi]的电压暂降时单位负荷造成的经济损失,则电力***内的电压暂降总经济损失可表示为
式中,NS表示***内的敏感负荷节点数,Nv表示所设置的电压暂降幅值区间数,L为***内的线路总数。
进一步地,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S40包括:
在一定的计算周期,计算配置方案的总成本,包括补偿设备的运行配置成本与电压暂降经济损失成本,优化配置模型的目标函数可表示为
式中,K21和K22分别为补偿设备的单位容量补偿成本和配置成本;NCOM为***中配置的补偿设备总数,优化配置模型的约束条件即为STATCOM的补偿电流约束和DVR的补偿容量约束,即
进一步地,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S50包括:
S501、可随机产生NP个候选配置方案作为初始抗体,以总成本为抗体亲和度算子进行个体评价;
S502、进行免疫进化操作,包括免疫选择、克隆、变异、抑制等;
(a)免疫选择:根据抗体亲和度和抗体浓度计算,选择优质抗体,使其活化;
(b)克隆:对活化抗体进行克隆复制,得到若干抗体副本;
(c)变异:对克隆得到的副本进行变异操作,使其发生亲和度突变;
(d)克隆抑制:对变异结果进行再选择,抑制亲和度低的抗体,保留亲和度高的变异结果;
S503、刷新种群,以新抗体代替种群中亲和度较低的抗体,形成新种群,判断是否满足收敛条件,若是,则收敛得到最优配置方案,否则转步骤S502。
本发明实施例提供的一种计及电压暂降经济损失不确定性的补偿设备优化配置方法,能够更准确的分析***的电压暂降造成的经济损失不确定性,相较于传统方法更符合实际情况,兼顾了STATCOM和DVR参与电压暂降补偿时的工作特性差异,提出了分布式补偿和集中式补偿系协同参与电压暂降治理的优化配置方法,充分发挥了不同设备的工作特性,优化了电压暂降改善能力和补偿方案的经济性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1是本发明实施例提供的一种计及电压暂降经济损失不确定性的补偿设备优化配置方法的流程示意图;
图2是本发明实施例中采用免疫算法求解优化配置方案的流程图。
图3是本发明实施例中采用免疫算法求解优化配置方案的收敛曲线。
图4是本发明实施例中采用最优配置方案后各节点电压暂降经济损失对比图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
请参阅附图1,为本发明实施例一提供的一种计及电压暂降经济损失不确定性的补偿设备优化配置方法的流程示意图,该方法具体包括如下步骤:
S10、根据电力***阻抗矩阵建立静止同步补偿器和动态电压恢复器的补偿模型,提出静止同步补偿器和动态电压恢复器在电力***不同位置接入时的补偿能力计算方法;
S20、结合不同故障类型下的补偿设备模型,建立电力***内任意位置发生故障时的电压暂降解析式;
S30、根据电压暂降解析式建立关于故障位置的电压暂降幅值分布的概率密度函数,以故障位置在线路上的概率分布反映电压暂降在不同暂降幅值的概率分布,计算不同电压区间的期望暂降频次及经济损失;
S40、以补偿设备的安装和运行成本及电压暂降经济损失最小为目标函数,以补偿设备容量限制为约束,建立优化配置模型,运用人工免疫算法进行优化求解;
S50、生成初始种群,进行种群免疫操作,包括选择、克隆、变异和克隆抑制,刷新种群,直至满足收敛条件,得到优化配置方案。
优选的,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S10包括:
S101、将静止同步补偿器(Static Synchronous CoMP1908580ensator,STATCOM)视为电力***中的一个可控电流源,基于阻抗矩阵的补偿模型可表示为:
式中,N为***内的节点数,Ish为STATCOM发出的补偿电流,Ik为故障电流;
(a)三相接地短路故障
根据式(1),STATCOM并联点(设为节点t)和故障点(设为节点k)的暂降电压可表示为:
并假设:
计算矩阵M中各元素,并假设:
可得到STATCOM的补偿电路Ish的计算式如下:
因此,若Vs≥kUt fcosβ,则通过STATCOM进行电压补偿,节点t的电压暂降可补偿至稳态电压水平Ut pre,其补偿电流可按式(7)计算得到;反之,若Vs<kUt fcosβ,则经补偿后的节点t电压幅值可表示为:
计算中,若补偿电流Ish达到补偿电流阈值Irate,则令
Ish=Irate (9)
这种情况下,补偿后节点t的电压暂降幅值式(7)计算得到;
(b)单相接地短路故障
对于不对称短路故障,应用对称分量法进行分析,假设A相为故障相,则有:
式中,Ztt0、Ztt1、Ztt2和Zkk0、Zkk1、Zkk2分别为节点t和节点k的零序、正序和负序自阻抗;Ztk 0、Ztk 1、Ztk 2分别为节点t和节点k之间的零序、正序和负序互阻抗,根据单相接地短路故障的边界条件,可以得到故障相补偿电流Ish a为:
式中,Ua,Ub和Uc分别为未补偿时的各相故障电压,Ish b和Ish c分别为非故障相的补偿电流,由于A相为故障相,故BC相补偿电流Ish b和Ish c均为0,根据式(1)则可得到经STATCOM补偿后的电压暂降幅值;
(c)两相接地短路故障
与单相接地短路故障的计算方法类似,根据式(1)和该故障类型的边界条件,即可得到故障相的补偿电流,在此不再赘述;
(b)两相相间短路故障
与接地短路故障不同,相间短路故障时不存在零序故障分量,将式(1)表示为序分量形式,可表示为:
进而根据两相相间短路故障的边界条件,矩阵M变换为:
则故障相的正序补偿电流Ish 1可表示为:
类似的利用相序变换可得到各相补偿电流与补偿电压;
S102、将动态电压调节器(Dynamic Voltage Regulator,DVR)视为电力***中的一个可控电压源,基于阻抗矩阵的补偿模型可表示为:
式中,UDVR为DVR设备输出的补偿电压,经DVR参与电压暂降补偿后的节点电压可表示为:
式中,Ut为补偿后的节点电压,α为补偿后的节点电压相角,γ为故障时的故障电压相角,θ为正常运行情况下的节点电压相角,DVR参与补偿电压电压暂降所需的所提供的无功功率QDVR可表示为:
本方法仅考虑DVR运行一般情况,即仅提供无功功率进行补偿,则有:
同样的,若Ut f>Ut precosθ,则故障后电压可恢复至正常运行状态下的电压水平,此时,则Ut=Ut pre,根据式(16)可计算得到DVR输出的补偿电压幅值和相角,若DVR参与电压补偿所提供的无功功率达到其容量上限Qrate,则经补偿后的节点电压可表示为:
式中,St为节点t所接负荷的视在功率。
优选的,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S20包括:
S201、假设故障点k位于线路m-n上的λ处,λ为故障位置在线路上的归一化距离,可表示为:
若故障点k不在线路首末节点,则将其视为虚拟节点,故障点k和其他任一节点t的互阻抗及其自阻抗可分别表示为
式中,zmn 012为线路m-n之间的零序、正序和负序线路阻抗;
S202、建立基于阻抗矩阵的不同故障类型下电压暂降方程,电压方程可表示为(假设A相为故障相)
式中,Ut f和Ut pre分别表示节点t的故障时的电压暂降幅值和稳态电压幅值;SLGF、LLF、DLGF和TPF分别表示单相接地故障、相间短路故障、相间接地短路故障和三相短路故障;
S203、根据步骤S10所述的补偿设备的补偿能力计算方法,推导节点经电压补偿后的电压暂降幅值,可表示为:
式中,UCOM为补偿设备的输出电压,根据补偿位置与补偿设备容量计算可得;Ut为故障时经补偿设备补偿后的电压水平。
优选的,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S30包括:
S301、根据上式,不同故障类型下电压暂降幅值是关于故障位置的函数,电压暂降幅值分布与线路上不同位置的故障率有关,电压暂降的概率密度函数可表示为:
式中,g(λ)为故障位置在线路的概率分布函数,将归一化故障位置用故障电压表示,反解电压方程,即可得到故障位置关于故障电压的微分,因此,当故障发生在线路j时,故障点位于线路j上区间[λlow,λup]的概率可用节点t的电压暂降幅值在区间[Vi-1,Vi]上的概率表示,即
Pj,t(Vi-1≤|V|≤Vi)=Pj,k(λlow≤|λ|≤λup) (29)
假设Nj,t表示统计的线路j发生故障时,节点t的电压暂降频次,则处于区间[Vi-1,Vi]的暂降电压频次可表示为
S302、建立不同电压区间的电压暂降随机经济损失模型,节点t发生幅值为[Vi-1,Vi]的电压暂降时造成的经济损失可表示为
式中K1 i表示发生幅值为[Vi-1,Vi]的电压暂降时单位负荷造成的经济损失,则电力***内的电压暂降总经济损失可表示为
式中,NS表示***内的敏感负荷节点数,Nv表示所设置的电压暂降幅值区间数,L为***内的线路总数。
优选的,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S40包括:
在一定的计算周期,计算配置方案的总成本,包括补偿设备的运行配置成本与电压暂降经济损失成本,优化配置模型的目标函数可表示为
式中,K21和K22分别为补偿设备的单位容量补偿成本和配置成本;NCOM为***中配置的补偿设备总数,优化配置模型的约束条件即为STATCOM的补偿电流约束和DVR的补偿容量约束,即
优选的,所述计及电压暂降经济损失不确定性的补偿设备优化配置方法中,所述步骤S50包括:
S501、可随机产生NP个候选配置方案作为初始抗体,以总成本为抗体亲和度算子进行个体评价;
S502、进行免疫进化操作,包括免疫选择、克隆、变异、抑制等;
(a)免疫选择:根据抗体亲和度和抗体浓度计算,选择优质抗体,使其活化;
(b)克隆:对活化抗体进行克隆复制,得到若干抗体副本;
(c)变异:对克隆得到的副本进行变异操作,使其发生亲和度突变;
(d)克隆抑制:对变异结果进行再选择,抑制亲和度低的抗体,保留亲和度高的变异结果;
S503、刷新种群,以新抗体代替种群中亲和度较低的抗体,形成新种群,判断是否满足收敛条件,若是,则收敛得到最优配置方案,否则转步骤S502。
利用免疫算法求解多类型补偿设备优化配置方案的流程如图2所示。经计算,免疫算法在第78次迭代时收敛,收敛曲线如图3所示,最优配置方案如表1所示。
表1最优配置方案
接入节点 | 接入类型 | 接入有功功率 |
6 | DVR | 100kVA |
12 | STATCOM | 18MVA |
31 | DVR | 160kVA |
28 | STATCOM | 12MVA |
图4所示为优化前后的节点电压暂降经济损失对比。
本发明实施例提供的一种计及电压暂降经济损失不确定性的补偿设备优化配置方法,利用STATCOM和DVR进行集中式和分布式补偿相结合的方式进行电压暂降综合治理,基于电压暂降分布的概率密度函数预估不同幅值区间的期望暂降频次,建立电压暂降经济损失的随机模型;综合考虑配置成本和电压暂降经济损失成本,利用免疫算法对一定时间周期内的分布式电源配置方案进行优化,改善配置方案的经济性和可靠性。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种计及电压暂降经济损失不确定性的补偿设备优化配置方法,其特征在于,所述方法包括:
S10、根据电力***阻抗矩阵建立静止同步补偿器和动态电压恢复器的补偿模型,提出静止同步补偿器和动态电压恢复器在电力***不同位置接入时的补偿能力计算方法;
S20、结合不同故障类型下的补偿设备模型,建立电力***内任意位置发生故障时的电压暂降解析式;
S30、根据电压暂降解析式建立关于故障位置的电压暂降幅值分布的概率密度函数,以故障位置在线路上的概率分布反映电压暂降在不同暂降幅值的概率分布,计算不同电压区间的期望暂降频次及经济损失;
S40、以补偿设备的安装和运行成本及电压暂降经济损失最小为目标函数,以补偿设备容量限制为约束,建立优化配置模型,运用人工免疫算法进行优化求解;
S50、生成初始种群,进行种群免疫操作,包括选择、克隆、变异和克隆抑制,刷新种群,直至满足收敛条件,得到优化配置方案;
所述步骤S10包括:
S101、将静止同步补偿器视为电力***中的一个可控电流源,基于阻抗矩阵的补偿模型可表示为:
式中,N为***内的节点数,Ish为STATCOM发出的补偿电流,Ik为故障电流;
(a)三相接地短路故障
根据式(1),STATCOM并联点设为节点t,故障点设为节点k的暂降电压可表示为:
并假设:
计算矩阵M中各元素,并假设:
可得到STATCOM的补偿电路Ish的计算式如下:
因此,若Vs≥kUt fcosβ,则通过STATCOM进行电压补偿,节点t的电压暂降可补偿至稳态电压水平Ut pre,其补偿电流可按式(7)计算得到;反之,若Vs<kUt fcosβ,则经补偿后的节点t电压幅值可表示为:
计算中,若补偿电流Ish达到补偿电流阈值Irate,则令
Ish=Irate (9)
这种情况下,补偿后节点t的电压暂降幅值式(7)计算得到;
(b)单相接地短路故障
对于不对称短路故障,应用对称分量法进行分析,假设A相为故障相,则有:
式中,和分别为节点t和节点k的零序、正序和负序自阻抗;Ztk 0、Ztk 1、Ztk 2分别为节点t和节点k之间的零序、正序和负序互阻抗,根据单相接地短路故障的边界条件,可以得到故障相补偿电流Ish a为:
式中,Ua,Ub和Uc分别为未补偿时的各相故障电压,Ish b和Ish c分别为非故障相的补偿电流,由于A相为故障相,故BC相补偿电流Ish b和Ish c均为0,根据式(1)则可得到经STATCOM补偿后的电压暂降幅值;
(c)两相接地短路故障
根据式(1)和该故障类型的边界条件,即可得到故障相的补偿电流;
(d)两相相间短路故障
与接地短路故障不同,相间短路故障时不存在零序故障分量,将式(1)表示为序分量形式,可表示为:
进而根据两相相间短路故障的边界条件,矩阵M变换为:
则故障相的正序补偿电流Ish 1可表示为:
类似的利用相序变换可得到各相补偿电流与补偿电压;
S102、将动态电压调节器视为电力***中的一个可控电压源,基于阻抗矩阵的补偿模型可表示为:
式中,UDVR为DVR设备输出的补偿电压,经DVR参与电压暂降补偿后的节点电压可表示为:
式中,Ut为补偿后的节点电压,α为补偿后的节点电压相角,γ为故障时的故障电压相角,θ为正常运行情况下的节点电压相角,DVR参与补偿电压暂降所需的所提供的无功功率QDVR可表示为:
本方法仅考虑DVR运行一般情况,即仅提供无功功率进行补偿,则有:
同样的,若Ut f>Ut precosθ,则故障后电压可恢复至正常运行状态下的电压水平,此时,则Ut=Ut pre,根据式(16)可计算得到DVR输出的补偿电压幅值和相角,若DVR参与电压补偿所提供的无功功率达到其容量上限Qrate,则经补偿后的节点电压可表示为:
式中,St为节点t所接负荷的视在功率。
2.根据权利要求1所述的计及电压暂降经济损失不确定性的补偿设备优化配置方法,其特征在于,所述步骤S20包括:
S201、假设故障点k位于线路m-n上的λ处,λ为故障位置在线路上的归一化距离,可表示为:
若故障点k不在线路首末节点,则将其视为虚拟节点,故障点k和其他任一节点t的互阻抗及其自阻抗可分别表示为
式中,zmn 012为线路m-n之间的零序、正序和负序线路阻抗;
S202、建立基于阻抗矩阵的不同故障类型下电压暂降方程,电压方程可表示为假设A相为故障相:
式中,Ut f和Ut pre分别表示节点t的故障时的电压暂降幅值和稳态电压幅值;SLGF、LLF、DLGF和TPF分别表示单相接地故障、相间短路故障、相间接地短路故障和三相短路故障;
S203、根据步骤S10所述的补偿设备的补偿能力计算方法,推导节点经电压补偿后的电压暂降幅值,可表示为:
式中,UCOM为补偿设备的输出电压,根据补偿位置与补偿设备容量计算可得;Ut为故障时经补偿设备补偿后的电压水平。
3.根据权利要求2所述的计及电压暂降经济损失不确定性的补偿设备优化配置方法,其特征在于,所述步骤S30包括:
S301、根据上式,不同故障类型下电压暂降幅值是关于故障位置的函数,电压暂降幅值分布与线路上不同位置的故障率有关,电压暂降的概率密度函数可表示为:
式中,g(λ)为故障位置在线路的概率分布函数,将归一化故障位置用故障电压表示,反解电压方程,即可得到故障位置关于故障电压的微分,因此,当故障发生在线路j时,故障点位于线路j上区间[λlow,λup]的概率可用节点t的电压暂降幅值在区间[Vi-1,Vi]上的概率表示,即
Pj,t(Vi-1≤|V|≤Vi)=Pj,k(λlow≤|λ|≤λup) (29)
假设Nj,t表示统计的线路j发生故障时,节点t的电压暂降频次,则处于区间[Vi-1,Vi]的暂降电压频次可表示为
S302、建立不同电压区间的电压暂降随机经济损失模型,节点t发生幅值为[Vi-1,Vi]的电压暂降时造成的经济损失可表示为
式中K1 i表示发生幅值为[Vi-1,Vi]的电压暂降时单位负荷造成的经济损失,则电力***内的电压暂降总经济损失可表示为
式中,NS表示***内的敏感负荷节点数,Nv表示所设置的电压暂降幅值区间数,L为***内的线路总数。
5.根据权利要求4所述的计及电压暂降经济损失不确定性的补偿设备优化配置方法,其特征在于,所述步骤S50包括:
S501、可随机产生NP个候选配置方案作为初始抗体,以总成本为抗体亲和度算子进行个体评价;
S502、进行免疫进化操作,包括免疫选择、克隆、变异、抑制;
(a)免疫选择:根据抗体亲和度和抗体浓度计算,选择优质抗体,使其活化;
(b)克隆:对活化抗体进行克隆复制,得到若干抗体副本;
(c)变异:对克隆得到的副本进行变异操作,使其发生亲和度突变;
(d)克隆抑制:对变异结果进行再选择,抑制亲和度低的抗体,保留亲和度高的变异结果;
S503、刷新种群,以新抗体代替种群中亲和度较低的抗体,形成新种群,判断是否满足收敛条件,若是,则收敛得到最优配置方案,否则转步骤S502。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910403118.1A CN109995050B (zh) | 2019-05-15 | 2019-05-15 | 计及电压暂降经济损失不确定性的补偿设备优化配置方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910403118.1A CN109995050B (zh) | 2019-05-15 | 2019-05-15 | 计及电压暂降经济损失不确定性的补偿设备优化配置方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109995050A CN109995050A (zh) | 2019-07-09 |
CN109995050B true CN109995050B (zh) | 2020-06-26 |
Family
ID=67136490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910403118.1A Active CN109995050B (zh) | 2019-05-15 | 2019-05-15 | 计及电压暂降经济损失不确定性的补偿设备优化配置方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109995050B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110535145B (zh) * | 2019-08-09 | 2020-07-14 | 广东电网有限责任公司 | 对称故障下含statcom电网电压暂降水平计算方法 |
CN110661294B (zh) * | 2019-10-11 | 2020-07-14 | 四川大学 | 基于自适应核密度估计的含有dfig的区间暂降频次估计方法 |
CN110687344B (zh) * | 2019-10-24 | 2021-11-12 | 南京南瑞继保电气有限公司 | 单相电压暂降检测方法及装置、电压恢复器、设备及介质 |
CN111308977B (zh) * | 2020-03-05 | 2021-03-23 | 四川大学 | 暂降及短时中断严重度评估与治理的工业过程建模方法 |
CN111641212B (zh) * | 2020-06-12 | 2022-01-25 | 国网重庆市电力公司电力科学研究院 | 一种储能优化方法、装置及可读存储介质 |
CN112488544A (zh) * | 2020-12-07 | 2021-03-12 | 深圳供电局有限公司 | 电压暂降补偿设备配置方法、装置和计算机设备 |
CN113363975B (zh) * | 2021-06-18 | 2022-09-13 | 国网四川省电力公司技能培训中心 | 一种考虑过程免疫时间的共享电压暂降补偿***及方法 |
CN117559402B (zh) * | 2023-11-13 | 2024-07-16 | 国网宁夏电力有限公司中卫供电公司 | 一种电网抗电压暂降性能评估方法、介质及*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109523165A (zh) * | 2018-11-19 | 2019-03-26 | 东南大学 | 一种基于有限信息的电压暂降状态估计方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101882237A (zh) * | 2010-05-14 | 2010-11-10 | 长沙理工大学 | 一种改进的免疫-粒子群优化算法 |
CN103795068A (zh) * | 2014-03-05 | 2014-05-14 | 广东电网公司电力调度控制中心 | 一种高压配网动态无功补偿装置容量的优化配置方法 |
CN104569702A (zh) * | 2015-02-13 | 2015-04-29 | 国家电网公司 | 一种统一电能质量控制器的测试***及其测试方法 |
CN105958478B (zh) * | 2016-05-23 | 2019-03-19 | 国网江苏省电力公司电力科学研究院 | 一种用于优质电力园区中电压暂降的综合治理方法 |
CN108898305A (zh) * | 2018-06-26 | 2018-11-27 | 国网山东省电力公司德州供电公司 | 主动配电网规划方法及其*** |
-
2019
- 2019-05-15 CN CN201910403118.1A patent/CN109995050B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109523165A (zh) * | 2018-11-19 | 2019-03-26 | 东南大学 | 一种基于有限信息的电压暂降状态估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109995050A (zh) | 2019-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109995050B (zh) | 计及电压暂降经济损失不确定性的补偿设备优化配置方法 | |
US7979239B2 (en) | System and method for monitoring and managing electrical power transmission and distribution networks | |
Mohiuddin et al. | Optimal distributed control of AC microgrids with coordinated voltage regulation and reactive power sharing | |
CN111327047B (zh) | 一种确定级联式电网多交流断面输电能力的方法及*** | |
Ramos et al. | Placement and sizing of utility-size battery energy storage systems to improve the stability of weak grids | |
Meng et al. | Offshore transmission network planning for wind integration considering AC and DC transmission options | |
CN105741016B (zh) | 一种用于中期电网规划的静态可靠性概率指标获取方法 | |
CN113239512B (zh) | 一种考虑韧性的交直流配电网规划方案筛选方法及*** | |
Zhao et al. | Optimal transmission switching to eliminate voltage violations during light-load periods using decomposition approach | |
Saadatmand et al. | Optimal coordinated tuning of power system stabilizers and wide‐area measurement‐based fractional‐order PID controller of large‐scale PV farms for LFO damping in smart grids | |
Sanni et al. | Strength assessment of electric power systems containing inverter-based distributed generation | |
Rupolo et al. | A new parallel and decomposition approach to solve the medium-and low-voltage planning of large-scale power distribution systems | |
CN104767207A (zh) | 基于发电机无功储备灵敏度的电压稳定预防控制方法 | |
CN112186764B (zh) | 一种配电网设备的接入优化方法、装置及电子设备 | |
Liu et al. | Sensitivity Region-Based Optimization for Maximizing Renewable Generation Hosting Capacity of an Islanded Microgrid | |
CN110994621B (zh) | 基于改进直流潮流算法的电网有功潮流调整方法及装置 | |
Van Dai et al. | Improving Power System Stability with Gramian Matrix‐Based Optimal Setting of a Single Series FACTS Device: Feasibility Study in Vietnamese Power System | |
Valencia | Expansion planning of joint medium-and low-voltage three-phase distribution networks considering the optimal integration of distributed energy resources | |
Dehghani-Arani et al. | Introduction a multi-objective function in unbalanced and unsymmetrical distribution networks for optimal placement and sizing of distributed generation units using NSGA-II | |
Pan et al. | A flexible black‐start network partitioning strategy considering subsystem recovery time balance | |
Behbahani et al. | Reconfiguration of distribution network for improving power quality indexes with flexible lexicography method | |
Margossian | Iterative state estimation with weight tuning and pseudo-measurement generation | |
Castillo et al. | Optimal location and size for various renewable distributed generators in distribution networks | |
Prasad et al. | A synchrophasor measurements based adaptive underfrequency load shedding scheme | |
Nazir et al. | ADMM based distributed load flow for radial power networks |
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 |