We present a simple numerical model that is used in conjunction with a systematic algorithm for parameter optimization to understand the three-dimensional stochastic intensity dynamics of stimulated Brillouin scattering in a two-mode optical fiber. The primary factors driving the complex dynamics appear to be thermal density fluctuations, transverse pump fluctuations, and asymmetric transverse mode fractions over the beam cross-section.