Loading [MathJax]/jax/output/SVG/jax.js
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

基于均值-方差理论的多阶段公交走廊设计

夏亮 江欣国 范英飞

陈维荣, 冉韵早, 韩莹, 李奇. 考虑两阶段电转气的区域综合能源系统优化调度[J]. 西南交通大学学报, 2023, 58(6): 1221-1230. doi: 10.3969/j.issn.0258-2724.20210605
引用本文: 夏亮, 江欣国, 范英飞. 基于均值-方差理论的多阶段公交走廊设计[J]. 西南交通大学学报, 2023, 58(6): 1294-1302. doi: 10.3969/j.issn.0258-2724.20210874
CHEN Weirong, RAN Yunzao, HAN Ying, LI Qi. Optimal Scheduling of Regional Integrated Energy Systems Under Two-Stage Power to Gas[J]. Journal of Southwest Jiaotong University, 2023, 58(6): 1221-1230. doi: 10.3969/j.issn.0258-2724.20210605
Citation: XIA Liang, JIANG Xinguo, FAN Yingfei. Phased Transit Service Design Based on Mean-Variance Theory[J]. Journal of Southwest Jiaotong University, 2023, 58(6): 1294-1302. doi: 10.3969/j.issn.0258-2724.20210874

基于均值-方差理论的多阶段公交走廊设计

doi: 10.3969/j.issn.0258-2724.20210874
基金项目: 国家自然科学基金(71771191)
详细信息
    作者简介:

    夏亮(1992—),男,博士研究生,研究方向为公共交通规划设计,E-mail:lxia@my.swjtu.edu.cn

    通讯作者:

    江欣国(1975—),男,教授,博士,研究方向为道路交通安全,E-mail:xjiang@swjtu.edu.cn

  • 中图分类号: U491

Phased Transit Service Design Based on Mean-Variance Theory

  • 摘要:

    随机性是公交出行需求的重要特性. 随机需求下的公交系统中,不同风险态度(中立风险、规避风险和寻求风险)的决策者会选择不同的设计方案,这在多阶段公交系统走廊设计中更为明显. 为研究不同决策风险态度下随机特性对公交服务设计的影响,首先,基于均值-方差理论,提出不同风险态度下的决策方法,并构建对应的多阶段公交走廊设计模型,同时优化不同阶段的公交线路长度、站点选址和发车间隔;之后,利用局部分解法和最优化理论,得出模型的解析解;最后,对不同设计策略(全覆盖设计、单阶段和多阶段的公交线路非全覆盖设计)和风险决策态度下的公交走廊设计进行对比分析. 结果表明:1) 公交线路多阶段设计具有比全覆盖设计和单阶段设计更低的期望系统成本,以及比单阶段设计更稳定的方案;2) 在不同设计策略中,决策风险态度对公交系统的参数配置(即公交线路长度、站点分布和发车间隔)有不同的影响,例如规避风险态度决策者在全覆盖设计中追求更密集的站点分布,而在单阶段设计中则会追求相对更低的站点密度,在多阶段设计中这就更为复杂. 该研究从政府管理者和运营商的角度,提供基于需求随机增长下公交系统多阶段设计的风险投资决策方法.

     

  • 能源是人类社会赖以生存和发展的物质基础. 针对能源型经济社会建设的目标,提出了开发新能源、减少环境污染、提高能源利用率的理念,区域综合能源系统已成为了近年研究热点[1-3]. 相较于传统单一能源系统,综合能源系统含有多种能源类型,不同能源互补耦合,基于能量梯级利用原则进行协调优化[4-5]. 目前,针对含有电转气(P2G)环节的综合能源系统已有较多研究:文献[6]建立了含P2G的电-热系统风电消纳优化运行模型,对比采用P2G设备前后和不同风电占比下风电利用率及效益;文献[7]分析了电、气、热需求响应及P2G、储能设备接入对区域综合能源系统运行优化的影响;文献[8]建立了电-气-热-冷联供的微型能源网模型,并在考虑P2G收益的基础上进行日前经济调度;文献[9]提出了一种考虑需求响应的含P2G 电-气综合能源系统优化调度模型. 上述研究中,P2G环节均采用电转天然气模型,利用P2G产生的天然气供给系统中气负荷使用. 实际上P2G环节从原理上分为2个阶段:一为电解水制氢气,二为氢气进一步与碳基原料反应合成甲烷[10]. 而电转氢气这一中间过程具有更高的转换效率[11]. 且随着氢燃料电池汽车等的推广,含有氢负荷的区域综合能源系统(RIES)日益受到重视[12]. 因此,本文建立“电-气-热-储-氢”区域综合能源系统,并综合考虑P2G的2阶段模型,将一阶段制得的氢气供给区内的氢负荷使用.

    随着新能源的日益发展,风电光伏发电场的建设装机容量不断增加,但由于其与电网建设不同步、调峰能力不匹配,导致弃风弃光现象较为明显[13-14]. 同时,区域综合能源系统耦合情况复杂,因此优化调度问题成为核心问题之一. 针对不同系统结构的优化调度问题,文献[15]对“以热定电”导致的弃风现象及运行成本较高问题,提出基于热电联供(CHP)灵活热电比的多目标优化调度方法;文献[16]提出了一种多目标电热综合能源系统低碳经济调度模型,采用多目标细菌群体驱药性算法进行求解,对系统风电消纳的效果进行了分析;文献[17]针对风电消纳问题进行研究,考虑电锅炉、电转气技术、储能设备,利用Mosek求解器求解. 文献[18-21]介绍了可转移负荷以及需求侧响应在系统能量管理方面的相关研究. 由于甲烷化反应器、电解槽、储氢罐和储气罐等装置的引入,本文“电-气-热-储-氢”RIES的拓扑结构与能量耦合情况将更为复杂,结合系统结构与运行模式的特殊性,需求侧响应影响以及建立优化调度模型问题有待深入研究.

    因此,本文针对电-气-热-储-氢耦合的区域综合能源系统,提出考虑两阶段P2G模型的RIES优化调度方法,以最小日运行成本为优化目标,充分考虑各方面约束条件,利用混合整数线性规划方法求解优化调度结果,并开展不同季节、是否含P2G环节、是否考虑需求侧响应影响因素下的技术经济分析,为RIES合理运行提供科学可靠的依据.

    考虑含P2G两阶段模型的“电-气-热-储-氢”区域综合能源系统可用图1结构描述,主要包括可再生能源发电系统、燃气轮机、燃气锅炉、电锅炉、电(气、热、氢)能源类型的负荷和储能设备、P2G环节.

    图  1  “电-气-热-储-氢”区域综合能源系统结构
    Figure  1.  Integrated energy system structure of electricity-gas-heat-storage-hydrogen coupling area
    1.1.1   传统设备模型

    1) 燃气锅炉

    燃气锅炉是以天然气为燃料,并将其转化为热能的设备. 其数学模型如式(1)所示.

    Ph,GB(t)=Pg,GB(t)ηGB, (1)

    式中:Ph,GB(t)、Pg,GB(t)分别为燃气锅炉的产热功率、天然气消耗功率,ηGB为燃气轮机的制热效率,t为时刻.

    2) 电锅炉

    电能经电锅炉转换成热能,其制热功率为

    Ph,EB(t)=Pe,EB(t)ηEB, (2)

    式中:Ph,EB(t)、Pe,EB(t)分别为电锅炉的产热功率、耗电功率,ηEB为电锅炉的制热效率.

    3) 燃气轮机

    燃气轮机能够将天然气的能量转化为电能及热能,其电、热输出功率与天然气输入功率间均存在近似线性的关系,数学模型如式(3)所示.

    {Pe,GT(t)=Pg,GT(t)ηe,GT,Ph,GT(t)=Pg,GT(t)ηh,GT,ηh,GT=1ηe,GTηloss,GT, (3)

    式中:Pg,GT(t)为燃气轮机消耗的天然气功率,Pe,GT(t)、Ph,GT(t)分别为时刻t燃气轮机的发电功率和产热功率,ηe,GTηh,GT分别为燃气轮机输出电能、热能的效率,ηloss,GT为热能自耗散率.

    1.1.2   电、热储能设备模型

    储能设备在系统运行过程中对多余能量进行储存,并在功率缺额时及时进行补偿,具有能量缓冲、平抑功率波动、保障系统稳定运行的特点,有助于实现系统的经济调度.

    电、热储能设备的充放能量为

    Ek(t)=Ek(t1)(1ηloss,k)+Δt(ηch,kPch,k(t)Pdis,k(t)/ηdis,k), (4)

    式中:k∈{e,h},e、h分别代表电、热储能设备;ηloss,k为自耗能率;Pch, k(t)、Pdis, k(t)分别为时刻t的设备充、放电功率;ηch,kηdis,k分别为储能的充、放能效率;∆t为步长,取1 h.

    1.1.3   P2G 2阶段模型

    本文构建了一种联合制氢与制甲烷的2阶段P2G模型,其流程如图2所示. 第1阶段为电转氢阶段,进行电解水反应,使用的电能主要是风力发电高峰期而电负荷相对较低时富余的电能,生成的H2中一部分注入储氢罐并供给氢燃料电池车负荷;另一部分富余H2进入第2阶段,即制甲烷阶段,进行甲烷化反应进一步制得CH4,注入储气罐并供应系统内的天然气负荷.

    图  2  P2G 2阶段模型
    Figure  2.  Two-stage P2G model

    1) 制氢阶段

    电解槽制氢量的数学模型如式(5)所示.

    VH2,EL(t)=0.0224PEL(t)ρH24νELMH2, (5)

    式中:$V_{{{\rm{H}}_2},{\rm{EL}}}^{}\left( t \right) $为时刻t电解槽制氢量(m3),PEL(t)为时刻t电解槽耗电功率,$\rho_{{\rm{H}}_2} $为H2密度(g/m3),νEL为电解氢速率,$M_{{\rm{H}}_2} $为H2摩尔质量(g/mol).

    由于储氢罐前、后压缩机的精细模型对系统优化运行的影响较小,因此,主要关注氢气的存储过程. 储氢设备数学模型表示为

    EHT(t)=EHT(t1)+(Ein,HT(t)Eout,HT(t)), (6)

    式中:EHT(t)为时刻t氢储存容量(kmol);Ein,HT(t)、Eout,HT(t)分别为时刻t存入、输出氢气的量(kmol).

    2) 制甲烷阶段

    氢气通过甲烷化反应器转化为天然气,储气罐模型同式(4),同时采用式(7)的固定效率简化模型描述甲烷化反应.

    Pg,M(t)=4ηMEH2,M(t)HLκ, (7)

    式中:Pg,M(t)为时刻t甲烷化反应的天然气功率(kW),$E_{{\rm{H}}_{2},{\rm{M}}} (t)$为时刻t输入甲烷化反应器的氢量(kmol),HL为天然气低热值(kW/m3),$ \kappa $为气体密度(kg/m3),ηM为运行效率.

    1.2.1   需求侧响应模型的目标函数

    采用激励型需求侧响应,通过调整可转移负荷的用电时间,使其与可再生能源发电在时序上更加贴合,进一步减少能源购买成本,提高系统经济性. 其目标函数为

    {minTt=1|PDG(t)Ldr(t)|,Ldr(t)=Lload(t)+Ldr_in(t)Ldr_out(t), (8)

    式中:T为调度周期,取T = 24 h,PDG(t)为时刻t可再生能源的输出功率,Lload(t)、Ldr(t)分别为时刻t响应前、后负荷大小,Ldr_in(t)、Ldr_out(t)分别为时刻t负荷转入量、转出量.

    1.2.2   需求侧响应模型的约束条件

    1) 负荷转移量约束

    在系统实际运行中,任一时刻下的实际负荷转移量均应不大于最大负荷转移量,且在一个调度周期内转出与转入的负荷总量应相等,以维持总负荷需求不变,故约束条件可表示为

    {mdr(t)Mdr(t),Ldr_in(T)=Ldr_out(T), (9)

    式中:mdr(t)、Mdr(t)分别为时刻t的负荷转移容量和最大负荷可转移容量.

    2) 负荷转移时段约束

    负荷转移只能在一个调度周期T内进行,即

    {tinT,toutT, (10)

    式中:tintout分别为负荷转入、转出时刻.

    针对电-气-热-氢耦合的区域综合能源系统的日前优化调度问题,从系统运营商角度出发,以日运行经济成本最优为目标. 目标函数F除包括购电成本Cbuye(t)、购气成本Cbuyg(t)外,还考虑弃风惩罚成本Ccurw(t)、弃光惩罚成本Ccurp(t),以体现弃风弃光现象对系统的影响,以及针对需求侧响应的补贴成本Cdr(t),如式(11)所示.

    F=minTt=1[Cbuye(t)+Cbuyg(t)+Ccurw(t)+Ccurp(t)+Cdr(t)]. (11)

    1) 系统功率平衡约束

    系统运行时需满足电、气、热、氢4种能源的能量平衡约束,即

    {Pe,net(t)+Pwt(t)+Ppv(t)+Pe,GT(t)+Pdis,e(t)=PEL(t)+Pe,EB(t)+Pch,e(t)+Le(t),Pg,net(t)+Pg,M(t)+Pg,dis(t)=Pg,GT(t)+Pg,GB(t)+Pg,ch(t)+Lg(t),Ph,GB(t)+Ph,GT(t)+Ph,EB(t)+Ph,dis(t)=Ph,ch(t)+Lh(t),Ein,HT(t)+EH2,M(t)=EH2,EL(t),2Eout,HT(t)=LH2(t), (12)

    式中:Le(t)、Lg(t)、Lh(t)、LH2(t)分别为电负荷、天然气负荷、热负荷、氢负荷功率;Pe,net(t)、Pwt(t)、Ppv(t)、Pe,GT(t)、Pe,EB(t)、PEL(t)、Pch,e(t)、Pdis,e(t)分别为购电功率、风力发电机功率、光伏系统功率、燃气轮机发电功率、电锅炉耗电功率、电解槽耗电功率、蓄电池充电功率和放电功率;Pg,net(t)、Pg,GT(t) (Pg,GB(t)/Pg,M(t))、Pg,ch(t)、Pg,dis(t)分别为购气功率、燃气轮机(燃气锅炉/甲烷化反应)耗气功率、储气设备充气功率和放气功率;Ph,GB(t)(Ph,GT(t)/Ph,EB(t))、Ph,ch(t)、Ph,dis(t)分别为燃气锅炉(燃气轮机/电锅炉)产热功率、储热设备充热功率和放热功率;$E_{{{\rm{H}}_2},{\rm{EL}}} $为甲烷化反应电解槽产氢量.

    2) 可再生能源出力约束

    光伏、风机在实际运行过程中的输出功率Pwt(t)、Ppv(t)均不能超过其预测输出功率值,故约束条件为

    {0Pwt(t)Pf,wt(t),0Ppv(t)Pf,pv(t), (13)

    式中:Pf,wt(t)、Pf,pv(t)分别为系统光伏、风机的预测最大输出功率.

    3) 电/燃气锅炉约束

    电锅炉、燃气锅炉在实际运行过程中,需满足额定功率约束及爬坡约束,即

    {0Pe,EB(t)Pe,EB,rated(t),Rg,EB,downΔtPe,EB(t)Pe,EB(t1)Rg,EB,upΔt, (14)
    {0Pg,GB(t)Pg,GB,rated(t),Rg,GB,downΔtPg,GB(t)Pg,GB(t1)Rg,GB,upΔt, (15)

    式中:Pe,EB,rated(t)、Pg,GB,rated(t)分别指系统中电锅炉、燃气锅炉的额定功率,Rg,EB,upRg,EB,down分别为电锅炉的爬坡功率上、下限,Rg,GB,upRg,GB,down分别为燃气锅炉的爬坡功率上、下限.

    4) 燃气轮机约束

    燃气轮机运行过程中需满足额定功率约束及爬坡约束,即

    {0Pg,GT(t)Pg,GT,rated(t),Rg,GT,downΔtPg,GT(t)Pg,GT(t1)Rg,GT,upΔt, (16)

    式中:Pg,GT,rated(t)为系统中燃气轮机的额定功率,Rg,GT,upRg,GT,down分别指燃气轮机的爬坡上、下限.

    5) 电解槽约束

    电解槽运行过程中需满足其上、下限功率约束,即

    uEL(t)PEL,minPEL(t)uEL(t)PEL,max, (17)

    式中:PEL,minPEL,max分别为电解槽运行功率的下、上限;uEL(t)为0-1变量,表示电解槽时刻t的启停状态,uEL(t) = 1即表示电解槽此时处于运行状态,uEL(t) = 0即表示电解槽此时处于停机状态.

    6) 甲烷化反应器约束

    uM(t)EH2,M,minEH2,M(t)uM(t)EH2,M,max, (18)

    式中:EH2,M,minEH2,M,max分别指甲烷化反应器耗氢量的下、上限;uM(t)为0-1变量,表示甲烷化反应器t时刻的启停状态,uM(t) = 1表示启动,反之表示停机.

    7) 储电/热/气设备

    电、热、气储能设备需满足上、下限约束,避免过度充放能影响设备寿命. 同时,为了给下一调度周期预留一定的调节裕量,使得在下一个调度周期开始时储能设备能够满足系统充放能要求,所以将运行一个调度周期T后的储能容量恢复到初始状态. 此外,考虑到同一时刻储能装置不能同时充放能,所以任一时刻充放能功率都必有一个为0. 故相关约束条件为

    {0Pch,κ(t)uκ(t)Pch,κ,max,0Pdis,κ(t)vκ(t)Pdis,κ,max,0uκ(t)+vκ(t)1, (19)
    Eκ,minEκ(t)Eκ,max, (20)
    Eκ(0)=Eκ(T), (21)

    式中:κ∈{e,h,g},g代表气设备;Pch,κ(t)、Pdis,κ(t)分别为储能设备κ充、放能功率;uκ(t)、vκ(t)均为0-1变量,表示时刻t储能设备κ的充、放能状态,当uκ(t)=0时,储能设备不充能,反之则充能;vκ(t)=0时,储能设备不放能,反之则放能;同一时刻储能设备不能同时充、放能;Eκ,minEκ,max分别为κ储能设备允许容量下、上限;Pch, κ,maxPdis, κ,max分别为储能设备κ充能和放能功率的最大值.

    8) 储氢罐

    储氢罐也需满足上、下限约束及调度周期T内充放氢总量相等,即

    {0Eout,HT(t)Eout,HT,max,0Ein,HT(t)Ein,HT,max, (22)
    EHT,minEHT(t)EHT,max, (23)
    EHT(0)=EHT(T), (24)

    式中:Eout,HT,maxEin,HT,max分别为各时刻储氢罐注入、输出的氢气最大值,EHT,minEHT,max分别为储氢罐储存容量上、下限,EHT(0)和EHT(T)分别为当前优化调度周期的初时刻储氢量与末时刻储氢量.

    9) 外部网络交互功率约束

    在本研究中,仅考虑从外部电网、天然气网购入能量,不考虑向外部电网售出能量. 系统与外部网络间联络线上的传输功率有如下关系:

    {0Pg,net(t)Pg,net,max(t),0Pe,net(t)Pe,net,max(t), (25)

    式中:Pg,net(t)、Pe,net(t)分别为购气功率和购电功率;Pg,net,max(t)、Pe,net,max(t)分别为购气和购电功率的最大值.

    电-气-热-储-氢耦合的区域综合能源系统中耦合了多种类型的设备,其优化变量众多,本质为混合整数规划(MILP)问题. 因此,搭建系统优化调度模型后,利用数学优化工具包YALMIP建模求解,实现系统内多种设备的优化调度.

    1) 输入数据:根据冬季、夏季典型日的风光资源条件,输入其风力、光伏发电预测值及其相应的各类型负荷数据;输入各设备模型的相关数据信息.

    2) 定义优化调度决策变量:定义一个决策变量矩阵P,第1~24列分别代表1 d的24个时段.

    3) 设定优化调度目标函数:日运行成本最优.

    4) 设定优化调度约束条件:结合各设备模型及系统功率平衡原则分别建立相应的约束条件,并输入微电网优化调度模型.

    5) 求解参数设置:设置求解算法,对决策变量矩阵P进行求解.

    6) 输出结果:通过优化算法求解输出满足系统所有约束条件下的最低日运行成本以及相应的系统最佳优化调度策略.

    本文采用的冬季典型日电、气、热、氢负荷数据如图3,夏、冬季典型日热负荷值比例取为1∶4. 分时电价与天然气价如图4所示. 风光出力预测数据如图5所示.

    图  3  冬季典型日负荷预测数据
    Figure  3.  Predicted data of normal daily load in winter
    图  4  电价与天然气价
    Figure  4.  Electricity and gas prices
    图  5  冬季典型日风光出力预测数据
    Figure  5.  Predicted normal output data of wind turbine and photovoltaic system in winter

    为研究P2G对系统的影响,本文按照是否含有P2G环节的冬、夏季节设置4个场景分析系统的优化调度结果,具体场景见表1.

    表  1  算例场景设置
    Table  1.  Scenario settings
    场景季节含电/热储能含气储能含 P2G
    1 × ×
    2
    3 × ×
    4
    下载: 导出CSV 
    | 显示表格

    加入P2G环节显著减少了系统日运行成本,提升了系统的经济性,不同场景下的系统运行成本见表2.

    表  2  不同场景下的系统运行成本
    Table  2.  Operating costs in different scenarios
    场景 购电成本 购气成本 购氢成本 弃风成本 弃光成本 总成本
    1 557.5 5863.6 1079.6 1099.4 56.6 8656.7
    2 401.0 5420.8 0 10.7 0 5832.5
    3 0 3025.6 1079.6 1426.6 207.7 5739.5
    4 0 2152.3 0 49.5 0 2201.8
    下载: 导出CSV 
    | 显示表格

    1) 购电购气成本

    分别对比场景1、2、3、4,购电购气成本均有不同程度的降低. 其中购气成本降幅更大,场景2较场景1减少442.8元,场景4较场景3减少873.3元,分别节省总成本的7.55%和28.86%. 对比可知,加入P2G环节后,系统将多余电能转化为天然气以满足气负荷,降低购气成本,尤其是在10:00—15:00 的风光发电高峰时段. 由于夏季热负荷需求降低,锅炉出力减少,更多的弃风量被P2G环节利用,所以场景4购气成本较场景2节约比例更大.

    2) 购氢成本

    场景1和场景3不含P2G环节,系统内氢负荷需求全部由加氢站购氢满足,氢负荷一天内的总量约为53.98 kg,则购氢成本固定为1079.6元,分别占场景1和场景3总日运行成本的12.47%和18.81%.

    3) 弃风弃光成本

    图35可知:夜间0:00—6:00电负荷需求较低,均小于1.0 MW,而此时风力发电出力最大可达3.5 MW,因此,夜间产生较大弃风量,系统经济性较差;日间电负荷需求于日间11:00—15:00明显降低,此时风电出力趋于稳定,且光伏出力处于高峰,故午间产生较严重的弃风弃光现象. 综上,场景1、3由于缺少P2G环节,都存在较严重的弃风现象和一定的弃光现象,而场景2、4加入P2G后,风光能源均得到良好利用,弃风弃光成本大幅降低.

    图6图7为各场景下风光利用情况. 对于场景1、3,午间时段光照较强,产生一定的弃光量,风机不出力;夜间低负荷时段,出现大量弃风,日弃风量分别为45.26%和58.73%. 场景2、4消除了弃光现象,仅在01:00—03:00和20:00—23:00时间段存在少量弃风,弃风率降低到0.44%和2.04%. 对比场景1、3,夏季热负荷需求较冬季更低,弃风弃光量更多,加入P2G环节后消纳弃风弃光效果更显著. 综上,加入P2G后电-气-氢-储-热RIES能有效消纳弃风弃光可再生能源,显著减少弃风弃光现象.

    图  6  风电出力及风机利用率
    Figure  6.  Optimal output of wind turbine and utilization
    图  7  光伏出力及光伏利用率
    Figure  7.  Optimal output of photovoltaic system and its utilization

    以系统日运行成本最小为目标求解得优化调度结果,以冬季场景为例进行分析,图89分别为场景1、2的电、气、热功率平衡图.

    图  8  场景1的功率平衡图
    Figure  8.  Balance of power in scenario 1
    图  9  场景2的功率平衡图
    Figure  9.  Balance of power in scenario 2

    图8可以看出:场景1在此时段内尽管电锅炉也满载运行,但由于缺少电解水制氢,仍无法消纳大量的弃风弃光量;只有蓄电池少量储放电对系统进行调节;由于场景1缺少制甲烷环节和储气罐,系统日购气量一直保持在较高值,以满足气负荷的需求. 观察图9可以看出:10:00—15:00时段,由于光伏出力达高峰且电负荷需求有所减少,电转气设备出力增加满足气负荷需求,购气量减少;热负荷需求降低,燃气锅炉出力减少,购电仅出现在16:00—21:00时间段内,是由于此时段内电负荷需求达到峰值,各设备无法满足其需求,故从外部购电来满足此功率缺额;场景2在01:00—16:00及22:00—24:00时段系统均未从外部购电,燃气轮机不出力,电锅炉满载运行,同时P2G环节出力制氢制甲烷. 由图5可知:风光出力之和在此时段内均保持较高水平,因此,存在富余电能供给电锅炉和P2G环节以减少弃风弃光量,提高系统经济性;燃气锅炉出力满足剩余热负荷需求,甲烷化反应补充供气. 各储能系统灵活充放能调节各能量需求,提高系统缓冲能力.

    图10为场景2、4电转气环节出力. 对比图10(a)和图5:场景2中电解槽产氢时段与场景1中的弃风弃光时段几乎重合,场景4与场景3亦然,即场景2、场景4中P2G环节均能为系统消纳多余可再生能源;由于储氢罐始末态一致,第1阶段电解水产氢输入储氢罐的氢量等于氢负荷日需求量. 由于夏季存在更多弃风弃光量,电转气制甲烷量更大,即图10(b)第2阶段甲烷化耗氢量更多,P2G环节运行时间略长.

    图  10  P2G环节出力
    Figure  10.  Optimal output of P2G

    场景5在场景2的基础上引入电力需求侧响应环节,图11为可再生能源出力曲线及需求响应前后电负荷曲线对比,表3为冬季各场景下的运行成本.

    图  11  冬季典型日电负荷需求响应前、后对比
    Figure  11.  Comparison of normal daily electric load demands in winter before and after introducing response
    表  3  冬季不同场景下的系统运行成本
    Table  3.  Operating costs in different scenarios in winter
    场景 购电成本 购气成本 购氢成本 弃风成本 弃光成本 补贴成本 总成本 节省成本/%
    1 557.5 5863.6 1079.6 1099.4 56.6 0 8656.7
    2 401.0 5420.8 0 10.7 0 0 5832.5 32.62
    5 22.6 5152.3 0 17.7 0 547.8 5740.4 33.69
    下载: 导出CSV 
    | 显示表格

    引入需求侧响应后,由于可转移负荷用电时段的改变,使得响应后负荷曲线更加平滑,且更加贴近可再生能源出力曲线,使得与外部电网的交互功率减少,大幅降低了购电成本;使得燃气轮机出力有所减少,购气成本降低;与场景2相比,弃风弃光成本仍处于较低水平. 因此,引入需求侧响应后,尽管新增了补贴成本,但总成本仍有所降低,相较场景2节省成本百分比进一步提升为33.69%.

    1) 引入P2G环节后,冬、夏2季节下的系统在购电购气成本、弃风弃光成本、购氢成本各方面都有所下降,使得系统总运行成本大幅降低,提升系统的经济性.

    2) 引入P2G环节后,配合各储能设备灵活充放能,系统能有效消纳弃风弃光可再生能源,显著减少弃风弃光现象. 夏季弃风弃光现象更严重,P2G环节运行时间更长,充分提高了系统能源利用效率.

    3) 考虑电负荷需求侧响应,通过负荷用电时间转移进行削峰填谷,使得与外部网络交互功率减少,通过减少能源购买成本,降低了系统的日运行成本.

  • 图 1  单中心公交走廊

    Figure 1.  Layout of monocentric corridor

    图 2  公交线路全覆盖设计的发车间隔

    Figure 2.  Headways of transit system with full-covered design

    图 3  公交线路单阶段设计的发车间隔

    Figure 3.  Headways with once-and-done design

    图 4  公交系统多阶段设计结果

    Figure 4.  Result of transit system under phased design

    表  1  实验中的参数值

    Table  1.   Value of parameters in experiment

    参数数值参数数值
    vw/(km·h−1 2 ζ 1
    vc/(km·h−1 60 k/(户·km−1 4600
    Td/(h·站−1 0.0125 $\bar {g} $/(户·年−1 0.03
    CV/(人·车−1 1800 σg/(户·年−1 0.01
    π1/(h·(km·d)−1 14256 + 475.2μ pg 0.7
    π2/ (h·(站·d)−1 7056 + 235.2μ κ 0.1
    π3/(h·km−1 2.2 Tw/(d·年−1 250
    π4/(h·h−1 101 + 5 μ N0/ 户 30000
    Tr/ (h·d−1 16 μ/ (h·h−1 20
    η/ (次·(户·d) −1 4
    下载: 导出CSV

    表  2  公交线路全覆盖设计结果

    Table  2.   Result of transit system with full-covered design

    风险决策
    态度
    公交线路/
    走廊长度/km
    平均站点
    密度/
    (站·km−1
    平均发
    车间隔/
    min
    期待用户
    通勤成本/
    (h·(户·d)−1
    期待运营
    商成本/
    (h·(户·d)−1
    期待系统
    成本/
    (h·(户·d)−1
    系统成本
    标准差/
    (h·(户·d)−1
    中立风险 (6.81, 14.93) 1.0937 2.1989 1.1280 0.8123 1.9403 0.0794
    规避风险 (6.81, 14.93) 1.1773 2.1721 1.1261 0.8154 1.9414 0.0761
    寻求风险 (6.81, 14.93) 0.9814 2.2252 1.152 0.8076 1.9595 0.1057
    下载: 导出CSV

    表  3  公交线路单阶段设计结果

    Table  3.   Result of transit system with once-and-done design

    风险决策
    态度
    线路
    长度/km
    走廊
    长度/km
    平均站
    点密度/
    (站·km−1
    平均发车
    间隔/min
    期待用户
    通勤成本/
    (h·(户·d)−1
    期待运营
    商成本/
    (h·(户·d)−1
    期待系统
    成本/
    (h·(户·d)−1
    系统成本
    标准差/
    (h·(户·d)−1
    中立风险 10.97 (6.81, 14.93) 1.1011 2.2979 1.208 0.8504 2.0584 0.1334
    规避风险 11.40 (6.81, 14.93) 1.0976 2.3131 1.1813 0.8827 2.064 0.1169
    寻求风险 10.48 (6.81, 14.93) 1.1099 2.2766 1.2499 0.8144 2.0644 0.1546
    下载: 导出CSV

    表  4  公交线路多阶段设计结果

    Table  4.   Result of transit system with phased design

    风险决策
    态度
    线路长度/km 走廊长
    度/km
    平均站
    点密度/
    (站·km−1
    平均发车
    间隔/min
    阶段
    数/个
    期待用户
    通勤成本/
    (h·
    (户·d)−1
    期待运营
    商成本/(h·
    (户·d)−1
    期待系统
    成本/
    (h·(户·d)−1
    系统成本
    标准差/(h·
    (户·d)−1
    初始 最终
    中立风险 6.82 12.31 (6.81, 14.93) 1.2080 2.0977 5 1.1539 0.7398 1.8936 0.1095
    规避风险 6.82 12.76 (6.81, 14.93) 1.2186 2.0930 5 1.1420 0.7580 1.9000 0.0947
    寻求风险 6.82 11.88 (6.81, 14.93) 1.2209 2.0916 5 1.1856 0.7147 1.9003 0.1352
    下载: 导出CSV
  • [1] FUSCO G, ALESSANDRINI A, COLOMBARONI C, 等. A model for transit design with choice of electric charging system[J]. Procedia - Social and Behavioral Sciences, 2013, 87: 234–249.

    FUSCO G X,ALESSANDRINI A,COLOMBARONI C,et al. A model for transit design with choice of electric charging system[J]. Procedia-Social and Behavioral Sciences,2013,87: 234-249.
    [2] FAN W B, RAN Y. Planning skip-stop services with schedule coordination[J]. Transportation Research Part E: Logistics and Transportation Review, 2021, 145: 102119.1-102119.19. doi: 10.1016/j.tre.2020.102119
    [3] 郭小乐. 高铁站接运公交线路设计与运营方案优化研究[D]. 北京: 北京交通大学, 2020.
    [4] LUO X L, FAN W B, JIANG Y S, et al. Optimal design of bus stop locations integrating continuum approximation and discrete models[J]. Journal of Advanced Transportation, 2020, 2020: 8872748.1-8872748.10.
    [5] 王殿海,熊满初,章立辉,等. 基于舒适度的公交线网优化方法[J]. 西南交通大学学报,2017,52(4): 789-795. doi: 10.3969/j.issn.0258-2724.2017.04.018

    WANG Dianhai, XIONG Manchu, ZHANG Lihui, et al. Optimization of public traffic network based on comfort improvement[J]. Journal of Southwest Jiaotong University, 2017, 52(4): 789-795. doi: 10.3969/j.issn.0258-2724.2017.04.018
    [6] 梁东,李远富,樊敏. 基于区间数距离的铁路线路方案决策模型与方法[J]. 西南交通大学学报,2019,54(4): 823-830. doi: 10.3969/j.issn.0258-2724.20180847

    LIANG Dong, LI Yuanfu, FAN Min. Decision making model and method based on distance measure between interval numbers in railway location[J]. Journal of Southwest Jiaotong University, 2019, 54(4): 823-830. doi: 10.3969/j.issn.0258-2724.20180847
    [7] MATISZIW T C, MURRAY A T, KIM C. Strategic route extension in transit networks[J]. European Journal of Operational Research, 2006, 171(2): 661-673. doi: 10.1016/j.ejor.2004.09.029
    [8] CHENG W C, SCHONFELD P. A method for optimizing the phased development of rail transit lines[J]. Urban Rail Transit, 2015, 1(4): 227-237. doi: 10.1007/s40864-015-0029-2
    [9] SUN Y S, SCHONFELD P, GUO Q W. Optimal extension of rail transit lines[J]. International Journal of Sustainable Transportation, 2018, 12(10): 753-769. doi: 10.1080/15568318.2018.1436730
    [10] XIA L, JIANG X G, FAN W B. Continuum approximation model for optimizing phased development of scalable transit corridors[C]//2019 IEEE Intelligent Transportation Systems Conference (ITSC), Auckland: IEEE, 2019: 852-857.
    [11] HASSANNAYEBI E, ZEGORDI S H, AMIN-NASERI M R, 等. Train timetabling at rapid rail transit lines: A robust multi-objective stochastic programming approach[J]. Operational Research, 2017, 17(2): 435–477.

    HASSANNAYEBI E,ZEGORDI S H,AMIN-NASERI M R,et al. Train timetabling at rapid rail transit lines: a robust multi-objective stochastic programming approach[J]. Operational Research,2017,17(2): 435-477.
    [12] HASSANNAYEBI E, ZEGORDI S H, AMIN-NASERI M R, 等. Demand-oriented timetable design for urban rail transit under stochastic demand[J]. Journal of Industrial and Systems Engineering, 2016, 9(3): 28–56.

    HASSANNAYEBI E,ZEGORDI S H,AMIN-NASERI M,et al. Demand-oriented timetable design for urban rail transit under stochastic demand[J]. Journal of Industrial and Systems Engineering,2016,9(3): 28-56.
    [13] HØYEM H, ODECK J. Optimal public transit frequency under stochastic demand and fixed vehicle size: application in the Norwegian car ferry sector[J]. Research in Transportation Economics, 2020, 82: 100878.1-100878.14. doi: 10.1016/j.retrec.2020.100878
    [14] HUANG Z F, REN G, LIU H X. Optimizing bus frequencies under uncertain demand: case study of the transit network in a developing city[J]. Mathematical Problems in Engineering, 2013, 2013: 1-10.
    [15] AN K, LO H K. Robust transit network design with stochastic demand considering development density[J]. Transportation Research Part B: Methodological, 2015, 81: 737-754. doi: 10.1016/j.trb.2015.05.019
    [16] LI M, XU J, WEI L, 等. Modeling a risk-based dynamic bus schedule problem under no-notice evacuation incorporated with dynamics of disaster, supply, and demand conditions[J]. Journal of Advanced Transportation, 2019, 2019: 1–17.

    LI M H,XU J L,WEI L Y,et al. Modeling a risk-based dynamic bus schedule problem under no-notice evacuation incorporated with dynamics of disaster, supply, and demand conditions[J]. Journal of Advanced Transportation,2019,2019: 1-17.
    [17] HADAS Y, SHNAIDERMAN M. Public-transit frequency setting using minimum-cost approach with stochastic demand and travel time[J]. Transportation Research Part B: Methodological, 2012, 46(8): 1068-1084. doi: 10.1016/j.trb.2012.02.010
    [18] AMIRIPOUR S M M, CEDER A , MOHAYMANY A S. Designing large-scale bus network with seasonal variations of demand[J]. Transportation Research Part C: Emerging Technologies, 2014, 48: 322-338. doi: 10.1016/j.trc.2014.08.017
    [19] 魏长钦,王伟智. 随机需求下定制公交站点及路径动态优化模型[J]. 福州大学学报(自然科学版),2020,48(1): 98-104.

    WEI Changqin, WANG Weizhi. Dynamic optimization model of customized bus stations and route under random demand[J]. Journal of Fuzhou University (Natural Science Edition), 2020, 48(1): 98-104.
    [20] YAN Y, MENG Q, WANG S, 等. Robust optimization model of schedule design for a fixed bus route[J]. Transportation Research Part C, 2012, 25: 113–121.

    YAN Y D,MENG Q,WANG S A,et al. Robust optimization model of schedule design for a fixed bus route[J]. Transportation Research Part C: Emerging Technologies,2012,25: 113-121.
    [21] MEDINA M, GIESEN R, MUÑOZ J C. Model for the optimal location of bus stops and its application to a public transport corridor in Santiago, Chile[J]. Transportation Research Record: Journal of the Transportation Research Board, 2013, 2352(1): 84-93. doi: 10.3141/2352-10
    [22] 李欢. 基于连续近似法的公交线路优化及实例分析[D]. 成都: 西南交通大学, 2018.
    [23] CHOI T M, LI D, YAN H M. Mean–variance analysis for the newsvendor problem[J]. IEEE Transactions on Systems, Man, and Cybernetics−Part A: Systems and Humans, 2008, 38(5): 1169-1180. doi: 10.1109/TSMCA.2008.2001057
    [24] FAN W B, MEI Y, GU W H. Optimal design of intersecting bimodal transit networks in a grid city[J]. Transportation Research Part B: Methodological, 2018, 111: 203-226. doi: 10.1016/j.trb.2018.03.007
    [25] GU W H, AMINI Z, CASSIDY M J. Exploring alternative service schemes for busy transit corridors[J]. Transportation Research Part B: Methodological, 2016, 93: 126-145. doi: 10.1016/j.trb.2016.07.010
    [26] DAGANZO C F. Structure of competitive transit networks[J]. Transportation Research Part B: Methodological, 2010, 44(4): 434-446. doi: 10.1016/j.trb.2009.11.001
    [27] LI Z C, LAM W H K, WONG S C, 等. Modeling the effects of integrated rail and property development on the design of rail line services in a linear monocentric city[J]. Transportation Research Part B: Methodological, 2012, 46(6): 710–728.

    LI Z C,LAM W H K,WONG S C,et al. Modeling the effects of integrated rail and property development on the design of rail line services in a linear monocentric city[J]. Transportation Research Part B: Methodological,2012,46(6): 710-728.
  • 加载中
图(4) / 表(4)
计量
  • 文章访问数:  314
  • HTML全文浏览量:  100
  • PDF下载量:  28
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-11-08
  • 修回日期:  2022-01-25
  • 网络出版日期:  2023-08-25
  • 刊出日期:  2022-04-01

目录

/

返回文章
返回