滕建標(biāo)，劉蘊(yùn)芳，周雯，李國敏，盧耀東

(1.環(huán)境保護(hù)部華南環(huán)境科學(xué)研究所，廣州 510655;2.中國科學(xué)院地質(zhì)與地球物理研究所，北京 100029;3.廣東省水文地質(zhì)一大隊(duì)，廣東湛江 524049)

大陸島是指地貌、地質(zhì)結(jié)構(gòu)與鄰近的大陸一致或有密切聯(lián)系的島嶼，按其形成的原因可分為構(gòu)造島和沖蝕島兩種［1］。湛江東海島島內(nèi)中、深層含水層與大陸(湛江市區(qū))通過湛江灣底部相連，相應(yīng)含水層水文地質(zhì)條件一致，地下水水力聯(lián)系密切，是典型的沖蝕型大陸島。

大陸島絕大多數(shù)是基巖島，由于海水入侵［2-3］、地層尖滅等原因，一般作為獨(dú)立的水文地質(zhì)單元進(jìn)行研究。國內(nèi)外不少學(xué)者對于這些島嶼作了相關(guān)的水文地質(zhì)研究。李國敏等［4］、LiuC［5］分別對廣西潿洲島、金門島作過水文地質(zhì)數(shù)值模擬研究。但是，對于與大陸地下水聯(lián)系密切的大陸島開展數(shù)值模擬研究仍較少。

數(shù)值模擬技術(shù)能夠真實(shí)地刻畫含水層系統(tǒng)復(fù)雜的水文地質(zhì)條件，在地下水資源的研究、評價和管理工作中被廣泛采用［6］。本文以廣東湛江東海島為例，對非獨(dú)立的沖蝕型大陸島地下水進(jìn)行地下水?dāng)?shù)值模擬研究，建立東海島地下水三維數(shù)值模型，模擬東海島地下水流場，對于沖蝕型大陸島數(shù)值模擬研究具有借鑒意義。

1 東海島水文地質(zhì)概況

東海島位于廣東省湛江市南部，面積約286km2。自20世紀(jì)60年代開始開發(fā)利用島內(nèi)中層承壓水作為水源以來，隨著本島及鄰區(qū)開采量的不斷增加，中、深層承壓含水層的水位也在逐年下降。據(jù)該區(qū)地下水動態(tài)長期觀測資料顯示，2007年島內(nèi)中層承壓水降幅最大的監(jiān)測孔水位標(biāo)高已達(dá)－9.71m，水位埋深比1966年降低了14.21m;降幅最大的深層承壓水監(jiān)測孔水位標(biāo)高為－7.80m，比1966年下降了15.80 m［7］。隨著東海島地下水水資源供需矛盾日益突出，有必要對研究區(qū)地下水水資源進(jìn)行定量研究，為島內(nèi)地下水資源開發(fā)利用提供指導(dǎo)。

東海島平均海平面高程為2.2m，地貌類型以北海組平原為主，其次為濱海平原環(huán)島分布，湛江組臺地僅局部分布，玄武巖臺地分布于東部龍水嶺一帶。

東海島位于東海斷凹陷盆地內(nèi)，盆底主要由白堊系地層構(gòu)成，地表出露地層為第四系，其下沉積厚約1000m的新生代砂質(zhì)和泥質(zhì)巖類互層。島內(nèi)地下水系統(tǒng)主要由松散巖類孔隙含水巖組構(gòu)成［8］。含水層按照埋深、水力特征和開采條件可分為三層，如圖1、圖2所示。

(1)淺層水:全島各地均有分布，含水層埋深一般5～30m，巖性主要為粗砂、中砂、細(xì)砂等，厚度3～10m不等。島內(nèi)淺層水接受大氣降水的垂向入滲補(bǔ)給，一部分垂向越流補(bǔ)給下伏含水層，一部分沿著其弱透水層底板潛流入海。

(2)中層承壓水:分布于全島各地，巖性為湛江組的粗砂、礫砂及中細(xì)砂，總厚度50～90m。中層水與淺層水之間以10～25m厚的湛江組黏土相隔。中層承壓水主要接受淺層水的越流補(bǔ)給和南側(cè)雷州灣的側(cè)向補(bǔ)給，受到湛江市區(qū)大量開采形成降落漏斗的影響，地下水總體上由南向北徑流。

(3)深層承壓水:分布于全島各地，巖性為下洋組的礫石、礫砂及中粗砂。深層水和中層水之間以下洋組黏土層相隔。含水層頂板埋深205～332m。地下水位標(biāo)高為－2～－8m。該層為島內(nèi)供水的主要開采層位之一，富水性較好，接受上層含水層的越流補(bǔ)給和西南側(cè)側(cè)向來水補(bǔ)給，向湛江市區(qū)漏斗區(qū)排泄。

圖1 東海島水文地質(zhì)剖面示意圖Fig.1 Hydrogeologic profile of Donghai Island

圖2 2008年東海島地下水流場實(shí)測圖Fig.2 Groundwater contour of Donghai Island in 2008

2 地下水?dāng)?shù)值模型建立

2.1 水文地質(zhì)條件概化

根據(jù)東海島水文地質(zhì)條件，島內(nèi)中砂、粗砂和礫砂等砂性土富水性較好，賦存有豐富的地下水，為區(qū)內(nèi)主要含水層;黏土、砂質(zhì)黏土等黏性土富水性和透水性均較差，為相對隔水層。

通過分析評價區(qū)實(shí)際水文地質(zhì)條件，建立如下水文地質(zhì)概念模型:非均質(zhì)各向異性;上邊界為降水補(bǔ)給、蒸發(fā)和井排泄邊界;下邊界為隔水邊界;模型第一層的四周邊界以海岸帶為邊界，概化為定水頭邊界;模型其他層的四周邊界由海岸帶向四周延伸2km，并將北部和南部邊界概化為流量邊界，通過modflow的井文件實(shí)現(xiàn)。中層水概化為2層含水層，淺、深層含水層分別概化為1層，即4層含水層。

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

上述水文地質(zhì)概念模型可用如下數(shù)學(xué)模型描述:

式中，Ω為地下水滲流區(qū)域;S1為模型的第一類邊界;S2為模型的第二類邊界;kxx、kyy、kzz分別為x、y、z主方向的滲透系數(shù)，m/s;w為源匯項(xiàng)，包括降水入滲補(bǔ)給、蒸發(fā)、井的抽水量和泉的排泄量，m3/s;μs為貯水率，1/m;H0(x，y，z)為初始地下水水頭函數(shù)，m;H1(x，y，z)為第一類邊界地下水水頭函數(shù)，m;q(x，y，z，t)為第二類邊界單位面積流量函數(shù)，m3/s。

2.3 模型構(gòu)建

根據(jù)水文地質(zhì)概念模型和數(shù)學(xué)模型，利用國際通用的visuamodflow軟件求解。模型計(jì)算單元為100m×100m矩形網(wǎng)格，垂向上按實(shí)際鉆孔資料進(jìn)行分層。最終剖分共計(jì)500列、392行、7層，共計(jì)592557個有效計(jì)算單元，如圖3所示。

為了較準(zhǔn)確地刻畫評價區(qū)水文地質(zhì)條件，本模型依據(jù)水文地質(zhì)圖及抽水試驗(yàn)結(jié)果，將評價區(qū)分成了14個參數(shù)分區(qū)，同時通過計(jì)算水位和實(shí)際水位擬合分析，反復(fù)調(diào)參，得到最終的含水層參數(shù)。

圖3 模型三維視圖(垂向放大30倍)Fig.3 Model of 3D view(Verticalmagnified30times)

3 模型識別與檢驗(yàn)

3.1 地下水水位動態(tài)擬合

根據(jù)水文地質(zhì)模型所建立的數(shù)學(xué)模型，必須反映實(shí)際流場的特點(diǎn)，因此，必須對數(shù)學(xué)模型進(jìn)行校正，即校正其參數(shù)以及邊界條件等是否能確切地反映計(jì)算區(qū)的實(shí)際水文地質(zhì)條件。本次模擬依據(jù)2007年的18個長觀孔地下水水位資料來對模型進(jìn)行校正，共歷時365d，源匯項(xiàng)包括降水、蒸發(fā)、人工開采等。對模型求解后得到在給定水文地質(zhì)參數(shù)和各均衡項(xiàng)條件下地下水位動態(tài)擬合曲線如圖4、圖5、圖6所示。通過擬合長觀孔水位動態(tài)，識別水文地質(zhì)參數(shù)、邊界和其他均衡項(xiàng)。由于參數(shù)分區(qū)和參數(shù)初值選取客觀地反映了研究區(qū)的水文地質(zhì)條件，經(jīng)細(xì)致調(diào)參后，模型識別取得了較好的效果。

圖4 2007年淺層地下水水位擬合曲線Fig.4 The curve fitting of shallow groundwater in 2007

3.2 模型驗(yàn)證

圖5 2007年中層地下水水位擬合曲線Fig.5 The curve fitting of middle groundwater in 2007

圖6 2007年深層地下水水位擬合曲線Fig.6 The curve fitting of deep groundwater in 2007

為進(jìn)一步驗(yàn)證所建立的數(shù)學(xué)模型和模型參數(shù)的可靠性，利用2008—2009年實(shí)測地下水水位動態(tài)觀測資料對數(shù)學(xué)模型進(jìn)行了檢驗(yàn)。采用模型識別后的參數(shù)來計(jì)算各長觀測孔水位動態(tài)，模擬計(jì)算得到的水位與實(shí)測水位比較結(jié)果如圖7、圖8、圖9所示。由圖可知，實(shí)測值與計(jì)算值擬合程度較高，差別不大，基本達(dá)到模型精度要求，表明模型正確，可以進(jìn)行預(yù)測。

圖7 2008年淺層地下水水位擬合曲線Fig.7 The curve fitting of shallow groundwater in 2008

圖8 2008年中層地下水水位擬合曲線Fig.8 The curve fitting of middle groundwater in 2008

圖9 2008年深層地下水水位擬合曲線Fig.9 The curve fitting of deep groundwater in 2008

圖10、圖11、圖12為東海島2009年6月淺層、中層和深層流場模擬圖，對比東海島的等水位線圖(圖3)可知，模擬結(jié)果與實(shí)際流場基本一致，模型可用來進(jìn)行后續(xù)的水量水質(zhì)預(yù)報。

圖10 2009年6月東海島淺層地下水等水位圖Fig.10 Calculated of shallow groundwater contour，Donghai Island in June 2009

4 結(jié)論

圖11 2009年6月東海島中層地下水等水位圖Fig.11 Calculated of middle groundwater contour，Donghai Is land in June 2009

圖12 2009年6月東海島深層地下水等水位圖Fig.12 Calculated of deep groundwater contour，Donghai Island in June 2009

在分析東海島水文地質(zhì)條件的基礎(chǔ)上，構(gòu)建了東海島地下水水流三維數(shù)值模型。本次主要工作內(nèi)容包括:①水文地質(zhì)條件分析與模型構(gòu)建;②針對沖蝕型大陸島不同含水層與大陸之間的水文地質(zhì)關(guān)系而采取的模型概化措施;③地下水水流三維數(shù)值模擬檢驗(yàn)。模型的識別驗(yàn)證結(jié)果表明，所建數(shù)學(xué)模型能夠較好地反映東海島水文地質(zhì)條件，可用做東海島未來的水質(zhì)水位預(yù)報和地下水資源管理。

［1］劉錫清.關(guān)于海洋島嶼的成因類型問題［J］.海洋地質(zhì)動態(tài)，2000，16(8):1-2.

［2］楊巧鳳，李文鵬，王瑞久.深圳沿海表層海水與地表水關(guān)系的地球化學(xué)分析［J］.水文地質(zhì)工程地質(zhì)，2010，37(1):45-49.

［3］楊巧鳳，李文鵬，王瑞文.深圳大沙河河水與海水關(guān)系的水文化學(xué)識別［J］.水文地質(zhì)工程地質(zhì)，2010，37(2):32-36.

［4］成建梅，李國敏，陳崇希.濱海、海島海水入侵?jǐn)?shù)值模擬研究——以山東煙臺市和廣西潿洲島為例［M］.武漢:中國地質(zhì)大學(xué)出版社，2004:112-117.

［5］LiuC，LinC，JangC，etal.Sustainable groundwater management in Kinmen Island［J］.Hydrological Processes，2006，20(20):4363-4372.

［6］薛禹群.中國地下水?dāng)?shù)值模擬的現(xiàn)狀與展望［J］.高校地質(zhì)學(xué)報，2010，16(1):1-6．

［7］廣東省地質(zhì)勘查局水文工程地質(zhì)一大隊(duì)．廣東省湛江市區(qū)2007年地下水動態(tài)監(jiān)測年鑒［R］．湛江:廣東省地質(zhì)勘查局水文工程地質(zhì)一大隊(duì)，2008．

［8］盧耀東．湛江經(jīng)濟(jì)技術(shù)開發(fā)區(qū)東海島新區(qū)水文地質(zhì)特征報告［R］．湛江:廣東省地質(zhì)局水文工程地質(zhì)一大隊(duì)，2007．