基于INVEST模型的黑水河流域水質(zhì)風(fēng)險評估

徐凌云,丁 玲

(上?？睖y設(shè)計研究院有限公司,上海 200434)

隨著人口擴(kuò)張和城市發(fā)展伴隨而來的生態(tài)系統(tǒng)的破壞，使得人們對生態(tài)系統(tǒng)服務(wù)重要性的認(rèn)識與日俱增，與此同時生態(tài)系統(tǒng)服務(wù)功能價值的評估已成為了生態(tài)學(xué)和生態(tài)經(jīng)濟(jì)學(xué)研究的熱點，生態(tài)經(jīng)濟(jì)學(xué)家根據(jù)生態(tài)系統(tǒng)提供的各項價值并結(jié)合環(huán)境經(jīng)濟(jì)學(xué)，逐步確立了一些評估生態(tài)系統(tǒng)服務(wù)功能的模型。目前常用的生態(tài)系統(tǒng)服務(wù)功能價值評估模型有INVEST模型、生態(tài)足跡模型、CITYgreen模型、EVR模型(Ecological Value at Risk)，其中，生態(tài)足跡模型[1-3]能客觀反映人類對生態(tài)系統(tǒng)(生態(tài)基礎(chǔ)設(shè)施)的需求與供給之間的矛盾，指示自然資源的壓力狀態(tài)，但缺點在于產(chǎn)量因子存在偏差，計算結(jié)果有高估地區(qū)生態(tài)狀態(tài)的可能，而且只涉及自然資源，對于人類可持續(xù)發(fā)展的其他方面難以測算；CITYgreen模型[4-5]可用于城市森林(小區(qū)域如一個公園、一個社區(qū)范圍，大區(qū)域如整個縣、整個市區(qū))的結(jié)構(gòu)分析與生態(tài)效益評價，同時將結(jié)果以報告形式輸出，其基于GIS軟件ArcView開發(fā)，對遙感影像的質(zhì)量要求較高。EVR模型(Ecological Value at Risk)[6-7]將生態(tài)服務(wù)價值的定量化與生態(tài)風(fēng)險分析的數(shù)學(xué)模型相結(jié)合，可以進(jìn)行基于生態(tài)服務(wù)價值的生態(tài)風(fēng)險分析研究，但模型對數(shù)據(jù)要求量較大，且模型是借鑒金融分析市場組合的VaR方法，設(shè)定了假設(shè)條件，可推廣性差。而InVEST模型[8-9]能較好地把握總體格局，體現(xiàn)人類活動對生境的威脅程度和影響范圍，使利益相關(guān)者權(quán)衡發(fā)展與保護(hù)的關(guān)系，尋求最佳策略。

InVEST(The Integrate Valuation of Ecosystem Services and Tradeoffs Tool)[10-11]是由美國斯坦福大學(xué)、世界自然基金(WWF)和大自然保護(hù)協(xié)會(TNC)于2007年聯(lián)合開發(fā)出的一種生態(tài)系統(tǒng)服務(wù)和交易的綜合評估模型。該模型基于GIS平臺，可以定量化評估各項生態(tài)服務(wù)功能，并以地圖的形式體現(xiàn)其空間分布特征，為自然資源管理提供最優(yōu)的決策支持。目前許多專家、學(xué)者把InVEST作為一款實用的分析工具，將其應(yīng)用于生態(tài)環(huán)境、水利水電、土地資源、水產(chǎn)養(yǎng)殖等不同的領(lǐng)域[12-13]。本研究以黑水河流域為例，介紹了如何利用INVEST模型評價黑水河流域營養(yǎng)鹽的空間分布特征，從不同的角度分析評價黑水河流域水質(zhì)威脅，并提出保護(hù)建議，為該流域生態(tài)和自然資源的保護(hù)提供決策依據(jù)。

1 研究區(qū)域與研究方法

1.1 流域概況

黑水河(圖1)是金沙江左岸一級支流，位于四川省涼山彝族自治州境內(nèi)，發(fā)源于昭覺縣瑪果梁子，自北向南流經(jīng)昭覺、普格、寧南3縣，于寧南縣東南部葫蘆口注入金沙江；流域面積為3 591 km2，全長173.0 km，天然落差1 931 m，平均比降1.105%，河口處多年平均流量為80 m3/s、徑流量25.25億m3。通過黑水河干流的棲息地修復(fù)與建設(shè)可發(fā)揮其流水生境對長江上游特有魚類的保護(hù)作用，并為白鶴灘庫區(qū)喜緩流和靜水生境但需流水刺激產(chǎn)卵的魚類提供適宜的水生生境。

黑水河流域淡水資源供應(yīng)是一項有助于社會發(fā)展的重要生態(tài)系統(tǒng)服務(wù)功能，保障了流域農(nóng)業(yè)灌溉、人口增加、工業(yè)發(fā)展和旅游活動等。目前，黑水河流域的工業(yè)開發(fā)強度較低，流域內(nèi)威脅水環(huán)境的主要因素有農(nóng)業(yè)面源、場鎮(zhèn)居民生活污水、采石挖沙及沿河工礦企業(yè)生產(chǎn)廢水、養(yǎng)殖廢水，以及地質(zhì)災(zāi)害引發(fā)的水環(huán)境威脅。為進(jìn)一步分析全流域氮、磷營養(yǎng)物空間分布情況，利用InVEST2.5.4中的產(chǎn)水量模塊和營養(yǎng)鹽截流(營養(yǎng)物沉積模塊)模塊模擬流域氮、磷營養(yǎng)物負(fù)荷分布，并在此基礎(chǔ)上對黑水河流域水質(zhì)威脅進(jìn)行分析評估，為黑水河干流的棲息地修復(fù)與建設(shè)提供決策依據(jù)。

1.2 研究方法

利用產(chǎn)水量模塊模擬黑水河流域的降雨徑流，然后將產(chǎn)水量空間分布導(dǎo)入營養(yǎng)鹽截流模塊，營養(yǎng)鹽截流模塊通過提取數(shù)字高程模型(DEM)的坡度，確定徑流路徑，通過分析匯流能力閾值提取出河網(wǎng)信息，并考慮了不同植被的截留能力，模擬了氮磷自一網(wǎng)格遷移至下一網(wǎng)格過程中的沉積與截留的損失量，以此來模擬氮磷營養(yǎng)物隨徑流輸入至河道的空間分布。

1.2.1 產(chǎn)水模塊

InVEST模型的產(chǎn)水模塊是一種基于Budyko曲線和年平均降水的水量平衡估算方法[14]，不考慮地表水與地下水的交互作用，認(rèn)為某柵格單元的降雨量減去沒有上游徑流補給時蒸散發(fā)后的水量，即為產(chǎn)水量。該模型首先要計算出柵格單元某種特定植被類型下的年平均產(chǎn)水量，由公式(1)計算獲得：

(1)

式中Yxj——柵格單元x中植被類型為j的年產(chǎn)水量；AETxj——柵格單元x中植被類型j的實際蒸散量；Px——柵格單元x的年平均產(chǎn)水量?？捎葿udyko曲線公式(2)近似得到：

(2)

式中Rxj——Bydyko的干燥指數(shù)；ωx——自然氣候-土壤性質(zhì)的非物理參數(shù)，分別由公式(3)和公式(4)計算得到。

(3)

式中Kxj——柵格單元x中土地利用覆被類型下的植被蒸散發(fā)系數(shù)；ETox——柵格單元x中的參考蒸散。

(4)

式中AWCx——土壤有效含水量，是由土壤質(zhì)地和有效土層厚度決定的；Z——季節(jié)性因子，代表季節(jié)性降雨分布和降水深度，不同的研究區(qū)域具有不同的適宜值。

1.2.2 營養(yǎng)鹽截留模塊

地表徑流中的氮磷通過植被緩沖帶后，其中部分氮磷通過植物吸收、微生物固定、硝化反硝化作用后被截留，部分氮磷通過土壤吸附而滲透到土層中、沉積礦化而被截留。營養(yǎng)物截留模塊基于改進(jìn)后的輸出系數(shù)法，考慮了氮磷營養(yǎng)物在輸移過程中由于截留和沉積作用而不可能完全進(jìn)入收納水體，以此來模擬氮磷營養(yǎng)物流入河川徑流中的空間分布：

ALVx=HSSx·polx

(5)

式中ALVx——柵格x調(diào)節(jié)的載荷值；polx——柵格x的輸出系數(shù)；HSSx——柵格x的水文敏感性得分值。

1.3 數(shù)據(jù)來源及處理

a) DEM數(shù)據(jù)。采用流域?qū)崪y30 m×30 m的DEM數(shù)據(jù)。DEM圖中可以看出，黑水河流域北高南低，地貌多以山嶺和河谷相間，河谷到分水嶺高差較大，在700～4 000 m左右，沿河支溝發(fā)育。

b) 土地利用。土地利用數(shù)據(jù)采用2016年度實測的流域1∶250 000土地利用圖。從土地利用圖可以看出，目前黑水河流域土地利用方式主要有坡耕地、耕地、梯田、林地、建設(shè)用地、裸地、水域及水利設(shè)施用地等，以林地和坡耕地占地面積最大，分別占整個流域面積的42.0%和47.9%。坡耕地和耕地主要沿黑水河、則木河、西洛河及其兩側(cè)支溝的河谷地帶分布。根據(jù)對黑水河流域的初步調(diào)查分析，威脅流域內(nèi)水文、水質(zhì)及自然生態(tài)系統(tǒng)的主要因素還包括居民點、采石挖沙、沿河企業(yè)、養(yǎng)殖等幾種方式，將這幾種類型更新至土地利用圖中進(jìn)行校正。

c) 土壤數(shù)據(jù)。土壤數(shù)據(jù)采用收集到的1∶100萬數(shù)字化土壤圖，包括土壤厚度、土壤顆粒組成、持水率等基本信息，土壤類型采用收集到的1∶250 000土壤類型空間分布圖。

d) 降雨與植物潛在蒸散量。降雨數(shù)據(jù)采用黑水河流域的多年平均年降雨量。根據(jù)InVEST用戶手冊上的建議和文獻(xiàn)查閱，發(fā)現(xiàn)修正的哈格里夫斯公式在計算年均潛在蒸散量上應(yīng)用較廣，采用該公式計算研究區(qū)域ET0。修正的哈格里夫斯公式如下：

ET0=0.0013×0.408×RA×(Tav+17)×

(TD-0.0123P)0.76

式中RA——地球外輻射(太陽總輻射)，MJ/(m2·d-1)；Tav——年均日最高氣溫和年均日最低氣溫的平均值，℃；TD——年均日最高氣溫和年均日最低氣溫的差值，單位是℃；P——月均降雨量，mm。

通過中國氣象科學(xué)數(shù)據(jù)共享平臺，得到寧南站多年日最高氣溫和年均日最低氣溫，計算出Tav和TD。查閱四川省太陽能輻射分布圖，得到研究區(qū)域的RA值，取地區(qū)均值5 760 MJ/(m2·a-1)。根據(jù)哈格里夫斯公式最終計算得出黑水河流域植物潛在蒸散量的空間范圍為605～635 mm/a。

e) 模型參數(shù)取值。InVEST模型參數(shù)是基于不同土地利用類型的參數(shù)集合，主要包括植被最大根深、蒸散發(fā)系數(shù)及各個營養(yǎng)鹽的輸出系數(shù)、截留能力。其中植被最大根深數(shù)據(jù)參考Canadell等對于全球范圍內(nèi)植被最大根深的研究[15]，并進(jìn)行了率定；蒸散系數(shù)參考InVEST用戶指南，并進(jìn)行了率定；而不同土地利用類型的氮磷輸出系數(shù)和氮磷去除效率參考了國內(nèi)外文獻(xiàn)和三峽庫區(qū)農(nóng)業(yè)面源污染負(fù)荷[16]，并進(jìn)行了率定。率定結(jié)果具體見表1。

表1 模型參數(shù)取值

2 結(jié)果與分析

2.1 產(chǎn)水量模塊的校驗與營養(yǎng)鹽輸出空間分布

應(yīng)用INVEST模型分別對流域產(chǎn)水和耗水進(jìn)行模擬，將產(chǎn)水量和耗水量的差作為流域?qū)嶋H出水量。采用寧南水文站多年平均流量對模擬結(jié)果進(jìn)行驗證(表2)。根據(jù)寧南水文站代表年多年年均流量為69.1 m3/s，年徑流量為21.79億m3，經(jīng)過對產(chǎn)水量模塊的反復(fù)調(diào)試，模擬出的寧南水文站多年平均徑流量為20.76億m3，誤差為4.7%，與實際情況吻合較好。

將校驗后的產(chǎn)水量(mm)柵格圖導(dǎo)入營養(yǎng)鹽截留模塊，可以提取黑水河流域出口處的氮、磷營養(yǎng)鹽濃度數(shù)據(jù)，將黑水河流域出口處的總氮和總磷的濃度值與2016年實測豐水期、平水期、枯水期及年均值進(jìn)行對比，對比發(fā)現(xiàn)模擬計算出的黑水河河口處的總氮和總磷年均濃度相對誤差分別為3%和20%(表3)。

表2 寧南水文站多年平均實測徑流量與計算徑流量對比

表3 黑水河河口處營養(yǎng)鹽濃度計算值與實測值對比

2.2 黑水河流域營養(yǎng)鹽負(fù)荷水平

根據(jù)INVEST模型模擬出2016年全年黑水河流域氮、磷營養(yǎng)鹽負(fù)荷空間分布及各個子流域的氮、磷營養(yǎng)鹽負(fù)荷空間分布，見圖2。

從圖2、3中可以看出，營養(yǎng)鹽的負(fù)荷輸出差異十分明顯，其中山嶺地區(qū)營養(yǎng)物輸出量小，而河谷地區(qū)輸出量大，這和區(qū)域的土地利用類型分布基本一致。山嶺地區(qū)土地利用類型主要為林地，人工開發(fā)強度低，面源污染產(chǎn)生的氮磷負(fù)荷量小，加上林地植被繁茂，對氮、磷污染物的截留率高，這些因素大大降低了山嶺地區(qū)的氮磷輸出。

各子流域氮、磷輸出負(fù)荷差異也較大，黑水河和則木河沿線所在子流域輸出負(fù)荷量較大，流域北部及黑水河?xùn)|西兩側(cè)的山嶺地區(qū)所在的子流域輸出負(fù)荷量較小，各子流域中氮、磷單位面積年輸出負(fù)荷最高值分別為26.17 kg/ha、2.62 kg/ha。

2.3 黑水河流域營養(yǎng)鹽濃度水平

為更直觀的反映黑水河流域氮、磷營養(yǎng)物質(zhì)的現(xiàn)狀水平及水環(huán)境類別，通過對各子流域氮磷輸出總量和實際產(chǎn)水量分布圖疊加完成單位換算，得出各子流域氮、磷平均濃度空間分布，見圖3。

從圖3可以看出，各子流域氮、磷平均濃度變化較大，總氮濃度在0.12～3.4 mg/L左右，總磷濃度在0.01～0.34 mg/L左右。根據(jù)GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》對黑水河流域氮磷濃度水平進(jìn)行評價，總氮和總磷濃度滿足各類水質(zhì)標(biāo)準(zhǔn)的子流域面積占流域總面積的比例見表4?？傮w來說，大部分子流域水環(huán)境質(zhì)量較好，少數(shù)子流域水體中氮、磷含量超過三類水標(biāo)準(zhǔn)，而超標(biāo)區(qū)域主要集中在黑水河干流、西洛河干流和則木河、西洛河交匯處所在的小流域。

表4 營養(yǎng)鹽濃度滿足相應(yīng)水質(zhì)標(biāo)準(zhǔn)的子流域面積百分比

3 結(jié)果討論

對黑水河流域水質(zhì)主要威脅源進(jìn)行分析和討論。

a) 從不同土地利用類型的營養(yǎng)鹽平均輸出量角度分析。根據(jù)氮磷營養(yǎng)物輸出結(jié)果，氮磷營養(yǎng)鹽輸出量與土地利用類型密切關(guān)聯(lián)。將黑水河流域氮磷營養(yǎng)鹽空間分布圖與土地利用圖進(jìn)行疊加統(tǒng)計，得出不同土地利用類型的氮、磷平均輸出量，不同土地利用類型的氮、磷平均輸出量及貢獻(xiàn)率見表5。

根據(jù)表5的統(tǒng)計結(jié)果，養(yǎng)殖和采石采砂的氮、磷平均輸出量較高，水域的其次。養(yǎng)殖以及采石采砂的高氮磷輸出量與其自身產(chǎn)生的污染負(fù)荷較高有關(guān)；水體由于基本位于流域出水口，接受上游輸入的污染負(fù)荷高，其自身截留能力差，因而輸出的氮磷負(fù)荷較高；耕地、坡耕地和梯田的氮、磷輸出量相近，要略低，這與其本身的截留能力較強有關(guān)；林地的截留能力相比耕地更強，而產(chǎn)生的負(fù)荷也更小，所以氮、磷輸出量最??；建設(shè)用地、企業(yè)、居民的氮、磷輸出量與其負(fù)荷產(chǎn)生量和截留能力密切相關(guān)；裸地截留能力最差，但是由于其污染物負(fù)荷產(chǎn)生量也小，因此氮、磷輸出量也不高。

表5 不同土地利用類型氮、磷平均輸出量及貢獻(xiàn)率

b) 從超過三類水標(biāo)準(zhǔn)的子流域營養(yǎng)鹽平均輸出量角度分析。對氮、磷濃度超過三類水標(biāo)準(zhǔn)的子流域各土地利用類型對應(yīng)的單位污染負(fù)荷貢獻(xiàn)率進(jìn)行統(tǒng)計，總氮濃度超三類的子流域有12個，總氮濃度超三類的子流域氮的平均輸出量的主要威脅因子貢獻(xiàn)率見表6?？偭诐舛瘸惖淖恿饔蛴?個，總磷濃度超三類的子流域磷的平均輸出量的威脅因子貢獻(xiàn)率見表7。從兩張表中可以看出，采石采砂、養(yǎng)殖、建設(shè)用地、裸地、坡耕地等的單位面積氮磷的高輸出是造成子流域氮磷濃度超標(biāo)的主要威脅源，如子流域146、251單位面積氮磷輸出量最大的是建設(shè)用地，子流域84、175單位面積氮磷輸出量最大的是采石采砂，子流域183、207單位面積氮磷輸出量最大的是養(yǎng)殖。

表6 總氮濃度超三類的子流域單位面積氮輸出量的主要威脅因子及其貢獻(xiàn)率

表7 總磷濃度超三類的子流域單位面積磷輸出量的主要威脅因子及其貢獻(xiàn)率

4 結(jié)論

綜合上述分析，黑水河流域氮磷營養(yǎng)負(fù)荷的高低與流域土地利用方式和開發(fā)強度密切相關(guān)，對氮磷濃度超過三類水標(biāo)準(zhǔn)的子流域分析發(fā)現(xiàn)，城鎮(zhèn)建設(shè)(人口)、養(yǎng)殖、農(nóng)業(yè)面源、采石采砂等是影響這些子流域水質(zhì)狀況的主要威脅源，應(yīng)該積極采取相應(yīng)保護(hù)措施，改善這些子流域的水質(zhì)狀況，如加強城鎮(zhèn)基礎(chǔ)設(shè)施建設(shè)、應(yīng)用農(nóng)業(yè)最佳管理措施，禁止干流河道采砂等，同時由于坡耕地的氮磷輸出總量的貢獻(xiàn)率較大，也需要對氮磷濃度超過三類水標(biāo)準(zhǔn)的子流域中的坡耕地采取相應(yīng)的措施，如將坡耕地改造為梯田。

總體而言，InVEST模型較好地模擬了黑水河流域產(chǎn)水量和氮磷營養(yǎng)鹽負(fù)荷，并且給出了不同土地利用類型的氮磷營養(yǎng)鹽負(fù)荷的主要威脅因子及其貢獻(xiàn)率，以及根據(jù)水質(zhì)標(biāo)準(zhǔn)給出了威脅程度及貢獻(xiàn)率較高的區(qū)域(子流域)。