Formula sheet

Tailings Water Balance and Seepage Control Formula Sheet

Mining engineering formula sheet for tailings storage facility water balance, pond volume, freeboard, beach geometry, seepage, drain response, water-quality load, monitoring trends, uncertainty, and release checks.

This formula sheet collects first-pass relationships used to review tailings storage facility water balance, pond storage, seepage control, drain response, freeboard, water-quality load, monitoring trends, uncertainty, and operating release evidence.

Use these equations to make tailings and mine-closure calculations traceable. They do not replace site-specific geotechnical design, consequence classification, climate modelling, dam safety review, approved operating manuals, emergency action plans, competent-person signoff, or regulatory requirements.

Before calculating, state the facility state: active deposition, temporary shutdown, storm operation, upset water inventory, construction raise, closure transition, or post-closure monitoring. A water balance that is acceptable during normal operation can be unsafe during an extreme rainfall sequence, blocked reclaim, drain malfunction, or restricted decant condition.

Symbols and Basis

Use SI units unless a site standard requires otherwise. Keep dry solids, wet slurry, pond water, void water, rainfall, seepage, and reclaim streams on separate bases.

SymbolMeaningCommon unit
m_sdry tailings solids mass\text{t} or \text{kg}
\dot{m}_sdry tailings solids rate\text{t/day}
\rho_ddeposited dry density\text{t/m}^3
\rho_wwater density\text{kg/m}^3
w_sslurry solids mass fractiondimensionless
Spond or system water storage\text{m}^3
Qvolumetric flow rate\text{m}^3/\text{s} or \text{m}^3/\text{day}
P_dprecipitation depth\text{m}
E_devaporation depth\text{m}
Aarea normal to flow or water surface area\text{m}^2
Khydraulic conductivity\text{m/s}
ihydraulic gradientdimensionless
F_bfreeboard\text{m}
Cconcentration\text{mg/L}
U_{95}approximate expanded uncertainty at 95 percent confidenceengineering unit

State whether water volumes are measured, modelled, estimated from stage-storage curves, or inferred from pump runtime. Tailings water balances are often wrong because one stream is reported as dry solids, another as wet tonnes, and another as pond level without a consistent reference boundary.

Tailings Solids and Deposited Volume

Dry tailings solids from ore feed and product yield:

m_t=m_{ore}(1-y_p)

where y_p is the product or concentrate mass yield.

Dry tailings rate:

\dot{m}_t=\dot{m}_{ore}(1-y_p)

Deposited solids volume:

\displaystyle V_{solids}=\frac{m_t}{\rho_d}

Average facility fill rate:

\displaystyle \dot{V}_{fill}=\frac{\dot{m}_t}{\rho_d}

The deposited dry density is not a universal material property. It depends on particle size, mineralogy, thickening, deposition method, beaching, consolidation, desiccation, entrained water, and time. Use measured beach and density data when the calculation supports a raise schedule, freeboard decision, or closure storage statement.

Slurry Water and Entrained Water

Slurry mass from dry solids and solids mass fraction:

\displaystyle m_{slurry}=\frac{m_s}{w_s}

Water mass in slurry:

\displaystyle m_w=m_s\left(\frac{1-w_s}{w_s}\right)

Approximate slurry volume:

\displaystyle V_{slurry}=\frac{m_s}{\rho_s}+\frac{m_w}{\rho_w}

where \rho_s is the particle density. If the site reports slurry density \rho_{slurry} directly, use:

\displaystyle V_{slurry}=\frac{m_{slurry}}{\rho_{slurry}}

The engineering check is not only pumpability. The water carried with the solids affects pond rise, reclaim demand, seepage load, deposition beaching, and closure consolidation.

Discrete TSF Water Balance

For a chosen time step:

S_{k+1}=S_k+\left(Q_{proc}+Q_{rain}+Q_{runon}+Q_{gw}-Q_{reclaim}-Q_{evap}-Q_{seep}-Q_{dis}\right)\Delta t

where the sign convention treats inflows to the TSF water inventory as positive. Use the same boundary for every term. For example, if seepage is collected and pumped back into the pond, it may leave one control volume but not the mine-wide water circuit.

Rainfall on pond surface:

V_{rain}=P_d A_{pond}

Runoff from contributing catchment:

V_{runon}=C_r P_d A_c

Evaporation loss:

V_{evap}=E_d A_{pond}

Reclaim volume from pump flow:

V_{reclaim}=Q_{reclaim}\Delta t

Use a stage-storage curve rather than a flat-area approximation when pond geometry changes significantly with elevation. For screening over a small level change:

\displaystyle \Delta z_{pond}\approx\frac{\Delta S}{A_{pond}}

Water-balance validation should reconcile pond survey, flow meters, pump runtime, rainfall gauges, seepage collections, water-quality sampling, and process plant make-up water. A model that closes numerically but disagrees with pond level trends is not release evidence.

Pond Freeboard and Storage Margin

Available freeboard:

F_b=z_{crest}-z_{pond}

Freeboard margin above the required operating minimum:

M_F=F_b-F_{req}

Approximate remaining storage above the pond before reaching the freeboard limit:

V_{margin}\approx A_{pond}M_F

This approximation is only a quick screen. A formal check should use the surveyed stage-storage relationship, wave run-up basis, storm inflow design, spillway or decant capacity, settlement allowance, construction tolerance, beach condition, and consequence classification.

Do not treat positive freeboard as proof that seepage risk is acceptable. Overtopping risk and internal erosion risk can move in different directions.

Beach Slope and Pond Offset

Approximate beach elevation drop over a deposition distance:

\Delta z_b=S_b L_b

where S_b is beach slope as vertical drop per horizontal length and L_b is horizontal distance. Pond offset from the embankment should be checked against the operating plan:

M_{offset}=L_{pond}-L_{min}

where L_{pond} is actual pond setback and L_{min} is the required minimum setback.

Beach slope is operational evidence, not just geometry. It changes with solids concentration, particle size distribution, flow rate, deposition point, spigot operation, thickener performance, wind, ice, and traffic. A flatter-than-assumed beach can move the pond toward the embankment and reduce seepage and freeboard margin.

Seepage by Darcy Flow

Darcy flow through a saturated control zone:

Q=KiA

Hydraulic gradient:

\displaystyle i=\frac{\Delta h}{L}

Hydraulic conductivity inferred from measured seepage:

\displaystyle K=\frac{Q}{iA}

Convert seepage from cubic metres per second to cubic metres per day:

Q_{\text{m}^3/\text{day}}=86400Q_{\text{m}^3/\text{s}}

Darcy flow assumes laminar flow through a representative porous medium. It is not enough for cracked zones, open defects, preferential foundation windows, internal erosion, concentrated pipe flow, or unsaturated transient flow unless those mechanisms have been explicitly modelled.

Exit Gradient and Piping Screen

Simplified exit gradient:

\displaystyle i_e=\frac{h_e}{L_e}

Critical hydraulic gradient for a quick condition screen:

\displaystyle i_{cr}=\frac{G_s-1}{1+e}

Exit-gradient safety factor:

\displaystyle FS_i=\frac{i_{cr}}{i_e}

where G_s is soil solids specific gravity and e is void ratio. The required safety factor depends on consequence, material susceptibility, filter compatibility, monitoring reliability, and design basis.

This screen should be interpreted with seepage clarity, suspended solids, piezometer trend, drain function, filter compatibility, deformation, cracking, pond position, and emergency action triggers. A calculated safety factor slightly above a minimum does not justify continued operation when seepage is turbid or rapidly increasing.

Drain and Collection Response

Drain response ratio:

\displaystyle R_d=\frac{Q_{meas}}{Q_{expected}}

Seepage collection balance:

Q_{seep,total}=Q_{drain}+Q_{toe}+Q_{uncollected}

Pond-source seepage fraction for a simplified water-balance audit:

\displaystyle f_{seep}=\frac{Q_{seep,total}}{Q_{pond,out}+Q_{seep,total}}

Use these as operating indicators, not as proof of hydraulic performance. A low drain response ratio may indicate clogging, measurement error, lower hydraulic gradient, flow bypass, frozen drains, sedimentation, damaged outlets, or a changed pond boundary. The response must be checked against piezometers and field inspection.

Water-Quality Load

Daily constituent load:

\displaystyle L=\frac{CQ}{1000}

where C is in \text{mg/L}, Q is in \text{m}^3/\text{day}, and L is in \text{kg/day}.

Mixing concentration for two water streams:

\displaystyle C_m=\frac{C_1Q_1+C_2Q_2}{Q_1+Q_2}

Removal efficiency:

\displaystyle \eta_R=\frac{C_{in}-C_{out}}{C_{in}}

Water-quality calculations must state whether the sample is filtered, unfiltered, dissolved, total, flow-proportional, grab, composite, field-preserved, or laboratory-corrected. Tailings seepage decisions often depend on turbidity, suspended solids, pH, sulfate, salinity, metals, process reagents, and oxidation products, not on flow alone.

Pumping Power for Reclaim or Pump-Back

Hydraulic power:

P_h=\rho_w gQH

Pump shaft or electrical input power estimate:

\displaystyle P_{in}=\frac{\rho_w gQH}{\eta}

where Q is in \text{m}^3/\text{s}, H is total dynamic head, and \eta is combined efficiency. Pump-back reliability should also check standby pumps, power availability, pipe capacity, winter operation, siltation, controls, instrumentation, and access during storm conditions.

Monitoring Trend and Time to Trigger

Linear rate of change for a monitored head, pond level, or seepage flow:

\displaystyle r=\frac{x_2-x_1}{t_2-t_1}

Estimated time to trigger:

\displaystyle t_{trigger}=\frac{x_{trigger}-x_{now}}{r}

Use this only when the trend is physically credible over the forecast period. Piezometric response can be nonlinear after rainfall, deposition changes, drain blockage, pond relocation, consolidation, thaw, or construction activity.

Uncertainty and Guard Bands

Approximate combined standard uncertainty:

u_c=\sqrt{u_1^2+u_2^2+\cdots+u_n^2}

Approximate expanded uncertainty:

U_{95}\approx2u_c

Guard-banded margin to an upper limit:

M_{GB}=x_{limit}-x_{meas}-U_{95}

For a lower limit, such as minimum freeboard:

M_{GB}=x_{meas}-x_{limit}-U_{95}

If the guard-banded margin is negative, do not treat the result as compliant without additional evidence. Tailings decisions should account for survey error, flow-meter uncertainty, rainfall spatial variation, stage-storage error, seepage measurement error, piezometer datum checks, and sampling variability.

Worked Example 1: Storm Pond Rise and Freeboard

A TSF pond starts at elevation 121.80\ \text{m} with crest elevation 125.00\ \text{m}. Required operating freeboard is 2.00\ \text{m}. During a two-day wet period:

QuantityValue
process water to pond18{,}000\ \text{m}^3/\text{day}
reclaim pumping20{,}000\ \text{m}^3/\text{day}
evaporation plus seepage loss1{,}800\ \text{m}^3/\text{day}
event rainfall depth0.095\ \text{m}
pond area145{,}000\ \text{m}^2
contributing catchment260{,}000\ \text{m}^2
runoff coefficient0.45

Rainfall on the pond:

V_{rain}=0.095(145{,}000)=13{,}775\ \text{m}^3

Runon:

V_{runon}=0.45(0.095)(260{,}000)=11{,}115\ \text{m}^3

Process inflow over two days:

V_{proc}=18{,}000(2)=36{,}000\ \text{m}^3

Reclaim and losses:

V_{out}=20{,}000(2)+1{,}800(2)=43{,}600\ \text{m}^3

Net storage increase:

\Delta S=36{,}000+13{,}775+11{,}115-43{,}600=17{,}290\ \text{m}^3

Approximate pond rise:

\displaystyle \Delta z_{pond}=\frac{17{,}290}{145{,}000}=0.119\ \text{m}

Final pond elevation:

z_{pond,final}=121.80+0.119=121.919\ \text{m}

Final freeboard:

F_b=125.00-121.919=3.081\ \text{m}

Freeboard margin:

M_F=3.081-2.00=1.081\ \text{m}

Engineering Comment

The event does not consume the simplified freeboard requirement, but the pond rises by about 0.12\ \text{m}. The release decision should still check forecast rainfall, stage-storage accuracy, pump availability, decant restrictions, beach condition, pond setback, and whether the seepage response changes after the event. This is an operating screen, not a dam-safety approval.

Worked Example 2: Seepage Flow and Exit-Gradient Screen

A low-permeability control zone has K=3.0\times10^{-7}\ \text{m/s}, average gradient i=0.42, and seepage area A=4200\ \text{m}^2. At the downstream toe, a local piezometric head of 3.1\ \text{m} is estimated over an exit path length of 5.5\ \text{m}. The exit-zone material has G_s=2.68 and e=0.78. Expected drain flow is 110\ \text{m}^3/\text{day}; measured drain flow is 65\ \text{m}^3/\text{day}.

Darcy seepage:

Q=KiA=(3.0\times10^{-7})(0.42)(4200)=5.29\times10^{-4}\ \text{m}^3/\text{s}

Daily flow:

Q=5.29\times10^{-4}(86400)=45.7\ \text{m}^3/\text{day}

Exit gradient:

\displaystyle i_e=\frac{3.1}{5.5}=0.564

Critical gradient:

\displaystyle i_{cr}=\frac{2.68-1}{1+0.78}=0.944

Exit-gradient safety factor:

\displaystyle FS_i=\frac{0.944}{0.564}=1.67

Drain response ratio:

\displaystyle R_d=\frac{65}{110}=0.59

Engineering Comment

The simplified exit-gradient safety factor is above a screening value of 1.5, but the drain response ratio is low. That means the decision should not be “continue as normal.” The engineering response is to inspect the drain outlet, check for blockage or bypass, validate the piezometer datum, compare seepage clarity and suspended solids, and increase monitoring until the hydraulic path is understood. If seepage becomes turbid, flow rises quickly, or piezometers continue to climb, the facility should move to the relevant trigger-action response level even if the simple gradient screen has not failed.

Common Mistakes

  1. Mixing dry tonnes, wet tonnes, slurry volume, and pond water without a declared basis.
  2. Using average annual rainfall to justify short-term storm freeboard.
  3. Treating reclaimed water as removed from risk when pump capacity, power supply, and access are not reliable.
  4. Using Darcy flow for a cracked, eroded, or preferential path without checking whether the porous-medium assumption still applies.
  5. Reviewing freeboard without pond offset, wave run-up, beach slope, stage-storage uncertainty, and decant restrictions.
  6. Treating clear seepage, turbid seepage, drain flow, and piezometer head as independent observations instead of one hydraulic system.
  7. Using a single laboratory hydraulic conductivity value without field calibration or sensitivity checks.
  8. Declaring a facility acceptable from a single calculated safety factor while trigger-action response evidence is incomplete.

Validation Checklist

A tailings water-balance and seepage-control calculation is ready for review only when it states:

  1. control-volume boundary and sign convention;
  2. dry-solids, slurry, pond, seepage, reclaim, rainfall, evaporation, and discharge bases;
  3. source of stage-storage, beach, density, and climate data;
  4. seepage measurement method, hydraulic-gradient basis, and drain condition;
  5. freeboard requirement, pond setback requirement, and operating trigger level;
  6. uncertainty or sensitivity for dominant measured inputs;
  7. water-quality constituents relevant to the site geochemistry;
  8. monitoring evidence needed before continuing deposition, raising the embankment, drawing down the pond, or closing an action item.

The practical release question is: does the calculation explain the observed pond level, seepage flow, drain response, piezometer trend, and water-quality evidence with enough margin for the current operating mode? If not, the result is a prompt for investigation, not a release package.

REF

See also