陈文溆乐， 向 月， 彭光博， 刘友波， 刘俊勇

(四川大学 电气工程学院， 成都 610065)

中国提出力争在2030年前实现碳达峰，努力争取在2060年前实现碳中和.从二氧化碳的排放源分析，与能源有关的碳排放占比近九成，其中电力的碳排放又占到八成，能源供给侧变革成为碳中和技术的主线[1-3].因此，研究双碳目标下未来电力系统中能源供给侧结构如何演化十分必要.

系统动力学(System Dynamics, SD)于1956年创立，在20世纪70年代被引进中国.通过对系统论、信息论、控制论进行有机结合，SD能够定性和定量地分析多变量、多重反馈、高阶次、高度非线性的变量关系，是解决复杂系统中多元信息反馈问题的有效方法之一，在解决电力系统相关领域的问题方面也体现出了一定的优越性[4-5].

基于中国输配电价政策，文献[6]研究了输配电价政策对电网企业现金流的影响，运用SD理论建立了企业投资优化决策模型，并以某市为例对电力市场改革后电网企业的投资决策提出建议和政策含义.文献[7]搭建了分布式光伏电能共享与配电网层面交互模型的SD模型，分析了不同分布式光伏安装率和上网电价下的共享价值变化.文献[8]通过系统动力学模型生成多种经济发展场景来模拟供给侧结构性改革的不确定性，实现中长期负荷外推预测.

系统动力学在电价分析、电力投资经济性、负荷预测方面已有较多研究，但在能源规划方面的应用较少，尤其是电力系统能源供给侧改革中少有应用.在“碳中和、碳达峰”背景下，通过对比其他国家，围绕技术链和产业链，文献[9]提出中国能源高质量发展的3个显性途径和1个隐性途径，建议将关键核心技术纳入国家科技计划.文献[10]调研国内外智能电网的背景、发展现状、关键技术，对我国建设智能电网的关键技术进行分析.文献[11]将产能投资和设备资产动态演变解耦，建立能源替代模型，揭示了天然气消费替代带来的能源供应安全冲击和能源设备资产搁浅的形成机理.基于我国30个省份的面板数据，文献[12]评估了碳交易市场规模的碳减排效应.文献[13]以发电市场碳减排政策为研究对象，构建了基于SD理论的发电行业市场节能减排政策影响模型，分析不同政策情景下，对碳减排及电力市场的影响.已有研究中对能源发展路径、智能电网主要技术、电能替代、减碳效应等进行了分析，但限于研究时未考虑电力系统中供给侧能源演化发展与碳排放之间的相互影响关系，无法定量分析供给侧能源发展路径中各因素对碳减排的贡献程度.

针对以往研究的不足，本文在长时间尺度下探讨不同情景碳排放可能发生的变化，从电力系统发电结构转型的角度出发，提出一种基于系统动力学的能源网发展演化模型，综合考虑双碳目标的政策指导、技术发展应用、电力市场潜力等影响因素.

1 系统动力学模型构建

SD首先以实际观测到的初始系统信息为基础，然后再辅以计算机仿真模拟技术帮助分析者进行路径优化从而预测系统行为变化，具体系统动力学建模步骤如图1所示.

图1 系统动力学模型建模步骤Fig.1 Modeling steps of system dynamics model

2 SD模型构建

2.1 分析提炼关键因素

基于SD建模原理，根据《能源生产和消费革命战略(2016—2030)》[14]对能源供给侧结构性改革作出的要求，非化石能源在一次能源消费中的比重、非化石能源发电量占比、电力系统碳排放总量都是重要考虑因素.人口数量、全社会用电量、国民生产总值的变化会引起电力需求的变化.电力需求决定总发电量，结合装机容量、投资力度、发电消纳水平，也会对新能源发电量起到一定的影响作用.绿证数量与新能源发电量有对应关系，其规模和交易直接造成发电结构的变化，进而影响电力碳排放.

总之，当电力碳排放呈增长趋势时，在双碳背景下，国家会对行业减排制定相应的政策措施.搭建图2所示的因果回路图，包括电力需求模块、电力技术模块、发电消纳模块、绿证市场模块以及能源消耗和碳排放模块.

图2 因果回路图Fig.2 Causal loop diagram

2.2 模型主要关系及参数说明

2.2.1政策投资模块 政策投资模块模拟政府投资行为，主要影响新能源新增装机容量，动力学演化模型如图3所示.

图3 政策投资模块Fig.3 Policy investment module

It,b=It,pkt,i

(1)

(2)

Ct,inc=kt,b-cIt,b+Et,ndkt,e-cPt0,gp/(Pt,gptt,c)

(3)

(4)

式中：It,b、kt,i为投资变化和投资变化率；It,p为投资力度；It0,p为仿真初始年限的投资力度；投资力度会促进新能源新增装机，具体由式(3)给出，其表明新能源新增装机容量Ct,inc不仅与投资变化有对应关系，还受电力需求Et,nd、当年绿证价格Pt,gp、新能源平均利用时间tt,c的影响，kt,b-c、kt,e-c为相应系数.式(4)给出新能源装机容量Ct,c计算式，由仿真初始年限新能源装机容量Ct0,c和新能源新增装机容量Ct,inc计算得到.

2.2.2电力需求模块 用电需求是推动电力碳排放的正向因素，用电需求主要受经济发展水平、人口、全国发电企业平均上网电价、产业结构等因素的影响.用电需求和经济发展水平呈现高度相关性，特别是第二产业和第三产业的用电量和相应产业GDP增加值的相关性在99%以上[15].人口数量的稳步增加也会对电力消费及碳排放产生较大影响，电力需求子模块的动力学演化模型如图4所示.

图4 电力需求模块Fig.4 Power demand module

(9)

式中：Pt为仿真年限人口数量；Pt0为仿真初始年限人口数量；pt,inc为人口增长率；Gt为国民生产总值；Gt0为仿真初始年限GDP；gt,inc为GDP增长率.式(6)、(7)分别给出仿真年限人口数量和国民生产总值计算式，Pt,G为人均收入.Et0,h、kh,inc分别为仿真初始年限的全行业用电量、行业用电增长率.全行业用电量Et,h用式(8)表示，总电力需求用式(9)表示，n1、n2、n3、n4分别为人均收入Pt,G、全国发电企业平均上网电价Pt,avg、第二产业产值Pt,s以及全行业用电量Et,h的权重系数.

2.2.3电力技术模块 电力技术子模块主要考虑发电新技术、清洁能源装机和发电占比、发电设备利用时间、线损率等因素.碳捕获、利用与封存(CCUS)技术通过降低源端发电排放物实现减排目标.多种清洁能源发电对电力碳排放起到一定的抑制作用.终端电力需求一定的条件下，提高发电设备的利用时间，发电量和碳排放能得到有效降低.系统动力学推演模型如图5所示.

图5 电力技术模块Fig.5 Power technology module

Qt,pro=p1Et,nd+p2Pt,avg+p3Et,l

(10)

qt,c=k1ct,c+k2It,m

(11)

ct,c=Ct,c/(Ct,c+Ct,f)

(12)

Qt,cpro=Qt,proqt,c

(13)

Qt,fpro=Qt,pro(1-qt,c)

(14)

Qt,coal=(Qt,fpro+Qt,x)ke-c

(15)

tt,f=Qt,fpro/Ct,f

(16)

tt,c=Qt,cpro/Ct,c

(17)

式中：Qt,pro为总发电量，受电力需求Et,nd、平均上网电价Pt,avg和线损率Et,l的影响；p1、p1、p1为相关系数；新能源发电量所占比重qt,c主要受到装机容量占比ct,c和投资力度的影响，k1、k2为相应相关系数.装机容量占比又由新能源装机容量Ct,c和火电企业装机容量Ct,f直接决定.新能源发电量Qt,cpro、火电企业发电量Qt,fpro由式(13)、(14)算出.Qt, coal为火电企业的煤耗量；Qt,x为火电企业参与辅助服务电量；Ct,f为火电企业装机容量.火电企业机组发电量和参与辅助服务的电量决定火电企业的煤耗量，ke-c为电煤转换系数，关系式如式(15)所示.式(16)、(17)分别给出火电企业设备平均利用时间tt,f、新能源平均利用时间tt,c的计算式.

2.2.4绿证市场模块 绿证市场模块的动力学演化模型如图6所示，主要考虑煤耗成本、煤价、绿证收益等因素.煤耗成本和煤价对电力碳排放起到一定的抑制作用，其动力学演化模型耦合在电力需求子模块及能源消耗和碳排放子模块中.绿证是衡量清洁能源发电量的指标，新能源发电量越多，对电力碳排放的抑制作用越强.在市场中销售绿电，新能源发电企业得到相应的收益，也会促进新能源机组的发展.

图6 绿证市场模块Fig.6 Power market module

Pt,gp=Pt0+Pt,ex

(18)

Pt,ex=(qt,pg-qt,sg)/qt,sg

(19)

Nt,ger=Nt0,ger+Gt,l-Gt,grid

(20)

Nt,grid=Nt0,grid+Gt,grid-Gt,up

溧阳南山片区,地处溧阳、宜兴、广德三地交汇处,自然条件较好,周边旅游资源丰富,是旅游发展“三山两湖一团城”体系中的重点发展景区(详见图2)。

(21)

Gt, l=Qt,cpro/kcpro-l

(22)

Gt,up=Et,ndPq

(23)

式中：Pt,gp为绿证价格变化，由绿证初始价格Pt0和绿证超额需求Pt,ex决定；qt,pg为预计购买量；qt,sg为预计销售量；Nt,ger为绿证厂商所持有的绿证数量；Nt,grid为电网所持有的绿证数量.超额需求是预计购买量和预计销售量的差值函数，如式(19)所示.式(20)、(21)给出绿证厂商和电网所持有的绿证数量，其中间数量关系由式(22)、(23)给出，Nt0,ger、Nt0,grid为仿真初始年限绿证厂商、电网持有的绿证数量；Gt,l为新能源发电量所对应的绿证数量，也就是绿证厂商获得的绿证数量；Gt,grid、Gt,up分别为电网获得的绿证和上交的绿证；kcpro-l为新能源发电量和绿证数量的转换关系；Pq为配额比例.

2.2.5发电消纳模块 发电消纳模块的动力学演化模型如图7所示.储能和火电机组灵活性改造是消纳清洁能源的重要方式.在新能源资源丰富时，储能的装机容量直接影响弃新能源的多少，火电机组减少出力也能最大程度减少弃新能源量.在新能源出力较低时，储能放电和增大火电机组出力都是缓解新能源出力不够的有效措施.

图7 发电消纳模块Fig.7 Power generation consumption module

(30)

式中：Qt,x为消纳量；Qt0,x为仿真初始年限消纳量；kt,x为消纳水平变化率.在本模型中，消纳水平变化率主要受煤电消纳因子ct,x和储能消纳因子st,x的影响.煤电消纳水平通过火电机组灵活性指数kt,fe体现，储能消纳水平通过储能装机容量Ct,s和充放电效率ηt,s体现，Qt,x、Qt,s为弃风弃光量；Ct0,s为仿真初始年限储能装机容量；kt,s为储能变化率.

2.2.6能源消耗和碳排放模块 能源消耗和碳排放子模块的动力学演化模型如图8所示.能源利用效率高低直接通过发电煤耗体现，发电能源结构和发电技术又是造成发电煤耗过高的原因，煤耗降低对碳排放起到直接抑制作用.

图8 能源消耗和电力碳排放模块Fig.8 Energy consumption and power carbon emission modules

Qt,fos=Qt,coalλc-f

(31)

Ct,fos=Qt,fosSt,foskfos-c

(32)

(33)

(34)

Pt,rawQt,raw=Pt,staQt,sta

(35)

式中：Qt,fos为化石能源消费量，与燃煤消耗量的关系由(31)给出，其直接影响化石能源碳排放量Ct,fos，如式(32)所示.λc-f为燃煤消耗量与化石能源消耗量的转换系数；St,fos、kfos-c分别为化石能源对标准煤的折算系数以及化石能源的碳排放系数.电力系统碳排放主要由化石能源碳排放决定，式(33)给出电力系统二氧化碳排放量的计算式，Ct0为仿真初始年限碳排放量.Ct,en为碳排放强度，指每单位国民生产总值所带来的二氧化碳排放量[16]，反映碳排放量和经济之间的关系.本文中能源价格指原煤价格，Pt,raw、Pt,sta表示原煤和标煤的价格；Qt,raw、Qt,sta表示原煤量和标煤量，《中华人民共和国价格法释义》[17]中规定计划电煤采用政府指导价格，在社会资源一定的情况下，原煤价格上涨时，标煤量会减少，关系如式(35)所示，会间接减少碳排放，进而影响碳排放强度.

3 算例仿真

3.1 算例基础数据

随着电力系统复杂性和不确定性的增强，系统环境日益复杂，本文采用情景预测分析方法[18-19]对未来发展情况进行分析和预测.通过分析系统特定因素之间的互相联系，产生可能出现的情景模式.

为全面模拟我国电力系统转型过程中可能出现的碳排放路径，本文设置基准情景、低碳情景、超低碳情景3种情景进行模拟，具体情况将在后文中详细叙述.

为验证前文所构建模型在各种情景下的可行性和有效性，调研中国电力企业协会数据和中国统计年鉴得到关键参数，如表1所示.调用Vensim软件，仿真分析投资力度、能源侧结构变化、关键技术发展等因素对发电结构、碳排放量的影响.仿真初始时间设为2020年1月，结束时间为2030年12月，共120个月，时间步长为1个月.

表1 参数初始值Tab.1 Initial value of parameters

3.2 仿真结果及分析

3.2.1基准情景 基准情景考虑的发展驱动是投资力度.未来十年我国新能源新增装机容量要达到7×108kW，平均投资变化率为9.15%.人口、全社会用电量、电力行业技术等指标延续历史发展的自然规律，不考虑碳排放对行业发展的约束，政府通过现行的环境排放标准版和有限的经济调控手段引导电力行业投资，仿真结果如图9所示.

图9 基准情景下电力系统指标变化Fig.9 Changes of power system index under base scenarios

仿真结果表明即使政府大力投资新能源，由于新能源发电技术成本较高，累计装机容量、装机占比、发电量比例增长缓慢.除此之外，为满足电力需求，火电企业燃煤发电量逐年增加，电力系统碳排放近似呈现线性增长趋势.

3.2.2低碳情景 低碳情景考虑我国环境与经济的协调、低碳可持续发展等因素.在基准情景的投资力度下，分别考虑发电技术发展(低碳情景1)和市场监督(低碳情景2)对新能源消纳量的影响，进而影响电力系统转型过程中的碳排放.

低碳情景1主要考虑火电机组的灵活性改造和储能规模化发展对弃新能源的影响，仿真结果如图10所示.

图10 低碳情景1下电力系统指标变化Fig.10 Changes in power system indicators under low-carbon scenarios 1

绿证市场是我国电力市场的发展趋势，绿电的交易不仅会影响新能源厂商的收入，间接影响装机容量，也是促进消纳新能源的重要方式之一，低碳情景2的仿真结果如图11所示.

图11 低碳情景2下电力系统指标变化Fig.11 Changes in power system indicators under low-carbon scenarios 2

在基准情景下，提倡低耗能和节能减排，火电机组灵活性改造和市场消纳机制投入示范运行或进入商业应用，带动发电效率提升.由于低碳和节能减排政策的影响，绿证市场快速发展，新能源机组逐渐增多，燃煤机组的增长速度受到约束，燃煤机组发电量略受影响，燃煤机组仍然是主要发电形式.此种情况下，弃新能源的现象未得到解决，新能源的发电占比增长相对缓慢.

3.2.3超低碳情景 超低碳情景主要考虑电力新技术的进步和市场消纳机制的协同作用.强化电力行业技术进步，建设坚强电网，协调新能源并网问题，新能源得到快速发展，增加低碳调度运行方式的灵活性.结合基准情景、低碳情景，得到图12所示仿真结果对比图.

图12 3种情景电力系统主要指标对比Fig.12 Comparison of main indexes of power systems under three scenarios

基于不同发展路径，图12给出未来十年电力系统主要指标的演化趋势.通过绿证市场交易绿电，绿电厂商收益增加，促进新能源装机增加，间接影响燃煤机组装机，低碳情景2对新能源装机容量的促进作用最明显.低碳情景1主要考虑灵活性改造的技术发展，“新能源+火电机组”的协调发展得到定量分析，灵活性改造降低火电机组的出力下限，利于火电机组参与辅助服务，促进新能源机组出力，新能源发电占比显著增加.燃煤机组发电逐步被新能源发电替代，其发电量也随之逐渐减少.大约在2027年，新能源发电量显著增加，超低碳情景对电力系统新能源的消纳优势更加明显，多种消纳方式协同作用能增强电力系统的稳定性.从碳排放变化角度看，除基准情景外，所模拟的情景均对碳排放均有一定的抑制作用.在仿真年限内，低碳情景1大约在2027年达到碳排放变化的最大值，超低碳情景大约在2028年达到碳排放变化的最大值，低碳情景2的碳排放变化最慢，但呈现持续增长的趋势.

4 结语

基于系统动力学，构建了双碳目标下电力系统供给侧形态发展的推演模型.政策投资力度是电力系统转型的基础，清洁能源消纳是关键.通过仿真模拟，验证了模型的可行性和有效性.结果表明：在电力系统供给侧转型初期，发电新技术、市场消纳能力都会限制能源发电结构的发展，此时投资力度对新能源累计装机容量、装机占比及发电量的促进作用不明显，碳排放近似线性增长.随着消纳技术和消纳市场的不断成熟，火电机组灵活性增强、储能规模化发展、绿证市场趋于成熟，对新能源机组出力促进作用增强，新能源发电成为主要发电形式，有效降低电力碳排放，助力双碳目标的实现.