JPS6034430A - フロー測定方法および装置 - Google Patents
フロー測定方法および装置Info
- Publication number
- JPS6034430A JPS6034430A JP59078722A JP7872284A JPS6034430A JP S6034430 A JPS6034430 A JP S6034430A JP 59078722 A JP59078722 A JP 59078722A JP 7872284 A JP7872284 A JP 7872284A JP S6034430 A JPS6034430 A JP S6034430A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- flow
- cross
- function
- fluid
- 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
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/704—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
- G01F1/708—Measuring the time taken to traverse a fixed distance
- G01F1/712—Measuring the time taken to traverse a fixed distance using auto-correlation or cross-correlation detection means
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Fluid Mechanics (AREA)
- Biomedical Technology (AREA)
- Hematology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (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)
- Measuring Volume Flow (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
[発明の技術分野]
不発明龜、フロー測定に関するもので必ジ、とくにビー
ボ血液フロー測定に関するものでるる。
ボ血液フロー測定に関するものでるる。
この明細−書で使用する[フロー」という用語は、流体
の密度と6tI、体の速度との積で表される流体の動き
の測定値である。血液フロー測定Ci種々の医学的用途
をもつ重要な診断技術でろる。主な用途は正確なフロー
・レート測定が極めて有用となる診…1用の分析と心臓
の諸状態の処置におけるものである。
の密度と6tI、体の速度との積で表される流体の動き
の測定値である。血液フロー測定Ci種々の医学的用途
をもつ重要な診断技術でろる。主な用途は正確なフロー
・レート測定が極めて有用となる診…1用の分析と心臓
の諸状態の処置におけるものである。
70−を測定するため、とくに血管系における非圧縮性
流体としての血液の体積束すなわち平均フロー・レー)
k測定するために種々の技術が開発されている。そのよ
うな技術の一例が、殺菌した指示染料すなわち標識染料
全正確に測って心臓血管系中に断続的に、または連続し
て注入し、それ力)ら下流の測定点において濃度を検出
するという、指示薬稀釈技術でろる。フロー・レート6
:質量保存の原理を基にして計η、できる。しかし、そ
の手法は、系についていくつかの仮定を必要とし、しか
もその実施には時間が非′話にかかる。
流体としての血液の体積束すなわち平均フロー・レー)
k測定するために種々の技術が開発されている。そのよ
うな技術の一例が、殺菌した指示染料すなわち標識染料
全正確に測って心臓血管系中に断続的に、または連続し
て注入し、それ力)ら下流の測定点において濃度を検出
するという、指示薬稀釈技術でろる。フロー・レート6
:質量保存の原理を基にして計η、できる。しかし、そ
の手法は、系についていくつかの仮定を必要とし、しか
もその実施には時間が非′話にかかる。
フローを測定するために、測定場15iに挿入される測
定器装備のカテーテルを用いるという熱稀釈(tber
mod i 1ution)技術が知られている。スワ
ンおよびガング(Swan & Ganz)により開発
された4ルーメン・カテーテル(four −1urn
encatheter) では、高渦tたは低温の生理
的食塩水のようなある量の熱を、測定点の上流側に注入
する。測定点の下流側における温度の変化全温度計によ
り測定する。温度の時間的変化を解析することによりフ
ローが計算される。染料注入技術のように、このスワン
ーガンズ手法でも、殺菌された外部物質を用い、その動
作は本質的に断続的である。
定器装備のカテーテルを用いるという熱稀釈(tber
mod i 1ution)技術が知られている。スワ
ンおよびガング(Swan & Ganz)により開発
された4ルーメン・カテーテル(four −1urn
encatheter) では、高渦tたは低温の生理
的食塩水のようなある量の熱を、測定点の上流側に注入
する。測定点の下流側における温度の変化全温度計によ
り測定する。温度の時間的変化を解析することによりフ
ローが計算される。染料注入技術のように、このスワン
ーガンズ手法でも、殺菌された外部物質を用い、その動
作は本質的に断続的である。
閉じた系における流体フローを分散関数(disper
sion function) から計算できる。その
分散関数は、一般に、インパルス注入による指示物の投
与の結果として得られるインパルス応答関数でを)る。
sion function) から計算できる。その
分散関数は、一般に、インパルス注入による指示物の投
与の結果として得られるインパルス応答関数でを)る。
ステップ応答測定および正弦波応答測定のような他の決
定技術も、流体フロー全計算するために用いられている
。
定技術も、流体フロー全計算するために用いられている
。
上記の在米技術はともに質量保存の原理を基にしている
。それらの技術の全ての性質は基本的には決定論的なも
のであって、必ずしも真ではないことがあり得る仮定に
依存している。たとえば、指示(標識)物理量のフロー
が全流体のフローを表すと仮定し、および基準の点と測
定場所との間を流れるのに要する時間の分布は全流体の
要素と指示物理量の要素とに対して同じであると仮定し
ている。その上、指示物理量の体積は1つの入口ボート
と1つの出口ポートとの間で一定で必ると仮定し、測定
期間中は出口ボートη・ら入口ボートへの環流はないと
仮定している。更に、測定場所におけるるる量の指示物
理量の平均到達時間により平均流れ時間を計算できると
仮定している。
。それらの技術の全ての性質は基本的には決定論的なも
のであって、必ずしも真ではないことがあり得る仮定に
依存している。たとえば、指示(標識)物理量のフロー
が全流体のフローを表すと仮定し、および基準の点と測
定場所との間を流れるのに要する時間の分布は全流体の
要素と指示物理量の要素とに対して同じであると仮定し
ている。その上、指示物理量の体積は1つの入口ボート
と1つの出口ポートとの間で一定で必ると仮定し、測定
期間中は出口ボートη・ら入口ボートへの環流はないと
仮定している。更に、測定場所におけるるる量の指示物
理量の平均到達時間により平均流れ時間を計算できると
仮定している。
非決定論的信号励振を行うことにより、電子装置のイン
パルス応答全測定づ゛るための技哲が信号処理の分野に
おいて知られている。それらの技術ケよ、一般に通信の
用途においてり)み従来用いられていたものでめる。
パルス応答全測定づ゛るための技哲が信号処理の分野に
おいて知られている。それらの技術ケよ、一般に通信の
用途においてり)み従来用いられていたものでめる。
必要なものはフロー測定技術でめり、とくに、決定論的
測定技術に伴う不確定性かはとんとない液フロー測定技
術でるる。
測定技術に伴う不確定性かはとんとない液フロー測定技
術でるる。
し発明の概要〕
本発明によれば、確率的励振信号(stochastl
c eXeltajlon 81gnal ) k系の
入口に加え、その結果としてその系の下流側の出口に測
定加能な出力信号を発生させることKより、フローとく
に血管系におけるような不均一なフロー、が測定される
。励振信号と出力信号との相互相関(cross −c
orrelation) fc”解析1−ることにより
70−・レート’t−得ることができる。質量保存の原
理により較正を行うことができ、それにより量的なパラ
メータが濃度型パラメータに関連させられる。確率的励
振信号は白色雑音の性質金有するから、希望の情報を得
るために簡単にされた拡散スペクトル検出技術と信号抽
出技術を用いることができる。ある特定の実施例におい
ては、指示物理量の形で注入された信号は温度丑たはそ
れに類似のパラメータに変化を与えて、測定中の系の分
散関数に従って検出点まで流体のフローで伝えられる効
果を生じさせる。擬似ランダム2進シーケンス(pse
udo −random binary 5equen
ce)の形の好適な確率的入力信号が、検出点で受けら
れた信号と相互相関させられて、相互相関関数を表す出
力信号を得る。その出力信号力・らフロー・レートを計
算できる。とくに、相互相関関数があたかも系の分散関
数であるかのようにその相互相関関数が積分されて、そ
れ力)らフローを計算できるモーメント値を得られる。
c eXeltajlon 81gnal ) k系の
入口に加え、その結果としてその系の下流側の出口に測
定加能な出力信号を発生させることKより、フローとく
に血管系におけるような不均一なフロー、が測定される
。励振信号と出力信号との相互相関(cross −c
orrelation) fc”解析1−ることにより
70−・レート’t−得ることができる。質量保存の原
理により較正を行うことができ、それにより量的なパラ
メータが濃度型パラメータに関連させられる。確率的励
振信号は白色雑音の性質金有するから、希望の情報を得
るために簡単にされた拡散スペクトル検出技術と信号抽
出技術を用いることができる。ある特定の実施例におい
ては、指示物理量の形で注入された信号は温度丑たはそ
れに類似のパラメータに変化を与えて、測定中の系の分
散関数に従って検出点まで流体のフローで伝えられる効
果を生じさせる。擬似ランダム2進シーケンス(pse
udo −random binary 5equen
ce)の形の好適な確率的入力信号が、検出点で受けら
れた信号と相互相関させられて、相互相関関数を表す出
力信号を得る。その出力信号力・らフロー・レートを計
算できる。とくに、相互相関関数があたかも系の分散関
数であるかのようにその相互相関関数が積分されて、そ
れ力)らフローを計算できるモーメント値を得られる。
本発明の目的は、フローに関連する分散関数の誘導とい
う問題に対する確率的解決技術を得ることでるる。
う問題に対する確率的解決技術を得ることでるる。
本発明の別の目的は、インパルス応答を用いる技術、ス
テップ応答を用いる技術、または正弦波励振応答を用い
る技術のような、分散関数を得るための決定論的手法に
関連する仮定において本質的に生ずる誤差を避けること
である。
テップ応答を用いる技術、または正弦波励振応答を用い
る技術のような、分散関数を得るための決定論的手法に
関連する仮定において本質的に生ずる誤差を避けること
である。
本発明の更に別の目的は、ビーポ血液フロー測定のため
の技術を得ることでろる。
の技術を得ることでろる。
本発明全理解するためには、流れを測定する技術につい
て調べることが助けとなる。
て調べることが助けとなる。
質量の保存により、系に入る一ffiqすなわち制御量
は、フロー・レー)Fと時間に依存する濃度関数c(t
)との積の時間積分に等しい。すなわち、q= J″、
F−c(t) d t (1)分散関数というのは、
ぞJる濃度関数C(t)の2つの点の間の移動時間の分
布を記述する関数である。
は、フロー・レー)Fと時間に依存する濃度関数c(t
)との積の時間積分に等しい。すなわち、q= J″、
F−c(t) d t (1)分散関数というのは、
ぞJる濃度関数C(t)の2つの点の間の移動時間の分
布を記述する関数である。
すなわち、
ここeこh(t)は、再循環がないと仮定した時の分散
関数でるる。
関数でるる。
血管系における分散関数は、サンプリング場所の間にお
ける血管部分の正規化されたユニット・インパルス応答
である。その分散関数は、変形、分散および遅延全記述
するために使用でき、かつフロー全計算するためにも使
用できる。/ことえば周知の原理eこよジ、ある流体の
70−・レートは、ろる母集団の全ての要素の到達時間
で体積を除したものに等し、いことをボすことができる
。丁なわら、 q F = v/ t = −(3) Σh(t) ここに、■は体積、tは指示物理量の平均到達時間でめ
る。
ける血管部分の正規化されたユニット・インパルス応答
である。その分散関数は、変形、分散および遅延全記述
するために使用でき、かつフロー全計算するためにも使
用できる。/ことえば周知の原理eこよジ、ある流体の
70−・レートは、ろる母集団の全ての要素の到達時間
で体積を除したものに等し、いことをボすことができる
。丁なわら、 q F = v/ t = −(3) Σh(t) ここに、■は体積、tは指示物理量の平均到達時間でめ
る。
この平均到達時間は、るる到達点におけるある体積の母
集団分布全記述する関数の第1のモーメントすなわち積
分を表す。この例においては、分散関数h(t)は母集
団分布でるる。したがって、分散関数をめることにより
、その他の仮定に関連してフローを計算できる。
集団分布全記述する関数の第1のモーメントすなわち積
分を表す。この例においては、分散関数h(t)は母集
団分布でるる。したがって、分散関数をめることにより
、その他の仮定に関連してフローを計算できる。
本発明によれば、拡散スペクトルまたは既知の確率的励
振信号’t6る系の入口点に刀lえ、その系の出口点に
おいて時間に依存する応答を検出し、励振信号を応答信
号に相互相関させて分散関数を構成することにより、フ
ローに対する分1枚関数力玉求められる。そして、較正
ノくラメータを参考にして、フローを計算できる。
振信号’t6る系の入口点に刀lえ、その系の出口点に
おいて時間に依存する応答を検出し、励振信号を応答信
号に相互相関させて分散関数を構成することにより、フ
ローに対する分1枚関数力玉求められる。そして、較正
ノくラメータを参考にして、フローを計算できる。
分散関数(これは伝達関数とも呼tまれる)hは、ウイ
ーt−=ホップ7 (Wi ener= 1(opf
)の式を用いて得ることがで、きる。そのウイーナー=
ホップ70式は、める系の入ブjとインノくルス応答す
なわち伝達関数との自己相関のコンボ1ツユ−ジョンが
、入力と出力との相互相関に等しいということを表すも
のでめる。すなわち、 φxy(t) =h(v)φX z (t −v )
d v (4)ここに φXy(t) (4a ) は入力と出力の相互相関で、1、 φXX(t)(4b) は入力の自己相関でろり、 h (v)(4c ) は分散関数、伝達関数またはインパルス応答関数でろる
。
ーt−=ホップ7 (Wi ener= 1(opf
)の式を用いて得ることがで、きる。そのウイーナー=
ホップ70式は、める系の入ブjとインノくルス応答す
なわち伝達関数との自己相関のコンボ1ツユ−ジョンが
、入力と出力との相互相関に等しいということを表すも
のでめる。すなわち、 φxy(t) =h(v)φX z (t −v )
d v (4)ここに φXy(t) (4a ) は入力と出力の相互相関で、1、 φXX(t)(4b) は入力の自己相関でろり、 h (v)(4c ) は分散関数、伝達関数またはインパルス応答関数でろる
。
(4)式は、任意の自己相関に対して真でめる一般式て
める。自己相関関数を適切に選択することにより、すな
わち、自己相関関数が白色雑音であるように作られてい
るものとすると、ウイーナー=ホツプフの式は簡単では
めるが重要な彫りy(t) = K h(t) + B
(5)となる。ここに、Kは比例定数、Bはオフセッ
トでろる。この簡単な形は、零に等しい時間値において
のみ自己相関値を有するという性質の結果である。
める。自己相関関数を適切に選択することにより、すな
わち、自己相関関数が白色雑音であるように作られてい
るものとすると、ウイーナー=ホツプフの式は簡単では
めるが重要な彫りy(t) = K h(t) + B
(5)となる。ここに、Kは比例定数、Bはオフセッ
トでろる。この簡単な形は、零に等しい時間値において
のみ自己相関値を有するという性質の結果である。
この簡単にさ%fしたウイーナー=ホツプフの式(式(
5))は、そのインパルス応答h (t)がめられてい
る系でろって、かつろる有限の時間で相互相関関数と自
己相関関数の適切な評価を行える系につ武 いての、適切な帯域幅を有する擬似ランダム雑音源によ
り実現できる。適当な励振源は、測定開始前の少くとも
安定化シーケンス期間中に連続して加えられる2進最大
長シーケンスである。その2進最大長シーケンスは、次
のような性質を一督する一連の2進数でろる。すなわち
、n個の個5j1」時間間隔に対して、そのシーケンス
はN個の個、MIJ時間間隔すなわちデジットで形成さ
れる。ここに、N=2−1 (6) でるり、 (N+1)/2 (7) 論理l状態と、 (N−1)/2 (8) 論理0状態とを含む。2進最大長シーケンスの状態と特
性とについての別の情報力(デービイス(W。
5))は、そのインパルス応答h (t)がめられてい
る系でろって、かつろる有限の時間で相互相関関数と自
己相関関数の適切な評価を行える系につ武 いての、適切な帯域幅を有する擬似ランダム雑音源によ
り実現できる。適当な励振源は、測定開始前の少くとも
安定化シーケンス期間中に連続して加えられる2進最大
長シーケンスである。その2進最大長シーケンスは、次
のような性質を一督する一連の2進数でろる。すなわち
、n個の個5j1」時間間隔に対して、そのシーケンス
はN個の個、MIJ時間間隔すなわちデジットで形成さ
れる。ここに、N=2−1 (6) でるり、 (N+1)/2 (7) 論理l状態と、 (N−1)/2 (8) 論理0状態とを含む。2進最大長シーケンスの状態と特
性とについての別の情報力(デービイス(W。
D、 T、 Davies )著[自己j轟応III御
c7) fc メ(7)系識別(System Ide
ntificationfor 5eif −Adap
tive Control ) J (NY: Joh
nWiley & 5ons、 1970)に見られる
。
c7) fc メ(7)系識別(System Ide
ntificationfor 5eif −Adap
tive Control ) J (NY: Joh
nWiley & 5ons、 1970)に見られる
。
系のインパルス応答関数h(t、) k得るために相互
相関させる目的でるる系への入力として2進最大長シー
ケンスが用いらする場合には、(4)入力)らのコンボ
リューション積分は次の形となる。
相関させる目的でるる系への入力として2進最大長シー
ケンスが用いらする場合には、(4)入力)らのコンボ
リューション積分は次の形となる。
2進シーケンスにおいては、乗算の代りに加算金使用で
きるから、実時間演算のために相互相関過程を簡単にで
きる。
きるから、実時間演算のために相互相関過程を簡単にで
きる。
系のインパルス応答関数h(t) ?適切に識別するf
cめに(r−t、、シーケンス長の周期T=N△を全、
ぞの糸のインパルスがほぼ零まで減衰するのに要する時
間ニジ長くしなければならない。さもないと、スペクト
ルの折9曲りが起き、そのために不正確な系関数が生ず
ることになる。
cめに(r−t、、シーケンス長の周期T=N△を全、
ぞの糸のインパルスがほぼ零まで減衰するのに要する時
間ニジ長くしなければならない。さもないと、スペクト
ルの折9曲りが起き、そのために不正確な系関数が生ず
ることになる。
同じ種類の制約が周波数領域に対してもめては凍る。
流体のフローとくに血液のフローの計算は、系の関数h
(t)の時間積分が定義により1に等しいことに注意
して、(9)式から行うことができる。個別の形におい
ては、系の関数h(t)は個別値h(kΔt)に一致す
る。
(t)の時間積分が定義により1に等しいことに注意
して、(9)式から行うことができる。個別の形におい
ては、系の関数h(t)は個別値h(kΔt)に一致す
る。
ユニット励振を仮定すると、相互相関関数φxyは係数
aによジ換算′きれる。ここに、aは濱1j定すべきフ
ロー系の中枢側(入口)月?−トにおける温度上昇また
は指示物理量の増加に等しい。Pffi入カパフカパワ
一時間当りの熱、または単位時間当りの指示物理量の注
入′@:)、Cを媒体(すな増)ち、血液の比熱または
媒体中の指示物理量の量の上ヒ密度、Fをフロm−レー
トとして、 P (lo) 8−]下 を認めると、系に与えられる2進最大長シーケンス励振
の相互相関関数の値を測定づ−ることVこよシフロー・
シー1−Ffz決定できることが示さねる。
aによジ換算′きれる。ここに、aは濱1j定すべきフ
ロー系の中枢側(入口)月?−トにおける温度上昇また
は指示物理量の増加に等しい。Pffi入カパフカパワ
一時間当りの熱、または単位時間当りの指示物理量の注
入′@:)、Cを媒体(すな増)ち、血液の比熱または
媒体中の指示物理量の量の上ヒ密度、Fをフロm−レー
トとして、 P (lo) 8−]下 を認めると、系に与えられる2進最大長シーケンス励振
の相互相関関数の値を測定づ−ることVこよシフロー・
シー1−Ffz決定できることが示さねる。
づ−なわぢ、
ΣφXy(kΔL)
1(=1
(11)式は、フローを相互相関関数に亘接関連づける
ものでるるから、非常に重要な式でるる。この例では、
相互相関関数は、擬似ランタ゛ム雑音自己相関関数との
インパルス応答関数のコンボリューションである。選択
された入力により信号処理を行うことができ、それによ
り糸の入口における信号は系の出口において検出された
信号と相互相関させられる。その結果を積分、具体的に
は加え合わせることができ、その積分の逆数に定数を乗
することによりフロー・レートを表す値が得られる。そ
の値が本発明に従ってめられるパラメータである。
ものでるるから、非常に重要な式でるる。この例では、
相互相関関数は、擬似ランタ゛ム雑音自己相関関数との
インパルス応答関数のコンボリューションである。選択
された入力により信号処理を行うことができ、それによ
り糸の入口における信号は系の出口において検出された
信号と相互相関させられる。その結果を積分、具体的に
は加え合わせることができ、その積分の逆数に定数を乗
することによりフロー・レートを表す値が得られる。そ
の値が本発明に従ってめられるパラメータである。
本発明の特定の実施例においてケよ、指示物理量は温度
レベルの間で変化できる加熱要素により加えられる熱で
ある。本発明の別の実施例においては、指示物理量は、
流体注入シーケンスにおいてオン−オフの変化をさせて
注入できる駆動要素により注入される染料でるる。
レベルの間で変化できる加熱要素により加えられる熱で
ある。本発明の別の実施例においては、指示物理量は、
流体注入シーケンスにおいてオン−オフの変化をさせて
注入できる駆動要素により注入される染料でるる。
以下、図面を参照して本発明の詳細な説明する。
図には、血管系のような流体系12の中枢側ボート14
と末梢側ボート16の間を流れる流体のフローを測定す
るための測定装置10の代表的な実施例が示されている
。熱のような指示物理量Q(t)だけ加熱要素18によ
って中枢(IIIボート14に導入される。温度のよう
な指示物理量の濃度の変1Lが末梢側ボート16におい
てトランスデユーサ20により検出される。るる特定の
実施例においては、加熱要素18は、フローにより与力
ふれる風船を先端部にと9つけら+’Lfc肺動脈カテ
ーテル(flow−directed balloon
−tipped pulmonary artery
catheter)にと9つり゛られた精密級の高抵
抗ニクロム加熱線とすること75玉でき、トランスデユ
ーサ20としてほそ・D同じカテーテルの内部にとりつ
けられたサーミスタで構成できる。そのカテーテルの構
造はスワンーガンズの熱稀釈カテーテルの構造に類似す
る。
と末梢側ボート16の間を流れる流体のフローを測定す
るための測定装置10の代表的な実施例が示されている
。熱のような指示物理量Q(t)だけ加熱要素18によ
って中枢(IIIボート14に導入される。温度のよう
な指示物理量の濃度の変1Lが末梢側ボート16におい
てトランスデユーサ20により検出される。るる特定の
実施例においては、加熱要素18は、フローにより与力
ふれる風船を先端部にと9つけら+’Lfc肺動脈カテ
ーテル(flow−directed balloon
−tipped pulmonary artery
catheter)にと9つり゛られた精密級の高抵
抗ニクロム加熱線とすること75玉でき、トランスデユ
ーサ20としてほそ・D同じカテーテルの内部にとりつ
けられたサーミスタで構成できる。そのカテーテルの構
造はスワンーガンズの熱稀釈カテーテルの構造に類似す
る。
加熱要素18は、約1o o ミIJ秒より短い時間中
に零状態から約2〜10ワツトの連続電力状態へ電流を
迅速に切換えることができる電子的増幅器22のような
電源により、加熱づれる。この増幅器22は確率的発振
器24により励振される。
に零状態から約2〜10ワツトの連続電力状態へ電流を
迅速に切換えることができる電子的増幅器22のような
電源により、加熱づれる。この増幅器22は確率的発振
器24により励振される。
その確率的発振器24は白色雑音に似た擬(以ランダム
2進パルス列またはランダム2ン篤長シーケンスを発生
する。
2進パルス列またはランダム2ン篤長シーケンスを発生
する。
トランスデユーサ2oから信号を受けるために筒利得検
出増幅器26が結合される。その検出増幅器26φj;
、センサ出カナなゎぢサーミスタ出力のマイクロボルト
a度の変化全検出できる。その変化は末梢側ポート16
における指示物理量濃度の小さな変1は示すものである
。演算処理の便宜のため、検出増幅器26の出方はアナ
ログ−デジタル変換器(ADC)28によりデジタル形
式に変換される。
出増幅器26が結合される。その検出増幅器26φj;
、センサ出カナなゎぢサーミスタ出力のマイクロボルト
a度の変化全検出できる。その変化は末梢側ポート16
における指示物理量濃度の小さな変1は示すものである
。演算処理の便宜のため、検出増幅器26の出方はアナ
ログ−デジタル変換器(ADC)28によりデジタル形
式に変換される。
確率的発振器24により発生さitた信号X(△t)と
AL)02B の出方端子に生じた信号y(Δt〕は相
互相関器30へ与えられる。その相互相関器は並列接続
されている一群の遅延器31と、乗算器33と、積分器
35とを備える。それらの積分器35trJ、、イ=号
y全乗ぜられた、N個の遅延させられた人力2進最大長
シーケンスを、並列に処理する。、相互相関器3oの出
方は、N個の信号φxy(k△t)である。相互相関器
3oは、相互相関関数φxy(kΔt)の全てのもの全
同時に計算できるように構成される。したがって、相互
相関器30の出力は、直列時間から並列時間に変換され
た所定ノシーケンス長Tの並列的にサンプリングした信
号と拷えることができる。シーケンス長Tは、インパル
ス応答関数h(k△t) (k−’=1・・・。
AL)02B の出方端子に生じた信号y(Δt〕は相
互相関器30へ与えられる。その相互相関器は並列接続
されている一群の遅延器31と、乗算器33と、積分器
35とを備える。それらの積分器35trJ、、イ=号
y全乗ぜられた、N個の遅延させられた人力2進最大長
シーケンスを、並列に処理する。、相互相関器3oの出
方は、N個の信号φxy(k△t)である。相互相関器
3oは、相互相関関数φxy(kΔt)の全てのもの全
同時に計算できるように構成される。したがって、相互
相関器30の出力は、直列時間から並列時間に変換され
た所定ノシーケンス長Tの並列的にサンプリングした信
号と拷えることができる。シーケンス長Tは、インパル
ス応答関数h(k△t) (k−’=1・・・。
N)がほぼ零まで減少できるように十分に長く選択され
る。
る。
相互相関器30の出力は、加算器32へ与えられる。こ
の加算器tよ、相互相関関数φXyの全てのものを加え
合わせることにより、前記した(11)式の分母の値を
得る積分器として動作する。
の加算器tよ、相互相関関数φXyの全てのものを加え
合わせることにより、前記した(11)式の分母の値を
得る積分器として動作する。
分散関数のための基準線を得ることが通常望ましい。そ
のために、相互相関器N をオフセットして片寄りを除
去する。とくvご、1回またはそれ以上の動作にわたっ
て相互相関器30の出力を入力として受け、(5)式に
示す基準線オフセット係数Bを出力として生ずる応答関
数正規化器34が設けられる。この応答関数正規化器は
、まず応答関数h(k△t)が1つの減衰する指数関数
過程を表すと仮定することにより、(5)式を基にして
Bを計算するように動作する。それから、応答関数正規
化器34は、インパルス応答関数の以前のピーク部分を
最も良い指数関数曲線に適合させるために、最小平均二
乗誤差技術(1east means 5quared
error technique)’fc用いる。最
も良く適合いで指数関数曲線に最も良く適合するそのイ
ンパルス応答関数は、基準線値Bを計算フロるために用
いられる 加算器32は、並列相互相関関数の和の形の出力をメモ
リ36に与える。そのメモリには多数回の試験?行って
得られたデータが格納される。メモリ36はフロー・レ
ート計算器38に結合される。フロー・レート計算器3
8は、(11)式の簡単な代数計算ケ行い、メモリ36
から相互相関関数をと9出し、オフセット係数を応答関
数正規化器34からと9出し、パワーPの値と、比AC
の値と、シーケンス長Nの値とを外部の、典型的には手
動の、入カ装fit40から得る。フロー・レート計算
器38の出方端子は適当なフロー・レート指示器39に
接続される。このフロー・レート指示器は、2ボート系
でろると仮定して、中枢側ポート14と末梢側ポート1
6を有する系を流れるフロー・レートの希望の値を表示
する。
のために、相互相関器N をオフセットして片寄りを除
去する。とくvご、1回またはそれ以上の動作にわたっ
て相互相関器30の出力を入力として受け、(5)式に
示す基準線オフセット係数Bを出力として生ずる応答関
数正規化器34が設けられる。この応答関数正規化器は
、まず応答関数h(k△t)が1つの減衰する指数関数
過程を表すと仮定することにより、(5)式を基にして
Bを計算するように動作する。それから、応答関数正規
化器34は、インパルス応答関数の以前のピーク部分を
最も良い指数関数曲線に適合させるために、最小平均二
乗誤差技術(1east means 5quared
error technique)’fc用いる。最
も良く適合いで指数関数曲線に最も良く適合するそのイ
ンパルス応答関数は、基準線値Bを計算フロるために用
いられる 加算器32は、並列相互相関関数の和の形の出力をメモ
リ36に与える。そのメモリには多数回の試験?行って
得られたデータが格納される。メモリ36はフロー・レ
ート計算器38に結合される。フロー・レート計算器3
8は、(11)式の簡単な代数計算ケ行い、メモリ36
から相互相関関数をと9出し、オフセット係数を応答関
数正規化器34からと9出し、パワーPの値と、比AC
の値と、シーケンス長Nの値とを外部の、典型的には手
動の、入カ装fit40から得る。フロー・レート計算
器38の出方端子は適当なフロー・レート指示器39に
接続される。このフロー・レート指示器は、2ボート系
でろると仮定して、中枢側ポート14と末梢側ポート1
6を有する系を流れるフロー・レートの希望の値を表示
する。
以上説明した本発明の装E<’c作って試験看:行った
。血管系を表す試@部分は、混合室と、右心臓・肺動脈
に数学的VC類似する寸つ丁ぐな部分とで構成された血
管系で、ちる。体積ンシ)@接測定され、フィッシャー
拳アンドΦポーター(Fisher &Porter)
社の10A11027A型流量計全用いて実際の流量を
測定した。注入量を脈動しない流れ、すなわち、一定の
流れの下に測定した。冷水と、一定水頭タンクと、−犀
抵抗値り)逃し弁とを用いて、諸条件を一定に保持した
。相互相関器をシミ:L シ/ −i−するためNi
PDP−11ミニコンピユータをプログラムした。(た
とえばヒユーレット・パラカード社(I−Iewlet
t Packerd Company ) により製造
・販売され1いる装置のような市販の相互相関器も使用
できる。) この試験の結果、測定しグ、−フロー・レートb″、広
い実験したフロー・レートの範囲にわたって、真のフロ
ー・レートとめまり大きくは異ならないことか判明して
いる。偏差は、装置の構成部品が理想的なものでないこ
と、雑音を避けられないCと、量子化誤差が生ずること
、相互相関関数φxy(1)がインパルス応答関数h(
t)に正比例するという仮定に起因する二次的な影響な
どにより生ずるものでるる。系の固有の雑音が高次の成
分の太きさより大きいために、相互相関関数と系の関数
が等しいという理論的仮定は正当化される。
。血管系を表す試@部分は、混合室と、右心臓・肺動脈
に数学的VC類似する寸つ丁ぐな部分とで構成された血
管系で、ちる。体積ンシ)@接測定され、フィッシャー
拳アンドΦポーター(Fisher &Porter)
社の10A11027A型流量計全用いて実際の流量を
測定した。注入量を脈動しない流れ、すなわち、一定の
流れの下に測定した。冷水と、一定水頭タンクと、−犀
抵抗値り)逃し弁とを用いて、諸条件を一定に保持した
。相互相関器をシミ:L シ/ −i−するためNi
PDP−11ミニコンピユータをプログラムした。(た
とえばヒユーレット・パラカード社(I−Iewlet
t Packerd Company ) により製造
・販売され1いる装置のような市販の相互相関器も使用
できる。) この試験の結果、測定しグ、−フロー・レートb″、広
い実験したフロー・レートの範囲にわたって、真のフロ
ー・レートとめまり大きくは異ならないことか判明して
いる。偏差は、装置の構成部品が理想的なものでないこ
と、雑音を避けられないCと、量子化誤差が生ずること
、相互相関関数φxy(1)がインパルス応答関数h(
t)に正比例するという仮定に起因する二次的な影響な
どにより生ずるものでるる。系の固有の雑音が高次の成
分の太きさより大きいために、相互相関関数と系の関数
が等しいという理論的仮定は正当化される。
動物実験も行った。熱励振要素として力ロ熱要素をと9
つけるタメに、スワン=ガンズ型のフローにより4かれ
る肺動脈力テーブルを改装した。麻酔をかけrc犬に大
腿動脈を通してカテーテル全挿入した。肺動脈の中にサ
ーミスタ検出器を挿入した。内部大静脈(1nteri
or vena cava)内に加熱要素を挿入した。
つけるタメに、スワン=ガンズ型のフローにより4かれ
る肺動脈力テーブルを改装した。麻酔をかけrc犬に大
腿動脈を通してカテーテル全挿入した。肺動脈の中にサ
ーミスタ検出器を挿入した。内部大静脈(1nteri
or vena cava)内に加熱要素を挿入した。
そうすると、確率的入力により7JI′I熱要素を加熱
し、時間を基準にして出力温度を監視し、本発明に従っ
て相互相関を計算することにより心臓の出力を決定した
。
し、時間を基準にして出力温度を監視し、本発明に従っ
て相互相関を計算することにより心臓の出力を決定した
。
本発明に従って作られた実用的な装置に誤差が生ずるこ
とについて少し述べることVcjる。第1に、本発明は
、理論的に完全な入力2進最大長シ一ケンス関数信号、
すなわち、ただ2Qの状態が存在し、それら2つの状態
の間のスイッチング時間が無限に短いと仮定する関数信
号、を仮定する。
とについて少し述べることVcjる。第1に、本発明は
、理論的に完全な入力2進最大長シ一ケンス関数信号、
すなわち、ただ2Qの状態が存在し、それら2つの状態
の間のスイッチング時間が無限に短いと仮定する関数信
号、を仮定する。
突際には、状態遷移時間を要求される。その状態遷移時
間は自己相関関数に誤差を生じさせると考えられる。し
かし、クロック時間と比較してスイッチング時間が短い
とすると、励振時には大きな誤差が生ずることはない。
間は自己相関関数に誤差を生じさせると考えられる。し
かし、クロック時間と比較してスイッチング時間が短い
とすると、励振時には大きな誤差が生ずることはない。
しかし、加熱要素の時定数の結果として無視でキナい立
上り時間と立下9時間が存在する。し〃)シ、加熱要素
は01は等しい温度上昇時間と温度低下時間を示す力・
ら、修正を行うことは不要でろる。
上り時間と立下9時間が存在する。し〃)シ、加熱要素
は01は等しい温度上昇時間と温度低下時間を示す力・
ら、修正を行うことは不要でろる。
同様に、出カドランスデューサは有限の時定数を示す。
試作した装置で用いたサーミスタの助定数は0.1秒台
でめった。しかし、この時定V、は0.5秒と1秒の間
のクロック時間間隔と比較して十分に短いから、データ
は(互とんど影響を受けなかった。時定数がサンプル時
間と比較して長くさ4ている力・ら、高い周波数の情報
が除去される。系のインパルス応答関数h(t) は比
較的少い高周波数成分を有する応答でろると仮定してお
り、かつ高い周波数成分が系のインパルス応答の一部で
はない二次の影響る表子から、出カドランスデューサの
時定数がクロック時間より長くなると、望ましくない高
い周波数成分を無くす効果的かつ望ましい低域フィルタ
関数が得られる。
でめった。しかし、この時定V、は0.5秒と1秒の間
のクロック時間間隔と比較して十分に短いから、データ
は(互とんど影響を受けなかった。時定数がサンプル時
間と比較して長くさ4ている力・ら、高い周波数の情報
が除去される。系のインパルス応答関数h(t) は比
較的少い高周波数成分を有する応答でろると仮定してお
り、かつ高い周波数成分が系のインパルス応答の一部で
はない二次の影響る表子から、出カドランスデューサの
時定数がクロック時間より長くなると、望ましくない高
い周波数成分を無くす効果的かつ望ましい低域フィルタ
関数が得られる。
定常状態利得といわゆる白色雑音が、系にオフセット誤
差を生じさせることが予測される。本発明の装置により
実行される計算は、面積またはモーメンIf得るために
系のインパルス応答関数を積分1゛ることに等しいから
、インパルス応答関数h(t)のための真の基準線を得
る要素を設けなけれ。
差を生じさせることが予測される。本発明の装置により
実行される計算は、面積またはモーメンIf得るために
系のインパルス応答関数を積分1゛ることに等しいから
、インパルス応答関数h(t)のための真の基準線を得
る要素を設けなけれ。
ばならない。前記した応答関数正規化器34のような要
素により曲線に適合する指数関数を使用する技術が最も
効果的な技術であることが判明している。
素により曲線に適合する指数関数を使用する技術が最も
効果的な技術であることが判明している。
この装置によるフロー測定の精度は、測定すべき系にお
けるフローの変動に左右される刀・ら、正確な平均フロ
ー・レートを得ることは困難となる。
けるフローの変動に左右される刀・ら、正確な平均フロ
ー・レートを得ることは困難となる。
この目的のために、データを得るための測定を、何回か
実行するのが一般的でβる。1つの標本に対する出力フ
ロー・レー+−存生ずるためにフロー・レート計算器全
使用でき、廿たは何回かの測定にわたってフロー・レー
トを平均するために希望に応じてフロー・レート計算器
を使用できる。この方法は、平均フロー・レートヲ−層
正確に決定するために積分時間を長くすることか、1つ
の標本についての測定を長く行うことシこ大体等しい。
実行するのが一般的でβる。1つの標本に対する出力フ
ロー・レー+−存生ずるためにフロー・レート計算器全
使用でき、廿たは何回かの測定にわたってフロー・レー
トを平均するために希望に応じてフロー・レート計算器
を使用できる。この方法は、平均フロー・レートヲ−層
正確に決定するために積分時間を長くすることか、1つ
の標本についての測定を長く行うことシこ大体等しい。
ッ、上、特定の実施例について本発明全説明したが、そ
の他の実施例も当業者には明らかでめろう。
の他の実施例も当業者には明らかでめろう。
たとえば、インライン流量計励低の挿入に適さない装置
において、直接的流体フロー測定が実用的でない場合に
は、常に本発明を用いることができる。可変熱入力以外
の手段、たとえば、2進最大長シーケンス励振パターン
の形の指示染料および放射性トレーサのような挿入さ′
#′1.なし1励振要素を使用できる。それらの指示体
をよ、それらに固有の欠点を示すが、同様に良く機能で
きる。本発明の装置の潜在的な用途は、多ボート流体系
の解析に才・・けるものである。予測さ才りる各出口ボ
ートにトランスデユーサを設けることができ、入力ボー
トにおける平均フロー・レートヲより正確に知るために
信号出力をいくつ刀島のディメンジョンで相互相関させ
ることができる。更に、本発明は自動連続心臓フロー測
定のために使用でき、心臓の鎖状態の処置における潜在
的に重装な診断手段となる。
において、直接的流体フロー測定が実用的でない場合に
は、常に本発明を用いることができる。可変熱入力以外
の手段、たとえば、2進最大長シーケンス励振パターン
の形の指示染料および放射性トレーサのような挿入さ′
#′1.なし1励振要素を使用できる。それらの指示体
をよ、それらに固有の欠点を示すが、同様に良く機能で
きる。本発明の装置の潜在的な用途は、多ボート流体系
の解析に才・・けるものである。予測さ才りる各出口ボ
ートにトランスデユーサを設けることができ、入力ボー
トにおける平均フロー・レートヲより正確に知るために
信号出力をいくつ刀島のディメンジョンで相互相関させ
ることができる。更に、本発明は自動連続心臓フロー測
定のために使用でき、心臓の鎖状態の処置における潜在
的に重装な診断手段となる。
添附図面は本発明に従って70−を測定する装置のブロ
ック図である。 10・・・・測定装置、12・・・・流体系、18・・
・・〃口前要素、20・・・・トランスデユーサ、22
・・・・増幅器、24・・・・確率的発振器、26・・
′・・検出増幅器、28・・・曽アナログーデジタル変
換器、30・・・・相互相関器、31・・・會遅延器、
32・・・・メモリ、33・・・・乗算器、34・・φ
・応答関数正規化器、35佛・・・積分器、36・・轡
・メモリ、38ψ・・・フロー吻し−)計算!、39・
・・・フロー・レート指示器。 特許出願人 ザ・ボード・オブ・トラステイーズ・オプ
φザ・リーランド・スクンフオー トーシュニア・ユニヴアーシティ 代理人 山川政樹(ほか2名)
ック図である。 10・・・・測定装置、12・・・・流体系、18・・
・・〃口前要素、20・・・・トランスデユーサ、22
・・・・増幅器、24・・・・確率的発振器、26・・
′・・検出増幅器、28・・・曽アナログーデジタル変
換器、30・・・・相互相関器、31・・・會遅延器、
32・・・・メモリ、33・・・・乗算器、34・・φ
・応答関数正規化器、35佛・・・積分器、36・・轡
・メモリ、38ψ・・・フロー吻し−)計算!、39・
・・・フロー・レート指示器。 特許出願人 ザ・ボード・オブ・トラステイーズ・オプ
φザ・リーランド・スクンフオー トーシュニア・ユニヴアーシティ 代理人 山川政樹(ほか2名)
Claims (1)
- 【特許請求の範囲】 (1,)、Mの入口点と系の出口点の間の流体のフロー
を測足フーる方法でるって、既知の確率的励振信号の形
の指示物理量を系の前fjL入口点において流体に加え
る過程と、 系の前記出口点において前記励振信号Vこ対する前記流
体の時間に依存する応答を検出して応答信号を得る過程
と、 前記励振信号をntJ記応答信号に相互相関させて相互
相関関数信号全構成する過程と、 υiJ it 4fj互al関関数信号および前記励振
4g号の指示物理の最の情報からフロー・レートを表子
イば号をとり出す過程と、 フローを表す前記信号を信号利用装置へ与える過程と、 全1+iiiえることを特徴とTるフロー測定方法。 (2)特i、V(′請求の範囲第1項記載の方法であっ
て、前記記とり出し過程は、Fを70−金表す信号、P
を単位時間当りの指示物理量の注入量、Cを前記流体中
における指示物理量の比畝匹、I’l前記確率的励振信
号にお(ハ)る間隔の数、l\tfa:クロック間隔、
φxYを相互相関関数、■<を・整数として、式 に従ってフロー全表す信号荀発生1゛る過程を含むこと
を特徴と1−るフロー測定方法。 (3)%許請求の範囲第1項lた社第2項記載の方法で
めって、前記確率的励振43号は白色雑廿源を近似する
ことを特徴とするフロー測定方法。 (4)特許請求の範囲第1項、第2項、または、第3項
記載の方法でろって、前記確率的励振信号は2進最大長
シーケンスであることを特徴とするフロー測定方法。 (5)特許請求の範囲第1項、第2項、第3項、−!た
は、第4項記載の方法でろって、前記加える過程は、前
記流体に温度変化を与えるために熱を加える過程を含む
ことを特徴とするフロー測定方法。 (6)特許請求の範囲第1項、第2項、第3項、また&
J:、第4項記載の方法であって、前記加える過程は、
前記流体中に、指示物理量を、その濃度が変化するよう
にして入れる過程を含むことを特徴とするフロー611
j定方法。 (7)特許請求の範囲第1項、第2項、第3項、第4項
、第5項、または、第6項記載の方法であって、ntJ
記相互相関関数信号を指数減衰関数として近似すること
によジ、前記相互相関過程のための基準線を設定する過
程と、前記基準線として前記指数減衰関数の漸近線に対
応する値全選択する過程とを含むCとを特徴とするフロ
ー測定方法。 (8)フロー・シー1測定すべき系の入口ボートに既知
の硲率的励振信号の形で指示物理量を加える斤めの要素
と、 前記系の出口ポートに配性され、前記系の前記入口ボー
トにおける前記励振信号に応答して前記系の流体パラメ
ータの変化を検出し、検出された信号を発生するように
動作する要素と、前記励振信号と検出された前記信号を
相互相関させて、前記系の特性インノクルス応答にほぼ
i線型に関連させられる相互相関関数信号を得るための
要素と、 複数の相互相関関数信号全積分して相互相関関数信号の
和を発生するための要素と、 相互相関関数信号の前記和に応答して、測定されたフロ
ー・レートを表す系の出力S号音発生するための要素と
、 全備えるフロー測定装置。 (9)%許請求の範囲第8項記載の装置であって、前記
信号を与える要素は、白色雑音源を近似する擬似ランダ
ム信号を発生するように動作することを特徴とするフロ
ー測定装置。 (10)特許請求の範囲第9項記載の装置でろって、前
記信号を与える要素は、2進最大長シーケンスの形で前
記確率的励振信号を発生するように動作することを特徴
とするフロー測定装置。 (11)特許請求の範囲第8項、第9項、または、第1
O項記載の装置であって、前記系を較正するための要素
を含み、この較正要素は、前記相互相関関数信号を指数
減衰関数と(7て近似するkめの要素と、前記指数減衰
関数の漸近線に対応する値を基準線と[7て選択するた
めの要素とを含むことを%徴とするフロー測定装部−8
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US487334 | 1983-04-21 | ||
US06/487,334 US4507974A (en) | 1983-04-21 | 1983-04-21 | Method and apparatus for measuring flow |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS6034430A true JPS6034430A (ja) | 1985-02-22 |
JPH0527051B2 JPH0527051B2 (ja) | 1993-04-20 |
Family
ID=23935317
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP59078722A Granted JPS6034430A (ja) | 1983-04-21 | 1984-04-20 | フロー測定方法および装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US4507974A (ja) |
EP (1) | EP0126931B1 (ja) |
JP (1) | JPS6034430A (ja) |
CA (1) | CA1210509A (ja) |
DE (1) | DE3483626D1 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63216536A (ja) * | 1987-03-05 | 1988-09-08 | テルモ株式会社 | 心拍出量の測定装置 |
JPS63216539A (ja) * | 1987-03-05 | 1988-09-08 | テルモ株式会社 | 心拍出量の相対的変化の連続測定装置 |
JPS63216537A (ja) * | 1987-03-05 | 1988-09-08 | テルモ株式会社 | 心拍出量の測定装置 |
JPS63216540A (ja) * | 1987-03-05 | 1988-09-08 | テルモ株式会社 | 心拍出量の測定装置 |
JPS63286165A (ja) * | 1987-05-01 | 1988-11-22 | プロダクト・イノウベイション・ホールディングズ・リミテッド | 医学注入装置 |
US7688098B2 (en) | 2007-03-08 | 2010-03-30 | Nec Corporation | Power supply noise measuring circuit and power supply noise measuring method |
Families Citing this family (70)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4933633A (en) * | 1981-06-09 | 1990-06-12 | Adec, Inc. | Computer controlled energy monitoring system |
DE3642771A1 (de) * | 1986-12-15 | 1988-06-23 | Kuipers Ulrich | Verfahren und vorrichtung zur messung der messgroesse eines messobjekts |
US4785823A (en) * | 1987-07-21 | 1988-11-22 | Robert F. Shaw | Methods and apparatus for performing in vivo blood thermodilution procedures |
DE3725312A1 (de) * | 1987-07-30 | 1989-02-09 | Jiri Hokynar | Steuergeraet fuer fluidfluss |
US5146414A (en) * | 1990-04-18 | 1992-09-08 | Interflo Medical, Inc. | Method and apparatus for continuously measuring volumetric flow |
US5720293A (en) * | 1991-01-29 | 1998-02-24 | Baxter International Inc. | Diagnostic catheter with memory |
US5553622A (en) * | 1991-01-29 | 1996-09-10 | Mckown; Russell C. | System and method for controlling the temperature of a catheter-mounted heater |
US6387052B1 (en) | 1991-01-29 | 2002-05-14 | Edwards Lifesciences Corporation | Thermodilution catheter having a safe, flexible heating element |
WO1992015239A1 (en) * | 1991-02-04 | 1992-09-17 | Kensey Nash Corporation | Apparatus and method for determining viscosity of the blood of a living being |
WO1992015241A1 (en) * | 1991-03-04 | 1992-09-17 | Kensey Nash Corporation | Apparatus and method for determining in vivo the instantaneous velocity of the flow of blood in a living being |
WO1992015240A1 (en) * | 1991-03-04 | 1992-09-17 | Kensey Nash Corporation | Determining the instantaneous shear stress on the wall of a blood vessel |
US5217019A (en) * | 1991-12-27 | 1993-06-08 | Abbott Laboratories | Apparatus and method for continuously monitoring cardiac output |
US5261411A (en) * | 1991-12-27 | 1993-11-16 | Abbott Laboratories | Thermal drift correction while continuously monitoring cardiac output |
US5305760A (en) * | 1992-02-07 | 1994-04-26 | Interflo Medical Inc. | Method for rejecting electrical interference from physiological measurements |
US5357967A (en) * | 1993-06-04 | 1994-10-25 | Baxter International Inc. | Method and apparatus for measuring flow using frequency-dispersive techniques |
US5363856A (en) * | 1993-08-13 | 1994-11-15 | Abbott Laboratories | Correcting thermal drift in cardiac output determination |
US5797398A (en) * | 1993-08-13 | 1998-08-25 | Thermal Technologies, Inc. | Method and apparatus for measuring continuous blood flow at low power |
US5423323A (en) * | 1993-08-30 | 1995-06-13 | Rocky Mountain Research, Inc. | System for calculating compliance and cardiac hemodynamic parameters |
CN1068429C (zh) * | 1993-12-23 | 2001-07-11 | 安德雷斯和霍瑟·弗罗泰克有限公司 | 卡式超声容积流量计 |
US5636638A (en) * | 1994-06-29 | 1997-06-10 | Baxter International Inc. | Electrical power amplifier for continuous cardiac output monitoring |
US5650571A (en) * | 1995-03-13 | 1997-07-22 | Freud; Paul J. | Low power signal processing and measurement apparatus |
US5687733A (en) * | 1995-10-26 | 1997-11-18 | Baxter International Inc. | System and method for estimating cardiac output |
US5620002A (en) * | 1995-12-22 | 1997-04-15 | Abbott Critical Care Systems | Method for correcting thermal drift in cardiac output determination |
US5954659A (en) * | 1996-02-16 | 1999-09-21 | Curley; Michael G. | Method and apparatus for invasively measuring cardiac output by detecting temperature differences of blood heated by radiation |
US5886576A (en) * | 1996-10-16 | 1999-03-23 | Baxter International Inc. | Electrical power amplifier for continuous cardiac output monitoring |
US5928155A (en) * | 1997-01-24 | 1999-07-27 | Cardiox Corporation | Cardiac output measurement with metabolizable analyte containing fluid |
US6169965B1 (en) * | 1997-12-31 | 2001-01-02 | Honeywell International Inc. | Fluid property and flow sensing via a common frequency generator and FFT |
US6223593B1 (en) | 1997-12-31 | 2001-05-01 | Honeywell International Inc. | Self-oscillating fluid sensor |
US6299583B1 (en) | 1998-03-17 | 2001-10-09 | Cardiox Corporation | Monitoring total circulating blood volume and cardiac output |
US6045512A (en) * | 1998-06-09 | 2000-04-04 | Baxter International Inc. | System and method for continuous estimation and display of cardiac ejection fraction and end diastolic volume |
US6393894B1 (en) | 1999-07-27 | 2002-05-28 | Honeywell International Inc. | Gas sensor with phased heaters for increased sensitivity |
ATE227022T1 (de) * | 1999-09-27 | 2002-11-15 | Mettler Toledo Gmbh | Verfahren und vorrichtung für die thermische analyse eines materials |
ATE265177T1 (de) | 1999-10-28 | 2004-05-15 | Pulsion Medical Sys Ag | Vorrichtung, computersystem und computerprogramm zur bestimmung eines kardiovasculären parameters |
US6371923B1 (en) | 1999-12-07 | 2002-04-16 | Edwards Lifesciences Corporation | Time-domain system and method for relaxation measurement and estimation of indicator dilution for continuous estimation and display of cardiac ejection fraction and end diastolic volume |
US6502459B1 (en) | 2000-09-01 | 2003-01-07 | Honeywell International Inc. | Microsensor for measuring velocity and angular direction of an incoming air stream |
ES2220607T3 (es) | 2001-03-01 | 2004-12-16 | Pulsion Medical Systems Ag | Aparato, programa informatico y conjunto de cateter venoso central para control hemodinamico. |
US6746408B2 (en) * | 2001-05-29 | 2004-06-08 | Transonic Systems Inc. | Method of blood flow measurement in arterio-venous hemodialysis shunts by indicator dilution |
US7131950B1 (en) * | 2002-09-23 | 2006-11-07 | E.P. Limited | System and method for noise reduction in thermodilution for cardiac measurement |
US20040224422A1 (en) * | 2002-09-27 | 2004-11-11 | Ulrich Bonne | Phased micro analyzer III, IIIA |
US20040223882A1 (en) * | 2002-09-27 | 2004-11-11 | Ulrich Bonne | Micro-plasma sensor system |
US20050063865A1 (en) * | 2002-09-27 | 2005-03-24 | Ulrich Bonne | Phased VII micro fluid analyzer having a modular structure |
WO2004029602A2 (en) * | 2002-09-27 | 2004-04-08 | Honeywell International Inc. | Phased micro analyser |
US7104112B2 (en) * | 2002-09-27 | 2006-09-12 | Honeywell International Inc. | Phased micro analyzer IV |
US7530257B2 (en) | 2002-09-27 | 2009-05-12 | Honeywell International Inc. | Phased micro analyzer VIII |
US7000452B2 (en) * | 2002-09-27 | 2006-02-21 | Honeywell International Inc. | Phased micro fluid analyzer |
US7494326B2 (en) * | 2003-12-31 | 2009-02-24 | Honeywell International Inc. | Micro ion pump |
US7452333B2 (en) * | 2003-12-05 | 2008-11-18 | Edwards Lifesciences Corporation | Arterial pressure-based, automatic determination of a cardiovascular parameter |
US7422562B2 (en) * | 2003-12-05 | 2008-09-09 | Edwards Lifesciences | Real-time measurement of ventricular stroke volume variations by continuous arterial pulse contour analysis |
US7220230B2 (en) * | 2003-12-05 | 2007-05-22 | Edwards Lifesciences Corporation | Pressure-based system and method for determining cardiac stroke volume |
US9029028B2 (en) | 2003-12-29 | 2015-05-12 | Honeywell International Inc. | Hydrogen and electrical power generator |
US20050142035A1 (en) * | 2003-12-31 | 2005-06-30 | Ulrich Bonne | Micro-discharge sensor system |
US7632235B1 (en) | 2004-11-22 | 2009-12-15 | Pacesetter, Inc. | System and method for measuring cardiac output via thermal dilution using an implantable medical device with an external ultrasound power delivery system |
US7651466B2 (en) | 2005-04-13 | 2010-01-26 | Edwards Lifesciences Corporation | Pulse contour method and apparatus for continuous assessment of a cardiovascular parameter |
US7502109B2 (en) * | 2005-05-17 | 2009-03-10 | Honeywell International Inc. | Optical micro-spectrometer |
US7578167B2 (en) * | 2005-05-17 | 2009-08-25 | Honeywell International Inc. | Three-wafer channel structure for a fluid analyzer |
EP1904814A1 (de) * | 2005-07-21 | 2008-04-02 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Messeinrichtung zum medienunabhängigen messen einer strömungsgeschwindigkeit eines mediums |
US7424366B2 (en) * | 2005-08-27 | 2008-09-09 | Schlumberger Technology Corporation | Time-of-flight stochastic correlation measurements |
US20070167861A1 (en) * | 2005-11-29 | 2007-07-19 | Lopez George A | Cardiac output measurement devices and methods |
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 |
WO2008144525A1 (en) * | 2007-05-16 | 2008-11-27 | Massachusetts Institute Of Technology | System and method for prediction and detection of circulatory shock |
WO2008144490A1 (en) * | 2007-05-16 | 2008-11-27 | Massachusetts Instutute Of Technology | Systems and methods for model-based estimation of cardiac ejection fraction, cardiac contractility, and ventricular end-diastolic volume |
US8282564B2 (en) | 2007-05-16 | 2012-10-09 | Massachusetts Institute Of Technology | Systems and methods for model-based estimation of cardiac output and total peripheral resistance |
EP2062528B1 (en) * | 2007-11-23 | 2012-03-14 | Pulsion Medical Systems AG | Central venous sensor assembly for measuring physiological data for cardiac output determination and method of determining cardiac output |
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 |
DE102008026708B4 (de) | 2008-06-04 | 2014-01-23 | Iprm Intellectual Property Rights Management Ag | Vorrichtung zur Bestimmung des Blutvolumens und/oder Blutvolumenstroms und Verfahren zum Betreiben derselben |
CN102834047A (zh) | 2010-01-29 | 2012-12-19 | 爱德华兹生命科学公司 | 在心血管参数确定中不规则心动周期影响的排除 |
WO2015006295A1 (en) | 2013-07-08 | 2015-01-15 | Edwards Lifesciences Corporation | Determination of a hemodynamic parameter |
EP3346913A4 (en) | 2015-09-25 | 2019-01-16 | C.R. Bard Inc. | CATHETER ASSEMBLY COMPRISING MONITORING CAPABILITIES |
DE102017128576A1 (de) | 2017-12-01 | 2019-06-06 | Ilmsens Gmbh | Vorrichtung zur Vitalitätsüberwachung |
WO2021142142A1 (en) | 2020-01-07 | 2021-07-15 | Bard Access Systems, Inc. | Diagnostic systems and methods including temperature-sensing vascular devices |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3214973A (en) * | 1963-02-14 | 1965-11-02 | Robert A Lester | Acoustic flow meter for measuring very slow fluid flow |
US3543578A (en) * | 1968-02-26 | 1970-12-01 | Universal Oil Prod Co | Flow metering system |
US3842670A (en) * | 1969-03-28 | 1974-10-22 | Nat Res Dev | Improvements in measuring the velocity of gases |
US3820393A (en) * | 1970-12-09 | 1974-06-28 | Gunigle R Mc | Signal velocity measurement system |
FR2291731A1 (fr) * | 1974-11-20 | 1976-06-18 | Obermajer Wladimir | Procede et dispositif de mesure de la capacite cardiaque |
US4244026A (en) * | 1978-11-06 | 1981-01-06 | General Electric Company | Velocity measuring correlation sonar |
-
1983
- 1983-04-21 US US06/487,334 patent/US4507974A/en not_active Expired - Lifetime
-
1984
- 1984-04-12 EP EP84104121A patent/EP0126931B1/en not_active Expired
- 1984-04-12 DE DE8484104121T patent/DE3483626D1/de not_active Expired - Lifetime
- 1984-04-16 CA CA000452128A patent/CA1210509A/en not_active Expired
- 1984-04-20 JP JP59078722A patent/JPS6034430A/ja active Granted
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63216536A (ja) * | 1987-03-05 | 1988-09-08 | テルモ株式会社 | 心拍出量の測定装置 |
JPS63216539A (ja) * | 1987-03-05 | 1988-09-08 | テルモ株式会社 | 心拍出量の相対的変化の連続測定装置 |
JPS63216537A (ja) * | 1987-03-05 | 1988-09-08 | テルモ株式会社 | 心拍出量の測定装置 |
JPS63216540A (ja) * | 1987-03-05 | 1988-09-08 | テルモ株式会社 | 心拍出量の測定装置 |
JPH0467854B2 (ja) * | 1987-03-05 | 1992-10-29 | Terumo Corp | |
JPH0467856B2 (ja) * | 1987-03-05 | 1992-10-29 | Terumo Corp | |
JPH0467853B2 (ja) * | 1987-03-05 | 1992-10-29 | Terumo Corp | |
JPH0467855B2 (ja) * | 1987-03-05 | 1992-10-29 | Terumo Corp | |
JPS63286165A (ja) * | 1987-05-01 | 1988-11-22 | プロダクト・イノウベイション・ホールディングズ・リミテッド | 医学注入装置 |
US7688098B2 (en) | 2007-03-08 | 2010-03-30 | Nec Corporation | Power supply noise measuring circuit and power supply noise measuring method |
Also Published As
Publication number | Publication date |
---|---|
EP0126931A3 (en) | 1986-01-08 |
DE3483626D1 (de) | 1991-01-03 |
EP0126931A2 (en) | 1984-12-05 |
US4507974A (en) | 1985-04-02 |
EP0126931B1 (en) | 1990-11-22 |
CA1210509A (en) | 1986-08-26 |
JPH0527051B2 (ja) | 1993-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPS6034430A (ja) | フロー測定方法および装置 | |
EP1125548B1 (en) | Triggered flow measurement | |
EP0525124B1 (en) | Method and apparatus for continuously measuring volumetric flow | |
JP3701285B2 (ja) | 心拍出量を評価するためのシステム | |
EP1076511B1 (en) | System for measuring flow in narrow channels comprising a temperature and a pressure sensor | |
Yelderman | Continuous measurement of cardiac output with the use of stochastic system identification techniques | |
US8784327B2 (en) | Method and system for obtaining dimension related information for a flow channel | |
KR100416894B1 (ko) | 심장 박출 계수 및 말기 확장기 용적의 연속적인 추정 및디스플레이를 위한 시스템 및 방법 | |
RU2179408C2 (ru) | Способ индивидуального определения в организме функции эластичности и непрерывного определения системного кровотока живого существа | |
US7890269B2 (en) | Method and device for acoustic measurement of the specific density of a gaseous or liquid medium | |
JPH0296636A (ja) | 粒径測定装置 | |
EP0701413B1 (en) | Method and apparatus for measuring flow using frequency-dispersive techniques | |
US4230126A (en) | Apparatus and method for measuring extravascular lung water | |
US20040088123A1 (en) | Method for real-time monitoring of cardiac output and blood flow in arteries and apparatus for implementing the same | |
Battocletti et al. | Analysis of a nuclear magnetic resonance blood flowmeter for pulsatile flow | |
Poots et al. | Comparison of CW Doppler ultrasound spectra with the spectra derived from a flow visualization model | |
Mellander et al. | Venous blood flow recorded with an isothermal flowmeter | |
Berbich et al. | Non-linear analysis of the arterial pulsatile flow: assessment of a model allowing a non-invasive ultrasonic functional exploration | |
EP0015294B1 (en) | Apparatus for measuring extravascular lung water | |
Zarinabad et al. | Modelling parameter role on accuracy of cardiac perfusion quantification | |
JP2002022508A (ja) | サンプリング方法および異常脈動検出方法 | |
Barnett et al. | Measurement of instantaneous blood flow | |
Melbin | A General Approach for Computation and Interpretation of Dilution Curves | |
De | Theoretical investigation of estimation of steady and pulsatile blood flow and blood vessel cross section by CW NMR excitation | |
Baker | Cardiac Blood Flow Measurement: A Component Of The Comprehensive Cardiac Examination |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EXPY | Cancellation because of completion of term |