Providing water with adequate quality to users is one of the main concerns for water utilities. In most countries, this is ensured through the introduction of disinfectants, such as chlorine, which are subjected to decay over time, with consequent loss of disinfection action and the possible formation of harmful by-products. In this context, water quality models can be a useful tool to support management and, thus, ensure sufficient standards in all network points, but most of these models require the input of reaction parameters which could be difficult to obtain based on the information available to water utilities, especially in the case of complex water distribution networks (WDNs) supplied by more than one source. This study proposes a pragmatic, interval-number-based method to model chlorine decay in complex WDNs by relying on the use of the network hydraulic model and the results of trace analysis, which are exploited to obtain overall reaction rates. The method is applied to the case of a real WDN supplied by water sources with different qualitative features. The results obtained highlight that the method can help water utilities in the identification of overall water quality parameters.