水利學報2002年7月SHUILIXUEBAO文章編號: 0559-9350 (2002) 07-0053-04排放口污水處理率的優(yōu)化計算薛聯(lián)青'，陸桂華'《1.河海大學水環(huán)院江蘇甫京210098)摘要:依據(jù)文中建立的河道水質(zhì)模型，以水環(huán)境目標為約束，污水處理費用作為排放口污水處理控制方案優(yōu)選的評價依據(jù)，提出了排放口污水處理綜合規(guī)劃方程的拉格朗日乘子求解方法，并應用于淮河流城某區(qū)城的污水治理規(guī)劃。計算結(jié)果表明該方法較傳統(tǒng)規(guī)劃方程的數(shù)值求解及懲罰函數(shù)法操作簡便，對央策者進行排放口污水處理控制方案的初選具有指導意義。關鍵詞:污染控制;排放口;解析方法;優(yōu)化;數(shù)值計算.中圈分類號: X522文獻標識碼: A由于工業(yè)與城市經(jīng)濟的高度發(fā)展，工業(yè)廢水和生活污水日益增加，大量的廢棄物流人江河朔泊，造成了嚴重的水體污染，河流的環(huán)境承載力下降。與此同時各污水處理廠的污水處理費用也逐年增高。如何合理的配置污水處理廠，減少整個區(qū)域處理運行費用，就要求我們用系統(tǒng)分析的方法，在河.流水質(zhì)狀況模擬的基礎上，進行排放口優(yōu)化處理計算。排放口優(yōu)化處理問題是水污染控制規(guī)劃計算中常見的一種，但在實際應用中，對污水處理率的優(yōu)化計算通常可形成凸規(guī)劃問題，尤其當有多約束條件時，若采用簡單的解析算法求解，即可避免繁雜的數(shù)值計算。同時，在-定程度上也避免了數(shù)值迭代計算誤差所引起的處理效率的不合理現(xiàn)象。1處理率優(yōu)化計算 的解析方法設河道的設計流量為Qo，上游來水中的污染物濃度為Co， 河道斷面水質(zhì)標準為C',第i個排放口排放的污水流量為qi，污水初始濃度為l，處理率為n，則相應的污水處理費用為Z=aqf+ar q7。當河流中水質(zhì)變化符合Streter-Phelps 模型!"1所描述的規(guī)律時，可建立排放口污水處理率的最優(yōu)化模型表示如下:minZ= 2z= b。+ bini+b2≌+...+ b,公.s.t.a1179 + a27 + "aIn7n≤c -d= cia217: + az2 +"*a2n7n≤c'-dx=ci(1)a.s7i + 0272 + .'.m7n≤c'-d, = c;式中:b,(i=0,1,.，n)與q，an,ar及β值有關;ag(i=1,2,.，n;j=1,2,.，n)與Qo，q., 4;及耗氧系數(shù)k有關; d, 是與ag相關的常數(shù)項; C* 為監(jiān)測斷面的水質(zhì)約束; Z為所有排放口污水處理的總費用。式(1)中的目標函數(shù)是二次，約柬條件呈線性，因此屬于非線性規(guī)劃中的凸規(guī)劃類型["。收稿日期: 20-191作者藺介:薛聯(lián)青(1973-),女，新疆石河子人，博士研究生，主要從事水資源規(guī)劃和環(huán)境保護的研究。-53-中國煤化工TTTTJMYHCNMHG.1.1懲罰函數(shù)法構(gòu)造懲罰函數(shù)[|如下:φ= b。+ 2bp+ r{Z[max(0,aun +ap+-+an.-ci)}(2)可將式(1) 的有約束的規(guī)劃問題變成無約束的規(guī)劃問題。在式(2) 中，y稱為懲罰因子，為一取值任意大的正數(shù)。由式(2) 知，只要有一個約束條件不滿足，φ將趨于無窮大，而在約束條件都滿足的情況下，φ值與z值相同，故使ψ有最小值的處理率組合便是使z有最小值，且又滿足約束條件的處理率組合。對于式(2)中的辦分別求導，并考慮導數(shù)為零及不滿足約東的情況得:里= 2b,7 +2rZan(anη +oa72 +. +a7.-ci) =0呢= 2b.272 +27 Za(auq +a2+.+awη.-ci)=0(3)器=26. +22_a(an+-+on-ci)=0整理后并注意到γ→∞，得: Sam + 2aua72+-+ 2auan7 = Zcian; Sana7i +二咖7+--+ Zaawn. = ..a..... Eana7+ 2anaB+-+ Ean = Ecia.Z Sana" Zanan式(1)是一個線性方程組，其解為:q。= -;A =Zana Z略- 2oa....... =l Sanan Sa-.J「EaEana. Zcian". Sana.Eanaa 2哈-. Zcian". Zan .則由上述構(gòu)造的懲罰函數(shù)即可求解排放口污水處理的效率l Sanan Zan."" 2cian...2a2J和污水處理費用。1.2 拉格朗日乘子法 由于在水污染控制中，通常選擇具有代表性，對河流水質(zhì)起主要影響作用的污染因子作為約束。鑒于淮河污染主要表現(xiàn)為有機污染，因此本文選取化學耗氧量COD≤15mg/1 (以.COD.依度表征)和溶解氧DC≥5ng/l作為河流水質(zhì)的衡量指標和約束條件!1，以整個河段上所有排放口污水處理的費用之和為目標函數(shù),采用拉格朗日乘子法，將有約束的極值問題轉(zhuǎn)化為無約束極值進行求解。該方法較前文中的懲罰函數(shù)法及數(shù)值求解方法簡便,并易于操作。根據(jù)河流COD和DO的響應矩陣U和V,及河流各斷面的COD、DO約柬向量L'和0°,并設排放口最大允許輸人河流各斷面的COD濃度向量為L,河段中的COD濃度的變化用S-P模型所描述。UL+m≤L°則在滿足約柬T的條件下， 根據(jù)方程(1), 建立以》為變量的污水處理的費用函數(shù)z，并且求出η. 使Z為最小，即滿足minZ=bo+br亦+bz≠+...+ b.好.以排放口的污水處理率為變量,建立各個斷面的約束條件。第1控制斷面: c,= Colo+Lh9i一知≤c'; C"。= Coe%。第2控制斷面: C,=Q.+gC, (00+9) +6292 (1- 2) .: C*; C'= Ce-1. ....第i控制斷面:C=Qo+g:+92C".. (+..+..) +6.9. (1-N)≤c;C'= C..eQo+q1+q2+..+ qs將約束條件轉(zhuǎn)換為表達式: h(》)=0 (i=1, 2, .，n)的形式。處理效率7是可以用排放口54 一中國煤化工MHCNMHG“-----T蔓u污水的初始濃度l所表示。建立費用函數(shù)的拉格朗日函數(shù)L (η, 2) = C+2,hj (η) + λrhz (η)+..+λ.h. (η)。其穩(wěn)定點(η", λ")應滿足下列方程:VqL(n",i")= VJ(n")+A; Vh(η°)+. + λiVh,(q")+.. + a:Vh.(7") = 0.若穩(wěn)定點存在，則穩(wěn)定點即為所求費用函數(shù)的最優(yōu)解。2排放口污水處理的優(yōu)化計算根據(jù)淮河流域某區(qū)域的地形、地貌特點以及原有及擬建排放口的情況，沿河段設置5個集中的排放口。擬規(guī)劃的5個排放口在預測年的排污量q分別為0.255、0.939、1.098、 7.529和7. 665m'/a,排放口污水初始COD依度L (i=I, 2, .5)分別為347.8、36.17、 28.73、 88.34 和140.1mg/1。溶解氧Do濃度分別為6.17、7.60、 9.12、4.52 和4.94ng/1。為了計算簡便，本文僅采用COD與DO作為主要約束因子，選擇來水保證率為P=75%，流量Q=55m'1s,河流流速v=5.5km/d。根據(jù)河流各排放口的位置進行概化，計算水流流經(jīng)各河段的時間分別為: rp.=1.0 (d)，T.2=2.5 (d); T3=1.5 (d)， T4=2.7 (d)， T4_s=0.5 (d)。同時其非線性的污水處理費用目標函數(shù)為: Z;=200q." + 1000qi*n.式中: 2.為污水處理的費用; q為污水處理的規(guī)模; y為污水處理的效率。由此可確定各排放口的費用函數(shù)分別為z =67.0+335.1積，z= 190.2+ 950.9站，z, =215.5+ 1077.7樂，44= 1005.6+ 5028.0m，Zs = 1020.0 +5100.0ξ .由于費用函數(shù)是非線性函數(shù)，因此在進行數(shù)值求解時，為了使線性函數(shù)能較好地代表原函數(shù)，通常采用分段線性化方法，即將污水處理效率分成幾個區(qū)間，每一區(qū)間用一線性函數(shù)代表原函數(shù)。不失一般性，假定在0≤η<1的區(qū)間里，對費用函數(shù)分三段實行線性化處理，分別計算各段的線性斜率"，效率分級為: 0≤劑≤0.3, 0.3≤取≤0.85, 0.85≤q≤1 (i=1, .5)，對所建立的規(guī)劃模型分別用數(shù)值及解析方法進行污水處理效率和處理費用的優(yōu)化計算。2.1用敷值方法求解模型若淮河淮南段河流水質(zhì)監(jiān)測斷面均采用國家三級標準作為約束控制，則河段上如果設置s個污水排放口，則計算機編程進行數(shù)值求解可得s個排放口污水處理率及其各自的污水處理費用，計算結(jié)果如表1所示。累計各個排放口的污水處理費用，可計算得當設有5個污水處理廠時，所設定的各排放口污水處理的總費用約為5919.75萬元。表」5個捶放口的污水處理費用計算標號液量/ (Vr)膿度1 (my/L)處理率(%}費用/萬元斷面COD/ (mg/L)斷面DO/ (nm/L)0.2550347.800.000.09.455.860.939036.07 .0.5.786.431.098028.734.567.087.529088.34 .49.493261.335.647.687.6650140.10 .29.982658.4215.001.35，2.2規(guī)劃模型的解析求解分析設各排放口初始污水的 COD濃度為1,并且認為河段中的COD濃度的變化可用-級動力學反應方程描述，則可以建立排放口處理費用的目標函數(shù): Z= 2z = 200q%*+ 00qm..目標函數(shù)及約柬條件中的排放口污水處理效率》可以排放口污水初始濃度l表示，即處理后輸人河流各斷面COD濃度為L= (1-η) I,采用拉格朗日乘子法進行求解，可求得5個集中排放口的污水處理效率為0.065, 0.019, 0.023, 0.458, 0.248， 總的污水處理費用約為5552萬元，與懲罰函數(shù)法計算結(jié)果相一致，并與數(shù)值求解計算結(jié)果相接近。鑒于第四與第五排放口相距較近，合并后新廠的污水處理流量為15. 194m'/s,排放污水的初始CoD濃度為114. 45mg/L,溶解氧濃度DO為4.731m/L.合并計算得污水處理費用函數(shù)為z. = 1763.455 -中國煤化工YHCNMHG+8817弱，同法可求得各排放口污水處理的效率及費用，并將數(shù)值計算結(jié)果列于表2。由結(jié)果可看出，當僅設4個污水處理廠時的污水處理的大致費用為3733.08萬元。較5個排放口的污水處理費用大大降低。表2 4個排放口的污水處理費用計算標號流量/ (ts)依度/ (mg/L)處理事(%)費用1萬元斷面COD/ (mg/L)斷面DO/ (ng/L)0.2550347.800.009.425.840.939036.070.5.726.411.098028.734.507.06I5. 1940114.45 .44.493733.0815.007.472.3結(jié)果分析若對模型 的求解采用上節(jié)所介紹的懲罰函數(shù)法以及拉格朗日乘子法進行求解，各排放口的條件和約束條件均不變?？捎嬎愕贸? (1)對于5個排放口時，污水處理效率分別為: 0.065,0.019，0.023, 0.458, 0.248，總的污水處理費用大約為5552萬元。(2) 對于4個排放口時，計算得處理效率分別為: 0.066, 0.019, 0.025, 0.42, 污水處理的總費用大約為3749萬元。4結(jié)論通過以上實例計算表明，拉格朗日乘子法與數(shù)值求解所求得的結(jié)果基本相近,并與懲罰函數(shù)法計算結(jié)果相一致。由于在數(shù)值求解中對費用函數(shù)進行了線性化近似，并且在計算選代中產(chǎn)生一定的誤差，由此導致了與解析求解的結(jié)果的偏差。計算表明，文中拉格朗日乘子法具有計算簡便的特點，尤其在做初步規(guī)劃時，可以較為方便地做出估算，對選擇方案具有指導意義。同時在對多約束條件的選擇上具有可選擇性，并且在有多約束時，拉格朗日的乘子法求解較傳統(tǒng)的數(shù)值求解及懲罰函數(shù)法易于操作，解算簡便，特別有利于規(guī)劃人員的設計計算。參考文獻:[1]付國偉，程聲通.水污染控制系統(tǒng)規(guī)劃[M].北京:猜華大學出版社，1985.[2]黃平，排放口污水處理宰優(yōu)化計算的解析方法[J].環(huán)境科學與技術(shù)，1997.[3]李維錚，等.運籌學[M].北京:清華大學出版杜，1982.[4]傅國偉.河灌水質(zhì)數(shù)學模型及其模擬計算[M]. 北京:中國環(huán)境科學出版社，1987.Optimization of treatment rate for discharged urban wastewaterXUE Lian-qiog' ，LU Gui-hua'(1. Hohai Unieriy, Narnjing 210098， China)Abstract: On the basis of river water model established, the Lagrange muliplier method for solving thecomprehensive planning equation of treatment rate for dischanged urban water is proposed. In this methodthe water environment is 8el as the objective constraint and the co8t of treatmnent is referred 88 the 88-ment slandard for optimization of pollution control. The method is applied to the wastewaler dischargeplanning of a river section in Huaihe River. It's indicated thal the method is suitable for the contralscheme selection to decision-makers .Key words: pllution control; wastewaler dischange outlet; treatment rate; optimization一56一中國煤化工MHCNMHG