(57)【要約】振動呼吸療法デバイス（100）は、呼気がデバイスを通過する間に開口（10）を開閉するロッカーアーム（12）上の弁要素（11）を有することにより、音を発生させる。マイクロフォン（21）を含む携帯電話（20）は、発生した音を拾い上げると共に、音信号を音エネルギー信号に変換してその周波数を測定するためのプロセッサ（22）を含む。プロセッサ（22）は、検出された周波数に固定倍率を乗じ、その積に固定定数を加算することによってデバイス（10）内の圧力測定値を計算する。圧力データは、患者によるデバイスの使用を監視するために使用される。【選択図】図３
2021525617.tif 000002
【特許請求の範囲】
【請求項１】
通過する呼吸流に対して交互抵抗を与えるよう構成されたデバイス（100）と、該デバイスが生成した音を感知するよう構成されたセンサユニット（20）とを備える呼吸療法装置であって、
前記センサユニット（20）が、前記デバイス（100）内で生成された圧力測定値を音周波数から計算するよう構成されたプロセッサ（22）を含むことを特徴とする装置。
【請求項２】
請求項１に記載の装置であって、前記プロセッサ（22）が、前記検出された周波数に固定倍率を乗じ、該積に固定定数を加算することによって前記圧力測定値を計算するよう構成されていることを特徴とする装置。
【請求項３】
請求項２に記載の装置であって、前記倍率が、約1.119であることを特徴とする装置。
【請求項４】
請求項２又は３に記載の装置であって、前記定数が、約-4.659であることを特徴とする装置。
【請求項５】
請求項１〜４の何れか一項に記載の装置であって、前記プロセッサ（22）が、音信号を音エネルギー信号に変換すると共に、該音エネルギー信号を前記圧力測定値の計算に使用することを特徴とする装置。
【請求項６】
請求項１〜５の何れか一項に記載の装置であって、前記センサユニット（20）が、可聴音に応答するマイクロフォン（21）を含むことを特徴とする装置。
【請求項７】
請求項１〜６の何れか一項に記載の装置であって、前記デバイス（100）が、呼気が前記デバイスを通過する間に開口（10）を開閉するロッカーアーム（12）上の弁要素（11）を有することを特徴とする装置。
【請求項８】
請求項１〜７の何れか一項に記載の装置であって、前記センサユニット（20）が、前記デバイス（100）上に取り付けられているのではなく、前記デバイスから分離されていることを特徴とする装置。
【請求項９】
請求項８に記載の装置であって、前記センサユニットが、マイクロフォン（21）及び画面（24）を含む携帯電話（20）によって提供され、前記プロセッサ（22）も、前記電話によって提供されていることを特徴とする装置。
【請求項１０】
呼吸療法デバイス（100）から発せられる音を感知するよう構成されたセンサユニット（20）であって、
前記センサユニット（20）が、前記音周波数から前記デバイス（100）内の圧力測定値を計算するよう構成されたプロセッサ（22）を含むことを特徴とするセンサユニット。
【請求項１１】
デバイスを通過する呼吸に対して振動抵抗を与えるタイプの前記呼吸療法デバイス（100）の使用を評価する方法であって、
前記デバイス（100）の使用によって生じると共に、空気を介して伝達される圧力波を監視するステップと、前記圧力波の周波数測定値を計算するステップと、前記デバイス（100）内の圧力測定値を前記周波数測定値から計算するステップとを含むことを特徴とする方法。
【請求項１２】
請求項１１に記載の方法であって、前記圧力測定値を、前記検出した周波数に固定倍率を乗じ、該積に固定定数を加算することによって計算することを特徴とする方法。
【請求項１３】
請求項１２に記載の方法であって、前記倍率が、約1.119であり、前記定数が、約-4.659であることを特徴とする方法。
【請求項１４】
請求項１１〜１３に記載の方法であって、音信号を音エネルギー信号に変換すると共に、前記音エネルギー信号を前記圧力測定値の計算に使用することを含むことを特徴とする方法。
【請求項１５】
請求項１１〜１４の何れか一項に記載の方法で使用するための装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、通過する呼吸流に対して交互抵抗を与えるよう構成されたデバイスと、そのデバイスが生成した音を感知するよう構成されたセンサユニットとを備えるタイプの呼吸療法装置に関する。
【背景技術】
【０００２】
呼気陽圧（PEP：positive expiratory pressure）装置、即ち、デバイスを通過する呼気に抵抗を与える装置は、様々な呼吸障害、例えば、慢性閉塞性肺疾患、気管支炎、嚢胞性線維症、並びに無気肺を患っている患者の治療を補助するのに広く使用されている。より最近では、流れに対して交互抵抗を与えるこのような装置は、特に有効であることが明らかになっている。このような装置の一例としては、スミス・メディカル社によってアカペラ（スミス・メディカル社の登録商標）の商標名の下で販売されており、特許文献１（米国特許第6581598号明細書）、特許文献２（米国特許第6776159号明細書）、特許文献３（米国特許第7059324号明細書）、並びに特許文献４（米国特許第7699054号明細書）に記載されている。他の振動呼吸療法装置としては、例えば、セイヤー（Thayer）社製の「Quake」、モナハン「Monaghan」社製の「AeroPEP」、スミス・メディカル社製の「TheraPEP」、並びにパーカッショネア（Percussionaire）社製の「IPV Percussionator」が入手可能である。フィリップス社製の「CoughAssist」のような他の装置も入手可能である。呼吸療法装置は逆に、流れに対して吸気中に交互抵抗を与えることができる。特許文献５（国際公開第14202923号パンフレット）に記載の装置は、療法デバイスの使用によって生成された音を監視すると共に、装置内の振動弁における動作周波数の指標を提供する。
【０００３】
これら装置を効果的に使用するには、所定の間隔で定期的に使用しなければならない。慢性疾患の場合、患者は、症状の緩和を持続的に維持するために残りの人生にわたって装置を毎日使用する必要がある。
【０００４】
これら装置は極めて効果的ではあるが、使用者はしばしば、装置を定期的かつ所定の頻度で使用するのを怠ることがある。装置の使用記録を残すことは極めて困難であり、特に患者が装置を自宅で使用している場合に困難である。臨床医は、患者の状態の悪化が指示通りに装置の使用を怠ったからなのか、又は他の要因によるものであるかを把握していないことが多い。使用を監視する目的で、又は患者にフィードバックを提供して正しい使用を促す目的で、療法デバイス内の圧力を監視することは有用であり得る。圧力を監視するために、療法デバイス内に圧力センサを取り付けることが以前に提案されている。しかしながら、この提案は、通常は使い捨てとして想定されるデバイスのコスト増につながると共に、デバイス内に電力が必要になるため、電気的安全性の懸念や、電池が使用されている場合には廃棄上の問題につながる可能性がある。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】米国特許第6581598号明細書
【特許文献２】米国特許第6776159号明細書
【特許文献３】米国特許第7059324号明細書
【特許文献４】米国特許第7699054号明細書
【特許文献５】国際公開第14202923号パンフレット
【発明の概要】
【０００６】
本発明の目的は、代替的な呼吸療法装置及び呼吸療法装置の使用を評価する方法を提供することである。
【０００７】
本発明の一態様によれば、上述したタイプの呼吸療法装置が提供され、その呼吸療法装置は、センサユニットが、音の周波数によって装置内で生成された圧力測定値を計算するよう構成されたプロセッサを含むことを特徴とする。
【０００８】
プロセッサは、好適には、検出された周波数に固定倍率を乗じ、その積に固定定数を加算することによって圧力測定値を計算するよう構成されている。この場合、倍率は、約1.119とすることが可能であり、定数は、約-4.659とすることが可能である。プロセッサは、好適には、音信号を音エネルギー信号に変換すると共に、その音エネルギー信号を圧力測定値の計算に使用する。センサユニットは、好適には、可聴音に応答するマイクロフォンを含む。デバイスは、呼気がデバイスを通過する間に開口を開閉するロッカーアーム上の弁要素を有することができる。センサユニットは、好適には、デバイス上に取り付けられているのではなく、デバイスから分離されている。センサユニットは、マイクロフォン及び画面を含む携帯電話によって提供され、プロセッサも、その電話によって提供されている。
【０００９】
本発明の他の一態様によれば、呼吸療法デバイスから発せられる音を感知するよう構成されたセンサユニットが提供され、そのセンサユニットは、音周波数からデバイス内の圧力測定値を計算するよう構成されたプロセッサを含むことを特徴とする。
【００１０】
本発明の更なる態様によれば、デバイスを通過する呼吸に対して振動抵抗を与えるタイプの呼吸療法デバイスの使用を評価する方法が提供され、その方法は、デバイスの使用によって生じると共に、空気を介して伝達される圧力波を監視するステップと、圧力波の周波数測定値を計算するステップと、デバイス内の圧力測定値を周波数測定値から計算するステップとを含むことを特徴とする。
【００１１】
圧力測定値は、検出した周波数に固定倍率を乗じ、その積に固定定数を加算することによって計算することができる。倍率は、約1.119とすることが可能であり、定数は、約-4.659とすることが可能である。本発明の方法は、好適には、音信号を音エネルギー信号に変換すると共に、音エネルギー信号を圧力測定値の計算に使用することを含む。
【００１２】
本発明の第４態様によれば、上述した更なる態様に係る方法において使用するための装置が提供される。
【００１３】
以下、振動PEPデバイスを含む装置を添付図面に基づいて例示的に説明する。
【図面の簡単な説明】
【００１４】
【図１】装置を示す分解図である。
【図２】一回の呼吸で生成された音の出力を示す説明図である。
【図３】使用中の装置を示す説明図である。
【図４】時間に対する音の振幅を示すグラフである。
【図５】図４におけるトレースを、拡大した時間にわたって示すグラフである。
【図６】図５におけるグラフを、音エネルギーに変換された状態で示すグラフである。
【図７】図６におけるグラフを、フィルタリング後の状態で示すグラフである。
【図８】周波数と圧力との関係を示すグラフである。
【図９】療法セッション中の携帯電話の画面を示す説明図である。
【図１０】一週間にわたる療法歴の表示が異なる、使用後の携帯電話の代替的な画面を示す説明図である。
【図１１】一週間にわたる療法歴の表示が異なる、使用後の携帯電話の代替的な
【発明を実施するための形態】
【００１５】
まず図１につき説明すると、デバイス100は、外側ハウジング２内に収容したロッカーアセンブリ１を有し、ハウジング２は、ほぼ半円筒形状の上側部分３及び下側部分４を有する。デバイスは、円形断面を有する調整可能ダイヤル５によって組み立てが完了する。ロッカーアセンブリ１は、一端に呼吸入口７を、他端に一方向弁（図示せず）が設けられた吸気入口８を含む空気流チューブ６を有し、その一方向弁により、空気流チューブ６内に空気が流入するが、吸気入口からの空気の流出は阻止される。空気流チューブ６は、非線形プロファイルの出口開口10を有し、その出口開口10は、ロッカーアーム12上に取り付けられた円錐状弁要素11によって開閉され、ロッカーアーム12はその長さの中間位置で横方向軸線周りに旋回する。空気流チューブ６及びハウジング２により、ロッカーアーム12が取り付けられる構造がもたらされる。ロッカーアーム12の呼吸入口７から遠い方の端部において、ロッカーアーム12は、調整可能な支持フレーム14上に取り付けられた永久磁石（図示せず）によって生ずる磁場と相互作用する鉄ピン13を保持している。この磁石構成は、患者がデバイスを通して呼吸していないときには、ロッカーアーム12の遠位端部が押し下げられ、従って弁要素11も出口開口10に封止係合するよう下方に押し下げられるものである。支持フレーム14の一端におけるカムフォロワ突出部15は、ダイヤル５におけるカム溝孔16内に位置し、ダイヤルを回転することにより、支持フレーム14が磁石と一緒に上下動可能であり、これにより鉄ピン３と相互作用する磁場強度が変化する。ダイヤル５により、使用者に対する療法の有用性を最大化するよう、デバイスを通過する空気流の動作周波数及び抵抗が調整可能である。
【００１６】
患者が呼吸入口７から空気を吸い込むと、空気が吸気入口８を通過して空気流チューブ６に沿って呼吸入口まで引き込まれる。患者が息を吐き出すと、吸気入口８における一方向弁が閉鎖し、この経路に沿っていかなる空気の流出も阻止される。その代わりに、呼気圧がロッカーアーム12における弁要素11の下側に加わり、弁要素11を磁石吸着力に抗して開口10から持ち上げ、これにより空気が大気に流出することができる。開口10は、ロッカーアーム12の遠位端部が持ち上がるにつれて有効排出面積が増加する非線形輪郭を有し、従ってアームが復帰下降して開口を閉じることが可能である。使用者が十分な呼気圧を加え続ける限り、開口10の開閉に合わせてロッカーアーム12が上昇及び下降を繰り返し、これによりデバイスを通過する呼気流に対して振動的、交互的、又は揺動的な抵抗をもたらす。ロッカーアーム12のこの振動運動は、療法デバイス２の外部に聞こえるガタガタ音を発生させる。デバイスの構造及び動作に関する更なる情報は特許文献１に記載されているが、その内容は本発明の理解にとって必須ではない。
【００１７】
ここまで説明したように、装置は従来のものである。
【００１８】
本発明の装置は、上述したデバイス100と、空気を介して伝達されると共に、デバイスの使用で生じる圧力波に応答するセンサ手段20とを備える。この圧力波は、典型的には、音響周波数内にあり、従来のマイクロフォンで感知できる可聴音を生成する。
【００１９】
図２は、所望の療法効果を得るために、患者がデバイス100を通して息を吐く様子を示す。ロッカーアーム12の振動運動は、周囲の空気を介して、圧力波として伝達される可聴音を生成する。これは、図２の右側に示す一回の呼気におけるトレース「Ｔ」によって表されている。図２に明らかなように、感知された音は、ロッカーアーム12における振動の周波数に依存する周波数で、正と負のピークを急速かつ交互に繰り返しており、呼吸の開始に際して最大値まで増加すると共に、呼吸の終了に際して徐々に減少する平均振幅を有する。
【００２０】
図３は、デバイス100と、携帯電話20の形態としたセンサ手段とを示す。この場合、携帯電話20には、内蔵されたマイクロフォン21の出力を処理するよう構成されたアプリがプログラムされている。これにより、電話20は、デバイス２とは別個ではあるが、使用時にデバイス近傍に配置されるスタンドアロン型の音響センサとして機能する。電話20内のマイクロフォン21は、メモリを含むプロセッサ22に接続され、表示画面24に出力するようになっている。マイクロフォン21は、可聴範囲の音に応答する。プロセッサ22にプログラムされたアプリは、マイクロフォン21からの信号の周波数を判定するよう構成された周波数検出アルゴリズムを含む。
【００２１】
作動中の療法デバイスにおけるフラッター周波数を測定するよう構成された周波数検出アルゴリズムは、時間領域における自己相関に基づいており、音の振幅ではなく音のエネルギーを具体化する、フィルタリングされた信号で動作する。音エネルギーは音振幅及び音周波数に依存し、波構造をより明確に表すことが分かっている。
【００２２】
図４は、療法デバイスから検出された音声信号を示し、図５は、その信号をより短時間にわたって拡大した状態で示す。図示するように、周波数又はピッチパターンは識別可能ではあるが、その上に重ね合わされた高周波信号によって混同されている。図６は、拡大された信号を単純な振幅から音エネルギーに変換する効果を示しており、波構造が遥かに明確である。図７に示すように、これをフィルタリングし、容易に識別可能なピーク及びトラフを有する平滑化された信号を得ることができる。
【００２３】
アルゴリズムでは、精度及び計算速度に等しく重点が置かれている。例えば、C++で実施した場合、20秒の音声フィードが僅か50ミリ秒で処理可能であり、リアルタイム処理に関して400倍のマージンが得られる。アルゴリズムが携帯電話などで実行される場合でも、このマージンは十分過ぎるほどである。アルゴリズムは、音声データのインクリメンタル・チャンク、即ち通常は0.1又は0.2ミリ秒ごとに最新の音声キャプチャの新たなチャンクで動作する。アルゴリズムの動作範囲は、10〜40 Hzに設定されており、換言すれば、この範囲の振動周波数を検出することを目的としている。
【００２４】
周波数検出アルゴリズムに関する更なる詳細については、以下に記載する。この場合、ピッチという用語は、付加的な高周波ノイズを除去した後の療法デバイス内における弁の振動周波数を意味するものと理解されたい。
【００２５】
・  ピッチの目標範囲が次の通り規定される：[flower fupper]=[10 Hz 40 Hz]
・  FS=48 kHz又は44.1 kHzの音声信号は、（980 Hzまで）45分の１にダウンサンプリングされる。ダウンサンプリング率は、周波数範囲の中央で所望のピッチ精度を2.5％に設定することで得られる
・  生の音声信号から、平滑化されたバージョンが減算される
・  平滑化されたバージョンは、176の音声サンプルポイントを有する最小フォルマント周波数（これは療法デバイス内で生成される典型的な固有周波数であり、250 Hz以上から開始）の長さに等しい近傍値を平均化することによって計算される
・  減算後、二乗を計算する
・  連続する45のサンプルのセットを平均化することにより、信号がダウンサンプリングされる（オーバーラップはなし）
・  ダウンサンプリングされたこの信号は、音声エネルギーと称される
・  音声エネルギーの平滑化されたバージョンは、sigma=0.2*45*fupper/Fsのガウシアンフィルタによる畳み込みによって計算される。これは、ガウシアンのFWHMが、求められる最高ピッチ周波数の周期の0.4倍であることを意味する。これ以上大きくした場合、高域（40 Hz）におけるピッチのシグネチャがぼやけ始める。
・  音声のチャンクが利用可能になると、以前の移動平均を追跡するパイプラインで処理される
・  平滑化された信号において、自己相関は効率的な方法で計算される。即ち、
・  ピッチが検出されると、次の自己相関は、（最後に検出されたピッチ）かつ前のポイントよりも１波長先の位置で計算される
・  自己相関の計算は、最も可能性の高い位置に集中するよう、ピッチの目標範囲に維持される
・  自己相関の時間分解能は、最初の計算では粗く、その後、最大値付近でより細かい分解能が得られる
・  相関が0.6を超える場合、自己相関関数のargmaxを選択する。即ち、
・  ピッチの移動平均値を、argmax値に対応するピッチで更新する
・  次の反復時に検索ウィンドウを絞り込むために、移動平均値の導関数を更新する
・  ピッチを検出するためのラン・レングスを更新する。即ち、
・  ラン・レングスを０に初期化する
・  ピッチが明確に識別される毎に、ラン・レングスが１だけ増加する
・  ただし、ラン・レングスが最大値に達した場合は除く。その最大値は、端数のない減衰率（１-１/decayRate）として与えられる。この場合、減衰率は0.8に設定され、最大ラン率は５である
・  検出された自己相関値（>0.6）から次の値への飛躍において、ピッチに関するランニング平均値が期待値を超える場合、ラン率が同等の数だけ減少される
・  混合値を次の通り計算する：mix=min（[decayRate１-１/runlength]）
・  周期（１/pitch）、振幅、導関数に関して移動平均値を更新する
・  MovingAve=movingAve*mix+newValue*（１-mix）
・  ピッチ又は周波数を次の通り計算する：Pitch=１/（movingAvePeriod）
【００２６】
ピッチ又は周波数検出アルゴリズム及びピッチ/圧力モデルは、ベクトルや文字列などのコアC++ライブラリだけを使用してC++で実施することができる。アルゴリズムが音声を0.2秒のインクリメンタル・チャンクで処理する場合、約50ミリ秒で20秒の音声サンプルが実行される。これにより、400倍のパフォーマンスヘッドルームがもたらされ、携帯電話又は同様の機器のコードに関して十分なヘッドルームが得られる。
【００２７】
広範な実験データから、流れに対して振動抵抗を与えるロッカーアーム12及び弁要素11の振動周波数は、空気流チューブ６内で測定された圧力と線形的かつ密接に関連していることが判明している。図８は、このことを示す。線形モデルを使用した場合、9993のデータポイントに対して0.886のr2が生成され、極めて良好にフィットする。抵抗ダイヤル５の設定にかかわらず、同様の結果が適用されることが判明している。
【００２８】
検出された周波数と圧力との関係は、次の通り表すことができる：P=A+B×f
ここに、Pは、cm単位での水柱圧力、fは周波数、Aは定数、Bは定数倍率である。より具体的には、上述した呼吸療法デバイス２に関しては、P= -4.659+1.119×fであることが判明している。
【００２９】
定数及び倍率Ａ及びＢは、振動抵抗療法デバイスにおいて異なる場合があることを理解されるだろうが、一般的には同じ形態の関係が存在すると想定される。
【００３０】
このことの重要性は、デバイス内に専用の圧力センサを必要とすることなく、検出された動作周波数からのみ、療法デバイス内の圧力に極めて近い指標を得ることができることを意味する。
【００３１】
このように、療法デバイスに修正を施すことなく、圧力の指標を正確に導き出すことが可能である。圧力センサを療法デバイス内に取り付ける必要がないため、コストを最小限に抑えることができる。また、療法デバイスへの電力供給や療法デバイス内における電力供給が不要である。これにより、デバイス内に電気的安全性のための絶縁物を提供するために特別な設備の必要性がなくなり、デバイスに又はデバイスから延在するケーブルの必要性がなくなる。更に、デバイス内の電池は、安全な保管時間を超えないように特別な対策を必要とする輸送安全上の問題を生じさせると共に、廃棄上の様々な問題を生じさせるが、デバイス内におけるその電池の必要性がなくなる。
【００３２】
使用に際して、電話又は他のセンサユニット20は、デバイス100とは別個かつ遠隔に配置されるが、マイクロフォン21の可聴範囲内にあるようにするためにデバイスに対して十分に近傍に位置している。この場合、センサユニットは、デバイスと直接接触する必要はない。電話20内の療法監視アプリは、探された後に起動される。使用者が療法セッションを開始すると、デバイス100が音波を発し始め、その音波が電話20内のマイクロフォン21によって拾われる。電話20内のプロセッサ22は、周波数検出アルゴリズムを実行して周波数値を識別して次式、即ち-4.659+1.119×fを適用することにより、cm単位での水柱圧力測定値が導き出される。この圧力測定値は、特定の使用者にとって適切な目標値及び限界値を表す、プロセッサ22内に保存されている事前プログラムされた圧力値と比較される。その後、圧力は、図９に示すように、療法セッション中に電話20の画面24上に、数値表示及びグラフ表示と組み合わせて表示される。圧力は、典型的には、小さ過ぎる、大き過ぎる、又は適性であることを示す異なる色のスケールに沿って表示される。療法セッションの終了後、電話画面24は、療法歴を表示するために切り替わる。図10及び図11は、２つの異なる代替的な履歴表示を示す。この場合、典型的には、平均呼気圧、呼気時間、一定期間（例えば毎週）にわたる使用回数、並びに療法セッション中の呼気回数が表されている。
【００３３】
上述したように、療法デバイス100におけるダイヤル５の設定は、デバイスを通して流れる周波数及び抵抗に影響を及ぼす。これは、有益な効果を最大限に得るために使用者によって設定される。電話20内のアプリは、使用者又は臨床医によるダイヤル設定の手動入力を受信するよう構成することができる。
【００３４】
上述したセンサ手段は、携帯電話20内のプログラム又はアプリで提供されるが、代替的に、適切にプログラムされたプロセッサを有する専用のセンサユニットで提供されてもよい。そのようなセンサユニットは、例えば、デバイス上でセンサを支持するクリップ又はストラップによって療法デバイスに取り付けることができる。センサ手段は、汎用コンピュータ内のプログラムアプリケーションによって提供されるような多機能ユニットにより、コンピュータに内蔵されたマイクロフォンを使用して、又は別箇のプラグインマイクロフォンを使用して提供することができる。プログラムアプリケーションは、特徴的な周波数の音を感知することなく、所定時間経過後に自動停止するよう構成することができる。
【００３５】
センサユニットで得られた情報を多くの異なる方法で表現し、従って使用者及び臨床医に最適に提供できることが理解されよう。
【００３６】
本発明の装置は、音信号を生成する任意かつ従来の呼吸療法装置と共に使用することができる。本発明の療法装置は、噴霧療法やエアゾール薬の投与など、他の治療法と組み合わせることができる。
【００３７】
本発明は、既存かつ従来の療法デバイスを使用している人に、その使用に関する有用なデータの提供を可能にする。これにより、使用者は、指示された治療プログラムをどの程度遵守しているかをより認識可能であり、従って最大の効果を得られるようデバイスの使用を適宜変更することができる。臨床医も、患者の順守状況をチェックすることが可能であり、従って患者の症状悪化が順守欠如によるものか否か、又は他の治療が必要か否かを特定することができる。
【図１】
2021525617.tif 000003
【図２】
2021525617.tif 000004
【図３】
2021525617.tif 000005
【図４】
2021525617.tif 000006
【図５】
2021525617.tif 000007
【図６】
2021525617.tif 000008
【図７】
2021525617.tif 000009
【図８】
2021525617.tif 000010
【図９】
2021525617.tif 000011
【図１０】
2021525617.tif 000012
【図１１】
2021525617.tif 000013
【国際調査報告】
2021525617.tif 000014
2021525617.tif 000015
2021525617.tif 000016