EP1022984A1 - Bestimmung der leberfunktion anhand einer plasmaverschwinderate - Google Patents

Bestimmung der leberfunktion anhand einer plasmaverschwinderate

Info

Publication number
EP1022984A1
EP1022984A1 EP99942720A EP99942720A EP1022984A1 EP 1022984 A1 EP1022984 A1 EP 1022984A1 EP 99942720 A EP99942720 A EP 99942720A EP 99942720 A EP99942720 A EP 99942720A EP 1022984 A1 EP1022984 A1 EP 1022984A1
Authority
EP
European Patent Office
Prior art keywords
dye
concentration
pdr
measurement
rez
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
Application number
EP99942720A
Other languages
English (en)
French (fr)
Other versions
EP1022984B1 (de
Inventor
Andreas Dr. Med. Hoeft
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Publication of EP1022984A1 publication Critical patent/EP1022984A1/de
Application granted granted Critical
Publication of EP1022984B1 publication Critical patent/EP1022984B1/de
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

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
    • A61B5/0275—Measuring blood flow using tracers, e.g. dye dilution
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/42—Detecting, measuring or recording for evaluating the gastrointestinal, the endocrine or the exocrine systems
    • A61B5/4222—Evaluating particular parts, e.g. particular organs
    • A61B5/4244—Evaluating particular parts, e.g. particular organs liver

Definitions

  • the invention relates to a method for determining liver function using a plasma shrinkage rate according to the preamble of claim 1 and a device for determining liver function using a plasma shrinkage rate according to the preamble of claim 7.
  • liver function is an important parameter in intensive care medicine and crucial for the prognosis of seriously ill patients.
  • liver function is being routinely assessed in intensive care medicine using various laboratory sizes that characterize the synthetic performance and the elimination performance of the liver.
  • the disadvantage of these laboratory sizes is that if the liver function fails, they only change to abnormal values with a considerable latency over time, so that a liver function disorder becomes evident only after days.
  • the indocyanine green is usually injected intravenously as a bolus and, following the bolus injection, at least 2, usually several, blood samples are taken at intervals of several minutes over a period of at least 15 minutes.
  • the elimination time constant can be calculated from the drop in the dye concentration in the blood samples. This procedure is described in but rarely used in the clinic, because it is still too time-consuming due to the laboratory analyzes.
  • the invention has for its object to provide a method and a device with which the liver function can be determined non-invasively and the measurement result is available more quickly.
  • a decrease in a dye concentration in the blood is measured, which is caused by the breakdown of the dye by the liver.
  • a suitable dye such as. B. indocyanine green
  • injection of a suitable dye such as. B. indocyanine green
  • the second maximum arises from recirculation, i.e. before the concentration wave decays on a first pass, a second pass already takes place.
  • the mean circulation transit time previously given mtt. can be calculated from a circulatory transport function g (t) that describes the transport behavior of the circulatory system.
  • the mean cycle transit time mtt circ is then according to the formula
  • the circulation transport function g (t) can be calculated from the measured dye concentration with the aid of an iterative, non-linear adaptation method.
  • a model function ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
  • g (t) a m g m (t) + a m + ⁇ g m + 1 (t) + ... + a n g n (t)
  • c (t) form the concentration-time profile of the indicator dye
  • c D0 us (t) the first portion of the dye concentration-time profile directly passing the measurement site
  • c rez (t) a recirculating portion of the dye Concentration-time course
  • k the elimination fraction of the dye eliminated by the liver.
  • optical measurement of the dye concentration resulting in the blood circulation can alternatively be carried out by fiber-optic measurement in a central vessel or non-invasively by means of light transmission or reflection measurement of light irradiated at suitable body locations, in particular on the fingers, earlobes, bridge of the nose, cheek mucosa or forehead.
  • 1 is a graphical representation of the concentration-time profile of an indicator dye
  • 2 is a semi-logarithmic representation of the concentration-time curve according to FIG. 1
  • Fig. 3 is a graphical representation of the circulatory transport function
  • Fig. 4 is a block diagram of a device for measuring and determining the plasma shrinkage rate.
  • 1 shows a graphical representation of the concentration-time curve of an indicator substance after a bolus injection.
  • sensors are arranged to record the dye concentration-time curve, which are connected to an evaluation circuit and a computer.
  • An example of a device according to the invention is shown in FIG. 4.
  • the use of a reference dye is cheap but not absolutely necessary.
  • the red blood pigment (hemoglobin) can be used as a reference dye.
  • one light transmitter 10 operates on a wavelength of 800 nanometers and the other light transmitter 12 operates on a wavelength of 900 nanometers.
  • the light receivers 14 and 16 can be designed so that they preferentially filter out these wavelengths.
  • the absorption behavior of the dyes can be measured both by light transmission or by light reflection measurement. Suitable places on the body are the fingers, the earlobes, the bridge of the nose, the cheek mucosa or the forehead.
  • the signals of the light intensities received by the light receivers 14 and 16 are evaluated by the evaluation circuit 18 and fed to a computer 20.
  • This computer 20 first preprocesses the pulsatile components of the light intensities I indpuls (t) and I r ⁇ £ puls (t), for example by forming the quotient:
  • the relative dye concentration as a function of time can be determined from these preprocessed signals, also by forming the quotient:
  • a non-linear adaptation algorithm from c (t) is used to determine a circulatory transport function g (t), which describes the transport behavior of the circulatory system, and a parameter k, which determines a fractional re-entry rate of the Represents dye after each circulation passage.
  • the following is based on a model in which a transport route of the injected indicator dye ICG as a bolus injection first leads through the lungs. Behind the lungs, the transport route is divided into a first, fast compartment, a second, slow compartment and an elimination of the dye by the liver. The first and second compartments then lead back to the lungs.
  • the concentration of the indicator dye at the exit of the lungs follows the function:
  • c (t) form the concentration-time profile of the indicator dye
  • k the fractional re-entry rate of the dye after partial elimi- liver liver.
  • the factor k is therefore always less than 1.
  • the compartments of the circulatory model are used for g (t).
  • Preliminary investigations showed that for dye dilution curves, which were measured on patients, one or two compartments have to be used for g (t), depending on the accuracy that is sought. This results in the general model function
  • g (t) a m g m (t) + a m + 1 g m + 1 (t) + ... + a n g n (t)
  • g (t) a g ⁇ t) + a2g2 ( fc ) -
  • Distribution functions can be optimized according to the method of the smallest quadratic deviation.
  • a transport function results as it was calculated as an example for the data contained in FIG. 1 for two compartments and is shown in FIG. 3.
  • the resulting transport function g (t) is composed of the transport function gl (t) for the first compartment and g2 (t) for the second compartment.
  • the mean cycle transit time (mtt.), which is the first moment of g (t), is also calculated from the transport function g (t): mt
  • the plasma shrinkage rate (PDR) of the ICG dye then results as:
  • the time constant K ⁇ 1xm of the elimination phase is determined in the linear part of the curve.
  • the PDR is:

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physiology (AREA)
  • Hematology (AREA)
  • Cardiology (AREA)
  • Gastroenterology & Hepatology (AREA)
  • Endocrinology (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

Die Leberfunktion wird nach Injektion eines Indikatorfarbstoffes in den Blutkreislauf, optische Messung der im Blutkreislauf resultierenden Farbstoffkonzentration und Berechnung einer Plasmaverschwindungrate aus dem Farbstoff-Konzentrations-Zeit-Verlauf des Indikatorfarbstoffes bestimmt. Dabei wird aus der Messung der Farbstoffkonzentration c(t) eine mittlere Kreislauftransitzeit mttcirc berechnet. Außerdem wird ein Parameter k ermittelt, der eine fraktionelle Wiedereintrittsrate des Indikatorfarbstoffes nach jeder Kreislaufpassage darstellt. Die Plasmaverschwindungrate PDR ergibt sich dann nach der Formel: PDR = (1-k)/(k.mttcirc).

Description

BESTIMMUNGDES LEBERFUNKTIONANHAND EINER PLASMAVERSCHWINDERATE
Die Erfindung betrifft ein Verfahren zur Bestimmung der Leberfunktion anhand einer Plasmaverschwinderate nach dem Oberbegriff des Anspruchs 1 und eine Vorrichtung zur Bestimmung der Leberfunktion anhand einer Plasmaverschwinderate nach dem Oberbegriff des Anspruchs 7.
Die Leberfunktion ist ein wichtiger Parameter in der Intensivmedizin und für die Prognose von schwerstkranken Patienten entscheidend. Zur Zeit wird die Leberfunktion anhand von verschiedenen Laborgrößen, die die synthetische Leistung und die Eliminationsleistung der Leber charakterisieren, in der Routine in der Intensivmedizin erfaßt. Der Nachteil dieser Laborgrößen ist jedoch, daß sie bei Ausfall der Leberfunktion erst mit einer erheblichen zeitlichen Latenz auf krankhafte Werte verändert sind, so daß eine Leberfunktionsstörung erst nach Tagen evident wird.
Eine Möglichkeit die Leberfunktion, zumindest bezüglich ihrer Eliminationsleistung sofort zu evaluieren, besteht in der Applikation von Indikatorstoffen, die durch die Leber eliminiert werden und der Bestimmung der Eliminationszeitkonstante dieser Indikatoren. Ein hierfür gängiger Indikator ist der Farbstoff Indocyaningrün. Üblicherweise wird das In- docyaningrün als Bolus intravenös injiziert und im Anschluß an die Bolusinjektion wird in Abständen von mehreren Minuten über eine Zeitspanne von mindestens 15 Minuten mindestens 2, meistens mehrere Blutproben entnommen. Aus dem Abfall der Farbstoffkonzentration in den Blutproben kann die Eliminati- onszeitkonstante berechnet werden. Dieses Verfahren wird in der Klinik jedoch selten angewendet, da es wegen der Laboranalysen immer noch zu zeitaufwendig ist.
Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren und eine Vorrichtung zu schaffen mit denen die Leberfunktion nichtinvasiv bestimmt werden kann und das Meßergebnis schneller zur Verfügung steht.
Diese Aufgabe wird bei einem Verfahren nach dem Oberbegriff des Anspruchs 1 durch die Merkmale dieses Anspruchs gelöst und bei einer Vorrichtung nach dem Oberbegriff des Anspruchs 7 durch die Merkmale jenes Anspruchs gelöst.
Auch bei dem erfindungsgemäßen Verfahren und der Vorrichtung wird eine Abnahme einer Farbstoffkonzentration im Blut gemessen, die durch Abbau des Farbstoffes durch die Leber zustandekommt. Nach Injektion eines geeigneten Farbstoffes, wie z. B. Indocyaningrün, in die Blutbahn ergibt sich bei normaler Leberfunktion an einer Meßstelle des Körpers ein charakteristischer Farbstoff-Konzentrations-Zeit-Verlauf . Hier zeigt sich zunächst ein erstes Maximum der Farbstoff- Konzentration und nach einem vorübergehenden Abklingen ein zweites Maximum der Farbstoff-Konzentration. Das zweite Maximum entsteht durch Rezirkulation, d.h. vor Abklingen der Konzentrationswelle bei einem ersten Durchlauf erfolgt bereits ein zweiter Durchlauf .
Aus der Messung des zeitlichen Verlaufs der Farbstoffkonzentration c(t) wird eine mittlere Kreislauftransitzeit mttoirc berechnet und außerdem ein Parameter k ermittelt, der eine fraktionelle Wiedereintrittsrate des Indikatorfarbstoffes nach jeder Kreislaufpassage darstellt. Mit diesen Werten kann dann die Plasmaverschwinderate (PDR) nach der Formel PDR = (l-k)/(k • mttcirc) berechnet werden. Das Ergebnis liegt bereits nach wenigen Rezirkulationszyklen vor.
Die zuvor angegebene mittlere Kreislauftransitzeit mtt . kann aus einer Kreislauftransportfunktion g(t) berechnet werden, die das Transportverhalten des Kreislaufs beschreibt. Die mittlere Kreislauftransitzeit mttcirc wird dann nach der Formel
mtt;„ = erhalten.
Die Berechnung der Kreislauftransportfunktion g(t) kann aus der gemessenen Farbstoffkonzentration mit Hilfe eines iterativen, nichtlinearen Anpassungsverfahrens durchgeführt werden. Hier wird unter Vorgabe einer Modellfunktion
g(t) = amgm(t) + am+ιgm+1(t) + ... + angn(t)
mit m=\ *m = 1,
bei der die einzelnen Kompartimente amgm durch linksschiefe Verteilungsfunktionen beschrieben sind, wiederholt eine rekursive Faltung gemäß den Gleichungen
c(t) = cbolus(t) + crez(t)
und crez ( t) = k * j" g ( t-u) * c (u) du
durchgeführt. Dabei werden die Parameter k, am sowie die Parameter der Verteilungsfunktionen nach dem Verfahren der kleinsten quadratischen Abweichung optimiert, wobei wenigstens ein Kompartiment aχg]_(t) angesetzt wird. In der Formel bilden c(t) den Konzentrations-Zeit-Verlauf des Indikatorfarbstoffes, cD0 us(t) der den Meßort direkt passierende erste Anteil des Farbstoff-Konzentrations-Zeit-Verlaufs , crez(t) einen rezirkulierender Anteil des Farbstoff- Konzentrations-Zeit-Verlaufs und k die Eliminationsfraktion des durch die Leber eliminierten Farbstoffes .
Für eine größere Genauigkeit können auch zwei Kompartimente (aι_g]_(t)) und (a2 2( )) angesetzt werden.
Die optische Messung der im Blutkreislauf resultierenden Farbstoffkonzentration kann alternativ durch fiberoptische Messung in einem zentralen Gefäß vorgenommen werden oder nicht-invasiv mittels Licht-Transmissions- oder Reflexionsmessung eingestrahlten Lichtes an geeigneten Körperorten, insbesondere am Finger, Ohrläppchen, Nasenrücken, Wangenschleimhaut oder Stirn.
Nachfolgend wird die Erfindung anhand der Zeichnung erläutert. In dieser zeigen:
Fig. 1 eine grafische Darstellung des Konzentrations-Zeit- Verlaufs eines Indikatorfarbstoffes, Fig. 2 eine halblogarith ische Darstellung des Konzentrations-Zeit-Verlaufs nach Fig. 1,
Fig. 3 eine grafische Darstellung der Kreislauftransportfunktion und
Fig. 4 ein Blockschaltbild einer Vorrichtung zur Messung und Bestimmung der Plasmaverschwinderate.
Fig. 1 zeigt eine grafische Darstellung des Konzentrations- Zeit-Verlaufs eines Indikator arbstoffes nach einer Bolusin- jektion.
Das nach dem ersten Maximum folgende zweite Maximum entsteht durch Rezirkulation, d.h. vor Abklingen der Konzentrationswelle bei einem ersten Durchlauf erfolgt bereits ein zweiter Durchlauf .
An einer Meßstelle des Körpers sind zur Erfassung des Farbstoff-Konzentrations-Zeit-Verlaufs Sensoren angeordnet, die mit einer Auswerteschaltung und einem Rechner verbunden sind. Ein Beispiel für eine erfindungsgemäße Vorrichtung zeigt Fig. 4. Dort sind zwei Lichtsender 10 und 12 und zwei Lichtempfänger 14 und 16 vorgesehen, die über eine Auswerteschaltung 18 mit einem Rechner 20 verbunden sind. Es sei davon ausgegangen, daß für die Messung der Leberfunktion als Indikatorfarbstoff Indocyaningrün verwendet wird.
Die Verwendung eines Referenzfarbstoffes ist günstig aber nicht zwingend notwendig. Als Referenzfarbstoff kann der rote Blutfarbstoff (Hämoglobin) verwendet werden. In diesem Fall arbeitet der eine Lichtsender 10 auf einer Wellenlänge von 800 Nanometern und der andere Lichtsender 12 auf einer Wellenlänge von 900 Nanometern. Die Lichtempfänger 14 und 16 können so ausgeführt sein, daß sie auch bevorzugt diese Wellenlängen herausfiltern. Alternativ kann aber auch zwischen zwei Lichtsendern 10 und 12 umgeschaltet werden und man kommt in diesem Fall mit nur einem Lichtempfänger 14 aus, der aber für beide Wellenlängen ausgelegt sein muß.
Die Messung des Absorptionsverhaltens der Farbstoffe kann sowohl durch Licht-Transmissions- oder auch durch Licht- Reflexionsmessung erfolgen. Geeignete Körperorte sind die Finger, die Ohrläppchen, der Nasenrücken, die Wangenschleimhaut oder die Stirn.
Die von den Lichtempfängern 14 und 16 empfangenen Signale der Lichtintensitäten werden von der Auswerteschaltung 18 ausgewertet und einem Rechner 20 zugeführt. In diesem Rechner 20 erfolgt zunächst eine Vorverarbeitung der pulsatilen Anteile der Lichtintensitäten Iindpuls(t) und Irβ£puls(t), beispielsweise durch Quotientenbildung:
Iin p ls ( ) = Iina ( t-dt ) / Iind ( t+dt ) ι«fpui.( ) = ι„.(t-dt) / ιre£(t+dt)
Aus diesen vorverarbeiteten Signalen kann in einem weiteren Schritt bereits die relative Farbstoff-Konzentration als Funktion der Zeit bestimmt werden, und zwar ebenfalls durch Quotientenbildung :
c(t) ind puls ref puls
Es bedarf hier keines Kalibrationsfaktors, da für die Bestimmung der Leberfunktion nur ein lineares Signal benötigt wird. Das Meßprinzip beruht nicht auf der Bestimmung absoluter Konzentrationswerte, sondern nur auf der Bestimmung von Zeitkonstanten .
Bei der weiteren Auswertung der Farbstoff-Konzentration als Funktion der Zeit wird anschließend mittels nichtlinearer Anpassungsalgorithmen aus c(t) eine Kreislauftransportf nk- tion g(t), die das Transportverhalten des Kreislaufs beschreibt, sowie ein Parameter k ermittelt, der eine fraktionelle Wiedereintrittsrate des Farbstoffes nach jeder Kreislaufpassage darstellt.
Im folgenden legt man ein Modell zugrunde, bei dem ein Transportweg des injizierten Indikatorfarbstoffes ICG als Bolusinjektion zuerst durch die Lunge führt. Hinter der Lunge teilt sich der Transportweg in ein erstes, schnelles Kompartiment, ein zweites langsames Kompartiment und in eine Elimination des Farbstoffes durch die Leber. Das erste und zweite Kompartiment führen anschließend wieder zur Lunge zurück. Die Konzentration des Indikatorfarbstoffes am Ausgang der Lunge folgt der Funktion:
c(t) = cboιus(t) + crez(t)
-rez (t) = k * J g(t-u) * c(u) du
Hierbei bilden c(t) den Konzentrations-Zeit-Verlauf des Indikatorfarbstoffes, Ck0τ_us(t) der den Meßort direkt passierende erste Anteil des Farbstoff-Konzentrations-Zeit- Verlaufs, crez(t) einen rezirkulierender Anteil des Farbstoff-Konzentrations-Zeit-Verlaufs und k die fraktionelle Wiedereintrittsrate des Farbstoffes nach teilweiser Elimi- nierung durch die Leber. Der Faktor k ist daher stets kleiner als 1.
Der Transportprozeß wird durch ein Faltungintegral beschrieben. Rezirkulation in dieser Darstellungsweise bedeutet, daß das Ergebnis der Faltung crez(t) gleichzeitig auch die Eingangsfunktion c(t) beeinflußt. Die Rezirkulation des Indikators führt somit zu einem Zusammenhang, der im Prinzip folgendermaßen beschrieben ist:
c(t) = f [k, g(t) , c(t)]
Für g(t) werden die Kompartimente des Kreislaufmodells angesetzt. Voruntersuchungen ergaben, daß bei Farbstoffverdün- nungskurven, die an Patienten gemessen worden sind, für g(t) ein bis zwei Kompartimente anzusetzen sind, je nachdem, welche Genauigkeit angestrebt wird. Damit ergibt sich für die allgemeine Modellfunktion
g(t) = amgm(t) + am+1gm+1(t) + ... + angn(t)
bei der die einzelnen Kompartimente amgm durch linksschiefe
Verteilungsfunktionen beschrieben sind. Für zwei Kompartimente gilt der spezielle Ausdruck:
g(t) = a gη t) + a2g2(fc)-
Die Berechnung der Kreislauftransportfunktionen aus gemessenen Farbstoffkurven wird durch den Rechner mit Hilfe eines iterativen, nicht-linearen Anpassungsverf hrens durchgeführt werden, bei dem unter Vorgabe der Modellfunktion
wiederholt eine rekursive Faltung gemäß den Gleichungen
c(t) = cbolus(t) + crez(t)
und
-rez (t) = k J g(t- u) c (u) du
durchgeführt wird, wobei k, aι_, a2 und die Parameter der
Verteilungsfunktionen nach dem Verfahren der kleinsten quadratischen Abweichung optimiert werden.
Nach Durchführung dieser Rechenschritte ergibt sich eine Transportfunktion, wie sie als Beispiel für die in Fig. 1 enthaltenen Daten für zwei Kompartimente errechnet wurde und in Fig. 3 dargestellt ist. Die resultierende Transportfunktion g(t) setzt sich aus der Transportfunktion gl(t) für das erste Kompartiment und g2 (t) für das zweite Kompartiment zu- sammem.
Aus der Transportfunktion g(t) wird des weiteren die mittlere Kreislauftransitzeit berechnet (mtt . ) , die das erste Mo- ment von g(t) ist: m t
Die Plasmaverschwinderate (PDR) des Farbstoffes ICG ergibt sich sodann als:
PDR = (l-k)/(k • mtteιre)
Alternativ kann auch so vorgegangen werden, daß bei der halblogarithmischen Darstellung der Konzentrations-Zeit- Kurve nach Fig. 2 in dem linearen Teil der Kurve die Zeitkonstante Kβlxm der Eliminationsphase bestimmt wird. In diesem Falle ist die PDR:
PDR = -Kelιm /x 100 x 60 (für die Konvertierung in Einhei- ten/%/min. )

Claims

P a t e n t a n s p r ü c h e
1. Verfahren zur Bestimmung der Leberfunktion anhand einer Plasmaverschwindungrate (PDR) nach Injektion eines Indikatorfarbstoffes in den Blutkreislauf, optische Messung der im Blutkreislauf resultierenden Farbstoffkonzentration und Berechnung der Plasmaverschwinderate aus dem Farbstoff- Konzentrations-Zeitverlauf des Indikatorfarbstoffes, dadurch gekennzeichnet , daß aus der Messung der Farbstoffkonzentration (c(t)) eine mittlere Kreislauftransitzeit (mttoιrc.) berechnet wird, daß ein Parameter (k) ermittelt wird, der eine fraktionelle Wiedereintrittsrate des Indikatorfarbstoffes nach jeder Kreislaufpassage darstellt, und daß die Plasmaverschwindungrate (PDR) nach der Formel
PDR = (l-k)/(k • mttcirc) berechnet wird.
2. Verfahren nach Anspruch 1, dadurch gekennzeichnet , daß die mittlere Kreislauftransitzeit (mttcirc) aus einer Kreislauftransportfunktion (g(t)) berechnet wird, die das Transportverhalten des Kreislaufs beschreibt, nach der Formel
mttcirc
3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daß die Berechnung der Kreislauftransportfunktion (g(t)) aus der gemessenen Farbstoffkonzentration mit Hilfe eines iterativen, nichtlinearen Anpassungsverfahrens durchgeführt wird, bei dem unter Vorgabe einer Modellfunktion
g(t) = amgm(t) + am+1gm+1(t) + ... + angn(t)
bei der die einzelnen Kompartimente amgm durch linksschiefe Verteilungsfunktionen beschrieben sind, wiederholt eine rekursive Faltung gemäß den Gleichungen
c(t) = cboιus(t) + crez(t)
und
crez (t) = k * J g(t-u) * c(u) du
durchgeführt wird und die Parameter k, &m sowie die Parameter der Verteilungsfunktionen nach dem Verfahren der kleinsten quadratischen Abweichung optimiert werden, wobei wenigstens ein Kompartiment aι_g]_(t) angesetzt wird und c(t) den Konzentrations-Zeit-Verlauf des Indikatorfarbstoffes, cιDO]_us(t) der den Meßort direkt passierende erste Anteil des
Farbstoff-Konzentrations-Zeit-Verlaufs, crez(t) einen rezirkulierender Anteil des Farbstoff-Konzentrations-Zeit- Verlaufs und k die Eliminationsfraktion des durch die Leber eliminierten Farbstoffes bilden.
4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, daß zwei Kompartimente (a]_g^(t)) und (a2g2 (t) ) angesetzt werden.
5. Verf hren nach einem der Ansprüche 1 bis 4 , dadurch gekennzeichnet, daß die optische Messung der im Blutkreislauf resultierenden Farbstoffkonzentration durch fiberoptische Messung in einem zentralen Gefäß vorgenommen wird.
6. Verfahren nach einem oder mehreren der Ansprüche 1 bis 4, dadurch gekennzeichnet, daß die optische Messung der im Blutkreislauf resultierenden Farbstoffkonzentration nicht- invasiv mitteis Licht-Transmissions- oder Reflexionsmessung eingestrahlten Lichtes an geeigneten Körperorten, insbesondere am Finger, Ohrläppchen, Nasenrücken, Wangenschleimhaut oder Stirn gemessen wird.
7. Vorrichtung zur Bestimmung der Leberfunktion anhand einer Plasmaverschwindungrate (PDR) mit einem optischen Meßsensor (10) und einem damit verbundenen Rechner (12), wobei mittels des Meßsensor (12) nach Injektion eines Indikatorfarbstoffes in den Blutkreislauf die im Blutkreislauf resultierende Farbstoffkonzentration gemessen wird und mittels des Rechners (12) die Plasmaverschwindungrate aus dem Farbstoff-Konzentrations-Zeitverlauf des Indikatorfarbstoffes bestimmbar ist, dadurch gekennzeichnet , daß der Rechner (12) so gesteuert ist, daß aus der Messung der Farbstoffkonzen- tration (c(t)) eine mittlere Kreislauftransitzeit (mttcirc) berechnet wird, daß ein Parameter (k) ermittelt wird, der eine fraktionelle Wiedereintrittsrate des Indikatorfarbstoffes nach jeder Kreislaufpassage darstellt, und daß die Plasmaverschwindungrate (PDR) nach der Formel
PDR = (l-k)/(k • mtt . ) berechnet wird.
8. Vorrichtung nach Anspruch 7, dadurch gekennzeichnet, daß der Rechner (12) ferner so gesteuert ist, daß die mittlere Kreislauftransitzeit (mttoxrc) aus einer Kreislauftransportfunktion (g(t) )berechnet wird, die das Transportverhalten des Kreislaufs beschreibt, nach der Formel
9 . Vorrichtung nach Anspruch 8, dadurch gekennzeichnet, daß der Rechner (12) so gesteuert ist, daß die Berechnung der Kreislauftransportfunktion (g(t)) aus der gemessenen Farbstoffkonzentration mit Hilfe eines iterativen, nichtlinearen Anpassungsverfahrens durchgeführt wird, bei dem unter Vorgabe einer Modellfunktion
g(t) = amgm{t) + am+ιgm+ι(t) + ... + angn(t)
bei der die einzelnen Kompartimente amgm durch linksschiefe Verteilungsfunktionen beschrieben sind, wiederholt eine rekursive Faltung gemäß den Gleichungen
c(t) = cbolus^) + crez (t)
und -rez ( t ] = k * f g ( t-u) * c (u ) du
durchgeführt wird und die Parameter k, am sowie die Parameter der Verteilungsfunktionen nach dem Verfahren der kleinsten quadratischen Abweichung optimiert werden, wobei wenigstens ein Kompartiment aι_g]_(t) angesetzt wird und c(t) den Konzentrations-Zeit-Verlauf des Indikatorfarbstoffes, cD0ιus(t) der den Meßort direkt passierende erste Anteil des
Farbstoff-Konzentrations-Zeit-Verlaufs, crez(t) einen rezirkulierender Anteil des Farbstoff-Konzentrations-Zeit- Verlaufs und k die Eliminationsfraktion des durch die Leber eliminierten Farbstoffes bilden.
10. Vorrichtung nach Anspruch 9, dadurch gekennzeichnet, daß der Rechner (12) so gesteuert ist, daß zwei Kompartimente (aj_gι_(t)) und (a2g2( )) angesetzt werden.
11. Vorrichtung nach einem der Ansprüche 7 bis 10, dadurch gekennzeichnet , daß- der optische Meßsensor (10) ein fiberoptischer Katheter ist, der zur Messung in einem zentralen Gefäß angeordnet wird.
12. Vorrichtung nach einem der Ansprüche 7 bis 10, dadurch gekennzeichnet , daß der optische Meßsensor (10) einen Lichtsender und Lichtempfänger zur nicht-invasiven Licht- Transmissions- oder Reflexionsmessung eingestrahlten Lichtes an geeigneten Körperorten, insbesondere am Finger, Ohrläppchen, Nasenrücken, Wangenschleimhaut oder Stirn umfaßt.
EP99942720A 1998-07-01 1999-07-01 Bestimmung der leberfunktion anhand einer plasmaverschwinderate Expired - Lifetime EP1022984B1 (de)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
DE19829273 1998-07-01
DE19829273 1998-07-01
PCT/DE1999/001899 WO2000001297A1 (de) 1998-07-01 1999-07-01 Bestimmung der leberfunktion anhand einer plasmaverschwinderate

Publications (2)

Publication Number Publication Date
EP1022984A1 true EP1022984A1 (de) 2000-08-02
EP1022984B1 EP1022984B1 (de) 2004-04-14

Family

ID=7872560

Family Applications (1)

Application Number Title Priority Date Filing Date
EP99942720A Expired - Lifetime EP1022984B1 (de) 1998-07-01 1999-07-01 Bestimmung der leberfunktion anhand einer plasmaverschwinderate

Country Status (5)

Country Link
US (1) US6640129B1 (de)
EP (1) EP1022984B1 (de)
DE (1) DE59909175D1 (de)
ES (1) ES2220105T3 (de)
WO (1) WO2000001297A1 (de)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10143995A1 (de) 2001-09-07 2003-04-03 Pulsion Medical Sys Ag System und Computerprogramm zur Bestimmung von Kreislaufgrößen eines Patienten
GB0808777D0 (en) 2008-05-15 2008-06-18 Norgine Bv Prognostic method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0657216B2 (ja) 1988-09-14 1994-08-03 住友電気工業株式会社 肝機能検査装置
JPH02309929A (ja) * 1989-05-24 1990-12-25 Sumitomo Electric Ind Ltd 肝機能検査装置
EP0403683A1 (de) 1989-06-23 1990-12-27 Siemens Aktiengesellschaft Anordnung zur Untersuchung eines flüssigen Messmediums
DE4130931C2 (de) * 1991-09-13 1994-05-19 Hoeft Andreas Verfahren und Vorrichtung zum Ermitteln des zirkulierenden Blutvolumens
GB9716962D0 (en) * 1997-08-12 1997-10-15 Univ Birmingham Liver function test

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO0001297A1 *

Also Published As

Publication number Publication date
DE59909175D1 (de) 2004-05-19
US6640129B1 (en) 2003-10-28
WO2000001297A1 (de) 2000-01-13
ES2220105T3 (es) 2004-12-01
EP1022984B1 (de) 2004-04-14

Similar Documents

Publication Publication Date Title
EP1144024B1 (de) Vorrichtung zur bestimmung des verteilungsvolumens eines blutinhaltsstoffes während einer blutbehandlung
DE3938662C2 (de)
DE3382674T2 (de) Sauerstoffmesseinrichtung mit sauerstoffmesssonde.
DE3851251T2 (de) Rückstrahlungsgerät zur Bestimmung des Sauerstoffgehaltes.
EP1835950B1 (de) Vorrichtung und verfahren zur erkennung von komplikationen wàhrend einer extrakorporalen blutbehandlung
DE69820367T2 (de) Bestimmung von Luftwegwiderstand und Lungennachgiebigkeit
EP0845273B1 (de) Vorrichtung zur in-vivo-Bestimmung von Parametern der Hämodialyse
DE68921947T2 (de) Vorrichtung zum Untersuchen der Leberfunktion.
DE68909635T2 (de) Vorrichtung zur Messung des Metabolismus eines Gewebes.
EP0249825A1 (de) Vorrichtung zur Impedanzmessung an Körpergeweben
DE2741981B2 (de) Vorrichtung zum Bestimmen der Sauerstoffsättigung im Blut
WO1998008434A1 (de) Verfahren und anordnung zur nicht invasiven bestimmung des zerebralen blutflusses mittels nah-infrarot-spektroskopie
DE102004005086A1 (de) Vorrichtung zum Messen der Konzentration einer lichtabsorbierenden Substanz in Blut
DE102010034626A1 (de) Vorrichtung zur extrakorporalen Blutbehandlung
EP0793113B1 (de) MR-Verfahren mit reduzierten Bewegungsartefakten
DE68911686T2 (de) Vorrichtung und verfahren zur messung der aktivität von radioaktiven mustern, die mehrere radioaktive isotope enthalten, ohne separate bestimmung des löschniveaus.
DE69029023T2 (de) Kalibrierungsverfahren für einen Impulsantwort-Durchflussmesser
DE4130931C2 (de) Verfahren und Vorrichtung zum Ermitteln des zirkulierenden Blutvolumens
DE102006058332A1 (de) Verfahren zur Steuerung einer Messung eines Magnetresonanzgeräts anhand eines Elektrokardiogramm-Signals
DE60029754T2 (de) Vorrichtung zur bestimmung des blutwiederumlaufes in einer vaskulären zugangseinrichtung
DE3887639T2 (de) Vorrichtung zum Untersuchen der Leberfunktion.
EP1022984A1 (de) Bestimmung der leberfunktion anhand einer plasmaverschwinderate
DE3787466T2 (de) Gerät zur untersuchung der leberfunktion.
DE2754334C2 (de) Vorrichtung zur Bestimmung des Blutdrucks
DE3209810A1 (de) Verfahren zum messen der magnetischen kernresonanz fuer die nmr-tomographie

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20000114

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LI LU MC NL PT SE

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): DE ES FR GB IT

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REF Corresponds to:

Ref document number: 59909175

Country of ref document: DE

Date of ref document: 20040519

Kind code of ref document: P

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: GERMAN

GBT Gb: translation of ep patent filed (gb section 77(6)(a)/1977)

Effective date: 20040701

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2220105

Country of ref document: ES

Kind code of ref document: T3

Ref country code: IE

Ref legal event code: FD4D

ET Fr: translation filed
PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20050117

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: IT

Payment date: 20110728

Year of fee payment: 13

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20130215

Year of fee payment: 14

Ref country code: DE

Payment date: 20130124

Year of fee payment: 14

Ref country code: GB

Payment date: 20130125

Year of fee payment: 14

Ref country code: ES

Payment date: 20130129

Year of fee payment: 14

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20130701

REG Reference to a national code

Ref country code: FR

Ref legal event code: ST

Effective date: 20140331

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20140201

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130701

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 59909175

Country of ref document: DE

Effective date: 20140201

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130701

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130731

REG Reference to a national code

Ref country code: ES

Ref legal event code: FD2A

Effective date: 20150710

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20130702