《1. 引言》
1. 引言
二氧化碳(CO2 )是电厂排放的主要温室气体。基于胺溶液的燃烧后碳捕集由于可以通过对传统火电厂的简单改造而实现，是一种前景广阔的大范围碳捕集技术[1]。而类属伯胺的乙醇胺(MEA)则是采用溶液类碳捕集法最可行的策略，因为相比仲胺和叔胺，它与二氧化碳的反应速度更快[2]。之前的研究聚焦在特定碳捕集水平下，碳捕集过程的优化运行[1,3–7]。然而，乙醇胺溶液的再生能耗高、价格贵，而且碳配额市场下每季度拍卖的结算价格又会变化。固定碳捕集过程中碳捕集水平的运行方式并不经济。文献[8,9]已经对不同二氧化碳排放价格条件下需对碳捕集水平进行优化有所阐释。不过这些二氧化碳排放的定价机制类似碳税[10]。对于灵活的市场导向的碳配额贸易机制[11–13]，决策者需确定向碳排放市场提出竞标的碳配额数量。因此，我们有必要为配有碳捕集设备的电厂设计一种统一的竞标和运行策略来最大化电厂整个寿命周期中的利润。我们运用了Sarsa时间差分算法来为装配碳捕集设备的特定电厂的决策者寻找竞标和运行策略。竞标与运行的关系是通过随时间变化的碳配额市场下决策者的持有账户建立的[14]。对于寻找到的策略的优劣，可由获得的贴现累计利润来评估，即电厂的贴现现金流[9]。本文将按如下顺序介绍。在第2节，在考虑到各季度碳配额拍卖市场的前提下，我们会建立一个应用了碳捕集技术的电厂的利润模型，并依此讨论、提出基于溶液进行碳捕集的燃煤电厂的竞标和运行问题。在第3节，Sarsa时间差分算法会被引入并用于寻找前述系统的最优解。在第4节，通过案例，我们会解释Sarsa时间差分算法可找到一个统一的竞标和运行策略，且该策略可以最大化特定电厂的利润。本文末尾会给出相应结论。
《2. 问题形成》
2. 问题形成
在本节，我们将建立一个配有碳捕集过程的燃煤电厂利润模型以及一个简化后的温室气体排放贸易系统。此后，在碳排放贸易系统下，我们会给出电厂寿命周期中基于贴现累计利润的目标函数。
《2.1. 基于乙醇胺的碳捕集模型的建立》
2.1. 基于乙醇胺的碳捕集模型的建立
基于乙醇胺的碳捕集过程模型可由Aspen Plus® 构建[15]。模型的物理特性由电解质非随机双液法(electrolytenon-random two liquid，eNRTL)计算。模型本身由来自碳捕集试验场的实验数据验证[16]。验证后的模型将按比例放大，匹配额定容量650 MW的亚临界燃煤电厂排放的废气。图1是基于乙醇胺的碳捕集过程的流程图[6,17]。图中共有两座吸收塔(Absorber)用于废气中二氧化碳的吸收[18]。表1列出了吸收塔和汽提塔(Stripper)的参数。贫碳MEA溶液由分流器Splitter 1分为相等的两个部分，随后分别送入两座吸收塔的顶部。同时，来自电厂的废气由分流器Splitter 2等分并注入吸收塔底部。在吸收塔中，废气中的二氧化碳将自发地与乙醇胺溶液反应。吸收了大部分二氧化碳的气相会排入空气而富含碳的液相则被抽送至换热器，之后输送至汽提塔内。在汽提塔中，二氧化碳从富碳MEA溶液中分解出来，同时重新生成贫碳MEA溶液，离开汽提塔底部。为了使重新生成的贫碳MEA溶液的温度满足两座吸收塔进液口特定的温度目标，新生成的溶液需由换热器及下游的冷却器冷却。此外，在循环前，流失的MEA和水分由混合器Mixer 3补充。最后，贫碳MEA溶液会被送回至吸收塔，用于持续的二氧化碳吸收。通过汽提塔的冷凝器(Condenser)，分解出的高浓度二氧化碳产品可被压缩和运输。
图1中，基于乙醇胺的燃烧后碳捕集过程受控于四个控制回路。类似的控制方案在文献[3,17]中已有讨论。相应地，在Aspen Plus® 建立的稳态模型中，我们提出以下设定：①通过改变冷凝器热负荷，汽提塔最上层的温度设定在35℃；②通过改变冷却器热负荷，吸收塔顶层贫碳MEA溶液的温度设定在40℃；③贫碳MEA溶液的贫碳胺荷载比(lean loading)设置在
左右[18,19]；④通过改变贫碳MEA溶液的流速，碳捕集水平可以设置为来自离散值集合{50%，60%，70%，80%，90%}中的某一特定数值。值得注意的是，实际中的leanloading测量困难，重沸器(reboiler)的温度会由其热负荷决定并以此体现lean loading的大致情况。
《图1》
图1. 燃烧后碳捕集过程[6,17]。
《表1》
表1 吸收塔和汽提塔参数
在给出了吸收塔、汽提塔的规格(表1)及废气、贫碳MEA溶液的基本输入设定(表2)后，我们进一步改变贫碳MEA溶液的摩尔流速和lean loading，意图在最低重沸器热负荷下，实现特定的碳捕集水平。表3摘录了各个设定下的运行表现，从中我们可以确定在第t季度不同碳捕集水平ct 对应的最优的重沸器热负荷Qreb (ct )。
《表2》
表2 燃烧后碳捕集物质流
《表3》
表3 不同设定下燃烧后碳捕集过程的运行性能
《2.2. 亚临界燃煤电厂的利润模型》
2.2. 亚临界燃煤电厂的利润模型
燃煤电厂消耗了全世界范围内绝大部分能源并且在所有电力生产系统中排放了最多的二氧化碳，因此十分重要[20]。这一节中，我们为配有碳捕集的亚临界燃煤电厂构建了利润模型。我们也会同时建立电厂运行设定与电力输出之间的关系。依据美国能源情报署(EIA)的报告[21]，电厂第t季度的成本Ct (USD·qtr–1 )可由如下公式计算：
其中，FOM为每一季度固定的运行和维修(OM)成本；VOMt 为每一季度变化的运行和维修成本；Ft 为每一季度的燃料成本；Bt 为每一季度碳配额竞标成本。根据加利福尼亚州碳配额拍卖机制[11]，竞标成本可定义为
其中，vt 为第t季度拍卖的碳配额结算价格(USD·allowance–1 )；wt 为决策者第t季度通过拍卖赢取的碳配额吨数(t·allowance–1 )。配额1 t将允许电厂排放1 t的二氧化碳。式(1)的其他项则定义如下：
其中，β和δ分别为固定和变化的运行维修成本系数；f为单位燃料成本；Pn 为电厂额定容量。这些变量及其单位都列在表4中。此外，Et 是电力输出(kW·h·qtr–1 )，Ht 是燃料消耗(GJ·qtr–1 )。燃煤电厂电力生产的收益可以表示为
其中，λ为电价[USD·(kW·h)–1 ]并在表4中给出。之后，第t季度的利润可作如下推导:
式中，利润Pt 是随变量Et ，Ht ，wt 和vt 变化的，所以我们将其记作Pt =P(Et , Ht , wt , vt )。
《表4》
表4 配有碳捕集的电厂的参数[21]
在本文中，电力生产和碳捕集的总燃料消耗Ht 假定为常量。据此，主电厂的电力输出和碳捕集设施的捕集能耗是要相互权衡的。对于如表4所设定的额定容量为Pn =650 000 kW的燃煤电厂，1个季度的燃料消耗可按下式计算：
式中的“2190”代表1个季度的小时数。每一季度的碳捕集能耗Ut (GJ·qtr–1 )需由下式约束：
Ut 应是第2.1节中讨论的重沸器热负荷的函数，即
结合式(9)和(10)，第t季度的电力输出Et 可以推导为
以上算式表明碳捕集水平ct 可以唯一确定电力输出Et 。电厂每一季度的利润[公式(7)]可以简化为
总体说来，对于某一配有碳捕集技术的特定电厂，假定燃料消耗H固定以及燃料和电力价格(f和l)持续不变，电厂第t季度的利润Pt 可以由碳捕集水平ct 、碳拍卖市场中获得的碳配额数量wt 和单位碳配额的结算价格vt唯一确定。
《2.3. 碳配额市场》
2.3. 碳配额市场
在第2.2节，虽然一个季度的利润P t 已经被完全定义，但我们只讨论了Et 和Ht 两个自由度。另外两个自由度的变量wt 和vt 会受到碳配额市场条件的影响。只有当市场包含的或选择加入的所有实体(如发电公司)向拍卖操作员上交了竞标选项(包括竞标数量q和竞标价格p)后，每一季度的市场条件才能被完全定义。决策者相关的实体的竞标数量和竞标价格分别记为q0,t 和p0,t ；所有其他实体的竞标数量和竞标价格分别记为qi,t 和pi,t 。其中，i∈I，I={1, 2, 3, …, I}，是除去决策者实体外的、包含在碳配额市场下的实体的集合。操作员之后会实行如下的密封竞标拍卖机制[14]。
在一个季度中，若某一实体或竞标人的竞标违反了购买限制、持有限制或竞标担保的相关规定，拍卖操作员会因其不合格而拒绝该竞标。接下来，会依据竞标价格按降序考虑所有竞标人的合格的竞标方案。从出价最高的竞标开始，各个竞标价格下呈递的竞标方案都会依竞标人的竞标数量出售等量的碳配额，直到满足如下任意一条件：一是所有碳交易市场下的碳配额A(单位：allowance)售完；二是下一个竞标人的竞标价格低于拍卖的竞标底价g t (单位：USD·allowance–1 )[11]。如果拍卖的碳配额售完，则结算价格是出售了碳配额的最后一个竞标方案的竞标价格；如果结算价格等于拍卖底价，则出售的碳配额应该等于所有在拍卖底价之上的各竞标方案共同累计的竞标数量。拍卖操作员之后可以计算出每一个竞标人或实体赢得的碳配额(如决策者赢得的碳配额在本文中记为wt )、总共售出的碳配额ut 和对应于所有实体的一个统一结算价格vt ，即
在式(13)、式(14)或式(15)中，决策者只能决定他自身的竞标数量q0,t 和竞标价格p0,t 。但利用其他实体的历史竞标数据，决策者可以对其他公司的竞标选项进行估计。例如，在第2.4节，其他实体的各个竞标选项是由概率来表示的。
在本文中，为了判断一个竞标方案是否合格，我们只考虑电厂决策者的持有限制hl (单位：allowance)。为了简单，不再考虑购买限制和竞标担保。持有限制是包含在碳配额市场下某一实体账户的持有上限。如果任意实体的竞标数量会潜在地导致其持有账户里的碳配额超过持有限制，那么这一呈递的竞标方案将会因不合格而受到拍卖操作员的拒绝。本文中的持有账户只拥有加利福尼亚州出台的法规中多个账户的类似功能[11]。我们假设持有账户的碳配额有以下约束：
其中，ht+1 为第t+1季度拍卖开始时持有账户的碳配额。值得注意的是，如果赢得的碳配额wt 低于电厂碳排放et ，额外的碳配额将从持有账户中上缴；如果赢得的碳配额高于碳排放，则剩余的赢取的碳配额将保留至持有账户。在第t季度前的所有季度中，累积在持有账户的总碳配额被记为ht 。式(17)的不等式表明持有账户的碳配额不能耗尽，否则，没有缴纳碳配额的过度排放将产生额外的罚金。根据加利福尼亚州法规[11]，对于某实体产生二氧化碳排放并未及时上缴碳配额的那部分，作为惩罚，新的遵守义务将是上缴更多的碳配额，其数量四倍于本应上缴的碳配额数量。因未及时上缴而导致的该义务的履行需要额外的竞标和贸易机制。为了简单，不同于加利福尼亚州法规中规定的惩罚，我们假定未及时上缴时，每吨过度排放二氧化碳需交的罚金为320 USD。因此，式(17)可以认为是一个软边界。另外，式(16)表明对于任意t，决策者只能呈递不会潜在地导致ht+1 超过之前提到的hl 的竞标数量q0,t 。变量et 代表电厂第t季度的碳排放(单位：t·qtr–1 )，由下式决定：
ct 碳捕集水平是基于溶液的碳捕集过程的运行参数，而wt 和q0,t 是在碳配额市场下的竞标参数。式(16)和式(17)中的不等式表明了竞标和运行的潜在关系。
《2.4. 目标构建》
2.4. 目标构建
在第2.2节，利润[式(12)]可由碳捕集水平ct 、决策者赢得的碳配额数量wt 和结算价格vt 表示。碳捕集水平ct 可以由决策者任意确定，而赢得的碳配额数量wt 和结算价格vt 必须由式(13)和式(14)中出现的所有实体的竞标选项决定。如果所有其他实体已经呈递了他们的竞标选项(如任意
下的pi,t 和qi,t )，配有碳捕集电厂的决策者只需决定自己运行方案ct 和竞标方案(q0,t , p0,t )，就能估计式(12)中相应的利润。可以定义如下统一的决策者的行为是
其中，
是在状态st 下的一个离散的行为集合，并且进一步认为是不随状态变化的集合A。决策者的电厂只知道自己的竞标数量q0,t 和价格p0,t ；任意下的其余竞标人的qi,t 和pi,t 必须由决策者通过先验知识获得。在本文中，其他实体的竞标数量和价格被提前假定为受到上一季度拍卖市场结算价格、出售碳配额的影响。一个相似的状态选择方案在文献[22]中有所讨论。因此，第t季度的状态st 可按下式记为
其中，因为持有账户的碳配额ht 必须充足[式(17)]但不能潜在地超过持有上线hl [式(16)]，所以被认为是一个状态变量。此外，我们趋向于最大化所关心的电厂寿命下的贴现累计利润。因此，时间t也作为式(20)中st 的一个状态。不同时间状态表示的是电厂寿命的不同阶段，在各个阶段决策者可以设定不同的行为方案。若实体的竞标数量集合和竞标价格集合分为
和
则决策者可估计任意竞标人选择某竞标选项的概率
需对任意 成立。虽然每个季度中实体会选择自己的竞标选项，对于数量和价格的竞标选项的集合(即和)在本文中假定为时不变且不随状态而变化。之后，我们建立以下马尔可夫决策过程。在特定的状态st =s下，决策者会执行行为at =a。在定义了联合概率
后，所有其他竞标人将如式(22)所设定的那样选择他们自己的竞标选项，使得下一季度的状态st+1 = (vt , ut ,ht+1 , t+1) = s′ 在执行行为at = a后可以被唯一确定。进一步地，奖励值函数rt+1 是由状态st 过渡至st+1 时基于式(12)导出，即
注意到“奖励值”是在强化学习(RL)框架下定义的专业术语。从物理意义上来说，第t+1季度的奖励值rt+1就是电厂在第t季度的利润Pt [式(12)]。因为t可以是任意时刻，决策者可以递归地获得有限时间区间内的奖励值序列st , at , st+1 , rt+1 , at+1 , st+2 , rt+2 , …, at+N–1 , st+N , rt+N ，称为一段竞标和运行事件。变量N表示配有基于乙醇胺溶液碳捕集过程的电厂的寿命。对于∀k∈{0, 1, …, N–1}，目标函数可以构建为
并服从
其中，Vp(st )为在方针p下关于s t 的状态值函数；rt+k+1 为由状态s t+k 过渡至s t+k+1 的奖励值；γ为贴现率；Eπ {·}为贴现奖励值序列在策略p下的期望。对于决策者，一个随机方针的概率可以写为p(st+k , at+k )，其中，在各个状态下任意行为出现的概率应该事先由决策者确定，目的是最大化电厂寿命下的贴现累计利润。我们先考虑随机方针是因为最优方针是由基于强化学习的Sarsa算法搜索得到的。最后，随机方针将逐渐演变为一个可应用的确定优化方针。此外，式(26)和式(27)可由式(16)和式(17)分别获得。
《3. Sarsa 算法：介绍与应用》
3. Sarsa 算法：介绍与应用
基于强化学习的Sarsa时间差分算法可在Maltlab® 中实现，找到在第2节中所定义问题的最优方针或策略。由于该算法的自适应及无模型特性，它可使电厂利润最大化。一个最初的优化方针可以根据第2节中的仿真环境自行找到；进一步的策略调整可以在agent与真实环境互动后实现。Sarsa时间差分算法比起动态规划法计算时间更少，而比起另一种基本的强化学习法(Q-学习法)[23]，收敛性要更好。然而，如果调节参数值(如ε)的计划给定序列不恰当，Sarsa时间差分算法经常会找到一些较差的方针。参数ε是在行为集合A内的等可能寻找行为的概率。在ε-贪心方针下，各个状态下的应对方针都包含两种模式。一种是直接选择该状态下的贪心行为。该贪心行为是当前状态下决策者认为最好的行为，且无穷迭代后将收敛至最优行为。选择该行为模式的概率是1–ε。另一种是为了探索当前状态下各个行为的优劣，即是否存在更好的行为。选择该模式的概率为ε。若决策者随机地进入了探索模式，行为集合A内各个行为都有被采用的可能，各行为被采用的概率都是ε/nA 。nA 是集合A包括的行为的数量，将在式(32)给出。
为了设计Sarsa时间差分算法，我们根据式(24)定义了一个最优行为值函数，即
该式对于所有
成立。其中，Qn 为在方针n下的行为值函数。以此定义的最优方针为
然而，在最优方针并未获取的前提下，最优的行为值Q* 是未知的。根据文献[23,24]，行为值函数的迭代方案[式(30)]可以确保：行为值函数Qt+k+1 (s, a)在k → ∞且无穷次访问所有状态
和所有行为
后，将收敛至Qπ (s, a)。该迭代方案记为
其中，a为当前方针π下的行为。如果根据式(30)得到的Qπ (s, a)的估计值为
改进方针可通过下式实现：
其中，
为使方针 π' （s ，a） 比π(s, a)有所提升的贪心行为；
为行为集合
中包含的行为的数量；ε为从行为集合中均匀选取任意行为的概率。除了贪心行为，其他所有行为称为探索行为。探索行为可以保证找到每个状态对应的全局最大状态值Vπ [式(24)]而非局部最大值。通过设定新的方针π←π′，式(30)至式(32)构成了一组可重复迭代无穷多次的行为值迭代算法。这一算法可以用来获取最优方针π* 。但使用前，a和ε都需设置计划值序列。学习率a起初要很大来实现对Q(s, a)[式(30)]的快速初始化，但最终减小以确保行为值的收敛。虽然a的计划值序列存在理论条件，但它们很少被实际应用[23]。探索行为集合的概率ε在探索初期等于1，但会逐渐降低至0，因为我们最终将推出一个确定方针。表5给出应用了ε-贪心方针的Sarsa时间差分算法的伪代码。
《表5》
表5 应用ε-贪心方针的Sarsa时间差分算法
Sarsa时间差分算法是无模型的在线算法，因此可以通过与环境(即真实碳配额市场)的交互来实现。然而，在本文中，对于其他实体的竞标选项及相应概率的估计是必要的，这些可以用来形成一个模拟碳配额市场。这样的先验知识可由其他电厂的历史竞标数据来获得。如果历史竞标数据无从获取，通过统计分析，历史的市场条件仍可以用来识别状态之间相互过渡的概率[22]。在此基础上，我们可以得到一个基于强化学习的Sarsa时间差分算法(表5)的初始方针。这样的好处是与真实碳拍卖市场的交互次数减少。文献[23]给出了一个统一的规划和学习视角，同时考虑了仿真模型和真实环境。
《4. 结果与讨论》
4. 结果与讨论
在案例研究中，共有8个包含在市场下的实体，标记为0，1，2，3，4，5，6，7。实体0是运行燃煤电厂的决策者，其设定如表4所示。这一实体被假定为我们自己的公司，并且想要最大化该公司下一家电厂的贴现累计利润。决策者将利用Sarsa时间差分算法来寻找各个状态下合适的竞标和运行行为。实体1运行着一家与实体0在设定上完全相同的电厂，但采用了不同的竞标和运行策略。该策略将会在第4.3节介绍。其他实体(即实体2~7)的竞标策略是事先定义且被假定由决策者的仿真环境来预先估计。对于式(24)中的目标函数，初始时刻是t=0，而相应的时间范围是N=100季度(即电厂的寿命是25年)，因此k的取值范围是k∈{0, 1, 2, …,99}。由以上定义推知，任意随时间变化的量可由k来索引，如sk ，ak ，rk+1 。拥有25年寿命的电厂的年度贴现率设置为8%[9]，所以季度制的贴现率为γ = 1/(1 +8%) 0.25 ≈ 0.98。持有限制h l 的制定需基于年度碳配额预算[11]。然而，年度碳配额预算是计划制定的，且可能每年各不相同。为了简单起见，本文认为hl 是一个常数，即6×106 allowances。在表5中，μ是8段事件，a由1/20逐渐变至1/200，ε由1逐渐变至0.1。变量a和ε都在执行方针改进后改变。
式(20)的状态变量需要聚合为离散的子集来缓解状态空间的维数灾难；该方法又被称为状态聚合[22,23]。状态聚合可由如下方式实现：结算价格和出售碳配额将一起被考虑，因为当一个处于某个特定的区间时，另一个也应被约束在某个特定值。例如，当碳配额拍卖中的出售碳配额uk-1 小于总体拍卖的碳配额A = 1 500 000 t，则结算价格必须等于竞标底价g，即表6中is =1，2，3对应行的状态。同样地，时刻k和持有账户的碳配额hk 分别聚合并记录在表6、表7中。根据这两个表格，原始的状态空间
被离散为8 × 5 × 14 = 560个聚合状态。行为变量[式(19)]则归为两部分。其一是运行部分，包括5种决策者燃煤电厂可设定的碳捕集水平，即表3中的
={50%, 60%, 70%, 80%, 90%}；其二是16种可行的竞标选项，包括竞标数量和价格，即(q0 , p0 )∈
。类似于其他实体的竞标数量和竞标价格集合(即
)，我们只考虑时不变的、不随状态变化的竞标选项集 。由此，对于决策者的各个聚合状态，共有5 × 16 = 80种不同的行为方式。接下来的内容里，我们将在决策者实际采用了某个特定行为时提及该行为的具体值。为简便起见，由和 集合给出的80中不同行为并未在本文列出。
《表6 》
表6 结算价格和售出碳配额对(vk-1 , uk-1 )的子集及时刻k的子集分布
《表7》
表7 持有账户碳配额的子集分布
《4.1. Sarsa 时间差分算法的收敛性》
4.1. Sarsa 时间差分算法的收敛性
在这一节，我们提供了特定状态某个行为的行为值的收敛特性。由于状态都已聚合，我们只考虑在表6和表7中各个状态聚合后的子集，而非实际任意k下的状态值sk 。图2展示了某一状态行为对(s, a)的行为值Q(s, a)的收敛情况。其中，状态s被归为元组(is , js , vs ) = (7, 7,4)，而行为a的索引为ia =61。索引ia =61对应的行为是a =(300 000, 14.5, 27)。该行为事先设定在离散行为集合
中。这一状态行为对应的Q(s, a)通过式(30)的迭代将收敛至式(29)中的最优行为值Q*(s, a)。由于表5中的算法可以访问到所有的状态行为对(s, a)，因此对于任意
和
，表5的迭代会得到对应的Q*(s, a)的值，从而各个状态的最优行为可由式(29)确定。如前所述，每个状态都有如图3所示的80种行为值。基于这些行为值，我们可以得出行为索引是ia = 61的行为有最大的Q值，所以它给出了在该状态下最好的行为方式。最优方针因此可以通过寻找每个聚合状态下的最大行为值找到。
《图2》
图2. 典型状态行为对在is = 5，js = 7，vs = 4，ia = 6处的收敛性。
《图3》
图3. 在特定状态is = 5，js = 7，vs = 4处所有行为的行为值。
《4.2. Sarsa 时间差分算法的性能》
4.2. Sarsa 时间差分算法的性能
在这一节，通过与大多数相关文献中设定固定碳捕集水平的运行方式对比，我们展示了考虑到随时间变化碳捕集水平的Sarsa时间差分算法可以获取更多的贴现累计利润。k=0时刻的初始竞标底价为12.73 USD·allowance–1 [11]。另外，年度竞标底价的增长率τ被引入，用来提高每年的竞标底价。年度竞标底价增长率可以刺激新兴碳捕集和埋存技术的进步。某一结算价格的案例已显示在图4中。从中不仅可以看到结算价格在整个时间范围内(即100个季度)的波动，而且由于5%的年度增长率，结算价格出现了计划增长。该年度增长率同样设置于加利福尼亚和魁北克的联合温室气体拍卖机制中[11,14]。
《图4》
图4. 在年增长率τ = 5%和初始竞标底价g =12.73 USD·allowance–1 下的结算价格。
为了展现Sarsa时间差分算法的自适应性，计划的碳配额竞标底价的年度增长率τ被分别假定为0%、5%、10%和15%。如果特定的年度增长率固定在如图5所示的τ = 0%，则除了竞争对手(competitor，即Entity1)所示曲线，四组不同的奖励值序列表明了实体0的决策者采用四组不同竞标和运行策略的获利情况。其中一组的竞标运行策略是用考虑到随时间变化碳捕集水平的Sarsa时间差分算法找到的。对于其余策略：运行方面，固定了碳捕集水平(即碳捕集水平在相关的一段竞标和运行事件中被设置为50%、70%或90%)；竞标方面，事先定义了各聚合状态下所有行为中各竞标选项出现的概率。对于固定了碳捕集水平的策略，可能的竞标选项来自于竞标选项集合
这一集合也是基于Sarsa的统一的竞标和运行策略集合。之前的奖励值序列显示了以季度为基础的电厂寿命下的利润。通过计算特定奖励值序列的贴现和，我们可以获得某个特定竞标和运行策略的贴现累计利润。例如，根据图5，我们得出在年度竞标底价增长率τ为0时各个策略对应的贴现累计利润。
《图5》
图5. 在年增长率τ = 0%和持有账户初始碳配额h 0 = 0.05 × 106 allowances下不同竞标和运行策略的奖励值。
类似地，如图6至图8所示，我们可以得到在持有账户初始碳配额为h0 = 0.05×106 和τ由5%变化至15%时的贴现累计利润。在持有账户特定初始碳配额为h0 =0.05×106 t的前提下，不同竞标底价增长率对应的贴现累计利润如图9所示。
《图6》
图6. 在年增长率τ = 5%和持有账户初始碳配额h 0 = 0.05 × 106 allowances下不同竞标和运行策略的奖励值。
《图7》
图7. 在年增长率τ = 10%和持有账户初始碳配额h 0 = 0.05 × 106allowances下不同竞标和运行策略的奖励值。
《图8》
图8. 在年增长率τ = 15%和持有账户初始碳配额h 0 = 0.05 × 106allowances下不同竞标和运行策略的奖励值。
《图9》
图9. 在持有账户初始碳配额h0 = 0.05 × 106allowances下的贴现累计利润。
对于以其他持有账户初始碳配额数量为前提的贴现累计利润展示在图10和图11中。可以推知，无论采用基于固定碳捕集水平方案的固定碳捕集水平设置为何值，Sarsa时间差分法找到的统一、灵活的运行和竞标策略都表现得更为出色。
《图10》
图10. 在持有账户初始碳配额h0 = 3 × 106allowances下的贴现累计利润。
《图11》
图11. 在持有账户初始碳配额h0 = 5 × 106allowances下的贴现累计利润。
《4.3. 与碳配额市场下另一实体的比较》
4.3. 与碳配额市场下另一实体的比较
在同一碳配额市场下，我们将比较采用了Sarsa时间差分算法的决策者和其竞争对手(实体1)的获利情况。对于该竞争对手，其电厂的所有设置都假定与实体0的决策者的相同。而在关于运行和竞标策略的问题上，实体1将其碳捕集水平固定在60%，而它的竞标选项则从与实体0选用的相同集合
中独立选取。我们假设实体1在竞标选项的选择行为上可由实体0的决策者根据玻尔兹曼分布近似，即
其中，y和z为可选的竞标选项的索引；nb 为总的竞标选项的数量(对于竞标选项集 ，nb =16)；Pr(y)记为选择以y为索引的竞标选项的概率；ξ为分布的温度常数。从式(33)中可以看出，ξ很大时，所有可能的竞标选项的选择接近于等概率。为了简单，本例中的温度常数为ξ=1。变量ω代表由y或z索引的各个选项的权重。在本文有关仿真中，所有权重中最大的被定义为常值，即ωmax = nb 。假定第13个竞标选项被赋予最大的权重，即ω(y = 13) =ωmax = 16，则以索引y = 13为中心，所有竞标选项的权重依据离该中心的远近而指定为：4，5，6，7，8，9，10，11，12，13，14，15，16，15，14，13。根据这些权重，选择某个可能竞标选项的概率可以由式(33)事先定义。实际中，决策者需获得该竞争对手的历史竞标数据或历史市场条件来识别这些权重。
决策者持有账户的碳配额hk 和竞争对手的碳配额数量分别如图12和图13所示。两个实体的贴现累计利润则可根据图5和图6得到并记录在图14中。图中，决策者在不同竞标底价的年度增长率下，都获得了更多的贴现累计利润。这些表明，决策者通过Sarsa时间差分法得到的竞标和运行策略要比处在同一碳配额市场下的竞争对手的策略更好。
《图12》
图12. 不同竞标底价下采用Sarsa时间差分策略的决策者持有账户的碳配额数量。
《图13》
图13. 不同竞标底价下竞争对手持有账户的碳配额数量。
《图14》
图14. 在持有账户初始碳配额h0 = 0.05 × 106allowances下实体0的决策者和实体1的贴现累计利润。
《5. 结论》
5. 结论
本文为燃煤电厂设计了由Sarsa时间差分算法得到的统一的竞标和运行策略。当考虑了来自碳捕集水平集合的随时间变化的灵活捕集水平和一个共同设计的竞标选项集合后，提出的策略要优于运行上固定碳捕集水平、竞标上独立设计的策略。Sarsa时间差分算法能够最大化在不同碳配额市场条件下电厂的贴现累计利润，如竞标底价采用不同的年增长率或是持有账户拥有不同的初始碳配额。此外，通过与另一个固定其碳捕集水平并以玻尔兹曼随机分布独立设计竞标策略的电厂比较，运用了Sarsa时间差分算法的决策者在碳配额市场下更具竞争力。
《Compliance with ethics guidelines》
Compliance with ethics guidelines
Ziang Li, Zhengtao Ding, and Meihong Wang declare that they have no conflict of interest or financial conflicts to disclose.
京公网安备 11010502051620号
