《1 引言》
1 引言
随着我国电力事业的大发展, 变电所供电半径的优选, 是一个有实际价值的问题。目前, 已建立了静态和动态的优化数学模型
《2 Fuzzy几何规划模型的建立》
2 Fuzzy几何规划模型的建立
考虑城市110 kV变电所直接降压为10 kV向用户供电, 采用年费用法建立模型。设
1) 供电范围是以110 kV变电所为中心, R为半径的圆内, 且10 kV配电网为辐射形结构;
2) 不考虑投资过程, 投资费用是在一定负荷水平下发生的年运行费μ;
3) 投资回收期确定为抵偿年限N (8~10年) , 即总投资在抵偿年限内回收。
4) 整个电网覆盖面上负荷密度均匀。
于是, 用年费用法建立的静态模型如下:
其中, F为年费用, Z为总投资费用。
单位容量的年费用函数可表示为
根据文献
为110 kV变电所建设投资费。
为10 kV中压配电网主干线建设投资费。其中:M=Scosφ /Pav为10kV中压配电网出线回路数, L=ER (km) 为10 kV中压配电网出线每回路线路长度, E为地形修正系数, Sl为110 kV变电所全部的10 kV中压配电网主干线导线总截面积 (mm2) , l是导线的长度 (km) 。
为110 kV变电所和10 kV中压配电网线路的运行费用的不变部分 (大修、小修、折旧费) 与总投资的正比例函数。
为全年10 kV线路损失费。
而110 kV变电所变压器的损失费用为
其中, γ为变压器的台数。
将式 (2) ~ (6) 代入式 (1) , 得
式中:H为对总投资提取的年运行费的年提取系数;S为变电所容量 (kVA) ;a1为投资中与变电所容量无关部分的系数 (元) ;a2为10 kV线路投资中与导线截面无关的每公里的投资 (元/km) ;b1为投资中与变电所容量有关部分的系数 (元/kVA) ;b2为10 kV线路投资中与导线截面有关部分的系数 (元/km·mm2) ;σ为平均负荷密度 (kW/km2) ;Kc为容截比 (110 kV变电所宜取2.2~2.5) ; Pav为每一回路线路平均负荷 (kW) ;UN为中压配电网额定电压 (本文为10 kV) ;J为导线经济电流密度 (A/mm2) ;ρ为导线电阻率 (Ω·mm2/km) ;τ为年平均最大负荷损耗小时数;C0为损耗电度成本电费 (元/kW·h) ;KFe为变压器的铁损系数 (≈0.0085 kW/ (kVA) 3/4) ;τb为变压器的铜损耗等值小时数;Kcu为变压器的铜损系数 (≈0.055 kW/ (kVA) 3/4) ;SNL为每台变压器的额定负荷 (kVA) ;SN为每台变压器额定容量 (kVA) ;静态模型的目标函数确定为单位容量的年费用最小, 即min F0 。
由于负荷、投资过程和电价等因素是不确定的, 含有大量的随机现象和模糊现象。然而, 经典的模型忽略这些重要的现象, 这会使模型失真, 且其适应性将受到很大的限制。这里将目标函数F0作模糊化处理, 使单位容量的年费用模糊最小, 即
《图1》
但变电所的容量是非负的, 即S>0。F0是关于S的指数多项式函数
《图2》
《图3》
即为Fuzzy正项几何规划模型。其中目标函数F0的系数既可为普通正实数, 也可为在某一区间内变化的正Fuzzy数。
获得变电所平均供电半径R, 单位面积上变电所的个数nb和变电所的平均个数Nb。这里Q为城市电网供电区域面积 (km2) 。
《3 模型的求解》
3 模型的求解
N=10, cosφ=0.9, τ =2 400 h, E =1.3, C0 = 0.06元/ kW·h, H = 8.0 %, ρ =31.5 Ω·mm2/km, UN=10 kV, J =1.15 A/mm2, a1=510 000元, b1=50元/ kVA, a2 = 18 640元/ km, b2 = 89元 / km·mm2, Pav=4 000 kW, Kc = 2.2, Q = 49.03 km2, KFe =0.008 5 kW/ (kVA) 3/4, Kcu= 0.055 kW (kVA) 3/4。
将它们代入式 (1) , 便得
代入式 (8) , 即为
《图4》
由于误差, 例如变电所和中压配电线的投资费用系数的波动及电价随时间变化的影响等因素, 目标函数具有一定的伸缩性。而σ是一个参变数, 单位为kW/km2, 经测定宜取:1 566.3, 3 234.9, 4 877.6, 5 915, 8 300, 11 044。
本文取定参数σ=5 915 kW/km2, 且以式 (12) 为例来讨论它的两种情形的求解:模型Ⅰ, 目标函数在某允许域内伸缩;模型Ⅱ, 目标函数的系数在某一区间内波动。
称
为原几何规划及其对偶规划的困难度
《3.1模型Ⅰ》
3.1模型Ⅰ
规划模型为
《图5》
1) 解确定性几何规划
它的对偶规划为
最优解为w1 = 1/3, w2 = 2/3;目标最优值为d0=3.9571。从而得到确定性几何规划式 (14) 的近似最优解为S0≈69 595 kVA。
2) 解确定性几何规划
其中t1=1 725是根据试验, 且由公式
《图6》
确定, 这里k=1/λ, λ为在区间[0, 1]中适当选取的数。
求解规划式 (15) 得, d1=3.960 9。
3) 解确定性几何规划
其中Z0=d1-d0=0.003 8。将其化为
它的项数J=4, 变量个数m=2, 困难度D=4-2-1 , 对偶规划为
式 (17) 的一般解为
这是一维线性约束优化问题。如取w02=2/3, 得
对偶规划的目标函数为
利用式子
即得α=0.337, 从而
获得式 (13) 的最小值为3.960 9。利用公式
即得到模型Ⅰ的近似最优解为S= 69 629, 近似Fuzzy最优值为3.960 9。由式 (9) ~ (11) , 得:
当σ取5 915, 110 kV变电所经济容量S=69 629∈[69 530, 69 728] 时, 经济供电半径R≈1.31 km。该市需建设110 kV变电所计算值9.1个。
《3.2模型Ⅱ》
3.2模型Ⅱ
假定模型中函数的系数为在区间内变化的Fuzzy数。设式 (12) 的目标函数的第二个系数
由式 (19) 有φ (c2) ≥α, 则得
此为含Fuzzy系数的正项几何规划, 即为模型Ⅱ, 它的对偶形式为
解得w0 =1/3, w1 = 2/3。当取α=0.8时, 有d≈3.688 7。利用模型Ⅰ的求解步骤3中的式 (18) , 获得式 (20) 的最优解为S=74 660.45 。S和d亦即为模型Ⅱ的近似最优解和最优值。且当取α=0.8时, 由式 (9) ~ (11) , 得R=1.35, Nb≈8.56。
《4 结果比较分析》
4 结果比较分析
将静态模型、动态模型
Table 1 Comparison of S, R and Nb
《表1》
σ=5 915 kW·km-2 | S/kVA | R/km | Nb |
静态 | 66 100 | 1.27 | 9.68 |
动 态 | 71 320 | 1.32 | 8.96 |
几何规划 | 69 595 | 1.27 | 9.2 |
模型Ⅰ | 69 629 | 1.31 | 9.09 |
模型Ⅱ | 74 660.45 | 1.35 | 8.56 |
对于不同的σj可计算出110 kV变电所的经济容量Sj, 经济供电半径Rj和变电所的个数Nb 。
实际上, 该市在2001年前将规划建设6~10个110 kV变电所, 上述结果与实际比较符合, 可应用于城市电网规划。但从表1中看到, 动态模型优于静态模型, Fuzzy几何规划模型Ⅰ和Ⅱ明显优于静态模型和几何规划模型, 甚至比动态模型更优。同时, 对于高负荷水平, 表1中变电所选的容量偏小, 变电所的个数偏多, 这主要是没有考虑负荷的变化增长过程及变电所建设的投资过程 (土地征用费及配电线的建设费等的增长) , 使得经典模型及其方法在城市电网建设中执行起来比较困难。但用模型Ⅰ和Ⅱ考虑目标的伸缩、系数的波动 (当然, 约束条件的伸缩与其系数的波动亦可作类似的讨论, 在此略去) 。这样模型可包含更多的信息, 给人们的决策留有充分的余地。决策者可根据实际需求进行调整。例如, 可选择经济容量偏大, 变电所个数偏少的最优方案, 以期达到最佳的效果。
如果再考虑变电所和中压配电线路的投资费用系数, 以及电价随时间变化等因素的影响, 其效果将会更好些。
《5 结束语》
5 结束语
本文建立的模型, 对于确定的阈值α和城区的不同平均负荷密度σj , 在变电所允许的经济容量Sj下, 可以很快求出平均供电半径Rj 。当Fuzzy几何规划的困难度D>0时, 可用搜索法求相应的形如式 (16) 、式 (17) 的参数几何规划最优解。这样, 利用Fuzzy几何规划, 建立一个系统控制模型, 根据实际需求调整其参数α, σj 和Sj (j=1, …, r) , 使之在变电所的最大经济容量之下, 求得最经济的供电半径, 而使年费用最小。