Nous élaborons et analysons mathématiquement des approximations numériques par des méthodes de type volumes finis de solutions faibles de systèmes hyperboliques pour des écoulements géophysiques. Dans une première partie nous approchons les solutions du système de la magnétohydrodynamique en faible épaisseur avec un fond plat. Nous développons un schéma de type Godunov utilisant un solveur de Riemann approché défini via une méthode de relaxation. Des expressions explicites sont établies pour les vitesses de relaxation, qui permettent d'obtenir un schéma satisfaisant un ensemble de bonnes propriétés de consistance et de stabilité. Il conserve la masse, préserve la positivité de la hauteur de fluide, vérifie une inégalité d'entropie discrète, résout les discontinuités de contact même résonantes, donne des vitesses de propagations contrôlées par les données initiales. Des tests numériques sont effectués, validant les résultats théoriques énoncés. Dans une seconde partie nous approchons les solutions du système de la magnétohydrodynamique en faible épaisseur avec fond variable. Nous développons un schéma équilibre pour certains états stationnaires au repos. Nous utilisons la méthode de reconstruction hydrostatique, avec des états reconstruits pour la hauteur d'eau et les composantes du champ magnétique. Nous trouvons des termes correctifs pour les flux numériques par rapport au cadre habituel, et nous prouvons que le schéma obtenu préserve la positivité de la hauteur d'eau, vérifie une inégalité d'entropie semi-discrète et est consistant. Des tests numériques sont effectués, validant les résultats théoriques. Dans une troisième partie nous établissons la convergence d'un schéma cinétique avec reconstruction hydrostatique pour le système de Saint-Venant avec topographie. De nouvelles estimations sur le gradient des solutions approchées sont obtenues par l'analyse de la dissipation d'énergie. La convergence est obtenue par la méthode de compacité par compensation, sous des hypothèses sur les données initiales et la régularité du fond / We build and analyze mathematically numerical approximations by finite volume methods of weak solutions to hyperbolic systems for geophysical flows. In a first part we approximate the solutions of the shallow water magneto hydrodynamics system with flat bottom. We develop a Godunov scheme using an approximate Riemann solver defined via a relaxation method. Explicit formulas are established for the relaxation speeds, that lead to a scheme satisfying good properties of consistency and stability. It preserves mass, positivity of the fluid height, satisfies a discrete entropy inequality, resolves contact discontinuities, and involves propagation speeds controlled by the initial data. Several numerical tests are performed, endorsing the theoretical results. In a second part we approximate the solutions of the shallow water magneto hydrodynamics system with non-flat bottom. We develop a well-balanced scheme for several steady states at rest. We use the hydrostatic reconstruction method, with reconstructed states for the fluid height and the magnetic field. We get some new corrective terms for the numerical fluxes with respect to the classical framework, and we prove that the obtained scheme preserves the positivity of height, satisfies a semi-discrete entropy inequality, and is consistent. Several numerical tests are presented, endorsing the theoretical results. In a third part we prove the convergence of a kinetic scheme with hydrostatic reconstruction for the Saint-Venant system with topography. Some new estimates on the gradient of approximate solutions are established, by the analysis of energy dissipation. The convergence is obtained by the compensated compactness method, under some hypotheses concerning the initial data and the regularity of the topography