The beam spreading, spectral degree of coherence and evolution behavior of the intensity profile of partially coherent Lorentz-Gaussian elliptical multi-Gaussian correlation Schell-Model beam propagating in anisotropic atmospheric turbulence are studied. The analytical expressions of cross-spectral density function, as well as root mean square (rms) beam width, are derived based on the extended Huygens-Fresnel principle and the relationship of Lorentz distribution and Hermitian Gaussian function. With the increases of propagation distance, the elliptical beam will first evolve into Gaussian beam and then change back to elliptical beam. In anisotropic atmospheric turbulence, the influence of the inner scale of turbulence on the spectral degree of coherence and rms beam width are obviously greater than that of the outer scale. For Lorentz-Gaussian elliptical multi-Gaussian correlated Schell-Model beams, better propagation performance was found in anisotropic atmospheric turbulence with larger anisotropic factor and smaller inner scale.