Flowback rate-transient analysis with spontaneous imbibition effects

Abstract

Analysis of flowback data, gathered immediately after fracture stimulation, can be performed to understand the fluid flow physics, investigate flow regimes, and obtain early estimates of fracture properties. During a hydraulic fracturing treatment, significant amounts of fracturing fluid will leak-off from the fractures into the reservoir due to Darcy flow, capillary, osmotic and electrostatic forces. Capillary invasion of fluids into the reservoir can cause a loss in gas relative permeability, leading to an altered zone near the fracture-matrix interface, therefore impeding the flow of hydrocarbons into the fracture. Due to this phenomenon and other fluid transport mechanisms, a simple application of Darcy's law might not be adequate for describing the fluid flow physics when solid-liquid interaction is significant. To overcome some of the above limitations, spontaneous imbibition effects are modeled at the fracture/matrix interface during the flowback period in this study.This paper presents a semi-analytical model for analyzing two-phase water and gas flowback data, when spontaneous imbibition occurs. This model was developed by solving the fracture and reservoir matrix flow equations simultaneously. The effects of fracture and reservoir matrix pressure gradients on gas and water influx at the fracture-matrix interface are accounted for in order to evaluate the reservoir matrix hydrocarbon influx. The proposed model accounted for spontaneous imbibition driven by capillary forces by quantifying the fluid influx due to capillary processes and adding it to the mass flow equations. Further, capillary pressure effects were incorporated into the PVT properties of matrix pseudovariables. The average phase pressures in the fracture and matrix were calculated iteratively using a modified material balance approach.The proposed semi-analytical model was successfully verified using fully-numerical simulation data. Practical application of the proposed model was then demonstrated using production data from a multi-fractured horizontal well.