In this paper, we develop and analyze a mixed finite element method for a nonlinear, higher-order model describing nonlinear wave phenomena and exhibiting important conservation properties. A central goal of our approach is to ensure that these properties are preserved at the discrete level while avoiding the challenges typically encountered when constructing finite element subspaces of $H^2$ as would be required in a standard continuous Galerkin framework. At the continuous level, we establish well-posedness and characterize the solution through energy laws and mass conservation. For the semi-discrete formulation, we derive error estimates in various Bôchner spaces. Furthermore, we establish that the implicit fully discrete scheme is well-posed, converges with optimal order and consistent with both mass conservation and an entropy dissipation law. Finally, we confirm the theoretical findings and conservation properties on a set of benchmark problems.