Dans cet exposé, nous présentons un nouveau modèle de type Boussinesq, dont le but est de pouvoir propager correctement les vagues et leurs cinématiques sur des domaines étendus allant jusqu'à plusieurs kilomètres au large. La construction du modèle repose sur trois idées: la première est de formuler le problème en fonction d'un opérateur de Dirichlet-Neumann exprimé, non pas à la surface libre de manière classique, mais au niveau de la surface du fluide au repos, de manière à travailler sur un opérateur statique. La seconde idée est de chercher une approximation de cet opérateur au moyen de séries de Taylor tronquées et d'approximants de Padé. La troisième est enfin d'utiliser une décompositionartificielle du fluide en deux couches de même densité, de manière à diminuer l'ordre des dérivées du problème. Le modèle final comprend ainsi quatre équations (en 2DH) ne faisant intervenir que des dérivées secondes au maximum, et nous montrons via une analyse linéaire et des simulations numériques non-linéaires que le modèle permet de propager des vagues avec précision jusqu'en eaux profondes.