1 Introduction
Two parameter classes may be used to characterize clay sediment properties. The first one concerns the geotechnical index properties linked to the particles’ mineralogical nature and can be represented, for example, by the Atterberg limits, or by their equivalent void ratios, , where γs and γw are the unit weight of the solids and of the water, respectively. The second relates to the mechanical parameters (e or wnat) that express the geometric arrangement of particles and their evolution. Particles may be bonded due to cementation during the sedimentation and consolidation processes [2,7]. In this case, the clay is non-remoulded and the link between the particles is an important parameter in the phenomenological behaviour [6], in particular during an oedometric loading in a compressibility study [8,9]. When there is no bond between the particles or when the bond has been destroyed (by a high loading for example), the soil is described as remoulded. It is important to study the behaviour of material that has been remoulded and reconstituted in the laboratory in normally consolidated and overconsolidated conditions [5,14] because the ‘intrinsic’ mechanical properties can then be deduced, making it possible to define a fixed reference framework in order to explain the behaviour of natural, non-remoulded clay. Biarez and Favre [3,4], and later Burland [8], introduced the concept of intrinsic compression properties in order to describe, in one normalized plane, the behaviour of one-dimensional, normally consolidated clay that had been remoulded and reconsolidated (at wsat = 1.5 wL). This paper proposes a comparative analysis of these two approaches. It considers two experimental curves, the first one pertaining to a remoulded reconstituted clay (the Kaolinite P300) which has been the subject of several laboratory studies [14], and the other one to a natural clay taken from deepwater sediments (at a depth of 700 m) in the Gulf of Guinea [13].
2 The Biarez & Favre (B&F) model
The soil – a granular and discontinuous medium (DM) – may be seen as a virtual continuous medium (CM), using, for statistical connections, the logical framework shown in Table 1. The ‘mechanical properties of soils’ thus appear as the integration of the grain properties (called ‘nature of grains’) in their spatial configuration (called ‘arrangement of grains’) and their mechanical boundary conditions (which can be represented by the consolidation stress tensor). In the case of remoulded clays, Favre [11,12] shows that the Atterberg limits (‘nature of grain’ parameters) allow us to explain the mechanical properties. In mineral clays, these properties are strongly linked by the following relation:
(1) |
Logical framework for statistical connections
Tableau 1 Cadre logique pour les liaisons statistiques
In this way, the model proposed in Fig. 1a, which, plotted with the Casagrande A line, defines the mineralogical nature of the clay.
On the other hand, a large number of results collected by Favre [11,12] show that the 6.5- and 1000-kPa loadings on the oedometric path correspond, on average, to wL and wP, respectively (or eL and eP, considering ), see relation (2) and Fig. 1b. This leads first to equations (3) or (6) for the compressibility index Cc, then to relation (4) when using (1) in (3).
(2) |
(3) |
(4) |
(5) |
(6) |
By comparison with relation (5), proposed by Skempton [16] for remoulded clays, we underestimate by 10% the compression index Cc for wL = 40% and by 5% for wL = 70%. Nevertheless, the simple model represented by relations (2) and (4) was retained for the oedometric test, and hence, by relation (1), it uses only one Atterberg limit, wL. It is a straight line in the or plane, passing through the two points given by relation (2), with a slope defined by equation (4). A single abacus was then obtained in equation (7), being expressed in kPa and :
(7) |
The ‘arrangement of grains’ therefore appears to be directly linked to the consolidation stress tensor, through the mineralogical properties. This connection corresponds in the plane to a straight line called NCRS (Normally Consolidated Remoulded Simplified, see Fig. 1c). Thus, a normally consolidated remoulded clay will be represented by points located on the NCRS, the overconsolidated remoulded clay by points below the NCRS. In this case, the overconsolidation of the material only depends on the clay loading history, expressed by the overconsolidation ratio . Here is the maximum effective stress and the overburden effective stress.
3 Burland's model
Burland [8] collected and analyzed the oedometric compressibility of several clay sediments remoulded and reconstituted in the laboratory and with liquid limit wL varying from 25% to 159%. The curves obtained in the () plane are slightly concave when is located between 10 and 1000 kPa (Fig. 2a). The intrinsic compressibility index is therefore introduced (equation (8)) as the difference between , two characteristic quantities corresponding to the consolidation stresses at 100 and 1000 kPa, respectively. is also the slope of the oedometric compression curve (linear in this interval of stresses), the quantity {−[log(100)−log(1000)]} being equal to 1. With the transformation of the variable defined in (9), and by introducing the normalized void index Iv, the Burland model requires all experimental curves to pass through two fixed points in the () plane, corresponding to and for and , respectively.
(8) |
(9) |
The ICL (Intrinsic Compression Line) curve, corresponding to equation (10) and represented by Fig. 2b, can then be deduced from a statistical model based on the experimental curves passing through the two specific points (, Iv = 0) and (, Iv = –1). is expressed in kPa.
(10) |
The two particular parameters for the ‘arrangement of grains’ are, on the other hand, connected by equation (11) to eL, the void ratio corresponding to the liquid limit, and the only parameter for the ‘nature of grains’.
(11) |
3.1 Comment
Although both approaches appear to be based on a similar method, there is a fundamental difference in the manner in which the intrinsic characteristics of the ‘arrangement of grains’ were defined on the oedometric path. Burland does this using two experimental values () for the ‘arrangement of grains’. On the other hand, the B&F method, by introducing IL, used two parameters of ‘nature of grains’ (eL and eP). However, both models link the ‘nature of grains’ parameters to the ‘arrangement of grains’ ones, using statistical correlations under oedometric loading. This corresponds to relation (2) in the B&F model and to relation (11) in Burland's model. An examination of Figs. 2a and 3a reveals a difference in the Cc (B&F approach) and (Burland approach) slopes, in particular between and , when the liquid limit is high (wL = 100% to 160%). On the other hand, for the lower wL values, Cc and are similar, and fit well the experimental data given by Burland [8] for different values of wL (see Fig. 3a).
4 Comparison of models and mixed model
4.1 B&F's model in Burland's space ()
According to correlation (7), we have:
(12) |
Using Eqs. (9) and (12), to express B&F's model as a function of Iv, we get the NCRS equation (13) in Burland's space.
(13) |
(14) |
(15) |
(16) |
Equation (16) is a straight line in the () plane, passing through two specific points (Iv = 0, ) and (Iv = −1, ), represented by the NCRS line in Fig. 4a. A representation of both the NCRS and ICL lines on Burland's plane shows that the two curves converge in a vertical stress range between 40 and 1000 kPa. Normalization by means of the Iv parameter and assumption are, of course, fundamental to this result, because they force the NCRS to have the (Iv = −1, ) point on the ICL line. However, in the stress range lower than 40 kPa, the two lines diverge slightly due to the concavity of the curves (Fig. 2a), which is taken into account in the ICL expression. Both approaches may thus be summarized in Table 2, with the indices for ‘B&F in Burland space’ being deduced from both interpolation (15) and relation (1), considering and .
Characteristics of the models
Tableau 2 Caractéristiques des modèles
Burland's model | B&F's model | B&F's mixed model | |
Equation | |||
Characteristic points | |||
Compression index |
4.2 B&F's ‘mixed’ model represented in Burland's space ()
In the linear relation (13) resulting from B&F, are replaced by expressions (11) proposed by Burland for the ICL expression. This gives the ‘mixed’ model (17):
(17) |
Here ζ and χ are the parameters given by equation (18) as function of the eL and eP limits, defined by their experimental values:
(18) |
The ‘mixed’ model described above is applied to the two kinds of clay and are represented in the () plane: the Kaolinite P300 (Fig. 4a, NCRSa(k)) whose Atterberg's limits are wL = 40% and wP = 20%, and the deepwater sediments from the Gulf of Guinea (Fig. 4a, NCRSa(m-s)), whose wL and wP limits vary according to depth (therefore according to too) from 110% to 160% for wL and from 30% to 90% for wP [13]. Since the assumption of is not considered, the NCRSa line is not forced to pass through a point on the ICL line. The NCRSa line, however, will depend on the eL limit. When equation (17) is applied to the kaolinite (through the eL value), we obtain a curve termed NCRSa(k) located very close to the ICL line, whatever the stress domain is. Both correlations are thus very similar, and compatible with the experimental curve, on the oedometric path, of the kaolinite represented in Fig. 4a in the same plane (the kaolinite having first been consolidated under one-dimensional conditions under an axial stress of 140 kPa). On the other hand, the NCRSa(m-s) curve calculated for the marine sediments, through their eL values, diverges from the ICL and the NCRS. The principal difference between the ICL and the NCRS, on the one hand, and NCRSa(m-s), on the other hand (whereas the calculation of the ‘mixed model’ NCRSa was focused on the remoulded state condition) seems to be due mainly to the fact that the ‘mixed model’ does not take into account the hypothesis of Biarez's model, expressed by . Two experimental curves (at z = 9.5 m and z = 16.5 m depths) of the natural non-remoulded marine sediments were plotted in Fig. 4a to compare the two models, in particular at the end of oedometric loading, when the initial structure was destroyed.
The experimental path evolves first from one point located above the ICL curve (sensitive clays) towards a domain that exhibits a ‘fictitious consolidation stress’ and a supplementary cohesion on the oedometric curve. The latter is due to the beginning of the degradation of the cementation between particles developed in situ during the sedimentation and the consolidation processes. Afterwards, the cementation is gradually destroyed (destructuring according to Leroueil et al. [15], and the experimental curves tend to converge at high stresses (up to 1000 kPa) towards the ICL curve of the remoulded state of the material. Clays structuring has been quantified by Cotecchia and Chandler [9,10] using the sensitivity parameter St (Skempton and Northey [17]) as a normalizing parameter for the structure. Thus, the evolution of St may be used to describe the destructuration phenomenon (destructuration law of Baudet and Ho [1]).
4.2.1 Remark
Burland's relation (10) expressed in the (IL – log σ’v) plane using equation (11) becomes equation (19) as a function of the eL and eP limits, ζ and χ being defined by relation (18):
(19) |
Since the ICL depends on eL (Fig. 4b), it is transformed in the () plane by a bundle of parallel curves. These curves, for low wL values, are located in the upper part of the bundle, whereas for high values of wL, they are located in the lower part. In this plane, the similarity between the two models can be clearly seen, which is also in accordance with the experimental curve for the low wL values (case of the kaolinite ICL(k)). As previously in Burland's plane, the ICL(m-s) of marine clay, located in the high wL domain, diverges from the NCRS. The experimental oedometric loading curves (at z = 9.5 m and z = 16.5 m) converge towards the ICL(m-s) in the domain of high stresses, which shows that Burland's correlation is closer to the experimental results for these high values of wL than B&F's one.
B&F's assumption that cannot be introduced in Burland's model, because it would force the ICL(m-s) to converge to the NCRS at (IL = 0, ) point; in this case the curve no longer reflects the experimental observation.
5 Conclusion
Two intrinsic compressibility concepts proposed for the remoulded reconstituted clays are discussed in this paper. The NCRS, resulting from the Biarez and Favre's correlation [2,3], is based on the liquid index parameter IL = (w–wP)/(wL/wP). The ICL, resulting from the Burland's correlation [8], is based on a normalized void index parameter .
The ‘arrangement of grains’ in Biarez and Favre is defined through wL and wP (or eL and eP) – two ‘nature parameters’ – whereas in Burland, this is done by means of the two ‘arrangement parameters’ . However, the two approaches have in common the fact that they propose statistical correlations between the ‘nature of grains’ parameters and the ‘arrangement of grains’ ones for the normalized oedometric tests.
When the NCRS is represented in Burland's plane (Iv,), the difference between the two statistical models is small and the two approaches appear to be similar. The ‘mixed’ model, however, shows that their behaviour is different and that the similarity observed with the NCRS is due to the assumption , which imposes the (IL = 0, ) point as a common point to both curves. If the ICL and the NCRS are, indeed, close when wL is low, they nevertheless diverge for the high values of wL. This result is obtained by representing the two curves in B&F's plane (IL,), which appears to be the more explicative plane. This study shows that Burland's ICL line reflects better the experimental results and represents the compressibility for a large range of wL. On the other hand, the NCRS in the B&F model is less well adapted to the high values of wL, but remains valid for representing compressibility in the case of medium and low values of wL. The NCRS also has the advantage of being simpler to use because it is based on two current ‘nature parameters’: the Atterberg's limits wL and wP.
Acknowledgements
This work was conducted in the framework of the CLAROM Deepwater sediments project – on the site B Stacor core from the Gulf of Guinea. Project partners were: IFP, IFREMER, Fugro-France, Saipem-SA, Stolt Offshore, Technip, and Total.