[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP2008096328A - Charging rate estimation method, charging rate estimating device, and secondary battery power supply system - Google Patents

Charging rate estimation method, charging rate estimating device, and secondary battery power supply system Download PDF

Info

Publication number
JP2008096328A
JP2008096328A JP2006279478A JP2006279478A JP2008096328A JP 2008096328 A JP2008096328 A JP 2008096328A JP 2006279478 A JP2006279478 A JP 2006279478A JP 2006279478 A JP2006279478 A JP 2006279478A JP 2008096328 A JP2008096328 A JP 2008096328A
Authority
JP
Japan
Prior art keywords
voltage
secondary battery
charging rate
adjustment parameter
equation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2006279478A
Other languages
Japanese (ja)
Other versions
JP4785056B2 (en
Inventor
Noriyasu Iwane
典靖 岩根
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.)
Furukawa Electric Co Ltd
Original Assignee
Furukawa Electric Co Ltd
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 Furukawa Electric Co Ltd filed Critical Furukawa Electric Co Ltd
Priority to JP2006279478A priority Critical patent/JP4785056B2/en
Publication of JP2008096328A publication Critical patent/JP2008096328A/en
Application granted granted Critical
Publication of JP4785056B2 publication Critical patent/JP4785056B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10Energy storage using batteries

Landscapes

  • Secondary Cells (AREA)
  • Tests Of Electric Status Of Batteries (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a charging rate estimation method, a charging rate estimation device, and a secondary battery power supply system, capable of stably determining the optimum values of adjustment parameters of an approximate expression that allows determination of high-accuracy approximation for the open-circuit voltage of a secondary battery, in a short time, and for the estimation of the charging rate. <P>SOLUTION: In step S104, a voltage drop amount ΔV in data acquisition period is calculated. In step S105, coefficient β is calculated from the ΔV, calculated in step S104 and Ci input from a storage part 140. In step S106, initial value Ai0 of an adjustment parameter Ai is calculated by using the coefficient β calculated in the step S105, and initial value Vc0 of an adjustment parameter Vc is calculated. By using the initial values Ai0, Vc0 of the adjustment parameters thus determined and Bi0 read from the storage part 140, advance is made to next step S107, where the adjustment parameters are optimized by the method of least squares. <P>COPYRIGHT: (C)2008,JPO&INPIT

Description

本発明は、負荷に電力を供給する二次電池の充電率を推定する充電率推定方法及び充電率推定装置の技術分野に関するものである。   The present invention relates to a technical field of a charging rate estimation method and a charging rate estimation device that estimate a charging rate of a secondary battery that supplies power to a load.

従来から、自動車等に搭載される鉛蓄電池等の二次電池に関し、残存する充電率を正確に知ることが要請されている。一般に、二次電池においては充電率と開回路電圧との間に相関があるため、開回路電圧を求めることにより充電率を推定することができる。   2. Description of the Related Art Conventionally, it is required to accurately know the remaining charge rate for secondary batteries such as lead storage batteries mounted on automobiles and the like. Generally, in a secondary battery, since there is a correlation between the charging rate and the open circuit voltage, the charging rate can be estimated by obtaining the open circuit voltage.

しかし、二次電池の開回路電圧は充電又は放電を行っていない状態で行う必要がある上に、充電又は放電の終了後に開回路電圧が安定するまでには長い時間を要している。そのため、所定の条件下で二次電池の開回路電圧を短時間内に測定し、この電圧測定値を用いて開回路電圧の時間特性を近似する関数を決定し、決定された関数に基づき開回路電圧の収束値を求める方法が種々提案されている(例えば、特許文献1 〜3)。   However, the open circuit voltage of the secondary battery needs to be performed in a state where charging or discharging is not performed, and it takes a long time for the open circuit voltage to become stable after completion of charging or discharging. Therefore, the open circuit voltage of the secondary battery is measured within a short time under a predetermined condition, and a function that approximates the time characteristic of the open circuit voltage is determined using the measured voltage value, and the open circuit voltage is opened based on the determined function. Various methods for obtaining the convergence value of the circuit voltage have been proposed (for example, Patent Documents 1 to 3).

上記従来の方法で二次電池の開回路電圧の収束値を求める場合、その精度は近似計算に用いる関数の精度に依存する。例えば、多項式関数や対数関数などの時間特性を持つ関数に基づき開回路電圧の収束値を計算する方法が一般的である。しかし、これらの関数は、二次電池の開回路電圧の時間特性を高い精度で近似することは困難であり、求めた開回路電圧の収束値の誤差が大きくなってしまう。   When the convergence value of the open circuit voltage of the secondary battery is obtained by the conventional method, the accuracy depends on the accuracy of the function used for the approximate calculation. For example, a general method is to calculate the convergence value of the open circuit voltage based on a function having time characteristics such as a polynomial function or a logarithmic function. However, it is difficult for these functions to approximate the time characteristic of the open circuit voltage of the secondary battery with high accuracy, and an error in the convergence value of the obtained open circuit voltage becomes large.

これに対し特許文献4では、二次電池の開回路電圧の近似計算に4次以上の指数減衰関数を用いており、これにより開回路電圧を高い精度で近似できるようにし、充電率を高精度に推定できるようにしている。
特開平7−98367号 特開2002−234408号 特開2003−75518号 特開2005−43339号
On the other hand, in Patent Document 4, an exponential decay function of the fourth order or higher is used for the approximate calculation of the open circuit voltage of the secondary battery, so that the open circuit voltage can be approximated with high accuracy, and the charging rate is highly accurate. Can be estimated.
JP-A-7-98367 JP 2002-234408 A JP 2003-75518 A JP 2005-43339 A

しかしながら、上記従来の二次電池の充電率推定方法では、次のような問題があった。開回路電圧の近似計算に用いる近似式は、開回路電圧の推定値の精度を高めるために調整パラメータを有しており、所定の条件下で二次電池の開回路電圧を測定して取得した電圧測定値を用いてこの調整パラメータを最適に決定していた。   However, the conventional secondary battery charging rate estimation method has the following problems. The approximate expression used for the approximate calculation of the open circuit voltage has an adjustment parameter to improve the accuracy of the estimated value of the open circuit voltage, and is obtained by measuring the open circuit voltage of the secondary battery under a predetermined condition. This tuning parameter was optimally determined using voltage measurements.

しかし、二次電池の開回路電圧の近似計算に4次以上の指数減衰関数のような複雑な非線形関数を用いる場合には、電圧測定値を用いて調整パラメータを最適化する際、最小二乗法の最適解をガウス-ニュートン法やレーベンベルグ-マルカート法のような逐次演算を解くか、カルマン-フィルタ演算のようなフィルタリング演算を行なう必要があり、何れの場合でも調整パラメータの初期値を設定する必要があり、その初期値によっては、計算が発散して調整パラメータの最適解値が得られない可能性や、得られたとしても非常に長い時間がかかってしまうといった問題があった。   However, when a complex nonlinear function such as a fourth or higher order exponential decay function is used for the approximate calculation of the open circuit voltage of the secondary battery, the least square method is used when optimizing the adjustment parameter using the voltage measurement value. It is necessary to solve the sequential solution such as Gauss-Newton method or Levenberg-Marquardt method or the filtering operation such as Kalman-filter operation. Depending on the initial value, there is a possibility that calculation may diverge and the optimum solution value of the adjustment parameter may not be obtained, and even if it is obtained, it takes a very long time.

そこで、本発明はこれらの問題を解決するためになされたものであり、二次電池の開回路電圧を高精度で近似できる近似式の調整パラメータの最適値を短時間でかつ安定的に決定して充電率を推定することができる充電率推定方法、充電率推定装置及び二次電池電源システムを提供することを目的とする。   Therefore, the present invention has been made to solve these problems, and the optimum value of the adjustment parameter of the approximate expression that can approximate the open circuit voltage of the secondary battery with high accuracy is determined in a short time and stably. It is an object of the present invention to provide a charging rate estimation method, a charging rate estimation device, and a secondary battery power supply system that can estimate a charging rate.

この発明の充電率推定方法の第1の態様は、充放電終了後の二次電池の開回路電圧(OCV;Open Circuit Voltage)の時間変化を、調整パラメータを有する電圧特性式で近似し、前記電圧特性式から前記二次電池安定時の開回路電圧を算出して前記二次電池の充電率を推定する充電率推定方法であって、前記二次電池の充放電終了後の所定のデータ取得期間に、前記二次電池の電圧測定値を複数個取得し、前記データ取得期間における前記電圧測定値の変化量に依存する所定の関係式及び事前に設定された所定の算出方法で前記調整パラメータの初期値を設定し、前記調整パラメータを前記初期値から漸次修正しながら前記電圧特性式で算出した前記開回路電圧と前記電圧測定値との誤差が最小となるときの前記調整パラメータの最適値を推定することを特徴とする。   According to a first aspect of the charging rate estimation method of the present invention, the time variation of the open circuit voltage (OCV) of the secondary battery after the end of charge / discharge is approximated by a voltage characteristic equation having an adjustment parameter, A charging rate estimation method for calculating an open circuit voltage when the secondary battery is stable from a voltage characteristic equation to estimate a charging rate of the secondary battery, and obtaining predetermined data after completion of charging and discharging of the secondary battery In the period, a plurality of voltage measurement values of the secondary battery are acquired, and the adjustment parameter is calculated using a predetermined relational expression that depends on a change amount of the voltage measurement value in the data acquisition period and a predetermined calculation method set in advance. The adjustment parameter is optimal when the error between the open circuit voltage calculated by the voltage characteristic equation and the voltage measurement value is minimized while the adjustment parameter is gradually corrected from the initial value. Estimate It is characterized in.

この発明の充電率推定方法の他の態様は、前記調整パラメータをAi(i=1〜n)、Bi(i=1〜n)及びVcとし、nを4以上の整数としたとき、前記電圧特定式が、
[数9]
V(t)=A1・exp(B1・t))+A2・exp(B2・t)+・・・
+An・exp(Bn・t)+Vc
で表され、前記データ取得期間初期の前記電圧測定値をVa、前記データ取得期間終了直前の前記電圧測定値をVbとしたとき、事前に設定された定数α及びCi(i=1〜n)を用いて
[数10]
(C1+C2+・・・+Cn)・β=(Va−Vb)・α
を満たすようパラメータβを決定し、
[数11]
Ai0=Ci・β (i=1〜n)
で算出されるAi0を前記調整パラメータAiの前記初期値とし、
[数12]
Vc0=Va−(Va−Vb)・α
で算出されるVc0を前記調整パラメータVcの前記初期値とすることを特徴とする。
In another aspect of the charging rate estimation method of the present invention, when the adjustment parameters are Ai (i = 1 to n), Bi (i = 1 to n) and Vc, and n is an integer of 4 or more, the voltage The specific formula is
[Equation 9]
V (t) = A1 · exp (B1 · t)) + A2 · exp (B2 · t) +.
+ An · exp (Bn · t) + Vc
And the constants α and Ci set in advance (i = 1 to n), where Va is the voltage measurement value at the beginning of the data acquisition period, and Vb is the voltage measurement value immediately before the end of the data acquisition period. Using
[Equation 10]
(C1 + C2 +... + Cn) · β = (Va−Vb) · α
Parameter β to satisfy
[Equation 11]
Ai0 = Ci · β (i = 1 to n)
Ai0 calculated in step S is used as the initial value of the adjustment parameter Ai,
[Equation 12]
Vc0 = Va− (Va−Vb) · α
Vc0 calculated by the above is used as the initial value of the adjustment parameter Vc.

この発明の充電率推定方法の他の態様は、前記調整パラメータBi(i=1〜n)の初期値が、事前に設定された定数であることを特徴とする。   Another aspect of the charging rate estimation method of the present invention is characterized in that an initial value of the adjustment parameter Bi (i = 1 to n) is a constant set in advance.

この発明の充電率推定方法の他の態様は、前記調整パラメータBi(i=1〜n)の初期値を前記二次電池の温度から算出する温度関数が事前に設定され、前記データ取得期間に前記温度を測定して温度測定値を取得し、前記温度測定値を前記温度関数に代入して前記調整パラメータBi(i=1〜n)の初期値を決定することを特徴とする。   In another aspect of the charging rate estimation method of the present invention, a temperature function for calculating an initial value of the adjustment parameter Bi (i = 1 to n) from a temperature of the secondary battery is set in advance, and the data acquisition period A temperature measurement value is obtained by measuring the temperature, and the initial value of the adjustment parameter Bi (i = 1 to n) is determined by substituting the temperature measurement value into the temperature function.

この発明の充電率推定方法の他の態様は、前記Va又は前記Vbの少なくともいずれか一方が、2以上の前記電圧測定値を平均して算出した平均値であることを特徴とする。   Another aspect of the charging rate estimation method of the present invention is characterized in that at least one of the Va and the Vb is an average value calculated by averaging two or more voltage measurement values.

この発明の充電率推定方法の他の態様は、前記二次電池安定時の前記開回路電圧が、前記調整パラメータの前記最適値を決定したときのVcで与えられることを特徴とする。   Another aspect of the charging rate estimation method of the present invention is characterized in that the open circuit voltage when the secondary battery is stable is given by Vc when the optimum value of the adjustment parameter is determined.

この発明の充電率推定方法の他の態様は、前記二次電池安定時の前記開回路電圧が、前記調整パラメータの前記最適値を決定した関数式に予め定めた所定時間を代入して算出することを特徴とする。   According to another aspect of the charging rate estimation method of the present invention, the open circuit voltage when the secondary battery is stable is calculated by substituting a predetermined time for a function equation that determines the optimum value of the adjustment parameter. It is characterized by that.

この発明の充電率推定装置の第1の態様は、二次電池の充電率を推定する充電率推定装置であって、前記二次電池の電圧を測定する電圧センサと、調整パラメータを有する電圧特性式を用いて前記充電率を推定するための演算を実行制御する制御部と、前記電圧センサで測定された電圧測定値と、前記制御部による前記演算に用いられるデータとを記憶する記憶部と、を備え、前記制御部が、前記二次電池の充放電終了後の所定のデータ取得期間に、前記電圧センサから前記二次電池の電圧測定値を複数個取得し、前記データ取得期間における前記電圧測定値の変化量に依存する所定の関係式及び事前に設定された所定の算出方法で前記調整パラメータの初期値を設定し、前記調整パラメータを前記初期値から漸次修正しながら前記電圧特性式で算出した前記開回路電圧と前記電圧測定値との誤差が最小となるときの前記調整パラメータの最適値を推定し、前記調整パラメータの最適値を用いた前記電圧特性式から前記二次電池安定時の前記開回路電圧を算出して前記充電率を推定することを特徴とする。   A first aspect of the charging rate estimation device according to the present invention is a charging rate estimation device that estimates a charging rate of a secondary battery, a voltage sensor that measures the voltage of the secondary battery, and a voltage characteristic having an adjustment parameter. A control unit that executes and controls calculation for estimating the charging rate using an equation; a storage unit that stores a voltage measurement value measured by the voltage sensor; and data used for the calculation by the control unit; The control unit acquires a plurality of voltage measurement values of the secondary battery from the voltage sensor in a predetermined data acquisition period after the end of charging and discharging of the secondary battery, and the control unit in the data acquisition period An initial value of the adjustment parameter is set by a predetermined relational expression that depends on a change amount of the voltage measurement value and a predetermined calculation method set in advance, and the voltage characteristic is gradually corrected from the initial value. Estimating the optimum value of the adjustment parameter when the error between the open circuit voltage calculated in step 1 and the voltage measurement value is minimum, and stabilizing the secondary battery from the voltage characteristic equation using the optimum value of the adjustment parameter The open circuit voltage at the time is calculated to estimate the charging rate.

この発明の充電率推定装置の第2の態様は、前記制御部が、前記調整パラメータをAi、Bi(i=1〜n)及びVcとし、nを4以上の整数としたとき、次式
[数13]
V(t)=A1・exp(B1・t))+A2・exp(B2・t)+・・・
+An・exp(Bn・t)+Vc
を前記電圧特定式として用い、前記データ取得期間初期の前記電圧測定値をVa、前記データ取得期間終了直前の前記電圧測定値をVbとしたとき、事前に設定された定数α及びCi(i=1〜n)を用いて
[数14]
(C1+C2+・・・+Cn)・β=(Va−Vb)・α
を満たすようパラメータβを決定し、
[数15]
Ai0=Ci・β (i=1〜n)
で算出されるAi0を前記調整パラメータAiの前記初期値とし、
[数16]
Vc0=Va−(Va−Vb)・α
で算出されるVc0を前記調整パラメータVcの前記初期値として用いていることを特徴とする。
According to a second aspect of the charging rate estimation device of the present invention, when the control unit sets the adjustment parameters to Ai, Bi (i = 1 to n) and Vc, and n is an integer of 4 or more,
[Equation 13]
V (t) = A1 · exp (B1 · t)) + A2 · exp (B2 · t) +.
+ An · exp (Bn · t) + Vc
Is used as the voltage specifying formula, Va is the voltage measurement value at the initial stage of the data acquisition period, and Vb is the voltage measurement value immediately before the end of the data acquisition period. 1-n)
[Equation 14]
(C1 + C2 +... + Cn) · β = (Va−Vb) · α
Parameter β to satisfy
[Equation 15]
Ai0 = Ci · β (i = 1 to n)
Ai0 calculated in step S is used as the initial value of the adjustment parameter Ai,
[Equation 16]
Vc0 = Va− (Va−Vb) · α
Is used as the initial value of the adjustment parameter Vc.

この発明の二次電池電源システムの第1の態様は、第1の態様又は第2の態様に記載の充電率推定装置を備えることを特徴する。   The 1st aspect of the secondary battery power supply system of this invention is provided with the charging rate estimation apparatus as described in a 1st aspect or a 2nd aspect.

本発明によれば、二次電池の開回路電圧の時間変化を、調整パラメータを有する電圧特性式で近似し、該調整パラメータの初期値を好適に設定するようにしたことにより、該調整パラメータの最適値を短時間でかつ安定的に決定することが可能となる。そして、最適化された調整パラメータを用いて二次電池安定時の開回路電圧を高精度で推定でき、これにより充電率を高精度に推定することができる充電率推定方法、充電率推定装置及び二次電池電源システムを提供することが可能となる。   According to the present invention, the time variation of the open circuit voltage of the secondary battery is approximated by the voltage characteristic equation having the adjustment parameter, and the initial value of the adjustment parameter is suitably set, so that the adjustment parameter The optimum value can be determined stably in a short time. Then, the open circuit voltage when the secondary battery is stable can be estimated with high accuracy using the optimized adjustment parameter, and thereby the charge rate estimation method, the charge rate estimation device, and the charge rate can be estimated with high accuracy. A secondary battery power supply system can be provided.

図面を参照して本発明の好ましい実施の形態における充電率推定方法、充電率推定装置及び二次電池電源システムの構成について詳細に説明する。なお、同一機能を有する各構成部については、図示及び説明簡略化のため、同一符号を付して示す。以下では、自動車等の車両に搭載される二次電池に対して本発明を適用する場合を説明する。   Configurations of a charging rate estimation method, a charging rate estimation device, and a secondary battery power supply system according to a preferred embodiment of the present invention will be described in detail with reference to the drawings. In addition, about each structural part which has the same function, the same code | symbol is attached | subjected and shown for simplification of illustration and description. Below, the case where this invention is applied with respect to the secondary battery mounted in vehicles, such as a motor vehicle, is demonstrated.

本発明の実施の形態に係る二次電池電源システムを、図2を用いて説明する。図2は、本実施形態に係る二次電池電源システム100の概略構成を示すブロック図である。二次電池電源システム100は、二次電池110と、二次電池110の電圧を測定する電圧センサ120と、電圧センサ120から電圧測定値を入力して二次電池110を制御している制御部130と、制御部130で行われる制御に用いられるデータ等を記憶する記憶部140と、制御部130からの制御指令に基づいて二次電池110を充電するための充電回路150とを備えている。   A secondary battery power supply system according to an embodiment of the present invention will be described with reference to FIG. FIG. 2 is a block diagram showing a schematic configuration of the secondary battery power supply system 100 according to the present embodiment. The secondary battery power supply system 100 includes a secondary battery 110, a voltage sensor 120 that measures the voltage of the secondary battery 110, and a control unit that controls the secondary battery 110 by inputting a voltage measurement value from the voltage sensor 120. 130, a storage unit 140 that stores data used for control performed by the control unit 130, and a charging circuit 150 for charging the secondary battery 110 based on a control command from the control unit 130. .

上記のように構成された二次電池電源システム100には、モータ等の負荷10が接続されており、二次電池110から負荷10に電力が供給される構成となっている。自動車等の車両の場合には、二次電池110として例えば鉛蓄電池が用いられている。   The secondary battery power supply system 100 configured as described above is connected to a load 10 such as a motor, and power is supplied from the secondary battery 110 to the load 10. In the case of a vehicle such as an automobile, for example, a lead storage battery is used as the secondary battery 110.

制御部130は、演算処理を行うCPU等で構成されており、二次電池電源システム100の全体の動作を制御するとともに、所定のタイミングで後述の充電率推定のための演算処理を実行している。制御部130に接続された記憶140は、制御プログラム等の各種プログラムを事前に記憶させておくためのROMや、制御部130の処理に用いられる各種データを記憶させるためのRAM等を備えている。   The control unit 130 is configured by a CPU or the like that performs arithmetic processing, controls the overall operation of the secondary battery power supply system 100, and executes arithmetic processing for charging rate estimation described later at a predetermined timing. Yes. The storage 140 connected to the control unit 130 includes a ROM for storing various programs such as a control program in advance, a RAM for storing various data used for processing of the control unit 130, and the like. .

本発明の充電率推定装置は、少なくとも後述の充電率推定のための演算処理を実行する制御部と、記憶部140及び電圧センサ120とで構成することができる。図2に示す実施形態では、充電率推定の演算処理を行う制御部を、二次電池電源システム100の全体の動作を制御している制御部130に含める構成としている。   The charging rate estimation apparatus according to the present invention can be configured with at least a control unit that executes arithmetic processing for charging rate estimation described later, a storage unit 140, and a voltage sensor 120. In the embodiment shown in FIG. 2, a control unit that performs calculation processing for charging rate estimation is included in the control unit 130 that controls the overall operation of the secondary battery power supply system 100.

次に、本実施形態に係る二次電池電源システム100で用いられる二次電池110の充電率推定方法について説明する。一般に、二次電池の充電率は、二次電池が安定しているときの開回路電圧と強い相関関係があることから、安定時の開回路電圧が求まるとこれから二次電池の充電率を推定することが可能となる。   Next, a charging rate estimation method for the secondary battery 110 used in the secondary battery power supply system 100 according to the present embodiment will be described. In general, the charging rate of the secondary battery has a strong correlation with the open circuit voltage when the secondary battery is stable, so when the open circuit voltage at the stable time is obtained, the charging rate of the secondary battery is estimated from this It becomes possible to do.

しかし、特に車両に搭載された二次電池では、充放電が頻繁に繰り返されているため、電圧が安定していないのが通常の状態である。充放電が行われると二次電池の内部で分極が発生し、これが充放電終了後に徐々に解消されていく。この分極が解消されて二次電池が安定するまでに、通常十数時間から数日という極めて長い時間を要する。そのため、二次電池が十分に安定する前に再び充放電されることが多く、二次電池の安定時の開回路電圧を測定するのが極めて困難であった。   However, particularly in a secondary battery mounted on a vehicle, since charging and discharging are frequently repeated, it is normal that the voltage is not stable. When charging / discharging is performed, polarization occurs inside the secondary battery, and this is gradually eliminated after completion of charging / discharging. It usually takes an extremely long time of several tens of hours to several days until the polarization is eliminated and the secondary battery is stabilized. Therefore, the secondary battery is often charged and discharged again before it becomes sufficiently stable, and it is extremely difficult to measure the open circuit voltage when the secondary battery is stable.

本実施形態では、制御部130において開回路電圧の時間的な変動を高い精度で近似できる電圧特性式を用いることで、二次電池110の安定時の開回路電圧を高精度で推定できるようにしている。電圧特性式として、4次以上の指数減衰関数を含む近似式を用いることで、開回路電圧の時間変動を高精度に推定することができる。本実施形態では、一実施例として下記の4次の指数減衰関数を含む近似式を用いている。
[数17]
V(t)=A1・exp(B1・t))+A2・exp(B2・t)+
A3・exp(B3・t))+A4・exp(B4・t)+Vc (式1)
In the present embodiment, the control unit 130 uses a voltage characteristic equation that can approximate the temporal variation of the open circuit voltage with high accuracy, so that the open circuit voltage when the secondary battery 110 is stable can be estimated with high accuracy. ing. By using an approximate expression including a fourth-order or higher exponential decay function as the voltage characteristic expression, it is possible to estimate the time fluctuation of the open circuit voltage with high accuracy. In this embodiment, an approximate expression including the following fourth-order exponential decay function is used as an example.
[Equation 17]
V (t) = A1 · exp (B1 · t)) + A2 · exp (B2 · t) +
A3 * exp (B3 * t)) + A4 * exp (B4 * t) + Vc (Formula 1)

上記(式1)において、A1〜A4、B1〜B4(以下ではそれぞれをAi、Biと示す)、及びVcは、V(t)が開回路電圧を高精度に推定できるよう調整するための調整パラメータ(フィッティングパラメータ)である。これらの調整パラメータは、所定の時刻tにおける電圧測定値とV(t)との誤差が最小となるように調整される。   In the above (Formula 1), A1 to A4, B1 to B4 (hereinafter referred to as Ai and Bi) and Vc are adjustments for adjusting V (t) so that the open circuit voltage can be estimated with high accuracy. It is a parameter (fitting parameter). These adjustment parameters are adjusted so that the error between the voltage measurement value at a predetermined time t and V (t) is minimized.

二次電池110の充放電終了後の開回路電圧の時間変動の一実施例として、二次電池110を充電した後60分間にわたる開回路電圧の時間変動を測定した例を図3に示す。同図では、充電終了後10分間の電圧測定値を●印で示し、それ以降の電圧測定値を○印で示している。   FIG. 3 shows an example of measuring the time fluctuation of the open circuit voltage over 60 minutes after charging the secondary battery 110 as an example of the time fluctuation of the open circuit voltage after the end of charging and discharging of the secondary battery 110. In the figure, the voltage measurement values for 10 minutes after the end of charging are indicated by ● and the voltage measurement values thereafter are indicated by ○.

充電終了後10分間の電圧測定値(●印)を、(式1)に示す4次の指数減衰関数を有する電圧特性式でフィッティングしたときと、対数関数を用いてフィッティングしたときの結果を図3に併せて示している。それぞれを最適にフィッティングした結果は、下記の通りである。
[数18]
V(t)=1.80933exp(-t/4.65331)+0.2895exp(-t/0.41691)
−0.90055exp(-t/4.65129)+0.9exp(-t/0.004)+13.35703 (式2)
V(t)=-0.2517ln(t)+14.072 (式3)
Fig. 2 shows the results of fitting the voltage measurement value (● mark) for 10 minutes after the end of charging with the voltage characteristic equation having the fourth-order exponential decay function shown in (Equation 1) and fitting using the logarithmic function. 3 is also shown. The results of optimal fitting of each are as follows.
[Equation 18]
V (t) = 1.80933exp (-t / 4.65331) + 0.2895exp (-t / 0.41691)
-0.90055exp (-t / 4.65129) + 0.9exp (-t / 0.004) +13.35703 (Formula 2)
V (t) =-0.2517ln (t) +14.072 (Formula 3)

上記のフィッティング式と充電終了後10分経過以降の電圧測定値(○印)とを比較すると明らかなように、対数関数を用いてフィッティングした(式3)では時間の経過とともに誤差が大きくなっていくのに対し、4次の指数減衰関数を有する電圧特性式でフィッティングした(式2)の場合には、充電終了後10分経過以降も電圧測定値と極めてよい一致を示しており、開回路電圧の時間変動を高精度に推定できていることがわかる。   As is obvious from comparing the above fitting equation with the measured voltage value (marked with ○) after 10 minutes from the end of charging, the error is increased with the passage of time in the case of fitting using the logarithmic function (Equation 3). On the other hand, in the case of fitting with a voltage characteristic equation having a fourth-order exponential decay function (Equation 2), even after 10 minutes from the end of charging, it shows a very good agreement with the measured voltage value. It can be seen that the time variation of the voltage can be estimated with high accuracy.

上記のように、(式1)に示す4次の指数減衰関数を有する電圧特性式を用いて二次電池110の開回路電圧を高精度に推定できるようにするためには、充放電終了後に二次電池110の開回路電圧を所定のデータ取得期間測定し、この電圧測定値を用いて(式1)を最適にフィッティングする必要がある。また、制御部130で充電率を推定して制御等に直ちに利用できるよう、上記フィッティングを計算が発散しないように、また短時間で行えるようにする必要がある。   As described above, in order to be able to estimate the open circuit voltage of the secondary battery 110 with high accuracy using the voltage characteristic equation having the fourth-order exponential decay function shown in (Equation 1), It is necessary to measure the open circuit voltage of the secondary battery 110 for a predetermined data acquisition period and optimally fit (Equation 1) using this voltage measurement value. In addition, the fitting needs to be performed in a short time so that the calculation does not diverge so that the controller 130 can estimate the charging rate and immediately use it for control or the like.

(式1)の電圧特性式を電圧測定値に短時間で高精度にフィッティングさせるためには、(式1)に用いられている調整パラメータの初期値を適切に選択して用いることが重要である。中でも、調整パラメータA1〜A4及びVcは、二次電池110の充電状態や劣化度等によって変化することから、調整パラメータA1〜A4及びVcの初期値を適切に設定することが、フィッティングを高速にかつ安定的に行う上で効果が大きい。以下では、本実施形態で用いる充電率推定方法において、調整パラメータの初期値を好適に設定している方法を、図4を用いて説明する。   In order to fit the voltage characteristic equation of (Equation 1) to the voltage measurement value with high accuracy in a short time, it is important to appropriately select and use the initial values of the adjustment parameters used in (Equation 1). is there. Among them, the adjustment parameters A1 to A4 and Vc change depending on the state of charge of the secondary battery 110, the degree of deterioration, and the like. Therefore, appropriately setting the initial values of the adjustment parameters A1 to A4 and Vc speeds up the fitting. In addition, the effect is large in performing stably. Hereinafter, a method of suitably setting the initial value of the adjustment parameter in the charging rate estimation method used in the present embodiment will be described with reference to FIG.

図4は、(式1)の調整パラメータの初期値を設定する方法を説明するための二次電池110の開回路電圧の変化の例を示す図である。同図において、○印は電圧センサ120から取得した電圧測定値21を示しており、二次電池110の充放電終了後のデータ取得期間に取得されたものである。データ取得期間初期の電圧をVa、データ取得期間終了直前の電圧をVbとし、その間の電圧変化をΔVとしたとき、ΔVは次式のように表わされる。
[数19]
ΔV=Va−Vb (式4)
FIG. 4 is a diagram illustrating an example of a change in the open circuit voltage of the secondary battery 110 for explaining a method of setting the initial value of the adjustment parameter of (Equation 1). In the figure, a circle indicates a voltage measurement value 21 acquired from the voltage sensor 120, which is acquired during a data acquisition period after the end of charging / discharging of the secondary battery 110. When the voltage at the beginning of the data acquisition period is Va, the voltage immediately before the end of the data acquisition period is Vb, and the voltage change during that period is ΔV, ΔV is expressed as the following equation.
[Equation 19]
ΔV = Va−Vb (Formula 4)

なお、Vaは、データ取得期間初期の複数の電圧測定値を平均したものであってもよく、例えば図4に示す21aの電圧測定値を平均してVaとしてよい。同様に、Vbは、データ取得期間終了直前の複数の電圧測定値を平均したものであってもよく、例えば図4に示す21bの電圧測定値を平均してVbとしてよい。   Va may be an average of a plurality of voltage measurement values at the beginning of the data acquisition period. For example, Va may be obtained by averaging the voltage measurement values of 21a shown in FIG. Similarly, Vb may be an average of a plurality of voltage measurement values immediately before the end of the data acquisition period. For example, the voltage measurement value of 21b shown in FIG.

上記のΔVを用いて、(式1)における定数項Vcの初期値(Vc0とする)を次式用いて設定している。
[数20]
Vc0=Va−α・ΔV (式5)
ここで、αは事前に実験等を行って設定した定数である。二次電池110が安定状態にあるときの開回路電圧は、(式1)において時間tに依存する指数減衰関数が0となったときのVcに一致する。よって、(式5)は、安定時の開回路電圧の初期値を設定する式となる。
Using the above ΔV, the initial value (constant Vc0) of the constant term Vc in (Expression 1) is set using the following expression.
[Equation 20]
Vc0 = Va−α · ΔV (Formula 5)
Here, α is a constant set through experiments and the like in advance. The open circuit voltage when the secondary battery 110 is in a stable state coincides with Vc when the exponential decay function depending on the time t becomes 0 in (Equation 1). Therefore, (Expression 5) is an expression for setting an initial value of the open circuit voltage at the time of stabilization.

また、調整パラメータA1〜A4に対しては、それぞれの初期値の基準となる値(以下では基準値という)を事前に決定して記憶部140に記憶させている。すなわち、Ai(i=1〜n)の初期値をAi0とし、上記の基準値をCiとしたとき、
[数21]
Ai0=Ci・β (式6)
によりAiの初期値Ai0を決定するようにしている。基準値をCiは、調整パラメータAiの初期値としてその比例関係を事前に決めておくようにしたものである。
In addition, for the adjustment parameters A1 to A4, values serving as references for the initial values (hereinafter referred to as reference values) are determined in advance and stored in the storage unit 140. That is, when the initial value of Ai (i = 1 to n) is Ai0 and the reference value is Ci,
[Equation 21]
Ai0 = Ci · β (Formula 6)
Thus, the initial value Ai0 of Ai is determined. The reference value Ci is an initial value of the adjustment parameter Ai and its proportional relationship is determined in advance.

以下では、(式6)の係数βを決定する方法を説明する。(式1)において、時間t=0としたとき次式が成り立つ。
[数22]
V(0)=A1+A2+A3+A4+Vc (式7)
ここで、V(o)はデータ取得開始時の開回路電圧であるから、
[数23]
V(0)=Va (式8)
となる。
Hereinafter, a method of determining the coefficient β in (Expression 6) will be described. In (Expression 1), the following expression is established when time t = 0.
[Equation 22]
V (0) = A1 + A2 + A3 + A4 + Vc (Formula 7)
Here, V (o) is an open circuit voltage at the start of data acquisition.
[Equation 23]
V (0) = Va (Formula 8)
It becomes.

(式7)において、Ai及びVcとしてそれぞれの初期値を用い、さらに(式5)及び(式8)を用いることで、
[数24]
Va=A10+A20+A30+A40+Vc0
=A10+A20+A30+A40+(Va−α・ΔV) (式9)
となる。上式より、
[数25]
A10+A20+A30+A40=α・ΔV (式10)
が成り立つ。
In (Expression 7), by using respective initial values as Ai and Vc, and further using (Expression 5) and (Expression 8),
[Equation 24]
Va = A10 + A20 + A30 + A40 + Vc0
= A10 + A20 + A30 + A40 + (Va-α · ΔV) (Formula 9)
It becomes. From the above formula,
[Equation 25]
A10 + A20 + A30 + A40 = α · ΔV (Formula 10)
Holds.

さらに、(式6)と(式10)とから
[数26]
(C1+C2+C3+C4)・β=α・ΔV (式11)
となり、これから係数βは次式で決定される。
[数27]
β=α・ΔV/(C1+C2+C3+C4) (式12)
(式12)において、ΔVは電圧測定値から算出され、α及びCiは事前に決定されて記憶部140に記憶されているデータである。よって、(式12)よりβが決定され、(式6)よりAiの初期値Ai0が設定される。
Furthermore, from (Equation 6) and (Equation 10)
[Equation 26]
(C1 + C2 + C3 + C4) · β = α · ΔV (Formula 11)
From this, the coefficient β is determined by the following equation.
[Equation 27]
β = α · ΔV / (C1 + C2 + C3 + C4) (Formula 12)
In (Expression 12), ΔV is calculated from the voltage measurement value, and α and Ci are data determined in advance and stored in the storage unit 140. Therefore, β is determined from (Expression 12), and the initial value Ai0 of Ai is set from (Expression 6).

調整パラメータAiの初期値Ai0及びVcの初期値Vc0が上記のようにして決定されるのに対し、別の調整パラメータであるBi(i=1〜n)は、初期値を事前に決定して記憶部140に記憶させておき、これを読み込んで用いることができる。 The initial value Ai0 of the adjustment parameter Ai and the initial value Vc0 of Vc are determined as described above, whereas Bi (i = 1 to n), which is another adjustment parameter, determines the initial value in advance. The data can be stored in the storage unit 140 and read and used.

次に、本実施形態の制御部130において実行される二次電池110の充電率推定方法について、その処理の流れを図1に示す流れ図を用いて説明する。図1に示す演算処理は、主に制御部130が記憶部140に保持される制御プログラムに基づいて実行する処理であり、二次電池電源システム100における二次電池110の充放電終了後に所定のタイミングで実行開始される。   Next, the process flow estimation method of the secondary battery 110 executed in the control unit 130 of the present embodiment will be described with reference to the flowchart shown in FIG. The calculation process shown in FIG. 1 is a process that is mainly executed by the control unit 130 based on a control program held in the storage unit 140, and has a predetermined value after the secondary battery 110 is charged and discharged in the secondary battery power supply system 100. Execution starts at the timing.

制御部130において、充放電終了後の所定のタイミングで充電率の推定処理が要求されると、まずステップS101において、演算処理に必要なパラメータの設定が行われる。演算処理に必要なパラメータとしては、電圧センサ120から電圧測定値を取得する際のサンプリング間隔ΔTs、サンプル取得数Ns、及び(式1)の電圧特性式の調整パラメータの初期値を設定するための初期化パラメータがある。サンプリング間隔ΔTs及びサンプル取得数Nsは、事前に決定された値を用いればよく、例えばΔTs=10(秒)、Ns=60(個)等の値を用いることができる。   When the control unit 130 requests charging rate estimation processing at a predetermined timing after completion of charging / discharging, first, in step S101, parameters necessary for calculation processing are set. As parameters necessary for the arithmetic processing, the sampling interval ΔTs when acquiring the voltage measurement value from the voltage sensor 120, the number Ns of sample acquisitions, and the initial value of the adjustment parameter of the voltage characteristic equation of (Equation 1) are set. There are initialization parameters. As the sampling interval ΔTs and the sample acquisition number Ns, values determined in advance may be used. For example, values such as ΔTs = 10 (seconds) and Ns = 60 (pieces) may be used.

また、初期化パラメータとしては、(式5)で用いたパラメータα、(式6)で用いた基準値Ci、及び調整パラメータBiの初期値Bi0がある。パラメータαは、実験等によって最適な値を事前に決定して用いることができる。基準値Ciについても、調整パラメータA1〜A4の相互の最適な比例関係を事前に評価し、この比例関係に対応した値をCiとすることができる。同様にして調整パラメータBiの初期値Bi0も、事前に好適な値を決定しておくことができる。これらは、事前に記憶部140に記憶させておく。   The initialization parameters include the parameter α used in (Expression 5), the reference value Ci used in (Expression 6), and the initial value Bi0 of the adjustment parameter Bi. As the parameter α, an optimum value can be determined in advance by experiments or the like. As for the reference value Ci, the optimum proportional relationship between the adjustment parameters A1 to A4 can be evaluated in advance, and a value corresponding to this proportional relationship can be set as Ci. Similarly, a suitable value can be determined in advance for the initial value Bi0 of the adjustment parameter Bi. These are stored in the storage unit 140 in advance.

ステップS102では、上記のサンプリング間隔ΔTsでサンプル取得数Nsに達するまで電圧センサ120から電圧測定値を取得する。取得した電圧測定値を、ここではVm(j)(j=1〜Ns)とし、これを記憶部140に順次記憶させていく。   In step S102, voltage measurement values are acquired from the voltage sensor 120 until the number of sample acquisitions Ns is reached at the sampling interval ΔTs. Here, the acquired voltage measurement value is Vm (j) (j = 1 to Ns), and this is sequentially stored in the storage unit 140.

ステップS103では、ステップS102におけるデータ取得開始初期の電圧測定値Vaと終了直前の電圧測定値Vbとを算出する。Vaとして取得開始後最初の電圧測定値Vm(1)を設定することができる。あるいは、取得開始初期の複数個(Maとする)の電圧測定値Vm(j)(j=1、Ma)の平均値をVaとしてもよい。同様にして、VbとしてVm(Ns)を設定することができる。あるいは、取得終了直前の複数個の電圧測定値の平均値をVbとしてもよい。   In step S103, the voltage measurement value Va at the beginning of data acquisition in step S102 and the voltage measurement value Vb immediately before the end of data acquisition are calculated. The first voltage measurement value Vm (1) after the start of acquisition can be set as Va. Alternatively, Va may be an average value of a plurality of voltage measurement values Vm (j) (j = 1, Ma) at the beginning of acquisition. Similarly, Vm (Ns) can be set as Vb. Alternatively, an average value of a plurality of voltage measurement values immediately before the end of acquisition may be Vb.

ステップS104では、データ取得期間における電圧低下量ΔVを(式4)に基づいて算出する。また、ステップS105では、ステップS104で算出したΔVと記憶部140から入力したCiとから、(式12)に基づいて係数βを算出する。   In step S104, the voltage drop amount ΔV in the data acquisition period is calculated based on (Equation 4). In step S105, the coefficient β is calculated based on (Equation 12) from ΔV calculated in step S104 and Ci input from the storage unit 140.

ステップS106では、ステップS105で算出した係数βを用いて(式6)に基づいて調整パラメータAiの初期値Ai0を、また(式5)に基づいて調整パラメータVcの初期値Vc0を、それぞれ算出する。以上で決定された調整パラメータの初期値Ai0とVC0、及び記憶部140から読み込んだBi0を用いて、次のステップS107の調整パラメータの最適化に進む。本実施例では最小二乗法による最適化を前提に述べるが、カルマン-フィルタ演算等を用いても全く問題ない。   In step S106, the initial value Ai0 of the adjustment parameter Ai is calculated based on (Equation 6) and the initial value Vc0 of the adjustment parameter Vc is calculated based on (Equation 5) using the coefficient β calculated in step S105. . Using the adjustment parameter initial values Ai0 and VC0 determined as described above and Bi0 read from the storage unit 140, the process proceeds to optimization parameter adjustment in the next step S107. Although the present embodiment is described on the premise of optimization by the least square method, there is no problem even if Kalman-filter operation or the like is used.

ステップS107では、(式1)の電圧特性式とステップS102で電圧センサ120から取得した電圧測定値Vm(j)とで表される次式の誤差分散が最小となるよう、調整パラメータをステップS106で求めた初期値から逐次更新していく。
[数28]
S=Σ[V(j)−Vm(j)] (式13)
In step S107, the adjustment parameter is set to step S106 so that the error variance of the following expression expressed by the voltage characteristic expression of (Expression 1) and the voltage measurement value Vm (j) acquired from the voltage sensor 120 in step S102 is minimized. Update sequentially from the initial value obtained in step 1.
[Equation 28]
S = Σ [V (j) −Vm (j)] 2 (Formula 13)

ステップS108では、ステップS107の最小二乗法によって求まった調整パラメータの最適値をAi、Bi、Vcとし、これを記憶部140に記憶させ、後に電圧特性式(式1)を算出する際にこれを用いる。   In step S108, the optimum values of the adjustment parameters obtained by the least square method in step S107 are Ai, Bi, and Vc, which are stored in the storage unit 140, and this is used when the voltage characteristic equation (Equation 1) is calculated later. Use.

最後のステップS109では、最適化された調整パラメータVcを用いて、二次電池110の充電率を推定する。電圧特性式(式1)を用いて安定時の開回路電圧Voを算出する場合、(式1)においてtを大きくすると指数減衰項は0に収束することから、安定時の開回路電圧Voは(式1)の定数項Vcに等しくなる。このように、二次電池110の電圧特性式として(式1)を用いた場合には、安定時の開回路電圧Voが(式1)の定数項Vcから直ちに知ることができる。
また別の方法として、最適化された調整パラメータを代入した電圧特性式(式1)に予め定めた所定の電池電圧の安定所要時間を代入することによって安定時の開回路電圧Voを求めることもできる。
In the last step S109, the charging rate of the secondary battery 110 is estimated using the optimized adjustment parameter Vc. When calculating the stable open circuit voltage Vo using the voltage characteristic equation (Equation 1), if t is increased in (Equation 1), the exponential decay term converges to 0. Therefore, the stable open circuit voltage Vo is It becomes equal to the constant term Vc of (Formula 1). Thus, when (Formula 1) is used as the voltage characteristic expression of the secondary battery 110, the open circuit voltage Vo at the time of stability can be immediately known from the constant term Vc of (Formula 1).
As another method, the stable open circuit voltage Vo can be obtained by substituting a predetermined required time for stabilization of a predetermined battery voltage into the voltage characteristic equation (Equation 1) into which the optimized adjustment parameter is substituted. it can.

以上説明したように、電圧特性式(式1)が二次電池110の開回路電圧の変化を高精度に推定できることから、上記で推定された充電率の信頼性が高く、これをユーザへのアラーム出力や負荷制御等に利用することが可能となる。   As described above, since the voltage characteristic equation (Equation 1) can accurately estimate the change in the open circuit voltage of the secondary battery 110, the reliability of the charging rate estimated above is high. It can be used for alarm output and load control.

本発明の別の実施形態に係る二次電池電源システムの概略構成を図5に示す。図5に示す本実施形態の二次電池電源システム200では、二次電池110の温度を測定するための温度センサ210を追加しており、温度センサ210で測定した電池温度を制御部130に取り込むように構成している。   FIG. 5 shows a schematic configuration of a secondary battery power supply system according to another embodiment of the present invention. In the secondary battery power supply system 200 of this embodiment shown in FIG. 5, a temperature sensor 210 for measuring the temperature of the secondary battery 110 is added, and the battery temperature measured by the temperature sensor 210 is taken into the control unit 130. It is configured as follows.

調整パラメータBi(i=1〜n)の初期値は、第1の実施形態では固定値として記憶部140に事前に記憶させたものを用いていた。しかしながら、調整パラメータBiを二次電池110の温度によって変更させるようにすることで、開回路電圧をより高精度に推定することが可能となる。従って、調整パラメータBiの初期値についても、これを二次電池110の温度に依存して決定するようにすることで、最小二乗法による調整パラメータの最適化をさらに高速化させることが可能となる。   As the initial value of the adjustment parameter Bi (i = 1 to n), the initial value stored in the storage unit 140 as a fixed value is used in the first embodiment. However, by changing the adjustment parameter Bi according to the temperature of the secondary battery 110, the open circuit voltage can be estimated with higher accuracy. Therefore, by determining the initial value of the adjustment parameter Bi depending on the temperature of the secondary battery 110, the optimization of the adjustment parameter by the least square method can be further accelerated. .

本実施形態では、調整パラメータBiの初期値Bi0を二次電池110の温度から算出する温度関数を事前に決定して記憶部140に記憶させておく。そして、温度センサ210から取得した測定温度をもとに、上記の温度関数から調整パラメータBiの初期値Bi0を算出する。   In the present embodiment, a temperature function for calculating the initial value Bi0 of the adjustment parameter Bi from the temperature of the secondary battery 110 is determined in advance and stored in the storage unit 140. Then, based on the measured temperature acquired from the temperature sensor 210, the initial value Bi0 of the adjustment parameter Bi is calculated from the above temperature function.

本実施形態において調整パラメータBiの初期値Bi0を二次電池110の温度から算出するようにする場合には、図1に示す制御部130での処理の流れが下記のように変更される。まず、ステップS101では、初期値Bi0に代えて上記の温度関数を記憶部140から入力する。   In the present embodiment, when the initial value Bi0 of the adjustment parameter Bi is calculated from the temperature of the secondary battery 110, the flow of processing in the control unit 130 shown in FIG. 1 is changed as follows. First, in step S101, the above temperature function is input from the storage unit 140 in place of the initial value Bi0.

また、ステップS102では、電圧測定値の取得に加えて、温度センサ210から二次電池110の温度測定値を取得する。さらに、ステップS106では、Ai0、Vc0の算出に加えて、上記の温度関数にステップS102で取得した温度測定値を代入することで、初期値Bi0を算出する。上記の温度関数に代入する温度測定値として、温度センサ210から取得した全ての温度測定値を平均して用いてもよいし、いずれか1つを選択して用いてもよい。   In step S102, in addition to obtaining the voltage measurement value, the temperature measurement value of the secondary battery 110 is obtained from the temperature sensor 210. Further, in step S106, in addition to calculating Ai0 and Vc0, the initial value Bi0 is calculated by substituting the temperature measurement value acquired in step S102 into the above temperature function. As temperature measurement values to be substituted into the above temperature function, all temperature measurement values acquired from the temperature sensor 210 may be averaged and used, or any one of them may be selected and used.

上記のように、本実施形態では調整パラメータBiの初期値を二次電池110の温度測定値から算出して用いるようにしていることから、初期値の調整パラメータを用いた電圧特性式(式1)は、電圧測定値により近い開回路電圧を算出するようになる。よって、本実施形態の初期値を用いることで、最小二乗法をさらに高速に収束させて調整パラメータの最適値を決定することが可能となる。   As described above, in the present embodiment, the initial value of the adjustment parameter Bi is calculated from the temperature measurement value of the secondary battery 110 and used, so that the voltage characteristic equation (formula 1) using the adjustment parameter of the initial value is used. ) Calculates an open circuit voltage closer to the voltage measurement. Therefore, by using the initial value of the present embodiment, it is possible to converge the least square method at a higher speed and determine the optimum value of the adjustment parameter.

上記実施形態では、電圧特性式(式1)が指数減衰関数を4項含む例を用いて説明したが、4項に限らずそれ以上の指数減衰関数を含むようにしてもよい。また、二次電池の充放電終了後に行う電圧測定値等の取得個数やサンプリング間隔等についても、二次電池の特性や劣化度等に合わせて適宜変更することができる。さらに、事前に設定された定数であるαやCiについても、二次電池の種類や温度、劣化度等によって選択あるいは算出するようにしてもよい。   In the above-described embodiment, the voltage characteristic formula (Formula 1) has been described using an example including four terms of the exponential decay function. However, the present invention is not limited to four terms and may include exponential decay functions of more terms. In addition, the number of voltage measurement values obtained after the completion of charging / discharging of the secondary battery, the sampling interval, and the like can be changed as appropriate according to the characteristics and the degree of deterioration of the secondary battery. Furthermore, α and Ci, which are constants set in advance, may be selected or calculated according to the type, temperature, deterioration degree, and the like of the secondary battery.

なお、本実施形態では、車両に搭載される車両用二次電池の充電率を推定する構成を備えた二次電池電源システムについて説明したが、本発明は車両用途の二次電池に限られることなく、一般的な二次電池を搭載した各種装置に対して広く適用することができる。   In addition, although this embodiment demonstrated the secondary battery power supply system provided with the structure which estimates the charging rate of the secondary battery for vehicles mounted in a vehicle, this invention is limited to the secondary battery for vehicles use. The present invention can be widely applied to various devices equipped with a general secondary battery.

本発明の二次電池の充電率推定方法を説明するためのフローチャートである。3 is a flowchart for explaining a charging rate estimation method for a secondary battery according to the present invention. 本発明の1実施形態に係る二次電池電源システムの概略構成を示すブロック図である。It is a block diagram which shows schematic structure of the secondary battery power supply system which concerns on 1 embodiment of this invention. 二次電池の充電後60分間にわたる開回路電圧の時間変動を測定した例を示す図である。It is a figure which shows the example which measured the time fluctuation of the open circuit voltage over 60 minutes after charge of a secondary battery. 調整パラメータの初期値を設定する方法を説明するための二次電池の開回路電圧の変化の例を示す図である。It is a figure which shows the example of the change of the open circuit voltage of a secondary battery for demonstrating the method to set the initial value of an adjustment parameter. 本発明の別の実施形態に係る二次電池電源システムの概略構成を示すブロック図である。It is a block diagram which shows schematic structure of the secondary battery power supply system which concerns on another embodiment of this invention.

符号の説明Explanation of symbols

10・・・負荷
100、200・・・二次電池電源システム
110・・・二次電池
120・・・電圧センサ
130・・・制御部
140・・・記憶部
150・・・充電回路
210・・・温度センサ
DESCRIPTION OF SYMBOLS 10 ... Load 100, 200 ... Secondary battery power supply system 110 ... Secondary battery 120 ... Voltage sensor 130 ... Control part 140 ... Memory | storage part 150 ... Charging circuit 210 ...・ Temperature sensor

Claims (10)

充放電終了後の二次電池の開回路電圧(OCV;Open Circuit Voltage)の時間変化を、調整パラメータを有する電圧特性式で近似し、前記電圧特性式から前記二次電池安定時の開回路電圧を算出して前記二次電池の充電率を推定する充電率推定方法であって、
前記二次電池の充放電終了後の所定のデータ取得期間に、前記二次電池の電圧測定値を複数個取得し、
前記データ取得期間における前記電圧測定値の変化量に依存する所定の関係式及び事前に設定された所定の算出方法で前記調整パラメータの初期値を設定し、
前記調整パラメータを前記初期値から漸次修正しながら前記電圧特性式で算出した前記開回路電圧と前記電圧測定値との誤差が最小となるときの前記調整パラメータの最適値を推定する
ことを特徴とする充電率推定方法。
The time variation of the open circuit voltage (OCV) of the secondary battery after the end of charging / discharging is approximated by a voltage characteristic equation having an adjustment parameter, and the open circuit voltage when the secondary battery is stable from the voltage characteristic equation. Is a charging rate estimation method for estimating the charging rate of the secondary battery,
In a predetermined data acquisition period after the end of charging / discharging of the secondary battery, a plurality of voltage measurement values of the secondary battery are acquired,
Setting an initial value of the adjustment parameter with a predetermined relational expression depending on a change amount of the voltage measurement value in the data acquisition period and a predetermined calculation method set in advance;
Estimating an optimum value of the adjustment parameter when an error between the open circuit voltage calculated by the voltage characteristic equation and the voltage measurement value is minimized while gradually adjusting the adjustment parameter from the initial value. How to estimate the charging rate.
前記調整パラメータをAi(i=1〜n)、Bi(i=1〜n)及びVcとし、nを4以上の整数としたとき、前記電圧特定式は時間をtとして、
[数1]
V(t)=A1・exp(B1・t))+A2・exp(B2・t)+・・・
+An・exp(Bn・t)+Vc
で表され、
前記データ取得期間初期の前記電圧測定値をVa、前記データ取得期間終了直前の前記電圧測定値をVbとしたとき、事前に設定された定数α及びCi(i=1〜n)を用いて
[数2]
(C1+C2+・・・+Cn)・β=(Va−Vb)・α
を満たすようパラメータβを決定し、
[数3]
Ai0=Ci・β (i=1〜n)
で算出されるAi0を前記調整パラメータAiの前記初期値とし、
[数4]
Vc0=Va−(Va−Vb)・α
で算出されるVc0を前記調整パラメータVcの前記初期値とする
ことを特徴とする請求項1に記載の充電率推定方法。
When the adjustment parameter is Ai (i = 1 to n), Bi (i = 1 to n), and Vc, and n is an integer of 4 or more, the voltage specifying expression is time t.
[Equation 1]
V (t) = A1 · exp (B1 · t)) + A2 · exp (B2 · t) +.
+ An · exp (Bn · t) + Vc
Represented by
When the voltage measurement value at the initial stage of the data acquisition period is Va and the voltage measurement value immediately before the end of the data acquisition period is Vb, constants α and Ci (i = 1 to n) set in advance are used.
[Equation 2]
(C1 + C2 +... + Cn) · β = (Va−Vb) · α
Parameter β to satisfy
[Equation 3]
Ai0 = Ci · β (i = 1 to n)
Ai0 calculated in step S is used as the initial value of the adjustment parameter Ai,
[Equation 4]
Vc0 = Va− (Va−Vb) · α
The charging rate estimation method according to claim 1, wherein Vc0 calculated in step 1 is used as the initial value of the adjustment parameter Vc.
前記調整パラメータBi(i=1〜n)の初期値は、事前に設定された定数である
ことを特徴とする請求項2に記載の充電率推定方法。
The charging rate estimation method according to claim 2, wherein an initial value of the adjustment parameter Bi (i = 1 to n) is a constant set in advance.
前記調整パラメータBi(i=1〜n)の初期値を前記二次電池の温度から算出する温度関数が事前に設定され、
前記データ取得期間に前記温度を測定して温度測定値を取得し、
前記温度測定値を前記温度関数に代入して前記調整パラメータBi(i=1〜n)の初期値を決定する
ことを特徴とする請求項2に記載の充電率推定方法。
A temperature function for calculating an initial value of the adjustment parameter Bi (i = 1 to n) from the temperature of the secondary battery is set in advance,
Measuring the temperature during the data acquisition period to obtain a temperature measurement value;
The charging rate estimation method according to claim 2, wherein an initial value of the adjustment parameter Bi (i = 1 to n) is determined by substituting the temperature measurement value into the temperature function.
前記Va又は前記Vbの少なくともいずれか一方は、2以上の前記電圧測定値を平均して算出した平均値である
ことを特徴とする請求項2に記載の充電率推定方法。
The charging rate estimation method according to claim 2, wherein at least one of the Va and the Vb is an average value calculated by averaging two or more voltage measurement values.
前記二次電池安定時の前記開回路電圧は、前記調整パラメータの前記最適値を決定したときのVcで与えられる
ことを特徴とする請求項2に記載の充電率推定方法。
The charging rate estimation method according to claim 2, wherein the open circuit voltage when the secondary battery is stable is given by Vc when the optimum value of the adjustment parameter is determined.
前記二次電池安定時の前記開回路電圧は、前記調整パラメータの前記最適値を決定した関数式に予め定めた所定時間を代入して算出する
ことを特徴とする請求項2に記載の充電率推定方法。
3. The charging rate according to claim 2, wherein the open circuit voltage when the secondary battery is stable is calculated by substituting a predetermined time for a function formula that determines the optimum value of the adjustment parameter. Estimation method.
二次電池の充電率を推定する充電率推定装置であって、
前記二次電池の電圧を測定する電圧センサと、
調整パラメータを有する電圧特性式を用いて前記充電率を推定するための演算を実行制御する制御部と、
前記電圧センサで測定された電圧測定値と、前記制御部による前記演算に用いられるデータとを記憶する記憶部と、を備え、
前記制御部は、前記二次電池の充放電終了後の所定のデータ取得期間に、前記電圧センサから前記二次電池の電圧測定値を複数個取得し、前記データ取得期間における前記電圧測定値の変化量に依存する所定の関係式及び事前に設定された所定の算出方法で前記調整パラメータの初期値を設定し、前記調整パラメータを前記初期値から漸次修正しながら前記電圧特性式で算出した前記開回路電圧と前記電圧測定値との誤差が最小となるときの前記調整パラメータの最適値を推定し、前記調整パラメータの最適値を用いた前記電圧特性式から前記二次電池安定時の前記開回路電圧を算出して前記充電率を推定する
ことを特徴とする充電率推定装置。
A charging rate estimation device for estimating a charging rate of a secondary battery,
A voltage sensor for measuring a voltage of the secondary battery;
A control unit that executes and controls calculation for estimating the charging rate using a voltage characteristic equation having an adjustment parameter;
A storage unit for storing a voltage measurement value measured by the voltage sensor and data used for the calculation by the control unit;
The control unit acquires a plurality of voltage measurement values of the secondary battery from the voltage sensor during a predetermined data acquisition period after the end of charging and discharging of the secondary battery, and sets the voltage measurement value in the data acquisition period. The initial value of the adjustment parameter is set by a predetermined relational expression that depends on the amount of change and a predetermined calculation method set in advance, and the adjustment parameter is calculated from the initial value and the voltage characteristic formula is used to calculate the initial value. The optimum value of the adjustment parameter when the error between the open circuit voltage and the measured voltage value is minimized is estimated, and the A charging rate estimation apparatus characterized by calculating a circuit voltage and estimating the charging rate.
前記制御部は、前記調整パラメータをAi、Bi(i=1〜n)及びVcとし、nを4以上の整数としたとき、次式
[数5]
V(t)=A1・exp(B1・t))+A2・exp(B2・t)+・・・
+An・exp(Bn・t)+Vc
を前記電圧特定式として用い、前記データ取得期間初期の前記電圧測定値をVa、前記データ取得期間終了直前の前記電圧測定値をVbとしたとき、事前に設定された定数α及びCi(i=1〜n)を用いて
[数6]
(C1+C2+・・・+Cn)・β=(Va−Vb)・α
を満たすようパラメータβを決定し、
[数7]
Ai0=Ci・β (i=1〜n)
で算出されるAi0を前記調整パラメータAiの前記初期値とし、
[数8]
Vc0=Va−(Va−Vb)・α
で算出されるVc0を前記調整パラメータVcの前記初期値として用いている
ことを特徴とする請求項7に記載の充電率推定装置。
When the control parameter is Ai, Bi (i = 1 to n) and Vc, and n is an integer of 4 or more,
[Equation 5]
V (t) = A1 · exp (B1 · t)) + A2 · exp (B2 · t) +.
+ An · exp (Bn · t) + Vc
Is used as the voltage specifying formula, Va is the voltage measurement value at the initial stage of the data acquisition period, and Vb is the voltage measurement value immediately before the end of the data acquisition period. 1-n)
[Equation 6]
(C1 + C2 +... + Cn) · β = (Va−Vb) · α
Parameter β to satisfy
[Equation 7]
Ai0 = Ci · β (i = 1 to n)
Ai0 calculated in step S is used as the initial value of the adjustment parameter Ai,
[Equation 8]
Vc0 = Va− (Va−Vb) · α
The charging rate estimation apparatus according to claim 7, wherein Vc0 calculated in step (b) is used as the initial value of the adjustment parameter Vc.
請求項7又は請求項8に記載の充電率推定装置を備える
ことを特徴する二次電池電源システム。
A secondary battery power supply system comprising the charge rate estimation device according to claim 7 or 8.
JP2006279478A 2006-10-13 2006-10-13 CHARGE RATE ESTIMATION METHOD, CHARGE RATE ESTIMATION DEVICE, AND SECONDARY BATTERY POWER SUPPLY SYSTEM Active JP4785056B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006279478A JP4785056B2 (en) 2006-10-13 2006-10-13 CHARGE RATE ESTIMATION METHOD, CHARGE RATE ESTIMATION DEVICE, AND SECONDARY BATTERY POWER SUPPLY SYSTEM

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006279478A JP4785056B2 (en) 2006-10-13 2006-10-13 CHARGE RATE ESTIMATION METHOD, CHARGE RATE ESTIMATION DEVICE, AND SECONDARY BATTERY POWER SUPPLY SYSTEM

Publications (2)

Publication Number Publication Date
JP2008096328A true JP2008096328A (en) 2008-04-24
JP4785056B2 JP4785056B2 (en) 2011-10-05

Family

ID=39379314

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006279478A Active JP4785056B2 (en) 2006-10-13 2006-10-13 CHARGE RATE ESTIMATION METHOD, CHARGE RATE ESTIMATION DEVICE, AND SECONDARY BATTERY POWER SUPPLY SYSTEM

Country Status (1)

Country Link
JP (1) JP4785056B2 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011122164A1 (en) 2010-03-30 2011-10-06 古河電気工業株式会社 State-of-charge estimation method, state-of-charge estimation device, and secondary-battery power system
WO2012066643A1 (en) * 2010-11-16 2012-05-24 古河電気工業株式会社 Method and apparatus for detecting state of electrical storage device
JP2013034339A (en) * 2011-08-03 2013-02-14 Toyota Motor Corp Vehicle control device
WO2013183480A1 (en) 2012-06-05 2013-12-12 株式会社豊田自動織機 State-of-charge estimation method and state-of-charge estimation device
JP5541391B1 (en) * 2013-03-28 2014-07-09 株式会社豊田自動織機 Open circuit voltage estimation method and vehicle
JP2014241656A (en) * 2013-06-11 2014-12-25 株式会社豊田自動織機 Charging start control method of secondary battery, and charging device
WO2015145616A1 (en) * 2014-03-26 2015-10-01 株式会社日立製作所 Lithium ion secondary battery control device and control method, and lithium ion secondary battery module

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004530880A (en) * 2001-05-02 2004-10-07 ハネウェル・インターナショナル・インコーポレーテッド Method of determining state of charge of battery by measuring open circuit voltage
JP2005043339A (en) * 2003-07-09 2005-02-17 Furukawa Electric Co Ltd:The Method for estimating charging rate, device for estimating charging rate, battery system and battery system for vehicle
JP2005100969A (en) * 2003-08-22 2005-04-14 Furukawa Electric Co Ltd:The Method and device for measuring internal impedance of secondary battery, and power supply system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004530880A (en) * 2001-05-02 2004-10-07 ハネウェル・インターナショナル・インコーポレーテッド Method of determining state of charge of battery by measuring open circuit voltage
JP2005043339A (en) * 2003-07-09 2005-02-17 Furukawa Electric Co Ltd:The Method for estimating charging rate, device for estimating charging rate, battery system and battery system for vehicle
JP2005100969A (en) * 2003-08-22 2005-04-14 Furukawa Electric Co Ltd:The Method and device for measuring internal impedance of secondary battery, and power supply system

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8903667B2 (en) 2010-03-30 2014-12-02 Furukawa Electric Co., Ltd. State-of-charge estimation method, state-of-charge estimation device, and secondary-battery power system
JP2011209237A (en) * 2010-03-30 2011-10-20 Furukawa Electric Co Ltd:The State-of-charge estimation method and device, and secondary-battery power system
CN102803980A (en) * 2010-03-30 2012-11-28 古河电气工业株式会社 State-of-charge estimation method, state-of-charge estimation device, and secondary-battery power system
EP2555007A4 (en) * 2010-03-30 2017-07-26 Furukawa Electric Co., Ltd. State-of-charge estimation method, state-of-charge estimation device, and secondary-battery power system
WO2011122164A1 (en) 2010-03-30 2011-10-06 古河電気工業株式会社 State-of-charge estimation method, state-of-charge estimation device, and secondary-battery power system
WO2012066643A1 (en) * 2010-11-16 2012-05-24 古河電気工業株式会社 Method and apparatus for detecting state of electrical storage device
CN102893169A (en) * 2010-11-16 2013-01-23 古河电气工业株式会社 Method and apparatus for detecting state of electrical storage device
JP5823292B2 (en) * 2010-11-16 2015-11-25 古河電気工業株式会社 Method and apparatus for detecting state of power storage device
JP2013034339A (en) * 2011-08-03 2013-02-14 Toyota Motor Corp Vehicle control device
CN104335058A (en) * 2012-06-05 2015-02-04 株式会社丰田自动织机 State-of-charge estimation method and state-of-charge estimation device
US20150153420A1 (en) * 2012-06-05 2015-06-04 Kabushiki Kaisha Toyota Jidoshokki a corporation State-of-charge estimation method and state-of-charge estimation device
EP2860539A4 (en) * 2012-06-05 2016-01-27 Toyota Jidoshokki Kk State-of-charge estimation method and state-of-charge estimation device
JPWO2013183480A1 (en) * 2012-06-05 2016-01-28 株式会社豊田自動織機 Charging rate estimation method and charging rate estimation device
WO2013183480A1 (en) 2012-06-05 2013-12-12 株式会社豊田自動織機 State-of-charge estimation method and state-of-charge estimation device
JP2014190929A (en) * 2013-03-28 2014-10-06 Toyota Industries Corp Open circuit voltage estimation method and vehicle
WO2014155829A1 (en) * 2013-03-28 2014-10-02 株式会社豊田自動織機 Open-circuit voltage estimation method and vehicle
JP5541391B1 (en) * 2013-03-28 2014-07-09 株式会社豊田自動織機 Open circuit voltage estimation method and vehicle
JP2014241656A (en) * 2013-06-11 2014-12-25 株式会社豊田自動織機 Charging start control method of secondary battery, and charging device
WO2015145616A1 (en) * 2014-03-26 2015-10-01 株式会社日立製作所 Lithium ion secondary battery control device and control method, and lithium ion secondary battery module

Also Published As

Publication number Publication date
JP4785056B2 (en) 2011-10-05

Similar Documents

Publication Publication Date Title
US8903667B2 (en) State-of-charge estimation method, state-of-charge estimation device, and secondary-battery power system
EP3745151B1 (en) Method and apparatus for correcting soc, battery management system and storage medium
JP4785056B2 (en) CHARGE RATE ESTIMATION METHOD, CHARGE RATE ESTIMATION DEVICE, AND SECONDARY BATTERY POWER SUPPLY SYSTEM
JP4015128B2 (en) CHARGE RATE ESTIMATION METHOD, CHARGE RATE ESTIMATION DEVICE, BATTERY SYSTEM, AND VEHICLE BATTERY SYSTEM
JP4865523B2 (en) Battery charge rate estimation method, battery charge rate estimation device, and battery power supply system
JP4532416B2 (en) Battery discharge capability determination method, battery discharge capability determination device, and power supply system
US7733063B2 (en) Apparatus for calculating quantity indicating charged state of on-vehicle battery
US9594122B2 (en) Estimating state of charge (SOC) and uncertainty from relaxing voltage measurements in a battery
US20190137573A1 (en) Power storage apparatus and controlling method for the same
US20130325379A1 (en) Internal resistance estimation device and method of estimating internal resistance
US20090248334A1 (en) Method for estimating the charge of a motor vehicle battery
US20090134843A1 (en) Method and apparatus for detecting internal electric state of in-vehicle secondary battery
JP2007292778A (en) Method for estimating state of battery charge
WO2013183480A1 (en) State-of-charge estimation method and state-of-charge estimation device
JP2023116546A (en) Charging rate estimation device
JP2007179968A (en) Battery status control device
JP6895786B2 (en) Secondary battery status detection device and secondary battery status detection method
JP5129029B2 (en) Open voltage value estimation method and open voltage value estimation apparatus
JP4865627B2 (en) Battery remaining capacity estimation method, battery remaining capacity estimation apparatus, and battery power supply system
TWI758760B (en) Battery control device and battery capacity estimation method
JP5367320B2 (en) Secondary battery state detection method, state detection device, and secondary battery power supply system
JP2012500980A (en) Method and controller for determining the state of charge of a battery
CN109633451B (en) Energy storage system self-track parameter calibration method and SOC estimation method
WO2012046285A1 (en) Battery status estimation method and power supply system
EP3923396B1 (en) Chargeable battery state detection device and chargeable battery state detection method

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20100921

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101004

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20101119

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: 20110701

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110707

R151 Written notification of patent or utility model registration

Ref document number: 4785056

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140722

Year of fee payment: 3

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350