張浩，尚大軍，樊霄雁，劉智勇

（1.河北工業(yè)大學化工學院，天津300130；2.中國石油大慶石化分公司，黑龍江大慶163711；3.河北工業(yè)大學海洋科學與工程學院，天津300130）

單雜質(zhì)整體水網(wǎng)絡(luò)設(shè)計新方法

張浩1，尚大軍2，樊霄雁3，劉智勇3

（1.河北工業(yè)大學化工學院，天津300130；2.中國石油大慶石化分公司，黑龍江大慶163711；3.河北工業(yè)大學海洋科學與工程學院，天津300130）

研究了單雜質(zhì)整體水網(wǎng)絡(luò)設(shè)計.首先分析了各子網(wǎng)絡(luò)之間的關(guān)系，指出具有再生單元的用水子系統(tǒng)與廢水處理子系統(tǒng)之間的聯(lián)系是高濃度廢水.在此基礎(chǔ)上將設(shè)計過程分為2步：首先設(shè)計包含再生單元用水子網(wǎng)絡(luò)，得到廢水流量和濃度，第2步以廢水處理量最小為目標設(shè)計廢水處理網(wǎng)絡(luò).在具有再生單元的子系統(tǒng)中，再生過程模型分為固定雜質(zhì)移除率（RR）及固定再生后濃度2種形式.對于固定移除率模型問題，首先估算出再生濃度，通過幾次迭代即可得到再生子網(wǎng)絡(luò)設(shè)計結(jié)果；對于固定再生后濃度模型問題無需迭代，只需1次計算即可得到再生子網(wǎng)絡(luò)設(shè)計結(jié)果.最后，通過2個實例來說明本文方法.結(jié)果表明，本文方法的設(shè)計過程簡單且有效.

水系統(tǒng)集成；單雜質(zhì)；整體水網(wǎng)絡(luò)；固定移除率；固定再生后濃度；廢水處理

0 引言

工業(yè)節(jié)水非常重要.傳統(tǒng)水網(wǎng)絡(luò)設(shè)計通常將用水子網(wǎng)絡(luò)和廢水處理子網(wǎng)絡(luò)分開研究.而包含用水子系統(tǒng)、再生子系統(tǒng)及污水處理子系統(tǒng)的整體水網(wǎng)絡(luò)，全面考慮廢水回用[1]、再生再利用/循環(huán)利用[2-3]、廢水處理[4-5]諸過程，可以實現(xiàn)最大程度上的節(jié)水減排，所以整體水網(wǎng)絡(luò)研究對于工業(yè)節(jié)水具有重要意義.上述子系統(tǒng)的設(shè)計均較為較復雜，而上述子系統(tǒng)構(gòu)成的整體水網(wǎng)絡(luò)的設(shè)計則更為復雜.

目前水網(wǎng)絡(luò)的研究方法主要有夾點法與數(shù)學規(guī)劃法.Takama等[6]首次采用數(shù)學規(guī)劃法對煉油廠中水優(yōu)化分配進行研究.Wang和Smith[1]提出水夾點技術(shù)研究廢水最小化問題，通過構(gòu)造極限濃度復合曲線確定水夾點的位置，由此得到新鮮水最小量.Doyle和Smith[7]采用數(shù)學規(guī)劃法研究了多雜質(zhì)水網(wǎng)絡(luò)的回用問題.對于廢水處理網(wǎng)絡(luò)，Wang和Smith[8]開發(fā)了圖形化方法確定廢水處理夾點，并由此得到廢水處理量目標值.Galan等[9]采用連續(xù)性松弛解法研究了分布式廢水處理網(wǎng)絡(luò).目前文獻中對用水網(wǎng)絡(luò)及廢水處理網(wǎng)絡(luò)的研究還有很多，限于篇幅，本文主要討論整體水網(wǎng)絡(luò)設(shè)計.

已知文獻報道中，對再生單元為固定出口濃度模型的簡單總體水網(wǎng)絡(luò)主要用夾點法，如Ng等[10-11]采用夾點圖與水級聯(lián)分析技術(shù)設(shè)計單雜質(zhì)整體水網(wǎng)絡(luò).Deng等[12]提出PGA（Process-based Graphical Approach）法設(shè)計固定雜質(zhì)負荷與固定流量整體水網(wǎng)絡(luò).Parand等[13]提出ECTA（Extended Composite Table Algorithm）法分別求解總水網(wǎng)絡(luò)中新鮮水消耗量、再生水流量、廢水流量和廢水濃度.對于復雜的再生單元為固定出口濃度模型的總體水網(wǎng)絡(luò)的設(shè)計，許多學者采用數(shù)學規(guī)劃法，如Putra和Amminudin[14]提出了一種由MILP和NLP構(gòu)成的2步法設(shè)計整體水網(wǎng)絡(luò).Tudor和Lavric[15]采用GA（Genetic Algorithm）法求解有進、出口濃度限制的整體水網(wǎng)絡(luò)數(shù)學模型.Ahmetovic和Grossmann[16]基于各類過程的超結(jié)構(gòu)模型得到非凸NLP和MINLP模型，推導出變量的嚴格邊界，由此得到整體水網(wǎng)絡(luò)全局最優(yōu)解.Kim[17]對總水網(wǎng)絡(luò)集成進行了系統(tǒng)的分析，采用圖形方法和數(shù)學優(yōu)化法結(jié)合的策略考察了新鮮水費用、污水處理費用、管道費用和廢水排放費用之間的關(guān)系.Khor等[18]采用全局優(yōu)化求解器對固定流量整體水網(wǎng)絡(luò)集成問題進行了研究.

上述文獻研究再生單元為固定再生后濃度模型的整體水網(wǎng)絡(luò)，而對于RR模型的水網(wǎng)絡(luò)設(shè)計方法并不多.Kuo和Smith[19]改進圖解法[1]確定RR模型用水網(wǎng)絡(luò)的最小新鮮水消耗量，但過程較為繁瑣.Bandyopad hyay等[20]采用了源水流復合曲線（Source Composite Curve）法設(shè)計RR模型的單雜質(zhì)水網(wǎng)絡(luò).Pan等[21]提出迭代法設(shè)計考慮再生再利用/循環(huán)過程的單雜質(zhì)用水網(wǎng)絡(luò)，對再生單元為RR模型與固定再生后濃度模型的水網(wǎng)絡(luò)設(shè)計都能有效解決.已知文獻研究對再生單元為RR模型的整體水網(wǎng)絡(luò)研究尚未見報道.

本文首先分析了整體水網(wǎng)絡(luò)中各子系統(tǒng)之間的聯(lián)系，指出具有再生單元的用水子系統(tǒng)與廢水處理子系統(tǒng)之間的聯(lián)系是高濃度廢水.在此基礎(chǔ)上將單雜質(zhì)整體水網(wǎng)絡(luò)的設(shè)計分為2個步驟：1）首先設(shè)計包含再生單元的用水網(wǎng)絡(luò)，同時得出待處理廢水流的量和濃度；2）在上述基礎(chǔ)上設(shè)計廢水處理網(wǎng)絡(luò).這樣就可以簡化設(shè)計過程，由設(shè)計實例可知，本文方法設(shè)計步驟簡單且有效.

1 整體水網(wǎng)絡(luò)分析

整體水網(wǎng)絡(luò)由用水子系統(tǒng)、再生子系統(tǒng)及污水處理子系統(tǒng)3部分組成，如圖1所示.

圖1 整體水網(wǎng)絡(luò)各部分間的關(guān)系Fig.1 Interactions between elements of a total water network

從用水子系統(tǒng)出來的未被利用的低濃度水流進入再生子系統(tǒng)，再生后水流返回用水子系統(tǒng)重新利用.而用水子系統(tǒng)中未被利用的高濃度水流則進入廢水處理系統(tǒng)進行處理，滿足環(huán)境排放要求后排放.

由圖1可見，包含再生過程的用水子網(wǎng)絡(luò)（圖1中虛線框內(nèi)部分）與廢水處理子網(wǎng)絡(luò)之間的聯(lián)系是高濃度廢水.這樣可將整體水網(wǎng)絡(luò)設(shè)計分為2步，即先設(shè)計包含再生單元的用水網(wǎng)絡(luò)，由此可得到待處理的廢水流的流量及濃度.然后以所得廢水流數(shù)據(jù)為基礎(chǔ)設(shè)計廢水處理子網(wǎng)絡(luò).采用上述設(shè)計策略可以簡化總體設(shè)計難度.設(shè)計過程如圖2所示.

2 含再生再利用/循環(huán)利用過程的用水網(wǎng)絡(luò)設(shè)計

Liu等[22]指出包含再生單元的水網(wǎng)絡(luò)與只考慮再利用的水網(wǎng)絡(luò)的區(qū)別在于前者增加了1條（或幾條）再生水流，如果確定了再生水流的濃度與流量，并把再生水流作為新增源水流加入到只包含再利用過程的水網(wǎng)絡(luò)當中，即可得到包含再生單元的水網(wǎng)絡(luò)所需數(shù)據(jù)，然后可采用對只考慮再利用水網(wǎng)絡(luò)提出的設(shè)計方法進行設(shè)計[22].

對再生單元采用固定雜質(zhì)移除率模型的水網(wǎng)絡(luò)，劉智勇及其合作者[21]在Liu等[22]工作的基礎(chǔ)上提出迭代方法設(shè)計包含再生單元的用水網(wǎng)絡(luò)，首先估算出再生水流濃度，假設(shè)再生水流的量足夠大，將再生水流作為新增源水流加入到只考慮再利用的水網(wǎng)絡(luò)中，即可形成具有再生單元的用水網(wǎng)絡(luò)，經(jīng)過幾次迭代，當兩次相鄰迭代再生濃度的差小于某一給定值ε，即可得到最終設(shè)計.對再生單元采用固定出口濃度模型時，采用Liu等[23]提出的多水源法，只需1次迭代即可得到具有再生單元的網(wǎng)絡(luò)設(shè)計.包含再生單元的用水網(wǎng)絡(luò)設(shè)計完成后，可得到廢水流的流量和濃度，然后進行廢水處理子網(wǎng)絡(luò)設(shè)計.

圖2 單雜質(zhì)總水網(wǎng)絡(luò)設(shè)計過程Fig.2 Design procedure for a total water network of single contaminant

3 廢水處理網(wǎng)絡(luò)設(shè)計

從上述用水網(wǎng)絡(luò)出來的未被利用的廢水可分為2部分，滿足環(huán)境排放要求的水流可以直接排放，高濃度的水流則需要經(jīng)過處理單元處理.

McLaughlin等[24]指出，雜質(zhì)負荷一定的情況下，水流濃度越低則處理量就越高，這樣會使處理費用增大.本文采用時靜[25]提出的單雜質(zhì)分布式水處理方法設(shè)計廢水處理網(wǎng)絡(luò)，首先將廢水流按照濃度降低的順序排列，然后確定廢水流中夾點水流處理量，通過合理設(shè)計處理過程，最大程度地減少不必要的水流混合來降低系統(tǒng)的總處理量.

4 實例研究

本節(jié)研究了2個實例.

4.1 例1

本例數(shù)據(jù)見表1，數(shù)據(jù)取自文獻[26]，規(guī)定再生單元固定移除率（RR值）為95%，處理單元固定移除率（RR值）為99%，環(huán)境極限排放濃度為50 g/t.

表1 例1數(shù)據(jù)Tab.1 The data of example 1

設(shè)計過程如下.

4.1.1 采用迭代法[21]設(shè)計包含再生單元的用水網(wǎng)絡(luò)

由文獻的設(shè)計步驟可得：首先估算再生濃度，將各源水流合并，計算再生濃度為12.45 g/t，再生水作為新增源水流加入水網(wǎng)絡(luò)中用于第1次迭代過程.執(zhí)行過程中入口極限濃度較小的過程優(yōu)先執(zhí)行，當入口極限濃度相同時，則出口濃度較小的過程優(yōu)先執(zhí)行.入口極限濃度最低為0 g/t的需求水流（D1，D2及D3）只能用新鮮水（FW）滿足，用量為49.17 t/h.

按照濃度順序，此時應該執(zhí)行P4過程.對于P4過程，因源水流S1、S2、S3濃度等于或大于P4過程出口極限濃度，所以不應回用源水流S1、S2、S3，而只能使用再生水Sreg，再生水量為29.61 t/h，P4過程出口濃度達到極限值.對于P5過程，可用源水流為S1、S2、S3（S1濃度最小優(yōu)先回用），根據(jù)文獻[21]方法計算得，源水流S1對于P5過程分配量為8.23 t/h，再生源水流Sreg對P5過程分配量為36.08 t/h，這樣過程入口與出口濃度均達到極限值.

類似地，對后續(xù)過程進行設(shè)計，直至所有過程執(zhí)行完成.第1次迭代得到再生水流量為272.83 t/h.按照高濃度源水流排放，低濃度源水流再生的規(guī)則（第1次迭代需再生處理的水流如表2所示），得到該過程再生濃度用于第2次迭代過程，直至第4次迭代再生水流濃度不再變化，每次迭代的結(jié)果見表3.通過比較可以看出，第3、4次的迭代結(jié)果相同，說明本例可以通過第3次迭代即可得到最終設(shè)計結(jié)果.

表2 例1第1次迭代需再生處理的水流數(shù)據(jù)Tab.2 The data of regeneration stream for iteration 1

表3 例1迭代結(jié)果Tab.3 Iteration results for example 1

最終的設(shè)計結(jié)果新鮮水消耗量與廢水排放量都為46.25 t/h，再生水量為276.73 t/h，再生水濃度為14.38 g/t.各過程水流分配量如表4所示.

表4 例1各過程水流分配數(shù)據(jù)Tab.4 The allocation data for example 1

4.1.2 廢水處理網(wǎng)絡(luò)設(shè)計

處理過程的脫除率為99%，環(huán)境極限濃度為50 g/t，環(huán)境允許排放的最大雜質(zhì)負載量為2.46 kg/h.

根據(jù)單雜質(zhì)分布式水處理方法[25]，把廢水流股按照濃度下降的順序進行排列，并計算出各水流雜質(zhì)負荷量（M），經(jīng)過處理后的剩余雜質(zhì)負荷量以及雜質(zhì)負荷累積量如表5所示.

表5 廢水流數(shù)據(jù)Tab.5 Wastewater data for example 1

根據(jù)文獻[25]，從表5的值可以看出，所以判斷出水流3為夾點水流.濃度濃度高于水流3的所有水流需要完全處理，水流3部分處理部分旁路，濃度低于水流3的水流則直接排放.由計算可知，夾點水流3的處理量為32.05 t/h，總處理量為49.17 t/h.最終的設(shè)計如圖3所示.

4.2 例2

本例數(shù)據(jù)見表6，數(shù)據(jù)取自文獻[11]，再生單元采用固定出口濃度模型，出口濃度為10 g/t，處理單元固定移除率（RR值）為90%，環(huán)境極限排放濃度為20 g/t.

假設(shè)新鮮水與再生水流量充足，再生水作為新增源水流加入水網(wǎng)絡(luò)中，將各水流按濃度升高的順序排列見表7.

圖3 例1最終網(wǎng)絡(luò)結(jié)構(gòu)Fig.3 Final design for example 1

表6 例2數(shù)據(jù)Tab.6 The data of example 2

設(shè)計過程如下.

4.2.1 采用多水源法[23]設(shè)計包含再生單元的用水網(wǎng)絡(luò)

濃度為0 g/t的需求水流（D1）只能用新鮮水來滿足，用量為120 t·h-1.

由文獻[23]設(shè)計步驟可得：對于第1個濃度為50 g/t的需求水流可以由高濃度源水流（HCS為100 g/t）和低濃度源水流（LCS為10 g/t）來滿足，其中高濃度源水流HCS用量為35.56 t/h，低濃度源水流LCS用量為44.44 t/h.類似地可以計算確定各個用水過程水流分配量，具體結(jié)果見表8.

由表中數(shù)據(jù)計算可得再生水流量為88.88 t/h，各水流剩余量分別為：100 g/t流量8.88 t/h，180 g/t流量35 t/h，250 g/t流量85 t/h.低濃度水流可以用于再生，所以再生水流由濃度為100 g/t、180 g/t的源水流及45 t/h濃度為250 g/t的水流組成.剩余40 t/h濃度為250 g/t的水流進入處理單元.

表7 例2數(shù)據(jù)Tab.7 The data of example 2

表8 例2各過程水流分配數(shù)據(jù)Tab.8 The allocation data for example 2

4.2.2 廢水處理網(wǎng)絡(luò)設(shè)計

進入廢水處理單元的只有一股水流，根據(jù)質(zhì)量守恒確定第一個處理單元處理量10 t/h，第二個處理單元處理量39.89 t/h，最終網(wǎng)絡(luò)結(jié)構(gòu)見圖4.

圖4 例2最終網(wǎng)絡(luò)結(jié)構(gòu)Fig.4 Final design for Example 2

5 結(jié)論

本文提出了設(shè)計單雜質(zhì)整體水網(wǎng)絡(luò)的新方法.首先分析了整體水網(wǎng)絡(luò)中各子系統(tǒng)之間的關(guān)系，在此基礎(chǔ)上將設(shè)計過程分為2步：先設(shè)計包含再生單元的用水網(wǎng)絡(luò)，再設(shè)計廢水處理網(wǎng)絡(luò).本文研究的再生系統(tǒng)包括固定再生出口濃度模型及固定雜質(zhì)移除率模型2種情況.其中包含固定雜質(zhì)移除率再生過程的整體水網(wǎng)絡(luò)尚未見諸文獻報道.對文獻實例研究結(jié)果表明：本文方法的設(shè)計過程簡單有效，具有明確的工程意義.

[1]Wang Y P，Smith R.Wastewater minimisation[J].Chemical Engineering Science，1994，49（94）：981-1006.

[2]Alva-Argáez A，Kokossis A C，Smith R.Wastewater minimisation of industrial systems using an integrated approach[J].Computers and Chemical Engineering，1998，22（12）：S741-S744.

[3]Olesen S G，Polley G T.A simple methodology for the design of water networks handling single contaminants[J].Chemical Engineering Research and Design，1997，75（4）：420-426.

[4]Kuo W C J，Smith R.Designing for the interactions between water-use and effluent treatment[J].Chemical Engineering Research and Design，1998，76（3）：287-301.

[5]Wang Y P，Smith R.Wastewater minimization with flowrate constraints[J].Chemical Engineering Research and Design，1995，73：889-904.

[6]Takama N，Kuriyama T，Shiroko K，et al.Optimal water allocation in a petroleum refinery[J].Computers and Chemical Engineering，1980，4（4）：251-258.

[7]Doyle S J，Smith R.Targeting water reuse with multiple contaminants[J].Process safety and environmental protection，1997，75（3）：181-189.

[8]Wang Y P，Smith R.Design of distributed effluent treatment systems[J].Chemical Engineering Science，1994，49（18）：3127-3145.

[9]Galan B，Grossmann I E.Optimal design of distributed wastewater treatment networks[J].Industrial and Engineering Chemistry Research，1998，37（10）：4036-4048.

[10]Ng D K S，F(xiàn)oo D C Y，Tan R R.Targeting for total water network 1.waste stream identification[J].Industrial and Engineering Chemistry Research，2007，46（26）：9107-9113.

[11]Ng D K S，F(xiàn)oo D C Y，Tan R R.Targeting for total water network.2.waste treatment targeting and interactions with water system elements[J].Industrial and Engineering Chemistry Research，2007，46（46）：9114-9125.

[12]Deng C，F(xiàn)eng X，Ng D K S，et al.Process-based graphical approach for simultaneous targeting and design of water network[J].Aiche Journal，2011，57（57）：3085-3104.

[13]Parand R，Hong M Y，Pareek V，et al.Use of pinch concept to optimize the total water regeneration network[J].Industrial and Engineering Chemistry Research，2014，53（8）：3222-3235.

[14]Putra Z A，Amminudin K A.Two-step optimization approach for design of a total water system[J].Industrial and Engineering Chemistry Research，2008，47（16）：6045-6057.

[15]Tudor R，Lavric V.Optimization of total networks of water-using and treatment units by genetic algorithms[J].Industrial and Engineering Chemistry Research，2010，49（8）：3715-3731.

[16]Ahmetovi E，Grossmann I E.Global superstructure optimization for the design of integrated process water networks[J].Aiche Journal，2011，57（2）：434-457.

[17]Kim J K.System analysis of total water systems for water minimization[J].Chemical Engineering Journal，2012，193-194（12）：304-317.

[18]Cheng S K，Chachuat B，Shah N.Fixed-flowrate total water network synthesis under uncertainty with risk management[J].Journal of Cleaner Production，2014，77：79-93.

[19]Kuo W C J，Smith R.Design of water-using systems involving regeneration[J].Process Safety and Environmental Protection，1998，76（2）：94-114.

[20]Bandyopadhyay S，Cormos C C.Water management in process industries incorporating regeneration and recycle through a single treatment Unit[J].Industrial and Engineering Chemistry Research，2008，47（4）：1111-1119.

[21]潘春暉，時靜，劉智勇.采用迭代方法設(shè)計具有再生再利用/循環(huán)的單雜質(zhì)用水系統(tǒng)[J].化工學報，2011，62（5）：1339-1344.

[22]Liu Z Y，Li Y M，Liu Z H，et al.A simple method for design of water-using networks with multiple contaminants involving regeneration reuse[J]. Aiche Journal，2010，55（55）：1628-1633.

[23]劉智勇，李艷梅，張廣林，等.Simultaneously designing and targeting for networks with multiple resources of different Qqualities[J].中國化學工程學報：英文版，2009，17（3）：445-453.

[24]McLaughlin L.Develop an effective wastewater treatment strategy[J].Chemical Engineering Progress，1992，88（9）：35-42.

[25]時靜.分布式廢水處理體系的設(shè)計研究[D].天津：河北工業(yè)大學，2011.

[26]Savelski M J，Bagajewicz M J.Algorithmic procedure to design water utilization systems featuring a single contaminant in process plants[J].Chemical Engineering Science，2001，56（5）：1897-1911.

[責任編輯 田豐]

A new method for design of total water networks of single contaminant

ZHANG Hao1,SHANG Dajun2,FAN Xiaoyan3,LIU Zhiyong3
（1.School of Chemical Engineering and Techology，Hebei University of Technology，Tianjin 300130，China；2.Refinery of Daqing Petrochemical Company，Petrochina，Heilongjiang Daqing 163711，China；3.School of Marine Science and Engineering，Hebei University of Technology，Tianjin 300130，China）

This paper studies the design of total water networks of single contaminant.Based on the analysis of interactions between the sub-networks of the system,the link between water using sub-network involving regeneration unit and effluent treatment sub-network is the high concentration wastewater.A two-stage design procedure was proposed for a design of total water networks.In the first stage,the waterusing sub-network involving regeneration unit was designed,and the concentration and the amounts of effluent streams were obtained.In the second stage,the wastewater treatment subnetwork was designed by minimizing unnecessary stream mixing.Regeneration units with models offixed removal ratio (RR)and fixed regeneration concentration can be used in the networks.For the model with fixed RR,the initial concentration(s)of the Sregshould be estimated first.The final design can be obtained after several iterations.For the model of fixed regenerated concentration,the final design can be obtained with one iteration.Two examples were chosen to illustrate the presented method.Results have shown that the design procedure proposed in this paper is simple and effective.

water system integration;single contaminant;total water network;fixed removal ratio;fixed regeneration concentration;wastewater treatment

TQ085

A

1007-2373（2017）02-0085-07

10.14081/j.cnki.hgdxb.2017.02.015

2017-02-18

國家自然科學基金（21176057）；國家重點基礎(chǔ)研究發(fā)展計劃（2012CB720305）

張浩（1990-），男，碩士研究生.通訊作者：劉智勇（1956-），男，教授，博士生導師.