国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

沖蝕型大陸島地下水?dāng)?shù)值模擬研究——以東海島為例

2015-10-23 02:47:20滕建標(biāo)劉蘊(yùn)芳周雯李國(guó)敏盧耀東
關(guān)鍵詞:湛江中層水文地質(zhì)

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

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

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

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

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

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

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

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

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

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

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

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

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

圖2 2008年?yáng)|海島地下水流場(chǎng)實(shí)測(cè)圖Fig.2 Groundwater contour of Donghai Island in 2008

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

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

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

通過(guò)分析評(píng)價(jià)區(qū)實(shí)際水文地質(zhì)條件,建立如下水文地質(zhì)概念模型:非均質(zhì)各向異性;上邊界為降水補(bǔ)給、蒸發(fā)和井排泄邊界;下邊界為隔水邊界;模型第一層的四周邊界以海岸帶為邊界,概化為定水頭邊界;模型其他層的四周邊界由海岸帶向四周延伸2km,并將北部和南部邊界概化為流量邊界,通過(guò)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é)模型,利用國(guó)際通用的visuamodflow軟件求解。模型計(jì)算單元為100m×100m矩形網(wǎng)格,垂向上按實(shí)際鉆孔資料進(jìn)行分層。最終剖分共計(jì)500列、392行、7層,共計(jì)592557個(gè)有效計(jì)算單元,如圖3所示。

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

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

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

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

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

圖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í)測(cè)地下水水位動(dòng)態(tài)觀測(cè)資料對(duì)數(shù)學(xué)模型進(jìn)行了檢驗(yàn)。采用模型識(shí)別后的參數(shù)來(lái)計(jì)算各長(zhǎng)觀測(cè)孔水位動(dòng)態(tài),模擬計(jì)算得到的水位與實(shí)測(cè)水位比較結(jié)果如圖7、圖8、圖9所示。由圖可知,實(shí)測(cè)值與計(jì)算值擬合程度較高,差別不大,基本達(dá)到模型精度要求,表明模型正確,可以進(jìn)行預(yù)測(cè)。

圖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月淺層、中層和深層流場(chǎng)模擬圖,對(duì)比東海島的等水位線圖(圖3)可知,模擬結(jié)果與實(shí)際流場(chǎng)基本一致,模型可用來(lái)進(jìn)行后續(xù)的水量水質(zhì)預(yù)報(bào)。

圖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)建;②針對(duì)沖蝕型大陸島不同含水層與大陸之間的水文地質(zhì)關(guān)系而采取的模型概化措施;③地下水水流三維數(shù)值模擬檢驗(yàn)。模型的識(shí)別驗(yàn)證結(jié)果表明,所建數(shù)學(xué)模型能夠較好地反映東海島水文地質(zhì)條件,可用做東海島未來(lái)的水質(zhì)水位預(yù)報(bào)和地下水資源管理。

[1]劉錫清.關(guān)于海洋島嶼的成因類型問(wèn)題[J].海洋地質(zhì)動(dòng)態(tài),2000,16(8):1-2.

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

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

[4]成建梅,李國(guó)敏,陳崇希.濱海、海島海水入侵?jǐn)?shù)值模擬研究——以山東煙臺(tái)市和廣西潿洲島為例[M].武漢:中國(guó)地質(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]薛禹群.中國(guó)地下水?dāng)?shù)值模擬的現(xiàn)狀與展望[J].高校地質(zhì)學(xué)報(bào),2010,16(1):1-6.

[7]廣東省地質(zhì)勘查局水文工程地質(zhì)一大隊(duì).廣東省湛江市區(qū)2007年地下水動(dòng)態(tài)監(jiān)測(cè)年鑒[R].湛江:廣東省地質(zhì)勘查局水文工程地質(zhì)一大隊(duì),2008.

[8]盧耀東.湛江經(jīng)濟(jì)技術(shù)開(kāi)發(fā)區(qū)東海島新區(qū)水文地質(zhì)特征報(bào)告[R].湛江:廣東省地質(zhì)局水文工程地質(zhì)一大隊(duì),2007.

猜你喜歡
湛江中層水文地質(zhì)
湛江美術(shù)
湛江鼓舞角逐“山花獎(jiǎng)”
漫步湛江
黃河之聲(2022年12期)2022-09-27 14:48:00
基于抽水試驗(yàn)計(jì)算水文地質(zhì)參數(shù)的分析與確定
基于GPRS實(shí)現(xiàn)自動(dòng)化水文地質(zhì)抽水試驗(yàn)
基于轉(zhuǎn)移概率的三維水文地質(zhì)結(jié)構(gòu)劃分
雙麥克斯韋分布下極區(qū)中層塵埃粒子帶電研究
水文地質(zhì)在工程地質(zhì)勘察中的有效應(yīng)用
寫(xiě)湛江、畫(huà)湛江大型書(shū)畫(huà)基層采風(fēng)作品選
企業(yè)中層管理團(tuán)隊(duì)在運(yùn)營(yíng)中的困境與作用
牡丹江市| 永康市| 启东市| 石嘴山市| 财经| 吉隆县| 扶风县| 沙坪坝区| 庆安县| 神池县| 长宁县| 深州市| 乌苏市| 柘荣县| 遂宁市| 林芝县| 时尚| 兴城市| 武鸣县| 武功县| 吉木乃县| 视频| 宁强县| 道孚县| 宣威市| 新邵县| 枝江市| 荥阳市| 黑山县| 涞源县| 兴业县| 长泰县| 休宁县| 伽师县| 汝阳县| 牡丹江市| 韩城市| 屏南县| 新安县| 隆子县| 荔波县|