计及低碳排放和需求侧管理的主动配电网双层规划

1、计及低碳排放和需求侧管理的主动配电网双层规划建设新型电力系统，利用清洁可再生能源发电，促进电网低碳发展，是配电网规划工作面临的新使命和新课题。风电、光电等间歇性能源大量并网，使配电网从传统被动单向辐射状的供电网络转变为主动控制功率双向流动的有源网络，即主动配电网络（active distribution netwok，ADN），同时也给配电网规划带来复杂性与不确定性。在电力行业中，碳排放被认为来自上游生产环节，需求导致生产，碳排放的源头归因于用电负荷1。因此，应用需求侧管理技术不仅可以削峰填谷、降低负荷需求，还能够促进配电网低碳发展。可中断负荷作为一种灵活的需求响应资源，如何针对其特性进行建模

2、，并融入到不确定环境下的配电网规划模型中，是目前规划工作面临的挑战。为应对不确定性问题，目前配电网规划主要有区间规划2、模糊规划3，4、概率建模5、场景分析法6，7、风险评估法8和鲁棒优化法9，10。区间规划适用于不确定因素变化小的场景，该方法计算量小，缺点是只能获得潮流区间值，精度比较差；模糊规划适用于不确定因素变化较大的场景，但是只能获得潮流大致分布；概率建模适用于精度要求较高的场景；场景分析法适用于组合场景较少的规划，存在典型场景获取困难的问题；风险评估法针对约束条件不苛刻的场景，存在风险如何量化和置信度如何选取的问题；鲁棒优化针对不确定变量形式多样的场景，计算效率高，但是优化结果较为保

3、守。本文提出一种计及低碳排放和需求侧管理的主动配电网分层规划方法建立的双层模型，上层规划模型考虑包括碳排放成本在内的配电公司各种经济成本之和，即以年综合费用最小化为目标，在此基础上，下层规划以分布式电源（distributed generation，DG）有功出力切除量和可中断负荷中断量最小为目标对每个场景进行优化，采用二进制量子粒子群算法结合原对偶内点法进行求解。考虑ADN中不确定因素如风速、光照强度和负荷，建立相应的概率模型，利用拉丁超立方抽样技术生成配电网全年运行序列。本文在上层规划模型总成本中加入碳排放成本，从而体现所建配电网兼顾运行经济性和排放低碳性；将DG有功出力切除量和可中断负荷

4、中断量联合进行优化，降低峰值负荷，优化用电方式，使DG出力和用户负荷需求匹配得更加紧密，有利于配电网的稳定运行。1 分布式风电、光伏发电及负荷描述方法1.1风电、光伏发电和负荷模型根据已有大量的观测数据和理论模拟研究可知，规划阶段的风速服从两参数韦布尔分布，光照强度服从贝塔分布，负荷服从高斯分布，其概率密度函数分别如下：f(v)=(kck)v(k1)e(vc)k（1）f(r)=(+)()()(rrmax)(1)(1rrmax)(1)（2）f(P)=1(2)exp(P)222（3）式中：v为实际风速；c为尺度参数，反映地方风资源的平均风速；k为形状参数，为无量纲量；r为实际观测光照强度；rmax

5、为最大光照强度，由安装光伏节点位置具体光照情况确定；和为beta分布的2个正值参数；P为电力网有功负荷；为其期望值；2为方差。风电与风速、光伏与光照强度的近似函数关系如下：PW=00vvciPW,nvvcivnvcivcivvnPW,nvnvvco0vcov（4）PP=PP,nrrmaxrrmaxPP,nrmaxr（5）式中：PW，n和PP，n分别为风机和光伏的额定输出功率；vci、vn、vco分别为风机的切入风速、额定风速和切出风速。1.2可中断负荷可中断负荷是一种基于激励的需求响应资源11，在高峰用电负荷加大或紧急状况下，电网自动向用户发出削减负荷的需求侧信号（demand respons

6、e，DR），用户响应后对用电量进行削减或中断以配合电力运营公司的需求响应策略，而配电公司支付用户一定补偿。可中断负荷经济补偿方式主要有2种：固定支出和停电后高赔偿。前者供电公司每年需要支付固定的补偿费用，后者则是在实际停电后按合约上的可中断负荷电价进行赔偿。本文研究的可中断负荷补偿机制为后者，即实施中断后按照负荷中断量及可中断负荷电价给予用户停电补偿。2 主动配电网上下双层规划模型主动配电网规划不仅涉及到DG的选址定容及网络馈线的架设，还需要统筹考虑需求响应和削减DG有功出力、调节有载变压器抽头等管理措施，因此，将ADN规划转换成一个双层规划问题，上层规划预设网架结构和DG容量位置，为下层决策

7、提供初始条件。下层决策对DG有功出力和可中断负荷量联合优化用电方式，维持配电网稳定运行。2.1上层规划模型在上层规划中以配电网综合费用Ctotal最小为优化目标，考虑供电线路的投资及运行维护等年值费用Cline、网损年费用Closs、DG安装及运行维护等年值费用Cdg、购电费用Cenbuy，并在此基础上考虑配电网碳排放成本Ccarbon、主动管理费用CAM和需求侧管理费用CDSM。目标函数为minCtotal：Ctotal=Cline+Cdg+Cenbuy+Closs+Ccarbon+CAM+CDSM（6）1）配电网碳排放成本Ccarbon包括碳排放权交易收益和碳税12：Ccarbon=t=1

8、8760k=1mPk,t+iD|Ii,t|2Rlii=1n(Pdg,i,tPcurt,i,t)tCai=1n(Pdg,i,tPcurt,i,t)tCb(7)式中：Pk，t为负荷点k在t时刻的负荷；Ii，t为线路i在t时刻的电流；D为配电网线路集合；R和li分别为电阻率和线路i的长度；Pdg，i，t和Pcurt，i，t分别为分布式电源i在t时刻的有功出力和削减功率；Ca为单位碳排放的税收；Cb为碳交易价格，即可再生能源发电减排所获收益；为单位电量碳排放量。2）主动管理费用CAM：CAM=t=18760i=1nFAMPdg,i,tt（8）式中：FAM为分布式电源i单位发电量的主动管理成本。3）需求

9、侧管理费用CDSM为配电公司给予用户削减负荷的停电补偿：CDSM=t=18760i=1nFDSMPDSM,i,tt（9）式中：FDSM为可中断负荷电价；PDSM，i，t为第i个被控负荷在t时刻的负荷削减量。4）其他各类费用：Cline=iDr0(1+r0)T1(1+r0)T11clli+ulxi（10）Cdg=h=1nr0(1+r0)T2(1+r0)T21cdPdgh+udh（11）Cenbuy=t=18760cePen,tt（12）Closs=t=18760ceiD|Ii,t|2rlit（13）式中：T1、T2分别为线路和DG的全寿命周期；r0为配电网建设投资贴现率；cl为配电网线路工程的单

10、位长度造价，万元/km；ul为线路运行维护费用；cd为DG单位容量造价；Pdg.h、ud.h分别为DG装机容量和运行维护成本；xi为ADN规划候选线路，为0-1变量；Pen，t为第t时刻配电网在上级电网的购电量；ce为购电单价。2.2上层规划约束条件上层规划包括以下约束条件：1）DG安装节点容量约束：PmindghPdghPmaxdgh（14）2）辐射连通状约束：n=m+1（15）式中：n、m分别为系统节点数和支路数。3）DG渗透率约束为减小ADN节点电压、线路电流越限风险，规定配电网接入DG的容量满足：max（16）式中：为DG的实际渗透率；max为DG的最大渗透率。4）可靠性约束：RASA

11、IRmin（17）RASAI=Ntotal8760UsiNsiNtotal8760（18）式中：RASAI为ADN平均供电可用率；Usi为每年线路分段故障停电的持续时间；Ntotal为ADN用户总数；Nsi为受影响用户数。2.3下层规划模型下层模型从数学角度看属于有约束条件的非线性规划问题，以DG有功出力切除量和可中断负荷中断量最小为优化目标，其表达式为mint=18760iNdgPcurt,i,t+iNDSMPDSM,i,t（19）式中：Ndg为有功负荷切用户数；NDSM为可中断负荷用户数。2.4下层规划约束条件下层规划包括以下约束条件：1）潮流约束条件：Pi=Uij=1nUj(Gijcos

12、ij+Bijsinij)Qi=Uij=1nUj(GijsinijBijcosij) （20）式中：Pi、Qi为节点i注入的有功和无功功率；Ui、Uj为节点i、j的电压幅值；Gij、Bij、ij为节点i和j之间的电导、电纳和相角。2）节点电压幅值上下限约束：UminiUiUmaxi（21）式中：Umini为允许的最低节点电压幅值；Umaxi为允许的最高节点电压幅值。3）配电线路电流Ij约束：IjImaxj（22）式中：Imaxj为配电线路安全运行电流上限值。4）DG有功出力切除量约束：Pmincurt,iPcurt,i,tPmaxcurt,i（23）5）需求侧可中断负荷弹性约束：0PDSM,i,

13、tPmaxDSM,i,t（24）6）无功电源出力约束：QminciQci,tQmaxci（25）式中：Qci,t为电源节点t时刻注入的电网无功；Qmaxci、Qminci为电源节点注入电网无功的上限值和下限值。7）有载调压变压器分接头调节约束：TminkTk,tTmaxk（26）式中：Tmaxk、Tmink为k节点的有载调压变压器分接头的最大和最小档位值。3 模型求解3.1量子粒子群算法为避免粒子群优化（particle swarm opi-timization，PSO）算法在进化过程中过早收敛和陷入局部最优解，Sun等13将量子力学的相关概念引入粒子进化过程，提出量子粒子群算法，提高了算法收

14、敛速度、种群多样性及全局寻优能力，其基本公式如下：pi=pbi+(1)gb（27）mb=i=1MpbiM（28）r(t)=1,rand0.51,rand0.5（29）x(t+1)=pi+r(t)|mbx(t)|In(1u)（30）式中：mb为个体平均最优；pi为每个粒子的局部吸引子；和u为（0，1）之间的随机数；pbi为每个粒子寻优所取得局部最优值；gb为种群历史最优值；为收缩扩张系数，线性减小可以加快粒子群算法的收敛速度；M为种群数量。3.2二进制量子粒子群算法ADN规划属于多离散变量的非线性混合整数规划问题，本文选择二进制（0或1）编码的量子粒子群算法14，分别进行线路编码和DG编码。线路

15、编码部分采用n位二进制数模拟线路投建状态和导线截面选择，假设被选导线型号有3种：LGJ-185、LGJ-240和LGJ-300，则编码取值范围为00，01，10，11，依次代表线路不架设及架设上述3种导线型号中的1种。DG编码部分采用n位二进制数来表示每个候选节点上DG的安装容量，即有2n种可能容量取值，当编码为0时表示不接入DG，其余情况表示接入相应DG的容量，通过DG编码实现DG类型、位置及接入容量的同时优化。量子粒子群算法更新公式如下：bi=dH(pi,mb)In(1/)（31）qmi=bi/dH1（32）x(i,t)=1,qmi00,qmi0（33）式中：bi为变异吸引子；dH（pi，

16、mb）为pi和mb2个粒子间的海明距离，即对应位置上不同字符的个数；qmi为变异概率。3.3规划流程ADN规划流程具体步骤如下：1）输入ADN线路、负荷、DG参数，利用拉丁超立方抽样技术15对风速、光照强度、负荷概率密度函数抽样，根据风机和光伏发电函数得到DG出力样本矩阵。2）初始化粒子群算法参数，对ADN线路和DG决策变量进行编码，生成满足约束条件的初始种群。3）更新量子粒子位置，产生子种群。4）采用原始对偶内点法对下层最优潮流问题进行求解，并将下层规划最优值代入到上层规划中，计算种群个体的适应度并排序，优胜劣汰。5）判断有关的收敛条件是否满足，如果不满足，则返回到第3步，满足输出最优解。4

17、 算例分析以1个待扩展的25节点配电网16为例进行仿真计算，如图1所示。设置负荷节点电压允许偏差为额定电压的5%，节点6、10、18为可中断负荷节点，其他节点为常规节点。线路造价11.9万元/km，电阻为0.145 /km，电抗为0.114 /km，维修时间和故障率分别为2 h和0.12 次/（kma）。待选接入风机的节点为5、7、8、11，待选接入光伏的节点为14、19、24，单台风机和光伏的额定功率均为100 kW，考虑到配电网实际运行情况，各节点最多允许接入4台DG，DG的单位投资费用为8 000 元/kW，运行费用为0.2 元/kW。风速分布形态服从二参数Weibul分布函数，形状参数

18、k为2.2、尺度参数c为8.21。风机切入时，额定和切出速度分别为4、14和25 m/s。太阳辐射强度参数：a=0.65，b=0.85。rmax=1 000 W/m2，ce=0.38 元/（kWh），可中断负荷单位电量中断补偿价格为1.2元/（kWh），t=1 h，=0.78 kg/（kWh），Ca=0.045 元/kg，Cb=0.065 元/kg，采样数为8 760。图125节点配电网Fig.125-node distribution network本文设计了3种不同方案进行ADN规划，并进行了对比分析：方案1采取控制DG的发出功率和可中断负荷2种主动管理措施；方案2只考虑控制DG的发出功率

19、，不考虑应用可中断负荷；方案3既不应用可中断负荷，也不对DG进行主动管理，即“安装便遗忘”。3种方案得到的最优规划方案结构如图2所示。表1为3种规划方案的各相关费用，表2为3种规划方案的DG位置和容量对比。图23种规划模型优化供电方案结果对比Fig.2Comparison of optimization results of power supply schemes of three planning models表1配电网规划费用Table 1Distribution network planning costs方案规划费用/万元ClineCdgClossCenbuyCcarbonCAMCD

20、SMCtotal1367.13341.2445.45326.10121.6734.5125.741 261.842448.11280.2973.34352.30178.4325.1201 357.593476.71240.2989.74398.98218.78001 424.50表23种规划方案DG位置和容量对比Table 2Comparison of DG location and capacity of 3 planning schemes方案风力发电光伏发电位置容量/kW位置容量/kW15,7,8,11200,300,400,20224,19,24200,100,30025,7,8,11

21、100,300,200,10014,19,24200,0,20225,7,8,11100,100,100,014,19,24100,0,100由表1、2可知：1）方案1年综合费用Ctotal最小，比方案2少95.75 万元，降低了7.05%，比方案3少162.66 万元，降低了11.4%。说明基于主动管理模式和需求侧管理的ADN规划方法的规划结果整体较优。2）方案1线路投资费用Cline分别比方案2、3少80.98、109.58 万元，这说明考虑需求侧响应和对DG采用主动管理手段可以减少和延缓电网固定资产投资，为配电公司创造效益。3）方案2的DG安装和维护费用Cdg比方案3高16.65%，装机容量比方案3高600 kW，说明采取主动管理措施提高了配电网对分布式电源的接纳能力，方案1的Cdg比方案2高21.75%，即