During fracturing operations in subterranean formations for the production of hydrocarbons, fractures are generated in a zone of interest by pumping a fluid at relatively high pressures that create and extend fractures in the rock. Induced formation fractures may then form channels that permit the return of hydrocarbon fluids to the wellbore. While pressurized fluids may be used for the fracturing process, injected fluids may also serve additional functions, including transporting particulate proppants that prevent the closure of induced and existing fractures. In order to adjust the solids transport properties of fracturing fluids, fluid pumping rates may be adjusted to increase transport capacity. Further, rheology may be modified using various viscosifiers and flow improvers that increase suspension capacity of the fracturing treatments and control the extent that fractures are propagated from the wellbore. Moreover, fracturing fluid design may also strike a balance between the requirements for the suspension properties for various solids and maintaining a fluid viscosity that is within the capacity of the pumping infrastructure at the well site.
This summary is provided to introduce a selection of concepts that are described further below in the detailed description. This summary is not intended to identify key or essential features of the claimed subject matter, nor is it intended to be used as an aid in limiting the scope of the claimed subject matter.
In one aspect, embodiments of the present disclosure may be directed to methods that include using geometrical and mechanical properties to generate a model of at least one interval of a wellbore; designing a fluid pumping schedule for a wellbore treatment fluid system; simulating flow of a fluid system containing solids using the model of at least one interval of the wellbore, wherein simulating comprises determining the flow rate distribution within the wellbore and modeling the settling and resuspension of the solids within the flow of the fluid system; and updating the fluid pumping schedule based on the output determined from the simulated flow of the fluid system.
In another aspect, embodiments of the present disclosure may be directed to methods that include using geometrical and mechanical properties to generate a model of at least one interval of a wellbore; simulating flow of a fluid system containing at least one of solid and additive within at least one interval of a wellbore, wherein simulating comprises determining the flow rate distribution within the wellbore and modeling settling and resuspension of the solids within the fluid system; determining the critical velocity of the fluid system for various concentrations of at least one of solid and additive; and selecting the concentration of at least one of solid and additive based on the critical velocities determined.
Other aspects and advantages of the disclosure will be apparent from the following description and the appended claims.
Embodiments disclosed herein are directed to the design of treatment pumping schedules used to deliver fluids and solids into a wellbore with control over solid settling. In one or more embodiments, methods in accordance with the present disclosure may be used to design fracturing treatments and optimize delivery of proppant for operations in hydrocarbon-bearing formations. As used herein, “fracturing” and “fracturing treatment” may be used to refer to types of fracturing operations, including the initial fracturing of one or more intervals in a wellbore, the treatment of existing fractures, and refracturing operations.
Multi-stage hydraulic fracturing technology may be used for commercial hydrocarbons production in unconventional resources to increase cumulative production and life-time of the well without drilling new wells. In addition, refracturing techniques may be employed in primary and secondary recovery operations to increase hydrocarbon production. Hydrocarbon extraction from a reservoir may involve injecting treatment fluids containing specialized additives such as proppants that are transported through the horizontal sections of a wellbore and delivered to identified intervals (or clusters) to prop open formation fractures created or reopened by the force exerted from the injected fluids. With particular respect to
During fracturing operations, flow rates may approach zero as injected fluids distribute among one or more clusters and into the fracture network, despite high rates of fluid delivery. Leak off of fluids into fracture in the formation may lead to the accumulation of solid particles at various intervals in the wellbore and the formation of blockages and plugs in some cases. Low flow rates and poor rheology of the transport fluid may make delivery of treatment fluids and solids to distant clusters and at long horizontal distances unmanageable.
Maintaining flow rates above a critical flow rate threshold, the rate at which solids begin to settle out and accumulate, or the “critical velocity,” becomes substantial factor in the design of treatment pumping schedules. For example, reduced fluid rates are design considerations for refracturing operations in which one or more intervals within the wellbore are perforated and connected to a large number of fractures created by the first fracturing treatment and considerable fluid loss and diversion is expected.
Methods in accordance with the present disclosure include the design of fluid treatments and pumping schedules that minimize the risk of solids settling and generation of barriers to adequate wellbore flow and plugging. In some embodiments, methods may include multi-stage fracturing and/or refracturing operations that control and/or minimize the settling of solids in a wellbore.
In one or more embodiments, the design of fluid treatments may include generating one or more models describing multicomponent slurry flow in a wellbore. Models in accordance with the present disclosure may provide measures for preventing solids settling that may include: pumping rate adjustments that ensure that the flow velocity is high enough to reduce settling; the addition of viscosifying agents to injected fluids; the addition of insoluble materials; and energizing the fluid with a gaseous phase such as nitrogen and carbon dioxide.
The following section describes an example of a workflow in accordance with the present disclosure that may be used to generate an optimized fluid treatment pumping schedule.
Treatment Design Workflow
Methods in accordance with the present disclosure may incorporate a workflow to design a wellbore treatment. With particular respect to
Pumping schedules generated may minimize the risk of negative outcomes during fracturing and refracturing operations such as wellbore plugging. Methods in accordance with the present disclosure may also be used as a part of workflow for multi-criteria optimization operations such as production simulation and calculation of net present value of a selected reservoir. In some embodiments, risk assessments may be based on simulation results obtained in a model of slurry flow in the wellbore.
Models in accordance with the present disclosure may be described as lists of input and output parameters that are discussed in further detail below.
Input Data
Inputs for models may include a pumping schedule that incorporates a list of stages characterized by fluid pumping rate and time, and concentrations of fluid components. Other model inputs include geometrical and mechanical properties of the wellbore and piping in contact with the pumped fluids. Geometrical properties may include hole diameter, interval lengths, depth, spatial orientation, and the like. Relevant mechanical properties may include standard geological and/or mechanical characterizations such as stresses, Young's modulus, Poisson's ratio, and the like.
Other model inputs include structural descriptors of the wellbore and downhole equipment in previously fractured or unconsolidated formations, which may include properties of perforated intervals such as cluster spacing, perforation diameters, perforation numbers, and the like. Models in accordance with the present disclosure may also account for properties of various fluid components of the anticipated fluid to be injected such as properties of any solids including mean particle diameter, specific gravity, overall particle geometry, and the like; and properties of fluids such as rheology, additive content, specific gravity, temperature, and the like.
Output Data
Models in accordance with the present disclosure may generate one or more outputs that may be used to generate a final pumping schedule. In one or more embodiments, outputs may include the cross-sectional area occupied by solids precipitated from an injected fluid, such as solid particle bed height, at each point of wellbore trajectory, and may also include the flow rates at each point of a wellbore trajectory and at all perforations. In some embodiments, models may calculate the defined as the critical velocity of a given fluid system, which is defined as the flow velocity below which the solids settle forming a bank at the wellbore bottom. Models in accordance with the present disclosure may also be constructed to account for the coupled problems of calculating flow rates distribution, and modeling of solids settling and resuspension.
In one or more embodiments, construction of models that indicate conditions within the wellbore and/or during simulated fluid flow may involve the use of 3D modeling techniques. For example, models and simulations of inflow into fractures may involve the use of Khristianovich-Geertsma-de Klerk (KGD), Perkins-Kern-Nordgren (PKN), radial, pseudo 3D, planar 3D, or full 3D approaches. Simulation of proppant transport may involve the use of commercially available software packages such as ANSYS® CFD or Star-CCM+.
In one or more embodiments, fluid treatments in accordance with the present disclosure may mitigate solids settling through a number of mechanisms including adjusting pumping rates of fluids into the wellbore and incorporating various additives into injected treatment fluids. Each of these mechanisms will be discussed in turn below.
Pumping Rate Adjustment
Methods in accordance with the present disclosure may incorporate one or more multicomponent slurry flow models to calculate the pumping rate in one or more stages of a fluid treatment pumping schedule. During fluid flow, a slurry composition containing solids may undergo various stages of solids transport. With particular respect to
In some embodiments, decreasing fluid flow velocity may alter the flow of a solids-containing fluid from a heterogeneous flow to “sliding bed” flow regime. To avoid this phenomenon, transition velocity may be held at a rate higher than the critical velocity. Flow models used to determine critical velocity include Turian's model as shown in Eq. 1, where Vcrit is the critical velocity, φ is the volume fraction of solids in the flow, dp is the particle diameter, D is the pipe diameter, s is the ratio of solid density to fluid density, μ is the viscosity of the fluid, and ρ is particle density.
Turian model represented by Eq. 1 was created for Newtonian solid-laden fluids. For non-Newtonian fluids, no calibrated models currently exist. However, additives may be combined with the fluid system to shift the critical velocity to lower values.
In one or more embodiments, the ability of different materials to decrease critical velocity may be determined empirically through testing of various additive concentrations using a flow loop apparatus or other similar testing device. In some embodiments, the critical velocity of various treatment fluids may be analyzed using a flow loop setup that combines fluids and solids and measures fluid flow and solids transport properties. With particular respect to
In some embodiments, testing methods in accordance with the present disclosure may include mixing a fluid composition such as water, slickwater, or linear gels with the selected additives for 5 to 30 minutes in a hopper achieve the proper level of dispersion. During the dispersion stage, materials may be added to the fluid composition using the by-pass line 104. The hopper 402 may then be charged with one or more solids to achieve a final circulating concentration of 0.1 to 10 pounds of proppant added (ppa) in the pumped fluid. In some embodiments, testing runs may be performed with a solid additive concentration within the range of 5 to 10 ppa. The flow loop apparatus is then operated at a pumping regime that provides fluid flow approximating the critical velocity rate through the main line. The rate of fluid flow is then decreased incrementally and the critical rate at which sand begins to settle is recorded along with the height of any generated proppant dune.
In addition to measuring the critical velocity of the system, solids settling may be studied by a number of additional techniques that include calculating the velocity distribution along the lateral and amount of proppant occupying a given cluster. In order to validate critical velocity, methods in accordance with the present disclosure may incorporate a control run in which the fluid composition is circulated throughout the flow loop 406 without additives to establish a baseline critical velocity for a fluid system.
In one or more embodiments, the concentration of additive that decreases the critical velocity of the fluid through the piping is determined and recorded. With particular respect to
During testing of a fluid system, the additive concentration, and effectively the critical velocity of the fluid, may then be adjusted using this method to fall within the desired operational parameters. The model is then executed with the newly determined critical velocity and estimates are generated for cross-sectional areas of the system that will experience solids settling. If the outcome is less than desired, the concentration of flow modifiers is altered or differing solid materials are selected, and the method is repeated until a satisfactory treatment composition and pumping schedule is generated.
Additive Materials
In one or more embodiments, additives that may be combined with fluids in accordance include solids and other materials that are often regarded as proppants. The geometry of additives may include materials such as fibers, slices, flakes, fines, colloidal particles, and pulps. Additive materials in accordance the present disclosure may have varying geometries including spherical, oblate, or fibrous and may vary in diameter (or length for materials having an asymmetrical aspect ratio) from 10 nm to 10 mm.
In some embodiments, additives may be rigid or flexible, straight or curved, and single or multi-component. Additives in accordance with the present disclosure may be be any material, such as comminuted plant materials, synthetic polymers such as polyester, polyaramide, polyamide, novoloid or novoloid-type polymers, fibrillated synthetic organic polymers, ceramics, inorganics, metals, carbon, glass, natural polymer, and mixtures thereof. In one or more embodiments, additives may include cellulose, polylactic acid polyester, polyglycolic acid polyester, polyvinyl alcohol, polyamide, and the like.
In one or more embodiments, fluid systems in accordance with the present disclosure may be modified using one or more viscosifying agents. Viscosifying agents may include polymers such as hydroxyethylcellulose, carboxymethylguar, carboxymethylhydroxypropyl guar, and various viscosoelastic surfactants. In some embodiments, polymeric viscosifying agents may be applied to a fluid system as a linear additive or as crosslinked and/or graft polymers. Viscosifying additives may also include clays such as bentonite and attapulgite, or any combination thereof.
In one or more embodiments, flow models may be developed for various compositions of matter including fluids systems containing solids, polymers, or gasses. For example, fluid flow models may be applied to energized fluids such as foaming fluids and stabilized foamed fracturing fluids.
Modeling Fluid Flow and Solids Transport in the Formation
In one or more embodiments, methods in accordance with the present disclosure may incorporate one or more models of solids flow in a wellbore to design pumping schedules for fracturing treatments. In some embodiments, models may account for combinations of problems encountered in the design of pumping schedules, including the problems of (1) calculating flow rate distribution, and (2) modeling solids settling and resuspension. The selected problems will be discussed in turn in the following sections.
Problem 1: Calculation of Flow Rates Distribution
With particular respect to
In one or more embodiments, algorithms for designing fluid treatments may also account for other flow-modifying effects that may include casing friction, hydrostatic forces, perforation friction, and fracture presence. Other known pressure drops in the system may be added to the model in some embodiments.
The casing friction in a wellbore element may be determined in some embodiments by Eq. 2, where c is the friction pressure coefficient, d is the friction pressure constant, qwb is the flow rate contributed by the wellbore, and ΔL is the length of the wellbore element.
ΔPfrict=c·qWBd·ΔL (2)
The hydrostatic pressure drop is given by Eq. 3, where Δzi is the height difference between reference point A and an i-th perforated interval, g is the local acceleration due to gravity, and ρslurry is the density of the slurry composition.
ΔPhydr,i=ρslurry·g·Δzi (3)
The perforation friction is given by Eq. 4, where qi is the fluid flow into the i-th fracture, mi is the number of perforations at i-th cluster, dperf is the diameter of perforation, and Cd is a dimensionless discharge coefficient.
Eqs. 2-4 couple an analytical model of a constant height PKN fracture to a set of equations describing conservation of slurry volume and pressure continuity. These equations provide coupling between multiple fractures via the wellbore. The net pressure corresponds to a Perkins-Kern-Nordgren (PKN) fracture for the zero leak-off case, where ν is Poisson's ratio, μ is the fluid viscosity, E is Young's modulus, h is the fracture height, and t is time since fracture initiation.
Variable qfrac is defined as shown in Eq. 6, wherein qi is the fluid flow into the i-th fracture.
While one example is shown above, any system simulating a hydraulic fracture could be used in methods in accordance with the present disclosure to calculate the net pressure caused by presence of the fracture.
One assumption in the above algorithm is the absence of fluid transfer between the fractures within the formation, which implies that the fractures are coupled through the wellbore. With this assumption the conservation of fluid volume may be described the equality shown in Eq. 7 stating that the total pumping rate Q is equal to the sum of each fluid flow into the i-th fracture qi.
Q=Σ
i=1
n
q
i (7)
Continuity of pressure implies that the sum of all pressure changes is equal to the reference pressure for any path connecting a fracture tip of the i-th fracture with the reference point A as described in Eq. 8 and as shown in
P
ref=Σk=1iΔPfrick,k+ΔPperf,i−ΔPhydr,i,i=1,n (8)
In one or more embodiments, the solution to the system of non-linear equations Eq. 7 and Eq. 8, where qi, . . . , qn, Pref are unknowns, may be found by the application of a Newton-like method for nonlinear system of equations.
Problem 2: Modeling of Solids Settling and Resuspension
Another component of models in accordance with the present disclosure is modeling the process of solid suspension in fluid flows. In order to study the phenomenon, an experimental flow loop is constructed from a horizontal pipe with n perforated intervals, which allow for fluid communication with a formation. The flow of slurry in the wellbore splits in multiple streams, an immobile bank 704 and a flowing slurry 702, as shown in
From this model, the average slurry velocity may be calculated as shown in Eq. 9, where Vslurry is the velocity of the slurry, qwb is the flow rate in the wellbore element, Spipe is the cross-sectional area of the pipe, and Sbank is the cross-sectional area of the bank created by solids settling.
In one or more embodiments, bank area Sbank may be calculated by determining if Vslurry is greater than the minimal transport velocity V1, at which point Sbank may be assumed to be zero. If Vslurry is less than V1, the area of bank may increase until Vslurry is equivalent to Vi, which may be described by Eq. 10.
In one or more embodiments, Eqs. 11 and 12 may be used to calculate the minimum flow rate required to transport solid particles. In some embodiments, the velocity of course particles in a fluid system may be described by the Durand equation as shown in Eq. 11, where VDurand is the minimum transport velocity, F is an empirical constant that varies between 0.4 and 1.5, s is the ratio of the densities of the particle component and the fluid component, g is the acceleration due to gravity, and D is the wellbore diameter.
V
Durand
=F√{square root over (2·g·(s−1)·D)} (11)
The Wasp equation shown in Eq. 12 may also be used as a correction to Eq. 11 that takes into account the particles' size, where dp is the diameter of the solid particles.
V
Wasp
=F√{square root over (2·g·(s−1)·D)}·(dp)1/6 (12)
In one or more embodiments, the problems of flow rate calculation as distributed throughout the fractures and solid settling from flowing fluids may be combined in a coupled Algorithm. For example, a two-way coupled algorithm may be used to solve problems 1 and 2 iteratively (as discussed in
The following examples, methods in accordance with the present disclosure are used to determine the settling rate of solids for selected systems.
In the first example, a method in accordance with the present disclosure is applied to overcome wellbore sanding, in a two cluster interval located in a formation containing layers of varying stresses. The wellbore under study had a horizontal section with an inner diameter (ID) of 5 inches. The lateral section of the well was 300 ft in length with 2 perforation clusters equally spaced, with each perforation cluster contained 6 perforations. The perforation diameter (dperf) was 0.38 inches, with a discharge coefficient (Cd) of 0.85. The formation was characterized as having a Young's modulus (E) of 4 Mpsi, and a Possion's ratio (ν) of 0.2. The height of both fractures is 50 ft. The stimulation uses slickwater with dynamic viscosity (μ) of 5 cP and 100 mesh sand having an average particle diameter (dp) of 0.15 mm with a ratio of particle and fluid densities (s) of 2.65.
The minimum transport velocity in this case is 3.24 ft/s, which corresponds to flow rate 4.7 barrel per minute (BPM) in the wellbore. For this example, the pumping rate into the formation for the samples varied from 0 BPM up to 40 BPM. The closure stress at the first interval was 6000 psi and 6700 psi at the second interval. The results of the simulation are presented in Table 1.
As demonstrated in Table 1, the first cluster of the two cluster system is the only cluster open at pumping rates of less than 10 BPM. Further, when the pumping rate is below 5 BPM, the cross-sectional area occupied by sand is greater than zero, indicating a pumping regime that is unsafe. Pumping rates in the range of 5 to 9 BPM showed no evidence of sand deposition and were regarded as a safe pumping regime. Pumping rates in the range of 10 to 35 BPM were sufficient to open the second cluster, but also produced a sand bank between first and second cluster intervals and were regarded as unsafe. On the other hand, pumping rates greater that 40 BPM failed to produce sand bank deposition, indicating a safe pumping regime.
In the next example, the developed algorithm is used to solve another set of well conditions within a horizontal section having an ID of 5 inches. The lateral was 500 feet in length with 10 perforation clusters equally spaced and a total number of perforations of 60. Diameter dperf was 0.38 inches under a Cd of 0.85. The formation under study possessed an E of 5 Mpsi and a ν of 0.2. The height of all the fractures was set to 50 ft. The injected fluid system was slickwater with a μ of 2 cP and contained 20/40 mesh sand having a dp of 0.5 mm. During testing, the pumping rate varied from 10 BPM up to 90 BPM or 1.5 bpm per perforation. The minimum transport velocity measured was 7 ft/s, which corresponds to a flow rate of 10 BPM in the wellbore.
With particular respect to
Although the preceding description has been described herein with reference to particular means, materials and embodiments, it is not intended to be limited to the particulars disclosed herein; rather, it extends to all functionally equivalent structures, methods and uses, such as are within the scope of the appended claims. In the claims, means-plus-function clauses are intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures. Thus, although a nail and a screw may not be structural equivalents in that a nail employs a cylindrical surface to secure wooden parts together, whereas a screw employs a helical surface, in the environment of fastening wooden parts, a nail and a screw may be equivalent structures.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/RU2016/000290 | 5/16/2016 | WO | 00 |