TH E BACKGROUN D OF TH E I NVENTION AN D P RIOR ART

The present invention relates generally to automatic determination of at least one attribute of a liquid substance. More particu- larly the invention relates to a monitoring unit according to the preamble of claim 1 and a method according to the preamble of claim 1 2. The invention also relates to a computer program according to claim 21 and a computer readable medium according to clai m 22. There are many areas where efficient and reliable handling of liquid substances is important, for instance in connection with production and processing of foodstuffs. Moreover, automatic milking solutions are becoming increasingly efficient and sophisticated. Today, there is also a strong demand for flexible and animal- friendly milk production. For example, so-called milking robots have been introduced, which enable animals to autonomously decide when they are to be milked. Milking robots are advantageous because they reduce the amount of manual labor involved in milking. Milking robots are also desirable from an animal health point-of-view, since thereby it is uncomplicated to extract milk more frequently than by applying the existing alternative solutions, and in general, high-frequency milking vouches for a good udder health. However, milking robots may be somewhat problematic because these machines are often operated without any human operator being present. This, in turn, renders a safe and reliable operation highly important. One aspect of such an operation is that problems related to unsatisfying milk quality must be resolved automatically. For example it is necessary to automatically detect unacceptably high concentrations certain cons- tituents in the milk, such that adequate measures can be taken promptly. US 4,678,995 describes an apparatus and method for determining the presence of specific substances in a sample, e.g. of tissue, by nuclear magnetic resonance (NM R) and then obtaining a high-resolution NM R image of these substances. US 7, 1 41 ,978 discloses a solution for analyzing fluid samples by means of an NM R probe comprising multiple N M R detection sites.

WO 00/39578 reveals a method and apparatus for estimation of a cell count in a body fluid. Here, a mid-I R spectrum is recorded and the spectral information therein serves as a basis for deriving a number of cells in the body fluid . This document also briefly mentions the possibility to employ NM R spectroscopy.

US 2009/0278685 describes methods and systems for calibration of radio-frequency identification (RFI D) sensors involving impedance measuring by determining a complex impedance spectrum, a phase angle and/or a magnitude of the impedance. Inter alia, the document mentions that the RFI D sensors may be adapted to measure physical , chemical and biological parameters; each sensor can have a digital I D and be calibrated to ac- curately react to a parameter of interest through changes in measurements of the sensor's complex impedance.

US 7,21 9,024 discloses a system, method and program product for determining in-place engineering properties, such as density and moisture content of certain engineering materials. A data- base, material model and sensor model are also shown. The document generally relates to material analysis, and to the field of impedance spectroscopy, and the determination of engineering properties of a material from the response to electromagnetic probing in a defined frequency spectru m. US 6,51 1 ,851 reveals a method for identifying a change in the composition of a liquid . Here, a time varying electrical or elect- romagnetical input signal is applied to the liquid in a range of frequencies encompassing a resonant frequency of an electrical circuit comprising the liquid . An impedance quantity is measured of said electrical circuit by means of the output signal as a function of the frequency of the time varying input signal in said range of frequencies. A resonant frequency of said electrical cir- cuit is determined. After a change in the composition of the liquid, variation in the impedance quantity is measured at or near the previously determined resonant frequency of said electrical circuit. The variation in the impedance quantity is then related to the change in the composition . WO 2008/076453 describes a flex fuel sensor, which is deployed in conjunction with a fuel transfer line, or at the bottom/side of a fuel tank. A radio-frequency signal at a constant frequency may be generated across a resonant circuit, which comprises an inductor and a printed circuit board trace capacitor, capacitor pla- tes, semi cylindrical capacitor plates, or the like. Electromagnetic radiation is propagated into the passing fuel in the transfer pipe. The conductivity and dielectric properties of the fuel change the capacitance of the trace capacitor plates. These changes are proportional to the ethanol/alcohol content of the fuel , and are preferably detected by a microcontroller, or the like, and then transmitted to a flex fuel vehicle engine management system .

PROBLEMS ASSOCIATE D WITH TH E PI ROR ART

Hence, diverse solutions are known for determining various pro- perties of different kinds of entities (e.g. represented by tissue pieces, fluid samples or RFI D sensors) . However, there is yet no straightforward solution for unobtrusive online monitoring of the characteristics a flowing liquid.

SUM MARY OF TH E I NVENTION

The object of the present invention is to alleviate the above problem, and thus offer an efficient solution for determining a measure representing an attribute of a flowing liquid substance, for instance in the form of a liquid food product.

According to one aspect of the invention, the object is achieved by the initially described monitoring unit, wherein the unit includes a dielectric conduit, a measuring circuitry and a processor. The dielectric conduit is configured to transport a flow of the liquid substance through the monitoring unit. The measuring circuitry includes a coil member encircling the dielectric conduit. The measuring circuitry is further configured to register an impedance measure of the liquid substance when flowing through the dielectric conduit. The processor is configured to derive the attribute of the liquid substance based on the registered impedance measure.

This monitoring unit is advantageous because in addition to allowing continuous monitoring in real time, the unit is straightfor- ward to install , uncomplicated to operate, and since no consumables are required, essentially maintenance free.

According to one preferred embodiment of this aspect of the invention, the determined attribute contains a first number of different parameters and the first number is larger than or equal to two. The monitoring unit further includes at least one dielectric conduit and a second number of measuring circuitries. The at least one dielectric conduit is configured to transport the flow of the liquid substance through the monitoring unit. The second nu mber, in turn, is larger than or equal to the first nu mber. Further, each measuring circuitry includes a coil member encircling one of the at least one dielectric conduit, and each measuring circuitry is configured to register a respective impedance measure of the liquid substance when flowing through said one of the at least one dielectric conduit in response an electro- magnetic field applied to the liquid substance via the at least one coil member. Here, the processor is configured to derive each of the first number of parameters based on the registered impedance measures. Such a design is advantageous because it provides convenient monitoring of multiple constituents and/or characteristics of the liquid substance.

According to another preferred embodiment of this aspect of the invention, the coil member in each measuring circuitry is configured to apply an electromagnetic field having a particular spectral range. The total number of different spectral ranges of electromagnetic energy applied to the liquid substance via the coil members is here larger than or equal to the first number. Preferably, each electromagnetic field has a relatively narrow spectral range with a distinct center frequency located in an in- terval from 1 M Hz to 1 500 M Hz. Thus, for example one constituent/characteristic per center frequency of the liquid substance can be determined.

According to yet another preferred embodiment of this aspect of the invention, the processor is associated with a data bank of coefficients describing relationships between impedance measures of the liquid substance and the different parameters. The processor is here configured to derive said parameters by applying analytic calculations and/or numerical methods to equations describing relationships between the impedance measures and the different parameters. Thereby, a comparatively large nu mber of constituents and/or characteristics of the liquid substance can be determined in a relatively uncomplicated manner.

According to still another preferred embodi ment of this aspect of the invention, the processor is associated with a lookup table describing relationships between impedance measures of the liquid substance and the different parameters of said attribute, and the processor is configured to derive said parameters from said lookup table, either directly or by interpolating between values therein. Hence, a comparatively large number of consti- tuents and/or characteristics of the liquid substance can be determined in a very processing efficient manner.

According to a further preferred embodi ment of this aspect of the invention, the determined attribute of the liquid substance in- eludes one or more of the parameters: water content, a concentration of sodium, the total concentration of solids, pH level and electrical conductivity. Consequently, the proposed monitoring unit is generally well suited for online checking of liquid food pro- ducts.

According to another preferred embodiment of this aspect of the invention, the liquid substance is presumed to be milk, and the determined attribute includes one or more of the parameters: a concentration of lactose, a concentration of fat, a concentration of protein , a concentration of urea and a concentration of somatic cells. Further preferably, parameters reflecting a concentration of macrophages, leucocytes and/or polymorphonuclear leukocytes may also be determined . Hence, the monitoring unit is especially adapted for employment in conjunction with a milking installation, for example to prevent mixing of poor-quality milk with good milk.

According to still a further preferred embodiment of this aspect of the invention , each measuring circuitry includes a power source and an interface. The power source is configured to ge- nerate electric energy to the coil member, such that in response thereto the coil member produces an electromagnetic field having a particular spectral range. The interface is configured to enable the processor to register the impedance measure of the measuring circuitry. This renders the overall design exceptional- ly uncomplicated.

According to yet another preferred embodiment of this aspect of the invention, the monitoring unit includes at least two parallel dielectric conduits, each of which is configured to transport a fraction of the flow of liquid substance. Each dielectric conduit is also encircled by at least one coil member of a respective one of said measuring circuitries. Thus, the monitoring unit can be made physically compact.

According to another aspect of the invention, the object is achie- ved by the method described initially, wherein a flow of the liquid substance is transported through a dielectric conduit. When the liquid substance is flowing there through , the liquid substance is subjected to the electromagnetic field, an impe- dance measure is registered, and based thereon said attribute is derived . The advantages of this method, as well as the preferred embodiments thereof, are apparent from the discussion above with reference to the proposed monitoring unit.

According to a further aspect of the invention the object is achieved by a computer program, which is directly loadable into the memory of a computer, and includes software adapted to implement the method proposed above when said program is run on a computer.

According to another aspect of the invention the object is achie- ved by a computer readable medium , having a program recorded thereon, where the program is to control a computer to perform the method proposed above when the program is loaded into the computer.

Further advantages, beneficial features and applications of the present invention will be apparent from the following description and the dependent claims.

BRI EF DESCRI PTION OF TH E D RAWI NGS

The invention is now to be explained more closely by means of preferred embodiments, which are disclosed as examples, and with reference to the attached drawings.

Figure 1 shows a perspective view of a monitoring unit according to a first embodiment of the invention ;

Figure 2 shows a schematic side view of a measuring circuitry according to one embodiment of the invention ;

Figure 3 shows graphs representing a respective measured voltage as a function of frequency for each of five different milk samples;

Figure 4 schematically illustrates a monitoring unit according to a second embodiment of the invention ; and

Figure 5 illustrates, by means of a flow diagram , the general method according to the invention .

DESCRI PTION OF PREFERRE D EMBODI M ENTS OF TH E I NVENTION

Methods for detecting constituents in liquids, primary water based solutions, may use the change in water's permittivity (ε>80) at electromagnetic wave frequencies in the range 1 M Hz to 200 M Hz. These methods are strongly dependent on the concentrations of impurities in the water. Since knowing only the dielectric property is insufficient to draw a conclusion regarding the concentration levels of individual elements, the present invention proposes that dielectric spectroscopy (DS) be applied.

In general , DS (which may also be referred to as impedance spectroscopy) measures the dielectric properties of a medium as a function of frequency. DS is based on the interaction of an external field with the electric dipole moment of a liquid sample. This interaction is often expressed by permittivity.

There are a number of different dielectric mechanisms connected to the way a studied medium reacts to an applied electromagne- tic field. Each dielectric mechanism is centered on its characteristic frequency, which is the reciprocal of the characteristic time of the process. Normally, the dielectric mechanisms are divided into relaxation processes and resonance processes, and the most common dielectric mechanisms are: dipole relaxation (in the M Hz-to-GHz range) , ionic relaxation (in the kHz range) and dielectric relaxation (in the kHz-to-GHz range) . The working frequencies for the present invention are in the M Hz range, and therefore dipole relaxation and dielectric relaxation are primarily involved.

Dipole relaxation originates from permanent and induced dipoles aligning to an electric field. The orientation of the dipoles' pola- rization is disturbed by thermal noise (which misaligns the dipole vectors from the overall direction of the field) , and the time needed for the dipoles to relax is determined by the local viscosity of the liquid medium. This renders dipole relaxation heavily dependent on temperature and chemical environment. The pola- rization of a dielectric material is a competition between torques due to the imposed electric field, which tends to align the molecules; and collisions, which tend to destroy the alignment.

Dielectric relaxation as a whole is the result of the movement of dipoles, i.e. dipole relaxation. Dielectric relaxation is the momen- tary delay, or lag , in the dielectric constant of a material . This is usually caused by the delay in the molecular polarization with respect to a changing electric field in a dielectric medium.

Ionic relaxation as a whole is the result of electric charges moving due to an applied alternating field. The time lag between electrical field and polarization implies an irreversible degradation of free energy. In physics, dielectric relaxation typically refers to the relaxation response of a dielectric medium to an external electric field of microwave frequencies. Such a relaxation is often described in terms of permittivity as a function of frequency.

The theory of DS gives us an understanding that an electromag- netical field can influence different molecules, and even cells, at a variety of frequencies being unique to each type of molecule or cell . This phenomenon allows quantifying a majority of the com- ponents and contaminants in a liquid substance.

We now refer to Figure 1 , which shows a perspective view of a monitoring unit 1 00 for determining an attribute P of a liquid substance according to a first embodiment of the invention. The monitoring unit 1 00 includes measurement means, a dielectric conduit 1 1 0 and a processor 1 30. The measurement means, which here include coil members 1 21 , 1 22 and 1 23, are configured to subject a liquid substance flowing F through the dielectric conduit 1 1 0 to an electromagnetic field. Under influence of the electromagnetic field, the measurement means are further configured to measure at least one electromagnetic property of the liquid substance, namely register impedance measures z^ ) , z(f ₂ ) and z(f ₃ ) respectively. To achieve the electromagnetic influence and enable said measurements, each of the coil members 1 21 , 1 22 and 1 23 includes at least one electric conductor loop encircling the dielectric conduit 1 1 0. Preferably, to avoid interference the coil members 1 21 , 1 22 and 1 23 are also electromagnetically screened from one another. The processor 1 30 is configured to derive the above-mentioned attribute P based on the registered i mpedance measures z^ ) , z(f ₂ ) and z(f ₃ ) . If the attribute P only includes one parameter it is sufficient that the measurement means contains a single measuring circuitry. Figure 2 shows an example of such a design, which will be used to explain operating principle of the monitoring unit according to the invention .

In Figure 2, a schematic measuring circuitry and dielectric conduit 1 1 0 are shown in a side view. As can be seen the coil member 1 20 here contains six loops of an electric conductor around the dielectric conduit 1 1 0.

A power source 21 0 is connected to the coil member 1 20 and configured to deliver electric energy to the coil member 1 20. Consequently, in response to the electric energy, the coil member 1 20 produces an electromagnetic field B(f) , which has a par- ticular spectral range f and extends into the dielectric conduit 1 1 0. Depending on the characteristics of the liquid substance flowing F through dielectric conduit 1 1 0, the electromagnetic field B(f) will experience particular impedance. The measuring circuitry further includes an interface, here represented by sym- bolic resistive and capacitive components R and C respectively, towards the processor 1 30 through which interface the processor 1 30 is enabled to register an impedance measure z(f) reflecting the impedance experienced by the electromagnetic field B(f) . Thereby, the processor 1 30 may determine the attribute P of the liquid substance, as well as detect any variations in the characteristics relating to this attribute P .

Returning now to Figure 1 , we see that the attribute P contains a first number of different parameters denoted p _{1 ;} p ₂ and p ₃ res- pectively. The first number is larger than or equal to two. In this example, the first number is three.

The parameters p ₂ and p ₃ of the attribute P may represent various characteristics and/or constituents of the liquid substance, such as water content, a concentration of sodium , the total concentration of solids, pH level and electrical conductivity. Provided that the liquid substance is milk, the parameters p _{1 ;} p ₂ and p ₃ may further reflect a concentration of lactose, a concentration of fat, a concentration of protein , a concentration of urea and/or a concentration of somatic cells.

In any case, a second number of measuring circuitries are included in the monitoring unit 1 00, where the second number is larger than or equal to the first number. In other words, at least as many measuring circuitries are included as the number of parameters p _{1 ;} p ₂ and p ₃ that shall be determined . This relationship will be explained in further detail below. In the embodiment of Figure 2, the attribute P contains three parameters p, , p ₂ and p ₃ , and therefore three coil members 1 21 , 1 22 and 1 23 respectively encircle the dielectric conduit 1 1 0. I n response the electromagnetic fields applied to the liquid substance via the coil members 1 21 , 1 22 and 1 23, each measuring circuitry is configured to register a respective i mpedance measure z^ ) , z(f ₂ ) and z(f ₃ ) of the liquid substance when flowing F through the dielectric conduit 1 1 0. It is should be noted that there is typically not a direct relationship between a particular impedance measure z{ ), z(f ₂ ) or z(f ₃ ) and a given parameter p,, p ₂ or p ₃ . Instead, inter alia due to various interference and combinatory effects, the relations between the impedance measure z{\^), z(f ₂ ) and z(f ₃ ) and the parameters p ₂ and p ₃ are best described by an equation system. To avoid underdetermination, the number of impedance measures is larger than or equal to the number of parameters. The equation system, in turn, may be fairly complex and can for instance include non-linear components. Nevertheless, based on the registered impedance measures z^), z(f ₂ ) and z(f ₃ ) and said relations, the processor 130 is configured to derive each of the parameters p,, p ₂ and p ₃ .

To illustrate how the relations between the impedance measures z^), z(f ₂ ) and z(f ₃ ) and the parameters p ₂ and p ₃ may be determined we will now refer to a simple linear example. For d if - ferent concentrations of constituents in the liquid substance corresponding to each parameter p ₂ and p ₃ , a respective impedance measure is registered, for instance via impedance spectroscopy (IS). Here, the energy storage and energy dissipation of the liquid substance is probed over a range of frequencies, say from 1 MHz to 1500 MHz, by means of impedance measurements. These measurements may be carried out using a so-called potentiostat instrument. Since essentially all physio-chemical systems (e.g. solutions, organic molecules and biological tissue) possess such energy storage and dissipation properties it is possible to characterize a given physio-chemical system through its impedance response to a particular electromagnetic frequency. Consequently, before the monitoring unit 100 is used for measurement purposes, liquid substance is transported through the unit with various known concentrations p^ p ₂ and p ₃ of constituents therein, and for each concentration an impedance measure z(f i ), z(f ₂ ) and z(f ₃ ) is registered while the liquid substance is subjected to an electromagnetic field having a selected frequency (or spectral range) f ₂ and f ₃ respectively, which has been found to be particularly impedance sensitive with respect to the constituent in question. This can be described as:

or more compact: Z = K P .

Thus, the attribute P , representing the parameters p, , p ₂ and p ₃ , can be derived as:

P = Κ ^"1 ·Ζ which may also be expressed :

or in general terms:

1

(adj(K)) - Z .

det(K)

According to the invention, the total number of different spectral ranges of electromagnetic energy applied to the liquid substance via the coil members 1 21 , 1 22 and 1 23 is larger than or equal to the first number. Of course, for large numbers of parameters p, of the attribute P, fairly extensive calculations may be required. This is especially true if (in contrast to the above example) the relationships between parameters p, and the impedance measures Z j are non-linear. For relatively few parameters p, and uncomplicated relationships, however, a lookup table can be a useful tool. For somewhat more complex equation systems, analytic calculations may prove to be more efficient, whereas determining the parameters p, in especially intricate cases may require various numerical methods. According to one embodiment of the invention, the processor 130 is associated with a data bank of coefficients describing the relationships between impedance measures z{\^), z(f ₂ ) and z(f ₃ ) of the liquid substance and the different parameters p,, p ₂ and p ₃ of the attribute P. Here, the processor 130 is configured to derive the parameters p,, p ₂ and p ₃ by applying analytic calculations and/or numerical methods to equations describing relationships between the impedance measures z{\^), z(f ₂ ) and z(f ₃ ) of the liquid substance and the different parameters p ₂ and p ₃ .

According to another embodiment of the invention, the processor 130 is associated with a lookup table describing relationships between impedance measures z^), z(f ₂ ) and z(f ₃ ) of the liquid substance and the different parameters p,, p ₂ and p ₃ of the attribute P. The processor 130 is here configured to derive the parameters p,, p ₂ and p ₃ from said lookup table, either directly or by interpolating between values therein.

In any case, to further illustrate how to obtain the above-mentioned matrix K ^"1 needed to derive the attribute P and any parame- ters p, thereof, we will now discuss a general example.

Let us assume that n designated coil members, or sensors, are to be used, and that the below linear system of equations describe the relationships between n different parameters p,, p ₂ , p _n of a liquid substance to be determined, coefficients b ₂ , b _n represent respective weight coefficients of the n sensors, C c ₂ , c _n designate respective constants for the n sensors, and z(f _! ), z(f ₂ ),..., z(f _n ) are the respective readings of the n sensors. z(f _! ) = a ₁ b ₁₁ +a ₂ b ₁ 2+a ₃ b ₁ 3+ ... +a _n 1 _n +C 1

z(f ₂ ) = a _! b ₂₁ +a ₂ b ₂₂ +a ₃ b ₂₃ + ... +a _n b _2n +C ₂ z(f _n ) = a ₁ bni+a ₂ bn2+a3 n3+ ■■■ +a To determine the n weight coefficients b ₂ , b _n and the respective constant for each sensor, i.e. c c ₂ , c _n it is necessary to study n+1 different samples with known composition of the liquid substance. Thus, for the first sensor, the following system of linear equations shall be solved: z(f 1 ) = a b +a^b^+a^b^ +a _1n b _1n +c ₁ z(f ₂ ) = a ₂ i b- _\ 1 +a ₂ 2bi2+a ₂ 3bi3+ +a _2n bi _n +Ci z(f _n+1 ) = 3( _η +ΐ)ΐ b-i 1 + a( _n +i)2bi ₂ +a(n+i)3bi3+ ... +a _(n+1)n b _1n +c ₁ where:

a , a ₁₂ , a _1n , a _2n , a _(n+1)n are n known portions of components, or parameters, in n+1 respective samples of the liquid substance; b , b ₁₂ , b _1n are n weight coefficients of the n known portions of components, or parameters, for the first sensor to be determined;

C! is the constant of the first sensor to be determined; and z{ ), z(f ₂ ), z(f _n+1 ) are the respective readings of the first sensor for the n+1 known samples.

Analogously, for the second sensor, the following system of linear equations should be solved: ^z (f 1 ) = 3i 1 b ₂ i +a ₁₂ b ₂₂ +a ₁₃ b ₂₃ + +a _1n b _2n +c ₂ z(f ₂ ) = a ₂ i b ₂ i +a ₂₂ b ₂₂ +a ₂ 3b ₂ 3+

···

z(f _n+ i) = 3( _n+ i )i b ₂ i +3( _n+ i ₎₂ b ₂₂ +3( _n+ i )3b ₂ 3+ ... +a( _n+ i _)n b _2n +c ₂ where:

a , a ₁₂ , a _1n , a _2n , a _(n+1)n are n known portions of components, or parameters, in n+1 respective samples of the liquid substance; b ₂ i . b ₂₂ , b _2n are n weight coefficients of the n known portions of components, or parameters, for the second sensor to be determined; c ₂ is the constant of the second sensor to be determined; and z{ ), z(f ₂ ), z(f _n+1 ) are the respective readings of the second sensor for the n+1 known samples.

Similarly, for the n:th sensor, the following system of linear equations is to be solved: ^z (f 1 ) = ¾i 1 b _n i +a ₁₂ b _n2 +a ₁₃ b _n 3+ +a _1n b _nn +c _n z(f ₂ ) = a ₂₁ b _n1 +a ₂₂ b _n2 +a ₂₃ b _n 3+ +a _2n b _nn +c _n ^z (fn + l) = ^(n + 1)1 b _n i +a(n ₊ i ₎₂ b _n2 + a(n ₊ i)3b _n 3+ ... +Ά ( _{n +} i ) _n b _nn +C _n where:

a , a ₁₂ , a _1n , a _2n , a _(n+1)n are n known portions of compo- nents, or parameters, in n+1 respective samples of the liquid substance; b _n1 , b _n2 , b _nn are n weight coefficients of the n known portions of components, or parameters, for the n:th sensor to be determined; c _n is the constant of the n:th sensor to be determined; and z{ ), z(f ₂ ), z(f _n+1 ) are the respective readings of the n:th sensor for the n + 1 known samples.

Preferably, according to embodiments of the present invention, a computer program is used for calculating the attribute P, i.e. the parameters p,, such as various component concentrations in the liquid substance. The computer program may be based on the so-called Gaussian elimination method. Thereby, first the method is used for solving systems of linear equations in which unknown factors are the weight coefficients by, and the known fac- tors are the samples data a _nm and the sensor readings z(f _k ). Se- cond, the method is used to verify the calculated weight coefficients b . In this case, unknown factors are the parameters p, , and known factors are the weight coefficients b and the sensor readings z(f _k ) . According to one embodiment of the present invention , the liquid substance is represented by milk, e.g. from cows. By using a standard potentiostat, areas may be identified where an electrical field influences specific milk components. Figure 4 shows a respective graph over five different milk samples scanned in a frequency range from 1 0 M Hz to 1 70 M Hz. Each graph represents a measured voltage over a sensor coil as a function of frequency for a given milk sample. Hence, the voltage reflects any impedance variations in the milk in response to the variation in frequency. As can be seen, there is a substantial overlap among the graphs. However, there are also frequency areas where the graphs split in amplitude. These are areas where phenomena of impedance spectroscopy take place, and consequently the influence of each individual constituent in the milk can be defined. In this case, a frequency area f _a was found especially interesting. To identify the most suitable frequencies, two milk samples were used with the minimum and maximu m concentrations of a target constituent. Hence, to identify the best frequency for monitoring fat content, milk samples were used which had 0% fat (skim milk) and 3,2% fat (pasteurized milk) respectively. Analogously, two milk samples were used which both had 3,2% fat (pasteurized milk and farm milk respectively) , but which had great difference in somatic cell count (SCC) to define an appropriate frequency where SCC has the strongest response; and finally, two milk samples were used that contained 3,2% pasteurized lactose free milk and 3,2% pasteurized regular milk to find a frequency where lactose has the strongest response.

Using the potentiostat, in this case it was determined that at a first frequency f i , around 1 1 7 M Hz, sufficient sensitivity was obtained for measuring lactose; at a second frequency f ₂ , around 1 26 M Hz, sufficient sensitivity was obtained for measuring SCC; and at a third frequency f ₃ , around 1 31 M Hz, sufficient sensitivity was obtained for measuring milk fat. Based on this information, in turn , measuring circuitries were designed , which are optimized for registering parameters around the identified frequen- cies f 1 , f ₂ and f ₃ respectively.

Table 1 below shows an example of impedance measurements on milk samples in respect of fat content, lactose content, content of solids non-fat (SN F) , direct microscopic somatic cell count (DMSCC x 1 000) , milk urea nitrogen (M UN) content, true protein content (i .e. compensated for non-nitrogen protein) , content of other solids and the total content of solids. Here, the deviations are relatively small for all parameters, except for the direct microscopic somatic cell count. Therefore, for all parameters other than the direct microscopic somatic cell count, the de- pendence of milk impedance on them can be assumed to be linear. On the other hand, among the samples the difference in the direct microscopic somatic cell count runs up to tenfold nu mbers. Consequently, the dependence of milk impedance on this parameter cannot be assumed linear. Instead, to take proper account of the milk impedance behavior on the direct microscopic somatic cell count a (DMSCC x 1 000) ² colu mn has been added to Table 1 .

Table 1 Table 2 below represents the result of measurements on milk samples made by a set of measuring circuitries, which each includes at least one coil member of the proposed type. Each measuring circuit here measured a resonance amplitude and re- sonance frequency respectively. According to the invention, both these values can be used for weight coefficient determination.

Table 2

As mentioned earlier, weight coefficients are preferably calculated and verified for each coil member used. In this example, using all 1 2 samples, the following data were obtained .

Coil #2 Coil #3

Coef #0 -269.972 Coef #0 82.7843

Coef #1 1 494.87 Coef #1 -783.27

Coef #2 -870.284 Coef #2 365.625

Coef #3 0.1 27788 Coef #3 0.61 1 944

Coef #4 32.4042 Coef #4 -9.53335

Coef #5 -0.00032 Coef #5 -0.00037

Coef #6 1 0457.7 Coef #6 6782.92 Coil #1 (frequency) Coil #3 (frequency)

Coef #0 -0.18802 Coef #0 -0.02558

Coef #1 1.35308 Coef #1 -0.06935

Coef #2 -0.91157 Coef #2 0.1982

Coef #3 -0.00014 Coef #3 8.12E-06

Coef #4 0.030247 Coef #4 0.007431

Coef #5 -7.93E-08 Coef #5 -4.65E-08

Coef #6 149.208 Coef #6 35.405

Coil #2 (frequency) Coil #4 (frequency)

Coef #0 -0.0711 Coef #0 0.01687

Coef #1 0.44809 Coef #1 -0.15115

Coef #2 -0.22186 Coef #2 0.109708

Coef #3 0.000268 Coef #3 -6.03E-06

Coef #4 0.003958 Coef #4 -0.00196

Coef #5 -3.83E-07 Coef #5 1.26E-08

Coef #6 91.7391 Coef #6 8.97722

Table 3 below shows the values of Table 1 after correction according to the above weight coefficients #0, #1, #2, #3, #4, #5 and #6 for the coil members #1 , #2, #3 and #4.

Table 3 Two different strategies are proposed to check the consistency of the weight coefficients. The consistency can be proved by calculating weight coefficients of selected parameters of the liquid substance (e.g. the above milk components) using diffe- rent sets of samples. Then, if the coefficients are sufficiently close, it is safe to assu me that these coefficients are "correct", and may therefore be used for real/live measurements on a liquid substance.

An alternative way to prove the consistency of the components' weight coefficients is to use the coefficients for calculating a nu mber of parameters of the liquid substance for samples which were not used in order to determine the coefficients, and then compare the calculated parameters with the measured data.

If the liquid substance is milk, SCC is one of the most important parameters to characterize the quality of this liquid food product. Today, various methods based on photolu minescence are employed for counting the number of somatic cells in milk samples. However, said methods cannot easily distinguish between different types of leukocytes. Namely, in milk there are 4 to 6 major groups of somatic cells, which differ in terms function, size and electrical polarization. Typically, in uninfected milk with a SCC less than 300,000 we find the following concentrations of somatic cells : macrophages 60%, leucocytes 28% and polymorphonuclear leukocytes (PMN) 5-25%. Nevertheless, if the milk originates from an animal having mastitis, the SCC is normally above 400,000, and here the concentration of PM N may be as high as 90% (usually in a range from 33% to 49%) . Hence, by studying parameters related to leukocytes a more accurate measure of the animal health can be obtained than by monitoring SCC alone. Different leukocytes have different polarization frequencies, and may thus be individually identified via the present invention.

Figure 4 schematically illustrates a monitoring unit 1 00 according to an embodiment of the invention, wherein the dielectric conduit 110 has been spilt into two separate branches 111 and 112 respectively. This enables a relatively compact design of the monitoring unit 100. Naturally, any number of parallel branches is here conceivable if further compactness is desired. In any case, each branch 111 and 112 is configured to transport a fraction of the flow F, and each branch is encircled by at least one coil member 121, 122, 123 and 124 respective 125, 126, 127 and 128, which coil member is included in a measuring circuitry. Preferably, for sensitivity reasons the electromagnetic field produced by each coil member has a relatively narrow spectral range with a distinct center frequency that is located in an interval from 1 MHz to 1500 MHz. Thus, based on the coil members 121, 122, 123, 124, 125, 126, 127 and 128 up to eight parameters of the attribute P can be determined by the proces- sor 130.

It is also preferable if the processor 130 includes, or is associated with, a computer readable medium M, e.g. in the form of a memory module, such that the processor 130 has access to the contents of this medium M. Furthermore, a program is recorded in the computer readable medium M, and the program is adapted to make the processor 130 control the process described above, as well as the embodiments thereof further elaborated on below, when the program is run on the processor.

In order to sum up, we will now describe the general method ac- cording to the invention with reference to the flow diagram in Figure 5.

In a first step 510, a flow of liquid substance is transported through at least one dielectric circuit. In parallel there with, a step 520 subjects the liquid substance to at least one electro- magnetic field having a respective spectral range. A step 530 parallel with steps 510 and 520, registers at least one impedance measure, or more precisely registers a number of impedance measures equal to the number of electromagnetic fields applied in step 520. Then, based on the registered impedance measures, a step 540 determines at least one parameter characterizing an attribute of the liquid substance. More precisely, step 540 determines up to as many parameters as the number of impedance measures registered in step 530.

Thereafter, the procedure either ends, or loops back to steps 51 0 to 530 for continued measuring and updating of the at least one parameter.

All of the process steps, as well as any sub-sequence of steps, described with reference to Figure 5 above may be controlled by means of a programmed computer apparatus. Moreover, although the embodiments of the invention described above with reference to the drawings comprise computer apparatus and processes performed in computer apparatus, the invention thus also extends to computer programs, particularly computer programs on or in a carrier, adapted for putting the invention into practice. The program may be in the form of source code, object code, a code intermediate source and object code such as in partially compiled form, or in any other form suitable for use in the implementation of the process according to the invention. The program may either be a part of an operating system, or be a separate application. The carrier may be any entity or device capable of carrying the program. For example, the carrier may comprise a storage medium, such as a Flash memory, a ROM (Read Only Memory) , for example a DVD (Digital Video/Versatile Disk) , a CD (Compact Disc) or a semiconductor ROM , an EP- ROM (Erasable Programmable Read-Only Memory) , an EEPROM (Electrically Erasable Programmable Read-Only Memory) , or a magnetic recording medium, for example a floppy disc or hard disc. Further, the carrier may be a transmissible carrier such as an electrical or optical signal which may be conveyed via electrical or optical cable or by radio or by other means. When the program is embodied in a signal which may be conveyed directly by a cable or other device or means, the carrier may be consti- tuted by such cable or device or means. Alternatively, the carrier may be an integrated circuit in which the program is embedded, the integrated circuit being adapted for performing, or for use in the performance of, the relevant processes.

Although the invention has been described in relation to the analysis of milk, it is equally well adapted for analyzing other liquids, especially other liquid food products. The parameters to be analyzed will in those cases of course be different than the ones described in this specification.

Although the invention is advantageous in connection with cow milking, the invention is equally well adapted for implementation in milking machines for any other kind of mammals, such as goats, sheep or buffaloes.

The term "comprises/comprising" when used in this specification is taken to specify the presence of stated features, integers, steps or components. However, the term does not preclude the presence or addition of one or more additional features, integers, steps or components or groups thereof.

The invention is not restricted to the described embodiments in the figures, but may be varied freely within the scope of the claims.