Brevet FR2987128A1 — publié le 20130823

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ée 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 CO2. 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 connaitre 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, tel 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 oedomé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 part 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 Cement 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-10 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 15 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 20 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. 25 BREVE DESCRIPTION DES DESSINS La figure 1 représente le logigramme du procédé selon l'invention. 30 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 Cement 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. 35 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 (I), 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 P1 et une température T1; (B) déterminer le degré d'hydratation a(t) en fonction de valeurs souhaitées de finesse (I). du système cimentaire, de pression Po et/ou de température To ; (C) déterminer la composition du système cimentaire en fonction du temps C(t) et en fonction de valeurs souhaitées de finesse (I)n 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 (1)n 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 oeuvre 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 C25, l'aluminate tricalcique de formule 3CaO.A1203 (Aluminate), de symbole C3A, l'aluminoferrite 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-SH, l'hydroxyde de calcium CH, le trisulfoaluminate de calcium hydraté TSA (Ettringite), le monosulfate de calcium hydraté AFm, l'aluminoferrite de calcium hydraté, et leurs associations. Les phases initiales telles que les silicates C3S, C25 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'aluminoferrite 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 Phase Fraction Possibilité de Morphologies Dimensions minéralogique volumique (%) cristallisation C-S-H 50-70 Très peu Non résolu 10-100 nm CH 15-20 Très bien Solide hexagonal 10-100 nm Ettringite 1-5 Bien Solide de forme des 10x0,5 nm aiguilles Monosulfoaluminate 1-5 bien Solide sous forme des plans lx1x0,1 nm hexagonaux 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 l'Ettringite 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.("), 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 100nm (Constantinides(1°)), 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 C3 S ; - du silicate bicalcique de formule 2CaO.Si02 (Belite), de symbole C25 ; - de l'aluminate tricalcique de formule 3CaO.A1203 (Aluminate), de symbole C3A; - du ferro-aluminate tétracalcique de formule 4CaO.A1203 Fe2O3 (Ferrite), de symbole C4F; - du sulfate dihydraté de calcium CaSO4.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.1MPa. 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. '5)). 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 (Na2SO4, K2504, CaSO4.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=oo) 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 Avramim, 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 cc*. 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(4) ; Fuji et Kondo(5)). Dans le cas du processus d'hydratation d'un ciment Portland, tel que représenté à la figure 2, le degré d'hydratation seuil cc* 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 Avramim, 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): 2C3S +10,6H ->C34 - S - H8 2,6CH 2C2S + 8,6H C34 - S2 - H8 0,6CH (I) 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 C25, 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 C25. 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 (I). du système cimentaire, de pression P. et/ou de température T.. 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 oedométrique Kv, ou le module de Lamé X. 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 K E 2t, y K, G 2G G 3A+2G - + 2G - A + 3 A +G 2(À, +G) - K - A - 3K - 2A 3 9K -(KA) 2 3K - A 3K - A - 9KG 2G 3K -2G K - K +4G - 3K +G 3 2(3K +G) 3 EG - GE-2G E G4G - E - -1 3(3G- E) 3G- E 2G 3G- E - - 3K3K-E 3K-E 3K3K +E 3KE 9K - E 6K 9K - E 9K -E Al+ v (1+v)(1-2v) Â - - Al-v A1-2v 3v v y 2v 1+v 2G(1+v) v - 1-v - 2G 2G 2G 3(1-2v) 1-2v 1-2v - 3K (1- 2v) 3K y - 3K 1-v 1- 2v 3K 1+v 1+v 2(1+v) E _ Ev _ E (1-v) E 3(1-2v) (1+v)(1- 2v) (1+v)(1-2v) 2(1+v) 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. 3y. - (1+ )- bB 3 - 2(1+ )- bB 3v d - (1+ 2v d)- bB 3 - (1- 2yd )- bB Ea - 3,(1- bB,) E - 3 - 11+ y. )- bB E - 3 -(1- 2yd )- bB K d = (1- bB) - K. K - d u 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 25 é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 = Ab-axial / 4e.,,, y = -4E',, / 4eaxial , avec An - - axial l'incrément de contrainte axiale, A -Eaxial l'incrément de déformation axiale et A -Eradial 30 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ù va - vu - 20 3Ed les déformations sont réversibles, par la formule suivante : K =46somopeAE volumique avec AGisotrope 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 oedométrique, K, est une constante élastique classiquement déterminée lors d'un essai de compression oedomé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 = 46,,, / 4e.a, , avec An -- axial l'incrément de contrainte axiale et Asaxial 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 : 46' = 46 - b - Pi), où a' est la contrainte effective, a 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 4/3, =B - 46 , dans laquelle 46 est la variation de contrainte isotrope , et APp 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 (I). 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 des 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 enfonction 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 Cement 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 1 : 2 V dyn 2 ( (HI) Vs Vp ) - 2 ( Vs 1 V P Module de Young dynamique : Edyn = densité x Vp(1+ v dyn ). (1-2 . vdyn) (1- vdyn (IV) 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, T1). 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 : US5859811, US5763773, US5357481, US5168470, US5001676, US4813028, US4779236, US4255798, US3401773, US2538114. Il est également évident que le procédé selon l'invention peut être mis en oeuvre 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 T1. Pour cela, le procédé selon l'invention utilise une corrélation empirique 110 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 ) a - (vo, - vo ) (V) avec Vo 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 Vs, 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 (I). du système cimentaire, de pression P. et/ou de température T.. 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 cc* é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 (1) 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 cc* 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 cc* 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 Tx ; - et le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau rx est fonction de la finesse (I) 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 ix s'exprime selon l'équation suivante : _ n,+ 1 Rx ( 4Rx ( (VI) Tx (T ,a) X Tx (T0,(0)exp T T 0 ) T T 0 avec (I)0 la finesse d'un ciment de référence, de préférence (I)0 = 3600 cm2/g, To la température initiale et Po la pression initiales au temps t=0 du processus d'hydratation, R la constante de gaz parfaits, 4E' l'énergie d'activation, 4V' le volume d'activation pour la nucléation et la croissance des hydrates lors du premier stade du processus d'hydratation, n,, une constante, déterminée expérimentalement dont les valeurs sont publiées dans la littérature (Bernard et al.(12)). Sous-étape 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 Avramim. 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 : a Inix, a x x=1 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, C,AF ). m'' 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: dax Ax (a) (VIII) dt Tx L'affinité chimique A' (a) de la phase initiale réactive X vérifie l'équation suivante : Ax (a) = (1-(a -a0x))-[-ln (1- (a -a0x))11n 1(IX) 10 où n, et aox sont des constantes (Bernard et al.(12)). de déterminées de expérimentalement la température. 15 l'homme de l'art caractéristique associé ( AE x à la réaction la phase et de AV x initiale et connues par Le temps de la finesse (I) du ciment, la pression P P réactive X avec l'eau Tx est fonction vérifie l'équation suivante : 0 Au cours de la sous-étape B-i, ( 1 1 T T (VI) T (T ,a)) _ T 0 ) 0 )) n,+1 (1)0 XTx )exp (1)0 R R (I) Dans laquelle (Do est la finesse d'un ciment de référence, de préférence (1)0=3600cm2/g, T0 la température initiale et Po la pression initiales au temps t=0 du processus d'hydratation, R la constante de gaz parfaits, égale à 8,314 J°K-1mo1-1, 4E' l'énergie d'activation, 4Vx le 20 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 à -27x10-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 25 la cinétique d'hydratation de ce dernier sont intégrés dans cette équation (V/). 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-étape B-ii 30 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 35 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 :5 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 Al(a) (XI) dt L'affinité chimique vérifie l'équation suivante : A(0 = (1- 02" Comme précédemment, on peut exprimer le temps caractéristique Tx associé à la réaction de la phase initiale réactive X avec l'eau en fonction de la finesse (I) du ciment, de la pression et de la température. rx vérifie l'équation (V/) précédemment explicitée. Au cours de ce deuxième stade, le temps caractéristique de référence 2x (T0,(1)0) est contrôlé par la diffusion des ions à travers la couche des hydrates autour des grains de ciment (Bernard et al. '2) : ix (T("(1)°) R2 (XIII) 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 4V' 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 cc, 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. 23 (X) (1-a *)1/3 -(1-a )1/335 Etape C 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 (I). du système cimentaire, de pression P, et/ou de température T., à 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 connaitre, à 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 (I). du système cimentaire, de pression P. et/ou de température T.. 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 C25, l'aluminate tricalcique C3A, l'aluminoferrite 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'aluminoferrite 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,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 West calculé selon l'équation (XIV) : Vx (t) Vx° (1-a, (t)) (XIV) avec T7,,° le volume initial de la phase initiale réactive X, et a \ (t) 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) : v(t)=Ivyx - ax (t) (XV) avec VxY 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 Vi,' West calculé selon l'équation (XVI) : (t) = -ax (t) (XVI) avec V,v° le volume initial d'eau dans le système cimentaire, Vwx 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 5 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) = T7-s-H (et * -(e)( * (t))x Hia * (t)i) (XVII) 10 VHD (t) = (7-s-H (e)( (t) - a *) x Hia (t) - a *.1 (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. 15 La densité de C-S-H HD et C-S-H LD est exprimée par: 20 P HD HD) - Pg HD - PW P LD LD) - Pg +0 LD - PW 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 OHD et OLD 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 25 température, via le degré d'hydratation seuil cc* 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é 30 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 à 35 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 (XX/) permettant de calculer les fractions 5 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 T1, sous couvert que la différence de porosité capillaire soit connue. fLD 2 fLD 1 ( * ( 1 - a 1* ^ 1* 1 - a 2 a 2 a2 al+1\ 1 - HD2 - LD) \- HD1 - LD) 10 fHD 2 fHD1 (XXI) * a2 * ( * 1a 2 (1 ail 1-a1 1 - 1 HD1l-TLD) ± 1 0HD 2 1 OLD ) J LD20 LD f HD20 HD2 f LD10LD f HD10HD1 - cp1 cp2 15 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% à une température de à 60°C. 20 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 (I). du 25 système cimentaire, de pression P. et/ou de température T., à 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 multicompositionnelle 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 30 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 35 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 Zaoui(14), 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. '6)) 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 Hill(17) 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(18) 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 : C hOM 3K dh OM J + 2 G hOM K (XXII) Ô Ô , 61k6j/ +6116,k Où Jijk/ = 11. k" ; K - I J ; I= . (Sli dénote le delta Kronecker, K dhmn est le 3 2 module d'incompressibilité drainé, Gh'n 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 : n K dho. = frk rArv (XXIII) r =1 (XXIV) n Ghom = f r o. Ad J r r =1 Dans lesquelles fr, kr, et g, 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 sont données par les équations (XXV) et (XXVI) suivantes : (1+64 (krIk0-1))-1 Av = (XXV) fr +a0Vfr I k0 -1))-1 r Ard (1+ po(gr 1 g 0 -Or (XXVI) r(1+ f30(gr I go -1))-1 r 3k, Avec a0 - Po 504 + 4g0) 6,(k, +2g 0), ko le module d'incompressibilité de la 3k, + 4g, matrice de référence, et go le module de cisaillement de la matrice de référence. 30 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: 1 n-1 b ho. 1_1f, (A)V, (XXVII) r=1 Le signe (A1),, dénote la valeur moyenne volumique de A dans le volume V. Le 10 module de Biot homogénéisé est écrit sous la forme (XXVIII) suivante : 1 n-1 fr Cr-1 : -1 :(A) v,) (XXVIII) Nhoen r =1 Où cr est le tenseur de rigidité de la phase r qui est définit par la formule (XXIX) suivante : 15 cr= 3/fr.' + 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 20 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', , csir qui sont respectivement le tenseur de rigidité, le 25 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.(")) : 30 n-1 bhom =1- (A)v, : (1- bl r =1 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 (XXX/) suivante (Ulm et al.(")) : 35 1 n-1 \- 1 1 fr(f c') :1-1:(A)v,):(1-brI)+ r=1 (XXXI) NI 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.0 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.1 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.1, 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 al(12)). 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 . Kxh°111 = (1- Ox)ksAlelx (XXXII) Gx111 = (1- Ox) gsAe,x (XXXIII) 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.0 (« 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 (XXV/) décrites ci- dessus, avec ko=k, 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 (Ghabezloo"). 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 : K;in = AL)Kle AILD cp f HDKIHr A1HD cp f CH k CH A1CH cp f CK CK CK A1CK cp (XX)M7) G;in =VLD f AL)G1L'°Din AIL D cp f HD HDD cp+ iCHgCHACv H cp+ iCKgCKACvK cp Les tenseurs de localisation sont déterminés selon les équations (XXV) et (XXV/) décrites ci-dessus, avec ko=Kch: n , et go =Gchp'n (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 oeuvre 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 Nom chimique Formule chimique Symbole Pourcentage de masse (%) Tricalcium silicate (Alite) 3CaO.Si02 C3 S 55 Dicalcium silicate (Belite) 2CaO.Si02 C2 S 18 Tricalcium aluminate 3CaO.A1203 C3 A 10 Tetracalcium aluminoferrite (Ferrite) 4CaO.A1203 Fe2O3 C4AF 8 Calcium sulfate dihydrate (gypsum) CaSO4.2H20 CSH2 6 La composition initiale Co du système cimentaire ceml est donnée dans le tableau 4 ci-dessous. Elle correspond à un rapport eau sur ciment égal à environ 044 (w/c en anglais). TABLEAU 4 Composant Quantité [gr] Ciment 914,12 Eau 396,13 Agent anti-mousse 7,31 (D047) Dispersant 11,06 (D604AM) Agent anti-sédimentation 1,37 (D153) 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 cem 1 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 !amis. 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 a(t) à PI Tl, et à différentes valeurs de finesse On du système cimentaire, de pression Pn et de température Tn 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 T1, tel que décrit dans le test UCA ci-dessus. Voc 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 a = (T7,, - Vo )i (T7c,, - Vo ) Ce degré d'hydratation en fonction du temps a(t) est déterminé pour une pression Pi, et une température T1, avec Pi = 0.3 ou 20 MPa et T1 = 25 °C. La valeur de Tx ) 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. TABLEAU 5 Clinker I', To nx+1 cco AExIR (X) C3S 13,24 1,76 0,02 4800 C2S 72,01 1,00 0 2500 C3A 59 0,90 0,04 5500 C4AF 24,68 2,34 0,4 4200 La valeur de Tx ) utilisée dans la formule (V/) 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 Clinker C3S C2S C3A C4AF tx (h) 3492 564939 1439 3433 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 4V exprimée dans la formule (V/). TABLEAU 7 T °C) 7 13 15 25 30 40 60 Pression 0,3 0,3 40 0,3 20 40 0,3 40 0,3 20 (MPa) a' 0,68 0,70 0,68 0,72 0,72 0,72 0,60 0,60 0,40 0,40 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'Avramim. 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%, 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, CS-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.40 TABLEAU 8 C3S C2S C3A C4AF CH Globule C-S-H k [GPa] 112,5 116,7 120,8 104,2 32,5 23,0 g [GP] 51,9 53,8 55,8 48,1 14,6 17,0 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. u) ; 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 : E statique = 0.83 x Edynamique 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 1000h (42 jours) est de 6,1GPa. 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/Gdmesures » 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 ages", 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 concrete 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., 'A model for the microstructure of calcium silicate hydrate in cernent paste', Cernent and Concrete 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/ Concrete Science and Engineering, 37, p. 59-70 (2004). (10) Constantinides G., "Invariant mechanical properties of Calcium-Silicate-Hydrates (CS-H) in Cernent-Based materials: instrumented nanoindentation and microporomechanical modelling", PhD thesis, Massachusetts institute of Technology (2006). (11) Jennings H. M., "Refinements to colloid model of C-S-H in cernent: CM-II", Cernent and Concrete 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 cernent-based materials", Cernent and Concrete Research, 33, p. 1293-1309 (2003). (13) Ghabezloo S., "Association of macroscopic laboratory testing and micromechanics modelling for the evaluation of the poroelastic parameters of a hardened cernent paste", Cernent and Concrete 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 concrete", Third edition, 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 media",J. Mech. Phys. Solids, 43, p. 1919-1951 (1995). (17) Hill, R. "The essential structure of constitutive laws for metal composites and polycrystals", J. Mech. Phys. Solids, 15, p. 79-95. (1967). (18) Eshelby, J.D. "The determination of the elastic field of an ellipsoidal inclusion", Proceedings of the Royal Society of London, 241, p. 376-392 (1957). (19) Ulm, F.-J., Constantinides, G., Heukamp, F.H. "Is concrete 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 cement paste", Cement and Concrete 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). REVENDICATIONS1. Procédé de détermination de paramètres mécaniques d'un système cimentaire de composition initiale Co et de finesse (I), 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 P1 et une température T1; (B) déterminer le degré d'hydratation a(t) en fonction de valeurs souhaitées de finesse (Do du système cimentaire, de pression Po et/ou de température To ; (C) déterminer la composition du système cimentaire en fonction du temps C(t) et en fonction de valeurs souhaitées de finesse (Do du système cimentaire, de pression Pn et/ou de température To, à 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 (I)n du système cimentaire, de pression Po et/ou de température To, à 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 =(T7:, -T70)1(T70, -Vo), avec Vo et V, 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 rs, 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 (1) 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 A (a) de la phase initiale réactive X, cette affinité chimique A (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 iX ;- et le temps caractéristique associé à la réaction de la phase initiale réactive X avec l'eau r, est fonction de la finesse (I) 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 -c x s'exprime selon l'équation suivante : ( ( \I 4E 1 1 AV ( P P Tx (T, (1)) ± n+1 (1)0 xTx 7,-, yo (1)0)exp x Avec (Do la finesse d'un ciment de référence, de préférence (Do = 3600 cm2/g, To la température initiale et Po la pression initiales au temps t=0 du processus d'hydratation, R la constante de gaz parfaits, 4E' l'énergie d'activation, 4V' le volume d'activation pour la nucléation et la croissance des hydrates lors du premier stade du processus d'hydratation et n,, 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'aluminoferrite 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'aluminoferrite 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é Vw° 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 V, West calculé selon l'équation (XIV) : vx (t), vx° (i-ax (t)) (XIV) x o (I) R T T R T T T° ) o ) (VI)avec 17,,° le volume initial de la phase initiale réactive X, et a \ (t) 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) : v(t)=Ivyx - ax (t) avec VxY le volume occupé par la phase hydratée Y 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 Vi,' (t) est calculé selon l'équation (XVI) : (t) = -ax (t) (XVI) avec 17,4) le volume initial d'eau dans le système cimentaire, Vwx 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 oeuvre le procédé selon l'une quelconque de revendications 1 à 16.