⛽ Comprehensive Drilling Calculator

180 Formulas · Field / SI / Metric · Hydrostatics · Well Control · Hydraulics · Fluids · Mechanical · Directional

Hydrostatics & Fundamental Pressure 22 Formulas

1. Pressure Gradient

psi/ft = MW(ppg) × 0.052kPa/m = ρ(kg/m³) ÷ 102kg/cm²/m = ρ(kg/l) ÷ 10
—

Field: psi/ft · SI: kPa/m · Metric: kg/cm²/m

2. Hydrostatic Pressure

HP(psi) = MW × 0.052 × TVD(ft)HP(kPa) = ρ(kg/m³) × 0.00981 × TVD(m)HP(bar) = ρ(kg/l) × 0.0981 × TVD(m)
—

3. Equivalent Mud Weight from LOT/FIT

EMW(ppg) = MW_LOT + LOT_psi ÷ 0.052 ÷ TVD_shoeEMW(kg/m³) = ρ_LOT + LOT_kPa × 102 ÷ TVD_shoeEMW(kg/l) = ρ_LOT + LOT_bar × 10.2 ÷ TVD_shoe
—

4. Specific Gravity from Mud Density

SG = MW(ppg) ÷ 8.33SG = ρ(kg/m³) ÷ 1000SG = ρ(kg/l)
—

5. Formation Pressure

Pf(psi) = Ph_DS(psi) + SIDPP(psi)Pf(kPa) = Ph_DS(kPa) + SIDPP(kPa)Pf(bar) = Ph_DS(bar) + SIDPP(bar)
—

Also: MW required = (Pf ÷ 0.052 ÷ TVD) ppg

6. Bottom-Hole Pressure — Shut-In

BHP = Ph_DS + SIDPP = Ph_Ann + SICPBHP(kPa) = Ph_DS + SIDPPBHP(bar) = Ph_DS + SIDPP
—

Verify: Ph_Ann + SICP should equal Ph_DS + SIDPP

7. BHP While Circulating

BHP_circ = Ph_Ann + APL (psi)BHP_circ = Ph_Ann + APL (kPa)BHP_circ = Ph_Ann + APL (bar)
—

8. BHP — Killing Well (Subsea)

BHP_kill = Ph_Ann + P_Ann + CLFL (psi)BHP_kill(kPa) = Ph_Ann + P_Ann + CLFLBHP_kill(bar) = Ph_Ann + P_Ann + CLFL
—

9. Equivalent Circulating Density (ECD)

ECD(ppg) = MW + APL ÷ 0.052 ÷ TVDECD(kg/m³) = ρ + APL(kPa) × 102 ÷ TVD(m)ECD(kg/l) = ρ + APL(bar) × 10.2 ÷ TVD(m)
—

10. Formation Temperature

FT(°F) = Surface_Temp + (0.012°F/ft × TVD)FT(°C) = Surface_Temp + (0.022°C/m × TVD)FT(°C) = Surface_Temp + (0.022°C/m × TVD)
—

11. Temperature Conversions

°R = °F + 460 · °K = °C + 273°K = °C + 273 · °F = (°C × 9/5) + 32°K = °C + 273
—

12. Crude Oil Density from °API

SG = 141.5 ÷ (131.5 + °API) · Density(ppg) = SG × 8.33
—

13. Gas Pressure at Depth (Charles/Boyle)

P_d = P_s × e^(0.2085 × Gg × D ÷ (T_avg + 460))P_d = P_s × e^(0.03415 × Gg × D_m ÷ (T_avg + 273))P_d = P_s × e^(0.03415 × Gg × D_m ÷ (T_avg + 273))
—

14. Gas Gradient

Gas Grad (psi/ft) = (P_deep − P_shallow) ÷ Depth_intervalGas Grad (kPa/m) = (P_deep − P_shallow) ÷ Interval(m)Gas Grad (bar/m) = (P_deep − P_shallow) ÷ Interval(m)
—

15. Gas / Fluid Gradient Calculator

G(psi/ft) = SG × 0.4335G(kPa/m) = SG × 9.81G(bar/m) = SG × 0.0981
—

16. Influx Type Identification

d_inf(ppg) = MW − (SICP−SIDPP) ÷ (0.052 × h_i)d_inf(kg/m³) = ρ − (SICP−SIDPP) ÷ (0.00981 × h_i)d_inf(kg/l) = ρ − (SICP−SIDPP) ÷ (0.0981 × h_i)
—

17. Gas Volume — Boyle-Charles Combined

V₂ = (P₁ × V₁ × T₂) ÷ (P₂ × T₁) · T in absolute (°R or K)
—

18. Pressure Gradient from Any Fluid Density

G(psi/ft) = ρ(ppg) × 0.052G(kPa/m) = ρ(kg/m³) × 0.00981G(bar/m) = ρ(kg/l) × 0.0981
—

19. Hydrostatic Pressure Difference

ΔHP(psi) = (MW₁ − MW₂) × 0.052 × TVDΔHP(kPa) = (ρ₁ − ρ₂) × 0.00981 × TVDΔHP(bar) = (ρ₁ − ρ₂) × 0.0981 × TVD
—

20. Mud Density Unit Converter

ppg × 119.83 = kg/m³ · ppg ÷ 8.33 = SG · ppg × 7.48 = pcf
—

21. Pressure Unit Converter

psi × 6.895 = kPa · psi × 0.06895 = bar · psi × 0.0703 = kg/cm²
—

22. MAASP — Max Allowable Annular Surface Pressure

MAASP(psi) = (FG − MW) × 0.052 × TVD_shoeMAASP(kPa) = (FG − ρ) × 0.00981 × TVD_shoeMAASP(bar) = (FG − ρ) × 0.0981 × TVD_shoe
—

MAASP is the maximum pressure allowed on the annulus at surface before fracturing the formation at the shoe

Well Control & Kick Management 25 Formulas

1. Kill Mud Weight (KMW)

KMW(ppg) = SIDPP(psi) ÷ (0.052 × TVD) + OMWKMW(kg/m³) = SIDPP(kPa) ÷ (0.00981 × TVD) + OMWKMW(kg/l) = SIDPP(bar) ÷ (0.0981 × TVD) + OMW
—

Metric alt: KMW(kg/m³) = (SIDPP(psi) × 10 ÷ TVD(m)) + OMW(kg/l×1000)

2. Initial Circulating Pressure (ICP)

ICP = SIDPP + SCRP (Slow Circulating Rate Pressure)
—

SCRP must be recorded at the kill rate (typically 25–50 spm) BEFORE the kick occurs

3. Final Circulating Pressure (FCP)

FCP = SCRP × (KMW ÷ OMW)
—

4. Estimated Influx Density

d_inf(ppg) = MW − (SICP − SIDPP) ÷ (0.052 × h_i)d_inf(kg/m³) = ρ − (SICP − SIDPP) ÷ (0.00981 × h_i)d_inf(kg/l) = ρ − (SICP − SIDPP) ÷ (0.0981 × h_i)
—

Gas: <3 ppg · Oil: 3–7 ppg · Saltwater: 8.3–9.0 ppg

5. Kick Tolerance

KT(ppg) = (MAASP ÷ 0.052 ÷ TVD) − (MW × h_i ÷ TVD)KT = [MAASP(kPa) ÷ (0.00981 × TVD)] − [ρ × h_i ÷ TVD]KT = [MAASP(bar) ÷ (0.0981 × TVD)] − [ρ × h_i ÷ TVD]
—

6. Influx Height

h_i(ft) = Pit Gain(bbl) ÷ Ann Capacity(bbl/ft)h_i(m) = Gain(m³) ÷ Ann Cap(m³/m)h_i(m) = Gain(m³) ÷ Ann Cap(m³/m)
—

7. Boyle's Law — Gas Expansion

P₁ × V₁ = P₂ × V₂ → V₂ = (P₁ × V₁) ÷ P₂
—

Also: P₂ = P₁V₁÷V₂. Critical for gas migration and kick handling

8. Gas Percolation / Migration Rate

Mig Rate(ft/hr) = ΔP(psi/hr) ÷ Mud Gradient(psi/ft)Mig Rate(m/hr) = ΔP(kPa/hr) ÷ Mud Grad(kPa/m)Mig Rate(m/hr) = ΔP(bar/hr) ÷ Mud Grad(bar/m)
—

Typical gas migration: 100–1000 ft/hr. Monitor SICP rise rate.

9. U-Tube Balance Equation

SIDPP + Ph_mud(DS) = SICP + Ph_mud(ann) + Ph_gas(ann)
—

10. Driller's Method — ICP & FCP

ICP = SIDPP + SCRP · FCP = SCRP × (KMW÷OMW) · Max BHP = const
—

11. Tapered String — Interim Circulating Pressure

P_interim = ICP − (Strks_pumped ÷ Strks_DS) × (ICP − FCP)
—

12. Strokes to Pump & Spot a Pill

Strokes_pump = Vol_pill(bbl) ÷ PO(bbl/stk) · Strokes_chase = chase_vol ÷ PO
—

13. Volumetric Well Control Method

ΔV(bbl) = ΔP(psi) ÷ (MW×0.052) · Allow P rise, bleed ΔV, maintain BHP
—

14. Lubricate & Bleed — Volume per Step

Vol_lub(bbl) = ΔP_target(psi) ÷ (MW×0.052) × Ann_cap(bbl/ft)
—

15. Riser Margin (Offshore)

RM(ppg) = (ρ_sea×0.052×D_sea − ρ_mud×0.052×(D_sea−TVD)) ÷ (0.052×TVD)RM(kg/m³) = (ρ_sw×D_sea×0.00981 − ρ_mud×0.00981×(D_sea−TVD)) ÷ (0.00981×TVD)RM(kg/l) = same with 0.0981
—

16. Max Pressure at Shoe (Driller's — Gas Kick)

P_shoe_max = Ph_mud_above_shoe + Ph_gas + Csg_pressure
—

17. Choke Line Friction Loss (CLFL) — Subsea

CLFL(psi) ≈ K × (ρ × Q²) — use pressure gauge readingsCLFL(kPa) ≈ measured at kill rate on choke lineCLFL(bar) ≈ measured at kill rate
—

CLFL = SPR × (Kill_Rate/Full_Rate)² where SPR = slow pump rate pressure

18. BOP Accumulator Sizing (API 16D)

Usable Vol(gal) = Close_Vol × P_precharge ÷ (P_min − P_precharge)Usable Vol(L) = Close_Vol × P_pre ÷ (P_min − P_pre)Usable Vol(L) = Close_Vol × P_pre ÷ (P_min − P_pre)
—

API 16D: Total accumulator = Usable Vol ÷ 0.5 (50% usable factor). Min precharge = 200 psi above seabed hydrostatic.

19. Formation Pressure from DST / Wireline

Pf(psi) = FBHP or = MW×0.052×TVD + SIDPPPf(kPa) = ρ×0.00981×TVD + SIDPPPf(bar) = ρ×0.0981×TVD + SIDPP
—

20. Additional Mud Weight for Kill (Over-balance)

ΔMW(ppg) = SIDPP ÷ (0.052 × TVD)Δρ(kg/m³) = SIDPP(kPa) ÷ (0.00981 × TVD)Δρ(kg/l) = SIDPP(bar) ÷ (0.0981 × TVD)
—

21. A-Factor — Gas Migration in Annulus

A = BHP − [Mud_Grad × (D − X) − P_influx]A = BHP − [Mud_Grad × (D − X) − P_influx] (kPa)A = BHP − [Mud_Grad × (D − X) − P_influx] (bar)
—

22. New MAASP After Mud Weight Increase

MAASP_new = (FG − KMW) × 0.052 × TVD_shoeMAASP_new = (FG − KMW) × 0.00981 × TVD_shoeMAASP_new = (FG − KMW) × 0.0981 × TVD_shoe
—

23. Pressure to Break Circulation (Gel Strength)

P_gel(psi) = (YP × L) ÷ (300 × D) + (PV × V × L) ÷ (90000 × D²)P_gel(kPa) = (τ_gel(Pa) × L) ÷ (300 × D) × 47.88P_gel(bar) = P_gel_psi × 0.06895
—

24. Gas Migration & Shut-In Pressure Rise

Mig_Rate(ft/hr) = ΔP_rate(psi/hr) ÷ MW(ppg) ÷ 0.052Mig_Rate(m/hr) = ΔP_rate(kPa/hr) ÷ (ρ × 0.00981)Mig_Rate(m/hr) = ΔP_rate(bar/hr) ÷ (ρ × 0.0981)
—

Typical gas migration: 100–1500 ft/hr. Monitor SICP trend every 15 min when shut in.

25. Pressure at Shoe — Gas Migrating to Shoe

P_shoe = P_surf + Ph_mud_above_shoe + Ph_gas_at_shoeP_shoe = P_surf + Ph_mud + Ph_gas (kPa)P_shoe = P_surf + Ph_mud + Ph_gas (bar)
—
Pump, Volume & Hydraulics 30 Formulas

1. Triplex Pump Output

PO(bbl/stk) = 0.000243 × Liner_ID²(in) × Stroke_Length(in)PO(L/stk) = 0.000785 × Liner_ID²(cm) × Stroke_Length(cm)PO(L/stk) = 0.000785 × Liner_ID²(cm) × Stroke(cm)
—

2. Duplex Pump Output

PO(bbl/stk) = 0.000162 × (2×Bore² − Rod²) × StrokePO(L/stk) = 0.000524 × (2×Bore² − Rod²) × Stroke (cm)PO(L/stk) same formula
—

3. Flow Rate from Pump SPM

Q(bbl/min) = PO(bbl/stk) × SPMQ(L/min) = PO(L/stk) × SPMQ(L/min) = PO(L/stk) × SPM
—

4. Annular Velocity (AV)

AV(ft/min) = (24.51 × Q_bbl/min) ÷ (Dh² − Dp²) in²AV(m/min) = Q(L/min) ÷ [π(Dh²−Dp²)÷4×1000] (Dh,Dp in mm)AV(m/min) = Q(m³/min) ÷ [π(Dh²−Dp²)÷4] (Dh,Dp in m)
—

Minimum AV for hole cleaning: ~100 ft/min (0.5 m/s) in vertical, higher in deviated wells

5. Hydraulic Horsepower (HHP)

HHP = (P_psi × Q_bbl/min) ÷ 1714HHP = (P_kPa × Q_L/min) ÷ 44,742HHP = (P_bar × Q_L/min) ÷ 447.4
—

6. HHP at Bit (HHPb)

HHPb = (Pb_psi × Q_bbl/min) ÷ 1714HHPb = (Pb_kPa × Q_L/min) ÷ 44,742HHPb = (Pb_bar × Q_L/min) ÷ 447.4
—

Target: ≥2.5 HHP/in² of bit area for PDC; ≥1.0 for rollercone

7. Bit Nozzle Pressure Drop

Pb(psi) = 156.5 × MW(ppg) × Q²(bbl/min) ÷ TFA²(in²)Pb(kPa) = 0.8 × ρ(kg/m³) × v_n²(m/s)Pb(bar) = Pb_kPa ÷ 100
—

Field TFA(in²): count nozzle sizes in 32nds. Three 12/32" nozzles → TFA = 3×(12/32)²×π/4 ≈ 0.332 in²

8. Nozzle Jet Velocity & Impact Force

v_n(ft/s) = Q(bbl/min) × 0.3208 ÷ TFA(in²)v_n(m/s) = Q(m³/s) ÷ TFA(m²)v_n(m/s) = Q(m³/s) ÷ TFA(m²)
—

9. Pipe Internal Capacity

Cap(bbl/ft) = ID²(in) ÷ 1029.4Cap(L/m) = π × ID²(mm) ÷ 4000Cap(m³/m) = π × ID²(m) ÷ 4
—

10. Annular Capacity

Ann Cap(bbl/ft) = (Dh² − Dp²)(in²) ÷ 1029.4Ann Cap(L/m) = π(Dh²−Dp²)(mm²) ÷ 4000Ann Cap(m³/m) = π(Dh²−Dp²)(m²) ÷ 4
—

11. Surface to Bit — Lag Strokes

Lag Strokes = DS_Volume ÷ Pump Output
—

12. Bit to Surface — Return Lag Strokes

Return Strokes = Annulus Volume ÷ Pump Output
—

13. Pump Pressure vs. SPM (Roughneck Formula)

P_new = P_old × (SPM_new ÷ SPM_old)²
—

14. Pump Pressure vs. Mud Weight Change

P_new = P_old × (MW_new ÷ MW_old)
—

15. Rectangular Pit Capacity

Vol(bbl) = L(ft) × W(ft) × D(ft) ÷ 5.615Vol(m³) = L × W × D (m)Vol(m³) = L × W × D (m)
—

16. Cylindrical Tank Capacity

Vol(bbl) = π × r²(ft) × H(ft) ÷ 5.615Vol(m³) = π × r²(m) × H(m)Vol(m³) = π × r²(m) × H(m)
—

17. Annular Pressure Loss (from ECD)

APL(psi) = (ECD − MW) × 0.052 × TVDAPL(kPa) = (ECD − ρ) × 0.00981 × TVDAPL(bar) = (ECD − ρ) × 0.0981 × TVD
—

18. Cuttings Slip Velocity & Transport Ratio

V_slip(ft/min) ≈ 0.45×[(65.4−MW)÷MW]^0.5 × d_cut(in)^0.667 ÷ μ_app^0.333V_slip(m/min) = V_slip(ft/min) × 0.3048same
—

19. Shear Rate in Pipe & Annulus

γ_pipe = 96×V(ft/s) ÷ ID(in) · γ_ann = 144×V(ft/s) ÷ (Dh−Dp)(in)γ_pipe = 8V(m/s) ÷ D(m) · γ_ann = 12V(m/s) ÷ (Dh−Dp)(m)same as SI
—

20. Reynolds Number

Re = 928 × ρ(ppg) × V(ft/s) × D(in) ÷ μ(cP)Re = ρ(kg/m³) × V(m/s) × D(m) ÷ μ(Pa·s)same as SI
—

Laminar: Re < 2100 · Transitional: 2100–3000 · Turbulent: >3000

21. Fanning Friction Factor

f_lam = 16÷Re (laminar) · f_turb = 0.0791÷Re^0.25 (Blasius, turbulent)
—

22. Power Law n and K Indices

n = 3.32 × log(θ600 ÷ θ300) · K = θ300 ÷ (511)^n
—

23. Plastic Viscosity & Yield Point (Bingham Plastic)

PV(cP) = θ600 − θ300 · YP(lb/100ft²) = θ300 − PV
—

24. Trip Volume Monitoring

Expected Fill(bbl) = Pipe Displacement(bbl/ft) × Stands Pulled × ft/standExpected Fill(m³) = Pipe Displacement(m³/m) × m pulledsame
—

25. Laminar Pressure Loss — Pipe (Bingham Plastic)

ΔP(psi) = PV×V×L ÷ (90000×ID²) + YP×L ÷ (225×ID)ΔP(kPa) = μ_p×v×L ÷ (3000×D²) + τ_y×L ÷ (6×D)ΔP(bar) = ΔP_kPa ÷ 100
—

26. Laminar Pressure Loss — Annulus (Bingham)

ΔP(psi) = PV×V×L ÷ (60000×(Dh−Dp)²) + YP×L ÷ (200×(Dh−Dp))ΔP(kPa) = μ_p×v×L ÷ (2000×(Dh−Dp)²) + τ_y×L ÷ (6×(Dh−Dp))ΔP(bar) = ΔP_kPa ÷ 100
—

27. Critical Velocity — Laminar/Turbulent Transition

Vc = [97×PV + √(97²×PV² + 8.2×D²×YP×ρ)] ÷ (6.7×ρ×D)Vc(m/s) = Vc_ft/min × 0.00508same
—

28. System HHP & Parasitic Pressure Loss

HHP_sys = P_surf × Q ÷ 1714 · HHP_parasitic = (P_surf − Pb) × Q ÷ 1714kW_sys = P(kPa) × Q(m³/s) ÷ 1000kW = P(bar) × Q(L/s) ÷ 10
—

Optimal: 65% of system HHP at bit. Parasitic loss should be 35% for efficient drilling.

29. Fluid Velocity from Flow Rate & Area

V(ft/s) = Q(bbl/min) × 0.3208 ÷ A(in²)V(m/s) = Q(m³/s) ÷ A(m²)V(m/s) = Q(L/s) ÷ (A(cm²) × 0.0001)
—

30. Minimum Flow Rate for Hole Cleaning

Q_min(bbl/min) = AV_min(ft/min) × (Dh²−Dp²) ÷ 24.51Q_min(L/min) = AV_min(m/min) × π(Dh²−Dp²)÷4000 (mm)Q_min(m³/min) = AV_min(m/min) × π(Dh²−Dp²)÷4 (m)
—

Vertical: 100–120 ft/min min AV · Deviated 30–60°: 150–180 ft/min · Horizontal: 200+ ft/min

Mud & Drilling Fluids Engineering 28 Formulas

1. Barite Required to Raise Mud Weight

Barite(lb/bbl) = 1470 × (MW_new − MW_old) ÷ (35 − MW_new)Barite(kg/m³) = 4200 × (ρ_new − ρ_old) ÷ (4200 − ρ_new)Barite(kg/m³) = 4200 × (ρ_new − ρ_old) ÷ (4200 − ρ_new)
—

Barite SG = 4.2 (35 ppg). Formula valid up to ~19 ppg

2. Mixing Two Fluids — Volume Balance

(V₁ × ρ₁) + (V₂ × ρ₂) = V_F × ρ_F → V₂ = (V_F × ρ_F − V₁ × ρ₁) ÷ ρ₂
—

3. Drill Pipe Slug Design

Dry Pipe (ft) = Slug_Vol × (ρ_slug − ρ_mud) × 1029.4 ÷ DP_ID²Dry Pipe (m) = Slug_Vol × (ρ_slug − ρ_mud) × 1000 ÷ Cap_m³/msame
—

4. Plastic Viscosity & Yield Point

PV(cP) = θ600 − θ300 · YP(lb/100ft²) = θ300 − PV · AV(cP) = PV + YP
—

5. Average Specific Gravity of Mud

ASG = (%vol_water × SG_water + %vol_oil × SG_oil + %vol_solids × SG_solids) ÷ 100
—

6. Low Gravity Solids (LGS) Content

LGS(%vol) = [%vol_solids × (4.2 − ASG)] ÷ (4.2 − 2.6)
—

LGS (SG 2.6) includes bentonite + drilled solids. HGS = Barite (SG 4.2)

7. Oil/Water Ratio Adjustment

Oil volume to add: V_oil = V_total × (%oil_target − %oil_current) ÷ 100
—

8. Bentonite Content Determination

Bent(%vol) = (Mud_CEC − Shale_CEC × %drill_sol) ÷ Bent_CEC
—

9. Complete Mud Solids Analysis

Drilled Solids = %LGS − %Bentonite · Barite = %Suspended − %LGS
—

10. Saltwater Content from Chlorides

SW(%vol) ≈ Cl_ppm ÷ (NaCl_sat_ppm × 0.01) — approx for NaCl
—

Sat. NaCl ≈ 315,000 ppm. 3.5% seawater ≈ 19,000 ppm Cl

11. Dilution Rate to Remove Drilled Solids

V_dilute = V_pit × (%LGS_curr − %LGS_target) ÷ (%LGS_target − %LGS_diluent)
—

12. Centrifuge Evaluation — Underflow Density

ρ_underflow ≈ (Mass_in − Mass_overflow) ÷ Volume_underflow
—

13. Barite for Weighted Cement Slurry

Barite(lb/sk) = [Target_density(ppg) − Base_density] × 94 ÷ (35 − Target)
—

14. Cement Slurry Weight Check

Slurry ρ(ppg) = Total Weight ÷ Total Volume (gal)
—

15. Shear Stress (Power Law & Bingham)

Power Law: τ = K × γⁿ · Bingham: τ = YP + PV×γ÷478.8
—

16. Apparent Viscosity at Shear Rate

AV(cP) = τ(lb/100ft²) × 478.8 ÷ γ(s⁻¹) · Or: AV = PV + YP (at 511 s⁻¹)
—

17. API Filter Loss & Spurt Loss

V_filtrate = 2 × C_w × √t · Spurt Loss = V_30min − 2×(V_7.5min − V_1min)
—

API spec: fresh water mud <15 mL/30min · OBM <2 mL/30min

18. Mud pH & Alkalinity Analysis

Excess lime (lb/bbl) = 0.26 × (Pm − 0.5×Pf) · Carbonate = 0.026×(2×Pf − Pm) · Bicarbonate = 0.012×(Mf − 10×Pm)
—

19. Methylene Blue Test (MBT) — Bentonite

Bentonite(lb/bbl) = MBT_reading × 5 · Bentonite(%vol) = lb/bbl ÷ 9.1
—

API grade bentonite: 22–32 lb/bbl typical · High yield: 10–15 lb/bbl · Native clay: 3–5 lb/bbl

20. Barite Volume % from Solids Analysis

%Barite = %Suspended_Solids − %LGS · Suspended = 100 − %oil − %water
—

21. Chloride Analysis & Salt Content

NaCl(ppm) = Cl⁻(ppm) × 1.65 · CaCl₂(ppm) = Cl⁻(ppm) × 1.56 · KCl(ppm) = Cl⁻(ppm) × 1.91
—

22. Lime Treatment for CO₂ / H₂S

Lime(lb/bbl) = CO₂(mg/L) ÷ 1000 × 0.74 · H₂S(lb/bbl) = H₂S(ppm) × Vol × 0.000034
—

23. Soda Ash for Ca²⁺ Treatment

Soda Ash(lb/bbl) = Ca²⁺(ppm) × 0.00094 · Total(lbs) = lb/bbl × Volume(bbl)
—

Target Ca²⁺: <200 ppm for WBM · High Ca²⁺ causes flocculation and viscosity increase

24. KCl Concentration for Shale Inhibition

KCl(lb/bbl) = %KCl × 3.5 · KCl(kg/m³) = %KCl × 10 · Activity of water: Aw = 1 − 0.00345×%KCl
—

25. Mass Rate Balance — Solids Removal Pit

Q_B(lb/min) = Q_m×ρ_m − Q_c×ρ_c − Q_w×ρ_w · (Vm=mud in, Vc=centrifuge out, Vw=water in)
—

26. Additive Concentration per Cement Sack

%Additive = W_add(lb) ÷ 94 × 100 · lb/sk = %add × 0.94 · gal/sk = V_add(gal)/sk
—

27. Turbulent Pressure Loss (Pipe)

ΔP(psi) = 7.7×10⁻⁵ × f × ρ⁰·⁸ × PV⁰·² × Q¹·⁸ × L ÷ D⁴·⁸ΔP(kPa) = ΔP_psi × 6.895ΔP(bar) = ΔP_psi × 0.06895
—

28. LCM Pill & Loss Circulation Treatment

Pill Vol(bbl) = Ann_Cap(bbl/ft) × Pill_Length(ft) · Barite to weight pill = standard formula
—
Mechanical & Engineering 22 Formulas

1. Buoyancy Factor

BF = 1 − (MW_ppg ÷ 65.44) · W_buoy = W_air × BFBF = 1 − (ρ_kg/m³ ÷ 7850) · (Steel density = 7850 kg/m³)BF = 1 − (ρ_kg/l ÷ 7.85)
—

2. Ton-Mile — Round Trip

TM = D×(L+D)×W + D×(M+0.5×C) ÷ 5,280,000TM = D(m)×(L+D)×W(kN) ÷ 1,000,000same as SI
—

Replace line at: 1000–2000 TM (new bit), 2000–4000 TM (no bit), 500 TM (coring)

3. Free Point — Stuck Pipe

FPD(ft) = (W_pipe_lb/ft × ΔL_stretch_in × 2,500,000) ÷ (Pull_lbs × 12)FPD(m) = (W_pipe × ΔL × E × A) ÷ Pull_forcesame formula
—

Free Point Constant = pipe weight (lb/ft). Stretch formula: FPD = (FPC × in_stretch × 2,500,000) ÷ Pull_lbs

4. D-Exponent (Drilling Rate Analysis)

d = log(R ÷ 60N) ÷ log(12W ÷ 10⁶D)
—

Normalized: dc = d × (ρ_normal ÷ ρ_actual). Decreasing dc_exp → abnormal pore pressure

5. Hookload & Overpull

HL_fluid = W_air × BF · Overpull = HL_measured − HL_expectedHL(kN) = W_air(kN) × BFHL(tonnes) = W_air(tonnes) × BF
—

6. Casing Stretch / Length Change

ΔL_total = ΔL_piston + ΔL_balloon + ΔL_buckle + ΔL_temp
—

E_steel = 30×10⁶ psi (207 GPa) · β_steel = 6.9×10⁻⁶ /°F (12.4×10⁻⁶ /°C)

7. Cement Volume — Annular Fill

Ann Fill(bbl/ft) = (Dh² − Csg_OD²) ÷ 1029.4 · Sacks = Vol ÷ YieldAnn Fill(m³/m) = π(Dh² − OD²) ÷ 4 (Dh,OD in m)same
—

8. Balanced Cement Plug

Plug Length(ft) = Plug_Vol(bbl) ÷ [(DS_Cap + Ann_Cap)(bbl/ft)]Plug Length(m) = Plug_Vol(m³) ÷ [(DS_Cap + Ann_Cap)(m³/m)]same
—

9. Differential Hydrostatic — Cementing

ΔP(psi) = (ρ_cement − ρ_mud) × 0.052 × TVD · Check: must not exceed MAASPΔP(kPa) = (ρ_cem − ρ_mud) × 0.00981 × TVDΔP(bar) = (ρ_cem − ρ_mud) × 0.0981 × TVD
—

10. Drill Collar Length for WOB

DC_Length(ft) = WOB(lbs) ÷ (BF × DC_Weight_lb/ft)DC_Length(m) = WOB(kN) ÷ (BF × DC_Weight_kN/m)same
—

Add 15–25% safety factor — run only 85% of DC in compression

11. Matthews & Kelly Fracture Gradient

FG(psi/ft) = Pore_Pressure(psi/ft) + Ki × σ_matrix(psi/ft) ÷ TVDFG(kPa/m) = Pp(kPa/m) + Ki × σ(kPa/m) ÷ TVDsame ratio
—

Ki ≈ 0.3 shallow · 0.5 at 5000ft · 0.8–0.9 deep. Check against LOT/FIT data.

12. Piston Effect — Casing Length Change ΔL₁

ΔL₁(in) = ΔFa(lb) × L(ft) × 12 ÷ (A(in²) × 30×10⁶)ΔL₁(mm) = ΔFa(N) × L(m) × 1000 ÷ (A(m²) × 207×10⁹)same as SI
—

E_steel = 30×10⁶ psi (207 GPa). Positive = elongation, negative = shortening.

13. Ballooning Effect — Length Change ΔL₂

ΔL₂(in) = −2ν × (ΔP_i×r_i² − ΔP_o×r_o²) ÷ (r_o²−r_i²) × L×12 ÷ EΔL₂(mm) = −2ν × ΔP_avg × L × 1000 ÷ E(Pa)same as SI
—

ν (Poisson's ratio for steel) = 0.3. Negative sign = shortening when pressure increases inside.

14. Temperature Effect — Length Change ΔL₄

ΔL₄(in) = β × L(ft) × ΔT(°F) × 12 · β_steel = 6.9×10⁻⁶ /°FΔL₄(mm) = β × L(m) × ΔT(°C) × 1000 · β = 12.4×10⁻⁶ /°Csame as SI
—

15. Buckling Effect — Length Change ΔL₃

ΔL₃(in) = −(r_n)² × ΔP_avg × L(ft) × 12 ÷ (8 × E × I / ID²) (simplified)ΔL₃(mm) ≈ −Fa_buckling × L ÷ (2 × E × A)same
—

16. Total Casing Length Change (All Effects)

ΔL_total = ΔL₁(piston) + ΔL₂(balloon) + ΔL₃(buckle) + ΔL₄(temp)
—

17. Casing Yield Strength at Temperature

Yield_T = Yield_20°C × [1 − 0.02 × (T − 200) ÷ 50] for T > 200°F
—

Reduce ~2% per 50°F above 200°F (93°C). Applies to L-80, P-110, Q-125 grades.

18. Mix Water for Cement — Volume & WCR

Water(gal/sk) = WCR × 94 ÷ 8.33 · WCR = water weight ÷ cement weight
—

Class A: WCR=0.46 (5.2 gal/sk) · Class G: WCR=0.44 (5.0 gal/sk) · Class H: WCR=0.38 (4.3 gal/sk)

19. Cement Slurry Yield (ft³/sk)

Yield(ft³/sk) = (94÷Density_cement + Water(gal/sk)÷7.481 + Additive_vol) ÷ 1 · 1 ft³ = 7.481 gal
—

Class G typical yield: 1.15 ft³/sk (5.0 gal/sk water) = 1.06 bbl/sk

20. Drill String Torque & Drag

Drag(lbs) = μ × N × ΔL · Torque(ft-lbs) = μ × N × r_string · N = Normal force on formationDrag(kN) = μ × N × ΔLsame
—

21. Neutral Point — BHA in Compression

NP_depth(ft) = TD − (WOB ÷ BF ÷ DC_wt_ppf)NP_depth(m) = TD − (WOB_kN ÷ (BF × DC_wt_kNm))same
—

Keep neutral point inside drill collars. DP should always be in tension.

22. Critical Buckling Force — Sinusoidal

F_crit(lbs) = 2 × √(E×I×W_buoy÷r_c) · r_c = (Dh−DP_OD)÷2F_crit(N) = 2 × √(E×I×W_buoy÷r_c)same
—

If WOB > F_crit, sinusoidal buckling occurs. If WOB > 2×F_crit, helical buckling.

Directional Drilling & Survey Calculations 13 Formulas

1. Dogleg Severity (DLS)

DLS = (100÷MD_interval) × arccos[cos(I₂−I₁) − sin(I₁)×sin(I₂)×(1−cos(Az₂−Az₁))]
—

Max DLS: PDC 3°/100ft · Motor 8°/100ft · Rotary BHA 1.5°/100ft

2. Minimum Radius of Curvature

R(ft) = 5730 ÷ DLS(°/100ft)R(m) = 1746 ÷ DLS(°/30m)R(m) = 1746 ÷ DLS(°/30m)
—

3. Minimum Curvature Method (North/East/Vertical)

ΔNorth = (MD/2) × (sinI₁cosAz₁ + sinI₂cosAz₂) × RF
ΔEast = (MD/2) × (sinI₁sinAz₁ + sinI₂sinAz₂) × RF
ΔVert = (MD/2) × (cosI₁ + cosI₂) × RF
—

4. Average Angle Method Survey

ΔN = MD × sin(AvgI) × cos(AvgAz) · ΔE = MD × sin(AvgI) × sin(AvgAz) · ΔV = MD × cos(AvgI)
—

5. Departure & Horizontal Distance

Departure(CD) = sin(I) × Build_Distance · Horizontal = ΔN² + ΔE²
—

6. Build Rate to Hit Target TVD/Departure

BUR(°/100ft) = [arctan(Dep÷TVD_interval)] × (100 ÷ MD_build)
—

7. Well Planning — KOP, Build, Tangent Distances

R = 5730÷BUR · TVD_build = R×sin(I) · MD_build = R×I(rad) · Tangent_TVD = (TVD_target − TVD_kop − TVD_build) ÷ cos(I)
—

8. Radius of Curvature Survey Method

ΔN = [300÷(I₂−I₁)] × [sin(I₂)−sin(I₁)] × [sin(Az₂)−sin(Az₁)] ÷ (Az₂−Az₁)
—

9. Effective WOB in Deviated Well

WOB_eff = WOB_surface × cos(I) − FF × W_BHA × sin(I)
—

10. Casing Setting Depth Design

D_shoe(ft) = (P_kick_psi) ÷ (FG_shoe × 0.052) · Check: FG_shoe(ppg) > KMW_nextD_shoe(m) = P_kick(kPa) ÷ (FG × 0.00981)D_shoe(m) = P_kick(bar) ÷ (FG × 0.0981)
—

Shoe must be deep enough that FG_shoe exceeds KMW of the next section plus safety margin.

11. Riser Margin — Deep Water

RM(ppg) = (ρ_sw × D_water − ρ_mud × (D_water − TVD)) ÷ TVDRM(kg/m³) = (ρ_sw×D_w − ρ_mud×(D_w−TVD)) ÷ TVDsame ratio
—

If riser is disconnected, mud weight must be increased by riser margin to maintain wellbore stability.

12. Tool Face Angle — Motor Orientation

Gravity TF: clockwise from high side · Magnetic TF: clockwise from magnetic north (0–360°)
—

Gravity TF used when inclination >5°. Magnetic TF used when inclination <5° (near vertical).

13. Closure Distance & Azimuth

CD = √(ΔN² + ΔE²) · CL_Az(°) = atan2(ΔE, ΔN) × 180/π (corrected to 0–360°)
—

Bonus: Wellbore Tortuosity Index

Tortuosity = Σ(DLS_i × MD_interval_i) ÷ Total_MD · Units: °/100ft
—

Low tortuosity: <1°/100ft · Moderate: 1–3°/100ft · High (affects T&D): >3°/100ft

Additional Calculations 40 Formulas · follows the Field / SI / Metric switch

Merged from the Lapeyrouse and Pressure calculators. Inputs and results use the selected unit system; the formulas themselves are shown in field units.

Pressure & Mud Weight

Required Mud Weight from Pressure

MW (ppg) = Pressure (psi) ÷ (0.052 × TVD ft)
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Lost Returns Pressure Change

ΔP (psi) = (MW_new − MW_old) × 0.052 × TVD
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Pore Pressure & Fracture Gradient

Overburden Gradient from Bulk Density

OBG (psi/ft) = 0.433 × ρ_bulk (g/cc); ppg = 8.33 × ρ
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Pore Pressure from d-exponent (Ratio Method)

PP = MW_normal × (dc_normal ÷ dc_observed)
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Pore Pressure from d-exponent (Eaton)

PP = OBG − (OBG − PP_n) × (dc_observed ÷ dc_normal)^1.2
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Pore Pressure from Resistivity (Eaton)

PP = OBG − (OBG − PP_n) × (R_observed ÷ R_normal)^1.2
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Effective Stress

σ_eff = Total stress (OBP) − Pore pressure
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Fracture Gradient (Eaton / Poisson)

FG = [ν ÷ (1 − ν)] × (OBG − PP) + PP
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Fracture Pressure at Depth

P_frac (psi) = FG (ppg) × 0.052 × TVD
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Mud Weight Window

Window = FG − PP (safe MW between PP and FG)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Overbalance & Margin Check

Overbalance = MW − PP; Margin to fracture = FG − MW
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Poisson's Ratio from Depth (Empirical)

ν ≈ min(0.125 + TVD ÷ 100,000 , 0.48) (Gulf Coast)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Well Control

Choke Pressure During Kill

P_choke = ICP − (strokes pumped ÷ strokes to bit) × (ICP − FCP)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Kick Tolerance (Subsea BOP)

KT = FP_shoe − MW×0.052×Shoe − Riser margin loss
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Hydraulics

HSI — Hydraulic Horsepower per Area of Bit

HSI = HHPb ÷ (π/4 × D_bit²)
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Standpipe Pressure Budget

P_s = P_surface + P_drillstring + P_bit + P_annulus
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Optimum Bit Pressure Drop

ΔP_bit,opt = % × P_standpipe (≈65% for soft formations)
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Annular Pressure Loss (Simplified)

APL ≈ MW × AV × L ÷ (300 × (D_h − D_p)²)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Surge & Swab Pressure (Approx.)

ΔP ≈ 0.4 × MW × V_trip² ÷ 25.8
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Volumes & Displacement

Displacement Volume — Pipe Steel

V (bbl) = (OD² − ID²) ÷ 1029.4 × Length
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Strokes to Displace / Fill

Strokes = Volume ÷ Pump output
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Dilution to Reduce Mud Weight

V_water = V × (MW_old − MW_new) ÷ (MW_new − 8.33)
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Cementing

Sacks of Cement Required

Sacks = V_bbl × 5.615 ÷ Yield × Excess factor
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Displacement Volume (Cementing)

V (bbl) = ID² ÷ 1029.4 × Depth
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Waiting on Cement (WOC) Estimate

Rule-of-thumb WOC time from temperature at depth
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Casing Design

Burst Pressure Rating (Barlow)

P_burst = 0.875 × 2 × Ym × t ÷ OD
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Burst Pressure on Casing (Net Load)

P_burst load = Formation/internal pressure − outside mud hydrostatic
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Collapse Pressure (Yield Strength)

P_col = Ym × 2(D/t − 1) ÷ (D/t)²
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Collapse Pressure on Casing (Net Load)

P_col load = (ρ_out − ρ_in) × 0.052 × Depth
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Axial Tension Capacity (Yield)

F_y = Ym × 0.7854 × (OD² − ID²)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Casing Design Safety Factors

SF = Rating ÷ Load (min: burst 1.1 · collapse 1.0 · tension 1.6)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Triaxial (Von Mises) Stress

σ_VME = √(σ_a² − σ_a·σ_h + σ_h²)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Drill String & Pipe

Drill Pipe Tensile Safety Factor

SF = Pipe yield strength ÷ Maximum hookload
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Maximum Allowable Overpull

Overpull_max = Yield ÷ SF − Current hookload
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Maximum Pipe Length (Tension Limit)

L_max = (Ym × As − W_other) ÷ (w × BF)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Drill String Weight in Mud (Hookload)

W_buoyed = (w_dp·L_dp + w_dc·L_dc) × BF, BF = (65.5 − MW) ÷ 65.5
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Drill String Stretch

ΔL = W × L ÷ (As × E), E = 30×10⁶ psi
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Bending Stress in Pipe

σ_b = E × OD × α ÷ (2 × 12), α = DLS/100 in rad/ft (≈ 218 × OD × DLS)
—

Follows the unit switch · formula shown in field units · from Pressure Calc

Ton-Miles — Setting Casing

TM = W_c × D × BF × 2 ÷ (2000 × 5280)
—

Follows the unit switch · formula shown in field units · from Lapeyrouse

Directional

Turn Rate (Azimuth Change)

Turn rate = |Az₂ − Az₁| ÷ ΔMD × 100
—

Follows the unit switch · formula shown in field units · from Pressure Calc