The aim of this paper is to introduce a new wavelet method for presenting approximate solutions of multitype variable-order (VO) fractional partial differential equations arising from the modeling of phenomena. In specific, this paper focuses on the numerical solution of the VO-fractional mobile-immobile advection-dispersion equation, Klein Gordon equation and Burgers equation. These equations are converted into a system of algebraic equations with the assistance of the bivariate Genocchi wavelet functions, their operational matrices, and the variable-order fractional Caputo derivative operator. Also, we present a new technique to get the operational matrix of integration and VO-fractional derivative. The modified operational matrices for solving the proposed equations are powerful and effective. So that, the accuracy of these matrices directly affects the implementation process. Finally, we consider numerical examples to confirm the superiority of the scheme, and for each example, exhibit the results through graphs and tables.