Stochastic approximation (SA) algorithms are proposed to identify a multi-input and multi-output (MIMO) Wiener system, in which the system input is taken to be a sequence of independent and identically distributed (i.i.d.) Gaussian random vectors u(k) is an element of N(0, I). The algorithm for identifying the nonlinear part is designed with multi-variable kernel functions. Under suitable conditions, we show that the estimates of the coefficients of the linear subsystem and of the values of the nonlinear function converge to the respective true values with probability one.