(57)【要約】      （修正有）【課題】動作条件下で電流を直接検出することなくバッテリパックまたはセルの電流及び充電状態を推定する【解決手段】バッテリパックまたはセルの特性評価電圧及び特性評価電流の測定された時間的傾向に関連する、バッテリパックまたはセルの特性評価データを取得するステップと、複数の特性評価温度値のそれぞれの値に関連する、特性評価データを取得するステップと、バッテリパックまたはセルの、温度及び充電状態に応じて、バッテリパックまたはセルの動作モデルの複数のパラメータを決定するために特性評価データを処理するステップと、バッテリパックまたはセルが動作条件にある間に、バッテリパックまたはセルの動作電圧及び動作温度を測定するステップと、測定された動作電圧、測定された動作温度及び複数のパラメータに基づく、動作モデルによって、バッテリパックまたはセルの動作電流を推定するステップと、を備える【選択図】図１
2019032299.tif 000002
【特許請求の範囲】
【請求項１】
バッテリパックまたはセルによって分配される動作電流（Ｉ）を推定する方法であって、
−バッテリパックまたはセルの特性評価電圧（Ｖｍ）及び特性評価電流（Ｉｍ）の測定された時間的傾向に関連する、バッテリパックまたはセルの特性評価データを取得するステップであって、それぞれの時間的傾向は、複数の特性評価温度値（Ｔｋ）のそれぞれの値に関連する、特性評価データを取得するステップと、
−バッテリパックまたはセルの、温度及び充電状態「ｓｔａｔｅ  ｏｆ  ｃｈａｒｇｅ」（ＳＯＣ）に応じて、バッテリパックまたはセルの動作モデルの複数のパラメータ（Ｐ）を決定するために前記特性評価データを処理するステップと、
−バッテリパックまたはセルが動作条件にある間に、バッテリパックまたはセルの動作電圧（Ｖ）及び動作温度（Ｔ）を測定するステップと、
−測定された動作電圧（Ｖ）、測定された動作温度（Ｔ）及び前記複数のパラメータ（Ｐ）に基づく、前記動作モデルによって、バッテリパックまたはセルの動作電流（Ｉ）を推定するステップと、を備え、
ここで、処理するステップは、前記複数の特性評価温度値（Ｔｋ）のそれぞれの値のために実行される、次のステージ、
−特性評価温度（Ｔｋ）に関連する特性評価電圧（Ｖｍ）及び特性評価電流（Ｉｍ）のそれぞれの時間的傾向に沿って複数の時間観測窓（Ｗｉ）を特定し、それぞれの特性評価電流（Ｉｍｉ）及び特性評価電圧（Ｖｍｉ）値を検知するステップと、
−それぞれの観測窓（Ｗｉ）において、バッテリパックまたはセルの充電状態（ＳＯＣｉ）のそれぞれの値を計算するステップと、
−それぞれの観測窓（Ｗｉ）において、パラメータ（Ｐ）のそれぞれの名目上の値に応じた、バッテリパックまたはセルの動作モデルによって、それぞれの推定した電圧（ＶＡＢｉ）を計算し、時間観測窓（Ｗｉ）の推定された電圧（ＶＡＢｉ）と特性評価電圧（Ｖｍｉ）の間の差によって決まるそれぞれの誤差関数（Ｅｉ）を決定するステップと、
−それぞれの観測窓（Ｗｉ）において、前記誤差関数（Ｅｉ）を最小化することによって、バッテリパックまたはセルのモデルの複数のパラメータ（Ｐ）のそれぞれについての実際の値（Ｐｉ）を計算するステップと、
−バッテリパックまたはセルの温度と充電状態に応じた前記複数のパラメータ（Ｐ）を得るために、それぞれの充電状態（ＳＯＣｉ）及び特性評価温度（Ｔｋ）と、観測窓（Ｗｉ）において計算される、パラメータの実際の値（Ｐｉ）を関連付けるステップと、を備える方法。
【請求項２】
バッテリパックまたはセルの動作モデルは、電気回路モデルであり、モデルパラメータは、電気回路パラメータを備える、請求項１に記載の方法。
【請求項３】
処理ステップはさらに、
−それぞれの観測窓（Ｗｉ）において、バッテリパックまたはセルのそれぞれの負荷なし電圧（Ｖｏｃｉ）値を決定ステップと、
−観測窓（Ｗｉ）で得られた複数の負荷なし電圧（Ｖｏｃｉ）及び充電状態（ＳＯＣｉ）値に基づいて、バッテリパックまたはセルの負荷なし電圧（Ｖｏｃ）及び充電状態（ＳＯＣ）の関係を定義するステップと、
ここで、それぞれの観測窓（Ｗｉ）において、パラメータの実際の値（Ｐｉ）を計算するステップは、それぞれの負荷なし電圧（Ｖｏｃｉ）値と、バッテリパックまたはセルの負荷なし電圧及び充電状態（ＳＯＣ）の前記関係も考慮するステップと、を備える請求項２に記載の方法。
【請求項４】
複数の観測窓（Ｗｉ）を特定するステップは、特性評価電流（Ｉｍ）が過渡電流現象の完全な消耗時に、ゼロであり、その結果として、負荷のない電圧（Ｖｏｃｉ）は、特性評価電圧（Ｖｍｉ）に対応する、時間窓を特定するステップを備える、請求項３に記載の方法。
【請求項５】
バッテリパックまたはセルの充電状態（ＳＯＣｉ）のそれぞれの値を計算するステップは、
−特性評価電流の時間的傾向に基づいて、バッテリパックまたはセルの充電状態（ＳＯＣｉ）のそれぞれの値を計算する、請求項１乃至４のいずれか１項に記載の方法。
【請求項６】
バッテリパックまたはセルの電気回路モデルは、負荷のない電圧（Ｖｏｃ）発生器、バッテリまたはセルの内部抵抗（Ｒ０）及び１以上の回路グループの直列を備え、それぞれは、それぞれの回路グループ抵抗及びそれぞれの回路グループ容量の並列を備え、
回路パラメータＰの群は、負荷のない電圧（Ｖｏｃ）、バッテリまたはセルの内部抵抗（Ｒ０）、１以上の回路グループ抵抗（Ｒ１、Ｒ２）、及び１以上の回路グループ容量（Ｃ１、Ｃ２）を備える、請求項２乃至５のいずれか１項に記載の方法。
【請求項７】
バッテリパックまたはセルの電気回路モデルは、２つの回路群を備え、回路パラメータＰの群は、第１の回路群抵抗（Ｒ１）、第１の回路群容量（Ｃ１）、第２の回路群抵抗（Ｒ２）、第２の回路群容量（Ｃ２）を備える、請求項６に記載の方法。
【請求項８】
バッテリパックまたはセルの動作モデルは、訓練された予測アルゴリズムを備え、パラメータ（Ｐ）は、訓練された予測アルゴリズムのパラメータであり、
複数のパラメータ（Ｐ）を決定するステップは、取得されるバッテリパックまたはセルの特性評価データによる予測アルゴリズムを訓練するステップを備える、請求項１に記載の方法。
【請求項９】
バッテリパックまたはセルの特性評価データを取得するステップは、
−動作条件下でバッテリパックまたはセルを使う前に、それぞれの特性評価温度（Ｔｋ）値において、バッテリパックまたはセルの特性評価電圧（Ｖｍ）の及び特性評価電流（Ｉｍ）の時間的傾向を測定するステップ、または、
−バッテリパックまたはセルの動作条件下の使用の前に処理される特性評価手順から特性評価データを取得するステップ、を備える請求項１乃至８のいずれか１項に記載の方法。
【請求項１０】
処理するステップはさらに、
−動作条件下のバッテリパックまたはセルを容易に利用できる電気保存手段で、入力として、温度と充電状態値を受け取り、パラメータ（Ｐ）のそれぞれの値を提供するように構成された、ルックアップテーブルとして、バッテリパックまたはセルの動作モデルのパラメータ（Ｐ）を保存するステップ、を備える請求項１乃至９のいずれか１項に記載の方法。
【請求項１１】
バッテリパックまたはセルの動作電流を推定するステップは、前記保存ルックアップテーブルにアクセスするように構成され、入力として、バッテリパックまたはセルに組み込まれたそれぞれのセンサから、測定された動作電圧（Ｖ）及び動作温度（Ｔ）を受け取るように構成された、バッテリパックまたはセルと関連する、電子プロセッサによって処理される、請求項１０に記載の方法。
【請求項１２】
さらに
−バッテリパックまたはセルの熱的モデルを定義するステップと、
−少なくとも１つの測定時間間隔において、前記間隔で測定される、バッテリパックまたはバッテリセルの電圧（Ｖ）及び温度（Ｔｃ）の値に基づいて、及びバッテリパックまたはバッテリセルの前記温度（Ｔｃ）に対応するパラメータ（Ｐ）に基づいて、バッテリパックまたはセルの動作モデルによって、浪費した電力（Ｐｅ）を推定するステップと、
−バッテリパックまたはバッテリセルの前記温度値（Ｔｃ）に基づいて、環境温度値（Ｔｅ）に基づいて、ここで両方の温度値は、前記間隔で測定され、及びバッテリパックまたはバッテリセルの前記温度（Ｔｃ）に対応するパラメータ（Ｐ）に基づいて、バッテリパックまたはセルの熱的モデルによって、浪費した熱出力（Ｐｄ）を推定するステップと、
−推定された、浪費した電力（Ｐｅ）と浪費した熱出力（Ｐｄ）の差（ΔＰ）を決定するステップと、
−劣化及び／または経年劣化現象を考慮するパラメータ（Ｐ）の訂正した値を得るように、前記推定した差（ΔＰ）によって、パラメータを訂正するステップと、を備える請求項１乃至１１のいずれか１項に記載の方法。
【請求項１３】
バッテリパックまたはセルの熱的モデルは、１以上の熱的セルモデルからなり、
それぞれの熱的セルモデルは、熱的セル容量（ＣＴ）、熱的セル抵抗（ＲＴ）、セル温度（Ｔｃ）及び環境温度（Ｔｅ）によって特性評価される、請求項１２に記載の方法。
【請求項１４】
訂正するステップは、
−電力と熱出力の差（ΔＰ）により、訂正因子（Ｇ）を決定するステップと、
−バッテリまたはセルの内部抵抗（Ｒ０）の更新された推定を得るために、前記訂正因子（Ｇ）を、内部抵抗（Ｒ０）の値に乗じるステップと、
−バッテリまたはセルの内部抵抗（Ｒ０）の更新された推定による、電力（Ｐｅ）を再計算するステップと、
−前記差（ΔＰ）があらかじめ決められたしきい値以下に減るまで、決定すること、乗じること、及び再計算することの前記ステップを反復するステップと、
−前記反復の終わりに得られた更新された推定をバッテリまたはセルの内部抵抗（Ｒ０）値と見なすステップと、を備える請求項１２または１３のいずれか１項に記載の方法。
【請求項１５】
さらに、
−動作条件下でバッテリパックまたはセルを使用する前に、異なる温度値及び使用時間に応じて、バッテリまたはセルの内部抵抗（Ｒ０）の値の変化の特性を評価するステップと、
−内部抵抗の特性評価データを保存するステップと、
−方法が適用されたとき、バッテリパックまたはセルの動作条件下で、バッテリパックまたはセルの実際の使用の経過時間を考慮して、バッテリまたはセルの内部抵抗（Ｒ０）の前記特性評価データによって、バッテリまたはセルの内部抵抗（Ｒ０）の訂正値を使用するステップと、を備える請求項６に記載の方法。
【請求項１６】
バッテリパックまたはセルの充電状態（ＳＯＣ）を推定する方法であって、
−請求項１乃至１５のいずれか１項に記載の、バッテリパックまたはセルの動作時間に沿って、バッテリパックまたはセルにより分配される動作電流（Ｉ）を推定するステップと、
−バッテリパックまたはセルの動作開始時から現在時まで動作する、動作電流（Ｉ）の推定された時間的傾向に基づき、初期充電状態（ＳＯＣ０）からバッテリパックまたはセルの充電状態（ＳＯＣ）を推定するステップと、を備える方法。
【請求項１７】
一体化した電流センサを除外しても、バッテリパックまたはセルを制御するように構成されたバッテリパックまたはセルの制御装置であって、
装置は、バッテリパックまたはセルの測定された動作電圧（Ｖ）値とバッテリパックまたはセルの測定された動作温度（Ｔ）値を受け取ように構成され、
装置は、請求項１乃至１５のいずれか１項に記載された、バッテリパックまたはセルの動作電流（Ｉ）を推定する方法を実行する及び／または請求項１６に記載されたバッテリパックまたはセルの充電状態（ＳＯＣ）を推定する方法を実行するように構成された電子処理手段を備える、バッテリパックまたはセル制御装置。
【請求項１８】
一体化された電流センサを備えることなく、動作するように構成されたバッテリパックまたはセルであって、
バッテリパックまたはセルの動作電圧（Ｖ）を検知するように構成された、電圧センサと、
バッテリパックまたはセルの動作温度（Ｔ）を検知するように構成された、温度センサと、
請求項１７に記載の制御装置と、を備えるバッテリパックまたはセル。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電気器具を制御する、電子方法及びシステムの技術的部門に関連する。
【０００２】
特に、本発明は、動作条件下で電流を直接検知することなく、バッテリ、または自動車のバッテリパックまたは、同様に、バッテリセル（単にセル「ｃｅｌｌ」としてこの明細書で定義される）の、動作電流及び充電状態を推定する方法に関連する。
【０００３】
さらに、本発明は、前記方法を実行する、バッテリパックまたはセル、電流センサを備える必要なく、前記方法の長所によって、動作電流及び充電状態をモニタできる自動車のバッテリパックまたはセルの制御装置に関連する。
【背景技術】
【０００４】
電気化学的なセル（特に、自動車分野において、電気またはハイブリッド駆動自動車の動作システム及びバッテリという観点において）から構成されるバッテリ、バッテリパックまたはセルの電圧及び電流を知ることは、バッテリ、バッテリパックまたはセルの正しい管理及び制御のとても重要な情報である。
【０００５】
バッテリパックまたはセルの複数のセルの電圧、電流、温度の取得は、セルのパフォーマンス及び耐久性の低下、また、極端な場合、重大な結果として、バッテリパックまたはセルの火災または爆発まで至る、熱暴走「ｔｈｅｒｍａｌ  ｒｕｎａｗａｙ」として知られる現象をもたらす、２次物理化学プロセスのきっかけを避けるように、セルそれ自身の充電状態を知り、バッテリパックまたはセルの正しい管理及び保護するために、必要とされる。
【０００６】
さらに、いくつかのモジュールまたは全てのバッテリパックまたはセルに共通の、電流を知ることは、セルの充電状態（ＳＯＣ）を決定するために役立つ。
【０００７】
共通の背景技術の解決方法において、バッテリまたはセルの電流を知ることは、通常、バッテリパックまたはセルそれ自身に設置された、シャントタイプの物理電流センサまたはホール効果ベースのセンサを使って、そのような電流の直接検知によって得られる。
【０００８】
シャントタイプセンサについて言えば、電流は、バッテリまたはセル電流が循環する既知の値の抵抗の電圧を測定することによって、決定される。
【０００９】
ホール効果ベースのセンサについて言えば、電流は、バッテリまたはセルの電流によって横切られ、電磁場を受ける導体に発生するホール電圧を測定することによって、決定される。
【００１０】
前述の両方の場合は、バッテリパックそれ自身またはセルの複雑さとコストの増加の不利益を必然的に伴う、複数の電流センサを備える、バッテリパックまたはセルを提供する必要性を想定している。
【００１１】
それゆえ、電流センサによって、そのような電流の直接の取得に関わらず、バッテリまたはセルの電流（すなわち、「バッテリまたはセル電流」）を推定する方法を特定する必要が生じる。
【００１２】
このアプローチは、しかしながら、実際に使用できるように十分に正確な電流推定を得るさらなる技術的な問題に遭遇する。実際に、必要とされた正確さは、（例えば５％に相当またはそれより小さい）制限されたパーセンテージの範囲内で推定の誤差を最小化することを提供する。
【００１３】
そのような正確なバッテリまたはセルの電流の推定を得ることは、実際のバッテリまたはセルの電流が、得られた瞬間において、その瞬間に検知されたバッテリまたはセルの動作電圧に依存するだけでなく、複雑な振る舞いで、バッテリまたはセルの温度と充電状態（今後、その過程及び寿命）にも依存するので、全く容易でない。それゆえ、バッテリまたはセルの電圧及び温度を知ることによって、必要とされる正確さで、バッテリまたはセル電流を得ることができる、単一の分析式はない。
【００１４】
さらに、バッテリまたはセル電流と他のバッテリまたはセルの電気的、熱的量の関係は、劣化及び経年劣化にも依存するので、さらに、そのような態様を適切に考慮する推定方法を有する必要性が、さらに感じられる。
【００１５】
既知の解決法は、この関連で、上記で言及された必要を満足しない。
【００１６】
前述を参照して、主に自動車応用の分野で、特に電気またはハイブリッドの自動車のバッテリパックとの関連で、バッテリまたはセルの電流及び充電状態を推定し、バッテリまたはセルの電流及び充電状態に関して十分正確な情報を提供する一方で、動作条件下で、バッテリまたはセルの電流を直接取得する必要を避ける方法の必要性が感じられる。
【発明の概要】
【００１７】
少なくとも部分的に、背景技術を参照して上記に示される欠点を克服すること、特に考えられる技術部門で感じられる上述の必要を満たすことができるような、電流それ自身の直接の取得を必要とすることなく、バッテリまたはセルの電流及び充電状態を十分に正確に推定することを提供することができる、バッテリパックまたはセルの動作電流及び充電状態を推定する方法を提供することが、本発明の目的である。
【００１８】
特に、本発明は、シャントまたはホール効果タイプの物理電流センサを用いる必要なく、バッテリまたはセルが電流によって横切られるとき、バッテリまたはセルの端部において電圧のみを測定することによって、バッテリまたはセルの電流を決定する方法を提供することが目的である。それによって、バッテリまたはセルの端部において、電圧を直接測定することのみで、バッテリパックまたはセルを完全に管理することが可能であり、それにより、温度及び経年劣化を含む、異なる変化に応じて、電圧測定とその内部抵抗の特性評価から電流を決定する。
【００１９】
そのような目的は、請求項１による電流を推定する方法で達成される。
【００２０】
そのような方法のさらなる実施形態は、請求項２−１５によって定義される。
【００２１】
請求項１６によるバッテリパックまたはセルの充電状態を推定する方法が、本発明の別の目的である。
【００２２】
本発明の目的は、また前述の方法を実行することができる、請求項１７で定義される、バッテリパックまたはセルの制御装置及び、請求項１８で定義される、電流センサなしで、前述の制御器具によって制御されるバッテリパックまたはセルでもある。
【図面の簡単な説明】
【００２３】
本発明による電流を推定する方法及び充電状態を推定する方法、及び関連する制御装置及びバッテリパックまたはセルのさらなる特徴及び利点は、添付した図を参照して、示す及び制限しない実施例を手段として得られる、好ましい実施形態の続く記載から明白になるであろう。
【００２４】
【図１】図式的に、本発明による推定方法の実施形態を、ブロック図によって、示す。
【図２】本発明による方法の実施形態に適合された、バッテリまたはセルの電気モデルの図を示す。
【図３】図式的に、方法の実施形態で用いられた、プロセスアルゴリズムを、ブロック図によって、示す。
【図４Ａ】方法の入力データを構成することができる、バッテリパックまたはセルの、電圧の測定された時間的傾向を示す。
【図４Ｂ】方法の入力データを構成することができる、バッテリパックまたはセルの、特性評価電流の測定された時間的傾向を示す。
【図５Ａ】図４Ａのグラフの一部を拡大した詳細を示す。
【図５Ｂ】図４Ｂのグラフの一部を拡大した詳細を示す。
【図６】バッテリまたはセルの充電状態ＳＯＣの推定された時間的傾向を示す。
【図７】バッテリまたはセルの負荷なし電圧Ｖｏｃの推定された時間的傾向を示す。
【図８】バッテリまたはセルの充電状態及び負荷なし電圧（関数）の間の推定された関係の実施例を示す。
【図９】バッテリまたはセルの充電状態及び負荷なし電圧の間の逆の（逆関数）推定された関係の実施例を示す。
【図１０】バッテリまたはセルの端部において測定された電圧と本発明の方法によって推定された対応する電圧の比較のグラフを示す。
【図１１】方法の典型的な実施形態で推定された、充電状態に応じたバッテリまたはセル（方法で使われるバッテリまたはセルの可能なモデルに関連する）の回路パラメータの値を示す。
【図１２】本発明による方法の実施形態に適合された、バッテリまたはセルの熱的モデルの図を示す。
【図１３】図式的に、方法の実施形態で用いられた、劣化及び／または経年劣化に応じた、パラメータを訂正及び更新した、アルゴリズムを、ブロック図によって、示す。
【図１４】発明の方法で推定された、バッテリまたはセルの測定された電流から導き出す充電状態とバッテリまたはセルの電流から導き出す充電状態の比較のグラフを示す。
【図１５】本発明の方法で推定された測定されたバッテリまたはセルの電流及びバッテリまたはセルの電流の比較のグラフを示す。
【発明を実施するための形態】
【００２５】
図１−１５を参照して、バッテリまたはセルによって分配された動作電流Ｉを推定する方法が、これから記載される。
【００２６】
方法は、第１に、バッテリパックまたはセルの特性評価電圧Ｖｍ及び特性評価電流Ｉｍの測定された時間的傾向に関連する、バッテリパックまたはセルの特性評価データを取得するステップを備える。前記時間的傾向のそれぞれは、複数の特性評価温度値Ｔｋのそれぞれの値に関連付けられる。
【００２７】
さらにバッテリパックまたはセルの温度及び充電状態ＳＯＣに応じた、バッテリパックまたはセルの動作モデルの複数のパラメータＰを決定するために前記の特性評価データを処理する方法が提供される。
【００２８】
その後、バッテリパックまたはセルが動作条件下にある間、測定された動作電圧Ｖ、測定された動作温度Ｔ及び前述の複数のパラメータＰ（前にバッテリパックまたはセルの温度及び充電状態に応じて決定されたように）に基づいて、前記の動作モデルによって、バッテリパックまたはセルの動作電圧Ｖ及び動作温度Ｔを測定し、バッテリパックまたはセルの動作電流Ｉを推定する方法が提供される。
【００２９】
動作電流Ｉは、バッテリパックまたはセルの瞬間的な動作電圧Ｖ、瞬間的な温度Ｔが検知されたとき、バッテリパックまたはセルに流れる瞬間的な電流であることに注目すべきである。そのような動作電流Ｉは、それゆえ、前記のパラメータＴ、Ｖ及びＰに基づいて推定される。
【００３０】
パラメータＰは、次に、すでに上記で記述されたように、下記で記載されるように、バッテリパックまたはセルの温度及び充電状態ＳＯＣに基づいて、さらに処理することによって、得られる。
【００３１】
本明細書で記載された方法において、処理するステップは、前記の複数の特性評価温度値Ｔｋのそれぞれの値で処理される、いくつかのステージを備える。
【００３２】
特に、処理するステップは、特性評価温度Ｔｋに関連した、特性評価電圧Ｖｍ及び特性評価電流Ｉｍのそれぞれの時間的傾向に沿った複数の時間観測窓Ｗｉを特定することと、それぞれ特性評価電流Ｉｍｉと特性評価電圧Ｖｍｉ値を検知することを含む。
【００３３】
それゆえ、それぞれの観測窓Ｗｉにおいて、バッテリパックまたはセルの充電状態ＳＯＣｉのそれぞれの値を計算する方法が提供される。
【００３４】
さらに、それぞれの観測窓Ｗｉにおいて、パラメータＰのそれぞれの名目上の値に応じて、バッテリパックまたはセルの動作モデルによって、それぞれ推定された電圧ＶＡＢｉを計算し、時間観測窓Ｗｉの推定された電圧ＶＡＢｉと特性評価電圧Ｖｍｉの差に依存するそれぞれの誤差関数Ｅｉを決定する方法が提供される。
【００３５】
さらに、それぞれの観測窓Ｗｉにおいて、前記の誤差関数Ｅｉを最小化することによって、バッテリパックまたはセルのモデルの複数のパラメータＰのそれぞれの実際の値Ｐｉを計算する方法が提供される。
【００３６】
最終的に、それぞれの充電状態ＳＯＣｉ及び特性評価温度Ｔｋと観測窓Ｗｉにおいて計算される、パラメータの実際の値Ｐｉを関連付けるステップは、バッテリパックまたはセルの温度及び充電状態ＳＯＣに応じた、前記の複数のパラメータＰを得るために、提供される。
【００３７】
方法の実施形態によると、さらに、それぞれの観測窓Ｗｉにおいて、バッテリパックまたはセルの負荷なし電圧Ｖｏｃｉのそれぞれを決定し、観測窓Ｗｉにおいて得られる複数の負荷なし電圧Ｖｏｃｉ及び充電状態ＳＯＣｉ値に基づいて、バッテリパックまたはセルの負荷なし電圧Ｖｏｃ及び充電状態ＳＯＣの関係を定義する、特性評価データを処理する前記のステップが提供される。
【００３８】
そのような場合において、それぞれの観測窓Ｗｉにおいて、パラメータの実際の値Ｐｉを計算するステップは、それぞれの負荷なし電圧Ｖｏｃｉ値及び決定されたバッテリパックまたはセルの負荷なし電圧Ｖｏｃ及び充電状態ＳＯＣの関係も考慮する。
【００３９】
特定の実施形態によると、複数の観測窓Ｗｉを特定するステップは、特性評価電流Ｉｍが過渡電流現象の完全な消耗時に、ゼロであり、その結果、負荷なし電圧Ｖｏｃｉが特性評価電圧Ｖｍｉに対応する、時間窓を特定することを備える。
【００４０】
特定の実施形態によると、バッテリパックまたはセルの充電状態の値を計算するステップは、特性評価電流の時間的傾向に基づく、バッテリパックまたはセルの充電状態（ＳＯＣｉ）のそれぞれの値を計算することを備える。
【００４１】
方法の異なる可能な実施形態によると、バッテリパックまたはセルの異なる動作モデルを使うことができる。
【００４２】
方法の典型的な実施形態によると、バッテリパックまたはセルの動作モデルは、電気回路モデルである。そのような場合において、モデルパラメータは、電気回路パラメータを含む。
【００４３】
図２で示された実施例において、バッテリパックまたはセルの電気回路モデルは、負荷のない電圧Ｖｏｃ発生器、１以上の回路グループの直列に配置されたバッテリまたはセルの内部抵抗Ｒ０を備え、それぞれは、回路グループ抵抗及び回路グループ容量の並列を備える。そのような場合において、回路パラメータＰの群は、負荷のない電圧Ｖｏｃ、バッテリまたはセルの内部抵抗Ｒ０、１以上の回路グループ抵抗Ｒ１、Ｒ２、及び１以上の回路グループ容量Ｃ１、Ｃ２を備える。
【００４４】
特に図２の実施例は、２次モデル（通常、正確で、同時に、相対的に単一のモデルを有するように使われる）、すなわち２つの回路グループを備えるモデルを示す。そのような場合において、回路パラメータＰの群は、バッテリまたはセルの内部抵抗Ｒ０及び負荷なし電圧Ｖｏｃに加えて、第１の回路群抵抗Ｒ１、第１の回路群容量Ｃ１、第２の回路群抵抗Ｒ２、第２の回路群容量Ｃ２を備える。
【００４５】
いくつかのパラメータ、偶数の非回路パラメータで特性評価される、他のタイプのモデル、偶数の非回路モデルを使用する方法が提供されることに注目すべきである。
【００４６】
１つの実施形態において、例えば、バッテリパックまたはセルの動作モデルは、訓練された予測アルゴリズムを備え、パラメータＰは、訓練された予測アルゴリズムのパラメータである。そのような場合において、複数のパラメータＰを決定するステップは、取得されるバッテリパックまたはセルの特性評価データによる予測アルゴリズムを訓練することを備える。
【００４７】
例えば、Ｎパラメータのニューラルネットワークを使用することができる。
【００４８】
電気回路モデル（例えば、図２に示されているもの）の実施例に戻ると、モデルパラメータ（Ｖｏｃ、Ｒ０、Ｒ１、Ｒ２、Ｃ１、Ｃ２）は、充電状態ＳＯＣと温度の関数であることに注目すべきである。
【００４９】
電圧、電流及び温度の関係は、次の式（１）で表される。
【００５０】
【数１】
2019032299.tif 000003
【００５１】
ここで、電圧ＶＡＢ（ｔ）は、バッテリセルまたはセルの終端の電圧を示し、Ｉ（ｔ）はセルを循環する電流である。
【００５２】
温度依存は、それぞれ温度に依存する、パラメータ（Ｖｏｃ、Ｒ０、Ｒ１、Ｒ２、Ｃ１、Ｃ２）に含まれるので、潜在している。
【００５３】
ジュール効果によって消散した電力とエネルギは、次の式（２）で表される。
【００５４】
【数２】
2019032299.tif 000004
【００５５】
すでに上記で記述されたように、電気パラメータは温度に依存し、さらに、バッテリセルまたはセルの充電状態、特に負荷なし電圧Ｖｏｃに依存する。
【００５６】
モデルを効果的に使うために、それゆえ、異なる充電状態ＳＯＣと温度Ｔ値に対するパラメータを特性評価する必要がある。
【００５７】
この目的を達成するために、特性評価データを処理するステップを、一般論として、すでに上記で記載されて、提供している。
【００５８】
さらなる詳細は、これによって、図２の電気回路モデルが使われている方法の実施形態を参照して、処理するステップに用いられているアルゴリズムに提供されている。
【００５９】
アルゴリズムを処理するそのような実施例は、図３に図式的に示されている。
【００６０】
温度を変えてあらかじめ定義され、設定されたＫに対応する、Ｋの時間的傾向、すなわち電圧及び電流傾向を有すると見なす。ｋ＝１に関して、第１の温度は設定され（Ｔｋ＝Ｔ１）、温度Ｔ１のバッテリまたはセルで取得される、電圧及び電流の時間的傾向が、選択される。
【００６１】
電圧と電流の時間的傾向は、例えば、バッテリまたはセルの、「オフライン」で特性を評価するステップで取得される。それゆえ、そのような電圧とそのような電流は、それぞれ上記で定義された、特性評価電圧Ｖｍ及び特性評価電流Ｉｍに対応する。
【００６２】
特性評価電圧Ｖｍと特性評価電流Ｉｍの時間的傾向の実施例は、温度２８℃、２０アンペアアウア（Ａｈ）チタン酸リチウム（ＬＴＯ）バッテリセルまたはセルについて言えば、それぞれ、図４Ａ及び４Ｂに示される。
【００６３】
Ｖｍ及びＩｍの前記時間的傾向において、電流がゼロであり、例えば、セルそれ自身の挙動が決定されると考えることができることなど、時間間隔の間、ゼロ値で残る、すなわち、回路モデルの容量Ｃ１及びＣ２がかなり放電された状態で、時間観測窓Ｗｉは、特定される。図４Ａ及び４Ｂにおいて、観測窓Ｗｉの終わりは、「Ｘ」マークによって示される。
【００６４】
第１の時間観測窓の拡大が、特性評価電圧Ｖｍに関して、図５Ａで、特性評価電流Ｉｍに関して、図５Ｂで示される。
【００６５】
次のステージは、測定された電流Ｉｍに基づいて、充電状態（ＳＯＣ）を計算することである。充電状態ＳＯＣを、例えば、次の式（４）に基づいて、いわゆる「クーロンカウンタ」技術で、計算することができる。
【００６６】
【数３】
2019032299.tif 000005
【００６７】
ここで、Ｑｎｏｍは、１クーロンの放電割合で推定された、セルの名目上の容量（示されている実施例において、２０Ａｈの値を有する）である。
【００６８】
時間の関数として計算されている、充電状態ＳＯＣは、図６に示される。
【００６９】
その後、負荷なし電圧Ｖｏｃも、時間の関数として計算される。この目的を達成するために、観測窓Ｗｉの最初と終わりの点は、電流がゼロであり、それゆえ、容量Ｃ１及びＣ２の端部の電圧が無視できると考慮できる、「変化しない」点であることに注目すべきである。それゆえ、そのような点において、負荷なし電圧Ｖｏｃは、測定された電圧Ｖｍと一致する。他の点において、負荷なし電圧Ｖｏｃの計算は、前に推定された充電状態ＳＯＣの関数の線形補間として実行されることができる。
【００７０】
時間の関数として計算された負荷なし電圧Ｖｏｃは、図７に示される。
【００７１】
前記情報に基づいて、負荷なし電圧Ｖｏｃの関数として充電状態ＳＯＣ、すなわち、図が図８に示されるＳＯＣ＝ｆ（Ｖｏｃ）、を計算することが可能である。
【００７２】
同様に、逆関数を負荷なし電圧Ｖｏｃを充電状態ＳＯＣと結び付けて、すなわち図が図９に示されるＶｏｃ＝ｆ−１（ＳＯＣ）で得ることができる。
【００７３】
ＶｏｃとＳＯＣの関係を参照して、第ｉの観測窓Ｗｉが設定された時点で、次の関係が、有効である。
【００７４】
【数４】
2019032299.tif 000006
【００７５】
ここで、下付文字１及び２は、それぞれ時間観測窓の最初及び終わりの点を示す。より小さいΔＳＯＣ＝ＳＯＣ２−ＳＯＣ１はそのような近似関係をより正確にする。
【００７６】
現時点で、残っている回路パラメータ（Ｒ０、Ｒ１、Ｒ２、Ｃ１、Ｃ２）を得ることができる：それぞれの観測窓Ｗｉにおいて、モデルは、それぞれ測定された電流Ｉｍｉで促され、前記５つの回路パラメータの組み合わせが、測定された電圧Ｖｍとモデルによって推定された電圧ＶＡＢの誤差を最小化するように、決定される。
【００７７】
そのような誤差の最小化プロセスの終わりに、推定された電圧ＶＡＢは、図１０（図のスケールにおいて、２つの電圧は、実際には区別ができない）に示されるように、測定された特性評価電圧Ｖｍを高い精度で、再現する。
【００７８】
特性評価温度Ｔｋ（特定の場合において、２８℃）についての処理手順の終わりに、モデルのすべてのパラメータは、図１１の図に示されるように、充電状態ＳＯＣの関数として得られる。
【００７９】
処理手順は複数の特性評価温度Ｔｋで反復される。理想的には、手順は、いくつもの温度値で実行されることができる。
【００８０】
全手順の終わりに、それゆえ、バッテリパックまたはセルのモデルの回路パラメータのそれぞれの値は、充電状態ＳＯＣの値の関数として、及び異なる動作温度値の関数として入手可能である。
【００８１】
処理手順の結果は、「ルックアップテーブル」の形で、集められ、保存することができる。
【００８２】
実際には、方法の実施形態によると、処理はさらに電子保存手段（動作条件でバッテリパックまたはセルにアクセス可能な）に、入力として温度と充電状態を受け取り、出力としてパラメータＰのそれぞれの値を提供するように構成された、ルックアップテーブルとして、バッテリパックまたはセルの動作モデルのパラメータＰを保存するステップを備える。
【００８３】
そのような場合において、方法の実施形態によると、バッテリパックまたはセルの動作電流を推定するステップは、前記保存されたルックアップテーブルにアクセスするように構成され、入力として、バッテリパックまたはセルに組み込まれたそれぞれのセンサから生じる測定された動作電圧Ｖと動作温度Ｔ値を受け取るように構成された、バッテリパックまたはセルに関連する電子プロセッサによって、実行される。
【００８４】
本明細書において、バッテリパックまたはセルの特性評価データを取得するステップは、広い意味で、限定されない実施例を手段として、ここから下記で記載される実施形態に記載されているそれらを含む、多くの異なる方法で実行されることができる、前の特性を評価する段階に続いて、そのようなデータを利用可能にするように、理解されるに違いない。
【００８５】
方法の可能な実施形態において、バッテリパックまたはセルの特性評価データを取得するステップは、動作条件下でバッテリパックまたはセルを使用する前に、それぞれの特性評価温度Ｔｋ値における、バッテリパックまたはセルの特性評価電圧Ｖｍ及び特性評価電流Ｉｍの時間的傾向を測定することを備える。
【００８６】
方法の別の可能な実施形態において、バッテリパックまたはセルの特性評価データを取得するステップは、バッテリパックまたはセルの動作条件下で使用する前に、実行される特性を評価する手順で特性評価データを取得することを備える。
【００８７】
方法のさらなる可能な実施形態において、特性評価データを取得するステップは、バッテリパックまたはセルの入手可能なデータシートから特性評価データを導き出すことを備える。
【００８８】
バッテリパックまたはセルの動作モデルのパラメータが、前に記載されたように、バッテリまたはセルの温度と充電状態の両方を考慮して、決定されると、バッテリパックまたはセルの動作電流Ｉを推定することが可能である。
【００８９】
これは、例えば、次の方程式（５）に基づいて、なされることができる。
【００９０】
【数５】
2019032299.tif 000007
【００９１】
ここで、
【００９２】
【数６】
2019032299.tif 000008
である。
【００９３】
反復法の、次のステージ（ステップ：Ｔｓ）を離散化することで、次の式が得られる。
【００９４】
【数７】
2019032299.tif 000009
【００９５】
【数８】
2019032299.tif 000010
【００９６】
ここで、下付文字ＯＬＤをつけて、反復法の前のステージＴｓから得られたＩ、ＩＲ１、ＩＲ２、Ｃ１、Ｃ２、Ｒ１及びＲ２の値が示される。
【００９７】
方法の前に記載された実施形態は、適当な方法で、いかなる動作条件でも、そのような電流を直接検知する必要なく、バッテリパックまたはセルの電流の信頼性のある推定を得るように、バッテリパックまたはセルの充電状態及び動作温度を考慮することができる。
【００９８】
電流の実際の値からの推定の狂いのさらなる原因は、バッテリまたはセル及びその構成要素の劣化及び／または経年劣化に由来することがかもしれない。
【００９９】
実際には、回路モデルのパラメータが、バッテリまたはセルが新しいか特性を評価することが実行されたときのものである場合有効である。充電状態と温度のそのようなパラメータの依存は、バッテリまたはセルの劣化により、徐々に変化する。それゆえ、厳しく言うと、特性を評価することは、バッテリパックまたはセルがハイブリッドまたは電気自動車に取り付けられた後、これは簡単ではなく、または、実現可能ではないが、周期的に実行されるべきである。
【０１００】
そのような態様に対処し、及び経年劣化及び劣化に関連する要因も考慮して、モデルのパラメータを正しく推定しようとするために、方法のさらなる実施形態は、続くさらなるステップを提供する。
【０１０１】
第１に、バッテリパックまたはセルの熱的モデルが、前述の動作モデルに追加して、定義される。
【０１０２】
その後、少なくとも１つの測定時間間隔（ｔ０−ｔ１）において、消散した電力Ｐｅを、バッテリパックまたはバッテリセルの電圧Ｖと温度Ｔｃの前記間隔（ｔ０−ｔ１）で測定された値に基づいて、及びバッテリパックまたはバッテリセルの前記温度Ｔｃにおいて推定されるパラメータＰに基づいて、バッテリパックまたはセルの動作モデルによって推定する。
【０１０３】
さらに、同じ測定時間間隔（ｔ０−ｔ１）において、消散した熱出力Ｐｄを、バッテリパックまたはバッテリセルの温度Ｔｃの前記値に基づいて、及び前記間隔（ｔ０−ｔ１）で測定された環境温度Ｔｅ値に基づいて、及びバッテリパックまたはバッテリセルの温度Ｔｃにおいて推定されるパラメータＰに基づいて、バッテリパックまたはセルの熱的モデルによって推定する。
【０１０４】
その後、前に推定された消散した電力Ｐｅと消散した熱出力Ｐｄの差ΔＰを決定するためのステップを実行する。
【０１０５】
最終的に、そのような劣化及び／または経年劣化の現象を考慮するために、パラメータＰの更新値を得るように、前記推定された差ΔＰに基づいてパラメータＰを、修正する。
【０１０６】
方法の実施形態によると、バッテリパックまたはセルの熱的モデルは、それぞれの熱的セルモデルが、熱的セル容量ＣＴ、熱的セル抵抗ＲＴ、セル温度Ｔｃ及び環境温度Ｔｅによって特性評価される、１以上の熱的セルモデルから成り立って使われる。
【０１０７】
熱的セル容量ＣＴ及び熱的セル抵抗ＲＴパラメータは、セル及びバッテリパックまたはセルの特定の形態に特有である。
【０１０８】
特に、方法に使われる熱的モデルの実施例は、例えば、図式的に図１２に示されるものなど、１次熱的モデルである。
【０１０９】
そのような熱的モデルによると、次の関係（６）が有効である。
【０１１０】
【数９】
2019032299.tif 000011
【０１１１】
ここで、パラメータの意味するものは、すでに上記に記載されている。
【０１１２】
時間間隔（ｔ０−ｔ１）で消散した熱エネルギＥＴは次の式（７）によって計算されることができる。
【０１１３】
【数１０】
2019032299.tif 000012
【０１１４】
現時点で、与えられた時間間隔で消散したエネルギ及び／または熱及び電力を利用可能に有するので、パラメータの値は更新されることができる。
【０１１５】
特に、図１３のブロック図で示される実施形態によって、モデルパラメータを修正する前記ステップは、前に記載されたように、劣化を除外して決定される値から始まる、内部抵抗Ｒ０パラメータの修正を備える。
【０１１６】
この目的を達成するために、修正因子Ｇが、推定された電力と推定された熱出力の差ΔＰに基づいて、決定される。その後、劣化を除外して有効な、バッテリまたはセルの内部抵抗Ｒ０値は、バッテリまたはセルの内部抵抗Ｒ０の更新された推定を得るために、前記修正因子Ｇで乗じられる。
【０１１７】
その後、内部抵抗Ｒ０の更新された推定に基づいて、電力Ｐｅを再計算し、その後、差ΔＰがあらかじめ決められたしきい値未満に減らされるまで、決定することと、乗じることと、再計算する前記ステップを繰り返すステップが提供される。最終的に、前記反復法の終わりにおいて得られる更新された推定は、バッテリまたはセルの内部抵抗Ｒ０値として考えられる。
【０１１８】
また、劣化及び／または経年劣化の存在するバッテリまたはセルの電気動作パラメータを推定においても、目的とされる方法の代替の実施形態は、動作条件下でバッテリパックまたはセルを使う前に、異なる温度値と使用の時間に応じたバッテリまたはセルの内部抵抗Ｒ０の値の変化の特性評価をし、バッテリまたはセルの内部抵抗の特性評価データを保存するステップを備える。
【０１１９】
さらに、方法が、バッテリパックまたはセルの動作条件下で適用されたとき、バッテリまたはセルの内部抵抗Ｒ０の修正された値は、バッテリパックまたはセルの実際の使用の経過時間を考慮して、バッテリまたはセルの内部抵抗Ｒ０の前記特性評価データによって、用いられる。
【０１２０】
本発明は、また、バッテリパックまたはセルの充電状態ＳＯＣを推定する方法を備える。そのような方法は、上記で記載された実施形態のいずれか１つによって、バッテリパックまたはセルの動作時間に沿って、バッテリパックまたはセルによって分配される動作電流Ｉを推定するステップを備える。その後、バッテリパックまたはセルの動作開始の瞬間から現在の瞬間まで、動作する、動作電流Ｉの推定された時間的傾向に基づく、最初の充電状態ＳＯＣ０からバッテリパックまたはセルの充電状態ＳＯＣを推定する方法が提供される。
【０１２１】
特に、実施形態によると、前述の式はこの目的において使われることができる（４）。
【０１２２】
【数１１】
2019032299.tif 000013
【０１２３】
そのような方法は、比較が、物理センサによって、取得される電流を測定することによって推定される充電状態ＳＯＣＩｍｉｓと、物理電流センサを除外して、この方法に基づいて推定される充電状態ＳＯＣＩｅｓｔの間でなされる、図１４の実施例を手段として示されるように、３％と同じかそれより低い誤差で、充電状態を正確に推定することができる。
【０１２４】
図１５において、比較は、測定された電流Ｉｍｉｓと本発明の方法によって推定されたように推定された電流Ｉｅｓｔでなされる。注目されることができるように、５％と同じかそれより低い、高い正確性が得られる。
【０１２５】
本発明は、さらに一体化された電流センサを除外しても、バッテリパックまたはセルを制御するように構成されたバッテリパックまたはセル制御装置を備え、装置は、バッテリパックまたはセルの測定された動作電圧Ｖ値及びバッテリパックまたはセルの測定された動作温度Ｔ値を受け取るように構成されている。
【０１２６】
装置は、上記に記載された実施形態のいずれかによって、バッテリパックまたはセルの動作電流Ｉを推定する方法を実行、及び／または上記に記載された実施形態のいずれか１つによって、バッテリパックまたはセルの充電状態ＳＯＣを推定する方法を実行するように構成された電子処理手段を備える。
【０１２７】
本発明は、さらに、一体化された電流センサを備えることなく、動作するように構成されたバッテリパックまたはセルを備える。
【０１２８】
バッテリパックまたはセルは、バッテリパックまたはセルの動作電圧Ｖを検知するように構成された、電圧センサ、バッテリパックまたはセルの動作温度Ｔを検知するように構成された、温度センサ、前に記載されたものによる制御装置を備える。
【０１２９】
注目されることができるように、本発明の目的は、記載された推定する方法によって、そのそれぞれの動作特徴の長所によって、完全に達成される。
【０１３０】
実際には、バッテリまたはセルの電流及びバッテリまたはセルの充電状態を推定する方法は、適切なセンサによって、バッテリまたはセルの電流を直接取得する必要なく、バッテリまたはセルの電圧及び温度のみを「オンライン」での取得（すなわち、バッテリまたはセルは動作条件下である間）に基づく。
【０１３１】
これは、さらに、複雑さとコストを減らすことに関して、明らかな利点を導き出すバッテリパックまたはセルから電流センサを排除する（並びに対応する配線、電流信号を取得するハードウェアシステム及び関連のあるソフトウェア管理も排除する）ことが可能である。
【０１３２】
本方法が物理電流センサの存在下で用いられたとしても、電流の追加の間接の推定を得る方法を提供すること、この方法の使用は、ついでの利益として、安全の観点から有益な、合同性または妥当性の推定をする可能性をどちらにせよ提供するであろうことは注目すべきである（ＩＳＯ２６２６２規格）。
【０１３３】
さらに、前に詳細に記載されているように、バッテリまたはセルの電流及びバッテリまたはセルの充電状態を推定する方法は、実用的に、効果的に用いられるように、十分正確な推定を得ることができる。特に、記載された方法は、５％より低い誤差でバッテリパックまたはセルを流れる瞬間電流を推定し、３％より低い誤差で充電状態を推定することができる。
【０１３４】
同様な利点を、前記推定方法を実行する、本発明による制御装置、並びに装置及び制御システムを使って、電流センサのなくてもよい、そのような推定方法から利益を得る、本発明によるバッテリパックまたはセルを参照して特定することができる。
【０１３５】
当業者は、付随する要求を満足するために、続く請求項の範囲から逸脱することなく、上記に記載された電流を推定する方法、充電状態を推定する方法、制御装置及びバッテリパックまたはセルの実施形態の他の機能的な均等物に、要素を改良し、適合し、置き換えることができる。
【０１３６】
可能な実施形態に付属するように記載された特徴のそれぞれは、記載された他の実施形態から独立に達成されることができる。用語「備える」（ｃｏｍｐｒｉｓｉｎｇ）は、他の要素またはステップを排除することがなく、用語「１つの」（ａ）または「１つの」（ｏｎｅ）は複数を排除することがないことも留意すべきである。より明確に説明するために、さまざまな部分に適切なハイライトをする要求を許容するので、図はノンスケールで描かれている。
【図１】
2019032299.tif 000014
【図２】
2019032299.tif 000015
【図３】
2019032299.tif 000016
【図４Ａ】
2019032299.tif 000017
【図４Ｂ】
2019032299.tif 000018
【図５Ａ】
2019032299.tif 000019
【図５Ｂ】
2019032299.tif 000020
【図６】
2019032299.tif 000021
【図７】
2019032299.tif 000022
【図８】
2019032299.tif 000023
【図９】
2019032299.tif 000024
【図１０】
2019032299.tif 000025
【図１１】
2019032299.tif 000026
【図１２】
2019032299.tif 000027
【図１３】
2019032299.tif 000028
【図１４】
2019032299.tif 000029
【図１５】
2019032299.tif 000030
【外国語明細書】
2019032299.tif 000031
2019032299.tif 000032
2019032299.tif 000033
2019032299.tif 000034
2019032299.tif 000035
2019032299.tif 000036
2019032299.tif 000037
2019032299.tif 000038
2019032299.tif 000039
2019032299.tif 000040
2019032299.tif 000041
2019032299.tif 000042
2019032299.tif 000043
2019032299.tif 000044
2019032299.tif 000045
2019032299.tif 000046
2019032299.tif 000047
2019032299.tif 000048
2019032299.tif 000049
2019032299.tif 000050
2019032299.tif 000051
2019032299.tif 000052
2019032299.tif 000053
2019032299.tif 000054
2019032299.tif 000055
2019032299.tif 000056
2019032299.tif 000057
2019032299.tif 000058
2019032299.tif 000059
2019032299.tif 000060
2019032299.tif 000061
2019032299.tif 000062
2019032299.tif 000063
2019032299.tif 000064
2019032299.tif 000065
2019032299.tif 000066
2019032299.tif 000067
2019032299.tif 000068
2019032299.tif 000069