(57)【要約】【課題】加速器の運転および運転調整時における運転パターンのパラメータ調整が容易であるとともに、加速器を適切に制御できることを目的とする。【解決手段】加速空洞６２を駆動するための高周波信号のパターンデータを、複数の制御期間毎に、時間に対する一次関数を定義するデータで記憶する記憶部３１と、パターンデータに基づいて、制御データを生成する制御データ生成部３３と、制御データに基づいて、加速空洞６２を逐次制御する制御部（３０、５０、５１）と、を備え、制御データ生成部３３は、複数の制御期間に応じた制御データを、それぞれ時間に対する一次関数上のデータとして生成するとともに、生成した制御データのうち、第１の制御期間の制御データを時間軸上でつないだ第１の直線と、第２の制御期間の制御データを時間軸上でつないだ第２の直線と、の交点を含む所定期間の制御データを平滑化する。【選択図】図１
2014002986.tif 000002
【特許請求の範囲】
【請求項１】
加速器の加速空洞を制御するための高周波制御装置であって、
前記加速空洞を駆動するための高周波信号のパターンデータを、前記加速器の加速サイクルにおける、フラットボトム期間、加速期間、フラットトップ期間、減速期間を含む複数の制御期間毎に、時間に対する一次関数を定義するデータで記憶する記憶部と、
前記記憶部に記憶させたパターンデータを読み込み、読み込んだパターンデータに基づいて、前記加速空洞を逐次制御するための制御データを生成する制御データ生成部と、
前記制御データに基づいて、前記加速空洞を逐次制御する制御部と、を備え、
前記制御データ生成部は、前記複数の制御期間に応じた制御データを、それぞれ時間に対する一次関数上のデータとして生成するとともに、
生成した制御データのうち、第１の制御期間の制御データを時間軸上でつないだ第１の直線と、第２の制御期間の制御データを前記時間軸上でつないだ第２の直線と、の交点を含む所定期間の制御データを平滑化することを特徴とする高周波制御装置。
【請求項２】
前記平滑化した制御データを前記時間軸上でつないだ曲線が、前記第１の直線および前記第２の直線の少なくとも一方と接線を共有するように、前記平滑化を行うことを特徴とする請求項１に記載の高周波制御装置。
【請求項３】
前記制御データ生成部は、前記所定期間の制御データの時間当たりの設定回数が、他の期間の制御データの時間当たりの設定回数よりも多くなるように、前記制御データを生成することを特徴とする請求項１または２に記載の高周波制御装置。
【請求項４】
前記記憶部が記憶するデータは、前記制御期間毎の前記高周波信号の周波数と振幅の少なくともいずれかの初期値と、当該制御期間の長さであることを特徴とする請求項１ないし３のいずれか１項に記載の高周波制御装置。
【請求項５】
前記制御部は、
前記加速空洞に出力された高周波信号の位相、前記加速空洞に印加された電圧、および前記加速器内の粒子線の位置情報の少なくともいずれかのフィードバック信号を用いたフィードバック制御により、前記加速空洞の制御を行うことを特徴とする請求項１ないし４のいずれか１項に記載の高周波制御装置。
【請求項６】
前記制御部は、
前記フラットトップ期間において、前記フィードバック信号が予め定められている範囲を超えると、当該フィードバック信号を前記範囲内に入るように補正することを特徴とする請求項５に記載の高周波制御装置。
【請求項７】
粒子線を加速するための加速空洞が設けられた加速器と、
前記加速空洞を制御する請求項１ないし６のいずれか１項に記載の高周波制御装置と、
前記加速器によりエネルギーを調整されて供給された粒子線を、照射対象に応じた照射野に成形して前記照射対象へ照射する照射装置と、
を備えたことを特徴とする粒子線治療装置。
【請求項８】
前記照射装置は、スキャニング方式によって前記照射野を形成することを特徴とする請求項７に記載の粒子線治療装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、粒子線治療装置の線源であるシンクロトロンのような加速器用の高周波制御装置および粒子線治療装置に関するものである。
【背景技術】
【０００２】
粒子線治療は、患部組織に粒子線を照射してダメージを与えることで治療を行うものであり、広義の放射線治療のひとつである。しかし、陽子線や重イオン線等の粒子線は、γ線、Ｘ線といった他の放射線と異なり、粒子線のエネルギーによって線量付与の深度範囲を調整することができ、患部の立体形状に応じた線量付与が可能となる。そのため、とくに、粒子線治療装置用の加速器には、正確なエネルギーや軌道で粒子線を供給することが求められている。
【０００３】
加速器は、周回軌道を形成するための偏向電磁石、高周波を使って粒子線を加速する加速空洞、および粒子線が通過するための通路となる真空ダクトから構成される。偏向電磁石の磁場は粒子線の加速（エネルギー増大）に伴い、定められたパターンに従って変化する。同時に粒子線の周回周波数も変化するため、安定に加速するためには上述した加速空洞に印加される高周波信号も定められたパターンによって周波数を制御する必要がある。
【０００４】
このような制御を単純化するため、例えば、計算機であらかじめ作成した高周波信号の周波数リファレンス信号を記憶しておき、運転中に順次読みだして制御するＲＦ制御装置（例えば、特許文献１参照。）が提案されている。しかし、このような方式では、データ量が膨大となり、データ調整が容易でなくなるとともに、構成要素も多大となる。そこで制御に使用する波形データを加速器の運転パターンを平坦な安定領域と平坦な加速部分を持つ領域とに分け、領域ごとのデータを用いて制御するようにした高周波加速制御装置（例えば、特許文献２参照。）が提案されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０００−２３２０００号公報（段落００３１〜００５０、図１）
【特許文献２】特開２０１０−３５３８号公報（段落００１８〜００２４、図２〜図５）
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、上述したように、フラットな部分（特許文献２での平坦な安定領域に対応）と単純な勾配部分（同平坦な加速部分に対応）を単純に切り替えて制御した場合、周波数の変化率が急変する部分が生じ、突変と呼ばれる現象が発生して、加速器を適切に制御できない場合があることが分かった。
【０００７】
本発明は、上記のような課題を解決するためになされたもので、加速器の運転および運転調整時における運転パターンのパラメータ調整が容易であるとともに、加速器を適切に制御できる高周波制御装置、およびそれを用いて正確な治療が可能な粒子線治療装置を得ることである。
【課題を解決するための手段】
【０００８】
本発明の高周波制御装置は、加速器の加速空洞を制御するための高周波制御装置であって、前記加速空洞を駆動するための高周波信号のパターンデータを、前記加速器の加速サイクルにおける、フラットボトム期間、加速期間、フラットトップ期間、減速期間を含む複数の制御期間毎に、時間に対する一次関数を定義するデータで記憶する記憶部と、前記記憶部に記憶させたパターンデータを読み込み、読み込んだパターンデータに基づいて、前記加速空洞を逐次制御するための制御データを生成する制御データ生成部と、前記制御データに基づいて、前記加速空洞を逐次制御する制御部と、を備え、前記制御データ生成部は、前記複数の制御期間に応じた制御データを、それぞれ時間に対する一次関数上のデータとして生成するとともに、生成した制御データのうち、第１の制御期間の制御データを時間軸上でつないだ第１の直線と、第２の制御期間の制御データを前記時間軸上でつないだ第２の直線と、の交点を含む所定期間の制御データを平滑化することを特徴とする。
【発明の効果】
【０００９】
本発明の高周波制御装置によれば、１次関数を定義する単純なデータを用いて逐次制御するための制御データを生成するとともに、所定期間の制御データを平滑化するようにしたので、加速器運転および運転調整時における運転パターンのパラメータ調整が容易であるとともに、加速器を適切に制御できる。さらに、本高周波制御装置を用いることによって、正確な照射が可能な粒子線治療装置を得ることができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施の形態１にかかる高周波制御装置の構成を示すブロック図である。
【図２】本発明の実施の形態１にかかる高周波制御装置の制御における加速サイクルでの区分制御期間を説明するための図である。
【図３】本発明の実施の形態１にかかる高周波制御装置の制御におけるデータ構造を説明するための図である。
【図４】本発明の実施の形態１にかかる高周波制御装置の制御データを生成する際のカウント動作を説明するための図である。
【図５】本発明の実施の形態１にかかる高周波制御装置の制御データのスムージング処理動作を説明するための図である。
【図６】本発明の実施の形態２にかかる高周波制御装置の構成を示すブロック図である。
【図７】本発明の実施の形態２にかかる高周波制御装置の制御におけるフィードバック補正を説明するための図である。
【図８】本発明の実施の形態３にかかる高周波制御装置の構成を示すブロック図である。
【図９】本発明の実施の形態４にかかる高周波制御装置の構成を示すブロック図である。
【図１０】本発明の実施の形態５にかかる粒子線治療装置の構成を示す図である。
【発明を実施するための形態】
【００１１】
実施の形態１．
以下、本発明の実施の形態１にかかる高周波制御装置の構成について説明する。図１〜図５は本発明の実施の形態１にかかる高周波制御装置について説明するためのもので、図１は高周波制御装置の構成を説明するための制御対象である加速器を含めた高周波制御装置のブロック図、図２は加速器の加速サイクルで、制御を行うためのデータの切替区分である区分制御期間を説明するための図、図３は制御データを生成するために記憶（読込）される運転パターンごとのデータ構造を説明するための図、図４は加速サイクルでの区分制御期間ごとの制御データを生成するためのタイミング制御部におけるカウント動作を説明するための図、図５は区分制御期間の内、スムージング対象となる期間における制御データのスムージング処理動作を説明するための図である。以下、詳細に説明する。
【００１２】
高周波制御装置３は、図１に示すように、シンクロトロンのような環状加速器（以下、単に加速器６と称する。）には、主構成要素である加速空洞６２を制御するためのものである。加速器６は、周回軌道を形成するための偏向電磁石６１、高周波を使って粒子線を加速する前述の加速空洞６２、および粒子線が通過するための通路となる真空ダクト６３が設けられている。これ以外の機器としては、粒子線を前もって加速する前段加速器６６（後述）、加速された粒子線を真空ダクト６３に入射する入射装置６５、粒子線（ビーム）の位置等を計測するビームモニタ６４、粒子線を加速器６からユーザ側に送り出す出射装置６７などから構成される。そして、加速器６の偏向磁界の変化を示すフラットベース期間、加速期間、フラットトップ期間及び減速期間に対応する加速空洞６２を駆動するための高周波信号の高周波信号制御データは、予め計算機１９で作成されているとして説明する。
【００１３】
高周波制御装置３は、上述した高周波信号制御データを保存するメモリモジュール３１と、基準クロック３６と、Ｄ／Ａコンバータに内蔵されたＤＤＳ５０（ダイレクトデジタルシンセサイザ）と、振幅変調波生成器５１と、振幅変調波生成器５１とＤＤＳ５０の出力信号を合成する乗算器５２と、乗算器５２からの出力信号のノイズを低減するローパスフィルタ５３と、メモリモジュール３１から読み出した高周波信号制御データからＤＤＳ５０へ設定する周波数設定データと振幅変調波生成器５１に設定する振幅設定データを演算生成する波形演算制御部３３と、基準クロック３６からのクロック信号により波形演算制御部３３に加速空洞６２を駆動するための高周波信号のための運転パターンの制御タイミングを与えるタイミング制御部３２とを備えている。そして、イベントカウンタ４２が設けられたタイミング制御部３２と、データメモリ４５が内蔵された波形演算制御部３３とで、デジタル制御回路３０を構成している。
【００１４】
つぎに、加速器６の運転に伴う高周波制御装置３の制御について説明する。なお、ここでは、一般的な制御部分から説明し、その後に本発明の実施の形態にかかる高周波制御装置３の特徴的な制御について説明する。
【００１５】
加速器６の偏向電磁石６１の磁場は、後述するように、一定のパターンに従って変化する。このとき、同時に粒子線の周回周波数も変化するため、安定に加速するためには加速空洞２に印加される高周波信号も、図２の２段目に示すように、定められたパターンによって周波数を制御する必要がある。
【００１６】
この制御において、加速中に粒子線は、シンクロトロン振動と呼ばれる現象を起こすため、そのままでは粒子線に対して何らかの外乱があった場合にビームが不安定になる。この対策として図１に示すように、位相モニタ回路１４により、検出された位相制御フィードバック信号Ｓ２０と加速用の高周波信号の位相差に基づき高周波制御装置３にて周波数のフィードバック制御を行い、シンクロトロン振動を減衰する方策が通常とられている。
【００１７】
ここで、加速空洞６２の周波数が偏向電磁石６１の磁場で決まる最適周波数からずれた場合には、粒子線の位置が理想的な中心軌道からずれるため、ビームモニタ６４、および位置モニタ回路１５でこれを検出し、検出した位置制御フィードバック信号Ｓ２１を用いて高周波制御装置３でフィードバック制御することも行われる。
【００１８】
また、加速空洞６２に実際印加される高周波電圧を電圧モニタ回路１３にて検出し、検出した電圧制御フィードバック信号Ｓ２２に基づき高周波制御装置３にて高周波信号の振幅をフィードバック制御することも行われる。
【００１９】
これらの制御において、デジタル制御回路３０は、メモリモジュール３１からの高周波信号制御データ（パターンデータ）と前記フィードバック制御にてＤＤＳ５０と振幅変調波生成器５１へ入力する設定データを生成する。
【００２０】
また、高周波制御装置３の外部構成として、変調出力を増幅するパワーアンプ１２、上述した計算機１９の制御により、メモリモジュール３１のデータ読み出しタイミングを制御する図示しないタイミング発生器等を備えることがある。
【００２１】
次に、加速器６の運転パターンについてさらに説明する。加速器６では、粒子線の加速に伴い、図２の上段に示すように、一般的にフラットボトムと称される一定値を保つ期間、一定勾配で加速する加速期間、一般的にフラットトップと称される一定値を保つ期間、および一定勾配で減速する減速期間からなる４つの期間を１サイクルとして、周期的な制御が繰り返される。
【００２２】
フラットボトム期間は、加速器６の最低エネルギーの状態である。フラットボトムでは前段加速器６６（図１０参照）から粒子線を入射し、加速器６内に蓄積する。次にＯＦＦ状態にあった高周波信号をＯＮとし、電圧を増加させることにより粒子線が高周波信号の一定の位相にあつまるように捕獲を行う。
【００２３】
次に、加速期間では偏向電磁石６１の磁場変動に基づいて高周波信号の周波数と電圧を制御して粒子線の加速を行う。なお、加速の開始部分と終了部分は区別してスムージングと呼ばれ、後述するように、特別なパターンが必要となる。
【００２４】
そして、粒子線が所定のエネルギーに達した時点で加速を終了し、磁場が変化しないフラットトップの状態に入る。この期間はフラットトップ期間と呼ばれ、粒子線を加速器６から出射するのに適した条件を作るため、高周波信号の振幅、周波数などの微調整を行う。
【００２５】
次に、加速器６の磁場を最低値まで下げるように減速を行う減速期間に入る。減速の開始部分と終了部分も区別してスムージングと呼ばれている。減速期間が終了すると、加速器６は、再びフラットボトム期間に入り、電圧を初期状態に戻す。上記の運転パターンの周期を加速サイクルと呼んでいる。
【００２６】
次に、本実施の形態にかかる加速器用高周波制御装置３でのデータ設定および動作について説明する。上述したように、加速器６の運転開始前には、予め計画された全ての加速空洞駆動用の高周波信号制御データが、計算機１９によって設定されており、設定されたデータは、随時読込みできるよう、メモリモジュール３１に保存されている。そして、このデータの設定方法と、その設定に従っての制御動作が本発明の実施の形態にかかる加速器用高周波制御装置３の特徴となる。
【００２７】
高周波信号制御データは、図２に示す加速空洞駆動の運転パターンの１周期分を示す加速サイクルを、上述したフラットボトム、加速、フラットトップ、減速、フラットボトムと各期間を接続するスムージング期間の計９に分割した区分制御期間Ｔ１〜Ｔ９毎に設定する。設定内容は、区分制御期間毎の、時間間隔、高周波信号の周波数設定データの初期値（図２の周波数制御の段の○印）、高周波信号の振幅設定データの初期値（図２の振幅制御の段の○印）、その期間における設定値の分解能に対応する区間分割数、および、設定値の修正（スムージング）の有無で構成される。なお、区間分割数とは、区分制御期間毎のＤＤＳ５０用の周波数設定データと振幅変調波生成器５１用の振幅設定データのデータ設定回数である。
【００２８】
上述したデータは、図３に示す様に、運転パターン毎にフォルダ構成を持つように、メモリモジュール３１に保存されている。その構成はＦＩＦＯ（先入れ先出し）となっており、運転パターンの１番目からＮ番目までを波形演算制御部３３からの読み出し要求に従って順番にフォルダ情報を出力する。運転パターンのＮ番目とは加速器６の運転のための最終加速サイクルの設定データを意味する。
【００２９】
加速器６の制御は、計算機１９からタイミング信号である加速サイクル開始信号Ｓ０が出力されることにより、波形演算制御部３３がメモリモジュール３１に保存された高周波信号制御データを読み込むことで開始される。そして、波形演算制御部３３は、メモリモジュール３１から読み込んだ高周波信号制御データに基づいてデータ演算し、ＤＤＳ５０と振幅変調波生成器５１にデータ設定を行う。以下、さらに詳細に説明する。
【００３０】
計算機１９は加速サイクルを開始するため、まずタイミング制御部３２に最初のデータ設定トリガ信号Ｓ１ａを出力する。タイミング制御部３２は図２の下側の段に示したように、最初のデータ設定トリガ信号Ｓ１ａのみ波形演算制御部３３へそのまま信号出力し、２回目以降のデータ設定トリガ信号Ｓ１ｂは、各加速サイクルの最後の区分制御期間Ｔ９に同期して波形演算制御部３３へ出力される。
【００３１】
波形演算制御部３３は、データ設定トリガ信号Ｓ１ａ、Ｓ１ｂが入力される毎に加速サイクルにおける高周波信号制御データを順次メモリモジュール３１から読み込む処理を行う。波形演算制御部３３は読み込んだ高周波信号制御データから区分制御期間Ｔ１ないしＴ９毎の時間間隔データと設定分解数データをタイミング制御部３２内のイベントカウンタ４２に毎回プリセットする。
【００３２】
タイミング制御部３２のイベントカウンタ４２は、区分制御期間（Ｔ１ないしＴ９）毎のレジスタを持っている。そして、計算機１９から加速サイクル開始信号Ｓ０が入力されると、カウントアップを開始し、基準クロック３６のクロック信号を基準として、図４に示すように区分制御期間毎にプリセットされた時間間隔に到達するまでカウントを継続する。カウントがプリセット値に到達した時点でカウントオーバー信号Ｓ４を波形演算制御部３３へ出力する。
【００３３】
また、タイミング制御部３２は区分制御期間の区間分割数に合せて時間間隔の均等分割処理を行い、それぞれの分割された時間にカウント値が到達した時点で波形演算制御部３３へイベントトリガ信号Ｓ３を出力する。例えば、図４において区分制御期間Ｔ２の時間間隔が１５０ｍｓで区間分割数が１５の場合には、イベントトリガ信号Ｓ３は１０ｍｓ間隔（１５０ｍｓ／１５）で波形演算制御部３３へ出力される。
【００３４】
次に、波形演算制御部３３は、メモリモジュール３１から高周波信号制御データを読み込んだ際に区分制御期間毎の高周波信号の周波数設定データの初期値、振幅設定データの初期値、周波数設定データと振幅設定データの区間分割数からイベントトリガ信号Ｓ３に同期して設定するＤＤＳ５０用の周波数設定データと振幅変調波生成器５１用の振幅設定データのデータ変化量を予め演算し、内蔵データメモリ４５に区分制御期間毎の変化量を保存する。
【００３５】
区分制御期間毎のＤＤＳ５０用の周波数設定データのデータ変化量は、当該区分制御期間での周波数設定データの初期値と、次の区分制御期間での周波数設定データの初期値、および当該区分制御期間での区分分割数から演算する。例えば、図４において、区分制御期間Ｔ２の変化量は、区分制御期間Ｔ２での周波数設定データの初期値１０ｋＨｚと、区分制御期間Ｔ３での周波数設定データの初期値２５ｋＨｚの差分１５ｋＨｚを、区分分割数１５で除した１ｋＨｚとなる。振幅変調波生成器５１用の振幅設定データの変化量についても、同様に演算される。
【００３６】
波形演算制御部３３はタイミング制御部３２からイベントトリガ信号Ｓ３が入力される度に、ＤＤＳ５０用の周波数設定データと振幅変調波生成器５１用の振幅設定データにデータ変化量を加算していく。このとき、区分制御期間が加速の場合は変化量は正となり、一定勾配の増加、減速の場合は負となり一定勾配の減少となる。そして、フラットボトムおよびフラットトップでは、変化量は０となり、一定値となる。そして、タイミング制御部３２からカウントオーバー信号Ｓ４が入力されると内蔵データメモリ４５の区分制御期間番号を切り替えて設定するデータ変化量も変更する。つまり、加速サイクルにおける４大工程であるフラットボトム（Ｔ１、Ｔ９）、加速（Ｔ３）、フラットトップ（Ｔ５）、減速（Ｔ７）では、工程（区分制御期間）毎の設定値を当該区間の設定値と次の区間の設定値、および区分分割数の３つのデータの組み合わせた値で設定でき、その値は直線上にのる。
【００３７】
一方、上記のように直線上の設定値で制御する工程間を直接接続して加速サイクルを形成した場合、接続部分において変化率がステップ状に急変し、突変が生じて設定通りに制御できない場合がある。そこで、本実施の形態にかかる高周波制御装置３では、４大工程間のスムージング期間と呼ばれる期間（区分制御期間Ｔ２、Ｔ４、Ｔ６、Ｔ８）については、４大工程と同様に、直線状の設定値を計算した上で、計算した設定値を修正する修正処理（スムージング処理）を行うようにした。修正値は、当該区分制御期間の暫定設定値と、前後の工程における設定値とを比較して設定される。
【００３８】
スムージング処理対象となる期間（Ｔ２、Ｔ４、Ｔ６、Ｔ８）においても、前述した区分制御期間毎に演算で得られる一定の変化量に応じた直線状の設定値を計算する。このとき、図３のデータ形式中の修正有無の項が、「有」と設定されていると、計算した設定値を暫定設定値とし、暫定設定値に対して修正値による修正を行う。修正値の設定方法およびスムージング処理について、区分制御期間Ｔ２を例に、図５を用いて説明する。
【００３９】
図５に示すように、フラットボトム（区分制御期間Ｔ１）で設定された一定値の設定値（Ｌ（Ｔ１）：破線）と、スムージング対象の区分制御期間Ｔ２で設定された暫定設定値（Ｌ（Ｔ２）：破線）との接続点１−２では、変化率が急変している。また、暫定設定値Ｌ（Ｔ２）と、加速（区分制御期間Ｔ３）で設定された一定勾配の設定値（Ｌ（Ｔ３）：破線）との接続点２−３でも、変化率が急変している。そこで、各接続点において、変化率が急変しないよう、例えば、フラットボトムＴ１の設定値Ｌ（Ｔ１）、および加速Ｔ２の設定値Ｌ（Ｔ３）と接線を共有する曲線（実線）を描けるような設定分割数に応じた修正データを作成する。そして、計算した暫定設定値に作成した修正データを加算することで、区分制御期間Ｔ２の修正後の設定値Ｓ（Ｔ２）を得ることができる。これにより、工程間において変化率が急変することがなく、突変の発生を抑止することができる。
【００４０】
このような、区間分割を示す各イベントトリガ信号Ｓ３に対応した修正データは、予め計算機１９で計算しておき、メモリー等に保存しておくように、あるいは、波形演算制御部３３に設定するようにすれば、制御中に余分な演算をする必要がなくなる。また、修正データを適用するか否かは、上述したようにデータ中に設定していてもよいが、計算機１９から修正データを使用するための信号を出すようにしてもよい。
【００４１】
波形演算制御部３３はイベントトリガ信号Ｓ３に同期して修正データを順次ＤＤＳ５０用の周波数設定データと振幅変調波生成器５１用の振幅設定データに加算することでスムージング処理を行う。スムージング対象期間の内、区分制御期間Ｔ２とＴ８においては暫定設定値のマイナス側に修正データを反映し、区分制御期間Ｔ４とＴ６においては暫定設定値のプラス側に修正データを反映することで、変化率の急変を抑えた設定値のカーブを得ることができる。なお、突変を抑制するためのスムージングとして、変化率を急変させないように接線を共有する曲線を得る例について説明したが、これに限ることはない。加速器６のパラメータ設定（運転調整）等において、突変を生じないような設定値が得られるのであれば、各イベントトリガ信号Ｓ３に対応した修正データは適宜調整（生成）すればよい。
【００４２】
なお、上述した例では、各区分制御期間の制御値は、３つの値で設定しているが、区間分割数はデータ設定の分解能を示すものであって、設定値自体は、当該区分制御期間の設定値と次の区間の設定値のみで決まる。つまり、各区分制御期間の高周波信号制御データとして、時間に対する一次関数を定義するための加速サイクル内の各区分制御期間の時間と初期値を保存する。それに加えて、スムージング処理の有無、スムージング対象となる４大工程をつなぐ遷移期間（Ｔ２、Ｔ４、Ｔ６、Ｔ８）での修正データさえ保持していればよい。これだけのデータにより、遷移期間をスムージング処理でき、加速や減速時の高周波信号パターンの突変対策や不要ノイズの軽減を図り、加速と減速時のビーム捕捉の安定性を高めることができる。
【００４３】
一方、従来のように、工程に関係なく、一定間隔でＤＤＳと振幅変調波生成器に設定するデータを逐次演算する場合には、調整要素が多く、メモリモジュールに保存するデータ量も多い。即ち、加速器の立ち上げ調整をする時点では最適な運転方法やパラメータが事前にわかっておらず、実際にシステムを運転してみながらパラメータ調整を行う場合が多いため、加速器の各施工場所において最適な加速運転パターンのパラメータ調整の単純化が図れないという問題があった。そのため、加速器の運転中にリアルタイムで信号波形を変化させるなどの柔軟性が乏しく、粒子線の加速中にフィードバック情報を用いてリアルタイムに高周波出力信号波形を自動調整し加速器の安定化を図る機能が不足しているという問題点があった。
【００４４】
それに対して、例えば、４大工程ごとに生成した高周波信号パターンを保存しておくようにするだけでは、工程間の接続部分（スムージング期間）で、変化量の急変により、突変が生じる可能性があり、治療計画に基づく治療の前提となるエネルギーと軌道が正確な粒子線を出射することが困難となる。
【００４５】
しかしながら、本実施の形態にかかる高周波制御装置３では、加速器６の運転中の加速空洞駆動制御において高周波信号の周波数と振幅の設定パラメータを加速サイクルの区分制御期間Ｔ１ないしＴ９毎に演算して加速運転パターンを生成することができるため、計算機１９からメモリモジュール３１に入力する高周波信号制御データを単純化することが可能となる。また、運転試験時のパラメータ調整においても調整要素が少なく微調整ができるという柔軟な対応が可能となる。
【００４６】
なお、上記例では、遷移期間に対して、別途暫定設定値を設けて修正を行う例について説明したが、これに限ることはない。例えば、スムージング対象期間Ｔ２部分にフラットボトムＴ１の設定値Ｌ（Ｔ１）の直線と加速Ｔ３の設定値Ｌ（Ｔ３）の直線を延長させて暫定設定値を算出する。そして、設定値Ｌ（Ｔ１）の直線と加速Ｔ３の設定値Ｌ（Ｔ３）の直線の接続点を含む変化量が急変する部分をスムージングするように修正データを生成してもよい。
【００４７】
以上のように、本実施の形態１にかかる高周波制御装置３によれば、加速器６の加速空洞６２を制御するための高周波制御装置３であって、加速空洞６２を駆動するための高周波信号のパターンデータを、加速器６の加速サイクルにおける、フラットボトム期間Ｔ１，Ｔ９、加速期間Ｔ３、フラットトップ期間Ｔ５、減速期間Ｔ７を含む複数の制御期間（区分制御期間）毎に、時間に対する一次関数を定義するデータで記憶する記憶部として機能するメモリモジュール３１と、記憶部（メモリモジュール３１）に記憶させたパターンデータを読み込み、読み込んだパターンデータに基づいて、加速空洞６２を逐次制御するための制御データを生成する制御データ生成部として機能する波形演算制御部３３と、制御データに基づいて、加速空洞６２を逐次制御する制御部として機能するデジタル制御回路３０、ＤＳＳ５０、振幅変調波生成器５１と、を備え、制御データ生成部（波形演算制御部３３）は、複数の制御期間に応じた制御データを、それぞれ時間に対する一次関数上のデータとして生成するとともに、生成した制御データのうち、第１の制御期間の制御データを時間軸上でつないだ第１の直線と、第１の制御期間に続く第２の制御期間の制御データを時間軸上でつないだ第２の直線と、の交点を含む所定期間の制御データを平滑化するように構成したので、加速器６の運転および運転調整時における運転パターンのパラメータ調整が容易であるとともに、加速器６を適切に制御できる。
【００４８】
とくに、制御データ生成部（波形演算制御部３３）は、平滑化した制御データを時間軸上でつないだ曲線が、第１の直線および第２の直線の少なくとも一方と接線を共有するように、平滑化を行うようにしたので、制御データがスムーズに変化し、突変の発生を抑えて、加速器６を適切に制御することができる。
【００４９】
さらに、制御データ生成部（波形演算制御部３３）は、スムージング期間と呼ばれる所定期間（区分制御期間Ｔ２，Ｔ４，Ｔ６，Ｔ８）の制御データの時間当たりの設定回数が、他の期間（区分制御期間Ｔ１，Ｔ３，Ｔ５，Ｔ７，Ｔ９）の制御データの時間当たりの設定回数よりも多くなるように、制御データを生成するようにしたので、制御データの変化が一層スムーズになり、効果的に突変の発生を抑えて、加速器６を適切に制御することができる。
【００５０】
また、記憶部（メモリモジュール３１）に記憶させる一次関数を定義するデータは、制御期間（区分制御期間）毎の高周波信号の周波数と振幅の少なくともいずれかの初期値と、当該制御期間の時間であるようにしたので、単純で直観的なイメージに沿って容易にパラメータ調整ができる。
【００５１】
また、制御部（デジタル制御回路３０）は、加速空洞６２に出力された高周波信号の位相（Ｓ２０）、加速空洞６２に印加された電圧（Ｓ２２）、および加速器６内の粒子線の位置情報（Ｓ２１）の少なくともいずれかのフィードバック信号を用いたフィードバック制御により、加速空洞６２の制御を行うようにしたので、加速器６を適切に制御できる。
【００５２】
実施の形態２．
本実施の形態２では、実施の形態１で示した高周波制御装置に対し、フィードバック制御信号を補正する補正処理部を追加したものである。図６と図７は本発明の実施の形態２にかかる高周波制御装置について説明するためのもので、図６は高周波制御装置の構成を説明するためのブロック図、図７は補正部によるフィードバック補正を説明するための図である。なお、図６において、図１との共通部分である高周波制御装置の外側の部分については、記載を省略している。
【００５３】
図６に示すように、本実施の形態２にかかる高周波制御装置３では、波形演算制御部３３とフィードバック制御用のモニタ回路１３、１４、１５間に補正処理部３４を設けた。そして、加速器６の運転中の粒子線加速状態を示すフィードバックモニタ信号である位相制御フィードバック信号Ｓ２０、位置制御フィードバック信号Ｓ２１、電圧制御フィードバック信号Ｓ２２を補正処理部３４を介して波形演算制御部３３に入力されるようにした。
【００５４】
次に補正処理部３４の動作について説明する。
補正処理部３４は各フィードバックモニタ信号Ｓ２１、Ｓ２２、Ｓ２３で得られたパラメータを基に波形演算制御部３３が設定するＤＤＳ５０用の周波数設定データと振幅変調波生成器５１用の振幅設定データを補正する目的を持つ。補正処理部３４は加速サイクルのフラットトップの期間、各フィードバックモニタ信号の信号品位を確保するため平滑化処理によるノイズ除去を行う機能を持つ。
【００５５】
平滑化処理にはＮ回のサンプリングから得られた測定値を積算平均する処理を用い、サンプリング数（Ｎ）は、計算機１９から設定および変更が可能である。また、平滑化処理機能のオンオフは、波形演算制御部３３が出力する補正ゲート信号Ｓ７により制御され、例えば、補正ゲート信号Ｓ７がアサートされている期間だけ平滑化処理機能をイネーブルにする。平滑化処理機能がオンの場合、補正処理部３４は平滑化した各フィードバックモニタ信号から補正すべき対象が周波数であるのか振幅であるのか、あるいは両方であるのかの補正判定信号Ｓ６をタイミング制御部３２から入力されるイベントトリガ信号Ｓ３に同期して波形演算制御部３３へ通知する。
【００５６】
補正判定信号Ｓ６は４ビットの信号であり、２ビットは周波数用、２ビットは振幅用に分けられる。それぞれの２ビット値が００および１０では補正なし、０１では後述する下限しきい値に対応した補正が必要で、１１では後述する上限しきい値に対応した補正が必要であることの意味を持つ。なお、位相制御フィードバック信号Ｓ２０または位置制御フィードバック信号Ｓ２１に関しては周波数補正に用い、電圧制御フィードバック信号Ｓ２２に関しては振幅補正を用いる。
【００５７】
補正判定方法の具体例について図７を用いて説明する。図に示すように、フィードバック補正判定には位相制御フィードバック信号Ｓ２０、位置制御フィードバック信号Ｓ２１、電圧制御フィードバック信号Ｓ２２の各フィードバックモニタ信号に対応した補正パラメータとして上限しきい値（図中上限値と表記）、下限しきい値（同下限値）およびこれらの判定基準となるフィードバック基準値（同ＦＢ基準値）を用いる。これらの補正パラメータは計算機１９により、補正処理部３４に予め設定しておき変更することも可能である。
【００５８】
補正処理部３４は平滑化した位相制御フィードバック信号Ｓ２０、位置制御フィードバック信号Ｓ２１、電圧制御フィードバック信号Ｓ２２の各フィードバックモニタ信号のレベルとフィードバック基準値を比較し、その差分が上限しきい値または下限しきい値を超えないか監視する。しきい値の範囲を超えた場合は直ちに波形演算制御部３３へ補正判定信号Ｓ６で通知する。なお、補正処理部３４は計算機１９から起動と停止の設定が可能である。
【００５９】
次に、波形演算制御部３３は補正判定信号Ｓ６の４ビットの値から、補正対象が周波数か振幅かを知ることができる。そして該当する補正対象を補正トリガ信号Ｓ９のタイミングでＤＤＳ５０用の周波数設定データと振幅変調波生成器５１用の振幅設定データにフィードバック補正値を加算して補正処理を行う。フィードバック補正値は加速器６の運転開始前に計算機１９より予め周波数と振幅の変化量として設定される。
【００６０】
図７では、上限しきい値対応分のフィードバック補正値を補正Ａ、下限しきい値対応分のフィードバック補正値を補正Ｂとした場合である。図において、フィードバックモニタ信号（図中ＦＢ値と表記）のレベルが、点ＰＡにおいてフィードバック上限しきい値を超えたため、ＤＤＳ５０用の周波数設定データあるいは振幅変調波生成器５１用の振幅設定データに補正Ａを加えてフィードバックモニタ信号のレベルをフィードバック基準値に近付ける例を示している。また点ＰＢにおいて、フィードバックモニタ信号のレベルがフィードバック下限しきい値を超えたためＤＤＳ５０用の周波数設定データあるいは振幅変調波生成器５１用の振幅設定データに補正Ｂを加えてフィードバックモニタ信号のレベルをフィードバック基準値に近付ける例を示している。
【００６１】
こうして、加速器６の運転中の粒子線加速状態を監視し、しきい値を超えると直ちに補正処理部３４から波形演算制御部３３へ情報を与えて、各フィードバックモニタ信号がフィードバック基準値に近づくよう引き込むトラッキング制御を行うことで、安定した加速サイクルのフラットトップを得ることができる。これにより、例えば、図２で示した加速サイクルの区分制御期間Ｔ５であるフラットトップの期間において、加速器６の運転中の加速空洞駆動制御の調整をリアルタイムに行い、粒子線加速の安定化を図ることができる。また、各フィードバックモニタ信号の外乱ノイズの影響を低減することができる。また補正処理部３４の平滑化平均回数、フィードバック補正判定用パラメータ、フィードバック補正値の各パラメータを計算機１９で調整設定できることで、シンクロトロンの試験運転時に柔軟に対応ができ最適なパラメータ調整ができる効果がある。
【００６２】
以上のように、本実施の形態２にかかる高周波制御装置３によれば、制御部（デジタル制御回路３０）は、フラットトップ期間Ｔ５において、フィードバック信号が予め定められている範囲を超えると、当該フィードバック信号を範囲内に入るように補正する補正部（補正処理部３４）を設けるようにしたので、安定した加速サイクルのフラットトップを得ることができる。
【００６３】
実施の形態３．
本実施の形態３では、実施の形態１と実施の形態２で示したデジタルダイレクトシンセサイザの代わりに、Ｄ／Ａコンバータ内蔵のマルチ出力ＤＤＳとマルチプレクサを備えるようにしたものである。その他の構成については、上記実施の形態１や２と同様であるので説明を省略する。
【００６４】
図８は、本発明の実施の形態３にかかる高周波制御装置の構成を説明するためのブロック図であり、基本的には実施の形態２の説明で用いた図６に対応するものである。図に示すように、本実施の形態３においても、実施の形態１あるいは実施の形態２で記載したデジタル制御回路３０の機能を有する。そして、制御のための信号を出力する部分において、ＤＤＳ５０の代わりに、Ｄ／Ａコンバータ内蔵のマルチ出力ＤＤＳ５４とマルチプレクサ５５を具備するようにした。波形演算制御部３３によりマルチプレクサ５５の入力ポートを制御することで３チャンネルのマルチ出力ＤＤＳ５４の出力信号を切り替える構成である。
【００６５】
次に実施の形態３にかかる高周波制御装置の動作について説明する。
実施の形態２での周波数補正では、フィードバックモニタ信号（Ｓ２０、Ｓ２１、Ｓ２２）のレベルとフィードバック基準値との比較結果により、ＤＤＳ５０用の周波数設定データにフィードバック補正値を加算する処理を行った例を示した。本実施の形態３ではフィードバック補正値の加算に代えて、マルチ出力ＤＤＳ５４の信号出力を直接変更するようにした。
【００６６】
マルチ出力ＤＤＳ５４は、周波数設定値としてメモリモジュール３１から波形演算制御部３３が読み出した高周波信号制御データによって与えられる周波数設定値とそれを基準値として高めの周波数と低めの周波数の設定値が設定される。なお、高めの周波数と低めの周波数の設定値は基準値からの差分パラメータとして加速器６の運転開始前に計算機１９より予め波形演算制御部３３に設定される。
【００６７】
波形演算制御部３３は、補正判定信号Ｓ６の４ビットの値から周波数に関する２ビットの状態により、補正トリガ信号Ｓ９のタイミングでマルチプレクサ５５の入力ポートを周波数選択信号Ｓ８により切り替える。補正判定信号Ｓ６の２ビットが００または１０ならば基準値設定を選択、０１ならば高めの周波数設定を選択、１１ならば低めの周波数設定値を選択する。
【００６８】
こうして、加速器６の運転中の粒子線加速状態を監視し、しきい値を超えると直ちに補正処理部３４から波形演算制御部３３へ情報を与える。これにより、位相制御フィードバック信号Ｓ２０と位置制御フィードバック信号Ｓ２１の各フィードバックモニタ信号がフィードバック基準値に近づくよう引き込むトラッキング制御を行い、安定した加速サイクルのフラットトップを得ることができる。また高めの周波数と低めの周波数の設定値は計算機１９で調整設定できることで、加速器６の試験運転時に柔軟に対応ができ最適なパラメータ調整ができる効果がある。
【００６９】
なお、周波数選択信号Ｓ８を制御してマルチプレクサ５５の入力ポートを任意に切り替えることで乗算器５２に入力する信号を周波数変調することも可能である。
【００７０】
以上のように、本実施の形態３にかかる高周波制御装置３によれば、制御部（５１〜５５部分）は、フラットトップ期間Ｔ５において、フィードバック信号が予め定められている範囲を超えると、当該フィードバック信号を範囲内に入るように補正する補正部として機能するマルチ出力ＤＤＳ５４とマルチプレクサ５５を設けるようにしたので、安定した加速サイクルのフラットトップを得ることができる。
【００７１】
実施の形態４．
本実施の形態４では、上記実施の形態３で示した高周波制御装置のマルチ出力ＤＤＳ部分の構成を変更したものである。図９は、本発明の実施の形態４にかかる高周波制御装置の構成を説明するためのブロック図である。
【００７２】
図に示すように、本実施の形態４にかかる高周波制御装置３では、マルチ出力ＤＤＳ５４の代わりに、ＤＤＳ５０とＤ／Ａコンバータ５７と電圧可変移相器５６と加算器５８とを備えるようにしたものである。別の言い方をすれば、上記各実施の形態１ないし３に記載のデジタル制御回路３０の機能を有するとともに、ＤＤＳ５０とＤ／Ａコンバータ５７と電圧可変移相器５６と加算器５８およびマルチプレクサ５５を具備した構成である。
【００７３】
次に実施の形態４にかかる高周波制御装置の動作について説明する。
実施の形態３ではマルチ出力ＤＤＳ５４の信号出力を異なる周波数設定にしてマルチプレクサ５５で出力周波数の切り替えを行った例を示した。本実施の形態４では、メモリモジュール３１から波形演算制御部３３が読み出した高周波信号制御データによって与えられるＤＤＳ５０の出力周波数設定を基準として高周波信号を電圧可変移相器５６に入力する。そして、波形演算制御部３３によるＤ／Ａコンバータ５７の制御から得られる位相制御電圧レベルに比例して電圧可変移相器５６の出力移相量を設定することにより、位相を任意に調整する機能を有するようにした。また、ＤＤＳ５０の出力信号と電圧可変移相器５６の出力信号を加算器５８で合成し、両者の中間となる位相を作る機能も有する。
【００７４】
電圧可変移相器５６の出力移相量は、事前評価調整に基づき最適な値を決めておき、加速器６の運転開始前に計算機１９より予め波形演算制御部３３に設定される。
【００７５】
そして、波形演算制御部３３は補正判定信号Ｓ６の４ビットの値から周波数に関する２ビットの状態により補正トリガ信号Ｓ９のタイミングでマルチプレクサ５５の入力ポートを周波数選択信号Ｓ８により切り替える。補正判定信号Ｓ６の２ビットが００または１０ならばＤＤＳ５０の出力信号を選択、０１ならば電圧可変移相器５６の出力信号を選択、１１ならばＤＤＳ５０と電圧可変移相器５６の合成出力信号を選択する。
【００７６】
こうして、加速器６の運転中の粒子線加速状態を監視し、しきい値を超えると直ちに補正処理部３４から波形演算制御部３３へ情報を与える。これにより、位相制御フィードバック信号Ｓ２０と位置制御フィードバック信号Ｓ２１の各フィードバックモニタ信号がフィードバック基準値に近づくよう引き込むトラッキング制御を行い、安定した加速サイクルのフラットトップを得ることができる。また高めの周波数と低めの周波数の設定値は計算機１９で調整設定できることで、加速器６の試験運転時に柔軟に対応ができ最適なパラメータ調整ができる効果がある。
【００７７】
なお、周波数選択信号Ｓ８を制御してマルチプレクサ５５の入力ポートを任意に切り替えることで乗算器５２に入力する信号を位相変調することも可能である。
【００７８】
尚、上記実施の形態１ないし実施の形態４に係る高周波制御装置３は、その構成の少なくとも一部をプログラマブル半導体デバイスの一種であるＦＰＧＡ（Field-Programmable Gate Array）を用いて構成することもできる。また、上記各実施の形態にかかる高周波制御装置９は、適宜その構成を増減したり組み合わせを変更することができる。
【００７９】
以上のように、本実施の形態４にかかる高周波制御装置３によれば、制御部（５０、５２、５３、５５〜５８部分）は、フラットトップ期間Ｔ５において、フィードバック信号が予め定められている範囲を超えると、当該フィードバック信号を範囲内に入るように補正する補正部として機能するＤＤＳ５０、電圧可変位相器５６、マルチプレクサ５５、Ｄ／Ａコンバータ５７等を設けるようにしたので、安定した加速サイクルのフラットトップを得ることができる。
【００８０】
実施の形態５．
本実施の形態５にかかる粒子線治療装置は、加速器の加速空洞の制御に、上記各実施の形態にかかる高周波制御装置を用いたものである。図１０は本実施の形態５にかかる粒子線治療装置１全体の構成を示す図である。なお、実際の粒子線治療装置では、患者に粒子線を照射する照射装置を複数備えていることが一般的であるが、簡略化のため、ひとつの照射装置のみを表示している。
【００８１】
粒子線治療装置１は、図１０に示すように、粒子線の供給源として、シンクロトロンである加速器６と、患者Ｋに対して粒子線を最適な形態（照射野）で照射するための照射装置８と、加速器６と照射装置８とをつなぎ、加速器６から粒子線を照射装置８に輸送する輸送経路７と、加速空洞６２の制御を行う高周波制御装置３を含む上述した各系統を連携して制御する図示しない制御装置とを備えている。そして、本実施の形態５にかかる粒子線治療装置１に特徴的な構成は、加速器６に設けた加速空洞６２の制御に、上記各実施の形態１〜４で説明した高周波制御装置３を備えるようにしたことである。
【００８２】
加速器６は、粒子線が周回する軌道経路となる真空ダクト、前段加速器６６から供給された粒子線を軌道経路に導くための入射装置６５、粒子線が周回軌道に沿って周回するように、粒子線の軌道を偏向させるための偏向電磁石６１、粒子線を加速させる加速空洞６２、輸送経路７に出射するための出射装置６７等を備えている。なお、前段加速器６６は、図では簡略化のためにひとつの機器のように記載しているが、実際には、陽子、炭素（重粒子）等の荷電粒子（イオン）を発生させるイオン源（イオンビーム発生装置）と、発生された荷電粒子を初期加速する線形加速器系とを備えている。そして、前段加速器６６から周回軌道内に導かれた荷電粒子は、高周波数の電界で加速され、磁石で曲げられながら、光速の約６０％まで加速された粒子線となる。
【００８３】
加速器６により加速された粒子線は、ＨＥＢＴ（高エネルギービーム輸送：High Energy Beam Transport）系と称される輸送経路７へと出射される。輸送経路７は、粒子線の輸送経路となる真空ダクトと、粒子線の軌道を切替える切替装置である切替電磁石と、粒子線を所定角度に偏向する偏向電磁石とを備えている。そして加速器６により十分にエネルギーが与えられ、出射装置６７から取り出された粒子線は、切替電磁石で必要に応じて軌道を変えることで、複数の照射装置のうちの指定された照射装置へと導くことができる。
【００８４】
照射装置８は、輸送経路７から供給された粒子線を照射対象である患者Ｋの患部の大きさや深さに応じた照射野に成形して患部へ照射するものである。照射野を成形する方法は複数あるが、例えば、粒子線を走査させて照射野を形成するスキャニング照射法では、とくに供給された粒子線の軌道精度が形成する照射野の精度に大きく影響する。したがって、上述した各実施の形態１にかかる高周波制御装置３を用いることにより、パラメータ設定が容易で、突変等の影響が抑制され、正確な軌道とエネルギーの粒子線が供給されるので、治療計画等で設定した通りに照射野を形成することができ、周辺組織への影響を最低限にして効果的な治療を行うことができる。なお、照射装置等の具体例として、例えば、文献（三菱電機技報、Ｖｏｌ．６９、Ｎｏ２、（１９９５）、ｐ３４、上田和宏他「ＨＩＭＡＣ用治療・照射システム」）等に挙げられた構成のものを用いてもよい。
【００８５】
以上のように、本実施の形態５にかかる粒子線治療装置１によれば、粒子線を加速するための加速空洞６２が設けられた加速器６と、加速空洞６２を制御するための各実施の形態にかかる高周波制御装置３と、加速器６によりエネルギーを調整されて供給された粒子線を、照射対象に応じた照射野に成形して照射対象へ照射する照射装置８と、を備えるように構成したので、運転パラメータを柔軟に変更して容易に最適化できるとともに、突変等の影響を抑えて粒子線を自在に加速できるので、正確な照射が可能な粒子線治療装置１を得ることができる。
【符号の説明】
【００８６】
１：粒子線治療装置、
３：高周波制御装置、
６：加速器（シンクロトロン）、
７：輸送系、
８：照射装置、
３０：デジタル制御回路（制御部）、
３１：メモリモジュール（記憶部）、
３２：タイミング制御部、
３３：波形演算制御部（制御データ生成部）、
３４：補正処理部（補正部）、
５０：ダイレクトデジタルシンセサイザ（制御部）、
５１：振幅変調波生成器（制御部）、
５４：マルチ出力ダイレクトデジタルシンセサイザ（制御部（補正部））、
６２：加速空洞。
【図１】
2014002986.tif 000003
【図２】
2014002986.tif 000004
【図３】
2014002986.tif 000005
【図４】
2014002986.tif 000006
【図５】
2014002986.tif 000007
【図６】
2014002986.tif 000008
【図７】
2014002986.tif 000009
【図８】
2014002986.tif 000010
【図９】
2014002986.tif 000011
【図１０】
2014002986.tif 000012