1.本发明公开了一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，属于电力工程研究区域，尤其涉及基于考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法。
背景技术：2.我国西南地区与东南亚国家边境地区人口密度高，贸易和能源需求增长迅速，跨境贸易、旅游、加工业、农业等方面合作发展迅速，未来发展潜力巨大。能源互联互通是中国对外合作的重要组成部分，跨境能源合作已成为助推经济快速发展的重要环节，对于实现能源合作利益和命运共同，保障各国能源安全起着至关重要的作用。然而，能源互联互通作为中外合作的重要组成部分，面临着法律、政策、金融层面的跨境市场交易风险、双边调度机制差异较大、信息交流有限等问题。不完善的制度建设和技术手段难以保障跨境综合能源系统安全稳定运行，此外，还可能降低跨境各方参与积极性，从而影响跨境资源的优化配置和整体福利的提升。
3.现有的能源跨境合作以单一能源形式为主，多为电力及天然气管线的点对点传输，建设周期长、通道利用率低，缺乏通道双向供给能力、新能源就地消纳能力及资源互补利用能力。
技术实现要素：4.本发明针对我国与周边邻国能源市场机制不同、多边协同调度难度大等问题，提出了一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，用于确定现有能源跨境协同运行和新能源规模化利用的最优化规模；
5.一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，其特征在于，包括如下步骤：
6.s1、构建跨境综合能源系统cbies模型，所述跨境综合能源系统cbies包括燃气轮机、燃气锅炉、可再生能源机组、储能设备和电锅炉，其中可再生能源机组包括风力发电和光伏发电，储能设备包括蓄电池和蓄热槽；
7.s2、在步骤s1所述的跨境综合能源系统cbies模型基础上引入跨境绿证交易模型和碳交易模型，建立考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型；
8.s3、根据步骤s2中建立的考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型，使用求解器求解得出优化结果。
9.s11、所述步骤s1中构建的跨境综合能源系统cbies模型具体包括燃气轮机模型、燃气锅炉模型、可再生能源机组模型、储能设备模型和电锅炉模型，具体如下：
10.燃气轮机模型具体如下：
11.[0012][0013]
其中表示第i个国家综合能源系统t时刻燃气轮机的输出电功率，表示第i个国家综合能源系统t时刻燃气轮机的输出热功率，l
ng
表示天然气的低位热值，表示燃气轮机在t时刻的天然气消耗量，η
gt
表示燃气轮机的发电效率，n
gt
表示余热回收系数；
[0014]
燃气锅炉模型具体如下：
[0015][0016]
其中表示第i个国家综合能源系统t时刻燃气锅炉的输出热功率，l
ng
表示天然气的低位热值，表示燃气锅炉在t时刻的天然气消耗量，η
gb
表示燃气锅炉的发热效率；
[0017]
可再生能源机组模型包括风力发电模型和光伏发电模型，具体如下：
[0018]
风力发电模型：
[0019][0020]
其中为第i个国家综合能源系统t时刻风电功率，ρ为空气密度，a为风轮机叶片扫掠面积，v为风速，c
wt
为风能利用系数，是单位时间内风轮所吸收的风能与通过风轮旋转面的全部风能之比，λ
wt
为叶尖速率比；
[0021]
光伏发电模型：
[0022][0023]
ts＝ta+0.0138
·
(1+0.031ta)
·
(1-0.042v)
·gꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0024]
其中为第i个国家综合能源系统t时刻光伏发电功率，g为光照强度，ts为光伏电池表面温度，p
stc
、g
stc
、t
stc
分别为标准测试条件下的最大输出功率、光照强度和光伏电池表面温度，ε为光伏电池温度系数；ta为环境温度；v为风速；
[0025]
储能设备模型包括蓄电池模型和蓄热槽模型，具体如下：
[0026]
蓄电池模型：
[0027]
放电
[0028][0029]
充电
[0030][0031]
其中为第i个国家综合能源系统t+1时刻蓄电池中存储的电量，w
te,i
为第 i个国家综合能源系统t时刻蓄电池中存储的电量，和分别为第i个国家综合能源系统t时刻蓄电池放电功率和充电功率；和分别为蓄电池自身的放电效率和充电效率；和分别为蓄电池自身的放电损耗和充电损耗；
[0032]
蓄热槽模型：
[0033]
放热
[0034][0035]
充热
[0036][0037]
其中为第i个国家综合能源系统t+1时刻蓄热槽中存储的热能，w
th,i
为第i个国家综合能源系统t时刻蓄热槽中存储的热能，和分别为第i 个国家综合能源系统t时刻蓄热槽放热功率和充热功率；和分别为蓄热槽自身的放热效率和充热效率；和分别为蓄热槽自身的放热损耗和充热损耗；
[0038]
电锅炉模型，具体如下：
[0039][0040]
其中，为第i个国家综合能源系统t时刻电锅炉的制热功率；为第 i个国家综合能源系统t时刻电锅炉所需电功率，η
eb
为电锅炉的转换效率。
[0041]
s21、在步骤s1所述的跨境综合能源系统cbies模型基础上引入碳交易模型具体包括碳交易成本，具体如下：
[0042][0043]
其中为跨境综合能源系统的碳交易成本，n为综合能源系统个数，为第i个国家综合能源系统碳排放量，为第i个国家综合能源系统碳初始配额，为市场上的碳排放权价格；d为碳排放量区间长度；σ为各阶梯碳排放权价格增长幅度，每上升一个阶梯，碳排放权价格增长时，碳交易成本为负值，此时表示出售碳排放权获得收益；
[0044]
第i个国家综合能源系统碳排放量计算如下：
[0045][0046]
其中为跨境综合能源系统直接碳排放量，为外购电力的碳排放量，为外购热力的碳排放量；
[0047]
跨境综合能源系统直接碳排放量计算如下：
[0048][0049]
其中n为综合能源系统个数；t为小时数，取为24小时；是第i个国家综合能源系统中燃气轮机t时刻输出的电功率，是第i个国家综合能源系统中燃气锅炉t时刻输出的热功率；和分别为第i个国家综合能源系统中的燃气轮机发电效率和燃气锅炉发热效率；
[0050]
外购电力的碳排放量计算如下：
[0051][0052]
其中表示第i个国家综合能源系统外购电力消耗量，单位为万千瓦时；为电网基准线排放因子，表示消耗单位电力的间接排放量，本文取值 0.8367，单位为吨co2/千瓦时；
[0053]
外购热力的碳排放量计算如下：
[0054][0055]
其中表示第i个国家综合能源系统外购热力消耗量，单位为百万千焦；为热网基准线排放因子，表示消耗单位热力的间接排放量，本文取0.096，单位为吨co2/百万千焦；
[0056]
第i个国家综合能源系统碳初始配额计算如下：
[0057][0058]
其中为第i个国家综合能源系统燃气机组t时刻供电量，单位为mwh， p
e,b
为燃气机组的供电基准值，本文取值0.3791，单位为tco2/mwh，γh为燃气机组供热量修正系数，燃气机组供热量修正系数为1；为第i个国家综合能源系统燃气机组供热量，单位为gj，p
h,b
为燃气机组供热基准值，本文取值 0.0600，单位为tco2/gj；
[0059]
s22、在步骤s1所述的跨境综合能源系统cbies模型基础上引入跨境绿证交易模型具体包括跨境绿证交易成本，具体如下：
[0060][0061]
其中为第i个国家综合能源系统跨境绿证交易成本，n为综合能源系统个数，为第i个国家综合能源系统可再生能源实际消纳量，为第i个国家综合能源系统可再生能源配额量，和分别为购买和出售绿色证书的价格，λf为惩罚系数；
[0062]
第i个国家综合能源系统可再生能源配额量计算如下：
[0063]
[0064]
其中为第i个国家综合能源系统区域预测年消纳可再生能源量，为第i个国家综合能源系统区域预测年接入可再生能源量，为第i个国家综合能源系统区域预计全社会年用电量，n为综合能源系统个数；
[0065]
s23、所述步骤s2中建立考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型，具体如下：
[0066][0067]
其中β为跨境绿证交易与碳交易权折算系数，αz为交易价格调整系数，由可再生能源类型及当地燃气价格确定；为市场上的碳排放权价格，为出售绿色证书的价格；
[0068]
所述考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型包括系统运行总成本，具体如下：
[0069][0070]
其中ci、和分别为第i个国家综合能源系统运行总成本、燃料成本、购能成本、售能收益、系统维护成本，各国综合能源系统间交互成本，碳交易成本和绿证交易成本，单位元；
[0071]
第i个国家综合能源系统运行燃料成本：
[0072][0073]
其中c
fuel
为燃料费用，单位元，表示燃气轮机在t时刻的天然气消耗量，表示燃气锅炉在t时刻的天然气消耗量；
[0074]
第i个国家综合能源系统运行购能成本：
[0075][0076]
其中和分别为各国综合能源系统t时刻向电网购电费用和购热费用，单位元；和分别为第i个国家综合能源系统t时刻与电网交互的购电功率和与热网交互的购热功率，；
[0077]
第i个国家综合能源系统运行售能收益：
[0078][0079]
其中和分别为售电费用和售热费用，单位元；和分别为第i个国家综合能源系统t时刻与电网交互的售电功率和与热网交互的售热功率；
[0080]
第i个国家综合能源系统运行系统维护成本：
[0081][0082]
其中c
gt
、c
gb
、c
eb
、c
hs
、c
es
、c
pv
、c
wt
分别为燃气轮机、燃气锅炉、电锅炉、蓄热槽、蓄电池、光伏、风机的单位功率维护成本，单位元，为第i 个国家综合能源系统t时刻燃气轮
机的输出电功率，为第i个国家综合能源系统t时刻燃气锅炉的输出热功率，为第i个国家综合能源系统t时刻电锅炉的制热功率，为第i个国家综合能源系统t时刻蓄热槽放热功率或充热功率，为第i个国家综合能源系统t时刻蓄电池充电功率或放电功率，为第 i个国家综合能源系统t时刻光伏发电功率，为第i个国家综合能源系统t 时刻风电功率；
[0083]
各国综合能源系统间交互成本：
[0084][0085]
其中何分别为综合能源系统间购电何购热费用，为i、j两国综合能源系统间传输的电功率，为i、j两国综合能源系统间传输的热功率；
[0086]
所述考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型的约束条件具体包括各国综合能源系统内部功率约束、各国综合能源系统与电网交互功率约束、各国综合能源系统与热网交互功率约束、各国综合能源系统间电网交互功率约束和各国综合能源系统间热网交互功率约束具体如下：
[0087]
各国综合能源系统内部功率约束具体包括电功率平衡约束、热功率平衡约束、各国综合能源系统设备出力上下限约束、蓄电池功率约束和蓄热槽功率约束，具体如下：
[0088]
电功率平衡约束：
[0089][0090]
其中，为第i个国家综合能源系统t时刻与电网交互的购电功率，为第i个国家综合能源系统t时刻与电网交互的售电功率，为第i个国家综合能源系统t时刻燃气轮机的输出电功率，为i、j两国综合能源系统间电网传输的电功率损耗，为i、j两国综合能源系统间传输的电功率，为第i个国家综合能源系统t时刻电锅炉所需电功率，分别为第i个国家综合能源系统t时刻蓄电池充电功率、放电功率，为第i个国家综合能源系统t时刻光伏发电功率，为第i个国家综合能源系统t时刻风电功率，为第i个国家综合能源系统的电负荷；
[0091]
热功率平衡约束：
[0092][0093]
其中为第i个国家综合能源系统t时刻与热网交互的购热功率，为第i个国家综合能源系统t时刻与热网交互的售热功率,为第i个国家综合能源系统t时刻燃气锅炉的输出热功率，为第i个国家综合能源系统t时刻电锅炉的制热功率，为i、j两国综合能源系统间热网传输的热功率损耗，为i、j两国综合能源系统间传输的热功率，分别为第i个国家综合能源系统t时刻蓄热槽放热功率、充热功率，为第i
个国家的热负荷；
[0094]
各国综合能源系统设备出力上下限约束：
[0095][0096]
其中为第i个国家综合能源系统中设备m的电功率，为第i个国家综合能源系统中设备m的电功率上下限；为设备m的电功率，为设备m的电功率，为第i个国家综合能源系统中设备m的热功率上下限；
[0097]
蓄电池功率约束：
[0098][0099][0100][0101][0102]
为第i个国家综合能源系统蓄电池容量，为最大充电倍率，为最大放电倍率，分别为第i个国家综合能源系统t时刻蓄电池充电功率、放电功率；为t时刻充能的状态位，为t时刻放能的状态位，是0-1变量，表示同一设备同一时刻充放能状态唯一，为第i个国家综合能源系统蓄电池的最大、最小储能量；
[0103]
蓄热槽功率约束：
[0104][0105][0106][0107][0108]
其中为第i个国家综合能源系统蓄热槽容量，为最大充热倍率，为最大放热倍率，分别为第i个国家综合能源系统t时刻蓄热槽放热功率、充热功率，为第i个国家综合能源系统蓄热槽的最大、最小储热量；
[0109]
各国综合能源系统与电、热网交互功率约束，具体如下：
[0110]
各国综合能源系统与电网交互功率约束：
[0111][0112]
其中为第i个国家综合能源系统与电网交互最大允许购电功率，为第i个国家综合能源系统与电网交互最大允许售电功率，为第i 个国家综合能源系统t时刻与电网交互的购电功率，为第i个国家综合能源系统t时刻与电网交互的售电功率；
[0113]
各国综合能源系统与热网交互功率约束：
[0114]
[0115]
其中为第i个国家综合能源系统与热网交互最大允许购电功率，为第i个国家综合能源系统与热网交互最大允许售电功率，为第i 个国家综合能源系统t时刻与热网交互的购热功率，为第i个国家综合能源系统t时刻与热网交互的售热功率；
[0116]
各国综合能源系统间电、热网交互功率约束，具体如下：
[0117]
各国综合能源系统间电网交互功率约束：
[0118][0119]
其中为i、j两国综合能源系统间传输电功率最大值，为i、j两国综合能源系统间传输热功率最大值，为i、j两国综合能源系统间传输的电功率，为i、j两国综合能源系统间传输的热功率；
[0120]
各国综合能源系统间热网交互功率约束：
[0121][0122]
其中为i、j两国综合能源系统间电网传输的电功率损耗最大值，为i、j两国综合能源系统间热网传输的热功率损耗最大值，为i、j 两国综合能源系统间电网传输的电功率损耗，为i、j两国综合能源系统间热网传输的热功率损耗；
[0123]
所述求解器为cplex求解器或gurobi求解器。
[0124]
本发明一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，首先构建跨境综合能源系统cbies，所述跨境综合能源系统cbies包括燃气轮机、燃气锅炉、可再生能源机组、储能设备、电锅炉，其中可再生能源机组包括风力发电和光伏发电，储能设备包括蓄电池和蓄热槽；其次提出考虑热网、电网损耗的跨境综合能源系统模型，在充分考虑跨境能量流动与转化的基础上，引入跨境绿证、碳交易，建立考虑碳-绿证联合交易的跨境综合能源系统协同优化模型；最后，通过所建立模型仿真求得协同运行优化结果。本发明提供的优化方法能够减少不同国家的综合能源系统运行成本，也提高了能源优化配置与可再生能源消纳量。
具体实施方式
[0125]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施案例对本发明进行深入地详细说明。应当理解，此处所描述的具体实施案例仅仅用以解释本发明，并不用于限定发明。
[0126]
本发明提供的一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，包括如下步骤：
[0127]
s1、构建跨境综合能源系统cbies模型，所述跨境综合能源系统cbies包括燃气轮机、燃气锅炉、可再生能源机组、储能设备和电锅炉，其中可再生能源机组包括风力发电和光伏发电，储能设备包括蓄电池和蓄热槽；
[0128]
s2、在步骤s1所述的跨境综合能源系统cbies模型基础上引入跨境绿证交易模型和碳交易模型，建立考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型；
[0129]
s3、根据步骤s2中建立的考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型，使用求解器求解得出优化结果。
[0130]
s11、所述步骤s1中构建的跨境综合能源系统cbies模型具体包括燃气轮机模型、燃气锅炉模型、可再生能源机组模型、储能设备模型和电锅炉模型，具体如下：
[0131]
燃气轮机模型具体如下：
[0132][0133][0134]
其中表示第i个国家综合能源系统t时刻燃气轮机的输出电功率，表示第i个国家综合能源系统t时刻燃气轮机的输出热功率，l
ng
表示天然气的低位热值，表示燃气轮机在t时刻的天然气消耗量，η
gt
表示燃气轮机的发电效率，n
gt
表示余热回收系数；
[0135]
燃气锅炉模型具体如下：
[0136][0137]
其中表示第i个国家综合能源系统t时刻燃气锅炉的输出热功率，l
ng
表示天然气的低位热值，表示燃气锅炉在t时刻的天然气消耗量，η
gb
表示燃气锅炉的发热效率；
[0138]
可再生能源机组模型包括风力发电模型和光伏发电模型，具体如下：
[0139]
风力发电模型：
[0140][0141]
其中为第i个国家综合能源系统t时刻风电功率，ρ为空气密度，a为风轮机叶片扫掠面积，v为风速，c
wt
为风能利用系数，是单位时间内风轮所吸收的风能与通过风轮旋转面的全部风能之比，λ
wt
为叶尖速率比；
[0142]
光伏发电模型：
[0143][0144]
ts＝ta+0.0138
·
(1+0.031ta)
·
(1-0.042v)
·gꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(47)
[0145]
其中为第i个国家综合能源系统t时刻光伏发电功率，g为光照强度，ts为光伏电池表面温度，p
stc
、g
stc
、t
stc
分别为标准测试条件下的最大输出功率、光照强度和光伏电池表面温度，ε为光伏电池温度系数；ta为环境温度；v为风速；
[0146]
储能设备模型包括蓄电池模型和蓄热槽模型，具体如下：
[0147]
蓄电池模型：
[0148]
放电
[0149][0150]
充电
[0151][0152]
其中为第i个国家综合能源系统t+1时刻蓄电池中存储的电量，为第 i个国家综合能源系统t时刻蓄电池中存储的电量，和分别为第i个国家综合能源系统t时刻蓄电池放电功率和充电功率；和分别为蓄电池自身的放电效率和充电效率；和分别为蓄电池自身的放电损耗和充电损耗；
[0153]
蓄热槽模型：
[0154]
放热
[0155][0156]
充热
[0157][0158]
其中为第i个国家综合能源系统t+1时刻蓄热槽中存储的热能，w
th,i
为第i个国家综合能源系统t时刻蓄热槽中存储的热能，和分别为第i 个国家综合能源系统t时刻蓄热槽放热功率和充热功率；和分别为蓄热槽自身的放热效率和充热效率；和分别为蓄热槽自身的放热损耗和充热损耗；
[0159]
电锅炉模型，具体如下：
[0160][0161]
其中，为第i个国家综合能源系统t时刻电锅炉的制热功率；为第 i个国家综合能源系统t时刻电锅炉所需电功率，η
eb
为电锅炉的转换效率。
[0162]
s21、在步骤s1所述的跨境综合能源系统cbies模型基础上引入碳交易模型具体包括碳交易成本，具体如下：
[0163][0164]
其中为跨境综合能源系统的碳交易成本，n为综合能源系统个数，为第i个国家综合能源系统碳排放量，为第i个国家综合能源系统碳初始配额，为市场上的碳
排放权价格；d为碳排放量区间长度；σ为各阶梯碳排放权价格增长幅度，每上升一个阶梯，碳排放权价格增长碳排放权价格增长时，碳交易成本为负值，此时表示出售碳排放权获得收益；
[0165]
第i个国家综合能源系统碳排放量计算如下：
[0166][0167]
其中为跨境综合能源系统直接碳排放量，为外购电力的碳排放量，为外购热力的碳排放量；
[0168]
跨境综合能源系统直接碳排放量计算如下：
[0169][0170]
其中n为综合能源系统个数；t为小时数，取为24小时；是第i个国家综合能源系统中燃气轮机t时刻输出的电功率，是第i个国家综合能源系统中燃气锅炉t时刻输出的热功率；和分别为第i个国家综合能源系统中的燃气轮机发电效率和燃气锅炉发热效率；
[0171]
外购电力的碳排放量计算如下：
[0172][0173]
其中表示第i个国家综合能源系统外购电力消耗量，单位为万千瓦时；为电网基准线排放因子，表示消耗单位电力的间接排放量，本文取值 0.8367，单位为吨co2/千瓦时；
[0174]
外购热力的碳排放量计算如下：
[0175][0176]
其中表示第i个国家综合能源系统外购热力消耗量，单位为百万千焦；为热网基准线排放因子，表示消耗单位热力的间接排放量，本文取0.096，单位为吨co2/百万千焦；
[0177]
第i个国家综合能源系统碳初始配额计算如下：
[0178][0179]
其中为第i个国家综合能源系统燃气机组t时刻供电量，单位为mwh， p
e,b
为燃气机组的供电基准值，本文取值0.3791，单位为tco2/mwh，γh为燃气机组供热量修正系数，燃气机组供热量修正系数为1；为第i个国家综合能源系统燃气机组供热量，单位为gj，p
h,b
为燃气机组供热基准值，本文取值 0.0600，单位为tco2/gj；
[0180]
s22、在步骤s1所述的跨境综合能源系统cbies模型基础上引入跨境绿证交易模型具体包括跨境绿证交易成本，具体如下：
[0181][0182]
其中为第i个国家综合能源系统跨境绿证交易成本，n为综合能源系统个数，为第i个国家综合能源系统可再生能源实际消纳量，为第i个国家综合能源系统可再生能源配额量，和分别为购买和出售绿色证书的价格，λf为惩罚系数；
[0183]
第i个国家综合能源系统可再生能源配额量计算如下：
[0184][0185]
其中为第i个国家综合能源系统区域预测年消纳可再生能源量，为第i个国家综合能源系统区域预测年接入可再生能源量，为第i个国家综合能源系统区域预计全社会年用电量，n为综合能源系统个数；
[0186]
s23、所述步骤s2中建立考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型，具体如下：
[0187][0188]
其中β为跨境绿证交易与碳交易权折算系数，αz为交易价格调整系数，由可再生能源类型及当地燃气价格确定；为市场上的碳排放权价格，为出售绿色证书的价格；
[0189]
所述考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型包括系统运行总成本，具体如下：
[0190][0191]
其中ci、和分别为第i个国家综合能源系统运行总成本、燃料成本、购能成本、售能收益、系统维护成本，各国综合能源系统间交互成本，碳交易成本和绿证交易成本，单位元；
[0192]
第i个国家综合能源系统运行燃料成本：
[0193][0194]
其中c
fuel
为燃料费用，单位元，表示燃气轮机在t时刻的天然气消耗量，表示燃气锅炉在t时刻的天然气消耗量；
[0195]
第i个国家综合能源系统运行购能成本：
[0196][0197]
其中和分别为各国综合能源系统t时刻向电网购电费用和购热费用，单位元；和分别为第i个国家综合能源系统t时刻与电网交互的购电功率和与热网交互的购热功率，；
[0198]
第i个国家综合能源系统运行售能收益：
[0199][0200]
其中和分别为售电费用和售热费用，单位元；和分别为第i个国家综合能源系统t时刻与电网交互的售电功率和与热网交互的售热功率；
[0201]
第i个国家综合能源系统运行系统维护成本：
[0202][0203]
其中c
gt
、c
gb
、c
eb
、c
hs
、c
es
、c
pv
、c
wt
分别为燃气轮机、燃气锅炉、电锅炉、蓄热槽、蓄电池、光伏、风机的单位功率维护成本，单位元，为第i 个国家综合能源系统t时刻燃气轮机的输出电功率，为第i个国家综合能源系统t时刻燃气锅炉的输出热功率，为第i个国家综合能源系统t时刻电锅炉的制热功率，为第i个国家综合能源系统t时刻蓄热槽放热功率或充热功率，为第i个国家综合能源系统t时刻蓄电池充电功率或放电功率，为第 i个国家综合能源系统t时刻光伏发电功率，为第i个国家综合能源系统t 时刻风电功率；
[0204]
各国综合能源系统间交互成本：
[0205][0206]
其中何分别为综合能源系统间购电何购热费用，为i、j两国综合能源系统间传输的电功率，为i、j两国综合能源系统间传输的热功率；
[0207]
所述考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型的约束条件具体包括各国综合能源系统内部功率约束、各国综合能源系统与电网交互功率约束、各国综合能源系统与热网交互功率约束、各国综合能源系统间电网交互功率约束和各国综合能源系统间热网交互功率约束具体如下：
[0208]
各国综合能源系统内部功率约束具体包括电功率平衡约束、热功率平衡约束、各国综合能源系统设备出力上下限约束、蓄电池功率约束和蓄热槽功率约束，具体如下：
[0209]
电功率平衡约束：
[0210][0211]
其中，为第i个国家综合能源系统t时刻与电网交互的购电功率，为第i个国家综合能源系统t时刻与电网交互的售电功率，为第i个国家综合能源系统t时刻燃气轮机的输出电功率，为i、j两国综合能源系统间电网传输的电功率损耗，为i、j两国综合能源系统间传输的电功率，为第i个国家综合能源系统t时刻电锅炉所需电功率，分别为第i个国家综合能源系统t时刻蓄电池充电功率、放电功率，为第i个国家综合能源系统t时刻光伏发电功率，为第i个国家综合能源系统t时刻风电功
率，l
ie
为第i个国家综合能源系统的电负荷；
[0212]
热功率平衡约束：
[0213][0214]
其中为第i个国家综合能源系统t时刻与热网交互的购热功率，为第i个国家综合能源系统t时刻与热网交互的售热功率,为第i个国家综合能源系统t时刻燃气锅炉的输出热功率，为第i个国家综合能源系统t时刻电锅炉的制热功率，为i、j两国综合能源系统间热网传输的热功率损耗，为i、j两国综合能源系统间传输的热功率，分别为第i个国家综合能源系统t时刻蓄热槽放热功率、充热功率，为第i个国家的热负荷；
[0215]
各国综合能源系统设备出力上下限约束：
[0216][0217]
其中为第i个国家综合能源系统中设备m的电功率，为第i个国家综合能源系统中设备m的电功率上下限；为设备m的电功率，为设备m的电功率，为第i个国家综合能源系统中设备m的热功率上下限；
[0218]
蓄电池功率约束：
[0219][0220][0221][0222][0223]
为第i个国家综合能源系统蓄电池容量，为最大充电倍率，为最大放电倍率，分别为第i个国家综合能源系统t时刻蓄电池充电功率、放电功率；为t时刻充能的状态位，为t时刻放能的状态位，是0-1变量，表示同一设备同一时刻充放能状态唯一，为第i个国家综合能源系统蓄电池的最大、最小储能量；
[0224]
蓄热槽功率约束：
[0225][0226][0227][0228][0229]
其中为第i个国家综合能源系统蓄热槽容量，为最大充热倍率，为最大放热倍率，分别为第i个国家综合能源系统t时刻蓄热槽放热功率、充热功率，为第i个国家综合能源系统蓄热槽的最大、最小储热量；
[0230]
各国综合能源系统与电、热网交互功率约束，具体如下：
[0231]
各国综合能源系统与电网交互功率约束：
[0232][0233]
其中为第i个国家综合能源系统与电网交互最大允许购电功率，为第i个国家综合能源系统与电网交互最大允许售电功率，为第i 个国家综合能源系统t时刻与电网交互的购电功率，为第i个国家综合能源系统t时刻与电网交互的售电功率；
[0234]
各国综合能源系统与热网交互功率约束：
[0235][0236]
其中为第i个国家综合能源系统与热网交互最大允许购电功率，为第i个国家综合能源系统与热网交互最大允许售电功率，为第i 个国家综合能源系统t时刻与热网交互的购热功率，为第i个国家综合能源系统t时刻与热网交互的售热功率；
[0237]
各国综合能源系统间电、热网交互功率约束，具体如下：
[0238]
各国综合能源系统间电网交互功率约束：
[0239][0240]
其中为i、j两国综合能源系统间传输电功率最大值，为i、j两国综合能源系统间传输热功率最大值，为i、j两国综合能源系统间传输的电功率，为i、j两国综合能源系统间传输的热功率；
[0241]
各国综合能源系统间热网交互功率约束：
[0242][0243]
其中为i、j两国综合能源系统间电网传输的电功率损耗最大值，为i、j两国综合能源系统间热网传输的热功率损耗最大值，为i、j 两国综合能源系统间电网传输的电功率损耗，为i、j两国综合能源系统间热网传输的热功率损耗；
[0244]
所述求解器为cplex求解器或gurobi求解器。
[0245]
实施例
[0246]
跨境综合能源系统cbies通过电网、供热管网进行跨境交互，并向各自国家综合能源系统的负荷供应电能、热能；跨境综合能源系统配备了微燃机、燃气锅炉、燃气轮机、电锅炉、储电、储热、光伏、风机和余热回收装置且不同国家分别设置参数，将构建好的模型在配备i7cpu和16gb ram的计算机上进行模拟，仿真平台为matlab2020b。
[0247]
(1)跨境综合能源系统的参数和仿真设置
[0248]
表1给出了仿真中使用的设备容量及参数，表2设置了跨境综合能源系统运行情景
的六组案例，前三组案例不考虑能源系统跨境交互，后三组案例考虑能源系统跨境交互，再考虑是否引入碳、绿证交易来对比跨境综合能源系统规划结果差异。
[0249]
表1设备容量及参数
[0250][0251]
表2跨境综合能源系统运行情景
[0252][0253]
case1、case2、case3三个国家综合能源系统独立运行，不通过电、热网在系统间传输能量；case4、case5、case6三个国家综合能源系统协同运行，通过电、热网在系统间传输能量，且当系统间传输能量时，购买能量的价格小于从电网购电的价格，两种方式都以三个综合能源系统总运行成本最小为目标函数。
[0254]
(2)优化分析
[0255]
六组案例的优化结果分别见表3和表4。
[0256]
表3不考虑跨境交互时各能源系统运行成本
[0257][0258]
由表3可知，与case1相比，考虑绿色证书交易的case2，各国综合能源系统的运行成本分别降低了57680、50260、41660元；与case2相比，考虑绿证及碳交易的case3，各国综
合能源系统的运行成本分别降低了46440、43970、43160 元。
[0259]
表4考虑跨境交互时各能源系统运行成本
[0260][0261]
由表4可知，与case4相比，考虑绿色证书交易的case5，各国综合能源系统的运行成本分别降低了57640、49560、52890元；与case5相比，考虑绿证及碳交易的case6，各国综合能源系统的运行成本分别降低了43230、42360、45623 元。
[0262]
从case1、4，case2、5，case3、6的对比中看出，在考虑系统跨境交互以后，a国综合能源系统运行成本分别下降了50580、65830、66410元，b国综合能源系统运行成本分别下降了50540、66120、77640元，c国综合能源系统运行成本分别下降了47330、64510、80103元。
[0263]
考虑碳交易与绿证对跨境对各国年运行成本及风光消纳量的影响，可以得到以下重要结论：
[0264]
a)在跨境地区建立碳、绿证交易机制的跨境综合能源系统可有效降低系统运行成本，引导各能源系统降低自身碳排放，使得跨境综合能源系统为实现收益最大优先考虑消纳可再生能源；
[0265]
b)考虑能源系统间的跨境交互可以调控国家间资源配置不平衡的问题，可有效减少电网购电量，提高使用天然气的比例与可再生能源消纳量，对于减少系统机组碳排放量具有重要意义，是实现跨境综合能源系统协同运行的重要方式。
[0266]
综上所述，仿真结果验证了所提出的在跨境综合能源系统cbies中进行考虑碳-绿证优化规划模型的有效性，结果表明，碳-绿证优化规划模型在跨境综合能源系统cbies中的设置只带来了很少的投资成本，但给系统运行带来了可观的收益，因此，本发明提供的一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法具有潜在的工程应用价值。
[0267]
以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
技术特征：1.一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，其特征在于，包括如下步骤：s1、构建跨境综合能源系统cbies模型，所述跨境综合能源系统cbies包括燃气轮机、燃气锅炉、可再生能源机组、储能设备和电锅炉，其中可再生能源机组包括风力发电和光伏发电，储能设备包括蓄电池和蓄热槽；s2、在步骤s1所述的跨境综合能源系统cbies模型基础上引入跨境绿证交易模型和碳交易模型，建立考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型；s3、根据步骤s2中建立的考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型，使用求解器求解得出优化结果。2.根据权利要求1所述的一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，其特征在于：s11、所述步骤s1中构建的跨境综合能源系统cbies模型具体包括燃气轮机模型、燃气锅炉模型、可再生能源机组模型、储能设备模型和电锅炉模型，具体如下：燃气轮机模型具体如下：燃气轮机模型具体如下：其中表示第i个国家综合能源系统t时刻燃气轮机的输出电功率，表示第i个国家综合能源系统t时刻燃气轮机的输出热功率，l
ng
表示天然气的低位热值，表示燃气轮机在t时刻的天然气消耗量，η
gt
表示燃气轮机的发电效率，n
gt
表示余热回收系数；燃气锅炉模型具体如下：其中表示第i个国家综合能源系统t时刻燃气锅炉的输出热功率，l
ng
表示天然气的低位热值，表示燃气锅炉在t时刻的天然气消耗量，η
gb
表示燃气锅炉的发热效率；可再生能源机组模型包括风力发电模型和光伏发电模型，具体如下：风力发电模型：其中为第i个国家综合能源系统t时刻风电功率，ρ为空气密度，a为风轮机叶片扫掠面积，v为风速，c
wt
为风能利用系数，是单位时间内风轮所吸收的风能与通过风轮旋转面的全部风能之比，λ
wt
为叶尖速率比；光伏发电模型：t
s
＝t
a
+0.0138
·
(1+0.031t
a
)
·
(1-0.042v)
·
g
ꢀꢀꢀꢀꢀ
(6)其中为第i个国家综合能源系统t时刻光伏发电功率，g为光照强度，t
s
为光伏电池表面温度，p
stc
、g
stc
、t
stc
分别为标准测试条件下的最大输出功率、光照强度和光伏电池表面温
度，ε为光伏电池温度系数；t
a
为环境温度；v为风速；储能设备模型包括蓄电池模型和蓄热槽模型，具体如下：蓄电池模型：放电充电其中为第i个国家综合能源系统t+1时刻蓄电池中存储的电量，w
te,i
为第i个国家综合能源系统t时刻蓄电池中存储的电量，和分别为第i个国家综合能源系统t时刻蓄电池放电功率和充电功率；和分别为蓄电池自身的放电效率和充电效率；和分别为蓄电池自身的放电损耗和充电损耗；蓄热槽模型：放热充热其中为第i个国家综合能源系统t+1时刻蓄热槽中存储的热能，w
th,i
为第i个国家综合能源系统t时刻蓄热槽中存储的热能，和分别为第i个国家综合能源系统t时刻蓄热槽放热功率和充热功率；和分别为蓄热槽自身的放热效率和充热效率；和分别为蓄热槽自身的放热损耗和充热损耗；电锅炉模型，具体如下：其中，为第i个国家综合能源系统t时刻电锅炉的制热功率；为第i个国家综合能源系统t时刻电锅炉所需电功率，η
eb
为电锅炉的转换效率。3.根据权利要求2所述的一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，其特征在于：s21、在步骤s1所述的跨境综合能源系统cbies模型基础上引入碳交易模型具体包括碳交易成本，具体如下：
其中为跨境综合能源系统的碳交易成本，n为综合能源系统个数，为第i个国家综合能源系统碳排放量，为第i个国家综合能源系统碳初始配额，为市场上的碳排放权价格；d为碳排放量区间长度；σ为各阶梯碳排放权价格增长幅度，每上升一个阶梯，碳排放权价格增长放权价格增长时，碳交易成本为负值，此时表示出售碳排放权获得收益；第i个国家综合能源系统碳排放量计算如下：其中为跨境综合能源系统直接碳排放量，为外购电力的碳排放量，为外购热力的碳排放量；跨境综合能源系统直接碳排放量计算如下：其中n为综合能源系统个数；t为小时数，取为24小时；是第i个国家综合能源系统中燃气轮机t时刻输出的电功率，是第i个国家综合能源系统中燃气锅炉t时刻输出的热功率；和分别为第i个国家综合能源系统中的燃气轮机发电效率和燃气锅炉发热效率；外购电力的碳排放量计算如下：其中表示第i个国家综合能源系统外购电力消耗量，单位为万千瓦时；为电网基准线排放因子，表示消耗单位电力的间接排放量，本文取值0.8367，单位为吨co2/千瓦时；外购热力的碳排放量计算如下：其中表示第i个国家综合能源系统外购热力消耗量，单位为百万千焦；为热网基准线排放因子，表示消耗单位热力的间接排放量，本文取0.096，单位为吨co2/百万千焦；
第i个国家综合能源系统碳初始配额计算如下：其中为第i个国家综合能源系统燃气机组t时刻供电量，单位为mwh，p
e,b
为燃气机组的供电基准值，本文取值0.3791，单位为tco2/mwh，γ
h
为燃气机组供热量修正系数，燃气机组供热量修正系数为1；为第i个国家综合能源系统燃气机组供热量，单位为gj，p
h,b
为燃气机组供热基准值，本文取值0.0600，单位为tco2/gj；s22、在步骤s1所述的跨境综合能源系统cbies模型基础上引入跨境绿证交易模型具体包括跨境绿证交易成本，具体如下：其中为第i个国家综合能源系统跨境绿证交易成本，n为综合能源系统个数，为第i个国家综合能源系统可再生能源实际消纳量，为第i个国家综合能源系统可再生能源配额量，和分别为购买和出售绿色证书的价格，λ
f
为惩罚系数；第i个国家综合能源系统可再生能源配额量计算如下：其中为第i个国家综合能源系统区域预测年消纳可再生能源量，为第i个国家综合能源系统区域预测年接入可再生能源量，为第i个国家综合能源系统区域预计全社会年用电量，n为综合能源系统个数；s23、所述步骤s2中建立考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型，具体如下：其中β为跨境绿证交易与碳交易权折算系数，α
z
为交易价格调整系数，由可再生能源类型及当地燃气价格确定；为市场上的碳排放权价格，为出售绿色证书的价格；所述考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型包括系统运行总成本，具体如下：其中c
i
、和分别为第i个国家综合能源系统运行总成本、燃料成本、购能成本、售能收益、系统维护成本，各国综合能源系统间交互成本，碳交易成本和绿证交易成本，单位元；第i个国家综合能源系统运行燃料成本：
其中c
fuel
为燃料费用，单位元，表示燃气轮机在t时刻的天然气消耗量，表示燃气锅炉在t时刻的天然气消耗量；第i个国家综合能源系统运行购能成本：其中和分别为各国综合能源系统t时刻向电网购电费用和购热费用，单位元；和分别为第i个国家综合能源系统t时刻与电网交互的购电功率和与热网交互的购热功率，；第i个国家综合能源系统运行售能收益：其中和分别为售电费用和售热费用，单位元；和分别为第i个国家综合能源系统t时刻与电网交互的售电功率和与热网交互的售热功率；第i个国家综合能源系统运行系统维护成本：其中c
gt
、c
gb
、c
eb
、c
hs
、c
es
、c
pv
、c
wt
分别为燃气轮机、燃气锅炉、电锅炉、蓄热槽、蓄电池、光伏、风机的单位功率维护成本，单位元，为第i个国家综合能源系统t时刻燃气轮机的输出电功率，为第i个国家综合能源系统t时刻燃气锅炉的输出热功率，为第i个国家综合能源系统t时刻电锅炉的制热功率，为第i个国家综合能源系统t时刻蓄热槽放热功率或充热功率，为第i个国家综合能源系统t时刻蓄电池充电功率或放电功率，为第i个国家综合能源系统t时刻光伏发电功率，为第i个国家综合能源系统t时刻风电功率；各国综合能源系统间交互成本：其中何分别为综合能源系统间购电何购热费用，为i、j两国综合能源系统间传输的电功率，为i、j两国综合能源系统间传输的热功率；所述考虑碳交易和跨境绿证交易的跨境综合能源系统协同运行优化模型的约束条件具体包括各国综合能源系统内部功率约束、各国综合能源系统与电网交互功率约束、各国综合能源系统与热网交互功率约束、各国综合能源系统间电网交互功率约束和各国综合能源系统间热网交互功率约束具体如下：各国综合能源系统内部功率约束具体包括电功率平衡约束、热功率平衡约束、各国综合能源系统设备出力上下限约束、蓄电池功率约束和蓄热槽功率约束，具体如下：
电功率平衡约束：其中，为第i个国家综合能源系统t时刻与电网交互的购电功率，为第i个国家综合能源系统t时刻与电网交互的售电功率，为第i个国家综合能源系统t时刻燃气轮机的输出电功率，为i、j两国综合能源系统间电网传输的电功率损耗，为i、j两国综合能源系统间传输的电功率，为第i个国家综合能源系统t时刻电锅炉所需电功率，分别为第i个国家综合能源系统t时刻蓄电池充电功率、放电功率，为第i个国家综合能源系统t时刻光伏发电功率，为第i个国家综合能源系统t时刻风电功率，为第i个国家综合能源系统的电负荷；热功率平衡约束：其中为第i个国家综合能源系统t时刻与热网交互的购热功率，为第i个国家综合能源系统t时刻与热网交互的售热功率,为第i个国家综合能源系统t时刻燃气锅炉的输出热功率，为第i个国家综合能源系统t时刻电锅炉的制热功率，为i、j两国综合能源系统间热网传输的热功率损耗，为i、j两国综合能源系统间传输的热功率，分别为第i个国家综合能源系统t时刻蓄热槽放热功率、充热功率，为第i个国家的热负荷；各国综合能源系统设备出力上下限约束：其中为第i个国家综合能源系统中设备m的电功率，为第i个国家综合能源系统中设备m的电功率上下限；为设备m的电功率，为设备m的电功率，为第i个国家综合能源系统中设备m的热功率上下限；蓄电池功率约束：蓄电池功率约束：蓄电池功率约束：蓄电池功率约束：蓄电池功率约束：为第i个国家综合能源系统蓄电池容量，为最大充电倍率，为最大放电倍率，分别为第i个国家综合能源系统t时刻蓄电池充电功率、放电功率；为t时刻充能的状态位，为t时刻放能的状态位，是0-1变量，表示同一设备同一时刻充放能状态唯一，为第i个国家综合能源系统蓄电池的最大、最小储能量；
蓄热槽功率约束：蓄热槽功率约束：蓄热槽功率约束：蓄热槽功率约束：其中为第i个国家综合能源系统蓄热槽容量，为最大充热倍率，为最大放热倍率，分别为第i个国家综合能源系统t时刻蓄热槽放热功率、充热功率，为第i个国家综合能源系统蓄热槽的最大、最小储热量；各国综合能源系统与电、热网交互功率约束，具体如下：各国综合能源系统与电网交互功率约束：其中为第i个国家综合能源系统与电网交互最大允许购电功率，为第i个国家综合能源系统与电网交互最大允许售电功率，为第i个国家综合能源系统t时刻与电网交互的购电功率，为第i个国家综合能源系统t时刻与电网交互的售电功率；各国综合能源系统与热网交互功率约束：其中为第i个国家综合能源系统与热网交互最大允许购电功率，为第i个国家综合能源系统与热网交互最大允许售电功率，为第i个国家综合能源系统t时刻与热网交互的购热功率，为第i个国家综合能源系统t时刻与热网交互的售热功率；各国综合能源系统间电、热网交互功率约束，具体如下：各国综合能源系统间电网交互功率约束：其中为i、j两国综合能源系统间传输电功率最大值，为i、j两国综合能源系统间传输热功率最大值，为i、j两国综合能源系统间传输的电功率，为i、j两国综合能源系统间传输的热功率；各国综合能源系统间热网交互功率约束：其中为i、j两国综合能源系统间电网传输的电功率损耗最大值，为i、j两国综合能源系统间热网传输的热功率损耗最大值，为i、j两国综合能源系统间电网传输的
电功率损耗，为i、j两国综合能源系统间热网传输的热功率损耗；4.根据权利要求3所述一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，其特征在于：所述求解器为cplex求解器或gurobi求解器。
技术总结本发明涉及一种考虑碳-绿证联合交易的跨境综合能源系统协同运行优化方法，首先构建跨境综合能源系统CBIES，跨境综合能源系统CBIES包括燃气轮机、燃气锅炉、可再生能源机组、储能设备、电锅炉，其中可再生能源机组包括风力发电和光伏发电，储能设备包括蓄电池和蓄热槽；其次提出考虑热网、电网损耗的跨境综合能源系统模型，在充分考虑跨境能量流动与转化的基础上，引入跨境绿证、碳交易，建立考虑碳-绿证联合交易的跨境综合能源系统协同优化模型；最后，通过所建立模型仿真求得协同运行优化结果。本发明提供的优化方法能够减少不同国家的综合能源系统运行成本，也提高了能源优化配置与可再生能源消纳量。与可再生能源消纳量。
技术研发人员：骆钊 杨林燕 刘鸿志 黎博文 聂灵峰 田肖 董晨鸣 罗蒙顺 吴谕侯 朱家祥
受保护的技术使用者：昆明理工大学
技术研发日：2022.06.25
技术公布日：2022/11/1
