紀軼群，郭高軒，2，李文忠，沈媛媛，許 亮，南英華，陸海燕

（1.北京市水文地質(zhì)工程地質(zhì)大隊，北京 100195；2.中國科學(xué)院大學(xué)地球科學(xué)學(xué)院，北京 100049；3.長慶油田分公司安全環(huán)保監(jiān)督部第九安全環(huán)保監(jiān)督站，內(nèi)蒙古烏審旗 017300）

0 引言

城市化進程加快、人口急劇增長，加之氣候的不確定性等因素使得北京市的供水矛盾日益嚴峻。平原區(qū)第四系地下水連年超采致使水位持續(xù)下降、漏斗面積不斷擴大[1～3]，地面沉降持續(xù)發(fā)展[4]。因此，以往被認為“戰(zhàn)備水”的巖溶水資源勘查于2011年正式啟動，項目是北京市于南水北調(diào)進京前開展的規(guī)模最大、手段最全的地下水資源勘查工程。

以往研究表明，北京地區(qū)的巖溶水不僅具有北方巖溶水的系統(tǒng)開放度高、發(fā)育規(guī)模大、系統(tǒng)資源要素構(gòu)成多的特點，而且各水資源要素間轉(zhuǎn)化關(guān)系復(fù)雜[5～6]。區(qū)域巖溶水資源量計算需要大量地理環(huán)境、基礎(chǔ)地質(zhì)、水文地質(zhì)等數(shù)據(jù)作為基礎(chǔ)，該數(shù)據(jù)呈現(xiàn)為海量、時空動態(tài)等特性，常規(guī)手段很難達到及時和有效的存儲、查詢、修改、建模、計算等工作。

地理信息系統(tǒng)(簡稱GIS)是一種集采集、存儲、管理、分析、顯示等功能于一體的計算機系統(tǒng)，是分析和處理海量空間數(shù)據(jù)的通用技術(shù)。GIS已被廣泛地應(yīng)用于資源調(diào)查、環(huán)境評估、水資源管理等領(lǐng)域[7～9]。因此，利用GIS技術(shù)，有效地對巖溶水資源勘查與開發(fā)及周圍地質(zhì)環(huán)境信息進行統(tǒng)一的數(shù)據(jù)庫管理，并在數(shù)據(jù)庫上層，構(gòu)建巖溶水資源量評價模型，可以在時間維和空間維上準確、高效的

計算和評價北京巖溶水資源量的分布。為水文地質(zhì)勘查、地下水的合理開采、高效利用、地質(zhì)環(huán)境監(jiān)測與保護、科學(xué)管理決策提供有力支持。

1 模塊框架設(shè)計

（1）均衡法介紹

水量均衡法是評價地下水資源的根本方法。方法原理嚴格，計算結(jié)果準確，適宜各個地區(qū)。應(yīng)用水量均衡法可以評價區(qū)域地下水補給資源量，初步確定可開采資源量。對于地下水系統(tǒng)來說，在補給和消耗的不平衡發(fā)展過程中，一個區(qū)域內(nèi)的含水層系統(tǒng)，在任一時間段（△t）的補給量和消耗量之差等于該區(qū)域中水體積的變化量，這就是水量均衡法的基本原理[10]。

根據(jù)基本原理可建立天然條件下水量均衡方程式：

Q補-Q排=Q儲或 Q補=Q排+Q儲

式中：Q補為計算期間內(nèi)地下水系統(tǒng)各補給量總和；Q排為計算期間內(nèi)地下水系統(tǒng)各排泄量總和；Q儲為計算期間內(nèi)地下水系統(tǒng)儲存量的變化量。

（2）模塊總體結(jié)構(gòu)設(shè)計

本計算模塊基于C/S（客戶端/服務(wù)器）技術(shù)體系結(jié)構(gòu)設(shè)計及采用COM組件開發(fā)，COM組件具有封裝性好，接口開放的特點，具體選用ESRI的ArcObject組件，基于C/S技術(shù)體系可以滿足“北京巖溶水資源計算模塊中”所需的資料量多、數(shù)據(jù)類型多、分析較復(fù)雜的特點，支持空間多元數(shù)據(jù)庫的特點。

根據(jù)水文地質(zhì)特點及區(qū)域巖溶水資源量計算評價的需要,北京巖溶水資源量計算模塊由空間數(shù)據(jù)處理系統(tǒng)、屬性數(shù)據(jù)處理系統(tǒng)、影像數(shù)據(jù)處理系統(tǒng)、區(qū)域資源補給項計算子模塊，區(qū)域資源排泄項計算子模塊、區(qū)域資源變化量子模塊、區(qū)域可采資源量評價模塊、資源量專題圖制作等功能組成。系統(tǒng)總體結(jié)構(gòu)（如圖1所示）采用了成熟的C/S結(jié)構(gòu)，方便對數(shù)據(jù)和功能的分類管理，同時安全性和穩(wěn)定性比較高。

圖1 模塊總體結(jié)構(gòu)

（3）模塊采用的關(guān)鍵技術(shù)

系統(tǒng)構(gòu)建采用ArcGIS平臺，ArcGIS是由美國ESRI公司研制開發(fā)的一個功能強大的通用性地理信息系統(tǒng)軟件，該軟件最顯著的特點是對地理信息本身的生成與處理能力較強，且支持拓撲結(jié)構(gòu)，可進行空間分析、緩沖區(qū)分析和疊置分析，對處理巨量數(shù)據(jù)有較大優(yōu)勢，加之其提供的面向?qū)ο蟪绦蛟O(shè)計開發(fā)工具AO、及VBA，使用戶能開發(fā)出自己需要的功能模塊，并嵌入到ArcGIS中；同時它還支持DLL（動態(tài)連接庫）和DDE（動態(tài)數(shù)據(jù)交換）技術(shù)，進一步提高了ArcGIS平臺的開放性[11]。因此，可以在此基礎(chǔ)上二次開發(fā)出適宜北京巖溶水資源量評價的專業(yè)化客戶軟件。

Geodatabase(空間數(shù)據(jù)庫)與城市水資源管理集成有著明顯的優(yōu)勢[12]。本次空間數(shù)據(jù)管理采用Geodatabase(空間數(shù)據(jù)庫)統(tǒng)一管理各種類型空間數(shù)據(jù)。Geodatabase是一種采用標準關(guān)系數(shù)據(jù)庫技術(shù)來表現(xiàn)地理信息的數(shù)據(jù)模型，支持在標準的數(shù)據(jù)庫管理系統(tǒng)(DBMS)表中存儲和管理地理信息，并且在很多方面增強了關(guān)系型數(shù)據(jù)庫。Geodatabase支持多種DBMS結(jié)構(gòu)和多用戶訪問，且大小可伸縮。多用戶Geodatabase通過ArcSDE支持多種數(shù)據(jù)庫平臺。

2 模塊數(shù)據(jù)庫構(gòu)成

北京巖溶水資源量均衡法評價模塊數(shù)據(jù)庫由巖溶水水資源空間數(shù)據(jù)庫（包括矢量空間數(shù)據(jù)庫和柵格空間數(shù)據(jù)庫）、屬性數(shù)據(jù)庫、影像數(shù)據(jù)庫、元數(shù)據(jù)庫構(gòu)成（如圖2，表1），各子庫間數(shù)據(jù)相互獨立，同時邏輯關(guān)系又相互交織在一起，共同為整個信息系統(tǒng)提供基礎(chǔ)數(shù)據(jù)?？臻g數(shù)據(jù)庫以WGS-84為基準坐標系統(tǒng)。

圖2 評價模塊數(shù)據(jù)庫結(jié)構(gòu)關(guān)系

空間數(shù)據(jù)庫由工作過程中建立的單要素圖層和綜合圖件組成。單要素圖層是反映巖溶區(qū)的基礎(chǔ)地理、地質(zhì)、水文地質(zhì)條件的圖形要素集，以點、線、面、體、柵格等形式存儲。它是成果圖件的基本組成單元。通過疊加地層圖層、構(gòu)造圖層、自然地理等圖層形成地質(zhì)圖；通過疊加富水性分區(qū)、巖溶水系統(tǒng)分區(qū)等圖層形成綜合水文地質(zhì)圖，因此，采取圖層疊加的方式組合成綜合圖件。

屬性數(shù)據(jù)庫包括一般事務(wù)性描述性屬性和圖形數(shù)據(jù)內(nèi)部屬性的擴展，在圖形內(nèi)部屬性結(jié)構(gòu)中，為了減小數(shù)據(jù)冗余，實現(xiàn)數(shù)據(jù)快速調(diào)用，內(nèi)部屬性只記錄了地理實體數(shù)據(jù)的基本信息，詳細的序列記錄建立在其他數(shù)據(jù)表中，例如觀測井分布圖中，內(nèi)部屬性信息只包括觀測井的統(tǒng)一編號，地理位置、坐標數(shù)據(jù)，而日常動態(tài)觀測數(shù)據(jù)則建立在一般性屬性表中，空間數(shù)據(jù)庫與屬性庫記錄通過統(tǒng)一編號進行關(guān)聯(lián)（如圖3），通過編程實現(xiàn)圖形和屬性間互動查詢。

影像數(shù)據(jù)庫包括各巖溶區(qū)航空遙感影像共9940km2。由于航片影像文件較大，為了便于后期使用，在入庫前對各航片采用影像金字塔技術(shù)，提高后期瀏覽效率。

表1 北京巖溶水資源量評價模塊數(shù)據(jù)庫組成和結(jié)構(gòu)

元數(shù)據(jù)庫，元數(shù)據(jù)是“關(guān)于數(shù)據(jù)的數(shù)據(jù)”，由于巖溶水資源量評價涉及眾多空間數(shù)據(jù)、屬性數(shù)據(jù)，需要對地理、地質(zhì)、社會數(shù)據(jù)集的內(nèi)容、質(zhì)量、表示方式、空間參考、管理方式以及數(shù)據(jù)集的其他特征進行描述，才能正確的管理、編輯、使用數(shù)據(jù)；它也是實現(xiàn)地理空間信息共享的核心標準之一。

3 模塊功能構(gòu)建

實現(xiàn)均衡法巖溶水資源計算，即每個子模塊應(yīng)根據(jù)建立的均衡方程式實現(xiàn)每一項的計算功能；首先按照北京地區(qū)巖溶水資源源匯項特點，構(gòu)建北京巖溶系統(tǒng)水均衡項，如公式1，后按照各源匯項進行子模塊的構(gòu)建。

式中：

補給項：

Q入——大氣降水入滲補給量（萬m3/a）

Q回——農(nóng)灌水回滲量（萬m3/a）

Q渠——地表水滲漏補給量（萬m3/a）

Q越——越流量（萬m3/ a）

Q側(cè)入——地下水側(cè)向流入量（萬m3/a）

消耗項：

Q開采——開采水量（萬m3/a）

Q泉——泉水流量（萬m3/a）

Q燕——地下水蒸發(fā)量（萬m3/a）

Q地——河水基流（萬m3/a）

Q第四系——向第四系頂托越流補給量（萬m3/a）

Q測出——地下水側(cè)向流出量（萬m3/a）

其它參數(shù)含義：

μ(s)——均衡區(qū)內(nèi)潛水水位變幅帶內(nèi)含水層給水度（μ）及承壓水彈性釋水系數(shù)（s）

F——均衡區(qū)面積（m2）

△H——均衡時段內(nèi)地下水位變幅（m）

△t——均衡時段（a）

3.1 補給量計算功能

（1）大氣降水入滲補給量

該模塊先進行區(qū)域降雨量計算，按照觀象臺位置數(shù)據(jù)生成泰森多邊形（Thiessen polygon），控制不同區(qū)域降水量。可溶巖大氣降水入滲計算使用Q入=F?P?α×10-4（大氣降水入滲補給量（萬m3/a）；F計算面積（m2）；P降水量（m）；α入滲系數(shù)）公式。模塊中可以量測點與點的距離和計算多邊形的面積F，可以按照巖溶水各等級系統(tǒng)來統(tǒng)計入滲量。

圖3 空間數(shù)據(jù)庫與屬性數(shù)據(jù)庫之間的鏈接

（2）農(nóng)業(yè)灌溉回滲量

該模塊實現(xiàn)農(nóng)業(yè)灌溉回滲量計算按照Q回=Q農(nóng)灌×β（Q回為農(nóng)灌回滲量（萬m3/a），Q農(nóng)灌地下水（或地表水）灌溉水量（萬m3/a），β地下水（或地表水）灌溉回歸系數(shù)））公式。灌溉面積可以在模塊中由土地利用數(shù)據(jù)獲取各系統(tǒng)裸露山區(qū)園地和耕地面積，保持了計算集成性。

（3）地表水滲漏補給量

根據(jù)河流不同測站的徑流量來計算此段河流的入滲量。如子模塊中可以查詢永定河上游雁翅流量站徑流量與下游三家店流量站徑流量，兩者相減，加入其他修正項后，可算出永定河從雁翅站到三家店段的永定河滲漏量。

（4）越流量

第四系地下水水位高于下伏巖溶水水位，存在著垂向交換量，以越流方式補給巖溶水，該子模塊實現(xiàn)計算通過V=KI×10-4=K（HA-HB）×10-4/M（V為單位面積越流量（萬m/d）；K為弱透水層垂向滲透系數(shù)（m/d）；I為兩個含水層之間的水力梯度；M為弱透水層厚度（m）；HA為上覆第四系含水層水位（m）；HB為下伏含水層水位（m））公式，在計算的數(shù)據(jù)結(jié)構(gòu)上可以抽象為不同柵格圖層的計算，即使用地圖代數(shù)進行多個柵格層之間的數(shù)學(xué)運算。

（5）地下水側(cè)向徑流量

該子模塊采用達西定律計算實現(xiàn)Q側(cè)入=K?I?H?B×365×10-4（Q側(cè)入側(cè)向流入量（萬m3/d），K含水層滲透系數(shù)（m/d），I水力坡度，H含水層厚度（m），B斷面寬度（m）），此模塊中提供巖溶水系統(tǒng)入流邊界的地質(zhì)剖面圖。因此，可以測量含水層寬度和厚度，方便計算不同入流斷面來水量，同時根據(jù)地下水水位等值線來計算水力梯度。

3.2 排泄量計算功能

（1）開采量

功能上實現(xiàn)基巖井及水源地集中井，按照不同時間段、不同行政區(qū)域、不同巖溶水系統(tǒng)，統(tǒng)計出開采量，支持Excel格式數(shù)據(jù)導(dǎo)入和導(dǎo)出。

（2）泉水流量

子模塊實現(xiàn)泉點位置添加、編輯、刪除，對流量表的查詢、修改并按照不同行政區(qū)域，不同巖溶水系統(tǒng)給出某個時間段的統(tǒng)計監(jiān)測泉流量。

（3）蒸發(fā)量

當巖溶水水埋深小于極限蒸發(fā)埋深時，地下水將產(chǎn)生蒸發(fā)，在功能實現(xiàn)上子模塊采用Q蒸=ε?F?10-4（Q蒸為地下水蒸發(fā)強度（萬m3/a），ε為蒸發(fā)系數(shù)，F(xiàn)為面積m2）；蒸發(fā)系數(shù)ε通過ε=ε水蒸?η?(1-h/h0）2，(ε水蒸為氣象站蒸發(fā)值（mm）;h為地下水水位埋深（m），h0為蒸發(fā)極限埋深（m），η為折算系數(shù))，在功能中讀取各氣象站蒸發(fā)量值，可以根據(jù)不同區(qū)域設(shè)定蒸發(fā)極限深度和折算系數(shù)。

（4）巖溶水向外排泄

巖溶水向外排泄包括巖溶水向河流排泄（河流基流），向第四系孔隙水排泄（第四系頂托越流補給），側(cè)向流出量。其中巖溶水向河流排泄功能實現(xiàn)與補給項中地表水滲漏補給量相同，第四系孔隙水排泄功能實現(xiàn)與補給項中越流補給量相同，側(cè)向流出量功能實現(xiàn)與補給項中地下水側(cè)向徑流量相同，在此不在贅述。

3.3 儲存變化量計算功能

儲存變化量功能按照公式△Q=μF△h（μ為給水度或者釋水系數(shù)，F(xiàn)為面積m2，△h為時段內(nèi)水頭變化量單位m）來計算，μ和△h可以支持各圖層參數(shù)分區(qū)柵格計算，從而計算出變化量與前面計算的總補給量減去總排泄量之差進行誤差分析，一般認為兩者計算總誤差在5%以內(nèi)時，各項量計算比較合理，結(jié)果可信度高。

3.4 可采資源量評價功能

地下水可開采量評價一般用地下水總補給量減去不可襲奪的排泄量來確定，不可襲奪的排泄項主要有河流的基流、向第四系排泄量、側(cè)向流出量等。在功能實現(xiàn)上子模塊可以設(shè)定不可襲奪項及其量的大小。對于不可襲奪排泄量的設(shè)定，可以根據(jù)地區(qū)水資源情況，評價出不同生態(tài)條件下的巖溶水開采量，如河流基流，泉流量排泄量僅為以往的一半時的區(qū)域可開采量。

3.5 計算結(jié)果的專題圖制作功能

專題圖就是根據(jù)相應(yīng)的專題數(shù)據(jù)和指標，采用一定的算法，對空間要素進行分類分級，并以一定的分類分級符號進行可視化展示，從而使用戶能夠以直觀的形式，了解某專題的空間分布和定性、定量特征[13]。北京巖溶水資源量計算結(jié)果的專題圖制作子模塊應(yīng)能夠為用戶提供最直觀的可視化表達。系統(tǒng)提供北京各巖溶水系統(tǒng)基于各種補給量、排泄量、資源量的簡單渲染專題圖、唯一值專題圖、點密度專題圖、數(shù)據(jù)分級專題圖、餅狀專題圖。例如，北京各巖溶系統(tǒng)巖溶水可采資源量柱狀專題圖，來定量渲染各區(qū)不同資源量，給使用者直觀感受同時便于理解。

4 模塊應(yīng)用實例

以西山巖溶水系統(tǒng)中的十渡—長溝子系統(tǒng)為例，使用模塊進行資源量計算，研究區(qū)在房山西南部面積為636.9km2。主要巖溶含水層為薊縣系霧迷山組（Jxw）和鐵嶺組（Jxt），拒馬河從西向東穿過研究區(qū)，同時有巖溶大泉分布其中。巖溶系統(tǒng)內(nèi)部均衡項多、數(shù)據(jù)量大，計算復(fù)雜。巖溶水補給來源由大氣降水入滲、水河渠滲漏量、農(nóng)業(yè)灌溉回歸、側(cè)向流入量構(gòu)成；巖溶水排泄項由開采量、泉流量、蒸發(fā)量、向地表排泄量、向第四系排泄量、側(cè)向流出量構(gòu)成。

圖4 模塊計算界面

表2 2011年十渡—長溝巖溶水系統(tǒng)水均衡表（萬m3）

在模塊中選擇該地區(qū)的2011年均衡要素和數(shù)據(jù)與圖層，設(shè)定合理水文地質(zhì)參數(shù)（計算界面如圖4）后得出2011年十渡－長溝巖溶水補給量和排泄量匯總以及儲存變化量見表2，兩者誤差在5%以內(nèi)，認為本次計算資源量結(jié)果較為合理。

5 結(jié)論

（1）本文采用組件GIS，結(jié)合空間數(shù)據(jù)庫（Geodatabase）技術(shù)，使用區(qū)域水均衡評價原理，構(gòu)建以均衡法為基礎(chǔ)資源量評價模塊。

（2）為了驗證模塊的可靠性，使用該模塊，結(jié)合人機交互方法對十渡—長溝巖溶水系統(tǒng)2011年資源進行了試算，結(jié)果表明區(qū)域補給量與排泄量差值與儲變量值誤差小于5%，認為本區(qū)資源量計算值較為合理。

（3）該模塊能夠滿足復(fù)雜水文地質(zhì)條件下，巖溶區(qū)域資源量計算要求，可以為北京巖溶水資源勘查工程中區(qū)域巖溶水資源量計算提供信息化技術(shù)手段，有著計算效率高，準確性好等特點。與此同時模塊對巖溶地質(zhì)條件和巖溶地下水資源量的相關(guān)數(shù)據(jù)進行空間信息化，提高各層次信息的可視程度、并全景式地顯示巖溶水資源量計算過程和時空分布特征；極大的提高地下水資源管理的效率和水平，因此建立該模塊具有很大的實際意義。

[1]張士鋒，孟秀敬等.北京市水資源與水量平衡研究[J].地理研究，2012，31（11）：1991～1996.

[2]陳 江，聶振龍等.北京市地下水資源利用現(xiàn)狀評價[J].南水北調(diào)與水利科技,2008, 6(6)：97～99.

[3]任志遠.北京市平原區(qū)地下水資源研究[D].吉林大學(xué)，2006：16～17.

[4]王 萍.北京市超量開采地下水引起的地面沉降研究[J].勘察科學(xué)技術(shù)，2004（05）：46～49.

[5]梁永平,韓行瑞等.中國北方巖溶地下水環(huán)境問題與保護[M].地質(zhì)出版社，2013：41～51.

[6]郭高軒,辛寶東等.北京巖溶水勘查開發(fā)的現(xiàn)狀與思考[J].南水北調(diào)與水利科技,2011,.9(2)：33～36.

[7]紀軼群,葉 超.北京平原區(qū)淺層地溫能空間信息系統(tǒng)設(shè)計與實現(xiàn)[J].首都師范大學(xué)學(xué)報,2009(01),54～57.

[8]諸云強，宮輝力，趙文吉等.基于組件技術(shù)的地理信息系統(tǒng)二次開發(fā)—以地下水資源空間分析系統(tǒng)為例[J].地理與地理信息科學(xué)，2003, 19(1)：16～19.

[9]武 強,鄒德禹,董東林等.塔里木盆地水資源開發(fā)管理的地理信息系[J].中國礦業(yè)大學(xué)學(xué)報，1999,28(1)：78～81.

[10]中國地質(zhì)調(diào)查局 主編.水文地質(zhì)手冊（第二版）[M].北京：地質(zhì)出版社2012：692～694.

[11]于 軍.基于ArcGIS平臺的蘇錫常地區(qū)地面沉降管理信息系統(tǒng)研究[D].吉林大學(xué)，2006：20～21.

[12]龍 瀛,杜鵬飛,趙東東等.基于Geodatabase的城市水資源管理系統(tǒng)[J].清華大學(xué)學(xué)報（自然科學(xué)版），2006, 46（9）：1560～1563.

[13]聶小波，吳北平，何保國.基于ArcGIS Engine的專題圖模塊的設(shè)計與實現(xiàn)[J].地理空間信,2006(1)：12～14.