Tsunarisque
| ,Chapitre 1. Reconstitution et modélisation du tsunami
L’analyse spatiale des dommages sur le bâti : contribution méthodologique et enseignements pour les futurs scénarios de risque tsunami
Résumé
The tsunami of the 26th of December 2004 caused an unprecedented disaster in the Indian Ocean. In Sumatra, a third of Banda Aceh area has been destroyed and 70,000 people died. The Tsunarisque Program brings new considerations about the tsunami dynamic and damages intensity in this urban area: an original method of damages spatial analysis is based on field surveys and photo-interpretations. The first result is a very accurate cartography of the tsunami breaking that is shown by a steep drop of the damaging gradient around 2,7 km from the coast. The second is a new “macro-tsunamic” intensity scale based on special typologies of buildings and damages. This analysis is complemented by damages curves that give the statistical relations between mean damages intensities and wave heights. These results allow developing application in tsunami potential losses modelling.
Texte intégral
1. Introduction
1Avec près de 70 000 morts (JICA, 2005) et des milliers de bâtiments détruits, l’agglomération de Banda Aceh (Sumatra, Indonésie) représente la zone de plus forte densité de dégâts liés à l’impact du tsunami du 26 décembre 2004 sur le pourtour de l’océan Indien. Au lendemain de la catastrophe, les premières évaluations satellitaires des dommages laissaient apparaître une agglomération dévastée à des degrés divers sur près des deux tiers de son territoire (SERTIT, 2004). Cette première analyse largement illustrée par les médias du monde entier permettait de prendre la mesure de l’exceptionnelle gravité de cet événement et de mobiliser l’aide d’urgence internationale. Elle fut suivie par les premières reconnaissances de terrain ponctuelles sur les dommages, notamment en Thaïlande (Ghosh et al., 2005 ; Srivichai et al., 2005 ; Thanawood C. et al., 2006) et en Indonésie (Borrero, 2005 ; Ghobarah et al., 2006). Dans un second temps (mars 2005), un premier diagnostic terrain complet des dommages aux infrastructures (réseaux vitaux et bâtiments) était réalisé par des équipes scientifiques japonaises. Ce travail a permis de dégager quelques bilans statistiques matériels et humains assortis de cartographies au 1/50 000epour préparer le schéma de reconstruction de la ville (JICA, 2005).
2Par la suite, aucune restitution cartographique détaillée des dommages à l’échelle locale (inférieure au 1/10 000e) ne semble avoir été menée à des fins scientifiques de reconstitution du phénomène et de ses effets. De telles approches spatiales basées sur un diagnostic terrain des pathologies du bâti sont pourtant courantes après de gros séismes destructeurs. Assurées habituellement par des spécialistes de la construction parasismique et des géologues, elles ont permis des avancées significatives dans la connaissance des mécanismes du phénomène et de ses effets, au moyen de cartes d’intensité macrosismique. Or, aucune démarche équivalente n’est décrite à ce jour dans la littérature scientifique pour les tsunamis.
3Notre recherche vise précisément à combler cette lacune. Elle privilégie une approche géographique et quantitative des dommages au bâti sur un secteur représentatif de Banda Aceh. Cela nécessite au préalable de définir un protocole méthodologique original mais reproductible, et de comprendre les modalités d’endommagement propres aux tsunamis. L’objectif est de mieux caractériser certains aspects du phénomène comme l’intensité de ses effets et la position de la zone de déferlement. On en tire également des bilans statistiques sur 892 bâtiments, à partir desquels sont ébauchées des courbes d’endommagement. Elles viendront alimenter les futurs scénarios de risque tsunami et de quantification des pertes associées, à l’instar de ce qui se pratique déjà pour les séismes. Les applications futures pourront porter sur des contextes géographiques proches, en termes d’occupation du littoral et de nature du bâti, en particulier dans les DOM-TOM.
2. Secteur d’étude et méthodologie
4Notre analyse s’est portée sur un secteur de 7 km2 compris entre le littoral de Uleelheue et la grande place de Kampung Baru (place de l’Avion) (Fig .1). Il s’agit d’un secteur représentatif de la diversité du bâti et des modalités d’endommagement. Il offre une forte densité de bâtiments avec en moyenne 885 unités au km2. D’ouest en est, ce secteur comprenait une portion littorale formée d’un lido (Uleelheue) et de lagunes aménagées en parcs à crevettes, des zones résidentielles entrecoupées de petites rizières et une zone d’habitat très concentré proche du centre-ville. Le tout est encore relié par une route principale en position diagonale nord-ouest-sud-est par rapport à notre cadre d’étude. L’est de la zone est traversé par un canal d’orientation subméridienne. Le choix du terrain a été défini avant tout sur la base d’un gradient décroissant nord-est-sud-ouest des hauteurs d’eau mesurées par rapport au sol, allant de 14,1 à 0 m. Il s’agit d’un secteur fortement éprouvé par le tsunami. Outre ses commodités d’accès, ce secteur présente l’intérêt d’être très bien couvert depuis l’espace : avant la catastrophe, par une scène du satellite Quick Bird du 23 juin 2004 (résolution 60 cm) ; après la catastrophe, par deux missions photographiques aériennes de janvier 2005 et juin 2005 (30 cm), une scène Quick Bird du 28 décembre 2004 (60 cm) et une scène Ikonos du 3 janvier 2005 (1 m).
5Un diagnostic systématique du bâti a été conduit sur le terrain du 16 au 25 janvier 2006, soit un peu plus d’un an après l’événement, en pleine phase de reconstruction. Néanmoins, les dommages étaient encore très visibles malgré quelques opérations de réparation de certains bâtiments. Par commodité de repérage, nous avons découpé notre secteur d’étude en mailles de 500 m de côté. Du fait de l’absence de cartes précises de la zone et des transformations du terrain (routes et chemins détruits) opérées par le tsunami, nous avons privilégié un repérage précis sur des fonds photographiques aériens de juin 2005, préalablement redressés et géoréférencés (orthophotos). La qualité de ces fonds (résolution 30 cm) a permis de nous passer d’un GPS et de pouvoir constituer la base d’un système d’information géographique très précis. Globalement, peu de transformations ont été constatées dans l’occupation du sol entre la période de notre mission et la date des photo aériennes utilisées (5 mois), si ce n’est l’apparition de nouvelles constructions.
6Après des premiers repérages de terrain dans la zone sélectionnée mais aussi en d’autres secteurs de la ville, nous avons établi une typologie simplifiée du bâti et des modes et niveaux d’endommagement du bâti visibles sur le terrain. Nous avons tenté d’établir en permanence une correspondance visuelle entre les observations de terrain et les clichés aériens. Sur le terrain, nous avons accordé une priorité d’analyse aux bâtiments partiellement détruits (1 230 unités). Ceux totalement rasés par le tsunami ou au contraire non affectés (angle SE) ont été identifiés par photo-interprétation (4 970 unités). Les bâtiments totalement disparus après déblaiement, ou ceux dont les dalles ont été recouvertes de gravats ou réutilisées pour reconstruire, ont été identifiés par photo-interprétation d’extraits (gratuits) des scènes Quick Bird d’après la catastrophe (source : http://www.digitalglobe.com).
7Le repérage de terrain a été complété par des photographies au sol permettant d’affiner le diagnostic en laboratoire. Ce diagnostic repose sur une analyse rapide et visuelle des dommages et ne nécessite pas de connaissances poussées en ingénierie des structures. L’attention s’est portée sur les atteintes aux éléments de la structure du bâtiment, leurs liaisons et les différents types de matériaux utilisés. Des codes correspondant aux classes de vulnérabilité des bâtiments et à leur niveau d’endommagement ont été portés sur chaque bâtiment identifié. L’ensemble de la base de données (6 200 bâtiments renseignés) a fait l’objet d’un SIG structuré dans le logiciel MapInfo® (Fig. 2). Les traitements cartographiques suivants ont été également assurés par un module d’analyse spatiale (Vertical Mapper®).
3. Description des modalités d’endommagement
8Bien que d’une magnitude (Mw) exceptionnelle de 9.3 (Stein et al. 2005), les effets du séisme ont été minimes à Banda Aceh, comparativement à ceux du tsunami. L’examen des dommages sur les bâtiments préservés par le tsunami le confirme, tout comme les témoignages des survivants collectés au cours de notre mission ou par d’autres équipes (Boen, 2006, p .S807), ou encore une séquence vidéo amateur tournée au moment des secousses.
9Le tsunami est avant tout un phénomène hydraulique de submersion marine caractérisé à terre par un écoulement à forte charge solide non newtonien (Chanson, 2005). L’action dommageable d’un tsunami se traduit sur le continent par plusieurs processus pouvant se succéder et/ou se surimposer dans le temps et l’espace. Chaque processus produit des mécanismes de ruine spécifiques qui renseignent à leur tour sur la magnitude du phénomène. Les reconnaissances de terrain ont permis de diagnostiquer 7 processus majeurs associés au passage d’un tsunami sur le milieu terrestre. Chacun peut être décrit par un ou plusieurs paramètres physiques (critères de magnitude) en théorie quantifiable(s) et modélisable(s). En lien avec la dynamique d’une vague de type plongeant (cas de la deuxième vague), on peut différencier ces processus élémentaires en nature et en nombre selon que l’on se trouve en amont (zone de transit ou de transfert), au sein, ou à l’aval (zone d’inondation) de la zone de déferlement (cassage de la vague). Les zones de transit et de déferlement se caractérisent par des écoulements turbulents de très forte énergie tandis que la zone d’inondation donne des écoulements laminaires beaucoup moins rapides. Ce découpage spatial simplifié et donc chronologique devient évidemment beaucoup plus complexe à modéliser ou à reconstituer lorsque plusieurs vagues se succèdent, comme ce fut le cas à Banda Aceh.
10Les processus dominants associés sont des actions hydrodynamiques temporaires impliquant des forces de pression et d’impact qui sont fonction de la vitesse de l’écoulement, de sa hauteur et de sa charge solide. Y sont également associés des effets d’érosion (abrasion, affouillement) et d’accumulation de débris divers, arrachés, charriés puis déposés par le tsunami. Les gros objets flottants comme les bateaux ou les voitures peuvent provoquer des effets d’impact et d’embâcle très dommageables pour les infrastructures, parfois à plusieurs kilomètres du rivage. On retrouve à Banda Aceh un bateau de 63 m de long et 20 m de large appartenant à la compagnie nationale d’électricité (PLN) et qui a été déposé à 2,6 km de sa zone de stationnement habituelle. Le retrait des eaux (backwash) joue un rôle dommageable mineur sur les infrastructures. La fin du processus laisse subsister des poches d’inondation plus ou moins permanentes d’eau salée préjudiciables pour certains matériaux. Un an après le passage du tsunami, certains secteurs affouillés, très anciennement marécageux (puis asséchés), ou affaissés, laissaient entrevoir des bâtiments avec les fondations dans l’eau.
11Toute tentative de mise au point de matrices ou de fonctions d’endommagement nécessite de pouvoir quantifier au préalable les processus d’endommagement avec des critères physiques simples débouchant sur des échelles de magnitude. La vitesse des écoulements et leur hauteur sont les deux paramètres physiques les plus accessibles et pertinents pour tenter de modéliser a priori ou a posteriori le potentiel de destruction d’un tsunami, faute de pouvoir reconstituer facilement le champ des pressions et d’en déduire les effets directs sur les infrastructures. Des auteurs mentionnent des pressions d’environ 6,8 t/m2 sur le front de mer de Kamala Beach (Phuket) en Thaïlande pour une hauteur de vague de 3 m (Srivichai et al., 2005). Les forces de pression exercées sont en effet fonction des masses en mouvement et donc des innombrables débris transportés, ce qui laisse supposer un champ des déformations sur le bâti très anisotrope et aléatoire. Les nombreux embâcles formés par ces mêmes débris et les effets d’abriliés aux bâtiments ne feraient que complexifier une éventuelle modélisation. À Banda Aceh, la vitesse de l’écoulement a pu être déterminée ponctuellement à partir de deux vidéos amateurs (source http://www.asiantsunamivideos.com). L’une tournée à environ 2,5 km de la côte dans le centre de Banda Aceh (près de la grande mosquée) montre un flot de débris de 1,5 m de hauteur dont le front avance à une vitesse de l’ordre de 1,5 à 1,7 m/sec (Chanson, 2005). Il s’agit pour Lavigne et al. (2008, cet ouvrage) de la première vague en fin de course. L’autre vidéo tournée à 8 h 50 en zone périurbaine au sud-ouest de l’agglomération à 2,4 km de la côte, soit près de 36 minutes après la deuxième vague, montre un écoulement beaucoup plus fluide mais très chargé, de 3 m de hauteur et dont la vitesse est environ 2 fois supérieure à la précédente.
12À côté de ces quelques valeurs ponctuelles de pression et de vitesse, l’autre grandeur exploitable spatialement pour caractériser la magnitude du tsunami est la hauteur atteinte par l’onde de submersion sur les terres. Elle est en théorie proportionnelle au carré de la vitesse. La mesure et la cartographie systématique des hauteurs de submersion ont fait l’objet de plusieurs campagnes, américaines, japonaises et françaises. La cartographie obtenue par l’équipe française dans le cadre du programme de recherche Tsunarisque offre un gradient de hauteur d’eau s’étendant de 0 à 15 m sur notre secteur d’étude. Rappelons qu’au contact des côtes, la vitesse de l’onde marine décroît tandis que son amplitude augmente pour satisfaire l’équation de conservation de l’énergie, jusqu’au moment où la vague déferle. Jusqu’à ce point de rupture, la hauteur des vagues est donc proportionnelle à la capacité de destruction du tsunami (son énergie potentielle) et constitue en cela un très bon indicateur de magnitude du phénomène. Les échelles de magnitude proposées successivement depuis 75 ans par Imamura & Iida, Soloviev, Abe, Hatori, Murty & Loomis sont basées sur des valeurs maximales de hauteurs de vague (Hmax en m) sur la côte (Papadopoulos & Imamura, 2001). Les valeurs maximales mesurées sur la côte dans notre secteur d’étude de Uleelheue (14 m) correspondent à une magnitude de 3 sur l’échelle d’Imamura & Iida (M = log2Hmax) et de 9,1 selon la classification de Abe (1981).
13La diversité des modes d’action du tsunami a entraîné des dommages variables sur le bâti, en nature et en intensité, qui vont de la légère dégradation à l’arasement total avec déblaiement. La première vague, probablement poussée par la deuxième, a semble-t-il surtout endommagé le bâti semi-permanent et léger (sederhana), comme le suggère la quantité de matériaux transportés et visibles sur la vidéo prise au niveau de la grande mosquée. En supposant que les effets dommageables majeurs ont été le fruit de la deuxième vague, on peut mettre en relation chaque processus avec les modes d’endommagement les plus fréquemment rencontrés sur le terrain. Les modes d’endommagement dominants sont liés aux effets d’impact et de poussées latérales exercés par le front du tsunami tout au long de la zone de transfert. Les bâtiments semblent avoir explosé au contact des vagues comme le confirment certains témoignages oraux établissant une analogie sonore entre l’arrivée du tsunami et un bombardement. Dans ces secteurs de hauteurs d’eau mesurées entre 6 et 14 m (jusqu’à environ 2,5 km de la côte selon une direction nord-ouest-sud-est) ne subsistent que les dalles des habitations. Les bâtiments présentaient pourtant le plus souvent une structure armée et un à deux étages. Les piliers ont été arrachés ou couchés dans le sens du courant et les débris disséminés probablement sur plusieurs kilomètres. On retrouve le plus souvent seulement les réservoirs à eau des salles de bain (ou mandi). Les secteurs intermédiaires de submersion comprise entre 3 et 6 m (entre 2,5 et 3,4 km de la côte environ) présentent des dommages beaucoup plus variables en intensité. Ces derniers sont surtout fonction de la résistance du bâti à des effets d’impacts ciblés. Les murs extérieurs sont souvent perforés et les constructions partiellement effondrées suite à des ruptures de poteaux. Au-delà d’une zone de déferlement supposée et toujours en direction du sud-est, les dommages sont à peine visibles un an après. Ils sont cependant présumés par les amoncellements de débris visibles sur les clichés aériens ou terrestres pris juste après le tsunami. Il s’agit de secteurs d’inondation comprise entre 0 et 3 m et concernés surtout par une submersion plus lente et donc par moins d’effets d’impact. Sur les bâtiments conservés, les dommages liés au séisme se distinguent facilement de ceux induits par le tsunami. Sur le terrain, ces différenciations ont été facilitées par les descriptions rapportées par une équipe d’ingénieurs du génie civil canadiens (Ghobarah et al., 2006). Rappelons que les séismes provoquent des oscillations du sol entraînant des dommages aux structures par efforts de balancement, de cisaillement et de flexion des structures, avec des pathologies caractéristiques (Zacek, 1996). Par ailleurs, dans les secteurs de submersion comprise entre 6 et 14 m, il semble logique que les bâtiments fragilisés par le séisme aient moins bien résisté au tsunami, mais sans que l’on puisse en apporter la démonstration.
4. Mise au point de typologies du bâti et des dommages
14Conformément à la méthodologie explicitée en début d’article, une analyse quantitative et spatiale des dommages au bâti a été entreprise. Pour cela nous avons établi, après plusieurs phases de validation sur le terrain, deux typologies simplifiées et facilement réutilisables. Elles concernent d’une part la nature du bâti et d’autre part les principaux modes et niveaux d’endommagement de ce bâti.
15Nous introduisons 5 classes de bâtiments de vulnérabilité décroissante (A, B, C, D, E) assez facilement distinguables sur le terrain et sur photos aériennes selon des critères de forme, de dimension et de nature de toiture (Pl. 1). Cette classification (Tab. 1) se base entre autres sur les observations de nature et de vulnérabilité effectuées par Ghobarah et al. (2006) et Boen (2006). Le type A est le moins résistant. Il s’agit de constructions très légères faites en bois et en tôles, parfois sur pilotis, comme c’était le cas en bordure des lagunes du secteur de Uleelheue. Le type B correspond à du bâti maçonné (briques) à structure légère comportant des poteaux de 20 × 20 cm renforcés par 4 fers lisses de 8 mm de diamètre aux quatre coins de l’édifice. C’est le type le plus représenté dans le secteur d’étude. Il s’agit le plus souvent d’habitations autoconstruites appartenant à des classes sociales moyennes. Le type C correspond toujours à une habitation individuelle du style villa, parfois de haut standing. La construction présente une géométrie en plan plus complexe et toujours un étage. La structure est renforcée et calculée. Le type D est celui des bâtiments collectifs (blocs) à usage d’habitation, de commerce ou de service public (enseignement, administration). La structure est en général moyennement renforcée pour la taille de ces bâtiments à plusieurs poteaux, de un à deux étages, mais de géométrie simple, qui dominent dans le centre-ville de Banda. Le type E caractérise également des bâtiments collectifs, le plus souvent à usage administratif, religieux (mosquées) ou privé (hôtels), mais de très bonne conception avec des structures en béton armé renforcées et calculées. Les nombreuses mosquées présentes sur le secteur d’étude (16) ont été regroupées dans la classe E précédente, compte tenu de la qualité de leur construction et de leur bonne transparence hydraulique. Sur les photos aériennes, des incertitudes de distinction peuvent se poser entre les classes A et B, ou B et C, lorsque les toitures sont respectivement semblables à A ou C. Une classe NI (non identifié) a été réservée aux très nombreux bâtiments rasés (4095), dont il est difficile de connaître la nature initiale a posteriori. La typologie retenue est assez adaptée à la classification du bâti dans le contexte de pays en développement des zones tropicales. Cela concerne aussi la plupart des territoires français d’outre-mer exposés à d’éventuels tsunamis (Martinique, Guadeloupe, Guyane, La Réunion, Mayotte, Tahiti, etc.).
16La typologie des dommages proposée est spécifique aux tsunamis et diffère de celles existant pour les séismes (Grünthal et al., 2001) ou les mouvements de terrain (Leone et al., 1996) par la nature des dommages décrits. Les dommages y sont hiérarchisés en 6 classes suivant une gravité croissante (D0 à D5) (Tab .2 et Pl. 4). Un seuil de dégât déterminant pour la stabilité de l’édifice et son éventuelle réoccupation est franchi en passant du niveau D2 à D3. Cela correspond à l’atteinte aux structures rigidifiant le bâtiment (poteaux verticaux et transversaux, planchers). Mis à part les destructions partielles ou totales (D4, D5) nettement identifiables par photo-interprétation, les autres classes nécessitent des observations in situ. La principale difficulté d’interprétation peut survenir avec le niveau D3 qui exige un examen plus détaillé des superstructures, sans pour autant requérir des compétences particulières en ingénierie. Compte tenu de la densité d’observations renseignées et géoréférencées, tant sur les dommages que sur le phénomène (hauteurs, directions, chronologie), il nous est apparu intéressant de parfaire notre analyse spatiale des dommages par un essai d’interpolation sous SIG.

PL .1. Illustration des différentes classes de vulnérabilité des bâtiments (photos F. Leone, janvier 2006).
5. Analyse spatiale des dommages
17L’interpolation des dommages suppose que ces derniers suivent une évolution continue dans l’espace. Elle repose sur les bâtiments de classes NI (entièrement détruits, D5) et B. Ce choix est déterminé par la recherche d’un comportement moyen homogène des structures pour les niveaux d’endommagement inférieurs à D5. De la sorte, les variations spatiales observables dans l’endommagement entre D0 et D4 seront logiquement et principalement liées à l’action et à la magnitude du phénomène plutôt qu’à la résistance du bâti (classe B à comportement homogène).

PL .2.Illustration des niveaux de dommages (D2 à D5) pour différentes classes de bâtiments (photos : F. Leone, janvier 2006).
18L’interpolation numérique a été réalisée par krigeage sous Vertical Mapper® (Gratton, 2002), à partir du jeu de données ponctuelles (5 693 unités) pour aboutir à un fichier raster renseigné de type Grid. Une vectorisation dans MapIn+fo® a ensuite permis d’extraire les contours des 5 niveaux d’endommagement.
19La carte obtenue montre très nettement un gradient d’endommagement diminuant d’ouest en est et qui s’estompe dans l’angle sud-est de la zone d’étude, dans des secteurs très densément urbanisés (Fig. 3). Toute la partie orientale de la zone, sur près de 2 km depuis l’ancien rivage de Uleelheue, a subi un endommagement maximum, avec une prédominance de l’arasement total des bâtiments (D5). Cela correspond à des hauteurs d’eau comprises entre 6 et 14 m. Des îlots d’endommagement plus faible (D2 à D4) apparaissent en se rapprochant de l’ouest. Ils se répartissent sur une bande clairsemée, d’orientation générale nord-est-sud-ouest, de 800 m de large. Ils révèlent des effets de site à mettre en relation avec une baisse d’énergie cinétique localisée de la vague, des effets d’abri (ou d’écran) provoqués par certains bâtiments ou bien plus probablement d’un meilleur comportement de certains bâtiments mieux construits. Des effets d’écran très ponctuels mais non cartographiés sont identifiables dans certains secteurs. C’est notamment ce que l’on a pu observer à proximité de la Grande Mosquée de Uleelheue qui a partiellement protégé une habitation de type B dans un secteur où la submersion a pourtant atteint près de 14 m de haut en front de mer. La zone clairsemée est relayée à l’est par une discontinuité très nette d’orientation générale NE-SO située à une distance moyenne d’environ 2,7 km depuis le littoral de Uleelheue. Le gradient d’endommagement y chute brutalement de D5 à D1 sur près de 100 m dans un secteur où les hauteurs d’eau passent soudainement de 6 à 2 m. Nous interprétons cette bande continue et étroite comme étant la zone de déferlement de la deuxième vague, donc d’un changement de dynamique et d’énergie de la vague avec ce que cela implique en termes de processus et de potentiel d’endommagement sur le bâti. Le niveau D2, facilement identifiable sur le terrain, est caractéristique de cette zone de déferlement où la submersion, plus basse et en perte d’énergie, est venue seulement défoncer les murs extérieurs et parfois les cloisons internes. Cette hypothèse est confirmé par la position des bateaux n° 4 et n° 5 qui se retrouvent précisément échoués sur cette zone de déferlement probable (Pl. 3). La position des bateaux n° 2 et n° 3 peut s’expliquer par un déferlement plus précoce dans un secteur qui n’a pas été correctement mis en évidence par l’interpolation, mais seulement esquissé. Enfin, plusieurs témoins localisent aussi un déferlement dans le secteur du bateau PLN (n° 4). La zone de déferlement est représentée sur la carte en se calant approximativement sur les limites de D5 (amont) et D2 (aval). Sa largeur est variable suivant les secteurs, ce qui renforce un tracé lobé qui peut s’expliquer par des différences de hauteur de vague. Cela semble être le cas dans le secteur B7 où l’on observe une extension plus large du déferlement qui coïncide avec des hauteurs d’eau sans doute un peu plus hautes qu’ailleurs à cette distance de la côte (autour de 6 m). La présence d’une route transversale dans ce secteur a pu également étaler le déferlement et ses effets sur les constructions. Deux autres secteurs présentent aussi une dilatation de la zone de déferlement, mais qui a été probablement détachée du front par l’interpolation sous forme de lobes isolés, à cheval sur les secteurs C6 et C7. Il existe là aussi un lien avec des hauteurs d’eau plus hautes (9 m) mesurées 500 m en amont, et, plus au sud, dans le secteur D6, mais sans démonstration équivalente. Par ailleurs, on remarque une très bonne corrélation spatiale entre la position du front de déferlement et les indicateurs d’orientation de la deuxième vague, relevés sur le terrain en amont de ce front par Lavigne et al. (2006). On constate un découplage de l’orientation des écoulements dans le secteur B5 en deux directions générales distinctes (130° N et 150° N), qui pourrait s’expliquer par un partage de la vague. Cette hypothèse suggère la présence, plus à l’ouest, d’une possible zone de déferlement révélée par des poches aval de moindre dommage, bien visibles dans le secteur C5, et compatible avec la position des bateaux n° 2 et n° 3. Au-delà de cette ou de ces zones de déferlement, l’endommagement est beaucoup plus faible (D0 à D1) et correspond surtout à des effets de submersion de faible hauteur (moins de 2 m) et de stagnation des eaux.
6. Proposition et validation d’une échelle d’intensité « macrotsunamique »
20La carte précédente des dommages interpolés n’intègre pas les observations sur les autres catégories de bâtiments. La prise en compte de ces informations complémentaires et utiles nécessite de mettre au point une échelle d’intensité des dommages intégrant les niveaux d’endommagement (D0, D1, D2, D3, D4, D5) et les classes de vulnérabilité du bâti (A, B, C, D, E). Cette combinaison empirique a pu se faire à partir de nos observations de terrain. Elle débouche sur une échelle d’intensité des dommages à 6 degrés (Tab. 3). Bien que basée sur des critères qualitatifs, mais clairement explicités, cette échelle d’intensité « macrotsunamique » (par analogie avec les échelles macrosismiques) permet de quantifier et de spatialiser l’endommagement sur un secteur géographique donné. En tant qu’échelle d’intensité, elle respecte la condition d’indépendance vis-à-vis des critères de magnitude du phénomène pour se baser uniquement sur les observations des dégâts. Sa pertinence est révélée par le fait qu’elle permet de restituer intégralement le gradient d’endommagement obtenu pour les seuls bâtiments B et surtout de retrouver précisément la zone de déferlement précédemment cartographiée (passage brutal de VI à I) (Fig .4). Mieux encore, cette échelle permet de révéler et de compléter la position de la zone de déferlement supposée à l’extrémité nord-est de la zone d’étude, c’est-à-dire dans un secteur qui présentait peu de bâtiments de classe B .Elle permet aussi d’affiner le découpage des secteurs à effets de site ainsi que le détail de la ligne de déferlement en confirmant en particulier sa position plus occidentale (et éventuellement plus précoce) dans la zone d’échouage des bateaux n° 2 et n° 3. Cette échelle présente donc un réel intérêt pour positionner une ligne de déferlement à partir de l’analyse fine des dommages sur des bâtiments de vulnérabilités différentes, et ce au même titre qu’une échelle macrosismique (EMS98 par exemple) permet de retrouver l’épicentre d’un séisme et de mettre en évidence une source (faille active) et des effets de site géologiques ou topographiques. Notre échelle est adaptée à une approche détaillée, rigoureuse et locale de l’endommagement.
TAB .3. Échelle d’intensité « macrotsunamique » basée sur les dommagesaux bâtiments (degrés I à VI).


FIG .4. Carte des intensités de dommages interpolées (tous types de bâtiments) (échelle macro-tsunamique).
21Elle diffère en cela des échelles d’intensité des dommages proposées par Papadopoulos et al. (2001) et leurs prédécesseurs (Sieberg, 1927 ; Ambraseys, 1962 ; Soloviev, 1970) qui reposent sur des critères beaucoup plus qualitatifs, subjectifs et variés (dommages, perceptions, effets environnementaux, etc.) et donc peu adaptés à une analyse spatiale fine. Enfin, un essai de corrélation entre l’intensité des dommages (ID), tous types de bâtiments confondus, et la magnitude du phénomène exprimée en hauteur de submersion par rapport au sol est présenté dans le paragraphe suivant. Il donne le tableau de correspondance suivant (Tab. 4) dans lequel nous avons également porté pour indication l’échelle d’intensité donnée par Papadopoulos et al. (2001).
TAB .4. Correspondances entre les intensités macrotsunamiques (I) et les hauteurs de submersion (H).

7. Élaboration de courbes d’endommagement
22L’élaboration de scénarios d’endommagement ou de pertes liés à un tsunami nécessite d’utiliser des matrices ou fonctions permettant d’établir la correspondance entre la magnitude du phénomène et les dommages attendus. En général, ces matrices ou fonctions d’endommagement s’expriment par la probabilité, pour un bâtiment de classe de vulnérabilité donnée, d’atteindre un certain niveau d’endommagement moyen : Dmoy. On peut aussi les exprimer par la probabilité, pour un stock de bâtiments d’une classe donnée, d’atteindre un certain pourcentage d’un niveau d’endommagement D donné. Les premières fonctions s’appliquent donc à l’échelle d’un bâtiment, les secondes à un volume de bâtiments. Les unes et les autres servent à quantifier les pertes économiques, du moment où les valeurs de remplacement du bâtiment ou du stock sont connues (Leone et al., 1996 ; Leone, 2004). Or de telles fonctions appliquées aux tsunamis sont absentes dans la littérature scientifique, contrairement à ce que l’on observe dans le champ de la sismologie et des scénarios sismiques (FEMA 310, 1998 ; Giovinazzi et al., 2003 ; RISK-UE, 2003 ; AFPS, 2005), celui des vents forts (Murlidharan et al., 1997 ; Khanduri et al., 2003 ; Stewart, 2003) ou des inondations (MOC, 1996 ; Dutta et al., 2003). Ce constat est probablement lié à la rareté des retours d’expériences pluridisciplinaires et détaillées sur les dommages provoqués par un tsunami.
23Nous avons tenté d’esquisser et d’illustrer ces deux familles de fonctions d’endommagement structurel (ou de vulnérabilité) en construisant les courbes s’y rapportant. Pour cela, nous avons échantillonné nos bâtiments, soit 892 unités, dans un rayon de 100 m autour de chaque point de mesure de hauteur d’eau par rapport au sol. Cela permet d’obtenir un bon taux de couverture de la zone tout en limitant au mieux les recouvrements (2,8 %).
24La première approche est appliquée aux bâtiments de classe B, les plus fréquents dans la zone. Le nuage de points obtenu donne une relation non linéaire de type logarithmique présentant une assez bonne qualité d’ajustement (R2 = 0,58). Nous nous sommes calés sur cette portion de courbe pour établir la courbe donnant le niveau d’endommagement moyen (Dmoy) en fonction de la hauteur d’eau mesurée au sol (H). Un ajustement manuel a été effectué aux limites n en se basant sur nos observations de terrain et de manière à borner la courbe par les niveaux D0 = 0 et D5 = 5. Entre ces points, l’évolution de Dmoy est considérée comme continue en fonction de la hauteur de submersion (Fig. 5). La même démarche appliquée à l’intensité moyenne des dommages (IDmoy), donnée par notre échelle, aboutit à une meilleure corrélation (R2 = 0,66). Elle donne une courbe d’allure semblable permettant théoriquement de retrouver l’intensité des dommages à partir de la magnitude du phénomène (H), et ce quelle que soit la classe des bâtiments analysés.
25La seconde approche suit le même principe de construction (R2 = 0,59). La courbe obtenue exprime le pourcentage de destruction totale (niveau D5) d’un stock de bâtiments (en %) en fonction de la hauteur de submersion (H). Elle est limitée par les valeurs 0 et 100 %. Cette courbe n’a qu’une valeur indicative car le bâti pris en compte pour D5 n’a pas été différencié. Par exemple, on obtiendrait probablement une évolution différente avec une plus grande proportion de bâtiments plus résistants et plus hauts (2 étages), notamment à partir des valeurs de H inférieures à 6 m. Compte tenu du fait que les pertes sont donc modulées par les classes de vulnérabilité, il conviendra de pousser l’analyse pour chaque classe de bâtiment après identification de la nature initiale des bâtiments NI (non identifiés) par photo-interprétation des images Quick Bird de juin 2004.
26Dans les deux approches, les autres corrélations ne sont pas significatives, hormis pour le niveau D1 qui présente une évolution globalement et logiquement inverse par rapport à D5. En vue d’une application opérationnelle dans le cadre de scénarios de perte, ces courbes et fonctions qui en découlent doivent être définies pour chaque classe de bâtiment. Cela suppose de multiplier de tels retours d’expériences et d’envisager des modélisations des interactions tsunami/structures. Par ailleurs, une très bonne corrélation linéaire (R2 = 0,91) apparaît entre les pertes humaines, exprimées en % de morts et disparus, et le pourcentage de destruction totale du bâti (niveau D5). Celle-ci est établie sur 16 quartiers (Desa) totalement ou en grande partie inclus dans notre zone d’étude et pour lesquels le JICA (2005) nous fournit des données sur les victimes. La droite de régression obtenue sur ces quelques points décrit une fonction de pertes humaines, aussi sommaire soit-elle.
27Enfin, le passage entre des niveaux de dommages structurels et des coûts monétaires directs et/ou indirects nécessite aussi de construire des fonctions de pertes économiques et de délais de réoccupation pour chaque niveau de dommage. Nos observations de terrain permettent seulement de dire que les bâtiments sont majoritairement abandonnés et reconstruits à partir du niveau D4 et que les niveaux D2 et D3 impliquent des travaux de remise en état pouvant aller respectivement de quelques jours à une année. Ceci dépend d’une multitude de facteurs d’ordre organisationnel, institutionnel, politico-administratif, socio-économique, culturel voire psychosociologique. Les travaux menés en complément de notre démarche par Vinet et al. (2006) apportent un éclairage très pertinent pour comprendre ce processus de réoccupation et de reconstruction après un tsunami dans un contexte de pays en développement.
8. Conclusion
28Ce travail est entièrement fondé sur l’analyse systématique des dommages au bâti à partir d’observations in situ et selon un protocole méthodologique original et reproductible. Il a conduit à la mise au point de typologies adaptées à une reconnaissance rapide de terrain et par photo-interprétation. La structuration et la densité spatiale des informations récoltées ont permis d’établir une cartographie fine du gradient d’endommagement et de révéler une zone de déferlement dont l’existence n’avait été évoquée que par quelques témoins. L’emprise spatiale de cette zone de déferlement a pu être affinée grâce à une échelle d’intensité macrotsunamique permettant d’intégrer l’ensemble du bâti endommagé. Les correspondances entre les 6 degrés de cette échelle et les hauteurs moyennes de submersion ont également été déduites.
29Cette analyse démontre tout l’intérêt d’une lecture géographique des dommages en complément des traditionnels relevés de hauteurs de submersion posttsunami pour reconstituer la dynamique et l’extension du phénomène. Elle souligne le rôle déterminant de cette dynamique sur les mécanismes de destruction et l’importance de la zone de déferlement en tant que ligne de discontinuité majeure dans le gradient d’endommagement. Les travaux menés jusqu’à présent sur la cartographie des dommages posttsunami par télédétection n’ont pas permis d’arriver à ce degré de précision ni de renseigner la dynamique du phénomène (JICA, 2005, vol .III ; Gruhier, 2005 ; SERTIT, 2004 ; Greidanus et al., 2005).

PL .3. Les principaux bateaux échoués dans ou à proximité de la zone d’étude (photos : F. Leone & F. Vinet, janvier 2006). voir Fig .4 pour situation des n° 3, 4 et 5.
30La pertinence des futurs scénarios ou modèles de pertes liées à un tsunami va donc largement dépendre de la capacité à simuler des hauteurs et vitesses d’eau sur les terres et donc la zone de déferlement de la vague principale. Ces scénarios nécessiteront dans un second temps d’intégrer des fonctions d’endommagement tsunami adaptées aux enjeux présents (bâtiments, routes, personnes, activités, etc.). Nous en avons ébauché certaines dans ce travail avec l’espoir de les améliorer et de les compléter à l’occasion de futurs retours d’expériences.
31Le seul modèle de mesure du risque tsunami décrit par la littérature n’est pour l’instant qu’au stade de développement en Grèce. Il est basé sur une approche multicritère de vulnérabilité et non sur des fonctions d’endommagement (Papathoma, 2003 et Dominey-Howes & Papathoma, 2006). Par ailleurs et en complément utile, Sugimoto et al. (2003) ont proposé un modèle de simulation des délais d’évacuation et de mise à l’abri de populations littorales au Japon en cas de tsunami. Ces scénarios, très utiles pour la prévention, méritent d’être développés dans des contextes socio-géographiques similaires à celui de Banda Aceh. Cela suppose d’intégrer des démarches, des outils et des données pluridisciplinaires, et de travailler à la lumière de retours d’expériences aussi riches que celui entrepris en Indonésie dans le cadre du programme Tsunarisque.
Table des illustrations
![]() | |
---|---|
Légende | FIG. 1. Situation de la zone d’étude au sein de la ville de Banda Aceh. |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-1.jpg |
Fichier | image/jpeg, 226k |
![]() | |
Légende | FIG .2. Carte des bâtiments diagnostiqués. |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-2.jpg |
Fichier | image/jpeg, 285k |
![]() | |
Titre | TAB .1. Classification du bâti en fonction d’une vulnérabilité décroissante. |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-3.jpg |
Fichier | image/jpeg, 1006k |
![]() | |
Titre | TAB .2. Classification des modes et niveaux d’endommagement du bâti suivant une gravité croissante. |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-4.jpg |
Fichier | image/jpeg, 731k |
![]() | |
Légende | PL .1. Illustration des différentes classes de vulnérabilité des bâtiments (photos F. Leone, janvier 2006). |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-5.jpg |
Fichier | image/jpeg, 285k |
![]() | |
Légende | PL .2.Illustration des niveaux de dommages (D2 à D5) pour différentes classes de bâtiments (photos : F. Leone, janvier 2006). |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-6.jpg |
Fichier | image/jpeg, 270k |
![]() | |
Légende | FIG .3. Carte des niveaux d’endommagement interpolés (bâtiments de classes B et NI). |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-7.jpg |
Fichier | image/jpeg, 147k |
![]() | |
Titre | TAB .3. Échelle d’intensité « macrotsunamique » basée sur les dommagesaux bâtiments (degrés I à VI). |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-8.jpg |
Fichier | image/jpeg, 62k |
![]() | |
Légende | FIG .4. Carte des intensités de dommages interpolées (tous types de bâtiments) (échelle macro-tsunamique). |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-9.jpg |
Fichier | image/jpeg, 154k |
![]() | |
Titre | TAB .4. Correspondances entre les intensités macrotsunamiques (I) et les hauteurs de submersion (H). |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-10.jpg |
Fichier | image/jpeg, 42k |
![]() | |
Légende | FIG. 5. Courbe d’endommagement moyen (bâtiments de classe B). |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-11.jpg |
Fichier | image/jpeg, 222k |
![]() | |
Légende | PL .3. Les principaux bateaux échoués dans ou à proximité de la zone d’étude (photos : F. Leone & F. Vinet, janvier 2006). voir Fig .4 pour situation des n° 3, 4 et 5. |
URL | http://books.openedition.org/psorbonne/docannexe/image/3800/img-12.jpg |
Fichier | image/jpeg, 198k |
© Éditions de la Sorbonne, 2011