The analytical expression for the cross spectral density function of partially coherent circular edge dislocation beams propagating in the deep dermis of mouse tissue is derived based on the generalized Huygens-Fresnel principle， the effects of the initial beam parameters （the beam wavelength λ and the number of circular edge dislocations n） and the propagation distance z on the normalized intensity distribution， phase evolution and propagation trajectory of the beam are investigated. The results show that the central intensity of the partially coherent circular edge dislocation beam with n dislocation number is the largest in the source plane， and 2n secondary peaks symmetrically distribute on both sides. With the increment of propagation distance， the intensity distribution gradually evolves from multi-peak to single-peak， the longer the wavelength is， the smaller the n is， the faster the intensity distribution evolution is. The more the number of dislocations is， the better the beam stability is. Additionally， in the source plane the radius of the innermost ring of n circular edge dislocations decreases as the number of dislocations increases. Owing to the combination of the biological tissue turbulence induction and diffraction effect， the n circular edge dislocation has split into n sets of coherent vortices whose topological charges are "+1" and "-1" from the beginning of the propagation， respectively. As the propagation proceeds， n sets of coherent vortices with topological charges of "+1" and "-1" will be generated. The bigger the wavelength is， the more the n is， the faster the evolution of the beam phase， the more the distribution of the coherent vortices tends to be concentrated， and finally all the coherent vortices are annihilated. The smaller the distance between the coherent vortices， the earlier the annihilation. The longer the wavelength is， the larger the value of n is， the earlier the annihilation of the initial coherent vortices pairs is， and the longer the propagation distance to the total annihilation is. The smaller the wavelength is， the larger the value of n is， the earlier the annihilation of the newly generated coherent vortices pairs is， and the longer the propagation distance to the total annihilation is.