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

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；

H₂S: hydrogen sulfide；

NH₃: ammonia；

CH₄: 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；

SO₄ ^2-: sulfate；

DS: dissolubility sulfide；

TS: total sulfur compound；

TCOD: total COD；

SCOD: resolvable chemical oxygen demand；

VSS: volatile suspended matter；

NH₄-N: ammonia nitrogen；

NOx-N: nitrate nitrogen；

Q: flow；

HRT: hydraulic detention time；

ASM model components related terms is explained:

S_O: soluble component includes dissolved oxygen；

S_F: easily biological-degradable dissolved organic matter；

S_I: inertia dissolved organic matter；

S_NH:Ammonium salt and ammonia nitrogen；

N₂: dinitrogen；

S_NO: nitrate nitrogen；

S_ALK: basicity；

S_SO4: sulfate；

S_H2S: Sulfide in water；

X_I: graininess inertia Organic substance；

X_S: degradable chemical at a slow speed；

X_H: Heterotroph amount；

X_STO: the intracellular storage product that different oxygen is biological；

X_A: nitrification is biological；

X_TS: float；

X_INA: inerting microorganism；

X_SRB: sulfate reducting bacteria；

X_NR-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 H₂S、NH₃、CH₄Etc..Wherein, hydrogen sulfide H₂S 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. H₂S 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.②H₂S has stench, spills into and stink can be caused in air to pollute.③H₂S 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 H₂The 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, H₂The 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 phase_FAerobic storage Tibetan, S_FAnaerobism storage, X_HAerobic growth, X_HAnaerobic growth, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_AAerobic respiration and X_AThe 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:

$\left\{\begin{array}{c}{\mathrm{\ϵ}}_l+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_{X_{i,j}^f}=1\\ {\mathrm{\ϵ}}_l=0.2+\left(\frac{0.5}{500*{0.25}^{\left({10}^3/\left({10}^6\·z+1\right)\right)}+1}\right)\\ {\mathrm{\ϵ}}_{X_{i,j}^f}=\frac{X_{i,j}^f}{{\mathrm{\ρ}}_i^f}\end{array}\right.$

In above formula, j represents biomembranous hierarchy number, and n represents biomembranous total hierarchy number, and z represents biomembrane depth of seam division, △Z_jRepresent 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, S_FAerobic storage, S_FAnaerobism storage, aerobic growth, anaerobism raw Long, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_AAerobic respiration, X_AAnaerobic respiration, SRB growth, SRB decline, heterotrophic bacteria inactivation, X_STOPassivation, 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 phase_FAerobic storage, S_FAnaerobism storage, X_HGood Oxide growth, X_HAnaerobic growth, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, nitrification are made With, X_AAerobic respiration and X_AAnaerobic 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 H₂S 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, S_FAerobic storage, S_FAnaerobism storage, X_HAerobic growth, X_HAnaerobic growth, X_HAerobic respiration, X_HAnaerobism exhale Suction, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_AAerobic respiration and X_AThe 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:

$\left\{\begin{array}{c}{\mathrm{\ϵ}}_l+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_{X_{i,j}^f}=1\\ {\mathrm{\ϵ}}_l=0.2+\left(\frac{0.5}{500*{0.25}^{\left({10}^3/\left({10}^6\·z+1\right)\right)}+1}\right)\\ {\mathrm{\ϵ}}_{X_{i,j}^f}=\frac{X_{i,j}^f}{{\mathrm{\ρ}}_i^f}\end{array}\right.$

In above formula, j represents biomembranous hierarchy number, and n represents biomembranous total hierarchy number, and z represents biomembrane depth of seam division, △Z_jRepresent 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, S_FAerobic storage, S_FAnaerobism storage, aerobic growth, anaerobism raw Long, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_AAerobic respiration, X_AAnaerobic respiration, SRB growth, SRB decline, heterotrophic bacteria inactivation, X_STOPassivation, 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 phase_FAerobic storage, S_FAnaerobism storage, X_HGood Oxide growth, X_HAnaerobic growth, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, nitrification are made With, X_AAerobic respiration and X_AAnaerobic 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 simultaneously₂S'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, SO₄ ^2-、DS、TS、 TCOD、SCOD、VSS、NH₄-N, NOx-N and H₂The 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 S_O、S_I、S_F、S_NH、N₂、S_NO、S_ALK；X is included with graininess component_I、X_S、X_H、X_STO、X_A、X_TS.And increase component X_INA、 X_SRB、S_SO4And S_H2SDeng, 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 hours₄ ^2-、DS、TS、TCOD、SCOD、TSS、VSS、NH₃-N、NO_xThe 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 MH17₂S (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, S_FAerobic storage, S_FAnaerobism storage, X_HAerobic growth, X_HAnaerobic growth, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_AAerobic respiration and X_AAnaerobic 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, i_NBMRepresent the nitrogen content in microorganism, i_NXIRepresent X_IIn 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:b_SRBX_SRB

Formula 15:

Formula 16:

Formula 17:

Formula 18:

Formula 19:

Formula 20:k_SO·S_O ^0.1·S_H2S

Formula 21:

Formula 22:k_La·(S_O,sat-S_O)

X in table 4₁,x₂,…x₂₀And y₁,y₂,…y₂₁Can draw according to COD and N conservation equation in each course of reaction:

$\underset{i=1}{\overset{21}{\Σ}}{v}_{k,i}\·l_{1,i}=0$

$\underset{i=1}{\overset{21}{\Σ}}{v}_{k,i}\·l_{2,i}=0$

As sought x1 and y1 of hydrolysis:

Only see COD conservation:

$\underset{i=1}{\overset{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:

$\underset{i=1}{\overset{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_{NSS}+y_1\·1+(-1\·i_{NBM})=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:

$k1=\frac{\[{\mathrm{HS}}^{-}\]\[H^{+}\]}{\[H_{2}S\]},k2=\frac{\[S^{2-}\]\[H^{+}\]}{\[HS\]}$

Wherein, k1 is H₂The 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；[H₂S] for not dissociating hydrogen sulfide molecule (H₂S) concentration；K2 is H₂Two grades of dissociation constants of S； [S^2-] for bearing the sulphion (S of bivalence^2-) 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 H₂S pollutes essential and plays a major role.Based on biomembrane for blow-off pipe H₂The 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 △ Z_j, every class the dissolved material S in each layering_i,jAnd the every class granulating material X in each layering_i,jAfter, Calculate in conjunction with following formula, and constantly adjust thickness in monolayer △ Z_j, every class the dissolved material S in each layering_i,jAnd it is each Every class granulating material X in layering_i,j, until obtaining dynamic biological membrane modle after meeting following formula:

$\left\{\begin{array}{c}{\mathrm{\ϵ}}_l+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_{X_{i,j}^f}=1\\ {\mathrm{\ϵ}}_l=0.2+\left(\frac{0.5}{500*{0.25}^{\left({10}^3/\left({10}^6\·z+1\right)\right)}+1}\right)\\ {\mathrm{\ϵ}}_{X_{i,j}^f}=\frac{X_{i,j}^f}{{\mathrm{\ρ}}_i^f}\end{array}\right.$

In above formula, j represents biomembranous hierarchy number, and n represents biomembranous total hierarchy number, and z represents biomembrane depth of seam division, △Z_jRepresent 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 △ Z_jAdjusting range be 0.5 × 10^-4～1.5 × 10^-4M, if △ is Z_jLess than 0.5*10^-4 Then it is incorporated to last layer, if △ is Z_jMore 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_{att,X_i^w}=k_{att}\·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^-1；Referring to the concentration of granulating material Xi in water body, unit is gCODm^-3, wherein subscript w represents water body；k_attIt 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:

$\left\{\begin{array}{c}{r}_{\det ,X_i}=k_{\det }{\mathrm{\μ}}_H{L_f}^2{\left(\frac{{\mathrm{\τ}}_w-{\mathrm{\τ}}_w^{\min }}{{\mathrm{\τ}}_w^{\min }}\right)}^{2.5}\tanh \left(\frac{M_i}{M_{TSS}}\right)\\ {\mathrm{\τ}}_w={\mathrm{\ρ}}_w\·g\·S_f\·R_h\\ r_{\det ,X_{i,j}}=-r_{\det ,X_i}\·X_{i,j}/\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{X}_{i,j}\end{array}\right..$

In above formula,It is particulate matter X_iPeel rate, unit is gCODm^-3d^-1；k_detBeing to peel off constant, unit is gm^-5；μ_HBeing the maximum growth rate of heterotrophic microorganism, unit is d^-1；L_fBeing biofilm thickness, unit is m；τ_wIt is shearing force, Unit is N m^-2；M_iIt is particulate matter X_iTS mass；And M_TSIt it is then monoblock biomembranous TS mass；S_fFriction gradient；R_hFor water Power radius, unit is m；ρ_wBeing the mass density of water, unit is gm^-3；And g is acceleration of gravity, unit is m²s^-1；It is Particulate matter X in biomembrane j layering_iPeel rate, unit is gCODm^-3d^-1；X_i,jIt is that in biomembrane j layering, particulate form is organic Thing X_iConcentration.

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:

$\left\{\begin{array}{c}{J}_{i,j}=-D_{f,i,j}\frac{\∂S_{i,j}}{\∂z_j}\\ D_f=D_w\·f_D\\ f_D=\left(1-\frac{0.43{\mathrm{\ρ}}_m^{0.92}}{11.19+0.27{\mathrm{\ρ}}_m^{0.99}}\right)\\ D_{f,i,j}=D_{w,i}\·\left(1-\frac{0.43{\mathrm{\ρ}}_{m,j}^{0.92}}{11.19+0.27{\mathrm{\ρ}}_{m,j}^{0.99}}\right)\end{array}\right.$

In above formula, J_i,jIt is the dissolved material S_iThe diffusion flux unit passing through biomembrane j layered boundary is gm^-2s^-1；S_i,j It is the dissolved material S_iIn the concentration of biomembrane j layer, unit is gm^-3；D_f,i,jIt is the dissolved material S_iExpansion in biomembrane j layer Dissipating coefficient, unit is m²s^-1；D_wBeing the diffusion coefficient in water body, unit is m²s^-1, and f_DIt 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:

D_w,O2=4.864 × 10^-13T²+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, S_FAerobic storage, S_FAnaerobism storage, aerobic growth, anaerobic growth, X_H's Aerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_AAerobic respiration, X_AAnaerobism Breathe, SRB grows, SRB becomes feeble and die, heterotrophic bacteria inactivation, X_STOPassivation, autotrophic bacteria inactivation, SRB inactivation, sulfide-oxidation, fermentation and Reoxygenation totally 21 microbial reaction formulas, and 13 kinds of components include S_O、S_F、S_I、S_NH、S_NO、S_SO4、S_H2S、X_I、X_S、X_H、X_STO、X_AWith X_SRB, 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, i_NBMIt is the nitrogen content in microorganism, i_NXIIt is X_IIn nitrogen content.

X in table 5₁,x₄,x₅And y₁,y₂,y₃Can draw according to COD conservation equation in each course of reaction following:

$\underset{i=1}{\overset{5}{\Σ}}{v}_{k,i}\·l_{1,i}=0$

$\underset{i=1}{\overset{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 certain_20℃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)₂S is diffused in gas phase, when it reaches balance at liquid-gas interface, uses Henry's law and unit conversion draw H in gas phase₂S pollutant levels, computing formula is as follows:

$\left\{\begin{array}{c}{P}^{*}=C/H\\ P^{*}V=nRT\\ Cg=M/V\end{array}\right.$

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 m³；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 Z_jLess than 0.5*10^-4Then it is incorporated to last layer；If △ is Z_jMore 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 K_SRB,O.Root It is predicted the H in the gentle phase of total sulfur compound in the sewage of model₂S simulation result contrasts with measured value, by model parameter μ_SRBIt is corrected to 5.6d^-1, K_SRB,OIt is corrected to 1.0gO₂m^-3, simulation result is the most identical with measured data, such as Fig. 5 and 6.Prediction mould Type can emulate sulfide and H₂The 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 MH17₂S gas meansigma methods is 12.9ppm, and the meansigma methods of simulation curve is 12.5ppm.Mould Type can simulate S^2-And H₂Outside the concentration change of S, additionally it is possible to prediction H₂The 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 H₂The 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 H₂S 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 H₂The 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%, H₂The 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 H₂The mean concentration of S.

3. prioritization scheme three, the nitrate dosage of increase by 100%, when dosage doubles, average H₂S concentration from 6.6 are reduced to 1.1ppm, reduce 88%, such as Figure 10.H in most cases₂The concentration of S reaches null value, but due to flow Change, however it remains H₂The situation of S peak value.Obviously result shows that dosage doubles, and still there will be certain time H₂S 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 H₂S concentration is minimized, and adds the most The H produced₂S concentration is the lowest, but for suppression H₂S spillover effect is little.Relatively control H₂S 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 H₂S gas is controlled effect emulation, even can Water quality such as COD, ammonia nitrogen, DO, SO to blow-off pipe₄ ^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 phase_FAerobic storage, S_FAnaerobism storage, X_HAerobic life Long, X_HAnaerobic growth, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_A Aerobic respiration and X_AThe 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:

$\left\{\begin{array}{c}{\mathrm{\ϵ}}_l+\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\ϵ}}_{X_{i,j}^f}=1\\ {\mathrm{\ϵ}}_l=0.2+\left(\frac{0.5}{500*{0.25}^{({10}^3/({10}^6\·z+1\left)\right)}+1}\right)\\ {\mathrm{\ϵ}}_{X_{i,j}^f}=\frac{X_{i,j}^f}{{\mathrm{\ρ}}_i^f}\end{array}\right.$

In above formula, j represents biomembranous hierarchy number, and n represents biomembranous total hierarchy number, and z represents biomembrane depth of seam division, △ Z_j 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, S_FAerobic storage, S_FAnaerobism storage, aerobic growth, anaerobic growth, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_AAerobic respiration, X_A's Anaerobic respiration, SRB growth, SRB decline, heterotrophic bacteria inactivation, X_STOPassivation, 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 phase_FAerobic storage, S_FAnaerobism storage, X_HAerobic life Long, X_HAnaerobic growth, X_HAerobic respiration, X_HAnaerobic respiration, X_STOAerobic respiration, X_STOAnaerobic respiration, Nitrification, X_A Aerobic respiration and X_AAnaerobic 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.