Extending the stochastic mean-field model by including pairing, an approach is proposed for describing evolutions of complex many-body systems in terms of an ensemble of Time-Dependent Hartree-Fock Bogoliubov trajectories which is determined by incorporating fluctuations in the initial state. Non-linear evolution of the initial fluctuations provides an approximate description of quantal correlations and fluctuations of collective observables. Since the initial-state fluctuations break the particle-number symmetry, the dynamical description in which pairing correlations play a crucial role is greatly improved as compare to the mean-field evolution. The approach is illustrated for a system of particles governed by a pairing Hamiltonian.