We analyze flowback production data of 502 multi-fractured horizontal oil and gas wells completed in the Montney Formation and 83 oil wells completed in the Duvernay and Horn River Formations. The main goal of this paper is to evaluate the possibility of distinguishing between formation and fracturing water based on the water-flowback response. We hypothesize that: 1) the slope of water-flowback harmonic decline (HD) profile is reversely proportional to formation water mobility, 2) the deviations from the unit slope on rate-normalized pressure (RNP) plots is proportional to the slope of HD, and 3) the slope of water-flowback HD correlates with the initial water saturation (Swi). To verify our hypothesis, we 1) classify the observed HD trends of water-flowback rate based on slopes, 2) construct RNP diagnostic plots of the studied wells, 3) analyze log data and estimate average Swi by using Archie equation (1952) for the studied wells. 4) investigate the effects of Swi on the water-flowback pattern. The results show that there are two distinct flowback patterns. The first flowback pattern shows sharp slope (>10-41/day) of water-flowback HD profile and relatively high slope values (0.64 to 0.984 kpa/m3) of the corresponding RNP plots. However, the second pattern shows very low slope of HD (<5 × 10-5 1/day), with some wells showing no significant decline of water rate through the entire flowback process, also relatively low slope values (0 to 0.23) of the corresponding RNP plots. Analysis of the log data shows a positive correlation between Swi and slope of water-flowback HD profile. We also found that the slopes are proportional to the slope of RNP. These results indicate that as Swi increases, slope of HD decreases and there is more deviation from the unit-slope on the RNP plots.