The saturated gain performance of the two-pump fiber optical parametric amplifier (2-P FOPA) in the presence of random dispersion fluctuation as well as fourth-order dispersion coefficient are investigated. The gain spectra of the 2-P FOPA is simulated by solving the coupled equations numerically, considering the dispersion fluctuations along the fiber as a stochastic process with a given standard deviation and correlation length. The gain performance of 2-P FOPA is analysed with the variation of correlation length, standard deviation and fourth-order dispersion coefficient, respectively. It shows that, the gain saturation curves of 2-P FOPA exhibit obvious differences in the presence of dispersion fluctuations and fourth-order dispersion coefficients. The peak gain reduces with the increase of the standard deviation and the reduction is higher in case of shorter correlation lengths, also the variations of 2-P FOPA saturation powers are affected by the standard deviation and correlation length significantly. In addition, the performance of the 2-P FOPA varies greatly while considering the fourth-order dispersion coefficient, where the peak gain remains the same, but the 3 dB bandwidth is wider with considering the fourth-order dispersion coefficient. The numerical analysis is meaningful for the research of 2-P FOPA saturated gain performance, which having attracted great research interests in the optical fiber communication systems.