JPWO2019150999A1 - 脈波伝播速度測定装置及びその方法 - Google Patents
脈波伝播速度測定装置及びその方法 Download PDFInfo
- Publication number
- JPWO2019150999A1 JPWO2019150999A1 JP2019569006A JP2019569006A JPWO2019150999A1 JP WO2019150999 A1 JPWO2019150999 A1 JP WO2019150999A1 JP 2019569006 A JP2019569006 A JP 2019569006A JP 2019569006 A JP2019569006 A JP 2019569006A JP WO2019150999 A1 JPWO2019150999 A1 JP WO2019150999A1
- Authority
- JP
- Japan
- Prior art keywords
- pulse wave
- velocity
- envelope
- wave
- piezoelectric vibration
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 18
- 238000012545 processing Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims description 22
- 210000004204 blood vessel Anatomy 0.000 claims description 16
- 230000000747 cardiac effect Effects 0.000 claims description 14
- 239000000284 extract Substances 0.000 claims description 5
- 238000000691 measurement method Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 abstract description 4
- 238000001914 filtration Methods 0.000 abstract 1
- 206010003210 Arteriosclerosis Diseases 0.000 description 4
- 208000011775 arteriosclerosis disease Diseases 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 210000001367 artery Anatomy 0.000 description 3
- 238000010009 beating Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 210000003127 knee Anatomy 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000002354 daily effect Effects 0.000 description 2
- 230000003203 everyday effect Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 206010005746 Blood pressure fluctuation Diseases 0.000 description 1
- 208000002251 Dissecting Aneurysm Diseases 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 206010002895 aortic dissection Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 206010008118 cerebral infarction Diseases 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000002837 heart atrium Anatomy 0.000 description 1
- 210000003709 heart valve Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000653 nervous system Anatomy 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 230000004218 vascular function Effects 0.000 description 1
- 210000000707 wrist Anatomy 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/026—Measuring blood flow
- A61B5/0285—Measuring or recording phase velocity of blood waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/291—Bioelectric electrodes therefor specially adapted for particular uses for electroencephalography [EEG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1102—Ballistocardiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
- A61B5/372—Analysis of electroencephalograms
- A61B5/374—Detecting the frequency distribution of signals, e.g. detecting delta, theta, alpha, beta or gamma waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
- A61B5/377—Electroencephalography [EEG] using evoked responses
- A61B5/38—Acoustic or auditory stimuli
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6887—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient mounted on external non-worn devices, e.g. non-medical devices
- A61B5/6891—Furniture
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- 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/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
Abstract
【課題】 圧電振動センサを人体に直接ではなく、間接的に取り付けた場合であっても、良好に脈波を検出し、その伝播速度を得る。【解決手段】 脈波や心弾波は、椅子１００の座面１０２に取り付けられた圧電振動センサ１１０に伝播し、電気信号に変換されて出力される。その後、ローパスデジタルフィルタ２０２Ｐ，バンドパスデジタルフィルタ２０２Ｂによるフィルタリング処理が施され、ローパスデジタルフィルタ２０２Ｐからは脈波ＧＰが出力され、バンドパスデジタルフィルタ２０２Ｂからは心弾波ＧＢが出力される。これら脈波ＧＰ及び心弾波ＧＢは、包絡線処理回路２１０Ｐ，２１０Ｂの絶対値回路２１２Ｐ，２１２Ｂ及びローパスフィルタ２１４Ｐ，２１４Ｂによる処理が施されて、包絡線が求められる。そして、求められた包絡線のピークの差に基づいて、脈波伝播速度ＰＷＶが演算される。【選択図】図２
Description
本発明は、心臓の拍動に伴う血管系内の血圧・体積の変化である脈波の伝播速度を測定する脈波伝播速度測定装置及びその方法に関する。
脈波に関する背景技術としては、例えば、下記特許文献１記載の「動脈硬化評価装置」がある。これは、入射波と反射波とを正確に分離して、個体差による動脈硬化度を精度よく評価できるようにしたもので、被験者の頸部において、圧電トランスデューサで動脈を伝わる脈波を変位信号として検出するとともに、超音波診断装置のプローブで動脈の血流速度を測定する。そして超音波診断装置のプローブで得られた血流速度を変位信号に変換して入射波を得た後、圧電トランスデューサによって検出された変位信号から入射波を差し引いて反射波を得、更に、入射波と反射波の振幅強度から生体の血管機能を評価するようにしたものである。
しかし、人体に直接圧電センサを取り付けることなく脈波を検出して健康状態を知ることができれば、利便性が向上し、より好都合である。例えば、椅子，自動車のシートなどに人が座ることで毎日脈波を検出し、循環器系，神経系などの健康状態を知ることができれば、昨今の高齢化社会では、極めて有益である。特に動脈硬化は、脳梗塞や大動脈解離といった重篤な疾患を引き起こす原因となるため、早期に発見することが重要である。最近の医療の分野では、この動脈硬化の指標として脈波伝播速度(PWV)が注目されており、人間ドックや健康診断においても、オプションとして測定されるようになってきている。
ところで、脈波伝播速度を測定する手法として、本件発明者らは、動脈が体表面近くに存在している人体の２ヶ所に圧電振動センサを固定し、これら２ヶ所で検出した脈波の時間差から脈波伝播速度を測定する手法を提案している。図７(A)にはその手法が示されており、指先と手首に圧電振動センサ１０，１２を取り付けている。圧電振動センサ１０，１２としては、例えば、国際公開第2016/167202号パンフレットに開示されている振動波形センサが好適な例の一つである。これら２つの圧電振動センサ１０，１２で検出した脈波の時間差を測定することで、脈波伝播速度を得ている。同図(B)には、２つの圧電振動センサ１０，１２によって検出された脈波信号が示されており、同図(C)には、同図(B)から得た脈波伝播速度が示されている。
このような脈波伝播速度を検出する圧電振動センサを、例えば、日常生活の中で毎日使用する椅子に取り付けた場合、椅子が硬質の材料でできているときは、人体で発生した脈波の振動が椅子内を音速で伝達してしまう。このため、違う場所で検出した脈波信号との間でクロストークが発生してしまうといった問題点がある。クロストークを低減する手段として、椅子にスリットを挿入するといった方法があり、相応の効果は期待できるが、椅子のデザイン上現実的ではない。
本発明は、かかる点に着目したもので、その目的は、圧電振動センサを人体に直接ではなく、間接的に取り付けた場合であっても、良好に脈波を検出し、その伝播速度を得ることである。他の目的は、少なくも１つの圧電振動センサのみであっても、良好に脈波を検出し、その伝播速度を得ることである。
本発明の脈波伝播速度測定装置は、人体が接触する面に取り付けられた圧電振動センサの出力振動波形に基づいて人の脈波伝播速度を得る脈波伝播速度測定装置であって、前記圧電振動センサの出力振動波形から、脈波を取り出す第１のフィルタ手段と、前記圧電振動センサの出力振動波形から、心弾波を取り出す第２のフィルタ手段と、前記第１のフィルタ手段で得られた脈波と、前記第２のフィルタ手段で得られた心弾波とを利用して、脈波伝播速度を演算する演算手段と、を備えたことを特徴とする。
主要な形態の一つによれば、前記第１のフィルタ手段は、前記圧電振動センサの出力振動波形から４Ｈｚ以下の周波数帯域成分を取り出し、前記第２のフィルタ手段は、前記圧電振動センサの出力振動波形から１０Ｈｚ以上３３Ｈｚ以下の周波数帯域成分を取り出すことを特徴とする。他の形態によれば、前記演算手段は、前記第１のフィルタ手段で得られた脈波の包絡線を得る第１の包絡線処理手段と、前記第２のフィルタ手段で得られた心弾波の包絡線を得る第２の包絡線処理手段と、前記第１及び第２の包絡線処理手段で得られた包絡線のピークとを利用して、前記脈波伝播速度を演算する伝播速度演算手段と、を備えたことを特徴とする。
他の形態によれば、前記伝播速度演算手段は、前記包絡線のピークから脈波と心弾波の時間差を求めるとともに、脈波が通る経路の血管長に対して、脈波伝播速度＝血管長／（脈波と心弾波の時間差）の演算を行うことを特徴とする。更に他の形態によれば、前記第１及び第２の包絡線処理手段は、絶対値回路とローパスフィルタによって構成されており、前記ローパスフィルタの遮断周波数を１．５Ｈｚ以上４Ｈｚ以下に設定したことを特徴とする。更に他の形態によれば、一つの圧電振動センサの出力振動波形から、前記脈波と心弾波を得ることを特徴とする。あるいは、異なる圧電振動センサの出力振動波形から、前記脈波と心弾波を得ることを特徴とする。
本発明の脈波伝播速度測定方法は、人体が接触する面に取り付けられた圧電振動センサの出力振動波形に基づいて人の脈波伝播速度を得る脈波伝播速度測定方法であって、前記圧電振動センサの出力振動波形から脈波を取り出す第１のステップと、前記圧電振動センサの出力振動波形から心弾波を取り出す第２のステップと、前記第１のステップで得られた脈波と、前記第２のステップで得られた心弾波とを利用して、脈波伝播速度を演算する第３のステップと、を備えたことを特徴とする。
主要な形態の一つによれば、前記第１のステップは、前記圧電振動センサの出力振動波形から４Ｈｚ以下の周波数帯域成分を取り出し、前記第２のステップは、前記圧電振動センサの出力振動波形から１０Ｈｚ以上３３Ｈｚ以下の周波数帯域成分を取り出すことを特徴とする。他の形態によれば、前記第３のステップは、前記第１のステップで得られた脈波の包絡線を得る第１の包絡線処理ステップと、前記第２のステップで得られた心弾波の包絡線を得る第２の包絡線処理ステップと、前記第１及び第２の包絡線処理ステップで得られた包絡線のピークとを利用して、前記脈波伝播速度を演算する伝播速度演算ステップと、を備えたことを特徴とする。
他の形態によれば、前記伝播速度演算ステップは、前記包絡線のピークから脈波と心弾波の時間差を求めるとともに、脈波が通る経路の血管長に対して、脈波伝播速度＝血管長／（脈波と心弾波の時間差）の演算を行うことを特徴とする。更に他の形態によれば、一つの圧電振動センサの出力振動波形から、前記脈波と心弾波を得ることを特徴とする。あるいは、異なる圧電振動センサの出力振動波形から、前記脈波と心弾波を得ることを特徴とする。本発明の前記及び他の目的，特徴，利点は、以下の詳細な説明及び添付図面から明瞭になろう。
本発明によれば、人体が接触する面に設置した圧電振動センサの出力振動波形から、フィルタ手段によって脈波及び心弾波をそれぞれ取り出し、それらを利用して脈波伝播速度を演算することとしたので、圧電振動センサを人体に直接取り付けることなく、良好に脈波伝播速度を得ることができる。
以下、本発明を実施するための最良の形態を、実施例に基づいて詳細に説明する。
最初に、図１を参照しながら、本発明の基本的な考え方について説明する。図１には、心拍ないし脈拍に関する主要な信号が示されている。まず、図１(A)は心電図で、心臓の電気的活動を体表面から波形として記録したものである。最も大きなＲ波は、心房からの刺激が刺激伝導系を通って、心室が興奮したときに起こる波形で、Ｒ波の頂点の間隔から心拍数を計算することができる。同図(B)は光電脈波で、心臓の拍動に応じて変化する血管の内圧の変化ないし容積の変化（外径の変化）である脈波を、ヘモグロビン濃度による透過光または反射光の減衰によって測定したものである。
同図(C)は心弾波で、心臓の物理的な動きに伴ってに生ずる微細な人体の振動を示し、心臓弁作動の振動が骨格等を通じて伝達されるもので、上述した脈波に比べて高速で伝達される。心弾波のピークは、同図(A)に示す心電波形や同図(B)に示す脈波に対して少し遅れて生ずる。同図(D)は心音波で、心臓の拍動に伴う弁膜の振動によって生じる音を電気信号に変換したものである。心臓の収縮時に生ずる第１音（Isound）と、拡張時に生ずる第２音（IIsound）が主なピークである。
これらのうち、心電波形及び光電脈波は周波数帯域が０〜４Ｈｚ程度であり、心弾波形及び心音波形は周波数帯域が１０〜４０Ｈｚである。従って、帯域分離を行うことで、心電波形ないし光電脈波と心弾波形ないし心音波形とを、区別することができる。本発明では、圧電振動センサを使用して得た信号波形に対して帯域分離を行って、脈波波形や心弾波形を分離するようにしている。
次に、図２〜図５も参照しながら、本発明の実施例１について説明する。図２には、本実施例における圧電振動センサの設置の様子が示されており、椅子１００の座面１０２の下側適宜位置に、圧電振動センサ１１０が設置されている。人１０４の血管が体表近くにあるほうが、脈波を測定しやすい。一方、血管は、太腿の裏側であって、膝に近い部位において、体表近くにある。そこで、圧電振動センサ１１０を座面１０２の先端側に設置することで、同図に矢印ＦＰで示すように、脈波が圧電振動センサ１１０に伝わる。一方、心弾波は、主として人１０４の骨を通じて伝播し、矢印ＦＢで示すように圧電振動センサ１１０に伝わる。
圧電振動センサ１１０には、図３(A)に示す測定回路２００が接続されており、圧電振動センサ１１０の振動波形出力側は、ローパスデジタルフィルタ２０２Ｐと、バンドパスデジタルフィルタ２０２Ｂが接続されている。両デジタルフィルタ２０２Ｐ，２０２Ｂの出力側は、包絡線処理回路２１０Ｐ，２１０Ｂにそれぞれ接続されている。包絡線処理回路２１０Ｐは、絶対値回路２１２Ｐとローパスフィルタ２１４Ｐによって構成されており、包絡線処理回路２１０Ｂは、絶対値回路２１２Ｂとローパスフィルタ２１４Ｂによって構成されている。これら包絡線処理回路２１０Ｐ，２１０Ｂの出力側は、それぞれ伝播速度演算部２２０に接続されている。
これらのうち、ローパスデジタルフィルタ２０２Ｐは、圧電振動センサ１１０の出力振動波形から、０〜４Ｈｚ（４Ｈｚ以下）の周波数帯域の成分である脈波成分を取り出すためのフィルタである。バンドパスデジタルフィルタ２０２Ｂは、圧電振動センサ１１０の出力振動波形から、１０〜３３Ｈｚの周波数帯域の成分である心弾波成分を取り出すためのフィルタである。サンプリング周波数は、例えば１０００Ｈｚとする。
包絡線処理回路２１０Ｐ，２１０Ｂは、絶対値回路２１２Ｐ，２１２Ｂによって入力信号の絶対値をとるとともに、ローパスフィルタ２１４Ｐ，２１４Ｂによるフィルタ処理を施して包絡線を得る回路である。伝播速度演算部２２０は、包絡線処理回路２１０Ｐ，２１０Ｂの出力に基づいて、脈波伝播速度ＰＷＶを演算する機能を備えており、
ａ，包絡線処理回２１０Ｐ，２１０Ｂの出力信号波形の包絡線からのピーク検出，
ｂ，脈波ＧＰの包絡線のピークと心弾波ＧＢの包絡線のピークとの差に基づく脈波伝播速度ＰＷＶの演算，
が行われるようになっている。
ａ，包絡線処理回２１０Ｐ，２１０Ｂの出力信号波形の包絡線からのピーク検出，
ｂ，脈波ＧＰの包絡線のピークと心弾波ＧＢの包絡線のピークとの差に基づく脈波伝播速度ＰＷＶの演算，
が行われるようになっている。
前記ｂの脈波伝播速度の演算は、例えば次のように行う。心弾波ＧＢの伝播速度は音速で、骨を伝播する場合１ｋｍ／ｓ程度であるのに対し、脈波ＧＰの伝播速度は１０ｍ／ｓ程度であり、両者を比較すると、心弾波ＧＢの伝播速度は、ほぼ無限大と考えることができる。してみると、心弾波ＧＢが圧電振動センサ１１０で検出された時点で脈波ＧＰが心臓を出発したと考えることができるので、心弾波ＧＢの検出から脈波ＧＰの検出に至る時間で脈波ＧＰが人体内の経路を伝播したと考えることができる。従って、脈波伝播速度ＰＷＶは、脈波ＧＰが通る経路である心臓から膝裏までの血管長に対し、脈波伝播速度(PWV)＝血管長／（脈波と心弾波の時間差）で表される。脈波と心弾波の時間差は、包絡線処理して得られた脈波ＧＰのピークと心弾波ＧＰの時間差となる。なお、血管長は、例えば、レントゲン写真等で得られている動脈血管配置図をベースに体表面をメジャーで測定するなどの方法で知ることができる。あるいは、複数人に関し、このようにして得た血管長と身長との相関関係を得ることで、身長から推定することが可能である。
次に、図４〜６も参照しながら、本実施例の動作について説明する。人１０４の心臓の拍動に基づいて生じた脈波や心弾波は、座面１０２や人１０４を通じて、座面１０２に取り付けられた圧電振動センサ１１０に伝播する（図２，矢印ＦＰ，ＦＢ参照）。圧電振動センサ１１０では、伝わった振動が電気信号に変換されて出力される。この振動波形に対しては、ローパスデジタルフィルタ２０２Ｐ，バンドパスデジタルフィルタ２０２Ｂによるフィルタリング処理が施され、ローパスデジタルフィルタ２０２Ｐからは脈波の帯域成分が出力され、バンドパスデジタルフィルタ２０２Ｂからは心弾波の帯域成分が出力される。
図４には、信号波形の一例が示されている。同図(A)は、座面１０２の比較的先端側に圧電振動センサ１１０を配置した場合であり、同図(B)は、座面１０２の比較的後端側に圧電振動センサ１１０を配置した場合である。ＧＰは脈波，ＧＢは心弾波である。これらのグラフを見ると、多少の相違はあるものの、いずれのセンサ配置であっても、脈波ＧＰ，心弾波ＧＢがそれぞれ測定されている。
圧電振動センサ１１０で得た脈波ＧＰ及び心弾波ＧＢは、包絡線処理回路２１０Ｐ，２１０Ｂにそれぞれ入力され、包絡線検出が行われる。すなわち、絶対値回路２１２Ｐ，２１２Ｂによって入力信号の絶対値が求められ、ローパスフィルタ２１４Ｐ，２１４Ｂによるフィルタ処理が施される。例えば、図５(A)に示す脈波ＧＰ，心弾波ＧＢに対して包絡線処理を行うと、同図(B)に示す包絡線ＧＰＥ，ＧＢＥがそれぞれ求められる。これらの包絡線ＧＰＥ，ＧＢＥの信号は、伝播速度演算部２２０に入力され、ここで、脈波ＧＰの包絡線ＧＰＥのピークＰＧＰＥと心弾波ＧＢの包絡線ＧＢＥのピークＰＧＢＥとの差が求められ、更に、その差に基づいて脈波伝播速度ＰＷＶの演算が行われる。
ところで、この場合において、包絡線処理回路２１０Ｐ，２１０Ｂのローパスフィルタ２１４Ｐ，２１４Ｂにおける遮断周波数をどのように設定するかにより、脈波伝播速度ＰＷＶの値も異なるようになる。上述した図５(B)は遮断周波数が２．５Ｈｚの場合であるが、これを１Ｈｚとすると図６(A)のようになり、６Ｈｚとすると同図(B)に示すようになる。このような遮断周波数と脈波伝播速度ＰＷＶとの関係を示すと、同図(C)に示すようになる。このグラフの結果からすると、遮断周波数が１．５Ｈｚ以下では、脈波ＧＰ，心弾波ＧＢのどちらの信号も信号強度が低下してしまい、これが原因で正しいピークの時間差値が得られない。またこれに対して、遮断周波数が４Ｈｚを超えてくると、心弾波ＧＢの包絡線処理が不十分になり、これまた正しい脈波伝播速度ＰＷＶを得ることができない。この理由により、遮断周波数は、１．５Ｈｚ以上４Ｈｚ以下に設定するのが好ましい。
このように、本実施例によれば、椅子１００に取り付けた圧電振動センサ１１０の出力振動波形から、脈波と心弾波を帯域抽出し、それらの包絡線のピークを利用して脈波伝播速度を演算することとしたので、
ａ，圧電振動センサを人体に直接取り付けることなく、良好に脈波伝播速度を得ることできる。
ｂ，１つの圧電振動センサのみで脈波と心弾波を同時に測定することができ、それらに基づいて良好に脈波伝播速度を得ることができる。
ａ，圧電振動センサを人体に直接取り付けることなく、良好に脈波伝播速度を得ることできる。
ｂ，１つの圧電振動センサのみで脈波と心弾波を同時に測定することができ、それらに基づいて良好に脈波伝播速度を得ることができる。
次に、図３(B)を参照しながら、実施例２について説明する。上述した実施例１では、一つの圧電振動センサ１１０を使用したが、本実施例では、同図に示すように、ローパスデジタルフィルタ２０２Ｐ，バンドパスデジタルフィルタ２０２Ｂのそれぞれに圧電振動センサ１１０Ｐ，１１０Ｂを接続している。このようにすることで、前記実施例１と比較して、圧電振動センサの数は増えるが、それぞれの圧電振動センサを、例えば、心弾波の感度が高い尾てい骨下付近と，脈波の感度が高い膝裏付近に設置することが可能になり、信号品質を更に向上させることができるという利点が生ずる。
なお、本発明は、上述した実施例に限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更を加え得ることができる。例えば、以下のものも含まれる。
(1)前記実施例では、人が日常生活の中で必ず使用する椅子の座面に圧電振動センサを取り付けたが、他に肘掛，背もたれ，ベッドや枕などの寝具等、人が接触する面であれば、各種のものに取り付けてよい。
(2)図３に示した回路構成も、例えばコンピュータプログラムによって演算を行うなど、同様の信号処理を行う各種の形態が考えられる。
(1)前記実施例では、人が日常生活の中で必ず使用する椅子の座面に圧電振動センサを取り付けたが、他に肘掛，背もたれ，ベッドや枕などの寝具等、人が接触する面であれば、各種のものに取り付けてよい。
(2)図３に示した回路構成も、例えばコンピュータプログラムによって演算を行うなど、同様の信号処理を行う各種の形態が考えられる。
本発明によれば、人体が接触する面に設置した圧電振動センサの出力振動波形から、フィルタ手段によって脈波及び心弾波をそれぞれ取り出し、それらを利用して脈波伝播速度を演算することとしたので、圧電振動センサを人体に直接取り付けることなく、良好に脈波伝播速度を得ることができ、医療の分野に好適である。
１０，１２：圧電振動センサ
１００：椅子
１０２：座面
１０４：人
１１０，１１０Ｐ，１１０Ｂ：圧電振動センサ
２００：測定回路
２０２Ｂ：バンドパスデジタルフィルタ
２０２Ｐ：ローパスデジタルフィルタ
２１０Ｐ，２１０Ｂ：包絡線処理回路
２１２Ｐ，２１２Ｂ：絶対値回路
２１４Ｐ，２１４Ｂ：ローパスフィルタ
２２０：伝播速度演算部
ＧＢ：心弾波
ＧＰ：脈波
ＧＰＥ，ＧＢＥ：包絡線
ＰＧＰＥ，ＰＧＢＥ：ピーク
ＰＷＶ：脈波伝播速度
１００：椅子
１０２：座面
１０４：人
１１０，１１０Ｐ，１１０Ｂ：圧電振動センサ
２００：測定回路
２０２Ｂ：バンドパスデジタルフィルタ
２０２Ｐ：ローパスデジタルフィルタ
２１０Ｐ，２１０Ｂ：包絡線処理回路
２１２Ｐ，２１２Ｂ：絶対値回路
２１４Ｐ，２１４Ｂ：ローパスフィルタ
２２０：伝播速度演算部
ＧＢ：心弾波
ＧＰ：脈波
ＧＰＥ，ＧＢＥ：包絡線
ＰＧＰＥ，ＰＧＢＥ：ピーク
ＰＷＶ：脈波伝播速度
Claims (13)
- 人体が接触する面に取り付けられた圧電振動センサの出力振動波形に基づいて人の脈波伝播速度を得る脈波伝播速度測定装置であって、
前記圧電振動センサの出力振動波形から、脈波を取り出す第１のフィルタ手段と、
前記圧電振動センサの出力振動波形から、心弾波を取り出す第２のフィルタ手段と、
前記第１のフィルタ手段で得られた脈波と、前記第２のフィルタ手段で得られた心弾波とを利用して、脈波伝播速度を演算する演算手段と、
を備えたことを特徴とする脈波伝播速度測定装置。 - 前記第１のフィルタ手段は、前記圧電振動センサの出力振動波形から４Ｈｚ以下の周波数帯域成分を取り出し、前記第２のフィルタ手段は、前記圧電振動センサの出力振動波形から１０Ｈｚ以上３３Ｈｚ以下の周波数帯域成分を取り出すことを特徴とする請求項１記載の脈波伝播速度測定装置。
- 前記演算手段は、
前記第１のフィルタ手段で得られた脈波の包絡線を得る第１の包絡線処理手段と、前記第２のフィルタ手段で得られた心弾波の包絡線を得る第２の包絡線処理手段と、
前記第１及び第２の包絡線処理手段で得られた包絡線のピークとを利用して、前記脈波伝播速度を演算する伝播速度演算手段と、
を備えたことを特徴とする請求項１又は２記載の脈波伝播速度測定装置。 - 前記伝播速度演算手段は、
前記包絡線のピークから脈波と心弾波の時間差を求めるとともに、脈波が通る経路の血管長に対して、脈波伝播速度＝血管長／（脈波と心弾波の時間差）の演算を行うことを特徴とする請求項３記載の脈波伝播速度測定装置。 - 前記第１及び第２の包絡線処理手段は、絶対値回路とローパスフィルタによって構成されており、
前記ローパスフィルタの遮断周波数を１．５Ｈｚ以上４Ｈｚ以下に設定したことを特徴とする請求項３又は４記載の脈波伝播速度測定装置。 - 一つの圧電振動センサの出力振動波形から、前記脈波と心弾波を得ることを特徴とする請求項１〜５のいずれか一項に記載の脈波伝播速度測定装置。
- 異なる圧電振動センサの出力振動波形から、前記脈波と心弾波を得ることを特徴とする請求項１〜５のいずれか一項に記載の脈波伝播速度測定装置。
- 人体が接触する面に取り付けられた圧電振動センサの出力振動波形に基づいて人の脈波伝播速度を得る脈波伝播速度測定方法であって、
前記圧電振動センサの出力振動波形から脈波を取り出す第１のステップと、
前記圧電振動センサの出力振動波形から心弾波を取り出す第２のステップと、
前記第１のステップで得られた脈波と、前記第２のステップで得られた心弾波とを利用して、脈波伝播速度を演算する第３のステップと、
を備えたことを特徴とする脈波伝播速度測定方法。 - 前記第１のステップは、前記圧電振動センサの出力振動波形から４Ｈｚ以下の周波数帯域成分を取り出し、前記第２のステップは、前記圧電振動センサの出力振動波形から１０Ｈｚ以上３３Ｈｚ以下の周波数帯域成分を取り出すことを特徴とする請求項８記載の脈波伝播速度測定方法。
- 前記第３のステップは、
前記第１のステップで得られた脈波の包絡線を得る第１の包絡線処理ステップと、前記第２のステップで得られた心弾波の包絡線を得る第２の包絡線処理ステップと、
前記第１及び第２の包絡線処理ステップで得られた包絡線のピークとを利用して、前記脈波伝播速度を演算する伝播速度演算ステップと、
を備えたことを特徴とする請求項８又は９記載の脈波伝播速度測定方法。 - 前記伝播速度演算ステップは、
前記包絡線のピークから脈波と心弾波の時間差を求めるとともに、脈波が通る経路の血管長に対して、脈波伝播速度＝血管長／（脈波と心弾波の時間差）の演算を行うことを特徴とする請求項１０記載の脈波伝播速度測定方法。 - 一つの圧電振動センサの出力振動波形から、前記脈波と心弾波を得ることを特徴とする請求項８〜１１のいずれか一項に記載の脈波伝播速度測定方法。
- 異なる圧電振動センサの出力振動波形から、前記脈波と心弾波を得ることを特徴とする請求項８〜１１のいずれか一項に記載の脈波伝播速度測定方法。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018015080 | 2018-01-31 | ||
| JP2018015080 | 2018-01-31 | ||
| PCT/JP2019/001568 WO2019150999A1 (ja) | 2018-01-31 | 2019-01-19 | 脈波伝播速度測定装置及びその方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2019150999A1 true JPWO2019150999A1 (ja) | 2020-11-19 |
| JP6874166B2 JP6874166B2 (ja) | 2021-05-19 |
Family
ID=67478142
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019569006A Active JP6874166B2 (ja) | 2018-01-31 | 2019-01-19 | 脈波伝播速度測定装置及びその方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20200359929A1 (ja) |
| JP (1) | JP6874166B2 (ja) |
| WO (1) | WO2019150999A1 (ja) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004305268A (ja) * | 2003-04-02 | 2004-11-04 | Colin Medical Technology Corp | 心音検出装置 |
| JP2011527589A (ja) * | 2008-07-11 | 2011-11-04 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 心弾動図解析方法及び装置 |
| JP2014507213A (ja) * | 2011-01-27 | 2014-03-27 | ザ ボード オブ トラスティーズ オブ ザ レランド スタンフォード ジュニア ユニバーシティー | 循環系を観察するためのシステム及び方法 |
| JP2016026516A (ja) * | 2013-12-07 | 2016-02-18 | 株式会社デルタツーリング | 音・振動情報収集機構、音・振動情報センシングシステム及びコンピュータプログラム |
| WO2017055670A1 (es) * | 2015-10-02 | 2017-04-06 | Universitat Politècnica De Catalunya | Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma |
| WO2017141976A1 (ja) * | 2016-02-15 | 2017-08-24 | ヘルスセンシング株式会社 | 睡眠状態測定装置及び方法、位相コヒーレンス算出装置、生体振動信号測定装置、ストレス状態測定装置及び睡眠状態測定装置及び心拍波形抽出方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4480785B1 (ja) * | 2009-04-23 | 2010-06-16 | メディカルトラスト株式会社 | 生体情報検出装置 |
| US11589758B2 (en) * | 2016-01-25 | 2023-02-28 | Fitbit, Inc. | Calibration of pulse-transit-time to blood pressure model using multiple physiological sensors and various methods for blood pressure variation |
-
2019
- 2019-01-19 US US16/963,092 patent/US20200359929A1/en active Pending
- 2019-01-19 JP JP2019569006A patent/JP6874166B2/ja active Active
- 2019-01-19 WO PCT/JP2019/001568 patent/WO2019150999A1/ja active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004305268A (ja) * | 2003-04-02 | 2004-11-04 | Colin Medical Technology Corp | 心音検出装置 |
| JP2011527589A (ja) * | 2008-07-11 | 2011-11-04 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 心弾動図解析方法及び装置 |
| JP2014507213A (ja) * | 2011-01-27 | 2014-03-27 | ザ ボード オブ トラスティーズ オブ ザ レランド スタンフォード ジュニア ユニバーシティー | 循環系を観察するためのシステム及び方法 |
| JP2016026516A (ja) * | 2013-12-07 | 2016-02-18 | 株式会社デルタツーリング | 音・振動情報収集機構、音・振動情報センシングシステム及びコンピュータプログラム |
| WO2017055670A1 (es) * | 2015-10-02 | 2017-04-06 | Universitat Politècnica De Catalunya | Método y aparato para estimar el tiempo de tránsito del pulso aórtico a partir de intervalos temporales medidos entre puntos fiduciales del balistocardiograma |
| WO2017141976A1 (ja) * | 2016-02-15 | 2017-08-24 | ヘルスセンシング株式会社 | 睡眠状態測定装置及び方法、位相コヒーレンス算出装置、生体振動信号測定装置、ストレス状態測定装置及び睡眠状態測定装置及び心拍波形抽出方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6874166B2 (ja) | 2021-05-19 |
| WO2019150999A1 (ja) | 2019-08-08 |
| US20200359929A1 (en) | 2020-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5016718B2 (ja) | 動脈硬化評価装置 | |
| Yang et al. | Pulse transit time measurement using seismocardiogram, photoplethysmogram, and acoustic recordings: Evaluation and comparison | |
| JP6131404B2 (ja) | 心臓の機能不全と異常を表す情報を測定するための方法及び器械 | |
| US20220378299A1 (en) | Noninvasive method for measuring sound frequencies created by vortices in a carotid artery, visualization of stenosis, and ablation means | |
| US6328698B1 (en) | Diagnostic system and method for coronary artery disease and others | |
| WO2009076566A1 (en) | Systems and methods for cardiac contractility analysis | |
| Saito et al. | Noninvasive assessment of arterial stiffness by pulse wave analysis | |
| CN107106118B (zh) | 检测重搏切迹的方法 | |
| US20190083038A1 (en) | System and method for detecting decoupling of an auscultatory sound sensor from a test-subject | |
| WO2008021118A2 (en) | Systems and methods for calibration of heart sounds | |
| CN110115574A (zh) | 心率监测的方法和装置 | |
| JP5328613B2 (ja) | 脈波伝播速度測定装置および脈波伝播速度測定プログラム | |
| JP7177443B2 (ja) | 生体振動信号検出装置 | |
| JP6129166B2 (ja) | 動脈の閉塞／再開を検出するための方法及び装置並びに収縮期血圧を測定するためのシステム | |
| CN105748038B (zh) | 一种冠心病风险指数的无损检测装置 | |
| JP6874166B2 (ja) | 脈波伝播速度測定装置及びその方法 | |
| Arathy et al. | An accelerometer probe for local pulse wave velocity measurement | |
| Jing et al. | Noninvasive acoustical analysis system of coronary heart disease | |
| US20210145398A1 (en) | Electronic stethoscope with enhanced features | |
| JP2019130035A (ja) | 脈波伝播速度測定装置 | |
| JP5016717B2 (ja) | 動脈硬化評価装置 | |
| Elgendi et al. | Measurement of aa intervals at rest in the second derivative plethysmogram | |
| JP2006271949A (ja) | 脈波から呼吸信号を検出する方法 | |
| JP5344294B2 (ja) | 血管壁硬軟度評価装置 | |
| KR20020013820A (ko) | 심음과 맥파를 이용한 맥파전달속도 측정시스템 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200420 |
|
| 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: 20210329 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210421 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6874166 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |