JP7081831B2 - 圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法 - Google Patents
圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法 Download PDFInfo
- Publication number
- JP7081831B2 JP7081831B2 JP2019510214A JP2019510214A JP7081831B2 JP 7081831 B2 JP7081831 B2 JP 7081831B2 JP 2019510214 A JP2019510214 A JP 2019510214A JP 2019510214 A JP2019510214 A JP 2019510214A JP 7081831 B2 JP7081831 B2 JP 7081831B2
- Authority
- JP
- Japan
- Prior art keywords
- vbrsa
- pressure
- pulse wave
- reflex
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000008035 nerve activity Effects 0.000 title claims description 180
- 230000002118 vasosympathetic effect Effects 0.000 title claims description 124
- 238000001514 detection method Methods 0.000 title claims description 63
- 230000011514 reflex Effects 0.000 claims description 118
- 230000002792 vascular Effects 0.000 claims description 83
- 230000002889 sympathetic effect Effects 0.000 claims description 66
- 230000006870 function Effects 0.000 claims description 60
- 238000004364 calculation method Methods 0.000 claims description 52
- 210000004204 blood vessel Anatomy 0.000 claims description 50
- 210000005036 nerve Anatomy 0.000 claims description 23
- 230000001174 ascending effect Effects 0.000 claims description 20
- 230000007423 decrease Effects 0.000 claims description 20
- 238000000034 method Methods 0.000 claims description 19
- 210000001367 artery Anatomy 0.000 claims description 13
- 238000000611 regression analysis Methods 0.000 claims description 8
- 230000000630 rising effect Effects 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims 3
- 230000036772 blood pressure Effects 0.000 description 35
- 238000006243 chemical reaction Methods 0.000 description 25
- 230000000747 cardiac effect Effects 0.000 description 21
- 238000005259 measurement Methods 0.000 description 17
- 230000035487 diastolic blood pressure Effects 0.000 description 14
- 230000001537 neural effect Effects 0.000 description 13
- 238000013500 data storage Methods 0.000 description 12
- 230000035488 systolic blood pressure Effects 0.000 description 10
- 230000007246 mechanism Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 9
- 230000035945 sensitivity Effects 0.000 description 9
- 210000001519 tissue Anatomy 0.000 description 8
- 201000010099 disease Diseases 0.000 description 7
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 7
- 238000011156 evaluation Methods 0.000 description 7
- 230000003340 mental effect Effects 0.000 description 6
- 239000008280 blood Substances 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 5
- 210000001715 carotid artery Anatomy 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000011084 recovery Methods 0.000 description 5
- 208000006011 Stroke Diseases 0.000 description 4
- 230000035581 baroreflex Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 230000010349 pulsation Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 208000024172 Cardiovascular disease Diseases 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 210000002820 sympathetic nervous system Anatomy 0.000 description 3
- 125000002066 L-histidyl group Chemical group [H]N1C([H])=NC(C([H])([H])[C@](C(=O)[*])([H])N([H])[H])=C1[H] 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 210000003403 autonomic nervous system Anatomy 0.000 description 2
- 210000000467 autonomic pathway Anatomy 0.000 description 2
- 238000009530 blood pressure measurement Methods 0.000 description 2
- 206010008118 cerebral infarction Diseases 0.000 description 2
- 208000026106 cerebrovascular disease Diseases 0.000 description 2
- 210000005069 ears Anatomy 0.000 description 2
- 230000000004 hemodynamic effect Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000001734 parasympathetic effect Effects 0.000 description 2
- 210000005037 parasympathetic nerve Anatomy 0.000 description 2
- 230000000541 pulsatile effect Effects 0.000 description 2
- 230000029865 regulation of blood pressure Effects 0.000 description 2
- 210000000707 wrist Anatomy 0.000 description 2
- 206010005746 Blood pressure fluctuation Diseases 0.000 description 1
- 206010007558 Cardiac failure chronic Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 206010047139 Vasoconstriction Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000005574 cross-species transmission Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 210000001002 parasympathetic nervous system Anatomy 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 210000000578 peripheral nerve Anatomy 0.000 description 1
- 238000013186 photoplethysmography Methods 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 210000001774 pressoreceptor Anatomy 0.000 description 1
- 230000008430 psychophysiology Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 206010042772 syncope Diseases 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 230000003845 vascular endothelial function Effects 0.000 description 1
- 230000006442 vascular tone Effects 0.000 description 1
- 230000025033 vasoconstriction Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 230000002618 waking effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B10/00—Other methods or instruments for diagnosis, e.g. instruments for taking a cell sample, for biopsy, for vaccination diagnosis; Sex determination; Ovulation-period determination; Throat striking implements
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0295—Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4029—Detecting, measuring or recording for evaluating the nervous system for evaluating the peripheral nervous systems
- A61B5/4035—Evaluating the autonomic nervous system
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Vascular Medicine (AREA)
- Psychiatry (AREA)
- Hematology (AREA)
- Educational Technology (AREA)
- Child & Adolescent Psychology (AREA)
- Developmental Disabilities (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Hospice & Palliative Care (AREA)
- Psychology (AREA)
- Social Psychology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
本発明は、血圧を一定に保つ役割を果たす圧反射機能に関わる血管の交感神経活動を非侵襲的に検出する圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法に関するものである。
従来、血圧は、心血管疾患や脳梗塞等における重要なリスクファクターであることが知られている。このため、血圧の状態を日常的にモニターすることにより、上述した心血管疾患や脳梗塞等を予防するための指針や、プレクリニカル段階での診断材料等として利用することができる。しかしながら、血圧を測定する行為は、身体への負担が大きく面倒な作業でもあるため、現実的には日常的または継続的に血圧をモニタリングすることは困難である。
そこで近年、血圧を一定の範囲に保とうとする自律神経系の圧反射機能を評価する研究が行われている。当該圧反射機能の低下が、脳卒中、高血圧症および失神等におけるリスクファクターと考えられているためである。ただし、血圧は、血行力学上、心拍出量（一回拍出量×心拍数）と血管抵抗との積として定義されるところ、圧反射においては、上記心拍数に対する心臓の交感神経活動の寄与は見られないため、上記心拍数は心臓の副交感神経活動によって制御され、上記血管抵抗は血管の交感神経活動によって制御されている。このため、前記圧反射機能には、心臓の副交感神経活動（心臓の神経性圧反射機能）と、血管の交感神経活動（血管の神経性圧反射機能）との双方が関わっている。
具体的には、図１５に示すように、心臓の神経性圧反射機能では、血圧の上昇により血管壁が伸展すると、頸動脈等に存在する圧受容体細胞が発火する。そして、当該発火信号を受けた脳が心臓の副交感神経活動を介して心拍数を減少させて（すなわち、心拍間隔または脈拍間隔を延長させて）心拍出量を下げ、血圧を下げようとする。一方、血管の神経性圧反射機能は、血管の弛緩により血圧が低下すると、所定値以下にならないよう、血管の交感神経活動を賦活して血管を緊張させ、血圧を上げようとする。
ただし、心臓の神経性圧反射機能と血管の神経性圧反射機能との作用バランスには個人差があり、心臓の神経性圧反射機能が優位に働く人もいれば、血管の神経性圧反射機能が優位に働く人もいる。また、同じ個人であっても、例えば安静にしている場合と運動した後とでは、上記作用バランスが変わることもある。したがって、圧反射機能を総合的に評価するには、心臓の神経性圧反射機能と血管の神経性圧反射機能とを別個独立に測定する必要がある。
心臓の神経性圧反射機能を測定する技術としては、当該機能を感度として非侵襲的に測定する方法が存在する。例えば、血管壁の伸展による血管径の増加や減少を超音波エコーにより毎拍測定し、これとともに生じる心拍ないし脈拍間隔の延長ないし短縮として現れる心臓の圧反射系列に基づいて検出する、いわゆるシーケンス法が知られている。また、本願発明者は、血管径ではなく、規準化脈波容積（非特許文献１）による血液容積を用いて、心臓の神経性圧反射感度を測定する方法を提案している（特許文献１）。
一方、血管の神経性圧反射機能を測定する技術としては、従来、図１６に示すように、拡張期血圧または血管抵抗と、いわゆるマイクロニューログラフィによって侵襲的に測定された血管の交感神経活動（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）を用いることにより、生体の組織成分に起因する誤差を除外し、同一人の測定結果のみならず、異なる個人間の測定結果も比較可能とし、広く一般的に圧反射性血管交感神経活動の動向を把握するという課題を解決するために、前記脈波データは、光電式容積脈波、規準化された前記光電式容積脈波である規準化脈波容積、または対数化された前記規準化脈波容積のうちのいずれかであってもよい。
また、本発明に係る圧反射性血管交感神経活動検出方法は、簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出するという課題を解決するために、生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出方法であって、前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ－１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出ステップを有する。
また、本発明に係る圧反射性血管交感神経活動検出方法の一態様として、圧反射性血管交感神経活動をＶＢＲＳＡ値として指標化し、圧反射性血管交感神経活動を簡易的かつ客観的に評価するという課題を解決するために、前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ－１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出ステップと、第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出ステップとを有していてもよい。
さらに、本発明の一態様として、圧反射性血管交感神経活動を評価するのに適したＶＢＲＳＡ値を算出するという課題を解決するために、前記ＶＢＲＳＡ値算出ステップでは、下記式（１）～（３）のいずれかを用いて前記ＶＢＲＳＡ値が算出されてもよい。
ＶＢＲＳＡ値＝｜ＰＶｔ－ＰＶｍ｜ ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）－ｌｎ（ＰＶｍ）｜ ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜ ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値
ＶＢＲＳＡ値＝｜ＰＶｔ－ＰＶｍ｜ ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）－ｌｎ（ＰＶｍ）｜ ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜ ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値
また、本発明の一態様として、上昇系列と下降系列とのそれぞれについてＶＢＲＳＡ値を個別に算出することで、疾患に応じて異なる傾向等を把握し、評価の精度を向上するという課題を解決するために、前記ＶＢＲＳＡ値算出ステップでは、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出してもよい。
さらに、本発明の一態様として、上昇系列および下降系列の合算値をＶＢＲＳＡ値として算出することで、ＶＢＲＳＡ系列の検出数が少なくても算出し、圧反射性血管交感神経活動を示す指標を分かりやすく表現するという課題を解決するために、前記ＶＢＲＳＡ値算出ステップでは、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出してもよい。
本発明によれば、簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出することができる。
本発明に係る圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法は、圧反射機能に関わる血管の交感神経活動である圧反射性血管交感神経活動（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）を用いることにより、生体の組織成分に起因する誤差を除外することができるため、同一人の測定結果のみならず、異なる個人間の測定結果も比較することができ、広く一般的に圧反射性血管交感神経活動の動向を把握することができる。
７．脈波検出手段２をウェアラブルに構成することで、日常生活においても、血圧を調節する自律神経機能を継続的にモニターすることができる。
８．ＶＢＲＳＡ値は、血圧などの他の測定値との関係を関数化することによって、血管の圧反射機能を示す指標（血管神経性圧反射感度）として使用することもできる。
９．脈波データからは、心臓の神経性圧反射感度を別途測定できるため、圧反射性血管交感神経活動との関係性を評価することで、被測定者のそれぞれについて血圧調節機能や血圧状態を幅広く調査・研究することができる。
１０．心電図データを用いることにより、脈波における拍間隔が読み取りにくい被測定者からも正確な拍間隔を取得することができる。
１１．圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）の検出数や単位時間当たりの検出頻度についても、医学的に有用な情報として提供できる可能性がある。
１．簡易的かつ非侵襲的に、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を個別に検出することができる。
２．圧反射性血管交感神経活動をＶＢＲＳＡ値として指標化し、圧反射性血管交感神経活動を簡易的かつ客観的に評価することができる。
３．上記式（１）～（３）のいずれかを用いることにより、圧反射性血管交感神経活動を評価するのに適したＶＢＲＳＡ値を算出することができる。
４．上昇系列と下降系列とのそれぞれについてＶＢＲＳＡ値を個別に算出することで、疾患に応じて異なる傾向等を把握でき、評価の精度を向上することができる。
５．上昇系列および下降系列の合算値をＶＢＲＳＡ値として算出することで、ＶＢＲＳＡ系列の検出数が少なくても算出でき、圧反射性血管交感神経活動を示す指標を分かりやすく表現することができる。
６．無次元の絶対量である規準化脈波容積（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）を使用しているが、必ずしもこれに限定されるものではない。すなわち、個人間で値を比較できない相対値としてのみ使用する場合には、光電式容積脈波の値をそのまま使用してもよい。
血圧は、後期高齢化社会の中で増加し続ける心血管疾患における重要なリスクファクターである。信頼性の高い圧反射性血管交感神経活動の測定を日常生活上にまで応用することで、血圧を調節する仕組みをより詳細にモニターし、調べることができる。したがって、慢性心不全を含む生活習慣病全体のプレクリニカルな段階での疾病予測、スクリーニング、診断、および長期の経過観察に役立ち、医療現場への波及効果は極めて広範である。
また、本発明に係る反射性血管交感神経活動測定装置１が、スマートフォンやタブレット端末等の携帯端末サイズに小型化されれば、子供からお年寄りまで手軽に扱うことができる健康機器として家庭用・個人用に広く普及することも考えられる。さらに、ウェアラブルな脈波検出手段２を組み合わせて使用すれば、日常生活におけるモニタリングも容易になるものと考えられる。
１ 圧反射性血管交感神経活動検出装置
１ａ 圧反射性血管交感神経活動検出プログラム
２ 脈波検出手段
３ 記憶手段
４ 演算処理手段
２１ 光電式センサ
２１ａ 発光部
２１ｂ 受光部
２２ 脈波アンプ
３１ プログラム記憶部
３２ 脈波データ記憶部
３３ 拍間隔記憶部
３４ 心電図データ記憶部
３５ ＶＢＲＳＡ系列記憶部
３６ ＶＢＲＳＡ値記憶部
４１ 脈波データ取得部
４２ 脈波データ変換部
４３ 拍間隔取得部
４４ ＶＢＲＳＡ系列検出部
４５ 理論値算出部
４６ ＶＢＲＳＡ値算出部
１ａ 圧反射性血管交感神経活動検出プログラム
２ 脈波検出手段
３ 記憶手段
４ 演算処理手段
２１ 光電式センサ
２１ａ 発光部
２１ｂ 受光部
２２ 脈波アンプ
３１ プログラム記憶部
３２ 脈波データ記憶部
３３ 拍間隔記憶部
３４ 心電図データ記憶部
３５ ＶＢＲＳＡ系列記憶部
３６ ＶＢＲＳＡ値記憶部
４１ 脈波データ取得部
４２ 脈波データ変換部
４３ 拍間隔取得部
４４ ＶＢＲＳＡ系列検出部
４５ 理論値算出部
４６ ＶＢＲＳＡ値算出部
Claims (18)
- 生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出装置であって、
前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ－１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出部を有する、圧反射性血管交感神経活動検出装置。 - 前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ－１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出部と、
第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出部と
を有する、請求項１に記載の圧反射性血管交感神経活動検出装置。 - 前記ＶＢＲＳＡ値算出部は、下記式（１）～（３）のいずれかを用いて前記ＶＢＲＳＡ値を算出する、請求項２に記載の圧反射性血管交感神経活動検出装置；
ＶＢＲＳＡ値＝｜ＰＶｔ－ＰＶｍ｜ ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）－ｌｎ（ＰＶｍ）｜ ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜ ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値 - 前記ＶＢＲＳＡ値算出部は、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出する、請求項２または請求項３に記載の圧反射性血管交感神経活動検出装置。
- 前記ＶＢＲＳＡ値算出部は、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出する、請求項４に記載の圧反射性血管交感神経活動検出装置。
- 前記脈波データは、光電式容積脈波、規準化された前記光電式容積脈波である規準化脈波容積、または対数化された前記規準化脈波容積のうちのいずれかである、請求項１から請求項５のいずれかに記載の圧反射性血管交感神経活動検出装置。
- 生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出プログラムであって、
前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ－１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出部としてコンピュータを機能させる、圧反射性血管交感神経活動測定プログラム。 - 前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ－１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出部と、
第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出部としてコンピュータを機能させる、請求項７に記載の圧反射性血管交感神経活動検出プログラム。 - 前記ＶＢＲＳＡ値算出部は、下記式（１）～（３）のいずれかを用いて前記ＶＢＲＳＡ値を算出する、請求項８に記載の圧反射性血管交感神経活動検出プログラム；
ＶＢＲＳＡ値＝｜ＰＶｔ－ＰＶｍ｜ ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）－ｌｎ（ＰＶｍ）｜ ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜ ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値 - 前記ＶＢＲＳＡ値算出部は、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出する、請求項８または請求項９に記載の圧反射性血管交感神経活動検出プログラム。
- 前記ＶＢＲＳＡ値算出部は、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出する、請求項１０に記載の圧反射性血管交感神経活動検出プログラム。
- 前記脈波データは、光電式容積脈波、規準化された前記光電式容積脈波である規準化脈波容積、または対数化された前記規準化脈波容積のうちのいずれかである、請求項７から請求項１１のいずれかに記載の圧反射性血管交感神経活動検出プログラム。
- 生体動脈の脈波データと、前記脈波データに対応する拍間隔とに基づいて、圧反射機能に伴う血管の交感神経活動である圧反射性血管交感神経活動を検出する圧反射性血管交感神経活動検出方法であって、
前記拍間隔がｎ（ｎは３以上の自然数）拍連続して増加または減少する系列のうち、前記拍間隔と前記脈波データとの相関係数が、第ｎ－１拍目までは任意の正の閾値より大きく、かつ、第ｎ拍目に前記相関係数が前記閾値以下となる系列を前記圧反射性血管交感神経活動を示す圧反射性血管交感神経活動系列（ＶＢＲＳＡ系列）として検出するＶＢＲＳＡ系列検出ステップを有する、圧反射性血管交感神経活動測定方法。 - 前記圧反射性血管交感神経活動系列における、ｎ拍分の前記拍間隔とｎ－１拍分の前記脈波データとの回帰分析によって、第ｎ拍目の前記脈波データの理論値を算出する理論値算出ステップと、
第ｎ拍目の前記脈波データの実測値と前記理論値とに基づいて、前記圧反射性血管交感神経活動を示す指標であるＶＢＲＳＡ値を算出するＶＢＲＳＡ値算出ステップとを有する、請求項１３に記載の圧反射性血管交感神経活動検出方法。 - 前記ＶＢＲＳＡ値算出ステップでは、下記式（１）～（３）のいずれかを用いて前記ＶＢＲＳＡ値が算出される、請求項１４に記載の圧反射性血管交感神経活動検出方法；
ＶＢＲＳＡ値＝｜ＰＶｔ－ＰＶｍ｜ ・・・式（１）
ＶＢＲＳＡ値＝｜ｌｎ（ＰＶｔ）－ｌｎ（ＰＶｍ）｜ ・・・式（２）
ＶＢＲＳＡ値＝｜ＰＶｔ／ＰＶｍ｜ ・・・式（３）
但し、各符号は以下を表す。
ＰＶｔ：脈波データの理論値
ＰＶｍ：脈波データの実測値 - 前記ＶＢＲＳＡ値算出ステップでは、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して増加する上昇系列と、前記ＶＢＲＳＡ系列のうち前記拍間隔がｎ拍連続して減少する下降系列とのそれぞれについて、前記ＶＢＲＳＡ値を個別に算出する、請求項１４または請求項１５に記載の圧反射性血管交感神経活動検出方法。
- 前記ＶＢＲＳＡ値算出ステップでは、前記上昇系列について算出された前記ＶＢＲＳＡ値の平均値と、前記下降系列について算出された前記ＶＢＲＳＡ値の平均値とを合算または平均化した値を前記ＶＢＲＳＡ値として算出する、請求項１６に記載の圧反射性血管交感神経活動検出方法。
- 前記脈波データは、光電式容積脈波、規準化された前記光電式容積脈波である規準化脈波容積、または対数化された前記規準化脈波容積のうちのいずれかである、請求項１３から請求項１７のいずれかに記載の圧反射性血管交感神経活動検出方法。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017067476 | 2017-03-30 | ||
| JP2017067476 | 2017-03-30 | ||
| PCT/JP2018/013525 WO2018181851A1 (ja) | 2017-03-30 | 2018-03-30 | 圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2018181851A1 JPWO2018181851A1 (ja) | 2020-02-20 |
| JP7081831B2 true JP7081831B2 (ja) | 2022-06-07 |
Family
ID=63676089
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019510214A Active JP7081831B2 (ja) | 2017-03-30 | 2018-03-30 | 圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US11517207B2 (ja) |
| EP (1) | EP3603495B1 (ja) |
| JP (1) | JP7081831B2 (ja) |
| WO (1) | WO2018181851A1 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA3116587A1 (en) * | 2018-10-24 | 2020-04-30 | Research Triangle Institute | Systems, methods, and devices for determining mild traumatic brain injuries |
| US20240261179A1 (en) * | 2023-02-03 | 2024-08-08 | Stryker Corporation | Enhancing patient monitoring with multiple sensors |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008086568A (ja) | 2006-10-02 | 2008-04-17 | Fukuda Denshi Co Ltd | 血圧反射機能測定装置 |
| JP2013202123A (ja) | 2012-03-27 | 2013-10-07 | Tohoku Univ | 自律神経機能測定装置 |
| WO2014157605A1 (ja) | 2013-03-29 | 2014-10-02 | 北海道公立大学法人札幌医科大学 | 神経性圧反射感度測定装置、神経性圧反射感度測定プログラムおよび神経性圧反射感度測定方法 |
-
2018
- 2018-03-30 EP EP18777468.2A patent/EP3603495B1/en active Active
- 2018-03-30 US US16/498,839 patent/US11517207B2/en active Active
- 2018-03-30 WO PCT/JP2018/013525 patent/WO2018181851A1/ja unknown
- 2018-03-30 JP JP2019510214A patent/JP7081831B2/ja active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008086568A (ja) | 2006-10-02 | 2008-04-17 | Fukuda Denshi Co Ltd | 血圧反射機能測定装置 |
| JP2013202123A (ja) | 2012-03-27 | 2013-10-07 | Tohoku Univ | 自律神経機能測定装置 |
| WO2014157605A1 (ja) | 2013-03-29 | 2014-10-02 | 北海道公立大学法人札幌医科大学 | 神経性圧反射感度測定装置、神経性圧反射感度測定プログラムおよび神経性圧反射感度測定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP3603495B1 (en) | 2024-10-23 |
| EP3603495A1 (en) | 2020-02-05 |
| WO2018181851A1 (ja) | 2018-10-04 |
| JPWO2018181851A1 (ja) | 2020-02-20 |
| EP3603495A4 (en) | 2021-01-06 |
| US11517207B2 (en) | 2022-12-06 |
| US20200113455A1 (en) | 2020-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11160459B2 (en) | Monitoring health status of people suffering from chronic diseases | |
| US10405791B2 (en) | Method and continuously wearable noninvasive apparatus for automatically detecting a stroke and other abnormal health conditions | |
| JP6130474B2 (ja) | 体重計装置及び脈波速度取得方法 | |
| RU2673379C2 (ru) | Система мониторинга и способ для мониторинга гемодинамического статуса субъекта | |
| JP5408751B2 (ja) | 自律神経機能測定装置 | |
| US20180199824A1 (en) | Device for measuring biological signals | |
| JP2011521702A (ja) | Ｃｏ２評価のための方法および装置 | |
| WO2013165474A1 (en) | Continuously wearable non-invasive apparatus for detecting abnormal health conditions | |
| JP6203737B2 (ja) | 指細小動脈拡張能検査方法、指細小動脈拡張能検査装置および指細小動脈拡張能検査プログラム | |
| Byfield et al. | Towards robust blood pressure estimation from pulse wave velocity measured by photoplethysmography sensors | |
| JP7081831B2 (ja) | 圧反射性血管交感神経活動検出装置、圧反射性血管交感神経活動検出プログラムおよび圧反射性血管交感神経活動検出方法 | |
| Foo et al. | Pulse transit time changes observed with different limb positions | |
| JP6297539B2 (ja) | 神経性圧反射感度測定装置、神経性圧反射感度測定プログラムおよび神経性圧反射感度測定方法 | |
| Foo et al. | Dual-channel photoplethysmography to monitor local changes in vascular stiffness | |
| JP6850493B2 (ja) | 指細小動脈拡張能連続検査プログラム、指細小動脈拡張能連続検査装置および指細小動脈拡張能連続検査方法 | |
| JP4453431B2 (ja) | 生体情報測定装置 | |
| Centonze | A Novel Approach For Blood Pressure Prediction Using Machine Learning Techniques | |
| US20230148884A1 (en) | Method and device for determining volemic status and vascular tone | |
| Seifer et al. | Monitoring devices for falls and syncope | |
| Huotari et al. | Merged Poincare plot based heart rate analysis | |
| Dvorak et al. | Pulse transition time in heamodynamic measurement: Methods, problems, hardware equipment and first results | |
| EP3581104A1 (en) | Method, device and computer program product for estimating a compliance of a blood vessel in a subject | |
| CN116439681A (zh) | 一种智能动态血压检测仪及血压检测方法 | |
| Etemadi et al. | Emerging Wearable Physiological Monitoring Technologies & Decision Aids for Health & Performance: Wearable ballistocardiogram and seismocardiogram systems for health and performance | |
| Dalla Pozza et al. | Calculating sympathovagal balance from heart rate variability: are there alternatives in adolescents? |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20190917 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210329 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20220427 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220519 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7081831 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |