Où poser un câble sous-marin n'est pas d'abord une question de profondeur, mais de stabilité du fond. Un tracé impeccable sur la carte bathymétrique peut traverser des sédiments qui se remettent en mouvement à chaque tempête, dénudent le câble et l'exposent à l'affouillement. Nous avons construit un outil qui chiffre ce risque cellule par cellule, puis en déduit le corridor de moindre coût entre deux atterrages.
Pourquoi la bathymétrie ne suffit pas
La profondeur dit où le câble peut passer physiquement, pas s'il y restera stable. Deux phénomènes gouvernent la tenue d'un câble posé ou ensouillé : la mobilité sédimentaire — la fraction du temps où le fond est mis en mouvement par la houle et les courants — et l'affouillement extrême, le creusement ponctuel lors des événements rares. Un sable fin peut être mobile 8 % de l'année là où un fond grossier voisin ne l'est que 0,6 % ; à quelques centaines de mètres près, le risque varie d'un ordre de grandeur. C'est cette information, absente d'une carte de profondeur, que l'atlas reconstruit.
La colonne vertébrale : des données hétérogènes sur une grille unique
Tout part de sources publiques reprojetées sur une seule grille de référence de 500 m (EPSG:3035). Chaque variable vient de sa source naturelle :
- Courants de fond — GLORYS12 (réanalyse Copernicus Marine), près du fond.
- États de mer (hauteur, période, vitesse orbitale au fond) — hindcast RESOURCECODE / HOMERE, horaire sur plus de 30 ans.
- Bathymétrie, nature de fond (classes de Folk), densité de trafic, zones d'exclusion — EMODnet.
Les extractions brutes sont mises en cache une fois puis réutilisées ; la physique lourde tourne aux nœuds de vagues du hindcast plutôt que sur les milliers de mailles, et seules les statistiques sont interpolées sur la grille. Le résultat est un cube unique où houle, courant, fond et contraintes s'alignent dans le temps et l'espace — la base que tout le reste vient lire.
Mesurer la mobilité : la contrainte de cisaillement au fond
La mobilité se déduit de la contrainte de cisaillement combinée houle + courant, comparée au seuil de mise en mouvement du sédiment. Nous appliquons la formulation de Soulsby (1997) pour la contrainte de fond, et le critère de Soulsby & Whitehouse pour la contrainte critique, calibré sur le d50 déduit de la nature de fond. La mobilité d'une maille est alors la fraction du temps où la contrainte dépasse ce seuil.
Sur notre tuile de démonstration (mer du Nord, janvier 2020), le fond est nettement dominé par la houle — le rapport contrainte de vagues sur contrainte de courant avoisine 8 — avec un maximum de cisaillement au 90ᵉ centile autour de 0,18 Pa et au 99ᵉ de 0,33 Pa. Chaque maille reçoit une classe (stable / intermédiaire / mobile) et un nombre d'événements par an, exportés en carte géoréférencée.
Chiffrer l'affouillement extrême par analyse de valeurs extrêmes
Le risque de dimensionnement ne vient pas de la moyenne mais des événements rares, que nous traitons par analyse de valeurs extrêmes (POT-GPD). À partir de la série de contraintes, la méthode « pics au-dessus d'un seuil » ajuste une loi de Pareto généralisée et produit les niveaux de retour à 1, 10, 50 et 100 ans, assortis d'un intervalle de confiance à 95 %.
Concrètement, l'atlas répond à la question de l'ingénieur d'ouvrage : quelle sollicitation le fond subit-il une fois par siècle à cet endroit, et avec quelle incertitude ? Une analyse de sensibilité (robustesse à ±1 classe de Folk sur le d50) accompagne chaque résultat pour ne pas masquer les hypothèses.
Du risque au tracé : le corridor de moindre coût
Les couches de risque sont fusionnées en une surface de coût, sur laquelle un algorithme A* trace le corridor optimal entre deux atterrages. Le coût de chaque maille combine, avec des poids ajustables, la mobilité, l'affouillement extrême, le risque d'ancrage (dérivé du trafic) et les exclusions réglementaires (Natura 2000, zones interdites). Là où un atterrage tombe sur une cellule barrée, il est ramené automatiquement sur la maille franchissable la plus proche.
L'outil ne propose pas un tracé unique mais trois variantes de Pareto — mini-risque, équilibré, mini-longueur — pour que la décision reste celle du client. Rejouer un jeu de poids recalcule le corridor en moins d'une seconde, ce qui en fait un support de discussion et non une boîte noire.
Une méthode agnostique de la zone
Changer de zone d'étude, c'est changer un fichier de configuration, pas le code. L'emprise, la résolution et les sources sont des paramètres ; rien d'une zone n'est codé en dur. La même chaîne, validée sur la mer du Nord, tourne aujourd'hui sur une tuile de la Manche — avec ses courants de marée horaires et ses contraintes propres — sans réécriture. C'est ce qui rend l'étude reproductible et transposable d'un projet à l'autre.
La physique et le solveur de tracé sont des modules isolés, testés unitairement et relus par un humain — pas des composants délégués à l'aveugle. C'est la condition pour signer un résultat d'ingénierie.
Ce que reçoit le client
Le livrable central est un atlas autoportant (un rapport HTML avec provenance, méthodologie, hypothèses et résultats) accompagné d'un service web : une API et un visualiseur cartographique en lecture seule, où l'on interroge un point (mobilité, niveau de retour à 100 ans, coût) ou l'on redessine le corridor en déplaçant des curseurs de poids. Les couches, stockées au format Zarr, se déploient aussi bien sur une machine que sur un stockage objet compatible S3.
Vous préparez la pose ou l'ensouillage d'un câble, d'un ombilical ou d'un pipeline ? Décrivez-nous votre zone et nous cadrons l'atlas de risque le plus sobre qui réponde à votre décision. Pour le choix des sources en amont, voir aussi notre guide des données océaniques.
Questions fréquentes
Quelles données faut-il pour évaluer la stabilité d’un câble sous-marin ?
Des courants de fond (GLORYS12), un hindcast de vagues longue durée (RESOURCECODE ou HOMERE), et la bathymétrie, la nature de fond, le trafic et les zones d’exclusion (EMODnet), le tout reprojeté sur une grille de référence commune.
Comment estime-t-on la mobilité des sédiments ?
On calcule la contrainte de cisaillement au fond combinée houle + courant (formulation de Soulsby), comparée au seuil de mise en mouvement du sédiment (Soulsby & Whitehouse). La mobilité est la fraction du temps où cette contrainte dépasse le seuil critique.
Peut-on appliquer la méthode à n’importe quelle zone ?
Oui. La chaîne est agnostique de la zone : l’emprise, la résolution et les sources sont des paramètres d’un fichier de configuration. Une nouvelle zone = un nouveau fichier de config et un run, sans réécriture du code.