CN112465271A - 一种面向储能平抑风电波动场景的储能电池选型方法 - Google Patents

一种面向储能平抑风电波动场景的储能电池选型方法 Download PDF

Info

Publication number
CN112465271A
CN112465271A CN202011490337.7A CN202011490337A CN112465271A CN 112465271 A CN112465271 A CN 112465271A CN 202011490337 A CN202011490337 A CN 202011490337A CN 112465271 A CN112465271 A CN 112465271A
Authority
CN
China
Prior art keywords
battery
energy storage
type
storage battery
weight
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.)
Granted
Application number
CN202011490337.7A
Other languages
English (en)
Other versions
CN112465271B (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.)
Beijing Lianzhi Huineng Technology Co ltd
Xinyuan Zhichu Energy Development Beijing Co ltd
North China University of Technology
Original Assignee
Beijing Lianzhi Huineng Technology Co ltd
North China University of Technology
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 Beijing Lianzhi Huineng Technology Co ltd, North China University of Technology filed Critical Beijing Lianzhi Huineng Technology Co ltd
Priority to CN202011490337.7A priority Critical patent/CN112465271B/zh
Publication of CN112465271A publication Critical patent/CN112465271A/zh
Application granted granted Critical
Publication of CN112465271B publication Critical patent/CN112465271B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/253Fusion techniques of extracted features
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/06Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
    • G06N3/061Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using biological neurons, e.g. biological neurons connected to an integrated circuit
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Economics (AREA)
  • Mathematical Physics (AREA)
  • Computational Linguistics (AREA)
  • Human Resources & Organizations (AREA)
  • Computing Systems (AREA)
  • Strategic Management (AREA)
  • Software Systems (AREA)
  • Neurology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Evolutionary Biology (AREA)
  • Marketing (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Secondary Cells (AREA)

Abstract

本发明涉及一种面向储能平抑风电波动场景的储能电池选型方法;选型方法包括如下步骤:设定储能电池类型Si,i=1,2,…,N,组成Si类型电池的电池特征的样本集Ki;使用SOM自组织映射算法,计算出第i类储能电池类型对应的各个储能电池特征的权重;使用层次分析法计算出在储能平抑风电波动场景下的储能电池四种特征的重要性权重,作为电池选型的依据;根据特征重要性权重,对储能平抑风电波动场景下所使用的储能电池类型进行筛选,以获得更高经济效益以及运行稳定性;本发明相对于现有技术的优点是:本发明利用层次分析法与自组织映射进行储能电池的选型,为电厂建设过程中的电池选型问题提供参考。

Description

一种面向储能平抑风电波动场景的储能电池选型方法
技术领域:
本发明涉及电池储能系统,进一步涉及一种面向储能平抑风电波动场景的储能电池选型方法。
背景技术:
近些年,新能源发电在电网中的占比逐渐上升,因为其发电的高效率,极低的成本以及对于环境的友好性质而受到大力投资。然而,新能源例如风力发电也有重大的缺点,那就是风力发电受到风力大小的影响。决定风力大小的因素很多,例如天气阴晴,温度,气象等等,其随机性与波动性较强,而这也导致了风电并网时导致电网波动,需要抑制此种波动以保持电网稳态运行。因此,风电站需要加装储能电池来平抑风电并网波动。
储能是一个新兴的领域,随着新能源发电形式的比重上升,储能在交通,家用等多个领域得到了发展。伴随着储能的发展,多种多样的储能电池面世,例如锂离子电池,铅酸电池等。对于不同化学原理的储能电池,具有各不相同的电池属性,例如电池寿命,充放电电压,充放电电流等。因此,在发电厂建立储能电站时,就需要进行储能电池的选型工作。由于需要考虑到储能电池的类型,不同种类电池的不同属性情况,电池选型工作往往极为复杂。而储能电池的选择是否合适直接关系到电厂的运行效率,运行稳定性,甚至在部分电厂储能担任调频工作,对电场的运行有更大的影响。
对电厂建设时储能电池类型的选择而言,此种储能电池对电厂的应用场景是否合适,是否满足电厂对电池某些属性,例如充放电电压,充放电电流,电池运行温度等等属性的要求,需要经过细致的分析后得出结果。
发明内容:
本发明是为了解决储能电池选型问题。具体技术方案如下:
一种面向储能平抑风电波动场景的储能电池选型方法,包括如下步骤:
步骤1:设定储能电池类型Si,i=1,2,…,N,在储能平抑风电波动场景下,基于当前N种类型储能电池数据,计算每个类型储能电池各自的使用数量在N类储能电池总使用数量中的占比,同时将占比作为权重WBi,其中WBi代表Si类型储能电池的权重;对运行中每个类型电池m个时刻的四类参数分别进行抽样,得到四个参数的四个m长度的抽样序列后,组成m个子样本,每一个子样本包括同一时刻的四个参数的抽样值,组成Si类型电池的电池特征的样本集Ki,Ki=(ki1,ki2,…kim),使用kij代表样本集Ki中第j个样本,kij=(pij1,pij2,pij3,pij4),pij1至pij4分别代表此电池样本的四种特征的大小;
步骤2:使用SOM自组织映射算法,将第i类型储能电池的样本集中的样本特征进行融合,得到第i类储能电池类型样本特征的融合向量,并根据第i类储能电池类型的样本特征的融合向量计算出第i类储能电池类型对应的各个储能电池特征的权重;
步骤3:将第i类储能电池类型的电池特征的权重向量作为层次分析法中电池特征的输入向量,使用层次分析法计算出在储能平抑风电波动场景下的储能电池四种特征的重要性权重,作为电池选型的依据;
步骤4:根据步骤3之中计算得出的储能平抑风电波动场景下的储能电池特征重要性权重,对储能平抑风电波动场景下所使用的储能电池类型进行筛选,以获得更高经济效益以及运行稳定性。
优选方案,一种面向储能平抑风电波动场景的储能电池选型方法包括如下步骤:
步骤1:设定储能电池类型Si,i=1,2,…,N,在储能平抑风电波动场景下,基于当前N种类型储能电池数据,计算每个类型储能电池各自的使用数量在N类储能电池总使用数量中的占比,同时将占比作为权重WBi,其中WBi代表Si类型储能电池的权重;对运行中每个类型电池m个时刻的四类参数分别进行抽样,得到四个参数的四个m长度的抽样序列后,组成m个子样本,每一个子样本包括同一时刻的四个参数的抽样值,组成Si类型电池的电池特征的样本集Ki,Ki=(ki1,ki2,…kim),使用kij代表样本集Ki中第j个样本,kij=(pij1,pij2,pij3,pij4),pij1至pij4分别代表此电池样本的四种特征的大小;
步骤2:使用SOM自组织映射算法,将第i类型储能电池的样本集中的样本特征进行融合,得到第i类储能电池类型样本特征的融合向量,并根据第i类储能电池类型的样本特征的融合向量计算出第i类储能电池类型对应的各个储能电池特征的权重;包括如下具体步骤:
步骤2.1:建立SOM神经网络模型,设定输入神经元为m个,输出神经元为Q个,为全部输出神经元赋予较小的初始权值向量,设定初始优胜邻域半径为R(0),设定初始学习率为α(0),设定迭代停止条件阈值F;
步骤2.2:使用不同类型储能电池的特征的样本集作为神经网络的输入向量,初始化全部输入向量与输出神经元初始权重向量;
步骤2.3:对于优胜神经元及其邻域神经元,使用公式计算并增益优胜神经元及其邻域神经元,使其对输入神经元的输入向量的欧氏距离缩短;
步骤2.4:不断重复上述迭代,直到学习率α(T)小于阈值F,停止迭代,得到对应于第i类储能电池类型的电池特征融合向量;
步骤2.5:计算电池特征的融合向量中各个特征对于电池特征的额定数据的误差并且取倒数,得到电池特征的权重向量,以电池特征的权重向量作为层次分析法的输入向量;
步骤3:将第i类储能电池类型的电池特征的权重向量作为层次分析法中电池特征的输入向量,使用层次分析法计算出在储能平抑风电波动场景下的储能电池四种特征的重要性权重,作为电池选型的依据;包括如下步骤:
步骤3.1:将四类电池类型对应的权重输入,依照此权重生成判断矩阵。将步骤2中使用SOM得出的四个电池特征的权重向量分别作为第i类储能电池类型的储能电池特征的输入向量,分别生成一个判断矩阵;
步骤3.2:得到电池类型的判断矩阵以及对应N种储能电池类型的四种储能电池特征的判断矩阵后,检验上述N+1个矩阵是否满足一致性要求,满足一致性要求则继续计算。
步骤3.3:根据N个不同电池类型所得出的N个矩阵的特征向量组合成一个特征矩阵,使用此特征矩阵与不同电池类型对应的权重计算,得到在储能平抑风电波动场景下的储能电池特征重要性权重;
步骤4:根据步骤3之中计算得出的储能平抑风电波动场景下的储能电池特征重要性权重,对储能平抑风电波动场景下所使用的储能电池类型进行筛选,以获得更高经济效益以及运行稳定性。
本发明相对于现有技术的优点是:本发明利用层次分析法与自组织映射进行储能电池的选型,对储能平抑风电波动场景下的储能电池属性的重要性程度进行研究,为电厂建设过程中的电池选型问题提供参考。
附图说明:
图1面向储能平抑风电波动场景的储能电池选型方法流程图。
图2是实施例步骤2中SOM自组织映射流程示意图。
图3是实施例步骤2中SOM结构示意图。
具体实施方式:
实施例:
以下结合附图对本发明作进一步详细说明。
一种面向储能平抑风电波动场景的储能电池选型方法,包括如下步骤:
步骤1:设定储能电池类型Si,i=1,2,3,4,分别代表在储能平抑风电波动场景下,磷酸铁锂电池、三元锂电池、钛酸锂电池、全钒液流电池;基于当前四种类型储能电池数据,计算每个类型储能电池各自的使用数量在四类储能电池总使用数量中的占比,同时将占比作为权重WBi,其中WBi代表Si类型储能电池的权重;对运行中每个类型电池m个时刻的四类参数分别进行抽样,得到四个参数的四个m长度的抽样序列后,组成m个子样本,每一个子样本包括同一时刻的四个参数的抽样值,组成Si类型电池的电池特征的样本集Ki,Ki=(ki1,ki2,…kim),使用kij代表样本集Ki中第j个样本,kij=(pij1,pij2,pij3,pij4),pij1至pij4分别代表此电池样本的四种特征的大小;
步骤2:使用步骤1中采集到的样本集Ki作为SOM自组织映射的输入样本集,使用SOM自组织映射对样本集进行融合,得到一个融合向量kfuse,i=(pfuse,i1,pfuse,i2,pfuse,i3,pfuse,i4);包括如下具体步骤:
步骤2.1:设定SOM输出层神经元nout,q,nout,q代表第q个输出神经元,使用Wq T=(Wq1 T,Wq2 T,Wq3 T,Wq4 T)代表第q个输出神经元的第T次迭代的权重向量,并设定输出神经元的权重向量初始值Wq 0=(Wq1 0,Wq2 0,Wq3 0,Wq4 0);设定初始优胜邻域半径R(0),设定初始学习率α(0)以及阈值F;建立SOM神经网络,输入层设定m个输入神经元nputin,j,输出层设定Q个输出神经元nout,q
步骤2.2:将第i个样本集Ki作为SOM模型的输入;每一个输入层神经元nputin,j输入一个子样本kij=(pij1,pij2,pij3,pij4),并且对全部输入神经元的输入向量kij以及输出神经元权重Wq T进行归一化,以某一输出神经元权重向量Wq T计算公式:
Figure BDA0002838904410000061
从j=1起始,计算第j个输入神经元nputin,j与每一个输出神经元权重向量Wq T的欧氏距离,此处欧氏距离计算公式:
Figure BDA0002838904410000062
公式中,Ljq T代表第j个输入神经元与第q个输出神经元之间的欧氏距离;
对于一个输入神经元来说,将其与所有的输出神经元进行欧氏距离的计算之后,得到了Q个欧氏距离,选出其中最小的欧氏距离,此欧氏距离对应的输出神经元称为第T次迭代的优胜神经元;
步骤2.3:对于优胜神经元,在其一定范围内的圆形区域称为优胜邻域,在优胜邻域内的神经元受到一次增益,使其在空间距离上更加靠近第j个输入神经元的输入向量;此增益的效果从优胜神经元向外部神经元逐渐减弱,计算公式如下:
Wq T+1=Wq T+η(T,LENTH)[kij-Wq T],
公式中,Wq T+1为增益过后的输出神经元的权值向量,作为第T+1次第q个输出神经元的权值向量,η(T,LENTH)为学习率,是一个关于迭代次数T以及优胜神经元与第q个输出神经元之间距离LENTH的函数,且随着迭代次数T的增加而不断减小,kij-Wq T为本次进行运算的输入神经元输入向量kij与第q个输出神经元权重向量的差向量;
步骤2.4:每一个输入神经元nputin,j进行一次步骤2.2与步骤2.3的计算过程并对获胜的输出经神元及其邻域神经元进行一次权重调整。全部输入神经元完成上述计算后为一次迭代;不断重复迭代过程直到学习率α(T)低于阈值F,结束迭代,完成第i个样本集Ki的数据融合,在输出层得到一个融合后的融合向量kfuse,i=(pfuse,i1,pfuse,i2,pfuse,i3,pfuse,i4);对四个电池类型Si的样本集Ki分别使用SOM模型进行特征融合,得到四个电池特征的融合向量kfuse,1,kfuse,2,kfuse,3,kfuse,4
步骤2.5:对于第i类储能电池类型Si,采集其出厂时四种电池特征的额定数据Ni1,Ni2,Ni3,Ni4,即第i个电池类型对应的电池特征的额定数据;第i种电池类型的第一种电池特征的额定数据Ni1对融合向量kfuse,i中的每个元素计算公式如下:
Figure BDA0002838904410000071
对电池特征的融合向量中的全部元素使用公式5进行计算后,再对全部元素取倒数,公式如下:
Figure BDA0002838904410000072
将对应于第i类电池类型的电池特征的权重向量Δkfuse.i作为层次分析法中电池类型Si对应的各个储能电池特征的输入;
步骤3:将电池特征的权重向量作为层次分析法中电池特征的初始权重,使用层次分析法计算出在储能平抑风电波动场景下储能电池四种特征的重要性程度,作为电池选型依据;包括如下具体步骤:
步骤3.1:将四类电池类型的权重向量WB=(WB1,WB2,WB3,WB4)输入层次分析法中,依照此权重向量生成判断矩阵A,具体生成公式如下:
Figure BDA0002838904410000081
将步骤2中使用SOM以及计算得出的对应于第i类储能电池类型的储能电池特征的权重向量Δkfuse.i作为层次分析法中的输入向量,每一类电池类型Si对应的权重向量Δkfuse.i均使用公式8进行计算,四类电池类型共得到四个判断矩阵Y1,Y2,Y3,Y4;其生成公式如下:
Figure BDA0002838904410000082
步骤3.2得到电池类型的判断矩阵A以及对应第i种电池类型的四种电池特征的判断矩阵Yi后,检验A矩阵以及Y1,Y2,Y3,Y4共五个判断矩阵是否满足一致性要求;包括如下具体步骤:
步骤3.2.1:对A矩阵以及Y1,Y2,Y3,Y4共五个判断矩阵分别计算其特征向量以及最大特征根,对判断矩阵A,将其最大特征根记为λA max,而对应于其最大特征根λA max的最大特征向量记为MλA max,对判断矩阵Yi,将其最大特征根记为λi,而对应于第i个电池类型的电池特征的权重向量kfuse,i的最大特征向量记为Mi
步骤3.2.2:对A矩阵以及Y1,Y2,Y3,Y4矩阵进行一致性校验,此处以A矩阵为例说明一致性校验计算步骤;
步骤3.2.2.1:计算一致性指标CI;CI使用如下公式计算:
Figure BDA0002838904410000091
公式中,λA max为对应判断矩阵A的最大特征根,n为判断矩阵A的维度,此处n=4;
步骤3.2.2.2:依靠专家的先验知识总结得出表格,通过查表可以确定随机一致性指标RI;
随机一致性指标数值表格
n 1 2 3 4 5 6 7 8 9
RI 0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45
步骤3.2.2.3:设定一致性比率CR,一致性比率CR使用如下公式进行计算:
Figure BDA0002838904410000092
若一致性比率小于0.1,则说明A矩阵以及Y1,Y2,Y3,Y4五个矩阵通过一致性校验,可以继续进行以下步骤;
步骤3.3:由Y1,Y2,Y3,Y4四个判断矩阵所计算得出的对应于各自最大特征根λ1,λ2,λ3,λ4的特征向量M1,M2,M3,M4,使用这四个特征向量组成一个特征矩阵Gcombo,计算公式如下:
Gcombo=(M1,M2,M3,M4),
使用特征矩阵Mcombo与对应A矩阵最大特征根λA max的特征向量MλA max计算得到重要性权重向量E=(e1,e2,e3,e4)。E矩阵的计算使用如下公式计算:
E=GCOMBO×MλA max
其中,e1即代表第一种电池特征的权重,此权重的大小代表在储能平抑风电波动的应用场景下,第一种电池特征的重要程度,权重最大到最小代表此种特征对储能平抑风电波动工况的重要性由大到小;
步骤4:根据储能平抑风电波动场景下储能电池四种特征的权重,对此种场景下使用的电池进行选型;权重越大的电池特征,说明其再储能平抑风电波动场景下的重要性越大;因此,选择权重最大的电池特征以及其次的电池特征均比较优异的电池类型。
本发明不局限于上述实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (3)

1.一种面向储能平抑风电波动场景的储能电池选型方法,其特征在于,包括如下步骤:
步骤1:设定储能电池类型Si,i=1,2,…,N,在储能平抑风电波动场景下,基于当前N种类型储能电池数据,计算每个类型储能电池各自的使用数量在N类储能电池总使用数量中的占比,同时将占比作为权重WBi,其中WBi代表Si类型储能电池的权重;对运行中每个类型电池m个时刻的四类参数分别进行抽样,得到四个参数的四个m长度的抽样序列后,组成m个子样本,每一个子样本包括同一时刻的四个参数的抽样值,组成Si类型电池的电池特征的样本集Ki,Ki=(ki1,ki2,…kim),使用kij代表样本集Ki中第j个样本,kij=(pij1,pij2,pij3,pij4),pij1至pij4分别代表此电池样本的四种特征的大小;
步骤2:使用SOM自组织映射算法,将第i类型储能电池的样本集中的样本特征进行融合,得到第i类储能电池类型样本特征的融合向量,并根据第i类储能电池类型的样本特征的融合向量计算出第i类储能电池类型对应的各个储能电池特征的权重;
步骤3:将第i类储能电池类型的电池特征的权重向量作为层次分析法中电池特征的输入向量,使用层次分析法计算出在储能平抑风电波动场景下的储能电池四种特征的重要性权重,作为电池选型的依据;
步骤4:根据步骤3之中计算得出的储能平抑风电波动场景下的储能电池特征重要性权重,对储能平抑风电波动场景下所使用的储能电池类型进行筛选,以获得更高经济效益以及运行稳定性。
2.根据权利要求1所述一种面向储能平抑风电波动场景的储能电池选型方法,其特征在于,包括如下步骤:
步骤1:设定储能电池类型Si,i=1,2,…,N,在储能平抑风电波动场景下,基于当前N种类型储能电池数据,计算每个类型储能电池各自的使用数量在N类储能电池总使用数量中的占比,同时将占比作为权重WBi,其中WBi代表Si类型储能电池的权重;对运行中每个类型电池m个时刻的四类参数分别进行抽样,得到四个参数的四个m长度的抽样序列后,组成m个子样本,每一个子样本包括同一时刻的四个参数的抽样值,组成Si类型电池的电池特征的样本集Ki,Ki=(ki1,ki2,…kim),使用kij代表样本集Ki中第j个样本,kij=(pij1,pij2,pij3,pij4),pij1至pij4分别代表此电池样本的四种特征的大小;
步骤2:使用SOM自组织映射算法,将第i类型储能电池的样本集中的样本特征进行融合,得到第i类储能电池类型样本特征的融合向量,并根据第i类储能电池类型的样本特征的融合向量计算出第i类储能电池类型对应的各个储能电池特征的权重;包括如下具体步骤:
步骤2.1:建立SOM神经网络模型,设定输入神经元为m个,输出神经元为Q个,为全部输出神经元赋予较小的初始权值向量,设定初始优胜邻域半径为R(0),设定初始学习率为α(0),设定迭代停止条件阈值F;
步骤2.2:使用不同类型储能电池的特征的样本集作为神经网络的输入向量,初始化全部输入向量与输出神经元初始权重向量;
步骤2.3:对于优胜神经元及其邻域神经元,使用公式计算并增益优胜神经元及其邻域神经元,使其对输入神经元的输入向量的欧氏距离缩短;
步骤2.4:不断重复上述迭代,直到学习率α(T)小于阈值F,停止迭代,得到对应于第i类储能电池类型的电池特征融合向量;
步骤2.5:计算电池特征的融合向量中各个特征对于电池特征的额定数据的误差并且取倒数,得到电池特征的权重向量,以电池特征的权重向量作为层次分析法的输入向量;
步骤3:将第i类储能电池类型的电池特征的权重向量作为层次分析法中电池特征的输入向量,使用层次分析法计算出在储能平抑风电波动场景下的储能电池四种特征的重要性权重,作为电池选型的依据;包括如下步骤:
步骤3.1:将四类电池类型对应的权重输入,依照此权重生成判断矩阵。将步骤2中使用SOM得出的四个电池特征的权重向量分别作为第i类储能电池类型的储能电池特征的输入向量,分别生成一个判断矩阵;
步骤3.2:得到电池类型的判断矩阵以及对应N种储能电池类型的四种储能电池特征的判断矩阵后,检验上述N+1个矩阵是否满足一致性要求,满足一致性要求则继续计算。
步骤3.3:根据N个不同电池类型所得出的N个矩阵的特征向量组合成一个特征矩阵,使用此特征矩阵与不同电池类型对应的权重计算,得到在储能平抑风电波动场景下的储能电池特征重要性权重;
步骤4:根据步骤3之中计算得出的储能平抑风电波动场景下的储能电池特征重要性权重,对储能平抑风电波动场景下所使用的储能电池类型进行筛选,以获得更高经济效益以及运行稳定性。
3.根据权利要求2所述一种面向储能平抑风电波动场景的储能电池选型方法,其特征在于,包括如下步骤:
步骤1:设定储能电池类型Si,i=1,2,3,4,分别代表在储能平抑风电波动场景下,磷酸铁锂电池、三元锂电池、钛酸锂电池、全钒液流电池;基于当前四种类型储能电池数据,计算每个类型储能电池各自的使用数量在四类储能电池总使用数量中的占比,同时将占比作为权重WBi,其中WBi代表Si类型储能电池的权重;对运行中每个类型电池m个时刻的四类参数分别进行抽样,得到四个参数的四个m长度的抽样序列后,组成m个子样本,每一个子样本包括同一时刻的四个参数的抽样值,组成Si类型电池的电池特征的样本集Ki,Ki=(ki1,ki2,…kim),使用kij代表样本集Ki中第j个样本,kij=(pij1,pij2,pij3,pij4),pij1至pij4分别代表此电池样本的四种特征的大小;
步骤2:使用步骤1中采集到的样本集Ki作为SOM自组织映射的输入样本集,使用SOM自组织映射对样本集进行融合,得到一个融合向量kfuse,i=(pfuse,i1,pfuse,i2,pfuse,i3,pfuse,i4);包括如下具体步骤:
步骤2.1:设定SOM输出层神经元nout,q,nout,q代表第q个输出神经元,使用Wq T=(Wq1 T,Wq2 T,Wq3 T,Wq4 T)代表第q个输出神经元的第T次迭代的权重向量,并设定输出神经元的权重向量初始值Wq 0=(Wq1 0,Wq2 0,Wq3 0,Wq4 0);设定初始优胜邻域半径R(0),设定初始学习率α(0)以及阈值F;建立SOM神经网络,输入层设定m个输入神经元nputin,j,输出层设定Q个输出神经元nout,q
步骤2.2:将第i个样本集Ki作为SOM模型的输入;每一个输入层神经元nputin,j输入一个子样本kij=(pij1,pij2,pij3,pij4),并且对全部输入神经元的输入向量kij以及输出神经元权重Wq T进行归一化,以某一输出神经元权重向量Wq T计算公式:
Figure FDA0002838904400000041
从j=1起始,计算第j个输入神经元nputin,j与每一个输出神经元权重向量Wq T的欧氏距离,此处欧氏距离计算公式:
Figure FDA0002838904400000042
公式中,Ljq T代表第j个输入神经元与第q个输出神经元之间的欧氏距离;
对于一个输入神经元来说,将其与所有的输出神经元进行欧氏距离的计算之后,得到了Q个欧氏距离,选出其中最小的欧氏距离,此欧氏距离对应的输出神经元称为第T次迭代的优胜神经元;
步骤2.3:对于优胜神经元,在其一定范围内的圆形区域称为优胜邻域,在优胜邻域内的神经元受到一次增益,使其在空间距离上更加靠近第j个输入神经元的输入向量;此增益的效果从优胜神经元向外部神经元逐渐减弱,计算公式如下:
Wq T+1=Wq T +η(T,LENTH)[kij-Wq T],
公式中,Wq T+1为增益过后的输出神经元的权值向量,作为第T+1次第q个输出神经元的权值向量,η(T,LENTH)为学习率,是一个关于迭代次数T以及优胜神经元与第q个输出神经元之间距离LENTH的函数,且随着迭代次数T的增加而不断减小,kij-Wq T为本次进行运算的输入神经元输入向量kij与第q个输出神经元权重向量的差向量;
步骤2.4:每一个输入神经元nputin,j进行一次步骤2.2与步骤2.3的计算过程并对获胜的输出经神元及其邻域神经元进行一次权重调整。全部输入神经元完成上述计算后为一次迭代;不断重复迭代过程直到学习率α(T)低于阈值F,结束迭代,完成第i个样本集Ki的数据融合,在输出层得到一个融合后的融合向量kfuse,i=(pfuse,i1,pfuse,i2,pfuse,i3,pfuse,i4);对四个电池类型Si的样本集Ki分别使用SOM模型进行特征融合,得到四个电池特征的融合向量kfuse,1,kfuse,2,kfuse,3,kfuse,4
步骤2.5:对于第i类储能电池类型Si,采集其出厂时四种电池特征的额定数据Ni1,Ni2,Ni3,Ni4,即第i个电池类型对应的电池特征的额定数据;第i种电池类型的第一种电池特征的额定数据Ni1对融合向量kfuse,i中的每个元素计算公式如下:
Figure FDA0002838904400000061
对电池特征的融合向量中的全部元素使用公式5进行计算后,再对全部元素取倒数,公式如下:
Figure FDA0002838904400000062
将对应于第i类电池类型的电池特征的权重向量Δkfuse.i作为层次分析法中电池类型Si对应的各个储能电池特征的输入;
步骤3:将电池特征的权重向量作为层次分析法中电池特征的初始权重,使用层次分析法计算出在储能平抑风电波动场景下储能电池四种特征的重要性程度,作为电池选型依据;包括如下具体步骤:
步骤3.1:将四类电池类型的权重向量WB=(WB1,WB2,WB3,WB4)输入层次分析法中,依照此权重向量生成判断矩阵A,具体生成公式如下:
Figure FDA0002838904400000063
将步骤2中使用SOM以及计算得出的对应于第i类储能电池类型的储能电池特征的权重向量Δkfuse.i作为层次分析法中的输入向量,每一类电池类型Si对应的权重向量Δkfuse.i均使用公式8进行计算,四类电池类型共得到四个判断矩阵Y1,Y2,Y3,Y4;其生成公式如下:
Figure FDA0002838904400000071
步骤3.2得到电池类型的判断矩阵A以及对应第i种电池类型的四种电池特征的判断矩阵Yi后,检验A矩阵以及Y1,Y2,Y3,Y4共五个判断矩阵是否满足一致性要求;包括如下具体步骤:
步骤3.2.1:对A矩阵以及Y1,Y2,Y3,Y4共五个判断矩阵分别计算其特征向量以及最大特征根,对判断矩阵A,将其最大特征根记为λAmax,而对应于其最大特征根λAmax的最大特征向量记为MλAmax,对判断矩阵Yi,将其最大特征根记为λi,而对应于第i个电池类型的电池特征的权重向量kfuse,i的最大特征向量记为Mi
步骤3.2.2:对A矩阵以及Y1,Y2,Y3,Y4矩阵进行一致性校验,此处以A矩阵为例说明一致性校验计算步骤;
步骤3.2.2.1:计算一致性指标CI;CI使用如下公式计算:
Figure FDA0002838904400000072
公式中,λAmax为对应判断矩阵A的最大特征根,n为判断矩阵A的维度,此处n=4;
步骤3.2.2.2:依靠专家的先验知识总结得出表格,通过查表可以确定随机一致性指标RI;
随机一致性指标数值表格
n 1 2 3 4 5 6 7 8 9 RI 0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45
步骤3.2.2.3:设定一致性比率CR,一致性比率CR使用如下公式进行计算:
Figure FDA0002838904400000081
若一致性比率小于0.1,则说明A矩阵以及Y1,Y2,Y3,Y4五个矩阵通过一致性校验,可以继续进行以下步骤;
步骤3.3:由Y1,Y2,Y3,Y4四个判断矩阵所计算得出的对应于各自最大特征根λ1,λ2,λ3,λ4的特征向量M1,M2,M3,M4,使用这四个特征向量组成一个特征矩阵Gcombo,计算公式如下:
Gcombo=(M1,M2,M3,M4),
使用特征矩阵Mcombo与对应A矩阵最大特征根λAmax的特征向量MλAmax计算得到重要性权重向量E=(e1,e2,e3,e4)。E矩阵的计算使用如下公式计算:
E=GCOMBO×MλAmax
其中,e1即代表第一种电池特征的权重,此权重的大小代表在储能平抑风电波动的应用场景下,第一种电池特征的重要程度,权重最大到最小代表此种特征对储能平抑风电波动工况的重要性由大到小;
步骤4:根据储能平抑风电波动场景下储能电池四种特征的权重,对此种场景下使用的电池进行选型;权重越大的电池特征,说明其再储能平抑风电波动场景下的重要性越大;因此,选择权重最大的电池特征以及其次的电池特征均比较优异的电池类型。
CN202011490337.7A 2020-12-16 2020-12-16 一种面向储能平抑风电波动场景的储能电池选型方法 Active CN112465271B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011490337.7A CN112465271B (zh) 2020-12-16 2020-12-16 一种面向储能平抑风电波动场景的储能电池选型方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011490337.7A CN112465271B (zh) 2020-12-16 2020-12-16 一种面向储能平抑风电波动场景的储能电池选型方法

Publications (2)

Publication Number Publication Date
CN112465271A true CN112465271A (zh) 2021-03-09
CN112465271B CN112465271B (zh) 2023-05-23

Family

ID=74803005

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011490337.7A Active CN112465271B (zh) 2020-12-16 2020-12-16 一种面向储能平抑风电波动场景的储能电池选型方法

Country Status (1)

Country Link
CN (1) CN112465271B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113219355A (zh) * 2021-03-29 2021-08-06 安徽江淮汽车集团股份有限公司 电池选型方法、装置、设备及存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104166790A (zh) * 2014-07-24 2014-11-26 广东电网公司电力科学研究院 基于逼近理想解理论的锂离子动力电池性能评价方法
CN104899459A (zh) * 2015-06-16 2015-09-09 北京亿利智慧能源科技有限公司 基于层次分析法的电池性能评价方法
CN105129109A (zh) * 2015-09-30 2015-12-09 北京航空航天大学 一种基于多重分形理论和自组织映射网络的飞机副翼作动器系统健康评估方法
CN105654175A (zh) * 2015-12-24 2016-06-08 北方民族大学 一种面向轴承制造企业的零件供应商多目标优选方法
CN107437135A (zh) * 2016-05-26 2017-12-05 中国电力科学研究院 一种新型储能选型方法
CN109800950A (zh) * 2018-12-17 2019-05-24 国家电网有限公司 基于层次分析法的梯次利用电池储能系统性能评估方法
CN110929791A (zh) * 2019-11-27 2020-03-27 北京交通大学 梯次利用电池的应用场景选择方法
CN111487532A (zh) * 2020-04-09 2020-08-04 北方工业大学 一种基于层次分析法和熵值法的退役电池筛选方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104166790A (zh) * 2014-07-24 2014-11-26 广东电网公司电力科学研究院 基于逼近理想解理论的锂离子动力电池性能评价方法
CN104899459A (zh) * 2015-06-16 2015-09-09 北京亿利智慧能源科技有限公司 基于层次分析法的电池性能评价方法
CN105129109A (zh) * 2015-09-30 2015-12-09 北京航空航天大学 一种基于多重分形理论和自组织映射网络的飞机副翼作动器系统健康评估方法
CN105654175A (zh) * 2015-12-24 2016-06-08 北方民族大学 一种面向轴承制造企业的零件供应商多目标优选方法
CN107437135A (zh) * 2016-05-26 2017-12-05 中国电力科学研究院 一种新型储能选型方法
CN109800950A (zh) * 2018-12-17 2019-05-24 国家电网有限公司 基于层次分析法的梯次利用电池储能系统性能评估方法
CN110929791A (zh) * 2019-11-27 2020-03-27 北京交通大学 梯次利用电池的应用场景选择方法
CN111487532A (zh) * 2020-04-09 2020-08-04 北方工业大学 一种基于层次分析法和熵值法的退役电池筛选方法及系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ALEXANDRE BARIN等: "Selection of storage energy technologies in a power quality scenario - the AHP and the fuzzy logic", 2009 35TH ANNUAL CONFERENCE OF IEEE INDUSTRIAL ELECTRONICS *
吴威;唐雨晨;叶荣;林章岁;江岳文;温步瀛;: "不同场景下基于AHP-TOPSIS退役电池梯次利用综合评价", 电网与清洁能源 *
李建林;马会萌;田春光;惠东;: "基于区间层次分析法的电化学储能选型方案", 高电压技术 *
王帅;尹忠东;郑重;王银顺;邹涵宇;严玉廷;: "基于电压曲线的退役电池模组分选方法", 中国电机工程学报 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113219355A (zh) * 2021-03-29 2021-08-06 安徽江淮汽车集团股份有限公司 电池选型方法、装置、设备及存储介质
CN113219355B (zh) * 2021-03-29 2022-04-08 安徽江淮汽车集团股份有限公司 电池选型方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN112465271B (zh) 2023-05-23

Similar Documents

Publication Publication Date Title
CN111157899B (zh) 一种基于模型融合思想对电池soc的估计方法
CN106532778B (zh) 一种计算分布式光伏并网最大准入容量的方法
CN110941929A (zh) 一种基于ARMA和Elman神经网络联合建模的电池健康状态评估方法
CN110429637B (zh) 一种概率静态电压稳定域的可视化方法
CN107681664A (zh) 一种基于内点法的分布式电源并网优化潮流计算方法
CN110635478B (zh) 一种基于单目标的新能源接入下的输电网规划的优化方法
CN110429636B (zh) 一种静态电压稳定故障筛选与排序的方法
CN110009141B (zh) 基于sdae特征提取和svm分类模型的爬坡事件预测方法及系统
Bagheri et al. Impacts of renewable energy sources by battery forecasting on smart power systems
CN115907339A (zh) 一种基于gmm的配电网光伏与电动汽车充电站协同规划方法
CN106777521B (zh) 基于双链量子遗传算法的发电机组涉网参数优化方法
CN112465271A (zh) 一种面向储能平抑风电波动场景的储能电池选型方法
CN108694475B (zh) 基于混合模型的短时间尺度光伏电池发电量预测方法
CN114094599A (zh) 一种多站融合参与调峰调频调压潜力评估方法及装置
Wang et al. Accurate solar cell modeling via genetic Neural network-based Meta-Heuristic algorithms
Sun et al. Short-term photovoltaic power prediction modeling based on AdaBoost algorithm and Elman
CN115204048A (zh) 一种光伏发电系统模型辨识系统和方法
CN115600503A (zh) 基于单纯形-模拟退火法的多结化合物电池模型参数优化方法
CN113690941A (zh) 配置光储充微网容量的方法、装置、设备及存储介质
Sun et al. Online model identification method of vanadium redox flow battery based on time-varying forgetting factor recursive least squares
CN110571791B (zh) 一种新能源接入下的输电网规划的优化配置方法
CN113659636A (zh) 分布不确定性条件下N-k分布鲁棒机组控制方法及装置
Singsathid et al. Parameter estimation of solar PV models using self-adaptive differential evolution with dynamic mutation and pheromone strategy
CN113807019A (zh) 基于改进场景分类和去粗粒化的mcmc风电模拟方法
Li et al. Development of control strategy to increase the lifetime of grid-connected Li-ion battery

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20220506

Address after: 100144 Beijing City, Shijingshan District Jin Yuan Zhuang Road No. 5

Applicant after: NORTH CHINA University OF TECHNOLOGY

Applicant after: Beijing Lianzhi Huineng Technology Co.,Ltd.

Applicant after: Xinyuan Zhichu energy development (Beijing) Co.,Ltd.

Address before: 100144 Beijing City, Shijingshan District Jin Yuan Zhuang Road No. 5

Applicant before: NORTH CHINA University OF TECHNOLOGY

Applicant before: Beijing Lianzhi Huineng Technology Co.,Ltd.

GR01 Patent grant
GR01 Patent grant