JP6277864B2 - Battery internal state estimation device - Google Patents
Battery internal state estimation device Download PDFInfo
- Publication number
- JP6277864B2 JP6277864B2 JP2014107775A JP2014107775A JP6277864B2 JP 6277864 B2 JP6277864 B2 JP 6277864B2 JP 2014107775 A JP2014107775 A JP 2014107775A JP 2014107775 A JP2014107775 A JP 2014107775A JP 6277864 B2 JP6277864 B2 JP 6277864B2
- Authority
- JP
- Japan
- Prior art keywords
- wavelet transform
- voltage
- internal impedance
- secondary battery
- current
- 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
- 238000004364 calculation method Methods 0.000 claims description 79
- 238000001514 detection method Methods 0.000 claims description 55
- 230000008859 change Effects 0.000 claims description 45
- 238000009792 diffusion process Methods 0.000 claims description 35
- 238000006243 chemical reaction Methods 0.000 claims description 30
- 238000010586 diagram Methods 0.000 claims description 27
- 230000006866 deterioration Effects 0.000 claims description 22
- 238000007599 discharging Methods 0.000 claims description 22
- 238000012937 correction Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 9
- 238000002485 combustion reaction Methods 0.000 claims description 7
- 230000000737 periodic effect Effects 0.000 claims description 5
- 239000002131 composite material Substances 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims 1
- 238000000034 method Methods 0.000 description 29
- 238000012545 processing Methods 0.000 description 17
- 230000006870 function Effects 0.000 description 13
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 7
- 229910001416 lithium ion Inorganic materials 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000000691 measurement method Methods 0.000 description 4
- 230000010287 polarization Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 235000009413 Ratibida columnifera Nutrition 0.000 description 1
- 241000510442 Ratibida peduncularis Species 0.000 description 1
- 239000011149 active material Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 229910052987 metal hydride Inorganic materials 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- PXHVJJICTQNCMI-UHFFFAOYSA-N nickel Substances [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 1
- -1 nickel metal hydride Chemical class 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- 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/10—Energy storage using batteries
Landscapes
- Battery Mounting, Suspending (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Tests Of Electric Status Of Batteries (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Secondary Cells (AREA)
Description
二次電池の内部インピーダンスを算出する電池内部状態推定装置に関する。 The present invention relates to a battery internal state estimation device that calculates internal impedance of a secondary battery.
二次電池を安全に運用するためには、二次電池の内部状態を把握する必要がある。二次電池の内部状態を示す指標として、二次電池の充電量（ＳＯＣ：Ｓｔａｔｅ ｏｆ Ｃｈａｒｇｅ）、充放電許容電力、内部抵抗、及び、劣化度などがある。 In order to operate the secondary battery safely, it is necessary to grasp the internal state of the secondary battery. As an index indicating the internal state of the secondary battery, there are a charge amount (SOC: State of Charge) of the secondary battery, allowable charge / discharge power, internal resistance, a degree of deterioration, and the like.
二次電池の内部状態を把握するためには、二次電池の内部で発生する電気的・化学的現象を捉える方法が必要となる。電圧と電流との応答は、電気的・化学的現象により変化する。電圧と電流との応答は二次電池の内部インピーダンスにあたる。この二次電池の内部インピーダンスに基づき、二次電池を等価回路でモデル化し、その等価回路を二次電池の特性計測や制御に用いる方法が知られている。 In order to grasp the internal state of the secondary battery, a method for capturing electrical and chemical phenomena occurring inside the secondary battery is required. The response between voltage and current varies with electrical and chemical phenomena. The response between voltage and current corresponds to the internal impedance of the secondary battery. A method is known in which a secondary battery is modeled by an equivalent circuit based on the internal impedance of the secondary battery, and the equivalent circuit is used for characteristic measurement and control of the secondary battery.
二次電池の内部インピーダンスは、複数の抵抗成分と複数の容量成分の組み合わせでモデル化できることが知られている。ある周波数における内部インピーダンスの実部及び虚部の絶対値が、その周波数における二次電池の内部抵抗に相当する。二次電池に入出力される電圧及び電流の周波数に応じて内部インピーダンスは変化するため、二次電池の内部抵抗は変化することになる。 It is known that the internal impedance of a secondary battery can be modeled by a combination of a plurality of resistance components and a plurality of capacity components. The absolute values of the real part and the imaginary part of the internal impedance at a certain frequency correspond to the internal resistance of the secondary battery at that frequency. Since the internal impedance changes according to the frequency of the voltage and current input / output to / from the secondary battery, the internal resistance of the secondary battery changes.
例えば、車両に搭載された二次電池について考える。車両搭載環境において、二次電池には、発電機及びモータなどの電気負荷が直流−交流変換装置を介して接続され、二次電池と発電機及び電気負荷との間で非周期的な充放電が行われる。この非周期的な充放電において、充放電電流の大きさ、及び、充放電の継続時間が変化するため、正弦波を印加してインピーダンスを計測する交流インピーダンス法は適用できない。また、刻々と二次電池の内部状態が変化する中で周波数特性を精度よく算出する解析手法が必要となる。 For example, consider a secondary battery mounted on a vehicle. In a vehicle-mounted environment, an electrical load such as a generator and a motor is connected to the secondary battery via a DC-AC converter, and aperiodic charge / discharge is performed between the secondary battery, the generator and the electrical load. Is done. In this non-periodic charge / discharge, since the magnitude of the charge / discharge current and the duration of charge / discharge change, the AC impedance method of measuring impedance by applying a sine wave cannot be applied. In addition, an analysis method for accurately calculating the frequency characteristics is required while the internal state of the secondary battery changes every moment.
車両搭載環境において得られる二次電池の電圧と電流から、二次電池の内部抵抗の周波数特性を算出する方法が提案されている。二次電池の充放電電流と端子間電圧についてウェーブレット変換を適用して、得られた電圧ウェーブレット変換係数と電流ウェーブレット変換係数の傾きから内部インピーダンスを算出し、内部インピーダンスの実部の大きさを内部抵抗とし、その周波数特性を算出する方法が提案されている（特許文献１）。 There has been proposed a method for calculating the frequency characteristics of the internal resistance of the secondary battery from the voltage and current of the secondary battery obtained in a vehicle-mounted environment. Applying wavelet transform to the charge / discharge current of the secondary battery and the voltage between terminals, the internal impedance is calculated from the slope of the obtained voltage wavelet transform coefficient and current wavelet transform coefficient, and the size of the real part of the internal impedance is A method of calculating the frequency characteristic of the resistor is proposed (Patent Document 1).
従来技術の方法では、車両が運転状態下である場合など、動作中のシステムにおいて得られる二次電池の電圧及び電流の波形から、二次電池のある周波数の内部抵抗を求めることができるとされている。しかしながら、従来技術の方法では、内部インピーダンスの実部の大きさを算出しているだけである。周波数情報と共に、内部インピーダンスの実部と虚部の両方あるいは実部と虚部の絶対値を求めないと、二次電池の各種現象の寄与度として、実部に関する複数の抵抗成分及び虚部に関する複数の容量成分を求めることができない。虚部がなければ、各種現象に対応する複数の抵抗成分を切り分けができない。また、虚部を含まないインピーダンス絶対値の周波数特性からは容量成分を求められない。以上より、各種現象の時定数が求まらない。つまり、動作中のシステムにおいて、二次電池の各種現象の寄与度を求めることができず、二次電池の内部状態を精度よく推定できない。 In the prior art method, the internal resistance at a certain frequency of the secondary battery can be obtained from the voltage and current waveforms of the secondary battery obtained in the operating system, such as when the vehicle is operating. ing. However, the prior art method only calculates the size of the real part of the internal impedance. If the absolute value of both the real part and the imaginary part of the internal impedance or the real part and the imaginary part is not obtained together with the frequency information, the contribution of various phenomena of the secondary battery is related to a plurality of resistance components and the imaginary part related to the real part. A plurality of capacitive components cannot be obtained. Without an imaginary part, a plurality of resistance components corresponding to various phenomena cannot be separated. Further, the capacitance component cannot be obtained from the frequency characteristic of the impedance absolute value not including the imaginary part. From the above, time constants of various phenomena cannot be obtained. That is, in the operating system, the contribution degree of various phenomena of the secondary battery cannot be obtained, and the internal state of the secondary battery cannot be accurately estimated.
本発明は、上記課題を解決するためになされたものであり、動作中のシステムにおいても、二次電池の内部状態の詳細を精度よく推定できる電池内部状態推定装置を提供することを目的とする。 The present invention has been made to solve the above-described problem, and an object of the present invention is to provide a battery internal state estimation device that can accurately estimate the details of the internal state of a secondary battery even in an operating system. .
本発明における電池内部状態推定装置は、二次電池（１０）の端子間電圧を検出する電圧検出手段（３１）及び前記二次電池に流れる充放電電流を検出する電流検出手段（３０）による検出値を時系列で記憶する記憶手段（５４）と、前記時系列で記憶された端子間電圧の検出値に対してウェーブレット変換を行うことで電圧ウェーブレット変換係数を算出するとともに、前記時系列で記憶された充放電電流の検出値に対してウェーブレット変換を行うことで電流ウェーブレット変換係数を算出するウェーブレット変換手段（６０）と、前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、周波数ごとの前記二次電池の内部インピーダンスの実部及び虚部を算出する内部インピーダンス算出手段（６１）と、を備えることを特徴とする。 The battery internal state estimation device according to the present invention is detected by voltage detection means (31) for detecting the voltage across the terminals of the secondary battery (10) and current detection means (30) for detecting the charge / discharge current flowing in the secondary battery. A storage means (54) for storing values in time series, and a voltage wavelet transform coefficient is calculated by performing wavelet transform on the detected value of the inter-terminal voltage stored in the time series, and is stored in the time series. Wavelet transform means (60) for calculating a current wavelet transform coefficient by performing wavelet transform on the detected charge / discharge current value, and a ratio between the voltage wavelet transform coefficient and the current wavelet transform coefficient for each frequency. Internal impedance calculation means (61) for calculating a real part and an imaginary part of the internal impedance of the secondary battery Characterized in that it obtain.
ウェーブレット変換は、フーリエ変換とは異なり、時間情報を失わず周波数情報と両方が得られる特徴があり、実部と虚部の両方の精度が期待できる。また、ウェーブレット変換は、変換を行う対象が非周期的な波形でも、各周波数成分を精度よく算出することができる。つまり、非周期的で複雑な充放電を行いながら、二次電池の内部状態として、各周波数における二次電池の内部インピーダンスの実部及び虚部を精度よく算出することができる。また、ウェーブレット変換を用いると、二次電池の充放電が実施され、その後に充放電が停止された際の緩やかな電圧変化、つまり、周期１つに対する断片的な波形からでも内部インピーダンスの実部及び虚部を算出できる。充放電を停止している際は、充電量変化がなく、二次電池内部の発熱もないため、低周波数域における内部インピーダンスの実部及び虚部を精度よく算出できる。このように、動作中のシステムにおいても、二次電池の内部状態の詳細を精度よく推定することができる。 Unlike the Fourier transform, the wavelet transform is characterized in that both frequency information can be obtained without losing time information, and the accuracy of both the real part and the imaginary part can be expected. In addition, the wavelet transform can accurately calculate each frequency component even if the object to be transformed is an aperiodic waveform. That is, the real part and the imaginary part of the internal impedance of the secondary battery at each frequency can be accurately calculated as the internal state of the secondary battery while performing aperiodic and complicated charge / discharge. In addition, when the wavelet transform is used, charging / discharging of the secondary battery is carried out, and then the gradual voltage change when charging / discharging is stopped, that is, the real part of the internal impedance even from a fragmentary waveform for one period. And the imaginary part can be calculated. When charging / discharging is stopped, there is no change in the charge amount, and there is no heat generation in the secondary battery, so that the real part and imaginary part of the internal impedance in the low frequency range can be calculated with high accuracy. Thus, even in the operating system, the details of the internal state of the secondary battery can be accurately estimated.
（第１実施形態）
本実施形態における組電池Ｂ及び電池内部状態推定装置としての制御装置５０を備える電気回路図を図１に示す。組電池Ｂは、リチウムイオン蓄電池からなる複数個の蓄電池１０（電池セル）が直列接続または並列接続された組電池であり、制御装置５０は、組電池Ｂを構成する蓄電池１０（電池セル）の内部状態をそれぞれ推定する。組電池Ｂ及び制御装置５０は、例えば、車両に搭載されているものとする。その車両は、内燃機関とモータジェネレータとを備えるハイブリッド車であるとする。
(First embodiment)
FIG. 1 shows an electric circuit diagram including the assembled battery B and the control device 50 as a battery internal state estimating device in the present embodiment. The assembled battery B is an assembled battery in which a plurality of storage batteries 10 (battery cells) made of lithium ion storage batteries are connected in series or in parallel, and the control device 50 includes the storage battery 10 (battery cells) constituting the assembled battery B. Each internal state is estimated. Assume that the assembled battery B and the control device 50 are mounted on a vehicle, for example. The vehicle is assumed to be a hybrid vehicle including an internal combustion engine and a motor generator.
蓄電池１０は、交流−直流変換を行うインバータ２１を介してモータジェネレータ２２に接続されている。モータジェネレータ２２が発電機として機能する場合、蓄電池１０はモータジェネレータ２２から供給される電力によって充電される。また、モータジェネレータ２２が動力源としてのモータとして機能する場合、蓄電池１０から放電することでモータジェネレータ２２に対して電力供給を行う。車両は、内燃機関での燃焼を停止させ駆動力の出力を停止し、モータジェネレータ２２をモータとして機能させるＥＶ（Ｅｌｅｃｔｒｉｃ Ｖｅｈｉｃｌｅ）走行モードを備える。また、内燃機関での燃焼を実施して駆動力を出力し、モータジェネレータ２２をモータ又は発電機として機能させるＨＶ（Ｈｙｂｒｉｄ Ｖｅｈｉｃｌｅ）走行モードを備える。
The
蓄電池１０は、電気負荷２０に接続され、電気負荷２０に対し電力を供給する。また、商用電源などから蓄電池１０に対して充電する事が可能である。また、蓄電池１０と電気負荷２０及びインバータ２１、商用電源との間にはリレースイッチ３３，３４が設置されていて、制御装置５０の指令により、蓄電池１０と主経路との接続及び切断を行う。蓄電池１０と電気負荷２０及びインバータ２１とを接続する経路上には蓄電池１０に流れる充放電電流を検出するための電流センサ３０が設けられている。蓄電池１０の端子間には、その端子間の電圧を検出するための電圧センサ３１が設けられている。また、蓄電池１０には温度センサ３２が設けられている。これら電流センサ３０（電流検出手段）、電圧センサ３１（電圧検出手段）及び温度センサ３２はそれぞれ充放電電流、端子間電圧及び電池温度に応じた検出信号を出力し、その検出信号は制御装置５０に入力される。
The
制御装置５０は、入力される検出信号に基づいて、電流センサ３０の検出値Ｉ、電圧センサ３１の検出値Ｖ及び温度センサ３２の検出値Ｔを取得する。制御装置５０は、これら取得した検出値Ｉ，Ｖ，Ｔに基づいて、蓄電池１０のＳＯＣ、劣化度、予測電力を算出する。
The control device 50 acquires the detection value I of the
続いて、制御装置５０で行われる本実施形態の演算方法について以下述べていく。図２に本実施形態の制御装置５０の機能を表す機能ブロック図を示す。制御装置５０は、内部インピーダンス処理部５１、ＲＣ成分算出部６２、ＳＯＣ処理部５２、劣化度処理部５３、及び、予測電力処理部５４から構成される。
Subsequently, the calculation method of the present embodiment performed by the control device 50 will be described below. FIG. 2 is a functional block diagram showing functions of the control device 50 of the present embodiment. The control device 50 includes an internal
まず、内部インピーダンス処理部５１について述べる。制御装置５０は、電流センサ３０及び電圧センサ３１によって検出された蓄電池１０の端子間電圧の検出値Ｖと、充放電電流の検出値Ｉを時系列で記憶する。そして、時系列で記憶された端子間電圧の検出値Ｖ（ｔ）と、充放電電流の検出値Ｉ（ｔ）をウェーブレット変換部６０でそれぞれウェーブレット変換する。
First, the internal
本発明で利用する技術の原理について述べる。波形ｆ（ｔ）の連続ウェーブレット変換は式（１）で定義される。 The principle of the technology used in the present invention will be described. The continuous wavelet transform of the waveform f (t) is defined by equation (1).
二次電池が車両に搭載され走行している場合の電圧と電流の波形に対して種々のマザーウェーブレットを適用した。その結果、モルレー関数を適用した場合と、ガウス関数を適用した場合とを比較すると、モルレー関数を適用した方が、得られたスペクトルの幅が狭く、二次電池の内部インピーダンスを高精度に算出できる。 Various mother wavelets were applied to the voltage and current waveforms when the secondary battery was mounted on a vehicle. As a result, comparing the case where the Morray function is applied and the case where the Gaussian function is applied, the obtained spectrum is narrower and the internal impedance of the secondary battery is calculated with higher accuracy when the Morray function is applied. it can.
離散ウェーブレット変換として、ハールのスケーリング関数による分解で得る方法がある。この方法では、式（２）について、ａ＝１／２＾ｊ，ｂ＝ｋ／２＾ｊの関係となる。
式（３）に、端子間電圧Ｖのウェーブレット変換係数ＷψＶと充放電電流Ｉのウェーブレット変換係数ＷψＩを示す。
As a discrete wavelet transform, there is a method obtained by decomposition using a Haar scaling function. In this method, the relationship of a = 1/2 ^ j and b = k / 2 ^ j is established for the expression (2).
Equation (3) shows the wavelet transform coefficient WψV of the inter-terminal voltage V and the wavelet transform coefficient WψI of the charge / discharge current I.
以上のように、ウェーブレット変換部６０によるウェーブレット変換により、蓄電池１０の周波数ｆ毎の内部インピーダンスＺの実部Ｒｅ（Ｚ（ｆ））及び虚部Ｉｍ（Ｚ（ｆ））を得ることができる。
As described above, the real part Re (Z (f)) and the imaginary part Im (Z (f)) of the internal impedance Z for each frequency f of the
次に、電圧補正部５６について述べる。図３に、蓄電池１０が持つ特性として、充電量ＳＯＣと開放端電圧ＯＣＶとの関係を表す概念図を示す。蓄電池１０が充放電されている場合の開放端電圧ＯＣＶの変化について着目する。蓄電池１０の特性として、充電量ＳＯＣと開放端電圧ＯＣＶとは一対一の関係があり、一般的には制御性を考慮して直線的な領域を使用し、充電量ＳＯＣの算出などに用いられる。蓄電池１０が充放電されると、充電量ＳＯＣが変化し、その変化に伴って開放端電圧ＯＣＶも変化する。
Next, the
仮に、ＳＯＣの変化による蓄電池１０の開放端電圧ＯＣＶの変化を考慮することなく、内部インピーダンスＺを算出した場合は、開放端電圧ＯＣＶの変化が内部インピーダンスＺの容量成分の一部に相当する振舞いとなる。つまり、蓄電池の内部インピーダンスとは異なる容量成分が誤って加算され、算出する内部インピーダンスＺの精度が低下する。そのため、内部インピーダンスＺを算出する際は、検出した端子間電圧の波形に含まれる開放端電圧ＯＣＶの変化を除いて、端子間電圧の基準を揃えるなどの処置が必要になる。
If the internal impedance Z is calculated without considering the change in the open circuit voltage OCV of the
そこで、本実施形態では、現在のＳＯＣにおけるＯＣＶをＯＣＶ算出部５７により算出し、ＳＯＣの変化による開放端電圧ＯＣＶの変化量ΔＯＣＶをΔＯＣＶ算出部５８により算出する（ΔＯＣＶ＝ＯＣＶｉｎｉ−ＯＣＶ）。なお、ＯＣＶは現在値、ＯＣＶｉｎｉは基準値であり、ＯＣＶｉｎｉは記憶部５５でＶａの記憶を開始する時のＯＣＶとする。そして、端子間電圧の検出値ＶからΔＯＣＶを減算部５９により減算することで補正した値Ｖａを記憶部５５により記憶する（Ｖａ＝Ｖ−ΔＯＣＶ）。そして、記憶部５５により時系列で記憶された補正後の電圧検出値Ｖａ（ｔ）をウェーブレット変換し、そのウェーブレット変換係数ＷψＶａを用いて内部インピーダンスＺの実部Ｒｅ（Ｚ）及び虚部Ｉｍ（Ｚ）を算出する構成とした。
Therefore, in the present embodiment, the OCV in the current SOC is calculated by the
次に、ＲＣ成分算出部６２について述べる。ＲＣ成分算出部６２は、インピーダンス算出部６１が算出した周波数ｆごとの内部インピーダンスＺの実部Ｒｅ（Ｚ（ｆ））及び虚部Ｉｍ（Ｚ（ｆ））に基づいて、内部インピーダンスＺの等価回路モデルにおける抵抗成分及び容量成分を算出する。
Next, the RC
蓄電池１０の等価回路を図４に示す。蓄電池１０の等価回路は、開放端電圧ＯＣＶ１１と内部インピーダンスＺのモデル１２とを直列に接続した構成とする。内部インピーダンスＺのモデル１２は、溶液中や電極の通電抵抗を表す直流抵抗１３、正極及び負極における電極界面反応を表す反応抵抗１４，１５、及び、活物質中や溶液中のイオン拡散を表す拡散抵抗１６を直列に接続した構成とする。
An equivalent circuit of the
拡散抵抗１６は、抵抗成分及び容量成分の並列接続体が複数個直列に接続した構成とする。抵抗成分と容量成分の並列接続による等価回路はフォスタ型等価回路と呼ばれる。この他に、拡散抵抗１６をカウエル型等価回路で表現する方法もある。尚、カウエル型等価回路表現とフォスタ型等価回路表現は、抵抗成分と容量成分の値は異なるが変換が可能であり、等価である。
The diffused
直流抵抗１３は、抵抗成分Ｒｓとして構成されている。反応抵抗１４は、抵抗成分Ｒ１及び容量成分Ｃ１が並列接続されて構成されている。反応抵抗１５は、抵抗成分Ｒ２及び容量成分Ｃ２が並列接続されて構成されている。
The
拡散抵抗１６は、複数の抵抗成分Ｒｗｉ及び複数の容量成分Ｃｗｉがそれぞれ並列接続されて構成されている。ただし、ｉは拡散抵抗１６を構成する抵抗成分と容量成分の並列接続体の個数であり、ｉ＝０，１，…，ｎの範囲をとる。
The diffused
直流抵抗１３、反応抵抗１４，１５、及び、拡散抵抗１６の時定数はそれぞれ異なるものとなっている。直流抵抗１３の時定数は０．００１秒であり、略０として扱うことができる。また、反応抵抗１４，１５の時定数は約０．１秒〜約１秒である。拡散抵抗１６の時定数は数十秒から数千秒のオーダーである。
The DC resistors 13,
民生電池（リチウムイオン二次電池）での劣化耐久試験において、直流抵抗１３、時定数の小さい順に反応抵抗１４，１５、及び拡散抵抗１６で分類した場合での各抵抗成及び各容量成分の推移について説明する。各抵抗成及び各容量成分の値は、試験前の値に比べて、劣化によりいずれも増加する。また、それぞれ異なった増加率（劣化度）を示す。このように、各抵抗成分及び各容量成分の劣化度が異なるので、劣化度を好適に見積もるためには各抵抗成分及び各容量成分を個別に算出する必要がある。特に、反応抵抗１５の抵抗成分Ｒ２及び容量成分Ｃ２は、他の成分に比べて、顕著に増加する結果がある。
Transition of each resistance component and each capacitance component when categorized by
ここで、車両での蓄電池１０の電流・電圧の周波数範囲について述べる。車両（車載電源システム）の動作中において、蓄電池１０には非周期的な充放電が行われる。この車両動作中における蓄電池１０の電流・電圧の変化の高周波数側は、インバータ２１のキャリア周波数が数ｋＨｚと、一般的な蓄電池の直流抵抗１３及び電極界面の反応抵抗１４，１５が取得できる周波数範囲にある。一方で、低周波数側は、拡散抵抗１６の周波数域について、反応抵抗１５近傍の低周波数域は取得できる。しかし、更に低い周波数域であり拡散抵抗１６の終点に位置する超低周波数域は、車両走行中には含まれず、超低周波数域の内部インピーダンスが取得できない場合がある。
Here, the frequency range of the current / voltage of the
また、拡散抵抗１６は複数の抵抗成分と容量成分の並列回路が重なり合う合成インピーダンスとなる。従って、１組の抵抗成分と容量成分の並列回路ずつに分離することはできず、あくまで、合成インピーダンスに合うように合わせ込む必要がある。以上を考慮して、以下の実施形態とした。
The diffused
ＲＣ成分算出部６２は、以下に説明するＲＣ成分算出方法１及びＲＣ成分算出方法２を用いて抵抗成分及び容量成分を算出する。ＲＣ成分算出方法１では、直流抵抗１３、電極界面の反応抵抗１４，１５を算出する。ＲＣ成分算出方法２では、拡散抵抗１６を算出する。
The RC
まず、ＲＣ成分算出方法１について述べる。ベクトル軌跡線図を用いた方法を適用する。ベクトル軌跡線図は、ボード線図に比べ、複雑な等価回路モデルにむいている。
First, the RC
図５に、内部インピーダンスＺの周波数特性を表すベクトル軌跡線図の概念図を示す。蓄電池１０の端子間電圧の変化の周波数ｆが低いほど、内部インピーダンスＺの実部Ｒｅ（Ｚ）は大きくなる。また、便宜上、虚部Ｉｍ（Ｚ）は負値を正側にとっている。図５では、説明のために周波数ｆに応じて領域Ｄ０〜Ｄ２に分割している。周波数ｆと時定数τの関係は、２πｆ＝１／τ＝１／（Ｒ・Ｃ）となる。周波数ｆがｆ≫１／（２π・τ１）の領域を領域Ｄ０とし、周波数ｆがｆａ１＝１／（２π・τ１）近傍の領域を領域Ｄ１とし、周波数ｆがｆａ２＝１／（２π・τ２）近傍の領域をＤ２とする。
FIG. 5 shows a conceptual diagram of a vector locus diagram representing the frequency characteristics of the internal impedance Z. The real part Re (Z) of the internal impedance Z increases as the frequency f of the change in the inter-terminal voltage of the
周波数ｆがｆ≫ｆａ１となる領域Ｄ０では、直流抵抗１３が内部インピーダンスＺとして寄与する。つまり、内部インピーダンスＺの実部は、Ｒｓとなる。また、周波数ｆがｆｂ１（ｆｂ１＞ｆａ１）となる領域Ｄ０と領域Ｄ１との境界の最下点Ｐ１（Ｒｅ（Ｚ）＝Ｒｓ，Ｉｍ（Ｚ）≒０）より周波数ｆが低下すると、内部インピーダンスＺの実部及び虚部が増加する。
In the region D0 where the frequency f is f >> fa1, the
周波数ｆがｆａ１近傍となる領域Ｄ１におけるベクトル軌跡線図では、円弧が形成されている。領域Ｄ１では、反応抵抗１４が内部インピーダンスＺとして寄与する。周波数ｆがｆ＝ｆａ１となると、領域Ｄ１における円弧状のベクトル軌跡の頂点Ｑ１となる。頂点Ｑ１において、Ｒｅ（Ｚ）＝Ｒｓ＋Ｒ１／２，Ｉｍ（Ｚ）＝Ｒ１／２である。また、周波数ｆがｆｂ２（ｆａ１＞ｆｂ２＞ｆａ２）となると、領域Ｄ１と領域Ｄ２との境界の最下点Ｐ２となる。最下点Ｐ２において、Ｒｅ（Ｚ）＝Ｒｓ＋Ｒ１，Ｉｍ（Ｚ）≒０である。
In the vector locus diagram in the region D1 where the frequency f is near fa1, an arc is formed. In the region D1, the
周波数ｆがｆａ２近傍となる領域Ｄ２におけるベクトル軌跡線図では、領域Ｄ２と同様に円弧が形成されている。領域Ｄ２では、反応抵抗１５が内部インピーダンスＺとして寄与する。周波数ｆがｆ＝ｆａ２となると、領域Ｄ２における円弧状のベクトル軌跡の頂点Ｑ２となる。頂点Ｑ２において、Ｒｅ（Ｚ）＝Ｒｓ＋Ｒ１＋Ｒ２／２，Ｉｍ（Ｚ）＝Ｒ２／２である。また、周波数ｆがｆｂ３（ｆａ２＞ｆｂ３＞ｆａ３）となると、最下点Ｐ３となる。最下点Ｐ３において、Ｒｅ（Ｚ）＝Ｒｓ＋Ｒ１＋Ｒ２，Ｉｍ（Ｚ）≒０である。
In the vector locus diagram in the region D2 where the frequency f is in the vicinity of fa2, an arc is formed as in the region D2. In the region D2, the
このように探索された最下点Ｐ１〜Ｐ３のそれぞれの実部に基づき、各抵抗成分Ｒｓ，Ｒ１，Ｒ２を算出することができる。具体的には、最下点Ｐ１の実部が抵抗成分Ｒｓに相当する。また、最下点Ｐ２の実部から最下点Ｐ１の実部を減算したものが抵抗成分Ｒ１に相当する。同様に、最下点Ｐ３の実部から最下点Ｐ２の実部を減算したものが抵抗成分Ｒ２に相当する。 Each resistance component Rs, R1, R2 can be calculated based on the respective real parts of the lowest points P1 to P3 searched in this way. Specifically, the real part of the lowest point P1 corresponds to the resistance component Rs. Further, a value obtained by subtracting the real part of the lowest point P1 from the real part of the lowest point P2 corresponds to the resistance component R1. Similarly, a value obtained by subtracting the real part of the lowest point P2 from the real part of the lowest point P3 corresponds to the resistance component R2.
また、探索された頂点Ｑ１，Ｑ２に基づき、容量成分Ｃ１，Ｃ２を算出することができる。具体的には、頂点Ｑ１の周波数ｆａ１は、１／（２π・Ｃ１・Ｒ１）と等しいため、周波数ｆａ１及び抵抗成分Ｒ１に基づいて、容量成分Ｃ１を算出できる（Ｃ１＝１／（２π・Ｒ１・ｆａ１））。同様に、頂点Ｑ２の周波数ｆａ２及び抵抗成分Ｒ２に基づいて、容量成分Ｃ２を算出できる（Ｃ２＝１／（２π・Ｒ２・ｆａ２））。容量成分Ｃ１，Ｃ２の算出において、上記最下点Ｐ１〜Ｐ３の実部に基づき算出された抵抗成分Ｒ１，Ｒ２を用いる構成とする。なお、抵抗成分Ｒ１，Ｒ２として予め定めた所定値を用いる構成としてもよい。 Further, the capacity components C1 and C2 can be calculated based on the searched vertices Q1 and Q2. Specifically, since the frequency fa1 of the vertex Q1 is equal to 1 / (2π · C1 · R1), the capacitance component C1 can be calculated based on the frequency fa1 and the resistance component R1 (C1 = 1 / (2π · R1). -Fa1)). Similarly, the capacitance component C2 can be calculated based on the frequency fa2 and the resistance component R2 of the vertex Q2 (C2 = 1 / (2π · R2 · fa2)). In the calculation of the capacitance components C1 and C2, the resistance components R1 and R2 calculated based on the real parts of the lowest points P1 to P3 are used. Note that a predetermined value may be used as the resistance components R1 and R2.
ここで、領域Ｄ２より低周波数側におけるベクトル軌跡線図では、拡散抵抗１６は複数の抵抗成分と容量成分の並列回路が重なり合う合成インピーダンスとなる。拡散抵抗については、円弧の最下点と頂点が判別できない場合が多いので、次に示すＲＣ成分算出方法２を用いて算出する。
Here, in the vector locus diagram on the lower frequency side than the region D2, the diffused
ＲＣ成分算出方法２について述べる。拡散抵抗１６の抵抗成分Ｒｗｉ及び容量成分Ｃｗｉを算出する方法を以下に示す。ただし、ｉは拡散抵抗１６を構成する抵抗成分と容量成分の並列接続体の個数であり、ｉ＝０，１，…，ｎの範囲をとる。例えば、ｉ＝４とする。拡散抵抗は時定数が大きく、拡散抵抗を精度よく算出するためには、その時定数に応じた長時間にわたって検出された検出値Ｉ，Ｖを用いる必要がある。長時間にわたって検出された検出値Ｉ，Ｖを用いると、その検出期間における充放電による開放端電圧ＯＣＶ変化や蓄電池の温度上昇が精度に大きく影響を与える。そこで、電池状態変化が少ない状況で端子間電圧Ｖや充放電電流Ｉを測定する必要がある。
The RC
本実施形態の内部インピーダンス処理部５１は、車両が動作状態から停止状態に移行する車両動作終了時点、及び、車両停止中における検出値Ｉ，Ｖを取得し、その検出値に基づいて内部インピーダンスの実部Ｒｅ（Ｚ）及び虚部Ｉｍ（Ｚ）を算出する。そして、その内部インピーダンスの実部Ｒｅ（Ｚ）及び虚部Ｉｍ（Ｚ）を用いて、ＲＣ成分算出部６２は、拡散抵抗１６の抵抗成分Ｒｗｉ及び容量成分Ｃｗｉを算出する。
The internal
具体的には、図６に示すように、蓄電池１０の充放電状態から充放電停止状態に切り替わる少し前から、充放電停止状態の電圧Ｖや電流Ｉを取得する。ここで、蓄電池１０の充放電状態は、車両の動作・停止状態に対応する。つまり、車両の動作状態において、蓄電池１０は充放電状態となり、車両の停止状態において、蓄電池１０は充放電停止状態となる。車両の停止状態において、蓄電池１０の端子間電圧Ｖは緩変化し、その電圧Ｖの緩変化が拡散抵抗１６による変化に相当する。そこで、電圧Ｖの緩変化する期間における検出値Ｉ，Ｖを用いることで拡散抵抗１６を精度よく算出することが可能になる。
Specifically, as shown in FIG. 6, the voltage V and the current I in the charge / discharge stop state are acquired slightly before the
なお、内部インピーダンスを算出する際に、分母の電流Ｉを０にしないために、充放電停止状態の切り替わる少し前の電流及び電圧の検出値Ｉ，Ｖを取得する。この時の充放電による開放端電圧変化については前述同様に電圧Ｖの補正を行う。また、拡散抵抗１６は時定数が数十秒以上と、直流抵抗１３及び反応抵抗１４，１５の時定数に比べて大きいため、検出値のサンプリング周期を直流抵抗１３及び反応抵抗１４，１５の算出時に比べて大きくしてもよい。
When calculating the internal impedance, the current and voltage detection values I and V just before the charge / discharge stop state is switched are acquired so that the denominator current I is not zero. As for the open-circuit voltage change due to charging / discharging at this time, the voltage V is corrected as described above. Further, since the
次に、拡散抵抗の抵抗成分と容量成分の算出方法を述べる。拡散抵抗１６は合成インピーダンスであるため、ＲＣ成分算出方法１で扱う最下点が虚部から離れていて、精度よく抵抗成分と容量成分をそれぞれ切り分けができない場合がある。そこで、インピーダンス算出部６１が算出した周波数ｆごとの内部インピーダンスＺの実部及び虚部を基準にして、等価回路モデルの合成インピーダンスと比較することで、拡散抵抗１６の抵抗成分Ｒｗｉ及び容量成分Ｃｗｉを算出する。
Next, a method for calculating the resistance component and the capacitance component of the diffused resistor will be described. Since the diffused
図７に拡散抵抗１６の抵抗成分Ｒｗｉ及び容量成分Ｃｗｉの算出処理のフローチャートを示す。また、内部インピーダンスＺ（ｆ）として、車両の走行終了時点（蓄電池が充放電状態から充放電停止状態になった時点）及び停車中（充放電停止中）における充電電流の検出値Ｉ及び端子間電圧の検出値Ｖからウェーブレット変換により算出されたものが用いられる。
FIG. 7 shows a flowchart of the calculation process of the resistance component Rwi and the capacitance component Cwi of the
まずは、ステップＳ１１において、制御装置５０に記憶した抵抗成分Ｒｗｉ及び容量成分Ｃｗｉの初期値を設定する。この初期値は、演算速度を速めるために、初回以降は前回値を適用する。ステップＳ１２において、測定値に相当する各周波数ｆにおける内部インピーダンスＺ（ｆ）をインピーダンス算出部６１から取得する。
First, in step S11, initial values of the resistance component Rwi and the capacitance component Cwi stored in the control device 50 are set. As the initial value, the previous value is applied after the first time in order to increase the calculation speed. In step S <b> 12, the internal impedance Z (f) at each frequency f corresponding to the measurement value is acquired from the
次に、ステップＳ１３において、抵抗成分Ｒｗｉ及び容量成分Ｃｗｉの初期値に基づいて、各周波数ｆにおける推定内部インピーダンスＺａ（ｆ）を算出する。次に、ステップＳ１４において、初期値に基づいて算出された推定内部インピーダンスＺａ（ｆ）と、インピーダンス算出部６１によって算出された内部インピーダンスＺ（ｆ）の偏差を算出する。
Next, in step S13, the estimated internal impedance Za (f) at each frequency f is calculated based on the initial values of the resistance component Rwi and the capacitance component Cwi. Next, in step S14, a deviation between the estimated internal impedance Za (f) calculated based on the initial value and the internal impedance Z (f) calculated by the
次に、ステップＳ１５において、ステップＳ１４で算出された偏差が所定値以下か否かを判定する。誤差が最小となるものを選択するようにしてもよい。ステップＳ１４で算出された偏差が所定値より大きい場合（Ｓ１５：ＮＯ）、ステップＳ１６において、その偏差が減少するように抵抗成分Ｒｗｉ及び容量成分Ｃｗｉを変更する。そして、再び、ステップＳ１３において、各周波数ｆにおける推定内部インピーダンスＺａ（ｆ）の算出処理を行う。ステップＳ１４で算出された偏差が所定値以下の場合（Ｓ１５：ＹＥＳ）、ステップＳ１７において、抵抗成分Ｒｗｉ及び容量成分Ｃｗｉの値を現在値として記憶し、処理を終了する。 Next, in step S15, it is determined whether or not the deviation calculated in step S14 is a predetermined value or less. You may make it select the thing with the smallest error. When the deviation calculated in step S14 is larger than the predetermined value (S15: NO), in step S16, the resistance component Rwi and the capacitance component Cwi are changed so that the deviation decreases. Then, in step S13, the estimated internal impedance Za (f) is calculated at each frequency f. If the deviation calculated in step S14 is equal to or smaller than the predetermined value (S15: YES), the values of the resistance component Rwi and the capacitance component Cwi are stored as current values in step S17, and the process is terminated.
上記の処理によれば、インピーダンス算出部６１により算出された内部インピーダンスＺ（ｆ）の実部及び虚部と、蓄電池１０における等価回路モデルを構成する抵抗成分Ｒｗｉ及び容量成分Ｃｗｉの仮値に基づき算出された推定内部インピーダンスＺａ（ｆ）の実部及び虚部とを比較し、その比較結果に基づいて、抵抗成分Ｒｗｉ及び容量成分Ｃｗｉを算出することができる。
According to the above processing, based on the real and imaginary parts of the internal impedance Z (f) calculated by the
図８に、車両用高容量二次電池（リチウムイオン二次電池）において、合成インピーダンスである拡散抵抗１６について、測定値を用いてウェーブレット変換により算出した内部インピーダンスＺ（ｆ）と、使用する抵抗成分及び容量成分の合わせ込みを行った推定内部インピーダンスＺａ（ｆ）とを示す。推定内部インピーダンスＺａ（ｆ）として、２組の抵抗成分及び容量成分とみなして算出した推定内部インピーダンスＺａ（２ＲＣモデル）、及び、拡散抵抗１６を４組の抵抗成分及び容量成分とみなして算出した推定内部インピーダンスＺａ（４ＲＣモデル）を示す。
FIG. 8 shows the internal impedance Z (f) calculated by wavelet transform using the measured values for the diffused
推定内部インピーダンスＺａ（４ＲＣモデル）は、推定内部インピーダンスＺａ（２ＲＣモデル）と比較して、インピーダンス算出部６１によって算出された内部インピーダンスＺを精度よく合わせ込むことができる。このようにして、合成インピーダンスである拡散抵抗１６のモデルの構成を選定することができ、また、複数の抵抗成分と複数の容量成分を算出することができる。
As compared with the estimated internal impedance Za (2RC model), the estimated internal impedance Za (4RC model) can match the internal impedance Z calculated by the
次に、ＳＯＣ処理部５２について述べる。図２に示す構成において、ＳＯＣ処理部５２のＯＣＶ算出部６３には、ＲＣ成分算出部６２から抵抗成分Ｒｓ，Ｒ１，Ｒ２，Ｒｗｉの値及び容量成分Ｃ１，Ｃ２，Ｃｗｉの値が入力される。また、ＯＣＶ算出部６３には、電流センサ３０から現在の充放電電流の検出値Ｉ、電圧センサ３１から現在の端子間電圧の検出値Ｖが入力される。ＯＣＶ算出部６３では、抵抗成分Ｒｓ，Ｒ１，Ｒ２，Ｒｗｉの値、容量成分Ｃ１，Ｃ２，Ｃｗｉの値に基づいて内部インピーダンスＺを形成し、現在の充放電電流の検出値Ｉ、及び、現在の端子間電圧の検出値Ｖに基づいて、蓄電池１０の開放端電圧ＯＣＶを算出する。
Next, the
具体的には、図４に示す等価回路モデルにおいて、蓄電池１０の開放端電圧ＯＣＶは、現在の端子間電圧の検出値Ｖ、現在の充放電電流の検出値Ｉに基づいて、下記の式（６）のように算出することができる。なお、ｓは微分演算子である。ｍは反応抵抗１４，１５のＲＣ並列回路の個数であり、例えばｍ＝２とする。ｎは拡散抵抗１６のＲＣ並列回路の個数であり、合成インピーダンスの再現性を考慮して、例えば、ｍ＝４とする。
Specifically, in the equivalent circuit model shown in FIG. 4, the open-circuit voltage OCV of the
ＳＯＣ処理部５２のＯＣＶ−ＳＯＣ変換部６４は、図３に示すＯＣＶ−ＳＯＣ特性をマップとして記憶しておき、ＯＣＶ算出部６３によって算出された開放端電圧ＯＣＶをＳＯＣに変換する。このようにしてＳＯＣ処理部５２によって算出されたＳＯＣを用いることで、蓄電池１０に対する好適な充放電制御を行うことが可能になる。
The OCV-
図９に、本発明の効果として、車両用高容量二次電池（リチウムイオン二次電池）において、車両走行時の波形を用いて算出したＳＯＣの時間変化を示す。但し、Ｓ１：ＳＯＣの真値、Ｓ２：本実施形態のＳＯＣ処理部５２により算出されるＳＯＣ、Ｓ３：本実施形態の直流抵抗１３及び反応抵抗１４，１５のみ（拡散抵抗１６を含まない等価回路モデル）を用いて算出されるＳＯＣ、Ｓ４:電流Ｉと電圧Ｖについて最小２乗法を適用し切片を開放端電圧ＯＣＶとして得られたＳＯＣ（ＩＶ測定法によるＳＯＣ）である。
FIG. 9 shows, as an effect of the present invention, the time change of the SOC calculated using the waveform when the vehicle travels in the high-capacity secondary battery for vehicles (lithium ion secondary battery). However, S1: True value of SOC, S2: SOC calculated by the
時刻Ｔ１〜Ｔ２において、車両はＥＶ走行モードで走行し、時刻Ｔ２以降において、車両はＨＶ走行モードで走行している。時刻Ｔ１〜Ｔ２におけるＥＶ走行モードでは、拡散抵抗１６を含まない等価回路モデルを用いて算出されるＳＯＣ（Ｓ３）及びＩＶ測定法を用いて算出されるＳＯＣ（Ｓ４）は、ＥＶ走行時において誤差が大きい。これは、ＥＶ走行モードに代表される高電流が充放電の一方に継続して流れる場合に、拡散現象が顕著となり、見かけ上の電圧低下を引き起こすためである。つまり、電圧Ｖを入力とするＳＯＣを算出する場合、ＥＶ走行モードにおいては拡散抵抗１６を考慮しなければならない。従って、拡散抵抗１６を含まない等価回路モデル及びＩＶ測定法はＥＶ走行モードには適さない。
From time T1 to T2, the vehicle travels in the EV travel mode, and after time T2, the vehicle travels in the HV travel mode. In the EV travel mode at times T1 to T2, the SOC (S3) calculated using the equivalent circuit model not including the diffused
一方、本実施形態のＳＯＣ処理部５２によって算出されるＳＯＣ（Ｓ２）は、ＥＶ走行モードにおいて、ＳＯＣの真値（Ｓ１）に精度よく追随する。また、ＥＶ走行モードからＨＶ走行モードに切り替わると、ＩＶ測定法によるＳＯＣ（Ｓ４）は拡散抵抗１６が解消しＳＯＣの真値（Ｓ１）に戻っていく。つまり、ＥＶ走行モードからＨＶ走行モードに切り替わった後の電圧変化から、拡散抵抗１６のインピーダンスを算出できることを示している。よって、車両で走行中でも拡散抵抗１６の成分が求まる。
On the other hand, the SOC (S2) calculated by the
次に、劣化度処理部５３について説明する。民生電池（リチウムイオン二次電池）での劣化耐久試験において、直流抵抗１３、時定数の小さい順に反応抵抗１４，１５、及び拡散抵抗１６で分類した場合での各抵抗成及び各容量成分の推移について説明する。各抵抗成及び各容量成分の値は、試験前の値に比べて、劣化によりいずれも増加する。また、それぞれ異なった増加率（劣化度）を示す。このように、各抵抗成分及び各容量成分の劣化度が異なるので、劣化度を好適に見積もるためには各抵抗成分及び各容量成分を個別に算出する必要がある。
Next, the deterioration degree processing unit 53 will be described. Transition of each resistance component and each capacitance component when categorized by
図２に示す構成において、劣化度処理部５３は、基準値算出部６５及び劣化度算出部６６を備える。基準値算出部６５では、ＲＣ成分の基準値として、無劣化状態での各抵抗成分の値Ｒｓｉｎｉ，Ｒ１ｉｎｉ，Ｒ２ｉｎｉ,Ｒｗｉｉｎｉ及び各容量成分の値Ｃ１ｉｎｉ，Ｃ２ｉｎｉ，Ｃｗｉｉｎｉの値を算出する。これらは蓄電池１０の温度Ｔと充電量ＳＯＣを引数とする。
In the configuration shown in FIG. 2, the degradation level processing unit 53 includes a reference
劣化度算出部６６には、現在の蓄電池１０の抵抗成分Ｒｓ，Ｒ１，Ｒ２，Ｒｗｉ及び容量成分Ｃ１，Ｃ２，ＣｗｉがＲＣ成分算出部６２から入力され、また、基準値算出部６５から基準抵抗成分及び基準容量成分が入力される。劣化度算出部６６は、抵抗成分と基準抵抗成分との比、及び、容量成分と基準容量成分との比として蓄電池１０の劣化度Ｋを算出する。
The degradation
次に、予測電力処理部５４について説明する。予測電圧算出部６８は、現在の抵抗成分Ｒと容量成分Ｃ、基準電流Ｉｕ、開放端電圧ＯＣＶ、式（６）に示す抵抗成分と容量成分の並列回路（内部インピーダンス）の印加電圧ΔＶｃを用いて、式（６）により、１秒後の予測電圧Ｖ（ｔ＋ｔ１）を推定する。予測電流に相当する基準電流Ｉｕは、例えば最大許容電流Ｉｍａｘとする。並列回路の印加電圧ΔＶｃは現在の容量成分の電荷量に応じて変化する電圧変化であり、式（６）の右辺第３項と第４項に対応する。例えば、ＳＯＣ処理部５２のＯＣＶ算出部６３の現在の値を用いて算出することができる。ｔ１秒後の予測電力Ｐ（ｔ＋ｔ１）＝ｔ１秒後の予測電圧Ｖ（ｔ＋ｔ１）×基準電流Ｉｕとして算出する。ｔ１秒後の予測電力Ｐ（ｔ＋ｔ１）と蓄電池１０の放電許容電力Ｗｏｕｔまたは充電許容電力Ｗｉｎと比較し、好適な制御方法を決定する。
Next, the predicted
以下、本実施形態における効果を述べる。 Hereinafter, effects in the present embodiment will be described.
ウェーブレット変換を用いて、蓄電池１０の周波数毎の内部インピーダンスＺの実部Ｒｅ（Ｚ（ｆ））と虚部Ｉｍ（Ｚ（ｆ））を算出できる。そして、その算出された周波数毎の内部インピーダンスＺの実部Ｒｅ（Ｚ（ｆ））と虚部Ｉｍ（Ｚ（ｆ））とに基づいて、複数の時定数にあたる複数の抵抗成分及び複数の容量成分を算出できる。
The real part Re (Z (f)) and imaginary part Im (Z (f)) of the internal impedance Z for each frequency of the
ウェーブレット変換は、変換を行う対象が非周期的な波形であっても周波数毎の内部インピーダンスＺを算出できる。これにより、車両に搭載された蓄電池の波形から内部インピーダンスＺが算出できる。特に、車両環境下での複雑な充放電の波形、充放電状態から充放電停止に切替わる点とその後の充放電停止中の緩電圧変化の波形、矩形状の波形において内部インピーダンスＺが算出できる。勿論、正弦波でも内部インピーダンスＺが算出できる。 The wavelet transform can calculate the internal impedance Z for each frequency even if the object to be transformed is an aperiodic waveform. Thereby, the internal impedance Z can be calculated from the waveform of the storage battery mounted on the vehicle. In particular, the internal impedance Z can be calculated in a complicated charging / discharging waveform under a vehicle environment, a point where switching from charging / discharging state to charging / discharging stop and a subsequent slow voltage change waveform during charging / discharging stopping, and a rectangular waveform. . Of course, the internal impedance Z can also be calculated with a sine wave.
内部インピーダンスＺを直流抵抗１３、反応抵抗１４，１５、及び、拡散抵抗１６の抵抗成分及び容量成分を算出できることで、例えば、蓄電池のＳＯＣ、予測電力、また、蓄電池１０の劣化度を精度よく算出できる。蓄電池内部で発生する個々の劣化モードの特定が成し得る。
The internal impedance Z can be calculated from the resistance component and capacity component of the
また、内部インピーダンスの実部Ｒｅ（Ｚ（ｆ））及び虚部Ｉｍ（Ｚ（ｆ））から得られることで、ベクトル軌跡線図に基づいて、各抵抗成分及び各容量成分を算出できる。 Further, by obtaining from the real part Re (Z (f)) and the imaginary part Im (Z (f)) of the internal impedance, each resistance component and each capacitance component can be calculated based on the vector locus diagram.
蓄電池１０の充放電に伴う開放端電圧ＯＣＶの変化が、内部インピーダンスＺの算出に与える影響を抑制すべく、端子間電圧の検出値Ｖを開放端電圧ＯＣＶの変化量ΔＯＣＶで補正することで、内部インピーダンスＺを精度よく算出できる。
By correcting the detected value V of the inter-terminal voltage with the change amount ΔOCV of the open-end voltage OCV in order to suppress the influence of the change in the open-end voltage OCV accompanying the charge / discharge of the
（その他の実施形態）
・上記実施形態では、内部インピーダンスＺを直流抵抗１３、反応抵抗１４，１５、拡散抵抗１６の直列接続体としてモデル化したが、これを変更してもよい。例えば、蓄電池の特性や周囲環境を制約し、反応抵抗１４，１５、拡散抵抗１６の発現を抑制する事ができれば、いずれかを省略する構成としてもよい。
(Other embodiments)
In the above embodiment, the internal impedance Z is modeled as a series connection body of the
・上記実施形態では、充放電状態から充放電停止に切替わる点とその後の充放電停止中の緩電圧変化の波形から拡散抵抗を算出していたが、これを変更してもよい。例えば、充放電中に限定して、平均電流が大きい充放電から平均電流が小さい充放電に切り替わる点とその後の平均電流が小さい充放電中における緩電圧変化の波形から拡散抵抗を算出してもよい。具体的には、ＥＶ走行からＨＶ走行に切り替わる点とその後のＨＶ走行において、電圧降下が回復する際に拡散抵抗を算出することができる。これにより、走行中において電圧変化が顕著な状態であるため、走行中に精度よく拡散抵抗を算出することができる。また、例えば、充電停止状態から充放電状態に切り替わる点とその後の充放電中の緩電圧変化の波形から拡散抵抗を算出してもよい。具体的には二次電池に対して急速充電を行う場合としてもよい。上記実施形態を含め、前述の波形において、直流抵抗、反応抵抗、及び拡散抵抗の少なくとも１つを算出するとしてもよい。 In the above embodiment, the diffusion resistance is calculated from the point of switching from the charge / discharge state to the charge / discharge stop and the waveform of the slow voltage change during the subsequent charge / discharge stop, but this may be changed. For example, the diffusion resistance can be calculated from the point of switching from charging / discharging with a large average current to charging / discharging with a small average current and the waveform of the slow voltage change during charging / discharging with a small average current, during charging / discharging. Good. Specifically, the diffusion resistance can be calculated when the voltage drop recovers at the point where the EV traveling is switched to the HV traveling and the subsequent HV traveling. Thereby, since the voltage change is remarkable during traveling, the diffusion resistance can be accurately calculated during traveling. Further, for example, the diffusion resistance may be calculated from the point at which the charge stop state is switched to the charge / discharge state and the waveform of the slow voltage change during the subsequent charge / discharge. Specifically, the secondary battery may be rapidly charged. In the above-described waveform including the above embodiment, at least one of direct current resistance, reaction resistance, and diffusion resistance may be calculated.
図９に示すように、車両用高容量二次電池（リチウムイオン二次電池）において、拡散抵抗１６を含む等価回路モデルを用いて算出したＳＯＣ（Ｓ２）と、拡散抵抗１６を含まない等価回路モデルを用いて算出したＳＯＣ（Ｓ３）との差は、時刻Ｔ１から増加し続け、時刻Ｔ２において最も大きくなる。そして、時刻Ｔ２の後、Ｓ２とＳ３との差が解消していく。これは、ＥＶ走行に伴って、蓄電池１０から大電流が出力され続けることで、拡散抵抗１６に生じる分極が大きくなり、その後、ＨＶ走行に切り替わることでその分極が解消されるためと考えられる。つまり、ＥＶ走行からＨＶ走行に切り替わった時点（時刻Ｔ２）及びその後拡散抵抗１６に生じる分極が解消するまでの間において、拡散抵抗１６による電圧Ｖへの影響が顕著となる。そこで、ＥＶ走行からＨＶ走行に切り替わった時点（時刻Ｔ２）及びその後拡散抵抗１６に生じる分極が解消するまでの間において検出された電流Ｉ及び電圧Ｖに基づいて、内部インピーダンスを算出する構成とすることで、拡散抵抗１６を精度よく算出することが可能になる。
As shown in FIG. 9, in a high capacity secondary battery for vehicles (lithium ion secondary battery), SOC (S2) calculated using an equivalent circuit model including a diffused
・上記実施形態では、蓄電池１０が搭載されている車両とはハイブリッド車を示したが、これを変更してもよい。例えば、蓄電池を含む電源システムを要するいずれの車両に適用してもよい。同様に、車両に限らず、蓄電池を含む電源システムを要するいずれのシステムに適用してもよい。
-In the said embodiment, although the vehicle in which the
・上記実施形態では、蓄電池１０としてリチウムイオン二次電池を用いたが、これを変更してもよい。例えば、ニッケル水素蓄電池を用いてもよい。
In the above embodiment, a lithium ion secondary battery is used as the
・上記実施形態では、充放電に伴う充電量ＳＯＣの変化による開放端電圧ＯＣＶの変化量ΔＯＣＶを算出し、その変化量ΔＯＣＶを用いて端子間電圧の検出値Ｖを補正する構成とした。これを変更してもよい。以下に別の補正方法として、内部インピーダンスＺに対して補正する方法を述べる。図３に示す充電量ＳＯＣと開放端電圧ＯＣＶとの関係を、内部インピーダンスの等価容量成分Ｃｏｃｖとして予め変換しておき、制御装置５０に記憶しておく。等価容量成分Ｃｏｃｖは、例えば、Ｃｏｃｖ＝Ａｆ・ΔＳＯＣ／ΔＯＣＶとして算出できる（ここで、Ａｆは蓄電池１０の満充電容量）。内部インピーダンスＺの算出に際し、開放端電圧の変化量ΔＯＣＶが含まれた状態で算出した内部インピーダンスＺに対して、補正項としてＣｏｃｖを減算（−１／ｓＣｏｃｖ）すればよい（但し、ｓは微分演算子）。これにより、蓄電池１０の内部インピーダンスＺとは異なる成分である充電量変化による開放端電圧変化の影響を取り除くことができる。この他に、等価回路モデルでＳＯＣを算出する際に、式（６）の等価回路モデル式に補正項として等価容量成分Ｃｏｃｖを減算（−１／ｓＣｏｃｖ）すればよい。図３に示す充電量ＳＯＣと開放端電圧ＯＣＶの関係は、蓄電池の劣化や蓄電池の温度変化に対して殆ど影響を受けず特性変化がないものとして扱えるため、等価容量成分Ｃｏｃｖは予め定数に変換して制御装置５０に記憶しておけばよい。また、等価容量成分Ｃｏｃｖは、ＳＯＣに依らない定数として制御装置５０に記憶してもよく、ＳＯＣをパラメータとするマップとして制御装置５０に記憶しておいてもよい。ＳＯＣをパラメータとするマップにより等価容量成分Ｃｏｃｖを算出する方法は、充電量ＳＯＣに基づいて内部インピーダンスＺの虚部を補正していると言える。これら等価容量成分Ｃｏｃｖによる内部インピーダンスの補正方法は、制御装置５０の記憶量や演算量が比較的少なくてすむ。
In the above embodiment, the change amount ΔOCV of the open circuit voltage OCV due to the change in the charge amount SOC accompanying charge / discharge is calculated, and the detected value V of the terminal voltage is corrected using the change amount ΔOCV. This may be changed. As another correction method, a method for correcting the internal impedance Z will be described below. The relationship between the charge amount SOC and the open-circuit voltage OCV shown in FIG. 3 is converted in advance as an equivalent capacitance component Cocv of internal impedance and stored in the control device 50. The equivalent capacity component Cocv can be calculated, for example, as Cocv = Af · ΔSOC / ΔOCV (where Af is the full charge capacity of the storage battery 10). When calculating the internal impedance Z, Cocv may be subtracted (−1 / sCocv) as a correction term with respect to the internal impedance Z calculated in a state in which the open-circuit voltage variation ΔOCV is included (where s is a differential) operator). Thereby, the influence of the open end voltage change by the charge amount change which is a component different from the internal impedance Z of the
・上記実施形態では、ＲＣ成分算出部６２で、ベクトル軌跡線図による抵抗成分及び容量成分の算出方法を用いていたが、これを変更してもよい。ベクトル軌跡線図に代えて、ボード線図による抵抗成分及び容量成分の算出方法を用いてもよい。ボード線図による算出方法を述べる。内部インピーダンスの実部と虚部の絶対値について周波数特性を求める。周波数特性の形状は、階段状になる。１段ずつの距離が直流抵抗１３、反応抵抗１４，１５の各抵抗成分にあたる。また、段の折れ点での周波数（折れ点周波数）と先に求めた各抵抗成分から各容量成分が算出できる。ＲＣ成分算出方法として、ベクトル軌跡線図での方法とボード線図での方法のいずれも、周波数毎の内部インピーダンスの実部と虚部が必要になる。
In the above embodiment, the RC
・直流抵抗１３、反応抵抗１４，１５、拡散抵抗１６について、ＲＣ成分算出部６２における方法を変更してもよい。例えば、ＲＣ成分算出方法１とＲＣ成分算出方法２の適用対象を変更してもよい。具体的には、全抵抗成分及び全容量成分について、ＲＣ成分算出方法１を用いて算出してもよい。また、全抵抗成分及び全容量成分について、ＲＣ成分算出方法２を用いて算出してもよい。
The method in the RC
・上記実施形態の内部状態推定装置は、蓄電池１０は、車両に搭載された状態において適用したが、これを変更してもよい。例えば、図１０に示すように蓄電池１０が電池特性計測装置２３に接続した状態に変更して、内部状態推定装置を適用してもよい。蓄電池１０が電池特性計測装置２３に接続した状態で、車両での蓄電池１０に対する非周期的な充放電を模した電流又は電圧の波形を電池特性計測装置から蓄電池に対して充放電を行ってもよい。この場合、内部状態推定装置は制御装置５０に変えて、電池特性計測装置２３が持つ制御装置に適用してもよい。これにより、蓄電池１０が車両に搭載された状態を模擬でき、走行中の内部インピーダンスの推移を手軽に計測できる。
-Although the
・ＳＯＣ処理部５２及び劣化度処理部５３、予測電力処理部５４の少なくともいずれかを省略する構成としてもよい。また、電圧補正部５６を省略する構成としてもよい。
-It is good also as a structure which abbreviate | omits at least any one of the
１０…蓄電池（二次電池）、３０…電流センサ、３１…電圧センサ、５４…記憶部（記憶手段）、６０…ウェーブレット変換部（ウェーブレット変換手段）、６１…インピーダンス算出部（内部インピーダンス算出手段）。
DESCRIPTION OF
Claims (20)
前記時系列で記憶された端子間電圧の検出値に対してウェーブレット変換を行うことで電圧ウェーブレット変換係数を算出するとともに、前記時系列で記憶された充放電電流の検出値に対してウェーブレット変換を行うことで電流ウェーブレット変換係数を算出するウェーブレット変換手段（６０）と、
前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、周波数ごとの前記二次電池の内部インピーダンスの実部及び虚部を算出する内部インピーダンス算出手段（６１）と、
を備え、
前記二次電池は、車両において走行用モータの電力源として搭載されるものであって、
前記二次電池には、電池特性計測装置から前記車両の動作状態及び停止状態での前記二次電池における非周期的な充放電を模した電流又は電圧が入力され、
前記ウェーブレット変換手段は、前記電池特性計測装置により前記非周期的な充放電を行う場合おいて、前記電圧検出手段及び前記電流検出手段によって検出された検出値に対してウェーブレット変換を行うことを特徴とする電池内部状態推定装置。 Storage means for storing the detection values by the voltage detection means (31) for detecting the voltage between the terminals of the secondary battery (10) and the current detection means (30) for detecting the charge / discharge current flowing in the secondary battery in time series ( 5 5 )
A voltage wavelet transform coefficient is calculated by performing wavelet transform on the detected value of the inter-terminal voltage stored in the time series, and wavelet transform is performed on the detected charge / discharge current value stored in the time series. Wavelet transform means (60) for calculating a current wavelet transform coefficient by performing,
An internal impedance calculating means (61) for calculating a real part and an imaginary part of the internal impedance of the secondary battery for each frequency according to a ratio of the voltage wavelet transform coefficient and the current wavelet transform coefficient;
Equipped with a,
The secondary battery is mounted as a power source of a driving motor in a vehicle,
The secondary battery is supplied with a current or voltage imitating non-periodic charge / discharge in the secondary battery in the operation state and stop state of the vehicle from a battery characteristic measurement device,
The wavelet transform unit performs wavelet transform on the detection values detected by the voltage detection unit and the current detection unit when the non-periodic charge / discharge is performed by the battery characteristic measurement device. A battery internal state estimating device.
前記時系列で記憶された端子間電圧の検出値に対してウェーブレット変換を行うことで電圧ウェーブレット変換係数を算出するとともに、前記時系列で記憶された充放電電流の検出値に対してウェーブレット変換を行うことで電流ウェーブレット変換係数を算出するウェーブレット変換手段（６０）と、
前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、周波数ごとの前記二次電池の内部インピーダンスの実部及び虚部を算出する内部インピーダンス算出手段（６１）と、
を備え、
前記二次電池は、車両において走行用モータの電力源として搭載されるものであって、
前記車両は、前記走行用モータ及び内燃機関を駆動力として用いるとともに、前記走行用モータのみを駆動力とするモータ走行と、前記走行用モータ及び内燃機関の両方を駆動力とするハイブリッド走行とを行い、
前記ウェーブレット変換手段は、前記車両において前記モータ走行から前記ハイブリッド走行への切り替えが行われた時点、及び、その切り替え後のハイブリット走行中において前記電圧検出手段及び前記電流検出手段によって検出された検出値に対してウェーブレット変換を行うことで前記電圧ウェーブレット変換係数及び前記電流ウェーブレット変換係数を算出し、
前記内部インピーダンス算出手段は、前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、前記二次電池の内部インピーダンスの実部及び虚部を算出することを特徴とする電池内部状態推定装置。 Storage means for storing the detection values by the voltage detection means (31) for detecting the voltage between the terminals of the secondary battery (10) and the current detection means (30) for detecting the charge / discharge current flowing in the secondary battery in time series ( 55) and
A voltage wavelet transform coefficient is calculated by performing wavelet transform on the detected value of the inter-terminal voltage stored in the time series, and wavelet transform is performed on the detected charge / discharge current value stored in the time series. Wavelet transform means (60) for calculating a current wavelet transform coefficient by performing,
An internal impedance calculating means (61) for calculating a real part and an imaginary part of the internal impedance of the secondary battery for each frequency according to a ratio of the voltage wavelet transform coefficient and the current wavelet transform coefficient;
With
The secondary battery is mounted as a power source of a driving motor in a vehicle,
The vehicle uses the traveling motor and the internal combustion engine as driving power, and motor traveling using only the traveling motor as driving force, and hybrid traveling using both the traveling motor and internal combustion engine as driving force. Done
The wavelet transform unit detects the detected value detected by the voltage detection unit and the current detection unit during the hybrid travel after the switching from the motor travel to the hybrid travel in the vehicle. The voltage wavelet transform coefficient and the current wavelet transform coefficient are calculated by performing wavelet transform on the
The internal impedance calculation unit calculates a real part and an imaginary part of the internal impedance of the secondary battery based on a ratio between the voltage wavelet transform coefficient and the current wavelet transform coefficient .
前記時系列で記憶された端子間電圧の検出値に対してウェーブレット変換を行うことで電圧ウェーブレット変換係数を算出するとともに、前記時系列で記憶された充放電電流の検出値に対してウェーブレット変換を行うことで電流ウェーブレット変換係数を算出するウェーブレット変換手段（６０）と、
前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、周波数ごとの前記二次電池の内部インピーダンスの実部及び虚部を算出する内部インピーダンス算出手段（６１）と、
を備え、
前記二次電池の充放電による充電量変化に基づいて、前記電圧検出手段による端子間電圧の検出値を補正する電圧補正手段（５６）を備えることを特徴とする電池内部状態推定装置。 Storage means for storing the detection values by the voltage detection means (31) for detecting the voltage between the terminals of the secondary battery (10) and the current detection means (30) for detecting the charge / discharge current flowing in the secondary battery in time series ( 55) and
A voltage wavelet transform coefficient is calculated by performing wavelet transform on the detected value of the inter-terminal voltage stored in the time series, and wavelet transform is performed on the detected charge / discharge current value stored in the time series. Wavelet transform means (60) for calculating a current wavelet transform coefficient by performing,
An internal impedance calculating means (61) for calculating a real part and an imaginary part of the internal impedance of the secondary battery for each frequency according to a ratio of the voltage wavelet transform coefficient and the current wavelet transform coefficient;
With
A battery internal state estimation device comprising: a voltage correction means (56) for correcting a detected value of the voltage between the terminals by the voltage detection means based on a change in charge amount due to charging / discharging of the secondary battery .
前記時系列で記憶された端子間電圧の検出値に対してウェーブレット変換を行うことで電圧ウェーブレット変換係数を算出するとともに、前記時系列で記憶された充放電電流の検出値に対してウェーブレット変換を行うことで電流ウェーブレット変換係数を算出するウェーブレット変換手段（６０）と、A voltage wavelet transform coefficient is calculated by performing wavelet transform on the detected value of the inter-terminal voltage stored in the time series, and wavelet transform is performed on the detected charge / discharge current value stored in the time series. Wavelet transform means (60) for calculating a current wavelet transform coefficient by performing,
前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、周波数ごとの前記二次電池の内部インピーダンスの実部及び虚部を算出する内部インピーダンス算出手段（６１）と、An internal impedance calculating means (61) for calculating a real part and an imaginary part of the internal impedance of the secondary battery for each frequency according to a ratio of the voltage wavelet transform coefficient and the current wavelet transform coefficient;
を備え、With
前記内部インピーダンス算出手段は、前記内部インピーダンスの算出に際し、前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により算出された前記内部インピーダンスの虚部を、前記二次電池の充放電による充電量変化に基づいて補正する補正手段を備えることを特徴とする電池内部状態推定装置。The internal impedance calculation means, when calculating the internal impedance, the imaginary part of the internal impedance calculated by the ratio of the voltage wavelet transform coefficient and the current wavelet transform coefficient, the charge amount due to charging and discharging of the secondary battery A battery internal state estimation device comprising correction means for correcting based on a change.
前記時系列で記憶された端子間電圧の検出値に対してウェーブレット変換を行うことで電圧ウェーブレット変換係数を算出するとともに、前記時系列で記憶された充放電電流の検出値に対してウェーブレット変換を行うことで電流ウェーブレット変換係数を算出するウェーブレット変換手段（６０）と、A voltage wavelet transform coefficient is calculated by performing wavelet transform on the detected value of the inter-terminal voltage stored in the time series, and wavelet transform is performed on the detected charge / discharge current value stored in the time series. Wavelet transform means (60) for calculating a current wavelet transform coefficient by performing,
前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、周波数ごとの前記二次電池の内部インピーダンスの実部及び虚部を算出する内部インピーダンス算出手段（６１）と、An internal impedance calculating means (61) for calculating a real part and an imaginary part of the internal impedance of the secondary battery for each frequency according to a ratio of the voltage wavelet transform coefficient and the current wavelet transform coefficient;
を備え、With
前記内部インピーダンス算出手段は、前記内部インピーダンスの算出に際し、前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により算出された前記内部インピーダンスの虚部を、前記二次電池の充放電による充電量変化に伴う前記二次電池の開放端電圧の変化を表す等価容量成分に基づいて補正する補正手段を備えることを特徴とする電池内部状態推定装置。The internal impedance calculation means, when calculating the internal impedance, the imaginary part of the internal impedance calculated by the ratio of the voltage wavelet transform coefficient and the current wavelet transform coefficient, the charge amount due to charging and discharging of the secondary battery A battery internal state estimation device comprising correction means for correcting based on an equivalent capacity component representing a change in open circuit voltage of the secondary battery accompanying a change.
前記ウェーブレット変換手段は、前記二次電池に対して非周期的な充放電を行う場合において、前記電圧検出手段及び前記電流検出手段によって検出された検出値に対してウェーブレット変換を行うことで前記電圧ウェーブレット変換係数及び前記電流ウェーブレット変換係数を算出し、
前記内部インピーダンス算出手段は、前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、前記二次電池の内部インピーダンスの実部及び虚部を算出することを特徴とする請求項２乃至８のいずれか１項に記載の電池内部状態推定装置。 The secondary battery is mounted as a power source of a driving motor in a vehicle, and charging / discharging is performed on the secondary battery according to a driving state of the vehicle,
The wavelet transform unit performs the wavelet transform on the detection value detected by the voltage detection unit and the current detection unit when performing non-periodic charge / discharge on the secondary battery. Calculating a wavelet transform coefficient and the current wavelet transform coefficient;
9. The internal impedance calculation unit calculates a real part and an imaginary part of the internal impedance of the secondary battery based on a ratio between the voltage wavelet transform coefficient and the current wavelet transform coefficient . The battery internal state estimation apparatus according to any one of the preceding claims.
前記ウェーブレット変換手段は、前記二次電池の充放電停止時点及び充放電停止中において、前記電圧検出手段及び前記電流検出手段によって検出された検出値に対してウェーブレット変換を行うことで前記電圧ウェーブレット変換係数及び前記電流ウェーブレット変換係数を算出し、
前記内部インピーダンス算出手段は、前記電圧ウェーブレット変換係数と前記電流ウェーブレット変換係数との比により、前記二次電池の内部インピーダンスの実部及び虚部を算出することを特徴とする請求項２乃至９のいずれか１項に記載の電池内部状態推定装置。 The secondary battery is mounted as a power source of a driving motor in a vehicle,
The wavelet transform unit performs the wavelet transform on the detection values detected by the voltage detection unit and the current detection unit during the charge / discharge stop time and the charge / discharge stop of the secondary battery, thereby performing the voltage wavelet transform. Calculating the coefficient and the current wavelet transform coefficient;
10. The internal impedance calculation unit calculates a real part and an imaginary part of the internal impedance of the secondary battery based on a ratio between the voltage wavelet transform coefficient and the current wavelet transform coefficient . The battery internal state estimation apparatus according to any one of the preceding claims.
前記直流抵抗、前記反応抵抗、及び、前記拡散抵抗はそれぞれ時定数が相互に異なるものであって、
前記成分算出手段は、前記内部インピーダンス算出手段により算出された前記内部インピーダンスの実部及び虚部に基づいて、前記直流抵抗、前記反応抵抗及び前記拡散抵抗をそれぞれ構成する抵抗成分及び容量成分の少なくともいずれかを算出することを特徴とする請求項１１に記載の電池内部状態推定装置。 The internal impedance comprises a direct current resistance representing a conduction resistance, a reaction resistance representing a reaction at an electrode interface, and a diffusion resistance representing ion diffusion,
The DC resistance, the reaction resistance, and the diffusion resistance have different time constants, respectively.
The component calculating means is based on a real part and an imaginary part of the internal impedance calculated by the internal impedance calculating means, and includes at least a resistance component and a capacitance component that respectively constitute the DC resistance, the reaction resistance, and the diffusion resistance. The battery internal state estimation device according to claim 11 , wherein either one is calculated.
前記電圧検出手段及び前記電流検出手段による検出値と、前記成分算出手段によって算出された前記抵抗成分及び前記容量成分とに基づいて、前記二次電池の開放端電圧を算出する開放端電圧算出手段（６３）を備えることを特徴とする請求項１１乃至１６のいずれか１項に記載の電池内部状態推定装置。 The component calculation means calculates the resistance component and the capacitance component,
An open-end voltage calculation means for calculating an open-end voltage of the secondary battery based on the detection value by the voltage detection means and the current detection means and the resistance component and the capacity component calculated by the component calculation means. The battery internal state estimation device according to any one of claims 11 to 16 , further comprising (63).
前記成分算出手段により算出された前記抵抗成分及び前記容量成分と、前記容量成分への電荷の蓄積により生じる電圧の現在値と、前記二次電池の開放端電圧の現在値と、将来の所定時刻において前記二次電池に流れる充放電電流の予測値である予測電流と、に基づいて、その将来の所定時刻における前記二次電池の端子間電圧の予測値である予測電圧を算出する予測電圧算出手段（６８）と、
前記予測電流及び前記予測電圧に基づいて、前記所定時刻における前記二次電池における充放電電力の予測値である予測電力を算出する予測電力算出手段（５４）を備えることを特徴とする請求項１１乃至１９のいずれか１項に記載の電池内部状態推定装置。 The component calculation means calculates the resistance component and the capacitance component,
The resistance component and the capacitance component calculated by the component calculation means, the current value of the voltage generated by the charge accumulation in the capacitance component, the current value of the open-circuit voltage of the secondary battery, and a predetermined future time A predicted voltage that calculates a predicted voltage that is a predicted value of the inter-terminal voltage of the secondary battery at a predetermined future time based on a predicted current that is a predicted value of a charge / discharge current flowing in the secondary battery Means (68);
Based on the predicted current and the predicted voltage, claim 11, characterized in that it comprises a predicted power calculating means (54) for calculating the predicted power is the predicted value of the charge-discharge electric power in the secondary battery in the predetermined time The battery internal state estimation device according to any one of 1 to 19 .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2014107775A JP6277864B2 (en) | 2014-05-26 | 2014-05-26 | Battery internal state estimation device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2014107775A JP6277864B2 (en) | 2014-05-26 | 2014-05-26 | Battery internal state estimation device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2015224876A JP2015224876A (en) | 2015-12-14 |
| JP6277864B2 true JP6277864B2 (en) | 2018-02-14 |
Family
ID=54841750
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2014107775A Active JP6277864B2 (en) | 2014-05-26 | 2014-05-26 | Battery internal state estimation device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6277864B2 (en) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6326452B2 (en) * | 2016-06-15 | 2018-05-16 | 本田技研工業株式会社 | Battery state estimation device and battery state estimation method |
| JP6790892B2 (en) * | 2017-02-16 | 2020-11-25 | 株式会社デンソー | Battery state estimation device and power supply device |
| JP6969464B2 (en) * | 2018-03-19 | 2021-11-24 | トヨタ自動車株式会社 | Deterioration state estimation method for secondary battery system and secondary battery |
| WO2019215786A1 (en) * | 2018-05-07 | 2019-11-14 | 三菱電機株式会社 | Cell degradation detection device and cell temperature estimation device |
| WO2020208762A1 (en) * | 2019-04-11 | 2020-10-15 | 三菱電機株式会社 | Device and method for assessing characteristics of storage battery |
| WO2021149774A1 (en) * | 2020-01-24 | 2021-07-29 | 株式会社デンソー | Battery measurement device |
| JP2021117221A (en) * | 2020-01-24 | 2021-08-10 | 株式会社デンソー | Battery measuring device |
| CN111624503B (en) * | 2020-04-26 | 2023-04-28 | 宁波普瑞均胜汽车电子有限公司 | Online estimation method for temperature of lithium ion battery |
| JP7405043B2 (en) * | 2020-09-02 | 2023-12-26 | トヨタ自動車株式会社 | Evaluation method for secondary batteries |
| KR102538419B1 (en) * | 2021-08-11 | 2023-05-31 | 한국생산기술연구원 | Method of measuring internal resistance of battery by current cut-off and method of measuring state of charge using same |
| JP2024008300A (en) * | 2022-07-07 | 2024-01-19 | 大和製罐株式会社 | Economical efficiency estimation device and economical efficiency estimation method using charge-discharge curve of secondary battery |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6556001B1 (en) * | 1999-10-12 | 2003-04-29 | Gerald Wiegand | Highly time resolved impedance spectroscopy |
| DE10012964A1 (en) * | 2000-03-16 | 2001-10-04 | Implex Hear Tech Ag | Device for operating re-chargeable electrical energy storage device selects current charging strategy for storage device depending on prognosis derived from model, storage device parameters |
| EP2472277A3 (en) * | 2003-06-27 | 2012-10-17 | The Furukawa Electric Co., Ltd. | Method and device for measuring secondary cell internal impedance and judging deterioration |
| JP5088081B2 (en) * | 2007-10-12 | 2012-12-05 | 富士通株式会社 | Battery measuring method and battery manufacturing method |
| JP2009106059A (en) * | 2007-10-23 | 2009-05-14 | Meidensha Corp | Battery simulator |
| JP5058959B2 (en) * | 2008-12-12 | 2012-10-24 | 愛三工業株式会社 | Secondary battery deterioration state diagnosis device |
| DE102009049589A1 (en) * | 2009-10-16 | 2011-04-21 | Bayerische Motoren Werke Aktiengesellschaft | Method for determining and / or predicting the maximum performance of a battery |
| JP2012083142A (en) * | 2010-10-07 | 2012-04-26 | Toshiba Corp | Calculation device for internal resistance of secondary battery |
| JP5307113B2 (en) * | 2010-12-20 | 2013-10-02 | 古河電気工業株式会社 | Full charge detection device and full charge detection method |
| JP6026120B2 (en) * | 2012-03-19 | 2016-11-16 | 株式会社東芝 | Secondary battery internal resistance calculation device and internal resistance calculation method thereof, secondary battery abnormality detection device and abnormality detection method thereof, secondary battery degradation estimation device and degradation estimation method thereof |
| JP5950728B2 (en) * | 2012-06-29 | 2016-07-13 | 日置電機株式会社 | Equivalent circuit analysis apparatus and equivalent circuit analysis method |
| JP6011135B2 (en) * | 2012-08-09 | 2016-10-19 | 株式会社デンソー | Power system |
| JP2014238302A (en) * | 2013-06-06 | 2014-12-18 | ニューロング精密工業株式会社 | Secondary battery impedance evaluation method and secondary battery impedance evaluation apparatus |
-
2014
- 2014-05-26 JP JP2014107775A patent/JP6277864B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2015224876A (en) | 2015-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6277864B2 (en) | Battery internal state estimation device | |
| CN110914696B (en) | Method and system for estimating battery open cell voltage, state of charge, and state of health during operation of a battery | |
| US10656210B2 (en) | Secondary battery state detection device and secondary battery state detection method | |
| US10641830B2 (en) | Battery's state of charge estimation apparatus | |
| US9685807B2 (en) | Battery control device | |
| JP6111275B2 (en) | Battery control device | |
| US20180321326A1 (en) | Method and system for estimating state of charge or depth of discharge of battery, and method and system for evaluating health of battery | |
| CN106133994B (en) | Secondary cell condition checkout gear and secondary cell condition detection method | |
| JP6197479B2 (en) | Power storage system and method for estimating full charge capacity of power storage device | |
| WO2014141796A1 (en) | Secondary cell state detection device and method for detecting secondary cell state | |
| JP2020504994A (en) | Hybrid battery charger / tester | |
| JP2006242880A (en) | Condition detector for power supply device, power supply device, and initial characteristic extractor used for power supply device | |
| KR20160120476A (en) | Apparatus and Method for caculating degradation degree | |
| JP2009081964A (en) | Power supply system for vehicle | |
| JP6382663B2 (en) | Battery state determination method and battery state determination device | |
| US10794961B2 (en) | Internal state estimating device | |
| JP2018096803A (en) | Internal resistance calculation device, method for calculating internal resistance, and internal resistance calculation program | |
| JP6183283B2 (en) | Parameter estimation device for equivalent circuit of secondary battery for vehicle | |
| JP2013088156A (en) | Secondary battery state detector and secondary battery state detection method | |
| JP2016161354A (en) | Deterioration detector and method for detecting deterioration | |
| Nejad et al. | Sensitivity of lumped parameter battery models to constituent parallel-RC element parameterisation error | |
| JP6668102B2 (en) | Deterioration detection device and deterioration detection method | |
| JP2016130658A (en) | Secondary battery control device | |
| JP6606153B2 (en) | Battery state estimation method and battery state estimation device | |
| JP2016161352A (en) | Deterioration detector and method for detecting deterioration |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160721 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170421 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170523 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170719 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20171219 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20180101 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 6277864 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |