The present invention relates to a method for measuring the moisture ratio of organic material, preferably wood, and apparatus herefor.
It is often necessary or desirable in various contexts to determine the moisture ratio of a material, thereby to obtain a measurement of its moisture content.
One such context in which this applies involves the treatment of wood pieces in the so-called sorting department of a saw mill, in which knowledge of the water content of the wood pieces is necessary in order to optimize the preceding drying stage thereof. The invention can also be applied when wishing to measure the moisture ratio of seed and wood chips.
The invention is not limited to any particular field of application, however, but relates to the measurement of the moisture ratio of organic material.
One method of measuring the moisture ratio of organic material is to produce an electric or electromagnetic field therein. The propagation or spread factor of the field is influenced chiefly by the water contained in the material. The water content is normally determined by measuring resistance or capacitance. These values are a measurement of a complex material-constant in the water-bearing material. The material constant is known as the complex, dielectric constant and comprises a true component and an imaginary component ε' and ε", where ε' is the dielectric constant and ε" is a measurement of the loss in the material, and where tan δ=ε"/ε'.
The moisture ratio is defined as ##EQU1## Wm is the weight of the material when wet, and Wd is the weight of the material when dry.
When the density of the material in its dry state is known and is constant, the moisture ratio can be determined simply by measuring the water content. However, the density of many organic materials, such as wood for example, varies. In order to obtain a true value of the moisture content, it is necessary to measure both the water content and the density of the material.
It is also known to determine the water content and density of organic materials with the aid of microwaves. Such a method is found described in the reference Kraszewski, A. "Microwave Aquametry - A review", Journal of Microwave Power, 15(4), 1980 p.p. 209-220. It can be said in general that only one transducer is needed when practising such a method. In order to be able to determine two values, namely water content and density, it is necessary to employ two parameters. In this respect there is measured the damping and the phase shift or change of a microwave signal transmitted through the material. These two values are related to the permitivity in accordance with the following relationship: ##EQU2## where f is frequency, ω is angular velocity, d is thickness, C is the speed of light, A is damping (dB), and F1 is phase shift (degrees).
Damping or attenuation A is a measurement of the water content and the phase shift F1 is a measurement of the density and the water content.
In the aforesaid Kraszewski reference it is given in formula 12 that A=F1 (Ww, Wd) and that F1=F2 (Ww, Wd) where Ww is the weight of the water in the material, Wd, and A and F1 have the aforesaid significance, and F1 and F2 represent different functions.
Kraszewski maintains that the moisture ratio can be determined independently of both the density and the thickness of the material. This does not apply in general to the aforesaid formula 12, which describes the model in a general form, but solely when using the special model incorporated in Kraszewski's formula 13, namely ##EQU3## V is the volume of the material and where k1 -k4 are constants.
However, it is not possible to measure the moisture ratio of an organic material with the aid of formula 13 without taking density and thickness into account, other than in the case of a material with which the damping (A) as a function of the water content per unit of area (Ww /Y) can be described with a single straight line irrespective of thickness, and with which the phase shift or change (F1) as a function of the water content per unit of area (Ww /Y) does not vary with thickness. The relationship A as a function of Ww /Y and F1 as a function of Ww /Y respectively is discussed hereinafter.
The problem of measuring with a high degree of accuracy the moisture ratio of organic material whose thickness and/or density varies is solved by the present invention.
Accordingly, the present invention relates to a method for measuring the moisture ratio of material, particularly of organic material, in which microwaves are transmitted to the material whose moisture ratio is to be determined; in which the transmitted microwaves are received subsequent to their transmission through said material; in which the damping (A) of the microwaves and their phase change (F1) upon transmission through the material is determined; and in which the so-called moisture ratio of the material is established with the aid of the determined damping (A) and phase-change values, the method being characterized by determining the term (Ww /Y) from the expression ##EQU4## where A is the measured damping and t is the thickness of the material, and where k1 and k3 are constants and are determined by a calibration process, and where the term Ww /Y is the weight of water per unit area; and by determining the term (Wd /Y) from the expression ##EQU5## where F1 is the measured phase change, in which expression the aforedetermined term Ww /Y is used, where t is the thickness of the material, the term Wd/Y is the weight of the material in a dry state per unit area, and where k4, k5, k6 are constants and determined by a calibration process; and by deriving the moisture content in a known manner through the ratio (Ww /Y)/Wd /Y).
The invention also relates to apparatus of the kind set forth in claim 4 and having the characteristic features disclosed therein.
The present invention is based on the understanding that the damping A is influenced by the thickness of the material, rather than the variable (Wd /Y), i.e. the weight of the material in a dry state per unit area, and that in the case of the majority of organic materials the phase change or shift (F1) varies as a function of the water content per unit area (Ww /Y) with the thickness.
By unit area Y above is meant the area perpendicular to the direction of propagation of the microwaves. The terms (Ww /Y) and (Ws /Y) are abbreviated from the respective expressions ##EQU6## where V is the volume of the material and t the thickness thereof.
In several areas where the invention can be applied, no advantage is gained by attempting to eliminate the dependency of the thickness of the material and therewith obtain a less accurate value of the moisture ratio. On the contrary, the thickness of the material is often known. Such is the case in saw mills for example, where a thickness measuring device is often on hand or can be readily installed. By measuring the thickness of the material, or when having knowledge of its thickness, the moisture ratio can be measured with a high degree of accuracy with the aid of the present invention, even in those cases where the density of the material varies. The invention therefore affords a particular advantage when desiring to measure the moisture ratio of wood. The density varies quite considerably between different wood species, such as between pine and spruce.
The present invention will now be described partially with reference to the separate figures of the accompanying drawings, of which
FIG. 1 illustrates damping per unit of thickness as a function of the weight of water in the material per unit of volume;
FIG. 2 illustrates damping as a function of the weight of water per unit area in respect of a material of thicknesses t1, t2 and t3 respectively;
FIG. 3 illustrates the phase change per unit of thickness as a function of the weight of water in the material per unit of volume;
FIG. 4 illustrates the phase change as a function of the weight of water per unit area in respect of a material having thicknesses t1, t2 and t3 and densities d1 and d2 ; and
FIG. 5 is a block schematic of apparatus or an arrangement for carrying out the method according to the invention.
According to the invention the damping A is determined by the relationship ##EQU7## where t is the thickness and k1, k2 and k3 are constants.
In the case of a microwave frequency of say 16 GHz, the constant k2 is so small that the term Wd /Y·k2 can be ignored.
There is then obtained ##EQU8##
According to the invention the phase change F1 is determined from the relationship ##EQU9## where k4, k5, k6 are constants.
With regard to the relationship (1) above, FIG. 1 shows the principle relationship between damping per unit of thickness A/t and the weight of the water in the material per unit of volume Ww /V.
Subsequent to passing a knee, the relationship is approximately linear. The function is approximately independent of the term Wd /Y. This non-dependency increases with an increase in microwave frequency. The relationship shown in FIG. 1 illustrates the conditions prevailing at a microwave frequency of 16 GHz. At higher frequencies, up to 20-24 GHz, the dependency on the term Wd /Y is practically zero. The term Wd /Y cannot be ignored at lower frequencies, for example 2450 MHz, where the dependency is greater.
From a general standpoint, the damping or attenuation is affected more precisely by the water contained in the material and the phase change is affected to a still larger extent by the density at higher frequencies. This influence has a maximum at 20-24 GHz.
Consequently, in accordance with a preferred embodiment there is used a microwave frequency above 12 HGz, and beneath 24 GHz, preferably 16 GHz.
FIG. 2 illustrates the relationship between the damping A as a function of the water weight per unit area (Ww /Y) in respect of a material having the thicknesses t1, t2, t3, i.e. the same relationship as that illustrated in FIG. 1 with the variables having been multiplied by t. FIG. 2 shows three curves corresponding to the curve shown in FIG. 1. All curves have a common inclination beneath the respective knees of each curve.
The physical explanation as to why a knee is obtainned at low water concentrations is because the water is bound more firmly to the organic material and does not therewith influence the microwaves to the same extent as the more loosely bound water molecules at higher water concentrations.
The relationship (2) according to the invention derives from the curve to the right of the knee. This section of the curve describes the damping A as a function of the water weight per unit area Ww /Y. The knee occurs at very low water concentrations. The knee referenced (Ww /V)k in FIG. 1 corresponds to the knee (Ww /Y)k in FIG. 2. According to the relationships in FIGS. 1 and 2 ##EQU10## Thus, when studying the function in respect of a material piece having the thickness t2 =2t1 the following relationship applies: ##EQU11##
It will be seen from FIG. 2 that the last term in the relationship (2) describes where respective curves cut the A-axis, i.e. at a value of k3 ·t1,2k3 ·t1 and 3·k3 ·t1 respectively.
Thus, the point of intersection with the A-axis depends on the thickness t and not on the term Wd /Y. In the aforecited relationship according to Kraszewski ##EQU12## it is said instead that damping A depends on ##EQU13##
In practice, however, there is found a relationship between the dependency of the damping on Wd /Y and on t. The major part of this relationship is proportional. A much better result is obtained when the thickness is used as a variable instead, due to the fact that physically it is the thickness t which influences the damping instead of Wd /Y.
Thus, the relationship (2) above is utilized in accordance with the invention in respect of the influence of the damping A. The thickness t is measured and inserted into the relationship or equation, as is also the value of the measured damping A. Thus, the sole unknown term (Ww /Y) can be obtained from the equation (2).
This is inserted in the equation (3) above together with the value of the measured thickness t.
As beforementioned, the invention is also based on the understanding that in the case of the majority of organic materials the phase change or shift (F1) is influenced as a function of the water content per unit area (Ww /Y) with the thickness, in a corresponding manner to that in which the damping (A) is affected. This influence has the same physical basis as that given above.
FIG. 3 illustrates the phase change per unit of thickness (F1/t) as a function of the water content per unit volume (Ww /V). As will be seen, each and every one of the curves has approximately the same appearance as the curve in FIG. 1, exhibiting a characteristic knee. It will also be seen that the density of the material affects the function, such that the curve representing the higher density is moved upwards parallel with the lower curve. The density d2 is higher than the density d1.
FIG. 4 illustrates the phase change (F1) as a function of the water content per surface unit (Ww /Y).
As will be seen from FIG. 4, curves corresponding to those illustrated in FIG. 2 are obtained, where a plurality of parallel curves extend from a first inclined curve section. As will be seen from FIG. 3, the position of the first inclined curve section along the F1-axis depends upon the density. Extending from this first curve section are mutually parallel curves whose positions are determined by the thickness t of the material. As will be seen from FIG. 4, the curves are displaced towards a higher value on Ww /Y in respect of a greater thickness, where t3 >t2 >t1.
FIG. 4 also indicates that the phase change F1 depends upon the thickness of material of the kind with which the relationship of F1/t to Ww /V exhibits a knee, i.e. at least the majority of organic materials.
The phase change F1 is measured and also inserted into the relationship or equation (3), whereafter the sole unknown term Wd /Y is obtained therefrom.
Because the term (Ww /Y) is determined by taking into account the thickness of the material and is thus far more accurate than would be the case if no attention had been paid to the thickness dependency, the value of the term Wd /Y taken from the relationship (3) affords a higher degree of accuracy.
By also realizing that the phase change or shift F1 is directly affected by the thickness, and therewith inserting the term k6 ·t, c.f. the relationship (3), a very accurate result is obtained with respect to the term Wd /Y.
Since the density of certain organic materials, such as dry wood, varies considerably between, for example, spruce and pine, it is essential that the term Wd /Y is determined very accurately, since it cannot be considered to be constant when applying the invention in a saw mill for example.
The constants k1, k3, k4, k5, k6 are determined with the aid of a conventional calibration process. ##EQU14## is obtained from the resultant values of the terms ##EQU15##
Because, as mentioned, the relationship (2) determines the water content as a function of the thickness of the material and the term (Ww /Y) is inserted in the relationship (3), by means of which term (Wd /Y) is determined as a function of the thickness, a correct value of the moisture ratio is obtained even when the thickness and/or density of the material varies.
Because a curve characteristic according to FIGS. 2 and 4 is formed from the common incline for organic material beneath the knee, the present invention can be applied with particular advantage for determining the moisture ratio of organic material where the thickness and/or density thereof varies from case to case or changes within a single piece of material.
A complicating factor when measuring the moisture ratio of organic material is that many organic materials swell, i.e. increase in volume, with increasing water content. Wood is a typical example of such material. In this case, it is preferred to compensate the dependency of the density on the water content.
According to one extremely important embodiment of the invention compensation is made for the dependency of the material density on the water content of said material. In this case the density of the material is measured at the prevailing water content with the aid of the relationship ##EQU16## where the term (Wdd /Y) is the term (Wd /Y) compensated for density variations due to the prevailing water content of the material.
The expression is a second degree polynomial. A relationship Wdd =f(Ww) is measured empirically for respective materials, wherewith the constants k7 and k8 are determined.
The values obtained from the relationships (2) and (3) in respect of the terms ##EQU17## are inserted into the relationship or equation (4).
The moisture ratio is then obtained in a manner analogous with the aforegoing, by utilizing the term Wdd /Y, i.e. ##EQU18##
The moisture ratio thus obtained provides an extremely accurate value of the water content of the material.
It will be evident from the aforegoing that the present invention solves the aforementioned difficulties of measuring the moisture ratio of organic material with a high degree of accuracy, when the thickness and/or the density of the material varies and, in addition, when the density varies with the water content of the material.
FIG. 5 illustrates schematically an arrangement for measuring the moisture ratio of materials, incorporating apparatus according to the present invention.
The arrangement 1 located within the region defined by the chain line is of a known kind and will therefore only be described in brief.
The arrangement 1 includes a microwave oscillator 2 which generates a signal of frequency f0. This signal is passed to a so-called 3 dB-hybrid 3, where the signal is divided between a reference channel 4 and a measuring channel 5. A second oscillator 6 generates a signal of frequency f1, which is substantially lower than the frequency f0.
According to one preferred embodiment of the invention, the frequency f0 exceeds 12 GHz and lies beneath 24 GHz, preferably 16 GHz. The frequency f1 is therefore suitably about 10 KHz. The frequency f0 is modulated in a modulator 7 with the frequency f1, so as to form the frequency f0 +f1. The signal formed is transmitted by means of a first microwave antenna 8 and is received on a second microwave antenna 9.
The material whose moisture ratio is to be determined, for example a woodpiece 10, is introduced between the pair of antennae 8,9. The signal received in the measuring channel 5 is mixed with the signal in the reference channel 4 in a mixer 11, such as to form a signal having the frequency f1, and is delivered through a conductor 12 to a signal processing circuit 13, to which a signal is also applied directly from the said second oscillator 6, through a conductor 14. The phase change or shift F1 and the damping A are determined in a known manner in the signal processing circuit 13, by comparing the signals on lines 12 and 14 respectively. A signal corresponding to the phase change F1 and a signal corresponding to the damping A is formed on the outputs 15,16 of the signal processing circuit.
In accordance with the present invention, the signals F1 and A respectively are passed through respective conductors 15 and 16 to a calculating unit which incorporates a microprocessor 17 or the like. The microprocessor 17 incorporates a memory 18, which suitably comprises a RAM-memory or the like. There is also provided a suitable means 19 of known kind for measuring the thickness of the material piece 10 whose moisture ratio is to be determined and/or an infeed unit 20, such as keyboard or the like, for feeding in the thickness of the material piece 10, this thickness being determined in some other way or being known. The keyboard 20 may be provided with a display 21 on which the thickness fed into the microprocessor is displayed.
The means 19 and the infeed unit 20 are connected to the memory 18 in which the pertinent thickness is stored.
According to the present invention the computer unit is adapted to calculate the value of the term (Ww /Y) in the expression ##EQU19## where A is the damping obtained on the conductor 16 and t is the thickness of the material obtained from the memory 18, and where k1 and k3 are constants, which are also found stored in the memory 18. The constants k1 and k3 are determined by a suitable conventional calibration process, whereafter the values obtained are introduced into the memory. The computer unit introduces the value of the term (Ww /Y) into the memory 18.
The computer unit 17 is also adapted to calculate the value of the term (Wd /Y) in the expression ##EQU20## where F1 is the phase change or shift obtained from the conductor 15 and t is the thickness of the material. In this case, the computer unit is adapted to collect the value of the term (Ww /Y) from the memory and insert said value into the expression for F1. The constants k4, k5, k6 are determined by a suitable conventional calibration process, whereafter the values obtained are introduced into the memory.
The constants, k1, k3, k4, k5 and k6 are determined prior to commencing the measuring process in question.
The computer unit 17 is adapted to then form the moisture ratio, by carrying out the calculation (Ww /Y)/(Wd /Y).
In accordance with one preferred embodiment of the invention, which can be applied to advantage when the material whose moisture ratio is to be determined has properties such that its density is a function of its water content, the computer unit is adapted, subsequent to determining the value of the terms (Ww /Y) and (Wd /Y) respectively to calculate the value of the term (Wdd /Y) from the expression ##EQU21## while utilizing the value of the terms (Ww /Y) and (Wd /Y), where the term (Wdd /Y) is the term (Wd /Y) compensated for density variations due to the prevailing water content of the material. The constants k7 and k8 are also determined by means of a suitable conventional calibration process and introduced into the memory 18. This calibration process thus includes an empirical determination of the aforesaid expression in respect of a pertinent material.
The computer unit 17 is adapted to then form the moisture ratio, by making the calculation
(W.sub.w /Y)/(W.sub.dd /Y).
The value of the moisture quotient (Ww /Y)/(Wd /Y) or, when applicable, (Ww /Y)/(Wdd /Y) is fed to a suitable visualizing device 22, such as a display unit, a writer, lamps or the like.
It will be understood that the block schematic illustrated in FIG. 5 and associated descriptive text, particularly with reference to the arrangement located within the region defined by chain lines, shall be considered an embodiment where the invention is applied. Any other array of apparatus or devices for calculating phase change and damping can be used together with the present invention.
Instead of the arrangement illustrated in the region defined by said chain lines, there can be used to advantage an arrangement described in the Swedish patent application No. 8403438-8.
Thus, the invention shall not be considered to be restricted to the aforedescribed embodiments and it will be understood that modifications can be made within the scope of the following claims.