In this paper, the notion of multiple vector-valued multiresolution analysis of space L-2(R-S, C-rxr) is introduced. A method for constructing biorthogonal multiple vector-valued wavelet packets in higher dimensions is presented and their properties is investigated by means of time-frequency analysis method, matrix theory and operator theory. Three biorthogonality formulas concerning these wavelet packets are obtained. Finally, new Riesz bases of space L-2(R-S, C-rxr) is obtained by constructing a series of subspaces of biorthogonal multiple vector-valued wavelet packets. (c) 2007 Elsevier Ltd. All rights reserved.