The SU(2) lattice gauge system with dynamical Wilson fermions is investigated by the variational cumulant expansion method with the novel criterion to determine the variational parameters. The agreement of calculated order parameter EP of the pure gauge sector with Monte Carlo data is much improved, especially, in the intermediate coupling region. The fermionic contributions are small to k6 hopping parameter expansion, but the general feature is in accordance with the Monte Carlo data qualitatively. The necessity of higher order calculations in the intermediate coupling region is discussed.