电动汽车充电负荷预测这两年一直是电网规划里的热门话题我最早接触这个方向是因为一个配电网扩容项目需要估算小区层面的负荷峰值。当时手头没有实测充电数据只有车辆保有量和用户出行统计而蒙特卡洛方法恰好能在缺乏实测数据的情况下通过“随机抽样统计聚合”把充电负荷的分布特征推出来。配合Matlab实现整个流程从建模到出图一天之内就能跑通。这篇文章就围绕“电动汽车充电负荷的蒙特卡洛预测方法”展开讲清楚建模思路、概率模型怎么定、代码怎么写以及我实际调试中踩过的坑。这个方法特别适合以下几类人参考做配电网规划或充电桩布点的工程师研究有序充电、需求响应的研究生以及刚接触蒙特卡洛模拟、想尽快上手Matlab仿真的学生。文章里的代码结构不复杂核心就是“生成大量随机车主行为逐台累加得到总负荷曲线”但细节处理不好结果很容易失真所以我会重点讲参数取值的来龙去脉和几个容易被忽略的坑。1. 为什么充电负荷预测要用蒙特卡洛1.1 充电负荷的核心难点不确定性单台电动汽车的充电行为在电网调度员眼里几乎是完全随机的车主几点插枪、出门前跑了多少公里、用快充还是慢充、电池剩多少电这些信息在事前都没法精确知道。但电网规划需要的不是某一台车的行为而是几百上千台车聚合后的总负荷曲线以及峰值出现在几点、峰值功率大概是多少范围。这个场景天然适合用概率方法处理。蒙特卡洛的核心思想很朴素既然单车的随机规律可以用概率分布描述那就按分布大量生成虚拟车主把每台车的充电曲线叠加起来模拟足够多次之后统计结果就会趋于稳定。背后是大数定律在起作用——样本量越大模拟得到的均值越接近真实期望。1.2 蒙特卡洛模拟用随机对抗随机实际做的时候整个流程可以拆成四步第一步建立每台车充电行为的概率模型比如日行驶里程服从什么分布、开始充电时间服从什么分布第二步用随机数生成器按这些分布抽取一辆车的参数计算它的充电功率曲线第三步把大量单车的曲线按时间轴叠加得到总负荷第四步重复多次模拟统计出负荷的均值、峰值、置信区间。这四步看着简单但每一步都有讲究。比如随机数生成器的质量直接影响结果稳定性抽样次数太少则曲线毛刺多、峰值不稳定。我在项目里一般单次模拟10000辆车整个蒙特卡洛循环重复200到500次既能保证收敛运行时间也控制在可接受范围内。1.3 与其他主流预测方法的对比做充电负荷预测还有其他路子。时间序列方法ARIMA、指数平滑需要较长时间的历史负荷数据适合已建成的充电站做短期预测但不适合新规划区域。神经网络类方法能拟合复杂非线性关系但需要大量高质量训练样本而且模型可解释性弱参数调整像黑盒。基于出行链的仿真方法更精细需要交通数据支撑建模复杂度高。相比之下蒙特卡洛方法的优势在于数据需求低只要有统计意义上的分布参数就能开工逻辑透明每个随机变量都能追溯到物理含义还能方便地输出概率区间而不只是一个单点预测值。缺点也很明显它对分布参数的敏感性高参数给得不合理结果会系统性偏移。所以后面我会花不少篇幅讲参数该怎么定。2. 建立随机变量的概率模型2.1 日行驶里程对数正态分布电动汽车每天的行驶里程是决定充电能量的第一关键变量。家用车的日均行驶里程统计下来通常呈正偏态分布大部分车每天跑几十公里少数车跑两百公里以上这种形态用对数正态分布拟合效果很好。Matlab里生成对数正态分布随机数的函数是lognrnd% 日行驶里程km对数正态分布 mu_m 3.20; % 均值参数 sigma_m 0.88; % 标准差参数 daily_mileage lognrnd(mu_m, sigma_m, N, 1);这里的mu_m和sigma_m不是最终里程的均值和标准差而是取对数后的均值和标准差。参考美国NHTS和国内一些出行调查数据私家车日行驶里程均值大约在30到50公里之间对应mu_m取3.0到3.3sigma_m取0.8到0.9比较常见。国内一线城市通勤距离偏长可以适当把mu_m调大一些。2.2 开始充电时刻晚高峰正态分布大部分私家车主回家后随手插枪充电所以开始充电时间呈现明显的晚高峰特征峰值集中在18点到20点。工程上常用正态分布描述均值取17.5到18.5标准差取3到4这样大部分开始充电时间落在16点到23点之间。% 开始充电时刻小时正态分布 mu_t 18.0; sigma_t 3.5; start_time normrnd(mu_t, sigma_t, N, 1);有一点要特别提醒normrnd生成的数值有可能小于0或大于24充电开始时间在物理上不可能出现这种值。处理办法是截断小于5点的按5点处理大于23点的按23点处理。截断会稍微改变分布形状但对总负荷曲线的影响不大。2.3 充电功率与电池容量参数充电功率取决于充电设施类型。居民小区慢充桩常见3.5kW和7kW两档公共快充桩功率在30kW到120kW不等。做小区配电网规划时通常按慢充为主考虑单台车功率取3.5kW或7kW。如果研究区域包含公共快充站就要按快充比例做混合建模。电池容量直接决定同样里程消耗多少电量。目前主流纯电动乘用车电池容量在40kWh到80kWh之间百公里电耗在12到18kWh。建模时可以取一个代表值比如电池容量60kWh百公里电耗15kWh也可以让电池容量在合理范围内随机分布模拟不同车型混跑的场景。2.4 充电时长的推导过程充电时长不是随便抽一个分布而是由行驶里程、能耗、电池容量和充电功率推导出来的% 充电能量kWh energy daily_mileage / 100 * e_consumption; % 充电时长小时 charge_hours energy / charge_power; % 充电结束时刻 end_time start_time charge_hours;这里默认假设车主每次回家都把电池充满实际中确实大部分慢充用户是回家就插枪直到充满。如果要做更精细的模型可以加一个起始SOC的分布然后充电时长按(1 - SOC_start) * battery_capacity / charge_power计算效果会更贴近真实场景。需要注意一个问题如果充电结束时间超过24点负荷曲线就需要折返到第二天凌晨。比如23点开始充3小时实际负荷分布在23点到次日2点。处理时可以用取模mod(end_time, 24)的方式也可以直接把时间轴扩展到48小时再截取当天部分。3. Matlab代码实现全过程3.1 初始化与环境准备代码结构我建议分成三个文件参数配置文件、核心模拟函数、主脚本。这样调参和复用都方便。参数配置把所有分布参数集中放改起来一目了然。% 参数配置 N 10000; % 模拟车辆数 num_sim 500; % 蒙特卡洛重复次数 charge_power 7; % 充电功率kW battery_capacity 60; % 电池容量kWh e_consumption 15; % 百公里电耗kWh/100km mu_m 3.20; sigma_m 0.88; % 日行驶里程对数正态参数 mu_t 18.0; sigma_t 3.5; % 开始充电时间正态参数 load_profile zeros(num_sim, 96); % 每条曲线96个点15分钟粒度时间粒度选15分钟还是1小时取决于规划精度需求。1小时粒度曲线平滑但峰值会被低估15分钟粒度能更好捕捉负荷尖峰。代码示例里我按15分钟粒度设置一天96个点。3.2 主循环抽样、聚合、累加核心模拟代码如下for s 1:num_sim daily_load zeros(1, 96); for i 1:N % 抽样 mileage lognrnd(mu_m, sigma_m); s_time normrnd(mu_t, sigma_t); s_time max(5, min(23, s_time)); % 截断 % 充电时长 energy mileage / 100 * e_consumption; c_hours energy / charge_power; e_time s_time c_hours; % 累加到15分钟粒度负荷曲线 idx_start floor(s_time * 4) 1; idx_end floor(e_time * 4) 1; idx_end min(idx_end, 96); daily_load(idx_start:idx_end) daily_load(idx_start:idx_end) charge_power; end load_profile(s, :) daily_load; end这段代码的逻辑是外层循环控制蒙特卡洛重复次数内层循环逐个生成车辆行为并累加。重点是索引换算——15分钟对应0.25小时所以用时刻乘以4取整就能定位到对应的时间段。idx_end做了上限保护防止超出数组边界。第一次跑的时候建议先设num_sim 5、N 2000快速验证流程确认没报错再加大规模。我见过不少新手一上来就设N 100000结果跑了几分钟还没结束还以为是死循环。3.3 结果可视化与置信区间绘制模拟完成后统计每条曲线的均值和分位数mean_load mean(load_profile, 1); p5 prctile(load_profile, 5, 1); p95 prctile(load_profile, 95, 1); t (0:95) * 0.25; figure; hold on; fill([t, fliplr(t)], [p5, fliplr(p95)], ... [0.85, 0.85, 0.95], EdgeColor, none); plot(t, mean_load, b-, LineWidth, 2); xlabel(时刻小时); ylabel(充电负荷kW); legend(5%-95%区间, 平均负荷);置信区间用fill填充半透明色带比单纯画多条模拟曲线清晰得多。实际运行后通常会看到双峰结构晚高峰充电开始后负荷爬升深夜很多车仍在充但新接入的车减少凌晨三四点达到峰值的概率最高早上陆续充满后负荷下降。如果做的是快充场景白天的充电峰会明显得多。3.4 向量化优化与运行效率上面双层循环有个问题模拟500次、每次10000辆车意味着要循环500万次Matlab跑起来会比较慢。我实际测试过大概要两三分钟。如果想提速可以把内层循环向量化。for s 1:num_sim mileage lognrnd(mu_m, sigma_m, N, 1); s_time normrnd(mu_t, sigma_t, N, 1); s_time max(5, min(23, s_time)); energy mileage / 100 * e_consumption; c_hours energy / charge_power; e_time s_time c_hours; for i 1:N idx_start floor(s_time(i) * 4) 1; idx_end min(floor(e_time(i) * 4) 1, 96); daily_load(idx_start:idx_end) daily_load(idx_start:idx_end) charge_power; end load_profile(s, :) daily_load; end这样改动只消除了随机抽样的循环开销累加部分的循环还在。要进一步优化可以用accumarray或稀疏矩阵做向量化累加但代码复杂度会上升。我的经验是对大多数研究和规划场景向量化到这一步已经够用没必要为缩短几十秒牺牲代码可读性。2016a之后的Matlab版本支持parfor并行循环把外层for换成parfor配合Parallel Computing Toolbox多核处理器能明显提速。要注意的是parfor里的随机数生成需要显式设置随机流否则每个工作进程可能生成相同序列parfor s 1:num_sim rng(s, twister); % 其余代码不变 end4. 结果分析与蒙特卡洛收敛性判断4.1 典型结果解读仿真结束后先看平均负荷曲线的整体形态。以10000辆私家车、7kW慢充为例晚高峰时段峰值负荷通常在3000到5000kW范围峰值出现时间大约在19点到21点之间。这里要注意“峰值出现时间”和“开始充电高峰时间”不一定重合——虽然大多数车主18点开始充电但充电功率恒定时负荷峰值取决于同时充电的车辆数而同时充电车辆数会在开始充电高峰后延后一段时间达到最大。如果发现曲线出现多个不合理的尖峰或者均值曲线明显抖动优先怀疑两个原因模拟次数太少或者抽样分布参数不合理。前者增加num_sim后者需要重新翻看分布假设。4.2 模拟次数怎么选收敛性验证方法怎么判断模拟次数够不够一个实用的做法是画出峰值负荷均值随模拟次数增加的变化曲线看它是否趋于平稳。peak_mean zeros(1, num_sim); for s 1:num_sim peak_mean(s) mean(max(load_profile(1:s, :), [], 2)); end plot(1:num_sim, peak_mean);曲线在前几十次可能波动明显随着次数增加逐渐收敛到某个稳定值。如果到500次还在持续漂移说明单次模拟的方差太大这时候增加模拟次数效果有限更有效的做法是增加单次模拟车辆数N。我在项目中常用的组合是N 10000、num_sim 300峰值均值的变化在1%以内满足工程精度。如果只是课程作业或方案初步估算N 5000、num_sim 100也能得到可用的结果。4.3 参数敏感性与多场景比较蒙特卡洛模型的输出对输入分布参数非常敏感。我做过的敏感性测试显示开始充电时间的均值每提前1小时晚高峰负荷峰值会显著降低因为更多车辆在电网低谷期充电日行驶里程的标准差增大会导致充电需求更分散负荷峰值下降但持续时间变长。因此实际应用时建议不只跑一组参数而是做多场景对比。比如基准场景正常通勤日私家车为主保守场景所有车主下班立即充电来的时间更集中乐观场景部分车主错峰充电开始充电时间更分散极端场景叠加高温天气空调使用导致百公里电耗上升每个场景跑完得到一组峰值的均值、标准差和95%置信区间规划时取保守场景的上限作为容量配置依据比单一预测值可靠得多。5. 常见问题与排查技巧实录5.1 充电负荷曲线出现“毛刺”或异常尖峰这是最常碰到的问题。毛刺通常有两个来源单次模拟车辆数N太小随机性没有充分平均掉加大N即可另一个来源是随机数发生器没有重置多次模拟结果几乎相同导致分位数区间很窄且不真实。另外检查时间索引换算是否出错。我调试时打印过前几辆车的idx_start和idx_end发现如果floor用成round同一个时刻会被分配到不同的时间槽造成阶梯状毛刺。始终用floor向下取整保证时间连续性。5.2 负荷曲线形状完全不符合预期比如晚高峰不明显或者全天负荷过于平坦。这种情况先检查开始充电时间的分布参数。如果sigma_t设得过大充电时间在一天内均匀铺开自然看不到晚高峰。还有一种可能是mu_t设得太早比如设成12点那么很多人中午就开始充电负荷曲线会在下午出现宽峰。建议先画出抽样分布的直方图确认一下histogram(lognrnd(mu_m, sigma_m, 100000, 1), 100);直方图能直观暴露参数设置错误。分布形态不对后面的模拟结果必然不对。5.3 程序运行时间过长如果模拟规模太大导致运行时间难以接受优先考虑这几招用向量化取代部分循环用parfor降低num_sim同时适当增加N因为N对结果稳定性的贡献通常比num_sim更直接最后实在不行再考虑把时间粒度从15分钟放宽到30分钟时间点从96个降到48个计算量直接减半。5.4 置信区间过宽没有参考价值95%置信区间太宽说明单次模拟的随机波动太大。原因是每辆车充电功率相同所有随机性都来自充电时间和时长极端情况下可能出现几千辆车同时充电的大尖峰。处理思路是增加充电功率的多样性比如一部分车3.5kW、一部分车7kW、少量车60kW快充让聚合曲线更平滑。6. 后续可以怎么扩展6.1 从无序充电到有序充电基础模型假设车主随心所欲地充电这个结果是“最坏情况”。如果想研究有序充电策略比如分时电价引导、智能调度控制充电功率可以在累加负荷的步骤加入控制逻辑。比如判断当前总负荷是否超过变压器容量上限超过则延迟部分车辆的充电开始时间或降低充电功率。6.2 与配电网潮流计算结合充电负荷预测的最终目的是评估配电网承载力。可以将蒙特卡洛输出的负荷曲线接入潮流计算程序比如Matpower逐时段计算节点电压、线路载流量找出重载或电压越限的时段和位置。这个扩展在工程上非常实用规划报告里“充电桩接入对配电网的影响分析”章节基本就是这个套路。6.3 多车型、多天气、多季节建模不同车型的电池容量差异很大冬季低温会让实际可用容量下降、百公里电耗上升节假日出行里程分布和通勤日完全不同。把这些因素分层建模可以让预测结果分季节、分场景输出而不是只给一条全年平均曲线。我后来做的版本就加入了温度修正系数冬季负荷峰值比夏季高10%到15%这个差异对配变容量选择是有决定性影响的。最后分享一个实际操作中的体会蒙特卡洛充电负荷预测这类模型真正决定结果质量的不是代码写得多漂亮而是对物理场景的理解深度。分布参数不是随便拍脑袋定的每一个都要能说清楚来源和适用条件。拿到一个预测结果先别急着往里套规划方案多问几句这个峰值在什么前提条件下成立参数变一变成分会怎么变想清楚这些问题模型才真正从“作业代码”变成了“决策工具”。