JP7167415B2 - 温度推定モデル決定方法および装置、温度推定モデルが適用されたバッテリ管理システム - Google Patents

温度推定モデル決定方法および装置、温度推定モデルが適用されたバッテリ管理システム Download PDF

Info

Publication number
JP7167415B2
JP7167415B2 JP2021520343A JP2021520343A JP7167415B2 JP 7167415 B2 JP7167415 B2 JP 7167415B2 JP 2021520343 A JP2021520343 A JP 2021520343A JP 2021520343 A JP2021520343 A JP 2021520343A JP 7167415 B2 JP7167415 B2 JP 7167415B2
Authority
JP
Japan
Prior art keywords
temperature
error
profile
parameter set
temperature estimation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2021520343A
Other languages
English (en)
Other versions
JP2022504862A (ja
Inventor
リー、ボヒュン
パク、ミンヒー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
LG Energy Solution Ltd
Original Assignee
LG Energy Solution 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 LG Energy Solution Ltd filed Critical LG Energy Solution Ltd
Publication of JP2022504862A publication Critical patent/JP2022504862A/ja
Application granted granted Critical
Publication of JP7167415B2 publication Critical patent/JP7167415B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/396Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • G01K7/16Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using resistive elements
    • G01K7/22Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using resistive elements the element being a non-linear resistance, e.g. thermistor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • G01K7/42Circuits effecting compensation of thermal inertia; Circuits for predicting the stationary value of a temperature
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • G01K7/42Circuits effecting compensation of thermal inertia; Circuits for predicting the stationary value of a temperature
    • G01K7/427Temperature calculation based on spatial modeling, e.g. spatial inter- or extrapolation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/3644Constructional arrangements
    • G01R31/3648Constructional arrangements comprising digital calculation means, e.g. for performing an algorithm
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/382Arrangements for monitoring battery or accumulator variables, e.g. SoC
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/385Arrangements for measuring battery or accumulator variables
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M10/00Secondary cells; Manufacture thereof
    • H01M10/42Methods or arrangements for servicing or maintenance of secondary cells or secondary half-cells
    • H01M10/425Structural combination with electronic components, e.g. electronic circuits integrated to the outside of the casing
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M10/00Secondary cells; Manufacture thereof
    • H01M10/42Methods or arrangements for servicing or maintenance of secondary cells or secondary half-cells
    • H01M10/48Accumulators combined with arrangements for measuring, testing or indicating the condition of cells, e.g. the level or density of the electrolyte
    • H01M10/486Accumulators combined with arrangements for measuring, testing or indicating the condition of cells, e.g. the level or density of the electrolyte for measuring temperature
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M10/00Secondary cells; Manufacture thereof
    • H01M10/42Methods or arrangements for servicing or maintenance of secondary cells or secondary half-cells
    • H01M10/425Structural combination with electronic components, e.g. electronic circuits integrated to the outside of the casing
    • H01M2010/4271Battery management systems including electronic circuits, e.g. control of current or voltage to keep battery in healthy state, cell balancing
    • 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

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • General Chemical & Material Sciences (AREA)
  • Nonlinear Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Secondary Cells (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Description

[関連出願との相互引用]
本出願は2019年8月29日付韓国特許出願第10-2019-0106742号に基づく優先権の利益を主張し、当該韓国特許出願の文献に開示されたすべての内容は本明細書の一部として含まれる。
本開示は温度推定モデル決定方法および装置、温度推定モデル決定方法によって決定された温度推定モデルが適用されたバッテリ管理システムに関する。
バッテリモジュールの充電時間を短縮するためには、バッテリモジュールを構成する複数のバッテリセル間の温度偏差をバッテリ管理システム(Battery Management System)で正確に感知し、感知した温度に基づいてバッテリモジュールに対する充電制御が行われなければならない。
しかし、バッテリ生産工程および費用の限界により、バッテリモジュール内には1-2個の温度センサのみが設けられ、バッテリモジュール内で最高温度および最低温度の地点を正確に検出できない。
本開示の実施形態によりバッテリの最高温度および最低温度の指定を推定できる温度推定モデル決定方法および装置、温度推定モデルが適用されたバッテリ管理システムを提供する。
発明の一特徴によるバッテリモジュール内の温度を推定する温度推定モデルを決定する方法は、温度推定のために選択されたモデル種類によるn個のパラメータおよびn個のパラメータそれぞれの区間数に基づいてすべてのパラメータセットを計算する段階、すべてのパラメータセットそれぞれに基づいたモデルに予備検証プロファイルの入力データを代入して温度を推定する段階、すべてのパラメータセットそれぞれに対して、予備検証プロファイルの測定された温度と推定された温度の間の差による第1誤差を算出する段階、すべてのパラメータセットに対する第1誤差それぞれが第1臨界値より小さいかを比較する段階、比較結果第1誤差のうちの第1臨界値より小さい第1誤差に対応するパラメータセットを第1パラメータセットとして抽出する段階、第1パラメータセットに基づいたモデルを用いて予備検証プロファイルの時系列的に変わる入力データにより推定された温度と予備検証プロファイルの測定された温度の間の誤差のうちの最も大きい誤差である最大誤差を算出する段階、最大誤差が第2臨界値より小さいかを比較する段階、比較結果第2臨界値より小さい最大誤差に対応する加重値を乗じた値と第1誤差のうちの第1パラメータセットに対応する第1誤差に対応する加重値を乗じた値を足した結果に基づいて第2パラメータセットを抽出する段階、第2パラメータセットに再検証プロファイルの入力データにより温度を推定する段階、再検証プロファイルの測定された温度と推定された温度の間の差による第2誤差を算出する段階、および第2誤差が第3臨界値より小さいと、第2パラメータセットに基づいたモデルを温度推定モデルに決定する段階を含み得る。
予備検証プロファイルおよび再検証プロファイルは、少なくとも時系列的に変わる温度センサの温度値、バッテリ電流、冷却水温度値、および空気温度値を入力データとして含み得る。
推定された温度は、推定されたバッテリモジュールの最高温度または最低温度を含み、測定された温度は、測定されたバッテリモジュールの最高温度または最低温度を含み得る。
第1誤差および第2誤差のうちの少なくとも一つは、二乗平均平方根誤差(Root Mean Square Error,RMSE)であり得る。
温度推定モデル決定方法は、第1誤差がすべて第1臨界値以上であれば、基準プロファイルを変更する段階をさらに含み、基準プロファイルは、温度推定モデルをビルドアップするための初期値として、温度推定に必要な因子および測定された温度を含み得る。
温度推定モデル決定方法は、第1パラメータセットすべての最大誤差が第2臨界値以上であれば、基準プロファイルを変更する段階をさらに含み得る。
温度推定モデル決定方法は、第2パラメータセットすべての第2誤差が第3臨界値以上であれば、基準プロファイルを変更する段階をさらに含み得る。
第1パラメータセットを抽出する段階は、第1誤差のうちの小さい順に上位i個に属する第1誤差に対応するパラメータセットを第1パラメータセットとして抽出する段階を含み得る。
第2パラメータセットを抽出する段階は、足した結果が小さい順に上位j個に属する足した結果に対応するパラメータセットを第2パラメータセットとして抽出する段階を含み得る。
本発明の他の特徴によるバッテリ管理システムは、温度推定モデル決定方法によって決定された温度推定モデルを含み得る。
本発明のまた他の特徴によるバッテリモジュール内の温度を推定する温度推定モデルを決定する装置は、温度推定のために選択されたモデル種類によるn個のパラメータおよびn個のパラメータそれぞれの区間数に基づいてすべてのパラメータセットを計算するパラメータセット生成部、すべてのパラメータセットそれぞれに基づいたモデルに予備検証プロファイルの入力データを代入して温度を推定するモデリング部、およびすべてのパラメータセットそれぞれに対して、予備検証プロファイルの測定された温度と推定された温度の間の差による第1誤差を算出する誤差計算部を含む。
モデリング部は、すべてのパラメータセットに対する第1誤差それぞれが第1臨界値より小さいかを比較し、比較結果第1誤差のうちの第1臨界値より小さい第1誤差に対応するパラメータセットを第1パラメータセットとして抽出する。
温度推定モデル決定装置は、第1パラメータセットに基づいたモデルを用いて予備検証プロファイルの時系列的に変わる入力データにより推定された温度と予備検証プロファイルの測定された温度の間の誤差のうちの最も大きい誤差である最大誤差を算出し、最大誤差が第2臨界値より小さいかを比較し、比較結果第2臨界値より小さい最大誤差に対応する加重値を乗じた値と第1誤差のうちの第1パラメータセットに対応する第1誤差に対応する加重値を乗じた値を足した結果に基づいて第2パラメータセットを抽出する最大誤差計算部、および第2パラメータセットに再検証プロファイルの入力データにより温度を推定し、再検証プロファイルの測定された温度と推定された温度の間の差による第2誤差を算出し、第2誤差が第3臨界値より小さいと、第2パラメータセットに基づいたモデルを温度推定モデルに決定する最終モデル決定部をさらに含み得る。
予備検証プロファイルおよび再検証プロファイルは、少なくとも時系列的に変わる温度センサの温度値、バッテリ電流、冷却水温度値、および空気温度値を入力データとして含み得る。
推定された温度は、推定されたバッテリモジュールの最高温度または最低温度を含み、測定された温度は、測定されたバッテリモジュールの最高温度または最低温度を含み得る。
第1誤差および第2誤差のうちの少なくとも一つは、二乗平均平方根誤差(Root Mean Square Error,RMSE)であり得る。
本開示の実施形態によりバッテリの最高温度および最低温度の指定を推定できる温度推定モデル決定方法および装置、温度推定モデルが適用されたバッテリ管理システムを提供する。
一実施形態による温度推定モデル決定装置を示す図である。 一実施形態による温度推定モデル決定方法を示すフローチャートである。 一実施形態による最高温度推定モデル決定方法を示すフローチャートである。 一実施形態による最低温度推定モデル決定方法を示すフローチャートである。 一実施形態による温度推定モデルが適用されたバッテリ管理システムを含むバッテリシステムを示す図である。
本開示の実施形態はシステム識別(system identification)方法に基づいてサーミスタ(thermistor)測定値、バッテリモジュールに流れる電流、冷却水および外気温度値を用いて実際のバッテリモジュール内の最高温度および最低温度を推定するロジックを導き出して最終的にサーミスタの限界を改善することができる。バッテリモジュール内の最高温度および最低温度推定によりバッテリ管理システムはより精巧な充電制御が可能である。
サーミスタ(thermistor)の温度値、バッテリ電流、冷却水温度値、および外気温度値に基づいたバッテリモジュール内の最高/最低温度の間の相関関係を推定する温度推定モデルを決定する方法およびそれを実現した装置に関する。
以下、添付する図面を参照して本明細書に開示された実施形態を詳細に説明するが、同一または似の構成要素は同一、類似の図面符号を付与し、これに対する重複する説明は省略する。以下の説明で使われる構成要素に対する接尾辞の「モジュール」および/または「部」は明細書作成の容易さだけを考慮して付与されたり混用されたりするものであり、それ自体が互いに区別される意味または役割を有するものではない。また、本明細書に開示された実施形態を説明することにおいて関連する公知技術に対する具体的な説明が本明細書に開示された実施形態の要旨を曖昧にすると判断される場合はその詳細な説明を省略する。また、添付する図面は本明細書に開示された実施形態を容易に理解できるようにするためのものであり、本明細書に開示された技術的思想は添付する図面によって制限されるものではなく、本発明の思想および技術範囲に含まれるすべての変更、均等物または代替物を含むものとして理解しなければならない。
第1、第2などのように序数を含む用語は多様な構成要素を説明するために使われるが、構成要素は用語によって限定されない。用語は一つの構成要素を他の構成要素から区別するための目的にのみ使用される。
ある構成要素が他の構成要素に「連結されて」いるか、または「接続されて」いると言及された時には、その他の構成要素に直接連結されているかまたは接続されていてもよいが、中間に他の構成要素が存在し得ると理解されなければならない。反面、ある構成要素が他の構成要素に「直接連結されて」いるか、「直接接続されて」いると言及された時には、中間に他の構成要素が存在しないと理解されなければならない。
本出願で、「含む」または「有する」などの用語は明細書上に記載された特徴、数字、段階、動作、構成要素、部品またはこれらを組み合わせたものが存在することを指定するためであり、一つまたはそれ以上の他の特徴や数字、段階、動作、構成要素、部品またはこれらを組み合わせたものなどの存在または付加の可能性をあらかじめ排除しないものとして理解されなければならない。
図1は一実施形態による温度推定モデル決定装置を示す図である。
図1に示すように、温度推定モデル決定装置1は、パラメータセット生成部10、モデリング部20、誤差計算部30、最大誤差計算部40、および最終モデル決定部50を含む。
図1に示す各構成は最高温度および最低温度推定モデル決定装置の構成を説明するための一例として、少なくとも二つの構成が合わさるか、一つの構成が少なくとも二つの構成に分離されるか、または追加的な構成が付加されるなどの変更が可能であり、図1に示す実施形態が発明を制限しない。
図2は一実施形態による温度推定モデル決定方法を示すフローチャートである。
図3は一実施形態による最高温度推定モデル決定方法を示すフローチャートである。
図4は一実施形態による最低温度推定モデル決定方法を示すフローチャートである。
以下、図1から図4を参照して一実施形態による温度推定モデル決定装置および方法について説明する。
先に、パラメータセット生成部10は実験データの入力を受け、基準プロファイル、予備検証プロファイル、および再検証プロファイルを選択する(S1)。基準プロファイルは温度推定モデルをビルドアップ(bulid)するために必要な初期値として、温度推定に必要な因子を含む入力データおよび測定された温度である出力データを含み得る。予備検証プロファイルは温度推定に必要なパラメータ個数とパラメータ範囲により決定されるすべてのパラメータセットのうちの温度推定モデルとして適したパラメータセットを決定するために利用されるデータを提供する。再検証プロファイルは温度推定モデルを最終的に検証するためのデータを提供する。
基準プロファイル、予備検証プロファイル、および再検証プロファイルそれぞれは時系列的に区分されて羅列された少なくともサーミスタの温度値、バッテリ電流、バッテリモジュールの冷却水温度値、バッテリモジュールの外気温度値、バッテリモジュール内の最高/最低温度に関するデータを含み得る。この時、サーミスタの温度値、バッテリ電流、バッテリモジュールの冷却水温度値、およびバッテリモジュールの外気温度値は温度推定モデルの入力データであり、バッテリモジュール内の最高温度または最低温度は温度推定モデルの出力データであり得る。
先に、温度推定モデル決定装置1は最高温度を推定するのか、最低温度を推定するのかを選択し(S2)、最高温度推定の場合、最高温度推定モデルを設定し(S3)、最低温度推定の場合、最低温度推定モデルを設定する(S4)。最高温度推定と最低温度推定それぞれのためのモデル設定は同様の方法により行われるので、以下、図3および図4に示す各段階を共に説明する。
パラメータセット生成部10は最高温度または最低温度推定のためのモデルの種類を選択し、選択されたモデル種類によるパラメータ個数および各パラメータの範囲区間数を乗じたn個のすべてのパラメータセットを計算する(S301,S401)。温度推定のためのモデルの種類はARXモデル(Auto Regressive eXternal model)および伝達関数モデル(transfer function model)などがあり得る。
パラメータセット生成部10で計算されたすべてのパラメータセットは予備検証プロファイルと共にモデリング部20に伝達され、モデリング部20はすべてのパラメータセットそれぞれに基づいたモデルに予備検証プロファイルの入力データを代入して最高温度または最低温度を推定する(S302,S402)。
誤差計算部30は予備検証プロファイルの出力データである最高温度または最低温度とモデリング部20の出力(OUT1~OUTn)である推定された最高温度または最低温度の間の二乗平均平方根誤差(Root Mean Square Error,RMSE)を算出する(S303,S403)。出力(OUT1-OUTn)それぞれは一つのパラメータセットによって定義されたモデルに対して、予備検証プロファイルに含まれた時系列的に区分されて羅列されたサーミスタの温度値、バッテリ電流、バッテリモジュールの冷却水温度値、およびバッテリモジュールの外気温度値それぞれを、入力をした時の推定された最高温度または最低温度を含み、誤差計算部30は予備検証プロファイルで入力データに対応する出力データと推定された最高温度または最低温度の間の差を二乗平均した結果を平方根処理して二乗平均平方根誤差(RMSE)を算出することができる。RMSEは推定結果と実際の測定結果の間の誤差を算出するための一例として、発明はこれに限定されるものではない。
モデリング部20は誤差計算部30からRMSEを受信し、最高温度または最低温度それぞれに対する臨界値TH2よりRMSEが小さいかを比較する(S304,S404)。
モデリング部20はRMSEが臨界値TH2より小さいパラメータセットとRMSEを保存し(S305,S405)、すべてのパラメータセットに対する段階S304または段階S404を完了した後、保存されたRMSEのうちの小さい順に上位i個に対応するi個の第1パラメータセットを抽出する(S307,S407)。i個は最大個数として、段階S304または段階S404で臨界値TH2よりRMSEが小さいパラメータセットの個数がi個以下である場合、該当パラメータセットを第1パラメータセットとして抽出する。
モデリング部20は段階S304および段階S404を満足するパラメータセットがないと判断すると(S306,S406)、パラメータセット生成部10に失敗信号(FS1)を伝送してこれを知らせ、パラメータセット生成部10は段階S1により基準プロファイルを変更する。この時、基準プロファイルの変更に応じて予備検証プロファイルおよび再検証プロファイルのうちの少なくとも一つが共に変更され得る。第1パラメータセットは予備検証プロファイルと共に最大誤差計算部40に伝送される。
最大誤差計算部40はi個の第1パラメータセットそれぞれに基づいたモデルを用いて時系列的に変わる推定最高温度または推定最低温度と予備検証プロファイルでの出力データの測定された最高温度または最低温度の間の最大誤差を算出する(S308,S408)。上で言及したように、予備検証プロファイルは時系列的に区分されて羅列された入力データおよび出力データを含み、最大誤差は時系列的に変わる入力データにより推定された最高温度または最低温度と入力データに対応する出力データすなわち、測定された最高温度または最低温度の間の誤差のうちの最も大きい誤差を意味する。
最大誤差計算部40は第1パラメータセットそれぞれの最大誤差と臨界値TH3より小さいかを比較する(S309,S409)。最大誤差計算部40は比較結果最大誤差が臨界値TH3より小さいパラメータセットに対して、誤差計算部30により段階S303またはS403で算出されたRMSEに対応する加重値を乗じた結果と最大誤差に対応する加重値を乗じた結果を足して加重誤差を算出してこれを保存する(S310,S410)。
最大誤差計算部40は段階S309および段階S408を満足するパラメータセットがないと判断すると(S311,S411)、パラメータセット生成部10に失敗信号(FS2)を伝送してこれを知らせ、パラメータセット生成部10は段階S1により基準プロファイルを変更する。この時、基準プロファイルの変更に応じて予備検証プロファイルおよび再検証プロファイルのうちの少なくとも一つが共に変更され得る。
最大誤差計算部40は段階S310またはS410で保存された加重誤差のうちの小さい順にj個に対応する第2パラメータセットを抽出する(S312,S412)。j個は最大個数として、段階S309または段階S409で臨界値TH3より最大誤差が小さいパラメータセットの個数がj個以下である場合、該当パラメータセットを第2パラメータセットとして抽出する。第2パラメータセットは再検証プロファイルと共に最終モデル決定部50に伝送される。
最終モデル決定部50は再検証プロファイルの入力データにより最高温度または最低温度を推定し、再検証プロファイルの出力データである測定された最高温度または最低温度と推定最高温度または推定最低温度の間のRMSEを算出する(S5)。
最終モデル決定部50はRMSEが臨界値TH1より小さいかを比較する(S6)。
段階S6で比較結果、RMSEが臨界値TH1より小さい時、最終モデル決定部50は該当パラメータセットに基づいたモデルを最終モデル(FOUT)に決定し得る。もし、j個の第2パラメータセットの少なくとも二つの第2パラメータセットに対するRMSEが臨界値TH1より小さい時、最も小さいRMSEに対応する第2パラメータセットを選択し得る。
最終モデル(FOUT)は最高温度に対する最高温度推定モデルおよび最低温度については最低温度推定モデルに区分され、バッテリ管理システム(Battery Management System)に挿入され得る(S7)。
図5は一実施形態による温度推定モデルが適用されたバッテリ管理システムを含むバッテリシステムを示す図である。
図5に示すように、バッテリシステム2は、バッテリ管理システム(Battery Management System,BMS,100)、バッテリモジュール200、少なくとも二つの温度センサ300,500、電流センサ400、および冷却部600を含む。
BMS100は先立って図1から図4で説明した実施形態によって決定された温度推定モデル110を含む。温度推定モデル110はバッテリモジュール200の最高温度および最低温度のうちの少なくとも一つを推定し得る。
バッテリモジュール200は複数のバッテリセルが直列/並列連結されて必要な電源を供給し得る。図5では、バッテリモジュール200が直列連結されている複数のバッテリセル(C1-Cn)を含み、バッテリシステム2の二つの出力端(+、-)の間に連結されていることを示しているが、図5に示す構成および構成間の連結関係は一例であり、発明はこれに限定されるものではない。
温度センサ300はバッテリモジュール200の温度を測定し、測定された温度を指示する情報をBMS100に伝送する。温度センサ300はサーミスタとして具現され得、サーミスタによって測定された値が測定された温度を指示する情報であり得る。図5ではバッテリモジュール200に付着した一つの温度センサ300が示されているが、発明はこれに限定されず、二つ以上の温度センサがバッテリモジュール200に付着してもよい。
電流センサ400はバッテリモジュール200の正極と出力端(+)の間に直列連結され、バッテリモジュール200に流れる電流を測定し、測定された電流を指示する情報をBMS100に伝送する。
温度センサ500は外気温度を測定し、測定された外気温度を指示する情報をBMS100に伝送する。
冷却部600はバッテリモジュール200の少なくとも一面に接触して位置し、バッテリモジュール200の温度を低くし得る。冷却部600の内部には冷却水があり、冷却部600は冷却水の温度を測定し、測定された冷却水温度を指示する情報をBMS100に伝送する。
BMS100に伝送されたサーミスタ測定値、バッテリに流れる電流、外気温度値、および冷却水温度値は温度推定モデル110に入力され、温度推定モデル110はバッテリモジュール200の最高温度および最低温度のうちの少なくとも一つを推定し得る。
上記実施形態により、バッテリモジュール内の最高温度および最低温度を正確に推定することができ、サーミスタの測定値にのみ依存してバッテリモジュールの温度を測定する従来技術の問題を克服することができる。
また、バッテリモジュール内の最高温度および最低温度を正確に推定することができ、バッテリモジュールの出力制御をより精巧に遂行できるバッテリ管理システムが提供される。
以上、本発明の実施形態について詳細に説明したが、本発明の権利範囲がこれに限定されるのではなく、本発明の属する分野における通常の知識を持った者が様々に変形および改良した形態も本発明の権利範囲に属する。

Claims (14)

  1. バッテリモジュール内の温度を推定する温度推定モデルを決定する温度推定モデル決定方法であって、
    温度推定のために選択されたモデル種類によるn個(nは自然数)のパラメータおよび前記n個のパラメータそれぞれの各パラメータ範囲区間数に基づいてすべてのパラメータセットを計算する段階と、
    温度推定に用いられる入力データと、前記バッテリモジュール内の測定温度とを含み、前記温度推定モデルの予備検証に用いられる予備検証プロファイルを取得する段階と、
    温度推定に用いられる入力データと、前記バッテリモジュール内の測定温度とを含み、前記温度推定モデルの再検証に用いられる再検証プロファイルを取得する段階と、
    すべての前記パラメータセットのそれぞれに対応するモデルに、前記予備検証プロファイルの前記入力データを代入して温度を推定する段階と、
    すべての前記パラメータセットのそれぞれに対して、前記予備検証プロファイルの前記温度と、前記推定された温度との差による第1誤差を算出する段階と、
    すべての前記パラメータセットに対する第1誤差のそれぞれと、第1臨界値とを比較する段階と、
    前記第1誤差のうち、前記第1臨界値より小さい第1誤差に対応するパラメータセットを第1パラメータセットとして抽出する段階と、
    前記第1パラメータセットに対応するモデルを用いて前記予備検証プロファイルの時系列的に変わる入力データにより推定された温度と、前記予備検証プロファイルの前記温度との差のうちの最も大きい差である最大誤差を算出する段階と、
    前記最大誤差と第2臨界値とを比較する段階と、
    前記第2臨界値より小さい前記最大誤差に対応する加重値を乗じた値と、前記第1パラメータセットに応じた第1誤差に対応する加重値を乗じた値を足した結果に基づいて第2パラメータセットを抽出する段階と、
    前記第2パラメータセットに対応するモデルに、前記再検証プロファイルの前記入力データを代入して温度を推定する段階と、
    前記再検証プロファイルの前記温度と、前記推定された温度との差による第2誤差を算出する段階と、
    前記第2誤差が第3臨界値より小さい場合、前記第2パラメータセットに対応するモデルを前記温度推定モデルに決定する段階と、を含む、
    温度推定モデル決定方法。
  2. 前記予備検証プロファイルおよび前記再検証プロファイルは、
    少なくとも時系列的に変わる温度センサの温度値、バッテリ電流、冷却水温度値、および空気温度値を前記入力データとして含む、請求項1に記載の温度推定モデル決定方法。
  3. 前記推定された温度は、推定された前記バッテリモジュールの最高温度または最低温度を含み、
    前記測温度は、測定された前記バッテリモジュールの最高温度または最低温度を含む、請求項1または2に記載の温度推定モデル決定方法。
  4. 前記第1誤差および前記第2誤差のうちの少なくとも一つは、二乗平均平方根誤差である、請求項1から3のいずれか一項に記載の温度推定モデル決定方法。
  5. 前記第1誤差がすべて前記第1臨界値以上であれば、基準プロファイルを変更する段階をさらに含み、
    前記基準プロファイルは、前記温度推定モデルをビルドアップするための初期値として、温度推定に必要な因子および測定された温度を含む、請求項1から4のいずれか一項に記載の温度推定モデル決定方法。
  6. 前記第1パラメータセットすべての前記最大誤差が前記第2臨界値以上であれば、前記基準プロファイルを変更する段階をさらに含む、請求項5に記載の温度推定モデル決定方法。
  7. 前記第2パラメータセットすべての前記第2誤差が前記第3臨界値以上であれば、前記基準プロファイルを変更する段階をさらに含む、請求項6に記載の温度推定モデル決定方法。
  8. 前記第1パラメータセットを抽出する段階は、
    前記第1誤差のうちの小さい順に上位i個に属する第1誤差に対応するパラメータセットを前記第1パラメータセットとして抽出する段階を含む、請求項1から7のいずれか一項に記載の温度推定モデル決定方法。
  9. 前記第2パラメータセットを抽出する段階は、
    前記足した結果が小さい順に上位j個(jは自然数)に属する足した結果に対応するパラメータセットを前記第2パラメータセットとして抽出する段階を含む、請求項1から8のいずれか一項に記載の温度推定モデル決定方法。
  10. バッテリモジュール内の温度を推定する温度推定モデルを決定する温度推定モデル決定装置であって、
    温度推定のために選択されたモデル種類によるn個(nは自然数)のパラメータおよび前記n個のパラメータのそれぞれの各パラメータ範囲区間数に基づいてすべてのパラメータセットを計算するパラメータセット生成部であって温度推定に用いられる入力データと、前記バッテリモジュール内の測定温度とを含み、前記温度推定モデルの予備検証に用いられる予備検証プロファイルを取得し、温度推定に用いられる入力データと、前記バッテリモジュール内の測定温度とを含み、前記温度推定モデルの再検証に用いられる再検証プロファイルを取得する、パラメータ生成部と、
    すべての前記パラメータセットのそれぞれに対応するモデルに予備検証プロファイルの前記入力データを代入して温度を推定するモデリング部と、
    すべての前記パラメータセットのそれぞれに対して、前記予備検証プロファイルの前記温度と、前記推定された温度との差による第1誤差を算出する誤差計算部と、を含み、
    前記モデリング部は、
    すべての前記パラメータセットに対する第1誤差のそれぞれと、第1臨界値とを比較し、前記第1誤差のうちの前記第1臨界値より小さい第1誤差に対応するパラメータセットを第1パラメータセットとして抽出し、
    前記第1パラメータセットに対応するモデルを用いて前記予備検証プロファイルの時系列的に変わる入力データにより推定された温度と、前記予備検証プロファイルの前記温度との差のうちの最も大きい差である最大誤差を算出し、前記最大誤差と、第2臨界値とを比較し、前記第2臨界値より小さい前記最大誤差に対応する加重値を乗じた値と、前記第1パラメータセットに応じた第1誤差に対応する加重値を乗じた値を足した結果に基づいて第2パラメータセットを抽出する最大誤差計算部と、
    前記第2パラメータセットに対応するモデルに、前記再検証プロファイルの前記入力データを代入して温度を推定し、前記再検証プロファイルの前記温度と、前記推定された温度との差による第2誤差を算出し、前記第2誤差が第3臨界値より小さい場合、前記第2パラメータセットに対応するモデルを前記温度推定モデルに決定する最終モデル決定部と、を有する、温度推定モデル決定装置。
  11. 前記予備検証プロファイルおよび前記再検証プロファイルは、
    少なくとも時系列的に変わる温度センサの温度値、バッテリ電流、冷却水温度値、および空気温度値を入力データとして含む、請求項10に記載の温度推定モデル決定装置。
  12. 前記推定された温度は、推定された前記バッテリモジュールの最高温度または最低温度を含み、
    前記測温度は、測定された前記バッテリモジュールの最高温度または最低温度を含む、請求項10または11に記載の温度推定モデル決定装置。
  13. 前記第1誤差および前記第2誤差のうちの少なくとも一つは、二乗平均平方根誤差である、請求項10から12のいずれか一項に記載の温度推定モデル決定装置。
  14. バッテリモジュール内の温度を推定する温度推定モデル決定方法によって決定された温度推定モデルを有し、
    前記温度推定モデル決定方法は、
    温度推定のために選択されたモデル種類によるn個(nは自然数)のパラメータおよび前記n個のパラメータのそれぞれの各パラメータ範囲区間数に基づいてすべてのパラメータセットを計算する段階と、
    温度推定に用いられる入力データと、前記バッテリモジュール内の測定温度とを含み、前記温度推定モデルの予備検証に用いられる予備検証プロファイルを取得する段階と
    温度推定に用いられる入力データと、前記バッテリモジュール内の測定温度とを含み、前記温度推定モデルの再検証に用いられる再検証プロファイルを取得する段階と
    すべての前記パラメータセットのそれぞれに対応するモデルに前記予備検証プロファイルの前記入力データを代入して温度を推定する段階と、
    すべての前記パラメータセットのそれぞれに対して、前記予備検証プロファイルの前記温度と、前記推定された温度との差による第1誤差を算出する段階と、
    すべての前記パラメータセットに対する第1誤差のそれぞれと、第1臨界値とを比較する段階と、
    前記第1誤差のうちの前記第1臨界値より小さい第1誤差に対応するパラメータセットを第1パラメータセットとして抽出する段階と、
    前記第1パラメータセットに対応するモデルを用いて前記予備検証プロファイルの時系列的に変わる入力データにより推定された温度と、前記予備検証プロファイルの前記温度との差のうちの最も大きい差である最大誤差を算出する段階と、
    前記最大誤差と、第2臨界値とを比較する段階と、
    前記第2臨界値より小さい前記最大誤差に対応する加重値を乗じた値と、前記第1パラメータセットに応じた第1誤差に対応する加重値を乗じた値を足した結果に基づいて第2パラメータセットを抽出する段階と、
    前記第2パラメータセットに対応するモデルに、前記再検証プロファイルの前記入力データを代入して温度を推定する段階と、
    前記再検証プロファイルの前記温度と、前記推定された温度との差による第2誤差を算出する段階と、
    前記第2誤差が第3臨界値より小さい場合、前記第2パラメータセットに対応するモデルを前記温度推定モデルに決定する段階と、を含む、
    バッテリ管理システム。
JP2021520343A 2019-08-29 2020-08-26 温度推定モデル決定方法および装置、温度推定モデルが適用されたバッテリ管理システム Active JP7167415B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
KR10-2019-0106742 2019-08-29
KR1020190106742A KR102634916B1 (ko) 2019-08-29 2019-08-29 온도 추정 모델 결정 방법 및 장치, 온도 추정 모델이 적용된 배터리 관리 시스템
PCT/KR2020/011375 WO2021040396A1 (ko) 2019-08-29 2020-08-26 온도 추정 모델 결정 방법 및 장치, 온도 추정 모델이 적용된 배터리 관리 시스템

Publications (2)

Publication Number Publication Date
JP2022504862A JP2022504862A (ja) 2022-01-13
JP7167415B2 true JP7167415B2 (ja) 2022-11-09

Family

ID=74685595

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021520343A Active JP7167415B2 (ja) 2019-08-29 2020-08-26 温度推定モデル決定方法および装置、温度推定モデルが適用されたバッテリ管理システム

Country Status (6)

Country Link
US (1) US11754630B2 (ja)
EP (1) EP3855196A4 (ja)
JP (1) JP7167415B2 (ja)
KR (1) KR102634916B1 (ja)
CN (1) CN113631937B (ja)
WO (1) WO2021040396A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113899527B (zh) * 2021-12-06 2022-03-01 中国空气动力研究与发展中心低速空气动力研究所 一种试验模型表面温度的修正方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012103220A (ja) 2010-11-15 2012-05-31 Honda Motor Co Ltd 電池温度検出装置および電池温度検出方法
CN102494778A (zh) 2011-11-14 2012-06-13 北京理工大学 一种基于人工神经网络的二次电池表面最高温度预测方法
JP2014524629A (ja) 2011-08-19 2014-09-22 ハートフォード スチーム ボイラー インスペクション アンド インシュアランス カンパニー 動的外れ値偏り低減システム及び方法
WO2016038658A1 (ja) 2014-09-08 2016-03-17 株式会社 東芝 電池パック、制御回路および制御方法
CN109987000A (zh) 2019-02-19 2019-07-09 中国第一汽车股份有限公司 一种动力电池温度预测***及方法
JP6555440B1 (ja) 2018-03-20 2019-08-07 株式会社Gsユアサ 異常要因判定装置、劣化判定装置、コンピュータプログラム、劣化判定方法及び異常要因判定方法

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6076964A (en) * 1998-11-11 2000-06-20 Chrysler Corporation Prediction of internal temperature of a battery using a non-linear dynamic model
JP3815218B2 (ja) 2000-12-20 2006-08-30 株式会社日立製作所 データ分析方法及び装置
JP4664581B2 (ja) 2003-08-27 2011-04-06 三菱重工業株式会社 燃料電池システム
EP2211335A1 (en) 2009-01-21 2010-07-28 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Apparatus, method and computer program for obtaining a parameter describing a variation of a signal characteristic of a signal
US8404572B2 (en) 2009-02-13 2013-03-26 Taiwan Semiconductor Manufacturing Co., Ltd Multi-zone temperature control for semiconductor wafer
US9111212B2 (en) 2011-08-19 2015-08-18 Hartford Steam Boiler Inspection And Insurance Company Dynamic outlier bias reduction system and method
KR100934410B1 (ko) 2009-04-29 2009-12-29 강원대학교산학협력단 근적외 분광분석법을 이용한 곡물 종실의 무게 측정방법
JP2013012617A (ja) 2011-06-30 2013-01-17 Sharp Corp 気相成長装置、気相成長装置の温度制御方法および温度制御プログラム
US10557840B2 (en) 2011-08-19 2020-02-11 Hartford Steam Boiler Inspection And Insurance Company System and method for performing industrial processes across facilities
PL2734977T3 (pl) * 2012-02-22 2016-02-29 Aselsan Elektronik Sanayi Ve Ticaret Anonim Sirketi System i sposób optymalizacji systemu śledzenia
KR101977729B1 (ko) 2012-11-26 2019-05-14 한국전자통신연구원 배터리 관리 시스템 알고리즘 최적화 방법 및 이를 위한 장치
EP3514700A1 (en) 2013-02-20 2019-07-24 Hartford Steam Boiler Inspection and Insurance Company Dynamic outlier bias reduction system and method
US10496927B2 (en) 2014-05-23 2019-12-03 DataRobot, Inc. Systems for time-series predictive data analytics, and related methods and apparatus
US10366346B2 (en) 2014-05-23 2019-07-30 DataRobot, Inc. Systems and techniques for determining the predictive value of a feature
US10558924B2 (en) 2014-05-23 2020-02-11 DataRobot, Inc. Systems for second-order predictive data analytics, and related methods and apparatus
DE112015002433T5 (de) 2014-05-23 2017-03-23 Datarobot Systeme und Techniken zur prädikativen Datenanalytik
KR101827108B1 (ko) 2016-05-04 2018-02-07 두산중공업 주식회사 플랜트 이상 감지 학습 시스템 및 방법
EP3488512A1 (en) 2016-07-22 2019-05-29 Eos Energy Storage, LLC Battery management system
GB2571651B (en) 2016-10-21 2022-09-21 Datarobot Inc Systems for predictive data analytics, and related methods and apparatus
JP7506967B2 (ja) 2018-03-06 2024-06-27 ナブテスコ株式会社 歯車ユニット
CN108535661B (zh) * 2018-05-18 2020-03-31 东北大学 一种基于模型误差谱的动力电池健康状态在线估计方法
WO2020246631A1 (ko) 2019-06-04 2020-12-10 엘지전자 주식회사 온도 예측 모델의 생성 장치 및 시뮬레이션 환경의 제공 방법

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012103220A (ja) 2010-11-15 2012-05-31 Honda Motor Co Ltd 電池温度検出装置および電池温度検出方法
JP2014524629A (ja) 2011-08-19 2014-09-22 ハートフォード スチーム ボイラー インスペクション アンド インシュアランス カンパニー 動的外れ値偏り低減システム及び方法
CN102494778A (zh) 2011-11-14 2012-06-13 北京理工大学 一种基于人工神经网络的二次电池表面最高温度预测方法
WO2016038658A1 (ja) 2014-09-08 2016-03-17 株式会社 東芝 電池パック、制御回路および制御方法
JP6555440B1 (ja) 2018-03-20 2019-08-07 株式会社Gsユアサ 異常要因判定装置、劣化判定装置、コンピュータプログラム、劣化判定方法及び異常要因判定方法
CN109987000A (zh) 2019-02-19 2019-07-09 中国第一汽车股份有限公司 一种动力电池温度预测***及方法

Also Published As

Publication number Publication date
US20220043065A1 (en) 2022-02-10
WO2021040396A1 (ko) 2021-03-04
JP2022504862A (ja) 2022-01-13
EP3855196A4 (en) 2022-02-23
US11754630B2 (en) 2023-09-12
CN113631937B (zh) 2023-07-18
EP3855196A1 (en) 2021-07-28
CN113631937A (zh) 2021-11-09
KR102634916B1 (ko) 2024-02-06
KR20210026238A (ko) 2021-03-10

Similar Documents

Publication Publication Date Title
CN105588667B (zh) 一种高精密热敏电阻温度计校准装置
CN100483146C (zh) 利用基于健康指标的存留寿命模型的自适应改变的状态来估计电池存留寿命的设备与方法
CN106291381B (zh) 一种联合估计动力电池***荷电状态与健康状态的方法
JP6509725B2 (ja) バッテリの充電状態の推定
US11527780B2 (en) Battery model and control application calibration systems and methods
CN105988086A (zh) 电池余量预测装置及电池组
JP7167415B2 (ja) 温度推定モデル決定方法および装置、温度推定モデルが適用されたバッテリ管理システム
CN111542760B (zh) 用于校正分流电阻器的电流值的***和方法
CN103884473B (zh) 一种水冷***的漏水判断方法
CN115828699B (zh) 功率半导体模块全生命周期结温预测方法、***及终端
CN109282502A (zh) 一种储水式热水器的剩余洗浴时间的预报方法
CN113919222A (zh) 一种电池组的内部温度在线计算方法
CN115158076A (zh) 计量误差评估方法、装置及计算机可读存储介质
CN113933734A (zh) 一种退役电池包内部单体参数的提取方法
CN110412472A (zh) 一种基于正态伽马滤波的电池荷电状态估计方法
CN110377941B (zh) 建立卫星温度场的罚盲似然克里金代理模型的方法
JP2008164469A (ja) サーミスタの短絡故障検出装置
CN111122189B (zh) 燃气设备的测试台位不确定度测定方法及热效率测定方法
CN109858187B (zh) 基于正则化方式建立卫星温度场的克里金代理模型的方法
CN114580214B (zh) 基于物联网的智能型蒸汽流量计校核测试***及方法
CN103003674A (zh) 用于确定在废气物料流量传感器上得到的总物料流量的方法
CN112213642B (zh) 一种并联电池组中单体连接断开的检测及定位方法
KR101511289B1 (ko) 서셉터 온도 균일도 평가 방법, 장치 및 시스템
CN209878979U (zh) 电流测试验证的试验***
CN116662763B (zh) 一种用于变频器温升测试的数据处理方法和***

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210415

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220511

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220524

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220817

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221007

R150 Certificate of patent or registration of utility model

Ref document number: 7167415

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150