城市工業(yè)COD總量優(yōu)化分配研究

來源：http://www.xfzhuw.com/ 作者：余氯檢測儀 時(shí)間：2019-10-11

　　摘要 總量控制是控制污染源發(fā)展趨勢、改善環(huán)境質(zhì)量、實(shí)現(xiàn)經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展的重要途徑,如何在適度公平的基礎(chǔ)上尋求環(huán)境、經(jīng)濟(jì)、技術(shù)、資源等整體最優(yōu)是環(huán)境科學(xué)領(lǐng)域的研究課題。本文以工業(yè)城市蘇州市為例,研究其工業(yè)化學(xué)需氧量的排放特征,以基尼系數(shù)法分行業(yè)對(duì)其工業(yè)化學(xué)需氧量排放量的公平性進(jìn)行評(píng)價(jià),并將總量控制與資源、社會(huì)和經(jīng)濟(jì)相聯(lián)系,以行業(yè)經(jīng)濟(jì)效益最大化和增加治理投資費(fèi)用最小化為目標(biāo),利用多目標(biāo)行業(yè)總量優(yōu)化分配模型對(duì)蘇州市的工業(yè)化學(xué)需氧量排放總量進(jìn)行優(yōu)化分配。研究結(jié)果表明,紡織業(yè)、化學(xué)原料及化學(xué)制品制造業(yè)、能源和水的生產(chǎn)與供應(yīng)業(yè)、造紙及木材加工、醫(yī)藥制造業(yè)等行業(yè)是蘇州市的化學(xué)需氧量重點(diǎn)排放行業(yè),經(jīng)優(yōu)化分配后,COD排放總量削減了10%,新鮮用水量減少了41.81%,行業(yè)年總產(chǎn)值增幅達(dá)到214.69%,資源和水環(huán)境容量在滿足一定的經(jīng)濟(jì)增長速度的條件下實(shí)現(xiàn)優(yōu)化配置,總量控制制度在市場經(jīng)濟(jì)體制下發(fā)揮出盡可能大的環(huán)境效益和經(jīng)濟(jì)效益。
中國論文網(wǎng) /7/view-3083670.htm
　　關(guān)鍵詞 化學(xué)需氧量(COD);公平性;基尼系數(shù);優(yōu)化分配
　　中圖分類號(hào) X321 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1002-2104(2010)03-0124-06 doi:10.3969/j.issn.1002-2104.2010.03.023
　　
　　總量控制是控制污染源發(fā)展趨勢、改善環(huán)境質(zhì)量、實(shí)現(xiàn)經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展的重要途徑,各個(gè)排污單位或污染源之間如何科學(xué)、合理地分配排放污染物的量,使其在市場經(jīng)濟(jì)體制下發(fā)揮出盡可能大的環(huán)境效益和經(jīng)濟(jì)效益,是總量控制的核心問題[1]。水污染物排放總量分配涉及到經(jīng)濟(jì)、社會(huì)、技術(shù)和環(huán)境等多種因素, 在每種因素中往往又包含了若干種定性和定量因子[2]。國外優(yōu)化分配的研究自20世紀(jì)60年代起步,80年代后各種分配模型的應(yīng)用得到很大改善,如開發(fā)了隨機(jī)動(dòng)態(tài)規(guī)劃費(fèi)用最小模型[3]、利用模糊多目標(biāo)優(yōu)化方法等[4-5],國內(nèi)也有一些總量控制公平性方面的研究,如最小處理費(fèi)用公平協(xié)調(diào)方法[6]、水污染物允許排放量公平分配方法[7-9]、多目標(biāo)行業(yè)總量優(yōu)化分配模型[10-11],但就我國目前國情而言,如何在適度公平的基礎(chǔ)上,尋求環(huán)境、經(jīng)濟(jì)、技術(shù)、資源等整體最優(yōu)仍是一個(gè)重要的研究課題。
　　本文以蘇州市為例,將總量控制與資源、社會(huì)和經(jīng)濟(jì)相聯(lián)系,研究該市工業(yè)化學(xué)需氧量(COD)的排放特征,以基尼系數(shù)法分行業(yè)對(duì)其工業(yè)COD排放量的公平性進(jìn)行評(píng)價(jià),同時(shí)考慮水資源量、水環(huán)境容量、社會(huì)行業(yè)優(yōu)化調(diào)整等多方面的約束,以行業(yè)經(jīng)濟(jì)效益最大化和增加治理投資費(fèi)用最小化為目標(biāo),參考利用多目標(biāo)行業(yè)總量優(yōu)化分配模型[10-11]對(duì)蘇州市的工業(yè)COD排放總量進(jìn)行優(yōu)化分配,力求實(shí)現(xiàn)盡可能大的環(huán)境效益和經(jīng)濟(jì)效益。
　　1 研究對(duì)象概況
　　蘇州市位于長江三角洲中部,江蘇省南部,東鄰上海,南接浙江,西傍無錫,北依長江。蘇州水系受滬寧鐵路的分割,分為南北兩大水系。北部地區(qū)由于地勢相對(duì)較高,水系主要以河網(wǎng)為主。而南部地區(qū)由于地勢低洼,水系以湖蕩為主。蘇州屬典型的江南水網(wǎng)化平原地區(qū),全市水面約?3 609 km²,?約占總面積的42.5%。蘇州擁有大小湖泊321個(gè),湖泊率約為8 %,是全國湖泊率的10倍。
　　2006年蘇州市三次產(chǎn)業(yè)比例為1.9∶65.4∶32.7,電子、紡織、交通設(shè)備、冶金、電氣、化工是支柱產(chǎn)業(yè),水污染點(diǎn)源以鋼鐵行業(yè)、紡織行業(yè)、化工行業(yè)、造紙業(yè)、印染行業(yè)等為主。工業(yè)企業(yè)廢水COD排放總量按行業(yè)分類統(tǒng)計(jì),全市排在前三位的行業(yè)為紡織業(yè)、化學(xué)原料及化學(xué)制品制造業(yè)和造紙及紙制品業(yè)。
　　2 研究基本原理
　　2.1 基尼系數(shù)
　　基尼系數(shù)的定義是在洛倫茲曲線的基礎(chǔ)上提出的,洛倫茲曲線如圖1。
　　用橫軸代表人口累積百分比,用縱橫代表收入累積百分比,實(shí)際收入分配曲線和收入分配絕對(duì)平等曲線之間的面積為A,實(shí)際收入分配曲線右下方面積為B,則基尼系數(shù)的定義式為:
　　G=SAS?(A+B)?(1)
　　
　　圖1 洛倫茲曲線
　　Fig.1 Lorenz Curve
　　楊占紅等:城市工業(yè)COD總量優(yōu)化分配研究中國人口•資源與環(huán)境 2010年 第3期基尼系數(shù)有多種計(jì)算方法,如直接計(jì)算法、回歸曲線法、人口等分法、城鄉(xiāng)分解法[12],葉禮奇提出了三角形面積法、弓形面積法[13],本文采用比較普遍使用的梯形面積法,將洛倫茨曲線下方的面積近似為若干梯形進(jìn)行計(jì)算,其公式如下:
　　G?ini?=1-?∑ni=1(Xi-X?i-1?)(Yi+Y?i+1?)?(2)
　　式中,Xi為評(píng)估指標(biāo)的累計(jì)比例,%;Yi為污染物的累計(jì)比例,%;i為分配對(duì)象數(shù)量。當(dāng)i=1時(shí),(X?i-1?,Y?i-1?)視為(0,0)。
　　2.2 多目標(biāo)行業(yè)總量優(yōu)化分配模型[10-11]
　　2.2.1 行業(yè)水環(huán)境經(jīng)濟(jì)綜合指數(shù)(R)
　　某行業(yè)的R可以用式描述:
　　R=F(ec,en,r)(3)
　　式中:ec表示該行業(yè)經(jīng)濟(jì)方面的特征,en表示該行業(yè)環(huán)境方面的特征,r表示該行業(yè)資源方面的特征,F表示一組復(fù)雜的函數(shù)關(guān)系。
　　假定共選擇m項(xiàng)指標(biāo)作為R的影響因子,先將每項(xiàng)因子的實(shí)際值與區(qū)域社會(huì)平均值相除,得到單因子社會(huì)發(fā)展度,再運(yùn)用數(shù)學(xué)方法調(diào)整m項(xiàng)因子對(duì)R的影響權(quán)重,綜合計(jì)算得到R值。在行業(yè)總量優(yōu)化分配模型中將R作為調(diào)整行業(yè)結(jié)構(gòu)這一關(guān)鍵約束條件來加以使用。
　　2.2.2 模型結(jié)構(gòu)
　　(1)目標(biāo)函數(shù)
　　1)行業(yè)年總產(chǎn)值總和S(10?4元•a-1?)最大:
　　MaxS=?∑ni=1X?iK?1i??(4)
　　式中:Xi―行業(yè)i控制污染物優(yōu)化排放量,t;
　　K?1i?―行業(yè)i萬元產(chǎn)值排放控制污染物量系數(shù),t/萬?元;?
　　i―總量控制區(qū)域內(nèi)行業(yè)類別。
　　2)行業(yè)新增處理污染物費(fèi)用總和G最小:本模型中假定生產(chǎn)規(guī)模保持不變,當(dāng)行業(yè)控制污染物現(xiàn)狀排放量Ci高于Xi時(shí),G為正;當(dāng)Ci低于Xi時(shí),G為負(fù)。本模型只考慮限制污染物排放量變化所產(chǎn)生的影響。
　　MinG=?∑Ni=1(Ci-Xi)PiK?2i??(5)
　　式中:Pi―行業(yè)i廢水中控制污染物平均濃度指標(biāo),t/萬t;
　　?K?2i?―行?業(yè)i處理單位體積廢水費(fèi)用系數(shù),萬元/?萬t。?
　　(2)約束條件
　　1)水環(huán)境質(zhì)量約束:為保證區(qū)域水環(huán)境質(zhì)量,各行業(yè)每年所排放控制污染物總量應(yīng)不大于本區(qū)域控制污染物的環(huán)境容量QR(t)。
　　?∑Ni=1Xi≤QR(6)??2)水資源合理利用約束:為保證當(dāng)?shù)厮Y源可持續(xù)利用,各行業(yè)每年所用新鮮水量總和應(yīng)不大于當(dāng)年本區(qū)域可供給行業(yè)生產(chǎn)的新鮮水資源量QG(萬t)。
　　?∑Ni=1XiK?1i?K?3i?≤QG?(7)
　　式中:K?3i?―行業(yè)i萬元產(chǎn)值新鮮用水量系數(shù),萬噸/萬元。
　　3)行業(yè)優(yōu)化調(diào)整約束:為保證區(qū)域行業(yè)布局結(jié)構(gòu)科學(xué)合理,水環(huán)境容量、水資源得到優(yōu)化配置,行業(yè)間的容量配額比應(yīng)與水環(huán)境經(jīng)濟(jì)綜合指數(shù)Ri其比例相近,相適應(yīng)程度通過上下限兩個(gè)調(diào)整系數(shù)控制。??K?min?•Ri∑Ni=1Ri≤Xi∑Ni=1Xi≤K?max?•Ri∑Ni=1Ri?(8)
　　式中:K?min?―行業(yè)優(yōu)化調(diào)整下限系數(shù);
　　K?max?―行業(yè)優(yōu)化調(diào)整上限系數(shù)。
　　4)行業(yè)分配容量非負(fù)約束:即:Xi> 0
　　(3)模型求解
　　為求解方便,本規(guī)劃將目標(biāo)函數(shù)由多目標(biāo)轉(zhuǎn)化為單目標(biāo),將非線性約束轉(zhuǎn)化為線性約束,轉(zhuǎn)化后模型形式如下:

　　?Max(S-G)=∑Ni=1XiK?1i?-∑Ni=1(Ci-Xi)PiK?2i?(9)
　　約束條件為:
　　∑Ni=1Xi≤QR
　　∑Ni=1XiK?1i?K?3i?≤QG
　　Xi≥K?min?•Ri∑Ni=1Ri•∑Ni=1Xi,i=1,2,L N
　　Xi≥K?miax?•Ri∑Ni=1Ri•∑Ni=1Xi,i=1,2,L N?
　　3 研究結(jié)果
　　3.1 工業(yè)COD排放特征
　　3.1.1 排放量與排放強(qiáng)度
　　根據(jù)《中國環(huán)境統(tǒng)計(jì)年報(bào)》、《蘇州市環(huán)境狀況公報(bào)》和《蘇州統(tǒng)計(jì)年鑒》對(duì)蘇州市1993-2006年間工業(yè)COD排放量和蘇州工業(yè)產(chǎn)值數(shù)據(jù)統(tǒng)計(jì),得出蘇州市近年來工業(yè)COD年排放量和排放強(qiáng)度的變化趨勢,分別見圖2和圖3。
　　
　　圖2 1993-2006年蘇州市工業(yè)COD排放量變化趨勢
　　Fig.2 The emission trend of industry COD in? Suzhou between 1993 and 2006
　　圖3 1993-2006年蘇州市工業(yè)COD排放強(qiáng)度變化趨勢
　　Fig.3 The emission intensity trend of industry COD in? Suzhou between 1993 and 2006從圖2中可以看出,1993-2006年14年間蘇州市工業(yè)COD年排放量平穩(wěn)中稍微增長,期間存在下降與上升的波動(dòng),在2005年排放量達(dá)到峰值6.71萬t。從圖3中可以看出,1993-2006年間蘇州市工業(yè)COD排放強(qiáng)度有了大幅下降,降低幅度達(dá)到85%,亦存在一定的升降波動(dòng)。
　　3.1.2 分行業(yè)排放量與放量強(qiáng)度
　　本研究中共調(diào)查了873家重點(diǎn)排污企業(yè),這些污染點(diǎn)源分屬于30余個(gè)子行業(yè)。為了便于宏觀層面上的管理,根據(jù)行業(yè)性質(zhì)經(jīng)合并得到17個(gè)行業(yè),分別為:能源和水的生產(chǎn)與供應(yīng)業(yè)、紡織業(yè)、造紙及木材加工、化學(xué)原料及化學(xué)制品制造業(yè)、食品加工制造業(yè)、醫(yī)藥制造業(yè)、冶金及金屬加工制造業(yè)、非金屬礦物制品業(yè)、電子及通訊設(shè)備、儀器儀表及文化辦公用機(jī)械制造業(yè)、交通運(yùn)輸設(shè)備制造業(yè)、工藝品及其他制造業(yè)、石油、煉焦及核燃料加工業(yè)、采選業(yè)、廢棄資源和廢舊材料回收加工業(yè)、研究與試驗(yàn)發(fā)展、公共設(shè)施管理業(yè)。各行業(yè)COD年排放量和萬元產(chǎn)值排放量情況見表1。
　　從表1可以看出,蘇州市各行業(yè)COD年排放量和排放強(qiáng)度差距較大,COD年排放量最高可相差20 653倍,排放強(qiáng)度最高相差3227倍,但單從行業(yè)年排放量和排放強(qiáng)度分析公平性是不盡合理的,下面采用基尼系數(shù)法對(duì)蘇州市各行業(yè)COD排放量的公平性進(jìn)行分析。
　　3.2 公平性評(píng)價(jià)
　　以各行業(yè)工業(yè)產(chǎn)值為評(píng)估指標(biāo),工業(yè)產(chǎn)值累積比例作為橫坐標(biāo),COD排放量的累積比例為縱坐標(biāo),以各行業(yè)工業(yè)生產(chǎn)總值-COD排放量繪制洛倫茨曲線,見圖4。
　　根據(jù)曲線圖應(yīng)用梯形面積法計(jì)算基尼系數(shù),通過計(jì)算得到各行業(yè)生產(chǎn)總值-COD排放量的基尼系數(shù)為0.652。
　　在經(jīng)濟(jì)學(xué)中,由于社會(huì)發(fā)展的局限性,人均收入的分配不可能平衡,因此基尼系數(shù)在0-0.2的可能性很小,基 尼系數(shù)的合理范圍為0.2-0.4。但在流域間環(huán)境基尼系
　　
　　圖4 蘇州市各行業(yè)生產(chǎn)總值-COD排放量洛倫茲曲線圖
　　Fig.4 Lorenz Curve of different industry total output ?value-COD emission in Suzhou city數(shù)的計(jì)算中,由于各流域間沒有水資源沖突,沒有不平等的前提,認(rèn)為基尼系數(shù)趨近于零是最合理的,因此本文將基尼系數(shù)的合理范圍界定為0-0.2。
　　根據(jù)以上計(jì)算結(jié)果可知,蘇州市工業(yè)生產(chǎn)總值指標(biāo)的基尼系數(shù)超過了合理范圍,主 要是由于紡織業(yè)、化學(xué)原料及化學(xué)制品制造業(yè)、能源和水的生產(chǎn)與供應(yīng)業(yè)、造紙及木材加工、醫(yī)藥制造業(yè)等行業(yè)污染物排放量過高,加之電子及通訊設(shè)備等低COD排放行業(yè)起步較晚,發(fā)展緩慢所致。由于行業(yè)基尼系數(shù)較高,需對(duì)指標(biāo)進(jìn)行調(diào)整,處于曲線頂點(diǎn)的行業(yè)是重點(diǎn)污染對(duì)象,是污染物削減排放的主要對(duì)象。
　　3.3 總量優(yōu)化分配
　　3.3.1 模型所需系數(shù)及邊界條件確定
　　(1)模型所需參數(shù)
　　將17個(gè)行業(yè)的污染點(diǎn)源數(shù)據(jù)進(jìn)行匯總,得到年總產(chǎn)值、廢水治理設(shè)施運(yùn)行費(fèi)用、年新鮮用水量、廢水排放量、COD排放量等數(shù)據(jù)的匯總結(jié)果;再經(jīng)過計(jì)算得到模型中所需的各行業(yè)系數(shù)值,包括 K1、K2、K3、C和P,見表1。
　　根據(jù)江蘇省2006年水資源公報(bào),蘇州市當(dāng)年水資源量約為值為19.98億t(QG)。按照“來水Ⅲ類,目標(biāo)Ⅳ類"標(biāo)準(zhǔn)控制測算水環(huán)境容量,得出蘇州市“十一五"期間COD最大允許排放量為129 500 t,扣除2006年城鎮(zhèn)生活COD排放量68 760.39 t,實(shí)際企業(yè)(占排污負(fù)荷約85%)所屬行業(yè)的COD為51 628.67 t(QR)。
　　行業(yè)優(yōu)化調(diào)整系數(shù)上下限:K?min?和K?max?分別取0.9和?1.1。?
　　(2)行業(yè)水環(huán)境經(jīng)濟(jì)綜合指數(shù)
　　采用BP網(wǎng)絡(luò)確定行業(yè)水環(huán)境經(jīng)濟(jì)綜合指數(shù)R,輸入因子為7個(gè),分別為:行業(yè)生產(chǎn)總值、工業(yè)用水總量、新鮮用水量、廢水設(shè)施治理運(yùn)行費(fèi)用、工業(yè)廢水排放量、COD排放量、氨氮排放量;輸出因子為行業(yè)水環(huán)境經(jīng)濟(jì)綜合指數(shù);利用Matlab神經(jīng)網(wǎng)絡(luò)工具箱,經(jīng)過反復(fù)多次試算,確定隱層層數(shù)為1,隱層神經(jīng)元5個(gè),隱層和輸出層傳遞函數(shù)分別為Tansig和Logsig,BP網(wǎng)絡(luò)的訓(xùn)練函數(shù)取Traingdx,誤差為0.001。
　　根據(jù)江蘇省《江蘇省“十一五"工業(yè)結(jié)構(gòu)調(diào)整和發(fā)展規(guī)劃綱要》,江蘇省在今后很長一段時(shí)期里將要大力發(fā)展扶植那些資源節(jié)約、低污染、高效益的重點(diǎn)行業(yè),改造或淘汰那些耗能高、污染重、效益低的行業(yè)。因此,本規(guī)劃在訓(xùn)練BP網(wǎng)絡(luò)的過程中吸納了上述規(guī)劃綱要中的觀點(diǎn),對(duì)于重點(diǎn)行業(yè)給予較高的R值。
　　BP網(wǎng)絡(luò)確定R值及其排名結(jié)果見表1。
　　表1 模型所需參數(shù)確定結(jié)果
　　Tab.1 The result of parameters in model
　　序 號(hào)?Serial ?number行 業(yè)
　　Industries排放強(qiáng)度
　　Emission ?intensity ?(t/10?4 Yuan)?K1處理單位?廢水費(fèi)用
　　Treat unit ?wastewater ?cost? (Yuan/t)
　　K2萬元產(chǎn)值?新鮮水耗
　　Fresh water ?consumption of? ten thousand? Yuan output value ? (t/10?4 Yuan)?K3COD排放量
　　COD ?emission
　　(t/a)
　　C廢水中?COD濃度
　　COD ?concentration? of wastewater ?(t/10?4 t)?P水環(huán)境經(jīng)濟(jì)
　　綜合指數(shù)
　　Water ?resource ?environment -?economic ?comprehensive ?coefficient
　　RR排名
　　R Ranking1能源和水的生產(chǎn)與供應(yīng)業(yè)0.003 2780.836 0 0.044 7546 900.982 60.720.476 8142紡織業(yè)0.007 7330.823 40.009 05018 589.155 80.900.100 0173造紙及木材加工0.004 2180.799 70.004 9016 380.913 90.970.883 312 4化學(xué)原料及化學(xué)制品制造業(yè)0.002 5501.217 00.002 8879 700.776 01.070 .347 5155食品加工制造業(yè)0.001 1822.100 40.001 2721 690.463 61.150.878 213 6醫(yī)藥制造業(yè)0.004 7400.609 30.004 9121 226.538 01.010.950 787冶金及金屬加工制造業(yè)0.000 4721.843 30.001 3904 455.157 20.710.1669168非金屬礦物制品業(yè)0.000 4041.370 20.000 908464.355 80 .720.947 599電子及通訊設(shè)備0.000 2062.933 00.000 3641 390.394 70.690.970 0310儀器儀表及文化辦公用機(jī)械制造業(yè)0.000 5556.945 10.000 5981 633.398 6 1.070.955 6711交通運(yùn)輸設(shè)備制造業(yè)0.000 2245.051 90.000 60145.011 4 0.870.960 3612工藝品及其他制造業(yè)0.000 3785.227 30.000 412548.621 1 1.160.964 4513石油、煉焦及核燃料加工業(yè)0.000 2002.222 20.000 2000 .900 1.001.000 0114采選業(yè)0.000 3595.210 00.000 435593.632 51.130.9 24 01115廢棄資源和廢舊材料回收加工業(yè)0.000 5599.722 20.000 7922.639 01.020.975 4216研究與試驗(yàn)發(fā)展0.505 3070.871 40.543 0171 809.000 0. 930.967 8417公共設(shè)施管理業(yè)0.645 6670.396 80.840 0001 937.000 0. 770.932 510

　　
　　
　　表2 模型優(yōu)化分配結(jié)果
　　Tab.2 Optimal distribution results by the model
　　序號(hào)?Serial ?number行業(yè)
　　Industries優(yōu)化前?Before optimal優(yōu)化后?after optimal工業(yè)總產(chǎn)值
　　Total industrial ?output value
　　(10?4 Yuan)新鮮用水量
　　Fresh water ?consumption
　　(10?4 t)COD排放量
　　COD emission
　　(t)工業(yè)總產(chǎn)值
　　Total industrial ?output value
　　(10?4 Yuan)新鮮用水量
　　Fresh water ?consumption
　　(10?4 t)COD排放量
　　COD emission
　　(t) 1能源和水的生產(chǎn)與供應(yīng)業(yè)2 105 228.3094 217.656 900.98504 343.9122 571.411 653.242紡織業(yè)2 403 986.6021 756.9318 589.1644 838.54405.79346.743造紙及木材加工1 512 801.507 413.736 380.91726 107.863 558.653 062.734化學(xué)原料及化學(xué)制品制造業(yè)3 804 201.7010 984.599 700.78472 513.381 364.151 204.915食品加工制造業(yè)1 430 629.001 820.181 690.462 576 175.453 276.903 045.046醫(yī)藥制造業(yè)258 771.201 270.981 226.54695 447.953 416.043 296.427冶金及金屬加工制造業(yè)9 446 324.6013 127.504 455.161 498 525.012 082.95707.308非金屬礦物制品業(yè)1 149 339.401 043.82464.369 939 110.299 024.714 015.409電子及通訊設(shè)備6 763 317.102 462.761 390.3919 955 112.027 263.664 110.7510儀器儀表及文化辦公用機(jī)械制造業(yè)2 943 821.401 761.171 633.407 035 452.354 207.203 904.6811交通運(yùn)輸設(shè)備制造業(yè)201 290.80120.9045.0118 168 060.4710 919.004 069.6512工藝品及其他制造業(yè)1 452 372.20598.35548.6210 812 224.584 454.644 087.0213石油、煉焦及核燃料加工業(yè)4 500.000.900.9021 189 448.844 237.894 237.8914采選業(yè)1 653 663.00719.25593.6310 907 549.154 744.783 915.8115廢棄資源和廢舊材料回收加工業(yè)4 722.003.742.646 050 209.654 791.773 382.0716研究與試驗(yàn)發(fā)展3 580.001 944.001 809.006 640.943 606.153 355.7217公共設(shè)施管理業(yè)3 000.002 520.001 937.005 007.724 206.483 233.32合 計(jì)35 141 548.80161 766.4557 368.94110 586 768.1094 132.1751 628.67??由表2可知,石油、煉焦及核燃料加工業(yè)、廢棄資源和廢舊材料回收加工業(yè)、電子及通訊設(shè)備等行業(yè)排名靠前,而相對(duì)的紡織、冶金及化工等行業(yè)排名居末。
　　3.3.2 模型求解及結(jié)果分析
　　在Lingo軟件中編寫程序求解上述模型,得出優(yōu)化分配結(jié)果見表2。
　　由表2可以看出,應(yīng)用模型優(yōu)化后,COD排放總量削減了10%,新鮮用水量減少了41.81%,行業(yè)年總產(chǎn)值總和大幅度增長,增幅達(dá)到214.69%,并且,改善水環(huán)境質(zhì)量所獲得的長期社會(huì)綜合效益遠(yuǎn)不止于此。從上述優(yōu)化結(jié)論里也可以看出,蘇州市在今后行業(yè)結(jié)構(gòu)調(diào)整、環(huán)境政策制定能夠得到以上結(jié)果,是因?yàn)樗h(huán)境容量、水資源量被重新優(yōu)化配置給各行業(yè),區(qū)域行業(yè)結(jié)構(gòu)更為合理。
　　4 結(jié) 論
　　(1)1993-2006年的14年間,蘇州市工業(yè)COD排放量和排放強(qiáng)度均呈現(xiàn)下降趨勢,紡織業(yè)、化學(xué)原料及化學(xué)制品制造業(yè)、能源和水的生產(chǎn)與供應(yīng)業(yè)、造紙及木材加工、醫(yī)藥制造業(yè)等行業(yè)是蘇州市的COD重點(diǎn)排放行業(yè);
　　(2)經(jīng)多目標(biāo)行業(yè)總量優(yōu)化分配模型優(yōu)化后,COD排放總量削減了10%,新鮮用水量減少了41.81%,行業(yè)年總產(chǎn)值增幅達(dá)到214.69%,本模型在既定的限制條件(包括水資源與其它資源要素的供給情況等)下,為蘇州市的行業(yè)水污染物容量總量控制提供了一種可借鑒的思路和方法。
　　 (編輯:溫武軍)
　　參考文獻(xiàn)(References)
　　[1]張玉清. 河流功能區(qū)水污染物容量總量控制的原理和方法[M]. 北京:中國環(huán)境科學(xué)出版社, 2001.Zhang Yuqing. The Principle and Method of Total Water Pollutant Amount Control in Water Functional Areas [M]. Beijing: China Environmental Science Press, 2001.
　　[2]李如忠, 錢家忠, 汪家權(quán). 水污染物允許排放總量分配方法研究[J]. 水利學(xué)報(bào), 2003, (5): 112-115.Li Ruzhong, Qian Jiazhong, Wang Jiaquan. Study on Distribution of Total Amount of Drainage Water Pollutant in a Region [J]. Journal of Hydraulic Engineering, 2003,(5): 112-115.
　　[3]Cardwell, H. and Ellis, H. Stochastic Dynamic Programming Models for Water Quality Management. Water Resources Research, 1993, 29(40):803-813.
　　[4]Wen, C?G and Lee, C?S. Fuzzy Programming Approach to Water Quality Management Proceedings of the National Science Council, Republic of China. Part A: Physical Science and Engineering [Proc. Natl. Sci. Counc. Rep. China Pt. A: Phys. Sci. Eng.], 1998,22(5):579-590.
　　[5]Sasikumar K; Mujumdar PP. Fuzzy Optimization Model for Water Quality Management of a River System. Journal of Water Resources Planning and Management?ASCE. 1998, 124(2): 79-88.

　　[6]徐華君, 徐百福. 污染物允許排放總量分配的公平協(xié)調(diào)思路與方法[J]. 新疆大學(xué)學(xué)報(bào)(自然科學(xué)版), 1996, 13(3): 86-89.Xu Huajun, Xu Baifu. The Thought Way and Method of Fair Coordination on the Total Control [J]. Journal of Xinjiang University(Natural Science Edition), 1996, 13(3):86-89.
　　[7]汪俊啟, 張穎. 總量控制中水污染物允許排放量公平分配研究[J]. 安慶師范學(xué)院學(xué)報(bào)(自然科學(xué)版), 2000, 6 (3): 37-40.Wang Junqi, Zhang Yin. Study on National Distribution of Permissible Discharge Quantity of Water Pollutants in Total Quantity Control [J]. Journal of Anqing Teachers College (Natural Science Edition), 2000, 6(3): 37-40.
　　[8]林高松, 李適宇, 江峰. 基于公平區(qū)間的污染物允許排放量分配方法[J]. 水利學(xué)報(bào), 2006, 37(1): 52-57.Lin Gaosong, Li Shiyu, Jiang Feng. Allocation Method of Allowable Waste Loads in River Based on the Principle of Equitable Interval [J]. Journal of Hydraulic Engineering, 2006, 37(1): 52-57.
　　[9]方秦華, 張珞平, 洪華生. 水污染負(fù)荷優(yōu)化分配研究[J]. 環(huán)境保護(hù), 2005, (13): 29-31.Fang Qinhua, Zhang Luoping, Hong Huasheng. Review on Optimization of Water Pollution Loads Allocation [J]. Environmental Protection, 2005, (13): 29-31.
　　[10]王亮, 張宏偉, 岳琳. 水污染物總量行業(yè)優(yōu)化分配模型研究[J]. 天津大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版), 2006, 8(1): 59-63.Wang Liang, Zhang Hongwei, Yue Lin. Optimal Industry Total Distribution Model of Water Pollutants [J]. Journal of Tianjin University (Social Sciences), 2006, 8(1): 59-63.
　　[11]王亮. 天津市重點(diǎn)水污染物容量總量控制研究[D]. 天津: 天津大學(xué), 2005. WangLiang. Study on Total Capacity Control for Main Water Pollutants of Tianjin [D]. Tianjin: Tianjin University, 2005.
　　[12]熊俊. 基尼系數(shù)四種估算方法的比較與選擇[J]. 商業(yè)研究, 2003, (23): 123-125.XiongJun. The Comparison and Selection of Four Estimation Methods by Gini Coefficient [J]. Commercial Research, 2003, (23): 123-125.
　　[13]葉禮奇. 基尼系數(shù)計(jì)算方法[J]. 中國統(tǒng)計(jì), 2003, (4): 58.Ye Liqi. The Gini Coefficient Calculation Methods [J]. China Statistics, 2003, (4): 58.?
　　Study on Optimal Distribution of Total COD of City Industry?
　　YANG Zhan?hong?1 LUO Hong?1 LV Lian?hong?1 XU Yong?li?2
　　(1. Chinese Research Academy of Environmental Sciences, Beijing 100012, China; ?2. Environmental Geology Perambulation Academy of Hebei Province, Shijiazhuang Hebei 050021,China)?
　　Abstract Total amount control is an important approach to control the pollution, improve the environment quality and realize sustainable development. How to find a optimal model to balance environment, economy, technology and resource under a relatively fair condition is a research subject in environmental science. This articlestudies the emission characteristics of industry COD in Suzhou city, and uses Gini coefficient to evaluate the equity of COD emission distribution of by. Meanwhile, this study relates the total amount control with resource, society and economic, aiming to maximize the economic benefits andminimize the investment cost of pollution control. The multi-object optimal industry total distribution model is used to redistribute the total amount of COD among industries of Suzhou. The results show that textile industry, chemical materials and chemical products industry, energy and water production and supply industry, papermaking and woo d processing industry and medicine industry are primary COD emission industries. After optimal distribution of total COD, the rate of COD emission decreased 10%, fresh water consumption decreased 41.81% and industry annual total output value increased 214.69%. Also, resources and water environmental capacity could realize the optimal allocation under the Correspondingspeed of economic growth. Total amount control may create the largest potential environmental and economic benefitsunder market economy system.
　　Key words Chemical Oxygen Demand(COD); equity; Gini coefficient; optimal distribution

——本文由豐臨科技整理發(fā)布，內(nèi)容供參考，如有侵權(quán)，請(qǐng)聯(lián)系刪除，謝謝！上海豐臨科技有限公司為你提供濁度儀(濁度計(jì))、在線濁度儀、余氯儀、余氯分析儀、工業(yè)在線pH計(jì)、cod測定儀、pH計(jì)等多種水質(zhì)檢測儀，水質(zhì)分析儀，歡迎您前來選購，豐臨科技竭誠為您服務(wù)！