The conventional water network synthesis approach greatly simplifies wastewater treatment units by using fixed recoveries, creating a gap for their applicability to industrial processes. This work describes a unifying approach combining various technologies capable of removing all the major types of contaminants through the use of more realistic models. The following improvements are made over the typical superstructure-based water network models. First, unit-specific shortcut models are developed in place of the fixed contaminant removal model to describe contaminant mass transfer in wastewater treatment units. Shortcut wastewater treatment cost functions are also incorporated into the model. In addition, uncertainty in mass load of contaminants is considered to account for the range of operating conditions. Furthermore, the superstructure is modified to accommodate realistic potential structures. We present a modified Lagrangean-based decomposition algorithm in order to solve the resulting nonconvex mixed-integer nonlinear programming (MINLP) problem efficiently. Several examples are presented to illustrate the effectiveness and limitations of the algorithm for obtaining the global optimal solutions.