左朝暉，李紹康，楊津津，袁穎，李翔*

1.中國環(huán)境科學(xué)研究院 2.河北地質(zhì)大學(xué)城市地質(zhì)與工程學(xué)院

頁巖氣是蘊藏于頁巖層可供開采的天然氣資源，相比于常規(guī)天然氣，頁巖氣具有開采壽命長和生產(chǎn)周期長的優(yōu)點[1]。中國的頁巖氣資源可采量約為115×1011～361×1011m3，近幾年頁巖氣年產(chǎn)量達到44.6×108m3，頁巖氣的大規(guī)模開發(fā)有助于優(yōu)化我國能源消費結(jié)構(gòu)[2]，同時促進頁巖氣開發(fā)區(qū)域經(jīng)濟發(fā)展。然而，頁巖氣開采過程中需要消耗大量淡水資源[3]，鉆井和壓裂過程所需水量通常達到10 000～24 000 m3，其中有30%～70%耗水被重復(fù)利用，其開采耗水十分嚴(yán)重[4-5]，若不協(xié)調(diào)頁巖氣開發(fā)與水資源之間的關(guān)系，可能導(dǎo)致區(qū)域水資源匱乏[6]。因此，有必要開展頁巖氣水資源承載力評價研究，為頁巖氣開發(fā)區(qū)域水資源可持續(xù)利用政策的制定提供借鑒。

目前，國外對水資源承載力的研究很少，多數(shù)研究將水資源承載力納入可持續(xù)發(fā)展研究范疇[7-9]。國內(nèi)對區(qū)域水資源承載力的研究主要是選取社會、水資源、經(jīng)濟、生態(tài)4個方面的指標(biāo)建立水資源承載力評價體系，運用模糊綜合評價法、系統(tǒng)動力學(xué)法、神經(jīng)網(wǎng)絡(luò)法、生態(tài)足跡法等進行評價[10-11]。在評價指標(biāo)體系方面，針對頁巖氣的耗水特點的區(qū)域水資源承載力研究較少，基于社會、水資源、經(jīng)濟、生態(tài)方面指標(biāo)的區(qū)域水資源承載力評價體系缺少客觀性與全面性，如何通過科學(xué)的模型分析區(qū)域頁巖氣開發(fā)過程中水資源承載力成為重要的研究課題。在評價方法方面，BP (back propagation)神經(jīng)網(wǎng)絡(luò)模型具有廣泛的適應(yīng)能力、學(xué)習(xí)能力和映射能力等特點，非線性數(shù)據(jù)處理能力十分強大，且不需要人工確定指標(biāo)權(quán)重，在多變量非線性系統(tǒng)的建模與預(yù)測方面可取得滿意的結(jié)果[12]。許莉等[13]提出了基于BP神經(jīng)網(wǎng)絡(luò)的水資源承載力綜合評價方法，并描述了方法的應(yīng)用過程；郭曉英等[14]將指標(biāo)作為BP神經(jīng)網(wǎng)絡(luò)的輸出，對廈門市的水資源承載力進行了動態(tài)分析；馮湘華等[15]使用模糊綜合評判模型、TOPSIS模型、徑向基神經(jīng)網(wǎng)絡(luò)模型完成了對淳化縣的水資源承載力綜合評價，得到神經(jīng)網(wǎng)絡(luò)與其他方法評價結(jié)果基本一致的結(jié)論；張彥等[16]在保定市水環(huán)境承載研究中將遺傳算法(genetic algorithm, GA)結(jié)合BP神經(jīng)網(wǎng)絡(luò)模型形成GA-BP神經(jīng)網(wǎng)絡(luò)評價模型，并對比了BP神經(jīng)網(wǎng)絡(luò)模型和GA-BP神經(jīng)網(wǎng)絡(luò)模型的可信度，指出GA-BP神經(jīng)網(wǎng)絡(luò)模型的精度更高且具有很強的泛化能力。

筆者以威遠縣為研究對象，在區(qū)域水資源承載力評價中社會、水資源、經(jīng)濟、生態(tài)指標(biāo)基礎(chǔ)上增加頁巖氣開采特征指標(biāo)，建立頁巖氣開采區(qū)域水資源承載力評價指標(biāo)體系。使用GA優(yōu)化BP神經(jīng)網(wǎng)絡(luò)評價模型，將GA-BP神經(jīng)網(wǎng)絡(luò)組合模型應(yīng)用于頁巖氣開發(fā)區(qū)域的水資源承載力評價，以期為頁巖氣開發(fā)區(qū)域的水資源承載力評價提供思路，同時為促進威遠縣的社會經(jīng)濟活動與水資源協(xié)調(diào)發(fā)展提供參考。

1 研究方法

1.1 BP神經(jīng)網(wǎng)絡(luò)基本原理

BP神經(jīng)網(wǎng)絡(luò)是按照誤差逆向傳播算法訓(xùn)練的多層前饋神經(jīng)網(wǎng)絡(luò)，是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)[17]。BP神經(jīng)網(wǎng)絡(luò)模型一般由輸入層、隱含層和輸出層組成。在運算時將歸一化后的指標(biāo)數(shù)據(jù)作為輸入層，將指標(biāo)期望結(jié)果作為輸出層，輸入層數(shù)據(jù)正向傳播至隱含層，通過傳遞函數(shù)處理后傳播至輸出層，當(dāng)測試輸出與期望結(jié)果之間存在誤差時，將誤差沿原傳播過程反向傳播并逐一調(diào)整各層連接權(quán)值(w1、w2)。經(jīng)過不斷迭代計算，當(dāng)誤差降至設(shè)定范圍時，得到最終的神經(jīng)網(wǎng)絡(luò)模型。BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)如圖1所示。

圖1 BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)示意Fig.1 Schematic diagram of BP neural network structure

1.2 GA基本原理

GA是一種基于自然選擇和遺傳理論的隨機優(yōu)化算法，其通過選擇、交叉和變異等操作對所有個體進行篩選，根據(jù)適應(yīng)度保留優(yōu)秀個體，淘汰適應(yīng)度差的個體，循環(huán)迭代至搜尋到最優(yōu)解[18]。由于其高效率的算法以及強大的全局尋優(yōu)能力，被廣泛應(yīng)用在組合優(yōu)化方面[19]。

1.3 GA-BP神經(jīng)網(wǎng)絡(luò)模型

雖然BP神經(jīng)網(wǎng)絡(luò)具有自學(xué)習(xí)和非線性映射的能力，但其算法收斂速度慢，對初始權(quán)值和閾值非常敏感，容易陷入局部最優(yōu)或出現(xiàn)過擬合現(xiàn)象。遺傳算法從群體出發(fā)，同時將多個個體一起比較，從而提高運算效率，同時具有全局搜索能力而不易陷入局部最優(yōu)。將GA與BP神經(jīng)網(wǎng)絡(luò)結(jié)合，使用GA全局尋優(yōu)的能力來優(yōu)化BP神經(jīng)網(wǎng)絡(luò)模型的初始權(quán)值和閾值，可以提高BP神經(jīng)網(wǎng)絡(luò)的收斂速度，從而使優(yōu)化后的網(wǎng)絡(luò)可更好地進行預(yù)測。具體流程見圖2。

圖2 GA-BP神經(jīng)網(wǎng)絡(luò)流程Fig.2 GA-BP neural network flowchart

2 水資源承載力評價與分析

2.1 研究區(qū)概況

威遠縣位于我國西南地區(qū)，地處四川盆地中南部，行政面積1 289 km2。屬亞熱帶暖濕季風(fēng)氣候，冬春干旱而夏秋多雨。地勢西北高，東南低，主要為低山和丘陵地貌。由于氣候與地形地貌影響，威遠縣降水時空分布不均，水資源開發(fā)利用率高達51%，年平均水資源總量為4.7億m3，但人均水資源占有量僅為600 m3左右，是四川省人均水資源占有量的14，屬于水資源貧乏地區(qū)。

我國已經(jīng)探明的頁巖氣大部分分布于西南地區(qū)，目前已有9個頁巖氣重點建產(chǎn)區(qū)，威遠縣大部分區(qū)域處于長寧—威遠國家級頁巖氣示范區(qū)內(nèi)。長寧—威遠示范區(qū)頁巖氣單井壓裂用水量在20 000 m3以上，其中僅有20%左右的壓裂用水得以循環(huán)使用，導(dǎo)致大部分水資源被注入地下。這使得威遠縣水資源情況更加復(fù)雜，進一步加劇社會、經(jīng)濟發(fā)展和水資源之間的矛盾。

2.2 評價指標(biāo)體系及其分級標(biāo)準(zhǔn)

受經(jīng)濟、社會和生態(tài)的影響與制約，水資源承載力評價是一個復(fù)雜的系統(tǒng)，如何確定評價體系中的指標(biāo)是承載力研究的基礎(chǔ)[20]。因此，基于研究區(qū)頁巖氣開發(fā)與生活、生產(chǎn)用水的矛盾，在區(qū)域水資源承載力研究[14,21]中結(jié)合頁巖氣井開發(fā)規(guī)模和用水特點，確定了由5個方面共16個指標(biāo)構(gòu)成的區(qū)域水資源承載力指標(biāo)體系，詳見表1。

表 1 研究區(qū)水資源承載力評價指標(biāo)體系

國內(nèi)外對于水資源承載力的評價鮮見統(tǒng)一標(biāo)準(zhǔn)，本研究在確定分級標(biāo)準(zhǔn)時依據(jù)以下原則：1)優(yōu)先采用已有國家標(biāo)準(zhǔn)或國際標(biāo)準(zhǔn)；2)盡量參照近似行業(yè)的設(shè)定標(biāo)準(zhǔn)；3)對于沒有參照標(biāo)準(zhǔn)的指標(biāo)，尋求專家經(jīng)驗或者采用等間距法、總分頻率法等進行劃分。依據(jù)上述原則以及GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》、GBT 50331—2002《城市居民生活用水量標(biāo)準(zhǔn)》和《頁巖氣發(fā)展規(guī)劃(2016—2020年)》等國家標(biāo)準(zhǔn)或規(guī)劃文本，將各指標(biāo)對水資源承載力的影響分為5個等級，為了方便區(qū)分研究區(qū)水資源承載力狀態(tài)，分別以“理想狀態(tài)”“可承載”“弱承載”“弱不可承載”和“不可承載”表示，各指標(biāo)對應(yīng)的分級標(biāo)準(zhǔn)見表2。其中水資源承載力指數(shù)為[0.8，1.0]，表明水資源承載力處于適宜人類活動的理想狀態(tài)，經(jīng)濟社會發(fā)展不會受到水資源阻礙，可以適當(dāng)加大發(fā)展力度；[0.6，0.8)，表明區(qū)域水資源承載力處在可以承載人類活動的狀態(tài)，經(jīng)濟社會發(fā)展與水資源的關(guān)系較為協(xié)調(diào)；[0.4，0.6)，表明水資源承載力處在較弱的狀態(tài)，社會經(jīng)濟發(fā)展受到水資源掣肘，需要保護水資源并進行資源分配結(jié)構(gòu)優(yōu)化；[0.2，0.4)，表明水資源承載力處于很弱的狀態(tài)，區(qū)域水資源問題較為突出，經(jīng)濟社會發(fā)展受到水資源阻礙，需要大力保護水資源并提高水資源供給；[0，0.2)，表明水資源承載力處于差的狀態(tài)，區(qū)域內(nèi)社會經(jīng)濟發(fā)展與水資源之間極不協(xié)調(diào)。

表 2 研究區(qū)水資源承載力指標(biāo)分級標(biāo)準(zhǔn)

2.3 GA-BP神經(jīng)網(wǎng)絡(luò)模型構(gòu)建

基于表2的分級標(biāo)準(zhǔn)，采用隨機技術(shù)公式[22]〔式(1)〕在每一等級中生成100組數(shù)據(jù)，其中隨機選出90組作為訓(xùn)練樣本數(shù)據(jù)，剩余10組作為測試樣本數(shù)據(jù)。由此，在5個等級區(qū)間中共計生成500組樣本，其中訓(xùn)練樣本450組，測試樣本50組。

(1)

為避免指標(biāo)體系中各指標(biāo)量綱不統(tǒng)一帶來的影響，采用最大最小方法在形成模型前對數(shù)據(jù)進行歸一化處理。公式如下：

(2)

式中：Xmax和Xmin為數(shù)據(jù)中最大值和最小值；Xij為數(shù)據(jù)的實際值。

選取3層神經(jīng)網(wǎng)絡(luò)模型，經(jīng)過反復(fù)訓(xùn)練設(shè)置隱含層神經(jīng)元數(shù)為6個，訓(xùn)練函數(shù)選用trainlm函數(shù)，模型精度設(shè)置為0.001。采用GA對所建立的網(wǎng)絡(luò)結(jié)構(gòu)進行優(yōu)化，經(jīng)過反復(fù)調(diào)試后確定GA參數(shù)，其中迭代次數(shù)為10，種群規(guī)模為10，交叉概率為0.3，變異概率為0.1，得到最佳適應(yīng)度曲線如圖3所示，最佳適應(yīng)度為1.431 9。將測試樣本數(shù)據(jù)代入訓(xùn)練好的模型，得到測試樣本的擬合值和原始值的線性回歸曲線(圖4)，相關(guān)系數(shù)為0.995 98，測試樣本最大相對誤差為6.50%，表明模型擬合具有良好的效果。

圖3 適應(yīng)度曲線Fig.3 Fitness curve

圖4 測試樣本回歸曲線Fig.4 Test sample regression curve

將相同網(wǎng)絡(luò)結(jié)構(gòu)的BP神經(jīng)網(wǎng)絡(luò)模型和GA-BP神經(jīng)網(wǎng)絡(luò)模型的測試樣本數(shù)據(jù)誤差進行對比，結(jié)果見圖5。由圖5可知，單一BP神經(jīng)網(wǎng)絡(luò)模型的評價結(jié)果波動明顯，在第40～50個樣本處理過程中出現(xiàn)較大誤差；GA-BP神經(jīng)網(wǎng)絡(luò)模型的評價結(jié)果未出現(xiàn)劇烈的波動，相比單一BP神經(jīng)網(wǎng)絡(luò)模型誤差更小，表明該模型在頁巖氣開發(fā)區(qū)域的水資源承載力評價中具有一定的適用性。

圖5 模型驗證誤差比較Fig.5 Comparison of validation errors between two models

2.4 評價結(jié)果與分析

評價指標(biāo)數(shù)據(jù)來源于歷年《內(nèi)江市水資源公報》《內(nèi)江市統(tǒng)計年鑒》和《威遠縣國民經(jīng)濟和社會發(fā)展統(tǒng)計公報》等，部分缺少數(shù)據(jù)采用插值法補全，具體數(shù)據(jù)見表3。選用訓(xùn)練好的BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，將研究區(qū)2014—2019年的評價指標(biāo)數(shù)據(jù)作為輸入層，將西南地區(qū)頁巖氣開發(fā)狀態(tài)下水資源承載力作為輸出層，對研究區(qū)2014—2019年水資源承載力進行評價，得到威遠縣頁巖氣開發(fā)區(qū)域的水資源承載力指數(shù)及承載力狀態(tài)，結(jié)果見表4。由表4可以看出，2014—2017年水資源承載力均處于可承載狀態(tài)，說明在現(xiàn)有頁巖氣開發(fā)狀態(tài)下，水資源相對充足，頁巖氣開發(fā)處于比較適合的狀態(tài)；2018年和2019年水資源承載力處于弱承載狀態(tài)，說明社會經(jīng)濟活動和水資源處于非協(xié)調(diào)狀態(tài)。從趨勢上來看，2014—2019年的水資源承載力雖然沒有達到不可承載的水平，但水資源承載力指數(shù)除2016年稍有回升外，均處于逐年下降趨勢。

進一步分析發(fā)現(xiàn)，2015年人均生活用水量、萬元工業(yè)增加值用水量、頁巖氣井群規(guī)模等均較2014年有所增加，而人均水資源量達到近幾年的最低水平，其他指標(biāo)變化不大，因此可以推斷水資源供給不足、生活用水劇烈增加和工業(yè)用水(包括頁巖氣開采用水)量增加3個方面的原因?qū)е滤Y源承載力劇烈波動；2016年人均水資源量達到近幾年最高水平，水資源供給充足，井群規(guī)模僅增加23口，人均生活用水量相比2015年減少45 m3，萬元工業(yè)增加值用水量減少，且水質(zhì)達標(biāo)率、植被覆蓋率、年污水處理率等指標(biāo)均有所增長，在這些因素的共同影響下，水資源承載力指數(shù)較2015年有所增長。由2014—2016年水資源承載力變化可以看出，其主要受到人均水資源量、人均生活用水量、萬元工業(yè)增加值用水量和井群規(guī)模等指標(biāo)的影響。根據(jù)頁巖氣井群規(guī)模、單井需水量、壓裂返排液回用率和壓裂液返排率4項指標(biāo)數(shù)據(jù)推算研究區(qū)頁巖氣總需水量、回用水量和新鮮水使用量(表5)，發(fā)現(xiàn)2017年之后雖然人均水資源量、人均生活用水量、萬元工業(yè)增加值用水量等主要指標(biāo)并未惡化，但頁巖氣井群規(guī)模急劇擴大，導(dǎo)致開發(fā)總需水量急劇增長，研究區(qū)水資源承載力指數(shù)下降明顯，這表明頁巖氣系統(tǒng)方面的耗水壓力變大是水資源承載力逐年下降的主要原因。

表3 2014—2019年評價指標(biāo)數(shù)據(jù)

表4 水資源承載力評價結(jié)果

綜上所述，研究區(qū)水資源承載力受到社會、水資源、生態(tài)等多方面的影響，其中主要影響因素為人均水資源量、萬元工業(yè)增加值用水量和井群規(guī)模。由于頁巖氣開發(fā)活動不斷增強，水資源承載力總體處

表5 頁巖氣開發(fā)所需水量數(shù)據(jù)

于下降趨勢，尤其在2017—2019年開發(fā)規(guī)?？焖贁U大時，水資源承載力明顯下降，但由于社會、水資源和生態(tài)三方面的指標(biāo)均比較穩(wěn)定，萬元工業(yè)增加值用水量、年污水處理率、單井需水量、壓裂返排液回用率等工業(yè)指標(biāo)負荷逐漸降低，區(qū)域水資源承載力并未出現(xiàn)不可承載的情況。因此，研究區(qū)在頁巖氣開發(fā)時要充分考慮頁巖氣耗水的影響，在井群規(guī)模維持現(xiàn)狀或繼續(xù)擴大的情況下，應(yīng)進一步加強水資源的節(jié)約和保護，同時調(diào)整產(chǎn)業(yè)結(jié)構(gòu)、建設(shè)節(jié)水型社會，通過節(jié)水、回用和暫時調(diào)整經(jīng)濟發(fā)展重心等手段不斷提高經(jīng)濟社會發(fā)展過程中的用水效率，保障經(jīng)濟社會的可持續(xù)發(fā)展。

3 結(jié)論

(1)在常規(guī)區(qū)域水資源承載力影響因素中充分考慮頁巖氣開發(fā)對區(qū)域水資源承載力的影響，增添頁巖氣開采耗水指標(biāo)，建立適合頁巖氣開發(fā)區(qū)域的水資源承載力評價體系。采用遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)組合模型中的權(quán)值和閾值進行評價，與單一的BP神經(jīng)網(wǎng)絡(luò)模型相比，組合模型收斂速度更快，誤差較小，評價結(jié)果更為穩(wěn)定。

(2)研究區(qū)水資源承載力總體處于波動下降狀態(tài)，主要受到人均水資源量、頁巖氣的井群規(guī)模和萬元工業(yè)增加值用水量等指標(biāo)影響。2014—2017年水資源承載力處于可承載水平，水資源與經(jīng)濟、社會、頁巖氣發(fā)展處于相對平衡狀態(tài)，區(qū)域水資源供給仍處于較協(xié)調(diào)的狀態(tài)；2018年和2019年研究區(qū)水資源承載力處于弱承載水平，頁巖氣開采規(guī)模急劇擴大改變了研究區(qū)的水資源承載力結(jié)構(gòu)，加劇了社會經(jīng)濟發(fā)展與水資源之間的不平衡。

(3)研究區(qū)萬元工業(yè)增加值用水量、年污水處理率、單井需水量、壓裂返排液回用率等指標(biāo)的持續(xù)好轉(zhuǎn)可以緩解頁巖氣大規(guī)模開發(fā)帶來的水資源壓力，因此，威遠縣可以在節(jié)約和保護水資源的基礎(chǔ)上，進一步降低頁巖氣井單井用水量，提高區(qū)域內(nèi)經(jīng)濟社會發(fā)展過程中的用水效率，通過重新優(yōu)化水資源分配結(jié)構(gòu)，推動經(jīng)濟社會發(fā)展與水資源承載能力相適應(yīng)。