本发明提出一种蓄热式电锅炉参与电力系统调峰的运行优化方法，属于灵活性资源参与调峰的运行优化领域。该方法首先利用环境温度预测值和环境温度实际值的历史数据通过经验概率分布函数转换后，构建转换后的环境温度预测值和环境温度实际值的联合概率分布函数，计算得到次日实际环境温度的经验分布；然后建立由目标函数和约束条件构成的蓄热式电锅炉参与电力系统调峰的运行优化模型，对该模型的目标函数进行分布鲁棒优化后对模型求解，得到各时段蓄热式电锅参与调峰时上报功率的最优解，最终得到蓄热式电锅炉参与调峰市场的投标功率曲线，优化完毕。本发明能够降低环境温度不确定性对蓄热式电锅炉参与调峰的影响，充分发挥蓄热式电锅炉的调峰能力。