In this paper the optimum position of optical phase conjugator (OPC) and the optimal dispersion coefficients of filler sections in WDM system with the conventional mid-span spectral inversion (MSSI) are numerically induced and then applied into $16{ imes}40$ Gbps WDM systems with 1,000km non zero - dispersion shifted fiber(NZ-DSF) in order to efficiently compensate for the distorted overall channels. It is confirmed that the compensation extents of overall WDM channels are more improved by applying the induced optimal parameters into WDM system than those in WDM system with the conventional MSSI. So it is expected to alternate with the forming method of the symmetrical distributions of power and local dispersion by applying these optimal parameters into the real WDM system, which generate a serious problem of applying the OPC into multi-channels WDM system if it is not formed. It is also confirmed that two optimal parameters depend on each other, but less related with the finding procedure. And, it will be possible to realize the flexible system design by applying the methods proposed in this paper into the real WDM system with OPC.