(57)【要約】【課題】  簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出することができる圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法を提供する。【解決手段】  生体動脈の脈波データと、脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出装置１であって、拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、拍間隔と脈波データとの相関係数が、第ｎ−１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に相関係数が閾値以下となる系列を圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列として検出するＶＢＲＳＡ系列検出部４４を有する。
2018181851.tif 000002
【特許請求の範囲】
【請求項１】
生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出装置であって、
前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ−１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出部を有する、圧反射性血管交感神経活動検出装置。
【請求項２】
前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ−１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出部と、
第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出部と
を有する、請求項１に記載の圧反射性血管交感神経活動検出装置。
【請求項３】
前記ＶＢＲＳＡ値算出部は、下記式（１）〜（３）のいずれかを用いて前記ＶＢＲＳＡ値を算出する、請求項２に記載の圧反射性血管交感神経活動検出装置；
ＶＢＲＳＡ値＝｜ＰＶｔ−ＰＶｍ｜  ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）−ｌｎ（ＰＶｍ）｜  ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜  ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値
【請求項４】
前記ＶＢＲＳＡ値算出部は、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出する、請求項２または請求項３に記載の圧反射性血管交感神経活動検出装置。
【請求項５】
前記ＶＢＲＳＡ値算出部は、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出する、請求項４に記載の圧反射性血管交感神経活動検出装置。
【請求項６】
前記脈波データは、光電式容積脈波、規準化された前記光電式容積脈波である規準化脈波容積、または対数化された前記規準化脈波容積のうちのいずれかである、請求項１から請求項５のいずれかに記載の圧反射性血管交感神経活動検出装置。
【請求項７】
生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出プログラムであって、
前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ−１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出部としてコンピュータを機能させる、圧反射性血管交感神経活動測定プログラム。
【請求項８】
前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ−１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出部と、
第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出部としてコンピュータを機能させる、請求項７に記載の圧反射性血管交感神経活動検出プログラム。
【請求項９】
前記ＶＢＲＳＡ値算出部は、下記式（１）〜（３）のいずれかを用いて前記ＶＢＲＳＡ値を算出する、請求項８に記載の圧反射性血管交感神経活動検出プログラム；
ＶＢＲＳＡ値＝｜ＰＶｔ−ＰＶｍ｜  ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）−ｌｎ（ＰＶｍ）｜  ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜  ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値
【請求項１０】
前記ＶＢＲＳＡ値算出部は、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出する、請求項８または請求項９に記載の圧反射性血管交感神経活動検出プログラム。
【請求項１１】
前記ＶＢＲＳＡ値算出部は、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出する、請求項１０に記載の圧反射性血管交感神経活動検出プログラム。
【請求項１２】
前記脈波データは、光電式容積脈波、規準化された前記光電式容積脈波である規準化脈波容積、または対数化された前記規準化脈波容積のうちのいずれかである、請求項７から請求項１１のいずれかに記載の圧反射性血管交感神経活動検出プログラム。
【請求項１３】
生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出方法であって、
前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ−１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出ステップを有する、圧反射性血管交感神経活動測定方法。
【請求項１４】
前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ−１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出ステップと、
第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出ステップとを有する、請求項１３に記載の圧反射性血管交感神経活動検出方法。
【請求項１５】
前記ＶＢＲＳＡ値算出ステップでは、下記式（１）〜（３）のいずれかを用いて前記ＶＢＲＳＡ値が算出される、請求項１４に記載の圧反射性血管交感神経活動検出方法；
ＶＢＲＳＡ値＝｜ＰＶｔ−ＰＶｍ｜  ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）−ｌｎ（ＰＶｍ）｜  ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜  ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値
【請求項１６】
前記ＶＢＲＳＡ値算出ステップでは、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出する、請求項１４または請求項１５に記載の圧反射性血管交感神経活動検出方法。
【請求項１７】
前記ＶＢＲＳＡ値算出ステップでは、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出する、請求項１６に記載の圧反射性血管交感神経活動検出方法。
【請求項１８】
前記脈波データは、光電式容積脈波、規準化された前記光電式容積脈波である規準化脈波容積、または対数化された前記規準化脈波容積のうちのいずれかである、請求項１３から請求項１７のいずれかに記載の圧反射性血管交感神経活動検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血圧を一定に保つ役割を果たす圧反射機能に関わる血管の交感神経活動を非侵襲的に検出する圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法に関するものである。
【背景技術】
【０００２】
従来、血圧は、心血管疾患や脳梗塞等における重要なリスクファクターであることが知られている。このため、血圧の状態を日常的にモニターすることにより、上述した心血管疾患や脳梗塞等を予防するための指針や、プレクリニカル段階での診断材料等として利用することができる。しかしながら、血圧を測定する行為は、身体への負担が大きく面倒な作業でもあるため、現実的には日常的または継続的に血圧をモニタリングすることは困難である。【０００３】
そこで近年、血圧を一定の範囲に保とうとする自律神経系の圧反射機能を評価する研究が行われている。当該圧反射機能の低下が、脳卒中、高血圧症および失神等におけるリスクファクターと考えられているためである。ただし、血圧は、血行力学上、心拍出量（一回拍出量×心拍数）と血管抵抗との積として定義されるところ、圧反射においては、上記心拍数に対する心臓の交感神経活動の寄与は見られないため、上記心拍数は心臓の副交感神経活動によって制御され、上記血管抵抗は血管の交感神経活動によって制御されている。このため、前記圧反射機能には、心臓の副交感神経活動（心臓の神経性圧反射機能）と、血管の交感神経活動（血管の神経性圧反射機能）との双方が関わっている。【０００４】
具体的には、図１５に示すように、心臓の神経性圧反射機能では、血圧の上昇により血管壁が伸展すると、頸動脈等に存在する圧受容体細胞が発火する。そして、当該発火信号を受けた脳が心臓の副交感神経活動を介して心拍数を減少させて（すなわち、心拍間隔または脈拍間隔を延長させて）心拍出量を下げ、血圧を下げようとする。一方、血管の神経性圧反射機能は、血管の弛緩により血圧が低下すると、所定値以下にならないよう、血管の交感神経活動を賦活して血管を緊張させ、血圧を上げようとする。【０００５】
ただし、心臓の神経性圧反射機能と血管の神経性圧反射機能との作用バランスには個人差があり、心臓の神経性圧反射機能が優位に働く人もいれば、血管の神経性圧反射機能が優位に働く人もいる。また、同じ個人であっても、例えば安静にしている場合と運動した後とでは、上記作用バランスが変わることもある。したがって、圧反射機能を総合的に評価するには、心臓の神経性圧反射機能と血管の神経性圧反射機能とを別個独立に測定する必要がある。【０００６】
心臓の神経性圧反射機能を測定する技術としては、当該機能を感度として非侵襲的に測定する方法が存在する。例えば、血管壁の伸展による血管径の増加や減少を超音波エコーにより毎拍測定し、これとともに生じる心拍ないし脈拍間隔の延長ないし短縮として現れる心臓の圧反射系列に基づいて検出する、いわゆるシーケンス法が知られている。また、本願発明者は、血管径ではなく、規準化脈波容積（非特許文献１）による血液容積を用いて、心臓の神経性圧反射感度を測定する方法を提案している（特許文献１）。【０００７】
一方、血管の神経性圧反射機能を測定する技術としては、従来、図１６に示すように、拡張期血圧または血管抵抗と、いわゆるマイクロニューログラフィによって侵襲的に測定された血管の交感神経活動（MSA：Muscle nerve Sympathetic Activity）との関数から導き出す方法しか存在していない（非特許文献２）。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】国際公開第２０１４／１５７６０５号
【非特許文献】
【０００９】
【非特許文献１】Y. Sawadaら,「Normalized pulse volume (NPV) derived photo-plethysmographically as a more valid measure of the finger vascular tone, International Journal of Psychophysiology, 2001, Vol.41」, p1-10, Elsevier B.V.
【非特許文献２】G. Sundlofら,「Human muscle nerve sympathetic activity at rest. Relationship to blood pressure and age, The Journal of Physiology, 1978, Vol.274」, p621-637, The Physiological Society.
【発明の概要】
【発明が解決しようとする課題】
【００１０】
上述したマイクロニューログラフィは、血管の交感神経活動を直接測定できるものの、末梢神経束内に針のような金属の微小電極を刺入しなければならず、侵襲的であるという問題がある。また、マイクロニューログラフィを実施するには、病院検査のように行動が制限されたり（安静等）、整備された施設（実験室や検査室等）や熟練の測定者を必要とするため、測定可能な時間、場所、シチュエーション等が極めて限定的であるという問題もある。【００１１】
本発明は、このような問題点を解決するためになされたものであって、簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出することができる圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法を提供することを目的としている。
【課題を解決するための手段】
【００１２】
本発明に係る圧反射性血管交感神経活動検出装置は、簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出するという課題を解決するために、生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出装置であって、前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ−１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出部を有する。また、本発明に係る圧反射性血管交感神経活動検出プログラムは、前記ＶＢＲＳＡ系列検出部としてコンピュータを機能させるものである。【００１３】
また、本発明に係る圧反射性血管交感神経活動検出装置の一態様として、圧反射性血管交感神経活動をＶＢＲＳＡ値として指標化し、圧反射性血管交感神経活動を簡易的かつ客観的に評価するという課題を解決するために、前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ−１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出部と、第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出部とを有していてもよい。また、本発明に係る圧反射性血管交感神経活動検出プログラムの一態様として、前記理論値算出部と、前記ＶＢＲＳＡ値算出部としてコンピュータを機能させてもよい。【００１４】
さらに、本発明の一態様として、圧反射性血管交感神経活動を評価するのに適したＶＢＲＳＡ値を算出するという課題を解決するために、前記ＶＢＲＳＡ値算出部は、下記式（１）〜（３）のいずれかを用いて前記ＶＢＲＳＡ値を算出してもよい。
ＶＢＲＳＡ値＝｜ＰＶｔ−ＰＶｍ｜  ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）−ｌｎ（ＰＶｍ）｜  ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜  ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値【００１５】
また、本発明の一態様として、上昇系列と下降系列とのそれぞれについてＶＢＲＳＡ値を個別に算出することで、疾患に応じて異なる傾向等を把握し、評価の精度を向上するという課題を解決するために、前記ＶＢＲＳＡ値算出部は、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出してもよい。【００１６】
さらに、本発明の一態様として、上昇系列および下降系列の合算値をＶＢＲＳＡ値として算出することで、ＶＢＲＳＡ系列の検出数が少なくても算出し、圧反射性血管交感神経活動を示す指標を分かりやすく表現するという課題を解決するために、前記ＶＢＲＳＡ値算出部は、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出してもよい。【００１７】
また、本発明の一態様として、無次元の絶対量である規準化脈波容積（NPV）や、対数化された規準化脈波容積（lnNPV）を用いることにより、生体の組織成分に起因する誤差を除外し、同一人の測定結果のみならず、異なる個人間の測定結果も比較可能とし、広く一般的に圧反射性血管交感神経活動の動向を把握するという課題を解決するために、前記脈波データは、光電式容積脈波、規準化された前記光電式容積脈波である規準化脈波容積、または対数化された前記規準化脈波容積のうちのいずれかであってもよい。【００１８】
また、本発明に係る圧反射性血管交感神経活動検出方法は、簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出するという課題を解決するために、生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出方法であって、前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ−１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出ステップを有する。【００１９】
また、本発明に係る圧反射性血管交感神経活動検出方法の一態様として、圧反射性血管交感神経活動をＶＢＲＳＡ値として指標化し、圧反射性血管交感神経活動を簡易的かつ客観的に評価するという課題を解決するために、前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ−１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出ステップと、第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出ステップとを有していてもよい。【００２０】
さらに、本発明の一態様として、圧反射性血管交感神経活動を評価するのに適したＶＢＲＳＡ値を算出するという課題を解決するために、前記ＶＢＲＳＡ値算出ステップでは、下記式（１）〜（３）のいずれかを用いて前記ＶＢＲＳＡ値が算出されてもよい。
ＶＢＲＳＡ値＝｜ＰＶｔ−ＰＶｍ｜  ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）−ｌｎ（ＰＶｍ）｜  ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜  ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値【００２１】
また、本発明の一態様として、上昇系列と下降系列とのそれぞれについてＶＢＲＳＡ値を個別に算出することで、疾患に応じて異なる傾向等を把握し、評価の精度を向上するという課題を解決するために、前記ＶＢＲＳＡ値算出ステップでは、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出してもよい。【００２２】
さらに、本発明の一態様として、上昇系列および下降系列の合算値をＶＢＲＳＡ値として算出することで、ＶＢＲＳＡ系列の検出数が少なくても算出し、圧反射性血管交感神経活動を示す指標を分かりやすく表現するという課題を解決するために、前記ＶＢＲＳＡ値算出ステップでは、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出してもよい。
【発明の効果】
【００２３】
本発明によれば、簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出することができる。
【図面の簡単な説明】
【００２４】
【図１】本発明に係る圧反射性血管交感神経活動検出装置および圧反射性血管交感神経活動検出プログラムの一実施形態を示すブロック図である。
【図２】（ａ）生体組織（静脈部分は省略）の光電式容積脈波を測定する様子を示す断面図、および（ｂ）光電式容積脈波の交流成分（PGac）と直流成分（PGdc）を説明する図である。
【図３】本発明に係る「拍間隔」および「脈波データに対応する拍間隔」を説明する図である。
【図４】（ａ）血圧（BP）、心拍間隔（RRI）および規準化脈波容積（NPV）の実測データ、および（ｂ）当該実測データにおける収縮期血圧（SBP）と規準化脈波容積（NPV）との散布図である。
【図５】図４（ａ）の実測データにおける、（ａ）収縮期血圧（SBP）と心拍間隔（RRI）との散布図、および（ｂ）規準化脈波容積（NPV）と心拍間隔（RRI）との散布図である。
【図６】実測データにおける心臓の神経性圧反射感度（Cardiac nBRS）および圧反射性血管交感神経活動系列（VBRSA系列）についての検出時間および測定結果を示す図である。
【図７】（ａ）心臓の神経性圧反射感度（Cardiac nBRS）の測定方法、および（ｂ）圧反射性血管交感神経活動（Vascular BRSA）の算出方法を示す図である。
【図８】本発明に係る圧反射性血管交感神経活動検出方法の一実施形態を示すフローチャート図である。
【図９】実施例１における個人の実測データを示すグラフである。
【図１０】図９の実測データにおける、ＶＢＲＳＡ値（Vascular BRSA）と頸動脈の血管抵抗（Carotid arterial resistance）との散布図である。
【図１１】実施例２において、精神的ストレスに対する反応機序の異なる群ごとに測定された、（ａ）平均血圧（MBP）、（ｂ）頸動脈の血管抵抗（Carotid arterial resistance）および（ｃ）ＶＢＲＳＡ値（Vascular BRSA）を示すグラフである。
【図１２】実施例２における、全被験者の実測データについてのＶＢＲＳＡ値（Vascular BRSA）と頸動脈の血管抵抗（Carotid arterial resistance）との相関図である。
【図１３】実施例３における日常生活上の実測データを示すグラフである。
【図１４】図１３の実測データのうち、血圧測定期間におけるＶＢＲＳＡ値と拡張期血圧（DBP）との相関図である。
【図１５】圧反射機能に関わる機序を説明する図である。
【図１６】従来のマイクロニューログラフィによる血管交感神経活動（MSA）と拡張期血圧との関係を示す図である。
【発明を実施するための形態】
【００２５】
本発明に係る圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法は、圧反射機能に関わる血管の交感神経活動である圧反射性血管交感神経活動（Vascular BRSA：Vascular Baroreflex-Related Sympathetic Activity）を検出するとともに、当該圧反射性血管交感神経活動を指標化して評価するのに好適なものである。
【００２６】
本願発明者は、上述した課題に鑑みて鋭意研究した結果、生体から非侵襲的に測定された脈波データのみを使用することにより、簡易的かつ非侵襲的に、信頼性の高い圧反射性血管交感神経活動を個別に検出できることを見出した。
【００２７】
以下、本発明に係る圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法の一実施形態について図面を用いて説明する。
【００２８】
図１に示すように、本実施形態の圧反射性血管交感神経活動検出装置１には、耳や指等における生体動脈の脈波を検出するための脈波検出手段２が接続されている。また、圧反射性血管交感神経活動検出装置１は、本実施形態の圧反射性血管交感神経活動検出プログラム１ａや各種のデータを記憶する記憶手段３と、この記憶手段３を制御するとともに各種のデータを取得して演算処理を実行する演算処理手段４とを有している。以下、各構成手段について詳細に説明する。
【００２９】
なお、本実施形態において、圧反射性血管交感神経活動検出装置１は、記憶手段３および演算処理手段４を備えたパーソナルコンピュータによって構成されており、別途、脈波検出手段２が周辺機器として接続されている。しかしながら、装置構成は上記に限定されるものではない。例えば、本実施形態の圧反射性血管交感神経活動検出プログラム１ａをアプリ化してスマートフォンやタブレット端末にインストールし、それらを本実施形態の圧反射性血管交感神経活動検出装置１として機能させてもよい。
【００３０】
脈波検出手段２は、光電式容積脈波記録法（PPG：Photo-Plethysmography）により、動脈における脈動性の容積変化を示す容積脈波を非侵襲的に検出するものである。本実施形態において、脈波検出手段２は、図１に示すように、生体に装着されて光量を検出する光電式センサ２１と、この光電式センサ２１からの出力信号を増幅して脈波データとして出力する脈波アンプ２２とを有している。
【００３１】
光電式センサ２１は、生体に固定されるＬＥＤ（Light Emitting Diode）等の発光部２１ａと、この発光部２１ａと対向する位置に生体を挟んで配置されるフォトダイオード等の受光部２１ｂとを備えている。そして、発光部２１ａから発光されて生体を透過した透過光量を受光部２１ｂで検出し、脈波アンプ２２から脈波データとして時系列で出力するようになっている。
【００３２】
なお、本実施形態において、脈波アンプ２２は、脈波データとして光電式容積脈波を出力している。しかしながら、この構成に限定されるものではなく、脈波アンプ２２は、光電式容積脈波を異なる値に変換してから出力してもよい。具体的には、脈波アンプ２２は、規準化された光電式容積脈波である規準化脈波容積（Normalized Pulse Volume：以下、「NPV」という場合がある）を脈波データとして出力してもよい。あるいは、脈波アンプ２２は、規準化脈波容積の自然対数（lnNPV）のように、対数化された規準化脈波容積を脈波データとして出力してもよい。なお、規準化脈波容積（NPV）の算出方法としては、光電式容積脈波を直流成分と交流成分とに分離するとともに、当該交流成分の振幅を同時刻における前記直流成分の平均値で除算する方法がある。しかしながら、当該方法に限定されるものではなく、算出結果としての規準化脈波容積（NPV）が等しいものであればよい。
【００３３】
また、本実施形態において、脈波検出手段２は、生体動脈の脈波データを耳から取得しているが、検出対象部位は耳に限定されるものではない。すなわち、脈波検出手段２によって生体動脈の脈波を検出できる部位であればよく、例えば指、鼻、掌、手首、腕、足等を検出対象部位とすることもできる。また、本実施形態において、受光部２１ｂとしては、フォトダイオードのような透過型の光センサを使用しているが、フォトトランジスタのような反射型の光センサを使用してもよい。
【００３４】
さらに、脈波検出手段２は、在宅環境等の日常生活においても、簡単かつ手軽に脈波をモニタリングできるように、ウェアラブルに構成されていることが好ましい。例えば、耳に取り付け可能なイヤフォンタイプやイヤリングタイプ、手首や上腕等に取り付け可能なブレスレットタイプ、指に取り付け可能な指輪タイプ等である。このような脈波検出手段２によれば、装着する際の負担が少なく、モニタリング中も何ら支障なく日常生活を過ごすことができる。
【００３５】
記憶手段３は、各種のデータを記憶するとともに、演算処理手段４が演算処理を行う際のワーキングエリアとして機能するものである。本実施形態において、記憶手段３は、ハードディスク、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、フラッシュメモリ等で構成されており、図１に示すように、プログラム記憶部３１と、脈波データ記憶部３２と、拍間隔記憶部３３と、心電図データ記憶部３４と、ＶＢＲＳＡ系列記憶部３５と、ＶＢＲＳＡ値記憶部３６とを有している。
【００３６】
プログラム記憶部３１には、本実施形態の圧反射性血管交感神経活動検出プログラム１ａがインストールされている。そして、この圧反射性血管交感神経活動検出プログラム１ａが演算処理手段４によって実行されることにより、パーソナルコンピュータやスマートフォン、タブレット端末等のコンピュータを後述する各構成部として機能させるようになっている。
【００３７】
なお、圧反射性血管交感神経活動検出プログラム１ａの利用形態は、上記構成に限られるものではない。例えば、ＣＤ−ＲＯＭやＵＳＢメモリ等のように、コンピュータで読み取り可能な非一時的な記録媒体に圧反射性血管交感神経活動検出プログラム１ａを記憶させておき、当該記録媒体から直接読み出して実行してもよい。また、外部サーバ等からクラウドコンピューティング方式やＡＳＰ（Application Service Provider）方式等で利用してもよい。
【００３８】
脈波データ記憶部３２は、後述する脈波データ取得部４１によって取得された脈波データまたは後述する脈波データ変換部４２によって変換された脈波データを時系列で記憶するものである。拍間隔記憶部３３は、後述する拍間隔取得部４３によって取得された拍間隔を時系列で記憶するものである。心電図データ記憶部３４は、心電計等の心電図取得手段（図示せず）によって取得された心電図データを時系列で記憶するものである。
【００３９】
また、ＶＢＲＳＡ系列記憶部３５は、後述するＶＢＲＳＡ系列検出部４４によって検出された圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）を記憶するものである。さらに、ＶＢＲＳＡ値記憶部３６は、後述するＶＢＲＳＡ値算出部４６によって算出された圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を記憶するものである。
【００４０】
演算処理手段４は、ＣＰＵ（Central Processing Unit）等から構成されており、圧反射性血管交感神経活動検出プログラム１ａを実行することにより、脈波データ取得部４１と、脈波データ変換部４２と、拍間隔取得部４３と、ＶＢＲＳＡ系列検出部４４と、理論値算出部４５と、ＶＢＲＳＡ値算出部４６として機能するようになっている。以下、各構成部についてより詳細に説明する。
【００４１】
脈波データ取得部４１は、脈波検出手段２から生体動脈の脈波データを取得するものである。本実施形態において、脈波データ取得部４１は、脈波検出手段２から光電式容積脈波を脈波データとして取得する。そして、当該光電式容積脈波を脈波データ変換部４２に提供し規準化脈波容積（NPV）に変換させるようになっている。
【００４２】
なお、脈波データ取得部４１は、上記構成に限定されるものではなく、脈波検出手段２側で予め光電式容積脈波が規準化脈波容積（NPV）や、対数化された規準化脈波容積（lnNPV）に変換されている場合、それらを脈波データとして取得してもよい。また、光電式容積脈波をそのまま脈波データとして使用してもよい。これらの場合、脈波データ取得部４１は、脈波検出手段２から取得した脈波データを脈波データ変換部４２に提供することなく、直接脈波データ記憶部３２に時系列で格納する。
【００４３】
また、本実施形態において、脈波データ取得部４１は、脈波検出手段２から直接脈波データを取得しているが、この構成に限定されるものではない。例えば、他の装置によって測定された脈波データについては、ＵＳＢメモリ等の記録媒体から取り込んでもよい。また、データサーバ等に保存されている脈波データについては、インターネット等を介してダウンロードしてもよい。
【００４４】
脈波データ変換部４２は、脈波データとしての光電式容積脈波を規準化脈波容積（NPV）や対数化された規準化脈波容積（lnNPV）等の異なる値に変換するものである。本発明において、圧反射性血管交感神経活動を検出するための脈波データとしては、光電式容積脈波よりも規準化脈波容積（NPV）や対数化された規準化脈波容積（lnNPV）を用いることが好ましい。以下、その理由について説明する。
【００４５】
まず、図２（ａ）に示すように、生体動脈における光電式容積脈波を測定する場合、ランバート・ベールの法則（Lambert-Beer's law）に従えば、光源からの光が、一層で成り立つ溶液中を通過する際の吸収は、溶液の濃度および通過する液層の厚さに比例する。したがって、図２（ｂ）に示すように、光電式容積脈波の交流成分（PGac）は、その直流成分（PGdc）における動脈の脈動成分を増幅することにより算出される。
【００４６】
しかしながら、図２（ａ）、（ｂ）に示すように、実際の生体には、光が生体を透過する際の液層として、動脈が拍動していないときの非拍動成分や、生体の組織にあたる固定成分等の個人差成分が存在する。このため、光電式容積脈波の交流成分（PGac）は、異なる個人間では比較できない相対値としてしか使用できない。そこで、異なる個人間でも比較可能な絶対値として、光電式容積脈波の交流成分（PGac）をその直流成分（PGdc）で除算することにより脈動に伴う血液容積を規準化した規準化脈波容積（NPV）が知られている（非特許文献１）。
【００４７】
当該規準化脈波容積（NPV）は、上記のとおり、同じ単位（mV）同士の除算により算出されるため、脈動に伴う血液容積を無次元化した絶対量といえる。したがって、ある個人について異なる日時に取得した測定結果のみならず、異なる個人について取得された測定結果も比較することができる。また、生体の組織成分に起因する誤差が除外されるというメリットを有する。したがって、本発明においては、圧反射性血管交感神経活動を検出するための脈波データとして、規準化脈波容積（NPV）や対数化された規準化脈波容積（lnNPV）を使用することが推奨される。
【００４８】
なお、規準化脈波容積（NPV）は、ΔＩ／Ｉによって評価されるところ、脈動変化分の透過光量（ΔＩ）は、脈波の交流成分の振幅によって決定される。また、生体（組織＋血液）の透過光量（Ｉ）は、同時刻における脈波の直流成分の平均値によって決定される。したがって、脈波データ変換部４２は、脈波データ取得部４１によって取得された脈波データ（光電式容積脈波）に基づいて、脈波の交流成分の振幅を同時刻における当該脈波の直流成分の平均値で除算してなる規準化脈波容積（NPV）を毎拍ごとに算出し、脈波データ記憶部３２に時系列で格納するようになっている。
【００４９】
また、脈波データ変換部４２は、毎拍ごとに算出した規準化脈波容積（NPV）をさらに対数化し、当該対数化された規準化脈波容積（lnNPV）を脈波データ記憶部３２に時系列で格納するようにしてもよい。一方、上記のとおり、脈波データ取得部４１が、脈波検出手段２側で予め変換された規準化脈波容積（NPV）や、対数化された規準化脈波容積（lnNPV）を脈波データとして取得する場合、脈波データ変換部４２は不要である。
【００５０】
拍間隔取得部４３は、脈波データに対応する拍間隔を取得するものである。本発明において、「拍間隔」とは、図３に示すように、脈波データにおける脈拍間隔（PI）、脈波データの交流成分が立ち上がる部分の間隔（RI：以下、「立ち上がり間隔」という）、および心電図データにおける心拍間隔（RR間隔）等のように、拍動の一拍分に相当する全ての時間間隔を含む概念である。
【００５１】
また、本発明において、「脈波データに対応する拍間隔」とは、脈波データが取得された時間帯と同じ時間帯における拍間隔のみならず、脈波データが取得された時間帯よりも一拍分遅れた拍間隔を含む概念である。具体的には、拍間隔が心拍間隔の場合、図３に示すように、脈波Ａ、Ｂ、Ｃに対して、心拍間隔ａ、ｂ、ｃのみならず、心拍間隔ｂ、ｃ、ｄも「脈波データに対応する拍間隔」に当たる。
【００５２】
したがって、本実施形態において、拍間隔取得部４３は、脈波データ記憶部３２に格納された脈波データを参照し、当該脈波データに対応する脈拍間隔や立ち上がり間隔を取得し、時系列で拍間隔記憶部３３に格納する。なお、本実施形態では、脈拍間隔や立ち上がり間隔を判別しにくい老人等にのために、別途、心電図データ記憶部３４に心電図データを用意している。この場合、拍間隔取得部４３は心電図データを参照し、脈波データに対応する心拍間隔を取得し、時系列で拍間隔記憶部３３に格納するようになっている。
【００５３】
ＶＢＲＳＡ系列検出部４４は、脈波データから圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）を検出するものである。この圧反射性血管交感神経活動系列は、心臓の神経性圧反射機能による影響を受けるものではなく、血管の神経性圧反射機能に関わる血管の交感神経系により影響を受けるものである。以下、脈波データのうち規準化脈波容積（NPV）からＶＢＲＳＡ系列を検出する方法について具体的に説明するが、他の脈波データについても同様に考えることができる。
【００５４】
上述したとおり、規準化脈波容積（NPV）は、異なる個人間での測定結果を比較でき、生体の組織成分に起因する誤差を除外するというメリットを有する。しかしながら、規準化脈波容積（NPV）は、心臓の圧反射機能（交感神経）を機序とする一回拍出量（血液容積の脈動変化）と、血管の圧反射機能（交感神経）を機序とする血管抵抗（血管収縮度）との双方から影響を受ける。
【００５５】
このため、図４（ａ），（ｂ）に示すように、収縮期血圧（SBP）と規準化脈波容積（NPV）とが正相関のループを繰り返すデータ系列については、心臓の圧反射機能が優位に機能しており、一回拍出量に依存する心臓圧反射優位系列と考えられる。一方、上記正相関のループから大きく逸脱するデータ系列については、規準化脈波容積（NPV）が全く異なる値を示すため、血管の圧反射機能が優位に機能しており、血管抵抗に依存する血管圧反射優位系列と考えられる。なお、図４（ｂ）において、心臓圧反射優位系列のデータは白丸（○）で示し、血管圧反射優位系列のデータは黒丸（●）で示す。
【００５６】
つまり、図４（ｂ）の一点鎖線で示すように、収縮期血圧（SBP）が同じ値であっても、心臓の圧反射機能と血管の圧反射機能のうち、どちらの機序が優位に働いているかに応じて、規準化脈波容積（NPV）の値が大きく異なる。このため、規準化脈波容積（NPV）をそのまま使用しても、指標としての独立性を欠くこととなる。
【００５７】
また、図５（ａ），（ｂ）に示すように、収縮期血圧（SBP）と拍間隔（心拍間隔RRI）との散布図、および規準化脈波容積（NPV）と拍間隔（心拍間隔RRI）との散布図においても、正相関のループを繰り返すデータ系列（図５中の白丸）から逸脱するデータ系列（図５中の黒丸）が存在する。したがって、規準化脈波容積（NPV）から血管の圧反射機能に関わる圧反射性血管交感神経活動を独立して検出するには、前記逸脱するデータ系列を検出できればよい。
【００５８】
そこで、図６の実測データにおける時間帯Ａに着目すると、図７（ａ）に示すように、規準化脈波容積（NPV）および拍間隔（脈拍間隔IBI）からなる３拍以上のデータ系列において正相関が成立しているため、心臓の神経性圧反射系列といえる（特許文献１）。なお、心臓の神経性圧反射感度（Cardiac nBRS：Cardiac neural Baroreflex Sensitivity）は、規準化脈波容積（NPV）と拍間隔（IBI）との相関を示す回帰直線の傾き（回帰係数）によって算出され、図７（ａ）に示す例では、１９．２（ms/unitless）である。一方、図６の実測データにおける時間帯Ｂでは、図７（ｂ）に示すように、拍間隔（IBI）が３拍連続して減少するのに対して、規準化脈波容積（NPV）が３拍目に増加している。
【００５９】
すなわち、心臓の神経性圧反射系列から予想される３拍目の規準化脈波容積（NPV）の理論値は、図７（ｂ）に示すように、２拍目の値よりも小さい値であるのに対し、３拍目の規準化脈波容積（NPV）の実測値は２拍目の値よりも大きく、理論値から逸脱している。このような系列は、心臓の神経性圧反射系列における副交感自律神経活動によるものではなく、血管の交感神経活動の影響によるものと考えられるため、本発明に係る圧反射性血管交感神経活動系列といえる。なお、圧反射性血管交感神経活動（Vascular BRSA）を示す指標である後述のＶＢＲＳＡ値は、例えば下記式（１）によれば、｜理論値−実測値｜によって算出され、図７（ｂ）に示す例では、４．８４（unitless）である。
【００６０】
以上より、本実施形態において、ＶＢＲＳＡ系列検出部４４は、脈波データ記憶部３２内の脈波データと拍間隔記憶部３３内の拍間隔とを参照し、拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列を探索する。そして、ＶＢＲＳＡ系列検出部４４は、当該系列のうち、拍間隔と脈波データとの相関係数が、第ｎ−１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に相関係数が閾値以下となる系列を圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出し、ＶＢＲＳＡ系列記憶部３５に時系列で記憶させるようになっている。
【００６１】
なお、本発明において、圧反射性血管交感神経活動系列として検出する連続拍数ｎを３以上の自然数とするのは、第ｎ−１拍目までは交感神経系の影響を受けておらず、第ｎ拍目から圧反射機能に伴う血管の交感神経活動の影響を受けていることの信頼性を担保するためである。また、本実施形態において、正の閾値としては、０．８５を採用しているが、この値に限定されるものではない。
【００６２】
理論値算出部４５は、心臓の神経性圧反射系列から予想される第ｎ拍目の脈波データの理論値を算出するものである。本実施形態において、理論値算出部４５は、ＶＢＲＳＡ系列記憶部３５に記憶された圧反射性血管交感神経活動系列のデータを参照する。そして当該データ系列における、ｎ拍分の拍間隔とｎ−１拍分の脈波データとの回帰分析によって推定される、第ｎ拍目の脈波データの理論値を算出する。
【００６３】
具体的には、理論値算出部４５は、まず、圧反射性血管交感神経活動系列における、ｎ−１拍分の拍間隔と脈波データについて回帰分析（直線回帰等）を行い、回帰直線を算出する。そして、理論値算出部４５は、当該回帰直線に第ｎ拍目の拍間隔を代入して得られる値を第ｎ拍目の脈波データの理論値として算出するようになっている。
【００６４】
ＶＢＲＳＡ値算出部４６は、圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するものである。本実施形態において、ＶＢＲＳＡ値算出部４６は、理論値算出部４５によって算出された第ｎ拍目の脈波データの理論値と、ＶＢＲＳＡ系列記憶部３５に格納された第ｎ拍目の脈波データの実測値とに基づいて、圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出する。
【００６５】
具体的には、ＶＢＲＳＡ値算出部４６は、下記式（１）〜（３）のいずれかを用いてＶＢＲＳＡ値を算出するようになっている。
ＶＢＲＳＡ値＝｜ＰＶｔ−ＰＶｍ｜  ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）−ｌｎ（ＰＶｍ）｜  ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜  ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値
【００６６】
なお、本発明において、実測値とは、脈波データとして光電式容積脈波を使用する場合は、実際に測定された光電式容積脈波の値をいうものとする。また、脈波データとして規準化脈波容積（NPV）や対数化された規準化脈波容積（lnNPV）を使用する場合は、光電式容積脈波から変換された後の値をいうものとする。また、ＶＢＲＳＡ値の算出方法は、上記式（１）〜（３）に限定されるものではなく、第ｎ拍目の脈波データに関する理論値と実測値とを用いて、圧反射性血管交感神経活動を評価するのに適した指標を算出うるものであれば、どのような式でもよい。
【００６７】
また、本実施形態において、ＶＢＲＳＡ値算出部４６は、ＶＢＲＳＡ値を算出するにあたり、同種のデータ系列ごとに分けて個別に算出してもよい。具体的には、ＶＢＲＳＡ系列のうち拍間隔がｎ拍連続して増加する上昇系列と、ＶＢＲＳＡ系列のうち拍間隔がｎ拍連続して減少する下降系列である。
【００６８】
上昇系列について算出したＶＢＲＳＡ値は、血圧を上げ過ぎないように働くものの指標となる。また、下降系列について算出したＶＢＲＳＡ値は、血圧を下げ過ぎないように働くものの指標となる。よって、これらの各ＶＢＲＳＡ値は、疾患（例えば糖尿病等）に応じて異なる可能性があるため、別々に算出することによって、評価の精度が向上するものと考えられる。
【００６９】
一方、ＶＢＲＳＡ値算出部４６は、上昇系列および下降系列の合算値をＶＢＲＳＡ値として算出してもよい。具体的には、所定の時間（１分間等）内に上昇系列について算出されたＶＢＲＳＡ値の平均値と、所定の時間（１分間等）内に下降系列について算出されたＶＢＲＳＡ値の平均値とを合算または平均化したものをＶＢＲＳＡ値としてもよい。これにより、ＶＢＲＳＡ系列の検出数が少なくても算出できるケースが増大するとともに、圧反射性血管交感神経活動に関する理解しやすい指標としても利用することができる。
【００７０】
つぎに、本実施形態の圧反射性血管交感神経活動検出プログラム１ａによって実行される圧反射性血管交感神経活動検出装置１の作用および圧反射性血管交感神経活動検出方法について、図８を参照しつつ説明する。
【００７１】
本実施形態の圧反射性血管交感神経活動検出装置１または圧反射性血管交感神経活動検出プログラム１ａを用いて圧反射性血管交感神経活動を測定する場合、まず、脈波データ取得部４１が脈波検出手段２から生体動脈の脈波データを取得する（ステップＳ１：脈波データ取得ステップ）。この脈波データは、脈波検出手段２の態様に応じて、簡易的かつ非侵襲的に測定されるため、被測定者にかかる負担が軽減される。また、ウェアラブルな脈波検出手段２であれば、日常的または継続的な測定も簡単になる。
【００７２】
つぎに、脈波データ変換部４２が、ステップＳ１で取得された脈波データとしての光電式容積脈波を規準化脈波容積（NPV）等に変換する（ステップＳ２：脈波データ変換ステップ）。これにより、変換後の脈波データにおいては、上述した非拍動成分や固定成分等の個人差成分が捨象される。このため、ある個人について異なる日時に取得した測定結果のみならず、異なる個人について取得された測定結果も比較可能となる。また、生体の組織成分に起因する誤差が除外される。
【００７３】
なお、本実施形態では、脈波データ取得部４１が光電式容積脈波を脈波データとして取得するため、圧反射性血管交感神経活動検出装置１側で脈波データを変換している（ステップＳ２）。しかしながら、上述したとおり、脈波検出手段２から規準化脈波容積（NPV）や、対数化された規準化脈波容積（lnNPV）を脈波データとして取得した場合、ステップＳ２はスキップされ、ステップＳ１からステップＳ３へジャンプする。
【００７４】
つづいて、拍間隔取得部４３が、脈波データに対応する拍間隔を取得すると（ステップＳ３：拍間隔取得ステップ）、ＶＢＲＳＡ系列検出部４４が、拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、拍間隔と脈波データとの相関係数が、第ｎ−１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に相関係数が閾値以下となる系列を圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出する（ステップＳ４：ＶＢＲＳＡ系列検出ステップ）。これにより、圧反射機能に関わる自律神経系のうち、心臓の副交感神経系ではなく、血管の交感神経系の影響を受ける圧反射性血管交感神経活動を示すデータ系列が別個独立に特定される。
【００７５】
つぎに、理論値算出部４５が、圧反射性血管交感神経活動系列における、ｎ拍分の拍間隔とｎ−１拍分の脈波データとの回帰分析によって、第ｎ拍目の脈波データの理論値を算出する（ステップＳ５：理論値算出ステップ）。これにより、仮に第ｎ拍目においても、心臓の神経性圧反射機能が働いていた場合に測定されたはずの脈波データが予想値として算出される。
【００７６】
最後に、ＶＢＲＳＡ値算出部４６が、ステップＳ５で算出された理論値と、第ｎ拍目の脈波データの実測値とに基づいて、ＶＢＲＳＡ値を算出する（ステップＳ６：ＶＢＲＳＡ値算出ステップ）。これにより、当該ＶＢＲＳＡ値を用いて、圧反射性血管交感神経活動を簡易的かつ客観的に評価することが可能となる。また、上記式（１）〜（３）のいずれかを用いることにより、圧反射性血管交感神経活動を評価するのに適したＶＢＲＳＡ値が算出される。
【００７７】
また、本実施形態において、ＶＢＲＳＡ値算出部４６は、上昇系列と下降系列とのそれぞれについて、ＶＢＲＳＡ値を個別に算出してもよい。これにより、両系列におけるＶＢＲＳＡ値が、疾患に応じて異なる傾向を示すようであれば、評価の精度が向上する。さらに、ＶＢＲＳＡ値算出部４６は、上昇系列および下降系列の合算値をＶＢＲＳＡ値として算出してもよい。これにより、ＶＢＲＳＡ系列の検出数が少なくても算出できるケースが増大する。また、圧反射性血管交感神経活動を示す指標が分かりやすい形で表現される。
【００７８】
以上のような本発明に係る本実施形態によれば、以下のような効果を奏する。
１．簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を個別に検出することができる。
２．圧反射性血管交感神経活動をＶＢＲＳＡ値として指標化し、圧反射性血管交感神経活動を簡易的かつ客観的に評価することができる。
３．上記式（１）〜（３）のいずれかを用いることにより、圧反射性血管交感神経活動を評価するのに適したＶＢＲＳＡ値を算出することができる。
４．上昇系列と下降系列とのそれぞれについてＶＢＲＳＡ値を個別に算出することで、疾患に応じて異なる傾向等を把握でき、評価の精度を向上することができる。
５．上昇系列および下降系列の合算値をＶＢＲＳＡ値として算出することで、ＶＢＲＳＡ系列の検出数が少なくても算出でき、圧反射性血管交感神経活動を示す指標を分かりやすく表現することができる。
６．無次元の絶対量である規準化脈波容積（NPV）や、対数化された規準化脈波容積（lnNPV）を用いることにより、生体の組織成分に起因する誤差を除外することができるため、同一人の測定結果のみならず、異なる個人間の測定結果も比較することができ、広く一般的に圧反射性血管交感神経活動の動向を把握することができる。
７．脈波検出手段２をウェアラブルに構成することで、日常生活においても、血圧を調節する自律神経機能を継続的にモニターすることができる。
８．ＶＢＲＳＡ値は、血圧などの他の測定値との関係を関数化することによって、血管の圧反射機能を示す指標（血管神経性圧反射感度）として使用することもできる。
９．脈波データからは、心臓の神経性圧反射感度を別途測定できるため、圧反射性血管交感神経活動との関係性を評価することで、被測定者のそれぞれについて血圧調節機能や血圧状態を幅広く調査・研究することができる。
１０．心電図データを用いることにより、脈波における拍間隔が読み取りにくい被測定者からも正確な拍間隔を取得することができる。
１１．圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）の検出数や単位時間当たりの検出頻度についても、医学的に有用な情報として提供できる可能性がある。
【００７９】
つぎに、本発明に係る圧反射性血管交感神経活動検出装置１、圧反射性血管交感神経活動検出プログラム１ａおよび圧反射性血管交感神経活動検出方法の具体的な実施例について説明する。
【実施例１】
【００８０】
本実施例１では、精神的ストレスによって変動する血行動態に基づいて、本発明に係る圧反射性血管交感神経活動と、血管の交感神経活動（血管の神経性圧反射機能）によって制御される血管抵抗との関連性を確認するための実験を行った。【００８１】
具体的には、仰臥位の被験者１名（２３歳、男性）に対し、血圧（BP）を測定するための血圧計（型名：MUB101、製造元：株式会社メディセンス）を指に装着した。また、光電式容積脈波を測定するための脈波検出手段２（型名：EBRS-01J、製造元：有限会社ニッチプロダクト）を被験者の耳に装着した。さらに、ドップラー法により頸動脈における心拍出量（Cardiac output）を測定するための超音波プローブ（型名：UNEXEF38G、製造元：株式会社ユネクス）を被験者の首に固定した。【００８２】
そして、上述した各データを９分間、時系列で同時測定した。当該９分間のうち、最初の３分間は被験者を安静（Rest）にした。また、次の３分間で、被験者に精神的ストレスを負荷する目的で難しい暗算を実施させた（Stress）。そして、最後の３分間は、再び被験者を安静にして回復を図った（Recovery）。【００８３】
そして、測定された心拍出量を血管内皮機能検査装置（型名：UNEXEF38G、製造元：株式会社ユネクス）により処理し、平均血圧（MBP）と心拍出量（Cardiac output）を算出するとともに、血行力学の式（平均血圧＝心拍出量×血管抵抗）に基づいて、頸動脈の血管抵抗（Carotid arterial resistance）を算出した。また、測定された光電式容積脈波を本発明に係る圧反射性血管交感神経活動検出装置１によって処理し、規準化脈波容積（NPV）、脈拍間隔（IBI）およびＶＢＲＳＡ値（Vascular BRSA）を算出した。この結果を図９に示す。【００８４】
図９は、被験者について得られた測定結果を時系列で示す図である。また、当該測定結果のうち、ＶＢＲＳＡ値（Vascular BRSA）と頸動脈の血管抵抗（Carotid arterial resistance）のそれぞれについて、１０秒毎の平均値をプロットした散布図を図１０に示す。なお、図１０において、白丸は安静時（Rest）、黒丸はストレス負荷時（Stress）、白三角は回復時（Recovery）のデータである。【００８５】
図１０に示すように、ストレス負荷時には、安静時に比較して脈拍間隔が半分程度に低減し、脈拍数としては２倍程度に増大している。また、ストレス負荷時には、安静時に比較して血管抵抗が増大し、血管が緊張して固くなっていることがわかる。このとき、ＶＢＲＳＡ値（Vascular BRSA）の動向は、血管抵抗が低い安静時には高い値を示し、血管がそれ以上弛緩しないように働いている。一方、血管抵抗が高いストレス負荷時には低い値を示し、あまり働いていない。このような結果は、従来のマイクロニューログラフィーによって測定された血管の交感神経活動と一致する。【００８６】
以上の本実施例１によれば、本発明に係る圧反射性血管交感神経活動と、血管の交感神経活動（血管の神経性圧反射機能）によって制御される血管抵抗との間には、理論通りの関係性があることが確認された。また、本発明に係るＶＢＲＳＡ値は、圧反射性血管交感神経活動を示す指標として利用できることが示された。
【実施例２】
【００８７】
本実施例２では、精神的ストレスに対して血圧を上昇させる反応機序が、心臓の圧反射機能または血管の圧反射機能のどちらが優位な人であっても、本発明に係る圧反射性血管交感神経活動が、血管抵抗と関連性を有するか否かを確認する実験を行った。【００８８】
具体的には、２０〜２５歳の男性９名および女性２名の各被験者に対して、実施例１と同様の測定および解析を実施し、上述した各種データを算出した。また、上述した精神的ストレスに対して血圧を上昇させる反応機序は、個人によってどちらが優位に作用するかが異なる。そこで、本実施例２では、主に心拍出量によって血圧を上昇させる（したがって、血管抵抗は低下する）心臓優位反応群と、主に血管抵抗によって血圧を上昇させる血管優位反応群とのいずれかに全被験者を分類した。【００８９】
そして、心臓優位反応群と血管優位反応群とのそれぞれについて、安静時（Rest）、ストレス負荷時（Stress）および回復時（Recovery）における、平均血圧（MBP）、頸動脈の血管抵抗（Carotid arterial resistance）およびＶＢＲＳＡ値（Vascular BRSA）の平均値と標準偏差をグラフ化した。その結果を図１１（ａ）〜（ｃ）に示す。なお、黒丸は血管優位反応群のデータであり、白丸は心臓優位反応群のデータを示している。【００９０】
また、全被験者の実測データを用いて作成された、ＶＢＲＳＡ値（Vascular BRSA）と頸動脈の血管抵抗（Carotid arterial resistance）との相関図を図１２に示す。本実施例２において、両者のピアソン相関係数は−０．６１であり、１％水準で有意であった。図１１および図１２に示すように、本発明に係るＶＢＲＳＡ値（Vascular BRSA）は、心臓優位反応群と血管優位反応群の双方において、頸動脈の血管抵抗（Carotid arterial resistance）との関連性が確認された。【００９１】
以上の本実施例２によれば、精神的ストレスに対して血圧を上昇させる反応機序が、心臓の圧反射機能または血管の圧反射機能のどちらが優位な人であっても、本発明に係る圧反射性血管交感神経活動が血管抵抗と関連性を有することが確認された。よって、本発明に係るＶＢＲＳＡ値は、反応機序に係る個人差によらず、圧反射性血管交感神経活動を評価しうる指標として利用できることが示された。
【実施例３】
【００９２】
本実施例３では、本発明に係るＶＢＲＳＡ値と、血管の弛緩による拡張期血圧との関連性に基づいて、ＶＢＲＳＡ値が反射性血管交感神経活動を評価する指標として信頼性を有するものであるか否かを確認する実験を行った。【００９３】
具体的には、被験者（２２歳、男性）に対し、血圧変動を測定するための携帯型自動血圧計（型名：TM-2431、製造元：株式会社エー・アンド・デイ）を上腕に装着した。また、光電式容積脈波を測定するための脈波検出手段２（型名：EBRS-01J、製造元：有限会社ニッチプロダクト）を被験者の耳に装着した。そして、上述した各データを就寝前１時間以上、就寝時、そして起床後１時間以上という長時間に渡り、時系列で同時測定した。【００９４】
そして、測定された血圧データを血圧解析ソフト（型名：TM-2430-15、製造元：株式会社エー・アンド・デイ）により処理し、拡張期血圧（DBP）を算出した。また、測定された光電式容積脈波を本発明に係る圧反射性血管交感神経活動検出装置１によって処理し、脈拍間隔（IBI）、規準化脈波容積（NPV）およびＶＢＲＳＡ値（Vascular BRSA）を算出した。その結果を図１３に示す。【００９５】
なお、被験者は、測定期間中における自身の行動事象（Event）について、その内容と時間を記録した。図１３において、Event1として記録されているのは、トイレに行った時刻である。また、Event2は就寝時刻であり、Event3は起床時刻である。本実施例３において、就寝時刻から起床時刻までの期間を就寝期（Sleep）とし、それ以外の期間を覚醒期（Awake）とした。【００９６】
一般的に、健康な若者においては、睡眠中、脈拍数が低下（脈拍間隔が増加）するとともに、血管の弛緩により血圧が低下することが想定されている。この点、本実施例３においても、図１３に示すように、就寝期（Sleep）においては、脈拍間隔（IBI）が増加するとともに、拡張期血圧（DBP）が低下しており、上記想定どおりの結果であることが確認された。【００９７】
また、図１４は、図１３に示す実測データのうち血圧測定時間における、本発明に係るＶＢＲＳＡ値の平均値（Vascular BRSA）と、拡張期血圧（DBP）との相関を示すグラフである。図１４に示すように、ＶＢＲＳＡ値と拡張期血圧との間には、高い相関が確認され、その相関係数は、−０．９４であった。なお、ＶＢＲＳＡ値と拡張期血圧との相関を示す回帰直線の傾きは、血管の圧反射機能を示す指標（血管神経性圧反射感度）として使用することもできる。【００９８】
以上の本実施例３によれば、本発明に係るＶＢＲＳＡ値は、日常生活上においても、血管の弛緩による拡張期血圧と対応した反射性血管交感神経活動を評価する指標として高い信頼性を有することが実証された。【００９９】
なお、本発明に係る圧反射性血管交感神経活動検出装置１、圧反射性血管交感神経活動検出プログラム１ａおよび圧反射性血管交感神経活動検出方法は、前述した実施形態および実施例に限定されるものではなく、適宜変更することができる。【０１００】
例えば、本発明に係るＶＢＲＳＡ値の数値範囲に対応させて、その評価を示す表示テーブル（例：良好、普通、悪い）を記憶手段３に格納しておき、ＶＢＲＳＡ値算出部４６によって算出されたＶＢＲＳＡ値に基づき、圧反射性血管交感神経活動に関する評価を表示手段や印刷手段から出力し、被測定者に提示するようにしてもよい。【０１０１】
また、上述した本実施形態では、圧反射性血管交感神経活動を検出するための脈波データとして、規準化脈波容積（NPV）を使用しているが、必ずしもこれに限定されるものではない。すなわち、個人間で値を比較できない相対値としてのみ使用する場合には、光電式容積脈波の値をそのまま使用してもよい。
【産業上の利用可能性】
【０１０２】
血圧は、後期高齢化社会の中で増加し続ける心血管疾患における重要なリスクファクターである。信頼性の高い圧反射性血管交感神経活動の測定を日常生活上にまで応用することで、血圧を調節する仕組みをより詳細にモニターし、調べることができる。したがって、慢性心不全を含む生活習慣病全体のプレクリニカルな段階での疾病予測、スクリーニング、診断、および長期の経過観察に役立ち、医療現場への波及効果は極めて広範である。【０１０３】
また、本発明に係る反射性血管交感神経活動測定装置１が、スマートフォンやタブレット端末等の携帯端末サイズに小型化されれば、子供からお年寄りまで手軽に扱うことができる健康機器として家庭用・個人用に広く普及することも考えられる。さらに、ウェアラブルな脈波検出手段２を組み合わせて使用すれば、日常生活におけるモニタリングも容易になるものと考えられる。
【符号の説明】
【０１０４】
１  圧反射性血管交感神経活動検出装置
１ａ  圧反射性血管交感神経活動検出プログラム
２  脈波検出手段
３  記憶手段
４  演算処理手段
２１  光電式センサ
２１ａ  発光部
２１ｂ  受光部
２２  脈波アンプ
３１  プログラム記憶部
３２  脈波データ記憶部
３３  拍間隔記憶部
３４  心電図データ記憶部
３５  ＶＢＲＳＡ系列記憶部
３６  ＶＢＲＳＡ値記憶部
４１  脈波データ取得部
４２  脈波データ変換部
４３  拍間隔取得部
４４  ＶＢＲＳＡ系列検出部
４５  理論値算出部
４６  ＶＢＲＳＡ値算出部
【図１】
2018181851.tif 000003
【図２】
2018181851.tif 000004
【図３】
2018181851.tif 000005
【図４】
2018181851.tif 000006
【図５】
2018181851.tif 000007
【図６】
2018181851.tif 000008
【図７】
2018181851.tif 000009
【図８】
2018181851.tif 000010
【図９】
2018181851.tif 000011
【図１０】
2018181851.tif 000012
【図１１】
2018181851.tif 000013
【図１２】
2018181851.tif 000014
【図１３】
2018181851.tif 000015
【図１４】
2018181851.tif 000016
【図１５】
2018181851.tif 000017
【図１６】
2018181851.tif 000018
【国際調査報告】
2018181851.tif 000019
2018181851.tif 000020