skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Dispersion of solids in fracturing flows of yield stress fluids
Solids dispersion is an important part of hydraulic fracturing, both in helping to understand phenomena such as tip screen-out and spreading of the pad, and in new process variations such as cyclic pumping of proppant. Whereas many frac fluids have low viscosity, e.g. slickwater, others transport proppant through increased viscosity. In this context, one method for influencing both dispersion and solids-carrying capacity is to use a yield stress fluid as the frac fluid. We propose a model framework for this scenario and analyse one of the simplifications. A key effect of including a yield stress is to focus high shear rates near the fracture walls. In typical fracturing flows this results in a large variation in shear rates across the fracture. In using shear-thinning viscous frac fluids, flows may vary significantly on the particle scale, from Stokesian behaviour to inertial behaviour across the width of the fracture. Equally, according to the flow rates, Hele-Shaw style models give way at higher Reynolds number to those in which inertia must be considered. We develop a model framework able to include this range of flows, while still representing a significant simplification over fully three-dimensional computations. In relatively straight fractures and for fluids of moderate rheology, this simplifies into a one-dimensional model that predicts the solids concentration along a streamline within the fracture. We use this model to make estimates of the streamwise dispersion in various relevant scenarios. This model framework also predicts the transverse distributions of the solid volume fraction and velocity profiles as well as their evolutions along the flow part.  more » « less
Award ID(s):
1554044 2050396
PAR ID:
10427377
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
830
ISSN:
0022-1120
Page Range / eLocation ID:
93 to 137
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Summary Slickwater fracturing has become one of the most leveraging completion technologies in unlocking hydrocarbon in unconventional reservoirs. In slickwater treatments, proppant transport becomes a big concern because of the inefficiency of low-viscosity fluids to suspend the particles. Many studies have been devoted to proppant transport experimentally and numerically. However, only a few focused on the proppant pumping schedules in slickwater fracturing. The impact of proppant schedules on well production remains unclear. The goal of our work is to simulate the proppant transport under real pumping schedules (multisize proppants and varying concentration) at the field scale and quantitatively evaluate the effects of proppant schedules on well production for slickwater fracturing. The workflow consists of three steps. First, a validated 3D multiphase particle-in-cell (MP-PIC) model has been used to simulate the proppant transport at real pumping schedules in a field-scale fracture (180-m length, 30-m height). Second, we applied a propped fracture conductivity model to calculate the distribution of propped fracture width, permeability, and fracture conductivity. In the last step, we incorporated the fracture geometry, propped fracture conductivity, and the estimated unpropped fracture conductivity into a reservoir simulation model to predict gas production. Based on the field designs of pumping schedules in slickwater treatments, we have generated four proppant schedules, in which 100-mesh and 40/70-mesh proppants were loaded successively with stair-stepped and incremental stages. The first three were used to study the effects of the mass percentages of the multisize proppants. From Schedules 1 through 3, the mass percentage of 100-mesh proppants is 30, 50, and 70%, respectively. Schedule 4 has the same proppant percentage as Schedule 2 but has a flush stage after slurry injection. The comparison between Schedules 2 and 4 enables us to evaluate the effect of the flush stage on well production. The results indicate that the proppant schedule has a significant influence on treatment performance. The schedule with a higher percentage of 100-mesh proppants has a longer proppant transport distance, a larger propped fracture area, but a lower propped fracture conductivity. Then, the reservoir simulation results show that both the small and large percentages of 100-mesh proppants cannot maximize well production because of the corresponding small propped area and low propped fracture conductivity. Schedule 2, with a median percentage (50%) of 100-mesh proppants, has the highest 1,000-day cumulative gas production. For Schedule 4, the flush stage significantly benefits the gas production by 8.2% because of a longer and more uniform proppant bed along the fracture. In this paper, for the first time, we provide both the qualitative explanation and quantitative evaluation for the impact of proppant pumping schedules on the performance of slickwater treatments at the field scale by using an integrated numerical simulation workflow, providing crucial insights for the design of proppant schedules in the field slickwater treatments. 
    more » « less
  2. Slickwater hydraulic fracturing is becoming a prevalent approach to economically recovering shale hydrocarbon. It is very important to understand the proppant’s transport behavior during slickwater hydraulic fracturing treatment for effective creation of a desired propped fracture geometry. The currently available models are either oversimplified or have been performed at limited length scales to avoid high computational requirements. Another limitation is that the currently available hydraulic fracturing simulators are developed using only single-sized proppant particles. Motivated by this, in this work, a computationally efficient, three-dimensional, multiphase particle-in-cell (MP-PIC) model was employed to simulate the multi-size proppant transport in a field-scale geometry using the Eulerian–Lagrangian framework. Instead of tracking each particle, groups of particles (called parcels) are tracked, which allows one to simulate the proppant transport in field-scale geometries at an affordable computational cost. Then, we found from our sensitivity study that pumping schedules significantly affect propped fracture surface area and average fracture conductivity, thereby influencing shale gas production. Motivated by these results, we propose an optimization framework using the MP-PIC model to design the multi-size proppant pumping schedule that maximizes shale gas production from unconventional reservoirs for given fracturing resources. 
    more » « less
  3. Abstract The viscosity of fluids and their dependence on shear rate, known as shear thinning, plays a critical role in applications ranging from lubricants and coatings to biomedical and food-processing industries. Traditional models such as the Carreau and Eyring theories offer competing explanations for shear-thinning behavior. The Carreau model attributes viscosity reduction to molecular distortions, while the Eyring model describes shear thinning as a stress-induced transition over an activation energy barrier. This work proposes an extended-Eyring model that incorporates stress-dependent activation volumes, bridging key aspects of both theories. In modifying transition-state theory by using an Evans-Polanyi perturbation analysis, we derive a generalized viscosity equation that accounts for the molecular-scale rearrangements governing fluid flow. The model is validated against computational and experimental data, including shear-thinning behavior of pure squalane and polyethylene oxide (PEO) aqueous solutions. Comparative analysis with Carreau-Yasuda and conventional Eyring models demonstrates excellent accuracy in predicting viscosity trends over a wide range of shear rates. The introduction of stress-dependent activation volumes provides a description of molecular exchange kinetics accounting for structural reorganization under shear. These findings offer a unified framework for modeling shear thinning and have broad implications for designing advanced lubricants, polymer solutions, and complex fluids with tailored flow properties. Graphical Abstract 
    more » « less
  4. We investigate the rheological behavior of athermal particle suspensions using experiments and theory. A generalized version of the homogenization estimates of Ponte Castañeda and Willis [J. Mech. Phys. Solids, 43(12), 1919–1951 (1995)] is presented for the effective viscosity of athermal suspensions accounting for additional microstructural features (e.g., polydispersity) via an empirical parameter, [Formula: see text]. For the case of identically sized spheres dispersed with statistical isotropy in a Newtonian fluid, the parameter [Formula: see text] is estimated from the results of Batchelor and Green [J. Fluid Mech. 56(2), 375–400 (1972)] for the Huggins coefficient. Predictions for the macroscopic viscosity are found to be in good agreement with measurements for monodisperse polymethyl methacrylate (PMMA) spheres in glycerol, as well as for the empirical Krieger–Dougherty equation for the shear viscosity. The proposed estimates have the added benefit that they can also be used to get information on the statistics of the stress and strain-rate fields in the fluid and particle phases. In addition, results for the effective shear viscosity are used in combination with the linear comparison method of Ponte Castañeda [J. Mech. Phys. Solids 39(1), 45–71 (1991)] to generate the corresponding estimates for the effective macroscopic behavior and field statistics of particle suspensions in (viscoplastic) yield stress fluids. Good agreement is also found between the theoretical estimates and experimental results for the effective yield and flow stress of suspensions with monodisperse PMMA spheres in Carbopol. Finally, it is argued that the results for the phase averages and fluctuations of the stress and strain-rate fields can be used to provide a physical interpretation for the parameter [Formula: see text] in terms of the polydispersity of the suspension and its implications for the percolation threshold. 
    more » « less
  5. Bacterial suspensions—a premier example of active fluids—show an unusual response to shear stresses. Instead of increasing the viscosity of the suspending fluid, the emergent collective motions of swimming bacteria can turn a suspension into a superfluid with zero apparent viscosity. Although the existence of active superfluids has been demonstrated in bulk rheological measurements, the microscopic origin and dynamics of such an exotic phase have not been experimentally probed. Here, using high-speed confocal rheometry, we study the dynamics of concentrated bacterial suspensions under simple planar shear. We find that bacterial superfluids under shear exhibit unusual symmetric shear bands, defying the conventional wisdom on shear banding of complex fluids, where the formation of steady shear bands necessarily breaks the symmetry of unsheared samples. We propose a simple hydrodynamic model based on the local stress balance and the ergodic sampling of nonequilibrium shear configurations, which quantitatively describes the observed symmetric shear-banding structure. The model also successfully predicts various interesting features of swarming vortices in stationary bacterial suspensions. Our study provides insights into the physical properties of collective swarming in active fluids and illustrates their profound influences on transport processes. 
    more » « less