Fluid migration in a fracture network plays an important role in the oil accumulation mechanism and hence is key to oil exploration. In this study, we build a model by combining one-dimensional (1D) Navier–Stokes equations, linear elastic equations, and energy equations, and validate the model by reproducing the thickness profile of a fluid-driven crack measured in an experiment. We employ this model to simulate the upward flow of viscous fluid in a single fracture during hydrocarbon migration. The simulation suggests that the parameters of both the fluid and the surrounding rock matrix, as well as the boundary condition imposed on the fracture outlet, affect the upward flow in the fracture. We then extend our model from the single fracture to the bifurcated fracture and the fracture network by maintaining homogeneous pressure and mass conservation at the connection of the channels. We find that the increase in network complexity leads to an increase in the inlet pressure gradient and inlet speed, and a decrease in the outlet pressure gradient and outlet speed. The effective area where the fluid is driven upward from the inlet to the outlet is offset toward the inlet. More importantly, the main novelty of our model is that it allows us to evaluate the effect of inconsistencies in individual branch parameters, such as matrix stiffness, permeability, temperature, and boundary conditions, on the overall upward flow of viscous fluid. Our results suggest that the heterogeneity enforces the greater impact on the closer branches.