CN105760695A - Method for predicting hydrogen sulfide output in drainage pipeline - Google Patents

Method for predicting hydrogen sulfide output in drainage pipeline Download PDF

Info

Publication number
CN105760695A
CN105760695A CN201610136984.5A CN201610136984A CN105760695A CN 105760695 A CN105760695 A CN 105760695A CN 201610136984 A CN201610136984 A CN 201610136984A CN 105760695 A CN105760695 A CN 105760695A
Authority
CN
China
Prior art keywords
hydrogen sulfide
drainage pipeline
biomembrane
parameter
aerobic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610136984.5A
Other languages
Chinese (zh)
Other versions
CN105760695B (en
Inventor
江峰
梁爽
梁振声
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China Normal University
Original Assignee
South China Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by South China Normal University filed Critical South China Normal University
Priority to CN201610136984.5A priority Critical patent/CN105760695B/en
Publication of CN105760695A publication Critical patent/CN105760695A/en
Application granted granted Critical
Publication of CN105760695B publication Critical patent/CN105760695B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Purification Treatments By Anaerobic Or Anaerobic And Aerobic Bacteria Or Animals (AREA)
  • Activated Sludge Processes (AREA)

Abstract

The invention discloses a method for predicting hydrogen sulfide output in a drainage pipeline. The method comprises the following steps: carrying out intensive sampling characteristic analysis on sewage in the drainage pipeline and simultaneously obtaining pipeline parameters of the drainage pipeline; (S2) building a water quality module, a dynamic biological membrane module and a liquid-gas diffusion module on the basis of ASM3 model; (S3) obtaining simulating parameters, input by a user, of a production mode and an initial value of a biological membrane; (S4) running a prediction model to obtain stable biological membrane parameters; (S5) calculating to obtain the concentration of sulfide and a hydrogen sulfide gas in effluent from the tail end of the drainage pipeline after each period of time; (S6) calibrating the prediction model; and (S7) carrying out real-time prediction on the output of hydrogen sulfide in the drainage pipeline by the calibrated prediction model. According to the method, biological dynamic change in the drainage pipeline can be simulated; the output of the hydrogen sulfide can be accurately predicted and obtained; and the method can be widely applied to the field of production prediction of pollution gas.

Description

The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline
Technical field
The dusty gas that the present invention relates to supply and drain water system produces prediction field, particularly relates to hydrogen sulfide in drainage pipeline The Forecasting Methodology of generation amount.
Background technology
Explanation of nouns:
COD: full name Chemical Oxygen Demand, COD, be chemically to measure in water sample to need The amount of oxidized reducing substances;
H2S: hydrogen sulfide;
NH3: ammonia;
CH4: methane;
SRB: sulfate reducting bacteria;
ASM: activated sludge series model (Activated Sludge Models:ASM1, ASM2, ASM2D and ASM3);
T: temperature;
VFA: VFA;
DO: dissolved oxygen;
SO4 2-: sulfate;
DS: dissolubility sulfide;
TS: total sulfur compound;
TCOD: total COD;
SCOD: resolvable chemical oxygen demand;
VSS: volatile suspended matter;
NH4-N: ammonia nitrogen;
NOx-N: nitrate nitrogen;
Q: flow;
HRT: hydraulic detention time;
ASM model components related terms is explained:
SO: soluble component includes dissolved oxygen;
SF: easily biological-degradable dissolved organic matter;
SI: inertia dissolved organic matter;
SNH:Ammonium salt and ammonia nitrogen;
N2: dinitrogen;
SNO: nitrate nitrogen;
SALK: basicity;
SSO4: sulfate;
SH2S: Sulfide in water;
XI: graininess inertia Organic substance;
XS: degradable chemical at a slow speed;
XH: Heterotroph amount;
XSTO: the intracellular storage product that different oxygen is biological;
XA: nitrification is biological;
XTS: float;
XINA: inerting microorganism;
XSRB: sulfate reducting bacteria;
XNR-SOB: denitrification desulfurizing bacteria.
Urban drainage pipe network is the requisite ingredient of urban infrastructure.In recent years.Along with quickly sending out of city Exhibition, the scale of drainage pipeline networks is also gradually increased, however during sewage transports, under biological chemistry action, some materials Can occur with harmful gaseous form, such as H2S、NH3、CH4Etc..Wherein, hydrogen sulfide H2S can cause water quality to be disliked as one Change, the smelly odorant pollutant of blackening, the most also there is toxicity, corrosivity, pipe-line maintenance etc. had threat.Hydrogen sulfide Threat includes: 1. H2S causes drainage pipeline to corrode.One research of houston, U.S.A Drainage Services Department shows, this city 70% impaired Drainage pipeline is caused by hydrogen sulfide corrosion;And Belgian Flanders, the loss caused due to hydrogen sulfide corrosion reaches 500 Ten thousand Euros/year, are about the 10% of this city annual sewage collecting processing cost;In China Shanghai City, the corrosion of hydrogen sulfide is also Cause one of main cause that drainage pipeline damages.②H2S has stench, spills into and stink can be caused in air to pollute.③H2S has Strong neurotoxicity, can cause people's severe acute pancreatits going into a coma rapidly and even die suddenly, in Shanghai, Zhejiang, Jiangsu, the edge such as Guangdong Haiti district the most a lot of generation sewer maintenance workman is poisoned to death event.
In general, life of urban resident sewage does not also contain sulfide, and the sulfide in blow-off pipe is that sewage is at pipe In road transportation, the SRB of the biomembrane internal layer of tube wall produce, and in dispersal events to air.Sewer hydrogen sulfide Generation be a complex process, it is made by temperature, organic concentration, sulfate concentration, biofilm structure, multiple-microorganism With, the impact of dissolved oxygen concentration, nitrate concentration, hydraulic detention time, effluent flow rate etc. factors.Meanwhile, city row Grid has the highest complexity, and water quality and quantity has unstability, therefore H2The generation position of S, concentration, time have very High uncertainty, and present the most day and the most night fluctuation and seasonal fluctuation feature.
Currently the control to sewer stench gaseous contaminant mainly take aeration aerating, add nitrate, sodium hydroxide, The methods such as iron salt, but owing to the water yield in drainage pipeline networks, variation water quality are big, H2The stench gaseous contaminants such as S at different time, no Very big with the concentration difference in sewage pipeline section, which results in reagent cost too high or control the phenomenon of poor efficiency of effect.
In order to realize the fixing quantity of sewer hydrogen sulfide, the generation of hydrogen sulfide in founding mathematical models, prediction sewer, Exhaust position and concentration are crucial.The Related Mathematical Models being predicted the generation amount of sewer hydrogen sulfide, currently mainly has The WATS model of Aalborg university of Denmark exploitation and the SewerX model of University of Queensland's exploitation, but above-mentioned model Existing defects, it is impossible to simulation sewer biodynamic change, does not embodies biomembrane and produces sewer hydrogen sulfide and conversion Material impact, thus be accordingly used in hydrogen sulfide product row's prediction has inconvenience more, it is impossible to prediction obtains the discharge capacity of hydrogen sulfide exactly.
Summary of the invention
In order to solve above-mentioned technical problem, it is an object of the invention to provide the prediction of hydrogen sulfide generation amount in drainage pipeline Method.
The technical solution adopted for the present invention to solve the technical problems is:
The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline, comprises the following steps:
S1, the sewage in drain pipe is carried out intensive sampling specificity analysis, record real-time influent quality parameter and water outlet Water quality parameter, obtains the pipe parameter of this drain pipe simultaneously;
S2, based on ASM3 model, set up water quality module, dynamic biological film module and liquid-gas diffusion module, formed hydrogen sulfide The forecast model of generation amount, and using the pipe parameter of influent quality water amount parameters and drain pipe as the input of this forecast model;
S3, the analog parameter obtaining the forecast model of user's input and biomembrane initial value;
S4, operation forecast model, use water quality module to carry out water quality parameter calculating, and use dynamic biological film module to carry out The simulation of biomembranous dynamic change calculations, physical action and biological chemistry action simulation, so stable arriving pipeline biomembrane After date, it is thus achieved that stable biomembrane parameter;
S5, using pipe parameter and real-time influent quality water amount parameters as the input of forecast model, run forecast model The concentration of material in calculating acquisition gas and water and biomembrane, and by setting time output drainage pipeline end at set intervals The sulfide of water outlet and sulfureted hydrogen gas concentration;
In the effluent quality water amount parameters that S6, the sulfureted hydrogen gas concentration in each moment obtained according to prediction are corresponding Sulfureted hydrogen gas concentration between difference, forecast model is calibrated;
Forecast model after S7, employing calibration carries out real-time estimate to the generation amount of hydrogen sulfide in drainage pipeline.
Further, described influent quality water amount parameters and effluent quality water amount parameters all include pH value, temperature and volatilization Property organic acid, dissolved oxygen, sulfate, dissolubility sulfide, total sulfur compound, total COD, resolvable chemical oxygen demand, wave The property sent out float, ammonia nitrogen, nitrate nitrogen and the concentration of gas phase hydrogen sulfide gas, the pipe parameter of described drain pipe includes pipeline stream State, pipe shape, pipe range, caliber, gradient, roughness and pipeline flow.
Further, water quality module described in described step S2 is used for calculating organic hydrolysis, S in aqueous phaseFAerobic storage Tibetan, SFAnaerobism storage, XHAerobic growth, XHAnaerobic growth, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XAAerobic respiration and XAThe microbial reaction speed of anaerobic respiration, trip temperature school of going forward side by side Just.
Further, in described step S2, dynamic biological film module is set up in the following manner:
Every class the dissolved material in biomembranous hierarchy number, thickness in monolayer, each layering under S21, acquisition original state And after the every class granulating material in each layering, calculate in conjunction with following formula, and constantly adjust thickness in monolayer, each layering In every class the dissolved material and every class granulating material in each layering, until obtaining dynamic biological film after meeting following formula Model:
ϵ l + Σ i = 1 n ϵ X i , j f = 1 ϵ l = 0.2 + ( 0.5 500 * 0.25 ( 10 3 / ( 10 6 · z + 1 ) ) + 1 ) ϵ X i , j f = X i , j f ρ i f
In above formula, j represents biomembranous hierarchy number, and n represents biomembranous total hierarchy number, and z represents biomembrane depth of seam division, △ZjRepresent biomembranous thickness in monolayer,Representing the i-th class granulating material in jth layering, i, j, z are positive integer, εlTable Show the volume fraction of liquid in jth layering,Represent the volume fraction that granulating material is shared in jth is layered,Represent Granulating materialCOD density;
S22, set up granulating material and exist at the adhesion formula within biomembrane, peel rate formula and the dissolved material Diffusion formula within biomembrane;
S23, set up biomembrane in organic hydrolysis, SFAerobic storage, SFAnaerobism storage, aerobic growth, anaerobism raw Long, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XAAerobic respiration, XAAnaerobic respiration, SRB growth, SRB decline, heterotrophic bacteria inactivation, XSTOPassivation, autotrophic bacteria inactivation, SRB inactivation, sulfide oxygen Change, ferment and the microbial reaction formula of reoxygenation, set up the reaction equation of sulfate and sulfide simultaneously.
Further, in described step S21, the adjusting range of biomembranous thickness in monolayer is 0.5 × 10-4~1.5 × 10-4m。
Further, described liquid-gas diffusion module is for according to Henry's law and unit conversion method, in conjunction with dynamic biological Film module calculates the dissolved hydrogen sulfide obtained, and calculates and obtains hydrogen sulfide in gas phase gas concentration.
Further, described step S7, including:
S71, acquisition pipe parameter and real-time influent quality water amount parameters, as the input of forecast model, obtain simultaneously and set Fixed simulated conditions parameter, and then it is predicted model initialization;
S72, calculate the hydraulic parameters in drain pipe according to pipe parameter;
The content of reoxygenation in S73, calculating aqueous phase;
S74, employing water quality module calculate organic hydrolysis, S in aqueous phaseFAerobic storage, SFAnaerobism storage, XHGood Oxide growth, XHAnaerobic growth, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, nitrification are made With, XAAerobic respiration and XAAnaerobic respiration biochemical reaction after obtain the reaction rate of each material, trip temperature of going forward side by side Correction;
S75, employing dynamic biological film module carry out biomembranous dynamic change calculations, physical action simulation and biochemistry After effect simulation calculates, it is thus achieved that the sulfide of drainage pipeline outlet water at tail end and sulfureted hydrogen gas concentration after at set intervals.
Further, described step S75, including:
S751, employing dynamic biological film module calculate each material in aqueous phase to carry out with biomembrane entry material adhering to, peeling off And the physical action reaction of diffusion;
S752, calculate the microbial reaction of each material in biomembrane and go forward side by side after trip temperature correction, it is thus achieved that sulfate and dissolving The total amount of state sulfide;
S753, according to Henry's law and unit conversion method, calculate the dissolved obtained in conjunction with dynamic biological film module Hydrogen sulfide, calculates and obtains hydrogen sulfide in gas phase gas concentration.
The invention has the beneficial effects as follows: the Forecasting Methodology of hydrogen sulfide generation amount in the drainage pipeline of the present invention, including: S1, Sewage in drain pipe is carried out intensive sampling specificity analysis, records real-time influent quality water amount parameters and the effluent quality water yield Parameter, obtains the pipe parameter of this drain pipe simultaneously;S2, based on ASM3 model, set up water quality module, dynamic biological film module and Liquid-gas diffusion module, forms the forecast model of hydrogen sulfide generation amount, and is joined by the pipeline of influent quality water amount parameters and drain pipe Number is as the input of this forecast model;S3, the analog parameter obtaining the forecast model of user's input and biomembrane initial value; S4, operation forecast model, use water quality module to carry out water quality parameter calculating, and it be biomembranous to use dynamic biological film module to carry out Dynamically change calculations, physical action simulation and biological chemistry action simulation, and then after arriving pipeline biomembrane stable phase, it is thus achieved that Stable biomembrane parameter;S5, using pipe parameter and real-time influent quality water amount parameters as the input of forecast model, calculate Obtain sulfide and the sulfureted hydrogen gas concentration of rear drainage pipeline outlet water at tail end at set intervals;S6, according to prediction obtain The difference between sulfureted hydrogen gas concentration in the effluent quality water amount parameters that the sulfureted hydrogen gas concentration in each moment is corresponding Value, calibrates forecast model;Forecast model after S7, employing calibration carries out reality to the generation amount of hydrogen sulfide in drainage pipeline Time prediction.This method can simulate the biodynamic change in drainage pipeline, and prediction accurately obtains the generation amount of hydrogen sulfide, On the basis of this method, it is possible to achieve efficient, low cost, exactly control municipal drainage pipeline H2S gas concentration, reduces pipeline Corrosion, extend pipeline life, reduce maintenance pipeline time hydrogen sulfide poisoning risk, avoid hydrogen sulfide stench to resident living do Disturb.
Accompanying drawing explanation
The invention will be further described with embodiment below in conjunction with the accompanying drawings.
Fig. 1 be the present invention drainage pipeline in hydrogen sulfide generation amount Forecasting Methodology a specific embodiment flow process signal Figure;
Fig. 2 be the present invention a specific embodiment in dynamic biological film module set up process schematic;
Fig. 3 be the present invention a specific embodiment in biochemical reaction schematic diagram in drainage pipeline;
Fig. 4 be the present invention a specific embodiment in run forecast model and carry out the detailed process of hydrogen sulfide Output forecast Figure;
Fig. 5 be the present invention a specific embodiment in predict the sulfureted hydrogen gas concentration of acquisition and the change curve of measured value Schematic diagram;
Fig. 6 be the present invention a specific embodiment in predict the change curve signal of the total sulfur compound of acquisition and measured value Figure;
Fig. 7 be the present invention a specific embodiment in use the change of forecast model simulation hydrogen sulfide emission under normal circumstances Change curve chart;
Fig. 8 be the present invention a specific embodiment in use the hydrogen sulfide row in the case of forecast model simulative optimization scheme one The change curve put;
Fig. 9 be the present invention a specific embodiment in use the hydrogen sulfide row in the case of forecast model simulative optimization scheme two The change curve put;
Figure 10 be the present invention a specific embodiment in use the hydrogen sulfide in the case of forecast model simulative optimization scheme three The change curve of discharge.
Detailed description of the invention
The invention provides the Forecasting Methodology of hydrogen sulfide generation amount in a kind of drainage pipeline, comprise the following steps:
S1, the sewage in drain pipe is carried out intensive sampling specificity analysis, record real-time influent quality parameter and water outlet Water quality parameter, obtains the pipe parameter of this drain pipe simultaneously;
S2, based on ASM3 model, set up water quality module, dynamic biological film module and liquid-gas diffusion module, formed hydrogen sulfide The forecast model of generation amount, and using the pipe parameter of influent quality water amount parameters and drain pipe as the input of this forecast model;
S3, the analog parameter obtaining the forecast model of user's input and biomembrane initial value;Detailed, wherein simulate ginseng Number includes: the step-length of duct segments simulation, the step of timesharing simulation, simulation total time, sampling interval, sampling number of times, biomembrane are steady Periodically;Biomembrane initial value includes: the concentration of biofilm thickness, all kinds of the dissolved and granulating material and ratio, biomembrane are close Degree;
S4, operation forecast model, use water quality module to carry out water quality parameter calculating, and use dynamic biological film module to carry out The simulation of biomembranous dynamic change calculations, physical action and biological chemistry action simulation, so stable arriving pipeline biomembrane After date, it is thus achieved that stable biomembrane parameter;
S5, using pipe parameter and real-time influent quality water amount parameters as the input of forecast model, run forecast model The concentration of material in calculating acquisition gas and water and biomembrane, and by setting time output drainage pipeline end at set intervals The sulfide of water outlet and sulfureted hydrogen gas concentration;
In the effluent quality water amount parameters that S6, the sulfureted hydrogen gas concentration in each moment obtained according to prediction are corresponding Sulfureted hydrogen gas concentration between difference, forecast model is calibrated;
Forecast model after S7, employing calibration carries out real-time estimate to the generation amount of hydrogen sulfide in drainage pipeline.
Being further used as preferred embodiment, described influent quality water amount parameters and effluent quality water amount parameters all include PH value, temperature and VFA, dissolved oxygen, sulfate, dissolubility sulfide, total sulfur compound, total COD, molten Solution property COD, volatile suspended matter, ammonia nitrogen, nitrate nitrogen and the concentration of gas phase hydrogen sulfide gas, the pipe of described drain pipe Road parameter includes pipeline fluidised form, pipe shape, pipe range, caliber, gradient, roughness and pipeline flow.
Being further used as preferred embodiment, described in described step S2, water quality module is used for calculating Organic substance in aqueous phase Hydrolysis, SFAerobic storage, SFAnaerobism storage, XHAerobic growth, XHAnaerobic growth, XHAerobic respiration, XHAnaerobism exhale Suction, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XAAerobic respiration and XAThe microbial reaction of anaerobic respiration Speed, trip temperature of going forward side by side corrects.
Being further used as preferred embodiment, in described step S2, dynamic biological film module is to set up in the following manner :
Every class the dissolved material in biomembranous hierarchy number, thickness in monolayer, each layering under S21, acquisition original state And after the every class granulating material in each layering, calculate in conjunction with following formula, and constantly adjust thickness in monolayer, each layering In every class the dissolved material and every class granulating material in each layering, until obtaining dynamic biological film after meeting following formula Model:
ϵ l + Σ i = 1 n ϵ X i , j f = 1 ϵ l = 0.2 + ( 0.5 500 * 0.25 ( 10 3 / ( 10 6 · z + 1 ) ) + 1 ) ϵ X i , j f = X i , j f ρ i f
In above formula, j represents biomembranous hierarchy number, and n represents biomembranous total hierarchy number, and z represents biomembrane depth of seam division, △ZjRepresent biomembranous thickness in monolayer,Representing the i-th class granulating material in jth layering, i, j, z are positive integer, εlTable Show the volume fraction of liquid in jth layering,Represent the volume fraction that granulating material is shared in jth is layered,Represent Granulating materialCOD density;
S22, set up granulating material and exist at the adhesion formula within biomembrane, peel rate formula and the dissolved material Diffusion formula within biomembrane;
S23, set up biomembrane in organic hydrolysis, SFAerobic storage, SFAnaerobism storage, aerobic growth, anaerobism raw Long, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XAAerobic respiration, XAAnaerobic respiration, SRB growth, SRB decline, heterotrophic bacteria inactivation, XSTOPassivation, autotrophic bacteria inactivation, SRB inactivation, sulfide oxygen Change, ferment and the microbial reaction formula of reoxygenation, set up the reaction equation of sulfate and sulfide simultaneously.
Being further used as preferred embodiment, in described step S21, the adjusting range of biomembranous thickness in monolayer is 0.5 ×10-4~1.5 × 10-4m。
Being further used as preferred embodiment, described liquid-gas diffusion module is for changing according to Henry's law and unit Calculation method, calculates the dissolved hydrogen sulfide obtained in conjunction with dynamic biological film module, calculates and obtains hydrogen sulfide in gas phase gas concentration.
It is further used as preferred embodiment, described step S7, including:
S71, acquisition pipe parameter and real-time influent quality water amount parameters, as the input of forecast model, obtain simultaneously and set Fixed simulated conditions parameter, and then it is predicted model initialization;
S72, calculate the hydraulic parameters in drain pipe according to pipe parameter;
The content of reoxygenation in S73, calculating aqueous phase;
S74, employing water quality module calculate organic hydrolysis, S in aqueous phaseFAerobic storage, SFAnaerobism storage, XHGood Oxide growth, XHAnaerobic growth, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, nitrification are made With, XAAerobic respiration and XAAnaerobic respiration biochemical reaction after obtain the reaction rate of each material, trip temperature of going forward side by side Correction;
S75, employing dynamic biological film module carry out biomembranous dynamic change calculations, physical action simulation and biochemistry After effect simulation calculates, it is thus achieved that the sulfide of drainage pipeline outlet water at tail end and sulfureted hydrogen gas concentration after at set intervals.
It is further used as preferred embodiment, described step S75, including:
S751, employing dynamic biological film module calculate each material in aqueous phase to carry out with biomembrane entry material adhering to, peeling off And the physical action reaction of diffusion;
S752, calculate the microbial reaction of each material in biomembrane and go forward side by side after trip temperature correction, it is thus achieved that sulfate and dissolving The total amount of state sulfide;
S753, according to Henry's law and unit conversion method, calculate the dissolved obtained in conjunction with dynamic biological film module Hydrogen sulfide, calculates and obtains hydrogen sulfide in gas phase gas concentration.
Below in conjunction with specific embodiment, the present invention is elaborated.
Embodiment one
The schematic flow diagram of reference Fig. 1, the Forecasting Methodology of hydrogen sulfide generation amount in a kind of drainage pipeline, including:
S1, the sewage in drain pipe is carried out intensive sampling specificity analysis, record real-time influent quality water amount parameters and Effluent quality water amount parameters, obtains the pipe parameter of this drain pipe simultaneously;Particularly as follows: every 2h gather respectively a water inlet water sample and Water outlet water sample, continuous 48h intensive sampling, wherein water inlet water sample is used for being predicted the calibration of model;Water outlet water sample is used for verifying The emulator of model, uses on-line instrument odalog 7000 to measure the H produced in water inlet water sample and water outlet water sample simultaneously2S's Gas concentration.When sampling, sample frequency and time can the most arbitrarily adjust.Influent quality water amount parameters and water outlet water Matter water amount parameters refers to each component ratio of waste water in drain pipe, all includes pH value, T and VFA, DO, SO4 2-、DS、TS、 TCOD、SCOD、VSS、NH4-N, NOx-N and H2The concentration of S gas.The present embodiment gathers the influent quality water amount parameters of acquisition Meansigma methods as shown in table 1 below:
Influent quality water amount parameters meansigma methods in table 1 drainage pipeline
The pipe parameter of drain pipe includes pipeline fluidised form, pipe shape, pipe range, caliber, gradient, roughness and pipeline stream Amount.Pipeline fluidised form refers to that in pipeline, the form of sewage stream is belonging to pressure current or gravity stream, and pipe shape can be circular, square Shape etc..The pipe parameter of the drain pipe of the present embodiment is as shown in table 2 below:
The pipe parameter of table 2 drain pipe
Project Parameter value Parameter declaration/unit explanation
Pipeline fluidised form 2 1 is gravity stream, and 2 is pressure current
Pipe shape 2 1 is square, and 2 is circular
Pipe range 700 Unit m
Caliber 0.7 Unit m
Roughness 0.0012 --
Guan Kuan 0 Unit m
Guan Gao 0 Unit m
Gradient 0.005 --
S2, based on ASM3 model, utilize fortran95 high-level language that ASM3 model is write and revised, set up water Matter module, dynamic biological film module and liquid-gas diffusion module, form the forecast model of hydrogen sulfide generation amount, and by influent quality The pipe parameter of water amount parameters and drain pipe is as the input of this forecast model;
According to No. 3 models of activated sludge (ASM3 model), blow-off pipe waste water quality component is divided into: soluble component bag Include dissolved oxygen SO、SI、SF、SNH、N2、SNO、SALK;X is included with graininess componentI、XS、XH、XSTO、XA、XTS.And increase component XINA、 XSRB、SSO4And SH2SDeng, build water quality module, dynamic biological film module and liquid-gas diffusion module.
Measure in blow-off pipe sewage according to " water and effluent monitoring analyze method " (fourth edition, China Environmental Science Press) Water quality index: field assay measure water quality index include VFA, DO, pH, T;Residue water sample delivers to experiment after sealing cold preservation Room, completed SO in 24 hours4 2-、DS、TS、TCOD、SCOD、TSS、VSS、NH3-N、NOxThe analysis of the water quality index such as-N is surveyed Fixed, wherein the dissolved sulfide and total sulfur compound water sample to be measured press national standard method and are used aluminium hydroxide, zinc acetate-sodium acetate solid in advance Fixed;Concentration of hydrogen sulfide H in artificial well MH172S (g) uses field instrumentation (OdaLog 7000) on-line determination, locate into Distance water surface 1m elevation plane.Drain discharge, obtains according to pressure feed pump working state recording.It addition, heterotrophic bacteria, autotrophic bacteria and The yield coefficient of sulfate reducting bacteria, the decline parameter such as coefficient and maximum growth rate use document " A biofilm model For prediction of pollutant transformationin sewers " in reference value, be shown in Table 3 in detail:
The value of table 3 prediction model parameters
It practice, hydrogen sulfide gas generation process relates to solid, liquid, gas three-phase, water quality module, dynamic biological film in drain pipe Module and liquid-gas diffusion module corresponding sewage phase, solid biologic film and the reaction of gas phase respectively.
(1) water quality module: characteristic based on blow-off pipe sewage, this water quality module is used for calculating organic water in aqueous phase Solution, SFAerobic storage, SFAnaerobism storage, XHAerobic growth, XHAnaerobic growth, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XAAerobic respiration and XAAnaerobic respiration totally 12 kinds of microorganisms Microbial reaction speed, trip temperature of going forward side by side correct.The reaction rate of 12 kinds of microorganisms is anti-by 12 corresponding microorganisms Formula is answered to carry out calculating, similar with ASM3 model, formula 1~12 in its concrete formula such as table 4.In table 4, displaying is Biochemistry in drainage pipeline and dynamical matrix, because the reason of form length, the wherein reaction rate of last string Concrete formula uses code name formula 1~formula 21 to represent, formula is listed under form, specific as follows:
Table 4 biochemistry and dynamical matrix
In upper table, iNBMRepresent the nitrogen content in microorganism, iNXIRepresent XIIn nitrogen content.
Wherein, formula 1~formula 21 are as follows:
Formula 1:
Formula 2:
Formula 3:
Formula 4:
Formula 5:Formula 6:
Formula 7:
Formula 8:
Formula 9:
Formula 10:
Formula 11:
Formula 12:
Formula 13:
Formula 14:bSRBXSRB
Formula 15:
Formula 16:
Formula 17:
Formula 18:
Formula 19:
Formula 20:kSO·SO 0.1·SH2S
Formula 21:
Formula 22:kLa·(SO,sat-SO)
X in table 41,x2,…x20And y1,y2,…y21Can draw according to COD and N conservation equation in each course of reaction:
Σ i = 1 21 v k , i · l 1 , i = 0
Σ i = 1 21 v k , i · l 2 , i = 0
As sought x1 and y1 of hydrolysis:
Only see COD conservation:
Σ i = 1 21 v k , i · l 1 , i = v 1 , 2 · l 1 , 2 + v 1 , 5 · l 1 , 2 + 1 , 5 · l 1 , 9 = x 1 · 1 + y 1 · ( - 4.57 ) + ( - 1 ) · 1 = 0
Only see N conservation:
Σ i = 1 21 v k , i · l 2 , i = v 1 , 2 · l 2 , 2 + v 1 , 5 · l 2 , 2 + 1 , 5 · l 2 , 9 = x 1 · i N S S + y 1 · 1 + ( - 1 · i N B M ) = 0
Two formulas can solve x1 and y1.
It addition, the sulfide in sewage forms the sulfur not dissociated being generated by the reaction of biomembrane module and being diffused in aqueous phase Compound, according to ionic equilibrium Equation for Calculating, forms the hydrogen sulfide of the dissolved, as follows:
k 1 = [ HS - ] [ H + ] [ H 2 S ] , k 2 = [ S 2 - ] [ H + ] [ H S ]
Wherein, k1 is H2The one-level dissociation constant of S;[HS-] it is sulfur hydrogen radical ion (HS-) ion concentration mg/L, [H+] it is H+Ion concentration mg/L of ion;[H2S] for not dissociating hydrogen sulfide molecule (H2S) concentration;K2 is H2Two grades of dissociation constants of S; [S2-] for bearing the sulphion (S of bivalence2-) ion concentration.
(2) dynamic biological film module: in blow-off pipe, microorganism particularly SRB has tremendous influence for the content of sulfide Effect, and SRB is mainly distributed in the biomembrane of tube wall the deep layer of (or bed mud), so biological film model is for blow-off pipe H2S pollutes essential and plays a major role.Based on biomembrane for blow-off pipe H2The tremendous influence degree that S produces, this is dynamically Biomembrane module includes biomembrane and the reaction of complicated physical biochemistry that can dynamically change.Specifically include:
1) biomembranous dynamic change: in order to embody biomembranous heterogeneity, needs to move biofilm density on vertical State change is simulated.At the density method of this model with reference to Horn pattern: density includes mass density and COD density two kinds. Certain graininess in being layered for a certain biomembrane is organic, and its COD density has maximum to limit, actual concentrations and maximum limit The ratio of value processed is the plot ratio of this material, and the volume fraction of various materials is 1 with the liquid volume fraction sum in this layering, Mass density then calculates according to TS.Owing to the COD density threshold limit value of various particulate matters is different, although plot ratio summation is protected Hold constant, but the change that in being layered along with biomembrane, various microorganisms are constituted, mass density can embody dynamic Vertical Difference.
Dynamic biological film module to set up process as follows: with reference to Fig. 2, obtain biomembranous hierarchy number j, list under original state Layer thickness △ Zj, every class the dissolved material S in each layeringi,jAnd the every class granulating material X in each layeringi,jAfter, Calculate in conjunction with following formula, and constantly adjust thickness in monolayer △ Zj, every class the dissolved material S in each layeringi,jAnd it is each Every class granulating material X in layeringi,j, until obtaining dynamic biological membrane modle after meeting following formula:
ϵ l + Σ i = 1 n ϵ X i , j f = 1 ϵ l = 0.2 + ( 0.5 500 * 0.25 ( 10 3 / ( 10 6 · z + 1 ) ) + 1 ) ϵ X i , j f = X i , j f ρ i f
In above formula, j represents biomembranous hierarchy number, and n represents biomembranous total hierarchy number, and z represents biomembrane depth of seam division, △ZjRepresent biomembranous thickness in monolayer,Representing the i-th class granulating material in jth layering, i, j, z are positive integer, εlTable Show the volume fraction of liquid in jth layering,Represent the volume fraction that granulating material is shared in jth is layered,Represent Granulating materialCOD density;
Biomembranous thickness in monolayer △ ZjAdjusting range be 0.5 × 10-4~1.5 × 10-4M, if △ is ZjLess than 0.5*10-4 Then it is incorporated to last layer, if △ is ZjMore than 1.5 × 10-4M is then divided into two-layer;
2) complicated physical biochemistry reaction: include physical action and biological chemistry action
Physical action, including adhering to, peeling off and the diffusion of the dissolved material
It is generally believed that the process that suspended particulate substance adheres to solid dielectric is mainly a physical process, hanging in water body Floating particles substrate concentration, hydraulics and surface of solids roughness have certain impact to this process, adhere to formula as follows:
r a t t , X i w = k a t t · X i w
In above formula,Being the granulating material Xi speed that adheres to the surface of solids from water body, unit is gCODm-3d-1Referring to the concentration of granulating material Xi in water body, unit is gCODm-3, wherein subscript w represents water body;kattIt is to adhere to speed Rate constant, unit is d-1
In pipeline, biomembrane peel rate is mainly affected by peel rate own and flow shear, uses below equation to enter Row emulation:
r det , X i = k det μ H L f 2 ( τ w - τ w min τ w min ) 2.5 tanh ( M i M T S S ) τ w = ρ w · g · S f · R h r det , X i , j = - r det , X i · X i , j / Σ j = 1 n X i , j .
In above formula,It is particulate matter XiPeel rate, unit is gCODm-3d-1;kdetBeing to peel off constant, unit is gm-5;μHBeing the maximum growth rate of heterotrophic microorganism, unit is d-1;LfBeing biofilm thickness, unit is m;τwIt is shearing force, Unit is N m-2;MiIt is particulate matter XiTS mass;And MTSIt it is then monoblock biomembranous TS mass;SfFriction gradient;RhFor water Power radius, unit is m;ρwBeing the mass density of water, unit is gm-3;And g is acceleration of gravity, unit is m2s-1It is Particulate matter X in biomembrane j layeringiPeel rate, unit is gCODm-3d-1;Xi,jIt is that in biomembrane j layering, particulate form is organic Thing XiConcentration.
Diffusion is the main path that the dissolved material enters within biomembrane.The equation using molecule diffusion is retouched State, it is assumed that this diffusion process defers to Fick First Law:
J i , j = - D f , i , j ∂ S i , j ∂ z j D f = D w · f D f D = ( 1 - 0.43 ρ m 0.92 11.19 + 0.27 ρ m 0.99 ) D f , i , j = D w , i · ( 1 - 0.43 ρ m , j 0.92 11.19 + 0.27 ρ m , j 0.99 )
In above formula, Ji,jIt is the dissolved material SiThe diffusion flux unit passing through biomembrane j layered boundary is gm-2s-1;Si,j It is the dissolved material SiIn the concentration of biomembrane j layer, unit is gm-3;Df,i,jIt is the dissolved material SiExpansion in biomembrane j layer Dissipating coefficient, unit is m2s-1;DwBeing the diffusion coefficient in water body, unit is m2s-1, and fDIt is the effective of solutes accumulation in biomembrane Coefficient;ρmFor biomembranous mass density, unit is gm-3
It addition, diffusion rate of dissolved oxygen coefficient is obtained by the following formula:
Dw,O2=4.864 × 10-13T2+2.880×10-11T+1.268×10-9
Wherein, T is water temperature, and unit is DEG C.
Biological chemistry action: according to the characteristic of drainage pipeline, different from aqueous biochemical chemical action, the life in biomembrane Thing chemical reaction is divided into organic hydrolysis, SFAerobic storage, SFAnaerobism storage, aerobic growth, anaerobic growth, XH's Aerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XAAerobic respiration, XAAnaerobism Breathe, SRB grows, SRB becomes feeble and die, heterotrophic bacteria inactivation, XSTOPassivation, autotrophic bacteria inactivation, SRB inactivation, sulfide-oxidation, fermentation and Reoxygenation totally 21 microbial reaction formulas, and 13 kinds of components include SO、SF、SI、SNH、SNO、SSO4、SH2S、XI、XS、XH、XSTO、XAWith XSRB, as shown in table 4 and Fig. 3.
Meanwhile, set up the reaction equation of sulfate and sulfide, as shown in table 5:
The biochemistry of table 5 sulfate and sulfide and dynamical matrix
Same, in upper table, iNBMIt is the nitrogen content in microorganism, iNXIIt is XIIn nitrogen content.
X in table 51,x4,x5And y1,y2,y3Can draw according to COD conservation equation in each course of reaction following:
Σ i = 1 5 v k , i · l 1 , i = 0
Σ i = 1 5 v k , i · l 2 , i = 0
Biochemical reaction is all relevant with temperature with diffusion coefficient, the temperature in actual sewer and modeling temperature, I.e. standard reaction temperature (20 DEG C) or typical coefficient temperature (25 DEG C) has difference, no matter aqueous phase module or biomembrane module is both needed to It it is modified, to use following formula to be modified:
In above formula, r is the emulation reaction rate of microorganism, r at a temperature of certain20℃It is the reaction rate under 20 degrees Celsius, θTFor Correction factor, in aqueous phase, the biological respinse factor takes 1.07, and in biomembrane, the biological respinse factor takes 1.03, reoxygenation reaction (formula 21) Take 1.024.
(3) by the dissolved H obtained in (1)2S is diffused in gas phase, when it reaches balance at liquid-gas interface, uses Henry's law and unit conversion draw H in gas phase2S pollutant levels, computing formula is as follows:
P * = C / H P * V = n R T C g = M / V
In above formula, P*Being dividing potential drop, unit is kpa;C is concentration in solution, and unit is mol/L;H is Henry's constant;V is gas Body volume, unit is m3;N is the amount of material, and unit is mol;R is gas constant;T is temperature, and unit is K;Cg is in gas phase Material concentration, unit is kg/m3;M is material mass, and unit is kg.
S3, the analog parameter of forecast model obtaining user's input and biomembrane initial value, defeated as forecast model Enter.Obvious, also it is the input of forecast model in conjunction with the pipe parameter of the drain pipe in table 2.The simulation obtained in the present embodiment Parameter is as shown in table 6.
Table 6 analog parameter
Parameter declaration Parameter value
The step-length of duct segments simulation, △ x, unit m 100
The step-length of pipeline timesharing simulation, △ t, unit s 3
Simulation total time, unit h 24
Sampling interval, unit h 2
Sampling number of times (is not counted in) for the first time, unit time 12
Biomembrane stable phase, unit h 90*24
S4, operation forecast model, use water quality module to carry out water quality parameter calculating, and use dynamic biological film module to carry out The simulation of biomembranous dynamic change calculations, physical action and biological chemistry action simulation, so stable arriving pipeline biomembrane After date, it is thus achieved that stable biomembrane parameter.Particularly as follows: operation forecast model, water quality module is used to carry out water quality parameter calculating, Material concentration after physics, chemistry, biological change in aqueous phase, uses liquid-gas diffusion module to carry out hydrogen sulfide in gas The simulation of concentration, then use dynamic biological film module to carry out biomembranous dynamic change calculations, physical action simulation and bioid Effect is simulated, and obtains material concentration after physics, chemistry, biological change in biomembrane, then carries out material expansion with aqueous phase Dissipate and exchange, through long-play, after the biomembrane stabilization time reaching to set, it is thus achieved that stable biomembrane parameter.This step A rapid specific embodiment as shown in Figure 4, specifically includes:
1.1, program starts;
1.2, data type and the distribution array space of all kinds of parameter are defined;
1.3, the pipeline of field pipes, waterpower and water quality data are read, including in the pipeline conditional parameter of table 2 and table 1 The water yield and water quality index.Read the simulated conditions parameter of table 6, it is assumed that a length of ttmax during total simulation, biomembrane stable phase duration Tini, pipeline divides xmax section, and biomembrane is that (every layer at 0.5*10 for the dynamic ymax layer changed-4With 1.5*10-4Between m);
1.4, judge whether step 1.3 parameter fully enters complete, the most then enter step 1.5, if it is not, then return to step Rapid 1.3;
1.5, simulation starts, and simulated time sets t=0, forms whole piece according to the parameter of step 1.3 with model inner equation The environment initial value of pipeline, including water quality and biomembrane initial value, completes model initialization;
1.6, water force: pipeline condition and pipeline flow according to table 2 calculate, and t time xth pipeline section is (during beginning this moment Wait x=1, i.e. first paragraph pipeline) the depth of water, flow velocity, Reynolds number (Re). Floyd's algorithm (Fr), shearing force etc. to biomembrane adhere to With the hydraulics with material impact that comes off.
1.7, reoxygenation calculates: calculating the resume speed of DO in aqueous phase, DO content has a strong impact on sulfide concentration;
1.8, there is formula (1)~the biochemical reaction of (12) in table 4 in material in aqueous phase, and substance for calculation is in aqueous phase Concentration, altogether imax kind material.
1.9, the temperature correction of biochemical reaction: the microbial reaction speed temperature influence of step 1.8, model is initial It is set as 20 DEG C, is required to carry out temperature adjustmemt.
1.10, imax kind material and the biomembrane entry material in aqueous phase exported by step 1.9 swap include molten Solve the diffusion of state material, the adhesion of particulate matter material and peeling.
1.11, after step 1.10 completes, the liquid-gas diffusion carrying out hydrogen sulfide respectively is anti-with the biochemistry within biomembrane Should and mass exchange.On the one hand, in the liquid-gas diffusion mainly aqueous phase of hydrogen sulfide, a sulfide part forms the H2S not dissociated, The H2S that do not dissociates forms H2S gas in Henry's law is diffused into gas phase;On the other hand, the biochemistry within biomembrane is anti- Should and mass exchange mainly imax kind material layering biomembrane between carry out such as biochemical reaction in table 4, through excess temperature After correction, the biomembrane of the material concentration of output last layer respectively and next layer carries out mass exchange, then judges this layer of biomembrane The need of adjusting (if △ is ZjLess than 0.5*10-4Then it is incorporated to last layer;If △ is ZjMore than 1.5*10-4Then divide two-layer, altogether ymax Layer).
1.12, next pipeline section, the i.e. material concentration of x+1 section: return to step 1.6 are calculated, until all pipeline section has been simulated, I.e. pipeline section xmax calculates complete.
1.13, calculating next simulation moment, i.e. the material concentration of t+ △ t: return to step 1.5, pipeline sets x=1 Restart to calculate.As t >=tini, i.e. represent that pipeline water quality and biomembrane condition are considered to meet with actual pipeline, Fang Kejin Row next step, simultaneously simulation duration tt at time be 0.
1.14, output result: dense object time, target pipeline section, target organism film layer and target substance as required Degree output, the present embodiment is concentration of hydrogen sulfide in sulfide in output end pipeline aqueous phase and gas phase.
1.15, calculate next to simulate moment, i.e. the material concentration of tt+ △ t: return to step 1.6, until tt >=ttmax.
1.16, EP (end of program).
S5, using pipe parameter and real-time influent quality water amount parameters as the input of forecast model, calculate obtain every The sulfide of drainage pipeline outlet water at tail end and sulfureted hydrogen gas concentration after a period of time.
In the effluent quality water amount parameters that S6, the sulfureted hydrogen gas concentration in each moment obtained according to prediction are corresponding Sulfureted hydrogen gas concentration between difference, forecast model is calibrated, the parameter of calibration mainly includes μSRBAnd KSRB,O.Root It is predicted the H in the gentle phase of total sulfur compound in the sewage of model2S simulation result contrasts with measured value, by model parameter μSRBIt is corrected to 5.6d-1, KSRB,OIt is corrected to 1.0gO2m-3, simulation result is the most identical with measured data, such as Fig. 5 and 6.Prediction mould Type can emulate sulfide and H2The variation tendency of S gaseous contaminant and concentration.Fig. 5 and Fig. 6 illustrates at 2012.7.11 extremely 2012.7.13 the H of period actual measurement MH172S gas meansigma methods is 12.9ppm, and the meansigma methods of simulation curve is 12.5ppm.Mould Type can simulate S2-And H2Outside the concentration change of S, additionally it is possible to prediction H2The time to peak of S gas overflowing.
Forecast model after S7, employing calibration carries out real-time estimate to the generation amount of hydrogen sulfide in drainage pipeline, and simulation changes Become the nitrate dosage control effect to hydrogen sulfide, provide effective, economic control program.
Scheme emulates: change nitrate dosage H2The production of S
It is previously noted this pipeline hydrogen sulfide contamination serious, although the airport parties concerned have added calcium nitrate and administered (proportion 1.4 Calcium nitrate solution, daytime 70L/h, night 84L/h), but hydrogen sulfide contamination problem is the severeest.Dirty for effectively reducing hydrogen sulfide Dye, the scheme of the different calcium nitrate dosage of modeling three kinds: be respectively and 1. reduce by 50% dosage;2. increase by 50% dispensing Amount;3. increase by 100% dosage.Result is as follows:
If Fig. 7 is the hydrogen sulfide gas simulation curve that the forecast model after using correction runs 180h continuously, putting down of simulation All H2S concentration is 6.6ppm, and peak value is 33ppm, in this, as the background value controlling hydrogen sulfide simulation case.
1. prioritization scheme one, the nitrate dosage of minimizing 50%, such as Fig. 8, analog result shows, if dosage reduces 50%, then H2The mean concentration of S rises to 19.5ppm from 6.6.Under this result shows that the existence of nitrate and dosage are to controlling Water channel sulfide plays an important role
2. prioritization scheme two, the nitrate dosage of increase by 50%, such as Fig. 9, simulation curve represents when dosage increases When 50%, H2The mean concentration of S is effectively reduced, and is reduced to 1.9ppm from 6.6, reduces 71%, and all peak values are dense Degree is below 30ppm.Although the dosage of increase by 50% yet suffers from the hydrogen sulfide gas peak condition of higher concentration, but can Efficiently control H2The mean concentration of S.
3. prioritization scheme three, the nitrate dosage of increase by 100%, when dosage doubles, average H2S concentration from 6.6 are reduced to 1.1ppm, reduce 88%, such as Figure 10.H in most cases2The concentration of S reaches null value, but due to flow Change, however it remains H2The situation of S peak value.Obviously result shows that dosage doubles, and still there will be certain time H2S overflows The phenomenon gone out.Compared with the dosage increasing by 50%, regulation effect slightly rises, but cost is greatly improved, and is its 1.3 times.
The dosage being increased nitrate by forecast model simulation makes average H2S concentration is minimized, and adds the most The H produced2S concentration is the lowest, but for suppression H2S spillover effect is little.Relatively control H2S effect understands with dosage, increases The dosage adding 50% can reach the most economical and control hydrogen sulfide effect.
It is the present invention preferably embodiment that above-mentioned nitrate controls the embodiment of municipal administration blow-off pipe hydrogen sulfide contamination, this Bright modes such as adding oxygen, alkali and iron salt of also can simulating is to Sulfide in water and H2S gas is controlled effect emulation, even can Water quality such as COD, ammonia nitrogen, DO, SO to blow-off pipe4 2-Emulate etc. material concentration.
It is above the preferably enforcement of the present invention is illustrated, but the invention is not limited to described enforcement Example, those of ordinary skill in the art it may also be made that all equivalent variations on the premise of spirit of the present invention or replacing Changing, modification or the replacement of these equivalents are all contained in the application claim limited range.

Claims (8)

1. the Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline, it is characterised in that comprise the following steps:
S1, the sewage in drain pipe is carried out intensive sampling specificity analysis, record real-time influent quality parameter and effluent quality Parameter, obtains the pipe parameter of this drain pipe simultaneously;
S2, based on ASM3 model, set up water quality module, dynamic biological film module and liquid-gas diffusion module, formed hydrogen sulfide produce The forecast model of amount, and using the pipe parameter of influent quality water amount parameters and drain pipe as the input of this forecast model;
S3, the analog parameter obtaining the forecast model of user's input and biomembrane initial value;
S4, operation forecast model, use water quality module to carry out water quality parameter calculating, and use dynamic biological film module to carry out biology The simulation of the dynamic change calculations of film, physical action and biological chemistry action simulation, and then after arriving pipeline biomembrane stable phase, Obtain stable biomembrane parameter;
S5, using pipe parameter and real-time influent quality water amount parameters as the input of forecast model, run forecast model and calculate The concentration of material in acquisition gas and water and biomembrane, and by setting time output drainage pipeline outlet water at tail end at set intervals Sulfide and sulfureted hydrogen gas concentration;
Sulfur in the effluent quality water amount parameters that S6, the sulfureted hydrogen gas concentration in each moment obtained according to prediction are corresponding Change the difference between hydrogen gas concentration, forecast model is calibrated;
Forecast model after S7, employing calibration carries out real-time estimate to the generation amount of hydrogen sulfide in drainage pipeline.
The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline the most according to claim 1, it is characterised in that described water inlet Water quality and quantity parameter and effluent quality water amount parameters all include pH value, temperature and VFA, dissolved oxygen, sulfate, Dissolubility sulfide, total sulfur compound, total COD, resolvable chemical oxygen demand, volatile suspended matter, ammonia nitrogen, nitrate nitrogen, With the concentration of gas phase hydrogen sulfide gas, the pipe parameter of described drain pipe includes pipeline fluidised form, pipe shape, pipe range, caliber, slope Fall, roughness and pipeline flow.
The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline the most according to claim 1, it is characterised in that described step Water quality module described in S2 is used for calculating organic hydrolysis, S in aqueous phaseFAerobic storage, SFAnaerobism storage, XHAerobic life Long, XHAnaerobic growth, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XA Aerobic respiration and XAThe microbial reaction speed of anaerobic respiration, trip temperature of going forward side by side corrects.
The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline the most according to claim 1, it is characterised in that described step In S2, dynamic biological film module is set up in the following manner:
S21, the every class the dissolved material obtained under original state in biomembranous hierarchy number, thickness in monolayer, each layering and After every class granulating material in each layering, calculate in conjunction with following formula, and constantly adjust in thickness in monolayer, each layering Every class granulating material in every class the dissolved material and each layering, until obtaining dynamic biological film mould after meeting following formula Type:
ϵ l + Σ i = 1 n ϵ X i , j f = 1 ϵ l = 0.2 + ( 0.5 500 * 0.25 ( 10 3 / ( 10 6 · z + 1 ) ) + 1 ) ϵ X i , j f = X i , j f ρ i f
In above formula, j represents biomembranous hierarchy number, and n represents biomembranous total hierarchy number, and z represents biomembrane depth of seam division, △ Zj Represent biomembranous thickness in monolayer,Representing the i-th class granulating material in jth layering, i, j, z are positive integer, εlRepresent the The volume fraction of liquid in j layering,Represent the volume fraction that granulating material is shared in jth is layered,Represent granule State materialCOD density;
S22, set up granulating material at the adhesion formula within biomembrane, peel rate formula and the dissolved material at biology Diffusion formula within film;
S23, set up biomembrane in organic hydrolysis, SFAerobic storage, SFAnaerobism storage, aerobic growth, anaerobic growth, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XAAerobic respiration, XA's Anaerobic respiration, SRB growth, SRB decline, heterotrophic bacteria inactivation, XSTOPassivation, autotrophic bacteria inactivation, SRB inactivation, sulfide-oxidation, send out The microbial reaction formula of ferment and reoxygenation, sets up the reaction equation of sulfate and sulfide simultaneously.
The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline the most according to claim 1, it is characterised in that described step In S21, the adjusting range of biomembranous thickness in monolayer is 0.5 × 10-4~1.5 × 10-4m。
The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline the most according to claim 1, it is characterised in that described liquid- Gas diffusion module for according to Henry's law and unit conversion method, calculating the dissolved obtained in conjunction with dynamic biological film module Hydrogen sulfide, calculates and obtains hydrogen sulfide in gas phase gas concentration.
The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline the most according to claim 1, it is characterised in that described step S7, including:
S71, acquisition pipe parameter and real-time influent quality water amount parameters, as the input of forecast model, obtain setting simultaneously Simulated conditions parameter, and then it is predicted model initialization;
S72, calculate the hydraulic parameters in drain pipe according to pipe parameter;
The content of reoxygenation in S73, calculating aqueous phase;
S74, employing water quality module calculate organic hydrolysis, S in aqueous phaseFAerobic storage, SFAnaerobism storage, XHAerobic life Long, XHAnaerobic growth, XHAerobic respiration, XHAnaerobic respiration, XSTOAerobic respiration, XSTOAnaerobic respiration, Nitrification, XA Aerobic respiration and XAAnaerobic respiration biochemical reaction after obtain the reaction rate of each material, trip temperature school of going forward side by side Just;
S75, employing dynamic biological film module carry out biomembranous dynamic change calculations, physical action simulation and biological chemistry action After simulation calculates, it is thus achieved that the sulfide of drainage pipeline outlet water at tail end and sulfureted hydrogen gas concentration after at set intervals.
The Forecasting Methodology of hydrogen sulfide generation amount in drainage pipeline the most according to claim 7, it is characterised in that described step S75, including:
S751, use dynamic biological film module calculate each material and biomembrane entry material in aqueous phase carry out adhering to, peel off and The physical action reaction of diffusion;
S752, calculate the microbial reaction of each material in biomembrane and go forward side by side after trip temperature correction, it is thus achieved that sulfate and the dissolved sulfur The total amount of compound;
S753, according to Henry's law and unit conversion method, calculate the dissolved sulfuration obtained in conjunction with dynamic biological film module Hydrogen, calculates and obtains hydrogen sulfide in gas phase gas concentration.
CN201610136984.5A 2016-03-10 2016-03-10 The prediction technique of hydrogen sulfide yield in drainage pipeline Active CN105760695B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610136984.5A CN105760695B (en) 2016-03-10 2016-03-10 The prediction technique of hydrogen sulfide yield in drainage pipeline

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610136984.5A CN105760695B (en) 2016-03-10 2016-03-10 The prediction technique of hydrogen sulfide yield in drainage pipeline

Publications (2)

Publication Number Publication Date
CN105760695A true CN105760695A (en) 2016-07-13
CN105760695B CN105760695B (en) 2018-08-31

Family

ID=56333041

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610136984.5A Active CN105760695B (en) 2016-03-10 2016-03-10 The prediction technique of hydrogen sulfide yield in drainage pipeline

Country Status (1)

Country Link
CN (1) CN105760695B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106339586A (en) * 2016-08-23 2017-01-18 环境保护部华南环境科学研究所 Tracing method and prediction method of mustard gas pollution
CN111354425A (en) * 2020-05-08 2020-06-30 广西壮族自治区农业科学院 Method for evaluating soil acidification effect driven by plant harvesting
CN111982740A (en) * 2020-06-30 2020-11-24 中国水利水电科学研究院 Analysis method for migration and diffusion of water elements based on long-film sediment
CN112632755A (en) * 2020-12-03 2021-04-09 中煤科工集团重庆研究院有限公司 Comprehensive judgment and identification method for risk of hydrogen sulfide disaster of abandoned oil well in coal mining area
CN115563862A (en) * 2022-09-27 2023-01-03 盐城工学院 Pipeline hydrogen sulfide prediction method based on optimized combined prediction model

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102436232A (en) * 2011-10-18 2012-05-02 中国石油化工股份有限公司 Equipment maintaining system for factories for purifying natural gas with high content of sulfur
CN104007041A (en) * 2014-04-17 2014-08-27 河南理工大学 Device for determining content of hydrogen sulfide in coal bed

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102436232A (en) * 2011-10-18 2012-05-02 中国石油化工股份有限公司 Equipment maintaining system for factories for purifying natural gas with high content of sulfur
CN104007041A (en) * 2014-04-17 2014-08-27 河南理工大学 Device for determining content of hydrogen sulfide in coal bed

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
FENG JIANG ET AL.: "A biofilm model for prediction of pollutant transformation in sewers", 《WATER RESEARCH》 *
K SHARMA ET AL.: "PREDICTING HYDROGEN SULFIDE FORMATION IN SEWERS: A NEW MODEL", 《WATER》 *
ORI LAHAV ET AL.: "A different approach for predicting H2S(g) emission rates in gravity sewers", 《WATER RESEARCH》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106339586A (en) * 2016-08-23 2017-01-18 环境保护部华南环境科学研究所 Tracing method and prediction method of mustard gas pollution
CN106339586B (en) * 2016-08-23 2019-01-18 环境保护部华南环境科学研究所 Mustard gas pollutes source tracing method and forecasting procedure
CN111354425A (en) * 2020-05-08 2020-06-30 广西壮族自治区农业科学院 Method for evaluating soil acidification effect driven by plant harvesting
CN111982740A (en) * 2020-06-30 2020-11-24 中国水利水电科学研究院 Analysis method for migration and diffusion of water elements based on long-film sediment
CN112632755A (en) * 2020-12-03 2021-04-09 中煤科工集团重庆研究院有限公司 Comprehensive judgment and identification method for risk of hydrogen sulfide disaster of abandoned oil well in coal mining area
CN112632755B (en) * 2020-12-03 2024-04-30 中煤科工集团重庆研究院有限公司 Comprehensive judgment method for hydrogen sulfide disaster risk of abandoned oil well in coal mining area
CN115563862A (en) * 2022-09-27 2023-01-03 盐城工学院 Pipeline hydrogen sulfide prediction method based on optimized combined prediction model

Also Published As

Publication number Publication date
CN105760695B (en) 2018-08-31

Similar Documents

Publication Publication Date Title
CN105760695B (en) The prediction technique of hydrogen sulfide yield in drainage pipeline
Hvitved-Jacobsen et al. A process and model concept for microbial wastewater transformations in gravity sewers
Liu et al. Sulfide and methane production in sewer sediments
Carrera-Chapela et al. Modeling the odor generation in WWTP: an integrated approach review
Zhang et al. Effects of physical and biochemical processes on the dissolved oxygen budget for the Pearl River Estuary during summer
Carrera et al. A review of sulfide emissions in sewer networks: overall approach and systemic modelling
Liang et al. Systematic evaluation of a dynamic sewer process model for prediction of odor formation and mitigation in large-scale pressurized sewers in Hong Kong
Chaudhury et al. Dissolved oxygen modeling of the Blackstone River (northeastern United States)
CN103792844B (en) A kind of sewage treatment control method based on orderly cluster
Chau A three-dimensional eutrophication modeling in Tolo Harbour
Kalyuzhnyi et al. Integrated mathematical model of UASB reactor for competition between sulphate reduction and methanogenesis
Lung et al. Eutrophication model for the Patuxent estuary: Advances in predictive capabilities
CN104418428A (en) Water treatment control apparatus
Vollertsen et al. Integrated design of sewers and wastewater treatment plants
Wang et al. Development of a new sediment flux model–Application in Chesapeake Bay
Heydarzadeh et al. Dissolved oxygen determination in sewers using flow hydraulic parameters as part of a physical-biological simulation model
Chau An unsteady three-dimensional eutrophication model in Tolo harbour, Hong Kong
Hua et al. Parameter-efficient bioclogging model: calibration and comparison with laboratory data
Jones et al. Towards a useful dynamic model of the anaerobic wastewater treatment process: a practical illustration of process identification
Chau et al. Two-layered, 2D unsteady eutrophication model in boundary-fitted coordinate system
Vollertsen et al. Sewer quality modeling–a dry weather approach
Calabro et al. In sewer processes: mathematical model development and sensitivity analysis
Gostelow et al. Integrated odour modelling for sewage treatment works
Mei et al. Estimation of water environment capacity: example as four basin in Shandong province, China
Sharma et al. Modelling of methane production and emissions

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant