FR 3088089 A1 20200508 FR 1860155 A 20181105 DOMAINE DE L'INVENTION L'invention concerne de manière générale le forage d'un tunnel, et plus particulièrement l'optimisation du forage par un tunnelier du type à pression de boue ou à densité variable en fonction de la nature du terrain en front de s taille. ARRIERE-PLAN TECHNOLOGIQUE On connaît des machines de creusement de tunnels, appelées tunneliers, qui comportent une structure mobile de grandes dimensions 10 consistant en une grande usine mobile à l'avant de laquelle est disposé un bouclier ayant une section compatible avec la future section en adéquation avec la forme définitive du tunnel (tunnel de section circulaire, bilobé, rectangulaire...) La partie antérieure du bouclier qui vient en contact avec le front de 15 taille pour réaliser la découpe de la formation géologique traversée par le tunnel comporte une tête de coupe supportant des outils de travail, notamment des molettes et couteaux, entraînée en rotation à une variable vitesse et couplée à une poussée que l'on adapte à la nature du terrain à creuser. 20 Il est connu du document WO 03/087537 un dispositif comprenant un détecteur d'usure, un détecteur de tours et un détecteur de charge placés sur le support de la mollette. On connaît également d'autres publications (EP 1 098 066, US 5 330 292, EP 1 253 287, JP 11 270283, JP 09 228 778) qui placent des capteurs d'acquisition des données sur l'arbre de la molette ou 25 sur l'arbre de la tête de coupe ou sur la tête de coupe. Il a également été proposé dans le document FR 2 874 959 de placer directement sur chaque molette concernée un ou plusieurs capteurs choisis parmi les capteurs aptes à fournir des signaux représentatifs de la force de la pénétration de la molette dans le terrain, de l'état de rotation de la molette, 30 de la position de la molette dans le terrain et de la température de la molette afin de fournir des informations plus précises sur la nature du terrain à creuser. Ce procédé s'avère très efficace dans la mesure où il s'affranchit 2 des déformations de l'arbre ou du support de la molette, de la tête de coupe ou de l'arbre de la tête de coupe, lors de la mesure des différents paramètres de forage. Quels que soient les procédés proposés, le pilote doit systématiquement déduire des valeurs mesurées par les différents capteurs les conditions imposées par le terrain traversé et ajuster différentes caractéristiques du tunnelier (vitesse d'avancement, couple, etc.) en fonction de ces valeurs. Cet ajustement doit cependant tenir compte du type de tunnelier utilisé, et nécessite donc un savoir-faire important de la part du pilote. De plus, la nature du terrain n'est généralement pas homogène, de sorte que les valeurs ainsi obtenues peuvent ne pas être représentatives de la nature réelle du terrain en front de taille. RESUME DE L'INVENTION Un objectif de l'invention est donc de proposer un procédé d'optimisation de caractéristiques d'un tunnelier, notamment d'un tunnelier de type à pression de boue ou à densité variable, qui permette d'optimiser en temps réel des caractéristiques du tunnelier afin d'améliorer le forage sans risquer d'endommager le tunnelier, quelle que soit la nature du terrain en front de taille ou encore le type de tunnelier utilisé. Pour cela, l'invention propose un procédé de détermination d'un modèle d'interactions terrain/machine lors d'un forage d'un terrain par un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD, ledit 25 procédé comprenant les étapes suivantes : - St : obtention d'un ensemble de paramètres caractéristiques d'un forage par au moins un tunnelier donné sur au moins un chantier de forage, notamment d'un tunnelier de type à pression de boue ou VD, - S2 : identification d'un ensemble de formules dépendant de tout ou 30 partie des paramètres de forage, - S3 : détermination, à partir des formules ainsi identifiées, d'un ensemble de variables, 3 - S4 : application d'un algorithme de classification non supervisée aux variables de sorte à obtenir des groupes d'individus identifiés selon un critère prédéfini dudit algorithme, - S5 : application d'un algorithme de classification supervisée aux 5 variables et aux groupes d'individus ainsi déterminés de sorte à obtenir un modèle d'interaction terrain/machine reliant lesdites variables auxdits groupes. Certaines caractéristiques préférées mais non limitatives du procédé 10 de détermination d'un modèle d'interactions décrit ci-dessus sont les suivantes, prises individuellement ou en combinaison : - la classification non supervisée comprend un algorithme de segmentation (clustering), par exemple un algorithme K-MEANS. - la classification supervisée comprend un algorithme des forêts 15 aléatoires (random forest). - les variables dépendent de l'un au moins des paramètres de forage suivants : un couple d'une roue de coupe du tunnelier donné, une vitesse de rotation de la roue de coupe du tunnelier donné, une vitesse d'avancement du tunnelier donné, une force de contact du tunnelier donné, une surface de zo la roue de coupe du tunnelier donné, un rayon de la roue de coupe du tunnelier donné, une pression de confinement à l'axe du tunnelier donné. - l'étape S4 de classification non supervisée détermine entre 8 et 10 groupes d'individus différents. - chaque variable décrit une tranche de terrain, telle qu'un anneau donné posé par l'au moins un tunnelier donné, sur l'au moins un chantier de forage et dans lequel les groupes d'individus déterminés à l'étape S4 correspondent à des groupes de tranches de terrain, notamment des groupes d'anneaux posés par l'au moins un tunnelier donné sur l'au moins un chantier de forage. 4 Selon un deuxième aspect, l'invention propose également un procédé d'optimisation de caractéristiques d'un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD, ledit procédé comprenant les étapes suivantes : s - SO : détermination d'un modèle d'interactions terrain/machine selon l'une des revendications 1 à 5, - 511 : mesure instantanée de l'ensemble des paramètres spécifiques de forage du tunnelier, - S13: détermination du groupe d'individus correspondant aux paramètres de forage mesurés à l'étape S11 à l'aide du modèle d'interactions terrain/machine, - S14: optimisation des caractéristiques du tunnelier en fonction du groupe d'individus ainsi déterminé. 15 Certaines caractéristiques préférées mais non limitatives du procédé d'optimisation décrit ci-dessus sont les suivantes, prises individuellement ou en combinaison : - le procédé d'optimisation comprend en outre, préalablement à l'étape S13, une étape de calcul de tout ou partie des variables déterminées à l'étape 20 S3 à partir des paramètres de forage mesurés à l'étape S11, et dans lequel, au cours de l'étape S13, le groupe d'individus est déterminé en appliquant le modèle d'interactions terrain/machine et aux variables ainsi déterminées et/ou lors du forage, le tunnelier pose successivement une pluralité 25 d'anneaux, et les étapes 511 à S14 sont réitérées lors de la pose de chaque anneau. Selon un troisième aspect, l'invention propose également un produit programme d'ordinateur comprenant des instructions de code pour 30 l'exécution d'un procédé de détermination d'un modèle d'interactions terrain/machine lors d'un forage d'un terrain par un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD comme décrit ci-dessus et/ou 5 d'un procédé d'optimisation de caractéristiques d'un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD comme décrit ci-dessus. Selon un quatrième aspect, l'invention propose un moyen de stockage s lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé de détermination d'un modèle d'interactions terrain/machine lors d'un forage d'un terrain par un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD comme décrit ci-dessus et/ou d'un procédé d'optimisation de caractéristiques d'un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD comme décrit ci-dessus. BREVE DESCRIPTION DES DESSINS D'autres caractéristiques, buts et avantages de la présente invention 15 apparaîtront mieux à la lecture de la description détaillée qui va suivre, et au regard des dessins annexés donnés à titre d'exemples non limitatifs et sur lesquels : La figure 1 est une vue en coupe partielle d'un exemple de réalisation d'un tunnelier conforme à un mode de réalisation de l'invention. 20 La figure 2 illustre de manière schématique un exemple de réalisation d'une unité de traitement adaptée pour exécuter le procédé d'optimisation selon l'invention, ainsi que les paramètres en entrée de ladite unité de traitement. La figure 3 est organigramme illustrant les étapes d'un exemple de 25 réalisation d'un procédé d'optimisation conforme à l'invention. DESCRIPTION DETAILLEE D'UN MODE DE REALISATION Afin d'optimiser des paramètres pour un tunnelier 1, notamment d'un tunnelier 1 de type à pression de boue ou à densité variable (variable density 30 en anglais, VD), lors d'un forage, l'invention propose de déterminer préalablement un modèle d'interactions terrain/machine (procédé SO). L'optimisation des caractéristiques du tunnelier 1 est alors réalisée à partir de 6 ce modèle d'interactions terrain/machine (procédé S10 d'optimisation), et non plus uniquement en fonction du type de terrain théorique en front de taille. Le procédé S10 part en effet du principe selon lequel des terrains de nature différente peuvent interagir de manière similaire avec la partie antérieure 3 5 du tunnelier 1 selon par exemple la profondeur où ces terrains se situent mais aussi selon l'état de la machine. Le procédé SO de détermination d'un modèle d'interactions terrain/machine et le procédé S10 d'optimisation peuvent être mis en oeuvre par une unité de traitement 4 adaptée pour exécuter le ou les procédés S10, par exemple un ordinateur (pouvant être logé dans le tunnelier 1 ou distant et connecté par voie filaire audit tunnelier 1) ou un serveur distant disposant de moyens de traitement. Ainsi, l'unité de traitement 4 peut par exemple comprendre une mémoire 5 dans laquelle sont stockées les instructions de 15 code pour l'exécution du ou des procédés et un calculateur 6 de type processeur, microprocesseur, microcontrôleur, etc., configuré pour exécuter lesdites instructions. Comme nous le verrons par la suite, cette unité de traitement 4 reçoit en entrée des paramètres spécifiques de forage, qui peuvent comprendre des mesures obtenues par un ou plusieurs capteurs C1 - 20 C5 fixés sur le tunnelier 1 et donnant une information sur le forage (pression de confinement à l'axe du tunnelier 1, couple de la roue 2 de coupe du tunnelier 1, vitesse d'avancement du tunnelier 1), ou des caractéristiques T1-T3 qui peuvent être entrées directement dans l'unité de traitement 4 et qui caractérisent le tunnelier 1 en lui-même (rayon de la roue 2 de coupe du 25 tunnelier 1, surface de la roue 2 de coupe du tunnelier 1, etc.) Les paramètres spécifiques sont donc soit issus de capteurs Cl -05 placés sur le tunnelier 1 et capables de mesurer une valeur instantanée du paramètre correspondant, soit liés au tunnelier 1 en lui-même et sont fixes (paramètres T1-T3 sur la figure 2). 30 7 Lors d'une première phase SO, le modèle d'interactions terrain/machine lors d'un forage d'un terrain par un tunnelier 1 est déterminé, de préférence dans le cas d'un tunnelier 1 du type à pression de boue ou VD. 5 A cet effet, au cours d'une première étape S1, un ensemble de paramètres spécifiques de forage par au moins un tunnelier 1 donné sur au moins un chantier de forage est obtenu. Le procédé s'appliquant préférentiellement aux tunneliers du type à pression de boue ou VD, les paramètres de forage sont de préférence 10 obtenus pour de tels tunneliers. Dans une forme de réalisation, de tels paramètres spécifiques sont obtenus pour différents types de chantiers et dont la nature du terrain diffère, afin que le modèle d'interaction terrain/machine soit le plus complet possible et puisse ensuite être appliqué à tout type de chantier. Dans une forme de 15 réalisation, ces paramètres spécifiques peuvent avoir été pré-enregistrés dans une base de données regroupant l'ensemble de ces paramètres spécifiques mesurés pour différents chantiers de forage. Lors de l'étape S1, les paramètres spécifiques sont alors obtenus en interrogeant ladite base de données. 20 Les paramètres de forage peuvent notamment comprendre tout ou partie des paramètres spécifiques suivants : des paramètres machines relevés par des capteurs C1-C4 de l'unité de traitement 4 du tunnelier 1 sur le chantier concerné, des caractéristiques reflétant des profils en long géologique, des mesures obtenues par des capteurs C1 -05 lors d'essais in 25 situ, des mesures obtenues par des capteurs C1-05 lors d'essais laboratoires, des capteurs Cl -05 placés sur des molettes instrumentées du tunnelier 1, des caractéristiques issues de rapports fournis par une station de traitement de boues du tunnelier 1, des caractéristiques obtenues par analyse de risque de creusement, des mesures de tassement en surface obtenus par 30 des capteurs dédiés, des mesures effectuées lors de la maintenance des molettes et couteaux sur la roue 2 de coupe, des mesures déterminant l'état d'usure des molettes et couteaux sur la roue 2 de coupe, des caractéristiques 8 reflétant la configuration du tunnelier 1 (rayon de la roue 2 de coupe, surface de la roue 2 de coupe, etc.). Plus précisément, au cours de l'étape S1, tout ou partie des paramètres suivants peuvent être obtenus, pour chaque chantier et chaque 5 tunnelier 1 examiné : un couple de la roue 2 de coupe du tunnelier 1 examiné, une vitesse de rotation de la roue 2 de coupe du tunnelier 1 examiné, une vitesse d'avancement du tunnelier 1 examiné, une surface de la roue 2 de coupe du tunnelier 1 examiné, un rayon de la roue 2 de coupe du tunnelier 1 examiné, une pression de confinement à l'axe du tunnelier 1 examiné. Le 10 couple est par exemple calculé à partir des mesures obtenues par plusieurs capteurs Cl mesurant le couple, par exemple à l'aide de convertisseurs de fréquence de chacun des moteurs, placés sur chacun des moteurs principaux de l'entrainement principal. La vitesse de rotation de la roue 2 de coupe est par exemple mesurée à l'aide d'un capteur C2 du type codeur angulaire, qui 15 peut être à lecture optique, placé sur le joint tournant de la roue 2 de coupe. La vitesse de d'avancement du tunnelier 1 peut être calculée à partir des mesures obtenues par plusieurs capteurs C3 d'élongation de vérins, par exemple des capteurs C2 de position à câble placés sur les différents vérins du tunnelier, puis dérivée du temps pour obtenir la vitesse. Un tel calcul est 20 réalisé en temps réel par l'unité de traitement 4. La force de contact est par exemple calculée à partir des mesures obtenues par plusieurs capteurs en fonction de la configuration du tunnelier, notamment en fonction de ses articulations et peut par exemple être mesurée à l'aide de capteurs C4 du type capteurs de pression placés sur les vérins hydrauliques de la tête de 25 coupe, sur les vérins des articulations et sur les vérins de poussée des anneaux. La pression de confinement à l'axe est par exemple mesurée à l'aide d'un capteur C5 du type capteur de pression, placé dans la chambre de confinement au plus près de l'axe de rotation de la roue de coupe. Ces paramètres spécifiques sont obtenus pour l'ensemble des 30 anneaux posés par le tunnelier 1 concerné sur le chantier. 9 Le cas échéant, une étape de tri des paramètres de forage peut également être réalisée, afin de ne pas tenir compte des valeurs aberrantes lors de la suite des étapes du procédé SO. Pour cela, les valeurs hors limite des différents capteurs C1-05 utilisés pour obtenir ces paramètres sont 5 écartées. Au cours d'une deuxième étape S2, un ensemble de formules dépendant de tout ou partie des paramètres de forage sont identifiées. Les formules peuvent notamment comprendre ou être établies à partir 10 de formules métiers utilisant tout ou partie des paramètres de forage afin de créer des combinaisons ayant un grand pouvoir statistique et qui sont susceptibles d'être interprétées sur le terrain. Les formules métier utilisées lors de cette étape peuvent être connues et correspondre par exemple à des formules d'énergie spécifiques d'excavation du milieu des forages et de 15 l'adimensionnement désiré afin de pouvoir généraliser les données ainsi obtenues à l'ensemble des chantiers. Le cas échéant, ces formules métier peuvent être ajustées de sorte à augmenter leur pouvoir statistique. Dans une forme de réalisation, chaque formule métier peut être une combinaison de tout ou partie des paramètres de forage suivants : le couple 20 de la roue 2 de coupe, la vitesse de rotation de la roue 2 de coupe, la vitesse d'avancement du tunnelier 1, la surface de la roue 2 de coupe, le rayon de la roue 2 de coupe, la pression de confinement à l'axe. Optionnellement, les formules ainsi identifiées peuvent être adimensionnées en fonction du type de tunnelier 1 utilisé afin que le modèle 25 soit applicable à tout type de tunnelier 1 et que les valeurs soient comparables les unes aux autres. A cet effet, les formules sont chacune divisées par le diamètre du tunnelier 1 utilisé sur le chantier concerné. Dans un mode de réalisation, les formules métier suivantes peuvent par exemple être établies : 30 - l'énergie de rotation SEErot, qui dépend du couple C, de la surface S et de la vitesse de rotation w de la roue 2 de coupe ainsi que de la vitesse d'avancement y du tunnelier 1. 10 - l'indice d'énergie de translation SEEtrans, qui dépend de la vitesse de rotation w et de la surface S de la roue 2 de coupe ainsi que de la force de contact (ou contact force en anglais) F (c'est-à-dire la force de contact entre la partie antérieure 3 du bouclier et le front de taille). 5 - le coefficient de frottement f, qui dépend du couple C et du rayon R de la roue 2 de coupe ainsi que de la contact force F. Par exemple, l'énergie en rotation SEErot (en J/m3) peut être déterminée comme suit : 10 SEErot = (C*0))/(v*S) L'indice d'énergie de translation SEEtrans (en J/m3/rn 1) peut être déterminée comme suit : SE Etrans = F/(v*S) où F est la force contact (en kN) 15 Le coefficient de frottement f (sans dimension) peut être déterminé comme suit : f = 3/2 * C/(F*R) 20 Au cours d'une troisième étape, un ensemble de variables peuvent alors être déterminées à partir des formules ainsi identifiées. Pour cela, les formules sont calculées en utilisant la valeur des paramètres de forage obtenus à l'étape S1. Le cas échéant, ces formules peuvent éventuellement être transformées de sorte à augmenter leur pouvoir statistique et/ou réduire 25 la quantité d'informations traitées. Par exemple, les variables peuvent être obtenues en déterminant, sur un anneau complet, la moyenne des formules appliquées aux paramètres de forage, leur maximum et/ou leur écart-type. En d'autres termes, les variables sont agrégées à l'échelle de l'anneau, par le calcul de la moyenne, du 30 maximum et/ou de l'écart-type pour l'ensemble des formules. Par exemple, dans le cas où les formules identifiées à l'étape S2 comprennent l'énergie en rotation SEErot, l'indice d'énergie en translation 11 SEEtrans et le frottement f, les variables peuvent être calculées en déterminant, pour chaque formule, la moyenne de la formule sur un anneau complet, et en déterminant en outre le maximum du frottement et son écart type pour chaque anneau complet. On obtient donc cinq variables en tout, s correspondant respectivement à la moyenne de l'énergie en rotation sur un anneau, la moyenne de l'indice d'énergie en translation sur un anneau, et la moyenne, le maximum et l'écart-type du frottement par anneau. On notera qu'une analyse en composantes principales (ACP) appliquée à ces cinq variables permet d'obtenir trois axes principaux suivant lesquels sont réparties les variables et présentant une variance de 99%. Cela confirme donc que l'information est résumée de manière effective par ces trois axes, qui permettent à la fois d'obtenir une opposition des variables et des formules qui fournissent une interprétation métier claire. 15 En particulier, le premier axe oppose le coefficient de frottement aux énergies en rotation et en translation. En particulier, lorsque les individus (c'est-à-dire les anneaux sur lesquels les variables ont été calculées) se situent dans les positifs de cet axe, cela implique qu'il y a beaucoup de frottements pour peu d'énergie dépensée. Ce phénomène se traduit au sens 20 métier par l'effet de colmatage de la roue 2 de coupe (par exemple l'effet de colmatage d'un terrain très argileux). A contrario, lorsque les individus se situent dans les négatifs de cet axe, ils ne subissent pas cet effet de colmatage. Ce premier axe reflète donc le colmatage de la machine par le terrain. 25 Le deuxième axe repose principalement sur les énergies en rotation et translation. Si les individus se situent dans les positifs de cet axe, cela implique une dépense d'énergie importante pour avancer et excaver le terrain, ce qui se traduit au sens métier par une dureté du terrain importante face à la roue 2 de coupe (comme c'est le cas pour la roche granitique). A 30 contrario, lorsque les individus se situent dans les négatifs de cet axe, le terrain est plutôt meuble de sorte que peu d'énergie est nécessaire pour avancer. Ce deuxième axe reflète donc la dureté du terrain. 12 Le troisième axe repose sur une friction forte mais stable en association avec l'énergie en rotation et en opposition avec l'énergie en translation. Si les individus se situent dans les positifs de cet axe, cela implique qu'il y a un frottement important et une forte consommation de s l'énergie en rotation, ce qui se traduit au sens métier par de l'hétérogénéité des terrains face à la roue 2 de coupe (comme c'est le cas par exemple pour une phase polylithologique telle qu'un terrain comprenant des lentilles d'argile ou des blocs de roche au sein d'un autre type de terrain - qui sera identifié à l'aide des autres axes). A contrario, lorsque les individus se situent dans les 10 négatifs de cet axe, cela se traduit par des à-coups. Ce troisième axe reflète donc l'homogénéité du terrain. Ces axes présentent donc chacun un sens métier permettant leur interprétation et facilitant ainsi l'application d'un algorithme de classification non supervisée aux variables pour obtenir un ensemble de groupes 15 d'individus (étape S4) et d'un algorithme de classification supervisée aux groupes d'individus et ces variables afin d'obtenir le modèle d'interactions terrain/machine (étape S5). Plus précisément, au cours d'une quatrième étape S4, on détermine 20 donc une classification non supervisée des anneaux décrits par ces variables de sorte à obtenir un ensemble de groupes d'anneaux (ou groupes d'individus) selon un critère prédéfini de l'algorithme. Cette étape de classification non supervisée a pour objectif d'identifier, dans le jeu des cinq variables appliquées à l'ensemble des anneaux 25 complets, un ou plusieurs groupes homogènes qui respectent le critère prédéfini. La classification non supervisée peut notamment comprendre un algorithme de segmentation (clustering), par exemple un algorithme KMEANS. Ceci n'est cependant pas limitatif, d'autres algorithmes de 30 classification non supervisée pouvant être mis en oeuvre sans sortir de la portée de la présente demande, tels que, à titre d'exemples non limitatifs, un algorithme de classification ascendante hiérarchique ou encore un algorithme 13 de DBSCAN (acronyme anglais de Density-Based Spatial Clustering of Applications with Noise, pour partitionnement basé sur la densité). De manière connue en soi, l'algorithme K-MEANS se base sur l'initialisation aléatoire de K points devenant des K-centres de gravités, puis s réalise l'attribution des points aux groupes d'individus par rapport à leur distance (minimale) aux k-centres. L'algorithme recalcule ensuite la position du k-centre de gravité dès qu'un nouveau point est attribué au groupe. Il s'agit donc d'un algorithme itératif qui s'arrête en fonction d'un critère prédéfini pouvant comprendre notamment un nombre d'itérations limite, une 10 stabilisation des k-centres de gravités, une attribution de l'ensemble des points à un groupe, etc. Lors de la mise en oeuvre de cette étape, des calculs de coefficients tels que le coefficient de silhouette peuvent donner un nombre de groupes d'individus optimal, correspondant à des groupes d'individus homogènes et 15 bien distants des uns des autres, en fonction des variables obtenues lors de la troisième étape S3 du procédé SO et des paramètres de forage obtenus à l'étape S1. Ce nombre de groupes d'individus peut, le cas échéant, être modifié afin d'obtenir un degré de précision et d'affinage plus adapté pour l'application métier du procédé SO, et notamment le besoin de détailler zo certaines interactions terrain/machine en vue d'optimiser le forage par le tunnelier 1. Typiquement, pour les paramètres de forage identifiés par la Demanderesse et les cinq variables décrites ci-dessus, l'algorithme KMEANS permet de segmenter les variables en six groupes d'individus 25 optimaux. La Demanderesse a toutefois augmenté le nombre de groupes d'individus à dix afin d'obtenir un affinage plus précis des différents groupes d'individus. Afin de permettre une interprétation métier de ces groupes d'individus, ceux-ci ont été placés sur les trois axes identifiés préalablement par l'ACP 30 (colmatage/dureté/homogénéité). Il ressort que trois des groupes d'individus représentent trois niveaux de dureté qui peuvent correspondre par exemple à du béton, du granite plus ou moins décomposé, galets, etc. Par ailleurs, 14 quatre des groupes d'individus représentent quatre niveaux intermédiaires entre la dureté et le colmatage qui peuvent correspondre par exemple à des terrains de type sable, alluvion, sablo-argileux, granite très décomposé. Enfin, trois des groupes d'individus représentent trois niveaux de colmatage s qui peuvent correspondre par exemple à différents types de terrain argileux. Au cours d'une cinquième étape S5, un algorithme de classification supervisée est appliqué aux variables et aux groupes d'individus déterminés à l'étape S4. lo Cette étape S5 de classification supervisée a pour objectif d'obtenir le modèle d'interaction terrain/machine reliant lesdites variables auxdits groupes d'individus. La classification supervisée peut notamment comprendre un algorithme des forêts aléatoires (random forest). Ceci n'est cependant pas 15 limitatif, d'autres algorithmes de classification supervisée pouvant être mis en oeuvre sans sortir de la portée de la présente demande, tels que, à titre d'exemples non limitatifs, un arbre de décision, une machine à vecteur de support ou encore une méthode des k-plus proches voisins. De manière connue en soi, l'algorithme random forest utilise un zo ensemble d'arbres de décisions auxquels n'est fournie qu'une vision parcellaire du problème (un échantillon pris aléatoirement tant en termes d'individus que de variables), ceci afin de prédire la classification. Chaque arbre va ainsi classer l'individu dans l'un des groupes d'individus déterminés à l'étape S4 en fonction de l'information qu'on lui aura fournie. Le random 25 forest fait ensuite le compte des décisions de chaque arbre et attribue le groupe d'individus qui aura reçu le plus de votes à l'individu (l'anneau). L'algorithme permet ainsi d'obtenir un modèle reliant les variables aux groupes d'individus déterminés à l'étape S4 et traduisant ainsi les interactions terrain/machine possibles lors d'un forage par un tunnelier 1. 30 Grâce au procédé SO, on a ainsi déterminé un modèle d'interaction terrain/machine reliant les variables aux groupes d'individus. Ce modèle est 15 adimensionné et peut donc être appliqué à tout tunnelier 1, quel que soit sa taille ou son type (à pression de boue ou VD), mais également à tout chantier, dans la mesure où l'ensemble des interactions en fonction de la nature des terrains ont de préférence été rencontrées dans les chantiers à partir 5 desquels les paramètres de forage ont été obtenus. On notera que le procédé SO permet de déterminer les interactions terrain/machine, et non le type de terrain rencontré par la partie antérieure 3 du bouclier. Deux terrains de natures différentes peuvent ainsi générer substantiellement la même interaction terrain/machine et appartenir au même lo groupe d'individus. Et c'est bien cette interaction qui est déterminée par le procédé SO et qui est pertinente pour optimiser le forage, et non la nature de chaque terrain. Au cours d'une deuxième phase, correspondant au forage d'un tunnel 15 par un tunnelier 1, notamment de type à pression de boue ou VD, des caractéristiques du tunnelier 1 peuvent alors être optimisées (étape 510) en appliquant le modèle d'interaction terrain/machine ainsi déterminé. On notera que cette optimisation 510 est de préférence réalisée pour chaque anneau posé par le tunnelier 1. 20 Pour cela, au cours d'une première étape 511, les paramètres spécifiques de forage sont obtenus instantanément, pour chaque anneau. Les paramètres spécifiques comprennent les paramètres de forages nécessaires à l'établissement des variables identifiés à l'étape S2, et peuvent 25 être soit entrés directement dans l'unité de traitement 4, soit mesurés instantanément par des capteurs Cl -05 du tunnelier 1. Afin que les données soient comparables, l'étalonnage des différents capteurs C1-05 utilisés est réalisé de manière similaire entre les capteurs utilisés pour l'obtention des paramètres de forage à l'étape S1 et les capteurs 30 utilisés pour la mesure instantanée desdits paramètres de forage à l'étape S11. 16 Au cours d'une deuxième étape S12, les variables qui ont été identifiées à l'étape S3 sont calculées pour chaque anneau. Ainsi, dans l'exemple détaillé ci-dessus, pour chaque anneau, on calcule la moyenne de l'énergie en rotation sur l'anneau, la moyenne de l'indice d'énergie en 5 translation sur l'anneau, et la moyenne, le maximum et l'écart-type du frottement pour l'anneau, soit cinq variables au total. Ceci n'est cependant pas limitatif, les variables pouvant être calculées pour des tranches de terrain, qui peuvent par exemple mesurer quelques centimètres de long, en lieu et place desdits anneaux. Dans ce qui suit, pas 10 souci de simplification, l'invention sera détaillée dans le cas où les variables sont calculées par anneau, bien que cela ne soit pas limitatif. Au cours d'une troisième étape S13, le modèle d'interactions terrain/machine est appliqué aux paramètres de forage de sorte à déterminer 15 le groupe d'individus correspondant aux paramètres de forage instantanés qui ont été mesurés pour cet anneau. Plus précisément, pendant l'étape Si 3, les variables déterminées à l'étape S12 sont calculées pour chaque anneau, le modèle d'interactions terrain/machine étant appliqué auxdites variables. Dans le cas où les variables sont calculées par tranche de terrain, le 20 modèle d'interaction terrain/machine est bien entendu appliqué auxdites tranches de terrain. On en déduit alors le groupe d'individus correspondant aux paramètres de forage instantanés du tunnelier pour cet anneau. Il est alors 25 possible de déterminer si des caractéristiques du tunnelier 1 peuvent être modifiées, compte-tenu du groupe d'individus, afin d'optimiser le forage du tunnel. Les caractéristiques du tunnelier 1 pouvant être modifiées en fonction des groupes d'individus comprennent notamment la vitesse d'avancement du 30 tunnelier 1, la contact force, la vitesse de rotation de la roue 2 de coupe du tunnelier 1 et le couple de la roue 2 de coupe du tunnelier 1. Dans le cas d'espèce, ces caractéristiques correspondent sensiblement aux paramètres 17 spécifiques mesurés par les différents capteurs Cl-05. Il ne s'agit ici que d'un exemple, les caractéristiques optimisées sur le tunnelier 1 pouvant être différentes des paramètres spécifiques utilisés dans le modèle d'interaction terrain/machine. 5 La modification des caractéristiques du tunnelier 1 peut être faite manuellement, par l'opérateur du tunnelier 1. Pour cela, dans une forme de réalisation, des statistiques regroupant, pour chaque groupe d'individus identifié, des informations relatives aux 10 caractéristiques du tunnelier 1, peuvent être réalisées. Ces statistiques peuvent par exemple comprendre une moyenne, un minimum, un maximum et/ou une valeur optimale de chacune de ces caractéristiques du tunnelier 1 au sein de chaque groupe. Par exemple, à chaque groupe d'individus peuvent être associés une moyenne, un minimum, un maximum et/ou une 15 valeur optimale de la vitesse d'avancement du tunnelier 1, de la contact force, de la vitesse de rotation de la roue 2 de coupe du tunnelier 1 et /ou du couple de la roue 2 de coupe du tunnelier 1. Ces statistiques peuvent alors être communiquées en temps réel à l'opérateur afin de l'assister dans sa prise de décision et éventuellement dans 20 la modification de l'une ou de plusieurs de ces caractéristiques du tunnelier 1. En variante, ces informations peuvent être utilisées afin d'automatiser ladite prise de décision. 25 On notera que l'invention concerne également un produit programme d'ordinateur comprenant des instructions de code pour l'exécution du procédé de détermination d'un modèle d'interactions terrain/machine ou du procédé d'optimisation de caractéristiques pour un tunnelier 1 conforme à l'invention, ainsi que des moyens de stockage lisibles par un équipement 30 informatique (par exemple un disque dur du calculateur 6) sur lequel on trouve ce produit programme d'ordinateur. REVENDICATIONS1. Procédé (SO) de détermination d'un modèle d'interactions terrain/machine lors d'un forage d'un terrain par un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD, ledit procédé comprenant les étapes suivantes : - Si : obtention d'un ensemble de paramètres caractéristiques d'un forage par au moins un tunnelier donné sur au moins un chantier de forage, 10 notamment d'un tunnelier de type à pression de boue ou VD, - S2 : identification d'un ensemble de formules dépendant de tout ou partie des paramètres de forage, - S3 : détermination, à partir des formules ainsi identifiées, d'un ensemble de variables, 15 - S4 : application d'un algorithme de classification non supervisée aux variables de sorte à obtenir des groupes d'individus identifiés selon un critère prédéfini dudit algorithme, - S5: application d'un algorithme de classification supervisée aux variables et aux groupes d'individus ainsi déterminés de sorte à obtenir un zo modèle d'interaction terrain/machine reliant lesdites variables auxdits groupes. 2. Procédé (SO) selon la revendication 1, dans lequel la classification non supervisée comprend un algorithme de segmentation (clustering), par 25 exemple un algorithme K-MEANS. 3. Procédé (SO) selon l'une des revendications fou 2, dans lequel la classification supervisée comprend un algorithme des forêts aléatoires (random forest). 30 4. Procédé (SO) selon l'une des revendications 1 à 3, dans lequel les variables dépendent de l'un au moins des paramètres de forage suivants : un 19 couple d'une roue de coupe du tunnelier donné, une vitesse de rotation de la roue de coupe du tunnelier donné, une vitesse d'avancement du tunnelier donné, une force de contact du tunnelier donné, une surface de la roue de coupe du tunnelier donné, un rayon de la roue de coupe du tunnelier donné, 5 une pression de confinement à l'axe du tunnelier donné. 5. Procédé (SO) selon l'une des revendications 1 à 4, dans lequel l'étape S4 de classification non supervisée détermine entre 8 et 10 groupes d'individus différents. 6. Procédé (SO) selon l'une des revendications 1 à 5, dans lequel chaque variable décrit une tranche de terrain, telle qu'un anneau donné posé par l'au moins un tunnelier donné, sur l'au moins un chantier de forage et dans lequel les groupes d'individus déterminés à l'étape S4 correspondent à des groupes de tranches de terrain, notamment des groupes d'anneaux posés par l'au moins un tunnelier donné sur l'au moins un chantier de forage. 7. Procédé (S10) d'optimisation de caractéristiques d'un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD, ledit procédé zo comprenant les étapes suivantes : - SO : détermination d'un modèle d'interactions terrain/machine selon l'une des revendications 1 à 5, - S11 : mesure instantanée de l'ensemble des paramètres spécifiques de forage du tunnelier, 25 - S13: détermination du groupe d'individus correspondant aux paramètres de forage mesurés à l'étape S11 à l'aide du modèle d'interactions terrain/machine, - S14: optimisation des caractéristiques du tunnelier en fonction du groupe d'individus ainsi déterminé. 30 8. Procédé (S10) selon la revendication 7, comprenant en outre, préalablement à l'étape S13, une étape (S12) de calcul de tout ou partie des variables déterminées à l'étape S3 à partir des paramètres de forage mesurés à l'étape S11, et dans lequel, au cours de l'étape S13, le groupe d'individus est déterminé en appliquant le modèle d'interactions terrain/machine et aux variables ainsi déterminées. 5 9. Procédé (S10) selon l'une des revendications 7 ou 8, dans lequel, lors du forage, le tunnelier pose successivement une pluralité d'anneaux, et les étapes 511 à S14 sont réitérées lors de la pose de chaque anneau. 10 10. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé (SO) de détermination d'un modèle d'interactions terrain/machine lors d'un forage d'un terrain par un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD, selon l'une des revendications 1 à 6 et/ou d'un procédé (S10) d'optimisation de 15 caractéristiques d'un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD, selon l'une des revendications 7 à 9. 11. Moyen de stockage lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des instructions de code 20 pour l'exécution d'un procédé (SO) de détermination d'un modèle d'interactions terrain/machine lors d'un forage d'un terrain par un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD, selon l'une des revendications 1 à 6 et/ou d'un procédé (S10) d'optimisation de caractéristiques d'un tunnelier, notamment d'un tunnelier de type à pression de boue ou VD, selon l'une des revendications 7 à 9.