A controlled rate of propagation of the fluid saturation front or thermal front is desired in may oil and gas and geothermal operations. Natural fractures and fractures created during hydraulic stimulation may have heterogeneous hydraulic properties resulting in uneven flow distributions, therefore leading to short-circuiting and breakthrough issues. The present invention relates to wellbores connected hydraulically by multiple fracture zones; methods are directed to control for even flow distribution among fractures, regardless of heterogeneities in fracture hydraulic properties, and to control propagation of saturation fronts and thermal fronts in subsurface reservoirs.