[发明专利]基于负荷预测和场景聚类的城市配电网规划方法

权利要求书

1.一种基于负荷预测和场景聚类的城市配电网规划方法，所述城市配电网中包含N个负荷节点，将所述负荷节点记为{L₁,L₂,…,L_i,…,L_N}，L_i表示第i个负荷节点，i＝1,2,…,N，其特征在于，所述城市配电网规划方法是按照以下步骤进行：

步骤1：采集城市历史H年的年总用电量数据并按照时间先后排序得到总用电量数据序列{x₁,x₂,…,x_h,…,x_H}，其中x_h代表城市第H-h年前的年总用电量负荷；h＝1,2,…,H；

步骤2：设置采样周期，并采集所述第i个负荷节点L_i在当前年份的第a个采样时刻的负荷数据P_i-a，从而得到第i个负荷节点L_i当前年份各采样时刻的负荷数据P_i＝{P_i-1,P_i-2…,P_i-a,…,P_i-A}，进而得到所有负荷节点在当前年份的各采样时刻的负荷数据P＝{P₁,P₂…,P_i,…,P_N}；a＝1,2,…,A，A对应一年内的总采样时刻；

步骤3：应用回归拟合法得到所述城市配电网未来F年的年总用电量预测序列表示应用回归拟合函数预测的第f个序列值，代表未来f年后的年总用电量的预测值；f＝1,2,…,F；

步骤4：利用式(1)得到第i个负荷节点L_i在未来第f年第a个采样时刻的负荷预测数据

从而得到第i个负荷节点L_i在未来第f年对应所有采样时刻的负荷数据进而得到城市所有负荷节点在未来第f年对应所有采样时刻的负荷预测数据

步骤5：利用k均值聚类算法对城市所有负荷节点在未来第f年对应所有采样时刻的负荷预测数据P^f进行聚类分析，得到聚类结果；

步骤5-1：基于所述城市所有负荷节点未来第f年对应所采样时刻的负荷预测数据，计算城市在未来第f年第a个采样时刻的总负荷预测数据进而获得城市未来A个采样时刻的总负荷预测数据并选择A个总负荷预测数据中的最大值与最小值分别记为与根据城市负荷分布特性，确定初始分类数M^f；并将f初值赋值为1；

步骤5-2：以步长将所述第i个负荷节点L_i的未来第f年城市负荷预测数据P_i^f分为M^f类，从而得到M^f个区间分组表示未来第f年的第m^f类区间分组；统计所述区间分组内的项数，记为

步骤5-3：定义为对应区间分组的第i个负荷节点L_i在未来第f年城市负荷预测数据的样本均值，并作为第i个负荷节点L_i的第m^f类区间分组的分类中心，从而获得第m^f类区间分组的均值初始向量

步骤5-4：采用式(2)计算所有负荷节点未来第f年的各采样时刻的负荷数据预测结果到所述第m^f类区间分组的分类中心的距离

步骤5-5：以距离最小为原则，调整未来第f年城市负荷预测数据的分类，获得调整后的M^f个区间分组重新统计调整后的各区间分组的项数重新计算调整后的第m^f类区间分组的类均值向量采用式(3)计算第i个负荷节点L_i的总均值从而得到总均值向量

步骤5-6：采用式(4)计算调整后的M^f类分组的类内离差矩阵W^f的元素以及调整后的M^f类分组的类间离差矩阵B^f的元素从而得到调整后的M^f类分组的类内离差矩阵W^f以及调整后的M^f类分组的类间离差矩阵B^f，并根据T^f＝W^f+B^f计算总离差矩阵T^f，i,j＝1,2,…,N；

步骤5-7：拟定置信区间下的χ²临界值，判断χ²检验的计算结果是否大于自由度n(M^f-1)，如果大于，则认为调整后的分组有效，将调整后的区间分组的类均值向量的各分类中心作为各负荷节点对应区间分组的典型负荷，进一步得到所有区间分组中各负荷节点的典型负荷；否则，认为调整后的分组无效，并将M^f+1赋值给M^f后，并返回步骤5-2顺序执行，重新聚类；

步骤5-8：将f+1赋值给f后，返回步骤5-1，从而依序对未来f年负荷预测数据进行聚类分析，直至f＝F；

步骤6：根据规划目标与决策变量列出城市配电网规划的目标函数，所述目标函数包括：城市规划期间电费收入C_Rev和网络综合费用；

所述网络综合费用包括：城市规划期间配电网投资建设成本C_inv、配电网系统运行成本C_ope、配电网系统维护成本C_Mn、配电网故障损失成本C_int和城市配电网评价指标对应的惩罚成本C_ind；

步骤6-1：以规划期间内城市配电网的净现值最大为目标函数，所述净现值为城市配电网电费收益与网络综合费用之差；并利用式(5)得到净现值C：

C＝C_Rev-(C_inv+C_ope+C_Mn+C_int+C_ind) (5)

步骤6-2：以时间T年为间隔，将城市配电网建设方案划分为K个阶段；

假设未来第f_k,t年表示城市配电网建设第k阶段的第t年，则所述城市配电网建设方案的第k阶段所包含的所有年份记为{f_k,1,f_k,2,…,f_k,t,…,f_k,T}；k＝1,2,…,K；t＝1,2,…,T；初始化k＝1；

步骤6-3：针对第k阶段配电网建设方案，统计城市第k阶段的线路建设集合Ψ_L(k)、线路类型建设集合Ψ′_L、变电站建设集合Ψ_B(k)和变电站类型建设集合Ψ′_B，从而利用式(6)计算第k阶段城市配电网投资建设成本

式(6)中，表示是否建设第s′类线路；若表示不建设第s′类线路，若表示建设第s′类线路，表示是否建设第m′类变电站，若表示不建设第m′类变电站，表示建设第m′类变电站；Len_s代表第s条待建线路的长度，Cap_m代表第m个待建变电站的容量；代表建设单位长度第s′类线路的成本；代表建立单位容量第m′类变电站的成本；

步骤6-4：根据城市未来第f_k,t年各负荷节点负荷预测结果，基于对应第类区间分组的典型负荷利用式(7)计算城市第k阶段电费收入

式(7)中，为未来第f_k,t年第类区间分组的项数，为第i个负荷节点L_i的负荷数据为时每千瓦小时电量买卖收益；t′为采样时间间隔；

步骤6-5：假设城市第k阶段配电网建设将在未来第f_k,T年完成，将城市未来第f_k-1,T年的配电网结构作为城市第k阶段配电网结构参考，确定城市第k阶段配电网结构的导纳矩阵Y^k，统计城市第k阶段配电网系统输电线路表示城市第k阶段配电网中第s_Lk条输电线路，s_Lk＝1,2,…,N_Lk；N_Lk表示城市第k阶段配电网所有输电线路；

统计所述城市第k阶段各变电站容量表示城市第k阶段配电网中第m_Bk个变电站；m_Bk＝1,2,…,N_Bk；N_Bk表示城市第k阶段配电网所有变电站；

利用式(8)所示的潮流计算方法计算城市未来第f_k,t年对应第类区间分组的各负荷节点电压幅值与负荷节点电压相角

式(8)中，与分别为导纳矩阵Y^k对应位置元素的实部与虚部；为所述城市未来第f_k,t年对应区间分组的各负荷节点的无功注入功率；

利用式(9)计算城市未来第f_k,t年对应区间分组的从第i个负荷节点经第s_Lk条输电线路流向第j个负荷节点的传输功率

式(9)中，i,j分别为连接第s_Lk条输电线路的两个负荷节点的编号，第i个负荷节点为第s_Lk条输电线路的首端节点；与为导纳矩阵Y^k对应位置的元素的共轭复数；

利用式(10)计算未来第f_kt年第s_Lk条输电线路对应第类区间分组的线路损耗功率

式(10)中，为城市未来第f_k,t年对应区间分组的从第j个负荷节点经第s_Lk条输电线路流向第i个负荷节点的传输功率；

根据式(11)确定城市第k阶段城市配电网系统线路损耗成本

式(11)中，λ_out为每千瓦小时电量损耗对应的金钱损失系数；为城市未来第f_k,t年所有线路的网损电量，并有：

利用式(12)计算城市第k阶段城市配电网系统变电站运行成本

式(13)中，λ_sub为发电站运行成本系数；

根据计算城市第k阶段城市配电网系统运行成本

步骤6-6：利用式(14)计算所述城市第k阶段配电网维护成本

式(14)中，λ_Mn为配电网维护成本系数；为第k′阶段的配电网投资成本；k′＝1,2,…,k；

步骤6-7：根据所述城市第k阶段配电网结构，设定所述城市第k阶段配电网中第s_Lk条输电线路的平均故障时间进一步得到所述城市第k阶段配电网中所有线路每年的平均故障时间根据城市第k阶段负荷预测数据及其聚类结果，确定在未来第f_kt年对应区间分组的N_Lk条输电线路的城市负荷供电功率不足期望值进一步获得未来第f_k,t年N_Lk条输电线路对应各区间分组城市负荷供电功率不足期望值；根据式(15)计算第k阶段所述城市故障损失费用

式(15)中，λ_int为所述城市在第k阶段规划过程中每Kwh电量不足对应的经济损失成本；为在未来第f_k,t年对应区间分组的第s_Lk条输电线路的城市负荷供电功率不足期望值；

步骤6-8：根据城市第k阶段城市配电网结构及未来第f_k,t年各负荷节点负荷预测结果，建立城市配电网评价指标体系；

1)线损率

根据式(16)计算未来第f_k,t年城市配电网线损率

2)轻载率、重载率、过载率

根据式(17)计算输电线路在未来第f_k,t年对应区间分组的负载率进一步计算所有线路在未来第f_k,t年区间分组的负载率

式(17)中，为未来第f_k,t年第s_Lk条输电线路的传输功率，为未来第f_k,t年第s_Lk条输电线路的额定功率；

将负载率小于20％的线路视为轻载线路，将负载率小于100％并超过80％的线路视为重载线路，将负载率大于100％的线路视为过载线路；

统计未来第f_k,t年城市对应区间分组的轻载、重载与过载线路数目，并分别记为

根据式(18)计算未来第f_k,t年城市对应区间分组的线路轻载率重载率与过载率

利用式(19)进一步得到第f_k,t年城市各线路的轻载率重载率与过载率

3)电压合格率

根据式(20)确定第f_k,t年城市电压偏差

式(20)中，V_R,i为所述城市第i个负荷节点的额定电压；将电压偏差小于额定电压的7％的负荷节点视为电压合格的负荷节点，统计未来第f_k,t年城市对应区间分组电压合格的负荷节点数

根据式(21)计算未来第f_kt年城市的电压合格率

4)N-1通过率

根据所述城市未来第f_k年的城市配电网结构及其负荷预测数据，确定未来第f_k年城市最大负荷运行方式；

在最大负荷运行方式下，变电站出现开关停运后，将全部负荷能在两次操作以内转移到其他线路供电的输电线路视为N-1合格线路，并统计N-1合格线路数目从而根据式(22)算未来第f_k年城市的N-1通过率

根据各评价指标对电力系统的重要程度确定各指标的权重，利用式(23)计算城市第k阶段配电网系统评价指标对应的惩罚成本

式(23)中，分别对应线损率、轻载率，重载率，过载率，电压合格率与N-1通过率的惩罚成本系数；

步骤7：根据城市第k阶段城市投资建设情况，设置约束项包括：1)配电网建设逻辑约束；2)配电网建设约束；3)系统供电不足期望值小于最大值约束；并根据当前阶段系统投资建设方案，调整下一阶段城市各负荷节点所有采样时刻的负荷预测数据；

步骤7-1：利用式(24)构建配电网建设逻辑约束：

步骤7-2：利用式(25)—式(27)构建配电网建设约束：

Contin＝{1，0} (27)

式(25)表示第k阶段建设过程中，每个待选线路位置只安装一种类型的线路；

式(26)表示第k阶段建设过程中，每个待新建变电站位置只建设一种类型的变电站；

式(27)表示网络连通性约束；若Contin＝1代表系统线路是联通的；若Contin＝0代表系统线路是不连通的；

步骤7-3：利用式(28)构建系统供电不足期望值小于最大值约束：

式(28)中，AENS为系统供电不足期望值的最大值；

步骤7-4：根据所述城市第k阶段城市投资建设情况，统计第k阶段与各负荷节点相连的线路和变电站建设情况，并根据式(29)计算第k+1阶段的第i个负荷节点对应区间分组的负荷预测数据进一步得到第k+1阶段的所有负荷节点对应区间分组的负荷预测数据

式(29)中，ψ′_L,i为第k阶段与第i个负荷节点L_i相连的新建线路集合，ψ′_B,i为第k阶段第i个负荷节点L_i供电的新建变电站集合；λ_L为反映线路投资对其相连负荷节点负荷影响的系数；λ_B为反映变电站投资对供电范围内负荷节点负荷影响的系数；

步骤7-5：将的值赋给k+1赋值给k后，返回步骤6-3执行直至k＝K，从而获得城市在整个配电网规划阶段内各规划阶段的配电网净现值{C¹,C²,…,C^k,…,C^K}；

步骤7-6：根据式(30)计算城市整个规划过程中的净现值C：

式(31)中，r_ann代表年利率；

步骤8：对所述城市配电网规划的目标函数进行建模，并利用内点法对模型进行求解，得到最优投资建设方案。