一种基于标准化矩阵模型的综合能源系统优化方法与流程

本发明属于综合能源系统设计规划领域，涉及一种综合能源优化方法，尤其是涉及一种面向结构、设计和运行的基于标准化矩阵模型的综合能源系统优化方法。

背景技术：

在ppp和新电改条件下，区域能源互联网的实质是多能互补基础上的综合能源服务。传统能源系统规划通常仅面对单一能源系统(如冷、热、电、天燃气)，不能优化配置各个能源系统的资源，整体能源资源利用效率低。针对这一问题，在解决分布式清洁能源的背景下，研究人员提出了综合能源系统协同规划的概念，即将冷、热、电、热水、天燃气等多种能源耦合，形成一个多种能源综合转换利用的物理系统，以充分发挥不同形式能源的互补与协同效应，在更大范围内实现各个能源系统的资源优化配置，提高系统运行灵活性的同时，提高分布式清洁能源消纳能力和系统综合能效与可靠性。国内外许多研究机构、能源企业已经开展了相关的项目探索，如国外的欧盟electra示范项目、欧盟e-dema项目、英国曼彻斯特示范工程、德国朗根费尔德示范工程、日本柏叶智慧城市、加拿大耶洛奈夫镇示范项目等与国内的延庆县示范工程、上海迪士尼度假区示范、崇明岛示范工程等。

能源综合服务系统的优化决策是一个复杂的系统问题，很大程度上取决于决策者偏好及其所确立的优化目标。从时间尺度上，能源综合服务系统的优化可以分为长期战略优化(战略规划)、中期战术优化(战术规划)及短期运行优化(运行规划)三个层次。任一层次均需实现该层次所规定的特定目标，同时该层次的决策将成为低层次的约束条件，如图1所示。

孔祥强等采用模式搜索和惩罚函数两种方法，以年运行费用最小为目标函数，对分布式热电联产系统的运行策略进行了优化分析。bischi等基于混合整数线性规划理论，以日运行费用最小为目标函数，考虑能源价格、运维费用、启停损耗、外界温度等相关因素，构建了分布式热电联产系统的运行规划模型。facci等分析了分布式热电联产系统常规“以热定电”和“以电定热”的模式，提出了系统运行优化的动态规划模型。黄子硕等分析了多能互补分布式能源系统综合能效的影响因素，并以夏季工况为例给出了几种典型配置下的综合能效水平。任洪波等论述了分布式热电联产系统设计优化研究进展，提出了由结构优化、设计优化和运行优化所构成的系统优化的层次型框架体系。相比于分布式热电联产系统，由于综合能源系统固有的多元性与复杂性，其系统优化研究将是当前及未来能源系统工程研究的主要课题之一。此外，尽管大量文章讨论了分布式热电联产系统的设计优化与运行优化研究，但很少有论文研究其标准化建模方法。基于调度系数，chicco等提出了一种为小规模冷、热、电三联供系统自动生成耦合矩阵的建模技术。由于涉及调度因子，所建立的综合能源系统规划模型是非线性的，该方法对于大型的综合能源系统而言太过复杂。almassalkhi等提出了一种能量“输入-存储-转换-存储-输出”结构的综合能源系统线性建模方法。然而，该建模方法只能处理固定结构形式的综合能源系统。当前，研究一种适用于各种结构形式的，且灵活自动的综合能源系统计算机标准化建模方法和适用于面向结构、设计和运行的综合优化方法具有重要意义。

技术实现要素：

本发明旨在解决采用传统综合能源系统优化方法所建立的综合能源系统规划模型是非线性的，对于大型的综合能源系统而言太过复杂，计算耗时长等问题，提出了一种基于标准化矩阵模型的综合能源系统优化方法。

本发明的目的可以通过以下技术方案来实现：

一种基于标准化矩阵模型的综合能源系统优化方法，其特征在于，该方法包括以下步骤：

1)建立综合能源系统标准化矩阵模型；

2)根据优化类型，缩放所述综合能源系统标准化矩阵模型，建立与所述优化类型对应的优化规划模型，采用智能优化算法获得优化结果。

进一步地，建立所述综合能源系统标准化矩阵模型时，以能源转换设备之间的能量流作为状态变量。

进一步地，所述综合能源系统标准化矩阵模型为：

式中，ec为供给电制冷设备的电量，hh为供给吸收式制冷设备的热量，fab为供给锅炉设备的燃料量，fchp为供给chp设备的燃料量，eh为供给电制热设备的电量，sh为光热、深层地热等产生的热量，ei为坚强主网架供电量，se为区域光伏产生的电量，we为区域风电产生的电量，vc,in为储冷设备充入的冷量，vc,out为储冷设备放出的冷量，vh,in为储热设备充入的热量，vh,out为储热设备放出的热量，ve,in为储电设备充入的电量，ve,out为储电设备放出的电量，为通过能源输送网络输送至用户端的总冷量，为通过能源输送网络输送至用户端的总热量，为通过能源输送网络输送至用户端的总电量，δec为冷量变化，δeh为热量变化，δee为电量变化，ηq、ηw分别为chp设备产热、产电效率，ηab为锅炉产热效率，ηc为电制冷设备的性能系数，ηh为电制热设备的性能系数，ηr为吸收式制冷设备的性能系数，ηch为充能效率，ηds为放能效率，下标e、c、h分别表示电、冷、热。

进一步地，所述优化类型包括综合能源系统结构优化、综合能源系统设计优化和综合能源系统运行优化。

进一步地，所述综合能源系统结构优化中，以最大化系统年综合能效为优化目标，优化确定系统组成单元及其相互联系，所述系统年综合能效综合考虑经济性指标、节能性指标和环境性指标。

进一步地，所述最大化系统年综合能效的优化函数表示为：

maxη＝f(α,β,x,y,z,u)

式中，变量x,α,β分别为燃气/燃油锅炉与chp设备消耗燃料量之比、通过电制冷设备提供的冷量与总冷量需求之比和通过电制热设备提供的热量与总热量需求之比，表示系统配置方案；变量y,z分别为总冷量需求与总电力需求之比和总热量需求与总电力需求之比，表示区域用户的总冷、热、电需求结构；u当地生物质、光伏、风电能源折合当量标准煤后与总电力需求之比，表示当地可再生能源利用情况。

进一步地，所述综合能源系统设计优化中，将综合能源系统设计成区域型分布式能源互联网络形式，优化获得系统组成单元的技术特性，所述技术特性包括机型、容量和台数。

进一步地，所述区域型分布式能源互联网络形式中，对于同类型能源转换设备，选择2台或2台以上设备并联运行。

进一步地，所述综合能源系统运行优化中，以区域最小化运维成本为优化目标，优化确定综合能源系统各设备单元最优的逐时运行策略。

进一步地，所述区域最小化运维成本的优化函数表示为：

式中，pi,t为第i种输入能源第t时段的价格，m为区域综合能源系统外购能源的种类，fi,t为第i种能源第t时段的外购量，t为时间周期。

与现有技术相比，本发明具有以如下有益效果：

1)本发明的综合能源系统标准化矩阵建模方法和优化方法简单高效，能根据不同综合能源系统应用场景，通过缩放综合能源系统标准化矩阵，快速建立相应的优化规划模型，在确保计算准确性的同时克服了以往非线性规划模型计算耗时长的缺点。

2)本发明所提出的综合能源系统标准化矩阵建模方法及优化方法能有效适用于区域综合能源系统包括结构、设计和运行的综合优化规划。

附图说明

图1为能源综合服务系统的优化决策过程示意图；

图2为一综合能源系统结构示意图；

图3为综合能源系统的标准化矩阵建模方法原理图；

图4为本发明建立的综合能源系统线性规划模型的物理结构示意图；

图5为综合能源系统结构优化示意图

图6为区域型分布式能源互联网络形式

图7为综合能源系统运行优化示意图

图8为系统最佳综合能效随区域总冷量需求与总电力需求的比值y的变化情况

图9为系统最佳综合能效随区域总热量需求与总电力需求的比值z的变化情况

图10为系统最优综合能效随区域总热/冷量需求与总电力需求的比值的变化情况

图11为某一区域典型冷、热、电三联供系统

图12为区域终端用户冷、热、电等逐时需求

图13为外购电的逐时电价

图14为区域系统各设备单元的最优逐时运行策略。

具体实施方式

下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。

本发明提供一种基于标准化矩阵模型的综合能源系统优化方法，包括以下步骤：1)建立综合能源系统标准化矩阵模型；2)根据优化类型，缩放所述综合能源系统标准化矩阵模型，建立与所述优化类型对应的优化规划模型，采用智能优化算法获得优化结果。

1、综合能源系统标准化矩阵建模

针对综合能源系统优化研究为目标，首先提出了一种综合能源系统标准化矩阵建模方法，从系统工程的视角出发，以能源转换设备之间的能量流作为状态变量来处理引入调度因子所导致的非线性问题，构建了综合能源系统线性规划模型。基于所提出的综合能源系统标准化矩阵模型，可用于开展综合能源系统包括结构、设计和运行的综合优化。

综合能源系统是一种集成多种能源输入、多种能源输出以及多种能源转换设备的系统，如图2所示。

综合能源系统是以坚强智能电网为基础的区域能源互联网，如图2所示。以系统工程的视角看，综合能源系统可由以下四个子系统组成：

①外部能源供应子系统：外部能源供应子系统支撑着整个综合能源系统的运行，供应整个系统运行需求的一次能源(天然气、燃油、煤)和二次能源(坚强主网架供电)，为整个系统运行提供能源保障。

②能源转换子系统：第一类是包括光伏发电(pv)、小型风力和小水力发电等的小规模可再生能源发电系统；第二类是以内燃机、燃气轮机、微燃机、燃料电池、斯特林发动机等为原动机的热电联产或冷、热、电三联产系统。第三类是辅助型能源转换系统，如燃气/油锅炉、储能设备等。通过各种方式将一次能源和二次能源高效的转化成终端用户需要的多种能源形式。

③能源输送网络：能源输送网络将各个能源转换子系统产生的能源高效合理的传输到每个需要不同形式能源的终端用户，包括电网、热网、冷网等。

④用户终端子系统：用户终端子系统最终将能源转换子系统产生的能源消耗的系统。在需求侧，将多个不同类型用户(如居民用户、工厂、商场、写字楼、医院等)融入同一供能体系，通过用户间互补、互动实现负荷平均化。

电力系统计算机标准化建模是电力系统自动化的重要基础，同样地，综合能源系统计算机标准化建模是未来多能源系统自动化的重要基础。

如图3所示，综合能源系统是一个多种形式能源输入vin，多种能源输出vout的双端口网络，本发明的中间的能量转换环节由耦合矩阵c表示：

vout＝cvin(1)

即

其中耦合矩阵c用以描述不同类型能源转换设备本身的转换特性。

各种能源转换设备的转换特性如下：

以内燃机、燃气轮机、微型燃机、燃料电池、斯特林发动机等为原动机的热电联产系统(chp)：

式中ηq、ηw分别为chp设备产热、产电效率。

燃气/油锅炉：

[vout]＝[ηab][vin](4)

式中ηab为锅炉产热效率，公式中“[]”表示向量。

电制冷设备：

[vout]＝[ηc][vin](5)

式中ηc为电制冷设备的性能系数。

电制热设备：

[vout]＝[ηh][vin](6)

式中ηh为电制热设备的性能系数。

吸收式制冷设备：

[vout]＝[ηr][vin](7)

式中ηr为吸收式制冷设备的性能系数。

储能设备：

储能(energystorage)除了充放能行为外，还包括内部储能状态(stateofcharge,soc)的变化：

其中ηch为充能效率，ηds为放能效率，δe为内部储能状态的变化。

储能设备第t时刻的储能状态为：

e(t)＝e(t-1)+δe(t)(9)

综合能源系统能源输送网络包括电网、热网、冷网、天然气网等，如图2所示。本发明提出的标准化矩阵建模方法通过以能源转换设备之间的能量流作为状态变量来处理引入调度因子所导致的非线性问题。由于任意两个能量流之间的关系是线性的，因此不需要引入调度因子作为决策变量，本发明的综合能源系统线性规划模型的物理结构如图4所示。

本发明的综合能源系统标准化矩阵模型可表示为：

式中ec为供给电制冷设备的电量；hh为供给吸收式制冷设备的热量；fab为供给锅炉设备的燃料量；fchp为供给热电联产(chp)系统的燃料量；eh为供给电制热设备的电量；sh为光热、深层地热等产生的热量；ei为坚强主网架供电量；se为区域光伏产生的电量；we为区域风电产生的电量；vc,in为储冷设备充入的冷量；vc,out为储冷设备放出的冷量；vh,in为储热设备充入的热量；vh,out为储热设备放出的热量；ve,in为储电设备充入的电量；ve,out为储电设备放出的电量；为通过能源输送网络输送至用户端的总冷量；为通过能源输送网络输送至用户端的总热量；为通过能源输送网络输送至用户端的总电量；δec为冷量变化；δeh为热量变化；δee为电量变化；ηq、ηw分别为chp设备产热、产电效率；ηab为锅炉产热效率；ηc为电制冷设备的性能系数；ηh为电制热设备的性能系数；ηr为吸收式制冷设备的性能系数；ηch为充能效率，ηds为放能效率，下标e、c、h分别表示电、冷、热。

2、综合能源系统结构优化

综合能源系统结构优化的目的是决定系统组成单元及其相互联系(能流结构和设备类型)。以系统年综合能效为优化目标，基于粒子群优化算法，构建了可用于开展区域综合能源系统最优结构规划的数学模型。

(1)优化目标

选取经济性、节能性和环保性3类指标作为优化目标。

1)经济性

选用系统年供热、冷、电总量作为经济性评价指标：

2)节能性

选取系统年一次能源消耗总量作为节能性评价指标：

式中θf为输入的化石燃料的折合标准煤系数，假设输入燃料为天燃气，取1；θe为国家电网所供电力的折合标准煤系数，国家电网火电占比75％时，取1.95。

3)环境性

综合能源系统的年co2排放量包括从天然气中co2的排放量及系统电力的co2排放量，计算式为：

式中为天然气的碳排放因子。

由上述可知，节能性指标和环境性指标只相差一个天然气的碳排放因子的乘积系数，成正比例关系，所以可以将环境性指标隐含在节能性指标中。因此，综合考虑经济性、节能性和环境性指标，采用系统年综合能效作为优化目标：

通过能源输送网络供给用户的总冷量需求

式中ec为用于电制冷的总电量需求，mj；hh为用于吸收式制冷的总热量需求，mj。

通过能源输送网络供给用户的总热量需求

式中fab为用于锅炉的总燃料量，mj；fchp为用于热电联产(chp)的总燃料量，mj；eh用于电制热的总电量需求，mj；sh为本地光热、深层地热等太阳能制热量，mj。

通过能源输送网络供给用户的总电量需求

式中ei为国家电网所供电量，mj；se为本地太阳能发电量，mj；we为本地风能发电量，mj。

设通过电制冷设备提供的冷量与总冷量需求之比：

设通过电制热设备(包括电加热和各种形式的电驱动热泵)提供的热量与总热量需求之比：

设系统中燃气/燃油锅炉与chp设备消耗燃料量的比值：

设系统总冷量需求与总电力需求的比值：

设系统总热量需求与总电力需求的比值：

设系统当地生物质、光伏、风电能源折合当量标准煤后与系统总电力需求的比值：

此时，可以得到无量纲化的系统年综合能效表达式：

式中表示通过chp设备消耗的燃料量与系统总电力需求的比值。

最终，以系统年综合能效为优化目标的优化函数为：

η＝f(α,β,x,y,z,u)(25)

式中变量x、α、β表示系统配置方案；变量y、z表示区域用户的总冷、热、电需求结构；u表示当地可再生能源利用情况。

综合能源系统结构优化的目的是为了确定系统组成单元及其相互联系(能流结构和设备类型)，如图5所示。建立优化目标后可以采用粒子群优化算法、遗传算法等全局优化算法进行求解。

3、综合能源系统设计优化

综合能源系统设计优化的目的是决定系统组成单元的技术特性(机型、容量和台数)。为提高综合能源系统供能可靠性，可以设计成区域型分布式能源互联网络形式，如图6所示，并且在满足区域用户的总冷、热、电量需求基础上，对于同类型能源转换设备，可以选择2台或2台以上设备并联运行，以此进一步提高系统综合能效与供能可靠性，实现区域综合能源系统的最优设计规划。

4、综合能源系统运行优化

综合能源系统运行优化的目的是决定系统组成单元的逐时运行策略。以区域综合能源系统运行优化研究为目标，并以最小化运维成本为优化目标函数，构建了可用于开展区域综合能源系统的最优运行规划的数学模型，来为指导综合能源系统各设备单元最优的逐时运行策略提供理论依据，如图7所示。

优化目标：最小化区域综合能源系统的运维成本，如式(26)：

其中pi,t为第i种输入能源第t时段的价格；m为区域综合能源系统外购能源的种类；fi,t为第i种能源第t时段的外购量。

本实施例基于上述优化方法实现综合能源系统结构优化和综合能源系统运行优化。

其中综合能源系统结构优化的具体实施案例如下：

在实施系统优化设计之前，作如下假设：

①区域用户的总冷量需求与总电力需求的比值取为y＝0.3；

②区域用户的总热量需求与总电力需求的比值取为z＝0.4；

③由式(24)可知，当地可再生能源利用与综合能效的关系比较确定，且当地可再生能源利用情况通常受本地自然环境资源和能源政策条件限制，这里假设取u＝0.2；

④当前能源转换设备技术水平见表1；

表1几种典型chp系统发电与产热参数

其余常用能源转换设备的技术水平基准值为ηh＝3，ηr＝1，ηc＝3.5，ηab＝0.90。

根据上述假设条件，系统的优化目标函数可以简化为：η＝f(α,β,x)。选取遗传算法(ga)、粒子群优化算法(pso)作为优化算法，相关参数设置如表2所示，用于搜索最优的年综合能效η。表3为系统优化结果，其中综合能源服务系统包含以燃料电池、内燃机、燃气轮机和汽轮机为原动机的热电联产系统分别表示为系统1、系统2、系统3和系统4。

表2pso算法和ga算法的相关参数

表3系统优化结果

由表3可知，相比ga算法，pso算法可以在同样种群规模和进化代数下更有效地搜索到最优目标结果。当区域用户冷/热/电需求、可再生能源资源利用情况和能源转换设备技术水平确定(如上述假设条件)时，综合能源服务系统配置方案为使用燃料电池为原动机的热电联产系统和吸收式制冷设备来供能，不使用电制冷/电制热设备、燃气/燃油锅炉时系统年综合能效最优，综合能效达83.5％。

基于区域用户侧冷热电需求变化情况，结合系统综合能效，综合能源系统结构优化计算结果如图8、图9、图10所示。

图10所示，为当u＝0.15，z＝0.1，y从0到1时，系统最优综合能效随区域总冷量需求与总电力需求的比值y的变化情况；图11所示，为当u＝0.15，y＝0.1，z从0到1时，系统最优综合能效随区域总热量需求与总电力需求的比值z的变化情况；图12所示，为当u＝0.15，y与z同时从0到1时，系统最优综合能效随区域总热/冷量需求与总电力需求的比值的变化情况。由上述三图可知，随着区域用户总冷量和/或总热量需求与总电力需求的比值增大，系统最优综合能效趋势为逐渐增大。

综合能源系统运行优化的具体实施案例如下：

采用某一区域典型冷、热、电三联供系统的应用案例来阐述。该系统由一个背压运行的chp设备，一个辅助锅炉(ab)，一个吸收式制冷设备(cerg)，及一个电制冷设备(warg)构成。采用本发明所提出的建模过程，其系统标准化矩阵模型的物理结构描述如图11所示。

该系统各设备单元的转换特性如表4所示。

表4系统各设备单元的转换特性

此时，综合能源系统线性规划矩阵模型(式(10))可缩放为：

其中，该区域终端用户冷、热、电等逐时需求如图12所示。

另外，天然气的气价为40euro/mwh且恒定。外购电的逐时电价如图13所示。

则优化目标of为：

区域能源供需平衡约束条件如下：

各设备单元转换特性约束条件如下：

0≤ηabfab,t≤400(32)

0≤ηqfchp,t≤160(33)

0≤ηcec,t≤300(34)

0≤ηrhh,t≤300(35)

采用matlab自带的fmincon求解函数对该区域综合能源系统线性规划模型进行求解，得到如图14所示的系统各设备单元的最优逐时运行策略。

图14为通过本发明所提出的综合能源系统运行优化方法得到的区域系统各设备单元在24小时内的最优逐时运行策略。在0～6小时阶段，由于电价较低，区域终端用户的电需求全部由外购电供应。此后，由于电价高于天然气价，且为满足区域终端用户的冷、热、电需求，chp设备、辅助锅炉(ab)、吸收式制冷设备(cerg)、电制冷设备(warg)及外购电方式协同运行。在21～24小时阶段，区域终端用户只有电需求，且电价回落，该区域终端用户的电需求全部由外购电供应。最终得到的最小化区域系统的运维成本为377.7725euro，计算总耗时0.367053s(采用4.0ghz双处理器的笔记本电脑)。

以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。