▶ 株式会社東芝の特許一覧 ▶ 東芝デバイス＆ストレージ株式会社の特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-12-15
(45)【発行日】2023-12-25
(54)【発明の名称】血圧測定装置、及び血圧測定方法
(51)【国際特許分類】
A61B 5/021 20060101AFI20231218BHJP
A61B 5/02 20060101ALI20231218BHJP
A61B 5/022 20060101ALI20231218BHJP
【ＦＩ】
A61B5/021
A61B5/02 310A
A61B5/022 400F
【請求項の数】
4
(21)【出願番号】P 2020026521
(22)【出願日】2020-02-19
(65)【公開番号】P2020189077
(43)【公開日】2020-11-26
【審査請求日】2022-03-14
(31)【優先権主張番号】P 2019035280
(32)【優先日】2019-02-28
(33)【優先権主張国・地域又は機関】JP
(31)【優先権主張番号】P 2019093999
(32)【優先日】2019-05-17
(33)【優先権主張国・地域又は機関】JP
(73)【特許権者】
【識別番号】000003078
【氏名又は名称】株式会社東芝
(73)【特許権者】
【識別番号】317011920
【氏名又は名称】東芝デバイス＆ストレージ株式会社
(74)【代理人】
【識別番号】100091487
【弁理士】
【氏名又は名称】中村 行孝
(74)【代理人】
【識別番号】100105153
【弁理士】
【氏名又は名称】朝倉 悟
(74)【代理人】
【識別番号】100107582
【弁理士】
【氏名又は名称】関根 毅
(74)【代理人】
【識別番号】100118843
【弁理士】
【氏名又は名称】赤岡 明
(74)【代理人】
【識別番号】100125151
【弁理士】
【氏名又は名称】新畠 弘之
(72)【発明者】
【氏名】川上 健
【審査官】蔵田 真彦
(56)【参考文献】
【文献】特開２００４－１５４２３１（ＪＰ，Ａ）
【文献】特開２０１６－１４６９５８（ＪＰ，Ａ）
【文献】特開２０１６－１１２２７７（ＪＰ，Ａ）
【文献】特開平１０－２９５６５７（ＪＰ，Ａ）
【文献】韓国登録実用新案第２０－０３８０９２８（ＫＲ，Ｙ１）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／０２－５／０３
(57)【特許請求の範囲】
【請求項1】
所定の周波数帯域の光信号を照射したときに被検者の体内で散乱された受光信号に基づいて、前記被検者の脈波を測定する測定部と、前記被検者の血圧を取得する血圧取得部と、を備え、
前記血圧取得部は、
前記測定された脈波の第１拍において、前記脈波を時間で１階微分した値が最大になる第１基準時の前記脈波の値から直流成分を減じた第１差分値と、前記脈波の最大ピークの時刻である第４基準時の脈波の値から前記脈波の直流成分を減じた第２差分値と、前記脈波の直流成分と、前記第４基準時から前記第１基準時を減じた第１差分時間と、前記脈波の振幅が前記第４基準時の位置から単調に減少した後のボトムの時刻である第２基準時から、前記第４基準時を減じた第３差分時間と、を生成する特徴点処理部と、
前記脈波の直流成分と、前記第１差分値と、前記第２差分値と、前記第１差分時間と、前記第３差分時間と、を少なくとも用いて拡張期血圧を演算する血圧演算部と、を有し、
前記血圧演算部は、前記第１差分値と、前記第２差分値と、前記第１差分時間と、前記脈波の直流成分と、前記第３差分時間とに、少なくとも基づく第１所定値と、
前記第１差分値と、前記脈波の直流成分との差分を前記第１基準時の１階微分値で除算した第２所定値と、に基づき、前記拡張期血圧を演算する、血圧測定装置。
【請求項2】
前記特徴点処理部は、前記第１差分値を前記第１基準時の１階微分値で除算した値に基づく第３基準時と前記第１基準時との第２差分時間を更に生成し、前記血圧演算部は、前記拡張期血圧の値と、前記第２所定値と、前記直流成分、前記第１差分値、及び前記第２差分時間に基づく第３値と、に基づき、収縮期血圧を演算する、請求項１に記載の血圧測定装置。
【請求項3】
前記特徴点処理部は、前記測定された脈波の第１拍は、数拍分のデータを平均化して演算する、請求項１又は２に記載の血圧測定装置。
【請求項4】
被検者に所定の周波数帯域の光信号を照射したときに被検者の体内で散乱されて受光された受光信号に基づいて、前記被検者の脈波を測定する測定工程と、拡張期血圧を取得する血圧取得工程と、を備え、
前記血圧取得工程は、
前記測定された脈波の第１拍において、前記脈波を時間で１階微分した値が最大になる第１基準時の前記脈波の値から直流成分を減じた第１差分値と、前記脈波の最大ピークの時刻である第４基準時の脈波の値から前記脈波の直流成分を減じた第２差分値と、前記脈波の直流成分と、前記第４基準時から前記第１基準時を減じた第１差分時間と、前記脈波の振幅が前記第４基準時の位置から単調に減少した後のボトムの時刻である第２基準時から、前記第４基準時を減じた第３差分時間と、を生成する特徴点処理工程と、
前記脈波の直流成分と、前記第１差分値と、前記第２差分値と、前記第１差分時間と、前記第３差分時間と、を少なくとも用いて拡張期血圧を演算する血圧演算工程と、
を有し、
前記血圧演算工程は、前記第１差分値と、前記第２差分値と、前記第１差分時間と、前記脈波の直流成分と、前記第３差分時間とに、少なくとも基づく第１所定値と、
前記第１差分値と、前記脈波の直流成分との差分を前記第１基準時の１階微分値で除算した第２所定値と、に基づき、前記拡張期血圧を演算する、血圧測定方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明の実施形態は、血圧測定装置、及び血圧測定方法に関する。
【背景技術】
【0002】
心拍数の変化に対応する動脈及び毛細血管の血液量の変化を測定することにより、心拍に伴う脈波を検出する光電式容積脈波計（ＰＰＧ：Ｐｈｏｔｏｐｌｅｔｈｙｓｍｏｇｒａｍ）センサが知られている。ＰＰＧセンサを用いて脈拍ごとに組織を通過する血液量に基づいて心拍数を検出する手法は、容積脈波（ＢＶＰ：Ｂｌｏｏｄ Ｖｏｌｕｍｅ Ｐｕｌｓｅ）測定と呼ばれている。
【0003】
容積脈波形状の特徴点に基づき血圧の推定を行う方法が一般に知られている。ところが、容積脈波の波形は、被検者の活動状態や精神状態によって変動し、容積脈波に乱れが生じる。容積脈波が乱れていると、特徴点等を正確に計測できなくなり、血圧の測定精度が悪くなる。
【先行技術文献】
【特許文献】
【0004】
【文献】特開２０００－２１７７９６号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
本発明が解決しようとする課題は、被検者の血圧を簡易かつ精度よく取得可能な血圧測定装置及び血圧測定方法を提供することである。
【課題を解決するための手段】
【0006】
本実施形態によれば、血圧測定装置は、測定部と、血圧取得部とを備える。測定部は、所定の周波数帯域の光信号を照射したときに被検者の体内で散乱された受光信号に基づいて、被検者の脈波を測定する。血圧取得部は、脈波を時間で１階微分した値が最大になる第１基準時から次の脈波の立ち上がる第２基準時までの期間内における第１期間の前記被検者の血液流量に対応する第１値と、被検者の血流抵抗に対応する第２値とに基づき、拡張期血圧を取得する。
【図面の簡単な説明】
【0007】
【図1】第１実施形態による血圧測定装置の概略構成を示すブロック図。
【図2】腕時計型の血圧測定装置の一例を示す図。
【図3】測定部により測定された容積脈波の一例を示す図。
【図4】円筒管で近似した血管を示す図。
【図5】第１実施形態に係る脈波の波形の一例とモデル化した円筒管の半径を示す図。
【図6】第１実施形態に係る血管の容積変化に伴う円筒管の半径の変化を模式的に示す図。
【図7A】点ＰＳから点ＰＥへ至る際の流量を模式的に示す図。
【図7B】点ＰＥから点ＰＨへ至る際の流量を模式的に示す図。
【図8】第１実施形態に係る特徴点処理部が取得する情報例を示す図。
【図9A】血圧が高めの被検者の測定データを示す図。
【図9B】血圧が低めの被検者の測定データを示す図。
【図10】血圧測定装置の処理を示すフローチャート。
【図11】第２実施形態に係る脈波の波形の一例とモデル化した円筒管の半径との関係を示す図。
【図12】第２実施形態に係る血管の容積変化に伴う円筒管の半径の変化を模式的に示す図。
【図13A】点ＰＳから点ＰＥを経て点ＰＨへ至る際の流量を模式的に示す図。
【図13B】点ＰＨから点ＰＤへ至る際の流量を模式的に示す図。
【図14】第２実施形態に係る特徴点処理部が取得する情報例を示す図。
【発明を実施するための形態】
【0008】
以下、図面を参照して本発明の実施形態を説明する。以下の実施形態では、血圧測定装置内の特徴的な構成および動作を中心に説明するが、血圧測定装置には以下の説明で省略した構成および動作が存在しうる。
【0009】
（第１実施形態）
図１は第１実施形態による血圧測定装置１の概略構成を示すブロック図である。血圧測定装置１は、測定部２と、血圧取得部４とを備えている。血圧測定装置１は、例えば図２に示すような腕時計型の生体測定装置６に組み込むことができる。なお、生体測定装置６は、上腕部、胸部などに配置しても良い。
図１は第１実施形態による血圧測定装置１の概略構成を示すブロック図である。血圧測定装置１は、測定部２と、血圧取得部４とを備えている。血圧測定装置１は、例えば図２に示すような腕時計型の生体測定装置６に組み込むことができる。なお、生体測定装置６は、上腕部、胸部などに配置しても良い。
【0010】
測定部２は、被検者の心拍数の変化に伴う動脈及び毛細血管の血液量の変化を測定することにより、心拍に伴う容積脈波の情報を取得する。以下では、容積脈波を単に脈波と呼ぶこともある。
【0011】
測定部２は、発光部２２と、受光部２４と、脈波生成部２６とを有する。発光部２２は、例えば、特定の波長帯域（緑色や近赤外帯域など）の光信号を発光するＬＥＤ（Ｌｉｇｈｔ Ｅｍｉｔｔｉｎｇ Ｄｅｖｉｃｅ）を有する。受光部２４は、発光部２２からの光信号が被検者の体内で吸収又は反射・散乱された後の信号を受光する。脈波生成部２６は、受光信号に基づいて、心拍の１拍ごとに脈波を生成する。
【0012】
光信号の発光量が変動すると、受光信号の受光量も変動する。このため、脈波生成部２６は、受光信号のＤＣ成分とＡＣ成分とに分離し、ＡＣ／ＤＣ比に基づいて脈波を生成する。このため、生成される脈波は、無次元のデータである。
【0013】
図３は、測定部２により測定された容積脈波の一例を示す図である。縦軸は脈波の値を示し、横軸は時間を示している。図３に示すように、脈波は一心拍毎に変動を繰り返す。ｉ拍目の脈波ｙｉは周期的に変動するＡＣ成分とＤＣ成分である
により構成される。
【0014】
血圧取得部４は、脈波に基づき、被検者の血圧を取得する。この血圧取得部４は、特徴点処理部４２と、血圧演算部４４とを有する。
【0015】
まず、図４乃至７Ｂに基づき、血圧取得部４で用いるモデルに関して説明する。図４は、モデル化された血管を示す図である。図４に示す血管のモデルは、半径ｒｉｓ、長さＬの円筒管で近似したものである。血圧変動は、心臓から拍出された血液により血管壁にかかる圧力の変動である。この血圧変動は、脈波ｙｉと連動する。
【0016】
円筒管の圧力差ΔＰ、流量Ｑ、および抵抗Ｒの関係は、ナビエ・ストークスの式から導かれ、（１）式で示される。
血圧取得部４は、円筒管のモデルに基づき、脈波ｙｉを用いて流量Ｑおよび抵抗Ｒに対応する値を演算し、被検者の血圧を取得する。なお、人の血圧は一般に、心臓の収縮期における血管内の最大圧力である収縮期血圧ＳＢＰ（Ｓｙｓｔｏｌｉｃ Ｂｌｏｏｄ Ｐｒｅｓｓｕｒｅ）、心臓の拡張期における血管内の最小圧力である拡張期血圧ＤＢＰ（Ｄｉａｓｔｏｌｉｃ Ｂｌｏｏｄ Ｐｒｅｓｓｕｒｅ）、拡張期血圧から収縮期血圧を減じた脈圧ＰＰ（Ｐｕｌｓｅ ｐｒｅｓｓｕｒｅ）を用いて評価される。
【数1】
【0017】
図５は、第１実施形態に係る脈波の波形の一例とモデル化した円筒管の半径との関係を示す図である。左図は、１拍分の正常な脈波の波形の一例を示す図である。縦軸は脈波の値を示し、横軸は時間を示している。右図は、モデル化した円筒管の半径を示す図である。血管の容積変化を、半径ｒｉｓと変化分であるΔｒｉｄで示している。すなわち、半径ｒｉｓは、点ＰＥにおける血管の半径であり、Δｒｉｄは、点ＰＥから点ＰＨへの容積の増加に係る半径の増分である。
【0018】
正常な脈波ｙｉは、振幅がボトムの位置（ｔ０）で開始し、振幅がほぼ単調に増加して最大ピーク（ｔ２）に達し、その後、振幅が単調に減少してボトムの位値（ｔ３）に達して終了する。ここで、添え字ｉは、容積脈波データにおいて、各一つのパルスを識別する番号である。すなわち、ｉ拍目の脈波に対応するデータを意味する。なお、本実施形態に係る演算では、一拍毎の演算を行っているが、これに限定されず、数拍分のデータを例えば平均化して演算してもよい。
【0019】
ｔ１は、ｔ０からｔ２までの間で、脈波ｙｉを時間で１階微分した値が最大になる時刻である。このｔ１は、後述する（４）式で示す粘弾性運動方程式の変位ｒ（ｔ）の平衡点に対応する。
【0020】
ｔｓは、（２）式で示すように、時刻ｔ１の脈波ｙｉの値ｙｉｓから直流成分
を減じた第１差分値を、最大になる１階微分値
で除算した値に基づく時刻である。ここで、１階微分値
は、例えば後述する（１５）式で演算される。
線Ｌｓは、点ＰＥの接線である。すなわち、線Ｌｓと横軸と水平なラインとの間の角度θにより演算されるｔａｎθは、この１階微分値
に対応する。ＰＳ、ＰＥ、ＰＨ、ＰＬは、それぞれ時刻ｔｓ、ｔ１、ｔ２、ｔ３に対応する点である。なお、本実施形態に係る時刻ｔ１が第１時刻に対応し、時刻ｔ２が第２時刻に対応し、時刻ｔ０又はｔｓが第３時刻に対応する。
【数2】
【0021】
図６は、点ＰＳから点ＰＬへの血管の容積変化に伴う円筒管の半径の変化を模式的に示す図である。縦軸は時間を示し、横軸は点ＰＳからの半径の変動分を示している。半径は時間の経過にしたがい、点ＰＳから点ＰＨまでは増加し、その後は減少する。
【0022】
血液流量を、モデル化された血管の容積変化に応じた半径ｒに基づく体積流量Ｑとして計量する。このため、単位時間あたりの半径の平均変化率（Δｒ／Δｔ）を用いて流量Ｑを定義できる。Δｔは、ｔの時間変化量を示し、Δｒは、Δｔにおける半径ｒの変化量を示す。
【0023】
図７Ａは、点ＰＳから点ＰＥへ至る際の流量ＱＳＥを模式的に示す図である。図７Ｂは、点ＰＥから点ＰＨへ至る際の流量ＱＥＨを模式的に示す図である。横軸は平均変化率（Δｒ／Δｔ）の二乗を示し、縦軸は長さＬとπの乗算値である。平均変化率ｍｉｓは、後述する（８）式で示される。図７Ａの平均変化率ｍｉｓは、点ＰＳから点ＰＥへ至る際の半径ｒｉｓを時刻ｔ１から時刻ｔｓを減算した値で除算した値である。図７Ｂの平均変化率ｍｉｄは、Δｒｉｄを時刻ｔ２から時刻ｔ１を減算した値で除算した値である。また、点ＰＨから点ＰＬへ至る際流量ＱＥＬ（図示せず）、は、点ＰＨから点ＰＬに至る時間Ｔｉｄ２と抵抗Ｒと血管のコンプライアンスＣと流量ＱＥＨを用いて、（３）式で得ることが可能である。（３）式において、ネイピア数の項は、２要素のウィンドケッセルモデルにおいて、収縮期後の圧力低下を表す方法として知られている。
【数3】
【0024】
血管壁の変位ｒ（ｔ）、すなわち半径ｒ（ｔ）の変位は脈波ｙｉの値に連動する。また、圧力を血管壁の変位ｒ（ｔ）で近似する事が可能であり、半径ｒ（ｔ）の変位は（４）式で示す粘弾性運動方程式と同等となる。すなわち、血管を円筒管で近似する場合、血圧を脈波ｙｉの情報に基づき演算可能となる。
【数4】
【0025】
ここで、ｂは粘性定数であり、ｋは弾性定数であり、血管の弾性を反映する。（４）式において左辺は、ニュートンの運動の第２法則における全体の力を示す。また、右辺の第１項は、減衰力を示し、第２項は、復元力を示し、第３項は、ウィンドケッセル効果による力を示している。変位ｒ（ｔ）の平衡位置は点ＰＥに対応する。
【0026】
心臓の収縮期間中、脈波ｙｉの立ち上がりから点ＰＥまでは、主に、ウィンドケッセル効果によって血管壁が変位する。一方で、点ＰＥから点ＰＨまでは、減衰力および復元力によって変位する。なお、点ＰＨでは、ウィンドケッセル効果および減衰力は無視できる。
【0027】
このため、心臓の収縮期間中、点ＰＥまでの半径ｒｉｓの拡張は、主にウィンドケッセル効果で生じる。従って、心臓の収縮期に生じる力、すなわち、収縮期血圧ＳＢＰは、流量ＱＳＥに反映されている。一方で、拡張期血圧ＤＢＰは、流量ＱＥＬに反映されている。流量ＱＥＬは、ウィンドケッセル効果による力の下限であるから、点ＰＥから点ＰＨへ至る際の流量ＱＥＨから得られる。なお、本実施形態に係る流量ＱＥＨが第１値に対応し、抵抗Ｒが第２値に対応し、流量ＱＳＥが第３値に対応する。
【0028】
そこで、本実施形態では、心拍ｉにおける収縮期血圧ＳＢＰｉを（５）式でモデル化する。Ｒｉは、流量ＱＳＥｉおよびＱＥＬｉの観測の際における末梢の循環抵抗を反映する値である。ここでａ、αは定数である。
【数5】
【0029】
また、拡張期血圧ＤＢＰｉを（６）式でモデル化する。
ここでｂ、βは定数である。定数ａ、α、ｂ、βなどは、血圧測定装置１での測定値と、血圧測定装置１と同時に測定した医療機器（例えば手首カフ式）で測定したデータが一致するように例えば最小二乗法などで演算することが可能である。一度、定数ａ、α、ｂ、βを決定すると、後の測定では、定数ａ、α、ｂ、βの演算は不要である。なお、脈圧ＰＰｉは、収縮期血圧ＳＢＰｉから拡張期血圧ＤＢＰｉを減算した値である。
【数6】
【0030】
以上が本実施形態に係る血圧取得部４が用いるモデルの説明であるが、以下に血圧取得部４の詳細な構成を説明する。図８は、脈波ｙｉから特徴点処理部４２が取得する情報例を示す図である。縦軸は脈波の値を示し、横軸は時間を示している。
【0031】
特徴点処理部４２は、ｔ０を立ち上がり時刻として検出し、ｔ２を最大ピークになる時刻として検出する。また、特徴点処理部４２は、ｔ０からｔ２までの間で、脈波を時間で１階微分した値が最大になるｔ１を算出する。
【0032】
また、特徴点処理部４２は、時刻ｔ１の脈波ｙｉの値から直流成分
を減じた第１差分値Δｙｉｓと、時刻ｔ２の脈波ｙｉの値から直流成分
を減じた第２差分値Δｙｉｈを演算する。
【0033】
特徴点処理部４２は、（７）式を用いて時間Ｔｉｓを演算する。Ｔｉｓは、時刻ｔ１の微分値に対応するｔａｎθで第１差分値Δｙｉｓを除算した時間である。そして、ｔ１からＴｉｓを減じた第２立ち上がり時刻ｔｓを演算する。Ｔｉｄ１は、時刻ｔ２からｔ１を減じた時間であり、Ｔｉｄ２は、時刻ｔ３からｔ４を減じた時間である。
立ち上がり時刻ｔ０における脈波ｙｉの形状は個人差が大きく、なだらかに変化する人と急峻に変化する人とが存在する。このため、立ち上がり時刻ｔ０と時刻ｔ１との差分値は、個人差により変動しやすくなる。これに対し、第２立ち上がり時刻ｔｓと時刻ｔ１との時間差Ｔｉｓは、個人差による変動が低減され、安定した値を示す。このため、流量の演算では、この時間差Ｔｉｓを用いる。なお、脈波ｙｉの形状によっては、時間Ｔｉｓを時刻ｔ０と時刻ｔ１との時間差として演算してもよい。これにより、計算を簡易化できる。
【数7】
【0034】
血圧演算部４４は、（３）、（６）式で示したように、脈波ｙｉを時間で１階微分した値が最大になる第１時刻ｔ１から脈波の最大ピークの第２時刻ｔ２までの被検者の血液流量に対応する流量ＱＥＨ（第１値）と、被検者の血流抵抗に対応する抵抗Ｒ（第２値）とに基づき、拡張期血圧ＤＢＰを取得する。すなわち、血圧演算部４４は、拡張期血圧ＤＢＰとして、流量ＱＥＨに基づく流量ＱＥＬと抵抗Ｒの乗算に所定係数ｂを乗算し、所定の定数βを更に加算した値を演算する。なお、抵抗Ｒは、後述する（１４）、（１５）式に基づき演算される。また、この血圧演算部４４は、（５）式を用いて脈波の立ち上がり第３時刻ｔｓ又はｔ０から時刻ｔ１までの血液流量に対応する流量ＱＳＥ（第３値）に更に基づき、収縮期血圧ＳＢＰを取得する。なお、本実施形態に係る第１時刻が第１基準時に対応し、第２時刻が第４基準時に対応し、第３時刻が第３基準時に対応する。
【0035】
より詳細には、血圧演算部４４は、ｒｉｓを（８）式に基づき演算し、平均変化率ｍｉｓを（９）式に基づき演算する。流量ＱＳＥを（１０）式に基づき演算する。ここで、
は、点ＰＥの血管容積に比例する。このように、血圧演算部４４は、（５）式及び（１０）式に基づき、収縮期血圧ＳＢＰｉを取得する。この際に演算する脈波ｙｉの第１差分値Δｙｉｓ、直流成分
、時間Ｔｉｓは、脈波ｙｉの変動に対しても安定的に且つ簡易に演算可能である。このため、収縮期血圧ＳＢＰｉを簡易かつ精度よく取得できる。Ｇは定数である。
【数8】
【数9】
【数10】
【0036】
血圧演算部４４は、Δｒｉｄを（１１）式に基づき演算し、平均変化率ｍｉｄを（１２）式に基づき演算する。そして、血圧演算部４４は、更に、点ＰＨから点ＰＬへ至る際の流量ＱＥＬを（２）式に基づき演算する。このように、血圧演算部４４は、（６）式及び（１３）式に基づき、拡張期血圧ＤＢＰｉを取得する。この際に演算する第１差分値Δｙｉｓ、第２差分値Δｙｉｈ、直流成分
、時間Ｔｉｓ、Ｔｉｄは、脈波ｙｉの変動に対しても安定的に且つ簡易に演算可能である。このため、収縮期血圧ＳＢＰｉを簡易かつ精度よく取得できる。
【数11】
【数12】
【数13】
【0037】
このように、第１値ＱＥＬ及び第３値ＱＳＥは、時刻ｔｓにおける被検者の血管容積に対応する
に基づく値である。また、血管容積は、第１差分値Δｙｉｓを直流成分
で除算した値に基づく値である。すなわち、第１値ＱＥＬは、第２差分値Δｙｉｈから第１差分値Δｙｉｓを減じて第１差分値Δｙｉｓで除算した値と、血管容積に対応する
の平方根と、の乗算に基づく値である。
【0038】
血圧演算部４４は、（１４）、（１５）式に基づき、被検者の血流抵抗に対応するＲｉを演算する。
ここで、ｆｓは脈波ｙｉのサンプリング周波数である。
【数14】
【数15】
【0039】
図９Ａは、血圧が高めの被検者の測定データを示す図である。図９Ｂは、血圧が低めの被検者の測定データを示す図である。縦軸は血圧を示し、横軸は時間を示す。また、菱形のマークが血圧測定装置１での測定値を示し、実線は比較対象用に測定した医療機器（手首カフ式）のデータを示す。本実施形態の血圧測定装置１で測定した値は、どちらの場合にも、比較対象用に測定したデータとよく一致する。
【0040】
図１０は血圧測定装置１の処理を示すフローチャートである。まず、測定部２は、被検者の脈波を取得する（ステップＳ１００）。続けて、特徴点処理部４２は、脈波ｙｉに基づき処理を行う。
【0041】
次に血圧演算部４４は、脈波ｙｉを時間で１階微分した値が最大になる時刻から脈波ｙｉの最大ピークの時刻までの脈波ｙｉ被検者の血液流量に対応する流量ＱＥＬ、被検者の血流抵抗に対応する抵抗Ｒ、および脈波の立ち上がり時刻から脈波を時間で１階微分した値が最大になる時刻までの血液流量に対応する流量ＱＳＥを演算する（ステップＳ１０２）。
【0042】
次に血圧演算部４４は、流量ＱＥＬ、抵抗Ｒ、流量ＱＳＥに基づき、拡張期血圧ＤＢＰ、収縮期血圧ＳＢＰ、脈圧ＰＰを算出する。（ステップＳ１０４）。血圧演算部４４は、全体処理を終了するか否かを判定し（ステップＳ１０６）、全体処理を終了する場合（ステップＳ１０６：ＹＥＳ）、全体処理を終了し、全体処理を終了しない場合（ステップＳ１０６：ＮＯ）、ステップＳ１００からの処理を繰り返す。
【0043】
このように、本実施形態では、（５）式及び（１０）式に基づき、収縮期血圧ＳＢＰｉを取得し、（６）式及び（１３）式に基づき、拡張期血圧ＤＢＰｉを取得するため、血圧を簡易かつ精度よく検出できる。
【0044】
（第２実施形態）
血圧測定装置１は、第１実施形態では収縮期血圧ＳＢＰを流量ＱＳＥ（図６）に基づき、演算していたのに対し、第２実施形態では収縮期血圧ＳＢＰを流量ＱＥＨにも基づき、演算する点で相違する。また、血圧測定装置１は、第１実施形態では拡張期血圧ＤＢＰを流量ＱＥＨ（図６）に基づき、演算していたいのに対し、第２実施形態では、拡張期血圧ＤＢＰを流量ＱＨＤに基づき、演算する点で相違する。以下では第１実施形態と異なる点に関して説明する。
血圧測定装置１は、第１実施形態では収縮期血圧ＳＢＰを流量ＱＳＥ（図６）に基づき、演算していたのに対し、第２実施形態では収縮期血圧ＳＢＰを流量ＱＥＨにも基づき、演算する点で相違する。また、血圧測定装置１は、第１実施形態では拡張期血圧ＤＢＰを流量ＱＥＨ（図６）に基づき、演算していたいのに対し、第２実施形態では、拡張期血圧ＤＢＰを流量ＱＨＤに基づき、演算する点で相違する。以下では第１実施形態と異なる点に関して説明する。
【0045】
第１実施形態に係る血圧測定装置１で測定した血圧は、一般的な人の収縮期血圧ＳＢＰ、及び拡張期血圧ＤＢＰと良く一致する。しかし、被検者の中には、脈波の特性が一般の人と異なる特性を有する人が存在することが分かってきた。第２実施形態に係る血圧測定装置１では、このような人にも対応可能とする。
【0046】
図１１は、第２実施形態に係る脈波の波形の一例とモデル化した円筒管の半径との関係を示す図である。図５と同様に、左図は、１拍分の正常な脈波の波形の一例を示す図である。縦軸は脈波の値を示し、横軸は時間を示している。右図は、モデル化した円筒管の半径を示す図である。血管の容積変化を、半径ｒｓｉと変化分であるΔｒｄｉで示している。点ＰＤは、点ＰＨと点ＰＬの間において点ＰＥと同じ容積脈波の値を示す点である。Ｉｄｃは、容積脈波の直流成分である。ここで、添え字ｉは、容積脈波データにおいて、各一つのパルスを識別する番号である。すなわち、ｉ拍目の脈波に対応するデータを意味する。なお、点ＰＤの時刻が第５基準時に対応する。
【0047】
図１２は、点ＰＳから点ＰＬへの血管の容積変化に伴う円筒管の半径の変化を模式的に示す図である。すなわち、図１２は一拍分の脈波に関する円筒管の半径の変化を示している。縦軸は時間を示し、横軸は点ＰＳからの半径の変動分を示している。半径は時間の経過にしたがい、点ＰＳから点ＰＨまでは増加し、その後は減少する。
【0048】
図１３Ａは、点ＰＳから点ＰＥを経て点ＰＨへ至る際の流量ＱＳを模式的に示す図である。図１３Ｂは、点ＰＨから点ＰＤへ至る際の流量ＱＨＤを模式的に示す図である。横軸は血管半径ｒの平均変化率（Δｒ／Δｔ）の二乗を示し、縦軸は長さＬとπの乗算値である。図１３Ａのｍｓｉは、点ＰＳから点ＰＥへ至る際の半径の平均変化率であり、ｍｄ１ｉは、点ＰＥから点ＰＨへ至る際の半径の平均変化率である。図１３Ｂのｍｄ２ｉは、点ＰＤから点ＰＨへ至る際の半径の平均変化率である。なお、本実施形態に係る流量ＱＨＤが第１値に対応し、抵抗Ｒが第２値に対応し、流量ＱＳが第３値に対応する。
【0049】
本実施形態では、収縮期血圧ＳＢＰを流量ＱＳにより演算する。心臓の収縮期間中の点ＰＥまでの血管径の拡張は、主にウィンドケッセル効果で生じる。そして、点ＰＥ以後から次第に復元力と減衰力が支配的になる。すなわち、本実施形態では、ウィンドケッセル効果に復元力が加えられた点ＰＥから点ＰＨまでの流量ＱＳＥにまで、収縮期に生じる力の範囲を拡張してモデル化している。被検者によっては、点ＰＥから点ＰＨまでの範囲にもウィンドケッセル効果がより強くでる人がいると考えられる。このような人を含めて測定をする場合には、流量ＱＳを用いると、収縮期血圧ＳＢＰの測定精度がより向上する。なお、流量ＱＳを用いても、一般的な人の収縮期血圧ＳＢＰの精度も低下しないことが実験的に検証されている。
【0050】
一方で、点ＰＥから点ＰＨまでの範囲にウィンドケッセル効果がより強くでる人の場合には、ウィンドケッセル効果が弱まる点は、点ＰＬ側にずれていると考えられる。拡張期血圧は、ウィンドケッセル効果による力の下限であるので、ウィンドケッセル効果が弱まる点を点ＰＨまでずらし、点ＰＨから点ＰＤまでの範囲の流量ＱＤを用いて拡張期血圧ＤＢＰをモデル化している。特に、流量ＱＤを流量ＱＨＤに基づき、演算する。なお、流量ＱＨＤを用いても、一般的な人の拡張期血圧ＤＢＰの精度も低下しないことが実験的に検証されている。
【0051】
以上が本実施形態に係る血圧取得部４が用いるモデルの説明であるが、以下に血圧取得部４の詳細な処理例を説明する。
図１４は、第２実施形態に係る特徴点処理部４２が脈波ｙｉから取得する情報例を示す図である。縦軸は脈波の値を示し、横軸は時間を示している。右図は、モデル化した円筒管の半径を示す図である。血管の容積変化を、半径ｒｓｉと変化分であるΔｒｄｉで示している。
図１４は、第２実施形態に係る特徴点処理部４２が脈波ｙｉから取得する情報例を示す図である。縦軸は脈波の値を示し、横軸は時間を示している。右図は、モデル化した円筒管の半径を示す図である。血管の容積変化を、半径ｒｓｉと変化分であるΔｒｄｉで示している。
【0052】
特徴点処理部４２は、時刻ｔ１の脈波ｙｉの値から直流成分Ｉｄｃを減じた第１差分値Δｙｓｉと、時刻ｔ２の脈波ｙｉの値から直流成分Ｉｄｃを減じた第２差分値Δｙｈｉを演算する。
また、特徴点処理部４２は、（１６）式を用いて時間Ｔｄ２ｉを演算する。Ｔｄ２ｉは、点ＰＤと点ＰＨ間の時間である。Ｔｓｉは、時刻ｔ１からｔ０を減じた時間であり、Ｔｄ１ｉは、時刻ｔ２からｔ１を減じた時間であり、Ｔｄ３ｉは、時刻ｔ３からｔ２を減じた時間である。すなわち、特徴点処理部４２は、脈波ｙｉの時刻ｔ２から時刻ｔ３までの期間内において時刻ｔ１における脈波ｙｉと同等の値を示す点ＰＤの時刻を第５基準時として取得し、時刻ｔ２と第５基準時との間の時間を時間Ｔｄ２ｉとして演算する。
また、特徴点処理部４２は、（１６）式を用いて時間Ｔｄ２ｉを演算する。Ｔｄ２ｉは、点ＰＤと点ＰＨ間の時間である。Ｔｓｉは、時刻ｔ１からｔ０を減じた時間であり、Ｔｄ１ｉは、時刻ｔ２からｔ１を減じた時間であり、Ｔｄ３ｉは、時刻ｔ３からｔ２を減じた時間である。すなわち、特徴点処理部４２は、脈波ｙｉの時刻ｔ２から時刻ｔ３までの期間内において時刻ｔ１における脈波ｙｉと同等の値を示す点ＰＤの時刻を第５基準時として取得し、時刻ｔ２と第５基準時との間の時間を時間Ｔｄ２ｉとして演算する。
【数16】
【0053】
点ＰＬに対応する容積を基準とすれば、Δｙｓｉ／Ｉｄｃは点ＰＥ、ＰＤにおける容積に比例し、同様にΔｙｈｉ／Ｉｄｃは点ＰＨにおける容積に比例する。Ｇは比例定数であり、Ｉｄｃは、脈波のＤＣ成分の値である。
円筒管の半径がｒｓｉからｒｓｉ＋Δｒｄｉに変化した場合、Δｒｄｉを点ＰＥにおける半径ｒｓｉを用いて、（１７）～（１９）式により演算可能である。
ここで、（１９）式の半径ｒｓｉを（１７）式を用いて、整理すると以下の（２０）、（２１）式に変形できる。血圧演算部４４は、（２０）、（２１）式を用いて半径ｒｓｉとΔｒｄｉを演算する。なお、Ｌはモデル円筒管の長さである。
円筒管の半径がｒｓｉからｒｓｉ＋Δｒｄｉに変化した場合、Δｒｄｉを点ＰＥにおける半径ｒｓｉを用いて、（１７）～（１９）式により演算可能である。
【数17】
【数18】
【数19】
【数20】
【数21】
【0054】
血圧演算部４４は、平均変化率ｍｓｉを（２２）式を用いて演算する。
また、血圧演算部４４は、平均変化率ｍｄ１ｉ及びｍｄ２ｉをそれぞれ、（２３）、（２４）式を用いて演算する。
【数22】
【数23】
【数24】
【0055】
血圧演算部４４は、流量ＱＳｉを平均変化率ｍｄ１ｉ及びｍｄ２ｉに基づき、（２５）式を用いて演算する。
【数25】
【0056】
血圧演算部４４は、抵抗Ｒｉを（２６）式を用いて演算する。ここで、Ｖｉは、モデル円筒管の体積であり、Ｖｉ（ｔ１）は、時間ｔ１でのモデル円筒管の体積である。すなわち、Ｉｄｃは、第１実施形態における
に対応する。このため、（２６）式は（１６）式と同等の値となる。
血圧演算部４４は、抵抗ＲｉとコンプライアンスＣを用いて、点ＰＨから点ＰＬへ至る際の流量ＱＤｉを（２７）式に基づき演算する。
抵抗Ｒdｉは、（２８）、（２９）式で得る。
【数26】
【数27】
【0057】
【数28】
【数29】
【数30】
【数31】
【0058】
このように、血圧演算部４４は、脈波ｙｉを時間で１階微分した値が最大になる第１時刻ｔ１から次の脈波の立ち上がる第４時刻ｔ３までの期間内における期間Ｔｄ２ｉ（第１期間）の被検者の血液流量に対応する流量ＱＨＤ（第１値）と、被検者の血流抵抗に対応するＲ（第２値）とに基づき、拡張期血圧ＤＢＰを取得する。また、血圧演算部４４は、脈波の立ち上がる第３時刻ｔ０から脈波の最大ピークの第２時刻ｔ２までの期間内における期間（Ｔｓｉ＋Ｔｄ１ｉ）（第２期間）の被検者の血液流量に対応する流量ＱＳ（第３値）に更に基づき、収縮期血圧を取得する。なお、本実施形態に係る第１時刻が第１基準時に対応し、第２時刻が第４基準時に対応し、第３時刻が第３基準時に対応し、第４時刻が第２基準時に対応する。
【0059】
このように、本実施形態では、（２７）式及び（３０）式に基づき、拡張期血圧ＤＢＰｉを取得し、（２５）式及び（３１）式に基づき、収縮期血圧ＳＢＰｉを取得するため、血圧を簡易かつ精度よく検出できる。
【0060】
血圧測定装置１の少なくとも一部は、ハードウェアで構成してもよいし、ソフトウェアで構成してもよい。ソフトウェアで構成する場合には、血圧測定装置１の少なくとも一部の機能を実現するプログラムをフレキシブルディスクやＣＤ－ＲＯＭ等の記録媒体に収納し、コンピュータに読み込ませて実行させてもよい。記録媒体は、磁気ディスクや光ディスク等の着脱可能なものに限定されず、ハードディスク装置やメモリなどの固定型の記録媒体でもよい。
【0061】
また、血圧測定装置１の少なくとも一部の機能を実現するプログラムを、インターネット等の通信回線（無線通信も含む）を介して頒布してもよい。さらに、同プログラムを暗号化したり、変調をかけたり、圧縮した状態で、インターネット等の有線回線や無線回線を介して、あるいは記録媒体に収納して頒布してもよい。
【0062】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【0063】
１：血圧測定装置、２：測定部、４：血圧取得部、６：生体測定装置。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7A】
【図7B】
【図8】
【図9A】
【図9B】
【図10】
【図11】
【図12】
【図13A】
【図13B】
【図14】