This study proposes a seismic reliability?based water distribution system (WDS) optimal design model that minimizes total costand maximizes seismic reliability. Here, seismic reliability is defined as the ratio of the available quantity of water to the required demandunder stochastic earthquake events. A new evaluation model is used to assess seismic reliability, while a multiobjective harmony search(MOHS) based on a ranking approach is used for optimization. The Anytown network was modified for the demonstration of the proposedmethod. First, this study performs the sensitivity analysis of MOHS parameter values [i.e., harmony search consideration rate (HMCR) andpitch adjustment rate (PAR)] to identify the best parameter set in a pipe-sizing problem of an Anytown network. Then, Pareto optimalsolutions with three different tank configurations are obtained and compared with respect to the final Pareto fronts and the system designs.For the sensitivity analysis, it reveals that higher PAR and lower HMCR values are also required to maintain high searchability in a multiobjective(MO) framework. In addition, Pareto-optimal solutions found for networks with tanks dominated those found for those withouttanks.