(57)【要約】【課題】細胞培養装置によって目的産物が産生されるか否か高い精度で評価できる細胞培養装置の評価方法及び評価プログラムを提供する。【解決手段】本発明に係る細胞培養装置の評価方法は、細胞培養を行う細胞培養装置を評価する方法であり、前記細胞培養装置内の細胞が培養時に受ける一つ又は複数の物理量を独立変数とする関数を定義する定義ステップＳ１と、輸送方程式を用いて一群の前記細胞に対する前記関数の累積積分分布の時間変化を算出する第１算出ステップＳ２と、前記関数の累積積分分布の時間変化を積分時間で除することで、前記関数の累積平均分布の時間変化を算出する第２算出ステップＳ３と、前記関数の累積平均分布の時間変化を前記細胞培養装置内の細胞への影響因子の一つとして前記細胞培養装置を評価する評価ステップＳ４と、をこの順で含む。【選択図】図１
2018161057.tif 000002
【特許請求の範囲】
【請求項１】
細胞培養を行う細胞培養装置を評価する方法であり、
前記細胞培養装置内の細胞が培養時に受ける一つ又は複数の物理量を独立変数とする関数を定義する定義ステップと、
輸送方程式を用いて一群の前記細胞に対する前記関数の累積積分分布の時間変化を算出する第１算出ステップと、
前記関数の累積積分分布の時間変化を積分時間で除することで、前記関数の累積平均分布の時間変化を算出する第２算出ステップと、
前記関数の累積平均分布の時間変化を前記細胞培養装置内の細胞への影響因子の一つとして前記細胞培養装置を評価する評価ステップと、
をこの順で含むことを特徴とする細胞培養装置の評価方法。
【請求項２】
請求項１において、
前記輸送方程式が、下記式（１）であることを特徴とする細胞培養装置の評価方法。
【数１】
2018161057.tif 000015
【請求項３】
細胞培養を行う細胞培養装置を評価するプログラムであり、
コンピュータを、
前記細胞培養装置内の細胞が培養時に受ける一つ又は複数の物理量を独立変数とする関数を定義する定義手段、
輸送方程式を用いて一群の前記細胞に対する前記関数の累積積分分布の時間変化を算出する第１算出手段、
前記関数の累積積分分布の時間変化を積分時間で除することで、前記関数の累積平均分布の時間変化を算出する第２算出手段、
前記関数の累積平均分布の時間変化を前記細胞培養装置内の細胞への影響因子の一つとして前記細胞培養装置を評価する評価手段、
として機能させることを特徴とする細胞培養装置の評価プログラム。
【請求項４】
請求項３において、
前記輸送方程式が、下記式（１）であることを特徴とする細胞培養装置の評価プログラム。
【数２】
2018161057.tif 000016
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、細胞培養装置の評価方法及び評価プログラムに関する。
【背景技術】
【０００２】
現在、医薬品として用いられている抗体をはじめとするタンパク質の多くは遺伝子組換え動物細胞を用いて製造されている。タンパク質は、一次構造までは遺伝子の塩基配列のみに依存するため微生物でも産生できるが、抗体などのタンパク質の多くは、ポリペプチド合成後に糖鎖で修飾される。これらのタンパク質はその一次構造のみでは生理活性が現れないものが多く、微生物では完全な合成が困難であり、二次、三次構造の構築と糖鎖修飾の機能を持つ動物細胞を用いる必要がある。
【０００３】
医薬品の生産に用いられている動物細胞培養の多くは、本来は付着依存性である組織細胞を剥離して順化することで単一細胞として液体培地中で浮遊させ、栄養源や酸素を供給しながら行っている。抗体を用いた抗体医薬品においては、培養液の単位容量当たりの抗体収量が小さいこと、培地が高価であること、一患者当たりの必要投与量が大きいことなどの理由により、非常に高価になっている。このため、細胞数密度が高く、生産収率の大きい大型細胞培養装置による培養が求められている。
【０００４】
大型細胞培養装置では、細胞が吸収する酸素を供給すると共に細胞が吐き出す二酸化炭素を除去するために、十分な量の通気と攪拌が必要である。一方で、過剰な攪拌による流体せん断応力と気泡通気は細胞にダメージを与え、細胞が死滅するという問題がある。これに対応する技術が、特許文献１及び非特許文献１に記載されている。
【０００５】
特許文献１には、接着依存性細胞を培養容器内で単層培養する際の培養細胞評価方法が記載されている。この方法は、接種した細胞が培養容器の底面に接着するまでの細胞接着期と、前記接着した細胞が細胞分裂を開始するまでの誘導期と、前記細胞分裂の開始後から培養容器の底面にほぼコンフルエント状態になるまでの対数増殖期と、前記コンフルエント状態になった後の定常期と、の４種のステージから選ばれる少なくとも１種のステージにおける細胞固有のパラメータを導出し、前記接着依存性細胞の単層培養過程を定量評価することを特徴とする旨記載されている。
【０００６】
非特許文献１には、動物細胞培養槽のスケールアップと最適化設計を目的とした細胞培養槽の数値シミュレーション技術が記載されている。なお、図７は、小型細胞培養装置に関する従来の数値シミュレーション技術を説明する説明図である。図８は、大型細胞培養装置に関する従来の数値シミュレーション技術を説明する説明図である。
【０００７】
非特許文献１には、図７及び図８に示すように、小型細胞培養槽及び大型細胞培養槽内のいずれにおいても、前記した数値シミュレーション技術を用いて、流体せん断応力による細胞死滅１０、混合性能１１及びガス交換性能１２で囲まれた運転可能領域１００を見出すことにより、細胞培養槽の設計最適化を図る旨記載されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００１−２２４３６６号公報
【非特許文献】
【０００９】
【非特許文献１】天野研、外４名、「抗体医薬品の生産性向上を支援する動物細胞培養槽シミュレーション技術」、日立評論、２００７年５月、Ｖｏｌ．８９，Ｎｏ．０５、ｐ．３４〜３７
【発明の概要】
【発明が解決しようとする課題】
【００１０】
実生産向けの大型細胞培養装置（２００Ｌから１０００００Ｌ規模）の適切な運転条件の評価にあたっては、事前に１ｍＬから１００Ｌ規模の小型細胞培養装置を作製し、当該小型細胞培養装置での細胞培養結果を基にして、運転条件を評価する。特許文献１に記載されている技術や非特許文献１に記載されている技術を適用すれば、細胞をより確実に増殖させることができる。
【００１１】
しかしながら、小型細胞培養装置での培養結果として得られる評価指標としての物理量のうち、特許文献１においては、細胞の増殖速度しか考慮できない。
また、非特許文献１においては、評価指標として取り入れるべき物理量として、流体せん断応力による細胞の死滅１０、混合性能１１、ガス交換性能１２、が挙げられているが、これらの物理量はいずれも、特定の時間間隔ごとに測定した、物理量の離散的な経時変化である。
【００１２】
本発明者が研究を進めたところ、細胞の増殖速度や物理量の離散的な経時変化では細胞増殖についての現象は説明したり、評価したりできるものの、細胞の遺伝子発現などの細胞の内部状態に関する現象は説明できないことが新たに分かった。つまり、従来の技術では、細胞が増殖しても、タンパク質などの目的産物が産生されない場合があることが分かった。
【００１３】
本発明は前記状況に鑑みてなされたものであり、細胞培養装置によって目的産物が産生されるか否か高い精度で評価できる細胞培養装置の評価方法及び評価プログラムを提供することを課題とする。
【課題を解決するための手段】
【００１４】
本発明者は、前記課題を解決するため鋭意研究した結果、細胞が受ける物理量の累積値、又は累積値の時間平均値の時間変化を、物理量をソース項とした輸送方程式から算出できることを見出した。そして、輸送方程式で算出した結果に基づいて、細胞が受ける物理量の累積値、又は累積値の時間平均値の時間変化を算出し、細胞培養装置内の細胞への影響因子の一つとして用いることで、細胞培養装置の運転条件を評価できることを見出した。本発明はこれらの知見に基づいて成されたものである。
【００１５】
前記課題を解決した本発明に係る細胞培養装置の評価方法は、細胞培養を行う細胞培養装置を評価する方法であり、前記細胞培養装置内の細胞が培養時に受ける一つ又は複数の物理量を独立変数とする関数を定義する定義ステップと、輸送方程式を用いて一群の前記細胞に対する前記関数の累積積分分布の時間変化を算出する第１算出ステップと、前記関数の累積積分分布の時間変化を積分時間で除することで、前記関数の累積平均分布の時間変化を算出する第２算出ステップと、前記関数の累積平均分布の時間変化を前記細胞培養装置内の細胞への影響因子の一つとして前記細胞培養装置を評価する評価ステップと、をこの順で含む。
【００１６】
また、本発明に係る細胞培養装置の評価プログラムは、細胞培養を行う細胞培養装置を評価するプログラムであり、コンピュータを、前記細胞培養装置内の細胞が培養時に受ける一つ又は複数の物理量を独立変数とする関数を定義する定義手段、輸送方程式を用いて一群の前記細胞に対する前記関数の累積積分分布の時間変化を算出する第１算出手段、前記関数の累積積分分布の時間変化を積分時間で除することで、前記関数の累積平均分布の時間変化を算出する第２算出手段、前記関数の累積平均分布の時間変化を前記細胞培養装置内の細胞への影響因子の一つとして前記細胞培養装置を評価する評価手段、として機能させる。
【発明の効果】
【００１７】
本発明に係る細胞培養装置の評価方法及び評価プログラムは、細胞培養装置によって目的産物が産生されるか否か高い精度で評価できる。
【図面の簡単な説明】
【００１８】
【図１】本発明の細胞培養装置の評価方法の一実施形態を説明するフローチャートである。
【図２】本発明の細胞培養装置の評価プログラムの一実施形態を説明するブロック図である。
【図３】物理量の存在比率を表すグラフである。
【図４】累積平均物理量の存在比率を表すグラフである。
【図５】累積平均物理量を評価指標の一つとして加えた場合における小型細胞培養装置の運転可能領域を説明する説明図である。図中、横軸は攪拌回転数を表し、縦軸は通気量を表している。
【図６】累積平均物理量を評価指標の一つとして加えた場合における大型細胞培養装置の運転可能領域を説明する説明図である。図中、横軸は攪拌回転数を表し、縦軸は通気量を表している。
【図７】小型細胞培養装置に関する従来の数値シミュレーション技術を説明する説明図である。図中、横軸は攪拌回転数を表し、縦軸は通気量を表している。
【図８】大型細胞培養装置に関する従来の数値シミュレーション技術を説明する説明図である。図中、横軸は攪拌回転数を表し、縦軸は通気量を表している。
【発明を実施するための形態】
【００１９】
以下、適宜図面を参照して本発明に係る細胞培養装置の評価方法及び評価プログラムの一実施形態について詳細に説明する。
【００２０】
［細胞培養装置の評価方法］
図１は、本発明の細胞培養装置の評価方法（以下、「本評価方法」と記載する）の一実施形態を説明するフローチャートである。
本評価方法は、細胞培養を行う細胞培養装置を評価する方法である。なお、「細胞培養装置を評価する」とは、細胞培養装置（細胞培養槽）の形状や、攪拌翼の形状、攪拌翼の回転数、通気量、ｐＨなどによる運転条件を評価するということである。
図１に示すように、本評価方法は、定義ステップＳ１と、第１算出ステップＳ２と、第２算出ステップＳ３と、評価ステップＳ４と、をこの順で含んでおり、この順に各ステップを行う。
【００２１】
ここで、本評価方法の評価対象となる細胞培養装置は、タンパク質などの目的産物を産生する細胞を培養する装置（槽・タンク）であればよく、容量の多少に関わりなくどのような装置も対象となる。細胞培養装置は、細胞培養を行う際に、攪拌翼の攪拌回転数、通気量、温度、気泡径、ｐＨ、表面張力などの運転条件を制御できる手段を有しているものが好ましい。各手段によって制御されるこれらの要素の物理量は、細胞培養装置の操作変数（独立変数）として扱うことができる。
【００２２】
本実施形態における目的産物（産生対象となる物質）としては、医薬品等の主原料となる物質医薬品等の主原料となる物質、例えば、抗体や酵素などのタンパク質や、低分子化合物及び高分子化合物等の生理活性物質などを挙げることができるが、これらに限定されない。
【００２３】
また、培養対象となる細胞としては、例えば、動物細胞、植物細胞、昆虫細胞、細菌、酵母、真菌や藻類などを挙げることができるが、これらに限定されない。本実施形態に係る細胞培養装置は、抗体や酵素などのタンパク質を産生する動物細胞を好適な培養対象とすることができる。
【００２４】
（定義ステップ）
定義ステップＳ１は、細胞培養装置内の細胞が培養時に受ける一つ又は複数の物理量を独立変数とする関数を定義する。
物理量としては、例えば、流体せん断応力、流速、動粘度、乱流エネルギー、乱流エネルギー散逸率、コルモゴロフスケール、酸素濃度、二酸化炭素濃度、培地成分の成分分率、ｐＨ、温度などが挙げられる。複数の物理量を独立変数とする場合、これらを複数組み合わせて得られたものを関数として定義することができる。
【００２５】
（第１算出ステップ）
第１算出ステップＳ２は、輸送方程式を用いて一群の細胞（細胞群）に対する前記関数の累積積分分布の時間変化を算出する。つまり、この第１算出ステップＳ２では、細胞が受けた物理量の累積値を輸送方程式によって計算する。これにより、前記物理量として、例えば、流体せん断応力を考慮した場合、個々の細胞の流体せん断応力累積値を直接計算することなく、特定の流体せん断応力累積値を持つ細胞が、全細胞のどの程度の割合存在するかという、存在比率を知ることができる。
【００２６】
なお、細胞が受けた物理量の累積計算は、本来であれば、個々の細胞を追跡することにより行うのが好ましいと考えられる。しかしながら、細胞培養においては一般的に１０７個／ｍＬオーダーの細胞が存在すると言われており、個々の細胞を追跡して累積計算を行うことは非常に困難である。これに対し、本評価方法では前記したように、細胞が受けた物理量の累積値を輸送方程式によって計算するので、細胞の培地中での動きを個々に追跡して累積計算する場合と比較して容易且つ迅速に計算することができる。
輸送方程式は、下記式（１）であるのが好ましい。
【００２７】
【数１】
2018161057.tif 000003
【００２８】
なお、一般的な輸送方程式は、下記式（２）の形をとり、分子拡散を表す拡散係数Ｄが存在するが、本実施形態で扱う物理量の累積値は、分子拡散する値ではないため、Ｄ＝０とした式（１）を用いる。
【００２９】
【数２】
2018161057.tif 000004
【００３０】
2018161057.tif 000005
【００３１】
【数３】
2018161057.tif 000006
【００３２】
2018161057.tif 000007
【００３３】
【数４】
2018161057.tif 000008
【００３４】
2018161057.tif 000009
【００３５】
2018161057.tif 000010
【００３６】
【数５】
2018161057.tif 000011
【００３７】
重み関数は任意のものを適用できるが、例えば、シグモイド関数（下記式（６）、（７））やヘヴィサイド関数（下記式（８）、（９））などを好適に適用できる。
【００３８】
【数６】
2018161057.tif 000012
【００３９】
2018161057.tif 000013
【００４０】
ここで、図７及び図８に示す従来技術による評価を参照して、図５及び図６に示す本評価方法による評価の優位性を説明する。
前述したように、実生産向けの大型細胞培養装置の適切な運転条件の評価にあたっては、事前に小型細胞培養装置を作製し、当該小型細胞培養装置での細胞培養結果を基にして、運転条件を評価している。つまり、小型細胞培養装置における各種の評価指標に対して、細胞培養可能な運転可能領域を調査している。
【００４１】
すなわち、図７に示すように、従来技術では、細胞培養装置が細胞培養に対して適切な能力を有しているか否かの評価指標として、流体せん断応力による細胞の死滅１０、混合性能１１、ガス交換性能１２の３つを挙げている。なお、図７では、細胞培養装置の操作変数を攪拌回転数と通気量の２種類としている。図７中、流体せん断応力による細胞の死滅１０に関しては運転可能領域１００の上限値、混合性能１１及びガス交換性能１２に関しては運転可能領域１００の下限値を表している。
【００４２】
そして、小型細胞培養装置について得られた実験式又は数値流体解析により求めた攪拌回転数と通気量を操作変数として、図８に示すように、実生産向けの大型細胞培養装置における流体せん断応力による細胞の死滅１０、混合性能１１、ガス交換性能１２からなる３つの評価指標に基づいた培養装置の運転可能領域１００を定め、これを調査している。なお、図８に示す実生産向け大型細胞培養装置の運転可能領域１００は、図７の小型細胞培養装置の運転可能領域１００に対応しており、異なる容量の装置で同様の細胞培養結果が得られる。ここで、細胞培養結果としては、一例として生細胞数などの物理量が挙げられる。
【００４３】
しかしながら、このようにして評価・調査された運転可能領域１００で細胞培養を行うと、前述したように、生細胞数は増えるが、目的とするタンパク質が発現しないことがある。これは、流体せん断力などの物理量（物理的な刺激）が細胞内の情報伝達系に作用して、目的とするタンパク質の発現を妨げてしまうことによるものと考えられる。
【００４４】
2018161057.tif 000014
【００４５】
［細胞培養装置の評価プログラム］
図２は、本発明の細胞培養装置の評価プログラム（以下、「本評価プログラム」と記載する）の一実施形態を説明するブロック図である。
本評価プログラムは、細胞培養を行う細胞培養装置を評価するプログラムである。図２に示すように、本評価プログラムは、コンピュータ５を用いて前記した本評価方法を実践するため、コンピュータ５を定義手段１、第１算出手段２、第２算出手段３、評価手段４として機能させる。
ここで、本評価プログラムにおける定義手段１、第１算出手段２、第２算出手段３、評価手段４は、それぞれ図１を参照して説明した本評価方法の定義ステップＳ１、第１算出ステップＳ２、第２算出ステップＳ３、評価ステップＳ４に対応するものであるため、各手段の意義や動作、作用等についての詳細な説明は省略する。
【００４６】
本評価プログラムは、ＣＤ−ＲＯＭ、フレキシブルディスク等のコンピュータ読み取り可能な記録媒体（不図示）に記録し、記録媒体駆動装置（不図示）によって、当該記録媒体から本評価プログラムを読み出して図示しないハードディスクドライブなどの記憶手段にインストールして実行することができる。
【００４７】
また、コンピュータ（クライアント）が通信ネットワークなどの通信手段を備えている場合、本評価プログラムが通信ネットワークを介して接続された他のコンピュータ（サーバ）に記憶され、当該コンピュータから通信ネットワークを介して本評価プログラムをダウンロードして実行させたり、サーバに記憶された本評価プログラムを実行させたりすることができる。
【００４８】
これらのようにして本評価プログラムによりコンピュータ５を機能させると、評価対象となる細胞培養装置が目的産物を産生できるか否か高い精度で評価できる。なお、この場合、数値解析した結果をサーバに備えられた記憶手段（不図示）に記憶することとしてもよい。
【実施例】
【００４９】
次に、実施例により本発明の内容を具体的に説明する。
本実施例では、細胞培養装置を用いて細胞を培養し、所定の培養期間における生細胞数とタンパク質を測定した。
【００５０】
細胞としてはマウスマウスハイブリドーマであるＣＲＬ−１６０６細胞（Ａｍｅｒｉｃａｎ  Ｔｙｐｅ  Ｃｕｌｔｕｒｅ  Ｃｏｌｌｅｃｔｉｏｎより購入）を用いた。この細胞は抗フィブロネクチン抗体を分泌する浮遊系の細胞である。
培養にはＩｓｃｏｖｅ’ｓ  Ｍｏｄｉｆｉｅｄ  Ｄｕｌｂｅｃｃｏ’ｓ  Ｍｅｄｉｕｍ（ＩＭＤＭ培地）に５％濃度になるようにＦｅｔａｌ  ｂｏｖｉｎｅ  ｓｅｒｕｍ（ＦＢＳ）を添加した培地を用いた。
【００５１】
細胞培養装置は、内径１５０ｍｍ、培養容積５Ｌのガラス製円筒型培養槽を用いた。この細胞培養装置には、加温用ラバーヒータ、マグネット駆動式攪拌翼、温度測定電極、ｐＨ電極、溶存酸素（ＤＯ）電極及びこれらを計測して調節する制御装置が接続されている。液中通気用に、平均細孔径１００μｍの焼結金属製スパージャを組み込んだ。
【００５２】
培養液のｐＨは、細胞培養装置の気相部に供給する混合ガス（空気、窒素、酸素、炭酸ガス）中の炭酸ガス濃度を増減させて自動で調節した。培養液の温度は３７℃に調節した。培養液が酸性化した場合には、濃度２％の水酸化ナトリウム溶液を注入して調節した。また、ＤＯ濃度については、気相部及び液中に供給するガスの酸素分圧を増減して調節した。攪拌翼の回転数は、攪拌翼を取り付けた駆動軸が直結されている駆動モータの出力を調節し、１００ｒｐｍで回転させた。
【００５３】
培養開始から１日に２〜３回の頻度で細胞培養装置内の培養液を無菌的にサンプリングし、生細胞数を測定した。
生細胞数の測定は、Ｂｅｃｋｍａｎ  Ｃｏｕｌｔｅｒ社製の生細胞計数計ＢｉｏＰｒｏｆｉｌｅ  １００Ｐｌｕｓを用い、トリパンブルーにより染色した死細胞と区別して画像処理により計数した。また、サンプリング液を二次元電気泳動にかけて、細胞内の全タンパク質を数千以上のスポットに分離した。そして、各スポット中のタンパク質をトリプシンで分解して質量分析で解析することでペプチドの質量や部分配列情報を得、ペプチドマスフィンガープリント法により細胞内のタンパク質を同定した。
【００５４】
その結果、攪拌回転数と通気量を操作することで、細胞に含まれるタンパク質、すなわち、細胞の培養過程で発現したタンパク質の種類及び比率が変化するような攪拌回転数と通気量の組合せが複数存在することが分かった。
【００５５】
検討の結果、これらのタンパク質発現変異点は、実施形態で説明した式（１）、（３）、（４）、（５）、（８）、（９）（但し、式（８）、（９）においては、ｃ＝０）を用いて算出した累積平均物理量の存在比率H(A,t)やそのスカラー値S(A,t)と相関があることが分かった。なお、この実施例では、物理量を流体せん断応力とし、当該流体せん断応力から前記式によって算出されるスカラー値として累積平均流体せん断応力を求めている。
【００５６】
そして、累積平均流体せん断応力によるタンパク質発現変異点１３を図７に重ねてプロットしたものが、前述した図５であり、図８に重ねてプロットしたものが、前述した図６である。
図７と図５及び図８と図６をそれぞれ比較すると、流体せん断応力による細胞の死滅１０、混合性能１１、ガス交換性能１２からなる３つの評価指標に基づいた培養装置の運転可能領域１００においては、生細胞数などでは、同一の結果を得ることができた。
【００５７】
一方、図５及び図６に示すように、タンパク質発現に関しては同一の結果を得るための運転可能領域２００があり、図７及び図８に示す運転可能領域１００に比べて狭くなっている。これは、タンパク質発現に関しては、流体せん断応力による細胞の死滅１０、混合性能１１、ガス交換性能１２からなる３つの評価指標に加えて、累積平均流体せん断応力を評価指標として考慮する必要があることを示している。
【００５８】
また、図５及び図６から、実生産向けの大型細胞培養装置において、小型細胞培養装置と同じタンパク質発現を実現するには、流体せん断応力による細胞の死滅１０、混合性能１１、ガス交換性能１２、累積平均流体せん断応力によるタンパク質発現変異点１３の４つのパラメータを一致させる必要があることが分かる。累積平均流体せん断応力によるタンパク質発現変異点１３の計算は本発明による一連の手順で実現できる。
【００５９】
以上、本発明に係る細胞培養装置の評価方法及び評価プログラムについて、実施形態及び実施例によって詳細に説明したが本発明の主旨はこれに限定されるものではなく、様々な変形例が含まれる。例えば、前記した実施形態は本発明を分かり易く説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施形態の構成の一部を他の実施形態の構成に置き換えることが可能であり、また、ある実施形態の構成に他の実施形態の構成を加えることも可能である。また、各実施形態の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
また、前記した各構成、機能、処理部、処理手段、制御手段等は、それらの一部又は全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリや、ハードディスク、ＳＳＤ（Solid State Drive）等の記録装置、又は、ＩＣカード、ＳＤカード、ＤＶＤ等の記録媒体に置くことができる。
また、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。
【符号の説明】
【００６０】
Ｓ１    定義ステップ
Ｓ２    第１算出ステップ
Ｓ３    第２算出ステップ
Ｓ４    評価ステップ
１      定義手段
２      第１算出手段
３      第２算出手段
４      評価手段
５      コンピュータ
【図１】
2018161057.tif 000017
【図２】
2018161057.tif 000018
【図３】
2018161057.tif 000019
【図４】
2018161057.tif 000020
【図５】
2018161057.tif 000021
【図６】
2018161057.tif 000022
【図７】
2018161057.tif 000023
【図８】
2018161057.tif 000024