CN107294087B - 一种含气象能源的综合能源系统典型场景集生成方法 - Google Patents

一种含气象能源的综合能源系统典型场景集生成方法 Download PDF

Info

Publication number
CN107294087B
CN107294087B CN201710490900.2A CN201710490900A CN107294087B CN 107294087 B CN107294087 B CN 107294087B CN 201710490900 A CN201710490900 A CN 201710490900A CN 107294087 B CN107294087 B CN 107294087B
Authority
CN
China
Prior art keywords
mentioned steps
value
meteorological
variable
distribution
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
Application number
CN201710490900.2A
Other languages
English (en)
Other versions
CN107294087A (zh
Inventor
孙宏斌
郭庆来
王彬
付学谦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201710490900.2A priority Critical patent/CN107294087B/zh
Publication of CN107294087A publication Critical patent/CN107294087A/zh
Application granted granted Critical
Publication of CN107294087B publication Critical patent/CN107294087B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种含气象能源的综合能源系统典型场景集生成方法,属于电力系统运行与分析技术领域。本发明对气象相关性进行建模,充分考虑了气象变量之间的相关性,所提算法可以显著提高气象变量随机生产模拟的准确性。采用基于瓦瑟斯坦距离的最优分位点理论对随机模拟的连续变量离散化,得到了数量较少的综合能源系统典型场景集合。本方法的最优分位点理论可以生成极端场景,采用离散变量,生成的典型场景集可以帮助运营商分析不确定性对综合能源系统运行和规划带来的问题,典型场景的分析结果可以为系统安全稳定运行提供支持。本发明准确分析了不同气象场景下的综合能源系统状态变量相关关系特征和运行特性,为运行规划提供有价值的信息。

Description

一种含气象能源的综合能源系统典型场景集生成方法
技术领域
本发明涉及一种含气象能源的综合能源系统典型场景集生成方法,属于电力系统运行与分析技术领域。
背景技术
随着能源和环境问题的日益严峻,为了提高能源的总体效率和可再生能源的消纳能力,对多种能源系统集成和互补融合的需求日益迫切。随着多种能源系统技术的发展,电力系统和供热供冷系统的耦合度不断增强,形成冷热电综合能源系统。含高渗透率气象能源的综合能源系统和传统电力系统之间有三个的显著差异:间歇性可再生能源和随机负荷比例日益提高、多种能源系统的耦合和联系日益紧密、天气变化对系统运行的多重影响越来越大。含高渗透率气象能源和气象敏感负荷的综合能源系统易受到复杂气象环境和地理环境的影响。例如,降雨对灌溉负荷影响很大;温度和湿度对降温负荷影响很大;风场的地理位置分布对风电出力特性影响很大。传统的方法采用数字模拟方法对电力系统的不确定性进行研究,通过抽样的方式实现不确定性变量的随机生产模拟。仿真中需要生成大量的随机变量样本,不利于场景分析的简便性并使得分析工作复杂和繁琐。需要研究简单和可行的综合能源系统典型气象场景集生成方法,进而准确分析不同气象场景下的综合能源系统状态变量相关关系特征和运行特性,为运行规划提供有价值的信息。
发明内容
本发明的目的是提出一种含气象能源的综合能源系统典型场景集生成方法,以克服现有电力系统不确定建模理论的不足之处,满足高渗透率气象能源的综合能源系统的分析、运行与控制要求。
本发明提出的含气象能源的综合能源系统典型场景集生成方法,包括以下步骤:
(1)从当地气象局获取供暖季节内的太阳辐射强度数据、环境温度数据和风速数据,太阳辐射强度、环境温度和风速共有M个变量;
(2)根据上述步骤(1)采集的气象变量数据,采用参数估计方法,分别计算太阳辐射强度的边缘累积概率分布函数、环境温度的边缘累积概率分布函数和风速的边缘累积概率分布函数;
(3)采用假设检验方法,分别对上述步骤(2)计算得到的边缘累积概率分布函数进行假设检验,若通过假设检验,则进行步骤(4),若未通过假设检验,则返回步骤(2);
(4)利用上述步骤(3)检验通过的各个气象变量的边缘概率分布公式,根据上述步骤(1)采集的气象变量数据,分别得到与各气象变量数据相对应的边缘累积概率值,利用所有气象变量的边缘累积概率值,得到气象变量的Copula函数,Copula函数的自变量为所有气象变量的边缘累积概率分布变量,共M个,因变量为Copula函数的概率累积分布变量,共1个;
(5)设定连续采样数量N,利用拉丁超立方体抽样方法,生成一个N×1维的数据序列,将该数据序列作为上述Copula函数概率累积分布变量的值;
(6)根据上述步骤(5)的Copula函数概率累积分布变量值和上述步骤(4)的Copula函数,计算得到N×M维的气象变量的边缘累积概率值分布矩阵,该边缘累积概率值分布矩阵的行数等于各气象变量的边缘累积概率值的个数,边缘累积概率值分布矩阵的列数等于上述步骤(1)的气象变量的个数;
(7)根据上述步骤(6)的边缘累积概率值分布矩阵和上述步骤(3)通过检验的各个气象变量的边缘概率分布函数,得到一个N×M维的气象变量模拟值的分布矩阵,该气象变量模拟值分布矩阵的行数等于各气象变量模拟值的个数,矩阵的列数等于上述步骤(1)的气象变量的个数;
(8)利用光伏发电系统的数学模型,根据上述步骤(7)的气象变量模拟值的分布矩阵中太阳辐射强度和环境温度,计算得到综合能源系统中的光伏发电功率;利用建筑物供暖热负荷的数学模型,根据上述步骤(7)的气象变量模拟值分布矩阵中的太阳辐射强度、环境温度和风速,计算得到综合能源系统中的建筑物供暖热负荷功率;
(9)利用热电联产系统的数学模型,根据上述步骤(8)的建筑物供暖热负荷功率,计算综合能源系统中的热电联产系统的发电功率;
(10)利用电力系统潮流计算方法,根据上述步骤(8)的光伏发电功率以及步骤(9)的热电联产系统的发电功率,分别计算得到与N组综合能源系统发电功率相对应的综合能源系统潮流计算结果;
(11)利用瓦瑟斯坦距离分位点方法,对上述步骤(10)的潮流计算结果中的任一连续变量进行离散化处理,得到综合能源系统的典型场景集,典型场景集包括多个非极端场景值和两个极端场景值,两个极端场景值的计算公式如下:
其中,z0为任一连续变量的最小值,zQ+1为任一连续变量的最大值,p0 d为z0的离散概率值,pQ+1 d为zQ+1的离散概率值,z1为任一连续变量的第1个分位点,zQ为任一连续变量的第Q个分位点,pc(.)是任一连续变量的概率密度分布函数。
本发明提出的本发明提出的含气象能源的综合能源系统典型场景集生成方法,其特点和效果是:本发明方法基于Copula函数理论和拉丁超立方体抽样理论对气象相关性进行建模,充分考虑了气象变量之间的相关性,本发明方法可以显著提高气象变量随机生产模拟的准确性。采用基于瓦瑟斯坦(Wasserstein)距离的最优分位点理论对随机模拟的连续变量离散化,得到了数量较少的综合能源系统典型场景集合。相比传统离散化理论,最优分位点理论可以生成极端场景,并大幅提升了计算精度。最优分位点理论采用离散变量,而不是采用连续变量抽样技术(如蒙特卡罗模拟),使本发明有别于其他已经有的研究成果。生成的典型场景集可以帮助运营商分析不确定性对综合能源系统运行和规划带来的问题,典型场景的分析结果可以为系统安全稳定运行提供支持。
具体实施方式
本发明提出的含气象能源的综合能源系统典型场景集生成方法,包括以下步骤:
(1)从当地气象局获取供暖季节内的太阳辐射强度数据、环境温度数据和风速数据,太阳辐射强度、环境温度和风速共有M个变量;
(2)根据上述步骤(1)采集的气象变量数据,采用参数估计方法,分别计算太阳辐射强度的边缘累积概率分布函数、环境温度的边缘累积概率分布函数和风速的边缘累积概率分布函数;
(3)采用假设检验方法,分别对上述步骤(2)计算得到的边缘累积概率分布函数进行假设检验,假设检验的方法可以有多种,当边缘累积概率分布函数是正态分布时,可以采用柯尔莫可洛夫-斯米洛夫检验(Kolmogorov–Smirnov test)方法检验;边缘累积概率分布函数是非正态分布时,采用贝叶斯理论进行假设检验(参考文献姜培华,范国良.几种非正态总体未知参数的贝叶斯假设检验问题[J].南通大学学报(自然科学版),2013,(01):82-86.)。若通过假设检验,则进行步骤(4),若未通过假设检验,则返回步骤(2);
(4)利用上述步骤(3)检验通过的各个气象变量的边缘概率分布公式,根据上述步骤(1)采集的气象变量数据,分别得到与各气象变量数据相对应的边缘累积概率值,利用所有气象变量的边缘累积概率值,得到气象变量的Copula函数,Copula函数的自变量为所有气象变量的边缘累积概率分布变量,共M个,因变量为Copula函数的概率累积分布变量,共1个;
(5)设定连续采样数量N,利用拉丁超立方体抽样方法,生成一个N×1维的数据序列,将该数据序列作为上述Copula函数概率累积分布变量的值;
(6)根据上述步骤(5)的Copula函数概率累积分布变量值和上述步骤(4)的Copula函数,计算得到N×M维的气象变量的边缘累积概率值分布矩阵,该边缘累积概率值分布矩阵的行数等于各气象变量的边缘累积概率值的个数,边缘累积概率值分布矩阵的列数等于上述步骤(1)的气象变量的个数;
(7)根据上述步骤(6)的边缘累积概率值分布矩阵和上述步骤(3)通过检验的各个气象变量的边缘概率分布函数,得到一个N×M维的气象变量模拟值的分布矩阵,该气象变量模拟值分布矩阵的行数等于各气象变量模拟值的个数,矩阵的列数等于上述步骤(1)的气象变量的个数;
(8)利用光伏发电系统的数学模型(光伏发电系统数学模型可以使用美国国家可再生能源实验室(NREL)研发的HOMER软件中的光伏发电系统数学模型),根据上述步骤(7)的气象变量模拟值的分布矩阵中太阳辐射强度和环境温度,计算得到综合能源系统中的光伏发电功率;利用建筑物供暖热负荷的数学模型(建筑物供暖热负荷的数学模型使用中华人民共和国国家标准《采暖通风与空调调节设计规范(GB50019-2003)》),根据上述步骤(7)的气象变量模拟值分布矩阵中的太阳辐射强度、环境温度和风速,计算得到综合能源系统中的建筑物供暖热负荷功率;
(9)利用热电联产系统的数学模型,(热电联产系统设备的生产厂家会提供具体的数学模型及运行参数),根据上述步骤(8)的建筑物供暖热负荷功率,计算综合能源系统中的热电联产系统的发电功率;
(10)利用电力系统潮流计算方法,(本发明的一个实施例中,使用了电力系统分析综合程序PSASP进行电力系统的潮流计算)根据上述步骤(8)的光伏发电功率以及步骤(9)的热电联产系统的发电功率,分别计算得到与N组综合能源系统发电功率相对应的综合能源系统潮流计算结果,综合能源系统潮流计算结果为多个连续变量,如多个功率变量和多个电压变量),其中N为潮流计算结果中的任一连续变量的值的个数;
(11)利用瓦瑟斯坦(Wasserstein)距离分位点方法,对上述步骤(10)的潮流计算结果中的任一连续变量进行离散化处理,得到综合能源系统的典型场景集,典型场景集包括多个非极端场景值和两个极端场景值,本发明的一个实施例中利用以下文献公开的内容:黎静华,韦化,莫东.含风电场最优潮流的Wait-and-See模型与最优渐近场景分析[J].中国电机工程学报,2012,(22):15-24.,计算得到非极端场景值,两个极端场景值的计算公式如下:
其中,z0为任一连续变量的最小值,zQ+1为任一连续变量的最大值,p0 d为z0的离散概率值,pQ+1 d为zQ+1的离散概率值,z1为任一连续变量的第1个分位点,zQ为任一连续变量的第Q个分位点,pc(.)是任一连续变量的概率密度分布函数。

Claims (1)

1.一种含气象能源的综合能源系统典型场景集生成方法,其特征在于该方法包括以下步骤:
(1)从当地气象局获取供暖季节内的太阳辐射强度数据、环境温度数据和风速数据,太阳辐射强度、环境温度和风速共有M个变量;
(2)根据上述步骤(1)采集的气象变量数据,采用参数估计方法,分别计算太阳辐射强度的边缘累积概率分布函数、环境温度的边缘累积概率分布函数和风速的边缘累积概率分布函数;
(3)采用假设检验方法,分别对上述步骤(2)计算得到的边缘累积概率分布函数进行假设检验,若通过假设检验,则进行步骤(4),若未通过假设检验,则返回步骤(2);
(4)利用上述步骤(3)检验通过的各个气象变量的边缘概率分布公式,根据上述步骤(1)采集的气象变量数据,分别得到与各气象变量数据相对应的边缘累积概率值,利用所有气象变量的边缘累积概率值,得到气象变量的Copula函数,Copula函数的自变量为所有气象变量的边缘累积概率分布变量,共M个,因变量为Copula函数的概率累积分布变量,共1个;
(5)设定连续采样数量N,利用拉丁超立方体抽样方法,生成一个N×1维的数据序列,将该数据序列作为上述Copula函数概率累积分布变量的值;
(6)根据上述步骤(5)的Copula函数概率累积分布变量值和上述步骤(4)的Copula函数,计算得到N×M维的气象变量的边缘累积概率值分布矩阵,该边缘累积概率值分布矩阵的行数等于各气象变量的边缘累积概率值的个数,边缘累积概率值分布矩阵的列数等于上述步骤(1)的气象变量的个数;
(7)根据上述步骤(6)的边缘累积概率值分布矩阵和上述步骤(3)通过检验的各个气象变量的边缘概率分布函数,得到一个N×M维的气象变量模拟值的分布矩阵,该气象变量模拟值分布矩阵的行数等于各气象变量模拟值的个数,矩阵的列数等于上述步骤(1)的气象变量的个数;
(8)利用光伏发电系统的数学模型,根据上述步骤(7)的气象变量模拟值的分布矩阵中太阳辐射强度和环境温度,计算得到综合能源系统中的光伏发电功率;利用建筑物供暖热负荷的数学模型,根据上述步骤(7)的气象变量模拟值分布矩阵中的太阳辐射强度、环境温度和风速,计算得到综合能源系统中的建筑物供暖热负荷功率;
(9)利用热电联产系统的数学模型,根据上述步骤(8)的建筑物供暖热负荷功率,计算综合能源系统中的热电联产系统的发电功率;
(10)利用电力系统潮流计算方法,根据上述步骤(8)的光伏发电功率以及步骤(9)的热电联产系统的发电功率,分别计算得到与N组综合能源系统发电功率相对应的综合能源系统潮流计算结果;
(11)利用瓦瑟斯坦距离分位点方法,对上述步骤(10)的潮流计算结果中的任一连续变量进行离散化处理,得到综合能源系统的典型场景集,典型场景集包括多个非极端场景值和两个极端场景值,两个极端场景值的计算公式如下:
其中,z0为任一连续变量的最小值,zQ+1为任一连续变量的最大值,p0 d为z0的离散概率值,pQ+1 d为zQ+1的离散概率值,z1为任一连续变量的第1个分位点,zQ为任一连续变量的第Q个分位点,pc(.)是任一连续变量的概率密度分布函数。
CN201710490900.2A 2017-06-23 2017-06-23 一种含气象能源的综合能源系统典型场景集生成方法 Active CN107294087B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710490900.2A CN107294087B (zh) 2017-06-23 2017-06-23 一种含气象能源的综合能源系统典型场景集生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710490900.2A CN107294087B (zh) 2017-06-23 2017-06-23 一种含气象能源的综合能源系统典型场景集生成方法

Publications (2)

Publication Number Publication Date
CN107294087A CN107294087A (zh) 2017-10-24
CN107294087B true CN107294087B (zh) 2019-08-16

Family

ID=60099644

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710490900.2A Active CN107294087B (zh) 2017-06-23 2017-06-23 一种含气象能源的综合能源系统典型场景集生成方法

Country Status (1)

Country Link
CN (1) CN107294087B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108959719B (zh) * 2018-06-08 2022-04-12 中国科学院电工研究所 交直流混合的分布式可再生能源系统测试场景模拟方法
CN109325621B (zh) * 2018-08-29 2022-04-22 华南理工大学 一种园区能源互联网两阶段优化调度控制方法
CN109378824B (zh) * 2018-11-23 2021-05-07 国网江苏省电力有限公司 基于光伏场景的协同多无功设备的电力系统电压控制方法
CN109902340B (zh) * 2019-01-20 2023-04-07 东北电力大学 一种计及复杂气象耦合特性的多源-荷联合场景生成方法
CN109871999B (zh) * 2019-02-27 2021-08-27 华北电力大学 一种综合能源系统的随机生产模拟方法及系统
CN110323740B (zh) * 2019-06-24 2020-08-21 中国农业大学 计及光伏出力和空调负荷动态相关性的经济调峰方法
CN111509768A (zh) * 2020-04-01 2020-08-07 杭州市电力设计院有限公司 一种综合能源系统的场景模拟方法及系统
CN111697572B (zh) * 2020-06-15 2021-09-17 西安交通大学 基于多阶段随机规划理论的电源及电力流结构优化方法
CN111682574B (zh) * 2020-06-18 2021-10-15 国网江苏省电力有限公司电力科学研究院 一种交直流混联系统运行场景辨识方法、存储介质和设备
CN117522012B (zh) * 2023-11-02 2024-07-16 长江水利委员会水文局 一种基于季节周期特性的径流场景生成方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105449713A (zh) * 2015-12-11 2016-03-30 天津大学 考虑分布式电源特性的有源配电网智能软开关规划方法
CN105790258A (zh) * 2016-03-15 2016-07-20 国电南瑞科技股份有限公司 基于正态Copula函数的拉丁超立方抽样法概率潮流计算方法
CN106227706A (zh) * 2016-07-25 2016-12-14 河海大学 一种多气候模式输出数据综合校正及不确定性评估方法
CN106340901A (zh) * 2016-09-27 2017-01-18 四川大学 一种基于惯性矩的最优分位点配电网可控资源配置方法
CN106655155A (zh) * 2016-10-12 2017-05-10 中国农业大学 一种考虑故障恢复时间不确定性的配电网故障恢复方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015203946A (ja) * 2014-04-14 2015-11-16 国立大学法人京都大学 ヒストグラムの重心を計算する方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105449713A (zh) * 2015-12-11 2016-03-30 天津大学 考虑分布式电源特性的有源配电网智能软开关规划方法
CN105790258A (zh) * 2016-03-15 2016-07-20 国电南瑞科技股份有限公司 基于正态Copula函数的拉丁超立方抽样法概率潮流计算方法
CN106227706A (zh) * 2016-07-25 2016-12-14 河海大学 一种多气候模式输出数据综合校正及不确定性评估方法
CN106340901A (zh) * 2016-09-27 2017-01-18 四川大学 一种基于惯性矩的最优分位点配电网可控资源配置方法
CN106655155A (zh) * 2016-10-12 2017-05-10 中国农业大学 一种考虑故障恢复时间不确定性的配电网故障恢复方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Optimal Siting and Sizing of Energy Storage System for Power Systems with Large-scale Wind Power Integration;Haoran Zhao et al.;《IEEE》;20151231;第1-6页

Also Published As

Publication number Publication date
CN107294087A (zh) 2017-10-24

Similar Documents

Publication Publication Date Title
CN107294087B (zh) 一种含气象能源的综合能源系统典型场景集生成方法
Ali et al. Statistical analysis of wind characteristics using Weibull and Rayleigh distributions in Deokjeok-do Island–Incheon, South Korea
Maia et al. Energy and exergy analysis of the airflow inside a solar chimney
Meschede et al. Classification of global island regarding the opportunity of using RES
Laoun et al. Modeling of solar photovoltaic-polymer electrolyte membrane electrolyzer direct coupling for hydrogen generation
Fu et al. Probabilistic power flow analysis considering the dependence between power and heat
Zhang et al. Study on the operation of small rooftop wind turbines and its effect on the wind environment in blocks
Abdullahi et al. Modelling and performance analysis of a silicon PV module
Kaleshwarwar et al. Assessment of the solar energy potential of diverse urban built forms in Nagpur, India
Salmi et al. Desing and simulation of an autonomous 12.6 kW solar plant in the Algeria’s M’sila region using PVsyst software
Herlambang et al. Numerical estimation of photovoltaic–electrolyzer system performance on the basis of a weather database
CN104239662A (zh) 一种定量评价区域风能资源开发潜力的方法
Namjoo et al. Exergy performance analysis of solar photovoltaic thermal (PV/T) air collectors in terms of exergy losses
Zhang et al. Multi-site measurement for energy application of small distributed wind farm in complex mountainous areas
Wang et al. Deciphering the sensitivity of urban canopy air temperature to anthropogenic heat flux with a forcing-feedback framework
Bansod et al. Solar chimney power plant-A review
Kuscu et al. The effect of flow rate on small solar chimney performance
Sainthiya et al. Thermal modeling and performance analysis of a hybrid photovoltaic/thermal system under combined surface water cooling in winter season: an experimental approach
Hakimi et al. Wind environment at a roof-mounted wind turbine on a peaked roof building
Khosravi et al. An artificial intelligence based-model for heat transfer modeling of 5G smart poles
Chaudhry et al. Computational analysis to factor wind into the design of an architectural environment
Roshan et al. Effects of climate change on wind energy production in Iran
Forest et al. An optimized platform for performance evaluation of solar battery chargers
Akkaya et al. Analysis and applicability of Mersin region wind speed data with artificial neural networks
Awino Design-integrated Urban Heat Island analysis tool and workflow: development and application

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