JP2019169104A - Analysis condition setting device, method, and program, and analysis device - Google Patents

Analysis condition setting device, method, and program, and analysis device Download PDF

Info

Publication number
JP2019169104A
JP2019169104A JP2018058684A JP2018058684A JP2019169104A JP 2019169104 A JP2019169104 A JP 2019169104A JP 2018058684 A JP2018058684 A JP 2018058684A JP 2018058684 A JP2018058684 A JP 2018058684A JP 2019169104 A JP2019169104 A JP 2019169104A
Authority
JP
Japan
Prior art keywords
time
analysis
series data
reference value
condition setting
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.)
Granted
Application number
JP2018058684A
Other languages
Japanese (ja)
Other versions
JP6884724B2 (en
Inventor
美香 田原
Mika Tawara
美香 田原
裕貴 鈴木
Yuki Suzuki
裕貴 鈴木
智 柳下
Satoshi Yagishita
智 柳下
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Toshiba Energy Systems and Solutions Corp
Original Assignee
Toshiba Corp
Toshiba Energy Systems and Solutions Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Toshiba Energy Systems and Solutions Corp filed Critical Toshiba Corp
Priority to JP2018058684A priority Critical patent/JP6884724B2/en
Publication of JP2019169104A publication Critical patent/JP2019169104A/en
Application granted granted Critical
Publication of JP6884724B2 publication Critical patent/JP6884724B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

To provide an analysis condition setting technique capable of reducing an analysis load.SOLUTION: An analysis condition setting device 10 at least includes: a register unit 14 for registering a previously set first reference value ΔW; a data acquisition unit 12 for acquiring a time series data group 11 of analytic values or the integrated quantity of them; a time section setting unit 15 for setting a time section having a time width on the temporal axis of the time series data group 11, on the basis of the differential value of the integrated quantity and the first reference value ΔW; and a plume division unit 16 for dividing the time series data group 11 into plumes on the basis of the set time section.SELECTED DRAWING: Figure 1

Description

本発明の実施形態は、時系列データ群の解析条件を設定する技術に関する。   Embodiments described herein relate generally to a technique for setting analysis conditions for a time-series data group.

原子力施設等で発生するあらゆる事故を対象として、その発生頻度と発生時の影響を定量評価し、その積であるリスクがどれほど小さいかで安全性の度合いを表現する確率論的リスク評価(PRA:Probabilistic Risk Assessment)が活用されている。原子力発電所の事故におけるPRAは、炉心損傷確率を解析するレベル1、格納容器破損確率及び放射性物質の施設外放出確率を解析するレベル2、健康影響確率を解析するレベル3に分けられる。   Probabilistic risk assessment (PRA) for all accidents that occur in nuclear facilities, etc., by quantitatively evaluating the frequency of occurrence and the impact at the time of occurrence, and expressing the degree of safety based on how small the product risk is Probabilistic Risk Assessment) is being used. PRA in an accident at a nuclear power plant is divided into level 1 for analyzing the core damage probability, level 2 for analyzing the containment vessel failure probability and the out-of-facility release probability of radioactive materials, and level 3 for analyzing the health effect probability.

このレベル3PRAにおける不確かさ解析は、原子力発電所の事故による公衆のリスクを評価するものであって、1つの事故シナリオに対し、数百ケースのソースターム(放射性物質の放出条件)が設定されている。そして、この不確かさ解析では、年間の気象条件の変化に対応させるために例えば1時間ごとのケースを考慮するとなると、8760ケースの気象シーケンスに対し計算が行われることとなる。このため、解析ケースの数が、数十万から数百万ケースに上る。   This uncertainty analysis in Level 3 PRA evaluates the public risk from accidents at nuclear power plants, and hundreds of cases of source terms (radioactive material release conditions) are set for one accident scenario. Yes. In this uncertainty analysis, for example, when an hourly case is considered in order to cope with a change in annual weather conditions, calculation is performed on a weather sequence of 8760 cases. For this reason, the number of analysis cases increases from hundreds of thousands to millions of cases.

さらに、レベル3PRAで使用する解析コードでは、ソースタームを複数のプルーム(放射性物質を含む雲)に分割し、各プルームについて大気中での移行・拡散を計算する。ソースタームは、事故シナリオによって様々であるため、プルームへの分割はソースタームの傾向を見ながら手動で行うか、単純に時間的に等間隔で分割するかのいずれかの方法が用いられている。   Furthermore, in the analysis code used in Level 3 PRA, the source term is divided into a plurality of plumes (clouds containing radioactive materials), and the migration / diffusion in the atmosphere is calculated for each plume. Since source terms vary depending on the accident scenario, division into plumes is either done manually while looking at the trends of the source terms, or simply divided at regular intervals in time. .

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,” 2016NUREG / CR-7155 “State-of-the-Art Reactor Consequence Analyzes Project, Uncertainty Analysis of the Unmitigated Long-Term Station Blackout of the Peach Bottom Atomic Power Station,” 2016

レベル2PRAで解析されレベル3PRAに受け渡されるソースタームは、放射性物質を物理的・化学的挙動からいくつかのグループに分類され、各グループの環境への放出量を時間ステップごとに表した時系列データ群で構成されたものが一般的である。   The source terms analyzed at Level 2PRA and delivered to Level 3PRA are classified into several groups based on physical and chemical behavior of radioactive substances, and the time series showing the release amount of each group into the environment at each time step. Generally composed of data groups.

レベル3PRAでは、このような一連の放射性物質の放出を表した時系列データ群を複数のプルームに分割し、各プルームについて大気中の移行・拡散を計算する。前述したように、レベル3PRAの不確かさ解析では、1つの事故シナリオに対し数百ケースのソースタームについて解析を実施する。したがって、数百のソースタームに対し、プルーム分割の操作を行う必要があり、ソースタームの傾向を見ながら、手動で分割位置を決定することは現実的ではない。   In level 3PRA, a time series data group representing such a series of radioactive substance release is divided into a plurality of plumes, and the migration / diffusion in the atmosphere is calculated for each plume. As described above, in the level 3 PRA uncertainty analysis, analysis is performed for several hundred cases of source terms for one accident scenario. Therefore, it is necessary to perform a plume division operation on several hundred source terms, and it is not realistic to manually determine the division position while observing the tendency of the source terms.

プルームを自動的に分割する方法として、時間的な等間隔でプルームを分割する方法がある。しかしこの方法では、プルーム分割の間隔を例えば気象データの単位時間である1時間に設定すると、放射性物質の放出期間が長期にわたる事故シナリオにおいては、プルームの数が増加してしまう。すると、プルームの数が増加した分だけ解析負荷、即ち解析にかかる時間及び必要なメモリ容量が増加してしまう。   As a method of automatically dividing the plume, there is a method of dividing the plume at regular time intervals. However, in this method, if the interval of the plume division is set to, for example, 1 hour which is a unit time of weather data, the number of plumes increases in an accident scenario in which the radioactive material release period is long. Then, the analysis load, that is, the time required for the analysis and the necessary memory capacity are increased by the increase in the number of plumes.

一方において、そのような時間的な等間隔でプルーム分割を粗く設定すると、プルームの数を減少させ解析負荷を低減させることができるが、解析結果の精度は低下し、公衆のリスクを正しく把握することが困難となる。さらに、一つの事故シナリオの中に複数の放出位置からの放射性物質の放出がある場合、放出位置毎にそれぞれプルームを設定することになり、その分だけプルームの数が増加してしまう。   On the other hand, if the plume division is set roughly at such regular intervals, the number of plumes can be reduced and the analysis load can be reduced. However, the accuracy of the analysis results is reduced and the public risk is correctly grasped. It becomes difficult. Furthermore, when there is a release of radioactive material from a plurality of release positions in one accident scenario, a plume is set for each release position, and the number of plumes increases accordingly.

本発明の実施形態はこのような事情を考慮してなされたもので、解析負荷を低減することのできる解析条件設定技術を提供することを目的とする。   The embodiment of the present invention has been made in view of such circumstances, and an object thereof is to provide an analysis condition setting technique capable of reducing an analysis load.

実施形態に係る解析条件設定装置において、予め設定された第1参照値を登録する登録部と、解析値又はその積算量の時系列データ群を取得するデータ取得部と、前記時系列データ群の時間軸において時間幅を有する時間区間を、前記積算量の差分値と前記第1参照値に基づいて設定する時間区間設定部と、設定された前記時間区間に基づいて、前記時系列データ群をプルームに分割するプルーム分割部と、を備えることを特徴とする。   In the analysis condition setting device according to the embodiment, a registration unit for registering a preset first reference value, a data acquisition unit for acquiring a time-series data group of analysis values or integrated values thereof, and the time-series data group A time interval setting unit for setting a time interval having a time width on the time axis based on the difference value of the integrated amount and the first reference value, and the time series data group based on the set time interval And a plume dividing unit that divides the plume.

実施形態に係る解析装置において、前記解析条件設定装置と、複数の前記プルームを単位として所定の解析コードで解析を実行する解析部と、前記解析の結果を出力する出力部と、を備えることを特徴とする。   The analysis apparatus according to the embodiment includes: the analysis condition setting device; an analysis unit that performs analysis with a predetermined analysis code in units of the plurality of plumes; and an output unit that outputs the analysis result. Features.

本発明の実施形態によれば、解析負荷を低減することのできる解析条件設定技術が提供される。   According to the embodiment of the present invention, an analysis condition setting technique capable of reducing an analysis load is provided.

本発明の各実施形態に係る解析条件設定装置及び解析装置のブロック図。The block diagram of the analysis condition setting apparatus and analysis apparatus which concern on each embodiment of this invention. 第1実施形態に係る解析条件設定装置で処理される時系列データ群の立ち上がり部分の詳細を示すグラフ。The graph which shows the detail of the rising part of the time series data group processed with the analysis condition setting apparatus which concerns on 1st Embodiment. 第1実施形態に係る解析条件設定装置で処理される時系列データ群の全体を示すグラフ。The graph which shows the whole time series data group processed with the analysis condition setting apparatus which concerns on 1st Embodiment. 第1実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャート。The flowchart explaining the analysis condition setting method or analysis condition setting program which concerns on 1st Embodiment. 第2実施形態に係る解析条件設定装置で処理される時系列データ群の立ち上がり部分の詳細を示すグラフ。The graph which shows the detail of the rising part of the time series data group processed with the analysis condition setting apparatus which concerns on 2nd Embodiment. 第2実施形態に係る解析条件設定装置で処理される時系列データ群の全体を示すグラフ。The graph which shows the whole time series data group processed with the analysis condition setting apparatus which concerns on 2nd Embodiment. 第2実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャート。The flowchart explaining the analysis condition setting method or analysis condition setting program which concerns on 2nd Embodiment. 第3実施形態に係る解析条件設定装置で処理される複数の時系列データ群を示すグラフ。The graph which shows the several time series data group processed with the analysis condition setting apparatus which concerns on 3rd Embodiment. 第3実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャート。The flowchart explaining the analysis condition setting method or analysis condition setting program which concerns on 3rd Embodiment. 第4実施形態に係る解析条件設定装置で処理される複数の時系列データ群を示すグラフ。The graph which shows the several time series data group processed with the analysis condition setting apparatus which concerns on 4th Embodiment.

以下、本発明の実施形態を添付図面に基づいて説明する。図1は、本発明の各実施形態に係る解析条件設定装置及び解析装置のブロック図である。   Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings. FIG. 1 is a block diagram of an analysis condition setting device and an analysis device according to each embodiment of the present invention.

(第1実施形態)
図2は第1実施形態に係る解析条件設定装置で処理される時系列データ群の立ち上がり部分の詳細を示すグラフである。図3は第1実施形態に係る解析条件設定装置で処理される時系列データ群の全体を示すグラフである。
(First embodiment)
FIG. 2 is a graph showing details of a rising portion of a time-series data group processed by the analysis condition setting device according to the first embodiment. FIG. 3 is a graph showing the entire time-series data group processed by the analysis condition setting device according to the first embodiment.

図1に示すように、第1実施形態に係る解析条件設定装置10は、予め設定された第1参照値ΔWを登録する登録部14と、解析値又はその積算量W(tr){W(t0),W(t1),…W(tR)}(図2,図3参照)の時系列データ群11を取得するデータ取得部12と、時系列データ群11の時間軸(図2,図3参照)において時間幅が調節された複数の時間区間Zg(Z1,Z2,…ZG)の少なくとも一つを、時間区間Zgの積算量の差分値(W(tm+n)−W(tm))(図4参照)と第1参照値ΔWに基づいて設定する時間区間設定部15と、設定された時間区間Zg(Z1,Z2,…ZG)に基づいて時系列データ群11を複数のプルームPg{W(tm),… W(tm+n-1)}(g=1〜G)(図4参照)に分割するプルーム分割部16と、を少なくとも備えている。 As shown in FIG. 1, the analysis condition setting device 10 according to the first embodiment includes a registration unit 14 that registers a first reference value ΔW that is set in advance, and an analysis value or an integrated value W (t r ) {W (t 0 ), W (t 1 ),... W (t R )} (see FIGS. 2 and 3), a data acquisition unit 12 that acquires the time series data group 11, and the time axis of the time series data group 11 ( 2, at least one, integrated quantity of the difference value of the time interval Z g in Figure 3 a plurality of time width in the reference) is adjusted time interval Z g (Z 1, Z 2 , ... Z G) (W ( t m + n ) −W (t m )) (see FIG. 4) and the time interval setting unit 15 set based on the first reference value ΔW, and the set time intervals Z g (Z 1 , Z 2 ,... Based on Z G ), the time series data group 11 is divided into a plurality of plumes P g {W (t m ),... W (t m + n−1 )} (g = 1 to G) (see FIG. 4). The plume division 16 I have at least.

登録部14に登録される第1参照値ΔWは、その値自身が予め設定されている場合以外に、演算部13において演算された値が設定される場合がある。この演算部13は、時間範囲がt0からtRである時系列データ群11における積算量の最大値W(t)(図3参照)を、指定した分割数で除算して第1参照値ΔWを演算する。これにより、積算量の最大値W(t)の大きさにかかわらず、時系列データ群11は、分割数に対応する数のプルームPg{W(tm),… W(tm+n-1)}(g=1〜G)(図4参照)に適切に分割される。 The first reference value ΔW registered in the registration unit 14 may be set to a value calculated by the calculation unit 13 other than the case where the value itself is set in advance. The calculation unit 13 divides the maximum value W (t R ) (see FIG. 3) of the integrated amount in the time series data group 11 in which the time range is t 0 to t R by the specified number of divisions and refers to the first. The value ΔW is calculated. Thus, regardless of the maximum value W (t R ) of the integrated amount, the time-series data group 11 has a number of plumes P g {W (t m ),... W (t m + ) corresponding to the number of divisions. n-1 )} (g = 1 to G) (see FIG. 4).

上述の通り、時間区間設定部15は、時系列データ群11の時間軸(図2,図3参照)において時間幅が調節された複数の時間区間Zg(Z1,Z2,…ZG)の少なくとも一つを、積算量の差分値(W(tm+n)−W(tm)など)(図4参照)と第1参照値ΔWに基づいて、すなわち時間区間Zgの積算量の差分値と第1参照値ΔWを参照して設定する。 As described above, the time interval setting unit 15 has a plurality of time intervals Z g (Z 1 , Z 2 ,... Z G ) whose time widths are adjusted on the time axis (see FIGS. 2 and 3) of the time series data group 11. ) Is calculated based on the difference value (W (t m + n ) −W (t m ) etc.) (see FIG. 4) and the first reference value ΔW, that is, the integration of the time interval Z g . It is set with reference to the amount difference value and the first reference value ΔW.

一般に解析値又はその積算量W(tr){W(t0),W(t1),…W(tR)}は時間的な離散データとして得られることが多いことから、後述するように、時間区分設定部15は、nを1ずつ増やしながら積算量の差分値を順次計算していき、差分値(W(tm+n)−W(tm))が第1参照値ΔWを超えた場合にW(tm)からW(tm+n-1)となる区間を時間区間ZgのプルームPg{W(tm),… W(tm+n-1)}(g=1〜G)として設定する。換言するとこの場合、時間区分設定部15は、少なくとも一つの時間区間Zgを、積分量の差分値(W(tm+n-1)−W(tm))が第1参照値ΔW以下の最大値となるように設定している。 In general, the analysis value or its integrated value W (t r ) {W (t 0 ), W (t 1 ),... W (t R )} is often obtained as temporal discrete data, and will be described later. In addition, the time segment setting unit 15 sequentially calculates the difference value of the integrated amount while increasing n by 1, and the difference value (W (t m + n ) −W (t m )) is the first reference value ΔW. If the interval exceeds W (t m ) and W (t m + n-1 ), the plume P g {W (t m ), ... W (t m + n-1 )} of the time interval Z g (G = 1 to G) is set. In other words, in this case, the time segment setting unit 15 has at least one time interval Z g , and the difference value (W (t m + n−1 ) −W (t m )) of the integration amount is equal to or less than the first reference value ΔW The maximum value is set.

さらに、解析値又はその積算量W(tr){W(t0),W(t1),…W(tR)}は時間的な離散データの量が十分に多い場合などには、第1参照値ΔWに予め定めた一定幅を持たせた構成としても構わない。 Furthermore, the analysis value or its integrated value W (t r ) {W (t 0 ), W (t 1 ),... W (t R )} is, for example, when the amount of temporal discrete data is sufficiently large. The first reference value ΔW may have a predetermined constant width.

図2に示した例では、時間的な離散データとして得られた解析値又はその積算量W(tr){W(t0),W(t1),…W(tR)}を内挿または外挿して連続データとしての近時曲線(関数)を定め、時間区分設定部15は時間区間Zgを、この連続データである近似曲線に基づいて、各々の時間区間Zgにおける積算量の差分値(W(tm+n)−W(tm))(図4参照)が予め定めた一定値である第1参照値ΔWとなるように設定している。 In the example shown in FIG. 2, the analysis value obtained as temporally discrete data or its integrated value W (t r ) {W (t 0 ), W (t 1 ),... W (t R )} is included. By inserting or extrapolating, a recent curve (function) as continuous data is determined, and the time segment setting unit 15 determines the time interval Z g based on the approximate curve that is the continuous data and the integrated amount in each time interval Z g . The difference value (W (t m + n ) −W (t m )) (see FIG. 4) is set to a first reference value ΔW that is a predetermined constant value.

レベル3PRAでは、放出点からの距離毎の早期被ばく線量、長期被ばく線量、被ばく経路(クラウドシャイン、グランドシャイン、吸入、摂取、再浮遊物質吸入)ごとの線量など、目的に応じて多数の項目の評価が可能である。また、気象条件の変動を考慮し、年間気象変化に対する平均値、50%値、95%値、最大値などの統計値を出力する。ソースターム条件と評価対象項目、気象統計値などによって、分割数の指定値と解析精度及びプルーム単体総数の関係が異なるため、ユーザーがソースタームの種類とレベル3PRAでの評価対象項目、気象統計値を選択すると、分割数と解析精度への影響及びプルーム単体数が出力されるデータベースを備え、ユーザーをサポートする。データベースはあらかじめ代表的ソースタームについて本発明の方法でプルーム分割を行い、レベル3PRAを実施した結果に基づき作成する。解析精度については、例えばプルーム単体の放出継続時間を気象データの単位時間である1時間で時間的に等分割したものをベースケースとし、分割数の値を変えて本発明の手法で分割した場合とベースケースとのレベル3PRA評価結果との誤差をデータベースとして保有しておく方法がある。これにより、ユーザーが解析の目的に応じ、解析精度と解析負荷のトレードオフの関係から適切なプルーム分割を行うための分割数を指定することができる。   In Level 3PRA, there are many items depending on the purpose, such as early exposure dose at each distance from the release point, long-term exposure dose, dose for each exposure route (cloud shine, ground shine, inhalation, ingestion, inhalation of resuspended substances). Evaluation is possible. In addition, taking into account fluctuations in weather conditions, statistical values such as an average value, 50% value, 95% value, and maximum value with respect to annual weather changes are output. Since the relationship between the specified value of the number of divisions, the analysis accuracy, and the total number of plumes differs depending on the source term conditions, evaluation target items, weather statistics, etc. When the is selected, a database that outputs the number of divisions, the effect on analysis accuracy, and the number of single plumes is provided to support the user. The database is created based on the result of performing the plume division on the representative source terms by the method of the present invention and performing the level 3 PRA in advance. For analysis accuracy, for example, when the duration of a single plume is equally divided in time by 1 hour which is the unit time of weather data, the base case is used, and the value of the number of divisions is changed and divided by the method of the present invention. There is a method of storing an error between the level 3PRA evaluation result of the base case and the base case as a database. Thereby, the user can designate the number of divisions for performing appropriate plume division from the relationship between the analysis accuracy and the analysis load according to the purpose of analysis.

データ取得部12で取得される時系列データ群11は、PRAのレベル2から引き渡されるソースタームである。このソースタームは、一般に放射性物質を物理的・化学的性質等の条件によっていくつかのグループに分類されており、それぞれのグループ毎に環境への放出量あるいは初期内蔵量に対する放出割合を時系列データ群11として出力したものである。第1実施形態では、放射性物質のグループが一つの場合を取り上げている。   The time series data group 11 acquired by the data acquisition unit 12 is a source term delivered from the PRA level 2. This source term is generally classified into several groups of radioactive substances according to physical and chemical properties, etc., and the release rate to the environment or the initial built-in amount for each group is time-series data. This is output as the group 11. In the first embodiment, the case where there is one group of radioactive substances is taken up.

この時系列データ群11は、単位時間当たりの変化量を示す解析値であったり、この解析値が予め時間積分された積算量W(tr){W(t0),W(t1),…W(tR)であったりする。時系列データ群11がそのような解析値である場合は、データ取得部12の後段において、適時、積分処理が行われる。 This time-series data group 11 is an analysis value indicating the amount of change per unit time, or an integrated amount W (t r ) {W (t 0 ), W (t 1 ) obtained by integrating the analysis value with time. , ... W (t R ). When the time series data group 11 has such an analysis value, an integration process is performed in a timely manner after the data acquisition unit 12.

時間区間設定部15で設定した複数の時間区間Zg(Z1,Z2,…ZG)の各々に基づいて、時系列データ群11が分割される。この時系列データ群11を分割して得られた複数のプルームPg{W(tm),… W(tm+n-1)}(g=1〜G)は、蓄積部17に蓄積される。なお表示部18において、複数のプルームPgに分割された時系列データ群11のグラフを表示させることができる。 The time series data group 11 is divided based on each of the plurality of time intervals Z g (Z 1 , Z 2 ,... Z G ) set by the time interval setting unit 15. A plurality of plumes P g {W (t m ),... W (t m + n-1 )} (g = 1 to G) obtained by dividing the time series data group 11 are accumulated in the accumulation unit 17. Is done. Note the display section 18 can display a graph of time series data group 11 is divided into a plurality of plumes P g.

さらに解析装置20は、上述した解析条件設定装置10と、複数のプルームPg(g=1〜G)を単位として所定の解析コードで解析を実行する解析部21と、この解析の結果を出力する出力部23と、を少なくとも備えている。ここで、解析部21とは、例えば1時間おきに情報提供される気象情報に基づいて、放出積算量の変化が時系列データ群11として与えられている放射性物質の拡散シミュレーションを行うものである。 Further, the analysis device 20 outputs the analysis result setting device 10, the analysis unit 21 that performs analysis with a predetermined analysis code in units of a plurality of plumes P g (g = 1 to G), and the result of the analysis. And at least an output unit 23. Here, the analysis unit 21 performs, for example, a diffusion simulation of a radioactive substance in which a change in the accumulated release amount is given as the time-series data group 11 based on weather information provided every other hour. .

図4は第1実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャートである。この図4に基づいて図2及び図3におけるデータ処理を説明する。
まず、予め設定された第1参照値ΔWが登録される(S11)。次に積算量W(tr){W(t0),W(t1),…W(tR)}の時系列データ群11を取得する(S12)。なお取得される時系列データ群11は、積算量W(tr){W(t0),W(t1),…W(tR)}を例示しているが、この積算量の微分値に相当する解析値である場合もある。
FIG. 4 is a flowchart for explaining the analysis condition setting method or the analysis condition setting program according to the first embodiment. The data processing in FIGS. 2 and 3 will be described with reference to FIG.
First, a preset first reference value ΔW is registered (S11). Next, the time series data group 11 of the integrated amount W (t r ) {W (t 0 ), W (t 1 ),... W (t R )} is acquired (S12). The acquired time series data group 11 exemplifies the integrated amount W (t r ) {W (t 0 ), W (t 1 ),... W (t R )}. It may be an analysis value corresponding to the value.

次に図2を参照しつつ、m=0(S13)、n=1(S14)として、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 No)。引き続き、m=0、n=2(S22)として、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 Yes)。その結果、[t0,t1]の時間区間Z1が設定され、この時間区間Z1に割り付けられている{W(t0),W(t1)}がプルームP1として分割される(S23)。 Next, referring to FIG. 2, assuming that m = 0 (S13) and n = 1 (S14), the difference between the first reference value ΔW and the integrated amount (W (t m + n ) −W (t m )) (S18 No). Subsequently, as m = 0 and n = 2 (S22), the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount (S18 Yes). As a result, the set time interval Z 1 of [t 0, t 1], is allocated to this time interval Z 1 {W (t 0) , W (t 1)} is split as plumes P 1 (S23).

引き続いて、m=2(S24 No,S25)、n=1(S14)として、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 No)。さらに、m=2、n=2(S22)として、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18Yes)。その結果、[t2,t3]の時間区間Z2が設定され、この時間区間Z2に割り付けられている{W(t2),W(t3)}がプルームP2として分割される(S23)。 Subsequently, as m = 2 (S24 No, S25) and n = 1 (S14), the first reference value ΔW and the difference value (W (t m + n ) −W (t m )) of the integrated amount are compared. (S18 No). Further, as m = 2 and n = 2 (S22), the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount (Yes in S18). As a result, the set time interval Z 2 of [t 2, t 3], are assigned to the time interval Z 2 {W (t 2) , W (t 3)} is split as plumes P 2 (S23).

引き続いて、m=4(S24 No,S25)、n=1,2として、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 No)。さらに、m=4、n=3(S22)として、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18Yes)。その結果、[t4,t6]の時間区間Z3が設定され、この時間区間Z3に割り付けられている{W(t4),W(t5),W(t6)}がプルームP3として分割される(S23)。 Subsequently, assuming that m = 4 (S24 No, S25) and n = 1, 2, the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount. It performs (S18 No). Further, as m = 4 and n = 3 (S22), the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount (Yes in S18). As a result, a time interval Z 3 of [t 4 , t 6 ] is set, and {W (t 4 ), W (t 5 ), W (t 6 )} assigned to this time interval Z 3 is a plume. It is split as P 3 (S23).

引き続いて、m=7(S24 No,S25)、n=1〜12として、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 No)。さらに、m=7、n=13(S22)として、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18Yes)。その結果、[t7,t19]の時間区間Z4が設定され、この時間区間Z4に割り付けられている{W(t7),…,W(t19)}がプルームP4として分割される(S23)。 Subsequently, assuming that m = 7 (S24 No, S25) and n = 1 to 12, the comparison between the first reference value ΔW and the difference value (W (t m + n ) −W (t m )) of the integrated amount is performed. It performs (S18 No). Further, as m = 7 and n = 13 (S22), the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount (S18 Yes). As a result, [t 7, t 19] time interval Z 4 of the set, are assigned to the time interval Z 4 {W (t 7) , ..., W (t 19)} is split as plumes P 4 (S23).

上述したフローが、最後の積分量W(tR)まで繰り返され、時系列データ群11は、複数のプルームPg(g=1〜G)に分割されることになる。これにより、プルーム分割が粗く設定されることなく、解析精度と解析負荷の適切な関係を維持することができる。 The above-described flow is repeated until the last integration amount W (t R ), and the time-series data group 11 is divided into a plurality of plumes P g (g = 1 to G). Thereby, an appropriate relationship between the analysis accuracy and the analysis load can be maintained without setting the plume division coarsely.

(第2実施形態)
次に図1、図5、図6、図7を参照して本発明における第2実施形態について説明する。なお、図1に示される第2実施形態の解析条件設定装置10において、第1実施形態と重複する構成については説明を省略する。
(Second Embodiment)
Next, a second embodiment of the present invention will be described with reference to FIG. 1, FIG. 5, FIG. 6, and FIG. In the analysis condition setting device 10 of the second embodiment shown in FIG. 1, the description of the same configuration as that of the first embodiment is omitted.

第2実施形態の解析条件設定装置10の登録部14には、第1参照値ΔWに加え、予め設定された第2参照値Δtmin及び第3参照値Δtmaxがさらに登録されている。そして、時間区間Zg(Z1,Z2,…ZG)は、第2参照値Δtminよりも小さくならないように分割され、さらに第3参照値Δtmaxよりも大きくならないように分割される。ここで、第2参照値Δtminの設定は、解析部21で利用する気象情報等のデータの更新間隔(例えば一時間程度)が一つの目安となる。 In addition to the first reference value ΔW, a preset second reference value Δt min and third reference value Δt max are further registered in the registration unit 14 of the analysis condition setting device 10 of the second embodiment. Then, the time interval Z g (Z 1 , Z 2 ,... Z G ) is divided so as not to be smaller than the second reference value Δt min , and further divided so as not to be larger than the third reference value Δt max. . Here, the setting of the second reference value Δt min is based on an update interval (for example, about one hour) of data such as weather information used by the analysis unit 21.

図5は第2実施形態に係る解析条件設定装置で処理される時系列データ群の立ち上がり部分の詳細を示すグラフである。図6は第2実施形態に係る解析条件設定装置で処理される時系列データ群の全体を示すグラフである。図7は第2実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャートである。この図7に基づいて図5及び図6におけるデータ処理を説明する。   FIG. 5 is a graph showing details of a rising portion of a time-series data group processed by the analysis condition setting device according to the second embodiment. FIG. 6 is a graph showing the entire time-series data group processed by the analysis condition setting device according to the second embodiment. FIG. 7 is a flowchart for explaining an analysis condition setting method or an analysis condition setting program according to the second embodiment. The data processing in FIGS. 5 and 6 will be described based on FIG.

まず、予め設定された第1参照値ΔW、第2参照値Δtmin及び第3参照値Δtmaxが登録される(S11)。次に積算量W(tr){W(t0),W(t1),…W(tR)}の時系列データ群11を取得する(S12)。なお取得される時系列データ群11は、積算量W(tr){W(t0),W(t1),…W(tR)}を例示しているが、この積算量の微分値に相当する解析値である場合もある。 First, a preset first reference value ΔW, second reference value Δt min and third reference value Δt max are registered (S11). Next, the time series data group 11 of the integrated amount W (t r ) {W (t 0 ), W (t 1 ),... W (t R )} is acquired (S12). The acquired time series data group 11 exemplifies the integrated amount W (t r ) {W (t 0 ), W (t 1 ),... W (t R )}. It may be an analysis value corresponding to the value.

次に図5を参照しつつ、m=0(S13)、n=1,2(S14,S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 No)。そして、m=0(S13)、n=3(S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 Yes)。 Next, referring to FIG. 5, assuming that m = 0 (S13) and n = 1, 2 (S14, S16), the comparison between the second reference value Δt min and the time width (t m + n −t m ) is performed. It performs (S15 No). Then, as m = 0 (S13) and n = 3 (S16), the second reference value Δt min is compared with the time width (t m + n −t m ) (S15 Yes).

m=0,n=3のまま引き続いて、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 Yes)。その結果、[t0,t2]の時間区間Z1が設定され、この時間区間Z1に割り付けられている{W(t0),W(t1),W(t2)}がプルームP1として分割される(S23)。 Subsequently, the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount while maintaining m = 0 and n = 3 (S18 Yes). As a result, a time interval Z 1 of [t 0 , t 2 ] is set, and {W (t 0 ), W (t 1 ), W (t 2 )} assigned to this time interval Z 1 is a plume. Divided as P 1 (S23).

引き続いて、m=3(S24 No,S25)、n=1,2(S14,S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 No)。そして、m=3のままn=3(S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 Yes)。 Subsequently, as m = 3 (S24 No, S25) and n = 1, 2 (S14, S16), the second reference value Δt min is compared with the time width (t m + n −t m ) (S15). No). Then, with m = 3, n = 3 (S16), and the second reference value Δt min is compared with the time width (t m + n −t m ) (S15 Yes).

m=3,n=3のまま引き続いて、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 Yes)。その結果、[t3,t5]の時間区間Z2が設定され、この時間区間Z2に割り付けられている{W(t3),W(t4),W(t5)}がプルームP2として分割される(S23)。 Subsequently, the first reference value ΔW is compared with the integrated value difference value (W (t m + n ) −W (t m )) while maintaining m = 3 and n = 3 (Yes in S18). As a result, a time interval Z 2 of [t 3 , t 5 ] is set, and {W (t 3 ), W (t 4 ), W (t 5 )} assigned to this time interval Z 2 is a plume. Divided as P 2 (S23).

引き続いて、m=6(S24 No,S25)、n=1(S14)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 No)。そして、m=6のままn=2(S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 Yes)。 Subsequently, as m = 6 (S24 No, S25) and n = 1 (S14), the second reference value Δt min is compared with the time width (t m + n −t m ) (S15 No). Then, with m = 6, n = 2 (S16), and the second reference value Δt min is compared with the time width (t m + n −t m ) (S15 Yes).

m=6,n=2のまま引き続いて、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 No)。さらに、m=6のままn=3に更新し(S21 No,S22)、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18Yes)。その結果、[t6,t8]の時間区間Z3が設定され、この時間区間Z3に割り付けられている{W(t6),W(t7),W(t8)}がプルームP3として分割される(S23)。 Subsequently, the first reference value ΔW and the difference value (W (t m + n ) −W (t m )) of the integrated amount are compared with m = 6 and n = 2 (No in S18). Further, while updating m = 6, it is updated to n = 3 (S21 No, S22), and the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount. Perform (S18 Yes). As a result, a time interval Z 3 of [t 6 , t 8 ] is set, and {W (t 6 ), W (t 7 ), W (t 8 )} assigned to this time interval Z 3 is a plume. It is split as P 3 (S23).

引き続いて、m=9(S24 No,S25)、n=1,2(S14,S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 No)。そして、m=9のままn=3(S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 Yes)。 Subsequently, as m = 9 (S24 No, S25) and n = 1, 2 (S14, S16), the second reference value Δt min is compared with the time width (t m + n −t m ) (S15). No). Then, with m = 9, n = 3 (S16), and the second reference value Δt min is compared with the time width (t m + n −t m ) (S15 Yes).

m=9のまま引き継いで、n=4〜6(S22)とし、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 No)。さらに、m=9のままn=7に更新し(S21 No,S22)、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18Yes)。その結果、[t9,t15]の時間区間Z4が設定され、この時間区間Z4に割り付けられている{W(t9),…,W(t15)}がプルームP4として分割される(S23)。 Taking m = 9 into succession, n = 4 to 6 (S22), and the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount ( S18 No). Furthermore, m = 9 is updated to n = 7 (S21 No, S22), and the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount. Perform (S18 Yes). As a result, [t 9, t 15] time interval Z 4 of the set, are assigned to the time interval Z 4 {W (t 9) , ..., W (t 15)} divided as plumes P 4 (S23).

引き続いて、m=16(S24 No,S25)、n=1,2(S14,S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 No)。そして、m=16のままn=3(S16)として、第2参照値Δtminと時間幅(tm+n−tm)との対比を行う(S15 Yes)。 Subsequently, as m = 16 (S24 No, S25) and n = 1, 2 (S14, S16), the second reference value Δt min is compared with the time width (t m + n −t m ) (S15). No). Then, with m = 16, n = 3 (S16), and the second reference value Δt min is compared with the time width (t m + n −t m ) (S15 Yes).

m=16のまま引き継いで、n=4〜9(S22)とし、第1参照値ΔWと積算量の差分値(W(tm+n)−W(tm))との対比を行う(S18 No,S21 No)。さらに、m=16のままn=10に更新すると、第3参照値Δtmaxと時間幅(tm+n−tm)との対比が行われる(S21 Yes)。その結果、[t16,t25]の時間区間Z5が設定され、この時間区間Z5に割り付けられている{W(t16),…,W(t25)}がプルームP5として分割される(S23)。 Taking m = 16 into succession, n = 4 to 9 (S22), and the first reference value ΔW is compared with the difference value (W (t m + n ) −W (t m )) of the integrated amount ( S18 No, S21 No). Furthermore, when updating to n = 10 with m = 16, the third reference value Δt max is compared with the time width (t m + n −t m ) (S21 Yes). As a result, [t 16, t 25] time interval Z 5 are set, are assigned to the time interval Z 5 {W (t 16) , ..., W (t 25)} is split as plumes P 5 (S23).

上述したフローが、最後の積分量W(tR)まで繰り返され、時系列データ群11は、複数のプルームPg(g=1〜G)に分割されることになる。これにより、プルーム分割が粗く設定されることなく、かつプルームの数が多すぎることもなく、解析精度と解析負荷の適切な関係を維持することができる。 The above-described flow is repeated until the last integration amount W (t R ), and the time-series data group 11 is divided into a plurality of plumes P g (g = 1 to G). As a result, it is possible to maintain an appropriate relationship between the analysis accuracy and the analysis load without setting the plume division coarsely and without having too many plumes.

(第3実施形態)
次に図1、図8及び図9を参照して本発明における第3実施形態について説明する。なお、図1に示される第3実施形態の解析条件設定装置10において、第1実施形態及び第2実施形態と重複する構成については説明を省略する。
(Third embodiment)
Next, a third embodiment of the present invention will be described with reference to FIG. 1, FIG. 8, and FIG. In the analysis condition setting device 10 of the third embodiment shown in FIG. 1, the description of the same configuration as that of the first embodiment and the second embodiment is omitted.

第3実施形態に係る解析条件設定装置10(図1)において、データ取得部12は、複数の時系列データ群W1(tr),W2(tr),…,WQ(tr)(図8参照)を取得する。そして、これら時系列データ群W1(tr),W2(tr),…,WQ(tr)の各々に設定された複数の時間区間Zg(Z1,Z2,…ZG)は、相互に対応するもの同士の時間幅が等しくなるように設定されている。さらにプルーム分割部16は、複数の時系列データ群W1(tr),W2(tr),…,WQ(tr)を、相互に対応する時間区間Zg(Z1,Z2,…ZG)ごとに、プルームPg(P1,P2,…PG)に分割する。 In the analysis condition setting device 10 (FIG. 1) according to the third embodiment, the data acquisition unit 12 includes a plurality of time series data groups W 1 (t r ), W 2 (t r ),..., W Q (t r ) (See FIG. 8). Then, these time-series data group W 1 (t r), W 2 (t r), ..., W Q (t r) , each set being a plurality of time intervals Z g of (Z 1, Z 2, ... Z G ) is set so that the time widths of the corresponding ones are equal. Further, the plume dividing unit 16 divides a plurality of time series data groups W 1 (t r ), W 2 (t r ),..., W Q (t r ) into mutually corresponding time intervals Z g (Z 1 , Z 2, ... for each Z G), the plume P g (P 1, P 2 , divides ... to P G).

図8は第3実施形態に係る解析条件設定装置で処理される複数の時系列データ群を示すグラフである。図9は第3実施形態に係る解析条件設定方法又は解析条件設定プログラムを説明するフローチャートである。この図9は、図8から追加された(S17)、(S19)、(S20)のフローに基づいて、異なる複数の時系列データ群W1(tr),W2(tr),…,WQ(tr)は、相互に対応する時間区間Zg(Z1,Z2,…ZG)が一致し、共通のプルームPg(P1,P2,…PG)に分割することができる。 FIG. 8 is a graph showing a plurality of time-series data groups processed by the analysis condition setting device according to the third embodiment. FIG. 9 is a flowchart illustrating an analysis condition setting method or an analysis condition setting program according to the third embodiment. 9 is based on the flow of (S17), (S19), (S20) added from FIG. 8, and a plurality of different time series data groups W 1 (t r ), W 2 (t r ),. , W Q (t r ) coincide with each other corresponding time intervals Z g (Z 1 , Z 2 ,... Z G ), and are divided into common plumes P g (P 1 , P 2 ,... P G ). can do.

複数の時系列データ群W1(tr),W2(tr),…,WQ(tr)を用いることにより、放射性物質がいくつかのグループに分類されている場合について模擬することができる。環境に放出された放射性物質による被ばく線量評価で重要な放射性物質のグループは、希ガスを代表するグループ、よう素を代表するグループ、セシウムを代表するグループである。 Simulate the case where radioactive substances are classified into several groups by using multiple time series data groups W 1 (t r ), W 2 (t r ), ..., W Q (t r ) Can do. The groups of radioactive materials that are important for radiation dose assessment by radioactive materials released to the environment are the group that represents rare gas, the group that represents iodine, and the group that represents cesium.

希ガスは放射性雲からの外部被ばくへの寄与が大きく、フィルターで除去されにくいことから、フィルターベント実施時において被ばく線量に支配的となる。よう素は、吸入した場合の内部被ばく線量に対する寄与が大きい。セシウムは地表沈着による汚染地面からの長期外部被ばく線量に対し支配的な核種である。これら3つのグループに着目してプルーム分割を行うことで、主な被ばく経路の線量を適切に評価する分割が可能となる。   Since noble gases contribute greatly to external exposure from radioactive clouds and are not easily removed by filters, they are dominant in the exposure dose during filter venting. Iodine contributes greatly to the internal exposure dose when inhaled. Cesium is the dominant nuclide for long-term external exposure doses from contaminated ground due to surface deposition. By performing plume division focusing on these three groups, it is possible to perform division that appropriately evaluates the dose of main exposure routes.

(第4実施形態)
図10は第4実施形態に係る解析条件設定装置で処理される複数の時系列データ群を示すグラフである。第4実施形態に係る解析装置20は、第1位置に由来する時系列データ群W(tr)に基づいて解析した第1解析結果と第2位置に由来する時系列データ群V(tr)に基づいて解析した第2解析結果とを蓄積する蓄積部22を備えている。さらに、これら第1解析結果及び第2解析結果を加算して出力させる加算部24を、さらに備えている。
(Fourth embodiment)
FIG. 10 is a graph showing a plurality of time-series data groups processed by the analysis condition setting device according to the fourth embodiment. The analysis device 20 according to the fourth embodiment includes a first analysis result analyzed based on the time series data group W (t r ) derived from the first position and a time series data group V (t r derived from the second position. ) Based on the second analysis result is stored. Further, an adder 24 that adds and outputs the first analysis result and the second analysis result is further provided.

第1位置に由来する時系列データ群W(tr)及び第2位置に由来する時系列データ群V(tr)は、開始時刻も異なるために、それぞれの体系で設定された時間区間Z1,Z2,…Z4及び時間区間X1,X2,X3に基づいてプルーム分割されている。しかし、それぞれの時系列データ群W(tr),V(tr)から導かれた第1解析結果及び第2解析結果は、時間軸をあわせて互いを加算して評価することができる。 Since the time series data group W (t r ) derived from the first position and the time series data group V (t r ) derived from the second position have different start times, the time interval Z set in each system is set. Plume is divided based on 1 , Z 2 ,... Z 4 and time intervals X 1 , X 2 , X 3 . However, the first analysis result and the second analysis result derived from the respective time series data groups W (t r ) and V (t r ) can be evaluated by adding each other along the time axis.

事故時の放射性物質放出の代表的なシナリオでは、事故進展に伴い、初めに建屋からの放射性物質の漏洩があり、その後、ベントスタックからフィルターベントが行われるケースがある。この場合には、放射性物質の放出位置が建屋とベントスタックの2か所となり、それぞれ、放射性物質の放出高さ、初期拡散幅が異なるため、レベル3PRAでは、これらを別々のプルームとして扱う。また、第1位置及び第2位置のそれぞれにおいてグループが一つの場合を取り上げているが、グループが複数の場合も適用することができる。   In a typical scenario of radioactive material release at the time of an accident, there is a case where the radioactive material leaks from the building first and then the filter vent is performed from the vent stack as the accident progresses. In this case, the radioactive material is released at two locations, the building and the vent stack, and the radioactive material emission height and the initial diffusion width are different from each other. In Level 3PRA, these are handled as separate plumes. Moreover, although the case where there is one group in each of the first position and the second position is taken up, the case where there are a plurality of groups can also be applied.

また、第4実施形態において、演算部13により第1参照値ΔWを演算する場合は、第1位置及び第2位置のそれぞれにおいて、異なる分割数を指定することができる。事故時の放射性物質放出の代表的なシナリオでは、事故進展に伴い、初めに建屋からの放射性物質の漏洩があり、その後、ベントスタックからフィルターベントが行われるケースがある。放出形態の違いにより、適切な分割数の値は異なるため、放出位置毎にこの分割数を指定することで、より最適なプルーム分割が可能となる。   Further, in the fourth embodiment, when the calculation unit 13 calculates the first reference value ΔW, it is possible to specify different division numbers in each of the first position and the second position. In a typical scenario of radioactive material release at the time of an accident, there is a case where the radioactive material leaks from the building first and then the filter vent is performed from the vent stack as the accident progresses. Since the value of the appropriate number of divisions differs depending on the emission form, it is possible to perform more optimal plume division by designating this number of divisions for each emission position.

以上述べた少なくともひとつの実施形態の解析条件設定装置によれば、積算量の差分値が参照値となるように時系列データ群をプルーム分割することにより、解析精度と解析負荷のトレードオフの関係を考慮しつつ、解析負荷を低減しかつ現実的な計算時間で評価することが可能となる。   According to the analysis condition setting device of at least one embodiment described above, the relationship between the analysis accuracy and the analysis load is traded off by plume dividing the time series data group so that the difference value of the integrated amount becomes the reference value. Thus, it is possible to reduce the analysis load and evaluate in a realistic calculation time.

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。また、解析条件設定装置の構成要素は、コンピュータのプロセッサで実現することも可能であり、解析条件設定プログラムにより動作させることが可能である。   Although several embodiments of the present invention have been described, these embodiments are presented by way of example and are not intended to limit the scope of the invention. These embodiments can be implemented in various other forms, and various omissions, replacements, changes, and combinations can be made without departing from the scope of the invention. These embodiments and their modifications are included in the scope and gist of the invention, and are also included in the invention described in the claims and the equivalents thereof. The components of the analysis condition setting device can be realized by a processor of a computer and can be operated by an analysis condition setting program.

以上説明した解析条件設定装置は、専用のチップ、FPGA(Field Programmable Gate Array)、GPU(Graphics Processing Unit)、又はCPU(Central Processing Unit)などのプロセッサを高集積化させた制御装置と、ROM(Read Only Memory)やRAM(Random Access Memory)などの記憶装置と、HDD(Hard Disk Drive)やSSD(Solid State Drive)などの外部記憶装置と、ディスプレイなどの表示装置と、マウスやキーボードなどの入力装置と、通信I/Fとを、備えており、通常のコンピュータを利用したハードウェア構成で実現できる。   The analysis condition setting device described above includes a control device in which a processor such as a dedicated chip, FPGA (Field Programmable Gate Array), GPU (Graphics Processing Unit), or CPU (Central Processing Unit) is highly integrated, and ROM ( Storage devices such as Read Only Memory (RAM) and Random Access Memory (RAM), external storage devices such as HDD (Hard Disk Drive) and SSD (Solid State Drive), display devices such as a display, and inputs such as a mouse and a keyboard The apparatus and the communication I / F are provided, and can be realized by a hardware configuration using a normal computer.

また解析条件設定装置で実行されるプログラムは、ROM等に予め組み込んで提供される。もしくは、このプログラムは、インストール可能な形式又は実行可能な形式のファイルでCD−ROM、CD−R、メモリカード、DVD、フレキシブルディスク(FD)等のコンピュータで読み取り可能な記憶媒体に記憶されて提供するようにしてもよい。   A program executed by the analysis condition setting device is provided by being incorporated in advance in a ROM or the like. Alternatively, this program is stored in a computer-readable storage medium such as a CD-ROM, CD-R, memory card, DVD, or flexible disk (FD) as an installable or executable file. You may make it do.

また、本実施形態に係る解析条件設定装置で実行されるプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせて提供するようにしてもよい。また、装置10は、構成要素の各機能を独立して発揮する別々のモジュールを、ネットワーク又は専用線で相互に接続し、組み合わせて構成することもできる。   The program executed by the analysis condition setting apparatus according to the present embodiment may be stored on a computer connected to a network such as the Internet and provided by being downloaded via the network. In addition, the apparatus 10 can be configured by combining separate modules that perform each function of the component elements independently by a network or a dedicated line.

10…解析条件設定装置、11…時系列データ群、12…データ取得部、13…演算部、14…登録部、15…時間区間設定部、16…プルーム分割部、17…蓄積部、18…表示部、20…解析装置、21…解析部、22…蓄積部、23…出力部、24…加算部、W(tr){W(t0),W(t1),…W(tR)}…積算量、Zg(Z1,Z2,…ZG)…時間区間、W(tm+n)−W(tm)…差分値、ΔW…第1参照値、Pg{W(tm),… W(tm+n-1)}…プルーム。 DESCRIPTION OF SYMBOLS 10 ... Analysis condition setting apparatus, 11 ... Time series data group, 12 ... Data acquisition part, 13 ... Operation part, 14 ... Registration part, 15 ... Time section setting part, 16 ... Plume division part, 17 ... Accumulation part, 18 ... Display unit, 20 ... analyzing device, 21 ... analyzing unit, 22 ... accumulating unit, 23 ... output unit, 24 ... adding unit, W (t r ) {W (t 0 ), W (t 1 ), ... W (t R )} ... integrated amount, Z g (Z 1 , Z 2 , ... Z G ) ... time interval, W (t m + n ) -W (t m ) ... difference value, ΔW ... first reference value, P g {W (t m ), ... W (t m + n-1 )} ... Plum.

Claims (7)

予め設定された第1参照値を登録する登録部と、
解析値又はその積算量の時系列データ群を取得するデータ取得部と、
前記時系列データ群の時間軸において時間幅を有する時間区間を、前記積算量の差分値と前記第1参照値に基づいて設定する時間区間設定部と、
設定された前記時間区間に基づいて、前記時系列データ群をプルームに分割するプルーム分割部と、を備えることを特徴とする解析条件設定装置。
A registration unit for registering a preset first reference value;
A data acquisition unit for acquiring a time-series data group of analysis values or integrated values thereof;
A time interval setting unit that sets a time interval having a time width on the time axis of the time-series data group based on the difference value of the integrated amount and the first reference value;
An analysis condition setting device comprising: a plume dividing unit that divides the time-series data group into plumes based on the set time interval.
請求項1に記載の解析条件設定装置において、
前記登録部に、予め設定された第2参照値がさらに登録されており、
前記時間区間は、前記第2参照値よりも小さくならないように分割されることを特徴とする解析条件設定装置。
In the analysis condition setting device according to claim 1,
A second reference value set in advance is further registered in the registration unit,
The analysis condition setting device, wherein the time interval is divided so as not to be smaller than the second reference value.
請求項1又は請求項2に記載の解析条件設定装置において、
前記登録部に、予め設定された第3参照値がさらに登録されており、
前記時間区間は、前記第3参照値よりも大きくならないように分割されることを特徴とする解析条件設定装置。
In the analysis condition setting device according to claim 1 or claim 2,
A third reference value set in advance is further registered in the registration unit,
The analysis condition setting device, wherein the time interval is divided so as not to be larger than the third reference value.
請求項1から請求項3のいずれか1項に記載の解析条件設定装置において、
前記データ取得部は、複数の前記時系列データ群を取得し、
前記時系列データ群の各々に設定された複数の時間区間は、相互に対応するもの同士の前記時間幅が等しくなるように設定されており、
プルーム分割部は、複数の前記時系列データ群を、相互に対応する時間区間ごとに、前記プルームに分割する、ことを特徴とする解析条件設定装置。
In the analysis condition setting device according to any one of claims 1 to 3,
The data acquisition unit acquires a plurality of the time series data groups,
The plurality of time intervals set in each of the time series data groups are set so that the time widths of the corresponding ones are equal to each other,
The plume dividing unit divides the plurality of time-series data groups into the plumes for each time interval corresponding to each other.
請求項1から請求項4のいずれか1項に記載の解析条件設定装置と、
複数の前記プルームを単位として所定の解析コードで解析を実行する解析部と、
前記解析の結果を出力する出力部と、を備えることを特徴とする解析装置。
The analysis condition setting device according to any one of claims 1 to 4,
An analysis unit that performs analysis with a predetermined analysis code in units of a plurality of the plumes;
And an output unit for outputting the result of the analysis.
予め設定された第1参照値を登録する登録ステップと、
解析値又はその積算量の時系列データ群を取得するデータ取得ステップと、
前記時系列データ群の時間軸において時間幅を有する時間区間を、前記積算量の差分値と前記第1参照値に基づいて設定する時間区間設定ステップと、
設定された前記時間区間に基づいて、前記時系列データ群をプルームに分割する分割ステップと、を含むことを特徴とする解析条件設定方法。
A registration step for registering a preset first reference value;
A data acquisition step for acquiring a time-series data group of the analysis value or its integrated amount;
A time interval setting step for setting a time interval having a time width on the time axis of the time series data group based on the difference value of the integrated amount and the first reference value;
A dividing step of dividing the time-series data group into plumes based on the set time interval.
コンピュータに、
予め設定された第1参照値を登録する登録ステップ、
解析値又はその積算量の時系列データ群を取得するデータ取得ステップ、
前記時系列データ群の時間軸において時間幅を有する時間区間を、前記積算量の差分値と前記第1参照値に基づいて設定する時間区間設定ステップ、
設定された前記時間区間に基づいて、前記時系列データ群を複数のプルームに分割する分割ステップ、を実行させることを特徴とする解析条件設定プログラム。
On the computer,
A registration step for registering a preset first reference value;
A data acquisition step for acquiring a time-series data group of the analysis value or the accumulated amount thereof,
A time interval setting step for setting a time interval having a time width on the time axis of the time series data group based on the difference value of the integrated amount and the first reference value;
An analysis condition setting program for executing a dividing step of dividing the time series data group into a plurality of plumes based on the set time interval.
JP2018058684A 2018-03-26 2018-03-26 Analysis condition setting device, method, and program, and analysis device Active JP6884724B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018058684A JP6884724B2 (en) 2018-03-26 2018-03-26 Analysis condition setting device, method, and program, and analysis device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018058684A JP6884724B2 (en) 2018-03-26 2018-03-26 Analysis condition setting device, method, and program, and analysis device

Publications (2)

Publication Number Publication Date
JP2019169104A true JP2019169104A (en) 2019-10-03
JP6884724B2 JP6884724B2 (en) 2021-06-09

Family

ID=68107435

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018058684A Active JP6884724B2 (en) 2018-03-26 2018-03-26 Analysis condition setting device, method, and program, and analysis device

Country Status (1)

Country Link
JP (1) JP6884724B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111931352A (en) * 2020-07-16 2020-11-13 嘉泰工程技术有限公司 Universal probability safety evaluation model conversion system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013170825A (en) * 2012-02-17 2013-09-02 Japan Atom Power Co Ltd:The Method for evaluating activated radioactivity and system for evaluating activated radioactivity
JP2014130117A (en) * 2012-12-28 2014-07-10 Mitsubishi Heavy Ind Ltd Diffusion state prediction device of diffuse matter and program
CN104933294A (en) * 2015-05-25 2015-09-23 苏州热工研究院有限公司 Method for constructing migration and propagation evaluation model of radioactive substance

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013170825A (en) * 2012-02-17 2013-09-02 Japan Atom Power Co Ltd:The Method for evaluating activated radioactivity and system for evaluating activated radioactivity
JP2014130117A (en) * 2012-12-28 2014-07-10 Mitsubishi Heavy Ind Ltd Diffusion state prediction device of diffuse matter and program
CN104933294A (en) * 2015-05-25 2015-09-23 苏州热工研究院有限公司 Method for constructing migration and propagation evaluation model of radioactive substance

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111931352A (en) * 2020-07-16 2020-11-13 嘉泰工程技术有限公司 Universal probability safety evaluation model conversion system

Also Published As

Publication number Publication date
JP6884724B2 (en) 2021-06-09

Similar Documents

Publication Publication Date Title
Zapartas et al. Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction
US11190535B2 (en) Methods and systems for inferring behavior and vulnerabilities from process models
Sanz-Forcada et al. A scenario of planet erosion by coronal radiation
US9189543B2 (en) Predicting service request breaches
CN107463904A (en) A kind of method and device for determining periods of events value
Ritenberga et al. Statistical modelling of non-stationary processes of atmospheric pollution from natural sources: example of birch pollen
CN113159362A (en) Forest fire grade prediction method and device
JP5717604B2 (en) Diffusion situation prediction system
Benito et al. Forecasting plant range collapse in a mediterranean hotspot: When dispersal uncertainties matter
JP6884724B2 (en) Analysis condition setting device, method, and program, and analysis device
US20160104125A1 (en) Computing system and method for providing information relating to maintenance actions
Garcia et al. Detection of period-spacing patterns due to the gravity modes of rotating dwarfs in the TESS southern continuous viewing zone
CN114708007A (en) Intelligent decomposition method and system for store sales plan
CN109325061A (en) A kind of method for processing abnormal data and device based on probability distribution
CN117035374B (en) Force cooperative scheduling method, system and medium for coping with emergency
Aguiar et al. Effect of chronological addition of records to species distribution maps: The case of Tonatia saurophila maresi (C hiroptera, P hyllostomidae) in S outh A merica
Shukla et al. An efficient Dijkstra’s semaphore based k-out-of-n Petri net voting model: An application to shutdown actuation system reliability
CN113419706A (en) Rapid random number generation method and system and inspection method and system thereof
CN108804640B (en) Data grouping method, device, storage medium and equipment based on maximized IV
CN109783646B (en) Data processing method and device
CN112907285A (en) Option data processing method based on self-adaptive grid and backward iteration integral
Ingram et al. Using early stage project data to predict change-proneness
US20140136161A1 (en) Precise simulation of progeny derived from recombining parents
JP6566499B2 (en) Shipment quantity prediction device, shipment quantity prediction method, and program
CN108399549B (en) Business object delivery method and device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200310

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210216

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210226

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210413

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210512

R150 Certificate of patent or registration of utility model

Ref document number: 6884724

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150