PRÉLIMINAIRES

I. Épistémologie et Enjeux Scientifiques du Domaine

La sismologie, initialement une science d’observation des catastrophes, a muté en une discipline d’imagerie quantitative de l’intérieur terrestre. Cette transformation, propulsée par la puissance de calcul et la théorie de la propagation des ondes, place le traitement du signal au cœur de la géophysique moderne. L’enjeu n’est plus seulement de localiser un séisme, mais de déchiffrer la signature mécanique des roches pour en extraire des informations sur les ressources (minérales, hydriques, énergétiques) et les risques (failles actives, stabilité du sol). Cette UE formalise cette transition épistémologique radicale.

II. Cartographie des Compétences et Transversalité

Cette unité d’enseignement forge une compétence unifiée à la croisée de trois piliers. Le traitement des images satellitaires et sismiques constitue le socle technique, permettant de passer de la donnée brute à une information intelligible. Cette information est ensuite mobilisée pour évaluer quantitativement les ressources et les risques, une compétence analytique cruciale pour l’ingénieur. Enfin, la modélisation et la digitalisation de cette information géographique dans un SIG achèvent le cycle, transformant l’analyse scientifique en un outil d’aide à la décision stratégique, directement exploitable par les planificateurs et les industriels.

III. Alignement Stratégique avec les Réalités Opérationnelles

La maîtrise de l’interprétation sismique quantitative répond à une demande critique du marché africain, notamment en RDC, où le sous-sol constitue un levier économique majeur. L’Ingénieur Géophysicien et Modélisateur devient l’architecte de la connaissance du sous-sol, indispensable aux secteurs minier et pétrolier. L’Expert en Télédétection spatiale et le Spécialiste SIG, armés de ces compétences, sont capables de contextualiser ces découvertes, de surveiller les impacts environnementaux et de planifier l’aménagement du territoire en intégrant les contraintes géologiques, des compétences vitales pour le développement durable des infrastructures.

Chapitre I. Fondations Physiques et Numériques du Signal Sismique

I.1 Propagation des Ondes et Modèle Convolutif

Ancrée dans la mécanique des milieux continus, la propagation des ondes sismiques (P, S, surface) constitue le langage physique que nous devons déchiffrer. Ce sous-chapitre formalise l’équation d’onde et établit le modèle convolutif comme pierre angulaire de la sismique réflexion. La trace sismique est ainsi conceptualisée comme la convolution de la réflectivité terrestre avec la signature de la source, filtrée par les effets de propagation. La maîtrise de ce modèle est le prérequis absolu pour toute action de traitement ou d’interprétation ultérieure.

I.2 Discrétisation du Signal et Transformée de Fourier

Sous l’angle du traitement numérique, le signal sismique continu est une abstraction. Sa capture par les géophones impose un échantillonnage, régi par le théorème de Shannon, dont la violation génère un repliement de spectre (aliasing) irréversible. Ce segment explore les implications pratiques de la discrétisation et introduit la Transformée de Fourier comme l’outil mathématique essentiel pour passer du domaine temporel au domaine fréquentiel. L’étudiant apprendra à analyser le contenu spectral d’une trace pour diagnostiquer la qualité du signal et préparer les étapes de filtrage.

I.3 Limites du Modèle Idéal : Atténuation, Bruit et Anisotropie

Critiquement, le modèle convolutif simple se heurte à la complexité du sous-sol réel. L’atténuation géométrique et intrinsèque dissipe l’énergie des ondes, modifiant leur amplitude et leur phase de manière non-linéaire, un effet particulièrement prononcé dans les terrains sédimentaires ou volcaniques. Ce module analyse les mécanismes physiques de l’atténuation et de l’anisotropie, ainsi que la nature du bruit sismique (cohérent et aléatoire). Comprendre ces limitations est fondamental pour ne pas interpréter un artefact de propagation comme une structure géologique réelle.

I.4 Application : Audit d’un Jeu de Données Sismiques Brut (Format SEG-Y)

Face à un fichier SEG-Y, la première action de l’ingénieur consiste à en évaluer la qualité. Cet exercice pratique guide l’étudiant dans le chargement d’un jeu de données sismiques 2D réel, issu du bassin côtier congolais, à l’aide de bibliothèques Python (telles qu’ObsPy ou Segyio). L’objectif est de réaliser un diagnostic technique : vérifier les paramètres d’acquisition dans l’en-tête, calculer le spectre d’amplitude moyen, visualiser les sections de tirs brutes et identifier les principaux types de bruits présents.

Chapitre II. Acquisition en Sismique Réflexion : De la Source au Récepteur

II.1 Géométrie d’Acquisition et Couverture Souterraine (CDP)

Héritage des explorations pétrolières, le concept de point-milieu-commun (CMP) ou point de réflexion commun (CDP) est une révolution géométrique. Il permet de sonder la même portion du sous-sol sous différents angles, une redondance essentielle pour augmenter drastiquement le rapport signal/bruit. Ce sous-chapitre dissèque les différentes géométries d’acquisition (2D, 3D, “split-spread”, “end-on”) et démontre mathématiquement comment le choix du dispositif en surface conditionne la qualité et la résolution de l’image finale en profondeur, notamment la couverture et le pli.

II.2 Mécanismes de la Source Sismique et Couplage au Sol

La nature de la source sismique contrôle la signature initiale de l’ondelette injectée dans le sol. Des sources impulsives (dynamite, marteau) aux sources vibratoires (Vibroseis), chaque technologie possède une signature spectrale, une directivité et une répétabilité qui lui sont propres. Ce segment analyse la physique de ces sources et le concept critique de couplage au sol, particulièrement complexe en environnements de forêt dense ou de sols latéritiques. Un mauvais couplage est la première source de dégradation irrécupérable du signal sismique.

I.3 Analyse Critique des Artefacts d’Acquisition

La perfection en acquisition est une illusion. L’empreinte d’acquisition (“acquisition footprint”), due à une illumination non-uniforme du sous-sol, et les ondes non-réfléchies (onde directe, air-wave, ground-roll) constituent des bruits cohérents puissants qui masquent les réflexions profondes. Ce module enseigne à reconnaître la signature de ces artefacts sur les enregistrements bruts. L’analyse de leur vitesse apparente et de leur contenu fréquentiel dans le domaine F-K (fréquence-nombre d’onde) est la clé pour concevoir des filtres efficaces lors du traitement.

II.4 Application : Conception d’une Campagne Frugale pour l’Hydrogéophysique

Pour répondre aux besoins en eau potable de la périphérie de Kinshasa, la cartographie des aquifères superficiels est une priorité. Ce cas d’étude impose la conception d’une campagne de sismique réflexion 2D à haute résolution et à faible coût. L’étudiant devra choisir une source adaptée (marteau instrumenté), définir l’espacement optimal des géophones pour éviter l’aliasing du ground-roll, et calculer le pli de couverture nécessaire pour imager des aquifères situés entre 30 et 100 mètres de profondeur, le tout avec un budget matériel contraint.

Chapitre III. Traitement du Signal Sismique : De la Donnée Brute à la Section Lisible

III.1 Correction des Effets Statiques et Géométriques

La topographie de surface et les variations d’épaisseur de la zone altérée (weathering layer) introduisent des décalages temporels qui distordent la cohérence des réflexions. Ce sous-chapitre aborde les corrections statiques, une étape cruciale et délicate, ainsi que la correction de “Normal Moveout” (NMO). Cette dernière compense la géométrie de l’acquisition pour transformer les hyperboles de réflexion en événements plats sur les gathers CDP, préparant le terrain pour l’étape de sommation (stack). Une mauvaise correction NMO dégrade la résolution verticale.

III.2 Outils de Rehaussement du Signal : Filtrage et Déconvolution

Au cœur du traitement sismique, la déconvolution vise à compresser l’ondelette sismique pour augmenter la résolution temporelle et à supprimer les réverbérations (multiples). Ce segment explore les différentes approches de déconvolution (prédictive, “spiking”) et leur fondement statistique. Parallèlement, les filtres fréquentiels (passe-bande) et spatiaux (F-K) sont déployés pour éliminer les bruits résiduels identifiés lors de l’analyse. L’étudiant apprendra à appliquer ces outils en séquence, chaque étape conditionnant le succès de la suivante dans une chaîne de traitement logique.

III.3 Limites de la Chaîne de Traitement Conventionnelle

La chaîne de traitement standard repose sur des hypothèses fortes, comme la propagation 1D pour la déconvolution ou une vitesse de NMO constante, qui sont souvent violées dans des contextes géologiques complexes. Ce module critique analyse les situations où cette approche échoue : présence de fortes hétérogénéités latérales, structures à fort pendage, ou présence de gaz superficiel. Reconnaître ces scénarios est vital pour éviter de produire une image sismique propre mais géologiquement fausse, et pour justifier le passage à des algorithmes plus avancés.

III.4 Application : Traitement d’une Ligne Sismique sur le Graben du Tanganyika

L’étudiant reçoit une ligne sismique 2D brute acquise sur le lac Tanganyika, caractérisée par une topographie de fond lacustre complexe, des structures de failles normales et une sédimentation rapide. La mission est de construire une chaîne de traitement complète dans un environnement logiciel open-source (ex: Seismic Unix). Il devra appliquer les corrections statiques flottantes, déterminer les vitesses de NMO par analyse de vitesse, appliquer une déconvolution et un filtre F-K pour atténuer les multiples de l’eau, et produire une section sommée interprétable.

Chapitre IV. Imagerie Sismique : Migration et Construction du Modèle de Vitesse

IV.1 Le Principe de Migration : Replacer l’Énergie Sismique

La section sommée (stack) n’est pas une image fidèle du sous-sol ; les événements ne sont pas à leur position réelle, notamment en présence de pendages. La migration est l’opération d’imagerie qui corrige cette distorsion en déplaçant l’énergie sismique des hyperboles de diffraction vers leurs apex. Ce sous-chapitre expose le fondement conceptuel de la migration, la présentant comme une solution inverse à l’équation d’onde. Elle transforme une section temporelle en une coupe géologique structurellement plus juste, augmentant la focalisation et la résolution latérale.

IV.2 Mécanismes de Migration : de Kirchhoff à la RTM

Les algorithmes de migration varient en complexité, en coût de calcul et en précision. Ce segment décortique les mécanismes de la migration de Kirchhoff, une approche intuitive basée sur la sommation le long des courbes de diffraction, et la compare aux méthodes basées sur l’équation d’onde comme la RTM (Reverse Time Migration). La RTM, bien que très coûteuse, est la seule capable d’imager correctement des structures extrêmement complexes comme les flancs de dômes de sel ou les réseaux de failles sous-volcaniques, fréquents dans le Rift Est-Africain.

IV.3 Le Rôle Central et la Critique du Modèle de Vitesse

Toute migration en profondeur est entièrement dépendante de la précision du modèle de vitesse de propagation des ondes. Une erreur de quelques pourcents sur le modèle de vitesse peut entraîner des erreurs de positionnement de plusieurs centaines de mètres en profondeur et une défocalisation complète de l’image. Ce module analyse de manière critique les méthodes de construction de ce modèle (“velocity picking”, tomographie) et les incertitudes associées. Le modèle de vitesse n’est pas une donnée, mais le résultat d’une interprétation itérative.

IV.4 Application : Migration en Temps d’une Structure Faillée au Kivu

À partir d’une section sommée d’une zone du Graben du Kivu, l’étudiant doit réaliser une migration en temps. La première étape consiste à l’analyse minutieuse des vitesses pour construire un modèle de vitesse lissé. Ensuite, en utilisant un logiciel de traitement, il appliquera un algorithme de migration (ex: Stolt ou Gazdag) pour transformer la section “stack” en une section migrée. La comparaison “avant/après” permettra de quantifier le gain en résolution et la correction du positionnement des failles normales, essentielles pour l’évaluation du risque sismique local.

Chapitre V. Interprétation Quantitative (IQ) : Extraction des Attributs Physiques

V.1 Au-delà de la Structure : le Concept d’Attribut Sismique

L’interprétation ne s’arrête pas au pointé des horizons et des failles. L’amplitude, la fréquence et la phase du signal sismique contiennent des informations quantitatives sur la lithologie, la porosité et le fluide contenu dans les roches. Ce sous-chapitre introduit le concept d’attribut sismique, qui consiste à calculer de nouvelles traces à partir des données sismiques pour mettre en évidence des caractéristiques géologiques subtiles. Des attributs comme l’enveloppe de la trace ou la fréquence instantanée transforment une image structurelle en une carte de propriétés physiques.

V.2 Mécanisme de l’AVO (Amplitude Versus Offset)

L’amplitude d’une réflexion n’est pas constante mais varie avec la distance source-récepteur (l’offset), un phénomène gouverné par les équations de Zoeppritz. L’analyse AVO (Amplitude Versus Offset) exploite cette variation pour discriminer les types de fluides (gaz, pétrole, eau) et les lithologies. Ce segment décortique la physique de l’AVO, les approximations (Shuey) utilisées en pratique et la classification des anomalies AVO. C’est l’outil principal de l’interprétation quantitative pour la détection directe d’hydrocarbures ou de réservoirs géothermiques.

V.3 Inversion Sismique : Limites et Non-Unicité

L’inversion sismique vise à transformer la section sismique (réflectivité) en un volume d’impédance acoustique, une propriété plus directement liée à la géologie. Cependant, ce problème inverse est mal posé et intrinsèquement non-unique : une infinité de modèles géologiques peut produire la même réponse sismique. Ce module aborde de front cette limitation fondamentale, en analysant comment l’intégration de données externes (logs de puits, modèles géologiques a priori) est nécessaire pour contraindre la solution et obtenir un résultat géologiquement plausible et non un simple artefact mathématique.

V.4 Application : Caractérisation d’un Réservoir Potentiel dans la Cuvette Centrale

Un prospect a été identifié sur une section sismique 3D de la Cuvette Centrale. La mission de l’étudiant est de mener une analyse quantitative pour évaluer son potentiel. Il devra d’abord réaliser une inversion sismique post-stack pour obtenir un cube d’impédance. Ensuite, il effectuera une analyse AVO sur les gathers pré-stack pour rechercher des anomalies indicatrices de la présence de gaz. La combinaison de ces deux approches permettra de proposer une cartographie du risque et de la probabilité de succès du prospect.

Chapitre VI. Synthèse Géospatiale : Intégration Sismique, Télédétection et SIG

VI.1 Le Principe de Fusion de Données Multi-Échelles

La connaissance géologique est renforcée lorsque des données de natures et d’échelles différentes sont confrontées. La sismique offre une vision en profondeur, tandis que la télédétection (images satellites optiques et radar) et la géologie de surface fournissent le contexte à grande échelle. Ce sous-chapitre établit les fondements théoriques de la fusion de données, montrant comment l’intégration de contraintes externes permet de valider ou d’invalider des hypothèses d’interprétation sismique et de réduire l’incertitude globale du modèle géologique.

VI.2 Outils d’Intégration : Co-visualisation et Corrélation Spatiale

Les Systèmes d’Information Géographique (SIG) sont la plateforme technique de cette synthèse. Ce segment se concentre sur les outils permettant de co-visualiser en 3D des horizons sismiques, des images satellites drapées sur un modèle numérique de terrain (MNT) et des données de puits. Au-delà de la visualisation, l’étudiant apprendra à utiliser des outils d’analyse spatiale pour corréler statistiquement l’expression en surface d’une faille (visible sur une image Sentinel) avec sa géométrie en profondeur (interprétée sur la sismique).

VI.3 Critique des Modèles Intégrés : Propagation des Incertitudes

L’intégration de multiples sources de données ne fait pas magiquement disparaître les incertitudes ; elle peut même les propager et les amplifier si le processus n’est pas rigoureusement contrôlé. Ce module adopte une posture critique en analysant les défis liés à l’harmonisation de données de résolutions, de précisions et de datums différents. La question centrale est de savoir comment construire un modèle géologique final qui quantifie et représente visuellement l’incertitude à chaque point, plutôt que de présenter une image faussement déterministe.

VI.4 Application : Modélisation du Risque de Glissement de Terrain à Bukavu

La ville de Bukavu est exposée à un risque élevé de glissements de terrain. Ce projet de synthèse finale charge l’étudiant de construire un modèle de risque intégré dans QGIS. Il devra combiner : 1) une carte des pentes dérivée d’un MNT radar, 2) une carte de lithologie interprétée à partir d’images multispectrales, 3) la proximité des failles actives identifiées par l’analyse de données sismiques locales et régionales. Le résultat est une carte de zonage du risque, un outil concret pour l’urbanisme et la protection civile.

ANNEXES

A. Maîtrise de GMT (Generic Mapping Tools)

Pour l’Ingénieur Géophysicien, la capacité à produire des cartes et des figures de qualité publication est non-négociable. GMT est une suite d’outils en ligne de commande, open-source et extrêmement puissante, qui permet d’automatiser la création de cartes géophysiques complexes (cartes bathymétriques, gravimétriques, sections sismiques géoréférencées). Son approche par script la rend idéale pour le travail reproductible et efficace, même sur des machines aux ressources limitées, un atout majeur pour un travail de terrain ou en laboratoire en RDC.

B. Déploiement de la Méthodologie InSAR

L’Expert en Télédétection spatiale doit maîtriser l’interférométrie radar (InSAR) pour quantifier les déformations de la surface terrestre avec une précision centimétrique. Cette annexe fournit un guide pratique pour télécharger les données gratuites du satellite Sentinel-1, les traiter avec des logiciels open-source (comme SNAP ou GMTSAR) pour générer des interférogrammes, et interpréter ces derniers pour mesurer la subsidence minière, le gonflement d’un volcan (comme le Nyiragongo) ou le mouvement d’une faille, liant directement la télédétection à la sismotectonique.

C. Manipulation du Format de Données SEED

Le format SEED (Standard for the Exchange of Earthquake Data) est la langue véhiculaire de la sismologie mondiale. Pour le Spécialiste SIG ou le Modélisateur, savoir extraire, lire et convertir ces données est fondamental pour intégrer les informations des réseaux sismologiques globaux dans des analyses locales. Cette annexe détaille la structure d’un fichier miniSEED, montre comment utiliser des outils (comme ObsPy) pour lire les formes d’onde, accéder aux métadonnées (stations, réponses instrumentales) et préparer les données pour des analyses de risque sismique ou des études de structure.

Lire aussi :

Cours de Télédétection Marine, master-2, TDM2231

Cours de Projet, master-1, PAP2122

Cours de Data mining, master-1, DAM2121

Cours de Algorithmique et programmation, licence-2, ALP1241

Cours de Stage professionnel, licence-3, CON1362

Cours de Techniques et structures, master-1, TST2121