CN111261903B - 一种基于模型的质子交换膜燃料电池阻抗在线估计方法 - Google Patents
一种基于模型的质子交换膜燃料电池阻抗在线估计方法 Download PDFInfo
- Publication number
- CN111261903B CN111261903B CN202010068829.0A CN202010068829A CN111261903B CN 111261903 B CN111261903 B CN 111261903B CN 202010068829 A CN202010068829 A CN 202010068829A CN 111261903 B CN111261903 B CN 111261903B
- Authority
- CN
- China
- Prior art keywords
- fuel cell
- model
- anode
- cathode
- water
- 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
- 239000000446 fuel Substances 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000012528 membrane Substances 0.000 title claims abstract description 34
- 239000002245 particle Substances 0.000 claims abstract description 38
- 230000002068 genetic effect Effects 0.000 claims abstract description 15
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 230000004044 response Effects 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 58
- 239000007789 gas Substances 0.000 claims description 33
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 17
- 239000001257 hydrogen Substances 0.000 claims description 15
- 229910052739 hydrogen Inorganic materials 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 13
- 238000012546 transfer Methods 0.000 claims description 13
- 239000007788 liquid Substances 0.000 claims description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 6
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000036571 hydration Effects 0.000 claims description 6
- 238000006703 hydration reaction Methods 0.000 claims description 6
- 239000001301 oxygen Substances 0.000 claims description 6
- 229910052760 oxygen Inorganic materials 0.000 claims description 6
- 239000008213 purified water Substances 0.000 claims description 6
- 230000036961 partial effect Effects 0.000 claims description 5
- 238000009792 diffusion process Methods 0.000 claims description 4
- 230000035772 mutation Effects 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 claims description 2
- 229910001882 dioxygen Inorganic materials 0.000 claims description 2
- ONUFESLQCSAYKA-UHFFFAOYSA-N iprodione Chemical compound O=C1N(C(=O)NC(C)C)CC(=O)N1C1=CC(Cl)=CC(Cl)=C1 ONUFESLQCSAYKA-UHFFFAOYSA-N 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 230000005284 excitation Effects 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000001228 spectrum Methods 0.000 abstract description 4
- 230000009467 reduction Effects 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 101100274110 Mycolicibacterium paratuberculosis (strain ATCC BAA-968 / K-10) groEL2 gene Proteins 0.000 description 1
- 101100390908 Pectobacterium carotovorum subsp. carotovorum fliN gene Proteins 0.000 description 1
- 101100457865 Rhodobacter capsulatus mopA gene Proteins 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005370 electroosmosis Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 101150077981 groEL gene Proteins 0.000 description 1
- 101150028210 groEL1 gene Proteins 0.000 description 1
- 235000003642 hunger Nutrition 0.000 description 1
- 238000001453 impedance spectrum Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 239000012495 reaction gas Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000037351 starvation Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04313—Processes for controlling fuel cells or fuel cell systems characterised by the detection or assessment of variables; characterised by the detection or assessment of failure or abnormal function
- H01M8/04537—Electric variables
- H01M8/04634—Other electric variables, e.g. resistance or impedance
- H01M8/04641—Other electric variables, e.g. resistance or impedance of the individual fuel cell
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04305—Modeling, demonstration models of fuel cells, e.g. for training purposes
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04992—Processes for controlling fuel cells or fuel cell systems characterised by the implementation of mathematical or computational algorithms, e.g. feedback control loops, fuzzy logic, neural networks or artificial intelligence
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/30—Hydrogen technology
- Y02E60/50—Fuel cells
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Chemical & Material Sciences (AREA)
- Electrochemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Chemical & Material Sciences (AREA)
- Sustainable Energy (AREA)
- Sustainable Development (AREA)
- Manufacturing & Machinery (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Medical Informatics (AREA)
- Fuzzy Systems (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Automation & Control Theory (AREA)
- Artificial Intelligence (AREA)
- Fuel Cell (AREA)
Abstract
本发明涉及一种基于模型的质子交换膜燃料电池阻抗在线估计方法，包括以下步骤：S1建立适用于实车控制器应用的动态燃料电池集总参数模型；S2采用遗传粒子群混合寻优算法对动态燃料电池集总参数模型进行参数辨识；S3将高频阶跃电流信号施加到辨识后的动态燃料电池集总参数模型中，得到高频电压响应数据；S4对时域数据进行加汉宁窗处理，并采用快速傅里叶变换对燃料电池的电压和电流进行时域到频域转换，基于变换后的电压和电流计算各频率下的燃料电池阻抗。与现有技术相比，本发明具有无需交流激励源、计算量少、速度快、减少频谱泄露等优点。
Description
技术领域
本发明涉及燃料电池技术领域，尤其是涉及一种基于模型的质子交换膜燃料电池阻抗在线估计方法。
背景技术
质子交换膜燃料电池(PEMFC)是一种以氢能为载体将化学能转化为电能的发电装置，具有反应温度低、动态响应速度快、反应效率高、功率密度高等优点，在交通领域具有广泛的应用前景。但是，燃料电池大规模商业化应用受到耐久性和可靠性制约，主要由供气不足引起的反应饥饿以及不当水管理引起的膜干和水淹造成。为进一步提升燃料电池性能和使用寿命，需要对燃料电池内部状态进行诊断并实时控制。电化学阻抗谱能够解析燃料电池内部不同电化学过程，从而可推断出燃料电池内部关键状态，如质子交换膜内水含量、多孔介质内部液态水含量、反应气体供给状态等，因此广泛用于燃料电池故障诊断。
中国专利CN103904348公开了一种用于诊断燃料电池堆阻抗的方法和系统。该发明主要是先合成不同频率的多个正弦信号，将合成的正弦信号施加于燃料电池堆，获取电压电流数据并进行傅里叶变换，从而计算阻抗。其中，合成多个正弦信号源是本发明方法的关键，因此严重依赖信号发生器和信号合成器，大大增加了硬件成本。
中国专利CN105699902公开了一种用于燃料电池诊断的阻抗测定装置及其方法，该发明主要是利用电流的变动状态对燃料电池组的电流及电压进行采样，并进行阻抗演算，但该方法主要针对一个或两个频率成分进行计算，频率范围宽度不够，高频阻抗获取受到硬件采样频率限制。
中国专利CN107482911公开了一种适用于氢燃料电池堆交流阻抗测试的DC/DC变换器，但该发明DC/DC的拓扑结构为传统boost电路，无法适用于大功率燃料电池系统。另外，中国专利CN109212431优化了DC/DC拓扑结构，但为了产生适用于燃料电池阻抗测试的激励源，硬件成本和复杂度都会增加。
因此，为了减小车载阻抗获取的成本和难度，有必要提出一种适用于实车控制器的宽频率带的在线阻抗计算方法。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于模型的质子交换膜燃料电池阻抗在线估计方法。
本发明的目的可以通过以下技术方案来实现：
一种基于模型的质子交换膜燃料电池阻抗在线估计方法，包括以下步骤：
S1建立适用于实车控制器应用的动态燃料电池集总参数模型；
S2采用遗传粒子群混合寻优算法对动态燃料电池集总参数模型进行参数辨识；
S3将高频阶跃电流信号施加到辨识后的动态燃料电池集总参数模型中，得到高频电压响应数据；
S4对时域数据进行加汉宁窗处理，并采用快速傅里叶变换对燃料电池的电压和电流进行时域到频域转换，基于变换后的电压和电流计算各频率下的燃料电池阻抗。
所述的步骤S1中，基于以下假设条件建立动态燃料电池集总参数模型，则有：
1)所有气体均视为理想气体；
2)内部气体组分的压力、温度和浓度均匀分布；
3)不考虑热力学动态过程，电池温度视为恒定不变。
所述的动态燃料电池集总参数模型包括用以描述阴极内部气体压力与外部操作条件耦合关系的阴极腔体模型、用以描述阳极内部气体压力与外部操作条件耦合关系的阳极腔体模型、用以描述膜内水传质过程的膜水合模型以及用以描述电压响应的电压电化学模型，所述的外部操作条件包括氢气质量流量、空气质量流量、阴阳极入口温度、阴阳极入口压力、阴阳极进口相对湿度、阴阳极出口压力和负载电流。
所述的阴极腔体模型具体为：
其中，
和mliq,ca分别为阴极流道内氧气、氮气、水蒸气和液态水的质量，
和
分别为各气体进入阴极腔体和排出阴极腔体的质量分数，Wca,in和Wca,out分别为进入阴极腔体和排出阴极腔体混合气的质量流量，Wliq,ca为排出阴极腔体的液态水质量，
为阴极侧水的相变率，β为膜内净水传质速率，Ist为负载电流，F为法拉第常数，
为氧气的摩尔质量，
为水的摩尔质量。
所述的阳极腔体模型具体为：
其中，
和yvap,anin分别为经过加湿器后进入电池前氢气质量分数和水蒸气质量分数，
mvap,an分别为阳极流道内氢气和水蒸气质量，Wan,in和Wan,out分别为进入阳极和排出阳极的混合气质量流量，
为氢气的摩尔质量，
为水的摩尔质量，mliq,an为阳极侧液态水质量，Wliq,an为排出阳极腔体的液态水质量，
为排出阳极腔体的氢气质量分数，yvap,anout为排出阳极水蒸气的质量分数，mphase,an为相变水质量，Ist为负载电流，F为法拉第常数，β为膜内净水传质速率。
所述的膜水合模型具体为：
其中，Acell为燃料电池的等效活性面积，F为法拉第常数，Ist为负载电流，β为膜内净水传质速率，Nv,mem为单位时间内穿透膜的水传递速率，nd为电子拖拽系数，Dw为浓差扩散系数，cv,ca和cv,an分别为阴阳极两侧水浓度，Lm为膜的厚度。
所述的电化学模型具体为：
其中，Vcell为模型估计输出电池电压，Tst为燃料电池温度，
为阳极侧氢气分压，
为阴极侧氧气分压，R为气体常数，F为法拉第常数，Lm为膜的厚度，Acell为电池的等效活性面积，λm为交换膜的含水量，θ1～θ10为待识别参数，Ist为负载电流。
所述的步骤S2中，混合遗传粒子群算法具体包括以下步骤：
定义种群大小、空间维度、遗传寻优操作和粒子群寻优操作的基本参数；
计算每个粒子的个体最优和全局最优，并记录每个粒子的初始位置；
对粒子进行粒子群算法操作，判断粒子适应度是否满足要求；
若不满足适应度要求，则对粒子进行选择、交叉、变异的遗传学操作；
将粒子再进行一次速度和位置更新，判断粒子是否满足适应度要求，若不满足则重复上述过程，直到搜索到问题最优解。
所述的粒子适应度定义为实测电压和模型输出电池电压的误差平方和。
所述的步骤S4中，各频率下的燃料电池阻抗Z的具体表达式为：
其中，|FFT(Vcell(t))|为快速傅里叶变换后的模型估计输出电池电压模，|FFT(Ist(t))|为快速傅里叶变换后的高频阶跃电流模，FFT为快速傅里叶计算，Vcell(t)为模型输出电压，Ist(t)为输入电流。
与现有技术相比，本发明具有以下优点：
一、本发明只需要事先对燃料电池模型进行辨识，不需要交流激励源装置，即可在线计算燃料电池不同频率下的阻抗，可实现车载应用并节省成本。
二、本发明采用的燃料电池模型能够描述外部操作条件对燃料电池内部状态的影响和电压的动态响应，计算量少，可嵌入于燃料电池控制器内。
三、本发明采用了混合遗传粒子群算法对燃料电池模型进行参数辨识，与传统粒子群算法相比能够减少陷入局部最优的可能，与传统遗传算法相比具有更快的计算速度。
四、本发明采用了汉宁窗对电压和电流数据进行处理，减少了频谱的泄露，采用快速傅里叶变换对燃料电池阻抗进行计算，与传统离散傅里叶变换计算相比计算速度更快。
附图说明
图1为本发明的方法流程图。
图2为本发明实例中所构建的燃料电池模型示意图。
图3为本发明实例中燃料电池模型参数辨识算法流程图。
图4为本发明实例模型辨识结果。
图5为本发明实例注入到模型的阶跃电流。
图6为本发明实例燃料电池阻抗计算结果。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
图1为本发明所构建的质子交换膜燃料电池阻抗在线估计方法，如图1所示，在本实例中，该方法包括以下步骤；
S1，建立适用于实车控制器应用的集总参数模型；
建立燃料电池模型之前，对模型做出如下假设；
1)所有气体均视为理想气体；
2)内部气体组分的压力、温度和浓度均匀分布；
3)不考虑热力学动态过程，电池温度视为恒定不变。
基于上述模型假设，建立燃料电池模型，图2为本发明实例中燃料电池模型的结构框图，模型包括电化学模型、阴极腔体模型、阳极腔体模型和膜水合模型。
阴极腔体模型主要基于质量守恒定律描述燃料电池阴极内部气体压力与外部操作条件耦合关系，其主要方程式如下：
其中，
mvap,ca和mliq,ca分别为阴极流道内氧气、氮气、水蒸气和液态水的质量；
yvap,in和
yvap,out分别为各气体进入阴极腔体和排出阴极腔体的质量分数；Wca,in和Wca,out分别为进入阴极腔体和排出阴极腔体混合气的质量流量，单位为kg/s；Wliq,ca为排出阴极腔体的液态水质量，
为阴极侧水的相变率；β为膜内净水传质速率。
阳极极腔体模型同样基于质量守恒定律描述燃料电池阳极内部气体压力与外部操作条件耦合关系，其主要方程式如下：
其中，
和yvap,anin分别为经过加湿器后进入电池前氢气质量分数和水蒸气质量分数，变量
mvap,an分别为阳极流道内氢气和水蒸气质量；Wan,in和Wan,out分别为进入阳极和排出阳极的混合气质量流量。
膜水合模型主要描述膜内水传质过程，通常包含两个主要机制：水分子被质子拖拽引起的“电渗透”现象；阴阳极两侧水浓度差引起的“反扩散”现象，主要由以下等式描述：
其中，Acell为燃料电池的等效活性面积；Nv,mem为单位时间内穿透膜的水传递速率(mol.s-1.cm-2)，基于上述两种主要机制，假设水浓度梯度对于交换膜厚度方向呈线性分布，因此水通量可表述为：
其中，nd为电子拖拽系数；Dw为浓差扩散系数；cv,ca和cv,an分别为阴阳极两侧水浓度；Lm为膜的厚度；Dw是膜含水量λm和膜温度Tm的函数，nd同样是膜含水量函数
电化学模型主要描述电池稳态电压输出。燃料电池由于受到极化作用，输出电压需要从理想电动势中减去极化作用造成的过电势，主要包括活化过电势、欧姆过电势以及浓差过电势，主要表达式如下：
其中，Tst为燃料电池温度，
为阳极侧氢气分压，
为阴极侧氧气分压，R为气体常数，F为法拉第常数，Lm为膜的厚度，Acell为电池的等效活性面积；λm为交换膜的含水量，θ1～θ10为待识别参数，Ist为负载电流。
S2，利用遗传粒子群混合寻优算法对模型进行参数辨识
辨识算法过程如图3所示，主要过程如下；
(1)种群初始化：定义种群大小、空间维度、遗传寻优操作和粒子群寻优操作的基本参数，并在一定范围内初始化粒子速度和位置；
(2)计算每个粒子的个体最优和全局最优，并记录每个粒子的初始位置；
(3)对粒子进行粒子群寻优操作，更新各粒子的位置和速度；
(4)判断更新后的粒子适应度是否满足要求，若满足适应度要求，则停止寻优；
(5)如不满足要求，则对粒子进行选择、交叉、变异操作；将经过遗传学操作后的粒子再进行一次速度和位置更新，看粒子是否满足要求，若不满足则重复上述过程，直到搜索到问题最优解。
本发明以新标欧洲循环测试(NEDC)工况数据为实例，将其工况下测得的外部操作条件作为模型输入，电池电压为模型输出。将实测电压和模型输出电压的误差平方差和定义为适应度，表达式如下：
其中，N为采样点个数，Vcell为实测电压，
为模型估计电压。
利用混合遗传粒子群优化算法进行参数辨识，直至适应度值满足一定要求。此实例模型辨识结果如图4所示，结果表明该模型及参数辨识结果能够有效跟踪实测电压。
S3，将高频阶跃电流信号施加到燃料电池模型中，得到高频电压响应数据
获取燃料电池阻抗需要不同频率的正弦激励，传统的阻抗测试设备难以用于实车在线测量，采用激励源装置方案会额外增加燃料电池系统成本，因此需要寻找一种无需激励源的在线阻抗估计方法。阶跃电流信号可认为是由无穷不同幅值不同频率的正弦信号组成，可由控制器处理器生成，与正弦信号发生器和信号合成器产生的信号等效，将其施加到步骤S2辨识后的动态燃料电池集总参数模型，即可获取电压动态响应。本实例模拟施加的阶跃电流信号如图5所示，阶跃变化的电流的变载幅值可为0.1A/cm2～1.2A/cm2电密，阶跃变载可为加载或减载，仿真电流的采样频率按照设定的目标分析频率设定。
S4，利用快速傅里叶变化估算燃料电池阻抗
一旦获取电流和电压的时域数据，需要对其进行时域到频率的转换。为了减小泄露，需要对数据进行加窗处理。选择窗函数的原则是：窗谱的主瓣窄而高，这样它可以有陡峭的通频带，从而保证谱分析有足够的分辨率；旁瓣幅值应小，且旁瓣随频率的增大要衰减得尽量快，以减小频谱泄露。但通常上述两点难以同时满足。窗口函数所确定的分辨率与精度是相互矛盾的，在相同信号样本长度的情况下，精度的提高往往会导致分辨率降低。因此，实际应用中进行窗函数的选用要综合考虑分辨率和精度问题，应根据信号性质特点以及需求进行具体分析。本实例选择汉宁窗对阶跃电流和电压的时域数据进加窗处理。
采用快速傅里叶变换加窗后的电压和电流进行快速傅里叶变换，得到对应的频域数据，利用不同频率下的电压模除以不同频率下的电流模，即可获得对应频率下的阻抗，具体表达式如下：
本实例计算的阻抗效果如图6所示，仿真结果表明该方法能够估计燃料电池阻抗。该模型可嵌入于燃料电池控制器中，可将阶跃电流施加到该模型中，通过阶跃电压响应即可计算燃料电池阻抗，可实现车载在线应用。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。显然，本行业的技术人员应该了解本发明的上述实施例仅是为清楚地说明本发明所做的举例，而并非是对本发明的实施方式的限定。本领域内的技术人员可以在本发明主要思想内做其他变化，这些所引申出的显而易见的变化，都应包含在本发明所要求保护的范围内。
Claims (4)
1.一种基于模型的质子交换膜燃料电池阻抗在线估计方法，其特征在于，包括以下步骤：
S1建立适用于实车控制器应用的动态燃料电池集总参数模型，所述的动态燃料电池集总参数模型包括用以描述阴极内部气体压力与外部操作条件耦合关系的阴极腔体模型、用以描述阳极内部气体压力与外部操作条件耦合关系的阳极腔体模型、用以描述膜内水传质过程的膜水合模型以及用以描述电压响应的电压电化学模型，所述的外部操作条件包括氢气质量流量、空气质量流量、阴阳极入口温度、阴阳极入口压力、阴阳极进口相对湿度、阴阳极出口压力和负载电流；
所述的阴极腔体模型具体为：
其中，
mvap,ca和mliq,ca分别为阴极流道内氧气、氮气、水蒸气和液态水的质量，
yvap,in和
yvap,out分别为各气体进入阴极腔体和排出阴极腔体的质量分数，Wca,in和Wca,out分别为进入阴极腔体和排出阴极腔体混合气的质量流量，Wliq,ca为排出阴极腔体的液态水质量，
为阴极侧水的相变率，β为膜内净水传质速率，Ist为负载电流，F为法拉第常数，
为氧气的摩尔质量，
为水的摩尔质量；
所述的阳极腔体模型具体为：
其中，
和yvap,anin分别为经过加湿器后进入电池前氢气质量分数和水蒸气质量分数，
mvap,an分别为阳极流道内氢气和水蒸气质量，Wan,in和Wan,out分别为进入阳极和排出阳极的混合气质量流量，
为氢气的摩尔质量，
为水的摩尔质量，mliq,an为阳极侧液态水质量，Wliq,an为排出阳极腔体的液态水质量，
为排出阳极腔体的氢气质量分数，yvap,anout为排出阳极水蒸气的质量分数，mphase,an为相变水质量，Ist为负载电流，F为法拉第常数，β为膜内净水传质速率；
所述的膜水合模型具体为：
其中，Acell为燃料电池的等效活性面积，F为法拉第常数，Ist为负载电流，β为膜内净水传质速率，Nv,mem为单位时间内穿透膜的水传递速率，nd为电子拖拽系数，Dw为浓差扩散系数，cv,ca和cv,an分别为阴阳极两侧水浓度，Lm为膜的厚度；
所述的电化学模型具体为：
其中，Vcell为模型估计输出电池电压，Tst为燃料电池温度，
为阳极侧氢气分压，
为阴极侧氧气分压，R为气体常数，F为法拉第常数，Lm为膜的厚度，Acell为电池的等效活性面积，λm为交换膜的含水量，θ1～θ10为待识别参数，Ist为负载电流；
S2采用遗传粒子群混合寻优算法对动态燃料电池集总参数模型进行参数辨识；
S3将高频阶跃电流信号施加到辨识后的动态燃料电池集总参数模型中，得到高频电压响应数据；
S4对时域数据进行加汉宁窗处理，并采用快速傅里叶变换对燃料电池的电压和电流进行时域到频域转换，基于变换后的电压和电流计算各频率下的燃料电池阻抗，各频率下的燃料电池阻抗Z的具体表达式为：
其中，|FFT(Vcell(t))|为快速傅里叶变换后的模型估计输出电池电压模，|FFT(Ist(t))|为快速傅里叶变换后的高频阶跃电流模，FFT为快速傅里叶计算，Vcell(t)为模型输出电压，Ist(t)为输入电流。
2.根据权利要求1所述的一种基于模型的质子交换膜燃料电池阻抗在线估计方法，其特征在于，所述的步骤S1中，基于以下假设条件建立动态燃料电池集总参数模型，则有：
1)所有气体均视为理想气体；
2)内部气体组分的压力、温度和浓度均匀分布；
3)不考虑热力学动态过程，电池温度视为恒定不变。
3.根据权利要求1所述的一种基于模型的质子交换膜燃料电池阻抗在线估计方法，其特征在于，所述的步骤S2中，混合遗传粒子群算法具体包括以下步骤：
定义种群大小、空间维度、遗传寻优操作和粒子群寻优操作的基本参数；
计算每个粒子的个体最优和全局最优，并记录每个粒子的初始位置；
对粒子进行粒子群算法操作，判断粒子适应度是否满足要求；
若不满足适应度要求，则对粒子进行选择、交叉、变异的遗传学操作；
将粒子再进行一次速度和位置更新，判断粒子是否满足适应度要求，若不满足则重复上述过程，直到搜索到问题最优解。
4.根据权利要求3所述的一种基于模型的质子交换膜燃料电池阻抗在线估计方法，其特征在于，所述的粒子适应度定义为实测电压和模型输出电池电压的误差平方和。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010068829.0A CN111261903B (zh) | 2020-01-21 | 2020-01-21 | 一种基于模型的质子交换膜燃料电池阻抗在线估计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010068829.0A CN111261903B (zh) | 2020-01-21 | 2020-01-21 | 一种基于模型的质子交换膜燃料电池阻抗在线估计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111261903A CN111261903A (zh) | 2020-06-09 |
| CN111261903B true CN111261903B (zh) | 2021-03-26 |
Family
ID=70947048
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010068829.0A Active CN111261903B (zh) | 2020-01-21 | 2020-01-21 | 一种基于模型的质子交换膜燃料电池阻抗在线估计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111261903B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111736082B (zh) * | 2020-08-24 | 2020-11-13 | 杭州华塑加达网络科技有限公司 | 一种交流阻抗谱测试方法、装置及介质 |
| CN112652794B (zh) * | 2020-12-18 | 2022-11-11 | 福州大学 | 一种利用时滞信息的阴极开放式燃料电池热管理系统及方法 |
| CN113340788A (zh) * | 2021-06-10 | 2021-09-03 | 同济大学 | 一种燃料电池气渗定量检测方法 |
| CN114361535B (zh) * | 2021-12-16 | 2023-08-04 | 同济大学 | 一种基于电化学阻抗谱的燃料电池氢渗量评估方法 |
| CN114492087B (zh) * | 2022-04-02 | 2022-07-19 | 国网浙江省电力有限公司电力科学研究院 | 氢储能电站质子交换膜燃料电池故障诊断方法及装置 |
| CN116207306B (zh) * | 2023-04-24 | 2023-08-01 | 中汽数据有限公司 | 燃料电池电拖拽系数确定方法 |
| CN116500340B (zh) * | 2023-05-18 | 2024-02-02 | 浙江蓝能氢能科技股份有限公司 | 一种电解制氢装置的阻抗测量方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105449241A (zh) * | 2015-12-18 | 2016-03-30 | 武汉理工大学 | 燃料电池在线电化学阻抗谱检测系统及方法 |
| CN106450388A (zh) * | 2016-09-26 | 2017-02-22 | 中国计量大学 | 一种水冷型燃料电池温度优化与控制方法 |
| CN110161421A (zh) * | 2019-05-22 | 2019-08-23 | 同济大学 | 一种在线重构设定频率范围内电池阻抗的方法 |
| CN110688808A (zh) * | 2019-09-24 | 2020-01-14 | 南通大学 | 一种动力电池模型的粒子群及lm优化混合迭代辨识方法 |
-
2020
- 2020-01-21 CN CN202010068829.0A patent/CN111261903B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105449241A (zh) * | 2015-12-18 | 2016-03-30 | 武汉理工大学 | 燃料电池在线电化学阻抗谱检测系统及方法 |
| CN106450388A (zh) * | 2016-09-26 | 2017-02-22 | 中国计量大学 | 一种水冷型燃料电池温度优化与控制方法 |
| CN110161421A (zh) * | 2019-05-22 | 2019-08-23 | 同济大学 | 一种在线重构设定频率范围内电池阻抗的方法 |
| CN110688808A (zh) * | 2019-09-24 | 2020-01-14 | 南通大学 | 一种动力电池模型的粒子群及lm优化混合迭代辨识方法 |
Non-Patent Citations (1)
| Title |
|---|
| 《小型燃料电池空气系统动态模型的仿真与分析》;熊力 等;《2007中国汽车工程学会年会论文集》;20071231;第181-186页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111261903A (zh) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111261903B (zh) | 一种基于模型的质子交换膜燃料电池阻抗在线估计方法 | |
| Yuan et al. | Understanding dynamic behavior of proton exchange membrane fuel cell in the view of internal dynamics based on impedance | |
| Hu et al. | A novel diagnostic methodology for fuel cell stack health: performance, consistency and uniformity | |
| Gu et al. | Data-based flooding fault diagnosis of proton exchange membrane fuel cell systems using LSTM networks | |
| CN111199122B (zh) | 基于多物理场的质子交换膜燃料电池水淹故障诊断方法 | |
| Wahdame et al. | Comparison between two PEM fuel cell durability tests performed at constant current and under solicitations linked to transport mission profile | |
| Hu et al. | Carbon corrosion induced fuel cell accelerated degradation warning: From mechanism to diagnosis | |
| CA2940150C (en) | Use of neural network and eis signal analysis to quantify h2 crossover in-situ in operating pem cells | |
| Sheng et al. | Energy management strategy based on health state for a PEMFC/Lithium-ion batteries hybrid power system | |
| CN113608131B (zh) | 一种动态机车工况下pemfc电堆衰退性能检测方法 | |
| CN109841879A (zh) | 燃料电池水含量估计系统、方法、计算机设备及存储介质 | |
| CN111063919B (zh) | 燃料电池阳极状态监测方法 | |
| Ma et al. | Water content diagnosis for proton exchange membrane fuel cell based on wavelet transformation | |
| CN117214746A (zh) | 一种燃料电池的交流阻抗测量系统 | |
| Chadha et al. | Study of new flow field geometries to enhance water redistribution and pressure head losses reduction within PEM fuel cell | |
| CN114264881A (zh) | 一种燃料电池阻抗在线监控方法及系统 | |
| Yuan et al. | Advanced online broadband impedance spectrum acquisition of fuel cells by S-transform | |
| CN114388855B (zh) | 一种燃料电池寿命预测方法 | |
| Zhao et al. | Inconsistency evaluation of vehicle-oriented fuel cell stacks based on electrochemical impedance under dynamic operating conditions | |
| Tang et al. | A novel multi-step investigation of in-plane heterogeneity for commercial-size fuel cells based on current distribution model and multi-point impedance method | |
| CN103852222B (zh) | 阳极渗漏位置检测 | |
| Choi et al. | Real-time detection of flooding in polymer electrolyte membrane fuel cells using high-frequency electrochemical impedance | |
| CN116231013A (zh) | 燃料电池故障嵌入模型的建立方法 | |
| CN111007404A (zh) | 一种基于关键频率点的燃料电池阻抗测量分析系统及方法 | |
| CN112993341A (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211215 Address after: 201821 room jt19773, building 4, zone B, No. 925, Yecheng Road, Jiading Industrial Zone, Jiading District, Shanghai Patentee after: Haizhuo power (Shanghai) Energy Technology Co.,Ltd. Address before: 200092 Siping Road 1239, Shanghai, Yangpu District Patentee before: TONGJI University |