Second order statistical moments of the phase fluctuations are obtained taking into account the boundary condition, diffraction effects and polarization coefficients of the ordinary and extraordinary waves. The variance and the correlation function are calculated for arbitrary 3D spectral function of electron density fluctuations containing both anisotropic Gaussian and power-law spectra; anisotropy coefficient and the orientation angle of elongated plasma irregularities. The phase scintillation index and the scintillation level are analyzed numerically. Maximum of the scintillation index for small-scale irregularities is in the interval 0.2-0.3 corresponding to the moderate scintillation intensity, within the weak-scatter regime. Splashes are revealed for different anisotropy factor of elongated largescale irregularities varying orientation angle with respect to the lines of force of geomagnetic field. Scintillation index is calculated for small-scale irregularities using the "frozenin" assumption and taking into account movement of rigid irregularities. Log-log plots of the power spectrum of the intensity fluctuations have the same minimums satisfying the "standard relationship" of scattered ordinary and extraordinary waves. It was shown that the normalized scintillation level growth in both non-fully-developed diffraction pattern and in transition zone increasing anisotropy factor. Rising orientation angle scintillation level decreases and splashes arises in fully developed scintillation region.