寧東能源化工基地水資源優(yōu)化配置研究

馮興平，江 濤，唐 林，李永紅*，范 磊，趙正宏，侯光才，錢 會

(1.西北農(nóng)林科技大學(xué)水土保持研究所，陜西楊凌 712100；2.長安大學(xué)環(huán)境科學(xué)與工程學(xué)院，陜西西安710054；3.西安市城市規(guī)劃設(shè)計研究院，陜西西安 710082；4.中國地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心，陜西西安 710054)

寧東能源化工基地位于寧夏回族自治區(qū)中東部，在靈武、鹽池、同心、紅寺堡地區(qū)，南北長130 km，東西寬50 km，面積約為3 500 km2，已探明煤炭儲量273×108t，遠(yuǎn)景儲量1 394.3×108t，是一個全國罕見的儲量大、煤質(zhì)好、地質(zhì)構(gòu)造簡單的整裝煤田，被列為國家13個重點(diǎn)開發(fā)的億噸級礦區(qū)之一。寧東能源化工基地位于鄂爾多斯盆地西緣靈鹽臺地之上，區(qū)內(nèi)水資源極為貧乏。基地供水來源單一，絕大部分供水來源于黃河引水，少量來源于吳忠平原地下水。黃河引水不但用水成本很高，而且地表水易受污染，輸水具有季節(jié)性，保障程度低。因此，如何高效利用水資源是寧東能源化工基地建設(shè)面臨的重要問題，受到社會各界的廣泛關(guān)注。目前文獻(xiàn)的研究，只是給出了基地各類水(黃河水、礦井疏干水、處理污水)配置數(shù)量，并提出礦井水、各類廢水處理利用及節(jié)約用水解決基地供水的措施[1-2]，但未涉及基地地表水與地下水的聯(lián)合運(yùn)用以及水資源優(yōu)化配置問題，因而前人的研究對基地合理高效開發(fā)利用水資源的實(shí)際指導(dǎo)作用尚有欠缺，有必要進(jìn)行深入研究。筆者依據(jù)“鄂爾多斯盆地(寧夏)能源基地地下水勘查”項(xiàng)目的地下水資源評價資料，結(jié)合寧東能源化工基地建設(shè)規(guī)劃對水資源的需求，提出了地表水和地下水聯(lián)合供水方案，并采用線性規(guī)劃模型對供水方案進(jìn)行了優(yōu)化研究，為基地水資源經(jīng)濟(jì)和高效利用提供基礎(chǔ)依據(jù)，也為干旱區(qū)工業(yè)基地經(jīng)濟(jì)合理利用水資源提供參考。

1 寧東能源化工基地產(chǎn)業(yè)規(guī)劃對水資源的需求及供水現(xiàn)狀

1.1寧東能源化工基地產(chǎn)業(yè)規(guī)劃按照寧東能源化工基地總體規(guī)劃，基地以煤炭、電力和煤化工為三大主導(dǎo)產(chǎn)業(yè)，輔助發(fā)展石油天然氣化工、精細(xì)化工、材料工業(yè)和其他服務(wù)業(yè)。煤礦區(qū)包括靈武、鴛鴦湖、橫城3個礦區(qū)，12個井田；電廠包括馬連臺電廠、靈州煤矸石綜合電廠、靈武電廠、鴛鴦湖電廠、水洞溝電廠、方家莊電廠和棗泉電廠等；3個工業(yè)園區(qū)包括煤化工園區(qū)、臨河綜合項(xiàng)目區(qū)和靈州綜合項(xiàng)目區(qū)。在寧夏東北部與內(nèi)蒙古交界地帶建設(shè)了平羅精細(xì)化工基地[3-4]。寧東能源化工基地地理位置及主要廠礦的分布見圖1。

圖1 寧東能源化工基地廠礦分布Fig.1 Distribution of factories and mines in Ningdong energy chemical base

1.2寧東基地對水資源的需求主要包括煤礦生產(chǎn)生活、電廠生產(chǎn)生活和煤化工項(xiàng)目對新鮮淡水的需求。根據(jù)《寧東基地總體規(guī)劃與建設(shè)綱要》，2020年總需水量達(dá)到3.43×108m3/a，其中煤礦生產(chǎn)生活用水0.16×108m3/a，電廠生產(chǎn)生活用水0.29×108m3/a，煤化工生產(chǎn)用水2.98×108m3/a，以煤化工用水量最大。環(huán)保綠化需水可通過礦坑廢水和工業(yè)廢水的處理或再循環(huán)利用解決。

1.3寧東基地引黃供水工程和應(yīng)急水源需求

1.3.1引黃供水工程。寧夏已經(jīng)按照基地規(guī)劃修建了2座引黃供水水庫：鴨子蕩水庫和劉家溝水庫。其中，鴨子蕩水庫位于基地北部，2020年供水量可以達(dá)到3.00×108m3/a以上；太陽山供水工程(劉家溝水庫)位于基地南部，供水量為0.36×108m3/a。這2座水庫的水源均來自黃河，通過專門的管道進(jìn)行輸送，水量的配額來自水權(quán)置換，即把農(nóng)業(yè)節(jié)約用水所產(chǎn)生的余額分配給工業(yè)用水[5]。

1.3.2應(yīng)急水源需求。應(yīng)急水源首先保障煤礦生產(chǎn)生活、電廠生產(chǎn)生活新鮮淡水的需求，即2020年煤礦生產(chǎn)生活用水0.16×108m3/a，以及電廠生產(chǎn)生活用水0.29×108m3/a，合計0.45×108m3/a。煤化工項(xiàng)目的用水需求暫且不納入應(yīng)急供水保障范圍。

2 研究資料與方法

筆者依據(jù)“鄂爾多斯盆地(寧夏)能源基地地下水勘查”項(xiàng)目的地下水資源評價資料(2014年)，通過統(tǒng)計和圖表方法分析可向?qū)帠|能源化工基地供水的各地下水水源地的可采資源量及基地周邊水源地應(yīng)急供水能力；利用拓?fù)鋱D方式表達(dá)地表水與地下水聯(lián)合供水的方案的基本內(nèi)容；采用線性規(guī)劃模型對聯(lián)合供水方案進(jìn)行優(yōu)化研究，確定供水方案中地下水的實(shí)際配額和供水路線。

3 結(jié)果與分析

3.1地下水源地及其供水潛力目前，已經(jīng)探明可以向?qū)帠|能源化工基地供水的地下水水源地有5處，各水源地的地下水可開采量統(tǒng)計如表1所示。備選應(yīng)急水源地的可開采量匯總見表2。將上述水源地的供水能力繪制成框圖(圖2)。

表1寧東能源化工基地水源地地下水可開采量統(tǒng)計

Table1StatisticsofgroundwaterrecoverableinNingdongenergychemicalbase

水源地Water head site現(xiàn)狀開采量Current exploitation quantity×108 m3/a工程可開采量Engineering recoverable quantity×108 m3/a水質(zhì)Water Quality駱駝井水源地Luotuojing water head site0.0180.046淡水靈武大泉水源地Lingwu Daquan water head site0.0900.270淡水靈武崇興水源地Lingwu Chongx-ing water head site0.0400.070淡水陶樂平原水源地Taole plain water head site0.0100.113淡水吳忠平原傍河水源地Wuzhong plain riverside water head site00.246淡水

表2 寧東基地備選應(yīng)急水源地的可開采量匯總

注：括號數(shù)據(jù)中為已有開采量(×108 m3/a) Note：The data in parentheses is the existing mining amount(×108 m3/a) 圖2 寧東能源基地周邊水源地應(yīng)急供水能力框圖Fig.2 Block diagram of emergency water supply capacity of water sources near Ningdong energy base

根據(jù)圖2分析，來自寧東能源化工基地西部的鹽池駱駝井地區(qū)可供水0.046×108m3/a，東部的靈武地區(qū)可供水0.465×108m3/a，而北部的陶樂地區(qū)可供水0.113×108m3/a。合計的工程可開采量為0.624 ×108m3/a，已經(jīng)能滿足寧東能源化工基地2020年應(yīng)急供水0.450×108m3/a的要求。更遠(yuǎn)期的應(yīng)急供水可能需要利用微咸水，要考慮對微咸水進(jìn)行處理。

3.2地表水與地下水聯(lián)合供水的方案按照寧東能源基地的總體規(guī)劃，能源化工企業(yè)的用水主要通過鴨子蕩水庫和劉家溝水庫提供[6]。引黃供水工程依賴于黃河及其附屬的水利設(shè)施，輸水具有季節(jié)性，需要修建水庫或大型蓄水池并處理泥沙沉淀，供水和維護(hù)運(yùn)營的成本很高，蒸發(fā)滲漏造成的水資源浪費(fèi)較大，而且地表水易受污染，保障程度低[7]。在探明了地下水水源地之后，通過適當(dāng)引入地下水的供水工程，可以減少蒸發(fā)損失，降低運(yùn)營成本，并起到應(yīng)急供水的作用。但是，僅靠現(xiàn)有的地下水水源地，還不能滿足能源化工基地的全部新鮮淡水需求。因此，采用地表水與地下水聯(lián)合供水的方案，將更有利于區(qū)域水資源的高效開發(fā)利用[8]。經(jīng)分析研究，提出了如圖3所示的寧東能源化工基地地表水與地下水聯(lián)合供水方案。方案的要點(diǎn)如下：

(1)陶樂水源地向平羅精細(xì)化工基地單獨(dú)供水。平羅精細(xì)化工基地在陶樂地區(qū)的北端，目前的供水措施是引黃河水，經(jīng)供水管道和泵站輸送到工業(yè)區(qū)。該基地規(guī)劃供水為2.0萬m3/d(0.073×108m3/a)，而陶樂傍河水源地的可開采量達(dá)到0.113×108m3/a，能夠滿足平羅精細(xì)化工基地的需求。

(2)修建輸水管道和配水站聯(lián)合駱駝井應(yīng)急水源地，向橫城工礦區(qū)和鴛鴦湖礦區(qū)供水。駱駝井應(yīng)急水源地的供水能力為0.046×108m3/a。

(3)在靈武地區(qū)，修建配水站和輸水管道聯(lián)合傍河3號水源地和鴨子蕩水庫的供水管道，向橫城工礦區(qū)供水。另外修建配水站聯(lián)合現(xiàn)有的崇興水源地和2號備選水源地聯(lián)合鴨子蕩水庫的供水管道，向靈武工礦區(qū)提供煤礦和電廠的生產(chǎn)生活用水。這種與鴨子蕩水庫聯(lián)合調(diào)度的方案中，地下水的可供水量達(dá)到0.276×108m3/a。鴨子蕩水庫則向煤化工基地提供煤化工項(xiàng)目的新鮮淡水需求，2020年設(shè)計供水能力為2.980×108m3/a，地下水可供水量占比很小。

(4)在靈武南部地區(qū)，將大泉水源地與1號備選水源地與劉家溝水庫的地表水進(jìn)行聯(lián)合調(diào)度，向積家井和馬家灘工礦區(qū)供水。馬家灘和積家井礦區(qū)2020年規(guī)劃需要0.054×108m3/a的新鮮淡水，低于劉家溝水庫的供水能力。大泉水源地和1號水源地的可供水總量為0.189×108m3/a，比馬家灘和積家井礦區(qū)需求量還大，因此剩余的水量可以向北分配到靈武工礦區(qū)。

圖3 寧東能源化工基地地表水與地下水聯(lián)合供水方案拓?fù)鋱DFig.3 Topological plan of joint water supply plan for surface water and underground water at Ningdong energy chemical base

3.3寧東能源化工基地供水方案的優(yōu)化以上的討論提出了寧東能源化工基地地表水與地下水聯(lián)合供水的方案，確定了周邊地下水水源地的供水線路和最大可能供水規(guī)模，但還存在很大的不確定性，即地下水的實(shí)際供水規(guī)模達(dá)到多大才是經(jīng)濟(jì)合理的方案。為了定量地確定方案中的供水配額，該研究通過線性規(guī)劃模型進(jìn)行分析。

3.3.1模型方程。寧東能源化工基地地表水與地下水聯(lián)合供水方案優(yōu)化的基本思路是用最低的經(jīng)濟(jì)成本達(dá)到供水目標(biāo)[9-10]。因?yàn)橐命S河水的供水方案是已經(jīng)有規(guī)劃的，工程建設(shè)已經(jīng)部分完成。因此，在該模型中不再考慮地表供水工程的基礎(chǔ)設(shè)施投資，而主要考慮地下水供水工程的基礎(chǔ)設(shè)施投資和運(yùn)營成本。

供水總量取決于需求，應(yīng)滿足以下關(guān)系式：

Q1+Q2=Qx

(1)

式中，Q1為地表水的有效供水量；Q2為地下水的供水量；Qx為需要達(dá)到的總供水量。設(shè)地表水引水量為Qs，考慮到蒸發(fā)和滲漏損失量，地表水有效供水量與引水量關(guān)系如下：

Q1=(1-η)Qs

(2)

式中，η為損失系數(shù)，按照鴨子蕩水庫的有關(guān)資料，η可取3.3%(源自“寧東能源化工基地一期供水工程水資源論證報告書”)。

對于損失的水量，屬于水資源的浪費(fèi)，引入以下罰函數(shù)：

(3)

式中，Jw為模型中引入的處罰成本；qs為為單方水的水權(quán)置換成本；es為處罰的倍增函數(shù)(即每多浪費(fèi)1 m3水需要多增加的單方成本)。根據(jù)寧夏地區(qū)水權(quán)置換的經(jīng)驗(yàn)，qs=2.0元/m3。es可取0.01元/萬m3，表示每多浪費(fèi)1 萬m3方水單價上漲0.01元。

地下水和地表水供水都需要滿足以下約束條件：

Q1≤QR；Q2≤QP

(4)

式中，QP為地下水的可開采量；QR為地表水庫的供水能力。

對于任一地下水水源地，開采量會引起當(dāng)?shù)毓┧蜕鷳B(tài)環(huán)境壓力，其水權(quán)置換的單方成本隨著實(shí)際開采量與可開采量的比例(環(huán)境壓力增大)增加，計算公式如下：

(5)

式中，qg為地下水的基礎(chǔ)水權(quán)置換單方成本(取2.0元/m3)。

總的經(jīng)濟(jì)成本按下式計算：

(6)

式中，t為模型考慮的投資回收周期(20年)；qsQ1為地表水的水權(quán)置換成本；qe為地下水供水工程的基礎(chǔ)設(shè)施投資單方成本[0.04元/(m3·km)]；Lg為輸水距離；ms為地表水供水工程的運(yùn)營配水成本(1.5元/m3)；mg為地下水供水工程的運(yùn)營配水成本(1.0元/m3)。

規(guī)劃模型的目標(biāo)函數(shù)如下：

Jx=minJx(Qs,Qg)

(7)

即聯(lián)合供水方案在20年內(nèi)的總成本達(dá)到最低。

3.3.2模型求解。對圖3中的供水方案可以分解幾個組成部分，分別求解各自的優(yōu)化模型，然后進(jìn)行匯總。

3.3.2.1陶樂地區(qū)優(yōu)化供水方案。陶樂地區(qū)平羅精細(xì)化工基地的需水量Qx為0.073×108m3/a，其中陶樂水源地的工程可開采量為QP為0.113×108m3/a，目前黃河三棵柳引水工程供水能力QR為0.055×108m3/a，輸水距離Lg為40 km。模型求解得到的結(jié)果如圖4所示。

圖4 陶樂地區(qū)供水方案經(jīng)濟(jì)成本隨地下水供水量的變化曲線Fig.4 Variation curve of economic cost of water supply scheme in Taole area

由圖4可以看出，隨著地下水的供水量增加，供水經(jīng)濟(jì)成本總體呈現(xiàn)下降趨勢。當(dāng)?shù)叵滤墓┧窟_(dá)到0.065×108m3/a時，供水成本降低到最低。陶樂傍河水源地的可開采量達(dá)到0.113×108m3/a，但只要開采0.065×108m3/a就可以滿足最優(yōu)供水目標(biāo)，這樣地表水只需要調(diào)度0.008×108m3/a。

3.3.2.2中部工礦區(qū)優(yōu)化供水方案。中部工礦區(qū)包括橫城工礦區(qū)、靈武工礦區(qū)和鴛鴦湖工礦區(qū)，聯(lián)合供水的總需水量Qx為0.396×108m3/a，地表供水工程為鴨子蕩水庫，除去向煤化工的供水外，剩余供水能力為0.256×108m3/a。地下水的水源地有4個：駱駝井應(yīng)急水源地、現(xiàn)有的靈武崇興水源地、3號備選水源地和2號備選水源地，總計可供水量達(dá)到0.352×108m3/a。首先進(jìn)行總量的優(yōu)化選擇，然后按照各個水源地的供水能力分配供水量，平均輸水距離按照60 km計，模型求解得到的結(jié)果如圖5所示。

圖5 中部工礦區(qū)供水方案經(jīng)濟(jì)成本隨地下水供水量的變化曲線Fig.5 Variation curve of economic cost of water supply scheme in central industrial and mining area

由圖5可以看出，隨著地下水的供水量增加，供水經(jīng)濟(jì)成本總體呈現(xiàn)上升趨勢，但當(dāng)?shù)叵滤墓┧窟_(dá)到0.175×108m3/a時，供水成本降低到最低。因中部礦區(qū)聯(lián)合供水的總需水量Qx為0.396×108m3/a，故最優(yōu)供水方案應(yīng)為地下水的供水總量采用0.175×108m3/a，地表水的供水量為0.245×108m3/a。在這些水源地中，崇興水源地已是靈武市水源地，可不作為寧東基地應(yīng)急供水的水源地。2號水源地與3號水源地按照供水能力和輸水距離進(jìn)行比例劃分，建議開采方案：駱駝井應(yīng)急水源地0.019×108m3/a(占可開采量的41%)；靈武平原2號備選水源地0.089×108m3/a(占可開采量的58%)；靈武平原3號備選水源地0.067×108m3/a(占可開采量的72%)。

3.3.2.3南部工礦區(qū)優(yōu)化供水方案。南部工礦區(qū)包括馬家灘工礦區(qū)和積家井工礦區(qū)，聯(lián)合供水的總需水量Qx為0.054×108m3/a，地表供水工程為劉家溝水庫，供水能力為0.061×108m3/a。地下水水源地為大泉鄉(xiāng)水源地和1號備選水源地，應(yīng)急供水量為0.189×108m3/a，輸水距離計為50 km。模型求解得到的結(jié)果如圖6所示。

圖6 南部工礦區(qū)供水方案經(jīng)濟(jì)成本隨地下水供水量的變化曲線Fig.6 Variation curve of economic cost of water supply scheme in the southern industrial and mining area

由圖6可以看出，采用地下水可以降低供水成本，且全部采用地下水時成本最低。靈武1號備選水源地的可開采量為0.170×108m3/a，超過了南部工礦區(qū)的需水量，因此可以考慮全部由1號備選水源地供水，應(yīng)急供水量0.054×108m3/a，只占可開采量的32%。

3.3.3聯(lián)合供水方案地下水的實(shí)際供水路線和供水配額。為實(shí)現(xiàn)用最低的經(jīng)濟(jì)成本達(dá)到供水目標(biāo)，通過對地表水與地下水聯(lián)合供水方案優(yōu)化，地下水的實(shí)際供水路線和供水配額見表3。

表3 聯(lián)合供水方案地下水的實(shí)際供水路線和供水配額

按總供水量0.294×108m3/a和最佳供水路線聯(lián)合配置地下水，地表水配額只需3.136×108m3/a，既能滿足寧東能源化工基地2020年規(guī)劃對水資源的需求，又能實(shí)現(xiàn)社會經(jīng)濟(jì)和環(huán)境效益的最大化。

4 結(jié)論

(1)該研究提出寧東能源化工基地地表水和地下水聯(lián)合供水的方案，并采用線性規(guī)劃模型對聯(lián)合供水的方案進(jìn)行優(yōu)化，根據(jù)優(yōu)化方案確定了地下水的實(shí)際供水路線和供水配額：陶樂傍河水源地可向平羅精細(xì)化工基地供水0.065×108m3/a；駱駝井應(yīng)急水源地向橫城和鴛鴦湖工礦區(qū)供水 0.019×108m3/a；吳忠平原2號和3號備選水源地向靈武工礦區(qū)供水0.089×108和0.067×108m3/a；大泉地區(qū)的1號水源地向南部的馬家灘和積家井工礦區(qū)供水0.054×108m3/a。

(2)按總供水量0.294×108m3/a和最佳供水路線聯(lián)合配置地下水，地表水配額只需3.136×108m3/a，既能滿足寧東能源化工基地2020年規(guī)劃對水資源需求，又能實(shí)現(xiàn)社會經(jīng)濟(jì)和環(huán)境效益的最大化。