杜 超, 肖長來, 呂 軍, 田浩然

1.松遼流域水資源保護(hù)局松遼水環(huán)境科學(xué)研究所, 長春 130021 2.吉林大學(xué)環(huán)境與資源學(xué)院, 長春 130021

地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)系統(tǒng)開發(fā)及應(yīng)用
——以下遼河平原為例

杜 超1, 肖長來2, 呂 軍1, 田浩然1

1.松遼流域水資源保護(hù)局松遼水環(huán)境科學(xué)研究所, 長春 130021 2.吉林大學(xué)環(huán)境與資源學(xué)院, 長春 130021

地下水?dāng)?shù)值模擬需要大量的監(jiān)測數(shù)據(jù)作為支撐，數(shù)據(jù)的錄入耗時(shí)耗力，存在出錯(cuò)的可能，且不能第一時(shí)間進(jìn)行實(shí)時(shí)快速預(yù)報(bào)。采用VB語言編程技術(shù),以識別和驗(yàn)證后的下遼河平原地下水流數(shù)值模型和地下水溶質(zhì)運(yùn)移數(shù)值模型為基礎(chǔ),重新編譯后的MF2K和MT3DMS為內(nèi)核計(jì)算程序,Microsoft SQL Server 2000為數(shù)據(jù)庫，開發(fā)了模型與地下水實(shí)時(shí)監(jiān)測系統(tǒng)的接口，集成建立了下遼河平原地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)模型，并開發(fā)了操作界面下的操作系統(tǒng)，實(shí)現(xiàn)了下遼河平原地下水水質(zhì)的實(shí)時(shí)預(yù)報(bào)功能。該系統(tǒng)將地下水實(shí)時(shí)監(jiān)控技術(shù)與地下水?dāng)?shù)值模擬技術(shù)結(jié)合，操作簡便、時(shí)效性強(qiáng)、出錯(cuò)概率小，能夠?qū)崟r(shí)獲取地下水水位、水質(zhì)監(jiān)測數(shù)據(jù)，對地下水流場和地下水溶質(zhì)濃度場的演化進(jìn)行快速預(yù)報(bào)。同時(shí)，本系統(tǒng)預(yù)留了更新接口，可根據(jù)逐步積累的地下水監(jiān)測數(shù)據(jù)和系統(tǒng)計(jì)算數(shù)據(jù)進(jìn)行進(jìn)一步驗(yàn)證，不斷提高系統(tǒng)預(yù)報(bào)的精度。

實(shí)時(shí)預(yù)報(bào)系統(tǒng)；數(shù)值模擬；地下水水質(zhì)；溶質(zhì)運(yùn)移；系統(tǒng)集成

下遼河平原人口眾多、工農(nóng)業(yè)發(fā)達(dá)、地下水超采嚴(yán)重，是我國人類活動最為強(qiáng)烈的區(qū)域之一[1-4]，地下水水質(zhì)在強(qiáng)烈的人類活動下也逐漸趨于惡化[5-6]。實(shí)時(shí)掌握地下水水質(zhì)狀況，并做出科學(xué)的模擬和預(yù)報(bào)，是對下遼河平原地下水資源管理工作者提出的必然要求。目前，地下水?dāng)?shù)值模擬技術(shù)已相對成熟，GMS、Visual Modflow等商業(yè)模擬軟件已被廣泛應(yīng)用于我國地下水研究和管理中[7-9]。但地下水?dāng)?shù)值模擬需要大量的監(jiān)測數(shù)據(jù)，數(shù)據(jù)的錄入耗時(shí)耗力，存在出錯(cuò)的可能，且不能第一時(shí)間進(jìn)行實(shí)時(shí)快速預(yù)報(bào)；同時(shí)要求操作者具有專業(yè)的水文地質(zhì)基礎(chǔ)知識作為技術(shù)支撐，不便于基層工作者開展工作。因此，將地下水實(shí)時(shí)監(jiān)控技術(shù)與地下水?dāng)?shù)值模擬技術(shù)結(jié)合，建立一套操作簡便、時(shí)效性強(qiáng)、出錯(cuò)率低的地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)系統(tǒng)，對于保障下遼河平原地下水水質(zhì)安全十分必要。

1 地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)原理

在地下水流數(shù)值模型及地下水溶質(zhì)運(yùn)移模型的基礎(chǔ)上，以Fortran修改編譯后的MF2K(Modflow2000)和MT3DMS為內(nèi)核計(jì)算程序，采用VB語言進(jìn)行系統(tǒng)集成，開發(fā)其與地下水實(shí)時(shí)監(jiān)測數(shù)據(jù)庫的接口，實(shí)現(xiàn)地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)功能。

1.1 內(nèi)核計(jì)算程序

MF2K和MT3DMS是分別用于模擬多孔介質(zhì)中三維地下水流和溶質(zhì)運(yùn)移的計(jì)算機(jī)程序，采用有限差分法進(jìn)行求解，是世界上兩大地下水?dāng)?shù)值模擬軟件GMS和Visual Modflow的水流模型和水質(zhì)模型內(nèi)核程序，已被廣泛應(yīng)用于世界各地的地下水?dāng)?shù)值模擬評價(jià)中。在近幾年剛剛完成的中國北方平原區(qū)地下水資源評價(jià)項(xiàng)目中，中國北方地區(qū)的主要平原區(qū)和盆地的地下水均進(jìn)行了數(shù)值模擬，采用的水流模型基本都為MF2K；MT3DMS已在實(shí)際應(yīng)用中被不斷驗(yàn)證，取得了良好的模擬效果。MF2K和MT3DMS是當(dāng)今世界上最為流行及可靠的地下水水質(zhì)模擬程序，源程序及源代碼完全公開[10-11]。

根據(jù)實(shí)時(shí)預(yù)報(bào)系統(tǒng)的需要，采用Fortran語言對MF2K和MT3DMS內(nèi)核計(jì)算程序進(jìn)行了部分修改及編譯，作為水流模型和溶質(zhì)運(yùn)移模型的計(jì)算內(nèi)核程序。

1.2 實(shí)時(shí)監(jiān)測數(shù)據(jù)讀取

地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)需要的實(shí)時(shí)數(shù)據(jù)主要包括降水、蒸發(fā)等實(shí)時(shí)氣象監(jiān)測數(shù)據(jù)，地下水位、地下水質(zhì)等實(shí)時(shí)地下水監(jiān)測數(shù)據(jù)，供水量、用水量等實(shí)時(shí)供用水監(jiān)測數(shù)據(jù)等。

系統(tǒng)采用VB中的ADO控件與地下水實(shí)時(shí)監(jiān)測數(shù)據(jù)庫連接，獲取實(shí)時(shí)監(jiān)測數(shù)據(jù)。

1.3 地下水流場實(shí)時(shí)預(yù)報(bào)

地下水流場實(shí)時(shí)預(yù)報(bào)模型的構(gòu)建基礎(chǔ)是利用MF2K建立的地下水流數(shù)值模型。獲取實(shí)時(shí)監(jiān)測數(shù)據(jù)或手動輸入數(shù)據(jù)后，系統(tǒng)對MF2K水流模型的輸入文件進(jìn)行修改，并對文件中的相關(guān)變量進(jìn)行實(shí)時(shí)更新，則可完成對實(shí)時(shí)預(yù)報(bào)模型的前處理過程；然后調(diào)用MF2K內(nèi)核程序進(jìn)行計(jì)算，對地下水流場進(jìn)行實(shí)時(shí)預(yù)報(bào)。地下水流場實(shí)時(shí)預(yù)報(bào)模型相關(guān)輸入文件見表1，運(yùn)行流程見圖1。

圖1 地下水流場實(shí)時(shí)預(yù)報(bào)模型運(yùn)行流程圖Fig. 1 Running process of real time groundwater flow prediction model

1.4 地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)

地下水流場實(shí)時(shí)預(yù)報(bào)模型的構(gòu)建基礎(chǔ)是利用MT3DMS建立的地下水溶質(zhì)運(yùn)移模型。由于水質(zhì)模型是依賴于水流模型的，因此，水質(zhì)模型的時(shí)空離散及計(jì)算方案不再單獨(dú)另設(shè)，而與水流模型同步。提取標(biāo)準(zhǔn)輸入文件后，根據(jù)應(yīng)力期的設(shè)置情況對標(biāo)準(zhǔn)輸入文件進(jìn)行實(shí)時(shí)更新，隨后鏈接水流模型MF2K生成的水位鏈接文件(.FTL)，調(diào)用MT3DMS模型進(jìn)行實(shí)時(shí)預(yù)報(bào)。地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)模型相關(guān)輸入文件見表2，運(yùn)行流程見圖2。

表1 地下水流場實(shí)時(shí)預(yù)報(bào)模型輸入文件及用途

Table 1 Input files and functions of real time groundwater flow prediction model

文件文件類型文件用途.namNAME文件定義程序運(yùn)行所需輸入文件和輸出文件的名稱和路徑.dis離散文件定義模型離散狀況及含水層、地下水流的性質(zhì).zon分區(qū)文件定義模型參數(shù)分區(qū).ba6基本屬性文件定義活動單元格、初始水位.oc輸出設(shè)置文件定義是否輸出結(jié)果與保存相關(guān)設(shè)置.lpf層特征文件定義滲透系數(shù)、給水度、儲水系數(shù)等水文地質(zhì)參數(shù).riv河流文件定義河流與地下水轉(zhuǎn)化過程中的參數(shù)及變量值.rch面狀補(bǔ)給率文件定義面狀補(bǔ)給源匯項(xiàng).wel井文件定義井流量.chd定水頭文件定義定水頭邊界.sip強(qiáng)隱式法文件定義強(qiáng)隱式求解法相關(guān)參數(shù)

表2 地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)模型輸入文件及用途

Table 2 Input files and functions of real time groundwater quality prediction model

文件文件類型文件用途.mtsNAME文件定義程序運(yùn)行所需輸入文件和輸出文件的名稱和路徑.btn基本溶質(zhì)運(yùn)移文件定義模型網(wǎng)格剖分、時(shí)間離散及初始濃度.adv對流文件定義對流過程求解的一些參數(shù)設(shè)置.dsp彌散文件定義彌散過程求解所需的一些參數(shù),如彌散度.gcgGCG文件定義GCG求解法的一些參數(shù)設(shè)置.lmt鏈接文件定義與MF2K模型的鏈接設(shè)置

圖2 地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)模型運(yùn)行流程圖Fig. 2 Running process of real time groundwater quality prediction model

2 下遼河平原概況

2.1 地理位置

下遼河平原地處遼河中下游，沿北東至南西方向分布在遼寧省中部地區(qū)，縱向長度約為240 km，寬度為120～140 km，平原區(qū)面積約21 929 km2。下遼河平原在行政區(qū)劃上包括了遼寧省中部城市群，是東北地區(qū)規(guī)模最大的區(qū)域一體化經(jīng)濟(jì)區(qū)，同時(shí)也是遼寧省的重要產(chǎn)糧區(qū)。

2.2 氣象水文

下遼河平原屬溫帶半濕潤、半干旱季風(fēng)氣候，四季區(qū)分明顯，雨季集中，日照時(shí)間長。多年平均氣溫為7.1～8.9℃，氣溫北高南低；降水量由南東向北西遞減，多年平均降水量為500～750 mm。

研究區(qū)內(nèi)水系發(fā)達(dá)，主要?jiǎng)澐譃檫|河--雙臺子河水系、渾河--太子河水系、大凌河--小凌河水系。

2.3 水文地質(zhì)

下遼河平原地下水資源較為豐富，平原上部廣泛分布第四系松散巖類孔隙含水層，第四系之下普遍分布上第三系館陶組和明化鎮(zhèn)組的砂礫巖、砂巖和泥巖互層的地層。本次研究對象為第四系松散巖類孔隙水，其含水層厚度大、分布穩(wěn)定、水量豐富、開采方便，是研究區(qū)內(nèi)最主要的地下水類型。

下遼河平原地下水受天然因素和人為因素的共同影響，從補(bǔ)給區(qū)、徑流區(qū)到排泄區(qū)形成了不同的地下水水化學(xué)類型，存在著較為明顯的水平分帶特征。地下水從補(bǔ)給區(qū)流向排泄區(qū)的過程，就是地下水化學(xué)組分的演化過程：降水?dāng)y帶的組分不斷減少，巖石當(dāng)中的組分不斷溶濾進(jìn)入地下水，礦化度逐漸升高；水化學(xué)類型由HCO3-Ca型逐漸轉(zhuǎn)變?yōu)镃l-Na型。

3 下遼河平原地下水流數(shù)值模型

3.1 水文地質(zhì)概念模型

本次研究針對的是下遼河平原第四系松散巖類含水層，含水層的厚度和巖性在研究區(qū)范圍內(nèi)均有不同程度的變化，因此將其概化為非均質(zhì)各向同性含水層。水流特征概化為準(zhǔn)三維非穩(wěn)定達(dá)西流。

研究區(qū)北、東、西三面為山前地帶，概化為二類流量邊界；研究區(qū)南部為濱海三角洲地區(qū)，在靠近入?？诘膫?cè)向邊界地下水位趨近于0，概化為一類水頭邊界。潛水含水層頂部為水量交換邊界；含水層底部為相對隔水邊界。

3.2 數(shù)學(xué)模型離散

根據(jù)研究區(qū)的水文地質(zhì)概念模型，建立第四系松散巖類孔隙含水層的非均質(zhì)各向同性非穩(wěn)定流的數(shù)學(xué)模型。

研究區(qū)總面積為21 929 km2，對其進(jìn)行矩形剖分，剖分為230列、221行，共30 808個(gè)網(wǎng)格單元，每個(gè)網(wǎng)格長1 000 m、寬1 000 m，面積為1 km2。采用2000年1月到2007年1月的地下水長觀資料進(jìn)行識別和驗(yàn)證，應(yīng)力期設(shè)置為1月，時(shí)間步長為15 d。故將整個(gè)模擬期離散為84個(gè)應(yīng)力期，每個(gè)應(yīng)力期劃分為2個(gè)時(shí)間段。

3.3 模型的識別和驗(yàn)證

將模擬期內(nèi)各源匯項(xiàng)及參數(shù)初值輸入模型后，運(yùn)行至模擬期末刻。將模型計(jì)算水位與實(shí)測水位進(jìn)行擬合，并對擬合誤差進(jìn)行統(tǒng)計(jì)。結(jié)果表明，水位擬合誤差小于0.5 m的結(jié)點(diǎn)數(shù)達(dá)到已知水位結(jié)點(diǎn)數(shù)的80 %以上，計(jì)算水位與實(shí)測水位擬合程度良好。典型站點(diǎn)地下水水位擬合情況見圖3。

圖3 地下水水位擬合圖Fig. 3 Fitting chart of groundwater level

模擬期涵蓋了豐、平、枯等不同代表性水文年，通過長時(shí)間序列地下水觀測數(shù)據(jù)識別和驗(yàn)證出來的模型具有較高的仿真性，能反映出下遼河平原第四系含水層的實(shí)際特征，通過反演所確定的水文地質(zhì)參數(shù)是可靠的。

4 下遼河平原地下水溶質(zhì)運(yùn)移模型

4.1 概念模型

礦化度反映水中含鹽量的多少，是區(qū)分咸水和淡水的重要標(biāo)志，也是地下水環(huán)境質(zhì)量評價(jià)中的一個(gè)重要指標(biāo)，因此選擇礦化度作為本次模擬的模擬因子。將礦化度處理為一種溶質(zhì)，其溶質(zhì)運(yùn)移過程符合對流-彌散原理，且彌散作用符合Fick定律[12-13]。

側(cè)向邊界均概化為濃度邊界。對南部入海邊界，取海水平均濃度為濃度邊界；對其他側(cè)向邊界，依據(jù)繪制的濃度等值線分別進(jìn)行賦值。對研究區(qū)內(nèi)的河流內(nèi)邊界，將其處理為地下水中礦化度的一個(gè)外源，其礦化度依據(jù)河水的平均成分進(jìn)行賦值。

垂向上，將大氣降水、地表灌溉用水?dāng)y帶的礦化度處理為地下水中相應(yīng)組分的外源；蒸發(fā)濃縮作用在潛水面發(fā)生，由模型根據(jù)蒸發(fā)量自動計(jì)算。含水層底部概化為隔水邊界。

4.2 數(shù)學(xué)模型

地下水準(zhǔn)三維非穩(wěn)定流溶質(zhì)運(yùn)移的可混溶對流-彌散數(shù)學(xué)模型為

地下水溶質(zhì)運(yùn)移對流-彌散模型的求解需要地下水流模型的支持。因此，對計(jì)算區(qū)的空間離散及時(shí)間離散沿用水流模型的設(shè)置。

4.3 模型的識別和驗(yàn)證

將源匯項(xiàng)、初始濃度場、初始參數(shù)輸入模型后，運(yùn)行模型并進(jìn)行擬合，直至擬合程度符合要求為止。礦化度擬合情況見圖4。

通過擬合結(jié)果可以看出，計(jì)算礦化度場與實(shí)測礦化度場的擬合程度較好，精度較高，說明所建立的模型能夠真實(shí)和客觀地反映研究區(qū)地下水的溶質(zhì)運(yùn)移規(guī)律，通過反演所確定的水質(zhì)模型參數(shù)是可靠的。

a. 實(shí)測礦化度場；b.計(jì)算礦化度場。圖4 地下水礦化度場擬合圖Fig. 4 Fitting chart of groundwater mineralization

5 下遼河平原地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)系統(tǒng)

5.1 系統(tǒng)總體結(jié)構(gòu)

本系統(tǒng)開發(fā)選擇以C/S結(jié)構(gòu)為基礎(chǔ)的類三層結(jié)構(gòu)，包括客戶層、類中間層和服務(wù)層。客戶層由本次編寫的程序來完成，是用戶直接操作的程序，絕大部分用戶的操作、數(shù)據(jù)的處理、功能的實(shí)現(xiàn)等都由它來完成，是類三層結(jié)構(gòu)中功能最強(qiáng)的一部分。類中間層由Microsoft SQL Server 2000數(shù)據(jù)庫的服務(wù)器端來完成，主要用于把客戶層所請求的操作轉(zhuǎn)換為服務(wù)層所能識別的語言，并進(jìn)行數(shù)據(jù)的完整性與規(guī)范性檢查。服務(wù)層所完成的功能主要是調(diào)用實(shí)時(shí)數(shù)據(jù)等，主要由Microsoft SQL Server 2000數(shù)據(jù)庫管理系統(tǒng)來實(shí)現(xiàn)[14]。

本系統(tǒng)包括地下水流實(shí)時(shí)預(yù)報(bào)模型和地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)模型，它們分別以MF2K和MT3DMS為計(jì)算內(nèi)核。在開發(fā)各模型的過程中各自獨(dú)立進(jìn)行，模型與模型之間通過鏈接文件實(shí)現(xiàn)對接。各模型通過連接實(shí)時(shí)監(jiān)控系統(tǒng)的數(shù)據(jù)庫獲取實(shí)時(shí)數(shù)據(jù)，也可不與數(shù)據(jù)庫相連，采用人工輸入的方式進(jìn)行操作。系統(tǒng)總體結(jié)構(gòu)見圖5。

圖5 下遼河平原地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)系統(tǒng)結(jié)構(gòu)圖Fig.5 Structure of real time groundwater quality prediction system in Lower Liaohe River Plain

5.2 系統(tǒng)開發(fā)模式

根據(jù)系統(tǒng)的總體結(jié)構(gòu)及功能設(shè)計(jì)，本次研究選擇模塊化的開發(fā)方式，即在系統(tǒng)總體結(jié)構(gòu)的框架下，單獨(dú)開發(fā)各特定功能的模塊化程序，以供其他模塊調(diào)用。每個(gè)模塊具有相對的獨(dú)立性，可完成各自模塊相應(yīng)的設(shè)計(jì)功能。采用模塊化封裝后的程序出錯(cuò)概率很小，提升了程序的穩(wěn)定性。

圖6 系統(tǒng)快速實(shí)時(shí)預(yù)報(bào)示意圖Fig. 6 Diagram of real time prediction of the system

本次系統(tǒng)的集成開發(fā)以及模型與地下水實(shí)時(shí)監(jiān)控系統(tǒng)對接的接口開發(fā)采用Visual Basic語言，對MF2K和MT3DMS內(nèi)核計(jì)算程序的修改及編譯采用Fortran語言。

5.3 系統(tǒng)整合及功能實(shí)現(xiàn)

在下遼河平原地下水流數(shù)值模型及地下水溶質(zhì)運(yùn)移模型的基礎(chǔ)上，以MF2K和MT3DMS為內(nèi)核計(jì)算程序，采用VB語言對其進(jìn)行了系統(tǒng)集成，并開發(fā)了其與地下水實(shí)時(shí)監(jiān)測數(shù)據(jù)庫的接口，建立了下遼河平原地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)系統(tǒng)，系統(tǒng)的主要功能如下：

1)地下水監(jiān)測信息的實(shí)時(shí)接收處理功能。能夠?qū)崟r(shí)接收通過自動監(jiān)測儀器獲得的監(jiān)測數(shù)據(jù)，如實(shí)時(shí)地下水位、實(shí)時(shí)地下水離子濃度等；聯(lián)機(jī)程序化輸入通過手持式移動數(shù)據(jù)采集器獲得的相關(guān)測站(點(diǎn))的監(jiān)測數(shù)據(jù)；交互式手工錄入通過人工觀測記錄獲得的監(jiān)測數(shù)據(jù)資料；程序化批量輸入歷史監(jiān)測(觀測)數(shù)據(jù)資料；聯(lián)系輸入通過網(wǎng)絡(luò)傳輸獲得的監(jiān)測信息。

2)地下水流場實(shí)時(shí)預(yù)報(bào)功能。根據(jù)實(shí)時(shí)監(jiān)測的水文氣象、地下水水位等信息，或通過手動輸入賦值，以及用戶設(shè)置的應(yīng)力期和時(shí)間步長，實(shí)時(shí)計(jì)算模擬期各源匯項(xiàng)，從而對研究區(qū)地下水流場進(jìn)行預(yù)報(bào)，為地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)提供地下水動力場數(shù)據(jù)支持。

3)地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)功能。根據(jù)實(shí)時(shí)監(jiān)測的地下水水質(zhì)數(shù)據(jù)，以及實(shí)時(shí)預(yù)報(bào)的地下水流場數(shù)據(jù)，對全區(qū)地下水水質(zhì)濃度場進(jìn)行實(shí)時(shí)預(yù)報(bào)。

4)實(shí)時(shí)模擬預(yù)報(bào)結(jié)果可視化功能。利用VB在系統(tǒng)內(nèi)部窗口中調(diào)用Surfer Automation功能，將實(shí)時(shí)預(yù)報(bào)的地下水流場和地下水水質(zhì)濃度場進(jìn)行圖形可視化。

5)模型的修改完善功能。隨著實(shí)際監(jiān)測值的不斷獲取，可利用模型計(jì)算值與監(jiān)測值進(jìn)一步進(jìn)行驗(yàn)證，不斷修改完善本系統(tǒng)。系統(tǒng)預(yù)留更新升級接口，只需替換程序中水量水質(zhì)模型的相關(guān)源文件即可完成。

5.4 系統(tǒng)應(yīng)用——快速預(yù)報(bào)實(shí)例

根據(jù)系統(tǒng)預(yù)設(shè)方案進(jìn)行下遼河平原地下水水質(zhì)實(shí)時(shí)快速預(yù)報(bào):

1)選擇平水年及模擬期為5 a的快速預(yù)報(bào)方案。

2)從地下水實(shí)時(shí)監(jiān)測系統(tǒng)數(shù)據(jù)庫中讀取當(dāng)前時(shí)間的地下水水位和水化學(xué)成分濃度值，并由模型在后臺自動插值生成初始流場和初始濃度場。

3)選擇預(yù)設(shè)源匯項(xiàng)為平水年的默認(rèn)設(shè)置，直接運(yùn)行水流模型，得到水質(zhì)模型所需的地下水流場。

4)在水流模型運(yùn)行完畢后，運(yùn)行水質(zhì)模型，得到預(yù)報(bào)結(jié)果，見圖6。

6 結(jié)語

采用VB語言和數(shù)據(jù)庫技術(shù)，基于MF2K和MT3DMS建立的下遼河平原地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)系統(tǒng)能夠?qū)崟r(shí)獲取地下水水位、水質(zhì)監(jiān)測數(shù)據(jù)，對地下水流場和地下水溶質(zhì)濃度場的演化進(jìn)行快速預(yù)報(bào)，實(shí)現(xiàn)了地下水實(shí)時(shí)監(jiān)控技術(shù)與地下水?dāng)?shù)值模擬技術(shù)的結(jié)合。與傳統(tǒng)的地下水?dāng)?shù)值模擬相比，使用該系統(tǒng)無需耗費(fèi)大量時(shí)間和精力進(jìn)行數(shù)據(jù)處理和錄入工作，減少了出錯(cuò)概率，系統(tǒng)操作簡便，時(shí)效性強(qiáng)，并可根據(jù)逐步積累的監(jiān)測數(shù)據(jù)和計(jì)算數(shù)據(jù)進(jìn)行進(jìn)一步驗(yàn)證，不斷提高預(yù)報(bào)精度。該系統(tǒng)可為下遼河平原地下水資源的合理開發(fā)利用及污染預(yù)警提供科學(xué)依據(jù)。在本系統(tǒng)的基礎(chǔ)上，建議進(jìn)一步開發(fā)集成水文地球化學(xué)運(yùn)移模型的實(shí)時(shí)預(yù)報(bào)系統(tǒng)，不斷完善系統(tǒng)功能，提高模型預(yù)報(bào)精度。

[1] 莫鳳珍, 潘明杰. 遼寧中部城市群水資源問題與對策[J]. 遼寧經(jīng)濟(jì), 2001(2): 15-16. Mo Fengzhen, Pan Mingjie. Problems and Measures of Water Resources in Liaoning Central Cities Group[J]. Liaoning Economy, 2001(2):15-16.

[2] 孫永平. 遼寧中部城市群功能關(guān)系優(yōu)化研究[D]. 長春: 東北師范大學(xué), 2006. Sun Yongping. Study on Optimization for Cities Group Function-Relation in Central Part of Liaoning Province[D]. Changchun: Northeast Normal University, 2006.

[3] 張遠(yuǎn), 鄭丙輝, 王西琴, 等. 遼河流域渾河、太子河生態(tài)需水量研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2007, 27(6): 937-943. Zhang Yuan, Zheng Binghui, Wang Xiqin, et al. Study of Ecological Instream of the Hun River and Taizi River in Liao River Basin[J]. Acta Scientiae Circumstance, 2007, 27(6): 937-943.

[4] 林嵐, 遲寶明, 戴長雷. 渾河太子河流域1956--2000年地下水動態(tài)特征分析[J]. 吉林水利, 2004 (10): 1-4. Lin Lan, Chi Baoming, Dai Changlei. Dynamic Feature Analysis for the Groundwater in Hunhetaizi River Basin from 1956 to 2000[J]. Jilin Water Resources, 2004 (10): 1-4.

[5] 單麗, 王宇麗. 遼寧省地下水水質(zhì)現(xiàn)狀變化趨勢[J]. 東北水利水電, 2008, 26(3): 68-70. Shan Li, Wang Yuli. Trend of Groundwater Change in Liaoning Province[J]. Water Resources & Hydropower of Northeast China, 2008, 26(3): 68-70.

[6] 楊小南, 李宇斌, 胡成, 等. 渾河、太子河污染物監(jiān)控?cái)嗝鎯?yōu)化分析[C]//2007中國環(huán)境科學(xué)學(xué)會學(xué)術(shù)年會優(yōu)秀論文集. 北京: 中國環(huán)境科學(xué)學(xué)會, 2007: 1561-1568. Yang Xiaonan, Li Yubin, Hu Cheng, et al. Optimiazation and Analysis of Monitoring Sections of Contaminants in Hun River and Taizi River[C]//Proceedings of 2007 Annual Scientific Conference of China Environmental Science Committee. Beijing: Chinese Society for Environmental Sciences, 2007: 1561-1568.

[7] 杜超. 雙城市地下水資源評價(jià)及可持續(xù)利用研究[D]. 長春: 吉林大學(xué), 2008. Du Chao. Research on the Assessment and Sustainable Utilization of Groundwater Resources in Shuangcheng City[D]. Changchun: Jilin University, 2008.

[8] 杜超, 肖長來, 劉秀亮, 等. 基于Bossel框架的雙城市地下水資源可持續(xù)利用評價(jià)[J]. 吉林大學(xué)學(xué)報(bào): 地球科學(xué)版, 2010, 40(2): 331-336. Du Chao, Xiao Changlai, Liu Xiuliang, et al. Sustainable Utilization Assessment of Groundwater Resources in Shuangcheng City Based on Bossel Frame[J]. Journal of Jilin University: Earth Science Edition, 2010, 40(2): 331-336.

[9] 劉昌軍, 趙華, 張順福, 等. 臺蘭河地下水庫輻射井抽水過程的非穩(wěn)定滲流場的有限元分析[J]. 吉林大學(xué)學(xué)報(bào): 地球科學(xué)版, 2013, 43(3): 922-930. Liu Changjun, Zhao Hua, Zhang Shunfu, et al. Finite Element Analysis on Unsteady Seepage Field of Groundwater Reservoir of Tailan River During the Pumping Water of the Radiation Well[J]. Journal of Jilin University: Earth Science Edition, 2013, 43(3): 922-930.

[10] Harbaugh A W, Banta E R, Hill M C, et al. MODFLOW-2000, the US Geologcial Survey Modular Groundwater Model-User Guide to Modularization Concepts and the Groundwater Flow Process[R]. London: US Geologcial Survey, 2000.

[11] Zheng Chunmiao. MT3DMS V5.2 Supplemental User’s Guide[R]. Tuscaloosa: University of Alabama, 2006.

[12] 趙勇勝. 彌散度及其在地下水污染模型中的作用[J]. 長春地質(zhì)學(xué)院學(xué)報(bào), 1992, 22(2): 208-211. Zhao Yongsheng. Dispersivity and Its Function on Groundwater Pollution Model[J]. Journal of Changchun University of Earth Sciences, 1992, 22(2): 208-211.

[13] 林麗蓉, 唐仲華. 地下水及溶質(zhì)運(yùn)移數(shù)值模擬系統(tǒng)[J]. 地質(zhì)科技情報(bào), 2003, 22(2): 103-106. Lin Lirong, Tang Zhonghua. Numerical Simulation System of Groundwater and Solute Migration[J]. Geological Science and Technology Information, 2003, 22(2): 103-106.

[14] 謝新民, 蔣云鐘, 閆繼軍, 等. 流域水資源實(shí)時(shí)監(jiān)控管理系統(tǒng)研究[J]. 水科學(xué)進(jìn)展, 2003, 14(3): 255-259. Xie Xinmin, Jiang Yunzhong, Yan Jijun, et al. Study on Real Time Monitoring and Management System for Water Resources in River Basin[J]. Advances in Water Science, 2003, 14(3): 255-259.

Development and Application of Real Time Groundwater Quality Prediction System: An Example in Lower Liaohe River Plain

Du Chao1, Xiao Changlai2, Lü Jun1, Tian Haoran1

1.SongliaoInstituteofWaterEnvironmentScience,SongliaoRiverBasinWaterResourcesProtectionBureau,Changchun130021,China2.CollegeofEnvironmentandResources,JilinUniversity,Changchun130021,China

Groundwater numerical simulation requires mass monitoring data as basic support. The input of data requires much time and energy, which may cause fault and could not predict groundwater quality immediately. Using VB language, based on calibrated and validated numerical model of groundwater flow and solute transport of Lower Liaohe River Plain, re-compiled MF2K and MT3DMS as calculation kernel programs, Microsoft SQL Server 2000 as database, the interface of real time groundwater monitoring system is established, integration of Lower Liaohe River Plain groundwater quality real-time forecast model is established, and the operation system is developed, achieving the function of real time groundwater quality prediction in Lower Liaohe River Plain. The system combines real-time monitoring technology with the groundwater numerical simulation technology. It has advantages of simple operation, strong timeliness and small error probability. The system can collect real time data of groundwater level and quality, and predict future groundwater level and groundwater solute concentration immediately. Meanwhile, update interface is pre-set to improve prediction accuracy according to accumulated groundwater monitoring data and calculated data.

real time prediction system; numerical simulation; groundwater quality; solute transport; system integration

10.13278/j.cnki.jjuese.201405206.

2013-12-19

國家水體污染控制與治理科技重大專項(xiàng)(2012ZX07601002-002)

杜超(1985--)，男，工程師, 博士，主要從事水資源保護(hù)方面的研究和工作，E-mail:cattle53@163.com。

10.13278/j.cnki.jjuese.201405206

P641

A

杜超, 肖長來, 呂軍，等.地下水水質(zhì)實(shí)時(shí)預(yù)報(bào)系統(tǒng)開發(fā)及應(yīng)用：以下遼河平原為例.吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2014,44(5):1625-1632.

Du Chao, Xiao Changlai, Lü Jun, et al.Development and Application of Real Time Groundwater Quality Prediction System: An Example in Lower Liaohe River Plain.Journal of Jilin University:Earth Science Edition,2014,44(5):1625-1632.doi:10.13278/j.cnki.jjuese.201405206.