CN104362642B - 提高交直流电网中长期电压稳定的动态无功备用优化方法 - Google Patents
提高交直流电网中长期电压稳定的动态无功备用优化方法 Download PDFInfo
- Publication number
- CN104362642B CN104362642B CN201410584184.0A CN201410584184A CN104362642B CN 104362642 B CN104362642 B CN 104362642B CN 201410584184 A CN201410584184 A CN 201410584184A CN 104362642 B CN104362642 B CN 104362642B
- Authority
- CN
- China
- Prior art keywords
- dynamic
- term
- compensation equipment
- voltage
- passive compensation
- 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
- 230000007774 longterm Effects 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000011105 stabilization Methods 0.000 title claims abstract description 36
- 230000006641 stabilisation Effects 0.000 title claims abstract description 33
- 238000005457 optimization Methods 0.000 claims abstract description 28
- 230000035945 sensitivity Effects 0.000 claims description 15
- 239000003990 capacitor Substances 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 240000002853 Nelumbo nucifera Species 0.000 claims description 7
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 7
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims description 6
- 230000002068 genetic effect Effects 0.000 claims description 6
- 239000005364 simax Substances 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 5
- XXXSILNSXNPGKG-ZHACJKMWSA-N Crotoxyphos Chemical compound COP(=O)(OC)O\C(C)=C\C(=O)OC(C)C1=CC=CC=C1 XXXSILNSXNPGKG-ZHACJKMWSA-N 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 abstract description 6
- 238000010276 construction Methods 0.000 abstract 1
- 230000010354 integration Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 5
- 230000001052 transient effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 206010068052 Mosaicism Diseases 0.000 description 1
- 208000032370 Secondary transmission Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- H02J1/00—Circuit arrangements for dc mains or dc distribution networks
-
- 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)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Control Of Electrical Variables (AREA)
Abstract
本发明提供一种提高交直流电网中长期电压稳定的动态无功备用优化方法,包括以下步骤:确定影响交直流电网中长期电压稳定的关键故障集合;调整动态无功补偿设备的无功出力,并计算动态无功补偿设备的灵敏度;对m个动态无功补偿设备进行排序,并计算动态无功补偿设备的权重系数;计算动态无功补偿设备备用容量,建立动态无功备用优化模型,并求解该动态无功备用优化模型。本发明为提高多直流落点电网中长期电压稳定水平提供了辅助决策支持,对提高大规模交直流电网中长期电压稳定裕度,建立送、受端之间畅通的电力传输通道,提升交直流输电通道输送能力,改善电网运行的经济性和电能质量,均具有重大意义。
Description
技术领域
本发明属于电力***技术领域,具体涉及一种提高交直流电网中长期电压稳定的动态无功备用优化方法。
背景技术
电压稳定问题被国内外学者重视以来,已经发展为多个研究分支,对电压稳定也有了明确合理的定义,经过多年研究努力,电力学者们已经在电压稳定问题的某些领域中取得了***的成果,如对静态电压稳定分析、动态电压稳定中的小干扰电压稳定及暂态电压稳定分析中,都已形成了一套较为完善的研究理论和分析方法,在电力***调度运行及监测控制等方面都发挥着不可替代的作用。然而目前国内对中长期电压稳定问题的研究尚不够深入,没有形成较为统一的认识,人们对中长期电压失稳机理及过程不能进行详细严谨的分析,因此,研究中长期电压稳定问题具有非常重要的理论意义。
电力***遭受大扰动后,由于负荷的电压灵敏性可能暂时保持电压稳定,然而电力***中很多影响电压稳定性的元件都存在慢动态动作过程,随着有载调压变压器分接头的换接,以及具有负荷恢复特性的元件功率恢复,经过一个较长的时间过程后,***仍然存在发生电压崩溃的可能,这就是中长期电压稳定问题,中长期电压稳定分析所研究时域范围为几分钟甚至几十分钟。负荷恢复特性对中长期电压失稳有极大影响,具有恢复特性的元件主要有感应电动机和恒温负荷,同时有载调压变压器分接头换接是造成负荷恢复的重要原因。由于上述3种动态元件的响应时间常数长短不一,从而形成快慢动态结合的中长期电压失稳过程。
当前,缺乏有效、快速、适应性强的电压稳定控制方法也是引发大停电事故的重要原因之一。我国虽然未曾发生由电压稳定问题引起的大停电事故,但随着“西电东送,南北互供”电力***联网格局的形成,负荷中心水平不断增长,大容量远距离输电不断增加,我国电力***的电压稳定性问题日益突出,发生电压失稳事故的几率也越来越大。由于电压稳定问题具有隐蔽性和突发性,事故期间难以察觉,一旦发生电压崩溃,在我国当前电网实际情况下,势必造成极其巨大的损失。因此,研究提高中长期电压稳定的动态无功备用优化问题,有效防止电压失稳和电压崩溃事故发生,具有重要的理论价值和实际意义。
发明内容
为了克服上述现有技术的不足,本发明提供一种提高交直流电网中长期电压稳定的动态无功备用优化方法,为提高多直流落点电网中长期电压稳定水平提供了辅助决策支持,对提高大规模交直流电网中长期电压稳定裕度,建立送、受端之间畅通的电力传输通道,提升交直流输电通道输送能力,改善电网运行的经济性和电能质量,均具有重大意义。
为了实现上述发明目的,本发明采取如下技术方案:
本发明提供一种提高交直流电网中长期电压稳定的动态无功备用优化方法,所述方法包括以下步骤:
步骤1:确定影响交直流电网中长期电压稳定的关键故障集合;
步骤2:调整动态无功补偿设备的无功出力,并计算动态无功补偿设备的灵敏度;
步骤3:对m个动态无功补偿设备进行排序,并计算动态无功补偿设备的权重系数;
步骤4:计算动态无功补偿设备备用容量,建立动态无功备用优化模型,并求解该动态无功备用优化模型。
所述步骤1中,对交直流电网进行故障扫描,计算负荷母线i的电压稳定裕度KMVSi,有:
其中,ZLi为负荷母线i处的负荷等值阻抗,ZTi为***戴维南等值阻抗;
选取KMVSi最小值为交直流电网的电压稳定裕度,记为KMVSI,根据交直流电网的电压稳定裕度值确定故障的严重情况,得到关键故障,从而得到关键故障集合。
所述步骤2中,动态无功补偿设备包括发电机、静止无功补偿器和静止同步补偿器。
所述步骤2具体包括以下步骤:
步骤2-1:分别调整各动态无功补偿设备的无功出力,并对关键故障再次进行时域仿真;
步骤2-2:在中长期时间尺度下,针对某故障l,计算动态无功补偿设备j的灵敏度SIl,j;
步骤2-3:在中长期时间尺度下,针对多个故障,计算动态无功补偿设备j的灵敏度SIj。
所述步骤2-2中,针对某故障l,动态无功补偿设备j的灵敏度SIl,j表示为:
其中,Qj0为动态无功补偿设备j的初始无功出力;ΔQj为调整动态无功补偿设备j的无功功率变化量;ΔQRj为调整动态无功补偿设备j的无功备用变化量;kMVSI,l(Qj0+ΔQj)为调整动态无功补偿设备j的无功出力后,在故障Fl下,交直流电网的负荷裕度值;kMVSI,l(Qj0)为调整动态无功补偿设备j的无功出力前,在故障Fl下,交直流电网的负荷裕度值。
所述步骤2-3中,针对多个故障,动态无功补偿设备j的灵敏度SIj表示为:
其中,Nl为关键故障总数。
所述步骤3具体包括以下步骤:
步骤3-1:根据SIj对m个动态无功补偿设备进行排序,SIj最大值表征该动态无功补偿设备对中长期电压稳定的贡献程度最大,贡献程度大的动态无功补偿设备留出更多无功备用量;
步骤3-2:以SIj最大值SImax为基准,归一化处理SIj,计算动态无功补偿设备的权重系数pj,有pj=SIj/|SImax|。
所述步骤4具体包括以下步骤:
步骤4-1:计算动态无功补偿设备的备用容量QRM;
步骤4-2:以提高QRM作为动态无功备用优化目标,建立动态无功备用优化模型;
步骤4-3:采用遗传算法求解该动态无功备用优化模型。
所述步骤4-1中,动态无功补偿设备的备用容量QRM表示为:
其中,Qgjmax为中长期电压稳定中动态无功补偿设备j的无功出力上限,Qgj为动态无功补偿设备j的当前无功出力。
所述步骤4-2中,动态无功备用优化模型的目标函数为:
动态无功备用优化模型的约束条件包括潮流方程约束和变量约束;所述变量约束为控制变量约束和状态变量约束;
(1)潮流方程约束:
在动态无功备用优化模型中,各个节点的有功出力和无功出力都满足以下潮流方程,有:
其中,PGi和QGi分别为电力***中发电机节点的有功出力和无功出力;PLi和QLi分别为负荷节点的有功出力和无功出力;QCi为节点的无功补偿容量;Gir和Bir分别为节点i、r之间的电导和电纳;Vi和Vr分别为节点i、r的电压;δir为节点i、r之间的电压相角差;n为节点总数;Pti(dc)和Qti(dc)分别为直流节点的有功输入和无功输入,分为以下两种情况:
1)节点i在整流侧换流母线上,Pti(dc)和Qti(dc)分别表示为:
其中,kp为换流器的极数;UdR为整流侧直流电压;Id为直流线路电流;KdR为整流侧换流变压器变比;b为每极的6脉波串联桥数;VR为整流侧的交流母线电压幅值;
2)节点i在逆变侧换流母线上,Pti(dc)和Qti(dc)分别表示为:
其中,UdI为逆变侧直流电压;KdI为逆变侧换流变压器变比;VI为逆变侧的交流母线电压幅值;
(2)控制变量约束:
其中,NG、NSVC、NSVG、NC、NT和Ndc分别为发电机节点数、静止无功补偿器节点数、静止同步补偿器节点数、并联电容器节点数、变压器可调分接头数和直流网络节点数;VGi为发电机节点的端电压,VGimin和VGimax分别为VGi的下限值和上限值;VSVCg为静止无功补偿器节点的端电压,VSVCgmin和VSVCgmax分别为VSVCg的下限值和上限值;VSVGh为静止同步补偿器节点的端电压,VSVGhmin和VSVGhmax分别为VSVGh下限值和上限值;QCu为并联电容器组的补偿容量,QCumin和QCumax分别为QCu下限值和上限值;Tk为变压器可调分接头,Tkmin和Tkmax分别为Tk下限值和上限值;Udl、Idm、Pdn和θdr分别为换流器控制电压、控制电流、控制功率以及控制角,Udlmin和Udlmax、Idmmin和Idmmax、Pdnmin和Pdnmax、θdrmin和θdrmax分别表示相应的下限值和上限值;
(3)状态变量约束:
其中,NL为负荷节点数;QGi为发电机节点无功出力,QGimin和QGimax分别为QGi的下限值和上限值;BSVCg为静止无功补偿器电纳,BSVCgmin和BSVCgmax分别为BSVCg的下限值和上限值;ISVGh为静止同步补偿器电流幅值,ISVGhmin和ISVGhmax分别为ISVGh的下限值和上限值;VLp为负荷节点电压幅值,VLpmin和VLpmax分别为VLp的下限值和上限值。
与现有技术相比,本发明的有益效果在于:
1.目前尚无适用于多馈入直流电网特征的提高中长期电压稳定的动态无功备用优化技术,本发明创新性地提出了一种适用于多馈入直流电网特征提高中长期电压稳定的动态无功备用优化方法;
2.与基于静态的传统无功备用优化方法相比,本方法详细考虑了***的动态特性,能够更加准确地确定动态无功补偿设备备用容量,为电网的优化运行提供基础;
3.通过时域仿真分析,可快捷、方便、准确地确定各无功源的参与因子,可应用于大规模电力***的动态无功备用优化,克服了传统电力***动态无功优化的算法只能应用于小***的缺点。
附图说明
图1是本发明实施例中提高交直流电网中长期电压稳定的动态无功备用优化方法流程图;
图2是本发明实施例中采用遗传算法求解动态无功备用优化模型流程图;
图3是本发明实施中3机10节点修正测试交直流***示意图;
图4是本发明实施例中发电机相对功角变化曲线图;
图5是本发明实施例中发电机2和发电机3的励磁电流曲线图;
图6是本发明实施例中节点9和节点10电压变化曲线图;
图7是本发明实施例中优化前后节点3(发电机G3机端)电压曲线图;
图8是本发明实施例中优化前后节点10电压曲线图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明提供一种提高交直流电网中长期电压稳定的动态无功备用优化方法,所述方法包括以下步骤:
步骤1:确定影响交直流电网中长期电压稳定的关键故障集合;
步骤2:调整动态无功补偿设备的无功出力,并计算动态无功补偿设备的灵敏度;
步骤3:对m个动态无功补偿设备进行排序,并计算动态无功补偿设备的权重系数;
步骤4:计算动态无功补偿设备备用容量,建立动态无功备用优化模型,并求解该动态无功备用优化模型。
所述步骤1中,对交直流电网进行故障扫描,计算负荷母线i的电压稳定裕度KMVSi,有:
其中,ZLi为负荷母线i处的负荷等值阻抗,ZTi为***戴维南等值阻抗;
选取KMVSi最小值为交直流电网的电压稳定裕度,记为KMVSI,根据交直流电网的电压稳定裕度值确定故障的严重情况,得到关键故障,从而得到关键故障集合。
所述步骤2中,动态无功补偿设备包括发电机、静止无功补偿器和静止同步补偿器。
所述步骤2具体包括以下步骤:
步骤2-1:分别调整各动态无功补偿设备的无功出力,并对关键故障再次进行时域仿真;
步骤2-2:在中长期时间尺度下,针对某故障l,计算动态无功补偿设备j的灵敏度SIl,j;
步骤2-3:在中长期时间尺度下,针对多个故障,计算动态无功补偿设备j的灵敏度SIj。
所述步骤2-2中,针对某故障l,动态无功补偿设备j的灵敏度SIl,j表示为:
其中,Qj0为动态无功补偿设备j的初始无功出力;ΔQj为调整动态无功补偿设备j的无功功率变化量;ΔQRj为调整动态无功补偿设备j的无功备用变化量;kMVSI,l(Qj0+ΔQj)为调整动态无功补偿设备j的无功出力后,在故障Fl下,交直流电网的负荷裕度值;kMVSI,l(Qj0)为调整动态无功补偿设备j的无功出力前,在故障Fl下,交直流电网的负荷裕度值。
所述步骤2-3中,针对多个故障,动态无功补偿设备j的灵敏度SIj表示为:
其中,Nl为关键故障总数。
所述步骤3具体包括以下步骤:
步骤3-1:根据SIj对m个动态无功补偿设备进行排序,SIj最大值表征该动态无功补偿设备对中长期电压稳定的贡献程度最大,贡献程度大的动态无功补偿设备留出更多无功备用量;
步骤3-2:以SIj最大值SImax为基准,归一化处理SIj,计算动态无功补偿设备的权重系数pj,有pj=SIj/|SImax|。
所述步骤4具体包括以下步骤:
步骤4-1:计算动态无功补偿设备的备用容量QRM;
步骤4-2:以提高QRM作为动态无功备用优化目标,建立动态无功备用优化模型;
步骤4-3:采用遗传算法求解该动态无功备用优化模型。
所述步骤4-1中,动态无功补偿设备的备用容量QRM表示为:
其中,Qgjmax为中长期电压稳定中动态无功补偿设备j的无功出力上限,Qgj为动态无功补偿设备j的当前无功出力。
所述步骤4-2中,动态无功备用优化模型的目标函数为:
动态无功备用优化模型的约束条件包括潮流方程约束和变量约束;所述变量约束为控制变量约束和状态变量约束;
(1)潮流方程约束:
在动态无功备用优化模型中,各个节点的有功出力和无功出力都满足以下潮流方程,有:
其中,PGi和QGi分别为电力***中发电机节点的有功出力和无功出力;PLi和QLi分别为负荷节点的有功出力和无功出力;QCi为节点的无功补偿容量;Gir和Bir分别为节点i、r之间的电导和电纳;Vi和Vr分别为节点i、r的电压;δir为节点i、r之间的电压相角差;n为节点总数;Pti(dc)和Qti(dc)分别为直流节点的有功输入和无功输入,分为以下两种情况:
1)节点i在整流侧换流母线上,Pti(dc)和Qti(dc)分别表示为:
其中,kp为换流器的极数;UdR为整流侧直流电压;Id为直流线路电流;KdR为整流侧换流变压器变比;b为每极的6脉波串联桥数;VR为整流侧的交流母线电压幅值;
2)节点i在逆变侧换流母线上,Pti(dc)和Qti(dc)分别表示为:
其中,UdI为逆变侧直流电压;KdI为逆变侧换流变压器变比;VI为逆变侧的交流母线电压幅值;
(2)控制变量约束:
其中,NG、NSVC、NSVG、NC、NT和Ndc分别为发电机节点数、静止无功补偿器节点数、静止同步补偿器节点数、并联电容器节点数、变压器可调分接头数和直流网络节点数;VGi为发电机节点的端电压,VGimin和VGimax分别为VGi的下限值和上限值;VSVCg为静止无功补偿器节点的端电压,VSVCgmin和VSVCgmax分别为VSVCg的下限值和上限值;VSVGh为静止同步补偿器节点的端电压,VSVGhmin和VSVGhmax分别为VSVGh下限值和上限值;QCu为并联电容器组的补偿容量,QCumin和QCumax分别为QCu下限值和上限值;Tk为变压器的变比,Tkmin和Tkmax分别为Tk下限值和上限值;Udl、Idm、Pdn和θdr分别为换流器控制电压、控制电流、控制功率以及控制角,Udlmin和Udlmax、Idmmin和Idmmax、Pdnmin和Pdnmax、θdrmin和θdrmax分别表示相应的下限值和上限值;
(3)状态变量约束:
其中,NL为负荷节点数;QGi为发电机节点无功出力,QGimin和QGimax分别为QGi的下限值和上限值;BSVCg为静止无功补偿器电纳,BSVCgmin和BSVCgmax分别为BSVCg的下限值和上限值;ISVGh为静止同步补偿器电流幅值,ISVGhmin和ISVGhmax分别为ISVGh的下限值和上限值;VLp为负荷节点电压幅值,VLpmin和VLpmax分别为VLp的下限值和上限值。
步骤4-3中,采用遗传算法求解该动态无功备用优化模型;
遗传算法的基本思想是,在某特定环境下的一群个体,由于环境限制,只有适应性强的可生存,而弱者被淘汰,它们适应环境的优良性状会遗传给后代。GA应用于无功备用优化问题时可以理解为:电力***下的一组初始潮流解,受各种约束条件约束,通过目标函数评价其优劣,评价值低的被抛弃,只有评价值高的有机会将其特征迭代至下一轮解,最后趋向最优。
具体过程如下:
(1)首先,根据下式随机产生第一代母体,有:
Xi=INT(RND(Ximax-Xiimn))+Ximin (11)
其中,RND为随机数,且0<RND<1;INT(*)为取整;
Xi若为VGi,则Ximax、Ximax分别表示发电机节点的端电压上下限;
Xi若为VSVCg,则Ximax、Ximax分别表示静止无功补偿器节点的端电压上下限;
Xi若为VSVGh,则Ximax、Ximax分别表示静止同步补偿器节点的端电压上下限;
Xi若为QCu,则Ximax、Ximax分别表示并联电容器组的补偿容量上下限;
Xi若为Tk,则Ximax、Ximax分别表示变压器的变比上下限。
式(11)使变量的约束方程可转化为整数变量的约束方程。如1+5×0.025%的变压器其YT的取值范围为1~11。编码采用二进制数,每五位顺序表示YVgi、YVsvcg、YVsvgh、YQcu、YTk的值:
H=[…,b5i-4,…,b5i,…,b5g-4,…b5g,…,b5h-4,…b5h,…,b5u-4,…,b5u,…,b5k-4,…,b5k,…] (12)
(2)对A中的每个个体根据式(13)进行解码,修改原始潮流数据中对应的值,再开始潮流计算,本发明的潮流计算程序采用的是N-R法;
式中:ΔVGi、ΔVSVCg、ΔVSVGh、ΔQCu、ΔTk对应变量有级调节单元值;
YVgi、YVsvcg、YVsvgh、YQcu、YTk代表控制变量开关位置的整数变量;
YVgi=1,表示第i个发电机节点端电压调到最大;
YVsvcg=1,表示第g个静止无功补偿器节点端电压调到最大;
YVsvgh=1,表示第h个静止同步补偿器节点端电压调到最大;
YQcu=1,表示第j个电容器投入一组电容量;
YTk=1,表示第k个变压器分接头置于变比最大位置;
(3)经过潮流计算,获得了各节点的电压、无功和动态无功备用容量等数据,并将其由大到小排序;
(4)依据适应值大小对各个体进行排序,保留亲和力大的个体组成个体群B,同时对B内的个体进行交叉变异操作,保留操作后整体适应值大的个体,组成个体群C;依据适应值大小对B、C进行排列组成个体群D;
(5)检查迭代结束条件,如果达到则结束,否则转下一步;
(6)随机产生一组新的个体群E,与D共同组成新一代迭代计算个体群F,转步骤(2),重新开始计算。
实施例
如图3所示,针对3机10节点***,500kV母线(Bus6)向负荷区域的两个负荷供电,其中的工业负荷(节点Bus7)通过OLTC变压器与500kV负荷母线连接,而居民负荷与商业负荷(节点Bus10)则通过两台OLTC变压器和一段代表次级输电***的阻抗接在500kV负荷母线。负荷地区有一台1600MVA的等值发电机(节点Bus3),并采用了大量的并联补偿装置,节点8上分别配置了容量为±240Mvar的静止无功补偿器(SVC)和容量为600Mvar的电容器组,该电容器组每组容量为100Mvar,共6组。两台远方的发电机通过4条500kV线路和1回双极直流输电线路向负荷区域输送功率。仿真所采用的主要模型:变压器(Bus9~Bus10)为OLTC变压器,其它分接头保持不变;节点Bus7上的负荷为恒功率模型,其它负荷为恒阻抗模型;发电机2和3(节点Bus2和Bus3)上的发电机有过励磁限制装置,发电机1(节点Bus1)为无穷大发电机。
对此***进行故障扫描,确定威胁***中长期电压稳定的关键故障集合。为了方便地说明TSI指标的有效性,本算例只考察最严重的N-1故障,故障形式为t=0.1s时节点5~节点6间的一回交流联络线在节点6侧发生三相永久性短路故障,故障后0.09s跳开故障线路节点6侧开关,0.1s跳开故障线路节点5侧开关。
图4表示发电机2、发电机3分别与发电机1之间相对功角摇摆曲线。由图4可见,这个扰动引起的初始快速暂态过程会很快消失,表明***是可以保持暂态功角稳定的,后续的中长期过程中功角虽有摇摆,但角度都比较小,表明***也可以保持中长期功角稳定。
图5中过励磁限制器的励磁电流表明了发电机2和发电机3电动势Eq的响应,这个电动势正比于励磁电流。如图所示,扰动之后,发电机2和发电机3的励磁电流会突然上升,如果超过了转子电流限制,就会启动过励磁限制器的反时间机制。扰动后,OLTC变压器的运行对发电机强加了一个非常重的无功需求。这个需求进一步恶化了转子过负荷,直到最终过励磁限制器被激励,致使励磁电流回到其额定值。注意,这个过励限制器是积分型的,以至于Eq被强迫到Eq lim。随后的分接头变换导致暂态励磁电流升高,这个升高的励磁电流很快被过励磁限制器所检测(如图5中发电机2、3励磁电流的最大值点),并给予校正。
图6给出了节点9和节点10的电压,即给负荷供电的OLTC的高压侧母线9电压和负荷侧节点10电压。由图可见,暂态过程中,节点9能够在0.87p.u.处稳定运行。OLTC变压器通过降低变比Tk,设法恢复负荷侧节点10电压。经过30秒的初始时间延迟后,OLTC转换器开始运行,约55秒,经过5次分接头调整之后,母线电压上升到0.915pu,非常接近事故前水平,发电机2和3的励磁电流输出也随之增加以满足***对无功功率的需求(图5)。但是到了347秒时,由于发电机3的过励磁限制装置开始动作,限制了其输出电流,使得该机的无功功率输出也随之下降,导致负荷节点10的电压再次下降。为保证电压,变压器分接头继续动作了18次。在442秒时,发电机2的过励磁限制装置也开始动作,***的无功功率缺口大大增加,导致了电压崩溃的发生。
在用本发明求出无功备用容量优化问题各控制变量后,利用时域仿真验证分析所提方法的有效性。
节点5~节点6间的一回交流联络线在节点6侧发生三相永久性短路故障,故障后0.09s跳开故障线路节点6侧开关,0.1s跳开故障线路节点5侧开关。图7和图8分别为发电机G3机端电压和节点10电压曲线,从图中可以看出,优化后***的中长期电压稳定性比优化前要好,这说明采用本发明提出的优化算法能有效提高电网的中长期电压稳定性。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (9)
1.提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述方法包括以下步骤:
步骤1:确定影响交直流电网中长期电压稳定的关键故障集合;
步骤2:调整动态无功补偿设备的无功出力,并计算动态无功补偿设备的灵敏度;
步骤3:对m个动态无功补偿设备进行排序,并计算动态无功补偿设备的权重系数;
步骤4:计算动态无功补偿设备备用容量,建立动态无功备用优化模型,并求解该动态无功备用优化模型;
所述步骤1中,对交直流电网进行故障扫描,计算负荷母线i的电压稳定裕度KMVSi,有:
其中,ZLi为负荷母线i处的负荷等值阻抗,ZTi为***戴维南等值阻抗;
选取KMVSi最小值为交直流电网的电压稳定裕度,记为KMVSI,根据交直流电网的电压稳定裕度值确定故障的严重情况,得到关键故障,从而得到关键故障集合。
2.根据权利要求1所述的提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述步骤2中,动态无功补偿设备包括发电机、静止无功补偿器和静止同步补偿器。
3.根据权利要求1所述的提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述步骤2具体包括以下步骤:
步骤2-1:分别调整各动态无功补偿设备的无功出力,并对关键故障再次进行时域仿真;
步骤2-2:在中长期时间尺度下,针对某故障l,计算动态无功补偿设备j的灵敏度SIl,j;
步骤2-3:在中长期时间尺度下,针对多个故障,计算动态无功补偿设备j的灵敏度SIj。
4.根据权利要求3所述的提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述步骤2-2中,针对某故障l,动态无功补偿设备j的灵敏度SIl,j表示为:
其中,Qj0为动态无功补偿设备j的初始无功出力;ΔQj为调整动态无功补偿设备j的无功功率变化量;ΔQRj为调整动态无功补偿设备j的无功备用变化量;kMVSI,l(Qj0+ΔQj)为调整动态无功补偿设备j的无功出力后,在故障l下,交直流电网的负荷裕度值;kMVSI,l(Qj0)为调整动态无功补偿设备j的无功出力前,在故障l下,交直流电网的负荷裕度值。
5.根据权利要求3所述的提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述步骤2-3中,针对多个故障,动态无功补偿设备j的灵敏度SIj表示为:
其中,Nl为关键故障总数。
6.根据权利要求3所述的提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述步骤3具体包括以下步骤:
步骤3-1:根据SIj对m个动态无功补偿设备进行排序,SIj最大值表征该动态无功补偿设备对中长期电压稳定的贡献程度最大,贡献程度大的动态无功补偿设备留出更多无功备用量;
步骤3-2:以SIj最大值SImax为基准,归一化处理SIj,计算动态无功补偿设备的权重系数pj,有pj=SIj/|SImax|。
7.根据权利要求1所述的提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述步骤4具体包括以下步骤:
步骤4-1:计算动态无功补偿设备的备用容量QRM;
步骤4-2:以提高QRM作为动态无功备用优化目标,建立动态无功备用优化模型;
步骤4-3:采用遗传算法求解该动态无功备用优化模型。
8.根据权利要求7所述的提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述步骤4-1中,动态无功补偿设备的备用容量QRM表示为:
其中,Qgjmax为中长期电压稳定中动态无功补偿设备j的无功出力上限,Qgj为动态无功补偿设备j的当前无功出力。
9.根据权利要求7所述的提高交直流电网中长期电压稳定的动态无功备用优化方法,其特征在于:所述步骤4-2中,动态无功备用优化模型的目标函数为:
动态无功备用优化模型的约束条件包括潮流方程约束和变量约束;所述变量约束为控制变量约束和状态变量约束;
(1)潮流方程约束:
在动态无功备用优化模型中,各个节点的有功出力和无功出力都满足以下潮流方程,有:
其中,PGi和QGi分别为电力***中发电机节点的有功出力和无功出力;PLi和QLi分别为负荷节点的有功出力和无功出力;QCi为节点的无功补偿容量;Gir和Bir分别为节点i、r之间的电导和电纳;Vi和Vr分别为节点i、r的电压;δir为节点i、r之间的电压相角差;n为节点总数;Pti(dc)和Qti(dc)分别为直流节点的有功输入和无功输入,分为以下两种情况:
1)节点i在整流侧换流母线上,Pti(dc)和Qti(dc)分别表示为:
其中,kp为换流器的极数;UdR为整流侧直流电压;Id为直流线路电流;KdR为整流侧换流变压器变比;b为每极的6脉波串联桥数;VR为整流侧的交流母线电压幅值;
2)节点i在逆变侧换流母线上,Pti(dc)和Qti(dc)分别表示为:
其中,UdI为逆变侧直流电压;KdI为逆变侧换流变压器变比;VI为逆变侧的交流母线电压幅值;
(2)控制变量约束:
其中,NG、NSVC、NSVG、NC、NT和Ndc分别为发电机节点数、静止无功补偿器节点数、静止同步补偿器节点数、并联电容器节点数、变压器可调分接头数和直流网络节点数;VGi为发电机节点的端电压,VGimin和VGimax分别为VGi的下限值和上限值;VSVCg为静止无功补偿器节点的端电压,VSVCgmin和VSVCgmax分别为VSVCg的下限值和上限值;VSVGh为静止同步补偿器节点的端电压,VSVGhmin和VSVGhmax分别为VSVGh下限值和上限值;QCu为并联电容器组的补偿容量,QCumin和QCumax分别为QCu下限值和上限值;Tk为变压器可调分接头,Tkmin和Tkmax分别为Tk下限值和上限值;Udl、Idm、Pdn和θdr分别为换流器控制电压、控制电流、控制功率以及控制角,Udlmin和Udlmax、Idmmin和Idmmax、Pdnmin和Pdnmax、θdrmin和θdrmax分别表示相应的下限值和上限值;
(3)状态变量约束:
其中,NL为负荷节点数;QGi为发电机节点无功出力,QGimin和QGimax分别为QGi的下限值和上限值;BSVCg为静止无功补偿器电纳,BSVCgmin和BSVCgmax分别为BSVCg的下限值和上限值;ISVGh为静止同步补偿器电流幅值,ISVGhmin和ISVGhmax分别为ISVGh的下限值和上限值;VLp为负荷节点电压幅值,VLpmin和VLpmax分别为VLp的下限值和上限值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410584184.0A CN104362642B (zh) | 2014-10-27 | 2014-10-27 | 提高交直流电网中长期电压稳定的动态无功备用优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410584184.0A CN104362642B (zh) | 2014-10-27 | 2014-10-27 | 提高交直流电网中长期电压稳定的动态无功备用优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104362642A CN104362642A (zh) | 2015-02-18 |
CN104362642B true CN104362642B (zh) | 2017-01-11 |
Family
ID=52529880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410584184.0A Active CN104362642B (zh) | 2014-10-27 | 2014-10-27 | 提高交直流电网中长期电压稳定的动态无功备用优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104362642B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105281342B (zh) * | 2015-10-26 | 2017-06-16 | 海南电网有限责任公司 | 优先保证动态无功裕度的光伏发电站无功分配方法 |
CN106099908B (zh) * | 2016-05-20 | 2022-02-01 | 中国电力科学研究院 | 一种受端电网中长期电压稳定性评估方法 |
CN108964010A (zh) * | 2017-05-19 | 2018-12-07 | 国网安徽省电力公司 | 一种确定电网设备对电网安全指标的灵敏度的方法和*** |
CN108090272B (zh) * | 2017-12-13 | 2020-11-17 | 广东电网有限责任公司电力科学研究院 | 一种模块化多电平换流器的建模仿真方法及装置 |
CN108879707A (zh) * | 2018-07-10 | 2018-11-23 | 福州大学 | 一种电力***无功补偿地点在线排序方法 |
CN111030196B (zh) * | 2019-12-17 | 2020-12-11 | 清华大学 | 一种基于动态灵敏度的受端电网动态无功储备优化方法 |
CN111193271B (zh) * | 2020-02-17 | 2021-11-12 | 苏州工业园区服务外包职业学院 | 电网无功优化方法、装置、计算机设备及存储介质 |
CN111682571B (zh) * | 2020-05-07 | 2021-11-02 | 山东大学 | 混合多馈入交直流混联***分级协调电压控制方法及*** |
CN112039074B (zh) * | 2020-09-25 | 2022-01-18 | 贵州电网有限责任公司 | 一种在线安全稳定紧急控制策略方式字生成方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5760860B2 (ja) * | 2011-08-30 | 2015-08-12 | 株式会社明電舎 | 無効電力補償装置 |
CN103825281B (zh) * | 2014-01-22 | 2016-02-24 | 清华大学 | 基于动态无功设备的电力***中中枢母线电压的控制方法 |
CN104037775B (zh) * | 2014-05-14 | 2016-01-27 | 浙江大学 | 一种电力***短期电压稳定控制方法 |
-
2014
- 2014-10-27 CN CN201410584184.0A patent/CN104362642B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104362642A (zh) | 2015-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104362642B (zh) | 提高交直流电网中长期电压稳定的动态无功备用优化方法 | |
CN103701140B (zh) | 提高交直流电网暂态电压稳定的动态无功备用优化方法 | |
Aamir et al. | Impact of synchronous condenser on the dynamic behavior of LCC-based UHVDC system hierarchically connected to AC system | |
Chaudhuri et al. | Adaptive droop control for effective power sharing in multi-terminal DC (MTDC) grids | |
CN104466984B (zh) | 一种提高直流换相安全水平的动态无功备用优化方法 | |
Byeon et al. | A research on the characteristics of fault current of DC distribution system and AC distribution system | |
Chaudhuri et al. | Modeling and stability analysis of MTDC grids for offshore wind farms: A case study on the North Sea benchmark system | |
Zhang et al. | The reactive power voltage control strategy of PV systems in low-voltage string lines | |
Zhang et al. | Generalized short circuit ratio for multi-infeed LCC-HVDC systems | |
Huang et al. | Improving photovoltaic and electric vehicle penetration in distribution grids with smart transformer | |
Isozaki et al. | On detection of cyber attacks against voltage control in distribution power grids | |
Gaigowal et al. | Distributed power flow controller using single phase DSSC to realize active power flow control through transmission line | |
Jayawardena et al. | Low-voltage ride-through characteristics of microgrids with distribution static synchronous compensator (DSTATCOM) | |
Aarathi et al. | Grid connected photovoltaic system with super capacitor energy storage and STATCOM for power system stability enhancement | |
Safitri et al. | Different techniques for simultaneouly increasing the penetration level of rooftop PVs in residential LV networks and improving voltage profile | |
CN106340906A (zh) | 基于改进遗传算法的交直流***低压减载优化方法 | |
Lin et al. | A linearized branch flow model considering line shunts for radial distribution systems and its application in Volt/VAr control | |
Zhao et al. | Optimal Configuration of ESS and SVG for the Coordinated Improvement of Power Quality in Low Voltage Distribution Network with high Penetration PV | |
Gandhar et al. | Application of SSSC for compensation assessment of interconnected power system | |
Shinde et al. | Investigation of effects of solar photovoltaic penetration in an IEEE 13-bus radial low-voltage distribution feeder system | |
Garces et al. | A voltage regulator based on matrix converter for smart grid applications | |
Rijesh et al. | Performance analysis of smart device—STATCOM for grid application | |
Safitri et al. | Coordination of single-phase rooftop PVs to regulate voltage profiles of unbalanced residential feeders | |
Bernal et al. | Reactive power fluctuations smoothing in optimal control of grid-connected PV systems | |
Dias et al. | Power electronics in the context of renewables, power quality and smart grids |
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 |