2023-136987 生体状態評価方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2023136987

(43)【公開日】2023-09-29

(54)【発明の名称】生体状態評価方法

(51)【国際特許分類】

   A61B 5/0537 20210101AFI20230922BHJP

【ＦＩ】

A61B5/0537 200

【審査請求】未請求

【請求項の数】5

【出願形態】ＯＬ

(21)【出願番号】P 2022042947

(22)【出願日】2022-03-17

(71)【出願人】

【識別番号】000219314

【氏名又は名称】東レエンジニアリング株式会社

(74)【代理人】

【識別番号】110001427

【氏名又は名称】弁理士法人前田特許事務所

(72)【発明者】

【氏名】坂上 友介

【テーマコード（参考）】

4C127

【Ｆターム（参考）】

4C127AA06

4C127EE01

4C127EE03

(57)【要約】

【課題】生体状態を評価する際の測定時間を短縮する。
【解決手段】２つの電極１０、２０を生体表面３０に配置し、２つの電極１０、２０を介して測定装置１と生体とで形成される電気回路Ｃを、第１の状態から第２の状態に変化させ、第１の状態から第２の状態に変化させたときに電気回路Ｃに生じる過渡応答ＲＴの情報を取得し、取得した過渡応答ＲＴの情報に基づいて生体の情報を評価する。
【選択図】図３

【特許請求の範囲】

【請求項1】

電極を有する測定装置を用いた生体状態評価方法であって、
前記電極を生体の表面に配置し、
前記電極を介して前記測定装置と前記生体とで形成される電気回路を、第１の状態から第２の状態に変化させ、
前記第１の状態から前記第２の状態に変化させたときに前記電気回路に生じる過渡応答の情報を取得し、
取得した前記過渡応答の情報に基づいて、前記生体の情報を評価することを特徴とする、生体状態評価方法。

【請求項2】

前記過渡応答の情報とは、該過渡応答の整定時間であることを特徴とする、請求項１に記載の生体状態評価方法。

【請求項3】

前記生体の情報とは、生体の筋疲労度であることを特徴とする、請求項１又は２に記載の生体状態評価方法。

【請求項4】

前記第１の状態から前記第２の状態に変化させることとは、第１の外部抵抗が並列接続された状態から第２の外部抵抗が並列接続された状態へと前記電気回路を切り替えることであることを特徴とする、請求項１から３のいずれか１つに記載の生体状態評価方法。

【請求項5】

前記第１の外部抵抗及び前記第２の外部抵抗のうちのいずれか一方の抵抗値が、無限大に設定されていることを特徴とする、請求項４に記載の生体状態評価方法。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、生体状態評価方法に関する。

【背景技術】

【0002】

生体状態を評価する方法として、生体表面に電極を配置して、その電極を介して形成される電気回路を用いた方法が知られている。

【0003】

例えば特許文献１には、生体状態評価方法の一例として、生体表面に少なくとも３つの電極を配置するとともに、それら各電極とグランド電位との間に、第１の外部抵抗を並列接続したときに各電極に生じる第１の電圧と、それら各電極とグランド電位との間に、第２の外部抵抗を並列接続したときに各電極に生じる第２の電圧とをそれぞれ測定することが開示されている。

【0004】

前記特許文献１によると、第１の電圧及び第２の電圧から３つの比を算出することで、生体内の信号源の位置を検出することができる。

【0005】

一方、特許文献２には、生体状態評価方法の別例として、生体表面に少なくとも２つの電極を配置して、その２つの電極間に第１の外部抵抗を並列接続したときに生じる第１の電圧と、第２の外部抵抗を並列接続したときに生じる第２の電圧とをそれぞれ測定することが開示されている。

【0006】

前記特許文献２によると、第１の電圧及び第２の電圧の測定結果から得られる電圧比に基づいて、生体表面下の筋肉部位における２つの電極間の生体インピーダンスを算出することができる。そして、算出された生体インピーダンスの時間変化に基づいて、筋肉部位における局所的な筋疲労を評価することができる。

【先行技術文献】

【特許文献】

【0007】

【特許文献1】特許第６４９２１０３号

【特許文献2】特開２０２０－１５７０３７号公報

【発明の概要】

【発明が解決しようとする課題】

【0008】

前記特許文献１及び２に記載された方法は、いずれも、第１の外部抵抗を接続した状態と、第２の外部抵抗を接続した状態とのそれぞれで測定された電位を用いるように構成されている。

【0009】

ところが、第２の外部抵抗が接続された状態での測定に際しては、第１の外部抵抗から第２の外部抵抗へと切り替えた場合に生じる過渡応答が収束した後に、電位を測定する必要があった。そのため、従来知られた方法では、外部抵抗の切り替え直後に電位を測定することができず、測定時間の短縮という点で性能向上の余地があった。

【0010】

本開示は、かかる点に鑑みてなされたものであり、その目的とするところは、生体状態を評価する際の測定時間を短縮することにある。

【課題を解決するための手段】

【0011】

本開示は、電極を有する測定装置を用いた生体状態評価方法に係る。この生体状態評価方法は、前記電極を生体の表面に配置し、前記電極を介して前記測定装置と前記生体とで形成される電気回路を、第１の状態から第２の状態に変化させ、前記第１の状態から前記第２の状態に変化させたときに前記電気回路に生じる過渡応答の情報を取得し、取得した前記過渡応答の情報に基づいて、前記生体の状態を評価することを特徴としている。

【0012】

前記構成によれば、過渡応答の収束後に電位を測定するのではなく、その過渡応答の情報それ自体を活用することで生体の情報を評価する。過渡応答の収束を待つ必要がない分、測定時間を短縮することができる。

【0013】

また、前記過渡応答の情報とは、該過渡応答の整定時間であってもよい。

【0014】

前記構成によれば、過渡応答の整定時間を取得することで、その整定時間に対応した時定数を算出することができる。そうして算出した時定数は、いわゆる抵抗値（レジスタンス）ばかりでなく、静電容量（キャパシタンス）の情報も含んでいる。したがって、抵抗値のみを用いた評価方法と比較して、より高精度に生体の状態を評価したり、より多方面からの評価を行ったりすることができるようになる。

【0015】

また、前記生体の情報とは、生体の筋疲労度であってもよい。

【0016】

また、前記第１の状態から前記第２の状態に変化させることとは、第１の外部抵抗が並列接続された状態から第２の外部抵抗が並列接続された状態へと前記電気回路を切り替えることであってもよい。

【0017】

前記構成によれば、電気回路の切り替えを利用して生体インピーダンス等の電気的なパラメータを算出し、その算出結果に基づいて生体の状態を評価することができるようになる。

【0018】

また、前記第１の外部抵抗及び前記第２の外部抵抗のうちのいずれか一方の抵抗値が、無限大に設定されていてもよい。

【0019】

前記構成によれば、第１の状態と第２の状態との間の切り替えを、電気回路が接続された状態と電気回路が遮断された状態との間の切り替えとみなすことができる。これにより、生体インピーダンス等、各種パラメータの計算式がよりシンプルなものとなり、生体の状態をより容易に評価することが可能になる。

【発明の効果】

【0020】

以上説明したように、本開示によれば、生体状態を評価する際の測定時間を短縮することができる。

【図面の簡単な説明】

【0021】

【図1】図１は、一実施形態における生体状態評価方法を模式的に示した図である。

【図2】図２は、図１に示した状態の等価回路図を示す。

【図3】図３は、第１の状態と第２の状態との切り替えに伴う過渡応答を例示するタイムチャートである。

【図4】図４（ａ）及び（ｂ）は、筋肉内の水分量と、生体インピーダンスの変化との関係を模式的に示したグラフである。

【発明を実施するための形態】

【0022】

以下、本開示の実施形態を図面に基づいて詳細に説明する。なお、以下の説明は例示である。また、以下の図面においては、説明の簡潔化のため、実質的に同一の機能を有する構成要素に同一の参照符号を付す。

【0023】

＜過渡応答を用いた筋疲労評価方法＞
図１は、本開示の一実施形態における、生体状態評価方法としての筋疲労評価方法を模式的に示した図である。

【0024】

図１に示す筋疲労評価方法は、図１に示す測定装置１を用いて実施することができる。この測定装置１は、２つの電極１０、２０と、電圧測定手段としての増幅器４０と、第１の外部抵抗Ｒg1と、第２の外部抵抗Ｒg2と、接続手段としてのスイッチＳＷとを有する。

【0025】

図１に示すように、生体の表面（生体表面）３０に、２つの電極１０、２０を所定の間隔を空けて配置する。ここでは、上腕部の表面に２つの電極１０、２０を貼り付けた例を示す。そして、電極１０、２０間に生じる電圧を、増幅器４０で増幅して測定する。また、２つの電極１０、２０の間に、第１の外部抵抗Ｒg1と、第２の外部抵抗Ｒg2とが並列に配置されている。そして、スイッチＳＷによって、電極１０、２０間に、第１の外部抵抗Ｒg1が並列接続された状態（第１の状態）と、第２の外部抵抗Ｒg2が並列接続された状態（第２の状態）とに切り替える。

【0026】

図２は、図１に示した状態の等価回路図を示す。

【0027】

ここで、Ｖbは、２つの電極１０、２０間の生体表面３０下にある筋肉部位（上腕筋）における筋電位である。この筋電位Ｖbは、上腕部３０を運動したとき、すなわち、上腕筋に負荷をかけたときに発生する。筋電位Ｖbを発生する信号源Ｓは、生体内における２つの電極１０、２０の間に介在している。

【0028】

なお、本明細書において、筋肉（部位）とは、上腕二頭筋、広背筋など、生理学上個別に区分された筋組織を指す。

【0029】

また、Ｚbは、生体内における２つの電極１０、２０の間の生体インピーダンスを示す。この生体インピーダンスＺbは、２つの電極１０、２０の間のレジスタンスＲbとキャパシタンスＣbとが並列接続されたＲＣ回路におけるインピーダンスとみなすことができる。また、Ｒinは、増幅器４０の入力抵抗を示す。

【0030】

図２に示すように、２つの電極１０、２０を介して、測定装置１と生体とで電気回路Ｃが形成される。具体的に、この電気回路Ｃは、生体インピーダンスＺbと、信号源Ｓと、２つの電極１０、２０と、第１の外部抵抗Ｒg1及び第２の外部抵抗Ｒg2と、スイッチＳＷと、入力抵抗Ｒinとによって形成されている。電気回路Ｃにおいて２つの電極１０、２０間に生じた電圧（特に、生体外で生じる電位差）は、増幅器４０で増幅されて、出力電圧Ｖoutとして計測される。

【0031】

本実施形態において、電気回路Ｃを第１の状態から第２の状態に変化させることとは、第１の外部抵抗Ｒg1が並列接続された状態から第２の外部抵抗Ｒg2が並列接続された状態へと電気回路Ｃを切り替えることに相当する。この切り替えは、前述のようにスイッチＳＷによって行われるようになっている。

【0032】

なお、スイッチＳＷによる切り替えは、所定周期で繰り返すように構成してもよい。そのように構成した場合、第１の外部抵抗Ｒg1及び第２の外部抵抗Ｒg2の間の切り替えは、第１の外部抵抗Ｒg1から第２の外部抵抗Ｒg2への切り替えと、第２の外部抵抗Ｒg2から第１の外部抵抗Ｒg1への切り替えとからなる工程を、前記所定周期で繰り返し行うことになる。

【0033】

また、第１の外部抵抗Ｒg1及び第２の外部抵抗Ｒg2のうちのいずれか一方は、その抵抗値が無限大に設定されていてもよい。例えば、第１の外部抵抗Ｒg1の抵抗値を有限の値とし、第２の外部抵抗Ｒg2の抵抗値を無限大とした場合、スイッチＳＷによる切り替えは、第１の外部抵抗Ｒg1によって電気回路Ｃが接続された状態と、第２の外部抵抗Ｒg2によって電気回路Ｃが遮断された状態との間の切り替えとみなすことができる。

【0034】

図３は、第１の状態と第２の状態との切り替えに伴う過渡応答を例示するタイムチャートである。

【0035】

従来知られた方法では、第１の外部抵抗Ｒg1を接続した状態（第１の状態）と、第２の外部抵抗Ｒg2を接続した状態（第２の状態）とのそれぞれで出力電圧Ｖoutを測定するとともに、各測定値の電圧比に基づいて、生体インピーダンスＺb、ひいては筋疲労度等の生体の状態を評価するように構成されていた。

【0036】

ここで、第２の状態での測定に際しては、第１の状態から第２の状態へと切り替える必要がある。しかし、そうした切り替えを行った直後、図３に示すような過渡応答ＲＴが生じることになる。

【0037】

したがって、従来知られた方法では、過渡応答ＲＴが収束した後に出力電圧Ｖoutを測定する必要があった。そのため、外部抵抗の切り替え直後に出力電圧Ｖoutを測定することができず、測定時間の短縮という点で性能向上の余地があった。

【0038】

こうした問題に対し、本実施形態に係る方法では、収束後の電圧ではなく、過渡応答を直に測定することで、生体インピーダンスＺb、ひいては生体の状態を評価するように構成されている。

【0039】

具体的に、電気回路Ｃにおいて、第２の外部抵抗Ｒg2が並列接続された状態（第２の状態）から第１の外部抵抗Ｒg1が並列接続された状態（第１の状態）へと変化させたときに２つの電極１０、２０間に生じる第１の電圧Ｖ_１は、下式（１）で与えられる。

【0040】

【数1】

【0041】

上式（１）において、ｔ_１は、第２の状態から第１の状態へ変化させてからの経過時間（スイッチＳＷによる切り替えからの経過時間）であり、τ_１は、第１の状態における時定数である。以下、τ_１を「第１の時定数」という。第１の時定数τ_１は、下式（２）で与えられる。

【0042】

【数2】

【0043】

第１の時定数τ_１は、第２の状態から第１の状態へ変化させたときに生じる過渡応答の整定時間に比例する。この整定時間は、過渡応答がピーク値に達してから、定常値の±ａ％（ａは実数）以内に収束するまでの所要時間を示す。この整定時間に基づいて、第１の時定数τ_１を算出することができる。

【0044】

また、電気回路Ｃにおいて第１の状態から第２の状態へと変化させたきに２つの電極１０、２０間に生じる第２の電圧Ｖ_２は、下式（３）で与えられる。

【0045】

【数3】

【0046】

上式（３）において、ｔ_２は、第１の状態から第２の状態へ変化させてからの経過時間（スイッチＳＷによる切り替えからの経過時間）であり、τ_２は、第２の状態における時定数である。以下、τ_２を「第２の時定数」という。第２の時定数τ_２は、下式（４）で与えられる。

【0047】

【数4】

【0048】

第２の時定数τ_１は、第１の状態から第２の状態へ変化させたときに生じる過渡応答の整定時間、つまり過渡応答がピーク値に達してから、定常値の±ｂ％（ｂは実数）以内に収束するまでの所要時間に比例する。したがって、この整定時間に基づいて、第２の時定数τ_２を算出することができる。

【0049】

式（２）と式（４）を組み合わせることで、生体インピーダンスＺbを構成するレジスタンスＲbとキャパシタンスＣbは、以下の式（５）及び（６）より求めることができる。

【0050】

【数5】

【0051】

【数6】

【0052】

式（５）及び式（６）に示すように、生体表面３０下にある筋肉部位（上腕筋）における２つの電極１０、２０間のレジスタンスＲb及びキャパシタンスＣbは、それぞれ、第１の外部抵抗Ｒg1及び第２の外部抵抗Ｒg2の値と、第１の時定数τ_１及び第２の時定数τ_２の値に基づいて算出することができる。

【0053】

そして、生体インピーダンスＺbを推定したい角周波数をωとすると、生体表面３０下にある筋肉部位（上腕筋）における２つの電極１０、２０間の生体インピーダンスＺbは、以下の式（７）より求めることができる。

【0054】

【数7】

【0055】

式（７）において、ｊは虚数単位を表している。ここで、例えば第２の外部抵抗Ｒg2の値を無限大に設定した場合、上式（７）は、下式（８）のように簡略化することができる。

【0056】

【数8】

【0057】

すなわち、本実施形態に係る生体インピーダンス測定方法では、まず、第１の状態及び第２の状態の間で電気回路Ｃを変化させたときに生じる過渡応答の情報を取得する。特に本実施形態では、過渡応答の情報として、第１の時定数τ_１に対応した過渡応答の整定時間と、第２の時定数τ_２に対応した過渡応答の整定時間を取得する。

【0058】

その後、取得した各整定時間に基づいて、第１の時定数τ_１と第２の時定数τ_２を算出するとともに、算出した第１の時定数τ_１及び第２の時定数τ_２と、第１の外部抵抗Ｒg1及び第２の外部抵抗Ｒg2の値を上式（７）に代入することで、生体表面３０下の筋肉部位における２つの電極１０、２０間の生体インピーダンスＺbを算出するようになっている。このように、本実施形態では、取得した過渡応答の情報に基づいて生体インピーダンスＺbを算出することができる。

【0059】

ところで、筋肉に負荷をかけて筋疲労が生じると、血中の乳酸濃度が増加することがよく知られているが、筋肉内の水分も増加することが分かっている。そのため、疲労した筋肉部位における生体インピーダンスを測定すると、平常時よりも生体インピーダンスが減少していることが予測できる。

【0060】

本願発明者らが鋭意検討を重ねた結果、得られた知見によると、腕の肘曲げ運動のように、筋肉に一定でない負荷をかけたときの生体インピーダンスの変化は、血中乳酸濃度の変化や、筋肉の厚みの変化と、強い相関関係があることが分かった。すなわち、特定の筋肉の疲労は、筋肉内の水分量の変化として捉えることができ、これにより、生体インピーダンスの変化が、生体の筋疲労度を反映する指標になり得ることが分かった。

【0061】

したがって、筋肉に負荷をかけて筋疲労度が生じると、筋肉内の水分量が変化し、筋肉内の水分量の変化を、生体インピーダンスの変化として捉えることによって、特定の筋肉の疲労度を評価することができる。

【0062】

そして、筋肉の疲労度を評価するための生体インピーダンスは、前述のように、第１の状態及び第２の状態の間の切り替えに際して取得した過渡応答の情報に基づいて算出することができる。すなわち、本実施形態では、取得した過渡応答の情報に基づいて、生体の状態としての筋疲労度を評価することができる。

【0063】

図４は、筋肉内の水分量の変化（図４（ａ））と、生体インピーダンスの変化（図４（ｂ））との関係を模式的に示したグラフである。図４（ａ）、（ｂ）に示すように、筋肉に負荷をかけて筋疲労が生じると、筋肉内の水分量が増加し（時刻ｔ_０～ｔ_１）、これに伴い、生体インピーダンスが減少する（時刻ｔ_０～ｔ_１）。そして、筋肉に負荷をかけるのを止めると、筋肉内の水分量が元の状態に戻り（時刻ｔ_１～ｔ_２）、これに伴い、生体インピーダンスも元の状態に戻る（時刻ｔ_１～ｔ_２）。すなわち、生体インピーダンスの時間変化を測定することにより、筋疲労が、時刻ｔ_０～ｔ_１の間で蓄積し、時刻ｔ_１～ｔ_２の間で回復することが分かる。これにより、筋肉に負荷をかけたときに生じる筋疲労度を評価することができる。

【0064】

ここで、時間変化のモニタ対象とする生体インピーダンスとしては、式（７）で表される生体インピーダンスＺbの絶対値を用いてもよいし、インピーダンスＺbの実部（レジスタンス）又は虚部（リアクタンス）の値を用いてもよいし、インピーダンスＺbの実部及び虚部の値を同時にモニタしてもよい。

【0065】

なお、図示は省略するが、運動していない場合にも、生体インピーダンスは、筋疲労以外の要因で、時間変化する場合がある。そのため、生体インピーダンスの時間変化に基づいて、筋疲労を評価する場合、筋疲労以外の要因による時間変化と区別する必要がある。

【0066】

通常、運動していない場合の生体インピーダンスの時間変化量（傾き）は、運動した場合の生体インピーダンスの時間変化量（傾き）よりも小さい。典型的には、運動前の生体インピーダンスに対して、運動直後の生体インピーダンスが所定の割合以上変化した場合は、筋疲労によるものと考えられる。

【0067】

一方、生体インピーダンスの時間変化が、前記所定の割合以下の場合は、筋疲労によるものではないと考えられる。従って、運動後に算出した生体インピーダンスの時間変化量が、所定の値以上になったときを、筋疲労と判断することによって、筋疲労を正しく評価することができる。

【0068】

このように、本実施形態に係る生体状態評価方法では、過渡応答ＲＴの収束後に電位を測定するのではなく、その過渡応答ＲＴの情報それ自体を活用することで生体の情報を評価する。過渡応答ＲＴの収束を待つ必要がない分、測定時間を短縮することができる。

【0069】

また、過渡応答の情報として整定時間を取得することで、その整定時間に対応した時定数を算出することができる。そうして算出した時定数は、いわゆる抵抗値（レジスタンス）ばかりでなく、静電容量（キャパシタンス）の情報も含んでいる。したがって、抵抗値のみを用いた評価方法と比較して、より高精度に生体の情報を評価したり、より多方面からの評価を行ったりすることができるようになる。

【0070】

また、式（７）及び式（８）に示すように、電気回路Ｃの切り替えを利用して生体インピーダンスＺｂを算出する際に、実部及び虚部の双方とも、従来用いられてきた電圧比を必要としない。これにより、測定時間の短縮を図る上で有利になる。

【0071】

また、一方の外部抵抗の抵抗値を無限大に設定することで、第１の状態と第２の状態との間の切り替えを、電気回路Ｃが接続された状態と電気回路Ｃが遮断された状態との間の切り替えとみなすことができる。これにより、生体インピーダンスＺbの算出に要する数式を、式（８）に示すように、よりシンプルなものにすることができる。そのことで、生体の状態をより容易に評価することが可能になる。

【0072】

＜他の実施形態＞
以上、本開示を好適な実施形態により説明してきたが、こうした記述は限定事項ではなく、もちろん、種々の改変が可能である。

【0073】

例えば、前記実施形態では、生体の情報として、当該生体の筋疲労度を評価するように構成されていたが、本開示における評価対象は、生体の筋疲労度には限定されない。生体の情報として、生体内の信号源の位置を評価したり、生体の体脂肪を評価したりしてもよい。

【0074】

例えば、特許文献１により開示されている生体内信号源位置検出方法は、生体内の信号源と電極との間の内部抵抗の抵抗値を用いたものであるが、その抵抗値として、過渡応答の情報によって得られる値（例えば、上式（５）を参照）を用いてもよい。

【0075】

また前記実施形態では、図３に示すような出力電圧Ｖoutの時間変化に基づいて過渡応答ＲＴの情報を取得するように構成されていたが、本開示は、そうした構成には限定されない。例えば、生体内から出力される電流の時間変化に基づいて過渡応答ＲＴの情報を取得してもよい。

【0076】

また、生体の体脂肪等を評価対象とした場合、生体表面に電極を配置して、その電極から生体に高周波電流を流すことによって、生体インピーダンスを測定してもよい。

【0077】

また、前記実施形態では、生体表面３０に２つの電極１０、２０を配置したが、生体表面３０にグランド電極を配置し、２つの電極１０、２０間に生じた電圧を、差動アンプ４０によって測定してもよい。この場合も、２つの電極１０、２０間における生体インピーダンスは、前記の式（７）により求めることができる。また、第１の電圧Ｖ_１、及び第２の電圧Ｖ_２は、その差分をとって差動アンプ４０で増幅されて測定されるため、外部からのノイズを除去することができる。これにより、生体インピーダンスＺｂをより精度よく測定することができる。

【0078】

また、前記実施形態では、筋繊維からなる筋肉部位（上腕筋）の筋疲労を評価するとき、２つの電極を筋繊維に沿って互いに近傍に配置したが、本開示は、そうした配置には限定されない。２つの電極を、筋繊維を取り囲むように、互いに近傍に配置してもよい。

【0079】

また、本開示は、電極を有する測定装置を用いた生体状態評価システムとして使用することもできる。すなわち、本開示に係る生体状態評価システムは、生体表面に所定の間隔を空けて配置される電極と、その電極を介して測定装置と生体とで形成される電気回路を第１の状態及び第２の状態の間で変化させる切り替え手段と、その切り替え手段により、第１の状態及び第２の状態の間で変化させたときに電気回路に生じる過渡応答の情報を取得する情報取得手段と、取得した過渡応答の情報に基づいて、生体の状態を評価する生体状態評価手段とを備えている。

【符号の説明】

【0080】

１ 測定装置
１０ 電極
２０ 電極
３０ 生体表面
Ｃ 電気回路
Ｒg1 第１の外部抵抗
Ｒg2 第２の外部抵抗

【図1】

【図2】

【図3】

【図4】