A Study of Proppant Transport With Fluid Flow in a Hydraulic Fracture

2018 ◽  
Vol 33 (04) ◽  
pp. 307-323 ◽  
Author(s):  
Christopher A. J. Blyton ◽  
Deepen P. Gala ◽  
Mukul M. Sharma
2015 ◽  
Vol 3 (3) ◽  
pp. ST43-ST53 ◽  
Author(s):  
Mehdi Mokhtari ◽  
Azra N. Tutuncu ◽  
Gregory N. Boitnott

Contrary to the assumption in cubic law, the surface of fractures has some degree of roughness, which impacts their fluid dynamics. Incorporating the effect of roughness can improve the simulation of fluid flow in fractures and faults, as well as proppant transport in hydraulic fracturing. To investigate the effect of roughness on the fluid flow, we created a fracture using the Brazilian test, and its roughness was measured using a laser profilometer. Experimental permeability measurements showed a reduction in permeability as the effective stress increased. However, the unmatching surfaces of the fracture prevented its complete mechanical closure. Numerical simulations of the fluid dynamics were conducted on the measured fracture geometry. We determined that the hydraulic fracture aperture is less than the mechanical fracture aperture and that there was anisotropy in the fracture permeability. The ratio of hydraulic fracture aperture to mechanical fracture aperture, as well as anisotropy in fracture permeability, increased when the fracture aperture decreased. The anisotropy in fracture permeability was 45% at the lowest simulated fracture aperture. Integrating the experimental and numerical data, we estimated the fracture porosity and fracture permeability.


2021 ◽  
Author(s):  
Seyhan Emre Gorucu ◽  
Vijay Shrivastava ◽  
Long X. Nghiem

Abstract An existing equation-of-state compositional simulator is extended to include proppant transport. The simulator determines the final location of the proppant after fracture closure, which allows the computation of the permeability along the hydraulic fracture. The simulation then continues until the end of the production. During hydraulic fracturing, proppant is injected in the reservoir along with water and additives like polymers. Hydraulic fracture gets created due to change in stress caused by the high injection pressure. Once the fracture opens, the bulk slurry moves along the hydraulic fracture. Proppant moves at a different speed than the bulk slurry and sinks down by gravity. While the proppant flows along the fracture, some of the slurry leaks off into the matrix. As the fracture closes after injection stops, the proppant becomes immobile. The immobilized proppant prevents the fracture from closing and thus keeps the permeability of the fracture high. All the above phenomena are modelled effectively in this new implementation. Coupled geomechanics simulation is used to model opening and closure of the fracture following geomechanics criteria. Proppant retardation, gravitational settling and fluid leak-off are modeled with the appropriate equations. The propped fracture permeability is a function of the concentration of immobilized proppant. The developed proppant simulation feature is computationally stable and efficient. The time step size during the settling adapts to the settling velocity of the proppants. It is found that the final location of the proppants is highly dependent on its volumetric concentration and slurry viscosity due to retardation and settling effects. As the location and the concentration of the proppants determine the final fracture permeability, the additional feature is expected to correctly identify the stimulated region. In this paper, the theory and the model formulation are presented along with a few key examples. The simulation can be used to design and optimize the amount of proppant and additives, injection timing, pressure, and well parameters required for successful hydraulic fracturing.


1986 ◽  
Vol 108 (2) ◽  
pp. 107-115 ◽  
Author(s):  
I. D. Palmer ◽  
C. T. Luiskutty

There is a pressing need to compare and evaluate hydraulic fracture models which are now being used by industry to predict variable fracture height. The fractures of concern here are vertical fractures which have a pronounced elongation in the direction of the payzone, i.e., there is a dominant one-dimensional fluid flow along the payzone direction. A summary is given of the modeling entailed in the basic ORU fracture model, which calculates fracture height as a function of distance from the wellbore in the case of a continuous sand bounded by zones of higher (but equal) minimum in-situ stress. The elastic parameters are assumed the same in each layer, and injected flow rates and fluid parameters are taken to be constant. Leak-off is included with spurt loss, as well as non-Newtonian flow. An advantage of the model is its small computer run time. Predictions for wellbore height and pressure from the ORU model are compared separately with the AMOCO and MIT pseudo-3D models. In one instance of high stress contrast the ORU wellbore pressure agrees fairly well with the AMOCO model, but the AMOCO wellbore height is greater by 32 percent. Comparison between the ORU and MIT models in two cases (also high stress contrast) indicates height disagreement at the wellbore by factors of 1.5–2.5 with the MIT model giving a lower height. Thus it appears there can be substantial discrepancies between all three models. Next we compare the ORU model results with six cases of elongated fractures from the TERRA-TEK fully-3D model. Although two of these cases are precluded due to anomolous discrepancies, the other four cases show reasonable agreement. We make a critical examination of assumptions that differ in all the models (e.g., the effective modulus-stiffness multiplier approximation in the AMOCO model, the effect of finite fluid flow in the vertical direction in the MIT model, and the effect of 2D flow and limited perforated height in the TERRA-TEK model). Suggestions are made for reconciling some of the discrepancies between the various models. For example, the ORU/AMOCO height discrepancy appears to be resolved; for other discrepancies we have no explanation. Our main conclusion is that the AMOCO, TERRA-TEK and ORU models for fracture height and bottomhole pressure are in reasonable agreement for highly elongated fractures. Despite the difficulties in understanding the different models, the comparisons herein are an encouraging first step towards normalizing these hydraulic fracture models.


2021 ◽  
Vol 56 (2) ◽  
pp. 164-177
Author(s):  
A. B. Kiselev ◽  
Li Kay-Zhui ◽  
N. N. Smirnov ◽  
D. A. Pestov

Minerals ◽  
2020 ◽  
Vol 10 (8) ◽  
pp. 657
Author(s):  
Chaojie Cheng ◽  
Harald Milsch

Fractures efficiently affect fluid flow in geological formations, and thereby determine mass and energy transport in reservoirs, which are not least exploited for economic resources. In this context, their response to mechanical and thermal changes, as well as fluid–rock interactions, is of paramount importance. In this study, a two-stage flow-through experiment was conducted on a pure quartz sandstone core of low matrix permeability, containing one single macroscopic tensile fracture. In the first short-term stage, the effects of mechanical and hydraulic aperture on pressure and temperature cycles were investigated. The purpose of the subsequent intermittent-flow long-term (140 days) stage was to constrain the evolution of the geometrical and hydraulic fracture properties resulting from pressure solution. Deionized water was used as the pore fluid, and permeability, as well as the effluent Si concentrations, were systematically measured. Overall, hydraulic aperture was shown to be significantly less affected by pressure, temperature and time, in comparison to mechanical aperture. During the long-term part of the experiment at 140 °C, the effluent Si concentrations likely reached a chemical equilibrium state within less than 8 days of stagnant flow, and exceeded the corresponding hydrostatic quartz solubility at this temperature. This implies that the pressure solution was active at the contacting fracture asperities, both at 140 °C and after cooling to 33 °C. The higher temperature yielded a higher dissolution rate and, consequently, a faster attainment of chemical equilibrium within the contact fluid. X-ray µCT observations evidenced a noticeable increase in fracture contact area ratio, which, in combination with theoretical considerations, implies a significant decrease in mechanical aperture. In contrast, the sample permeability, and thus the hydraulic fracture aperture, virtually did not vary. In conclusion, pressure solution-induced fracture aperture changes are affected by the degree of time-dependent variations in pore fluid composition. In contrast to the present case of a quasi-closed system with mostly stagnant flow, in an open system with continuous once-through fluid flow, the activity of the pressure solution may be amplified due to the persistent fluid-chemical nonequilibrium state, thus possibly enhancing aperture and fracture permeability changes.


SPE Journal ◽  
2019 ◽  
Vol 24 (05) ◽  
pp. 2292-2307 ◽  
Author(s):  
Jizhou Tang ◽  
Kan Wu ◽  
Lihua Zuo ◽  
Lizhi Xiao ◽  
Sijie Sun ◽  
...  

Summary Weak bedding planes (BPs) that exist in many tight oil formations and shale–gas formations might strongly affect fracture–height growth during hydraulic–fracturing treatment. Few of the hydraulic–fracture–propagation models developed for unconventional reservoirs are capable of quantitatively estimating the fracture–height containment or predicting the fracture geometry under the influence of multiple BPs. In this paper, we introduce a coupled 3D hydraulic–fracture–propagation model considering the effects of BPs. In this model, a fully 3D displacement–discontinuity method (3D DDM) is used to model the rock deformation. The advantage of this approach is that it addresses both the mechanical interaction between hydraulic fractures and weak BPs in 3D space and the physical mechanism of slippage along weak BPs. Fluid flow governed by a finite–difference methodology considers the flow in both vertical fractures and opening BPs. An iterative algorithm is used to couple fluid flow and rock deformation. Comparison between the developed model and the Perkins–Kern–Nordgren (PKN) model showed good agreement. I–shaped fracture geometry and crossing–shaped fracture geometry were analyzed in this paper. From numerical investigations, we found that BPs cannot be opened if the difference between overburden stress and minimum horizontal stress is large and only shear displacements exist along the BPs, which damage the planes and thus greatly amplify their hydraulic conductivity. Moreover, sensitivity studies investigate the impact on fracture propagation of parameters such as pumping rate (PR), fluid viscosity, and Young's modulus (YM). We investigated the fracture width near the junction between a vertical fracture and the BPs, the latter including the tensile opening of BPs and shear–displacement discontinuities (SDDs) along them. SDDs along BPs increase at the beginning and then decrease at a distance from the junction. The width near the junctions, the opening of BPs, and SDDs along the planes are directly proportional to PR. Because viscosity increases, the width at a junction increases as do the SDDs. YM greatly influences the opening of BPs at a junction and the SDDs along the BPs. This model estimates the fracture–width distribution and the SDDs along the BPs near junctions between the fracture tip and BPs and enables the assessment of the PR required to ensure that the fracture width at junctions and along intersected BPs is sufficient for proppant transport.


1982 ◽  
Vol 22 (03) ◽  
pp. 321-332 ◽  
Author(s):  
M.E. Hanson ◽  
G.D. Anderson ◽  
R.J. Shaffer ◽  
L.D. Thorson

Abstract We are conducting a U.S. DOE-funded research program aimed at understanding the hydraulic fracturing process, especially those phenomena and parameters that strongly affect or control fracture geometry. Our theoretical and experimental studies consistently confirm the well-known fact that in-situ stress has a primary effect on fracture geometry, and that fractures propagate perpendicular to the least principal stress. In addition, we find that frictional interfaces in reservoirs can affect fracturing. We also have quantified some effects on fracture geometry caused by frictional slippage along interfaces. We found that variation of friction along an interface can result in abrupt steps in the fracture path. These effects have been seen in the mineback of emplaced fractures and are demonstrated both theoretically and in the laboratory. Further experiments and calculations indicate possible control of fracture height by vertical change in horizontal stresses. Preliminary results from an analysis of fluid flow in small apertures are discussed also. Introduction Hydraulic fracturing and massive hydraulic fracturing (MHF) are the primary candidates for stimulating production from tight gas reservoirs. MHF can provide large drainage surfaces to produce gas from the low- permeability formation if the fracture surfaces remain in the productive parts of the reservoir. To determine whether it is possibleto contain these fractures in the productive formations andto design the treatment to accomplish this requires a much broader knowledge of the hydraulic fracturing process. Identification of the parameters controlling fracture geometry and the application of this information in designing and performing the hydraulic stimulation treatment is a principal technical problem. Additionally, current measurement technology may not be adequate to provide the required data. and new techniques may have to be devised. Lawrence Livermore Natl. Laboratory has been conducting a DOE-funded research program whose ultimate goal is to develop models that predict created hydraulic fracture geometry within the reservoir. Our approach has been to analyze the phenomenology of the fracturing process to son out and identify those parameters influencing hydraulic fracture geometry. Subsequent model development will incorporate this information. Current theoretical and stimulation design models are based primarily on conservation of mass and provide little insight into the fracturing process. Fracture geometry is implied in the application of these models. Additionally, pressure and flow initiation in the fractures and their interjection with the fracturing process is not predicted adequately with these models. We have reported previously on some rock-mechanics aspects of the fracturing process. For example, we have studied, theoretically and experimentally, pressurized fracture propagation in the neighborhood of material interfaces. Results of interface studies showed that natural fractures in the interfacial region negate any barrier effect when the fracture is propagating from a lower modulus material toward a higher modulus material. On the other hand, some fracture containment could occur when the fracture is propagating from a higher modulus into a lower modulus material. Effect of moduli changes on the in-situ stress field have to be taken into consideration to evaluate fracture containment by material interfaces. Some preliminary analyses have been performed to evaluate how stress changes when material properties change, but we have not evaluated this problem fully. SPEJ P. 321^


Sign in / Sign up

Export Citation Format

Share Document