JP7385696B2 - 測定装置、測定方法、及びプログラム - Google Patents
測定装置、測定方法、及びプログラム Download PDFInfo
- Publication number
- JP7385696B2 JP7385696B2 JP2022046656A JP2022046656A JP7385696B2 JP 7385696 B2 JP7385696 B2 JP 7385696B2 JP 2022046656 A JP2022046656 A JP 2022046656A JP 2022046656 A JP2022046656 A JP 2022046656A JP 7385696 B2 JP7385696 B2 JP 7385696B2
- Authority
- JP
- Japan
- Prior art keywords
- bode diagram
- impedance
- vertices
- vertex
- parallel circuits
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 25
- 238000010586 diagram Methods 0.000 claims description 101
- 239000000284 extract Substances 0.000 claims description 35
- 238000005259 measurement Methods 0.000 claims description 28
- 238000009792 diffusion process Methods 0.000 claims description 6
- 238000002847 impedance measurement Methods 0.000 description 15
- 238000000691 measurement method Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000003775 Density Functional Theory Methods 0.000 description 2
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000010365 information processing Effects 0.000 description 2
- 229910001416 lithium ion Inorganic materials 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000003990 capacitor Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000001453 impedance spectrum Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/389—Measuring internal impedance, internal conductance or related variables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/08—Measuring resistance by measuring both voltage and current
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
- G01R31/3842—Arrangements for monitoring battery or accumulator variables, e.g. SoC combining voltage and current measurements
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Tests Of Electric Status Of Batteries (AREA)
- Secondary Cells (AREA)
Description
本発明は、測定装置、測定方法、及びプログラムに関する。
近年、リチウムイオン二次電池等の二次電池が搭載された電動車両の普及が進んでいる。二次電池は、充放電の繰り返しとともに、又は時間の経過とともに劣化する。二次電池の劣化の進行度合いを診断する手法として、交流インピーダンス測定法が公知である。
交流インピーダンス測定法では、電極系に印加した交流電圧(電流)に対する応答電流(電圧)を測定し、その比がインピーダンスとなる。その表記法としては、ボード線図とナイキスト図がよく用いられる。ボード線図は、印加した交流電圧の周波数f(Hz)の対数を横軸に、インピーダンスの絶対値|Z|(Ω)の対数と位相差θ(deg)を、それぞれ縦軸にプロットした図である。|Z|は印加した交流電圧と応答電流の振幅の比、θは電圧と電流の位相の差である。ナイキスト図では、横軸にインピーダンスの実数部を、縦軸に虚数部をとりインピーダンスを複素平面上に表記する。実数部Re[Z](抵抗)は、式:Re[Z]=|Z|cosθにより、虚数部Im[Z](リアクタンス)は、式:Im[Z]=|Z|sinθにより、それぞれ計算できる。
特許文献1には、複数の二次電池を含む電池モジュールの特性に関する情報を処理する電池情報処理システムであって、前記電池モジュールの交流インピーダンスの測定結果を解析する解析装置と、前記交流インピーダンスを表現する等価回路モデルに含まれる複数の回路定数と、前記特性との間の相関関係が記憶された記憶装置とを備え、前記解析装置は、前記交流インピーダンスの実数成分に関するボード線図である第1の周波数特性図上と、前記交流インピーダンスの虚数成分に関するボード線図である第2の周波数特性図上とに、前記交流インピーダンスの測定結果をプロットし、前記第1の周波数特性図上のプロット結果に対するフィッティング処理により第1の多項式曲線を取得するとともに、前記第2の周波数特性図上のプロット結果に対するフィッティング処理により第2の多項式曲線を取得し、前記第1及び第2の多項式曲線を、前記交流インピーダンスの実数成分と虚数成分とに関するナイキスト図上のインピーダンス曲線に変換し、前記インピーダンス曲線から前記複数の回路定数を抽出し、前記相関関係を参照することによって、抽出された前記複数の回路定数から前記特性を評価する、電池情報処理システムが記載されている。
特許文献2には、複数個の電池セルが直列に接続され、実負荷を高周波域を含む負荷変動を生じる状態で駆動する電池モジュールをリアルタイムで測定監視する電池監視装置に用いられるインピーダンス測定装置であって、前記各セルの電圧波形データ及び電流波形データを離散フーリエ変換し、電圧波形データの離散フーリエ変換結果を電流波形データの離散フーリエ変換結果で除算することによりインピーダンスを演算するDFT演算部と、このDFT演算部で演算されたインピーダンスデータの特徴に基づいて最適な等価回路モデルを選択する回路モデル選択部と、この回路モデル選択部で選択された等価回路モデルにおいて定数フィッティングを行う回路定数推定演算部、とで構成されている電池インピーダンス測定装置が記載され、さらに前記電池インピーダンス測定装置において、前記回路モデル選択部は、等価回路モデルの具体的な構成について、1)ワールブルグ素子の有無、2)LR並列回路の有無、3)RC並列回路の段数の手順で順次決定することが記載されている。
交流インピーダンス測定では、交流の入力信号(電圧又は電流)に対する出力信号(電圧又は電流)の比|Z|と位相遅れθを測定する。入力信号の周波数fを変えることによって異なる律速過程の|Z|及びθが得られる。すなわち、(1)微弱な交流信号を入力する、(2)応答を時間場から周波数場へ高速フーリエ変換する、(3)出力/入力インピーダンス比|Z|及び位相遅れθを算出する(ボード線図)、(4)インピーダンス実部、虚部を複素平面上に表示する(ナイキスト図)、という手順で測定が行われる。インピーダンススペクトルは、横軸に実部(Z’)、縦軸に虚部(Z”)をプロットしたナイキスト図で表現される。二次電池の等価回路として、抵抗RとキャパシタCを並列接続したR/C並列回路が用いられる。二次電池の交流インピーダンス測定においては、測定対象の二次電池の等価回路を決める必要があるが、理想的には、1つのR/C並列回路が、ナイキスト図上で独立した半円弧として現れる。しかし、実際の二次電池では複数の円弧が重なり合って現れるため、正しく分離することが難しい。特に、他と比較して小さい円弧成分は目視確認が難しく、詳細な抵抗成分を見落としていた。すなわち、R/C並列回路で表現される円弧成分が何個あるかを見極めるのが困難であった。
本発明は、このような事情を考慮してなされたものであり、二次電池の交流インピーダンス測定において、等価回路におけるR/C並列回路数をより精確にカウントできる、測定装置、測定方法、及びプログラムを提供することを目的の一つとする。
この発明に係る測定装置、測定方法、及びプログラムは、以下の構成を採用した。
(1):この発明の一態様に係る測定装置は、インピーダンス演算器と、前記インピーダンス演算器と信号を入出力するための信号入出力インターフェイスと、を備え、前記インピーダンス演算器は次の(S1)~(S5)の処理を実行する測定装置である。
(S1)前記インピーダンス演算器は、横軸に周波数f(Hz)、縦軸に位相θ(度)をプロットしてボード線図を作成する。
(S2)前記インピーダンス演算器は、作成したボード線図上で頂点を抽出して、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングする。
(S3)前記インピーダンス演算器は、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる。
(S4)前記インピーダンス演算器が、ボード線図上で頂点を抽出しなくなればR/C並列回路数の同定を終了するが、頂点を抽出すれば、前記(S2)及び前記(S3)を繰り返す。ただし、頂点の位相(θ)の大きさは、予め設定した位相(θ)の大きさが予め定めたθ1以上である。
(S5)前記インピーダンス演算器が、R/C並列回路数の総カウント数を報告する。
ただし、(S1)において、前記インピーダンス演算器は、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、(S2)において、測定ノイズと頂点とを見分けるため、頂点の位相(θ)の大きさの下限値(θ0)を予め設定しておき、前記インピーダンス演算器は、θ値がθ0以上である頂点候補から頂点を抽出する。
さらに、ステップS4において、頂点の位相(θ)の大きさの下限値(θ1)を予め設定しておき、前記インピーダンス演算器は、θ値がθ1以上である頂点候補を頂点として抽出する。
(1):この発明の一態様に係る測定装置は、インピーダンス演算器と、前記インピーダンス演算器と信号を入出力するための信号入出力インターフェイスと、を備え、前記インピーダンス演算器は次の(S1)~(S5)の処理を実行する測定装置である。
(S1)前記インピーダンス演算器は、横軸に周波数f(Hz)、縦軸に位相θ(度)をプロットしてボード線図を作成する。
(S2)前記インピーダンス演算器は、作成したボード線図上で頂点を抽出して、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングする。
(S3)前記インピーダンス演算器は、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる。
(S4)前記インピーダンス演算器が、ボード線図上で頂点を抽出しなくなればR/C並列回路数の同定を終了するが、頂点を抽出すれば、前記(S2)及び前記(S3)を繰り返す。ただし、頂点の位相(θ)の大きさは、予め設定した位相(θ)の大きさが予め定めたθ1以上である。
(S5)前記インピーダンス演算器が、R/C並列回路数の総カウント数を報告する。
ただし、(S1)において、前記インピーダンス演算器は、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、(S2)において、測定ノイズと頂点とを見分けるため、頂点の位相(θ)の大きさの下限値(θ0)を予め設定しておき、前記インピーダンス演算器は、θ値がθ0以上である頂点候補から頂点を抽出する。
さらに、ステップS4において、頂点の位相(θ)の大きさの下限値(θ1)を予め設定しておき、前記インピーダンス演算器は、θ値がθ1以上である頂点候補を頂点として抽出する。
(2):この発明の一態様に係る測定方法は、
プロセッサが、横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成するステップS1と、
プロセッサが、作成したボード線図上で頂点を抽出して、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングするステップ2と、
プロセッサが、フィッティングしたR/C並列回路で表現できるボード線図成分を、作成したボード線図から減ずるステップS3と、
プロセッサが、ボード線図上で頂点を抽出しなければステップS5に進み、さらに頂点を抽出すれば、ステップS2及びステップS3を繰り返すように処理を選択するステップS4と、
プロセッサが、R/C並列回路数の総カウント数を報告して、R/C並列回路数の同定を終了するステップS5と
を含む。
ただし、ステップS1において、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、S2において、測定ノイズと頂点とを見分けるため、頂点の位相(θ)の大きさの下限値(θ0)を予め設定しておき、プロセッサが、θ値がθ0以上である頂点候補から頂点を抽出する。
さらに、ステップS4において、頂点の位相(θ)の大きさの下限値(θ1)を予め設定しておき、プロセッサが、θ値がθ1以上である頂点候補を頂点として抽出する。
プロセッサが、横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成するステップS1と、
プロセッサが、作成したボード線図上で頂点を抽出して、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングするステップ2と、
プロセッサが、フィッティングしたR/C並列回路で表現できるボード線図成分を、作成したボード線図から減ずるステップS3と、
プロセッサが、ボード線図上で頂点を抽出しなければステップS5に進み、さらに頂点を抽出すれば、ステップS2及びステップS3を繰り返すように処理を選択するステップS4と、
プロセッサが、R/C並列回路数の総カウント数を報告して、R/C並列回路数の同定を終了するステップS5と
を含む。
ただし、ステップS1において、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、S2において、測定ノイズと頂点とを見分けるため、頂点の位相(θ)の大きさの下限値(θ0)を予め設定しておき、プロセッサが、θ値がθ0以上である頂点候補から頂点を抽出する。
さらに、ステップS4において、頂点の位相(θ)の大きさの下限値(θ1)を予め設定しておき、プロセッサが、θ値がθ1以上である頂点候補を頂点として抽出する。
(3):この発明の一態様に係るプログラムは、
インピーダンス測定装置が備えるインピーダンス演算器に搭載されたプロセッサが、
横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成して、作成したボード線図上で頂点を抽出し、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングし、
フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じ、
ボード線図上で頂点を抽出しなくなればR/C並列回路の総カウント数を報告してR/C並列回路数の同定を終了し、頂点をさらに抽出すれば、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングし、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じることを繰り返す、
プログラムである。
インピーダンス測定装置が備えるインピーダンス演算器に搭載されたプロセッサが、
横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成して、作成したボード線図上で頂点を抽出し、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングし、
フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じ、
ボード線図上で頂点を抽出しなくなればR/C並列回路の総カウント数を報告してR/C並列回路数の同定を終了し、頂点をさらに抽出すれば、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングし、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じることを繰り返す、
プログラムである。
(1)~(3)の態様によれば、合成インピーダンスを表すボード線図では区別できなかった円弧成分を区別することが可能となり、等価回路におけるR/C並列回路数をより精確にカウントできるようになる。
以下、図面を参照し、本発明の測定装置、測定方法、及びプログラムの実施形態について説明する。
本発明において、二次電池は特に限定されないが、リチウムイオン二次電池が好ましい態様である。
本発明において、二次電池は特に限定されないが、リチウムイオン二次電池が好ましい態様である。
[測定装置]
本実施形態の測定装置100は、後述するインピーダンス演算器80及びインピーダンス演算器80と信号を入出力するための信号入出力インターフェイス(図示せず)を備える。インピーダンス演算器80は、また、交流電圧モニタ60及び交流電流モニタ70からの信号を入力するためのインターフェイスを備える。測定装置100はインピーダンス演算器80からなるものでもよいが、交流電圧モニタ60、交流電流モニタ、電流センサ12、掃引信号発生器50を統合していてもよい。
本実施形態の測定装置100は、後述するインピーダンス演算器80及びインピーダンス演算器80と信号を入出力するための信号入出力インターフェイス(図示せず)を備える。インピーダンス演算器80は、また、交流電圧モニタ60及び交流電流モニタ70からの信号を入力するためのインターフェイスを備える。測定装置100はインピーダンス演算器80からなるものでもよいが、交流電圧モニタ60、交流電流モニタ、電流センサ12、掃引信号発生器50を統合していてもよい。
測定装置100の使用形態を示す図1では、二次電池10と電流センサ12の直列回路の両端には、掃引信号発生器50が接続されている。掃引信号発生器50は、出力周波数が掃引変化する交流信号を、二次電池10と電流センサ12の直列回路に出力する。
交流電圧モニタ60は、二次電池10の両端の交流電圧を測定してインピーダンス演算器80に入力する。交流電流モニタ70は、電流センサ12に流れる交流電流を測定してインピーダンス演算器80に入力する。
インピーダンス演算器80は、掃引信号発生器50の出力信号の各周波数における交流電圧モニタ60の測定電圧と交流電流モニタ70の測定電流との比である二次電池10の複素インピーダンスを算出する。
インピーダンス演算器80におけるR/C並列回路数の同定は以下のように行われる。
(1)インピーダンス演算器80は、横軸に掃引信号発生器50の出力信号の周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成する(図2)。この際、インダクタンス(L)成分と拡散(例えば、CPE)成分のフィッティングを行い、測定結果のボード線図から除去する。
(1)インピーダンス演算器80は、横軸に掃引信号発生器50の出力信号の周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成する(図2)。この際、インダクタンス(L)成分と拡散(例えば、CPE)成分のフィッティングを行い、測定結果のボード線図から除去する。
(2)次に、インピーダンス演算器80は、作成したボード線図上で極大かつ最大の点(以下「頂点」ともいう。)T1を抽出し、R/C並列回路数として1をカウントする。インピーダンス演算器80は、点T1を中心にR/C並列回路で表すことができる山の形状をフィッティングする(図3)。なお、測定ノイズと極大かつ最大の点とを見分けるため、棄却するθ値(θ0)が設定されており、インピーダンス演算器80は、位相(θ)の値が予め定めたθ0以上である点から極大かつ最大の点を抽出する。
(3)インピーダンス演算器80は、次に、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる(図4)。すなわち、インピーダンス演算器80は、同じ周波数(Hz)の位相(θ)の値を引き算する。
(4)インピーダンス演算器80が、ボード線図上で極大かつ最大の点(頂点)を抽出しなくなればインピーダンス演算器80によるR/C並列回路数の同定を終了するが、まだ、極大かつ最大の点(頂点)を抽出すれば、(2)及び(3)を繰り返す。
ボード線図上の頂点の高さが予め設定した位相(θ)の大きさが予め定めたθ1(例えば、0.5)以下になったところでR/C並列回路数のカウント増加を止める。
ボード線図上の頂点の高さが予め設定した位相(θ)の大きさが予め定めたθ1(例えば、0.5)以下になったところでR/C並列回路数のカウント増加を止める。
(2-2)インピーダンス演算器80は、ボード線図上でさらに頂点T2を抽出し、R/C並列回路数として2をカウントする。頂点T2を中心にR/C並列回路で表すことができる山の形状をフィッティングする(図5)。
(3-2)インピーダンス演算器80は、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる(図6)。
(4-2)インピーダンス演算器80が、ボード線図上で極大かつ最大の点(頂点)を抽出しなくなればインピーダンス演算器80によるR/C並列回路数の同定を終了するが、まだ、極大かつ最大の点を抽出すれば、さらに(2)及び(3)を繰り返す。
(2-3)インピーダンス演算器80は、ボード線図上で頂点T3を抽出し、R/C並列回路数として3をカウントする。頂点T3を中心にR/C並列回路で表すことができる山の形状をフィッティングする(図7)。
(3-3)インピーダンス演算器80は、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる(図8)。
(4-3)インピーダンス演算器80が、ボード線図上で極大かつ最大の点(頂点)を抽出しなくなればインピーダンス演算器80によるR/C並列回路数の同定を終了するが、まだ、極大かつ最大の点を抽出すれば、さらに(2)及び(3)を繰り返す。
図8では、インピーダンス演算器80はさらに頂点を抽出しなくなるので、ここでR/C並列回路数の同定を終了する。
図8では、インピーダンス演算器80はさらに頂点を抽出しなくなるので、ここでR/C並列回路数の同定を終了する。
(5)インピーダンス演算器80は、R/C並列回路数の総カウント数を報告する。本実施形態のインピーダンス演算器80は、図2に示すボード線図から3つの頂点を抽出できる。すなわち、本発明のインピーダンス測定装置100は、図2に示すボード線図から、R/C並列回路数を3と同定することができる。
従来のインピーダンス測定装置によるR/C並列回路数の同定結果を図9に示す。従来のインピーダンス測定装置は、図2と同一のボード線図から、R/C並列回路数を2と同定することができる。
したがって、本発明のインピーダンス測定装置は、従来、合成されたインピーダンスデータでは抽出できず見落としていた円弧成分(R/C並列回路)を抽出することができ、より精確なインピーダンス測定が可能となる。
したがって、本発明のインピーダンス測定装置は、従来、合成されたインピーダンスデータでは抽出できず見落としていた円弧成分(R/C並列回路)を抽出することができ、より精確なインピーダンス測定が可能となる。
本実施形態のインピーダンス測定装置は、インピーダンス演算器及び信号入出力インターフェイスを備える。この構成要素は、例えば、CPU(Central Processing Unit)などのハードウェアプロセッサがプログラム(ソフトウェア)を実行することにより実現される。これらの構成要素のうち一部又は全部は、LSI(Large Scale Integration)やASIC(Application Specific Integrated Circuit)、FPGA(Field-Programmable Gate Array)、GPU(Graphics Processing Unit)などのハードウェア(回路部;circuitryを含む)によって実現されてもよいし、ソフトウェアとハードウェアの協働によって実現されてもよい。プログラムは、予めHDD(Hard Disk Drive)やフラッシュメモリなどの記憶装置(非一過性の記憶媒体を備える記憶装置)に格納されていてもよいし、DVDやCD-ROMなどの着脱可能な記憶媒体(非一過性の記憶媒体)に格納されており、記憶媒体がドライブ装置に装着されることでインストールされてもよい。
上記説明した実施形態は、以下のように表現することができる。
コンピュータによって読み込み可能な命令(computer-readable instructions)を格納する記憶媒体(storage medium)と、
前記記憶媒体に接続されたプロセッサと、を備え、
前記プロセッサは、前記コンピュータによって読み込み可能な命令を実行することにより(the processor executing the computer-readable instructions to:)
横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成し、
作成したボード線図上で極大かつ最大の点(以下「頂点」ともいう。)T1を抽出して、R/C並列回路数として1をカウントし、点T1を中心にR/C並列回路で表すことができる山の形状をフィッティングし、
フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じ、
ボード線図上で極大かつ最大の点(頂点)を抽出しなくなればR/C並列回路数の同定を終了してR/C並列回路数の総カウント数を報告し、さらに極大かつ最大の点(頂点)を抽出すれば、(ステップS2)及び(ステップS3)を繰り返す、
装置。
コンピュータによって読み込み可能な命令(computer-readable instructions)を格納する記憶媒体(storage medium)と、
前記記憶媒体に接続されたプロセッサと、を備え、
前記プロセッサは、前記コンピュータによって読み込み可能な命令を実行することにより(the processor executing the computer-readable instructions to:)
横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成し、
作成したボード線図上で極大かつ最大の点(以下「頂点」ともいう。)T1を抽出して、R/C並列回路数として1をカウントし、点T1を中心にR/C並列回路で表すことができる山の形状をフィッティングし、
フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じ、
ボード線図上で極大かつ最大の点(頂点)を抽出しなくなればR/C並列回路数の同定を終了してR/C並列回路数の総カウント数を報告し、さらに極大かつ最大の点(頂点)を抽出すれば、(ステップS2)及び(ステップS3)を繰り返す、
装置。
[測定方法]
本実施形態の測定方法は、以下のステップS1~S5を備える。
本実施形態の測定方法の説明では、図10に示すフローチャート、図1の構成図、及び図2~4のボード線図を参照する。また、交流電圧モニタ60及び交流電流モニタ70は、測定データをインピーダンス演算器80に入力しているものとする。
本実施形態の測定方法は、以下のステップS1~S5を備える。
本実施形態の測定方法の説明では、図10に示すフローチャート、図1の構成図、及び図2~4のボード線図を参照する。また、交流電圧モニタ60及び交流電流モニタ70は、測定データをインピーダンス演算器80に入力しているものとする。
(ステップS1)インピーダンス演算器80に搭載されたプロセッサ(以下、単に「プロセッサ」という。)が、横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成する(図2)。この際、インダクタンス(L)成分と拡散(例えば、CPE)成分のフィッティングを行い、測定結果のボード線図から除去する。
(ステップS2)プロセッサが、作成したボード線図上で極大かつ最大の点(以下「頂点」ともいう。)T1を抽出し、R/C並列回路数として+1をカウントする。プロセッサが、点T1を中心にR/C並列回路で表すことができる山の形状をフィッティングする(図3)。なお、測定ノイズと極大かつ最大の点とを見分けるため、棄却するθ値(θ0)が設定されており、プロセッサは、θ値がθ0以上である点から極大かつ最大の点(頂点)を抽出する。
(ステップS3)プロセッサが、次に、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる(図4)。すなわち、プロセッサが、同じ周波数(Hz)の位相(θ)の値を引き算する。
(ステップS4)プロセッサが、ボード線図上で極大かつ最大の点(頂点)を抽出しなくなればステップS5に進む。インピーダンス演算器80が、さらに極大かつ最大の点(頂点)を抽出すれば、(ステップS2)及び(ステップS3)を繰り返す。
(ステップS5)プロセッサが、R/C並列回路数の総カウント数を報告して、R/C並列回路数の同定を終了する。
[プログラム]
本実施形態のプログラムは、上述したインピーダンス演算器に搭載されたプロセッサに以下の(S1)~(S5)の処理を命じるプログラムである。
(S1)横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成する。
(S2)作成したボード線図上で頂点を抽出し、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングする。
(S3)フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じる。
(S4)ボード線図上で頂点を抽出しなくなればR/C並列回路の総カウント数を報告してR/C並列回路数の同定を終了し、頂点をさらに抽出すれば、(S2)及び(S3)の処理を繰り返す。
(S5)R/C並列回路数の総カウント数を報告する。
本実施形態のプログラムは、上述したインピーダンス演算器に搭載されたプロセッサに以下の(S1)~(S5)の処理を命じるプログラムである。
(S1)横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成する。
(S2)作成したボード線図上で頂点を抽出し、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングする。
(S3)フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じる。
(S4)ボード線図上で頂点を抽出しなくなればR/C並列回路の総カウント数を報告してR/C並列回路数の同定を終了し、頂点をさらに抽出すれば、(S2)及び(S3)の処理を繰り返す。
(S5)R/C並列回路数の総カウント数を報告する。
以上、本発明を実施するための形態について実施形態を用いて説明したが、本発明はこうした実施形態に何等限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々の変形及び置換を加えることができる。
10…二次電池
12…電流センサ
50…掃引信号発生器
60…交流電圧モニタ
70…交流電流モニタ
80…インピーダンス演算器
100…インピーダンス測定装置
12…電流センサ
50…掃引信号発生器
60…交流電圧モニタ
70…交流電流モニタ
80…インピーダンス演算器
100…インピーダンス測定装置
Claims (3)
- 二次電池の交流インピーダンスの測定装置であって、インピーダンス演算器と、前記インピーダンス演算器と前記インピーダンス演算器に接続された機器と信号を入出力するための信号入出力インターフェイスと、を備え、前記インピーダンス演算器は次の(S1)~(S5)の処理を実行する測定装置。
(S1)前記インピーダンス演算器は、横軸に周波数f(Hz)、縦軸に位相θ(度)をプロットしてボード線図を作成する。ボード線図の作成時に、前記インピーダンス演算器は、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
(S2)前記インピーダンス演算器は、作成したボード線図上で頂点を抽出して、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングする。ただし、測定ノイズと頂点とを見分けるため、棄却するθ値(θ0)が設定されており、前記インピーダンス演算器は、位相(θ)の値が予め定めたθ0以上である点から頂点を抽出する。
(S3)前記インピーダンス演算器は、フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる。
(S4)前記インピーダンス演算器が、ボード線図上で頂点を抽出しなくなればR/C並列回路数の同定を終了するが、頂点を抽出すれば、前記(S2)及び前記(S3)を繰り返す。ただし、頂点の位相(θ)の大きさは、予め設定した位相(θ)の大きさが予め定めたθ1以上である。
(S5)前記インピーダンス演算器が、R/C並列回路数の総カウント数を報告する。 - 二次電池の交流インピーダンスの測定方法であって、
プロセッサが、横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成するステップS1と、
プロセッサが、作成したボード線図上で頂点を抽出して、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングするステップ2と、
プロセッサが、フィッティングしたR/C並列回路で表現できるボード線図成分を、作成したボード線図から減ずるステップS3と、
プロセッサが、ボード線図上で頂点を抽出しなければステップS5に進み、さらに頂点を抽出すれば、ステップS2及びステップS3を繰り返すように処理を選択するステップS4と、
プロセッサが、R/C並列回路数の総カウント数を報告して、R/C並列回路数の同定を終了するステップS5と
を含む、測定方法。
ただし、ステップS1において、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、S2において、測定ノイズと頂点とを見分けるため、頂点の位相(θ)の大きさの下限値(θ0)を予め設定しておき、プロセッサが、θ値がθ0以上である頂点候補から頂点を抽出する。
さらに、ステップS4において、頂点の位相(θ)の大きさの下限値(θ1)を予め設定しておき、プロセッサが、θ値がθ1以上である頂点候補を頂点として抽出する。 - 二次電池の交流インピーダンスの測定のためのプログラムであって、インピーダンス演算器と前記インピーダンス演算器と前記インピーダンス演算器に接続された機器との間で信号を入出力するための信号入出力インターフェイスと、を備えるインピーダンス測定装置において、前記インピーダンス演算器に搭載されたプロセッサにお以下の(S1)~(S5)の処理を命じるプログラム。
(S1)横軸に周波数f(Hz)(対数軸)、縦軸に位相θ(度)をプロットしてボード線図を作成する。
(S2)作成したボード線図上で頂点を抽出し、R/C並列回路数として+1をカウントし、前記頂点を中心にR/C並列回路で表すことができる山の形状をフィッティングする。
(S3)フィッティングしたR/C並列回路で表現できるボード線図成分を、測定結果のボード線図から減じる。
(S4)ボード線図上で頂点を抽出しなくなればR/C並列回路の総カウント数を報告してR/C並列回路数の同定を終了し、頂点をさらに抽出すれば、(S2)及び(S3)の処理を繰り返す。
(S5)R/C並列回路数の総カウント数を報告する。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2022046656A JP7385696B2 (ja) | 2022-03-23 | 2022-03-23 | 測定装置、測定方法、及びプログラム |
CN202310099506.1A CN116804703A (zh) | 2022-03-23 | 2023-02-08 | 测定装置、测定方法以及存储介质 |
US18/113,629 US11977124B2 (en) | 2022-03-23 | 2023-02-24 | Measurement device, measurement method, and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2022046656A JP7385696B2 (ja) | 2022-03-23 | 2022-03-23 | 測定装置、測定方法、及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2023140688A JP2023140688A (ja) | 2023-10-05 |
JP7385696B2 true JP7385696B2 (ja) | 2023-11-22 |
Family
ID=88078657
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2022046656A Active JP7385696B2 (ja) | 2022-03-23 | 2022-03-23 | 測定装置、測定方法、及びプログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US11977124B2 (ja) |
JP (1) | JP7385696B2 (ja) |
CN (1) | CN116804703A (ja) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008522350A (ja) | 2004-11-29 | 2008-06-26 | ハイドロジェニクス コーポレイション | 電気化学電池の故障状態を検出し表示するためのシステムおよび方法 |
CN104914312A (zh) | 2015-06-18 | 2015-09-16 | 哈尔滨工业大学 | 一种计算交流阻抗谱弛豫时间分布的方法 |
JP2015197363A (ja) | 2014-04-01 | 2015-11-09 | プライムアースEvエナジー株式会社 | 電池状態判定方法及び電池状態判定装置 |
JP3207998U (ja) | 2016-10-03 | 2016-12-15 | 株式会社International Wise Laboratory | リファレンス回路および装置 |
WO2020261799A1 (ja) | 2019-06-27 | 2020-12-30 | ヌヴォトンテクノロジージャパン株式会社 | 電池管理回路、電池管理システムおよび電池管理ネットワーク |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5088081B2 (ja) * | 2007-10-12 | 2012-12-05 | 富士通株式会社 | 電池の測定方法及び電池の製造方法 |
JP2011232075A (ja) * | 2010-04-26 | 2011-11-17 | Yokogawa Electric Corp | 電池特性評価装置および電池特性評価方法 |
JP2012079582A (ja) | 2010-10-04 | 2012-04-19 | Nec Corp | 二次電池、二次電池の正極または負極の抵抗測定方法、および、二次電池の正極または負極の劣化程度診断方法 |
JP5278508B2 (ja) | 2011-07-25 | 2013-09-04 | 横河電機株式会社 | 電池劣化判定装置および方法 |
JP5589988B2 (ja) | 2011-07-28 | 2014-09-17 | 横河電機株式会社 | 電池インピーダンス測定装置 |
JP5912884B2 (ja) * | 2012-06-04 | 2016-04-27 | 日置電機株式会社 | 等価回路解析装置および等価回路解析方法 |
JP5950728B2 (ja) * | 2012-06-29 | 2016-07-13 | 日置電機株式会社 | 等価回路解析装置および等価回路解析方法 |
CN106897522B (zh) * | 2017-02-27 | 2023-08-01 | 长安大学 | 基于磷酸铁锂动力电池阻抗谱的复参数电路模型及方法 |
JP7070043B2 (ja) | 2018-04-26 | 2022-05-18 | トヨタ自動車株式会社 | 電池情報処理システム、電池モジュールの特性評価方法、および、組電池の製造方法 |
JP7070044B2 (ja) * | 2018-04-26 | 2022-05-18 | トヨタ自動車株式会社 | 電池情報処理システム、電池モジュールの満充電容量算出方法、電池モジュールの内部抵抗算出方法、および、組電池の製造方法 |
JP7020273B2 (ja) | 2018-04-26 | 2022-02-16 | トヨタ自動車株式会社 | 電池情報処理システム、組電池、電池モジュールの容量算出方法、および、組電池の製造方法 |
JP7070224B2 (ja) * | 2018-08-08 | 2022-05-18 | トヨタ自動車株式会社 | 電池情報処理システムおよび電池情報処理方法 |
CN112327172B (zh) * | 2020-11-30 | 2021-09-03 | 同济大学 | 一种基于弛豫时间分布的锂离子电池建模方法 |
CN112763545B (zh) * | 2020-12-30 | 2022-08-26 | 宁德新能源科技有限公司 | 锂离子电池eis的交流阻抗数据处理与解读的方法及电池测试设备 |
-
2022
- 2022-03-23 JP JP2022046656A patent/JP7385696B2/ja active Active
-
2023
- 2023-02-08 CN CN202310099506.1A patent/CN116804703A/zh active Pending
- 2023-02-24 US US18/113,629 patent/US11977124B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008522350A (ja) | 2004-11-29 | 2008-06-26 | ハイドロジェニクス コーポレイション | 電気化学電池の故障状態を検出し表示するためのシステムおよび方法 |
JP2015197363A (ja) | 2014-04-01 | 2015-11-09 | プライムアースEvエナジー株式会社 | 電池状態判定方法及び電池状態判定装置 |
CN104914312A (zh) | 2015-06-18 | 2015-09-16 | 哈尔滨工业大学 | 一种计算交流阻抗谱弛豫时间分布的方法 |
JP3207998U (ja) | 2016-10-03 | 2016-12-15 | 株式会社International Wise Laboratory | リファレンス回路および装置 |
WO2020261799A1 (ja) | 2019-06-27 | 2020-12-30 | ヌヴォトンテクノロジージャパン株式会社 | 電池管理回路、電池管理システムおよび電池管理ネットワーク |
Also Published As
Publication number | Publication date |
---|---|
US20230305071A1 (en) | 2023-09-28 |
US11977124B2 (en) | 2024-05-07 |
CN116804703A (zh) | 2023-09-26 |
JP2023140688A (ja) | 2023-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6991616B2 (ja) | 電流パルス法による電池診断装置及び電池診断方法 | |
US8762109B2 (en) | Crosstalk compensation in analysis of energy storage devices | |
US8909490B2 (en) | Battery state estimation device and battery state estimation method | |
Das et al. | Performance of a load-immune classifier for robust identification of minor faults in induction motor stator winding | |
JP4835757B2 (ja) | 電池特性評価装置 | |
Schulze et al. | Parameter identification of unsymmetrical transmission lines using fault records obtained from protective relays | |
WO2014054259A1 (ja) | バッテリのパラメータ推定装置およびその推定方法 | |
US9625503B2 (en) | Method and device for analysing the quality of the electrical energy in a three-phase electric network | |
US20130204555A1 (en) | Method and Apparatus for Electrically Locating a Fault in a Cable | |
CN107992665B (zh) | 一种特高压换流站交流滤波器在线故障诊断分析方法 | |
JP5924617B2 (ja) | 等価回路合成方法並びに装置、および回路診断方法 | |
US10209310B2 (en) | Rotary machine diagnostic system | |
CN112098756B (zh) | 电磁兼容问题定位方法、装置、设备及存储介质 | |
JP4087672B2 (ja) | インピダンススペクトルのパターンマッチング技法を用いた単位電池分類方法 | |
EP3012643A1 (en) | Method and apparatus for identifying causes for cable overcurrent | |
JP7385696B2 (ja) | 測定装置、測定方法、及びプログラム | |
Papadopoulos et al. | Online parameter identification and generic modeling derivation of a dynamic load model in distribution grids | |
Bohórquez et al. | One-ended fault location method based on machine learning models | |
JP7275589B2 (ja) | 周波数応答波形生成装置、異常診断装置、周波数応答波形生成方法、異常診断方法 | |
Yousof et al. | Study of transformer winding deformation by frequency response analysis | |
CN109901030B (zh) | 一种电抗器匝间绝缘状态监测方法、***及应用 | |
JP6931830B2 (ja) | インダクタンスの計測方法、計測システム、及びプログラム | |
EP4198532A1 (en) | A computer implemented method and a system for determining electromotive force curves of a battery | |
Cicilio et al. | Dynamic Composite Load Model Priority Placement Based on Electrical Centrality | |
CN110032758A (zh) | 计算电信号的能量的方法、装置和计算机存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20221129 |
|
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: 20231024 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20231110 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7385696 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |