基于GIS和AHP的南通市主城區地下水地源熱泵適宜性評價

周紅衛， 王　琦， 張沙莎， 楊　儀， 郭　華

(江蘇省地質環境勘查院，江蘇南京210002)

?

基于GIS和AHP的南通市主城區地下水地源熱泵適宜性評價

周紅衛， 王琦， 張沙莎， 楊儀， 郭華

(江蘇省地質環境勘查院，江蘇南京210002)

摘要：將南通市主城區的地質和水文地質條件、地下水動力場、溫度場與熱物性、地下水化學場和環境地質等自然因素作為層次分析法(AHP)模型的評價因子，應用AHP獲得各影響因子的權重值，并結合GIS的空間分析技術，將研究區地下水地源熱泵適宜性分為適宜、較適宜和不適宜3個大區，其中適宜區和較適宜區占總調查區的97.78%。經檢驗，評價方法適合該研究區。

關鍵詞：層次分析法(AHP)； 地下水地源熱泵； MapGIS空間分析； 適宜性分區； 江蘇南通

0引言

淺層地溫能開發利用方式主要為地埋管式和地下水式地源熱泵。對于地埋管地源熱泵方式，適宜性分區主要考慮巖土體特性、地下水分布和滲流情況等因素；對于地下水地源熱泵方式，則主要考慮含水層巖性、分布、埋深、厚度、富水性、滲透性、地下水動態變化、地質災害等因素(王貴玲等，2012)。

結合GIS的空間分析技術，按照層次分析法計算所得的分區標準，可將多個評價因子間交互計算后的屬性結果在GIS圖件上直觀表達。層次分析法(Analytic Hierarchy Process，AHP)(Saaty，1986)是一種簡便、靈活而又實用的多準則決策方法，適合于針對淺層地溫能開發利用適宜性評價(董殿偉等，2010;苑雷等，2013；楊澤等，2013;Chandio et al.，2013)。

南通市主城區地下水地源熱泵適宜性的分區工作，可為淺層地溫能資源量評價、開發利用規劃及政府管理提供可靠依據。

1分區模型的構建與評價

1.1研究區概況

南通市位于江蘇省東南隅、長江三角洲前緣，東臨黃海，南靠長江，瀕江臨海。研究區為南通市主城區，包括崇川區、港閘區、經濟技術開發區3個區，面積約470.0 km2，其中陸域面積約410.0 km2，水域面積約60.0 km2(圖1)。

圖1　研究區交通位置圖Fig.1　Map showing traffic location of the study area

中國廣大平原松散地層分布區，或具備地下換熱條件的基巖出露區，都具有開發利用淺層地溫能的條件(韓再生等，2007)。南通市位于長江三角洲平原，成陸時間不長，區內地形平坦，河渠縱橫，水網密布，長江沿岸屹立有5個規模較小的孤丘，高程約30～110 m。該區地下水水源充足，富水性分區明顯，為淺層地溫能的開發利用提供了天然地質條件。本次淺層地溫能調查研究的含水層自上而下為潛水層、第I承壓含水層、第II承壓含水層、第Ⅲ承壓含水層。區內各含水層均通過黏性土弱透水含水層發生水力聯系，其中第I承壓水含水層(組)由上更新統沖積、沖海積松散地層組成，分布廣泛，含水層頂板埋深40～60 m，富水性較好，單井涌水量1 000～5 000 m3/d。水質較復雜，受海水及后期地表水的補給淡化所影響，具有明顯的分帶性，其一般規律為由西向東，礦化度增加，水質變差。第Ⅲ承壓水含水層為南通市地下水主采層，經過數十年的大量或過量開采，主采層水位持續下降，改變了含水層組內部和含(隔)水層組間的原始相平衡狀況，水質呈咸化趨勢(單衛華等，2008)。

結合南通市主城區特有的地質、水文地質條件和地源熱泵資源潛力評價要求，優先選擇富水性較好的第Ⅰ承壓含水層作為評價含水層，第Ⅰ承壓含水層各要素分區圖和等值線圖作為參與適宜性分區計算的各評價因子的重要屬性。

1.2 AHP模型建立——評價因子的選擇

為了突出影響地源熱泵適宜性的主導因子，參考了國內外較為普遍適用的評價指標，主要有地質、水文地質條件，地下水動力場，地質體的溫度場、熱物性，地下水化學特征以及施工條件等(董殿偉等，2010；苑雷等，2013；楊澤等，2013)。

在前人研究的基礎上，通過物探測井，抽水回灌試驗等勘察手段，確定研究區地源熱泵適宜性分區主要受5個條件控制：地質與水文地質條件，地下水動力場，溫度場與熱物性，地下水化學場和環境地質。評價因子分列如下。

1.2.1地質、水文地質條件(1) 富水性：全區根據含水層特征及富水性劃分為4個富水區：水量豐富區，單井涌水量>3 000 m3/d；水量中等區，單井涌水量為1 000～3 000 m3/d；水量較小區，單井涌水量<1 000 m3/d；缺失區，無第Ⅰ承壓含水組區及基巖區。

(2) 咸水底界埋深。全區根據咸水底界埋深等值線與第I承壓含水層隔水底板的相對位置可分為2個區：咸水底界埋深>140 m，位于隔水底板下覆；咸水底界埋深120～140 m，局部隔水底板埋深小于咸水底界埋深。

(3) 含水層回灌能力。依據回灌實驗將研究區含水層回灌能力分為2個區：回灌較好區，采灌比0.35～0.65；回灌一般區，采灌比>0.65。

1.2.2地下水動力場(1) 地下水位埋深。根據水位統測，研究區地下水位埋深可劃分為3個區，即地下水位埋深<2、2～4、>4 m。

(2) 含水層滲透系數。根據抽水試驗，滲透系數分為3個區：中部滲透性一般區，滲透系數10～30 m/d；其他2個滲透性較好的區，滲透系數>30 m/d。1.2.3溫度場、熱物性(1) 地層溫度梯度。根據溫度測溫，地層溫度梯度可分為3個區：中部溫度梯度較大區，溫度梯度>3 ℃每100 m；其他2個地層溫度梯度較小區，溫度梯度<3 ℃每100 m。

(2) 熱傳導率。研究區不同地點，同一地層熱傳導率可分為2個區：熱傳導率較低區1.0～1.3 W/(m·K)，熱傳導率較高區1.3～1.9 W/(m·K)。

(3) 比熱容。根據研究區不同地點，同一地層比熱容可分為2個區：比熱容較低區1.0～1.2 kJ/(kg·K)，比熱容較高區1.2～1.52 kJ/(kg·K)。

1.2.4地下水化學場水質：第Ⅰ承壓含水組由于受到二次海侵的影響，礦化度由西向東逐漸增高，劃分為3個區：礦化度<1 g/L的淡水區，礦化度為1～3 g/L的微咸水區，礦化度>3 g/L的咸水區。

1.2.5環境地質根據多年的沉降監測數據，南通市主城區地面沉降率可劃分為3個區，即年沉降量>4、3～4、<3 mm。

地質、水文地質條件決定著淺層地溫能的賦存環境，含水層出水能力、回灌能力及抽回灌井比例等決定場地是否能使用地下水地源熱泵方式。地下水動力場需結合含水層厚度，地下水流動狀況，地下水超采狀況等確定地下水地源熱泵工程適宜性。水化學場分析主要針對城市地下水源地分布區及地下水資源的保護區劃，結合不同層位地下水水質資料，提出地下水地源熱泵工程適宜區域，確定防止地下水含水層污染的適宜開采回灌層位。

1.3影響因子的權重

1.3.1建立層次結構模型對南通市主城區水文地質勘察數據進行分析，根據適宜性分析目的性及評價因子的從屬關系，建立層次分析模型(圖2)。

1.3.2構造判斷矩陣在建立層次結構關系后，上下層次之間元素的隸屬關系被確定，根據1～9標度法，確定同一層次下元素間兩兩間相對重要程度，即權重，最終建立判斷矩陣(表1—表4)。

圖2　適宜性分區層次結構模型Fig.2　Hierarchical structure model of suitability zoning

AB1B2B3B4B5B113462B21/314/322/3B31/43/4121/2B41/61/21/211/3B51/23/2231

表2　B1-C層次判斷矩陣

表3　B2-C層次判斷矩陣

表4　B3-C層次判斷矩陣

(1)

(2)

(3)

A-B層判斷矩陣λmax=5.01，特征向量為(0.828 6，0.276 2，0.220 2，0.130 8，0.414 3)，CR=0.002 2<0.1，滿足一致性檢驗，歸一化得權重Wi=(0.443 1，0.147 7，0.117 7，0.069 9，0.221 5)。B-C層所有判斷矩陣CR均為0，滿足一致性檢驗，其中B1-C特征向量歸一化權重Wi=(0.562 5，0.312 5，0.125 0)，B2-C特征向量歸一化權重Wi=(0.200 0，0.800 0)，B3-C特征向量歸一化權重Wi=(0.466 7，0.333 3，0.200 0)。

1.3.4影響因子的權重各影響因子權重如表5所示。

1.4空間分析

本次適宜性區劃是通過將南通市主城區網格化，計算各單元格的適宜性指數，從而得到整個研究區的適宜性分區。為提高計算準確度，對南通市主城區陸地大約410 km2的范圍進行1 km × 1 km的網格剖分，參與計算評價的單元格個數為466個。

對評價因子分區圖進行矢量化，并對各評價因子的分區進行屬性賦值，通過MapGIS軟件的空間分析功能使各單元格獲得各評價因子的屬性，根據式(4)計算各單元格的適宜性分值。

(4)

式(4)中,Gi為適宜性分值，Wi為第i個因子的權重(表5)，Si為單元格第i個因子的賦值，m=7。

應用MapGIS軟件空間分析功能，對影響各適宜性分區的評價因子進行對應區屬性的區對區合并分析，再根據表5所列的各評價因子權重，對各適宜性分區迭加后的值進行權重計算。

對于?Ai∈D，檔案實體A的特征信息分屬1個或多個視圖，在構建檔案信息關聯分析模型之前，需要對來自不同數據源異構的檔案數據進行數據融合，通過對檔案信息進行特征提取形成關聯節點，將多源異構的檔案數據映射到不同視圖的特征關系圖中，形成多視圖下的檔案信息特征關聯圖。

表5　各影響因子的權重

2淺層地溫能開發利用適宜性分區

2.1適宜性分區標準

結合各評價因子分區套圖信息，將南通市主城區地下水地源熱泵適宜性分區標準分為不適宜、較適宜、適宜3類(表6)，并通過區屬性重新定義分區，獲得研究區淺層地溫能開發利用適宜性分區圖(圖3)。

表6　南通市主城區地下水地源熱泵適宜性分區標準

圖3　南通市主城區地下水地源熱泵適宜性分區圖Fig.3　Map showing suitability zoning of groundwater source heat pump system in the main urban area of Nantong

2.2適宜性分區

運用MapGIS對適宜性分區圖進行統計，得出南通市主城區地下水地源熱泵系統適宜區面積246.9 km2，較適宜區面積154.0 km2，占總調查面積的97.78%。適宜區水量豐富，滲透條件較好，便于熱量交換，有利于提高地埋管式地源熱泵的換熱效率。研究區地下水地源熱泵適宜區主要集中在第四紀長江古河道活動區，沉積物以河床相砂層為主，厚度較大，富水性好，回灌能力強，水文地質條件優越，非常適宜做地下水地源熱泵工程；較適宜區沉積物較復雜，主要由黏性土和砂層組成的雙層結構，富水性和回灌能力均不及適宜區；不適宜區狼山、劍山、軍山為孤山殘丘，長期遭受風化剝蝕，其余地區為長江河道現代沉積區，沉積物以黏性為主，砂層不發育，富水性差，不適合建立水源熱泵工程。結合目前南通市主城區已建成的地下水地源熱泵的實例，對本次適宜性分區進行檢驗，驗證了分區結果的準確性，證明本次工作評價體系、評價方法基本正確，評價結果可信，可以作為南通市主城區建立地下水地源熱泵的依據。

3結論

(1) 對江蘇省沿海地區淺層地溫能適宜性進行了分區，并獲得了南通市主城區的地下水地源熱泵適宜性分區結果，為南通地區淺層地溫能合理開發利用提供了科學依據，對其他地區的分區具有指導意義。

(3) 地下水地源熱泵的適宜性主要考慮的是良好的水文地質條件，應綜合考慮多種因素(包括非地質因素)來決定該地區的適宜性評價方式，如單井的出水量和回灌量。本次評價取用采灌比分區作為評價因子，研究區采灌比差別較小，對研究區分區評價影響不大。

參考文獻：

董殿偉，江劍，王立發，等.2010.北京市平原區地源熱泵系統適宜性分區評價.北京水務，(2):12-14.

DZ/T 0225—2009,淺層地熱能勘查評價規范.

韓再生，冉偉彥，佟紅兵，等. 2007.淺層地熱能勘查評價.中國地質，34(6):1115-1121.

江蘇省地質礦產局第一水文地質工程地質大隊.1989.南通市水文地質、工程地質、環境地質綜合勘察報告(1∶5萬).江蘇南京：江蘇省地質礦產局第一水文地質工程地質大隊.

江蘇省地質工程勘察院.2008. 南通市地下水資源評價報告.江蘇南京：江蘇省地質工程勘察院.

單衛華，邢衛兵. 2008.南通沿海地區主采層地下水咸化特征.地質學刊，32(4):286-291.

王貴玲，藺文靜，張薇. 2012. 我國主要城市淺層地溫能利用潛力評價.建筑科學，28(10):1-8.

苑雷，陳剛. 2013.改進層次分析法下的淺層地溫能適宜性評價.可再生能源，31(7):112-115.

楊澤，于慧明，都基眾. 2013. 基于層次分析法的沈陽市地下水地源熱泵適宜性研究.建筑節能，41(12):32-35.

CHANDIO I A， MATORI A N B， YUSOF K B W， et al. 2013.GIS-based analytic hierarchy process as a multicriteria decision analysis instrument: a review. Arabian Journal of Geosciences， 6(8): 3059-3066.SAATY T L. 1986. Axiomatic foundation of the analytic hierarchy process. Management Science， 32(7): 841-855.

GIS/AHP-based suitability assessment of groundwater source heat pump system in the main Nantong urban area

ZHOU Hongwei, WANG Qi, ZHANG Shasha, YANG Yi, GUO Hua

(Environmental Geology Exploration Institute of Jiangsu Province, Nanjing 210002, Jiangsu, China)

Abstract：This study chosed the geological and hydrological conditions, groundwater dynamic field, temperature field, thermal properties, water chemistry field and geological environment of the main urban area of Nantong as evaluation factors of the Analytic Hierarchy Process (AHP) model, and obtained the weight values of each factor according to the AHP method. Combining with the spatial analysis technology of GIS, we classified the suitability of groundwater source heat pump system in the study area into three types, i.e., suitable, fairly suitable and unsuitable types, with the suitable and fairly suitable areas accounting for 97.78% of the total study area. The further engineering inspection suggests that this evaluation method is suitable for the study area.

Keywords：Analytic Hierarchy Process(AHP); groundwater source heat pump; MapGIS spatial analysis; division of suitability; Nantong in Jiangsu

doi:10.3969/j.issn.1674-3636.2016.01.178

收稿日期：2015-03-09；修回日期：2015-05-18；編輯：陸李萍

基金項目：江蘇省地質勘查基金項目(蘇國土資[2012]999號)

作者簡介：周紅衛(1970—)，男，高級工程師，地下水科學與工程專業，長期從事工程地質、水文地質、環境地質及地熱地質研究工作，E-mail： zsscug@163.com

中圖分類號：P314.9

文獻標識碼：A

文章編號：1674-3636(2016)01-0178-05