JP4761829B2 - 軸方向ボイド率分布測定方法および収納装置収納前の燃料集合体中性子増倍率評価方法 - Google Patents

軸方向ボイド率分布測定方法および収納装置収納前の燃料集合体中性子増倍率評価方法 Download PDF

Info

Publication number
JP4761829B2
JP4761829B2 JP2005143884A JP2005143884A JP4761829B2 JP 4761829 B2 JP4761829 B2 JP 4761829B2 JP 2005143884 A JP2005143884 A JP 2005143884A JP 2005143884 A JP2005143884 A JP 2005143884A JP 4761829 B2 JP4761829 B2 JP 4761829B2
Authority
JP
Japan
Prior art keywords
distribution
axial
neutron
void ratio
burnup
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 - Fee Related
Application number
JP2005143884A
Other languages
English (en)
Other versions
JP2006322727A (ja
Inventor
精 植田
宏徳 熊埜御堂
偉司 三橋
司 菊池
研一 吉岡
智治 佐々木
宏司 平岩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP2005143884A priority Critical patent/JP4761829B2/ja
Priority to US11/434,195 priority patent/US20070076839A1/en
Priority to TW095117333A priority patent/TWI320189B/zh
Publication of JP2006322727A publication Critical patent/JP2006322727A/ja
Priority to US12/345,438 priority patent/US8401141B2/en
Application granted granted Critical
Publication of JP4761829B2 publication Critical patent/JP4761829B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21CNUCLEAR REACTORS
    • G21C1/00Reactor types
    • G21C1/04Thermal reactors ; Epithermal reactors
    • G21C1/06Heterogeneous reactors, i.e. in which fuel and moderator are separated
    • G21C1/08Heterogeneous reactors, i.e. in which fuel and moderator are separated moderator being highly pressurised, e.g. boiling water reactor, integral super-heat reactor, pressurised water reactor
    • G21C1/084Boiling water reactors
    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21CNUCLEAR REACTORS
    • G21C17/00Monitoring; Testing ; Maintaining
    • G21C17/10Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Plasma & Fusion (AREA)
  • General Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Description

本発明は、沸騰水型原子炉の炉心内の軸方向ボイド率分布測定方法および燃料集合体を収納装置に収納する前の燃料集合体中性子増倍率評価方法に関する。
沸騰水型原子炉(以下BWR、または単に原子炉という)では冷却水が沸騰し、軸方向、すなわち冷却水の流れの上流側から下流側(炉心の上部、燃料集合体の上部)に向かって、チャンネルボックス内の冷却水が通過し得る空間に占める気泡(ボイド)の割合(以下ボイド率という)が増大する。ボイド率の軸方向分布を、ここでは簡単にボイド率分布ともいうことにする。ボイド率はBWRにおいては、炉心の反応度や出力分布などの核的特性や冷却特性に大きな影響を与えるため、ボイド率分布を定量的に明らかにすることは極めて重要なテーマである。しかしながら、実用可能な測定技術が存在しないため、実際の商業用BWR炉心の中で測定された例は見当たらない。従来は炉心の外でモックアップ実験を行ない、その結果に基づいて理論モデルを構築し、そのモデルを用いて実機BWRのボイド率分布の設計が成されてきた。
ボイド率分布測定の必要性は古くから認識されており、実質的に発熱を伴わない臨界実験装置においては、模擬ボイドを用いてボイド率による核特性の変化を測定した例がある。この方法では、マンガン(Mn)を含むワイヤ(マンガニン線)を実験炉心内に挿入し、その一部に熱中性子吸収体であるカドミウム(Cd)を巻き付け、Mn55の(n,γ)反応における反応率のカドミ比(Cd比)を測定している。実験炉心の中にボイドを模擬するアルミニウム中空管を挿入し、中空管の大きさを変えてボイド率を変え、Cd比がボイド率と良好な相関関係にあることから、Cd比を測定すればボイド率を測定できるとしている。この方法の基本原理は、熱外中性子束と熱中性子束との比がボイド率と良好な関係にあることに着目した点で優れている。
本発明者の一人は前述の基本原理を用いながらCdを用いる必要の無い、より実用的な方法を、特許文献1および特許文献2において提案した。これらの文献では、燃料集合体内部、炉内計装管、固定型あるいは移動型の中性子検出器に融点の高いGd2O3 のような強い熱中性子吸収体、あるいはステンレス鋼材などのような弱い中性子吸収材を局所的に配置し、熱中性子束の歪みを局所的に引き起こすことによって熱中性子束と熱外中性子束とを分離し、両者の比からボイド率を求める方法である。
特開昭55-121195号公報 特開昭55-125489号公報 特開昭61-262689号公報
前述の方法では融点が低く、また実機原子炉構造材に悪影響を与える恐れがあるCdを使用する点で実機への応用は困難であった。特許文献1および特許文献2では、熱中性子束の歪みを人為的に引き起こさせる点で、いまだに実機への適用の見込みはたっていない。
本発明者らは最近、炉内計装管を活用する新たな方法を検討している。そこでは、熱外中性子束以上にボイド率に依存する度合いが高い高速中性子束と熱中性子束との比にも着目して、この比および熱外中性子束と熱中性子束との比を利用してボイド率を測定する案である。この方法では原子炉の運転を停止して行なわれる定期検査の際に特別に炉内計装管を取り出して計測する必要があった。計測に当たっては細くて正確な形状を保持することが容易でない計装管と検出器とを正確に位置決めする専用の保持具を準備する必要があった。さらに計装管は燃料集合体相互間の水ギャップに設置されているため、少なくとも計装管を取り囲む4体の燃料集合体の平均的なボイド率に対応し、かつ燃料集合体内部より検出感度が大幅に低下する位置であるという問題が残っていた。
以上の議論のように軸方向ボイド率分布はBWR炉心において非常に重要な要素であるにもかかわらず、いまだに実機商用炉における実測は実現していない。
ところで、照射済BWR燃料集合体を原子炉から取り出して水中で輸送容器(キャスク)や燃料貯蔵ラックのような燃料集合体収納装置に収納する場合を考えると、軸方向燃料有効部のうち、下端から2/3ないし3/4付近で、無限増倍率(あるいは反応度)が高くなる部位が発生する。出力運転中の軸方向出力分布を平坦にする、あるいは軸方向出力分布が炉心最大出力にならないようにする、などの設計条件と、軸方向ボイド率分布の影響により、上部(冷却材下流側)においてウランや生成したプルトニウムの燃焼遅れがあることと高い転換比によるプルトニウムの高い生成割合に起因して無限増倍率(あるいは反応度)が高くなる部位が発生するものである。
このような背景から、水が満たされた状態の燃料集合体収納装置においては中性子増倍率が高く、したがって未臨界度が浅くなる部位は、上述のように、軸方向燃料有効部のうち、下部から2/3ないし3/4付近に発生する。この部位は設計や運転条件によってほぼ決まる。したがって中性子増倍率が厳しくなる背景である軸方向ボイド率分布の影響を考慮しないと、安全な未臨界度評価、すなわち臨界安全性確保ができないことになる。従来は、設計余裕を過剰に取ることによりこの問題を回避してきた。
本発明は上述の課題を解決するために成されたものであり、第1の目的は軸方向ボイド率分布を測定する実現性の高い新しい方法を提供することであり、第2の目的は照射済(または使用済)BWR燃料集合体を燃料集合体収納装置に収納する際に、過度の設計裕度を取らなくても良いように、未臨界度が浅くなる(中性子増倍率が高く臨界に近付く)部位における燃料集合体の中性子増倍率を実験的に評価することができる方法を提供することである。
上記目的を達成するため、請求項1の発明は、原子炉内で中性子に照射された照射燃料集合体において、重核種が中性子捕獲を行って生成しかつ中性子を放出する放射性核種に起因する第1群の放射能と、核分裂で生成した後、中性子捕獲反応を受けておらずかつ所定の範囲で生成量が燃焼度に比例しかつガンマ線を放出する核分裂生成物に起因する第2群の放射能とを測定し、測定された第1群の放射能の強度である中性子放出率Aを燃焼度Eの指数(α)乗に比例するものとして、比例定数をaとして、A=a・Eαで表すとともに、測定される第2群の放射能の強度であるガンマ線強度Bを燃焼度Eに比例するものとして、比例定数をbとして、B=b・Eで表し、第2群の放射能強度の式の両辺をα乗して第1群の式との比をとって燃焼度の因子を消去して得られる式を、「0」を付して示す基準場所における値との比として求めた(式1)、すなわち
(式1): a/a0 =(A/A0)・(B0/B)α・(b/b0α
において、(b/b0)は1とみなし、(A/A0 )および(B0 /B)の測定値、および計算によって求めた指数αから(a/a0 )を求め、(a/a0 )が原子炉内部の軸方向ボイド率分布に依存する特性を利用して、別途求められた「a/a 0 とボイド率との相関関係」を用いて軸方向ボイド率分布を求めること、を特徴とする軸方向ボイド率分布測定方法である。
また、請求項の発明は、中性子放出率はキュリウム242を除く核種から放出される中性子放出率であり、核分裂生成物はセシウム137であることを特徴とする請求項記載の軸方向ボイド率分布測定方法である。
また、請求項の発明は、中性子放出率はキュリウム244から放出される中性子放出率であり、核分裂生成物はセシウム137またはセリウム144であることを特徴とする請求項記載の軸方向ボイド率分布測定方法である。
また、請求項の発明は、中性子放出率はキュリウム242から放出される中性子放出率であり、核分裂生成物はセシウム137またはセリウム144であることを特徴とする請求項記載の軸方向ボイド率分布測定方法である。
また、請求項の発明は、中性子放出率は、キュリウム242およびキュリウム244から放出される中性子放出率を除く中性子放出率であり、核分裂生成物はセシウム137またはセリウム144であることを特徴とする請求項記載の軸方向ボイド率分布測定方法である。
また、請求項6の発明は、原子炉内で中性子に照射された照射燃料集合体において、核分裂生成物がさらに中性子捕獲して生成しかつガンマ線を放出する放射性核種に起因する第1群の放射能と、核分裂で生成した後、中性子捕獲反応を受けておらずかつ所定の範囲で生成量が燃焼度に比例しかつガンマ線を放出する核分裂生成物に起因する第2群の放射能とを測定し、測定された第1群の放射能の強度であるガンマ線強度Aを燃焼度Eの指数(α)乗に比例するものとして、比例定数をaとして、A=a・Eαで表すとともに、測定される第2群の放射能の強度であるガンマ線強度Bを燃焼度Eに比例するものとして、比例定数をbとして、B=b・Eで表し、第2群の放射能強度の式の両辺をα乗して第1群の式との比をとって燃焼度の因子を消去して得られる式を、「0」を付して示す基準場所における値との比として求めた(式1)、すなわち
(式1): a/a0 =(A/A0)・(B0/B)α・(b/b0α
において、(b/b0)は1とみなし、(A/A0 )および(B0 /B)の測定値、および計算によって求めた指数αから(a/a0 )を求め、(a/a0 )が原子炉内部の軸方向ボイド率分布に依存する特性を利用して、別途求められた「a/a0とボイド率との相関関係」を用いて軸方向ボイド率分布を求めること、を特徴とする軸方向ボイド率分布測定方法である。
また、請求項の発明は、第1群の放射能を有する核種がセシウム134またはユーロピウム154であり、第2群の放射能を有する核種がセシウム137またはセリウム144であることを特徴とする請求項記載の軸方向ボイド率分布測定方法である。
また、請求項8の発明は、原子炉内で中性子に照射された照射燃料集合体から放出されている、重核種が中性子捕獲を行って生成する放射性核種に起因する中性子放出率を、所定の冷却時間を挟んで少なくとも2回測定し、半減期の違いを利用してキュリウム242から放出される第1の中性子放出率とその他核種から放出される第2の中性子放出率とに分類し、測定された第1の中性子放出率Aを燃焼度Eの指数(α)乗に比例するものとして、A=a・Eα(ただしaは比例定数)で表すとともに、測定される第2の中性子放出率Bを燃焼度Eの指数(β)乗に比例するものとして、比例定数をbとして、B=b・Eβで表し、「0」を付して示す基準場所における値との比として求めた(式2)、すなわち
(式2): (a/a0β・(b0/b)α=(A/A0β・(B0/B)α
において、A/A0 およびB0 /Bの測定値および計算によって求めた指数αとβから、(a/a0β・(b0/b)αを求め、別途求められた「(a/a0β・(b0/b)αとボイド率との相関関係」を用いて軸方向ボイド率分布を求めること、を特徴とする軸方向ボイド率分布測定方法である。
また、請求項9の発明は、ガンマ線強度γの軸方向分布と中性子束または中性子計数率φの軸方向分布を測定してそれぞれ基準場所における値γ0 およびφ0 に対する相対値γ/γ0 およびφ/φ0 を求め、一方では計算によって求めた指数αを用いて(γ0/γ )αを、他方では(式3)、すなわち
(式3): S/S0=(φ/φ0)・(1−k)/(1−k0
ただし、(S/S0 )は基準場所における値に対して規格化した中性子放出率の軸方向分布、k0基準場所における中性子増倍率、kはφ測定点における中性子増倍率
で表し、中性子増倍率分布(kおよびk0 )に計算値を用いて(式3)から(S/S0 )を求め、(式1)に対応する(式4)すなわち
(式4): a/a0=(S/S0)・(γ0/γ)α
によりa/a0 を算出し、別途計算で求めた「a/a0 とボイド率との相関」を示す校正曲線を用いて軸方向ボイド率分布を求めること、を特徴とする請求項1乃至5のいずれか記載の軸方向ボイド率分布測定方法である。
また、請求項10の発明は、軸方向ガンマ線分布測定値の平均値をこの申告値に規格化することによって軸方向燃焼度分布(E、E0 )を求め、一方、無限増倍率k を、少なくともボイド率に依存する係数を用いた燃焼度Eの1次式または2次式で表し、前記軸方向燃焼度分布と仮定した軸方向ボイド率分布とを用いて軸方向無限増倍率k 分布を求め、測定条件に対応する中性子輸送拡散計算により換算因子Fk を求めて、k 分布にFk をかけることにより軸方向中性子増倍率k分布へ換算し、一方、軸方向中性子束分布から(φ/φ0)を求め、これらと(式3)すなわち
(式3): S/S0=(φ/φ0)・(1−k)/(1−k0
ただし、(S/S0 )は基準場所における値に対して規格化した中性子放出率の軸方向分布、k0基準場所における中性子増倍率、kはφ測定点における中性子増倍率
とからS/S0 を求め、一方では、計算によって求めた指数αと軸方向ガンマ線分布から(γ0 /γ)αを求めて(式4)すなわち
(式4): a/a0=(S/S0)・(γ0/γ)α
により(a/a0 )を算出し、別途計算で求めた「a/a0 とボイド率との相関関係」を用いて軸方向ボイド率分布を求め、この分布が前記軸方向ボイド率分布仮定値または繰り返し計算における直前に求められた軸方向ボイド率分布と指定された範囲で一致していれば最終結果とし、一致しない場合には初期値として与えた軸方向ボイド率分布仮定値あるいは直前の繰り返し前のボイド率分布と置換することによって繰り返し計算を行ない、最終結果が得られて収束した値を軸方向ボイド率分布とすること、を特徴とする請求項1乃至5のいずれか記載の軸方向ボイド率分布測定方法である。
また、請求項11の発明は、基準場所における燃焼度E0 を中性子放出率法によって求め、軸方向ガンマ線分布測定値のうち基準場所における値γ0 を前記燃焼度E0 に規格化することによって軸方向燃焼度分布(E、E0 )を求め、一方、無限増倍率k を少なくともボイド率に依存する係数を用いた燃焼度Eの1次式または2次式で表し、前記燃焼度分布と仮定した軸方向ボイド率分布とを用いて軸方向無限増倍率k 分布を求め、測定条件に対応する中性子輸送拡散計算により換算因子Fk を求めてk 分布にFk をかけることにより軸方向中性子増倍率(k)分布へ換算し、一方、軸方向中性子束分布から(φ/φ0)を求め、これらと(式3)すなわち
(式3): S/S0=(φ/φ0)・(1−k)/(1−k0
ただし、S/S0 基準場所における値に対して規格化した中性子放出率の軸方向分布、k0基準場所における中性子増倍率、kはφ測定点における中性子増倍率
とからS/S0 を求め、
一方では、計算によって求めた係数αと軸方向ガンマ線分布から(γ0 /γ)αを求めて、(式4)すなわち
(式4): a/a0=(S/S0)・(γ0/γ)α
によりa/a0 を算出し、別途計算で求めた「a/a0 とボイド率との相関関係」を用いて軸方向ボイド率分布を求め、この分布が前記軸方向ボイド率分布仮定値または繰り返し計算における直前に求められた軸方向ボイド率分布と指定された範囲で一致していれば最終結果とし、一致しない場合には初期値として与えた軸方向ボイド率分布仮定値あるいは直前の繰り返し前のボイド率分布と置換することによって繰り返し計算を行ない、最終結果が得られて収束した値を軸方向ボイド率分布とすること、を特徴とする請求項1乃至5のいずれか記載の軸方向ボイド率分布測定方法である。
また、請求項12の発明は、燃料集合体の軸方向位置を変えながらガンマ線スペクトルを測定し、それぞれの高さ位置において セシウム134または ユーロピウム154と セシウム137または セリウム144のフォトピークを測定してガンマ線強度比G、すなわち Cs134/Cs137 比、 Eu154/Cs137 比、または Cs134/Ce144 比の 0 」を付して示す基準場所における値で規格化した軸方向分布G/G0 を求め、G/G0 値とボイド率との相関を示す校正曲線を用いて軸方向ボイド率分布を求めること、を特徴とする請求項6または7記載の軸方向ボイド率分布測定方法である。
また、請求項13の発明は、請求項1ないし12のいずれか一項に記載の軸方向ボイド率分布測定方法によって軸方向ボイド率分布を求め、照射燃料集合体の基準場所と未臨界度モニタ点に対して中性子測定を行なって、基準場所における値に「 0 」を付して表すものとして、基準場所における中性子増倍率k0 および中性子束または中性子計数率の比φ0 /φを求め、また、グロスガンマ線強度比[(γg )/(γg0 ]を測定し、別途求めたグロスガンマ線強度相対分布と燃焼度相対分布との相関関係(計算値)を用いてグロスガンマ線強度比[(γg )/(γg0 ]から燃焼度相対比(E/E0)へ換算し、計算で求めた指数αを用いて(E/E0αを求め、前記軸方向ボイド率分布と、「a/a0 とボイド率との相関関係」を用いてa/a0 の値を求め、(式5)、すなわち
(式5): k=1−(1−k0)・(φ0/φ)・(E/E0α・(a/a0
により未臨界度モニタ点の中性子増倍率kを求めること、を特徴とする燃料集合体収納装置への収納前の燃料集合体中性子増倍率評価方法である。
本発明によれば従来の測定技術を活用して燃料集合体自体の測定を行なうことにより、測定のための新たな治具を原子炉の中に設置することなく、実機沸騰水型原子炉における軸方向ボイド率分布を測定することができる。この結果、炉心に大きな影響を与えるボイドの影響が明らかになり、炉心および燃料設計法の信頼度を評価することができる。さらに、照射(使用済)BWR燃料集合体を収納装置に収納する際の臨界安全裕度を設計において過剰に取り入れる必要がなくなり、収納装置の合理化設計に寄与し、製造コスト低減を図ることができる。
本発明は、原子炉の中で中性子に照射された照射BWR燃料集合体(取り出し燃料および再装荷される燃料集合体)から放出されている放射線のうち、ウランやプルトニウムが中性子を捕獲して生成した中性子放出核種(放射能)からの中性子放出率、核分裂生成物(放射能)が放出するガンマ線強度、および核分裂生成物がさら中性子を捕獲して生成した核種(放射能)が放出するガンマ線強度に着目した新規な方法により、照射BWR燃料集合体の軸方向ボイド率(割合)分布を測定する実用性の高い方法であり、さらに軸方向ボイド率(割合)分布を取り入れた、燃料集合体の軸方向中性子増倍率を評価する方法である。
上述の特性を理解するために、先ず本発明が着目する現象について詳細に説明する。
<<中性子放出率測定対象核種の特性>>
中性子で照射され、原子炉から取り出された原子炉用燃料から自発的に放出される中性子は、冷却時間が例えば2ケ月以内というように特に短い場合には核分裂生成物のLa140 から放出される高いエネルギーのγ線が水素原子と光核反応を起こして放出する中性子が存在するが、本発明ではそれ以降を対象とする。その場合、半減期と中性子放出率の大きさを考慮して、半減期が 162.8日と短いキュリウム242(Cm242)、半減期が18.1年と比較的長いCm244 、および実用的にはほとんど減衰しない「その他核種」(Residuals) に分類するのが便利であり、本発明でもこの方式を踏襲する。
Cm242からの中性子放出率を(S2 )、 Cm244からの中性子放出率を(S4 )、「その他核種」からの中性子放出率を(SR )と記載することにする。(SR )に属する核種として燃焼計算コードORIGENでは、U234、U238、Pu238 、Pu239 、Pu240 、Pu241 、Pu242 、Am241 、Am243 、Cm243 、Cm246 、および Cf252を取り上げているが、通常実質的には Pu238から Am241の範囲が対象として差し支えない。
次に、図9を用いて上述の中性子放出核種の生成過程をウラン燃料の場合について説明する。σc は中性子エネルギーが速度換算で 2200m/sの熱中性子に対する中性子捕獲断面積(核分裂断面積を除く)であり、Ic は共鳴領域の捕獲断面積であって、無限希釈における捕獲の共鳴積分と呼ばれている。減速材のボイド率(割合)が高くなると中性子が減速されにくくなり熱中性子の割合が低下するので、その結果熱中性子捕獲率に対して相対的に共鳴領域の中性子捕獲率が高くなる。ここでは「CR=Ic /σc : capture-ratio 」を定義する。この比が大きいとその反応がボイド率(割合)に対して感度が高いことを示すと考えられる。参考値として共鳴領域に共鳴現象が生じないボロンB10 の場合の値(CR=0.449)も記した。
燃料に含まれる大量のU238は中性子捕獲に関するCR比は著しく大きく、したがってボイド率に大きく依存して捕獲反応が進みU239が生成する。U238の実質的なCR比は自己遮蔽現象によって大幅に低下するものの、大きいことに変りはない。U239は半減期23.5分でNp239に、そしてNp239は2.35日の半減期でPu239 へ崩壊する。Pu239 のCR値はU238を初め他の核種より大幅に小さく、B10 の2倍にも及ばないので、Pu239 が中性子を捕獲する場合のボイド率依存性は小さい。Pu239 は一方では核分裂(Fis )を起こすが、他方では中性子を捕獲してPu240 が生成する。
Pu240 は比較的大きなCR値を有するのでボイド率依存性は比較的高い。Pu240 は自発核分裂で中性子(自発核分裂中性子)を放出しその放出率(NER) は比較的高いが、一方では中性子を捕獲してPu241 が生成する。Pu241 は一方では核分裂し、他方では捕獲反応によりPu242 となるがCR値が小さいためボイド率依存性は小さい。また半減期14.4年で崩壊してAm241 となるものもある。Am241 はα崩壊し、α線は燃料に含まれている酸素と(α,n)反応を起こして中性子((α,n)中性子)を放出する。Am241 は一方では比較的弱いながらボイド率に依存した中性子捕獲反応によりAm242 となり、Am242 は半減期16時間で崩壊してCm242 となる。Cm242 は自発核分裂中性子と(α,n)中性子の放出率(S2 )が非常に高い。Cm242 はα崩壊すると(α,n)中性子放出率が比較的高いPu238 となる。
一方、Pu241 は低いボイド率依存の中性子捕獲反応によりPu242 となる。Pu242 は高いボイド率依存の中性子捕獲反応によりPu243 となり、5時間の半減期でAm243 となる。Am243 は高いボイド率依存の中性子捕獲反応によりAm244 となり、半減期10時間で崩壊してCm244 となる。Cm244 は半減期18.1年でα崩壊してPu240 に戻るが、Cm244 は主に自発核分裂による中性子放出率(S4 )が非常に高い。
以上の説明からわかるように、図9において太枠で示す核種(U238、Pu240、Pu242、Am243)の中性子捕獲反応はボイド率に依存し、その結果生成する核種の中性子放出率は照射中のボイド率の影響を受けている。Cm242 からの中性子放出率(S2 )は主にU238とPu240 の中性子捕獲反応のボイド率依存性によって、Cm244 からの中性子放出率(S4 )はさらにPu242 とAm243 の中性子捕獲反応のボイド率依存性によって特性が規定されることがわかる。(SR )に属するものは、主にU238とPu240 の中性子捕獲反応のボイド率依存性によって特性が規定されることがわかる。
<<γ線測定対象の核種の特性>>
セシウム137(Cs137)は核分裂した後、中性子捕獲反応を起こすこと無く直接生成し、半減期30.2年で崩壊する。その際662keVのγ線を放出するのでこれを測定する。Cs137 の生成量は広い範囲の燃焼度にわたって燃焼度に比例する。冷却時間が半年以内とか燃焼度が低い場合には、Cs134 、Zr95-Nb95 、その他が放出するγ線による非常に高いバックグランドに邪魔されて測定精度が低下する恐れがあるため、そのような場合にはCs137 の代わりにセリウム144(Ce144)などを用いる方が好適な場合がある。 Cs137生成時のボイド率依存はほとんど無い。核分裂核種による収率の変化も概して小さい。
Cs134 は核分裂で生成した安定核種の Cs133が中性子を捕獲して生成し、半減期2.06年で崩壊する。Cs134 の生成量は大まかには燃焼度の2乗に比例する。 Cs133はU238と類似して大きな共鳴捕獲が生じるため、 Cs134の生成量はボイド率に依存する。 Cs134がβ崩壊に続いて放出するγ線では、796keVと1.365MeVの測定が容易である。Cs133 の核分裂核種による収率の変化は小さい。
ユーロピウム154(Eu154)は核分裂で生成した安定核種の Eu153が中性子を捕獲して生成し、半減期 8.5年で崩壊する。 Eu153はU238と類似して大きな共鳴捕獲が生じるため、 Eu154の生成量はボイド率に依存する。Eu154 の生成量はおおよそ燃焼度の2乗に比例する。 Eu154がβ崩壊に続いて放出するγ線では、1.274MeVの測定が比較的容易である。Eu153 の核分裂核種による収率の変化は大きいので、燃料組成の影響を受けると共に、Eu153は中性子捕獲断面積が大きく、燃焼度が高い場合にEu154の生成量は若干飽和現象が現われる点にも留意が必要となる。
Ce144 は核分裂した後、中性子捕獲反応を起こすこと無く直接生成し、半減期 284日でβ崩壊してPr144 になり、半減期17.3分でβ崩壊した直後に放出する 2.186MeV のγ線を測定するのが特に容易である。Ce144 の生成量は燃焼度が低い場合には燃焼度に比例するが、半減期が短いため飽和現象が早々に現れる。核分裂核種による収率の変化は概して小さい。したがって生成時のボイド率依存はほとんど無い。半減期が 284日と短いため、本発明で使用する場合には照射期間が短く、すなわち燃料の燃焼度が余り進んでいない状態での利用が好適である。
なお、燃料の燃焼度が低く冷却時間が短い場合には半減期が短いZr95やNb95の利用も考慮に値する。
上述の諸特性がボイド率測定に活用できる根拠を次に詳細に説明する。
図10は、濃縮度4% のBWR燃料集合体をボイド率40%(ボイド割合0.40)で燃焼計算コードORIGENを用いて計算した結果を両対数のグラフでプロットしたものである。照射終了以降の期間、すなわち冷却期間(時間)は1年の場合である。縦軸は燃料金属1トン当たりの毎秒中性子放出率であり、自然対数表示では 10 5 を対数値でゼロとなるように規格化している。燃焼度は燃料金属1トン当たりの燃焼度である。それぞれの曲線はそれぞれの燃焼度範囲において直線になっている部分を見出すことができる。本発明では、直線部では中性子放出率は燃焼度の指数乗に比例することに着目する。勾配は相互に異なっており、この勾配の違いは上記「指数乗」の指数の値が異なっていることを示している。
図11は、図10における各曲線を燃焼度10GWd/tにおいて規格化したものであるが、ガンマ線放出核種の相対変化も付記した。この図に記載していないEu154 は Cs134とほぼ重なる曲線である。 Cs137は半減期が30年と長い(この図の60GWd/t にあたる照射期間は5年程度)ため、ほとんど直線となっている。Cs137 濃度はほとんど燃焼度に比例する(指数= 1.0)ため、他の曲線は燃焼度の指数が1より大きい特性を持っていることがわかる。
勾配が異なる2種の曲線(両対数で直線)の組み合わせにより、ボイド率の評価が可能か否かという点に着目し、前記計算値を用いて検討した例を図12に示す。
この図は冷却期間1年の照射燃料から放出されている Cm242を除く核種からの中性子放出率(S4R)の自然対数値を燃焼度の自然対数値に対してプロットしたもので、パラメータはボイド率(0、40、70 %)である。ボイド率依存計算はボイド率依存の断面積を用いてORIGENコードにより行なった。
燃焼度約18GWd/t から45GWd/t において良好な直線となっており、「 4R =a・E α 、(E:燃焼度、a:比例定数)」の形で近似できることがわかる。ボイド率が高いほど中性子放出率は高くなっている。0% ボイドにおいて燃焼度が20GWd/t から40GWd/t の間のα値は 4.9である。
ところで Cs137(濃度、あるいはγ線強度)は燃焼度に比例するため、「 Cs137=b・E、(b:比例定数)」と書くことができ、この式の両辺をα乗して比を取ると、燃焼度の項が消去できて、
a=S 4R ・(b/Cs137) α
と書くことができる。この式を基準のボイド率(ここでは便宜上0% とする)との比で表現すると、bの値が燃料集合体の軸方向で一定(b/b0 =1)の場合には、
a/a 0 =[S 4R /(S 4R ) 0 ]・[(Cs137) 0 /(Cs137)] α
の形で表現でき、右辺の2個の比の測定値から相対値(a/a0 )を求めることができることがわかる。一方この相対値とボイド率との関係は図12に示すような曲線における基準のボイド率に対する比として計算で求めることができるので、ボイド率の評価が可能となる。
この比は燃焼度に依存しないことが望ましいが、実際には緩やかに燃焼度に依存するので、燃焼度の概算値が要求される。放射能群の組み合わせによってはなるべく正確な燃焼度が必要な場合もあることが後述の例では出てくる。
図13は、前述の図12(請求項)に対応するもので、S4RとCs137 との比であるが、図11に示す2組の曲線をボイド率40% の場合について示している。この図から燃焼度20GWd/tないし50GWd/t の間で本発明を適用できることが、上述の説明から理解できる。
図14は、S4 、すなわち Cm244からの中性子放出率と Cs137との組み合わせを用いる場合の例である(請求項対応)。 Cm244の半減期は18.1年と長いので便利であるが、燃焼度が高いとS4Rの場合と同じであり、燃焼度が低い場合(例えば15GWd/t 以下)には、SR の影にすっかり隠れてしまうため好適ではないと予想できる。10〜20GWd/t付近では計算でSR の寄与分を差し引いて使用することになる。
図15は、S2 、すなわち Cm242からの中性子放出率と Cs137との組み合わせを用いる場合の例である(請求項対応)。 Cm242の半減期が短いので、燃焼度が高い場合(例えば20GWd/t 以上)で使用するのは好適ではないと予想できる。
図16は、SR 、すなわち Cm242と Cm244を除く中性子放出率と Cs137との組み合わせを用いる場合の例である(請求項対応)。測定によってS4 とSR を分離することはできないので、燃焼度が7GWd/t 程度以下でないと実用できないことが、図11から理解できる。
図17は、 Cs134と Cs137との組み合わせ、および Eu154と Cs137との組み合わせの例を示している(請求項対応)。Cs134の半減期が2年程度であるため、2年程度に対応する20GWd/t付近以下での適用が好適である。Eu154は半減期は長いが、捕獲断面積が大きいことから燃焼が進むと Eu154の一部には Eu155になるものがあるため、飽和現象が現れるので、やはり高い燃焼度への適用には補正が必要となる。
図18は、図15において Cs137の代わりに Ce144を用いる例である(請求項対応)。低い燃焼度では冷却時間は余り長くない方が好適であり、その場合 Cs137の測定精度は低下する恐れがあり、半減期は長くないが測定しやすい Ce144を用いる方が得策の場合が考えられる。
図19は、図16において Cs137の代わりに Ce144を用いる例である(請求項対応)。低い燃焼度では冷却時間は余り長くない方が好適であり、その場合 Cs137の測定精度は低下する恐れがあり、半減期は長くないが測定しやすい Ce144を用いる方が得策の場合が考えられるのは図17の場合と同様である。
図20は、図17において Cs137の代わりに Ce144を用いる例である(請求項対応)。Eu154 は燃焼度が低く、冷却時間が短い場合には測定精度が大幅に低下する恐れがあるため、この例では除外してある。
図21〜25は、対零ボイド比の(a/a0 )に相当する値の燃焼度依存性を示しており、値が大きいほどボイドに対する感度が高いことを示している。また燃焼度依存量として見た場合、平坦であれば本発明では燃焼度の情報は必要無いことであり、緩やかに変化している場合には燃焼度の概算値が必要であり、急激な変化をしている場合には燃焼度の正確な情報が必要であることを示している。
図21はS4Rを用いる場合で、燃焼度が約18GWd/t 以上で適当なボイドに対する感度があり、18〜30GWd/t 付近では燃焼度の情報はごく大まかな値であれば良いことが読み取れる。
図22はS4 を用いる場合である。低い燃焼度の場合にはSR の影に隠れて利用できないことを考慮すると、S4Rの場合と大差ない。
図23はS2 を用いる場合である。燃焼度が例えば20GWd/t 以下で用いることになるので、ある程度正確な燃焼度の値が要求されることがわかる。
図24はSR を用いる場合である。燃焼度か低いとき、冷却時間を変えた2回ないし2回以上の中性子測定からS2 とSR とを分離することができ、かつS4 の値は小さく計算で補正できるので便利であるが、ボイドに対する感度が低い点に難点がある。燃焼度の正確な情報は必要無い。
図25は Cs134を用いる場合である。感度の燃焼度依存が概して小さく、燃焼度情報はごく大まかで良い。感度が低い点が難点である。
なお、請求項記載の方法は図10および11におけるS2 とS4Rの勾配の違いを利用する方法であり、燃焼度が10GWd/t ないし7GWd/t以下という非常に低い燃焼度の場合に適用できる方法である。実施手順は第6の実施例において説明する。
以下、具体的な実施例を実施手順として説明する。手順に含まれている詳細な背景は上述したので手順中心に説明する。
[実施例1]
図1に示す手順は請求項1に対応する軸方向ボイド率分布測定方法である。第1群の放射能はその強度(A=S:中性子放出率)が燃焼度(E)の指数(α)乗に比例するものとして、S=a・E α (ただしaは比例定数)で表すことができる中性子束または中性子計数率(φ)と、第2群の放射能はその強度(B=γ:ガンマ線強度)を燃焼度(E)に比例するものとして、γ=g・E(ただしgは比例定数)で表すことができる核分裂生成物のCs137 或いは Ce144のガンマ線強度(γ)とを測定する。そして、ボイド率既知または評価が容易な位置、特にボイド率0% の位置などの基準点の値に対する相対値[(φ/φ0 )および(γ/γ0 )]を求め、(式1)を具体化した(式1a)、すなわち
(式1a): a/a 0 =(S/S 0 )・(γ 0 /γ) α ・(g/g 0 α
により(a/a0 )の値を求め、計算で求めた(a/a0 )の値とボイド率との相関を示す校正曲線を用いて軸方向ボイド率分布を求める。
ここで、比例定数のgは燃料集合体の軸方向でほとんど変化しない(g/g0 =1)が、燃料の設計によっては微妙に変化する場合がある。その場合には必要に応じて燃料棒に起因するガンマ線遮蔽効果を計算により補正することになる。(S/S0 )の値は(φ/φ0 )の値から求められる。厳密には(φ/φ0 )の値に中性子増倍率の軸方向変化の寄与分を補正することになるが、BWR燃料集合体ではそれを考慮しなくても(φ/φ0 −1)に比べてかなり小さく、増倍率変化効果はその10〜20% 程度であり、近似的に「S/S0=φ/φ0 」としても差支えない場合が多い。
ところで、BWR燃料集合体の有効部を慣例にならって24等分し、それぞれの区間をノードと呼び、下端(冷却水上流)から1、2、…、24ノードと呼ぶことにすると、ボイド率はノード2〜3ではほぼ0% 、上端付近の23〜24ノードでは70〜75% となることがわかっている。原子炉内の炉心高さ方向のボイド率分布形は原子炉の運転方法や燃料の燃焼度の進行に伴って変化する。本実施例では、S4Rと Cs137を用いて使用済燃料集合体について測定するときには、両者の半減期がそれぞれ長いことから、全燃焼期間の平均的なボイド率分布を測定することになる。後出の実施例では、測定対象の放射能の半減期が長い場合には燃焼期間の平均的な、測定対象の放射能の半減期が短い場合には炉停止直前におけるボイド率分布を測定することになる。
[実施例2]
図2に示す手順は請求項に対応する軸方向ボイド率分布測定方法であって、実施例1と異なる点として、(φ/φ0 )の値に中性子増倍率の軸方向変化の寄与分を設計計算で得られる標準的な中性子増倍率の分布を補正計算に用いる方法である。前述のように、BWR燃料集合体では中性子増倍率変化の影響は(φ/φ0 −1)に比べてかなり小さく、その10〜20% 程度であり、設計計算値を用いてもほとんど問題ない。中性子束または計数率(φ)は中性子増倍率(k)および中性子放出率(S)との間に、比例定数(c)を介して「φ=cS/(1−k)」の関係があることは良く知られており、この関係式を用いる。
比例定数(c)は、例えば「固定源計算法」(中性子源強度を与えて中性子輸送拡散現象を解き、目的の位置の中性子束を求める方法)を用い、組成と中性子源強度分布も共通で、仮想的に中性子増倍効果が無い場合(k=0)の中性子束(φNM:非増倍)と増倍効果が正規に存在する場合の中性子束(φM :増倍)との比から(k)の値を求め、その値を「φ=cS/(1−k)」の関係式に代入して求めることができる。Sの値は計算の入力値であるため既知量である。
[実施例3]
図3に示す手順は請求項10に対応する軸方向ボイド率分布測定方法であって、実施例2と異なる点は、(φ/φ0 )の値に中性子増倍率の軸方向変化の寄与分を、電力申告値として最も入手容易と考えられる「燃料集合体平均燃焼度」を用いて補正する点にある。しかしながら、燃料集合体の軸方向ボイド率分布を与えられる可能性はかなり低いと考えられるので、仮想したボイド率分布を測定値を用いて繰り返し計算により修正しながらボイド率分布を決定し、そのボイド率に対応する燃焼度依存の無限増倍率(k )を求め、その値から中性子増倍率(k)を求めるという新規な方法を考案するに至った。以下詳細にその手順を説明する。
軸方向の中性子束(または中性子計数率)分布とガンマ線強度分布は実施例1および2と同様に測定され、(γ/γ 0 α も同様に求めることができるが、(φ/φ0 )から(S/S0 )を単純に求めることはできない。そこで中性子増倍率変化を補正する因子が(φ/φ0 −1)に比べてかなり小さいため、中性子増倍率変化補正因子を繰り返し計算によって求めることができること(繰り返し計算で収束させることが可能であること)が見出された。
軸方向のガンマ線強度分布(特に Cs137分布)の平均値を電力申告値の燃料集合体平均燃焼度に規格化することにより、軸方向燃焼度分布(E)が求められる。中性子増倍率(k)は無限増倍率(k )に計算で求められる一定の値(Fk )を掛けて求めることができるが、(k )の値は燃焼中の冷却水ボイド率の影響を受ける。
図3の手順に関して述べる前に(k )と(Fk )について説明する。本実施例を実施するのは室温でボイドが無く、燃焼時は高温高圧でボイドがあるという点で混乱を招く恐れがあるからである。
高温高圧でボイド率を与え、燃料集合体の燃焼計算(正確な設計計算では中性子スペクトル計算を行ないながら燃焼計算も行なう。簡易計算では燃焼途中のスペクトル計算を行なわない場合もある。)を行ない、燃焼に伴う燃料組成の変化を求め、次に室温(冷温、冷態ともいう)において燃焼により変化した燃料組成を用いて中性子スペクトル計算を行ない、無限増倍率を求めると共に、中性子輸送拡散計算を行なうための群定数を求める。
前記固定源計算法を用いて、測定する体系に対して中性子輸送拡散計算を行なって、測定位置において観測される中性子増倍率(k)を求める。この値は通常用いられている実効増倍率とは水中に置かれたBWR燃料集合体の場合、定量的には僅かに異なっているが、本実施例では無視できる程度の差異しかない。測定位置が燃料集合体から2〜3cm以上離れると、水中では中性子増倍率と無限増倍率との比(Fk )は燃料組成の影響を受けないことが本発明者らの研究によってすでに明らかになっている。「測定体系の無限増倍率」が「ボイド率の影響を受ける」ということは、燃焼時の冷却水のボイド率によって組成が変わるため、無限増倍率も変化するということである。
以上の認識の下、図3の手順を次に説明する。無限増倍率(k )は燃焼度(E)が進むにつれてほぼ直線的に減少するが、燃焼度が高くなるとプルトニウムの生成による中性子増倍が現れてくるので、図3では燃焼度の2次式で近似している。すなわち、
=J 0 −J 1 ・E+J 2 ・E 2
ただし、1次式近似で問題ないことも多い。2次式で近似する場合の定数(J0 )および係数(J1 、J2 )はボイド率に依存するものであり、依存特性はボイド率(割合)を変えた計算によって求めることができる。燃料集合体の上下で初期濃縮度を若干変えた設計の場合には、これらの定数および係数に濃縮度依存性を含めることになる。
こうして計算により、無限増倍率と燃焼度に関する相関式がボイド依存量として表現される。したがって、燃焼度と仮定した(繰り返し計算の開始点であるため、i=1とする)ボイド率によって無限増倍率(k∞i)の軸方向分布が得られ、(Fk )を介して、軸方向中性子増倍率(ki )分布が得られる。この値と測定値(φ/φ0 )とから(S/S0 i の値が求められ、(a/a0i の軸方向分布が求められ、計算で求めたこの比とボイド率との関係から、仮定ボイド率(i=1)から修正された軸方向ボイド率(vi 、i=2)の分布が得られる。この修正された分布を修正される前の分布と比較し、許容範囲で一致していれば軸方向分布は求められたとし、一致が不十分な場合には修正された分布をk と燃焼度との相関式に代入する繰り返し計算を行なうことにより、収束した値を軸方向ボイド率分布とするものである。
[実施例4]
図4に示す手順は請求項11に対応する軸方向ボイド率分布測定方法であって、実施例3と異なる点は、(φ/φ0 )の値に中性子増倍率の軸方向変化の寄与分を、燃料から放出されている中性子放出率から、本発明者らが開発した「中性子放出率法」(特許文献3)と本質的に同じ方法を用いて基準点における燃焼度(E0 )を求めて利用する点である。他の部分は実施例3の場合と同じであるため、異なる部分を、図5(基準点に対する中性子放出率法測定手順)を用いて説明する。
中性子放出率法では、燃料集合体側面の中性子束(φ)を測定する。本測定では中性子束の絶対値が必要であり、例えば金箔放射化法、あるいはその方法などにより校正した中性子検出器を使用する方法がある。
中性子束または中性子計数率(φ)は中性子増倍率(k)および中性子放出率(S)との間に、比例定数(c)を介して「φ=c・S/(1−k)」の関係がある。比例定数(c)は前述のように計算で求めることができる。繰り返し計算の最初(i=1)に中性子増倍率の初期値(k1 )を与え、「Si =φi ・(1−ki )/c」(全中性子放出率)を求める。Cm242 からの中性子放出率の補正が必要な場合には、計算により、または半減期の違いを利用して冷却時間を異にする2回ないし2回以上測定することによりS2 成分を除去する。すなわち、「(S4Ri =Si −S2 」の演算を行なう。この値から計算で求められるボイド率0% に対するS4Rと燃焼度との相関を示す校正曲線を用いて(Ei )が求められる。
さらに、計算で求められたk とE、S4RとEとの相関を示す校正曲線を用いてk i が求められる。この値から前述の測定体系に対する換算因子(Fk )を用いて修正された中性子増倍率ki が求められる。修正された中性子増倍率と修正される前、または直前に修正されていた中性子増倍率との差異が許容範囲であれば最終値とし、一致が不十分な場合には初期値または前回の修正値を入れ替えて繰り返し計算を行ない、収束値を求める。
繰り返し計算が必要な理由は、超越関数が本手法の中に実質的に含まれているためであることがわかっており、実用的には3〜5回程度で収束する。中性子放出率法では、燃焼度の他に、中性子放出率、プルトニウム濃度、無限増倍率、中性子増倍率なども同時的に求めることができるが、図4の手順では燃焼度のみを利用している。
[実施例5]
図6に示す手順は請求項、12に対応する軸方向ボイド率分布測定方法であって、上述の実施例1〜4と異なり、ガンマ線スペクトル測定のみによって目的を達成するものであり、従来の十分確立された計測技術では通常Geなどの半導体検出器が用いられる。照射燃料の側面にガンマ線コリメータを設置してガンマ線を導き、半導体検出器などによりガンマ線スペクトルを測定する。これによって目的のCs137、Cs134 、Ce144 、Eu154 などから放出されているガンマ線の強度を求めることができる(ガンマ線スペクトル分析法)。
これらはどの時点で測定しても全てを精度良く測定できるものではないが、どのような時点や条件で測定できるかは専門家の間では周知である。図6ではCs134 とCs137 、Eu154 とCs137 、およびCs134 とCe144 との組み合わせを「比」の形とし、その軸方向分布(G/G0 )を求め、計算で求めた「(G/G0 )とボイド率との相関」を示す校正曲線を使用して軸方向ボイド率分布が求められる。
[実施例6]
図7に示す手順は請求項に対応する軸方向ボイド率分布測定方法であって、上述の実施例1〜5と異なり、中性子放出率測定のみによって、燃焼度が低い場合に好適に実施することができるものであって、Cm242 からの中性子放出率(S2 )とそれを除く核種からの中性子放出率(S4R)との比を求め、この比の軸方向分布を計算で求めた校正曲線を使用して軸方向ボイド率分布が求められる。実施の手順は実施例5の場合と酷似している。
具体的には、原子炉内で中性子に照射された照射燃料集合体から放出されている中性子捕獲に起因する中性子放出率を、所定の冷却時間を挟んで少なくとも2回測定する。そして、半減期の違いを利用してCm242 から放出されている第1の中性子放出率(S2 )と、その他核種から放出される第2の中性子放出率(S4R)とに分類する。測定された第1の中性子放出率(A=S2 )を燃焼度(E)の指数(α)乗に比例するものとして、A=a・E α (ただしaは比例定数)で表す。
さらに、測定される第2の中性子放出率(B=S4R)を燃焼度(E)の指数(β)乗に比例するものとして、B=b・E β (ただしbは比例定数)で表し、第1の中性子放出率の式の両辺をβ乗し、第2の中性子放出率の式の両辺をα乗して両式の比を取って燃焼度の因子を消去して得られる式を基準場所(「0」を付す)における値との比として求めた(式2)、すなわち
(式2): (a/a 0 β ・(b 0 /b) α =(A/A 0 β ・(B 0 /B) α
を得る。ここで、(a/a 0 β ・(b 0 /b) α が原子炉内部の軸方向ボイド率分布に依存する特性を利用して、(A/A0 )および(B0 /B)の測定から、別途求められた校正曲線を介して軸方向ボイド率分布を求める。
この方法は中性子測定のみで実施できる点で大きな特徴があるが、中性子放出率(S4R)の中でのCm244 からの中性子放出率が小さいことという条件が含まれている。すなわち、(a/a0 )と(b0/b)のボイド率分布に対する依存性が逆であるため、結果として(a/a 0 β ・(b 0 /b) α のボイド率分布依存性は相殺され、低い燃焼度の場合にしか適用できないという難点がある。また、Cm242を利用するため照射終了後の冷却時間はせいぜい2年以内でないと実施できない、ほかに、感度が燃焼度で大きく変化するため、燃焼度のデータが要求される、などという難点がある。
[実施例7]
図8は燃料集合体収納装置への収納前の燃料集合体の中性子増倍率を評価する方法に関するもので、本発明の第2の目的を達成する請求項14に対応するものである。照射済(使用済)BWR燃料集合体では、前述のように、軸方向燃料有効部のうち、下端から2/3ないし3/4付近で、無限増倍率が高くなる部位が存在する。この部位が収納装置に燃料集合体を収納した場合、未臨界度を最も厳しくする(臨界未満で臨界に最も近付く)ので、従来は過度の設計裕度を取って臨界安全性を確保する必要があったが、本実施例では燃料集合体ごとに前記部位を実測することができるので設計裕度を大幅に低減できる。
本実施例の内容は以下に述べるグロスガンマ線測定ではなく、通常装置が大きくデータ処理も複雑になるガンマ線スペクトル分析法でも当然実施できるが、実用性を考慮し、測定(評価)精度は若干低下する可能性はあるものの、装置がより小型で単純なグロスガンマ線測定法を用いる方法を提供するものである。
すなわち、本実施例では、照射燃料集合体の基準点「0」と未臨界度モニタ点に対して中性子測定を行なって(φ0 )および(φ)を求めると共に、グロスガンマ線強度比[(γg )/(γg0 ]を測定し、中性子測定により前記「中性子放出率法」の手法により、基準点の中性子増倍率(k0 )および中性子束または中性子計数率の比(φ0 /φ)を求める。そして、別途求めたグロスガンマ線強度相対分布と燃焼度相対分布との相関関係(計算値、分布形の差異は通常10% 程度以下)を用いて前記グロスガンマ線強度比[(γg )/(γg 0 から燃焼度相対比へ換算し、計算で求めた燃焼度(E)の指数(α)乗値を用いて(E/E 0 α を求める。ここで、前記諸請求項等を用いて同種の燃料集合体について複数体の測定を行うことなどにより、代表的なボイド率分布が求められたとき、ここの未臨界度モニタ点のボイド率を与えることができる。このボイド率から、「a/a0 とボイド率との相関関係」を逆に用いて(a/a0 )の値を求め、
(式5): k=1−(1−k 0 )・(φ 0 /φ)・(E/E 0 α ・(a/a 0
により、燃料集合体収納装置へ収納する前に未臨界度モニタ点の中性子増倍率を求めるものである。
なお、グロスガンマ線測定法は本発明者らによってすでに確立され、日本の商業用再処理工場で採用されている。また、照射済(使用済)加圧水型原子炉(PWR)用燃料集合体でも、ボイドに起因するものではなく水温、水密度、制御棒効果、等が原因で類似の現象が発生することがわかっているので、因子(a/a0 )にそれらの影響を計算により取り込むことによって、本実施例を実施することができる。
第1実施例の測定手順(請求項1対応)を示すフロー図。 第2実施例の測定手順(請求項10対応)を示すフロー図。 第3実施例の測定手順(請求項11対応)を示すフロー図。 第4実施例の測定手順(請求項12対応)を示すフロー図。 基準点に対する中性子放出率法測定手順を示すフロー図。 第5実施例の測定手順(請求項13対応)を示すフロー図。 第6実施例の測定手順(請求項9対応)を示すフロー図。 未臨界度モニタ点(軸方向高さ)の中性子増倍率簡易評価法(請求項14対応)を示すフロー図。 超ウラン核種の生成・消滅過程を説明するフロー図。 中性子放出率と燃焼度との関係を示すグラフ。 中性子放出率(相対値)および Cs134、 Cs137濃度(相対値)と燃焼度との関係を示すグラフ。 燃焼度依存中性子放出率(S4R)のボイド率依存性を示すグラフ。 4Rと Cs137との組み合わせの燃焼度依存性(請求項3対応)を示すグラフ。 4 と Cs137との組み合わせの燃焼度依存性(請求項4対応)を示すグラフ。 2 と Cs137との組み合わせの燃焼度依存性(請求項5対応)を示すグラフ。 はSR と Cs137との組み合わせの燃焼度依存性(請求項6対応)を示すグラフ。 Cs134と Cs137との組み合わせ、および Eu154と Cs137との組み合わせの燃焼度依存性(請求項7、8対応)を示すグラフ。 図15において Cs137の代りに Ce144を用いる例(請求項5対応)を示すグラフ。 図16において Cs137の代りに Ce144を用いる例(請求項6対応)を示すグラフ。 図17において Cs137の代りに Ce144を用いる例(請求項7、8対応)を示すグラフ。 4Rを用いる場合の対零ボイド比の(a/a0 )に相当する値の燃焼度依存性を示すグラフ。 4 を用いる場合の対零ボイド比の(a/a0 )に相当する値の燃焼度依存性を示すグラフ。 2 を用いる場合の対零ボイド比の(a/a0 )に相当する値の燃焼度依存性を示すグラフ。 R を用いる場合の対零ボイド比の(a/a0 )に相当する値の燃焼度依存性を示すグラフ。 Cs134を用いる場合の対零ボイド比の(a/a0 )に相当する値の燃焼度依存を示すグラフ。

Claims (13)

  1. 原子炉内で中性子に照射された照射燃料集合体において、重核種が中性子捕獲を行って生成しかつ中性子を放出する放射性核種に起因する第1群の放射能と、核分裂で生成した後、中性子捕獲反応を受けておらずかつ所定の範囲で生成量が燃焼度に比例しかつガンマ線を放出する核分裂生成物に起因する第2群の放射能とを測定し、
    測定された第1群の放射能の強度である中性子放出率Aを燃焼度Eの指数(α)乗に比例するものとして、比例定数をaとして、A=a・Eαで表すとともに、
    測定される第2群の放射能の強度であるガンマ線強度Bを燃焼度Eに比例するものとして、比例定数をbとして、B=b・Eで表し、
    第2群の放射能強度の式の両辺をα乗して第1群の式との比をとって燃焼度の因子を消去して得られる式を、「0」を付して示す基準場所における値との比として求めた(式1)、すなわち
    (式1): a/a0 =(A/A0)・(B0/B)α・(b/b0α
    において、
    (b/b0)は1とみなし、(A/A0 )および(B0 /B)の測定値、および計算によって求めた指数αから(a/a0 )を求め、(a/a0 )が原子炉内部の軸方向ボイド率分布に依存する特性を利用して、別途求められた「a/a0とボイド率との相関関係」を用いて軸方向ボイド率分布を求めること、
    を特徴とする軸方向ボイド率分布測定方法。
  2. 中性子放出率はキュリウム242を除く核種から放出される中性子放出率であり、核分裂生成物はセシウム137であることを特徴とする請求項1記載の軸方向ボイド率分布測定方法。
  3. 中性子放出率はキュリウム244から放出される中性子放出率であり、核分裂生成物はセシウム137またはセリウム144であることを特徴とする請求項1記載の軸方向ボイド率分布測定方法。
  4. 中性子放出率はキュリウム242から放出される中性子放出率であり、核分裂生成物はセシウム137またはセリウム144であることを特徴とする請求項1記載の軸方向ボイド率分布測定方法。
  5. 中性子放出率は、キュリウム242およびキュリウム244から放出される中性子放出率を除く中性子放出率であり、核分裂生成物はセシウム137またはセリウム144であることを特徴とする請求項1記載の軸方向ボイド率分布測定方法。
  6. 原子炉内で中性子に照射された照射燃料集合体において、核分裂生成物がさらに中性子捕獲して生成しかつガンマ線を放出する放射性核種に起因する第1群の放射能と、核分裂で生成した後、中性子捕獲反応を受けておらずかつ所定の範囲で生成量が燃焼度に比例しかつガンマ線を放出する核分裂生成物に起因する第2群の放射能とを測定し、
    測定された第1群の放射能の強度であるガンマ線強度Aを燃焼度Eの指数(α)乗に比例するものとして、比例定数をaとして、A=a・Eαで表すとともに、
    測定される第2群の放射能の強度であるガンマ線強度Bを燃焼度Eに比例するものとして、比例定数をbとして、B=b・Eで表し、
    第2群の放射能強度の式の両辺をα乗して第1群の式との比をとって燃焼度の因子を消去して得られる式を、「0」を付して示す基準場所における値との比として求めた(式1)、すなわち
    (式1): a/a0 =(A/A0)・(B0/B)α・(b/b0α
    において、
    (b/b0)は1とみなし、(A/A0 )および(B0 /B)の測定値、および計算によって求めた指数αから(a/a0 )を求め、(a/a0 )が原子炉内部の軸方向ボイド率分布に依存する特性を利用して、別途求められた「a/a0とボイド率との相関関係」を用いて軸方向ボイド率分布を求めること、
    を特徴とする軸方向ボイド率分布測定方法。
  7. 第1群の放射能を有する核種がセシウム134またはユーロピウム154であり、第2群の放射能を有する核種がセシウム137またはセリウム144であることを特徴とする請求項6記載の軸方向ボイド率分布測定方法。
  8. 原子炉内で中性子に照射された照射燃料集合体から放出されている、重核種が中性子捕獲を行って生成する放射性核種に起因する中性子放出率を、所定の冷却時間を挟んで少なくとも2回測定し、
    半減期の違いを利用してキュリウム242から放出される第1の中性子放出率とその他核種から放出される第2の中性子放出率とに分類し、
    測定された第1の中性子放出率Aを燃焼度Eの指数(α)乗に比例するものとして、A=a・Eα(ただしaは比例定数)で表すとともに、
    測定される第2の中性子放出率Bを燃焼度Eの指数(β)乗に比例するものとして、比例定数をbとして、B=b・Eβで表し、
    0」を付して示す基準場所における値との比として求めた(式2)、すなわち
    (式2): (a/a0β・(b0/b)α=(A/A0β・(B0/B)α
    において、
    A/A0 およびB0 /Bの測定値および計算によって求めた指数αとβから、(a/a0β・(b0/b)αを求め、別途求められた「(a/a0β・(b0/b)αとボイド率との相関関係」を用いて軸方向ボイド率分布を求めること、を特徴とする軸方向ボイド率分布測定方法。
  9. ガンマ線強度γの軸方向分布と中性子束または中性子計数率φの軸方向分布を測定してそれぞれ基準場所における値γ0 およびφ0 に対する相対値γ/γ0 およびφ/φ0 を求め、
    一方では計算によって求めた指数αを用いて(γ0/γ )αを、他方では(式3)、すなわち
    (式3): S/S0=(φ/φ0)・(1−k)/(1−k0
    ただし、(S/S0 )は基準場所における値に対して規格化した中性子放出率の軸方向分布、k0基準場所における中性子増倍率、kはφ測定点における中性子増倍率
    で表し、
    中性子増倍率分布(kおよびk0 )に計算値を用いて(式3)から(S/S0 )を求め、(式1)に対応する(式4)すなわち
    (式4): a/a0=(S/S0)・(γ0/γ)α
    によりa/a0 を算出し、
    別途計算で求めた「a/a0 とボイド率との相関」を示す校正曲線を用いて軸方向ボイド率分布を求めること、を特徴とする請求項1乃至5のいずれか記載の軸方向ボイド率分布測定方法。
  10. 燃料集合体平均の燃焼度が与えられた場合に、
    軸方向ガンマ線分布測定値の平均値をこの申告値に規格化することによって軸方向燃焼度分布(E、E0 )を求め、
    一方、無限増倍率k を、少なくともボイド率に依存する係数を用いた燃焼度Eの1次式または2次式で表し、前記軸方向燃焼度分布と仮定した軸方向ボイド率分布とを用いて軸方向無限増倍率k 分布を求め、
    測定条件に対応する中性子輸送拡散計算により換算因子Fk を求めて、k 分布にFk をかけることにより軸方向中性子増倍率k分布へ換算し、一方、軸方向中性子束分布から(φ/φ0)を求め、これらと(式3)すなわち
    (式3): S/S0=(φ/φ0)・(1−k)/(1−k0
    ただし、(S/S0 )は基準場所における値に対して規格化した中性子放出率の軸方向分布、k0基準場所における中性子増倍率、kはφ測定点における中性子増倍率
    とからS/S0 を求め、
    一方では、計算によって求めた指数αと軸方向ガンマ線分布から(γ0 /γ)αを求めて(式4)すなわち
    (式4): a/a0=(S/S0)・(γ0/γ)α
    により(a/a0 )を算出し、
    別途計算で求めた「a/a0 とボイド率との相関関係」を用いて軸方向ボイド率分布を求め、
    この分布が前記軸方向ボイド率分布仮定値または繰り返し計算における直前に求められた軸方向ボイド率分布と指定された範囲で一致していれば最終結果とし、一致しない場合には初期値として与えた軸方向ボイド率分布仮定値あるいは直前の繰り返し前のボイド率分布と置換することによって繰り返し計算を行ない、最終結果が得られて収束した値を軸方向ボイド率分布とすること、を特徴とする請求項1乃至5のいずれか記載の軸方向ボイド率分布測定方法。
  11. 基準場所における燃焼度E0 を中性子放出率法によって求め、軸方向ガンマ線分布測定値のうち基準場所における値γ0 を前記燃焼度E0 に規格化することによって軸方向燃焼度分布(E、E0 )を求め、
    一方、無限増倍率k を少なくともボイド率に依存する係数を用いた燃焼度Eの1次式または2次式で表し、
    前記燃焼度分布と仮定した軸方向ボイド率分布とを用いて軸方向無限増倍率k 分布を求め、
    測定条件に対応する中性子輸送拡散計算により換算因子Fk を求めてk 分布にFk をかけることにより軸方向中性子増倍率(k)分布へ換算し、一方、軸方向中性子束分布から(φ/φ0)を求め、
    これらと(式3)すなわち
    (式3): S/S0=(φ/φ0)・(1−k)/(1−k0
    ただし、S/S0 基準場所における値に対して規格化した中性子放出率の軸方向分布、k0基準場所における中性子増倍率、kはφ測定点における中性子増倍率
    とからS/S0 を求め、
    一方では、計算によって求めた係数αと軸方向ガンマ線分布から(γ0 /γ)αを求めて、(式4)すなわち
    (式4): a/a0=(S/S0)・(γ0/γ)α
    によりa/a0 を算出し、
    別途計算で求めた「a/a0 とボイド率との相関関係」を用いて軸方向ボイド率分布を求め、
    この分布が前記軸方向ボイド率分布仮定値または繰り返し計算における直前に求められた軸方向ボイド率分布と指定された範囲で一致していれば最終結果とし、一致しない場合には初期値として与えた軸方向ボイド率分布仮定値あるいは直前の繰り返し前のボイド率分布と置換することによって繰り返し計算を行ない、最終結果が得られて収束した値を軸方向ボイド率分布とすること、を特徴とする請求項1乃至5のいずれか記載の軸方向ボイド率分布測定方法。
  12. 燃料集合体の軸方向位置を変えながらガンマ線スペクトルを測定し、
    それぞれの高さ位置において セシウム134または ユーロピウム154と セシウム137または セリウム144のフォトピークを測定してガンマ線強度比G、すなわち Cs134/Cs137 比、 Eu154/Cs137 比、または Cs134/Ce144 比の 0 」を付して示す基準場所における値で規格化した軸方向分布G/G0 を求め、
    G/G0 値とボイド率との相関を示す校正曲線を用いて軸方向ボイド率分布を求めること、を特徴とする請求項6または7記載の軸方向ボイド率分布測定方法。
  13. 請求項1ないし12のいずれか一項に記載の軸方向ボイド率分布測定方法によって軸方向ボイド率分布を求め、
    照射燃料集合体の基準場所と未臨界度モニタ点に対して中性子測定を行なって、
    基準場所における値に「 0 」を付して表すものとして、
    基準場所における中性子増倍率k0 および中性子束または中性子計数率の比φ0 /φを求め、
    また、グロスガンマ線強度比[(γg )/(γg0 ]を測定し、
    別途求めたグロスガンマ線強度相対分布と燃焼度相対分布との相関関係(計算値)を用いてグロスガンマ線強度比[(γg )/(γg0 ]から燃焼度相対比(E/E0)へ換算し、
    算で求めた指数αを用いて(E/E0αを求め、
    前記軸方向ボイド率分布と、「a/a0 とボイド率との相関関係」を用いてa/a0 の値を求め、
    (式5)、すなわち
    (式5): k=1−(1−k0)・(φ0/φ)・(E/E0α・(a/a0
    により未臨界度モニタ点の中性子増倍率kを求めること、を特徴とする燃料集合体収納装置への収納前の燃料集合体中性子増倍率評価方法。
JP2005143884A 2005-05-17 2005-05-17 軸方向ボイド率分布測定方法および収納装置収納前の燃料集合体中性子増倍率評価方法 Expired - Fee Related JP4761829B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2005143884A JP4761829B2 (ja) 2005-05-17 2005-05-17 軸方向ボイド率分布測定方法および収納装置収納前の燃料集合体中性子増倍率評価方法
US11/434,195 US20070076839A1 (en) 2005-05-17 2006-05-16 Axial void fraction distribution measurement method and neutron multiplication factor evaluating method
TW095117333A TWI320189B (en) 2005-05-17 2006-05-16 Axial void fraction distribution measurement method and neutron multiplication factor evaluating method
US12/345,438 US8401141B2 (en) 2005-05-17 2008-12-29 Axial void fraction distribution measurement method and neutron multiplication factor evaluating method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005143884A JP4761829B2 (ja) 2005-05-17 2005-05-17 軸方向ボイド率分布測定方法および収納装置収納前の燃料集合体中性子増倍率評価方法

Publications (2)

Publication Number Publication Date
JP2006322727A JP2006322727A (ja) 2006-11-30
JP4761829B2 true JP4761829B2 (ja) 2011-08-31

Family

ID=37542533

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005143884A Expired - Fee Related JP4761829B2 (ja) 2005-05-17 2005-05-17 軸方向ボイド率分布測定方法および収納装置収納前の燃料集合体中性子増倍率評価方法

Country Status (3)

Country Link
US (2) US20070076839A1 (ja)
JP (1) JP4761829B2 (ja)
TW (1) TWI320189B (ja)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012122929A (ja) * 2010-12-10 2012-06-28 Toshiba Corp 原子炉燃料の非破壊燃焼度評価法および非破壊燃焼度評価装置
JP5752467B2 (ja) * 2011-04-07 2015-07-22 株式会社東芝 原子炉燃料非破壊燃焼度評価方法およびその装置
JP6257129B2 (ja) * 2012-07-30 2018-01-10 一般財団法人電力中央研究所 管理対象物の臨界管理方法
JP6038575B2 (ja) * 2012-09-27 2016-12-07 株式会社東芝 核燃料燃焼度評価装置、方法およびプログラム
RU2528401C1 (ru) * 2013-06-20 2014-09-20 Федеральное государственное бюджетное учреждение "Национальный исследовательский центр "Курчатовский институт" Способ измерения нейтронной мощности ядерного реактора в абсолютных единицах
JP6427015B2 (ja) * 2015-01-26 2018-11-21 原子燃料工業株式会社 原子炉長期停止時の炉心特性評価及びそのシステム
CN110739091B (zh) * 2019-10-22 2022-05-10 上海核工程研究设计院有限公司 一种核电厂次临界度直接测量系统及硬件平台和测量方法
TWI773001B (zh) * 2019-12-06 2022-08-01 美商西屋電器公司 使用釩中子探測器的方法和用具
CN113096841A (zh) * 2021-03-29 2021-07-09 中国科学院上海应用物理研究所 一种反应堆系统防核扩散性能评估方法
CN113566970B (zh) * 2021-06-22 2023-06-16 中国辐射防护研究院 一种Pu-238同位素热源的搜寻方法和装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5011037B1 (ja) * 1969-10-07 1975-04-26
JPS55121195A (en) * 1979-03-13 1980-09-18 Nippon Atomic Ind Group Co Boiling water reactor
JP3103361B2 (ja) * 1990-02-15 2000-10-30 株式会社東芝 原子燃料の燃焼度測定方法
US5315117A (en) * 1993-03-04 1994-05-24 Texaco Inc. Volume meter system
JP3651716B2 (ja) * 1996-06-17 2005-05-25 株式会社東芝 原子炉燃料の非破壊燃焼度評価法
JPH1039083A (ja) * 1996-07-18 1998-02-13 Toshiba Corp 炉内情報監視装置
JP3564239B2 (ja) * 1996-10-11 2004-09-08 東京電力株式会社 原子炉炉心監視方法
JPH11101890A (ja) * 1997-09-26 1999-04-13 Toshiba Corp 原子炉内出力監視装置
JP2003043183A (ja) * 2001-07-30 2003-02-13 Toshiba Corp 照射燃料の発熱率測定方法
JP4358026B2 (ja) * 2004-05-07 2009-11-04 株式会社東芝 沸騰水型原子炉のボイド率分布測定方法

Also Published As

Publication number Publication date
US20090252271A1 (en) 2009-10-08
TW200713333A (en) 2007-04-01
JP2006322727A (ja) 2006-11-30
TWI320189B (en) 2010-02-01
US8401141B2 (en) 2013-03-19
US20070076839A1 (en) 2007-04-05

Similar Documents

Publication Publication Date Title
JP4761829B2 (ja) 軸方向ボイド率分布測定方法および収納装置収納前の燃料集合体中性子増倍率評価方法
JP5546174B2 (ja) 放射性廃棄物の放射能濃度評価方法及び評価プログラム、並びに放射能濃度評価装置
JP5752467B2 (ja) 原子炉燃料非破壊燃焼度評価方法およびその装置
JP6038575B2 (ja) 核燃料燃焼度評価装置、方法およびプログラム
Bolind The use of the BIC set in the characterization of used nuclear fuel assemblies by nondestructive assay
Kuramoto et al. Absolute measurement of βeff based on Feynman-α experiments and the two-region model in the IPEN/MB-01 research reactor
JP4327522B2 (ja) 燃料物質内プルトニウム238の非破壊測定方法
JP3103361B2 (ja) 原子燃料の燃焼度測定方法
Blanchet et al. Qualification of a gamma-ray heating calculation scheme for the future Jules Horowitz material testing reactor (RJH)
JP2003043183A (ja) 照射燃料の発熱率測定方法
Mulyana et al. A gamma spectroscopy-based non-destructive approach for pebble bed reactor safeguards
JP3026455B2 (ja) 照射燃料集合体の燃焼度測定方法
JP2006138795A (ja) 中性子検出器感度校正方法および燃料収納体系の未臨界度評価方法
Husnayani et al. Kr-85m activity as burnup measurement indicator in a pebble bed reactor based on ORIGEN2. 1 Computer Simulation
Lundkvist Investigation of possible non-destructive assay (NDA) techniques for at the future Swedish encapsulation facility
Eigenbrodt Spent fuel measurements: passive neutron albedo reactivity (PNAR) and photon signatures
Chen et al. Non-destructive determinaton of burn-up by gamma-scanning: an assessment of 144 Ce/Pr as a fission monitor in CANDU fuels
JPH0426718B2 (ja)
JP5951538B2 (ja) 核燃料燃焼度評価装置、その評価方法およびそのプログラム
Pergreffi et al. Gamma dose rates from a spent UO2 fuel Assembly: Calculations vs measurements
Kaplan-Trahan et al. Spent Fuel Measurements
JP2023037880A (ja) 放射能評価方法及び崩壊熱評価方法
Phillips et al. Nondestructive verification with minimal movement of irradiated light-water-reactor fuel assemblies
JPS6316298A (ja) 使用済核燃料集合体の非破壊測定方法
Terremoto Gamma-ray spectroscopy on irradiated fuel rods

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20071010

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20091130

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100126

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100329

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110208

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110406

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20110420

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: 20110510

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110607

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140617

Year of fee payment: 3

R151 Written notification of patent or utility model registration

Ref document number: 4761829

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: 20140617

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees