CN110956329B - 基于分布式光伏与电动汽车时空分布的负荷预测方法 - Google Patents
基于分布式光伏与电动汽车时空分布的负荷预测方法 Download PDFInfo
- Publication number
- CN110956329B CN110956329B CN201911213446.1A CN201911213446A CN110956329B CN 110956329 B CN110956329 B CN 110956329B CN 201911213446 A CN201911213446 A CN 201911213446A CN 110956329 B CN110956329 B CN 110956329B
- Authority
- CN
- China
- Prior art keywords
- load
- time
- photovoltaic
- formula
- membership
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000011156 evaluation Methods 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 12
- 230000007704 transition Effects 0.000 claims abstract description 7
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 230000005611 electricity Effects 0.000 claims description 9
- 230000001174 ascending effect Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 238000012549 training Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 239000002245 particle Substances 0.000 claims description 2
- 238000011161 development Methods 0.000 description 11
- 238000012546 transfer Methods 0.000 description 5
- 230000006399 behavior Effects 0.000 description 4
- 238000012886 linear function Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000005286 illumination Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000004378 air conditioning Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种基于分布式光伏与电动汽车时空分布的负荷预测方法,该方法步骤为:S1:分别建立光伏时空负荷预测模型和电动汽车充电负荷时空预测模型;S2:采用层次分析法确定各项影响指标之间的权重、采用模糊综合评价法进行各项影响指标的综合评价,确定各项影响指标下的综合配比;S3:分布式光伏时间出力预测模型由LS‑SVM模型的预测结果进行相应规划区的配比计算即可得到各规划区光伏出力大小;S4:充电负荷时空预测模型由基于出行时间概率转移矩阵下的蒙特卡洛算法预测得出;S5:将上述预测结果与传统电力负荷进行叠加,获得城市配电网综合负荷的时空分布情况。本发明得到的空间负荷预测值能为城市配电网规划提供重要参考。
Description
技术领域
本发明属于城市配电网空间负荷预测技术领域,具体地说是一种基于分布式光伏与电动汽车时空分布的负荷预测方法。
背景技术
随着城市配电网规划的完善以及智能电网建设的不断推进,大量分布式能源、分布式储能、智能用电设备等多元主体不断接入配电侧电力结构中,给配电网空间负荷预测带来新形势下的严峻挑战。
根据城市空间发展规划,屋顶分布式光伏与公共建筑、居民住宅以及工业园区等基础设施天然契合,具有巨大的发展潜力;同时,电动汽车续航能力的提升与充电设施的不断完善极大地促进了电动汽车在城市的保有量。因此,大规模屋顶分布式光伏电源与电动汽车充电负荷的高比例渗透影响了城市电网源荷时空分布特性:时间分布上,光伏出力受天气、光照等因素影响,用户的出行习惯与电动汽车充电时间密切相关;空间分布上,不同规划区内建筑结构与地理位置的差异将造成屋顶光伏在空间出力分布上的随机性,而用户出行目的地的随机转移将影响各功能区的充电负荷差异。因此,如何在计及光伏和电动汽车时空随机分布的情况下实现有效的配电网空间负荷预测,是一个亟待解决的问题。
对于配电网空间负荷预测的研究,不仅局限于用地仿真法、负荷密度指标法、智能算法下传统用电负荷的空间预测,还基于交通出行矩阵、停车生成率等进行电动汽车空间负荷预测。对于光伏空间负荷预测的研究大多涉及光伏电站时空相关性下的并网规划。一部分文献结合光伏与电动汽车对负荷影响的历年数据,采用LS-SVM算法进行空间负荷预测。一部分考虑屋顶光伏安装结构与电动汽车停车需求概率预测远景年的空间负荷。然而,上述文献大多侧重于传统负荷在空间分布上的特性,所涉及的新能源规划大多只包含单一的光伏或电动汽车,对于充分挖掘光伏与电动汽车时空分布的具体影响因素与特性涉及较少。
发明内容
本发明的目的是针对现有技术存在的问题,提供一种基于分布式光伏与电动汽车时空分布的负荷预测方法;该方法在在考量屋顶分布式光伏与电动汽车的时空分布特性的基础上,为城市配电网的空间负荷预测提供新形势下的合理的参考。
本发明的目的是通过以下技术方案解决的:
一种基于分布式光伏与电动汽车时空分布的负荷预测方法,
S1:屋顶分布式光伏、电动汽车充电负荷和传统电力负荷在时空分布上具有一定的独立性,故分别建立光伏时空负荷预测模型和电动汽车充电负荷时空预测模型;
S2:光伏在空间上的出力主要以屋顶光伏的形式呈现,而不同规划区用地的建筑结构、地理位置等多项因素会对屋顶光伏出力造成很大程度的影响,故采用层次分析法确定对屋顶分布式光伏出力造成影响的各项影响指标之间的权重、采用模糊综合评价法进行各项影响指标的综合评价,最终确定各项影响指标下的综合配比;
S3:分布式光伏时间出力预测模型由LS-SVM模型的预测结果进行相应规划区的配比计算即可得到各规划区光伏出力大小;
S4:充电负荷时空预测模型由基于出行时间概率转移矩阵下的蒙特卡洛算法预测得出;
S5:综合光伏时空负荷预测模型和电动汽车充电负荷时空预测模型预测得到的结果,与传统电力负荷进行叠加,即可得城市配电网综合负荷的时空分布情况。
进一步,所述步骤S2光伏空间出力配比模型的构建中,对于层次分析法评价因素的确定与划分,本发明将城市规划区划分为居民住宅区、商业区和工业区,选取主要七项指标建立评价因素集U={u1,u2,u3,u4,u5,u6,u7},如表1所示:
表1评价因素集
在考量上述各项指标的基础之上构成配比评语集V={v1,v2,v3},即各规划区的配比评语越高,光伏出力的占比则越大,如表2所示:
表2评语集
进一步,所述步骤S2光伏空间出力配比模型的构建中,对于模糊综合评价法中单项因素隶属函数的建立,由于评语集以各规划区的分类来衡量配比权重,故针对各项因素指标采用统一的隶属度函数来评判各规划区相应指标的优劣。
屋顶可利用率隶属度函数如式(1)所示,城市建筑屋顶主要分为水平和倾斜两种,大多为南北朝向,而带东/北坡的倾斜屋顶不适合安装光伏面板,且不同规划区的屋顶管道、遮阳面、空调架等存在不同程度的影响,故屋顶可利用率是影响光伏出力的因素指标之一。其隶属度函数选用升半正态分布,隶属度越高重要性越大。
式(1)中,u1表示屋顶可利用率,C1(u1)表示相应的隶属度。
白天传统电力负荷密度的隶属度函数如式(2)所示。白天传统电力负荷密度越高意味着该规划区在白天的电力需求越密集,相应地可配一定程度的屋顶光伏出力来支撑传统电力负荷,故其隶属度函数选用升半正态分布。
式(2)中,u2表示白天传统电力负荷密度且单位为W/m2,C2(u2)表示相应的隶属度。
高层建筑占比的隶属度函数如式(3)所示。光伏组件单体面积大,屋顶光伏面板所处的建筑物高度越高,其风荷载越大,安装与维护的成本越高,不利于屋顶光伏的稳定发展,故规划区内高层建筑比的隶属度函数取降半正态分布,隶属度越高重要性越小。
式(3),u3表示高层建筑占比,C3(u3)表示相应的隶属度。
光伏上网电价的隶属度函数如式(4)所示。高电价在一定程度上能促进屋顶分布式光伏的发展,但超过一定幅值的电价可能对用户的积极性和普及率造成影响,阻碍屋顶光伏的发展,故光伏上网电价的隶属度函数取正态分布。
式(4),u4表示光伏上网电价且单位为元/kWh,C4(u4)表示相应的隶属度。
光伏财政补贴的隶属度函数如式(5)所示。各规划区不同的财政补贴政策有利于屋顶光伏的发展,而越高补贴下的提升,其发展的增长率会变缓,故光伏财政补贴隶属度函数取升半柯西分布。
式(5),u5表示光伏财政补贴且单位为万元,C5(u5)表示相应的隶属度。
热岛效应温度的隶属度函数如式(6)所示。随着城市化的不断发展,各规划区之间的热岛效应越来越明显,以不同程度的温度差异影响着屋顶光伏的发展。由于平均温度与光伏出力具有一次函数相关性,故其隶属度函数取升一次函数分布。
式(6)中,u6表示热岛效应温度且单位为℃,C6(u6)表示相应的隶属度。
交通流量等级隶属度分析如式(7)所示。交通流量等级的区别在一定程度上反映了不同规划区的经济流动情况,屋顶光伏的发展与等级之间存在一定的正比关系,故其隶属函数取降一次函数分布(1~10等级):
式(7)中,u7表示交通流量等级,C7(u7)表示相应的隶属度;
确定各因素指标对评语集的隶属度函数后,可形成隶属度函数矩阵R,如下所示:
式(8)中,上标H、B和I分别表示各因素指标在住宅区、商业区和工业区的隶属度。
进一步,所述步骤S2的光伏空间出力配比模型的构建中,对于模糊综合评价法中各因素指标权重的确定,由于各因素指标对评语集重要程度的关系复杂,难以直观地给出各项权重,故引入层次分析法中矩阵判断标度,通过两两因素指标之间的重要性来构成相关联的有序层次,减少一定程度的主观判断。矩阵判断标度如表3所示。
表3矩阵判断标度
进一步,所述步骤S2光伏空间出力配比模型的构建中,对于模糊综合评价集,为兼顾各规划区的配比,较好地依权重大小进行评价,本发明采用加权平均型模糊合成算子并对其进行归一化处理,计算总评价集B,如式(9)所示
式(9)中,Bk为第k个规划区的配比;El为第l个因素指标所占的权重;Rl,k为第l个因素指标对第k个规划区的评价隶属度。
进一步,所述步骤S3光伏时间出力预测模型的构建中,对于训练样本的选取与处理,由于光伏出力受季节和天气的影响,故选取同一季节同一天气情况下的数据作为历史数据样本,包括一天的平均温度、平均湿度、最高温度、最低温度和各时段的光伏出力值。为充分利用各样本间的相关性,选取以下特性指标组成模型训练样本:
1)前两天的样本日同一预测点负荷、平均湿度、最高温度、最低温度、平均温度。
2)前一天的样本日同一预测点负荷、平均湿度、最高温度、最低温度、平均温度。
3)预测日的平均湿度、最高温度、最低温度、平均温度。
其中,温度、湿度数据进行归一化处理,负荷数据进行对数归一化处理。
对于核函数的确定,径向基函数将样本数据非线性变换到高维空间中,能够处理输入、输出为非线性关系的情况,且表示形式简单,参数设置少,故采用径向基函数作为回归模型中的核函数:
对于模型参数的确定,LS-SVM(least squares support vectormachine,最小二乘支持向量机)模型在建模过程中要确定2个重要参数:正则化参数γ和径向基核函数宽度γ作为平衡因子,控制优化目标超出误差的惩罚程度;影响模型的预测性能。故本发明采用粒子群算法自动迭代优化选取γ和提高预测精度。适应度函数为预测结果的平均绝对百分误差比:
基于所建立的LS-SVM预测模型,整体的光伏时间出力预测流程如图2所示。
进一步,所述步骤S4电动汽车空间负荷预测模型的构建中,对于电动汽车时空分布的描述,电动汽车充电负荷与一般性负荷的区别在于其随机流动性所带来的负荷转移,是一个随机离散的过程,并具有一定程度的概率特性,故参考马尔科夫链的状态转移过程来描述电动汽车的时空分布特性。
电动汽车的出行与用户的目的地密切相关,依据上文对城市规划区的划分,任意一个时刻t,电动汽车共有四种状态:停在住宅区、停在商业区、停在工业区和行驶中,且不同状态转移之间存在一定的概率i,j∈{H,B,I,D}。如图3所示。H、B、I、D分别表示住宅区、商业区、工业区和行驶状态,ij任意二者组合表示由i→j的状态转移。
则任意一个时刻t,电动汽车从一个状态到另一个状态的转移概率矩阵为:
进一步,所述步骤S4电动汽车空间负荷预测模型的构建中,对于电动汽车充电负荷预测的方法,考虑凌晨时刻0的起始状态,电动汽车89%停在住宅区,4%停在商业区,2%停在工业区,5%处于行驶状态;考虑日行驶距离模型,私家车日行驶里程近似服从对数正态分布,故每时段下的行驶里程概率密度函数为:
式(13)中,dΔ为每段行驶里程且单位为km;μD为日行驶里程均值,取3.7km;σD为行驶里程标准差,取0.9。
考虑充电决策的设置,电动汽车用户多数情况下一天一充即可保证当日的出行需求,同时为考虑到电池寿命与低电量出行保护,电池SOC(state of charge,SOC)低于20%时需充电,其满足的条件如下:
式(14)中,Sini为初始SOC值;dΔ(t)为t时段内行驶的里程数且单位为km;Tch表示需充电时的时段;Dmax为最大行驶里程且单位为km。充电时满足的条件如下:
式(15)中,S(t)为时段t结束时刻的SOC值;pch为充电功率;Δt为一个负荷预测时段的长度;B为电池容量;Sexp为期望充到的SOC值。
时空分布下电动汽车充电负荷预测的整体流程基于蒙特卡洛的概率随机模拟的方法进行时空充电负荷的预测,如图4所示
本发明相比现有技术有如下优点:
本发明建立了完备的屋顶光伏空间出力评价体系,在充分考量各项因素指标的基础上实现了合理的算例分析;建立了光伏时间出力预测模型,在一定的数据背景下,可实现了较精确的预测结果。
本发明建立了电动汽车空间负荷预测模型,充分考虑了电动汽车时空转移特性,实现了合理的充电负荷预测。
本发明通过叠加整合得到不同规划区整体的空间负荷预测结果,可为城市配电网的合理规划提供一定的参考。
附图说明
附图1为本发明的方法流程图;
附图2为光伏时间出力预测流程图;
附图3为电动汽车时空转移分布图;
附图4为时空分布下电动汽车充电负荷预测流程图;
附图5为各因素指标权重图;
附图6为工业区夏季晴天光伏出力预测结果图;
附图7为夏季晴天住宅区的配电网日空间负荷预测结果图;
附图8为夏季晴天商业区的配电网日空间负荷预测结果图;
附图9为夏季晴天工业区的配电网日空间负荷预测结果图。
其中:
具体实施方式
下面结合附图与实施例对本发明作进一步的说明。
步骤1:
1)根据某地区实际经验,考察评价各因素指标之间的重要性关系,确定各因素指标判断矩阵,如表4所示。
表4各因素指标判断矩阵
2)根据判断标度矩阵,采用层次分析法可计算出各因素指标之间的权重E,如图5所示。
3)某城区各规划区之间的因素指标实际数据如表5所示。
表5各规划区因素指标实际数据
4)选取该地工业区夏季晴天的屋顶光伏相关数据作为LS-SVM的历史训练样本。
5)假设该地区私家车10万辆,电动汽车渗透率为5%,具体的电动汽车参数如表6所示。初始时刻0时大部分用户在家充电,故假设住宅区0时已充至90%,工商业已充至70%,且住宅区充电站为慢充,工商业为快充,其余参数选取一致。
表6电动汽车相关参数
步骤2:将表5数据逐一代入各因素指标的隶属度函数中,得到的隶属度矩阵如式(16)所示。
根据说明书中式(9)可计算出总评价集B=[0.2019 0.3338 0.4643],分别表示住宅区、商业区和工业区的光伏出力配比。
步骤3:基于上述数据,工业区夏季晴天光伏时间出力预测结果如图6所示。由图6可知,受限于光照的时间特性,光伏只能在白天出力,午后左右到达峰值,且含有较大的预测误差,其余时段的预测误差不大,整体上平均绝对百分误差比为6.6%,具有较高的预测精度。
步骤4:基于光伏空间配比与工业区光伏出力预测即可预测各规划区的光伏出力情况,并与各规划区电动汽车空间负荷和传统电力负荷叠加,可得到该城在夏季晴天配电网空间负荷预测曲线,如图7-9所示。
由图7-9可知,各规划区屋顶光伏出力与配比模型一一对应,住宅区光伏出力最小,工业区光伏出力最大,并对其白天负荷高峰起到不同程度的削峰作用。不同规划区电动汽车的充电行为呈现一定的差异性,住宅区用户一般在晚上到家后陆续开始充电,负荷逐步到达高峰,且整体负荷量大;商业区与工业的充电行为类似,白天频繁的随机移动带来充电负荷的累积,下午的负荷量高于上午,整体上工业区的负荷量小。用户的充电行为与其生活密切相关,可能与传统用电习惯相吻合,造成“峰上加峰”的情况,如住宅区的夜晚时段和商业区的下午时段,加剧了规划区负荷的峰谷差。
整体上,在计及屋顶分布式光伏出力和电动汽车充电负荷的基础上,各规划区的负荷变化与传统电力空间负荷有较大的差异,主要表现在光伏出力的削峰与电动汽车负荷的增峰引起的负荷波动。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内;本发明未涉及的技术均可通过现有技术加以实现。
Claims (3)
1.一种基于分布式光伏与电动汽车时空分布的负荷预测方法,其特征在于:该方法包括以下步骤:
S1:分别建立光伏时空负荷预测模型和电动汽车充电负荷时空预测模型;
S2:采用层次分析法确定对屋顶分布式光伏出力造成影响的各项影响指标之间的权重、采用模糊综合评价法进行各项影响指标的综合评价,最终确定各项影响指标下的综合配比;
S3:分布式光伏时间出力预测模型由LS-SVM模型的预测结果进行相应规划区的配比计算即可得到各规划区光伏出力大小;
S4:充电负荷时空预测模型由基于出行时间概率转移矩阵下的蒙特卡洛算法预测得出;
S5:综合分布式光伏时空负荷预测模型和电动汽车充电负荷时空预测模型预测得到的结果,与传统电力负荷进行叠加,即可得城市配电网综合负荷的时空分布情况;
所述步骤S2中,在分布式光伏的空间出力配比模型的构建中,对于层次分析法评价因素的确定与划分,将城市规划区划分为居民住宅区v1、商业区v2和工业区v3,采用屋顶可利用率u1、白天传统负荷密度u2、高层建筑比u3、光伏上网电价u4、光伏财政补贴u5、热岛效应u6和交通流量等级u7计七项指标建立评价因素集U={u1,u2,u3,u4,u5,u6,u7},在考量上述七项指标的基础之上构成配比评语集V={v1,v2,v3},即各规划区的配比评语越高,光伏出力的占比则越大;
所述步骤S2中,对于模糊综合评价法中单项因素隶属函数的建立,由于评语集以各规划区的分类来衡量配比权重,故针对各项因素指标采用统一的隶属度函数来评判各规划区相应指标的优劣;
其中屋顶可利用率u1、隶属度函数如式(1)所示,其隶属度函数选用升半正态分布,隶属度越高重要性越大,
式(1)中,u1表示屋顶可利用率,C1(u1)表示相应的隶属度;
白天传统负荷密度u2的隶属度函数如式(2)所示,隶属度函数选用升半正态分布,
式(2)中,u2表示白天传统电力负荷密度且单位为W/m2,C2(u2)表示相应的隶属度;
高层建筑比u3的隶属度函数如式(3)所示,隶属度函数取降半正态分布,隶属度越高重要性越小,
式(3),u3表示高层建筑占比,C3(u3)表示相应的隶属度;
光伏上网电价u4的隶属度函数如式(4)所示,隶属度函数取正态分布,
式(4),u4表示光伏上网电价且单位为元/kWh,C4(u4)表示相应的隶属度;
光伏财政补贴u5的隶属度函数如式(5)所示,隶属度函数取升半柯西分布,
式(5),u5表示光伏财政补贴且单位为万元,C5(u5)表示相应的隶属度;
热岛效应u6的隶属度函数如式(6)所示,其隶属度函数取升一次函数分布,
式(6)中,u6表示热岛效应温度且单位为℃,C6(u6)表示相应的隶属度;
交通流量等级u7的隶属度分析如式(7)所示,其隶属函数取降一次函数分布,
式(7)中,u7表示交通流量等级,C7(u7)表示相应的隶属度;
确定各因素指标对评语集的隶属度函数后,可形成隶属度函数矩阵R,如式(8)所示:
式(8)中,上标H、B和I分别表示各因素指标在住宅区、商业区和工业区的隶属度;
所述步骤S2中,对于模糊综合评价法中各因素指标权重的确定,引入层次分析法中矩阵判断标度,通过两两因素指标之间的重要性来构成相关联的有序层次,减少一定程度的主观判断;所述步骤S2中,对于模糊综合评价集,采用加权平均型模糊合成算子并对其进行归一化处理,计算总评价集B,如式(9)所示,
式(9)中,Bk为第k个规划区的配比;El为第l个因素指标所占的权重;Rl,k为第l个因素指标对第k个规划区的评价隶属度;
所述步骤S4中的电动汽车空间负荷预测模型中,电动汽车的出行与用户的目的地密切相关,任意一个时刻t,电动汽车共有四种状态:停在住宅区、停在商业区、停在工业区和行驶中,且不同状态转移之间存在一定的概率H、B、I、D分别表示住宅区、商业区、工业区和行驶状态,ij任意二者组合表示由i→j的状态转移,则任意一个时刻t,电动汽车从一个状态到另一个状态的转移概率矩阵为:
式(12)中,均以行驶目的地来表示行驶状态之间的连续,则为0,即不存在一个行驶状态到另一个行驶状态的过程;Pt trans每行概率之和为1。
2.根据权利要求1所述的基于分布式光伏与电动汽车时空分布的负荷预测方法,其特征在于:所述步骤S3中,在分布式光伏的时间出力预测模型的构建中,对于训练样本的选取与处理,由于光伏出力受季节和天气的影响,故选取同一季节同一天气情况下的数据作为历史数据样本,包括一天的平均温度、平均湿度、最高温度、最低温度和各时段的光伏出力值;
为充分利用各样本间的相关性,选取以下特性指标组成模型训练样本:1)前两天的样本日同一预测点负荷、平均湿度、最高温度、最低温度、平均温度,2)前一天的样本日同一预测点负荷、平均湿度、最高温度、最低温度、平均温度,3)预测日的平均湿度、最高温度、最低温度、平均温度;其中,温度、湿度数据进行归一化处理,负荷数据进行对数归一化处理;
对于核函数的确定,采用径向基函数作为回归模型中的核函数:
式(10)中,K(x,xi)为核函数,x是m维输入向量;xi是第i个径向基函数的中心;为径向基核函数宽度,决定了该函数围绕中心点的宽度;||x-xi||是向量x-xi的范数,表示x与xi之间的距离;
对于模型参数的确定,正则化参数γ作为平衡因子,控制优化目标超出误差的惩罚程度;影响模型的预测性能;采用粒子群算法自动迭代优化选取γ和提高预测精度;
适应度函数为预测结果的平均绝对百分误差比:
式(11)中,fadp为适应度函数,Tpv为预测的总时段数;Wt为第t个时段内实测的出力值;为第t个时段内预测的出力值。
3.根据权利要求1所述的基于分布式光伏与电动汽车时空分布的负荷预测方法,其特征在于:所述步骤S4中,对于电动汽车充电负荷预测的方法,考虑凌晨时刻0点的起始状态,电动汽车89%停在住宅区,4%停在商业区,2%停在工业区,5%处于行驶状态;考虑日行驶距离模型,私家车日行驶里程近似服从对数正态分布,故每时段下的行驶里程概率密度函数为:
式(13)中,dΔ为每段行驶里程且单位为km;μD为日行驶里程均值且单位为km;σD为行驶里程标准差;
设定电动汽车用户一天一充且电池SOC低于20%时需充电,其满足的条件如下:
式(14)中,Sini为初始SOC值;dΔ(t)为t时段内行驶的里程数且单位为km;Tch表示需充电时的时段;Dmax为最大行驶里程且单位为km;充电时满足的条件如下:
式(15)中,S(t)为时段t结束时刻的SOC值;pch为充电功率;Δt为一个负荷预测时段的长度;B为电池容量;Sexp为期望充到的SOC值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911213446.1A CN110956329B (zh) | 2019-12-02 | 2019-12-02 | 基于分布式光伏与电动汽车时空分布的负荷预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911213446.1A CN110956329B (zh) | 2019-12-02 | 2019-12-02 | 基于分布式光伏与电动汽车时空分布的负荷预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110956329A CN110956329A (zh) | 2020-04-03 |
CN110956329B true CN110956329B (zh) | 2023-04-25 |
Family
ID=69979306
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911213446.1A Active CN110956329B (zh) | 2019-12-02 | 2019-12-02 | 基于分布式光伏与电动汽车时空分布的负荷预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110956329B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112257907B (zh) * | 2020-09-30 | 2023-04-07 | 广东电网有限责任公司广州供电局 | 一种考虑电价敏感性的电动汽车负荷建模方法 |
CN111967804B (zh) * | 2020-10-23 | 2021-01-29 | 北京国新智电新能源科技有限责任公司 | 基于移动充电设备主动配电的调度*** |
CN113325701B (zh) * | 2021-04-21 | 2022-06-07 | 西安交通大学 | 含新能源的电动车与建筑空调***协同优化方法与装置 |
CN113642787A (zh) * | 2021-08-04 | 2021-11-12 | 国网江苏省电力有限公司南通供电分公司 | 基于配售电竞争态势下区域配电网网格化负荷区间预测方法 |
CN113971348B (zh) * | 2021-11-11 | 2024-04-02 | 国网山东省电力公司日照供电公司 | 一种基于综合能源***的源端功率短期预测方法 |
CN114465256B (zh) * | 2022-01-24 | 2024-06-21 | 东北电力大学 | 多节点电动汽车充电负荷联合对抗生成区间预测方法 |
CN114498632A (zh) * | 2022-01-27 | 2022-05-13 | 国网江苏省电力有限公司泰州供电分公司 | 基于新能源和充电设施的配电台区负荷预测方法及装置 |
CN114548171A (zh) * | 2022-02-22 | 2022-05-27 | 国网江苏省电力有限公司南通供电分公司 | 一种基于盲信号分离的配电台区电动汽车负荷识别方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514487A (zh) * | 2013-07-15 | 2014-01-15 | 国家电网公司 | 一种含分布式电源的配电网负荷预测方法 |
CN104993522A (zh) * | 2015-06-30 | 2015-10-21 | 中国电力科学研究院 | 一种基于mpc的主动配电网多时间尺度协调优化调度方法 |
CN107392400A (zh) * | 2017-09-04 | 2017-11-24 | 重庆大学 | 计及实时交通与温度的ev充电负荷时空分布预测方法 |
CN108062591A (zh) * | 2016-11-05 | 2018-05-22 | 南京理工大学 | 电动汽车充电负荷时空分布预测方法 |
CN108288112A (zh) * | 2018-01-30 | 2018-07-17 | 广西电网有限责任公司柳州供电局 | 基于用户出行模拟的区域电动汽车充电站负荷预测方法 |
CN108288231A (zh) * | 2018-01-19 | 2018-07-17 | 广东电网有限责任公司河源供电局 | 一种分布式光伏接入对配电台区负荷特性影响的评估方法 |
CN108510128A (zh) * | 2018-04-11 | 2018-09-07 | 华南理工大学广州学院 | 一种区域电动汽车充电负荷时空分布预测方法 |
CN109193636A (zh) * | 2018-10-08 | 2019-01-11 | 华东交通大学 | 一种基于分类不确定集的电力***环境经济鲁棒调度方法 |
CN109871981A (zh) * | 2019-01-10 | 2019-06-11 | 国家电网有限公司 | 一种计及分布式电源和电动汽车影响的负荷特性预测方法 |
CN109934403A (zh) * | 2019-03-12 | 2019-06-25 | 浙江大学 | 基于数学模型的电动汽车居民区域充电负荷分析预测方法 |
CN110210647A (zh) * | 2019-04-29 | 2019-09-06 | 国网江苏省电力有限公司电力科学研究院 | 一种分布式电源、储能与柔性负荷联合调度方法及装置 |
CN110350535A (zh) * | 2018-09-07 | 2019-10-18 | 国网甘肃省电力公司电力科学研究院 | 一种计及大规模风电和电动汽车并网的配电网电压稳定概率评估方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012162646A1 (en) * | 2011-05-26 | 2012-11-29 | Ice Energy, Inc. | System and method for improving grid efficiency utilizing statistical distribution control |
US10693294B2 (en) * | 2012-09-26 | 2020-06-23 | Stem, Inc. | System for optimizing the charging of electric vehicles using networked distributed energy storage systems |
-
2019
- 2019-12-02 CN CN201911213446.1A patent/CN110956329B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514487A (zh) * | 2013-07-15 | 2014-01-15 | 国家电网公司 | 一种含分布式电源的配电网负荷预测方法 |
CN104993522A (zh) * | 2015-06-30 | 2015-10-21 | 中国电力科学研究院 | 一种基于mpc的主动配电网多时间尺度协调优化调度方法 |
CN108062591A (zh) * | 2016-11-05 | 2018-05-22 | 南京理工大学 | 电动汽车充电负荷时空分布预测方法 |
CN107392400A (zh) * | 2017-09-04 | 2017-11-24 | 重庆大学 | 计及实时交通与温度的ev充电负荷时空分布预测方法 |
CN108288231A (zh) * | 2018-01-19 | 2018-07-17 | 广东电网有限责任公司河源供电局 | 一种分布式光伏接入对配电台区负荷特性影响的评估方法 |
CN108288112A (zh) * | 2018-01-30 | 2018-07-17 | 广西电网有限责任公司柳州供电局 | 基于用户出行模拟的区域电动汽车充电站负荷预测方法 |
CN108510128A (zh) * | 2018-04-11 | 2018-09-07 | 华南理工大学广州学院 | 一种区域电动汽车充电负荷时空分布预测方法 |
CN110350535A (zh) * | 2018-09-07 | 2019-10-18 | 国网甘肃省电力公司电力科学研究院 | 一种计及大规模风电和电动汽车并网的配电网电压稳定概率评估方法 |
CN109193636A (zh) * | 2018-10-08 | 2019-01-11 | 华东交通大学 | 一种基于分类不确定集的电力***环境经济鲁棒调度方法 |
CN109871981A (zh) * | 2019-01-10 | 2019-06-11 | 国家电网有限公司 | 一种计及分布式电源和电动汽车影响的负荷特性预测方法 |
CN109934403A (zh) * | 2019-03-12 | 2019-06-25 | 浙江大学 | 基于数学模型的电动汽车居民区域充电负荷分析预测方法 |
CN110210647A (zh) * | 2019-04-29 | 2019-09-06 | 国网江苏省电力有限公司电力科学研究院 | 一种分布式电源、储能与柔性负荷联合调度方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110956329A (zh) | 2020-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110956329B (zh) | 基于分布式光伏与电动汽车时空分布的负荷预测方法 | |
Zhang et al. | A review of urban energy systems at building cluster level incorporating renewable-energy-source (RES) envelope solutions | |
Jin et al. | Hierarchical microgrid energy management in an office building | |
Shepero et al. | Modeling of photovoltaic power generation and electric vehicles charging on city-scale: A review | |
Tian et al. | A morphology-based evaluation on block-scale solar potential for residential area in central China | |
Jiang et al. | Comprehensive assessment of wind resources and the low-carbon economy: An empirical study in the Alxa and Xilin Gol Leagues of inner Mongolia, China | |
CN110110912A (zh) | 一种光伏功率多模型区间预测方法 | |
Zhang et al. | Data-driven cooling, heating and electrical load prediction for building integrated with electric vehicles considering occupant travel behavior | |
Zhu et al. | Optimization of photovoltaic provision in a three-dimensional city using real-time electricity demand | |
CN102769298B (zh) | 一种太阳能并网发电功率预测方法及预测*** | |
Dewi et al. | Equitable, affordable, and deep decarbonization pathways for low-latitude developing cities by rooftop photovoltaics integrated with electric vehicles | |
CN117407939A (zh) | 一种基于数字孪生的光伏***设计模型***及其使用方法 | |
Rezaei Nasab et al. | Investigating the probability of designing net-zero energy buildings with consideration of electric vehicles and renewable energy | |
Díaz et al. | Key factors influencing urban wind energy: A case study from the Dominican Republic | |
Hu et al. | Exploring cost-effective strategies for emission reduction of public buildings in a life-cycle | |
CN114498632A (zh) | 基于新能源和充电设施的配电台区负荷预测方法及装置 | |
CN117094521A (zh) | 一种光伏电站智能选址方法 | |
CN115018376B (zh) | 一种考虑新型电力***特征的负荷调控优化方法 | |
CN113052450B (zh) | 一种适应电能替代发展战略的城市能源互联网规划方法 | |
CN115545429A (zh) | 基于cpss架构的多场景低压可再生资源负荷预测方法 | |
CN111191954A (zh) | 一种综合能源技术适应性评价方法 | |
Yuying et al. | Spatial Load Forecasting of Distribution Network Based on Artificial Intelligent Algorithm | |
Zhang et al. | Techno-environmental-economical performance of allocating multiple energy storage resources for multi-scale and multi-type urban forms towards low carbon district | |
Belloni et al. | A street lighting based on solar energy: optimization of a neural network for the productivity evaluation in different climate conditions | |
Feng et al. | Low Carbon Building Design Optimization Based on Intelligent Energy Management System |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |