エアカウンターSをいろいろ調べてみた

自分まとめ。エステー・エアカウンターSについて任意取り調べ中。
38

感度

線量計には感度があって、cpm/(μSv/h)という単位で表される。つまり1μSv/hの所で1分間に何本放射線を捉えられるか。

エアカウンターSは廉価だけあって感度はかなり低めだと思われるのだけど、どれくらいか調べてみたのが下のツイート。

ちなみに放射性物質はバカなので記憶を持っていない。自分がどの位長く放射性物質をやっているか奴らは知らない。だから放射性物質から放射線が出てくるタイミングはまったくランダム。しかつめらしくポアソン過程と言っても同じこと。

30秒なら30秒、ポアソン過程のものを繰り返しながめていると、出てくる放射線の数はポアソン分布なんつーものにしたがう。このポアソン分布だと放射線の数の平均と分散(値ひく平均を2乗したものの平均)が一致するという性質がある。そこから感度が調べられるさ。

  感度 = 指示値の平均 / (指示値の分散×測定時間)

畠山元彦 @MuiMuiZ

エアカウンターSきたさ。最初の30秒の値をひたすら書き留めて、ポアソンだと仮定してγ線に対する感度を求めてみたら 38.3 cpm/(μSv/h) でした。0.10μSv/hだと2分の計測でγ線7~8本。なかなかしんどそう。 http://t.co/7zz4orxC

2012-02-22 19:34:27

30 秒の測定だと値が飛び飛びなのだけど、それが放射線 1 本分の違いとしてそこから推測しても 40 cpm/(μSv/h) ほどとなる。

んで、大体、捉えた放射線の数で一回一回の誤差が決まる。

感度が低ければ誤差が大きい。単純に放射線の数を数え、時間あたりの平均を取るなら、これは単に捉えた放射線の本数 n で決まって、標準偏差は放射線の本数の平方根に反比例 1/√n する。例えば、100 本なら 1/√100 = 1/10 = 10%。

(ここで言う誤差は、放射線がランダムに出てくるため、あるときには多かったり、別のときには少なかったりすることによる誤差。たとえ十分に放射線を捉えていて、この意味での誤差が小さくても、μSv/h単位での値には依然±20%の誤差があるのだと仕様書はしている)

1 μSv/h なら 1 分測定してガンマ線が平均 40 本捉えられる。よって平均 100 本にするには 2.5 分。0.1 μSv/h のような低いところだと 1 分あたり平均 4 本なので、25 分測定する必要がある。でもそんなに長く測定して平均してくれるの?

値の平均方法

そういうわけで、精度よく調べたければ長く測定して捉える放射線の数を増やせばいいのだけれど、そもそもエアカウンターSが長い測定で捉えた数の平均をきっちり出してくれなければ意味がない。

でも、そんな平均を出していては値の変化に対応できないから、たぶんそんなことはしていない。どのくらい「していない」かを調べたのが、その次のツイートさ。

畠山元彦 @MuiMuiZ

エアカウンターSがどんなふうに値を平均化しているか調べてみたさ。0.67μSv/hのところにおいた後、0.18μSv/hのところへ移動させたときの値の変化。 http://t.co/j3up7Wom

2012-02-28 21:43:56
拡大
畠山元彦 @MuiMuiZ

(エアカウンターS)それで、ここからあるときの指示値がどれだけ前の値に影響されているかを調べられる。 http://t.co/TVQprUAR

2012-02-28 21:45:26
拡大
畠山元彦 @MuiMuiZ

(エアカウンターS)移動平均なら矩形、記憶なしに指数的に影響が減るなら指数関数になるはずだけど、なにか山型か台形のようになった。各々のパスを見るとすぐに変化するものと暫く前の値を引きずるものがあるので平均としての台形かもしれない。3分経つとほぼ影響はなくなる

2012-02-28 21:48:28
畠山元彦 @MuiMuiZ

(エアカウンターS)2分ぐらいの移動平均で、範囲が変わるのに最大1分ぐらいかかるのかもしれない。ともかく同じ箇所で3分以上連続して測定しても精度はよくならないことがわかる。精度を求めるなら2分に1回ぐらいリセットして平均するのが正解か

2012-02-28 21:50:53

そういうわけで、0.1 μSv/h のような低いところで測定ごとのバラツキを 1 割程度に抑えたければ、2~3 分の測定を 10 回ほどやって平均を取るのがよさそう。

分解されてる方がいたよ

うまぷ @umapu

@MuiMuiZ 使ってるMCUのRAMは512Byteみたいですね。どういう計算してるかわかりませんが、RAM容量以上は保持できないからそれ以下のサンプル数ってことでしょうか。http://t.co/KAmbyTUe

2012-02-28 21:54:35
畠山元彦 @MuiMuiZ

@umapu おぉ、ありがとう。512byteなんですか…。素直に指数関数にすれば、直前の値だけでいいはずなのだけど。

2012-02-28 22:01:58
うまぷ @umapu

@MuiMuiZ センサが4個繋がってるようなので、なんか難しい計算してるんですかね・・・

2012-02-28 22:12:33
畠山元彦 @MuiMuiZ

@umapu うーん、何か難しいことやる必然性はないとおもうのだけどな。

2012-02-28 22:19:08