EP1015682A1

EP1015682A1 - A method for the control of a groundwood pulping process

Info

Publication number: EP1015682A1
Application number: EP97928285A
Authority: EP
Inventors: Tom Forsman
Original assignee: Individual
Current assignee: Individual
Priority date: 1996-06-25
Filing date: 1997-06-24
Publication date: 2000-07-05
Anticipated expiration: 2017-06-24
Also published as: CA2294424A1; FI962626A0; WO1997049857A1; FI102975B; FI102975B1; DE69718683D1; FI962626A; EP1015682B1

Abstract

The invention relates to a method for the control of a groundwood pulping process. According to the invention the drainability (freeness, CF) of the pulp and another quantity Q characterising the pulp quality are measured, whereby the quantity Q preferably is the tearing resistance (tear index RI), and the measured values CFx and Qx are compared with the set points CF0 and Q0 of the corresponding quantities. The wood supply rate Vn or the wood supply pressure Fn, and the grinding stone's peripheral speed Vp are adjusted so that the sum (CFx - CF0)2 + (Qx - Q0)2 obtains its minimum value.

Description

A METHOD FOR THE CONTROL OF A GROUNDWOOD PULPING PROCESS

The present invention relates to a method for the control of a groundwood pulping process in order to achieve an optimal value for both the drainability of the pulp and for another characteristic of the pulp, preferably for the tearing resistance of the pulp.

In controlling the pulp grinding process one object is usually to have a constant drainability value or freeness (CF) of the pulp. The control is for instance made so that the wood supply pressure is kept constant, whereby the wood supply rate is allowed to vary. Alternatively the wood supply rate can be kept constant and the supply pressure is allowed to vary.

When only the CF value of the pulp is used as the measured variable to control the process this of course has a disadvantage in that the CF value will not provide all information about the other quality properties of the pulp, which can be characterised by many measured quantities, such as tearing resistance and tensile strength, light- scattering and opacity.

The Finnish patent FI 70438 proposes a method to control a groundwood pulp process with the aid of a new quantity, the plasticity of the wood, as the control parameter. A desired pulp property is obtained at a given (constant) peripheral speed of the grinding stone when the supply pressure and the wood supply rate is selected so that during otherwise constant operating conditions (constant wood quality, constant peripheral speed and sharpness of the grindstone) a plasticity value is obtained.

From tests which are partly published and which are summarised below, it is known that at a constant freeness it is possible to improve the strength characteristics of the pulp, particularly the tearing resistance, by reducing the peripheral speed of the grinding stone. According to an article by Jan-Anders Fagerhed, "Development of wood grinding", Paperi ja Puu - Paper and Timber 72 (1990) :7, the tearing resistance increases about 40 % at a grinding overpressure of 0 to 1 bar when the peripheral speed of the grinding stone is reduced from 30 m/s to 10 m/s. Correspondingly, the tearing resistance increases about 20 % at an overpressure of 2 bar, and about 8 % at an overpressure of 3 to 4 bar. The same article also discloses that the tensile strength (at a given freeness value) can be affected to a certain amount by the peripheral speed of the grinding stone, even if the effect is not as obvious as concerning the tearing resistance. However, the tensile strength increases about 35 % when the grinding is made at atmospheric pressure and the peripheral speed is reduced from 30 m/s to 10 m/s.

In the method presented according to FI 70438 there was not proposed any variation of the peripheral speed in order to obtain an improved tearing resistance in addition to the desired freeness value.

The object of the present invention is to control the pulping process so that an optimal pulp quality is obtained, in other words so that optimal values are obtained both for the CF value of the pulp and for another quantity characterising the quality of the pulp, such as the tearing resistance, which usually is stated as the tear index (RI). As a criterion one uses the minimum sum of the squares of the system deviation from the desired levels concerning these quantities.

The features of the invention are presented in claim 1.

The method can be used in common stone pulping without overpressure (so called stone groundwood or SGW pulp) as well as in so called overpressure pulping (pressure groundwood or PGW) .

In principle the pulping process can be controlled by two control variables, i.a. the wood supply rate (or power) and the peripheral speed of the grinding stone. The supply rate can keep the CF value of the pulp at a desired level, and the peripheral speed of the stone can keep another variable at a desired level. Thus it is possible to control the process by a multivariable method with two input signals and two output signals.

The control can be effected with the aid of a multivariable control algorithm or with two SISO loops (single input, single output) .

The CF value and the tear index of the pulp are kept on a desired level, and the sum of the deviations

(CF_X-CF₀)² + (RI_X-RI₀)²

is minimised, where CF₀ = freeness set point; CF_X = measured freeness value; RI₀ = tear index set point; and RI_X = measured value of the tear index.

The multivariable control algorithm can also be made adaptive in order to compensate for changes in the grinding stone* s sharpness with time.

The relation between the grinding stone' s sharpness and the properties of the mass has been earlier published (see for instance Georg v. Alftan, "Valmistusolojen vaikutus mekaanisen massan ominaisuuksiin" , in the textbook "Puukemia" , Waldemar Jensen, Helsinki 1967.

Measurement data which has been published by Jan-Anders Fagerhed (Development of wood grinding, Part 3 Effects of casing pressure and pulpstone speed, Paperi-Puu - Paper and Timber 72 (1990) :7, 680 - 686) and which is supplemented by previously unpublished material are presented below.

A list of the symbols used below:

m = mass flow kg/h

P = grinding overpressure bar

F_n = supply pressure N

V_n - supply rate rran/s

V_p = peripheral speed m/s

SER - specific energy requirement MWh/t

Tear = tear index mmNNrm /g

CFS = Canadian Standard Freeness ml

Results:

Table 1 : Po ' T = 80 °C +/- - 1 °C m P F„ v_n ^VP SER Tear_ CFS kg/h bar N mm/s m/s MWh/t mNm /g ml

0.97 0 180 0.56 30.0 1.90 2.90 68

1.97 0 200 0.71 30.0 1.52 3.00 120

1.60 0 265 0.85 30.0 1.37 2.80 146

1.85 0 240 1.05 30.0 1.26 2.90 157

0.84 0 185 0.56 20.0 1.58 3.85 75

1.17 0 320 0.64 20.0 1.38 3.80 110

1.47 0 290 0.80 20.0 1.23 3.40 110

1.57 0 355 0.92 20.0 1.07 3.15 180

0.66 0 280 0.36 10.1 1.44 3.75 90

0.92 0 380 0.50 10.0 1.29 4.20 100

1.12 0 500 0.58 9.9 1.14 4.35 150

1.23 0 465 0.69 10.0 1.01 4.20 170 Table 2 : P_: T = 95 ^UC +/ - 1 C m P F_n v_n V_P SER rear CFS kg/h bar N mm/s m/s MWh/t mNm /g ml

0.99 1.0 110 0.41 30.0 1.79 3.70 90

1.07 1.0 170 0.53 30.0 1.84 3.90 65

1.28 1.0 200 0.63 30.0 1.55 3.85 105

1.50 1.0 225 0.74 30.0 1.40 3.25 120

0.75 1.0 150 0.38 20.0 1.57 4.65 90

1.00 1.0 245 0.48 20.0 1.45 4.40 85

1.28 1.0 265 0.59 20.1 1.15 5.15 140

1.34 1.0 230 0.69 20.0 1.31 4.55 60

0.64 1.0 335 0.30 10.0 1.38 5.35 85

0.79 1.0 420 0.38 10.0 1.02 4.95 95

1.04 1.0 435 0.49 10.0 1.09 5.30 110

1.18 1.0 460 0.59 10.0 0.93 5.45 120

Table 3: P₂ T = 110 °c +/- 1 °C m P F_n v_n ^VP SER rear ( ZFS kg/h bar N mm/s m/s MWh/t mNm /g ml

0.94 2.0 110 0.51 30.0 1.61 4.55 120

1.28 2.0 210 0.62 30.0 1.39 5.05 130

1.66 2.0 200 0.76 30.0 1.06 4.80 220

1.88 2.0 195 0.94 30.0 1.18 4.50 175

0.81 2.0 80 0.41 20.0 1.34 5.40 100

0.88 2.0 210 0.51 20.0 1.20 5.10 145

1.35 2.0 310 0.61 20.0 1.45 5.25 135

1.44 2.0 220 0.69 20.0 1.67 4.70 95

0.57 2.0 285 0.28 10.0 1.44 5.85 75

0.73 2.0 355 0.38 10.0 1.24 5.55 160

1.01 2.0 425 0.49 9.9 1.09 5.10 195

1.21 1.9 475 0.59 10.0 0.95 6.05 255 Tabl<2 4: P₃ T = :L20 C +/- 1 °C m P F„ v„ v_P SER Tear CFS kg/h bar N mm/s m/s MWh/t ; mNm /qf ml

0.76 3.0 75 0.40 30.0 1.67 5.35 75

1.01 3.0 135 0.50 30.0 1.39 5.25 105

1.26 3.0 150 0.60 30.0 1.20 5.45 100

1.48 3.0 155 0.72 30.0 1.24 5.75 100

0.74 3.0 130 0.35 20.0 1.30 5.90 100

0.94 3.0 250 0.45 20.0 1.42 5.55 60

1.10 3.0 255 0.56 20.0 1.45 5.85 70

1.29 3.0 225 0.67 20.0 1.12 5.75 140

0.58 3.0 310 0.28 10.0 1.52 6.00 100

0.70 3.0 350 0.36 10.0 1.40 5.65 115

0.89 3.0 420 0.46 10.0 1.19 5.80 175

1.05 3.0 480 0.54 10.0 1.19 6.45 150

Table 5: P, T = 130 °C +/- 1 °C m P F_n V_n V_p SER Tear i CSF kg/h bar N mm/s m/s MWh/t mNm /g ml

0.77 4.0 95 0.40 30.0 1.71 5.35 70

0.95 4.0 120 0.50 30.1 1.60 4.95 65

1.05 4.0 145 0.58 30.0 1.25 5.30 120

1.26 4.0 165 0.67 30.1 1.09 5.00 155

0.64 4.0 120 0.33 20.0 1.06 5.75 130

0.81 4.0 205 0.42 20.0 1.45 5.60 85

1.00 4.0 185 0.52 20.0 1.35 5.50 100

1.23 4.0 190 0.62 20.0 1.11 5.45 135

0.48 4.0 265 0.25 10.0 1.61 5.55 80

0.60 4.0 365 0.33 10.0 1.34 5.40 155

0.79 4.0 375 0.42 10.0 0.22 6.10 180

1.01 4.0 385 0.53 10.0 0.97 5.90 230

Table 6: P₅ T = 140 °C m P F_n v„ ^VP SER rear CSF kg/h bar N mm/s m/s MWh/t mNm /g ml

0.80 5.0 175 0.39 30.1 1.64 5.30 80

0.98 5.0 165 0.50 30.0 1.28 5.70 95

1.21 5.0 125 0.60 30.0 1.02 5.40 215

1.29 5.0 160 0.69 30.1 1.19 5.75 125

0.70 5.0 180 0.33 20.0 1.56 5.65 65

0.85 5.0 140 0.42 20.0 1.13 5.35 120

0.93 5.0 155 0.51 20.0 1.19 5.70 120

1.19 5.0 225 0.60 20.0 1.03 5.35 145

0.45 5.0 215 0.25 10.0 1.51 5.65 65

0.62 5.0 320 0.32 10.0 1.41 6.45 150

0.41 5.0 210 0.21 10.0 1.49 4.85 75

0.77 5.0 270 0.42 10.0 1.11 6.10 210

The relation between quantities characterising the pulp properties (freeness, tear index) and the operating conditions of the process can be determined by regression analysis based on the measurement data presented above.

The results show that the mass flow can be kept rather constant despite the lower peripheral speeds because the supply pressure is increased.

The method according to the invention also reduces the specific energy consumption (SER).

Control methodics:

An adaptive (self-adjusting) control algorithm is presented below. The controller is a generalisation of the multivariable control algorithm of Astrδm and Wittenmark (1973) .

The process can be described by the equation below: y(t) + A_lY(t-l) + ... + A„y(t-n) =

= B₀u(t-k-l) + ... + B_n^uft-k-n) + e(t) +

+ Cιβ(t-l) + ... C_ne(t-n) (1)

where u is the input vector and y is the output vector, and {e(t)} is a sequence of independent evenly distributed random vectors with a mean value of zero and the covariance

E[e(t)e^τ(t) ] = R

The dimension of all vectors u, y and e is p, and the dimension of all matrices A_t, Bj. and C*_. is pxp. The matrix B₀ is non-singular.

Now we introduce the shift operator q^"1 defined as

q^_1(t) = y(t-l)

and the polynomial matrices

A(z) = I + AiZ + ... + A„zⁿ n-l B B((zz)) == BB₀₀ ++ BBiiZZ ++ ... + B B_nn_-_z

It is assumed that all zeros of B(z) are outside the unit circle. B₀ is non-singular. The system (1) can be written as

A(q^_1)y(t) - B(q-¹)u(t-k-l) + C(q^'1)e(t) (2)

In each sampling interval the adaptive algorithm performs an identification based on the least squares method according to the model presented below.

The obtained parameters are used for calculation of the control strategy. 7/49857 PC17FI97/00406

Estimation

The algorithm estimates the parameters for the model

y(t) + A(q-^χ)y(t) = Bfq'¹ )y(t-k-l) + 6(t) (3)

so that the error t=(t) is minimised according to the least squares.

In the model (3) k is selected as the dead time for the process (2), and A(z) and B(z) are pxp polynomial matrices according to

S(z) = B₀ + BjZ + ... + B_nBz^nB

First we assume that

B₀ = I and

B₀ = I

where B₀ is a matrix in the constant term of B(z) for the process (2) .

Now we introduce the column vectors

σ_t = [Ct_j_! . . .vϊ_j_p . ..OJ_j__j . . .Ot_j_p

nfl

>il Hp •βi ^•PLP J / i - l/ ^■P (4)

where α^ is the (i,j)^c element in the matrix A_k; βi_j ^k is the (ijj)¹¹ element in the matrix B_k, and so on. Then the column vector θi can be considered to contain the coefficients of the i^th row in the model (3). Further we introduce the row vector

φ(t-k-l) * [-y^x(t-k-l) -y^t-k-l-n*)

u (t-k-2) u (t-k-l-n_B) (5)

The i^th row in model (3) can be written as

€(t) = y_t(t) - Ui(t-k-l) - φ(t-k-l)θj.

According to the least squares criterion the vector θj_. at each moment N is calculated so that

V„(θ_t) = 1/N Σ €_t (t), i = 1, ... , p (6) i'\ is minimised. This results in a least squares estimation of each row in (2) based on data which is available at the moment N. When N is large, the initial values are of insignificant importance in (6). The criterion (6) can be written as

V_H(θ_t) = 1/N Σ Σ _yi(t) - Ui(t-k-l) - φ(t-k-l)θi]²

i = 1, ..., p (7)

The value θi which minimises (7) is given by the normal equations, see Astrδm and Eykhoff (1971).

A/ _τ t Σ φ (t-k-l)φ(t-k-l) ]θi(N) =

= Σ φ^τ(t-k-l)[y_t(t)-Ui(t-k-l)]

i = 1, ... p (8)

Control

At each moment t the control strategy is calculated from -1 -l ,

S( q ) u ( t ) = A( q ^* )y( t ) ( 9 )

where the polynomial matrices A(z) and B( z ) are obtained from the current value of the estimated parameters. The control strategy can be written as

u_t(t) = -φ(t)θ_t(t) i = 1, ..., p (10)

The parameters for the controller are the same as the estimated parameters. When we use

the strategy (10) can be written as

u^τ(t) = -φ(t)θ(t) (12)

✓v.

The estimated parameter vector β_L in (8) can be recursively calculated from

θ_i(t)=θ_i(t-l) + K(t-l)[y_i(t)-u_i(t-k-l)-φ(t-k-l)θ_i(t-l)]

K(t-l)=P(t-l)φ^T(t-k-l)[l+φ(t-k-l)P(t-l)φ^T(t-k-l)]^'1 (13)

P(t)=P(t-l)-K(t-l)[l+φ(t-k-l)P(t-l)φ^τ(t-k-l)]K^T(t-l)

P(t) is a normalised covariance matrix of the estimated parameters Q_t .

The initial values of P(t) are assumed to be the same for all parameter vectors θj,. The corresponding amplification vectors K(t-l) will also be the same for all estimators.

Sometimes it may be useful to introduce an exponential weighting for the parameter estimation. This can be done by modifying the criterion (6) to rV

Σ λ" ^ce^( t ) i = 1 , λ < 1 ( 14 ) f1

The last equation in (13) changes to

P(t)= l/λ{P(t-l)-K(t-l)[l+φ(t-k-l)P(t-l)φ^x(t-k-l)] x K^τ(t-1)] (15)

Another possibility is to use Kalman filters. The covariance matrix P(t) is supplemented by adding to it a matrix Rι instead of the division by λ.

Then the equation (15) will be

P(t) = P(t-l)-K(t-l)[l+φ(t-k-l)P(t-l)φ^T(t-k-l)] x K^T(t-l) + R.

It should be noted that the algorithm can be construed as a union of a plurality (here 2) of simple self-adjusting controllers. For instance the controller 2 controls the output signal y₂(t) by using the control variable u₂(t). yj(t-i) and u^t-l-i) (i > 0) can be used as feedforward signals. This means that the two simple self-adjusting controllers can operate in a cascade mode.

The possibilities for this feature strongly depend on the process properties regarding the model (2) and character of the minimum variance strategy. The multivariable self- adjusting control algorithm can in some circumstances result in the minimum variance, in other words when C(z) = I (the process interference is white noise).

Another possibility is an exclusively multivariable minimum variance control algorithm, which is not adaptive.

At a pulping overpressure of 0 to 2 bar the control of the tear index at lower peripheral speeds results in great advantages (40 % to 20 %). As the multivariable control algorithm also is adaptive, changed sharpness is taken into account by increasing the peripheral speed. During this the freeness can be freely selected.

At higher pulping overpressures the advantage is an improvement of about 10 % concerning the tear index, and the changes in sharpness can be controlled in the periods between sharpening actions. During these periods the freeness can be freely selected.

If the sharpening is not made with pressurised water or similar at regular intervals, then the sharpening is made at P-^ at the maximum power consumption.

Claims

1. A method for the control of a groundwood pulping process, whereby pulpwood logs are pressed against the periphery of a rotating grinding stone, the grinding stone is sprayed with water, and the generated fiber suspension, the pulp, is stored, characterised in that the drainability or the freeness CF of the pulp and another quantity Q characterising the pulp quality are measured, the measured values CF_X and Q_x are compared with the set points CF₀ and Q₀ of the corresponding quantities, and the wood supply rate V_n or the wood supply pressure F_n, and the grinding stone' s peripheral speed V_p are adjusted so that the sum (CF_X - CF₀)² + (Q_x - Qo)² obtains its minimum value.

2. A method according to claim 1, characterised in that the quantity Q is a measure of the tearing resistance of the pulp, for instance the tear index RI.

3. A method according to claim 1 or 2, characterised in that the control is effected with the aid of a multivariable control algorithm.

4. A method according to claim 3, characterised in that the multivariable control algorithm is adaptive in order to compensate for changes in the grinding stone' s sharpness with time.