Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP6736110B2 - 動脈血圧に関する脈波伝播時間の補正方法 - Google Patents
[go: Go Back, main page]

JP6736110B2 - 動脈血圧に関する脈波伝播時間の補正方法 - Google Patents

動脈血圧に関する脈波伝播時間の補正方法 Download PDF

Info

Publication number
JP6736110B2
JP6736110B2 JP2019541843A JP2019541843A JP6736110B2 JP 6736110 B2 JP6736110 B2 JP 6736110B2 JP 2019541843 A JP2019541843 A JP 2019541843A JP 2019541843 A JP2019541843 A JP 2019541843A JP 6736110 B2 JP6736110 B2 JP 6736110B2
Authority
JP
Japan
Prior art keywords
pulse wave
toe
blood pressure
correction
ear
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
Application number
JP2019541843A
Other languages
English (en)
Other versions
JP2019530560A (ja
JP2019530560A5 (ja
Inventor
妍 陳
妍 陳
瑜 陳
瑜 陳
Original Assignee
浙江脉聯医療設備有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Priority claimed from CN201611045054.5A external-priority patent/CN106580303B/zh
Priority claimed from CN201611046184.0A external-priority patent/CN106377238B/zh
Application filed by 浙江脉聯医療設備有限公司 filed Critical 浙江脉聯医療設備有限公司
Publication of JP2019530560A publication Critical patent/JP2019530560A/ja
Publication of JP2019530560A5 publication Critical patent/JP2019530560A5/ja
Application granted granted Critical
Publication of JP6736110B2 publication Critical patent/JP6736110B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026Measuring blood flow
    • A61B5/0285Measuring or recording phase velocity of blood waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/021Measuring pressure in heart or blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A61B5/02125Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave propagation time
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/01Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A61B5/02116Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave amplitude

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Hematology (AREA)
  • Vascular Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Description

本発明は、動脈血圧測定の技術分野に関し、具体的には、拡張期血圧及び収縮期血圧に関する脈波伝播時間の補正方法に関する。
動脈血圧は、循環器系の状態を反映し、臓器灌流を評定する主な指標の一つであり、周術期モニタリングの重要な生命徴候パラメータである。現在、周術期に一般的に用いられる血圧モニタリング方法は、侵襲的測定と非侵襲的測定に分けられる。侵襲的測定とは、専用パイプを体の循環器系に配置し、変換器を介して機械的位置エネルギーを電子信号に変換して、モニタリング装置に血圧変化をリアルタイムに表示する技術である。侵襲的測定方法は、1拍ごとの血圧を連続的で正確に測定できるが、引き起こされる危険や傷害は無視できない。非侵襲的測定で一般的に用いられる方法はカフ・オシロメトリック法であり、このような技術は操作しやすいだけでなく、精度も臨床において認められており、身体検査や周術期モニタリングにおいて幅広く用いられている。しかしながら、カフ・オシロメトリック法では3〜5分間おきに血圧を断続的に測定することができるが、動脈血圧の変化をリアルタイムに追跡することはできない。
このため、医学界では、1拍ごとの血圧を連続的で非侵襲的に測定する技術を提案しており、脈波伝播時間/速度(PTT/PWV)を利用して1拍ごとの血圧を連続的で非侵襲的に測定する方法は研究のホットスポットとなりつつある。当該測定方法では、1つまたは複数の光電子センサを一組の心電図電極と同期させることによって、体積脈波(PhotoPlethysmoGraphy PPG)及び心電図信号(ECG)を取得し、PPGとECGとの間の時間差または2つのPPG間の時間差を用いてPTT/PWVを算出し、PTT/PWVと血圧との間の関数関係を調べて数学モデルを確立し、測定可能なPTT/PWVを用いて血圧を推定する。多くの学術論文では、PTT/PWVを用いて1拍ごとの血圧を連続的で非侵襲的に測定する原理を開示しており、たとえば、Yan Chen、Changyun Wen、Guocai Tao、Min Bi、及びGuoqi Li著の『A Novel Modeling Methodology of the Relationship Between Blood Pressure and Pulse Wave Velocity』、Yan Chen、Changyun Wen、Guocai Tao、及びMin Bi著の『Continuous and Noninvasive Measurement of Systolic and Diastolic Blood Pressure by One Mathematical Model with the Same Model Parameters and Two Separate Pulse Wave Velocities』、Younhee Choi、Qiao Zhang、Seokbum Ko著の『Noninvasive cuffless blood pressure estimation using pulse transit time and Hilbert−Huang transform』、Zheng Y、Poon CC、Yan BP、Lau JY著の『Pulse Arrival Time Based Cuff−Less and 24−H Wearable Blood Pressure Monitoring and its Diagnostic Value in Hypertension』、Mukkamala R、Hahn JO、Inan OT、Mestha LK、Kim CS、To¨reyin H、Kyal S著の『Toward Ubiquitous Blood Pressure Monitoring via Pulse Transit Time: Theory and Practice』が挙げられる。特許文献1〜8など多数の特許には、PTT/PWVを用いて1拍ごとの血圧を連続的で非侵襲的に測定する具体的な実施方法または装置が開示されている。
PTT/PWVを用いた従来の血圧測定方法及び技術によれば、初期較正のために従来のカフ・オシロメトリック法を用いて1つまたは一組の血圧値を測定する必要があり、較正の理由は、PTT/PWVと血圧との相関関係が対象に依存し、すなわち、個人ごとにPTT/PWVと血圧との間に明確な関係があることにあり、較正の目的は、対象に適した数学モデルパラメータを決定することにある。
ところが、従来方法には、循環器系が外部からの干渉を受けない条件下しか適用できないという限りがある。干渉がない場合のみに、個体のPTTと血圧との関係は特定の関数と数学モデルによって説明できるほど強い規律性を有する。しかし、周術期では、輸液、薬物、外科手術、温度などの複雑な要素による影響下で患者の循環器系において、PTTが一連の変異が生じるため、変異が発生したPTT及び固有の数学モデルを用いて血圧を推定すると、大きな誤差が生じる。変異が発生したPTTと血圧との関係には所定の規律性を持たないので、PTTの変異に適応するように数学モデルパラメータを頻繁に較正しても、根本的な問題を解決できず、臨床測定による正確性やリアルタイム性への要件を満たすことができない。
中国特許出願公開第101229058号明細書 中国特許出願公開第102811659号明細書 中国特許第1127939号明細書 米国特許第5865755号明細書 米国特許第5857975号明細書 米国特許第5649543号明細書 米国特許第9364158号明細書 欧州特許第0413267号明細書
従来技術に存在する欠陥に鑑み、本発明は、臨床条件で輸血・輸液、血管作用薬、外科的介入などの原因により引き起こされる脈波伝播時間の変異を自己適応的に補正でき、精度の高い、拡張期血圧に関する脈波伝播時間PTTの補正方法を提供することを目的とする。
本発明はさらに、臨床条件で輸血・輸液、血管作用薬、外科的介入などの原因により引き起こされる脈波伝播時間の変異を自己適応的に補正でき、精度の高い、収縮期血圧に関する脈波伝播時間PTTの補正方法を提供することも目的とする。
第1態様では、拡張期血圧に関する脈波伝播時間の補正方法は、
各心周期内の耳脈波をリアルタイムに検出して分析し、耳脈波における大動脈弁閉鎖点の高さhsd、耳脈波の収縮期時間t(単位はミリ秒)、耳脈波の拡張期時間t(単位はミリ秒)、耳脈波の最大高さhmaxを含む耳脈波のデータを取得するステップS1と、
各心周期内の足指脈波をリアルタイムに検出して分析し、足指脈波の収縮期時間ts−toe(単位はミリ秒)、足指脈波の拡張期時間td−toe(単位はミリ秒)、足指脈波の最大高さhmax−toe、足指脈波の開始点から、ピークでの立ち上がりエッジの転換点と立ち下がりエッジの転換点との中点であるピーク中点までの時間tch−toe(単位はミリ秒)、足指脈波の開始点からピークの最高点までの時間tmax−toe(単位はミリ秒)を含む足指脈波のデータを取得するステップS2と、
拡張期血圧に関する脈波伝播時間Tを計算するステップS3であって、該Tが耳脈波の開始点から足指脈波の開始点までの時間差であり、hが耳脈波または足指脈波の縦軸方向での振幅であるステップS3と、
同一心周期にステップS1、S2で取得したデータを用いて、該心周期での補正変数を算出するステップS4と、
ステップS4で取得した心周期での補正変数に基づき、該心周期での補正行列を算出するステップS5と、
複数の心周期での補正行列を連続的に取得し、ステップS3で取得したTを補正するステップS6と、を含む。
好ましくは、前記ステップS5における補正行列は、下式である。
好ましくは、前記ステップS6は、具体的に、8つの心周期での補正行列を連続的に取得し、
とする方法で補正を行うことである。
好ましくは、第1補正変数bは、下式に従って算出される。
好ましくは、第2補正変数bは、下式に従って算出される。
好ましくは、第3補正変数bは、下式に従って算出される。
好ましくは、第4補正変数bは、下式に従って算出される。
好ましくは、第5補正変数bは、下式に従って算出される。
好ましくは、第6補正変数bは、下式に従って算出される。
好ましくは、第7補正変数bは、下式に従って算出される。
この拡張期血圧に関する脈波伝播時間の補正方法は、同一心周期内の耳脈波と足指脈波とをリアルタイムに検出し、拡張期血圧に関する脈波伝播時間を計算し、そして脈波の形態学的特徴に基づき補正変数を抽出し、補正行列を取得することにより、上記脈波伝播時間の変異を自己適応的に補正し、補正済みの伝播時間は従来の数学モデルに用いて、臨床条件で各心周期の拡張期血圧を連続的で正確に測定できる。
第2態様では、収縮期血圧に関する脈波伝播時間の補正方法は、
各心周期内の耳脈波をリアルタイムに検出して分析し、耳脈波における大動脈弁閉鎖点の高さhsd、耳脈波の収縮期時間t(単位はミリ秒)、耳脈波の拡張期時間t(単位はミリ秒)、耳脈波の最大高さhmaxを含む耳脈波のデータを取得するステップS1と、
各心周期内の足指脈波をリアルタイムに検出して分析し、足指脈波の収縮期時間ts−toe(単位はミリ秒)、足指脈波の拡張期時間td−toe(単位はミリ秒)、足指脈波の最大高さhmax−toe、足指脈波の開始点から、ピークでの立ち上がりエッジの転換点と立ち下がりエッジの転換点との中点であるピーク中点までの時間tch−toe(単位はミリ秒)、足指脈波の開始点からピークの最高点までの時間tmax−toe(単位はミリ秒)を含む足指脈波のデータを取得するステップS2と、
収縮期血圧に関する脈波伝播時間Tを計算するステップS3であって、該Tが耳脈波における大動脈弁閉鎖点から足指脈波における大動脈弁閉鎖点までの時間差であり、hが耳脈波または足指脈波の縦軸方向での振幅であるステップS3と、
同一心周期にステップS1、S2で取得したデータを用いて、該心周期での補正変数を算出するステップS4と、
ステップS4で取得した心周期での補正変数に基づき、該心周期での補正行列を算出するステップS5と、
複数の心周期での補正行列を連続的に取得し、ステップS3で取得したTを補正するステップS6と、を含む。
好ましくは、前記ステップS5における補正行列は、下式である。
好ましくは、前記ステップS6は、具体的に、8つの心周期での補正行列を連続的に取得し、
とする方法で補正を行うことである。
好ましくは、第1補正変数aは、下式に従って算出される。
好ましくは、第2補正変数aは、下式に従って算出される。
好ましくは、第3補正変数aは、下式に従って算出される。
好ましくは、第4補正変数aは、下式に従って算出される。
好ましくは、第5補正変数aは、下式に従って算出される。
好ましくは、第6補正変数aは、下式に従って算出される。
好ましくは、第7補正変数aは、下式に従って算出される。
該収縮期血圧に関する脈波伝播時間の補正方法は、同一心周期内の耳脈波と足指脈波をリアルタイムに検出し、収縮期血圧に関する脈波伝播時間を計算し、そして脈波の形態学的特徴に基づき補正変数を抽出し、補正行列を取得することにより、上記脈波伝播時間の変異を自己適応的に補正し、補正済みの伝播時間は従来の数学モデルに用いて、臨床条件で各心周期の拡張期血圧を連続的で正確に測定できる。
以下、本発明の技術案の実施形態を詳細に説明する。以下の実施形態は本発明の技術案をより明確に説明するためのものに過ぎず、従って、単なる例示であり、本発明の保護範囲を限定することを意図するものではない。
周術期PTTの変化は2クラスに分けられる。クラス1の変化は血圧変化によるPTT変化、クラス2の変化はPTTと血圧とが同期していない変化(両者の変化方向または変化量が通常の関数規則に合致しない)である。たとえば、血液量がわずかに足りない場合は、PTTが増加するが、体自身による末梢抵抗への調整のため、血圧はあまり変化しない。また、胸部や腹部の手術中におけるフックの使用はPTTに深刻な影響を与える恐れがあるが、血圧への影響は小さい。さらに、ノルエピネフリンは小動脈を強く収縮させて、その結果、血圧が著しく高くなるが、全身平均PTTへの影響は小さい。
PTTにクラス1の変化が発生した場合、PTTと血圧との関係は依然として特定の関数によって表現することができ、数学モデルを用いて血圧の変化を推定できる。PTTにクラス2の変化が発生した場合、通常の循環器系に基づく数学モデルを用いて血圧を推定すると、大きな誤差が生じる。このような誤差は、PTTを用いた血圧測定による原理的誤差であり、数学モデルパラメータの初期キャリブレーション及び定期的な較正によって解決することができない。個体間のPTTの差異と同一個体のPTTの変異は、性質が異なる二種類の問題であるため、異なる方法で解決しなければならない。このため、本発明では、脈波の形態変化に応じて複数種の変数を抽出して、PTTの各種のクラス2の変化を間接的に識別して自己適応的に補正し、上記原理的誤差を克服し、従来の数学モデルを組み合わせることによって、自己適応較正機能を有する血圧の連続的な非侵襲的測定方法を確立し、これにより、カフ・オシロメトリック法のような常法を用いて繰り返して較正することが不要になる。
脈波を検出する人体の位置として、耳と足指が好ましく、この2つの部位での脈波は、大動脈と末梢動脈の生理学的、病理学的情報を取得することができ、伝播経路において代表的なものである。パルス信号を検出するセンサは、好ましくは赤外フォトプレチスモグラフ(PPG)である。
耳と足指の脈波自体の形態変化及び2種の脈波間の形態の相対変化は、PTTのクラス2の変化及び人体の異なる部位間の血圧の差を識別するために豊富な情報を提供する。本発明は、数年間にわたって、大量の手術症例の侵襲的動脈血圧、耳と足指の脈波波形及びPTTを収集して分析し、2つの脈波自体及び相対形態変化に基づき複数種の変数を抽出し、各種の変数とPTTの異なるクラス2の変化との間の関係を研究し、かつ各種変数の適用範囲を定義する。
臨床応用では、PPTによる血圧の連続測定中に、脈波波形をリアルタイムに分析し変数を抽出し、変数が適用範囲にあるか否かに基づきPTTにクラス2の変化が発生したか否かを判断する。そして、適用変数の性質に従ってPTTのクラス2の変化の性質と程度を決定し、ある変数が適用範囲を超えると、PTTには対応したクラス2の変化が発生しないことを示し、この場合、該変数は適用できない。適用可能な複数種の変数を融合して、補正量を算出してPTTを補正すると、補正されたPTT/PWVを従来の数学モデルに用いて、血圧を正確に計算できる。
本発明は、脈波形態の最も重要で基本的な変化規律を限りのある変数で表現し、そして、これらの規律とPTTとの間の関係を研究している。以下に説明する脈波は、平面座標において、縦軸が振幅h、横軸が時間t、脈波開始点が座標原点である。
(実施形態1)
拡張期血圧に関する脈波伝播時間の補正方法は、以下のステップを含む。
S1)各心周期内の耳脈波をリアルタイムに検出して分析し、耳脈波における大動脈弁閉鎖点の高さhsd(すなわち、収縮期と拡張期の耳脈波に現れる境界部での高さ)、耳脈波の収縮期時間t(単位はミリ秒)、耳脈波の拡張期時間t(単位はミリ秒)、耳脈波の最大高さhmaxを含む耳脈波のデータを取得する。
S2)各心周期内の足指脈波をリアルタイムに検出して分析し、足指脈波の収縮期時間ts−toe(単位はミリ秒)、足指脈波の拡張期時間td−toe(単位はミリ秒)、足指脈波の最大高さhmax−toe、足指脈波の開始点からピーク中点までの時間tch−toe(単位はミリ秒)、足指脈波の開始点からピークまでの時間tmax−toe(単位はミリ秒)を含む足指脈波のデータを取得する。ピーク中点は、ピークでの立ち上がりエッジの転換点と立ち下がりエッジの転換点との中点である。ピーク中点の定義について、YAN CHEN、CHANGYUN WEN、GUOCAI TAO、及びMIN BI著の『Continuous and Noninvasive Measurement of Systolic and Diastolic Blood Pressure by One Mathematical Model with the Same Model Parameters and Two Separate Pulse Wave Velocities』を参照することができる。
S3)拡張期血圧に関する脈波伝播時間Tを計算する。Tの定義については、YAN CHEN、CHANGYUN WEN、GUOCAI TAO、及びMIN BI著の『Continuous and Noninvasive Measurement of Systolic and Diastolic Blood Pressure by One Mathematical Model with the Same Model Parameters and Two Separate Pulse Wave Velocities』を参照することができる。hは耳脈波または足指脈波の縦軸方向での振幅である。
S4)同一心周期にステップS1、S2で取得したデータを用いて、該心周期での補正変数を算出する。
S5)ステップS4で取得した心周期での補正変数に基づき、該心周期での補正行列を算出する。
S6)複数の心周期での補正行列を連続的に取得し、ステップS3で取得したTを補正する。
この方法は、同一心周期内の耳脈波と足指脈波をリアルタイムに検出し、拡張期血圧に関する脈波伝播時間を計算し、そして脈波の形態学的特徴に基づき補正変数を抽出し、補正行列を取得し、上記脈波伝播時間の変異を補正することができ、補正済みの伝播時間を従来の数学モデルに用いて、臨床条件で各心周期の拡張期血圧を連続的に測定できる。
(第1補正変数b
ステップS4で取得した補正変数は、低血圧状態での拡張期血圧に関する伝播時間Tのクラス2の変化の補正に用いられる第1補正変数bを含み、bの適用範囲はb>0であり、bが大きいほど血圧が低いことを示す。
sd−m−0は、hsdと耳脈波収縮期の平均高さとの比を表す。一部の症例は、低血圧状態で、脈波ピークが前傾した三角形として現れ、hsdが大きく減少し、ksd−m−0が小さくなるが、これは、大動脈収縮期の最後部分の波形が大きく低下することを示し、脈波伝播を促進する持続動力が不足であり、伝播時間が長くなると推定される。このような状態では、拡張期情報が不安定になり、使用に不適である。
1−b=74〜82であり、好ましくは78である。d1−2−b=98〜106であり、好ましくは102である。
脈波伝播を促進する持続動力が不足である場合、伝播時間Tが長くなり、bで補正する必要がある。すなわち、d1−b≦ksd−m−0≦d1−2−bである場合、b=(d1−2−b−ksd−m−0)×0.4とする。
脈波伝播を促進する持続動力が深刻に不足になる場合、伝播時間Tがはるかに長くなり、bの上限値で補正する必要がある。すなわち、ksd−m−0<d1−bである場合、b=24×0.4とする。
脈波伝播を促進する持続動力が十分である場合、Tを補正する必要がなく、bが適用できない。すなわち、ksd−m−0>d1−2−bである場合、b=0とする。
(第2補正変数b
ステップS4で取得した補正変数は、さらに、高血圧状態で拡張期血圧に関する伝播時間Tのクラス2の変化の補正に用いられる第2補正変数bを含み、bの適用範囲はb>0であり、bが大きいほど拡張期血圧が高いことを示す。
sd−m−tsは、hsdと耳脈波拡張期のt−2t部分での平均高さとの比を表し、脈波拡張期の変異を判断することに用いられる。たとえば、胸部や腹部手術中におけるフックの使用は大動脈受力の変化を引き起こし、耳脈波拡張期の波形を低くして、ksd−m−tsを増大させる。
sd−m−2は、hsdと耳脈波の0−2t部分での平均高さとの比率を表し、収縮期と一部の拡張期の波形情報を含み、主に、気管内挿管による心拍数と血圧の上昇のような高血圧状態に用いられる。高血圧状態では、耳脈波は正三角形または後方へ傾斜した三角形を呈し、hsdがはるかに上昇し、ksd−m−2が大きくなる。正常の血圧状態での波形に比べて、高血圧状態での波形の立ち上がりエッジの傾きが小さくなり、脈波伝播を促進する動力が不足になり、伝播時間Tが長くなる。
|ksd−m−0−ksd−m−ts|≧40かつ(ksd−m−0+ksd−m−ts)/2≧ksd−m−2
ならば、ksd−m=2×ksd−m−2−(ksd−m−0+ksd−m−ts)/2、
さもなくば、ksd−m=ksd−m−2
これは、耳脈波拡張期の波形が変化する場合、たとえば、胸部や腹部手術中におけるフックの使用により大動脈の受力が変化した結果、脈波拡張期の形態が大きく変化する場合、ksd−mを補正し、そうでない場合はksd−m=ksd−m−2とするということである。
2−b=1.33〜1.43であり、好ましくは1.38である。
sd−m>(d2−b+(age−14)/15/100)(ageは年齢)であり、拡張期血圧に対応する持続動力が不足になり、それに対応して伝播時間Tが長くなり、bで補正する必要がある場合、b=(ksd−m−(d2−b+(age−14)/15/100))×0.5とする。bの変化は脈波の立ち上がりエッジの傾きの変化に反比例し、ここで、0.5は比例係数である。
sd−m≦(d2−b+(age−14)/15/100)であり、拡張期血圧に対応する持続動力が十分であり、bが適用できない場合、b=0とする。
(第3補正変数b
ステップS3で取得した補正変数は、さらに、血液量が変化しまたはセンサの配置部位の体温が変化した状態でTを補正するための第3補正変数bを含む。
d−m−tdは、耳脈波拡張期の平均高さと最大高さhmaxとの比率である。手術前に患者が断食しかつ飲水量を減少させると、血液量が低下して、kd−m−tdが減少し、脈波伝播時間が長くなり、手術中に輸血・輸液により血液量が増加すると、kd−m−tdが大きくなり、伝播時間が短くなる。(訳注:「kd−m−td」の下付き部分は、「d−m−t」とすべきだが、便宜上、「t」の「d」の下付きを解除して示す。同様に、明細書中で、下付き文字を含むパラメータが下付き文字として用いられる場合、当該パラメータ中の下付きを適宜解除して示す。)
sd−m−ts≦d3−2である場合、耳脈波の拡張早期での波形が正常な範囲を超えるまで上昇することを表し、kd−m−tdを修正する必要があり、修正結果をkd−m−td−1として表記する。
d−m−td−1=kd−m−td−(d3−2−ksd−m−ts)×75/100であり、kd−m−td≦dである場合、耳脈波が干渉を受けたと判断でき、この場合、kd−m−td−1=dとする。d=0.02〜0.14であり、好ましくは0.08である。d3−2=1.21〜1.31であり、好ましくは1.26である。
d−m−td−toeは、足指脈波の拡張期での平均高さと最大高さhmax−toeとの比率であり、ts−toeは足指脈波において識別された心臓収縮期時間を表し、td−toeは足指脈波において識別された拡張期時間を表す。kd−m−td−toe≦dである場合、kd−m−td−toe=dとする。kd−m−td−toeとkd−m−tdの作用及び性質は同じである。
d−m−a=(kd−m−td−1+kd−m−td−toe)/2として、耳と足指の脈波の性質が同じである2つの変数を合併し、その平均値を脈波伝播時間を補正する変数としたとき、脈波拡張期の波形に変化が発生した場合に、kd−m−aを補正する。
|ksd−m−0−ksd−m−ts|≧40かつ(ksd−m−0+ksd−m−ts)/2≧ksd−m−2かつksd−m−ts≧d3−2である場合、
d−ma−a=(kd−m−td−1+kd−m−td−toe+(ksd−m−0+ksd−m−ts)/2−ksd−m−2)/2とする。
血液量が正常であるとともに、センサの配置部位の体温が正常である状態には、bが適用できない。すなわち、c<kd−m−a<cである場合、b=0とする。c=(d+(age−14)/8)/100、d=23〜35、好ましくは29、c=(d+(age−14)/8)/100、d=27〜39、好ましくは33である。
血圧が非常に低くまたは非常に高い状態では、拡張期情報が不安定になり、bが適用できない。すなわち、ksd−m−0<dまたはksd−m−2>dである場合、b=0とする。d=0.97〜1.03であり、好ましくは1.00である。d=1.52〜1.58であり、好ましくは1.55である。
正常血圧状態では、血液量が低下しまたはセンサの配置部位の体温が低下する場合、bは正の値の67%とされる。すなわち、ksd−m−0≧d+0.10かつksd−m−2≦dかつkd−m−a≦cである場合、b=(c−kd−m−a)×67/100とする。d=1.42〜1.48であり、好ましくは1.45である。
低血圧または高血圧状態では、血液量が低下しまたはセンサの配置部位の体温が低下する場合、bは正常血圧状態の数値の50%とされる。すなわち、下式を満たす場合、b=(c−kd−m−a)×50/100とする。
正常血圧状態では、血液量が増加しまたはセンサの配置部位の体温が上昇する場合、bは負の値の62%とされる。すなわち、ksd−m−0≧d+0.10かつksd−m−2≦dかつkd−m−a≧cである場合、b=(c−kd−m−a)×62/100とする。
低血圧または高血圧状態では、血液量が増加しまたはセンサの配置部位の体温が上昇する場合、bは正常血圧状態の負の値の45%とされる。すなわち、下式を満たす場合、b=(c−kd−m−a)×45/100とする。
(第4補正変数b
ステップS4で取得した補正変数は、さらに、末梢血管拡張により下肢血圧(橈骨動脈血圧に対する)の低下をもたらす場合にTを補正する第4補正変数bを含み、bの適用範囲はb>0であり、bが大きいほど、橈骨動脈血圧に対する下肢血圧の低下が大きいことを示す。
末梢血管の収縮と拡張は、足指脈波のピーク位置の時間軸における前後移動を引き起こす。tmax−toe≧tch−toeである場合、ks−t−toe=[tmax−toe+tch−toe+400]/[(ts−toe+200)×2]とし、そうでない場合、ks−t−toe=[tmax−toe+200]/[ts−toe+200]とする。
s−t−toeは、足指脈波開始点からピークまでの時間と収縮期時間との比率であり、200が調整係数である。ピークの最高点が中点を超えるまで後方へ移動し、すなわちtmax−toe≧tch−toeである場合、ks−t−toeを補正する。ks−t−toeの値が大きいとき、足指の血管が拡張して、下肢血圧が低下することを示す。すなわち、ks−t−toe>0.8である場合、b=ks−t−toe−0.8とする。ks−t−toe≦0.8であり、bが適用できない場合、b=0とする。
(第5補正変数b
ステップS4で取得した補正変数はさらに第5補正変数bを含み、bの作用及び性質はbと同じであり、橈骨動脈血圧に対して下肢血圧が低下した場合にTを補正する。
s−m−toeは、足指脈波の収縮期での平均高さと最大高さhmax−toeとの比率であり、ks−m−toeが非常に大きいと、足指脈波ピークが広くて穏やかであることを示し、足指血管が拡張して、橈骨動脈に対して下肢血圧が低下することを表す。
足指血管が拡張しない場合、bが適用できない。すなわち、ks−m−toe<dである場合、b=0とする。d=0.67〜0.73であり、好ましくは0.7である。
足指血管が拡張しかつ脈波ピークの最高点が中点を超えるまで後方へ移動する場合、bは正の値とされる。すなわち、ks−m−toe≧dかつks−t−toe≧0.8である場合、b=ks−m−toe−dとする。
足指血管が拡張するが、脈波ピークの最高点の位置が中点を超えない場合、bは正の値を取りかつ半分にする。すなわち、ks−m−toe≧dかつks−t−toe<0.8である場合、b=(ks−m−toe−d)/2とする。
(第6補正変数b
ステップS4で取得した補正変数は、さらに、2つの脈波面積の相対変化を表し、足指血管が拡張して橈骨動脈血圧に対する下肢血圧が低下する場合にTを補正することに用いられる第6補正変数bを含む。bの適用範囲はb>0である。
s−m−toe−earは、足指脈波の収縮期での面積と耳脈波の収縮期での面積との比率であり、100は調整係数であり、ks−m−toe−earとkts−toe−earの作用及び性質は同じである。
足指脈波の面積が耳脈波の面積より小さい場合、足指血管が相対的に拡張せず、bが適用できない。すなわち、ks−m−toe−ear<1.0である場合、b=0とする。
第1の前提条件下、足指の面積が耳面積よりはるかに大きく、足指血管が多く拡張する場合、Cは定数1.08を最大値として用意される。すなわち、ks−m−toe−ear>1.08である場合、C=1.08とする。
耳脈波の形態が正常であると、bは最大補正値を取る。すなわち、t>220かつksd−m−0>0.88である場合、b=C−1.0とする。
耳脈波において非常に鋭い前へ傾斜した三角形が現れるかまたは波形が非常に狭いと、耳脈波の形態が深刻に異常に変化したことを表し、このとき、2つの脈波間の相対変化が拡大されるため、補正値を減少して使用する必要があり、bは最大補正値の1/3とされる。すなわち、t<160またはksd−m−0<0.80である場合、b=(C−1.0)×0.34とする。
耳脈波の形態の変異がそれほど深刻ではない場合、bは最大補正値の2/3とされる。すなわち、160<t≦220または0.80<ksd−m−0≦0.88である場合、b=(C−1.0)×0.67とする。
第2前提条件下、足指の面積が耳の面積より大きく、足指血管の相対拡張がそれほど深刻ではない場合、Cは正変数として用意される。すなわち、1.0≦ks−m−toe−ear≦1.08である場合、C=ks−m−toe−ear−1.0とする。
耳脈波の形態が正常であると、bは正変数を補正値とする。すなわち、t>220かつksd−m−0>0.88である場合、b=Cとする。
耳脈波の形態が深刻に異常に変化した場合、脈波間の相対変化が拡大されるため、補正値を減少して使用する必要があり、bは正変数の1/3とされる。すなわち、t≦160またはksd−m−0≦0.80である場合、b=C×0.34とする。
耳脈波の変異がそれほど深刻ではない場合、bは正変数の2/3とされる。すなわち160<t≦220または0.80<ksd−m−0≦0.88である場合、b=C×0.67とする。
(第7補正変数b
ステップS4で取得した補正変数はさらに第7補正変数bを含み、bの作用及び性質はbと同じであり、bは2つの脈波収縮期の幅(収縮期時間)の相対変化を表す。
ts−toe−ear=(ts−toe+825)/(t+825)。kts−toe−earは、足指の脈波において識別される心臓収縮期時間と耳脈波において識別される収縮期時間との比率であり、825は調整係数である。kts−toe−earが大きくなると、足指血管が拡張して、橈骨動脈血圧に対して下肢血圧が低下することを示す。
足指血管が相対的に拡張しない場合、bが適用できない。すなわち、kts−toe−ear<1.0である場合、b=0とする。
第1前提条件下、足指の血管が比較的大きく拡張する場合、Cは定数1.08を最大値として用意される。すなわち、kts−toe−ear>1.08である場合、C=1.08とする。
耳脈波の形態が正常である場合、bは最大補正値とされる。すなわち、t>220かつksd−m−0>0.88である場合、b=C−1.0とする。
耳脈波の形態が深刻に変化した場合、脈波間の相対変化が拡大されるため、補正値を減少して使用する必要があり、bは最大補正値の1/3とされる。すなわち、t<160またはksd−m−0<0.80である場合、b=(C−1.0)×0.34とする。
耳脈波の形態の変化が深刻ではない場合、bは最大補正値の2/3とされる。すなわち、160<t≦220または0.80<ksd−m−0≦0.88である場合、b=(C−1.0)×0.67とする。
第2前提条件下、足指の幅が耳の幅より大きく、足指血管の相対拡張がそれほど深刻ではない場合、Cは正変数として用意される。すなわち、1.0≦kts−toe−ear≦1.08である場合、C=kts−toe−ear−1.0とする。
耳脈波の形態が正常である場合、bは正変数を補正値とする。すなわち、t>220かつksd−m−0>0.88である場合、b=Cとする。
耳脈波の形態が深刻に変化した場合、bは正変数の1/3とされる。すなわち、t≦160またはksd−m−0≦0.80である場合、b=C×0.34とする。
耳脈波の形態の変化がそれほど深刻ではない場合、bは正変数の2/3とされる。すなわち、160<t≦220または0.80<ksd−m−0≦0.88である場合、b=(C−1.0)×0.67とする。
ステップS5における補正行列は下式である。
ここで、b=0である場合、該bが適用できない。ステップS6は、具体的には、再帰方式で選択されるものであって1つの最新変数を算出するたびに最も古い変数が1つ削除される8つの心周期での補正行列を連続的に取得し、8つの心周期の変数の平均値で呼吸変異による干渉を解消する。補正方法は下式である。
(実施形態2)
収縮期血圧に関する脈波伝播時間の補正方法は、以下のステップを含む。
S1)各心周期内の耳脈波をリアルタイムに検出して分析し、耳脈波における大動脈弁閉鎖点の高さhsd(すなわち、収縮期と拡張期の耳脈波に現れる境界部での高さ)、耳脈波の収縮期時間t(単位はミリ秒)、耳脈波の拡張期時間t(単位はミリ秒)、耳脈波の最大高さhmaxを含む耳脈波のデータを取得する。
S2)各心周期内の足指脈波をリアルタイムに検出して分析し、足指脈波の収縮期時間ts−toe(単位はミリ秒)、足指脈波の拡張期時間td−toe(単位はミリ秒)、足指脈波の最大高さhmax−toe、足指脈波の開始点からピーク中点までの時間tch−toe(単位はミリ秒)、足指脈波の開始点からピークの最高点までの時間tmax−toe(単位はミリ秒)を含む足指脈波のデータを取得する。ピーク中点は、ピークでの立ち上がりエッジの転換点と立ち下がりエッジの転換点との中点である。ピーク中点の定義について、YAN CHEN、CHANGYUN WEN、GUOCAI TAO、及びMIN BI著の『Continuous and Noninvasive Measurement of Systolic and Diastolic Blood Pressure by One Mathematical Model with the Same Model Parameters and Two Separate Pulse Wave Velocities』を参照することができる。
S3)収縮期血圧に関する脈波伝播時間Tを計算する。Tの定義については、YAN CHEN、CHANGYUN WEN、GUOCAI TAO、及びMIN BI著の『Continuous and Noninvasive Measurement of Systolic and Diastolic Blood Pressure by One Mathematical Model with the Same Model Parameters and Two Separate Pulse Wave Velocities』を参照することができる。hは耳脈波または足指脈波の縦軸方向での振幅である。
S4)同一心周期にステップS1、S2で取得したデータを用いて、該心周期での補正変数を算出する。
S5)ステップS4で取得した心周期での補正変数に基づき、該心周期での補正行列を算出する。
S6)複数の心周期での補正行列を連続的に取得し、ステップS3で取得したTを補正する。
この方法は、同一心周期内の耳脈波と足指脈波をリアルタイムに検出し、収縮期血圧に関する脈波伝播時間を計算し、そして脈波の形態学的特徴に基づき補正変数を抽出し、補正行列を取得し、上記脈波伝播時間の変異を補正することができ、補正済みの伝播時間を従来の数学モデルに用いて、臨床条件で各心周期の収縮期血圧を連続的に測定できる。
(第1補正変数a
ステップS4で取得した補正変数は、低血圧状態での収縮期血圧に関する伝播時間Tのクラス2の変化の補正に用いられる第1補正変数aを含み、aの適用範囲はa>0であり、aが大きいほど血圧が低いことを示す。
sd−m−0は、hsdと耳脈波収縮期の平均高さとの比を表す。一部の症例は、低血圧状態で、脈波ピークが前傾した三角形として現れ、hsdが大きく減少し、ksd−m−0が小さくなるが、これは、大動脈収縮期の最後部分の波形が大きく低下することを示し、脈波伝播を促進する持続動力が不足であり、伝播時間が長くなると推定される。このような状態では、拡張期情報が不安定になり、使用に不適である。
=76〜84であり、好ましくは80である。d1−2=104〜112であり、好ましくは108である。
脈波伝播を促進する持続動力が不足である場合、伝播時間Tが長くなり、aで補正する必要がある。すなわち、d≦ksd−m−0≦d1−2である場合、a=(d1−2−ksd−m−0)×0.50とする。
脈波伝播を促進する持続動力が深刻に不足になる場合、伝播時間Tがはるかに長くなり、aの上限値で補正する必要がある。すなわち、ksd−m−0<dである場合、a=28×0.50とする。
脈波伝播を促進する持続動力が十分である場合、Tを補正する必要がなく、aが適用できない。すなわち、ksd−m−0>d1−2である場合、a=0とする。
(第2補正変数a
ステップS4で取得した補正変数は、さらに、高血圧状態または正常血圧状態から高血圧状態に変わる状態で収縮期血圧に関する伝播時間Tのクラス2の変化の補正に用いられる第2補正変数aを含み、aの適用範囲はa>0であり、aが大きいほど拡張期血圧が高いことを示す。
sd−m−tsはhsdと耳脈波拡張期のt−2t部分での平均高さとの比を表し、脈波拡張期の変異を判断することに用いられる。たとえば、胸部や腹部手術中におけるフックの使用は大動脈受力の変化を引き起こし、耳脈波拡張期の波形を低くして、ksd−m−tsを増大させる。
sd−m−2は、hsdと耳脈波の0−2t部分での平均高さとの比率を表し、収縮期と一部の拡張期の波形情報を含み、主に、気管内挿管による心拍数と血圧の上昇のような高血圧状態に用いられる。正常血圧状態から高血圧状態に変わる過程では、耳脈波は正三角形または後方へ傾斜した三角形を呈し、hsdははるかに上昇し、ksd−m−2は次第に大きくなる。高血圧状態では、耳脈波は正三角形または後方へ傾斜した三角形を呈し、hsdははるかに上昇し、ksd−m−2はかなり大きくなる。上記二つの三角形波形の頂部(最も高い血圧)は維持時間が短く、最も高い血圧に対応する促進する動力が不足になり、伝播時間Tが長くなる。
|ksd−m−0−ksd−m−ts|≧40かつ(ksd−m−0+ksd−m−ts)/2≧ksd−m−2
ならば、ksd−m=2×ksd−m−2−(ksd−m−0+ksd−m−ts)/2、
さもなくば、ksd−m=ksd−m−2
これは、耳脈波拡張期の波形が変化する場合、たとえば、胸部や腹部手術中におけるフックの使用により大動脈の受力が変化した結果、脈波拡張期の形態が大きく変化する場合、ksd−mを補正し、そうでない場合はksd−m=ksd−m−2とするということである。
=1.17〜1.27であり、好ましくは1.22である。
sd−m>(d+(age−14)/15/100)(ageは年齢で、age≧14)であり、最も高い血圧に対応する持続動力が不足になり、それに対応して伝播時間Tが長くなり、aで補正する必要がある場合、a=(ksd−m−(d+(age−14)/15/100))とする。
sd−m≦(d+(age−14)/15/100)であり、脈波ピークの部分が緩やかで、最も高い血圧に対応する持続動力が十分であり、aが適用できない場合、a=0とする。
(第3補正変数a
ステップS3で取得した補正変数は、さらに、血液量が変化しまたはセンサの配置部位の体温が変化した状態でTを補正するための第3補正変数aを含む。
d−m−tdは、耳脈波拡張期の平均高さと最大高さhmaxとの比率である。手術前に患者が断食しかつ飲水量を減少させると、血液量が低下して、kd−m−tdが減少し、脈波伝播時間が長くなり、手術中に輸血・輸液により血液量が増加すると、kd−m−tdが大きくなり、伝播時間が短くなる。
sd−m−ts≦d3−2である場合、耳脈波の拡張早期での波形が正常な範囲を超えるまで上昇することを表し、kd−m−tdを修正する必要があり、修正結果をkd−m−td−1として表記する。
d−m−td−1=kd−m−td−(d3−2−ksd−m−ts)×75/100であり、kd−m−td≦dである場合、耳脈波が干渉を受けたと判断でき、この場合kd−m−td−1=dとする。d=0.02〜0.14であり、好ましくは0.08である。d3−2=1.21〜1.31であり、好ましくは1.26である。
d−m−td−toeは、足指脈波の拡張期での平均高さと最大高さhmax−toeとの比率であり、ts−toeは足指脈波において識別された心臓収縮期時間を表し、td−toeは足指脈波において識別された拡張期時間を表す。kd−m−td−toe≦dである場合、kd−m−td−toe=dとする。kd−m−td−toeとkd−m−tdの作用及び性質は同じである。
d−m−a=(kd−m−td−1+kd−m−td−toe)/2として、耳と足指の脈波の性質が同じである2つの変数を合併し、その平均値を脈波伝播時間Tを補正する変数としたとき、脈波拡張期の波形に変化が発生した場合に、kd−m−aを補正する。
|ksd−m−0−ksd−m−ts|≧40かつ(ksd−m−0+ksd−m−ts)/2≧ksd−m−2かつksd−m−ts≧d3−2である場合、
d−ma−a=(kd−m−td−1+kd−m−td−toe+(ksd−m−0+ksd−m−ts)/2−ksd−m−2)/2とする。
血液量が正常であるとともに、センサの配置部位の体温が正常である状態には、aが適用できない。すなわち、c<kd−m−a<cである場合、a=0とする。c=(d+(age−14)/8)/100、d=23〜35、好ましくは29、c=(d+(age−14)/8)/100、d=27〜39、好ましくは33である。
血圧が非常に低くまたは非常に高い状態では、拡張期情報が不安定になり、aが適用できない。すなわち、ksd−m−0<dまたはksd−m−2>dである場合、a=0とする。d=0.97〜1.03であり、好ましくは1.00である。d=1.52〜1.58であり、好ましくは1.55である。
正常血圧状態では、血液量が低下しまたはセンサの配置部位の体温が低下する場合、aは正の値の67%とされる。すなわち、ksd−m−0≧d+0.10かつksd−m−2≦dかつkd−m−a≦cである場合、a=(c−kd−m−a)×67/100とする。d=1.42〜1.48であり、好ましくは1.45である。
低血圧または高血圧状態では、血液量が低下しまたはセンサの配置部位の体温が低下する場合、aは正常血圧状態の数値の50%とされる。すなわち、下式を満たす場合、a=(c−kd−m−a)×50/100とする。
正常血圧状態では、血液量が増加しまたはセンサの配置部位の体温が上昇する場合、aは負の値の62%とされる。すなわち、ksd−m−0≧d+0.10かつksd−m−2≦dかつkd−m−a≧cである場合、a=(c−kd−m−a)×62/100とする。
低血圧または高血圧状態では、血液量が増加しまたはセンサの配置部位の体温が上昇する場合、aは正常血圧状態の負の値の45%とされる。すなわち、下式を満たす場合、a=(c−kd−m−a)×45/100とする。
(第4補正変数a
ステップS4で取得した補正変数は、さらに、末梢血管拡張により下肢血圧(橈骨動脈血圧に対する)の低下をもたらす場合にTを補正する第4補正変数aを含み、aの適用範囲はa>0であり、aが大きいほど、橈骨動脈血圧に対する下肢血圧の低下が大きいことを示す。
末梢血管の収縮と拡張は、足指脈波のピーク位置の時間軸における前後移動を引き起こす。tmax−toe≧tch−toeである場合、ks−t−toe=[tmax−toe+tch−toe+400]/[(ts−toe+200)×2]とし、そうでない場合、ks−t−toe=[tmax−toe+200]/[ts−toe+200]とする。
s−t−toeは、足指脈波開始点からピークまでの時間と収縮期時間との比率であり、200が調整係数である。ピークの最高点が中点を超えるまで後方へ移動し、すなわちtmax−toe≧tch−toeである場合、ks−t−toeを補正する。ks−t−toeの値が大きいとき、足指の血管が拡張して、下肢血圧が低下することを示す。すなわち、ks−t−toe>0.8である場合、a=ks−t−toe−0.8とする。ks−t−toe≦0.8であり、aが適用できない場合、a=0とする。
(第5補正変数a
ステップS4で取得した補正変数はさらに第5補正変数aを含み、aの作用及び性質はaと同じであり、橈骨動脈血圧に対して下肢血圧が低下した場合にTを補正する。
s−m−toeは、足指脈波の収縮期での平均高さと最大高さhmax−toeとの比率であり、ks−m−toeが非常に大きいと、足指脈波ピークが広くて穏やかであることを示し、足指血管が拡張して、橈骨動脈に対して下肢血圧が低下することを表す。
足指血管が拡張しない場合、aが適用できない。すなわち、ks−m−toe<dである場合、a=0とする。d=0.67〜0.73であり、好ましくは0.7である。
足指血管が拡張しかつ脈波ピークの最高点が中点を超えるまで後方へ移動する場合、aは正の値とされる。すなわち、ks−m−toe≧dかつks−t−toe≧0.8である場合、b=ks−m−toe−dとする。
足指血管が拡張するが、脈波ピークの最高点の位置が中点を超えない場合、aは正の値を取りかつ半分にする。すなわち、ks−m−toe≧dかつks−t−toe<0.8である場合、a=(ks−m−toe−d)/2とする。
(第6補正変数a
ステップS4で取得した補正変数は、さらに、2つの脈波面積の相対変化を表し、足指血管が拡張して橈骨動脈血圧に対する下肢血圧が低下する場合にTを補正することに用いられる第6補正変数aを含む。aの適用範囲はa>0である。
s−m−toe−earは、足指脈波の収縮期での面積と耳脈波の収縮期での面積との比率であり、100は調整係数であり、ks−m−toe−earとkts−toe−earの作用及び性質は同じである。
足指脈波の面積が耳脈波の面積より小さい場合、足指血管が相対的に拡張せず、aが適用できない。すなわち、ks−m−toe−ear<1.0である場合、a=0とする。
第1の前提条件下、足指の面積が耳面積よりはるかに大きく、足指血管が多く拡張する場合、Cは定数1.08を最大値として用意される。すなわち、ks−m−toe−ear>1.08である場合、C=1.08とする。
耳脈波の形態が正常であると、aは最大補正値を取る。すなわち、t>220かつksd−m−0>0.88である場合、a=C−1.0とする。
耳脈波において非常に鋭い前へ傾斜した三角形が現れるかまたは波形が非常に狭いと、耳脈波の形態が深刻に異常に変化したことを表し、このとき、2つの脈波間の相対変化が拡大されるため、補正値を減少して使用する必要があり、aは最大補正値の1/3とされる。すなわち、t<160またはksd−m−0<0.80である場合、a=(C−1.0)×0.34とする。
耳脈波の形態の変異がそれほど深刻ではない場合、bは最大補正値の2/3とされる。すなわち、160<t≦220または0.80<ksd−m−0≦0.88である場合、a=(C−1.0)×0.67とする。
第2前提条件下、足指の面積が耳の面積より大きく、足指血管の相対拡張がそれほど深刻ではない場合、Cは正変数として用意される。すなわち、1.0≦ks−m−toe−ear≦1.08である場合、C=ks−m−toe−ear−1.0とする。
耳脈波の形態が正常であると、aは正変数を補正値とする。すなわち、t>220かつksd−m−0>0.88である場合、a=Cとする。
耳脈波の形態が深刻に異常に変化した場合、脈波間の相対変化が拡大されるため、補正値を減少して使用する必要があり、aは正変数の1/3とされる。すなわち、t≦160またはksd−m−0≦0.80である場合、a=C×0.34とする。
耳脈波の変異がそれほど深刻ではない場合、aは正変数の2/3とされる。すなわち160<t≦220または0.80<ksd−m−0≦0.88である場合、a=C×0.67とする。
(第7補正変数a
ステップS4で取得した補正変数はさらに第7補正変数aを含み、aの作用及び性質はaと同じであり、aは2つの脈波収縮期の幅(収縮期時間)の相対変化を表す。
ts−toe−ear=(ts−toe+825)/(t+825)。kts−toe−earは、足指の脈波において識別される心臓収縮期時間と耳脈波において識別される収縮期時間との比率であり、825は調整係数である。kts−toe−earが大きくなると、足指血管が拡張して、橈骨動脈血圧に対して下肢血圧が低下することを示す。
足指血管が相対的に拡張しない場合、aが適用できない。すなわち、kts−toe−ear<1.0である場合、a=0とする。
第1前提条件下、足指の血管が比較的大きく拡張する場合、Cは定数1.08を最大値として用意される。すなわち、kts−toe−ear>1.08である場合、C=1.08とする。
耳脈波の形態が正常である場合、aは最大補正値とされる。すなわち、t>220かつksd−m−0>0.88である場合、a=C−1.0である。
耳脈波の形態が深刻に変化した場合、脈波間の相対変化が拡大されるため、補正値を減少して使用する必要があり、aは最大補正値の1/3とされる。すなわち、t<160またはksd−m−0<0.80である場合、a=(C−1.0)×0.34とする。
耳脈波の形態の変化が深刻ではない場合、aは最大補正値の2/3とされる。すなわち、160<t≦220または0.80<ksd−m−0≦0.88である場合、a=(C−1.0)×0.67とする。
第2前提条件下、足指の幅が耳の幅より大きく、足指血管の相対拡張がそれほど深刻ではない場合、Cは正変数として用意される。すなわち、1.0≦kts−toe−ear≦1.08である場合、C=kts−toe−ear−1.0とする。
耳脈波の形態が正常である場合、aは正変数を補正値とする。すなわち、t>220かつksd−m−0>0.88である場合、a=Cとする。
耳脈波の形態が深刻に変化した場合、aは正変数の1/3とされる。すなわち、t≦160またはksd−m−0≦0.80である場合、a=C×0.34とする。
耳脈波の形態の変化がそれほど深刻ではない場合、aは正変数の2/3とされる。すなわち、160<t≦220または0.80<ksd−m−0≦0.88である場合、a=C×0.67とする。
ステップS5における補正行列は下式である。
ここで、a=0である場合、該aが適用できない。ステップS6は、具体的には、再帰方式で選択されるものであって1つの最新変数を算出するたびに最も旧い変数が1つ削除された8つの心周期での補正行列を連続的に取得し、8つの心周期の変数の平均値で呼吸変異による干渉を解消する。補正方法は下式である。
最後に、以上の各実施形態は本発明の技術案を説明するためのものに過ぎず、限定するものではない。本発明を上述の実施形態を参照して詳細に説明したが、当業者であれば、上述の各実施形態に記載の技術案を修正したりまたはその一部もしくは全ての技術的特徴に対して同等置換を行ったりすることができ、これら修正または置換によって対応する技術案の本質が本発明の各実施形態の技術案の範囲から脱逸することがなく、これら修正または置換は本発明の特許請求の範囲及び明細書の範囲に属すべきであることを理解できるだろう。
[付記]
[付記1]
各心周期内の耳脈波をリアルタイムに検出して分析し、耳脈波における大動脈弁閉鎖点の高さhsd、耳脈波の収縮期時間t(単位はミリ秒)、耳脈波の拡張期時間t(単位はミリ秒)、耳脈波の最大高さhmaxを含む耳脈波のデータを取得するステップS1と、
各心周期内の足指脈波をリアルタイムに検出して分析し、足指脈波の収縮期時間ts−toe(単位はミリ秒)、足指脈波の拡張期時間td−toe(単位はミリ秒)、足指脈波の最大高さhmax−toe、足指脈波の開始点から、ピークでの立ち上がりエッジの転換点と立ち下がりエッジの転換点との中点であるピーク中点までの時間tch−toe(単位はミリ秒)、足指脈波の開始点からピークの最高点までの時間tmax−toe(単位はミリ秒)を含む足指脈波のデータを取得するステップS2と、
拡張期血圧に関する脈波伝播時間Tを計算するステップS3であって、該Tが耳脈波の開始点から足指脈波の開始点までの時間差であり、hが耳脈波または足指脈波の縦軸方向での振幅であるステップS3と、
同一心周期にステップS1、S2で取得したデータを用いて、該心周期での補正変数を算出するステップS4と、
ステップS4で取得した心周期での補正変数に基づき、該心周期での補正行列を算出するステップS5と、
複数の心周期での補正行列を連続的に取得し、ステップS3で取得したTを補正するステップS6と、
を含む、
ことを特徴とする、拡張期血圧に関する脈波伝播時間の補正方法。
[付記2]
前記ステップS5における補正行列は、下式である、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記3]
前記ステップS6は、具体的には、8つの心周期での補正行列を連続的に取得し、
とする方法で補正を行うことである、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記4]
第1補正変数bは、下式に従って算出される、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記5]
第2補正変数bは、下式に従って算出される、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記6]
第3補正変数bは、下式に従って算出される、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記7]
第4補正変数bは、下式に従って算出される、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記8]
第5補正変数bは、下式に従って算出される、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記9]
第6補正変数bは、下式に従って算出される、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記10]
第7補正変数bは、下式に従って算出される、
付記1に記載の拡張期血圧に関する脈波伝播時間の補正方法。
[付記11]
各心周期内の耳脈波をリアルタイムに検出して分析し、耳脈波における大動脈弁閉鎖点の高さhsd、耳脈波の収縮期時間t(単位はミリ秒)、耳脈波の拡張期時間t(単位はミリ秒)、耳脈波の最大高さhmaxを含む耳脈波のデータを取得するステップS1と、
各心周期内の足指脈波をリアルタイムに検出して分析し、足指脈波の収縮期時間ts−toe(単位はミリ秒)、足指脈波の拡張期時間td−toe(単位はミリ秒)、足指脈波の最大高さhmax−toe、足指脈波の開始点から、ピークでの立ち上がりエッジの転換点と立ち下がりエッジの転換点との中点であるピーク中点までの時間tch−toe(単位はミリ秒)、足指脈波の開始点からピークの最高点までの時間tmax−toe(単位はミリ秒)を含む足指脈波のデータを取得するステップS2と、
収縮期血圧に関する脈波伝播時間Tを計算するステップS3であって、該Tが耳脈波の開始点から足指脈波の開始点までの時間差であり、hが耳脈波または足指脈波の縦軸方向での振幅であるステップS3と、
同一心周期にステップS1、S2で取得したデータを用いて、該心周期での補正変数を算出するステップS4と、
ステップS4で取得した心周期での補正変数に基づき、該心周期での補正行列を算出するステップS5と、
複数の心周期での補正行列を連続的に取得し、ステップS3で取得したTを補正するステップS6と、
を含む、
ことを特徴とする、収縮期血圧に関する脈波伝播時間の補正方法。
[付記12]
前記ステップS5における補正行列は、下式である、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。
[付記13]
前記ステップS6は、具体的に、8つの心周期での補正行列を連続的に取得し、
とする方法で補正を行うことである、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。
[付記14]
第1補正変数aは、下式に従って算出される、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。
[付記15]
第2補正変数aは、下式に従って算出される、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。
[付記16]
第3補正変数aは、下式に従って算出される、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。
[付記17]
第4補正変数aは、下式に従って算出される、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。
[付記18]
第5補正変数aは、下式に従って算出される、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。
[付記19]
第6補正変数aは、下式に従って算出される、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。
[付記20]
第7補正変数aは、下式に従って算出される、
付記11に記載の収縮期血圧に関する脈波伝播時間の補正方法。

Claims (9)

  1. 各心周期内の耳脈波をリアルタイムに検出して分析し、耳脈波における大動脈弁閉鎖点の高さhsd、耳脈波の収縮期時間t(単位はミリ秒)、耳脈波の拡張期時間t(単位はミリ秒)、耳脈波の最大高さhmaxを含む耳脈波のデータを取得するステップS1と、
    各心周期内の足指脈波をリアルタイムに検出して分析し、足指脈波の収縮期時間ts−toe(単位はミリ秒)、足指脈波の拡張期時間td−toe(単位はミリ秒)、足指脈波の最大高さhmax−toe、足指脈波の開始点から、ピークでの立ち上がりエッジの転換点と立ち下がりエッジの転換点との中点であるピーク中点までの時間tch−toe(単位はミリ秒)、足指脈波の開始点からピークの最高点までの時間tmax−toe(単位はミリ秒)を含む足指脈波のデータを取得するステップS2と、
    拡張期血圧に関する脈波伝播時間T又は収縮期血圧に関する脈波伝播時間Tを計算するステップS3であって、該Tが耳脈波の開始点から足指脈波の開始点までの時間差であり、該Tが耳脈波における大動脈弁閉鎖点から足指脈波における大動脈弁閉鎖点までの時間差であり、hが耳脈波または足指脈波の縦軸方向での振幅であるステップS3と、
    同一心周期にステップS1、S2で取得したデータを用いて、該心周期での補正変数を算出するステップS4と、
    ステップS4で取得した心周期での補正変数 〜b 又はa 〜a に基づき、該心周期での補正行列を算出するステップS5と、
    複数の心周期での補正行列を連続的に取得し、ステップS3で取得したT又はTを補正するステップS6と、
    を含
    前記ステップS4で計算した補正変数は、低血圧状態でT 又はT を補正するb とa と、高血圧状態でT を補正するb と、高血圧状態または正常血圧から高血圧状態に変わる状態でT を補正するa と、血液量が変化しまたはセンサの配置部位の体温が変化した状態でT 又はT を補正するb とa と、橈骨動脈血圧に対して下肢血圧が低下した場合にT 又はT を補正するb 、b 、b 、b とa 、a 、a 、a とを含み、
    前記ステップS3でT を計算するとき、前記ステップS5における補正行列は、下式であり、
    前記ステップS3でT を計算するとき、前記ステップS5における補正行列は、下式である、
    ことを特徴とする、動脈血圧に関する脈波伝播時間の補正方法。
  2. 前記ステップS6は、具体的には、8つの心周期での補正行列を連続的に取得し、
    前記ステップS3でTを計算するときは、
    とする方法で補正を行うことであり、
    前記ステップS3でTを計算するときは、
    とする方法で補正を行うことである、
    請求項に記載の動脈血圧に関する脈波伝播時間の補正方法。
  3. 前記ステップS3でTを計算するとき、第1補正変数bは、下式に従って算出され、
    前記ステップS3でTを計算するとき、第1補正変数aは、下式に従って算出される、
    請求項1又は2に記載の動脈血圧に関する脈波伝播時間の補正方法。
  4. 前記ステップS3でTを計算するとき、第2補正変数bは、下式に従って算出され、
    前記ステップS3でTを計算するとき、第2補正変数aは、下式に従って算出される、
    請求項1から3のいずれか1項に記載の動脈血圧に関する脈波伝播時間の補正方法。
  5. 前記ステップS3でTを計算するとき、第3補正変数bは、下式に従って算出され、
    前記ステップS3でTを計算するとき、第3補正変数aは、下式に従って算出される、
    請求項1から4のいずれか1項に記載の動脈血圧に関する脈波伝播時間の補正方法。
  6. 前記ステップS3でTを計算するとき、第4補正変数bは、下式に従って算出され、
    前記ステップS3でTを計算するとき、第4補正変数aは、下式に従って算出される、
    請求項1から5のいずれか1項に記載の動脈血圧に関する脈波伝播時間の補正方法。
  7. 前記ステップS3でTを計算するとき、第5補正変数bは、下式に従って算出され、
    前記ステップS3でTを計算するとき、第5補正変数aは、下式に従って算出される、
    請求項1から6のいずれか1項に記載の動脈血圧に関する脈波伝播時間の補正方法。
  8. 前記ステップS3でTを計算するとき、第6補正変数bは、下式に従って算出され、
    前記ステップS3でTを計算するとき、第6補正変数aは、下式に従って算出される、
    請求項1から7のいずれか1項に記載の動脈血圧に関する脈波伝播時間の補正方法。
  9. 前記ステップS3でTを計算するとき、第7補正変数bは、下式に従って算出され、
    前記ステップS3でTを計算するとき、第7補正変数aは、下式に従って算出される、
    請求項1から8のいずれか1項に記載の動脈血圧に関する脈波伝播時間の補正方法。
JP2019541843A 2016-11-22 2017-11-20 動脈血圧に関する脈波伝播時間の補正方法 Active JP6736110B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
CN201611045054.5A CN106580303B (zh) 2016-11-22 2016-11-22 与收缩压相关的脉搏波传播时间的校正方法
CN201611046184.0A CN106377238B (zh) 2016-11-22 2016-11-22 与舒张压相关的脉搏波传播时间的校正方法
CN201611046184.0 2016-11-22
CN201611045054.5 2016-11-22
PCT/CN2017/111799 WO2018095291A1 (zh) 2016-11-22 2017-11-20 与舒张压、收缩压相关的脉搏波传播时间的校正方法

Publications (3)

Publication Number Publication Date
JP2019530560A JP2019530560A (ja) 2019-10-24
JP2019530560A5 JP2019530560A5 (ja) 2019-12-05
JP6736110B2 true JP6736110B2 (ja) 2020-08-05

Family

ID=62194728

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019541843A Active JP6736110B2 (ja) 2016-11-22 2017-11-20 動脈血圧に関する脈波伝播時間の補正方法

Country Status (4)

Country Link
US (1) US20190246919A1 (ja)
EP (1) EP3545832B1 (ja)
JP (1) JP6736110B2 (ja)
WO (2) WO2018095083A1 (ja)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108272446B (zh) * 2018-01-30 2021-03-26 浙江大学 无创连续血压测量系统及其校准方法
CN110853755A (zh) * 2019-11-07 2020-02-28 山西中医药大学 一种更年期女性脉象特点的研究方法
CN114305358B (zh) * 2021-02-24 2023-04-14 心永(北京)科技有限公司 血压测量模型的标定方法、装置、计算机设备和存储介质
CN113040738B (zh) * 2021-03-29 2023-03-24 南京邮电大学 血压检测装置
CN113598724A (zh) * 2021-09-03 2021-11-05 上海市高血压研究所 用于评估心血管自主神经功能的装置
CN113598737A (zh) * 2021-09-08 2021-11-05 上海中医药大学 一种基于特征融合的血压预测方法
TW202348196A (zh) * 2022-01-28 2023-12-16 石明正 脈診的脈搏波量測裝置及其使用方法
CN114391822B (zh) * 2022-02-28 2024-12-10 维沃移动通信有限公司 数据修正方法、数据修正装置、电子设备和可读存储介质
CN116172528A (zh) * 2022-08-17 2023-05-30 深圳市奥极健康科技有限公司 一种基于误差补偿的脉搏波连续血压测量系统及方法
CN119856916B (zh) * 2025-03-25 2025-05-13 沈阳恒德医疗器械研发有限公司 一种心血管检测方法及系统

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3926800A1 (de) 1989-08-14 1991-02-21 Hoechst Ag Wasserdispergierbare granulate zur anwendung im pflanzenschutz
EP0443267A1 (en) * 1990-02-23 1991-08-28 Sentinel Monitoring, Inc. Method and apparatus for continuous non-invasive blood pressure monitoring
JP3318727B2 (ja) 1994-06-06 2002-08-26 日本光電工業株式会社 脈波伝播時間方式血圧計
CN1127939C (zh) 1995-07-28 2003-11-19 北京新兴生物医学工程研究发展中心 一种利用脉波测量动脉血压的方法
JP3778655B2 (ja) * 1996-07-01 2006-05-24 コーリンメディカルテクノロジー株式会社 血圧監視装置
US5865755A (en) 1996-10-11 1999-02-02 Dxtek, Inc. Method and apparatus for non-invasive, cuffless, continuous blood pressure determination
JPH10295657A (ja) * 1997-04-24 1998-11-10 Matsushita Electric Ind Co Ltd 血圧測定装置
JPH10314129A (ja) * 1997-05-23 1998-12-02 Nippon Colin Co Ltd 脈波伝播速度情報測定機能付き自動血圧測定装置
JP3330079B2 (ja) * 1998-05-12 2002-09-30 日本コーリン株式会社 非観血連続血圧推定装置
JP3781956B2 (ja) * 2000-08-28 2006-06-07 株式会社タニタ 脈波検出装置
CN100560019C (zh) 2007-01-26 2009-11-18 香港中文大学 脉搏波传输时间法测量动脉血压的初始校准装置
JP4971041B2 (ja) * 2007-06-11 2012-07-11 株式会社デンソー 血圧測定装置及びプログラム並びに記録媒体
US8602997B2 (en) 2007-06-12 2013-12-10 Sotera Wireless, Inc. Body-worn system for measuring continuous non-invasive blood pressure (cNIBP)
CN102258364A (zh) * 2010-05-31 2011-11-30 深圳市盛力康实业发展有限公司 一种脉搏波识别方法、系统及动脉功能检测仪器
CN102008296B (zh) * 2010-12-24 2013-09-04 吉林大学 基于脉搏波信号和心电信号测量动脉血压装置及测量方法
WO2012092303A1 (en) 2010-12-28 2012-07-05 Sotera Wireless, Inc. Body-worn system for continous, noninvasive measurement of cardiac output, stroke volume, cardiac power, and blood pressure
CN103385702B (zh) * 2013-07-26 2015-08-26 中国科学院深圳先进技术研究院 一种无创血压连续检测装置及方法
CN104382571B (zh) * 2014-10-28 2017-06-16 深圳市维亿魄科技有限公司 一种基于桡动脉脉搏波传导时间的测量血压方法及装置
CN106580303B (zh) * 2016-11-22 2018-03-06 浙江脉联医疗设备有限公司 与收缩压相关的脉搏波传播时间的校正方法
CN106377238B (zh) * 2016-11-22 2018-03-06 浙江脉联医疗设备有限公司 与舒张压相关的脉搏波传播时间的校正方法

Also Published As

Publication number Publication date
WO2018095083A1 (zh) 2018-05-31
WO2018095291A1 (zh) 2018-05-31
EP3545832B1 (en) 2024-03-20
EP3545832A1 (en) 2019-10-02
US20190246919A1 (en) 2019-08-15
EP3545832A4 (en) 2020-08-05
JP2019530560A (ja) 2019-10-24

Similar Documents

Publication Publication Date Title
JP6736110B2 (ja) 動脈血圧に関する脈波伝播時間の補正方法
Kim et al. Ballistocardiogram-based approach to cuffless blood pressure monitoring: Proof of concept and potential challenges
Chandrasekhar et al. PPG sensor contact pressure should be taken into account for cuff-less blood pressure measurement
Peter et al. A review of methods for non-invasive and continuous blood pressure monitoring: Pulse transit time method is promising?
Yoon et al. Cuff-less blood pressure estimation using pulse waveform analysis and pulse arrival time
EP3030145B1 (en) Monitoring system and method for monitoring the hemodynamic status of a subject
EP3551060B1 (en) Pulse wave velocity determination, for example for blood pressure monitoring
EP3551059B1 (en) An apparatus and method for determining a calibration parameter for a blood pressure measurement device
WO2008094598A2 (en) Method and system for determining cardiac function
EP3220812A2 (en) Pulse wave velocity, arterial compliance, and blood pressure
JP2011520502A (ja) 血圧カフを使って得た中心圧力波形の推定方法
WO2016179425A1 (en) A device and method for non-invasive left ventricular end diastolic pressure (lvedp) measurement
CN106618504A (zh) 一种具备自适应校准的高准确度连续无创血压测量装置
KR20190048878A (ko) 광학 센서를 이용한 혈압 측정 방법 및 장치
CN118175958A (zh) 用于校准血压替代以用于监测对象的血压的设备、系统和方法
CN106580303B (zh) 与收缩压相关的脉搏波传播时间的校正方法
CN117064354A (zh) 无袖带连续血压测量方法及系统
CN118175957A (zh) 用于校准血压替代以用于监测对象的血压的设备、系统和方法
Liu et al. Arteries stiffen with age, but can retain an ability to become more elastic with applied external cuff pressure
CN106377238B (zh) 与舒张压相关的脉搏波传播时间的校正方法
US20230293035A1 (en) Cardiac Cycle Analysis Method and System
Savkar et al. Study on techniques involved in tourniqueteless blood pressure measurement using PPG
Kaushik et al. A comprehensive study on blood pressure measurement techniques
Gupta Blood Pressure Monitoring
US20220378307A1 (en) Method for correcting pulse wavetransit time associated with diastolic blood pressure or systolic blood pressure

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190412

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190412

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200225

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200221

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200522

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: 20200623

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200708

R150 Certificate of patent or registration of utility model

Ref document number: 6736110

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250