乳化状態の液体の安定性を評価する方法
【課題】測定された時間ごとの濁度の変化を解析することによって、濁度の変化の半減期をもとめ、その半減期の大小によって、乳化状態の液体の安定性を評価することができる。
【解決手段】
本発明は、測定された時間ごとの光の透過度から濁度に算出し、時間ごとの濁度の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から濁度の変化の半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、つまりことによって、その乳化状態の液体の安定性を評価するのである。
【解決手段】
本発明は、測定された時間ごとの光の透過度から濁度に算出し、時間ごとの濁度の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から濁度の変化の半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、つまりことによって、その乳化状態の液体の安定性を評価するのである。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、測定された時間ごとの光の透過度を吸光度に、その吸光度から濁度に換算し、測定された時間（T）ごとの濁度（D）の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、乳化状態の液体の安定性を評価する方法の発明に関するものである。
【背景技術】
【０００２】
乳化状態の液体とは、水と油のように相互に交じり合わないで、分散質（油）が分散媒（水）中に液滴状に分散している液体のことである。したがって、その液体を長時間放置しておくと、最終的に上相と下相の二相に分離する。
従来、乳化状態の液体についての乳化状態の安定性を評価する方法としては、とある測定時間と決めた時間の間に乳化状態の液体が分離して、上相と下相の二相に分離したかどうか、目視で確認していた（特許文献１から１１参照）。それに付け加えて、上相と下相のそれぞれの長さを測定することによって評価されていた（特許文献１２参照）。この方法は多くの特許出願書類に記載されている方法である（特許文献１から１２参照）。
他の同じような方法としては、乳化状態の液体に、遠心分離機で遠心力をかけた後、上相と下相の二相に分離したかどうか、目視で確認していた（特許文献１３から１６参照）。それに付け加えて、上相と下相のそれぞれの長さを測定することによって評価されていた。
これらの方法では、上相と下相の二相に分離の有無だけの判断しかできなかった。その他にも、エマルションの粒径測定装置（特許文献１７から１８参照）と後方散乱変化率（特許文献１９から２１参照）などがある。
また他方では、分光光度計を用いて乳化状態の安定性を評価する方法がよく使われている。通常、乳化状態の液体では、分光光度計の装置から照射させる光が通過できないので、光の透過度が直接測定することができない。したがって、乳化状態の液体を分散媒で、ある一定の割合で希釈して光の透過度が測定されていた。その光の透過度を吸光度［吸光度＝−１×log(透過度)］に換算して、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、さらに、そこから乳化活性係数（乳化活性係数＝２×濁度÷界面活性剤の濃度÷分散質の体積分率）が算出され、その定量的な値である乳化活性係数をもとに乳化安定性が評価されていた（非特許文献１から２参照） 。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１１−４１５１２号公報
【特許文献２】特開２００９−２５６３７７号公報
【特許文献３】特開２００９−０９６９３６号公報
【特許文献４】特開２００８−２８０２５６号公報
【特許文献５】特開２００８−２６３７９０号公報
【特許文献６】特開２００７−２７５０２０号公報
【特許文献７】特開２００７−２０２５３４号公報
【特許文献８】特開２００５−３１４２８８号公報
【特許文献９】特開２００５−１３４２８９号公報
【特許文献１０】特開２００５−１１２８２３号公報
【特許文献１１】特開２００５−０８２５１７号公報
【特許文献１２】特開２００７−０６８４０７号公報
【特許文献１３】特開２００８−０５４５４０号公報
【特許文献１４】特開２００７−３２５５４４号公報
【特許文献１５】特開２００６−２１２０３９号公報
【特許文献１６】特開２００４−２０１６７２号公報
【特許文献１７】特開２０１０−１８９４６１号公報
【特許文献１８】特開２０１０−０４３１４６号公報
【特許文献１９】特開２００６−０９４８６７号公報
【特許文献２０】特開２００５−０２７６７４号公報
【特許文献２１】特開２００２−２３８４７４号公報
【非特許文献】
【０００４】
【非特許文献１】Kevin N. Pearce, John E. Kinsella著 題名「Emulsifying properties of proteins: evaluation of a turbidimetric technique」 学術雑誌名「Journal of Agricultural and Food Chemistry」、26巻 3号、ページ 716-723、The American Chemical Society、1978年
【非特許文献２】水谷令子著 題名「乳化性測定方法の検討」 学術雑誌名「鈴鹿短期大学紀要」 、10巻、 ページ9-17, 鈴鹿短期大学、 1990年
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし本来、乳化状態の液体を長時間放置しておくと、最終的に上相と下相の二相に分離する。乳化状態の液体は時間とともに変化していくのである。それにもかかわらず、乳化活性係数では、時間とともに変化していくことについては、なにも考慮されていないので、乳化状態の液体の安定性を評価する指標としては問題を含んでいる。
【０００６】
本発明は、測定された時間ごとの光の透過度を吸光度に、その吸光度から濁度に算出し、時間ごとの濁度の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から濁度の変化の半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、乳化状態の液体の安定性を評価する方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
発明の第１態様は、測定された時間ごとの光の透過度を測定し、時間ごとの光の透過度を時間ごとの吸光度に換算し、時間ごとの吸光度から時間ごとの濁度に換算し、測定時間（T）ごとの濁度（D）に、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝をあわせることによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価する方法であって、
乳化試料を作成する乳化試料作製手段ステップと、
該乳化試料作製手段ステップで作成された乳化試料をある一定の割合で希釈する乳化試料希釈手段ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料に光を照射する分光ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料を通過した光の透過度を検出する測光ステップと、
該測光ステップで検出された光の透過度と時間を時間ごとに記憶し、さらに、記憶された光の透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する記憶ステップと、
該記憶ステップで記憶された濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
測定時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝で、その算出濁度（D２）式があらかじめ記憶手段に記憶されており、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する算出濁度初期値設定ステップと、
該算出濁度初期値設定ステップで設定された初期値が代入された算出濁度（D２）式に、該記憶ステップで記憶された測定時間（T）を代入して、算出濁度（D２）を算出して、測定時間（T）と算出濁度（D2）を用いて、測定時間（T）を横軸に、算出濁度（Ｄ２）を縦軸にした算出濁度設定グラフを作成する算出濁度設定グラフ化ステップと、
測定グラフ化ステップで求められた測定化グラフと算出濁度設定グラフ化ステップで求められた算出濁度設定グラフの差の２乗が最小になるように算出濁度（D２）式の測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）を変化させて算出する濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、算出濁度初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶ステップに記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式と算出濁度初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、を有する。
【０００８】
本発明によれば、測定時間（T）ごとの濁度（D）に、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝をあわせることによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価することができる。
発明の第２態様は、測定された時間ごとの光の透過度を測定し、時間ごとの光の透過度を時間ごとの吸光度に換算し、時間ごとの吸光度から時間ごとの濁度に換算し、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定し、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示し、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、その乳化状態の液体の安定性を評価する方法であって、
乳化試料を作成する乳化試料作製手段ステップと、
該乳化試料作製手段ステップで作成された乳化試料をある一定の割合で希釈する乳化試料希釈手段ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料に光を照射する分光ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料を通過した光の透過度を検出する測光ステップと、
該測光ステップで検出された光の透過度と時間を時間ごとに記憶し、さらに、記憶された光の透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する記憶ステップと、
該記憶ステップで記憶された濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を初期値として設定する解析グラフ初期値設定ステップと、
該解析グラフ初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化する解析グラフ化ステップと、
該解析グラフ化ステップでグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲の縦軸のｌｎ（D−Ｄｆｉｎａｌ）と横軸の（T−Ｔｉｎｉｔｉａｌ）の右下がりの直線関係から、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、該濁度の変化速度算出ステップで算出された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに重ね書きをし、それに付け加えて、右下がりの直線関係の式つまり一次関数の式と濁度の変化速度算出ステップで設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、を有する。
【０００９】
本発明によれば、測定された時間（T）ごとの濁度（D）に対して、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定し、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示し、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価することができる。
【００１０】
本発明の第３態様は、第1および第２の態様のいずれかだけの態様をとるばかりでなく、第1および第２の態様を兼ね備えかつ実験者つまりユーザーの選択によって第１もしくは第２態様のいずれかを選択することができることも特徴とする
【発明の効果】
【００１１】
本発明によれば、測定された時間（T）ごとの濁度（D）の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価することができる
【図面の簡単な説明】
【００１２】
【図１】本発明の実施形態１における乳化状態の液体の安定性を評価する方法に用いられる乳化状態の液体の乳化状態安定性評価装置のブロック図である。
【図２】本発明の実施形態１における乳化状態の液体の安定性を評価する方法のフローチャートである。
【図３】本発明の実施形態１と実施形態２における測定化グラフを示す図である。
【図４】本発明の実施形態１と実施形態２における解析初期値としては、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定することを示す図である。
【図５】本発明の実施形態１における解析結果を示す図である。
【図６】本発明の実施形態２における乳化状態の液体の安定性を評価する方法に用いられる乳化状態の液体の乳化状態安定性評価装置のブロック図である。
【図７】本発明の実施形態２における乳化状態の液体の安定性を評価する方法のフローチャートである。
【図８】本発明の実施形態２における右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定することを示す図である。
【図９】本発明の実施形態２における解析結果を示す図である。
【発明を実施するための形態】
【００１３】
以下、本発明を実施形態に沿って詳しく説明する。
【実施形態１】
【００１４】
以下、本発明の乳化状態の液体の安定性を評価する方法の第1番目の一実施形態について図面を参照しながら説明をする。図１は、本発明の第1番目の一実施形態における乳化状態の液体の安定性を評価方法に用いられる乳化状態安定性評価装置のブロック図である。
図１に示すように、乳化状態安定性評価装置１は、乳化試料作製ならびに希釈部２、分光部３、コンピュータ部４を有する。
【００１５】
乳化試料作製ならびに希釈部２は、乳化試料作製手段５と乳化試料希釈手段６を有する。乳化試料作製手段５は、乳化試料を作製する手段である。分散質（油）と分散媒（水）と乳化剤の混合液をホモジナイザーや超音波撹拌装置などの撹拌装置で乳化試料が作成される。乳化試料希釈手段６は、乳化試料作製手段５で作成されて乳化試料を乳化希釈液で希釈する手段である。少量の乳化試料を乳化希釈液にゆっくりと混ぜて分散させる。
【００１６】
分光部３は、分光ステップ７、乳化希釈試料容器８、測光ステップ９を有する。分光ステップ７は光源と回折格子からなり、光を発生し設定された波長の光に分光し、乳化希釈試料容器８に照射するものである。石英もしくはプラスチック製の乳化希釈試料容器８中に分光ステップ９より照射された光が、乳化希釈試料の粒子によって散乱される。乳化希釈試料容器８を乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。測光ステップ９は、普通は光電子増倍管が用いられているが、フォトダイオードアレイ検出器でもよい。
【００１７】
分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出され、時間と関連付けて、後述するコンピュータ部４の記憶手段１０に記憶される。
コンピュータ部４は、記憶手段１０、測定グラフ化手段１１、算出濁度初期値設定手段１２、算出濁度設定グラフ化手段１３、濁度の変化速度算出手段１４、濁度変化の半減期算出手段１５、結果表示手段１６を有している。
【００１８】
記憶手段１０は、上述のとおり分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。分光ステップ７から照射された光の強度と乳化希釈試料の粒子に散乱されずに通過した光の強度から透過度（透過度＝乳化希釈試料の粒子に散乱されずに通過した光の強度÷分光ステップ７から照射された光の強度）が算出され時間と関連付けられて記憶し、時間ごとの透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する。またこればかりでなく記憶手段１０は、測定時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの濁度（D）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝があらかじめ記憶手段に記憶されている。
なお、前もって時間とその時間における吸光度、さらには時間とその時間における濁度などを測定しておき、別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込んでもよい。
【００１９】
測定グラフ化手段（時間ごとの測定された濁度のグラフ算出手段）１１は、記憶手段１０で記憶された濁度を縦軸に、記憶手段１０で記憶された時間を横軸にした測定グラフをもとめるものである。
【００２０】
算出濁度初期値設定手段１２は、算出濁度（D２）式があらかじめ記憶手段１０に記憶されており、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する。
【００２１】
算出濁度設定グラフ化手段１３は、算出濁度初期値設定手段１２により初期値が設定され代入された算出濁度（D２）式に、記憶手段１０によって記録された測定時間を代入して、算出濁度（D2）を算出して、その記録された測定時間を横軸に、算出濁度を縦軸として算出濁度設定グラフを求めるものである。これにより、初期値が代入された算出濁度設定グラフを求めることができる。
濁度の変化速度算出ステップ１４は、測定グラフ化手段１１で求められた測定化グラフと算出濁度設定グラフ化手段１３によって求められた算出濁度設定グラフの差の２乗が最小になるように、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変化させる。具体的には、まず測定化グラフ化手段１１で求められた測定化グラフと算出濁度設定グラフ化手段１３で求められた算出濁度設定グラフを比較する。そして、上記算出濁度初期値設定手段１２によって設定された測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を変化させ、測定グラフ化手段１１で求められた測定化グラフと算出濁度設定グラフ化手段１３によって求められた算出濁度設定グラフの差の２乗が最小になるように、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を求める。これにより、算出濁度（D２）式について、測定化グラフと近似するグラフを求めることができ、なおかつ測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を求めることができる。
【００２２】
濁度変化の半減期算出手段１５は、濁度の変化速度算出ステップ１４で算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出して、それを記憶手段１０に記憶する。
【００２３】
結果表示手段１６は、濁度の変化速度算出手段１４で算出された濁度の変化速度（K）と濁度変化の半減期算出手段１５で算出された濁度変化の半減期を記憶手段１０に記憶させ、算出濁度初期値設定手段で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出手段で算出された時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶手段１０に記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式と算出濁度初期値設定ステップ１２で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出手段で算出された濁度変化の半減期のそれぞれを表示する。
【００２４】
次に、本発明の乳化状態の液体の安定性を評価する方法について、図１と図２を参照しながら説明する。図２は本発明の一実施形態である実施形態１における乳化状態の液体の安定性を評価する方法のフローチャートである。
【００２５】
Ｓ１１において、乳化状態安定性評価装置の電源スイッチ（図示略）がONかが判断される。そして、Ｓ１１でYESと判断された場合には、Ｓ１２において装置初期設定画面（図示略）が表示される。ここで、Ｓ１３で照射される光の波長、乳化希釈試料容器８の厚さ（光路長）、何秒ごとに測定するのかの時間間隔、最大何時間まで測定するのかの最大時間をそれぞれ設定する。それぞれの設定値は乳化試料の状態によってことなり、測定するそれぞれの乳化試料によって調整する必要がある。基本的には、乳化希釈試料容器の厚さ（光路長）を長くするか、最大時間を長くすることで解決することが多いということは、発明者の経験から得られた知見である。
実施例として、乳化希釈試料容器８に照射される光の波長を５００ｎｍ、乳化希釈試料容器の厚さ（光路長）は１ｃｍ、何秒ごとに測定するのかの時間間隔は３０秒間隔、最大時間は１７１０秒とした。装置初期設定画面でそれぞれの値を設定すると、Ｓ１３において、設定された光の波長の光が照射される。あるいは、後述するＳ１７の測定開始ボタンがONとさせたときに光を照射してもよい。
【００２６】
S１４において、乳化試料作製操作に入る。乳化試料は、分散質（油）、分散媒（水）、乳化剤を撹拌装置（ホモジナイザー、超音波撹拌装置など）で撹拌することで得られる。この操作は、実験者の実験の目的と用途に応じてさまざまな乳化状態の乳化試料を調製することが想定される。
実施例として、分散質（油）として、コーンオイル、分散媒として水、乳化剤として牛血清アルブミン、ピーエイチ（ｐH）ならびにイオン強度調節剤としてリン酸水素二ナトリウムとリン酸二水素ナトリウムを用いた。水３００ｍLに、リン酸二水素ナトリウム０．０８６ｇとリン酸水素二ナトリウム３．５１３ｇを溶解させると、ピーエイチ（ｐH）が８．０、イオン強度が０．１のリン酸緩衝液３００ｍLが調製される。ここに牛血清アルブミン１．５ｇを溶解させると０．５％牛血清アルブミンリン酸緩衝溶液（ｐH８．０、イオン強度０．１）が３００ｍL調製される。
コーンオイルを３ｍLと０．５％牛血清アルブミンリン酸緩衝溶液２７ｍLの割合で溶液を調製すると分散質の割合が０．１の乳化前液ができる。
もともと乳化前液に溶存している空気を脱気する（乳化前液を密封状態にし、アスピレータで２０分脱気するとよい）、その理由は、自然界における水溶液は、もともと空気が溶存しており、目に見えない気泡が存在しているのである。したがって、この気泡を除去するためにこの操作をするのである。脱気後、撹拌装置（ホモジナイザー）で５分間撹拌して、白濁の乳化試料が得られる。
【００２７】
S１５において、乳化試料希釈操作に入る。S14において調製させた乳化試料をある一定の割合で希釈する操作である。この操作は、実験者の目的と用途に応じてさまざまな乳化希釈試料を調製することが想定される。実施例においては、上述の操作で得られた白濁の乳化試料７５マイクロリットルを採取し、０．１％ラウリル硫酸ナトリウム溶液（ラウリル硫酸ナトリウム０．１５ｇに水１５０ｍLを溶解させたもの）１５ｍLに入れ、手で軽く振って分散させたものを乳化希釈試料とした。
【００２８】
S１６において、乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作に入る。実験者が乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作をするのである。実施例としては、光路長１ｃｍの石英セルに乳化希釈試料をゆっくりと流しいれ、装置に備えつけた。なお、S１４→S１５→S１６の一連の操作を自動化してもよい。
【００２９】
S１６において、乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作を完了したら、S１７の測定開始ボタンをONにする操作に入る。実施例においては、測定開始ボタンをONにする。
なお実施例にかぎらず、ここで、Ｓ１２で行う乳化希釈試料容器８に照射される光の波長、乳化希釈試料容器の厚さ（光路長）、何秒ごとに測定するのかの時間間隔、最大何時間まで測定するのかの最大時間をそれぞれ設定してもよい。
【００３０】
S１８の測定開始ボタンをONと判断されると、S１９において、分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。S２０において、分光ステップ７から照射された光の強度と乳化希釈試料の粒子に散乱されずに通過した光の強度から透過度（透過度＝乳化希釈試料の粒子に散乱されずに通過した光の強度÷分光ステップ７から照射された光の強度）が算出され時間と関連付けられて記憶し、時間ごとの透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する。
なお、前もって時間とその時間における吸光度、さらには時間とその時間における濁度などを測定しておき、別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込んでもよい。
【００３１】
S２０において、S１９において算出し記憶された濁度を縦軸に、時間を横軸に、実験グラフを作成する。実施例としては、濁度を縦軸に、時間を横軸にした図３が得られる。
【００３２】
Ｓ２１において解析初期値設定におこなう。解析初期値としては、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する。実施例としては、図４を参照にしながら設定すると、解析開始時間（Ｔｉｎｉｔｉａｌ）は、曲線がなだらかに減少し始めている開始時間２４０秒を設定し、解析終了時間は、測定最終時間である１７１０秒を設定し、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）は、測定最終時間である１７１０秒の濁度である０．３９を設定し、濁度の変化速度（K）は、反応開始時間と反応開始時間における濁度、そこから３点後である90秒後の時間である３３０秒とその時の濁度からの傾きを求めて、濁度の変化速度（K）の仮の値として設定した。
Ｓ２１において初期値が設定され代入された算出濁度（D２）式に、Ｓ１９において記録された測定時間を代入して、算出濁度（D2）を算出して、その記録された測定時間を横軸に、算出濁度を縦軸として算出濁度設定グラフを求めるものである。これにより、初期値が代入された算出濁度設定グラフを求めることができる。
Ｓ２３においては、Ｓ２０で求められた測定化グラフとＳ２２で求められた算出濁度設定グラフの差の２乗が最小になるように、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変化させる。具体的には、まずＳ２０で求められた測定化グラフとＳ２２で求められた算出濁度設定グラフを比較する。そして、Ｓ２１で設定された測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を変化させ、Ｓ２０で求められた測定化グラフとＳ２２で求められた算出濁度設定グラフの差の２乗が最小になるように、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を求める。これにより、算出濁度（D２）式について、測定化グラフと近似するグラフを求めることができ、なおかつ測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を求めることができる。算出された濁度の変化速度（K）の数を記憶手段１０に記憶する。
【００３３】
Ｓ２４は、Ｓ２３で算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出して、それを記憶手段１０に記憶する。
【００３４】
Ｓ２５は、Ｓ２３で算出された濁度の変化速度（K）とＳ２４で算出された濁度変化の半減期を記憶手段１０に記憶させ、Ｓ２１で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出手段で算出された時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶手段１０に記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式とＳ２１で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）とＳ２３で算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）とＳ２４で算出された濁度変化の半減期のそれぞれを図５に表示する。
以上説明したように、本発明の実施形態１によれば、測定された時間ごとの濁度の変化を、時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの算出濁度（D２）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝にあわせることによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、乳化状態の液体の安定性を評価することができる。
【実施形態２】
【００３５】
以下、本発明の乳化状態の液体の安定性を評価する方法の第2番目の一実施形態について図面を参照しながら説明をする。図６は、本発明の第2番目の一実施形態における乳化状態の液体の安定性を評価方法に用いられる乳化状態安定性評価装置のブロック図である。
図６に示すように、乳化状態安定性評価装置１は、乳化試料作製ならびに希釈部２、分光部３、コンピュータ部４を有する。
【００３６】
乳化試料作製ならびに希釈部２は、乳化試料作製手段５と乳化試料希釈手段６を有する。乳化試料作製手段５は、乳化試料を作製する手段である。分散質（油）と分散媒（水）と乳化剤の混合液をホモジナイザーや超音波撹拌装置などの撹拌装置で乳化試料が作成される。乳化試料希釈手段６は、乳化試料作製手段５で作成されて乳化試料を乳化希釈液で希釈する手段である。少量の乳化試料を乳化希釈液にゆっくりと混ぜて分散させる。
【００３７】
分光部３は、分光ステップ７、乳化希釈試料容器８、測光ステップ９を有する。分光ステップ７は光源と回折格子からなり、光を発生し設定された波長の光に分光し、乳化希釈試料容器８に照射するものである。石英もしくはプラスチック製の乳化希釈試料容器８中に分光ステップ７より照射された光が、乳化希釈試料の粒子によって散乱される。乳化希釈試料容器８を乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。測光ステップ８は、普通は光電子増倍管が用いられているが、フォトダイオードアレイ検出器でもよい。
【００３８】
分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出され、時間と関連付けて、後述するコンピュータ部４の記憶手段１０に記憶される。
コンピュータ部４は、記憶手段１０、測定グラフ化手段１１、解析グラフ初期値設定手段１２、解析グラフ化手段１３、濁度の変化速度算出手段１４、濁度変化の半減期算出手段１５、結果表示手段１６を有している。
【００３９】
記憶手段１０は、上述のとおり分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。分光ステップ７から照射された光の強度と乳化希釈試料の粒子に散乱されずに通過した光の強度から透過度（透過度＝乳化希釈試料の粒子に散乱されずに通過した光の強度÷分光ステップ７から照射された光の強度）が算出され時間と関連付けられて記憶し、時間ごとの透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する。
なお、前もって時間とその時間における吸光度、さらには時間とその時間における濁度などを測定しておき、別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込んでもよい。
【００４０】
測定グラフ化手段（時間ごとの測定された濁度のグラフ算出手段）１１は、記憶手段１０で記憶された濁度を縦軸に、記憶手段１０で記憶された時間を横軸にした測定グラフをもとめるものである。
【００４１】
解析グラフ初期値設定手段１２は、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定する。
【００４２】
解析グラフ化手段１３は、解析グラフ初期値設定手段１２で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、解析グラフを作成する。
濁度の変化速度算出ステップ１４は、解析グラフ化手段１３でグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲の縦軸のｌｎ（D−Ｄｆｉｎａｌ）と横軸の（T−Ｔｉｎｉｔｉａｌ）の右下がりの直線関係から、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度を算出する。
【００４３】
濁度変化の半減期算出手段１５は、濁度の変化速度算出ステップ１４で算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出して、それを記憶手段１０に記憶する。
【００４４】
濁度の変化速度算出手段１４で算出された濁度の変化速度（K）と濁度変化の半減期算出手段１５で算出された濁度変化の半減期を記憶ステップに記憶させ、濁度の変化速度算出ステップで算出された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに重ね書きをし、それに付け加えて、右下がりの直線関係の式つまり一次関数の式と濁度の変化速度算出ステップ１４で設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示する。
【００４５】
次に、本発明の乳化状態の液体の安定性を評価する方法について、図６と図７を参照しながら説明する。図７は本発明の一実施形態である実施形態２における乳化状態の液体の安定性を評価する方法のフローチャートである。
【００４６】
Ｓ１１において、乳化状態安定性評価装置の電源スイッチ（図示略）がONかが判断される。そして、Ｓ１１でYESと判断された場合には、Ｓ１２において装置初期設定画面（図示略）が表示される。ここで、Ｓ１３で照射される光の波長、乳化希釈試料容器８の厚さ（光路長）、何秒ごとに測定するのかの時間間隔、最大何時間まで測定するのかの最大時間をそれぞれ設定する。それぞれの設定値は乳化試料の状態によってことなり、測定するそれぞれの乳化試料によって調整する必要がある。基本的には、乳化希釈試料容器の厚さ（光路長）を長くするか、最大時間を長くすることで解決することが多いということは、発明者の経験から得られた知見である。
実施例として、乳化希釈試料容器８に照射される光の波長を５００ｎｍ、乳化希釈試料容器の厚さ（光路長さ）は１ｃｍ、何秒ごとに測定するのかの時間間隔は３０秒間隔、最大時間は１７１０秒とした。装置初期設定画面でそれぞれの値を設定すると、Ｓ１３において、設定された光の波長の光が照射される。あるいは、後述するＳ１７の測定開始ボタンがONとさせたときに光を照射してもよい。
【００４７】
S１４において、乳化試料作製操作に入る。乳化試料は、分散質（油）、分散媒（水）、乳化剤を撹拌装置（ホモジナイザー、超音波撹拌装置など）で撹拌することで得られる。この操作は、実験者の実験の目的と用途に応じてさまざまな乳化状態の乳化試料を調製することが想定される。
実施例として、分散質（油）として、コーンオイル、分散媒として水、乳化剤として牛血清アルブミン、ピーエイチ（ｐH）ならびにイオン強度調節剤としてリン酸水素二ナトリウムとリン酸二水素ナトリウムを用いた。水３００ｍLに、リン酸二水素ナトリウム０．０８６ｇとリン酸水素二ナトリウム３．５１３ｇを溶解させると、ピーエイチ（ｐH）が８．０、イオン強度が０．１のリン酸緩衝液３００ｍLが調製される。ここに牛血清アルブミン１．５ｇを溶解させると０．５％牛血清アルブミンリン酸緩衝溶液（ｐH８．０、イオン強度０．１）が３００ｍL調製される。
コーンオイルを３ｍLと０．５％牛血清アルブミンリン酸緩衝溶液２７ｍLの割合で溶液を調製すると分散質の割合が０．１の乳化前液ができる。
もともと乳化前液に溶存している空気を脱気する（乳化前液を密封状態にし、アスピレータで２０分脱気するとよい）、その理由は、自然界における水溶液は、もともと空気が溶存しており、目に見えない気泡が存在しているのである。したがって、この気泡を除去するためにこの操作をするのである。脱気後、撹拌装置（ホモジナイザー）で５分間撹拌して、白濁の乳化試料が得られる。
【００４８】
S１５において、乳化試料希釈操作に入る。S14において調製させた乳化試料をある一定の割合で希釈する操作である。この操作は、実験者の目的と用途に応じてさまざまな乳化希釈試料を調製することが想定される。実施例においては、上述の操作で得られた白濁の乳化試料７５マイクロリットルを採取し、０．１％ラウリル硫酸ナトリウム溶液（ラウリル硫酸ナトリウム０．１５ｇに水１５０ｍLを溶解させたもの）１５ｍLに入れ、手で軽く振って分散させたものを乳化希釈試料とした。
【００４９】
S１６において、乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作に入る。実験者が乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作をするのである。実施例としては、光路長１ｃｍの石英セルに乳化希釈試料をゆっくりと流しいれ、装置に備えつけた。なお、S１４→S１５→S１６の一連の操作を自動化してもよい。
【００５０】
S１６において、乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作を完了したら、S１７の測定開始ボタンをONにする操作に入る。実施例においては、測定開始ボタンをONにする。
なお実施例にかぎらず、ここで、Ｓ１２で行う乳化希釈試料容器８に照射される光の波長、乳化希釈試料容器の厚さ（光路長）、何秒ごとに測定するのかの時間間隔、最大何時間まで測定するのかの最大時間をそれぞれ設定してもよい。
【００５１】
S１８の測定開始ボタンをONと判断されると、S１９において、分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。S２０において、分光ステップ７から照射された光の強度と乳化希釈試料の粒子に散乱されずに通過した光の強度から透過度（透過度＝乳化希釈試料の粒子に散乱されずに通過した光の強度÷分光ステップ７から照射された光の強度）が算出され時間と関連付けられて記憶し、時間ごとの透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する。
なお、前もって時間とその時間における吸光度、さらには時間とその時間における濁度などを測定しておき、別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込んでもよい。
【００５２】
S２０において、S１９において算出し記憶された濁度を縦軸に、時間を横軸に、測定グラフを作成する。実施例としては、濁度を縦軸に、時間を横軸にした図３が得られる。
【００５３】
Ｓ２１においては解析初期値設定におこなう。解析初期値としては、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定する。実施例としては、図４を参照にしながら設定すると、解析開始時間（Ｔｉｎｉｔｉａｌ）は、曲線がなだらかに減少し始めている開始時間２４０秒を設定し、解析終了時間は、測定最終時間である１７１０秒を設定し、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）は、測定最終時間である１７１０秒の濁度である０．３９と設定した。
Ｓ２２の解析式設定グラフ化においては、S２１で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、解析設定グラフを作成する。実施例においては、図８が得られる。
Ｓ２３は、Ｓ２２の解析式設定グラフ化においてグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定する。実施例においては、図８のとおり最初から最後から１０点前まで明確な右下がりの傾きをもつ直線関係を示すので最初から最後から１０点前までを選択した。
【００５４】
Ｓ２４においては、S２３で右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定すると、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度を算出して、それぞれを記憶手段１０に記憶する。
Ｓ２５においては、S２４で算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出して、それを記憶手段１０に記憶する。
S２６においては、記憶手段１０に記憶された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、結果の右下がりの直線関係の式つまり一次関数の式とし、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに直線として重ね書きをし、それに付け加えて、結果の右下がりの直線関係の式つまり一次関数の式とS２３で設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と濁度の変化速度（K）とS25で算出された濁度変化の半減期のそれぞれを表示する。実施例をしては、図９のとおり示す。
【００５５】
以上説明したように、本発明の実施形態２によれば、測定された時間（T）ごとの濁度（D）に対して、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定し、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示し、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、その半減期の大小によって、その乳化状態の液体の安定性を評価することができる。
本発明の実施形態１によれば、測定された時間ごとの濁度の変化を、時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの算出濁度（D２）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝にあわせることによって、濁度の変化速度（K）を求めることができる。しかし、発明者の経験則によれば、測定された濁度の変化がきれいな右下がりの曲線になることのではなく、乱れることもあるので、解析開始時間（Ｔｉｎｉｔｉａｌ）から解析終了時間の間における濁度の変化において算出濁度（D２）式の曲線にきちんとのらないことがある。このような場合においては、信頼性のある濁度の変化速度（K）が求められないことがある。しかしながら、本発明の実施形態２によれば、測定された時間（T）ごとの濁度（D）に対して、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定し、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示し、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって濁度の変化速度をもとめる方法では、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示す区間だけの計算で濁度の変化速度をもとめることができるので、実施形態１よりも実施形態２の方が実用的な方法であるといえる。
【００５６】
以上、これまで実施形態１と実施形態２を別々の形態の方法として説明してきたが、両方の形態を兼ね備えて、どちらか一方を選択できるようにしておいてもよい。図１の算出濁度初期値設定手段１２もしくは図６の解析グラフ初期値設定手段１２において、図１の実施形態１のブロック図における算出濁度設定グラフ化手段１３→濁度の変化速度算出手段１４→濁度変化の半減期算出手段１５→記憶手段１０→結果表示手段１６の流れか、もしくは、図６の実施形態２のブロック図における解析グラフ化手段１３→濁度の変化速度算出手段１４→濁度変化の半減期算出手段１５→記憶手段１０→結果表示手段１６の流れのいずれかを選択できるようにしてもよい。同様に、図２もしくは図７のS２１において、図２の実施形態１のフローチャートにおけるS２２→S２３→S２４→S２５の流れか、もしくは、図７の実施形態２のフローチャートにおけるS２２→S２３→S２４→S２５→S２６の流れのいずれかを選択できるようにしてもよい。
【００５７】
以上、本発明によれば、実施形態１と実施形態２のどちらか一方の形態だけでも、実施形態１と実施形態２の両方の形態を兼ね備えて、どちらか一方を選択できるようにしておいても、測定された時間（T）ごとの濁度（D）の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価することができる。
【符号の説明】
【００５８】
１ 乳化状態安定性評価装置
２ 乳化試料作製ならびに希釈部
３ 分光部
４ コンピュータ部
５ 乳化試料作製手段
６ 乳化試料希釈手段
７ 分光ステップ
８ 乳化希釈試料容器
９ 測光ステップ
１０ 記憶手段
１１ 測定グラフ化手段
１２ 算出濁度初期値設定手段
１３ 算出濁度設定グラフ化手段
１４ 濁度の変化速度算出手段
１５ 濁度変化の半減期算出手段
１６ 結果表示手段
１ 乳化状態安定性評価装置
２ 乳化試料作製ならびに希釈部
３ 分光部
４ コンピュータ部
５ 乳化試料作製手段
６ 乳化試料希釈手段
７ 分光ステップ
８ 乳化希釈試料容器
９ 測光ステップ
１０ 記憶手段
１１ 測定グラフ化手段
１２ 解析グラフ初期値設定手段
１３ 解析グラフ化手段
１４ 濁度の変化速度算出手段
１５ 濁度変化の半減期算出手段
１６ 結果表示手段
【技術分野】
【０００１】
本発明は、測定された時間ごとの光の透過度を吸光度に、その吸光度から濁度に換算し、測定された時間（T）ごとの濁度（D）の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、乳化状態の液体の安定性を評価する方法の発明に関するものである。
【背景技術】
【０００２】
乳化状態の液体とは、水と油のように相互に交じり合わないで、分散質（油）が分散媒（水）中に液滴状に分散している液体のことである。したがって、その液体を長時間放置しておくと、最終的に上相と下相の二相に分離する。
従来、乳化状態の液体についての乳化状態の安定性を評価する方法としては、とある測定時間と決めた時間の間に乳化状態の液体が分離して、上相と下相の二相に分離したかどうか、目視で確認していた（特許文献１から１１参照）。それに付け加えて、上相と下相のそれぞれの長さを測定することによって評価されていた（特許文献１２参照）。この方法は多くの特許出願書類に記載されている方法である（特許文献１から１２参照）。
他の同じような方法としては、乳化状態の液体に、遠心分離機で遠心力をかけた後、上相と下相の二相に分離したかどうか、目視で確認していた（特許文献１３から１６参照）。それに付け加えて、上相と下相のそれぞれの長さを測定することによって評価されていた。
これらの方法では、上相と下相の二相に分離の有無だけの判断しかできなかった。その他にも、エマルションの粒径測定装置（特許文献１７から１８参照）と後方散乱変化率（特許文献１９から２１参照）などがある。
また他方では、分光光度計を用いて乳化状態の安定性を評価する方法がよく使われている。通常、乳化状態の液体では、分光光度計の装置から照射させる光が通過できないので、光の透過度が直接測定することができない。したがって、乳化状態の液体を分散媒で、ある一定の割合で希釈して光の透過度が測定されていた。その光の透過度を吸光度［吸光度＝−１×log(透過度)］に換算して、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、さらに、そこから乳化活性係数（乳化活性係数＝２×濁度÷界面活性剤の濃度÷分散質の体積分率）が算出され、その定量的な値である乳化活性係数をもとに乳化安定性が評価されていた（非特許文献１から２参照） 。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１１−４１５１２号公報
【特許文献２】特開２００９−２５６３７７号公報
【特許文献３】特開２００９−０９６９３６号公報
【特許文献４】特開２００８−２８０２５６号公報
【特許文献５】特開２００８−２６３７９０号公報
【特許文献６】特開２００７−２７５０２０号公報
【特許文献７】特開２００７−２０２５３４号公報
【特許文献８】特開２００５−３１４２８８号公報
【特許文献９】特開２００５−１３４２８９号公報
【特許文献１０】特開２００５−１１２８２３号公報
【特許文献１１】特開２００５−０８２５１７号公報
【特許文献１２】特開２００７−０６８４０７号公報
【特許文献１３】特開２００８−０５４５４０号公報
【特許文献１４】特開２００７−３２５５４４号公報
【特許文献１５】特開２００６−２１２０３９号公報
【特許文献１６】特開２００４−２０１６７２号公報
【特許文献１７】特開２０１０−１８９４６１号公報
【特許文献１８】特開２０１０−０４３１４６号公報
【特許文献１９】特開２００６−０９４８６７号公報
【特許文献２０】特開２００５−０２７６７４号公報
【特許文献２１】特開２００２−２３８４７４号公報
【非特許文献】
【０００４】
【非特許文献１】Kevin N. Pearce, John E. Kinsella著 題名「Emulsifying properties of proteins: evaluation of a turbidimetric technique」 学術雑誌名「Journal of Agricultural and Food Chemistry」、26巻 3号、ページ 716-723、The American Chemical Society、1978年
【非特許文献２】水谷令子著 題名「乳化性測定方法の検討」 学術雑誌名「鈴鹿短期大学紀要」 、10巻、 ページ9-17, 鈴鹿短期大学、 1990年
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし本来、乳化状態の液体を長時間放置しておくと、最終的に上相と下相の二相に分離する。乳化状態の液体は時間とともに変化していくのである。それにもかかわらず、乳化活性係数では、時間とともに変化していくことについては、なにも考慮されていないので、乳化状態の液体の安定性を評価する指標としては問題を含んでいる。
【０００６】
本発明は、測定された時間ごとの光の透過度を吸光度に、その吸光度から濁度に算出し、時間ごとの濁度の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から濁度の変化の半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、乳化状態の液体の安定性を評価する方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
発明の第１態様は、測定された時間ごとの光の透過度を測定し、時間ごとの光の透過度を時間ごとの吸光度に換算し、時間ごとの吸光度から時間ごとの濁度に換算し、測定時間（T）ごとの濁度（D）に、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝をあわせることによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価する方法であって、
乳化試料を作成する乳化試料作製手段ステップと、
該乳化試料作製手段ステップで作成された乳化試料をある一定の割合で希釈する乳化試料希釈手段ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料に光を照射する分光ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料を通過した光の透過度を検出する測光ステップと、
該測光ステップで検出された光の透過度と時間を時間ごとに記憶し、さらに、記憶された光の透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する記憶ステップと、
該記憶ステップで記憶された濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
測定時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝で、その算出濁度（D２）式があらかじめ記憶手段に記憶されており、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する算出濁度初期値設定ステップと、
該算出濁度初期値設定ステップで設定された初期値が代入された算出濁度（D２）式に、該記憶ステップで記憶された測定時間（T）を代入して、算出濁度（D２）を算出して、測定時間（T）と算出濁度（D2）を用いて、測定時間（T）を横軸に、算出濁度（Ｄ２）を縦軸にした算出濁度設定グラフを作成する算出濁度設定グラフ化ステップと、
測定グラフ化ステップで求められた測定化グラフと算出濁度設定グラフ化ステップで求められた算出濁度設定グラフの差の２乗が最小になるように算出濁度（D２）式の測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）を変化させて算出する濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、算出濁度初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶ステップに記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式と算出濁度初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、を有する。
【０００８】
本発明によれば、測定時間（T）ごとの濁度（D）に、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝をあわせることによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価することができる。
発明の第２態様は、測定された時間ごとの光の透過度を測定し、時間ごとの光の透過度を時間ごとの吸光度に換算し、時間ごとの吸光度から時間ごとの濁度に換算し、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定し、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示し、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、その乳化状態の液体の安定性を評価する方法であって、
乳化試料を作成する乳化試料作製手段ステップと、
該乳化試料作製手段ステップで作成された乳化試料をある一定の割合で希釈する乳化試料希釈手段ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料に光を照射する分光ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料を通過した光の透過度を検出する測光ステップと、
該測光ステップで検出された光の透過度と時間を時間ごとに記憶し、さらに、記憶された光の透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する記憶ステップと、
該記憶ステップで記憶された濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を初期値として設定する解析グラフ初期値設定ステップと、
該解析グラフ初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化する解析グラフ化ステップと、
該解析グラフ化ステップでグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲の縦軸のｌｎ（D−Ｄｆｉｎａｌ）と横軸の（T−Ｔｉｎｉｔｉａｌ）の右下がりの直線関係から、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、該濁度の変化速度算出ステップで算出された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに重ね書きをし、それに付け加えて、右下がりの直線関係の式つまり一次関数の式と濁度の変化速度算出ステップで設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、を有する。
【０００９】
本発明によれば、測定された時間（T）ごとの濁度（D）に対して、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定し、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示し、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価することができる。
【００１０】
本発明の第３態様は、第1および第２の態様のいずれかだけの態様をとるばかりでなく、第1および第２の態様を兼ね備えかつ実験者つまりユーザーの選択によって第１もしくは第２態様のいずれかを選択することができることも特徴とする
【発明の効果】
【００１１】
本発明によれば、測定された時間（T）ごとの濁度（D）の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価することができる
【図面の簡単な説明】
【００１２】
【図１】本発明の実施形態１における乳化状態の液体の安定性を評価する方法に用いられる乳化状態の液体の乳化状態安定性評価装置のブロック図である。
【図２】本発明の実施形態１における乳化状態の液体の安定性を評価する方法のフローチャートである。
【図３】本発明の実施形態１と実施形態２における測定化グラフを示す図である。
【図４】本発明の実施形態１と実施形態２における解析初期値としては、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定することを示す図である。
【図５】本発明の実施形態１における解析結果を示す図である。
【図６】本発明の実施形態２における乳化状態の液体の安定性を評価する方法に用いられる乳化状態の液体の乳化状態安定性評価装置のブロック図である。
【図７】本発明の実施形態２における乳化状態の液体の安定性を評価する方法のフローチャートである。
【図８】本発明の実施形態２における右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定することを示す図である。
【図９】本発明の実施形態２における解析結果を示す図である。
【発明を実施するための形態】
【００１３】
以下、本発明を実施形態に沿って詳しく説明する。
【実施形態１】
【００１４】
以下、本発明の乳化状態の液体の安定性を評価する方法の第1番目の一実施形態について図面を参照しながら説明をする。図１は、本発明の第1番目の一実施形態における乳化状態の液体の安定性を評価方法に用いられる乳化状態安定性評価装置のブロック図である。
図１に示すように、乳化状態安定性評価装置１は、乳化試料作製ならびに希釈部２、分光部３、コンピュータ部４を有する。
【００１５】
乳化試料作製ならびに希釈部２は、乳化試料作製手段５と乳化試料希釈手段６を有する。乳化試料作製手段５は、乳化試料を作製する手段である。分散質（油）と分散媒（水）と乳化剤の混合液をホモジナイザーや超音波撹拌装置などの撹拌装置で乳化試料が作成される。乳化試料希釈手段６は、乳化試料作製手段５で作成されて乳化試料を乳化希釈液で希釈する手段である。少量の乳化試料を乳化希釈液にゆっくりと混ぜて分散させる。
【００１６】
分光部３は、分光ステップ７、乳化希釈試料容器８、測光ステップ９を有する。分光ステップ７は光源と回折格子からなり、光を発生し設定された波長の光に分光し、乳化希釈試料容器８に照射するものである。石英もしくはプラスチック製の乳化希釈試料容器８中に分光ステップ９より照射された光が、乳化希釈試料の粒子によって散乱される。乳化希釈試料容器８を乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。測光ステップ９は、普通は光電子増倍管が用いられているが、フォトダイオードアレイ検出器でもよい。
【００１７】
分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出され、時間と関連付けて、後述するコンピュータ部４の記憶手段１０に記憶される。
コンピュータ部４は、記憶手段１０、測定グラフ化手段１１、算出濁度初期値設定手段１２、算出濁度設定グラフ化手段１３、濁度の変化速度算出手段１４、濁度変化の半減期算出手段１５、結果表示手段１６を有している。
【００１８】
記憶手段１０は、上述のとおり分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。分光ステップ７から照射された光の強度と乳化希釈試料の粒子に散乱されずに通過した光の強度から透過度（透過度＝乳化希釈試料の粒子に散乱されずに通過した光の強度÷分光ステップ７から照射された光の強度）が算出され時間と関連付けられて記憶し、時間ごとの透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する。またこればかりでなく記憶手段１０は、測定時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの濁度（D）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝があらかじめ記憶手段に記憶されている。
なお、前もって時間とその時間における吸光度、さらには時間とその時間における濁度などを測定しておき、別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込んでもよい。
【００１９】
測定グラフ化手段（時間ごとの測定された濁度のグラフ算出手段）１１は、記憶手段１０で記憶された濁度を縦軸に、記憶手段１０で記憶された時間を横軸にした測定グラフをもとめるものである。
【００２０】
算出濁度初期値設定手段１２は、算出濁度（D２）式があらかじめ記憶手段１０に記憶されており、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する。
【００２１】
算出濁度設定グラフ化手段１３は、算出濁度初期値設定手段１２により初期値が設定され代入された算出濁度（D２）式に、記憶手段１０によって記録された測定時間を代入して、算出濁度（D2）を算出して、その記録された測定時間を横軸に、算出濁度を縦軸として算出濁度設定グラフを求めるものである。これにより、初期値が代入された算出濁度設定グラフを求めることができる。
濁度の変化速度算出ステップ１４は、測定グラフ化手段１１で求められた測定化グラフと算出濁度設定グラフ化手段１３によって求められた算出濁度設定グラフの差の２乗が最小になるように、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変化させる。具体的には、まず測定化グラフ化手段１１で求められた測定化グラフと算出濁度設定グラフ化手段１３で求められた算出濁度設定グラフを比較する。そして、上記算出濁度初期値設定手段１２によって設定された測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を変化させ、測定グラフ化手段１１で求められた測定化グラフと算出濁度設定グラフ化手段１３によって求められた算出濁度設定グラフの差の２乗が最小になるように、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を求める。これにより、算出濁度（D２）式について、測定化グラフと近似するグラフを求めることができ、なおかつ測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を求めることができる。
【００２２】
濁度変化の半減期算出手段１５は、濁度の変化速度算出ステップ１４で算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出して、それを記憶手段１０に記憶する。
【００２３】
結果表示手段１６は、濁度の変化速度算出手段１４で算出された濁度の変化速度（K）と濁度変化の半減期算出手段１５で算出された濁度変化の半減期を記憶手段１０に記憶させ、算出濁度初期値設定手段で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出手段で算出された時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶手段１０に記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式と算出濁度初期値設定ステップ１２で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出手段で算出された濁度変化の半減期のそれぞれを表示する。
【００２４】
次に、本発明の乳化状態の液体の安定性を評価する方法について、図１と図２を参照しながら説明する。図２は本発明の一実施形態である実施形態１における乳化状態の液体の安定性を評価する方法のフローチャートである。
【００２５】
Ｓ１１において、乳化状態安定性評価装置の電源スイッチ（図示略）がONかが判断される。そして、Ｓ１１でYESと判断された場合には、Ｓ１２において装置初期設定画面（図示略）が表示される。ここで、Ｓ１３で照射される光の波長、乳化希釈試料容器８の厚さ（光路長）、何秒ごとに測定するのかの時間間隔、最大何時間まで測定するのかの最大時間をそれぞれ設定する。それぞれの設定値は乳化試料の状態によってことなり、測定するそれぞれの乳化試料によって調整する必要がある。基本的には、乳化希釈試料容器の厚さ（光路長）を長くするか、最大時間を長くすることで解決することが多いということは、発明者の経験から得られた知見である。
実施例として、乳化希釈試料容器８に照射される光の波長を５００ｎｍ、乳化希釈試料容器の厚さ（光路長）は１ｃｍ、何秒ごとに測定するのかの時間間隔は３０秒間隔、最大時間は１７１０秒とした。装置初期設定画面でそれぞれの値を設定すると、Ｓ１３において、設定された光の波長の光が照射される。あるいは、後述するＳ１７の測定開始ボタンがONとさせたときに光を照射してもよい。
【００２６】
S１４において、乳化試料作製操作に入る。乳化試料は、分散質（油）、分散媒（水）、乳化剤を撹拌装置（ホモジナイザー、超音波撹拌装置など）で撹拌することで得られる。この操作は、実験者の実験の目的と用途に応じてさまざまな乳化状態の乳化試料を調製することが想定される。
実施例として、分散質（油）として、コーンオイル、分散媒として水、乳化剤として牛血清アルブミン、ピーエイチ（ｐH）ならびにイオン強度調節剤としてリン酸水素二ナトリウムとリン酸二水素ナトリウムを用いた。水３００ｍLに、リン酸二水素ナトリウム０．０８６ｇとリン酸水素二ナトリウム３．５１３ｇを溶解させると、ピーエイチ（ｐH）が８．０、イオン強度が０．１のリン酸緩衝液３００ｍLが調製される。ここに牛血清アルブミン１．５ｇを溶解させると０．５％牛血清アルブミンリン酸緩衝溶液（ｐH８．０、イオン強度０．１）が３００ｍL調製される。
コーンオイルを３ｍLと０．５％牛血清アルブミンリン酸緩衝溶液２７ｍLの割合で溶液を調製すると分散質の割合が０．１の乳化前液ができる。
もともと乳化前液に溶存している空気を脱気する（乳化前液を密封状態にし、アスピレータで２０分脱気するとよい）、その理由は、自然界における水溶液は、もともと空気が溶存しており、目に見えない気泡が存在しているのである。したがって、この気泡を除去するためにこの操作をするのである。脱気後、撹拌装置（ホモジナイザー）で５分間撹拌して、白濁の乳化試料が得られる。
【００２７】
S１５において、乳化試料希釈操作に入る。S14において調製させた乳化試料をある一定の割合で希釈する操作である。この操作は、実験者の目的と用途に応じてさまざまな乳化希釈試料を調製することが想定される。実施例においては、上述の操作で得られた白濁の乳化試料７５マイクロリットルを採取し、０．１％ラウリル硫酸ナトリウム溶液（ラウリル硫酸ナトリウム０．１５ｇに水１５０ｍLを溶解させたもの）１５ｍLに入れ、手で軽く振って分散させたものを乳化希釈試料とした。
【００２８】
S１６において、乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作に入る。実験者が乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作をするのである。実施例としては、光路長１ｃｍの石英セルに乳化希釈試料をゆっくりと流しいれ、装置に備えつけた。なお、S１４→S１５→S１６の一連の操作を自動化してもよい。
【００２９】
S１６において、乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作を完了したら、S１７の測定開始ボタンをONにする操作に入る。実施例においては、測定開始ボタンをONにする。
なお実施例にかぎらず、ここで、Ｓ１２で行う乳化希釈試料容器８に照射される光の波長、乳化希釈試料容器の厚さ（光路長）、何秒ごとに測定するのかの時間間隔、最大何時間まで測定するのかの最大時間をそれぞれ設定してもよい。
【００３０】
S１８の測定開始ボタンをONと判断されると、S１９において、分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。S２０において、分光ステップ７から照射された光の強度と乳化希釈試料の粒子に散乱されずに通過した光の強度から透過度（透過度＝乳化希釈試料の粒子に散乱されずに通過した光の強度÷分光ステップ７から照射された光の強度）が算出され時間と関連付けられて記憶し、時間ごとの透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する。
なお、前もって時間とその時間における吸光度、さらには時間とその時間における濁度などを測定しておき、別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込んでもよい。
【００３１】
S２０において、S１９において算出し記憶された濁度を縦軸に、時間を横軸に、実験グラフを作成する。実施例としては、濁度を縦軸に、時間を横軸にした図３が得られる。
【００３２】
Ｓ２１において解析初期値設定におこなう。解析初期値としては、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する。実施例としては、図４を参照にしながら設定すると、解析開始時間（Ｔｉｎｉｔｉａｌ）は、曲線がなだらかに減少し始めている開始時間２４０秒を設定し、解析終了時間は、測定最終時間である１７１０秒を設定し、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）は、測定最終時間である１７１０秒の濁度である０．３９を設定し、濁度の変化速度（K）は、反応開始時間と反応開始時間における濁度、そこから３点後である90秒後の時間である３３０秒とその時の濁度からの傾きを求めて、濁度の変化速度（K）の仮の値として設定した。
Ｓ２１において初期値が設定され代入された算出濁度（D２）式に、Ｓ１９において記録された測定時間を代入して、算出濁度（D2）を算出して、その記録された測定時間を横軸に、算出濁度を縦軸として算出濁度設定グラフを求めるものである。これにより、初期値が代入された算出濁度設定グラフを求めることができる。
Ｓ２３においては、Ｓ２０で求められた測定化グラフとＳ２２で求められた算出濁度設定グラフの差の２乗が最小になるように、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変化させる。具体的には、まずＳ２０で求められた測定化グラフとＳ２２で求められた算出濁度設定グラフを比較する。そして、Ｓ２１で設定された測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を変化させ、Ｓ２０で求められた測定化グラフとＳ２２で求められた算出濁度設定グラフの差の２乗が最小になるように、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を求める。これにより、算出濁度（D２）式について、測定化グラフと近似するグラフを求めることができ、なおかつ測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）の数を求めることができる。算出された濁度の変化速度（K）の数を記憶手段１０に記憶する。
【００３３】
Ｓ２４は、Ｓ２３で算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出して、それを記憶手段１０に記憶する。
【００３４】
Ｓ２５は、Ｓ２３で算出された濁度の変化速度（K）とＳ２４で算出された濁度変化の半減期を記憶手段１０に記憶させ、Ｓ２１で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出手段で算出された時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶手段１０に記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式とＳ２１で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）とＳ２３で算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）とＳ２４で算出された濁度変化の半減期のそれぞれを図５に表示する。
以上説明したように、本発明の実施形態１によれば、測定された時間ごとの濁度の変化を、時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの算出濁度（D２）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝にあわせることによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、乳化状態の液体の安定性を評価することができる。
【実施形態２】
【００３５】
以下、本発明の乳化状態の液体の安定性を評価する方法の第2番目の一実施形態について図面を参照しながら説明をする。図６は、本発明の第2番目の一実施形態における乳化状態の液体の安定性を評価方法に用いられる乳化状態安定性評価装置のブロック図である。
図６に示すように、乳化状態安定性評価装置１は、乳化試料作製ならびに希釈部２、分光部３、コンピュータ部４を有する。
【００３６】
乳化試料作製ならびに希釈部２は、乳化試料作製手段５と乳化試料希釈手段６を有する。乳化試料作製手段５は、乳化試料を作製する手段である。分散質（油）と分散媒（水）と乳化剤の混合液をホモジナイザーや超音波撹拌装置などの撹拌装置で乳化試料が作成される。乳化試料希釈手段６は、乳化試料作製手段５で作成されて乳化試料を乳化希釈液で希釈する手段である。少量の乳化試料を乳化希釈液にゆっくりと混ぜて分散させる。
【００３７】
分光部３は、分光ステップ７、乳化希釈試料容器８、測光ステップ９を有する。分光ステップ７は光源と回折格子からなり、光を発生し設定された波長の光に分光し、乳化希釈試料容器８に照射するものである。石英もしくはプラスチック製の乳化希釈試料容器８中に分光ステップ７より照射された光が、乳化希釈試料の粒子によって散乱される。乳化希釈試料容器８を乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。測光ステップ８は、普通は光電子増倍管が用いられているが、フォトダイオードアレイ検出器でもよい。
【００３８】
分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出され、時間と関連付けて、後述するコンピュータ部４の記憶手段１０に記憶される。
コンピュータ部４は、記憶手段１０、測定グラフ化手段１１、解析グラフ初期値設定手段１２、解析グラフ化手段１３、濁度の変化速度算出手段１４、濁度変化の半減期算出手段１５、結果表示手段１６を有している。
【００３９】
記憶手段１０は、上述のとおり分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。分光ステップ７から照射された光の強度と乳化希釈試料の粒子に散乱されずに通過した光の強度から透過度（透過度＝乳化希釈試料の粒子に散乱されずに通過した光の強度÷分光ステップ７から照射された光の強度）が算出され時間と関連付けられて記憶し、時間ごとの透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する。
なお、前もって時間とその時間における吸光度、さらには時間とその時間における濁度などを測定しておき、別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込んでもよい。
【００４０】
測定グラフ化手段（時間ごとの測定された濁度のグラフ算出手段）１１は、記憶手段１０で記憶された濁度を縦軸に、記憶手段１０で記憶された時間を横軸にした測定グラフをもとめるものである。
【００４１】
解析グラフ初期値設定手段１２は、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定する。
【００４２】
解析グラフ化手段１３は、解析グラフ初期値設定手段１２で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、解析グラフを作成する。
濁度の変化速度算出ステップ１４は、解析グラフ化手段１３でグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲の縦軸のｌｎ（D−Ｄｆｉｎａｌ）と横軸の（T−Ｔｉｎｉｔｉａｌ）の右下がりの直線関係から、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度を算出する。
【００４３】
濁度変化の半減期算出手段１５は、濁度の変化速度算出ステップ１４で算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出して、それを記憶手段１０に記憶する。
【００４４】
濁度の変化速度算出手段１４で算出された濁度の変化速度（K）と濁度変化の半減期算出手段１５で算出された濁度変化の半減期を記憶ステップに記憶させ、濁度の変化速度算出ステップで算出された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに重ね書きをし、それに付け加えて、右下がりの直線関係の式つまり一次関数の式と濁度の変化速度算出ステップ１４で設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示する。
【００４５】
次に、本発明の乳化状態の液体の安定性を評価する方法について、図６と図７を参照しながら説明する。図７は本発明の一実施形態である実施形態２における乳化状態の液体の安定性を評価する方法のフローチャートである。
【００４６】
Ｓ１１において、乳化状態安定性評価装置の電源スイッチ（図示略）がONかが判断される。そして、Ｓ１１でYESと判断された場合には、Ｓ１２において装置初期設定画面（図示略）が表示される。ここで、Ｓ１３で照射される光の波長、乳化希釈試料容器８の厚さ（光路長）、何秒ごとに測定するのかの時間間隔、最大何時間まで測定するのかの最大時間をそれぞれ設定する。それぞれの設定値は乳化試料の状態によってことなり、測定するそれぞれの乳化試料によって調整する必要がある。基本的には、乳化希釈試料容器の厚さ（光路長）を長くするか、最大時間を長くすることで解決することが多いということは、発明者の経験から得られた知見である。
実施例として、乳化希釈試料容器８に照射される光の波長を５００ｎｍ、乳化希釈試料容器の厚さ（光路長さ）は１ｃｍ、何秒ごとに測定するのかの時間間隔は３０秒間隔、最大時間は１７１０秒とした。装置初期設定画面でそれぞれの値を設定すると、Ｓ１３において、設定された光の波長の光が照射される。あるいは、後述するＳ１７の測定開始ボタンがONとさせたときに光を照射してもよい。
【００４７】
S１４において、乳化試料作製操作に入る。乳化試料は、分散質（油）、分散媒（水）、乳化剤を撹拌装置（ホモジナイザー、超音波撹拌装置など）で撹拌することで得られる。この操作は、実験者の実験の目的と用途に応じてさまざまな乳化状態の乳化試料を調製することが想定される。
実施例として、分散質（油）として、コーンオイル、分散媒として水、乳化剤として牛血清アルブミン、ピーエイチ（ｐH）ならびにイオン強度調節剤としてリン酸水素二ナトリウムとリン酸二水素ナトリウムを用いた。水３００ｍLに、リン酸二水素ナトリウム０．０８６ｇとリン酸水素二ナトリウム３．５１３ｇを溶解させると、ピーエイチ（ｐH）が８．０、イオン強度が０．１のリン酸緩衝液３００ｍLが調製される。ここに牛血清アルブミン１．５ｇを溶解させると０．５％牛血清アルブミンリン酸緩衝溶液（ｐH８．０、イオン強度０．１）が３００ｍL調製される。
コーンオイルを３ｍLと０．５％牛血清アルブミンリン酸緩衝溶液２７ｍLの割合で溶液を調製すると分散質の割合が０．１の乳化前液ができる。
もともと乳化前液に溶存している空気を脱気する（乳化前液を密封状態にし、アスピレータで２０分脱気するとよい）、その理由は、自然界における水溶液は、もともと空気が溶存しており、目に見えない気泡が存在しているのである。したがって、この気泡を除去するためにこの操作をするのである。脱気後、撹拌装置（ホモジナイザー）で５分間撹拌して、白濁の乳化試料が得られる。
【００４８】
S１５において、乳化試料希釈操作に入る。S14において調製させた乳化試料をある一定の割合で希釈する操作である。この操作は、実験者の目的と用途に応じてさまざまな乳化希釈試料を調製することが想定される。実施例においては、上述の操作で得られた白濁の乳化試料７５マイクロリットルを採取し、０．１％ラウリル硫酸ナトリウム溶液（ラウリル硫酸ナトリウム０．１５ｇに水１５０ｍLを溶解させたもの）１５ｍLに入れ、手で軽く振って分散させたものを乳化希釈試料とした。
【００４９】
S１６において、乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作に入る。実験者が乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作をするのである。実施例としては、光路長１ｃｍの石英セルに乳化希釈試料をゆっくりと流しいれ、装置に備えつけた。なお、S１４→S１５→S１６の一連の操作を自動化してもよい。
【００５０】
S１６において、乳化希釈試料を乳化状態安定性評価装置に備え付ける（セットする）操作を完了したら、S１７の測定開始ボタンをONにする操作に入る。実施例においては、測定開始ボタンをONにする。
なお実施例にかぎらず、ここで、Ｓ１２で行う乳化希釈試料容器８に照射される光の波長、乳化希釈試料容器の厚さ（光路長）、何秒ごとに測定するのかの時間間隔、最大何時間まで測定するのかの最大時間をそれぞれ設定してもよい。
【００５１】
S１８の測定開始ボタンをONと判断されると、S１９において、分光ステップ７から照射された光が、乳化希釈試料容器８中を通過し、乳化希釈試料の粒子に散乱されずに通過した光は時間ごとに測光ステップ９に検出される。S２０において、分光ステップ７から照射された光の強度と乳化希釈試料の粒子に散乱されずに通過した光の強度から透過度（透過度＝乳化希釈試料の粒子に散乱されずに通過した光の強度÷分光ステップ７から照射された光の強度）が算出され時間と関連付けられて記憶し、時間ごとの透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する。
なお、前もって時間とその時間における吸光度、さらには時間とその時間における濁度などを測定しておき、別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込んでもよい。
【００５２】
S２０において、S１９において算出し記憶された濁度を縦軸に、時間を横軸に、測定グラフを作成する。実施例としては、濁度を縦軸に、時間を横軸にした図３が得られる。
【００５３】
Ｓ２１においては解析初期値設定におこなう。解析初期値としては、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定する。実施例としては、図４を参照にしながら設定すると、解析開始時間（Ｔｉｎｉｔｉａｌ）は、曲線がなだらかに減少し始めている開始時間２４０秒を設定し、解析終了時間は、測定最終時間である１７１０秒を設定し、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）は、測定最終時間である１７１０秒の濁度である０．３９と設定した。
Ｓ２２の解析式設定グラフ化においては、S２１で設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、解析設定グラフを作成する。実施例においては、図８が得られる。
Ｓ２３は、Ｓ２２の解析式設定グラフ化においてグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定する。実施例においては、図８のとおり最初から最後から１０点前まで明確な右下がりの傾きをもつ直線関係を示すので最初から最後から１０点前までを選択した。
【００５４】
Ｓ２４においては、S２３で右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定すると、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度を算出して、それぞれを記憶手段１０に記憶する。
Ｓ２５においては、S２４で算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出して、それを記憶手段１０に記憶する。
S２６においては、記憶手段１０に記憶された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、結果の右下がりの直線関係の式つまり一次関数の式とし、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに直線として重ね書きをし、それに付け加えて、結果の右下がりの直線関係の式つまり一次関数の式とS２３で設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と濁度の変化速度（K）とS25で算出された濁度変化の半減期のそれぞれを表示する。実施例をしては、図９のとおり示す。
【００５５】
以上説明したように、本発明の実施形態２によれば、測定された時間（T）ごとの濁度（D）に対して、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定し、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示し、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、その半減期の大小によって、その乳化状態の液体の安定性を評価することができる。
本発明の実施形態１によれば、測定された時間ごとの濁度の変化を、時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの算出濁度（D２）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝にあわせることによって、濁度の変化速度（K）を求めることができる。しかし、発明者の経験則によれば、測定された濁度の変化がきれいな右下がりの曲線になることのではなく、乱れることもあるので、解析開始時間（Ｔｉｎｉｔｉａｌ）から解析終了時間の間における濁度の変化において算出濁度（D２）式の曲線にきちんとのらないことがある。このような場合においては、信頼性のある濁度の変化速度（K）が求められないことがある。しかしながら、本発明の実施形態２によれば、測定された時間（T）ごとの濁度（D）に対して、解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を設定し、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示し、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって濁度の変化速度をもとめる方法では、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化すると、右下がりの傾きをもつ直線関係を示す区間だけの計算で濁度の変化速度をもとめることができるので、実施形態１よりも実施形態２の方が実用的な方法であるといえる。
【００５６】
以上、これまで実施形態１と実施形態２を別々の形態の方法として説明してきたが、両方の形態を兼ね備えて、どちらか一方を選択できるようにしておいてもよい。図１の算出濁度初期値設定手段１２もしくは図６の解析グラフ初期値設定手段１２において、図１の実施形態１のブロック図における算出濁度設定グラフ化手段１３→濁度の変化速度算出手段１４→濁度変化の半減期算出手段１５→記憶手段１０→結果表示手段１６の流れか、もしくは、図６の実施形態２のブロック図における解析グラフ化手段１３→濁度の変化速度算出手段１４→濁度変化の半減期算出手段１５→記憶手段１０→結果表示手段１６の流れのいずれかを選択できるようにしてもよい。同様に、図２もしくは図７のS２１において、図２の実施形態１のフローチャートにおけるS２２→S２３→S２４→S２５の流れか、もしくは、図７の実施形態２のフローチャートにおけるS２２→S２３→S２４→S２５→S２６の流れのいずれかを選択できるようにしてもよい。
【００５７】
以上、本発明によれば、実施形態１と実施形態２のどちらか一方の形態だけでも、実施形態１と実施形態２の両方の形態を兼ね備えて、どちらか一方を選択できるようにしておいても、測定された時間（T）ごとの濁度（D）の変化を解析することによって、濁度の変化速度をもとめ、その濁度の変化速度から半減期を算出し、半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小から、乳化状態の液体の安定性を評価することができる。
【符号の説明】
【００５８】
１ 乳化状態安定性評価装置
２ 乳化試料作製ならびに希釈部
３ 分光部
４ コンピュータ部
５ 乳化試料作製手段
６ 乳化試料希釈手段
７ 分光ステップ
８ 乳化希釈試料容器
９ 測光ステップ
１０ 記憶手段
１１ 測定グラフ化手段
１２ 算出濁度初期値設定手段
１３ 算出濁度設定グラフ化手段
１４ 濁度の変化速度算出手段
１５ 濁度変化の半減期算出手段
１６ 結果表示手段
１ 乳化状態安定性評価装置
２ 乳化試料作製ならびに希釈部
３ 分光部
４ コンピュータ部
５ 乳化試料作製手段
６ 乳化試料希釈手段
７ 分光ステップ
８ 乳化希釈試料容器
９ 測光ステップ
１０ 記憶手段
１１ 測定グラフ化手段
１２ 解析グラフ初期値設定手段
１３ 解析グラフ化手段
１４ 濁度の変化速度算出手段
１５ 濁度変化の半減期算出手段
１６ 結果表示手段
【特許請求の範囲】
【請求項１】
乳化状態の液体の安定性を評価する方法であって、
乳化試料を作成する乳化試料作製手段ステップと、
該乳化試料作製手段ステップで作成された乳化試料をある一定の割合で希釈する乳化試料希釈手段ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料に光を照射する分光ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料を通過した光の透過度を検出する測光ステップと、
該測光ステップで検出された光の透過度と時間を時間ごとに記憶し、さらに、記憶された光の透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する記憶ステップと、
該記憶ステップで記憶された濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
測定時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの濁度（D）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝で、その算出濁度（D２）式があらかじめ記憶手段に記憶されており、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する算出濁度（D２）式初期値設定ステップと、
該算出濁度初期値設定ステップで設定された初期値が代入された算出濁度（D２）式に、該記憶ステップで記憶された測定時間（T）を代入して、算出濁度（D２）を算出して、測定時間（T）と算出濁度（D2）を用いて、測定時間（T）を横軸に、算出濁度（Ｄ２）を縦軸にした算出濁度設定グラフを作成する算出濁度設定グラフ化ステップと、
測定グラフ化ステップで求められた測定化グラフと算出濁度設定グラフ化ステップで求められた算出濁度グラフの差の２乗が最小になるように算出濁度（D２）式の測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を変化させて算出する濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、算出濁度初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶ステップに記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式と算出濁度（D２）式初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、
を有する、
半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、乳化状態の液体の安定性を評価する方法。
【請求項２】
請求項１記載の方法を具現化し実行できる装置。
【請求項３】
乳化状態の液体の安定性を評価するプログラムであって、
別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込む、データ読み込みステップと、
該読み込みステップで読み込んだ濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
測定時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの濁度（D）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝で、その算出濁度（D２）式があらかじめ記憶手段に記憶されており、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する算出濁度（D２）式初期値設定ステップと、
該算出濁度（D２）式初期値設定ステップで設定された初期値が代入された算出濁度（D２）式に、該記憶ステップで記憶された測定時間（T）を代入して、算出濁度（D２）を算出して、測定時間（T）と算出濁度（D2）を用いて、測定時間（T）を横軸に、算出濁度（Ｄ２）を縦軸にした算出濁度設定グラフを作成する算出濁度設定グラフ化ステップと、
測定グラフ化ステップで求められた測定化グラフと算出濁度設定グラフ化ステップで求められた算出濁度（D２）式グラフの差の２乗が最小になるように算出濁度（D２）式の測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を変化させて算出する濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、算出濁度（D２）式初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶ステップに記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式と算出濁度（D２）式初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、
を有し、
半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、その乳化状態の液体の安定性を評価する方法を具現化し実行できるプログラム。
【請求項４】
乳化状態の液体の安定性を評価する方法であって、
乳化試料を作成する乳化試料作製手段ステップと、
該乳化試料作製手段ステップで作成された乳化試料をある一定の割合で希釈する乳化試料希釈手段ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料に光を照射する分光ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料を通過した光の透過度を検出する測光ステップと、
該測光ステップで検出された光の透過度と時間を時間ごとに記憶し、さらに、記憶された光の透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する記憶ステップと、
該記憶ステップで記憶された濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を初期値として設定する解析グラフ初期値設定ステップと、
該解析グラフ初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化する解析グラフ化ステップと、
該解析グラフ化ステップでグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲の縦軸のｌｎ（D−Ｄｆｉｎａｌ）と横軸の（T−Ｔｉｎｉｔｉａｌ）の右下がりの直線関係から、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、該濁度の変化速度算出ステップで算出された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに重ね書きをし、それに付け加えて、右下がりの直線関係の式つまり一次関数の式と濁度の変化速度算出ステップで設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、
を有し、
半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、その乳化状態の液体の安定性を評価する方法。
【請求項５】
請求項４記載の方法を具現化し実行できる装置。
【請求項６】
乳化状態の液体の安定性を評価するプログラムであって、
別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込む、データ読み込みステップと、
該読み込みステップで読み込んだ濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を初期値として設定する解析グラフ初期値設定ステップと、
該解析グラフ初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化する解析グラフ化ステップと、
該解析グラフ化ステップでグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲の縦軸のｌｎ（D−Ｄｆｉｎａｌ）と横軸の（T−Ｔｉｎｉｔｉａｌ）の右下がりの直線関係から、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、該濁度の変化速度算出ステップで算出された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに重ね書きをし、それに付け加えて、右下がりの直線関係の式つまり一次関数の式と濁度の変化速度算出ステップで設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、
を有し、
半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、その乳化状態の液体の安定性を評価する方法を具現化し実行するプログラム。
【請求項６】
請求項１記載の方法と請求項４記載の方法の両方を兼ね備え、実験者あるいは使用者の選択により請求項１記載の方法と請求項４記載の方法のどちらか一方の方法を選択し実行することによって乳化状態の液体の安定性を評価する方法
【請求項７】
請求項６記載の方法を具現化し実行できる装置
【請求項８】
請求項３記載のプログラムと請求項６記載のプログラムの両方のプログラムを兼ね備え、実験者あるいは使用者の選択により請求項３記載のプログラムと請求項６記載のプログラムのどちらか一方のプログラムを選択し実行するプログラム。
【請求項１】
乳化状態の液体の安定性を評価する方法であって、
乳化試料を作成する乳化試料作製手段ステップと、
該乳化試料作製手段ステップで作成された乳化試料をある一定の割合で希釈する乳化試料希釈手段ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料に光を照射する分光ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料を通過した光の透過度を検出する測光ステップと、
該測光ステップで検出された光の透過度と時間を時間ごとに記憶し、さらに、記憶された光の透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する記憶ステップと、
該記憶ステップで記憶された濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
測定時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの濁度（D）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝で、その算出濁度（D２）式があらかじめ記憶手段に記憶されており、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する算出濁度（D２）式初期値設定ステップと、
該算出濁度初期値設定ステップで設定された初期値が代入された算出濁度（D２）式に、該記憶ステップで記憶された測定時間（T）を代入して、算出濁度（D２）を算出して、測定時間（T）と算出濁度（D2）を用いて、測定時間（T）を横軸に、算出濁度（Ｄ２）を縦軸にした算出濁度設定グラフを作成する算出濁度設定グラフ化ステップと、
測定グラフ化ステップで求められた測定化グラフと算出濁度設定グラフ化ステップで求められた算出濁度グラフの差の２乗が最小になるように算出濁度（D２）式の測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を変化させて算出する濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、算出濁度初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶ステップに記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式と算出濁度（D２）式初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、
を有する、
半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、乳化状態の液体の安定性を評価する方法。
【請求項２】
請求項１記載の方法を具現化し実行できる装置。
【請求項３】
乳化状態の液体の安定性を評価するプログラムであって、
別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込む、データ読み込みステップと、
該読み込みステップで読み込んだ濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
測定時間（T）、解析開始時間（Ｔｉｎｉｔｉａｌ）、測定時間ごとの濁度（D）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）、濁度の変化速度（K）をそれぞれ変数として表した算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝で、その算出濁度（D２）式があらかじめ記憶手段に記憶されており、測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を未知数として、解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）を初期値として設定し、それに付け加えて測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）の仮の値を初期値として設定する算出濁度（D２）式初期値設定ステップと、
該算出濁度（D２）式初期値設定ステップで設定された初期値が代入された算出濁度（D２）式に、該記憶ステップで記憶された測定時間（T）を代入して、算出濁度（D２）を算出して、測定時間（T）と算出濁度（D2）を用いて、測定時間（T）を横軸に、算出濁度（Ｄ２）を縦軸にした算出濁度設定グラフを作成する算出濁度設定グラフ化ステップと、
測定グラフ化ステップで求められた測定化グラフと算出濁度設定グラフ化ステップで求められた算出濁度（D２）式グラフの差の２乗が最小になるように算出濁度（D２）式の測定時間の最大の範囲における濁度が最小の値（Ｄｆｉｎａｌ）と濁度の変化速度（K）を変化させて算出する濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、算出濁度（D２）式初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析終了時間と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）、そして測定時間（T）を算出濁度（D２）式｛D２＝（Dｉｎｉｔｉａｌ−Ｄｆｉｎａｌ）×exp［（−１）×K×（T−Ｔｉｎｉｔｉａｌ）］＋Ｄｆｉｎａｌ、ここで、exp［ ］はexp関数を表し、自然対数の底つまりネイピア数（ｅ）の（−１）×K×（T−Ｔｉｎｉｔｉａｌ）乗を意味する。｝に代入し、測定時間ごとの算出濁度（D２）を算出し記憶ステップに記憶させ、記憶された算出濁度（D2）を縦軸に、測定時間を横軸にした算出濁度設定グラフと測定グラフを重ね書きし、それに付け加えて、算出濁度（D２）式と算出濁度（D２）式初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）と解析開始時間における濁度（Dｉｎｉｔｉａｌ）と濁度の変化速度算出ステップで算出された測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、
を有し、
半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、その乳化状態の液体の安定性を評価する方法を具現化し実行できるプログラム。
【請求項４】
乳化状態の液体の安定性を評価する方法であって、
乳化試料を作成する乳化試料作製手段ステップと、
該乳化試料作製手段ステップで作成された乳化試料をある一定の割合で希釈する乳化試料希釈手段ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料に光を照射する分光ステップと、
該乳化試料希釈手段ステップで希釈された乳化希釈試料を通過した光の透過度を検出する測光ステップと、
該測光ステップで検出された光の透過度と時間を時間ごとに記憶し、さらに、記憶された光の透過度から時間ごとの吸光度［吸光度＝−１×log(透過度)、ここで、logは常用対数をあらわし、透過度の常用対数に−１をかけたものが吸光度である］を算出し、さらに該吸光度から濁度（濁度＝２．３０３×吸光度÷液相の長さ）に算出して、透過度と吸光度と濁度と時間のそれぞれを時間ごとに記憶する記憶ステップと、
該記憶ステップで記憶された濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を初期値として設定する解析グラフ初期値設定ステップと、
該解析グラフ初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化する解析グラフ化ステップと、
該解析グラフ化ステップでグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲の縦軸のｌｎ（D−Ｄｆｉｎａｌ）と横軸の（T−Ｔｉｎｉｔｉａｌ）の右下がりの直線関係から、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、該濁度の変化速度算出ステップで算出された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに重ね書きをし、それに付け加えて、右下がりの直線関係の式つまり一次関数の式と濁度の変化速度算出ステップで設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、
を有し、
半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、その乳化状態の液体の安定性を評価する方法。
【請求項５】
請求項４記載の方法を具現化し実行できる装置。
【請求項６】
乳化状態の液体の安定性を評価するプログラムであって、
別の記憶媒体（CD、CD-ROM、USBメモリなど）に時間とその時間における吸光度、さらには時間とその時間における濁度などのデータが保存されており、別の記憶媒体に保存された時間とその時間における吸光度を読み込み、時間とその時間における濁度に換算する、もしくは時間とその時間における濁度のデータを読み込む、データ読み込みステップと、
該読み込みステップで読み込んだ濁度を縦軸に、時間を横軸にしたグラフを作成する測定グラフ化ステップと、
解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）を初期値として設定する解析グラフ初期値設定ステップと、
該解析グラフ初期値設定ステップで設定された解析開始時間（Ｔｉｎｉｔｉａｌ）、解析開始時間における濁度（Dｉｎｉｔｉａｌ）、測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）をもとに、測定された時間（T）から解析開始時間（Ｔｉｎｉｔｉａｌ）の差（T−Ｔｉｎｉｔｉａｌ）を、濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）をそれぞれ算出し、さらにその濁度（D）から測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）の差（D−Ｄｆｉｎａｌ）の自然対数つまりｌｎ（D−Ｄｆｉｎａｌ）をそれぞれ算出し記憶ステップに記憶させ、縦軸にｌｎ（D−Ｄｆｉｎａｌ）を、横軸に（T−Ｔｉｎｉｔｉａｌ）として、グラフ化する解析グラフ化ステップと、
該解析グラフ化ステップでグラフ化されたグラフに対して、右下がりの傾きをもつ直線関係になっている範囲つまり解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）を設定して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲の縦軸のｌｎ（D−Ｄｆｉｎａｌ）と横軸の（T−Ｔｉｎｉｔｉａｌ）の右下がりの直線関係から、線形の最小二乗法もちいて、右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B、ここで−１×Aは傾きでBは切片を表す）の傾き（−１×A）と切片Bをもとめ、傾き（−１×A）に−１をかけることによって、濁度の変化速度（K）をもとめる濁度の変化速度算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）から濁度変化の半減期（半減期＝０．６９３÷濁度の変化速度）を算出する濁度変化の半減期算出ステップと、
該濁度の変化速度算出ステップで算出された濁度の変化速度（K）と該濁度変化の半減期算出ステップで算出された濁度変化の半減期を記憶ステップに記憶させ、該濁度の変化速度算出ステップで算出された傾き（−１×A）と切片Bを右下がりの直線関係の式つまり一次関数の式（縦軸の値＝−１×A×横軸の値＋B）に代入して、その解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）の間の範囲における横軸に対する縦軸の値をそれぞれ算出し記憶ステップに記憶させ、解析グラフに重ね書きをし、それに付け加えて、右下がりの直線関係の式つまり一次関数の式と濁度の変化速度算出ステップで設定された解析計算開始時間（TAｉｎｉｔｉａｌ）と解析計算終了時間（TAｆｉｎａｌ）と測定時間の最大の範囲における濁度が最小の値を（Ｄｆｉｎａｌ）と濁度の変化速度（K）と濁度変化の半減期算出ステップで算出された濁度変化の半減期のそれぞれを表示する結果表示ステップと、
を有し、
半減期が小さければ乳化状態の液体の安定性が低く、半減期が大きければ乳化状態の液体の安定性が高いことを示す半減期の大小によって、その乳化状態の液体の安定性を評価する方法を具現化し実行するプログラム。
【請求項６】
請求項１記載の方法と請求項４記載の方法の両方を兼ね備え、実験者あるいは使用者の選択により請求項１記載の方法と請求項４記載の方法のどちらか一方の方法を選択し実行することによって乳化状態の液体の安定性を評価する方法
【請求項７】
請求項６記載の方法を具現化し実行できる装置
【請求項８】
請求項３記載のプログラムと請求項６記載のプログラムの両方のプログラムを兼ね備え、実験者あるいは使用者の選択により請求項３記載のプログラムと請求項６記載のプログラムのどちらか一方のプログラムを選択し実行するプログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【公開番号】特開２０１３−６８４８１（Ｐ２０１３−６８４８１Ａ）
【公開日】平成２５年４月１８日（２０１３．４．１８）
【国際特許分類】
【出願番号】特願２０１１−２０６３７６（Ｐ２０１１−２０６３７６）
【出願日】平成２３年９月２１日（２０１１．９．２１）
【特許番号】特許第４９８６１９６号（Ｐ４９８６１９６）
【特許公報発行日】平成２４年７月２５日（２０１２．７．２５）
【出願人】（５０９００６２７７）
【Ｆターム（参考）】
【公開日】平成２５年４月１８日（２０１３．４．１８）
【国際特許分類】
【出願日】平成２３年９月２１日（２０１１．９．２１）
【特許番号】特許第４９８６１９６号（Ｐ４９８６１９６）
【特許公報発行日】平成２４年７月２５日（２０１２．７．２５）
【出願人】（５０９００６２７７）
【Ｆターム（参考）】
[ Back to top ]