A numerical scheme based on hybrid central finite-volume and finite-difference method is presented to model Green-Naghdi water wave equations. The governing equations are reformulated into the conservative form, and the convective flux is estimated using a Godunov-type finite volume method while the remaining terms are discretized using finite difference method. To enhance the robustness of the model, a central-upwind flux evaluation and a well-balanced non- negative water depth construction are incorporated. Numerical tests demonstrate that present model has the advantages of stability preserving and numerical efficiency.