The coupled nonlinear Schr*idinger equations for two wavelength optical pulses in a birefringence fiber are given. The model for treating polarization mode dispersion in two channels WDM system is established. Based on this model, optical signal propagation behaviors in two channels WDM with polarization mode dispersion are numerically simulated. The influences of polarization mode dispersion on two channels WDM system are also analyzed.