環(huán)境數(shù)模課程設(shè)計報告說明書

2016"環(huán)境數(shù)學模型"課程設(shè)計說明書題目活性污泥系統(tǒng)生化反響器中底物降解與微生物增長數(shù)學模型的建立實驗方法與結(jié)果實驗方法工藝流程與反響器本設(shè)計采用的工藝流程如下列圖所示：圖2-1活性污泥系統(tǒng)工藝流程圖本設(shè)計工藝采用活性污泥法處理污水，工藝的主要反響器包括生化反響器和沉淀池。污水通過蠕動泵恒速加到生化反響器中，反響器活性污泥和污水在機械攪拌設(shè)備和鼓風曝氣設(shè)備的共同作用下充分接觸，并在氧氣充足的條件下進展反響。經(jīng)處理后，污泥混液通過管道自流到沉淀池中，在里面實現(xiàn)泥水別離。別離后的水通過溢流堰從周邊排出，直接被排放到下水道系統(tǒng)，沉淀下來的污泥則通過回流泵，全部被抽回進展回流。系統(tǒng)運行過程中，進出水流量、進水質(zhì)量、污水的停留時間、生化反響器的容積、機械攪拌設(shè)備轉(zhuǎn)軸轉(zhuǎn)速、鼓風曝氣裝置的曝氣風量氣速、污泥回流量等參數(shù)在系統(tǒng)運行的過程中都保持不變。待系統(tǒng)持續(xù)運行一周穩(wěn)定后再取樣進展分析。實驗的進水為實驗室配置的污水，污水分別以葡萄糖、尿素、磷酸二氫鉀為碳源、氮源和磷源，其中C：N：P=100：40：1〔濃度比〕，TOC含量為200mg/L。生化反響器污泥混液的容量為12L，污水停留時間為6h。系統(tǒng)運行時間為兩周，第一周是調(diào)適階段，第二周取樣測試，測得的數(shù)據(jù)作為建模的原始數(shù)據(jù)。表2-1污水中各營養(yǎng)物質(zhì)的含量營養(yǎng)物質(zhì)葡萄糖〔C〕尿素〔N〕磷酸二氫鉀〔P〕含量〔mg/L〕500〔200〕85〔40〕8.77〔1〕取樣方法每隔24h取一次樣，通過虹吸管取樣。每次取樣時，先取進水和出水水樣用于測水體的COD指標，其中進水直接取配得的污水溶液，出水取沉淀池上清液。取得的水樣過膜除去水中的懸浮固體和微生物，保存在5ml玻璃消解管中，并在4℃下冷藏保存。取完用于測COD的水樣后，全開污泥回流泵，將沉淀池中的污泥全部抽回生化反響器〔由于實驗裝置的原因，沉淀池排泥管易堵，污泥易積聚在沉淀池中，為更準確測定活性污泥的增長情況，在此實驗中將泥完全抽回后再測定〕，待攪拌均勻后，取5ml污泥混液于干凈、衡重的坩堝中，待用于測污泥混液的SS。分析方法本實驗一共分析進出水COD和污泥混液SS兩個指標。其中COD采用"水質(zhì)快速消解分光光度法"〔HJ/T399-2007〕方法進展分析，SS采用"水質(zhì)懸浮物的測定重量法"〔GB11901-89〕方法進展分析。準確取2ml經(jīng)過膜處理的水樣于5mlcod消解管中，以重鉻酸鉀為氧化劑，硫酸銀-濃硫酸為催化劑，硫酸汞為抗氯離子干擾劑，按一定比例與水樣混合均勻。將消解管放在COD消解儀中，在150℃條件下消解2h。待經(jīng)消解的溶液冷卻后，以空白樣為參比液，在COD分析儀上讀出待測水樣的COD值，記錄數(shù)據(jù)。將裝在已衡重稱重的坩堝中的污泥混液放在烘箱中，在105℃溫度下烘3h以上，保證污泥中的水分被充分除去。坩堝冷卻后衡重稱重，記錄干污泥的質(zhì)量，求得活性污泥的SS。實驗過程的所有樣品都設(shè)置兩個平行樣，最后結(jié)果取平行樣的算術(shù)平均值。實驗結(jié)果實驗數(shù)據(jù)實驗測得數(shù)據(jù)如下表：表2-2活性污泥系統(tǒng)水質(zhì)分析結(jié)果時間MLSS〔mg/L〕進水COD〔mg/L〕出水COD〔mg/L〕消耗COD〔mg/L〕DAAY216104530453DAY3226027375199DAY4258054766482DAY558804140414數(shù)據(jù)分析根據(jù)表中數(shù)據(jù)，混合液懸浮固體MLSS的濃度隨著活性污泥系統(tǒng)運行時間的延長而增大，對于同一活性污泥系統(tǒng)，MLSS濃度與活性污泥量成正比，MLSS的變化在一定程度上反響系統(tǒng)的活性污泥量的變化，系統(tǒng)中的活性污泥量隨時間的延長而呈增長趨勢。在本設(shè)計中，活性污泥系統(tǒng)生化反響器處于穩(wěn)定狀態(tài)，底物充足且濃度根本不變。微生物在營養(yǎng)物質(zhì)充足、曝氣充分的條件下得以快速生長繁殖。實驗測得的結(jié)果與現(xiàn)實規(guī)律一致。本設(shè)計配的營養(yǎng)液中葡萄糖濃度為0.5g/L，對應的理論COD值為533mg/L。實際測得的進水COD值普遍低于533mg/L，這是由于葡萄糖非常容易被微生物降解利用。實驗所用營養(yǎng)液暴露在空氣中，可以被空氣中的微生物利用，使得實際測得值偏低。另一個可能原因是實驗所用的COD消解儀由于使用年限較長，工作性能下降，不能消解充分，進而使得實驗結(jié)果偏低，且實驗結(jié)果偏差較大。進水中碳源葡萄糖為溶解性有機物，容易被微生物利用。進水COD含量不算高，停留時間長達6h，根據(jù)已有的經(jīng)歷，出水COD值應該為0或接近0。但表中DAY3和DAY4的出水COD高達66-75mg/L。這是由于人員操作失誤，在測COD時沒有進展空白樣的參比，由此致使出水COD值比其他時期的要高。數(shù)學建模模型假設(shè)與前提為了簡化實際問題，方便研究問題，在建立底物降解與微生物增長數(shù)學模型時對反響系統(tǒng)作出如下假設(shè)：生化反響器處于完全混合狀態(tài)，所有生物反響均在恒定溫度下進展；進水中的微生物濃度與曝氣池中的活性污泥微生物濃度相比很小，可以忽略；系統(tǒng)處于穩(wěn)定狀態(tài)；二沉池中沒有微生物的活動，沒有污泥累積，泥水別離良好；生化反響器曝氣充足，污泥混液氧含量為常數(shù)，不存在微生物的厭氧代；反響系統(tǒng)的PH維持在適宜微生物生長代的水平，PH為常數(shù)；反響系統(tǒng)的氮含量很低，與以溶解性可生化有機物的代相比，微生物的氨化作用、硝化作用、反硝化作用以及有機氮的水解作用水平很低，可以忽略。不考慮顆粒態(tài)有機物因吸附作用而從反響體系中被除去。顆粒態(tài)有機物質(zhì)的生物網(wǎng)捕瞬間完成。微生物的衰減與電子受體的形式無關(guān)。模型工藝流程與反響器為了使建立的數(shù)學模型更符合實際實驗情況，建立模型采用的工藝流程盡量接近實際情況?；谶@樣的考慮，本設(shè)計采用的模型工藝流程與反響器如下列圖所示：模型處于穩(wěn)定運行的理想狀態(tài)，進出水流量相等。污水從進水管流入生化反響器，經(jīng)生化反響器的活性污泥代處理后，進入沉淀池。污泥和處理水在沉淀池中完全別離，污泥全部回流到生化反響器中，回流流量為進水流量的50%。經(jīng)處理后的凈水從沉淀池中被排出系統(tǒng)。生化反響模型模型的建立是基于活性污泥1號模型〔ActivatedSludgeModelNo.Ⅰ，簡稱ASM1〕建立。ASM1采用了Dold等人1980年提出的死亡—再生理論對單級活性污泥系統(tǒng)的碳氧化、硝化和反硝化三種主要生物學過程中的相關(guān)速率進展了定量描述。模型中采用Monod比生長速率動力學來解釋自養(yǎng)菌或異養(yǎng)菌的生長，與生長速率有關(guān)的單個過程中各組分之間的數(shù)量關(guān)系用化學當量系數(shù)描述。為簡化單位的換算，模型對全部有機組分和生物體統(tǒng)一采用COD當量來表示。ASM1從呼吸過程中電子受體的角度將活性污泥體系劃分為8個過程，包括異養(yǎng)菌的好氧和厭氧生長、自養(yǎng)菌的好氧生長、異養(yǎng)菌和自養(yǎng)菌的衰減、可溶性有機氮的氨化以及網(wǎng)捕性有機物、有機氮的水解。本實驗中，在實驗假設(shè)前提下，生化反響器充分曝氣，不存在厭氧代，不考慮氮磷的影響，故僅存在異養(yǎng)菌的好氧生長、異養(yǎng)菌的衰減和網(wǎng)捕性有機物的水解三種反響過程。另外由于忽略生化反響器中PH的變化和氮對過程微生物生長的影響，本設(shè)計的模型方程式與ASM1存在差異，兩者的動力學方程式比照方下表。表3-1ASM1方程式與本實驗模型方程式比照反響過程ASM1方程式本實驗模型方程式異養(yǎng)菌的好氧生長速率異養(yǎng)菌的衰減速率網(wǎng)捕性有機物的水解速率變量與常數(shù)模型的變量本模型圍繞底物降解對微生物生長的影響建立，過程包括的變量有：溶解性底物濃度SS、異養(yǎng)菌生物量*BH、微生物衰減產(chǎn)生的顆粒性產(chǎn)物含量*p、慢速可生物降解基質(zhì)*S、系統(tǒng)的運行時間t。模型的常數(shù)對于ASM1，模擬的反響過程常數(shù)主要可分為三種類型：化學計量常數(shù)、動力學常數(shù)、反響器常數(shù)。對于本設(shè)計模型，涉及異養(yǎng)菌的好氧生長、異養(yǎng)菌的衰減和網(wǎng)捕性有機物的水解三種過程，設(shè)計的動力學常數(shù)包括：異養(yǎng)菌最大比增長速率μH、異養(yǎng)菌版飽和系數(shù)KS、異養(yǎng)菌的氧半飽和系數(shù)KOH、異養(yǎng)菌衰減系數(shù)bH、最大比水解速率Kh、慢速可生化降解基質(zhì)水解的半飽和系數(shù)K*，化學計量常數(shù)包括：生物量中可轉(zhuǎn)化為顆粒性產(chǎn)物的比例fp、異養(yǎng)菌產(chǎn)率YH，反響器常數(shù)包括反響器的體積V。另外，在實驗假設(shè)前提下，模型在常溫下運行，反響器充分曝氣，生化反響器中污泥混液的溶解氧含量也為定值，各常數(shù)參數(shù)的具體取值如下表。表3-2本設(shè)計模型常數(shù)的取值動力學常數(shù)μH(d-1)KS(mgCOD/L)KOH(mgCOD/L)bH(1/d)Kh(1/d)K*(g/g)6.020.00.100.402.210.15化學計量常數(shù)體積常數(shù)fp{g細胞殘渣/gCOD(生物體)}YH{gCOD(生物體)/gCOD(利用)}生化反響器/L沉淀池/L0.080.671224溶解氧含量So/(mgCOD/L)7.5進水流量Q/(L/d)48計算機實現(xiàn)建模軟件本設(shè)計采用AQUASIM2.0作為建模軟件，AQUASIM是一個被廣泛應用于水系統(tǒng)辨析與過程模擬的電腦程序，不同版本的AQUASIM的根本功能一致。啟動程序后，進入的AQUASIM界面如圖4-1所示。除了一般軟件常有的翻開文件、保存文件等功能外，AQUASIM的主要編輯功能包括：編輯系統(tǒng)、計算、敏感度分析、參數(shù)估值、查看計算結(jié)果、關(guān)閉對話框等，上圖中編號1-6對應的菜單欄按鈕依次對應控制上述六種編輯功能。圖4-1AQUASIM2.0程序界面編輯系統(tǒng)功能選按鈕1編輯系統(tǒng)(editsystem)，AQUASIM界面就會出現(xiàn)如圖4-1的四個窗口：編輯變量(EditVariation)窗口、編輯過程(EditProcesses)窗口、編輯組件(Editpartments)窗口和編輯聯(lián)結(jié)(EditLinks)窗口。這四個窗口編輯的容最終構(gòu)成了整個模型構(gòu)造的主要因素。圖4-2構(gòu)成模型的主要因素之間的關(guān)系圖編輯變量窗口編輯變量窗口用于編輯輸入模型的設(shè)計參數(shù)和設(shè)計變量。變量類型包括狀態(tài)變量(statevarible)、程序變量(programvarible)、常數(shù)變量(constantvarible)、真值數(shù)列(reallistvarible)、變值變量(variablelistvarible)、方程變量(formulavarible)、調(diào)查變量(probevarible)。在建模中用得比擬多的是狀態(tài)變量、常數(shù)變量、真值數(shù)列和方程變量。各種類型變量代表的含義如下表所示。表4-1變量類型及其代表的含義變量類型代表含義statevaribletype用于描述模型過程的狀態(tài)變量programvaribletype用于描述時間、空間坐標、流出量等程序變量constantvaribletype用于描述模型的常數(shù)變量，包括化學計量常數(shù)、動力學常數(shù)等reallistvaribletype用于輸入實測數(shù)據(jù)variablelistvaribletype用于給定的任意參數(shù)值的其他變量間的插，比方多維插值formulavaribletype描述表達方程式的變量probevaribletype用于描述在給定可行組件的可評價變量編輯過程窗口編輯過程窗口主要用于編輯模型反響的過程，包括反響的速率方程式、反響組分的化學計量關(guān)系。編寫反響速率方程所用的變量表達式必須與編輯變量窗口變量的表達式一樣，這樣才能將設(shè)置的變量與速率方程式聯(lián)系起來。化學計量數(shù)的編輯需滿足反響方程中反響組分間的反響關(guān)系。圖4-3編輯過程窗口界面編輯組件窗口編輯組件窗口用于編輯參與反響過程的反響器，如沉淀池、格柵、生化反響器等。通過此窗口，可編輯*一工藝流程組件的變量(variables)、發(fā)生的反響過程(processes)、組件的初始條件(initialdition)和輸入物質(zhì)(input)。除此之外還可設(shè)計反響器的容積。圖4-4編輯組件窗口界面編輯聯(lián)結(jié)窗口編輯聯(lián)結(jié)窗口用于編輯反響組件間的聯(lián)結(jié)關(guān)系，包括物質(zhì)的流向，和流動物質(zhì)的載荷。如圖4-5所示，左邊窗口表示流體方向由反響器流到沉淀池，而右邊窗口表示污泥從沉淀池中回流至反響器，回流的水流量與物質(zhì)量則如圖中輸入框所示。圖4-5編輯聯(lián)結(jié)窗口界面計算功能計算功能的作用在于設(shè)置計算的“步寬〞和“步數(shù)〞。假設(shè)模型運算是以時間為基準，則需設(shè)置相鄰時間的間隔，以及模型一共需要模擬的時長。一般而言，步寬越小模仿的效果就越接近現(xiàn)實，但是由于后面參數(shù)估值和敏感度分析的運算時間隨步寬變小而延長，因而步寬也不宜過小。在每次計算前需要先初始化，去除之前的運算紀錄，再進展正式的運算。敏感度分析功能由于設(shè)計的參數(shù)對模擬效果的影響存在差異，敏感度分析的作用在于通過對設(shè)計的模型預計算，分析選定設(shè)計參數(shù)對模型模擬的影響大小程度。分析結(jié)果最終以.sen文件的形式保存，設(shè)計參數(shù)的敏感度值越大，說明該參數(shù)對模型模擬的影響越大。敏感度分析的結(jié)果給參數(shù)估值時估值參數(shù)的選定提供參考。圖4-6敏感度分析功能窗口界面參數(shù)估值功能參數(shù)估值過程也稱模型的率定。由于通過AQUASIM軟件模擬出來的模型為在規(guī)定的設(shè)計條件下〔參數(shù)規(guī)定、過程規(guī)定、反響器及反響器間的聯(lián)系規(guī)定、運算方式規(guī)定等〕，滿足模型假設(shè)條件而建立起來的的理想模型，模型與實際過程存在一定差異，因此需要進展參數(shù)估值。所謂參數(shù)估值，是指在進展敏感度分析的前提下，或沒有進展敏感度分析，但清楚設(shè)計參數(shù)對模型的敏感程度的前提下，選定幾個敏感度較大設(shè)計參數(shù)，通過指定并激活模型的狀態(tài)變量(statevariable)與實際測得值的計量關(guān)系(fit)，進展計算機的迭代運算，并通過在指定圍修改設(shè)計參數(shù)，使參數(shù)估值后的模型更接近實際情況的過程。經(jīng)率定后的模型能更真實地反映實際情況，但有時可能還不能完全反映真實情況。這時可通過檢查并修改模型的設(shè)計條件，屢次進展模型的計算、敏感度分析和參數(shù)估值，最后才能得到模擬程度較高的設(shè)計模型。圖4-7參數(shù)估值功能窗口界面查看結(jié)果功能在實際操作時，查看結(jié)果功能常與“計算功能〞一起應用。在模型設(shè)計時，編輯好模型的根本因素和運算條件后進展運算，運算結(jié)果通過“查看結(jié)果功能〞查看。使用者可根據(jù)需要，自行設(shè)置需要查看的變量的結(jié)果，也可設(shè)定查看不同反響組件*狀態(tài)變量的運算結(jié)果，運算結(jié)果最終以圖線的形式表現(xiàn)。圖4-8查看計算結(jié)果功能界面關(guān)閉對話框功能關(guān)閉對話框功能用于關(guān)閉編輯系統(tǒng)功能時展開的四個窗口。建模過程建模思路與根本流程采用AQUASIM2.0建模，建模的根本思路和流程圖4-9所示。首先是按照模型設(shè)計時選定的設(shè)計參數(shù)、變量及模擬過程在AQUASIM設(shè)計模型，在通過設(shè)置計算的步寬、步長和計算的對象，進展設(shè)計計算。計算后選定參數(shù)進展敏感度分析，根據(jù)敏感度分析結(jié)果，選取幾個敏感度較大的參數(shù)進展模型率定〔即參數(shù)估值〕。調(diào)出率定后模型的圖線，檢查實驗模擬的過程規(guī)律是否與實際情況相符。對于模擬效果不符合實際規(guī)律的情況，需重新檢查對模型進展調(diào)適，直至得到與現(xiàn)實契合程度較高的模型圖線為止。圖4-9AQUASIM建模的根本流程設(shè)計模型設(shè)計變量本設(shè)計設(shè)計五種類型的變量：狀態(tài)變量、常數(shù)變量、真實數(shù)值變量、方程變量、程序變量。各種類型變量包含的模型設(shè)計參數(shù)情況如表4-2所示。常數(shù)變量率定的最小值一般為設(shè)計值的1/10，最大值一般為設(shè)計值的10倍，但異養(yǎng)菌衰減系數(shù)bH、異養(yǎng)菌產(chǎn)率系數(shù)YH除外，bH的最大率定值取1，YH的最大率定值取0.8。表4-2模型各參數(shù)類型包含的設(shè)計參數(shù)參數(shù)類型狀態(tài)變量常數(shù)變量真值變量方程變量程序變量包含的設(shè)計參數(shù)SS、*BH、*P、*SYH、bH、fp、Kh、KOH、KS、K*、μH、SS_ini、SS_input、*BH_iniSS_real、*BH_realQ、Sot說明：設(shè)計參數(shù)中，參數(shù)YH、bH、fp、Kh、KOH、KS、K*、μH、t、Q、So、SS、*BH、*P、*S具有一般意義。變量SS_ini、SS_input、*BH_ini、SS_real、*BH_real分別代指初始SS濃度、輸入SS濃度、初始*BH濃度、實驗測得出水COD數(shù)據(jù)、實驗測得MLSS*0.75數(shù)據(jù)。由于實驗過程采用溶解性易降解的葡萄糖為唯一碳源，實驗測得的COD值可近似等價于底物濃度SS含量。另外，由于對于同一系統(tǒng)，f=MLSS/MLVSS為一定值，因此將實驗測得的MLSS值乘上轉(zhuǎn)換系數(shù)f則可得知體系中微生物的含量*BH。由于在實驗過程中測得的進水COD值變化較大，與理論值差異較大，為保證輸入碳源能滿足微生物代需要，本設(shè)計進水和初始的SS值都取實驗測得的最大值547mgCOD/L，初始*BH值取423*0.75=313mgCOD/L?？紤]到在實驗過程中由于人員誤差、儀器誤差帶來的實驗結(jié)果誤差較大，結(jié)合理論知識和相關(guān)經(jīng)歷，本設(shè)計5d的出水COD值同一取0，而MLSS數(shù)據(jù)DAY2-4的實驗數(shù)據(jù)作為模型的1-3天的MLSS數(shù)據(jù)，轉(zhuǎn)化為*BH_real時MLSS還需乘上轉(zhuǎn)換系數(shù)f。表4-3實驗測得MLSS值與模型*BH_real的關(guān)系轉(zhuǎn)換實驗數(shù)據(jù)模型數(shù)據(jù)時間MLSS/(mgCOD/L)時間MLSS/(mgCOD/L)*BH_real/(mgCOD/L)DAY1420initial420313DAY21610DAY116101208DAY32260DAY222601965DAY42580DAY325801995設(shè)計過程設(shè)計過程包括異養(yǎng)菌的好氧生長、異養(yǎng)菌的衰減和網(wǎng)捕性有機物的水解三種反響過程，具體方程式和化學計量參數(shù)形式按3.3.和3.4.要求編輯即可。設(shè)計反響器設(shè)計的反響器包括生化反響器(reactor)和沉淀池(clarifier)兩局部，反響器具有固定的容積。按照模型的假設(shè)前提，微生物反響過程僅在生化反響器進展，沉淀池不涉及任意反響過程。兩個反響器的具體設(shè)計情況如表4-4所示。表4-4設(shè)計反響器的具體設(shè)計情況反響器狀態(tài)變量涉及的反響初始條件輸入條件生化反響器SS、*BH、*P、*S異養(yǎng)菌的好氧生長異養(yǎng)菌的衰減網(wǎng)捕性有機物的水解SS=SS_ini*BH=*BH_iniSS=Q*SS_input沉淀池SS、*BH、*P、*S///設(shè)計聯(lián)結(jié)對于本設(shè)計而言，聯(lián)結(jié)的對象為生化反響器和沉淀池。水從生化反響器流向沉淀池，污泥則從沉淀池回流至生化反響器。本設(shè)計回流泥水混合液流量設(shè)為0.5*Q。由于模型假設(shè)在沉淀池中泥水充分別離，因此固體*P、*S、*BH全部回流，而溶解性的SS僅有50%參與回流。設(shè)計計算本設(shè)計計算采用的計算步寬為0.1d，步數(shù)為100，則對應的計算時長為10d。計算對象為反響器底物濃度SS和微生物生物量*BH分別隨時間的變化。敏感度分析本設(shè)計中參與敏感度分析的變量包括：YH、bH、fp、Kh、KOH、KS、K*、μH。各參數(shù)對于底物濃度變化和微生物量變化過程的敏感度和誤差奉獻情況如表4-5所示。由表中數(shù)據(jù)可知，參數(shù)μH、KS、YH、bH、KOH對底物濃度變化的影響較大，而所有參數(shù)對生物量變化過程都有較大影響。表4-5參數(shù)對底物濃度變化與微生物量變化過程的敏感度和誤差奉獻情況底物濃度變化生物量變化參數(shù)敏感度(mgCOD/L)誤差奉獻(mgCOD/L)參數(shù)敏感度(mgCOD/L)誤差奉獻(mgCOD/L)μH31.835.305YH326.5487.3KS14.190.7095bH28.2570.62YH6.5239.736μH19.173.194bH4.3310.82KS8.4960.4248KOH1.6664.165Kh4.5622.064Kh0.087640.03966K*3.61824.12K*0.061240.4083fp2.3329.13fp0.019150.2393KOH1.2463.116模型率定本設(shè)計是將模型模擬的底物濃度隨時間的變化、微生物生物量隨時間的變化兩個過程，與實驗過程測得的底物濃度和微生物生物量變化進展率定。即將狀態(tài)變量SS、*BH分別與真值變量SS_real、*BH_real進展率定。根據(jù)敏感度分析的結(jié)果，YH、bH、fp、Kh、KOH、KS、K*、μH八個參數(shù)對過程的影響都較大，考慮到生物量中可轉(zhuǎn)化為顆粒性產(chǎn)物的比例fp對于活性污泥而言變化不大，因此不考慮f