Apparatus for determining a parameter indicative of the progress of an extracorporeal blood treatment

Abstract

An apparatus for extracorporeal treatment of blood (1) comprising a treatment unit, a blood withdrawal line, a blood return line, a preparation line and a spent dialysate line. A control unit (10) is configured to calculate values of a parameter relating to treatment effectiveness based on measures of the conductivity in the spent dialysate line. The value of the effectiveness parameter is calculated using one or more values representative of the conductivity in the spent dialysate line obtained relying on a mathematical model.

Claims

1. An apparatus for extracorporeal treatment of blood comprising: a preparation line having one end connected to an inlet of a secondary chamber of a treatment unit and configured to convey fresh treatment liquid to the secondary chamber, the fresh treatment liquid presenting a characteristic (Cd.sub.in) which is one selected from the group of: conductivity in the fresh treatment liquid, and concentration of at least one substance in the fresh treatment liquid; a spent dialysate line having one end connected to an outlet of said secondary chamber and configured to remove spent liquid from the secondary chamber, the spent liquid presenting a characteristic (Cd.sub.out) which is one selected from the group of: conductivity in the spent liquid, and concentration of at least one substance in the spent liquid; and a control unit configured for commanding execution of a task for determination of a parameter indicative of the effectiveness of the extracorporeal blood treatment, said task comprising the following steps: receiving at least one prescription baseline (Cd.sub.set) for the characteristic (Cd.sub.in) in the fresh treatment liquid; causing fresh treatment liquid to flow in the preparation line to the secondary chamber with the characteristic being at said prescription baseline (Cd.sub.set); causing spent liquid to flow out of the secondary chamber into the spent dialysate line; receiving at least one parametric mathematical model which puts into relation the characteristic (Cd.sub.in) in the fresh treatment liquid with the characteristic (Cd.sub.out) in the spent liquid, said parametric mathematical model presenting a prefixed number of free parameters; causing an upstream variation of the value of the characteristic (Cd.sub.in) in the fresh treatment liquid with respect to said prescription baseline (Cd.sub.set) thereby causing a corresponding and timely delayed downstream variation of the same characteristic (Cd.sub.out) in the spent liquid flowing in the spent dialysate line; measuring a plurality of values taken by a reference portion of said downstream variation of the characteristic (Cd.sub.out) in the spent liquid, said reference portion having duration shorter than the entire duration of the downstream variation; estimating said free parameters of the at least one parametric mathematical model by means of said reference portion measured values and identifying one single characteristic mathematical model which puts into relation the characteristic (Cd.sub.in) in the fresh treatment liquid with the characteristic (Cd.sub.out) in the spent liquid; and computing at least one value of a parameter indicative of the effectiveness of the extracorporeal blood treatment by using said characteristic mathematical model and one or more values taken by the characteristic (Cd.sub.in) in the fresh treatment liquid.

2. Apparatus according to claim 1, wherein the step of computing at least one value of a parameter indicative of the effectiveness of the extracorporeal blood treatment comprises: computing one or more significant values of said downstream variation of the characteristic (Cd.sub.out), said significant value or values of the downstream variation relating to a time subsequent to the duration of the reference portion and being computed by using said characteristic mathematical model; computing at least one value of a parameter indicative of the effectiveness of the extracorporeal blood treatment from said computed significant value or values and from one or more values taken by the characteristic (Cd.sub.in) in the fresh treatment liquid.

3. Apparatus according to claim 2, wherein said reference portion has a duration which is less than 50% compared to the duration of the entire downstream perturbation.

4. Apparatus according to claim 1, wherein the control unit is configured to set, or to allow setting of, the duration of the said reference portion.

5. Apparatus according to claim 1, wherein computing one or more significant values of said delayed variation of the characteristic (Cd.sub.out) comprises determining the value (Cd.sub.out(n)) of the characteristic (Cd.sub.out) in the spent liquid at time instant (n) by using as input to the characteristic mathematical model: the values of the characteristic (Cd.sub.in) in the fresh treatment liquid at a plurality of time instants (n1, n2, n3) preceding in time the time instant (n); or a mathematically calculated version of the characteristic (Cd.sub.in) in the fresh treatment liquid.

6. Apparatus according to claim 5, wherein computing one or more significant values of said delayed variation of the characteristic (Cd.sub.out) comprises determining the value (Cd.sub.out(n)) of the characteristic (Cd.sub.out) in the spent liquid at time instant (n) subsequent to said reference portion with the following recursive equation, which represents the parametric mathematical model in the time domain:
Cd.sub.out(n)=a.sub.0.Math.Cd.sub.in(n)+b.sub.1.Math.Cd.sub.out(n1)+b.sub.2.Math.Cd.sub.out(n2)+ . . . b.sub.m.Math.Cd.sub.out(nm), wherein: Cd.sub.out(n) is the calculated value of the outlet characteristic at time instant (n), Cd.sub.in(n) is the known value of the inlet characteristic at time instant (n), Cd.sub.out(n1), Cd.sub.out(n2), . . . , Cd.sub.out(nm) are values of the outlet characteristic at preceding time instants (n1, n2, . . . nm) prior to time instant (n) and recursively computed by means of the mathematical model, and a.sub.0, b.sub.1, b.sub.2, . . . , b.sub.m are constant parameters of the mathematical model, as estimated by using said measured values of the reference portion of the downstream variation.

7. Apparatus according to claim 6, wherein the at least one mathematical model is a time invariant linear (LTI) model; and whereinin the frequency domain and using the z-Transformthe mathematical model is described by a transfer function H(z) having at least one zero and at least one pole.

8. Apparatus according to claim 7, wherein the transfer function H(z) comprises from 3 to 5 poles and is described by one of the following:
H(z)=Cd.sub.out(z)/Cd.sub.in(z)=a.sub.0/(1b.sub.1.Math.z.sup.1b.sub.2.Math.z.sup.2b.sub.3.Math.z.sup.3b.sub.4.Math.z.sup.4b.sub.5.Math.z.sup.5),
H(z)=Cd.sub.out(z)/Cd.sub.in(z)=a.sub.0/(1b.sub.1.Math.z.sup.1b.sub.2.Math.z.sup.2b.sub.3.Math.z.sup.3b.sub.4.Math.z.sup.4),
H(z)=Cd.sub.out(z)/Cd.sub.in(z)=a.sub.0/(1b.sub.1.Math.z.sup.1b.sub.2.Math.z.sup.2b.sub.3.Math.z.sup.3),
H(z)=Cd.sub.out(z)/Cd.sub.in(z)=a.sub.0/(1b.sub.1.Math.z.sup.1b.sub.2.Math.z.sup.2), wherein a.sub.0, b.sub.1, b.sub.2, b.sub.3, b.sub.4, b.sub.5 are constant parameters of the model, as estimated by using said measured values of the reference portion of the downstream variation.

9. Apparatus according to claim 1 comprising: a memory connected to the control unit and storing one or more change setting procedures, each of said change setting procedures when executed by the control unit configuring the control unit to vary the value of the characteristic (Cd.sub.in) of the fresh dialysis liquid from said set value to a new set value, said control unit being further configured to: prevent execution of the change setting procedure(s) only until end of a measurement instant (T.sub.END.sub._.sub.MEAS) at which the measurement of said plurality of values of the reference portion of said downstream variation in the characteristic (Cd.sub.out) in the spent liquid has been completed, and allow execution of the change setting procedure(s) immediately after the end of the measurement instant (T.sub.END-MEAS).

10. Apparatus according to claim 9, wherein the control unit is further configured to: receive a total treatment time (T), consecutively repeating at time intervals said task for the determination of the parameter during the treatment time (T) such that a plurality of values indicative of said parameter are correspondingly determined, and allow execution of one or more change setting procedures between two consecutive tasks.

11. Apparatus according claim 1, wherein: varying the value of the characteristic (Cd.sub.in) in the fresh treatment liquid comprises imposing a change from a first inlet value (Cd.sub.in1) to a second inlet value (Cd.sub.in2), which is kept constant for a prefixed time interval, thereby causing a corresponding change of the characteristic (Cd.sub.out) in spent liquid from a respective first outlet value (Cd.sub.out1) to a respective second outlet value (Cd.sub.out2) defining said timely delayed downstream variation of the characteristic (Cd.sub.out), the reference portion of said downstream variation begins after the characteristic in the spent liquid changes from said first outlet value (Cd.sub.out1) and lasts a period during which the characteristic either continuously increases or decreases without reaching the second outlet value (Cd.sub.out2), and the control unit is further configured to: calculate the second outlet value (Cd.sub.out2) of the characteristic (Cd.sub.out) by using as input to the characteristic mathematical model the values of the characteristic (Cd.sub.in) in the fresh treatment liquid, or a mathematically calculated version of the characteristic (Cd.sub.in) in the fresh treatment liquid, and use the calculated second outlet value (Cd.sub.out2) as a significant value for the computation of at least one value of a parameter indicative of the effectiveness of the extracorporeal blood treatment.

12. Apparatus according to claim 1, wherein: varying the value of the characteristic (Cd.sub.in) in the fresh treatment liquid comprises changing from a first inlet value (Cd.sub.in1) to a second inlet value (Cd.sub.in2), and then changing to a third inlet value (Cd.sub.out3) thereby causing a corresponding time delayed downstream variation of the characteristic (Cd.sub.out) in the spent liquid comprising a change from a respective first outlet value (Cd.sub.out1) to a respective second outlet value (Cd.sub.out2) and then to a third out value (Cd.sub.out3); the reference portion of the downstream variation begins after the characteristic in the spent liquid changes from said first outlet value (Cd.sub.out1) and lasts a period shorter than a fraction of the duration of said downstream variation; and calculating, as a significant value at least the third outlet value (Cd.sub.out3) or both the second and the third outlet values (Cd.sub.out2, Cd.sub.out3) of the characteristic (Cd.sub.out) by using as input to the characteristic mathematical model the values of the characteristic (Cd.sub.in) in the fresh treatment liquid, or a mathematically calculated version of the characteristic (Cd.sub.in) in the fresh treatment liquid.

13. Apparatus according to claim 1, wherein: varying the value of the characteristic (Cd.sub.in) in the fresh treatment liquid comprises imposing an upstream perturbation in the characteristic (Cd.sub.in) of the fresh treatment liquid thereby causing a corresponding downstream perturbation of the characteristic (Cd.sub.out) in spent liquid, the reference portion of said downstream perturbation begins after the characteristic (Cd.sub.out) in the spent liquid changes from said first outlet value (Cd.sub.out1) and lasts a prefixed period shorter than 60% of the duration of the downstream perturbation, and the control unit is further configured to: extrapolate a plurality of significant values of the characteristic (Cd.sub.out), describing a remaining portion of the downstream perturbation consecutive to said reference portion, by using as input to the mathematical model the values of the characteristic (Cd.sub.in) in the fresh treatment liquid, or a mathematically calculated version of the characteristic (Cd.sub.in) in the fresh treatment liquid, obtain a calculated downstream perturbation from said extrapolated significant values, and compute at least one value of a parameter indicative of the effectiveness of the extracorporeal blood treatment by comparing the calculated downstream perturbation and the upstream perturbation.

14. Apparatus according to claim 1, wherein the reference portion of said downstream variation lasts a period during which the characteristic (Cd.sub.out) either continuously increases or decreases without reaching 80% of the second outlet value (Cd.sub.out2).

15. Apparatus according to claim 1, wherein said parameter comprises one selected from the group of: an effective dialysance for one or more substances of the treatment unit (D), an effective clearance for one or more substances of the treatment unit (K), a concentration of a substance in blood (Cb.sub.in) upstream the blood treatment unit, and a dialysis dose at time (t) after start of the treatment (K.Math.t/V).

16. Apparatus according to claim 15, wherein the parameter comprises the effective dialysance (D) and wherein each computed value (D.sub.k) of said parameter at each respective variation is obtained using the formula:
D.sub.k=(Qd+WLR).Math.[1(Cd.sub.out2Cd.sub.out1)]/(Cd.sub.in2Cd.sub.in1) where: Cd.sub.out1 is the first outlet value taken by the characteristic (Cd.sub.out) in the spent dialysate line downstream of the secondary chamber in response to the change of the characteristic (Cd.sub.in) in the preparation line to said first inlet value Cd.sub.in1, Cd.sub.out2 is the second value taken by the characteristic (Cd.sub.out) in the spent dialysate line downstream of the secondary chamber in response to the change of the characteristic (Cd.sub.in) in the preparation line at said second inlet value (Cd.sub.in2), Cd.sub.in1, Cd.sub.in2 are first and second inlet values taken by the characteristic (Cd.sub.in) in the preparation line upstream of the secondary chamber, Qd is the fresh treatment liquid flow rate in the preparation line, WLR is the weight loss rate of a patient under treatment, or the formula,
D.sub.K=(Qd+WLR)[1(2Cd.sub.out1Cd.sub.out2Cd.sub.out3)/(2Cd.sub.in1Cd.sub.in2Cd.sub.in3)] where: Cd.sub.out1 is the first outlet value taken by the characteristic in the spent dialysate line downstream of the secondary chamber in response to the change of the characteristic (Cd.sub.in) in the preparation line to said first inlet value (Cd.sub.in1), Cd.sub.out2 is the calculated second value which is representative of the value taken by the characteristic in the spent dialysate line downstream of the secondary chamber in response to the change of the characteristic (Cd.sub.in) in the preparation line from said first inlet value (Cd.sub.in1) to said second inlet value (Cd.sub.in2), Cd.sub.out3 is the calculated third value which is representative of the value taken by the characteristic in the spent dialysate line downstream of the secondary chamber in response to the change of the characteristic (Cd.sub.in) in the preparation line from said second inlet value (Cd.sub.in2) to said third inlet value (Cd.sub.in3), Cd.sub.in1, Cd.sub.in2, Cd.sub.in3 are first, second and third inlet values taken by the characteristic (Cd.sub.in) in the preparation line upstream of the secondary chamber, Qd is the fresh treatment liquid flow rate in the preparation line, and WLR is the weight loss rate of a patient under treatment.

17. Apparatus according to claim 1, wherein the control unit is configured to: determine a downstream curve (Cd.sub.out(t)) representative of the values taken by the characteristic (Cd.sub.out) in the spent liquid, determine an angular deviation of a downstream baseline of the downstream curve (Cd.sub.out(t)) with respect to the prescription baseline (Cd.sub.set), and compensate for said angular deviation by angularly rotating the downstream curve to obtain a corrected downstream curve (Cd.sub.out-correct(t)).

18. Apparatus according to claim 17, wherein the step of determining a downstream curve (Cd.sub.out(t)) representative of the values taken by the characteristic (Cd.sub.out) in the spent liquid comprises: using measured values across a time interval covering at least the entire downstream variation; or measuring values of the characteristic (Cd.sub.out) in the spent liquid until the end of said reference portion, estimating the free parameters of the parametric mathematical model to identify the characteristic mathematical model, and using said identified characteristic mathematical model to calculate the downstream curve (Cd.sub.out(t)).

19. Apparatus according to claim 1, wherein the control unit is configured to: determine a downstream curve (Cd.sub.out(t)) representative of the values taken by the characteristic (Cd.sub.out) in the spent liquid, analyze a frequency spectrum of the downstream curve (Cd.sub.out(t)), and filter out harmonics of said frequency spectrum of the downstream curve (Cd.sub.out(t)) lying at frequencies higher than a prefixed threshold to eliminate noise and undesired perturbations possibly present in the downstream curve and obtain a corrected downstream curve (Cd.sub.out-correct(t)).

20. Apparatus according to claim 19, wherein the step of determining a downstream curve (Cd.sub.out(t)) representative of the values taken by the characteristic (Cd.sub.out) in the spent liquid comprises: using measured values across a time interval covering at least the entire downstream variation; or measuring values of the characteristic (Cd.sub.out) in the spent liquid until the end of said reference portion, estimating the free parameters of the parametric mathematical model to identify the characteristic mathematical model, and using said identified characteristic mathematical model to calculate the downstream curve (Cd.sub.out(t)).

21. Apparatus according to claim 1 wherein the control unit is further configured to: store a plurality of mathematical models each of which puts into relation the characteristic (Cd.sub.in) in the fresh treatment liquid with the characteristic (Cd.sub.out) in the spent liquid; and select the mathematical model to be used for computing the at least one significant value of said downstream variation based at least on one selected from the group of the shape of said upstream variation and the type of blood treatment unit used by the apparatus.

22. An apparatus for extracorporeal treatment of blood comprising: a blood treatment unit having a primary chamber and a secondary chamber separated by a semi-permeable membrane; a preparation line having one end connected to an inlet of the secondary chamber of the blood treatment unit and configured to convey fresh treatment liquid to the secondary chamber, the fresh treatment liquid presenting a characteristic (Cd.sub.in) which is one selected the group of: conductivity in the fresh treatment liquid, and concentration of at least one substance in the fresh treatment liquid; a spent dialysate line having one end connected to an outlet of said secondary chamber and configured to remove spent liquid from the secondary chamber, the spent liquid presenting a characteristic (Cd.sub.out) which is one selected from the group of: conductivity in the spent liquid, and concentration of at least one substance in the spent liquid; and a control unit configured for commanding execution of a task for determination of a parameter indicative of the effectiveness of the extracorporeal blood treatment, said task comprising the following steps: receiving at least one prescription baseline (Cd.sub.set) for the characteristic (Cd.sub.in) in the fresh treatment liquid; causing fresh treatment liquid to flow in the preparation line to the secondary chamber with the characteristic being at said prescription baseline (Cd.sub.set); causing spent liquid to flow out of the secondary chamber into the spent dialysate line; causing an upstream variation of the value of the characteristic (Cd.sub.in) in the fresh treatment liquid with respect to said prescription baseline (Cd.sub.set) thereby causing a corresponding and time delayed downstream variation of the same characteristic (Cd.sub.out) in the spent liquid flowing in the spent dialysate line; measuring a plurality of values taken by a reference portion of said downstream variation of the characteristic (Cd.sub.out) in the spent liquid, determining a downstream curve (Cd.sub.out(t)) representative of the values taken by the characteristic (Cd.sub.out) in the spent liquid, determining an angular deviation of a downstream baseline of the downstream curve (Cd.sub.out(t)) with respect to the prescription baseline (Cd.sub.set), compensating for said angular deviation by angularly rotating the downstream curve to obtain a corrected downstream curve (Cd.sub.out-correct(t)), computing at least one value of a parameter indicative of the effectiveness of the extracorporeal blood treatment by using one or more values of said corrected downstream curve and one or more values taken by the characteristic (Cd.sub.in) in the fresh treatment liquid.

23. Apparatus according to claim 22, wherein said parameter comprises one selected from the group of: an effective dialysance for one or more substances of the treatment unit (D), an effective clearance for one or more substances of the treatment unit (K), a concentration of a substance in blood (Cb.sub.in) upstream the blood treatment unit, and a dialysis dose at time (t) after start of the treatment (K.Math.t/V).

24. Apparatus according to claim 22, wherein the step of determining a downstream curve (Cd.sub.out(t)) representative of the values taken by the characteristic (Cd.sub.out) in the spent liquid comprises using measured values across a time interval covering at least the entire downstream variation.

25. Apparatus according to claim 22, comprising: receiving at least one parametric mathematical model which puts into relation the characteristic (Cd.sub.in) in the fresh treatment liquid with the characteristic (Cd.sub.out) in the spent liquid, said parametric mathematical model presenting a prefixed number of free parameters, wherein the step of determining a downstream curve (Cd.sub.out(t)) representative of the values taken by the characteristic (Cd.sub.out) in the spent liquid comprises: measuring values of the characteristic (Cd.sub.out) in the spent liquid until the end of said reference portion, said reference portion having duration shorter than the entire duration of the downstream variation, estimating the free parameters of the parametric mathematical model to identify the characteristic mathematical model, and using said identified characteristic mathematical model to calculate the downstream curve (Cd.sub.out(t)).

Description

DESCRIPTION OF THE DRAWINGS

(1) Aspects of the invention are shown in the attached drawings, which are provided by way of non-limiting example, wherein:

(2) FIG. 1 shows a schematic diagram of a blood treatment apparatus according to one aspect of the invention;

(3) FIG. 2 shows a schematic diagram of an alternative embodiment of a blood treatment apparatus according to another aspect of the invention;

(4) FIG. 3 shows a schematic diagram of another alternative embodiment of a blood treatment apparatus according to a further aspect of the invention;

(5) FIG. 4 shows a conductivity (or concentration) vs. time diagram showing the conductivity (or concentration) profile in the fresh and in the spent dialysate line, according to another aspect of the invention;

(6) FIGS. 5 and 6 show a conductivity (or concentration) vs. time diagram in the fresh and in the spent dialysate line wherein the conductivity (or concentration) variation in the fresh dialysate line is in the form of a relatively long step;

(7) FIGS. 7 and 8 show a conductivity (or concentration) vs. time diagram in the fresh and in the spent dialysate line wherein the conductivity (or concentration) variation in the fresh dialysate line is in the form of a relatively short step;

(8) FIG. 9 shows a conductivity (or concentration) vs. time diagram in the fresh and in the spent dialysate line wherein the conductivity (or concentration) variation in the fresh dialysate line is in the form of pulse;

(9) FIG. 10 is a diagram showing the outlet conductivity (or concentration) vs. time and schematically illustrating an angular correction of the conductivity (or concentration) baseline;

(10) FIG. 11 is a schematic flowchart of a method according to one aspect of the invention;

(11) FIG. 12 shows a conductivity (mS.Math.100/cm) vs. time (seconds) diagram showing the real measured conductivity profile in the fresh and in the spent dialysate line in the case of a step conductivity variation in the fresh dialysate of 1 mS/cm.

(12) FIG. 13 is an enlarged view of the measured outlet conductivity profile of FIG. 12;

(13) FIG. 14 is an enlarged view of the measured outlet conductivity profile of FIG. 12 where the reference time T.sub.R is identified;

(14) FIG. 15 is an enlarged view of the measured outlet conductivity profile of FIG. 12 (dotted line) and of the calculated curve representing the outlet conductivity as determined using a mathematical model according to aspects of the invention.

DETAILED DESCRIPTION

(15) Non-limiting embodiments of an apparatus 1 for extracorporeal treatment of bloodwhich may implement innovative aspects of the inventionare shown in FIGS. 1 to 3. FIG. 1 is a more schematic representation of the extracorporeal blood treatment apparatus 1, while FIGS. 2 and 3 represent, in greater detail, two possible non limiting embodiments of the apparatus 1.

(16) The apparatus 1 may be configured to determine a parameter indicative of the effectiveness of the treatment delivered to a patient (here below also referred to as effectiveness parameter). The effectiveness parameter may be one of the following: an effective dialysance for one or more substances of the treatment unit (D), e.g. electrolyte or sodium clearance; an effective clearance for one or more substances of the treatment unit (K), e.g. urea clearance; a concentration of a substance in blood (Cbin) upstream the blood treatment unit, e.g. sodium concentration in the blood upstream the treatment unit; a dialysis dose delivered until a certain point in time after start of the treatment (K.Math.t/V), where K is clearance, t represents the time interval from start of treatment until the point in time, and V represents a reference volume characteristic of the patient.

(17) Note that a parameter proportional to one of the above parameters or known function of one or more of the above parameters may alternatively be used as effectiveness parameter.

(18) In below description and in FIGS. 1 to 3 same components are identified by same reference numerals.

(19) In FIG. 1 it is represented an apparatus for the extracorporeal treatment of blood 1 comprising a treatment unit 2 (such as an hemofilter, an ultrafilter, an hemodiafilter, a dialyzer, a plasmafilter and the like) having a primary chamber 3 and a secondary chamber 4 separated by a semi-permeable membrane 5; depending upon the treatment, the membrane of the filtration unit may be selected to have different properties and performances. A blood withdrawal line 6 is connected to an inlet of the primary chamber 3, and a blood return line 7 is connected to an outlet of the primary chamber 3. In use, the blood withdrawal line 6 and the blood return line 7 are connected to a needle or to a catheter or other access device (not shown) which is then placed in fluid communication with the patient vascular system, such that blood may be withdrawn through the blood withdrawal line, flown through the primary chamber and then returned to the patient's vascular system through the blood return line. An air separator, such as a bubble trap 8, may be present on the blood return line; moreover, a safety clamp 9 controlled by a control unit 10 may be present on the blood return line downstream the bubble trap 8. A bubble sensor 8a, for instance associated to the bubble trap 8 or coupled to a portion of the line 7 between bubble trap 8 and clamp 9 may be present: if present, the bubble sensor is connected to the control unit 10 and sends to the control unit 10 signals for the control unit 10 to cause closure of the clamp 9 in case one or more bubbles above certain safety thresholds are detected. The blood flow through the blood lines is controlled by a blood pump 11, for instance a peristaltic blood pump, acting either on the blood withdrawal line or on the blood return line.

(20) An operator may enter a set value for the blood flow rate Q.sub.B through a user interface 12 and the control unit 10, during treatment, is configured to control the blood pump based on the set blood flow rate. The control unit 10 may comprise a digital processor (CPU) and a memory (or memories), an analogical type circuit, or a combination thereof as explained in greater detail in below section dedicated to the control unit 10.

(21) An effluent fluid line or spent dialysate line 13 is connected, at one end, to an outlet of the secondary chamber 4 and, at its other end, to a waste which may be a discharge conduit or an effluent fluid container collecting the fluid extracted from the secondary chamber. A fresh dialysis fluid line 19 is connected to the inlet of the secondary chamber 4 and supplies fresh dialysate to from a source to said second chamber. Conductivity or concentration sensors 109, 109a are respectively positioned on the fresh dialysis fluid line 19 and on the spent dialysate line 13. Concentration or conductivity sensor 109 is configured for detecting the conductivity or the concentration for one substance of for a group of substancesidentified as Cd.sub.inin the fresh dialysis fluid line 19. Concentration or conductivity sensor 109a is configured for detecting the conductivity or the concentration for one substance of for a group of substancesidentified as Cd.sub.outin the spent dialysate line 19.

(22) FIG. 2 shows an apparatus 1 configured to deliver any one of treatments like ultrafiltration and hemodialysis and hemodiafiltration. The apparatus of FIG. 2 comprises all the features described above in connection with FIG. 1, which are identified in FIG. 2 with same reference numerals. Furthermore, in the apparatus of FIG. 2, other features of a possible embodiment of the apparatus 1 are schematically shown: an effluent fluid pump 17 that operates on the effluent fluid line under the control of control unit 10 to regulate the flow rate Q.sub.eff across the effluent fluid line. The apparatus 1 may also include an ultrafiltration line 25 branching off the effluent line 13 and provided with a respective ultrafiltration pump 27 also controlled by control unit 10. The embodiment of FIG. 2 presents a pre-dilution fluid line 15 connected to the blood withdrawal line: this line 15 supplies replacement fluid from an infusion fluid container 16 connected at one end of the pre-dilution fluid line. Although in FIG. 2 a container 16 is shown as the source of infusion fluid, this should not be interpreted in a limitative manner: indeed, the infusion fluid may also come from an on line preparation section 100 part of the apparatus 1. Note that alternatively to the pre-dilution fluid line the apparatus of FIG. 1 may include a post-dilution fluid line (not shown in FIG. 2) connecting an infusion fluid container to the blood return line Finally, as a further alternative (not shown in FIG. 2) the apparatus of FIG. 2 may include both a pre-dilution and a post infusion fluid line: in this case each infusion fluid line may be connected to a respective infusion fluid container or the two infusion fluid lines may receive infusion fluid from a same source of infusion fluid such as a same infusion fluid container. Once again, the source of infusion fluid may alternatively be an online preparation section part of the apparatus 1 (similar to the device 100 described herein below) supplying fluid to the post and/or pre dilution lines. Furthermore, an infusion pump 18 operates on the infusion line 15 to regulate the flow rate Q.sub.rep through the infusion line. Note that in case of two infusion lines (pre-dilution and post-dilution) each infusion line may be provided with a respective infusion pump. The apparatus of FIG. 2, which includes a dialysis fluid line 19 connected at one end with a water inlet and at its other end with the inlet of the secondary chamber 4 of the filtration unit for supplying fresh dialysis liquid to the secondary chamber 4; a dialysis fluid pump 21 is operative on the dialysis liquid fluid line 19 under the control of said control unit 10, to supply fluid from the dialysis liquid container to the secondary chamber at a flow rate Q.sub.dial. The dialysis fluid pump 21, the ultrafiltration pump 27, the concentrate pumps 105 and 108, the infusion fluid pump 15 and the effluent fluid pump 17 are operatively connected to the control unit 10 which controls the pumps as it will be in detail disclosed herein below. An initial tract of line 19 links the haemodialyser or hemodiafilter 2 to a device 100, for preparing the dialysis liquid, which also includes a further tract of said line 19. The device 100 comprises a main line 101, the upstream end of which is designed to be connected to a supply of running water. Connected to this main line 101 are a first secondary line 102 and a second secondary line 103. The first secondary line 102, which may be looped back onto the main line 101, is provided with a connector configured for fitting a container 104, such as a bag or cartridge or other container, containing sodium bicarbonate in granule form (alternatively a concentrate in liquid form may be used). Line 102 is furthermore equipped with a concentrate pump 105 for metering the sodium bicarbonate into the dialysis liquid: as shown in FIG. 7 the pump may be located downstream of the container 104. The operation of the pump 105 is determined by the comparison between 1) a conductivity set point value for the solution forming at the junction of the main line 101 and the first secondary line 102 and 2) the value of the conductivity of this mixture measured by means of a first conductivity probe 106 located in the main line 101 immediately downstream of the junction between the main line 101 and the first secondary line 102. The free end of the second secondary line 103 is intended to be immersed in a container 107 for a concentrated saline solution, e.g. containing sodium chloride, calcium chloride, magnesium chloride and potassium chloride, as well as acetic acid. The second secondary line 103 is equipped with a pump 108 for metering sodium into the dialysis liquid, the operation of which pump depends on the comparison between 1) a second conductivity set point value for the solution forming at the junction of the main line 101 and the second secondary line 103 and 2) the value of the conductivity of this solution measured by means of a second conductivity probe 109 located in the main line 12 immediately downstream of the junction between the main line 12 and the secondary line 103. Note that as an alternative, instead of conductivity sensors concentration sensors may in principle be used. Moreover, the specific nature of the concentrates contained in containers 104 and 107 may be varied depending upon the circumstances and of the type of dialysis fluid to be prepared.

(23) The control unit 10 is also connected to the user interface 12, for instance a graphic user interface, which receives operator's inputs and displays the apparatus outputs. For instance, the graphic user interface 12 may include a touch screen, a display screen and hard keys for entering user's inputs or a combination thereof.

(24) The embodiment of FIG. 3 shows an alternative apparatus 2 designed for delivering any one of treatments like hemodialysis and ultrafiltration. The apparatus of FIG. 3 includes the same components described for the apparatus of FIG. 1. In the apparatus shown in FIG. 3 the same components described for the embodiment of FIG. 2 are identified by same reference numerals and thus not described again. In practice, differently from the hemodiafiltration apparatus of FIG. 2, the apparatus of FIG. 3 does not present any infusion line.

(25) In each one of the above described embodiments, flow sensors 110, 111 (either of the volumetric or of the mass type) may be used to measure flow rate in each of the lines. Flow sensors are connected to the control unit 10. In the example of FIG. 2 where the infusion line 15 and the ultrafiltration line 25 lead to a respective bag 16, 23, scales may be used to detect the amount of fluid delivered or collected. For instance, the apparatus of FIG. 2 includes a first scale 33 operative for providing weight information W.sub.1 relative to the amount of the fluid collected in the ultrafiltration container 23 and a second scale 34 operative for providing weight information W.sub.2 relative to the amount of the fluid supplied from infusion container 16. In the embodiment of FIG. 3, the apparatus includes a first scale 33 operative for providing weight information W.sub.1 relative to the amount of the fluid collected in the ultrafiltration container 23. The scales are all connected to the control unit 10 and provide said weight information W.sub.i for the control unit 10 to determine the actual quantity of fluid in each container as well as the actual flow rate of fluid supplied by, or received in, each container. In the example of FIG. 1 there is no dedicated ultrafiltration line and the total amount of ultrafiltration is determined by the difference of the flow rates detected by sensors 110 and 111. The control unit 10 is configured to act on appropriate actuators, such as pumps, present on lines 13 and 19 andusing the information concerning the difference of flow rates as detected by sensors 110, 111to make sure that a prefixed patient fluid removal is achieved in the course of a treatment time T, as required by the prescription provided to the control unit 10, e.g. via user interface 12.

(26) In the example of FIGS. 2 and 3, in order to control the fluid balance between the quantity of fluid supplied to the secondary chamber 4 and the quantity of fluid extracted from the secondary chamber, the flow-meters 110, 111 positioned on the fresh dialysate line and on the waste line 13 provide the control unit 10 with signals indicative of the flow of fluid through the respective lines and the scale or scales provide weight information which allow the control unit 10 to derive the flow rate through the ultrafiltration line 25 and, if present, through the infusion line 15. The control unit 10 is configured to control at least pumps 17, 21 and 27 (in case of FIG. 2 also pump 18) to make sure that a prefixed patient fluid removal is achieved in the course of a treatment time T, as required by the prescription provided to the control unit 10, e.g. via user interface 12. Note that other fluid balance systems may be used: for instance in case the apparatus includes a container as source of fresh dialysis fluid and a container to collect waste, then scales may be used to detect the amount of fluid delivered or collected by each container and then inform the control unit 10 accordingly. As a further alternative, systems based on volumetric control may be used where the fresh dialysis liquid line 19 and the waste line 13 are connected to a balance chamber system assuring thatat each instantthe quantity of liquid flowing into line 19 is identical to the quantity of fluid exiting from line 13.

(27) From a structural point of view one or more, containers 104, 107, 14, 16, 23 may be disposable plastic containers. The blood lines 6, 7 lines and the filtration unit may also be plastic disposable components which may be mounted at the beginning of the treatment session and then disposed of at the end of the treatment session. Pumps, e.g. peristaltic pumps or positive displacement pumps, have been described as means for regulating fluid flow through each of the lines; however, it should be noted that other flow regulating means may alternatively be adopted such as for example valves or combinations of valves and pumps. The scales may comprise piezoelectric sensors, or strain gauges, or spring sensors, or any other type of transducer able to sense forces applied thereon. As already explained, the conductivity sensors may be replaced by concentration sensors.

(28) Determination of the Effectiveness Parameter

(29) As mentioned at the beginning of the detailed description, the apparatus 1 is capable of determining an effectiveness parameter. In this regard, the control unit 10 of the apparatus 1 is configured for commanding execution of a number of procedures including a task specifically devoted to the determination of the parameter indicative of the effectiveness of the extracorporeal blood treatment. The task devoted to determination of the effectiveness parameter comprises the steps described herein below.

(30) First, the control unit 10 is configured for receiving at least one prescription baseline Cd.sub.set for the characteristic Cd.sub.in in the fresh treatment liquid; the characteristic may be the concentration for one substance in the dialysis liquid (e.g. the sodium concentration, or the calcium concentration), or the concentration for a group of substances in the dialysis liquid (such as the electrolyte concentration) or the conductivity of the dialysis liquid. Furthermore, the set value for the prescription baseline may be either pre set in a memory connected to the control unit 10 or, alternatively, it may be entered by the user via user interface 12.

(31) Although the prescription baseline is frequently a constant value, it may alternatively comprise a time-variable value which changes during treatment according to a prefixed law.

(32) The control unit 10, acting on appropriate actuators such as pumps 21 and 17, causes circulation of dialysis fluid through lines 19 and 13 and through the secondary chamber 4 of the treatment unit 2. In greater detail, the control unit 10 is configured for causing fresh treatment liquid to flow in the preparation line 19 to the secondary chamber 4 with the characteristic being at said prescription baseline Cd.sub.set: the characteristic at the baseline value may for instance be achieved by appropriately controlling the concentrate pumps 105, 108 of the preparation section 100. Furthermore, the control unit 10 is configured for reading the value of the characteristic in the spent dialysis fluid using sensor 109a. Depending upon the case, sensor 109a may for instance be a conductivity sensor, or a concentration sensor (sensitive to one or more substances).

(33) The control unit 10 is also configured to receive at least one parametric mathematical model which puts into relation the characteristic Cd.sub.in in the fresh treatment liquid with the characteristic Cd.sub.out in the spent liquid. The parametric mathematical model, which mathematically describes the components interposed between the two sensors 109, 109a, may for instance be pre-stored in a memory connected to the control unit 10, or it may be transferred to said memory via user interface 12 or via other input means such as a data reader, or it may be remotely transmitted from a remote source. The parametric model mathematically models the portion of hydraulic circuit between the sensors 109 and 109a and presents a prefixed number of free parameters that are determined as described herein below in order to characterize the parametric mathematical model into one single model. In practice, the parametric mathematical model defines a family of mathematical models and is univocally characterized only once the parameters of the model are determined.

(34) In addition to command the circulation of dialysis liquid in lines 19 and 13, the control unit 10, e.g. by acting one or more concentrate pumps 105, 108, causes an upstream variation of the value of the characteristic Cd.sub.in in the fresh treatment liquid with respect to said prescription baseline Cd.sub.set and then re-establishes the characteristic Cd.sub.in in the fresh treatment liquid to said prescription baseline Cd.sub.set. Note that the alteration of the characteristic may be made using any means able to momentarily change the characteristic of the dialysis liquid, e.g. the conductivity or the concentration for one or more substances in the fresh dialysis fluid: for instance, a bolus pump configured to inject a predefined bolus of saline may be used for this purpose. The upstream variation causes a corresponding and timely delayed downstream variation of the same characteristic Cd.sub.out in the spent liquid flowing in the spent dialysate line: FIG. 4 schematically shows the time delay T.sub.D between the upstream variation and the downstream variation; the time delay which is also referred to as hydraulic delay depends upon the components such as tubing and second chamber interposed between the sensor 109 and the sensor 109a. The time delay T.sub.D between the upstream variation and the downstream variation is also shown in FIGS. 5, 6, 7, 8.

(35) In order to determine the parameters of the parametric mathematical model, the control unit 10 is configured to receive, e.g. from sensor 109a, measures of a plurality of values taken by a reference portion 200 of the downstream variation of the characteristic Cd.sub.out in the spent liquid. The measured values taken by the reference portion 200 of the variation in the characteristic Cd.sub.out may be measured by first identifying the initiation of a ramp-up or of a ramp-down portion of the downstream variation with respect to a respective baseline value of the same characteristic Cd.sub.out in the spent liquid, and then by measuring the plurality of values, as values taken by said ramp-up portion or ramp-down portion of said downstream variation. According to an aspect of the invention, the reference portion 200 which is used by the control unit 10 to characterize the mathematical model has a duration T.sub.R significantly shorter than the entire duration T.sub.v of the downstream variation: duration T.sub.R may be less than 70% and optionally less than 50% of duration T.sub.v. This is visible e.g. in FIGS. 5 and 6: FIG. 5 shows the duration T.sub.v of the entire downstream variation which certain conventional systems have to wait in order to calculate the effectiveness parameter, while FIG. 6 shows the much shorter interval T.sub.R necessary to characterize the mathematical model and then calculate the effectiveness parameter. More in detail, according to an aspect of the invention, the control unit 10 characterizes the mathematical model without having to wait for the entire interval T.sub.v by estimating the free parameters of the parametric mathematical model using measured values taken by the reference portion thereby identifying one single characteristic mathematical model using measured values taken during time interval T.sub.R which is much shorter than T.sub.v.

(36) Once the parameters of the model have been determined, the control unit 10 has the characteristic mathematical model and may compute the value of the effectiveness parameter supplying as input to the characteristic mathematical model one or more values taken by the characteristic Cd.sub.in in the fresh treatment liquid. In other words with use of the parametric mathematical model and with the characterization of the same by means of measured values taken by the characteristic Cd.sub.out during T.sub.R, it is possible to then calculate the effectiveness parameter with no need to take measures during the entire downstream variation, thus shortening the time during which control of the characteristic (e.g. concentration or conductivity of the dialysis liquid) should not be taken over by procedures other than the task for the determination of the effectiveness parameter. In other words the task for determining the effectiveness parameter should prevent execution of other procedures acting on the characteristic of the fresh dialysis liquid only until the end measurement instant T.sub.END.sub._.sub.MEAS represented in FIGS. 6, 8 and 9 which is the instant at which the measurement of said plurality of values of the reference portion of said downstream variation necessary for characterizing the model has been completed.

(37) In order to calculate the effectiveness parameter, the control unit 10 may for instance first compute at least one significant value of said downstream variation of the characteristic Cd.sub.out: the significant value of the downstream variation is a computed not measured value which, as shown in the example of FIG. 6, relates to a time subsequent to the duration of the reference portion, for instance it may represent an asymptotic value Cd.sub.out2 that the downstream variation would take after a relatively long time. This value is computed by using the characteristic mathematical model providing one or more real or set values representative of the upstream variation; once the significant value Cd.sub.out2 has been determined, the control unit 10 may compute at least one value of a parameter (D, Cb.sub.in, K, K.Math.t/V) indicative of the effectiveness of the extracorporeal blood treatment from said computed significant value and from one or more values taken by the characteristic Cd.sub.in in the fresh treatment liquid.

(38) The computation of the at least one significant value or directly of the effectiveness parameter comprises determining the value Cd.sub.out(n) of characteristic Cd.sub.out in the spent liquid at time instant (n) by using as input to the mathematical model: a) the measured values of characteristic Cd.sub.in in the fresh treatment liquid at a plurality of time instants (n1, n2, n3) preceding in time the time instant (n), as measured for instance by sensor 109; or b) a mathematically calculated version of characteristic Cd.sub.in in the fresh treatment liquid; in this second case the input is a set curve or a number of set values which are fed as input to the mathematical model.

(39) The mathematical modelfor instance a time invariant linear (LTI) modelmay be represented in the time domain by the following recursive equation:
y(n)=a.sub.0.Math.u(n)+b.sub.1.Math.y(n1)+b.sub.2.Math.y(n2)+ . . . b.sub.m.Math.y(nm),
Thus, the value Cd.sub.out(n) of characteristic Cd.sub.out in the spent liquid at time instant (n) subsequent to said reference portion is calculated with the following recursive equation:
Cd.sub.out(n)=a.sub.0.Math.Cd.sub.in(n)+b.sub.1.Math.Cd.sub.out(n1)+b.sub.2.Math.Cd.sub.out(n2)+ . . . b.sub.m.Math.Cd.sub.out(nm),
wherein:
Cd.sub.out(n) is the calculated value of the outlet characteristic at time instant (n),
Cd.sub.in(n) is the known value of the inlet characteristic at time instant (n),
Cd.sub.out(n1), Cd.sub.out(n2), Cd.sub.out(nm) are values of the outlet characteristic at preceding time instants (n1, n2, . . . nm) prior to time instant (n) and recursively computed by means of the mathematical model.
a.sub.0, b.sub.1, b.sub.2, . . . , b.sub.m are constant parameters that characterize the mathematical model, as estimated by using said measured values of the reference portion of the downstream variation.

(40) In the frequency domain and using the z-Transformthe mathematical model is described by a transfer function H(z) having at least one zero and at least one pole. In an embodiment, the transfer function H(z) comprises a plurality of poles, e.g. from 2 to 5 poles, and is described by one of the following:
H(z)=Cd.sub.out(z)/Cd.sub.in(z)=a.sub.0/(1b.sub.1.Math.z.sup.1b.sub.2.Math.z.sup.2b.sub.3.Math.z.sup.3b.sub.4.Math.z.sup.4b.sub.5.Math.z.sup.5),
H(z)=Cd.sub.out(z)/Cd.sub.in(z)=a.sub.0/(1b.sub.1.Math.z.sup.1b.sub.2.Math.z.sup.2b.sub.3.Math.z.sup.3b.sub.4.Math.z.sup.4),
H(z)=Cd.sub.out(z)/Cd.sub.in(z)=a.sub.0/(1b.sub.1.Math.z.sup.1b.sub.2.Math.z.sup.2b.sub.3.Math.z.sup.3),
H(z)=Cd.sub.out(z)/Cd.sub.in(z)=a.sub.0/(1b.sub.1.Math.z.sup.1b.sub.2.Math.z.sup.2),
wherein
a.sub.0, b.sub.1, b.sub.2, b.sub.3, b.sub.4, b.sub.5 are constant parameters of the model, as estimated by using said measured values of the reference portion of the downstream variation.

(41) FIGS. 5,6 and 7,8 respectively show two possible implementations of the invention. In FIGS. 5 and 6 the characteristic is altered from a first to a second value and kept at the second value for a relatively long time, while in FIGS. 7 and 8 the characteristic is kept at the second value for a relatively short time. More in detail, in the example of FIGS. 5 and 6, the value of the characteristic Cd.sub.in in the fresh treatment liquid is varied by imposing a change of the same from a first inlet value Cd.sub.in1 to a second inlet value Cd.sub.in2, which may be kept constant for a prefixed time interval of e.g. 3 to 10 minutes, thereby causing a corresponding change of the characteristic Cd.sub.out in spent liquid from a respective first outlet value Cd.sub.out1 to a respective second outlet value Cd.sub.out2 defining said timely delayed downstream variation of the characteristic Cd.sub.out. In the example of FIGS. 5 and 6, the reference portion of the downstream variation begins after the characteristic in the spent liquid changes from said first outlet value Cd.sub.out1 and lasts a periodfor instance prefixed period T.sub.Rduring which the characteristic either continuously increases or decreases without reaching the second outlet value Cd.sub.out2. In the example shown, during T.sub.R the characteristic Cd.sub.out does not reach a prefixed fraction, e.g. 80%, of the second outlet value Cd.sub.out2. Moreover, there is no need to wait until the real value Cd.sub.out2 is actually reached. Instead, the second outlet value Cd.sub.out2 of the characteristic Cd.sub.out is calculated by using as input to the characteristic mathematical model the values of characteristic Cd.sub.in in the fresh treatment liquid, or a mathematically calculated version of the characteristic Cd.sub.in in the fresh treatment liquid.

(42) Then, the calculated second outlet value Cd.sub.out2 is used as significant value for the computation of at least one value of a parameter (D, Cb.sub.m, K, K.Math.t/V) indicative of the effectiveness of the extracorporeal blood treatment. In accordance with an aspect, if the parameter comprises is effective dialysance D, each computed value D.sub.k of the dialysance each respective variation is obtained using the formula:
D.sub.k=(Qd+WLR).Math.[1(Cd.sub.out2Cd.sub.out1)]/(Cd.sub.in2Cd.sub.in1)
where:
Cd.sub.out1 is the first outlet value taken by the characteristic in the spent dialysate line downstream of the secondary chamber in response to the change of characteristic Cd.sub.in in the preparation line to said first inlet value Cd.sub.in1,
Cd.sub.out2 is the calculated second value (namely the significant value) which is representative of the value taken by the characteristic in the spent dialysate line downstream of the secondary chamber in response to the change of characteristic Cd.sub.in in the preparation line from said first inlet value Cd.sub.in1 to said second inlet value Cd.sub.in2,
Cd.sub.in2 are first and second inlet values taken by the characteristic (Cd) in the preparation line upstream of the secondary chamber,
Qd is the fresh treatment liquid flow rate in the preparation line,
WLR is the weight loss rate of a patient under treatment.

(43) In FIGS. 7, 8 the upstream perturbation is shorter and the value of the characteristic Cd.sub.in in the fresh treatment liquid is varied by imposing a change of the same from a first inlet value Cd.sub.in1 to a second inlet value Cd.sub.in2, which may optionally be kept constant for a prefixed time interval of e.g. 1 to 2 minutes, and then a further change to a third inlet value Cd.sub.out3 thereby causing a corresponding change of the characteristic Cd.sub.out in spent liquid from a respective first outlet value Cd.sub.out1 to a respective second outlet value Cd.sub.out2 and then back to a third value Cd.sub.out3 to define said timely delayed downstream variation of the characteristic Cd.sub.out. In the example of FIGS. 7 and 8, Cd.sub.in1 is equal or close to Cd.sub.in3. In case a short perturbation is used, the formula needed for the calculation of the effectiveness parameter requires more than simply the knowledge of one significant value such as Cd.sub.out2. In the example of FIGS. 7 and 8, the reference portion of the downstream variation begins after the characteristic in the spent liquid changes from said first outlet value Cd.sub.out1 and lasts a periodfor instance prefixed period T.sub.Rduring which the characteristic either continuously increases or decreases. During T.sub.R the characteristic Cd.sub.out may or may not reach the second outlet value Cd.sub.out2. In accordance with an aspect, there is no need to wait until Cd.sub.out reaches Cd.sub.out2 and returns to the baseline value Cd.sub.out3. Instead, the second and third Cd.sub.out2 and Cd.sub.out3 or at least the third outlet value Cd.sub.out3 of the characteristic Cd.sub.out are/is calculated by using as input to the characteristic mathematical model the values of characteristic Cd.sub.in in the fresh treatment liquid, or a mathematically calculated version of the characteristic Cd.sub.in in the fresh treatment liquid.

(44) Once the values Cd.sub.out1, Cd.sub.out2, Cd.sub.out3 have been calculated, the effectiveness parameter may be determined based on these calculated values and on one or more inlet values of the conductivity, e.g. Cd.sub.in1, Cd.sub.in2, Cd.sub.in3. For instance if dialysance is to be calculated, the following formula may be adopted:
D=(Qd+WLR)[1(2Cd.sub.out1Cd.sub.out2Cd.sub.out3)/(2Cd.sub.in1Cd.sub.in2Cd.sub.3)]
where:
Cd.sub.out1 is the first outlet value taken by the characteristic in the spent dialysate line downstream of the secondary chamber in response to the change of characteristic Cd.sub.in in the preparation line to said first inlet value Cd.sub.in1,
Cd.sub.out2 is the calculated second value (namely one of the significant values) which is representative of the value taken by the characteristic in the spent dialysate line downstream of the secondary chamber in response to the change of characteristic Cd.sub.in in the preparation line from said first inlet value Cd.sub.in1 to said second inlet value Cd.sub.in2,
Cd.sub.out3 is the calculated third value (namely one of the significant values) which is representative of the value taken by the characteristic in the spent dialysate line downstream of the secondary chamber in response to the change of characteristic Cd.sub.in in the preparation line from said second inlet value Cd.sub.in2 to said third inlet value Cd.sub.in3,
Cd.sub.in1, Cd.sub.in2, Cd.sub.in3 are first, second and third inlet values taken by the characteristic (Cd) in the preparation line upstream of the secondary chamber,
Qd is the fresh treatment liquid flow rate in the preparation line,
WLR is the weight loss rate of a patient under treatment.

(45) According to a further embodiment, see FIG. 9, varying the value of the characteristic Cd.sub.in in the fresh treatment liquid comprises imposing an upstream perturbation, which may be in the shape of a sinusoid or of a short peak, in the characteristic of the fresh treatment liquid thereby causing a corresponding downstream perturbation of the characteristic Cd.sub.out in spent liquid. The reference portion of said downstream perturbation begins after the characteristic in the spent liquid changes from said first outlet value Cd.sub.out1 and lasts a prefixed period shorter than a fraction, e.g. 60% or even 50%, of the duration of the entire downstream perturbation. The control unit 10 determines in this case a plurality, e.g. 10 or more, of significant values of the characteristic Cd.sub.out, describing a remaining portion of the downstream perturbation consecutive to said reference portion, by using as input to the mathematical model the values of characteristic Cd.sub.in in the fresh treatment liquid, or a mathematically calculated version of the characteristic Cd.sub.in in the fresh treatment liquid, thereby obtaining a calculated downstream perturbation from said extrapolated significant values. Then, using e.g. the formulas described in EP 0920877, the control unit computes at least one value of a parameter (D, Cb.sub.m, K, K.Math.t/V) indicative of the effectiveness of the extracorporeal blood treatment by comparing the calculated downstream perturbation and the upstream perturbation.

(46) In accordance with a further aspect of the invention, the control unit 10 may also be configured to determine the baseline of the downstream curve representative of the values Cd.sub.out(t) taken over time by said characteristic in the spent dialysate line downstream of the secondary chamber. The baseline of the downstream curve Cd.sub.out(t) may be determined using measured values of the characteristic Cd.sub.out in the spent liquid or using a calculated curve representative of the downstream variation which has been previously determined using the characteristic mathematical model. In this second option only measured values of the characteristic Cd.sub.out in the spent liquid during said reference portion are used for the determination of the free parameters to identify the characteristic mathematical model; then using said identified characteristic mathematical model, a downstream curve Cd.sub.out(t) representative of the values taken by the characteristic Cd.sub.out in the spent liquid is mathematically determined and the baseline thereof identified.

(47) The control unit may be configured to determine an angular deviation between the baseline of the downstream curve Cd.sub.out(t) with respect to the prescription baseline Cd.sub.set, and to compensate for said angular deviation by angularly rotating the downstream curve to obtain a corrected downstream curve Cd.sub.out-correct(t), as shown in a the enlarged representation of FIG. 10.

(48) According to a yet further aspect, the control unit 10 is configured to remove undesired noise and perturbation from the characteristic Cd.sub.out. In accordance with an aspect, the control unit may receive measured values of the characteristic Cd.sub.out in the spent liquid during said reference portion, estimate the free parameters of the parametric mathematical model to identify the characteristic mathematical model, determine a downstream curve Cd.sub.out(t) representative of the values taken by the characteristic Cd.sub.out in the spent liquid using said identified characteristic mathematical model, analyze a frequency spectrum of the downstream curve Cd.sub.out(t), filter out harmonics of said frequency spectrum of the downstream curve Cd.sub.out(t) lying at frequencies higher than a prefixed threshold to eliminate noise and undesired perturbations possibly present in the downstream curve and obtain a corrected downstream curve Cd.sub.out-correct(t).

(49) Although the above description referred to one single parametric mathematical model, the control unit 10 may further be configured for storing a plurality of mathematical models each of which puts into relation the characteristic (Cd.sub.in) in the fresh treatment liquid with the characteristic (Cd.sub.out) in the spent liquid. In this case the control unit may be configured for selecting the mathematical model to be used for computing the at least one significant value of said downstream variation based on certain factors such as for instance: the shape of the upstream variation (one mathematical model may be better suited for a long step perturbation while another model may more properly operate for a short sinusoidal change), the type of blood treatment unit used by the apparatus, whether or not particular hydraulic components are present in the circuit section between sensor 109 and sensor 109a.

(50) Aspects of the invention are also disclosed in FIG. 11, which shows a flowchart exemplifying a method for determining an effectiveness parameter. The method may be executed by the control unit 10 of any one of the apparatuses disclosed herein above or claimed in the appended claims.

(51) The method comprises the following steps. step 300: selection of the mathematical model; step 301: measurement of values of conductivity, or concentration, Cd.sub.out in the spent dialysate corresponding to a variation respectively in the conductivity, or in the concentration of at least one substance, Cd.sub.in made on the fresh dialysis liquid flowing upstream the blood treatment unit (step 301); the measures are taken during the reference time T.sub.R which is sensibly shorter than the duration of the downstream variation, as already explained herein above; step 302: characterization of mathematical model using the measured value(s) of Cd.sub.out taken during the reference time T.sub.R and identification of a single mathematical model; step 303: determination, using the mathematical model, of significant value(s) necessary for the calculation of the effectiveness parameter; the significant values may be one or more calculated conductivity or concentration values of the downstream variation at instants following the reference period (such as Cd.sub.out2 or Cd.sub.out2 and Cd.sub.out3) step 304: determination of effectiveness parameter using the calculated significant value or values. The calculation of the effectiveness parameter may be made using any one of the formulas described above.

Example

(52) Here below an example is described, with reference to FIGS. 12-15, showing use of a one-zero and three-pole mathematical model to mathematically calculate the entire downstream variation; it is relevant noticing that to characterize the model, only measured values relative to a reference portion of the downstream variation having relatively short duration compared to the duration of the entire downstream variation are used. The Example provided adopts an exemplifying mathematical model and makes reference to a step perturbation imposed in the liquid flowing upstream the blood treatment unit. Of course other mathematical models may be adopted and the upstream perturbation may be different from a step-shaped perturbation. Furthermore, the example makes reference to conductivity variations and corresponding measures: of course the same procedure may be adopted using variations, and corresponding measures, in the concentration of at least one substance in the dialysis liquid.

(53) Referring now to the diagram of FIG. 12, two curves are represented: a first curve represents the inlet conductivity Cd.sub.in and shows a step-shaped conductivity variation (1 mS/cm) that has been imposed to the conductivity of the dialysis liquid flowing upstream the blood treatment unit, while a second curve (below the first curve) represents the downstream conductivity Cd.sub.out and shows the corresponding variation in the conductivity of the spent dialysis liquid as a consequence of the step-shaped perturbation on the upstream conductivity Cd.sub.in. FIG. 13 is an enlarged view of FIG. 12 and focuses on the outlet conductivity: notice that the curve in FIG. 13 is obtained measuring the outlet conductivity values from time 700 s to time 950 s (i.e. 250 seconds). FIG. 13 shows the value of the outlet conductivity Cd.sub.out2 at time 910 which is regarded as the significant value of interest, necessary for the calculation of e.g. dialysance when using formula:
D.sub.k=(Qd+WLR).Math.[1(Cd.sub.out2Cd.sub.out1)]/(Cd.sub.in2Cd.sub.in1)

(54) According to one aspect of the invention, instead of measuring the conductivity values until time 950 s, measures are taken only during reference portion T.sub.R (please refer to FIG. 14) i.e. for the 100 seconds only

(55) Then, using the following a one-zero and three-pole model:

(56) $H (z) = \frac{a_{0}}{1 - b_{1} z^{- 1} - b_{2} z^{- 2} - b_{3} z^{- 3}}$

(57) The following parameters are estimated using the measured values of Cd.sub.out during reference portion T.sub.R: a0=0.004209932871 b1=2.905495405197 b2=2.815777778625 b3=0.910210132599
giving

(58) $H (z) = \frac{0.004209932871}{\begin{matrix} 1 - 2.905495405197 z^{- 1} + 2.815777778625 z^{- 2} - \\ 0.910210132599 z^{- 3} \end{matrix}}$

(59) By feeding an idealized unit step (i.e. a calculated step) of appropriate length (e.g. 200 to 300 s) to this model and by suitably adding the baseline value Cd.sub.out1 to the model output, we get a signal as shown in FIG. 15 (continuous line represents model output, while dotted line schematically represents measured Cd.sub.out), which closely approximates the behavior of the system in a time interval following the reference portion.

(60) The following table shows the measured versus computed values of Cd.sub.out in the neighborhood of time n=910 where the good match between measured and computed values can be seen.

(61) TABLE-US-00001 Time Cd.sub.out model (mS .Math. 100/cm) Cd.sub.out measured (mS .Math. 100/cm) 905 1358.567173 1359 906 1358.586262 1359 907 1358.604656 1359 908 1358.622398 1359 909 1358.639698 1359 910 1358.656872 1359 911 1358.674126 1359 912 1358.691517 1359 913 1358.709053 1359 914 1358.726700 1359 915 1358.744549 1359

(62) The calculated significant value Cd.sub.out2 at time 910 is 13.59 mS/cm is very close to the corresponding measured value (13.58656872 mS/cm). Thus, the dialysance calculation using the above formula and relying on the calculated value Cd.sub.out2 of 13.59 mS/cm will provide exactly the same result as when using a measured valued for Cd.sub.out2, while requiring actual measurements only during T.sub.R.

(63) Control Unit 10

(64) As already indicated the apparatus according to the invention makes use of at least one control unit 10. This control unit 10 may comprise a digital processor (CPU) with memory (or memories), an analogical type circuit, or a combination of one or more digital processing units with one or more analogical processing circuits. In the present description and in the claims it is indicated that the control unit 10 is configured or programmed to execute certain steps: this may be achieved in practice by any means which allow configuring or programming the control unit 10. For instance, in case of a control unit 10 comprising one or more CPUs, one or more programs are stored in an appropriate memory: the program or programs containing instructions which, when executed by the control unit 10, cause the control unit 10 to execute the steps described and/or claimed in connection with the control unit 10. Alternatively, if the control unit 10 is of an analogical type, then the circuitry of the control unit 10 is designed to include circuitry configured, in use, to process electric signals such as to execute the control unit 10 steps herein disclosed.

(65) While the invention has been described in connection with what is presently considered to be the most practical and preferred embodiments, it is to be understood that the invention is not to be limited to the disclosed embodiments, but on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and the scope of the appended claims.

Apparatus for determining a parameter indicative of the progress of an extracorporeal blood treatment

Assignee

Inventors

Cpc classification

Classification Explorer

A61M1/16

HUMAN NECESSITIES

Classification Explorer

A61M1/1617

HUMAN NECESSITIES

Classification Explorer

A61M2205/3393

HUMAN NECESSITIES

Classification Explorer

G05B13/04

PHYSICS

Classification Explorer

A61M2209/02

HUMAN NECESSITIES

Classification Explorer

G05D7/0617

PHYSICS

Classification Explorer

A61M2205/3317

HUMAN NECESSITIES

Classification Explorer

A61M1/367

HUMAN NECESSITIES

Classification Explorer

A61M1/1609

HUMAN NECESSITIES

Classification Explorer

A61M2205/3324

HUMAN NECESSITIES

Classification Explorer

A61M1/1613

HUMAN NECESSITIES

Classification Explorer

G06F17/10

PHYSICS

Classification Explorer

A61M1/1607

HUMAN NECESSITIES

Classification Explorer

A61M2205/3334

HUMAN NECESSITIES

International classification

Classification Explorer

A61M37/00

HUMAN NECESSITIES

Classification Explorer

G06F17/10

PHYSICS

Classification Explorer

G05D7/06

PHYSICS

Classification Explorer

G05B13/04

PHYSICS

Classification Explorer

A61M1/16

HUMAN NECESSITIES

Classification Explorer

A61M1/36

HUMAN NECESSITIES

Abstract

Claims

Description