Starting from any two compactly supported refinable functions in L-2(R) with dilation factor d,we show that it is always possible to construct 2d wavelet functions with compact support such that they generate a pair of dual d-wavelet frames in L-2(R). Moreover, the number of vanishing moments of each of these wavelet frames is equal to the approximation order of the dual MRA; this is the highest possible. In particular, when we consider symmetric refinable functions, the constructed dual wavelets are also symmetric or antisymmetric. As a consequence, for any compactly supported refinable function phi in L-2(R), it is possible to construct, explicitly and easily, wavelets that are finite linear combinations of translates phi(d . - k), and that generate a wavelet frame with an arbitrarily preassigned number of vanishing moments.We illustrate the general theory by examples of such pairs of dual wavelet frames derived from B-spline functions.