Abstract: The paper gives the coupling part differential equation system about ZD two-phase flow in discrete media under certain assumption. The permeabilities and fracture apertures in the equations are expanded into Karhunen-Loeve series. Other parameters are expanded into polynomial series of stochastic functions. The errors caused by the series trimming are analyzed.The eigenvalues of random variables can be solved when the basic unknowns are solved.