(57)【要約】【課題】係数行列Ｋの算出時間を短くすることによって、粒子径分布の算出時間を短縮する。【解決手段】測定対象である粒子に光を照射する光源２３と、前記光の照射により生じる回折／散乱光の光強度を検出する複数の光検出器２４と、前記各光検出器２４から出力された光強度信号を受信し、ベクトルｓがベクトルｑと係数行列Ｋとの積を含む所定式で表されることに基づいて前記粒子の粒子径分布を算出する演算部３とを具備し、前記演算部３が、前記係数行列Ｋの要素のうち少なくとも２つの要素を並列計算して算出するようにした。【選択図】図４
2018036277.tif 000002
【特許請求の範囲】
【請求項１】
測定対象である粒子に光を照射する光源と、前記光の照射により生じる回折／散乱光の光強度を検出する複数の光検出器と、前記各光検出器から出力された光強度信号を受信し、ベクトルｓがベクトルｑと係数行列Ｋとの積を含む所定式で表されることに基づいて前記粒子の粒子径分布を算出する演算部とを具備し、
前記演算部が、前記粒子の粒子径に依存する複数の第１パラメータ及び前記回折／散乱光の拡がり角に依存する複数の第２パラメータを用いて前記係数行列Ｋの要素を算出するものであり、
前記演算部が、前記係数行列Ｋの要素のうち少なくとも２つの要素を並列計算して算出することを特徴とする粒子径分布測定装置。
ここで、ベクトルｓは、各光検出器から出力された光強度信号の値から得られる回折／散乱光の拡がり角度毎の光強度分布を示すベクトル、ベクトルｑは、測定対象である粒子の粒子径分布を示すベクトル、係数行列Ｋは、ベクトルｑをベクトルｓに変換するための行列である。
【請求項２】
前記演算部が、互いに同じ数の要素を並列計算して算出する複数の演算実行部を有し、
第１の演算実行部が、互いに異なる複数の行それぞれに含まれる要素又は互いに異なる複数の列それぞれに含まれる要素を算出し、
第２の演算実行部が、前記互いに異なる複数の行又は前記互いに異なる複数の列それぞれに含まれ、前記第１の実行部により算出された要素とは別の要素を算出することを特徴とする請求項１記載の粒子径分布測定装置。
【請求項３】
前記第１の演算実行部が、１つの行又は１つの列に含まれるそれぞれの要素の算出し、
前記第２の演算実行部が、前記第１の演算実行部により算出された行とは別の１つの行又は前記第１の演算実行部により算出された列とは別の１つの列に含まれるそれぞれの要素の算出することを特徴とする請求項２記載の粒子径分布測定装置。
【請求項４】
測定対象である粒子に光を照射する光源と、前記光の照射により生じる回折／散乱光の光強度を検出する複数の光検出器とを具備する粒子径分布測定装置に搭載されるプログラムであって、
前記各光検出器から出力された光強度信号を受信し、ベクトルｓがベクトルｑと係数行列Ｋとの積を含む所定式で表されることに基づいて前記粒子の粒子径分布を算出する演算部としての機能をコンピュータに発揮させるものであり、
前記演算部が、前記粒子の粒子径に依存する複数の第１パラメータ及び前記回折／散乱光の拡がり角に依存する複数の第２パラメータを用いて前記係数行列Ｋの要素を算出するように構成されるとともに、
前記演算部が、前記係数行列Ｋの要素のうち少なくとも２つの要素を並列計算して算出するように構成されていることを特徴とするプログラム。
ここで、ベクトルｓは、各光検出器から出力された光強度信号の値から得られる回折／散乱光の拡がり角度毎の光強度分布を示すベクトル、ベクトルｑは、測定対象である粒子の粒子径分布を示すベクトル、係数行列Ｋは、ベクトルｑをベクトルｓに変換するための行列である。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、測定対象である粒子に光を照射し、その照射により生じる回折光や散乱光（以下、回折／散乱光とも言う）の光強度分布に基づいて粒子の粒子径分布を算出する粒子径分布測定装置に関するものである。
【背景技術】
【０００２】
この種の粒子径分布測定装置は、特許文献１に示されるように、回折／散乱光を複数の光検出器で検出し、各光検出器から出力された光強度信号の値により得られる光強度分布ベクトルを用いて、下記の式（１）に基づき粒子径分布を算出するようにしている。
【０００３】
【数１】
2018036277.tif 000003
ここで、ｓは、各光検出器から出力された光強度信号の値から得られる回折／散乱光の拡がり角度毎の光強度分布を示すベクトル、ｑは、測定対象である粒子の粒子径分布を示すベクトル、Ｋは、粒子径分布ベクトルを光強度分布ベクトルに変換するための係数行列である。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−１６４５３９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、係数行列Ｋは、粒子の屈折率等の物性や粒子径や光検出器の配置位置によって定まるものであるところ、上記の式（１）に基づき粒子径分布を算出するためには、まず係数行列Ｋを算出しておく必要がある。
【０００６】
しかしながら、従来は、係数行列Ｋの各要素を順に、ＭＩＥ散乱理論等に基づいて１つずつ求めているので、係数行列Ｋの算出にはかなりの時間を要する。
【０００７】
そこで本発明は、かかる課題を鑑みてなされたものであって、係数行列Ｋの算出時間を短くすることによって、粒子径分布の算出時間を短縮することをその主たる目的とするものである。
【課題を解決するための手段】
【０００８】
すなわち、本発明に係る粒子径分布測定装置は、測定対象である粒子に光を照射する光源と、前記光の照射により生じる回折／散乱光の光強度を検出する複数の光検出器と、前記各光検出器から出力された光強度信号を受信し、ベクトルｓがベクトルｑと係数行列Ｋとの積を含む所定式で表されることに基づいて前記粒子の粒子径分布を算出する演算部とを具備し、前記演算部が、前記係数行列Ｋの要素のうち少なくとも２つの要素を並列計算して算出することを特徴とする粒子径分布測定装置。
ここで、ベクトルｓは、各光検出器から出力された光強度信号の値から得られる回折／散乱光の拡がり角度毎の光強度分布を示すベクトル、ベクトルｑは、測定対象である粒子の粒子径分布を示すベクトル、係数行列Ｋは、ベクトルｑをベクトルｓに変換するための行列である。
なお、所定式としては、例えば上述した式（１）や、この式に例えばノイズ等を表す項を加えた式が挙げられる。
【０００９】
このようなものであれば、係数行列Ｋの要素を並列計算して算出しているので、各要素を１つずつ算出する場合に比べて、係数行列Ｋの算出時間を短くすることができ、ひいては粒子径分布を算出するまでの時間を短縮することができる。
【００１０】
具体的な実施態様としては、前記演算部が、前記係数行列Ｋの要素のうち１つの行に含まれる２以上の要素又は１つの列に含まれる２以上の要素を並列計算して算出するものが挙げられる。
【００１１】
また、係数行列Ｋの算出時間をより短くするためには、前記演算部が、前記係数行列Ｋの１つの要素を算出するために用いられる複数のパラメータの値を算出するとともに、これらの値のうち少なくとも１つを記憶して、前記係数行列Ｋの別の要素を算出するために用いることが好ましい。
【００１２】
ここで、係数行列Ｋは、下記の式（２）に表されるように、
【数２】
2018036277.tif 000004
である。なお、ｘは、光検出器の数、ｙは、測定対象となる粒子径範囲の分割数である。
この係数行列Ｋのうち、同一行に属する要素を算出する際には、回折／散乱光の拡がり角に依存するパラメータを共通して用いることができ、同一列に属する要素を算出する際には、粒子の屈折率及び粒子径に依存するパラメータを共通して用いることができる。
【００１３】
そこで、具体的な実施態様としては、前記演算部が、前記係数行列Ｋの要素のうちある行に含まれる要素を算出する際に記憶した複数の前記パラメータの値と、ある列に含まれる要素を算出する際に記憶した複数の前記パラメータの値とを用いて、これらの行及び列が交わる位置の要素を算出するものが挙げられる。
【００１４】
また、本発明に係る粒子径分布測定装置は、測定対象である粒子に光を照射する光源と、前記光の照射により生じる回折／散乱光の光強度を検出する複数の光検出器と、前記各光検出器から出力された光強度信号を受信し、ベクトルｓがベクトルｑと係数行列Ｋとの積を含む所定式で表されることに基づいて前記粒子の粒子径分布を算出する演算部とを具備し、前記演算部が、前記係数行列Ｋの要素のうち１つの要素を算出するために用いられる前記粒子の粒子径に依存する複数の第１パラメータ及び前記回折／散乱光の拡がり角に依存する複数の第２パラメータのうち少なくとも２つのパラメータを並列計算して算出することを特徴とするものである。
【００１５】
このようなものであれば、係数行列Ｋの要素のうち１つの要素を算出するために用いられる複数の第１パラメータ及び複数の第２パラメータのうち少なくとも２つのパラメータを並列計算しているので、従来に比べて係数行列Ｋを算出する時間を短くすることができ、粒子径分布を算出するまでの時間を短縮することができる。
【００１６】
係数行列Ｋの算出時間をより短くするためには、前記演算部が、２種類の前記第１パラメータ及び２種類の前記第２パラメータを並列計算して算出し、下記の式（３）に基づき係数行列Ｋの各要素を算出するものが好ましい。
【数３】
2018036277.tif 000005
ここで、ｋは、係数行列Ｋの要素の値、ｍは、測定対象である粒子の屈折率、αは、測定対象である粒子の粒子径に関連する値、θは、回折／散乱光の拡がり角、ａ及びｂは、粒子の屈折率及び粒子径に依存する第１パラメータ、π及びτは、回折／散乱光の拡がり角に依存する第２パラメータ、Ｎは、演算部がシグマ記号で表される和を演算する際の最終項を示す値である。
【発明の効果】
【００１７】
このように構成した本発明によれば、係数行列Ｋの算出時間を従来に比べて大幅に短くすることが可能となり、ひいては測定全体にかかる時間を短縮することができる。
【図面の簡単な説明】
【００１８】
【図１】本発明の一実施形態における粒子径分布測定装置を示す概略図。
【図２】同実施形態における演算部の機能構成を示す機能ブロック図。
【図３】同実施形態における係数行列算出部の機能構成を示す機能ブロック図。
【図４】同実施形態における係数行列算出部の算出手順を示す図。
【図５】同実施形態における要素算出部が要素を算出する手順を示すフローチャート。
【図６】同実施形態における粒子径分布測定装置の効果を説明する図。
【図７】その他の実施形態における係数行列算出部の算出手順を示す図。
【発明を実施するための形態】
【００１９】
以下に本発明に係る粒子径分布測定装置１の一実施形態について図面を参照して説明する。
【００２０】
本実施形態に係る粒子径分布測定装置１は、粒子に光を照射した際に生じる回折／散乱光の拡がり角に応じる光強度分布が、ＭＩＥ散乱理論から粒子径によって定まることを利用し、前記回折／散乱光を検出することによって粒子径分布を測定するようにしたものである。
【００２１】
粒子径分布測定装置１は、図１に模式的に示すように、装置本体２と演算部３とを備えている。
【００２２】
装置本体２は、粒子を分散させたサンプルを収容するセル２１と、そのセル２１内の粒子にレンズ２２を介してレーザ光を照射する光源２３たるレーザ装置と、レーザ光の照射により生じる回折／散乱光の光強度を拡がり角に応じて検出する複数の光検出器２４（Ａ）、２４（Ｂ）とを備えたものである。
なお、セル２１は、本実施形態では、バッチ式セルを用いているが、循環式セルを用いても構わない。
【００２３】
演算部３は、物理的に言えば、ＣＰＵ、メモリ、入出力インターフェース等を備えた汎用乃至専用のコンピュータであり、各光検出器２４（Ａ）、２４（Ｂ）から出力された光強度信号を受信して式（１）に基づき粒子径分布を算出するものである。
【００２４】
この演算部３は、前記メモリの所定領域に記憶させた所定プログラムにしたがって、ＣＰＵや周辺機器を協働させることにより、図２に示すように、係数行列Ｋを算出する係数行列算出部３１、各光強度信号の値から求められる光強度分布ベクトルｓと前記係数行列Ｋとを用いて式（１）に基づき粒子径分布を算出する粒子径分布算出部３２としての機能を少なくとも発揮するようにしたものである。
【００２５】
本実施形態では、係数行列算出部３１が特徴的であるので、以下に詳述する。
【００２６】
係数行列算出部３１は、測定対象である粒子の物性、粒子径、及び光検出器２４の配置位置に応じた係数行列Ｋを算出するものであり、本実施形態では、オペレータが入力した粒子の例えば屈折率の値を受け付けて、粒子径及び光検出器２４の配置位置に応じた係数行列Ｋの各要素を算出するように構成されている。
【００２７】
より詳細には、この係数行列算出部３１は、図３に示すように、式（３）に基づいて要素を算出する複数の要素算出部４と、前記メモリの所定領域に設定され、前記各要素算出部４の計算過程で算出される２種類の第１パラメータａ、ｂ及び２種類の第２パラメータπ、τの値を記憶する記憶部５としての機能を備えている。
なお、記憶部５は、本実施形態では、キャッシュメモリの所定領域に設定されたものである。
【００２８】
各要素算出部４は、図３に示すように、式（３）に基づく演算を実行する演算実行部４１と、演算実行部４１が演算に用いるパラメータを算出するパラメータ算出部４２としての機能を具備している。
【００２９】
しかして、本実施形態では、係数行列算出部３１が、上述した要素算出部４を係数行列Ｋの行の数（ｘ）と同数有しており、図４に示すように、これらの要素算出部４が、１つの列に含まれる要素を並列して計算し、この並列計算を第１列から第ｙ列まで順次進めていくことで全ての要素を算出するように構成されている。
なお、ここで言う並列計算とは、計算を開始するタイミングや終了するタイミング等が必ずしも一致している必要はなく、各要素算出部４が、多少のタイミングのずれを有して要素を算出している計算状態も含む。
また、これらの要素算出部４は、複数のパソコンにまたがらせて構成されていても良い。
【００３０】
続いて、要素算出部４が、例えばｉ行ｊ列に位置する要素ｋｉｊを算出する手順について、図３及び図５を参照しながら、要素算出部４における各部の動作説明を兼ねて詳述する。ただし、ｉは、１≦ｉ≦ｘの整数とし、ｊは、１≦ｊ≦ｙの整数とする。
なお、下記の説明では、オペレータにより係数行列算出部３１に粒子の屈折率ｍが入力されているものとする。
【００３１】
ここで、要素ｋｉｊは、測定対象である粒子径範囲を分割したうちｊ番目の範囲に属する単位粒子量の粒子によって生じる回折／散乱光が、複数の光検出器２４のうちｉ番目の光検出器２４に入射したときに検出される光強度を示している。したがって、この要素ｋｉｊを求めるうえで、式（３）における粒子径に関連する値α及び回折／散乱光の拡がり角θは定数となり、要素算出部４は、ｎのみを変数とした下記の式（３）’に基づき要素ｋｉｊを算出する。
【数４】
2018036277.tif 000006
【００３２】
ここで、説明の便宜上、
【数５】
2018036277.tif 000007
とおく。また、以下では、ｎ＝ｍの場合を例として説明する。
【００３３】
上述した式（３）’及び（４）において、本実施形態では、１種目の第１パラメータは、複数あり、それらは、それぞれａｊ（１）、ａｊ（２）、・・・、ａｊ（Ｎ）である。２種目の第１パラメータも同様に複数あり、それらは、それぞれｂｊ（１）、ｂｊ（２）、・・・、ｂｊ（Ｎ）である。
上述した式（３）’及び（４）において、本実施形態では、１種目の第２パラメータは、複数あり、それらは、それぞれπｉ（１）、πｉ（２）、・・・、πｉ（Ｎ）である。２種目の第２パラメータも同様に複数あり、それらは、それぞれτｉ（１）、τｉ（２）、・・・、τｉ（Ｎ）である。
【００３４】
始めに、演算実行部４１は、記憶部５にアクセスして（Ｓ１）、ｊ列又はｉ行に含まれる別の要素を算出する要素算出部４がパラメータａｊ（ｍ）、πｉ（ｍ）、ｂｊ（ｍ）、τｉ（ｍ）を既に算出しており、その値を記憶部５が記憶しているか確認する（Ｓ２１〜２４）。
【００３５】
Ｓ２１〜２４において、前述した別の要素を算出する要素算出部４が既に算出しているパラメータに関しては、演算実行部４１がそのパラメータの値を記憶部５から取得する。
【００３６】
Ｓ２１〜２４において、前述した別の要素を算出する要素算出部４がいずれも算出していないパラメータに関しては、演算実行部４１が、そのパラメータを算出するための算出信号をパラメータ算出部４２へ送信する。
【００３７】
算出信号を受信したパラメータ算出部４２は、その算出信号に対応するパラメータの値を算出して演算実行部４１へ送信するとともに、その値を記憶部５にも送信する（Ｓ３１〜３４）。
【００３８】
演算実行部４１は、記憶部５及びパラメータ算出部４２から取得した各パラメータの値を用いてｆ（ｍ）を算出し、その値をメモリの所定領域に記憶させる（Ｓ４）。
【００３９】
続いて、演算実行部４１は、ｆ（１）からｆ（ｍ）までの値を取得するとともに、それらの値を足し合わせ（Ｓ５）、その総和が一定の値に収束しているかを判断する（Ｓ６）。
なお、収束しているかの判断は、例えば、ｆ（ｍ−１）までの総和とｆ（ｍ）までの総和との比率が一定の値以下となるかによって判断する。
【００４０】
Ｓ６において、前記総和が収束していないと判断された場合は、ｎ＝ｍ＋１として（Ｓ７）、再びＳ１へ戻り、Ｓ６において収束されると判断されるまでステッ１からＳ５までを繰り返す。
【００４１】
Ｓ６において、前記総和が収束していると判断された場合は、演算実行部４１は、その収束値を要素ｋｉｊの値として粒子径算出部に出力する（Ｓ８）。
【００４２】
以上のように構成された本実施形態に係る粒子径分布測定装置１によれば、以下のような効果を得ることができる。
【００４３】
複数の要素算出部４が、ある列に含まれている全ての要素を並列計算して算出するとともに、ある要素算出部４が、計算過程で算出したパラメータの値を記憶部５に送信し、別の要素算出部４が、このパラメータの値を記憶部５から取得して要素を算出するために用いるので、従来のように１つずつ要素を算出する場合に比べて、係数行列Ｋの算出時間を大幅に短くすることができる。
【００４４】
以下にこの効果を詳述する。
【００４５】
複数の要素算出部４が、例えば第１列の要素を並列計算する場合、図６に示すように、２種類の第１パラメータａ１（ｎ）及びｂ１（ｎ）は、各要素を算出するために共通して用いられるパラメータである。したがって、複数の要素算出部４のうちいずれか１つが第１パラメータａ１（ｎ）又はｂ１（ｎ）を算出すれば、別の要素算出部４は、記憶部５からこの値を取得して要素の算出に用いることができる。
【００４６】
次に、複数の要素算出部４が、例えば第２列の要素を並列計算する場合、図６に示すように、各要素を算出するために用いられる２種類の第２パラメータπｉ（ｎ）及びτｉ（ｎ）は、第１列に含まれる要素を算出する過程で既に算出されている。したがって、各要素算出部４は、既に算出されている２種類の第２パラメータπｉ（ｎ）及びτｉ（ｎ）の値を記憶部５から取得し、各要素の算出に用いることができる。
このことは、第３列から第ｙ列の要素を算出する場合においても同様である。
【００４７】
したがって、前述したように、本実施形態に係る粒子径分布測定装置１によれば、従来に比べて、各要素を算出する時間が短くなり、係数行列Ｋの算出時間を大幅に短くすることができ、粒子径分布の算出時間を短縮することができる。
【００４８】
なお、本発明は前記実施形態に限定されるものではない。
【００４９】
例えば、前記実施形態では、１つの列に含まれる複数の要素を並列計算するように構成されていたが、図７上段に示すように、１つの行に含まれる複数の要素を並列計算するように構成されていても良い。
また、図７下段に示すように、斜めに位置する複数の要素を並列計算するように構成されていても良い。この場合も、複数の要素を並列計算して算出されたパラメータの値を、別の要素の算出に用いることが可能である。
【００５０】
さらに、各列の計算順序は自由に選択することができ、第ｙ列から第１列まで順次計算するように構成しても良く、このように構成することで、後方又は側方への回折／散乱光に対応する要素であって、Ｓ６における収束に比較的長い時間を要する要素を優先して計算することができる。
【００５１】
また、前記実施形態では、パラメータ算出部４２が、ａ、π、ｂ、τを順次算出するように構成されているが、これらのパラメータのうち少なくとも２種類のパラメータを並列計算して算出するように構成されていても良い。さらに、例えば１種目の第１パラメータａｊ（１）及び１種目の第１パラメータａｊ（２）のように、同じ種類のパラメータを並列計算して算出するように構成されていても良い。
このように構成することで、係数行列Ｋの算出時間をより短くすることができる。
【００５２】
さらに、前記実施形態では、パラメータ算出部４２が、算出したパラメータの値を全て記憶部５に送信していたが、予め定めた項数（例えばｎ＝１００）までに算出した値を記憶部５に送信し、それ以降（例えばｎ＝１０１以降）に算出した値は記憶部５に送信しないように構成しても良い。
【００５３】
粒子径分布算出部３２は、前記実施形態では、式（１）に基づき粒子径分布を算出するものであったが、この式（１）に例えばノイズ等を表す項を加えた式に基づいて粒子径分布を算出するものであっても良い。
【００５４】
また、式（１）のベクトルｓは、各光検出器２４から出力された光強度信号の値から得られる回折／散乱光の光強度分布を示すベクトルであったが、ベクトルｓを前記光強度信号の値そのものを要素とするベクトルとしても良い。
この場合、粒子径分布を示すベクトルｑをベクトルｓに変換すべく、係数行列Ｋの各要素が、式（２）から得られる値を、例えば各光検出器２４からセル２１までの距離等を用いてさらに変換した値であれば良い。
【００５５】
装置本体２に関して言えば、前記実施形態では、光源２３としてレーザ装置を用いていたが、光源２３として例えば無偏向な光を発するＬＥＤ素子を用いても良い。
この場合、要素算出部４は、光強度の垂直な偏向成分と水平な偏向成分とを表している下記の２つの式（５）及び（６）により得られるｋ１及びｋ２の平均値に基づいて各要素を算出するように構成されていれば良い。
【数６】
2018036277.tif 000008
なお、この式（５）は式（２）と同一であり、光源２３からの光が偏向を有する光であっても、無偏向な光であっても、この式は共通して利用することができる。
【００５６】
その他、本発明は前記各実施形態に限られないし、その各部分構成を組み合わせて
も良く、その趣旨を逸脱しない範囲で種々の変形が可能であるのは言うまでもない。
【符号の説明】
【００５７】
１    ・・・粒子径分布測定装置
２    ・・・装置本体
３    ・・・演算部
３１  ・・・係数行列算出部
４    ・・・要素算出部
４１  ・・・演算実行部
４２  ・・・パラメータ算出部
５    ・・・記憶部
【図１】
2018036277.tif 000009
【図２】
2018036277.tif 000010
【図３】
2018036277.tif 000011
【図４】
2018036277.tif 000012
【図５】
2018036277.tif 000013
【図６】
2018036277.tif 000014
【図７】
2018036277.tif 000015