(57)【要約】【課題】解析負荷を低減することのできる解析条件設定技術を提供する。【解決手段】解析条件設定装置１０は、予め設定された第１参照値ΔＷを登録する登録部１４と、解析値又はその積算量の時系列データ群１１を取得するデータ取得部１２と、時系列データ群１１の時間軸において時間幅を有する時間区間を積算量の差分値と第１参照値ΔＷに基づいて設定する時間区間設定部１５と、設定された時間区間に基づいて時系列データ群１１をプルームに分割するプルーム分割部１６と、を少なくとも備えている。【選択図】図１
2019169104.tif 000002
【特許請求の範囲】
【請求項１】
予め設定された第１参照値を登録する登録部と、
解析値又はその積算量の時系列データ群を取得するデータ取得部と、
前記時系列データ群の時間軸において時間幅を有する時間区間を、前記積算量の差分値と前記第１参照値に基づいて設定する時間区間設定部と、
設定された前記時間区間に基づいて、前記時系列データ群をプルームに分割するプルーム分割部と、を備えることを特徴とする解析条件設定装置。
【請求項２】
請求項１に記載の解析条件設定装置において、
前記登録部に、予め設定された第２参照値がさらに登録されており、
前記時間区間は、前記第２参照値よりも小さくならないように分割されることを特徴とする解析条件設定装置。
【請求項３】
請求項１又は請求項２に記載の解析条件設定装置において、
前記登録部に、予め設定された第３参照値がさらに登録されており、
前記時間区間は、前記第３参照値よりも大きくならないように分割されることを特徴とする解析条件設定装置。
【請求項４】
請求項１から請求項３のいずれか１項に記載の解析条件設定装置において、
前記データ取得部は、複数の前記時系列データ群を取得し、
前記時系列データ群の各々に設定された複数の時間区間は、相互に対応するもの同士の前記時間幅が等しくなるように設定されており、
プルーム分割部は、複数の前記時系列データ群を、相互に対応する時間区間ごとに、前記プルームに分割する、ことを特徴とする解析条件設定装置。
【請求項５】
請求項１から請求項４のいずれか１項に記載の解析条件設定装置と、
複数の前記プルームを単位として所定の解析コードで解析を実行する解析部と、
前記解析の結果を出力する出力部と、を備えることを特徴とする解析装置。
【請求項６】
予め設定された第１参照値を登録する登録ステップと、
解析値又はその積算量の時系列データ群を取得するデータ取得ステップと、
前記時系列データ群の時間軸において時間幅を有する時間区間を、前記積算量の差分値と前記第１参照値に基づいて設定する時間区間設定ステップと、
設定された前記時間区間に基づいて、前記時系列データ群をプルームに分割する分割ステップと、を含むことを特徴とする解析条件設定方法。
【請求項７】
コンピュータに、
予め設定された第１参照値を登録する登録ステップ、
解析値又はその積算量の時系列データ群を取得するデータ取得ステップ、
前記時系列データ群の時間軸において時間幅を有する時間区間を、前記積算量の差分値と前記第１参照値に基づいて設定する時間区間設定ステップ、
設定された前記時間区間に基づいて、前記時系列データ群を複数のプルームに分割する分割ステップ、を実行させることを特徴とする解析条件設定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、時系列データ群の解析条件を設定する技術に関する。
【背景技術】
【０００２】
原子力施設等で発生するあらゆる事故を対象として、その発生頻度と発生時の影響を定量評価し、その積であるリスクがどれほど小さいかで安全性の度合いを表現する確率論的リスク評価（ＰＲＡ：Probabilistic Risk Assessment）が活用されている。原子力発電所の事故におけるＰＲＡは、炉心損傷確率を解析するレベル１、格納容器破損確率及び放射性物質の施設外放出確率を解析するレベル２、健康影響確率を解析するレベル３に分けられる。
【０００３】
このレベル３ＰＲＡにおける不確かさ解析は、原子力発電所の事故による公衆のリスクを評価するものであって、１つの事故シナリオに対し、数百ケースのソースターム（放射性物質の放出条件）が設定されている。そして、この不確かさ解析では、年間の気象条件の変化に対応させるために例えば１時間ごとのケースを考慮するとなると、８７６０ケースの気象シーケンスに対し計算が行われることとなる。このため、解析ケースの数が、数十万から数百万ケースに上る。
【０００４】
さらに、レベル３ＰＲＡで使用する解析コードでは、ソースタームを複数のプルーム（放射性物質を含む雲）に分割し、各プルームについて大気中での移行・拡散を計算する。ソースタームは、事故シナリオによって様々であるため、プルームへの分割はソースタームの傾向を見ながら手動で行うか、単純に時間的に等間隔で分割するかのいずれかの方法が用いられている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】NUREG/CR-7155 “ State-of-the-Art Reactor Consequence Analyses Project, Uncertainty Analysis of the Unmitigated Long-Term Station Blackout of the Peach Bottom Atomic Power Station,” 2016
【発明の概要】
【発明が解決しようとする課題】
【０００６】
レベル２ＰＲＡで解析されレベル３ＰＲＡに受け渡されるソースタームは、放射性物質を物理的・化学的挙動からいくつかのグループに分類され、各グループの環境への放出量を時間ステップごとに表した時系列データ群で構成されたものが一般的である。
【０００７】
レベル３ＰＲＡでは、このような一連の放射性物質の放出を表した時系列データ群を複数のプルームに分割し、各プルームについて大気中の移行・拡散を計算する。前述したように、レベル３ＰＲＡの不確かさ解析では、１つの事故シナリオに対し数百ケースのソースタームについて解析を実施する。したがって、数百のソースタームに対し、プルーム分割の操作を行う必要があり、ソースタームの傾向を見ながら、手動で分割位置を決定することは現実的ではない。
【０００８】
プルームを自動的に分割する方法として、時間的な等間隔でプルームを分割する方法がある。しかしこの方法では、プルーム分割の間隔を例えば気象データの単位時間である１時間に設定すると、放射性物質の放出期間が長期にわたる事故シナリオにおいては、プルームの数が増加してしまう。すると、プルームの数が増加した分だけ解析負荷、即ち解析にかかる時間及び必要なメモリ容量が増加してしまう。
【０００９】
一方において、そのような時間的な等間隔でプルーム分割を粗く設定すると、プルームの数を減少させ解析負荷を低減させることができるが、解析結果の精度は低下し、公衆のリスクを正しく把握することが困難となる。さらに、一つの事故シナリオの中に複数の放出位置からの放射性物質の放出がある場合、放出位置毎にそれぞれプルームを設定することになり、その分だけプルームの数が増加してしまう。
【００１０】
本発明の実施形態はこのような事情を考慮してなされたもので、解析負荷を低減することのできる解析条件設定技術を提供することを目的とする。
【課題を解決するための手段】
【００１１】
実施形態に係る解析条件設定装置において、予め設定された第１参照値を登録する登録部と、解析値又はその積算量の時系列データ群を取得するデータ取得部と、前記時系列データ群の時間軸において時間幅を有する時間区間を、前記積算量の差分値と前記第１参照値に基づいて設定する時間区間設定部と、設定された前記時間区間に基づいて、前記時系列データ群をプルームに分割するプルーム分割部と、を備えることを特徴とする。
【００１２】
実施形態に係る解析装置において、前記解析条件設定装置と、複数の前記プルームを単位として所定の解析コードで解析を実行する解析部と、前記解析の結果を出力する出力部と、を備えることを特徴とする。
【発明の効果】
【００１３】
本発明の実施形態によれば、解析負荷を低減することのできる解析条件設定技術が提供される。
【図面の簡単な説明】
【００１４】
【図１】本発明の各実施形態に係る解析条件設定装置及び解析装置のブロック図。
【図２】第１実施形態に係る解析条件設定装置で処理される時系列データ群の立ち上がり部分の詳細を示すグラフ。
【図３】第１実施形態に係る解析条件設定装置で処理される時系列データ群の全体を示すグラフ。
【図４】第１実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャート。
【図５】第２実施形態に係る解析条件設定装置で処理される時系列データ群の立ち上がり部分の詳細を示すグラフ。
【図６】第２実施形態に係る解析条件設定装置で処理される時系列データ群の全体を示すグラフ。
【図７】第２実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャート。
【図８】第３実施形態に係る解析条件設定装置で処理される複数の時系列データ群を示すグラフ。
【図９】第３実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャート。
【図１０】第４実施形態に係る解析条件設定装置で処理される複数の時系列データ群を示すグラフ。
【発明を実施するための形態】
【００１５】
以下、本発明の実施形態を添付図面に基づいて説明する。図１は、本発明の各実施形態に係る解析条件設定装置及び解析装置のブロック図である。
【００１６】
（第１実施形態）
図２は第１実施形態に係る解析条件設定装置で処理される時系列データ群の立ち上がり部分の詳細を示すグラフである。図３は第１実施形態に係る解析条件設定装置で処理される時系列データ群の全体を示すグラフである。
【００１７】
図１に示すように、第１実施形態に係る解析条件設定装置１０は、予め設定された第１参照値ΔＷを登録する登録部１４と、解析値又はその積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝（図２，図３参照）の時系列データ群１１を取得するデータ取得部１２と、時系列データ群１１の時間軸（図２，図３参照）において時間幅が調節された複数の時間区間Ｚg（Ｚ1，Ｚ2，…ＺG）の少なくとも一つを、時間区間Ｚgの積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））（図４参照）と第１参照値ΔＷに基づいて設定する時間区間設定部１５と、設定された時間区間Ｚg（Ｚ1，Ｚ2，…ＺG）に基づいて時系列データ群１１を複数のプルームＰg｛Ｗ(ｔm),… Ｗ(ｔm+n-1)｝（ｇ＝１〜Ｇ）（図４参照）に分割するプルーム分割部１６と、を少なくとも備えている。
【００１８】
登録部１４に登録される第１参照値ΔＷは、その値自身が予め設定されている場合以外に、演算部１３において演算された値が設定される場合がある。この演算部１３は、時間範囲がｔ0からｔRである時系列データ群１１における積算量の最大値Ｗ（ｔＲ）（図３参照）を、指定した分割数で除算して第１参照値ΔＷを演算する。これにより、積算量の最大値Ｗ（ｔＲ）の大きさにかかわらず、時系列データ群１１は、分割数に対応する数のプルームＰg｛Ｗ(ｔm),… Ｗ(ｔm+n-1)｝（ｇ＝１〜Ｇ）（図４参照）に適切に分割される。
【００１９】
上述の通り、時間区間設定部１５は、時系列データ群１１の時間軸（図２，図３参照）において時間幅が調節された複数の時間区間Ｚg（Ｚ1，Ｚ2，…ＺG）の少なくとも一つを、積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm）など）（図４参照）と第１参照値ΔＷに基づいて、すなわち時間区間Ｚgの積算量の差分値と第１参照値ΔＷを参照して設定する。
【００２０】
一般に解析値又はその積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝は時間的な離散データとして得られることが多いことから、後述するように、時間区分設定部１５は、ｎを１ずつ増やしながら積算量の差分値を順次計算していき、差分値（Ｗ（ｔm+n）−Ｗ（ｔm））が第１参照値ΔＷを超えた場合にＷ（ｔm）からＷ（ｔm+n-1）となる区間を時間区間ＺgのプルームＰg｛Ｗ(ｔm),… Ｗ(ｔm+n-1)｝（ｇ＝１〜Ｇ）として設定する。換言するとこの場合、時間区分設定部１５は、少なくとも一つの時間区間Ｚgを、積分量の差分値（Ｗ（ｔm+n-1）−Ｗ（ｔm））が第１参照値ΔＷ以下の最大値となるように設定している。
【００２１】
さらに、解析値又はその積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝は時間的な離散データの量が十分に多い場合などには、第１参照値ΔＷに予め定めた一定幅を持たせた構成としても構わない。
【００２２】
図２に示した例では、時間的な離散データとして得られた解析値又はその積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝を内挿または外挿して連続データとしての近時曲線（関数）を定め、時間区分設定部１５は時間区間Ｚgを、この連続データである近似曲線に基づいて、各々の時間区間Ｚgにおける積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））（図４参照）が予め定めた一定値である第１参照値ΔＷとなるように設定している。
【００２３】
レベル３ＰＲＡでは、放出点からの距離毎の早期被ばく線量、長期被ばく線量、被ばく経路（クラウドシャイン、グランドシャイン、吸入、摂取、再浮遊物質吸入）ごとの線量など、目的に応じて多数の項目の評価が可能である。また、気象条件の変動を考慮し、年間気象変化に対する平均値、５０％値、９５％値、最大値などの統計値を出力する。ソースターム条件と評価対象項目、気象統計値などによって、分割数の指定値と解析精度及びプルーム単体総数の関係が異なるため、ユーザーがソースタームの種類とレベル３ＰＲＡでの評価対象項目、気象統計値を選択すると、分割数と解析精度への影響及びプルーム単体数が出力されるデータベースを備え、ユーザーをサポートする。データベースはあらかじめ代表的ソースタームについて本発明の方法でプルーム分割を行い、レベル３ＰＲＡを実施した結果に基づき作成する。解析精度については、例えばプルーム単体の放出継続時間を気象データの単位時間である１時間で時間的に等分割したものをベースケースとし、分割数の値を変えて本発明の手法で分割した場合とベースケースとのレベル３ＰＲＡ評価結果との誤差をデータベースとして保有しておく方法がある。これにより、ユーザーが解析の目的に応じ、解析精度と解析負荷のトレードオフの関係から適切なプルーム分割を行うための分割数を指定することができる。
【００２４】
データ取得部１２で取得される時系列データ群１１は、ＰＲＡのレベル２から引き渡されるソースタームである。このソースタームは、一般に放射性物質を物理的・化学的性質等の条件によっていくつかのグループに分類されており、それぞれのグループ毎に環境への放出量あるいは初期内蔵量に対する放出割合を時系列データ群１１として出力したものである。第１実施形態では、放射性物質のグループが一つの場合を取り上げている。
【００２５】
この時系列データ群１１は、単位時間当たりの変化量を示す解析値であったり、この解析値が予め時間積分された積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)であったりする。時系列データ群１１がそのような解析値である場合は、データ取得部１２の後段において、適時、積分処理が行われる。
【００２６】
時間区間設定部１５で設定した複数の時間区間Ｚg（Ｚ1，Ｚ2，…ＺG）の各々に基づいて、時系列データ群１１が分割される。この時系列データ群１１を分割して得られた複数のプルームＰg｛Ｗ(ｔm),… Ｗ(ｔm+n-1)｝（ｇ＝１〜Ｇ）は、蓄積部１７に蓄積される。なお表示部１８において、複数のプルームＰgに分割された時系列データ群１１のグラフを表示させることができる。
【００２７】
さらに解析装置２０は、上述した解析条件設定装置１０と、複数のプルームＰg（ｇ＝１〜Ｇ）を単位として所定の解析コードで解析を実行する解析部２１と、この解析の結果を出力する出力部２３と、を少なくとも備えている。ここで、解析部２１とは、例えば１時間おきに情報提供される気象情報に基づいて、放出積算量の変化が時系列データ群１１として与えられている放射性物質の拡散シミュレーションを行うものである。
【００２８】
図４は第１実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャートである。この図４に基づいて図２及び図３におけるデータ処理を説明する。
まず、予め設定された第１参照値ΔＷが登録される（Ｓ１１）。次に積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝の時系列データ群１１を取得する（Ｓ１２）。なお取得される時系列データ群１１は、積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝を例示しているが、この積算量の微分値に相当する解析値である場合もある。
【００２９】
次に図２を参照しつつ、ｍ＝０（Ｓ１３）、ｎ＝１（Ｓ１４）として、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｎｏ）。引き続き、ｍ＝０、ｎ＝２（Ｓ２２）として、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｙｅｓ）。その結果、[ｔ0,ｔ1]の時間区間Ｚ1が設定され、この時間区間Ｚ1に割り付けられている｛Ｗ(ｔ0),Ｗ(ｔ1)｝がプルームＰ1として分割される（Ｓ２３）。
【００３０】
引き続いて、ｍ＝２（Ｓ２４ Ｎｏ，Ｓ２５）、ｎ＝１（Ｓ１４）として、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｎｏ）。さらに、ｍ＝２、ｎ＝２（Ｓ２２）として、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８Ｙｅｓ）。その結果、[ｔ2,ｔ3]の時間区間Ｚ2が設定され、この時間区間Ｚ2に割り付けられている｛Ｗ(ｔ2),Ｗ(ｔ3)｝がプルームＰ2として分割される（Ｓ２３）。
【００３１】
引き続いて、ｍ＝４（Ｓ２４ Ｎｏ，Ｓ２５）、ｎ＝１，２として、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｎｏ）。さらに、ｍ＝４、ｎ＝３（Ｓ２２）として、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８Ｙｅｓ）。その結果、[ｔ4,ｔ6]の時間区間Ｚ3が設定され、この時間区間Ｚ3に割り付けられている｛Ｗ(ｔ4),Ｗ(ｔ5),Ｗ(ｔ6)｝がプルームＰ3として分割される（Ｓ２３）。
【００３２】
引き続いて、ｍ＝７（Ｓ２４ Ｎｏ，Ｓ２５）、ｎ＝１〜１２として、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｎｏ）。さらに、ｍ＝７、ｎ＝１３（Ｓ２２）として、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８Ｙｅｓ）。その結果、[ｔ7,ｔ19]の時間区間Ｚ4が設定され、この時間区間Ｚ4に割り付けられている｛Ｗ(ｔ7),…,Ｗ(ｔ19)｝がプルームＰ4として分割される（Ｓ２３）。
【００３３】
上述したフローが、最後の積分量Ｗ(ｔR)まで繰り返され、時系列データ群１１は、複数のプルームＰg（ｇ＝１〜Ｇ）に分割されることになる。これにより、プルーム分割が粗く設定されることなく、解析精度と解析負荷の適切な関係を維持することができる。
【００３４】
（第２実施形態）
次に図１、図５、図６、図７を参照して本発明における第２実施形態について説明する。なお、図１に示される第２実施形態の解析条件設定装置１０において、第１実施形態と重複する構成については説明を省略する。
【００３５】
第２実施形態の解析条件設定装置１０の登録部１４には、第１参照値ΔＷに加え、予め設定された第２参照値Δｔmin及び第３参照値Δｔmaxがさらに登録されている。そして、時間区間Ｚg（Ｚ1，Ｚ2，…ＺG）は、第２参照値Δｔminよりも小さくならないように分割され、さらに第３参照値Δｔmaxよりも大きくならないように分割される。ここで、第２参照値Δｔminの設定は、解析部２１で利用する気象情報等のデータの更新間隔（例えば一時間程度）が一つの目安となる。
【００３６】
図５は第２実施形態に係る解析条件設定装置で処理される時系列データ群の立ち上がり部分の詳細を示すグラフである。図６は第２実施形態に係る解析条件設定装置で処理される時系列データ群の全体を示すグラフである。図７は第２実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャートである。この図７に基づいて図５及び図６におけるデータ処理を説明する。
【００３７】
まず、予め設定された第１参照値ΔＷ、第２参照値Δｔmin及び第３参照値Δｔmaxが登録される（Ｓ１１）。次に積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝の時系列データ群１１を取得する（Ｓ１２）。なお取得される時系列データ群１１は、積算量Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝を例示しているが、この積算量の微分値に相当する解析値である場合もある。
【００３８】
次に図５を参照しつつ、ｍ＝０（Ｓ１３）、ｎ＝１，２（Ｓ１４，Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｎｏ）。そして、ｍ＝０（Ｓ１３）、ｎ＝３（Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｙｅｓ）。
【００３９】
ｍ＝０，ｎ＝３のまま引き続いて、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｙｅｓ）。その結果、[ｔ0,ｔ2]の時間区間Ｚ1が設定され、この時間区間Ｚ1に割り付けられている｛Ｗ(ｔ0),Ｗ(ｔ1),Ｗ(ｔ2)｝がプルームＰ1として分割される（Ｓ２３）。
【００４０】
引き続いて、ｍ＝３（Ｓ２４ Ｎｏ，Ｓ２５）、ｎ＝１，２（Ｓ１４，Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｎｏ）。そして、ｍ＝３のままｎ＝３（Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｙｅｓ）。
【００４１】
ｍ＝３，ｎ＝３のまま引き続いて、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｙｅｓ）。その結果、[ｔ3,ｔ5]の時間区間Ｚ2が設定され、この時間区間Ｚ2に割り付けられている｛Ｗ(ｔ3),Ｗ(ｔ4),Ｗ(ｔ5)｝がプルームＰ2として分割される（Ｓ２３）。
【００４２】
引き続いて、ｍ＝６（Ｓ２４ Ｎｏ，Ｓ２５）、ｎ＝１（Ｓ１４）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｎｏ）。そして、ｍ＝６のままｎ＝２（Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｙｅｓ）。
【００４３】
ｍ＝６，ｎ＝２のまま引き続いて、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｎｏ）。さらに、ｍ＝６のままｎ＝３に更新し（Ｓ２１ Ｎｏ，Ｓ２２）、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８Ｙｅｓ）。その結果、[ｔ6,ｔ8]の時間区間Ｚ3が設定され、この時間区間Ｚ3に割り付けられている｛Ｗ(ｔ6),Ｗ(ｔ7),Ｗ(ｔ8)｝がプルームＰ3として分割される（Ｓ２３）。
【００４４】
引き続いて、ｍ＝９（Ｓ２４ Ｎｏ，Ｓ２５）、ｎ＝１，２（Ｓ１４，Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｎｏ）。そして、ｍ＝９のままｎ＝３（Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｙｅｓ）。
【００４５】
ｍ＝９のまま引き継いで、ｎ＝４〜６（Ｓ２２）とし、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｎｏ）。さらに、ｍ＝９のままｎ＝７に更新し（Ｓ２１ Ｎｏ，Ｓ２２）、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８Ｙｅｓ）。その結果、[ｔ9,ｔ15]の時間区間Ｚ4が設定され、この時間区間Ｚ4に割り付けられている｛Ｗ(ｔ9),…,Ｗ(ｔ15)｝がプルームＰ4として分割される（Ｓ２３）。
【００４６】
引き続いて、ｍ＝１６（Ｓ２４ Ｎｏ，Ｓ２５）、ｎ＝１，２（Ｓ１４，Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｎｏ）。そして、ｍ＝１６のままｎ＝３（Ｓ１６）として、第２参照値Δｔminと時間幅（ｔm+n−ｔm）との対比を行う（Ｓ１５ Ｙｅｓ）。
【００４７】
ｍ＝１６のまま引き継いで、ｎ＝４〜９（Ｓ２２）とし、第１参照値ΔＷと積算量の差分値（Ｗ（ｔm+n）−Ｗ（ｔm））との対比を行う（Ｓ１８ Ｎｏ，Ｓ２１ Ｎｏ）。さらに、ｍ＝１６のままｎ＝１０に更新すると、第３参照値Δｔmaxと時間幅（ｔm+n−ｔm）との対比が行われる（Ｓ２１ Ｙｅｓ）。その結果、[ｔ16,ｔ25]の時間区間Ｚ5が設定され、この時間区間Ｚ5に割り付けられている｛Ｗ(ｔ16),…,Ｗ(ｔ25)｝がプルームＰ5として分割される（Ｓ２３）。
【００４８】
上述したフローが、最後の積分量Ｗ(ｔR)まで繰り返され、時系列データ群１１は、複数のプルームＰg（ｇ＝１〜Ｇ）に分割されることになる。これにより、プルーム分割が粗く設定されることなく、かつプルームの数が多すぎることもなく、解析精度と解析負荷の適切な関係を維持することができる。
【００４９】
（第３実施形態）
次に図１、図８及び図９を参照して本発明における第３実施形態について説明する。なお、図１に示される第３実施形態の解析条件設定装置１０において、第１実施形態及び第２実施形態と重複する構成については説明を省略する。
【００５０】
第３実施形態に係る解析条件設定装置１０（図１）において、データ取得部１２は、複数の時系列データ群Ｗ1(ｔr),Ｗ2(ｔr),…,ＷQ(ｔr)（図８参照）を取得する。そして、これら時系列データ群Ｗ1(ｔr),Ｗ2(ｔr),…,ＷQ(ｔr)の各々に設定された複数の時間区間Ｚg（Ｚ1，Ｚ2，…ＺG）は、相互に対応するもの同士の時間幅が等しくなるように設定されている。さらにプルーム分割部１６は、複数の時系列データ群Ｗ1(ｔr),Ｗ2(ｔr),…,ＷQ(ｔr)を、相互に対応する時間区間Ｚg（Ｚ1，Ｚ2，…ＺG）ごとに、プルームＰg（Ｐ1，Ｐ2，…ＰG）に分割する。
【００５１】
図８は第３実施形態に係る解析条件設定装置で処理される複数の時系列データ群を示すグラフである。図９は第３実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャートである。この図９は、図８から追加された（Ｓ１７）、（Ｓ１９）、（Ｓ２０）のフローに基づいて、異なる複数の時系列データ群Ｗ1(ｔr),Ｗ2(ｔr),…,ＷQ(ｔr)は、相互に対応する時間区間Ｚg（Ｚ1，Ｚ2，…ＺG）が一致し、共通のプルームＰg（Ｐ1，Ｐ2，…ＰG）に分割することができる。
【００５２】
複数の時系列データ群Ｗ1(ｔr),Ｗ2(ｔr),…,ＷQ(ｔr)を用いることにより、放射性物質がいくつかのグループに分類されている場合について模擬することができる。環境に放出された放射性物質による被ばく線量評価で重要な放射性物質のグループは、希ガスを代表するグループ、よう素を代表するグループ、セシウムを代表するグループである。
【００５３】
希ガスは放射性雲からの外部被ばくへの寄与が大きく、フィルターで除去されにくいことから、フィルターベント実施時において被ばく線量に支配的となる。よう素は、吸入した場合の内部被ばく線量に対する寄与が大きい。セシウムは地表沈着による汚染地面からの長期外部被ばく線量に対し支配的な核種である。これら３つのグループに着目してプルーム分割を行うことで、主な被ばく経路の線量を適切に評価する分割が可能となる。
【００５４】
（第４実施形態）
図１０は第４実施形態に係る解析条件設定装置で処理される複数の時系列データ群を示すグラフである。第４実施形態に係る解析装置２０は、第１位置に由来する時系列データ群Ｗ(ｔr)に基づいて解析した第１解析結果と第２位置に由来する時系列データ群Ｖ(ｔr)に基づいて解析した第２解析結果とを蓄積する蓄積部２２を備えている。さらに、これら第１解析結果及び第２解析結果を加算して出力させる加算部２４を、さらに備えている。
【００５５】
第１位置に由来する時系列データ群Ｗ(ｔr)及び第２位置に由来する時系列データ群Ｖ(ｔr)は、開始時刻も異なるために、それぞれの体系で設定された時間区間Ｚ1，Ｚ2，…Ｚ4及び時間区間Ｘ1，Ｘ2，Ｘ3に基づいてプルーム分割されている。しかし、それぞれの時系列データ群Ｗ(ｔr)，Ｖ(ｔr)から導かれた第１解析結果及び第２解析結果は、時間軸をあわせて互いを加算して評価することができる。
【００５６】
事故時の放射性物質放出の代表的なシナリオでは、事故進展に伴い、初めに建屋からの放射性物質の漏洩があり、その後、ベントスタックからフィルターベントが行われるケースがある。この場合には、放射性物質の放出位置が建屋とベントスタックの２か所となり、それぞれ、放射性物質の放出高さ、初期拡散幅が異なるため、レベル３ＰＲＡでは、これらを別々のプルームとして扱う。また、第１位置及び第２位置のそれぞれにおいてグループが一つの場合を取り上げているが、グループが複数の場合も適用することができる。
【００５７】
また、第４実施形態において、演算部１３により第１参照値ΔＷを演算する場合は、第１位置及び第２位置のそれぞれにおいて、異なる分割数を指定することができる。事故時の放射性物質放出の代表的なシナリオでは、事故進展に伴い、初めに建屋からの放射性物質の漏洩があり、その後、ベントスタックからフィルターベントが行われるケースがある。放出形態の違いにより、適切な分割数の値は異なるため、放出位置毎にこの分割数を指定することで、より最適なプルーム分割が可能となる。
【００５８】
以上述べた少なくともひとつの実施形態の解析条件設定装置によれば、積算量の差分値が参照値となるように時系列データ群をプルーム分割することにより、解析精度と解析負荷のトレードオフの関係を考慮しつつ、解析負荷を低減しかつ現実的な計算時間で評価することが可能となる。
【００５９】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。また、解析条件設定装置の構成要素は、コンピュータのプロセッサで実現することも可能であり、解析条件設定プログラムにより動作させることが可能である。
【００６０】
以上説明した解析条件設定装置は、専用のチップ、ＦＰＧＡ（Field Programmable Gate Array）、ＧＰＵ（Graphics Processing Unit）、又はＣＰＵ（Central Processing Unit）などのプロセッサを高集積化させた制御装置と、ＲＯＭ（Read Only Memory）やＲＡＭ（Random Access Memory）などの記憶装置と、ＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）などの外部記憶装置と、ディスプレイなどの表示装置と、マウスやキーボードなどの入力装置と、通信Ｉ／Ｆとを、備えており、通常のコンピュータを利用したハードウェア構成で実現できる。
【００６１】
また解析条件設定装置で実行されるプログラムは、ＲＯＭ等に予め組み込んで提供される。もしくは、このプログラムは、インストール可能な形式又は実行可能な形式のファイルでＣＤ−ＲＯＭ、ＣＤ−Ｒ、メモリカード、ＤＶＤ、フレキシブルディスク（ＦＤ）等のコンピュータで読み取り可能な記憶媒体に記憶されて提供するようにしてもよい。
【００６２】
また、本実施形態に係る解析条件設定装置で実行されるプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせて提供するようにしてもよい。また、装置１０は、構成要素の各機能を独立して発揮する別々のモジュールを、ネットワーク又は専用線で相互に接続し、組み合わせて構成することもできる。
【符号の説明】
【００６３】
１０…解析条件設定装置、１１…時系列データ群、１２…データ取得部、１３…演算部、１４…登録部、１５…時間区間設定部、１６…プルーム分割部、１７…蓄積部、１８…表示部、２０…解析装置、２１…解析部、２２…蓄積部、２３…出力部、２４…加算部、Ｗ(ｔr)｛Ｗ(ｔ0),Ｗ(ｔ1),…Ｗ(ｔR)｝…積算量、Ｚg（Ｚ1，Ｚ2，…ＺG）…時間区間、Ｗ（ｔm+n）−Ｗ（ｔm）…差分値、ΔＷ…第１参照値、Ｐg｛Ｗ(ｔm),… Ｗ(ｔm+n-1)｝…プルーム。
【図１】
2019169104.tif 000003
【図２】
2019169104.tif 000004
【図３】
2019169104.tif 000005
【図４】
2019169104.tif 000006
【図５】
2019169104.tif 000007
【図６】
2019169104.tif 000008
【図７】
2019169104.tif 000009
【図８】
2019169104.tif 000010
【図９】
2019169104.tif 000011
【図１０】
2019169104.tif 000012