FR 3140434 A1 20240405 FR 2210027 A 20220930 L’invention concerne un procédé de contrôle de géométrie (100) mis en œuvre par ordinateur. Le procédé comprend : - fournir (102) une pluralité de mesures (50) représentant chacune une position mesurée dans un repère en trois dimensions d’un point situé sur une structure à contrôler ; - déterminer (103) un système de coordonnées cylindriques (r, θ, z) présentant un axe vertical (19) parallèle à un axe dudit repère en trois dimensions ; - calculer (105), à partir des mesures (50), une pluralité de points de discrétisation (250) ; - calculer (106), à partir des points de discrétisation (250), des positions estimées (300) de N arêtes verticales (13) de la structure à contrôler ; - calculer (107), à partir des positions estimées (300) des N arêtes verticales (13), des positions estimées (400) de N médianes, chaque médiane étant une médiane d’un pan vertical de la structure à contrôler ; et - comparer (108) des distances entre des médianes diamétralement opposées avec une plage acceptable de diamètre (D). Figure pour l’abrégé : Fig. 4 Procédé de contrôle de géométrie mis en œuvre par ordinateur L'invention se rapporte à un procédé de contrôle de géométrie mis en œuvre par ordinateur pour contrôler la géométrie d’une structure, notamment pour contrôler la géométrie d’une structure porteuse pour une installation de stockage de gaz liquéfié. Plus particulièrement, la structure à contrôler peut comporter une paroi globalement verticale de forme globalement polygonale régulière. Arrière-plan technologique On connaît par le document US 8,550,276 B2 une installation de stockage de gaz liquéfié comprenant une paroi verticale et une paroi de fond, où ladite paroi de fond présente une pluralité de secteurs images les uns des autres par rotation, et où ladite paroi de fond a la forme d'un polygone régulier dont chaque côté correspond à un desdits secteurs. Une telle structure est avantageuse en cela qu’elle permet de réaliser chaque secteur avec des éléments identiques, ce qui diminue le nombre d’éléments différents à employer. En particulier, une grande partie de la paroi de fond est réalisée à l’aide d’éléments rectangulaires de dimensions identiques. Une autre installation de stockage de gaz liquéfié de ce type est connue par le document WO 2011/048300 A1. Dans ce document également, l’installation comprend une paroi verticale et une paroi de fond. La paroi verticale présente une pluralité de pans verticaux. La paroi de fond inclut une pluralité de pièces rectangulaires réparties en secteurs images les uns des autres par rotation, les bords des pièces rectangulaires d'un desdits secteurs étant respectivement parallèles et perpendiculaires à l'un desdits pans verticaux. Toutefois, contrairement au document US 8,550,276 B2, le nombre desdits pans verticaux est le double du nombre desdits secteurs. Le nombre de pans verticaux est par exemple choisi égal à 56. Comme cela décrit dans ce document, le fait de prévoir un nombre de pans verticaux qui est élevé, en particulier le double du nombre des secteurs, permet de limiter la quantité de matière nécessaire pour réaliser la structure porteuse devant recevoir la paroi verticale la paroi de fond, à volume de stockage égal. Dans ces deux documents, la structure porteuse est par exemple réalisée en béton. Résumé Certains aspects de l’invention partent du constat que la structure porteuse présente, en pratique, certains écarts dimensionnels avec la forme de polygone régulier idéalement prévue. De tels écarts dimensionnels peuvent rendre difficile la construction de l’installation. Il convient donc que les écarts dimensionnels n’excèdent pas des tolérances spécifiées à l’avance. Toutefois, si la structure porteuse est de dimensions importantes, vérifier manuellement que la structure porteuse tombe dans les tolérances spécifiées peut prendre un temps considérable, pendant lequel on ne peut pas continuer la construction de l’installation, puisqu’on ne sait pas si la structure porteuse est conforme ou non. Une idée à la base de l’invention consiste à fournir un procédé pour contrôler la géométrie d’une structure, qui est mis en œuvre par ordinateur à partir d’une pluralité de mesures représentant chacune une position mesurée en trois dimensions d’un point situé sur la structure à contrôler. L’invention propose ainsi un procédé de contrôle de géométrie mis en œuvre par ordinateur, comprenant : - spécifier un modèle de structure, par exemple représentant une structure porteuse pour une installation de stockage de gaz liquéfié, le modèle de structure comportant des attributs géométriques et des attributs dimensionnels, les attributs géométriques comportant un nombre N de côtés d’un polygone, N étant un nombre entier supérieur ou égal à 3, les attributs dimensionnels incluant une plage acceptable de diamètre ; - fournir une pluralité de premières mesures représentant chacune une position mesurée dans un repère en trois dimensions d’un point situé sur une structure à contrôler ; - déterminer un système de coordonnées cylindriques présentant un axe vertical parallèle à un axe dudit repère en trois dimensions, ledit axe vertical présentant une position qui minimise une dispersion desdites premières mesures selon une cordonnée radiale dans ledit système de coordonnées cylindriques ; - spécifier des premiers secteurs disjoints de l’espace dans le système de coordonnées cylindriques, lesdits premiers secteurs étant définis par des incréments de coordonnées azimutales et des incréments de coordonnées axiales dans ledit système de coordonnées cylindriques ; - calculer, à partir des premières mesures, une pluralité de premiers points de discrétisation, chaque premier point de discrétisation représentant une position moyenne de la structure à contrôler dans l’un desdits premiers secteurs ; - calculer, à partir des premiers points de discrétisation, des positions estimées de N arêtes verticales de la structure à contrôler ; - calculer, à partir des positions estimées des N arêtes verticales, des positions estimées de N médianes, chaque médiane étant une médiane d’un pan vertical de la structure à contrôler ; et - à partir des positions estimées des N médianes, comparer des distances entre des médianes diamétralement opposées avec ladite plage acceptable de diamètre. Avec le procédé de contrôle de géométrie tel que défini ci-dessus, un contrôle de la géométrie de la structure porteuse peut être réalisé à partir des premières mesures, de façon automatisée, c’est-à-dire sans intervention d’un utilisateur excepté pour spécifier le modèle de la structure porteuse. Le contrôle de géométrie est ainsi aisé à mettre en œuvre même lorsque la structure porteuse est de très grandes dimensions. Le nombre N représente un nombre de côtés d’un polygone devant servir de directrice à une paroi verticale de la structure à contrôler. La paroi verticale, si elle est conforme, sera proche d’une forme idéale composée de N pans verticaux séparés par N arêtes verticales et formant une surface cylindrique polygonale ayant un polygone régulier à N côtés comme directrice. En d’autres termes, la plage de diamètre peut représenter une distance acceptable entre des médianes respectives de deux pans verticaux diamétralement opposés, qui assure une ressemblance suffisante avec la forme idéale. Selon des modes de réalisation, le procédé de contrôle de géométrie peut comporter une ou plusieurs des caractéristiques suivantes. Selon un mode de réalisation, ledit axe dudit repère en trois dimensions est parallèle à la direction du champ de gravité terrestre à l’emplacement de la structure à contrôler. Selon un mode de réalisation, les attributs dimensionnels incluent en outre une plage acceptable de rayon, et le procédé de contrôle de géométrie comprend en outre comparer des distances entre lesdites positions estimées des N médianes et ledit axe vertical avec ladite plage acceptable de rayon. En d’autres termes, la plage acceptable de rayon peut représenter une distance acceptable entre les médianes et ledit axe vertical, qui assure une ressemblance suffisante avec la forme idéale. Selon un mode de réalisation, les attributs dimensionnels incluent en outre une plage acceptable de largeur de pan, et le procédé de contrôle de géométrie comprend en outre comparer des distances entre les positions estimées de deux arêtes verticales voisines avec ladite plage acceptable de largeur de pan. En d’autres termes, la plage acceptable de largeur de pan peut représenter une distance acceptable entre deux arêtes verticales voisines, qui assure une ressemblance suffisante avec la forme idéale. Selon un mode de réalisation, les attributs dimensionnels incluent en outre une plage acceptable de non-verticalité d’arête, et le procédé de contrôle de géométrie comprend en outre comparer des distances, perpendiculairement à l’axe vertical dudit système de coordonnées cylindriques, entre des points d’extrémité de chaque arête verticale avec ladite plage acceptable de non-verticalité d’arête. En d’autres termes, la plage acceptable de non-verticalité d’arête peut représenter une distance acceptable, perpendiculairement à un axe vertical de la structure à contrôler, entre deux points d’extrémité d’une arête verticale donnée, qui assure une ressemblance suffisante avec la forme idéale. Selon un mode de réalisation, les attributs dimensionnels incluent en outre une plage acceptable d’ovalisation, et le procédé de contrôle de géométrie comprend en outre comparer une différence de diamètre entre un plus grand cercle inscrit de la paroi verticale et un plus petit cercle circonscrit de la paroi verticale de la structure à contrôler avec ladite plage acceptable d’ovalisation. En d’autres termes, la plage acceptable de non-verticalité d’arête peut représenter une différence de diamètre acceptable entre un plus grand cercle inscrit de la paroi verticale et un plus petit cercle circonscrit de la paroi verticale, qui assure une ressemblance suffisante avec la forme idéale. Selon un mode de réalisation, les attributs dimensionnels incluent en outre une plage acceptable de déformation locale des pans verticaux, et le procédé de contrôle de géométrie comprend en outre comparer des déformations locales de la structure à contrôler avec ladite plage acceptable de déformation locale des pans verticaux. Selon un mode de réalisation, les attributs dimensionnels incluent en outre une plage acceptable de non-verticalité des pans verticaux, et le procédé de contrôle de géométrie comprend en outre : - sélectionner un sous-ensemble des premiers points de discrétisation, dans lequel sous-ensemble les coordonnées azimutales des premiers points de discrétisation sont égales entre elles ; - déterminer des écarts de coordonnée radiale entre les premiers points de discrétisation du sous-ensemble ; et - comparer les écarts de coordonnée radiale avec ladite plage acceptable de non-verticalité des pans verticaux. Selon un mode de réalisation, les écarts de coordonnée radiale sont des écarts de coordonnée radiale entre un premier point de discrétisation du sous-ensemble qui présente la coordonnée verticale la plus faible et les autres premiers points de discrétisation du sous-ensemble. Selon un mode de réalisation, les attributs géométriques comportent en outre une paroi de fond plane, les attributs dimensionnels incluent en outre une plage acceptable de non-planéité de la paroi de fond, et le procédé de contrôle de géométrie comporte en outre : - fournir une pluralité de deuxièmes mesures représentant chacune une position mesurée dans ledit repère en trois dimensions d’un point situé sur une paroi de fond de la structure à contrôler ; - associer aux deuxièmes mesures un système de coordonnées cartésiennes, ledit système de coordonnées cartésiennes présentant un premier axe, un deuxième axe et un troisième axe perpendiculaires entre eux, le troisième axe coïncidant avec ledit axe vertical dudit système de coordonnées cylindriques ; - spécifier des deuxièmes secteurs disjoints de l’espace dans ledit système de coordonnées cartésiennes, lesdits deuxièmes secteurs étant définis par des incréments de coordonnées selon le premier axe et le deuxième axe dudit système de coordonnées cartésiennes ; - calculer, à partir des deuxièmes mesures, une pluralité de deuxièmes points de discrétisation, chaque deuxième point de discrétisation représentant une position moyenne de la paroi de fond de la structure à contrôler dans l’un desdits deuxièmes secteurs ; - calculer un paramètre de non-planéité à partir des deuxièmes points de discrétisation, et comparer le paramètre de non-planéité avec ladite plage acceptable de non-planéité de la paroi de fond. Selon un mode de réalisation, les attributs dimensionnels incluent en outre une plage acceptable de déformation locale de la paroi de fond, et le procédé de contrôle de géométrie comprend en outre comparer des déformations locales de la paroi de fond de la structure à contrôler avec ladite plage acceptable de déformation locale de la paroi de fond. Selon un mode de réalisation, N est un nombre pair, plus particulièrement N est un nombre pair compris entre 8 et 56, plus particulièrement encore N = 8 ou N = 56. Selon un mode de réalisation, le procédé de contrôle de géométrie comprend en outre une étape de génération consistant à générer une représentation visuelle d’un résultat d’au moins une dite comparaison, ou de toutes lesdites comparaisons ou de certaines desdites comparaisons. Une telle représentation visuelle peut être affichée à un utilisateur ou bien enregistrée en vue d’un affichage ultérieur à un utilisateur, par exemple. Ceci permet à permet à l’utilisateur d’avoir des informations utilisables pour décider d’actions correctives à entreprendre sur la paroi verticale et/ou la paroi de fond de la structure à contrôler. Selon un mode de réalisation, le procédé de contrôle de géométrie comporte en outre une étape d’affichage consistant à afficher ladite représentation visuelle sur un dispositif d’affichage pour un utilisateur. Selon un mode de réalisation, le procédé de contrôle de géométrie comporte en outre l’étape d’attribuer une qualification de conformité à la structure à contrôler, la qualification de conformité étant choisie parmi conforme et non-conforme. Selon un mode de réalisation, la pluralité de premières mesures provient d’une acquisition de mesures au moyen d’un instrument de mesure dans la structure à contrôler. Selon un mode de réalisation, l’instrument de mesure comporte un appareil de télédétection par laser disposé dans un espace interne de la structure à contrôler. Selon un mode de réalisation, l’invention fournit en outre un programme d’ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, conduisent l’ordinateur à mettre en œuvre le procédé de contrôle de géométrie selon l’un quelconque des modes de réalisation décrits ci-dessus. Selon un mode de réalisation, l’invention fournit en outre un support de données lisible par ordinateur sur lequel ledit programme d’ordinateur est enregistré. Selon un mode de réalisation, l’invention fournit en outre un dispositif de contrôle de géométrie comprenant au moins un processeur et au moins une mémoire contenant un programme d’ordinateur, dans lequel ladite au moins une mémoire et le programme d’ordinateur sont configurés, avec ledit au moins un processeur, pour causer l’exécution du procédé de contrôle de géométrie selon l’un quelconque des modes de réalisation décrits ci-dessus par le dispositif de contrôle de géométrie. Selon un mode de réalisation, l’invention fournit en outre un procédé de contrôle de conformité pour contrôler une conformité géométrique d’une structure à contrôler, le procédé de contrôle de conformité comprenant : - une étape d’acquisition consistant à acquérir une pluralité de mesures représentant chacune une position en trois dimensions mesurée d’un point situé sur la structure à contrôler ; et - la mise en œuvre du procédé de contrôle de géométrie selon l’un quelconque des modes de réalisation décrits ci-dessus. Selon un mode de réalisation, l’étape d’acquisition est réalisée en plusieurs sous-étapes, l’appareil de télédétection par laser étant disposé à chaque sous-étape en une localisation différente dans l’espace interne de la structure à contrôler. Selon un mode de réalisation, pendant l’étape d’acquisition, les mesures sont associées à une mesure de la direction du champ de gravité terrestre à l’emplacement de la structure à contrôler. Selon un mode de réalisation, la structure à contrôler est une structure porteuse pour une installation de stockage de gaz liquéfié. Selon un mode de réalisation, la structure à contrôler est en béton. Selon un mode de réalisation, l’installation de stockage de gaz liquéfié est destinée à être installée à terre. Dans ce cas, la structure à contrôler peut être réalisée en béton. Selon un autre mode de réalisation, l’installation de stockage de gaz liquéfié est destinée à être installée à bord d’une structure flottante, telle qu’un navire. Dans ce cas, la structure à contrôler peut être une portion d’une double coque que présente la structure flottante. Dans un mode de réalisation, le gaz liquéfié est du GNL, à savoir un mélange à forte teneur en méthane stocké à une température d’environ -162°C à la pression atmosphérique. D’autres gaz liquéfiés peuvent aussi être envisagés, notamment l’éthane, le propane, le butane ou l’éthylène. Des gaz liquéfiés peuvent aussi être stockés sous pression, par exemple à une pression relative comprise entre 2 et 20 bars, et en particulier à une pression relative voisine de 2 bars. La cuve peut être réalisée selon différentes techniques, notamment sous la forme d’une cuve intégrée à membrane ou d’une cuve autoporteuse. En particulier, la paroi verticale de la cuve peut être réalisée par juxtaposition de de rangées verticales de modules de paroi isolants plans et de rangées verticales de modules de paroi isolants de coin comme décrit par exemple dans la demande internationale n°PCT/EP2022/057845. Brève description des figures L’invention sera mieux comprise, et d’autres buts, détails, caractéristiques et avantages de celle-ci apparaîtront plus clairement au cours de la description suivante de plusieurs modes de réalisation particuliers de l’invention, donnés uniquement à titre illustratif et non limitatif, en référence aux dessins annexés. La représente, en coupe perpendiculairement à son axe vertical, la forme d’une structure porteuse polygonale pour une installation de stockage de gaz liquéfié. La est une illustration de principe de la non-planéité d’un pan porteur vertical de la structure porteuse de la . La est un diagramme-blocs représentant les étapes d’un procédé de contrôle de conformité de la structure porteuse de la . La est une vue partielle et en perspective, depuis l’espace interne de la structure porteuse, illustrant une étape d’acquisition de mesures réalisée dans la structure porteuse. La est un schéma illustrant à titre d’explication un repère cartésien dans lequel sont exprimées les mesures obtenues lors de l’étape d’acquisition. La est un schéma analogue à la , illustrant divers attributs dimensionnels de la structure porteuse. La est un schéma analogue à la , illustrant à titre d’explication un paramètre d’ovalisation de la structure porteuse. La est un schéma illustrant à titre d’explication un paramètre de non-verticalité d’arête de la structure porteuse. La est un schéma illustrant à titre d’explication, ensemble avec la , le principe d’une étape de calcul utilisée pour exprimer les coordonnées des mesures obtenues lors de l’étape d’acquisition dans le repère cylindrique de la . La est un schéma illustrant à titre d’explication, ensemble avec la , le principe d’une étape de calcul utilisée pour exprimer les coordonnées des mesures obtenues lors de l’étape d’acquisition dans le repère cylindrique de la . La est un schéma illustrant à titre d’explication un repère cylindrique utilisé pour exprimer les coordonnées des mesures obtenues lors de l’étape d’acquisition. La est un schéma illustrant à titre d’explication, ensemble avec la et la , le principe d’une discrétisation réalisée sur les mesures de la paroi porteuse verticale. La est un schéma illustrant à titre d’explication, ensemble avec la et la , le principe d’une discrétisation réalisée sur les mesures de la paroi porteuse verticale. La est un schéma illustrant à titre d’explication, ensemble avec la et la , le principe d’une discrétisation réalisée sur les mesures de la paroi porteuse verticale. La est un schéma illustrant à titre d’explication le principe d’un calcul de positions estimées d’arêtes verticales et de médianes des pans verticaux de la paroi porteuse verticale. La est un schéma analogue à la et à la , illustrant à titre d’explication le principe d’un contrôle par le calcul de divers paramètres de la paroi porteuse verticale. La est un schéma illustrant à titre d’explication le principe d’un contrôle par le calcul de déformations locales de la paroi porteuse verticale. La est un schéma illustrant à titre d’explication le principe d’un contrôle par le calcul d’une verticalité de la paroi porteuse verticale. La est un schéma illustrant à titre d’explication le principe d’une discrétisation réalisée sur les mesures de la paroi porteuse de fond visible sur la . La est un graphe illustrant à titre d’explication un paramètre de non-planéité de la paroi porteuse de fond. La est un graphe illustrant à titre d’explication un autre paramètre de non-planéité de la paroi porteuse de fond. La est un schéma illustrant à titre d’explication le principe d’un contrôle par le calcul de déformations locales de la paroi porteuse de fond. La est un diagramme-blocs fonctionnel d’un dispositif de contrôle de géométrie pour mettre en œuvre le procédé de contrôle de conformité de la . Comme on l’a mentionné ci-dessus, l’invention s’intéresse à la réalisation d’une installation de stockage de gaz liquéfié, qui porte la référence 1 dans la description qui va suivre. Selon une variante, l’installation 1 est apte à stocker un gaz liquéfié, en particulier du gaz naturel liquéfié (GNL) à une température d’environ -162°C et à pression atmosphérique ou d’autres gaz liquéfiés. On décrit tout d’abord la structure porteuse 10. La structure porteuse 10 comprend au moins une paroi porteuse qui définit une cavité destinée à recevoir la cuve étanche 20. Dans un mode de réalisation, une paroi porteuse principale 12 présente une géométrie à peu près cylindrique qui entoure la cavité. Une telle paroi porteuse principale 12 peut en outre être fermée par une autre paroi porteuse à au moins une extrémité dans la direction directrice. Dans un mode de réalisation, une telle paroi porteuse principale 12 peut s’étendre entre une paroi porteuse de fond et une paroi porteuse de couvercle. L’installation 1 peut être prévue pour être située à terre. La paroi porteuse principale 12 est alors typiquement verticale, c’est-à-dire située dans un plan parallèle à la direction de l’accélération de pesanteur aux tolérances dimensionnelles près. La structure porteuse 10 est par exemple réalisée en béton. De façon non représentée sur les dessins, la paroi porteuse de fond peut être située au niveau du sol ou éventuellement sous le niveau du sol. De façon non représentée sur les dessins, à l’extrémité de la paroi porteuse principale 12 opposée à la paroi porteuse de fond, la structure porteuse 10 comprend une paroi porteuse de couvercle fermant l’espace interne délimité par la paroi porteuse de fond et la paroi porteuse verticale 12. Cette paroi porteuse de couvercle peut supporter divers équipements utilisables pour acheminer le produit liquide depuis ou vers cet espace interne. La paroi porteuse de fond et/ou la paroi porteuse de couvercle peuvent par exemple être planes. Toutefois, d’autres formes sont possibles pour la paroi porteuse de fond et la paroi porteuse de couvercle, notamment des formes de calotte sphérique. En alternative, l’installation 1 peut être prévue pour être installée à bord d’une structure flottante, telle qu’un navire. Dans ce cas, la structure porteuse 10 est une portion d’une double coque que présente la structure flottante. La paroi porteuse principale 12 peut éventuellement être non verticale, et même présenter une direction directrice perpendiculaire à la direction de l’accélération de pesanteur lorsque la structure flottante est au repos. Dans la suite, on considère plus particulièrement le cas d’une installation 1 située à terre et où la paroi porteuse principale 12 est verticale. On parlera ainsi dans la suite d’une paroi porteuse verticale 12. Il est néanmoins précisé que la description qui suit s’applique à une orientation quelconque de la paroi porteuse principale 12 par rapport à la direction de l’accélération de pesanteur. La est une vue schématique en coupe de la structure porteuse 10, prise perpendiculairement à un axe vertical 9 de la paroi porteuse verticale 12. La paroi porteuse verticale 12 est représentée en trait plein sur la . La paroi porteuse verticale 12 forme une surface cylindrique polygonale et a été typiquement construite par des techniques de génie civil. Ainsi, la paroi porteuse verticale 12 présente des pans verticaux 14 séparés les uns des autres par des arêtes 13. La cuve étanche 20 (en pointillés sur la ) est destinée à être installée dans l’espace interne 11 de la structure porteuse 10. La cuve 20 comporte une paroi périphérique verticale 22 destinée à être en regard de la paroi porteuse verticale 12. De façon non représentée sur les dessins, la cuve 20 comporte en outre une paroi de fond en regard de la paroi porteuse de fond et une paroi de couvercle en regard de la paroi porteuse de couvercle. La paroi périphérique verticale 22 peut être formée de rangées verticales de modules de paroi isolants plans et de rangées verticales de modules de paroi isolants de coin comme décrit par exemple dans la demande internationale n°PCT/EP2022/057845. La position et les orientations des pans verticaux 14 peuvent présenter des déviations dimensionnelles par rapport à une forme prévue de polygone régulier. En outre, comme représenté schématiquement sur la , chaque pan vertical 14 peut présenter des déviations dimensionnelles par rapport à une forme idéale plane 14P. Ces déviations dimensionnelles peuvent par exemple être dues à des tolérances dimensionnelles sur une construction en béton. De tels écarts dimensionnels peuvent rendre difficile la construction de l’installation 1. Il convient donc que les écarts dimensionnels n’excèdent pas des tolérances spécifiées à l’avance. Toutefois, si la structure porteuse 10 est de dimensions importantes, vérifier manuellement que la structure porteuse 10 tombe dans les tolérances spécifiées peut prendre un temps considérable, pendant lequel on ne peut pas continuer la construction de l’installation 1 puisqu’on ne sait pas si la structure porteuse 10 est conforme ou non. On décrit ci-après, en référence aux figures 3 à 16, un procédé 1000 de contrôle de conformité (ci-après désigné par « le procédé 1000 » par commodité) pour contrôler la conformité géométrique de la structure porteuse 1 par rapport à un modèle. Dans une première étape 1001 (cf. ) du procédé 1000, la structure porteuse 1 ayant été construite, on acquiert une pluralité de mesures 50 (cf. et ). Chaque mesure 50 représente une position en trois dimensions mesurée d’un point situé sur la paroi porteuse verticale 12. En principe, chaque mesure 50 peut être acquise par toute technique convenable de télédétection. Toutefois, il est particulièrement avantageux d’acquérir les mesures 50 au moyen d’un appareil de télédétection par laser (lidar) 91 disposé dans l’espace interne 11. Le principe de fonctionnement d’un tel appareil 91 est connu en tant que tel et n’est donc pas décrit en détail ici. Dans une variante d’exécution très simple, l’étape d’acquisition 1001 est réalisée en une seule fois, l’appareil 91 ayant été disposé en un seul point de l’espace interne 11. Toutefois, il est préférable que l’étape d’acquisition 1001 soit réalisée en plusieurs sous-étapes, l’appareil 91 étant disposé dans en une localisation différente dans l’espace interne 11 à chaque sous-étape. Ceci peut permettre d’obtenir une acquisition avec une résolution satisfaisante sur l’ensemble de la paroi porteuse verticale 12. Bien entendu, dans ce cas, les mesures 50 obtenues à chaque sous-étape sont post-traitées de telle sorte qu’elles soient toutes exprimées dans un même repère de l’espace. En tout état de cause, à l’issue de l’étape d’acquisition 1001, on dispose d’un grand nombre de mesures 50, exprimées dans un même repère de l’espace. Selon un mode de réalisation, comme représenté schématiquement sur la , les mesures 50 sont exprimées dans un même repère orthogonal cartésien (x 1 , y 1 , z 1 ), dont l’origine A 1 est fixée arbitrairement. La troisième coordonnée z 1 est de préférence exprimée le long d’un axe correspondant précisément à la direction G (cf. et ) du champ de gravité terrestre local, c’est-à-dire du champ de gravité terrestre à l’emplacement de la structure porteuse 10. La direction G du champ de gravité terrestre local peut avoir été acquise pendant l’étape 1001 au moyen d’un capteur 92 approprié, par exemple intégré à l’appareil 91 (cf. ). On précise que bien que les figures 4 et 5 représentent à titre d’explication quelques mesures 50, le résultat de l’étape d’acquisition 1001 peut comprendre un très grand nombre de mesures 50, de l’ordre de 10 7 à 10 9 mesures 50 dans le cas d’une structure porteuse 10 présentant des dimensions internes de l’ordre de 10 1 mètres. Le nombre de mesures 50 dépend, comme cela est connu, des paramètres de balayage de l’appareil 91 et du nombre de sous-étapes de l’étape d’acquisition 1001. Après l’étape d’acquisition 1001, on met en œuvre un procédé de contrôle de géométrie 100 (ci-après désigné par « le procédé 100 »), dont les étapes sont décrites ci-après. Selon un mode de réalisation, le procédé de contrôle de géométrie 100 est mis en œuvre par un ordinateur exécutant un programme d’ordinateur approprié. Toutefois, le procédé 100 peut être mis en œuvre par toute combinaison appropriée de matériel et de logiciel, y compris dans un environnement informatique distribué. Le procédé 100 comprend une première étape 101 consistant à spécifier un modèle auquel on va comparer la structure porteuse 10. Dans le modèle, il est spécifié le nombre N des côtés du polygone qui doit servir de directrice à la paroi porteuse verticale 12. N est un nombre entier supérieur ou égal à 3. Par exemple, N est un nombre pair. Plus particulièrement, N est par exemple un nombre pair compris entre 8 et 56. Dans un mode de réalisation particulier illustré sur les dessins, N est égal à 8. Dans un autre mode de réalisation particulier, N est égal à 56. Dans ce qui suit, le rayon Q, le diamètre D, la largeur de pan W, l’ovalisation OV et la non-verticalité d’arête NV désignent des grandeurs physiques de la structure porteuse 10. Le modèle inclut des attributs dimensionnels, dont notamment des plages de valeurs acceptables pour chacune de ces grandeurs. Les figures 6 à 8 illustrent ces grandeurs dans la structure porteuse 10. La illustre à titre d’explication le rayon Q et le diamètre D. Comme la , la est une vue en coupe de la structure porteuse 10, prise perpendiculairement à un axe vertical de la paroi porteuse verticale 12. L’axe vertical porte la référence 9 sur la . À chaque pan vertical 14, on peut associer un point médian 14A situé dans le plan de coupe de la . Le rayon Q est la distance entre le point médian 14A et l’axe vertical 9 dans ce plan de coupe. Le diamètre D est la distance entre les points médians 14A de deux pans verticaux 14 diamétralement opposés. La illustre aussi à titre d’explication la largeur de pan W. Un pan vertical 14 étant considéré, la largeur de pan W associée à ce pan vertical est la distance, dans le plan de coupe de la , entre les deux arêtes verticales 13 délimitant le pan 14. La illustre à titre d’explication l’ovalisation OV. Comme la et la , la est une vue en coupe de la structure porteuse 10, prise perpendiculairement à l’axe vertical 9. Compte tenu d’un contour réel 12R de la paroi porteuse verticale 12 dans le plan de coupe de la , on peut définir un plus grand cercle inscrit IC dans ce contour réel 12R et un plus petit cercle circonscrit CC du contour réel 12R. Le plus grand cercle inscrit IC présente un diamètre noté D_inner et le plus petit cercle circonscrit CC présente un diamètre D_outer. L’ovalisation OV est la différence entre ces deux diamètres, c’est-à-dire que OV = D_outer – D_inner. La illustre à titre d’explication la non-verticalité d’arête NV. La est une vue schématique de face d’une arête 13 séparant deux pans verticaux 14 (non représentés sur la ). Étant donné deux points d’extrémité 13-1 et 13-2 de l’arête 13, la non-verticalité d’arête NV est la valeur absolue d’une distance, perpendiculairement à l’axe vertical 9, entre les points d’extrémité 13-1 et 13-2. On précise que les points d’extrémité 13-1 et 13-2 peuvent être des points de l’arête 13 dans deux plans de référence, espacés parallèlement à l’axe vertical 9, dans la paroi porteuse verticale 12. La plage acceptable pour chacune des grandeurs mentionnées ci-dessus peut être spécifiée sous la forme d’un intervalle de grandeurs numériques. En variante, une plage acceptable peut être spécifiée sous la forme d’une valeur de référence V et d’une tolérance T sur cette valeur de référence V ; la plage acceptable est alors l’intervalle compris entre V-T et V+T, cet intervalle pouvant être ouvert, fermé ou semi-fermé, selon les besoins de l’implémentation. En ce qui concerne le diamètre D, le rayon Q, la largeur de pan W, la valeur de référence est une valeur non nulle spécifiée à l’avance compte tenu des dimensions et proportions souhaitées de la structure porteuse 10, et les tolérances sont des valeurs non nulles spécifiées à l’avance compte tenu des tolérances acceptables sur ces dimensions et proportions souhaitées de la structure porteuse 10. En ce qui concerne l’ovalisation OV et la non-verticalité d’arête NV, la valeur de référence est nulle, et la tolérance est une valeur non nulle spécifiée à l’avance. Les attributs géométriques et/ou les attributs dimensionnels du modèle peuvent être spécifiées par un utilisateur, par exemple au moyen d’une interface utilisateur implémentée sur ordinateur. En variante, le modèle peut aussi avoir été spécifié à l’avance sous forme d’un fichier de paramètres, éventuellement éditable par l’utilisateur. Le procédé 100 comprend en outre une étape 102 dans laquelle on fournit les mesures 50 décrites ci-dessus, réalisées sur la structure porteuse 10 à contrôler. Concrètement, les mesures 50 sont fournies sous la forme d’un fichier informatique, sur un support de données ou via une interface réseau. Par exemple le fichier informatique est fourni sous le format de fichier .pts bien connu dans le domaine des appareils de télédétection par laser. Le procédé 100 comprend en outre une étape 103 dans laquelle on détermine, à partir des mesures 50, un système de coordonnées cylindriques (r, θ, z). Comme cela est connu, définir un système de coordonnées cylindriques suppose de définir un axe vertical le long duquel on mesure une coordonnée z et des coordonnées polaires r, θ perpendiculairement à cet axe vertical. Ainsi, dans la mise en œuvre de l’étape 103, il convient d’abord de définir un axe vertical le long duquel on mesure la coordonnée z. Dans ce but, le repère orthogonal cartésien (x 1 , y 1 , z 1 ) mentionné ci-dessus dans lequel sont exprimées les mesures 50 est défini de telle sorte que son axe vertical Z 1 soit parallèle à la direction G du champ de gravité terrestre local. L’axe vertical 19 du système de coordonnées cylindrique est choisi parallèle à l’axe vertical Z 1 . On comprend bien que l’axe vertical 19 est ainsi parallèle ou quasiment parallèle à l’axe vertical 9 de la structure porteuse 10 à contrôler. En outre, la position de l’axe vertical 19 présente une position qui minimise une dispersion selon la coordonnée radiale r des mesures 50. On comprend bien que l’axe vertical 19 est ainsi très proche de l’axe vertical 9 de la structure porteuse 10 à contrôler. En outre, comme on va le détailler ci-après, les arêtes 13 sont ainsi précisément les endroits où existent des maxima locaux de la coordonnée radiale r. Diverses méthodes mathématiques sont appropriées pour trouver une position de l’axe vertical 19 qui minimise une dispersion selon la coordonnée radiale r des mesures 50. Dans un exemple simple d’exécution, comme représenté schématiquement sur la , on sélectionne plusieurs (ici, quatre) sous-ensembles T 1 , T 2 , T 3 , T 4 de mesures 50, les sous-ensembles T 1 , T 2 , T 3 , T 4 étant définis par des incréments Δz 1 de coordonnées verticales z 1 ; et on calcule pour chaque sous-ensemble T 1 , T 2 , T 3 , T 4 un cercle d’interpolation C 1 , C 2 , C 3 , C 4 (cf. ) associé et perpendiculaire à l’axe vertical Z 1 , puis on calcule un barycentre C des centres des cercles C 1 , C 2 , C 3 , C 4 . L’axe vertical 19 est défini comme passant par ce barycentre C et étant parallèle à l’axe vertical Z 1 . Ensuite, les coordonnées cartésiennes des mesures 50 dans le repère cartésien (x 1 , y 1 , z 1 ) (cf. ) sont converties par le calcul en coordonnées cylindriques (r, θ, z) (cf. ) selon les relations mathématiques connues. Le choix de la direction dans laquelle θ = 0° peut être ou bien fait arbitrairement, ou bien être fait en fonction de la position d’un repère tracé sur l’un des pans 14 ou d'un élément structurel préalablement choisi sur l’un des pans 14. Le sens positif de l’axe vertical 19 est choisi de préférence contraire au sens de l’accélération de la gravité terrestre. Le procédé 100 comprend en outre une étape 104 dans laquelle on spécifie des secteurs disjoints 200 (cf. ) de l’espace dans le système de coordonnées cylindriques (r, θ, z). Concrètement, ces secteurs 200 sont définis par un incrément δθ de coordonnées θ (cf. et ) et par un incrément δz de coordonnées z (cf. et ). Comme les grandeurs numériques du modèle décrit ci-dessus en rapport à l’étape 101, les incréments δθ et δz peuvent être spécifiés par un utilisateur, par exemple au moyen d’une interface utilisateur implémentée sur ordinateur. En variante, le modèle peut aussi avoir été spécifié à l’avance sous forme d’un fichier de paramètres, éventuellement éditable par l’utilisateur. Les incréments δθ et δz peuvent être spécifiés a priori grâce à une connaissance de la valeur de référence Q 0 du rayon Q décrit ci-dessus. Par exemple les incréments δθ et δz peuvent être spécifiés de telle sorte qu’une superficie de la paroi porteuse verticale 12 incluse dans un secteur 200 donné soit comprise entre 1 cm 2 et 10 cm 2 . Le procédé 100 comprend en outre une étape 105 dans laquelle, à partir des mesures 50, on calcule des points de discrétisation 250, dont un est représenté à titre d’illustration sur la . Plus concrètement, il est calculé un point de discrétisation 250 par secteur 200. À chaque secteur 200, il est associé une plage de coordonnées cylindriques {Θ ± δθ/2, Z ± δz/2}. Le point de discrétisation 250 est le centre du secteur 200. Ses coordonnées cylindriques sont {R, Θ, Z}, où R est la moyenne des coordonnées r de toutes les mesures 50 se situant dans le secteur 200. Le procédé 100 comprend en outre une étape 106 dans laquelle, à partir des points de discrétisation 250 calculés à l’étape 105, on calcule des positions estimées 300 des N arêtes verticales 13. Ce calcul peut être effectué d’un grand nombre de manières. Avantageusement, ce calcul consiste à rechercher, grâce à la connaissance de la valeur de référence W 0 de la largeur de pan W décrite ci-dessus, un ensemble de N droites verticales 13 approximativement espacées de W 0 et dont les positions correspondent à N maxima locaux de r. On comprendra mieux le principe d’un tel calcul en se référant à la . Sur cette figure, on a représenté très schématiquement, dans le système de coordonnées cylindriques (r, θ, z), quelques-uns des points de discrétisation 250 calculés à l’étape 105, ensemble avec des positions estimées 300 des N arêtes verticales 13. Le calcul consiste à rechercher une position d’un ensemble de N droites verticales 300 sous contrainte des deux critères suivants : - (i) chacune des N droites 300 est suffisamment voisine de maxima locaux de r en fonction de θ ; - (ii) chaque droite 300 est séparée d’une droite 300 voisine d’une distance (perpendiculairement à l’axe vertical 19 du système de coordonnées cylindriques (r, θ, z)) qui est incluse dans l’intervalle [W 0 - W 2 ; W 0 + W 2 ], où W 2 est une valeur non nulle. On notera que W 2 n’est pas nécessairement égale à la valeur de la tolérance sur W 0 ; au contraire, W 2 peut par exemple être choisie égale à 2,5 fois ou 3,0 fois cette tolérance. Le critère (ii) peut alors permettre d’empêcher qu’une dépression locale sur l’un des pans verticaux 14 soit considérée à tort comme une arête verticale 13. On notera également que les droites 300 ne passent pas nécessairement par des points de discrétisation 250, et peuvent même ne passer par aucun point de discrétisation 250. En tout état de cause, à l’issue de l’étape 106, on dispose des positions estimées 300 des arêtes verticales 13. Le procédé 100 comprend en outre une étape 107 dans laquelle, à partir des positions estimées 300 des arêtes verticales 13, on calcule des positions estimées de N médianes 500 des pans verticaux 14 de la structure porteuse 10 à contrôler. Dans un exemple d’exécution, toujours en référence à la , cette étape 107 comprend calculer les positions de points 400 qui sont équidistants, perpendiculairement à l’axe vertical 19, des droites 300. De préférence, afin de simplifier le calcul, les positions des points 400 sont calculées à la même coordonnée z que des points de discrétisation 250. S’il s’avère qu’à cette coordonnée z, un point 400 ne coïncide pas avec un point de discrétisation 250, la coordonnée r de ce point 400 est calculée par interpolation (par exemple par interpolation linéaire) à partir des coordonnées de points de discrétisation 250 voisins. Les positions des médianes 500 sont ensuite calculées par interpolation, à partir des positions des points 400. Le procédé 100 comprend en outre une étape 108 dans laquelle, à partir des positions estimées des médianes calculées à l’étape 107, on compare des distances entre des médianes diamétralement opposées avec la plage acceptable de diamètre D. En référence à la , ce calcul peut consister à calculer une distance D c entre deux points 400 diamétralement opposés situés dans un même plan vertical (autrement dit situés à la même coordonnée z), et à contrôler que cette distance D c tombe dans la plage acceptable de diamètre D. Certaines ou tous des comparaisons énumérées ci-après peuvent également être effectuées lors de l'étape 108. Notamment, en référence à la , à partir des positions estimées des médianes 500 calculées à l’étape 107, et plus spécifiquement des positions des points 400, on compare les distances Q c entre l’axe vertical 19 et chacun des N points 400 avec la plage acceptable de rayon Q. En outre, toujours en référence à la , on calcule les distances entre les positions estimées 300 des arêtes verticales 13 calculées à l’étape 106, et on compare ces distances avec la plage acceptable de largeur de pan W. En outre, en revenant à la , on calcule les positions et les diamètres D_inner et D_outer du plus grand cercle inscrit IC et du plus petit cercle circonscrit CC, on calcule la valeur de l’ovalisation OV = D_outer – D_inner, et on compare cette valeur calculée avec la plage acceptable d’ovalisation OV. En outre, à partir des droites 300 calculées à l’étape 106, on calcule les positions de deux points d’extrémité situés sur les droites 300, on calcule la valeur absolue de la distance NV, perpendiculairement à l’axe vertical 19, entre ces deux points d’extrémité, et on compare NV avec la plage acceptable de non-verticalité d’arête. En outre, on compare des déformations locales des pans porteurs verticaux 14 avec une plage acceptable de déformation locale des pans porteurs verticaux 14, qui a été préalablement spécifiée à l’étape 101. Cette comparaison va être expliquée plus en détail en référence à la . Pour un secteur 200 donné, à partir de toutes les mesures 50 incluses dans ce secteur 200, il est calculé un plan d’interpolation locale 600. Ensuite, il est associé à ce plan 600 un repère cartésien local (a, b, c) qui est orthogonal et défini de la façon suivante : l’axe (c) est perpendiculaire au plan 600 ; l’axe (b) est parallèle à l’axe vertical 19 du système de coordonnées cylindriques (r, θ, z) ; et l’axe (a) est orthogonal aux axes (b) et (c). Ensuite, les coordonnées des mesures 50 sont exprimées dans ce repère cartésien local (a, b, c) selon les relations mathématiques connues. Ensuite, une valeur σa ayant été spécifiée, on recherche dans le secteur 200 les mesures 50 dont les coordonnées (a) sont incluses dans un intervalle de largeur σa (cf. ) ; parmi ces mesures 50, on recherche les deux mesures 50 les plus éloignées du plan 600, perpendiculairement à l’axe (c), l’une dans la direction positive de l’axe (c) et l’autre dans la direction négative de l’axe (c) ; on calcule la différence σc entre les coordonnées suivant l’axe (c) de ces deux mesures 50 ; et on compare cette différence σc avec la plage acceptable de déformation locale. Ce processus est répété jusqu’à épuisement du secteur 200, en décalant à chaque fois l’intervalle de largeur σa d’un incrément fixe ; autrement dit, l’intervalle de largeur σa est utilisé comme fenêtre glissante. Les opérations qui précèdent sont répétées pour chacun des secteurs 200. En outre, on calcule un écart de coordonnée radiale entre tous les points de discrétisation 250 ayant une même coordonnée azimutale Θ et le plus bas d'entre eux et on compare cet écart avec une plage acceptable de non-verticalité des pans porteurs verticaux 14, qui a été préalablement spécifiée à l’étape 101. Cette comparaison va être expliquée plus en détail en référence à la . Parmi tous les points de discrétisation 250 ayant une coordonnée Θ donnée, on obtient la coordonnée radiale Rb du point de discrétisation 250b qui présente la plus petite coordonnée Zb suivant l’axe vertical 19. Ensuite, pour chaque point de discrétisation 250 ayant la même coordonnée Θ, on calcule l’écart de coordonnée radiale ΔR = R – Rb, et on compare ΔR avec la plage acceptable de non-verticalité des pans porteurs verticaux 14. Les opérations qui précèdent sont répétées pour chaque coordonnée Θ entre 0 et 360°. En variante, l’écart ΔR peut n’être calculé que pour certains points de discrétisation 250 ayant une coordonnée Θ donnée. En variante, l’écart de coordonnée radiale ΔR peut être calculé avec la coordonnée radiale d’un autre point de discrétisation 250 que le point de discrétisation 250b. Le procédé 100 peut optionnellement comprendre une étape 109 de génération d’au moins une représentation visuelle de la conformité de la structure porteuse des résultats des contrôles effectués à l’étape 108. Par exemple cette représentation visuelle peut comprendre ou consister en des graphes indiquant par des couleurs les régions de la paroi porteuse verticale 12 où les divers paramètres calculés à l’étape 108 ne tombent pas dans les plages acceptables spécifiées à l’étape 101. Pour cela, plus concrètement, pour chaque secteur 200 de la paroi porteuse verticale 12 et pour chaque paramètre calculé à l’étape 108, une variable booléenne est calculée, et prend la valeur NON (pour non-conforme) si le paramètre ne tombe pas dans la plage acceptable associée spécifiée à l’étape 101, et OUI (pour conforme) sinon. En complément ou en alternative, la représentation visuelle peut comprendre ou consister en des graphes indiquant par des couleurs les divers paramètres calculés à l’étape 108 sur une représentation graphique de la paroi porteuse verticale 12. Dans une étape 1010 du procédé 1000, la ou les représentations visuelles générées à l’étape 109 peuvent être affichées à un utilisateur, ce qui permet à celui-ci d’avoir des informations utilisables pour décider d’actions correctives à entreprendre sur la paroi porteuse verticale 12. Le procédé peut aussi aboutir à générer une qualification de conformité choisie parmi conforme et non-conforme, par exemple une variable booléenne. Dans un exemple, la variable booléenne peut prendre la valeur OUI (pour conforme) si un pourcentage de secteurs 200 non-conformes de la paroi porteuse verticale 12 est inférieur à un pourcentage prédéterminé, et prendre la valeur NON (pour non-conforme) sinon. Jusqu’ici, on n’a décrit qu’un contrôle de géométrie de la paroi porteuse verticale 12. Toutefois, comme mentionné précédemment, la structure porteuse 10 peut également comporter une paroi porteuse de fond. On décrit ci-après une variante du procédé 100 qui est applicable dans le cas où cette paroi porteuse de fond 15 (cf. ) est plane. Dans cette variante, à l’étape d’acquisition 1001, l’appareil 91 acquiert en outre pluralité de mesures 51 (cf. ), où chaque mesure 51 représente une position en trois dimensions mesurée d’un point situé sur la paroi porteuse de fond 15 de la structure porteuse 10 à contrôler. À l’étape 102, ces mesures 51 sont fournies en même temps que les mesures 50. Les mesures 51 sont exprimées dans le même repère orthogonal cartésien (x 1 , y 1 , z 1 ) que les mesures 50. À l’étape 103, les coordonnées cartésiennes des mesures 51 dans le repère cartésien (x 1 , y 1 , z 1 ) (cf. ) sont converties par le calcul, selon les relations mathématiques connues, en coordonnées cartésiennes dans un système de coordonnées cartésiennes (x, y, z) dont l’axe vertical z coïncide avec l’axe vertical 19 du système de coordonnées cylindriques pour les mesures 50. À l’étape 104, on spécifie en outre des secteurs disjoints 201 de l’espace dans le système de coordonnées cartésiennes (x, y, z). Concrètement, ces secteurs 201 sont définis par un incrément δx de coordonnées x et par un incrément δy de coordonnées y (cf. ). Comme les grandeurs numériques du modèle décrit ci-dessus en rapport à l’étape 101, les incréments δx et δy peuvent être spécifiés par un utilisateur, par exemple au moyen d’une interface utilisateur implémentée sur ordinateur. En variante, le modèle peut aussi avoir été spécifié à l’avance sous forme d’un fichier de paramètres, éventuellement éditable par l’utilisateur. Les incréments δx et δy peuvent être spécifiés a priori grâce à une connaissance de la valeur de référence Q 0 du rayon Q décrit ci-dessus. Par exemple les incréments δx et δy peuvent être spécifiés de telle sorte qu’une superficie de la paroi porteuse de fond 15 incluse dans un secteur 201 donné soit comprise entre 5 cm 2 et 50 cm 2 .À l’étape 105, toujours en référence à la , il est calculé un point de discrétisation 251 par secteur 201. À chaque secteur 201, il est associé une plage de coordonnées cartésiennes {X ± δx/2, Y ± δy/2}. Le point de discrétisation 251 est le centre du secteur 201. Ses coordonnées cartésiennes sont {X, Y, Z}, où Z est la moyenne des coordonnées z de toutes les mesures 51 se situant dans le secteur 201. À l’étape 108, on calcule une dispersion Δz, suivant l’axe vertical z, des mesures 51 ; autrement dit, on calcule la différence entre la coordonnée z la plus élevée et la coordonnée z la plus faible des mesures 51. En référence à la , qui est un graphe représentant une distribution statistique des coordonnées z des mesures 51 dans un exemple, on comprend que plus Δz est élevé, plus haute est la plus haute aspérité de la paroi porteuse de fond 15. Δz est donc un paramètre qui quantifie la non-planéité de la paroi porteuse de fond 15. Toujours à l’étape 108, on compare Δz avec une plage acceptable de non-planéité, qui a été préalablement spécifiée à l’étape 101. À l’étape 108, on calcule également un plan d’interpolation global 601 à partir de toutes les mesures 51 (cf. ). Ensuite, on recherche les deux mesures 51 les plus éloignées du plan 601, parallèlement à l’axe vertical z, l’une dans la direction positive de l’axe vertical z et l’autre dans la direction négative de l’axe z ; et on calcule la différence Δd entre les coordonnées z de ces deux mesures 51. En référence à la , qui est un graphe représentant une distribution statistique des coordonnées z des mesures 51 dans un autre exemple, on comprend que de cette manière, Δd quantifie la non-planéité de la paroi porteuse de fond 15 d’une façon différente de Δz. Toujours à l’étape 108, on compare Δd avec une plage acceptable de non-planéité, qui peut ou non être identique à celle à laquelle on compare Δz. En variante, une seule des dispersions Δz ou Δd peut être calculée. De façon avantageuse, la dispersion Δz peut n’être calculée que sur les mesures 51 qui sont situées suffisamment loin des pans porteurs verticaux 14. Plus concrètement, grâce à la connaissance a priori de la valeur de référence Q 0 du rayon Q, on peut exclure du calcul de la dispersion Δz les mesures 51 dont la coordonnée cylindrique r vérifie r > R 0 – R 2 où R 2 est un seuil prédéterminé. Il en est avantageusement de même pour la dispersion Δd. Ceci peut permettre d’éviter que l’évaluation de la non-planéité de la paroi porteuse de fond 15 tienne compte du raccordement entre la paroi porteuse de fond 15 et les pans porteurs verticaux 14, plus particulièrement lorsque ce raccordement présente la forme d’un chanfrein. Le plan d’interpolation global 601 peut aussi être utilisé pour contrôler des déformations locales de la paroi porteuse de fond 15, comme on va le décrire en référence à la . Pour un secteur 201 donné, une valeur σx ayant été spécifiée, on recherche les mesures 51 dont les coordonnées (a) sont incluses dans un intervalle de largeur σx (cf. ) ; parmi ces mesures 51, on recherche dans le secteur 201 les deux mesures 51 les plus éloignées du plan 601, parallèlement à l’axe vertical z, l’une dans la direction positive de l’axe vertical z et l’autre dans la direction négative de l’axe z ; on calcule une différence σz entre les coordonnées suivant l’axe vertical z de ces deux mesures 51 ; et on compare cette différence σz avec une plage acceptable de déformation locale de la paroi de fond, qui a été préalablement spécifiée à l’étape 101. Ce processus est répété jusqu’à épuisement du secteur 201, en décalant à chaque fois l’intervalle de largeur σx d’un incrément fixe ; autrement dit, l’intervalle de largeur σx est utilisé comme fenêtre glissante. Les opérations qui précèdent sont répétées pour chacun des secteurs 201. Les principes décrits ci-dessus pour une paroi porteuse de fond 15 plane sont également applicables à une paroi de couvercle plane de la structure porteuse 10. Les procédés de contrôle de géométrie décrits ci-dessus peuvent être appliqués à toute structure de génie civil ou toute structure mécano-soudée devant être construite en conformité avec le modèle spécifié. Les procédés de contrôle de géométrie décrits ci-dessus peuvent être mis en œuvre au moyen d’un dispositif de contrôle de géométrie 3000 (cf. ). Le dispositif de contrôle de géométrie 3000 comprend au moins un processeur 3001 et au moins une mémoire 3002 contenant un programme d’ordinateur ; la au moins une mémoire 3002 et le programme d’ordinateur sont configurés, avec ledit au moins un processeur 3001, pour causer l’exécution du procédé de contrôle de géométrie 100 décrit précédemment. Le dispositif de contrôle de géométrie 3000 peut être réalisé sous différentes formes, de manière unitaire ou distribuée, au moyen de composants matériels et/ou logiciels. Des composants matériels utilisables sont les circuits intégrés spécifiques ASIC, les réseaux logiques programmables FPGA ou les microprocesseurs. Des composants logiciels peuvent être écrits dans différents langages de programmation, par exemple C, C++, Java (marque déposée) ou VHDL. Cette liste n’est pas exhaustive. Bien que l'invention ait été décrite en liaison avec plusieurs modes de réalisation particuliers, il est bien évident qu'elle n'y est nullement limitée et qu'elle comprend tous les équivalents techniques des moyens décrits ainsi que leurs combinaisons si celles-ci entrent dans le cadre de l'invention. L’usage du verbe « comporter », « comprendre » ou « inclure » et de ses formes conjuguées n’exclut pas la présence d’autres éléments ou d’autres étapes que ceux énoncés dans une revendication. Dans les revendications, tout signe de référence entre parenthèses ne saurait être interprété comme une limitation de la revendication. Procédé de contrôle de géométrie (100) mis en œuvre par ordinateur, comprenant : - spécifier (101) un modèle de structure (10), le modèle de structure (10) comportant des attributs géométriques et des attributs dimensionnels, les attributs géométriques comportant un nombre N de côtés d’un polygone, N étant un nombre entier supérieur ou égal à 3, les attributs dimensionnels incluant une plage acceptable de diamètre (D) ; - fournir (102) une pluralité de premières mesures (50) représentant chacune une position mesurée dans un repère en trois dimensions (x 1 , y 1 , z 1 ) d’un point situé sur une structure à contrôler ; - déterminer (103) un système de coordonnées cylindriques (r, θ, z) présentant un axe vertical (19) parallèle à un axe dudit repère en trois dimensions, ledit axe vertical (19) présentant une position qui minimise une dispersion desdites premières mesures (50) selon une cordonnée radiale (r) dans ledit système de coordonnées cylindriques ; - spécifier (104) des premiers secteurs disjoints (200) de l’espace dans le système de coordonnées cylindriques (r, θ, z), lesdits premiers secteurs (200) étant définis par des incréments de coordonnées azimutales (δθ) et des incréments de coordonnées axiales (δz) dans ledit système de coordonnées cylindriques (r, θ, z) ; - calculer (105), à partir des premières mesures (50), une pluralité de premiers points de discrétisation (250), chaque premier point de discrétisation (250) représentant une position moyenne de la structure à contrôler dans l’un desdits premiers secteurs (200) ; - calculer (106), à partir des premiers points de discrétisation (250), des positions estimées (300) de N arêtes verticales (13) de la structure à contrôler ; - calculer (107), à partir des positions estimées (300) des N arêtes verticales (13), des positions estimées (400) de N médianes, chaque médiane étant une médiane d’un pan vertical de la structure à contrôler ; et - à partir des positions estimées des N médianes, comparer (108) des distances entre des médianes diamétralement opposées avec ladite plage acceptable de diamètre (D). Procédé de contrôle de géométrie (100) selon la revendication 1, dans lequel les attributs dimensionnels incluent en outre une plage acceptable de rayon (Q), et dans lequel le procédé de contrôle de géométrie (100) comprend en outre comparer (108) des distances entre lesdites positions estimées (400) des N médianes et ledit axe vertical (19) avec ladite plage acceptable de rayon. Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 2, dans lequel les attributs dimensionnels incluent en outre une plage acceptable de largeur de pan (W), et dans lequel le procédé de contrôle de géométrie (100) comprend en outre comparer (108) des distances entre les positions estimées (300) de deux arêtes verticales voisines (13) avec ladite plage acceptable de largeur de pan. Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 3, dans lequel les attributs dimensionnels incluent en outre une plage acceptable de non-verticalité d’arête (NV), et dans lequel le procédé de contrôle de géométrie (100) comprend en outre comparer (108) des distances, perpendiculairement à l’axe vertical (19) dudit système de coordonnées cylindriques, entre des points d’extrémité de chaque arête verticale (13) avec ladite plage acceptable de non-verticalité d’arête (NV). Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 4, dans lequel les attributs dimensionnels incluent en outre une plage acceptable d’ovalisation (OV), et dans lequel le procédé de contrôle de géométrie (100) comprend en outre comparer (108) une différence de diamètre entre un plus grand cercle inscrit (CI) de la paroi verticale (12) et un plus petit cercle circonscrit (CC) de la paroi verticale (12) de la structure à contrôler avec ladite plage acceptable d’ovalisation. Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 5, dans lequel les attributs dimensionnels incluent en outre une plage acceptable de déformation locale des pans verticaux (14), et dans lequel le procédé de contrôle de géométrie (100) comprend en outre comparer des déformations locales de la structure à contrôler avec ladite plage acceptable de déformation locale des pans verticaux (14). Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 6, dans lequel les attributs dimensionnels incluent en outre une plage acceptable de non-verticalité des pans verticaux (14), et dans lequel le procédé de contrôle de géométrie (100) comprend en outre : - sélectionner un sous-ensemble des premiers points de discrétisation (250), dans lequel sous-ensemble les coordonnées azimutales (Θ) des premiers points de discrétisation (250) sont égales entre elles ; - déterminer des écarts (ΔR) de coordonnée radiale (R) entre les premiers points de discrétisation (250) du sous-ensemble ; et - comparer les écarts (ΔR) de coordonnée radiale (R) avec ladite plage acceptable de non-verticalité des pans verticaux (14). Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 7, dans lequel les attributs géométriques comportent en outre une paroi de fond plane, les attributs dimensionnels incluent en outre une plage acceptable de non-planéité de la paroi de fond, et le procédé de contrôle de géométrie (100) comporte en outre : - fournir (102) une pluralité de deuxièmes mesures (51) représentant chacune une position mesurée dans ledit repère en trois dimensions (x 1 , y 1 , z 1 ) d’un point situé sur une paroi de fond (15) de la structure à contrôler ; - associer (103) aux deuxièmes mesures (51) un système de coordonnées cartésiennes (x, y, z), ledit système de coordonnées cartésiennes présentant un premier axe, un deuxième axe et un troisième axe perpendiculaires entre eux, le troisième axe coïncidant avec ledit axe vertical (19) dudit système de coordonnées cylindriques ; - spécifier (104) des deuxièmes secteurs disjoints (201) de l’espace dans ledit système de coordonnées cartésiennes, lesdits deuxièmes secteurs (201) étant définis par des incréments (δx, δy) de coordonnées selon le premier axe et le deuxième axe dudit système de coordonnées cartésiennes (x, y, z) ; - calculer (105), à partir des deuxièmes mesures (51), une pluralité de deuxièmes points de discrétisation (251), chaque deuxième point de discrétisation (251) représentant une position moyenne de la paroi de fond (15) de la structure à contrôler dans l’un desdits deuxièmes secteurs (201) ; - calculer (108) un paramètre de non-planéité (Δz ; Δd) à partir des deuxièmes points de discrétisation (251), et comparer (108) le paramètre de non-planéité (Δz ; Δd) avec ladite plage acceptable de non-planéité de la paroi de fond. Procédé de contrôle de géométrie (100) selon la revendication 8, dans lequel les attributs dimensionnels incluent en outre une plage acceptable de déformation locale de la paroi de fond, et dans lequel le procédé de contrôle de géométrie (100) comprend en outre comparer des déformations locales de la paroi de fond de la structure à contrôler avec ladite plage acceptable de déformation locale de la paroi de fond. Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 9, dans lequel N est un nombre pair, plus particulièrement N est un nombre pair compris entre 8 et 56, plus particulièrement encore N = 8 ou N = 56. Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 10, comprenant en outre une étape de génération consistant à générer (109) une représentation visuelle d’un résultat d’au moins une dite comparaison (108). Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 11, comportant en outre l’étape d’attribuer une qualification de conformité à la structure à contrôler, la qualification de conformité étant choisie parmi conforme et non-conforme. Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 12, dans lequel la pluralité de premières mesures (50) provient d’une acquisition de mesures (1001) au moyen d’un instrument de mesure dans la structure à contrôler, et dans lequel l’instrument de mesure comporte un appareil de télédétection par laser (91) disposé dans un espace interne de la structure à contrôler. Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 13, dans lequel la structure à contrôler est une structure porteuse pour une installation de stockage de gaz liquéfié (1). Procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 14, dans lequel la structure à contrôler est en béton. Programme d’ordinateur comprenant des instructions qui, lorsque le programme d’ordinateur est exécuté par un ordinateur, conduisent l’ordinateur à exécuter le procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 15. Dispositif de contrôle de géométrie (3000) comprenant au moins un processeur (3001) et au moins une mémoire (3002) contenant un programme d’ordinateur, dans lequel ladite au moins une mémoire (3002) et le programme d’ordinateur sont configurés, avec ledit au moins un processeur (3001), pour causer l’exécution du procédé de contrôle de géométrie (100) selon l’une quelconque des revendications 1 à 15 par le dispositif de contrôle de géométrie (3000).