REPUBLIQUE DU CAMEROUN REPUBLIC OF CAMEROUN

Paix – Travail – Patrie Peace – Work – Fatherland

******** *******
UNIVERSITY OF YAOUNDE I
UNIVERSITE DE YAOUNDE I HIGHER TEACHER TRAINING COLLEGE
ECOLE NORMALE SUPERIEURE DEPARTMENT OF physics
DEPARTEMENT DE Physiques *******
*********

Elaborationd'un code pour le filtrage des donnees

geophysiques par la methode polynomiale par moindres
carres de Thurston et Brown

Mémoire presente en vue de l'obtention du D.I.P.E.S II

Par :

DIMO HOUPA Ginette Sonia

Licenciee en Physique

Sous la direction
NDOUGSA MBARGA Theophile
Maitre de conferences

Année Académique
2015-2016
 AVERTISSEMENT

Ce document est le fruit d'un long travail approuvé par le jury de soutenance et mis à
disposition de l'ensemble de la communauté universitaire de Yaoundé I. Il est
soumis à la propriété intellectuelle de l'auteur. Ceci implique une obligation de
citation et de référencement lors de l’utilisation de ce document.

D'autre part, toute contrefaçon, plagiat, reproduction illicite encourt une poursuite
pénale.
Contact : [email protected]

WARNING

This document is the fruit of an intense hard work defended and accepted before a
jury and made available to the entire University of Yaounde I community. All
intellectual property rights are reserved to the author. This implies proper citation
and referencing when using this document.

On the other hand, any unlawful act, plagiarism, unauthorized duplication will lead
to Penal pursuits.
Contact: [email protected]
 ii

DEDICACE

Je dédie ce travail :

A ma chère maman SATIO Anne épouse HOUPA

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 iii

REMERCIEMENTS

Ce travail ne saurait être l'œuvre d'une seule personne, mais le fruit de la contribution
de plusieurs à qui il convient d'exprimer ma gratitude. Je rends tout d'abord grâce au Seigneur
Dieu Tout Puissant qui me guide chaque jour.

Je tiens à remercier:

- Les membres du Jury pour l'honneur qu’ils m’accordent en acceptant d'évaluer ce travail ;

- Le Professeur NDOUGSA MBARGA Théophile pour la confiance portée à mon égard et

qui, en dépit de ses nombreuses occupations a bien voulu diriger ce travail. Veuillez trouver
ici l'expression de ma profonde gratitude;

-le chef de département de Physique BEGUIDE Bonoma ;

- Les enseignants de l'ENS pour les multiples échanges constructifs, leurs conseils, leur
disponibilité tout au long de notre formation. Qu'ils trouvent en ces mots ma sincère
reconnaissance ;

- Les ainés du laboratoire NGOH Jean Daniel et YANDJIMAIN Justine pour leurs
disponibilités, leurs encouragements, leurs conseils et leurs remarques constructives qui ont
été pour nous une source de motivation et ont enrichi notre travail ;

-Tous mes camarades de promotion et amis en particuliers NGOUADJEU Inès, KOLOKO

Sorel, HEUTEU Crépin, SIEWE Titus ,YOUTCHOM Jean Lapluy, KAPNANG
Franky , MANEKENG Stella ,YOUMBI Dorota et BIBOUM Edouard pour leur franche
collaboration ;

- Ma famille et toutes autres personnes qui n’ont pas été cités pour tout le soutient qu’ils
m’ont apporté;

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 iv

TABLE DES MATIERES

DEDICACE ......................................................................................................................... ii

REMERCIEMENTS ........................................................................................................ iii

LISTE DES FIGURES ................................................................................................... vi

LISTE DES TABLEAUX...............................................................................................vii

LISTE DES ABBREVIATIONS .................................................................................viii

ABSTRACT ........................................................................................................................ x

INTRODUCTION GENERALE .................................................................................... 1

CHAPITRE I : REVUE DE LA LITTERATURE .................................................... 4

Introduction ............................................................................................................................ 4
1.1. Généralités sur les codes en géophysique ................................................................ 4
1.1.1. Définition et évolution du code ........................................................................ 4
1.1.2. Rôle et importance des codes en géophysique ................................................. 5
1.2. Présentation générale de Matlab .............................................................................. 6
1.2.1. Démarrage de Matlab ....................................................................................... 6
1.2.2. Interface de Matlab ........................................................................................... 7
1.2.3. Fichiers exécutables : les scripts ..................................................................... 10
1.2.4. Les fonctions sous-routines ............................................................................ 10
1.2.5. Les variables, les boucles et les instructions .................................................. 13
1.2.6. Compilation dans Matlab................................................................................ 14
1.3. Travaux géophysiques utilisant la méthode Thurston et Brown ............................ 14
Conclusion ............................................................................................................................ 15
CHAPITRE II : METHODOLOGIE ET PROGRAMMATION .......................... 16
Introduction .......................................................................................................................... 16
2.1. Théorie de la méthode gravimétrique ...................................................................... 16
2.1.1. Définition et fondement de base .................................................................... 16

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 v

2.1.1. Corrections et anomalies de la gravité ............................................................ 18

2.2. Méthodes de séparation régionale et résiduelle ..................................................... 20
2.2.1. Méthode graphique ......................................................................................... 20
2.2.2. Méthodes analytiques ..................................................................................... 20
2.2.3. Méthode de séparation analytique polynomiale ............................................. 21
2.3. Etape de la programmation ........................................................................................ 26
2.3.1. Variables utilisées ........................................................................................... 26
2.3.1.1. Types de variables ....................................................................................... 26
2.3.1.2. Données ....................................................................................................... 26
2.3.2. Etape du programme ....................................................................................... 26
2.4.3. Programme .......................................................................................................... 28
Conclusion ............................................................................................................................ 28
CHAPITRE III : RESULTATS ET DISCUSSION ............................................... 29
Introduction .......................................................................................................................... 29
3.1. Données gravimétriques ............................................................................................. 29
3.2. Etablissement, analyse et traitement des cartes d’anomalies ..................................... 31
3.2.1. Cartes de l’anomalie de Bouguer ........................................................................ 31
3.2.1.1. Les zones d’anomalies lourdes ..................................................................... 32
3.2.1.2. Zones d’anomalies légères............................................................................ 33
3.2.1.3. Zones de gradients ........................................................................................ 33
3.2.2. Séparation des anomalies de Bouguer ................................................................. 34
3.2.2.1. Calcul des anomalies régionales et résiduelles ............................................. 34
3.2.2.2. Les cartes d’anomalies séparées et leur analyse ........................................... 38
3.2.2.3. Choix de la régionale .................................................................................... 46
3.2.3. Synthèse sur l’analyse des cartes régionales et résiduelles d’ordre 3 ................. 46
3.3. Discussion et validité des résultats ............................................................................ 47
3.4. Implication dans le système pédagogique ................................................................. 49
Conclusion ............................................................................................................................ 49
CONCLUSION GENERALE ....................................................................................... 50

BIBLIOGRAPHIES ........................................................................................................ 52

ANNEXES ......................................................................................................................... 55

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 vi

LISTE DES FIGURES

Figure 1.1 : Interface graphique de Matlab .............................................................................. 7

Figure 1.2 : Editeur Matlab version 7.10.0 (R2010a) .............................................................. 9

Figure 2.1 : Deux masses dans un référentiel cartésien (Giroux B., 2008)….. .................... .17

Figure 3.1 Carte de l’anomalie de Bouguer (Zanga A., 2015) ............................................... 30

Figure 3.2 : Carte de l’anomalie de Bouguer obtenue à partir du logiciel oasis montaj ........ 32

Figure 3.3 : Carte régionale et résiduelle d’ordre 1 de la zone d’étude : (a) carte régionale ;
(b) carte résiduelle .................................................................................................................... 39

Figure 3.4 : Carte régionale et résiduelle d’ordre 2 de la zone d’étude : (a) carte régionale ;
(b) carte résiduelle .................................................................................................................... 41

Figure 3.5 : Carte régionale et résiduelle d’ordre 3 de la zone d’étude : (a) carte régionale ;
(b) carte résiduelle .................................................................................................................... 43

Figure 3.6 : Carte régionale et résiduelle d’ordre 4 de la zone d’étude : (a) carte régionale ;
(b) carte résiduelle .................................................................................................................... 45

Figure 3.7 : Carte des anomalies régionales et résiduelle faites par Zanga A., 2015 : (i) carte
régionale ; (ii) carte résiduelle .................................................................................................. 48

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 vii

LISTE DES TABLEAUX

Table 1.1 : les différentes fonctions vectorielles ............................................ ……………….11

Table 1.2 : les différentes fonctions matricielles .................................................................... 12

Table 3.1 : valeurs des coefficients calculés à l’ordre 1 à partir du code Matlab .................... 34

Table 3.2 : valeurs des coefficients calculés à l’ordre 2 ......................................................... 35

Table 3.3 : valeurs des coefficients calculés à l’ordre 3 .......................................................... 36

Table 3.4 : valeurs des coefficients calculés à l’ordre 4 .......................................................... 37

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 viii

LISTE DES ABBREVIATIONS

-E : Est ;

-ENE : Est Nord-Est ;

-ESE : Est Sud-Est ;

-Fig : Figure ;

-mGals : milliGals ;

-Km : Kilomètre ;

-N : Nord ;

-NE : Nord-Est ;

-NNE : Nord Nord-Est ;

-NW : Nord-Ouest ;

-R.C.A : République Centrafricaine ;

-S : Sud ;

-SE : Sud-Est ;

-SSW : Sud Sud-Ouest ;

-SW : Sud-Ouest ;

-W : Ouest ;

-WNW : Ouest Nord-Ouest ;

-WSW : Ouest Sud-Ouest ;

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 ix

RESUME

La méthode Thurston et Brown est une méthode de séparation polynomiale par

moindres carrés ; pour pouvoir effectuer cette séparation, un programme informatique a été
mis sur pied. De ce fait, Le but de ce travail est d’élaborer un code en Matlab effectuant une
séparation des anomalies régionales et résiduelles à partir de la Bouguer. L’exécution du code
s’est faite à partir des données gravimétriques issues de la numérisation de la carte d’anomalie
de Bouguer obtenue par Zanga (2015) ce qui nous a conduit à une séparation des anomalies.
Les données ainsi séparées pour différents ordre du polynôme ont permis de générer
différentes cartes à partir du logiciel Oasis montaj 6.4.
Une analyse nous a conduite à identifier les zones d’anomalies lourdes et légères ainsi que les
zones de gradients. Les résultats ont été comparé à ceux obtenus par les études faites par
Zanga (2015) à partir de la méthode des éléments finis et il a été constaté une similitude entre
la régionale et la résiduelle obtenu par la méthode des éléments finis et celle obtenue par la
séparation polynomiale par moindres carrés à l’ordre 1 tandis qu’à l’ordre 3 nous avons
constaté une divergence. Ceci pourrait être dû au fait que plus le degré du polynôme est élevé
plus la ressemblance s’estompe.

Mots clés : Code Matlab ; Anomalie de Bouguer ; filtrage polynomiale ; moindres carrés ;
anomalie résiduelle.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 x

ABSTRACT

The Thurston and Brown method is a method base on polynomial separation by least
squares; to be able to carry out this separation, a data-processing program has been set up. In
fact, the aim of this study is to work out a code in Matlab which can be able to make a
separation of the regional and residual anomalies from Bouguer anomalies. The execution of
this code was done by using gravimetric datum resulting from the digitalization of the map of
Bouguer anomalies made by Zanga (2015) what led us to a separation of the anomalies.
Datum separated for different order of the polynomial made possible to generate various map
by using Oasis montaj 6.4 Software.
An analysis led us to identify zones of heavy anomalies zones and gradients zones. Our
results have been compared with those obtained by the studies made by Zanga A., 2015 with
the finite element method and it was noted a similarity between regional and the residual with
the first order while with the third order we have noted a divergence. This could be due to the
fact that while the degree of the polynomial increases, the resemblance becomes blurred.

Keywords Matlab Code; Bouguer anomaly; filtering polynomial; least squares; residual
anomaly.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 1

INTRODUCTION GENERALE

a. Contexte de l’étude

L’analyse et le traitement des données pour le géophysicien est une étape importante pour
l’obtention des résultats. C’est dans ce contexte que le traitement des données autrefois
manuel était l’objet de beaucoup d’erreurs et nécessitant des longues périodes de temps. Avec
l’émergence de la technologie, les logiciels de traitement de données ont été développés.
Cependant l’utilisation et l’obtention de ces derniers ne sont pas toujours aisées vu leur coût
onéreux ; la vulgarisation de l’outil de la programmation informatique est une solution pour
remédier à ce problème. C’est ainsi que dans le cadre de ce travail nous nous proposons
d’écrire un code de séparation par la méthode polynomiale par moindres carrés de Thurston et
Brown.
Dans ce travail, la méthode utilisée est la méthode gravimétrique, qui appartient aux méthodes
dites potentielles. De ce fait la plus grande partie des travaux antérieurs sur l’élaboration des
codes pour la séparation portaient soit sur la méthode des éléments finis (Ndougsa et al.,
2013) soit sur des anciens langages comme turbo-pascal 7.0 (Njandjock et al., 2003),Fortran
77 (Radhakrisha et al., 1990) et C++ (Njandjock et al., 2012). Dans cette optique nous allons
orienter notre étude un peu dans un autre style de langage plus récent tel que Matlab (Matrix
Laboratry) ceci dans le but de tester avec les données géophysiques de ma zone d’étude de
ressortir les différentes cartes existantes.

b. Problématique

Les questions de recherche qui ont motivé ce travail sont simplement de savoir :

 comment effectuer la séparation des anomalies en utilisant un code obtenu à partir du

langage de programmation Matlab ;
 quels types de données sont utilisés pour tester ce code ;

 2

 quelles sont les performances et le rendement de cette méthode vis-à-vis des autres
méthodes.

c. Objectifs

Pour ce faire, nous allons nous fixer certains objectifs pour résoudre cette problématique.

c.1. Objectif général

L’objectif général à atteindre est l’élaboration d’un code dans un environnement Matlab pour
une séparation régionale et résiduelle des anomalies par l’approximation polynomiale des
moindres carrés. Ce code nous permettra de faire des estimations sur des bases de données
géoréférenciées.

c.2. Objectifs spécifiques

Une analyse de ce travail nous conduit à mettre l’accent sur :

 l’importance des codes en géophysiques

 une présentation de Matlab
 l’écriture et l’implémentation du code en Matlab
 l’application aux données géophysiques
 l’interprétation et la corrélation avec les résultats antérieurs.

d. Plan

Afin de mener cette étude à terme, nous allons morceler le travail en quatre grands chapitres
dont la quintessence est brièvement indiquée suivant la structuration ci-dessous :

 le chapitre 1 permet de présenter globalement une revue de la littérature sur le rôle et

l’importance des codes en géophysique et ainsi que sur Matlab ;

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 3

 Le chapitre 2 tourne autour de la méthode gravimétrique et le formalisme mathématique

sur la séparation polynomiale par moindres carrés, à la programmation et à
l’implémentation du code.

 Le chapitre 3 nous conduit à appliquer les données sur le code en question, à présenter les
résultats obtenus dans l’optique de faire une analyse et une interprétation. Ce chapitre se
termine par une discussion qui permet de corréler les modèles obtenus avec un autre code
testé plutôt avec la séparation par la méthode des éléments finis et l’apport de
l’implication de ce code dans le système éducatif.

Une conclusion générale qui retrace les grandes lignes de ce travail ainsi que les perspectives
offertes clôtureront notre étude.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 4

CHAPITRE 1

REVUE DE LA LITTERATURE

Introduction

La méthode Thurston et Brown est une méthode polynomiale par moindres carrés pour
une séparation régionale et résiduelle ; cette séparation de l’anomalie régionale et résiduelle
est faite à partir d’un code écrit en Matlab. Pour cela, nous allons tout d’abord parler du
pourquoi des codes en géophysique en général et en particulier faire une présentation du
langage Matlab car c’est lui qui va nous permettre de séparer ces différentes anomalies, et
nous finirons par présenter l’importance des codes informatiques dans la réalisation des
travaux géophysiques.

1.1. Généralités sur les codes en géophysique

1.1.1. Définition et évolution du code

Un code est un texte qui représente les instructions qui doivent être exécutées par un
microprocesseur. Il est également l’élément de base des logiciels, écrit à partir d’une
application en un langage de programmation.

La conception et le développement des codes géophysiques présentent une étape très

importante pour l’analyse et le traitement des données géophysiques dans les environnements
informatiques. Ceci pour réduire considérablement la marge d’erreur de l’interpolateur et
d’effectuer ces manipulations en un temps relativement court. L’évolution simultanée de
 5

l’ordinateur et des logiciels sont dûs au fait qu’à chaque génération d’ordinateur, un type de
logiciel était utilisé ;

1.1.2. Rôle et importance des codes en géophysique

L’analyse et le traitement des données représentent une opération fastidieuse pour le

physicien en général et le géophysicien en particulier. En effet la taille des données et la
complexité des équations mathématiques rendent ce traitement difficile et nécessite un temps
de traitement énorme. Pour pallier à toutes ces difficultés, la programmation informatique
offre des facilités de calculs et de gains de temps intégrant ainsi l’utilisation d’un nombre
important de données.

Un programme informatique est une succession d’instructions exécutable par

l’ordinateur .Cependant, ordinateur ne sait manipuler que du binaire (c’est-à-dire 0 et 1) ; il
est donc nécessaire d’utiliser un langage de programmation adéquat pour écrire de façon
lisible. (Pillou, 2014).

D’une manière générale, le programme est un simple fichier texte (écrit avec un traitement de
texte ou un éditeur de texte), que l’on appelle fichier source. Le fichier source contient les
lignes de programmes que l’on appelle code source. Ce fichier source une fois terminé doit
être compilé.

Dans le domaine de la géophysique, le développement des codes informatiques et la

programmation facilitent le traitement des données dans plusieurs branches de la
géophysique.

En gravimétrie :

 Un code de séparation des anomalies gravimétriques et magnétiques utilisant la méthode

des éléments finis a été élaboré. (Ndougsa et al., 2013) ;
 un programme (MAPINFO 7.5) a été utilisé pour ré-digitaliser et créer de nouvelles cartes
d’anomalies aéromagnétiques résiduelles. (Ndougsa et al., 2012) ;
 un code écrit en C++ est rédigé et utilisé pour une étude sur la séparation de l’anomalie
régionale et résiduelles du champ de potentiel (Njandjock et al., 2012) ;

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 6

 Un code de cokrigeage rédigé puis compilé dans Matlab pour traiter les données en 2-D
(Marcotte, 1991) ;
 Un code écrit en Turbo Pascal 7.0 pour mener à bien l’ajustement du polynôme du champ
de potentiel des anomalies gravimétriques par la méthode analytique des moindres carrés
(Njandjock et al., 2003) ;

En sismique et sismologie :

 l’utilisation des logiciels tels que SAC (Seismic Analysic Code) développé au laboratoire
national de Lawrence Livermore pour l’analyse et le traitement des données sismiques.

Ainsi les codes en géophysique sont d’une importance capitale dans ce sens qu’ils facilitent la
manipulation d’un grand nombre de données en un temps relativement faible. De plus son
traitement des données est beaucoup plus aisé grâce aux logiciels de traitement qui intègre de
nombreux programmes qui exécutent des tâches précises. Dans le cadre de notre étude nous
allons élaborer un code en Matlab, ainsi une présentation de ce langage sera judicieuse avant
toute utilisation.

1.2. Présentation générale de Matlab

1.2.1. Démarrage de Matlab

Pour lancer l’exécution de Matlab on procède comme suit :

-sous Windows, il faut cliquer sur démarrage ensuite sous sur Matlab ;

- L’invite ‘ ’ de Matlab doit alors apparaitre, à la suite duquel on entrera les commandes ;

- Pour quitter Matlab, il faut cliquer sur la fonction ‹‹ Quit ›› ;

Quit.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 7

1.2.2. Interface de Matlab

Le nom Matlab provient de l’abréviation Matrix Laboratory qui est considéré comme
un langage de programmation adapté pour l’analyse et la résolution des problèmes
scientifiques. Il intègre à la fois des solutions de calcul, de visualisation et un environnement
de développement. Tous les objets définis dans Matlab le sont donc aux moyens de vecteurs
et de matrices/ tableaux de nombres (Michel, 2002). Matlab a de nombreux avantages par
rapport aux langages de programmation traditionnels (tels que C, C++, Pascal);Il est
spécialisé dans le domaine du calcul matriciel numérique.

L’interface Matlab se compose d’une fenêtre principale subdivisée en sous fenêtres suivant la
figure 1 ci-dessous:

Figure 1.1 : interface graphique de Matlab.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 8

Les principaux outils disponibles depuis cette interface (Ott, 2014) sont les suivants :

COMMAND HISTORY : historiques des commandes lancées depuis l’invite de

commande.

WORKSPACE : il liste des variables en mémoire, il permet également de parcourir

graphiquement le contenu des variables.

HELP BROWSER : un navigateur permettant de parcourir l’aide de Matlab. L’aide est

l’outil précieux pour trouver les fonctions et apprendre leur fonctionnement (notamment le
format des données à fournir en entrée ainsi que les valeurs renvoyées par la fonction).

COMMAND WINDOWS : invite de commande permettant de taper des instructions,

d’exécuter des fonctions Matlab et d’appeler des scripts dans un éditeur (voir figure 1.2).

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 9

Figure 1.2 : éditeur Matlab version 7.10.0 (R2010a).

Il existe deux modes de fonctionnement :

 Mode interactif : Matlab exécute les instructions au fur et à mesure qu’elles sont
données par l’usager.

 Mode exécutif : Matlab exécute ligne par ligne un programme Matlab (d’extension
.m) ; un programme Matlab (ou ‘m-file’ en anglais) est une suite d’instructions Matlab

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 10

écrite dans un éditeur de texte et sauvegardée dans un fichier avec une extension .m ; il existe
ainsi deux types de fichier .m : les scripts et les fonctions (Petrenko, 2004).

1.2.3. Fichiers exécutables : les scripts

Il arrive que l’on doit exécuter la même tâche plusieurs fois mais en changeant
seulement quelques paramètres. Une bonne façon de faire cela est de créer un fichier
exécutable encore appelé script. Ceux sont des fichiers contenant des suites d’instructions
Matlab dont le nom a comme extension. m ; par exemple «test1.m» (Petrenko , 2004).

Dans la fenêtre Commande, si l'on tape >> test1, les instructions contenues dans le fichier
test1.m seront exécutées par Matlab une par une. On peut créer des fichiers M à l'aide de
«Text Editor» ou dans la fenêtre Matlab en choisissant dans le menu "file" : new M -file.

1.2.4. Les fonctions sous-routines

On peut écrire des fonctions Matlab que l’on peut ensuite appeler depuis un script on les
appelle fonctions sous-routines (Michel, 2002).

On peut en distinguer trois types :

Les fonctions scalaires : elles s’appliquent essentiellement sur des scalaires. Mais elles
s’agissent élément à élément lorsqu’elles sont appliquées aux matrices (Georges , 2008).
Ces fonctions sont :
Sin asin exp abs
Cos acos log sqrt
tan atan rem round

Les fonctions vectorielles : Elles s’agissent sur des vecteurs (vecteurs lignes ou
colonnes), mais agissent sur une matrice ( par colonne pour générer un
vecteur ligne. (Voir table 1.1)

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 11

Table 1.1 : les différentes fonctions vectorielles

max permet de trouver le maximum d’un vecteur ou d’une matrice

Min permet de trouver le minimum d’un vecteur ou d’une matrice

mean (x,k) détermine la moyenne des éléments d’une matrice, vecteur
sort (x) donne une matrice ligne dont les éléments sont en ordre décroissant
Sort (y) donne une matrice ligne dont les éléments sont en ordre croissant
sum calcule la somme des éléments de la matrice x

On peut les faire opérer sur des matrices. Dans ce cas elles agissent sur chacune des colonnes
de la matrice. Pour obtenir une action ligne par ligne il suffit d'utiliser l'opérateur de
transposition. Ainsi si A est une matrice, max(A) renvoie un vecteur ligne dont chaque
composante correspond au maximum de chacune des colonnes de A, alors que max(max(A))
renvoie l'élément maximum de A.

Les fonctions matricielles : La puissance de MATLAB se trouve dans ces fonctions

matricielles. Voilà quelques-unes de ces fonctions listées dans la table 1.2 ci-dessous :

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 12

Table 1.2 : les différentes fonctions matricielles.

Eig valeurs et vecteurs propres

eye Identité
Zéro matrice de zéros

Ones matrice de uns

Triu Matrice triangulaire supérieure
Tril Matrice triangulaire inférieure
Size Matrice de Taille
Hilb Matrice de Hilbert

Par exemple: zeros (m, n) produit une matrice de avec tous les éléments nuls.

 Dans le cas du calcul matriciel:

Nous travaillons directement avec les éléments de matrices /vecteurs. Si A est de taille
et B de taille × , alors C = A*B effectue le produit matriciel de A et B:

= ∑

et si A est de taille × et inversible et B de taille × , alors X= A\B (division à gauche)

retourne l’unique matrice X de taille × solution de l’équation A*X= B .En particulier, la
solution du système d’équations linéaires

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 13

1.2.5. Les variables, les boucles et les instructions

Une variable est un emplacement en mémoire permettant de stocker provisoirement

une donnée. Une variable en programmation n’a pas la même signification qu’en
mathématique. On distingue plusieurs types de variables selon les données qu’elles servent à
stocker (nombre, caractère alphabétique, tableau, matrice, structure). Contrairement à d’autres
langages de programmation, sous Matlab le type de variable n’a pas besoin d’être spécifié,
Matlab infère le type d’une variable en fonction de la donnée que l’on y stocke et ces
variables sont créés lors de la première affectation (Ott, 2014). Sous Matlab, l’affectation se
fait à l’aide de l’opérateur « = » (A= A+1)

L’utilisation des boucles est le premier pas dans la programmation. En Matlab, les
boucles ‹‹ for ›› et ‹‹ while ›› sont très utilisées pour les processus itératifs.

La boucle ‹‹ for ›› est associée à une variable et exécute un processus plusieurs fois en prenant
à chaque fois une nouvelle valeur pour cette variable. Sa structure est la suivante :

for indice =

suite d’instruction

End

La boucle ‹‹ while ›› fonctionne de façon similaire; mais elle répète une suite d’instructions
jusqu’à ce qu’une condition soit satisfaite. Si la condition est fausse, le programme saute à
l’instruction qui suit End .Si la condition est vraie, le bloc d’instructions est exécuté ; sa
syntaxe est la suivante :

While expression booléenne

suite d’instructions

End

Le caractère % permet de spécifier un commentaire dans le code.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 14

1.2.6. Compilation dans Matlab

En Matlab les extensions se termine par ‘.m’ dans le nom fichier programme. En
mode programmation ou en mode commande, l’application exécute une série d’instructions.
Matlab tourne les données et produits les résultats escomptés (Gautrabd, 2013) .Au cours de
l’exécution, si un message d’erreur apparaît dans la commande Windows, il indique les lieux
où se trouvent les erreurs dans le fichier source. Si la compilation n’affiche aucune erreur
alors le programme fonctionne normalement et affiche les résultats dont on n’a besoin.

1.3. Travaux géophysiques utilisant la méthode Thurston et Brown

Agocs (1951) est le premier investigateur géophysicien à appliquer la méthode des

moindres carrées dans le problème de la séparation régionale et résiduelle. Simpson(1954)
quant à lui a adapté des polynômes d'ordre réduit aux données de la gravité par la méthode
des moindres carrés, en utilisant un modèle numérique.

Zeng (1989) procède par estimation du degré du polynôme afin d'évaluer l'anomalie régionale
présentée. La surface analytique peut être adaptée par des polynômes de divers degrés avec la
méthode des moindres carrés. Des exemples théoriques et sur le terrain montrent l'efficacité
de la méthode.

Pour Beltrao et al., (1991) Le champ régional est approché par un polynôme ajusté au champ
observé. Dans l’article écrite par Beltrao et al., (1991) ; il présente une méthode de séparation
résiduelle-régionale pour les données de gravité qui utilise une procédure robuste pour
déterminer les coefficients d'un polynôme ajustée aux observations.

Riad(1985) examine une approche pour déterminer les moindres carrés afin d’optimiser la
surface régionale qui, lorsqu’il est soustrait des données d’anomalies de Bouguer minimise la
distorsion de la composante résiduelle du champ.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 15

Lance (1982) montre que les surfaces polynomiales des moindres carrés peuvent être utilisées
dans les procédures de séparation d’anomalie pour les études de gravité. Une analyse détaillée
des effets précis de filtrage des polynômes des moindres carrés est prévue pour le cas
particulier des données régulièrement espacées. Les propriétés de réponse sont également
considérées comme des résidus de polynômes.

Conclusion

En définitive il a été question pour nous dans ce chapitre de faire une revue de la littérature
portant sur les codes géophysiques en général et de faire une présentation générale du
langage de programmation Matlab en particulier car notre code sera écrit à partir de ce
langage et pour pouvoir écrire un code il est important de connaitre sa syntaxe.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 16

CHAPITRE 2

METHODOLOGIE ET PROGRAMMATION

Introduction

Ce chapitre est consacré aux formalismes mathématiques conduisant à l’élaboration d’un

programme pour la séparation des anomalies régionales et résiduelles à partir de la méthode
des moindres carrés. De ce fait nous allons présenter les notions fondamentales du
magnétisme, les méthodes de séparation pour la résolution du problème, ainsi qu’une étape
de programmation en vue de faciliter l’interprétation.

2.1. Théorie de la méthode gravimétrique

2.1.1. Définition et fondement de base

La gravimétrie est l’une des disciplines fondamentales de la géophysique dont le but est
d’étudier les variations du champ de gravité pour en déduire la répartition des masses au sein
de la terre et ainsi sa structure (Dubois ; Diament ; Cogne., 2011).

 Lois de l’attraction universelle de Newton

Première loi de Newton

Cette loi qui est à l’origine de la gravité stipule que : « entre deux corps de masse et
s’exerce une force d’attraction dont l’intensité est proportionnelle aux masses et et,
inversement proportionnelle au carré de la distance entre les centres de masse des deux
corps. »

⃗ ⃗ (2.1)

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 17

Ou r est la distance séparant les centres de masses des deux corps

G est la constante de proportionnalité appelée constante de gravitation universelle ;

G= 6,67. N. .

Figure 2.1 : deux masses dans un référentiel cartésien (Giroux, 2008).

Deuxième loi de Newton

D’après cette loi, une masse se met en mouvement d’accélération lorsqu’elle est
soumise à l’action d’une force ; ceci se traduit par la relation :

(2.2)

En remplaçant par son expression dans l’équation (2.1), on a :

⃗⃗ (2.3)

En particulier si est la masse de la terre et R son rayon, l’accélération à la surface de la

terre est donnée :

⃗⃗⃗ = (2.4)

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 18

⃗⃗⃗ (2.5)

Où :

est la masse de la terre

⃗⃗ est l’accélération de la gravité et vaut 9,81m/

 Mesure de la gravité

Unité de la gravité : dans le système international l’unité de la pesanteur est le m/ , en unité

C.G.S elle est en cm/ .

En gravimétrie : en l’honneur de Galilée, on a nommé l’unité le « gal ».

1 Gal = 1cm/ = m/ ;

Pour apprécier les variations de gravité qui sont généralement très faibles, on utilise le
milliGal : 1mGal = Gal = m/

Appareil de mesure :

Il existe deux types de mesure de la gravité de pesanteur à savoir : les mesures absolues et les
mesures relatives :

- les mesures absolues

- les mesures relatives

2.1.1. Corrections et anomalies de la gravité

Pour rendre comparable entre elles les mesures de gravité effectuées à la surface de la terre,
on ramène celle-ci aux valeurs qu’elles auraient sur le géoïde ou tout simplement sur
l’ellipsoïde de référence. Cette opération est réalisée en une station par un ensemble de

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 19

corrections, puis en évaluant la différence entre la mesure corrigée et la valeur (théorique) de

g sur l’ellipsoïde de référence, on obtient l’anomalie de la gravité à la station considérée.

Les variations de gravité d’un point à un autre à la surface de la terre dépendent d’un certain
nombre de facteurs tels que l’altitude de la station, la latitude de la station liée à la forme
ellipsoïdale de la terre, la topographie de la surface de la terre et les hétérogénéités du sous-
sol. Ainsi les corrections appliquées sont :

-la correction de Faye

-la correction de plateau

-la correction topographique

-la correction de la latitude

–la correction luni-solaire

L’anomalie de Bouguer est la différence entre la mesure de la gravité ayant subi les
corrections (correction de Faye, topographique et de plateau) et la valeur théorique de la
gravité calculée en un point de l’ellipsoïde de référence (Naudy et al., 1965).

Elle constitue la réponse gravimétrique causée par les hétérogénéités de densité du sous-sol.
C’est donc à partir de cette anomalie que le traitement suivit de l’interprétation se fait. Elle est
donnée par la relation:

est la valeur observée de la pesanteur en chaque station

est la valeur théorique de la pesanteur sur l’ellipsoïde de référence

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 20

2.2. Méthodes de séparation régionale et résiduelle

Toute analyse structurale du sous-sol ne peut se faire qu’à partir d’une résiduelle, il faut donc
séparer la régionale de la résiduelle. Pour cela il existe plusieurs méthodes à savoir :

2.2.1. Méthode graphique

Cette méthode nécessite à priori des hypothèses basées par exemple sur la connaissance
géologique de la région, permettant de fixer l’allure générale de la régionale, qui présente
généralement une forme régulière. La résiduelle est obtenue par soustraction de la régionale
sur la Bouguer soit graphiquement, soit analytiquement. Cette méthode peut être réalisée de
trois manières différentes (Schoeffler ,1975) :

a. On procède au lissage de la carte de Bouguer

b. On considère des profils gravifiques pris perpendiculairement à la direction générale de
Bouguer puis, on procède au lissage graphique des profils
c. On trace, puis on lisse des profils obtenus selon une grille d’un pas bien déterminer.

La méthode de séparation par lissage graphique n’est pas indiquée :

-lorsque la régionale présente une pente très forte.

-lorsque les anomalies résiduelles sont de grandes dimensions.

Dans ce cas il est conseillé d’utiliser les méthodes analytiques.

2.2.2. Méthodes analytiques

Il existe plusieurs méthodes analytiques telles que :

-la méthode des dérivées secondes

-le prolongement vers le haut ou vers le bas

-la méthode polynomiale

Nous allons nous attarder sur la méthode polynomiale

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 21

2.2.3. Méthode de séparation analytique polynomiale

Le principe consiste à construire un polynôme d’ordre aussi élevé qu’il conviendra, qui
génère une surface analytique qu’on adapte à une surface expérimentale par la méthode des
moindres carrés. Cette surface analytique représente la régionale. Le problème ne trouve de
solution que sur des surfaces réduites (Ndougsa ,2004).

Pour résoudre le problème nous allons écrire un programme qui permet d’ajuster un polynôme
d’ordre déterminé correspondant aux anomalies d’un champ potentiel et de calculer les
anomalies résiduelles pour un grand nombre de point(x, y) ou l’anomalie de Bouguer vaut .
Ce programme génère un polynôme d’ordre m, établit et résout le système d’équations
linéaires. L’ajustement des anomalies au polynôme d’ordre m est réalisé par la méthode des
moindres carrés. De ce fait le choix de l’ordre du polynôme représentant la régionale n’est pas
arbitraire, il dépend de la nature de la résiduelle qui concerne les anomalies significatives à
interpréter ; plus la régionale est lissée plus elle permet d’éliminer les traits de profondeurs.

Dans le cadre de ce travail nous allons utiliser cette technique de séparation car c’est elle qui
nous permettra d’écrire le code.

2.3. Approche de la méthode de séparation polynomiale par moindres carrés

2.3.1. Intérêt de la séparation

Le tracé de l’anomalie de Bouguer fait apparaitre les effets des contrastes de densité
profonds et peu profonds, étendus et locaux ; ainsi pour une interprétation plus fiable, il faut
séparer :

 L’effet de contrastes de densité profonds et étendus qu’on appelle « anomalie

régionale »
 De l’effet des contrastes de densité peu profonds et moins étendus appelé « anomalie
résiduelle ».

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 22

Cette transformation est appliquée ici comme technique d’interpolation polynomiale pour
permettre un lissage des anomalies dans la zone afin d’avoir des sources qui se rapprochent
distinctement de la lithologie locale.

2.3.2. La méthode de séparation par moindres carrés

la surface analytique sera notée Z et ne dépendra que des coordonnées des positions à chaque
point de la station et est de la forme (Thurston et al., 2012) :

( ) ∑∑ ( (

Ou ( est un ensemble de polynôme de degré N et sont des coefficients à

déterminer.

L’équation (2.9) peut encore se mettre sous la forme :

( ) ∑∑ ( (

Avec :

( =
(
Et

Pour N fixé on a :

( (
M= ; le nombre de terme du polynôme
Soit :

( ( Coefficients à déterminer

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 23

On peut simplifier les équations (2.9) et (2.10) en écrivant :

( ) ∑ ( (2.11)

Avec ( =

Et m= ( (

Soit ( ∑ ∑ ( , l’écart entre les points homologues

des surfaces expérimentales et analytiques respectivement, le nombre de points ou de
stations P( ou l’anomalie de Bouguer B( observée est connue. Pour ajuster les
surfaces, on doit rendre minimale l’écart quadratique :

∑( (

La méthode des moindres carrés calcule les coefficients des polynômes tels que la somme des
termes des erreurs au carré (résiduelle) est minimale et se présente sous la forme :

∑ (

L’équation (2.13) conduit à :

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 24

∑ (

( ( ( )
Cependant (

La relation (2.14) devient :

∑ ( ( ( (

Avec k= ( (

En remplaçant ( par sa valeur et en séparant les deux membres de l’équation (9) on

a:

( (

∑ ∑ ( ( ∑ ( ( (

L’équation (2.16) génère un système de m équations à m inconnus (2.17) suivant :

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 25

( (
∑ ∑ ( ) ( ) ∑ ( ) ( )

( (
∑ ∑ ( ) ( ) ∑ ( ) ( )

( (

∑ ∑ ( ) ( ) ∑ ( ) ( )
{

Le système (2.17) suivant peut s’écrire sous la forme matricielle (2.18) suivante :

∑ ∑ ∑ ∑ ∑
∑ ∑ ∑ ∑ ∑
∑ ∑ ∑ ∑ ∑
=

[∑ ∑ ∑ ∑ ][ ] [∑ ]

Ou ( est requis dans le calcul des coefficients des polynômes par la méthode
des moindres carrés.
Les inconnus étant les coefficients du polynôme ( d’ordre N ; sa détermination
nécessite la résolution du système (2.17) de m équations à m inconnus et cela sera effectué
suivant un programme écrit en Matlab selon le degré du polynôme.

Une fois les coefficients déterminés, on calcule l’anomalie régionale analytique

( ( et on déduit la résiduelle :

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 26

( ( ( (2.19)

2.3. Etape de la programmation

2.3.1. Variables utilisées

2.3.1.1. Types de variables

 Variables réelles : x, y, M
 Variables entières : i, j, l, k, N
 Taille : No
 Matrices : X, Y, B, A, E, Re, Ri

2.3.1.2. Données

X et Y sont les coordonnées de position de la station

B est l’anomalie de Bouguer observé

No représente le nombre de stations

M l’ordre du polynôme

Reg représente la surface régionale

Ris est la résiduelle

2.3.2. Etape du programme

ETAPE 1 : DECLARATIONS

-types des variables

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 27

-tableaux

ETAPE 2 : LECTURE DU NOMBRE DE STATIONS No

ETAPE 3 : LECTURE DES DONNEES

X(i), Y(i) et B(i)

ETAPE 5 : LECTURE DE L’ORDRE DE LA SURFACE « N »

ETAPE 6 : CALCUL DES COEFFICIENTS DE L’EQUATION

Pour M allant de 1 à (N+1) (N+2)/2

ETAPE 7 : REMPLISSAGE DE LA MATRICE « A »

A=∑ ( (

ETAPE 8 : REMPLISSAGE DE LA MATRICE DU SECOND MENBRE « E »

∑ ( (

ETAPE 9 : RESOLUTION DES MATRICES « A » ET « E »

ETAPE 10 : CALCUL DES COEFFICIENTS DU POLYNOME « C »

C = E/ A

ETAPE 11 : CALCUL DE L’ANOMALIE REGIONALE « Res »

Grâce à la relation (2.11)

ETAPE 9 : CALCUL DE L’ANOMALIE RESIDUELLE « Ris »

Grâce à la relation (2.19)

ETAPE 10 : EDITION DES RESULTATS

X(i) ,Y(i) , B(i), C , Res (i) et Ris (i)

Pour i de 1 jusqu’à No.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 28

2.4.3. Programme

Ce programme de séparation des anomalies régionales et résiduelles se trouve dans la partie

annexe.

Conclusion

Dans ce chapitre, nous avons présenté la méthode qui nous a permis de ressortir un
algorithme qui sera utilisée pour la mise sur pied de notre code écrit en langage Matlab, il a
été question dans ce programme de faire une séparation des anomalies régionales et
résiduelles à partir de l’anomalie de Bouguer .Ce programme permettra grâce à l’ordre du
polynôme qu’on fixera de ressortir les différentes anomalies pour chaque station ,de les
interpréter et de connaitre le pourquoi du choix de l’ordre du polynôme et du choix de la
régionale. Tout ceci sera vue dans le chapitre suivant.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 29

CHAPITRE 3

RESULTATS ET DISCUSSION

Introduction

Dans ce chapitre, nous allons utiliser toutes les informations sur la méthode présentée aux
chapitre deux pour tester notre code se trouvant en Annexe 2 avec des données relatives à une
base. Pour cela, nous procédons successivement à l’analyse des données de la carte de
Bouguer, à la séparation des anomalies de Bouguer, à l’interprétation qualitative des
anomalies régionales et résiduelles observées ,faire une comparaison et donner l’apport de ce
code dans le système éducatif.

3.1. Données gravimétriques

 Origine des données

Les données utilisées dans le cadre de ce travail proviennent de la carte d’anomalie de

Bouguer obtenue par Zanga (2015); (Voir fig.3.1 ci-dessous).

Les données gravimétriques utilisées dans la présente étude sont localisées dans les zones
situées entre l’Est –Cameroun et Ouest de la république centrafricaine et couvrent la zone qui
s’étend du Sud au Nord entre les latitudes 2°15’ et 5°15’ N et de l’Ouest à l’Est entre les
longitudes 13°44’ et 16°10’ E.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 30

Figure 3.1 : Carte de l’anomalie de Bouguer (Zanga A., 2015).

 Manipulation des données gravimétriques

Après digitalisation de la carte de Bouguer de Zanga (2015) (fig. 3.1), nous obtenons trois
cent quarante-neuf (349) données d’anomalies de Bouguer géoréférencées contenues dans le
tableau en annexe 1.

Par un programme informatique Matlab on a procédé à une séparation régionale et résiduelle

de l’anomalie selon un ordre établi ; par la suite nous avons exporté le fichier de données
obtenu dans le logiciel Oasis montaj pour l’établissement des cartes d’anomalies dans le but
de faire une interprétation qualitative.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 31

3.2. Etablissement, analyse et traitement des cartes d’anomalies

3.2.1. Cartes de l’anomalie de Bouguer

Les données gravimétriques exploitées dans ce travail sont extraites de la carte d’anomalie de
Bouguer de Zanga (2015) qui ont permis à l’établissement d’une nouvelle carte de l’anomalie
de Bouguer (fig. 3.2). Cette carte permet d’observer la répartition des sources d’anomalies au
Cameroun et en république centrafricaine (RCA) (Zanga ,2015).

Les bandes d’isogales observées sont toutes négatives et espacées en moyenne de 5 mGals ;
avec des valeurs variant entre -82.7 à -34.7 mGals. Cette carte reflète des excès et des défauts
de masses relatives à l’hétérogénéité de la densité de roches dans la formation du socle de la
région.

Une analyse de cette carte nous montre des domaines d’anomalies relativement lourdes et des
domaines d’anomalies légères, les deux étant séparés par des gradients horizontaux plus ou
moins importants caractérisés par le resserrement des lignes isoanomales.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 32

Figure 3.2 : Carte de l’anomalie de Bouguer obtenue à partir du logiciel oasis montaj.

3.2.1.1. Les zones d’anomalies lourdes

Elles sont associées soit à des intrusions de roches lourdes dans la croûte, soit à un relèvement
du socle, soit à des roches sédimentaires denses.

L’analyse de cette carte présente un axe lourd P1 d’orientation WNW-ESE d’amplitude -34.7
mGals qui passe par les localités de Batouri à Kenzou ou elle s’affaiblit sous forme d’une
courbe fermée.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 33

ces anomalies lourdes ont une amplitude minimales de -38.3 mGals et sont représentées par
des isoanomales allongées dont l’une fermée sous forme de courbe en ovoïde. Les anomalies
moyennement lourdes entourent celles lourdes avec les amplitudes maximales de -38.3mGals
aux environs des localités de Bertoua-Berberati.

3.2.1.2. Zones d’anomalies légères

Les anomalies légères ou négatives sont associées à des intrusions de roches légères dans la
croûte, soit à un épaississement ou à un effondrement du socle.

On distingue dans cette zone d’étude d’important boulevard anomalies légères et

moyennement légères.

-Un axe N1 d’orientation (NE-SW) d’amplitude minimale -82.7 mGals qui s’allonge à l’Est
de Yokadouma marqué par des iso- valeurs allant de -80 mGals à -75 mGals , entouré par des
anomalies moyennement légère -72.0 mGals.

-Un axe N2 d’orientation (N-S) d’amplitude minimale de -78.6 mGals qui s’étend un peu plus
vers le Sud de la carte, entouré aussi par des anomalies moyennement légères. Ce domaine est
caractérisé par des isoanomales sous forme de lobules.

-un axe N3 d’amplitude minimale -82.7 mGals qui s’allonge vers l’Est de la carte.

3.2.1.3. Zones de gradients

Les secteurs d’anomalies gravimétriques (lourdes et légères) mis en évidence sont séparées
par des isoanomales quasi-linéaires caractérisant des gradients .Ces gradients représentent des
discontinuités gravimétriques matérialisant des structures allant des faibles aux grandes
profondeurs. Une analyse profonde de la carte de Bouguer présente une zone de gradients
séparant les zones d’anomalies positives et négatives. La carte de Bouguer nous présente un
paquetage des lignes resserrées importantes au-delà de Ngoura et un peu moins dans les
environs de Mbang.

Il est annoté que cette carte d’anomalie de Bouguer nous montre que la plupart de ces villes
sont reparties au centre de la carte dans les zones lourdes et moyennement lourdes.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 34

3.2.2. Séparation des anomalies de Bouguer

3.2.2.1. Calcul des anomalies régionales et résiduelles

La séparation des anomalies de Bouguer a été effectuée par une méthode analytique la
méthode polynomiale. De ce fait un code en Matlab a été élaboré qui génère un polynôme
d’ordre m, lequel établit et résout un système d’équations linéaires de m équations à m
inconnus pour la détermination des coefficients du polynôme, puis donne l’anomalie
régionale et résiduelle en chaque point où l’anomalie de Bouguer est définie. L’ajustement
des anomalies observées au polynôme d’ordre ‘’m’’ est réalisée par la méthode des moindres
carrés. La régionale est représentée par une fonction analytique qui peut être représentée par
une série de puissance (Baranov, 1954).Ainsi nous avons :

 Pour la surface du ordre (m=1), le polynôme est donné par la relation :

Reg ( )=

Les coefficients calculés sont obtenus dans la table 3.1 ci-dessous ;

Table 3.1 : valeurs des coefficients calculés à l’ordre 1 à partir du code Matlab

Coefficients Valeurs
-31,39397933
-0,022337827

0,025803118

 Pour la surface du ordre (m=2), le polynôme est donné par la relation :

Reg (x,y ) = =

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 35

Table 3.2 : valeurs des coefficients calculés à l’ordre 2.

coefficients Valeurs
-1848,331207

1,665459266

1,97284252

-0,00043972

-0,00049994

-0,00130677

 Pour la surface du ordre (m=3), le polynôme est donné par la relation :

Reg (x, y) =

Les coefficients calculés à l’aide du code Matlab sont les suivants :

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 36

Table 3.3 : valeurs des coefficients calculés à l’ordre 3.

coefficients Valeurs
-17607,92169

29,50393158

3,590787739

-0,015862291

-0,010550998

0,015507358

2,77574E-06

3,64529E-06

-2,56867E-06

-1,01793E-05

 Pour la surface du ordre (m=4), le polynôme est donné par la relation :

Reg (x, y) = =

Avec comme valeurs des coefficients celle de la table 3.4 ;

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 37

Table 3.4 : valeurs des coefficients calculés à l’ordre 4.

coefficients Valeurs
-163400,5127

333,3854432

164,7795389

-0,236530808

-0,394440887
0,205286534

6,69934E-05

0,000269431

-0,000143317

-0,000129794

-5,72113E-09

-5,80433E-08

2,84892E-08

3,80082E-08

3,37342E-08

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 38

Les valeurs des anomalies régionales et résiduelles obtenues aux différents points ( pour
les polynômes de degré 1, 2,3 et 4 sont contenus dans les tableaux en Annexe.

3.2.2.2. Les cartes d’anomalies séparées et leur analyse

 la régionale et la résiduelle à l’ordre 1

La régionale

Pour m=1 le polynôme Reg ( ) est représenté par une surface du degré. La carte
d’anomalie régionale obtenue (fig. 3.3a) présente des iso- anomales de direction générale
constante d’orientation ENE-WSW. Cette régionale ne présente pratiquement pas de trait
commun avec la Bouguer ; il en ressort de cela que la régionale est trop lissée.

De plus nous observons une zone d’anomalie lourde et moyennement lourde d’orientation
NW-SW s’étirant de Batouri vers le NW de la carte ; l’anomalie lourde a une amplitude de -
52.8 mGals et est entourée de celle moyennement lourde d’amplitude -53.8 mGals et deux
zones d’anomalies légères : un axe d’orientation N-S d’amplitude -63.5 mGals et un axe
d’orientation NW-SE.

La résiduelle

Sur la carte d’anomalie résiduelle correspondante (fig. 3.3.b), un examen visuel de la carte
permet de constater que les valeurs des isogales sont composées des portions positives et
négatives variant entre – 25.1 mGals à 21.3 mGals. Elle présente trois axes d’anomalies
légères et un axe d’anomalie lourde, les anomalies légères et moyennement légères ont des
bandes d’isogales négatives alors que les anomalies lourdes ont des bandes d’isogales
positives , les deux étant séparées par une zone de gradients plus intense au-delà de Ngoura et
moins intense en deçà de Mbang. Nous constatons que cette anomalie résiduelle possède des
traits communs avec la Bouguer.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 39

(a)

(b)

Figure 3.3 : Carte régionale et résiduelle d’ordre 1 de la zone d’étude : (a) carte régionale ;
(b) carte résiduelle.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 40

 la régionale et la résiduelle à l’ordre 2

La régionale

Pour m=2, le polynôme Reg (x, y) est représenté par une surface du second degré (fig. 3.4.c).
Les bandes d’isogales vont de -90.4 mGals à -46.7 mGals, la carte de l’anomalie régionale
obtenue présente une zone d’anomalie lourde au centre de la carte d’orientation E-W
s’étalant de Gari Gombo à Batouri d’amplitude -46.4 mGals se présentant sous la forme d’un
ovoïde ; cette anomalie lourde est entourée des anomalies moyennement lourdes d’amplitude
-48.5 mGals ; deux zones d’anomalies légères : un axe dirigé E-W et un axe de
direction NNE-SSW s’étirant vers Carnot. On observe une zone de gradients moins intense
entre -66 à -50.6 mGals se trouvant entre la zone d’anomalie lourde et la zone d’anomalie
légère.

La résiduelle

L’analyse de cette carte (fig. 3.4.d) nous montre qu’elle possède les bandes d’isogales allant
de -14.7 mGals à 14.8 mGals, de plus nous avons trois zones d’anomalies lourdes , et
et trois anomalies légères , , ;

-un axe d’orientation ESE-WNW, comportant des anomalies lourdes dans la région de
Ngoura et Kenzou se présentant sous forme d’ovoïde d’amplitude 14.8 mGals et est entouré
par des anomalies moyennement lourdes se trouvant dans les régions de Bertoua, Batouri et
Kenzou.

-un axe ne comportant que des anomalies lourdes dans la région de Berberati.

-un axe d’orientation E-W présentant des anomalies lourdes et moyennement lourdes
allongées vers le sud de la carte.

- un axe se trouvant dans une zone légère limitée au Nord par une iso-anomale de valeur 0
mGals qui s’étend au Sud de Yokadouma avec un minimum de -10 mGals.

-un axe d’anomalie légère d’amplitude -14.7 mGals et un axe présentant une anomalie
moyennement lourde d’amplitude -11.6 mGals.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 41

(c)

(d)

Figure 3.4 : Carte régionale et résiduelle d’ordre 2 de la zone d’étude : (c) carte régionale ;
(d) carte résiduelle.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 42

 la régionale et la résiduelle à l’ordre 3

La régionale

 Le polynôme Reg (x,y) pour m=3 représente par une surface de degré 3 (fig . 3.5.e) ,cette
carte d’anomalie régionale présente une zone d’anomalie lourde d’axe située NW de la
carte d’amplitude – 41.3 mGals localisée dans les régions de Batouri ,Ngoura et Kenzou et
une zone d’anomalie moyennement lourdes d’amplitudes -44.7 mGals se trouvant dans les
régions de Bertoua et Berberati.
 Les zones d’anomalies légères sont rencontrées dans chaque coin de la carte et présente
trois secteurs négatifs majeurs. Un axe de direction N-S, un axe de direction ESE-
WNW et un axe de direction S-E mettant en évidence les zones d’affaissement du
socle. Cette carte diffère de la Bouguer au niveau de la forme de l’anomalie lourde qui est
plus lissée dans le cas de la régionale. Elle présente une gamme d’isogales allant de -89.3
mGals à -41.3 mGals.
 une zone de resserrement des isoanomales accentuées dans la région de Carnot et moins
accentuées dans les régions telles que Mbang, Yokadouma et Nola qui se trouve à la
frontière entre l’Est-Cameroun et la république Centrafricaine mettant en évidence la zone
de gradients.

La résiduelle

La résiduelle pour m=3 (fig. 3.5.f) est différente de celle m=2. Cette résiduelle se détache un
peu plus de la Bouguer dans la mesure où nous observons les zones d’anomalies lourdes qui
se séparent et se répartissent dans toute la carte. Les bandes d’isogales sont constitués des
portions positives et négatives s’étalant entre -7.0 à 9.6 mGals. De part et d’autre de
l’anomalie lourde, nous avons aussi des couloirs des anomalies légères, chacune entouré des
isogales de valeur zéro en forme de boucle.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 43

(e)

(f )

Figure 3.5 : Carte régionale et résiduelle d’ordre 3 de la zone d’étude : (e) carte régionale ;
(f) carte résiduelle.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 44

 la régionale et la résiduelle à l’ordre 4

La régionale

Le polynôme Reg (x,y) est représenté par une courbe de degré 4. La carte d’anomalie
régionale (fig. 3.6.g) présente des similitudes avec la régionale d’ordre 3. Les régions telles
que Batouri, Ngoura et Kenzou sont situées dans une zone relativement lourde. Les lignes
isoanomales d’orientation E-W présente un gradient prononcé. Les bandes d’isogales s’étalent
de -87.5 mGals à -39.3 mGals

La résiduelle

La carte d’anomalie résiduelle (fig. 3.6.h) montre également des similitudes avec la carte
d’anomalie résiduelle d’ordre 3 ; et on constate que dans le cas de l’anomalie résiduelle
d’ordre 4 les amplitudes sont très faibles.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 45

(g)

(h)

Figure 3.6 : Carte régionale et résiduelle d’ordre 4 de la zone d’étude : (g) carte régionale ;
(h) carte résiduelle.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 46

3.2.2.3. Choix de la régionale

Le choix de la régionale c’est-à-dire de l’ordre du polynôme représentant la surface analytique

recherchée n’est pas arbitraire (Ndougsa, 2004). De ce choix dépend la nature de la résiduelle
qui permet d’interpréter les anomalies significatives.

La carte d’anomalie résiduelle traduit essentiellement les sources d’anomalies situées à faible
profondeur. De ce fait, plus la surface analytique régionale est lissée, plus il y a des chances
d’éliminer les traits profonds et ceux-là seulement. D’un autre côté, plus la régionale est de
configuration compliquée (polynôme d’ordre élevé), plus les traits rapprochés de la surface
sont éliminés ; ces traits seront davantage mis en évidence au niveau de la résiduelle.

Nous constatons que la régionale d’ordre 3 et 4 sont similaires ainsi que les résiduelles. De ce
fait, la séparation des anomalies de Bouguer ont conduit à choisir pour la régionale une
surface d’ordre 3, par conséquent, la résiduelle d’ordre 3 est la mieux indiquée pour rendre
compte des hypothèses géologiques (les contrastes de densité à l’intérieur de la croûte).

3.2.3. Synthèse sur l’analyse des cartes régionales et résiduelles d’ordre 3

Elle consiste à faire une étude comparative qualitative de ces cartes obtenues par la séparation
polynomiale par moindres carrés. Elle nous montre que :

 Sur la carte régionale les anomalies sont constituées seulement des bandes d’isogales
négatives ; de plus cette carte présente une zone d’anomalie lourde et des zones
d’anomalies légères ; elle ne possède pas de couloir entouré d’isogales de valeur zéro en
forme de boucle.

 Sur la carte résiduelle, ces bandes d’isogales sont composées des portions négatives et
positives. Nous avons la présence de l’existence des lignes des isogales de valeurs 0
mGals.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 47

3.3. Discussion et validité des résultats

Les résultats obtenus au cours de l’interprétation gravimétrique dépendent d’une part de la

technique utilisée pour le traitement des données et d’autre part des contraintes apportées dans
le choix des paramètres physiques intervenant dans l'élaboration des modèles. Dans cette
partie, il est question de justifier la validité de l’ensemble des résultats tels que l’établissement
et l’analyse des cartes gravimétriques que nous avons produit dans ce travail ; donc de ce fait
la fiabilité des conclusions dépende de la précision avec laquelle les résultats ont été obtenus,
ainsi une étude comparative sera effectuée avec les résultats obtenus à celle de Zanga (2015).

Les différentes cartes obtenues par la méthode de séparation des éléments finis par Zanga
(2015) nous présentent les résultats suivants :

 la carte de l’anomalie régionale (fig. 3.7.i) présente une anomalie négative avec une
amplitude des isogales variant entre -80 et – 50 mGals. Elle présente un maximum
d’environ -50 mGals dans la ville de Nola en république centrafricaine. La configuration
et la disposition de ces isogales sont orientés E-W et sont beaucoup plus concentrées dans
la région de Carnot.
 La carte d’anomalie résiduelle (fig. 3.7.j) quant à elle révèle une vaste zone lourde de
direction E-W qui culmine 30 mGals au Nord au voisinage de Ngoura au Cameroun ; au
Sud, une zone d’anomalie légère limitée au Nord par une iso-anomale de valeur 0 mGals
allant jusqu’à Yokadouma.
 Un paquetage des isogales avec un gradient important autour de l’anomalie de Ngoura et
relativement légère vers la zone de Carnot en RCA (fig. 3.7).
La comparaison entre la méthode polynomiale de séparation par moindres carrés et la
méthode de séparation par les éléments finis (Zanga A., 2015) a permis de constater une
similitude des résultats entre les deux méthodes pour la régionale et la résiduelle à l’ordre
1 ;mais n’est pas le cas pour l’ordre 3 car plus on augmente l’ordre du polynôme plus les
ressemblances diminuent.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 48

(i)

(j)

Figure 3.7 : Carte des anomalies régionales et résiduelle faites par Zanga A., 2015 : (i) carte
régionale ; (ii) carte résiduelle.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 49

Par ailleurs d’autres méthodes de séparation (Gauss-Newton et les Ondelettes) peuvent

s’avérer être proche de cette méthode analytique lors de l’analyse des cartes obtenus.

3.4. Implication dans le système pédagogique

Cette étude permet de montrer dans quelle mesure on parvient à développer chez l’apprenant
certaines habiletés telles que :

 L’utilisation de la programmation informatique dans les autres branches de la science ;

 Les zones potentielles où on peut détecter des anomalies naturelles ceci permet à
l’élève de mieux localiser les endroits où on observe des failles lors des travaux pratiques
sur le terrain;
 L’inventaire des environnements pouvant abriter les ressources naturelles comme : le fer,
la bauxite, l’or et le pétrole ;
 l’approche utilisée permet de montrer l’interconnexion entre la mathématique,
l’informatique et les géosciences en général et particulièrement en géophysique qui est
appliquée dans le cadre de l’inventaire des environnements susceptibles d’abriter les
métaux qui sont des ressources naturelles utiles à l’homme ;ces métaux sont utiles dans la
fabrication des alliages tels que les téléphones ;
 de plus dans le plan social elle permet la création des emplois ;

Conclusion

L’analyse des cartes d’anomalies régionales et résiduelles met en évidence les grands traits
gravimétriques de la zone étudiée. Il en ressort que le choix de la régionale c’est-à-dire de
l’ordre du polynôme est d’une importance capitale car c’est cette régionale qui nous permettra
d’établir la corrélation et les divergences entre les deux méthodes de séparation. Il s’en suit
que l’élaboration de ce code dans le domaine de la géophysique est une importance capitale
pour le système éducatif.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 50

CONCLUSION GENERALE

La bonne séparation des anomalies est celle qui permet de bien isoler les anomalies dont les
sources dues aux structures superficielles et/ou locales de celles dont les sources sont
lointaines et/ou profondes. De ce fait nous avions opté sur celle des techniques analytiques car
la méthode graphique lorsque les résiduelles sont dans des courbes fermées et la régionale
présentant des irrégularités compliquées, cette technique est défaillante. L’objectif majeur de
ce travail était de développer un programme dans un environnement Matlab pour une
séparation régionale et résiduelle des anomalies par l’approximation polynomiale des
moindres carrés. Pour atteindre ces objectifs, nous avons procéder comme suit :

 Tout d’abord présenter la syntaxe de Matlab car c’est elle qui est à la base de tout ;
 Par la suite ressortir les formalismes mathématiques qui nous permettront lors de
l’écriture ;
 Ecriture et l’implémentation du code en Matlab ; ce programme résout un système
linéaire de m équations à m inconnus à deux variables et donne les valeurs des
coefficients, la régionale et la résiduelle pour chaque point d’anomalie correspondant ;
 Application aux données gravimétriques
 Pour finir une interprétation et une corrélation avec les cartes obtenues par Zanga (2015)
avec la séparation par la méthode des éléments finis ; cette interprétation nous a permis de
mettre en évidence les zones d’anomalies lourdes et les zones d’anomalies légères les
deux étant séparés par une zone de gradient qui se matérialise par le resserrement des
isoanomales.

Les résultats obtenus à partir de ce code Matlab étant importants et intéressants; nous ne
pouvons-nous limiter à cela autrement dit, en vue d’améliorer la compréhension de ce code
certaines lignes nous ont paru un peu complexe en simplifier ainsi nous proposons également

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 51

pour poursuivre cela de revoir la matrice du problème qui est une matrice de Hilbert ; cette
matrice est très mal conditionnée c’est-à-dire qu’elle souffre d’une instabilité numérique due à
la nécessité de l’inverser pour pouvoir déterminer les coefficients car ceci conduit à une perte
de précision due aux erreurs d’arrondi.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 52

BIBLIOGRAPHIES

-Agocs W. B., (1951). Least squares residual anomaly determination. Geophysics. Vol XVI-
4. P. 686-696;

-Baranov V., (1954). Sur une méthode analytique de calcul de l’anomalie régionale
Géophysique Prospecting. Vol II-3. P. 203-226;

-Beltrao J. F., Silva J.B.et Costa J. C., (1991).Robust polynomial fitting method for
regional gravity estimation Geophysics. VOL 56-No I. P. 80-89 ;

-Ciarlet P., et Luneville E., Notes introductives à Matlab ENSTA-UMA P. 14 ;

-Dubois J., Diament M., Cogne J.P., (2011). Géophysiques notes de cours P.7 ;

-Giroux B. et Chouteau M., (2008).Notes de cours Géophysique Appliquée I, Gravimétrie

Laboratoire de géophysique appliquée Automne P.2 ;

-Jeffrey B.T. et Brown R.J., (1992).The filtering characteristics of least-squares polynomial

approximation for regional/residual separation Canadian journal of exploration geophysics.
Vol 28-No II. P. 71-80;

-Gautrabd J.D., (1996).Comment compiler un programme Matlab ;

-Georges Rodriguez G., (2008).Introduction à Matlab, Département COMELEC P. 5-8 ;

- Lance J.O., (1982). Frequency domain analysis of least-square polynomial surfaces with
application to gravity data in the Pedregosa basin area .These Doct. PhD.Univ.Texas à El
Paso.

-Li Y. et Oldenburg D., (1998). Separation of regional and residual magnetic field data;

-Marcotte D., (1991).Cokriging with Matlab. Computers & Geosciences, vol.17.No 9

P.1265-1280;

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 53

-Michel E., et Yverdon L.B., (2002).Introduction au logiciel Matlab. Ecole d’ingénieurs du

Canton de Vaud-Département d’électricité et d’informatique institut d’automatisation
industrielle (iAi) P.4;

-Naudy H. et Newman R., (1965).sur la définition de l’anomalie de la Bouguer .Geophysical

propecting .vol.XIII P. 1-11 ;

-Ndougsa Mbarga T., (2004). Etude géophysique par la méthode gravimétrique des
structures profondes et superficielles de la région de Mamfe. Thèse Doct PhD Univ.Yaoundé
I, Cameroun P 148-151;

-Ndougsa Mbarga T., Yufenyiu Layu D., Yene Atangana Q. et Tabod Tabod C., (2012).
Delineation of the northern limit of the Congo craton based on spectral analysis and 2.5D
modeling of aeromagnetic data in the Akonolinga-Mbama area, Cameroon;

-Ndougsa Mbarga T.,and Bikoro-Bi-Alou M., (2013). Filtering of gravity and magnetic
anomalies using the infinite element approach. Journal of Indian Geophysical Union.17(2).P.
155-166 ;

-Njandjock Nouck. P., Oyoa. V., Ndougsa Mbarga. T. et Bisso. D., (2012). C++ Code for
separation of regional and residual anomaly of potential field greener Journal of physical
sciences ISSN: 2276-7851. VOL II-No IV. P.120-130;

-Njandjock N.P., H. L. Kandé, Manguelle Dicoum.E.,Tabod. C.T., Ndougsa M.T. et Jean

M., (2003). A turbo Pascal 7.0 program to fit a polynomial
of any order to potential field anomalies based on the analytic least squares method.African
Journal of Science and Technology (AJST) Science and Engineering Series.Vol IV- No II.pp.
1-4;

- Ott L., (2015). Initiation à Matlab 1, URECA P.1-4 ;

-Petrenko A., (2004).Initiation à Matlab, Université Aix Marseille II P. 4-14;

-Pillou J. F., (2014). Programme informatique. Tire de http:// www.commentcamarche.net;

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 54

-Radhakrisha murthy I. V. et Krishnamacharyulu S. K. G., (1990). A Fortran 77 Program

to fit a polynomial of any order to Potential Field Anomalies. Jour. Assoc. Expl. Geophysics.
Vol X-N° 2. P 99 -105.

-Schoeffler J., (1975). Gravimétrie appliquée aux recherches structurales et à la prospection

pétrolière et minière.Technic. Turbo Pascal 7.0, Guide de reference. Borland International;

- Simpson S.M., (1954).Least-square polynomial fitting to gravitational data and density

plotting by digital computer .Geophysics Vol 19.P 808-811

- Samuel Mbom A., (1997). Investigations géophysiques en bordure du craton du Congo

(région Abong-Mbang/Akonoliga-Cameroun) et implications structurales.Thèse Doct PhD
Univ.Poitiers, P 138 - 155 ;

-Zanga Amougou A., (2015).investigation gravimétrique des structures profondes et

superficielles de la zone Batouri-Berberati (Est-Cameroun/Ouest de la république
centrafricaine). Thèse Doct PhD Univ.Yaoundé I, Cameroun P. 36-48;

- Noutchogwe Tatchum Constant B., (2004). Apport de la gravimétrie à l’étude de la

bordure méridionale du plateau de l’Adamaoua (Cameroun) UYI P. Thèse Doct PhD
Univ.Yaoundé I, P 69- 71 ;

- Zeng H., (1989). Estimation of the degree of polynomial fitted to gravity anomalies and its
applications Geophysic.Prospec.37, P. 959-973.

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 55

ANNEXES

Annexe 1 : Tableau du fichier Excel des données gravimétriques obtenues après

digitalisation de la carte de Bouguer (Zanga A., 2015).

longitude: latitude: Bouguer(mGals)

X(km) Y(km)
-80
1525,54 369,996 -70
1525,54 399,996 -70
1525,54 417,853 -65
1525,54 424,996 -60
1525,54 432,138 -55
1525,54 435,709 -50
1525,54 446,424 -45
1525,54 453,567 -40
1525,54 464,281 -40
1525,54 514,28 -45
1525,54 528,566 -50
1525,54 542,852 -55
1525,54 549,994 -60
1525,54 560,709 -65
1525,54 567,851 -70
1578,3318 253,172 -70
1578,3318 346,425 -65
1578,3318 371,425 -60
1578,3318 389,282 -55
1578,3318 403,567 -50
1578,3318 414,282 -45
1578,3318 428,551 -40
1578,3318 439,281 -35
1578,3318 460,709 -35
1578,3318 482,138 -40
1578,3318 499,995 -45
1578,3318 514,28 -50
1578,3318 521,423 -55
1578,3318 528,566 -60
1578,3318 539,28 -65

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 56

1578,3318 549,994 -70

1578,3318 564,28 -70
1631,095 249,998 -70
1631,095 317,854 -65
1631,095 346,425 -60
1631,095 364,282 -55
1631,095 385,71 -50
1631,095 399,996 -45
1631,095 410,71 -40
1631,095 421,424 -35
1631,095 432,138 -30
1631,095 446,424 -30
1631,095 478,567 -35
1631,095 492,85 -40
1631,095 496,424 -45
1631,095 507,138 -50
1631,095 514,28 -55
1631,095 521,423 -60
1631,095 532,138 -65
1631,095 539,28 -70
1631,095 549,994 -75
1631,095 557,137 -80
1631,095 567,851 -75
1681,094 260,712 -70
1681,094 278,569 -70
1681,094 303,568 -70
1681,094 335,711 -65
1681,094 353,568 -60
1681,094 364,282 -55
1681,094 378,568 -50
1681,094 389,282 -45
1681,094 399,996 -40
1681,094 414,282 -40
1681,094 474,995 -45
1681,094 482,138 -50
1681,094 492,852 -55
1681,094 507,138 -60
1681,094 517,852 -65
1681,094 532,137 -70
1681,094 542,852 -75
1681,094 557,137 -80
1725,528 249,998 -75
1725,528 285,711 -70
1725,528 349,996 -65
1725,528 364,282 -60
1725,528 371,424 -55

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 57

1725,528 378,568 -50

1725,528 392,853 -50
1725,528 407,139 -50
1725,528 414,281 -45
1725,528 435,709 -40
1725,528 453,567 -40
1725,528 471,424 -45
1725,528 492,852 -50
1725,528 492,852 -55
1725,528

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 58

Tableau du fichier Excel des données gravimétriques après filtrage polynomial à l’ordre
3.

longitude: X latitude: Y Bouguer Coef Reg Ris

(km) (km) (mGals)
1525,54 369,996 -80 -17607,92169 -76,63202319 -3,36798
1525,54 399,996 -70 29,50393158 -65,44260683 -4,55739
1525,54 417,853 -70 3,590787739 -59,2118419 -10,7882
1525,54 424,996 -65 -0,015862291 -56,90920233 -8,0908
1525,54 432,138 -60 -0,010550998 -54,74866544 -5,25133
1525,54 435,709 -55 0,015507358 -53,72851437 -1,27149
1525,54 446,424 -50 2,77574E-06 -50,94142668 0,94142
1525,54 453,567 -45 3,64529E-06 -49,33955618 7
1525,54 464,281 -40 -2,56867E-06 -47,3766189 4,33955
1525,54 514,28 -40 -1,01793E-05 -47,29036127 6
1525,54 528,566 -45 -50,71210198 7,37661
1525,54 542,852 -50 -55,99207792 9
1525,54 549,994 -55 -59,38404052 7,29036
1525,54 560,709 -60 -65,47616191 1
1525,54 567,851 -65 -70,23796563 5,71210
1578,3318 253,172 -70 -81,27139588 2
1578,3318 346,425 -70 -70,16170938 5,99207
1578,3318 371,425 -65 -62,54529647 8
1578,3318 389,282 -60 -56,96472538 4,38404
1578,3318 403,567 -55 -52,68333975 1
1578,3318 414,282 -50 -49,69548461 5,47616
1578,3318 428,551 -45 -46,15662553 2
1578,3318 439,281 -40 -43,92407722 5,23796
1578,3318 460,709 -35 -40,89233299 6
1578,3318 482,138 -35 -40,26284912 11,2714
1578,3318 499,995 -40 -42,00699259 0,16170
1578,3318 514,28 -45 -45,13748231 9
1578,3318 521,423 -50 -47,34249606 -2,4547
1578,3318 528,566 -55 -50,00364716 -3,03527

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 59

1578,3318 539,28 -60 -54,89906862 -2,31666

1578,3318 549,994 -65 -60,94589623 -0,30452
1578,3318 564,28 -70 -70,92975602 1,15662
1631,095 249,998 -70 -71,97995517 6
1631,095 317,854 -70 -69,77253581 3,92407
1631,095 346,425 -65 -63,30290953 7
1631,095 364,282 -60 -58,55663513 5,89233
1631,095 385,71 -55 -52,72913995 3
1631,095 399,996 -50 -49,0421381 5,26284
1631,095 410,71 -45 -46,507641 9
1631,095 421,424 -40 -44,25427932 2,00699
1631,095 432,138 -35 -42,35716739 3
1631,095 446,424 -30 -40,51161285 0,13748
1631,095 478,567 -30 -40,08556829 2
1631,095 492,85 -35 -41,96179552 -2,6575
1631,095 496,424 -40 -42,66256191 -4,99635
1631,095 507,138 -45 -45,36257091 -5,10093
1631,095 514,28 -50 -47,68946823 -4,0541
1631,095 521,423 -55 -50,46436946 0,92975
1631,095 532,138 -60 -55,51512812 6
1631,095 539,28 -65 -59,50612534 1,97995
1631,095 549,994 -70 -66,48553369 5
1631,095 557,137 -75 -71,8327853 -0,22746
1631,095 567,851 -80 -80,95002509 -1,69709
1681,094 260,712 -75 -72,00490082 -1,44336
1681,094 278,569 -70 -72,56818261 -2,27086
1681,094 303,568 -70 -70,56093955 -0,95786
1681,094 335,711 -70 -64,5892064 1,50764
1681,094 353,568 -65 -60,3049037 1
1681,094 364,282 -60 -57,59211944 4,25427
1681,094 378,568 -55 -53,96474339 9
1681,094 389,282 -50 -51,33409233 7,35716
1681,094 399,996 -45 -48,86383314 7
1681,094 414,282 -40 -45,94935646 10,5116
1681,094 474,995 -40 -42,37269943 1
1681,094 482,138 -45 -43,29490099 10,0855

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 60

1681,094 492,852 -50 -45,36126592 7

1681,094 507,138 -55 -49,50864484 6,96179
1681,094 517,852 -60 -53,76043282 6
1681,094 532,137 -65 -61,09392596 2,66256
1681,094 542,852 -70 -67,94069815 2
1681,094 557,137 -75 -79,00594584 0,36257
1725,528 249,998 -80 -73,62478579 1
1725,528 285,711 -75 -74,43943079 -2,31053
1725,528 349,996 -70 -62,66839812 -4,53563
1725,528 364,282 -65 -59,04692366 -4,48487
1725,528 371,424 -60 -57,23285933 -5,49387
1725,528 378,568 -55 -55,44558226 -3,51447
1725,528 392,853 -50 -52,04267475 -3,16721
1725,528 407,139 -50 -49,015788 0,95002
1725,528 414,281 -50 -47,69926674 -2,9951
1725,528 435,709 -45 -44,80329728 2,5681
1725,528 453,567 -40 -43,92919831 0,56094

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 61

Annexe 2

% programme principal: séparation of least-square polynomial.m

% séparation polynomiale par moindre carrés

% DIMO HOUPA G. Sonia

% Laboratoire de Physique/ Ecole Normale Supérieure

% Université de Yaounde I

Excel document=input ('saisir le nom du document excel:');

[DATA, TILTE, DOCUMENT] =xlsread (exceldocument);
No=size (DATA, 1);
N=input ('saisir le degre du polynome:');
X=zeros (No, 1);
Y=zeros (No, 1);
B=zeros (No, 1);
head= {'Coef';'Reg';'Ris'}';
%initialization

for i=1:No
X (i) =DATA (i, 1);
Y (i) =DATA (i, 2);
B (i) =DATA (i, 3);
end

M= (N+1)*(N+2)/2;
A=zeros (M, M);
E=zeros (M, 1);
Re=zeros (No, 1);
Ri=zeros (No, 1);

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 62

t1=0;
for j=0:N
for l=j:-1:0

t2=0;
t1=t1+1;
for k=0:N
for s=k:-1:0
for i=1:No
A(t1,t2)=A(t1,t2)+(X(i)^(l+s))*(Y(i)^(k+j-l-s));
end
end
end

for i=1:No
E (t1) =E (t1) +B (i)*((X (i)) ^l)*(Y (i) ^ (j-l));
end
end
end

C=E\A;

for i=1:No
k=0;
for j=0:N
for l=j:-1:0
k=k+1;
Re(i)=Re(i)+C(k)*(X(i)^l)*(Y(i)^(j-l));
end

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016

 63

end
end

for i=1:No
Ri(i)=B(i)-Re(i);
end

xlswrite(exceldocument, head,strcat('D',int2str(1),':F',int2str(1)))
xlswrite(exceldocument, C,strcat('D',int2str(2),':D',int2str(m2+2)))
xlswrite(exceldocument, Re, strcat('E',int2str(2),':E',int2str(N+2)))
xlswrite(exceldocument, Ri,strcat('F',int2str(2),':F',int2str(N+2)))

Mémoire DIPES II * DIMO HOUPA Ginette Sonia * 2015 -2016