Hydraulic fracturing enables the commercial development of unconventional resources in shales and tight formations. The conductivity and complexity of created fractures are critically dependent on the rheology of fracking fluid and the mechanics properties of rocks. Literatures show that both the rheology of fracturing fluid and fracture propagation dynamics are affected by the temperature of fracturing fluid. Neglecting the temperature transient behaviour may defeat the purpose of fracturing optimization during fracture initiation, propagation, and sand packing. The objective of this paper is to investigate the impact of temperature on fracturing design by studying the transient temperature behaviour across a complex wellbore using numerical modelling by coupling a finite difference heat transfer model with a dynamic fracture propagation model. The study results show that with the injection of cold fracturing fluid, hydraulic fracture propagation is decelerated, and production prediction is thus lessened compared with the case ignoring temperature effect. For multistage fractured wells, fracture geometry enlarges along the fluid flow direction in a horizontal segment. This potentially lowers the cost of hydraulic fracturing designs.
The exploration of oil, gas, and geothermal reservoirs has been getting deeper and deeper in these decades [
Research on temperature impacts of fracking is rare compared with pressure. Meyer [
In this paper, a numerical solution to wellbore fluid heat transient transfer is proposed based on Jiang’s transient model and used to predict the temperature behaviour during the hydraulic fracturing process. The temperature effects on fluid viscosity and thermal stresses are included into a dynamic fracture propagation model to investigate coupled temperature impacts on fracture growth by analysing the fracture geometry in deep reservoir well fracturing cases.
Heat transfer during a hydraulic fracturing process occurs when cold fracturing fluid interacts with a high-temperature wellbore and formation through convection and conduction. To include more complexity, in this case, we define a well drilled in a deep reservoir into two separate segments: the vertical and horizontal segments. The structures of segments are shown in Figures
Structure of the vertical segment of a complex wellbore.
Structure of the horizontal segment of a wellbore.
For a completed well before fracking, we write the energy balance equation coupled with Fourier’s law in the form of
And the continuity equation is
We apply the spatial and time finite difference to these governing equations, and for different flow regions shown in Figure
The factors
To connect the vertical and horizontal segments in the case of horizontal wells, a boundary condition between horizontal piece inlet temperature
We use a symmetrical initial temperature distribution around the vertical segment of a wellbore along the radial direction, while heat flux of the upper and lower sides is evaluated separately along the horizontal piece. This is calculated explicitly, and fluid temperature is updated at the start of every time step.
During the hydraulic fracturing process, a desirable viscosity is required for proppant transportation along a fracture. Fluid viscosity is also reported to have significant impact on fracture propagation [
Stress state changes when near-wellbore field temperature distribution varies, which significantly affects hydraulic fracture propagation. We use a revised model from Gogoi [
We use a fracture geometry model to couple the temperature-influenced parameter for analysing the impact on fracture geometry and production forecast. A few assumptions are made before setting up the model:
Constant fluid leak-off coefficient Hydraulic fractures fully penetrate the formation in the vertical direction The formation is under a normal stress regime. Fracture propagates in the normal to minimum horizontal stress direction
The models of fracture width and length are shown as follows:
To account for the change in the stress state, more implementations are made to further calibrate the fracture length and width.
The fracture half-length is a function of pore pressure difference
The correlation for fracture width is derived and shown in Appendix
A workflow of step-by-step modelling and implementation of a hydraulic fracture is shown as follows:
Model the temperature behaviour at initial reservoir condition for vertical and horizontal segments Calculate fracking fluid viscosity change based on the temperature distribution at the first pumping time step Calculate the pore pressure and normal stress state of formation based on data from step 2 Calculate the thermal stress distribution and superpose with normal stress of formation using revised Gogoi’s method for this pumping time step Model the fracture half-length and width propagation at this pumping time step Repeat steps 2-5 for every pumping time step until the end of the pumping schedule
A sample case from field data is used to validate the model. A vertical well is drilled at a depth of 7530 meters and open-hole completed and fractured at a perforation depth of 7430 meters using brine. The averaged corresponding properties are shown in Tables 1 and 2 in Appendix
The transient temperature profile of formation in the radial direction at the bottomhole is simulated using the heat transfer model, and results are shown in Figure
Temperature profile with time at the formation.
Brine viscosity is also affected by this process. While the wellbore and surrounding formation are being cooled, the heat is transferred into fracking fluids, heating them up. We calibrated the fracking fluid viscosity using the transient heat profile and noticed a decrease from 1 to 0.067 cP during the process while the fracking fluid is being heated.
Thermal-induced stress during pumping is calculated and superposed into the original stress state profiles, shown in Figure
Thermal-induced stress during fracturing.
Coupling all the temperature impacts described in the above sections, we modelled the dynamic growth of a hydraulic fracture with respect to transient heat behaviour. A comparison between the isothermal fracture half-length and the transient heat model fracture half-length is shown in Figure
Comparison of fracture half-length between the isothermal fracture model and the transient heat fracture propagation model.
On the other hand, the impact on fracture width is much smaller. We only noticed a change of about 0.016 mm in fracture width, which is negligible compared with the total fracture width. This may result from the shape of hydraulic fractures, where the fracture width itself is already small enough.
Simulation has been made by using the different fracture geometries with and without transient heat transfer. We compare the gas production with true field data, and the results are shown in Figure
Comparison of production predictions between the isothermal fracture model and the transient heat fracture propagation model. Dots represent history measurements.
In this specific case, we can easily infer that the isothermal fracture geometry failed to predict the production at a good accuracy in deep reservoirs. It tends to overestimate the fracture length and thus overestimate the total production. On the other hand, the transient heat model gives a better estimation of production trends, inferring that the method we use is validated.
Horizontal wells tend to have a much longer horizontal segment compared with the vertical piece. When we have a long horizontal well, fracking fluids have longer travel time inside the wellbore and cool the formation in a slightly different manner. In this case, the formation temperatures at the inlet of the horizontal segment and the end of the horizontal wellbore are different, causing difference in multistage hydraulic fracturing cases.
We added a horizontal wellbore of 2000 meters to the vertical case and apply a 5-stage hydraulic fracturing process. Stage 1 is at the ending of the horizontal segment while stage 5 is at the inlet position. The simulated temperature profile for the formation near the wellbore along the horizontal segment is shown in Figure
Formation temperature profile at the wellbore along horizontal segments at different pumping times.
We notice that the formation is being cooled down during pumping, and the end of the horizontal wellbore tends to have a much higher temperature than the inlet. Once we take this into consideration, we can model the multistage averaged fracture half-length at different locations, shown in Figure
Multistage hydraulic fracturing with respect to transient temperature behaviour.
We notice that in this case, the hydraulic fractures at different stages (and thus different distances) have shown their relationship with temperature. The 1st-stage fractures at the ending of the horizontal wellbore have the longest half-length, and the 5th-stage fractures at the inlet of the horizontal segment are the shortest. This trend cannot be found if we use isothermal fracture conditions, for which this difference around 15 m will also have impact on production predictions.
Based on the previous analysis and discussions, the following conclusions are drawn:
Temperature has shown its impact on different mechanisms during the hydraulic fracking process, including changing the fracking fluid viscosity and formation stress distributions. These changes affect the half-length of hydraulic fractures and therefore have an impact on production predictions. Current fracture models with isothermal assumption may be valid in conventional reservoirs with low reservoir temperature and small geothermal gradients but become inaccurate in deep reservoirs as the temperature impacts become nonnegligible By coupling the mechanisms into numerical simulations, considering these temperature impacts tends to reduce simulated fracture half-length compared with isothermal fracture conditions in deep reservoirs and thus reduce production prediction. This temperature impact is significant and cannot be ignored. On the other hand, the change of fracture width is much smaller and still negligible For multistage hydraulic fractures, the temperature behaviour along the horizontal wellbore tends to reduce the hydraulic fracture half-length with respect to different stage locations. As multistage fracturing usually starts from the ending of horizontal segments and moves along to the inlet, hydraulic fracture half-length also tends to decrease along this process. As the horizontal wellbore is usually much longer than the vertical piece, this impact on fracture length cannot be neglected
The
For an open-hole completed production well which is meant to be fractured, there are 3 different regions to be analysed: pipeflow, casing and cementing, and formation:
For the pipeflow region (subscript 1),
For the casing and cementing region (subscript 2),
For the formation region (subscript 3),
The derivation of fracture width calibration with respect to stress-state variation is shown below.
For a propagating fracture, the fracture width can be described by a separation [
The net stress effect is approximated as being purely elastic. Detournay et al. [
The net pressure poroelastic effect can be described by Boone and Detournay [
Boone and Detournay [
Table
Reservoir properties and fracking schedule for the case study.
Property names and units | Values |
---|---|
Well depth (m) | 7530 |
Liquid rate (m3/s) | 0.09 |
Perforation depth (m) | 7430 |
Formation thickness (m) | 50 |
Casing OD (m) | 0.178 |
Liquid injection temperature (°C) | 25 |
Surface temperature (°C) | 25 |
Geothermal gradient (°C | 0.018 |
Formation porosity | 0.076 |
Poisson ratio | 0.28 |
Overburden density (kg/m3) | 1.797 |
Young’s modulus (GPa) | 31.9 |
Geomechanics heterogeneity | 2 |
Leak-off coefficient (m/min0.5) | 0.004 |
Poroelastic coefficient | 0.2 |
Rock permeability (nD) | 400 |
Liquid Knudsen salinity (ppt) | 35 |
Thermal properties for the case study.
Heat conductivity (W/m°C) | Specific heat capacity (J/kg°C) | Density (kg/m3) | Viscosity (Pa·s) | |
---|---|---|---|---|
Formation | 2.20 | 920.0 | 2640 | Not applied |
Casing & cement | 43.33 | 418.7 | 8048.0 | Not applied |
Fluid | 0.586 | 4002.0 | 1000.0 | 0.0011 |
The data used to support the findings of this study are included in the article.
The authors declare that there is no conflict of interest regarding the publication of this paper.