主权项 |
一种基于粒子群优化与LCC的电力变压器检修决策方法,其特征在于:包括以下具体步骤:步骤1:初始化LCC两次大修时间:第一次和第二次大修时间t<sub>1</sub>﹑t<sub>2</sub>的初值取值范围为分别为[5,30]和[15,30];步骤2:计算LCC初值y1:y1=W1+W2+W3+W4+W5 (1)其中,W1为初始投资成本、W2为运行成本、W3为断电成本、W4为报废成本、W5为维修成本;所述初始投资成本W1的估计方法为:W1=W<sub>gz</sub>+W<sub>az</sub>+W<sub>zy</sub>+W<sub>yg</sub>+W<sub>ts</sub> (2)其中,W<sub>gz</sub>代表设备购置费,W<sub>az</sub>代表安装工程费,W<sub>zy</sub>代表征用土地费,W<sub>yg</sub>代表员工培训费,W<sub>ts</sub>代表调试费用;所述运行成本W<sub>2</sub>的估计方法为:<img file="re-FDA0001181871930000015.GIF" wi="1264" he="68" />其中,a为电价,单位为元/kWh;p<sub>0</sub>为变压器空载损耗,单位为kW,β为平均负载率,p<sub>k</sub>为变压器有载损耗,单位为kW;i为折现率;t<sub>n</sub>为名义役龄,t<sub>a</sub>为等效役龄,其计算方法为:<img file="re-FDA0001181871930000011.GIF" wi="1346" he="221" />其中,α<sub>3</sub>为大修的役龄回退因子,t<sub>n</sub>为名义役龄。断电成本W3的估计方法为:<img file="re-FDA0001181871930000013.GIF" wi="1620" he="71" />式中:b为单位电量售电利润,元/kWh;S<sub>N</sub>为变压器的额定容量,kVA;β为平均负载率,kW;<img file="re-FDA0001181871930000014.GIF" wi="113" he="41" />为平均功率因数;T为设备年单位故障中断供电时间,h;λ′(t<sub>n</sub>)为设备年平均故障数,次/年;RC为单位时间平均修复成本,万元/小时;MTTR为设备平均修复时间,小时;λ'(t<sub>n</sub>)为故障率函数,计算方法为:<img file="re-FDA0001181871930000012.GIF" wi="1429" he="249" />报废成本W4的估计方法为:W4=(0.32×W<sub>az</sub>‑0.05×W1)×(1+i)<sup>‑40</sup> (7)其中,W<sub>az</sub>为安装工程费,W1为初始投资成本,i为折现率;维修成本W5的估计方法为:W5=M<sub>l</sub>+M<sub>s</sub> (8)其中,M<sub>s</sub>为变压器预防性维修费用,M<sub>l</sub>为两次大修费用之和:<img file="re-FDA0001181871930000021.GIF" wi="1232" he="135" />所述变压器预防性维修费用M<sub>s</sub>的估计方法为:<img file="re-FDA0001181871930000022.GIF" wi="1350" he="470" />其中,β为预防性维修系数;t<sub>n</sub>为名义役龄;W1为初始投资成本;步骤3:采用粒子群优化方法求解变压器全寿命周期成本估计模型,确定最佳维修年限,由以下具体步骤组成:步骤3‑1:初始化变异速度,步骤2计算的LCC,作为最优变压器全寿命周期成本的初值;步骤3‑2:计算相应的变压器全寿命周期成本;V<sub>ik</sub>(n+1)=wV<sub>ik</sub>(n)+c<sub>1</sub>·rand·[(5+12*rand)‑t<sub>ik</sub>(n)]+c<sub>2</sub>·rand()·[‑(5+12*rand)+t<sub>ik</sub>(n)]t<sub>ik</sub>(n+1)=t<sub>ik</sub>(n)+V<sub>ik</sub>(n+1) (11)式中V<sub>ik</sub>(n)为第n次变异速度,V<sub>ik</sub>(n+1)为第n+1次变异速度;最大变异速度为12;rand()为0到1之间的随机数;t<sub>ik</sub>(n)为第n两次大修时间;惯性因子w=0.6;加速系数c<sub>1</sub>=c<sub>2</sub>=2;步骤3‑3:判断所述变压器全寿命周期成本是否小于LCC的当前值,如果是,转向步骤3‑4;否则转向步骤3‑6;步骤3‑4:更新第一次和第二次大修时间t<sub>1</sub>﹑t<sub>2</sub>、变异速度与LCC的当前值;步骤3‑5:判断是否满足算法终止条件,如果是,转向步骤3‑6,否则转向步骤3‑2;步骤3‑6:输出第一次和第二次大修时间t<sub>1</sub>﹑t<sub>2</sub>与LCC的当前值。 |