JP3958069B2 - 放射線測定装置 - Google Patents
放射線測定装置 Download PDFInfo
- Publication number
- JP3958069B2 JP3958069B2 JP2002046788A JP2002046788A JP3958069B2 JP 3958069 B2 JP3958069 B2 JP 3958069B2 JP 2002046788 A JP2002046788 A JP 2002046788A JP 2002046788 A JP2002046788 A JP 2002046788A JP 3958069 B2 JP3958069 B2 JP 3958069B2
- Authority
- JP
- Japan
- Prior art keywords
- pulse
- value
- power
- output
- 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.)
- Expired - Lifetime
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/17—Circuit arrangements not adapted to a particular type of detector
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、放射線測定装置に関する。より詳細には、本発明は、放射線測定装置により放射線測定を行う場合に広レンジの監視に用いられるパルス計測およびキャンベル計測に関し、そのディジタル信号処理による耐ノイズ性の改善に関する。
【０００２】
【従来の技術】
広いレンジの放射線を測定するためにパルス計測およびキャンベル計測が併用されている。パルス計測は、一般に放射線センサからのパルス数をカウントするものであり、そのパルスが重なりあって計数できなくなった場合に、キャンベル計測が行われる。たとえば、原子力の出力監視では、原子炉圧力容器内部に６〜１０本の起動領域センサ（ＳＲＮＭセンサ）および１００〜２００本の局所出力領域センサ（ＬＰＲＭセンサ）を設置し、それらの出力をそれぞれ起動領域モニタおよび出力領域モニタによって測定し、原子炉の出力を約１１桁の監視幅で監視している。
【０００３】
このうち起動領域モニタは、原子炉出力の低い領域、つまり原子炉出力が定格出力の１０−９％から１０−４％まではＳＲＮＭセンサの出力パルスの個数を計数すること（以後、パルス計測と呼ぶ）により出力を監視する。一方、原子炉出力が高い領域、つまり原子炉出力が１０−５％〜１０％では、センサの出力パルスの重なりにより生じる揺らぎのパワーを測定すること、すなわちキャンベル計測により原子炉の出力を監視する。
【０００４】
以下、原子炉起動監視装置におけるパルス計測およびキャンベル計測の先行技術例（特開２０００−１６２３６６号公報参照）について図１８を参照して説明する。
【０００５】
図１８に示した原子炉監視装置は、原子炉内の中性子に感応して中性子量に対応する電気パルスを出力するＳＲＮＭセンサ１と、その出力パルス信号を増幅し整形するアナログ式プリアンプ２と、ＳＲＮＭセンサ１の出力パルス幅より狭い時間間隔でサンプリングしディジタル値に変換するＡ／Ｄ変換器３と、このサンプリングされたサンプルデータからパルス数を計数し、原子炉の低い領域での出力レベル値に変換するパルス計数手段（ＰＣ）２３と、測定精度を向上させるためにＡ／Ｄ変換器３のサンプル値を積算する積算手段２４と、その積算値を２乗してパワーを演算するパワー演算手段２５と、そのパワーを平均化する加算平均手段２６と、パルス計数手段２３の計数結果および加算平均手段２６の演算結果に基づいて原子炉起動時の出力を連続監視する原子炉出力評価手段２７とから構成されている。
【０００６】
このように構成されたディジタル式原子炉起動監視装置においては、ＳＲＮＭセンサ１の出力パルスの形状をプリアンプ２で増幅整形し、増幅整形されたパルスをＡ／Ｄ変換器３が高速にサンプリングし、そのサンプル値に対してある一つの、または、複数のロジック演算を行い、その演算結果がそれぞれ予め設定した範囲内にある場合にセンサ出力パルスとしてパルス計数手段２３で計数する。
【０００７】
一方、同じサンプル値を積算手段２４で積算処理し、キャンベル計測に必要なサンプリングレートまで落とすとともに、ダイナミックレンジを稼ぐために積算し、等価ビット数の改善を計る。その結果をパワー演算手段２５でバンドパスフィルタ処理を行ったのち、その２乗和を演算し、その演算結果を加算平均手段２６で平均し、キャンベル出力値を算出する。これらパルス計数値とキャンベル出力値は原子炉出力評価手段２７で評価され、原子炉出力として表示される。
【０００８】
この構成では、パルス計数手段２３によりパルスの波高のみではなく、パルス幅の情報も取り入れた識別により広いパルス幅のノイズを除去し、センサ出力パルスのみを計数することができる。つまり、図１８の原子炉起動監視装置では、たとえば、１００ｎｓパルス幅のＳＲＮＭセンサ１の出力を２５ｎｓ間隔でサンプリングする。このサンプリングによるサンプルデータのうち、パルス幅に相当するｋ−３番めからｋ番めまでの４つのサンプル値[S(k), S(k-1), S(k-2), S(k-3)]を用い、パルス立ち上がり部相当のサンプル値S(k-3)、立ち下がり部相当のサンプル値S(k)と、その間の２つのサンプル値S(k-2)およびS(k-1)を用いて下記のような演算を行い、その演算結果(Out(k))をパルス識別の指標とし、それが予め設定したレベルの範囲内にある場合、中性子パルスとして計数する。
【０００９】
この演算を行うことにより、ＳＲＮＭセンサ１の出力パルスとほぼ同じパルス幅の信号のみを計数することが可能となる。つまり、大きなサージ状のノイズに信号パルスが重畳された場合でも、パルスのグランドレベルを差し引くことにより正確な測定値を計数することが可能となる。
【００１０】
なお、このようなセンサパルス形状に合致したケースを検出できる指標を複数設けて、それらの論理積（ＡＮＤ）をとることによって、この識別性能はさらに向上する。これにより、最もＳＲＮＭに混入しやすいとされる数μｓ幅のサージノイズが混入しても、それをほぼ完全に除去し、１００ｎｓ幅のセンサパルスのみを計数可能とすることができる。
【００１１】
一方、キャンベル計測に関しては、パワー演算手段２５により周波数帯域を制限し、２乗平均値を演算する。この構成では、ソフトウエアにより周波数帯域を設定できるため、ノイズが測定帯域と同等の周波数となった場合、測定帯域をソフトウエア上で変更することによりノイズの誘導を低減することができる。
【００１２】
【発明が解決しようとする課題】
しかし、上述の先行技術による原子炉起動監視装置においては、以下に述べるような課題が存在している。
【００１３】
第１の課題は、両極性ノイズの低減の問題である。つまり、数μｓ幅のサージノイズとセンサパルスとが重なった場合、不感時間を設けずにその重なったセンサパルスを計数するには、差分を用いてパルス波高値相当の値を算出する必要がある。しかし、この差分をとった場合、センサ出力などの正のみ、または負のみの単極性のパルスは、従来のグランドレベルからの波高値と同等の識別レベルとなるが、回路からのホワイトノイズなどの両極性のパルスに対しては、そのピーク間の電圧をパルス識別レベルと識別する。このため、従来のグランドレベルからの波高で識別する場合に比べ、約２倍の識別レベルが必要となる。したがって、センサ出力のみを計数するのに必要な識別レベルが、従来の方式より２倍程度必要となり、センサ信号／ホワイトノイズ比（ＳＮ比）が悪くなる。
【００１４】
第２の課題は、キャンベル計測における耐ノイズ性の向上である。従来のモータ等からのノイズ試験の結果、数μｓ幅のサージノイズが原子炉起動監視装置に誘導しやすいことが分かっている。パルス計測に関しては、前記のディジタル演算によるパルス識別によりこのサージノイズの低減が可能である。しかし、キャンベル計測においては、測定帯域をセンサのパルス形状から選定した数百Ｈｚから１ＭＨｚ以下の周波数帯域に設定しているが、前記の先行技術例では、この測定帯域をシフトすることにより誘導ノイズを除去している。しかし、最も誘導しやすいノイズの周波数が、ほぼ測定帯域の範囲になっており、完全に除去することは難しく、また、センサ感度がわずかに変化するため、これらを補正する必要がある。
【００１５】
一方、一般に線量当量を求める計測装置は、γ線エネルギーに対する感度特性を人体に対するエネルギー吸収特性と等しくなるように、センサの放射線入射窓、反応体積などが最適化されている。しかし、γ線の入射方向によって感度特性が異なる場合があり、感度特性を正確に一致させることは難しい。また、人体に対する線量当量を正確に換算するには、人体の部位によってエネルギー吸収特性が異なるため、γ線エネルギーに対するセンサ感度と一致させただけの計測装置のみでは、人体の各部位に対する正確な評価は難しい。さらに、γ線以外の中性子、たとえばβ線などが混在している場合、センサ構造で感度を補正しているセンサでは、物質への吸収特性が大きく異なるこれらの放射線に対しては評価できないため、それぞれ専用の測定系が必要となる。従来、これらの課題を解決するために、γ線のエネルギースペクトルを求め、その値から換算する手法が提案され、実用化されている。しかし、この手法ではパルス波高を用いたエネルギー情報の取得を元にしており、パルスのパイルアップが生じる条件では、エネルギー情報の取得が難しくなり、精度が悪くなる。つまり、センサの出力パルス幅にもよるが、通常エネルギー測定を行う場合約１×１０５ｃｐｓ程度が計測の上限となる。計測下限は、応答要求を満たす計数とする必要があるが、１ｃｐｓと仮定しても、計測レンジは約５桁程度となり、さらに広レンジの線量を連続で行える測定方法が望まれている。
【００１６】
特開平３−１８３９８３号公報には、センサを二重構造にすることにより測定の精度を改善することが記載されている。この手法では、前記したパルス計測のパイルアップの影響は、電流計測を行うことによって回避されているが、この手法に関しても、センサ構造または処理内容が複雑なため、その簡素化が望まれていた。
【００１７】
本発明はかかる従来技術に対処してなされたものであり、ディジタル処理を用いたパルス計測およびキャンベル計測において、耐ノイズ性を改善する放射線計測装置を提供することを目的とする。
【００１８】
【課題を解決するための手段】
上記目的を達成するため、請求項１に係る発明の放射線測定装置は、放射線に感応してパルス信号を出力する放射線検出手段と、この放射線検出手段の出力パルス信号を整形しサンプル値に変換するＡ／Ｄ変換手段と、サンプル値のｎ乗値相当（ｎ＞＝２）を演算し放射線のパルス信号のみを識別するｎ乗パルス識別手段と、このｎ乗パルス識別手段によって識別された放射線のパルス信号を計数する計数手段とを備え、しかも前記ｎ乗パルス識別手段は、前後ｎ個のサンプル値を掛け算することにより前記ｎ乗値相当を演算するか、又は同じサンプル値をｎ乗した後に前後２個のｎ乗値を平均化処理することにより前記ｎ乗値相当を演算するものであることを特徴とする。
【００１９】
請求項２に係る発明は、請求項１に記載の放射線測定装置において、ｎ乗パルス識別手段は、放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、パルス立ち上がり部およびパルス立ち下がり部に対応する両サンプリング点におけるサンプル値の和のｎ乗値を、両サンプリング点の間に位置する両サンプリング点におけるサンプル値の和のｎ乗値から差し引くことにより出力パルスの波高に相当する差を演算し、その差から所望のパルスのみを識別する差分識別手段を備えていることを特徴とする。
【００２０】
請求項３に係る発明は、請求項１に記載の放射線測定装置において、ｎ乗パルス識別手段は、放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、前半２つのサンプリング点におけるサンプル値の差のｎ乗値と、後半２つのサンプリング点におけるサンプル値の差のｎ乗値との和を演算し、その和に基づいて所望のパルスのみを識別する差分識別手段を備えていることを特徴とする。
【００２１】
請求項４に係る発明は、請求項１に記載の放射線測定装置において、ｎ乗パルス識別手段は、放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、各サンプル値のｎ乗値の積算値と放射線検出手段の出力パルスの幅に相当するサンプル値の波高との比率を算出し、その比率から放射線検出手段の出力パルスを識別する波高・パワー識別手段を備え、計数手段は波高・パワー識別手段によって識別されたパルスを計数することを特徴とする。
【００２２】
請求項５に係る発明は、請求項１に記載の放射線測定装置において、正および負の両極性の入力信号を単極性の信号に変換する単極性変換手段をさらに具備し、かつｎ乗パルス識別手段は、単極性変換手段によって得られた単極性の信号を偶数ｎを用いてｎ乗演算を行うｎ乗演算手段と、このｎ乗演算手段の出力を平均化する平滑手段とを備えていることを特徴とする。
【００２３】
請求項６に係る発明の放射線測定装置は、放射線に感応してパルス信号を出力する放射線検出手段と、この放射線検出手段の出力パルスを整形し周波数帯域を制限する帯域制限手段と、この帯域制限手段の各出力パルスのｎ乗値（ｎ＞＝２）を算出するｎ乗値演算手段と、このｎ乗値演算手段の出力信号を所定の時間幅で平滑化する第１の平滑化手段と、この第１の平滑化手段の出力の大小関係を評価し、その評価結果から特定のデータを除くデータ除去平均化手段と、このデータ除去平均化手段からの出力をさらに平滑化する第２の平滑化手段と、この第２の平滑化手段の出力を放射線強度に換算する変換手段とを具備する。
【００２４】
請求項７に係る発明は、請求項６に記載の放射線測定装置において、第１の平滑化手段の平滑時間幅を外来ノイズのパルス幅以上に設定し、データ除去平均化手段の平滑化区間を外来ノイズパルスの到来する時間間隔以下に設定し、計測されたノイズ波形によりこれら平滑時間幅とデータ除去数を調整するノイズ特性評価手段をさらに具備したことを特徴とする。
【００２５】
請求項８に係る発明は、請求項６に記載の放射線測定装置において、データ除去平均化手段は、第１の平滑化手段の過去の出力と現在の出力とを比較し、現在の出力値が第１の平滑化手段の揺らぎから過去の出力値をもとに事前に評価したしきい値以上に増加していた場合、第１の平滑化手段の現在の出力を除去し、予め評価によって求めた評価比率を過去の出力値に乗じた値に置き換えるデータ制限手段を有することを特徴とする。
【００２６】
請求項９に係る発明の放射線測定装置は、ｎは３以上の奇数であり、ｎ乗値演算手段は帯域制限手段の各出力パルスに対して奇数次のモーメントを算出することを特徴とする。
【００３０】
【発明の実施の形態】
＜第１の実施の形態＞
本発明の第１の実施の形態による放射線測定装置を、図１を参照して説明する。
【００３１】
図１に示す放射線測定装置は、放射線に感応し放射線量に応じた出力パルスを発生するＳＲＮＭセンサ１と、その出力パルスを増幅するプリアンプ２Ａと、プリアンプ２Ａの出力パルスをそのパルス幅より狭い時間間隔でサンプリングしサンプル値を得るＡ／Ｄ変換器３と、そのサンプル値をｎ乗しＳＲＮＭセンサ１の出力パルス幅に相当するｎ乗値を用いて信号を識別するｎ乗パルス識別装置４と、このｎ乗パルス識別装置４により識別されたパルスを計数するパルス計数装置（ＰＣ）５とを備えている。ＳＲＮＭセンサ１は、パルス出力が得られる核***センサであるが、これは同様のパルス出力が得られる電離箱等に置換することもできる。
【００３２】
このように構成された放射線測定装置のＳＲＮＭセンサ１に中性子が入射し、内部で核***すると図２（ａ）に示すような電気パルスが出力される。このＳＲＮＭセンサ１の出力パルスの時間幅は１００ｎｓ程度である。この出力パルスは、プリアンプ２Ａに入力され、増幅される。プリアンプ２Ａは、ＳＲＮＭセンサ１へ動作電圧を印加する機能をも有する。プリアンプ２Ａの出力パルスは、Ａ／Ｄ変換器３に入力され、図２（ａ）で点により示されるようなサンプリング時間間隔でサンプリングされ、ディジタル化される。このサンプリング時間間隔は短いほど波形情報を多く抽出でき、他の誤差信号、たとえば外来ノイズによる信号を除いてセンサ出力パルスのみを正確に計数することが可能となる。Ａ／Ｄ変換器３では、サンプリング前に、サンプリング定理で必要とされる周波数帯域に制限するバンドパスフィルタ処理も行う。
【００３３】
Ａ／Ｄ変換器３から出力されるサンプル値は、ｎ乗パルス識別装置４に内蔵されるｎ乗演算部に入力され、ここでｎ乗演算される。つまり、各サンプル値をそれぞれｎ乗演算するか、または、前後ｎ個のサンプル値を掛け算する。ここで、ｎ＞＝２とする。同じサンプル値をｎ乗演算した場合は、前後２個のｎ乗値を平均化処理する。たとえば、８個でパルス波形をサンプリングした場合は、おのおのの２乗値を演算した後、前後２個の演算値を平均化処理し、結果的に４つの演算値を得る。この場合、移動平均処理を行い、８個のサンプル値を得るようにしてもよい。
【００３４】
例えばｎ＝２として、２乗演算を行った場合のサンプル値の変化の様子を図２（ｂ）に示す。このようにサンプル値をｎ乗値に変換することにより、パルス波形の波高と回路によるノイズ成分の比率は、従来の比率のｎ乗倍に改善することができる。
【００３５】
ただし、単純にサンプル値ごとにｎ乗演算した値をパルス識別に用いても、サンプル値をそのままパルス識別に用いた場合と同じ識別性能となる。そこで、ディジタル演算でｎ乗する場合は、すでに述べたように前後のデータを掛け合わせるか、ｎ乗後に前後のデータを平均化する処理を加えることが必ず必要となる。この２乗演算された結果を、予め設定しておいた識別レベルと比較することにより、ＳＲＮＭセンサ１の出力パルスと回路ノイズとの識別が容易になる。この放射線測定装置によりＳＲＮＭセンサ１のパルスとして識別されたものは、パルス計数装置６でパルス発生率に換算され、最終的にはＳＲＮＭセンサ１の位置での中性子束レベルに換算される。
【００３６】
この実施の形態によれば、回路ノイズレベルと同程度の波高を有するパルスを従来の差分を用いた手法よりも良好に識別することができる。図３（ａ）にホワイトノイズと、それと同程度のパルス波高を有するセンサ出力パルスをＡ／Ｄ変換器３でサンプリングした場合の一例を示す。この図３（ａ）では、垂直破線で囲まれた時間４．２０×１０３ｎｓの近傍においてセンサ出力パルスが発生した場合を示しており、この区間で回路ノイズとセンサパルスとが重畳している。このデータに対して、従来のパルス波高によるノイズ識別を行った場合、図３（ａ）中の水平破線を識別レベルLbとすると、両パルスとも計数してしまい、正確な測定をすることはできない。また、上記の差分を用いたパルス計数方式では、図３（ａ）中、垂直破線で囲まれた範囲における正負両ピーク間の電圧がパルスの波高と認識される。したがって、正と負の両極性で生ずる回路ノイズは、従来の０Ｖからの波高値で識別していたものに比べ、２倍の電圧まで識別レベルをあげないと識別することができない。一方、加算平均化することにより、両極性のノイズは正負の信号がキャンセルされて平均化され、単極性のセンサ出力はもともと正のみ又は負のみの信号であるため、このようにキャンセルされることはないが、図３（ｂ）に示すように、パルス幅が広くなる。パルス幅が広くなった場合、センサ出力のパルス数の多い条件ではパルスが重なりパルスを計数できなくなり、パルス計測の計数上限が低下するという問題がある。
【００３７】
本実施の形態のｎ乗パルス識別装置４によるｎ乗演算処理後のサンプル値を図３（ｃ）に示す。回路ノイズ相当のホワイトノイズ（左側）の識別レベルとセンサ出力（右側）の識別レベルの比率は、従来の差分を用いた場合より大きく改善することができる。また、パルス幅も、平均化のみの処理をした場合の図３（ｂ）に比べ広がることもなく、これにより、パルス計数の上限を悪化させることなくパルス計測を可能とすることができる。特に、本実施の形態においてｎが奇数の場合、両極性ノイズの符号を維持して平均化することができるため、両極性信号と単極性の信号を良好なＳＮ比（信号雑音比）をもって識別することができる。
【００３８】
以上により、パルス計測において、出力パルスの識別レベルを設定する際に、回路ノイズまたはα線ノイズを除去する識別レベルを低く設定することができ、その分、センサ出力パルスが小さい場合でも計測感度を落とすことなくパルス計測を行うことができる。
【００３９】
＜第２の実施の形態＞
次に、本発明の第２の実施の形態を図４に基づいて説明する。
【００４０】
図４の放射線測定装置のｎ乗パルス識別装置４は、積分識別装置６、差分識別装置７、および波高・パワー識別装置８を備えている。Ａ／Ｄ変換器３の出力パルスは積分識別装置６および差分識別装置７に入力され、両者の出力は波高・パワー識別装置８に入力される。波高・パワー識別装置８の出力信号はパルス計数装置５に入力される。ここで、積分識別装置６は、第１の実施の形態で説明したパルス識別方式によりパルスを識別する。
【００４１】
以下、本実施の形態における差分識別装置７が行う処理の第１実施例について説明する。この実施例では次の原理により識別を行う。
【００４２】
Ａ／Ｄ変換器３の出力の最大値Top(k)およびボトム値Bottom(k)をそれぞれ、
Top(k) = +b×S(k-2) + c×S(k-1)
Bottom(k) = a×S(k-3) + d×S(k)
とすると、波高値High(k)は、先行技術での式（１）を簡略化し、
High(k) ＝ ＋Top(k) - Bottom(k)
とすることができる。
【００４３】
本実施の形態では、まず最大値Top(k)の２乗値とボトム値Bottom(k)の２乗値の差を求める。すなわち、
ここで、通常の回路ノイズのみに信号パルスが重畳したケース（以下、ケース１とする）では、
Top(k) >> +Bottom(k)
であるため、
X = High(k) × Top(k) …（３）
一方、非常に大きなサージノイズに信号パルスが重畳した極端なケース（以下、ケース２とする）では、近似値として、
Top(k) = Bottom(k)
であるため、式（２）から、
X = High(k) × (2×Top(k)) …（４）
となる。式（３）および式（４）から、
X／Top(k) = α × High(k) …（５）
（ここで、ケース１ではα＝１、ケース２ではα＝２）
となり、X／Top(k)はほぼパルス波高High(k)の一次関数となる。つまり、サージノイズ上にＳＲＮＭセンサの出力パルスが重畳した場合でも、この値によって識別することにより、数μｓの周期のサージノイズ上に重なった数百ｎｓのＳＲＮＭセンサ出力パルスを識別し、計数することが可能となる。
【００４４】
このように本第１実施例によれば、ｎ乗値の差を用いることによって、仮にセンサパルスよりもパルス幅の広い外来ノイズが誘導されても、その影響を低減することができる。
【００４５】
次に、本実施の形態における差分識別装置７の処理の第２実施例について説明する。式（１）を、
D1(k) = c×S(k-1) - d×S(k) …（６）
D2(k) = a×S(k-3) - b×S(k-2) …（７）
とすると、波高値High(k)は、
High(k) = +D1(k) + D2(k) …（８）
となる。式（８）の右辺各項のｎ乗の和Yは、
Y = D1(k)ｎ + D2(k)ｎ
となり、これは近似値として、
Y = High(k)ｎ
すなわち、
Y−ｎ = High(k) …（９）
となる。よって、この場合は、近似値としてY−ｎを指標に用いることにより、パルス識別が可能である。なお、ｎ＝１のときは、式（１）に相当する。
【００４６】
以上、差分識別装置７により、サージノイズ上にＳＲＮＭセンサ１の出力パルスが重畳した場合でも、差分によって求めた指標を用いることにより、数μｓの周期のサージノイズ上に重なった数百ｎｓのＳＲＮＭセンサ１のパルスを識別し、計測することが可能となる。
【００４７】
図５に積分識別装置６および差分識別装置７の識別性能の差異をまとめて示す。積分識別装置６は、ホワイトノイズなど両極性ノイズで、センサ出力のパルス幅より短い周期のノイズを低減するのに有効であり、一般に回路ノイズ、センサのα線ノイズなど、センサ出力より信号（波高値）が小さいノイズをさらに低減するのに有効である。一方、差分識別装置７は、センサ１の出力パルス幅より長い周期のノイズに対して有効であり、一般に外部からの数μｓパルス幅の誘導ノイズを除去するのに有効である。したがって、両者のロジックを最適に調整し、両者がともに成立した場合のみパルスを計数することにより、センサ出力のみをより正確に計数することができる。
【００４８】
以上のように、本第２実施例によれば、波高値に相当するサンプル値の差のｎ乗値を用いることにより、第１実施例と同様に、センサパルスよりもパルス幅の広い外来ノイズが誘導されても、その影響を低減し、より高精度の放射線測定を行うことができる。
【００４９】
次に、本実施の形態の第３実施例として、波高・パワー識別装置８について説明する。波高・パワー識別装置８は、積分識別装置６からパルスの積分値を、差分識別装置７からパルスの波高値相当の値を受け取る。この両者の比率（積分値／波高値）は、センサ出力パルスの場合、ほぼパルス幅相当となり一定の値を示す。一方、高周波成分を含むホワイトノイズは、波高値がセンサ出力と同等の場合でも積分値が小さいため、比率は小さくなる。また、パルス幅の広いサージノイズは、積分値は大きく波高値は小さいため、その比率はセンサ出力より大きくなる。したがって、この比率が一定範囲内である場合に、パルスを計数することにより、これらのノイズの影響を低減することができる。
【００５０】
以上により、パルス計測において、サージ状の外来ノイズが誘導された場合でも、その周期が数μｓとセンサ出力パルス幅１００ｎｓより長い場合、そのサージ状ノイズの影響を除去し、そのノイズに重畳したパルスも計数することができる。かくして、本第３実施例によれば、本実施の形態の第１実施例または第２実施例で述べた差分によるパルス計数と、第１の実施の形態で述べたｎ乗値によるパルス計数を併用し、その比率を用いることにより両者の特性を同時に実現し、より正確なパルス計測が可能となる。
【００５１】
＜第３の実施の形態＞
次に、本発明の第３の実施の形態を図６に基づいて説明する。
【００５２】
図６に示す放射線測定装置は、図１に示した放射線測定装置において、Ａ／Ｄ変換器３とｎ乗パルス識別装置４との間に、両極性の信号をパルスの主成分である極性のみの信号に変換する単極性変換手段９を介挿したことを特徴とするものである。なお、ＳＲＮＭセンサ１はパルス出力が得られる核***検出器であるが、これは他のパルス出力が得られる電離箱等の放射線検出器であってもよい。
【００５３】
このように構成された本実施の形態における単極性変換手段９の機能について図７を参照して説明する。図７（ａ）は、図２（ａ）に示されるような検出器出力を２次微分処理した場合のパルス波形の一例を示す。２次微分処理は、プリアンプ２Ａにおいてアナログ回路で処理する場合と、図２の波形をＡ／Ｄ変換器３でサンプリングし、ディジタル処理で２次微分演算を行う場合が想定されるが、今回は両者とも想定している。この図７（ａ）をｎ乗パルス識別装置４において偶数ｎを用いてｎ乗演算後、平均処理した結果を図７（ｂ）に示す。両極性が単極性に変化するために、パルス幅が広くなる。
【００５４】
そこで、単極性変換手段９は、図７（ａ）のパルスの主成分は負極性であるため（図７（ａ）の下方が負極性）、正極性の成分をゼロまたはそれに近い数値に置き換え、図７（ｃ）の波形に変換する。この波形を、ｎ乗パルス識別装置４に内蔵されるｎ乗演算手段によって、偶数ｎを用いてｎ乗演算した後、たとえば２乗演算した後、ｎ乗パルス識別装置４に内蔵される平滑手段によって平均化処理した結果を図７（ｄ）に示す。これから分かるように、パルス幅の広がりを、単極性変換手段９を用いない図７（ｂ）に比べ狭くすることができる。
【００５５】
かくして本実施の形態によれば、パルス計測において、偶数ｎを用いてｎ乗演算を行ってもパルス幅が広がることなく、その結果、パルスの重なりによって１個のパルスを計数できなくなる、パルスパイルアップによる誤計数を低減することができる。
【００５６】
＜第４の実施の形態＞
次に、本発明の第４の実施の形態を図８に基づいて説明する。
【００５７】
ＳＲＮＭセンサ１の出力は、帯域制限機能を有するプリアンプ２Ｂで増幅・帯域制限され、Ａ／Ｄ変換器３に入力される。Ａ／Ｄ変換器３の出力は、バンドパスフィルタ（ＢＰＦ）１０に入力され、特定の周波数帯域に制限される。通常、バンドパスフィルタ１０の周波数帯域はＳＲＮＭセンサ１の出力特性によって決定されるが、たとえば１００ｋＨｚから４００ｋＨｚまでの範囲の周波数成分のみを通過するディジタルフィルタとして構成することができる。この測定帯域は複数個設定してもよいが、本実施の形態では一つの帯域に特定した単純な場合の構成で代表させる。
【００５８】
ディジタルフィルタ処理を行う際、入力のサンプル値は、バンドパスフィルタ１０の出力周波数帯域に対して適切なサンプリング間隔およびビット数に調整される。このバンドパスフィルタ１０の出力はｎ乗値演算装置１１に入力され、ｎ次モーメント値に変換される。
【００５９】
本実施の形態では、ｎ乗値演算装置１１の出力側に、第１平滑装置１２Ａ、中間値の平均を取るデータ除去平均化装置（ＤＥＡ）１３、および第２平滑装置１２Ｂが配設され、その平滑個数とデータ除去個数を外来ノイズのサージノイズ幅と到来周期およびＭＳＶ計測に必要なデータ除去可能率を元に調整することによって、外来ノイズの影響を受けないＭＳＶ計測が可能となる。これによってさらに信頼性の高い計測を実現することができる。
【００６０】
さらに、データ除去平均化装置１３において、最大値のみを平均化することにより、低計数率時にもＭＳＶ指示値を得ることができ、プリアンプ２Ｂの回路ノイズ等で決まっているＭＳＶ計測の測定下限を拡大することができる。
【００６１】
なお、本実施の形態においては、ｎを奇数に設定した場合、サージ状のノイズは正負両極性のノイズであるため、平均化処理することによって正負の値が相殺され、単極性のセンサ出力と同じ極性の値を出力として選択することにより、サージ状ノイズの影響を低減することができる。
【００６２】
以下に述べる実施の形態では、ｎ次モーメント演算として最も簡単な２乗値に代表させて説明する。つまり、バンドパスフィルタ９のｋ番めの出力サンプル値をS(k)とし、ｎ乗値演算装置において、
Out1(k) = S(k) × S(k) …（１０）
の演算が行われるものとして説明する。
【００６３】
この出力Out1(k)は、第１平滑装置１２Ａに入力され、ある個数分について平均化される。この際、バンドパスフィルタ１０の出力の平均値がオフセットを有する場合は、サンプル値S(k)の平均値も同時に算出し、その２乗値を算出し、２乗演算結果Out1(k)から引き算する。つまり、ｎ個平均化する場合は、
Out2(k2) = (ΣOut1(k))/n − {(Σs(k))/n }２ …（１１）
を演算する。ただし、Σはｎ個のサンプル値を加算することを示すものとする。
【００６４】
この際、加算するｎの個数は、想定される外来ノイズのパルス幅に相当する個数分に設定され、また、この個数はディジタル演算のしやすいように２のｎ乗個に設定される。つまり、図９に示すようなノイズを想定した場合、すなわちパルス幅２０μｓのサージ状パルスが２ｍｓ周期で到来するパルス状ノイズを想定した場合、１μｓ周期でバンドパスフィルタ９の出力が得られるものとして、２０μｓ以上のデータを平均化する。ただし、ディジタル演算で行う場合、ビットシフトによる割り算が行えるように２のｎ乗個に加算数を設定した方が演算に有利であるため、このケースでは３２（＝２５）個加算する。出力間隔ｋ２は、バンドパスフィルタ９の出力間隔１μｓに対して、その３２倍である３２μｓ周期となる。
【００６５】
この第１平滑装置１２Ａの出力は、データ除去平均化装置１３に入力される。データ除去平均化装置１３では、３２μｓ周期の第１平滑装置１２Ａの出力を、特定個数ごとに区切って、その大小関係を比較し、特定のデータを除く処理を行う。図９に示すようにサージノイズの除去を想定した場合は、以下の処理を行う。図９に示すサージ状パルスの到来間隔は、約２ｍｓである。よって、２ｍｓ以下の間隔でデータを区切り、その中でサージ状パルスをサンプリングして得たパルスデータのみを除けば、残りはノイズのないデータのパワー値に換算することができる。つまり、この区間の大小関係を比較し、値の大きいものから２個以上を除けば、このサージ状ノイズの影響を除去することができる。ここで最低の除去数を２個に設定したのは、タイミングによっては２０μｓ幅のサージノイズが２個分の第１平滑装置１２Ａの出力信号に混入する可能性があるためである。また、大きい値を除いた場合、平均値を保持するために、小さい値も除いた最大の個数と同数、つまり、２個分の小さい値のデータを除く必要がある。このように大きい値と小さい値を２個ずつ除いた残りのデータを平均化する。図９のノイズの場合は、３６個分の第１平滑装置１２Ａの出力をひと区切りとすると、約１．１６ｍｓの平滑区間となり、その中に含まれるサージパルスは最大２個であるため、大きいもの２個と小さいもの２個を除き、残りの３２個を平均化する。この場合も、除いた後の個数が２のｎ乗個となるように、加算数を決定する。
【００６６】
また、サージ状パルスの到来間隔が短くなったケースでは、データ除去平均化装置１３での加算数を少なくする場合と、加算平均区間内でのサージ状パルスの到来個数を評価し、その個数の２倍分の最大・最小データを除く処理を行えば、ノイズを除去することができる。
【００６７】
ただし、このようにデータを除く場合に問題となるのは、データの除去率である。つまり、ＳＲＮＭセンサ１の出力パルス数が少ない場合、データを除去することによって、平均したパワーつまりＭＳＶ値が実際の値より低めに表示される。そのパルス数が充分な場合、そのランダム性からある程度のデータ除去は誤差範囲内におさまる。したがって、その測定精度から得られる除去限度率を確保するようにデータ除去を行えば、測定に悪影響を与えることなしに、間欠的に誘導されるサージノイズによる突発的な指示値の上昇を、あるしきい値の範囲内の値に制限することができる。
【００６８】
データ除去平均化装置１３の他の機能として、最大値のみを選択することにより、ＳＲＮＭセンサ１のパルス発生率が低い場合にも、ＭＳＶ指示値を得ることが可能となる。つまり、低いパルス発生率の場合、ひとつの平滑化区間においてセンサパルスが１個も到来しない時間が多く存在する。したがって、この到来しない区間を平滑化に含めないようにすることにより、パルス発生率が低い場合もＭＳＶ指示値を得ることが可能となり、ＭＳＶ計測の下限を延長させることができる。これは、パルス計測とＭＳＶ計測を併用した計測手法であるということができる。ただし、この場合、パルス発生率に対するＭＳＶ指示値の直線性を補正するために、データの除去割合をパラメータとした補正関数を用いて補正する必要がある。
【００６９】
次に、データ除去平均化装置１３内でデータを除去する処理の基準に関する一実施例を以下説明する。模擬中性子パルスの入力に対し第１平滑装置１２Ａにおいて３２μｓで平滑化した場合のシミュレーション結果の一例を図１０に示す。図１０の実線は図の縦左軸に対応し、平滑後のＭＳＶ指示値の平均値Ｓを示している。これから分かるように、パルス計数率がＭＳＶ計測範囲内すなわち１×１０5ｃｐｓ以上においては、概ね、パルス計数率と平滑後のＭＳＶ指示値すなわち２乗電圧とは、ほぼ比例して推移する。他方、パルス計数率が１×１０5ｃｐｓ〜２×１０5ｃｐｓの範囲では回路ノイズの影響で比例関係にはなっていない。
【００７０】
また、３２μｓでのＭＳＶ指示値の平均値Ｓと標準偏差σに対して、
Ｘ＝（Ｓ＋６σ）／Ｓ
で表される指標をパルス発生率ごとに評価したものを、図１０の縦右軸に対応する破線で示す。ＸはＭＳＶ計測下限で最大５．３程度となり、同様の評価で１×１０6ｃｐｓ以上では２〜３程度で推移する。
【００７１】
一方、他の評価方法として、バンドパスフィルタ１０で１００〜４００ｋＨｚの帯域制限処理を施して第１平滑装置１２Ａで模擬中性子パルスを３２μｓで平滑化したときのＭＳＶ計測値の揺らぎ率Ｙ、すなわち標準偏差σと平均値Ｓとの比率：Ｙ＝σ／Ｓ の推移を理論式で評価した場合のグラフを図１１に示す。図から分かるように、ＭＳＶ計測範囲（計測下限１×１０5ｃｐｓ）におけるＹの値は０．４以下であり、ＭＳＶ計測範囲の下限で揺らぎは最大となる。
【００７２】
これらの図から分かるように、波形の揺らぎは、特にＭＳＶ計測範囲の下限である約１×１０5ｃｐｓ程度において大きい。したがって、特にこの揺らぎの大きい計数率：１×１０5ｃｐｓ程度における、上述したＸあるいはＹにより例示される揺らぎの度合いを示す指標を、第１平滑装置１２Ａの積分時間の条件で事前に評価し、得られた評価値を超えたものはノイズと判定することとする。
【００７３】
すなわち、第１平滑装置１２Ａでパルスを３２μｓで平滑化した場合、第１平滑装置１２Ａで以前に平滑処理された結果や、あるいは予め模擬中性子パルスを用いて評価した結果を用いて得られる評価指標および評価基準と、平滑化した結果との比較を行い、この評価指標に基づいてパルスを評価したときに、評価基準を超えるような場合、例えば評価指標が前回の平滑処理結果の８倍と規定されるしきい値を超えた場合、この部分を異常値と判定する。
【００７４】
図１２は、ノイズを含むパルス波形を模式的に示したものである。通常のＭＳＶ値は破線で囲まれる範囲を揺らいでいる。この範囲を超えるものはノイズである。このしきい値は、その平均のＭＳＶ指示値によって変化するが、もっとも変動の大きいＭＳＶの計測下限での評価値を用いれば、全ＭＳＶ計測範囲において保守的な評価として適用可能である。
【００７５】
また、ここでのデータ除去平均化装置１３においては、ノイズと判定された部分のデータについては、このデータ自体を除去する場合と、このデータを上述したしきい値の範囲内の値に置き換える場合とがある。前者の場合、ノイズをほぼ完全に除去することができるが、上述したデータ除去率の範囲で除去する必要があるためデータ除去率を評価する必要となる。逆に後者の場合、完全なノイズ除去はできないが、データ除去率の評価が不要であるという特徴がある。
【００７６】
以下、後者の場合、すなわちノイズ部分の値を置換する場合の一例を詳述する。図１３は、ＭＳＶ指示値の揺らぎの範囲と平均値の変化幅の関係を示すグラフである。図中の実線はパルスを32μｓで平滑化したときのＭＳＶ信号出力、破線は中性子束の最大変化を示している。図中一点破線で示された、上述の方法により評価されるＭＳＶ計測の揺らぎの変化幅は、破線で示された、本来計測すべき中性子束の変化率に比べて十分大きい。よって、この評価されたＭＳＶ計測の揺らぎの変化幅を超えた場合を、異常値と判定することとする。
【００７７】
この揺らぎの変化幅すなわち最大揺らぎを超えるノイズが誘導された場合のパルスのＭＳＶ指示値の推移の一例を図１４（ａ）に示す。図中矢印Ａは、ＭＳＶ計測の揺らぎの変化幅、例えば前回の平滑処理結果の８倍と規定されたしきい値を示す。このＡを超えている部分は異常値と判定され、この異常値を、正常値として、本来の中性子束の最大変化率を前回のサンプリング値に乗じた値に置き換えることとする。置換処理の結果を図１４（ｂ）に示す。ここで、３２μｓ幅での想定される中性子束の最大変化率は、例えば1.03倍程度であり、上述したＭＳＶ計測の揺らぎ範囲よりも十分小さいものである。よってここでの置換処理は、異常値を前回値の1.03倍に置き換えるものとしている。このように、異常値が検出された場合の変化率も、平滑区間内での監視すべき中性子束の最大変化率をあらかじめ評価しておくことにより、計測の時間応答性を損なうことなく、異常値のみを除去することができる。
【００７８】
次に、データ除去平均化装置１３の出力は、第２平滑装置１２Ｂに入力され、測定値の揺らぎが必要とされる測定精度を満たし、また、応答要求を満たす範囲となるように、平均化処理される。この結果は、ＭＳＶ中性子評価装置１４に入力され、測定されたＭＳＶ値が中性子束の値に換算される。
【００７９】
また、本実施の形態では、ノイズ特性評価装置１５を設けるのが好適である。ここでは、ノイズ波形の特性であるサージ状パルスのパルス幅と到来周期の最小値を評価し、ＭＳＶ計測の測定範囲において必要となる最大データ除去率の範囲内において、第１平滑装置１２Ａで平均化のために用いるデータ数、データ除去平均化装置１３のデータ平均区間と、データ除去数、および、第２平滑装置１２Ｂの平滑化フィルタの時定数を設定する。
【００８０】
このような構成により、ＭＳＶ計測において、サージ状の外来ノイズが誘導された場合でも、そのサージ状のパルス幅と、その到来周期を評価し、ＭＳＶ計測の許容データ除去率を満たす範囲内でサージノイズデータを除去することにより、間欠的なサージノイズを完全に除去することができる。
【００８１】
＜参考例＞
次に、本発明の参考例を図１５に基づいて説明する。
【００８２】
図１５に示す放射線測定装置では、放射線センサとして、常温半導体であるＣｄＴｅを用いたＣｄＴｅセンサ１６を用いている。放射線センサとしては、エネルギー情報の取得が可能な、ＮａＩ等のシンチレーションセンサと光電子増倍管を組み合わせたもの、半導体センサであるＧｅセンサなども同様に用いることが可能である。ＣｄＴｅセンサ１６の出力は、チャージアンプ（ＣＡ）１７に入力される。チャージアンプ１７は、入力パルスの電荷を積分し、その電荷量に応じたパルス波高を有するパルスに変換し出力する。なお、チャージアンプ１７はＣｄＴｅセンサ１６に対しその動作電圧を供給する。チャージアンプ１７の出力は、一般に放射線のエネルギー測定に用いられるパイルアプリジェクション回路、ポールゼロキャンセル回路などで波形整形処理が施され、ＭＳＶ計測装置１８、電流計測装置（ＣＤ）１９、およびパルスカウンタ（ＰＣ）２０に入力される。ＭＳＶ計測装置１８では、周波数帯域の制限を行った後、ｎ乗の平均化演算が行われ、ＭＳＶ計測値つまり２次モーメント値に変換される。電流計測装置１９では、１次モーメント値である平均電流値が測定される。また、パルス計測装置２０では、パルス数の計数が行われる。これらＭＳＶ計測値、電流計測値、およびパルス計数値は、それぞれエネルギー評価装置２１に入力され、ＭＳＶ値とパルス計数の比率、または、ＭＳＶ値と直流電流値の比率（１次と２次のモーメント値の比率）をもとに平均放射線エネルギーを評価する。この平均のエネルギー値と、前記の測定値は、線量評価装置２２に入力され、照射線量、または、物質での吸収線量、または、人体へのリスク率を含めた線量当量に換算される。
【００８３】
エネルギー評価装置２１における放射線の平均エネルギーの推定方法を以下に説明する。
【００８４】
チャージアンプ１７の出力は、ＣｄＴｅセンサ１６内で放射線が吸収されたエネルギーに比例した波高値を有するパルスである。したがって、その反応の起こる確率をＮとし、吸収エネルギーをｑとすると、ＭＳＶ値、パルス計数値、および電流値は、以下の式で近似することができる。
【００８５】
ＭＳＶ値 ： k1×q２×N
ｎ次モーメント値 ： kn×qｎ×N
パルス計数値 ： k2×N
直流電流値（１次モーメント値） ： k0×q×N
ただし、k0、k1、k2、knはそれぞれ補正係数である。
【００８６】
よって、これらの比率は、
ＭＳＶ値／パルス計数値＝(k1／k2)×q２（一般に、kn×qｎ）
ＭＳＶ値／直流電流値＝(k1／k0)×q
ｎ次／ｎ’次モーメント値 ＝ kn×q（ｎ−ｎ’）／kn’
となり、これらの補正係数k0，k1，k2，kn等を予め評価し、これら測定値の比率を用いることにより、結晶中での吸収エネルギーの推定が可能となる。
【００８７】
図１６は、ＣｄＴｅセンサ１６によって種々の核種からの放射線、つまり、エネルギーの異なる放射線、を測定した場合の、パルス計数値およびＭＳＶ値（縦軸）を、市販の放射線サーベイメータで測定した線量（横軸）に対してプロットしたものである。一般に、サーベイメータ等は、放射線のエネルギーに対する感度特性を、放射線のエネルギーに対する線量当量の評価曲線に一致するように遮蔽、または、内部の補正係数を調整している。つまり、エネルギーの低い放射線に対しては、１個の放射線が入射した場合の結晶中での吸収エネルギーが小さいためパルス計数率が大きくなるが、１個の放射線の線量は低い値を示す。逆に、エネルギーの高い放射線に対しては、パルス計数率は低くなるが、１個の放射線で発生する電荷量が大きいため線量は大きくなる。したがって、測定しているパルス計数値または電流値の感度を線量当量への寄与率と同じになるように遮蔽等で調整する。図１６は、これらの感度補正を行っていないため、低いエネルギーの放射線に対してはパルス計数値が高くなっており、線量に対してパルス計数値およびＭＳＶ計測値はランダムになっている。
【００８８】
しかし、これをＭＳＶ値／パルス計数値の線量に対する特性としてプロットすると、図１７に示すように線量当量に対して単調な特性となる。したがって、この特性を事前に評価しておくことによって、両者の比率から線量当量への換算が可能である。また、同様に、このＭＳＶ値／パルス計数値の比率は、入射エネルギーに対しても単調な特性となり、この特性も評価しておくことにより、平均の入射放射線エネルギーを推定することができる。この場合、放射線のエネルギーに対する人体の各部位の吸収特性を用いることにより、各部位で吸収線量をより正確に評価することが可能となる。
【００８９】
さらに、パルス計測で線量評価を行う場合、パルスの波高分布を測定し、入射エネルギーのエネルギー情報を取得し、線量に換算する方法と、前記のセンサ自体の構造を工夫し、パルス計測の感度を線量応答特性と同一にする場合がある。また、電流値で線量を評価する方法としては、後者のセンサ構造を工夫し、感度応答を調整する方法がある。したがって、これら周知の手法と本参考例の線量評価方法を併用することにより、さらに精度のよい線量評価を行うことができる。つまり、たとえば、センサの感度特性を単独にある程度調整し、その後に本発明の補正関数を評価して用いることにより、さらに正確な線量評価が可能となる。
【００９０】
さらに、パルスがパイルアップし、パルス計数の数え落としが生じるような高い計数率では、波高情報を取得し線量に換算する前者の方法では正確な線量は評価できない。また、遮蔽を設けた後者の場合も、数え落とし分を補正する必要がある。しかし、本参考例のようにパルス計測とＭＳＶ計測を同時に実施することにより、パルス計測がパイルアップによって飽和した場合でも、ＭＳＶ計測によって線量の評価が可能となり、広いレンジの測定が可能となる。ただし、この場合の平均エネルギーの推定は、パルス計測のパイルアップの効果を補正する必要があるが、ある程度センサ自体の感度を線量応答に近似させておくことにより、その誤差は無視できる範囲に抑えることができる。
【００９１】
また、ＭＳＶ値／パルス計数値の比率と同様に、電流値／パルス計数値の比率と、上述した本参考例によるＭＳＶ値／電流値の比率を用いた場合も、同様の補正関数を求めておくことにより平均放射線エネルギーの推定が可能となる。
【００９２】
かくして、本参考例によれば、線量評価において、ｎ次モーメント値とパルス計測を併用することにより、平均入射エネルギーの推定から、より正確な線量評価を行うことができる。また、パルスがパイルアップし、計数誤差が生じる条件においても、ｎ次モーメントどうしの比率を用いることにより、同様に平均エネルギーの推定が可能となり、従来よりも広い測定レンジで正確な線量評価を行うことができる。
【００９３】
さらに、本参考例の変形例を以下説明する。ここでは、上述したＭＳＶ値の代わりに、１乗値（平均電流）、２乗値、３乗値、…、ｎ乗値をそれぞれ算出し、それぞれの補正関数を算出しておき、この補正関数の逆マトリックスを解くことにより、エネルギー分布の推定が可能となる。
【００９４】
つまり、各ｎ次のモーメントの計測値は、
x1 = a1[1:n]×E[n:1] (電流計測に相当)
x2 = a2[1:n]×E[n:1] (ＭＳＶ計測に相当)
x3 = a3[1:n]×E[n:1]
・・・
xn = an[1:n]×E[n:1]
xn: ｎ次モーメント値［スカラー量］
an: 応答マトリックス[１行ｎ列行列]
E: エネルギー分布[ｎ行１列行列]
となるから、行列の乗算として、
X[n:1] = A[n:n] × E[n:1]
と表すことができる。この行列Ａの逆行列Ａ−１を用いて、放射線のエネルギー分布は、
E[n:1] = A-1[n:n] × X[n:1]
として求めることができる。ただし、この１乗からｎ乗までのモーメント計測は、必要とするエネルギーバンド幅に相当する個数分を前記の式から選択すればよく、１次モーメントである平均電流値を除くことにより、すべて交流の計測手段のみで構成することができる。
【００９５】
以上のように、ＭＳＶ計測とパルス計測または電流計測を組み合せ、その比率から平均放射線エネルギーを推定し、線量に換算することができる。これは、従来のセンサ構造を用いて感度を補正していた手法と併用することにより、より線量応答に近い特性を容易に実現することができる。また、従来のパルス波高を求めて線量を算出していた手法に比べ、波高を選別する必要がなく、簡単な構成で実現が可能である。さらに、複数のｎ次モーメント値を用い、その応答関数から放射線のエネルギーを再構成することができ、パルス計測が困難な高計数率においても放射線のエネルギー分布の測定が可能となり、その情報からより正確な線量を評価することができる。かくして、これらを単独に、または従来の線量評価手法と併用することにより、より広いレンジの線量を、より正確に、一括して監視することの可能な放射線測定装置を提供することができる。
【００９６】
【発明の効果】
以上説明したように、本発明の放射線測定装置によれば、パルス波形のｎ乗値を算出し、その値でパルス識別を行うことによって、信号レベルの小さい両極性の回路ノイズおよびセンサのα線ノイズを低減し、センサパルスに対するこれらノイズの割合を低減することができる。これにより、従来は回路ノイズに紛れていたセンサ信号を測定することが可能となる。
【図面の簡単な説明】
【図１】 本発明による放射線測定装置の第１の実施の形態を示すブロック図である。
【図２】 （ａ）はＳＲＮＭセンサ出力パルスとサンプリングの例を示す説明図、（ｂ）はＳＲＮＭセンサ出力の２乗値を示す説明図である。
【図３】 （ａ）〜（ｃ）はｎ乗パルス計測手段の作用を説明するためのグラフである。
【図４】 本発明による放射線測定装置の第２の実施の形態を示すブロック図である。
【図５】 積分識別装置と差分識別装置の特徴を比較する図表である。
【図６】 本発明による放射線測定装置の第３の実施の形態を示すブロック図である。
【図７】 （ａ）〜（ｄ）は単極性変換手段の機能を説明するためのグラフである。
【図８】 本発明による放射線測定装置の第４の実施の形態を示すブロック図である。
【図９】 ノイズ波形例を示す波形図である。
【図１０】 第４の実施の形態において模擬中性子パルスを入力したときの第１平滑装置のシミュレーション結果と指標Ｘの推移を示すグラフである。
【図１１】 第４の実施の形態において第１平滑装置の出力の揺らぎ率Ｙを評価した結果を示すグラフである。
【図１２】 ＭＳＶ揺らぎとノイズの関係を説明するための波形図である。
【図１３】 ＭＳＶ指示値の揺らぎの範囲と平均値の変化幅の関係を示すグラフである。
【図１４】 （ａ）は最大揺らぎを超えるノイズが誘導された場合のパルスのＭＳＶ指示値の推移の一例を示すグラフ、（ｂ）は（ａ）を第４の実施の形態のデータ除去平均化手段によって処理した結果を示すグラフである。
【図１５】 本発明の参考例を示すブロック図である。
【図１６】 線量に対するＣｄＴｅセンサのパルス計数およびＭＳＶ指示値を示すグラフである。
【図１７】 線量に対するＣｄＴｅセンサのＭＳＶ指示値とパルス計数値の比率を示すグラフである。
【図１８】 従来のディジタル式原子炉起動監視装置のブロック図である。
【符号の説明】
１ ＳＲＮＭセンサ
２ アナログ式増幅器
２Ａ プリアンプ
２Ｂ プリアンプ
３ Ａ／Ｄ変換器
４ ｎ乗パルス識別装置
５ パルス計数装置（ＰＣ）
６ 積分識別装置
７ 差分識別装置
８ 波高・パワー識別装置
９ 単極性変換手段
１０ バンドパスフィルタ（ＢＰＦ）
１１ ｎ乗演算装置
１２Ａ 第１平滑装置
１２Ｂ 第２平滑装置
１３ データ除去平均化装置（ＤＥＡ）
１４ ＭＳＶ中性子評価装置
１５ ノイズ特性評価装置
１６ ＣｄＴｅセンサ
１７ チャージアンプ
１８ ＭＳＶ計測装置
１９ 電流計測装置
２０ パルス計数装置
２１ エネルギー評価装置
２２ 線量評価装置
２３ パルス計数手段
２４ 和演算手段
２５ パワー演算手段
２６ 加算平均手段
２７ 原子炉出力評価手段
【発明の属する技術分野】
本発明は、放射線測定装置に関する。より詳細には、本発明は、放射線測定装置により放射線測定を行う場合に広レンジの監視に用いられるパルス計測およびキャンベル計測に関し、そのディジタル信号処理による耐ノイズ性の改善に関する。
【０００２】
【従来の技術】
広いレンジの放射線を測定するためにパルス計測およびキャンベル計測が併用されている。パルス計測は、一般に放射線センサからのパルス数をカウントするものであり、そのパルスが重なりあって計数できなくなった場合に、キャンベル計測が行われる。たとえば、原子力の出力監視では、原子炉圧力容器内部に６〜１０本の起動領域センサ（ＳＲＮＭセンサ）および１００〜２００本の局所出力領域センサ（ＬＰＲＭセンサ）を設置し、それらの出力をそれぞれ起動領域モニタおよび出力領域モニタによって測定し、原子炉の出力を約１１桁の監視幅で監視している。
【０００３】
このうち起動領域モニタは、原子炉出力の低い領域、つまり原子炉出力が定格出力の１０−９％から１０−４％まではＳＲＮＭセンサの出力パルスの個数を計数すること（以後、パルス計測と呼ぶ）により出力を監視する。一方、原子炉出力が高い領域、つまり原子炉出力が１０−５％〜１０％では、センサの出力パルスの重なりにより生じる揺らぎのパワーを測定すること、すなわちキャンベル計測により原子炉の出力を監視する。
【０００４】
以下、原子炉起動監視装置におけるパルス計測およびキャンベル計測の先行技術例（特開２０００−１６２３６６号公報参照）について図１８を参照して説明する。
【０００５】
図１８に示した原子炉監視装置は、原子炉内の中性子に感応して中性子量に対応する電気パルスを出力するＳＲＮＭセンサ１と、その出力パルス信号を増幅し整形するアナログ式プリアンプ２と、ＳＲＮＭセンサ１の出力パルス幅より狭い時間間隔でサンプリングしディジタル値に変換するＡ／Ｄ変換器３と、このサンプリングされたサンプルデータからパルス数を計数し、原子炉の低い領域での出力レベル値に変換するパルス計数手段（ＰＣ）２３と、測定精度を向上させるためにＡ／Ｄ変換器３のサンプル値を積算する積算手段２４と、その積算値を２乗してパワーを演算するパワー演算手段２５と、そのパワーを平均化する加算平均手段２６と、パルス計数手段２３の計数結果および加算平均手段２６の演算結果に基づいて原子炉起動時の出力を連続監視する原子炉出力評価手段２７とから構成されている。
【０００６】
このように構成されたディジタル式原子炉起動監視装置においては、ＳＲＮＭセンサ１の出力パルスの形状をプリアンプ２で増幅整形し、増幅整形されたパルスをＡ／Ｄ変換器３が高速にサンプリングし、そのサンプル値に対してある一つの、または、複数のロジック演算を行い、その演算結果がそれぞれ予め設定した範囲内にある場合にセンサ出力パルスとしてパルス計数手段２３で計数する。
【０００７】
一方、同じサンプル値を積算手段２４で積算処理し、キャンベル計測に必要なサンプリングレートまで落とすとともに、ダイナミックレンジを稼ぐために積算し、等価ビット数の改善を計る。その結果をパワー演算手段２５でバンドパスフィルタ処理を行ったのち、その２乗和を演算し、その演算結果を加算平均手段２６で平均し、キャンベル出力値を算出する。これらパルス計数値とキャンベル出力値は原子炉出力評価手段２７で評価され、原子炉出力として表示される。
【０００８】
この構成では、パルス計数手段２３によりパルスの波高のみではなく、パルス幅の情報も取り入れた識別により広いパルス幅のノイズを除去し、センサ出力パルスのみを計数することができる。つまり、図１８の原子炉起動監視装置では、たとえば、１００ｎｓパルス幅のＳＲＮＭセンサ１の出力を２５ｎｓ間隔でサンプリングする。このサンプリングによるサンプルデータのうち、パルス幅に相当するｋ−３番めからｋ番めまでの４つのサンプル値[S(k), S(k-1), S(k-2), S(k-3)]を用い、パルス立ち上がり部相当のサンプル値S(k-3)、立ち下がり部相当のサンプル値S(k)と、その間の２つのサンプル値S(k-2)およびS(k-1)を用いて下記のような演算を行い、その演算結果(Out(k))をパルス識別の指標とし、それが予め設定したレベルの範囲内にある場合、中性子パルスとして計数する。
【０００９】
【００１０】
なお、このようなセンサパルス形状に合致したケースを検出できる指標を複数設けて、それらの論理積（ＡＮＤ）をとることによって、この識別性能はさらに向上する。これにより、最もＳＲＮＭに混入しやすいとされる数μｓ幅のサージノイズが混入しても、それをほぼ完全に除去し、１００ｎｓ幅のセンサパルスのみを計数可能とすることができる。
【００１１】
一方、キャンベル計測に関しては、パワー演算手段２５により周波数帯域を制限し、２乗平均値を演算する。この構成では、ソフトウエアにより周波数帯域を設定できるため、ノイズが測定帯域と同等の周波数となった場合、測定帯域をソフトウエア上で変更することによりノイズの誘導を低減することができる。
【００１２】
【発明が解決しようとする課題】
しかし、上述の先行技術による原子炉起動監視装置においては、以下に述べるような課題が存在している。
【００１３】
第１の課題は、両極性ノイズの低減の問題である。つまり、数μｓ幅のサージノイズとセンサパルスとが重なった場合、不感時間を設けずにその重なったセンサパルスを計数するには、差分を用いてパルス波高値相当の値を算出する必要がある。しかし、この差分をとった場合、センサ出力などの正のみ、または負のみの単極性のパルスは、従来のグランドレベルからの波高値と同等の識別レベルとなるが、回路からのホワイトノイズなどの両極性のパルスに対しては、そのピーク間の電圧をパルス識別レベルと識別する。このため、従来のグランドレベルからの波高で識別する場合に比べ、約２倍の識別レベルが必要となる。したがって、センサ出力のみを計数するのに必要な識別レベルが、従来の方式より２倍程度必要となり、センサ信号／ホワイトノイズ比（ＳＮ比）が悪くなる。
【００１４】
第２の課題は、キャンベル計測における耐ノイズ性の向上である。従来のモータ等からのノイズ試験の結果、数μｓ幅のサージノイズが原子炉起動監視装置に誘導しやすいことが分かっている。パルス計測に関しては、前記のディジタル演算によるパルス識別によりこのサージノイズの低減が可能である。しかし、キャンベル計測においては、測定帯域をセンサのパルス形状から選定した数百Ｈｚから１ＭＨｚ以下の周波数帯域に設定しているが、前記の先行技術例では、この測定帯域をシフトすることにより誘導ノイズを除去している。しかし、最も誘導しやすいノイズの周波数が、ほぼ測定帯域の範囲になっており、完全に除去することは難しく、また、センサ感度がわずかに変化するため、これらを補正する必要がある。
【００１５】
一方、一般に線量当量を求める計測装置は、γ線エネルギーに対する感度特性を人体に対するエネルギー吸収特性と等しくなるように、センサの放射線入射窓、反応体積などが最適化されている。しかし、γ線の入射方向によって感度特性が異なる場合があり、感度特性を正確に一致させることは難しい。また、人体に対する線量当量を正確に換算するには、人体の部位によってエネルギー吸収特性が異なるため、γ線エネルギーに対するセンサ感度と一致させただけの計測装置のみでは、人体の各部位に対する正確な評価は難しい。さらに、γ線以外の中性子、たとえばβ線などが混在している場合、センサ構造で感度を補正しているセンサでは、物質への吸収特性が大きく異なるこれらの放射線に対しては評価できないため、それぞれ専用の測定系が必要となる。従来、これらの課題を解決するために、γ線のエネルギースペクトルを求め、その値から換算する手法が提案され、実用化されている。しかし、この手法ではパルス波高を用いたエネルギー情報の取得を元にしており、パルスのパイルアップが生じる条件では、エネルギー情報の取得が難しくなり、精度が悪くなる。つまり、センサの出力パルス幅にもよるが、通常エネルギー測定を行う場合約１×１０５ｃｐｓ程度が計測の上限となる。計測下限は、応答要求を満たす計数とする必要があるが、１ｃｐｓと仮定しても、計測レンジは約５桁程度となり、さらに広レンジの線量を連続で行える測定方法が望まれている。
【００１６】
特開平３−１８３９８３号公報には、センサを二重構造にすることにより測定の精度を改善することが記載されている。この手法では、前記したパルス計測のパイルアップの影響は、電流計測を行うことによって回避されているが、この手法に関しても、センサ構造または処理内容が複雑なため、その簡素化が望まれていた。
【００１７】
本発明はかかる従来技術に対処してなされたものであり、ディジタル処理を用いたパルス計測およびキャンベル計測において、耐ノイズ性を改善する放射線計測装置を提供することを目的とする。
【００１８】
【課題を解決するための手段】
上記目的を達成するため、請求項１に係る発明の放射線測定装置は、放射線に感応してパルス信号を出力する放射線検出手段と、この放射線検出手段の出力パルス信号を整形しサンプル値に変換するＡ／Ｄ変換手段と、サンプル値のｎ乗値相当（ｎ＞＝２）を演算し放射線のパルス信号のみを識別するｎ乗パルス識別手段と、このｎ乗パルス識別手段によって識別された放射線のパルス信号を計数する計数手段とを備え、しかも前記ｎ乗パルス識別手段は、前後ｎ個のサンプル値を掛け算することにより前記ｎ乗値相当を演算するか、又は同じサンプル値をｎ乗した後に前後２個のｎ乗値を平均化処理することにより前記ｎ乗値相当を演算するものであることを特徴とする。
【００１９】
請求項２に係る発明は、請求項１に記載の放射線測定装置において、ｎ乗パルス識別手段は、放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、パルス立ち上がり部およびパルス立ち下がり部に対応する両サンプリング点におけるサンプル値の和のｎ乗値を、両サンプリング点の間に位置する両サンプリング点におけるサンプル値の和のｎ乗値から差し引くことにより出力パルスの波高に相当する差を演算し、その差から所望のパルスのみを識別する差分識別手段を備えていることを特徴とする。
【００２０】
請求項３に係る発明は、請求項１に記載の放射線測定装置において、ｎ乗パルス識別手段は、放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、前半２つのサンプリング点におけるサンプル値の差のｎ乗値と、後半２つのサンプリング点におけるサンプル値の差のｎ乗値との和を演算し、その和に基づいて所望のパルスのみを識別する差分識別手段を備えていることを特徴とする。
【００２１】
請求項４に係る発明は、請求項１に記載の放射線測定装置において、ｎ乗パルス識別手段は、放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、各サンプル値のｎ乗値の積算値と放射線検出手段の出力パルスの幅に相当するサンプル値の波高との比率を算出し、その比率から放射線検出手段の出力パルスを識別する波高・パワー識別手段を備え、計数手段は波高・パワー識別手段によって識別されたパルスを計数することを特徴とする。
【００２２】
請求項５に係る発明は、請求項１に記載の放射線測定装置において、正および負の両極性の入力信号を単極性の信号に変換する単極性変換手段をさらに具備し、かつｎ乗パルス識別手段は、単極性変換手段によって得られた単極性の信号を偶数ｎを用いてｎ乗演算を行うｎ乗演算手段と、このｎ乗演算手段の出力を平均化する平滑手段とを備えていることを特徴とする。
【００２３】
請求項６に係る発明の放射線測定装置は、放射線に感応してパルス信号を出力する放射線検出手段と、この放射線検出手段の出力パルスを整形し周波数帯域を制限する帯域制限手段と、この帯域制限手段の各出力パルスのｎ乗値（ｎ＞＝２）を算出するｎ乗値演算手段と、このｎ乗値演算手段の出力信号を所定の時間幅で平滑化する第１の平滑化手段と、この第１の平滑化手段の出力の大小関係を評価し、その評価結果から特定のデータを除くデータ除去平均化手段と、このデータ除去平均化手段からの出力をさらに平滑化する第２の平滑化手段と、この第２の平滑化手段の出力を放射線強度に換算する変換手段とを具備する。
【００２４】
請求項７に係る発明は、請求項６に記載の放射線測定装置において、第１の平滑化手段の平滑時間幅を外来ノイズのパルス幅以上に設定し、データ除去平均化手段の平滑化区間を外来ノイズパルスの到来する時間間隔以下に設定し、計測されたノイズ波形によりこれら平滑時間幅とデータ除去数を調整するノイズ特性評価手段をさらに具備したことを特徴とする。
【００２５】
請求項８に係る発明は、請求項６に記載の放射線測定装置において、データ除去平均化手段は、第１の平滑化手段の過去の出力と現在の出力とを比較し、現在の出力値が第１の平滑化手段の揺らぎから過去の出力値をもとに事前に評価したしきい値以上に増加していた場合、第１の平滑化手段の現在の出力を除去し、予め評価によって求めた評価比率を過去の出力値に乗じた値に置き換えるデータ制限手段を有することを特徴とする。
【００２６】
請求項９に係る発明の放射線測定装置は、ｎは３以上の奇数であり、ｎ乗値演算手段は帯域制限手段の各出力パルスに対して奇数次のモーメントを算出することを特徴とする。
【００３０】
【発明の実施の形態】
＜第１の実施の形態＞
本発明の第１の実施の形態による放射線測定装置を、図１を参照して説明する。
【００３１】
図１に示す放射線測定装置は、放射線に感応し放射線量に応じた出力パルスを発生するＳＲＮＭセンサ１と、その出力パルスを増幅するプリアンプ２Ａと、プリアンプ２Ａの出力パルスをそのパルス幅より狭い時間間隔でサンプリングしサンプル値を得るＡ／Ｄ変換器３と、そのサンプル値をｎ乗しＳＲＮＭセンサ１の出力パルス幅に相当するｎ乗値を用いて信号を識別するｎ乗パルス識別装置４と、このｎ乗パルス識別装置４により識別されたパルスを計数するパルス計数装置（ＰＣ）５とを備えている。ＳＲＮＭセンサ１は、パルス出力が得られる核***センサであるが、これは同様のパルス出力が得られる電離箱等に置換することもできる。
【００３２】
このように構成された放射線測定装置のＳＲＮＭセンサ１に中性子が入射し、内部で核***すると図２（ａ）に示すような電気パルスが出力される。このＳＲＮＭセンサ１の出力パルスの時間幅は１００ｎｓ程度である。この出力パルスは、プリアンプ２Ａに入力され、増幅される。プリアンプ２Ａは、ＳＲＮＭセンサ１へ動作電圧を印加する機能をも有する。プリアンプ２Ａの出力パルスは、Ａ／Ｄ変換器３に入力され、図２（ａ）で点により示されるようなサンプリング時間間隔でサンプリングされ、ディジタル化される。このサンプリング時間間隔は短いほど波形情報を多く抽出でき、他の誤差信号、たとえば外来ノイズによる信号を除いてセンサ出力パルスのみを正確に計数することが可能となる。Ａ／Ｄ変換器３では、サンプリング前に、サンプリング定理で必要とされる周波数帯域に制限するバンドパスフィルタ処理も行う。
【００３３】
Ａ／Ｄ変換器３から出力されるサンプル値は、ｎ乗パルス識別装置４に内蔵されるｎ乗演算部に入力され、ここでｎ乗演算される。つまり、各サンプル値をそれぞれｎ乗演算するか、または、前後ｎ個のサンプル値を掛け算する。ここで、ｎ＞＝２とする。同じサンプル値をｎ乗演算した場合は、前後２個のｎ乗値を平均化処理する。たとえば、８個でパルス波形をサンプリングした場合は、おのおのの２乗値を演算した後、前後２個の演算値を平均化処理し、結果的に４つの演算値を得る。この場合、移動平均処理を行い、８個のサンプル値を得るようにしてもよい。
【００３４】
例えばｎ＝２として、２乗演算を行った場合のサンプル値の変化の様子を図２（ｂ）に示す。このようにサンプル値をｎ乗値に変換することにより、パルス波形の波高と回路によるノイズ成分の比率は、従来の比率のｎ乗倍に改善することができる。
【００３５】
ただし、単純にサンプル値ごとにｎ乗演算した値をパルス識別に用いても、サンプル値をそのままパルス識別に用いた場合と同じ識別性能となる。そこで、ディジタル演算でｎ乗する場合は、すでに述べたように前後のデータを掛け合わせるか、ｎ乗後に前後のデータを平均化する処理を加えることが必ず必要となる。この２乗演算された結果を、予め設定しておいた識別レベルと比較することにより、ＳＲＮＭセンサ１の出力パルスと回路ノイズとの識別が容易になる。この放射線測定装置によりＳＲＮＭセンサ１のパルスとして識別されたものは、パルス計数装置６でパルス発生率に換算され、最終的にはＳＲＮＭセンサ１の位置での中性子束レベルに換算される。
【００３６】
この実施の形態によれば、回路ノイズレベルと同程度の波高を有するパルスを従来の差分を用いた手法よりも良好に識別することができる。図３（ａ）にホワイトノイズと、それと同程度のパルス波高を有するセンサ出力パルスをＡ／Ｄ変換器３でサンプリングした場合の一例を示す。この図３（ａ）では、垂直破線で囲まれた時間４．２０×１０３ｎｓの近傍においてセンサ出力パルスが発生した場合を示しており、この区間で回路ノイズとセンサパルスとが重畳している。このデータに対して、従来のパルス波高によるノイズ識別を行った場合、図３（ａ）中の水平破線を識別レベルLbとすると、両パルスとも計数してしまい、正確な測定をすることはできない。また、上記の差分を用いたパルス計数方式では、図３（ａ）中、垂直破線で囲まれた範囲における正負両ピーク間の電圧がパルスの波高と認識される。したがって、正と負の両極性で生ずる回路ノイズは、従来の０Ｖからの波高値で識別していたものに比べ、２倍の電圧まで識別レベルをあげないと識別することができない。一方、加算平均化することにより、両極性のノイズは正負の信号がキャンセルされて平均化され、単極性のセンサ出力はもともと正のみ又は負のみの信号であるため、このようにキャンセルされることはないが、図３（ｂ）に示すように、パルス幅が広くなる。パルス幅が広くなった場合、センサ出力のパルス数の多い条件ではパルスが重なりパルスを計数できなくなり、パルス計測の計数上限が低下するという問題がある。
【００３７】
本実施の形態のｎ乗パルス識別装置４によるｎ乗演算処理後のサンプル値を図３（ｃ）に示す。回路ノイズ相当のホワイトノイズ（左側）の識別レベルとセンサ出力（右側）の識別レベルの比率は、従来の差分を用いた場合より大きく改善することができる。また、パルス幅も、平均化のみの処理をした場合の図３（ｂ）に比べ広がることもなく、これにより、パルス計数の上限を悪化させることなくパルス計測を可能とすることができる。特に、本実施の形態においてｎが奇数の場合、両極性ノイズの符号を維持して平均化することができるため、両極性信号と単極性の信号を良好なＳＮ比（信号雑音比）をもって識別することができる。
【００３８】
以上により、パルス計測において、出力パルスの識別レベルを設定する際に、回路ノイズまたはα線ノイズを除去する識別レベルを低く設定することができ、その分、センサ出力パルスが小さい場合でも計測感度を落とすことなくパルス計測を行うことができる。
【００３９】
＜第２の実施の形態＞
次に、本発明の第２の実施の形態を図４に基づいて説明する。
【００４０】
図４の放射線測定装置のｎ乗パルス識別装置４は、積分識別装置６、差分識別装置７、および波高・パワー識別装置８を備えている。Ａ／Ｄ変換器３の出力パルスは積分識別装置６および差分識別装置７に入力され、両者の出力は波高・パワー識別装置８に入力される。波高・パワー識別装置８の出力信号はパルス計数装置５に入力される。ここで、積分識別装置６は、第１の実施の形態で説明したパルス識別方式によりパルスを識別する。
【００４１】
以下、本実施の形態における差分識別装置７が行う処理の第１実施例について説明する。この実施例では次の原理により識別を行う。
【００４２】
Ａ／Ｄ変換器３の出力の最大値Top(k)およびボトム値Bottom(k)をそれぞれ、
Top(k) = +b×S(k-2) + c×S(k-1)
Bottom(k) = a×S(k-3) + d×S(k)
とすると、波高値High(k)は、先行技術での式（１）を簡略化し、
High(k) ＝ ＋Top(k) - Bottom(k)
とすることができる。
【００４３】
本実施の形態では、まず最大値Top(k)の２乗値とボトム値Bottom(k)の２乗値の差を求める。すなわち、
Top(k) >> +Bottom(k)
であるため、
X = High(k) × Top(k) …（３）
一方、非常に大きなサージノイズに信号パルスが重畳した極端なケース（以下、ケース２とする）では、近似値として、
Top(k) = Bottom(k)
であるため、式（２）から、
X = High(k) × (2×Top(k)) …（４）
となる。式（３）および式（４）から、
X／Top(k) = α × High(k) …（５）
（ここで、ケース１ではα＝１、ケース２ではα＝２）
となり、X／Top(k)はほぼパルス波高High(k)の一次関数となる。つまり、サージノイズ上にＳＲＮＭセンサの出力パルスが重畳した場合でも、この値によって識別することにより、数μｓの周期のサージノイズ上に重なった数百ｎｓのＳＲＮＭセンサ出力パルスを識別し、計数することが可能となる。
【００４４】
このように本第１実施例によれば、ｎ乗値の差を用いることによって、仮にセンサパルスよりもパルス幅の広い外来ノイズが誘導されても、その影響を低減することができる。
【００４５】
次に、本実施の形態における差分識別装置７の処理の第２実施例について説明する。式（１）を、
D1(k) = c×S(k-1) - d×S(k) …（６）
D2(k) = a×S(k-3) - b×S(k-2) …（７）
とすると、波高値High(k)は、
High(k) = +D1(k) + D2(k) …（８）
となる。式（８）の右辺各項のｎ乗の和Yは、
Y = D1(k)ｎ + D2(k)ｎ
となり、これは近似値として、
Y = High(k)ｎ
すなわち、
Y−ｎ = High(k) …（９）
となる。よって、この場合は、近似値としてY−ｎを指標に用いることにより、パルス識別が可能である。なお、ｎ＝１のときは、式（１）に相当する。
【００４６】
以上、差分識別装置７により、サージノイズ上にＳＲＮＭセンサ１の出力パルスが重畳した場合でも、差分によって求めた指標を用いることにより、数μｓの周期のサージノイズ上に重なった数百ｎｓのＳＲＮＭセンサ１のパルスを識別し、計測することが可能となる。
【００４７】
図５に積分識別装置６および差分識別装置７の識別性能の差異をまとめて示す。積分識別装置６は、ホワイトノイズなど両極性ノイズで、センサ出力のパルス幅より短い周期のノイズを低減するのに有効であり、一般に回路ノイズ、センサのα線ノイズなど、センサ出力より信号（波高値）が小さいノイズをさらに低減するのに有効である。一方、差分識別装置７は、センサ１の出力パルス幅より長い周期のノイズに対して有効であり、一般に外部からの数μｓパルス幅の誘導ノイズを除去するのに有効である。したがって、両者のロジックを最適に調整し、両者がともに成立した場合のみパルスを計数することにより、センサ出力のみをより正確に計数することができる。
【００４８】
以上のように、本第２実施例によれば、波高値に相当するサンプル値の差のｎ乗値を用いることにより、第１実施例と同様に、センサパルスよりもパルス幅の広い外来ノイズが誘導されても、その影響を低減し、より高精度の放射線測定を行うことができる。
【００４９】
次に、本実施の形態の第３実施例として、波高・パワー識別装置８について説明する。波高・パワー識別装置８は、積分識別装置６からパルスの積分値を、差分識別装置７からパルスの波高値相当の値を受け取る。この両者の比率（積分値／波高値）は、センサ出力パルスの場合、ほぼパルス幅相当となり一定の値を示す。一方、高周波成分を含むホワイトノイズは、波高値がセンサ出力と同等の場合でも積分値が小さいため、比率は小さくなる。また、パルス幅の広いサージノイズは、積分値は大きく波高値は小さいため、その比率はセンサ出力より大きくなる。したがって、この比率が一定範囲内である場合に、パルスを計数することにより、これらのノイズの影響を低減することができる。
【００５０】
以上により、パルス計測において、サージ状の外来ノイズが誘導された場合でも、その周期が数μｓとセンサ出力パルス幅１００ｎｓより長い場合、そのサージ状ノイズの影響を除去し、そのノイズに重畳したパルスも計数することができる。かくして、本第３実施例によれば、本実施の形態の第１実施例または第２実施例で述べた差分によるパルス計数と、第１の実施の形態で述べたｎ乗値によるパルス計数を併用し、その比率を用いることにより両者の特性を同時に実現し、より正確なパルス計測が可能となる。
【００５１】
＜第３の実施の形態＞
次に、本発明の第３の実施の形態を図６に基づいて説明する。
【００５２】
図６に示す放射線測定装置は、図１に示した放射線測定装置において、Ａ／Ｄ変換器３とｎ乗パルス識別装置４との間に、両極性の信号をパルスの主成分である極性のみの信号に変換する単極性変換手段９を介挿したことを特徴とするものである。なお、ＳＲＮＭセンサ１はパルス出力が得られる核***検出器であるが、これは他のパルス出力が得られる電離箱等の放射線検出器であってもよい。
【００５３】
このように構成された本実施の形態における単極性変換手段９の機能について図７を参照して説明する。図７（ａ）は、図２（ａ）に示されるような検出器出力を２次微分処理した場合のパルス波形の一例を示す。２次微分処理は、プリアンプ２Ａにおいてアナログ回路で処理する場合と、図２の波形をＡ／Ｄ変換器３でサンプリングし、ディジタル処理で２次微分演算を行う場合が想定されるが、今回は両者とも想定している。この図７（ａ）をｎ乗パルス識別装置４において偶数ｎを用いてｎ乗演算後、平均処理した結果を図７（ｂ）に示す。両極性が単極性に変化するために、パルス幅が広くなる。
【００５４】
そこで、単極性変換手段９は、図７（ａ）のパルスの主成分は負極性であるため（図７（ａ）の下方が負極性）、正極性の成分をゼロまたはそれに近い数値に置き換え、図７（ｃ）の波形に変換する。この波形を、ｎ乗パルス識別装置４に内蔵されるｎ乗演算手段によって、偶数ｎを用いてｎ乗演算した後、たとえば２乗演算した後、ｎ乗パルス識別装置４に内蔵される平滑手段によって平均化処理した結果を図７（ｄ）に示す。これから分かるように、パルス幅の広がりを、単極性変換手段９を用いない図７（ｂ）に比べ狭くすることができる。
【００５５】
かくして本実施の形態によれば、パルス計測において、偶数ｎを用いてｎ乗演算を行ってもパルス幅が広がることなく、その結果、パルスの重なりによって１個のパルスを計数できなくなる、パルスパイルアップによる誤計数を低減することができる。
【００５６】
＜第４の実施の形態＞
次に、本発明の第４の実施の形態を図８に基づいて説明する。
【００５７】
ＳＲＮＭセンサ１の出力は、帯域制限機能を有するプリアンプ２Ｂで増幅・帯域制限され、Ａ／Ｄ変換器３に入力される。Ａ／Ｄ変換器３の出力は、バンドパスフィルタ（ＢＰＦ）１０に入力され、特定の周波数帯域に制限される。通常、バンドパスフィルタ１０の周波数帯域はＳＲＮＭセンサ１の出力特性によって決定されるが、たとえば１００ｋＨｚから４００ｋＨｚまでの範囲の周波数成分のみを通過するディジタルフィルタとして構成することができる。この測定帯域は複数個設定してもよいが、本実施の形態では一つの帯域に特定した単純な場合の構成で代表させる。
【００５８】
ディジタルフィルタ処理を行う際、入力のサンプル値は、バンドパスフィルタ１０の出力周波数帯域に対して適切なサンプリング間隔およびビット数に調整される。このバンドパスフィルタ１０の出力はｎ乗値演算装置１１に入力され、ｎ次モーメント値に変換される。
【００５９】
本実施の形態では、ｎ乗値演算装置１１の出力側に、第１平滑装置１２Ａ、中間値の平均を取るデータ除去平均化装置（ＤＥＡ）１３、および第２平滑装置１２Ｂが配設され、その平滑個数とデータ除去個数を外来ノイズのサージノイズ幅と到来周期およびＭＳＶ計測に必要なデータ除去可能率を元に調整することによって、外来ノイズの影響を受けないＭＳＶ計測が可能となる。これによってさらに信頼性の高い計測を実現することができる。
【００６０】
さらに、データ除去平均化装置１３において、最大値のみを平均化することにより、低計数率時にもＭＳＶ指示値を得ることができ、プリアンプ２Ｂの回路ノイズ等で決まっているＭＳＶ計測の測定下限を拡大することができる。
【００６１】
なお、本実施の形態においては、ｎを奇数に設定した場合、サージ状のノイズは正負両極性のノイズであるため、平均化処理することによって正負の値が相殺され、単極性のセンサ出力と同じ極性の値を出力として選択することにより、サージ状ノイズの影響を低減することができる。
【００６２】
以下に述べる実施の形態では、ｎ次モーメント演算として最も簡単な２乗値に代表させて説明する。つまり、バンドパスフィルタ９のｋ番めの出力サンプル値をS(k)とし、ｎ乗値演算装置において、
Out1(k) = S(k) × S(k) …（１０）
の演算が行われるものとして説明する。
【００６３】
この出力Out1(k)は、第１平滑装置１２Ａに入力され、ある個数分について平均化される。この際、バンドパスフィルタ１０の出力の平均値がオフセットを有する場合は、サンプル値S(k)の平均値も同時に算出し、その２乗値を算出し、２乗演算結果Out1(k)から引き算する。つまり、ｎ個平均化する場合は、
Out2(k2) = (ΣOut1(k))/n − {(Σs(k))/n }２ …（１１）
を演算する。ただし、Σはｎ個のサンプル値を加算することを示すものとする。
【００６４】
この際、加算するｎの個数は、想定される外来ノイズのパルス幅に相当する個数分に設定され、また、この個数はディジタル演算のしやすいように２のｎ乗個に設定される。つまり、図９に示すようなノイズを想定した場合、すなわちパルス幅２０μｓのサージ状パルスが２ｍｓ周期で到来するパルス状ノイズを想定した場合、１μｓ周期でバンドパスフィルタ９の出力が得られるものとして、２０μｓ以上のデータを平均化する。ただし、ディジタル演算で行う場合、ビットシフトによる割り算が行えるように２のｎ乗個に加算数を設定した方が演算に有利であるため、このケースでは３２（＝２５）個加算する。出力間隔ｋ２は、バンドパスフィルタ９の出力間隔１μｓに対して、その３２倍である３２μｓ周期となる。
【００６５】
この第１平滑装置１２Ａの出力は、データ除去平均化装置１３に入力される。データ除去平均化装置１３では、３２μｓ周期の第１平滑装置１２Ａの出力を、特定個数ごとに区切って、その大小関係を比較し、特定のデータを除く処理を行う。図９に示すようにサージノイズの除去を想定した場合は、以下の処理を行う。図９に示すサージ状パルスの到来間隔は、約２ｍｓである。よって、２ｍｓ以下の間隔でデータを区切り、その中でサージ状パルスをサンプリングして得たパルスデータのみを除けば、残りはノイズのないデータのパワー値に換算することができる。つまり、この区間の大小関係を比較し、値の大きいものから２個以上を除けば、このサージ状ノイズの影響を除去することができる。ここで最低の除去数を２個に設定したのは、タイミングによっては２０μｓ幅のサージノイズが２個分の第１平滑装置１２Ａの出力信号に混入する可能性があるためである。また、大きい値を除いた場合、平均値を保持するために、小さい値も除いた最大の個数と同数、つまり、２個分の小さい値のデータを除く必要がある。このように大きい値と小さい値を２個ずつ除いた残りのデータを平均化する。図９のノイズの場合は、３６個分の第１平滑装置１２Ａの出力をひと区切りとすると、約１．１６ｍｓの平滑区間となり、その中に含まれるサージパルスは最大２個であるため、大きいもの２個と小さいもの２個を除き、残りの３２個を平均化する。この場合も、除いた後の個数が２のｎ乗個となるように、加算数を決定する。
【００６６】
また、サージ状パルスの到来間隔が短くなったケースでは、データ除去平均化装置１３での加算数を少なくする場合と、加算平均区間内でのサージ状パルスの到来個数を評価し、その個数の２倍分の最大・最小データを除く処理を行えば、ノイズを除去することができる。
【００６７】
ただし、このようにデータを除く場合に問題となるのは、データの除去率である。つまり、ＳＲＮＭセンサ１の出力パルス数が少ない場合、データを除去することによって、平均したパワーつまりＭＳＶ値が実際の値より低めに表示される。そのパルス数が充分な場合、そのランダム性からある程度のデータ除去は誤差範囲内におさまる。したがって、その測定精度から得られる除去限度率を確保するようにデータ除去を行えば、測定に悪影響を与えることなしに、間欠的に誘導されるサージノイズによる突発的な指示値の上昇を、あるしきい値の範囲内の値に制限することができる。
【００６８】
データ除去平均化装置１３の他の機能として、最大値のみを選択することにより、ＳＲＮＭセンサ１のパルス発生率が低い場合にも、ＭＳＶ指示値を得ることが可能となる。つまり、低いパルス発生率の場合、ひとつの平滑化区間においてセンサパルスが１個も到来しない時間が多く存在する。したがって、この到来しない区間を平滑化に含めないようにすることにより、パルス発生率が低い場合もＭＳＶ指示値を得ることが可能となり、ＭＳＶ計測の下限を延長させることができる。これは、パルス計測とＭＳＶ計測を併用した計測手法であるということができる。ただし、この場合、パルス発生率に対するＭＳＶ指示値の直線性を補正するために、データの除去割合をパラメータとした補正関数を用いて補正する必要がある。
【００６９】
次に、データ除去平均化装置１３内でデータを除去する処理の基準に関する一実施例を以下説明する。模擬中性子パルスの入力に対し第１平滑装置１２Ａにおいて３２μｓで平滑化した場合のシミュレーション結果の一例を図１０に示す。図１０の実線は図の縦左軸に対応し、平滑後のＭＳＶ指示値の平均値Ｓを示している。これから分かるように、パルス計数率がＭＳＶ計測範囲内すなわち１×１０5ｃｐｓ以上においては、概ね、パルス計数率と平滑後のＭＳＶ指示値すなわち２乗電圧とは、ほぼ比例して推移する。他方、パルス計数率が１×１０5ｃｐｓ〜２×１０5ｃｐｓの範囲では回路ノイズの影響で比例関係にはなっていない。
【００７０】
また、３２μｓでのＭＳＶ指示値の平均値Ｓと標準偏差σに対して、
Ｘ＝（Ｓ＋６σ）／Ｓ
で表される指標をパルス発生率ごとに評価したものを、図１０の縦右軸に対応する破線で示す。ＸはＭＳＶ計測下限で最大５．３程度となり、同様の評価で１×１０6ｃｐｓ以上では２〜３程度で推移する。
【００７１】
一方、他の評価方法として、バンドパスフィルタ１０で１００〜４００ｋＨｚの帯域制限処理を施して第１平滑装置１２Ａで模擬中性子パルスを３２μｓで平滑化したときのＭＳＶ計測値の揺らぎ率Ｙ、すなわち標準偏差σと平均値Ｓとの比率：Ｙ＝σ／Ｓ の推移を理論式で評価した場合のグラフを図１１に示す。図から分かるように、ＭＳＶ計測範囲（計測下限１×１０5ｃｐｓ）におけるＹの値は０．４以下であり、ＭＳＶ計測範囲の下限で揺らぎは最大となる。
【００７２】
これらの図から分かるように、波形の揺らぎは、特にＭＳＶ計測範囲の下限である約１×１０5ｃｐｓ程度において大きい。したがって、特にこの揺らぎの大きい計数率：１×１０5ｃｐｓ程度における、上述したＸあるいはＹにより例示される揺らぎの度合いを示す指標を、第１平滑装置１２Ａの積分時間の条件で事前に評価し、得られた評価値を超えたものはノイズと判定することとする。
【００７３】
すなわち、第１平滑装置１２Ａでパルスを３２μｓで平滑化した場合、第１平滑装置１２Ａで以前に平滑処理された結果や、あるいは予め模擬中性子パルスを用いて評価した結果を用いて得られる評価指標および評価基準と、平滑化した結果との比較を行い、この評価指標に基づいてパルスを評価したときに、評価基準を超えるような場合、例えば評価指標が前回の平滑処理結果の８倍と規定されるしきい値を超えた場合、この部分を異常値と判定する。
【００７４】
図１２は、ノイズを含むパルス波形を模式的に示したものである。通常のＭＳＶ値は破線で囲まれる範囲を揺らいでいる。この範囲を超えるものはノイズである。このしきい値は、その平均のＭＳＶ指示値によって変化するが、もっとも変動の大きいＭＳＶの計測下限での評価値を用いれば、全ＭＳＶ計測範囲において保守的な評価として適用可能である。
【００７５】
また、ここでのデータ除去平均化装置１３においては、ノイズと判定された部分のデータについては、このデータ自体を除去する場合と、このデータを上述したしきい値の範囲内の値に置き換える場合とがある。前者の場合、ノイズをほぼ完全に除去することができるが、上述したデータ除去率の範囲で除去する必要があるためデータ除去率を評価する必要となる。逆に後者の場合、完全なノイズ除去はできないが、データ除去率の評価が不要であるという特徴がある。
【００７６】
以下、後者の場合、すなわちノイズ部分の値を置換する場合の一例を詳述する。図１３は、ＭＳＶ指示値の揺らぎの範囲と平均値の変化幅の関係を示すグラフである。図中の実線はパルスを32μｓで平滑化したときのＭＳＶ信号出力、破線は中性子束の最大変化を示している。図中一点破線で示された、上述の方法により評価されるＭＳＶ計測の揺らぎの変化幅は、破線で示された、本来計測すべき中性子束の変化率に比べて十分大きい。よって、この評価されたＭＳＶ計測の揺らぎの変化幅を超えた場合を、異常値と判定することとする。
【００７７】
この揺らぎの変化幅すなわち最大揺らぎを超えるノイズが誘導された場合のパルスのＭＳＶ指示値の推移の一例を図１４（ａ）に示す。図中矢印Ａは、ＭＳＶ計測の揺らぎの変化幅、例えば前回の平滑処理結果の８倍と規定されたしきい値を示す。このＡを超えている部分は異常値と判定され、この異常値を、正常値として、本来の中性子束の最大変化率を前回のサンプリング値に乗じた値に置き換えることとする。置換処理の結果を図１４（ｂ）に示す。ここで、３２μｓ幅での想定される中性子束の最大変化率は、例えば1.03倍程度であり、上述したＭＳＶ計測の揺らぎ範囲よりも十分小さいものである。よってここでの置換処理は、異常値を前回値の1.03倍に置き換えるものとしている。このように、異常値が検出された場合の変化率も、平滑区間内での監視すべき中性子束の最大変化率をあらかじめ評価しておくことにより、計測の時間応答性を損なうことなく、異常値のみを除去することができる。
【００７８】
次に、データ除去平均化装置１３の出力は、第２平滑装置１２Ｂに入力され、測定値の揺らぎが必要とされる測定精度を満たし、また、応答要求を満たす範囲となるように、平均化処理される。この結果は、ＭＳＶ中性子評価装置１４に入力され、測定されたＭＳＶ値が中性子束の値に換算される。
【００７９】
また、本実施の形態では、ノイズ特性評価装置１５を設けるのが好適である。ここでは、ノイズ波形の特性であるサージ状パルスのパルス幅と到来周期の最小値を評価し、ＭＳＶ計測の測定範囲において必要となる最大データ除去率の範囲内において、第１平滑装置１２Ａで平均化のために用いるデータ数、データ除去平均化装置１３のデータ平均区間と、データ除去数、および、第２平滑装置１２Ｂの平滑化フィルタの時定数を設定する。
【００８０】
このような構成により、ＭＳＶ計測において、サージ状の外来ノイズが誘導された場合でも、そのサージ状のパルス幅と、その到来周期を評価し、ＭＳＶ計測の許容データ除去率を満たす範囲内でサージノイズデータを除去することにより、間欠的なサージノイズを完全に除去することができる。
【００８１】
＜参考例＞
次に、本発明の参考例を図１５に基づいて説明する。
【００８２】
図１５に示す放射線測定装置では、放射線センサとして、常温半導体であるＣｄＴｅを用いたＣｄＴｅセンサ１６を用いている。放射線センサとしては、エネルギー情報の取得が可能な、ＮａＩ等のシンチレーションセンサと光電子増倍管を組み合わせたもの、半導体センサであるＧｅセンサなども同様に用いることが可能である。ＣｄＴｅセンサ１６の出力は、チャージアンプ（ＣＡ）１７に入力される。チャージアンプ１７は、入力パルスの電荷を積分し、その電荷量に応じたパルス波高を有するパルスに変換し出力する。なお、チャージアンプ１７はＣｄＴｅセンサ１６に対しその動作電圧を供給する。チャージアンプ１７の出力は、一般に放射線のエネルギー測定に用いられるパイルアプリジェクション回路、ポールゼロキャンセル回路などで波形整形処理が施され、ＭＳＶ計測装置１８、電流計測装置（ＣＤ）１９、およびパルスカウンタ（ＰＣ）２０に入力される。ＭＳＶ計測装置１８では、周波数帯域の制限を行った後、ｎ乗の平均化演算が行われ、ＭＳＶ計測値つまり２次モーメント値に変換される。電流計測装置１９では、１次モーメント値である平均電流値が測定される。また、パルス計測装置２０では、パルス数の計数が行われる。これらＭＳＶ計測値、電流計測値、およびパルス計数値は、それぞれエネルギー評価装置２１に入力され、ＭＳＶ値とパルス計数の比率、または、ＭＳＶ値と直流電流値の比率（１次と２次のモーメント値の比率）をもとに平均放射線エネルギーを評価する。この平均のエネルギー値と、前記の測定値は、線量評価装置２２に入力され、照射線量、または、物質での吸収線量、または、人体へのリスク率を含めた線量当量に換算される。
【００８３】
エネルギー評価装置２１における放射線の平均エネルギーの推定方法を以下に説明する。
【００８４】
チャージアンプ１７の出力は、ＣｄＴｅセンサ１６内で放射線が吸収されたエネルギーに比例した波高値を有するパルスである。したがって、その反応の起こる確率をＮとし、吸収エネルギーをｑとすると、ＭＳＶ値、パルス計数値、および電流値は、以下の式で近似することができる。
【００８５】
ＭＳＶ値 ： k1×q２×N
ｎ次モーメント値 ： kn×qｎ×N
パルス計数値 ： k2×N
直流電流値（１次モーメント値） ： k0×q×N
ただし、k0、k1、k2、knはそれぞれ補正係数である。
【００８６】
よって、これらの比率は、
ＭＳＶ値／パルス計数値＝(k1／k2)×q２（一般に、kn×qｎ）
ＭＳＶ値／直流電流値＝(k1／k0)×q
ｎ次／ｎ’次モーメント値 ＝ kn×q（ｎ−ｎ’）／kn’
となり、これらの補正係数k0，k1，k2，kn等を予め評価し、これら測定値の比率を用いることにより、結晶中での吸収エネルギーの推定が可能となる。
【００８７】
図１６は、ＣｄＴｅセンサ１６によって種々の核種からの放射線、つまり、エネルギーの異なる放射線、を測定した場合の、パルス計数値およびＭＳＶ値（縦軸）を、市販の放射線サーベイメータで測定した線量（横軸）に対してプロットしたものである。一般に、サーベイメータ等は、放射線のエネルギーに対する感度特性を、放射線のエネルギーに対する線量当量の評価曲線に一致するように遮蔽、または、内部の補正係数を調整している。つまり、エネルギーの低い放射線に対しては、１個の放射線が入射した場合の結晶中での吸収エネルギーが小さいためパルス計数率が大きくなるが、１個の放射線の線量は低い値を示す。逆に、エネルギーの高い放射線に対しては、パルス計数率は低くなるが、１個の放射線で発生する電荷量が大きいため線量は大きくなる。したがって、測定しているパルス計数値または電流値の感度を線量当量への寄与率と同じになるように遮蔽等で調整する。図１６は、これらの感度補正を行っていないため、低いエネルギーの放射線に対してはパルス計数値が高くなっており、線量に対してパルス計数値およびＭＳＶ計測値はランダムになっている。
【００８８】
しかし、これをＭＳＶ値／パルス計数値の線量に対する特性としてプロットすると、図１７に示すように線量当量に対して単調な特性となる。したがって、この特性を事前に評価しておくことによって、両者の比率から線量当量への換算が可能である。また、同様に、このＭＳＶ値／パルス計数値の比率は、入射エネルギーに対しても単調な特性となり、この特性も評価しておくことにより、平均の入射放射線エネルギーを推定することができる。この場合、放射線のエネルギーに対する人体の各部位の吸収特性を用いることにより、各部位で吸収線量をより正確に評価することが可能となる。
【００８９】
さらに、パルス計測で線量評価を行う場合、パルスの波高分布を測定し、入射エネルギーのエネルギー情報を取得し、線量に換算する方法と、前記のセンサ自体の構造を工夫し、パルス計測の感度を線量応答特性と同一にする場合がある。また、電流値で線量を評価する方法としては、後者のセンサ構造を工夫し、感度応答を調整する方法がある。したがって、これら周知の手法と本参考例の線量評価方法を併用することにより、さらに精度のよい線量評価を行うことができる。つまり、たとえば、センサの感度特性を単独にある程度調整し、その後に本発明の補正関数を評価して用いることにより、さらに正確な線量評価が可能となる。
【００９０】
さらに、パルスがパイルアップし、パルス計数の数え落としが生じるような高い計数率では、波高情報を取得し線量に換算する前者の方法では正確な線量は評価できない。また、遮蔽を設けた後者の場合も、数え落とし分を補正する必要がある。しかし、本参考例のようにパルス計測とＭＳＶ計測を同時に実施することにより、パルス計測がパイルアップによって飽和した場合でも、ＭＳＶ計測によって線量の評価が可能となり、広いレンジの測定が可能となる。ただし、この場合の平均エネルギーの推定は、パルス計測のパイルアップの効果を補正する必要があるが、ある程度センサ自体の感度を線量応答に近似させておくことにより、その誤差は無視できる範囲に抑えることができる。
【００９１】
また、ＭＳＶ値／パルス計数値の比率と同様に、電流値／パルス計数値の比率と、上述した本参考例によるＭＳＶ値／電流値の比率を用いた場合も、同様の補正関数を求めておくことにより平均放射線エネルギーの推定が可能となる。
【００９２】
かくして、本参考例によれば、線量評価において、ｎ次モーメント値とパルス計測を併用することにより、平均入射エネルギーの推定から、より正確な線量評価を行うことができる。また、パルスがパイルアップし、計数誤差が生じる条件においても、ｎ次モーメントどうしの比率を用いることにより、同様に平均エネルギーの推定が可能となり、従来よりも広い測定レンジで正確な線量評価を行うことができる。
【００９３】
さらに、本参考例の変形例を以下説明する。ここでは、上述したＭＳＶ値の代わりに、１乗値（平均電流）、２乗値、３乗値、…、ｎ乗値をそれぞれ算出し、それぞれの補正関数を算出しておき、この補正関数の逆マトリックスを解くことにより、エネルギー分布の推定が可能となる。
【００９４】
つまり、各ｎ次のモーメントの計測値は、
x1 = a1[1:n]×E[n:1] (電流計測に相当)
x2 = a2[1:n]×E[n:1] (ＭＳＶ計測に相当)
x3 = a3[1:n]×E[n:1]
・・・
xn = an[1:n]×E[n:1]
xn: ｎ次モーメント値［スカラー量］
an: 応答マトリックス[１行ｎ列行列]
E: エネルギー分布[ｎ行１列行列]
となるから、行列の乗算として、
X[n:1] = A[n:n] × E[n:1]
と表すことができる。この行列Ａの逆行列Ａ−１を用いて、放射線のエネルギー分布は、
E[n:1] = A-1[n:n] × X[n:1]
として求めることができる。ただし、この１乗からｎ乗までのモーメント計測は、必要とするエネルギーバンド幅に相当する個数分を前記の式から選択すればよく、１次モーメントである平均電流値を除くことにより、すべて交流の計測手段のみで構成することができる。
【００９５】
以上のように、ＭＳＶ計測とパルス計測または電流計測を組み合せ、その比率から平均放射線エネルギーを推定し、線量に換算することができる。これは、従来のセンサ構造を用いて感度を補正していた手法と併用することにより、より線量応答に近い特性を容易に実現することができる。また、従来のパルス波高を求めて線量を算出していた手法に比べ、波高を選別する必要がなく、簡単な構成で実現が可能である。さらに、複数のｎ次モーメント値を用い、その応答関数から放射線のエネルギーを再構成することができ、パルス計測が困難な高計数率においても放射線のエネルギー分布の測定が可能となり、その情報からより正確な線量を評価することができる。かくして、これらを単独に、または従来の線量評価手法と併用することにより、より広いレンジの線量を、より正確に、一括して監視することの可能な放射線測定装置を提供することができる。
【００９６】
【発明の効果】
以上説明したように、本発明の放射線測定装置によれば、パルス波形のｎ乗値を算出し、その値でパルス識別を行うことによって、信号レベルの小さい両極性の回路ノイズおよびセンサのα線ノイズを低減し、センサパルスに対するこれらノイズの割合を低減することができる。これにより、従来は回路ノイズに紛れていたセンサ信号を測定することが可能となる。
【図面の簡単な説明】
【図１】 本発明による放射線測定装置の第１の実施の形態を示すブロック図である。
【図２】 （ａ）はＳＲＮＭセンサ出力パルスとサンプリングの例を示す説明図、（ｂ）はＳＲＮＭセンサ出力の２乗値を示す説明図である。
【図３】 （ａ）〜（ｃ）はｎ乗パルス計測手段の作用を説明するためのグラフである。
【図４】 本発明による放射線測定装置の第２の実施の形態を示すブロック図である。
【図５】 積分識別装置と差分識別装置の特徴を比較する図表である。
【図６】 本発明による放射線測定装置の第３の実施の形態を示すブロック図である。
【図７】 （ａ）〜（ｄ）は単極性変換手段の機能を説明するためのグラフである。
【図８】 本発明による放射線測定装置の第４の実施の形態を示すブロック図である。
【図９】 ノイズ波形例を示す波形図である。
【図１０】 第４の実施の形態において模擬中性子パルスを入力したときの第１平滑装置のシミュレーション結果と指標Ｘの推移を示すグラフである。
【図１１】 第４の実施の形態において第１平滑装置の出力の揺らぎ率Ｙを評価した結果を示すグラフである。
【図１２】 ＭＳＶ揺らぎとノイズの関係を説明するための波形図である。
【図１３】 ＭＳＶ指示値の揺らぎの範囲と平均値の変化幅の関係を示すグラフである。
【図１４】 （ａ）は最大揺らぎを超えるノイズが誘導された場合のパルスのＭＳＶ指示値の推移の一例を示すグラフ、（ｂ）は（ａ）を第４の実施の形態のデータ除去平均化手段によって処理した結果を示すグラフである。
【図１５】 本発明の参考例を示すブロック図である。
【図１６】 線量に対するＣｄＴｅセンサのパルス計数およびＭＳＶ指示値を示すグラフである。
【図１７】 線量に対するＣｄＴｅセンサのＭＳＶ指示値とパルス計数値の比率を示すグラフである。
【図１８】 従来のディジタル式原子炉起動監視装置のブロック図である。
【符号の説明】
１ ＳＲＮＭセンサ
２ アナログ式増幅器
２Ａ プリアンプ
２Ｂ プリアンプ
３ Ａ／Ｄ変換器
４ ｎ乗パルス識別装置
５ パルス計数装置（ＰＣ）
６ 積分識別装置
７ 差分識別装置
８ 波高・パワー識別装置
９ 単極性変換手段
１０ バンドパスフィルタ（ＢＰＦ）
１１ ｎ乗演算装置
１２Ａ 第１平滑装置
１２Ｂ 第２平滑装置
１３ データ除去平均化装置（ＤＥＡ）
１４ ＭＳＶ中性子評価装置
１５ ノイズ特性評価装置
１６ ＣｄＴｅセンサ
１７ チャージアンプ
１８ ＭＳＶ計測装置
１９ 電流計測装置
２０ パルス計数装置
２１ エネルギー評価装置
２２ 線量評価装置
２３ パルス計数手段
２４ 和演算手段
２５ パワー演算手段
２６ 加算平均手段
２７ 原子炉出力評価手段
Claims (9)
- 放射線に感応してパルス信号を出力する放射線検出手段と、この放射線検出手段の出力パルス信号を整形しサンプル値に変換するＡ／Ｄ変換手段と、前記サンプル値のｎ乗値相当（ｎ＞＝２）を演算し前記放射線のパルス信号のみを識別するｎ乗パルス識別手段と、このｎ乗パルス識別手段によって識別された放射線のパルス信号を計数する計数手段とを備え、しかも前記ｎ乗パルス識別手段は、前後ｎ個のサンプル値を掛け算することにより前記ｎ乗値相当を演算するか、又は同じサンプル値をｎ乗した後に前後２個のｎ乗値を平均化処理することにより前記ｎ乗値相当を演算するものであることを特徴とする放射線測定装置。
- 前記ｎ乗パルス識別手段は、前記放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、パルス立ち上がり部およびパルス立ち下がり部に対応する両サンプリング点におけるサンプル値の和のｎ乗値を、前記両サンプリング点の間に位置する両サンプリング点におけるサンプル値の和のｎ乗値から差し引くことにより前記出力パルスの波高に相当する差を演算し、その差から所望のパルスのみを識別する差分識別手段を備えていることを特徴とする請求項１記載の放射線測定装置。
- 前記ｎ乗パルス識別手段は、前記放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、前半２つのサンプリング点におけるサンプル値の差のｎ乗値と、後半２つのサンプリング点におけるサンプル値の差のｎ乗値との和を演算し、その和に基づいて所望のパルスのみを識別する差分識別手段を備えていることを特徴とする請求項１記載の放射線測定装置。
- 前記ｎ乗パルス識別手段は、前記放射線検出手段の出力パルスのパルス幅に相当する連続する４つのサンプリング点のサンプル値を用い、前記各サンプル値のｎ乗値の積算値と前記放射線検出手段の出力パルスの幅に相当するサンプル値の波高との比率を算出し、その比率から前記放射線検出手段の出力パルスを識別する波高・パワー識別手段を備え、前記計数手段は前記波高・パワー識別手段によって識別されたパルスを計数することを特徴とする請求項１記載の放射線測定装置。
- 正および負の両極性の入力信号を単極性の信号に変換する単極性変換手段をさらに具備し、かつ前記ｎ乗パルス識別手段は、前記単極性変換手段によって得られた単極性の信号を偶数ｎを用いてｎ乗演算を行うｎ乗演算手段と、このｎ乗演算手段の出力を平均化する平滑手段とを備えていることを特徴とする請求項１記載の放射線測定装置。
- 放射線に感応してパルス信号を出力する放射線検出手段と、この放射線検出手段の出力パルスを整形し周波数帯域を制限する帯域制限手段と、この帯域制限手段の各出力パルスのｎ乗値（ｎ＞＝２）を算出するｎ乗値演算手段と、このｎ乗値演算手段の出力信号を所定の時間幅で平滑化する第１の平滑化手段と、この第１の平滑化手段の出力の大小関係を評価し、その評価結果から特定のデータを除くデータ除去平均化手段と、このデータ除去平均化手段からの出力をさらに平滑化する第２の平滑化手段と、この第２の平滑化手段の出力を放射線強度に換算する変換手段とを具備した放射線測定装置。
- 前記第１の平滑化手段の平滑時間幅を外来ノイズのパルス幅以上に設定し、前記データ除去平均化手段の平滑化区間を外来ノイズパルスの到来する時間間隔以下に設定し、計測されたノイズ波形によりこれら平滑時間幅とデータ除去数を調整するノイズ特性評価手段をさらに具備したことを特徴とする請求項６記載の放射線測定装置。
- 前記データ除去平均化手段は、前記第１の平滑化手段の過去の出力と現在の出力とを比較し、現在の出力値が前記第１の平滑化手段の揺らぎから過去の出力値をもとに事前に評価したしきい値以上に増加していた場合、前記第１の平滑化手段の現在の出力を除去し、予め評価によって求めた評価比率を過去の出力値に乗じた値に置き換えるデータ制限手段を有することを特徴とする請求項６記載の放射線測定装置。
- 前記ｎは３以上の奇数であり、前記ｎ乗値演算手段は前記帯域制限手段の各出力パルスに対して奇数次のモーメントを算出することを特徴とする請求項６記載の放射線測定装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002046788A JP3958069B2 (ja) | 2001-03-28 | 2002-02-22 | 放射線測定装置 |
| US10/106,895 US6836523B2 (en) | 2001-03-28 | 2002-03-27 | Radiation measurement device |
| US10/844,350 US7139354B2 (en) | 2001-03-28 | 2004-05-13 | Radiation measurement device |
| US11/348,343 US7333585B2 (en) | 2001-03-28 | 2006-02-07 | Radiation measurement device |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001093306 | 2001-03-28 | ||
| JP2001-93306 | 2001-03-28 | ||
| JP2002046788A JP3958069B2 (ja) | 2001-03-28 | 2002-02-22 | 放射線測定装置 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007100635A Division JP4417972B2 (ja) | 2001-03-28 | 2007-04-06 | 放射線測定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002357692A JP2002357692A (ja) | 2002-12-13 |
| JP3958069B2 true JP3958069B2 (ja) | 2007-08-15 |
Family
ID=26612395
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002046788A Expired - Lifetime JP3958069B2 (ja) | 2001-03-28 | 2002-02-22 | 放射線測定装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (3) | US6836523B2 (ja) |
| JP (1) | JP3958069B2 (ja) |
Families Citing this family (44)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3741006B2 (ja) * | 2001-08-08 | 2006-02-01 | 株式会社日立製作所 | 荷電粒子測定装置およびその測定方法 |
| JP4334256B2 (ja) * | 2003-03-26 | 2009-09-30 | 株式会社東芝 | 放射線出力監視装置 |
| ATE526880T1 (de) * | 2003-05-14 | 2011-10-15 | Koninkl Philips Electronics Nv | Verfahren und apparat zur verbesserten strahlendetektion |
| US8358052B2 (en) * | 2003-06-30 | 2013-01-22 | Piedek Technical Laboratory | Unit having resonator, oscillator having unit and electronic apparatus having unit |
| JP4659492B2 (ja) * | 2005-03-18 | 2011-03-30 | 日立Ｇｅニュークリア・エナジー株式会社 | 起動領域モニタシステム用ノイズ監視システム |
| US7180055B1 (en) * | 2005-08-09 | 2007-02-20 | Sandia Corporation | Nuclear sensor signal processing circuit |
| US7609204B2 (en) * | 2005-08-30 | 2009-10-27 | Honeywell International Inc. | System and method for dynamically estimating output variances for carrier-smoothing filters |
| US7294840B2 (en) * | 2006-03-01 | 2007-11-13 | Nucsafe, Inc. | Apparatus and method for reducing microphonic susceptibility in a radiation detector |
| DE102007034982B4 (de) * | 2007-07-26 | 2016-07-07 | Siemens Healthcare Gmbh | Verfahren zum Betreiben eines getakteten, zählenden Röntgenstrahlendetektors |
| US8039787B2 (en) * | 2007-08-03 | 2011-10-18 | Pulsetor, Llc | Digital pulse processor slope correction |
| US7741609B2 (en) * | 2007-08-03 | 2010-06-22 | Pulsetor, Llc | Adapting a high-performance pulse processor to an existing spectrometry system |
| EP2183579B1 (en) * | 2007-08-03 | 2018-09-12 | Judith B. Mott Revocable Trust Dated June 6, 2006, As Amended And Restated | Pileup rejection in an energy-dispersive radiation spectrometry system |
| US8253942B2 (en) * | 2007-09-27 | 2012-08-28 | Scott Technologies, Inc. | Optical gas detector |
| US7807973B2 (en) * | 2008-08-01 | 2010-10-05 | Pulsetor, Llc | Pileup rejection in an energy-dispersive radiation spectrometry system |
| EP2410920A1 (en) * | 2009-03-26 | 2012-02-01 | Koninklijke Philips Electronics N.V. | Spectral imaging |
| WO2010119358A2 (en) * | 2009-04-16 | 2010-10-21 | Koninklijke Philips Electronics N.V. | Spectral imaging |
| JP5336934B2 (ja) * | 2009-06-05 | 2013-11-06 | 株式会社東芝 | 広域中性子束監視システムおよび検出器特性評価装置 |
| JP5419670B2 (ja) * | 2009-12-14 | 2014-02-19 | 三菱電機株式会社 | 放射線測定装置及びその診断方法 |
| CN102110439B (zh) * | 2009-12-24 | 2012-07-25 | 联咏科技股份有限公司 | 频谱噪声检测方法 |
| FR2956496B1 (fr) * | 2010-02-17 | 2012-03-09 | Commissariat Energie Atomique | Procede de mesure en ligne de rayonnements ionisants |
| US9081102B2 (en) * | 2010-05-21 | 2015-07-14 | Lewis Ronald Carroll | Apparatus for extending a scintillation detector's dynamic range |
| CN102455430A (zh) * | 2010-10-27 | 2012-05-16 | 成都理工大学 | 一种核脉冲信号的数字化成形方法 |
| JP5535100B2 (ja) * | 2011-02-03 | 2014-07-02 | 三菱電機株式会社 | 炉外核計装装置 |
| US8930155B2 (en) | 2011-11-15 | 2015-01-06 | Fuji Electrict Co., Ltd. | Pulse processing device and radiation measuring device |
| CN102540240B (zh) * | 2012-01-12 | 2014-05-07 | 深圳大学 | 一种脉冲幅度分析电路及脉冲幅度分析器 |
| JP5773949B2 (ja) * | 2012-06-06 | 2015-09-02 | 三菱電機株式会社 | 放射線監視装置 |
| CN103995279B (zh) * | 2013-02-18 | 2017-03-22 | 中国科学技术大学 | 用于宽动态范围中子通量监测的动态标定***及方法 |
| JP5753551B2 (ja) * | 2013-04-25 | 2015-07-22 | 日立アロカメディカル株式会社 | 放射線測定装置 |
| CN103705258B (zh) * | 2013-11-19 | 2015-06-24 | 江苏康众数字医疗设备有限公司 | 一种成像设备的自动曝光控制方法及装置 |
| JP6080778B2 (ja) * | 2014-01-14 | 2017-02-15 | 三菱電機株式会社 | 放射線監視装置 |
| JP6325343B2 (ja) * | 2014-05-26 | 2018-05-16 | 株式会社日立製作所 | 放射線測定装置 |
| US9847215B2 (en) * | 2014-11-08 | 2017-12-19 | Jefferson Science Associates, Llc | Method for detecting and distinguishing between specific types of environmental radiation using a high pressure ionization chamber with pulse-mode readout |
| JP6522926B2 (ja) * | 2014-11-11 | 2019-05-29 | 株式会社東芝 | 中性子測定装置および中性子測定方法 |
| JP6441062B2 (ja) * | 2014-12-17 | 2018-12-19 | 株式会社東芝 | 中性子測定装置および中性子測定方法 |
| KR101733092B1 (ko) | 2015-01-07 | 2017-05-08 | 한국해양과학기술원 | 휴대용 방사능 검출기의 시간 영역에서의 최고값 검출 장치 및 그 방법 |
| JP6628701B2 (ja) * | 2016-08-05 | 2020-01-15 | 三菱電機株式会社 | 放射線測定装置 |
| CN109791211B (zh) | 2016-10-21 | 2020-08-25 | 珀金埃尔默健康科学公司 | 用改进的事件类型区别进行辐射检测的***及方法 |
| JP6748588B2 (ja) * | 2017-02-01 | 2020-09-02 | キヤノン株式会社 | 放射線撮影装置、放射線撮影方法およびプログラム |
| JP2019052854A (ja) * | 2017-09-12 | 2019-04-04 | 株式会社東芝 | 放射線検出装置および放射線検出方法 |
| KR20200054231A (ko) | 2017-09-15 | 2020-05-19 | 퍼킨엘머 헬스 사이언시즈, 아이엔씨. | 전자 테스트 소스를 사용하여 신틸레이션 이벤트들을 에뮬레이팅하기 위한 시스템들 및 방법들 |
| EP3567404A1 (en) | 2018-05-09 | 2019-11-13 | Target Systemelektronik GmbH & Co. KG | Method and device for the measurement of high dose rates of ionizing radiation |
| CN110276441B (zh) * | 2019-05-20 | 2021-05-04 | 成都理工大学 | 一种基于深度学习的梯形重叠核脉冲估计方法 |
| CN110161548A (zh) * | 2019-06-27 | 2019-08-23 | 南华大学 | 一种剂量率探测的方法、装置以及剂量率探测器 |
| CN112906851B (zh) * | 2021-03-24 | 2023-09-19 | 中国兵器装备集团自动化研究所有限公司 | 一种基于分段比例跳变阈值判断的计数方法 |
Family Cites Families (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NL278375A (ja) * | 1961-05-12 | |||
| US3894912A (en) * | 1972-06-21 | 1975-07-15 | Us Energy | Determination of parameters of a nuclear reactor through noise measurements |
| DE2634590A1 (de) * | 1976-07-31 | 1978-02-02 | Hartmann & Braun Ag | Verfahren zur kontinuierlichen funktionsueberwachung von kernstrahlungsmesskanaelen |
| JPS5921111A (ja) * | 1982-07-26 | 1984-02-03 | Toshiba Corp | ワイドレンジモニタ装置 |
| US4623508A (en) * | 1984-02-15 | 1986-11-18 | Reuter-Stokes, Inc. | Wide range flux monitor assembly |
| US4920548A (en) | 1988-09-28 | 1990-04-24 | Westinghouse Electric Corp. | Source range neutron flux count rate system incorporating method and apparatus for eliminating noise from pulse signal |
| US4931646A (en) * | 1989-03-17 | 1990-06-05 | The United States Of America As Represented By The Secretary Of The Army | Remote multichannel coincident nuclear detector and characterization system |
| US5067090A (en) * | 1989-05-01 | 1991-11-19 | Schlumberger Technology Corporation | Nuclear spectroscopy method and apparatus for digital pulse height analysis |
| JP2580811B2 (ja) | 1989-12-13 | 1997-02-12 | 富士電機株式会社 | １▲ｃｍ▼深部線量当量測定用放射線測定器 |
| JP2698680B2 (ja) * | 1990-02-01 | 1998-01-19 | 株式会社東芝 | 原子炉出力監視装置 |
| US5164895A (en) * | 1990-03-16 | 1992-11-17 | Westinghouse Electric Corp. | Neutron flux mapping system for nuclear reactors |
| US5076998A (en) * | 1990-06-21 | 1991-12-31 | Westinghouse Electric Corp. | Monitoring of low frequency pulse rate |
| US5114665A (en) * | 1990-07-25 | 1992-05-19 | The Babcock & Wilcox Company | Autonormalizing reactimeter |
| JPH0545490A (ja) * | 1991-08-21 | 1993-02-23 | Toshiba Corp | 原子炉出力監視装置 |
| JP2877609B2 (ja) * | 1992-03-30 | 1999-03-31 | 株式会社東芝 | 原子力発電設備用起動領域モニタシステム |
| JPH06112772A (ja) | 1992-09-25 | 1994-04-22 | Hitachi Ltd | フィルタ演算方法、フィルタ演算装置、サージ除去装置 |
| JP3183983B2 (ja) | 1993-02-23 | 2001-07-09 | ミサワホーム株式会社 | 建物ユニットの骨組み構造 |
| JP3226698B2 (ja) | 1994-02-24 | 2001-11-05 | 東芝エレベータ株式会社 | エレベータの防犯装置 |
| US5532944A (en) * | 1994-07-28 | 1996-07-02 | Sorrento Electronics, Inc. | Multi-channel analysis system and method using digital signal processing |
| US6194726B1 (en) * | 1994-12-23 | 2001-02-27 | Digirad Corporation | Semiconductor radiation detector with downconversion element |
| JP3567045B2 (ja) * | 1996-04-03 | 2004-09-15 | 株式会社東芝 | 原子炉出力監視装置 |
| JP3517331B2 (ja) | 1997-04-07 | 2004-04-12 | 三菱電機株式会社 | ディジタル計数率計 |
| CN100555470C (zh) * | 1997-12-01 | 2009-10-28 | 东芝株式会社 | 原子反应堆输出功率的监控装置 |
| JP4299927B2 (ja) * | 1998-08-31 | 2009-07-22 | 株式会社東芝 | 中性子束計測装置 |
| JP3936107B2 (ja) | 1998-09-21 | 2007-06-27 | 株式会社東芝 | 原子炉起動監視装置および炉内出力監視方法 |
| JP2001021654A (ja) | 1999-07-09 | 2001-01-26 | Mitsubishi Electric Corp | 放射線計測システム |
-
2002
- 2002-02-22 JP JP2002046788A patent/JP3958069B2/ja not_active Expired - Lifetime
- 2002-03-27 US US10/106,895 patent/US6836523B2/en not_active Expired - Lifetime
-
2004
- 2004-05-13 US US10/844,350 patent/US7139354B2/en not_active Expired - Lifetime
-
2006
- 2006-02-07 US US11/348,343 patent/US7333585B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US20040206909A1 (en) | 2004-10-21 |
| JP2002357692A (ja) | 2002-12-13 |
| US7139354B2 (en) | 2006-11-21 |
| US7333585B2 (en) | 2008-02-19 |
| US20020146087A1 (en) | 2002-10-10 |
| US20060126776A1 (en) | 2006-06-15 |
| US6836523B2 (en) | 2004-12-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3958069B2 (ja) | 放射線測定装置 | |
| RU2414724C2 (ru) | Способ и устройство для спектральной компьютерной томографии | |
| JP5171891B2 (ja) | 放射線測定装置 | |
| EP3637147A1 (en) | Gain correction apparatus and method for scintillation detector | |
| US9029769B2 (en) | Dose rate measuring apparatus | |
| JP4417972B2 (ja) | 放射線測定装置 | |
| JP2004108796A (ja) | 放射線測定装置 | |
| JP2003513250A (ja) | 差分訂正法と差分訂正装置 | |
| JP3807652B2 (ja) | 放射線測定装置及び方法 | |
| CN108196293B (zh) | 一种基于闪烁体探测器剂量率检测方法 | |
| US10353081B2 (en) | Gamma system count loss correction with virtual pulse injection | |
| EP3719541A1 (en) | Guard efficiency compensation system | |
| JP2001194460A (ja) | 放射線モニタ | |
| JP2001116844A (ja) | 放射線測定装置 | |
| JP7183206B2 (ja) | 放射能検査装置 | |
| JP2005049144A (ja) | 放射線計測方法 | |
| RU2729600C1 (ru) | Способ диагностирования стабильности работы устройства с коронным счетчиком для измерения нейтронных потоков в присутствии гамма-излучения | |
| JP2019120656A (ja) | 放射線計測装置及び放射線計測方法 | |
| JP2004138511A (ja) | 中性子測定装置 | |
| JPH0921880A (ja) | 放射線測定装置 | |
| JP3728220B2 (ja) | 比例計数管型中性子検出器のγ線感度試験方法 | |
| JP3563922B2 (ja) | 放射線測定装置 | |
| US10330797B2 (en) | Method for dead time determination in a gamma camera and a system for accomplishing the same | |
| JPH03189586A (ja) | 放射能測定装置 | |
| JP2020148689A (ja) | 検出装置及びスペクトル生成装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040618 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070209 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070406 |
|
| 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: 20070427 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070509 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 3958069 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110518 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110518 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120518 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120518 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130518 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130518 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140518 Year of fee payment: 7 |
|
| EXPY | Cancellation because of completion of term |