At last, the problem can be solved with the general algebraic modeling system (GAMS) to acquire the optimal design scheme of water utilization networks.