JP2011106970A - 誤差推定装置および誤差推定プログラム - Google Patents

Abstract

【課題】シミュレーション結果に含まれるパラメータに起因する誤差を定量的に推定する。
【解決手段】誤差推定装置１０は、線形結合定数演算部１５と誤差演算部１６とを備える。線形結合定数演算部１５は、シミュレーションの入力値の不確かさの割合の関係を示す共分散行列をＷとし、目的体系および実験体系のシミュレーションへの入力値の感度係数ベクトルをそれぞれＳ_ＲおよびＳ_ｉとしたときにＳ_ｉを線形結合したＳ_ｋの線形結合定数α_ｉを、Ｓ_ｋ ^ＴＷＳ_ｋ＝Ｓ_Ｒ ^ＴＷＳ_Ｒを満足しかつＳ_ＲとＳ_ｋとのなす角が最小になるように求める。誤差演算部１６は、実験で測定された物理量Ｒの測定値に対するシミュレーションで得られた物理量Ｒの計算値のシミュレーションのモデルへの入力値に関する相対誤差をＥ_Ｐｉとしたときに、α_ｉを重みとしてＥ_Ｐｉの値を合成してシミュレーションの結果得られる物理量に含まれる相対誤差を推定する。
【選択図】図１

Description

本発明は、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に含まれる誤差を、その対象を模擬した実験の結果を用いて推定する誤差推定装置および誤差推定プログラムに関する。

工業分野で製作される製品は、実際に使用する前に目的とする性能が発揮できるかを確認する必要がある。最もよく用いられる方法は、最終的な製品と殆ど同じものを試作して、その試作品の性能を確かめる方法で、実証試験と呼ばれている。たとえば家電製品でも自動車でも鉄道列車でも試作品を製作して、その性能を確認し品質を保証して最終的な製品として出荷している。

一方、特定の工業分野では、製作する製品の大きさが巨大で、最終的な製品の性能を試作によって確認することが極めて困難である場合がある。また、製作する製品が非常に高価であるか製作数が唯一つあるいは少数であるか、経済的あるいは他のいくつかの理由によって試作品を製作することが合理性をもたない場合もある。

たとえば、実際の原子力発電所や核物質に関係する施設では、これらの発電所や施設を試作することは経済的にも物理的にもほとんど不可能である。理由の一つとして、使用する核物質を自由に用意することができないことがある。使用済み核燃料の貯蔵施設を建設する場合、予め貯蔵する使用済み核燃料と同じものを用意しておいて施設の性能を測定することは現実的には不可能である。原子力発電所を建設する際も、予め別の燃料集合体を用意しておいてその核燃料によって炉心の核的性能を確認することは、経済的な面からも成立性がほとんど無い。

そこで、原子力分野では実物を用いた性能の確認ではなく、計算による性能の確認が産業分野の誕生時点から行われてきた。また他の分野、特に高額な製品である航空機やロケット、大型船舶などの設計でも、製品のコストを下げ、製作時間を短縮するために設計段階で計算によって最終的な性能を把握することが、ますます一般的になってきている。

予測計算には、一般的にコンピュータを用いる。製品が従う物理理論に基づいて作成された計算機プログラムによってその「製品」の性能を把握する技術は、計算物理、計算実験あるいは計算機シミュレーション（コンピュータシミュレーション）などと呼ばれ、大きく発達した技術分野になっている。

原子力発電所を建設する場合、設計段階において、実際の燃料集合体を用意して原子炉の炉心の核的性能を確認することは不可能であるので、コンピュータシミュレーションで性能を確認する。また原子炉の炉心に装荷する製作する燃料集合体の体数や燃料集合体に含まれる核物質の濃度や分布についてもコンピュータシュミレーションに基づいて決定される。原子力分野ではこのようなコンピュータシュミレーション技術は不可欠である。

コンピュータシュミレーションは、数値計算の積み重ねであって、必ず計算誤差が付きまとう。計算の確からしさを計算精度と呼ぶが、計算精度を把握することとは設計の信頼性・安全性の向上、経済性の改善、設計の合理化のために非常に重要である。計算精度を正確に把握することとは、得られた計算値に伴っている計算誤差を正確に把握することと同じ意味である。

次に、計算誤差の種類について説明する。

計算で生じる誤差の要因は幾つかの種類に分けられる。計算手法の誤差、数値計算の誤差および計算に用いる数値（入力値）に起因する誤差である。

計算に用いる物理モデルは、あくまで自然界の物理現象を数学的にモデル化したものであり、自然界がこの物理モデルに従って変化しているわけではない。加えて、物理モデルをコンピュータで計算できるように近似・変形した場合にも誤差が生じる。これらが計算手法の誤差である。

原子力の核計算において、中性子の運動の基礎方程式は、ボルツマンの輸送方程式によって記述される。ボルツマンの輸送方程式をコンピュータで解く場合、何らかの近似を用いて解くことになる。原子炉の炉心計算では、ボルツマンの輸送方程式を近似した拡散方程式が広く利用されている。この拡散方程式は、ボルツマンの輸送方程式を中性子の運動の角度方向の情報を無視して近似したものである。よって中性子の運動方向が均一でない体系では、拡散方程式による計算で誤差が大きくなるおそれがある。加えて、拡散方程式で使用される拡散係数は、密度の薄い媒質に対して適切な数値を設定することができない場合があり、そのため誤差を生ずる。計算対象の体系が小さい場合、すなわち中性子の洩れの割合が大きい小型の炉心を拡散理論で計算すれば、洩れを過大評価して臨界固有値を過小評価する傾向があり、計算誤差が大きくなる傾向があるといわれている。

このように、計算を行なう際、計算に用いる理論やモデルが適切であるか、計算手法で生じている誤差の検討が必要である。

数値計算の誤差は、コンピュータで計算する過程の四則演算で生じる誤差である。数値計算の誤差として代表的なものには、たとえば丸め誤差とか、切り上げ・切り捨て誤差がある。ただし、数値計算の誤差は、現在の工夫された数値計算技術とコンピュータのハードウェアの発達によって十分無視できるレベルになっており、一般に他の誤差と比較して十分小さく、誤差の主たる支配因子ではない。

誤差の要因として一番大きいと考えられるものは、計算に使用する数値に起因する誤差である。計算に用いられる数値は、寸法、体系を構成する物質、その物質の数密度などの計算対象に固有な数値と、計算対象それぞれには依存せず、計算に共通に使用される数値がある。たとえば、水の温度と圧力が決まれば水の密度は決まるが、その密度の値は計算体系には依存せず共通に使用される値である。計算に用いられるこれら全ての数値を、以下パラメータと呼ぶことにする。またパラメータによって計算結果に生ずる誤差を、以下パラメータ誤差と呼ぶことにする。

一般的に用いられる物理定数も真の値ではなく、誤差が含まれている。ただし、物理定数の精度は他の数値に比べて遥かに高く、有効数字も６桁を越えるものが多くある。このため、一般的に、物理定数の誤差が原子力施設の設計や建設、燃料集合体のなどの設計、製作で問題となることはないと判断できる。

一方、これまで核計算の中で最も重要と判断されてきたパラメータは「核データライブラリ」に関する数値であり、中性子反応断面積、崩壊定数、収率、遅発中性子割合などがそれにあたる。「核データライブラリ」は、核計算に直接的に関わる数値群で、これら数値の変化が核特性を示す計算値に与える影響が大きい。また、もともと「核データライブラリ」は、全ての値が測定によって正確に確認・決定されたものではなく、理論計算によって定められた数値も含まれる。このため、「核データライブラリ」に含まれる誤差は、他のパラメータ誤差よりも大きいと判断される。したがって、核計算では、「核データライブラリ」のパラメータ誤差が重要視される。

パラメータ誤差のなかで次に重要なものは、原子数密度の誤差である。原子力分野において、運転後の燃料集合体などについては、燃焼計算によって原子数密度を得る。このため、その後の計算で原子数密度を用いるときには、既にこの原子数密度に計算誤差が含まれていると判断される。

原子力産業では、扱う製品の構成元素、元素の原子数密度などの組成を高い精度で確認・測定することが困難あるいは非現実的な場合がある。燃料集合体内の核種組成や核種の原子数密度は、次にその燃料集合体を使用する場合の反応度の決定、保管・輸送の際の臨界安全性や放射線量を定量化するために極めて重要な数値である。しかし、核種の組成や原子数密度の測定のために燃料集合体を炉心から取り出すこと、さらに、破壊して分析測定することは現実的ではなく、また、高い放射線レベルの燃料集合体１体ずつを外部から正確に測定することも技術的にも経済的にも成立性が乏しい。

そこで燃料集合体の核種組成や原子数密度は、過去燃料集合体が置かれていた環境をできるかぎり正確に把握して、計算機プログラムによって計算する。これを燃焼計算と呼ぶ。このようにして求めた原子数密度には、計算誤差が含まれる。この原子数密度もこの数値を入力とした核計算の計算結果に与える影響が大きいものであり、計算で生じる誤差の原因を考察する際に重要である。

実際の原子力発電所に関する計算を行なう場合には、原子炉の中で移動する流体の流量や温度、材料温度は、計算に必要な入力点数について精度の高い数値が得られないことがあり、経験的な数値を仮定し入力して計算を行なう。炉心の中に水と蒸気の二層流が存在する沸騰水型原子炉（ＢＷＲ）の炉心計算では、流量の誤差によって生ずる計算誤差は有意な値である考えられている。

原子力産業の初めから、核物質の臨界性という原子力に特有の現象を確認するために、実験装置が作られ、利用されてきた。そのひとつが、臨界実験装置である。臨界実験装置は、一般的に、大気圧下、常温（室温）で運転・稼動できるように設計された装置である。

臨界実験装置は小型の原子炉であって、実際の原子炉で使用するウランやプルトニウムなどの核物質を使用して臨界状態を達成する。しかし、装置が非常に小型であるために熱をほとんど出さない。原子炉のミニチュア版といった装置であり、臨界状態を実現できることから、臨界実験装置と呼ばれている。臨界実験装置を用いた実験は、臨界実験と呼ばれる。

臨界実験は、計算の誤差を減らすために役立てられてきた。物理的に非常に単純化され、簡略化された条件で、体系を組み上げて、同時に精度の高い測定データを取得することが臨界実験の目的である。物理的に非常に単純で簡略化した体系では、形状や組成に関わる計算入力パラメータの誤差を減らすことができ、誤差の少ない測定値は計算値との比較を容易にする。過去から現在に至るまで「計算手法に起因する誤差」と「核データライブラリに起因する誤差」を明らかにすることを目的に、臨界実験によって得られた多くの測定値と計算値が比較されて、その結果、誤差要因が特定され誤差が定量化されてきた。

加えて、臨界実験の測定値を計算値がよく再現できていれば、計算に使った手法、計算に使った主として核データライブラリなどのパラメータの品質が高いと判断される。品質が高いと判断されれば、同じ手法とパラメータの組み合わせで目的とする体系の設計を行っても良いという品質保証・判断基準となってきた。

臨界実験の測定値で一番重要視されるものは、臨界になった条件である。そこで、臨界になった諸条件を、データとして正確に取得する。原子炉物理の言葉では、臨界になった条件を「臨界質量」という言葉で表すことがあり、これは核物質種類、質量のみならず、体系の幾何形状、温度や核物質の組成や質量一式の正確な数値を指す。

次に重要とされるものは、臨界実験装置が臨界になったときの核分裂反応の空間分布である。燃料棒を組み合わせて構成された臨界実験装置では、燃料棒から放出される放射線を測定して、放射線量の比によって核分裂反応の分布を測定することが多い。これは核分裂反応と放出される放射線の量は比例関係にあると考えられるからである。なお、臨界実験装置に、たとえば小型核分裂電離箱などの特別の放射線測定器を挿入して、目的とする位置での放射線の量を測定することも多い。この場合も、測定値の比によって中性子束の分布などが求められる。

目的とする測定値が得られれば、次に臨界実験の体系や実験条件を入力として計算を行なう。計算で求めた値を臨界実験で得られた臨界量と核分裂反応の分布、中性子束の分布などの測定値と比較することで、計算全体の品質が把握できる。すなわち計算値と測定値が測定誤差の範囲や許容できる程度で一致していれば、臨界実験を計算した計算機プログラム（計算手法）と核データライブラリを中心として計算に用いたパラメータが十分な品質を有しているという根拠になる。その後、同じ計算機プログラム（計算手法）と核データライブラリを用いて最終的に目的とする原子炉や原子炉施設の設計の計算に適用できるという判断がなされる。

一方、計算値と測定値に有意な差が認められる場合は、計算手法や核データライブラリの問題点や改良すべき点を特定し、計算値と測定値の一致が改善されるように改良がなされる。このように、臨界実験は、計算機（プログラム）と核データライブラリの品質の確認や保証、あるいは計算手法や核データライブラリの改良に寄与してきた。

臨界実験で得られた数値（測定値）にも誤差が含まれる。その誤差は、測定誤差と呼ばれる。測定誤差には、それぞれの測定の際に偶然発生する統計誤差（ランダム誤差、統計誤差）と、測定に用いた計測器や手法に伴う誤差（系統誤差）が含まれる。統計誤差は、確率分布に従い、測定の回数を増やせば誤差の割合は減少する性質がある。系統誤差は、測定方法そのものに付随した誤差であるので、測定ごとに常に発生し、系統誤差の割合は測定の回数によって変化することはない。なお、測定誤差の割合は、計算誤差の割合よりも小さいと考えられる。

これまで臨界実験で得られたデータを利用する方法として、一般的に次の二つの応用がなされてきた。補正因子（バイアス）法と、断面積アジャストメントである。

補正因子（バイアス）法では、計算によって生じる誤差を減らすことを目的としている。特定の計算手法（計算機プログラム）と特定の核データライブラリ（計算で用いるパラメータ）を使用して、目的とする体系の設計計算を行なうものとする。特定の物理量に注目して、臨界実験で得られた計算値の相対誤差の割合が、目的とする体系の計算値の相対誤差の割合と同じと仮定する。

次に、臨界実験の計算値と測定値とから、補正因子（バイアス）を計算する。目的とする体系の計算値に、この補正因子を乗じれば目的とする体系の計算値の相対誤差を取り除くことが可能になり、目的とする体系の計算精度を向上させることができる。

臨界実験の真値をＴ_Ｅ、臨界実験での計算値をＣ_Ｅ＝Ｔ_Ｅ（１＋Δｃ_ｅ）、臨界実験での測定値をＲ_Ｅ＝Ｔ_Ｅ（１＋Δｅ）とする。ここで、Δｃ_ｅは臨界実験での計算誤差、Δｅは測定誤差である。

ｆ＝Ｒ_Ｅ／Ｃ_Ｅ＝Ｔ_Ｅ（１＋Δｅ）／Ｔ_Ｅ（１＋Δｃ_ｅ）＝（１＋Δｅ）／（１＋Δｃ_ｅ）であるから、このｆを目的とする計算値Ｃ_ｄ＝Ｔ_ｄ（１＋Δｃ_ｄ）に乗ずれば、一般に実験の測定誤差Δｅは計算誤差Δｃ_ｅよりも小さいので、Δｃ_ｄとΔｃ_ｅとが同じ程度であるならば、
Ｒ_ｄ＝Ｃ_ｄ×ｆ
＝Ｔ_ｄ（１＋Δｃ_ｄ）×（１＋Δｅ）／（１＋Δｃ_ｅ）
≒Ｔ_ｄ（１＋Δｅ＋Δｃ_ｄ−Δｃ_ｅ）
≒Ｔ_ｄ（１＋Δｅ）
となり、目的とする体系の計算値Ｒ_ｄから計算による誤差を減らした値を推定することができる。これが補正因子（バイアス）法である。なお、Ｔ_ｄは真値、Δｃ_ｄは目的とする計算値に含まれる計算誤差である。

補正因子（バイアス）法では、目的とする体系において計算で生じる相対誤差を減らすことを目的としており、その誤差の原因を取り除くことは考察しない。臨界実験の計算値と測定値を利用して目的とする体系の計算精度を向上させる別の手法として、断面積アジャストメントがある。計算誤差の要因を核データライブラリ（断面積ライブラリ）と考えて、臨界実験での計算値と測定値とが一致するように、核データライブラリの値を調整・修正（アジャストメント）し、その調整を行った核データを用いて、目的とする体系の計算を行ない、計算精度を向上させるという考え方である。

すなわち、複数の臨界実験で得られた測定値と計算値の相対差の合計（通常は二乗和）を最小にするように、核データライブラリ（断面積）の数値を調整する手法である。核特性に対するそれぞれの断面積の感度係数を用いて、測定値と計算値の相対差の二乗和を最小にするように最終的に断面積を調整することを目的とする手法である。行われる数学的演算は、行列形式の最小二乗法であり、必然的に補正を行う核データライブラリは多群形式のライブラリである。よって、計算も中性子エネルギーを離散化した多群計算を出発点としており、連続エネルギーモンテカルロ計算などに、直接的には適用できない。

この手法の優れた点は、数学的な式の扱いが最小二乗法に基づいているので、感度係数を用いて断面積の値を調整する場合の数学的取り扱いが非常に明快なことである。計算手法や計算機プログラムで生じる誤差が無視できると考えられる場合で、設計計算に関わる不確かさ（誤差）が計算に使用しているパラメータの不確かさにのみ起因しており、パラメータの中でも断面積ライブラリが誤差の最も大きな要因ならば、この手法の適用性は高いといえる。

一方、この手法の問題点は、計算誤差（不確かさ）の原因を断面積ライブラリに限定していることで、たとえば他の未知の原因によって臨界実験での計算値と測定値の差が生じている場合は、的外れな操作を行っていることである。また、仮に断面積ライブラリが計算値と測定値の差の支配因子であっても、臨界実験（あるいはこれまで知りえた比較例）で計算値と測定値が一致することはあくまで必要条件であって、そのように調整した断面積ライブラリを使用すれば目的とする設計体系の誤差（不確かさ）を必ず小さくできるという根拠や理由（十分条件）は、示されていない。

特許文献１には、主に流体の計算に際して、幾何学的な寸法誤差・公差、境界条件など入力値の持つ誤差、数値解析誤差について、計算値に含まれる誤差を、感度係数を用いて評価する手法が開示されている。また、非特許文献１には、目的とする体系に対して最終的な設計値を得る際に、計算値にバイアスを乗じ、計算値を補正する手法が開示されている。この手法では、実験で得られた測定値と計算値との比を組み合わせることで、一般的なバイアスを得る。

特許文献１に記載された手法は、あくまで計算値だけを利用しており、実験で得られる測定値を有効的に利用する手段がない。また、非特許文献１に記載された手法は、数学的手法が記載されているものの、指数関数を利用するため演算が極めて複雑であり各実験に必要な係数を一意的に簡便に計算できない。さらに、基本的問題として計算のアルゴリズムに、なぜ指数関数を用いるべきかの物理学的な理由付けがないので、理論の基盤が十分でない。

特開２００８−２１７１３９号公報

Teruhiko KUGO、他２名、"Theoretical Study on New Bias Factor Methods to Effectively Use Critical Experiments for Improvement of Prediction Accuracy of Neutronic Characteristics"、Journal of NUCLEAR SCIENCE and TECHNOLOGY、Vol.44, No.12、Page 1509-1517、２００７年

原子力分野や、試作品を製作して最終的な性能を確認することが現実的ではない分野では、特定の設計手法に基づき、特定のパラメータを入力として目的とする体系の着目する物理量を計算して設計値としている。この計算で用いる方程式を、ここではモデル方程式と呼ぶ。設計計算の段階で計算値である設計値に含まれる数値の誤差を定量評価することは、設計精度を確認するために極めて重要な作業である。

設計精度を確認するために、原子力分野では臨界実験が行われてきた。実験において着目する物理量を測定し、設計する際に使用するのと同じ計算手法（計算機プログラム）と物理定数など（パラメータ）を用いて計算を行ない、計算値と測定値の差によって設計計算での誤差を推測する。あるいは計算値と測定値の差がある許容限界の範囲ならば、この方法で設計計算を行なっても問題はないという判断がなされてきた。

臨界実験は、従来、理論や計算手法の誤差と、パラメータによる誤差とを定量化することを目的として実施されてきた。ただし、臨界実験で得られる情報について、これまで物理学的・数学的に詳しい検討がなされることは少なく、目的とする体系と比較して重要な特徴が再現されている実験であれば有用であるという判断に基づいて実験がなされてきた。

臨界実験の測定値と計算値とは、補正因子（バイアス）法や断面積アジャストメントに利用されてきた。よって、臨界実験で得られた計算誤差の情報を組み合わせて、目的の体系の計算誤差にする手法は提供されていなかった。ここで情報とは、目的とする体系を、特定の設計手法、特定のパラメータで計算したときに得られる数値（設計予想値）に対して、誤差がどの程度含まれるかを確認するための情報である。

実験結果に対して定量的、数学的に詳しい検討がなされなかったため、実験を行う際に、着目する物理量の性質を明らかにするには実験においてどの項目を重要視すべきかの判断ができないので、最善の実験を行ったかどうか判断できないという問題がある。

また、目的とする体系の状態を常に広い範囲について実験で模擬できるとは限らない。たとえば目的とする体系で使用する物質について実験施設で十分な量を用意できない場合、実験ではその物質を使用している空間領域は実験体系の一部になり、このような実験は部分模擬実験と呼ばれる。部分模擬実験で得られた情報で目的とする体系の着目する物理量の誤差を評価する手法は知られていない。

さらに、複数の実験を行った際、どの実験が優れていたかを定量的に比較する方法、複数の実験の情報を組み合わせて設計計算の誤差を推定する手法は知られていない。その他、過去に実施した実験の情報を有効利用したり、実験に必要な資源（人的、時間的、経済的資源）の最適化を行う有効な手法がないのが現状である。

そこで、本発明は、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に含まれるパラメータに起因する誤差を、その対象を模擬した実験の結果を用いて定量的に推定することを目的とする。

上述の目的を達成するため、本発明は、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に含まれる誤差を、その対象を模擬した実験の結果を用いて推定する誤差推定装置において、ｉを１から実験の総数ｎまでの自然数としたときにｉ番目の実験についての前記モデルを用いたシミュレーションで得られたある物理量Ｒの計算値の当該実験で測定された当該物理量Ｒの測定値に対する前記モデルへの入力値に関する相対誤差Ｅ_Ｐｉを記憶する相対誤差記憶部と、シミュレーションに用いる前記モデルへの入力値の不確かさの割合の関係を示す共分散行列Ｗを記憶する共分散行列記憶部と、前記対象について前記モデルを用いたシミュレーションの結果のそのモデルへの入力値の単位変化に対する変化量を示す目的体系感度係数ベクトルＳ_Ｒを記憶する目的体系感度係数ベクトル記憶部と、ｉ番目の実験の体系について前記モデルを用いたシミュレーションした結果のそのモデルへの入力値の単位変化に対する変化量を表す実験体系感度係数ベクトルＳ_ｉを記憶する実験体系感度係数ベクトル記憶部と、前記実験体系感度係数ベクトルＳ_ｉを線形結合した線形結合ベクトルＳ_ｋの線形結合定数α_ｉを、Ｓ_ｋ ^ＴＷＳ_ｋ＝Ｓ_Ｒ ^ＴＷＳ_Ｒを満足しかつ前記目的体系感度係数ベクトルＳ_Ｒと前記線形結合ベクトルＳ_ｋとのなす角が最小になるように求める線形結合定数演算部と、前記線形結合定数α_ｉを重みとして前記相対誤差Ｅ_Ｐｉの値を合成して前記モデルを用いたシミュレーションの結果得られる前記対象の物理量Ｚに含まれる相対誤差ΔＺ／Ｚを推定する誤差演算部と、を有することを特徴とする。

また、本発明は、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に含まれる誤差を、その対象を模擬した実験の結果を用いて推定する誤差推定プログラムにおいて、コンピュータに、ｉを１から実験の総数ｎまでの自然数としたときにｉ番目の実験についての前記モデルを用いたシミュレーションで得られたある物理量Ｒの計算値の当該実験で測定された当該物理量Ｒの測定値に対する前記モデルへの入力値に関する相対誤差Ｅ_Ｐｉを記憶する相対誤差記憶機能と、シミュレーションに用いる前記モデルへの入力値の不確かさの割合の関係を示す共分散行列Ｗを記憶する共分散行列記憶機能と、前記対象について前記モデルを用いたシミュレーションの結果のそのモデルへの入力値の単位変化に対する変化量を示す目的体系感度係数ベクトルＳ_Ｒを記憶する目的体系感度係数ベクトル記憶機能と、ｉ番目の実験の体系について前記モデルを用いたシミュレーションした結果のそのモデルへの入力値の単位変化に対する変化量を表す実験体系感度係数ベクトルＳ_ｉを記憶する実験体系感度係数ベクトル記憶機能と、前記実験体系感度係数ベクトルＳ_ｉを線形結合した線形結合ベクトルＳ_ｋの線形結合定数α_ｉを、Ｓ_ｋ ^ＴＷＳ_ｋ＝Ｓ_Ｒ ^ＴＷＳ_Ｒを満足しかつ前記目的体系感度係数ベクトルＳ_Ｒと前記線形結合ベクトルＳ_ｋとのなす角が最小になるように求める線形結合定数演算機能と、前記線形結合定数α_ｉを重みとして前記相対誤差Ｅ_Ｐｉの値を合成して前記モデルを用いたシミュレーションの結果得られる前記対象の物理量Ｚに含まれる相対誤差ΔＺ／Ｚを推定する誤差演算機能と、を実現させることを特徴とする。

本発明によれば、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に含まれるパラメータに起因する誤差を、その対象を模擬した実験の結果を用いて定量的に推定できる。

本発明に係る計算誤差推定装置の一実施の形態におけるブロック図である。 本発明に係る計算誤差推定装置の一実施の形態におけるコンピュータのブロック図である。 本発明に係る計算誤差推定装置の一実施の形態における誤差推定方法のフローチャートである。

本発明に係る計算誤差推定装置の実施の形態を、図面を参照して説明する。なお、以下の説明は単なる例示であり、本発明はこれらに限定されない。

図１は、本発明に係る計算誤差推定装置の一実施の形態におけるブロック図である。図２は、本実施の形態におけるコンピュータのブロック図である。

計算誤差推定装置１０は、相対誤差記憶部１１と、共分散行列記憶部１２と、目的体系感度係数ベクトル記憶部１３と、実験体系感度係数ベクトル記憶部１４と、線形結合定数演算部１５と、誤差演算部１６と、を有している。この誤差推定装置１０は、コンピュータ２１上に構築することができる。

コンピュータ２１では、ＣＰＵ（Central Processing Unit）２２、ＲＡＭ（Random Access Memory）２３、ＲＯＭ（Read Only Memory）２４、入力制御部２６、表示制御部２９、ハードディスク制御部３１などがバス３２を介して接続されている。入力制御部２６には、キーボードやマウスなどの入力装置２５が接続されている。表示制御部２９には、液晶ディスプレイなどの表示装置２８が接続されている。ハードディスク制御部３１には、ハードディスク３０が接続されている。

ＣＰＵ２は、計算誤差を推定する誤差推定処理、およびその他の様々な演算処理や制御処理を行う。ＲＡＭ２３は、ＣＰＵ２２が処理を行う際に一時的にデータを記憶する。ＲＯＭ２４あるいはハードディスク３０は、測定値を介して得られる相対誤差、共分散（誤差）行列、計算して得られた感度係数などをＣＰＵ２２が処理するために必要なデータやプログラムを長期的に記憶する。演算開始の指示などコンピュータ２１への入力は、入力装置２５を介して行われる。推定された誤差などの必要な情報は、表示装置２８に表示される。

誤差推定装置１０は、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に付随する計算誤差を推定する。ここで、計算誤差推定装置１０が計算誤差を推定する体系を「目的の体系」あるいは「設計体系」と呼ぶこととする。この目的の体系とは、たとえば原子炉などの製品である。

計算誤差を推定するときに、誤差推定装置１０は、ｎ個の実験によって得られた結果を用いる。ここで、実験とは、目的の体系の少なくとも一部を模擬しておこなうものであり、実験に用いた装置などの体系を「実験体系」と呼ぶこととする。本実施の形態では、模擬実験として臨界実験を行う場合について説明する。

相対誤差記憶部１１は、ｉ番目の実験についてのシミュレーションで得られたある物理量Ｒの計算値の当該実験で測定された当該物理量Ｒの測定値に対するシミュレーションのモデルへの入力値に関する相対誤差Ｅ_Ｐｉを記憶する。ここで、ｉは、１から実験の総数ｎまでの自然数である。

共分散行列記憶部１２は、シミュレーションに用いるモデルへの入力値の不確かさの割合の関係を示す共分散行列Ｗを記憶する。目的体系感度係数ベクトル記憶部１３は、目的の体系についてのシミュレーションの結果の、そのシミュレーションに用いたモデルへの入力値の単位変化に対する変化量を示す目的体系感度係数ベクトルＳ_Ｒを記憶する。実験体系感度係数ベクトル記憶部１４は、ｉ番目の実験の体系についてのシミュレーションの結果の、そのシミュレーションに用いたモデルへの入力値の単位変化に対する変化量を表す実験体系感度係数ベクトルＳ_ｉを記憶する。

線形結合定数演算部１５は、実験体系感度係数ベクトルＳ_ｉを線形結合した線形結合ベクトルＳ_ｋの線形結合定数α_ｉを、Ｓ_ｋ ^ＴＷＳ_ｋがＳ_Ｒ ^ＴＷＳ_Ｒと等しくかつ目的体系感度係数ベクトルＳ_Ｒと線形結合ベクトルＳ_ｋとのなす角が最小になるように求める。誤差演算部１６は、線形結合定数α_ｉを重みとして相対誤差Ｅ_Ｐｉの値を合成してシミュレーションの結果得られる目的の体系の物理量Ｚに含まれる相対誤差ΔＺ／Ｚを推定する。

図３は、本実施の形態における誤差推定方法のフローチャートである。

まず、臨界実験などの模擬実験を行う（Ｓ１０１）。着目する物理量Ｒに対してこのとき得られた測定値をＲ_Ｅとする。

次に、あるモデルを用いてシミュレーション計算を行う（Ｓ１０２）。このとき得られた計算値をＲ_Ｃとする。

計算値と測定値の相対差は、
（Ｒ_Ｃ−Ｒ_Ｅ）／Ｒ_Ｅ＝（Ｒ_Ｃ／Ｒ_Ｅ）−１
と表わされる。この値は、
「パラメータによる相対誤差」＋「計算手法による相対誤差」−「測定誤差」
にほぼ等しい。なお、ここでの測定誤差は相対誤差である。

この値から「計算手法による相対誤差」と「測定誤差」を取り除く。「計算手法による相対誤差」および「測定誤差」は、いずれも通常、正負が明らかではないから、二乗して、
｛（Ｒ_Ｃ／Ｒ_Ｅ）−１｝^２−「計算手法による相対誤差」^２−「相対測定誤差」^２
を実験体系でのパラメータによって生じている相対誤差Ｅ_Ｐｉの二乗とみなす。なお、この演算は正確な評価ではないが、便宜上の評価である。計算手法による相対誤差が明確でないときには、０とおいてもよい。

このようにして得られたｉ番目の実験体系でのパラメータによって生じている相対誤差Ｅ_Ｐｉは、誤差推定装置１０に取り込まれて、相対誤差記憶部１１に記憶される（Ｓ１０３）。あるいは、誤差推定装置１０は、Ｒ_ＥとＲ_Ｃとを入力されて実験体系でのパラメータによって生じている相対誤差Ｅ_Ｐｉを求めてもよい。ステップＳ１０１ないしステップＳ１０３は、実験の総数ｎ回繰り返される。

次に、目的とする体系と実験体系とで使用する同じ構造の入力値の不確かさを表す共分散（誤差）行列を設定する。目的とする体系と実験体系で共通に使用するこの共分散（誤差）行列をＷとする。Ｗの対角成分ｗ_ｉｉにはそれぞれのパラメータの相対誤差を二乗した値Ｅ_Ｐｉ ^２が収められており、非対角成分ｗ_ｉｊにはパラメータｉとパラメータｊの間の相対誤差の積Ｅ_ＰｉＥ_Ｐｊが収められている。

ここで、この共分散（誤差）行列Ｗの各成分の大きさは、互いの大きさの関係が正しければ良く、絶対値は問題ではない。このＷは、記憶装置に記憶される（Ｓ１０４）。なお、ここでは、目的とする体系の入力に関しての不確かさを表した共分散（誤差）行列とそれぞれの実験についての体系の入力に関しての不確かさを表した共分散（誤差）行列が同じ行列という仮定をしている。

一般に、モデルへの入力値の不確かさの割合を表す共分散（誤差）行列は、必ずしも明確に定義されない。この場合は、単位行列を共分散（誤差）行列としてもよい。

着目する物理量をＲ、このＲを計算するときに用いるパラメータをｘ_ｉ（ｉ＝１，２，３，…，ｎ）とする。ｘ_ｉに関するＲの感度係数Ｓは

で定義される。すなわち、感度係数Ｓパラメータｘ_ｉの変化割合の単位量あたりに生じる物理量Ｒの変化割合を示している。着目する物理量について目的とする体系、実験体系に関してそれぞれ各パラメータに対しての感度係数を求める。感度係数はパラメータの数に対するベクトル量となる。

目的とする体系の感度係数ベクトルＳ_Ｒは目的体系感度係数ベクトル記憶部１３に、ｉ番目の実験体系の感度係数ベクトルＳ_ｉとは実験体系感度係数ベクトル記憶部１４にそれぞれ記憶される（Ｓ１０５）。これらの感度係数ベクトルは、誤差推定装置１０内で計算してもよいし、外部から入力されてもよい。

次に、目的とする体系の感度係数ベクトルＳ_Ｒに対してｉ番目の実験体系の感度係数ベクトルＳ_ｉを用いて、（１）式を満足する感度係数ベクトルの線形結合（線形結合ベクトル）Ｓ_ｋを作る（Ｓ１０６）。

ここで、α_ｉは任意定数で、線形結合定数と呼ぶこととする。

線形結合定数演算部１５は、線形結合定数α_ｉを、（２）式を満足し、かつ目的体系感度係数ベクトルＳ_Ｒと線形結合ベクトルＳ_ｋとのなす角θが最小になるように求める（Ｓ１０７）。

Ｓ_ｋ ^ＴＷＳ_ｋ＝Ｓ_Ｒ ^ＴＷＳ_Ｒ …（２）
ここで、Ｔはベクトルや行列の転置を示す記号である。この（２）式は、相対誤差の総量が同じであるという条件を表している。

目的体系感度係数ベクトルＳ_Ｒと線形結合ベクトルＳ_ｋとのなす角θが最小になるように求める方法はいくつか想定される。たとえば、（Ｓ_Ｒ−Ｓ_ｋ）^ＴＷ（Ｓ_Ｒ−Ｓ_ｋ）の絶対値を最小とすることにより、目的体系感度係数ベクトルＳ_Ｒと線形結合ベクトルＳ_ｋとのなす角θを最小にすることができる。（Ｓ_Ｒ−Ｓ_ｋ）^ＴＷ（Ｓ_Ｒ−Ｓ_ｋ）の値が零（０）となった場合、すなわち、
（Ｓ_Ｒ−Ｓ_ｋ）^ＴＷ（Ｓ_Ｒ−Ｓ_ｋ）＝０ …（３）
を満足する場合、自動的に、
Ｓ_ｋ ^ＴＷＳ_Ｒ／｛（Ｓ_ｋ ^ＴＷＳ_ｋ）^１／２（Ｓ_Ｒ ^ＴＷＳ_Ｒ）^１／２｝＝１…（４）
となる。この（４）式は、模擬性評価因子が１となっていることを示している。（４）式の左辺の分子は、共分散行列Ｗを考慮した線形結合ベクトルＳ_ｋと目的体系感度係数ベクトルＳ_Ｒとの内積を示している。（４）式の左辺の分母は、共分散行列Ｗを考慮した線形結合ベクトルＳ_ｋと目的体系感度係数ベクトルＳ_Ｒのベクトルの大きさになっている。

つまり、模擬性評価因子とは、目的体系感度係数ベクトルＳ_Ｒと線形結合ベクトルＳ_ｋとのなす角をθとしたときの、ｃｏｓθを表している。この値が１のとき、すなわちｃｏｓθ＝１のときに、θ＝０となり、線形結合ベクトルＳ_ｋと目的体系感度係数ベクトルＳ_Ｒとが重なった状態を表現している。

（Ｓ_Ｒ−Ｓ_ｋ）^ＴＷ（Ｓ_Ｒ−Ｓ_ｋ）の絶対値を最小とする線形結合定数α_ｉは、たとえばLagrangeの未定定数法を用いて求めることができる。（２）式の値をｄとおき、この式を束縛条件とし、Lagrangeの未定定数λを用いて、以下の式Ｌを作る。

Ｌ＝｛（Ｓ_Ｒ−Ｓ_ｋ）^ＴＷ（Ｓ_Ｒ−Ｓ_ｋ）｝＋λ｛Ｓ_ｋ ^ＴＷＳ_ｋ−ｄ｝
Ｌ＝｛２ｄ−２（Ｓ_ｋ ^ＴＷＳ_Ｒ）｝＋λ｛Ｓ_ｋ ^ＴＷＳ_ｋ−ｄ｝ …（５）
この（５）式に関して極値条件から、

として、α_ｉの条件式（連立方程式）を作成する。この条件式から、λとα_ｉ（ｉ＝１，２，…，ｎ）を決定する。

今、ｒ_ｉｊ＝Ｓ_ｉ ^ＴＷＳ_ｊ、ｒ_Ｒｉ＝Ｓ_ｉ ^ＴＷＳ_ｊとすれば、Lagrange未定定数λを用いて、λとα_ｉ（ｉ＝１，２，…，ｎ）とを解く方程式は以下のように書ける。

Ｒλα＝ｒ …（６）
ここで、ＲはＲ_ｉｊで構成される行列、αはα_ｉで構成される列ベクトル、ｒはｒ_Ｒｉで構成される列ベクトルである。

例として、Ｒが４行４列の行列である場合について説明する。この場合、上述の方程式は、

となる。この式を解いて、λα_ｉ（ｉ＝１，２，…，ｎ）を求める。すなわち、次式により、λα_ｉ（ｉ＝１，２，…，ｎ）を算出する。

加えて、α_ｉは、次式が成り立つように規格化される。

このようにして、λとα_ｉ（ｉ＝１，２，…，ｎ）が決定される。

なお、λとα_ｉ（ｉ＝１，２，…，ｎ）は、正負で与えられる。原理的には、２ｄ−２（Ｓ_ｋ ^ＴＷＳ_Ｒ）が小さくなるほうの符号を選ぶことになる。しかし、新たに定義される模擬性評価因子がλに等しいことがわかっており、模擬性評価因子は１に近い値であるべきなので、λを負の値とすることは目的と合致しない。そこで、λとしては常に正の値を選択すればよいことになる。

次に、誤差演算部１６は、目的の体系の着目する物理量Ｚの誤差ΔＺに対して、計算に用いたパラメータに起因する相対誤差の二乗（ΔＺ／Ｚ）^２を、線形結合定数α_ｉを用いて求める（Ｓ１０８）。この相対誤差の二乗（ΔＺ／Ｚ）^２は、次式で与えられる。

ここで、
（Ｓ_ｉ ^ＴＷＳ_ｉ）＝Ｅ_Ｐｉ ^２ …（１１）
とおき、
Ｓ_ｉ ^ＴＷＳ_ｊ／｛（Ｓ_ｉ ^ＴＷＳ_ｉ）^１／２（Ｓ_ｊ ^ＴＷＳ_ｊ）^１／２｝＝ＣＯＲ_ｉｊ
とおけば、ＣＯＲ_ｉｊは、ｉ番目の実験とｊ番目の実験との間の相関を表す。Ｅ_Ｐｉ ^２は、実験と計算で確認された計算値のパラメータに関する相対誤差の二乗である。

（１０）式および（１１）式から、
Ｓ_ｉ ^ＴＷＳ_ｊ＝ＣＯＲ_ｉｊ×｛（Ｓ_ｉ ^ＴＷＳ_ｉ）^１／２（Ｓ_ｊ ^ＴＷＳ_ｊ）^１／２｝
＝ＣＯＲ_ｉｊ×｛（Ｅ_Ｐｉ ^２）^１／２（Ｅ_Ｐｊ ^２）^１／２｝
＝ＣＯＲ_ｉｊ×｛（Ｅ_Ｐｉ ^２）（Ｅ_Ｐｊ ^２）｝^１／２ …（１２）
となる。よって、（１０）式および（１２）式から、目的の体系の相対誤差の二乗は、

と表わされる。このように、本実施の形態では、臨界実験などの模擬実験で得られたパラメータに起因する相対誤差を組み合わせて、目的の体系の計算値に関してパラメータに起因する相対誤差を評価することになる。

ｉ番目の実験とｊ番目の実験との間に相関がない場合には、ＣＯＲ_ｉｊは０となる。このとき、目的の体系の相対誤差の二乗は、

となる。ｉ番目の実験とｊ番目の実験との間に相関がほとんどない場合にも、（１４）式を用いて目的体系の相対誤差を求めてもよい。

目的とする体系（製品）を、モデル方程式を用いた計算をして設計する際、その計算値に含まれる誤差は、製品に品質の保証に直接関係する内容であり、この誤差を評価しておくこと極めて重要である。計算値に含まれる誤差は、計算理論やモデル方程式事態の誤差も含まれるが、設計の際に用いた数値（入力値）に含まれる誤差によって生ずる割合が最も大きい。特に原子力の分野では、この入力値として用いられる核データライブラリの値などの物理定数に不確かさがあり、材質や寸法以外の値から生ずる誤差が無視できない。そこで、これらの計算全体の品質の確認のために、実験が行われ、計算値と測定値との差から目的とする製品の計算誤差を確認したり、見積もっている。

実験を実施する際、その実験が目的とする体系をどこまで再現しているのかという価値判断が必要となる。従来、その判断には材質や寸法が似かよっているとか、目的の体系の一部分をそのまま再現すれば良いという判断がなされてきた。

本実施の形態では、その類似性を表す数値を模擬性（評価）因子として、次のように定義している。

実験と目的とする体系のそれぞれの計算値に対して、入力値のそれぞれが単位変化した場合に計算値が変化する割合を求める。これらの割合は、感度係数と呼ばれる数値である。感度係数は、入力値の種類数だけ存在するから、実験でも目的とする体系でも感度係数は入力値の種類を次元とするベクトルとなる。

模擬性（評価）因子は、実験と目的とする体系での感度係数ベクトルがどれだけ近いかを表す数値である。一般に、模擬性（評価）因子は、実験と目的とする体系の二つの感度係数ベクトルが構成する角度θに対して、ｃｏｓθを表す数値となっており、統計学では相関係数と呼ばれる数値である。これらのベクトルが重なっていれば、角度θは零（０）であるので、ｃｏｓθは１となる。模擬性（評価）因子は０から１の実数である。

本実施の形態では、実験での測定値に対する計算誤差を用いて、目的とする体系での誤差を推定している。このとき、誤差は、模擬性（評価）因子を最大するようにして推定されている。さらに、複数の実験を組み合わせた場合にも、模擬性（評価）因子を最大するようにして推定されている。

考え方としては、実験での感度係数ベクトルＳ_ｉを一次結合して、目的の体系の感度係数ベクトルＳ_Ｒにできるだけ近いものを作るという手法である。本実施の形態では、このとき入力値の不確かさの相対割合を表す共分散（誤差）行列Ｗを介してその演算を行うので、ベクトルだけの演算ではない。しかし、考え方としては、実験の感度係数ベクトルを線形和して目的の感度係数ベクトルに近づけるという手法である。そのとき得られた重み係数を掛けて実験で得られた誤差の和をとり、目的とする体系の誤差を推定する。この手法で、実験で得られた計算誤差を合成する場合、模擬性（評価）因子を最大にする条件下で合成されているため最適な値が得られる。

線形結合した感度係数ベクトルＳ_ｋを、目的とする体系の感度係数ベクトルＳ_Ｒにできる限り近づけるために、（ａ）ベクトルの長さ（大きさ）を同じにして、同時に、（ｂ）２つのベクトルのなす角度をできるかぎり０に近づけるように重み係数を決定することにより、模擬性（評価）因子を最大にしている。

線形結合定数α_ｉを重みとして、モデルへの入力値に関する相対誤差Ｅ_Ｐｉを合成して製品の性能確認のための数値解（計算値）に含まれる誤差を推定すれば、実験から得られた計算誤差の情報を最も有効に利用することができる。このように、本実施の形態では、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に含まれるパラメータに起因する誤差を、その対象を模擬した実験の結果を用いて定量的に推定することができる。

また、複数個の実験の誤差の情報を組み合わせて、目標の体系の誤差の情報を推定することができる。このため、目標の体系を一度に模擬した実験ができない場合でも、本質的な条件を分割して実験を行うことで、目標の体系の誤差の情報を推定することが可能となる。

さらに、複数の実験を組み合わせて、目的とする体系の計算誤差の割合を定量評価する手法をことができるため、実験の組み合わせ、すなわち個々の実験内容を最適化できる。たとえば、最終的に目的とする体系には３つの際立った特徴がある場合について考える。実験では実験の制約上、この３つの特徴を全て同時に実現できない場合、１つずつの特徴を取り入れた別々の３ケースの実験を行う。これら３つの実験の計算値と測定値の３つの飲み合わせを全て利用して、目的とする体系の着目する物理量に生ずる計算誤差を評価できる。

実験のデータを組み合わせることができるので、現在の実験データと過去の実験データを組み合わせて、実験の実施ケースを減らし、時間的にも経済的にも資源を節約することができる。なお、同じ実験装置で得られたデータのみを組み合わせるのではなくて、理論的には別の実験で得られたデータとの組み合わせも可能である。つまり、部分実験から最終的な製品・施設の設計精度を確認することができる。

また、この手法は、線型代数の簡単な数学的操作と行列演算しか用いないので、式の形式が非常に単純である。手法や数値の処理方法が単純なために、適用範囲が広い。

この手法で扱う誤差行列は、計算の入力となった全てのパラメータについて定義される。このため、臨界実験を考えると、たとえば断面積に起因する誤差のみでなく、数密度による誤差、物理定数（核分裂収率、崩壊定数、動特性パラメータなど：断面積はこのグループ）、長さ、温度、時間等の誤差も含めて一般的な取扱が可能である。

最終的に実験で得られたパラメータに関わる誤差を合成するため、計算のアルゴリズムで用いる相対誤差行列（共分散行列）は、相対誤差の絶対値を正しく評価するものではなく、相対誤差の相対値の関係を表したものや相関行列で計算の処理を行うことができる。つまり、相対誤差の割合の関係が得られれば評価可能であり、誤差の要因について互いの関係が数値的な割合で得られればよい。

これまで重要視されてこなかった原子数密度の計算誤差に起因する核特性の誤差についても評価が可能で、使用済みの燃料集合体の臨界安全性や、保管、輸送の評価に適用可能である。

着目する物理量Ｒが、たとえば分布系のベクトル量であっても、同様の方法で、パラメータに起因する誤差を定量的に推定することができる。模擬実験で得られるベクトル量の測定値をｂ_Ｅとする。また、モデルを用いてシミュレーションした結果得られるベクトル量の計算値をｂ_ｃとする。ｂ_Ｅおよびｂ_ｃはいずれもベクトルである。この場合には、物理量Ｒの測定値と計算値との比（Ｒ_Ｃ／Ｒ_Ｅ）は、次式によって求めればよい。

Ｒ_Ｃ／Ｒ_Ｅ＝｛（ｂ_ｃ・ｂ_Ｅ）／（｜ｂ_ｃ｜｜ｂ_Ｅ｜）｝
ここで、（ｂ_ｃ・ｂ_Ｅ）は内積を、｜ｂ_ｃ｜はｂ_ｃの大きさを表す。

このようにして得られた物理量Ｒの測定値と計算値との比（Ｒ_Ｃ／Ｒ_Ｅ）を用いて、実験体系のパラメータによって生じた相対誤差の二乗Ｅ_Ｐｉを求めることがきる。この相対誤差の二乗Ｅ_Ｐｉを用いて、上述の演算を行い、パラメータに起因する誤差を定量的に推定することができる。つまり、扱える物理量は、たとえば臨界固有値などのスカラー量だけではなく、分布系（ベクトル量）にも適用可能である。

また、目的とする体系と実験体系とで異なる共分散（誤差）行列を用いてもよい。目的とする体系の共分散行列をＷ_Ｒ、実験体系での共分散行列をＷ_Ｅとする。ただし、Ｗ_ＲとＷ_Ｅの構造は同じで、この場合、ふたつの行列の互いの（ｉ，ｊ）成分は同じ単位で定義されるものとする。すなわち、互いの（ｉ，ｊ）成分は、大きさの比較が直接行えるものとする。

この場合、（２）式および（３）式は、
Ｓ_ｋ ^ＴＷ_ＲＳ_ｋ＝Ｓ_Ｒ ^ＴＷ_ＲＳ_Ｒ …（２’）
および、
（Ｓ_Ｒ−Ｓ_ｋ）^ＴＷ_Ｒ（Ｓ_Ｒ−Ｓ_ｋ）＝０ …（３’）
となる。これらの式を用いて、線形結合定数α_ｉを求める。

さらに、次式が成り立つように、定数ｔ_ｉ（ｉ＝１，２，…，ｎ）を決定する。

ｔ_ｉＳ_ｉ ^ＴＷ_ＥＳ_ｉ＝Ｓ_ｉ ^ＴＷ_ＲＳ_ｉ
定数ｔ_ｉ（ｉ＝１，２，…，ｎ）の決定は、たとえば線形結合定数演算部１５で行えばよい。このようにして得られた定数ｔ_ｉを用いて、目的の体系の相対誤差の二乗は、

と求められる。

このように、目的とする体系の入力に関しての不確かさを表した共分散（誤差）行列とそれぞれの実験についての体系の入力に関しての不確かさを表した共分散（誤差）行列が同じ行列という仮定をせずに、この共分散行列が異なる場合でも処理すべき式を一般化できて計算誤差を推定することができる。

１０…計算誤差推定装置、１１…相対誤差記憶部、１２…共分散行列記憶部、１３…目的体系感度係数ベクトル記憶部、１４…実験体系感度係数ベクトル記憶部、１５…線形結合定数演算部、１６…誤差演算部、２１…コンピュータ、２２…ＣＰＵ、２３…ＲＡＭ、２４…ＲＯＭ、２５…入力装置、２６…入力制御部、２８…表示装置、２９…表示制御部、３０…ハードディスク、３１…ハードディスク制御部、３２…バス

Claims (5)

1. 対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に含まれる誤差を、その対象を模擬した実験の結果を用いて推定する誤差推定装置において、
   ｉを１から実験の総数ｎまでの自然数としたときにｉ番目の実験についての前記モデルを用いたシミュレーションで得られたある物理量Ｒの計算値の当該実験で測定された当該物理量Ｒの測定値に対する前記モデルへの入力値に関する相対誤差Ｅ_Ｐｉを記憶する相対誤差記憶部と、
   シミュレーションに用いる前記モデルへの入力値の不確かさの割合の関係を示す共分散行列Ｗを記憶する共分散行列記憶部と、
   前記対象について前記モデルを用いたシミュレーションの結果のそのモデルへの入力値の単位変化に対する変化量を示す目的体系感度係数ベクトルＳ_Ｒを記憶する目的体系感度係数ベクトル記憶部と、
   ｉ番目の実験の体系について前記モデルを用いたシミュレーションした結果のそのモデルへの入力値の単位変化に対する変化量を表す実験体系感度係数ベクトルＳ_ｉを記憶する実験体系感度係数ベクトル記憶部と、
   前記実験体系感度係数ベクトルＳ_ｉを線形結合した線形結合ベクトルＳ_ｋの線形結合定数α_ｉを、Ｓ_ｋ ^ＴＷＳ_ｋ＝Ｓ_Ｒ ^ＴＷＳ_Ｒを満足しかつ前記目的体系感度係数ベクトルＳ_Ｒと前記線形結合ベクトルＳ_ｋとのなす角が最小になるように求める線形結合定数演算部と、
   前記線形結合定数α_ｉを重みとして前記相対誤差Ｅ_Ｐｉの値を合成して前記モデルを用いたシミュレーションの結果得られる前記対象の物理量Ｚに含まれる相対誤差ΔＺ／Ｚを推定する誤差演算部と、
   を有することを特徴とする誤差推定装置。
2. 前記線形結合定数演算部は、（Ｓ_Ｒ−Ｓ_ｋ）^ＴＷ（Ｓ_Ｒ−Ｓ_ｋ）の絶対値を最小とすることにより前記目的体系感度係数ベクトルＳ_Ｒと前記線形結合ベクトルＳ_ｋとのなす角を最小にすることを特徴とする請求項１に記載の誤差推定装置。
3. 前記誤差演算部は、Ｓ_ｉ ^ＴＷＳ_ｊ／｛（Ｓ_ｉ ^ＴＷＳ_ｉ）^１／２（Ｓ_ｊ ^ＴＷＳ_ｊ）^１／２｝＝ＣＯＲ_ｉｊとして、
   

   

   により前記相対誤差ΔＺ／Ｚを推定することを特徴とする請求項１または請求項２に記載の誤差推定装置。
4. 前記共分散行列Ｗは、前記対象での共分散行列Ｗ_Ｒと前記実験の体系での共分散行列Ｗ_Ｅとを含み、
   前記線形結合定数演算部は、Ｓ_ｋ ^ＴＷ_ＲＳ_ｋ＝Ｓ_Ｒ ^ＴＷ_ＲＳ_Ｒを満足する線形結合定数α_ｉを求め、さらに、ｔ_ｉＳ_ｉ ^ＴＷ_ＥＳ_ｉ＝Ｓ_ｉ ^ＴＷ_ＲＳ_ｉを満足する定数ｔ_ｉを求め、
   前記誤差演算部は、Ｓ_ｉ ^ＴＷ_ＥＳ_ｊ／｛（Ｓ_ｉ ^ＴＷ_ＥＳ_ｉ）^１／２（Ｓ_ｊ ^ＴＷ_ＥＳ_ｊ）^１／２｝＝ＣＯＲ_ｉｊとして、
   

   

   により前記相対誤差ΔＺ／Ｚを推定することを特徴とする請求項１または請求項２に記載の誤差推定装置。
5. 対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果に含まれる誤差を、その対象を模擬した実験の結果を用いて推定する誤差推定プログラムにおいて、コンピュータに、
   ｉを１から実験の総数ｎまでの自然数としたときにｉ番目の実験についての前記モデルを用いたシミュレーションで得られたある物理量Ｒの計算値の当該実験で測定された当該物理量Ｒの測定値に対する前記モデルへの入力値に関する相対誤差Ｅ_Ｐｉを記憶する相対誤差記憶機能と、
   シミュレーションに用いる前記モデルへの入力値の不確かさの割合の関係を示す共分散行列Ｗを記憶する共分散行列記憶機能と、
   前記対象について前記モデルを用いたシミュレーションの結果のそのモデルへの入力値の単位変化に対する変化量を示す目的体系感度係数ベクトルＳ_Ｒを記憶する目的体系感度係数ベクトル記憶機能と、
   ｉ番目の実験の体系について前記モデルを用いたシミュレーションした結果のそのモデルへの入力値の単位変化に対する変化量を表す実験体系感度係数ベクトルＳ_ｉを記憶する実験体系感度係数ベクトル記憶機能と、
   前記実験体系感度係数ベクトルＳ_ｉを線形結合した線形結合ベクトルＳ_ｋの線形結合定数α_ｉを、Ｓ_ｋ ^ＴＷＳ_ｋ＝Ｓ_Ｒ ^ＴＷＳ_Ｒを満足しかつ前記目的体系感度係数ベクトルＳ_Ｒと前記線形結合ベクトルＳ_ｋとのなす角が最小になるように求める線形結合定数演算機能と、
   前記線形結合定数α_ｉを重みとして前記相対誤差Ｅ_Ｐｉの値を合成して前記モデルを用いたシミュレーションの結果得られる前記対象の物理量Ｚに含まれる相対誤差ΔＺ／Ｚを推定する誤差演算機能と、
   を実現させることを特徴とする誤差推定プログラム。

Images