REPUBLIQUE DU CAMEROUN REPUBLIC OF CAMEROON
PAIX-TRAVAIL-PATRIE PEACE-WORK-FATHERLAND
******** ******** ******** ********
1UNIVERSITE DE DSCHANG UNIVERSITY OF DSCHANG
******** ******** ******** ********
Ecole Doctorale Post Graduate School
******** ********
DSCHANG SCHOOL OF SCIENCES AND TECHNOLOGY
DEPARTEMENT DE PHYSIQUE/ DEPARTMENT OF PHYSICS
Unité de Recherche de Mécanique et de Modélisations des Systèmes Physiques
(UR-2MSP)
Modélisation de la source de l’anomalie Gravimétrique de l’axe Poli-Boki
(Nord Cameroun)
Projet de recherche
Master II en Physique
Option : Mécanique - Energétique
Présenté Par :
MOTEDE DARI DAMAL
Licencié en Physique
Matricule : CM-UDS-18SCI2274
Sous la direction de :
KOUMETIO Fidèle
Maitre de Conférence
Université de Dschang
Année académique : 2020-2021
INTRODUCTION
OBJECTIFS
MATERIELS
METHODES
TRAVAIL A FAIRE
RESULTATS ESCOMPTES
COMPETANCES REQUISES ET ACQUISES
CHRONOGRAMME 1
16/02/2026 2
La géophysique est la science qui utilise les principes de la
Physique à l'étude de la Terre.
Gravimétrie est une science qui mesure, analyse et étudie les
variations du champ gravitationnel terrestre dues au contraste de densité
entre les roches du sous-sol en vue de déterminer la forme et la
constitution de la terre.
Notons qu’une bonne interprétation des données permet de
déterminer les zones de discontinuités du sous-sol. Ceci décrivant les
facteurs géologiques et la tectonique de la zone d’étude. 2
16/02/2026 3
Les études gravimétriques faites par Poudjom-Djomani (1993), Poudjom-
Djomani et Diament (1997) et Noutchogwe (2004) soulignent un soulèvement du
manteau supérieur-Garoua, résultant de la croûte continentale anormalement
mince. L’amincissement de la croûte a été interprété comme l’effet de l’extension
lithosphérique.
Les études de Mouzong et al., (2014) révèlent une estimation de l’épaisseur
des sédiments dans le bassin de Garoua dans l’intervalle 4,4-8,9 km. La carte
d’anomalie résiduelle et l’analyse spectrale montre que l’épaisseur des sédiments
dans le bassin décroit lorsqu’on va de l’Est vers l’Ouest.
Eyike et Ebbing (2015) utilisent un modèle gravimétrique régional 3D afin
de décrire la structure lithosphérique. D’après ces auteurs, la croûte se serait 3
amincie due à une remontée de l’asthénosphère.
16/02/2026
Notre motivation est de contribuer à la connaissance du sous-
sol, à indiquer des zones potentiellement riches en minerais, afin
d’orienter les futures campagnes de recherche minière.²
4
16/02/2026 5
Les travaux antérieurs réalisé dans la région du nord Cameroun
montrent que la zone couvrant Poli-Boki apparait insuffisamment
explorée. Nous pouvons donc nous poser la question, comment
obtenir un model permettant de caractériser la source de l’anomalie
gravimétrique de l’axe Poli-Boki ?
5
16/02/2026 6
Les objectifs principaux de notre travail sont de caractériser la
source de l’anomalie de l’axe Poli-Boki, déterminer le contraste de
densité, la profondeur et la géométrie de la structure perturbatrice.
Pour atteindre ces objectifs principaux on s’est fixé quelques
objectifs spécifiques :
Faire une étude sur les récents résultats ;
Déterminer le(s) types d’anomalies ;
Maîtriser la théorie de la méthode gravimétrique ;
D’étudier les différentes méthodes d’interprétation gravimétriques.
6
16/02/2026 7
Logiciel surfer12.0 ;
7
16/02/2026 8
La dérive correspond à une variation de la
mesure au cours du temps indépendamment des
variations éventuelles de la pesanteur.
𝑽𝟐 −𝑽𝟏
𝑻𝑫 =
𝑻𝟐 −𝑻𝟏
(1)
Si TD est positif, les mesures ont été
surestimées et la correction est négative ;
Si TD est négatif, les mesures ont été sous-
estimées et la correction est positive. 8
16/02/2026 9
En supposant qu’entre la surface de
référence et le point de mesure P, il n’existe pas
de masse.
𝑮𝑴 𝑮𝑴 𝟏 𝟐 𝒉 −𝟐
𝒈= = = 𝒈𝟎 𝟏 + (2)
𝒉+𝑹 𝟐 𝑹𝟐 𝟏+𝒉 𝑹
𝑹
𝑮𝑴
Avec 𝒈𝟎 =
𝑹𝟐
𝟐𝒉
𝚫𝒈𝑭 = 𝒈 − 𝒈𝟎 = − 𝒈𝟎
𝑹
(3)
9
16/02/2026 10
Entre le point de mesure et la surface
de existe une masse.
𝑷𝑨 𝑷𝑨
𝒅𝒈𝟏 = 𝑮𝒅𝒎 = 𝑮𝒅𝒎 (4)
𝑨𝑷𝟑 𝒓𝟑
𝑷𝑨 𝒄𝒐𝒔 𝜽
𝒅𝒈 = 𝒅𝒈𝟏 . 𝒖 = 𝑮𝒅𝒎 .𝒖 = 𝑮𝒅𝒎 (5)
𝒓𝟑 𝒓𝟐
𝒅𝒔 𝒄𝒐𝒔 𝜽 𝒅𝒔 𝒄𝒐𝒔 𝜽
𝒅𝒈 = 𝑮𝝆𝒅𝒛 = 𝑮𝝆𝒅𝒛𝒅𝜴 𝒐𝒖 𝒅𝜴 = (6)
𝒓𝟐 𝒓𝟐
𝜟𝒈 = 𝑮𝝆𝒉𝜴 (7)
𝜟𝒈𝑷 = 𝒈 𝑷 − 𝒈 𝑶 = 𝟐𝝅𝑮𝝆𝒉 ≈ 𝟎, 𝟎𝟒𝟏𝟗𝟏𝝆𝒉 (8) 10
16/02/2026 11
𝒈 = 𝒈𝟎 𝟏 + 𝒌𝟏 𝐬𝐢𝐧𝟐 𝝋 − 𝒌𝟐 𝐬𝐢𝐧𝟐 𝟐𝝋 (9)
𝒈 = 𝟗𝟕𝟖𝟎𝟑𝟐. 𝟕 𝟏 + 𝟎. 𝟎𝟎𝟓𝟑𝟎𝟐𝟒 𝐬𝐢𝐧𝟐 𝝋 − 𝟎. 𝟎𝟎𝟎𝟎𝟎𝟓𝟖𝐬𝐢𝐧𝟐 𝟐𝝋 (10)
𝒅𝒈 𝟏 𝒅𝒈
= (11)
𝒅𝒔 𝑹𝒆 𝒅𝝓
𝒈 ≈ 𝒈𝒆 𝟏 + 𝒌𝟏 𝐬𝐢𝐧𝟐 ∅ (12)
𝚫𝒈𝑳
≈ 𝟎, 𝟖𝟏 𝐬𝐢𝐧 𝟐∅ 𝒎𝑮𝒂𝒍/𝒌𝒎 (13)
𝚫𝒔
11
16/02/2026 12
Cette correction tient compte des
irrégularités de la surface du sol à proximité
de la station de mesure.
Pour faire la correction de Bouguer,
on enlève l'attraction d'une tranche de terrain
d'épaisseur h qu’on peut qualifier de tranche
de Bouguer.
𝒛𝟎
𝚫𝒈𝒕 = 𝑽
𝑮𝝆 𝟐 𝟐 𝟐 𝒅𝒗 (14)
𝒙𝟎 +𝒚𝟎 +𝒛𝟎
12
16/02/2026 13
𝒈𝒕𝒉 = 𝒈𝟎 + 𝚫𝒈𝑭 + 𝚫𝒈𝑷 + 𝚫𝒈𝒕 (15)
𝒂𝑩 = 𝒈𝒐𝒃𝒔 − 𝒈𝒕𝒉 = 𝒈𝒐𝒃𝒔 − 𝒈𝟎 − 𝚫𝒈𝑭 + 𝚫𝒈𝑷 + 𝚫𝒈𝒕 (16)
13
16/02/2026 14
𝑩 = 𝑹é𝒈 + 𝑹é𝒔 (17)
𝑵 𝒋
𝑭 𝑿𝒊 ; 𝒀𝒊 = 𝑩𝒊 + 𝒋=𝟏 𝑳=𝟎 𝑩𝒎 𝑨𝒎 𝑿𝒊 + 𝒀𝒊 (18)
𝒋−𝑳 𝒋 𝒋+𝟑
Avec 𝑨𝒎 = 𝑿𝑳𝒊 𝒀𝒊 ; 𝒎= −𝑳+𝟏
𝟐
𝜺 = 𝒈 𝑿𝒊 ; 𝒀𝒊 − 𝑭 𝑿𝒊 ; 𝒀𝒊 (19)
𝑵𝟎 𝟐
𝑬= 𝒊=𝟏 𝜺𝒊 (20)
𝝏𝑬
= 𝟎 avec 𝟏 ≤ 𝑲 ≤ 𝑴 (21)
𝝏𝑩𝑲
𝑹é𝒈(𝒙𝒊 ; 𝒚𝒊 ) = 𝑭(𝒙𝒊 𝒚𝒊 ) (22)
𝑹é𝒔 𝑿𝒊 ; 𝒀𝒊 = 𝒈 𝑿𝒊 ; 𝒀𝒊 − 𝑹é𝒈(𝑿𝒊 ; 𝒀𝒊 ) (23) 15
15
Cas d’un profil
𝒂𝟎 ∞
𝒇 𝒙 = + 𝒏 𝒂𝒏 𝐜𝐨𝐬 𝒎𝒘𝒙 + 𝒃𝒏 𝐬𝐢𝐧 𝒎𝒘𝒙 (24)
𝟐
𝒍 𝑳
𝟐 𝟐 𝒊𝒘𝒙 𝟐 𝟐
𝒇 𝒙 = 𝒍 𝑭(𝒘)𝒆 ) 𝒅𝒘 𝑭 𝒘 = 𝒇 𝒙 𝒆−𝒊𝒘𝒙 𝒅𝒙 (25)
𝑻 − 𝑻 −𝑳
𝟐 𝟐
𝟐𝝅𝒌𝒏
𝟏 𝑵 −
𝑭𝒌 𝒘 = 𝒏=𝟏 𝒇𝒏 𝒆 𝑵 (26)
𝑵
𝑬𝒌 𝒘 = 𝑭𝒌 (𝒘) 𝟐 (27)
𝟐𝝅
𝑬𝒌 = 𝒆−𝟒𝝅𝒉𝒌 d’où 𝒌 = (28)
𝒉
𝚫(𝒍𝒏𝑬)
𝒉= (29)
𝟒𝝅𝚫𝑲
16
16/02/2026 16
Cas d’une grille
𝑵𝒙 −𝒙 𝑵𝒚 −𝟏 𝟐𝝅 𝟐𝝅
𝑮 𝒏𝒙, 𝒏𝒚 = 𝒊=𝟎 𝒋=𝟎 𝑷𝒏𝒙,𝒏𝒚 𝐜𝐨𝐬 𝑵𝒙 𝑿 + 𝑵𝒚 𝒀 + 𝑸𝒏𝒙,𝒏𝒚 𝐬𝐢𝐧 𝑵𝒙 𝑿 + 𝑵𝒚 𝒀 (30)
𝑵𝒙 𝑵𝒚
𝟐 𝟐
𝑲𝒏 = 𝑷𝒏𝒙,𝒏𝒚 + 𝑸𝒏𝒙,𝒏𝒚 (31)
𝒏= 𝒏𝟐𝒙 + 𝒏𝟐𝒚 (32)
17
16/02/2026 17
Le concept de corps idéal de Parker(1974-1975)
𝒈= 𝑮𝒊 𝒓 𝝆 𝒓 𝒅𝒗 𝒊 = 𝟏, 𝟐, 𝟑 … 𝑵 (33)
𝑵 𝑵
𝒈= 𝒋=𝟏 𝒌=𝟏 𝑮𝒊𝒋𝒌 𝝆𝒋𝒌 (34)
𝑵
𝝆𝟎 𝒑𝒐𝒖𝒓 𝒊=𝟎 𝜶𝒊 𝑮𝒊 > 𝟎
𝝆 𝒓 = 𝑵 (35)
𝟎 𝒑𝒐𝒖𝒓 𝒊=𝟎 𝜶𝒊 𝑮𝒊 ≪ 𝟎
18
16/02/2026 18
Les méthodes directes sont utilisées lorsqu'on a une
connaissance de la géologie de la région d'étude et donc de la
structure qui crée l'anomalie observée.
Modélisation 2D;
Modélisation 2,5D;
Modélisation 3D.
19
16/02/2026 19
Notre travail sera articulé comme suit :
Maîtriser la géologie et la tectonique de la zone couvrant Poli, Boki et leurs environs ainsi que
les résultats des études géophysiques déjà effectuées dans cette région ;
Maîtriser la théorie de la méthode gravimétrique ;
Effectuer la séparation résiduelle-régionale des anomalies gravimétriques de la région d’étude ;
Représenter la carte d’anomalies de Bouguer et les cartes d’anomalies résiduelles ;
Analyser ces cartes et choisir la carte résiduelle appropriée pour l’interprétation quantitative de
l’anomalie de l’axe Poli-Boki ;
Faire l’analyse spectrale en vue de déterminer les profondeurs d’interfaces (méthode
indirecte) ;
Déterminer le contraste de densité du corps intrusif par la méthode inverse en tenant compte
des résultats de l’analyse spectrale (recherche de la solution du corps idéal) ;
Faire la modélisation en tenant compte des résultats des méthodes inverse et indirecte. 20
16/02/2026 20
Nous pensons être à mesure de déterminer les différents
paramètres (profondeur, géométrie du corps, contraste de
densité) nous permettant d’obtenir un modèle qui permettra de
caractériser la source de l’anomalie dans cette zone d’étude. Les
analyses intégrant les valeurs des profondeurs et des densités
permettront de préciser les implications dans l’orientation des
futures recherches minières dans cette zone.
21
16/02/2026
Pour mener à bien ce travail nous devons avoir une
connaissance suffisante sur :
La théorie de la méthode gravimétrique ;
Les méthodes de prospection gravimétrique (interprétation
par séparation régionale-résiduelle) ;
Les logiciels surfer12.0.
22
16/02/2026
A la fin de ce travail nous serons capables d’interpréter et de
modéliser la source d’anomalie à partir des données collectées
dans une zone d’étude.
23
16/02/2026
CHRONOGRAMME DU TRAVAIL
Mois
Novembre Décembre Janvier Février Mars Avril Mai juin
Activités
Synthèse
bibliographique
Théorie de la
méthode
gravimétrique
Revue de la
littérature
Résultat et
discussion
Introduction
générale et
conclusion
générale plus
perspectives
Rédaction
Dépôts et
soutenance
24
16/02/2026
MERCI POUR VOTRE AIMABLE ATTENTION