METHOD FOR REDUCING CONCENTRATION OF MICROORGANISM-DERIVED DISSOLVED ORGANIC NITROGEN IN WASTEWATER

Abstract

A method for reducing mDON concentration in wastewater, including a) acquiring a kinetics associated with production and consumption of a mDON of an activated sludge system, and importing a kinetic expression of the mDON into a conventional activated sludge model No. 1 (ASM1) to build a kinetic equation for the mDON; b) inputting component variables, parameter variables, model matrices, process rate equation and operating parameters of a predictive model into a simulation software to build an ASM-mDON model; c) inputting initial values of the component variables and the parameter variables into the simulation software for model initialization; d) acquiring initial mDON kinetic and sensitivity analysis results, selecting corresponding parameters, calibrating kinetic and stoichiometric parameters of the ASM-mDON model using a parameter estimation function of the simulation software; and e) replacing the initial values of the ASM-mDON model with optimal values obtained in d).

Claims

1. A method for reducing mDON concentration in wastewater, the method comprising: a) integrating, by a data processor, a kinetic expression for microbial-derived dissolved organic nitrogen (mDON) into a conventional activated sludge model No. 1 (ASM1) to build an mDON predictive model; b) inputting into software, by the data processor, component variables, parameter variables, model matrices, process rate equations, and operational parameters of the mDON predictive model to build an ASM-mDON predictive model; c) inputting initial values, by the data processor, for the component variables and the parameter variables to initialize the ASM-mDON predictive model in software; d) running, by the data processor, a preliminary simulation of mDON kinetics; conducting, by the data processor, a sensitivity analysis; and selecting corresponding parameters that impact the predictions of the ASM-mDON predictive model based on the results from the preliminary simulation and the sensitivity analysis, e) calibrating, by parameter estimation feature in software, via the data processor, dynamic parameters and stoichiometric parameters in the ASM-mDON predictive model to obtain optimal values, thereby enhancing the accuracy of mDON concentration prediction; f) replacing, by the data processor, the initial values for the dynamic parameters and the stoichiometric parameters in the ASM-mDON predictive model with the optimal values obtained from the calibration to optimize the ASM-mDON predictive model; g) connecting a programmable logic controller via data signals to a dissolved oxygen sensor, a level sensor, a pH sensor, a plurality of actuators, and the data processor; collecting, by the programmable logic controller, the component variables and the parameter variables from the dissolved oxygen sensor, the level sensor, and the pH sensor; and transmitting, by the programmable logic controller, the collected data to the data processor; h) inputting, by the data processor, the collected component variables and the parameter variables into the ASM-mDON predictive model, and outputting an mDON concentration; i) inputting variations in operational parameters of an activated sludge system into the ASM-mDON predictive model to determine the changes in mDON concentration under different operational conditions; and identifying an optimal combination of the operational parameters that effectively reduce mDON production based on the model predictions; and j) returning, by the data processor, the optimal combination of the operational parameters to the programmable logic controller; adjusting, by the programmable logic controller, the operational parameters of the plurality of actuators, thereby ensuring effective reduction in mDON concentration during wastewater treatment.

2. The method of claim 1, wherein the activated sludge system comprises a fully mixed steady state activated sludge; the activated sludge has a sludge age of 5-30 days, and a concentration of 2000-5000 mg/L.

3. The method of claim 1, wherein the ASM-mDON model is used for study of the mDON released by microorganisms in the activated sludge system, and the model comprises: seven components: heterotrophic bacteria X.sub.H, autotrophic bacteria X.sub.A, inert particles X.sub.I, nitrate nitrogen S.sub.NO, ammonia nitrogen S.sub.NH, microorganism-derived dissolved organic nitrogen S.sub.DON, dissolved oxygen S.sub.O; five reaction processes: a growth process and an endogenous respiration process of heterotrophic bacteria using ammonium chloride as a substrate; a growth process and an endogenous respiration process of autotrophic bacteria using ammonium chloride as a substrate; and an ammonization process of mDON; and eighteen parameters: maximum specific growth rate {circumflex over ()}.sub.H of heterotrophic bacteria, yield coefficient Y.sub.H of heterotrophic bacteria, attenuation coefficient b.sub.H of heterotrophic bacteria, half-saturation constant K.sub.H,NH for ammonia nitrogen of heterotrophic bacteria, half-saturation constant K.sub.H,O for dissolved oxygen of heterotrophic bacteria, maximum specific growth rate {circumflex over ()}.sub.A of autotrophic bacteria, substrate utilization ratio f.sub.H,DON of heterotrophic bacteria converting the substrate into the mDON, yield coefficient Y.sub.A of autotrophic bacteria, attenuation coefficient b.sub.A of autotrophic bacteria, half-saturation constant K.sub.A,NH for ammonia nitrogen of autotrophic bacteria, half-saturation constant K.sub.A,O for dissolved oxygen of autotrophic bacteria, substrate utilization ratio f.sub.A,DON of autotrophic bacteria converting the substrate into the mDON, proportion of nitrogen i.sub.XB in an organism, proportion of nitrogen i.sub.XP in the product of the organism, substrate utilization ratio f.sub.NO of autotrophic bacteria converting the substrate into the nitrate nitrogen, proportion of inert particles f.sub.1 yielded in the organism, ammonification rate k.sub.a, and half-saturation constant K.sub.H,DON for mDON.

4. The method of claim 3, wherein change rates of the seven components of the ASM-mDON model satisfy with the following formulas: $\begin{array}{cc}{X}_H\text{:}\frac{{\mathrm{dX}}_H}{dt}={\stackrel{}{}}_H{M}_{H,\mathrm{NH}}\left(t\right)M_{H,O}\left(t\right)X_H\left(t\right)-b_H{M}_{H,O}\left(t\right)X_H\left(t\right)& \left(1\right)\end{array}$ $\begin{array}{cc}{X}_A\text{:}\frac{{\mathrm{dX}}_A}{\mathrm{dt}}={\stackrel{}{}}_A{M}_{A,\mathrm{NH}}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)-b_A{M}_{A,O}\left(t\right)X_A\left(t\right)& \left(2\right)\end{array}$ $\begin{array}{cc}{S}_{\mathrm{NH}}\text{:}\frac{d{S}_{NH}}{dt}=-\left(\frac{f_{H,\mathrm{DON}}}{Y_H}+i_{\mathrm{XB}}\right){\stackrel{}{}}_H{M}_{H,\mathrm{NH}}\left(t\right)M_{H,O}\left(t\right)X_H\left(t\right)-\left(\frac{f_{A,\mathrm{DON}}+f_{\mathrm{NO}}}{Y_A}+i_{\mathrm{XB}}\right){\stackrel{}{}}_A{M}_{A,NH}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)+k_a{M}_{H,\mathrm{DON}}\left(t\right)X_H\left(t\right)& \left(3\right)\end{array}$ $\begin{array}{cc}{S}_{\mathrm{DON}}\text{:}\frac{d{S}_{DON}}{dt}=\frac{f_{H,\mathrm{DON}}}{Y_H}{\stackrel{}{}}_H{M}_{H,\mathrm{NH}}\left(t\right)M_{H,O}\left(t\right)X_H\left(t\right)+\frac{f_{A,\mathrm{DON}}}{Y_A}{\stackrel{}{}}_A{M}_{A,NH}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)-k_a{M}_{H,\mathrm{DON}}\left(t\right)X_H\left(t\right)& \left(4\right)\end{array}$ $\begin{array}{cc}{S}_{\mathrm{NO}}\text{:}\frac{{\mathrm{dS}}_{\mathrm{NO}}}{\mathrm{dt}}=\frac{f_{\mathrm{NO}}}{Y_A}{\stackrel{}{}}_A{M}_{A,\mathrm{NH}}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)& \left(5\right)\end{array}$ $\begin{array}{cc}{X}_I:\frac{{\mathrm{dX}}_I}{\mathrm{dt}}=f_I{b}_H{M}_{H,O}\left(t\right)X_H\left(t\right)+f_I{b}_A{M}_{A,O}\left(t\right)X_A\left(t\right)& \left(6\right)\end{array}$ $\begin{array}{cc}{S}_O:\frac{{\mathrm{dS}}_O}{\mathrm{dt}}=k_L\left(S_O^{*}-S_O\right)-\left(1-\frac{2.86f_{H,\mathrm{DON}}}{Y_H}\right){\stackrel{}{}}_H{M}_{H,\mathrm{NH}}\left(t\right)M_{H,O}\left(t\right)X_H\left(t\right)-\left(1-\frac{2\text{.86}{f}_{A,\mathrm{DON}}}{Y_A}-\frac{4\text{.57}{f}_{\mathrm{NO}}}{Y_A}\right){\stackrel{}{}}_A{M}_{A,\mathrm{NH}}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)+\left(i_{\mathrm{XB}}-f_I{i}_{\mathrm{XP}}\right)b_H{M}_{H,O}\left(t\right)X_H\left(t\right)+\left(i_{\mathrm{XB}}-f_I{i}_{\mathrm{XP}}\right)b_A{M}_{A,O}\left(t\right)X_A\left(t\right)& \left(7\right)\end{array}$ M.sub.H,NH (t) is a Monod term determined by the substrate for the heterotrophic bacteria; M.sub.A,NH (t) is a Monod term determined by the substrate for the autotrophic bacteria; M.sub.H,O (t) is a Monod term determined by the dissolved oxygen for the heterotrophic bacteria; M.sub.A,O (t) is a Monod term determined by the dissolved oxygen for the autotrophic bacteria; M.sub.H,DON(t) is a Monod term determined by the mDON in the heterotrophic bacteria; k.sub.L is an exchange rate between a gas phase and a liquid phase; and S.sub.O* is a maximum solubility of oxygen.

5. The method of claim 3, wherein the mDON in wastewater is calculated using the following kinetic equation: $\begin{array}{cc}\frac{{\mathrm{dS}}_{DON}}{\mathrm{dt}}=\frac{f_{H,\mathrm{DON}}}{Y_H}{\stackrel{}{}}_H{M}_{H,\mathrm{NH}}\left(t\right)M_{H,O}\left(t\right)X_H\left(t\right)+\frac{f_{A,\mathrm{DON}}}{Y_A}{\stackrel{}{}}_A{M}_{A,\mathrm{NH}}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)-k_a{M}_{H,\mathrm{DON}}\left(t\right)X_H\left(t\right).& \left(8\right)\end{array}$

6. The method of claim 3, wherein a single-step size of the AMS-mDON model is 0.1, and a total response time for the predictive model is a product of a calculation capacity and the single-step size.

7. The method of claim 1, wherein in (g), collecting the component variables and the parameter variables comprises filtering an influent sample from a wastewater treatment plant using a membrane filter; measuring chemical oxygen demand (COD), concentrations of total nitrogen, nitrate nitrogen, nitrite nitrogen, ammonia nitrogen, and dissolved organic nitrogen of the influent sample filtered, respectively; and measuring yield coefficient Y.sub.H of heterotrophic bacteria, attenuation coefficient b.sub.H heterotrophic bacteria, and maximum specific growth rate {circumflex over ()}.sub.H of heterotrophic bacteria for the activated sludge.

8. The method of claim 7, wherein the wastewater treatment plant operates at an ambient temperature ranging from 15 to 25 C., and an influent pH thereof is 6.0-8.0.

9. The method of claim 7, wherein the concentration of the dissolved organic nitrogen is a difference between concentrations of total nitrogen and ammonia nitrogen, nitrate nitrogen and nitrite nitrogen; the concentration of the total nitrogen is measured by using potassium persulfate oxidation-ion chromatography, or potassium persulfate oxidation-ultraviolet spectrophotometry; the concentration of the ammonia nitrogen is measured by using salicylic acid-hypochlorite spectrophotometry; the concentration of the nitrate nitrogen is measured by using the ion chromatography or ultraviolet-visible spectrophotometry; the concentration of the nitrite nitrogen is measured by using ion chromatography or N-(1-naphthyl)-ethylenediamine spectrophotometry; and the COD is measured by using potassium dichromate method or rapid digestion method.

10. The method of claim 1, wherein in (i), the variations in operational parameters satisfy the following conditions: hydraulic retention time (HRT) between 3 hours and 12 hours, activated sludge age between 5 days and 30 days, activated sludge concentration between 2000 mg/L and 5000 mg/L, carbon-to-nitrogen ratio (C/N) greater than 6, and daily nitrogen load per kilogram of activated sludge less than 0.07 kg TN.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0036] FIG. 1 is a graft showing the predictive result of the concentration of the mDON according to one embodiment of the disclosure; and

[0037] FIG. 2 is a graft showing the predictive result of the concentration of the mDON according to verification embodiment of the disclosure.

DETAILED DESCRIPTION

[0038] To further illustrate the disclosure, embodiments detailing a method for reducing mDON concentration in wastewater are described below. It should be noted that the following embodiments are intended to describe and not to limit the disclosure.

Example 1

[0039] The example was a simulation of the operation of a laboratory-scale sequencing batch reactor (SBR) for the treatment of activated sludge. As a raw material, the wastewater containing particular compositions (excluding dissolved organic nitrogen) was prepared to support the growth of microorganisms in the activated sludge. The prepared wastewater contained the following compositions: 30030 mg/L COD, 205 mg/L total nitrogen, and 3.50.5 mg/L total phosphorus. The operating parameters of the sequencing batch reactor: the effective volume of 2 L, the operating cycle of 6 h, and the hydraulic retention time of 12 h, and the activated sludge age of 20 d. The operating mode of the sequencing batch reactor at ambient temperature of 25 C.: the inlet valve opens and the influent was filled in, followed by mixing and aeration for 300 min. The mixed liquor was sedimented for 50 min and the supernatant was drained out of the sequencing batch reactor. The activated sludge concentration is within the range of 2000200 mg/L, and a pH value of 7.50.5. The term hydraulic retention time, as used herein, refers to the average time that wastewater remains in the SBR. The term activated sludge age, as used herein, refers to the average time a particle of suspended solids remains in an activated sludge system. The term activated sludge concentration, as used herein, refers to the concentration of suspended solids in the activated sludge system. The term daily nitrogen load, as used herein, refers to the amount of nitrogen that is introduced into the activated sludge system. Within the wastewater treatment system, a programmable logic controller is connected via data signals to a dissolved oxygen sensor, a level sensor, a pH sensor, a plurality of actuators, and a data processor. The dissolved oxygen sensor is configured to measure the concentration of dissolved oxygen in the prepared wastewater. The pH sensor is configured to measure the pH level of the activated sludge in the SBR. The level sensor is configured to continuously monitor the liquid level within various stages of the wastewater treatment process. The plurality of actuators are configured to control various processes such as flow, level, aeration, and sludge handling. The data processor is configured to build an ASM-mDON model and analyzes the data transmitted from the programmable logic controller. The plurality of actuators include, but are not limited to, an influent pump, an effluent pump, a sludge discharge valve, an aeration blower, and a sludge recirculation pump. The programmable logic controller collects the component variables and the parameter variables from the dissolved oxygen sensor, the level sensor, and the pH sensor, and then transmits the collected data to the data processor. The data processor analyzes the component variables and the parameter variables and returns an optimal combination of the operational parameters to the programmable logic controller. Upon receiving the optimal operational variables, the programmable logic controller manages the plurality of actuators by sending specific commands to start, stop, or adjust their operating levels, thereby ensuring effective reduction in mDON concentration during wastewater treatment. A method for reducing mDON concentration in wastewater comprises:

1. Building the ASM-mDON Model:

[0040] The inputs to the ASM-mDON model comprise: X.sub.H (concentration of heterotrophic bacteria), X.sub.A (concentration of autotrophic bacteria), X.sub.I(concentration of inert particles), S.sub.NO(concentration of nitrate nitrogen), S.sub.NH (concentration of ammonia nitrogen), S.sub.DON(concentration of the mDON), S.sub.O(concentration of dissolved oxygen), which were state variables; the model matrix corresponding to the process rate equation for the components were inputted into the reaction process included in the software, thereby building a simulation of the sequencing batch reactor for the treatment of activated sludge.

[0041] The simulation model was the simplified ASM-mDON model as shown in Table 1:

TABLE-US-00001 TABLE 1 Process rate equations for components Components Process rate equation No. X.sub.H [00003]$\frac{{\mathrm{dX}}_H}{\mathrm{dt}}={\hat{}}_H{M}_{H,\mathrm{NH}}\left(t\right)M_{H,O}\left(t\right)X_H\left(t\right)-b_H{M}_{H,O}\left(t\right)X_H(t($ (1) X.sub.A [00004]$\frac{{\mathrm{dX}}_A}{\mathrm{dt}}={\hat{}}_A{M}_{A,\mathrm{NH}}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)-b_A{M}_{A,O}\left(t\right)X_A\left(t\right)$ (2) S.sub.NH [00005]$\frac{{\mathrm{dS}}_{\mathrm{NH}}}{\mathrm{dt}}=-\left(\frac{f_{H,\mathrm{DON}}}{Y_H}+i_{\mathrm{XB}}\right){\hat{}}_H{M}_{H,\mathrm{NH}}\left(t\right)M_{H,O}\left(t\right)X_H\left(t\right)-\left(\frac{f_{A,\mathrm{DON}}+f_{\mathrm{NO}}}{Y_A}+i_{\mathrm{XB}}\right){\hat{}}_A{M}_{A,\mathrm{NH}}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)+k_a{M}_{H,\mathrm{DON}}\left(t\right)X_H\left(t\right)$ (3) S.sub.DON [00006]$\frac{{\mathrm{dS}}_{\mathrm{DON}}}{\mathrm{dt}}=\frac{f_{H,\mathrm{DON}}}{Y_H}{\hat{}}_H{M}_{H,\mathrm{NH}}\left(t\right)M_{H,O}\left(t\right)X_H\left(t\right)+\frac{f_{A,\mathrm{DON}}}{Y_A}{\hat{}}_A{M}_{A,\mathrm{NH}}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)-k_a{M}_{H,\mathrm{DON}}\left(t\right)X_H\left(t\right)$ (4) S.sub.NO [00007]$\frac{{\mathrm{dS}}_{\mathrm{NO}}}{\mathrm{dt}}=\frac{f_{\mathrm{NO}}}{Y_A}{\hat{}}_A{M}_{A,\mathrm{NH}}\left(t\right)M_{A,O}\left(t\right)X_A\left(t\right)$ (5) X.sub.I [00008]$\frac{{\mathrm{dX}}_I}{\mathrm{dt}}=f_I{b}_H{M}_{H,O}\left(t\right)X_H\left(t\right)+f_I{b}_A{M}_{A,O}\left(t\right)X_A\left(t\right)$ (6)

[0042] where M.sub.H,NH(t) is the Monod term determined by the substrate for the heterotrophic bacteria; M.sub.A,NH (t) is the Monod term determined by the substrate for the autotrophic bacteria; M.sub.H,O (t) is the Monod term determined by the dissolved oxygen for the heterotrophic bacteria; M.sub.A,O (t) is the Monod term determined by the dissolved oxygen for the autotrophic bacteria; M.sub.H,DON (t) is the Monod term determined by the mDON in the heterotrophic bacteria; k.sub.L is the exchange rate between the gas phase and the liquid phase; S.sub.O* is the maximum solubility of oxygen.

[0043] The model matrix corresponding to the process rate equation for the components were shown in Table 2:

TABLE-US-00002 TABLE 2 Matrix imported into reaction process included in software Reaction process X.sub.H X.sub.A S.sub.NH S.sub.DON S.sub.NO X.sub.I S.sub.O Process rate equation Growth of heterotrophic 1 [00009]$-\frac{f_{H,\mathrm{DON}}}{Y_H}-i_{\mathrm{XB}}$ [00010]$\frac{f_{H,\mathrm{DON}}}{Y_H}$ [00011]$1-\frac{2\text{.86}{f}_{H,\mathrm{DON}}}{Y_H}$ {circumflex over ()}.sub.H M.sub.H, NH(t) M.sub.H,O(t)X.sub.H(t) bacteria Growth of autotrophic 1 [00012]$-\frac{f_{A,\mathrm{DON}}+f_I}{Y_A}-i_{\mathrm{XB}}$ [00013]$\frac{f_{A,\mathrm{DON}}}{Y_A}$ [00014]$\frac{f_{\mathrm{NO}}}{Y_A}$ [00015]$1-\frac{2\text{.86}{f}_{A,\mathrm{DON}}}{Y_A}-\frac{4\text{.57}{f}_{\mathrm{NO}}}{Y_A}$ {circumflex over ()}.sub.AM.sub.A, NH(t)M.sub.A, O(t)X.sub.A(t) bacteria Endogenous 1 f.sub.I b.sub.HM.sub.H, O(t)X.sub.H(t) respiration of heterotrophic bacteria Endogenous 1 f.sub.I b.sub.AM.sub.A, O(t)X.sub.A(t) respiration of autotrophic bacteria Ammonization 1 1 k.sub.aM.sub.H,DON(t)X.sub.H(t)

2. Determining the Components of Influent and the Values of Model Parameters of the ASM-mDON Model:

[0044] 50 mL of influent was sampled directly from the sequencing batch reactor, and filtered using a cellulose acetate membrane filter having pore size of 0.45 m. The concentration of total nitrogen was 20 mg/L measured by using potassium persulfate oxidation-ion chromatography, or potassium persulfate oxidation-ultraviolet spectrophotometry; the concentration of ammonia nitrogen was 20 mg/L measured by using salicylic acid-hypochlorite spectrophotometry; the nitrate nitrogen was 0 mg/L measured by using the ion chromatography or ultraviolet-visible spectrophotometry; the nitrite nitrogen was 0 mg/L measured by using ion chromatography or N-(1-naphthyl)-ethylenediamine spectrophotometry; and the COD was 300 mg/L measured by using potassium dichromate method or rapid digestion method. The concentration of the dissolved organic nitrogen was 0 mg/L, that is, the difference between the sum of total nitrogen and ammonia nitrogen and the sum of nitrate nitrogen and nitrite nitrogen.

[0045] The default values for kinetics and stoichiometric parameters of the conventional model, and the water quality parameters determined in 1) were used for simulation of parameters in relation to the mDON yielded in the wastewater treatment process; the simulation parameters were set as follows: {circumflex over ()}.sub.H, 0.8 h.sup.1; Y.sub.H, 0.67 mg(COD)/mg(N); b.sub.H, 0.62 h.sup.1; K.sub.H,NH, 0.05 mg(N)/L; K.sub.H,O, 0.2 mg(N)/L; {circumflex over ()}.sub.A, 0.3 h.sup.1; f.sub.H,DON, 0.04; Y.sub.A, 3.4 mg(COD)/mg(N); b.sub.A, 0.15 h.sup.1; K.sub.A,NH, 5 mg(N)/L; K.sub.A,O, 0.4 mg(N)/L; f.sub.A,DON, 0.04; i.sub.XB, 0.07 mg(N)/mg(COD); i.sub.XP, 0.03 mg(N)/mg(COD); f.sub.NO, 0.8; f.sub.I, 0.2; k.sub.a, 0.04 L/(mg(N) d); K.sub.H,DON, 1.5 mg(N)/L; where H was the maximum specific growth rate of heterotrophic bacteria, Y.sub.H was the yield coefficient of heterotrophic bacteria, b.sub.H was the attenuation coefficient of heterotrophic bacteria, K.sub.H,NH was the half-saturation constant for ammonia nitrogen of heterotrophic bacteria, K.sub.H,O was the half-saturation constant for dissolved oxygen of heterotrophic bacteria, {circumflex over ()}.sub.A was the maximum specific growth rate of autotrophic bacteria, f.sub.H,DON was the substrate utilization ratio of heterotrophic bacteria converting the substrates into the mDON, Y.sub.A was the yield coefficient of autotrophic bacteria, b.sub.A was the attenuation coefficient of autotrophic bacteria, K.sub.A,NH was the half-saturation constant for ammonia nitrogen of autotrophic bacteria, K.sub.A,O was the half-saturation constant for dissolved oxygen of autotrophic bacteria, f.sub.A,DON was the substrate utilization ratio of autotrophic bacteria converting the substrates into the mDON, i.sub.XB was the proportion of nitrogen in the organism, i.sub.XP was the proportion of nitrogen in the product of the organism, f.sub.NO was the substrate utilization ratio of autotrophic bacteria converting the substrates into the nitrate nitrogen, f.sub.I was the rate of inert particles yielded in the organism, k.sub.a was the ammonification rate, K.sub.H,DON was the half-saturation constant for mDON.

3. Predicting the Concentration of mDON in Wastewater:

[0046] The components of influent and values of model parameters determined in 2) were fed into the software for modeling mDON to predict the concentration of the mDON in wastewater; where the single-step size was 0.1, and the calculation capacity was 60 steps, and the simulation process was based on the mDON participating in the biochemical reactions. The model-predicted result was shown in FIG. 1.

Example 2

[0047] The example was the same as Example 1, except for the influent from the municipal wastewater treatment plant A. The operating parameters of the sequencing batch reactor: the influent temperature was 15 C., and hydraulic retention time was 8 h, activated sludge age was 20 d. The influent contains the following compositions: COD 96.2-120.6 mg/L, total nitrogen 23.7-29.1 mg/L, total phosphorus 2.0-3.5 mg/L, pH 7.4-8.0, and inert particles 3000-3200 mg/L. The influent (of greater than 200 mL) in the biological treatment process (i.e. oxidation ditch) and the activated sludge (of greater than 50 mL) were sampled for analysis of the components of the influent, as well as parameter estimation. The influent sample was then filtered using a cellulose acetate membrane filter having pore size of 0.45 m. The COD was measured by using potassium dichromate method or rapid digestion method; the concentration of total nitrogen was measured by using potassium persulfate oxidation-ion chromatography, or potassium persulfate oxidation-ultraviolet spectrophotometry; the concentration of ammonia nitrogen was measured by using salicylic acid-hypochlorite spectrophotometry; the nitrate nitrogen was measured by using the ion chromatography or ultraviolet-visible spectrophotometry; the nitrite nitrogen was measured by using ion chromatography or N-(1-naphthyl)-ethylenediamine spectrophotometry; the concentration of dissolved organic nitrogen was the difference between the sum of total nitrogen and ammonia nitrogen and the sum of nitrate nitrogen and nitrite nitrogen. According to the parameters determined in 1), the initial values of the yield coefficient (Y.sub.H) of heterotrophic bacteria, the attenuation coefficient of heterotrophic bacteria, and the maximum specific growth rate ({circumflex over ()}.sub.H) of heterotrophic bacteria were 0.26 mgCOD/mgN, 0.09 h.sup.1, and 1.0 h.sup.1, respectively.

[0048] Building the ASM-mDON model and finding the optimal parameter values i.sub.XB by parameter estimation: 0.07 mg(N)/mg(COD); k.sub.a, 0.04 L/(mg(N) d); {circumflex over ()}.sub.H, 1.0 h.sup.1; Y.sub.H, 0.30 mg(COD)/mg(N); b.sub.H, 0.05 h.sup.1; K.sub.H,NH, 0.05 mg(N)/L; K.sub.H,O, 0.2 mg(N)/L; {circumflex over ()}.sub.A, 0.3 h.sup.1; f.sub.H,DON, 0.04; Y.sub.A, 3.4 mg(COD)/mg(N); b.sub.A, 0.15 h.sup.1; K.sub.A,NH, 5 mg(N)/L; K.sub.A,O, 0.4 mg(N)/L; f.sub.A,DON, 0.04; i.sub.XP, 0.03 mg(N)/mg(COD); f.sub.NO, 0.8; f.sub.I, 0.2; K.sub.H,DON, 1.5 mg(N)/L.

[0049] Predicting the concentration of the mDON in wastewater: the components of influent and values of model parameters determined in 2) were fed into the software for modeling mDON to predict the concentration of the mDON in wastewater; where the single-step size was 0.1, and the calculation capacity was 240 steps. The model-predicted concentration of the mDON yielded in the oxidation ditch was 2.32 mg/L.

Example 3

[0050] The example was the same as Example 2, except for the influent coming from the municipal wastewater treatment plant A was sampled at different times. The operating parameters of the sequencing batch reactor: the influent temperature was 20 C., and hydraulic retention time was 8 h, activated sludge age was 5 d. The influent contains the following compositions: COD 96.2-120.6 mg/L, total nitrogen 23.7-29.1 mg/L, total phosphorus 2.0-3.5 mg/L, pH 7.4-8.0, and inert particles 3000-3200 mg/L. The influent (of greater than 200 mL) in the biological treatment process (i.e. oxidation ditch) and the activated sludge (of greater than 50 mL) were sampled for analysis of the components of the influent and parameter estimation. The influent sample was then filtered using a cellulose acetate membrane filter having pore size of 0.45 m. COD was measured by using potassium dichromate method or rapid digestion method; the concentration of total nitrogen was measured by using potassium persulfate oxidation-ion chromatography, or potassium persulfate oxidation-ultraviolet spectrophotometry; the concentration of ammonia nitrogen was measured by using salicylic acid-hypochlorite spectrophotometry; the nitrate nitrogen was measured by using the ion chromatography or ultraviolet-visible spectrophotometry; the nitrite nitrogen was measured by using ion chromatography or N-(1-naphthyl)-ethylenediamine spectrophotometry; the concentration of dissolved organic nitrogen was the difference between the sum of total nitrogen and ammonia nitrogen and the sum of nitrate nitrogen and nitrite nitrogen. According to the parameters determined in 1), the initial values of the yield coefficient (Y.sub.H) of heterotrophic bacteria, the attenuation coefficient of heterotrophic bacteria, and the maximum specific growth rate ({circumflex over ()}.sub.H) of heterotrophic bacteria were 0.26 mgCOD/mgN, 0.09 h.sup.1, and 1.0 h.sup.1, respectively.

[0051] Building the ASM-mDON model and finding the optimal parameter values i.sub.XB by parameter estimation: 0.07 mg(N)/mg(COD); k.sub.a, 0.04 L/(mg(N) d); {circumflex over ()}.sub.H, 1.0 h.sup.1; Y.sub.H, 0.30 mg(COD)/mg(N); b.sub.H, 0.05 h.sup.1; K.sub.H,NH, 0.05 mg(N)/L; K.sub.H,O, 0.2 mg(N)/L; {circumflex over ()}.sub.A, 0.3 h.sup.1; f.sub.H,DON, 0.04; Y.sub.A, 3.4 mg(COD)/mg(N); b.sub.A, 0.15 h.sup.1; K.sub.A,NH, 5 mg(N)/L; K.sub.A,O, 0.4 mg(N)/L; f.sub.A,DON, 0.04; i.sub.XP, 0.03 mg(N)/mg(COD); f.sub.NO, 0.8; f.sub.I, 0.2; K.sub.H,DON, 1.5 mg(N)/L.

[0052] Predicting the concentration of mDON in wastewater: the components of influent and values of model parameters determined in 2) were fed into the software for modeling mDON to predict the concentration of the mDON in wastewater; where the single-step size was 0.1, and the calculation capacity was 240 steps. The model-predicted concentration of the mDON yielded in the oxidation ditch was 1.89 mg/L.

Example 4

[0053] The example was the same as Example 1, except for the influent from the municipal wastewater treatment plant B. The operating parameters of the sequencing batch reactor: the influent temperature was 20 C., and hydraulic retention time was 6 h, activated sludge age was 30 days. The influent contained the following compositions: COD 130.9 mg/L, total nitrogen 25.1 mg/L, total phosphorus 5.1 mg/L, pH 7.2, and inert particles 3000-3200 mg/L. The influent (of greater than 200 mL) in the biological treatment process (i.e. oxidation ditch) and the activated sludge (of greater than 50 mL) were sampled for analysis of the components of the influent and parameter estimation. The influent sample was then filtered using a cellulose acetate membrane filter having pore size of 0.45 m. COD was measured by using potassium dichromate method or rapid digestion method; the concentration of total nitrogen was measured by using potassium persulfate oxidation-ion chromatography, or potassium persulfate oxidation-ultraviolet spectrophotometry; the concentration of ammonia nitrogen was measured by using salicylic acid-hypochlorite spectrophotometry; the nitrate nitrogen was measured by using the ion chromatography or ultraviolet-visible spectrophotometry; the nitrite nitrogen was measured by using ion chromatography or N-(1-naphthyl)-ethylenediamine spectrophotometry; the concentration of dissolved organic nitrogen was the difference between the sum of total nitrogen and ammonia nitrogen and the sum of nitrate nitrogen and nitrite nitrogen. According to the parameters determined in 1), the initial values of the yield coefficient (Y.sub.H) of heterotrophic bacteria, the attenuation coefficient of heterotrophic bacteria, and the maximum specific growth rate ({circumflex over ()}.sub.H) of heterotrophic bacteria were 0.2 mgCOD/mgN, 0.05 h.sup.1, and 0.3 h.sup.1, respectively.

[0054] Building the ASM-mDON model and finding the optimal parameter values i.sub.XB by parameter estimation: 0.07 mg(N)/mg(COD); k.sub.a, 0.04 L/(mg(N) d); {circumflex over ()}.sub.H, 0.33 h.sup.1; Y.sub.H, 0.32 mg(COD)/mg(N); b.sub.H, 0.05 h.sup.1; K.sub.H,NH, 0.05 mg(N)/L; K.sub.H,O, 0.2 mg(N)/L; {circumflex over ()}.sub.A, 0.3 h.sup.1; f.sub.H,DON, 0.04; Y.sub.A, 3.0 mg(COD)/mg(N); b.sub.A, 0.15 h.sup.1; K.sub.A,NH, 5 mg(N)/L; K.sub.A,O, 0.4 mg(N)/L; f.sub.A,DON, 0.04; i.sub.XP, 0.03 mg(N)/mg(COD); f.sub.NO, 0.8; f.sub.I, 0.2; K.sub.H,DON, 1.5 mg(N)/L.

[0055] Predicting the concentration of mDON in wastewater: the components of influent and values of model parameters determined in 2) were fed into the software for modeling mDON to predict the concentration of the mDON in wastewater; where the single-step size was 0.1, and the calculation capacity was 60 steps. The model-predicted concentration of the mDON yielded in the oxidation ditch was 4.31 mg/L.

Verification Example

[0056] The Verification Example involves comparing the predicted values of mDON obtained through the ASM-mDON prediction model with actual measurements from a sequencing batch reactor (SBR), using the procedures detailed in Example 1 to ensure accuracy of the predictions of the ASM-mDON prediction model. Water samples were taken every 0.5 hours during a 5-hour reaction period within a single operational cycle of the SBR. The pre-treatment and measurement methods for the water samples followed the procedures detailed in Example 1. As illustrated in FIG. 2, the measured mDON production in the SBR over one cycle closely matched the simulated curve, staying within an acceptable error margin. Additionally, the ASM-mDON prediction model was utilized for further comparative validation. The variations in the operational parameters of the SBR satisfied the following conditions: hydraulic retention time (HRT) between 3 hours and 12 hours, activated sludge age between 5 days and 30 days, activated sludge concentration between 2000 mg/L and 5000 mg/L, carbon-to-nitrogen ratio (C/N) greater than 6, and daily nitrogen load per kilogram of activated sludge less than 0.07 kg TN. Using an orthogonal design, each operational parameter was set at five different levels. The different settings were input into the ASM-mDON prediction model to figure out which combination of the operational parameters produced the lowest concentration of mDON. The optimal combination of the operational parameters for minimizing the concentration of mDON in the SBR was as follows: a HRT of 6.4 hours, an activated sludge age of 20.1 days, an activated sludge concentration of 4333 mg/L, a C/N ratio of 23.6, and a daily nitrogen load of 0.04 kg TN per kilogram of activated sludge. Under the conditions, the concentration of mDON was reduced to 0.52 mg/L. The optimal combination of the operational parameters results in a 30% cost reduction compared to achieving similar mDON reductions through the addition of a carbon source.

[0057] It will be obvious to those skilled in the art that changes and modifications may be made, and therefore, the aim in the appended claims is to cover all such changes and modifications.