JPH08502360A - 軽油のセタン価の予測方法 - Google Patents
軽油のセタン価の予測方法Info
- Publication number
- JPH08502360A JPH08502360A JP6509633A JP50963394A JPH08502360A JP H08502360 A JPH08502360 A JP H08502360A JP 6509633 A JP6509633 A JP 6509633A JP 50963394 A JP50963394 A JP 50963394A JP H08502360 A JPH08502360 A JP H08502360A
- Authority
- JP
- Japan
- Prior art keywords
- data
- cetane number
- network
- cetane
- gas oil
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- CIWBSHSKHKDKBQ-JLAZNSOCSA-N Ascorbic acid Chemical compound OC[C@H](O)[C@H]1OC(=O)C(O)=C1O CIWBSHSKHKDKBQ-JLAZNSOCSA-N 0.000 title claims abstract description 110
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000001228 spectrum Methods 0.000 claims abstract description 39
- 238000013528 artificial neural network Methods 0.000 claims abstract description 37
- 230000003595 spectral effect Effects 0.000 claims abstract description 28
- 238000005259 measurement Methods 0.000 claims abstract description 24
- 238000010521 absorption reaction Methods 0.000 claims abstract description 20
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 15
- 238000007619 statistical method Methods 0.000 claims abstract description 8
- 239000003921 oil Substances 0.000 claims description 69
- 230000008569 process Effects 0.000 claims description 12
- 238000012417 linear regression Methods 0.000 claims description 5
- 238000000513 principal component analysis Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010561 standard procedure Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 abstract description 26
- 238000002835 absorbance Methods 0.000 abstract 1
- 239000000446 fuel Substances 0.000 description 21
- QPUYECUOLPXSFR-UHFFFAOYSA-N 1-methylnaphthalene Chemical compound C1=CC=C2C(C)=CC=CC2=C1 QPUYECUOLPXSFR-UHFFFAOYSA-N 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 239000002283 diesel fuel Substances 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- TVMXDCGIABBOFY-UHFFFAOYSA-N octane Chemical compound CCCCCCCC TVMXDCGIABBOFY-UHFFFAOYSA-N 0.000 description 3
- 238000007405 data analysis Methods 0.000 description 2
- 238000004821 distillation Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- DCAYPVUWAIABOU-UHFFFAOYSA-N hexadecane Chemical compound CCCCCCCCCCCCCCCC DCAYPVUWAIABOU-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 238000003909 pattern recognition Methods 0.000 description 2
- 238000004886 process control Methods 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 241000282326 Felis catus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 238000012864 cross contamination Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000275 quality assurance Methods 0.000 description 1
- 238000012113 quantitative test Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/26—Oils; Viscous liquids; Paints; Inks
- G01N33/28—Oils, i.e. hydrocarbon liquids
- G01N33/2829—Mixtures of fuels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3577—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Spectrometry And Color Measurement (AREA)
- Sampling And Sample Adjustment (AREA)
- Investigating Or Analysing Materials By The Use Of Chemical Reactions (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
- Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
Abstract
(57)【要約】
軽油のセタン価を予測する方法であって、a)多様な給源から得た多数の軽油のセットのI.R.スペクトルを測定するステップ;b)前記スペクトル域内で一定範囲の波数を選択し、多数の当該波長を吸収データに変換し、前記吸収データを、多変量統計解析またはニューラルネットワークへの入力として使用するステップ;c)多変量統計解析またはニューラルネットワークを使用してスペクトルデータを分析するステップ;d)従来の測定によって軽油のセタン価を決定するステップ;e)測定I.R.データと従来の測定によるセタン価データとを含むトレーニングデータセットを選択し、得られた吸収値をセタン値と相関させ、予測データセットを生成するステップ;f)前記データを、セタン価が未知の軽油の同じ条件下で得た赤外線スペクトルに適用し、未知の軽油のセタン価を得るステップを実施する方法。
Description
【発明の詳細な説明】
軽油のセタン価の予測方法
本発明は、(近)赤外線((N.)I.R.)スペクトルをセタン価と相関さ
せることにより軽油のセタン価を予測する方法に係わる。当業者には公知のよう
に、有機化合物は赤外線域(約1〜300μm)に固有のスペクトルフィンガー
プリントを有する。
最近の市場では、バルク燃料に対して、製品の品質及び性能により重きが置か
れる傾向がある。全ての販売キャンペーンは、顧客に、その燃料が他よりも著し
く優れた性能を表わすことを納得させるよう図っており、この考え方は続くと予
想される。このような方策には、精製からブレンディング、そして最終的に顧客
に配給されるまでの十分な品質保証による支援が必要である。この条件は、燃料
の性能を測定する品質モニター装置を使用することにより満足されるのが最善で
ある。
今のところ、製油所を離れた後、燃料がまだ規格範囲内にあるかどうかを見る
定量検査が継続的に行われることはなく、透明度や臭い、表示された密度を有す
るかどうかは、手作業で検査するしかない。最近の進歩した燃料の高度に
複雑な特性を考えると、これは十分ではない。燃料の性能規格を実際に測定する
品質モニター装置パッケージは有益であり、そうすれば、燃料の劣化、他の燃料
との相互汚染(cross contamination)、燃料組成の変化を流通経路において検
出し、指摘し得る。このことにより、顧客が所望のものを確実に入手するだけで
なく、流通経路で起こり得る問題を発見し、対処し得る。
軽油などの燃料の重要な性能尺度として、ディーゼル油の着火性を表わす量で
あるセタン価がある。セタン価測定用のエンジンが存在する。
自動車用ガソリンのオクタン価測定と同様に標準燃料が使用される。かかる標
準燃料はn−セタン及びα−メチルナフタレンである。着火性は、試験エンジン
において着火性能が燃料のそれと一致するα−メチルナフタレンとのブレンド中
のセタンの容量%であるセタン価で表わされる。実際にはセタンエンジンはめっ
たに使用されず、セタンは通常は他の測定値、例えば密度及び蒸留性状から「セ
タン指数」として計算される。しかしながらセタン指数は、将来のディーゼル燃
料には適当でない可能性があり、セタン価に代わる測定値が必要とされている。
また、軽油のセタン価を近赤外線(N.I.R.)スペクトルの関数として決
定することが既に提案されている(例えば欧州特許出願第0,304,232号
明細書参照)。
しかしながら、これまでのところ(近赤外線及び中赤外線を含む)I.R.ス
ペクトル(有利には0.78〜30μmの波長を適用する)から軽油のセタン価
を予測することはできなかった。
従って本発明の目的は、未知の軽油のセタン価を予測する方法であって、広範
囲のセタン価を有する混合生成物及びプロセスストリーム(process stream)の
セタン価を予測する方法を提供することである。
本発明は、軽油のセタン価を予測する方法であって、
a)多様な給源から得た多数の軽油のI.R.スペクトルを測定するステップ
;
b)前記スペクトル域内で一定範囲の波数を選択し、多数の該当波長を吸収デ
ータに変換し、前記吸収データを、多変量統計解析またはニューラル(neural)
ネットワークへの入力として使用するステップ;
c)多変量統計解析またはニューラルネットワークを使用してスペクトルデー
タを分析するステップ;
d)従来の測定によって軽油のセタン価を決定するステップ;
e)測定されたI.R.データと従来の測定によるセタン価データとを含むト
レーニングデータセットを選択し、得られた吸収値をセタン価と相関させ、予測
データセットを生成するステップ;及び
f)前記データを、セタン価が未知の軽油の同じ条件下で得た赤外線スペクト
ルに適用し、未知の軽油のセタン価を得るステップ
からなる方法を提供する。
前述したように、本発明は、軽油のセタン価をその赤外線スペクトルと相関さ
せるという原理に基づいている。この原理自体は公知であり(例えば欧州特許出
願第0,304,232号明細書及び欧州特許出願第0,285,251号明細
書)、詳細には記載しない。
本発明によれば、多様な給源から得た多数の軽油(有利には少なくとも100
種)の赤外線スペクトルが測定される。異なった種々の2次変換プロセスによる
成分があっても良い。
このあとの任意の統計的予測ツールの汎用性及び適用性
を決定することになるので、多様性は重要である。
使用するスペクトル域は、最終的な完全ブレンド軽油に対しては9000〜4
500波数である。セタン価に影響するので、セタン着火性向上剤濃度の測定も
必要である(1600〜1700波数、有利には1630波数)。リアルタイム
プロセス制御用の変換プロセスストリームで使用する場合にも、スペクトル域は
9000〜4500波数である。
次いで、従来の試験エンジン測定によりセタン価を決定すると共に、例えば部
分最小二乗法(Partial Least Squares)、多重線形回帰、既約階数回帰(Reduc
ed Rank Regression)、主成分解析などのそれ自体は公知の多変量統計解析、ま
たはニューラルネットワークを使用し、スペクトルを分析する。
次いで、セタン価が未知の軽油の同じ条件下で得た赤外線スペクトルに適用し
得る予測データを生成する。
予測の標準誤差はセタン価で0−4未満が達成された。
ニューラルネットワーク自体の一般理論及び一般操作は当業者には公知である
ので、詳細は記載しない。
一般にニューラルネットワークは、学習期間内に入力変
数と出力変数との相関が探索されるシステムであると定義される。この学習期間
内に十分な例を与えられると、ニューラルネットワークは任意の入力に対して相
関出力を生成することができる。ニューラルネットワークは例えばパターン認識
問題に用途を有する。
当業者には、ニューラルネットワークは、(脳のニューロンと類似の)各々が
重み付けされ且つ他の層のエレメントに接続されている処理エレメントの層から
なることが理解される。ネットワークは、間違いなく適確なデータを用いてトレ
ーニングされる間に、エレメント間の重みを調整することでパターンを学習する
。
有利な学習アルゴリズムに従って、トレーニング誤差、即ち実際の結果と予測
結果との差が、トレーニングパターンからフィードバックを受け取らない隠れ層
に、ネットワークを介して逆伝送される。相互接続の重みは誤差を最小化するよ
うに誤差の方向で僅かずつ調整され、トレーニングデータがもう一度実行される
。誤差が容認可能レベル、通常は初期測定値が再現可能なレベルに達するまでこ
れが何回も繰り返される。
前述したように、多変量統計データ解析法はそれ自体公
知であるが、ここでは明瞭化のため、主成分解析を使用する赤外線スペクトルデ
ータの分析を詳述する。
上記予測において、トレーニングセットは、そのスペクトル及び「エンジン測
定」セタン価を、あとで未知の軽油のセタン価を決定するために使用される予測
データの分析及び生成に使用する軽油と定義される。このトレーニングセットの
選択は、データの汎用性及び予測性能に厳密に影響するが故に明らかに重要であ
る。トレーニングセットの軽油に対応するスペクトルセットに関するデータ解析
は以下のように行う:
1.セットの平均スペクトルを生成し、個々のスペクトルと平均との差を計算
する。
2.平均スペクトルは約5000データポイント程度となるので、1セット1
00の軽油の解析の問題は極めて難しい。取扱い可能な問題変数までデータを削
減する方法が必要とされている。差スペクトルを同時に解析することにより人工
スペクトルセットを計算し、データ変動を人工スペクトルによって説明すること
ができる。ここでは差スペクトルは、軽油セットの平均スペクトルと個々の軽油
のスペクトルとの吸収の算術差と定義される。人工スペクトル
は主成分として公知であり、適正な割合で平均スペクトルに加算されたときは、
個々の軽油スペクトルを再生する。例えば、典型的な主成分数は10であり、各
主成分の適当な因子(contribution)を平均に加算することにより、元の軽油ス
ペクトルを再生することができる。各主成分に帰する因子は主成分スコア(scor
e)として公知であり、スコアは各軽油に対して異なる。
3.軽油スペクトルの情報を、それぞれ5000データポイントから10主成
分スコアに削減する。多重線形回帰のような手法を使用し、これらのスコアをエ
ンジン測定セタン価と相関させ得る。
4.このあとの未知の軽油の分析は、スペクトルを測定し、該スペクトルを、
トレーニングセットの軽油のスペクトルを表わすのに前に使用した主成分によっ
て表わすことを含む。この結果、未知の軽油の主成分スコアセットが得られる。
先の多重線形回帰の相関係数を、未知の軽油スペクトルの主成分スコアに適用し
、未知の軽油のセタン価を得る。
5.ニューラルネットワーク法による分析は、測定する波長の数を物理的に減
らすことによりデータを削減するこ
とから、上記とは異なる。このデータ削減は以下のように行われる:トレーニン
グ軽油セットにおいて主成分解析を行い、各スペクトルデータポイントのセタン
価との相関に対する相対重要性を表わす「固有スペクトル」を生成する。次いで
、スペクトル測定値を典型的には5〜10の個別の波長に単純化する。
1つの波長は、計測ドリフトを補正するための伝送基準として使用することが
有利である。
前記基準によって補正された残りの波長は吸収データに変換される。これは対
数的に行うことができ、データは、各波長に対する所定の区間内で数学的に目盛
りを割り当てられる。いずれかの燃料、より多くはプロセスストリームに予想さ
れる最終値(extreme values)を使用し、試験すべき燃料に対して目盛りを設定
し得る、各波長における容認可能な吸収域を与え得る。吸収値は、ニューラルネ
ットワークへの入力として使用される。ニューラルネットワークはパターン認識
タスクを実行するが、ここでの用途は明らかに、セタン価に関するスペクトルデ
ータのパターンを認識することである。ニューラルネットワークを使用して未知
の軽油のセタン価を予測する前に、ネットワークを
「トレーニング」する必要がある。トレーニング方法は、サイズ及びアーキテク
チャーが該認識問題に適合するように設計されたネットワークに対して「トレー
ニングセット」を選択することである。トレーニングフェーズの間、ニューラル
ネットワークには各軽油の赤外線データとエンジン測定セタン価とが繰返し与え
られ、関係が「学習」される。「学習」方法は、赤外線データを与える都度ニュ
ーラルネットワークによってセタン価が予測され、それをエンジン価と比較し、
差をフィードバックとして使用してニューラルネットワークをその次の予測に対
して補正するというものである。一旦ニューラルネットワークが関係を「学習」
したならば、データセットは別のトレーニングセットと、「学習」フェーズでは
使用されない妥当性検査セットとに分割されるべきである。ネットワークが十分
に「学習」し、予測誤差が元来のエンジンによるセタン価決定における不確実性
に匹敵するまで、このプロセスが継続される。一旦このフェーズが完了したなら
ば、トレーニングの間に使用した波長における赤外線スペクトルデータを測定す
れば、ニューラルネットワークを使用して未知の軽油のセタン価を予測し得る。
特に、入力と既知の出力、即ち軽油の赤外線とその関連セタン価データとを繰
り返し与えて両者の関係を学習させ、標準的方法によって測定された実際の関連
セタン価データに対する予測値の性能をトレーニングデータにおいてモニターし
、吸収値とセタン価とを相関させ、前記学習期間後に調整ネットワークの相互接
続重み及びバイアスの値のセットを生成し、ネットワークアルゴリズムを使用し
、セタン価が未知の軽油の同じ条件下で得た赤外線スペクトルに対して調整値を
適用することにより、前記予測データを生成する。
使用するネットワークは、例えば第1層に4つの入力ノード、入力と出力の間
の第2層に2つの隠れノード、更に第3層に1つの出力ノードを含む3層アーキ
テクチャーを有するのが好ましい。これは(4,2,1)ネットワークと称され
ている。スペクトルデータは入力ノードに入力として与えられ、そこでは生成物
品質情報が出力である。
当業者には公知のように、ノードは所定の相互接続重みを有し、更にはバイア
スされ得る。
ネットワークの重み及びバイアスを記憶し、それを使用して測定された赤外線
吸収値を含む入力データを分析し、
そのパターンを軽油のセタン価と相関させることができる。即ち、赤外線データ
からセタン価を表わすネットワークアルゴリズムを使用する予測において、トレ
ーニングされたあと妥当性検査データセットによって試験された重要なパラメー
ターは、ノード間の相互接続の重みと、隠れノード及び出力ノードにおけるバイ
アスである。
かかるパラメーターは、将来の燃料サンプルのセタン価分析用ネットワークア
ルゴリズムにおいて問合せ及び実行され得る。
複数の出力に対しては、各出力に対してニューラルネットワークアルゴリズム
が実行される。実行は、マイクロプロセッサーチップ上のソフトウェアコードに
よるものであり、従って、容易に再プログラムし得るネットワークパラメーター
の変更に対して融通性がある。
適用されるネットワークアーキテクチャーの各層内に存在するノードの正確な
数、場合によっては実際の層の数は可変であることは当業者には理解されよう。
2〜5層を適用することが有利である。
本発明によれば、入力層のノード数は3〜10、(1つまたはそれ以上の)隠
れ層のノード数は1〜10、出力層
のノード数は1〜3であるのが有利である。特に(3,5,1)、(6,6,3
)及び(6,6,6,3)ネットワークを適用し得る。
以下、本発明の方法を実施例によって説明する。実施例A
種々の製油所から収集した20種のブレンド軽油からなるトレーニングセット
を、赤外線スペクトル域9000〜4500波数内で測定した。「セタンエンジ
ン」を使用するASTMD613に従い、全てのサンプルのセタン価を測定した
。トレーニングセットを上述の方法で分析した。主成分を使用して「人工スペク
トル」によってデータを表わす分析を実施し、最初の3つの主成分スペクトルを
6400〜4800波数域内で得た。この実施例に使用したトレーニングセット
には5つの主成分が必要であり、エンジン測定セタン価と、赤外線予測セタン価
との相関0.94、標準セタン価予測誤差0.4未満という結果が示された。上
記分析に代え、セタン価予測用ニューラルネットワークを使用するならば、まず
、第1のトレーニングセットに対して主成分を分析し、重要な個別の波長を同定
し得る適正スペクトルを生成する。次に使用されるニューラルネット
ワーク方法は、先の分析で同定された波長においてトレーニングセット赤外線デ
ータを測定する。実施例B
赤外線セタン価測定の適用可能範囲を示すため、20〜60の広範囲のセタン
価を有する11の軽油ブレンド成分、例えばlight cat cracked cycle oil、直
留分(straightrun distillate)のトレーニングセットを使用した。赤外線スペ
クトルは9000〜4500波数域で測定した。3つの主成分を使用する予測デ
ータを上述の方法において生成した。セタン価の予測値と測定値の相関係数0.
99が得られた。相互妥当性検査から、未知の燃料における予測性能が優れてい
ることが示された。この実施例のデータセットは、0.9より優れた標準予測誤
差をもつ結果を与えた。この性能は先のサンプルより低く、トレーニングセット
がより小さかったこと、及び燃料のセタン価がより広範囲であったことを反映し
ている。より大きなトレーニングセットで、この数字をサンプルAに匹敵するレ
ベルに向上し得るであろう。
上記説明及び添付の図面から、本発明の種々の変形が当業者には明らかである
。このような変形も請求の範囲内に
含まれるものとする。
【手続補正書】特許法第184条の8
【提出日】1994年9月21日
【補正内容】
自動車用ガソリンのオクタン価測定と同様に標準燃料が使用される。かかる標
準燃料はn−セタン及びα−メチルナフタレンである。着火性は、試験エンジン
において着火性能が燃料のそれと一致するα−メチルナフタレンとのブレンド中
のセタンの容量%であるセタン価で表わされる。実際にはセタンエンジンはめっ
たに使用されず、セタンは通常は他の測定値、例えば密度及び蒸留性状から「セ
タン指数」として計算される。しかしながらセタン指数は、将来のディーゼル燃
料には適当でない可能性があり、セタン価に代わる測定値が必要とされている。
また、軽油のセタン価を近赤外線(N.I.R.)スペクトルの関数として決
定することが既に提案されている(例えば欧州特許出願第0,304,232号
明細書参照)。
しかしながら、これまでのところ(近赤外線及び中赤外線を含む)I.R.ス
ペクトル(有利には0.78〜30μmの波長を適用する)から軽油のセタン価
を予測することはできなかった。
従って本発明の目的は、未知の軽油のセタン価を予測する方法であって、広範
囲のセタン価を有する混合生成物及びプロセスストリーム(process stream)の
セタン価を予
測する方法を提供することである。
本発明は、軽油のセタン価を予測する方法であって、
a)複数の軽油セットのI.R.スペクトルを測定するステップ;
b)前記スペクトル域内で一定範囲の波数を選択し、多数の該当波長を吸収デ
ータに変換し、前記吸収データをニューラル(neural)ネットワークへの入力と
して使用するステップ
を含み、更に、
c)ニューラルネットワークを使用してスペクトルデータを分析するステップ
であって、該ニューラルネットワークの層の数が2〜5であり、層が3つ以上の
場合には、入力層のノード数が3〜10であり、(1つまたはそれ以上の)隠れ
層のノード数が1〜10であり、出力層のノード数が1〜3であるステップ;
d)従来の測定によって軽油のセタン価を決定するステップ;
e)測定されたI.R.データと従来の測定によるセタン価データとを含むト
レーニングデータセットを選択し、得られた吸収値をセタン価と相関させ、予測
データセットを
生成するステップ;及び
f)前記データを、セタン価が未知の軽油の同じ条件下で得た赤外線スペクト
ルに適用し、未知の軽油のセタン価を得るステップ
を含むことを特徴とする方法を提供する。
前述したように、本発明は、軽油のセタン価をその赤外線スペクトルと相関さ
せるという原理に基づいている。この原理自体は公知であり(例えば欧州特許出
願第0,304,232号明細書及び欧州特許出願第0,285,251号明細
書)、詳細には記載しない。技術文献RD327135は、ガソリンのオクタン
価を決定するニューラルネットワークを開示しているが、本発明のセタン価を決
定する特定のニューラルネットワークはこれまでに示されていない。
本発明によれば、多様な給源から得た多数の軽油(有利には少なくとも100
種)の赤外線スペクトルが測定される。異なった種々の2次変換プロセスによる
成分があっても良い。
このあとの任意の統計的予測ツールの汎用性及び適用性を決定することになる
ので、多様性は重要である。
使用するスペクトル域は、最終的な完全ブレンド軽油に対しては9000〜4
500波数である。セタン価に影響するので、セタン着火性向上剤濃度の測定も
必要である(1600〜1700波数、有利には1630波数)。
リアルタイムプロセス制御用の変換プロセスストリームで使用する場合にも、
スペクトル域は9000〜4500波数である。
次いで、従来の試験エンジン測定によりセタン価を決定すると共に、例えば部
分最小二乗法(Partial Least Squares)、多重線形回帰、既約階数回帰(Reduc
ed Rank Regression)、主成分解析などのそれ自体は公知の多変量統計解析、ま
たはニューラルネットワークを使用し、スペクトルを分析する。
次いで、セタン価が未知の軽油の同じ条件下で得た赤外線スペクトルに適用し
得る予測データを生成する。
請求の範囲
1.軽油のセタン価を予測する方法であって、
a)複数の軽油セットのI.R.スペクトルを測定するステップ;
b)前記スペクトル域内で一定範囲の波数を選択し、多数の当該波長を吸収デ
ータに変換し、前記吸収データをニューラルネットワークへの入力として使用す
るステップを含み、更に、
c)ニューラルネットワークを使用してスペクトルデータを分析するステップ
であって、該ニューラルネットワークの層の数が2〜5であり、層が3つ以上の
場合には、入力層のノード数が3〜10であり、(1つまたはそれ以上の)隠れ
層のノード数が1〜10であり、出力層のノード数が1〜3であるステップ;
d)従来の測定によって軽油のセタン価を決定するステップ;
e)測定I.R.データと従来の測定によるセタン価データとを含むトレーニ
ングデータセットを選択し、得られた吸収値をセタン価と相関させ、予測データ
セットを生成するステップ;及び
f)前記データを、セタン価が未知の軽油の同じ条件下で得た赤外線スペクト
ルに適用し、未知の軽油のセタン価を得るステップ
を含むことを特徴とする方法。
2.前記スペクトル域が波長0.78〜30μmである請求項1に記載の方法。
3.前記軽油セットが少なくとも100である請求項1または2に記載の方法。
4.前記軽油セットが、異なった種々の2次変換プロセスによる成分を含む請求
項1から3のいずれか一項に記載の方法。
5.最終的な完全ブレンド軽油に対し、使用するスペクトル域が9000〜45
00波数であり、セタン着火性向上剤濃度を、スペクトル域1600〜1700
波数で測定する請求項1から4のいずれか一項に記載の方法。
6.広範囲のセタン価を有する変換プロセスストリームに対し、前記スペクトル
域が9000〜4500である請求項1から4のいずれか一項に記載の方法。
7.入力と既知の出力、即ち軽油の近赤外線データとそのセタン価とを繰り返し
与えて両者の関係を学習させ、標準
的方法によって測定した実際のセタン価データに対する予測値の性能をトレーニ
ングデータにおいてモニターし、吸収値とセタン価とを相関させるステップ;
前記学習期間後に、調整値としてネットワークの相互接続重みとバイアスから
なる値のセットを生成するステップ;
ネットワークアルゴリズムを使用し、セタン価が未知の軽油の同じ条件下で得
た赤外線スペクトルに前記調整値を適用するステップ
を含む、ニューラルネットワークを適用する請求項1から6のいずれか一項に記
載の方法。
8.前記ネットワークが、4つの入力ノード、2つの隠れノード及び1つの出力
ノードを有する((4,2,1)ネットワーク)請求項1から7のいずれか一項
に記載の方法。
9.前記ネットワークが(3,5,1)ネットワークである請求項1から7のい
ずれか一項に記載の方法。
10.前記ネットワークが(6,6,3)ネットワークである請求項1から7の
いずれか一項に記載の方法。
11.前記ネットワークが(6,6,6,3)ネットワークである請求項1から
7のいずれか一項に記載の方法。
─────────────────────────────────────────────────────
フロントページの続き
(81)指定国 EP(AT,BE,CH,DE,
DK,ES,FR,GB,GR,IE,IT,LU,M
C,NL,PT,SE),OA(BF,BJ,CF,CG
,CI,CM,GA,GN,ML,MR,NE,SN,
TD,TG),AT,AU,BB,BG,BR,BY,
CA,CH,CZ,DE,DK,ES,FI,GB,H
U,JP,KP,KR,KZ,LK,LU,MG,MN
,MW,NL,NO,NZ,PL,PT,RO,RU,
SD,SE,SK,UA,VN
(72)発明者 ギルクリスト,キヤサリン・アン
イギリス国、チエシヤー・シー・エイチ・
1・3・エス・エイチ、チエスター、イン
ス、プール・レーン(番地なし)
(72)発明者 トルチヤード,ジヨン・マイケル
イギリス国、チエシヤー・シー・エイチ・
1・3・エス・エイチ、チエスター、イン
ス、プール・レーン(番地なし)
Claims (1)
- 【特許請求の範囲】 1.軽油のセタン価を予測する方法であって、 a)多様な給源から得た多数の軽油のセットのI.R.スペクトルを測定する ステップ; b)前記スペクトル域内で一定範囲の波数を選択し、多数の当該波長を吸収デ ータに変換し、前記吸収データを、多変量統計解析またはニューラルネットワー クへの入力として使用するステップ; c)多変量統計解析またはニューラルネットワークを使用してスペクトルを分 析するステップ; d)従来の測定によって軽油のセタン価を決定するステップ; e)測定I.R.データと従来の測定によるセタン価データとを含むトレーニ ングデータセットを選択し、得られた吸収値をセタン価と相関させ、予測データ セットを生成するステップ;更に f)前記データを、セタン価が未知の軽油の同じ条件下で得た赤外線スペクト ルに適用し、未知の軽油のセタン価を得るステップ からなる方法。 2.前記スペクトル域が波長0.78〜30μmである請求項1に記載の方法。 3.前記軽油セットが少なくとも100である請求項1または2に記載の方法。 4.前記軽油セットが、異なった種々の2次変換プロセスによる成分を含む請求 項1から3のいずれか一項に記載の方法。 5.最終的な完全ブレンド軽油に対し、使用するスペクトル域が9000〜45 00波数であり、セタン着火性向上剤濃度を、スペクトル域1600〜1700 波数で測定する請求項1から4のいずれか一項に記載の方法。 6.広範囲のセタン価を有する変換プロセスストリームに対し、前記スペクトル 域が9000〜4500である請求項1から4のいずれか一項に記載の方法。 7.入力と既知の出力、即ち軽油の近赤外線データとそのセタン価とを繰り返し 与えて両者の関係を学習させ、標準的方法によって測定した実際のセタン価デー タに対する予測値の性能をトレーニングデータにおいてモニターし、吸収値とセ タン価とを相関させるステップ; 前記学習期間後に、調整値としてネットワークの相互接 続重みとバイアスからなる値のセットを生成するステップ; ネットワークアルゴリズムを使用し、セタン価が未知の軽油の同じ条件下で得 た赤外線スペクトルに前記調整値を適用するステップ を含む、ニューラルネットワークを適用する請求項1から6のいずれか一項に記 載の方法。 8.前記ニューラルネットワークの層の数が2〜5である請求項7に記載の方法 。 9.適用される前記ニューラルネットワークが3層または4層アーキテクチャー を有する請求項1から7のいずれか一項に記載の方法。 10.入力層のノード数が3〜10であり、(1つまたはそれ以上の)隠れ層の ノード数が1〜10であり、出力層のノード数が1〜3である請求項9に記載の 方法。 11.前記ネットワークが、4つの入力ノード、2つの隠れノード及び1つの出 力ノードを有する((4,2,1)ネットワーク)請求項9または10に記載の 方法。 12.前記ネットワークが(3,5,1)ネットワークである請求項9または1 0に記載の方法。 13.前記ネットワークが(6,6,3)ネットワークである請求項9または1 0に記載の方法。 14.前記ネットワークが(6,6,6,3)ネットワークである請求項9また は10に記載の方法。 15.前記予測データセットが、I.R.スペクトルの主成分スコア及びそれら の多重線形回帰係数であり、適用される多変量統計解析が主成分解析である請求 項1から6のいずれか一項に記載の方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP92309405.6 | 1992-10-15 | ||
EP92309405 | 1992-10-15 | ||
PCT/EP1993/002835 WO1994009362A1 (en) | 1992-10-15 | 1993-10-12 | Method for prediction of cetane numbers of gasoils |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH08502360A true JPH08502360A (ja) | 1996-03-12 |
Family
ID=8211520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP6509633A Pending JPH08502360A (ja) | 1992-10-15 | 1993-10-12 | 軽油のセタン価の予測方法 |
Country Status (20)
Country | Link |
---|---|
EP (1) | EP0671003B1 (ja) |
JP (1) | JPH08502360A (ja) |
KR (1) | KR950703735A (ja) |
AT (1) | ATE148561T1 (ja) |
AU (1) | AU677609B2 (ja) |
BR (1) | BR9307250A (ja) |
CA (1) | CA2146985A1 (ja) |
DE (1) | DE69307926T2 (ja) |
DK (1) | DK0671003T3 (ja) |
ES (1) | ES2098063T3 (ja) |
FI (1) | FI951720A0 (ja) |
GR (1) | GR3022931T3 (ja) |
HK (1) | HK91297A (ja) |
MY (1) | MY110065A (ja) |
NO (1) | NO951416L (ja) |
NZ (1) | NZ256694A (ja) |
SG (1) | SG45395A1 (ja) |
TR (1) | TR27871A (ja) |
WO (1) | WO1994009362A1 (ja) |
ZA (1) | ZA937593B (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004132921A (ja) * | 2002-10-15 | 2004-04-30 | Yukihiro Toyoda | 性状予測方法 |
JP2004514902A (ja) * | 2000-11-20 | 2004-05-20 | ザ、プロクター、エンド、ギャンブル、カンパニー | ポリマーの予測方法 |
JP2007183242A (ja) * | 2005-12-29 | 2007-07-19 | Ind Technol Res Inst | 移動式燃料分析装置および該移動式燃料分析装置を用いた燃料の品質を測定する方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5572030A (en) * | 1994-04-22 | 1996-11-05 | Intevep, S.A. | Method for determining parameter of hydrocarbon |
GB2301897B (en) * | 1995-06-08 | 1999-05-26 | Univ Wales Aberystwyth The | Composition analysis |
AR003845A1 (es) * | 1995-10-18 | 1998-09-09 | Shell Int Research | Un método para predecir una propiedad física de un residuo de petróleo crudo, de un petróleo combustible (fuel oil) residual o de un material bituminoso. |
AR003846A1 (es) * | 1995-10-18 | 1998-09-09 | Shell Int Research | Una celula de transmision adecuada para el uso en un dispositivo para medir espectros infrarrojos (proximos) de un material hidrocarbonaceo, un espectrometro que la comprende, uso del mismo, metodos para predecir una propiedad fisica de dicho material con dicha celula o espectrometro, y un procedimiento para preparar una composición de betum usando dicho metodo con dicho espectometro |
DE102007042507B4 (de) | 2007-09-07 | 2019-12-12 | HELLA GmbH & Co. KGaA | Verfahren zur Bestimmung von Inhaltsstoffen von Flüssigkeiten, insbesondere von Motoröl |
CN117589711A (zh) * | 2023-11-21 | 2024-02-23 | 天津市产品质量监督检测技术研究院 | 一种汽油中非常规添加剂检测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2041801T3 (es) * | 1987-08-18 | 1993-12-01 | Bp Oil International Limited | Metodo para la determinacion directa de propiedades fisicas de productos hidrocarbonados. |
US5349188A (en) * | 1990-04-09 | 1994-09-20 | Ashland Oil, Inc. | Near infrared analysis of piano constituents and octane number of hydrocarbons |
-
1993
- 1993-10-12 WO PCT/EP1993/002835 patent/WO1994009362A1/en active IP Right Grant
- 1993-10-12 AU AU51514/93A patent/AU677609B2/en not_active Ceased
- 1993-10-12 SG SG1996005908A patent/SG45395A1/en unknown
- 1993-10-12 DK DK93922561.1T patent/DK0671003T3/da active
- 1993-10-12 JP JP6509633A patent/JPH08502360A/ja active Pending
- 1993-10-12 CA CA002146985A patent/CA2146985A1/en not_active Abandoned
- 1993-10-12 BR BR9307250A patent/BR9307250A/pt active Search and Examination
- 1993-10-12 DE DE69307926T patent/DE69307926T2/de not_active Expired - Fee Related
- 1993-10-12 KR KR1019950701442A patent/KR950703735A/ko not_active Application Discontinuation
- 1993-10-12 ES ES93922561T patent/ES2098063T3/es not_active Expired - Lifetime
- 1993-10-12 NZ NZ256694A patent/NZ256694A/en unknown
- 1993-10-12 EP EP93922561A patent/EP0671003B1/en not_active Expired - Lifetime
- 1993-10-12 AT AT93922561T patent/ATE148561T1/de not_active IP Right Cessation
- 1993-10-13 TR TR00917/93A patent/TR27871A/xx unknown
- 1993-10-13 ZA ZA937593A patent/ZA937593B/xx unknown
- 1993-10-13 MY MYPI93002099A patent/MY110065A/en unknown
-
1995
- 1995-04-11 NO NO951416A patent/NO951416L/no unknown
- 1995-04-11 FI FI951720A patent/FI951720A0/fi unknown
-
1997
- 1997-03-26 GR GR970400618T patent/GR3022931T3/el unknown
- 1997-06-26 HK HK91297A patent/HK91297A/xx not_active IP Right Cessation
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004514902A (ja) * | 2000-11-20 | 2004-05-20 | ザ、プロクター、エンド、ギャンブル、カンパニー | ポリマーの予測方法 |
JP2004132921A (ja) * | 2002-10-15 | 2004-04-30 | Yukihiro Toyoda | 性状予測方法 |
JP2007183242A (ja) * | 2005-12-29 | 2007-07-19 | Ind Technol Res Inst | 移動式燃料分析装置および該移動式燃料分析装置を用いた燃料の品質を測定する方法 |
Also Published As
Publication number | Publication date |
---|---|
KR950703735A (ko) | 1995-09-20 |
TR27871A (tr) | 1995-10-11 |
EP0671003B1 (en) | 1997-01-29 |
DE69307926T2 (de) | 1997-05-15 |
DE69307926D1 (de) | 1997-03-13 |
NO951416D0 (no) | 1995-04-11 |
FI951720A (fi) | 1995-04-11 |
NO951416L (no) | 1995-04-11 |
GR3022931T3 (en) | 1997-06-30 |
NZ256694A (en) | 1996-02-27 |
CA2146985A1 (en) | 1994-04-28 |
SG45395A1 (en) | 1998-01-16 |
AU5151493A (en) | 1994-05-09 |
FI951720A0 (fi) | 1995-04-11 |
ES2098063T3 (es) | 1997-04-16 |
BR9307250A (pt) | 1999-08-24 |
AU677609B2 (en) | 1997-05-01 |
MY110065A (en) | 1997-12-31 |
WO1994009362A1 (en) | 1994-04-28 |
DK0671003T3 (da) | 1997-08-18 |
EP0671003A1 (en) | 1995-09-13 |
ZA937593B (en) | 1994-05-03 |
HK91297A (en) | 1997-08-01 |
ATE148561T1 (de) | 1997-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5572030A (en) | Method for determining parameter of hydrocarbon | |
EP0742901B1 (en) | Cracking property determination | |
AU676854B2 (en) | An apparatus for fuel quality monitoring | |
JPH08505944A (ja) | 炭化水素製品の物理特性データの予測方法 | |
JPH08502360A (ja) | 軽油のセタン価の予測方法 | |
Xie et al. | Calibration transfer as a data reconstruction problem | |
EP3861320B1 (en) | Systems and methods for implicit chemical resolution of vacuum gas oils and fit quality determination | |
Amat-Tosello et al. | Contribution of external parameter orthogonalisation for calibration transfer in short waves—Near infrared spectroscopy application to gasoline quality | |
CN102539370B (zh) | 一种多组分气体傅立叶变换红外光谱在线分析滤波方法 | |
CN116380816A (zh) | 一种煤质在线检测系统故障判定方法、程序以及存储介质 | |
CN114428067A (zh) | 一种预测汽油辛烷值的方法 | |
NL1003058C2 (nl) | Werkwijze voor het bepalen van parameters van koolwaterstof. | |
de Freitas et al. | Contributions of multivariate statistics in oil and gas industry | |
Gemperline | Rugged spectroscopic calibration using neural networks | |
BIAKTLUANGA et al. | IDENTIFICATION AND MEASUREMENT OF ADULTERATION IN MIZORAM GASOLINE USING FTIR SPECTROSCOPY AND CHEMOMETRIC METHODS | |
MXPA96001605A (en) | Method to determine hydrocarbon parameters |