Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7570970B2 - Acquisition processing device, processing method, processing program, and measurement system - Google Patents
[go: Go Back, main page]

JP7570970B2 - Acquisition processing device, processing method, processing program, and measurement system - Google Patents

Acquisition processing device, processing method, processing program, and measurement system Download PDF

Info

Publication number
JP7570970B2
JP7570970B2 JP2021087400A JP2021087400A JP7570970B2 JP 7570970 B2 JP7570970 B2 JP 7570970B2 JP 2021087400 A JP2021087400 A JP 2021087400A JP 2021087400 A JP2021087400 A JP 2021087400A JP 7570970 B2 JP7570970 B2 JP 7570970B2
Authority
JP
Japan
Prior art keywords
component
value
measured
impedance
real
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
JP2021087400A
Other languages
Japanese (ja)
Other versions
JP2022180739A (en
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.)
Hioki EE Corp
Original Assignee
Hioki EE Corp
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 Hioki EE Corp filed Critical Hioki EE Corp
Priority to JP2021087400A priority Critical patent/JP7570970B2/en
Publication of JP2022180739A publication Critical patent/JP2022180739A/en
Application granted granted Critical
Publication of JP7570970B2 publication Critical patent/JP7570970B2/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

  • Measurement Of Resistance Or Impedance (AREA)
  • Tests Of Electric Status Of Batteries (AREA)

Description

本発明は、蓄電デバイスのインピーダンスを測定する取得処理装置、処理方法、処理プログラム、及び測定システムに関する。 The present invention relates to an acquisition processing device, a processing method, a processing program, and a measurement system for measuring the impedance of an electricity storage device.

特許文献1には、電池を測定対象とした等価回路解析装置が開示されている。この等価回路解析装置は、電池の等価回路を用いた等価回路解析を行う。 Patent document 1 discloses an equivalent circuit analysis device that measures a battery. This equivalent circuit analysis device performs equivalent circuit analysis using the equivalent circuit of a battery.

特開2014-10037号公報JP 2014-10037 A

このような装置において、非線形アルゴリズムで等価回路解析を行う場合、特性を合わせ込む際に、どのパラメータを固定するかによって結果が変わる虞がある。このため、残差値が少なくなったとしても、その値が電池の特性を示すものであることを証明することができない。 When performing equivalent circuit analysis using a nonlinear algorithm on such a device, there is a risk that the results will vary depending on which parameters are fixed when fitting the characteristics. Therefore, even if the residual value is reduced, it cannot be proven that the value represents the battery characteristics.

本発明は、上記の問題点に鑑みてなされたものであり、測定対象となる蓄電デバイスの特性を示す測定結果の取得を可能とすることを目的とする。 The present invention was made in consideration of the above problems, and aims to make it possible to obtain measurement results that indicate the characteristics of the electricity storage device being measured.

本発明のある態様の取得処理装置は、蓄電デバイスのインピーダンスを取得処理する取得処理装置である。取得処理装置は、所定の周波数帯域において計測装置を用いて三点以上の周波数で測定された測定対象成分のインピーダンスを取得する取得手段を含む。取得処理装置は、前記取得手段で取得した前記三点以上のインピーダンスの実数成分又は虚数成分の極値を求める決定手段を含む。取得処理装置は、前記三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて前記測定対象成分の予測値を算出する算出手段を含む。取得処理装置は、前記極値が、前記予測値と同じ又は前記予測値よりも前記測定対象成分のインピーダンスに近づくように突出する場合には、前記極値を前記測定対象成分に関する測定値として確定する確定手段を含む。 An acquisition processing device according to one aspect of the present invention is an acquisition processing device that acquires and processes the impedance of an electric storage device. The acquisition processing device includes an acquisition means that acquires the impedance of a measurement target component measured at three or more frequencies using a measuring device in a predetermined frequency band. The acquisition processing device includes a determination means that determines the extreme values of the real or imaginary components of the impedance at the three or more points acquired by the acquisition means. The acquisition processing device includes a calculation means that calculates a predicted value of the measurement target component based on an approximation formula that represents the relationship between the real and imaginary components of the impedance at the three or more points. The acquisition processing device includes a determination means that, when the extreme value protrudes so as to be the same as the predicted value or closer to the impedance of the measurement target component than the predicted value, determines the extreme value as the measured value for the measurement target component.

この態様において、測定で取得した三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて測定対象成分の予測値を算出する。そして、三点以上の実数成分又は虚数成分の極値が、近似式に基づく予測値と同じ又は予測値よりも測定対象成分のインピーダンスに近づくように突出する場合に、極値を測定対象成分に関する測定値として確定する。 In this embodiment, a predicted value of the component to be measured is calculated based on an approximation formula that expresses the relationship between the real and imaginary components of the impedance at three or more points obtained by measurement. Then, when the extreme values of the real or imaginary components at three or more points stand out so as to be the same as the predicted value based on the approximation formula or to be closer to the impedance of the component to be measured than the predicted value, the extreme values are determined as the measured value for the component to be measured.

そして、この測定値は、実際に測定された実測値によって定められており、この実測値は、蓄電デバイスの特性を示す値である。 This measurement value is determined based on an actual measurement, and this actual measurement value indicates the characteristics of the energy storage device.

したがって、測定対象となる蓄電デバイスの特性を示す測定結果を測定値として取得することが可能となる。 Therefore, it is possible to obtain measurement results that indicate the characteristics of the energy storage device being measured as measured values.

図1は、一実施形態に係る測定システムを示すブロック図である。FIG. 1 is a block diagram illustrating a measurement system according to one embodiment. 図2は、一実施形態に係る測定システムの計測装置を示すブロック図である。FIG. 2 is a block diagram illustrating a measurement device of a measurement system according to an embodiment. 図3は、一実施形態に係る測定システムの取得処理装置を示すブロック図である。FIG. 3 is a block diagram illustrating an acquisition processing device of a measurement system according to an embodiment. 図4は、一実施形態に係る取得処理装置を示す機能ブロック図である。FIG. 4 is a functional block diagram illustrating an acquisition processing device according to an embodiment. 図5は、一実施形態に係る測定システムで蓄電デバイスのインピーダンスを測定する状態を示す等価回路である。FIG. 5 is an equivalent circuit showing a state in which the impedance of an electricity storage device is measured by a measurement system according to an embodiment. 図6は、一実施形態に係る測定システムで測定する蓄電デバイスの各周波数でのインピーダンスを複素平面に示したコールコールプロットを示す説明図である。FIG. 6 is an explanatory diagram showing a Cole-Cole plot in which the impedance at each frequency of an electricity storage device measured by a measurement system according to one embodiment is plotted on a complex plane. 図7は、一実施形態に係る測定システムで測定する蓄電デバイスのコールコールプロットのおける溶液抵抗成分の部分を示す説明図である。FIG. 7 is an explanatory diagram showing a part of the solution resistance component in a Cole-Cole plot of an electricity storage device measured by a measurement system according to one embodiment. 図8は、一実施形態に係る測定システムで測定する蓄電デバイスのコールコールプロットにおける電気二重層容量成分の部分を示す説明図である。FIG. 8 is an explanatory diagram showing a portion of the electric double layer capacitance component in a Cole-Cole plot of an electricity storage device measured by a measurement system according to one embodiment. 図9は、一実施形態に係る測定システムで測定する蓄電デバイスのコールコールプロットにおける反応抵抗成分の部分を示す説明図である。FIG. 9 is an explanatory diagram showing a reaction resistance component portion in a Cole-Cole plot of an electricity storage device measured by a measurement system according to one embodiment. 図10は、一実施形態に係る測定処理の動作を示すフローチャートである。FIG. 10 is a flowchart showing the operation of the measurement process according to an embodiment. 図11は、一実施形態に係る測定処理で予測値を予測する様子を示す説明図である。FIG. 11 is an explanatory diagram showing how a predicted value is predicted in a measurement process according to an embodiment. 図12は、一実施形態に係る周波数算出処理の動作を示すフローチャートである。FIG. 12 is a flowchart showing the operation of the frequency calculation process according to an embodiment.

以下、添付図面を参照しながら本発明の各実施形態について説明する。 Each embodiment of the present invention will be described below with reference to the attached drawings.

<実施形態>
以下、添付図面を参照しながら一実施形態について説明する。図1は、一実施形態に係る測定システム10を示すブロック図である。図2は、一実施形態に係る測定システム10の計測装置12を示すブロック図である。図3は、一実施形態に係る測定システム10の取得処理装置14を示すブロック図である。
<Embodiment>
An embodiment will be described below with reference to the accompanying drawings. Fig. 1 is a block diagram showing a measurement system 10 according to an embodiment. Fig. 2 is a block diagram showing a measurement device 12 of the measurement system 10 according to an embodiment. Fig. 3 is a block diagram showing an acquisition processing device 14 of the measurement system 10 according to an embodiment.

測定システム10は、図1に示すように、蓄電デバイス16のインピーダンスを測定するシステムである。 As shown in FIG. 1, the measurement system 10 is a system that measures the impedance of the power storage device 16.

ここで、蓄電デバイス16は、充放電可能に構成された蓄電池であり、電気二重層キャパシタなどのコンデンサ型の蓄電素子を含む。蓄電デバイス16は、複数の素電池が並列、直列又は直並列に接続された組電池でもよく、単電池であってもよい。蓄電デバイス16は、内部抵抗を有し、化学反応によって直流電圧を出力する。 Here, the power storage device 16 is a storage battery configured to be capable of being charged and discharged, and includes a capacitor-type power storage element such as an electric double layer capacitor. The power storage device 16 may be a battery pack in which multiple cells are connected in parallel, series, or series-parallel, or may be a single cell. The power storage device 16 has an internal resistance and outputs a DC voltage through a chemical reaction.

(ハードウエア構成)
測定システム10は、蓄電デバイス16のインピーダンスを測定するための複数の装置で構成されている。具体的に測定システム10は、蓄電デバイス16のインピーダンスを測定する計測装置12と、計測装置12での測定結果を利用して測定対象成分のインピーダンスを求める取得処理装置14とを備えている。
(Hardware configuration)
The measurement system 10 is composed of a plurality of devices for measuring the impedance of the electricity storage device 16. Specifically, the measurement system 10 includes a measuring device 12 that measures the impedance of the electricity storage device 16, and an acquisition processing device 14 that uses the measurement results from the measuring device 12 to determine the impedance of the measurement target component.

計測装置12は、蓄電デバイス16のインピーダンスを測定する。取得処理装置14は、計測装置12で測定したインピーダンスを処理して測定対象成分の測定値を確定する。 The measuring device 12 measures the impedance of the power storage device 16. The acquisition processing device 14 processes the impedance measured by the measuring device 12 to determine the measurement value of the component to be measured.

具体的に説明すると取得処理装置14は、計測装置12において、各周波数で測定したインピーダンスの実数成分と虚数成分とを複素平面にプロットしてコールコールプロットを形成する。これにより、取得処理装置14は、計測装置12での測定結果を利用して測定対象成分のインピーダンスを求める。 Specifically, the acquisition processing device 14 plots the real and imaginary components of the impedance measured at each frequency in the measurement device 12 on a complex plane to form a Cole-Cole plot. In this way, the acquisition processing device 14 uses the measurement results from the measurement device 12 to determine the impedance of the component to be measured.

測定対象成分としては、例えば、電解液抵抗を検出できる溶液抵抗成分、及び電気二重層容量を検出できる電気二重層容量成分が挙げられる。 Examples of components to be measured include a solution resistance component that can detect the electrolyte resistance, and an electric double layer capacitance component that can detect the electric double layer capacitance.

[計測装置]
図2に示す計測装置12は、図1に示す定電流源20と、電圧検出部22と、処理部24と、通信部26とを備える。定電流源20は、交流電流供給部と電流検出部とを備える。
[Measuring equipment]
The measuring device 12 shown in Fig. 2 includes the constant current source 20, the voltage detection unit 22, the processing unit 24, and the communication unit 26 shown in Fig. 1. The constant current source 20 includes an AC current supply unit and a current detection unit.

通信部26は、取得処理装置14と通信を行い、取得処理装置14から受信した測定条件などを処理部24へ引き渡す。また、通信部26は、処理部24から受けた測定結果などを取得処理装置14に引き渡す。 The communication unit 26 communicates with the acquisition processing device 14 and passes on the measurement conditions and other information received from the acquisition processing device 14 to the processing unit 24. The communication unit 26 also passes on the measurement results and other information received from the processing unit 24 to the acquisition processing device 14.

処理部24は、通信部26より受けた測定条件などに従って定電流源20の交流電流供給部を作動する。また、処理部24は、定電流源20の電流検出部で測定した電流値を受け取る。 The processing unit 24 operates the AC current supply unit of the constant current source 20 according to the measurement conditions received from the communication unit 26. The processing unit 24 also receives the current value measured by the current detection unit of the constant current source 20.

また、処理部24は、定電流源20の交流電流供給部からの出力と、電圧検出部22で検出した電圧の変化と、電流検出部で測定した電流値とに基づいて、インピーダンスを算出する。そして、処理部24は、インピーダンスの虚数成分及び実数成分を測定結果として通信部26に出力する。 The processing unit 24 also calculates the impedance based on the output from the AC current supply unit of the constant current source 20, the voltage change detected by the voltage detection unit 22, and the current value measured by the current detection unit. The processing unit 24 then outputs the imaginary and real components of the impedance to the communication unit 26 as the measurement results.

定電流源20の交流電流供給部は、処理部24からの測定条件に応じた大きさ及び周波数の交流電流を生成する。また、定電流源20の交流電流供給部は、生成した交流電流を蓄電デバイス16の正極端子16A及び負極端子16Bに印加する。定電流源20の電流検出部は、蓄電デバイス16に流れる電流を検出する。 The AC current supply unit of the constant current source 20 generates an AC current with a magnitude and frequency according to the measurement conditions from the processing unit 24. The AC current supply unit of the constant current source 20 also applies the generated AC current to the positive terminal 16A and the negative terminal 16B of the power storage device 16. The current detection unit of the constant current source 20 detects the current flowing through the power storage device 16.

電圧検出部22は、定電流源20から蓄電デバイス16に供給した交流電流に応じて蓄電デバイス16の正極端子16A及び負極端子16Bに生ずる電圧を検出する。また、電圧検出部22は、検出した電圧を処理部24へ出力する。 The voltage detection unit 22 detects the voltage generated at the positive terminal 16A and the negative terminal 16B of the power storage device 16 in response to the AC current supplied from the constant current source 20 to the power storage device 16. The voltage detection unit 22 also outputs the detected voltage to the processing unit 24.

[取得処理装置]
取得処理装置14は、図3に示すように、コンピュータを構成するプロセッサ30を中心に構成されている。プロセッサ30には、通信部32と、記憶部34と、入力部36と、表示部38と、報知部40と、時計部42と、外部通信部44とが接続されている。
[Acquisition and processing device]
3, the acquisition processing device 14 is configured with a processor 30 constituting a computer as a core component. A communication unit 32, a storage unit 34, an input unit 36, a display unit 38, a notification unit 40, a clock unit 42, and an external communication unit 44 are connected to the processor 30.

通信部32は、計測装置12の通信部26と通信を行う。通信は、有線通信又は無線通信で行われる。通信部32は、プロセッサ30から送られた測定条件などを計測装置12の通信部26に送信する。また、通信部32は、計測装置12の通信部26から受けた測定結果などを受信してプロセッサ30に引き渡す。 The communication unit 32 communicates with the communication unit 26 of the measuring device 12. The communication is performed by wired communication or wireless communication. The communication unit 32 transmits the measurement conditions and the like sent from the processor 30 to the communication unit 26 of the measuring device 12. The communication unit 32 also receives the measurement results and the like received from the communication unit 26 of the measuring device 12 and passes them over to the processor 30.

記憶部34は、プロセッサ30によってデータを読み出し可能に記憶する。記憶部34には、取得処理装置14の動作を制御する処理プログラムが格納される。記憶部34は、取得処理装置14の機能を実現する処理プログラムを格納する記憶媒体として機能する。 The memory unit 34 stores data so that it can be read by the processor 30. The memory unit 34 stores a processing program that controls the operation of the acquisition processing device 14. The memory unit 34 functions as a storage medium that stores a processing program that realizes the functions of the acquisition processing device 14.

記憶部34は、不揮発性メモリ(ROM:Read Only Memory)、及び揮発性メモリ(RAM:Random Access Memory)などにより構成される。 The memory unit 34 is composed of non-volatile memory (ROM: Read Only Memory) and volatile memory (RAM: Random Access Memory), etc.

また、記憶部34は、処理プログラムで使用するデータが読み出し可能に記憶される。 The memory unit 34 also stores data used by the processing program in a readable manner.

具体的に説明すると、記憶部34には、測定対象成分のインピーダンスを測定するための周波数帯域を示す情報が記憶されている。また、記憶部34は、測定対象成分の測定を開始する際に入力された周波数、及び測定対象成分の測定値を演算する過程で使用するデータなどを一時的に記憶する。さらに、記憶部34は、測定結果を記憶する。 Specifically, the memory unit 34 stores information indicating the frequency band for measuring the impedance of the component to be measured. The memory unit 34 also temporarily stores the frequency input when starting the measurement of the component to be measured, and data used in the process of calculating the measurement value of the component to be measured. The memory unit 34 also stores the measurement results.

入力部36は、利用者が入力したデータをプロセッサ30に送る。入力部36は、利用者の入力操作を受け付ける入力インターフェースとして機能する。入力部36は、一例として、複数の操作ボタン及び数字ボタン、又はタッチパネルで構成される。 The input unit 36 sends data input by the user to the processor 30. The input unit 36 functions as an input interface that accepts input operations by the user. As an example, the input unit 36 is composed of a plurality of operation buttons and numeric buttons, or a touch panel.

表示部38は、プロセッサ30からのデータに従って表示を行う。一例として、表示部38は、測定結果等を表示する。また、取得処理装置14が測定結果に基づいて、蓄電デバイス16をランク分けする機能を有する場合、表示部38は、ランク分けされたランクを表示する。表示するランクとしては、蓄電デバイス16の良不良を示す良否、又は蓄電デバイス16を用途ごとに分類する為の品質が挙げられる。 The display unit 38 performs display according to data from the processor 30. As an example, the display unit 38 displays measurement results, etc. Furthermore, if the acquisition processing device 14 has a function of ranking the power storage device 16 based on the measurement results, the display unit 38 displays the ranked rank. The displayed rank may be good or bad, indicating whether the power storage device 16 is good or bad, or quality, for classifying the power storage device 16 according to its use.

なお、蓄電デバイス16の出力特性の良否は、内部インピーダンスの状態に基づいて判断する。具体的には、内部インピーダンスが所定の基準よりも低いものを良品とする。 The quality of the output characteristics of the power storage device 16 is determined based on the state of the internal impedance. Specifically, a device with an internal impedance lower than a predetermined standard is considered to be a good product.

表示部38を構成する装置としては、発光ダイオード又はLCD(Liquid Crystal Display)等の表示パネルが挙げられる。本実施形態の表示部38は、一例として、液晶パネルで構成される。 The device constituting the display unit 38 may be a display panel such as a light-emitting diode or an LCD (Liquid Crystal Display). As an example, the display unit 38 in this embodiment is composed of a liquid crystal panel.

報知部40は、プロセッサ30からのデータに従って報知を行う。報知部40は、案内又は警告音等を利用者に音で報知する。音で報知する装置としては、圧電ブザー又はスピーカなどが挙げられる。本実施形態の報知部40は、一例として、スピーカで構成される。 The notification unit 40 issues notifications according to data from the processor 30. The notification unit 40 notifies the user by sound, such as guidance or warning sounds. Examples of devices that notify by sound include a piezoelectric buzzer or a speaker. In this embodiment, the notification unit 40 is, as an example, configured with a speaker.

時計部42は、現在の年月日及び時刻を示すととともに時間を測定する。時計部42は、現在の年月日及び時刻をプロセッサ30に出力する。時計部42が示す年月日及び時刻は、一例として、測定対象成分の測定値として確定した値に関連付けられて、記憶部34に記憶される。 The clock unit 42 indicates the current date and time and measures time. The clock unit 42 outputs the current date and time to the processor 30. The date and time indicated by the clock unit 42 are stored in the memory unit 34, for example, in association with a value determined as the measurement value of the component to be measured.

外部通信部44は、プロセッサ30と外部装置との間でデータの送受信を可能とする。外部通信部44は、データを送受信するためのインターフェースを構成する。外部通信部44は、USB(universal serial bus)、Bluetooth(登録商標)、無線LANなどを用いて通信を行うハードウエアで構成される。 The external communication unit 44 enables data transmission and reception between the processor 30 and an external device. The external communication unit 44 constitutes an interface for transmitting and receiving data. The external communication unit 44 is composed of hardware that performs communication using a USB (universal serial bus), Bluetooth (registered trademark), wireless LAN, etc.

前述した処理プログラムが外部装置から供給される場合、外部通信部44は、外部装置から処理プログラムを受信してプロセッサ30に送る。プロセッサ30は、受信した処理プログラムを記憶部34に格納する。 When the above-mentioned processing program is supplied from an external device, the external communication unit 44 receives the processing program from the external device and sends it to the processor 30. The processor 30 stores the received processing program in the memory unit 34.

プロセッサ30は、一例として、中央演算処理装置(CPU:Central Processing Unit)で構成される。プロセッサ30は、記憶部34に格納された処理プログラムを読み出すとともに、読み出した処理プログラムに従って動作する。これにより、プロセッサ30は、取得処理装置14の各部を制御して処理方法を実施する。 The processor 30 is, for example, configured with a central processing unit (CPU). The processor 30 reads a processing program stored in the memory unit 34 and operates according to the read processing program. In this way, the processor 30 controls each part of the acquisition processing device 14 to carry out the processing method.

また、プロセッサ30は、測定結果等を表示部38で表示したり、報知部40から報知したり、外部通信部44を介して外部装置へ送信したりする。 The processor 30 also displays the measurement results, etc. on the display unit 38, notifies them from the notification unit 40, and transmits them to an external device via the external communication unit 44.

(機能ブロック)
図4は、一実施形態に係る取得処理装置14の機能構成を示す機能ブロック図である。
(Function block)
FIG. 4 is a functional block diagram showing the functional configuration of the acquisition processing device 14 according to an embodiment.

取得処理装置14は、図4に示すように、取得部50と、決定部51と、算出部52と、確定部54とを備える。取得処理装置14における各部50、51、52、54の機能は、プロセッサ30が記憶部34から処理プログラムとして読み出したソフトウエアプログラムを実行することで実現される。 As shown in FIG. 4, the acquisition processing device 14 includes an acquisition unit 50, a determination unit 51, a calculation unit 52, and a confirmation unit 54. The functions of the units 50, 51, 52, and 54 in the acquisition processing device 14 are realized by the processor 30 executing a software program read from the storage unit 34 as a processing program.

[取得部]
取得手段として機能する取得部50は、測定対象成分のインピーダンスを測定するための所定の周波数帯域において、三点以上の周波数で測定されたインピーダンスを計測装置12から取得する。測定対象成分は、一例として蓄電デバイス16の溶液抵抗成分を含み、測定対象成分についての詳細は後述する。
[Acquisition Department]
The acquiring unit 50, which functions as an acquiring means, acquires impedance measured at three or more frequencies in a predetermined frequency band for measuring the impedance of the measurement target component from the measuring device 12. The measurement target component includes, as an example, a solution resistance component of the electricity storage device 16, and the measurement target component will be described in detail later.

また、取得部50は、算出部52で予測した対象成分の予測値が、計測装置12より取得した三点以上のインピーダンスの実数成分又は虚数成分の極値よりも測定対象成分のインピーダンスに近づくように突出するか否かを判断する。そして、予測値が極値よりも測定対象成分のインピーダンスに近づくように突出する場合には、周波数帯域のうち予測値に対応する予測周波数で測定された追加インピーダンスを取得する。 The acquisition unit 50 also determines whether the predicted value of the target component predicted by the calculation unit 52 is closer to the impedance of the measurement target component than the extreme values of the real or imaginary components of the impedance at three or more points acquired by the measurement device 12. If the predicted value is closer to the impedance of the measurement target component than the extreme values, an additional impedance measured at a predicted frequency in the frequency band that corresponds to the predicted value is acquired.

言い換えると、取得部50は、計測装置12より取得した値が、算出部52で予測した測定対象成分の予測値よりも測定対象成分のインピーダンスから遠いか否かを判断する。そして、取得部50は、取得した値が測定対象成分のインピーダンスよりも遠い場合には、上記周波数帯域のうち予測値に対応する周波数である予測周波数で測定した追加インピーダンスを測定結果として新たに取得する。 In other words, the acquisition unit 50 determines whether the value acquired from the measuring device 12 is farther from the impedance of the measured component than the predicted value of the measured component predicted by the calculation unit 52. If the acquired value is farther from the impedance of the measured component, the acquisition unit 50 newly acquires, as the measurement result, an additional impedance measured at a predicted frequency that is a frequency in the above frequency band that corresponds to the predicted value.

予測値に対応する予測周波数は、予測値の実数成分から求められる実数成分周波数と予測値の虚数成分から求められる虚数成分周波数とに基づいて算出される。 The predicted frequency corresponding to the predicted value is calculated based on the real component frequency obtained from the real component of the predicted value and the imaginary component frequency obtained from the imaginary component of the predicted value.

この予測周波数は、一例として、取得部50で取得した実数成分の最大値と最小値との差分及び虚数成分の最大値と最小値との差分に基づいて実数成分周波数及び虚数成分周波数に重みづけすることにより得られた値である加重平均値とする。 As an example, this predicted frequency is a weighted average value obtained by weighting the real component frequency and the imaginary component frequency based on the difference between the maximum and minimum values of the real component and the difference between the maximum and minimum values of the imaginary component acquired by the acquisition unit 50.

ここで、計測装置12より取得した値が予測値よりも測定対象成分のインピーダンスから遠いか否かは、測定対象成分に基づいて把握することできる。 Here, whether the value obtained from the measuring device 12 is farther from the impedance of the component being measured than the predicted value can be determined based on the component being measured.

具体例を挙げて説明すると、測定対象成分が電解液抵抗を含む溶液抵抗成分の場合、電解液抵抗は、測定したインピーダンスの実数成分のうち最も小さい値と相関がある。そして、測定したインピーダンスにおいて実数成分の最も小さい値が電解液抵抗を示すことが知られている。 To give a specific example, when the component to be measured is a solution resistance component that includes electrolyte resistance, the electrolyte resistance correlates with the smallest value among the real components of the measured impedance. It is known that the smallest value of the real component in the measured impedance indicates the electrolyte resistance.

このため、測定対象成分が電解液抵抗を含む溶液抵抗成分の場合において、測定したインピーダンスの実数成分が予測値のインピーダンスの実数成分よりも大きいときは、測定した値である極値が予測値よりも測定対象成分のインピーダンスから遠いと判断することができる。また、測定したインピーダンスの実数成分が予測値のインピーダンスの実数成分よりも小さいときは、測定した値である極値が予測値よりも測定対象成分のインピーダンスに近いと判断することができる。 Therefore, when the component to be measured is a solution resistance component including electrolyte resistance, if the real component of the measured impedance is larger than the real component of the predicted impedance, it can be determined that the extreme value, which is the measured value, is farther from the impedance of the component to be measured than the predicted value. Also, if the real component of the measured impedance is smaller than the real component of the predicted impedance, it can be determined that the extreme value, which is the measured value, is closer to the impedance of the component to be measured than the predicted value.

[決定部]
決定手段として機能する決定部51は、取得部50で取得した三点以上のインピーダンスの実数成分又は虚数成分の極値を求める。
[Decision section]
The determination unit 51 functioning as a determination means determines extreme values of real or imaginary components of the impedance at three or more points acquired by the acquisition unit 50 .

ここで、極値とは、取得した三点以上のインピーダンスの実数成分又は虚数成分の値のうち、測定対象成分のインピーダンスに最も近い値を示す。なお、測定対象成分のインピーダンスに近いか否かは、前述した方法により判断することができる。 Here, the extreme value refers to the value of the real or imaginary components of the impedance obtained at three or more points that is closest to the impedance of the component to be measured. Whether or not a value is close to the impedance of the component to be measured can be determined by the method described above.

また、極値は、三点以上のインピーダンスの実数成分又は虚数成分を複素平面にプロットして曲線を形成した場合、この曲線において、最も曲線の突出方向に位置した値と言い換えることができる。 In addition, when the real or imaginary components of impedance at three or more points are plotted on a complex plane to form a curve, the extreme value can be said to be the value located in the direction in which the curve protrudes most.

[算出部]
算出手段として機能する算出部52は、取得部50が取得した三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて測定対象成分の予測値を算出する。一例として、測定対象成分が電解液抵抗を含む溶液抵抗成分の場合、測定対象成分の予測値は、実数成分の値とする。
[Calculation section]
The calculation unit 52, which functions as a calculation means, calculates a predicted value of the measurement target component based on an approximation formula expressing the relationship between the real components and imaginary components of the impedance at three or more points acquired by the acquisition unit 50. As an example, when the measurement target component is a solution resistance component including an electrolyte resistance, the predicted value of the measurement target component is the value of the real component.

また、算出部52は、三点以上のインピーダンスの実数成分及び虚数成分の関係を二次関数で表した近似式の頂点を測定対象成分の予測値として算出する。 In addition, the calculation unit 52 calculates the apex of an approximation equation that expresses the relationship between the real and imaginary components of the impedance at three or more points as a quadratic function as a predicted value of the component to be measured.

具体的に算出部52は、取得部50で取得したインピーダンスの実数成分及び虚数成分の各値から最小二乗法によって二次曲線を算出し、二次曲線が突出方向に最も突出した曲点に基づいて予測値を定める。 Specifically, the calculation unit 52 calculates a quadratic curve from the real and imaginary components of the impedance acquired by the acquisition unit 50 using the least squares method, and determines a predicted value based on the inflection point where the quadratic curve is most prominent in the protruding direction.

算出部52は、取得部50で追加インピーダンスを取得した場合、追加インピーダンスの実数成分及び虚数成分と、追加インピーダンスよりも前に測定した三点以上のインピーダンスの実数成分及び虚数成分とに基づいて近似式を求める。そして、取得部50は、求めた近似式に基づいて、測定対象成分の予測値を新たに算出する。 When the acquisition unit 50 acquires an additional impedance, the calculation unit 52 obtains an approximation formula based on the real and imaginary components of the additional impedance and the real and imaginary components of the impedances at three or more points measured before the additional impedance. The acquisition unit 50 then calculates a new predicted value of the measurement target component based on the obtained approximation formula.

[確定部]
確定手段である確定部54は、決定部51で求めた極値が、前述の予測値と同じ又は予測値よりも測定対象成分のインピーダンスに近づくように突出する場合、極値を測定対象成分に関する測定値として確定する。
[Determined part]
The determination unit 54, which is a determination means, determines the extreme value as the measurement value for the measured component when the extreme value obtained by the decision unit 51 stands out so as to be the same as the predicted value or closer to the impedance of the measured component than the predicted value.

言い換えると確定部54は、三点以上の実数成分又は虚数成分の各値のうち測定対象成分に最も近い値が、予測値と同じ又は予測値よりも測定対象成分に近い場合に、その値を測定対象成分に関する値、すなわち最終結果として確定する。 In other words, if the value closest to the measured component among the three or more real or imaginary component values is the same as the predicted value or is closer to the measured component than the predicted value, the determination unit 54 determines that value as the value related to the measured component, i.e., the final result.

これにより、測定対象成分のインピーダンスは、蓄電デバイスを測定した値に基づいて定められる。 As a result, the impedance of the component to be measured is determined based on the value measured on the energy storage device.

取得部50で取得した追加インピーダンスを用いて測定対象成分の予測値を新たに算出した場合、確定部54は、次のようにして測定対象成分に関する値を確定する。すなわち、確定部54は、測定対象成分に最も近い測定値が、算出部52が新たに算出した予測値と同じ又は予測値よりも測定対象成分に近い場合に、測定対象成分に関する値として測定値を確定する。 When a new predicted value of the component to be measured is calculated using the additional impedance acquired by the acquisition unit 50, the determination unit 54 determines the value related to the component to be measured as follows. That is, when the measured value closest to the component to be measured is the same as the predicted value newly calculated by the calculation unit 52 or is closer to the component to be measured than the predicted value, the determination unit 54 determines the measured value as the value related to the component to be measured.

ここで、測定対象成分に最も近い測定値は、予測周波数で測定された追加インピーダンスを含めた総ての値において、測定対象成分のインピーダンスに最も近い値とする。 Here, the measurement value closest to the component to be measured is the value that is closest to the impedance of the component to be measured among all values, including the additional impedance measured at the predicted frequency.

一例として、測定対象成分が電解液抵抗を含む溶液抵抗成分の場合、測定対象成分の予測値は、実数成分の値とする。この場合、確定部54は、三点以上の実数成分の各値のうち最も小さい値が、予測値と同じ又は予測値よりも小さい場合に、測定対象成分に関する値として測定値を確定する。 As an example, when the component to be measured is a solution resistance component including electrolyte resistance, the predicted value of the component to be measured is the value of the real component. In this case, the determination unit 54 determines the measured value as the value related to the component to be measured when the smallest value among the values of the three or more real components is the same as or smaller than the predicted value.

これにより、電解液抵抗を含む溶液抵抗成分を示す値が測定対象成分に関する値として確定される。 This allows the value indicating the solution resistance component, including the electrolyte resistance, to be determined as the value relating to the component being measured.

ここで、極値が予測値よりも測定対象成分のインピーダンスから遠い場合、インピーダンスの測定、予測値の算出、及び測定した値を加えた極値と予測値との比較を繰り返すが、この繰返回数は、予め指定することも可能である。 If the extreme value is farther from the impedance of the component being measured than the predicted value, the process of measuring the impedance, calculating the predicted value, and comparing the extreme value to which the measured value is added with the predicted value is repeated, but the number of repetitions can also be specified in advance.

(等価回路)
図5は、一実施形態に係る測定システム10で蓄電デバイス16のインピーダンスを測定する状態を示す等価回路60である。
(Equivalent circuit)
FIG. 5 is an equivalent circuit 60 illustrating a state in which the impedance of the power storage device 16 is measured by the measurement system 10 according to one embodiment.

計測装置12のケーブルを蓄電デバイス16の各端子16A、16Bに接続した状態での蓄電デバイス16のインピーダンスは、ケーブル抵抗成分62と、溶液抵抗成分64と、反応抵抗成分66とが直列接続された等価回路60によって表すことができる。 The impedance of the energy storage device 16 when the cable of the measuring device 12 is connected to each terminal 16A, 16B of the energy storage device 16 can be represented by an equivalent circuit 60 in which a cable resistance component 62, a solution resistance component 64, and a reaction resistance component 66 are connected in series.

ケーブル抵抗成分62は、互いに並列接続されたインダクタンスLiと抵抗Riとで示される。インダクタンスLiは、ケーブルが有するインダクタンスを示し、抵抗Riは、ケーブルが有する抵抗を示す。 The cable resistance component 62 is represented by an inductance Li and a resistance Ri connected in parallel to each other. The inductance Li represents the inductance of the cable, and the resistance Ri represents the resistance of the cable.

溶液抵抗成分64は、蓄電デバイス16の電極間の電解液の溶液抵抗を示す。溶液抵抗成分64は、蓄電デバイス16の電解液抵抗及び電極のタブの溶接抵抗からなる電解液抵抗Rsolを有する。 The solution resistance component 64 represents the solution resistance of the electrolyte between the electrodes of the electricity storage device 16. The solution resistance component 64 has an electrolyte resistance Rsol consisting of the electrolyte resistance of the electricity storage device 16 and the welding resistance of the electrode tabs.

反応抵抗成分66は、正極の電極で生ずる正極反応抵抗68と、負極の電極で生ずる負極反応抵抗70とで構成される。正極反応抵抗68と負極反応抵抗70とは、等価回路60において直列接続される。 The reaction resistance component 66 is composed of a positive electrode reaction resistance 68 generated at the positive electrode and a negative electrode reaction resistance 70 generated at the negative electrode. The positive electrode reaction resistance 68 and the negative electrode reaction resistance 70 are connected in series in the equivalent circuit 60.

正極反応抵抗68は、並列接続された電気二重層容量C1と反応抵抗R1とを有する。電気二重層容量C1は、正極の電極と電荷液との間に生ずる電気容量を示し、反応抵抗R1は、正極の電極と電解液との間に生ずる電荷移動抵抗を示す。 The positive electrode reaction resistor 68 has an electric double layer capacitance C1 and a reaction resistance R1 connected in parallel. The electric double layer capacitance C1 indicates the electric capacitance occurring between the positive electrode and the electrolytic solution, and the reaction resistance R1 indicates the charge transfer resistance occurring between the positive electrode and the electrolytic solution.

負極反応抵抗70は、互いに並列接続された電気二重層容量C2と反応抵抗R2とを有する。電気二重層容量C2は、負極の電極と電荷液との間に生ずる電気容量を示し、反応抵抗R2は、負極の電極と電解液との間に生ずる電荷移動抵抗を示す。 The negative electrode reaction resistor 70 has an electric double layer capacitance C2 and a reaction resistance R2 connected in parallel to each other. The electric double layer capacitance C2 indicates the electric capacitance occurring between the negative electrode and the electrolytic solution, and the reaction resistance R2 indicates the charge transfer resistance occurring between the negative electrode and the electrolytic solution.

図6は、一実施形態に係る測定システム10で測定する蓄電デバイス16の各周波数でのインピーダンスを複素平面に示したコールコールプロット80を示す説明図である。図7は、一実施形態に係る測定システム10で測定された蓄電デバイス16のコールコールプロット80のおける溶液抵抗成分64の部分を示す説明図である。 Figure 6 is an explanatory diagram showing a Cole-Cole plot 80 in which the impedance at each frequency of the electricity storage device 16 measured by the measurement system 10 according to one embodiment is plotted on a complex plane. Figure 7 is an explanatory diagram showing the solution resistance component 64 in the Cole-Cole plot 80 of the electricity storage device 16 measured by the measurement system 10 according to one embodiment.

図8は、一実施形態に係る測定システム10で測定された蓄電デバイス16のコールコールプロット80における電気二重層容量成分72の部分を示す説明図である。図9は、一実施形態に係る測定システム10で測定する蓄電デバイス16のコールコールプロット80における反応抵抗成分74の部分を示す説明図である。 Figure 8 is an explanatory diagram showing the electric double layer capacitance component 72 in the Cole-Cole plot 80 of the power storage device 16 measured by the measurement system 10 according to one embodiment. Figure 9 is an explanatory diagram showing the reaction resistance component 74 in the Cole-Cole plot 80 of the power storage device 16 measured by the measurement system 10 according to one embodiment.

図6に示すように、コールコールプロット80(ナイキストプロットともいう)は、各周波数における蓄電デバイス16のインピーダンスを複素平面にプロットしたものである。複素平面の横軸(R)は、インピーダンスの実数成分を示し、縦軸(X)は、インピーダンスの虚数成分を示す。 As shown in FIG. 6, a Cole-Cole plot 80 (also called a Nyquist plot) is a plot of the impedance of the power storage device 16 at each frequency on a complex plane. The horizontal axis (R) of the complex plane indicates the real component of the impedance, and the vertical axis (X) indicates the imaginary component of the impedance.

このコールコールプロット80は、図6及び図7に示すように、電解液抵抗Rsolがインピーダンスの主要因となる領域を有する。電解液抵抗Rsolがインピーダンスの主要因となる領域は、第一周波数帯域に現われ、本実施形態の蓄電デバイス16の場合、第一周波数帯域は、一例として、数100Hz以上数10kHz未満の帯域である。 As shown in Figures 6 and 7, this Cole-Cole plot 80 has a region where electrolyte resistance Rsol is the main factor in impedance. The region where electrolyte resistance Rsol is the main factor in impedance appears in a first frequency band, and in the case of the electricity storage device 16 of this embodiment, the first frequency band is, for example, a band between several hundred Hz and several tens of kHz.

また、コールコールプロット80は、図6及び図8に示すように、電気二重層容量C1、C2の影響が表れるインピーダンスの領域を有する。電気二重層容量C1、C2の影響が表れるインピーダンスの領域は、第二周波数帯域に現われ、本実施形態の蓄電デバイス16の場合、第二周波数帯域は、一例として、数10Hz以上数100Hz未満の帯域である。 The Cole-Cole plot 80 also has an impedance region in which the effects of the electric double layer capacitances C1 and C2 appear, as shown in Figures 6 and 8. The impedance region in which the effects of the electric double layer capacitances C1 and C2 appear appears in the second frequency band, and in the case of the power storage device 16 of this embodiment, the second frequency band is, for example, a band of several tens of Hz or more and less than several hundreds of Hz.

さらに、コールコールプロット80は、図6及び図9に示すように、反応抵抗R1、R2及び溶液の拡散領域のインピーダンスの領域を有する。反応抵抗R1、R2及び拡散領域のインピーダンスの領域は、第三周波数帯域に現われ、本実施形態の蓄電デバイス16の場合、第三周波数帯域は、一例として、0.01Hz以上数10Hz未満の帯域である。 Furthermore, as shown in Figures 6 and 9, the Cole-Cole plot 80 has regions of reaction resistances R1, R2 and the impedance of the diffusion region of the solution. The reaction resistances R1, R2 and the impedance region of the diffusion region appear in a third frequency band, and in the case of the power storage device 16 of this embodiment, the third frequency band is, for example, a band of 0.01 Hz or more and less than several tens of Hz.

本実施形態において、測定対象成分とは、電解液抵抗Rsolを示す溶液抵抗成分64と、電気二重層容量C1、C2を示す電気二重層容量成分72のことである。 In this embodiment, the components to be measured are the solution resistance component 64 indicating the electrolyte resistance Rsol and the electric double layer capacitance component 72 indicating the electric double layer capacitances C1 and C2.

ここで、反応抵抗R1、R2及び拡散領域の良否は、一例として、高周波エリアで検出された電解液抵抗Rsolのインピーダンスと低周波エリアの反応抵抗R1、R2のインピーダンスとに基づいて評価することができる。具体的に説明すると、コールコールプロットにおいて、高周波エリアで検出された電解液抵抗Rsolのインピーダンスの位置と低周波エリアの反応抵抗R1、R2のインピーダンスの位置とを結んだ直線の長さ等によって評価することができる。 Here, the quality of the reaction resistances R1, R2 and the diffusion region can be evaluated based on, for example, the impedance of the electrolyte resistance Rsol detected in the high frequency area and the impedance of the reaction resistances R1, R2 in the low frequency area. To be more specific, in the Cole-Cole plot, the quality can be evaluated based on the length of the straight line connecting the position of the impedance of the electrolyte resistance Rsol detected in the high frequency area and the position of the impedance of the reaction resistances R1, R2 in the low frequency area.

このため、反応抵抗R1、R2及び拡散領域を示す反応抵抗成分74も、測定対象成分のインピーダンスとすることができる。 Therefore, the reaction resistances R1, R2 and the reaction resistance component 74 indicating the diffusion region can also be considered as the impedance of the component to be measured.

第一周波数帯域は、電解液抵抗Rsolを測定するための周波数帯域を示し、第二周波数帯域は、電気二重層容量C1、C2を測定するための周波数帯域を示す。また、第三周波数帯域は、反応抵抗R1、R2及び拡散領域のインピーダンスを測定するための周波数帯域を示す。 The first frequency band indicates a frequency band for measuring the electrolyte resistance Rsol, the second frequency band indicates a frequency band for measuring the electric double layer capacitances C1 and C2, and the third frequency band indicates a frequency band for measuring the reaction resistances R1 and R2 and the impedance of the diffusion region.

(動作説明)
次に、取得処理装置14の動作を、図10から図12を用いるとともに、取得処理装置14のプロセッサ30が実行する処理手順に従って説明する。なお、測定システム10は、取得処理装置14の動作に伴って機能する。
(Operation description)
Next, the operation of the acquisition processing device 14 will be described with reference to Figures 10 to 12 and in accordance with the processing procedure executed by the processor 30 of the acquisition processing device 14. Note that the measurement system 10 functions in conjunction with the operation of the acquisition processing device 14.

図10は、取得処理装置14の動作の一例を示すフローチャートである。図10には、蓄電デバイス16のインピーダンスの測定対象成分を測定するための測定処理が示されている。 Figure 10 is a flowchart showing an example of the operation of the acquisition processing device 14. Figure 10 shows a measurement process for measuring the measurement target component of the impedance of the power storage device 16.

ここで、この測定処理を実行する前には、測定処理に必要となるパラメータを設定するための設定処理が行われる。この設定処理では、蓄電デバイス16において測定対象とする測定対象成分が利用者によって選択されている。また、設定処理では、選択された測定対象成分に対応する周波数帯域において三点の周波数が利用者によって入力されているものとする。 Before executing this measurement process, a setting process is performed to set the parameters required for the measurement process. In this setting process, the measurement target component to be measured in the power storage device 16 is selected by the user. In addition, in the setting process, it is assumed that the user inputs three frequencies in the frequency band corresponding to the selected measurement target component.

一例を挙げて説明すると、本実施形態に係る蓄電デバイス16において、測定対象成分として電解液抵抗Rsolが選択されている。また、測定する周波数として、第一周波数帯域(数100Hz以上数10kHz未満の帯域)において、三点の周波数が入力されている。ここで、測定対象成分は電解液抵抗Rsolであるため、測定対象成分は、インピーダンスの実数成分の値で示される。 To explain by way of example, in the electricity storage device 16 according to this embodiment, the electrolyte resistance Rsol is selected as the component to be measured. In addition, three frequencies are input in the first frequency band (a band of several hundred Hz or more and less than several tens of kHz) as the frequencies to be measured. Here, since the component to be measured is the electrolyte resistance Rsol, the component to be measured is represented by the value of the real component of the impedance.

なお、本実施形態では、第一周波数帯域の周波数での測定結果を用いて測定対象成分である電解液抵抗Rsolを測定する場合について説明するが、これに限定されるものではない。例えば、他の周波数帯域の各周波数での測定結果を用いることで、他の測定対象成分を測定することができる。 In this embodiment, the electrolyte resistance Rsol, which is the component to be measured, is measured using the measurement results at a frequency in the first frequency band, but this is not limited to the above. For example, other components to be measured can be measured by using the measurement results at each frequency in other frequency bands.

取得処理装置14のプロセッサ30は、記憶部34に記憶された処理プログラムの測定処理を実行する。すると、プロセッサ30は、測定対象成分を測定するための周波数帯域において、予め入力された三点の周波数で測定されたインピーダンスを取得し、取得した各インピーダンスの実数成分及び虚数成分を得る(ステップS1)。 The processor 30 of the acquisition processing device 14 executes the measurement process of the processing program stored in the memory unit 34. The processor 30 then acquires impedances measured at three pre-input frequencies in the frequency band for measuring the measurement target components, and obtains the real and imaginary components of each acquired impedance (step S1).

具体的に説明すると、取得処理装置14のプロセッサ30は、通信部32を介して計測装置12と通信する。この通信により、プロセッサ30は、計測装置12に対して、入力された三点の周波数で蓄電デバイス16のインピーダンスの測定を実施させるとともに、測定した各インピーダンスを計測装置12から取得する。 Specifically, the processor 30 of the acquisition processing device 14 communicates with the measurement device 12 via the communication unit 32. Through this communication, the processor 30 causes the measurement device 12 to measure the impedance of the power storage device 16 at the three input frequencies, and acquires each measured impedance from the measurement device 12.

そして、プロセッサ30は、計測装置12から各インピーダンスの実数成分及び虚数成分を取得して記憶部34に記憶する。 Then, the processor 30 acquires the real and imaginary components of each impedance from the measuring device 12 and stores them in the memory unit 34.

ここで、実施形態では、測定した各インピーダンスの実数成分及び虚数成分を計測装置12から取得する場合について説明するが、これに限定されるものではない。例えば、計測装置12から出力されたインピーダンスを、取得処理装置14のプロセッサ30で実数成分と虚数成分とに分解し、その実数成分及び虚数成分を記憶部34に記憶してもよい。 Here, in the embodiment, a case will be described in which the real and imaginary components of each measured impedance are acquired from the measuring device 12, but this is not limited to this. For example, the impedance output from the measuring device 12 may be decomposed into real and imaginary components by the processor 30 of the acquisition processing device 14, and the real and imaginary components may be stored in the memory unit 34.

図11は、一実施形態に係る測定処理で予測値102を予測する様子を示す説明図である。図11には、三点の周波数で測定したインピーダンスとして、第一実測インピーダンス90、第二実測インピーダンス92、及び第三実測インピーダンス94が、複素平面に示されている。 Figure 11 is an explanatory diagram showing how a predicted value 102 is predicted in a measurement process according to one embodiment. In Figure 11, a first measured impedance 90, a second measured impedance 92, and a third measured impedance 94 are shown on a complex plane as impedances measured at three frequencies.

プロセッサ30は、図10に示したように、各実測インピーダンス90、92、94の実数成分及び虚数成分の各値のうち測定対象成分に最も近い値を極値とする(ステップS2)。 As shown in FIG. 10, the processor 30 determines the extreme value among the real and imaginary components of each measured impedance 90, 92, 94 that is closest to the component to be measured (step S2).

本実施形態において、測定対象成分のインピーダンスは電解液抵抗Rsolである。このため、プロセッサ30は、各実測インピーダンス90、92、94のうち実数成分が最も小さい第二実測インピーダンス92の実数成分を極値として定める。 In this embodiment, the impedance of the component to be measured is the electrolyte resistance Rsol. Therefore, the processor 30 determines the real component of the second measured impedance 92, which has the smallest real component among the measured impedances 90, 92, and 94, as the extreme value.

そして、プロセッサ30は、記憶部34に予め設定された「カウンタ」を「0」とし(ステップS2-1)、「カウンタ」がn未満であるか否かを判断する(ステップS2-2)。なお、nは、予め定めらえた値であり、例えば「3」等の数字で構成される。また、このnは、一例として、入力部からの入力によって予め定めてもよい。 Then, the processor 30 sets the "counter" preset in the memory unit 34 to "0" (step S2-1), and determines whether the "counter" is less than n (step S2-2). Note that n is a predetermined value, and is composed of a number such as "3". Also, as one example, this n may be predetermined by input from the input unit.

現時点において、「カウンタ」は、「0」なので、ステップS2-2の判断では、「カウンタ」はn未満であると判断し、ステップS3へ移行する。 At this point, the "counter" is "0", so in step S2-2 it is determined that the "counter" is less than n, and the process moves to step S3.

ステップS3において、プロセッサ30は、取得した各実測インピーダンス90、92、94から二次曲線100を算出する(ステップS3)。これにより、各実測インピーダンス90、92、94の実数成分及び虚数成分の関係が二次関数で表される。 In step S3, the processor 30 calculates a quadratic curve 100 from the acquired measured impedances 90, 92, and 94 (step S3). As a result, the relationship between the real and imaginary components of the measured impedances 90, 92, and 94 is expressed as a quadratic function.

具体的に説明すると、プロセッサ30は、図11に示したように、各実測インピーダンス90、92、94から一例として最小二乗法を用いて二次関数を求め、この二次関数が示す二次曲線100を取得する。なお、最小二乗法以外の近似法を用いて関数を求めてもよい。 Specifically, as shown in FIG. 11, the processor 30 calculates a quadratic function from each measured impedance 90, 92, 94 using, for example, the least squares method, and obtains the quadratic curve 100 represented by this quadratic function. Note that the function may be calculated using an approximation method other than the least squares method.

そして、プロセッサ30は、取得した二次曲線100から測定対象成分の予測値102を算出する(ステップS4)。この二次曲線100は、二次関数を表す近似式を示し、二次関数の頂点が予測値102として算出される。二次関数の頂点は、二次曲線100が突出方向に最も突出した曲点を示し、この曲点に基づいて予測値102が定められる。 Then, the processor 30 calculates the predicted value 102 of the measurement target component from the acquired quadratic curve 100 (step S4). This quadratic curve 100 indicates an approximation formula representing a quadratic function, and the apex of the quadratic function is calculated as the predicted value 102. The apex of the quadratic function indicates the inflection point where the quadratic curve 100 is most prominent in the protruding direction, and the predicted value 102 is determined based on this inflection point.

具体的に説明すると、二次曲線100は、Y=aX+bX+cの式で示される。この式を変形すると、Y=F(x)=a(X+b/2a)+(-b+4ac)/4aとなる。 More specifically, the quadratic curve 100 is expressed by the formula Y=aX 2 +bX+c. This formula can be rearranged to give Y=F(x)=a(X+b/2a) 2 +(-b 2 +4ac)/4a.

ここで、本実施形態では、測定対象成分として電解液抵抗Rsolが選択されている。このため、この二次曲線100は、便宜上、横軸をリアクタンス、縦軸を抵抗とするとともに、縦軸方向に凸となる曲線を示す式によって表す。 In this embodiment, the electrolyte resistance Rsol is selected as the component to be measured. For this reason, for convenience, the quadratic curve 100 is expressed by an equation that shows a curve that is convex in the vertical axis direction, with the horizontal axis representing reactance and the vertical axis representing resistance.

この変形した式から頂点に対応する予測値102の虚数成分の値Xtは、Xt=(-b/2a)となる。実数成分の値Rtは、Rt={(-b+4ac)/4a}となる。実数成分のうち最も小さい値は{(-b+4ac)/4a}となる。最も小さい実数成分は、頂点を示すので、この値を予測値102とすることができる。 From this transformed equation, the value Xt of the imaginary component of the predicted value 102 corresponding to the apex is Xt = (-b/2a). The value Rt of the real component is Rt = {(-b 2 + 4ac)/4a}. The smallest value of the real components is {(-b 2 + 4ac)/4a}. Since the smallest real component indicates the apex, this value can be used as the predicted value 102.

このように、虚数成分を変数とする二次関数F(x)の頂点に基づいて予測値102が定められる。 In this way, the predicted value 102 is determined based on the peak of the quadratic function F(x) with the imaginary component as a variable.

なお、本実施形態では、近似式として二次関数を例に挙げて説明するが、これに限定されるものではない。例えば、四点の周波数でインピーダンスを測定した場合、一方の隣接する二点を結ぶ直線と、他方の隣接する二点を結ぶ直線とを示す二つの近似式を想定する。この場合、二本の直線が交差する点を予測値102とすることができる。 In this embodiment, a quadratic function is used as an example of an approximation equation, but the present invention is not limited to this. For example, when impedance is measured at four frequencies, two approximation equations are assumed, one of which shows a straight line connecting two adjacent points, and the other of which shows a straight line connecting two adjacent points. In this case, the point where the two straight lines intersect can be set as the predicted value 102.

そして、プロセッサ30は、図10に示したように、各実測インピーダンス90、92、94の最小値(極値:以下同じ)が、予測値102よりも大きいか否かを判断する(ステップS5)。ステップS5において、各実測インピーダンス90、92、94の最小値が予測値102と同じ又は予測値102よりも小さいと判断した場合、各実測インピーダンス90、92、94の最小値を、測定対象成分に関する測定値として確定する(ステップS6)。そして、測定処理を終了する。 Then, as shown in FIG. 10, the processor 30 judges whether the minimum value (extreme value: the same below) of each measured impedance 90, 92, 94 is greater than the predicted value 102 (step S5). If it is judged in step S5 that the minimum value of each measured impedance 90, 92, 94 is equal to or smaller than the predicted value 102, the processor 30 determines the minimum value of each measured impedance 90, 92, 94 as the measured value for the component to be measured (step S6). Then, the measurement process ends.

これにより、ステップS1で取得した各実測インピーダンス90、92、94の実数成分又は虚数成分の各値のうち測定対象成分のインピーダンスに最も近い値が、予測値102と同じ又は予測値102よりも測定対象成分に近い場合に、測定対象成分に関する値として測定値が確定される。 As a result, if the value closest to the impedance of the component to be measured among the real or imaginary components of the measured impedances 90, 92, and 94 acquired in step S1 is the same as the predicted value 102 or is closer to the component to be measured than the predicted value 102, the measured value is determined as the value related to the component to be measured.

具体的に説明すると、測定対象成分は、溶液抵抗成分64を含み、測定対象成分の予測値102は、実数成分の値である。このため、取得した各実測インピーダンス90、92、94の実数成分の各値のうち最も小さい値が予測値102と同じ又は予測値102よりも小さい場合に、測定対象成分に関する値として測定値が確定される。 Specifically, the component to be measured includes a solution resistance component 64, and the predicted value 102 of the component to be measured is a real component value. Therefore, when the smallest value among the real component values of each of the acquired actual impedances 90, 92, 94 is the same as or smaller than the predicted value 102, the measured value is determined as the value related to the component to be measured.

なお、測定対象成分は溶液抵抗成分64を含む。このため、測定対象成分に関する値として確定された測定値は、測定対象成分である溶液抵抗成分64のインピーダンスを示す。 The components to be measured include the solution resistance component 64. Therefore, the measured value determined as the value related to the components to be measured indicates the impedance of the solution resistance component 64, which is the component to be measured.

なお、取得処理装置14は、他の処理において、確定したインピーダンスを表示部38に表示してもよい。また、取得処理装置14は、確定したインピーダンスを、外部通信部44を介して、他の装置に送信してもよい。 In addition, the acquisition processing device 14 may display the determined impedance on the display unit 38 in other processing. In addition, the acquisition processing device 14 may transmit the determined impedance to another device via the external communication unit 44.

また、取得処理装置14は、他の処理において、確定したインピーダンスを用いて蓄電デバイス16の特性を評価してもよい。具体的に説明すると、取得処理装置14は、確定したインピーダンスに基づいて、蓄電デバイス16をランク分けし、その結果を表示部38に表示したり、蓄電デバイス16の品質を表示したりしてもよい。 The acquisition processing device 14 may also use the determined impedance to evaluate the characteristics of the power storage device 16 in other processes. Specifically, the acquisition processing device 14 may rank the power storage device 16 based on the determined impedance and display the result on the display unit 38 or display the quality of the power storage device 16.

ステップS5において、プロセッサ30は、各実測インピーダンス90、92、94の最小値が予測値102よりも大きいと判断した場合、周波数算出処理を実行する(ステップS7)。 If in step S5 the processor 30 determines that the minimum value of each measured impedance 90, 92, 94 is greater than the predicted value 102, it executes a frequency calculation process (step S7).

このステップS7は、ステップS4で求めた予測値102が最小の値であるか否かの検証に用いたり、信頼度区間などの統計処理に用いたりするために予測値102における予測周波数Faを算出する。 This step S7 calculates the predicted frequency Fa of the predicted value 102 in order to verify whether the predicted value 102 obtained in step S4 is a minimum value or to use it for statistical processing such as confidence intervals.

図12は、一実施形態に係る周波数算出処理の処理手順を示すフローチャートであり、図12を用いて周波数算出処理の動作を説明する。 Figure 12 is a flowchart showing the processing steps of the frequency calculation process according to one embodiment, and the operation of the frequency calculation process will be explained using Figure 12.

周波数算出処理は、予測値102の予測周波数Faを算出する処理である。周波数算出処理において、プロセッサ30は、取得した各実測インピーダンス90、92、94の実数成分の値R1、R2、R3と、各値R1、R2、R3を取得した各周波数F1、F2、F3と、に基づいて一般多項式を算出する。そして、プロセッサ30は、算出した多項式を実数成分多項式とする(ステップSB1)。なお、実数成分多項式の次数は選択可能とする。 The frequency calculation process is a process for calculating the predicted frequency Fa of the predicted value 102. In the frequency calculation process, the processor 30 calculates a general polynomial based on the real component values R1, R2, R3 of the acquired measured impedances 90, 92, 94 and the frequencies F1, F2, F3 at which the values R1, R2, R3 were acquired. The processor 30 then sets the calculated polynomial as a real component polynomial (step SB1). Note that the degree of the real component polynomial is selectable.

そして、プロセッサ30は、算出した実数成分多項式に、ステップS4求めた実数成分の値Rtを代入して、実数成分の値Rtとなる実数成分の周波数Frを算出する(ステップSB2)。 Then, the processor 30 substitutes the real component value Rt determined in step S4 into the calculated real component polynomial to calculate the frequency Fr of the real component that results in the real component value Rt (step SB2).

また、プロセッサ30は、取得した各実測インピーダンス90、92、94の虚数成分の値X1、X2、X3と、各値X1、X2、X3を取得した各周波数F1、F2、F3と、によって一般多項式を算出する。そして、プロセッサ30は、算出した多項式を虚数成分多項式とする(ステップSB3)。なお、虚数成分多項式の次数は選択可能とする。 The processor 30 also calculates a general polynomial using the imaginary component values X1, X2, and X3 of the acquired measured impedances 90, 92, and 94, and the frequencies F1, F2, and F3 at which the values X1, X2, and X3 were acquired. The processor 30 then sets the calculated polynomial as an imaginary component polynomial (step SB3). Note that the degree of the imaginary component polynomial is selectable.

そして、プロセッサ30は、算出した虚数成分多項式に、ステップS3で求めた虚数成分の値Xtを代入して、虚数成分の値Xtとなる虚数成分の周波数Fxを算出する(ステップSB4)。 Then, the processor 30 substitutes the imaginary component value Xt calculated in step S3 into the calculated imaginary component polynomial to calculate the frequency Fx of the imaginary component that results in the imaginary component value Xt (step SB4).

次に、プロセッサ30は、測定した各実測インピーダンス90、92、94の実数成分の値の最大値から最小値を減算して実数成分の差分ΔRを取得する(ステップSB5)。また、プロセッサ30は、測定した各実測インピーダンス90、92、94の虚数成分の値の最大値から最小値を減算して虚数成分の差分ΔXを取得する(ステップSB6)。 Next, the processor 30 subtracts the minimum value from the maximum value of the real component of each of the measured impedances 90, 92, and 94 to obtain the difference in the real components ΔR (step SB5). The processor 30 also subtracts the minimum value from the maximum value of the imaginary component of each of the measured impedances 90, 92, and 94 to obtain the difference in the imaginary components ΔX (step SB6).

そして、プロセッサ30は、実数成分の周波数Frに実数成分の差分ΔRを乗算して(F’r=Fr×ΔR)、Frに重みづけを行う(ステップSB7)。また、プロセッサ30は、虚数成分の周波数Fxに虚数成分の差分ΔXを乗算して(F’x=Fx×ΔX)、Fxに重みづけを行う(ステップSB8)。 Then, the processor 30 multiplies the frequency of the real component Fr by the difference ΔR between the real components (F'r = Fr x ΔR) to weight Fr (step SB7). The processor 30 also multiplies the frequency of the imaginary component Fx by the difference ΔX between the imaginary components (F'x = Fx x ΔX) to weight Fx (step SB8).

プロセッサ30は、重みづけを行った実数成分の周波数F’rと虚数成分の周波数F’xとの加算値を、実数成分の差分ΔRに虚数成分の差分ΔXを加えた加算値で除算して、加重平均化された予測周波数Faを算出し(ステップSB9)、測定処理へ戻る。予測周波数Faの計算式を次に示す。 The processor 30 divides the sum of the weighted real component frequency F'r and the weighted imaginary component frequency F'x by the sum of the real component difference ΔR and the imaginary component difference ΔX to calculate the weighted average predicted frequency Fa (step SB9), and returns to the measurement process. The calculation formula for the predicted frequency Fa is as follows:

Fa=(F’r+F’x)/(ΔR+ΔX) Fa=(F'r+F'x)/(ΔR+ΔX)

これにより、予測値102の実数成分から求められる実数成分の周波数F’rと予測値102の虚数成分から求められる虚数成分の周波数F’xとに基づいて、予測周波数Faが算出される。なお、Faは、FxとFrとの平均値でもよい。 As a result, the predicted frequency Fa is calculated based on the frequency F'r of the real component calculated from the real component of the predicted value 102 and the frequency F'x of the imaginary component calculated from the imaginary component of the predicted value 102. Note that Fa may be the average value of Fx and Fr.

予測周波数Faは、実数成分の周波数F’r及び虚数成分の周波数F’xに、ステップS1で取得した実数成分の最大値と最小値との差分ΔR及び虚数成分の最大値と最小値との差分ΔXに基づいて重みづけをした加重平均値である。 The predicted frequency Fa is a weighted average of the real component frequency F'r and the imaginary component frequency F'x, weighted based on the difference ΔR between the maximum and minimum real component values and the difference ΔX between the maximum and minimum imaginary component values obtained in step S1.

測定処理において、図10に示したように、プロセッサ30は、ステップS2で求めた極値の実数成分のR値における周波数と、最小分解能の桁を切り捨て等で処理した予測周波数Faと比較して、実測可能であるか否かを判断する(ステップS8)。 In the measurement process, as shown in FIG. 10, the processor 30 compares the frequency at the R value of the real component of the extreme value obtained in step S2 with the predicted frequency Fa processed by truncating the minimum resolution digits, etc., to determine whether or not it is actually measurable (step S8).

ステップS8の判断において、ステップS2で求めた極値の実数成分のR値における周波数と、最小分解能の桁を切り捨て等で処理した予測周波数Faとが等しい場合、新たな値を得られない。このため、現時点における各実測インピーダンス90、92、94の最小値を測定対象成分のインピーダンスとして確定し(ステップS6)、測定処理を終了する。 In the judgment of step S8, if the frequency at the R value of the real component of the extreme value obtained in step S2 is equal to the predicted frequency Fa processed by truncating the minimum resolution digits, etc., a new value cannot be obtained. Therefore, the minimum value of each of the currently measured impedances 90, 92, 94 is determined as the impedance of the component to be measured (step S6), and the measurement process ends.

ステップS8の判断において、ステップS2で求めた極値の実数成分のR値における周波数と、最小分解能の桁を切り捨て等で処理した予測周波数Faとが異なる場合、算出した予測周波数Faで新たな値を得ることができる。このため、予測周波数Faで蓄電デバイス16を再度測定して追加インピーダンスを取得する(ステップS9)。なお、インピーダンスの取得方法は、ステップS1と同じ手順で行う。 In the judgment of step S8, if the frequency at the R value of the real component of the extreme value obtained in step S2 differs from the predicted frequency Fa processed by truncating the minimum resolution digits, etc., a new value can be obtained at the calculated predicted frequency Fa. Therefore, the power storage device 16 is measured again at the predicted frequency Fa to obtain additional impedance (step S9). The impedance is obtained using the same procedure as in step S1.

これにより、測定値が予測値102よりも測定対象成分のインピーダンスから遠い場合には、周波数帯域のうち予測値102に対応する予測周波数Faで測定されたインピーダンスが追加インピーダンスとして取得される。 As a result, when the measured value is farther from the impedance of the component to be measured than the predicted value 102, the impedance measured at the predicted frequency Fa in the frequency band that corresponds to the predicted value 102 is obtained as the additional impedance.

そして、プロセッサ30は、各実測インピーダンス90、92、94の最小値が、予測周波数Faで測定した追加インピーダンスよりも大きいか否かを判断する(ステップS10)。 Then, the processor 30 determines whether the minimum value of each measured impedance 90, 92, 94 is greater than the additional impedance measured at the predicted frequency Fa (step S10).

ステップS10において、プロセッサ30は、各実測インピーダンス90、92、94の最小値が予測周波数Faで取得した追加インピーダンスと同じ又は追加インピーダンスよりも小さいと判断した場合、ステップS6へ分岐する。ステップS6において、プロセッサ30は、各実測インピーダンス90、92、94の最小値を測定対象成分に関する測定値として確定して測定処理を終了する(ステップS6)。 In step S10, if the processor 30 determines that the minimum value of each measured impedance 90, 92, 94 is equal to or smaller than the additional impedance acquired at the predicted frequency Fa, the processor 30 branches to step S6. In step S6, the processor 30 determines the minimum value of each measured impedance 90, 92, 94 as the measured value for the component to be measured, and ends the measurement process (step S6).

一方、ステップS10において、プロセッサ30は、各実測インピーダンス90、92、94の最小値が、予測周波数Faで測定した追加インピーダンスよりも大きいと判断した場合、追加インピーダンスを最小値とする(ステップS11)。 On the other hand, in step S10, if the processor 30 determines that the minimum value of each measured impedance 90, 92, 94 is greater than the additional impedance measured at the predicted frequency Fa, the additional impedance is set to the minimum value (step S11).

そして、プロセッサ30は、「カウンタ」に「1」を加えて、ステップS2-2へ移行する。 Then, the processor 30 adds "1" to the "counter" and proceeds to step S2-2.

ステップS2-2において、プロセッサ30は、「カウンタ」がn未満であるか否かを判断する(ステップS2-2)。ステップS2-2において、プロセッサ30は、「カウンタ」がnであると判断した場合、現時点での全実測インピーダンスの最小値を測定対象成分に関する測定値として確定して測定処理を終了する(ステップS6)。 In step S2-2, the processor 30 determines whether the "counter" is less than n (step S2-2). If the processor 30 determines in step S2-2 that the "counter" is n, it determines the minimum value of all the actual measured impedances at this point in time as the measured value for the component to be measured, and ends the measurement process (step S6).

これにより、ステップS3からステップS12の繰り返し回数をn回に制限することができる。 This allows the number of repetitions of steps S3 to S12 to be limited to n.

一方、ステップS2-2において、プロセッサ30は、「カウンタ」がn未満であると判断した場合、ステップS3へ移行する。 On the other hand, if in step S2-2, the processor 30 determines that the "counter" is less than n, it proceeds to step S3.

ステップS3において、プロセッサ30は、追加インピーダンスが各実測インピーダンス90、92、94に加えられた全実測インピーダンスを用いた前述した方法で二次曲線を算出する(ステップS3)。 In step S3, the processor 30 calculates a quadratic curve using the above-described method using the total measured impedance with the additional impedance added to each measured impedance 90, 92, 94 (step S3).

また、プロセッサ30は、取得した二次曲線から測定対象成分の予測値を新たに算出する(ステップS4)。 The processor 30 also calculates new predicted values of the components to be measured from the acquired quadratic curve (step S4).

これにより、追加インピーダンスの実数成分及び虚数成分と各実測インピーダンス90、92、94の実数成分及び虚数成分とに基づいて、近似式である二次曲線を求め、この近似式に基づいて、測定対象成分の予測値を新たに算出する。 This allows a quadratic curve, which is an approximation equation, to be calculated based on the real and imaginary components of the additional impedance and the real and imaginary components of each measured impedance 90, 92, and 94, and a new predicted value of the component to be measured is calculated based on this approximation equation.

そして、プロセッサ30は、全実測インピーダンスの最小値が、新たに算出された予測値よりも大きいか否かを判断する(ステップS5)。 Then, the processor 30 determines whether the minimum value of all measured impedances is greater than the newly calculated predicted value (step S5).

ステップS5において、全実測インピーダンスの最小値が新たに算出された予測値と同じ又は予測値よりも小さいと判断した場合、全実測インピーダンスの最小値を測定対象成分に関する値として確定して測定処理を終了する(ステップS6)。 If it is determined in step S5 that the minimum value of all measured impedances is equal to or smaller than the newly calculated predicted value, the minimum value of all measured impedances is determined as the value related to the component to be measured, and the measurement process is terminated (step S6).

これにより、全実測インピーダンスのうち測定対象成分に最も近い値が、ステップS4で新たに算出した予測値と同じ又は予測値よりも測定対象成分のインピーダンスに近い場合に、測定対象成分に関する値として測定値が確定される。 As a result, if the value closest to the component to be measured among all the measured impedances is the same as the predicted value newly calculated in step S4 or is closer to the impedance of the component to be measured than the predicted value, the measured value is determined as the value related to the component to be measured.

また、ステップS5において、全実測インピーダンスの最小値が、新たに算出した予測値よりも大きいと判断した場合、ステップS7へ移行し、測定対象成分のインピーダンスが確定するまで、各ステップを繰り返す。 Also, if it is determined in step S5 that the minimum value of all measured impedances is greater than the newly calculated predicted value, the process proceeds to step S7, and each step is repeated until the impedance of the component to be measured is determined.

(作用及び効果)
次に、本実施形態による作用効果について説明する。
(Action and Effects)
Next, the effects of this embodiment will be described.

本実施形態における取得処理装置14は、蓄電デバイス16のインピーダンスを測定する取得処理装置14である。取得処理装置14は、所定の周波数帯域において計測装置12を用いて三点以上の周波数で測定された測定対象成分のインピーダンスを取得する取得部50を含む。取得処理装置14は、取得部50で取得した三点以上のインピーダンスの実数成分又は虚数成分の最小値である極値を求める決定部51を含む。取得処理装置14は、三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて測定対象成分の予測値を算出する算出部52を含む。取得処理装置14は、極値が、予測値と同じ又は予測値よりも測定対象成分のインピーダンスに近づくように突出する場合には、極値を測定対象成分に関する測定値として確定する確定部54を含む。 The acquisition processing device 14 in this embodiment is an acquisition processing device 14 that measures the impedance of the power storage device 16. The acquisition processing device 14 includes an acquisition unit 50 that acquires the impedance of the measurement target component measured at three or more frequencies using the measurement device 12 in a predetermined frequency band. The acquisition processing device 14 includes a determination unit 51 that determines an extreme value that is the minimum value of the real component or imaginary component of the impedance at three or more points acquired by the acquisition unit 50. The acquisition processing device 14 includes a calculation unit 52 that calculates a predicted value of the measurement target component based on an approximation formula that represents the relationship between the real component and the imaginary component of the impedance at three or more points. The acquisition processing device 14 includes a determination unit 54 that determines the extreme value as a measured value related to the measurement target component when the extreme value protrudes so as to be the same as the predicted value or closer to the impedance of the measurement target component than the predicted value.

この構成によれば、測定して取得した三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて測定対象成分の予測値を算出する。そして、三点以上の実数成分又は虚数成分の極値が、近似式に基づく予測値と同じ又は予測値よりも測定対象成分のインピーダンスに近づくように突出する場合に、極値を測定対象成分に関する測定値として確定する。 According to this configuration, a predicted value of the component to be measured is calculated based on an approximation formula that expresses the relationship between the real and imaginary components of the impedance at three or more points obtained by measurement. Then, when the extreme values of the real or imaginary components at three or more points stand out so as to be the same as the predicted value based on the approximation formula or to be closer to the impedance of the component to be measured than the predicted value, the extreme values are determined as the measured value for the component to be measured.

そして、この測定値は、実際に測定された実測値によって定められており、この実測値は、蓄電デバイスの特性を示す値である。 This measurement value is determined based on an actual measurement, and this actual measurement value indicates the characteristics of the energy storage device.

したがって、測定対象となる蓄電デバイスの特性を示す測定結果を測定値として取得することが可能となる。 Therefore, it is possible to obtain measurement results that indicate the characteristics of the energy storage device being measured as measured values.

具体的に説明すると、非線形アルゴリズムで等価回路解析を行う場合、特性を合わせ込む際に、どのパラメータを固定するかによって結果が変わる虞がある。このため、残差値が少なくなったとしても、その値が蓄電デバイスの特性を示すものであることを証明することができない。 To be more specific, when performing equivalent circuit analysis using a nonlinear algorithm, there is a risk that the results will change depending on which parameters are fixed when fitting the characteristics. Therefore, even if the residual value is reduced, it cannot be proven that the value represents the characteristics of the energy storage device.

しかし、本実施形態では、予測値102よりも測定対象成分に近い実測値を、測定対象成分に関する測定値として確定する。このため、実際に測定された値を用いて測定対象成分に関する測定値を定めることができる。 However, in this embodiment, the actual measured value that is closer to the component to be measured than the predicted value 102 is determined as the measured value for the component to be measured. Therefore, the measured value for the component to be measured can be determined using the actually measured value.

また、この構成によれば、取得した三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて測定対象成分の予測値102を算出する。そして、取得した実数成分又は虚数成分の各値のうち測定対象成分に最も近い値が予測値102と同じ又は予測値102よりも測定対象成分のインピーダンスに近い場合に、測定対象成分に関する値として測定値を確定する。 In addition, according to this configuration, a predicted value 102 of the component to be measured is calculated based on an approximation formula that expresses the relationship between the real and imaginary components of the impedance at three or more points that have been acquired. Then, if the value closest to the component to be measured among the acquired real or imaginary components is the same as the predicted value 102 or is closer to the impedance of the component to be measured than the predicted value 102, the measured value is determined as the value related to the component to be measured.

これにより、少なくとも三点のインピーダンスを処理することで、測定対象成分のインピーダンスを取得することができる。このため、多数の測定結果を演算処理しなければならない場合と比較して、プロセッサ30による演算処理の負担軽減が可能となる。 As a result, the impedance of the component to be measured can be obtained by processing the impedance at at least three points. This reduces the burden of calculation processing on the processor 30 compared to when a large number of measurement results must be processed.

また、多数の測定結果を要する場合と比較して、測定時間を短縮することができる。さらに、多数の測定結果の入力を要する場合と比較して、取得処理装置14への測定結果の入力処理が軽減される。これにより、測定対象成分に関する値を効率よく算出することができる。 In addition, the measurement time can be shortened compared to when a large number of measurement results are required. Furthermore, the input process of the measurement results to the acquisition processing device 14 is reduced compared to when a large number of measurement results are required to be input. This allows values related to the components to be measured to be calculated efficiently.

また、等価回路解析を行う場合のように、高度な技術を要する等価回路の設計及び初期パラメータ値の設定などが不要となる。これにより、熟練を要しないものであっても、測定対象成分の値を安定的に取得することができる。 In addition, unlike equivalent circuit analysis, there is no need to design an equivalent circuit or set initial parameter values, which require advanced techniques. This makes it possible to stably obtain values for the components to be measured, even if one does not require advanced skills.

また、本実施形態の処理方法は、蓄電デバイス16のインピーダンスを測定する処理方法である。処理方法は、所定の周波数帯域において三点以上の周波数で測定された測定対象成分のインピーダンスを取得する取得工程を含む。処理方法は、取得工程で取得した三点以上のインピーダンスの実数成分又は虚数成分の極値を求める決定工程を含む。処理方法は、三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて測定対象成分の予測値を算出する算出工程を含む。処理方法は、極値が、予測値と同じ又は予測値よりも測定対象成分のインピーダンスに近づくように突出する場合には、極値を測定対象成分に関する測定値として確定する確定工程を含む。 The processing method of this embodiment is a processing method for measuring the impedance of the power storage device 16. The processing method includes an acquisition step of acquiring the impedance of the measurement target component measured at three or more frequencies in a predetermined frequency band. The processing method includes a determination step of determining the extreme values of the real or imaginary components of the impedance at three or more points acquired in the acquisition step. The processing method includes a calculation step of calculating a predicted value of the measurement target component based on an approximation formula expressing the relationship between the real and imaginary components of the impedance at three or more points. The processing method includes a determination step of determining the extreme value as the measured value for the measurement target component when the extreme value protrudes so as to be the same as the predicted value or closer to the impedance of the measurement target component than the predicted value.

さらに、本実施形態の処理プログラムは、プロセッサ30を取得処理装置14の各手段である各部50、52、54として機能させる。 Furthermore, the processing program of this embodiment causes the processor 30 to function as each of the units 50, 52, and 54, which are each means of the acquisition processing device 14.

そして、本実施形態の測定システム10は、蓄電デバイス16のインピーダンスを測定するための複数の装置を有する測定システムである。測定システム10は、所定の周波数帯域において三点以上の周波数で測定された測定対象成分のインピーダンスを取得する取得部50と、取得部50で取得した三点以上のインピーダンスの実数成分又は虚数成分の極値を求める決定部51とを備える。測定システム10は、三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて測定対象成分の予測値を算出する算出部52を備える。測定システム10は、極値が、予測値と同じ又は予測値よりも測定対象成分のインピーダンスに近づくように突出する場合には、極値を測定対象成分に関する測定値として確定する確定部54を備える。 The measurement system 10 of this embodiment is a measurement system having multiple devices for measuring the impedance of the power storage device 16. The measurement system 10 includes an acquisition unit 50 that acquires the impedance of the measurement target component measured at three or more frequencies in a predetermined frequency band, and a determination unit 51 that determines the extreme values of the real or imaginary components of the impedance at the three or more points acquired by the acquisition unit 50. The measurement system 10 includes a calculation unit 52 that calculates a predicted value of the measurement target component based on an approximation formula that represents the relationship between the real and imaginary components of the impedance at the three or more points. The measurement system 10 includes a determination unit 54 that determines the extreme value as the measured value for the measurement target component when the extreme value protrudes so as to be the same as the predicted value or closer to the impedance of the measurement target component than the predicted value.

これらの処理方法、処理プログラム、及び測定システム10においても、前述した作用効果を奏することができる。 These processing methods, processing programs, and measurement system 10 can also achieve the effects described above.

また、本実施形態の取得処理装置14において、算出手段である算出部52は、三点以上のインピーダンスの実数成分及び虚数成分の関係を二次関数で表した近似式の頂点を予測値102として算出する。 In addition, in the acquisition processing device 14 of this embodiment, the calculation unit 52, which is a calculation means, calculates the apex of an approximation equation that expresses the relationship between the real components and the imaginary components of the impedance at three or more points as a quadratic function as a predicted value 102.

この構成によれば、予測値102を二次関数の頂点に基づいて算出することができる。このため、予測値102が収束するまで近似式を繰り返し演算しなければならない場合と比較して、プロセッサ30による演算処理の負担をさらに軽減することができる。 With this configuration, the predicted value 102 can be calculated based on the vertex of the quadratic function. Therefore, the burden of calculation processing by the processor 30 can be further reduced compared to a case where the approximation formula must be repeatedly calculated until the predicted value 102 converges.

さらに、本実施形態の取得処理装置14において、算出手段である算出部52は、取得手段である取得部50で取得した実数成分及び虚数成分の各値に基づいて二次曲線を算出し、二次曲線が突出方向に最も突出した曲点に基づいて予測値102を定める。 Furthermore, in the acquisition processing device 14 of this embodiment, the calculation unit 52, which is a calculation means, calculates a quadratic curve based on the values of the real components and imaginary components acquired by the acquisition unit 50, which is an acquisition means, and determines a predicted value 102 based on the inflection point where the quadratic curve protrudes most in the protruding direction.

この構成によれば、二次曲線において最も突出した曲点を特定することにより、この関数を頂点として予測値102を正確に定めることができる。 With this configuration, by identifying the most prominent inflection point in the quadratic curve, it is possible to accurately determine the predicted value 102 with this function as the apex.

また、本実施形態の取得処理装置14において、取得手段である取得部50は、予測値102が極値よりも測定対象成分のインピーダンスに近づくように突出する場合には、周波数帯域のうち予測値102に対応する予測周波数で測定された追加インピーダンスを取得する。算出手段である算出部52は、追加インピーダンスの実数成分及び虚数成分と三点以上のインピーダンスの実数成分及び虚数成分とに基づいて近似式を求め、近似式に基づいて測定対象成分の予測値を新たに算出する。確定手段である確定部54は、測定対象成分のインピーダンスに最も近い極値が、算出部52が新たに算出した予測値と同じ又は予測値よりも測定対象成分のインピーダンスに近い場合に、極値を測定値として確定する。 In addition, in the acquisition processing device 14 of this embodiment, the acquisition unit 50, which is an acquisition means, acquires the additional impedance measured at a predicted frequency in the frequency band corresponding to the predicted value 102 when the predicted value 102 protrudes closer to the impedance of the measurement target component than the extreme value. The calculation unit 52, which is a calculation means, finds an approximation formula based on the real and imaginary components of the additional impedance and the real and imaginary components of the impedance at three or more points, and calculates a new predicted value of the measurement target component based on the approximation formula. The determination unit 54, which is a determination means, determines the extreme value as the measured value when the extreme value closest to the impedance of the measurement target component is the same as the newly calculated predicted value by the calculation unit 52 or is closer to the impedance of the measurement target component than the predicted value.

この構成によれば、測定値が予測値102よりも測定対象成分のインピーダンスから遠い場合、予測値102に対応する予測周波数Faで測定された追加インピーダンスを追加して新たな予測値102で測定対象成分に関する値を確定することができる。 With this configuration, if the measured value is farther from the impedance of the component to be measured than the predicted value 102, an additional impedance measured at the predicted frequency Fa corresponding to the predicted value 102 can be added to determine the value related to the component to be measured with the new predicted value 102.

また、本実施形態の取得処理装置14において、取得手段である取得部50は、予測値102の実数成分から求められる実数成分周波数Frと予測値102の虚数成分から求められる虚数成分周波数Fxとに基づいて予測周波数Faを算出する。 In addition, in the acquisition processing device 14 of this embodiment, the acquisition unit 50, which is an acquisition means, calculates the predicted frequency Fa based on the real component frequency Fr obtained from the real component of the predicted value 102 and the imaginary component frequency Fx obtained from the imaginary component of the predicted value 102.

この構成によれば、予測値102の実数成分から求められる実数成分周波数Frと予測値102の虚数成分から求められる虚数成分周波数Fxと二つの周波数を用いて予測周波数Faを算出する。このため、一つの周波数を用いる場合と比べて、予測周波数Faの誤差範囲を狭くすることができる。 According to this configuration, the predicted frequency Fa is calculated using two frequencies: the real component frequency Fr calculated from the real component of the predicted value 102, and the imaginary component frequency Fx calculated from the imaginary component of the predicted value 102. Therefore, the error range of the predicted frequency Fa can be narrowed compared to when one frequency is used.

よって、予測周波数Faで測定される追加インピーダンスを、測定対象成分のインピーダンスの真値に近づけることが可能となる。 As a result, it is possible to bring the additional impedance measured at the predicted frequency Fa closer to the true value of the impedance of the component being measured.

さらに、本実施形態の取得処理装置14において、予測周波数Faは、実数成分周波数Frと虚数成分周波数Fxの加重平均値である。加重平均値は、実数成分周波数Fr及び虚数成分周波数Fxに、取得手段である取得部50で取得した実数成分の最大値と最小値との差分ΔR及び虚数成分の最大値と最小値との差分ΔXに基づいて重みづけした値である。 Furthermore, in the acquisition processing device 14 of this embodiment, the predicted frequency Fa is a weighted average value of the real component frequency Fr and the imaginary component frequency Fx. The weighted average value is a value obtained by weighting the real component frequency Fr and the imaginary component frequency Fx based on the difference ΔR between the maximum and minimum values of the real component acquired by the acquisition unit 50, which is an acquisition means, and the difference ΔX between the maximum and minimum values of the imaginary component.

この構成によれば、異なる周波数における実数成分の値の変化量と虚数成分の値の変化量とを考慮して予測周波数Faを求めることで、計測装置12による測定値のばらつきの影響を減らすことが可能となる。 With this configuration, the predicted frequency Fa is calculated by taking into account the amount of change in the value of the real component and the amount of change in the value of the imaginary component at different frequencies, making it possible to reduce the effects of variability in the measured values from the measuring device 12.

一例を挙げて説明すると、電解液抵抗を算出する周波数帯域において、実数成分の差分ΔRよりも虚数成分の差分ΔXが大きくなる。このため、虚数成分の値が実数成分の値よりも正確に測定された可能性が高い。よって、予測周波数Faを、前述した加重平均値とすることで、計測装置12による測定値のばらつきの影響を抑制することが可能となる。 To give an example, in the frequency band for calculating the electrolyte resistance, the difference ΔX of the imaginary components is larger than the difference ΔR of the real components. Therefore, it is highly likely that the value of the imaginary components is measured more accurately than the value of the real components. Therefore, by setting the predicted frequency Fa to the weighted average value described above, it is possible to suppress the effects of variability in the measured values by the measuring device 12.

また、本実施形態の取得処理装置14において、蓄電デバイス16の測定対象成分は、溶液抵抗成分64を含み、測定対象成分の予測値102は、実数成分の値である。確定手段である確定部54は、三点以上のインピーダンスの実数成分の各値のうち最も小さい測定値が、予測値102と同じ又は予測値102よりも小さい場合には、測定対象成分に関する値として測定値を確定する。 In addition, in the acquisition processing device 14 of this embodiment, the measurement target component of the electricity storage device 16 includes a solution resistance component 64, and the predicted value 102 of the measurement target component is a real component value. When the smallest measurement value among the real component values of the impedance at three or more points is the same as or smaller than the predicted value 102, the determination unit 54, which is a determination means, determines the measurement value as the value related to the measurement target component.

この構成によれば、溶液抵抗成分64に関する値を精度よく取得することができる。 This configuration allows the value of the solution resistance component 64 to be obtained with high accuracy.

ここで、溶液抵抗成分64を測定できる領域は、他の領域と比較して、蓄電デバイス16を構成する等価回路において、隣合う回路素子の影響を受けにくい。このため、溶液抵抗成分に関する値の取得に適している。 The area where the solution resistance component 64 can be measured is less affected by adjacent circuit elements in the equivalent circuit that constitutes the electricity storage device 16 than other areas. Therefore, it is suitable for obtaining values related to the solution resistance component.

また、測定対象成分に関する値として確定した測定値は、高周波帯で測定している。このため、ケーブル抵抗成分62及び反応抵抗成分66が測定値に与える影響は相対的に小さくなる。これにより、測定対象成分に関する値として確定した測定値が溶液抵抗成分64のインピーダンスである確度は高い。 In addition, the measured value determined as the value related to the component to be measured is measured in the high frequency band. Therefore, the influence of the cable resistance component 62 and the reaction resistance component 66 on the measured value is relatively small. As a result, there is a high degree of certainty that the measured value determined as the value related to the component to be measured is the impedance of the solution resistance component 64.

したがって、蓄電デバイス16の電解液の成分及び濃度などの均一性の評価、又は蓄電デバイス16の電極タブの溶接の評価を、より正確に行うことが可能となる。 This makes it possible to more accurately evaluate the uniformity of the electrolyte components and concentration of the electricity storage device 16, or to evaluate the welding of the electrode tabs of the electricity storage device 16.

以上、本発明の実施形態について説明したが、上記実施形態は本発明の適用例の一部を示したに過ぎず、本発明の技術的範囲を上記実施形態の具体的構成に限定する趣旨ではない。 Although the embodiments of the present invention have been described above, the above embodiments merely show some of the application examples of the present invention, and are not intended to limit the technical scope of the present invention to the specific configurations of the above embodiments.

10 測定システム
12 計測装置
14 取得処理装置
16 蓄電デバイス
30 プロセッサ
34 記憶部
50 取得部
52 算出部
54 確定部
64 溶液抵抗成分
66 反応抵抗成分
72 電気二重層容量成分
92 第二実測インピーダンス
100 二次曲線
102 予測値
C1、C2 電気二重層容量
Fa 予測周波数
Rsol 電解液抵抗
ΔR、ΔX 差分
REFERENCE SIGNS LIST 10 Measurement system 12 Measurement device 14 Acquisition processing device 16 Electricity storage device 30 Processor 34 Memory unit 50 Acquisition unit 52 Calculation unit 54 Determination unit 64 Solution resistance component 66 Reaction resistance component 72 Electric double layer capacitance component 92 Second measured impedance 100 Quadratic curve 102 Predicted value C1, C2 Electric double layer capacitance Fa Predicted frequency Rsol Electrolyte resistance ΔR, ΔX Difference

Claims (10)

蓄電デバイスのインピーダンスを取得処理する取得処理装置であって、
所定の周波数帯域において計測装置を用いて三点以上の周波数で測定された測定対象成分のインピーダンスを取得する取得手段と、
前記取得手段で取得した前記三点以上のインピーダンスの実数成分又は虚数成分の極値を求める決定手段と、
前記三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて前記測定対象成分の予測値を算出する算出手段と、
前記極値が、前記予測値と同じ又は前記予測値よりも前記測定対象成分のインピーダンスに近づくように突出する場合には、前記極値を前記測定対象成分に関する測定値として確定する確定手段と、
を含む取得処理装置。
An acquisition processing device that acquires and processes impedance of a power storage device,
An acquisition means for acquiring impedance of a measurement target component measured at three or more frequencies using a measurement device in a predetermined frequency band;
A determination means for determining an extreme value of a real component or an imaginary component of the impedance at the three or more points acquired by the acquisition means;
a calculation means for calculating a predicted value of the measurement target component based on an approximation formula expressing a relationship between real components and imaginary components of the impedance at the three or more points;
a determination means for determining, when the extreme value is equal to or closer to the impedance of the measured component than the predicted value, the extreme value as a measured value related to the measured component;
an acquisition processing device including:
請求項1に記載の取得処理装置であって、
前記算出手段は、前記取得手段によって取得された前記三点以上のインピーダンスの実数成分及び虚数成分の関係を二次関数で表した前記近似式の頂点を前記予測値として算出する、
取得処理装置。
The acquisition processing device according to claim 1 ,
The calculation means calculates, as the predicted value, a peak of the approximation equation obtained by expressing a relationship between real components and imaginary components of the impedances at the three or more points acquired by the acquisition means as a quadratic function.
Acquisition processing device.
請求項2に記載の取得処理装置であって、
前記算出手段は、前記取得手段によって取得された前記三点以上のインピーダンスの実数成分及び虚数成分の各値に基づいて二次曲線を算出し、前記二次曲線が突出方向に最も突出した曲点に基づいて前記予測値を定める、
取得処理装置。
The acquisition processing device according to claim 2 ,
The calculation means calculates a quadratic curve based on each value of a real component and an imaginary component of the impedance at the three or more points acquired by the acquisition means, and determines the predicted value based on a bending point at which the quadratic curve is most prominent in a protruding direction.
Acquisition processing device.
請求項1から請求項3のいずれか一項に記載の取得処理装置であって、
前記取得手段は、前記予測値が前記極値よりも前記測定対象成分のインピーダンスに近づくように突出する場合には、前記周波数帯域のうち前記予測値に対応する予測周波数で測定された追加インピーダンスを取得し、
前記算出手段は、前記追加インピーダンスの実数成分及び虚数成分と前記三点以上のインピーダンスの実数成分及び虚数成分とに基づいて前記近似式を求め、当該近似式に基づいて前記測定対象成分の前記予測値を新たに算出し、
前記確定手段は、前記測定対象成分のインピーダンスに最も近い前記極値が、前記算出手段が新たに算出した前記予測値と同じ又は前記予測値よりも前記測定対象成分のインピーダンスに近い場合に、前記極値を前記測定値として確定する、
取得処理装置。
The acquisition processing device according to any one of claims 1 to 3,
When the predicted value is closer to the impedance of the measurement target component than the extreme value, the acquisition means acquires an additional impedance measured at a predicted frequency corresponding to the predicted value within the frequency band;
the calculation means determines the approximation formula based on the real and imaginary components of the additional impedance and the real and imaginary components of the impedances at the three or more points, and newly calculates the predicted value of the measurement target component based on the approximation formula;
The determination means determines the extreme value as the measured value when the extreme value closest to the impedance of the measured component is the same as the predicted value newly calculated by the calculation means or is closer to the impedance of the measured component than the predicted value.
Acquisition processing device.
請求項4に記載の取得処理装置であって、
前記取得手段は、前記予測値の実数成分から求められる実数成分周波数と前記予測値の虚数成分から求められる虚数成分周波数とに基づいて前記予測周波数を算出する、
取得処理装置。
The acquisition processing device according to claim 4 ,
the obtaining means calculates the predicted frequency based on a real component frequency obtained from a real component of the predicted value and an imaginary component frequency obtained from an imaginary component of the predicted value.
Acquisition processing device.
請求項5に記載の取得処理装置であって、
前記予測周波数は、前記実数成分周波数及び前記虚数成分周波数に、前記三点以上のインピーダンスの前記実数成分の最大値と最小値との差分及び前記虚数成分の最大値と最小値との差分に基づいて重みづけをした加重平均値である、
取得処理装置。
The acquisition processing device according to claim 5 ,
The predicted frequency is a weighted average value obtained by weighting the real component frequency and the imaginary component frequency based on a difference between a maximum value and a minimum value of the real component of the impedance at the three or more points and a difference between a maximum value and a minimum value of the imaginary component.
Acquisition processing device.
請求項1から請求項6のいずれか一項に記載の取得処理装置であって、
前記測定対象成分は、溶液抵抗成分を含み、
前記測定対象成分の前記予測値は、前記実数成分の値であり、
前記確定手段は、前記三点以上の前記実数成分の各値のうち最も小さい値を前記極値とし、前記極値が前記予測値と同じ又は前記予測値よりも小さい場合に、前記極値を前記測定値として確定する、
取得処理装置。
The acquisition processing device according to any one of claims 1 to 6,
The component to be measured includes a solution resistance component,
the predicted value of the measured component is the value of the real component;
the determination means determines the smallest value among the values of the real components of the three or more points as the extreme value, and determines the extreme value as the measurement value when the extreme value is equal to or smaller than the predicted value.
Acquisition processing device.
蓄電デバイスのインピーダンスを測定する処理方法であって、
所定の周波数帯域において三点以上の周波数で測定された測定対象成分のインピーダンスを取得する取得工程と、
前記取得工程で取得した前記三点以上のインピーダンスの実数成分又は虚数成分の極値を求める決定工程と、
前記三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて前記測定対象成分の予測値を算出する算出工程と、
前記極値が、前記予測値と同じ又は前記予測値よりも前記測定対象成分のインピーダンスに近づくように突出する場合には、前記極値を前記測定対象成分に関する測定値として確定する確定工程と、
を含む処理方法。
A processing method for measuring impedance of a power storage device, comprising:
an acquisition step of acquiring impedance of a measurement target component measured at three or more frequencies in a predetermined frequency band;
a determination step of determining extreme values of real components or imaginary components of the impedances at the three or more points acquired in the acquisition step;
a calculation step of calculating a predicted value of the measurement target component based on an approximation formula expressing a relationship between real components and imaginary components of the impedance at the three or more points;
a determination step of determining, when the extreme value is the same as the predicted value or is closer to the impedance of the measured component than the predicted value, the extreme value as a measured value related to the measured component;
A processing method comprising:
プロセッサを、請求項1から請求項7のいずれか一項に記載の取得処理装置の各手段として機能させる処理プログラム。 A processing program that causes a processor to function as each of the means of the acquisition processing device according to any one of claims 1 to 7. 蓄電デバイスのインピーダンスを測定するための複数の装置を有する測定システムであって、
所定の周波数帯域において三点以上の周波数で測定された測定対象成分のインピーダンスを取得する取得部と、
前記取得部で取得した前記三点以上のインピーダンスの実数成分又は虚数成分の極値を求める決定部と、
前記三点以上のインピーダンスの実数成分及び虚数成分の関係を表した近似式に基づいて前記測定対象成分の予測値を算出する算出部と、
前記極値が、前記予測値と同じ又は前記予測値よりも前記測定対象成分のインピーダンスに近づくように突出する場合には、前記極値を前記測定対象成分に関する測定値として確定する確定部と、
を備える測定システム。
A measurement system having a plurality of devices for measuring impedance of a power storage device,
an acquisition unit that acquires impedance of a measurement target component measured at three or more frequencies in a predetermined frequency band;
A determination unit that determines extreme values of real components or imaginary components of the impedances at the three or more points acquired by the acquisition unit;
a calculation unit that calculates a predicted value of the measurement target component based on an approximation formula that represents a relationship between real components and imaginary components of the impedance at the three or more points;
a determination unit that determines the extreme value as a measurement value related to the measured component when the extreme value is the same as the predicted value or is closer to the impedance of the measured component than the predicted value;
A measurement system comprising:
JP2021087400A 2021-05-25 2021-05-25 Acquisition processing device, processing method, processing program, and measurement system Active JP7570970B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021087400A JP7570970B2 (en) 2021-05-25 2021-05-25 Acquisition processing device, processing method, processing program, and measurement system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021087400A JP7570970B2 (en) 2021-05-25 2021-05-25 Acquisition processing device, processing method, processing program, and measurement system

Publications (2)

Publication Number Publication Date
JP2022180739A JP2022180739A (en) 2022-12-07
JP7570970B2 true JP7570970B2 (en) 2024-10-22

Family

ID=84327524

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021087400A Active JP7570970B2 (en) 2021-05-25 2021-05-25 Acquisition processing device, processing method, processing program, and measurement system

Country Status (1)

Country Link
JP (1) JP7570970B2 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001235525A (en) 2000-02-22 2001-08-31 Furukawa Electric Co Ltd:The Lead storage battery characteristic evaluation method and lead storage battery characteristic evaluation device
WO2012073997A1 (en) 2010-11-30 2012-06-07 本田技研工業株式会社 Degradation estimation device for secondary battery
JP2013170862A (en) 2012-02-20 2013-09-02 Hioki Ee Corp Battery selector and battery selection method
US20210123981A1 (en) 2019-10-23 2021-04-29 NOVUM engineerING GmbH Estimating a Battery State of an Electrochemical Battery

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6046418B2 (en) * 2012-08-27 2016-12-14 日置電機株式会社 Equivalent circuit analysis apparatus and equivalent circuit analysis method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001235525A (en) 2000-02-22 2001-08-31 Furukawa Electric Co Ltd:The Lead storage battery characteristic evaluation method and lead storage battery characteristic evaluation device
WO2012073997A1 (en) 2010-11-30 2012-06-07 本田技研工業株式会社 Degradation estimation device for secondary battery
JP2013170862A (en) 2012-02-20 2013-09-02 Hioki Ee Corp Battery selector and battery selection method
US20210123981A1 (en) 2019-10-23 2021-04-29 NOVUM engineerING GmbH Estimating a Battery State of an Electrochemical Battery

Also Published As

Publication number Publication date
JP2022180739A (en) 2022-12-07

Similar Documents

Publication Publication Date Title
US6208147B1 (en) Method of and apparatus for measuring battery capacity by impedance spectrum analysis
JP6145824B2 (en) Battery diagnostic device and battery diagnostic method
EP4024069A1 (en) Method and apparatus for estimating battery temperature, electronic device, and storage medium
JP7502905B2 (en) Battery degradation diagnosis system and diagnosis processing device
JP2019191030A (en) Battery information processing system, battery pack, battery module characteristics evaluation method, and battery pack manufacturing method
EP3812783B1 (en) Estimating a battery state from electrical impedance measurements using convolutional neural network means
KR102649551B1 (en) Battery electric flow test method and system
JP7588479B2 (en) State measurement device, state measurement method, and program
CN107734488A (en) Bluetooth equipment quality detection device and method
CN115469200A (en) Insulation test method, device, equipment and storage medium
CN114325278A (en) Air gap tolerance voltage measuring and calculating method and device, computer equipment and medium
JP7570970B2 (en) Acquisition processing device, processing method, processing program, and measurement system
CN112698212B (en) Estimating battery state from gradients of electrical impedance measurements
CN117321413B (en) Battery water content detection method and device, electronic equipment and storage medium
EP3812782A1 (en) Estimating a temperature of an electrochemical battery
CN108414860A (en) Method for testing performance, device and the readable storage medium storing program for executing of high frequency cable
JP4833788B2 (en) Battery voltage prediction method, program, state monitoring device, and power supply system
US20220349944A1 (en) Battery performance evaluation device, electronic apparatus, charger, and battery performance evaluation method
US20240426929A1 (en) Battery pack management system and battery pack management method
CN113425280B (en) Identification method, measurement device and storage medium for abnormal measurement of impedance of human body
Křivík et al. Impedance methods for SoC determination of lead acid battery cell
CN114285118A (en) Battery parameter correction system, method and storage medium
CN113155901A (en) Uniformity acquisition device, uniformity acquisition method, and program
KR102956623B1 (en) Apparatus and method for diagnosing battery
CN106169944B (en) Automatic parameter adjusting system and automatic parameter adjusting method thereof

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240214

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240930

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241009

R150 Certificate of patent or registration of utility model

Ref document number: 7570970

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150