• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    硇洲島海水入侵?jǐn)?shù)值模擬

    2016-05-10 05:50:45段梅
    關(guān)鍵詞:數(shù)值模擬

    段梅

    (廣東省地質(zhì)局第四地質(zhì)大隊(duì),湛江 524049)

    ?

    硇洲島海水入侵?jǐn)?shù)值模擬

    段梅

    (廣東省地質(zhì)局第四地質(zhì)大隊(duì),湛江524049)

    摘要:通過對硇洲島地質(zhì)和水文地質(zhì)條件的綜合分析,建立了本研究區(qū)三維地下水?dāng)?shù)值模型,并獲取了相關(guān)含水層參數(shù)。從而建立海水入侵?jǐn)?shù)值模型,對硇洲島在不同開采條件下的海水入侵趨勢進(jìn)行了模擬預(yù)測,結(jié)果表明,減少開采量能減緩硇洲島海水入侵;在滿足當(dāng)?shù)毓┧枨蟮那疤嵯?,減少淺層水開采量,在遠(yuǎn)離沿海區(qū)域的中部地區(qū)適當(dāng)增加中層水開采量的開采方案優(yōu)于現(xiàn)狀開采條件下的方案。本研究也為減緩海水入侵提供了相關(guān)科學(xué)依據(jù)。

    關(guān)鍵詞:硇洲島;海水入侵;數(shù)值模擬;模型預(yù)測

    海水入侵是指由于濱海地區(qū)地下水動(dòng)力條件或水文地質(zhì)條件發(fā)生變化,引起海水或高礦化咸水向陸地淡水含水層運(yùn)移而發(fā)生的水體侵入的過程和現(xiàn)象。海水入侵早在19世紀(jì)就引起了人們的注意,Badon Ghyben和Herzberg分別于1889年和1901年獨(dú)立的提出了著名的計(jì)算咸淡水界面的Ghyben-Herzbeg公式,開創(chuàng)了海水入侵的研究序幕[1]。100多年來,海水入侵經(jīng)歷了從理論假設(shè)到合理概化,從理論模型、室內(nèi)實(shí)驗(yàn)?zāi)P偷綌?shù)值模型等多個(gè)階段的發(fā)展過程。Glover[2]、Cooper[3]、Henry[4]等研究了荷蘭、以色列、美國等地的咸淡水混合帶或過渡帶,用可混溶流體的對流-彌散方程代替不可混溶的鋒面表達(dá)式,沿這一方向經(jīng)過多年研究形成了對流-彌散模型。Gupta 和 Yapa等[5]運(yùn)用對流-彌散模型研究泰國曼谷附近的海水入侵。國內(nèi)關(guān)于海水入侵的研究起步較晚,李國敏等[6]利用人工彌散加權(quán)方法建立三維有限元模型,對廣西北海潿洲島的海水入侵進(jìn)行了研究。盧薇等[7]運(yùn)用FEFLOW軟件,利用數(shù)值法對珠江口東岸地區(qū)的海水入侵進(jìn)行了數(shù)值模擬研究。梁池生等[8]采用數(shù)值模擬方法建立的管理模型對雷州半島防海水入侵進(jìn)行了定量優(yōu)化研究,獲得了不同區(qū)域在給定的地下水水位降深約束下的允許開采量,為雷州半島地區(qū)防海水方案制定提供了重要依據(jù)。任雅嫻[9]通過對硇洲島地下水水質(zhì)狀況分析,認(rèn)為海水入侵是造成地下水咸化的主要原因。羅強(qiáng)等[10]對湛江市硇洲島的水文地質(zhì)條件和地下水質(zhì)變咸的現(xiàn)象進(jìn)行了分析,發(fā)現(xiàn)當(dāng)?shù)睾K肭质艹块_采和區(qū)域水文地質(zhì)條件改變兩個(gè)因素的共同影響,提出了及早涵養(yǎng)和恢復(fù)淺層水,謹(jǐn)慎、合理開采中、深層承壓水和調(diào)整產(chǎn)業(yè)發(fā)展規(guī)劃等防治對策。以上對硇洲島地區(qū)海水入侵的相關(guān)研究已表明當(dāng)?shù)睾K肭脂F(xiàn)象已相當(dāng)明顯,并有加劇趨勢。為了減輕海水入侵造成的危害,本研究從硇洲島地質(zhì)、水文地質(zhì)條件和地下滲流力學(xué)的基本原理出發(fā),選取海水中最主要的穩(wěn)定常量元素氯離子為特征指標(biāo),利用數(shù)值法對硇洲島海水入侵進(jìn)行了模擬研究,預(yù)測不同開采條件下的海水入侵趨勢,為了防治海水入侵及減緩其擴(kuò)張趨勢提供依據(jù)。

    1研究區(qū)概況及水文地質(zhì)條件

    硇洲島是雷州半島近岸一個(gè)獨(dú)立火山島嶼,位于湛江市東南方向,面積約56 km2。地面標(biāo)高10~50 m,地形沿制高點(diǎn)向四周輻射緩降。島上地表水系不發(fā)育,居民生活用水及灌溉用水、工業(yè)用水均依賴地下水。淺層火山巖孔洞裂隙水和中層承壓水是島內(nèi)的主要供水水源。20世紀(jì)90年代以來,香蕉種植和高位養(yǎng)蝦成為當(dāng)?shù)氐闹еa(chǎn)業(yè)并不斷發(fā)展,導(dǎo)致地下水開采強(qiáng)度增加,地下水水位下降,地下水水質(zhì)變咸。尤其是近來來,海水入侵范圍擴(kuò)大的速度大大增加,已經(jīng)對當(dāng)?shù)鼐用竦纳a(chǎn)生活造成了比較嚴(yán)重的影響。

    硇洲島位于雷瓊斷陷自流盆地東北部,自第三紀(jì)以來沉積了一套厚愈千米的松散-半固結(jié)沉積物,由于后期火山噴發(fā)導(dǎo)致地表普遍覆蓋著一層數(shù)米~數(shù)十米的火山巖被。其中,由砂性土和孔洞裂隙發(fā)育的火山巖組成的含水層與相對隔水的粘性土層相互疊置,構(gòu)成了多層結(jié)構(gòu)含水層組。根據(jù)地下水埋藏條件、水力特征以及本島的實(shí)際情況,島內(nèi)的地下水自上而下依次可分為第一含水巖組即淺層潛水-微承壓水(埋藏深度一般≤30 m)、第二含水巖組即中層承壓水(埋藏深度一般在30~200 m)和第三含水巖組即深層承壓水(埋藏深度一般>200 m)三大含水層組,三層含水巖組之間各有一套厚20~30 m粘土層,為相對隔水層。

    第一層地下水由火山巖孔洞裂隙水及湛江組上部微承壓孔隙水共同組成,受降水補(bǔ)給,地下水流向由中心向四周輻射,最終以泉排泄或匯入大?;蛘呦挠谌斯ら_采和越流補(bǔ)給中層地下水;富水性中東部火山口富水性最好,向四周依次降低,主要受含水層厚度控制。第二層地下水賦存于湛江組的砂礫石中,富水性較豐富,地下水位埋深最深47 m,火山口附近一般大于27 m,島嶼邊緣一般小于8 m,水位隨開采量增加逐年下降。第三層地下水賦存于下洋組的礫石、砂土中,地下水由于埋深過大,供水意義極小。第二、三含水層與陸地同層含水層相連,第二、三層地下水除主要接受上層水越流補(bǔ)給外,還接受來自大陸方向側(cè)向補(bǔ)給及通過火山口垂向補(bǔ)給,但其側(cè)向補(bǔ)給微弱,故與第一層地下水具有相似的流向和水力特征。

    硇洲島為一個(gè)相對獨(dú)立的水文地質(zhì)單元,研究區(qū)源匯項(xiàng)主要為降雨補(bǔ)給及人工開采。由于地處亞熱帶氣候區(qū),蒸發(fā)強(qiáng)烈,灌溉用水一般用于蒸發(fā),灌溉回灌與蒸發(fā)項(xiàng)相抵消。本研究區(qū)第一與第二含水層之間以及第二與第三含水層之間存在弱透水層粘土,火山口也是各含水層之間的很好連接通道,使三層含水層地下水之間存在水力聯(lián)系。故可將研究區(qū)概化為三層結(jié)構(gòu)的非穩(wěn)定的地下水三維流動(dòng)系統(tǒng)。

    2海水入侵?jǐn)?shù)學(xué)模型

    本次海水入侵的數(shù)學(xué)模型包括水流模型和氯離子遷移模型兩部分。首先要求解水流數(shù)學(xué)模型,然后將流速代入氯離子運(yùn)移模型求解氯離子隨時(shí)間與空間的變化規(guī)律。

    在研究區(qū),影響地下水壓力分布的地表水主要是海水,由于海水的密度大于淡水的密度,因此,在研究海水入侵地下水時(shí),若考慮了海水的密度,則計(jì)算結(jié)果更加合理。描述海水入侵的數(shù)學(xué)模型,變密度水流方程為

    式中,n為有效孔隙度;q為源匯項(xiàng),表示單位體積多孔介質(zhì)源或者匯的流量;ej為重力方向單位矢量的第j個(gè)分量(e1=e2=0,e3=1);ρs為海水密度;ρ為地下水密度。

    描述海水入侵鹽分運(yùn)移的對流-彌散方程為

    (1)

    (2)

    地下水水動(dòng)力彌散參數(shù)是不確定性程度很強(qiáng)的參數(shù),受很多因素影響,在實(shí)際研究中,通常將水動(dòng)力彌散系數(shù)表示為

    式中,u是合速度;ux、uy和uz分別是x、y和z軸方向上u的分速度;αT和αL分別是縱向彌散度和橫向彌散度;Dx0、Dy0和Dz0分別是分子擴(kuò)散系數(shù)。

    采用有限差分法求解上述三層結(jié)構(gòu)的地下水準(zhǔn)三維非穩(wěn)定流定解問題。在數(shù)值模擬計(jì)算中,我們采用波前法求解最終需要求解的線性代數(shù)方程組。

    3地下水?dāng)?shù)值模型

    3.1網(wǎng)格剖分及模擬期確定

    根據(jù)硇洲島33個(gè)鉆孔及地質(zhì)圖、剖面圖[11]建立其概化含水層三維結(jié)構(gòu)模型。模型共分5層,共17 670個(gè)單元格,單元格最小100 m×110 m,最大200 m×200 m。研究區(qū)自1992年開始大量開采地下水。在此之前,地下水位基本處于天然的穩(wěn)定狀態(tài)。故選取模擬期為1992年1月1日至2010年1月1日。模型共18個(gè)應(yīng)力期,一個(gè)應(yīng)力期為1 a。每個(gè)應(yīng)力期劃分12個(gè)時(shí)間步長。

    3.2邊界條件

    第一含水層被海水切割,海岸線是其定水頭邊界,因此將海水作為第一含水層的水頭邊界,由于模型中以月為最小時(shí)間步長,因此在本次模擬中忽略潮汐作用的影響。第一層地下水邊界取定水頭邊界,水頭值為0 m。第二、三含水層向海域延伸,根據(jù)“等效排泄邊界”的原理,把其邊界選在距海岸線5~15 km處,為定水頭邊界[12,13]。根據(jù)2010年野外調(diào)查中研究區(qū)邊界上三口第二層地下水開采井的水位分別為-0.05 m、-0.2 m和-0.8 m,與0 m海面高程相差不過1 m。因此設(shè)第二層地下水邊界為定水頭邊界,水頭值同淺層取0 m。第三層地下水埋深一般大于200 m,供水意義小,對第一、二層地下水影響較小,故設(shè)其邊界條件與第二層地下水一致。

    3.3初始條件和參數(shù)分區(qū)的確定

    由于硇洲島地區(qū)缺乏初始水位及參數(shù)分區(qū)等資料,我們先建立該區(qū)的穩(wěn)定流模型。利用1981年廣東省地質(zhì)圖測繪大隊(duì)繪制的1∶20萬綜合水文地質(zhì)圖中硇洲島部分第二層地下水5 m等水位標(biāo)高線[14]進(jìn)行模型初步擬合,并且根據(jù)含水層的巖性、富水性的不同,分別把第一含水層和第二含水層的水文地質(zhì)參數(shù)分別分為3個(gè)參數(shù)區(qū),而第一和第二及第二和第三含水層之間的弱透水層及第三含水層不分區(qū),即各自只有1個(gè)參數(shù)區(qū)。擬合結(jié)果見圖1,水位擬合的總體趨勢較好,說明采用的參數(shù)分區(qū)較為合理。利用該識(shí)別后穩(wěn)定流模型模擬出水位結(jié)果即作為地下水流模型的初始水位。

    圖1 第二含水層水位標(biāo)高等值線實(shí)測值與計(jì)算值對比圖

    3.4源匯項(xiàng)

    研究區(qū)源匯項(xiàng)主要包括降雨入滲和人工開采。源項(xiàng)為降雨補(bǔ)給量,由降雨入滲系數(shù)和降雨量資料計(jì)算得出。匯項(xiàng)中淺層水開采主要用于生活和農(nóng)業(yè)灌溉,開采分散,開采井密度大,成面狀開采。淺層面狀反補(bǔ)給量是根據(jù)研究區(qū)的農(nóng)業(yè)灌溉用水、水產(chǎn)養(yǎng)殖和居民生活用水這三項(xiàng)相加,再除以研究區(qū)總面積得到的一個(gè)反補(bǔ)給強(qiáng)度。中層開采井抽水量利用模型擬合機(jī)井2008年抽水量對比得出。

    3.5模型擬合及結(jié)果分析

    模型的識(shí)別擬合是利用研究區(qū)內(nèi)特征點(diǎn)、線、面、體上已知的地下水動(dòng)態(tài)觀測資料與模型在該點(diǎn)上的計(jì)算值進(jìn)行對比,以此判斷數(shù)值模型對地質(zhì)體的仿真程度,通常要進(jìn)行反復(fù)的修改參數(shù)和調(diào)整某些源匯項(xiàng)才能達(dá)到較為理想的擬合結(jié)果。本次利用2010年野外調(diào)查收集的20口機(jī)井的靜水位標(biāo)高及成井時(shí)間進(jìn)行模型擬合。擬合結(jié)果見圖2。

    圖2 第二層地下水水位標(biāo)高實(shí)測值與計(jì)算值對比圖(黑色為實(shí)測值)

    我們把水位擬合的誤差和浮動(dòng)百分比計(jì)算出來,如表1所示。從表1可以看出,所有機(jī)井?dāng)M合的誤差值最大為1 m,大多數(shù)在0.05 m以下。其中14口機(jī)井其擬合浮動(dòng)百分比都在5%以下,擬合效果較好。個(gè)別雖然浮動(dòng)大于100%,但其誤差值很小,最大僅為0.62 m。通過綜合分析水位擬合曲線和各種信息可知,所建立的地下水流數(shù)值模型基本達(dá)到模型精度要求,可用于預(yù)測計(jì)算。

    表1第二層地下水水位標(biāo)高實(shí)測值與計(jì)算值對比表

    孔號成井時(shí)間實(shí)測值/m計(jì)算值/m誤差/m浮動(dòng)百分比/%z40022004-2.58-2.450.135.14z40052007-3-2.670.3311.03z40071991-5.3-5.250.050.86z40082005-0.050.070.12243.37z40112005-0.2-0.530.33166.53z40122002-0.8-0.800.38z40131999-3.72-3.690.030.91z40141974-7.16-7.140.020.3z40152003-2.5-2.450.051.92z40162005-2.75-2.720.031.1z40171985-10.9-9.919.14z40192009-10-9.870.131.27z40212001-3.5-3.60.12.99z40221968-7.55-7.60.050.7z40251969-12.5-12.010.493.95z40272005-4.4-4.350.051.05z40281995-5.5-5.410.091.55z40311989-2.5-2.560.062.48z40322009-3.18-3.210.031.06z40331991-0.57-1.190.62109.49

    4硇洲島海水中氯離子遷移數(shù)值模擬

    氯離子遷移模型是利用水流計(jì)算結(jié)果進(jìn)行溶質(zhì)運(yùn)移的耦合計(jì)算的。故其網(wǎng)格剖分、邊界條件及模擬期都等同于硇洲島地下水流動(dòng)模型中的設(shè)定。

    根據(jù)海水中氯離子的平均濃度,在模型中設(shè)邊界的氯離子的補(bǔ)給濃度為20 000 mg/l,區(qū)域背景值假定為0 mg/l[15]。研究區(qū)所收集的資料中沒有關(guān)于進(jìn)行溶質(zhì)運(yùn)移計(jì)算所需要的參數(shù)資料,參考相關(guān)資料,設(shè)定縱向彌散度為20 m,橫向/縱向彌散度比率為0.2,垂向/縱向彌散度比率為0.2[16-18]。

    由于未收集到數(shù)據(jù)暫未擬合,但在本模型參數(shù)標(biāo)準(zhǔn)下,對氯離子運(yùn)移進(jìn)行研究還是有一定的指導(dǎo)意義的。模擬結(jié)果與實(shí)際情況是基本吻合的。模擬結(jié)果表明,以氯離子濃度為200 mg/l作為判斷海水入侵的標(biāo)準(zhǔn),硇洲島第一含水層沿海1 km以內(nèi)地區(qū)幾乎全部受到了海水入侵的影響。第二含水層在西部和北部沿海也已逐步受到淺層水咸化的影響。該結(jié)果與陳苑[19]2011年對硇洲島地下水咸化與海水入侵調(diào)查得出硇洲島內(nèi)地勢較為低緩的西南、西北與東北部沿海區(qū)域已經(jīng)出現(xiàn)海水入侵,海水最大縱深滲透距離達(dá) 1 000 m 左右,入侵區(qū)內(nèi) 30~80 m 深度范圍的地下含水層水質(zhì)普遍發(fā)生不同程度咸化的結(jié)論是相符的,這表明本次建立的模擬模型是合理的。

    5海水入侵預(yù)測

    5.1現(xiàn)狀開采條件下2020年海水入侵預(yù)測

    保持開采量及邊界條件不變,降雨補(bǔ)給量取1992年~2009年平均值0.004 0 m/d,模擬時(shí)段為2010年~2020年,共10個(gè)應(yīng)力期,每個(gè)應(yīng)力期12個(gè)時(shí)間步長。模擬結(jié)果表明淺層含水層受海水入侵影響的范圍在這10 a有所擴(kuò)展,平面擴(kuò)展速率約20~40 m/a。由于第一含水層的地形優(yōu)勢,除沿海地區(qū)海水向淡水滲流,其余大部分地區(qū)地下水流向仍保持由高地勢地區(qū)向低地勢地區(qū)滲流。所以第一含水層的海水入侵速率相對來說有所降低。

    第二含水層受海水入侵影響的范圍急劇增加,幾乎占了全區(qū)面積的三分之一。2009年底第二含水層僅西部和北部沿海地區(qū)受到海水入侵的影響,而經(jīng)過10 a的發(fā)展,幾乎所有沿海地區(qū)全部發(fā)生了海水入侵。第二含水層是硇洲島地區(qū)的主要供水層位。假如維持現(xiàn)狀開采,僅僅10 a第二含水層就將面臨著十分嚴(yán)重的海水入侵這一地質(zhì)環(huán)境問題。

    5.2優(yōu)化開采條件下2020年海水入侵預(yù)測

    本次分別預(yù)測了以下兩種開采方案下,2020年氯離子的濃度分布特征:方案1,第一含水層開采量減小而第二含水層開采量不變;方案2,第一含水層開采量減小但第二含水層開采量在遠(yuǎn)離沿海區(qū)域的中部地區(qū)增加以保證總開采量穩(wěn)定。

    方案1模擬結(jié)果與維持現(xiàn)開采條件下的模擬結(jié)果進(jìn)行對比,第一含水層水位有所上升,海水入侵的范圍減??;第二含水層仍不可避免的受到了海水入侵的影響,擴(kuò)展范圍比現(xiàn)狀開采條件下略有縮小。

    方案2模擬結(jié)果表明增加第二層地下水開采量沒有導(dǎo)致第一含水層海水入侵的進(jìn)一步發(fā)展,第二含水層受海水入侵的影響范圍相對于方案1略有增加,證明了第二含水層的海水入侵程度在一段時(shí)間內(nèi)會(huì)無可避免的加重,必須采取除減小開采量以外更有效的措施來阻止海水入侵程度的發(fā)展。

    為了更好對比上訴三種開采方案,選取z4022井所處單元格在這三種開采條件下中該單元格中氯離子濃度隨時(shí)間變化的特征,來具體分析三種開采條件下對中層含水層海水入侵程度的影響。如圖3所示。

    圖3不同開采條件下z4022所處單元格氯離子濃度隨時(shí)間變化曲線

    從圖3中可以看出,雖然方案2中增加了中層含水層的開采量,但z4022所處單元格中氯離子濃度卻低于維持現(xiàn)狀開采條件下的氯離子濃度。這說明了在不減少開采量的情況下,合理的開采布局是可以在一定程度上抑制海水入侵的發(fā)展的。

    6結(jié)論

    海水入侵對硇洲島的社會(huì)和經(jīng)濟(jì)發(fā)展會(huì)帶來不利影響,因此本文根據(jù)硇洲島水文地質(zhì)的特點(diǎn)和鉆孔資料,通過水文地質(zhì)條件概化建立硇洲島的三維穩(wěn)定流模型及瞬變流模型到氯離子的遷移模型,運(yùn)用收集到的地下水位資料對模型進(jìn)行識(shí)別和校驗(yàn),運(yùn)用識(shí)別后的模型模擬預(yù)測不同開采條件下海水入侵趨勢。結(jié)果表明,減少開采量能減緩海水入侵;在滿足當(dāng)?shù)氐墓┧枨蟮那疤嵯?,減小淺層水開采量,在遠(yuǎn)離沿海區(qū)域的中部地區(qū)適當(dāng)增加中層水開采量的開采方案優(yōu)于現(xiàn)狀開采條件下的方案。

    由于硇洲島專門的海水入侵勘探工作量較少,缺乏長觀資料和邊界資料及試驗(yàn)數(shù)據(jù),現(xiàn)有的機(jī)井觀測點(diǎn)分布不均,本研究僅為初步模擬研究。

    特別感謝張莉在本文寫作過程中提出寶貴建議和給以的幫助!

    參考文獻(xiàn)

    [1]蔡祖煌,馬鳳山.海水入侵的基本理論及其在入侵發(fā)展預(yù)測中的應(yīng)用[J].中國地質(zhì)災(zāi)害與防治學(xué)報(bào),1996,7(3):1-9.

    [2]Glover R E.The pattern of fresh water flow in a coastal aquifer. J Geophys Res, 1959, 64(4):457-459.

    [3]Cooper H H, F A Kohout, H R Henry, and R E Glover.Sea water in coastal aquifers.U.S.Geol. Surv. Water Supply Paper, 1964, 1614-C:70-84.

    [4]Henry H R.Effect of dispersion on salt encroachment in coastal aquifers.U.S.Geol. Surv., Water Supply Paper. 1964, 1613-C: 70-84.

    [5]Oude Essink G H P, Boekelman R H.Problems with large-scale modeling of salt water intrusion in 3D. 14thSalt Water Intrusion Meeting, 1996.

    [6]李國敏,陳崇希,沈照理,等.潿洲島海水入侵模擬[J].水文地質(zhì)工程地質(zhì),1995,(5):1-5.

    [7]盧薇,朱照宇,劉衛(wèi)平.基于FEFLOW的海水入侵?jǐn)?shù)值模擬[J].地下水,2010,32(3):19-21.

    [8]梁池生,周訓(xùn),陳明佑.雷州半島防海水入侵最優(yōu)地下水開采方案探討[J].現(xiàn)代地質(zhì),2002,16(4):429-434.

    [9]任雅嫻.硇洲島地下水咸化成因及變化趨勢分析[J].廣東水利水電,2013,(5):23-27.

    [10]羅強(qiáng),蘇肇漢.硇洲島地下水海水入侵的特征分析及治理對策[J],地質(zhì)災(zāi)害與環(huán)境保護(hù),2007,18(2):28-32.

    [11]廣東省地質(zhì)礦產(chǎn)局.湛江綜合區(qū)域地質(zhì)調(diào)查報(bào)告(1:50000)[R].1989.

    [12]陳崇希,林敏,舒本媛,等.濱海承壓含水層地下水的等效排泄邊界——以北海禾塘水源地為例[J].水文地質(zhì)工程地質(zhì),1990,17(4):2-4.

    [13]成建梅,黃丹紅,胡進(jìn)武.海水入侵模擬理論與方法研究進(jìn)展[J].水資源保護(hù),2004,20(2):3-8.

    [14]廣東省地質(zhì)局.雷州半島區(qū)域水文地質(zhì)普查報(bào)告(1:200000)[R].1981.

    [15]付旭.泉州市典型區(qū)域海水入侵實(shí)驗(yàn)與數(shù)值模擬研究[D].華僑大學(xué),2013.

    [16]孫訥正.地下水水質(zhì)的數(shù)學(xué)模擬(四)水動(dòng)力彌散方程的數(shù)值解法[J].水文地質(zhì)工程地質(zhì),1982,4(1):49-55.

    [17]鄒立芝,潘俊,楊昌兵,等.含水層水力參數(shù)的尺度效應(yīng)研究現(xiàn)狀[J].長春地質(zhì)學(xué)院學(xué)報(bào),1994,24(1):66-69.

    [18]張宇峰,張雪英,徐炎華,等.土壤中水動(dòng)力彌散系數(shù)的研究進(jìn)展[J].環(huán)境污染治理技術(shù)與設(shè)備,2003,4(7):8-12.

    [19]陳苑.硇洲島地下水水質(zhì)咸化和海水入侵主要成因及對策[J].大眾科技,2011,(4):114-115.

    NUMERICAL SIMULATION OF SEAWATER INTRUSION IN NAOZHOU ISLAND

    DUAN Mei

    (The Fourth Geology ofGuangdong Geological Bureau ,Zhanjiang Guangdong524049,China)

    Abstract:Through the comprehensive analysis of geological and hydrogeological conditions of the NaoZhou island, we have established the three-dimensional numerical model of groundwater and obtained the aquifer parameters in the study area. The trend of seawater intrusion of the NaoZhou island under different mining conditions have been simulated and predicted through the seawater intrusion numerical model. The results showed that reducing production could slow down the seawater intrusion of the NaoZhou island; and on the premise of meet the demand of the local water supply, the solution scheme that reduces shallow water production, and increases the middle water production appropriately which is away from the central part of the coastal areas, is superior to the present situation of mining conditions. This study also provided related scientific basis for slowing down seawater intrusion.

    Key words:Naozhou Island; seawater intrusion; numerical simulation; model prediction

    作者簡介:段梅(1985-)女,工程師,環(huán)境工程專業(yè),主要從事水文、工程和環(huán)境地質(zhì)研究工作。E-mail:41548316@qq.com

    中圖分類號:P736;P641

    文獻(xiàn)標(biāo)識(shí)碼:A

    收稿日期:2015-08-20改回日期:2015-11-17

    文章編號:1006-4362(2016)01-0108-05

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    国内久久婷婷六月综合欲色啪| 国产成人影院久久av| 亚洲国产精品成人综合色| 色播亚洲综合网| av女优亚洲男人天堂| 久9热在线精品视频| 国产成人福利小说| 99国产极品粉嫩在线观看| 欧美乱色亚洲激情| 真人一进一出gif抽搐免费| 国产亚洲欧美在线一区二区| 国产真实伦视频高清在线观看 | 成人永久免费在线观看视频| 欧美高清性xxxxhd video| 久久久久国内视频| 久久精品国产自在天天线| 亚洲经典国产精华液单 | 91av网一区二区| 少妇的逼好多水| а√天堂www在线а√下载| 国产精华一区二区三区| 婷婷丁香在线五月| 少妇的逼水好多| 三级国产精品欧美在线观看| 国产精品久久久久久精品电影| 亚洲av日韩精品久久久久久密| 中国美女看黄片| bbb黄色大片| 国产色婷婷99| 精品国产亚洲在线| 男人舔女人下体高潮全视频| or卡值多少钱| 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人| 欧美日韩综合久久久久久 | 91午夜精品亚洲一区二区三区 | 中文字幕熟女人妻在线| 久久6这里有精品| 日韩有码中文字幕| av在线老鸭窝| 一级黄色大片毛片| 国内揄拍国产精品人妻在线| 日韩欧美免费精品| 三级毛片av免费| 国产精品电影一区二区三区| 亚洲成人免费电影在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产黄a三级三级三级人| 亚洲国产欧美人成| 精品久久久久久久末码| 欧美色欧美亚洲另类二区| 国产单亲对白刺激| 色视频www国产| av女优亚洲男人天堂| 午夜精品一区二区三区免费看| 成人高潮视频无遮挡免费网站| 不卡一级毛片| 国产精品免费一区二区三区在线| 亚洲精品久久国产高清桃花| 亚洲avbb在线观看| netflix在线观看网站| 网址你懂的国产日韩在线| 国内少妇人妻偷人精品xxx网站| 男女下面进入的视频免费午夜| 日韩欧美一区二区三区在线观看| 精品免费久久久久久久清纯| 成人亚洲精品av一区二区| 中文字幕久久专区| 久久精品91蜜桃| 丁香欧美五月| 男女床上黄色一级片免费看| 99riav亚洲国产免费| 99热6这里只有精品| 亚洲久久久久久中文字幕| 男女那种视频在线观看| 亚洲不卡免费看| 网址你懂的国产日韩在线| 国产高清视频在线播放一区| 久久欧美精品欧美久久欧美| 香蕉av资源在线| 国产高清激情床上av| 久久久久久久午夜电影| 99久久无色码亚洲精品果冻| 一本一本综合久久| 内地一区二区视频在线| 精品国产亚洲在线| 亚洲精品成人久久久久久| 90打野战视频偷拍视频| 国产成年人精品一区二区| 国产精品久久久久久精品电影| 少妇裸体淫交视频免费看高清| 日韩欧美精品免费久久 | 欧美日韩综合久久久久久 | 在线a可以看的网站| www.色视频.com| 午夜福利在线在线| 精品免费久久久久久久清纯| 国产人妻一区二区三区在| 国产激情偷乱视频一区二区| 午夜福利高清视频| 97超视频在线观看视频| 久久久色成人| 日韩大尺度精品在线看网址| 婷婷精品国产亚洲av| 亚洲精品在线观看二区| 国产野战对白在线观看| 国产真实伦视频高清在线观看 | 男女视频在线观看网站免费| 亚洲av成人精品一区久久| 久久草成人影院| 两人在一起打扑克的视频| 午夜激情欧美在线| 欧美成狂野欧美在线观看| 亚洲欧美清纯卡通| 国产成人影院久久av| 嫁个100分男人电影在线观看| 性插视频无遮挡在线免费观看| 嫩草影院新地址| 久久6这里有精品| 一区二区三区免费毛片| 欧美黑人巨大hd| 禁无遮挡网站| 国产又黄又爽又无遮挡在线| 国产亚洲精品久久久久久毛片| avwww免费| 精品乱码久久久久久99久播| 88av欧美| 简卡轻食公司| 免费观看精品视频网站| 欧美一区二区精品小视频在线| 五月玫瑰六月丁香| 欧美色视频一区免费| 九色国产91popny在线| 亚洲av二区三区四区| 18禁黄网站禁片免费观看直播| or卡值多少钱| 亚洲人成电影免费在线| 99精品在免费线老司机午夜| 亚洲中文字幕日韩| 亚洲熟妇中文字幕五十中出| 久9热在线精品视频| av在线老鸭窝| 日韩欧美 国产精品| 亚洲精品成人久久久久久| 欧美一区二区亚洲| 国产中年淑女户外野战色| 国产av麻豆久久久久久久| 国产乱人伦免费视频| 精品国产三级普通话版| 最近中文字幕高清免费大全6 | 色综合站精品国产| 禁无遮挡网站| 午夜久久久久精精品| 成人国产综合亚洲| 97人妻精品一区二区三区麻豆| av在线天堂中文字幕| av女优亚洲男人天堂| 国产大屁股一区二区在线视频| 国产精品电影一区二区三区| 一二三四社区在线视频社区8| 免费看美女性在线毛片视频| 日本黄色片子视频| 韩国av一区二区三区四区| 国产精品久久久久久亚洲av鲁大| 床上黄色一级片| 久久午夜亚洲精品久久| 国产欧美日韩精品一区二区| 久久精品人妻少妇| a级毛片免费高清观看在线播放| 欧美在线一区亚洲| 可以在线观看的亚洲视频| 国产 一区 欧美 日韩| 亚洲国产精品久久男人天堂| 欧美高清成人免费视频www| 91字幕亚洲| 国产私拍福利视频在线观看| 精品欧美国产一区二区三| 日日摸夜夜添夜夜添av毛片 | av在线观看视频网站免费| 一个人观看的视频www高清免费观看| 首页视频小说图片口味搜索| 成年版毛片免费区| 国产国拍精品亚洲av在线观看| 午夜精品久久久久久毛片777| 国产大屁股一区二区在线视频| 嫩草影视91久久| 九九热线精品视视频播放| 一进一出抽搐动态| 国产伦在线观看视频一区| 国产三级中文精品| 舔av片在线| 国产成人欧美在线观看| 欧美不卡视频在线免费观看| 999久久久精品免费观看国产| 亚洲内射少妇av| 久久九九热精品免费| 国产在线精品亚洲第一网站| 国产精品乱码一区二三区的特点| 99在线人妻在线中文字幕| 欧美乱妇无乱码| 亚洲国产精品久久男人天堂| 国产极品精品免费视频能看的| 深爱激情五月婷婷| 中文字幕人妻熟人妻熟丝袜美| 精品福利观看| 宅男免费午夜| 国内精品久久久久久久电影| 国产亚洲精品av在线| 好男人电影高清在线观看| 99久久精品热视频| 国产白丝娇喘喷水9色精品| 国产大屁股一区二区在线视频| 国产成人啪精品午夜网站| 一本精品99久久精品77| 亚洲欧美日韩高清在线视频| 精品国内亚洲2022精品成人| 中文字幕av在线有码专区| 国产白丝娇喘喷水9色精品| 亚洲专区中文字幕在线| 男女做爰动态图高潮gif福利片| 一二三四社区在线视频社区8| 舔av片在线| 亚洲精品久久国产高清桃花| 国产美女午夜福利| 亚洲人成网站在线播| 久久人妻av系列| 国产精品99久久久久久久久| 一区二区三区激情视频| 婷婷色综合大香蕉| 99在线人妻在线中文字幕| 亚洲欧美日韩无卡精品| 亚洲成av人片在线播放无| 欧美激情久久久久久爽电影| 老熟妇仑乱视频hdxx| 最新中文字幕久久久久| 少妇被粗大猛烈的视频| 国产蜜桃级精品一区二区三区| 少妇被粗大猛烈的视频| 久久久久久久久久黄片| 欧美性猛交╳xxx乱大交人| 亚洲国产高清在线一区二区三| 国产av在哪里看| 亚洲国产精品sss在线观看| www.999成人在线观看| 免费高清视频大片| 欧美一区二区国产精品久久精品| 美女xxoo啪啪120秒动态图 | 国产精品亚洲av一区麻豆| a级一级毛片免费在线观看| 日韩欧美 国产精品| 噜噜噜噜噜久久久久久91| aaaaa片日本免费| 亚洲国产日韩欧美精品在线观看| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久电影中文字幕| 99视频精品全部免费 在线| 十八禁网站免费在线| 国产成人aa在线观看| 夜夜躁狠狠躁天天躁| 一进一出抽搐gif免费好疼| 免费看美女性在线毛片视频| 国产精品久久久久久久久免 | 国内精品美女久久久久久| 久久久久久久午夜电影| 亚洲av电影不卡..在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲av成人av| 欧美成人免费av一区二区三区| 婷婷亚洲欧美| 久久精品影院6| 蜜桃亚洲精品一区二区三区| 在线观看免费视频日本深夜| 亚洲最大成人av| 精品日产1卡2卡| 熟女电影av网| 宅男免费午夜| 国产欧美日韩精品一区二区| 宅男免费午夜| 性欧美人与动物交配| 内射极品少妇av片p| 亚洲乱码一区二区免费版| 国产伦人伦偷精品视频| 久久人妻av系列| 91久久精品电影网| 嫩草影院入口| 欧美激情国产日韩精品一区| 热99在线观看视频| 99在线视频只有这里精品首页| 可以在线观看毛片的网站| 国产精品亚洲一级av第二区| 99久久久亚洲精品蜜臀av| 日韩免费av在线播放| 亚洲片人在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产精品日韩av在线免费观看| 我要看日韩黄色一级片| 免费看美女性在线毛片视频| 亚洲成av人片免费观看| 在线看三级毛片| 欧美一级a爱片免费观看看| 激情在线观看视频在线高清| 欧美不卡视频在线免费观看| 国产三级中文精品| 国产高清激情床上av| 男女之事视频高清在线观看| 男女之事视频高清在线观看| 国产精品美女特级片免费视频播放器| 国产 一区 欧美 日韩| 欧美区成人在线视频| 成人国产一区最新在线观看| 久久香蕉精品热| 午夜免费男女啪啪视频观看 | 亚洲片人在线观看| 免费看美女性在线毛片视频| 精品人妻熟女av久视频| 欧美一级a爱片免费观看看| 国产视频一区二区在线看| 一区福利在线观看| 女生性感内裤真人,穿戴方法视频| 成年人黄色毛片网站| 国产精品电影一区二区三区| 国产一区二区三区视频了| 美女免费视频网站| 成人午夜高清在线视频| 亚洲美女搞黄在线观看 | 婷婷精品国产亚洲av| 老女人水多毛片| 亚洲国产精品合色在线| 国产精品乱码一区二三区的特点| 久久欧美精品欧美久久欧美| 美女黄网站色视频| 久久久久久久精品吃奶| aaaaa片日本免费| 脱女人内裤的视频| 国产午夜福利久久久久久| 国产精品av视频在线免费观看| а√天堂www在线а√下载| 两人在一起打扑克的视频| 极品教师在线视频| 国产亚洲精品综合一区在线观看| 999久久久精品免费观看国产| 久久婷婷人人爽人人干人人爱| 国产高清有码在线观看视频| 精品熟女少妇八av免费久了| 精品国内亚洲2022精品成人| 国产乱人视频| 毛片女人毛片| 国产伦一二天堂av在线观看| 国产精品,欧美在线| 免费在线观看亚洲国产| 国产亚洲精品久久久久久毛片| 综合色av麻豆| 国产 一区 欧美 日韩| 精品久久久久久,| 一区二区三区四区激情视频 | 欧美三级亚洲精品| 91在线精品国自产拍蜜月| 精品久久久久久久久久免费视频| 国产精品av视频在线免费观看| 亚洲av电影不卡..在线观看| 亚洲av成人不卡在线观看播放网| 午夜福利视频1000在线观看| 99国产精品一区二区三区| 久久久久久久久大av| 一本一本综合久久| 欧美+日韩+精品| 国产一区二区在线av高清观看| www.熟女人妻精品国产| 亚洲片人在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线美女| 毛片一级片免费看久久久久 | 脱女人内裤的视频| 婷婷精品国产亚洲av在线| 51国产日韩欧美| 国产精品1区2区在线观看.| 色综合站精品国产| 夜夜躁狠狠躁天天躁| 长腿黑丝高跟| 天堂影院成人在线观看| av国产免费在线观看| 国产真实乱freesex| 成人国产综合亚洲| 超碰av人人做人人爽久久| 日日夜夜操网爽| 久久婷婷人人爽人人干人人爱| 97超级碰碰碰精品色视频在线观看| 一区二区三区激情视频| 简卡轻食公司| 欧美日本亚洲视频在线播放| 欧美另类亚洲清纯唯美| 亚洲综合色惰| 精品一区二区三区av网在线观看| 精品久久久久久久久av| 免费一级毛片在线播放高清视频| 90打野战视频偷拍视频| 国产探花极品一区二区| 国产精品久久视频播放| 国产一区二区三区视频了| 简卡轻食公司| 成人特级黄色片久久久久久久| 女人十人毛片免费观看3o分钟| 天堂av国产一区二区熟女人妻| 国产免费av片在线观看野外av| 99热这里只有是精品在线观看 | 亚洲av一区综合| 久久人妻av系列| 中文字幕人成人乱码亚洲影| x7x7x7水蜜桃| 国产一级毛片七仙女欲春2| 免费观看人在逋| 99久久成人亚洲精品观看| 两性午夜刺激爽爽歪歪视频在线观看| or卡值多少钱| 脱女人内裤的视频| 亚洲国产精品999在线| 国产成人aa在线观看| 亚洲成人免费电影在线观看| 最近中文字幕高清免费大全6 | 国产野战对白在线观看| 成年人黄色毛片网站| 精品福利观看| 亚洲aⅴ乱码一区二区在线播放| 九色国产91popny在线| 久久人人精品亚洲av| 好男人在线观看高清免费视频| 精品福利观看| 国产伦人伦偷精品视频| 中文字幕av在线有码专区| 日本精品一区二区三区蜜桃| 欧美高清性xxxxhd video| 日韩欧美 国产精品| 中亚洲国语对白在线视频| 亚洲,欧美,日韩| 欧美一区二区国产精品久久精品| 久久国产精品人妻蜜桃| av福利片在线观看| 国产91精品成人一区二区三区| 麻豆国产av国片精品| 一本综合久久免费| 国产精品av视频在线免费观看| 91在线精品国自产拍蜜月| 级片在线观看| 成人特级黄色片久久久久久久| 热99在线观看视频| 久久久国产成人免费| 免费在线观看日本一区| 久久欧美精品欧美久久欧美| 搡老妇女老女人老熟妇| 亚洲av中文字字幕乱码综合| 日韩欧美三级三区| 99久久精品一区二区三区| 在线免费观看的www视频| 在线看三级毛片| 搡老岳熟女国产| 国产精品久久久久久久久免 | 超碰av人人做人人爽久久| 免费人成在线观看视频色| 午夜a级毛片| 1024手机看黄色片| 亚洲va日本ⅴa欧美va伊人久久| 最新中文字幕久久久久| 国产精华一区二区三区| 成人精品一区二区免费| 亚洲乱码一区二区免费版| 97热精品久久久久久| 国产伦精品一区二区三区四那| 五月伊人婷婷丁香| 国产欧美日韩一区二区精品| 亚洲一区二区三区不卡视频| 国产精品av视频在线免费观看| bbb黄色大片| 三级男女做爰猛烈吃奶摸视频| 亚洲熟妇熟女久久| 婷婷丁香在线五月| 两个人的视频大全免费| 制服丝袜大香蕉在线| 老女人水多毛片| 婷婷精品国产亚洲av在线| 国产黄色小视频在线观看| 乱码一卡2卡4卡精品| 又粗又爽又猛毛片免费看| 69av精品久久久久久| 日本免费a在线| 极品教师在线视频| 久久香蕉精品热| 国产高清视频在线观看网站| 亚洲成a人片在线一区二区| 白带黄色成豆腐渣| 日本精品一区二区三区蜜桃| 99热6这里只有精品| 白带黄色成豆腐渣| 精品久久国产蜜桃| 免费在线观看亚洲国产| 波野结衣二区三区在线| 亚洲精品粉嫩美女一区| 欧美3d第一页| 小说图片视频综合网站| 亚洲精品粉嫩美女一区| 两人在一起打扑克的视频| 国产精品亚洲一级av第二区| 亚洲18禁久久av| 婷婷丁香在线五月| 12—13女人毛片做爰片一| av在线天堂中文字幕| 十八禁国产超污无遮挡网站| 国产精品久久久久久久电影| 欧美激情在线99| 天天一区二区日本电影三级| 欧美性猛交╳xxx乱大交人| 小蜜桃在线观看免费完整版高清| 99久久精品热视频| 亚洲最大成人手机在线| 两人在一起打扑克的视频| 丁香六月欧美| 精华霜和精华液先用哪个| a级毛片免费高清观看在线播放| 久久久久久久亚洲中文字幕 | 好看av亚洲va欧美ⅴa在| 中文字幕人成人乱码亚洲影| 国产人妻一区二区三区在| 亚洲精品一卡2卡三卡4卡5卡| 麻豆国产97在线/欧美| 亚洲精品乱码久久久v下载方式| 欧美又色又爽又黄视频| 嫁个100分男人电影在线观看| 色5月婷婷丁香| 搡老岳熟女国产| 成年人黄色毛片网站| 亚洲av五月六月丁香网| 久久人人爽人人爽人人片va | 99久久无色码亚洲精品果冻| 中文资源天堂在线| 国产成年人精品一区二区| 久久精品国产亚洲av天美| 91午夜精品亚洲一区二区三区 | 日韩欧美国产一区二区入口| 亚洲五月天丁香| 精品人妻熟女av久视频| 亚洲av一区综合| 女生性感内裤真人,穿戴方法视频| 成人永久免费在线观看视频| 欧美乱妇无乱码| 美女免费视频网站| 窝窝影院91人妻| 一个人观看的视频www高清免费观看| 午夜福利高清视频| 国产久久久一区二区三区| 国产精品日韩av在线免费观看| 亚洲av二区三区四区| 精品无人区乱码1区二区| 香蕉av资源在线| 国产伦在线观看视频一区| 亚洲成av人片免费观看| 97人妻精品一区二区三区麻豆| 亚洲三级黄色毛片| 啪啪无遮挡十八禁网站| 久久欧美精品欧美久久欧美| 国产又黄又爽又无遮挡在线| 俄罗斯特黄特色一大片| 亚洲电影在线观看av| 亚洲第一欧美日韩一区二区三区| 亚洲熟妇熟女久久| 黄色日韩在线| 亚洲人成伊人成综合网2020| 日韩欧美三级三区| 国内精品一区二区在线观看| 亚洲人与动物交配视频| 我的女老师完整版在线观看| 欧美乱色亚洲激情| 99精品久久久久人妻精品| 国内揄拍国产精品人妻在线| 人人妻人人看人人澡| 夜夜爽天天搞| 91字幕亚洲| 乱码一卡2卡4卡精品| 男女之事视频高清在线观看| ponron亚洲| 日韩欧美在线二视频| 久久久精品欧美日韩精品| 久久香蕉精品热| 亚洲第一欧美日韩一区二区三区| 国产三级黄色录像| 麻豆成人av在线观看| 最新在线观看一区二区三区| 十八禁网站免费在线| h日本视频在线播放| 中文字幕精品亚洲无线码一区| 美女黄网站色视频| 欧美日韩瑟瑟在线播放| 狂野欧美白嫩少妇大欣赏| 色吧在线观看| 免费看光身美女| 51国产日韩欧美| 欧美黄色片欧美黄色片| 99久久九九国产精品国产免费| 美女 人体艺术 gogo| 久久久久久久久久黄片| 欧美不卡视频在线免费观看| 在线看三级毛片| 久久人人爽人人爽人人片va | 欧美一区二区国产精品久久精品| 国产精品亚洲一级av第二区| 婷婷丁香在线五月| 亚洲最大成人av| 日本黄大片高清| 日本在线视频免费播放| 午夜亚洲福利在线播放| 久久精品人妻少妇| 国产高清三级在线| 久久亚洲真实| 亚洲精品一区av在线观看| 啦啦啦韩国在线观看视频| 亚洲五月天丁香|