This invention relates to a method and equipment for real-time or on-line measurement of intact wood or pulp fibers, more particularly for measurement of a fiber for the microfibril angle and the path difference, a parameter proportional to the cell wall thickness. The method and equipment can be modified for real-time or on-line measurement of other birefringent samples, including retardation films and waveplates.
Wood or pulp fibers are closely related to paper properties. The increasing demand on high quality paper products requires optimal and more efficient use of the available wood fiber resources. Fibers' properties vary widely and they are different from fiber to fiber even within a tree. Research and analytical tools for measuring the properties of single pulp fibers are essential for a better utilization of available wood resources. The basic characteristics of a fiber include the fiber's length, width, shape, microfibril angle (MFA) and cell wall thickness (CWT). On-line measurement equipment is commercially available only for the fiber length, width and shape. The MFA and CWT are difficult to measure due to the fiber's two-wall structure.
A wood fiber is made of a primary wall enveloped in lignin to form the middle lamella and a secondary wall comprising three secondary layers, called S1, S2 and S3 layers (e.g. ref. Preston, R. D. (1974), The physical biology of plant cell walls, Chapman and Hall Ltd.). All the three secondary layers are concentric and composed of cellulosic microfibrils, embedded in an amorphous matrix of hemicelluloses and lignin. The outer secondary layer S1 and the inner secondary layer S3 are very thin, and their microfibrils are wound almost transversely to the fiber axis. The middle layer S2 contains the majority of the cell-wall material (80-95%) (see Page, D. H. (1969), Journal of Microscopy 90, 137-143 and Prud'homme, R. E. and Noah, J. (1975), Wood and fiber 6, 282-289) and it is widely accepted that a wood or pulp fiber can be approximated with the S2 layer. The microfibrils of the S2 layer trace a steep spiral around the fiber axis with the microfibrils of the layer's front and back walls crossed. The angle between the fibrillar direction and the fiber axis is termed the microfibril angle or S2 fibril angle of the fiber.
The experimental evidences provided by a number of investigations (e.g. see Mark, R. E. and Gillis, P. P (1973), Tappi 56, 164-167; El-Hosseiny, F. and Page, D. H. (1975), Fibre Science and Technology 8, 21-30; Page, D. H., El-Hosseiny, F., Winkler, K., and Lancaster, A. P. S. (1977), Tappi 60, 114-117; and Page, D. H. and El-Hosseiny, F. (1983), Journal Pulp and Paper Science, TR 99-100) indicate that the MFA φ is closely related to the mechanical properties of the fibers, such as the strength, the elastic modulus and the shrinkage. The CWT is a parameter often involved in pulp fiber measurement. For instance, the CWT is related to the fiber flexibility, strength and collapsibility (e.g. Dinwwodie, J. M. (1965), Tappi 48, 440-447; Horn, R. A. (1974), USDA For. Serv. Res. Pap. FPL 242, Horn, R. A. (1978), USDA For. Serv. Res. Pap. FPL 312 and Jang, H. F. and Seth, R. S. (1998), Conference proceeding. 84th annual meeting, Canadian Pulp and Paper Association, Montrel 27-30 Jan. 205-212) and it is also related to the surface quality and optical properties of paper.
To determine the MFA and/or the CWT, methods have been developed and used, for example the striation observation, angle of the slit pits, iodine staining, X-ray diffraction, confocal laser scanning microscopy (CLSM) and imaging ellipsometry. The first three techniques are tedious and only applicable to some wood species. The X-ray diffraction method generally is suitable for giving a measure of the mean microfibril angle of a piece of wood consisting of a few hundred fibers. Additionally, the application of the X-ray technique relies heavily on fiber geometry that is uncertain (see e.g. Prud'homme R. E. and Noah, J. (1975), Wood and fiber 6, 282-289).
The CLSM can be used for optically sectioning a wood or pulp fiber and generating its cross-sectional image. With the help of image analysis, the CLSM can determine fiber's transverse dimensions, including the CWT (e.g. Jang, H. F., Robertson, A. G. and Seth, R. S. (1992), Journal of Materials Science, 27, 6391-6400). To generate a cross-sectional image, a fiber is optically scanned in the cross-sectional direction while the fiber is stepped in the perpendicular direction and cross-sectional image is reconstructed from a series such line scans. A fiber for measurement by this technique needs to be pretreated, for example dyed with fluorochrome dye, and oriented to be perpendicular to the scanning direction. By combining the optical sectioning ability of the confocal microscope with the difluorescence of fluorochromadyed cellulose, also the MFA can be measured (Jang, H. F. (1998), Journal of Pulp and Paper Science, 24, 224-230). This method needs not only pretreatment of fibers (dyed with special fluorochromes) but also rotation of the incident polarized light.
Imaging ellipsometry enables nondestructive determination of both the MFA and the path difference PD or phase retardation Δ, a parameter proportional to the CWT, without sample pretreatment (Ye, C. and Sundström, M. O. PCT Patent Application (1996), WO9610168). The path difference PD or phase retardation Δ is proportional to the cell wall thickness d as described by PD=Δλ/2π=d(n2−n1) or Δ=2πd(n2−n1)/λ, where n2−n1 is the birefringence of the wall material and λ the light wavelength. In most applications, it is important to know the distribution of fibers' wall thicknesses instead of their absolute values. In these cases, the absolute quantity is not necessary and a parameter like the path difference PD, which is proportional to the CWT, can directly be used for replacing the CWT, as experimentally demonstrated (Ye C, Räty, J., Nyblom, I., Hyvärinen, H. and Moss, P. (2001), Nordic Pulp & Paper Research Journal. 2(16), 143-148). Furthermore, the CWT can be determined from the PD with a proper calibration procedure.
In an early work of imaging ellipsometry, a theoretical model describing the two-wall structure of a pulp fiber was established and a multiple-wavelength method was developed based thereon for measuring both the retardation Δ and the MFA of single pulp fibers (Ye, C, Sundström, M. O., and Remes, K. (1994), Appl. Opt. 33, 6626-6637 and Ye, C. and Sundström, M. O. PCT Patent Application (1996), WO9610168). According to the new model, the two opposite cell walls (S2 layer) of a fiber is optically equivalent to two identical linear retarders arranged in series with their axes symmetrical around the fiber axis. The two retarders have the same phase retardation Δ, which is proportional to the thickness of the fiber walls and the birefringence of the wall material, and their orientation angles have the same value as the microfibril angle φ of the fiber, but with opposite signs.
The multiple-wavelength method enabled non-destructive measurement of wood or pulp fibers for the MFA φ and phase retardation Δ without sample pretreatment. The method employs a plane polariscope and it measures light intensities of a fiber sample at different wavelengths (at least two) by rotating the analyzer of the plane polariscope. With the intensity data obtained, intermediate results for φ and at each wavelength are calculated, from which the measurement results of φ and at each wavelength are determined according to the criteria that the MFA φ is a constant parameter, whereas the phase retardation Δ is a function of the light wavelength λ as described by Δ=2πd(n2−n1)/λ, where d is the thickness of fiber's cell walls and n2−n1 is the birefringence of the wall material. Then an optimal estimation for φ and Δ is obtained by comparing the measured results as a function of the light wavelength with the theoretical description according to the least squares principle. However, the fiber sample has to be aligned to a certain orientation. Due to this limitation, the measurement speed is restricted to be further increased.
Later a more advanced method (Ye, C. (1999) Appl. Opt. 38, 1975-1985) was developed based on the Mueller-matrix ellipsometry. The Mueller-matrix method permits non-destructive determination of the MFA and Δ of a wood or pulp fiber oriented arbitrarily by measuring the Mueller matrix of the fiber at one wavelength. This method has all the advantages of the multiple-wavelength method but not subject to the same limitation. Based on this method, a more powerful research tool (Ye C, Räty, J., Nyblom, I., Hyvärinen, H. and Moss, P. (2001), Nordic Pulp & Paper Research Journal. 2(16), 143-148) for characterization of pulp fivers was constructed, which allows a semi-automatic measurement of single wood or pulp fibers for the MFA and Δ with the measurement speed significantly enhanced. However, the Mueller-matrix method still needs the fiber sample keeping stationary during measurement and requires sequentially acquiring images from the fiber sample.
Most recently, another ellipsometric method (Jang, H. F. (2005), US Patent Application, US 2005/9122514 A1) is reported for measurement of the microfibril angle MFA and the phase retardation Δ. The method of Jang employs a circular polariscope (Theocaria, P. S. and Gdoutos, E. E., Matrix Theory of Photoelasticity, Springer-Verlag, New York, 1979, pages 117-123), which creates an image determined only by the sample's properties related to polarized light. With a circular polariscope it is not necessary to align the sample because the equipment is insensitive to the sample's orientation. This method is also a multiple-wavelength method and it measures light intensities of a fiber sample created with the circular polariscope at several well-separated wavelengths (at least three). The microfibril angle MFA and phase retardation Δ are determined from the measured intensity data by fitting the data with the theoretical description.
Theoretically it is feasible to automate the method of Jang for real-time or on-line measurement of wood or pulp fiber properties with the circular polariscope adapted to be able to simultaneously create and detect multiple images of the sample at individual wavelengths. As known, the light intensity emergent from a fiber, no matter in a plane or circular polariscope, is non-linearly related to the fiber's properties, the incident light and the light wavelength. The non-linear relationship between the intensity data and the unknowns to be determined implies that light intensity measurement at two wavelengths is not enough to determine two unknown parameters of the sample, e.g. the microfibril angle and phase retardation, in all cases, in which fibers have different cell wall thicknesses ranging to cover all possible values for the phase retardation Δ. In fact, a simulation calculation shows that ambiguous results can occur when determining φ and Δ even in case of measuring light intensities at three wavelengths. This is a liability for the method of Jang. To guarantee the measurement results reliable, the number of individual wavelengths, at which light intensity needs to measure, should be higher or much higher than stated by Jang. For a non real-time measurement, this is not a big problem. However, it will be completely different if multiple images need to be simultaneously created and detected to realize real-time measurement especially when the number of the wavelengths increases. In addition, the light intensity incident on the fiber sample is an additional unknown parameter, which needs to be determined to measure the microfibril angle and phase retardation. A practical system meeting the requirements above is not only complicated and expensive but also technically not desired because for example the light intensity of the multiple images will be further reduced with increasing number of wavelengths.
For pulp fiber characterization, valuable fiber quality information needs to be gathered from a representative sample of thousands of individual fibers so that a real-time measurement method is the most desirable, capable of measuring a great mass of fibers in sufficiently short time to provide more reliable results and statistical analysis. Furthermore, because of largely automated pulping process it is more important that a real-time measurement method can be used or adapted for use under the on-line condition to measure moving fibers to provide on-line feedback information of the fiber quality for pulp evaluation and controlling the production process. As described above, however, the methods so far available for determination of the MFA φ and the CWT or PD are either limited for use in laboratories or restricted for the liabilities.
It is therefore an object of the present invention to provide a method and equipment capable of real-time measurement of wood or pulp fibers for both the microfibril angle and the path difference without sample pretreatment.
It is another object of the invention to provide a better solution for the real-time measurement of wood or pulp fibers, which enables measurement of a fiber or a moving fiber oriented arbitrarily for the microfibril angle and phase retardation by acquiring only one image from the fiber and capable of generating the reliable measurement results by over-determining the parameters to avoid possible ambiguous data in signal processing.
The present invention provides a method developed based on the spectroscopic ellipsometry as a solution for this challenging task. The method of the invention employs a circular polariscope in combination with a line spectral camera with the former generating a polarizing micrograph of a pulp fiber under test, which is insensitive to the fiber's orientation, and the latter completing a real-time spectral analysis of the micrograph generated.
A circular polariscope (Theocaria, P. S. and Gdoutos, E. E., Matrix Theory of Photoelasticity, Springer-Verlag, New York, 1979, pages 117-123) has an optical arrangement comprising two achromatic quarter-wave retarders inserted between and oriented perpendicular or parallel to each other and at 45° to a pair of parallel or perpendicular polarizers with the fibers for measurement, preferably immersed and distributed in suspension, embraced by the retarders. The image of a wood or pulp fiber created by the arrangement is independent of the fiber's orientation and it is formed or determined only by the fiber's MFA φ and phase retardation Δ in addition to the incident light intensity. With this property, the equipment of the invention enables measurement of a fiber oriented arbitrarily.
Another aspect of the method of the invention is to simultaneously measure all the light intensities emergent from a fiber sample in a continuous spectral range, from which a real-time spectral analysis of the fiber image can be carried out and the unknown parameters φ and Δ can be over-determined to avoid any possible ambiguous results and to improve the measurement accuracy. For this purpose, the method of the invention uses a line spectral camera placed behind the circular polariscope. The light emergent from the exit polarizer of the polariscope is scanned by a line spectral camera, which preferably is an ImSpector (http://www.specim.fi/) followed by a CCD camera. The ImSpector captures a line image of the fiber's image and disperses light from the line image into a continuous spectrum, which is detected by the CCD camera, so that a real-time spectroscopic analysis is feasible. As the MFA φ is a constant parameter, while the retardation α is a function of the light wavelength λ, as described by Δ=2πd(n2−n1)/λ, where d is the thickness of fiber's cell walls and n2−n1 is the birefringence of the wall material, a spectroscopic analysis of the invention based on the least-squares principle results in an optimal estimation for φ and Δ.
A further aspect of the invention is to determine the incident light intensity without using additional component. This is desired or needed to carry out a spectral analysis of the sample image as mentioned above so that the unknown parameters φ and Δ can be determined by acquiring and using only one image created from the circular polariscope. Because the incident light intensity of a fiber sample in a practical system is also a function of the light wavelength λ, it is necessary to compensate the effect caused by the spectral transmission of the incident light intensity to measure the fiber's φ and Δ. In accordance with the invention, when a fiber is measured an image part is captured by the spectral camera, which contains the fiber's segment selected for measurement and a neighboring background image part without fiber. The CCD camera behind the ImSpector detects and outputs the spectrum I[Δ(λ),φ] of the selected fiber segment and the spectrum I0(λ) of the background image part as well. The background image part is near the fiber segment to measure and it can be used to approximately describe the light intensity transmitted by the equipment and detected by the CCD camera at the position of the fiber segment in case the fiber is absent. The spectrum I[Δ(λ),φ] is normalized with I0(λ) and the normalized spectrum I[Δ(λ),φ]/I0(λ), which is independent of the spectral transmission of the equipment or the incident light intensity, is used for the real-time spectroscopic analysis, from which an optimal estimation for the fiber's φ and Δ can be generated based on the least square principle.
As only one image from the fiber sample is needed, the method of the present invention exerts no restriction on the measurement speed and does not need special equipment having complicated structure for simultaneously creating and detecting multiple images from the sample. In addition, as experimentally demonstrated this exclusive feature of the invention allows measurement of a moving fiber oriented arbitrarily for φ and Δ, an assignment required for fiber measurement under the on-line condition.
The equipment of the invention is developed based on the invention and it works as a spectroscopic imaging ellipsometer without moving part capable of measuring wood or pulp fibers oriented arbitrarily by acquiring only one image from the fibers. Besides the wood pulp fiber, the method is suitable for measurement of other cellulose fibers. In addition, the method and equipment of the invention can be used or adapted for use for real-time or on-line measurement of ordinary birefringent samples including retardation films and waveplates. The principle, advantages and features of this invention will become more apparent from the following description in conjunction with the accompanying drawings.
A wood or pulp fiber under investigation in the present invention is assumed to be a flattened cylinder in shape and its two walls are of the same thickness and their microfibrils lay symmetrically around the fiber axis. It is also assumed that a narrow region in the middle of the cell walls is examined so that the effect caused by light scattering at the fiber wall edges can be ignored. Under these conditions, the structure of a pulp fiber can be described according to the two-wall model (Ye, C. and Sundström, M. O. PCT Patent Application (1996), WO9610168 and Ye, C, Sundström, M. O., and Remes, K. (1994), Appl. Opt. 33, 6626-6637) with its opposite cell walls approximated by two identical linear retarders in series with their axes symmetrical around the fiber axis. The two retarders have the same relative retardation Δ, which is proportional to the thickness of the fiber walls and the birefringence of the wall material, and their orientation angles have the same value as the microfibril angle φ of the fiber, but with opposite signs. The properties of a wood or pulp fiber related to the polarized light can be described by using the Mueller-matrix formulation (e.g. refer Theocaris, P. S. and Gdoutos, E. E. (1979), Matrix Theory of Photoelasticity, Springer-Verlag Berlin, and Kliger, D. S., Lewis, J. W. and Randall, C. E. (1990), Polarized light in optics and spectroscopy, Academic Press, Harcourt Brace Jovanovich). As calculated (Ye, C, Räty, J., Nyblom, I., Hyvärinen, H. and Moss, P. (2001), Nordic Pulp and Paper Research Journal. 16, 143-148), the Mueller matrix T of a wood or pulp fiber having the retardation Δ and microfibril angle φ when it is in an optical system with the fiber's axis oriented at an angle θ related to a chosen reference axis can be expressed by
According to Equations (1-3), the matrix element t44=1−2cos22φsin2Δ is a function of the retardation Δ and MFA φ, but independent of the fiber's orientation angle θ. The method of the present invention uses a circular polariscope (Theocaria, P. S. and Gdoutos, E. E., Matrix Theory of Photoelasticity, Springer-Verlag, New York, 1979, pages 117-123) that allows transmission of the matrix element t44 through and filters out the all other matrix elements so that it is insensitive to the orientation of a fiber and thus can be used for its real-time measurement. The equipment of the present invention uses the optical arrangement of a circular polariscope to acquire the required measurement information from the fiber sample in combination with a proper spectral analysis.
where I0 is the light intensity transmitted by the equipment and detected by the spectral camera 4 when the fiber 14 is absent in sample unit 10. In principle, the condenser 9 can be placed at any position between the light source 1 and the fiber 14 while the objective 11 may be anywhere between the fiber 14 and the beamsplitter 3.
According to the invention, an optimal estimation of φ and Δ can be generated with a proper spectroscopic analysis of the light intensity I=I[Δ(λ),φ], because the MFA φ is a constant parameter, while the retardation Δ is a function of the light wavelength λ as described by Δ=2πd(n2−n1)/λ. The light beam emergent from the polarizer 13 is divided into two component beams by the beamsplitter 3. One of the split beams reaches CCD camera 5 and the other one is scanned by the spectral camera 4. Only for the purpose of measuring a fiber's φ and Δ, it is not necessary to split the light beam into two component beams. The image of the fiber 14 generated by CCD camera 5 is desirable to serve for controlling and monitoring the measurement procedure and it is interfaced or outputted to image-processing unit 6 and can additionally be used for determining the other parameters of the fiber such as the length, width and shape. The spectral camera 4 typically is a line spectral camera, consisting of an ImSpector 15 followed by a CCD camera 16. The ImSpector 15 is a grating-based imaging spectrograph device and it captures a line image of the fiber's polarizing micrograph and disperses light from the line image into spectrums, which are detected by the CCD camera 16, so that a real-time spectroscopic analysis is feasible.
For a practical system, not only the incident light intensity but also the transmission of each optical component and the response of a used detector as well is a function of the light wavelength X. This means that the light intensity I0 transmitted by the empty equipment is also a function of the light wavelength, i.e. I0=I0(λ), which is contributed by the spectral transmission of all the components in the equipment and the spectral response of the detector to be used in addition to the incident light. Thus, it is necessary to compensate the effect caused by the light intensity I0(λ) in order to determine the parameters φ and Δ based on Equation (4). In accordance with the invention, when a fiber segment is measured, a neighboring background image part without fiber is selected and also scanned by the ImSpector 15. The CCD camera 16 detects and outputs both the spectrum I=I[Δ(λ),φ] of the selected fiber segment and the spectrum I0=I0(λ) of the background image part, which is of the same size as that of the selected fiber segment. Because the background image part I0(λ) is near the fiber segment to be measured and it can be used to approximately describe the light intensity transmitted by the equipment and detected by the CCD camera 16 at the position of the fiber segment in case the fiber is absent. The spectrums I[Δ(λ),φ] and I0(λ) are sent to the image-processing unit 6, which is interfaced to a computer, where the spectrums I[Δ(λ),φ] and I0(λ) are digitized. With the help of proper software, the obtained data are further processed in the computer and the spectrum I[Δ(λ),φ] of the fiber 14 is normalized with I0(λ). The normalized spectrum T=I[Δ(λ),φ]/I0(λ) is determined only by the parameters φ and Δ and from Equation (4) it can be calculated as given by
T=1−cos2 2φ sin2 Δ(λ) (5)
The spectral analysis of the present invention is based on Equation (5) and the least-squares principle. A fitting curve is calculated according to Equation (5) and compared with the measured spectrum T. Estimates may be given for the parameters φ and Δ of the fitting curve as their starting values. Then the measurement results or the optimal estimation for φ and Δ can be generated by varying the estimates for φ and Δ in accordance with the least-squares principle until the sum of the squares of the measured spectral transmission T with respect to the fitting curve is minimized.
The solution provided by the present invention for real-time measurement of intact wood or pulp fibers for the microfibril angle φ and the phase retardation Δ mainly has two exclusive features. As described above, the spectrum T=I[Δ(λ),φ]/I0(λ) of a fiber segment is created and measured in a continuous spectral range, containing the data at all wavelengths in the equipment's spectral range, so that the parameters φ and Δ are over-determined. In this way, any ambiguous measurement results possibly occurred when determining φ and Δ due to insufficient spectral data can be avoided and the measurement accuracy can additionally be improved. Another feature of the invention is that it needs acquiring only one image from the sample and thus it dispenses with special equipment for simultaneously creating and detecting multiple images, which are expensive and technically not desirable.
The equipment of the invention works as a spectroscopic imaging ellipsometer without moving part and it measures a pulp fiber oriented arbitrarily by acquiring only one image from the fiber. Due to this feature, the equipment can be employed for measuring moving fibers if a high-speed CCD camera is used. In addition, it is feasible to measure several fibers simultaneously. It is obvious that the method and equipment of the invention can further be used or adapted to be used for measurement of pulp fibers under on-line conditions. For this purpose, the sample unit 10 of the equipment can be a capillary (e.g. Kajaani FS-200) or a flowing cuvette. The fibers to be measured will be guided for sequentially passing through the capillary or flowing cuvette for measurement.
The equipment of the invention can be modified for real-time or on-line measurement of other birefringent samples such as retardation films and waveplates in addition to other birefringent particles. Because a retardation film or waveplate is not a micro sample as wood fibers, the condenser 9 and the objective 11 are not necessary and thus can be removed. However due to this reason, the method described above for determining the light intensity I0(λ) or I0 cannot directly be applied. A retardation film or waveplate can be considered as a special fiber sample with an imaginary microfibril angle φ=0° and the retardation of its cell walls equal to half the retardation of the retardation film or waveplate.
I
c
=I
0 sin2 Δ(λ), (6)
which is complementary to I[Δ(λ)] such that I[Δ(λ)]+Ic[Δ(λ)]=I0(λ). The spectrums I[Δ(λ)] and I0(λ) are sent to the image-processing unit 6, where the spectrum I[Δ(λ)] of the sample 18 is normalized with I0(λ) and the normalized spectrum T=I[Δ(λ)]/(I[Δ(λ)]+Ic[Δ(λ)]) is used for determining the retardation of sample 18 according to the method of least-squares principle described above.
To test the present invention, a polarizing microscope (Leica DM RX HC, Leica Microsystems Wetzlar GmbH, Wetzlar Germany) was reconstructed with additional optical components added in accordance with the arrangement in
Experiments and measurements were carried out, in which single pulp fibers were measured by using the constructed equipment. To better test the method and the equipment's real-time measurement capability, a fiber for measurement was repeatedly measured after it was rotated to different orientations. A fiber for measurement was first oriented with the fiber's axis parallel or approximately parallel to the polarizers' axes (θ=0°), i.e. the reference axis, and it was then rotated to new orientations from θ=0° to θ=180° with an increment of 22.5°. The present invention will be explained in more detail with reference to the following Examples, which are a small part of the results obtained in the measurements.
The first example is a measured pine kraft pulp fiber. As an example,
The fiber of
The second example was a birch kraft pulp fiber. As an example,
The measurement results of the fiber segment 35 in
Number | Date | Country | Kind |
---|---|---|---|
20060715 | Aug 2006 | FI | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/FI2007/000193 | 7/24/2007 | WO | 00 | 7/20/2009 |