FR 2734315 A1 19961122 FR 9505825 A 19950515 1 2734315 La présente invention concerne une méthode permettant de déterminer les conditions de forage d'un outil comportant plusieurs taillants en interaction avec une roche. La méthode comporte l'utilisation d'un modèle de foration fondé sur le couplage des effets de destruction de la roche par les taillants et des effets de l'évacuation des déblais par un fluide. L'invention s'applique de préférence à l'étude du phénomène de bourrage d'un outil du type PDC. Le bourrage est un dysfonctionnement fréquemment observé par le foreur et 1 0 très néfaste car il peut diminuer très fortement la vitesse d'avancement du forage et parfois même, dans certains terrains, annihiler irréversiblement les effets de foration. Plusieurs travaux ont déjà été présentés, mais aucun ne prend en compte le dégagement de matière comme le fait la représentation modélisée dans la présente méthode. Les principaux travaux sont cités dans la liste des références ciincluse. Ainsi, la présente invention concerne une méthode permettant d'améliorer les performances de forage dans laquelle on met en oeuvre un modèle de foration. Le modèle prend en compte les effets de destruction d'une roche par au moins un taillant fixé sur un corps d'outil entraîné en rotation et les effets d'évacuation des débris de roche par un fluide, en calculant un bilan matière à partir: - de la production de débris de roche par le taillant ayant pénétré dans la roche d'une profondeur 6, - d'un lit de débris recouvrant la roche sous une épaisseur 1, - d'une lame de fluide d'épaisseur h comprise entre le lit de débris et le corps, la lame fluide ayant une concentration c en débris, - de paramètres de commandes, - de paramètres d'environnement. 2 2734315 Par la méthode, on détermine les conditions de forage en fonction de la réponse du modèle pour des valeurs déterminées desdits paramètres. L'un au moins des paramètres: poids sur l'outil, vitesse de rotation de l'outil et débit de fluide, peut être un paramètre de commande. Dans le modèle, la portance W de l'outil peut être décomposée en une composante solide Ws et une composante hydraulique Wh fonction notamment de la lame fluide. On peut considérer un étalement granulométrique des débris distribués selon une loi normale fonction de la profondeur 6 de l'entaille, de moyenne g lié à la ductilité de la roche et d'une dispersion caractérisée par l'écart-type a. 1 0 Le bilan matière solide B(t) peut être tel que B(t)=B+(t)- B-(t), dans lequel B+(t) est un terme de production de débris dépendant de È et correspondant au rythme de destruction de la roche, et B-(t) est un terme d'expulsion dépendant de I et h. On peut représenter le forage comme un système dynamique comportant, dans la 1 5 représentation interne classique par variables d'état x, des entrées u qui seront celles d'un système de commandes "poids sur outil", vitesse de rotation des tiges, puissance hydraulique, système soumis également à des perturbations non contrôlables v associées à la variabilité des propriétés des roches. Par le moyen du présent modèle, le système est observé grâce aux variables de sortie y qui peuvent être, entre autres, le couple au niveau de l'outil, la vitesse de pénétration dans l'axe &u trou, des indicateurs liés au niveau vibratoire comme l'élargissement de diamètre du trou, des indicateurs d'usure des taillants de la tête de forage, malheureusement difficiles à concevoir, l'ensemble de ces variables de sortie pouvant être perturbé par un bruit w. L'optimisation de la foration peut ainsi être la recherche d'une stratégie de commande qui permette au foreur: - d'éviter les risques relatifs à des écueils localisés, par exemple, liés à des intercalations de roche très dures ou, à l'autre extrême, susceptibles d'entraîner le bourrage de l'outil, 3 2734315 -d'avoir une stratégie cohérente à l'échelle du forage: par exemple, la détermination du nombre et de la durée d'utilisation optimaux des outils de forage, ou la nécessité d'une adaptation de la conduite du forage au fur et à mesure de l'usure des taillants. Il est clair également que la présente méthode peut permettre d'aider à une détermination de la structure des outils de forage: par exemple, forme et implantation des taillants, détermination des écoulement hydrauliques dans le voisinage de la destruction de la roche. Les références suivantes, peuvent servir d'illustration de l'arrière-plan 1 0 technologique du domaine concerné ainsi que de compléments à la description de la présente invention. Andersen E.E. and Azar J.J., 1990, "PDC performance under simulated borehole conditions" SPE 20412, New Orleans sept. 1990. Cheatham C.A. and Nahm J.J., 1990, "Bit balling in water-reactive shale during 1 5 full-scale drilling rate tests" IADC/SPE n 19926, Houston. Déliac E.P., 1986, "Optimisation des machines d'abattage à pics" Doctoral Dissertation, U. Paris 6 ed by ENSMP/CGES France. Detournay E. and Atkinson C., 1991, "Influence of pore pressure on the drilling response of PDC bits", Rock Mechanics as a Multidisciplinary Science, Roegiers (ed.), Rotterdam. Detournay E. and Defourny P., 1992, "A Phenomenological Model for the Drilling Action of Drag Bits", Int. J. Rock Mech. Min. Sci. & Geomech. Abstr. Vol. 29, No. 1, p13-23. Falconer I.G., Burgess T.M. and Sheppard M.C., 1988, "Separating Bit and Lithology Effects from Drilling, Mechanics Data", IADC/SPE Drilling Conference, Dallas, feb 28-march 2, 1988. Garnier A.J. and van Lingen N.H., 1958, "Phenomena affecting drilling rates at depth "SPE fall meeting, Houston. 4 2734315 Glowka D.A., 1985, "Implications of Thermal Wear Phenomena for PDC Bit Design and Operation", 60th Annual Technical Conference and Exhibition of the Society of Petroleum Engineers in Las Vegas, sept 22-25, 1985, SPE 14222. Karasawa H. and Misawa S., 1992, "Development of New PDC Bits for Drilling of Geothermal Wells - Part 1: Laboratory Testing", Journal of Energy Resources Technology, dec 1992, vol. 114 p 323. Pessier R.C. and Fear M.J., 1992, "Quantifying common drilling problems with mechanical specific energy and a bit specific coefficient of sliding friction" SPE 24584. Pessier R.C., Fear M.J., and Wells M.R., 1994, "Different shales dictate 1 0 fundamentally different strategies in hydraulics, bit selection, and operating practices". Pierry J. and Charlier R., 1994, "Finite element modelling of shear band localization and application to rock cutting by a PDC tool" SPE/ISRM Eurock Conference, Delft. Putot C., 1995,"Un modèle de foration prenant en compte les effets de destruction 1 5 de la roche et d'évacuation des déblais", 2e Colloque national en calcul des structures, Giens. Sellami H., 1987, "Etude des pics usés, application aux machines d'abattage" Doctoral Dissertation ENSMP/CGES France. Sellami H., Fairhurst C., Déliac E. and Delbast B., 1989, "The Role of in-situ Stresses and Mud Pressure on the Penetration of PDC bits" Rock at Great Depth, Maury & Fourmaintraux eds, Rotterdam 1989. Sinor A. and Warren T.M., 1989, "Drag Bit Wear Model", SPE Drilling Engineering, June 1989, p 128. Sinor A., Warren T.M., Behr S.M., Wells M.R. and Powers J.R., 1992, "Development of an anti-whirl core bit", SPE 24587. Wardlaw H.W.R., 1971, "Optimization of Rotary Drilling Parameters" PHD Dissertation, U. of Texas. 2734315 Warren T.M. and Winters W.J., 1986, "Laboratory Study of Diamond-Bit Hydraulic Lift", SPE Drilling Engineering, aug 1986. Warren T.M., 1987 "Penetration-Rate Performance of Roller-Cone Bits", SPE Drilling Engineering, march 1987. Warren T.M. and Armagost W.K. "Laboratory drilling performance of PDC bits" SPE Drilling Engineering, june 1989. Warren T.M. and Sinor A., "Drag-bit performance modeling", SPE Drilling Engineering, june 1989. Wells R., "Dynamics of rock-chip removal by turbulent jetting" SPE Drilling 1 0 Engineering, june 1989. Zijsling D.H. "Single cutter testing: a key for PDC bit development" SPE 16529 Offshore Europe Aberdeen, 1987. La présente invention sera mieux comprise à la description suivante, illustrée par 1 5 les figures ci-annexées, parmi lesquelles: - Les figures 1A et lB représentent le modèle physique dans les conditions initiales et en cours d'évolution à l'instant t. - La figure 2 représente la courbe d'équilibre obtenue dans le cas d'une application particulière du modèle selon l'invention. Le modèle présenté ci-après est un modèle d'évolution non linéaire avec, dans une première variante, trois variables indépendantes supposées caractériser complètement l'état du système de forage. Il s'agit en fait d'un modèle de taillant dit "local" dont le fonctionnement suffit à décrire, dans cette variante, une moyenne du comportement global de l'outil de forage. Sur la figure lB, le taillant est en interaction avec la roche vierge 2 et la pénétration actuelle S constitue une première variable d'état. La figure 1A montre les conditions initiales 6 2734315 o le taillant de hauteur H, fixé sur un corps 3, a pénétré de la profondeur S, dans la roche. Des études spécifiques sont conduites par ailleurs sur le processus de coupe qui montrent la difficulté de la prise en compte et la diversité des modes de représentation: indépendance plus ou moins garantie des effets de coupe et de butée, lien non nécessairement biunivoque de la pénétration et de l'effort normal, justifié par la théorie de la plasticité, influence de reprises successives (écrouissage). L'hypothèse retenue dans ce travail consiste à se contenter d'un lien biunivoque entre effort normal s'exerçant sur le taillant et enfoncement. Soit Ws l'effort vertical dit 1i 0 "solide" associé à cette pénétration. Le lien de Ws avec 8 sera explicité ci-après. Chacun des NC taillants équivalents constituant l'outil produit des copeaux de roche et cette production instantanée, supposée proportionnelle à S est partiellement évacuée dans l'espace annulaire, partiellement stockée au voisinage immédiat du taillant sous forme 1 5 d'un lit de débris dont l'épaisseur actuelle est la seconde variable d'état de notre formulation, baptisée l; ce lit de débris est supposé tapisser uniformément le front rocheux. L'espace résiduel entre le corps de l'outil et le lit de débris permet l'évacuation des copeaux de roche. Cette évacuation est rendue difficile lorsque l'espace résiduel est restreint; on notera h l'épaisseur de lame fluide, évidemment liée à la hauteur totale H du taillant à l'état neuf par la relation: H=h+6+l+7 o y est la hauteur de lame usée, grandeur lentement évolutive que nous considérerons en fait comme paramètre. L'évacuation est également entravée lorsque la viscosité équivalente de la suspension est augmentée, en raison de l'accroissement de concentration en particules solides. Ces deux effets sont traduits par la relation suivante: N3= D4 WH= h3 7 2734315 On pourra se référer à l'article de Jordaan I.J., Maes M.A.and J.P. Nadreau, 1988, "The crushing and clearing of ice in fast spherical indentation tests", Offshore Mechanics and Arctic Engineering, Houston. La troisième variable d'état s'introduit aussi tout naturellement: il peut s'agir de la concentration c de la suspension mais l'on fera plutôt choix de la viscosité dynamique "équivalente" associée r7 ou de la viscosité cinématique équivalente v (à distinguer de la viscosité v0 du fluide proprement dit). 1 0 Comme il a été annoncé, on définit comme commandes les grandeurs sur lesquelles une intervention est possible ou souhaitable; ce sont essentiellement: * le poids sur outil W * la vitesse de rotation N 1 5 * le débitfluide ou la puissance hydraulique; en fait, dans le présent modèle, la vitesse de l'écoulement en sortie de duse v". Dans le présent exemple, ces grandeurs sont supposées constantes et donc assimilables aux nombreux paramètres du problème. On pourra néanmoins envisager la réponse du système à une perturbation de cette commande, et considérer divers types de régulation associées à la variabilité des propriétés des roches. Dans le présent modèle, l'analyse de la décomposition du poids sur l'outil repose sur le principe de séparation entre une composante classique dite solide Ws justiciable des formules de représentation habituelles, et une portance hydraulique W,, qui s'accroît considérablement lorsque l'épaisseur h de la tranche fluide diminue et la viscosité équivalente 7/ s'accroît; on écrit: W = Ws + Wn 8 2734315 La composante solide Ws est explicitée selon l'article de Kuru E. and Wojtanowsicz A.K., 1988, "A Method for Detecting In-Situ PDC Dull and Lithology Change", IADC/SPE Drilling Conference, Dallas, feb 28, march 2, 1988. A, l'aire de butée de chaque taillant au stade d'usure y Ac(â) l'aire de coupe lorsque l'usure est y, la pénétration solide S Sp et Sc les résistances de la roche, respectivement en compression et au cisaillement Nc nombre de taillants i 0 DB le diamètre de l'outil a et g+ caractéristiques liées à l'interface outil/roche On écrit: WNA =S + A() Sc(sina+,u+ cosa) 1 5 La composante hydraulique est explicitée selon l'article de Jordaan I.J., Maes M.A.and J.P. Nadreau, 1988, "The crushing and clearing of ice in fast spherical indentation tests", Offshore Mechanics and Arctic Engineering, Houston. Wu = N6hD 3 W-h3 11 viscosité (dynamique) équivalente de la suspension boue plus particules solides. La circulation contrariée du fluide de forage (enrichi en particules) et en particulier la perte de charge au front de l'outil sont des indicateurs de cet effet de portance. La présente invention décrit également un modèle de rupture de roche intégré au modèle de foration. Il s'agit d'un modèle de représentation avec un schéma idéalisé de copeau parallélépipédique d'épaisseur 6 et d'aire carrée, de côté mDc, o D, est le diamètre 9 2734315 hydraulique considéré pour l'évacuation. Malgré la simplicité de cette géométrie, on considère qu'il est important de prendre en compte un étalement granulométrique. On considère ainsi une distribution gaussienne de tailles Dc qui tient compte: * de la profondeur actuelle de l'entaille 6 * de la ductilité de la roche exprimée au travers du paramètreu = E(Dc) / 6 * d'une dispersion caractérisée par l'écart-type a. 1 0 E ( D.) exprime la moyenne de la distribution en tailles et /u reflète le degré de ductilité de la roche rompue dans les conditions du forage, caractéristique supposée indépendante de 6; m>l est un paramètre liant diamètre hydraulique et géométrie; on supposera souvent m = 1. 1 5 Plutôt que la variable Do, on pourra préférer introduire le nombre n de copeaux libérés par chacun des N, taillants d'un outil de diamètre DB au cours d'une révolution, de sorte que: 7r DB n = N - 4m2N- []-- Dans l'article "A Dynamic Model for Rotary Rock Drilling", Journal of Energy Resources Technology, june 1982, vol 104 p 108, des auteurs Eronini I. E., Somerton W.H. and Auslander D.M., 1982, on considère, pour un outil tricône, un modèle d'évacuation de copeau qui est reproduit ici avec, cependant, l'introduction d'une granulométrie "étalée". L'expression des efforts hydrodynamiques s'exerçant sur le copeau de roche délimité par la rupture, utilisée dans le présent modèle, est également décrit dans l'article ci- dessus. 2734315 Les fondements du modèle sont les suivants: Pour libérer le copeau, il faut en premier lieu vaincre l'effet de rétention dû à la différence de pression existant entre pression de boue et pression de pore dont l'effet est considérable devant celui de la gravité; l'effort associé est supposé vaincu par l'effet de portance seul FL (L=lift) dont l'expression est présentée en annexe 1. La constante de temps VL du processus est extrêmement courte et donc négligée vis- à-vis de celle associée à l'effet d'entraînement proprement dit (FD et CD:; D = drag). Le copeau est ensuite accéléré depuis la position o il est conceptuellement sorti de son logement sous l'effet de l'effort de traînée FD jusqu'à l'espace annulaire; Soit òe le poids propre du copeau de roche de taille courante De et (o. l'effort de succion s'exerçant sur ce fragment pour le retenir; la condition d'évacuation s'écrit: avec un modèle de représentation de c dû à Eronini (1982), dont le détail n'est pas reproduit ici, condensé grâce au paramètre A, en fonction notamment de la présence d'un cake dont la perméabilité est supposée connue. FL >1>+ AP îDo 6p, g pc masse volumique des particules solides. En pratique, le terme 1 est tout à fait négligeable devant le second. Seules les particules caractérisées par Dc 1 1 2734315 La position de D par rapport à la courbe granulométrique conditionne la proportion de particules "évacuées" par rapport à celles "produites". On supposera que les distributions sont normales; la répartition de tailles Dc e n fonction de a dépend, certes, de la ductilité de la roche mais on a supposé qu'il n'y a pas d'effet d'échelle et que seule la distribution Dc / 5 est à caractériser. On montre que: ò Le seuil de décollement est d'autant plus élevé que l'épaisseur È est plus 1 0 petite. * Le fractionnement en un nombre de copeaux élevé (roche ductile avec /u faible) favorise le décollement et donc les possibilités d'évacuation. * L'augmentation du débit (à travers la vitesse en sortie de duse v") et de la viscosité stimulent évidemment l'évacuation. L'écriture du bilan de masse est effectuée comme suit: Supposant un instant que la granulométrie ne soit pas étalée. On écrit alors: dVR / dt=O si FL w, o r = rD puisque l'accélération du copeau s'effectue principalement sous l'effet des efforts de traînée. Vf est le volume élémentaire du copeau et N, le nombre de sites de production, autrement dit, le nombre de taillants. VR, homogène à un volume par unité de temps, est le débit d'évacuation solide. Le rythme de production solide (volume par unité de temps) doit être supposé égal à: 4 N 12 2734315 ce qui donne un bilan de progression, exprimé cette fois en unité de longueur par unité de temps: s=Nc5- /VR sA- 2 4r D2 si ce bilan est positif, il y a accumulation de débris et enrichissement de la suspension; si le bilan est négatif, les conclusions sont inversées en présence d'un fond enrichi de matières solides; sinon, l'évacuation est parfaite et il n'y a pas lieu de se poser le présent problème. Nous utilisons dans le présent modèle une granulométrie distribuée selon la loi 1 0 normale; plus précisément Dc / 6 est supposée distribuée selon une loi normale de moyenne u et d'écart-type ar. Il en résulte un facteur de minoration Z (calculé en annexe 2) multiplicateur de NoVf I /D fonction, comme il a été dit plus haut, du décalage entre D,', taille de copeau réalisant exactement l'équilibrage, et la distribution. On écrit donc: (Do VR =Zce i 5 ò0 est le poids de copeau dont la taille est D (pour l'épaisseur 6) Le cheminement du calcul est présenté en annexe 1. Il permet d'évaluer successivement, pour le copeau de taille courante DÀ * l'effort de portance FL * l'effort de traînée FD et le temps caractéristique associé rD Le bilan se présente alors sous la forme: + = N6 s- = moZ(8) o mo est défini plus loin. 13 2734315 Le terme d'évacuation ne dépend de 6 que par l'intermédiaire de Z et est conditionné, à technologie fixée, par: * la vitesse vn * la viscosité de la boue * la pression de rétention XP * la densité de boue essentiellement, la densité de copeaux accessoirement. ds Pour éviter la notation dérivée du type d= s, nous notons B(È) le bilan, homogène à une accumulation (longueur) par unité de temps. En fait, on opère dans ce qui suit deux modifications: (i) la première est une modification de pure forme consistant, par souci d'homogénéité, à rendre adimensionnelle en la remplaçant par y, = 6 / 6. On note encore B le bilan adimensionnel, homogène à l'inverse d'un temps, de sorte que: B(y,) = Ny, - mn.x(y,) (ii) la seconde est réalisée pour rendre compte assez correctement du phénomène de bourrage notamment. Elle consiste à reconnaître la dépendance du terme d'expulsion vis-à-vis des variables d'état I et h. Il nous est apparu assez commode, dans un premier temps, de rendre compte du phénomène en faisant dépendre le terme d'expulsion de la seule variable adimensionnelle y3 = 1/l o de sorte que: B(y, y3) = Nyl -ml (Y3)(Y) o la dépendance m (y3) est explicitée en annexe 3. 1 4 2734315 En toute rigueur, le terme d'expulsion dépend également visiblement de l'épaisseur résiduelle actuelle de la lame fluide, soit h, qui est plutôt considérée comme paramètre dans l'annexe 3. En définitive, le bilan matière solide comporte un terme de production B+ correspondant au rythme de destruction de la roche et un terme d'expulsion B-. Pour ce qui est de la dépendance vis-à-vis des variables d'état yI, Y2, y3, on a fait le choix suivant: B+(t) = B+(yl) roche détruite 1 0 B-(t) = B-(yl, y3) roche expulsée B(t) = B+(t) - B-(t) y, (t) = 6(t) / 8o entaille réduite dans la roche vierge y2(t) = Log h(t)/ho viscosité équivalente de la suspension y3(t) = 1(t) / o épaisseur réduite du lit de débris B (y, Y3) = mofv, (Y3)ZX,/ (YI) = mofv(z)X/,,(x) mo "jauge" de l'évacuation, norme du terme d'expulsion X^,/(x) dépendance, dite principale, vis-à-vis de la pénétration (yI); issue de la fonction de répartition de la loi normale f(z) modulation de l'expulsion selon l'épaisseur du lit de débris (Y3) ( 9\9 ad 4 4 N d p v. v. --ai-Il Jrm. = aB a c,)oJ(,;,p)4 : z m8 DB P(RP)4 ad, ac, al coefficients utilisés dans la formulation hydrodynamique et dont les valeurs peuvent être trouvées dans l'article d'Eronini d diamètre de duse; vn vitesse de fluide en sortie de duse DB diamètre de l'outil Pm, Pc masses volumiques de la boue et de la roche respectivement XP effet de rétention par pression différentielle au travers du copeau 1 5 2734315 x et z sont des variables associées respectivement à y, et y3 permettant une écriture explicite (Annexes 2 et 3). Avant réduction à trois variables d'état, le problème comporte a priori cinq variables dont trois de type géométrique: o, I, h respectivement profondeur d'entaille dans la roche vierge, épaisseur de lit de débris et épaisseur de la tranche fluide. (y hauteur usée de lame est une variable d'évolution lente en comparaison de celles qui seront étudiées dans ce problème; elle intervient donc ici comme paramètre); puis deux variables d'état de type concentration de la suspension; c la concentration, v la viscosité cinématique 1 0 "équivalente" associée (à distinguer de la viscosité vo du fluide de forage proprement dit). Les équations d'évolution résultent de l'écriture: * d'une équation de conservation de la somme des épaisseurs des différentes i 5 tranches qui, exprimée sous forme différentielle sur les variables dimensionnelles r, l, h s'énonce: d6 + dl + dh = O * d'une écriture de répartition du bilan de matière B(yl, y3) ou taux d'accumulation ds/dt entre contributions partielles dues à: (i) épaississement du lit de débris (dl) (ii) augmentation de concentration de la suspension (h dc) (iii) diminution d'épaisseur de la lame fluide (c dh) de sorte que: dl + cdh + hdc = ds * de la loi de commande W = cte = Ws + W. 16 2734315 L'expression de Ws et WH mentionnée plus haut permet d'expliciter, toujours sous forme différentielle, cette commande très particulière. On utilise la notation condensée: ws,â = " dW= ScNc (sin a + it+ cos aod et on suppose que cette grandeur est invariante avec &. La relation différentielle s'écrit donc: dW = O soit: d + dl 3dh + WS, d = 0 ô i 7hWH * de relations de comportement simplifiées: On énonce dans ce qui suit deux relations différentielles, dépendant des seuls i 5 paramètres a et b, liant concentration de la suspension c, viscosité équivalente 17 et épaisseur de tranche fluide h. Ce sont: i dh + ah dl =0 q dc - b(1-c) dl = 0 Ecriture des équations d'évolution: On pose: 1 + WS,ô K=- WH 1 +3a Ce facteur sera noté K(yl,Y2) dans la présentation finale. La manipulation des cinq relations conduit à la réduction suivante: -ds = d6{1 + (a + b)(1 - c)hK} 1 7 2734315 Les cinq variables d'état évoluent ainsi selon le schéma, o X désigne pour simplifier le vecteur état et u la commande, très élémentaire ici: 6= f(X,u) = -(1 + ahK)f(X,u) - = -Kf(X,u) r7 h aK f(X, u) h c = _bKf(X,u) 1-c avec: 1 0 f(X,u) -."Xu) 1 + h(1 - c)(a + b)K Nous explicitons d'une part la similitude des trois dernières relations et utilisons les formes adimensionnelles: A= = A =Y, o0 1 5 E = Log 17 E = Y2 r10 L=-o L=y3 h H = -Log- F = Log(1 - c) 18 2734315 Les équations différentielles prennent alors la forme réduite à trois variables indépendantes seulement puisque l'on a manifestement: H(t) = aE(t) F(t) = bE(t) Les équations d'évolution se présentent alors sous la forme très particulière: I = F1(YlY2,Y3) 1 0 Y2 =-K(y1,y2)F1(y,,Y2,Y3) 3= -(1 + afe-Y' K(yl,y2))F (y,,y2,y3) o: B(y1 Y3) i + (a + b) h- e-(a+b)Y2 K(Yl y2) do 1 5 avec: Ws 6e-(1+3a)y2 1+ ' K(yy2)= WH (1 + 3a)y2 K(yl, Y2) caractérise l'aptitude, compte tenu du bilan B, à canaliser les dépôts sur le lit de débris; K est une forme explicite des paramètres. A titre d'illustration, voici un exemple cohérent de valeurs ayant permis la résolution du cas représenté figure 2. Les simulations ont consisté à varier l'entrée ro, profondeur d'entaille initiale en l'absence de lit de débris (représentative du poids sur outil en conditions de dégagement idéales). Le résultat du calcul est 5*, entaille à l'équilibre - une fois le transitoire écoulé - et qui conditionne la vitesse de pénétration stabilisée. L'efficacité de la pénétration peut devenir 19 2734315 nulle, passé un certain seuil en poids, fonction des paramètres du problème (et ceci correspond au seuil de bourrage). On juge du degré d'efficacité de la foration en comparant les effets de portance "solide" et "hydraulique". La forme, ici très particulière, des équations d'évolution conduit à une convergence monotone de 6 vers sa valeur d'équilibre 6* alors que, intuitivement, on s'attend à des fluctuations (voir commentaires en annexe 4). La liste qui suit concerne donc les entrées du modèle nécessaires à l'identification i 0 du cas. Pour en faciliter la lecture, un classement de ces entrées a été réalisé. -Paramètres de commande: 6o pénétration initiale dans la roche vierge (lien avec le poids sur outil WOB) (varié dans la plage 0 à 1,26 mm) 1 5 N vitesse de rotation, supposée invariable (N = 0,7 tour/s) v, vitesse du jet fluide, en sortie de duse (lien avec le débit de boue Q (v" = 50ms-'))-Paramètres liés à l'outil DB diamètre de l'outil. (DB = 0,2 m) d diamètre de duse. (d = 0,01 m) Nc nombre de taillants; autant de "sites" producteurs de copeaux, autant de supports pour la reprise de l'effort vertical. (Nc = 81) - Paramètre lié au taillant: H hauteur efficace de taillant. (H = 2,65 mm) Le paramètre conditionne la répartition initiale H = 5o + ho. 2734315 - Paramètres liés à l'interface taillant/roche Ay aire caractéristique pour la représentation de l'effort vertical, (fonction de l'usure y). (Ay= 1 mm2) A terme proportionnel à la pénétration 5 représentatif de l'effort de coupe. ( Ac, = 5 mm, soit 5 mm2 de variation d'aire par mm de pénétration) ac et A+ angle caractéristique de coupe; coefficient de frottement, le choix: sin ac + g+cosa = 1 a été fait Sc résistance de "coupe" (cisaillement). (Sc = 500 MPa) Sp résistance en "butée" (compression). (Sp = 500 MPa) Paramètres liés au copeau de roche Pc masse volumique du copeau. (Pc = 2500 kg.m'3) - Paramètres liés à la découpe: Iu élancement moyen de copeau illustrant le degré de fragilité de la coupe I élevé, rupture fragile; IX bas, rupture ductile (I = 2) a resserrement granulométrique de la distribution (écart-type) (o = 0,5) - Paramètres liés à l'expulsion: ! coefficient intervenant dans la définition du diamètre hydraulique v 0 21 2734315 vo viscosité cinématique de la boue; à distinguer de la "viscosité équivalente" caractérisant la suspension notamment pour l'effet de portance hydraulique vo = 10-3m2s- 1 (viscosité dynamique r/=1,25 Pa.s) Paramètres constitutifs liant entre eux certains paramètres d'évolution au niveau des lois d'interface a pour le lien entre viscosité équivalente et épaisseur de lame fluide (a = 1) b pour le lien entre viscosité équivalente et concentration de la suspension (b = 1) La courbe représentée figure 2 est donc la traduction du comportement de l'outil de forage en termes d'efficacité pour ce choix particulier des 23 paramètres. On remarquera notamment une séparation en deux zones représentatives chacune d'un type de fonctionnement tout à fait caractéristique. 1 5 Zone 1: les possibilités d'évacuation des déblais par l'hydraulique sont excédentaires; dans ces conditions, seules interviennent pour limiter les performances en termes de vitesses de pénétration les caractéristiques techniques usuelles liant poids sur outil (WOB) et vitesse de pénétration (ROP). Les cas représentatifs de la zone 1 sont caractérisés évidemment par 8o = '*, puisque le lit de débris ne peut se reconstituer de manière durable. Zone 2: les conditions de dégagement deviennent progressivement de plus en plus défavorables vis-à-vis des conditions de production de roche, avec l'augmentation du poids sur l'outil (équivalente à l'augmentation de 6o). La reprise de ce poids se fait de plus en plus sous forme de portance hydraulique WH due à des conditions graduellement plus difficiles d'expulsion du fluide de forage enrichi de particules (pertes de charge croissantes) au détriment de l'effort vertical solide WS affecté au travail efficace de désagrégation de la roche vierge. A titre d'exemple, et pour compléter l'illustration du cas représenté figure 2, l'effort vertical correspondant à une pénétration de 50 = 0,65 mm de chacun des taillants (A), 22 2734315 compte tenu des caractéristiques de la roche est 172 kN; pour une pénétration 5o = 1 mm, le poids sur outil associé est 243 kN, la contribution hydraulique WH à l'équilibre commence à devenir très significative puisqu'elle avoisine 50 kN. Le seuil de bourrage 80SEUIL = 1,26 mm (dans le cas traité) (B) correspond à la condition d'application du poids sur outil WOB = 295 kN, lequel conduit irrémédiablement en 4,5 secondes vers un colmatage complet de l'espace compris entre le corps d'outil et la formation: le bilan massique production/expulsion de roche est devenu si défavorable qu'aucune possibilité "d'équilibre dynamique" (avec â', pénétration non nulle) n'existe. Il est clair que la détermination de la valeur de 80 au point (C) de la figure 2 donne le point de fonctionnement optimal pour les conditions paramétriques données. En effet, le sommet de la courbe en cloche représente la vitesse d'avancement la plus grande, donc le meilleur rendement de l'outil de forage. ANNEXE 1 EXPRESSION DES EFFORTS HYDRODYNAMIQUES S'EXER ANT SUR LE COPEAU ET CONSTANTES DE TEMPS ASSOCIÉES AUX MÉCANISMES CORRESPONDANTS Effort de portance: évaluation de la taille caractéristique seuil DC La formule de base due à Eronini exprime de manière semi-empirique l'effet de portance s'exerçant sur une particule de diamètre hydraulique DC au voisinage d'un outil de diamètre DB lorsque la vitesse de fluide est v" en sortie de duse de diamètre d. FL = aLp(Dc2d2v Dc dvM- 2 L=aLPm.DB, DB V0 23 2734315 Pm est la masse volumique de boue, Do sa viscosité et aL une constante de proportionnalité. Le copeau de taille DC a un poids propre: oe = m2D.26pcg de sorte que: FL _ aL Pm 1 d2V(dv1 2 0) 1 /2D3/2 ( co0 gm2 p D'2DB26, Vo) 1 0 Plutôt que la variable DC on pourra préférer introduire le nombre n de copeaux libérés par chacune des N, lames au cours d'une révolution: Ir fD2 4m2N, DC L'expression précédente est alors modifiée en: 1 5 FL = aL d2v2 p. NV4n"/4 1 dv, 2 F(J)l4TgDB P, m ô1 v FL=1+ AP AP o 6Pcg p;cg Ona: Nl'4n'/4 ( / 4)14 gDBM p RP dv. 2 m3/2 aL d2v2 P. PI v, m aL, d2v Pm P--gk---j 24 2734315 On montre, après quelques calculs, que: Do = a2 dvn 4dvn - 2 4 -1 2 DB m t8 DB)tV)tA Plus la taille seuil D au dessous de laquelle il y a libération du copeau est * importante et plus le débit solide potentiel est élevé. Effort de trainée: estimation du temps caractéristique d'expulsion et du terme d'évacuation L'effet de traînée FD est évalué par Eronini selon une formule analogue à celle décrivant l'effet de portance. Le temps caractéristique associé est tel que: FD ()'J oV = adV ID V ( DB) g est une estimation de la vitesse moyenne de circulation du fluide au front donnée par Eronini. Ainsi: 1 1 d2V2pm 1 1 (dv 4 v-D C gD Pc m: a ô VO) g Ce temps caractéristique est indépendant de la taille de particule Dc. 2734315 Le bilan solide s'écrit finalement, en explicitant la définition du terme d'expulsion: o VIR = XNC Pg rD 9 9 s=N V aD4 4 Nó d) pldv.) s=N6-7 = =Nô -X- aL--1 > La Y 4DB a,, X m DB PC VO (P)4 ANNEXE 2 CALCUL DU FACTEUR DE MINORATION X Remarque préliminaire Soit Fg,c la fonction de la loi normale de moyenne g et d'écart-type a; il est incorrect d'écrire: ( f) L'écriture supposerait une répartition en nombre relatif de copeaux alors que nous recherchons une pondération en volume: g: granulométrie moyenne imposée par le type de découpe y = E(D) : 6épaisseur de coupe actuelle D : taille seuil pour l'hydraulique (indépendante de 5). 26 2734315 Calcul de X Il faut rapporter le volume correspondant à des copeaux de taille distribuée DC au volume de matière dégagé si la granulométrie était assimilable à une distribution de Dirac sur D,, copeaux dont le volume élémentaire est alors 0>0 Pcg Le volume élémentaire du copeau de taille hydrodynamique DC est: o (JDc) 2 pcg OO de sorte que, pour une granulométrie étalée de moyenne g et d'écart type a: _, - -_ 2_e 2 t2dt ( 2 avec: t= Dc DC= D D S S (Il aurait été plus séduisant d'adopter une loi log normale de manière à ne pas avoir à envisager les t négatifs). D'o: 2 DC D ( X5=, D ) _.y 25e2 t2dt 27 2734315 Une intégration par parties donne, après quelques calculs: Dc X(6)- [ avec: _Xe Y 1) + [(x)- x +2) e -i) (I(x) fonction de répartition de la loi normale centrée réduite. (x) dépend uniquement du rapport de la moyenne sur l'écart-type de la a distribution des tailles g est lié à la ductilité de la roche 1 0 ca est lié au resserrement granulométrique. (g est une variable sans dimension reflétant en quelque sorte la fragilité moyenne de la roche sous forme d'un "élancement" du fragment; on peut également établir un lien entre cette caractéristique et l'angle de frottement d'un matériau). -L E(D,) D-o est le seuil, tenant compte essentiellement des conditions de pression, au- dessus duquel les particules n'ont pas à être comptabilisées dans le bilan d'évacuation. Remarques * Si les conditions (viscosité uo et densité Pm de boue, pression différentielle P) n'évoluent pas, le seuil D est invariant. * Si, en outre, la roche ne change pas de caractéristiques g et C, la courbe X(o) est invariante. 28 2734315 ANNEXE 3 FONCTION MODULATRICE SELON L'ÉPAISSEUR DU LIT DE DÉBRIS Il est apparu nécessaire de moduler le terme d'expulsion selon la variable épaisseur du lit de débris pour tenir compte des mécanismes suivants dits A et B; le mécanisme A étant dominant lorsque la lame fluide h est importante, le mécanisme B prévalant lorsque cette lame devient étroite. * le mécanisme A est caractérisé par une diminution de l'effet de rétention du copeau en présence d'un lit de débris qui consacre notamment l'augmentation des performances à l'équilibre en termes de vitesse de pénétration ROP lorsqu'on augmente le poids sur l'outil WOB (3 * le mécanisme B correspond à une expulsion du copeau rendue difficile en raison de l'étroitesse de la lame fluide (h) vis-à-vis de la taille du copeau (5) (È > h). * la transition (8 = h) correspond à une compétition des mécanismes A et B. Formellement, le bilan d'évacuation est écrit: B-(y ',Y3) = mof (Y3)XZ,(YI) 29 2734315 o v est un index dit de sensibilité au bourrage intervenant comme paramètre dans la fonction f, retenue: ft (z) = 1 + Z1ez -V fV(z)=f f (y3) avec ho- =Z o ho et o0 sont les paramètres associés aux valeurs initiales des variables h et 5. ANNEXE 4 LINEARISATION AU VOISINAGE DE L'ÉQUILIBRE RAPPEL ET NOTATIONS Un système dont les variables d'état sont définies par le vecteur x et les paramètres principaux regroupés dans le vecteur X est décrit par des équations d'évolution de la forme: X=F(X,A) Il est souvent examiné au voisinage d'une position d'équilibre Xo, satisfaisant donc: F(Xo,)=O 2734315 On introduit les coordonnées locales autour de Xo de sorte que: X = Xo + et le système différentiel initial se présente alors sous la forme dite linéarisée [si l'on fait abstraction du terme résiduel f (,X)]. = A(A)ó + f(,;) o A(X) = Ai;. x=x et f(ó,2) = o("1f2) contient des termes de degré supérieur ou égal à 2. Il découle des théorèmes de stabilité de Liapunov que: (i) Xo est stable si les valeurs propres de A(X) ont des parties réelles négatives. (ii) Xo est instable si au moins une valeur propre a une partie réelle positive. (iii) le cas critique est celui o la partie réelle de l'une ou plusieurs valeurs propres s'annule alors que les autres valeurs propres gardent leur partie réelle négative. Application: écriture du jacobien A(R)= dFi= 9Xj/x=Xo 31 2734315 Le système différentiel se présente sous la forme assez particulière: y, =FI(yIy2,y3) Y2 = -K(y,,y2)F1(y1,y2,y3) Y, = -[I + a ho K(yl,y2)eY2]F:(y',Y2,Y3) avec: -B(y'Yy3) l. (Yl 1 + (a + b) ho K(y, y2)e-(a+b)y2 0o A l'équilibre: B(y1,y3) =<> Fl(yl,Y2,Y3)= o Le jacobien s'écrit: K- F d aF, dy, oy2 dy3 do F K dF, A(X)= --F-_K - --_ Fl-K-AF' - K AF -d1, Yl Y2 oY3 a31 a32 a33 32 2734315 Tenant compte des relations traduisant l'équilibre, la matrice jacobienne se présente sous la forme simplifiée suivante: avec = 1 B 1 + (a + b) h Ke-(a+b)y J F. dB a,, oB :al -=_ al2 =o FI aB B oy3 a =Ka, a23 = Ka13 a31 = -( + a gKe- a32 = O a33 = -(1 + a h Ke-y' al3 Ceci se traduit au voisinage de l'équilibre par les écritures: 4-d Y + 0(l 3) (12 Y2 =-g,, + O(1ó1) ÉY = -K h1 + a ( Ke-a y+) soit L* = 1 + a- K e7aY2 33 2734315 o les symboles * rappellent que le calcul doit être effectué pour un triplet (Yl,Y2,Y3) satisfaisant la condition de l'équilibre. Alors: =Fi d B _ dB) décrit le comportement linéarisé au voisinage de l'équilibre, l'allure est de type exponentiel (convergent) si: i- --12-">0 F( B -L' a Y3 La condition: dOB dBB B L--= O définit, conjointement à B(y1,y3) = O, la condition de bifurcation caractérisant le bourrage. La résolution complète et explicite des relations de bifurcation peut s'effectuer si l'on tire parti du fait que, lorsque la transition vers le bourrage se produit, K(y1,y2) est proche de 1. On montre alors que, en variables x et z (voir annexes 2 et 3), la condition de bifurcation se présente sous la forme explicite z fonction de x: 34 2734315 Condition d'équilibre: (B = 0) kLy (z),. (1±) e-zZ(x)=- (T> Condition de bifurcation (dB = 0) X2 e 2 (X+8) -2 =1+(1-) (J xZ. (X) v +z ho - x + a o oa avec k=- M60 On remarquera que la condition de bifurcation est indépendante de N. 1 0 REMARQUE Dans un cas moins particulier que celui-ci, o, notamment, la proportionnalité supposée des grandeurs E = Log - H = -Log F = - Log(1 x) implique un rio h'o' comportement linéarisé à une dimension, on aurait: =A 4= (1,2,(23) Z(s) = déterminant(sl - A) avec 3 valeurs propres complexes ou réelles et une convergence vers l'équilibre d'un type différent de celui supposé ici. 2734315 REVENDICATIONS 1) Méthode permettant d'améliorer les performances de forage dans laquelle on met en oeuvre un modèle de foration, caractérisée en ce que ledit modèle prend en compte les effets de destruction (2) d'une roche par au moins un taillant (1) fixé sur un corps d'outil (3) entraîné en rotation et les effets d'évacuation des débris de roche par un fluide, en calculant un bilan matière à partir: - de la production de débris de roche par le taillant ayant pénétré dans la roche d'une profondeur 6, - d'un lit de débris recouvrant ladite roche sous une épaisseur I, - d'une lame de fluide d'épaisseur h comprise entre ledit lit de débris et ledit corps, ladite lame fluide ayant une concentration c en débris, 1 5 - de paramètres de commandes, - de paramètres d'environnement, et en ce que l'on détermine les conditions de forage en fonction de la réponse dudit modèle pour des valeurs déterminées desdits paramètres. 2) Méthode selon la revendication 1, caractérisée en ce que l'un au moins desdits paramètres: poids sur l'outil, vitesse de rotation de l'outil et débit de fluide, est un paramètre de commande. 3) Méthode selon l'une des revendications précédentes, caractérisée en ce que, dans ledit modèle, la portance W de l'outil est décomposée en une composante solide Ws et une composante hydraulique W. fonction notamment de la lame fluide. 36 2734315 4) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on considère un étalement granulométrique des débris distribués selon une loi normale fonction de la profondeur t de l'entaille, de moyenne g lié à la ductilité de la roche et d'une dispersion caractérisée par l'écart-type a. ) Méthode selon l'une des revendications précédentes, caractérisée en ce que ledit bilan matière solide B(t) est tel que B(t)=B+(t)-B (t), dans lequel B+ (t) est un terme de production de débris dépendant de 3 et correspondant au rythme de destruction de la roche, et B-(t) est un terme d'expulsion dépendant de I et h.