CN115660187B - 一种低碳小镇地源热泵容量优化配置方法 - Google Patents
一种低碳小镇地源热泵容量优化配置方法 Download PDFInfo
- Publication number
- CN115660187B CN115660187B CN202211363685.7A CN202211363685A CN115660187B CN 115660187 B CN115660187 B CN 115660187B CN 202211363685 A CN202211363685 A CN 202211363685A CN 115660187 B CN115660187 B CN 115660187B
- Authority
- CN
- China
- Prior art keywords
- heat pump
- ground source
- source heat
- power
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 45
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 38
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 10
- 238000003064 k means clustering Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 239000003245 coal Substances 0.000 claims description 19
- 238000005338 heat storage Methods 0.000 claims description 18
- 238000012423 maintenance Methods 0.000 claims description 13
- 238000010248 power generation Methods 0.000 claims description 12
- 238000009434 installation Methods 0.000 claims description 11
- 238000005315 distribution function Methods 0.000 claims description 10
- 230000008878 coupling Effects 0.000 claims description 9
- 238000010168 coupling process Methods 0.000 claims description 9
- 238000005859 coupling reaction Methods 0.000 claims description 9
- 230000035772 mutation Effects 0.000 claims description 9
- 239000000446 fuel Substances 0.000 claims description 8
- 230000014509 gene expression Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 238000010438 heat treatment Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- RZGZMLICFFEUIQ-UHFFFAOYSA-N 5-[(1-phenylcyclohexyl)amino]pentanoic acid Chemical compound C=1C=CC=CC=1C1(NCCCCC(=O)O)CCCCC1 RZGZMLICFFEUIQ-UHFFFAOYSA-N 0.000 claims description 6
- 230000020169 heat generation Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 5
- 210000000349 chromosome Anatomy 0.000 claims description 3
- 238000009833 condensation Methods 0.000 claims description 3
- 230000005494 condensation Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 2
- 241000039077 Copula Species 0.000 abstract description 2
- 230000001737 promoting effect Effects 0.000 abstract 2
- 238000010276 construction Methods 0.000 abstract 1
- 230000010485 coping Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 abstract 1
- 241000711969 Chandipura virus Species 0.000 description 20
- 208000015951 Cytophagic histiocytic panniculitis Diseases 0.000 description 20
- 230000007547 defect Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000000540 analysis of variance Methods 0.000 description 2
- 230000003749 cleanliness Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
Landscapes
- Wind Motors (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种低碳小镇地源热泵容量优化配置方法，结合k均值聚类与Copula函数生成电源侧和负荷侧的联合场景，通过F统计确定最优典型场景，以综合成本为优化目标建立优化配置模型，为应对风电的不确定性而引入区间数学，并利用量子衍生算法对所建立的优化模型进行求解，获取地源热泵最优配置区间。本发明为促进低碳建设，推广地热能的利用，提出计及不确定性的低碳小镇地源热泵容量优化配置方法，通过北方某低碳小镇的仿真算例验证了所提方法可有效降低年综合成本，同时减少弃风量。
Description
技术领域
本发明涉及综合能源系统领域，具体涉及一种低碳小镇地源热泵容量优化配置方法。
背景技术
近年来，风电因其清洁、可再生、装机灵活等优势，成为低碳小镇能源利用的重中之重。然而，风电的波动性和间歇性使得可靠地获取风电变得困难，因此，需要对低碳小镇的多种能源及其转换设备进行优化配置。
传统含风电力系统的优化中，一般采取随机优化和鲁棒优化两类方法来应对风电的不确定性问题。其中，随机优化以获得风电精确分布为前提，而鲁棒优化则不需要准确的概率分布函数，但由于鲁棒优化存在求解结果较保守的缺点，对系统的经济性和环保性可能产生不利影响。
发明内容
为了解决现有发明的不足，本文优化低碳小镇的经济性同时促进弃风消纳，解决低碳小镇地源热泵容量优化配置问题，为此，本发明提供了一种低碳小镇地源热泵容量优化配置方法。
所采用的技术方案如下：
一种低碳小镇地源热泵容量优化配置方法，包括如下步骤：
整理低碳小镇以往风电出力及负荷变化数据；
采用k-means聚类与Frank Copula函数构建电源侧和负荷侧的联合场景；
通过伪F-统计确定最优典型场景；
以综合成本为优化目标，建立优化配置目标函数；
引入区间数学，并利用量子衍生算法对所建立的目标函数进行求解，确定地源热泵最优配置方案。
进一步地，所述的采用k-means聚类与Frank Copula函数构建电源侧和负荷侧的联合场景，其具体方法如下：
步骤1、通过整理分析过往风电出力以及负荷变化的数据，求出在规划年限下的风电以及负荷的边际概率密度分布；
步骤2、采用Frank Copula函数表示风电功率与负荷之间关系；
步骤3、利用相关矩阵图得到树型结构，构造风电-负荷联合分布函数；
步骤4、对所构造的风电-负荷联合分布函数进行分时段采样，得到风电-负荷联合采样结果；对所构造的风电-负荷联合分布函数进行逆变换，得到每个时段的风电-负荷输出，最终形成考虑风电与负荷相关性的随机场景；
步骤5、采用K-means聚类算法对步骤4得到的风电与负荷相关性的随机场景进行聚类，通过多维PFS判别函数确定最佳聚类个数，聚类中心为典型风电-负荷曲线组合，得到典型的风电电源侧和负荷侧的联合场景。
所述的通过伪F-统计确定最优典型场景，其具体方法是：
对于P(P≥1)维变量的样本，定义一个伪F统计比率，如式(1)所示。
式中：分别为矩阵/>的迹；
m为样本数；k为类数；
分别为P维变量样本的类间和类内散布矩阵，其表达式如式(2)所示：
式中：xj为第j个样本向量；
vi为第i类Ci的聚类中心；
μij的表达式如式(3)所示：
对于式(1)，表征了类内距离，/>表征了类间距离，聚类数目k的增加，类内的距离呈负相关，类间的距离呈正相关；(m-k)/(k-1)与k呈负相关；
结合上述式(1)(2)(3)，求得PFS为最大值时的最大聚类数K，即得到最优典型场景。
所述的建立优化配置目标函数是以综合能源系统中地源热泵的投资成本与系统年运行成本之和最低，其表达式为：
式中：ξPl为年综合成本，ξCap为年投资成本，ξOpe为系统运行成本,D为冬季供暖期天数，一般取150天；
其中的年投资成本ξCap为：
ξCap＝ξInc+ξFmc (5)
其中，ξInc为年初始投资成本，ξFmc为地源热泵固定维护成本；
CBP为地源热泵的安装容量；为地源热泵的单位容量成本；
RBP为资金收回系数；r为贴现率，取5％；l为灵活性资源生命周期，与具体灵活性资源相关；
系统运行成本ζOpe为：
其中，ξFuel为燃料成本，ξOm为运行维护成本，为弃风惩罚成本；
ωcoal为消耗煤的价格,T表示时间被分成的总数量，n表示CHP机组数量，为第i台CHP机组燃煤在t时刻的消耗量，其用发电功率的另一种形式表示：
其中，ai，bi，ci为CHP机组i的耗煤比例数，为t时刻CHP机组i的发出的电能功率；
根据泵式CHP机组的运行原理，若t时刻第i台CHP机组纯凝工况下的发电功率为随着抽气量的增加，CHP机组的发电量/>和热输出/>有如下关系：
其中，Cv为电-热耦合参数；
将式(9)代入式(8)，即可得到t时刻第i台CHP机组的煤耗量和发电功率与供热功率/>的关系：
其中，是第i台CHP机组的煤耗系数，可由ai、bi、ci和Cv计算得出；
运行维护成本ξOm为：
其中，分别为t时刻CHP机组、电锅炉和地源热泵的输出功率，/>分别为CHP机组、电锅炉和地源热泵单位输出功率的运行维护成本；
弃风惩罚成本为：
其中，是弃风惩罚成本系数，m代表风电机组数，/>是t时刻风机的弃风量，其等于风电预测值与风电实际发电量的实际值的差。
所述的优化配置目标函数中还包括优化配置约束条件，所述的优化配置约束条件中包括供能平衡约束、设备运行约束和热泵容量约束。
优选地，所述的供能平衡约束为：
其中，表示在t时刻第j台风机实际的发电量，/>表示在t时刻第i台CHP机组实际的发电量，/>和/>分别表示在t时刻第h台电锅炉和第k台和地源热泵的耗电量，PL,t为t时刻系统中的电负荷，/>表示在t时刻第i台CHP机组实际的产热量，/>和/>分别表示在t时刻第h台电锅炉和第k台和地源热泵的产热量，QL,t为t时段内系统运行时的热量值，/>为t时段内系统运行时的热负荷。
所述的设备运行约束包括CHP机组出力约束、地源热泵处理约束和储热罐热量约束：
所述的CHP机组出力约束为：
其中：PCHP,t和QCHP,t表示t时刻CHP机组的电功率和热功率；PCHP 和表示CHP机组的电能功率的下限和上限，QCHP 和/>表示CHP机组的热功率的下限和上限；
Cv，PCHP,D，Cm，PCHP,C和PCHP,A为CHP机组的热电耦合比例数，表示CHP机组的电出力和热出力的热电耦合联系；
所述的热泵处理约束为：
其中，QBP地源热泵热出力，地源热泵热出力最大值；
所述的储热罐热量约束为：
其中，分别为蓄热罐在运行调度周期的初始时刻t0、结束时刻tN存储的热量，/> γHS 分别为蓄热罐运行中状态的上、下限,CHS为蓄热罐的安装容量。
优选地，所述的地源热泵容量约束为：
其中， CBP 为地源热泵的安装容量上限和下限，根据实际情况选取。
进一步地，所述的引入区间数学，并利用量子衍生算法对所建立的目标函数进行求解的方法是：
步骤1、输入典型日的风电功率数据、电负荷数据、热负荷数据以及各灵活性资源的参数、待定容灵活性资源的容量范围进行遗传计算的基本数据，设定好合理的种群数量、迭代次数计算基本参数；
步骤2、将目标函数中待求解的变量进行染色体实数编码，设定迭代次数g＝1，随机生成的初始种群数量为N；
步骤3、对于种群中的每个个体，将该个体对应的配置方案传递至运行层，确定运行层运行模型，在满足运行层运行约束条件的情况下，得到该配置方案下的最优运行成本；
步骤4、将运行层的最优值结果返回至目标函数的投资部分进行迭代，计算出每一个体对应的适应度值；
步骤5、保留种群中得到的最优个体，同时记录最优个体所对应的适应度值；
步骤6、进行种群的选择、交叉以及变异等操作，进而生成新种群，经过迭代选择、交叉和变异后的种群个体可能不满足所设定的范围，故需要进行检验操作，若满足则进行后续步骤，不满足则需重新进行选择、交叉和变异；
步骤7、判断此时是否达到了最大迭代次数G，如果不满足就使g＝g+1；返回步骤3，如果满足则跳出循环，输出最优个体的记录，得到低碳小镇地源热泵的最优配置方案。
本发明技术方案具有如下优点：
本发明利用区间数学表征风电出力的不确定性，利用k均值聚类与Copula函数生成电源侧和负荷侧的联合场景，通过F-统计确定最优典型场景，以综合成本为优化目标建立优化配置目标函数，引入区间数学并利用量子衍生算法对所建立的目标函数进行求解以进一步提高优化精度，从而解决了低碳小镇地源热泵容量优化配置问题,且有效降低年运行成本，本发明计提方法与现有鲁棒优化方法相比，可更有效的消纳弃风。
附图说明
为了更清楚地说明本发明具体实施方式，下面将对具体实施方式中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明所提供的低碳小镇地源热泵容量优化配置方法流程图；
图2是本发明所提供的典型日风电区间出力图示；
其中：中间曲线是风出力曲线，上下两根曲线分别为不确定性区间的上限和下限；
图3是本发明所提供的典型日用电负荷曲线图示；
图4a-图4d是本发明所提供的地源热泵容量配置中典型日功率平衡图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明提供了低碳小镇地源热泵容量优化配置方法，整理低碳小镇以往风电出力及负荷变化数据；采用k-means聚类与Frank Copula函数构建电源侧和负荷侧的联合场景；通过伪F-统计确定最优典型场景；以综合成本为优化目标，建立优化配置目标函数；引入区间数学，并利用量子衍生算法对所建立的目标函数进行求解，确定地源热泵最优配置方案。
在低碳小镇中，不同时间段、不同空间点处的风速和负荷均表现为非高斯分布，因此，存在源荷双侧的不确定性。风电和负荷的场景选择将决定低碳小镇能源容量优化配置方案的合理性，为兼顾源-荷双侧的不确定性，本发明基于k-means聚类和F-统计的源荷典型场景选择方法。首先，对小镇的以往的负荷及风电出力数据进行处理，剔除噪声数据，之后利用Frank Copula函数生成风电-负荷联合概率密度函数，对抽样结果采用F-统计确定最优聚类个数，并利用聚类方法获取具有代表性的典型场景。
负荷时序变化作为随机变量一定程度上受到了风速以及温度等环境因素的影响，因此可以从风电出力以及负荷曲线中得出非线性相关性。
本发明利用Frank Copula函数来构建负荷之间的非线性相关模型，并利用K-means聚类得到典型的风电及负荷典型场景组合。联合场景生成详细步骤如下：
步骤1、通过整理分析过往风电出力以及负荷变化的数据，可求出在规划年限下的风电以及负荷的边际概率密度分布；
步骤2、由于负荷和风电间有负相关，本发明采用Frank Copula函数来表示风电功率与负荷之间关系；
步骤3、根据所用的上述函数，利用相关矩阵图来得到树型结构，构造风电-负荷联合分布函数；
步骤4、对所构造的风电-负荷联合分布函数进行分时段采样，得到风电-负荷联合采样结果；对所构造的风电-负荷联合分布函数进行逆变换，得到每个时段的风电-负荷输出，最终形成考虑风电与负荷相关性的随机场景；
步骤5、采用K-means聚类算法对步骤4得到的风电与负荷相关性的随机场景进行聚类，通过多维PFS判别函数确定最佳聚类个数，聚类中心为典型风电-负荷曲线组合，得到典型的风电电源侧和负荷侧的联合场景。
本发明通过伪F-统计确定最优典型场景，其采用的下述多维PFS判别函数。
伪F-统计(Pseudo F-statistics.PFS)是来自方差分析(Analysis ofVariance.ANOVA)领域的一个统计量。对于P(P≥1)维变量的样本，定义一个“伪F统计比率”，如式(1)所示。
式中：分别为矩阵/>的迹；m为样本数；k为类数；/>分别为P维变量样本的类间和类内散布矩阵，表达式如式(2)所示。
式中：xj为第j个样本向量；vi为第i类Ci的聚类中心；μij的表达式为
对于式(1)，表征了类内距离，/>表征了类间距离。聚类数目k的增加，类内的距离呈负相关，类间的距离呈正相关。此外，(m-k)/(k-1)与k呈负相关。故伴随k的增加，PFS的值呈先增大后减小的趋势，在某个k值时PFS可能达到最大值，这是最优的簇数。也就是说，找到最佳聚类数k等价于找到最大PFS值。
本发明中低碳小镇地源热泵容量优化配置方法如下：
风电因其清洁、可再生、装机灵活等优势，成为低碳小镇能源利用的重中之重。然而，风电的波动性和间歇性使得可靠地获取风电变得困难。在传统的风电系统优化中，一般使用两种方法来处理风电的不确定性：随机优化和鲁棒优化。其中随机优化依赖于风能的准确分布，而鲁棒优化又存在保守性的问题。因此，本发明将既不依赖不确定参数的精确分布，且可以突出不确定参数对系统的影响的区间优化引入低碳小镇的优化配置方案，以应对典型日优化运行计算中风电不确定性的影响。通过区间数学将风电出力区间化，分别以优化运行结果的上限和下限对地源热泵容量进行优化，二者结果的均值作为最终配置结果。具体描述如下。
本发明考虑地源热泵的投资成本和系统运行成本，不考虑CHP机组、风电场、储热罐的投资成本，以综合能源系统中地源热泵的投资成本与系统年运行成本之和最低为目标函数，如式(4)所示：
式中：ξPl为年综合成本，ξCap为年投资成本，ξOpe为系统运行成本,D为冬季供暖期天数，一般取150天。
其中，年投资成本ξCap由两部分组成，即年初始投资成本ξInc和地源热泵固定维护成本ξFmc。
ξCap＝ξInc+ξFmc (5)
年初始投资成本是地源热泵容量相关的单位投资成本，具体计算如下：
其中，CBP为地源热泵的安装容量，为地源热泵的单位容量成本；RBP为资金收回系数；r为贴现率，本发明取5％；l为灵活性资源生命周期，与具体灵活性资源相关。
系统运行成本ζOpe为：
其中，ξFuel为燃料成本，ξOm为运行维护成本，为弃风惩罚成本。
燃料成本ξFuel：
其中，ωcoal为消耗煤的价格,T表示时间被分成的总数量，n表示CHP的数量，为第i台CHP机组燃煤在t时刻的消耗量，其可以用发电功率的另一种形式表示：
其中，ai，bi，ci为CHP机组i的耗煤比例数；为t时刻CHP机组i的发出的电能功率。根据泵式CHP机组的运行原理，若t时刻第i台CHP机组纯凝工况下的发电功率为/>随着抽气量的增加，CHP机组的发电量/>和热输出/>有如下关系：
其中，Cv为电-热耦合参数。将式(9)代入式(8)，即可得到t时刻第i台CHP机组的煤耗量和发电功率与供热功率/>的关系
其中，是第i台CHP机组的煤耗系数，可由ai、bi、ci和Cv计算得出。
运行维护成本ξOm：
其中，分别为t时刻热电联产机组、电锅炉和地源热泵的输出功率，/>分别为热电联产机组、电锅炉和地源热泵单位输出功率的运行维护成本。
弃风惩罚成本
由于风电具有随机性和不可控的特点，为进一步增加系统对风电的消纳，在运行层目标函数中加入弃风惩罚成本，如式(14)
其中，是弃风惩罚成本系数，m代表风电机组数，/>是t时刻风机的弃风量，其等于风电预测值与风电的实际发电量的实际值的差。
所设置的优化配置约束条件如下：
(1)供能平衡约束：
其中，表示在t时刻第j台风机实际的发电量，/>表示在t时刻第i台CHP机组实际的发电量，/>和/>分别表示在t时刻第h台电锅炉和第k台和地源热泵的耗电量，PL,t为t时刻系统中的电负荷，/>表示在t时刻第i台CHP机组实际的产热量，/>和/>分别表示在t时刻第h台电锅炉和第k台和地源热泵的产热量，QL,t为t时段内系统运行时的热量值，/>为t时段内系统运行时的热负荷，/>为t时段第g台内蓄热罐的做功(吸热为负，放热为正)。
(2)设备运行约束：
CHP机组出力约束如下：
其中：PCHP,t和QCHP,t表示t时刻CHP的电功率和热功率；PCHP 和表示CHP电能功率的下限和上限，QCHP 和/>表示CHP热功率的下限和上限。Cv，PCHP,D，Cm，PCHP,C和PCHP,A为CHP的热电耦合比例数，表示CHP电出力和热出力的热电耦合联系。
热泵处理约束如下：
其中，QBP地源热泵热出力，地源热泵热出力最大值。
储热罐热量约束如下：
其中，和分别为蓄热罐在运行调度周期的初始时刻t0、结束时刻tN存储的热量，/> γHS 分别为蓄热罐运行中状态的上、下限。CHS分别为蓄热罐的安装容量。
(3)热泵容量约束
其中， CBP 为热泵的安装容量上限和下限，根据实际情况选取。
本发明优化配置方案求解步骤如下：
步骤1、输入典型日的风电功率数据、电负荷数据、热负荷数据以及各灵活性资源的参数、待定容灵活性资源的容量范围等进行遗传计算的基本数据，设定好合理的种群数量、迭代次数等计算基本参数；
步骤2、将目标函数中待求解的变量进行染色体实数编码，设定迭代次数g＝1，随机生成的初始种群数量为N；
步骤3、对于种群中的每个个体，将该个体对应的配置方案传递至运行层，确定运行层运行模型，在满足运行层运行约束条件的情况下，得到该配置方案下的最优运行成本；
步骤4、将运行层的最优值结果返回至目标函数的投资部分进行迭代，计算出每一个体对应的适应度值；
步骤5、保留种群中得到的最优个体，同时记录最优个体所对应的适应度值；
步骤6、进行种群的选择、交叉以及变异等操作，进而生成新种群，经过迭代选择、交叉和变异后的种群个体可能不满足所设定的范围，故需要进行检验操作，若满足则进行后续步骤，不满足则需重新进行选择、交叉和变异；
步骤7、判断此时是否达到了最大迭代次数G，如果不满足就使g＝g+1；返回步骤3，如果满足则跳出循环，输出最优个体的记录，得到低碳小镇地源热泵的最优配置方案。
实施例
本发明涉及的北方某低碳能源小镇能源系统包含CHP机组1台，风电场容量为100MW，储热罐1台。选取一个自然年中冬季供暖150天的负荷数据，地源热泵设备费用：2000元/kW；运维费用0.0097元/kW·h；效率：3.5；运行年限：20年。煤炭价格为330元/吨，弃风惩罚价格为200元/MWh。其它设备参数如表1～表3所示。
表1CHP机组参数
表2CHP机组煤耗系数
表3储热罐参数
设置聚类数范围为3～10，通过本发明地源热泵容量优化配置方法进行聚类有效性函数的评估，最佳聚类数为4，该种情况下的典型风电-负荷典型组合如图2及图3所示，其中风电被赋予20％的不确定区间。对于所选的150天的风电-负荷场景，分别划分到4种典型场景的天数分别为35、42、10、63，因此在计算运行成本时，每个典型日成本乘以对应天数，相加之后即为总的运行成本。
设置3种场景，分别是不考虑不确定性的双层优化(方法1)；考虑运行优化的鲁棒优化(方法2)以及本发明所提计及源荷双侧不确定性的区间优化方法(方法3)。设置风电出力和负荷的不确定区间均为20％，对比3种场景下系统运行成本、弃风指标和总成本的对比，如表4所示，其中各结果均取4组典型风电出力-负荷曲线的结果的平均值。
表4不同方法优化配置结果对比
由表4可见，本发明所提方法与其它两种方法相比，可有效降低年运行成本，分别减少了750.598万元和300.608万元。此外，方法2和方法3均远低于不考虑不确定性的方法1，而本发明所提方法与采用鲁棒优化的方法相比可更有效的消纳弃风，这是因为鲁棒优化方法存在结果偏保守的缺点。图4a至图4d为20％不确定水平下4种典型日组合的区间优化功率平衡图，其中典型日3对热泵的利用较低，这是因为该种场景下全天风电出力较小，由典型场景划分实验可知，属于这类场景的只有10天。表5为区间优化和鲁棒优化的运行成本对比。
表5 10％～30％不确定度下运行成本对比
表5显示了在各种风电场不确定条件下，区间优化得到的总系统成本区间和鲁棒优化得到的总系统成本。由表5中的区域优化成本区间即可发现,风电不确定性越大,其总成本区间也越大。这就表明了,风电的不确定性将直接影响综合能源系统及风电的经济运行。表5中的总成本范围信息可以帮助系统调度器迅速地得到在特定风能区间情况下的系统总成本范围信息,并给出了风电不确定性对系统优化结果所影响的具体数据。另外,表5中利用鲁棒控制优化得到的总系统成本范围类似于通过区域设计得到的总成本范围的最大值。所以,与区域设计比较,利用鲁棒控制优化解的结论是比较保守的。
本发明未述及之处均适用于现有技术。
显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明的保护范围之中。
Claims (8)
1.一种低碳小镇地源热泵容量优化配置方法，其特征在于，
整理低碳小镇以往风电出力及负荷变化数据；
采用k-means聚类与Frank Copula函数构建电源侧和负荷侧的联合场景，其具体方法：
步骤1、通过整理分析过往风电出力以及负荷变化的数据，求出在规划年限下的风电以及负荷的边际概率密度分布；
步骤2、采用Frank Copula函数表示风电功率与负荷之间关系；
步骤3、利用相关矩阵图得到树型结构，构造风电-负荷联合分布函数；
步骤4、对所构造的风电-负荷联合分布函数进行分时段采样，得到风电-负荷联合采样结果；对所构造的风电-负荷联合分布函数进行逆变换，得到每个时段的风电-负荷输出，最终形成考虑风电与负荷相关性的随机场景；
步骤5、采用K-means聚类算法对步骤4得到的风电与负荷相关性的随机场景进行聚类，通过多维PFS判别函数确定最佳聚类个数，聚类中心为典型风电-负荷曲线组合，得到典型的风电电源侧和负荷侧的联合场景；
通过伪F-统计确定最优典型场景；
以综合成本为优化目标，建立优化配置目标函数；
引入区间数学，并利用量子衍生算法对所建立的目标函数进行求解，确定地源热泵最优配置方案。
2.根据权利要求1所述的低碳小镇地源热泵容量优化配置方法，其特征在于，所述的通过伪F-统计确定最优典型场景，其具体方法是：
对于P(P≥1)维变量的样本，定义一个伪F统计比率，如式(1)所示：
式中：分别为矩阵/>的迹；
m为样本数；k为类数；
分别为P维变量样本的类间和类内散布矩阵，其表达式如式(2)所示：
式中：xj为第j个样本向量；
vi为第i类Ci的聚类中心；
μij的表达式如式(3)所示：
对于式(1)，表征了类内距离，/>表征了类间距离，聚类数目k的增加，类内的距离呈负相关，类间的距离呈正相关；(m-k)/(k-1)与k呈负相关；
结合上述式(1)(2)(3)，求得PFS为最大值时的最大聚类数K，即得到最优典型场景。
3.根据权利要求1所述的低碳小镇地源热泵容量优化配置方法，其特征在于，所述的建立优化配置目标函数是以综合能源系统中地源热泵的投资成本与系统年运行成本之和最低，其表达式为：
式中：ξPl为年综合成本，ξCap为年投资成本，ξOpe为系统运行成本,D为冬季供暖期天数，D为150天；
其中的年投资成本ξCap为：
ξCap＝ξInc+ξFmc (5)
其中，ξInc为年初始投资成本，ξFmc为地源热泵固定维护成本；
CBP为地源热泵的安装容量；为地源热泵的单位容量成本；
RBP为资金收回系数；r为贴现率，取5％；l为灵活性资源生命周期，与具体灵活性资源相关；
系统运行成本ζOpe为：
其中，ξFuel为燃料成本，ξOm为运行维护成本，为弃风惩罚成本；
ωcoal为消耗煤的价格,T表示时间被分成的总数量，n表示CHP机组数量，为第i台CHP机组燃煤在t时刻的消耗量，其用发电功率的另一种形式表示：
其中，ai，bi，ci为CHP机组i的耗煤比例数，为t时刻CHP机组i的发出的电能功率；
根据泵式CHP机组的运行原理，若t时刻第i台CHP机组纯凝工况下的发电功率为随着抽气量的增加，CHP机组的发电量/>和热输出/>有如下关系：
其中，Cv为电-热耦合参数；
将式(9)代入式(8)，即可得到t时刻第i台CHP机组的煤耗量和发电功率与供热功率/>的关系：
其中，是第i台CHP机组的煤耗系数，可由ai、bi、ci和Cv计算得出；
运行维护成本ξOm为：
其中，分别为t时刻CHP机组、电锅炉和地源热泵的输出功率，分别为CHP机组、电锅炉和地源热泵单位输出功率的运行维护成本；
弃风惩罚成本为：
其中，是弃风惩罚成本系数，m代表风电机组数，/>是t时刻风机的弃风量，其等于风电预测值与风电实际发电量的实际值的差。
4.根据权利要求1所述的低碳小镇地源热泵容量优化配置方法，其特征在于，所述的优化配置目标函数中还包括优化配置约束条件，所述的优化配置约束条件中包括供能平衡约束、设备运行约束和热泵容量约束。
5.根据权利要求4所述的低碳小镇地源热泵容量优化配置方法，其特征在于，所述的供能平衡约束为：
其中，表示在t时刻第j台风机实际的发电量，/>表示在t时刻第i台CHP机组实际的发电量，/>和/>分别表示在t时刻第h台电锅炉和第k台和地源热泵的耗电量，PL,t为t时刻系统中的电负荷，/>表示在t时刻第i台CHP机组实际的产热量，/>和分别表示在t时刻第h台电锅炉和第k台和地源热泵的产热量，QL,t为t时段内系统运行时的热量值，/>为t时段内系统运行时的热负荷。
6.根据权利要求4所述的低碳小镇地源热泵容量优化配置方法，其特征在于，所述的设备运行约束包括CHP机组出力约束、地源热泵处理约束和储热罐热量约束：
所述的CHP机组出力约束为：
其中：PCHP,t和QCHP,t表示t时刻CHP机组的电功率和热功率；PCHP 和表示CHP机组的电能功率的下限和上限，QCHP 和/>表示CHP机组的热功率的下限和上限；
Cv，PCHP,D，Cm，PCHP,C和PCHP,A为CHP机组的热电耦合比例数，表示CHP机组的电出力和热出力的热电耦合联系；
所述的热泵处理约束为：
其中，QBP地源热泵热出力，地源热泵热出力最大值；
所述的储热罐热量约束为：
其中，分别为蓄热罐在运行调度周期的初始时刻t0、结束时刻tN存储的热量，/> γHS 分别为蓄热罐运行中状态的上、下限,CHS为蓄热罐的安装容量。
7.根据权利要求4所述的低碳小镇地源热泵容量优化配置方法，其特征在于，所述的地源热泵容量约束为：
其中， CBP 为地源热泵的安装容量上限和下限，根据实际情况选取。
8.根据权利要求3-7任一所述的低碳小镇地源热泵容量优化配置方法，其特征在于，所述的引入区间数学，并利用量子衍生算法对所建立的目标函数进行求解的方法是：
步骤1、输入典型日的风电功率数据、电负荷数据、热负荷数据以及各灵活性资源的参数、待定容灵活性资源的容量范围进行遗传计算的基本数据，设定好合理的种群数量、迭代次数计算基本参数；
步骤2、将目标函数中待求解的变量进行染色体实数编码，设定迭代次数g＝1，随机生成的初始种群数量为N；
步骤3、对于种群中的每个个体，将该个体对应的配置方案传递至运行层，确定运行层运行模型，在满足运行层运行约束条件的情况下，得到该配置方案下的最优运行成本；
步骤4、将运行层的最优值结果返回至目标函数的投资部分进行迭代，计算出每一个体对应的适应度值；
步骤5、保留种群中得到的最优个体，同时记录最优个体所对应的适应度值；
步骤6、进行种群的选择、交叉以及变异等操作，进而生成新种群，经过迭代选择、交叉和变异后的种群个体可能不满足所设定的范围，故需要进行检验操作，若满足则进行后续步骤，不满足则需重新进行选择、交叉和变异；
步骤7、判断此时是否达到了最大迭代次数G，如果不满足就使g＝g+1；返回步骤3，如果满足则跳出循环，输出最优个体的记录，得到低碳小镇地源热泵的最优配置方案。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211363685.7A CN115660187B (zh) | 2022-11-02 | 2022-11-02 | 一种低碳小镇地源热泵容量优化配置方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211363685.7A CN115660187B (zh) | 2022-11-02 | 2022-11-02 | 一种低碳小镇地源热泵容量优化配置方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115660187A CN115660187A (zh) | 2023-01-31 |
| CN115660187B true CN115660187B (zh) | 2024-04-30 |
Family
ID=84996084
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211363685.7A Active CN115660187B (zh) | 2022-11-02 | 2022-11-02 | 一种低碳小镇地源热泵容量优化配置方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115660187B (zh) |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6021402A (en) * | 1997-06-05 | 2000-02-01 | International Business Machines Corporaiton | Risk management system for electric utilities |
| JP2005223963A (ja) * | 2004-02-03 | 2005-08-18 | Toho Gas Co Ltd | コージェネレーションシステムの運転制御システム |
| GB201409030D0 (en) * | 2014-05-21 | 2014-07-02 | Stephens David H | Energy efficient heat pumps for warm dry homes |
| KR20190023791A (ko) * | 2017-08-30 | 2019-03-08 | 인천대학교 산학협력단 | 신재생에너지와 부하의 불확실성을 고려한 마이크로그리드의 가상등가부하변동밴드 적용에 의한 단순화된 강인최적운영 시스템 및 그 방법 |
| CN109978404A (zh) * | 2019-04-04 | 2019-07-05 | 国网经济技术研究院有限公司 | 计及不确定性因素与场景削减的输电通道扩展规划方法 |
| CN110348605A (zh) * | 2019-06-12 | 2019-10-18 | 广东工业大学 | 基于改进遗传算法的微电网经济运行优化方法 |
| CN110380444A (zh) * | 2019-06-26 | 2019-10-25 | 东南大学 | 一种基于变结构Copula的多场景下分散式风电有序接入电网的容量规划方法 |
| CN111144668A (zh) * | 2020-01-19 | 2020-05-12 | 河海大学 | 考虑场景模拟的综合能源系统随机优化模型的建立方法 |
| CN112234655A (zh) * | 2020-07-30 | 2021-01-15 | 国网甘肃省电力公司电力科学研究院 | 一种送端电网安全稳定的源网协调调峰优化方法 |
| CN112434915A (zh) * | 2020-11-09 | 2021-03-02 | 沈阳工程学院 | 面向弃风消纳的区域综合能源系统灵活性优化配置方法 |
| CN113078684A (zh) * | 2020-12-09 | 2021-07-06 | 国网辽宁省电力有限公司锦州供电公司 | 一种基于双层优化的区域能源社区规划方法 |
| CN113343478A (zh) * | 2021-06-24 | 2021-09-03 | 东北电力大学 | 考虑不确定性和需求响应的独立微网容量优化配置方法 |
| WO2021203738A1 (zh) * | 2020-04-08 | 2021-10-14 | 华北电力大学 | 考虑需求侧资源分层分散控制的配电系统可靠性计算方法 |
| CN113592133A (zh) * | 2021-05-06 | 2021-11-02 | 深圳第三代半导体研究院 | 一种能量枢纽优化配置方法及系统 |
| CN113987744A (zh) * | 2021-09-15 | 2022-01-28 | 国网吉林省电力有限公司松原供电公司 | 一种考虑风电不确定性的综合能源系统储能优化方法 |
| CN114254799A (zh) * | 2021-11-05 | 2022-03-29 | 国网浙江省电力有限公司嘉兴供电公司 | 一种碳电协同的消纳方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9058443B2 (en) * | 2012-07-17 | 2015-06-16 | International Business Machines Corporation | Planning economic energy dispatch in electrical grid under uncertainty |
-
2022
- 2022-11-02 CN CN202211363685.7A patent/CN115660187B/zh active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6021402A (en) * | 1997-06-05 | 2000-02-01 | International Business Machines Corporaiton | Risk management system for electric utilities |
| JP2005223963A (ja) * | 2004-02-03 | 2005-08-18 | Toho Gas Co Ltd | コージェネレーションシステムの運転制御システム |
| GB201409030D0 (en) * | 2014-05-21 | 2014-07-02 | Stephens David H | Energy efficient heat pumps for warm dry homes |
| KR20190023791A (ko) * | 2017-08-30 | 2019-03-08 | 인천대학교 산학협력단 | 신재생에너지와 부하의 불확실성을 고려한 마이크로그리드의 가상등가부하변동밴드 적용에 의한 단순화된 강인최적운영 시스템 및 그 방법 |
| CN109978404A (zh) * | 2019-04-04 | 2019-07-05 | 国网经济技术研究院有限公司 | 计及不确定性因素与场景削减的输电通道扩展规划方法 |
| CN110348605A (zh) * | 2019-06-12 | 2019-10-18 | 广东工业大学 | 基于改进遗传算法的微电网经济运行优化方法 |
| CN110380444A (zh) * | 2019-06-26 | 2019-10-25 | 东南大学 | 一种基于变结构Copula的多场景下分散式风电有序接入电网的容量规划方法 |
| CN111144668A (zh) * | 2020-01-19 | 2020-05-12 | 河海大学 | 考虑场景模拟的综合能源系统随机优化模型的建立方法 |
| WO2021203738A1 (zh) * | 2020-04-08 | 2021-10-14 | 华北电力大学 | 考虑需求侧资源分层分散控制的配电系统可靠性计算方法 |
| CN112234655A (zh) * | 2020-07-30 | 2021-01-15 | 国网甘肃省电力公司电力科学研究院 | 一种送端电网安全稳定的源网协调调峰优化方法 |
| CN112434915A (zh) * | 2020-11-09 | 2021-03-02 | 沈阳工程学院 | 面向弃风消纳的区域综合能源系统灵活性优化配置方法 |
| CN113078684A (zh) * | 2020-12-09 | 2021-07-06 | 国网辽宁省电力有限公司锦州供电公司 | 一种基于双层优化的区域能源社区规划方法 |
| CN113592133A (zh) * | 2021-05-06 | 2021-11-02 | 深圳第三代半导体研究院 | 一种能量枢纽优化配置方法及系统 |
| CN113343478A (zh) * | 2021-06-24 | 2021-09-03 | 东北电力大学 | 考虑不确定性和需求响应的独立微网容量优化配置方法 |
| CN113987744A (zh) * | 2021-09-15 | 2022-01-28 | 国网吉林省电力有限公司松原供电公司 | 一种考虑风电不确定性的综合能源系统储能优化方法 |
| CN114254799A (zh) * | 2021-11-05 | 2022-03-29 | 国网浙江省电力有限公司嘉兴供电公司 | 一种碳电协同的消纳方法 |
Non-Patent Citations (6)
| Title |
|---|
| Power to gas: addressing renewable curtailment by converting to hydrogen;xiaohe yan 等;Frontiers in Energy;20180928;第12卷;560–568 * |
| Wind-electricity-heat Correlation and Potential Uncertainty Analysis Based on Copula Function;Tie Li等;Journal of Physics: Conference Series;20200530;1-6 * |
| 基于改进k-means聚类的风电功率典型场景在日前调度中的应用;廖攀峰;齐军;孙绥;智李;薛冬;;电工材料;20200220(第01期);全文 * |
| 基于物元可拓模型的特色小镇能源系统综合评价;黄伟;杨子力;柳思岐;;现代电力;20201010(第05期);全文 * |
| 考虑复杂级联断面的多能源协调优化技术及应用;褚云龙等;电网技术;20201005;第44卷(第10期);3656-3662 * |
| 风光发电功率与用电负荷联合预测方法研究;王函;中国博士学位论文全文数据库 工程科技Ⅱ辑;20220115(第1期);C042-152 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115660187A (zh) | 2023-01-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xu et al. | Data-driven configuration optimization of an off-grid wind/PV/hydrogen system based on modified NSGA-II and CRITIC-TOPSIS | |
| Zhang et al. | A review of urban energy systems at building cluster level incorporating renewable-energy-source (RES) envelope solutions | |
| Yan et al. | Stochastic multi-scenario optimization for a hybrid combined cooling, heating and power system considering multi-criteria | |
| CN109103926A (zh) | 基于多辐照特性年气象场景的光伏发电接纳能力计算方法 | |
| CN112290592B (zh) | 一种风光储联合发电系统的容量优化规划方法、系统及可读存储介质 | |
| CN111754037B (zh) | 区域终端集成供能系统长期负荷混合预测方法 | |
| CN112862253B (zh) | 考虑风电不确定性的电-气-热储能协同规划方法及系统 | |
| CN112994115A (zh) | 一种基于wgan场景模拟和时序生产模拟的新能源容量配置方法 | |
| CN115455835A (zh) | 一种含可再生能源的多热源联网供热系统优化运行方法 | |
| CN117272850B (zh) | 配电网安全运行调度弹性空间分析方法 | |
| CN114066315A (zh) | 一种适应多元源荷接入的配电网规划系统 | |
| CN112072643A (zh) | 一种基于深度确定性梯度策略的光-蓄系统在线调度方法 | |
| Wang et al. | Multi-attribute decision analysis for optimal design of park-level integrated energy systems based on load characteristics | |
| CN115759610A (zh) | 一种电力系统源网荷储协同的多目标规划方法及其应用 | |
| CN115115193A (zh) | 一种工业园区低碳分析及优化规划方法 | |
| Li et al. | Optimal design and operation method of integrated energy system based on stochastic scenario iteration and energy storage attribute evaluation | |
| CN113887809A (zh) | 双碳目标下配电网供需平衡方法、系统、介质及计算设备 | |
| CN113988471A (zh) | 一种微电网运行多目标优化方法 | |
| Wang et al. | Research on short‐term and mid‐long term optimal dispatch of multi‐energy complementary power generation system | |
| CN115660187B (zh) | 一种低碳小镇地源热泵容量优化配置方法 | |
| CN116128154A (zh) | 一种农业园区综合能源系统的能源优化配置方法及装置 | |
| CN111126675A (zh) | 多能互补微网系统优化方法 | |
| CN114662922B (zh) | 考虑光伏不确定性的居民需求响应潜力评估方法及系统 | |
| CN116244927A (zh) | 一种零碳岛屿能源系统优化方法、装置、设备和存储介质 | |
| Chen et al. | Refined binary particle swarm optimization and application in power system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |