特許 6306456 診断装置、及び診断方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6306456

(24)【登録日】2018年3月16日

(45)【発行日】2018年4月4日

(54)【発明の名称】診断装置、及び診断方法

(51)【国際特許分類】

   G01N 23/05 20060101AFI20180326BHJP

   G01T 3/00 20060101ALI20180326BHJP

   G01T 1/36 20060101ALI20180326BHJP

   G01T 1/17 20060101ALI20180326BHJP

   G21K 5/02 20060101ALI20180326BHJP

   G21K 1/00 20060101ALI20180326BHJP

   H01M 10/48 20060101ALI20180326BHJP

【ＦＩ】

   G01N23/05

   G01T3/00 Z

   G01T1/36 A

   G01T1/17 H

   G01T3/00 H

   G21K5/02 N

   G21K1/00 N

   H01M10/48 P

【請求項の数】8

【全頁数】15

(21)【出願番号】特願2014-142069(P2014-142069)

(22)【出願日】2014年7月10日

(65)【公開番号】特開2016-17901(P2016-17901A)

(43)【公開日】2016年2月1日

【審査請求日】2016年8月5日

(73)【特許権者】

【識別番号】504173471

【氏名又は名称】国立大学法人北海道大学

(73)【特許権者】

【識別番号】000003207

【氏名又は名称】トヨタ自動車株式会社

(74)【代理人】

【識別番号】100103894

【弁理士】

【氏名又は名称】家入 健

(72)【発明者】

【氏名】加美山 隆

(72)【発明者】

【氏名】佐藤 博隆

(72)【発明者】

【氏名】吉澤 章博

(72)【発明者】

【氏名】三和田 靖彦

(72)【発明者】

【氏名】石井 博行

(72)【発明者】

【氏名】松本 清市

【審査官】 佐藤 仁美

(56)【参考文献】

【文献】 特開２００２−３４０８２１（ＪＰ，Ａ）

【文献】 特開２０１１−１８５７２２（ＪＰ，Ａ）

【文献】 特開２００８−１４５１９６（ＪＰ，Ａ）

【文献】 国際公開第２０１０／０５２９５０（ＷＯ，Ａ１）

【文献】 欧州特許出願公開第００６８２２６８（ＥＰ，Ａ１）

【文献】 佐藤博隆，ブラッグエッジ透過分光法による結晶組織構造情報のイメージング，RADIOISOTOPES，日本，２０１２年，Vol.61 (2012)，449-462

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｎ ２３／００−２３／２２７、

Ｇ０１Ｔ １／００− ７／１２、

Ｇ２１Ｋ １／００− ３／００、 ５／００− ７／００

(57)【特許請求の範囲】

【請求項1】

電池を透過した中性子を検出することで、電池を診断する診断方法であって、
電池を透過した中性子のスペクトルを測定するステップと、
前記電池の電圧と負極材料のブラッグエッジとの相関を示す検量線を参照して、前記中性子のスペクトルから前記電池の電圧を求めて、電池を診断するステップと、を備えた電池の診断方法。

【請求項2】

前記電池の設計データを基に、中性子輸送計算のシミュレーションを行うことで、前記電池で散乱して、検出器に到達する散乱中性子の量を算出し、
前記検出器と前記電池との距離を変えてシミュレーションした時の前記散乱中性子の量の算出結果に基づいて、前記検出器と前記電池との間の最適距離を求め、
前記電池を最適距離に応じた位置に配置して、前記中性子のスペクトルを測定する請求項１に記載の診断方法。

【請求項3】

前記電池を所定の充電状態に設定して、前記中性子のスペクトルを測定することで、前記電池の充電状態の２次元分布を診断する請求項１、又は２に記載の診断方法。

【請求項4】

前記電池の充電状態を変えて、前記中性子のスペクトルを測定することで、前記電池の内部抵抗の２次元分布を診断する請求項１〜３のいずれか１項に記載の診断方法。

【請求項5】

中性子を電池に向けて出射させる中性子発生装置と、
複数の画素を有し、前記電池を透過した中性子を検出する検出器と、
前記画素毎に前記中性子のスペクトルを測定する処理装置と、を備え、
前記電池の電圧と負極材料のブラッグエッジとの相関を示す検量線を参照して、前記中性子のスペクトルから前記電池の電圧を求めて、前記電池を診断する電池の診断装置。

【請求項6】

前記電池の設計データを基に、中性子輸送計算のシミュレーションを行うことで、前記電池で散乱される散乱中性子の量を算出し、
前記検出器と前記電池との距離を変えてシミュレーションした時の前記散乱中性子の量の算出結果に基づいて、前記検出器と前記電池との間の最適距離を求め、
前記電池を最適距離に応じた位置に配置して、前記中性子のスペクトルを測定する請求項５に記載の診断装置。

【請求項7】

前記電池を充放電する充放電装置をさらに備え、
前記充放電装置によって前記電池を所定の充電状態に設定して、前記中性子のスペクトルを測定することで、前記電池の充電状態の２次元分布を診断する請求項５、又は６に記載の診断装置。

【請求項8】

前記電池の充電状態を変えて、前記中性子のスペクトルを測定することで、前記電池の内部抵抗の２次元分布を診断する請求項５〜７のいずれか１項に記載の診断装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、電池の診断装置、及び診断方法に関するものである。

【背景技術】

【0002】

特許文献１には、Ｘ線、ガンマ線、又は中性子を用いた電池特性評価装置が開示されている。この電池特性評価装置では、Ｘ線、ガンマ線、又は中性子を多方向から走査して電池に照射している。そして、電池を透過したＸ線、ガンマ線、又は中性子を検出して、検出結果から活物質分布を求めている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２００４−１８１４６１号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

特許文献１の評価装置では、Ｘ線、ガンマ線、又は中性子の透過量に基づいて、断層像を得ているため、電池の充電状態や劣化状態の情報を得ることができないという問題がある。

【0005】

本発明は、上記の問題点を鑑みてなされたものであり、電池の充電状態や劣化状態等を適切に診断することができる診断装置、及び診断方法を提供することを目的とする。

【課題を解決するための手段】

【0006】

本発明の一態様に係る診断方法は、電池を透過した中性子を検出することで、電池を診断する診断方法であって、電池を透過した中性子のスペクトルを測定するステップと、前記電池の電圧と負極材料のブラッグエッジとの相関を示す検量線を参照して、前記中性子のスペクトルから前記電池の電圧を求めて、電池を診断するステップと、を備えたものである。これにより、電池の充電状態や劣化状態等を適切に診断することができる。
上記の診断方法では、前記電池の設計データを基に、中性子輸送計算のシミュレーションを行うことで、前記電池で散乱して、前記検出器に到達する散乱中性子の量を算出し、前記検出器と前記電池との距離を変えてシミュレーションした時の前記散乱中性子の量の算出結果に基づいて、前記検出器と前記電池との間の最適距離を求め、前記電池を最適距離に応じた位置に配置して、前記中性子のスペクトルを測定するようにしてもよい。これにより、散乱中性子によるノイズの影響を軽減することができ、適切な診断を行うことができる。
上記の診断方法では、前記電池を所定の充電状態に設定して、前記中性子のスペクトルを測定することで、前記電池の充電状態の２次元分布を診断するようにしてもよい。これにより、電池の充電状態の２次元分布を適切に診断することができる。
上記の診断方法では、前記電池の充電状態を変えて、前記中性子のスペクトルを測定することで、前記電池の内部抵抗の２次元分布を診断するようにしてもよい。これにより、電池の内部抵抗の２次元分布を適切に診断することができる。

【0007】

本発明の一態様に係る診断装置は、中性子を電池に向けて出射させる中性子発生装置と、複数の画素を有し、前記電池を透過した中性子を検出する検出器と、前記画素毎に前記中性子のスペクトルを測定する処理装置と、を備え、前記電池の電圧と負極材料のブラッグエッジとの相関を示す検量線を参照して、前記中性子のスペクトルから前記電池の電圧を求めて、前記電池を診断するものである。これにより、電池の充電状態や劣化状態等を適切に診断することができる。
上記の診断装置では、前記電池の設計データを基に、中性子輸送計算のシミュレーションを行うことで、前記電池で散乱される散乱中性子の量を算出し、前記検出器と前記電池との距離を変えてシミュレーションした時の前記散乱中性子の量の算出結果に基づいて、前記検出器と前記電池との間の最適距離を求め、前記電池を最適距離に応じた位置に配置して、前記中性子のスペクトルを測定するようにしてもよい。これにより、散乱中性子によるノイズの影響を軽減することができ、適切な診断を行うことができる。
上記の診断装置では、前記電池を充放電する充放電装置をさらに備え、前記充放電装置によって前記電池を所定の充電状態に設定して、前記中性子のスペクトルを測定することで、前記電池の充電状態の２次元分布を診断するようにしてもよい。これにより、電池の充電状態の２次元分布を適切に診断することができる。
上記の診断装置では、前記電池の充電状態を変えて、前記中性子のスペクトルを測定することで、前記電池の内部抵抗の２次元分布を診断するようにしてもよい。これにより、電池の内部抵抗の２次元分布を適切に診断することができる。

【発明の効果】

【0008】

本発明によれば、電池の充電状態や劣化状態等を適切に診断することができる診断装置、及び診断方法を提供することができる。

【図面の簡単な説明】

【0009】

【図1】本実施形態にかかる診断装置の構成を模式的に示す図である。

【図2】診断装置の検出器の画素構成を模式的に示す図である。

【図3】本実施形態に係る診断装置において、電池と検出器の適切な距離を設定するための方法を示すフローチャートである。

【図4】中性子輸送計算による散乱中性子のイメージ図である。

【図5】距離Ｄ１における、散乱中性子と透過中性子を示す図である。

【図6】距離Ｄ２における、散乱中性子と透過中性子を示す図である。

【図7】最適距離を求める処理の一例を示すフローチャートである。

【図8】中性子源における中性子の放出時間と放出数を示すグラフである。

【図9】本実施形態１に係る診断方法を示すフローチャートである。

【図10】異なる充電状態における中性子スペクトルを示す図である。

【図11】ブラッグエッジと電圧の相関を示すグラフである。

【図12】本実施形態２０に係る診断方法を示すフローチャートである。

【図13】電池の内部抵抗分布が均一な場合における電位分布を示す図である。

【図14】電池の内部抵抗分布にばらつきがある場合における電位分布を示す図である。

【発明を実施するための形態】

【0010】

以下、本発明に係る処理方法、処理装置の実施形態を、図面に基づいて詳細に説明する。但し、本発明が以下の実施形態に限定される訳ではない。また、説明を明確にするため、以下の記載及び図面は、適宜、簡略化されている。各図において同一の符号は実質的に同じ構成を示す。

【0011】

本実施形態にかかる診断装置は、リチウムイオン電池などの二次電池を診断するものである。電池の充電状態や劣化状態を診断することで、電池の高性能化を図ることができる。特に、車に搭載される大型電池では、２次元で電池内部状態を非破壊で把握することが、電池の高性能化に重要である。

【0012】

従来の計測方法としては、Ｘ線回折法、中性子回折法、ＸＡＦＳ（X-rays Absorption Fine Structure）などがある。しかしながら、回折法では１度に狭い範囲からの１点しか計測することができない。また、ＸＡＦＳでは、Ｘ線の透過量が少なく、大型の電池には適用が困難である。

【0013】

そこで、本実施形態では、中性子透過法を用いて、電池を診断している。具体的には、電池を透過した中性子を２次元アレイ検出器で検出している。そして、２次元アレイ検出器の画素毎に中性子のスペクトルを測定している。さらに、電池の負極の構造変化によるブラッグエッジパターンの変化と充電状態を検量線で対応付ける方法を用いている。すなわち、既知の検量線を参照して、ブラッグエッジから、充電状態の２次元分布を推定している。こうすることで、大型の電池であっても、レントゲン撮影のように２次元で電池状態を診断することができる。

【0014】

また、本実施形態では、電池と検出器との間の距離を最適化している。こうすることで、水素等の軽元素を多く含有する電池特有の散乱ノイズを低減することができる。よって、適切に電池の充電状態や劣化状態を診断することができる。

【0015】

さらに、電池１の充放電状態を変えて、透過中性子のスペクトルを測定することも可能である。こうすることで、電池１の内部抵抗の２次元分布を診断することができる。

【0016】

実施の形態１
本実施の形態にかかる診断装置について、図１を用いて説明する。図１は、診断装置の構成を模式的に示す図である。診断装置は、中性子発生装置１１、検出器１２、処理装置１３、及び充放電装置１５を備えている。

【0017】

中性子発生装置１１は、中性子Ｎ１を電池１に向けて出射させる。中性子発生装置１は、例えば、加速器、及び中性子源（ターゲット）を備えている。加速器は、電子や陽子等の荷電粒子を加速する。そして、加速した荷電粒子のパルスをタングステン、タンタル、水銀などの中性子源に照射することで、パルス状の中性子Ｎ１が発生する。中性子源で発生した中性子をコリメートすることで、入射中性子Ｎ１を得ることができる。また、中性子発生装置１１は、中性子を減速するモデレータ等を備えていてもよい。もちろん、中性子を発生するものであれば、中性子発生装置１１の構成は限定されるものではない。例えば、放射性同位体（ＲＩ）や原子炉からの中性子を、チョッパによってパルス状にしてもよい。中性子発生装置１１は、パルス中性子の発生タイミングに応じたトリガ信号を検出器１２に出力している。

【0018】

中性子発生装置１１で発生した中性子Ｎ１は、診断対象である電池１に入射する。中性子Ｎ１の一部は、電池１を透過して、透過中性子Ｎ２となる。電池１を透過した透過中性子Ｎ２は検出器１２に入射する。検出器１２は、入射した中性子数に応じた検出信号を処理装置１３に出力する。図２に示すように、検出器１２は複数の画素１２ａを備えており、各画素１２ａにおいて、透過中性子Ｎ２が検出される。例えば、検出器１２は、アレイ状に配列された複数の画素１２ａを備えた２次元アレイ検出器である。検出器１２は、例えば、２５６ｃｈ（１６×１６）のＬｉ−Ｇｌａｓｓ検出器である。なお、図１に示すように、検出器１２と電池１とは、距離Ｄだけ離間して配置されている。

【0019】

充放電装置１５は、電池１を充放電して、所望の充電状態に設定する。すなわち、電池１を所定の充電状態として、診断を行うことができる。電池１は、例えば、負極が炭素（黒鉛）であり、正極がリチウム金属酸化物であるリチウムイオン二次電池である。正極としては、例えば、Ｌｉ（ＮｉＣｏＭｎ）Ｏ_２等の三元系材料を用いることができる。もちろん、電池１の材料については特に限定されるものではない。

【0020】

処理装置１３は、例えば、パーソナルコンピュータ等の情報処理装置である。処理装置１３は、プロセッサやメモリ等のハードウェアを有しており、メモリ等に格納されたソフトウェアプログラムにしたがって処理を実行する。

【0021】

処理装置１３は、検出器１２からの検出信号に基づいて、スペクトル測定を行う。例えば、処理装置１３は、トリガ信号と検出信号の時間差に基づいて、ＴＯＦ（Time Of Flight）測定を行う。すなわち、処理装置１３は、中性子発生装置１１から検出器１２までの中性子の飛行時間を計測する。そして、処理装置１３は、距離Ｄと飛行時間とによって、透過中性子Ｎ２のエネルギー（波長）を測定する。これにより、中性子のスペクトル測定を行うことができる。処理装置１３は、スペクトル測定を画素１２ａ毎に行う。すなわち、検出器１２の画素毎にスペクトルが測定される。これにより、透過中性子Ｎ２のスペクトルの２次元分布を１度に測定することができる。

【0022】

さらに、処理装置１３は、電池１の電圧とブラッグエッジとの相関を示す検量線を記憶している。これにより、特定の充電状態における電池状態を２次元的に診断することができる。ここでは、負極のブラッグエッジに着目して、診断を行っている。なお、ブラッグエッジとは、吸収断面積が急激に変化する波長である。例えば、ブラッグエッジでは、中性子が回折するため、電池１を透過する中性子が少なくなる。よって、中性子スペクトルから、ブラッグエッジを画素毎に求めることで、画素に対応する位置での電池１の電圧が推定される。電池１の電圧は充電状態を示すため、電池１の充電状態の２次元分布を診断することができる。なお、透過中性子Ｎ２のスペクトルを用いた診断手法については、後述する。

【0023】

ＴＯＦ法によって、透過中性子Ｎ２の波長（エネルギー）を測定している。ＴＯＦ測定において、電池１内において散乱した散乱中性子が検出器１２で検出されると、ノイズになってしまう。例えば、電池１の特定の位置を透過した透過中性子は、検出器１２の１つの画素１２ａで検出される。一方、電池１内で中性子が散乱すると、散乱中性子は、ランダムな方向に散乱される。電池１内の特定位置に中性子が入射した場合に、散乱するか透過するかで検出画素が異なってしまう。すなわち、特定位置で散乱した散乱中性子は、特定位置を透過した透過中性子が検出される画素と異なる画素で検出される。換言すると、特定位置からの散乱中性子は、検出器１２の様々な画素１２ａで検出されるため、中性子の２次元透過像がぼけてしまう。

【0024】

さらに、中性子が電池１内の元素と衝突して散乱すると、検出器１２に到達する角度が変わる。したがって、透過中性子Ｎ２と散乱中性子とが同じ波長であったとしても、飛行時間が変わってしまう。このため、透過中性子Ｎ２のＴＯＦ測定において、散乱中性子がノイズとなってしまう。本実施の形態では、散乱中性子の影響を軽減するために、検出器１２と電池１との間の距離Ｄが最適化されている。すなわち、電池１内の原子で散乱した散乱中性子の影響が小さくなるように、距離Ｄが設定されている。

【0025】

距離Ｄを最適化する方法について、図３を用いて説明する。図３は、距離Ｄの最適化方法を示すフローチャートである。まず、診断する電池１の厚み、形状、及び元素数密度の情報に基づいて、中性子輸送計算の体系を構築する（Ｓ１０１）。すなわち、実際の電池１の設計データ等を用いて、元素の３次元レイアウトを構築する。そして、中性子輸送計算によって、電池１からの散乱中性子量を距離毎に計算する（Ｓ１０２）。すなわち、前記電池の設計データを基に、中性子輸送計算を行って、散乱中性子の量を算出する。電池１から検出器１２までの距離Ｄを変えて、検出器１２に到達する散乱中性子の量を算出する。そして、散乱中性子量がサチュレートする距離を求めて最適距離とする（Ｓ１０３）。なお、中性子輸送計算のシミュレーションでは、中性子発生装置１１と電池１との距離を一定として、検出器１２の位置のみを変えている。

【0026】

中性子輸送計算により散乱中性子の量を求めると図４に示すようになる。図４では、中性子源１１ａから出射した中性子が電池１で散乱する場合において、散乱中性子の強度が濃淡で示されている。電池１内の元素で散乱した散乱中性子は、図４に示すように、任意の方向に散乱する。一般に、単位面積当たりの散乱中性子の量は、散乱位置からの距離の２乗に反比例する。検出器１２に到達する散乱中性子の量は、距離Ｄの２乗に反比例する。一方、単位面積当たりの透過中性子Ｎ２の量は、中性子発生装置１１から検出器１２の距離の２乗に反比例する。中性子発生装置１１と電池１との距離は、検出器１２と電池１との距離に比べて、十分に長くなっている。例えば、中性子発生装置１１と電池１との距離は数ｍ以上であり、検出器１２と電池１との距離Ｄは１ｍ以下である。したがって、透過中性子Ｎ２の量は、距離Ｄを大きくしたとしても、あまり変化しない。

【0027】

中性子輸送計算では、評価位置に検出器１２を配置した場合に、検出器１２に入射する散乱中性子の量を求める。そして、図４の矢印のように評価位置を離していき、散乱中性子の量を計算する。上記のように、電池１から離れ行くほど、評価位置に入射する散乱中性子の量が減少していく。図５のように、距離Ｄ１において、検出器１２に到達する散乱中性子Ｎ３の量を計算する。さらに、図６のように、距離Ｄ１よりも長い距離Ｄ２として、検出器１２に到達する散乱中性子Ｎ３の量を計算する。距離を所定の値ずつ変えていき、検出器１２に到達する散乱中性子Ｎ３の量を計算していく。散乱中性子Ｎ３の量の算出は、中性子粒子輸送計算の体系を基に行われる。

【0028】

図５のように、最適距離よりも短い距離では、電池１で散乱した散乱中性子Ｎ３の多くが検出器１２に入射する。したがって、散乱中性子Ｎ３によるノイズ成分が大きくなる。処理装置１３が散乱中性子によるノイズを誤検出してしまうため、スペクトル測定の精度が劣化する。一方、最適距離よりも長い距離では、検出器１２に入射する散乱中性子Ｎ３の量が低下するが、検出器１２に入射する透過中性子Ｎ２の量が低下する。よって、検出器１２の検出結果にぼけが発生してしまう。最適距離では、散乱中性子Ｎ３によるノイズが低減して、ぼけも最小限となる。このように、本実施形態では、中性子散乱の多い水素を多く含む電池１を中性子透過法で評価するために、最適位置をシミュレーションにより予め算出しておく。

【0029】

なお、散乱中性子の量は電池１のサイズや含有物によって最適位置が異なる。診断対象の電池１の構成が変わるごとに、最適位置を求めることが好ましい。設計データ毎に、中性子の粒子輸送計算を行うことで、電池１のサイズや含有物が異なっても、最適な測定距離での測定が可能となる。

【0030】

このように、電池１の設計データを基に、中性子輸送計算のシミュレーションを行う。こうすることで、電池１内の元素で散乱して、検出器１２に到達する散乱中性子の量を算出することができる。検出器１２と電池１との距離を変えてシミュレーションした時の散乱中性子の量の算出結果に基づいて、検出器１２と電池１との間の最適距離を求める。すなわち、電池１から検出器１２を離していき、散乱中性子の量が飽和する評価位置を最適位置として設定する。そして、電池１を最適距離に応じた位置に配置して、中性子のスペクトルを測定する。

【0031】

さらに、最適距離を求める方法の具体例について、図７を用いて説明する。図７は、最適距離を求める処理を示すフローチャートであり、図３のフローをより具体的に示すものである。

【0032】

まず、電池１の厚さ、形状、及び元素数密度等の情報をインプットしたシミュレーションレイアウトを作成する（Ｓ２０１）。これにより、電池１の構造に応じたシミュレーションレイアウトが構築される。ここでは、電池１に含まれる元素、その元素の配置、及び元素数密度を示す設計データによって、シミュレーションレイアウトが作成される。

【0033】

次に、中性子源１１ａの情報を入力したシミュレーションレイアウトを作成する（Ｓ２０２）。ここでは、発生する中性子のエネルギー分布と放出時間分布等が中性子源１１ａの情報として入力される。中性子発生装置１１で発生する中性子に応じたシミュレーションレイアウトが構築される。このようにして、Ｓ２０１，及びＳ２０２によって、中性子輸送計算に基づくシミュレーションレイアウトが構築される。

【0034】

中性子源１１ａからの中性子の放出時間の分布例を図８に示す。図８には、低エネルギーの中性子と高エネルギーの中性子の放出時間が示されている。図８において、横軸が放出時間、縦軸が放出される中性子数となっている。図８では、高エネルギー中性子のデータを実線で示し、低エネルギー中性子のデータを破線で示している。

【0035】

図８に示すように、低エネルギーの中性子の放出時間のピークは、高エネルギーの中性子の放出時間のピークよりも遅くなっている。高エネルギーの中性子は、放出時間が短くなっている。低エネルギーの中性子は、高エネルギーの中性子に遅れて放出される。

【0036】

次に、中性子源１１ａと評価位置の関係から、６．７オングストローム（０．６７ｎｍ）付近となる到達時間を計算して、その到達時間を評価時間として設定する（Ｓ２０３）。中性子の波長によって、速度が決まるため、処理装置１３が中性子源１１ａから評価位置までの到達時間を導出する。なお、６．７オングストロームは、負極材料である炭素のブラッグエッジである。中性子源１１ａから評価位置までの距離から、波長６．７オングストロームの中性子が、中性子源１１ａから評価位置に到達する到達時間を処理装置１３が算出する。そして、処理装置１３は、この到達時間を評価時間に設定する。なお、評価時間は、評価位置に応じて変化する。したがって、シミュレーションする評価位置毎に、評価時間を算出する。

【0037】

そして、中性子源１１ａから直接評価位置に届いた透過中性子の数をＳとする（Ｓ２０４）。すなわち、電池１内で散乱せずに透過して評価位置に到達した透過中性子の数をＳとする。電池１で散乱して評価位置に届いた散乱中性子の数をｎとする（Ｓ２０５）。

【0038】

モンテカルロシミュレーションを実施して、評価位置を電池１から遠ざけながら、Ｓ、及びｎの値を導出する（Ｓ２０６）。例えば、電池１を構成する元素に対して、波長毎の吸収断面積、散乱断面積等は既知となっている。したがって、電池１の構成、及び入射中性子のスペクトルが既知であれば、中性子の散乱確率を求めることができる。Ｓ２０１、及びＳ２０２で作成したシミュレーションレイアウトを基に、処理装置１３がシミュレーションを実行する。このように、中性子輸送計算の体系に基づいて、評価位置毎に、透過中性子の数Ｓと、散乱中性子の数ｎを算出する。例えば、公知の中性子輸送計算のモンテカルロコードを用いて、Ｓとｎとを算出することができる。さらに、Ｓ２０３で設定した評価時間近辺に、評価位置に到達する中性子を、Ｓ及びｎとして導出する。こうすることで、より適切に最適位置として設定することができる。なお、シミュレーションでは、中性子発生装置１１と電池１との距離を一定としている。

【0039】

そして、Ｓ／ｎ^１／２＝２となる評価位置を最適位置とする（Ｓ２０７）。このようにすることで、透過中性子を散乱中性子のノイズ成分から十分に判別して、スペクトル測定することができる。すなわち、ノイズの影響を低減することができ、適切にスペクトルを測定することができる。もちろん、Ｓ／ｎ^１／２の値は２に限らず、２以上となっていてもよい。評価位置を電池１から遠ざけていき、Ｓ／ｎ^１／２の値が所定の閾値以上となった位置を最適位置とすることができる。

【0040】

次に、図９を参照して、本実施形態にかかる診断方法について説明する。図９は、本実施形態に係る診断方法を示すフローチャートである。まず、上記したように、中性子輸送計算のシミュレーションによって、電池１と検出器１２との間の最適距離を導出する（Ｓ３０１）。これにより、中性子発生装置１１に対する電池１及び検出器１２の配置が決定する。

【0041】

充放電装置１５によって電池１を充放電して、測定したい電圧となるように、電池１の電池電位を設定する（Ｓ３０２）。具体的には、充放電装置１５が電池の充電状態を所望の値とする。負極の回折面付近を測定領域として設定する（Ｓ３０３）。負極が炭素（Ｃ）の場合、黒鉛結晶の（００２）面由来の回折のピーク波長が６．７オングストロームである。黒鉛結晶の（００２）面の面間隔は０．３３５４ｎｍであるため、回折のピーク波長は、約６．７オングストロームとなる。したがって、中性子のスペクトル測定領域を６．７オングストローム付近に設定する。この波長は、炭素のブラッグエッジである。このようにすることで、スペクトル測定する波長範囲を適切に限定することができ、容易に測定することができる。

【0042】

次に、電池１がない状態で、入射中性子Ｎ１を波長毎に測定する（Ｓ３０４）。これにより、入射中性子Ｎ１のスペクトルが測定される。電池１がない状態での中性子の測定結果をＡとする。次に、電池１と検出器１２が最適距離となるように、電池１を設置して、電池１を透過した透過中性子Ｎ２を波長毎に測定する（Ｓ３０５）。これにより、電池１を透過した透過中性子Ｎ２のスペクトルが測定される。電池１を透過した中性子の測定結果をＢとする。もちろん、Ｓ３０４とＳ３０５とで、中性子発生装置１１に対する検出器１２の位置は同じとなっている。

【0043】

ステップＳ３０４、及びＳ３０５では、上記の通り、処理装置１３がＴＯＦ法によって、検出器１２で検出した中性子のスペクトルを測定している。なお、電池１がない状態での中性子スペクトルの測定は、いずれのタイミングで行われていてもよい。例えば、Ｓ３０５の後に、Ｓ３０４が実施されていてもよい。すなわち、Ｓ３０４、及びＳ３０５において、検出器１２と中性子発生装置１１との距離が同じであれば、いずれの測定を先に行ってもよい。ステップＳ３０４、及びＳ３０５のスペクトル測定は、検出器１２の画素１２ａ毎に行われている。

【0044】

次に測定結果Ａを測定結果Ｂで割って、補正する（Ｓ３０６）。処理装置１３は、Ａ／Ｂを算出して、波長毎の透過率とする。これにより、透過率のスペクトルを求めることができる。図１０は、透過率のスペクトル測定結果を示す。図１０において、横軸が中性子の波長を示し、縦軸が透過率を示している。図１０では、充電０％の電池１と、充電１００％の電池１でのスペクトル測定結果を示している。また、図１０は、負極のブラッグエッジ（６．７オングストローム）付近のスペクトルデータであり、検出器１２の１つの画素１２ａでの測定データを示している。

【0045】

図１０に示すように、電池１の充電状態によって、透過率のスペクトルが変化する。ブラッグエッジでは、中性子が回折するため、透過率が低くなる。充電０％のブラッグエッジλ_０と、充電１００％のブラッグエッジλ_１００は異なっている。すなわち、充電状態に応じて、ブラッグエッジが異なる波長となる。これは、電池１が充電されるにつれて、負極を構成する炭素の元素間にリチウムイオンが吸蔵されていくためであると考えられる。炭素元素間にリチウムイオンが吸蔵されていくと、炭素元素の間隔が広がる。電池１の充電状態に応じて面間隔が変わり、ブラッグエッジが変化する。処理装置１３は、スペクトル解析を行い、ブラッグエッジを導出する。

【0046】

例えば、処理装置１３は、６．７オングストローム付近の透過率のスペクトルデータに対してフィッティングを行って、透過率が極小となる波長を求める。透過率が極小となる波長が、ブラッグエッジの測定値となる。処理装置１３は、検出器１２の画素１２ａ毎に、ブラッグエッジ付近の測定領域において、透過率が極小となる波長を求める。

【0047】

そして、スペクトルの解析結果を、検量線と比較して、波長を電池電圧に変換する（Ｓ３０８）。検量線は、電池１の電圧と、ブラッグエッジとの相関を示すデータである。処理装置１３のメモリには、検量線のデータが予め格納されている。例えば、検量線は、図１１に示すようになる。図１１において、横軸が電池１の電圧を示し、縦軸が負極材料由来のブラッグエッジを示している。電池１の電圧が高くなるほど、負極に吸蔵されるリチウムイオンの数が増える。このため、電圧が高いほど、炭素元素間の間隔が広くなり、ブラッグエッジが長くなる。

【0048】

なお、検量線は、実際の測定によって求めることができる。例えば、電池１の基準サンプルを実際に計測する。具体的には、電池１の基準サンプルを作成して、Ｘ線回折や中性子回折を行う。回折ピークの波長がブラッグエッジに対応する。基準サンプルの充電状態を変えて回折ピークを測定することで、電池１の電圧とブラッグエッジを対応付けることができる。よって、検量線を求めることができる。

【0049】

あるいは、電池１の負極のみを取り出して、負極のみの状態でＸ線回折や中性子回折で回折ピークを求める。そして、充電状態を変えて、負極のみの状態における回折ピークを求める。これにより、負極の電位と回折ピークを対応付けることができる。さらに、負極と正極とを組み合わせた状態で、Ｘ線回折や中性子回折で回折ピークを求める。そして、負極と正極の電位を変化させながら、回折ピークを測定する。こうすることで、電池１の電圧とブラッグエッジを対応付けることができる。よって、検量線を求めることができる。

【0050】

処理装置１３は、検量線を参照して、ブラッグエッジの測定値から電池１の電圧を求める。すなわち、処理装置１３は、透過率が極小となる波長をブラッグエッジの測定値として求め、図１１に示すような検量線から、ブラッグエッジの測定値を電池１の電圧に変換する。画素毎に求められたブラッグエッジの測定値が電池１の局所的な電圧を示すことになる。

【0051】

そして、処理装置１３は、電池１の電圧の測定結果をマッピングして、電池状態を診断する（Ｓ３０９）。Ｓ３０７では、ブラッグエッジの測定値が検出器１２の画素毎に求められている。Ｓ３０８において、処理装置１３は、画素毎に、透過率が極小となる波長を電池１の電圧に変換している。処理装置１３は、電池１の電圧を画素毎に求めて、電圧を画素の値として示した２次元マップを作成する。こうすることで、電池１の充電状態を２次元分布で表示することができる。すなわち、電池１の電圧（充電状態）の２次元分布を求めることができる。こうすることで、電池１の充電状態のばらつきや劣化を診断することができる。

【0052】

各画素１２ａでのブラッグエッジの測定値が電池１の局所的な電圧を示すことになる。電圧が高い画素に対応する箇所ほど、より充電されており、ＳＯＣ（State Of Charge）が高くなる。充放電装置１５が電池を所定の充電状態として、スペクトル測定を行うと、電池１の充電状態や劣化状態等の２次元分布を診断することができる。これにより、電池１の内部状態を診断することができる。

【0053】

上記したように、本実施形態では、アレイ状に画素が配列された検出器１２を用いて、透過中性子を検出している。こうすることで、複数の画素での検出を同時に行うことができる。したがって、電池１を透過した透過中性子のスペクトルの２次元分布を簡便に測定することができる。また、検出器１２を電池１と同程度のサイズとすることで、１度の測定で、電池１全体の診断が可能となる。また、透過中性子を検出しているため、大型の電池１の診断を容易に行うことができる。すなわち、Ｘ線回折や中性子回折のように１点ずつ測定しなくでもよいため、大型の電池１を短時間で測定することができる。
本実施の形態に係る診断方法により、電池１を診断することで、電池１の高性能化を図ることができる。例えば、電池１の診断を行い、その診断結果に基づいて、電池１の製造工程の条件出し等を行う。これにより、バラツキの少ない電池１や高性能の電池１を作製することできるようになる。

【0054】

本実施形態では、充放電装置１５によって電池１をある充電状態にして、スペクトル測定を行っている。こうすることで、ある特定の充電状態における電池１の内部状態を面で診断することができる。すなわち、負極の充電状態や劣化状態を２次元で診断することができる。よって、電池１の充電状態や劣化状態等を適切に診断することができる診断することができる。このように電池１を診断していくことで、電池１の高性能化に資することができる。また、非破壊で診断することができるため、診断結果に基づいて、電池１の良否判断を行ってもよい。

【0055】

さらに、本実施の形態では、電池１と検出器１２との距離が最適化されている。すなわち、シミュレーションによって算出された最適距離に電池１及び検出器１２を設置して、スペクトル測定を行っている。こうすることで、散乱中性子によるノイズの影響を低減することができるため、高いＳＮ比での測定が可能になる。よって、適切に電池１を診断することができる。特に、中性子の透過率は、１％以下であり、非常に小さい。このため、電池１と検出器１２との距離を最適化することで、測定に対するノイズの影響を低減することができる。

【0056】

実施の形態２．
次に、実施の形態２に係る診断方法について、図１２を用いて説明する。図１２は、実施の形態２に係る診断方法を示すフローチャートである。本実施形態では、電池１の充電状態を変えて、中性子のスペクトルを測定している。こうすることで、電池１の内部抵抗の２次元分布を診断することができる。なお、本実施の形態にかかる診断装置は実施の形態１と同様であるため、説明を省略する。診断方法の基本的な手順についても、実施の形態１と共通であるため、適宜説明を省略する。例えば、Ｓ４０１、Ｓ４０３〜Ｓ４０８までのステップは、それぞれＳ３０１、Ｓ３０３〜Ｓ３０８と同様であるため、詳細な説明を省略する。

【0057】

まず、実施の形態１のＳ３０１と同様に、電池１と検出器１２の間の最適距離を導出する（Ｓ４０１）。次に、充放電装置１５によって、電池１の充電と放電を繰り返す（Ｓ４０２）。そして、実施の形態１のＳ３０３〜Ｓ３０８と同様に、測定、及び補正等を実施する（Ｓ４０３〜Ｓ４０８）。このように、本実施の形態では、充放電装置１５が電池１を連続して充放電しながら、スペクトル測定を行っている。こうすることで、充電状態を変えて、電池１の電圧を導出することができる。

【0058】

次に、充放電装置１５から入力された入力電流量と、電池１の電極面積から、予測される面積当りの電位を推定する（Ｓ４０９）。例えば、充放電装置１５から入力される電流をモニタする。そして、充放電装置１５からの入力電流量に基づいて、電池１の電位を推定する。そして、変換した電圧値と推定した電位から、場所ごとの電極の内部抵抗の値を推定して、劣化を診断する（Ｓ４１０）。

【0059】

図１３、図１４を用いて、内部抵抗の診断について説明する。図１３は、内部抵抗が均一な電池１を示し、図１４は、内部抵抗が不均一な電池１を示している。図１３、図１４では、電位、及び内部抵抗の高低を濃淡で示している。

【0060】

図１３に示すように、内部抵抗が均一な電池１では、充放電を繰り返し行っても、電位が均一になる。一方、内部抵抗に差がある電池１では、図１４に示すように充放電の際に内部抵抗値に応じて、電位にばらつきが生じる。すなわち、内部抵抗が高い箇所では、充電速度が遅くなるため、充電中の電位が低くなる。一方、内部抵抗が低い箇所では、充電速度が速くなるため、充電中の電位が高くなる。このように、あるタイミングでの電圧と電位から、電池１の内部抵抗の２次元分布を測定することができる。

【0061】

このように、本実施の形態では、充電状態を変えて、中性子のスペクトルを測定している。こうすることで、電池１の内部抵抗の２次元分布を診断することができる。これにより、内部抵抗のばらつきを低減することができ、より高性能の電池１を作製することが可能になる。なお、１つの電池１に対して、実施の形態１、２の診断を実施してもよい。

【0062】

負極の材料は炭素（黒鉛）に限られるものでなく、結晶構造を有し、充電状態の変化に応じて結晶面間隔が変化するものであればよい。

【0063】

なお、本発明は上記実施の形態に限られたものではなく、趣旨を逸脱しない範囲で適宜変更することが可能である。

【符号の説明】

【0064】

１ 電池
１１ 中性子発生装置
１２ 検出器
１３ 処理装置
１５ 充放電装置

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】