Herschel-Bulkley Fluid | Oil Authority

A Herschel-Bulkley fluid is a non-Newtonian fluid model that combines a yield stress (a minimum shear stress that must be exceeded before the material begins to flow) with a power-law relationship between shear stress and shear rate above that yield point — described by the equation τ = τ₀ + K × γⁿ, where τ is shear stress, τ₀ is the yield stress, K is the consistency index, γ is the shear rate, and n is the flow behavior index — making it the most general and accurate three-parameter rheological model used in drilling fluid engineering to characterize the behavior of weighted muds, polymer-enhanced fluids, and cement slurries that exhibit both a yield stress threshold and power-law flow characteristics that neither the simpler Bingham plastic nor pure power-law models can adequately describe alone.

Key Takeaways

The three parameters of the Herschel-Bulkley model each capture a physically distinct aspect of drilling fluid behavior: the yield stress τ₀ represents the gel strength that holds cuttings in suspension when circulation stops (critical for preventing barite sag and cuttings settling in inclined wellbores), the consistency index K governs viscosity at intermediate shear rates typical of annular flow, and the flow behavior index n characterizes whether the fluid shear-thins (n less than 1, typical of polymer and weighted muds where viscosity decreases with increasing shear rate) or shear-thickens (n greater than 1, rare in drilling fluids) — together providing a complete description of fluid behavior across the full shear rate range from near-zero in the annulus to high rates at the bit.
The Herschel-Bulkley model is more accurate than the Bingham plastic model (τ = τ₀ + μₚ × γ, a two-parameter model that assumes linear flow above yield stress) for most modern drilling fluids because real weighted muds do not exhibit linear shear stress-shear rate behavior — the power-law exponent n is typically 0.3 to 0.7 for polymer-weighted muds, meaning viscosity decreases significantly with shear rate; using Bingham plastic parameters for hydraulics calculations in these fluids overpredicts pressure losses at low shear rates and underpredicts at high shear rates, leading to incorrect equivalent circulating density calculations and wellbore stability assessments.
Herschel-Bulkley parameters are determined from six-speed viscometer readings (3, 6, 100, 200, 300, 600 RPM) on a Fann VG meter by fitting the three-parameter model to the full shear stress-shear rate dataset using least-squares regression or the API standard three-point fit procedure — unlike Bingham plastic which uses only the 300 and 600 RPM readings, or power-law which ignores yield stress, the Herschel-Bulkley fit uses all six readings to capture the complete rheological profile, particularly the critical low-shear behavior at 3 and 6 RPM that governs cuttings transport in horizontal and deviated wells where gravitational settling is the dominant concern.
Hydraulic pressure loss calculations for Herschel-Bulkley fluids in the annulus and drillstring require iterative numerical solutions because the governing equations (modified Navier-Stokes equations with yield stress and power-law viscosity) do not have closed-form analytical solutions for annular geometry — drilling engineering software (Landmark WELLPLAN, Halliburton WellPlan, Baker Hughes ADVANTAGE) uses numerical integration of the velocity profile across the annulus cross-section at each depth increment to compute friction pressure loss, ECD, and cuttings transport efficiency, making Herschel-Bulkley hydraulics computationally more demanding than Bingham plastic but significantly more accurate for modern fluid systems.
Herschel-Bulkley rheology is particularly critical for cement slurry design in challenging wellbores — oil well cement slurries exhibit yield stress and power-law flow behavior that determines the pump pressure required to displace the slurry to TD, the minimum displacement velocity needed to achieve turbulent flow in the annulus for good mud displacement, and the settling resistance of cement during the thickening time period; API Specification 10A and ISO 10426 cement testing standards include Herschel-Bulkley characterization as the recommended rheological model for cements because Bingham plastic parameters systematically underestimate displacement pressures in deviated casing strings where plug-to-turbulent flow transitions govern cement job quality.

Fast Facts

The Herschel-Bulkley model was published by Winslow Herschel and Ronald Bulkley in 1926 in a paper describing the flow of rubber compounds, predating its widespread application to drilling fluids by several decades. The model gained prominence in petroleum engineering through the work of researchers at the University of Tulsa and oil company research laboratories in the 1970s and 1980s who demonstrated that Bingham plastic parameters systematically mispredicted annular pressure losses in polymer-weighted muds. Today, the Herschel-Bulkley model is the default rheological model in all major drilling engineering software packages and is specified in API RP 13D (Rheology and Hydraulics of Oil-Well Drilling Fluids) as the recommended model for drilling fluid hydraulics calculations when all six viscometer readings are available.

What Is a Herschel-Bulkley Fluid?

Every fluid that flows obeys some relationship between the force applied to make it flow (shear stress) and how fast it moves in response (shear rate). Water and simple oils follow a linear relationship — double the force, double the flow rate — and are called Newtonian fluids. Drilling fluids, cement slurries, and most industrial suspensions are far more complex. They exhibit a yield stress, meaning they behave like a solid at rest and only begin flowing when a threshold force is exceeded. Above that threshold, their viscosity changes with shear rate in a non-linear way.

The Herschel-Bulkley model captures this complexity in a single equation: τ = τ₀ + K × γⁿ. The yield stress τ₀ is what prevents cuttings from settling when the pumps are shut down. The power-law term K × γⁿ describes how the fluid flows once mobilized — for most drilling muds, n is less than 1, meaning the fluid becomes less viscous as shear rate increases (shear-thinning behavior), which is exactly what engineers want: high viscosity at low shear rates in the annulus to carry cuttings, but low viscosity at high shear rates at the bit to minimize hydraulic horsepower requirements.

The genius of the model is that it encompasses the two simpler models as special cases. When n = 1, Herschel-Bulkley reduces to the Bingham plastic model. When τ₀ = 0, it reduces to the pure power-law model. For real drilling fluids that have both yield stress and non-linear flow behavior, neither simplification is adequate, and the full three-parameter Herschel-Bulkley model is required for accurate hydraulics prediction.

Herschel-Bulkley Fluid Applications in Drilling Engineering

Equivalent circulating density (ECD) calculation is the primary application of Herschel-Bulkley rheology in drilling — the friction pressure loss in the annulus adds to the hydrostatic pressure of the mud column to give the ECD, which must be managed between the pore pressure (fractures formation if ECD exceeds fracture gradient) and the collapse pressure (wellbore unstable if ECD falls below minimum stress). Herschel-Bulkley-based ECD calculations are more accurate than Bingham plastic calculations for weighted polymer muds, particularly in horizontal and extended-reach wells where the annular velocity varies substantially along the wellbore and the low-shear yield stress behavior dominates in low-flow-rate sections.

Cuttings transport modeling in deviated and horizontal wells requires accurate characterization of low-shear rheology — in wells deviated beyond 45 degrees, cuttings form a mobile bed on the low side of the wellbore that moves intermittently, and the critical velocity to re-suspend settled cuttings depends on the yield stress and low-shear viscosity of the drilling fluid. Herschel-Bulkley models predict this critical resuspension velocity more accurately than Bingham plastic models because the low-RPM viscometer readings that capture yield stress and near-yield behavior are explicitly included in the three-parameter fit, rather than being represented by a single yield stress parameter that may not capture the shape of the flow curve in the critical 3-100 RPM range.

Barite sag prediction in highly deviated wells uses Herschel-Bulkley parameters to assess whether the yield stress is sufficient to suspend barite (density approximately 4.3 g/cc) against gravity when circulation is interrupted for connections, surveys, or wiper trips. The critical yield stress required to prevent barite settling scales with the particle diameter and density difference, and comparing the measured Herschel-Bulkley yield stress τ₀ against the critical threshold provides a quantitative assessment of sag risk that guides fluid conditioning and weighting agent selection decisions before drilling high-angle sections.

Herschel-Bulkley Fluid Applications Across International Jurisdictions

Canada (AER / WCSB): WCSB horizontal Montney and Duvernay wells use Herschel-Bulkley-characterized synthetic oil-based muds for the 3,000 to 5,000 meter horizontal laterals where ECD management within narrow pressure windows between pore pressure and fracture gradient is critical. AER Directive 008 well licensing requirements implicitly require that drilling programs demonstrate the drilling fluid system can maintain wellbore stability throughout the planned trajectory, and Herschel-Bulkley hydraulics models are the basis for the ECD calculations submitted in drilling programs for HPHT or narrow-margin wells on the WCSB. Suncor, ConocoPhillips, and Shell Canada report using Herschel-Bulkley-based hydraulics software for Montney horizontal development well design.

United States (API / BSEE): API RP 13D (Rheology and Hydraulics of Oil-Well Drilling Fluids) is the foundational US standard for drilling fluid rheology and explicitly recommends the Herschel-Bulkley model for all hydraulics calculations when six-speed viscometer readings are available. Gulf of Mexico deepwater drilling programs use Herschel-Bulkley-characterized synthetic muds where the narrow ECD window in deepwater (between the weak deepwater sediment fracture gradient and the pore pressure) requires the highest accuracy hydraulics models available. BSEE well control regulations for HPHT wells require documentation of the drilling fluid hydraulics analysis that forms the basis for casing design and BOP selection, and Herschel-Bulkley models are standard in the hydraulics software used for these submissions.

Norway (Sodir / NORSOK): NORSOK D-010 well integrity standards for NCS drilling operations specify that hydraulics calculations must use a rheological model that accurately represents the full shear rate behavior of the drilling fluid — in practice this means Herschel-Bulkley for all modern polymer and oil-based muds used on the NCS. Equinor's internal drilling engineering standards reference Herschel-Bulkley hydraulics as the minimum acceptable approach for well designs with ECD margins less than 0.5 ppg. The NCS's deepwater frontier wells (Barents Sea, Norwegian Trench) use Herschel-Bulkley-based real-time hydraulics monitoring to manage ECD within the narrow pressure window between shallow gas formation pore pressures and weak overburden fracture gradients.

Middle East (Saudi Aramco): Saudi Aramco's maximum reservoir contact horizontal wells in the Arab Formation use Herschel-Bulkley-characterized non-aqueous drilling fluids selected for ECD management in the long horizontal sections where the hydrostatic pressure gradient must be maintained above minimum horizontal stress to prevent wellbore breakout while remaining below the natural fracture opening pressure that would cause lost circulation. Aramco's drilling engineering department uses Herschel-Bulkley parameters measured on-site at field conditions (temperature-corrected using the HPHT viscometer) rather than surface measurements, because Arab Formation well temperatures of 80 to 120°C significantly alter the rheological parameters relative to surface values and Bingham plastic interpolations fail to capture the temperature dependence of yield stress and consistency in weighted PHPA-polymer muds.