CN106355297A - 一种基于电解铝负荷特性的电网降损优化方法 - Google Patents
一种基于电解铝负荷特性的电网降损优化方法 Download PDFInfo
- Publication number
- CN106355297A CN106355297A CN201610888439.1A CN201610888439A CN106355297A CN 106355297 A CN106355297 A CN 106355297A CN 201610888439 A CN201610888439 A CN 201610888439A CN 106355297 A CN106355297 A CN 106355297A
- Authority
- CN
- China
- Prior art keywords
- load
- represent
- electrolytic aluminium
- power
- particle
- 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
Links
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 title claims abstract description 81
- 229910052782 aluminium Inorganic materials 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005457 optimization Methods 0.000 title claims abstract description 15
- 230000009467 reduction Effects 0.000 title abstract description 5
- 239000002245 particle Substances 0.000 claims abstract description 55
- 238000013178 mathematical model Methods 0.000 claims abstract description 8
- 239000004411 aluminium Substances 0.000 claims description 72
- 230000004308 accommodation Effects 0.000 claims description 13
- 230000005611 electricity Effects 0.000 claims description 12
- 230000009194 climbing Effects 0.000 claims description 11
- AZDRQVAHHNSJOQ-UHFFFAOYSA-N alumane Chemical group [AlH3] AZDRQVAHHNSJOQ-UHFFFAOYSA-N 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims description 3
- 230000001960 triggered effect Effects 0.000 claims description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 2
- 240000002853 Nelumbo nucifera Species 0.000 claims description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 2
- 238000002347 injection Methods 0.000 claims 1
- 239000007924 injection Substances 0.000 claims 1
- 238000005868 electrolysis reaction Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- PNEYBMLMFCGWSK-UHFFFAOYSA-N aluminium oxide Inorganic materials [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于电解铝负荷特性的电网降损优化方法,其特征在于,包括如下步骤:步骤1、对电解铝负荷的物理特性进行分析,并结合历史数据得到电解铝负荷的可调节特性;步骤2、在电解铝负荷可调节特性的基础上,建立电解铝负荷调节特性数学模型;步骤3、建立电网降损优化模型,包括目标函数和约束条件;步骤4、利用粒子群算法对电解铝负荷的功率和常规电源的出力进行优化;步骤5、输出优化结果。解决新能源并网电网损耗增加的问题。
Description
技术领域
本发明涉及一种基于电解铝负荷特性的电网降损优化方法。
背景技术
进入21世纪后,全球范围内的能源短缺趋势越发明显,以开发新能源和新技术相结合为基础,建立一个新型的清洁、安全、可持续能源系统成为当前能源发展的趋势。然而随着我国新能源产业的快速发展,大规模新能源基地也面临着电网消纳水平有限、外送通道狭窄、外送能力受阻等一系列技术困境。而且大规模的新能源远距离输送增加了网损,降低了地区电网及主网的电压质量,这又进一步增加了网损。此外,由于新能源与常规电源的逆向分布特点,在常规电源调峰时电网损耗巨大,新能源的随机波动导致输电通道电压大幅波动,也会增加输电通道的损耗。
现有技术中,第一种是利用数据挖掘技术,根据风光出力特性采用不同场景数据划分方法对风光资源丰富地区的实际风电、光电出力数据进行场景划分,在划分的典型场景数据下对风光出力互补耦合特性进行分析,研究合成出力跟踪系统负荷的机制特性及提高出力预测精度问题,提出耦合度和跟踪负荷度计算方法。但是该方法只考虑了不同类型新能源之间的耦合,无法解决大规模新能源无法消纳的问题。
第二种是通过优化常规火电机组出力,实现动态经济调度。该方法在优化模型中引入了正、负旋转备用约束,以应对风电功率预测误差给系统调度带来的影响。主要针对优化模型的不可微和多峰值特性,采用光滑化处理技术,并提出粒子群内点混合优化策略。但是该方法只考虑到火电机组的动态经济调度来减轻风电给系统调度带来的影响,没有考虑风电给电网损耗带来的影响。
第三种是基于需求响应消纳风电的积极效用,综合考虑不确定性调度决策的特点,分别建立了基于实时电价的考虑不确定性的风电消纳确定性和随机机组组合模型,并建立了经济调度模型评估2类机组组合决策的效果,但该方法未对需求侧资源进行详细的分析,也未提出针对性的方法。
发明内容
针对上述问题,本发明提供一种基于电解铝负荷特性的电网降损优化方法,解决新能源并网电网损耗增加的问题。
为实现上述技术目的,达到上述技术效果,本发明通过以下技术方案实现:
一种基于电解铝负荷特性的电网降损优化方法,其特征在于,包括如下步骤:
步骤1、对电解铝负荷的物理特性进行分析,并结合历史数据得到电解铝负荷的可调节特性;
步骤2、在电解铝负荷可调节特性的基础上,建立电解铝负荷调节特性数学模型,如下式所示:
Pmin≤P≤Pmax
cmin≤ci≤cmax
式中,P代表负荷功率,P0代表额定负荷或调节前的负荷,t代表时间,ton代表负荷调节的触发时间,Ru代表增负荷爬坡率,Pt代表t时刻的负荷功率,tret代表返回时间或是下一次调节时间,Rd代表减负荷爬坡率,Pmin代表负荷调节范围最小值,Pmax代表负荷调节范围最大值,ci代表负荷产量,T代表时间总长,pi,t代表负荷在t时刻消耗的有功功率,ηi代表单位产品耗电功率,cmin代表负荷产量范围最小值,cmax代表负荷产量范围最大值;
步骤3、建立电网降损优化模型,包括目标函数和约束条件,其中:
取电网损耗作为目标函数:
式中,Pi,Qi,Ui分别为支路i首端的有功功率、无功功率及电压,RLi为线路阻抗,Pj,Qj,Uj分别为支路j首端的有功功率、无功无功功率及电压,RTj为变压器阻抗,N为线路总数,M为变压器总数,T为时间总长,f为目标函数;
约束条件包括:
式中,∑PG、∑PL、∑PLOSS分别为系统的总发电出力、总负荷和总网损,其中总发电出力包括风力发电出力和常规电源出力;PLmin代表总负荷调节范围最小值,PLmax代表总负荷调节范围最大值,PL为总负荷功率;ΔPmaxdown为单次调节最大减负荷量,用负数表示,ΔPmaxup为单次调节最大增负荷量,Padj(w)为w调节时刻的负荷功率,Padj(w-1)为(w-1)调节时刻的负荷功率;PGmin代表常规电源调节范围最小值,PGmax代表常规电源调节范围最大值,PG为发电机出力功率;
步骤4、利用粒子群算法对电解铝负荷的功率和常规电源的出力进行优化;
步骤5、输出优化结果。
优选,步骤4具体包括如下步骤:
步骤401、对电解铝负荷的功率和常规电源的出力进行编码,编码的值即为电解铝负荷向电网注入的功率和常规电源发出的有功功率;
步骤402、取电解铝负荷和常规电源节点总数α为粒子群维数,对β个粒子进行编码,得α*β的矩阵变量;同时,对每个粒子的速度进行初始化,将个体的历史最优zbest设为当前位置,群体中最优者为当前的gbest;
步骤403、由电解铝负荷的功率和常规电源的出力编码组成的种群为自变量、网损值为因变量,并将网损值作为适应度函数,通过潮流计算,得到各个粒子的适应度函数值;
步骤404、
如果该粒子当前的适应度函数值比历史最优值好,那么历史最优值将被当前位置所替代;
如果该粒子的历史最优值比全局最优值好,那么全局最优值将被该粒子的历史最优值所替代;
步骤405、对每个粒子按照下列公式对速度和位置进行更新:
vi(k+1)=ψvi(k)+c1r1(pτ-xi(k))+c2r2(pg-xi(k))
xi(k+1)=xi(k)+vi(k+1)
式中,pτ表示个体极值,pg表示全体极值,ψ,c1,c2是常数,vi(k),vi(k+1)分别是第i个粒子第k步、k+1步迭代时当前粒子的速度,xi(k)、xi(k+1)分别是第i个粒子第k步、k+1步迭代时的位置,r1是随机数参数,r2是随机参数;
步骤406、迭代次数加1,如果未达到最大迭代次数,则返回至步骤403,否则输出结果并结束。
本发明的有益效果是:
本发明在分析了电解铝负荷调节特性的基础上,构建了电解铝负荷调节特性的数学模型,并在电解铝负荷调节特性数学模型的基础上,以电网损耗为优化目标,考虑电解铝负荷调节和常规电源出力的约束,建立了基于电解铝负荷调节特性的电网降损优化模型。通过粒子群算法对本发明提出的电网降损优化模型进行优化,得到了优化结果,并实例证明了通过对本发明得到的基于电解铝负荷调节特性的电网降损优化模型进行优化,能够很好的降低电网的损耗。
附图说明
图1是本发明典型的电解铝日负荷曲线;
图2是本发明典型电解铝负荷调节曲线示意图;
图3是本发明一种基于电解铝负荷特性的电网降损优化方法的流程图;
图4是本发明实施例中甘肃电网的结构示意图;
图5是本发明常规电源调节的降损优化结果;
图6是本发明考虑电解铝负荷调节的降损优化结果。
具体实施方式
下面结合附图和具体的实施例对本发明技术方案作进一步的详细描述,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
一种基于电解铝负荷特性的电网降损优化方法,包括如下步骤:
步骤1、对电解铝负荷的物理特性进行分析,并结合历史数据得到电解铝负荷的可调节特性:
目前,电解铝工业生产采用冰晶石-氧化铝熔盐电解法,通入强大的直流电后,在高温下,在电解槽内的两极上进行电化学反应,电解出铝金属。铝电解生产过程依靠不断地供给电解槽直流电能,工艺生产要求恒定的直流电流,若直流电流大幅度波动,会显著影响电解槽运行状态(尤其是热平衡)的平稳性,并干扰铝电解的正常控制过程,从而显著影响铝电解槽的技术经济指标。因此,对电解铝负荷进行调节时,要在电解槽能够承受的范围内压低负荷生产,通过控制可控硅的导通角调节电解铝负荷电流大小,降低系列电流运行。
电解铝负荷在正常情况下,负荷容量变化较小,负荷曲线较为平稳,如附图1所示。一般的,电解铝负荷通常有多条馈线,分别连接一台整流机组,每台机组均可在一定范围内进行调节。由现有数据可知,电解铝负荷调节范围为-10%~10%。通过步骤1获得可调节特性,包括负荷功率的调节范围,负荷调节的单次调节范围。
步骤2、在电解铝负荷可调节特性的基础上,建立电解铝负荷调节特性数学模型:
在电解铝负荷可调状态下,负荷开始参与调节时,负荷按照爬坡率从额定值开始变化,直至设定值为止。从额定值往下调时负荷爬坡率称为减负荷爬坡率,往上调节时称为增负荷爬坡率,并且不同机组的电解铝负荷爬坡率存在差异。
为不影响电解铝负荷的正常生产,要保证电解铝负荷的产量在允许波动范围内,负荷调节设定值需限制在可调范围内。以单次调节为例,假定由额定值调节到某设定值再返回到额定值,可得典型的电解铝负荷调节曲线如附图2所示。由上述调节曲线可建立电解铝负荷调节特性的数学模型,具体如下:
Pmin≤P≤Pmax
cmin≤ci≤cmax
式中,P代表负荷功率,P0代表额定负荷或调节前的负荷,t代表时间,ton代表负荷调节的触发时间,Ru代表增负荷爬坡率,Pt代表t时刻的负荷功率,tret代表返回时间或是下一次调节时间,Rd代表减负荷爬坡率,Pmin代表负荷调节范围最小值,Pmax代表负荷调节范围最大值,ci代表负荷产量,T代表时间总长,pi,t代表负荷在t时刻消耗的有功功率,ηi代表单位产品耗电功率,cmin代表负荷产量范围最小值,cmax代表负荷产量范围最大值。
步骤3、建立电网降损优化模型,包括目标函数和约束条件,其中:
取电网损耗作为目标函数:
式中,Pi,Qi,Ui分别为支路i首端的有功功率、无功功率及电压,RLi为线路阻抗,Pj,Qj,Uj分别为支路j首端的有功功率、无功无功功率及电压,RTj为变压器阻抗,N为线路总数,M为变压器总数,T为时间总长,f为目标函数;
分别考虑电网有功功率的约束、电解铝负荷的调节特性、负荷调节的爬坡速率约束和常规电源出力约束,约束条件包括:
式中,∑PG、∑PL、∑PLOSS分别为系统的总发电出力、总负荷和总网损,其中总发电出力包括风力发电出力和常规电源出力;PLmin代表总负荷调节范围最小值,PLmax代表总负荷调节范围最大值,PL为总负荷功率;ΔPmaxdown为单次调节最大减负荷量,用负数表示,ΔPmaxup为单次调节最大增负荷量,Padj(w)为w调节时刻的负荷功率,Padj(w-1)为(w-1)调节时刻的负荷功率;PGmin代表常规电源调节范围最小值,PGmax代表常规电源调节范围最大值,PG为发电机出力功率。
步骤4、利用粒子群算法对电解铝负荷的功率和常规电源(火力发电机,水力发电机等)的出力进行优化:
优选,如图3所示,具体包括如下步骤:
步骤401、对电解铝负荷的功率和常规电源的出力进行编码,编码的值即为电解铝负荷向电网注入的功率和常规电源发出的有功功率;
步骤402、取电解铝负荷和常规电源节点总数α为粒子群维数,对β个粒子进行编码,得α*β的矩阵变量;同时,对每个粒子的速度进行初始化,将个体的历史最优zbest(zbest即算法上次计算结果的pBest)设为当前位置,群体中最优者为当前的gbest;
粒子群算法初始化为一群随机粒子(随机解),然后通过迭代找到最优解,在每一次迭代中,粒子通过跟踪两个“极值”来更新自己:第一个就是粒子本身所找到的最优解,这个解叫做个体极值pBest;另一个极值是整个种群目前找到的最优解,这个极值是全局极值gBest。一般的,依次运用随机函数对β个粒子进行编码,同时,运用随机函数对每个粒子的速度进行初始化。
步骤403、由电解铝负荷的功率和常规电源的出力编码组成的种群为自变量、网损值为因变量,并将网损值作为适应度函数,通过潮流计算,得到各个粒子的适应度函数值;
步骤404、
如果该粒子当前的适应度函数值比历史最优值好,那么历史最优值将被当前位置所替代;
如果该粒子的历史最优值比全局最优值好,那么全局最优值将被该粒子的历史最优值所替代;
步骤405、对每个粒子按照下列公式对速度和位置进行更新:
vi(k+1)=ψvi(k)+c1r1(pτ-xi(k))+c2r2(pg-xi(k))
xi(k+1)=xi(k)+vi(k+1)
式中,pτ表示个体极值,pg表示全体极值,ψ,c1,c2是常数,vi(k),vi(k+1)分别是第i个粒子第k步、k+1步迭代时当前粒子的速度,xi(k)、xi(k+1)分别是第i个粒子第k步、k+1步迭代时的位置,r1是随机参数,r2是随机参数。
步骤406、迭代次数加1,如果未达到最大迭代次数,则返回至步骤403,否则输出结果并结束。
步骤5、输出优化结果。
如图4所示,以甘肃电网为基础,其中有2台风电机组,2台风电机组并网功率分别为60MW和70MW,并分别通过节点1和9并网,节点6和节点11所连发电机均为常规电源,节点7接入电解铝负荷,其有功功率为310MW。利用粒子群算法对电解铝负荷功率和常规电源出力进行优化,具体步骤如下:
第一步,对电解铝负荷的功率和常规电源的出力进行编码,编码的值即为电解铝负荷向电网注入的功率和常规电源发出的有功功率,根据给定的负荷确定编码的取值区间。设定粒子群规模为N=50,惯性权重系数的最大值和最小值分别为:ωmax=0.9,ωmin=0.1,加速系数c=2,迭代次数为50等参数。
第二步,取电解铝负荷和常规电源节点总数3为粒子群维数,依次运用随机函数对50个粒子进行编码,可得3*50的矩阵变量,代表了粒子群算法中的种群。同时,运用随机函数对每个粒子的速度进行初始化,将个体的历史最优zbest设为当前位置,群体中最优者为当前的gbest。
第三步,由电解铝负荷功率和常规电源出力编码组成的种群为自变量,网损值为因变量,并将网损值作为适应度函数,通过潮流计算,得到各个粒子的适应度函数值。
第四步,如果该粒子当前的适应度函数值比历史最优值好,那么历史最优值将被当前位置所替代。
第五步,如果该粒子的历史最优值比全局最优值好,那么全局最优值将被该粒子的历史最优值所替代。
第六步,对每个粒子进行粒子的速度和位置更新。每次更新后,检查速度是否超过最大值Vmax=2或是最小值Vmim=-2,若超过最大值,则将速度取值为2,若速度小于最小值,则将速度取值为-2,其他情况速度的取值不变。
第七步,更新粒子的飞行对粒子进行更新,并重新计算适应度函数值。
第八步,若达到最大迭代次数,则停止计算;否则返回第六步。
由上述步骤可得优化结果,即电解铝负荷的功率及常规能源的出力的值。
若仅对常规能源进行调节,则在第二步中仅对常规能源进行编码,可得优化结果如附图5所示,若对常规能源与电解铝负荷同时调节,则可得优化结果如附图6所示。
电解铝负荷具有调节特性,装机容量大,绝大部分采用330/750kV高压直供电模式。基于电解铝负荷的可调节特性,使其与新能源出力配合,将有助于消纳新能源、降低电网损耗。本发明结合电解铝负荷历史数据,对电解铝负荷物理特性进行分析,说明电解铝负荷具有可调节特性,能够配合新能源出力,在一定程度上解决新能源消纳的问题,并降低由新能源并网带来的电网损耗的增加。且电解铝负荷在参与调节的过程中受到一定的约束,因此,深入分析电解铝负荷的调节特性对于建立电解铝负荷调节特性数学模型具有十分重要的意义。
本发明在分析了电解铝负荷调节特性的基础上,构建了电解铝负荷调节特性的数学模型,并在电解铝负荷调节特性数学模型的基础上,以电网损耗为优化目标,考虑电解铝负荷调节和常规电源出力的约束,建立了基于电解铝负荷调节特性的电网降损优化模型。通过粒子群算法对本发明提出的电网降损优化模型进行优化,得到了优化结果,并实例证明了通过对本发明得到的基于电解铝负荷调节特性的电网降损优化模型进行优化,最大程度的降低了电网损耗。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或者等效流程变换,或者直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (3)
1.一种基于电解铝负荷特性的电网降损优化方法,其特征在于,包括如下步骤:
步骤1、对电解铝负荷的物理特性进行分析,并结合历史数据得到电解铝负荷的可调节特性;
步骤2、在电解铝负荷可调节特性的基础上,建立电解铝负荷调节特性数学模型,如下式所示:
Pmin≤P≤Pmax
cmin≤ci≤cmax
式中,P代表负荷功率,P0代表额定负荷或调节前的负荷,t代表时间,ton代表负荷调节的触发时间,Ru代表增负荷爬坡率,Pt代表t时刻负荷功率,tret代表返回时间或是下一次调节时间,Rd代表减负荷爬坡率,Pmin代表负荷调节范围最小值,Pmax代表负荷调节范围最大值,ci代表负荷产量,T代表时间总长,pi,t代表负荷在t时刻消耗的有功功率,ηi代表单位产品耗电功率,cmin代表负荷产量范围最小值,cmax代表负荷产量范围最大值;
步骤3、建立电网降损优化模型,包括目标函数和约束条件,其中:
取电网损耗作为目标函数:
式中,Pi,Qi,Ui分别为支路i首端的有功功率、无功功率及电压,RLi为线路阻抗,Pj,Qj,Uj分别为支路j首端的有功功率、无功无功功率及电压,RTj为变压器阻抗,N为线路总数,M为变压器总数,T为时间总长,f为目标函数;约束条件包括:
式中,ΣPG、ΣPL、ΣPLOSS分别为系统的总发电出力、总负荷和总网损,其中总发电出力包括风力发电出力和常规电源出力;PLmin代表总负荷调节范围最小值,PLmax代表总负荷调节范围最大值,PL为总负荷功率;ΔPmaxdown为单次调节最大减负荷量,用负数表示,ΔPmaxup为单次调节最大增负荷量,Padj(w)为w调节时刻的负荷功率,Padj(w-1)为(w-1)调节时刻的负荷功率;PGmin代表常规电源调节范围最小值,PGmax代表常规电源调节范围最大值,PG为发电机出力功率;
步骤4、利用粒子群算法对电解铝负荷的功率和常规电源的出力进行优化;
步骤5、输出优化结果。
2.根据权利要求1所述的一种基于电解铝负荷特性的电网降损优化方法,其特征在于,步骤4具体包括如下步骤:
步骤401、对电解铝负荷的功率和常规电源的出力进行编码,编码的值即为电解铝负荷向电网注入的功率和常规电源发出的有功功率;
步骤402、取电解铝负荷和常规电源节点总数α为粒子群维数,对β个粒子进行编码,得α*β的矩阵变量;同时,对每个粒子的速度进行初始化,将个体的历史最优zbest设为当前位置,群体中最优者为当前的gbest;
步骤403、由电解铝负荷的功率和常规电源的出力编码组成的种群为自变量、网损值为因变量,并将网损值作为适应度函数,通过潮流计算,得到各个粒子的适应度函数值;
步骤404、
如果该粒子当前的适应度函数值比历史最优值好,那么历史最优值将被当前位置所替代;
如果该粒子的历史最优值比全局最优值好,那么全局最优值将被该粒子的历史最优值所替代;
步骤405、对每个粒子按照下列公式对速度和位置进行更新:
vi(k+1)=ψvi(k)+c1r1(pτ-xi(k))+c2r2(pg-xi(k))
xi(k+1)=xi(k)+vi(k+1)
式中,pτ表示个体极值,pg表示全体极值,ψ,c1,c2是常数,vi(k),vi(k+1)分别是第i个粒子第k步、k+1步迭代时当前粒子的速度,xi(k)、xi(k+1)分别是第i个粒子第k步、k+1步迭代时的位置,r1是随机参数,r2是随机参数;
步骤406、迭代次数加1,如果未达到最大迭代次数,则返回至步骤403,否则输出结果并结束。
3.根据权利要求2所述的一种基于电解铝负荷特性的电网降损优化方法,其特征在于,步骤402中,依次运用随机函数对β个粒子进行编码,同时,运用随机函数对每个粒子的速度进行初始化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610888439.1A CN106355297B (zh) | 2016-10-11 | 2016-10-11 | 一种基于电解铝负荷特性的电网降损优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610888439.1A CN106355297B (zh) | 2016-10-11 | 2016-10-11 | 一种基于电解铝负荷特性的电网降损优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106355297A true CN106355297A (zh) | 2017-01-25 |
CN106355297B CN106355297B (zh) | 2019-08-23 |
Family
ID=57866188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610888439.1A Active CN106355297B (zh) | 2016-10-11 | 2016-10-11 | 一种基于电解铝负荷特性的电网降损优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106355297B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108808688A (zh) * | 2018-07-04 | 2018-11-13 | 云南电网有限责任公司 | 一种基于电解铝控制特性的稳控切负荷方法 |
CN111737873A (zh) * | 2020-06-24 | 2020-10-02 | 广东电网有限责任公司电力调度控制中心 | 基于电解铝负荷提供辅助服务的机组组合方法及相关装置 |
CN112417792A (zh) * | 2020-11-23 | 2021-02-26 | 云南电网有限责任公司 | 一种基于饱和电抗器可控边界的电解铝外特性建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101354423A (zh) * | 2008-09-16 | 2009-01-28 | 华北电力大学 | 基于实测的冲击负荷建模系统及方法 |
CN102376026A (zh) * | 2011-10-31 | 2012-03-14 | 冶金自动化研究设计院 | 工业企业用电负荷优化系统 |
CN102851703A (zh) * | 2012-10-09 | 2013-01-02 | 国电南瑞科技股份有限公司 | 一种电解铝供电系统智能自动调压优化控制方法 |
CN105870915A (zh) * | 2016-04-06 | 2016-08-17 | 清华大学 | 用于电力系统控制的电解铝和自备机组鲁棒等值建模方法 |
-
2016
- 2016-10-11 CN CN201610888439.1A patent/CN106355297B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101354423A (zh) * | 2008-09-16 | 2009-01-28 | 华北电力大学 | 基于实测的冲击负荷建模系统及方法 |
CN102376026A (zh) * | 2011-10-31 | 2012-03-14 | 冶金自动化研究设计院 | 工业企业用电负荷优化系统 |
CN102851703A (zh) * | 2012-10-09 | 2013-01-02 | 国电南瑞科技股份有限公司 | 一种电解铝供电系统智能自动调压优化控制方法 |
CN105870915A (zh) * | 2016-04-06 | 2016-08-17 | 清华大学 | 用于电力系统控制的电解铝和自备机组鲁棒等值建模方法 |
Non-Patent Citations (1)
Title |
---|
张学兵: "电网降损节能办法探讨", 《管理观察》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108808688A (zh) * | 2018-07-04 | 2018-11-13 | 云南电网有限责任公司 | 一种基于电解铝控制特性的稳控切负荷方法 |
CN108808688B (zh) * | 2018-07-04 | 2021-07-02 | 云南电网有限责任公司 | 一种基于电解铝控制特性的稳控切负荷方法 |
CN111737873A (zh) * | 2020-06-24 | 2020-10-02 | 广东电网有限责任公司电力调度控制中心 | 基于电解铝负荷提供辅助服务的机组组合方法及相关装置 |
CN112417792A (zh) * | 2020-11-23 | 2021-02-26 | 云南电网有限责任公司 | 一种基于饱和电抗器可控边界的电解铝外特性建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106355297B (zh) | 2019-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109066805B (zh) | 一种跨区互联电网发输电系统动态调度优化方法 | |
CN107301472B (zh) | 基于场景分析法和电压调节策略的分布式光伏规划方法 | |
CN108092324B (zh) | 一种风电参与调峰调频的agc控制系统和控制方法 | |
CN109256799B (zh) | 一种基于样本熵的新能源电力系统优化调度方法 | |
CN107292766B (zh) | 面向风电消纳的电力系统调峰手段经济性评估方法与系统 | |
CN110620397B (zh) | 一种高比例可再生能源电力系统的调峰平衡评估方法 | |
CN107959285B (zh) | 一种区域间日前发输电计划的优化方法、系统及装置 | |
CN107480896A (zh) | 一种支持多分区的新能源受端电网消纳能力评估方法 | |
CN103997039A (zh) | 基于概率区间预测的计及风电接纳的旋转备用区间预测方法 | |
CN106355297A (zh) | 一种基于电解铝负荷特性的电网降损优化方法 | |
CN113241803A (zh) | 一种基于新能源消纳的储能调度方法及计算机介质 | |
CN114996908B (zh) | 考虑智能软开关接入的主动配电网扩展规划方法及系统 | |
CN115345380A (zh) | 一种基于人工智能的新能源消纳电力调度方法 | |
CN112072643A (zh) | 一种基于深度确定性梯度策略的光-蓄系统在线调度方法 | |
CN105305485A (zh) | 一种消纳大规模间歇性能源的安全约束经济调度方法 | |
CN105305463B (zh) | 计及光伏发电和谐波污染的基于随机潮流的无功优化方法 | |
CN110429591A (zh) | 一种基于电力系统时序耦合性的输电网利用率评估方法 | |
CN105279575B (zh) | 基于发电预测的多能源主体分布式博弈优化方法 | |
CN116073448B (zh) | 一种基于低碳效益的配电系统源网荷储协同调峰方法 | |
CN116544951A (zh) | 一种水电调峰灵活性量化与调度方法 | |
CN105976046B (zh) | 一种考虑需求侧管理的低碳电网规划方法 | |
CN111767621A (zh) | 一种基于知识迁移q学习算法的多能源系统优化调度方法 | |
CN115271244A (zh) | 一种基于两阶段分布鲁棒优化的梯级水电站短期调峰模型 | |
CN115549118A (zh) | 配电网无功补偿设备出力平滑指数预测方法 | |
CN108039739A (zh) | 一种主动配电网动态随机经济调度方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |