(57)【要約】  最適なパイプスチル運転の決定のための方法であって、沸点範囲の画分に分けられる原油流れを、パイプスチル内へ供給する工程、予測される沸点範囲の画分の生産高を決定するために、原油流れの仮想分析を行なう工程、予測される沸点範囲の画分の生産高を、パイプスチルからの実際の沸点範囲画分の生産高と比較して、これら画分の生産高間の差を決定する工程、および画分の生産高間の差をパイプスチルの運転と関連付ける工程を含む方法。【選択図】図１
2008510877.tif 000002
【特許請求の範囲】
【請求項１】
最適なパイプスチル運転を決定する方法であって、
ａ）原油流れを前記パイプスチルに供給し、そこで前記原油流れを沸点範囲の画分に分ける工程；
ｂ）予測される沸点範囲の画分の生産高を決定するために、前記原油流れの仮想分析を行なう工程；
ｃ）予測される沸点範囲の画分の生産高を、パイプスチルからの実際の沸点範囲画分の生産高と比較して、これら画分の生産高間の差を決定する工程；および
ｄ）前記画分の生産高間の差をパイプスチルの運転と関連付ける工程
を含むことを特徴とする方法。
【請求項２】
前記仮想分析は、前記原油流れのＦＴ−ＩＲスペクトルを用いて行なわれることを特徴とする請求項１に記載の方法。
【請求項３】
前記仮想分析は、前記原油流れのＦＴ−ＩＲスペクトルおよびＡＰＩ比重を用いて行なわれることを特徴とする請求項１に記載の方法。
【請求項４】
仮想分析は、原油流れのＦＴ−ＩＲスペクトル、ＡＰＩ比重および粘度を用いて行なわれることを特徴とする請求項１に記載の方法。
【請求項５】
仮想分析は、パイプスチルへの入口ラインから採取された原油のサンプルに対して行なわれることを特徴とする請求項１に記載の方法。
【請求項６】
仮想分析は、パイプスチルへの入口ラインに供給する良く混合された原油タンクから採取された原油のサンプルに対して行なわれることを特徴とする請求項１に記載の方法。
【請求項７】
仮想分析は、パイプスチルへの入口ラインに供給する不均質な原油タンクから採取された全レベルの原油サンプルに対して行なわれることを特徴とする請求項１に記載の方法。
【請求項８】
仮想分析は、パイプスチルへの入口ラインに供給する二つ以上の良く混合された原油タンクから採取されたサンプルに対して行なわれることを特徴とする請求項１に記載の方法。
【請求項９】
パイプスチル運転における偏差の源としての、原油の変動を除去する工程を更に含むことを特徴とする請求項１に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はパイプスチルに供給される特定の原油流れに対して、原油のパイプスチルが最適に作動しているかどうかを決定する方法に関する。
【背景技術】
【０００２】
全ての原油はそれらの沸点範囲の画分において量の変化をきたし、各画分は存在する特定の分子種により決定される異なった物理的特性を有するであろう。これらの二つの要因、体積と物理的性質の組合せは原油の全体的品質を決定し、該材料の価格を決定する重要な要素である。原油の品質は該原油が加工される時に、精製所の運用設定を定義するためにも用いられる。
【０００３】
石油化学工業において、原油の品質は伝統的に原油分析法を用いて評価されていた。原油が分析試験されるとき、それは二段階で蒸留される。ＡＴＳＭ  Ｄ２８９２のような方法（ＡＴＳＭ規格の年鑑、５．０１〜５．０３巻、米国材料試験協会（Ａｍｅｒｉｃａｎ  Ｓｏｃｉｅｔｙ  ｆｏｒ  Ｔｅｓｔｉｎｇ  ａｎｄ  Ｍａｔｅｒｉａｌｓ）、フィラデルフィア、ペンシルベニア州（ＰＡ）を参照）が約６５０°Ｆ（３４３℃）未満で沸騰する蒸留液カットを分離するために使用される。この蒸留からの残留物は更に、６５０°Ｆから約１０００〜１０５４°Ｆ（３４３℃から５３８〜５６８℃）の範囲をカバーする蒸留カット分および真空残渣油カットを生み出すために、ＡＴＳＭ  Ｄ５２３６のような方法を用いて蒸留される。少なくともＬＰＧ（初留点６８°Ｆ以下）、ＬＳＲ（６８〜１５５°Ｆ）、ナフサ（１５５〜３５０°Ｆ）、灯油（３５０〜５００°Ｆ）、ディーゼル油（５００〜６５０°Ｆ）、減圧軽油（６５０°Ｆから１０００〜１０５４°Ｆ）、および真空残渣油（１０００〜１０５４°Ｆ＋）を含む、代表的な製品または個別供給に対応するカット分は一般的に分離される。各々の蒸留カットは次に元素的、分子的、物理的および／または性能特性に関して解析される。実施される特定の解析はカットの代表的な性質に依存する。これらの解析から由来するデータは一般に、あらゆる所望の蒸留範囲に関する原油の品質を評価するために数学的に操作されることができる、電子データベースに格納されるであろう。商品原油分析ライブラリーは、いずれもアスペンテック社（Ａｓｐｅｎｔｅｃｈ  Ｉｎｃ．）と同じように、データを操作するためのツールを備えた、ハバリー・システム社（Ｈａｖｅｒｌｙ  Ｓｙｓｔｅｍｓ  Ｉｎｃ．）およびＨＰＩコンサルタンツ社（ＨＰＩ  Ｃｏｎｓｕｌｔａｎｔｓ  Ｉｎｃ．）から利用可能である。分析データはクルード・クォリティー社（Ｃｒｕｄｅ  Ｑｕａｌｉｔｙ  Ｉｎｃ．）、シェル石油会社（Ｓｈｅｌｌ  Ｏｉｌ  Ｃｏｍｐａｎｙ）、およびスタットオイル（Ｓｔａｔｏｉｌ）により発行されている。特性対蒸留温度のデータは一般に滑らかな曲線に合致し、それは次にあらゆる所望の蒸留カットに関する特性を評価するために使用することができる。
【０００４】
詳細な原油分析は完了まで数週間から数ヶ月かかることがある。その結果、ビジネスの決定をするため、および計画、管理および作業を最適化するために使用される分析データは、一般的に最近購入、販売または加工されている貨物からではなく、むしろ古いデータとなる。分析は操業に重要な影響を持ち得る貨物間の変動を考慮しない。非特許文献１は、油の品質変動の原因、精製所の作業に対するそれらの影響、および原油品質モニタリングに用いるための改善された解析技術の必要性について論じている。ワグスパックは時間に対する、およびその輸送期間中の原油の変動性の原因を、製造用タンクの老化、相対的な分野の生産割合の変化、集油システムにおける原油の混合、バッチ・インターフェースに対するパイプラインの劣化、汚染、および共通仕様の原油流の中への著しく異なる品質の油流の注入としてリストアップしている。そのような変動は、原油価格およびそこから生産され得る製品において大きな変化を生じる可能性がある。
【０００５】
パイプスチルとも呼ばれる原油精製ユニットは、原油を異なる沸点温度（カットポイント）における、それらの構成する沸点範囲の部分別に分け、それらは次に別の精製プロセスユニットへの、または精製された石油製品へのブレンド用の供給材料となる。それぞれのカットポイントは、各沸点範囲の部分において入手出来ると見込まれる材料の量と同様に、経済的要因から決定される。精製所の作業は最新の原油分析を用いて、プラント操業の高度化された数学的モデルにより決定されるため、最高の価値の油流および製品を最大限に回収するために最適化される。
【０００６】
最適な作業からの逸脱は高価となる可能性を有し、ユニットはそれらを作業の目標値内に保つため絶えずモニターされる。逸脱が観察されると、プラントのスタッフはそれらが修正されるように、内在する原因を理解するよう試みる。これらの逸脱には多くの考えられる原因がある。蒸留塔内部および／または関連する熱交換器設備の汚損、タワー内部の機械的損傷、および故障した器具類等の機械的問題を含むことがある。逸脱はまた誤った制御設定により発生し得る。逸脱の根本的原因を識別することは困難で時間のかかる仕事であるかも知れない。解析を複雑にするのは、実験室での分析を用いて最適な作業が決定される間に、納入された原油の品質が分析において規定されたものから、時には大きく逸脱し得ることである。更に、原料油流れはしばしば異なる原油の混合物であり、該混合物の各原油の正確な割合は高い精度では知り得ないことがある。プラントのスタッフは逸脱が次善のプラント運転によるものか、あるいは原油品質および／または原油混合の構成における正常な変動の結果であるかを決定しなければならない。この不確実性は、連続した次善の運転においてもたらされる内在した操業上の問題を修正することに対して、遅れまたは不作為を生じ得る。観察された逸脱に対して内在する原因としての原油品質を確認または排除する能力は、それゆえ問題の解決を加速させることができる。
【０００７】
【特許文献１】米国特許第６，６６２，１１６Ｂ２号明細書
【特許文献２】米国特許第５，１２１，３３７号明細書
【非特許文献１】Ｋ．Ｇ．ワグスパック（炭化水素の加工、７７巻９号、１９９８年、特集記事）（Ｋ．Ｇ．Ｗａｇｕｅｓｐａｃｋ（Ｈｙｄｒｏｃａｒｂｏｎ  Ｐｒｏｃｅｓｓｉｎｇ，７７（９），１９９８  Ｆｅａｔｕｒｅ  Ａｒｔｉｃｌｅ））
【非特許文献２】Ｃ．Ｌ．ローソンおよびＲ．Ｊ．ハンソン（最小二乗法問題の解、サイアム、１９９５年）（Ｃ．Ｌ．Ｌａｗｓｏｎ  ａｎｄ  Ｒ．Ｊ．Ｈａｎｓｏｎ  （Ｓｏｌｖｉｎｇ  Ｌｅａｓｔ  Ｓｑｕａｒｅｓ  Ｐｒｏｂｌｅｍｓ，  ＳＩＡＭ，１９９５））
【非特許文献３】Ｒ．ブロおよびＳ．Ｄｅ  ジョング（ケモメトリックス・ジャーナル、１１巻、３９３〜４０１頁、１９９７年）（Ｒ．Ｂｒｏ  ａｎｄ  Ｓ．Ｄｅ  Ｊｏｎｇ  （Ｊｏｕｒｎａｌ  ｏｆ  Ｃｈｅｍｏｍｅｔｒｉｃｓ，Ｖｏｌ．１１，３９３−４０１，１９９７））
【発明の開示】
【発明が解決しようとする課題】
【０００８】
原油精製ユニットの最適な運転は、最適値からの逸脱を識別できるように、該ユニットに供給される原油の品質をプラントのスタッフが正確に説明することを必要とする。これはそのとき逸脱に関して内在する原因が適切に調査され、必要に応じて作業の変更を通じて修正されることを可能にする。原油の品質は一般に実験室での分析を原油に対して行なうことにより決定される。原油の品質は非常に変りやすく、実験室での分析により貨物を日常的に測定することは、それらの比較的高価な費用および実験室での分析を実施するためにかかる時間のために実際的ではない。原油ユニットに供給される材料から期待される実際の生産高を正確に説明することの不可能さは、解析に不確実性を加え、正しくない結論をもたらす可能性がある。
【課題を解決するための手段】
【０００９】
本発明は実験室での分析の代わりに、原油に対して仮想分析を実施することにより、与えられた原油流れに対してパイプスチルが最適に作動しているかどうかを決定する方法である。これは、ここに参考として盛り込まれている特許文献１に記載されている多変量解析データが原油に対して得られることを必要とする。特許文献１の方法は、これからは仮想分析（Ｖｉｒｔｕａｌ  Ａｓｓａｙ）と呼ばれよう。従って、本発明はパイプスチルに対して何らかの物理的または機械的試験を実施する前に、パイプスチルの「健康度」を決定することを可能にする。
【００１０】
仮想分析は、この偏差解析において用いられる原油の品質を正確かつ迅速に決定するための方法を用意することにより、この不確実性を克服する。特許文献１に記載された仮想分析は、分光法などの多変量解析技術または、多変量解析技術と検査の組合せを用いて未知の材料を解析するための方法である。そのような検査は大量の材料に対して安価で容易に出来る、物理的または化学的特性の測定であり、ＡＰＩ（米国石油協会）比重もしくは比重、および粘度を含むが、それに限定されない。未知の材料はその多変量解析データ（例えばスペクトル）、またはその多変量解析データと検査を、同じタイプの基準材料に関する多変量解析データ、または多変量解析データと検査データを含むデータベースと比較することにより解析される。該比較は未知の材料について含んでいる多変量解析データ、または含んでいる多変量解析データと検査に適合する基準材料の部分集合の混合を計算するために行なわれる。基準材料の計算された混合は、次に基準材料および既知の混合関係性の測定された化学的、物理的および性能特性を用いて、未知の材料の追加の化学的、物理的または性能上の特性を予測するために用いられる。
【００１１】
原油ユニットの健康モニタリングに仮想分析を適用する際に、解析者が仮想分析予測は健康度の決定に対して妥当であるかどうか決定することを可能にするため、該適用に関する仮想分析デ−タベース（ライブラリー）を最適化することが必要である。該最適化は検査の重み付けの設定および、解析の範囲を定義するための、予測される原油の特性が健康度決定における使用に対して十分であると見なされる、ティア１フィット（Ｔｉｅｒ  １  ｆｉｔｓ）と称される適合品質カットオフを含む。
【発明を実施するための最良の形態】
【００１２】
原油ユニットのプロセスのモニタリングは精製最適化モデル、日程計画の適用、または分析の伝達ツールを用いて、実際の生産高および鍵となる品質と、予測されたそれらとを比較する。生産高および品質を予測するための適切なソフトウェア・パッケージは、マスプロ社により販売されているアドバンス・リファイナリー・モデリングシステム（ｔｈｅ  Ａｄｖａｎｃｅ  Ｒｅｆｉｎｅｒｙ  Ｍｏｄｅｌｉｎｇ  Ｓｙｓｔｅｍ  ｓｏｌｄ  ｂｙ  ＭａｔｈＰｒｏ，Ｉｎｃ．）、アスペンテックにより販売されているオリオンＴＭおよびＰＩＭＳＴＭ（ｔｈｅ  ＯＲＩＯＮＴＭ  ａｎｄ  ＰＩＭＳＴＭ  ｓｏｌｄ  ｂｙ  ＡｓｐｅｎＴｅｃｈ）、およびＨＰＩコンサルタンツ社により販売されているアッセイ・シミュレータ（Ａｓｓａｙ  Ｓｉｍｕｌａｔｏｒ  ｓｏｌｄ  ｂｙ  ＨＰＩ  Ｃｏｎｓｕｌｔａｎｔｓ，  Ｉｎｃ．）を含むが、それらには限定されない。多くの石油会社は類似の「社内」システムを有する。次に、逸脱はそれらが適切に構成されていない実際のユニット運転、装置の問題、または単に期待と異なる供給物の品質に起因するかどうかを決定するために調査される。
【００１３】
仮想分析なしでは、プラントは供給物の品質が逸脱の原因かどうか決定するためのデータを殆ど持たないであろう。実験室での分析は比較的高価であり、完了するためには数週間から数ヶ月を要する。それゆえ受け取った原油の各貨物またはバッチに対して実験室での分析を実施することは実際的でない。プラントのスタッフは従って、原油が数ヶ月または数年前に実施され得た、最新の実験室での分析に基づくような期待された品質から異なるかどうかを決定するために、受け取った貨物のＡＰＩ比重のような、容易に測定できる特性に依存する。しかしながら、ＡＰＩ比重における差異は原油が変化したことを示し得る一方で、それは原油の特定の産出パターンがどのように変化したかを決定するには不十分である。この不確実性は精製スタッフが、原油の品質はプラントの逸脱による原因らしいと結論づけることをもたらし、それによって他の内在する原因を効果的に覆い隠し、未解決に至らせる可能性がある。仮想分析は与えられた原油ユニットの供給物の産出パターンを、高い精度で迅速かつ安価に決定する能力を備えることにより、この不確実性を取り除く。
【００１４】
仮想分析は、この偏差解析において用いられる原油の品質を正確かつ迅速に決定する方法を備えることにより、この不確実性を克服する。特許文献１に記述されている仮想分析は、分光法のような多変量解析技術、または多変量解析技術と検査の組合せを用いて未知の材料を解析する方法である。そのような検査は大量の材料に対する実験室での分析に比べて容易かつ安価に出来る、物理的または化学的特性測定であり、ＡＰＩ比重もしくは比重および粘度を含むが、それに限定されない。未知の材料はその多変量解析データ（例えばスペクトル）、またはその多変量解析データと検査を、同じタイプの基準材料に関する多変量解析データ、または多変量解析データと検査データを含むデータベースと比較することにより解析される。該比較は未知の材料について含んでいる多変量解析データ、または含んでいる多変量解析データと検査に適合する基準材料の部分集合の混合を計算するために行なわれる。基準材料の計算された混合は、次に基準材料および既知の混合関係性の測定された化学的、物理的および性能上の特性を用いて、未知の材料の追加の化学的、物理的または性能上の特性を予測するために用いられる。
【００１５】
仮想分析は７０００〜４００ｃｍ−１スペクトル域において、ＦＴ−ＭＩＲスペクトルデータを利用することが望ましい。スペクトルはスペクトル望ましくはＣａＦ２ウィンドウを有する０．２５ｍｍの公称パスレングス・セルを用いて集められる。不連続のスペクトル域は、吸光度が分光計の線形応答範囲を超えるデータ、およびスペクトル変動および従って情報量が小さい領域を避けるために、一般にこのスペクトルから選択される。スペクトルデータは、ブラウン（Ｂｒｏｗｎ）において述べられた直交化手順（特許文献２および特許文献１）を用いて、余分な信号に関して修正される。これらの余分な信号はサンプルの有機物成分から発生しない信号を表わし、基線の変動、水蒸気などの分光計パージ混入物質に起因する吸収、および原油の場合、未精製のサンプル中に溶解または分散している水分を含む。修正されたスペクトルデータは望ましくは検査データで増大される。検査データは直線的に混合できる形に変換され、重み付けされ、そして分光ベクトルの終端に連結される。例えば、ＡＰＩ比重は比重に変換され、粘度は粘度混合ナンバーに変換されるであろう。解析された未知の原油に関して増大され、修正された分光ベクトルは、望ましくは高速非負最小二乗アルゴリズム（Ｆａｓｔ  ＮｏｎＮｅｇａｔｉｖｅ  Ｌｅａｓｔ  Ｓｑｕａｒｅｓ  ａｌｇｏｒｉｔｈｍ）を用いて、基準原油に関して増大され修正された分光ベクトルの直線的な組合せとして適合させられる。一つの適切なアルゴリズムは非特許文献２に記述されている。一つの望ましいアルゴリズムは非特許文献３に記述されている。高速非負最小二乗アルゴリズムは、混合のための係数が１に十分近い値に合計されるまで、増大されたベクトルのスペクトル部分のスケールを調整する、反復アルゴリズムの中で使用され得る。該解析は指示された割合で加えられた時に、その増大された分光ベクトルが、解析された未知の原油に関する増大された分光ベクトルと最も近く一致する、仮想混合（Ｖｉｒｔｕａｌ  Ｂｌｅｎｄ）と称される基準原油の処方を生み出す。仮想分析は、それらに関して同じ指示された割合での基準原油が既知の混合関係を用いる、分析データを混合することにより生み出される。予測は材料の本当の混合に関して品質を計算するために設計されたソフトウェアを用いて行なわれ得る。これらの「混合」計算を行なうことができるソフトウェアは、ハバリー・システム社（Ｈａｖｅｒｌｙ  Ｓｙｓｔｅｍｓ  Ｉｎｃ．）、ＨＰＩコンサルタンツ社（ＨＰＩ  Ｃｏｎｓｕｌｔａｎｔｓ  Ｉｎｃ．）、およびアスペンテック社（Ａｓｐｅｎｔｅｃｈ  Ｉｎｃ．）から市販されている。多くの石油会社は類似の「社内」システムを有する。
【００１６】
統計用ツールは適合品質を評価するために用いられ、それにより分析予測の期待される品質を評価するために用いられる。様々な統計が、解析された未知の原油に関する増大され、修正された分光ベクトルと、基準原油に関する増大され、修正された分光ベクトルの直線的な組合せとの間の一致を測定するために使用され得る。一旦、そのような統計は適合品質比率（Ｆｉｔ  Ｑｕａｌｉｔｙ  Ｒａｔｉｏ）と呼ばれる。該適合品質比率は次の手順により計算される。
【００１７】
ステップ１：Ｒ２の計算
仮想混合が増大されないスペクトルデータを用いて計算される場合、Ｒ２は〔１〕により計算される。
【数１】
2008510877.tif 000003
ｘｕは解析された未知の原油に関する修正された分光ベクトルであり、
【数２】
2008510877.tif 000004
は「仮想混合」に関する計算された修正スペクトルである。ｆは分光ベクトル当たりの周波数ポイントの数であり、ｃは混合に関するゼロ以外の係数の数である。Ｓは係数の合計が１になるように、スペクトルデータを測るために用いられる反復して決定される因子である。転置は上付き文字ｔにより表示される。
【００１８】
ＡＰＩ比重および粘度を増大された分光ベクトルが使用される場合、Ｒ２は次のように計算される。
【数３】
2008510877.tif 000005
λｕ（ａｐｉ）およびλｕ（Ｖｉｓｃ）はＡＰＩおよび粘度の容積測定で混合可能な形であり、ｗＡＰＩおよびｗＶｉｓｃは二つの検査に対する重み付け因子である。
【数４】
2008510877.tif 000006
および
【数５】
2008510877.tif 000007
は仮想混合に基づく計算されたＡＰＩおよび粘度の、推定された混合可能な形である。ＡＰＩおよび粘度の容積測定で混合可能な形が解析において別個に使用される場合には、Ｒ２に関して同様の表現が用いられる。
【００１９】
ステップ２
適合品質ＦＱは次のように計算される。
【数６】
2008510877.tif 000008
ｃは仮想混合における構成要素に関するゼロ以外の係数の数である。指数εは適合品質がＲ２のみに依存するように、ゼロに設定可能であるが、０．２５のオーダーの値に設定することが望ましい。
【００２０】
ステップ３
適合品質比率ＦＱＲは次のように計算される。
【数７】
2008510877.tif 000009
ＦＱＣは適合品質カットオフである。ＦＱＣは、ＦＱＲ≦１．０を伴う解析が意図された適用に対して妥当な精度の予測を生じるように選択される。ＦＱＣに対する値が、解析においてどの検査が用いられたかによって異なるであろうことに注意されたい。ＦＱＲ≦１．０に対する解析はティア１解析と称される。増大されたベクトルにおける検査に関する重み付けは、ティア１フィットに対する所望の精度を達成するために調整される。ＦＱＣを調整する手順および検査に関する重み付けは、付録１において議論される。
【００２１】
仮想分析の解析は図１に示される図に従って行なわれることが望ましい。未知の材料に対するＡＰＩ比重および粘度が測定されたと想定すると、解析図式はポイント１において出発する。ユーザーは解析において使用される基準の特定のセットを供給してもよい。適合は付録１に記述されている三つのステップに従って行なわれる。ＦＴ−ＩＲのみに基づく適合（ステップ１）およびＦＴ−ＩＲとＡＰＩに基づく適合（ステップ２）が計算されるが、それらはこのポイントでは評価されない。ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく適合がティア１フィットを生じる場合、解析は完了し結果が報告される。
【００２２】
ポイント１での解析がティア１フィットを生じない場合、次に該プロセスはポイント２に進む。基準セットは最初に選択された原油と同じ原油等級の、全ての基準を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告される。
【００２３】
ポイント２での解析がティア１フィットを生じない場合、次に該プロセスはポイント３に進む。基準セットは最初に選択された原油と同じ場所から来た全ての基準を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告される。
【００２４】
ポイント３での解析がティア１フィットを生じない場合、次に該プロセスはポイント４に進む。基準セットは最初に選択された原油と同じ地域から来た全ての基準を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告される。
【００２５】
ポイント４での解析がティア１フィットを生じない場合、次に該プロセスはポイント５に進む。基準セットは全ての基準原油を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告される。
【００２６】
ポイント５での解析がティア１フィットを生じない場合、次に該プロセスはポイント６に進む。基準セットは全ての基準原油および混入物質を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告され、そして該サンプルは品質が落ちていると報告される。混入物質が最大許容水準を超えない場合、分析結果は依然として計算されてもよく、信頼区間（Ｃｏｎｆｉｄｅｎｃｅ  Ｉｎｔｅｒｖａｌｓ）が適合ＦＱＲに基づき見積もられる。混入物質が許容水準を超えた場合、その結果はＦＱＲにより示されたものよりも精度が低くなり得る。
【００２７】
ポイント６での解析がティア１フィットを生じない場合、次にＦＴ−ＩＲおよび（各ポイントにおけるステップ２からの）ＡＰＩ比重に基づく適合は、これらのいずれかがティア１フィットを生じるかどうか決定するために調査される。選択された基準に関する適合は最初に調査される（ポイント７）。この解析がティア１フィットを生じた場合、解析は完了し結果が報告される。そうでない場合、該プロセスはポイント８に続き、ＦＴ−ＩＲおよびＡＰＩ比重を用いて選択された原油と同じ等級の原油に基づく適合が調査される。該プロセスはポイント９（同じ場所の原油）、ポイント１０（同じ地域の原油）、ポイント１１（全ての原油）、およびポイント１２（全ての原油と混入物質）に関して適合のチェックを続け、ティア１フィットが見出された場合は中止され、さもなければ継続される。ティア１フィットがＦＴ−ＩＲおよびＡＰＩ比重を用いて見出されない場合、（各ポイントにおけるステップ１からの）ＦＴ−ＩＲのみの適合が、ポイント１３（選択された基準）、ポイント１４（同じ等級）、ポイント１５（同じ場所）、ポイント１６（同じ地域）、ポイント１７（全ての原油）、およびポイント１８（全ての原油と混入物質）に関する適合をチェックしながら調査され、ティア１フィットが見出された場合は中止され、さもなければ継続される。
【００２８】
ティア１フィットが見出されない場合、最高のＦＱＲ値を生み出す解析が選択され、報告される。ＦＱＲ値が１．５以下の場合、その結果はティア２フィットとして報告される。さもなければ、それは不適合として報告される。
【００２９】
粘度のデータが得られない場合、該解析図式はポイント７において出発し、上記に論じられているように続けられるであろう。粘度もＡＰＩ比重も得られなかった場合、該解析図式はポイント１５において出発し、上記に論じられているように続けられるであろう。
【実施例】
【００３０】
実施例１
例１は仮想分析ライブラリーが、パイプスチルの健康モニタリングに使用するためにどのように生み出され、最適化されるかを実証する。計算および最適化の方法論は付録１において与えられている。
【００３１】
仮想分析ライブラリーは、特許文献１および付録１に記述されている方法論を用いた５０４原油に関するＦＴ−ＭＩＲスペクトルを使用して生み出される。４６８５．２〜３４５０．０、２２３８．０〜１５４９．５、および１３４０．３〜１０４５．２ｃｍ−１の領域が使用される。ルジャンドル多項式が各領域において基線の変動を修正するために用いられる。五次（四次）多項式が高周波領域で用いられ、四次（三次）多項式がその他の二つの領域で用いられる。水蒸気および、原油中に拡散した液状水に関して修正がまた生じる。液状水のスペクトルを発生させるために用いられる差スペクトルは、修正ベクトルを生み出す前に、それらのノイズレベルを減らすため平滑化される。合計１７の直交修正ベクトルが１３の多項式、二つの水蒸気修正、および二つの液状水修正を含めて生み出される。５０４原油に関するスペクトルは１７の修正ベクトルに直交させられる。スペクトル変数は容積測定で混合可能な検査により増大される。ＡＰＩ比重は比重に変換され、４０℃における粘度は粘度混合指数に変換される。容積測定で混合可能な検査データは以下に論じられるように重み付けされる。解析の間、未知の原油に関するスペクトルは１７の修正に直交させられ、同じく重み付けされた、容積測定で混合可能な検査で増大され、これらの５０４基準原油に関する増大されたスペクトルの負でない直線的な組合せとして解析される。
【００３２】
比重は再現性により分割され、重み付けパラメータで乗算されることによって重み付けされる。粘度測定に関する再現性は相対的に７％と想定される。粘度混合指数に関する再現性は、粘度混合指数に対して±３．５％で解析され、二つの計算された指標間で絶対的な差を取るサンプルの粘度を変換することにより計算される。該粘度混合指数はこの計算された再現性により分割され、重み付けパラメータを乗じられる。
【００３３】
検査データに関するＦＱＣ値および重み付けパラメータは交差検証法を用いて決定される。各々の５０４原油はライブラリーから取り出され、それがあたかも５０３の残った基準を用いた未知の原油であるかのように解析される。該プロセスは、一つは放置された原油と同じ等級の基準を用いて、一つは放置された原油と同じ場所からの基準を用いて、一つは放置された原油と同じ地域からの基準を用いて、そして一つはライブラリー中の全ての原油を用いて、各原油が解析されるまで５０４回繰り返される。各解析に関して「仮想混合」が計算され、「仮想分析」が予測される。仮想分析予測結果は選択された特性に関して、測定された湿式分析のデータと比較される。パイプスチル健康モニタリングに関して、様々な蒸留カットに対する容積割合生産高は付録１において論じられている手順を用いて予測され、そしてＦＱＣおよび重み付けパラメータを設定するために使用される。
【００３４】
交差検証法は所望の性能が達成されるまで、ＦＱＣおよび重み付けパラメータに関する様々な値を用いて繰り返される。この例として、目標性能は平均生産高予測がティア１解析（ＦＱＲ≦１．０での解析）の時に対して、体積百分率９０％で１．５倍以内であった。交差検証の標準誤差は各々の蒸留カットに関して計算され、適切な統計値ｔで乗算され、そして平均される。検査のための重み付けは、ティア１解析に関するＡＰＩ比重および粘度についての予測誤差が、検査測定の再現性と同等であるように、個々に調整される。この例におけるライブラリーに関して、０．００７９８９のＦＱＣ値とＡＰＩ比重に関する１．４の重み付けパラメータ、および粘度３．２が表１に示すデータを生み出すために使用された。５０４の原油のうち２６０がティア１解析を実現した。
【００３５】
意図された適用によって、様々な蒸留カットおよび追加の特性がＦＱＣの設定において使用可能である。様々な確率レベルもまた統計値ｔを選定するために使用できる。例えば、９５％確率レベルが使用される場合、より少ないティア１フィットしか得られず、しかしＶＡ予測と湿式分析測定の間のより近い一致が達成されるであろう。
【００３６】
【表１】
2008510877.tif 000010
【００３７】
次の例は原油における生産高の変化を決定するために、何故ＡＰＩの変化だけでは不十分であるかを例示している。
【００３８】
実施例２
この例では、同じ等級の原油の二つの貨物が、ＡＰＩ対最新の実験室での分析において著しい差異を示す（０．５を超える差の数字は実験室での再現性の範囲外と考えられ、それゆえ大幅である）。生産高のパターンは仮想分析により決定された。双方の場合とも仮想分析はティア１フィットであったため、統計的に実験室での蒸留と等価である。
貨物１は０．９だけ増加した。
貨物２は１．４だけ増加した。
【００３９】
主な沸点範囲カットに対する生産高を以下に示す。１．５％を超える生産高の差は実験室での再現性の範囲外と考えられ、それゆえ大幅である。
【００４０】
【表２】
2008510877.tif 000011
【００４１】
生産高から見られるように、ＡＰＩ増加から予測できる変化への設定パターンはない。仮想分析を用いることなしには、これらの期待される生産高は定義されないままであろう。仮想分析はそれがなければ可能でない、ユニットからの生産高との比較に関する基礎を提供する。
【００４２】
実施例３
この例においては、各々が実験室での分析に対して（統計的に大きい）０．６ＡＰＩの偏差を有する、別の原油の三つの貨物が解析された。仮想分析は各サンプルに対して実施された。データから見られるように、ＡＰＩが変化しても原油の内部生産高は統計的に顕著な仕方で変わってはいない。仮想分析なしでは、精製所のスタッフは実験室での分析データから予測されたものに対する何らかの実際の生産高偏差は、原油における変化に起因すると結論付けたかも知れない。仮想分析を用いることは、これが誤った仮定らしいことを示す。
【００４３】
【表３】
2008510877.tif 000012
【００４４】
実施例４
これは多量のデータが極めて変わりやすい原油に関して記録された更に極端な例である。データから見られるように、この原油はＡＰＩにおいて広範囲の変動を示し、生じる生産高の変化もまた顕著な変動を示す。ＶＡを用いて、精製所のスタッフは原油ユニットから期待される生産高を追跡することができ、それによって性能をモニターすることができる。ＶＡなしでは、実際の性能と予測された性能の間の差のどれ位が原油品質の変化に起因するかを決定することが困難であろう。
【００４５】
【表４】
2008510877.tif 000013
【００４６】
プラントの利用
以下はどのように仮想分析情報がプラントにおいて運転をモニター／トラブルシューティングするために使用できるかの例である。
・或る精製所は、特定の原油を加工する間に期待された生産量における不足を経験した。仮想分析を用いて、彼らは原油が変化し、そして潤滑油の生産高が原油の結果として低くなり、プラントの運転の結果ではないことを決定することが出来た。
・別の精製所は、常圧パイプスチルからのナフサの生産高において減少を経験した。パイプスチルへの供給物は仮想分析により解析された。その結果、該供給物は実験室での分析に対して大きくは変化していないことが確認され、それは生産量減少の原因として、ユニットの性能を更に調査するように彼らを促した。
・三番目の精製所は、パイプスチルの標準から離れた軽油生産量の予期しない増加を経験した。それらの仮想分析結果を用い、そしてこれらを四番目の精製所からの仮想分析結果と共に確認し、彼らは原油が大きく変っていないことを確認した。これは、それらがプラントの運転スタッフに、増加した軽油生産量の原因はユニットの運転にあることを確信させるのを可能にし、それによって問題がより早く修正されることを可能にした。
【００４７】
作業プロセス
仮想分析を用いて原油ユニットへの供給の正確な表示を発展させることは、原油運搬プロセスにおいて様々なポイントで行なわれる品質測定を利用することができる。仮想分析は下記の原油サンプルに対して実施されてもよい。
・船から精製所への最初の陸揚げの際に採取されるサンプル。
・パイプラインで送られたバッチから精製所または外部の貯蔵タンクへのサンプル。
・タンク内に貯蔵された原油の混合物または前の等級からの最後の残りを含み得る、精製所原油タンクからのサンプル。
・移送ライン上の入口から原油ユニットへのサンプル。
【００４８】
これらの測定結果が作業プロセスにどのように用いられるかは、プラント構成と同様に原油ユニットに供給される混合物の複雑さにより決定される。
【００４９】
仮想分析の品質は変化する可能性があり、それはその結果の品質の内部的に生み出された尺度の条件による。
【００５０】
ティア１フィットの結果は、生産高予測の質における実験室での分析と統計的に等価であり、プラントのスタッフにより原油ユニットの性能予測に使用されるであろう。
【００５１】
方法１−原油ユニット入口移送ラインからの原油測定
サンプルが原油ユニット入口移送ラインからの場合、結果として生じる仮想分析はユニットに供給される原油の実際の品質を表わすものとみなされる。仮想分析情報はその時、原油ユニットが所定の最適な運転包絡線の中で作動しているかどうかを決定するために、実際の原油ユニットの生産高と比較する目的で直接的に使われることができる。予測された運転と観察された運転の間の如何なる偏差も、運転の差異により説明可能であり、実際の原油品質に対する予測された原油品質における未知の偏差としてではない。
【００５２】
方法２−タンクからの原油測定
タンクから採取された原油サンプルよりの仮想分析は、
Ａ）タンクが良く混合されているとみなされ、サンプルがタンクでの混合物全体を表わしている場合、および
Ｂ）全レベルのタンクでのサンプルがＡＳＴＭ  Ｄ４０５７の手順に従って採取されている場合に直接用いられることができる。
【００５３】
一旦タンクでの品質が仮想分析を用いて決定されると、結果として生じた生産高および品質は、原油ユニットが方法１に類似した所定の最適な運転包絡線の中で作動しているかどうかを決定するために使用可能である。
【００５４】
方法３−タンクの混合物からの原油測定
方法２において考慮されるステップは、原油ユニットが単一のタンクから供給される場合にのみ有効である。二つ以上のタンクが原油ユニットに供給するために使用される場合、それらのタンクから結果として生じる原油の混合物を決定しなければならない。これは方法２における仮想分析により計算されるタンクの構成を用いた容量分析計算から行なうことができる。
【００５５】
これら三つの方法のステップは、沸点範囲の画分に分けられる原油を、パイプスチル内へ供給し、予測される沸点範囲の画分を決定するために原油の仮想分析を実施し、予測される沸点範囲の画分を分離された沸点範囲の画分と比較して、二つの画分間の差異を決定し、その差異をパイプスチルの運転と相互に関連付けることを含む。パイプスチルの運転は、その時パイプスチルの出力を予測される出力へと一致させるために修正され得る。
【００５６】
単純化のために、予測された生産高と測定された生産高の間の差は、それらが湿式分析蒸留の見積もられた再現性１．５％を超えた場合に、一般的に大きいと考えられるであろう。しかしながら、パイプスチル性能のより詳細な評価のために、特定のカットに関する予測された生産高と測定された生産高の間の差は、そのカットに関する予測される不確実性（表１のｔ×ＳＥＣ）と、または付録１に記載の手順に従って計算された各生産高予測に関する信頼区間と比較されることができる。
【００５７】
原油ユニットのプロセス・モニタリングは実際の生産高および鍵となる品質を、精製所のＯＲＭモデル、スケジュール適用、または分析の伝達ツールを用いて予測されたものと比較する。次に逸脱は適切に構成されていない実際のユニット運転によるものか、装置の問題か、または単に期待されたものと異なる供給物の品質によるものかを決定するために調査される。
【００５８】
付録１
特許文献１の好適な実施形態において、原油に関する分析データを予測するためにＦＴ−ＩＲスペクトルがＡＰＩ比重および粘度との組合せで使用される。未知の原油のＦＴ−ＩＲスペクトルは検査データで増大され、基準原油に関して増大されたＦＴ−ＩＲスペクトルの直線的な組合せとして適合する。この特許文献１の好適な実施形態は、数学的に〔１〕式で表わされる。
【数８】
2008510877.tif 000014
式中、
【数９】
2008510877.tif 000015
である。
ｘｕは未知の原油に関するＦＴ−ＩＲを含む列ベクトルであり、Ｘは基準原油のＦＴ−ＩＲスペクトルのマトリックスである。ＦＴ−ＩＲスペクトルは、それらが体積ベースで混合されるように原油の一定体積に対して測定される。ｘｕとＸの双方は特許文献１に記載されている修正に対して直交させられていてもよい。ｘｕは列の下方に二つの追加要素、ｗＡＰＩλｕ（ＡＰＩ）およびｗＶｉｓｃλｕ（Ｖｉｓｃ）を加えることにより増大される。λｕ（ａｐｉ）およびλｕ（ｖｉｓｃ）は未知の材料に関するＡＰＩ比重検査および粘度検査の、容積測定で混合可能なバージョンであり、Λ（ＡＰＩ）およびΛ（ｖｉｓｃ）は基準原油に関する、対応する容積測定で混合可能な検査である。ｗＡＰＩおよびｗＶｉｓｃは二つの検査に関する重み付け因子である。
【数１０】
2008510877.tif 000016
および
【数１１】
2008510877.tif 000017
の値は係数ｃｕを伴う、計算された直線的な組合せに基づくスペクトルおよび検査の推定値である。直線的な組合せは非負の最小二乗アルゴリズムを用いて計算されることが望ましい。
【００５９】
特許文献１において、λｕ（ｖｉｓｃ）およびΛ（ｖｉｓｃ）を計算する際に用いられる粘度データは同じ温度で測定される必要があり、
【数１２】
2008510877.tif 000018
の関係を用いて、粘度混合数（Ｖｉｓｃｏｓｉｔｙ  Ｂｌｅｎｄｉｎｇ  Ｎｕｍｂｅｒ）に変換される。
【００６０】
１．５ｃＳｔを超える粘度に対し、パラメータｃは０．６〜０．８の範囲内にある。１．５未満の粘度に対し、ｃは一般に粘度の関数として表わされる。ｃに関する適切な関数は、
ｃ＝０．０９８８６５ｖ４−０．４９９１５ｖ３＋０．９９０６７ｖ２−０．９６３１８ｖ＋０．９９９８８    〔３〕
で与えられる。
【００６１】
特許文献１およびこの特許の目的のため、パラメータａは０に設定され、パラメータｂは１に設定される。粘度が重量ベースで混合すると想定される場合、〔１３〕から計算されるＶＢＮは、容積測定で混合可能な数を得るために材料の比重を乗じられるであろう。容積測定で混合可能な数を得るために用いられる方法は、分析予測を生み出す詳細な解析からのデータを操作するプログラムで使用される方法と合致するように選ばれよう。
【００６２】
基準原油に関する粘度データが、未知の材料に関して粘度が測定される温度において得られない場合、方程式〔１〕は直接適用出来ない。
【００６３】
原油に関して、ＡＳＴＭ  Ｄ３４１（ＡＳＴＭ規格の年鑑、５．０１〜５．０３巻、米国材料試験協会、フィラデルフィア、ペンシルベニア州（ＰＡ）を参照）は、粘度の温度依存性を記述している。この関係を表わす代わりの方法は〔４〕により与えられる。
ＶＢＮ（Ｔ）＝ｌｏｇ（ｌｏｇ（ν（Ｔ）＋ｃ））＝Ａ＋ＢｌｏｇＴ      〔４〕
【００６４】
Ｔは℃または°Ｒにおける絶対温度である。パラメータＡおよびＢは二つ以上の温度において測定された粘度に関する適合〔４〕に基づいて計算される。
【００６５】
未知の材料の粘度が、基準原油に関して粘度データが測定された温度において測定されない場合、二つの代替が適用できる。最初に、未知の粘度が測定された温度におけるＶｒｅｆｅｒｅｎｃｅｓを計算するため、基準原油に関する粘度データに対して方程式〔４〕が適用可能である。基準に関する計算された粘度は、次にΛ（ｖｉｓｃ）を計算するために用いられ、そして方程式〔１〕が適用される。代わりに、〔２〕における勾配ＢがＦＴ−ＩＲスペクトルの解析、またはＦＴ−ＩＲスペクトルとＡＰＩ比重の解析に基づいて推定でき、そして共通の基準温度における未知の材料の粘度を推定するため、Ｂが測定された粘度と組み合わせて使用され得る。
【００６６】
次のアルゴリズム的方法は、未知の材料の解析に対して利点を提供することが見出されている。
ステップ１：
ステップ１において、検査データは使用されない。
【数１３】
2008510877.tif 000019
式中、
【数１４】
2008510877.tif 000020
である
【００６７】
方程式〔４〕は未知の材料のＦＴ−ＩＲスペクトルと合致する直線的な組合せを計算するために、増大されないスペクトルデータに適用される。非負の最小二乗アルゴリズムが、望ましくは係数ｃｓｔｅｐ１を計算するために用いられる。該係数の合計が計算され、スケーリング因数ｓが合計の逆数として計算される。該係数はスケーリング因数により測られる。Ｒ２の値は〔６〕を用いて計算される。
【数１５】
2008510877.tif 000021
ｆは分光ベクトルｘｕ内のポイントの数であり、ｃは適合からのゼロ以外の係数の数である。他の適合度統計がＲ２の代わりに使用可能である。
【００６８】
ステップ２：
ステップ２において、ステップ１からの測られたスペクトルは、ベクトル
【数１６】
2008510877.tif 000022
を形成するために、ＡＰＩ比重データ（即ち比重）の容積測定で混合可能なバージョンを用いて増大される。増大されたベクトル
【数１７】
2008510877.tif 000023
の推定は、ステップ１からの係数、および方程式〔１ｂ〕の関係より計算される。最初のＲ２値は〔７〕を用いて計算される。
【数１８】
2008510877.tif 000024
【数１９】
2008510877.tif 000025
は、ベクトル
【数２０】
2008510877.tif 000026
と同じ長さのベクトルであり、その全ての要素はベクトル
【数２１】
2008510877.tif 000027
における要素の平均である。
【００６９】
測られ、増大された分光ベクトルは次に
【数２２】
2008510877.tif 000028
を用いて適合され、ここで
【数２３】
2008510877.tif 000029
である。
【００７０】
望ましくは非負の最小二乗適合から計算された係数、ｃｓｔｅｐ２は合計され、そして新たなスケーリング因数ｓは、前のスケーリング因数の合計回数の逆数として計算される。該係数は一体として合計するために測られ、増大された分光ベクトルの推定
【数２４】
2008510877.tif 000030
はこれらの正規化された係数および〔８ｂ〕に基づいて再計算される。Ｒ２の値は〔７〕および新たなスケーリング因数を用いて再度計算される。新たなＲ２の値が前の値よりも大きい場合、新しい適合は受け入れられる。方程式〔８〕は新たに計算されたスケーリング因数を用いて再度適用される。本プロセスは計算されたＲ２の値において更なる増加が得られなくなるまで続く。
【００７１】
ステップ３  粘度混合数の使用
単一の固定された温度において測定された粘度に基づく粘度混合数が使用されるべき場合、ステップ３において、最良のＲ２値を与えたステップ２からの、測られ、増大された分光ベクトルは、ベクトル
【数２５】
2008510877.tif 000031
を形成するために、粘度データの容積測定で混合可能なバージョンを用いて更に増大される。増大されたベクトル
【数２６】
2008510877.tif 000032
の推定はｃｓｔｅｐ２、および方程式〔１ｂ〕の関係を用いて計算される。最初のＲ２値は〔９〕を用いて計算される。
【数２７】
2008510877.tif 000033
【数２８】
2008510877.tif 000034
は、
【数２９】
2008510877.tif 000035
と同じ長さのベクトルであり、その要素は
【数３０】
2008510877.tif 000036
における要素の平均である。
【００７２】
測られ、増大された分光ベクトルは次に
【数３１】
2008510877.tif 000037
を用いて適合され、ここで
【数３２】
2008510877.tif 000038
である。
【００７３】
望ましくは非負の最小二乗適合から計算された係数、ｃｓｔｅｐ３は合計され、そして新たなスケーリング因数ｓは、前のスケーリング因数の合計回数の逆数として計算される。該係数は一体として合計するために測られ、増大された分光ベクトルの推定
【数３３】
2008510877.tif 000039
はこれらの正規化された係数および〔１０ｂ〕に基づいて再計算される。Ｒ２の値は〔９〕および新たなスケーリング因数を用いて再度計算される。新たなＲ２の値が前の値よりも大きい場合、新しい適合は受け入れられる。方程式〔１０ａ〕および〔１０ｂ〕は新たに計算されたスケーリング因数を用いて再度適用される。本プロセスは計算されたＲ２の値において更なる増加が得られなくなるまで続く。基準原油の「仮想混合」は最後のｃｓｔｅｐ３係数に基づいて計算され、分析特性は特許文献１に記載された既知の混合関係を用いて予測される。
【００７４】
ステップ２でＡＰＩ比重が得られない場合
ステップ２においてＡＰＩ比重が得られない場合、ステップ１からの測られたスペクトルは、ベクトル
【数３４】
2008510877.tif 000040
を形成するために、粘度データの容積測定で混合可能なバージョンを用いて増大される。増大されたベクトルの推定、
【数３５】
2008510877.tif 000041
はステップ１からの係数、および方程式〔１ｂ〕の関係より計算される。最初のＲ２値は〔１１〕を用いて計算される。
【数３６】
2008510877.tif 000042
【数３７】
2008510877.tif 000043
は、
【数３８】
2008510877.tif 000044
と同じ長さのベクトルであり、その要素は
【数３９】
2008510877.tif 000045
における要素の平均である。
【００７５】
測られ、増大された分光ベクトルは次に
【数４０】
2008510877.tif 000046
を用いて適合され、ここで
【数４１】
2008510877.tif 000047
である。
【００７６】
望ましくは非負の最小二乗適合から計算された係数、ｃｓｔｅｐ２は合計され、そして新たなスケーリング因数ｓは、前のスケーリング因数の合計回数の逆数として計算される。該係数は一体として合計するために測られ、増大された分光ベクトルの推定
【数４２】
2008510877.tif 000048
はこれらの正規化された係数および〔１２ｂ〕に基づいて再計算される。Ｒ２の値は〔１１〕および新たなスケーリング因数を用いて再度計算される。新たなＲ２の値が前の値よりも大きい場合、新しい適合は受け入れられる。方程式〔１２ａ〕および〔１２ｂ〕は新たに計算されたスケーリング因数を用いて再度適用される。本プロセスは計算されたＲ２の値において更なる増加が得られなくなるまで続く。基準原油の「仮想混合」は最後のｃｓｔｅｐ２係数に基づいて計算され、分析特性は特許文献１に記載された既知の混合関係を用いて予測される。
【００７７】
ステップ３代案
上記のステップ３において基準に関する粘度データは、未知の材料に関する粘度が測定される温度において知られるか、または計算可能でなければならない。代わりに、粘度／温度勾配Ｂが推定でき、そして基準原油に関する粘度データが知られている固定の温度における粘度を計算するために使用され得る。
【００７８】
未知の
【数４３】
2008510877.tif 000049
に関する粘度／温度勾配は、ステップ２からの係数ｃｓｔｅｐ２を用いて、基準原油の粘度／温度勾配の混合として推定される。勾配が重量ベースで混合される場合、ｃｓｔｅｐ２係数は基準の比重を用いてそれらの対応する重量割合に変換される。推定された勾配
【数４４】
2008510877.tif 000050
未知の材料に関する粘度νｕ、および粘度が測定される温度Ｔｕが、〔１３〕の関係を用いて固定の温度Ｔｆにおける粘度、νｕ（Ｔｆ）を計算するために使用される。
【数４５】
2008510877.tif 000051
νｕ（Ｔｆ）の値は
【数４６】
2008510877.tif 000052
において使用するための、容積測定で混合可能な粘度値λｕの計算に用いられる。その都度新たな係数ｃｓｔｅｐ３が計算され、勾配
【数４７】
2008510877.tif 000053
は新しい混合に基づいて再び推定され、方程式〔９〕を通じて新たなＲ２の計算に用いられる、νｕ（Ｔｆ）およびλｕの新しい値を計算するために用いられる。
【００７９】
ＡＰＩ比重が得られない場合のステップ２代案
ＡＰＩ比重が得られない場合、νｕ（Ｔｆ）の計算における粘度／温度勾配を推定するために、ステップ３代案の下で上記に説明された手順が係数ｃｓｔｅｐ１を用いて適用される。
【００８０】
追加検査データの組込み
ＡＰＩ比重および粘度に加えた他の検査がオプションとして計算で使用され得る。これらの検査のためのデータの、容積測定で混合可能な形は、増大されたベクトル
【数４８】
2008510877.tif 000054
を形成するため、粘度データと共にステップ２において増大されたベクトル内に含まれる。計算はそのとき上記のように進む。計算の各ステップにおいて、追加検査の予測は〔１４〕により与えられる。
【数４９】
2008510877.tif 000055
【００８１】
含まれる可能性のある他の検査は、硫黄、窒素、および酸価を含むが、それに限定されない。Ｒ２の値は次のように計算されよう。
【数５０】
2008510877.tif 000056
ｉは用いられた検査の数である。
【００８２】
容積測定で混合可能な粘度
ＡＰＩ比重の容積測定で混合可能なバージョンは比重である。ＡＰＩ比重が現在の発明への入力として用いられる場合、それは使用前に比重に変換される。粘度データもまた容積測定で混合可能な形に変換される。
【００８３】
特許文献１は粘度を混合可能な形に変換するために使用できる幾つかの方法を説明している。現在の発明はまた粘度混合指数（ＶＢＩ）の使用を規定している。ＶＢＩは２１０°Ｆにおける粘度に基づく。基準原油に関し、２１０°Ｆにおける粘度は二つ以上の温度で測定された粘度、および方程式〔４〕と〔１３〕の適用に基づいて計算される。未知の材料に関しては、代替のステップ３において用いられるＴｆ値が２１０°Ｆとして選ばれる。粘度混合指数は方程式〔４〕により、２１０°Ｆでの粘度と関連する。
【数５１】
2008510877.tif 000057
【００８４】
与えられた粘度に対応するＶＢＩ値は、ＭＡＴＬＡＢ（登録商標）（マスワークス社（Ｍａｔｈｗｏｒｋｓ，Ｉｎｃ．））におけるｆｍｉｎｂｎｄ関数のような、標準スカラー線形関数最小化ルーチンを用いて、〔１０〕から見出される。
【００８５】
検査データの重み付け
上記のアルゴリズムにおいてステップ２および３で使用される検査データは、特許文献１に記載のように重み付けされる。特に、重量Ｗは〔１７〕の形を有する。
【数５２】
2008510877.tif 000058
【００８６】
Ｒは解析された未知の材料に関するレベルで計算された検査データの再現性である。εはＸにおける修正された基準スペクトルの、ポイント当たりの平均分散である。０．２〜０．２５ｍｍセル内で収集された原油のスペクトルに関して、εは０．００５と想定できる。αは調整可能なパラメータである。αはステップ２および３からの検査データの予測に関する、所望の誤差分布を得るために選ばれる。
【００８７】
粘度データの大きさは温度と共に変化するため、ステップ３または代替ステップ２における適合への、その貢献度もまた変化するであろう。従って、重み付けのための調整可能なパラメータは、粘度データを異なる温度において用いる時には、比較出来る結果を得るように調整されなければならない。二つ以上の検査が適合の中に含まれる時の検査データ間の相互作用のため、全ての重み付けは粘度測定温度Ｔに依存するであろう。
【数５３】
2008510877.tif 000059
【００８８】
αの値は各々の基準原油がＸから取り出され、未知の原油ｘｕとして扱われる、交差検証解析を用いて各粘度測定温度において決定される。
【００８９】
予測品質
異なる検査入力、または異なる基準の集合を用いた予測は異なるであろう。検査データはそれが幾つかの分析データの予測を改善する場合にのみ、解析の中に含まれる。しかしながら、異なる検査入力および／または異なる基準の集合を用いて行なわれた予測の品質を比較できることは、有用である。実験室における適用に関して、そのような比較は検査データの品質のチェックとして使用可能である。オンラインでの適用に関して、検査データを生み出すために用いられるアナライザは、一時的に故障またはメンテナンスのために利用不可能になることがあり、検査データの欠如がどのように予測の品質に影響するかを知ることが望ましい。
【００９０】
検査データの異なる部分集合を用いて行なわれる予測を比較する目的で、該予測データの全体の品質を表わす単一の品質パラメータを有することが望ましい。予測可能な多数の分析特性を前提とすると、全ての可能な予測の品質を表わすことは実際的でない。しかしながら鍵となる特性の集合に関しては、単一の品質パラメータが定義可能である。
【００９１】
適合品質（ＦＱ）は〔１９〕により定義される。
【数５４】
2008510877.tif 000060
ｆ（ｃ，ｆ，ｉ）は適合におけるゼロ以外の係数の数ｃ、スペクトルポイントの数ｆ、および使用される検査の数ｉの関数である。原油分析データの予測への本発明の適用のために、一つの妥当な関数は次の形で見出される。
【数５５】
2008510877.tif 000061
【００９２】
ε指数は望ましくは０．２５のオーダーである。ＦＱは該計算の中の各ステップにおいてＲ２の値から計算される。適合品質カットオフ（ＦＱＣＩＲ）は該計算のステップ１からの結果、即ちＦＴ−ＩＲスペクトルのみに基づく解析のために定義される。ＦＱＣＩＲは幾つかの最小性能基準に基づいて選択される。適合品質比率は次に〔１６〕により定義される。
【数５６】
2008510877.tif 000062
【００９３】
該アルゴリズムにおけるステップ２および３に関して、ＦＱＣＩＲ，ＡＰＩカットオフおよびＦＱＣＩＲ，ＡＰＩ，Ｖｉｓｃカットオフもまた定義される。これらのカットオフは、異なる入力を用いてなされる予測の精度に、出来る限り近く合致するように設計された最適化手順により決定される。該カットオフはＦＱＲＩＲ，ＡＰＩおよびＦＱＲＩＲ，ＡＰＩ，Ｖｉｓｃを定義するために使用される。
【００９４】
これらのＦＱＲ値は、比較されるべき異なる検査入力および異なる基準の小集合を用いて、解析がなされることを可能にする所望の品質パラメータである。一般に、より低いＦＱＲ値を生み出す解析は、一般により正確な予測を生み出すと期待できる。同様に、同じＦＱＲの適合を生み出す異なる検査入力、または異なる基準小集合を用いてなされる二つの解析は、類似した精度の分析予測を生じることが期待される。
【００９５】
ＦＱＣＩＲ，ＡＰＩおよびＦＱＣＩＲ，ＡＰＩ，Ｖｉｓｃの値はまた性能基準に基づいて設定される。分析特性の限界的な集合が選択される。ステップ２（ＦＴ−ＩＲおよびＡＰＩ比重）およびステップ３（ＦＴ−ＩＲ、ＡＰＩ比重、および粘度）からの分析予測に関して、１以下のＦＱＲ値を有するサンプルに関する予測が、ステップ１（ＦＴ−ＩＲのみ）から得られるものと同程度であるようにＦＱＣ値が選ばれる。検査のための重み付けは検査に関する予測誤差が、それらの試験方法に関して予期される誤差と合致するように、同時に調整される。ＦＱＣ値および検査の重み付けは、標準最適化手順を用いて調整されることができる。
【００９６】
１以下のＦＱＲ値を生み出す解析は、ティア１フィットと呼ばれる。１よりも大きく、１．５以下のＦＱＲ値を生み出す解析はティア２フィットと呼ばれる。
【００９７】
信頼区間
特定の分析予測がプロセス適用に関して妥当かどうかを決定する際に、該予測の不確実性の推定を用意することは有用である。信頼区間は、未知の材料に関する予測特性と、未知の材料が基準の解析に従った場合に得られたであろう値との間の、予期される一致を表わす。各特性に関する信頼区間は、ＦＱＲの関数として推定される。
信頼区間の一般的な形は：
【数５７】
2008510877.tif 000063
で表わされる。
ｆ（Ｅｒｅｆ）は基準特性測定における誤差の関数である。ｔは選択された確率レベルおよびＣＩ計算における自由度の数に関するｔ統計値である。ｓは一旦ＦＱＲおよび基準特性誤差依存度が取り除かれた場合の、予測残差の標準偏差である。
【００９８】
原油分析データの予測への本発明の適用に関して、信頼区間の次のような形が予測誤差の有用な推定を提供するために見出されている。
絶対誤差ＣＩ：
【数５８】
2008510877.tif 000064
相対誤差ＣＩ：
【数５９】
2008510877.tif 000065
ａおよびｂは、基準データの交差検証解析の間に得られる誤差分布に適合するために計算されたパラメータである。ｙは測定された分析特性であり、そして
【数６０】
2008510877.tif 000066
は対応する予測特性である。いずれのＣＩが適用されるかは、参照方法の誤差特性に依存する。参照方法誤差が特性レベルから独立していると予想される特性データに対しては絶対誤差ＣＩが使用され、パラメータｂはゼロである。参照方法誤差が特性レベルに直接比例すると予想される特性データに対しては相対誤差ＣＩが使用される。参照方法誤差が特性レベルに依存するが、直接比例しないと予想される特性データに対しては絶対誤差ＣＩが使用され、ａとｂの両方はゼロでないことが可能である。
【００９９】
該適合に含まれる検査データに関して、信頼区間は僅かに異なる形をとる。
検査に関する絶対誤差ＣＩ：
【数６１】
2008510877.tif 000067
検査に関する相対誤差ＣＩ：
【数６２】
2008510877.tif 000068
方程式〔２５〕は参照方法誤差が特性レベルから独立しているＡＰＩ比重のような検査に適用される。方程式〔２６〕は参照方法誤差が特性レベルに直接比例する粘度のような検査に適用される。
【０１００】
基準部分集合を用いた解析
現在の発明が原油分析データの予測のため原油の解析に適用されるとき、結果として生じる適合および予測の品質が妥当である限り、解析において使用される基準を、解析されている未知の原油と最も類似の原油に限定することが望ましい。様々なサイズの部分集合が未知の原油に対するそれらの類似性に基づいて試験されることができる。原油に関して、次の部分集合の定義が有用であることが見出されている。
【０１０１】
【表５】
2008510877.tif 000069
【０１０２】
未知の原油の解析中に、より小さな部分集合を用いてティア１フィットが得られた場合、次の利点が実現される。
・解析により生み出された仮想混合はより少ない構成要素を有し、分析特性データの計算を単純化および迅速化するであろう。
・多変量解析または検査測定により直接検知されない痕跡レベルの構成要素に関する分析予測は改善され得る。
・解析はエンドユーザー（精製者）がより慣れている可能性のある、原油の仮想混合に基づく。
【０１０３】
部分集合は地理的情報の代わりに地球化学的情報にも基づく可能性がある。プロセスの流れへの適用に対し、部分集合はサンプルのプロセス履歴に基づき得る。
【０１０４】
解析されたサンプルが混合物の場合、部分集合は該混合物において予想される原油の構成要素として、等級、場所、および地域のサンプルから成ることがある。
【０１０５】
混入物質
解析において用いられる基準は、解析されているサンプル中に観察されることがある共通の混入物質を含み得る。一般に、そのような混入物質は未知の材料中に通常存在すると予想されない材料であり、それらは多変量解析測定により検出でき、識別できる。アセトンは恐らく原油サンプリング容器の汚染に起因する、幾らかの原油のＦＴ−ＩＲスペクトルにおいて観察された混入物質の一例である。
【０１０６】
混入物質に関する基準スペクトルは一般的に差により生成される。或る原油サンプルは故意に汚染される。汚染されていない原油のスペクトルは、混入物質のスペクトルを生成するために、故意に汚染されたサンプルのスペクトルから差し引かれる。次に差のスペクトルは純粋な材料を表わすために拡大される。例えば、混入物質が０．１％加えられる場合、差のスペクトルは１０００倍拡大されるであろう。
【０１０７】
混入物質は、ティア１フィットが原油だけを基準として用いて得られなかった時のみ、解析における基準として試験される。ティア１フィットが混入物質なしで得られなかった時に、基準としての混入物質の含有がティア１フィットを生じる場合、そのサンプルは汚染されていると想定される。
【０１０８】
検査データは仮想混合に関し、混入物質を含めて、また除外して計算される。計算された検査データにおける変化が、検査測定方法の再現性の二分の一よりも大きい場合、サンプルは正確に解析するには余りにも汚染されていると考えられる。計算された検査データにおける変化が、検査測定方法の再現性の二分の一よりも小さい場合、混入物質のない仮想混合に基づく分析結果は、サンプルの正確な表示だと想定される。
【０１０９】
代わりに、最大限許容できる汚染レベルは、代表的な原油サンプルに関して上記の判定基準に基づき設定可能である。計算された汚染レベルが、この最大許容可能なレベルを越えた場合、サンプルは正確に解析するには余りにも汚染されていると考えられる。原油内のアセトンに対し、ミディアムＡＰＩ原油に関して算定された４〜５％の粘度変化に基づいて、０．２５％レベルの最大許容汚染レベルを用いることができる。
【０１１０】
基準として使用される各混入物質に対して、最大許容レベルが設定される。混入物質の計算されたレベルが許容レベルよりも低い場合、分析予測は依然として継続することができ、そして不確実性が適合品質比率に基づき推定される。この最大許容レベルを超えると、分析予測は混入物質の存在のために正確さが減少し得る。
【０１１１】
複数の混入物質が基準として使用される場合、最大の組合せレベルが設定されてもよい。組合せ汚染レベルが最大の組合せレベルよりも低い場合、分析予測は依然として継続することができ、そして不確実性が適合品質比率に基づき推定される。この最大組合せレベルを超えると、分析予測は混入物質の存在のために正確さが減少し得る。
【０１１２】
解析図式
〔１９〕における関数ｆ（ｃ，ｆ，ｉ）が１に近い（例えば、〔２０〕のεの値がゼロに近い）場合、より多くの構成要素が混合物に加えられるに従ってＦＱは減少する傾向となり、基準のより大きな部分集合についてなされる解析は、より小さいＦＱ値を生じる傾向になるであろう。この場合、原油分析データの予測に対する本発明の適用に関して、図１に表わされた「第一のティア１フィット」図式は、妥当な予測品質を生み出すことが見出されている。単純化のため、ＦＴ−ＩＲのみ、ＦＴ−ＩＲとＡＰＩ、またはＦＴ−ＩＲ、ＡＰＩ、および粘度に基づく解析のみが示されている。ＦＴ−ＩＲおよび粘度に関する解析もまた用いられたならば、別の列が図中の図式に加えられたであろう。
【０１１３】
未知の材料に関するＡＰＩ比重および粘度が測定されたと仮定すると、解析図式はポイント１において出発する。ユーザーは解析において使用されるべき基準の特定の集合を供給してもよい。適合はここで上記に説明された三つのステップに従って行なわれる。ＦＴ−ＩＲのみに基づく適合（ステップ１）、およびＦＴ−ＩＲとＡＰＩに基づく適合（ステップ２）が計算されるが、それらはこのポイントでは評価されない。ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく適合がティア１フィットを生じる場合、解析は完了し結果が報告される。
【０１１４】
ポイント１での解析がティア１フィットを生じない場合、次に該プロセスはポイント２に進む。基準セットは最初に選択された原油と同じ原油等級の、全ての基準を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告される。
【０１１５】
ポイント２での解析がティア１フィットを生じない場合、次に該プロセスはポイント３に進む。基準セットは最初に選択された原油と同じ場所から来た全ての基準を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告される。
【０１１６】
ポイント３での解析がティア１フィットを生じない場合、次に該プロセスはポイント４に進む。基準セットは最初に選択された原油と同じ地域から来た全ての基準を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告される。
【０１１７】
ポイント４での解析がティア１フィットを生じない場合、次に該プロセスはポイント５に進む。基準セットは全ての基準原油を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告される。
【０１１８】
ポイント５での解析がティア１フィットを生じない場合、次に該プロセスはポイント６に進む。基準セットは全ての基準原油および混入物質を含むように拡張される。三ステップの解析は再び行なわれ、ＦＴ−ＩＲ、ＡＰＩ比重、および粘度に基づく解析が調査される。この解析がティア１フィットを生じる場合、解析は完了し結果が報告され、そして該サンプルは品質が落ちていると報告される。混入物質が最大許容水準を超えない場合、分析結果は依然として計算されてもよく、信頼区間（Ｃｏｎｆｉｄｅｎｃｅ  Ｉｎｔｅｒｖａｌｓ）が適合ＦＱＲに基づき見積もられる。混入物質が許容水準を超えた場合、その結果はＦＱＲにより示されたものよりも精度が低くなり得る。
【０１１９】
ポイント６での解析がティア１フィットを生じない場合、次にＦＴ−ＩＲおよび（各ポイントにおけるステップ２からの）ＡＰＩ比重に基づく適合は、これらのいずれかがティア１フィットを生じるかどうか決定するために調査される。選択された基準に関する適合は最初に調査される（ポイント７）。この解析がティア１フィットを生じた場合、解析は完了し結果が報告される。そうでない場合、該プロセスはポイント８に続き、ＦＴ−ＩＲおよびＡＰＩ比重を用いて選択された原油と同じ等級の原油に基づく適合が調査される。該プロセスはポイント９（同じ場所の原油）、ポイント１０（同じ地域の原油）、ポイント１１（全ての原油）、およびポイント１２（全ての原油および混入物質）に関して適合のチェックを続け、ティア１フィットが見出された場合は中止され、さもなければ継続される。ティア１フィットがＦＴ−ＩＲおよびＡＰＩ比重を用いて見出されない場合、（各ポイントにおけるステップ１からの）ＦＴ−ＩＲのみの適合が、ポイント１３（選択された基準）、ポイント１４（同じ等級）、ポイント１５（同じ場所）、ポイント１６（同じ地域）、ポイント１７（全ての原油）、およびポイント１８（全ての原油および混入物質）に関する適合をチェックしながら調査され、ティア１フィットが見出された場合は中止され、さもなければ継続される。
【０１２０】
ティア１フィットが見出されない場合、最高のＦＱＲ値を生み出す解析が選択され、報告される。ＦＱＲ値が１．５以下の場合、その結果はティア２フィットとして報告される。さもなければ、それは適合不成功として報告される。
【０１２１】
粘度のデータが得られない場合、該解析図式はポイント７において出発し、上記に論じられているように続けられるであろう。粘度もＡＰＩ比重も得られなかった場合、該解析図式はポイント１５において出発し、上記に論じられているように続けられるであろう。
【０１２２】
〔１９〕における関数ｆ（ｃ，ｆ，ｉ）が１に近くない（例えば、〔２０〕のεの値が例えば０．２５である）場合、より多くの構成要素が混合物に加えられるに従ってＦＱは必ずしも減少せず、基準のより大きな部分集合についてなされる解析は、より小さいＦＱ値を生じない可能性がある。この場合、原油分析データの予測に対する本発明の適用に関して、「最良適合」図式がより妥当な予測品質を生み出し得る。
【０１２３】
ＡＰＩ比重および粘度データの両方が得られる場合、図１における列１の解析１〜６は評価され、最小のＦＱＲを生じる解析が最良適合として選択される。最良適合に対するＦＱＲが１未満の場合、解析は完了し結果が報告される。
【０１２４】
ＡＰＩ比重および粘度を用いて得られた最良適合がティア１フィットでない場合、図１における列２の解析７〜１２が評価され、最小のＦＱＲを生じる解析が最良適合として選択される。最良適合に対するＦＱＲが１未満の場合、解析は完了し結果が報告される。
【０１２５】
ＡＰＩ比重を用いて得られた最良適合がティア１フィットでない場合、図１における列３の解析１３〜１８が評価され、最小のＦＱＲを生じる解析が最良適合として選択される。最良適合に対するＦＱＲが１未満の場合、解析は完了し結果が報告される。
【０１２６】
いずれの解析もティア１フィットを生じない場合、最小のＦＱＲを生じる解析が最良適合として選択され、報告される。ＦＱＲが１．５未満の場合、その結果はティア２フィットとして報告され、そうでなければ不適合として報告される。
【０１２７】
ライブラリー交差検証
基準ライブラリーの性能を評価し最適化するために、交差検証法が用いられる。一つの反復法において、一つの基準がライブラリーから除去され、あたかもそれが未知であるかのように解析される。該基準は次にライブラリーに戻される。この手順は各基準が除外され、そして一度解析されるまで繰り返される。
【０１２８】
交差検証法は解析図式内のいずれのポイントでもシミュレートされることができる。従って例えば、交差検証は検査入力としてＡＰＩ比重および粘度の両方を用いて行なうことができ、そして基準が除外されているのと同じ場所からの基準のみを用いて行なうことができる（ポイント３のシミュレーション）。
【０１２９】
基準ライブラリー最適化
与えられたＦＱＲに関する解析が、検査入力または基準の部分集合選定にかかわらず同程度の分析予測を生じるために、ＦＱＣ値および検査の重み付けを注意深く最適化することが必要である。この最適化は次の方法で完遂できる。
【０１３０】
ＦＴ−ＩＲのみの解析に関して：
Ｉ．  最小限の性能が設定される。
ＩＩ．  ＦＴ−ＩＲのみを用いて行なわれる解析に対し、解析図式中のポイント１３〜１７をシミュレートするために交差検証解析が行なわれる。これらのポイントに対する結果は組み合わされ、適合品質（ＦＱ）が各々の結果に対して計算される。選択された分析特性は各適合に基づいて予測される。
ＩＩＩ．  該結果は適合品質（ＦＱ）を増すために選別される。
ＩＶ．  順々に、各ＦＱ値は仮のＦＱＣとして選択され、そして仮のＦＱＣ値が計算される。各原油に関して、いずれのポイント（１３〜１７）で解析が終了したであろうかについて決定が行なわれる。これらの停止ポイントに対応する結果が収集され、そして分析予測のための統計値が計算される。これらの結果はこの仮のＦＱＣに関する反復結果と呼ばれる。
Ｖ．  最小限の性能基準に合う最大のＦＱ値がＦＱＣＩＲとして選択される。
ＶＩ．  ステップＩＶからの反復結果は、示されたＦＱＣについての解析から得られたであろう結果の代表である。
ＦＴ−ＩＲおよび検査を用いた解析に関して：
ＶＩＩ．  予測がＦＴ−ＩＲだけの解析からのものと合致すべき、分析特性の集合が選択される。
ＶＩＩＩ．  検査データへの適合の基準が設定される。
ＩＸ．  検査の重み付けのために初期評価がなされる。
Ｘ．  ポイント１〜５または７〜１１をシミュレートするために、交差検証解析が行なわれる。これらのポイントに対する結果は組み合わされ、適合品質（ＦＱ）が各々の結果に対して計算される。選択された分析特性は各適合に基づいて予測される。
ＸＩ．  該結果は適合品質（ＦＱ）を増すために選別される。
ＸＩＩ．  順々に、各ＦＱ値は仮のＦＱＣとして選択され、そして仮のＦＱＣ値が計算される。各原油に関して、いずれのポイント（１〜５または７〜１１）で解析が終了したであろうかについて決定が行なわれる。これらの停止ポイントに対応する結果が収集され、そして分析予測のための統計値が計算される。これらの結果はこの仮のＦＱＣに関する反復結果と呼ばれる。
ＸＩＩＩ．  ＦＴ−ＩＲおよび検査を用いてなされた分析予測に関する統計は、ＦＴ−ＩＲのみに基づくものと比較される。それに対する予測が同等である最大のＦＱ値は仮のＦＱＣＩＲ，ＡＰＩまたはＦＱＣＩＲ，ＡＰＩ，ｖｉｓｃとして選択される。
ＸＩＶ．  検査データへの適合は統計的に調査され、確立された基準と比較される。統計値が確立された基準に合致する場合、仮のＦＱＣＩＲ，ＡＰＩまたはＦＱＣＩＲ，ＡＰＩ，ｖｉｓｃの値は受け入れられる。そうでなければ、検査の重み付けは調整され、９〜１３が繰り返される。
ＸＶ．  ステップＸＩＩからの反復結果は、示されたＦＱＣおよび検査の重み付けについての解析から得られたであろう結果の代表である。
【０１３１】
様々な統計的手段がライブラリーの性能を評価するため、および検査への適合を評価するために使用できる。これらは以下を含むが、それに限定されない。
・ティア１フィットに関する分析特性の予測のための交差検証の標準誤差。ｔ（ｐ，ｎ）は確率レベルｐおよび自由度ｎに関するｔ統計値である。その合計はティア１フィットを生じるｎのサンプルにわたって計算される。
【数６３】
2008510877.tif 000070
・ＦＱＲ＝１での信頼区間。
・それに対する予測と測定された特性の間の差が該測定の再現性より小さい、ティア１フィットに関する予測の割合。
【０１３２】
基準原油が混入物質を含まないため、ステップ６、１２、および１８に対する適合はライブラリー最適化に含まれないことに注意されたい。
【０１３３】
信頼区間の計算
適合に含まれる検査に関し、信頼区間（ＣＩ）はＦＱＲに関してのみ定義される。次の手順は、含まれる検査に関する信頼区間の計算に用いられる。
・検査（例えばＡＰＩ比重）に関する絶対誤差ＣＩ。
―  上記のステップＸＶからのｎ個の反復結果の各々について、適合から予測される検査と、入力（測定された）検査値
【数６４】
2008510877.tif 000071
との間の差を計算する。
―  ｄｉを
【数６５】
2008510877.tif 000072
で割る。
―  これらの測られた結果の二乗平均平方根、
【数６６】
2008510877.tif 000073
を計算する。
―  所望の確率レベルおよび自由度ｎに対するｔ値を計算する。
―  信頼区間はそのとき方程式〔２５〕で与えられる。
・検査（例えば粘度）に関する相対誤差ＣＩ。
―  上記のステップＸＶからのｎ個の反復結果の各々について、適合から予測される検査と、入力（測定された）検査値
【数６７】
2008510877.tif 000074
との間の相対的差異を計算する。
―  ｒｉを
【数６８】
2008510877.tif 000075
で割る。
―  これらの測られた結果の二乗平均平方根、
【数６９】
2008510877.tif 000076
を計算する。
―  所望の確率レベルおよび自由度ｎに対するｔ値を計算する。
―  信頼区間はそのとき方程式〔２６〕で与えられる。
・分析予測に関する絶対誤差
―  パラメータａおよびｂの推定は交差検証解析からの全ての結果（ポイント１〜５、ポイント７〜１１、またはポイント１３〜１７）を用いてなされる。
―  交差検証解析からのｍ個の結果の各々について、分析特性値の予測と測定の間の差ｄｉ、即ち
【数７０】
2008510877.tif 000077
を計算する。
―  ａおよびｂの初期推定に関して、ｍ個の結果の各々について
【数７１】
2008510877.tif 000078
を計算する。
―  各々の結果について、比率
【数７２】
2008510877.tif 000079
を計算する。
―  ｍ比率の分布に関して、分布の正規性の尺度である統計値を計算する。そのような統計値はアンダーソン・ダーリング（Ａｎｄｅｒｓｏｎ−Ｄａｒｌｉｎｇ）統計、リリーフォース（Ｌｉｌｌｉｅｆｏｒｓ）統計、ジャルク・ベラ（Ｊａｒｑｕｅ−Ｂｅｒａ）統計、またはコルモゴロフ・スミルノフ（Ｋｏｌｍｏｇｏｒｏｖ−Ｓｍｉｒｎｏｖ）統計を含むが、それらに限定されない。ａおよびｂの値は計算された正規性統計に基づく分布の正規性を最大化するために調整される。アンダーソン・ダーリング統計に関して、これは統計値を最小化するためのａおよびｂの調整を含む。
―  ｎ個の反復結果の各々について、分析特性値の予測と測定の間の差ｄｉ、即ち
【数７３】
2008510877.tif 000080
を計算する。
―  上記で決定されたａおよびｂの値を用いて、ｎ個の反復結果の各々について
【数７４】
2008510877.tif 000081
を計算する。
―  測られた差異の二乗平均平方根、
【数７５】
2008510877.tif 000082
を計算する。
―  所望の確率レベルおよび自由度ｎに対するｔ統計値を計算する。
―  信頼区間はそのとき方程式〔２３〕で与えられる。
―  基準特性測定の再現性がレベルと無関係である場合、パラメータｂはゼロに設定されてもよく、パラメータａのみが調整される。
―  別の、更に複雑な表現がｆ（Ｅｒｅｆ）の代わりに置き換わって、上記と同じやり方で最適化され得る。例えば、公開された再現性による方法に関して、ｆ（Ｅｒｅｆ）は公開された再現性と同じ関数形式で表わし得る。
・分析予測に関する相対誤差
―  パラメータａの推定は、交差検証解析からの全ての結果（ポイント１〜５、ポイント７〜１１、またはポイント１３〜１７）を用いてなされる。
―  交差検証解析からのｍ個の結果の各々について、分析特性値の予測と測定の間の相対的差異ｒｉ、即ち
【数７６】
2008510877.tif 000083
を計算する。
―  ａおよびｂの初期推定に関して、ｍ個の結果の各々について
【数７７】
2008510877.tif 000084
を計算する。
―  各々の結果について、比率
【数７８】
2008510877.tif 000085
を計算する。
―  ｍ比率の分布に関して、分布の正規性の尺度である統計値を計算する。そのような統計値はアンダーソン・ダーリング（Ａｎｄｅｒｓｏｎ−Ｄａｒｌｉｎｇ）統計、リリーフォース（Ｌｉｌｌｉｅｆｏｒｓ）統計、ジャルク・ベラ（Ｊａｒｑｕｅ−Ｂｅｒａ）統計、またはコルモゴロフ・スミルノフ（Ｋｏｌｍｏｇｏｒｏｖ−Ｓｍｉｒｎｏｖ）統計を含むが、それらに限定されない。ａおよびｂの値は計算された正規性統計に基づく分布の正規性を最大化するために調整される。アンダーソン・ダーリング統計に関して、これは統計値を最小化するためのａおよびｂの調整を含む。
―  ｎ個の反復結果の各々について、分析特性値の予測と測定の間の相対的差異ｒｉ、即ち
【数７９】
2008510877.tif 000086
を計算する。
―  上記で決定されたａおよびｂの値を用いて、ｎ個の反復結果の各々について
【数８０】
2008510877.tif 000087
を計算する。
―  測られた差異の二乗平均平方根、
【数８１】
2008510877.tif 000088
を計算する。
―  所望の確率レベルおよび自由度ｎに対するｔ統計値を計算する。
―  信頼区間はそのとき方程式〔２３〕で与えられる。
―  基準特性測定の再現性がレベルと無関係である場合、パラメータｂはゼロに設定されてもよく、パラメータａのみが調整される。
―  別の、更に複雑な表現がｆ（Ｅｒｅｆ）の代わりに置き換わって、上記と同じやり方で最適化され得る。例えば、公開された再現性による方法に関して、ｆ（Ｅｒｅｆ）は公開された再現性と同じ関数形式で表わし得る。
【図面の簡単な説明】
【０１３４】
【図１】予測の原油分析データの概要を示す。
【図１】
2008510877.tif 000089
【国際調査報告】
2008510877.tif 000090