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

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

Info

Publication number
CN112465271B
CN112465271B CN202011490337.7A CN202011490337A CN112465271B CN 112465271 B CN112465271 B CN 112465271B CN 202011490337 A CN202011490337 A CN 202011490337A CN 112465271 B CN112465271 B CN 112465271B
Authority
CN
China
Prior art keywords
battery
energy storage
neuron
type
input
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
CN202011490337.7A
Other languages
English (en)
Other versions
CN112465271A (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
Xinyuan Zhichu Energy Development Beijing 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, Xinyuan Zhichu Energy Development Beijing 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 GDA0004150863410000061
从j=1起始,计算第j个输入神经元nputin,j与每一个输出神经元权重向量Wq T的欧氏距离,此处欧氏距离计算公式:
Figure GDA0004150863410000062
公式中,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 GDA0004150863410000071
对电池特征的融合向量中的全部元素使用上面
Figure GDA0004150863410000072
公式进行计算后,再对全部元素取倒数,公式如下:
Figure GDA0004150863410000073
将对应于第i类电池类型的电池特征的权重向量Δkfuse.i作为层次分析法中电池类型Si对应的各个储能电池特征的输入;
步骤3:将电池特征的权重向量作为层次分析法中电池特征的初始权重,使用层次分析法计算出在储能平抑风电波动场景下储能电池四种特征的重要性程度,作为电池选型依据;包括如下具体步骤:
步骤3.1:将四类电池类型的权重向量WB=(WB1,WB2,WB3,WB4)输入层次分析法中,依照此权重向量生成判断矩阵A,具体生成公式如下:
Figure GDA0004150863410000081
将步骤2中使用SOM以及计算得出的对应于第i类储能电池类型的储能电池特征的权重向量Δkfuse.i作为层次分析法中的输入向量,四类电池类型共得到四个判断矩阵Y1,Y2,Y3,Y4;其生成公式如下:
Figure GDA0004150863410000082
步骤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 GDA0004150863410000091
公式中,λ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 GDA0004150863410000092
若一致性比率小于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),
使用特征矩阵Gcombo与对应A矩阵最大特征根λAmax的特征向量MλAmax计算得到重要性权重向量E=(e1,e2,e3,e4)。E矩阵的计算使用如下公式计算:
E=GCOMBO×MλAmax
其中,e1即代表第一种电池特征的权重,此权重的大小代表在储能平抑风电波动的应用场景下,第一种电池特征的重要程度,权重最大到最小代表此种特征对储能平抑风电波动工况的重要性由大到小;
步骤4:根据储能平抑风电波动场景下储能电池四种特征的权重,对此种场景下使用的电池进行选型;权重越大的电池特征,说明其再储能平抑风电波动场景下的重要性越大;因此,选择权重最大的电池特征以及其次的电池特征均比较优异的电池类型。
本发明不局限于上述实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (1)

1.一种面向储能平抑风电波动场景的储能电池选型方法,其特征在于,包括如下步骤:
步骤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 FDA0004150863400000021
从j=1起始,计算第j个输入神经元nputin,j与每一个输出神经元权重向量Wq T的欧氏距离,此处欧氏距离计算公式:
Figure FDA0004150863400000022
公式中,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,pfusei2,pfuse,i3,pfuse,i4);对四个电池类型Si的样本集Ki分别使用SOM模型进行特征融合,得到四个电池特征的融合向量kfuse,1,kfuse,2,kfuse3,kfuse,4
步骤2.5:对于第i类储能电池类型Si,采集其出厂时四种电池特征的额定数据Ni1,Ni2,Ni3,Ni4,即第i个电池类型对应的电池特征的额定数据;第i种电池类型的第一种电池特征的额定数据Ni1对融合向量kfuse,i中的每个元素计算公式如下:
Figure FDA0004150863400000031
对电池特征的融合向量中的全部元素使用公式
Figure FDA0004150863400000032
进行计算后,再对全部元素取倒数,公式如下:
Figure FDA0004150863400000033
将对应于第i类电池类型的电池特征的权重向量Δkfuse.i作为层次分析法中电池类型Si对应的各个储能电池特征的输入;
步骤3:将电池特征的权重向量作为层次分析法中电池特征的初始权重,使用层次分析法计算出在储能平抑风电波动场景下储能电池四种特征的重要性程度,作为电池选型依据;包括如下具体步骤:
步骤3.1:将四类电池类型的权重向量WB=(WB1,WB2,WB3,WB4)输入层次分析法中,依照此权重向量生成判断矩阵A,具体生成公式如下:
Figure FDA0004150863400000041
将步骤2中使用SOM以及计算得出的对应于第i类储能电池类型的储能电池特征的权重向量Δkfuse.i作为层次分析法中的输入向量,四类电池类型共得到四个判断矩阵Y1,Y2,Y3,Y4;其生成公式如下:
Figure FDA0004150863400000042
步骤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 FDA0004150863400000051
公式中,λ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 FDA0004150863400000052
若一致性比率小于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),
使用特征矩阵Gcombo与对应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 CN112465271A (zh) 2021-03-09
CN112465271B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113219355B (zh) * 2021-03-29 2022-04-08 安徽江淮汽车集团股份有限公司 电池选型方法、装置、设备及存储介质
CN117933764B (zh) * 2024-03-22 2024-08-20 江西师范大学 基于可持续发展指数的区域生态安全格局构建方法及系统
CN118381466B (zh) * 2024-06-21 2024-09-17 江苏协航能源科技有限公司 基于大数据的光伏设备运行监测方法

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
Selection of storage energy technologies in a power quality scenario - the AHP and the fuzzy logic;Alexandre Barin等;2009 35th Annual Conference of IEEE Industrial Electronics;第3615-3620页 *
不同场景下基于AHP-TOPSIS退役电池梯次利用综合评价;吴威;唐雨晨;叶荣;林章岁;江岳文;温步瀛;;电网与清洁能源(第04期);第115-122页 *
基于区间层次分析法的电化学储能选型方案;李建林;马会萌;田春光;惠东;;高电压技术(第09期);第2707-2714页 *
基于电压曲线的退役电池模组分选方法;王帅;尹忠东;郑重;王银顺;邹涵宇;严玉廷;;中国电机工程学报(第08期);第2691-2705页 *

Also Published As

Publication number Publication date
CN112465271A (zh) 2021-03-09

Similar Documents

Publication Publication Date Title
CN112465271B (zh) 一种面向储能平抑风电波动场景的储能电池选型方法
CN106532778B (zh) 一种计算分布式光伏并网最大准入容量的方法
CN108205717A (zh) 一种光伏发电功率多时间尺度预测方法
CN111008728A (zh) 一种用于分布式光伏发电系统短期出力的预测方法
CN105207253B (zh) 考虑风电及频率不确定性的agc随机动态优化调度方法
CN104732300B (zh) 一种基于模糊分区理论的神经网络风功率短期预测方法
CN109599872B (zh) 基于堆栈降噪自动编码器的电力系统概率潮流计算方法
CN108960491A (zh) 基于rbf神经网络的光伏发电量预测方法
CN110009141B (zh) 基于sdae特征提取和svm分类模型的爬坡事件预测方法及系统
CN110110434B (zh) 一种概率潮流深度神经网络计算的初始化方法
CN110212551B (zh) 基于卷积神经网络的微网无功自动控制方法
CN106777521B (zh) 基于双链量子遗传算法的发电机组涉网参数优化方法
CN111105005B (zh) 一种风电功率预测方法
CN118134290B (zh) 一种基于改进anfis的光伏阵列运行状态评估及故障诊断方法
CN108694475B (zh) 基于混合模型的短时间尺度光伏电池发电量预测方法
CN107330248A (zh) 一种基于改进型神经网络的短期风功率预测方法
CN112418504B (zh) 一种基于混合变量选择优化深度信念网络风速预测方法
CN114094599A (zh) 一种多站融合参与调峰调频调压潜力评估方法及装置
CN117875752A (zh) 基于自组织映射决策树的电力系统灵活运行域评估方法
CN117117827A (zh) 一种基于卷积神经网络的新型配电网状态估计方法
CN115864491A (zh) 一种处于拓扑不明状态台区的电压灵敏度拟合方法
CN114781244A (zh) 一种风电场内分群与参数优化方法
CN112531725B (zh) 一种对静止无功发生器的参数进行识别的方法及系统
Cho et al. A variable step size incremental conductance MPPT of a photovoltaic system using DC-DC converter with direct control scheme
CN104850914B (zh) 一种基于特征建模的新能源发电量预测方法

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

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.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant