(57)【特許請求の範囲】
【請求項１】
心電図センサ、心音図センサ、脈波センサ及び超音波センサを含む生体信号測定センサ部と、前記生体信号測定センサ部の各センサに連結されて測定された生体信号の受信を受けて信号処理する生体信号受信及び処理部とで構成された生体信号測定計、及び
前記生体信号受信及び処理部と連結されて互いに通信して測定データの受信を受けて脳血管を分析するための生物力学的指標を算出する主処理部と、前記主処理部に連結されて使用者の制御命令の入力を受ける入力部と、前記主処理部に連結されて算出された結果を表示する出力部とで構成された分析指標算出計を含んで構成され、
前記脈波センサは、カフ（Ｃｕｆｆ）脈波センサ、頚動脈脈波センサ及び大腿動脈脈波センサのうち何れか一つであり、
前記生体信号測定計は、前記生体信号測定センサ部からＥＣＧ、ＰＣＧ及びＡＰＧ波形を同時に得て、
前記主処理部は、前記生体信号測定計から獲得した生体信号により得た脳血管の圧力曲線、前記圧力曲線に関する収縮期と弛緩期の面積及び前記脳血管の血流量を含む基礎情報により生物力学的指標を算出して、
前記主処理部は、
前記入力器により基礎情報の入力を受け、前記生体信号測定計を通して前記生体信号の受信を受ける第１段階、
前記受信された生体信号の波形を分析し、前記分析された波形資料をもとに脳血管の圧力曲線、前記圧力曲線に関する収縮期と弛緩期の面積及び前記脳血管の血流量を得る第２段階、及び
前記得た圧力曲線、圧力曲線の面積及び血流量と前記入力された基礎情報をもとに前記脳血管の順応性Ｃ及び抵抗Ｒを含む生物力学的指標を算出し、脳血管の分析結果を表示する第３段階が含まれるようにプログラムされることを特徴とする脳血管分析装置。
【請求項２】
前記生体信号受信及び処理部は、
前記生体信号測定センサ部から受信される生体信号を処理し、前記主処理部に測定データを伝送するように制御する制御器（ｍｉｃｒｏｃｏｎｔｒｏｌｌｅｒ）、
前記制御器の制御信号により前記心電図センサ、心音図センサ、脈波センサ及び超音波センサから受信される生体信号を選択する多重信号選択器、
前記多重信号選択器により選択された生体信号を前記制御器の制御信号によって雑音除去又は増幅も調節する雑音除去及び信号増幅器、
前記雑音除去及び信号増幅器を通過した生体信号を受け、前記入力部の制御命令又は前記主処理部に内蔵されたプログラムの制御命令が、前記制御器を通して必要な生体信号の選択を受けるようにする信号切換器、
前記信号切換器で選択された生体信号を前記制御器の制御信号によってサンプリング（ｓａｍｐｌｉｎｇ）してホールディング（ｈｏｌｄｉｎｇ）する標本維持器、及び
前記標本維持器を通してホールディング（ｈｏｌｄｉｎｇ）された生体信号を、前記制御器の制御信号によってデジタル信号に変えて前記制御器に送るＡ／Ｄ変換器を含んで構成されたことを特徴とする請求項１に記載の脳血管分析装置。
【請求項３】
前記頸動脈脈波センサと前記大腿動脈脈波センサは同一の圧力センサであり、
前記カフ（Ｃｕｆｆ）脈波センサは、カフ血圧計に圧力センサが更に付着されたことを特徴とする請求項１又は請求項２に記載の脳血管分析装置。
【請求項４】
前記カフ（Ｃｕｆｆ）脈波センサは、前記カフ血圧計の空気袋と連結されたゴム管に支路管を形成し、前記支路管の出口にアダプタを装着し、前記アダプタを前記頸動脈脈波センサ又は前記大腿動脈脈波センサと同一の構造のセンサ開放溝に装着したことを特徴とする請求項３に記載の脳血管分析装置。
【請求項５】
前記第２段階の前記脳血管は、左、右側後大脳動脈を含み、
前記第３段階における前記生物力学的指標中の左、右側後大脳動脈の順応性（Ｃｐ１、Ｃｐ２）及び左、右側後大脳動脈の抵抗（Ｒｐ１、Ｒｐ２）は、下記数式によりそれぞれ算出されることを特徴とする請求項１に記載の脳血管分析装置。
左側後大脳動脈の順応性
【数１】
0005628180.tif 000081
右側後大脳動脈の順応性
【数２】
0005628180.tif 000082
左側後大脳動脈の抵抗
【数３】
0005628180.tif 000083
及び
右側後大脳動脈の抵抗
【数４】
0005628180.tif 000084
前記数式１乃至４において、Ｐｐ１ｓは左側後大脳動脈の収縮期血圧、Ｐｐ１ｄは左側後大脳動脈の弛緩期血圧、Ｐｐ２ｓは右側後大脳動脈の収縮期血圧、Ｐｐ２ｄは右側後大脳動脈の弛緩期血圧、Ａｐ１ｓは左側後大脳動脈圧力曲線の収縮期面積、Ａｐ１ｄは左側後大脳動脈圧力曲線の弛緩期面積、Ａｐ２ｓは右側後大脳動脈圧力曲線の収縮期面積、Ａｐ２ｄは右側後大脳動脈圧力曲線の弛緩期面積、Ｓｐ１は左側後大脳動脈の血流量、Ｓｐ２は右側後大脳動脈の血流量である。
【請求項６】
前記第２段階の前記脳血管は、左、右側前大脳動脈を含み、
前記第３段階における前記生物力学的指標中の左、右側前大脳動脈の順応性（Ｃａ１、Ｃａ２）及び左、右側前大脳動脈の抵抗（Ｒａ１、Ｒａ２）は、下記数式によりそれぞれ算出されることを特徴とする請求項１に記載の脳血管分析装置。
左側前大脳動脈の順応性
【数５】
0005628180.tif 000085
右側前大脳動脈の順応性
【数６】
0005628180.tif 000086
左側前大脳動脈の抵抗
【数７】
0005628180.tif 000087
及び
右側前大脳動脈の抵抗
【数８】
0005628180.tif 000088
前記数式５乃至８において、Ｐａ１ｓは左側前大脳動脈の収縮期血圧、Ｐａ１ｄは左側前大脳動脈の弛緩期血圧、Ｐａ２ｓは右側前大脳動脈の収縮期血圧、Ｐａ２ｄは右側前大脳動脈の弛緩期血圧、Ａａ１ｓは左側前大脳動脈圧力曲線の収縮期面積、Ａａ１ｄは左側前大脳動脈圧力曲線の弛緩期面積、Ａａ２ｓは右側前大脳動脈圧力曲線の収縮期面積、Ａａ２ｄは右側前大脳動脈圧力曲線の弛緩期面積、Ｓａ１は左側前大脳動脈の血流量、Ｓａ２は右側前大脳動脈の血流量、Ｋは臨床係数である。
【請求項７】
前記第２段階の前記脳血管は、左、右側中大脳動脈を含み、
前記第３段階における前記生物力学的指標中の左、右側中大脳動脈の順応性（Ｃｍ１、Ｃｍ２）及び左、右側中大脳動脈の抵抗（Ｒｍ１、Ｒｍ２）は、下記数式によりそれぞれ算出されることを特徴とする請求項１に記載の脳血管分析装置。
左側中大脳動脈の順応性
【数９】
0005628180.tif 000089
右側中大脳動脈の順応性
【数１０】
0005628180.tif 000090
左側中大脳動脈の抵抗
【数１１】
0005628180.tif 000091
及び
右側中大脳動脈の抵抗
【数１２】
0005628180.tif 000092
前記数式９乃至１２において、Ｐｍ１ｓは左側中大脳動脈の収縮期血圧、Ｐｍ１ｄは左側中大脳動脈の弛緩期血圧、Ｐｍ２ｓは右側中大脳動脈の収縮期血圧、Ｐｍ２ｄは右側中大脳動脈の弛緩期血圧、Ａｍ１ｓは左側中大脳動脈圧力曲線の収縮期面積、Ａｍ１ｄは左側中大脳動脈圧力曲線の弛緩期面積、Ａｍ２ｓは右側中大脳動脈圧力曲線の収縮期面積、Ａｍ２ｄは右側中大脳動脈圧力曲線の弛緩期面積、Ｓｍ１は左側中大脳動脈の血流量、Ｓｍ２は右側中大脳動脈の血流量、Ｓａ１は左側前大脳動脈の血流量、Ｓａ２は右側前大脳動脈の血流量、Ｋは臨床係数である。
【請求項８】
前記主処理部は、前記第３段階で算出された前記脳血管の順応性Ｃ及び抵抗Ｒを前記出力部にＣ−Ｒ状態図（Ｃｈａｒｔ）上の点で表示されるように制御することを特徴とする請求項１乃至請求項７の何れか１項に記載の脳血管分析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血管の生物力学的特性と血流特性を測定し、脳血管の疾患を分析する血管分析システムに関し、より詳しくは、脳血管系を内頚動脈枝、前大脳血管枝、中大脳血管枝、後大脳血管枝、椎骨動脈枝及び脳底動脈枝からなる系に区分し、これを分析して脳血管系の血管枝ごとで血管の気質的変化を示す血管の弾性係数を求め、脳血管系の血流特性、気質的変化と機能的変化とを同時に示す血管の順応性と血流抵抗、血流量を測定し、脳血栓を始めとする様々な脳血管の難治の病を早期診する脳血管分析置に関する。
【背景技術】
【０００２】
今臨床においては、脳血管の疾患を早期診断するために、超音波ドップラー計を用いている。しかし、超音波ドップラー計は、血管特性を測定することができないため、脳血管の疾患を早期診断するには、一定の制限を有している。
【０００３】
今まで開発された脳血管の疾患の分析機としては、血管造影機、ＭＲＡ、ＦＭＲＩ、ＳＰＥＴ、ＴＣＤ、ＴＥＥ、ＴＴＥ、ＱＦＭ、ＣＶＤ等がある。
【０００４】
この中で、血管造影機は、血管自体の病的変化を直感的に観察することができる長所はあるが、造影剤の注入のための血管侵襲的な操作が必須でかつ操作が複雑である。
ＭＲＡやＦＭＲＩは、血管造影機の欠陥を克服した体系ではあるが、製作原価と診断費用が高価であるため、特定な病棟でのみ利用することができるという欠陥がある。
【０００５】
特に、ＭＲＡ、ＦＭＲＩ、ＳＰＥＴは、互いに多少の差があるが、これらの設備による検査は、血管分布、血流特性、低血流区域等を確認することができる反面、血管の特性を確認することができない。
【０００６】
超音波血流測定機（ultrasonic Quanti tative Flow Measurement System：ＱＦＭ）と脳血管特性測定機（ＣＶＤ）は、安い費用で頚動脈の血流量と、中大脳血管と前大脳血管の順応性を求めることができる。
【０００７】
しかし、血管の特性を特徴化する血管の気質的状態及び機能的状態を評価するには、血管に流れる血流量や血管壁に作用する血圧のような血流状態と関連する情報よりは、血管自体の気質的状態と機能的状態を反映する血管の弾性係数、血管の順応性と血流抵抗等の情報を認知することが必要である。
【０００８】
しかし、血管の気質状態を反映する脳血管系の血管枝ごとで血管の弾性係数と血管の順応性と直径、血流抵抗、血流量を測定することは、極めて難しい問題である。
【０００９】
これは、脳血管系の構造が複雑であり、脳血管系で血管枝は互いに異なる力学的挙動をするのにも起因するが、生きている人体で血管の気質的変化を示す指標である弾性係数を測定するには、生体内の血管を採集して引っ張り器で引っ張ることが最も正確な方法ではあるが、そのようなことは不可能であるためである。
【００１０】
２００２年に米国のカルフォニア綜合大学教授らであるWerner G、Marifan C、Tonny M、Jeffrey C等の複数の学者たちがヒトの脳血管の生物力学的特性を研究し、‘Mechanical and Failure Propertis of Human Cerebral Blood Vesseies’に脳血管の特性に関する論文を発表したことがある。
【００１１】
しかし、これらの結果は、ヒトの血管を採集して実験したものであり、診断を行うのに用いることのできる指標としては利用することができない。
【００１２】
脳血管系の血流量と順応性、弾性係数、血流抵抗等を間接的に求めるための研究は多様に行われていた。
【００１３】
１９９７〜２００４年の間米国学者であるBiedma、Haoliu、Cwako shiを始めとする多くの学者たちが脳血管の順応性と抵抗を求める問題を提案した。
【００１４】
しかし、以上の研究結果は、脳血管で血圧、血管の順応性と弾性係数、血流抵抗、血流量の分布に関する一般的な開示のみを盛り込んでいるだけであり、臨床に直接適用することができる結果は得られなかった。
【００１５】
臨床に適用する目的で、２００６年に丁光宏（Ding Guanghong）教授（上海▲復▼旦大学部校長）は、脳血管の個別的な枝の血流量を求める設備であるＫＦ−３０００を開発している。
【００１６】
ＫＦ−３０００は、脳血管系で血流特性を求めるのに革新的な成果を上げたが、血管の特性までは把握することができず、ＴＣＤを脳血管の疾病を早期診断するまでは発展させることができなかった。
【００１７】
丁光宏博士の研究とは異なり、日本ＣＯＬＬＩＮは、脳血管の疾病を早期診断するのに決定的意義を持つ脳血管系の血流特性と血管特性を評価する脳血流量計ＱＦＭ−２０００Ｘを提案し、類似する装置として、特許出願された‘脳血管のパラメーター測定装置及び方法’（An apparatus measuring parameter of cerebrovascular and method thereof）に基づくＣＶＤ−１０００も提案されている。
【００１８】
臨床に適用する目的で開発された超音波血流測定機（ultrasonic Quanti tative Flow Measurement System：ＱＦＭ）と、特許出願された‘脳血管のパラメーター測定装置及び方法’による構成は、安い費用で内頚動脈の血流量と中大脳血管と、前大脳血管の順応性を求めることにより、脳血管疾病を早期診断することができる可能性を示した。
【００１９】
しかし、超音波血流測定機（ＱＦＭ）と上記特許の構成は、基本測定指標として血管の順応性と血流指標とを選定することにより、血管特性の２つの基本特性である血管の気質的変化と機能的変化を区分して究明することはできなかった。
【００２０】
特に、超音波血流測定機（ＱＦＭ−２０００Ｘ）と脳血管特性測定機（ＣＶＤ−１０００）は、脳血管の順応性と抵抗を求めるには、様々な欠点を示していた。
【００２１】
ＱＦＭ−２０００Ｘは、脳血管系を左右に分離して内頚動脈に入り込む血流量を脳血流量と仮定し、左側脳血管系と右側脳血管系を評価する順応性Ｃと抵抗Ｒを求めた。
【００２２】
従って、脳の個別的血管枝を評価することができなかった。
【００２３】
また、これらの構成は、脳血管の順応性と抵抗を求めるために、血圧波形を圧力脈波に、血流波形を超音波ドップラーに与え、両波形が一致するまで校正してＣ、Ｒを求めたが、Ｃ、Ｒの解は、変動振幅が大きく、また曲線近似が実際現象と大きく差が出る欠点があった。
【００２４】
実際、超音波ドップラーにより血流波形を測るとき、水平面の誤差のために測定誤差が非常に大きい。ところが、そのような測定誤差を持つ波形を置いて両波形を近似させることは事実と合わず再現性が少ない。また、曲線対曲線の近似は、若干の波形変動にも両波形の近似時に大きい誤差を引き起こす。
【００２５】
従って、ＱＦＭ−２０００Ｘにおいて、順応性Ｃや抵抗Ｒは測定するヒトによってほぼ１０〜１００倍までの差をもたらして再現性がなく、臨床指標として用いるのが困難である。
【００２６】
特許出願された脳血管のパラメーター測定装置及び方法の構成は、脳血管系を模型化して分析するとき、脳に入り込む血流量を求めずに、脳に流れる血流量が心拍出量のｋ倍になるものと仮定することにより、正確な臨床指標を得ることができなかった。特許出願された脳血管のパラメーター測定装置及び方法の構成においては、次のように開示している。
【００２７】
内頚動脈の断面積が８０〜９０％減っても内頚動脈に流れ込む血流量には変化がないので、これに基づいて内頚動脈に流れ込む血流量は、
Ｑｃ=ＫｃＳｖ
として計算することができる。ここで、Ｓｖは心拍出量、Ｋｃは比例定数である。
【００２８】
しかし、上記のような仮定は、医療診断機器として十分な担保になり得ない。
また、特許出願された脳血管のパラメーター測定装置及び方法の構成は、前大脳血管と中大脳血管の順応性と抵抗が一定の比率で分けられる仮定することにより、疾病診断において正確度を落とした。
【００２９】
特に、ＱＦＭ−２０００Ｘと特許出願された脳血管のパラメーター測定装置及び方法は、後大脳血管の順応性と抵抗を求める方法を提示することができなかった。
【００３０】
このように、ＱＦＭ−２０００Ｘと上記特許された脳血管のパラメーター測定装置及び方法は、脳血管の特性を評価するために、弾性係数は求めることができずに、順応性と抵抗のみを求めたが、それ自体も多くの欠陥を持っている。
【００３１】
従って、脳血管系の弾性係数、順応性、抵抗、血流量を正確に求める新たな解法に関する要求が存在する。
【発明の概要】
【発明が解決しようとする課題】
【００３２】
本発明は、前述の従来技術の問題点を解消するために導き出されたものであり、脳血管の疾患を診断するために、脳血管の生物力学的特性と血流特性を測定するように心電図、心音図、脈波信号及び超音波ドップラー信号を基礎資料として脳血管系を生物力学的に分析することにより、脳血管系の血管枝ごとに弾性係数、血管の順応性と血流抵抗、血流量を測定し、脳血栓を始めとする様々な脳血管の難治の病を早期診断することができる脳血管分析装置を提供することを目的とする。
【課題を解決するための手段】
【００３３】
前述の目的を達成するために、本発明による脳血管分析装置は、心電図センサ、心音図センサ、脈波センサ及び超音波センサを含む生体信号測定センサ部と、生体信号測定センサ部の各センサに連結されて測定された生体信号の受信を受けて信号処理する生体信号受信及び処理部とで構成された生体信号測定計、及び生体信号受信及び処理部と連結されて互いに通信して測定データの受信を受けて脳血管を分析するための生物力学的指標を算出する主処理部と、主処理部に連結されて使用者の制御命令の入力を受ける入力部と、主処理部に連結されて算出された結果を表示する出力部とで構成された分析指標算出計を含んで構成され、主処理部は、生体信号測定計から獲得した生体信号により得た脳血管の圧力曲線、圧力曲線に関する収縮期と弛緩期の面積及び脳血管の血流量を含む基礎情報により生物力学的指標を算出することを第１の特徴とする。
【００３４】
そして、生体信号受信及び処理部は、生体信号測定センサ部から受信される生体信号を処理し、主処理部に測定データを伝送するように制御する制御器（microcontroller）、制御器の制御信号により心電図センサ、心音図センサ、脈波センサ及び超音波センサから受信される生体信号を選択する多重信号選択器、多重信号選択器により選択された生体信号を制御器の制御信号によって雑音除去又は増幅も調節する雑音除去及び信号増幅器、雑音除去及び信号増幅器を通過した生体信号を受け、入力部の制御命令又は主処理部に内蔵されたプログラムの制御命令が、制御器を通して必要な生体信号の選択を受けるようにする信号切換器、信号切換器で選択された生体信号を制御器の制御信号によってサンプリング（sampling）してホールディング（holding）する標本維持器、及び標本維持器を通してホールディング（holding）された生体信号を、制御器の制御信号によってデジタル信号に変えて制御器に送るＡ／Ｄ変換器を含んで構成されたことを第２の特徴とする。
【００３５】
そして、脈波センサは、カフ（Ｃｕｆｆ）脈波センサ、頚動脈脈波センサ及び大腿動脈脈波センサのうち何れか一つであり、生体信号測定計は、生体信号測定センサ部からＥＣＧ、ＰＣＧ及びＡＰＧ波形を同時に得ることを第３の特徴とする。
【００３６】
そして、頚動脈脈波センサと大腿動脈脈波センサは同一の圧力センサであり、カフ（Ｃｕｆｆ）脈波センサは、カフ血圧計に圧力センサが更に付着されたことを第４の特徴とする。
【００３７】
そして、カフ（Ｃｕｆｆ）脈波センサは、カフ血圧計の空気袋と連結されたゴム管に支路管を形成し、支路管の出口にアダプタを装着し、アダプタを頚動脈脈波センサ又は大腿動脈脈波センサと同一の構造のセンサ開放溝に装着したことを第５の特徴とする。
【００３８】
そして、主処理部は、入力部により基礎情報の入力を受け、生体信号測定計を通して生体信号の受信を受ける第１段階、受信された生体信号の波形を分析し、分析された波形資料をもとに脳血管の圧力曲線、圧力曲線に関する収縮期と弛緩期の面積及び脳血管の血流量を得る第２段階、及び得た圧力曲線、圧力曲線の面積及び血流量と入力された基礎情報をもとに脳血管の順応性Ｃ及び抵抗Ｒを含む生物力学的指標を算出し、脳血管の分析結果を表示する第３段階が含まれるようにプログラムされたことを第６の特徴とする。
【００３９】
そして、第２段階の脳血管は、左、右側後大脳動脈を含み、第３段階における生物力学的指標中の左、右側後大脳動脈の順応性（Ｃｐ１、Ｃｐ２）及び左、右側後大脳動脈の抵抗（Ｒｐ１、Ｒｐ２）は、所定の数式によりそれぞれ算出されたことを第７の特徴とする。
【００４０】
そして、第２段階の脳血管は、左、右側前大脳動脈を含み、第３段階における生物力学的指標中の左、右側前大脳動脈の順応性（Ｃａ１、Ｃａ２）及び左、右側前大脳動脈の抵抗（Ｒａ１、Ｒａ２）は、所定の数式によりそれぞれ算出されたことを第８の特徴とする。
【００４１】
そして、第２段階の脳血管は、左、右側中大脳動脈を含み、第３段階における生物力学的指標中の左、右側中大脳動脈の順応性（Ｃａ１、Ｃａ２）及び左、右側中大脳動脈の抵抗（Ｒａ１、Ｒａ２）は、所定の数式によりそれぞれ算出されたことを第９の特徴とする。
【００４２】
そして、主処理部は、第３段階で算出された脳血管の順応性Ｃ及び抵抗Ｒを出力部にＣ−Ｒ状態図（Ｃｈａｒｔ）上の点で表示されるように制御することを第１０の特徴とする。
【発明の効果】
【００４３】
本発明は、脳血管の個別的血管枝の気質的変化を反映する弾性係数を分析し、脳血管系の血流特性と脳血管の個別的血管枝の気質的変化と機能的変化を反映する血流量、血管の順応性、血流抵抗を算出することにより、脳血管の疾患の発生の危険性を早期に診断することができる。
【図面の簡単な説明】
【００４４】
【図１】本発明による脳血管分析装置の一実施例による全体システムのブロック図である。
【図２】図１の生体信号受信及び処理部の構成と信号流れを概念的に示すブロック図である。
【図３】図１の脈波センサの一種であるカフ脈波センサの構成を示す一面図及び要部分解斜視図である。
【図４】脳血管の血管枝の連結状態を概念的に示すウィリス輪（Willi's circle）である。
【図５】内頚動脈で前大脳動脈と中大脳動脈が分岐される場合であると仮定して模型化させた図４による脳血管の模型化回路図である。
【図６】内頚動脈と中大脳動脈を１つの血管枝と仮定して模型化させた図４による脳血管の模型化回路図である。
【図７】図１の主処理部の一例示的作業図である。
【図８】図１の主処理部による分析結果の一例としてＣ−Ｒ状態図（Ｃｈａｒｔ）を示す一例示図である。
【発明を実施するための形態】
【００４５】
以下、本発明が属する技術分野における通常の知識を有する者が容易に実施することができるように、本発明の実施例について添付した図面を参考として詳細に説明する。しかし、本発明は、様々な相違する形態に具現することができ、ここで説明する実施例に限定されない。
【００４６】
図１は本発明による脳血管分析装置の一実施例による全体システムのブロック図であり、図２は図１の生体信号受信及び処理部の構成と信号流れを概念的に示すブロック図であり、図３は図１の脈波センサの一種であるカフ脈波センサの構成を示す一面図及び要部分解斜視図であり、図４は脳血管の血管枝の連結状態を概念的に示すウィリス輪（Willi's circle）である。図５は内頚動脈で前大脳動脈と中大脳動脈が分岐される場合であると仮定して模型化させた図４による脳血管の模型化回路図であり、図６は内頚動脈と中大脳動脈を１つの血管枝と仮定して模型化させた図４による脳血管の模型化回路図であり、図７は図１の主処理部の一例示的作業図であり、図８は図１の主処理部による分析結果の一例としてＣ−Ｒ状態図（Ｃｈａｒｔ）を示す一例示図である。
【００４７】
本発明による脳血管分析装置の一実施例は、基本的に図１のように、心電図センサ１２２、心音図センサ１２４、脈波センサ１２６及び超音波センサ１２８を含む生体信号測定センサ部１２０と、生体信号測定センサ部１２０の各センサに連結されて測定された生体信号の受信を受けて信号処理する生体信号受信及び処理部１４０とで構成された生体信号測定計１００、及び生体信号受信及び処理部１４０と連結されて互いに通信して測定データの受信を受けて脳血管を分析するための生物力学的指標を算出する主処理部２１０と、主処理部２１０に連結されて使用者の制御命令の入力を受ける入力部２２０と、主処理部２１０に連結されて算出された結果を示す出力部２４０とで構成された分析指標算出計２００を含んで構成され、主処理部２１０は生体信号測定計１００から得られた生体信号により得た脳血管の圧力曲線、圧力曲線に関する収縮期と弛緩期の面積及び脳血管の血流量を含む基礎情報により生物力学的指標を算出することを特徴とする。
【００４８】
ここで、心電図センサ１２２は、少なくいとも３つ以上の電極で構成され、心電図（Electrocardiogram：ＥＣＧ）波形を得るためのものであり、心音図センサと共に脳血管圧力曲線の特徴点（収縮期の開始点、収縮期の最高点、切痕点、弛緩期の最高点、弛緩期の終了点）を把握するためのものである。
【００４９】
心音図センサ１２４は、心臓の弁膜が開閉するときに出る音を感知するためにマイクロフォンで構成され、これを通して心音図（Phonocardiogram：ＰＣＧ）波形を得て脳血管圧力曲線の特徴点を把握するためのものである。
【００５０】
脈波センサ１２６は、脈動による脈波を感知してＡＰＧ（Accelerated Plethysmogram：ＡＰＧ）波形を得るためのものであり、圧電素子で構成された圧力センサを用いることができるが、脈動を感知することができれば、これに限定されない。
【００５１】
本実施例による脈波センサ１２６は、脳血管の頻度スペクトル情報を得るためのカフ（Ｃｕｆｆ）脈波センサ、左、右側頚動脈の脈波を直接測定して脳血管の確率密度スペクトル情報を得るための頚動脈センサ、大腿動脈の脈波を直接測定して脈波伝達速度（Plse Wave Velocity：ＰＷＶ）等を求めるための大腿動脈センサのうち何れか１つを指す。
【００５２】
ここで、頚動脈脈波センサと大腿動脈脈波センサは同一種類の圧力センサであり得、カフ（Ｃｕｆｆ）脈波センサはカフ血圧計に圧力センサが更に付着されたもので構成することができる。
【００５３】
図３は、カフ（Ｃｕｆｆ）脈波センサの具体的構成の例を示すものであり、これによると、従来カフ血圧計１０の空気袋１３と連結されたゴム管（１４又は１７）に支路管２１を形成し、支路管の出口にアダプタ２１を装着し、アダプタ２１を頚動脈脈波センサ又は大腿動脈脈波センサと同一の構造のセンサ（例えば、圧力センサ、３４)の開放溝３２に装着して用いることができる。
【００５４】
超音波センサ１２８は、超音波をよく通過させることができる頭蓋骨の位置に探触子（probe）を置いて超音波を発射させて反射波を分析することにより、脳血管の血流速度及び血流量を測定するためのものであり、頭蓋骨貫通ドップラー（Trans-Carnial Doppler：ＴＣＤ）とも言う。
【００５５】
以上のように、心電図センサ１２２、心音図センサ１２４、脈波センサ１２６及び超音波センサ１２８は、別個の生体信号を感知するための生体信号測定センサ部１２０の必須構成となり、生体信号測定センサ部１２０と連結される生体信号受信及び処理部１４０が内蔵された装置には、少なくとも４つの連結端子が備えられる。
【００５６】
そして、生体信号受信及び処理部１４０は、図２のように、生体信号測定センサ部１２０から受信される生体信号を処理して主処理部２１０に測定データを伝送するように制御する制御器（microcontroller：マイコン、１４８）、制御器１４８の制御信号により心電図センサ１２２、心音図センサ１２４、脈波センサ１２６及び超音波センサ１２８から受信される生体信号を選択する多重信号選択器１４１、多重信号選択器１４１により選択された生体信号を制御器１４８の制御信号によって雑音除去又は増幅も調節する雑音除去及び信号増幅器１４２、雑音除去及び信号増幅器１４２を通過した生体信号を受けて入力部２２０の制御命令又は主処理部２１０に内蔵されたプログラムの制御命令が、制御器１４８を通して必要な生体信号を選択されるようにする信号切換器１４３、信号切換器１４３で選択された生体信号を制御器１４８の制御信号によってサンプリング（sampling）してホールディング（holding）する標本維持器１４４、及び標本維持器１４４を通してホールディング（holding）された生体信号を制御器１４８の制御信号によってデジタル信号に変えて制御器１４８に送るＡ／Ｄ変換器１４５を含んで構成される。
【００５７】
ここで、多重信号選択器１４１は、心電図センサ１２２、心音図センサ１２４、脈波センサ１２６及び超音波センサ１２８により同時測定して同時入力されるとき、これらの測定信号を順次選別して処理するためのものであり、雑音除去及び信号増幅器１４２は、得られた生体信号で各種雑音を除去して標準波とし、患者（被験者）によって増幅度を調節することができるように備えられる。
【００５８】
上記のように、生体信号受信及び処理部１４０は、生体信号測定計１００に含まれるように備えられるのが好ましいが、回路設計によって後述の主処理部２１０と一体として構成することもできる。
【００５９】
次に、生体信号測定計１００で得られて処理された生体信号は分析指標算出計２００に伝達され、分析指標算出計２００で脳血管圧力曲線を算出し、脳血管圧力曲線の面積及び血流量の情報を利用して生物力学的指標を算出するようになる。
【００６０】
生体信号受信及び処理部１４０が、図１のように、主処理部２１０と離れて構成されるときは、両者間の所定の通信手段（例えば、ＲＳ−２３２Ｃ）によりデータをやり取りするようになる。
【００６１】
主処理部２１０は、内蔵型メモリ部又は外付型メモリ部に保存されたプログラムによって生体信号受信及び処理部１４０から伝達された測定データを処理して脳血管を分析するための生物力学的指標を算出する核心装置であり、コンピュータの中央処理装置に当たる。
【００６２】
ここで、脳血管を分析するための生物力学的指標は、各脳血管枝の血流量、順応性、血流抵抗、動脈硬化度及び血流速度を言う。
【００６３】
先ず、各生物力学的指標が、本明細書で用いる定義と他の指標間の関係について簡単に説明する。
【００６４】
血流量は、脳血管枝に沿って流れる血液量を言うものであって単位はｍＬであり、時間の関数で表現する場合はＱ又はＱ（ｔ）であり、一定時間流れた血液量（Ｑの時間積分量）で表現する場合はＳで表記する。血流量は、一般的に脳血管枝の長さ方向に離隔された２ヶ所の血圧差（Ｐ−Ｐｖ）に正比例し、血流抵抗（Ｒ）に反比例する。血流量が少ないと、それによる虚血症状等が現れるようになる。
【００６５】
順応性（Complience）は、単位体積の血管に単位力を与えたときに起こる体積変化を言うものであって単位はｍＬ／ｍｍＨｇであり、簡略にＣと表記する。Ｃが小さいとは血管壁が硬化されるか、又は収縮されることを示し、反対にＣが大きいとは血管壁が柔軟であるか、拡張型痙攣が起こることを意味する。
【００６６】
血流抵抗（Resistance）は、脳血管枝に沿って流れる血液が受ける抵抗を言うものであって単位はｍｍＨｇ／Ｌであり、簡略にＲで表記される。Ｒは近似的に脳血管枝の長さ方向に離隔された２ヶ所の血圧差（Ｐ−Ｐｖ）と血流量（Ｑ）の比で決定される。
【００６７】
動脈硬化度（Ａｓｃ）は、血管を単位長さだけ変形させるために、どの程度力を与えるかを示す指標、即ち血管の硬化度を示す指標であり、血管の気質的変化を反映するものであって単位はＫｇ／ｃｍ２であり、一般的に弾性波伝播速度の２乗に比例する。
【００６８】
最後に、血流速度（Ｖ）は、脳血管枝に沿って流れる血液の速度であって単位はｃｍ／ｓであり、これは主に超音波センサ１２８により測定される。一方、脈波伝達速度（Plse Wave Velocity：ＰＷＶ）は、頚動脈と大腿動脈で脈波記録方法により測定したものであり、大動脈の弾力状態を反映する。血管壁が硬くなるほど速くなるが、特に動脈硬化性変化が激しいほど血流速度又は脈波伝達速度が速くなる。
【００６９】
また、各生物力学的指標を示す文字において、下付文字ａは前大脳動脈４１、ｂは脳底動脈４６、ｃは内頚動脈４２、ｄは心臓の弛緩期、ｍは中大脳動脈４３、ｐは後大脳動脈４５、ｓは心臓の収縮期、ｖは椎骨動脈４８、ａｃは前大脳交通動脈４０、ｐｃは後大脳交通動脈４４、１は左側、２は右側をそれぞれ指す。
【００７０】
一方、主処理部２１０には、使用者の制御命令の入力を受ける入力部２２０と、主処理部で算出された結果を表示する出力部２４０とがそれぞれ連結される。
【００７１】
ここで、出力部２４０は、プリンタだけでなく、モニタを通した画面出力部も含む。従って、図１に図示した映像処理部２３０は画面出力部に内蔵するようになる。
【００７２】
そして、入力部２２０は、通常のキーボード、マウスだけでなく、画面出力部（モニタ）に備えられたタッチ入力手段も含む。
【００７３】
上記のような構成で核心的な部分は、主処理部２１０の制御により生体信号を測定、分析してこれをもとに所定の数式により各生物力学的指標を算出させることにあるので、以下においては、これについて詳細に説明する。
【００７４】
主処理部２１０の制御は、基本的に入力部２２０により基礎情報（被検者の血圧、身長、体重、人種等）の入力を受け、生体信号測定計１００を通して生体信号の受信を受けた後（第１段階）、受信された生体信号の波形を分析し、分析された波形資料をもとに脳血管の圧力曲線、圧力曲線に関する収縮期と弛緩期の面積及び脳血管の血流量を得るステップ（第２段階）を経て、得た圧力曲線、圧力曲線の面積及び血流量と入力された基礎情報をもとに脳血管の順応性Ｃ及び抵抗Ｒを含む生物力学的指標を算出し、出力部２４０に脳血管の分析結果を表示するステップ（第３段階）が含まれるようにプログラムされたことを特徴として多様に実施することができる。
【００７５】
脳血管は、図４に示された脳血管枝を言うものであり、本実施例によると、後述の所定の数式により、左、右側の後大脳動脈４５、左、右側前大脳動脈４１及び左、右側中大脳動脈４３を含む各脳血管枝の生物力学的指標を自動的に計算し、出力部２４０に図８のようなＣ−Ｒ状態図（Ｃｈａｒｔ）等でその結果を表示するようになる。
【００７６】
後述のように、主処理部２１０により得られた測定データら（Ｑｃ１、Ｑｃ２、Ｑｖ１、Ｑｖ２）から複数個の方程式を連立して解くようにすることで、各脳血管枝の圧力曲線（Ｐａ１、Ｐａ２、Ｐｃ１、Ｐｃ２、Ｐｐ１、Ｐｐ２、Ｐｖ１、Ｐｖ２）を得て、各脳血管枝の生物力学的指標を計算するようになる。
【００７７】
その例を検討すると、左、右側の後大脳動脈の順応性（Ｃｐ１、Ｃｐ２）及び左、右側の後大脳動脈の抵抗（Ｒｐ１、Ｒｐ２）は、下記数式１〜４によりそれぞれ算出する。
【００７８】
左側後大脳動脈の順応性
【００７９】
【数１】
0005628180.tif 000002
【００８０】
右側後大脳動脈の順応性
【００８１】
【数２】
0005628180.tif 000003
【００８２】
左側後大脳動脈の抵抗
【００８３】
【数３】
0005628180.tif 000004
【００８４】
及び
右側後大脳動脈の抵抗
【００８５】
【数４】
0005628180.tif 000005
【００８６】
であり
数式１乃至４において、Ｐｐ１ｓは左側後大脳動脈の収縮期血圧、Ｐｐ１ｄは左側後大脳動脈の弛緩期血圧、Ｐｐ２ｓは右側後大脳動脈の収縮期血圧、Ｐｐ２ｄは右側後大脳動脈の弛緩期血圧、Ａｐ１ｓは左側後大脳動脈圧力曲線の収縮期面積、Ａｐ１ｄは左側後大脳動脈圧力曲線の弛緩期面積、Ａｐ２ｓは右側後大脳動脈圧力曲線の収縮期面積、Ａｐ２ｄは右側後大脳動脈圧力曲線の弛緩期面積、Ｓｐ１は左側後大脳動脈の血流量、Ｓｐ２は右側後大脳動脈の血流量である。
【００８７】
また、左、右側前大脳動脈の順応性（Ｃａ１、Ｃａ２）及び左、右側前大脳動脈の抵抗（Ｒａ１、Ｒａ２）は、下記数式５〜８によりそれぞれ算出される。
【００８８】
左側前大脳動脈の順応性
【００８９】
【数５】
0005628180.tif 000006
【００９０】
右側前大脳動脈の順応性
【００９１】
【数６】
0005628180.tif 000007
【００９２】
左側前大脳動脈の抵抗
【００９３】
【数７】
0005628180.tif 000008
【００９４】
及び
右側前大脳動脈の抵抗
【００９５】
【数８】
0005628180.tif 000009
【００９６】
であり
数式５乃至８において、Ｐａ１ｓは左側前大脳動脈の収縮期血圧、Ｐａ１ｄは左側前大脳動脈の弛緩期血圧、Ｐａ２ｓは右側前大脳動脈の収縮期血圧、Ｐａ２ｄは右側前大脳動脈の弛緩期血圧、Ａａ１ｓは左側前大脳動脈圧力曲線の収縮期面積、Ａａ１ｄは左側前大脳動脈圧力曲線の弛緩期面積、Ａａ２ｓは右側前大脳動脈圧力曲線の収縮期面積、Ａａ２ｄは右側前大脳動脈圧力曲線の弛緩期面積、Ｓａ１は左側前大脳動脈の血流量、Ｓａ２は右側前大脳動脈の血流量、Ｋは臨床係数である。
【００９７】
そして、左、右側中大脳動脈の順応性（Ｃｍ１、Ｃｍ２）及び左、右側中大脳動脈の抵抗（Ｒｍ１、Ｒｍ２）は、下記数式９〜１２によりそれぞれ算出される。
【００９８】
左側中大脳動脈の順応性
【００９９】
【数９】
0005628180.tif 000010
【０１００】
右側中大脳動脈の順応性
【０１０１】
【数１０】
0005628180.tif 000011
【０１０２】
左側中大脳動脈の抵抗
【０１０３】
【数１１】
0005628180.tif 000012
【０１０４】
及び
右側中大脳動脈の抵抗
【０１０５】
【数１２】
0005628180.tif 000013
【０１０６】
であり
数式９乃至１２において、Ｐｍ１ｓは左側中大脳動脈の収縮期血圧、Ｐｍ１ｄは左側中大脳動脈の弛緩期血圧、Ｐｍ１ｓは右側中大脳動脈の収縮期血圧、Ｐｍ２ｄは右側中大脳動脈の弛緩期血圧、Ａｍ１ｓは左側中大脳動脈圧力曲線の収縮期面積、Ａｍ１ｄは左側中大脳動脈圧力曲線の弛緩期面積、Ａｍ２ｓは右側中大脳動脈圧力曲線の収縮期面積、Ａｍ２ｄは右側中大脳動脈圧力曲線の弛緩期面積、Ｓｍ１は左側中大脳動脈の血流量、Ｓｍ２は右側中大脳動脈の血流量、Ｋは臨床係数である。
【０１０７】
一方、主処理部２１０は、第３段階で算出された脳血管の各枝の順応性Ｃ及び抵抗Ｒを出力部２４０にＣ−Ｒ状態図（Ｃｈａｒｔ）上の点で表示することができる。
【０１０８】
図８において、Ｃ−Ｒ状態図（Ｃｈａｒｔ）の領域区画は、多様な臨床結果によって精密度を高めることができるように、再区画することができるのは当然であるが、臨床結果の一例として各領域を定義すると、下記のとおりである。
【０１０９】
領域１は重症度脳血管の起始部攣縮であると判断される区域であり、領域２は血管狭窄の移行が始まる区域であり、領域３は血管狭窄が移行される区域であり、領域４は脳動脈硬化及び狭窄区域であり、領域５は起始部攣縮が疑われる区域であり、領域６は正常若しくは攣縮移行区域であり、領域７は正常区域であり、領域８と領域１０は正常若しくは自覚症状によって脳血管痙攣であると診断することができる区域であり、領域９は起始部攣縮であると判断される区域であり、領域１１は脳血管痙攣又は脳血管拡張剤の服用状態が疑われる区域である。
【０１１０】
以下においては、実施例を裏付ける関連理論及び臨床資料を補充的に提示する。
【０１１１】
本発明おいては、主処理部２１０が測定されたデータを利用して複数個の血圧と血流量に関する方程式を連立して解くように演算を進めることになる。
【０１１２】
主処理部２１０が脳血管状態を評価するための演算をするには、脳血管系が複雑であるために、先ず脳血管系を簡略化しなければならない。
【０１１３】
ここで、脳血管系の分析のために、進めた実験資料を紹介すると、表１のとおりである。
【０１１４】
【表１】
0005628180.tif 000014
表１と図４から分かるように、脳血管系は、内頚動脈、椎骨動脈、中大脳動脈、前大脳動脈、後大脳動脈を基本系として構成することができる（表１から分かるように、順応性は他の動脈らに比べて１／１００〜１／１００００程度であるので無視する）。
【０１１５】
図４と表１の実験結果から、脳血管系は、内頚動脈枝、前大脳動脈枝、中大脳動脈枝、後大脳動脈枝、椎骨動脈枝及び脳底動脈枝からなっているので、動脈枝ごとにウインドケッセル（Windkesell）模型を対応させ、血流特性を考慮して弾性鋼らに連結して多弾性鋼とし、弾性鋼ごとを１つの血流管路として分析すると、図４で表示された脳血管系を完璧に分析したものとなる。
【０１１６】
先ず、弾性鋼ごとを連結させて多弾性鋼とする場合をみると、血流特性を考慮して２つに模型化することができる。
【０１１７】
動脈枝ごとにウインドケッセルの模型を対応させ、血流特性を考慮して弾性鋼らを連結して多弾性鋼とする方法は、図５のように、内頚動脈で前大脳動脈と中大脳動脈が分岐する場合と仮定して多弾性鋼を構成する場合と、内頚動脈と中大脳動脈を１つの血管枝と仮定し、即ち内頚動脈＋中大脳動脈を１つの弾性鋼と仮定して多弾性鋼を構成する場合とに分けることができる。
【０１１８】
この場合に準じ、図４の脳血管系は、内頚動脈４２、椎骨動脈４８、中大脳動脈４３、前大脳動脈４１、後大脳動脈４５をそれぞれ１つの血流管路と仮定して考察することができる。
【０１１９】
図５で換算された脳血管系の血管枝らで血管の順応性と抵抗、弾性係数を求めるには、次のような問題を解明しなければならない。
【０１２０】
先ず、左右椎骨動脈４６と内頚動脈４２の入口で超音波により描かれた曲線と脈波曲線８個を与え、前大脳、中大脳、後大脳の血管枝の左右側の順応性と抵抗、弾性係数等１８個を求めることができないため、医学的に妥当でありながらも、生物力学では矛盾のない仮定を与え、図４を簡略化しなければならない。
【０１２１】
この問題を解明するために、表１の実験資料を検討すると、実験資料から後大脳動脈の順応性Ｃｐと後交通動脈の順応性Ｃｐｃの比は２５：１、前大脳動脈の順応性Ｃａと前大脳交通動脈の順応性Ｃａｃの比は４０：１である。
【０１２２】
また、脳血管に関する実験資料で得られたＣ、ＲとＰ１＝Ｐ３は、平均血圧１０３ｍｍＨｇ、Ｐ２＝Ｐ４の平均血圧を１０５ｍｍＨｇと与えて計算を進めると、Ｑｐｃ、Ｑａｃの血流量はＱａ、Ｑｐ、Ｑｍの１／３００程度しかない。
【０１２３】
以上の事実から前大脳交通動脈と後大脳交通動脈を力学的な固定端と見ることができる。
【０１２４】
従って、交通動脈を無視して動脈枝らの順応性と抵抗を求めることができる。
【０１２５】
このような仮定は、医学的にも妥当である。事実上前大脳交通動脈と後大脳交通動脈においては、脳血栓はそれほど起こらず、脳出血がたまに起こるので血圧と血流量のみ求め、順応性と抵抗を求めなくても診断においては不便がない。
【０１２６】
これにより、脳血管系は交通動脈で分割されたものと考察することができる。
【０１２７】
先ず、後大脳動脈を分析する問題について検討する。
【０１２８】
この問題は、脳底動脈で左右後大脳動脈が分岐される点でＰｐ１＝Ｐｐ２と左右Ｑｐ１、Ｑｐ２を与え、後大脳動脈をそれぞれ１つの弾性鋼と仮定して左右後大脳動脈の順応性と抵抗を求める問題と、後大脳動脈を単純管路と仮定して血流を考察する問題に帰する。
【０１２９】
後大脳動脈で順応性、血流抵抗、血流量を求めるのに最も難問題は、左右後大脳動脈が分岐される点でＰｐ１＝Ｐｐ２であるとき、左右後大脳動脈に流れる血流量（Ｑｐ１、Ｑｐ２）を求める問題である。椎骨動脈は、脳底動脈に合流したのち後大脳動脈で分岐される。
【０１３０】
脳底動脈で左右後大脳動脈に流れる血流量の特性を解明するために実験した結果、左右側Ｋｖ＊椎骨動脈＋Ｋｃ＊内頚動脈の血流量比と左右側後大脳動脈の血流量比との間には、高い相関関係があることを示した。
【０１３１】
その実験過程は、次のとおりである。
【０１３２】
脳底動脈で左右後大脳動脈に流れる血流量の特性を鶏鳴するために、総５０名の男たちを標本として実験を進めた。
【０１３３】
平均動脈圧は１１８ｍｍＨｇ〜１３２ｍｍＨｇ、脳底動脈の入口における血流量は５．２〜７．８ｍＬ／ｓ、総心動周期は５４１個を取った。実験結果は、下記表２のとおりである。
【０１３４】
【表２】
0005628180.tif 000015
【０１３５】
上記実験結果、左右側椎骨動脈＋内頚動脈の血流量比（ＳＶ１＊／ＳＶ２＊）と左右側後大脳動脈の血流量比（ＳＰ１／ＳＰ２）との間には、下記数式１３、１４のように表現される高い相関関係があることが分かった。
【０１３６】
【数１３】
0005628180.tif 000016
【０１３７】
【数１４】
0005628180.tif 000017
【０１３８】
上記実験結果を理論的に考察してみると、左右内頚動脈＋椎骨動脈の血流量が左右後大脳血管の血流量と相関がある理由は、次のとおり説明することができる。
【０１３９】
脳底動脈が後大脳動脈と連結される部位においては、シグマ効果が生じて前下小脳動脈と後下小動脈、上小脳血管は、脳底動脈で血液を吸い込む‘吸入点’のように挙動する。また、内頚動脈で流れる血液流れは、後大脳交通動脈を通して後大脳動脈の血液流れを補充する。その量は後大脳に流れる総血流量の３０〜３８％まである。一方、脳に流入される血液は、流体力学的速度境界層‘引き’により軸対称の大田幹を保存しようとする特徴を有している。
【０１４０】
このような現象を考慮し、椎骨動脈から脳底動脈を通して流れる血液の流れ量は、左右後大脳動脈に必ずしも同様に流れない。
【０１４１】
η＝Ｓ＊Ｖ１／Ｓ＊Ｖ２を左右側Ｋｖ＊椎骨動脈の血流量＋Ｋｃ＊頚動脈の血流量の比とし（ここで、ＫｖとＫｃはそれぞれ実験定数であって、Ｋｖ＝０．１３１〜０．１５２、Ｋｃ＝０．７３〜０．８２である）、ζ＝ＳＰ１／ＳＰ２を左、右側後大脳動脈の血流量の比とすると、数式１３は下記数式１５のように表現される。
【０１４２】
【数１５】
0005628180.tif 000018
【０１４３】
従って、後大脳動脈に流れる血流量は、椎骨動脈の血流量の８０〜８５％と内頚動脈で２９％〜３２％補充されるという事情を考慮し、下記数式１６のように求めることができる。
【０１４４】
【数１６】
0005628180.tif 000019
【０１４５】
Ｓｐ２＝１．２４（Ｑｖ１＋Ｑｖ２）−Ｓｐ１  （数式１７）
次に、左右後大脳動脈が分岐する点でＰｐ１＝Ｐｐ２を求める。
Ｐｐ１＝Ｐｐ２＝Ｐ３−Ｒｖ３＊Ｑ３−Ｒｂ＊Ｑｂ  （数式１８）
数式１８において、Ｒｖ３、ＲｂはPoisenilleの公式により、
【０１４６】
【数１７】
0005628180.tif 000020
のように求める。
【０１４７】
のように求める。
【０１４８】
上記式において、λは動脈の長さ、Ｄは直径、μは血液の粘稠度である。
【０１４９】
Ｄは流体力学で広く知られているプランク公式により求める。
【０１５０】
Ｑ、Ｐを求めた後、主処理部２１０は、血管の順応性Ｃｐｓ、ＣｐｄとＲｐ、Ｒｂを求める。ここで、後大脳動脈は、脳底動脈から左右に分岐されるので、後大脳交通動脈を固定端とする単弾性鋼として取り扱うことができる（図５参照）。
【０１５１】
一方、脳血管は、攣縮作用と痙攣が激しいので、後大脳動脈の模型化方程式は、血管が収縮するときと血管が拡張するときとに分けて考察する。
【０１５２】
【数１８】
0005628180.tif 000021
【０１５３】
【数１９】
0005628180.tif 000022
【０１５４】
実験資料によると、１７０〜１８０ｍｍＨｇ血圧で収縮期血管の順応性と拡張期血管の順応性とは同一である。ところが、後大脳動脈内圧は普通３４〜１１０ｍｍＨｇ程度であるので、後大脳動脈で血管の収縮時の順応性と拡張時の順応性との大きさは同一である。
従って、Ｃｐｓ＝Ｃｐｄ＝Ｃである。
【０１５５】
数式１９と２０は、Ｐ、Ｑ、Ｃ、Ｒの間の関係である。Ｃ、Ｒを求める算法は、血圧曲線Ｐと血流量曲線Ｑとが一致するまでＣ、Ｒを調整する代わりに、血圧曲線Ｐの面積と血流量Ｓとの間の関数関係によりＣ、Ｒを求める。
【０１５６】
面積対面積の間の関数関係を得ると、再現性のあるＣ、Ｒを得ることができる。
【０１５７】
このような目的から、数式１９と２０をそれぞれ積分して足し引きすると、下記数式２１のようになる。
【０１５８】
【数２０】
0005628180.tif 000023
【０１５９】
ここで、Ｓｖ＊は後大脳動脈に流れる１回拍出時の血流量であり、Ｐｓは収縮期血圧であり、Ｐｄは弛緩期（拡張期）血圧であり、Ａｓは収縮期の間の血圧曲線Ｐの面積であり、Ａｄは弛緩期の間の血圧曲線Ｐの面積である。
【０１６０】
従って、数式２１から左、右側後大脳動脈の順応性（Ｃｐ１、Ｃｐ２）及び抵抗（Ｒｐ１、Ｒｐ２）を求めると、下記のとおりである。
【０１６１】
左側後大脳動脈の順応性
【０１６２】
【数２１】
0005628180.tif 000024
【０１６３】
右側後大動脈の順応性
【０１６４】
【数２２】
0005628180.tif 000025
【０１６５】
左側後大脳動脈の抵抗
【０１６６】
【数２３】
0005628180.tif 000026
【０１６７】
及び
右側後大脳動脈の抵抗
【０１６８】
【数２４】
0005628180.tif 000027
【０１６９】
数式１乃至４において、Ｐｐ１ｓは左側後大脳動脈の収縮期血圧、Ｐｐ１ｄは左側後大脳動脈の弛緩期血圧、Ｐｐ２ｓは右側後大脳動脈の収縮期血圧、Ｐｐ２ｄは右側後大脳動脈の弛緩期血圧、Ａｐ１ｓは左側後大脳動脈圧力曲線の収縮期の面積、Ａｐ１ｄは左側後大脳動脈圧力曲線の弛緩期の面積、Ａｐ２ｓは右側後大脳動脈圧力曲線の収縮期の面積、Ａｐ２ｄは右側後大脳動脈圧力曲線の弛緩期の面積、Ｓｐ１は左側後大脳動脈の血流量、Ｓｐ２は右側後大脳動脈の血流量である。
【０１７０】
次に、後大脳動脈を弾性管に血液が流れる単純管路と仮定し、血液が流れる弾性管路で有弾性体の問題を解き、脳血管の気質的変化と機能的変化を区分する問題を解決する（図５参照）。
【０１７１】
１つの血管の弾性管で有弾性体の連続方程式と運動方程式を考察すると、
【０１７２】
【数２５】
0005628180.tif 000028
【０１７３】
【数２６】
0005628180.tif 000029
【０１７４】
であり、式２２、２３において、Ａは血管の断面積、Ｕは血流速度、Ｐは血圧である。
【０１７５】
【数２７】
0005628180.tif 000030
【０１７６】
数式２４において、Ｙは血管半径、μは粘度、τｗは接線応力である。
【０１７７】
【数２８】
0005628180.tif 000031
【０１７８】
数式２５において、Ｆは平均血流速度、ａは弾性波伝播速度である。
【０１７９】
【数２９】
0005628180.tif 000032
【０１８０】
数式２２乃至２６を考慮してまとめると、
【０１８１】
【数３０】
0005628180.tif 000033
【０１８２】
【数３１】
0005628180.tif 000034
【０１８３】
数式２７、２８において、Ｐは血圧曲線、μは粘度、Ａは血管の断面積、ρは血液密度である。
【０１８４】
ここで、数式２８において、
【０１８５】
【数３２】
0005628180.tif 000035
【０１８６】
を無視してＸによって積分すると、
【０１８７】
【数３３】
0005628180.tif 000036
【０１８８】
これにより、単弾性管において、
【０１８９】
【数３４】
0005628180.tif 000037
【０１９０】
が成立する。
【０１９１】
数式３０から分かるように、脳血管系で起こる血圧変動、脳血管の痙攣発作、脳血管の攣縮、薬物作用等による補充的内圧の作用は、脳血管の断面積の変化をもたらす。血管の順応性と抵抗から分かるように、順応性と抵抗が血圧変動、痙攣発作、攣縮、薬物作用等により激しい変動をもたらすのは、脳血管系の断面積の変化のためである。
【０１９２】
ところが、弾性係数Ｅは、弾性波伝播速度で表示されるので、脳血管で血圧変動、痙攣発作、攣縮、薬物作用等に関係なく、脳血管の気質的変化を代表することができる。
【０１９３】
一方、Moensu Kortewegによると、
【０１９４】
【数３５】
0005628180.tif 000038
【０１９５】
であるので、弾性係数Ｅ＝ρ（ｄ／ｈ）ＰＷＶ２である。
【０１９６】
従って、Ｃ、ＲからＡを消去し、μとＰＷＶの相関式から脳血管の動脈硬化度Ａｓｃを得る。
【０１９７】
【数３６】
0005628180.tif 000039
【０１９８】
数式３１において、Ｓ＝ｆ（ＰＷＶ）、Ｋ３は臨床で得られた係数である。
【０１９９】
次に、上記で提起した脳血管特性と血流特性を反映する指標を臨床に適用するには、内頚動脈及び椎骨動脈に流れる血流量を求める問題を解明しなければならない。
【０２００】
本発明においては、内頚動脈波曲線と超音波ドップラーにより求めた血流量曲線とが一致するまでＣ、Ｒを校正してＣ、Ｒを求める代わりに、超音波ドップラーにより得た頚椎動脈の血流量を内頚動脈波曲線の面積で表示されたＣ、Ｒ公式に代入する方式によりＣ、Ｒを求める。
【０２０１】
この方法は、脳血流特性を測定するときに経験する超音波ドップラーの弱点である水平面誤差がＣ、Ｒを求めるのに何の支障も与えないということを言う。
【０２０２】
言い換えれば、現在超音波測定技術により測定した血流量が臨床で無難に用いられているので、超音波測定技術により測定した血流量のみ求めると、Ｃ、Ｒ、Ａｓｃも臨床に無難に用いることができるということを示す。
【０２０３】
次に、交通動脈を無視した内頚動脈を分析する問題について検討する。
【０２０４】
この場合も上記の場合のように、前大脳動脈と中大脳動脈とが分岐する点で血圧Ｐａ＝Ｐｍ＝Ｐ、この点で血流量をそれぞれＱａ、Ｑｍとし、前大脳動脈と中大脳動脈をそれぞれ１つの弾性鋼と仮定し、次のような弾性鋼方程式を立てることができる。
【０２０５】
血管収縮時と血管拡張時との順応性は同一であるので、
【０２０６】
【数３７】
0005628180.tif 000040
【０２０７】
【数３８】
0005628180.tif 000041
【０２０８】
【数３９】
0005628180.tif 000042
【０２０９】
である。
【０２１０】
上記右側血管抵抗Ｒ２は５．９×１０４ｄｙｎ・Ｓ／ｃｍ５程度であり、左側血管抵抗Ｒ１は３４００ｄｙｎ・Ｓ／ｃｍ５であるので、Ｒ１／Ｒ２＝０である。
【０２１１】
従って、
【０２１２】
【数４０】
0005628180.tif 000043
【０２１３】
【数４１】
0005628180.tif 000044
【０２１４】
【数４２】
0005628180.tif 000045
【０２１５】
上記数式３１、３５、３６のＲ、Ｃ、Ａｓｃは、血流量が分からないため診断に用いることができない。
【０２１６】
従って、診断に用いることができるＲ、Ｃ、Ａｓｃを求めるために、Ｐｍ＝Ｐａの条件でＱｍとＱａを求めなければならない。
【０２１７】
この問題を解くために、内頚動脈が中大脳動脈と１つに連結されており、頚動脈＋中大脳動脈枝で前大脳動脈が分岐されたものと仮定し、頚動脈系を再び模型化する（図６参照）。
【０２１８】
ここで、中大脳動脈の血管の順応性をＣｍ、前大脳の血管の順応性をＣａ、中大脳動脈における血流抵抗をＲｍ、前大脳動脈のウィリス輪部門における血流抵抗をＲａ、前大脳動脈の残りの部門における血流抵抗をＲｔ、中大脳における血圧をＰｍ、前大脳における血圧をＰａ、中大脳動脈における血流量をＱｍ、前大脳動脈における血流量をＱａ、内頚動脈でウィリス輪の連結点における血圧をＰ、静脈圧はＰｖとしよう。
【０２１９】
そして、前大脳動脈と中大脳動脈の総順応性をＣとすると、
【０２２０】
【数４３】
0005628180.tif 000046
【０２２１】
である。
【０２２２】
ここで、図６から内頚動脈ＰＷＶを１次近似で中大脳動脈のＰＷＶとすると、内頚動脈＋中大脳動脈の順応性Ｃｃｍを求めることができる。
従って、
【０２２３】
【数４４】
0005628180.tif 000047
【０２２４】
【数４５】
0005628180.tif 000048
【０２２５】
【数４６】
0005628180.tif 000049
【０２２６】
以上の結果を総合すると、左、右側前大脳動脈の順応性（Ｃａ１、Ｃａ２）及び左、右側前大脳動脈の抵抗（Ｒａ１、Ｒａ２）は、下記数式５〜８によりそれぞれ算出することがでるようになる。
【０２２７】
左側前大脳動脈の順応性
【０２２８】
【数４７】
0005628180.tif 000050
【０２２９】
右側前大脳動脈の順応性
【０２３０】
【数４８】
0005628180.tif 000051
【０２３１】
左側前大脳動脈の抵抗
【０２３２】
【数４９】
0005628180.tif 000052
【０２３３】
及び
右側前大脳動脈の抵抗
【０２３４】
【数５０】
0005628180.tif 000053
【０２３５】
上記数式５乃至８において、Ｐａ１ｓは左側前大脳動脈の収縮期血圧、Ｐａ１ｄは左側前大脳動脈の弛緩期血圧、Ｐａ２ｓは右側前大脳動脈の収縮期血圧、Ｐａ２ｄは右側前大脳動脈の弛緩期血圧、Ａａ１ｓは左側前大脳動脈圧力曲線の収縮期面積、Ａａ１ｄは左側前大脳動脈圧力曲線の弛緩期面積、Ａａ２ｓは右側前大脳動脈圧力曲線の収縮期面積、Ａａ２ｄは右側前大脳動脈圧力曲線の弛緩期面積、Ｓａ１は左側前大脳動脈の血流量、Ｓａ２は右側前大脳動脈の血流量、Ｋは臨床係数である。
【０２３６】
また、左、右側中大脳動脈の順応性（Ｃｍ１、Ｃｍ２）及び左、右側中大脳動脈の抵抗（Ｒｍ１、Ｒｍ２）は、下記数式９〜１２によりそれぞれ算出することができるようになる。
【０２３７】
左側中大脳動脈の順応性
【０２３８】
【数５１】
0005628180.tif 000054
【０２３９】
右側中大脳動脈の順応性
【０２４０】
【数５２】
0005628180.tif 000055
【０２４１】
左側中大脳動脈の抵抗
【０２４２】
【数５３】
0005628180.tif 000056
【０２４３】
及び
右側中大脳動脈の抵抗
【０２４４】
【数５４】
0005628180.tif 000057
【０２４５】
上記数式９乃至１２において、Ｐｍ１ｓは左側中大脳動脈の収縮期血圧、Ｐｍ１ｄは左側中大脳動脈の弛緩期血圧、Ｐｍ２ｓは右側中大脳動脈の収縮期血圧、Ｐｍ２ｄは右側中大脳動脈の弛緩期血圧、Ａｍ１ｓは左側中大脳動脈圧力曲線の収縮期面積、Ａｍ１ｄは左側中大脳動脈圧力曲線の弛緩期面積、Ａｍ２ｓは右側中大脳動脈圧力曲線の収縮期面積、Ａｍ２ｄは右側中大脳動脈圧力曲線の弛緩期面積、Ｓｍ１は左側中大脳動脈の血流量、Ｓｍ２は右側中大脳動脈の血流量、Ｋは臨床係数である。
【０２４６】
次に、頚動脈と椎骨動脈において、血流量と頚動脈で平均血圧Ｐａｖが与えられた条件で脳血管系の血管枝ごとで血圧と血流量を求める。
【０２４７】
このとき、血管ごとの抵抗Ｒは、上記で得た値を与える。
【０２４８】
臨床において、血圧は、抵抗と血流量の積であると定義する。
【０２４９】
従って、上記内容をもとに本発明による処理部２１０で血圧と血流量を求めるための連立方程式を作成すると、
【０２５０】
【数５５】
0005628180.tif 000058
【０２５１】
【数５６】
0005628180.tif 000059
【０２５２】
【数５７】
0005628180.tif 000060
【０２５３】
【数５８】
0005628180.tif 000061
【０２５４】
【数５９】
0005628180.tif 000062
【０２５５】
【数６０】
0005628180.tif 000063
【０２５６】
【数６１】
0005628180.tif 000064
【０２５７】
【数６２】
0005628180.tif 000065
【０２５８】
【数６３】
0005628180.tif 000066
【０２５９】
【数６４】
0005628180.tif 000067
【０２６０】
【数６５】
0005628180.tif 000068
【０２６１】
【数６６】
0005628180.tif 000069
【０２６２】
【数６７】
0005628180.tif 000070
【０２６３】
【数６８】
0005628180.tif 000071
【０２６４】
【数６９】
0005628180.tif 000072
【０２６５】
【数７０】
0005628180.tif 000073
【０２６６】
【数７１】
0005628180.tif 000074
【０２６７】
【数７２】
0005628180.tif 000075
【０２６８】
【数７３】
0005628180.tif 000076
【０２６９】
【数７４】
0005628180.tif 000077
【０２７０】
【数７５】
0005628180.tif 000078
【０２７１】
【数７６】
0005628180.tif 000079
【０２７２】
【数７７】
0005628180.tif 000080
【０２７３】
である。
【０２７４】
以上の方程式において、未知数は、Ｑｐ１１、Ｑｐ１２、Ｑｐ２１、Ｑｐ２２、Ｑａ１１、Ｑａ１２、Ｑａ２１、Ｑａ２２、Ｑｍ１、Ｑｍ２、Ｑｐｃ１、Ｑｐｃ２、Ｑａｃ、Ｐａ１、Ｐａ２、Ｐｃ１、Ｐｃ２、Ｐｖ１、Ｐｖ２、Ｐｐ１、Ｐｐ２、Ｐ３、Ｐ４、Ｒｂであり、
Ｐ１、Ｐ２、Ｑｖ１、Ｑｖ２、Ｑｃ１、Ｑｃ２が分かるとすると、上記連立方程式は閉ざされるようになり、未知数を求めることができるようになる。
【産業上の利用可能性】
【０２７５】
本発明の脳血管分析装置によると、脳血管の個別的血管枝の気質的変化を反映する弾性係数を分析し、脳血管系の血流特性と脳血管の個別的血管枝の気質的変化と機能的変化を反映する血流量、血管の順応性、血流抵抗を算出することにより、脳血管疾患の発生危険性を早期に診断することができるものであるので、産業上利用可能性が非常に高い。
【符号の説明】
【０２７６】
１０：カフ血圧計
１１：カフ
１２：接着手段（ベルクロ）
１３：内蔵された空気袋
１４、１７、１８：ゴム管
１５：空気排出バルブ
１６：送気口
２０：アダプタ
２１：支路管
２２：支路管装着部
２４：蓋
２６：アダプタ突出連結部
３０：脈波センサ（圧力センサ）
３１：通気孔
３２：開放溝
３４：ハウジング胴体
３６：センシングリード線
４０：前大脳交通動脈
４１：前大脳動脈
４２：内頚動脈
４３：中大脳動脈
４４：後大脳交通動脈
４５：後大脳動脈
４６：脳底動脈
４７：前下小脳動脈
４８：椎骨動脈
４９：後下小脳動脈
【図１】
0005628180.tif 000093
【図２】
0005628180.tif 000094
【図３】
0005628180.tif 000095
【図４】
0005628180.tif 000096
【図５】
0005628180.tif 000097
【図６】
0005628180.tif 000098
【図７】
0005628180.tif 000099
【図８】
0005628180.tif 000100