JP2001208848A - Radiation mesuring system and method - Google Patents
Radiation mesuring system and methodInfo
- Publication number
- JP2001208848A JP2001208848A JP2000019152A JP2000019152A JP2001208848A JP 2001208848 A JP2001208848 A JP 2001208848A JP 2000019152 A JP2000019152 A JP 2000019152A JP 2000019152 A JP2000019152 A JP 2000019152A JP 2001208848 A JP2001208848 A JP 2001208848A
- Authority
- JP
- Japan
- Prior art keywords
- positions
- measured
- dose rate
- value
- radiation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Measurement Of Radiation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Instructional Devices (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は放射線測定システム
及び方法に関し、特に地理上における線量率分布の作成
に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a radiation measurement system and method, and more particularly, to a method of generating a geographic dose rate distribution.
【0002】[0002]
【従来の技術及びその課題】地理上の線量率の分布を得
るためには、所定間隔をもった多数の地点において線量
率を実測し、その多数の実測値を地図上に表現するのが
理想的である。しかし、固定設置型のモニタリングポス
トの設置数を増大させると、システムコストが飛躍的に
増大し、また設置場所の確保や設置条件などの制約があ
るために多くの困難が伴う。一方、放射線測定器を自動
車などの移動体に搭載し、走行を行いながらあるいは駐
車位置を移動させて、線量率の測定を逐次的に行うこと
も可能である。しかし、その場合には、自動車が走行で
きる箇所でしか線量率の測定を行うことができないとい
う問題がある。いずれにしても、有限かつ非均等的に設
定された実測位置から得られる有限個の実測値をもっ
て、信頼性の高い線量率の分布を描くことが必要とな
る。2. Description of the Related Art In order to obtain a geographical distribution of dose rates, it is ideal to measure the dose rates at a number of points at predetermined intervals and to represent the measured values on a map. It is a target. However, if the number of fixed installation type monitoring posts is increased, the system cost is dramatically increased, and there are many difficulties due to restrictions on securing an installation place and installation conditions. On the other hand, it is also possible to mount the radiation measuring instrument on a moving body such as an automobile and measure the dose rate sequentially while traveling or moving the parking position. However, in that case, there is a problem that the dose rate can be measured only in a place where the automobile can travel. In any case, it is necessary to draw a highly reliable dose rate distribution using a finite number of measured values obtained from finite and non-uniformly set measured positions.
【0003】なお、特開平8−334563号公報には
緊急時における放射線モニタリングシステムについて開
示されている。[0003] Japanese Patent Laid-Open No. 8-334563 discloses a radiation monitoring system in an emergency.
【0004】本発明は、上記従来の課題に鑑みなされた
ものであり、その目的は、有限の実測値をもって線量率
などの放射線情報の分布を作成することにある。The present invention has been made in view of the above-mentioned conventional problems, and has as its object to create a distribution of radiation information such as a dose rate using a finite measured value.
【0005】[0005]
【課題を解決するための手段】(1)上記目的を達成す
るために、本発明は、地理上の複数の実測位置において
放射線の線量率測定を行って実測値を求める測定手段
と、前記複数の実測位置についての実測値に基づいて、
地理上の複数の補間位置について線量率の推定値を推定
する推定手段と、前記複数の実測位置についての実測値
と前記複数の補間位置についての推定値とに基づいて、
地理上における線量率分布を形成する形成手段と、を含
むことを特徴とする。(1) In order to achieve the above object, the present invention provides a measuring means for measuring a radiation dose rate at a plurality of geographically measured positions to obtain a measured value; Based on the measured value of the measured position of
Estimating means for estimating the estimated value of the dose rate for a plurality of geographical interpolation positions, based on the measured values for the plurality of measured positions and the estimated values for the plurality of interpolation positions,
Forming means for forming a dose rate distribution in geography.
【0006】上記構成によれば、複数の実測値に基づい
て、複数の補間位置について推定値が推定され、それら
の実測値及び推定値に基づいて、地理上における線量率
分布を形成できる。よって、実際に測定を行っていない
位置についても線量率を合理的に推定できるという利点
がある。特に、各補間位置ごとに、複数の実測定位置か
ら距離を考慮し、それに基づいて線量率を推定するのが
望ましい。補間位置は固定的に設定してもよいし、ユー
ザー指定されるようにしてもよい。その補間位置の密度
や個数などをユーザー指定できるように構成するのが望
ましい。According to the above configuration, estimated values are estimated for a plurality of interpolation positions based on a plurality of actually measured values, and a geographical dose rate distribution can be formed based on the actually measured values and the estimated values. Therefore, there is an advantage that the dose rate can be reasonably estimated even at a position where measurement is not actually performed. In particular, for each interpolation position, it is desirable to consider the distance from a plurality of actual measurement positions and estimate the dose rate based on the distance. The interpolation position may be fixedly set or may be specified by the user. It is desirable that the density and the number of the interpolation positions be specified by the user.
【0007】望ましくは、前記測定手段は、前記複数の
実測位置ごとに固定設置された複数のモニタリング装置
を含む。あるいは、望ましくは、前記測定手段は、移動
体に搭載された放射線測定器及び移動***置検出器を含
む。Preferably, the measuring means includes a plurality of monitoring devices fixedly installed for each of the plurality of actually measured positions. Alternatively, preferably, the measuring means includes a radiation measuring instrument and a moving object position detector mounted on the moving object.
【0008】望ましくは、前記線量率分布は、等高線表
現及び色付け表現の少なくとも一方を利用して表現され
る。それらの表示形態をユーザー選択できるように構成
するのが望ましい。また、地理上の表示エリアを自在設
定できるようにしてもよい。更に、分布表示が合成され
る地図上に、地名、道路、河川などの地理情報を表すよ
うにするのが望ましい。加えて、原子力発電所、核燃料
施設などの原子力施設など明示するようにしてもよい。
望ましくは、前記複数の補間位置の密度を可変設定する
手段を含む。[0008] Preferably, the dose rate distribution is represented using at least one of a contour line expression and a coloring expression. It is desirable to configure such a display form that the user can select it. The geographic display area may be freely set. Furthermore, it is desirable that geographical information such as place names, roads, and rivers be displayed on a map on which the distribution display is synthesized. In addition, a nuclear facility such as a nuclear power plant or a nuclear fuel facility may be specified.
Preferably, the apparatus further includes means for variably setting the density of the plurality of interpolation positions.
【0009】(2)また、上記目的を達成するために、
本発明は、地理上の複数の実測位置において放射線の線
量率測定を行って実測値を求める測定手段と、空間統計
学的手法に従って、前記複数の実測位置についての実測
値に基づいて、地理上の複数の補間位置について線量率
の推定値を推定する推定手段と、前記複数の実測位置に
ついての実測値と前記複数の補間位置についての推定値
とに基づいて、地理上における線量率分布を形成する形
成手段と、を含むことを特徴とする。(2) To achieve the above object,
The present invention provides a measuring means for measuring a radiation dose rate at a plurality of actually measured positions on a geographical site to obtain an actually measured value, and a geographical method based on the actually measured values for the plurality of actually measured positions according to a spatial statistical technique. Estimating means for estimating an estimated value of the dose rate for a plurality of interpolated positions; and forming a geographical dose rate distribution based on the measured values for the plurality of actually measured positions and the estimated values for the plurality of interpolated positions. And forming means.
【0010】上記構成によれば、空間統計学的手法によ
って実測値に関する統計的な傾向から推定値が推定され
るので当該推定値の信頼性を高めることが可能となる。
この場合、複数の実測位置で囲まれる領域内に補間位置
が設定されるように実測位置を決定するのが望ましく、
またそれを考慮して、できる限り実測位置が分散するよ
うにするのが望ましい。According to the above configuration, since the estimated value is estimated from the statistical tendency of the actually measured value by the spatial statistical technique, the reliability of the estimated value can be improved.
In this case, it is desirable to determine the measured position so that the interpolation position is set in a region surrounded by a plurality of measured positions,
In consideration of this, it is desirable that the measured positions are dispersed as much as possible.
【0011】望ましくは、前記推定手段は、前記複数の
実測位置についての実測値と前記複数の実測位置の相互
間距離とに基づいて、各距離ごとの経験バリオグラムを
演算する手段と、前記各距離ごとの経験バリオグラムに
対し、理論バリオグラムを定義する理論バリオグラム関
数がフィッティングするように、その理論バリオグラム
関数に含まれるパラメータの値を決定する手段と、前記
決定されたパラメータの値によって理論共分散関数を決
定し、前記各補間位置ごとに、前記理論共分散関数に基
づいて前記複数の実測値が当該補間位置へ与える影響を
重み値として求める手段と、前記各補間位置ごとに、前
記各実測値に対してそれに対応する重み値を乗算した結
果を加算し、これにより前記推定値を演算する手段と、
を含む。ここで、望ましくは、前記理論バリオグラム関
数は、球形モデル、指数型モデル又はガウス型モデルの
1つを基礎とする関数である。Preferably, the estimating means calculates an empirical variogram for each of the distances based on the measured values of the plurality of measured positions and the distance between the plurality of measured positions. Means for determining a value of a parameter included in the theoretical variogram function so that a theoretical variogram function defining the theoretical variogram is fitted to each of the empirical variograms; Determining, for each of the interpolation positions, means for determining, as a weight value, the influence of the plurality of measured values on the interpolation position based on the theoretical covariance function, and for each of the interpolation positions, Means for calculating the estimated value by adding the result of multiplying the corresponding weight value, thereby calculating the estimated value;
including. Preferably, the theoretical variogram function is a function based on one of a spherical model, an exponential model, and a Gaussian model.
【0012】望ましくは、前記推定手段は、理論バリオ
グラムを定義する理論バリオグラム関数に含まれるパラ
メータの値を格納する手段と、前記パラメータの値によ
って理論共分散関数を決定し、前記各補間位置ごとに、
前記理論共分散関数に基づいて前記複数の実測値が当該
補間位置へ与える影響を重み値として求める手段と、前
記各補間位置ごとに、前記各実測値に対してそれに対応
する重み値を乗算した結果を加算し、これによって前記
推定値を演算する手段と、を含む。Preferably, the estimating means stores a value of a parameter included in a theoretical variogram function defining a theoretical variogram, and determines a theoretical covariance function based on the value of the parameter. ,
Means for obtaining, as a weight value, the influence of the plurality of measured values on the interpolation position based on the theoretical covariance function, and for each of the interpolation positions, each measured value is multiplied by a corresponding weight value. Means for summing the results and thereby calculating said estimate.
【0013】(3)また、上記目的を達成するために、
本発明は、地理上の複数の実測位置において放射線の線
量率測定を行って実測値を求める工程と、空間統計学的
手法に従って、前記複数の実測位置についての実測値に
基づいて、地理上の複数の補間位置について線量率の推
定値を推定する工程と、前記複数の実測位置についての
実測値と前記複数の補間位置についての推定値とに基づ
いて、地理上における線量率分布を形成する工程と、を
含むことを特徴とする。上記の推定工程及び形成工程を
実質的にソフトウエアにより実現するようにしてもよ
い。(3) To achieve the above object,
The present invention provides a step of obtaining a measured value by performing radiation dose rate measurement at a plurality of measured positions on a geographical basis, and according to a spatial statistical technique, based on the measured values for the plurality of measured positions, Estimating a dose rate estimate for a plurality of interpolated positions; and forming a geographical dose rate distribution based on the measured values for the plurality of measured positions and the estimated values for the plurality of interpolated positions. And characterized in that: The above estimation step and formation step may be substantially realized by software.
【0014】(4)推定方法の説明以下に、有限個の実
測値(実測線量率)から補間位置の推定値(推定線量
率)を算出する方法について、図9を用いて説明する。
この手法は、空間統計学の地理統計データ解析手法を応
用したものである。(4) Description of estimation method A method of calculating an estimated value (estimated dose rate) of an interpolation position from a finite number of actually measured values (actually measured dose rates) will be described below with reference to FIG.
This method applies the geostatistical data analysis method of spatial statistics.
【0015】実測値の収集を前提とし、まず、S201
では、後述の各距離ごとに「経験バリオグラム」が演算
される。ここで、経験バリオグラムγ(h)は下記のよ
うに定義される。Assuming that actual measurement values are collected, first, in S201
Then, an "experimental variogram" is calculated for each distance described below. Here, the empirical variogram γ (h) is defined as follows.
【0016】[0016]
【数1】 この経験バリオグラムγ(h)は、2点間の距離ごとの
統計的な傾向を表すものである。つまり、γ(h)は距
離hだけ離れた2点間の線量率の差の2乗平均の2分の
1、すなわち差分の2乗平均の2分の1に相当する。(Equation 1) This empirical variogram γ (h) represents a statistical tendency for each distance between two points. That is, γ (h) corresponds to one half of the root mean square of the difference between the dose rates between two points separated by the distance h, that is, one half of the root mean square of the difference.
【0017】例えば、実測定位置が3点の場合、各点の
線量率値をそれぞれZ1、Z2、Z3として、For example, when the actual measurement positions are three points, the dose rate values at each point are Z 1 , Z 2 , and Z 3 , respectively.
【数2】 となる。この経験バリオグラムは、距離の関数として、
一般的に図7に示す関数形を有する。ここで、ナゲット
(nugget)は、微視規模変動あるいは測定誤差に相当し、
バリオグラムγ(h)のh=0における値で定義され
る。シル(sill)は、hを無限大にした場合におけるγ
(h)での確率場の分散を表すものである。レンジ(ran
ge)は、自己相関がなくなる距離に相当する。(Equation 2) Becomes This empirical variogram, as a function of distance,
Generally, it has the function form shown in FIG. Where the nuggets
(nugget) is equivalent to microscopic scale fluctuation or measurement error,
It is defined by the value at h = 0 of the variogram γ (h). Sill is γ when h is infinite.
It represents the variance of the random field in (h). Range (ran
ge) is equivalent to the distance at which the autocorrelation disappears.
【0018】上記のS201において、実用的な規則と
して考慮すべき点が2つある。In the above S201, there are two points to be considered as a practical rule.
【0019】a)バリオグラムは、ペアの数が30以上
になる距離hのみで計算すべき b)バリオグラムの信頼できる距離はh < D/2 ここで、Dはデータの2点間の距離の内で最大になる距
離である。もちろん、上記条件は測定条件や要求精度な
どに応じて異なる。なお、実測位置は、理想的には地図
上を網羅しているのが望ましく、そうでない場合でもあ
る程度広く分散しているのが望ましい。A) The variogram should be calculated only at the distance h at which the number of pairs is 30 or more. B) The reliable distance of the variogram is h <D / 2, where D is the distance between two points in the data. Is the maximum distance. Of course, the above conditions differ depending on measurement conditions, required accuracy, and the like. Note that the measured positions ideally desirably cover the map, and even if not, it is desirable that the measured positions be widely dispersed to some extent.
【0020】S202では、上記のようにして求められ
た経験バリオグラムを利用して、「理論バリオグラム」
を規定する所定の関数についてのフィッティングが行わ
れる。その関数としては、各種のものがあるが、ここで
は、代表的な理論バリオグラムγ(h)として球型モデ
ルを例にとって説明する。なお、理論バリオグラムに
は、この球形モデルの他に、指数型モデル、ガウス型モ
デルなどがある。In S202, the "theoretical variogram" is obtained by using the empirical variogram obtained as described above.
Is performed for a predetermined function that defines Although there are various functions as the function, a spherical model will be described as a typical theoretical variogram γ (h). The theoretical variogram includes an exponential model, a Gaussian model, and the like in addition to the spherical model.
【0021】球形モデルの理論バリオグラムは以下のよ
うに定義される。The theoretical variogram of the spherical model is defined as follows.
【0022】[0022]
【数3】 図8に示すような座標系上においてS201で求められ
た経験バリオグラムをプロットし、各プロット点に対し
て、理論バリオグラム関数がフィッテイングするよう
に、当該理論バリオグラムが有する各パラメータC1、
C2、aの各値を最適化する(S203参照)。これら
のパラメータの値を特定するのは、以下に説明する理論
共分散関数に含まれるパラメータの値を決定するためで
ある。(Equation 3) The empirical variogram obtained in S201 is plotted on a coordinate system as shown in FIG. 8, and each parameter C1 of the theoretical variogram is set so that the theoretical variogram function is fitted to each plot point.
Each value of C2 and a is optimized (see S203). The values of these parameters are specified to determine the values of the parameters included in the theoretical covariance function described below.
【0023】S204では、理論バリオグラムに対応す
る理論共分散関数ρ(h)がまず決定される。つまり、
各パラメータC1、C2、aの値が代入され、以下のよ
うに当該関数が特定される。In S204, the theoretical covariance function ρ (h) corresponding to the theoretical variogram is first determined. That is,
The values of the parameters C1, C2, and a are substituted, and the function is specified as follows.
【0024】[0024]
【数4】 S205では、この理論共分散関数を用いて、各補間位
置に対する各実測値の影響度が重みWiとして計算され
る。具体的には、(Equation 4) In S205, using the theoretical covariance function, the degree of influence of each measured value on each interpolation position is calculated as a weight Wi. In particular,
【数5】 として、以下の条件を満たすように、Wiが計算され
る。ここで、W=R-1×Dである。(Equation 5) Is calculated so that the following condition is satisfied. Here, W = R −1 × D.
【0025】[0025]
【数6】 ここで、ρijは定点iと測定点j間(距離h)における
ρ(h)の値である。添え字番号の0は補間位置を表
し、1からnは実測定位置ないし実測値を表す。また、
uはラグランジュ変数といわれ、それは上記条件を満た
すように適宜決定される。(Equation 6) Here, ρ ij is the value of ρ (h) between the fixed point i and the measurement point j (distance h). The subscript number 0 represents an interpolation position, and 1 to n represent an actual measurement position or an actual measurement value. Also,
u is called a Lagrange variable, which is appropriately determined so as to satisfy the above conditions.
【0026】S206では、重みWiを使用し、各実測位
置での実測値をZi、補間位置の推定値をZ0として、下記
のようにして、推定値が演算される。[0026] In S206, using the weight W i, the measured value at each measured position Z i, the estimated value of the interpolation position as Z 0, as follows, the estimated value is calculated.
【0027】[0027]
【数7】 なお、以上のようにして求められた複数の推定値と複数
の実測値によって地理上の線量率分布が所定の表現形態
によって描かれる。(Equation 7) The geographical dose rate distribution is drawn in a predetermined expression form based on the plurality of estimated values and the plurality of actually measured values obtained as described above.
【0028】[0028]
【発明の実施の形態】以下、本発明の好適な実施形態を
図面に基づいて説明する。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Preferred embodiments of the present invention will be described below with reference to the drawings.
【0029】図1及び図2には、本発明に係る放射線測
定システムの好適な実施形態が示されている。FIGS. 1 and 2 show a preferred embodiment of the radiation measuring system according to the present invention.
【0030】図1において、この図1に示す構成例で
は、車両としての移動体8にデータ処理部10などの構
成が搭載されている。データ処理部10は例えばパーソ
ナルコンピュータで構成され、それには放射線測定器1
2及び地理上の位置を検出するGPS14が接続されて
いる。In FIG. 1, in the configuration example shown in FIG. 1, a structure such as a data processing unit 10 is mounted on a moving body 8 as a vehicle. The data processing unit 10 is composed of, for example, a personal computer,
2 and a GPS 14 for detecting a geographical position are connected.
【0031】放射線検出器12は、半導体検出器やシン
チレーション検出器などによって構成されるものであ
り、データ処理部10に対して検出値が出力される。デ
ータ処理部10はその検出値を処理することによって線
量率などの演算を行う機能を有している。この場合、G
PS14から出力される地理上の位置情報と共に線量率
が管理される。それらのデータはメモリ16上に格納さ
れ、必要に応じてそのメモリ16からデータが読み出さ
れて所定の演算が実行される。その演算結果は表示器1
8に表示される。本実施形態においては、特に表示器1
8上において後に具体的に示す地理上の線量率分布が表
示されている。The radiation detector 12 includes a semiconductor detector, a scintillation detector, and the like, and outputs a detection value to the data processing unit 10. The data processing unit 10 has a function of calculating the dose rate or the like by processing the detected value. In this case, G
The dose rate is managed together with the geographical position information output from the PS 14. These data are stored in the memory 16, and the data is read out from the memory 16 as necessary, and a predetermined operation is executed. The calculation result is displayed on display 1
8 is displayed. In the present embodiment, in particular, the display 1
8 shows a geographical dose rate distribution that will be specifically described later.
【0032】図2に示す構成例では、放射線測定システ
ムが複数個のモニタリング装置20と、ホスト装置26
とで構成されている。モニタリング装置20は地理上の
所定箇所に固定的に設置されるものであり、そのモニタ
リング装置20は放射線検出器22及び通信部24を有
している。放射線検出器22は、上記の放射線検出器1
2と同様に放射線の検出を行うものであり、その検出値
は通信部24を介してホスト装置26へ送信されてい
る。その通信は有線あるいは無線方式によって実現され
ている。ここで、各モニタリング装置20の位置があら
かじめ求められており、モニタリング装置20から出力
される検出データ(線量率データ)には位置データが付
加されている。In the configuration example shown in FIG. 2, the radiation measurement system includes a plurality of monitoring devices 20 and a host device 26.
It is composed of The monitoring device 20 is fixedly installed at a predetermined geographical location. The monitoring device 20 has a radiation detector 22 and a communication unit 24. The radiation detector 22 is the radiation detector 1 described above.
The detection of the radiation is performed in the same manner as in No. 2, and the detected value is transmitted to the host device 26 via the communication unit 24. The communication is realized by a wired or wireless system. Here, the position of each monitoring device 20 is determined in advance, and the position data is added to the detection data (dose rate data) output from the monitoring device 20.
【0033】ホスト装置26においては、通信部28に
よって各モニタリング装置からの送信データが受信さ
れ、それらの受信データがデータ処理部32によって処
理される。メモリ30上には上記のメモリ16と同様に
線量率データ及び位置データが相互に関連づけて格納さ
れる。データ処理部32はそのメモリ30上から必要な
データを読み出して所定の演算を実行し、その演算結果
を表示器34に出力する。その表示器34上には地理上
における線量率分布が表示される。In the host device 26, transmission data from each monitoring device is received by the communication unit 28, and the received data is processed by the data processing unit 32. The dose rate data and the position data are stored in the memory 30 in association with each other, similarly to the memory 16 described above. The data processing section 32 reads necessary data from the memory 30 and executes a predetermined operation, and outputs the operation result to the display 34. The display 34 displays a geographical dose rate distribution.
【0034】なお、上記の図1に示した構成例では、移
動体8上にデータ処理部10を搭載したが、データ処理
部10をホスト装置に設け、移動体8からのデータ転送
によってデータ処理を行うようにしてもよい。Although the data processing unit 10 is mounted on the mobile unit 8 in the configuration example shown in FIG. 1, the data processing unit 10 is provided in the host device, and the data processing unit 10 performs data processing by transferring data from the mobile unit 8. May be performed.
【0035】図4には、本実施形態に係る線量率分布を
描くためのソフトウエアによって作成された画像例が示
されている。ここで、画面内には地図が表示されてお
り、その地図上には移動体の走行経路上にマークとして
の測定箇所102が表されている。ちなみに、そのウイ
ンドウの上側には各種のアイコンが表示され、線量率の
単位を選択するためのボタン103が含まれる。また、
ウインドウの左側には地図の縮尺を選択するためのボタ
ン100が含まれる。FIG. 4 shows an example of an image created by software for drawing a dose rate distribution according to the present embodiment. Here, a map is displayed on the screen, and a measurement point 102 as a mark is displayed on the traveling route of the moving object on the map. Incidentally, various icons are displayed on the upper side of the window, and a button 103 for selecting a unit of the dose rate is included. Also,
The left side of the window includes a button 100 for selecting the scale of the map.
【0036】なお、図4には、一部の測定点のみが示さ
れているが、主要な各道路に沿って走行測定あるいは停
車測定を実行し、できるだけ広い範囲にわたって測定点
が分散するようにその測定を行うのが望ましい。FIG. 4 shows only a part of the measurement points. However, the traveling measurement or the stop measurement is performed along each of the main roads, and the measurement points are dispersed over as wide a range as possible. It is desirable to make that measurement.
【0037】以上のような実測定の結果に基づいて、図
5及び図6に示すような放射線分布が表示される。図5
に示す表示形態においては、いわゆる等高線表示104
によって放射線の線量率の分布が表されている。また、
図6に示す表示例においては、いわゆるカラー表示ある
いはカラーメッシュ表示106などの表示形態によって
放射線の線量率分布が表されている。ちなみに、このよ
うな表示を行う場合に、線量率分布と地図とを合成して
表示するのが望ましく、このような形態によれば、地図
上の各位置との関係において線量率を把握することが可
能となる。Based on the results of the actual measurement as described above, the radiation distribution as shown in FIGS. 5 and 6 is displayed. FIG.
In the display form shown in FIG.
Represents the distribution of the radiation dose rate. Also,
In the display example shown in FIG. 6, the dose rate distribution of the radiation is represented by a display form such as a so-called color display or a color mesh display 106. Incidentally, when such a display is performed, it is desirable to combine and display the dose rate distribution and the map. According to such a form, it is necessary to grasp the dose rate in relation to each position on the map. Becomes possible.
【0038】図3には、図1及び図2に示したデータ処
理部10,32によって実行される処理がフローチャー
トとして示されている。FIG. 3 is a flowchart showing the processing executed by the data processing units 10 and 32 shown in FIGS.
【0039】まずS101では、上述したように実測位
置において、実測値が測定される。この場合、各実測位
置についての位置データも同時に収集される。S102
では、実測位置が存在していない領域に補間位置が定め
られ、その補間位置において線量率の推定値が推定され
る。そして、S103では、図5あるいは図6に示した
ような分布マップが作成され、そしてS104において
当該分布マップが表示される。First, in S101, an actual measurement value is measured at the actual measurement position as described above. In this case, position data for each measured position is also collected at the same time. S102
In, an interpolation position is determined in a region where no actual measurement position exists, and an estimated value of the dose rate is estimated at the interpolation position. Then, in S103, a distribution map as shown in FIG. 5 or FIG. 6 is created, and in S104, the distribution map is displayed.
【0040】上記のS102の工程において推定値を推
定する場合には、上述した図9に示す手法を適用するの
が望ましい。すなわち、実測値に基づいて経験バリオグ
ラムを演算し、その演算結果に基づいて理論バリオグラ
ムを表す関数のフィッティングを実行し、そのフィッテ
ィングによりパラメータの値を特定し、そのようなパラ
メータの値によって規定される所定の関数から各実測値
ごとの重みを算出し、そのような重みに基づいて推定値
を推定する。もちろん、この手法以外の手法を利用して
推定値の推定を行ってもよいが、上記の空間統計的理論
に基づいた手法によれば、より診断性の高い線量率の推
定を行えるという利点がある。When estimating the estimated value in the step S102, it is desirable to apply the method shown in FIG. 9 described above. That is, an empirical variogram is calculated based on actual measurement values, a function representing a theoretical variogram is performed based on the calculation result, a parameter value is specified by the fitting, and the parameter value is specified by such a parameter value. A weight for each measured value is calculated from a predetermined function, and an estimated value is estimated based on such a weight. Of course, the estimation value may be estimated using a method other than this method. However, according to the method based on the above-mentioned spatial statistical theory, there is an advantage that the dose rate can be estimated with higher diagnostic performance. is there.
【0041】[0041]
【発明の効果】以上説明したように、本発明によれば、
地理上における信頼性の高い線量率分布を作成できると
いう利点がある。As described above, according to the present invention,
There is an advantage that a reliable geographic dose rate distribution can be created.
【図1】 本発明に係るシステムの構成例を示す図であ
る。FIG. 1 is a diagram showing a configuration example of a system according to the present invention.
【図2】 本発明に係るシステムの他の構成例を示す図
である。FIG. 2 is a diagram showing another configuration example of the system according to the present invention.
【図3】 データ処理部の動作内容を示すフローチャー
トである。FIG. 3 is a flowchart showing an operation content of a data processing unit.
【図4】 地理上における測定点を示す図である。FIG. 4 is a diagram showing measurement points in geography.
【図5】 等高線表示による線量率分布を示す図であ
る。FIG. 5 is a diagram showing a dose rate distribution by contour display.
【図6】 カラーメッシュ表示による線量率分布を示す
図である。FIG. 6 is a diagram showing a dose rate distribution by color mesh display.
【図7】 経験バリオグラムの関数系を示す図である。FIG. 7 is a diagram showing a functional system of an empirical variogram.
【図8】 理論バリオグラムの関数系を示す図である。FIG. 8 is a diagram showing a functional system of a theoretical variogram.
【図9】 統計的手法を基礎とする線量率の推定方法を
示すフローチャートである。FIG. 9 is a flowchart showing a method of estimating a dose rate based on a statistical method.
8 移動体、10 データ処理部、12 放射線検出
器、14 GPS、16メモリ、18 表示器、20
モニタリング装置、22 放射線検出器、24通信部、
26 ホスト装置、28 通信部、30 メモリ、32
データ処理部、34 表示器。8 moving body, 10 data processing unit, 12 radiation detector, 14 GPS, 16 memory, 18 display, 20
Monitoring device, 22 radiation detector, 24 communication unit,
26 host device, 28 communication unit, 30 memory, 32
Data processing unit, 34 display.
───────────────────────────────────────────────────── フロントページの続き Fターム(参考) 2C032 HB22 HC26 2G088 AA06 EE10 EE11 FF17 KK20 KK24 KK27 KK32 KK35 LL13 MM02 MM04 5J062 CC07 HH05 9A001 GG14 GZ05 JJ11 JJ72 ──────────────────────────────────────────────────続 き Continued on the front page F term (reference) 2C032 HB22 HC26 2G088 AA06 EE10 EE11 FF17 KK20 KK24 KK27 KK32 KK35 LL13 MM02 MM04 5J062 CC07 HH05 9A001 GG14 GZ05 JJ11 JJ72
Claims (11)
の線量率測定を行って実測値を求める測定手段と、 前記複数の実測位置についての実測値に基づいて、地理
上の複数の補間位置について線量率の推定値を推定する
推定手段と、 前記複数の実測位置についての実測値と前記複数の補間
位置についての推定値とに基づいて、地理上における線
量率分布を形成する形成手段と、 を含むことを特徴とする放射線測定システム。A measuring means for measuring a radiation dose rate at a plurality of geographically measured positions to obtain an actually measured value; and a plurality of geographically interpolated positions based on the measured values at the plurality of measured positions. Estimating means for estimating the estimated value of the dose rate, and forming means for forming a geographical dose rate distribution based on the measured values for the plurality of measured positions and the estimated values for the plurality of interpolation positions. A radiation measurement system comprising:
れた複数のモニタリング装置を含むことを特徴とする放
射線測定システム。2. The radiation measurement system according to claim 1, wherein said measurement means includes a plurality of monitoring devices fixedly installed for each of said plurality of measurement positions.
移動***置検出器を含むことを特徴とする放射線測定シ
ステム。3. The radiation measuring system according to claim 1, wherein said measuring means includes a radiation measuring instrument and a moving object position detector mounted on the moving object.
とも一方を利用して表現されることを特徴とする放射線
測定システム。4. The radiation measurement system according to claim 1, wherein the dose rate distribution is represented by using at least one of a contour line expression and a coloring expression.
とを特徴とする放射線測定システム。5. The radiation measurement system according to claim 4, further comprising: means for variably setting a density of the plurality of interpolation positions.
の線量率測定を行って実測値を求める測定手段と、 空間統計学的手法に従って、前記複数の実測位置につい
ての実測値に基づいて、地理上の複数の補間位置につい
て線量率の推定値を推定する推定手段と、 前記複数の実測位置についての実測値と前記複数の補間
位置についての推定値とに基づいて、地理上における線
量率分布を形成する形成手段と、 を含むことを特徴とする放射線測定システム。6. A measuring means for measuring a radiation dose rate at a plurality of geographically measured positions to obtain a measured value, and a geographical location based on the measured values at the plurality of measured positions according to a spatial statistical technique. Estimating means for estimating the estimated value of the dose rate for the plurality of interpolation positions above, Based on the measured values for the plurality of measured positions and the estimated values for the plurality of interpolation positions, the geographical dose rate distribution A radiation measuring system, comprising: forming means for forming.
位置の相互間距離とに基づいて、各距離ごとの経験バリ
オグラムを演算する手段と、 前記各距離ごとの経験バリオグラムに対し、理論バリオ
グラムを定義する理論バリオグラム関数がフィッティン
グするように、その理論バリオグラム関数に含まれるパ
ラメータの値を決定する手段と、 前記決定されたパラメータの値によって理論共分散関数
を決定し、前記各補間位置ごとに、前記理論共分散関数
に基づいて前記複数の実測値が当該補間位置へ与える影
響を重み値として求める手段と、 前記各補間位置ごとに、前記各実測値に対してそれに対
応する重み値を乗算した結果を加算し、これにより前記
推定値を演算する手段と、 を含むことを特徴とする放射線測定システム。7. The system according to claim 6, wherein the estimating means calculates an empirical variogram for each distance based on the measured values of the plurality of measured positions and the distance between the plurality of measured positions. Means for determining, for the empirical variogram for each distance, a value of a parameter included in the theoretical variogram function so that a theoretical variogram function defining the theoretical variogram is fitted; and Means for determining a theoretical covariance function according to the value, and for each of the interpolation positions, a means for determining, as a weight value, an influence of the plurality of actually measured values on the interpolation position based on the theoretical covariance function; And the result obtained by multiplying each of the measured values by the corresponding weight value is added thereto, thereby calculating the estimated value. Radiation measurement system characterized in that it comprises a stage, a.
ル又はガウス型モデルの1つを基礎とする関数であるこ
とを特徴とする放射線測定システム。8. The radiation measurement system according to claim 7, wherein the theoretical variogram function is a function based on one of a spherical model, an exponential model, and a Gaussian model.
まれるパラメータの値を格納する手段と、 前記パラメータの値によって理論共分散関数を決定し、
前記各補間位置ごとに、前記理論共分散関数に基づいて
前記複数の実測値が当該補間位置へ与える影響を重み値
として求める手段と、 前記各補間位置ごとに、前記各実測値に対してそれに対
応する重み値を乗算した結果を加算し、これによって前
記推定値を演算する手段と、 を含むことを特徴とする放射線測定システム。9. The system according to claim 1, wherein the estimating unit stores a value of a parameter included in a theoretical variogram function that defines a theoretical variogram, and determines a theoretical covariance function based on the value of the parameter. ,
For each of the interpolation positions, means for determining, as a weight value, the influence of the plurality of measured values on the interpolation position based on the theoretical covariance function, for each of the interpolation positions, Means for adding the result of multiplying by the corresponding weight value and calculating the estimated value based on the result.
線の線量率測定を行って実測値を求める工程と、 空間統計学的手法に従って、前記複数の実測位置につい
ての実測値に基づいて、地理上の複数の補間位置につい
て線量率の推定値を推定する工程と、 前記複数の実測位置についての実測値と前記複数の補間
位置についての推定値とに基づいて、地理上における線
量率分布を形成する工程と、 を含むことを特徴とする放射線測定方法。10. A step of measuring a radiation dose rate at a plurality of geographically measured positions to obtain a measured value; and a method of calculating a geographical position based on the measured values at the plurality of measured positions according to a spatial statistical technique. Estimating an estimated value of the dose rate for the plurality of interpolation positions, and forming a geographical dose rate distribution based on the measured values of the plurality of actually measured positions and the estimated values of the plurality of interpolation positions. A radiation measurement method, comprising:
グラムを格納した媒体であって、 前記プログラムは、 地理上の複数の実測位置において放射線の線量率測定を
行って得られた実測値を入力する機能と、 空間統計学的手法に従って、前記複数の実測位置につい
ての実測値に基づいて、地理上の複数の補間位置につい
て線量率の推定値を推定する機能と、 前記複数の実測位置についての実測値と前記複数の補間
位置についての推定値とに基づいて、地理上における線
量率分布を形成する機能と、 を含むことを特徴とするプログラム記録媒体。11. A medium storing a program used in a radiation measurement system, the program having a function of inputting a measured value obtained by measuring a radiation dose rate at a plurality of geographically measured positions. And a function of estimating an estimated value of the dose rate for a plurality of geographically interpolated positions based on the actually measured values of the plurality of actually measured positions, according to a spatial statistical method, and an actual measured value of the plurality of actually measured positions. And a function of forming a geographical dose rate distribution based on the estimated values of the plurality of interpolation positions.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000019152A JP2001208848A (en) | 2000-01-27 | 2000-01-27 | Radiation mesuring system and method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000019152A JP2001208848A (en) | 2000-01-27 | 2000-01-27 | Radiation mesuring system and method |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2001208848A true JP2001208848A (en) | 2001-08-03 |
Family
ID=18545903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000019152A Pending JP2001208848A (en) | 2000-01-27 | 2000-01-27 | Radiation mesuring system and method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2001208848A (en) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004205375A (en) * | 2002-12-25 | 2004-07-22 | Tokyo Electric Power Co Inc:The | Movable outdoor radiation measuring system |
JP2004264085A (en) * | 2003-02-28 | 2004-09-24 | Aloka Co Ltd | Radiation monitoring device |
WO2006090634A1 (en) * | 2005-02-22 | 2006-08-31 | National Institute Of Radiological Sciences | Method for predicting radiation detection value and prediction response radiation detector |
US7723696B2 (en) | 2005-02-22 | 2010-05-25 | National Institute Of Radiological Sciences | Radiation detected value forecasting method and forecast responsive radiation detector |
JP2010266304A (en) * | 2009-05-14 | 2010-11-25 | Genden Information System Co Ltd | System and method for monitoring of environmental radiation in emergency |
JP2011528116A (en) * | 2008-07-14 | 2011-11-10 | イーツェーイクス・テクノロジーズ・ゲーエムベーハー | Detector system for radiometry |
JP2012504749A (en) * | 2008-07-14 | 2012-02-23 | イーツェーイクス・テクノロジーズ・ゲーエムベーハー | Detector system with position determination system |
JP2013032926A (en) * | 2011-08-01 | 2013-02-14 | Kyoto Univ | Method for measuring radiation dose rate and method for creating radiation dose rate map |
JP2013092438A (en) * | 2011-10-26 | 2013-05-16 | Sumitomo Mitsui Construction Co Ltd | Method of removing radioactive material-contaminated soil, and system of removing and managing radioactive material-contained soil |
KR101308937B1 (en) | 2012-01-16 | 2013-09-23 | 한국원자력안전기술원 | System and method for car-borne survey of radioactivity |
JP2014095677A (en) * | 2012-11-12 | 2014-05-22 | Masahiko Tsujita | Cleaning method for radioactive material and cleaning apparatus used for the same |
JP2014145700A (en) * | 2013-01-30 | 2014-08-14 | Japan Atomic Energy Agency | Method for measuring deposition quantity of radioactive cesium in vicinity of nuclear power facility |
JP2014145628A (en) * | 2013-01-28 | 2014-08-14 | Mitsubishi Heavy Ind Ltd | System and method for creating radioactive substance distribution map |
JP2015526702A (en) * | 2012-06-12 | 2015-09-10 | アンスティテュ ド ラディオプロテクシオン エ ド スルテ ニュクレエール | How to map photon distribution in real time on site |
JP2016075710A (en) * | 2016-01-26 | 2016-05-12 | 大成建設株式会社 | Estimation method for amount of space radiation and entrance and exit management system |
JP2017016490A (en) * | 2015-07-03 | 2017-01-19 | 富士電機株式会社 | Monitoring system |
KR20180031984A (en) * | 2016-09-21 | 2018-03-29 | 한국원자력연구원 | Disaster preventer trainning system, portable terminal for trainning disaster preventer and method for operating disaster preventer trainning system |
JP2018087767A (en) * | 2016-11-29 | 2018-06-07 | 地方独立行政法人東京都立産業技術研究センター | Radiation dosimetry device, radiation dosimetry method and radiation dosimetry system |
KR20180070944A (en) * | 2016-12-19 | 2018-06-27 | (주)인스페이스 | Method for prediction the solar proton particle |
CN108701338A (en) * | 2017-01-13 | 2018-10-23 | 拉德搜公司 | The emergent in-situ measurement training method of radiation based on virtual accident |
-
2000
- 2000-01-27 JP JP2000019152A patent/JP2001208848A/en active Pending
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004205375A (en) * | 2002-12-25 | 2004-07-22 | Tokyo Electric Power Co Inc:The | Movable outdoor radiation measuring system |
JP2004264085A (en) * | 2003-02-28 | 2004-09-24 | Aloka Co Ltd | Radiation monitoring device |
US7723696B2 (en) | 2005-02-22 | 2010-05-25 | National Institute Of Radiological Sciences | Radiation detected value forecasting method and forecast responsive radiation detector |
JPWO2006090634A1 (en) * | 2005-02-22 | 2008-07-24 | 独立行政法人放射線医学総合研究所 | Prediction method of radiation detection value and prediction response type radiation detector |
WO2006090634A1 (en) * | 2005-02-22 | 2006-08-31 | National Institute Of Radiological Sciences | Method for predicting radiation detection value and prediction response radiation detector |
JP2011528116A (en) * | 2008-07-14 | 2011-11-10 | イーツェーイクス・テクノロジーズ・ゲーエムベーハー | Detector system for radiometry |
JP2012504749A (en) * | 2008-07-14 | 2012-02-23 | イーツェーイクス・テクノロジーズ・ゲーエムベーハー | Detector system with position determination system |
JP2010266304A (en) * | 2009-05-14 | 2010-11-25 | Genden Information System Co Ltd | System and method for monitoring of environmental radiation in emergency |
JP2013032926A (en) * | 2011-08-01 | 2013-02-14 | Kyoto Univ | Method for measuring radiation dose rate and method for creating radiation dose rate map |
JP2013092438A (en) * | 2011-10-26 | 2013-05-16 | Sumitomo Mitsui Construction Co Ltd | Method of removing radioactive material-contaminated soil, and system of removing and managing radioactive material-contained soil |
KR101308937B1 (en) | 2012-01-16 | 2013-09-23 | 한국원자력안전기술원 | System and method for car-borne survey of radioactivity |
US9678228B2 (en) | 2012-06-12 | 2017-06-13 | Institut De Radioprotection Et De Surete Nuclaire | Method of real-time mapping of a distribution of photons in a site |
JP2015526702A (en) * | 2012-06-12 | 2015-09-10 | アンスティテュ ド ラディオプロテクシオン エ ド スルテ ニュクレエール | How to map photon distribution in real time on site |
JP2014095677A (en) * | 2012-11-12 | 2014-05-22 | Masahiko Tsujita | Cleaning method for radioactive material and cleaning apparatus used for the same |
JP2014145628A (en) * | 2013-01-28 | 2014-08-14 | Mitsubishi Heavy Ind Ltd | System and method for creating radioactive substance distribution map |
US9990751B2 (en) | 2013-01-28 | 2018-06-05 | Mitsubishi Heavy Industries, Ltd. | Radioactive substance distribution map producing system and method of producing radioactive substance distribution map |
JP2014145700A (en) * | 2013-01-30 | 2014-08-14 | Japan Atomic Energy Agency | Method for measuring deposition quantity of radioactive cesium in vicinity of nuclear power facility |
JP2017016490A (en) * | 2015-07-03 | 2017-01-19 | 富士電機株式会社 | Monitoring system |
JP2016075710A (en) * | 2016-01-26 | 2016-05-12 | 大成建設株式会社 | Estimation method for amount of space radiation and entrance and exit management system |
KR20180031984A (en) * | 2016-09-21 | 2018-03-29 | 한국원자력연구원 | Disaster preventer trainning system, portable terminal for trainning disaster preventer and method for operating disaster preventer trainning system |
KR101938550B1 (en) * | 2016-09-21 | 2019-04-10 | 한국원자력연구원 | Disaster preventer trainning system, portable terminal for trainning disaster preventer and method for operating disaster preventer trainning system |
JP2018087767A (en) * | 2016-11-29 | 2018-06-07 | 地方独立行政法人東京都立産業技術研究センター | Radiation dosimetry device, radiation dosimetry method and radiation dosimetry system |
KR20180070944A (en) * | 2016-12-19 | 2018-06-27 | (주)인스페이스 | Method for prediction the solar proton particle |
KR101890494B1 (en) * | 2016-12-19 | 2018-08-21 | (주)인스페이스 | Method for prediction the solar proton particle |
CN108701338A (en) * | 2017-01-13 | 2018-10-23 | 拉德搜公司 | The emergent in-situ measurement training method of radiation based on virtual accident |
JP2019507893A (en) * | 2017-01-13 | 2019-03-22 | ラドサーチ カンパニー | On-site exploration training method of virtual accident base radiation emergency |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2001208848A (en) | Radiation mesuring system and method | |
Ikechukwu et al. | Accuracy assessment and comparative analysis of IDW, spline and kriging in spatial interpolation of landform (topography): an experimental study | |
US10514696B2 (en) | Navigation driving metric | |
US8874367B2 (en) | Method for estimating and displaying range of a vehicle | |
US7233859B2 (en) | Method and device for planning a trajectory | |
US8152366B2 (en) | Estimation of subsurface thermal structure using sea surface height and sea surface temperature | |
US20180010936A1 (en) | River discharge and depth estimation | |
JPWO2018181974A1 (en) | Judgment device, judgment method, and program | |
CN106341789A (en) | A mobile computer atmospheric barometric pressure system | |
JP6159594B2 (en) | Observation value processing apparatus, observation value processing method, and observation value processing program | |
CN102984714B (en) | A kind of method for detecting area based on radio broadcast coverage data and equipment | |
CN112762935B (en) | Positioning method and system based on hull attitude monitoring | |
KR101599613B1 (en) | Positioning apparatus and positioning method | |
KR20190124625A (en) | Method, apparatus, and system for analyzing vertical distribution of particle material using satellite sensor and surface weather observation | |
JP2019174191A (en) | Data structure, information transmitting device, control method, program, and storage medium | |
JP5871536B2 (en) | Position specifying device, position specifying method, and program | |
WO2019188820A1 (en) | Information transmission device, data structure, control method, program, and storage medium | |
KR20010008320A (en) | Apparatus and method of satellite navigation using sequential kalman filter algorithm | |
JP2016057438A (en) | Travel evaluation device, travel estimation method, and travel estimation program | |
JP4210763B2 (en) | Method and apparatus for continuous positioning of moving body using both wireless LAN positioning and GPS positioning, and continuous positioning program for moving body | |
JP3589186B2 (en) | Onboard marine forecasting device | |
CN118225180B (en) | Meteorological station headroom environment analysis method, equipment and storage medium | |
KR102450271B1 (en) | Apparatus and method for estimating a signal propagation model or estimating the location of a base station | |
JP2020160014A (en) | Position estimation device, position estimation method and position estimation program | |
JP2019174675A (en) | Data structure, map data generator, control method, program, and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20031210 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040106 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040308 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20040308 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20040406 |