Les méthodes géophysiques aéroportées sont très utilisées pour la prospection du sous-sol à l'échelle régionale, permettant ainsi de couvrir de grandes surfaces en particulier les zones difficiles d'accès. Le sujet de thèse concerne le développement de techniques d'interprétation des données géophysiques pour le problème des ressources naturelles et de l'environnement. La première partie de cette thèse concerne le développement d'une méthode de calcul direct des composantes des champs gravimétrique et magnétique à partir d'une structure (corps) géologique modélisé(e) par une grille à base de tétraèdres, ce qui permet de représenter des modèles géologiques très complexes particulièrement en présence de zones faillées et hétérogènes avec un nombre d'éléments optimal. Plusieurs techniques d'inversions utilisent des contraintes mathématiques pour la résolution du problème inverse en modélisation. Ces contraintes permettent de réduire le nombre de modèles possibles. Cependant les solutions proposées appelées aussi « le modèle le plus probable » présentent des solutions lisses, ce qui est loin de représenter la réalité géologique. Pour éluder ce problème, les deuxièmes et troisièmes parties de la thèse proposent des améliorations majeures du processus d'inversion par l'utilisation des méthodes géostatistiques telles que la Simulation Gaussienne Séquentielle ou la Co-Simulation dans le cas d'une inversion conjointe afin d'estimer les probabilités a posteriori des modèles simulés. La quatrième partie de ce mémoire présente une alternative à la simulation de plusieurs variables. L'apprentissage du réseau de neurones supervisé par un certain nombre de points permet d'établir une relation entre les différentes variables