This article explains how to analyze hydrogen sulfide formation (sometimes referred to as septicity) in a sewer hydraulic model.
Hydrogen sulfide forms in wastewater collection systems when anoxic conditions exist and results in corrosion, odor and toxicity problems. The rate of formation depends primarily on the strength of the wastewater as characterized by BOD concentration (Biochemical Oxygen Demand) and rate constants for the reactions which are dependent on temperature. Hydrogen sulfide in most commonly a problem in systems with long detention times in warm weather. The hydrogen sulfide calculations assume that sulfate concentration is not limiting.
The calculations are based on the Pomeroy-Parkhurst equation for predicting H2S concentrations in a sewer network; they can be used to evaluate the build-up or decay of the H2S concentration in a sewer system. The sewer system can be gravity and pressure combined system consisting of conduits, ponds, channels, wet-wells, pumps, pressure pipes and junctions.
Note that the calculated H2S produced by the model represents the H2S concentration in the aqueous phase. The amount that escapes to the gas phase depends on the turbulence, and various coefficients in Henry’s Law.
The Pomeroy-Parkhurst equation used for the calculation was developed on steady flow conditions. Since most models using SWMM solver are of dynamic modeling, some assumptions must be made for proper application of the equation. The primary assumption is that the solver will use the average hydraulic condition for the H2S calculations. The solver determines the average pipe flow by the total flow volume that has passed through the pipe for the simulation duration, the node H2S/BOD mixing calculations will also be based on the average flows of incoming and outlet pipes.
For more background information, see the Help topic "Hydrogen-Sulfide (H2S) Modeling" and the video further down.
To analyze hydrogen sulfide formation in a SewerGEMS model, the SWMM solver must be used and various data must be input into the manholes and pipes.
The primary assumption is that the solver will use the average hydraulic condition for the H2S calculations. The solver determines the average pipe flow by the total flow volume that has passed through the pipe for the simulation duration, the node H2S/BOD mixing calculations will also be based on the average flows of incoming and outlet pipes. The pipe detention time is determined by the pipe length and average flow velocity; the node detention time is determined by the average node volume and average outlet flow.
As a result of these assumptions the model provides a steady result set for the network which represents the expected average changes.
1) Select the Explicit (SWMM) solver in the Calculation Options or Compute Center
2) Select "True" for the calculation option "Run Hydrogen Sulfide Analysis." A reasonable model run time and a small reporting time step increment (such as 0.05 hours) should be set.
3) For each manhole, enter:
H2S (Local Inflow)BOD (local inflow)
Usually, there is very little H2S in the inflow but the presence of BOD (biochemical oxygen demand), along with sulfate, leads to the formation of H2S. For domestic wastewater, the BOD concentration may be on the order of 200 to 300 mg/L while for industrial waste, it can vary widely.
4) For each pipe, enter:
TemperatureH2S Flux CoefficientH2S Loss Factor
5) For wetwells, enter:
H2S (Local Inflow)TemperatureReaction rateBOD (local inflow)6) Compute the model and view the computed H2S and BOD concentration in the "Results (H2S)" section of the properties. The results will not change with time. This is because the average flow conditions are applied to the H2S calculations.
H2S forms in link elements (conduits, pressure pipes and channels) and is a function of temperature. The temperature default value is 20 deg C. The H2S flux concentration is on the order of 0.0003/hr while the H2S loss rate is on the order of 0.9. These values need to be calibrated for local conditions as there can be a great deal of variation. For pressure pipe and inverted siphons, there is no loss coefficient as H2S gas cannot leave the liquid phase. For wet wells, the Reaction Rate should be on the order of 0.2/day.
See also : Troubleshooting "N/A" results in H2S Modeling for pressure pipes