In-Situ Frost Heave Loads in Artificially Frozen Ground For Tunnelling
In-Situ Frost Heave Loads in Artificially Frozen Ground For Tunnelling
In-Situ Frost Heave Loads in Artificially Frozen Ground For Tunnelling
R.H.B. Rijkers TNO-NITG / Netherlands Institute of Applied Geosciences B. Hemmen GeoDelft N.M. Naaktgeboren Ministry of Transport and Public Works H. Weigl Kombinatie Middelplaat Westerschelde ABSTRACT: In the Netherlands the Westerscheldetunnel will be constructed with two tunnel bore machines and 26 connecting cross passages. The cross passages are built between the two main bored tunnels by using artificial ground freezing. During the construction phase of the cross passages an extensive monitoring program is carried out: (1) ground stress conditions and ground deformation and (2) deformation of the tunnel constructions by measuring temperatures, frost heave loads and deformation. Measurements are taken in three different orientations: perpendicular, parallel and vertical to the cross passage with stress -monitoring stations and extensometers. Significant differences have been observed in loads and deformations, which are explained by different soil conditions. A significant variation in ground stresses has been monitored between the perpendicular and parallel direction with respect to the orientation of the freezing tubes. All separate phases of the construction activities of the c ross passages (freezing, excavation, lining, thawing) are recognised clearly in the data. 1 INTRODUCTION AND OBJECTIVES For the construction of the Westerscheldetunnel in the Netherlands 26 cross passages (or connections) are built between two parallel bored tunnels with artificial ground freezing (AGF) techniques (Hass et al., 2000). Special in -situ measurements of stress and deformation are carried out during the freezing period for the first two cross passages (DV1 and DV2) in order to evaluate the effects of the ground freezing technique. Due to the increase of volume of water to ice with 9% and ice lens formation frost heave loads are expected. The monitoring program includes measurements of stress and deformation development in the soil and the tunnel constructions. The objectives of this study are (1) in -situ soil monitoring during period of artificial soil freezing , (2) evaluation of the frost heave loads to the tunnel due to soil freezing and (3) increase of knowledge on use of ground freezing techniques in soft soils for tunnelling. 2 WESTERSCHELDETUNNEL 2.1 Artificial ground freezing The Westerscheldetunnel is located in the southwestern part of the Netherlands below the estuarine river that connects to the port of Antwerpen (Belgium). In order to increase safety at the Westerscheldetunnel the two lanes in each direction are designed as separate tunnel tubes. These two parallel tunnel tubes are bored with inner diameters of 10.10m. The total length of each bored tunnel tube is 6.6km. To increase safety in the tunnel two bored tubes are connected with 26 cross connections at every 250m. The cross passages are built from the inner side of the main tunnel tubes by using artificial ground freezing in order to make underground excavations possible. Artificial ground freezing techniques are used for both sandy soils and Boom clay because stiffness and permeability of all soil units in the trajectory is not considered to be sufficient for open excavation . 2.2 Construction of cross passages The 26 cross passages will be constructed in different ground conditions. 9 cross passages are made in sand, 11 in the stiff clay and 6 in mixed (sandy) soils. The cross passage has an elliptical reinforced concrete lining with an thickness of 400mm, height of 2.75m and width of 2.5m. Outside these lining there is temporary shotcrete of 300mm thickness. In the main tunnel tubes five specials rings have been designed with a previously fixed and unchangeable sequence. In the middle ring steel segments replace two concrete segments. In these special enlarged steel segment the safety doors to the cross connection are constructed (Fig. 1). 22 freezing pipes are bored trough special watertight
connections in the specially prepared ring segments (Fig. 2). On completion of soil freezing (a minimum thickness of 2m frozen ground) the excavation is started by the NATM-method. During the last phase the concrete lining of the cross passage is built and the freezing machine is stopped eventually for thawing the soil. During the excavation period the emergency door can be closed when an unsafe situation occurs.
Frozen ring
Train
order to avoid high pressure within the frozen body after closure of the frozen body (cylindrical formed), a draining possibility was installed near the axis of the cross passage. With this draining system water pressure in the ice-enclosed area was released in a controlled manner (an increase of water pressure was measured in sandy soil at DV1). In order to safeguard a watertight connection between ice body and the tubings, approx. 26 temperature sensors, distributed in a plane, are distributed in the western tube. With these sensors it was possible to monitor the ice connection by temperature of the frozen soil. In addition, a 10cm thick insulation is installed at the interior side of the tunnel tube. This insulation provided a faster connection of the ice body and warrants a longer remaining ice connection in case of a failure of Body of frozen soil approx. 2m the freezing system.
Kabelkanaal
2.2 Ground freezing machine Ground freezing is performed as brine glaciation (NaCl2) since a nitrogen glaciation in the tunnel was excluded for safety reasons. A freezing unit was developed for soil freezing at temperatures of the freezing tubes at approx. 38C. This unit has a min. ammonia content of approx. 80kg for a performance of 95kW and numerous safety features. In a circular pattern 22 freezing tubes are bored with a distance of about 1.0m to each other and parallel to the cross way axis (Fig. 2). The freezin g tubes are connected to the main pipe by means of a flexible pipe. The complete piping system includes approx. 3m brine. When starting the soil freezing system, the brine temperature at every brine head on the starting side is measured in order to check whether the freezing tube works properly and whether the freezing capacity is distributed evenly. 2.3 Safety valve The spatial and temporal development of the complete freezing body was checked permanently. The freezing body has a target thickness of > 2.00m which was reached after approx. 26 days in sand at DV1 and after approx. 46 days in clay at DV2. Two diagonally arranged temperature monitoring devices controlled the formation of the frozen soil body. In
outer lining = shotcrete (250 mm) inner lining = concrete cast in -situ
= freezing tube
3 SOIL CONDITIONS 3.1 Soil profile and geotechnical characteristics Soil freezing was carried out at cross passages DV1 and DV2 in sand and clay (Table 1). Figure 3 gives a soil profile based on several geological borings, cone penetration tests and geophysical bore hole measurements. Soils of three geotechnical units will be frozen for the first two cross passages: Z1, BK1
and BK2. Boom clay units (BK1 and BK2) dip slightly towards the north and are covered with young coastal sands (Z1). Note the large difference in plasticity and permeability between soil types at DV1 and DV2 (Table 2). The top of the Boom clay varies strongly due to erosion (Figure 3). The axis off cross passage DV1 is at 17,5m-NAP and is constructed in loose sand of Holocene age (Z1 unit). The axis of cross passage DV2 is located at a depth of 28,7mN A P a n d i s c o n s t r u c t e d i n o v e rconsolidated clay of Oligocene age (BK1 and BK2
Table 1
unit). Both the tops of the BK2 unit and the Z1 unit are characterised with a salt contents of 2000-13000 Clmg/l. At both cross passages no significant ground water flow was recognised. 3.2 Frost susceptibility Accordingly international standards of ISSMFE frost heave susceptibility of unit H1 (sand) is classified as negligible and the overconsolidated clay of unit BK1 and BK2 as medium to strong Frost heave susceptibility is the tendency of the soil to expand during soil freezing - due to the volumetrical
Soil properties of geotechnical units Z1, BK1 and BK2 Z1 BK1 BK2
soil properties
lithology clay contents ( mean) sand contents ( mean) water contents w wet dry porosity p permeability k vertical permeability k horizontal k vert . / k hor. W p plasticity limit W l liquidity limit I p plasticity index shear strength torvane salt contents
silty CLAY 81.2 18.8 25.3 19.4 ~ 50 1.5E-9 2.4E-11 2.7E-10 5.8E-8 ~ 30 30 91 61 ~ 220 2000-13000
silty CLAY
with sand layers
Table 3 In-situ total stress and water pressure u at cross passages DV1 and DV2 measured by SM1 and SM4 (stress monitoring stations) before artificial ground freezing. The values between brackets are theoretically values based on depth and volumetrical unit weights of the soil column; K 0 =0.5.
vertical [kN/m 2 ] cross passage DV1 (SM 1) depth axis at 17,5 m-NAP cross passage DV2 (SM 4) depth axis at 28,7 m-NAP 204 (309) 310 (607)
expansion of water to ice with 9% - and the growth of ice lenses. Note the difference in BKJ1 and BK2 in clay contents. 3.3 Natural stress conditions In-situ stress conditions of unfrozen soil have been measured with stress monitoring stations and have been compared with calculated theoretically soil stresses (Table 3). Calculated stress conditions are based on soil columns and volumetrical unit weights. Measured and calculated horizontal
stresses match very well (Table 3). The measured vertical stresses are significantly lower than calculated values, probably due to the installation procedure of the monitoring equipment and remoulding of the soil.
Figure 3
Soil profile with the position and depth of the cross passage DV1 and DV2. Soils of units Z1, BK1 and BK2 have been frozen during the construction of cross passages DV1 and DV2.
4 GEOMECHANICAL BEHAVIOUR OF FROZEN SOIL In sand the water contents shall freeze almost instantaneously due to the high percentage of free water. Clay shows a different behaviour in frozen state due to the strongly bonding forces between water and clay minerals. In clay a percentage of water will be unfrozen at temperatures of 35C. The grain size distribution also controls consequent differences in frost heave and strength development. 4.1 Frost heave and ice lenses In freezing soil the formation water can segregate into ice lenses before it is frozen (Konrad & Morgenstern, 1981). As the result of growing ice lenses a water pressure gradient will develop that is depending on grain size and permeability of the soil (Penner, 1986). The thickness of ice lenses can range from < 10m to several cm. The freezing expansion of clay is partly caused by ice lenses. The growth of the ice lenses is enhanced by the (slow) migration of water at this frozen front. That is the reason why soils with low permeability such as silts and clays are frost heave susceptible and (drained) sand is not. The process of migration of water due to freezing temperatures is referred as cryosuction. Konrad & Morgenstern (1981) have defined experimentally a linear relation between the frost heave rate h of a soil, the segregation potential SP 0 and the temperature gradient T:
h = SP 0 . grad T
Figure 4 Frost heave rate vs temperature gradient of BK1 / BK2 soil samples.
Three-dimensional frost heave behaviour of Boom clay has been experimentally investigated and reported (Rijkers et al, 2000). 4.2 Creep Creep theoretically decreases frost heave stresses that have developed in an initial stage in a groundfreezing project. Assur (1963) states that creep of frozen soil is defined by temperature, time, load and material characteristics of soil. Creep processes are expected to occur in the frozen soils bodies at DV1 and DV2. 5 MONITORING PROGRAM 5.1 Monitoring equipment In order to observe changes in stress, temperature and deformation a special monitoring program was designed. This program consists of water pressure recorders, stress monitoring stations, spade cells, temperature recorders, e x t e n s o m e t e r s a n d inclinometers. The stress monitoring stations consist of three flat plates that record water pressure and temperature. The three plates can measure stresses directed parallel, perpendicular and vertically with respect to the freezing tubes. The inclinometers have also been installed, but results are not discussed in this paper.
This linear relation between temperature gradient and frost heave rate is experimentally been established with samples of the Boom clay from the tunnel trajectory (Fig. 4). Resulting SP0values smaller than 0.5m2/C resulted in a low category of frost susceptibility. Following experiences in civil engineering with AGF works frost heaves were expected in the direction perpendicular on the freezing front or isotherms and not in the direction parallel to the axis the freezing tubes.
Table 4 Characteristics and orientation of monitoring data of cross passage DV1 and DV2. cross passage DV1 data type
tangential stress axial stress radial stress
orientation
vertical parallel to freezing tubes perpendicular to freezing tubes vertical
station
SM1 SM1 SM1 SM1 EX2 VW1
Extensometers are located at the top of the each cross passage and consist of four packers at different depth levels. The packers are fixed in the soil and connected to the ground surface with bars. At the surface the shortening of the bars (distance between the packer and the surface) is measured. Water pressure recorders are situated on most of the instrument. During the complete period of ground freezing and thawing at cross passage DV1 and DV2 data was recorded. Due to the large amount of recorded data and different data types, it is not in the scope of this paper to address all measurements. We restrict this paper to the highlights of stress monitoring stations (SM1 and SM4; Figs 6-7) and extensometers (EX2 and EX3; Figs 8-9) at cross connections DV1 and DV2. In table 4 the reported data types are summarised. 6 STRESS AND DEFORMATION 6.1 Monitoring results Data records of stress monitoring, extensometers and temperature data are given as time series in Figures 6-9. Stress is given in these graphs as absolute values. In the upper parts of the graphs main construction phases of the cross connections are given, such as freezing, excavation and covering. 6.2 Stress deviations from SM1 at DV1 When the freezing process starts the frozen front has not immediately reached the monitoring instruments (Fig. 6). Radial and axial stresses in the surrounding unfrozen soil (sand) decrease. From the moment where the freezing front reaches the monitoring instrument (temperature < 0 C), the instrument is part of the frozen soil. The frozen soil is expanding due to the continuing freezing process. Inside the frozen soil radial, axial and tangential stresses increase. Consequently the cylindrical frozen soil body is tightening. The first construction phase is to excavate the heart of the cylindrical frozen soil body. Support to the frozen soil from inside is then removed. Axial and tangential stresses increase, while radial stress decreases. Shortly after excavations the frozen soil rises slightly in temperature, because the excavated area is in contact with relatively hot open air and hydration heat of the shotcrete. Temperature alternations cause deviations in the stiffness of the frozen soil and other equilibrium of stress distribution over the total cylindrical frozen soil. It leads to axial and tangential stress increase and radial stress decrease at the point of the stress monitoring stations. The freezing process is continui n g a n d t h e temperature in the total cylindrical frozen soil body decreases further. These actions provide radial and
axial support from the inside. Radial soil stress increases; the frozen soil reacts against the shotcrete. Axial and tangential stresses decrease. The freezing apparatus starts working with a lower capacity (33%). Temperature is increasing all over the frozen soil in a fast rate. Stress deviations also occur at a fast rate and stresses are redistributed over the frozen soil. The decrease in stiffness (due to temperature rise) leads to a decrease in radial and axial stresses. An overall radial stress release leads to an overall tangential stress increase. After the fast increase the temperature stabilises causing a stabilisation of all stresses. Finally, ground freezing stops and temperatures start to increase in the direction of the freezing point. The same process as freezing with a lower capacity recurs. There is one difference. Radial stress increases immediately after shutting the freezin g apparatus. This reaction is probably caused by the redistribution of stresses over the frozen soil, with this difference that the total distribution of stiffness over the frozen soil is different as it was during the previous period of temperature increase. The final tangential stress increased with respect to the starting value. This is a result of the installation procedure of stress monitoring stations. The vertical stress after installation is less as the original vertical in-situ stress (Table 3). Probably because soil has been excavated via a pulse boring to install this type of instrument. 6.3 Water pressure deviations from SM1 at DV1 There are no large water pressure deviations while freezing in sand. The maximum amount of water pressure decrease is approximately -15kPa. The sand is well drained and the decrease is a result of the chase of water during volume increase of frozen water. 6.4 Soil deformations from EX2 at DV1 In order to make installation of the instruments possible, soil is improved at both cross passages at ground level. Besides an extra cinder level is deposited. The extra load causes settlement to occur. In the beginning before the period of freezing this settlement is visible as a shortening of the packer bars. After a while when the soil around the freezing tubes freezes and expands, packer EX2-1 reacts and shortening takes place (Fig. 8). Packer EX2-2 to EX2-4 reacts in a later stadium. The frozen front progresses slowly. In the period of excavation the shortening drops (most for packer EX2-1 and respectively less for packer EX2-2 to EX2-4). The excavation causes a radial release and therefore a contraction of the cylindrical frozen soil body. After this period deviations in deformation are low. Packer EX2-1 to EX2-4 show similar fluctuations
conformable with radial stress observed at SM1. Further decrease in temperature causes a small expansion. The freezing apparatus starts working with a lower capacity. Packer EX2-1 shows a small decrease in shortening followed by an increase. Packers EX2-2 to EX2-4 stabilise and increase. The freezing process continued although at a lower level. The frozen soil is still expanding. Thawing leads to an overall decrease in shortening. The residual shortening is the result of the surface settlement for packer EX2-2 to EX2-4. Packer EX21 has a larger residual settlement. In fact this packer was the only one of EX2 packers inside the frozen soil. 6.5 Stress deviations from SM4 at DV2 All stresses increase strongly from the beginning after a frozen soil body (clay) starts forming (Fig. 7). Also tangential stress increases. The reason for this strong frost heave is the typical undrained character of clay. Water pressure increases with the same amount as the tangential stress. Radial and axial stresses increase conformable expectations, the freezing soil expands causing a stress increase in the surrounding unfrozen soil due to freezing of the water contents of the clay and to ice lens formation. From the moment the temperature drops under the freezing point, water pressure decreases in a fast rate to just above the vacuum pressure. Radial, axial and tangential stresses react and show a slight decrease. The water pressure is not as dominant any more as most of the water starts to freeze. Radial stress increases again to about a level where all three stresses start to decrease. This decrease in stress is due to creep or the volumetric contraction of ice at lower temperatures. Creep occurs at a certain level of temperature and stresses. The crystalline structure of the frozen water is not stable anymore and the soil more or less liquefies. At cross passage DV1 stresses in all principal directions increase in this period. The stress rate is much lower. Excavation leads to a radial release. Axial stress increases for the support from the heart is removed. The cylindrical frozen soil body contracts and tangential stress increases. Just after excavation the periods of cover and freezing with lower capacity begins. Temperature increases for a period of time and decreases again after a while. Radial stress increases slightly while the freezing process is continuing, only at a lower rate. Tangential stress increases fast. These phenomena also take place at cross passage 1 where a temperature increase leads to a fast tangential stress increase, caused by the redistribution of stresses and stiffness deviations in the cylindrical frozen soil body. The axial stress decreases also due to the redistribution of stresses over the in the meantime covered cross connection. Finally,
thawing does the same to the tangential stress as freezing with lower capacity. Tangential stress increases rapidly and due to a combination of creep and thawing the tangential stress is reduced to its original starting level. The axial and vertical stresses decrease by thawing. The total expansion stops and stresses decrease. The water pressure finds its original level at a temperature just below the freezing point. 6.6 Deformations from EX3 at DV2 The same soil improvement and cinder installation takes place at cross passage DV2. The same amount of settlement is visible in the beginning before the freezing period (Fig. 9). From the moment the freezing starts, expansion of the frozen soil body takes place. The difference between cross passage DV1 and cross passage DV2 is clear. The maximum shortening at DV1 is approximately 10 mm at packer EX3-1. The maximum shortening at DV2 is approximately 50 mm at packer EX3-1. Excavation leads to a decrease in shortening. The cylindrical frozen soil body shrinks. From this moment temperature increases and, but freezing continuous at a lower level. A smooth small expansion is the result. Thawing leads to the final decrease in shortening and all packers return to their original situation including the surface settlement. The residual shortening is large compared to DV1. 6.7 Drainage of the frozen soil and insulation measures The average temperature evolution was the same as calculated for both cross passages. However, at the connection between the main tunnel tube and the frozen ring the temperature drops slower than expected. At the steel segments it is clear that the conduction of heat is high that proper insulation was necessary. At cross connection DV1 dry ice was eventually used to maintain low temperatures at the connection between the steel segments and the frozen ring. For drainage during soil freezing a borehole was used to control the water pressure inside the frozen body. Water expands when it is frozen, so in an undrained situation water pressure rises when the frozen ring is watertight (Fig. 10, point B). In the first part the increasing water pressure (Fig. 10, point A) gives an indication that the ring of frozen soil is watertight. In our case at DV1 monitored temperatures through the tunnel segments gives that soil was not frozen directly behind the segments. Another test was conducted with opening and closing the drainage tap (Fig. 10). Because the water didnt stop we concluded that the frozen ring was not watertight and after one day using dry-ice a significant rise of the water pressure occurs that proved the water tightness of the frozen ring.
Stress monitoring station 1 (SM1) Absolute ground- and waterpressure versus Temperature and Date
330
Air release Water release T = 0 C Freezing decreased capacity Cover Excavation Freezing T = 0 C
15
310
10
290 5
Stress [kPa]
250
-SM1-5 temperature
-5
230
190
-15
170
SM1-4 waterpressure
-20
150 4/9 4/23 5/7 5/21 6/4 6/18 7/2 7/16 7/30 8/13 8/27 9/10 9/24 10/8 10/22 11/5 11/19 12/3 12/17 12/31 1/14
-25
Date
15
1400
T = 0 C Cover Excavation
T = 0 C
10
Freezing
1200
1000
0
Temperature = 0 C
400
-SM4-5 temperature
200
-5
SM4-4 waterpressure
0 5/26 6/9 6/23 7/7 7/21 8/4 8/18 9/1 9/15 9/29 10/13 10/27 11/10 11/24 12/8 12/22 1/5 1/19 2/2 2/16 3/2 3/16 3/30 4/13 -10
Date
Temperature [C]
Stress [kPa]
Temperature [C]
Temperature = 0 C
15 Water release
13
T = 0 C Cover Excavation
12
11 5 10
Temperature = 0 C
Shortening [mm]
-VW1-1 temperature
7 -5
EX2-1
-10
EX2-2
3
-15
EX2-3 EX2-4
-20
0 4/9 4/23 5/7 5/21 6/4 6/18 7/2 7/16 7/30 8/13 8/27 9/10 9/24 10/8 10/22 11/5 11/19 12/3 12/17 12/31 1/14
-25
Date
10 -15
0 5/26 6/9 6/23 7/7 7/21 8/4 8/18 9/1 9/15 9/29 10/13 10/27 11/10 11/24 12/8 12/22 1/5 1/19 2/2 2/16 3/2 3/16 3/30 Date
-20 4/13
Temperature [C]
Shortening [mm]
Temperature [C]
2.5
Waterpressure [bar]
1.5
0.5
0 5/6 5/13 5/20 5/27 6/3 date 6/10 6/17 6/24 7/1 7/8
Frost heave pressure is defined as the increase of ground pressure due to soil freezing. The absolute measurements of frost heave pressures at DV1 and DV2 are higher (see Figs. 6-7). At cross passage DV2 (in overconsolidated clay) the maximum offset of 20mm was registered. Also corresponding measured frost heave values are significantly higher. In the direction parallel to the freezing tubes 680kPa of frost heave stress has been observed that eventually caused 20mm offset of the tunnel segments (SM3 SM4 //; Fig. 11). 8 CONCLUSIONS Frost heave loads and soil deformations have successfully been measured in (frozen) soil bodies of sand and clay during an artificial ground freezing activities. Frost heave pressures differ significantly between frozen sand and frozen clay bodies. All phases of construction of the cross passage are recognized clearly in the data. 8. ACKNOWLEDGEMENTS We thank the organizations of KMW and COB for the permission publishing this paper. We also thank Andr Koers and Roel van de Kraan of TNO for preparing the drawings. 9. REFERENCES Assur, A. (1963) Creep of frozen soils. Proc. 1st Permafrost Int. In: Conf. Nat. Acad. Sci. Pucbl. Vol. 1287, pp. 339 Hass, H, R. Jagow-Klaff & H.W. Seidel (2000) Use of brine freezing for the construction of the traverse galleries. In: ISGF2000 - Frost action in soils (ed. J.F. Thimus), Louvain-la-Neuve, pp. 289-294 Kofoed, N. & S.R. Doran (1996) Storbaelt Eastern Railway Tunnel, Denmark ground freezing for tunnels and cross passages. In: Geotechnical Aspects of Underground Construction in Soft ground, Balkema Rotterdam 1996, pp. 391-398 Konrad, J.M. & N.R. Morgenstern (1981) The segregation potential of a freezing soil, Cndn Geotech. J. 18, pp. 482-491 Penner, E. & T. Ueda (1977) The dependence of frost heaving on load application preliminary results. In: Proc. Int. Symp. Frost Action in Soils, Lulea , Vol. 1, pp. 92-101 Rijkers, R, M. Naaktgeboren, H. Ct, J.-F. Thimus, M. Visschedijk (2000) Soil behaviour during artificial freezing Part 2: Application in bored tunnels, In: ISGF2000 - Frost action in soils (ed. J.F. Thimus), Louvain-la-Neuve, pp. 237-242
At the cross passage DV2 it was not possible to drain the inner part of the frozen ring because of Boom clay with very low permeability. However, the drainage pipe was kept open to reduce ground pressure on the main tunnel segments. 7 FROST HEAVE LOADS AND TUNNEL DEFORMATION During freezing activities tunnel convergence and offset of tunnel segments was measured (in the bored tunnel tube). Deformation and convergence of the main tunnel construction has been monitored during the freezing period. At cross passage DV1 (in sand) a maximum offset of tunnel segments of 4mm was observed. The deformation of the western tunnel at DV1 has been caused by the increase of ground pressure due to frost heave of 50kPa measured in the direction parallel to the freezing tubes (SM1 SM2 //; Fig. 11).
Figure 11 Measured (maximum) frost heave pressures at DV1 and DV2 versus permeability (adapted from Kofoed & Doran, 1996).