1. Pressure Gradient
psi/ft = MW(ppg) × 0.052kPa/m = ρ(kg/m³) ÷ 102kg/cm²/m = ρ(kg/l) ÷ 10Field: 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_shoe4. 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 + SIDPPVerify: 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 + CLFL9. 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 + 27312. Crude Oil Density from °API
SG = 141.5 ÷ (131.5 + °API) · Density(ppg) = SG × 8.3313. 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.098116. 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.098119. Hydrostatic Pressure Difference
ΔHP(psi) = (MW₁ − MW₂) × 0.052 × TVDΔHP(kPa) = (ρ₁ − ρ₂) × 0.00981 × TVDΔHP(bar) = (ρ₁ − ρ₂) × 0.0981 × TVD20. Mud Density Unit Converter
ppg × 119.83 = kg/m³ · ppg ÷ 8.33 = SG · ppg × 7.48 = pcf21. 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_shoeMAASP is the maximum pressure allowed on the annulus at surface before fracturing the formation at the shoe
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) + OMWMetric 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 = const11. 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 ÷ PO13. Volumetric Well Control Method
ΔV(bbl) = ΔP(psi) ÷ (MW×0.052) · Allow P rise, bleed ΔV, maintain BHP14. 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.098116. Max Pressure at Shoe (Driller's — Gas Kick)
P_shoe_max = Ph_mud_above_shoe + Ph_gas + Csg_pressure17. 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 rateCLFL = 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 + SIDPP20. 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_shoe23. 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.0689524. 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)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 formula3. Flow Rate from Pump SPM
Q(bbl/min) = PO(bbl/stk) × SPMQ(L/min) = PO(L/stk) × SPMQ(L/min) = PO(L/stk) × SPM4. 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.46. HHP at Bit (HHPb)
HHPb = (Pb_psi × Q_bbl/min) ÷ 1714HHPb = (Pb_kPa × Q_L/min) ÷ 44,742HHPb = (Pb_bar × Q_L/min) ÷ 447.4Target: ≥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 ÷ 100Field 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) ÷ 410. Annular Capacity
Ann Cap(bbl/ft) = (Dh² − Dp²)(in²) ÷ 1029.4Ann Cap(L/m) = π(Dh²−Dp²)(mm²) ÷ 4000Ann Cap(m³/m) = π(Dh²−Dp²)(m²) ÷ 411. Surface to Bit — Lag Strokes
Lag Strokes = DS_Volume ÷ Pump Output12. Bit to Surface — Return Lag Strokes
Return Strokes = Annulus Volume ÷ Pump Output13. 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 × TVD18. 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.3048same19. 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 SI20. Reynolds Number
Re = 928 × ρ(ppg) × V(ft/s) × D(in) ÷ μ(cP)Re = ρ(kg/m³) × V(m/s) × D(m) ÷ μ(Pa·s)same as SILaminar: 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)^n23. Plastic Viscosity & Yield Point (Bingham Plastic)
PV(cP) = θ600 − θ300 · YP(lb/100ft²) = θ300 − PV24. Trip Volume Monitoring
Expected Fill(bbl) = Pipe Displacement(bbl/ft) × Stands Pulled × ft/standExpected Fill(m³) = Pipe Displacement(m³/m) × m pulledsame25. 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 ÷ 10026. 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 ÷ 10027. Critical Velocity — Laminar/Turbulent Transition
Vc = [97×PV + √(97²×PV² + 8.2×D²×YP×ρ)] ÷ (6.7×ρ×D)Vc(m/s) = Vc_ft/min × 0.00508same28. 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) ÷ 10Optimal: 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
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³/msame4. Plastic Viscosity & Yield Point
PV(cP) = θ600 − θ300 · YP(lb/100ft²) = θ300 − PV · AV(cP) = PV + YP5. Average Specific Gravity of Mud
ASG = (%vol_water × SG_water + %vol_oil × SG_oil + %vol_solids × SG_solids) ÷ 1006. 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) ÷ 1008. Bentonite Content Determination
Bent(%vol) = (Mud_CEC − Shale_CEC × %drill_sol) ÷ Bent_CEC9. Complete Mud Solids Analysis
Drilled Solids = %LGS − %Bentonite · Barite = %Suspended − %LGS10. Saltwater Content from Chlorides
SW(%vol) ≈ Cl_ppm ÷ (NaCl_sat_ppm × 0.01) — approx for NaClSat. 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_underflow13. 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.816. 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.1API 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 − %water21. Chloride Analysis & Salt Content
NaCl(ppm) = Cl⁻(ppm) × 1.65 · CaCl₂(ppm) = Cl⁻(ppm) × 1.56 · KCl(ppm) = Cl⁻(ppm) × 1.9122. Lime Treatment for CO₂ / H₂S
Lime(lb/bbl) = CO₂(mg/L) ÷ 1000 × 0.74 · H₂S(lb/bbl) = H₂S(ppm) × Vol × 0.00003423. 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×%KCl25. 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)/sk27. 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.0689528. LCM Pill & Loss Circulation Treatment
Pill Vol(bbl) = Ann_Cap(bbl/ft) × Pill_Length(ft) · Barite to weight pill = standard formula1. 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 SIReplace 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 formulaFree 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) × BF6. Casing Stretch / Length Change
ΔL_total = ΔL_piston + ΔL_balloon + ΔL_buckle + ΔL_tempE_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)same8. 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)]same9. 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 × TVD10. 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)sameAdd 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 ratioKi ≈ 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 SIE_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 SI15. 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)same16. 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°FReduce ~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 weightClass 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 galClass 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 × ΔLsame21. Neutral Point — BHA in Compression
NP_depth(ft) = TD − (WOB ÷ BF ÷ DC_wt_ppf)NP_depth(m) = TD − (WOB_kN ÷ (BF × DC_wt_kNm))sameKeep 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)sameIf WOB > F_crit, sinusoidal buckling occurs. If WOB > 2×F_crit, helical buckling.
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 ratioIf 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: °/100ftLow tortuosity: <1°/100ft · Moderate: 1–3°/100ft · High (affects T&D): >3°/100ft
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 × TVDFollows 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.2Follows 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.2Follows the unit switch · formula shown in field units · from Pressure Calc
Effective Stress
σ_eff = Total stress (OBP) − Pore pressureFollows the unit switch · formula shown in field units · from Pressure Calc
Fracture Gradient (Eaton / Poisson)
FG = [ν ÷ (1 − ν)] × (OBG − PP) + PPFollows the unit switch · formula shown in field units · from Lapeyrouse
Fracture Pressure at Depth
P_frac (psi) = FG (ppg) × 0.052 × TVDFollows 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 − MWFollows 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 lossFollows 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_annulusFollows 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.8Follows the unit switch · formula shown in field units · from Pressure Calc
Volumes & Displacement
Displacement Volume — Pipe Steel
V (bbl) = (OD² − ID²) ÷ 1029.4 × LengthFollows the unit switch · formula shown in field units · from Lapeyrouse
Strokes to Displace / Fill
Strokes = Volume ÷ Pump outputFollows 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 factorFollows the unit switch · formula shown in field units · from Lapeyrouse
Displacement Volume (Cementing)
V (bbl) = ID² ÷ 1029.4 × DepthFollows the unit switch · formula shown in field units · from Lapeyrouse
Waiting on Cement (WOC) Estimate
Rule-of-thumb WOC time from temperature at depthFollows the unit switch · formula shown in field units · from Lapeyrouse
Casing Design
Burst Pressure Rating (Barlow)
P_burst = 0.875 × 2 × Ym × t ÷ ODFollows 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 hydrostaticFollows 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 × DepthFollows 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 hookloadFollows the unit switch · formula shown in field units · from Lapeyrouse
Maximum Allowable Overpull
Overpull_max = Yield ÷ SF − Current hookloadFollows 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.5Follows the unit switch · formula shown in field units · from Pressure Calc
Drill String Stretch
ΔL = W × L ÷ (As × E), E = 30×10⁶ psiFollows 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 × 100Follows the unit switch · formula shown in field units · from Pressure Calc