WO2013124588A1 - Methode de caracterisation du comportement mecanique de ciments - Google Patents

Methode de caracterisation du comportement mecanique de ciments Download PDF

Info

Publication number
WO2013124588A1
WO2013124588A1 PCT/FR2013/050355 FR2013050355W WO2013124588A1 WO 2013124588 A1 WO2013124588 A1 WO 2013124588A1 FR 2013050355 W FR2013050355 W FR 2013050355W WO 2013124588 A1 WO2013124588 A1 WO 2013124588A1
Authority
WO
WIPO (PCT)
Prior art keywords
hydration
cement
initial
function
cementitious system
Prior art date
Application number
PCT/FR2013/050355
Other languages
English (en)
Inventor
Axel-Pierre BOIS
André Garnier
Jean-Benoît LAUDET
Manh-Huyen VU
Siavash GHABEZLOO
Jean SULEM
Original Assignee
Total Sa
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Total Sa filed Critical Total Sa
Priority to DK13710496.4T priority Critical patent/DK2817615T3/en
Priority to EP13710496.4A priority patent/EP2817615B1/fr
Priority to SG11201405034YA priority patent/SG11201405034YA/en
Priority to BR112014020577-9A priority patent/BR112014020577B1/pt
Priority to RU2014138108A priority patent/RU2628032C2/ru
Priority to US14/379,549 priority patent/US9664648B2/en
Priority to CA2864513A priority patent/CA2864513C/fr
Priority to CN201380019523.4A priority patent/CN104335037A/zh
Publication of WO2013124588A1 publication Critical patent/WO2013124588A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/38—Concrete; ceramics; glass; bricks
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/02—Analysing fluids
    • G01N29/024—Analysing fluids by measuring propagation velocity or propagation time of acoustic waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/38—Concrete; ceramics; glass; bricks
    • G01N33/383—Concrete, cement
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/01—Indexing codes associated with the measuring variable
    • G01N2291/011—Velocity or travel time
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/022—Liquids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/023—Solids
    • G01N2291/0232—Glass, ceramics, concrete or stone
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/024—Mixtures
    • G01N2291/02441—Liquids in porous solids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/025—Change of phase or condition
    • G01N2291/0251—Solidification, icing, curing composites, polymerisation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/028—Material parameters
    • G01N2291/02827—Elastic parameters, strength or force

Definitions

  • the field of the invention is that of the characterization of the mechanical behavior of cements, in particular cements used during operations for cementing boreholes, for example in the context of the exploration and production of hydrocarbons.
  • the invention relates to a method for determining mechanical parameters of a cementitious system of given composition, as a function of time, and as a function of the fineness of the cement system, the pressure and / or the temperature.
  • the cementing of a casing in a petroleum well consists in placing a cement sheath in the annular portion between the extrados of the casing and the wall of the hole, the hole being able to be constituted by another casing or by rock.
  • a cement plug can also be created in the well in the context of drilling or completion operations, or at the time of definitive capping thereof. This sheath or cement plug has a key role in the stability and insulation of oil wells.
  • the sheath or cement plug is obtained by pumping a slag or cement slurry made from cement, water and admixtures.
  • This cement slurry is in a liquid state when pumped.
  • the hydration of the cement particles leads the liquid slag to a solid state, characterized by the existence of a skeleton and pores, thus forming a porous medium.
  • the cement sheath is exposed, during the life of the well and throughout its length, to various mechanical and thermal stresses (in situ conditions) coming from operations carried out on the surface and in the well (pressure tests, change of mud, stimulations cold and hot, production reserves ...) or phenomena arising directly in the subsoil (tank compaction, earthquakes ...), until its abandonment, or even beyond.
  • the cement plug is also subjected to various mechanical and thermal stresses. These stresses can damage the constituent material of the sheath or the cement plug, degrade its mechanical and hydraulic properties and, consequently, modify its contribution to the stability and the tightness of the well.
  • the knowledge of the mechanical behavior of the cement under in situ conditions and of the evolution of this behavior over time is essential for the analysis of the operation of the well during its drilling, its exploitation, and its definitive closure. For example, it is important to know the mechanical behavior of the cement to ensure the tightness of the well during these operations, during the storage and sequestration of fluids in underground reservoirs, such as greenhouse gases, such as C0 2 .
  • This knowledge of the mechanical behavior of cement and its evolution over time makes it possible to simulate, by means of analytical-digital models, the life of the sheaths or cement plugs. In situ conditions affect the formation process of hardened cement, and consequently its mechanical behavior.
  • the main mechanical parameters determined from these measurements are rupture parameters, such as compressive strength or tensile strength (uniaxial or triaxial), and deformability parameters, in particular elastic parameters of the material such as the Young's modulus or the Poisson's ratio.
  • a first category of techniques covers static mechanical tests (compressive loading / uniaxial traction, triaxial, odometric test etc.) on samples that are hardened inside a mold in aging benches, in pressure and temperature, then unloaded to be positioned in a measuring device. These tests require the samples to be brought back to atmospheric pressure and to ambient temperature, in order to equip them with measurement sensors, to position them under a press, then to carry out the actual tests, after a possible delivery of the samples under temperature conditions. and pressure similar to those encountered in the well. These loading / unloading cycles can not only damage the samples prior to measurement, but also disturb the measurement of the characteristics of said samples.
  • this first category of techniques does not allow testing of very young cements, because it is necessary to wait until the cement sample can be removed from its mold to be tested.
  • these techniques allow only a one-off measurement over time, and do not offer the possibility of accessing the evolution of mechanical parameters over time, depending on the hydration of the cement.
  • a second category of techniques covers dynamic tests based on a measurement of ultrasonic wave propagation, and does not include back samples at ambient conditions. These techniques, however, are of limited value because of their indirect nature.
  • static parameters which are the parameters used to simulate the behavior of cladding or cement plugs during their lifetime, are accessible only indirectly: they are calculated from dynamic parameters, determined from the measurement of ultrasonic waves, using correlation formulas; these formulas are themselves obtained from static tests, which may be tainted by error, or even not to cover the field of application of the tested materials.
  • static tests are generally performed on cubic samples of cement under certain pressure and temperature conditions that do not necessarily correspond to the conditions of the dynamic tests.
  • Ultrasonic Cernent Analyzer UCA type measurements are part of these techniques.
  • a UCA type measurement involves measuring the velocity of the compression waves and using a correlation to evaluate the compressive strength of the cement.
  • a third category of techniques includes some proposals for static mechanical tests without a loading / unloading step as mentioned above which necessitated bringing the cement sample back to pressure and ambient temperature before making a measurement.
  • EP 1541987 discloses a system in which a cement composition is poured into a bone-shaped mold, where the sample is aged under temperature and under pressure, and where uniaxial tension loading is performed until on breaking the sample, without having to unload the sample.
  • this method does not make it possible to carry out the measurements under in situ conditions, the pressure being able to be exerted only on two faces of the sample, the other faces undergoing a condition of loading by reaction of the mold and not by application of constraints in situ conditions.
  • the measures are therefore biased.
  • only tensile tests are possible, but these are biased as regards the measurement of elastic constants compared to compression tests, due to the appearance of microcracks which invalidate the elasticity hypothesis.
  • the range of determination of the elastic parameters is therefore very small.
  • the geometry used is unconventional.
  • US Pat. No. 7,549,320 corresponds to a technique of the same type, with a variation of the loading technology.
  • Sample taking is carried out, as in US Pat. No. 7,089,816, inside a flexible membrane.
  • the rigid containment enclosure surrounding the flexible membrane is compartmentalized by active sealing devices, and allows the application of a differential pressure to the sample, by the injection of fluids.
  • the flexible membrane used in this device can also be semi-permeable, thus limiting the deformation of the sample during its setting, but disrupting the hydration process of the cement sample.
  • US Pat. No. 7,552,648 describes yet another variant, in which a fluid is injected into the sample itself, which is porous, in order to obtain the desired pressure. Then a tensile test is carried out. No compression test is provided, and the supply of external fluid does not correctly simulate the water exchange under in situ conditions, as is also the case for the technique according to US Pat. No. 7,549,320.
  • this type of test does not make it possible to measure on cements of very young age, that is to say from the very first hours of hydration. nor to follow the evolution of the mechanical parameters measured over time, according to the hydration of the cement. None of the three categories of techniques described makes it possible to measure, without artefacts, the mechanical properties of a cement during setting.
  • the invention aims to respond to a need for simple characterization of the mechanical behavior of a cementitious system, through the determination of static mechanical properties and hydro-mechanical properties expressed as a function of time, and as a function of different corresponding pressures and temperatures. conditions encountered along a wellbore.
  • Another object of the invention is to provide an easily industrializable method for characterizing the mechanical behavior of a cementitious system, which uses simple and standardizable measurements.
  • Another object of the invention is to provide a method for estimating the evolution of the mechanical parameters of a cement system over time, from the beginning to the end of its hydration.
  • the invention also aims to provide an estimate of static mechanical parameters of the cement, which are the parameters used to reliably simulate the mechanical behavior of the sheaths or cement plugs during their lifetime.
  • FIG. 1 represents the logic diagram of the method according to the invention.
  • Figure 2 is a diagram illustrating the different stages of the hydration process of a Portland cement.
  • FIG. 3 shows an example of an Ultrasonic Cernent Analyzer (UCA) type device making it possible to measure the speed of compressional waves in a cement sample used in the process according to the invention.
  • UCA Ultrasonic Cernent Analyzer
  • FIG. 4 is a diagram showing the velocity of compressional waves as a function of time V p (t) measured during a UCA-type test on a cementitious system comprising Portland cement, as well as the degree of hydration as a function of time a (t) of the same cement system measured during a calorimetric test.
  • FIG. 5 is a diagram illustrating the linear relationship between the velocity of compressional waves as a function of time V p (t) and the degree of hydration as a function of time a (t) as represented in FIG.
  • FIGS. 6 to 8 are diagrams comparing the degree of hydration as a function of time a (t) determined by the process according to the invention for a given cementitious system, and the degree of hydration as a function of time a (t) for this same system, determined experimentally.
  • the degrees of hydration are represented for different values of pressure and temperature.
  • a (t) is shown for a low pressure of 0.3 MPa and at different temperatures (7 ° C, 13 ° C, 25 ° C, 40 ° C and 60 ° C).
  • a (t) is shown at higher pressures, depending on the following pressure / temperature torques: (40 MPa, 30 ° C.), (40 MPa, 15 ° C.), (20 MPa, 30 ° C.).
  • a (t) is represented for the following two pressure / temperature pairs: (0.3 MPa, 60 ° C.), (20 MPa, 60 ° C.),
  • FIG. 9 is a diagram illustrating the composition of a given cementitious system as a function of the degree of hydration of the system, determined by the process according to the invention.
  • the volume fraction of the various constituents of the cementitious system as a function of the degree of hydration of the cementitious system is given, for a fixed pressure and temperature torque.
  • Figures 10 to 13 are diagrams comparing the mechanical parameters as a function of time determined by the method according to the invention, for a given cementitious system, those experimentally determined and from the literature (Boumiz et al 1 - *, and Haecher al 2 - *) for this same cement system.
  • the present invention provides a method for determining mechanical parameters of a cementitious system of initial composition Co and fineness ⁇ , as a function of time, and as a function of the fineness of the cement system, pressure and / or temperature. This process comprises the following steps:
  • step B determining the degree of hydration a (t) as a function of desired values of fineness ⁇ ⁇ of the cementitious system, of pressure P n and / or of temperature T n ;
  • C determining the composition of the cementitious system as a function of time C (t) and as a function of desired values of fineness ⁇ ⁇ of the cementitious system, of pressure P n and / or of temperature T n , from the degree of hydration a (t) determined in step B;
  • the method makes it possible to characterize the mechanical behavior of a cementitious system used as a sheath or cement plug in a well, preferably a petroleum well.
  • the method according to the invention is not limited to this use, and can be implemented in the context of any well construction comprising a sheath or a cement plug for which the characterization of the mechanical behavior of the cementitious system used is desired, for example in the field of gas storage or geothermal energy.
  • a cementitious system within the meaning of the present invention corresponds to a cement slurry or a hardened cement.
  • cement slurry is meant a fluid composition (liquid, pasty, granular, etc.) based on cement and water, capable of hardening, that is to say likely to change to a solid or almost solid state. over time (undergoing a hardening step).
  • cement grout is used as a synonym for cement slag.
  • the hardening essentially corresponds to the hydration (or setting) of the cementitious and water-based composition.
  • Hardened cement is understood to mean a cement slag which has been hardened, and in which a porous skeleton has begun to develop, which can give it mechanical strength.
  • a hardened cement is not necessarily a completely hardened cement, that is to say whose hydration is complete.
  • a fully hardened cement is a cement whose hydration is 100% complete.
  • static mechanical parameter is meant a parameter which is not obtained from the measurement of wave speeds.
  • Static parameters are the parameters used to reliably simulate the mechanical behavior of cladding or cement plugs during their lifetime, unlike dynamic parameters. These latter, generally derived from sonic / ultrasonic wave measurement methods, differ from the static parameters due to the dependence of the mechanical properties of the cement on the frequency of the waves used and the size of the cement samples tested.
  • fineness of a cement is meant the grinding fineness of a cement expressed by its specific surface (developed surface per unit mass, expressed in cm 2 / g), and measured using a device of the type Blaine according to standard NF EN 196-6.
  • compressional wave or "P-wave” is understood to mean sonic or ultrasonic waves created by the mechanical vibration of a fluid or solid support and propagated by virtue of the elasticity of the surrounding medium in the form of longitudinal waves.
  • reactive initial phase is meant an initial phase of the cement of the cement system, able to hydrate to give a hydrated phase, also called hydrate in the present description.
  • the process according to the present invention applies to all the cements whose hydration is through a series of exothermic chemical reactions between its components and water, for example Portland cements, aluminous cements, and other cements defined in accordance with EN-196-2, EN-196-6, EN196-7, ISO 3310-1, ISO 13500, and ISO 10426-1 to ISO 10426-6 for the petroleum industry .
  • the cementitious system according to the present invention may comprise adjuvants, making it possible to obtain the desired properties:
  • Accelerators for example: calcium chloride
  • the objectives of which are to reduce the induction period of the cementitious system, and possibly to accelerate the development of uniaxial compressive strength
  • ⁇ retarders eg lignosulphonates, which aim to increase the induction period of the cement system
  • Extenders eg bentonite, hollow beads, which aim to reduce the density of cement slag and / or reduce the concentration of cement or increase the concentration of water;
  • ⁇ weights for example: hematite which aim to increase the density of cement slag
  • ⁇ dispersants for example lignin sulphonates
  • lignin sulphonates which are intended to reduce the viscosity of cement slurry in the cement grains deflocculant
  • ⁇ filtrate controllers (example: latex), which aim to limit the transfer of water cement slurry to the external system by building a cake little or not permeable against the formation;
  • ⁇ loss controllers eg gilsonite, which aim to limit cement losses in highly permeable areas
  • This initial composition Co comprises the initial phases of the cement, among which at least one reactive initial phase which can be hydrated, and the initial water, according to given fractions, for example expressed in volume.
  • the hydration of the cementitious system is advanced, and the composition of the cement C (t) is different from the initial composition C.
  • the composition of the cement C (t) comprises m phases, among which the initial phases of the cement in different proportions of Co, at least one hydrated phase resulting from the hydration of at least one reactive initial phase of the cement, and water.
  • the cement of the cementitious system comprises at least one initial reactive phase X, preferably chosen from the group consisting of tricalcium silicate of formula 3CaO.SiO 2 (Alite), of symbol C 3 S dicalcium silicate of formula 2CaO.SiO 2 (Belite), of C 2 S symbol, tricalcium aluminate of formula 3CaO.Al 2 O 3 (Aluminate), of C 3 A symbol, alumino ferrite tetracalcic C 4 F, and their associations.
  • X preferably chosen from the group consisting of tricalcium silicate of formula 3CaO.SiO 2 (Alite), of symbol C 3 S dicalcium silicate of formula 2CaO.SiO 2 (Belite), of C 2 S symbol, tricalcium aluminate of formula 3CaO.Al 2 O 3 (Aluminate), of C 3 A symbol, alumino ferrite tetracalcic C 4 F, and their associations.
  • the cement of the cementitious system comprises at least one hydrated phase resulting from the hydration of at least one reactive initial phase of the cement, preferably chosen from the group consisting of hydrated calcium silicate CS-H, calcium hydroxide CH, calcium trisulfoaluminate hydrate TSA (Ettringite), hydrated calcium monosulfate AFm, alumino calcium ferrite hydrate, and their combinations.
  • the initial phases such as silicates C 3 S, C 2 S and aluminates are reactive phases that react with the water of the cementitious system to give hydrated phases, by means of exothermic chemical reactions.
  • the silicates react with water to form hydrated calcium silicates CSH, and calcium hydroxide CH (Portlandite of formula Ca (OH) 2 ).
  • the aluminates react with water and gypsum, or with water alone to form the following hydrated aluminates: TSA hydrated calcium trisulfoaluminate (Ettringite), hydrated calcium monosulfate AFm and hydrated calcium alumino ferrite.
  • the microstructure of the cement system, especially hardened cement is very complex. The volume fractions, the possibility of crystallization, morphology and phase sizes of hydrated hardened
  • Porlandite CH is formed of broad hexagonal crystals which are surrounded by the CSH phase, which is in the form of gels, and occupies about 20% of the volume of the hardened cement. Monosulfoaluminates and PEttringite play a minor role in the structure of hardened cement and generally represent a volume fraction of 15% to 20%.
  • the CSH phase is the major product of hardened cement.
  • the CSH phase is a porous phase, of amorphous and colloidal structure, and has a varied chemical composition. It occupies between 50% and 70% of the total volume of hardened cement, which gives the properties of this phase a dominant role in the macroscopic properties of the cementitious system.
  • the CSH phase exists in two forms: the "high density” phase CSH HD (also referred to as “HD” in the present description) and the “low density” phase CSH LD (also referred to as “LD” in the present description).
  • CSH HD also referred to as “HD” in the present description
  • CSH LD also referred to as "LD” in the present description
  • Each of these two forms has, according to Jennings et al. (8 ' 9) , a common unit called globule.
  • the cementitious system comprises Portland cement, which is one of the most common cements, including the following initial phases used (see standards EN-196-2, EN-196-6 , EN196-7, ISO 3310-1, ISO 13500, and ISO 10426-1 to ISO 10426-6 for initial phase contents, depending on the grade and grade of the cement):
  • the silicate phases and aluminates are derived from a milled clinker, the main constituent of Portland cement.
  • the process according to the invention is not limited to the cement system comprising a Portland cement.
  • the method can be adapted for the determination of mechanical parameters of various cementitious systems, having chemical compositions and hydration reactions different from those of a cementitious system comprising a Portland cement, without departing from the scope of the present invention. 'invention.
  • the hydration process of the cementitious system takes into account different phases of the hydration process of a cementitious system for the determination of at least one mechanical parameter thereof.
  • the hydration process can be followed by the evolution of the heat flux created by the exothermic hydration reactions of the initial reactive phases of cement with water.
  • the amount of heat released by the exothermic reactions during the hydration of the cementitious system can be measured using a calorimeter.
  • Two types of calorimeter are generally used. The first type corresponds to quasi-adiabatic calorimeters, in which the reaction takes place inside a perfectly thermally insulated enclosure. This type of calorimeter is not suitable for measuring samples of cement systems under high pressures, and the operating pressure limit is about 0.1 MPa.
  • a second type corresponds to the isothermal calorimeters, with which the tests are carried out in a space of the order of cubic centimeter, and where the temperature is controlled so that the temperature of the space does not vary. The materials used for this type of measuring device make it possible to work with higher pressures, of the order of Mega Pascal.
  • the hydration of the cement system comprises in this order:
  • a dormant period also called an induction period, during which an increase in the concentration of C 3 A and of OH " ions occurs;
  • FIG. 2 illustrates the hydration process in the case of Portland cement (Rixom et al (15) ).
  • the hydration process is usually divided into 5 steps. Each step corresponds to a trend of evolution of the heat caused by the hydration reactions:
  • Period 1 an initial hydration period proceeds for about 15 minutes, during which the more soluble chemical components (Na 2 S0 4, K 2 S0 4, CaS0 4 .xH 2 0) of the cement dissolves, by increasing the pH of the cement slurry.
  • this third period corresponds to an acceleration of the hydration, which lasts approximately 4 hours. During this one, the heat of hydration evolves very quickly.
  • the hydrates formed come into contact, which forms the strength of the dough.
  • Period 4 a period of slowdown then takes place and lasts a few days. The evolution of the heat of hydration slows down.
  • Period 5 a cure period occurs for several months after the beginning of the setting of the cement system, and the evolution of the amount of hydrates is gradually stabilized.
  • the kinetics of hydration vary during the hydration process. Among the periods of the hydration process, a first stage where the hydration is governed mainly by a nucleation and growth process, as described for example by Avrami (3) , is conventionally recognized.
  • a second stage of the hydration process of the cementitious system is characterized by a hydration governed mainly by an ion diffusion process.
  • This second stage begins when the degree of hydration has reached a value of a degree of hydration threshold a *.
  • a degree of hydration threshold a * In the first stage, governed by the nucleation and growth phenomena, a layer of hydrates is created around the grains of cement, and is relatively permeable, which allows the diffusion of the ions and the progress of the chemical reactions.
  • a threshold degree of hydration is reached, a decrease in ion diffusion occurs, related to the thickening of the hydrate layer around the cement grain and the lower permeability of this layer.
  • CSH hydrated phases and CH are the main constituents of hardened cement, and are primarily responsible for the strength of hardened cement.
  • the hydration kinetics of the C 3 S phase is faster than that of the C 2 S phase, and the amount of CH created by the hydration of C 3 S is approximately three times greater than that created by hydration. C 2 S.
  • phase C 3 S plays a major role in the process of hydration and in the development of the mechanical behavior of the cementitious system, introducing the cement paste resistance during the first days.
  • the method takes into account the two reactions (I) described above to simulate the hydration process.
  • the aim of the invention is to characterize the mechanical behavior of a cementitious system, and in particular to determine certain mechanical parameters and their evolution over time, under in situ conditions, which can then be used, for example, in modeling the behavior of sheaths or cement plugs from boreholes during their lifetime.
  • the method according to the invention makes it possible to determine at least one mechanical parameter of the cementitious system as a function of time, for desired values of fineness ⁇ ⁇ of the cementitious system, of pressure P n and / or of temperature T n .
  • the method makes it possible to determine at least one mechanical parameter which is chosen from among the static deformability parameters, preferably the static elastic parameters such as static Young's modulus E, the static Poisson's ratio v , the incompressibility module K, the shear modulus G, and their combinations, such as for example the odometric module K v , or the lamé ⁇ module.
  • the static elastic parameters such as static Young's modulus E, the static Poisson's ratio v , the incompressibility module K, the shear modulus G, and their combinations, such as for example the odometric module K v , or the lamé ⁇ module.
  • two static elastic parameters are determined.
  • At least one mechanical parameter determined according to step D is a parameter under drained conditions.
  • static deformability parameters preferably static elastic parameters
  • a flow of pore fluids for example liquids
  • an undrained system it is assumed no exchange of pore fluids between the cement system and its environment.
  • the value of the mechanical parameters of the cement system, in particular the static deformability parameters will be different depending on whether the system is drained or not drained.
  • the method furthermore comprises the determination of at least one hydro-mechanical coupling parameter, such as the Biot b coefficient or the Skempton B coefficient, and their combinations.
  • hydro-mechanical coupling parameters such as the Biot b coefficient or the Skempton B coefficient
  • the knowledge of such hydro-mechanical coupling parameters makes it possible in particular to determine the static deformability parameters, preferably the static elastic parameters, under undrained conditions from the knowledge of the static elastic parameters, under drained conditions.
  • These parameters can, for example, be used to simulate the behavior of a cement cladding whose pores are saturated with water when the loading speed is significantly faster than the rate of diffusion of the pore pressures.
  • the set of formulas (II) used to pass from static elastic properties drained (index d) to static undrained properties (index u), and vice versa, is given below.
  • the shear modulus G is generally not measured directly, although it can be measured theoretically, and is evaluated from the knowledge of two other elastic parameters.
  • the method according to the invention makes it possible to determine the static deformability parameters, which have the advantage of being those used in the modeling of shafts or plugs of wellbore cement. These static parameters should not be confused with those called “dynamic”, which are evaluated from the knowledge of the speed of the compression and shear waves.
  • the hydro-mechanical coupling parameters are for example the Biot coefficient
  • ⁇ is the isotropic stress variation
  • ⁇ ⁇ is the pore pressure variation
  • these parameters can be determined for the cement system from an early age, which allows a better modeling of the mechanical behavior of the cement system. It may be important to know the behavior of the early stage cement in the context of drilling operations, for example for operations that must be carried out quickly after pumping cement slurry into the well, such as thermal logs allowing locate the top of the cement column, additional drilling operations once the casing is laid, the perforation of the production zones, or to evaluate the state of stress in the cement sheath once the cement has hardened.
  • FIG. 1 schematically represents the method according to the invention, which comprises steps A to D described below.
  • the cementitious system tested by the process according to the invention has a known initial composition Co and a fineness ⁇ .
  • Another initial parameter used by the process according to the invention is the speed of the compressional waves in the cementitious system as a function of time V p (t). The knowledge of these initial parameters makes it possible, thanks to the method according to the invention, to determine at least one mechanical parameter of the cementitious system.
  • the method further comprises an initial step of measuring the compression wave velocity as a function of time V p (t) in a sample of the cement system.
  • a measurement can, for example, be carried out according to a method of Ultrasonic Cernent Analyzer UCA type, such as for example described in US Pat. No. 4,259,868.
  • UCA-type methods are well known methods for evaluating the quality of a cementitious system, and are part of the second category of measurement methods described above, making it possible to characterize the mechanical behavior of a cementitious system.
  • FIG. 3 illustrates a UCA type measuring device, as described in the patent
  • Such a measuring device comprises a measuring cell 34 intended to receive the cement slurry and to be pressurized.
  • This cell 34 is pressurized in an autoclave 33 through its insertion into an opening 34a and a connection with a pressurizing conduit 34b.
  • the pressure and the temperature of the autoclave, transmitted to the cell 34, are controlled by inverters (37a, 36) and control screens (37, 35).
  • Ultrasonic transducers (not shown) are coupled to the measuring cell 34, and allow the emission and reception of ultrasonic waves passing through the cementitious system sample.
  • a control computer 31 is connected to the autoclave and comprises a keyboard 38 for inputting data and control functions into the measuring device, and a digital plotter 32 for displaying the measured wave velocity and the resistance uniaxial compression calculated in real time.
  • Step A consists in determining the degree of hydration of the cementitious system as a function of time a (t), from the knowledge of the compression wave velocity V p (t) of the cementitious system in question, measured in a sample of the cement system at a pressure Pi and a temperature ⁇ .
  • the method according to the invention uses an empirical correlation 1 10 which connects a (t) to V p (t).
  • the degree of hydration of the cementitious system as a function of time a (t) is calculated from V p (t) in a linear relation.
  • the degree of hydration of the cementitious system as a function of time a (t) is calculated from V p (t) according to the following linear relationship:
  • V p and Vo are known at the beginning of step A.
  • the value of V ⁇ can be known experimentally with the same measurement as that allowing to know V p (t), and corresponds to a sample in which the hydration is complete. Knowing V p (t), V ⁇ is for example calculated using a linear regression.
  • the value of V ⁇ is derived from a model established from preliminary experimental data from UCA type tests, making it possible to predict this value for any cementitious system. According to this embodiment, it is not necessary to carry out a measurement of V ⁇ for the cementitious system that is the subject of the process according to the invention.
  • V ⁇ is equal to approximately 3980 m / s.
  • the method according to the invention comprises a step B consisting in determining the degree of hydration a (t), resulting from stage A, as a function of desired values of fineness ⁇ ⁇ of the cementitious system, of pressure P n and / or temperature T n .
  • this determination uses a kinetic model 120 of the hydration of the cement system.
  • this kinetic model 120 comprises two stages corresponding to two substeps of step B: B-1 and B-ii.
  • the hydration process comprises a first stage where hydration is governed mainly by a nucleation and growth process, and a second stage where hydration is governed mainly by a diffusion process ions, said second stage beginning when the degree of hydration has reached a threshold degree of hydration value a *, this threshold value a * being a function of the temperature.
  • Step B comprises the following substeps:
  • Each of the substeps B-i and B-ii takes into account the fineness ⁇ of the cement system, the pressure and the temperature for the determination of the degree of hydration a (t).
  • the method therefore advantageously takes into account the effects of temperature and pressure to simulate the evolution of the degree of hydration of the cement system.
  • the value of degree of hydration threshold a * is evaluated by minimizing the difference between a (t), determined using the kinetic model 120, and a (t), determined experimentally to from the velocity of the compression waves, for different temperatures, and at a constant pressure, so as to take into account the variation of a * as a function of the temperature in step B.
  • the experimental determination of a (t ) is carried out at different temperatures, and at a constant pressure less than or equal to 1 MPa, preferably less than or equal to 0.5 MPa, more preferably less than or equal to 0.3 MPa. At such pressures, it can be considered that the effect of the activation volume is negligible.
  • the cementitious system of initial composition Co comprises a cement and water, the cement comprising at least one reactive initial phase X;
  • the degree of hydration a (t) determined in step B corresponds to the weighted average of the degrees of hydration of each of the initial reactive X phases of the cement;
  • the degree of hydration of each of the initial reactive phases X of the cement is a function of the ratio between the chemical affinity A x (a) of the initial reactive phase X, this chemical affinity A x (a) controlling the rate of evolution hydration of the initial reactive phase X, and the characteristic time associated with the reaction of the initial reactive phase X with the water ⁇ x ;
  • the characteristic time associated with the reaction of the initial reactive phase X with the water ⁇ ⁇ is a function of the fineness ⁇ of the cement, the pressure and the temperature.
  • This substep Bi it is considered the first stage of evolution of the degree of hydration or hydration is governed primarily by a process of nucleation and growth, such as for example described by Avrami 1 ⁇ .
  • This first stage comprises the phases of acceleration and slowing down of the hydration process in the case of a Portland cement.
  • This substep Bi consists of determining the degree of hydration of the cementitious system from the degrees of hydration of each of the initial reactive phases X of the cement according to the following relationship:
  • N is the number of initial reactive phases X.
  • X would be one of the four main constituents of clinker (C 3 S, C 2 S, C 3 A, C 4 AF).
  • m x ' is the mass fraction of component X.
  • the degree of hydration of the initial reactive phase a x is defined as the ratio between the amount of this reacted phase and the initial quantity of this phase, and verifies the following equation:
  • a x (a) (1 - (a - a 0x )) ⁇ [- ln (l - (a - a 0x )
  • n x and cio x are experimentally determined constants and known to those skilled in the art (Bernard et al 12 - *).
  • the characteristic time associated with the reaction of the initial reactive phase X with the water ⁇ ⁇ is a function of the fineness ⁇ of the cement, the pressure and the temperature.
  • ⁇ ⁇ verifies the following equation:
  • R the constant of perfect gases, equal to 8,314 J ° K-lmol-1
  • ⁇ ⁇ activation energy AV X activation volume for nucleation and hydrate growth during the growth and deceleration phases of the hydration process.
  • hydration equal to -27 ⁇ 10 -6 m 3 / mol
  • n x is an experimentally determined constant known to those skilled in the art (Bernard et al (12) ).
  • sub-step B-ii we consider the second stage of evolution of the degree of hydration, where hydration is governed mainly by an ion diffusion process.
  • This second stage includes the cure phase in the case of Portland cement.
  • the value of the radius of the cement grain surrounded by a layer of hydrates is called r * when the degree of hydration of the cement reaches a threshold value a *.
  • the value of the radius decreases with the progress of the hydration, and its rate of decrease can be written according to a diffusion constant D (cm 2 / h) according to the following equation: dr _ D
  • a (t) can be expressed as a function of the chemical affinity A (a).
  • a (a) satisfy the following equation: da _ A ⁇ a)
  • the kinetic model 120 used in the process according to the invention advantageously takes into account the effect of the fineness of the cement system, the temperature and the pressure to simulate the results of the evolution of the degree of hydration of the cement system.
  • the activation volume AV X is advantageously considered to be identical for each of the initial reactive phases.
  • the value of degree of hydration threshold a * is evaluated by minimizing the difference between a (t), determined according to step A, and a (t), determined experimentally for the cement system, for different temperatures, and a constant pressure.
  • the pressure is less than or equal to 1 MPa, preferably less than or equal to 0.5 MPa, even more preferably less than or equal to 0.3 MPa. At this low pressure, it is indeed possible to neglect the effect of the activation volume.
  • the pressure does not affect the degree of hydration threshold a *, which varies only as a function of the temperature.
  • a * which varies only as a function of the temperature.
  • the diffusion phase starts earlier at a higher hydration temperature, and the higher temperature causes the increase in the long-term CSH HD ratio. term.
  • This step consists in determining the composition of the cementitious system as a function of time C (t) and as a function of desired values of fineness ⁇ ⁇ of the cementitious system, of pressure P n , and / or of temperature T n , from the degree of hydration a (t) determined in step B.
  • the method according to the invention uses a compositional model 130, which allows to know, at any degree of hydration of the cement system, the content of the system in each phase considered.
  • the volume fraction, by mass or molar, preferably by volume, of each phase of the cementitious system is determined.
  • the knowledge of a (t), for desired values of fineness of the cement system, of pressure and / or temperature makes it possible to estimate the fraction, preferably volume, of the initial phases of the cementitious system (initial reactive phases cement and initial water consumed), as a function of time.
  • the fraction for example volume fraction, of the hydrated phases of the cementitious system is calculated at a given time. Indeed, the chemical equations make it possible to determine the proportion in numbers of moles of the reactants and the hydrates.
  • the volume fraction is then determined by knowing the density and the molar mass of the hydrates.
  • the fractions, preferably volume fractions, of the hydrated phases of the cementitious system, in particular the C-S-H and CH phases are determined.
  • the fraction, preferably volume, aluminates is determined to be the difference between the initial proportion, preferably the initial total volume, and the calculated proportion, preferably the calculated volume, of all other phases.
  • compositional model 130 takes into account the temperature at which the hydration is carried out, in particular as regards the determination of the fractions, preferably volume, of the hydrated phases, in particular of the C-S-H phase.
  • composition of the cementitious system C (t) is determined as a function of time, for desired values of fineness ⁇ ⁇ of the cementitious system, pressure P n and / or temperature T n .
  • the cementitious system comprises m phases, among which:
  • At least one reactive initial phase X preferably chosen from the group consisting of C 3 S tricalcium silicate, C 2 S dicalcium silicate, C 3 A tricalcium aluminate, C 4 F alumino ferrite tetracalcium, and combinations thereof;
  • At least one hydrated phase Y resulting from the hydration of at least one reactive initial phase X preferably chosen from the group consisting of calcium silicate hydrate CSH, calcium hydroxide CH, calcium trisulfoaluminate hydrate TSA hydrated calcium monosulfate AFm, alumino calcium ferrite hydrate, and their combinations; - some water;
  • the initial composition of the cementitious system Co comprising a determined initial volume of water V ° and at least one initial reactive phase X.
  • Stage C comprises estimating the composition of the cementitious system as a function of time C (t) by determining the mole fraction, by mass or volume, of the m phases of the cementitious system.
  • V x (t) V x ° (l - a x (t)) (XIV) with V x the initial volume of the initial reactive phase X, and x (i) the degree of hydration of the reactive initial phase X as a function of time;
  • V w (t) the volume of water as a function of time
  • V w (t) the volume of water as a function of time
  • V w (t) the volume of water as a function of time
  • the method according to the invention can also take into account other phases, such as incoming water, diffusing environmental medium to the pores of the cement system, or at least one gaseous phase.
  • the method also makes it possible to determine a fraction, for example a volume fraction, of incoming water in the cementitious system as a function of time, and a fraction, for example a volume fraction, of gas as a function of time.
  • the fractions of the incoming water and / or gaseous phases are determined based on the phenomenon of chemical shrinkage occurring during the hydration of the cement system.
  • the volume of the hydrated phase CSH is determined as a function of the temperature at which hydration occurs.
  • the fraction of the CSH LD and CSH HD phases is determined differently depending on whether the hydration is carried out at room temperature (20 ° C ⁇ 5 ° C), or at a temperature above or below room temperature.
  • the fractions, preferably volume fractions, of CSH LD and CSH HD can be expressed according to the following equations (Bernard et al., 12 - *):
  • V HD ⁇ t) ⁇ V _ s _ H (a '(t) - a *) x K [a' (t) - a * ⁇ (XVIII) with H the function of Heaviside, and i denotes each of the reactions having the hydrated phase as the product of the reaction.
  • PLD (1 " ⁇ ) - Pg + ⁇ t> LD ⁇ Pw (XX) with p g the density of the globule, preferably equal to 2.65 g / cm 3 , p w the density of the water, equal to 1 g / cm 3 ), and and LD the respective porosities of CSH HD and CSH LD.
  • the determination of the composition C (t) takes into account the influence of the temperature, via the degree of hydration threshold a * depending on the temperature , on the fractions of CSH LD and CSH HD, and possibly on the fraction of the incoming water.
  • the fractions, preferably volume fractions, of CSH LD and CSH HD can be expressed by the same equations (XVII) and (XVIII) above, which take into account a value of the porosity of the CSH HD phase different from that at ambient temperature.
  • the determination of the porosity of the HD HSC at a given temperature is obtained by assuming that 1) the amount of globules (CSH LD and CSH HD) depends only on the degree of hydration, 2) the volume of the globules is constant for one given hydration temperature, 3) the density of CSH globules is virtually unaffected by temperature, 4) the total porosity is not affected by temperature, and 5) the total pore volume measured by changing the temperature.
  • the oven at 105 ° C is considered to be the sum of pore volumes from globules, CSH LD, CSH HD and capillary volume.
  • the volume fraction of CSH HD and its density are greater (porosity of CSH HD lower) and the volume fraction of CSH LD is lower than at room temperature.
  • the macroporosity of the cementitious system also increases with a higher temperature, for the same degree of hydration.
  • the difference in capillary porosity is, for example, taken equal to 6% between the case at a temperature of 60 ° C. and that at a temperature of 20 ° C., based on mercury porosimeter measurements.
  • the porosities of C-S-H HD and C-S-H LD which are, for example, respectively 24% and 37% at an ambient temperature of 20 ° C., are respectively 17% and 37% at a temperature of 60 ° C.
  • This step consists in determining at least one mechanical parameter of the cementitious system as a function of time, and as a function of the desired values of fineness ⁇ ⁇ of the cementitious system, of pressure P n and / or of temperature T n , from the composition of the C cement system (t) determined in step C.
  • the method according to the invention uses a multi-composition analysis method 140, based on homogenization techniques as used in micromechanics.
  • This multi-composition analysis method 140 takes into account a multi-scale model of the cementitious system for determining, using homogenization techniques, at least one mechanical parameter of the cementitious system as a function of time, for a given value of fineness. of the cementitious system, of pressure and / or of temperature, knowing the composition of the cementitious system C (t) determined in stage C, and knowing the evolution of the cementitious components of the cementitious system obtained from the model of the kinetics of hydration used in step A.
  • the multi-scale model includes at least one elementary scale representing CSH solid phase globules and a macroscopic scale of the cementitious system, preferably a multi-scale model comprising three scales.
  • the physico-chemical heterogeneity of the cementitious system, in particular hardened cement is represented by a multi-scale model.
  • a local scale includes elementary units called globules that represent the CSH solid hydrate phase. This local scale represents the smallest scale of the multi-scale model used in this step.
  • a macroscopic scale represents the macroscopic scale cement system. Intermediate scales can be used.
  • the incompressibility modulus k s and the shear modulus g s of the globules constitute the two unknowns of the multiscale model.
  • micromechanics in material mechanics, consists in estimating the mechanical properties of a material by considering a problem at the limits of a representative elementary volume (VER). according to there are three stages for homogenization: description, localization and homogenization.
  • the VER In the first step, the VER must be well chosen to represent the macroscopic behavior of the material. Thus, to incorporate all the mechanical and geometrical information of the heterogeneous phases of the material, the VER is preferably sufficiently large. Moreover, in the context of continuous media, this volume is preferably small enough to describe the continuity of the macroscopic structure.
  • the interaction of the phases in the VER is considered.
  • Various well-known homogenization schemes are used.
  • the "diluted” scheme does not consider any interaction between inclusions in the VER, and can only be used for low volume fractions of the inclusions.
  • the "Mori-Tanaka” scheme which takes into account the interaction between inclusions, is commonly used for a microstructure of matrix composites and inclusions.
  • the difference between the "Mori-Tanaka” scheme and the "self-consistent” scheme is that the latter takes the homogenized material as the matrix by solving a system of nonlinear equations for the estimation of the actual mechanical properties of the material.
  • the elastic and hydromechanical coupling parameters can be expressed as a function of the strain localization tensors.
  • the tensor of the rigidity of the material is expressed by:
  • G hom is the shear modulus.
  • Zaoui (14) for an isotropic case with n spherical inclusions, the actual modulus of drained incompressibility and the shear modulus are calculated by the following expressions (XXIII) and (XXIV):
  • the material considered comprises a single porous phase (pores) and n-1 solid phase.
  • pores porous phase
  • n-1 solid phase the Biot coefficient tensor of a material having a single porous phase
  • the sign (A) denotes the average volume value of A in volume V.
  • the homogenized Biot modulus is written in the following form (XXVIII): f r c; 1 : (l - 1: ( ⁇ )) (XXVIII)
  • the material is a multi-scale porous material, for example two scales for the sake of simplification.
  • the generalization to N successive scales is then done without difficulty.
  • the following configuration is taken into account: the material comprises q porous phases (solid and small pores), pore volume (large pores) and nq-1 solid phases.
  • the homogenization consists in carrying out two steps I and II: the first step II consists in homogenizing separately q porous phases in order to obtain the effective parameters of each phase of step I (c r , which are respectively the stiffness tensor, the Biot coefficient tensor, the Biot modulus, and the solid phase stiffness tensor of the r phase in step I) as presented in the first case above.
  • the second step II consists of homogenizing q porous phases, a volume of large pores and nq-1 solid phases.
  • the homogenized Biot tensor in the second step is expressed by the following equation (XXX) (Ulm et al (19) ):
  • the heterogeneity of the cementitious system in step D manifests itself at three scales (Jennings ⁇ 11 - *):
  • a first scale Ech.l corresponds to the CSH LD and CSH HD phases, in the form of globules, having a characteristic length of between approximately 10 ⁇ 9 - 10 ⁇ 8 m, and further comprises the CSH gel pores.
  • a second Ech.2 scale corresponds to the cementitious system having a characteristic length greater than about 10 "8 m, and comprises the steps CSH LD, HD CSH, CH, the volume of capillary pores, aluminates and reactive initial phases (the clinker grains).
  • the matrix formed by solid globules, and on the other hand the pores in the CSH LD and CSH HD phases.
  • the matrix occupies for example 63% of volume for the case of the CSH phase LD, and it occupies for example a proportion greater than 76%, depending on the temperature, volume for the CSH HD phase.
  • the Mori-Tanaka scheme (Ghabezloo ⁇ 13 ⁇ ) is chosen, with the solid phase playing the role of the reference medium. This Mori-Tanaka scheme is suitable for materials whose solid phase is dominant and for which there are interactions between the particles (Bernard et al . 12 ).
  • G (I -) 3 ⁇ 4 pooan)
  • X denotes the CSH phase LD or the CSH HD phase
  • k s and g s are respectively the incompressibility modulus and the shear modulus of the globules which represent the solid hydrated phase CSH at l Ech.O elementary scale ("s" as solid).
  • the "self-consistent" scheme is used on this scale. This one allows to take into account the threshold of percolation, and is appropriate to estimate the elastic properties of the cement system at this scale (Ghabezloo 1 '13 ⁇ ).
  • This second Ech.2 scale comprises six phases: CSH LD (LD), CSH HD (HD), CH, the volume of capillary pores (cp), aluminates (AL), and the initial reactive phases of cement, which are example clinker grains (CK).
  • CSH LD LD
  • CSH HD HD
  • CH the volume of capillary pores
  • cp aluminates
  • CK clinker grains
  • the phase of the aluminates solid phase
  • the expressions of the homogenized elastic modules are as follows:
  • the parameters on this Ech.2 scale are, for example, the elastic properties of CH and the initial reactive phases of the cement (for example the four main components of the clinker).
  • the invention relates to a computer program product for keeping in a memory of a unit of a processor or on a removable memory medium suitable for cooperating with said processor unit, the product- program comprising instructions for implementing the method according to the invention
  • the exemplified cementitious system below includes Class G Portland cement, the composition of which is given in Table 3. TABLE 3
  • the initial composition Co of the cement system ceml is given in Table 4 below. It corresponds to a water-cement ratio equal to about 044 (w / c in English).
  • the manufacture of the cement slurry is carried out as follows, by mixing in 5 steps:
  • the pressure is first increased to the desired value (for example 0.3 MPa, 20 MPa or 40 MPa) and the temperature is then increased up temperature of the test in 30 minutes.
  • FIG. 4 shows the velocity of compression waves as a function of time V p (t) obtained.
  • Figure 4 shows the degree of hydration resulting from this calorimetric test.
  • the curves a (t) and V p (t) merge almost completely, which shows that a linear relationship exists between the evolution over time of the hydration degree a (t) of the cement system, measured by calorimetry. , and the evolution over time of the compression wave velocity V p (t).
  • FIG. 5 shows that the relation between a (t) and V p (t) is almost linear.
  • Steps A and B of the process determination of a) Pu Tu and different fineness values ⁇ "of the cement system, pressure P" and temperature T "
  • the velocity of compressional waves as a function of time V p (t) is measured in a sample of the ceml cementitious system, at a pressure Pi and a temperature T ls as described in the UCA test above.
  • FIGS. 6 to 8 which compare experimental data resulting from UCA tests as described above, and the data resulting from the determination of the degree of hydration according to steps A and B of the process which is the subject of the invention (curves in dotted lines, indicated as “simulation” in the legend), show the good predictive ability of steps A and B of the method which is the subject of the invention.
  • Step C of the process determination of the composition C (t)
  • stages 1 and 2 of the hydration process including the initial and dormant hydration phases, are not taken into account.
  • Stages 3 and 4 including the acceleration and deceleration phases of hydration, are taken into account and simulated according to the Avrami 1 nucleation and growth theory.
  • Stage 5 comprising the curing phase is also taken into account, and is considered to be controlled by diffusion phenomena.
  • FIG. 9 represents the evolution of the volume fractions of the different phases of the cementitious system cem1 as a function of the degree of hydration, for a hydration temperature of 20 ° C. and at atmospheric pressure.
  • Table 8 presents the incompressibility modules k and shear modulus g of the different solid components taken into account during the multiscale analysis, according to step D of the process: the four phases of the clinker of the ceml system, the phase hydrated CH and the globules representing the CSH solid phase of the cement system ceml. All these values were obtained from the literature, except for the cell modules that were obtained by calibration with the Bourissai test data (21) .
  • Multi-scale analysis includes two homogenizations:
  • the first homogenization considers, on the one hand, the matrix formed by the solid CSH globules, and on the other hand the pores in the CSH LD and CSH HD phases. It makes it possible to calculate the incompressibility and shear modules of the "porous" CSH LD and CSH HD.
  • the characteristic length is about 10 "9 - 10 " 8 m.
  • the matrix occupies 63% of volume for the case of the CSH phase LD, and it occupies a proportion greater than 76% o, depending on the temperature, volume for the CSH HD phase;
  • the second homogenization considers, on the one hand, the solid phases (CSH LD, C-HSH, CH, aluminates and initial reactive phases of the clinker) and, on the other hand, the capillary porosity. It makes it possible to calculate the incompressibility and shearing modules of the cementitious system and subsequently all of the elastic and hydro-mechanical parameters.
  • the characteristic length is about 10 "8 m. TABLE 8
  • FIGS. 10 to 13 show a good correspondence between the mechanical parameters determined according to the method of the invention, and the experimental data obtained in the literature (Boumiz et al 1 - *; Haecker et al 2 - *), thus validating the process according to the invention.
  • Data from the literature has been corrected to transform dynamic modules into static modules using the following formula: dynamic static
  • the measured data ("E-measurements” in the legend, cross-shaped symbols) and determined by the method according to the invention of the Young's modulus (E) are compared for a temperature of cure of 60 ° C ("Ehom” in the legend).
  • the incompressibility (K) and shear (G) modules determined by the method according to the invention are also represented ("Khom” and "Ghom” in the legend).
  • the shear modulus of the cement paste hydrated at 60 ° C. at the age of 1000 h (42 days) is 6.1 gpa.
  • the difference between this value and the measured value (5.9GPa) is 3.3%.

Abstract

La présente invention concerne un procédé pour la détermination de paramètres mécaniques d'un système cimentaire, en fonction du temps, et en fonction de la finesse du système cimentaire, de la pression et/ou de la température, représentatives des conditions in situ rencontrées dans les puits de forage. La composition initiale du système cimentaire, sa finesse Φ et la vitesse des ondes de compression en fonction du temps Vp(t) sont les seules données d'entrée du procédé. Celui-ci comprend : ■ une étape A d'estimation du degré d'hydratation du système cimentaire en fonction du temps α(t) à partir de Vp(t), à une pression Pl et une température Tl; ■ une étape B dans laquelle α(t) est déterminé en fonction de valeurs souhaitées de finesse Φn du système cimentaire, de pression Pn et/ou de température Tn; ■ une étape C dans laquelle la composition du système cimentaire en fonction du temps C(t) et en fonction de valeurs souhaitées de finesse Φn du système cimentaire, de pression Pn et/ou de température Tn à partir de α(t) déterminé à l'étape B; ■ et une étape D de détermination d'au moins un paramètre mécanique du système cimentaire en fonction du temps et en fonction des valeurs souhaitées de finesse Φn du système cimentaire, de pression Pn et/ou de température Tn, à partir de C(t) déterminée à l'étape C. Selon le procédé de l'invention, ces paramètres peuvent être déterminés dès le plus jeune âge du système cimentaire. En particulier, les paramètres de déformabilité statique et des paramètres de couplage hydro -mécaniques sont déterminés par le procédé selon l'invention.

Description

METHODE DE CARACTERISATION DU COMPORTEMENT MECANIQUE DE
CIMENTS
DOMAINE DE L'INVENTION
Le domaine de l'invention est celui de la caractérisation du comportement mécanique des ciments, en particulier des ciments utilisés lors des opérations de cimentation des puits de forage, par exemple dans le cadre de l'exploration et la production d'hydrocarbures. En particulier, l'invention concerne un procédé de détermination de paramètres mécaniques d'un système cimentaire de composition donnée, en fonction du temps, et en fonction de la finesse du système cimentaire, de la pression et/ou de la température.
ETAT DE LA TECHNIQUE
La cimentation d'un cuvelage dans un puits pétrolier consiste à placer une gaine de ciment dans l'annulaire entre l'extrados du cuvelage et la paroi du trou, le trou pouvant être constitué par un autre cuvelage ou par de la roche. Un bouchon de ciment peut également être créé dans le puits dans le cadre d'opérations de forage ou de complétion, ou au moment du bouchage définitif de celui-ci. Cette gaine ou ce bouchon de ciment ont un rôle primordial dans la stabilité et l'isolation des puits pétroliers.
La gaine ou le bouchon de ciment sont obtenus par pompage d'un laitier ou coulis de ciment fabriqué à partir de ciment, d'eau et d'adjuvants. Ce laitier de ciment est à l'état liquide lorsqu'il est pompé. L'hydratation des particules de ciment conduit le laitier liquide vers un état solide, caractérisé par l'existence d'un squelette et de pores, formant ainsi un milieu poreux.
La gaine de ciment est exposée, durant la vie du puits et sur toute sa longueur, à diverses sollicitations mécaniques et thermiques (conditions in situ) provenant d'opérations menées en surface et dans le puits (tests en pression, changement de boue, stimulations froides et chaudes, production des réserves...) ou de phénomènes prenant naissance directement dans le sous-sol (compaction du réservoir, séismes...), et ce jusqu'à son abandon, voire au-delà. Le bouchon de ciment est également soumis à diverses sollicitations mécaniques et thermiques. Ces sollicitations peuvent endommager le matériau constitutif de la gaine ou du bouchon de ciment, dégrader ses propriétés mécaniques et hydrauliques et, par conséquent, modifier sa contribution à la stabilité et à l'étanchéité du puits.
La connaissance du comportement mécanique du ciment dans les conditions in situ et de l'évolution de ce comportement au cours du temps est essentielle pour l'analyse du fonctionnement du puits pendant son forage, son exploitation, et son bouchage définitif. Par exemple, il est important de connaître le comportement mécanique du ciment pour garantir l'étanchéité du puits lors de ces opérations, lors du stockage et de la séquestration de fluides dans des réservoirs souterrains, tels que des gaz à effet de serre, comme le C02. Cette connaissance du comportement mécanique du ciment et de son évolution dans le temps, permet notamment de simuler, par le biais de modèles analytico-numériques, la vie des gaines ou des bouchons de ciment. Les conditions in situ influent sur le processus de formation du ciment durci, et par conséquent sur son comportement mécanique. Les ciments obtenus par durcissement de laitiers dans les conditions très spécifiques que sont celles rencontrées dans les puits sont en général très différents des ciments du même type formés dans des conditions ambiantes (c'est-à-dire à l'air ambiant et à pression atmosphérique). Dans les cas extrêmes, pour des puits de 6500 mètres de profondeur, forés avec des boues de densité égale à environ 2, la pression de fond peut atteindre plus de 130 MPa et la température plus de 250°C.
Il existe donc un besoin de connaître les paramètres physiques et mécaniques du ciment et leur évolution dans le temps, à différentes conditions de pression et température, analogues à celles rencontrées le long du puits, et plus particulièrement aux endroits les plus sollicités en pression et/ou en température. Il est à ce titre intéressant d'évaluer ces paramètres après la mise en place de la gaine ou du bouchon de ciment, c'est-à-dire dès le plus jeune âge du ciment lorsque l'hydratation du ciment commence à former une structure, aussi appelée squelette, et ce jusqu'à un temps infini correspondant à un ciment complètement hydraté, dans lequel les processus d'hydratation sont achevés (stabilisation de la composition du ciment, notamment de la quantité d'hydrates).
De nombreuses techniques de mesure sur des échantillons de ciment ont été proposées pour caractériser le comportement mécanique de tels matériaux. Les principaux paramètres mécaniques déterminés à partir de ces mesures sont des paramètres de rupture, tels que la résistance à la compression ou la résistance à la traction (uniaxiale ou triaxiale), et des paramètres de déformabilité, en particulier des paramètres élastiques du matériau tels que le module de Young ou le coefficient de Poisson.
Une première catégorie de techniques recouvre les essais mécaniques statiques (chargement en compression/traction uniaxiale, triaxiale, essai œdométrique etc.) sur des échantillons qui sont durcis à l'intérieur d'un moule dans des bancs de vieillissement, en pression et température, puis qui sont déchargés afin d'être positionnés dans un appareil de mesure. Ces essais nécessitent de ramener les échantillons à pression atmosphérique et à température ambiante, afin de les équiper de capteurs de mesure, de les positionner sous une presse, puis de réaliser les essais proprement dits, après une éventuelle remise des échantillons dans des conditions de température et de pression analogues à celles rencontrées dans le puits. Ces cycles de chargement/déchargement peuvent non seulement endommager les échantillons préalablement à la mesure, mais aussi perturber la mesure des caractéristiques desdits échantillons. En outre, cette première catégorie de techniques ne permet pas de tester les ciments de très jeune âge, car il est nécessaire d'attendre que l'échantillon de ciment puisse être retiré de son moule afin d'être testé. Enfin, de par leur nature destructive, ces techniques ne permettent qu'une mesure ponctuelle dans le temps, et n'offrent pas la possibilité d'accéder à l'évolution des paramètres mécaniques dans le temps, en fonction de l'hydratation du ciment.
Une deuxième catégorie de techniques recouvre les essais dynamiques fondés sur une mesure de la propagation d'ondes ultrasoniques, et ne comprenant pas de retour des échantillons à des conditions ambiantes. Ces techniques ont cependant un intérêt limité de par leur caractère indirect. En particulier, les paramètres statiques, qui sont les paramètres utilisés pour simuler le comportement des gaines ou des bouchons de ciment durant leur vie, ne sont accessibles qu'indirectement : ils sont calculés à partir des paramètres dynamiques, déterminés à partir de la mesure des ondes ultrasoniques, en utilisant des formules de corrélation ; ces formules sont elles-mêmes obtenues à partir d'essais statiques, qui peuvent être entachés d'erreur, voire ne pas couvrir le domaine d'application des matériaux testés. En effet, ces essais statiques sont généralement réalisés sur des échantillons cubiques de ciment dans certaines conditions de pression et de température qui ne correspondent pas nécessairement aux conditions des essais dynamiques. Les mesures de type Ultrasonic Cernent Analyser UCA font partie de ces techniques. Une mesure de type UCA consiste à mesurer la vitesse des ondes de compression et à utiliser une corrélation pour évaluer la résistance à la compression du ciment.
Une troisième catégorie de techniques comprend quelques propositions d'essais mécaniques statiques sans étape de chargement/déchargement telle que mentionnée ci-dessus qui nécessitée de ramener l'échantillon de ciment à pression et température ambiante avant d'effectuer une mesure.
Ainsi, le document EP 1541987 décrit un système dans lequel une composition de ciment est coulée dans un moule ayant une forme d'os, où l'échantillon est vieilli sous température et sous pression, et où un chargement en traction uniaxiale est effectué jusqu'à la rupture de l'échantillon, sans avoir à décharger l'échantillon. Toutefois, cette méthode ne permet pas d'effectuer les mesures en conditions in situ, la pression ne pouvant être exercée que sur deux faces de l'échantillon, les autres faces subissant une condition de chargement par réaction du moule et non par application de contraintes en conditions in situ. Les mesures sont donc biaisées. En outre, seuls des tests de traction sont possibles, or ceux-ci sont biaisés en ce qui concerne la mesure des constantes élastiques par rapport aux tests de compression, du fait de l'apparition de microfissures qui invalident l'hypothèse d'élasticité. La plage de détermination des paramètres élastiques est donc très réduite. De plus il n'est pas possible de mesurer les paramètres de rupture en compression, ni de mesurer les paramètres au jeune âge du ciment. Enfin, la géométrie utilisée est non conventionnelle.
Le document US 7,621,186 décrit une variante du système précédent, adoptant une géométrie de type tronconique. Il souffre donc des mêmes inconvénients. Le document WO 2007/020435 propose une technique qui consiste à faire prendre la composition de ciment dans un annulaire situé entre deux tubes concentriques, puis à faire varier les pressions en intrados du tube intérieur et/ou en extrados du tube extérieur tout en mesurant les déformations induites. Cette technique ne permet pas de réaliser des essais de compression axiale. En outre, cette technique présente l'inconvénient de se fonder sur une mesure en champ de contraintes hétérogène (en élasticité, les champs de contraintes et déformation dans un cylindre creux varient en 1/r2). Ainsi la mesure des propriétés élastiques de l'échantillon est très imprécise (très sensible aux erreurs), tout comme celle des propriétés d'endommagement et de rupture de l'échantillon.
Le document US 7,089,816 décrit une technique qui consiste à faire prendre la composition de ciment dans une enveloppe cylindrique, constituée d'une membrane déformable et de deux pistons, et placée dans une enceinte de confinement. Il est ensuite directement procédé aux essais mécaniques en appliquant une pression de confinement au travers de la membrane et un chargement axial par l'intermédiaire des pistons, comme pour une cellule triaxiale classique. Un inconvénient de cette technique est que l'utilisation d'une membrane souple pour la prise du ciment aboutit à un échantillon de forme irrégulière à l'issue de la prise. Du fait des variations de volume associées à la prise, il se produit en effet des instabilités de Taylor conduisant l'échantillon à perdre sa géométrie initiale. Aussi, cette technique ne permet pas d'effecteur des mesures en accord avec les procédures existantes, l'hydratation du ciment n'étant pas correctement reproduite.
Le document US 7,549,320 correspond à une technique du même type, avec une variation de la technologie de chargement. La prise de l'échantillon s'effectue, comme dans US 7,089,816, à l'intérieur d'une membrane flexible. L'enceinte de confinement rigide entourant la membrane flexible est compartimentée grâce à des dispositifs d'étanchéité actifs, et permet l'application d'une pression différentielle à l'échantillon, par l'injection de fluides. La membrane flexible utilisée dans ce dispositif peut en outre être semi-perméable, limitant ainsi la déformation de l'échantillon durant sa prise, mais perturbant le processus d'hydratation de l'échantillon de ciment.
Le document US 7,552,648 décrit encore une autre variante, dans laquelle on injecte un fluide dans l'échantillon lui-même, qui est poreux, afin d'obtenir la pression souhaitée. On effectue ensuite un essai de traction. Aucun test de compression n'est prévu, et l'apport de fluide extérieur ne simule pas correctement les échanges hydriques en conditions in situ, comme c'est aussi le cas pour la technique selon le document US 7,549,320.
Outre les inconvénients mentionnés ci-dessus pour cette troisième catégorie de techniques, ce type d'essais ne permet pas d'effectuer des mesures sur des ciments de très jeune âge, c'est-à-dire dès les premières heures de l'hydratation, ni de suivre l'évolution des paramètres mécaniques mesurés dans le temps, en fonction de l'hydratation du ciment. Aucune des trois catégories de techniques décrites ne permet de mesurer sans artéfacts les propriétés mécaniques d'un ciment en cours de prise.
OBJECTIFS DE L'INVENTION
Il existe par conséquent un besoin de disposer d'un procédé de détermination des paramètres mécaniques d'un système cimentaire et de leur évolution dans le temps, dans des conditions in situ, qui ne comporte pas les inconvénients des techniques de mesure décrites ci- dessus.
L'invention vise à répondre à un besoin de caractérisation simple du comportement mécanique d'un système cimentaire, à travers la détermination de propriétés mécaniques statiques et de propriétés hydro-mécaniques exprimées en fonction du temps, et en fonction de différentes pressions et températures correspondant aux conditions rencontrées le long d'un puits de forage.
Un autre objectif de l'invention consiste à fournir un procédé facilement industrialisable pour caractériser le comportement mécanique d'un système cimentaire, qui utilise des mesures simples et normalisables.
Un autre objectif de l'invention consiste à fournir un procédé permettant d'estimer l'évolution des paramètres mécaniques d'un système cimentaire au cours du temps, du début à la fin de son hydratation.
L'invention vise également à fournir une estimation de paramètres mécaniques statiques du ciment, qui constituent les paramètres utilisés pour simuler de manière fiable le comportement mécanique des gaines ou des bouchons de ciment durant leur vie.
BREVE DESCRIPTION DES DESSINS
La figure 1 représente le logigramme du procédé selon l'invention.
La figure 2 est un schéma illustrant les différents stades du processus d'hydratation d'un ciment Portland.
La figure 3 présente un exemple de dispositif de type Ultrasonic Cernent Analyser (UCA) permettant de mesurer la vitesse des ondes de compression dans un échantillon de ciment utilisée dans le procédé selon l'invention.
La figure 4 est un diagramme représentant la vitesse des ondes de compression en fonction du temps Vp(t) mesurée lors d'un essai de type UCA sur un système cimentaire comprenant du ciment Portland, ainsi que le degré d'hydratation en fonction du temps a(t) du même système cimentaire mesuré lors d'un essai calorimétrique. La figure 5 est un diagramme illustrant la relation linéaire entre la vitesse des ondes de compression en fonction du temps Vp(t) et le degré d'hydratation en fonction du temps a(t) tels que représentés à figure 3.
Les figures 6 à 8 sont des diagrammes comparant le degré d'hydratation en fonction du temps a(t) déterminé par le procédé selon l'invention pour un système cimentaire donné, et le degré d'hydratation en fonction du temps a(t) pour ce même système, déterminé expérimentalement. Les degrés d'hydratation sont représentés pour différentes valeurs de pression et de température. Dans la figure 6, a(t) est représenté pour une pression faible de 0.3 MPa et à différentes températures (7°C, 13°C, 25°C, 40°C et 60°C). Dans la figure 7, a(t) est représenté à des pressions plus élevées, selon les couples pression/température suivants : (40MPa, 30°C), (40MPa, 15°C), (20MPa, 30°C). Sur la figure 8, a(t) est représenté pour les deux couples pression/température suivants : (0.3MPa, 60°C), (20MPa, 60°C),
La figure 9 est un diagramme illustrant la composition d'un système cimentaire donné en fonction du degré d'hydratation du système, déterminé par le procédé selon l'invention. En particulier, la fraction volumique des différents constituants du système cimentaire en fonction du degré d'hydratation du système cimentaire est donnée, pour un couple de pression et de température fixé.
Les figures 10 à 13 sont des diagrammes comparant les paramètres mécaniques en fonction du temps déterminés par le procédé selon l'invention, pour un système cimentaire donné, à ceux déterminés expérimentalement et issus de la littérature (Boumiz et al 1-*, Haecher et al 2-*) pour ce même système cimentaire.
DESCRIPTION DETAILLEE DE L'INVENTION
Dans la description de l'invention et les exemples particuliers de l'invention qui suivent, il est fait référence aux dessins annexés.
Pour atteindre au moins l'un des objectifs susvisés, parmi d'autres, la présente invention propose un procédé de détermination de paramètres mécaniques d'un système cimentaire de composition initiale Co et de finesse Φ, en fonction du temps, et en fonction de la finesse du système cimentaire, de la pression et/ou de la température. Ce procédé comprend les étapes suivantes :
(A) déterminer le degré d'hydratation du système cimentaire en fonction du temps a(t) à partir de la vitesse des ondes de compression en fonction du temps Vp(t) mesurée dans un échantillon du système cimentaire, à une pression Pi et une température ΤΊ ;
(B) déterminer le degré d'hydratation a(t) en fonction de valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn ; (C) déterminer la composition du système cimentaire en fonction du temps C(t) et en fonction de valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn, à partir du degré d'hydratation a (t) déterminé à l'étape B ;
(D) déterminer au moins un paramètre mécanique du système cimentaire en fonction du temps, et en fonction de valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn, à partir de la composition du système cimentaire C(t) déterminée à l'étape C.
Selon un mode de réalisation de l'invention, le procédé permet de caractériser le comportement mécanique d'un système cimentaire utilisé en tant que gaine ou bouchon de ciment dans un puits, de préférence un puits pétrolier. Le procédé selon l'invention ne se limite pas à cette utilisation, et peut être mis en œuvre dans le cadre de toute construction de puits comprenant une gaine ou un bouchon de ciment pour lequel/laquelle la caractérisation du comportement mécanique du système cimentaire utilisé est souhaitée, par exemple dans le domaine du stockage de gaz ou de la géothermie.
Définitions
Dans tout le présent exposé, sauf indication expresse contraire, un singulier doit être interprété comme un pluriel et réciproquement.
Les définitions qui suivent sont données à titre d'exemples, pour faciliter l'interprétation du présent exposé.
Un système cimentaire au sens de la présente invention correspond à un laitier de ciment ou un ciment durci.
Par laitier de ciment, on entend une composition fluide (liquide, pâteuse, granulaire...) à base de ciment et d'eau, apte à durcir, c'est-à-dire susceptible de passer à un état solide ou quasiment solide au cours du temps (en subissant une étape de durcissement). L'expression « coulis de ciment » est utilisée en tant que synonyme de laitier de ciment. Ainsi, le durcissement correspond dans ce cas essentiellement à l'hydratation (ou prise) de la composition à base de ciment et d'eau.
Par ciment durci, on entend un laitier de ciment qui a subit un durcissement, et dans lequel a commencé à se développer un squelette poreux, pouvant lui conférer une résistance mécanique. Au sens de l'invention, un ciment durci n'est pas nécessairement un ciment complètement durci, c'est-à-dire dont l'hydratation est achevée. Un ciment complètement durci est un ciment dont l'hydratation est complète à 100%.
Par paramètre mécanique statique, on entend un paramètre qui n'est pas obtenu à partir de la mesure de vitesses d'ondes. Les paramètres statiques constituent les paramètres utilisés pour simuler de manière fiable le comportement mécanique des gaines ou des bouchons de ciment durant leur vie, à l'inverse des paramètres dynamiques. Ces derniers, généralement issus de méthodes de mesures d'ondes soniques/ultrasoniques, divergent des paramètres statiques, en raison de la dépendance des propriétés mécaniques du ciment à la fréquence des ondes utilisées et de la dimension des échantillons de ciment testés.
Par finesse d'un ciment, on entend la finesse de broyage d'un ciment exprimée par sa surface spécifique (surface développée par unité de masse, exprimée en cm2/g), et mesurée à l'aide d'un appareil de type Blaine selon la norme NF EN 196-6.
Par ondes de compression ou ondes P, on entend des ondes soniques ou ultrasoniques crées par la vibration mécanique d'un support fluide ou solide et propagée grâce à l'élasticité du milieu environnant sous forme d'ondes longitudinales.
Par les termes initial (aux), initiale(s), utilisés par exemple en référence à la composition du système cimentaire, à la température, la pression ou tout autre paramètre physico-chimique dont il est question dans la présente description, il est fait référence au stade initial du système cimentaire sous forme d'un laitier de ciment, et dans lequel les réactions d'hydratation n'ont pas ou quasiment pas débutées, correspondant au temps initial t=0.
Par phase initiale réactive, on entend une phase initiale du ciment du système cimentaire, apte à s'hydrater pour donner une phase hydratée, aussi appelée hydrate dans la présente description.
Le système cimentaire
Le procédé selon la présente invention s'applique à l'ensemble des ciments dont l'hydratation se fait au travers d'une série de réactions chimiques exothermiques entre ses composants et l'eau, par exemple les ciments Portland, les ciments alumineux, et autres ciments définis selon les normes EN-196-2, EN-196-6, EN196-7, ISO 3310-1, ISO 13500, et les normes ISO 10426-1 à ISO 10426-6 en ce qui concerne l'industrie pétrolière.
Le système cimentaire selon la présente invention peut comprendre des adjuvants, permettant d'obtenir les propriétés désirées :
■ des accélérateurs (exemple : le chlorure de calcium), qui ont pour objectifs de réduire la période d'induction du système cimentaire, et éventuellement d'accélérer le développement de la résistance à la compression uniaxiale ;
■ des retardateurs (exemple : les lignosulfonates), qui ont pour objectif d'augmenter la période d'induction du système cimentaire ;
■ des extendeurs (exemple : la bentonite, des billes creuses), qui ont pour objectif de diminuer la densité du laitier de ciment et/ou réduire la concentration en ciment ou augmenter la concentration d'eau ;
■ des alourdisseurs (exemple : l'hématite) qui ont pour objectif d'augmenter la densité du laitier de ciment ; ■ des dispersants (exemple : les lignosulfonates), qui ont pour objectif de diminuer la viscosité du laitier de ciment en défloculant les grains de ciment ;
■ des contrôleurs de filtrat (exemple : le latex), qui ont pour objectif de limiter les transferts d'eau du laitier de ciment vers le système extérieur par construction d'un cake peu ou pas perméable contre la formation ;
■ des contrôleurs de pertes (exemple : la gilsonite), qui ont pour objectif de limiter les pertes du ciment dans les zones très perméables ;
■ du sel ;
■ des additifs spéciaux, comme des anti-mousses, des fibres, des agents d'expansion, des agents anti-gaz, des agents anti-sédimentation, ou des agents donnant des ciments résilients.
Le système cimentaire selon la présente invention comprend une composition initiale Co à l'instant t = 0. Cette composition initiale Co comprend les phases initiales du ciment, parmi lesquelles au moins une phase initiale réactive pouvant s'hydrater, et l'eau initiale, selon des fractions données, par exemples exprimées en volume. A l'instant t, l'hydratation du système cimentaire est avancée, et la composition du ciment C(t) est différente de la composition initiale Co. La composition du ciment C(t) comprend m phases, parmi lesquelles les phases initiales du ciment dans des proportions différentes de Co, au moins une phase hydratée issue de l'hydratation d'au moins une phase initiale réactive du ciment, et de l'eau.
Selon un mode de réalisation de l'invention, le ciment du système cimentaire comprend au moins une phase initiale réactive X, de préférence choisie dans le groupe constitué par le silicate tricalcique de formule 3CaO.Si02 (Alite), de symbole C3S, le silicate bicalcique de formule 2CaO.Si02 (Belite), de symbole C2S, Paluminate tricalcique de formule 3CaO.Al203 (Aluminate), de symbole C3A, l'alumino ferrite tétracalcique C4F, et leurs associations. Selon ce mode de réalisation, le ciment du système cimentaire comprend au moins une phase hydratée issue de l'hydratation d'au moins une phase initiale réactive du ciment, de préférence choisie dans le groupe constitué par le silicate de calcium hydraté C-S- H, l'hydroxyde de calcium CH, le trisulfoaluminate de calcium hydraté TSA (Ettringite), le monosulfate de calcium hydraté AFm, l'alumino ferrite de calcium hydraté, et leurs associations.
Les phases initiales telles que les silicates C3S, C2S et les aluminates sont des phases réactives qui réagissent avec l'eau du système cimentaire pour donner des phases hydratées, par le biais de réactions chimiques exothermiques. Les silicates réagissent avec l'eau pour former des silicates de calcium hydratés C-S-H, et de l'hydroxyde de calcium CH (Portlandite de formule Ca(OH)2). Les aluminates réagissent avec l'eau et le gypse, ou avec l'eau seule pour former les aluminates hydratés suivants : le trisulfoaluminate de calcium hydraté TSA (Ettringite), le monosulfate de calcium hydraté AFm et l'alumino ferrite de calcium hydraté. La microstructure du système cimentaire, en particulier du ciment durci, est très complexe. Les fractions volumiques, la possibilité de cristallisation, la morphologie et les dimensions des phases hydratées du ciment durci sont présentées dans le tableau 1 ci-dessous, relatif à un ciment Portland (Bernard et al. ^12-*).
TABLEAU 1
Figure imgf000012_0001
La Porlandite CH est formée de cristaux hexagonaux larges qui sont entourés par la phase C-S-H, qui est sous forme de gels, et occupe environ 20% du volume du ciment durci. Les monosulfoaluminates et PEttringite jouent un rôle mineur dans la structure du ciment durci et représentent généralement une fraction volumique de 15% à 20%. La phase C-S-H constitue le produit majeur du ciment durci. La phase C-S-H est une phase poreuse, de structure amorphe et colloïdale, et présente une composition chimique variée. Elle occupe entre 50%> et 70%> du volume total du ciment durci, ce qui donne aux propriétés de cette phase un rôle dominant dans les propriétés macroscopiques du système cimentaire. La phase C-S-H existe sous deux formes : la phase "haute densité" C-S-H HD (aussi dénommée « HD » dans la présente description) et la phase "faible densité" C-S-H LD (aussi dénommée « LD » dans la présente description). Chacune de ces deux formes possède, selon Jennings et al.(8'9), une unité commune appelée globule. La différence entre ces deux types de C-S-H consiste en l'arrangement des globules, résultant en différentes porosités selon le type de C-S-H : environ 24% pour le C-S-H HD, avec une taille de pores comprise entre environ lOnm et lOOnm (Constantinides<:i0)), et environ 37% pour le C-S-H LD, avec une taille de pores comprise entre environ 5nm et 50nm
Selon un mode de réalisation de l'invention, le système cimentaire comprend du ciment Portland, qui constitue l'un des ciments les plus courants, comprenant les phases initiales suivantes utilisés (voir les normes EN-196-2, EN-196-6, EN196-7, ISO 3310-1, ISO 13500, et les normes ISO 10426-1 à ISO 10426-6 pour les teneurs en phases initiales, dépendant de la classe et du grade du ciment) :
- du silicate tricalcique de formule 3CaO.Si02 (Alite), de symbole C3S ;
du silicate bicalcique de formule 2CaO.Si02 (Belite), de symbole C2S ; de Paluminate tricalcique de formule 3CaO.Al203 (Aluminate), de symbole C3A ; du ferro-aluminate tétracalcique de formule 4CaO.Al203 Fe203 (Ferrite), de symbole C4F ;
du sulfate dihydraté de calcium CaS04.2H20 (gypse).
Les phases silicatées et les aluminates sont issus d'un clinker broyé, principal constituant du ciment Portland.
Il est bien entendu que le procédé selon l'invention ne se limite pas au système cimentaire comprenant un ciment Portland. En particulier, le procédé peut être adapté pour la détermination de paramètres mécaniques de systèmes cimentaires variés, présentant des compositions chimiques et des réactions d'hydratation différentes de celles d'un système cimentaire comprenant un ciment Portland, sans pour autant sortir du cadre de l'invention.
Le processus d'hydratation du système cimentaire Le procédé selon l'invention prend en compte différentes phases du processus d'hydratation d'un système cimentaire pour la détermination d'au moins un paramètre mécanique de celui-ci. Le processus d'hydratation peut être suivi par l'évolution du flux de chaleur crée par les réactions exothermiques d'hydratation des phases initiales réactives du ciment avec l'eau.
La quantité de chaleur dégagée par les réactions exothermiques lors de l'hydratation du système cimentaire peut être mesurée en utilisant un calorimètre. Deux types de calorimètre sont généralement utilisés. Le premier type correspond aux calorimètres quasi adiabatiques, dans lesquels la réaction a lieu à l'intérieur d'une enceinte parfaitement isolée thermiquement. Ce type de calorimètre n'est pas adapté à la mesure d'échantillons de systèmes cimentaires sous pressions élevées, et la pression limite d'utilisation est d'environ 0,1 MPa. Un deuxième type correspond aux calorimètres isothermes, avec lesquels les essais sont réalisés dans un espace de l'ordre du centimètre cube, et où la température est contrôlée afin que la température de l'espace ne varie pas. Les matériaux utilisés pour ce type d'appareil de mesure permettent de travailler avec des pressions plus élevées, de l'ordre du Méga Pascal.
En général, plusieurs stades peuvent être identifiés au cours du processus d'hydratation d'un système cimentaire. Classiquement, l'hydratation du système cimentaire comprend dans cet ordre :
- une période de dissolution, au cours de laquelle se produit la dissolution de l'Ettringite et la formation de gel CSH ;
- une période dormante, dite aussi période d'induction, au cours de laquelle se produit une augmentation de la concentration en C3A et en ions OH" ;
- une période de formation rapide de C-S-H et de Portlandite CH, dite aussi période d'accélération ; une période de ralentissement de la formation de C-S-H et de Portlandite CH ; et de formation éventuelle de monosulfoaluminate ;
et une période de cure, dite aussi période de diffusion, où les réactions chimiques sont contrôlées par des phénomènes de diffusion.
La figure 2 illustre le processus d'hydratation dans le cas d'un ciment Portland (Rixom et al.(15)). Dans ce cas, le processus d'hydratation est généralement divisé en 5 étapes. Chaque étape correspond à une tendance d'évolution de la chaleur provoquée par les réactions d'hydratation :
- Période 1 : une période d'hydratation initiale se déroule pendant environ 15 minutes, durant laquelle les composants chimiques les plus solubles (Na2S04, K2S04, CaS04.xH20) du ciment se dissolvent, faisant augmenter le pH du laitier de ciment.
- Période 2 : une période dormante s'ensuit, pendant environ 4 heures. La chaleur d'hydratation n'évolue presque pas. La résistance de la pâte de ciment n'est pas encore créée.
- Période 3 : cette troisième période correspond à une accélération de l'hydratation, qui dure environ 4 heures. Au cours de celle-ci, la chaleur d'hydratation évolue très rapidement.
Les hydrates formés se mettent en contact, ce qui forme la résistance de la pâte.
- Période 4 : une période de ralentissement s'opère ensuite et dure quelques jours. L'évolution de la chaleur d'hydratation ralentit.
- Période 5 : une période de cure se manifeste pendant plusieurs mois après le début de la prise du système cimentaire, et l'évolution de la quantité des hydrates se stabilise progressivement.
L'évolution de l'hydratation du système cimentaire dans le temps est classiquement décrite par l'évolution du degré d'hydratation dans le temps a(t), qui est un paramètre intrinsèque, et dont la valeur varie de 0 à 1. Celui-ci, à un instant t, peut être évalué comme le rapport de la quantité de chaleur Q(t) produite par les réactions chimiques à l'instant t sur la quantité totale de chaleur Q(t=∞) produite lorsque les réactions chimiques ont toutes eu lieu. La cinétique d'hydratation varie au cours du processus d'hydratation. Parmi les périodes du processus d'hydratation, un premier stade où l'hydratation est gouvernée principalement par un processus de nucléation et de croissance, tel que décrit par exemple par Avrami(3), est classiquement reconnu. Un deuxième stade du processus d'hydratation du système cimentaire est caractérisé par une hydratation gouvernée principalement par un processus diffusion des ions. Ce deuxième stade commence lorsque le degré d'hydratation a atteint une valeur d'un degré d'hydratation seuil a*. Lors du premier stade, gouverné par les phénomènes de nucléation et de croissance, une couche d'hydrates se crée autour des grains de ciment, et est relativement perméable, ce qui permet la diffusion des ions et l'avancement des réactions chimiques. Lorsqu'un degré d'hydratation seuil a* est atteint, une diminution de la diffusion des ions se produit, liée à l'épaississement de la couche d'hydrates autour du grain de ciment et à la baisse de perméabilité de cette couche. C'est alors la diffusion des ions qui contrôle l'avancement des réactions (Kondo et Kodama ' ; Fuji et Kondo '). Dans le cas du processus d'hydratation d'un ciment Portland, tel que représenté à la figure 2, le degré d'hydratation seuil a* marque le passage de la période 4 à la période 5. Les périodes 3 et 4 correspondent aux phases de nucléation et de croissance selon Avrami(3), et la période 5 (période de cure) est principalement gouvernée par le processus de diffusion des ions.
Le silicate tricalcique C3S et le silicate bicalcique C2S réagissent avec l'eau selon les équations suivantes, issues de Jennings et al.(6):
2C,S + 10,6H→C, 4 - S - H„+ 2,6CH
3 3'4 8 (I)
2C2S + 8,6H→ C3 4 - S2 - H8 + 0,6CH
Les phases hydratées C-S-H et le CH forment les principaux constituants du ciment durci, et sont principalement responsables de la solidité du ciment durci. La cinétique d'hydratation de la phase C3S est plus rapide que celle de la phase C2S, et la quantité de CH créée par l'hydratation du C3S est environ trois fois plus grande que celle crée par l'hydratation de C2S. Ainsi, la phase C3S joue un rôle majeur dans le processus d'hydratation et dans le développement du comportement mécanique du système cimentaire, en initiant la résistance de la pâte de ciment pendant les premiers jours.
Selon un mode de réalisation de l'invention, le procédé prend en compte les deux réactions (I) décrites ci-dessus pour simuler le processus d'hydratation.
Nature des paramètres mécaniques déterminés par le procédé selon l'invention
L'invention vise à caractériser le comportement mécanique d'un système cimentaire, et notamment à déterminer certains paramètres mécaniques et leur évolution dans le temps, dans des conditions in situ, qui pourront ensuite être utilisés, par exemple, dans la modélisation du comportement des gaines ou des bouchons de ciment des puits de forage durant leur vie. Le procédé selon l'invention permet de déterminer au moins un paramètre mécanique du système cimentaire en fonction du temps, pour des valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn.
Selon un mode de réalisation de l'invention, le procédé permet de déterminer au moins un paramètre mécanique qui est choisi parmi les paramètres de déformabilité statiques, de préférence les paramètres élastiques statiques tel que module de Young statique E, le coefficient de Poisson statique v, le module d'incompressibilité K, le module de cisaillement G, et leurs combinaisons, comme par exemple le module œdométrique Kv, ou le module de Lamé λ. De manière préférée, deux paramètres élastiques statiques sont déterminés. Avantageusement, il est possible de décrire les propriétés élastiques du système cimentaire, déterminées de manière unique par la combinaison de deux paramètres élastiques statiques, en utilisant des formules bien connues en élasticité, telles que décrites dans le tableau 2 ci- dessous.
TABLEAU 2
Figure imgf000016_0001
Selon un mode de réalisation de l'invention, au moins un paramètre mécanique déterminé selon l'étape D, en particulier choisi parmi les paramètres de déformabilité statiques, de préférence les paramètres élastiques statiques, est un paramètre en conditions drainées. Dans un système drainé, il est fait l'hypothèse d'une circulation des fluides de pores, par exemple des liquides, entre le système cimentaire et son environnement (encaissant) afin de conserver la pression de pores constante. A l'opposé, dans un système non drainé, il est fait l'hypothèse d'aucun échange de fluides de pores entre le système cimentaire et son environnement. La valeur des paramètres mécaniques du système cimentaire, en particulier des paramètres de déformabilité statique, sera différente selon que le système est drainé ou non drainé.
Selon un autre mode de réalisation, le procédé comprend en outre la détermination d'au moins un paramètre de couplage hydro-mécanique, tel que le coefficient de Biot b ou le coefficient de Skempton B, et leurs combinaisons. La connaissance de tels paramètres de couplage hydro-mécanique permet notamment de déterminer les paramètres de déformabilité statiques, de préférence les paramètres élastiques statiques, en conditions non drainées à partir de la connaissance des paramètres élastiques statiques, en conditions drainées. Ces paramètres peuvent, par exemple, être utilisés pour simuler le comportement d'une gaine de ciment dont les pores sont saturés en eau lorsque la vitesse de chargement est nettement plus rapide que la vitesse de diffusion des pressions de pores. L'ensemble des formules (II) utilisées pour passer des propriétés élastiques statiques drainées (indice d) à celles statiques non-drainées (indice u), et réciproquement, est donné ci-dessous.
Figure imgf000017_0001
2,vd - {\ + 2vd )- bB
v„
- {\ - 2vd )- bB
Figure imgf000017_0002
3E,
3 - (l - 2vd )- bB (l - bB)- Ku
K - Kd
1 - bB
Le module de Young E et le coefficient de Poisson v sont des constantes élastiques qui sont classiquement déterminées lors d'un essai de compression uniaxiale ou triaxiale où un échantillon cylindrique d'élancement 2 est soumis à une augmentation de la contrainte axiale alors qu'aucun confinement n'est appliqué et que la pression de pores est constante. Ces coefficients sont alors définis classiquement, dans la zone où les déformations sont réversibles, par les formules suivantes : E = ΑσαχΜ I ΑεαχΜ , v = -Asradial I Αεαχία1 , avec Δσαχίαι l'incrément de contrainte axiale, Δεαχ l'incrément de déformation axiale et Asradiai l'incrément de déformation radiale.
Le module d'incompressibilité K est une constante élastique classiquement déterminée lors d'un essai de compression isotrope où un échantillon est soumis à une augmentation de contrainte axiale égale à l'augmentation de confinement (dite contrainte isotrope) alors que la pression de pores est constante. Ce coefficient est alors défini classiquement, dans la zone où les déformations sont réversibles, par la formule suivante : K = Δσ^ I Asvolumjque , avec
Δσί80ιΓορε l'incrément de contrainte isotrope et AsvoiUmique l'incrément de déformation volumique.
Le module de cisaillement G n'est généralement pas mesuré directement, bien qu'il puisse l'être théoriquement, et est évalué à partir de la connaissance de deux autres paramètres élastiques.
Le module œdométrique, Kv est une constante élastique classiquement déterminée lors d'un essai de compression œdométrique où un échantillon cylindrique est soumis à une augmentation de la contrainte axiale alors qu'aucun déplacement radial n'est autorisé et que la pression de pores est constante. Ce coefficient est alors défini classiquement, dans la zone où les déformations sont réversibles, par la formule suivante : Kv = ΔσαχΜ I Δεαχία1 , avec Δσαχ l'incrément de contrainte axiale et Δεαχ l'incrément de déformation axiale.
Le procédé selon l'invention permet de déterminer les paramètres de déformabilité statiques, qui ont pour avantage d'être ceux utilisés dans la modélisation des gaines ou des bouchons de ciment des puits de forage. Ces paramètres statiques ne doivent pas être confondus avec ceux dits « dynamiques », qui sont évalués à partir de la connaissance de la vitesse des ondes de compression et de cisaillement.
Les paramètres de couplage hydro-mécanique sont par exemple le coefficient de Biot
(b) et le coefficient de Skempton (B). Le coefficient de Biot (b) est le coefficient qui permet de définir la contrainte effective, c'est-à-dire la contrainte qui induit la déformation d'un matériau dans le cas où la pression de pores n'est pas nulle, et s'exprime dans la relation suivante : Δσ'= Δσ - b · P , où σ' est la contrainte effective, σ est la contrainte totale, et Pp est la pression de pores. Le coefficient de Skempton (B) est le coefficient qui permet de calculer la variation de pression due à une variation de contrainte isotrope sans expulsion de fluide pour un échantillon saturé, et s'exprime selon la relation APp = B■ Ασ , dans laquelle
Δσ est la variation de contrainte isotrope , et ΔΡΡ est la variation de pression de pores.
Selon le procédé de l'invention, ces paramètres peuvent être déterminés pour le système cimentaire dès son plus jeune âge, ce qui permet notamment une meilleure modélisation du comportement mécanique du système cimentaire. Il peut être important de connaître le comportement du ciment au stade précoce dans le cadre des opérations de forage, par exemple pour des opérations qui doivent pouvoir être effectuées rapidement après le pompage du laitier de ciment dans le puits, telles que des logs thermiques permettant de localiser le sommet de la colonne de ciment, des opérations de forage supplémentaires une fois le tubage posé, la perforation des zones de productions, ou encore pour évaluer l'état de contraintes dans la gaine de ciment une fois le ciment durci. Le procédé
L'organigramme de la figure 1 représente schématiquement le procédé selon l'invention, qui comprend les étapes A à D décrites ci-dessous.
Le système cimentaire testé par le procédé selon l'invention présente une composition initiale Co connue et une finesse Φ. Un autre paramètre initial utilisé par le procédé selon l'invention est la vitesse des ondes de compression dans le système cimentaire en fonction du temps Vp(t). La connaissance de ces paramètres initiaux permet, grâce au procédé selon l'invention, de déterminer au moins un paramètre mécanique du système cimentaire.
Le procédé selon l'invention ne se limite pas à la combinaison des étapes A à D, selon l'ordre cité, et des étapes intermédiaires peuvent être réalisées, sans pour autant sortir du cadre de l'invention.
Etape de mesure de la vitesse des ondes de_ compression en Jonction du temps du système cimentaire
Selon un mode de réalisation de l'invention, le procédé comprend en outre une étape initiale de mesure de la vitesse des ondes de compression en fonction du temps Vp(t) dans un échantillon du système cimentaire. Une telle mesure peut, par exemple, être réalisée selon une méthode de type Ultrasonic Cernent Analyser UCA, telle que par exemple décrite dans le brevet US 4,259,868. Les méthodes de type UCA sont des méthodes bien connues pour évaluer la qualité d'un système cimentaire, et font partie de la deuxième catégorie de méthodes de mesure décrite plus haut, permettant de caractériser le comportement mécanique d'un système cimentaire. Il s'agit de méthodes non destructives permettant d'analyser des échantillons de ciment, sous forme coulis de ciment, en utilisant la mesure de la vitesse d'ondes ultrasoniques, notamment de compression, traversant l'échantillon, sans nécessiter de ramener les échantillons dans les conditions environnementales (pression et température) du laboratoire préalablement à la mesure. Une amélioration des méthodes de type UCA permet en outre de mesurer les vitesses des ondes de cisaillement, et sont désignées en tant que mesures MPRO (Reddy et al.(7)). La vitesse des ondes de compression et de cisaillement en fonction du temps ainsi enregistrée permet de calculer d'une part la résistance à la compression uniaxiale, reliée à la vitesse des ondes de compression par une fonction de corrélation, et d'autre part les paramètres élastiques dynamiques, à partir des relations suivantes : Coefficient de poisson dynamique :
Figure imgf000019_0001
Module de Young dynamique : l - V yn ) La figure 3 illustre un appareil de mesure de type UCA, tel que décrit dans le brevet
US 4,259,868. Un tel dispositif de mesure comprend une cellule de mesure 34 destinée à recevoir le coulis de ciment et à être mise sous pression. Cette cellule 34 est mise sous pression dans un autoclave 33 par le biais de son insertion dans une ouverture 34a et d'une connexion avec un conduit de mise sous pression 34b. La pression et la température de l'autoclave, transmises à la cellule 34, sont contrôlées par des variateurs (37a, 36) et des écrans de contrôle (37, 35). Des transducteurs ultrasoniques (non représentés) sont couplés à la cellule de mesure 34, et permettent l'émission et la réception d'ondes ultrasoniques traversant l'échantillon du système cimentaire. Un ordinateur de contrôle 31 est relié à l'autoclave et comporte un clavier 38 permettant l'entrée de données et les fonctions de commande dans le dispositif de mesure, et un traceur numérique 32 permettant de visualiser la vitesse des ondes mesurées ainsi que la résistance à la compression uniaxiale calculée en temps réel.
Le protocole de mesure selon une méthode de type UCA ou MPRO suit le schéma suivant :
· préparation d'un volume de laitier de ciment ;
• remplissage de la cellule de mesure ;
• mise en conditions de l'échantillon du système cimentaire en appliquant une rampe d'augmentation de pression et de température pour atteindre les conditions in situ (Pi, ΤΊ). Une fois ces conditions atteintes, la pression et la température restent constantes ;
· mesure, en fonction du temps, de :
o la vitesse des ondes de compression, et de vitesse des ondes de cisaillement dans le cas d'une mesure de type MPRO ;
o de la pression ; et
o de la température ;
et utilisation d'une fonction de corrélation pour évaluer la résistance à la compression uniaxiale du système cimentaire en fonction du temps, et utilisation des relations (III) et (IV) décrites ci-dessus, pour évaluer le coefficient de poisson dynamique et le module de Young dynamique en fonction du temps dans le cas d'une mesure MPRO ;
· une fois l'essai terminé, démontage de l'éprouvette et observation de l'échantillon de ciment. Il est entendu que l'invention ne se limite pas à la mesure de Vp(t) selon les méthodes de type UCA et MPRO décrites ici, et que l'homme du métier pourra choisir toute méthode permettant la mesure de la vitesse des ondes de compression d'un système cimentaire. Des exemples de telles méthodes sont par exemple décrits dans les documents suivants : US585981 1 , US5763773, US5357481 , US5168470, US5001676, US4813028, US4779236, US4255798, US3401773, US25381 14. Il est également évident que le procédé selon l'invention peut être mis en œuvre sans avoir à effectuer cette étape de mesure de Vp(t), si cette donnée est connue par ailleurs pour le système cimentaire étudié.
Un avantage à utiliser la mesure de la vitesse des ondes de compression est que celle- ci constitue une mesure simple et normalisable, contribuant de ce fait à fournir un procédé facilement industrialisable. Le protocole expérimental utilisé pour déterminer au besoin Vp(t) est simple. Seule la connaissance de la vitesse des ondes de compression est nécessaire pour déterminer le comportement mécanique d'un ciment selon le procédé objet de la présente invention.
Etape A
L'étape A consiste à déterminer le degré d'hydratation du système cimentaire en fonction du temps a(t), à partir de la connaissance de la vitesse des ondes de compression Vp(t) du système cimentaire en question, mesurée dans un échantillon du système cimentaire, à une pression Pi et une température ΤΊ .
Pour cela, le procédé selon l'invention utilise une corrélation empirique 1 10 qui relie a(t) à Vp(t).
Selon un mode de réalisation de l'invention, le degré d'hydratation du système cimentaire en fonction du temps a(t) est calculé à partir de Vp(t) selon une relation linéaire. De manière préférée, le degré d'hydratation du système cimentaire en fonction du temps a(t) est calculé à partir de Vp(t) selon la relation linéaire suivante :
= (Vp - Vo ) v)
{y. - y. ) m avec V0 et Vp correspondant respectivement à la vitesse des ondes de compression mesurée dans l'échantillon du système cimentaire au temps t = 0 et au temps t, et V∞ correspondant à la vitesse des ondes de compression dans l'échantillon du système cimentaire complètement hydraté.
Les valeurs de Vp et Vo sont connues au début de l'étape A. La valeur de V∞ peut être connue expérimentalement avec la même mesure que celle permettant de connaître Vp(t), et correspond à un échantillon dans lequel l'hydratation est complète. Connaissant Vp(t), V∞ est par exemple calculé à l'aide d'une régression linéaire. Selon un mode de réalisation de l'invention, la valeur de V∞ est issue d'un modèle établi à partir de données expérimentales préliminaires issues d'essais de type UCA, permettant de prédire cette valeur pour tout système cimentaire. Selon ce mode de réalisation, il n'est pas nécessaire de réaliser une mesure de V∞ pour le système cimentaire objet du procédé selon l'invention.
Selon un mode de réalisation de l'invention, V∞ est égal à environ 3980 m/s.
Etape B
Le procédé selon l'invention comprend une étape B consistant à déterminer le degré d'hydratation a(t), issu de l'étape A, en fonction de valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn.
Cette détermination utilise un modèle cinétique 120 de l'hydratation du système cimentaire. Selon un mode préféré de l'invention, ce modèle cinétique 120 comprend deux stades correspondant à deux sous-étapes de l'étape B : B-i et B-ii.
Selon un mode de réalisation de l'invention, le processus d'hydratation comprend un premier stade où l'hydratation est gouvernée principalement par un processus de nucléation et de croissance, et un deuxième stade où l'hydratation est gouvernée principalement par un processus diffusion des ions, ledit deuxième stade commençant lorsque le degré d'hydratation a atteint une valeur de degré d'hydratation seuil a*, cette valeur seuil a* étant fonction de la température. L'étape B comprend les sous-étapes suivantes :
- (B-i) la détermination du degré d'hydratation a(t) pendant le premier stade du processus d'hydratation du système cimentaire;
- (B-ii) la détermination du degré d'hydratation a(t) pendant le deuxième stade du processus d'hydratation du système cimentaire.
Chacune des sous-étapes B-i et B-ii prend en compte la finesse Φ du système cimentaire, la pression et la température pour la détermination du degré d'hydratation a(t).
Le procédé prend donc avantageusement en compte les effets de la température et de la pression pour simuler l'évolution du degré d'hydratation du système cimentaire.
Selon un mode de réalisation de l'invention, la valeur de degré d'hydratation seuil a* est évaluée en minimisant l'écart entre a(t), déterminé en utilisant le modèle cinétique 120, et a(t), déterminé expérimentalement à partir de la vitesse des ondes de compression, pour différentes températures, et à une pression constante, de manière à prendre en compte la variation de a* en fonction de la température dans l'étape B. Avantageusement, la détermination expérimentale de a(t) est effectuée à différentes températures, et à une pression constante inférieure ou égale à 1 MPa, de préférence inférieure ou égale à 0.5 MPa, encore plus préférentiellement inférieure ou égale à 0.3 MPa. A de telles pressions, il peut être considéré que l'effet du volume d'activation est négligeable.
De manière préférée : - le système cimentaire de composition initiale Co comprend un ciment et de l'eau, le ciment comprenant au moins une phase initiale réactive X ;
- le degré d'hydratation a(t) déterminé à l'étape B correspond à la moyenne pondérée des degrés d'hydratation de chacune des phases initiales réactives X du ciment ;
- le degré d'hydratation de chacune des phases initiales réactives X du ciment est fonction du rapport entre l'affinité chimique Ax{a) de la phase initiale réactive X, cette affinité chimique Ax{a) contrôlant le taux d'évolution de l'hydratation de la phase initiale réactive X, et le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau τ x ;
- et le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau τχ est fonction de la finesse Φ du ciment, de la pression et de la température.
Avantageusement, le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau τχ s'exprime selon l'équation suivante :
Figure imgf000023_0001
avec Φο la finesse d'un ciment de référence, de préférence Φο = 3600 cm2/g, T0 la température initiale et P0 la pression initiales au temps t=0 du processus d'hydratation, R la constante de gaz parfaits, ΔΕΧ l'énergie d'activation, AVX le volume d'activation pour la nucléation et la croissance des hydrates lors du premier stade du processus d'hydratation, nx une constante, déterminée expérimentalement dont les valeurs sont publiées dans la littérature (Bernard et al.(12)).
Sous-étap_e_ B-i
Lors de la sous-étape B-i, on considère le premier stade d'évolution du degré d'hydratation, où l'hydratation est gouvernée principalement par un processus de nucléation et de croissance, tel que par exemple décrit par Avrami1^. Ce premier stade comprend les phases d'accélération et de ralentissement du processus d'hydratation dans le cas d'un ciment Portland. Cette sous-étape B-i consiste à déterminer le degré d'hydratation du système cimentaire à partir des degrés d'hydratation de chacune des phases initiales réactives X du ciment selon la relation suivante :
N
a = ∑ x'- x (VII)
x=l où N est le nombre de phases initiales réactives X. Par exemple, pour le clinker, X serait un des quatre constituants principaux du clinker ( C3S, C2S , C3A, C4AF). mx' est la fraction massique du constituant X. Le degré d'hydratation de la phase initiale réactive ax est défini comme le rapport entre la quantité de cette phase ayant réagie et la quantité initiale de cette phase, et vérifie l'équation suivante:
^ = W (VIII)
dî τχ
L'affinité chimique Ax (a ) de la phase initiale réactive X vérifie l'équation suivante : Ax (a ) = (1 - (a - a0x )) · [- ln (l - (a - a0x ))]' ^ (IX) où nx et ciox sont des constantes déterminées expérimentalement et connues par l'homme de l'art (Bernard et al 12-*).
Le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau τχ est fonction de la finesse Φ du ciment, de la pression et de la température. Au cours de la sous-étape B-i, τχ vérifie l'équation suivante :
Figure imgf000024_0001
Dans laquelle Φο est la finesse d'un ciment de référence, de préférence <ï>o=3600cm2/g, T0 la température initiale et P0 la pression initiales au temps t=0 du processus d'hydratation, R la constante de gaz parfaits, égale à 8,314 J°K-lmol-l , ΔΕΧ l'énergie d'activation, AVX le volume d'activation pour la nucléation et la croissance des hydrates lors des phases de croissance et de ralentissement du processus d'hydratation, égal à -27x 10"6 m3/mol, et nx est une constante déterminée expérimentalement et connue par l'homme de l'art (Bernard et al.(12)).
Les effets de la finesse du système cimentaire, de la pression et de la température sur la cinétique d'hydratation de ce dernier sont intégrés dans cette équation (VI). Cette équation établit que la vitesse d'évolution d'hydratation d'un système cimentaire croît avec la finesse, la pression et la température.
Sous-étap_e_ Bzii
Lors de la sous-étape B-ii, on considère le deuxième stade d'évolution du degré d'hydratation, où l'hydratation est gouvernée principalement par un processus de diffusion des ions. Ce deuxième stade comprend la phase de cure dans le cas d'un ciment Portland.
On appelle r* la valeur du rayon du grain de ciment entouré d'une couche d'hydrates lorsque le degré d'hydratation du ciment atteint une valeur de seuil a*. La valeur du rayon diminue avec l'avancement de l'hydratation, et son taux de diminution peut être écrit en fonction d'une constante de diffusion D (cm2/h) selon l'équation suivante : dr _ D
dt r * -r (χ)
Comme précédemment, on peut exprimer le degré d'hydratation a(t) en fonction de l'affinité chimique A(a). a(t) et A(a) vérifient l'équation suivante : da _ A{a)
(XI)
dt τ
L'affinité chimique vérifie l'équation suivante :
Figure imgf000025_0001
Comme précédemment, on peut exprimer le temps caractéristique τχ associé à la réaction de la phase initiale réactive X avec l'eau en fonction de la finesse Φ du ciment, de la pression et de la température. τχ vérifie l'équation (VI) précédemment explicitée.
Au cours de ce deuxième stade, le temps caractéristique de référence τχ {Τ0, Φ0 ) est contrôlé par la diffusion des ions à travers la couche des hydrates autour des grains de ciment (Bernard et al.(12)) : τ, (Γ0> Φο ) = ^ (ΧΙΠ)
Le modèle cinétique 120 utilisé dans le procédé selon l'invention prend avantageusement en compte l'effet de la finesse du système cimentaire, de la température et de la pression pour simuler les résultats de l'évolution du degré d'hydratation du système cimentaire. Le volume d'activation AVX est avantageusement considéré comme étant identique pour chacune des phases initiales réactives. La valeur de degré d'hydratation seuil a* est évaluée en minimisant l'écart entre a(t), déterminé selon l'étape A, et a(t), déterminé expérimentalement pour le système cimentaire, pour différentes températures, et à une pression constante. Avantageusement, la pression est inférieure ou égale à 1 MPa, de préférence inférieure ou égale à 0.5 MPa, encore plus préférentiellement inférieure ou égale à 0.3 MPa. A cette faible pression, il est en effet possible de négliger l'effet du volume d'activation.
Selon un mode de réalisation de l'invention, la pression n'affecte pas le degré d'hydratation seuil a*, qui varie seulement en fonction de la température. En particulier, au- delà d'une certaine valeur de la température d'hydratation, la phase de diffusion commence plus tôt à une température d'hydratation plus élevée, et la température plus élevée provoque l'augmentation du taux de C-S-H HD à long terme. Etage
Cette étape consiste à déterminer la composition du système cimentaire en fonction du temps C(t) et en fonction de valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn, et/ou de température Tn, à partir du degré d'hydratation a(t) déterminé à l'étape B.
Pour cela, le procédé selon l'invention utilise un modèle compositionnel 130, qui permet de connaître, à tout degré d'hydratation du système cimentaire, la teneur du système en chaque phase considérée. En particulier, la fraction volumique, massique ou molaire, de préférence volumique, de chaque phase du système cimentaire est déterminée.
En premier lieu, la connaissance de a(t), pour des valeurs souhaitées de finesse du système cimentaire, de pression et/ou de température permet d'estimer la fraction, de préférence volumique, des phases initiales du système cimentaire (phases réactives initiales du ciment et eau initiale consommée), en fonction du temps. Dans un second temps, à partir des relations entre la proportion molaire des réactants et des hydrates des réactions chimiques de l'hydratation du système cimentaire, la fraction, par exemple volumique, des phases hydratées du système cimentaire est calculée à un temps donné. En effet, les équations chimiques permettent de déterminer la proportion en nombres de moles des réactants et des hydrates. La fraction volumique est alors déterminée en connaissant la densité et la masse molaire des hydrates. Ainsi, les fractions, de préférence volumiques, des phases hydratées du système cimentaire, en particulier les phases C-S-H et de CH, sont déterminées. Enfin, la fraction, de préférence volumique, des aluminates est déterminée comme étant la différence entre la proportion initiale, de préférence le volume total initial, et la proportion calculée, de préférence le volume calculé, de toutes les autres phases.
Avantageusement, le modèle compositionnel 130 prend en compte la température à laquelle s'effectue l'hydratation, notamment en ce qui concerne la détermination des fractions, de préférence volumiques, des phases hydratées, en particulier de la phase C-S-H.
Lors de cette étape, la composition du système cimentaire C(t) est déterminée en fonction du temps, pour des valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn.
Dans le procédé selon l'invention, le système cimentaire comprend m phases, parmi lesquelles :
- au moins une phase initiale réactive X, de préférence choisie dans le groupe constitué par le silicate tricalcique C3S, le silicate dicalcique C2S, Paluminate tricalcique C3A, l'alumino ferrite tétracalcique C4F, et leurs associations ;
- au moins une phase hydratée Y issue de l'hydratation d'au moins une phase initiale réactive X, de préférence choisie dans le groupe constitué par le silicate de calcium hydraté C-S-H, l'hydroxyde de calcium CH, le trisulfoaluminate de calcium hydraté TSA, le monosulfate de calcium hydraté AFm, l'alumino ferrite de calcium hydraté, et leurs associations ; - de l'eau;
la composition initiale du système cimentaire Co comprenant un volume initial d'eau déterminé V° et au moins une phase initiale réactive X.
L'étape C comprend l'estimation de la composition du système cimentaire en fonction du temps C(t) par la détermination de la fraction molaire, massique ou volumique des m phases du système cimentaire.
Selon un mode de réalisation de l'invention :
le volume d'au moins une phase initiale réactive X en fonction du temps Vx (t) est calculé selon l'équation (XIV) :
Vx (t) = Vx° (l - ax (t)) (XIV) avec Vx le volume initial de la phase initiale réactive X, et ax (i) le degré d'hydratation de la phase initiale réactive X en fonction du temps ;
le volume d'au moins une phase hydratée en fonction du temps est calculé selon l'équation (XV) : ( =∑*7 - « ( w avec V* le volume occupé par la phase hydratée Y formée par la phase réactive X et ax (t) le degré d'hydratation de la phase réactive X en fonction du temps ;
le volume de l'eau en fonction du temps Vw (t) est calculé selon l'équation (XVI) : (ή = ν -∑ν - αχ (ή (XVI) avec Vl le volume initial d'eau dans le système cimentaire, V^ le volume de l'eau consommée par la phase X, ax (t) le degré d'hydratation de X en fonction du temps.
Le procédé selon l'invention peut en outre prendre en compte d'autres phases, telle que de l'eau entrante, diffusant du milieu environnemental vers les pores du système cimentaire, ou au moins une phase gazeuse. Dans ce cas, le procédé permet en outre de déterminer une fraction, par exemple volumique, d'eau entrante dans le système cimentaire en fonction du temps, et une fraction, par exemple volumique, de gaz en fonction du temps. Pour cela, les fractions de l'eau entrante et/ou des phases gazeuses sont déterminées en se basant sur le phénomène de retrait chimique se produisant lors de l'hydratation du système cimentaire.
Selon un mode de réalisation de l'invention, le volume de la phase hydratée C-S-H est déterminé en fonction de la température à laquelle se produit l'hydratation. En particulier, la fraction des phases C-S-H LD et C-S-H HD est déterminée différemment selon que l'hydratation s'effectue à température ambiante (20°C±5°C), ou à température supérieure ou inférieure à la température ambiante.
Dans le cas où l'hydratation se produit à température ambiante, les fractions, de préférence volumiques, de C-S-H LD et C-S-H HD peuvent s'exprimer selon les équations suivantes (Bernard et al. ^12-*) :
VLD{t)
Figure imgf000028_0001
(XVII)
VHD{t) =∑V _s_H (a'(t)- a *)x K[a'(t)- a *\ (XVIII) avec H la fonction de Heaviside, et i désigne chacune des réactions ayant la phase hydratée comme produit de la réaction.
La densité de C-S-H HD et C-S-H LD est exprimée par:
PHD = (! " <t>HD ) · Pg + <t>HD■ Pw (XIX)
PLD = (l " Φΐθ ) - Pg + <t>LD■ Pw (XX) avec pg la densité du globule, de préférence égale à 2,65 g/cm3, pw la densité de l'eau, égale à 1 g/cm3), et et LD les porosités respectives de C-S-H HD et C-S-H LD.
Dans le cas où l'hydratation se produit à une température supérieure à la température ambiante, la détermination de la composition C(t) prend en compte l'influence de la température, via le degré d'hydratation seuil a* dépendant de la température, sur les fractions de C-S-H LD et C-S-H HD, et éventuellement sur la fraction de l'eau entrante. Les fractions, de préférence volumiques, de C-S-H LD et C-S-H HD peuvent s'exprimer par les mêmes équations (XVII) et (XVIII) ci-dessus, qui prennent en compte une valeur de la porosité de la phase C-S-H HD différente de celle à température ambiante. La détermination de la porosité du C-S-H HD à une température donnée est obtenue en supposant que 1) la quantité de globules (C-S-H LD et de C-S-H HD) ne dépend que du degré d'hydratation, 2) le volume des globules est constant pour une température d'hydratation donnée, 3) la densité des globules de C-S-H n'est quasiment pas modifiée par la température, 4) la porosité totale n'est pas modifiée par la température et 5) le volume total des pores mesuré par passage à l'étuve à 105°C est considéré comme étant la somme des volumes des pores provenant des globules, du C-S-H LD, du C-S-H HD et du volume des capillaires.
En effet, dans le cas d'une température supérieure à la température ambiante, la fraction volumique de C-S-H HD et sa densité sont plus importantes (porosité de C-S-H HD plus faible) et la fraction volumique de C-S-H LD est plus faible qu'à température ambiante. La macroporosité du système cimentaire augmente également avec une température plus élevée, pour un même degré d'hydratation.
On obtient ainsi un système d'équation (XXI) permettant de calculer les fractions volumiques en C-S-H LD et C-S-H HD ainsi que la porosité du C-S-H HD à une température T2 en fonction de ces mêmes paramètres à une température ΤΊ, sous couvert que la différence de porosité capillaire soit connue.
Figure imgf000029_0001
La différence de porosité capillaire est par exemple prise égale à 6% entre le cas à température de 60°C et celui à température de 20°C, basé sur des mesures au porosimètre à mercure. Ainsi, les porosités de C-S-H HD et de C-S-H LD qui sont par exemple respectivement de 24% et 37% à une température ambiante de 20°C, sont respectivement de 17%) et 37%o à une température de à 60°C.
Etape D
Cette étape consiste à déterminer au moins un paramètre mécanique du système cimentaire en fonction du temps, et en fonction des valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn, à partir de la composition du système cimentaire C(t) déterminée à l'étape C.
Pour cela, le procédé selon l'invention utilise une méthode d'analyse multi- compositionnelle 140, basée sur des techniques d'homogénéisation telles qu'utilisées en micromécanique. Cette méthode d'analyse multi-compositionnelle 140 prend en compte un modèle multi-échelle du système cimentaire pour déterminer, à partir de techniques d'homogénéisation, au moins un paramètre mécanique du système cimentaire en fonction du temps, pour une valeur donnée de finesse du système cimentaire, de pression et/ou de température, connaissant la composition du système cimentaire C(t) déterminée à l'étape C, et connaissant l'évolution des composants du ciment du système cimentaire obtenue à partir du modèle de la cinétique d'hydratation utilisé à l'étape A. Le modèle multi-échelle comprend au moins une échelle élémentaire représentant des globules de la phase solide C-S-H et une échelle macroscopique du système cimentaire, de préférence un modèle multi-échelle comprenant trois échelles.
Lors de cette étape D, l'hétérogénéité physico-chimique du système cimentaire, en particulier du ciment durci, est représentée par un modèle multi-échelle. Une échelle locale comprend des unités élémentaires appelées globules qui représentant la phase hydratée solide C-S-H. Cette échelle locale représente la plus petite échelle du modèle multi-échelle utilisé lors de cette étape. Une échelle macroscopique représente le système cimentaire à l'échelle macroscopique. Des échelles intermédiaires peuvent être utilisées. Le module d'incompressibilité ks et le module de cisaillement gs des globules constituent les deux inconnues du modèle multi-échelle.
Classiquement, la micromécanique, en mécanique des matériaux, consiste à estimer les propriétés mécaniques d'un matériau en considérant un problème aux limites d'un volume élémentaire représentatif (VER). Selon
Figure imgf000030_0001
il y a trois étapes pour l'homogénéisation : description, localisation et homogénéisation.
Dans la première étape, le VER doit être bien choisi pour représenter le comportement macroscopique du matériau. Ainsi, pour incorporer toute l'information mécanique et géométrique des phases hétérogènes du matériau, le VER est de préférence suffisamment grand. De plus, dans le contexte des milieux continus, ce volume est de préférence suffisamment petit pour décrire la continuité de la structure macroscopique.
Dans la deuxième étape, l'interaction des phases dans le VER est considérée. Différents schémas d'homogénéisation bien connus sont utilisés. Le schéma « dilué » ne considère aucune interaction entre les inclusions dans le VER, et ne peut être utilisé que pour des fractions volumiques faibles des inclusions. Le schéma « Mori-Tanaka », qui prend en compte l'interaction entre les inclusions, est couramment utilisé pour une microstructure de composites à matrice et inclusions. La différence entre le schéma « Mori-Tanaka » et le schéma « auto-cohérent » est que ce dernier prend le matériau homogénéisé comme la matrice en résolvant un système des équations non-linéaires pour l'estimation des propriétés mécaniques effectives du matériau. Pour le cas d'un matériau composite à matrice et inclusions, le schéma plus avancé de « Ponte-Castaneda et Willis » (Pont Castenada et al 16-*) prend en outre en compte l'information de l'orientation et de la distribution spatiale des inclusions dans le VER. Après avoir choisi le schéma d'estimation, les relations entre les tenseurs locaux de déformations et de contraintes et les tenseurs macroscopiques de déformations et de contraintes sont établies dans cette étape. Ces relations forment un tenseur appelé « tenseur de localisation des déformations ». Dans les cas où les conditions aux limites sont homogènes, le lemme de Hilr07^ permet de démontrer que la déformation et la contrainte macroscopiques sont respectivement égales à la moyenne du champ des déformations et des contraintes microscopiques respectivement. Pour des inclusions ellipsoïdales élastiques linéaires homogènes plongées dans un milieu solide élastique linéaire homogène, Eshelby ' a démontré que la déformation dans les inclusions est homogène. Enfin, le tenseur de localisation des déformations est déterminé en fonction des caractéristiques de la matrice et des inclusions.
Dans l'étape d'homogénéisation, en combinant des équations du comportement du matériau au niveau local et macroscopique, les paramètres élastiques et de couplage hydromécanique peuvent être exprimés en fonction des tenseurs de localisation de déformations. Pour un cas drainé, le tenseur de la rigidité du matériau est exprimé par :
Chom = 3K omJ + 2GhomK (XXII) où im = LJ ; K = I - J ; 1= ik β a }k . δβ dénote le delta Kronecker, ^om est le
3 2
module d'incompressibilité drainé, Ghom est le module de cisaillement. Selon Zaoui(14), pour un cas isotrope avec n inclusions sphériques, le module d'incompressibilité drainé et le module de cisaillement effectifs sont calculés par les expressions (XXIII) et (XXIV) suivantes :
K =∑frkrAr v (XXIII) Ghom =∑frgrAr J (XXIV)
dans lesquelles fr, kr, et gr sont respectivement la fraction volumique, le module d'incompressibilité et le module de cisaillement de la phase r. Pour le cas où toutes les inclusions sont sphériques, le tenseur de localisation de déformation volumique A[ et le tenseur de localisation de déformation déviatorique Ad r sont données par les équations (XXV) et (XXVI) suivantes :
Figure imgf000031_0001
avec aa =—— , Ba = ^fr + , ko le module d'incompressibilité de la
3k0 + 4g0 o 5(3£0 + 4g0 ) '
matrice de référence, et go le module de cisaillement de la matrice de référence. Pour calculer le tenseur de rigidité de la phase solide, le tenseur de coefficient de Biot et le module Biot homogénéisés, deux cas sont considérés :
Dans le premier cas, le matériau considéré comprend une phase poreuse unique (les pores) et n-1 phase solide. D'après Ulm et al.(19), le tenseur de coefficient de Biot d'un matériau ayant une seule phase poreuse est donné par l'équation (XXVII) suivante:
Figure imgf000032_0001
Le signe (A) dénote la valeur moyenne volumique de A dans le volume V. Le module de Biot homogénéisé est écrit sous la forme (XXVIII) suivante : frc;1 : (l - 1 :(Α) ) (XXVIII)
Figure imgf000032_0002
où cr est le tenseur de rigidité de la phase r qui est définit par la formule (XXIX) suivante : cr = 3krJ + 2grK (XXIX)
Dans le deuxième cas, le matériau est un matériau poreux multi-échelle, par exemple deux échelles par souci de simplification. La généralisation à N échelles successives se fait alors sans difficulté. La configuration suivante est prise en compte : le matériau comprend q phases poreuses (solide et les petits pores), un volume de pore (larges pores) et n-q-1 phases solides. L'homogénéisation consiste à réaliser deux étapes I et II : la première étape II consiste à homogénéiser séparément q phases poreuses pour obtenir les paramètres effectifs de chaque phase de l'étape I ( cr ,
Figure imgf000032_0003
qui sont respectivement le tenseur de rigidité, le tenseur de coefficient de Biot, le module de Biot, et le tenseur de rigidité de la phase solide de la phase r dans l'étape I) comme présenté dans le premier cas ci-dessus. La deuxième étape II consiste à homogénéiser q phases poreuses, un volume des larges pores et n-q-1 phases solides. Le tenseur Biot homogénéisé dans la deuxième étape est exprimé par l'équation (XXX) suivante (Ulm et al.(19)) :
Figure imgf000032_0004
où br est le coefficient de Biot de la phase r. Il faut noter que le coefficient de Biot des phases solides est égal à 0. Le module de Biot effectif est donné par l'expression (XXXI) suivante (Ulm et al 19-*) : -hom - 1 {A)VR ): (l - ftr')+ } (XXXI)
Figure imgf000033_0001
Selon un mode de réalisation, l'hétérogénéité du système cimentaire à l'étape D se manifeste à trois échelles (Jennings ^11-*) :
- Une échelle élémentaire Ech.O dans laquelle les globules ont une longueur caractéristique de l'ordre de 10~9m, et sont considérés comme représentant la phase solide de C-S-H.
- Une première échelle Ech.l correspond aux phases C-S-H LD et C-S-H HD, sous forme de globules, ayant une longueur caractéristique comprise entre environ 10~9 - 10~8m, et comprend en outre les pores de gel C-S-H.
- Une deuxième échelle Ech.2 correspond au système cimentaire ayant une longueur caractéristique supérieure à environ 10"8m, et comprend les phases C-S-H LD, C-S-H HD, CH, le volume des pores capillaires, les aluminates et les phases initiales réactives (les grains du clinker).
Au sein de la première échelle Ech.l, deux inclusions sont considérées : d'une part la matrice formée par des globules solides, et d'autre part les pores dans les phases C-S-H LD et C-S-H HD. La matrice occupe par exemple 63% de volume pour le cas de la phase C-S-H LD, et elle occupe par exemple une proportion supérieure à 76%, dépendant de la température, de volume pour la phase C-S-H HD. Le schéma de Mori-Tanaka (Ghabezloo^13^) est choisi, avec la phase solide jouant le rôle le milieu de référence. Ce schéma de Mori- Tanaka est approprié pour les matériaux dont la phase solide est dominante et pour lesquels il y a des interactions entre les particules (Bernard et a^l2)). Les paramètres mécaniques effectifs, par exemple le module d'incompressibilité et le module de cisaillement, de la phase C-S-H, sous ses deux formes C-S-H LD et C-S-H HD, peuvent être déterminées selon des équations (XXXII) et (XXXIII) suivantes, présentées dans les travaux de Ghabezloo(13' 20
= (i - ) I (XXXII)
G = (I - )¾ pooan) Où X désigne la phase C-S-H LD ou la phase C-S-H HD, ks et gs sont respectivement le module d'incompressibilité et le module de cisaillement des globules qui représentent la phase hydratée solide C-S-H à l'échelle élémentaire Ech.O (« s » comme solide).
En supposant une géométrie sphérique pour toutes les phases, les tenseurs de localisation de déformations sont calculés selon les équations (XXV) et (XXVI) décrites ci- dessus, avec ko=ks et go=gs (ko est le module d'incompressibilité de la matrice de référence, et go est le module de cisaillement de la matrice de référence). Au sein de la deuxième échelle Ech. 2, l'évolution de la fraction volumique des phases du système cimentaire au cours de l'hydratation est prise en compte. Avantageusement, le schéma « auto-cohérent » est utilisé à cette échelle. Celui-ci permet de prendre en compte le seuil de percolation, et est approprié pour estimer les propriétés élastiques du système cimentaire à cette échelle (Ghabezloo1'13^). Cette deuxième échelle Ech.2 comprend six phases : C-S-H LD (LD), C-S-H HD (HD), CH, le volume des pores capillaires (cp), des aluminates (AL), et les phases initiales réactives du ciment, qui sont par exemple les grains de clinker (CK). A cette échelle Ech.2, la phase des aluminates (phase solide) est supposée avoir des propriétés élastiques similaires à celles de C-S-H LD, qui ne sont pas modifiées par la température d'hydratation. Les expressions des modules élastiques homogénéisés sont les suivantes :
KCP = (flD + /AL )KLD ALD,cp + fHDKHD AHD CP + fCH ^CH ACH,cp + f CK ^ CK ACK ,cp (XXXIV) = (fLD + /AL ^LD 1 ALD,cp + fHD GHD AHD,cp + fcH g CH ACH ,cp + fcK g CK ACK ,cp (XXXV)
Les tenseurs de localisation sont déterminés selon les équations (XXV) et (XXVI) décrites ci-dessus, avec k0 = K^°M , et g0 = G^°M (ko est le module d'incompressibilité de la matrice de référence, et go est le module de cisaillement de la matrice de référence).
Les paramètres à cette échelle Ech.2 sont par exemple les propriétés élastiques de CH et des phases initiales réactives du ciment (par exemple les quatre composants principaux du clinker).
Programme informatique Selon un autre de ses aspects, l'invention concerne un produit programme informatique permettant de conserver dans une mémoire d'une unité d'un processeur ou sur un support de mémoire amovible approprié pour coopérer avec ladite unité du processeur, le produit- programme comprenant des instructions pour mettre en œuvre le procédé selon l'invention
EXEMPLES
L'invention sera mise en évidence avec ses avantages apparents au travers de l'exemple ci-dessous, donné à titre non limitatif.
Le système cimentaire ceml
Le système cimentaire ceml donné en exemple ci-dessous comprend un ciment Portland de classe G, dont la composition est donnée dans le tableau 3. TABLEAU 3
Figure imgf000035_0001
La composition initiale Co du système cimentaire ceml est donnée dans le tableau 4 -dessous. Elle correspond à un rapport eau sur ciment égal à environ 044 (w/c en anglais).
TABLEAU 4
Figure imgf000035_0002
La fabrication du laitier de ciment est effectuée de la manière suivante, par un mixage en 5 étapes:
· mélange de l'eau distillée avec un agent anti-sédimentation (D153);
• mixage pendant au moins 5 minutes avec une vitesse d'environ 4200 tours par minute ;
introduction des deux autres adjuvants : un dispersant et un agent anti-mousse (D604AM et D47) ;
· malaxage à haute vitesse (une vitesse de 4200 tours par minute) pendant environ 15 secondes en ajoutant du ciment ;
• augmentation de la vitesse à 12000 tours par minute en mixant pendant 35 secondes. Test UCA sur un échantillon du système cimentaire ceml
Le principe de la méthode et un dispositif de type UCA ont été décrits plus haut et s'appliquent ici. Après l'installation de la cellule UCA dans le système de génération de pression et de température, la pression est tout d'abord augmentée à la valeur souhaitée (par exemple 0,3 MPa, 20 MPa ou 40 MPa) et la température est ensuite augmentée jusqu'à température de l'essai en 30 mn. Les mesures de temps de transit ont été réalisées dès le début de l'essai, à t=0 et T=25°C.
La figure 4 montre la vitesse des ondes de compression en fonction du temps Vp(t) obtenue.
Essai calorimétrique sur un échantillon du système cimentaire ceml
Le principe d'un essai calorimétrique a été décrit plus haut et s'applique ici. Le système cimentaire ceml est testé en réalisant un essai calorimétrique isotherme à 25°C.
La figure 4 montre le degré d'hydratation issu de cet essai calorimétrique. Les courbes a(t) et Vp(t) se confondent presque complètement, ce qui montre qu'une relation linéaire existe entre l'évolution au cours du temps du degré d'hydratation a(t) du système cimentaire, mesuré par calorimétrie, et l'évolution au cours du temps de la vitesse des ondes de compression Vp(t).
La figure 5 montre que la relation entre a(t) et Vp(t) est quasi linéaire. La régression linéaire s'exprime ainsi : V = 2500 a + 1479.
Essais uniaxiaux sur un échantillon du système cimentaire ceml
Des essais uniaxiaux permettant d'obtenir des mesures des paramètres mécaniques, en particulier de déformabilité, sont effectués sur le système ceml. Afin de mesurer l'effet de la température sur l'évolution des propriétés élastiques du système cimentaire en cours d'hydratation, les blocs cylindriques d'une longueur de 250 mm et d'un diamètre de 100 mm ont été curés à 20°C et 60°C à la pression atmosphérique dans l'eau saturée en chaux. Après deux jours de cure, ces blocs ont été carottés et sciés pour obtenir des échantillons cylindriques d'une longueur d'environ 100 mm et d'un diamètre d'environ 40 mm. Ces échantillons ont ensuite été conservés dans une solution neutre (pH=13). Des essais uniaxiaux ont été effectués aux âges suivants : 3 jours, 4 jours, 7 jours, 14 jours et 35 jours. Les essais uniaxiaux ont été effectués à l'aide d'une presse d'une capacité de 50 tonnes. La vitesse de déplacement uniaxial utilisée est de 0,3 μητ/s. Trois capteurs LVDT (Linear variable Differential Transformer) ont été employés pour mesurer la déformation axiale. Le module de Young est évalué à partir de cycles de décharge - recharge de OMPa à 18MPa de contrainte axiale.
Etapes A et B du procédé : détermination de ait) à Pu Tu et à différentes valeurs de finesse Φ„ du système cimentaire, de pression P„ et de température T„
La vitesse des ondes de compression en fonction du temps Vp(t) est mesurée dans un échantillon du système cimentaire ceml, à une pression Pi et une température Tl s tel que décrit dans le test UCA ci-dessus.
V∞ est déterminé à partir du test UCA, et correspond à 37979 m/s. Connaissant Vp(t) (voir figure 4), le degré d'hydratation en fonction du temps a(t) est déterminé en utilisant la relation linéaire CL = Vp ~ VQ )/ V∞— V0 ) .
Ce degré d'hydratation en fonction du temps a(t) est déterminé pour une pression Pl s et une température Tls avec Pi = 0.3 ou 20 MPa et Ti = 25 °C.
La valeur άετχ (Τ0 , Φ0 ) utilisée dans la formule (VI) dans le cas du premier stade du processus d'hydratation, gouverné principalement par les phénomènes de nucléation et de croissance, est calculée, en utilisant les valeurs des constantes connues par l'homme de l'art (Bernard et al.(12)). Ces valeurs sont reprises dans le tableau 5. Pour chaque phase initiale réactive (C3S, C2S, C3A, C4AF), l'affinité chimique est calculée en utilisant les valeurs du tableau 5.
T ABLEAl J 5
Figure imgf000037_0001
La valeur άετχ (Γ0 ,Φ0) utilisée dans la formule (VI) dans le cas du deuxième stade du processus d'hydratation du système cimentaire, gouverné principalement par le phénomène de diffusion (période de cure dans le cas de ceml), est calculée à partir des valeurs données dans le tableau 6 suivant :
TABLEAU 6
Figure imgf000037_0003
La valeur du degré d'hydratation seuil a* pour le système cimentaire ceml est donnée dans le tableau 7 ci-dessous, a* est évalué en minimisant l'écart entre les courbes déterminées à partir du procédé selon l'invention, et les courbes expérimentales pour différentes températures sous 0,3MPa car, pour cette faible pression, on peut négliger l'effet de la variable ΔΥ exprimée dans la formule (VI).
TABLEAU 7
Figure imgf000037_0002
Les figures 6 à 8, qui comparent des données expérimentales issues d'essais UCA tels que décrit plus haut, et les données issues de la détermination du degré d'hydratation selon les étapes A et B du procédé objet de l'invention (courbes en pointillés, indiquée en tant que « simulation » dans la légende), montrent la bonne capacité prédictive des étapes A et B du procédé objet de l'invention.
Etape C du procédé : détermination de la composition C(t))
Dans cette étape, les stades 1 et 2 du processus d'hydratation, comprenant les phases d'hydratation initiale et dormante, ne sont pas pris en compte. Les stades 3 et 4, comprenant les phases d'accélération et de ralentissement de l'hydratation, sont pris en compte et simulées selon par la théorie de la nucléation et croissance d'Avrami1^. Le stade 5 comprenant la phase de cure est également pris en compte, et est considéré comme étant contrôlé par des phénomènes de diffusion.
La figure 9 représente l'évolution des fractions volumiques des différentes phases du système cimentaire ceml en fonction du degré d'hydratation, pour une température d'hydratation de 20°C et à pression atmosphérique.
Etape D du procédé
Le tableau 8 présente les modules d'incompressibilité k et module de cisaillement g des différents composants solides pris en compte lors de l'analyse multi-échelle, selon l'étape D du procédé : les quatre phases du clinker du système ceml, la phase hydratée CH et les globules représentant la phase solide C-S-H du système cimentaire ceml. Toutes ces valeurs ont été obtenues à partir de la littérature, sauf en ce qui concerne les modules des globules qui ont été obtenus par calage avec les données d'essais réalisés par Bourissai(21).
L'analyse multi-échelle inclut deux homogénéisations :
- la première homogénéisation considère d'une part la matrice formée par les globules de C-S-H solides, et d'autre part les pores dans les phases C-S-H LD et C-S-H HD. Elle permet de calculer les modules d'incompressibilité et de cisaillement des C-S-H LD et C-S-H HD « poreux ». La longueur caractéristique est d'environ 10"9 - 10"8m. La matrice occupe 63% de volume pour le cas de la phase C-S-H LD, et elle occupe une proportion supérieure à 76%o, dépendant de la température, de volume pour la phase C-S-H HD ;
- la seconde homogénéisation considère d'une part les phases solides (C-S-H LD, C- S-H HD, CH, aluminates et phases initiales réactives du clinker) et d'autre part la porosité capillaire. Elle permet de calculer les modules d'incompressibilité et de cisaillement du système cimentaire et par la suite l'ensemble des paramètres élastiques et hydro -mécaniques. La longueur caractéristique est d'environ 10"8 m. TABLEAU 8
Figure imgf000039_0001
Les figures 10 à 13 montrent une bonne correspondance entre les paramètres mécaniques déterminés selon le procédé de l'invention, et les données expérimentales obtenues dans la littérature (Boumiz et al 1-* ; Haecker et al 2-*), validant ainsi le procédé selon l'invention. Les données issues de la littérature ont été corrigées, afin de transformer les modules dynamiques en modules statiques, à l'aide de la formule suivante : statique dynamique
Sur la figure 10, par exemple, sont comparées les données mesurées (« E-Mesures » dans la légende, symboles en forme de croix) et déterminées par le procédé selon l'invention du module de Young (E), pour une température de cure de 60°C (« Ehom » dans la légende). Les modules d'incompressibilité (K) et de cisaillement (G) déterminés par le procédé selon l'invention sont également représentés (« Khom » et « Ghom » dans la légende). Le module de cisaillement de la pâte de ciment hydratée à 60°C, à l'âge de lOOOh (42 jours) est de 6,lGPa. La différence entre cette valeur et la valeur mesurée (5,9GPa) est de 3,3%.
Sur la figure 11 , on observe également une bonne correspondance entre les données expérimentales (représentées pas des symboles en croix et en losange, et indiquées en tant que « E-mesures » dans la légende) et les paramètres mécaniques déterminés selon l'invention (courbes en trait plein et en pointillé, et indiqués en tant que « Ehom » dans la légende). Cette figure illustre en outre l'effet de la température d'hydratation sur l'évolution des propriétés élastiques du système cimentaire : la valeur du module de Young à long terme est plus grande pour l'hydratation à 20°C que pour celle à 60°C.
Sur la figure 12, on observe une très bonne correspondance entre les données mesurées
(représentées pas des symboles en carré et en triangle, et indiquées en tant que « Ed/Gd- mesures » dans la légende) et prédites selon l'invention (courbes en trait plein et en pointillé, et indiqués en tant que « Ed/Gd-calcul » dans la légende), tant pour le module de Young que pour celui de cisaillement.
Sur la figure 13, on observe également une bonne correspondance entre les données expérimentales (représentées par les différents symboles) et les paramètres mécaniques déterminés selon l'invention (courbes en trait plein). Cette figure illustre en outre l'effet du rapport eau sur ciment sur l'évolution des propriétés élastiques du système cimentaire : la valeur du module de Young à long terme est plus grande pour un rapport eau sur ciment faible. Liste de références bibliographiques citées dans la description :
(1) Boumiz A., Veet C, Cohen Tenoudjit F., "Mechanical properties of cernent pastes and mortars at early âges", Advanced Cernent based materials, 3, p.94-106 (1996).
(2) Haecker C.J., Garboczi E.J., Bohn J.W., Sun Z., Shah S.P., Voigt T., "Modeling the linear elastic properties of Portland cernent paste", Cernent concrète research, 35, p. 1948-1960 (2005).
(3) Avrami, M., "Kinetics of phase change, Journal of Chemical Physics", 7, p. 1103- 1124, 9, p.177-184 (1939-1940).
(4) Kondo R., Kodama, M., "On the hydration kinetics of cernent, Semento Gijutsu Nenpo", 21, p. 77-828 (1967).
(5) Fuji, K., Kondo, W., "Kinetics of the hydration of tricalcium silicate", Journal of the American Ceramic Society, 57 (11), p. 492-497 (1974).
(6) Jennings H. M., Tennis P.D., "Model for the developing microstructure in Portland
Cernent Pastes", J. Am. Ceram. Soc, 77 (12), p. 3161-3172 (1994).
(7) Reddy, B.R., Santra, A., McMechan, D., Gray, D., Brenneis, C, Dunn, R., "Cernent mechanical property measurements under wellbore conditions", SPE 95921 (2005).
(8) Jennings H. M., Ά model for the microstructure of calcium silicate hydrate in cernent paste', Cernent and Concrète Research, 30, p. 101-116 (2000).
(9) Jennings H. M., "Colloid model of C-S-H and implications to the problem of creep and shrinkage", Materials and Structures/ Concrète Science and Engineering, 37, p. 59-70 (2004).
(10) Constantinides G., "Invariant mechanical properties of Calcium-Silicate-Hydrates (CS-H) in Cement-Based materials: instrumented nanoindentation and microporomechanical modelling", PhD thesis, Massachusetts institute of Technology (2006).
(11) Jennings H. M., "Refînements to colloid model of C-S-H in cernent: CM-II", Cernent and Concrète Research, 38, p. 275-289 (2008).
(12) Bernard, O., Ulm, F.-J., Lemarchand, E. "A multiscale micromechanics-hydration model for the early-age elastic properties of cement-based materials", Cernent and Concrète Research, 33, p. 1293-1309 (2003).
(13) Ghabezloo S., "Association of macroscopic laboratory testing and micromechanics modelling for the évaluation of the poroelastic parameters of a hardened cernent paste", Cernent and Concrète research, 40 (8), p. 1197-1210 (2010).
(14) Zaoui A., "Continuum micromechanics: survey", Journal of Engineering Mechanics, 128 (8), p. 808-816 (2002).
(15) Rixom, R., Mailvaganam, N. "Chemical admixtures for concrète", Third édition, Spons Architecture Price Book, 456 p. (1999). (16) Pont Castenada, P., Willis, R. "The effect of spatial distribution on the effective behaviour of composite materials and cracked média", J. Mech. Phys. Solids, 43, p. 1919-1951 (1995).
(17) Hill, R. "The essential structure of constitutive laws for métal composites and polycrystals", J. Mech. Phys. Solids, 15, p. 79-95. (1967).
(18) Eshelby, J.D. "The détermination of the elastic fîeld of an ellipsoïdal inclusion", Proceedings of the Royal Society of London, 241, p. 376-392 (1957).
(19) Ulm, F.-J., Constantinides, G., Heukamp, F.H. "Is concrète a poromechanics material? - A multiscale investigation of poroelastic properties", Materials and Structures, 37 (265), p. 43-58 (2004).
(20) Ghabezloo S., "Micromechanics analysis of thermal expansion and thermal pressurization of a hardened cernent paste", Cernent and Concrète research, 41 (5), p. 520-532 (2011).
(21) Bourissai, M. "Comportement thermo-chimio-hydro-mécanique d'un ciment pétrolier au très jeune âge en conditions de prise HP/HT. Approche expérimentale et analyse par changement d'échelle", Thèse de doctorat, Université Paris Est, 246 p. (2010).

Claims

REVENDICATIONS
1. Procédé de détermination de paramètres mécaniques d'un système cimentaire de composition initiale Co et de finesse Φ, en fonction du temps, et en fonction de la finesse du système cimentaire, de la pression et/ou de la température, comprenant les étapes suivantes:
(A) déterminer le degré d'hydratation du système cimentaire en fonction du temps a(t) à partir de la vitesse des ondes de compression en fonction du temps Vp(t) mesurée dans un échantillon du système cimentaire, à une pression Pi et une température Ti;
(B) déterminer le degré d'hydratation a(t) en fonction de valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn ;
(C) déterminer la composition du système cimentaire en fonction du temps C(t) et en fonction de valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn, à partir du degré d'hydratation a (t) déterminé à l'étape B ;
(D) déterminer au moins un paramètre mécanique du système cimentaire en fonction du temps, et en fonction des valeurs souhaitées de finesse Φη du système cimentaire, de pression Pn et/ou de température Tn, à partir de la composition du système cimentaire C(t) déterminée à l'étape C.
2. Procédé selon la revendication 1 , comprenant en outre une étape initiale de mesure de la vitesse des ondes de compression en fonction du temps Vp(t) dans un échantillon du système cimentaire.
3. Procédé selon la revendication 1 ou 2, dans lequel ledit paramètre mécanique est choisi parmi les paramètres de déformabilité statiques, de préférence les paramètres élastiques statiques, tel que module de Young statique E, le coefficient de Poisson statique v, le module d'incompressibilité K, le module de cisaillement G, et leurs combinaisons.
4. Procédé selon la revendication 3, dans lequel deux paramètres élastiques statiques sont déterminés.
5. Procédé selon la revendication 3 ou 4, comprenant en outre la détermination d'un paramètre mécanique choisi parmi les paramètres de couplage hydro-mécanique, tel que le coefficient de Biot ou le coefficient de Skempton, et leurs combinaisons.
6. Procédé selon l'une des revendications précédentes, dans lequel le degré d'hydratation du système cimentaire en fonction du temps a(t) est calculé à partir de Vp(t) selon une relation linéaire.
7. Procédé selon la revendication 6, dans lequel le degré d'hydratation du système cimentaire en fonction du temps a(t) est calculé à partir de Vp(t) selon la relation a = (VP - V0 )/ (V∞ - V0 ) , avec V0 et Vp correspondant respectivement à la vitesse des ondes de compression mesurée dans l'échantillon du système cimentaire au temps t = 0 et au temps t, et V∞ correspondant à la vitesse des ondes de compression dans l'échantillon du système cimentaire complètement hydraté.
8. Procédé selon l'une des revendications précédentes, dans lequel le processus d'hydratation comprend un premier stade où l'hydratation est gouvernée principalement par un processus de nucléation et de croissance, et un deuxième stade où l'hydratation est gouvernée principalement par un processus diffusion des ions, ledit deuxième stade commençant lorsque le degré d'hydratation a atteint une valeur de degré d'hydratation seuil a*, cette valeur seuil a* étant fonction de la température, et dans lequel l'étape B comprend les sous-étapes suivantes :
(B-i) la détermination du degré d'hydratation a(t) pendant le premier stade du processus d'hydratation du système cimentaire ;
(B-ii) la détermination du degré d'hydratation a(t) pendant le deuxième stade du processus d'hydratation du système cimentaire ;
chacune des sous-étapes B-i et B-ii prenant en compte la finesse Φ du système cimentaire, la pression et la température pour la détermination du degré d'hydratation a(t).
9. Procédé selon la revendication 8, dans lequel la valeur de degré d'hydratation seuil a* est évaluée en minimisant l'écart entre a(t) déterminé en utilisant le modèle cinétique 120 et a(t) déterminé expérimentalement à partir de la vitesse des ondes de compression, pour différentes températures, et à une pression constante, de manière à prendre en compte la variation de a* en fonction de la température dans l'étape B.
10. Procédé selon la revendication 9, dans lequel :
- le système cimentaire de composition initiale Co comprend un ciment et de l'eau, le ciment comprenant au moins une phase initiale réactive X ;
- le degré d'hydratation a(t) déterminé à l'étape B correspond à la moyenne pondérée des degrés d'hydratation de chacune des phases initiales réactives X du ciment ;
- le degré d'hydratation de chacune des phases initiales réactives X du ciment est fonction du rapport entre l'affinité chimique t (° de la phase initiale réactive X, cette affinité chimique t (° contrôlant le taux d'évolution de l'hydratation de la phase initiale réactive X, et le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau τχ ; - et le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau τχ est fonction de la finesse Φ du ciment, de la pression et de la température.
11. Procédé selon la revendication 10, dans lequel le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau τ x s'exprime selon l'équation suivante :
Figure imgf000044_0001
avec Φο la finesse d'un ciment de référence, de préférence Φο = 3600 cm2/g, T0 la température initiale et P0 la pression initiales au temps t=0 du processus d'hydratation, R la constante de gaz parfaits, ΔΕΧ l'énergie d'activation, AVX le volume d'activation pour la nucléation et la croissance des hydrates lors du premier stade du processus d'hydratation et nx une constante.
12. Procédé selon l'une des revendications précédentes, dans lequel le système cimentaire comprend m phases, parmi lesquelles :
- au moins une phase initiale réactive X, de préférence choisie dans le groupe constitué par le silicate tricalcique C3S, le silicate dicalcique C2S, l'aluminate tricalcique C3A, l'alumino ferrite tétracalcique C4F, et leurs associations;
- au moins une phase hydratée Y issue de l'hydratation d'au moins une phase initiale réactive X, de préférence choisie dans le groupe constitué par le silicate de calcium hydraté C-S-H, l'hydroxyde de calcium CH, le trisulfoaluminate de calcium hydraté TSA, le monosulfate de calcium hydraté AFm, l'alumino ferrite de calcium hydraté, et leurs associations ;
- de l'eau;
la composition initiale du système cimentaire Co comprenant un volume initial d'eau déterminé V° et au moins une phase initiale réactive X,
et dans lequel l'étape C comprend l'estimation de la composition du système cimentaire en fonction du temps C(t) par la détermination de la fraction molaire, massique ou volumique des m phases du système cimentaire.
13. Procédé selon la revendication 12, dans lequel :
- le volume d'au moins une phase initiale réactive X en fonction du temps Fx (t) est calculé selon l'équation (XIV) :
^ ( = ^ (1 - 0, ( ) (xiv) avec Vx le volume initial de la phase initiale réactive X, et ax (i) le degré d'hydratation de la phase initiale réactive X en fonction du temps ;
- le volume d'au moins une phase hydratée en fonction du temps est calculé selon l'équation (XV) :
Figure imgf000045_0001
avec V* le volume occupé par la phase hydratée 7 formée par la phase réactive X dans un volume élémentaire représentatif, et ax (t) le degré d'hydratation de la phase réactive X en fonction du temps ;
- le volume de l'eau en fonction du temps Vw (t) est calculé selon l'équation (XVI) : (ή = ν -∑ν - αχ (ή (XVI) avec Vl le volume initial d'eau dans le système cimentaire, V^ le volume de l'eau consommée par la phase X, ax (t) le degré d'hydratation de X en fonction du temps ;
14. Procédé selon l'une des revendications précédentes, dans lequel l'étape D est réalisée selon une méthode d'analyse multi-compositionnelle (140), cette méthode d'analyse multi-compositionnelle (140) :
- prenant en compte un modèle multi-échelles du système cimentaire, ledit modèle multi-échelles comprenant au moins une échelle élémentaire représentant des globules de la phase solide C-S-H et une échelle macroscopique du système cimentaire, de préférence un modèle multi-échelles comprenant trois échelles ;
- et permettant la détermination d'au moins un paramètre mécanique du système cimentaire en fonction du temps, pour une valeur donnée de finesse du système cimentaire, de pression et/ou de température, en utilisant des techniques d'homogénéisation, de préférence le schéma de Mori-Tanaka et le schéma auto-cohérent, connaissant la composition du système cimentaire C(t) déterminée à l'étape C, et connaissant l'évolution des composants du ciment du système cimentaire obtenue à partir du modèle de la cinétique d'hydratation utilisé à l'étape A.
15. Procédé selon l'une des revendications précédentes, dans lequel le système cimentaire comprend du ciment Portland.
16. Procédé selon l'une des revendications précédentes, permettant de caractériser le comportement mécanique d'un système cimentaire utilisé en tant que gaine ou bouchon de ciment dans un puits, de préférence un puits pétrolier.
17. Produit programme informatique permettant de conserver dans une mémoire d'une unité d'un processeur ou sur un support de mémoire amovible approprié pour coopérer avec ladite unité du processeur, le produit-programme comprenant des instructions pour mettre en œuvre le procédé selon l'une quelconque de revendications 1 à 16.
PCT/FR2013/050355 2012-02-22 2013-02-21 Methode de caracterisation du comportement mecanique de ciments WO2013124588A1 (fr)

Priority Applications (8)

Application Number Priority Date Filing Date Title
DK13710496.4T DK2817615T3 (en) 2012-02-22 2013-02-21 Method for characterizing cement mechanical conditions
EP13710496.4A EP2817615B1 (fr) 2012-02-22 2013-02-21 Méthode de caractérisation du comportement mécanique de ciments
SG11201405034YA SG11201405034YA (en) 2012-02-22 2013-02-21 Method for characterising the mechanical behaviour of cements
BR112014020577-9A BR112014020577B1 (pt) 2012-02-22 2013-02-21 processo de determinação de parâmetros mecânicos de um sistema cimentário e suporte de memória apropriado para cooperar com um processador
RU2014138108A RU2628032C2 (ru) 2012-02-22 2013-02-21 Способ определения механических свойств цементов
US14/379,549 US9664648B2 (en) 2012-02-22 2013-02-21 Method for characterising the mechanical behaviour of cements
CA2864513A CA2864513C (fr) 2012-02-22 2013-02-21 Methode de caracterisation du comportement mecanique de ciments
CN201380019523.4A CN104335037A (zh) 2012-02-22 2013-02-21 水泥力学性能的表征方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1251612A FR2987128B1 (fr) 2012-02-22 2012-02-22 Methode de caracterisation du comportement mecanique de ciments
FR1251612 2012-02-22

Publications (1)

Publication Number Publication Date
WO2013124588A1 true WO2013124588A1 (fr) 2013-08-29

Family

ID=47901209

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR2013/050355 WO2013124588A1 (fr) 2012-02-22 2013-02-21 Methode de caracterisation du comportement mecanique de ciments

Country Status (10)

Country Link
US (1) US9664648B2 (fr)
EP (1) EP2817615B1 (fr)
CN (1) CN104335037A (fr)
BR (1) BR112014020577B1 (fr)
CA (1) CA2864513C (fr)
DK (1) DK2817615T3 (fr)
FR (1) FR2987128B1 (fr)
RU (1) RU2628032C2 (fr)
SG (1) SG11201405034YA (fr)
WO (1) WO2013124588A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017538131A (ja) * 2014-12-17 2017-12-21 ベンタナ メディカル システムズ， インコーポレイテッド 真の拡散定数の取得

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140321240A1 (en) * 2013-04-26 2014-10-30 Siemens Medical Solutions Usa, Inc. Elastography for cement integrity inspection
CN105259335B (zh) * 2015-10-26 2018-03-23 中国建筑材料科学研究总院 水泥熟料性能的评价方法及煅烧工艺优化方法
US11061011B2 (en) * 2016-04-15 2021-07-13 Total Se Method for determining a time window in which a casing pressure test can be performed without damaging a cement sheath
WO2017178860A1 (fr) 2016-04-15 2017-10-19 Total Sa Procédé de détermination de paramètre de plasticité de pâte de ciment hydratant
CN107575210A (zh) 2016-11-30 2018-01-12 中国石油天然气股份有限公司 一种油气井固井水泥石弹性的评价方法
WO2018156122A1 (fr) 2017-02-22 2018-08-30 Halliburton Energy Services, Inc. Régulation de la chaleur d'hydratation par caractérisation de constituants cimentaires
CN108519317B (zh) * 2018-04-24 2023-10-31 北部湾大学 直接拉伸荷载下的岩石应力-渗流耦合试验装置
CN109142530A (zh) * 2018-08-30 2019-01-04 华东理工大学 基于信息融合和椭圆成像的压力容器缺陷导波定位方法
CN109903819A (zh) * 2019-01-15 2019-06-18 沈阳理工大学 一种硅酸盐水泥复合胶凝材料水化动力参数的分析方法
US20210151133A1 (en) * 2019-04-05 2021-05-20 Halliburton Energy Services, Inc. Method For Designing For Temperature Sensitivity Of Hydration Of Cement Slurry
US11386247B2 (en) 2019-04-05 2022-07-12 Halliburton Energy Services, Inc. Method for designing for induction time using supplementary cementitious materials
WO2021067502A1 (fr) * 2019-09-30 2021-04-08 Saudi Arabian Oil Company Systèmes et procédés de calorimétrie isothermique à haute température et à haute pression
CN112751499B (zh) * 2019-10-29 2022-04-08 苏州慧闻纳米科技有限公司 基于微机电系统的摩擦纳米发电机及其制备方法
CN112782038B (zh) * 2019-11-07 2023-07-25 中国石油化工股份有限公司 一种用于泡沫水泥的声密度特性的评价装置及评价方法
CN111208281B (zh) * 2020-01-17 2022-06-24 西南石油大学 水泥水化过程中液塑固三态体积变化的测试装置及方法
CN111462839B (zh) * 2020-04-21 2023-10-13 浙江水利水电学院 一种硬化水泥砂浆热膨胀系数的多尺度预测方法
CN113722969B (zh) * 2021-08-12 2022-06-14 武汉理工大学 一种硅酸三钙水化全过程性能指标确定方法及其计算分析仪
CN114778301B (zh) * 2022-06-17 2022-09-06 煤炭科学研究总院有限公司 一种岩体稳定性的判定方法、装置及电子设备
CN116466066B (zh) * 2023-04-20 2023-10-31 深圳大学 压缩浇筑纤维混凝土中纤维最佳掺量及长细比的确定方法

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2538114A (en) 1944-10-17 1951-01-16 Bell Telephone Labor Inc Thickness measurement
US3401773A (en) 1967-12-04 1968-09-17 Schlumberger Technology Survey Method and apparatus for cement logging of cased boreholes
US4255798A (en) 1978-05-30 1981-03-10 Schlumberger Technology Corp. Method and apparatus for acoustically investigating a casing and cement bond in a borehole
US4259868A (en) 1979-10-01 1981-04-07 Halliburton Company Method and apparatus for nondestructive testing of cement
US4779236A (en) 1986-07-28 1988-10-18 Amoco Corporation Acoustic well logging method and system
US4813028A (en) 1987-07-07 1989-03-14 Schlumberger Technology Corporation Acoustic well logging method and apparatus
US5001676A (en) 1990-04-27 1991-03-19 Mobil Oil Corporation Acoustic borehole logging
US5168470A (en) 1992-01-23 1992-12-01 Mobil Oil Corporation Apparatus for rotating a transducer assembly of a borehole logging tool in a deviated borehole
US5357481A (en) 1992-11-04 1994-10-18 Western Atlas International, Inc. Borehole logging tool
US5412990A (en) * 1992-03-20 1995-05-09 Schlumberger Technology Corporation Method for measuring cement thickening times
US5763773A (en) 1996-09-20 1998-06-09 Halliburton Energy Services, Inc. Rotating multi-parameter bond tool
US5859811A (en) 1996-02-29 1999-01-12 Schlumberger Technology Corporation Method of analyzing waveforms
EP1541987A2 (fr) 2003-12-12 2005-06-15 Bj Services Company Dispositif d'essai et méthode de dérivation du module d'Young des rapports effort/contrainte de tension
US7089816B2 (en) 2004-01-13 2006-08-15 Halliburton Energy Services, Inc. Method and apparatus for testing cement slurries
WO2007020435A1 (fr) 2005-08-18 2007-02-22 Halliburton Energy Services, Inc. Appareil et procede permettant de determiner des proprietes mecaniques de ciment destine a un puits de forage
US7549320B2 (en) 2007-01-11 2009-06-23 Halliburton Energy Services, Inc. Measuring cement properties
US7552648B2 (en) 2007-09-28 2009-06-30 Halliburton Energy Services, Inc. Measuring mechanical properties
US7621186B2 (en) 2007-01-31 2009-11-24 Halliburton Energy Services, Inc. Testing mechanical properties

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU612167A1 (ru) * 1976-08-01 1978-06-25 Всесоюзный Научно-Исследовательский И Проектный Институт "Теплопроект" Способ определени тепловыделени при гидратации образца исследуемого вещества
SU1474548A1 (ru) * 1987-03-20 1989-04-23 Московский Инженерно-Строительный Институт Им.В.В.Куйбышева Способ определени степени гидратации цемента
RU2008673C1 (ru) * 1991-09-23 1994-02-28 Уральский государственный технический университет Способ определения гидрофобности портландцемента
US20020112540A1 (en) * 2000-12-20 2002-08-22 Schlumberger Technology Corporation Acoustic method for estimating mechanical properties of a material and apparatus therefor
US6941819B1 (en) * 2001-09-28 2005-09-13 Chandler Instruments Company L.L.C. Apparatus and method for determining the dynamic mechanical properties of a cement sample
EP1726947A1 (fr) * 2005-04-20 2006-11-29 Sika Technology AG appareil et procédé pour la determination ultrasonore de module d'élasticité dynamique
RU2370766C1 (ru) * 2008-01-17 2009-10-20 Евгений Сергеевич Шитиков Способ оценки активности цементов в водоцементных системах
FR2965925B1 (fr) * 2010-10-12 2015-06-19 Total Sa Mesure des proprietes d'echantillons de compositions durcissantes a pression elevee

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2538114A (en) 1944-10-17 1951-01-16 Bell Telephone Labor Inc Thickness measurement
US3401773A (en) 1967-12-04 1968-09-17 Schlumberger Technology Survey Method and apparatus for cement logging of cased boreholes
US4255798A (en) 1978-05-30 1981-03-10 Schlumberger Technology Corp. Method and apparatus for acoustically investigating a casing and cement bond in a borehole
US4259868A (en) 1979-10-01 1981-04-07 Halliburton Company Method and apparatus for nondestructive testing of cement
US4779236A (en) 1986-07-28 1988-10-18 Amoco Corporation Acoustic well logging method and system
US4813028A (en) 1987-07-07 1989-03-14 Schlumberger Technology Corporation Acoustic well logging method and apparatus
US5001676A (en) 1990-04-27 1991-03-19 Mobil Oil Corporation Acoustic borehole logging
US5168470A (en) 1992-01-23 1992-12-01 Mobil Oil Corporation Apparatus for rotating a transducer assembly of a borehole logging tool in a deviated borehole
US5412990A (en) * 1992-03-20 1995-05-09 Schlumberger Technology Corporation Method for measuring cement thickening times
US5357481A (en) 1992-11-04 1994-10-18 Western Atlas International, Inc. Borehole logging tool
US5859811A (en) 1996-02-29 1999-01-12 Schlumberger Technology Corporation Method of analyzing waveforms
US5763773A (en) 1996-09-20 1998-06-09 Halliburton Energy Services, Inc. Rotating multi-parameter bond tool
EP1541987A2 (fr) 2003-12-12 2005-06-15 Bj Services Company Dispositif d'essai et méthode de dérivation du module d'Young des rapports effort/contrainte de tension
US7089816B2 (en) 2004-01-13 2006-08-15 Halliburton Energy Services, Inc. Method and apparatus for testing cement slurries
WO2007020435A1 (fr) 2005-08-18 2007-02-22 Halliburton Energy Services, Inc. Appareil et procede permettant de determiner des proprietes mecaniques de ciment destine a un puits de forage
US7549320B2 (en) 2007-01-11 2009-06-23 Halliburton Energy Services, Inc. Measuring cement properties
US7621186B2 (en) 2007-01-31 2009-11-24 Halliburton Energy Services, Inc. Testing mechanical properties
US7552648B2 (en) 2007-09-28 2009-06-30 Halliburton Energy Services, Inc. Measuring mechanical properties

Non-Patent Citations (22)

* Cited by examiner, † Cited by third party
Title
AVRAMI, M., KINETICS OF PHASE CHANGE, JOURNAL OF CHEMICAL PHYSICS, vol. 7, 9, 1939, pages 1103 - 1124,177-184
BERNARD, O.; ULM, F.-J.; LEMARCHAND, E.: "A multiscale micromechanics-hydration model for the early-age elastic properties of cement-based materials", CEMENT AND CONCRETE RESEARCH, vol. 33, 2003, pages 1293 - 1309
BOUMIZ A.; VEET C.; COHEN TENOUDJIT F.: "Mechanical properties of cement pastes and mortars at early ages", ADVANCED CEMENT BASED MATERIALS, vol. 3, 1996, pages 94 - 106
BOURISSAI, M.: "Comportement thermo-chimio-hydro-mécanique d'un ciment pétrolier au très jeune âge en conditions de prise HP/HT. Approche expérimentale et analyse par changement d'échelle", THÈSE DE DOCTORAT, 2010, pages 246
CONSTANTINIDES G.: "Invariant mechanical properties of Calcium-Silicate-Hydrates (CS-H) in Cement-Based materials: instrumented nanoindentation and microporomechanical modelling", PHD THESIS, 2006
ESHELBY, J.D.: "The determination of the elastic field of an ellipsoidal inclusion", PROCEEDINGS OF THE ROYAL SOCIETY OF LONDON, vol. 241, 1957, pages 376 - 392
FUJI, K.; KONDO, W.: "Kinetics of the hydration of tricalcium silicate", JOURNAL OF THE AMERICAN CERAMIC SOCIETY, vol. 57, no. 11, 1974, pages 492 - 497
GHABEZLOO S.: "Association of macroscopic laboratory testing and micromechanics modelling for the evaluation of the poroelastic parameters of a hardened cement paste", CEMENT AND CONCRETE RESEARCH, vol. 40, no. 8, 2010, pages 1197 - 1210, XP027103113
GHABEZLOO S.: "Micromechanics analysis of thermal expansion and thermal pressurization of a hardened cement paste", CEMENT AND CONCRETE RESEARCH, vol. 41, no. 5, 2011, pages 520 - 532, XP028156751, DOI: doi:10.1016/j.cemconres.2011.01.023
H LEE: "Ultrasonic in-situ monitoring of setting process of high-performance concrete", CEMENT AND CONCRETE RESEARCH, vol. 34, no. 4, 1 April 2004 (2004-04-01), pages 631 - 640, XP055006415, ISSN: 0008-8846, DOI: 10.1016/j.cemconres.2003.10.012 *
HAECKER C.J.; GARBOCZI E.J.; BOHN J.W.; SUN Z.; SHAH S.P.; VOIGT T.: "Modeling the linear elastic properties of Portland cement paste", CEMENT CONCRETE RESEARCH, vol. 35, 2005, pages 1948 - 1960, XP005090407, DOI: doi:10.1016/j.cemconres.2005.05.001
HILL, R.: "The essential structure of constitutive laws for metal composites and polycrystals", J. MECH. PHYS. SOLIDS, vol. 15, 1967, pages 79 - 95, XP024469378, DOI: doi:10.1016/0022-5096(67)90018-X
JENNINGS H. M.: "A model for the microstructure of calcium silicate hydrate in cement paste", CEMENT AND CONCRETE RESEARCH, vol. 30, 2000, pages 101 - 116
JENNINGS H. M.: "Colloid model of C-S-H and implications to the problem of creep and shrinkage", MATERIALS AND STRUCTURES/ CONCRETE SCIENCE AND ENGINEERING, vol. 37, 2004, pages 59 - 70
JENNINGS H. M.: "Refinements to colloid model of C-S-H in cement: CM-II", CEMENT AND CONCRETE RESEARCH, vol. 38, 2008, pages 275 - 289, XP022432803, DOI: doi:10.1016/j.cemconres.2007.10.006
JENNINGS H. M.; TENNIS P.D.: "Model for the developing microstructure in Portland Cement Pastes", J. AM. CERAM. SOC., vol. 77, no. 12, 1994, pages 3161 - 3172
KONDO R.; KODAMA, M., ON THE HYDRATION KINETICS OF CEMENT, SEMENTO GIJUTSU NENPO, vol. 21, 1967, pages 77 - 828
PONT CASTENADA, P.; WILLIS, R.: "The effect of spatial distribution on the effective behaviour of composite materials and cracked media", J. MECH. PHYS. SOLIDS, vol. 43, 1995, pages 1919 - 1951
REDDY, B.R.; SANTRA, A.; MCMECHAN, D.; GRAY, D.; BRENNEIS, C.; DUNN, R.: "Cement mechanical property measurements under wellbore conditions", SPE, 2005, pages 95921
RIXOM, R.; MAILVAGANAM, N.: "Spons Architecture Price Book", 1999, article "Chemical admixtures for concrete", pages: 456
ULM, F.-J.; CONSTANTINIDES, G.; HEUKAMP, F.H.: "Is concrete a poromechanics material? - A multiscale investigation of poroelastic properties", MATERIALS AND STRUCTURES, vol. 37, no. 265, 2004, pages 43 - 58
ZAOUI A.: "Continuum micromechanics: survey", JOURNAL OF ENGINEERING MECHANICS, vol. 128, no. 8, 2002, pages 808 - 816

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017538131A (ja) * 2014-12-17 2017-12-21 ベンタナ メディカル システムズ， インコーポレイテッド 真の拡散定数の取得
US10620037B2 (en) 2014-12-17 2020-04-14 Ventana Medical Systems, Inc. Obtaining true diffusivity constant

Also Published As

Publication number Publication date
CA2864513C (fr) 2020-04-28
BR112014020577B1 (pt) 2020-12-22
EP2817615A1 (fr) 2014-12-31
CN104335037A (zh) 2015-02-04
US9664648B2 (en) 2017-05-30
FR2987128B1 (fr) 2014-02-28
DK2817615T3 (en) 2018-08-06
SG11201405034YA (en) 2014-10-30
US20150033862A1 (en) 2015-02-05
RU2628032C2 (ru) 2017-08-14
CA2864513A1 (fr) 2013-08-29
FR2987128A1 (fr) 2013-08-23
EP2817615B1 (fr) 2018-05-02
RU2014138108A (ru) 2016-04-10

Similar Documents

Publication Publication Date Title
EP2817615B1 (fr) Méthode de caractérisation du comportement mécanique de ciments
Khalifeh et al. Geopolymers as an alternative for oil well cementing applications: A review of advantages and concerns
CN101539566B (zh) 混凝土早期毛细管负压的测试方法
Rabbi et al. Effect of curing stress and period on the mechanical properties of cement-mixed sand
Karakosta et al. Pore structure evolution and strength development of G-type elastic oil well cement. A combined 1H NMR and ultrasonic study
Ye et al. Experimental study on ultrasonic pulse velocity evaluation of the microstructure of cementitious material at early age.
Romanowski et al. Investigations on oilwell cement strength response to ultrasonic measurements in the presence of additives
Moradi-Marani et al. Effect of the temperature on the nonlinear acoustic behavior of reinforced concrete using dynamic acoustoelastic method of time shift
Anya et al. A novel apparatus and method for lab-scale study of wellbore integrity using CT imaging and analysis
Sun et al. Experimental Study of the Effect of Curing Temperature and Pressure on the Property Evolution of Oil Well Cement
Samudio Modelling of an oil well cement paste from early age to hardened state: hydration kinetics and poromechanical behaviour
Barría et al. Cement with bacterial nanocellulose cured at reservoir temperature: Mechanical performance in the context of CO2 geological storage
Elkatatny Improved carbonation resistance and durability of Saudi Class G oil well cement sheath in CO2 rich environments using laponite
Kongsukprasert et al. Small strain stiffness and non-linear stress-strain behaviour of cement-mixed gravelly soil
Sun et al. Early age properties of Portland cement pastes investigated with ultrasonic shear waves and numerical simulation
Kelessidis et al. Comprehensive assessment of additive and class g cement properties affecting rheology, fluid loss, setting time and long term characteristics of elastic cements
Akgün et al. Influence of degree of saturation on the borehole sealing performance of an expansive cement grout
Akintola et al. Performance Evaluation of Local Material Rice Husk Ash under Downhole Conditions with the Addition of Basic Oil Well Additives Antifoam, Fluid Loss, Dispersant and Retarder on Oil Well Cementing
Agofack et al. Effect of loading history on mechanical properties of hardened oilwell cement
Todorovic et al. Effect of Long-Term Aging in Carbonated Brine on Mechanical Properties of a Novel Cement System with an Expandable Agent
Alexa et al. Experimental Determination of Early Shrinkage of Alkali-Activated Slag
Sleiman et al. A numerical microscopically informed upscale approach for analyzing the reliability of testing method for concrete resistance to freeze-thaw
Haghighat-Nejad Flexible Geopolymer for Oil and Gas Well Cementing: An Experimental Study
Khial et al. Assessment of the aggressive agents' penetration into concrete by non-destructive techniques
Zhang et al. Deterioration of the performance of cement mortar lining due to sulfate attack

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13710496

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2013710496

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2864513

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 14379549

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2014138108

Country of ref document: RU

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112014020577

Country of ref document: BR

ENP Entry into the national phase

Ref document number: 112014020577

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20140821