1.一种计及故障风险及成本的配电网无功规划优化模型,其特征在于,所述模型的建模步骤包括:
步骤1:将配电网每小时负荷曲线的有功功率损耗P、无功功率损耗Q作为特征参数,采用模糊聚类方法将配电线路负荷分为最大、一般、最小三个负荷水平;
步骤2:在确定安装节点的前提下,以各节点的无功补偿容量为优化变量,建立考虑故障风险影响下的无功补偿设备的全寿命周期成本模型CLCC,所述全寿命周期成本模型由投资成本CIC、运行成本COC、检修维护成本CMC、运行条件X下对应的故障风险Risk形成的故障风险成本R、退役处置成本CDC构成:
CLCC=CIC+COC+CMC+Risk(R|X)+CDC
并以CLCC最小为改进的遗传算法的目标函数;
步骤3:根据步骤2求得的最大负荷和最小负荷运行方式下的补偿容量来确定无功补偿设备的安装容量和投切档数;
所述无功补偿设备的全寿命周期成本模型CLCC模型的等式约束为:
(1)CIC为投资成本:
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>I</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>c</mi>
</msub>
</munderover>
<msub>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>c</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>c</mi>
</msub>
<msub>
<mi>N</mi>
<mi>c</mi>
</msub>
</mrow>
式中:Qcj为最大负荷水平下节点j处的无功补偿总容量,单位为kVar;Ccj为节点j处补偿电容器单位容量价格,单位为万元/kVar;Kc为每个节点安装补偿设备的设备费用和固定投资费用,单位为万元/节点;Nc为补偿节点总数;
(2)COC为运行成本:
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>O</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>i</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<msub>
<mi>K</mi>
<mi>e</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>T</mi>
<mi>N</mi>
</msub>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>max</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中:T为使用年限,Ke为电价,单位为元/kWh;TN为系统的年运行小时数;μi1、μi2、μi3分别代表一负荷到最大、一般、最小负荷水平的隶属度;Pmax、Pnor、Pmin分别代表最大、一般和最小三种负荷水平下补偿后系统功率损耗,其中:
Pmax=PZ1+K(QZ1-Qp1)
Pnor=PZ2+K(QZ2-Qp2)
Pmin=PZ3+K(QZ3-Qp3)
式中:K为无功经济当量,Qp1、Qp2、Qp3分别代表最大、一般、最小负荷水平下的无功补偿容量;PZ1、PZ2、PZ3代表最大、一般、最小负荷水平下的有功功率损耗典型向量,QZ1、QZ2、QZ3代表最大、一般、最小负荷水平下的无功功率损耗典型向量;
(3)CMC为检修维护成本:
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>M</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>i</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<msub>
<mi>C</mi>
<mi>M</mi>
</msub>
</mrow>
式中:CM为年检修维护成本;
(4)Risk(R|X)为在运行条件X下对应的故障风险值:
<mrow>
<mi>R</mi>
<mi>i</mi>
<mi>s</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>|</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msubsup>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
</msubsup>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>|</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&theta;</mi>
<mi>d</mi>
<mi>t</mi>
</mrow>
式中:t0为初始时间,t1为终止时间;θ为条件X下,设备的过载热点温度;θ0为过载热点温度初值,θ1为过载热点温度终值;Pr(θ|X)是考虑观测变量X后,变量θ的概率密度函数;H(t|θ)为无功补偿设备在热点温度达到θ时发生故障的概率;R(θ)为无功补偿设备在过载热点温度θ的故障风险成本;Risk(R|X)中R为无功补偿设备的故障风险成本;
(5)CDC为退役处置成本:
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>i</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>S</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中:CSC为报废成本,CR为无功补偿设备的残值;
(6)各节点的功率平衡条件:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<msubsup>
<mi>U</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>cos&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>sin&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>;</mo>
</mrow>
<mrow>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<msubsup>
<mi>U</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>sin&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>cos&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>;</mo>
</mrow>
式中:k=1,2,3代表最大、一般、最小负荷水平;PGi、QGi为加在节点i的有功功率、无功功率;PDi、QDi为节点i负荷的有功功率、无功功率;Qci为节点i无功补偿装置补偿的无功功率;Ui、Uj分别为节点i和节点j的电压;Gij、Bij分别为节点i、j之间导纳矩阵的电导和电纳;θij为节点i、j之间的相角差;yi=1表示节点i为无功补偿点。