JP3293621B2 - 血液拍出量を決定するための装置 - Google Patents

血液拍出量を決定するための装置

Info

Publication number
JP3293621B2
JP3293621B2 JP50509092A JP50509092A JP3293621B2 JP 3293621 B2 JP3293621 B2 JP 3293621B2 JP 50509092 A JP50509092 A JP 50509092A JP 50509092 A JP50509092 A JP 50509092A JP 3293621 B2 JP3293621 B2 JP 3293621B2
Authority
JP
Japan
Prior art keywords
pressure
aorta
compliance
relationship
flow
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.)
Expired - Lifetime
Application number
JP50509092A
Other languages
English (en)
Other versions
JPH06505651A (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
Application filed by ネーデルランドセ・オルガニザテイエ・フール・テゲパスト−ナトウールベテンシヤツペリーク・オンデルツエク・テイエヌオー filed Critical ネーデルランドセ・オルガニザテイエ・フール・テゲパスト−ナトウールベテンシヤツペリーク・オンデルツエク・テイエヌオー
Publication of JPH06505651A publication Critical patent/JPH06505651A/ja
Application granted granted Critical
Publication of JP3293621B2 publication Critical patent/JP3293621B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

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 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/026Measuring blood flow
    • A61B5/029Measuring or recording blood output from the heart, e.g. minute volume

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Hematology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Physiology (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

【発明の詳細な説明】 本発明は、人の心臓の血液拍出量を決定するための装
置に関する。
血流内に低温液体を注入しかつ規則的な間隔でその下
流を測定する熱希釈(thermodilution)判定の支援によ
り1回拍出量及び心臓血液拍出量を算出することが実際
上知られている。
別の公知の方法においては、人の大動脈において測定
された圧力信号p(t)より1回拍出量Vs及び心臓血液
拍出量又は1分間心臓拍出量Qを決定するために、脈拍
輪郭法又は処方が使用される。これに関連して、1回拍
出量は、1収縮又は1拍動に心臓から排出される血液の
拍出量である。典型的な量は例えば70cm3であるが、こ
の量は拍動毎に変化しうる。心臓血液拍出量は単位時間
1分間における心臓が送り出す血液の量である。これに
関連して、典型的な値は5リットル/分であり、この量
は典型的には2から30リットル/分の間で変わりうる。
過去においては、測定された大動脈血圧信号p(t)
と平均大動脈流量q(t)との間に比較的簡単な関係が
仮定された。これに関連して、圧力は比較的容易にうま
く測定できるが流量は測定できないので、出発点は圧力
信号である。実際は、かかる相互関係又は「モデル」
は、最初に仮定されたよりも相当に複雑であることが見
いだされている。特に、このことは、このいわゆる脈拍
輪郭計算が重篤患者の場合に失敗の生ずる病院において
発生する。
特に、普通、大動脈を信号コンプライアンスとして考
えるいわゆる「ウインドケセル(Windkessel)」又はエ
ア・レシーバー・モデルが考えられる(図1b)。この場
合、ウインドケセル−液体及び気体で部分的に満たされ
た容器−であって、流れのサージを吸収し気泡中のサー
ジを緩衝でき、かつ拍動ポンプの下流に組み込まれたウ
インドケセルが意味される。心臓ポンプより排出された
血液の量は大動脈又はウインドケセルに大部分受け入れ
られ、一部は大動脈から分岐した種々の器官の末梢欠陥
部を通って流れ去る。心臓ポンプが血液を排出しない時
期、即ち格調期Tdには、ウインドケセルからの排出流量
が末梢血管部への供給を続ける。このとき、大動脈は気
体を少しも含まずかつ大動脈壁は弾性的であり、この弾
性的な血管壁ウインドケセルと同じ機能を果たすであろ
う。
この概念は無視された種々の要因を含む。第1に、大
動脈は幅よりも長さが大きい。実際上、心臓により作ら
れる圧力波は、大動脈の端部に達するのに0.1から0.3秒
の時間を要し、次いで心臓に戻るためにほぼ同じ時間が
必要であり、一方、心臓の全排出期(収縮時又は収縮
期、Ts)は0.2から0.4秒のみ継続する。排出の開始時に
は、排出流の緩衝用としてコンプライアンスの小部分だ
けが利用可能である。しかし、更により重要なことは、
圧力液の伝播又は移動時間の2倍後に、圧力液振幅が最
初の2倍にはなるが拍出量の2倍を利用できることであ
る。更に、一般に行われる第2の無視項目はウインドケ
セルが線形であることであり、これは圧力を2倍にする
ことにより蓄積される1回拍出量を2倍にできると言う
ことができる。実際に、ジー・ジェー・ランゲウーテル
ズ(G.J.Langewouters)他のジャーナル・オブ・バイオ
メカニックス(J.Biomechanics)17、ページ425−435
(1984)の論文より、何年も前から、圧力pと拍出量V
との間の関係はアークタンジェント関数に従った強度の
非線形である(図3a)ことが知られている。更に、前記
アークタンジェント関数は、患者の年齢及び性別に依存
する。最後に、大動脈から末梢血管部への流出は1カ所
に集中せず大動脈に沿って分布する。言い換えれば、大
動脈の始点における増加圧力は、これが末梢血管部への
種々の分岐点に達する以前にいくらかの時間を要し、増
加した大動脈圧力のため実際に流出が増加する。
要約すれば、これらの無視は、大動脈が線形の挙動を
しないこと、及び大動脈に沿った波の移動による時間の
要約があるということである。
しかし、実際には、別のモデル、いわゆる伝達経路モ
デル(図1a)が想定される。この近似においては、大動
脈は液体(血液)で満たされ、かつ終端部に集中した末
梢血管部への流出のある「適切な」長さの一様な弾性管
として考えられる。かかる管は二つの特性を持つ。管の
始点において心臓により作られた圧力波は有限速度vp
管の終点に移動し、そこで一部が反射されて管の始点に
戻る。反射波がここに到着するまでに時間がかかるの
で、心臓の送出期、収縮期は既に終わり、大動脈弁は閉
じられ、結局、反射波は閉じられた弁により抑止される
ので、心臓は反射波から何の障害も受けない。「適切な
長さ」とは、波の伝播速度を与える管が、反射圧力波の
直ぐには戻らないことを保証する長さであることであ
る。伝播速度vpは、大動脈の断面積A、大動脈の単位長
さ当たりのコンプライアンスC'及び血液の比重ρによ
り、最初に、次式のように決定される: 伝播速度に加えて、第2の特性は特性インピーダンス
Zoである。このインピーダンスは、圧力波の振幅と大動
脈における随伴流の波の振幅との間の比率を反映し、か
つ断面積A及びコンプライアンスC'にも影響される。実
際上、特性インピーダンスは次式で与えられる: 圧力p(t)との流量q(t)との間の比率が一定であ
りかつ式q(t)=p(t)/ZOで与えられることは、
少なくとも反射圧力(及び流量)の波が戻るときまで、
大動脈の始点における圧力液と流量波とは形状が同じで
あることを意味する。
このモデルは、1回拍出量の計算にはあまり重要では
ないが、ある状況下では計算をあまりにも不正確とする
種々の近似を含む。第1に、このモデルでは、実際に
は、大動脈は圧力pと拍出量V、従って一定長さの大動
脈の場合の断面積A=V/lとの間の関係に関して線形的
な挙動をしない。単位長さ当りコンプライアンスC'は、
通常の圧力poにおいては断面積の圧力に関する導関数、
即ちC'=[dA/dp]であるので、Aは圧力の上昇とと
もに増加するであろうが、同時に、大動脈が前述のアー
クタンジェント関数(図3b)に従って圧力上昇により膨
張しにくくなるのでC'は減少するであろう。特定の圧力
領域においては、これら2個の要因はZOについての式に
おいて互いに補償し、後者はZOについては見かけ上一定
となるであろう。より高い圧力においてはC'の減少がA
の増加よりも速やかに生じ、特性インピーダンスZOはA
とC'との積の逆数の平方根に従って徐々に増加し始める
であろう。従って、この非線形は役に立つが、しかしウ
インドケセル・モデルほどではない。第2の近似は、大
動脈を一様とすることである。実際は、大動脈は断面が
一様でないだけでなく、血管壁の膨張性も一様でない。
大動脈断面積は末梢に向かって減少する。最初は、大動
脈の分岐動脈は、圧力波及び流量波がその分岐に達する
と直ちに全断面積に寄与し始めるので前記減少は補償さ
れる。しかし、更に血管壁は末梢に向かって膨張性が減
少する。これは、それ自体がコンプライアンスC'の減少
を明示する。特性インピーダンスの式中のA及びC'の両
者は、大動脈を通じて末梢に向かって減少するため、zo
は増加するであろう。このことは、モデルに仮定された
大動脈管の終点における末梢部の抵抗の単一反射の問題
はないが、大動脈全体に沿い始めは小さいが収縮期に後
者を相当に増加させる反射が分布することを意味する。
収縮期は短いので、心臓は、収縮時に大動脈の始点にお
いて(低い)インピーダンスに合うだけであるが、長く
続いた収縮では心臓は(高い)末梢部インピーダンスも
経験する。これは、血液の瞬間排出量が増加する反射イ
ンピーダンスにより徐々に抑止されるので、これが収縮
の進行と共に次第に小さくなることを意味する。
第3の近似は次の通りである。血液が圧力の増加によ
り大動脈内に圧入された後に、次の拡張期の圧力が直ち
にゼロにならないことは勿論であるが、大動脈は末梢抵
抗を経て徐々に空虚になるであろう。これに続いて、次
の拍動が幾らか増加した圧力に対して押し出され、事実
上、この二つの圧力が互いに重ねられるであろう。拡張
圧力は、収縮時の内向きの流れと拡張時の外向きの流れ
とが平衡に達するまで、各拍動毎に常に増加するであろ
う。これがウインドケセルの機能を補償する。
最後に、流出が大動脈の終点に集中せずその全長に亙
って分布されることもまた真実である。
纏めれば、これらの無視は、大動脈が非線形の性質を
有しかつ一様として考えられないことの結果である。
上のモデルの概念は、両者とも、その目的がモデルを
使用して収縮中の時間の関数として大動脈圧力から大動
脈流量のパターンを算出することであるならば適切でな
い。伝送線路モデルにより見いだされる曲線は持続時間
の短い収縮という条件下の適宜量の実際の流量曲線に似
ているだけである。ウインドケセルモデルは、実際、流
量曲線の計算もできない。実際は、両方のモデルでは、
心臓の1回拍出量のみ、即ち1回の収縮について積分さ
れた流量曲線が圧力曲線から作られるだけである。
本発明の目的は、上述の問題を無くし、かつウインド
ケセル特性及び伝送線路特性の両者を考慮した1回拍出
量及び心臓血液拍出量の前記決定のための装置を提供す
ることである。このことは、圧力p(t)から流量q
(t)を算出しかつ収縮期に亙る流量を積分することに
より、最初に説明された形式の方法で達成される。この
場合、大動脈は、ウインドケセルコンプライアンス、及
びアークタンジェント関係としての大動脈における圧力
/流量に関係で補足された伝送線路として考えられる。
この場合は、非線形の大動脈特性が重要でありこの方法
に含まれるが、もし流量曲線が、収縮期中、連続的に計
算され表示できるならば、この方法の諸段階の精度の検
証のために有利である。この場合、流量曲線の積分より
1回拍出量が得られる。
本発明は、図面を参照した実施例に基づきより詳細に
説明されるであろう。図面において、 図1a及び1bは、公知の2種の脈拍計数器モデル、即ち
ウインドケセルモデル及び均一伝送線路モデルを示す。
図2a及び2bは、それぞれ、ウインドケセル機能により
補足された伝送線路モデルに基づく本発明による脈拍計
数器モデルの末梢部流出のある場合とない場合とを示
す。
図3aは、40歳男子の大動脈に対する様式化された圧力
/断面積及び単位長さ当たりの圧力/コンプライアンス
の関係のグラフを示す。
図3bは、初期値として動脈の最終拡張期圧力を使用し
た線形化圧力の計算の線図的表現を示す。
図4a及び4bは、それぞれ、脈拍計数器モデルの支援に
より計算された入力圧力信号p(t)及び計算流量曲線
q(t)を示す。
図5は、熱希釈決定数に対して描かれた6名の心臓切
開手術患者の脈拍計数器(線)及び熱希釈(o)心臓血
液拍出量のグラフを多数示す。
本発明は、伝送線路モデル、特に圧力波及び流量波が
固定要因、少なくとも末梢部から心臓に反射が戻るより
以前の特性インピーダンスZoと関係する伝送線路特性を
仮定する。このモデルは早期の反射による心臓からの外
向き流の漸進的な抑止及び拡張期ウインドケセル作用を
欠如している。本発明により、これら二つの欠点は、伝
送線路モデルに1個の要素手段、大動脈ウインドケセル
コンプライアンスCwを付加することにより無くされる
(図2a参照)。結局、収縮期においては、流量は、最初
はZoにより決定され、次いで、流入量q(t)によりウ
インドケセルレシーバー内に作り上げられたpw(t)の
圧力上昇により送出期中に段々と抑止される。
別な省略は、末梢部流出が伝送線路の終点に集中する
と考えることである。本発明による相当末梢部抵抗RP
ウインドケセルコンプライアンスCWとの新しい組合せの
結果、後者は、(低域濾波されたのもではあるが)平均
圧力と連結され、従ってより良く現実を表す。図2b参
照。このモデルについては、内向き流量q(t)から外
向き流量 が差し引かれ、従って内向き流の抑止は幾分かゆっくり
と生ずる。
本発明による改良されたモデルは線形の時間不変要素
を含む。実際は、圧力/断面積の関係がアークタンジェ
ントに従って変化する強度の非線形大動脈に関心があ
る、本発明により、2種のバージョンが可能である。第
1のバージョンにおいては、この方法で記録された圧力
曲線がまず線形化され、次いで線形脈拍輪郭モデルに適
用される。第2のバージョンにおいては、この方法で測
定された圧力は所与と考えられるが、脈拍輪郭モデルに
は圧力依存要素が設けられる。
線形化された圧力曲線モデル: この第1のバージョンにおいては、ある量の血液が心
臓により大動脈内に圧送され、かつこれに続いて圧力が
上昇するので、圧力曲線は上昇すると仮定する。圧力の
上昇の大きさは大動脈の圧力/拍出量の関係に依存す
る。高圧においては、少しの拍出量増加が大きな圧力上
昇をもたらし、一方、低圧においては同じ拍出量はより
小さな圧力上昇をもたらすであろう(アークタンジェン
ト関係)。しかし、既知であると仮定された圧力/拍出
量の関係の手段により、大動脈が線形の圧力/拍出量の
関係を持ったときに存在したであろう圧力曲線を算出す
ることが可能である。この方法は簡単な方法でグラフ的
にこれを示すことができる(図3b参照)。各拍動に対し
て、圧力曲線の放出期の始点及び組み合わせられた最終
の拡張期圧力pdが決定される。前記圧力において、非線
形の大動脈圧力/拍出量の関係、又は大動脈の長さを一
定としたときの圧力/断面積関係は、曲線のタンジェン
トの手段により線形化され、圧力曲線上で得られる圧力
サンプル毎に相当線形圧力が計算される。初期圧力を70
mmHgと仮定すると、圧力サンプル100mmHgのこの例で
は、相当線形圧力はほぼ92mmHgとなる。次いで、この圧
力が脈拍輪郭モデルに加えられる。このモデルは、適用
可能な圧力/断面積の関係から算出された特性インピー
ダンス及びウインドケセルコンプライアンスについて、
関係の拍動に対する最終拡張期圧力と組み合わせられた
値を持つ。
非線形モデル: この第2のバージョンにおいては、圧力曲線は、上と
同じく、そのウインドケセル特性だけでなくその総ての
特性を与えられた非線形大動脈内への血液の外向き流の
結果として上昇する。このとき、圧力曲線には線形修正
は適用されず、そのまま受け入れられる。しかし、関係
の拍動の拡張期中は、特性インピーダンス及びウインド
ケセルコンプライアンスの値は、いかなるときもこれが
存在するように圧力値に従って変化し、かつ上述の式に
従って非線形アークタンジェントの圧力/断面積の関係
から計算される。圧力サンプルが得られたときはいつで
も、組み合わせられた特性インピーダンス及びウインド
ケセルコンプライアンス、並びに流れにおけるこれらの
変化された値の効果が算出される。
アークタンジュント関係の使用: ランゲウーテルズ他によれば、大動脈断面積と大動脈
圧力との間の関係は、次式を用いてこれを特定すること
がてきる: A(p)=Am[0.5+(1/π.)tan-1{(p−po)/p1}] (1) ここに、A(p)は圧力pの関数としての断面積、Am
非常に高圧力における限界断面積、そしてp0とp1とはそ
れぞれ関係式の変曲点における圧力と変曲点における関
係式の傾斜である。典型的な曲線は図3a参照。
圧力の関数としての単位長さ当たりのコンプライアン
スC'はpに関しての微分の結果として以下の通りであ
る: ここに、Cm=Amp1はp=p0における最大コンプライ
アンスである。パラメーターAm、p0及びp1の値はランゲ
ウーテルズ他の研究から知られる。特に、これらのパラ
メーターは人の性別及び年齢に依存する。更に、Amの値
には約20%の分散がある。従って、約20%の標準偏差を
許容すると考えない限り、より正確な較正なしに絶対的
な意味における脈拍輪郭モデルにおける要素の値を決定
することはできない。
Z0及びvpは、この方法で見いだされたA及びC'の支援
で容易に計算できる。大動脈のウインドケセルコンプラ
イアンスを計算するためには、更に1個のパラメータ
ー、即ち大動脈の有効長さを必要とする。ウインドケセ
ルコンプライアンスは、単位長さ当たりのコンプライア
ンスと有効長さとの積、CW=C'である。成人についての
良好な値は60cmである。
末梢部抵抗の値: 図2bの脈拍輪郭モデルは、並列な総ての血管排出床の
全末梢部抵抗Rpを含む。後者の値は最初は未知数であ
る。最初に1回拍出量Vsを算出した後で、心拍数f,即
ち、1分間当たりの心臓の拍動数は、f=60/Tとして、
収縮時間と拡張期間との和T=Ts+Tdから見いだせる。
心臓血液拍出量は心拍数と関係の拍動に対する1回拍出
量との積として計算することができる; Q=fVs (3) 平均大動脈圧力も、ある時間にわたる圧力の積分をそ
の時間で割ったものとして知られる: 次いで、末梢部抵抗の値は容易に次のようになる: RP=pm/Q (5) 末梢部抵抗は緩慢にしか変化しないので、今の拍動から
算出されたその値は、これを次の拍動の流量計算に使用
することができる。第1の拍動に対しては、合理的な初
期値の推定を行わねばならないであろう。収縮期におけ
る流量値は末梢部抵抗値には強く依存しないので、この
方法で急速な収束が得られる。
脈拍輪郭モデルを使用した流量の計算: 適宜の時間における流量qは、シミュレーションを使
用してこれを計算することができる。これに関連して、
2bのモデルの挙動を記述する式が作られる。大動脈の流
量は次式に従う: q(t)=1/Z0)・{(p(t)−pw(t)} 末梢部抵抗を通る流量は次式に従う:qR(t)=pw/
RP;また、エアレシーバーコンプライアンスCWは次式に
従う: VW(t)=V0+∫{(q(t)−qR(t)}dt,ここ
に、V0はCWにおける最初の拍出量である。最後に、エア
レシーバー内の圧力は拍出量/圧力の関係pw(t)=f
(VW(t))に従う。
先に述べたように、1回拍出量VSは収縮期中の流量を
積分することにより流量曲線から見いだすことができ、
更に式(3)、(4)及び(5)に示されるように、心
拍数fは60を各拍動時間で割った商として、また心臓血
液拍出量Qは1回拍出量VSと心拍数fとの積として、そ
の末梢部抵抗RPは平均圧力pmを心臓血液拍出量Qで割っ
た商として見いだすことができる。
62歳の男子患者についての測定された圧力曲線及びこ
れから算出された流量曲線が図4に示される。各拍動の
前に、図4の圧力曲線の上昇端の始点がパターン認識ア
ルゴリズムを用いて決定される。このとき、曲線の拡張
期圧力pdが記録され、更に後者により、特性インピーダ
ンス及びウインドケセルコンプライアンスの正確な値が
関数発生器を使用して見いだされ、モデルに適用され
る。圧力は、アークタンジェント圧力/断面積の関係を
使用して線形化され、モデルに加えられる。アークタン
ジェント圧力/断面積の関係を決定するパラメーターの
正確な値は、ランゲウーテルズ他による回帰式に従い、
年齢n歳の女子及び男子についてそれぞれ次の通りであ
る。
女子 男子 Am=4.12 Am=5.62 p0=72−0.89n p0=76−0.89n p1=57−0.44n p1=57−0.44n ここに、Amはcm2で、圧力pはmmHgで、そして年齢n
は年で表す。
明らかに、これらの計算はデジタルプロセッサーによ
り有利に遂行することができる。この目的に対しては、
連続圧力曲線を規則的間隔を有する一連のサンプルに置
換し、得られた各サンプル毎に計算を繰り返すことが必
要である。従って、流量曲線はこの時間間隔の各瞬間に
利用可能となる。好ましいサンプル間隔10msである。高
速コンピューターが利用できれば精度の限界の改良だけ
でなく、より短い5msの時間間隔が可能である。1回拍
出量の計算のための流量曲線の積分は、収縮期における
流量のサンプルの足し算で置換される。
新しい方法の精度: この新しい方法は心臓切開手術を受けた6名の患者で
評価された。心臓血液拍出量についての基準値は、呼吸
サイクルに亙って正確に規則的に分布された冷液の注入
後に少なくとも4回、熱希釈(TD)定量を繰り返すこと
により手術中の何回かの瞬間に測定された。1回の熱希
釈定量はあまりにも不正確であった。特別な特徴は、呼
吸サイクルに亙り規則的に広がるように4回の定量が注
入されることであり、その結果、この4個の値を平均す
ることは、実際の心臓血液拍出量の推定量の精度を係数
4で改良することを意味する。この改良は機械的吸入法
により達成しただけであり、従ってもし患者が自発呼吸
するならば達成しない。図5に記されたように、自発呼
吸は第1の定量中に、即ちsにおいて常に生じ、続く総
ての定量中、患者は麻酔をかけられ人工通気される。
この新奇な脈拍輪郭法により算出された心臓血液拍出
量は比較のために同時に記録される。その結果が図5に
各患者後とに別々に示される。患者番号8910から9212及
び8915から8917についての同時定量の数、平均誤差、誤
差の標準偏差(SD)及び誤差範囲が表1に示され、かつ
グループで平均された。6番目の患者は3個の定量値に
つき大きな誤差がある点に注意されたい。これらは、全
心臓−長バイパスの期間後に行われ、この間、心臓は既
に鼓動し血液を圧送していたが心肺機械のポンプはなお
毎分約1リットルで支援していた。この流量は、実際は
熱希釈法で測定されたのもであり脈拍輪郭法によったも
のではなく、この矛盾を生じた。
比較のために、心臓血液拍出量はウェセリング他のUS
3841313による初期の方法によっても計算された。この
方法は文献において最も正確であるとして知られる。こ
の場合の平均誤差の標準偏差は、新奇な方法の8%に対
して12%であった。
概説すれば、この新方法での最後の2列に示された平
均誤差は、旧方法によるよりも計数 以上小さかった。
手術室又は集中治療室における患者の集中看視中は、
大動脈の動脈圧力を記録することは稀に可能であるだけ
である。患者に対する便利と安全の理由で、通常は橈骨
又は上腕又は指の動脈の圧力が非侵略的に記録されるで
あろう。この圧力は変形され大動脈圧力と比較される。
即ち圧力曲線の形状が変えられる。これは次の事実によ
る。即ち、大動脈において作られた圧力曲線が動脈系統
を経て、より末梢に置かれた上腕、脈拍又は指に伝わる
こと、及びその伝達経路はある種の特性を有することで
あり、この伝達経路の特性中で最も重要なものはその狭
帯特性とその低周波共振である。多数の患者のこれら特
性を測定することは可能である。これに関連し、発明者
は、伝達特性は平均値を中心とした比較的小さな範囲に
あることを見いだした。この事例では、濾波器内の反共
振により共振特性を修正することが可能であり、その結
果、心臓の左半分の脈拍輪郭1回拍出量を正確に決定す
ることが実際に可能である。良好な近似のためには、こ
の目的で1個の2次濾波器を使用できる。実際上は、デ
ジタルのII R(不定期間インパルス反応)2次又は3次
の濾波器を使用し、平均共振を周波数8Hzまで正確に修
正する。
この場合、一連の入力信号サンプルxnに基づいて一連
の出力信号サンプルynが次式を使用して計算される: ここに、c及びdは一定の濾波器定数である。2次濾波
器に対しては、次式が使用される。
C0=0.077 d1=1.484 d2=−0.561
フロントページの続き (56)参考文献 特開 昭58−216034(JP,A) 特開 平2−1230(JP,A) W.H.GUIER et al., Beat−by−Beat Strok e Volume from Aort ic−Pulse−Pressure Analysis,IEEE Tran sactions on BIOMED ICAL ENGINEERING,V ol.21/No.4,285−292 G.J.LANGEWOUTERS et al.,THE STATIC ELASTIC PROPERTIES OF 45 HUMAN THORAC IC AND 20 ABDOMINAL AORTAS IN VITRO A ND THE PARAMETE,Jo unal of Biomechani cs,Vol.17/No.6,425−435 (58)調査した分野(Int.Cl.7,DB名) A61B 5/00 - 5/0295

Claims (9)

    (57)【特許請求の範囲】
  1. 【請求項1】対象の動脈血圧を測定して、収縮及び拡張
    期からなる対応する大動脈血流圧信号を生成するセンサ
    ーと、 大動脈における非線型圧力−量関係のアークタンジエン
    ト関係を提供する、ウインドケセルコンプライヤンスで
    補足された伝送路としての、大動脈のシミュレーション
    モデルと、 上記モデル内に圧力を導入し、時間の関数として、そこ
    から出る流れを計算する導入手段と、 収縮期にわたって流れを積分し、拍出量を提供する積分
    手段と を具備することを特徴とする対象の脈拍状の大動脈血流
    圧信号から血液拍出量を決定するための装置。
  2. 【請求項2】アークタンジェント圧力−容積関係が、 A(p)=Am[0.5+(1/π)・tan-1{(p−pO)/
    p1}] ここで、 A(p)は圧力pの関数としての断面積、 Amはごく高圧における限界断面積、 p0はアークタンジェント関係の変曲点における圧力、 p1は変曲点におけるアークタンジェント関係の傾斜 に従う、大動脈圧力/断面積から、固定大動脈長さで、
    導かれる請求項1の装置。
  3. 【請求項3】測定された圧力がシミュレーションモデル
    に連続的に導入され、 各拍動毎にアークタンジエントの関係に従って、伝送線
    路の特性インピーダンス及びウインドケセルコンプライ
    アンスが、ウインドケセルコンプライアンスの圧力に、
    連続的に適用される請求項2の装置。
  4. 【請求項4】大動脈の線形圧力/容積関係の基礎で圧力
    が線形され、 最終拡張圧力と関連する伝送線路の特性インピーダンス
    の値及びウインドケセルコンプライアンスを維持して、
    上記圧力がパルス輪郭モデルに導入される請求項2の装
    置。
  5. 【請求項5】流量q(t)が q(t)=(1/Z0)・{(p(t)−pw(t)} ここで、 pwはコンプライアンス圧力、 ρは血液の比重、 C'(p)=Cm/{1+(p−p02/p1 2} ここで、 Cm=Am/πp1はp=p0における最大コンプライアンス の式により計算される請求項4の装置。
  6. 【請求項6】心臓血液拍出量Qが1回拍出量と心拍数と
    の積として計算される請求項1の装置。
  7. 【請求項7】Rp=pm/Q の式による末梢部抵抗RPが計算され、ここで、 pm=平均圧力=(1/T)・∫p(t)dt であり、 qRp(t)=pw/Rp の式で、末梢部抵抗RPから流量が計算されるる請求項6
    の装置。
  8. 【請求項8】橈骨又は上腕の動脈血圧を測定する手段を
    備えており、伝達により生じた圧力曲線の共振歪を修正
    するために圧力信号に濾波作業が行われる請求項1の装
    置。
  9. 【請求項9】指で、非侵略的に血圧を測定する手段を備
    えており、 伝達により生じた圧力曲線の共振歪を修正するために圧
    力信号に濾波作業が行われる請求項1の装置。
JP50509092A 1991-01-29 1992-01-28 血液拍出量を決定するための装置 Expired - Lifetime JP3293621B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
NL9100150A NL9100150A (nl) 1991-01-29 1991-01-29 Werkwijze voor het bepalen van het slagvolume en het hartminuutvolume van het menselijk hart.
NL9100150 1991-01-29
PCT/NL1992/000022 WO1992012669A1 (en) 1991-01-29 1992-01-28 Method of determining the stroke volume and the cardiac output of the human heart

Publications (2)

Publication Number Publication Date
JPH06505651A JPH06505651A (ja) 1994-06-30
JP3293621B2 true JP3293621B2 (ja) 2002-06-17

Family

ID=19858791

Family Applications (1)

Application Number Title Priority Date Filing Date
JP50509092A Expired - Lifetime JP3293621B2 (ja) 1991-01-29 1992-01-28 血液拍出量を決定するための装置

Country Status (10)

Country Link
US (1) US5400793A (ja)
EP (1) EP0569506B1 (ja)
JP (1) JP3293621B2 (ja)
AT (1) ATE128837T1 (ja)
CA (1) CA2101531C (ja)
DE (1) DE69205417T2 (ja)
DK (1) DK0569506T3 (ja)
ES (1) ES2080489T3 (ja)
NL (1) NL9100150A (ja)
WO (1) WO1992012669A1 (ja)

Families Citing this family (70)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5584298A (en) * 1993-10-25 1996-12-17 Kabal; John Noninvasive hemodynamic analyzer alterable to a continuous invasive hemodynamic monitor
US5836884A (en) 1993-12-17 1998-11-17 Pulse Metric, Inc. Method for diagnosing, monitoring and treating hypertension and other cardiac problems
US5535753A (en) * 1994-10-04 1996-07-16 Rutgers University Apparatus and methods for the noninvasive measurement of cardiovascular system parameters
WO1997016114A1 (fr) * 1995-11-01 1997-05-09 Seiko Epson Corporation Appareil pour mesurer l'etat d'un organisme vivant
US6569103B2 (en) 1995-12-22 2003-05-27 Arrow International Investment Corp. Device for determining a characteristic point in the cardiac cycle
US6010457A (en) * 1996-03-15 2000-01-04 Pmv Medical Pty Ltd Non-invasive determination of aortic flow velocity waveforms
US5971934A (en) * 1996-10-04 1999-10-26 Trustees Of The University Of Pennsylvania Noninvasive method and apparatus for determining cardiac output
GB9714550D0 (en) 1997-07-10 1997-09-17 Lidco Ltd Improved method and apparatus for the measurement of cardiac output
DE19814371A1 (de) * 1998-03-31 1999-10-14 Pulsion Verwaltungs Gmbh & Co Verfahren zur in-vivo Bestimmung der Compliance-Funktion und des systemischen Blutflusses eines Lebewesens und Vorrichtung zur Durchführung der Verfahren
US6117087A (en) * 1998-04-01 2000-09-12 Massachusetts Institute Of Technology Method and apparatus for noninvasive assessment of a subject's cardiovascular system
US6986744B1 (en) 1999-02-02 2006-01-17 Transonic Systems, Inc. Method and apparatus for determining blood flow during a vascular corrective procedure
US6315735B1 (en) 1999-03-31 2001-11-13 Pulsion Medical Systems Ag Devices for in-vivo determination of the compliance function and the systemic blood flow of a living being
IT1315206B1 (it) * 1999-04-27 2003-02-03 Salvatore Romano Metodo e apparato per la misura della portata cardiaca.
US6939303B2 (en) * 1999-07-23 2005-09-06 Provecto Tecnico R-2000, Ca Method and system of determining cardiac contractility
US6440078B1 (en) 1999-07-23 2002-08-27 Roberto Curiel Method and system of determining cardiac contractility
AU776098B2 (en) * 1999-11-24 2004-08-26 Edwards Lifesciences Corporation Method and apparatus for determining cardiac output or total peripheral resistance
AUPQ420599A0 (en) * 1999-11-24 1999-12-16 Duncan Campbell Patents Pty Ltd Method and apparatus for determining cardiac output or total peripheral resistance
US6510337B1 (en) * 1999-11-26 2003-01-21 Koninklijke Philips Electronics, N.V. Multi-phase cardiac imager
NL1016320C2 (nl) 2000-10-03 2002-04-04 Jozef Reinier Cornelis Jansen Inrichting voor het aansturen van hartondersteunende apparaten.
US7011631B2 (en) * 2003-01-21 2006-03-14 Hemonix, Inc. Noninvasive method of measuring blood density and hematocrit
US7815578B2 (en) * 2003-02-10 2010-10-19 Massachusetts Institute Of Technology Methods and apparatus for determining cardiac output
FR2860980B1 (fr) * 2003-10-16 2005-12-30 Ela Medical Sa Dispositif medical implantable actif comprenant des moyens d'ajustement de la frequence maximale de stimulation ventriculaire en fonction de l'etat hemodynamique du patient
US7452333B2 (en) 2003-12-05 2008-11-18 Edwards Lifesciences Corporation Arterial pressure-based, automatic determination of a cardiovascular parameter
US7220230B2 (en) 2003-12-05 2007-05-22 Edwards Lifesciences Corporation Pressure-based system and method for determining cardiac stroke volume
US7422562B2 (en) 2003-12-05 2008-09-09 Edwards Lifesciences Real-time measurement of ventricular stroke volume variations by continuous arterial pulse contour analysis
US7813808B1 (en) 2004-11-24 2010-10-12 Remon Medical Technologies Ltd Implanted sensor system with optimized operational and sensing parameters
US20060167361A1 (en) * 2005-01-27 2006-07-27 Bennett Tommy D Method and apparatus for continuous pulse contour cardiac output
US7935062B2 (en) * 2005-01-27 2011-05-03 Medtronic, Inc. Derivation of flow contour from pressure waveform
US7651466B2 (en) * 2005-04-13 2010-01-26 Edwards Lifesciences Corporation Pulse contour method and apparatus for continuous assessment of a cardiovascular parameter
DE602006007989D1 (de) * 2005-05-19 2009-09-03 Jansen Jozef Reinier C T dem patienten in beziehung stehenden parameter zur patientenüberwachung
US7742815B2 (en) * 2005-09-09 2010-06-22 Cardiac Pacemakers, Inc. Using implanted sensors for feedback control of implanted medical devices
AT502922B1 (de) * 2005-11-08 2011-07-15 Arc Seibersdorf Res Gmbh Verfahren und vorrichtung zur bestimmung des herzzeitvolumens
US20070142727A1 (en) * 2005-12-15 2007-06-21 Cardiac Pacemakers, Inc. System and method for analyzing cardiovascular pressure measurements made within a human body
US8343061B2 (en) 2006-03-15 2013-01-01 Board Of Trustees Of Michigan State University Method and apparatus for determining central aortic pressure waveform
US8905939B2 (en) 2006-07-13 2014-12-09 Edwards Lifesciences Corporation Method and apparatus for continuous assessment of a cardiovascular parameter using the arterial pulse pressure propagation time and waveform
US7955268B2 (en) * 2006-07-21 2011-06-07 Cardiac Pacemakers, Inc. Multiple sensor deployment
WO2008144404A1 (en) * 2007-05-16 2008-11-27 Massachusetts Instutute Of Technology Systems and methods for model-based estimation of cardiac output and total peripheral resistance
WO2008144525A1 (en) 2007-05-16 2008-11-27 Massachusetts Institute Of Technology System and method for prediction and detection of circulatory shock
US8235910B2 (en) * 2007-05-16 2012-08-07 Parlikar Tushar A Systems and methods for model-based estimation of cardiac ejection fraction, cardiac contractility, and ventricular end-diastolic volume
US8465435B2 (en) * 2007-07-20 2013-06-18 Bmeye B.V. Method, a system and a computer program product for determining a beat-to-beat stroke volume and/or a cardiac output
DE102007057235A1 (de) * 2007-11-28 2009-06-04 Iprm Intellectual Property Rights Management Ag Verfahren und Vorrichtung zur Bestimmung des Endes des systolischen Teils eines Druckverlaufs
US20090270739A1 (en) * 2008-01-30 2009-10-29 Edwards Lifesciences Corporation Real-time detection of vascular conditions of a subject using arterial pressure waveform analysis
WO2009102613A2 (en) * 2008-02-11 2009-08-20 Cardiac Pacemakers, Inc. Methods of monitoring hemodynamic status for ryhthm discrimination within the heart
WO2009102640A1 (en) 2008-02-12 2009-08-20 Cardiac Pacemakers, Inc. Systems and methods for controlling wireless signal transfers between ultrasound-enabled medical devices
DE102008008840A1 (de) * 2008-02-13 2009-09-24 Up Management Gmbh Verfahren und Vorrichtung zur nicht-invasiven Messung von dynamischen Herz-Lungen Interaktionsparametern
CA2713675C (en) * 2008-02-15 2017-08-22 Assistance Publique-Hopitaux De Paris Device and process for calculating new indices of arterial stiffness, and/or for stroke volume monitoring
KR101006198B1 (ko) * 2008-04-16 2011-01-12 ㈜비엔컴 블루투스를 이용한 핸즈프리 시스템
US10251566B2 (en) 2008-05-27 2019-04-09 Board Of Trustees Of Michigan State University Methods and apparatus for determining a central aortic pressure waveform from a peripheral artery pressure waveform
US20100152547A1 (en) 2008-07-02 2010-06-17 Sterling Bernhard B Method and system for determining cardiac performance
WO2010027652A1 (en) * 2008-08-26 2010-03-11 Cardiac Pacemakers, Inc. Cardiac output estimation using pulmonary artery pressure
EP2334230A1 (en) * 2008-10-10 2011-06-22 Cardiac Pacemakers, Inc. Systems and methods for determining cardiac output using pulmonary artery pressure measurements
GB0821084D0 (en) * 2008-11-18 2008-12-24 King S College London Apparatus and method
US8632470B2 (en) 2008-11-19 2014-01-21 Cardiac Pacemakers, Inc. Assessment of pulmonary vascular resistance via pulmonary artery pressure
US20100324378A1 (en) * 2009-06-17 2010-12-23 Tran Binh C Physiologic signal monitoring using ultrasound signals from implanted devices
US20110040194A1 (en) * 2009-08-12 2011-02-17 Rajendra Prasad Jadiyappa Method and system for determining cardiac parameters
CN102791192B (zh) * 2009-10-29 2015-02-25 Cn体系药物技术有限公司 增强和分析来自连续无创血压设备的信号的装置
US20130053664A1 (en) 2010-01-29 2013-02-28 Edwards Lifesciences Corporation Elimination of the effects of irregular cardiac cycles in the determination of cardiovascular parameters
US9782089B2 (en) 2010-09-08 2017-10-10 Siemens Healthcare Gmbh Worksheet system for determining measured patient values for use in clinical assessment and calculations
WO2013052896A1 (en) 2011-10-07 2013-04-11 Edwards Lifesciences Corporation Detecting a vasoactive agent in the bloodstream
CN102499658B (zh) * 2011-11-08 2013-12-11 中国科学院深圳先进技术研究院 中心血压波形重建模型及重建方法
EP3019075B1 (en) 2013-07-08 2020-03-25 Edwards Lifesciences Corporation Determination of a hemodynamic parameter
WO2015130705A1 (en) 2014-02-25 2015-09-03 Icu Medical, Inc. Patient monitoring system with gatekeeper signal
WO2017070120A1 (en) 2015-10-19 2017-04-27 Icu Medical, Inc. Hemodynamic monitoring system with detachable display unit
CN105496360B (zh) * 2015-12-03 2018-06-19 北京大学人民医院 一种辅助判断全人群系统血管顺应性的试剂盒
CN105361886B (zh) * 2015-12-03 2018-06-19 北京大学人民医院 一种辅助判断亚健康人群系统血管顺应性的试剂盒
DE102015015784A1 (de) * 2015-12-08 2017-06-08 WSH Engineering Services GbR (vertretungsberechtigter Gesellschafter: Hieronymi, Andreas, 61440 Oberursel) Verfahren zur Erfassung von arteriellen Parametern eines Menschen und Vorrichtung zur Durchführung des Verfahrens
DE102016001710B4 (de) 2016-02-15 2022-08-25 Fresenius Medical Care Deutschland Gmbh Gerät zur extrakorporalen Blutbehandlung mit einer Auswerte- und Steuereinheit
EP3427651A4 (en) * 2016-04-15 2019-12-11 Omron Corporation DEVICE AND SYSTEM FOR ANALYZING BIOLOGICAL INFORMATION AND PROGRAM THEREOF
CN111493855B (zh) * 2020-04-21 2023-01-06 重庆理工大学 个体化心输出量的无创测量系统与方法
GB2622396A (en) * 2022-09-14 2024-03-20 Carelight Ltd Real time opto-physiological monitoring method and system

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4632125A (en) * 1984-01-13 1986-12-30 American Hospital Supply Corp. Right heart ejection fraction and cardiac output catheter
JPS62207435A (ja) * 1986-03-07 1987-09-11 テルモ株式会社 心拍出量測定用カテ−テル
US5211177A (en) * 1990-12-28 1993-05-18 Regents Of The University Of Minnesota Vascular impedance measurement instrument
US5183051A (en) * 1991-01-14 1993-02-02 Jonathan Kraidin Means and apparatus for continuously determining cardiac output in a subject

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
G.J.LANGEWOUTERS et al.,THE STATIC ELASTIC PROPERTIES OF 45 HUMAN THORACIC AND 20 ABDOMINAL AORTAS IN VITRO AND THE PARAMETE,Jounal of Biomechanics,Vol.17/No.6,425−435
W.H.GUIER et al.,Beat−by−Beat Stroke Volume from Aortic−Pulse−Pressure Analysis,IEEE Transactions on BIOMEDICAL ENGINEERING,Vol.21/No.4,285−292

Also Published As

Publication number Publication date
DK0569506T3 (da) 1995-12-04
CA2101531C (en) 2003-01-14
EP0569506A1 (en) 1993-11-18
NL9100150A (nl) 1992-08-17
EP0569506B1 (en) 1995-10-11
DE69205417D1 (de) 1995-11-16
US5400793A (en) 1995-03-28
WO1992012669A1 (en) 1992-08-06
DE69205417T2 (de) 1996-05-30
ATE128837T1 (de) 1995-10-15
ES2080489T3 (es) 1996-02-01
JPH06505651A (ja) 1994-06-30
CA2101531A1 (en) 1992-07-30

Similar Documents

Publication Publication Date Title
JP3293621B2 (ja) 血液拍出量を決定するための装置
Raines et al. A computer simulation of arterial dynamics in the human leg
Forster et al. Oscillometric determination of diastolic, mean and systolic blood pressure—a numerical model
US8282564B2 (en) Systems and methods for model-based estimation of cardiac output and total peripheral resistance
US8905939B2 (en) Method and apparatus for continuous assessment of a cardiovascular parameter using the arterial pulse pressure propagation time and waveform
US20080015451A1 (en) Method and Apparatus for Continuous Assessment of a Cardiovascular Parameter Using the Arterial Pulse Pressure Propagation Time and Waveform
WO2009139646A1 (en) Method for estimating a central pressure waveform obtained with a blood pressure cuff
JP2002253519A (ja) 血液量測定方法及び生体信号モニタ装置
WO2008019207A2 (en) Method and apparatus for continuous assessment of a cardiovascular parameter using the arterial pulse pressure propagation time and waveform
EP3185767A1 (en) Method for oscillatory non-invasive blood pressure (nibp) measurement and control unit for an nibp apparatus
US11678808B2 (en) Haemodynamic data estimation apparatus and method of use
JP2000501327A (ja) 超音波エコーグラフィーによって動脈のコンプライアンス及び血圧を決定する方法及び装置
US10342437B2 (en) Detection of progressive central hypovolemia
US20100262022A1 (en) Detection of Progressive Central Hypovolemia using the System of the present invention with Pulse-Decomposition Analysis (PDA)
KR100951777B1 (ko) 혈액의 점도를 고려한 심장 모니터링 장치
US20240008748A1 (en) Method and apparatus for determining information about an arterial property of a subject
KR100585848B1 (ko) 말초혈관 혈류 측정을 이용한 비관혈적 혈압 측정 시스템
WO2008109185A2 (en) A noninvasive method to determine characteristics of the heart
JP2004105682A (ja) 自動脈波変動計測装置
Kolli et al. Optimization of balloon obstruction for simulating equivalent pressure drop in physiological stenoses
Karnaukhov et al. Developing a Pulse Wave Model Using the Lumped-Element Method
Gupta Blood Pressure Monitoring
RU2371084C1 (ru) Способ определения общей эластичности периферического сосудистого русла магистральных артерий
Yu et al. Oscillometric measurement of arterial pulse pressure for patients supported by a rotary blood pump
CN115645734A (zh) 体外反搏治疗冠心病的效果评价方法、系统及存储介质

Legal Events

Date Code Title Description
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

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080405

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090405

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100405

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110405

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120405

Year of fee payment: 10

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120405

Year of fee payment: 10