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

    鉛鋅礦廢石淋濾液在地下水潛水含水層中的遷移模擬研究
    ——基于有限差分法

    2018-07-04 06:57:24伍陽黃鑫張靜楊利均
    四川有色金屬 2018年2期
    關(guān)鍵詞:居民點(diǎn)運(yùn)移河流

    伍陽,黃鑫,張靜,楊利均

    (1.四川景星環(huán)境科技有限公司,四川成都 610037;2.四川農(nóng)業(yè)大學(xué),四川成都 611130)

    地下水埋藏在土壤、巖石的孔隙、裂隙和溶隙中,以各種不同形式的形態(tài)存在[1],全世界的地下水估計(jì)總量多達(dá)1.5億立方公里,幾乎占地球總水量的十分之一。目前,我國(guó)約有70%人口以地下水為主要飲用水源,全國(guó)95%以上的農(nóng)村人飲用地下水,全國(guó)40%的耕地使用地下水灌溉[2-3]。它為全球的絕大部分人口提供著水源且支持著農(nóng)業(yè)和工業(yè)活動(dòng),是一種十分重要的自然資源。然而地下水水質(zhì)在水環(huán)境中十分脆弱,一旦受到污染,不僅危害時(shí)間長(zhǎng),并且治理難度很大,因此對(duì)地下水污染的預(yù)防重要性遠(yuǎn)大于治理。

    金屬礦產(chǎn)的開發(fā)帶來巨大經(jīng)濟(jì)與能源效益,但也會(huì)出現(xiàn)一系列對(duì)環(huán)境的負(fù)面影響。其中,廢石場(chǎng)淋濾液對(duì)周圍地下水環(huán)境的影響尤為突出。鉛鋅尾礦在降雨作用下形成含有重金屬元素的淋濾液,下滲到地下水環(huán)境中隨著地下水遷移。它在污染地下水環(huán)境的同時(shí),帶來飲用水源污染風(fēng)險(xiǎn)和對(duì)居民健康的危害。通過數(shù)值模擬手段定量評(píng)估污染物在地下水環(huán)境中的遷移規(guī)律,對(duì)加強(qiáng)礦區(qū)監(jiān)測(cè)管理,以及周圍水環(huán)境的保護(hù)有著重要意義。盧曉華[4]通過建立地下水溶質(zhì)運(yùn)移模型,對(duì)廢水處理車間發(fā)生泄漏情況下地下水重金屬污染情景進(jìn)行了數(shù)值模擬與環(huán)境影響預(yù)測(cè)評(píng)價(jià),對(duì)環(huán)保部門從源頭上防治和管理潛在的地下水重金屬污染、保護(hù)地下水環(huán)境具有重要的意義;馮書順等[5]采用數(shù)值方法對(duì)礦坑井涌水量進(jìn)行了預(yù)測(cè)和模擬;鄧鼎興[6]采用有限元法對(duì)礦山地下水污染的分區(qū)治理效果進(jìn)行了數(shù)值模擬檢驗(yàn),結(jié)果表明數(shù)值模擬預(yù)測(cè)結(jié)果與實(shí)際監(jiān)測(cè)值較為接近,可為類似礦山地下水污染治理提供借鑒。憑借數(shù)值方法,在研究二維以及三維的地下水滲流和污染物的基本遷移的基礎(chǔ)之上,國(guó)內(nèi)外還出現(xiàn)了對(duì)于考慮污染物的化學(xué)反應(yīng)、生物降解、吸附、衰減等多過程的污染物遷移轉(zhuǎn)化研究[7-9],以及地下水與地表水耦合[10]、地面沉降與地下水滲流耦合等多場(chǎng)耦合模型[11]。

    本文以西部某鉛鋅礦為研究區(qū)域,將研究區(qū)域合理概化后,建立地下水水流模型與溶質(zhì)運(yùn)移模型。以有限差分方法上游加權(quán)法為基礎(chǔ),應(yīng)用MATLAB軟件模擬了廢石場(chǎng)在非正常情況下特征污染物Pb離子的遷移結(jié)果。模擬結(jié)果可以為庫區(qū)地下水環(huán)境評(píng)價(jià)和治理提供科學(xué)依據(jù)。

    1 研究區(qū)域概況

    鉛鋅礦所在的東、北及西側(cè)都以山脊為界,南側(cè)以斯弄多支溝為界,地形起伏,總體呈現(xiàn)西高東低,從西北到東南向海拔高度逐漸減小。該地區(qū)受西南支氣流、西南季風(fēng)等大氣環(huán)流的控制而具有高原亞寒帶半干旱季風(fēng)氣候的特點(diǎn),具有日照充足,空氣干燥稀薄,多大風(fēng),夏秋多雨,干濕季節(jié)分明的主要?dú)夂蛱卣?。最大年降水量?05.7mm,年平均降水量為352.38mm,最大日降雨量24.2mm,降水主要集中在6~9月,占全年降水的87%;年蒸發(fā)量1938.4mm。

    鉛鋅礦位于該區(qū)域支溝內(nèi),礦區(qū)地形坡度大,切割較強(qiáng)烈,溝谷發(fā)育且縱坡降大,第四系覆蓋嚴(yán)重,水文地質(zhì)邊界條件簡(jiǎn)單,構(gòu)成晶屑凝灰?guī)r、角礫凝灰?guī)r、溶結(jié)凝灰?guī)r與巖漿巖裂隙含水為主的水文地質(zhì)單元,屬以裂隙含水層為主的礦床。松散堆積物呈條帶狀分布在溝谷平坦部位,形成匯水谷地;溝谷斜坡上部基巖裸露,中、下部覆蓋殘坡積層。巖石類型以晶屑凝灰?guī)r、角礫凝灰?guī)r、溶結(jié)凝灰?guī)r為主,少量地帶出露花崗斑巖、閃長(zhǎng)玢巖。因此,區(qū)內(nèi)地下水主要賦存在基巖裂隙中,溝谷為其主要排泄區(qū),地下水與地表水分水嶺基本一致,地下水徑流途徑短,分散排泄,水量貧乏。地表水補(bǔ)給源以大氣降水和冰雪融水為主,降水多沿地表徑流,少量沿表土、基巖裂隙入滲,最終在排泄基準(zhǔn)面溢出。

    主要水化學(xué)類型包括Ca-Mg-HCO3,pH值介于8.19~8.38之間,地下水礦化度介于414.45~422.62mg/L,屬中性偏弱堿性,質(zhì)軟到低礦化淡水,水質(zhì)狀況良好。

    2 滲流與溶質(zhì)遷移數(shù)學(xué)模型的建立

    2.1 水文地質(zhì)概化模型

    圖1 水文模型概念圖Fig.1 Concept of hydrological model

    研究區(qū)內(nèi),在天然情況下地下水由西北向東南流動(dòng),與河流1平行的兩側(cè)分水嶺的垂向上幾乎無地下水交換,故兩側(cè)山脊線可概化為零流量邊界。西北處的為海拔較高的地形起伏較大的區(qū)域,故選取西北的天然分水嶺作為隔水邊界線。河流1與河流2所處第四系松散層,加之居民處的水井開采對(duì)邊界水位幾乎無影響,故概化為定水頭邊界。研究區(qū)頂部常年接受降水補(bǔ)給,可概化為潛水面邊界。研究區(qū)底部主要為為礦體底板,巖層的透水性弱,故概化為隔水邊界[12]。

    2.2 潛水含水層滲流控制方程

    因?yàn)楸镜貐^(qū)巖層傾角平緩,可以把地下水當(dāng)作平面二維問題。此外,由于受地形起伏控制研究區(qū)地下水流場(chǎng)形態(tài)較為復(fù)雜,地下水以泉水及溢出面形式出露,為此本次計(jì)算采用穩(wěn)定流建模。根據(jù)達(dá)西定律,對(duì)其建立起二維穩(wěn)定流形式的控制方程[13]:

    其中Kx為在橫向上的滲透系數(shù),Ky為在縱向上的滲透系數(shù),其單位為m/d。

    2.3 污染物遷移控制方程

    污染物的遷移主要受到兩方面的驅(qū)動(dòng):一是濃度梯度引起的溶質(zhì)分子的擴(kuò)散,其擴(kuò)散通量與濃度差的關(guān)系滿足費(fèi)克定律;二是由于地下水流動(dòng)引起的機(jī)械彌散。溶質(zhì)運(yùn)移控制方程[14]:

    其中θ為含水層有效孔隙率;Dx為橫向水動(dòng)力彌散系數(shù);Dy為縱向水動(dòng)力彌散系數(shù);Rd為阻滯因子;Vx代表在橫向上的地下水流速;Vy代表在縱向上的地下水流速。流速Vx或Vy由達(dá)西定律得到:

    水動(dòng)力彌散系數(shù)包括機(jī)械彌散系數(shù)與分子擴(kuò)散系數(shù):

    其中V為地下水實(shí)際流速,Vx與Vy分別為地下水速度在橫向和縱向上的速度分量,Dm為分子擴(kuò)散系數(shù)。

    2.4 初始條件與邊界條件

    在本研究中,地下水滲流場(chǎng)模型中海需要確定的邊界條件有第一類邊界與第二類邊界。

    地下水運(yùn)移的定解條件:

    (1)第一類邊界條件:地下水滲流方程中的第一類邊界為定水頭邊界,本研究區(qū)域中的河流可視為定水頭邊界:

    (2)第二類邊界條件:地下水滲流方程中第二類邊界為定流量邊界,本研究區(qū)域中的分水嶺為隔水邊界,即零流量邊界:

    溶質(zhì)運(yùn)移的定解條件包括了初始條件和邊界條件:初始條件:

    邊界條件:

    (3)第一類邊界條件:廢石浸出液入滲為連續(xù)注入式,地下水污染溶質(zhì)遷移模擬公式參考《多孔介質(zhì)滲流與污染物遷移數(shù)學(xué)模型》(仵彥卿編著)一書中平面連續(xù)注入示蹤劑,一維流動(dòng)二維彌散模型,污染源連續(xù)注入含水層,污染物濃度在各個(gè)時(shí)刻均為常數(shù):

    (4)第二類邊界條件:污染物溶質(zhì)運(yùn)移的第二類邊界為定擴(kuò)散通量邊界,此邊界與地下水滲流的第二類邊界線重合,則此邊界為零擴(kuò)散通量邊界:

    3 數(shù)值模擬及結(jié)果

    3.1 有限差分?jǐn)?shù)值解法

    有限差分法是一種將微分方程的控制方程及對(duì)應(yīng)定界條件離散,利用迭代法求解出微分方程數(shù)值解的方法。地下水流動(dòng)模型與污染物運(yùn)移模型是耦合的系統(tǒng)模型,但在穩(wěn)定流的情況下,可先求出穩(wěn)定流條件下的流場(chǎng)分布,再根據(jù)研究區(qū)域內(nèi)的各個(gè)點(diǎn)水頭求出速度場(chǎng),最后根據(jù)離散化后的溶質(zhì)運(yùn)移差分方程迭代求出污染物隨時(shí)空變化的濃度值。

    本文研究區(qū)域?yàn)楦鶕?jù)自然邊界線劃分的不規(guī)則區(qū)域。采用矩形網(wǎng)格剖分方法,共剖分了12334個(gè)節(jié)點(diǎn),每個(gè)矩形網(wǎng)格大小為30×30m2。由于研究區(qū)域被劃分為正方形網(wǎng)格且含水層各向同性,即Δx=Δy,Kx=Ky對(duì)于地下水滲流方程,其離散后差分格式為[14-15]:

    采用上游加權(quán)法,得到地下水溶質(zhì)運(yùn)移方程的顯式差分方程為:

    3.2 定解條件的確定

    如圖1所示,地下水水流模型的隔水邊界在整個(gè)研究區(qū)域的北側(cè)、西側(cè)和南側(cè),定水頭邊界為內(nèi)部河流與東側(cè)邊界處河流邊界。地下水水流模型的隔水邊界為污染物溶質(zhì)運(yùn)移的零通量邊界。根據(jù)2016年地勘資料對(duì)研究區(qū)域的污染物檢測(cè)值作為初始濃度即CO=0.02mg/L。

    3.3 參數(shù)選取

    在服務(wù)期滿(11a)前,認(rèn)為污染源為連續(xù)注入式,即淋濾液通過地表和包氣帶到達(dá)地下水系統(tǒng)后的濃度為常數(shù)。其源濃度Cq計(jì)算按公式計(jì)算[14]:

    其中sourceC為污染物下滲漏量[14]:

    其中Cp為污染物實(shí)驗(yàn)最大浸出液濃度(mg/L),α降水入滲補(bǔ)給系數(shù),A為廢石場(chǎng)滲水面積(m2),I為降水量(mm)。

    廢石場(chǎng)中的廢水滲透穿過庫底后,進(jìn)入第四系的潛水含水層,漂石含量小于2%,卵石含量40%~55%,卵石最大直徑10cm;混雜粘土與砂礫,不均勻。根據(jù)項(xiàng)目水文地質(zhì)勘查報(bào)告,評(píng)價(jià)范圍內(nèi)坡積層孔隙水的滲透系數(shù)為5.4m/d。廢石場(chǎng)所在區(qū)域的潛水含水層都是第四系坡積物、洪積物或沖積物,有效孔隙度取值0.10。采取類比的方法,確定含水層的縱向彌散度為0.3(m2/d)。根據(jù)經(jīng)驗(yàn)公式取橫向彌散度為縱向彌散度的0.1倍,即橫向彌散度為0.03(m2/d)[16]。將研究區(qū)域已有水文資料中的參數(shù)作為初始值賦給模型,通過實(shí)測(cè)的地下水水文位數(shù)據(jù)與模擬的水文數(shù)據(jù)進(jìn)行對(duì)比驗(yàn)證,采用經(jīng)過識(shí)別與驗(yàn)證水位后的參數(shù)作為模型參數(shù)來預(yù)測(cè)地下水流場(chǎng)與溶質(zhì)運(yùn)移。驗(yàn)證后的參數(shù)取值見表1。

    表1 模擬水文地質(zhì)參數(shù)及取值Tab.1 Chemical characteristics of the tested swine manure before and after leaching

    3.4 結(jié)果分析

    3.4.1 地下水流場(chǎng)分析

    以海拔高度表示研究區(qū)域內(nèi)各點(diǎn)的水位高度,穩(wěn)定流情況下的地下水的水位等值線圖見圖2。

    研究區(qū)的地下水流場(chǎng)受該地區(qū)的地形控制,一致地呈現(xiàn)出西高東低,由西北向東南遞減的規(guī)律。內(nèi)部河流的前半段為河流溝谷,地下水除自西北向東南向流動(dòng),還會(huì)向內(nèi)部河流形成的溝谷匯集流動(dòng),實(shí)現(xiàn)地下水補(bǔ)給地表河流。內(nèi)部河流后半段的沖積平層,即圖2由紅色虛線圍成三角區(qū)域,形成了一個(gè)微小的分水嶺,它的分流作用使水流不僅僅沿著內(nèi)部河流向邊界河流匯集。此外,邊界處的河流流向?yàn)樽阅舷虮?,匯集內(nèi)部河流后向下游流動(dòng),在研究區(qū)域內(nèi)水流最終以邊界處的河流下游為出口。

    圖2 水頭等值線圖Fig.2 Contour line of water head

    研究區(qū)域內(nèi)的水力梯度不大,變化不劇烈,但由于所處區(qū)域的巖層滲透系數(shù)較大,會(huì)由于較強(qiáng)的對(duì)流作用使得污染物能夠較為快速地隨著地下水遷移。另一方面,水流主要流向邊界河流,最終沿邊界河流的下游方向遷移,使得該區(qū)域地下水污染物的不易累積。

    3.4.2 污染物遷移規(guī)律分析

    分別取在1000天、11年(服務(wù)期滿)、16年(服務(wù)期滿后5年)、30年后這幾個(gè)時(shí)間節(jié)點(diǎn),得到各時(shí)間節(jié)點(diǎn)的污染物濃度分布見圖3。

    在服務(wù)期滿前,廢石場(chǎng)淋濾液被簡(jiǎn)化為濃度恒定的污染源。污染物在地下水系統(tǒng)中主要受對(duì)流作用遷移,沿著地下水主要流動(dòng)方向由西北向東南方向遷移,同時(shí)具有向內(nèi)部河流形成的河谷匯集的趨勢(shì)。根據(jù)對(duì)流場(chǎng)的分析可知,由于靠近邊界河流處的扇形沖積平層的分水嶺作用,污染物在隨著地下水遷移到內(nèi)部河流下游時(shí)將被分散到河流兩側(cè)。同時(shí)在彌散作用下,在河流2的上游地帶部分污染物擴(kuò)散到了河流2的對(duì)側(cè)。由于流速較大同時(shí)縱向彌散度較小,所形成的污染羽形狀較為扁平。在離污染源最遠(yuǎn)處的濃度最小,在污染羽中心且靠近源強(qiáng)處的濃度最大,呈現(xiàn)出中心高四周低的規(guī)律。

    在服務(wù)期滿后,廢石場(chǎng)下滲淋濾液量減少為零。污染物繼續(xù)沿著地下水主要流動(dòng)方向遷移,在服務(wù)期滿前后污染物遷移到達(dá)居民點(diǎn),并且在河流1的下游的污染帶寬度逐漸擴(kuò)大。污染羽仍然呈現(xiàn)出中心濃度高而邊緣濃度低的特征。污染物匯入河流1后,隨著河流1從南向北向下游遷移,使得研究區(qū)域內(nèi)的污染物不斷減小,所形成的污染羽長(zhǎng)度逐漸減小,中心也逐步向河流2的下游移動(dòng)。在30年后,大部分污染物已經(jīng)流出研究區(qū)域,居民點(diǎn)也不再處于污染羽的范圍內(nèi)。污染物在研究區(qū)域內(nèi)總體快速遷移而不滯留。

    3.4.3 對(duì)居民點(diǎn)的影響

    圖3 污染物在各時(shí)期的濃度分布圖Fig.3 Concentration of Pollution in different times

    居民點(diǎn)內(nèi)與整個(gè)研究區(qū)域內(nèi)污染物的最大濃度隨天數(shù)變化規(guī)律曲線如圖4所示。總體來看,研究區(qū)和居民點(diǎn)的最大污染物濃度變化規(guī)律較為一致,最大濃度隨天數(shù)變化曲線形狀呈現(xiàn)為寬峰。最初研究區(qū)域內(nèi)的最大濃度為廢石場(chǎng)下滲處的源強(qiáng)濃度,隨著時(shí)間增加最大濃度也增加,在5500天之后區(qū)域內(nèi)的最大濃度達(dá)到最大的0.2mg/L,并保持穩(wěn)定。在8700天左右,區(qū)域內(nèi)的最大濃度開始下降,到10000天左右降低到0.05mg/L度,隨著時(shí)間增加最大濃度也增加,在5500天之后區(qū)域內(nèi)的最大濃度達(dá)到最大的0.2mg/L,并保持穩(wěn)定。在8700天左右,區(qū)域內(nèi)的最大濃度開始下降,到10000天左右降低到0.05mg/L以下。由于污染物到達(dá)居民點(diǎn)需要一定的時(shí)間,所以0-3900天居民點(diǎn)內(nèi)未檢測(cè)到污染物。3900天后居民點(diǎn)內(nèi)最大濃度開始隨時(shí)間上升,在6000天左右達(dá)到最大值。居民點(diǎn)內(nèi)保持該最大濃度到8700天后開始下降,到11000天左右,居民點(diǎn)內(nèi)無法檢出污染物。

    圖4 研究區(qū)與居民區(qū)的最大鉛離子濃度變化圖Fig.4 Maximum concentration of pollution in research area and residence area

    居民點(diǎn)鉛離子的最大檢出濃度為0.0386mg/L。根據(jù)GB/T14848-93Ⅲ類水體標(biāo)準(zhǔn),Pb≤0.05mg/L。居民點(diǎn)內(nèi)的鉛離子濃度未超標(biāo)。以0.003mg/L的檢出限值為標(biāo)準(zhǔn),污染物到達(dá)居民點(diǎn)需要3910天左右。而在研究區(qū)域內(nèi),大部分時(shí)間鉛離子濃度超過了0.05mg/L,且最大濃度達(dá)到了0.2mg/L。

    4 結(jié)論

    本文根據(jù)西藏某鉛鋅礦的水文地質(zhì)條件,建立了地下水滲流與廢石淋濾液在地下水系統(tǒng)中運(yùn)移的數(shù)學(xué)模型,數(shù)值結(jié)果合理可靠,對(duì)于定量化研究地下水污染防治提供理論依據(jù)。

    (1)研究區(qū)域內(nèi)的水位與高程基本吻合,地下水流場(chǎng)受地形控制,自西北向東南的流動(dòng)為地下水的主要滲流方向,地下水沿河流1向下游流動(dòng),與河流2匯集后向河流2下游流動(dòng),有利于污染物的擴(kuò)散而不累積。

    (2)污染物沿著地下水主要流動(dòng)方向由西北向東南方向遷移,同時(shí)具有向內(nèi)部河流形成的河谷匯集的趨勢(shì),當(dāng)污染物到達(dá)河流2后將沿河流2下游方向遷移離開研究區(qū)域。整體來看,運(yùn)移過程中所形成的污染羽形狀較為扁平。在離污染源最遠(yuǎn)處的濃度最小,在污染羽中心且靠近源強(qiáng)處的濃度最大,呈現(xiàn)出中心高四周低的規(guī)律。

    (3)污染物到達(dá)居民點(diǎn)需要約3910天,污染物的最大檢出濃度為0.0386mg/L,未超過地下水Ⅲ類水質(zhì)量標(biāo)準(zhǔn)限值,因此模擬結(jié)果表明淋濾液對(duì)居民點(diǎn)沒有造成風(fēng)險(xiǎn)和危害。但在5000-8000天的時(shí)段內(nèi),需要加強(qiáng)對(duì)敏感點(diǎn)的監(jiān)測(cè),提前預(yù)防可能出現(xiàn)的風(fēng)險(xiǎn)并進(jìn)行實(shí)時(shí)相應(yīng)防治措施。

    [1]焦甜,王軍霞,唐仲華,等.地下水溶質(zhì)運(yùn)移方程有限差分格式的實(shí)證研究[J].安全與環(huán)境工程,2016,23(3):17-23.

    [2]孟展.西藏斯弄多鉛鋅銀礦床礦物學(xué)特征研究及地質(zhì)意義[D].成都理工大學(xué),2016:23-39.

    [3]陳華君,胡若鵬,鄧競(jìng)超.淋濾條件下鉛鋅尾礦中重金屬對(duì)環(huán)境影響研究[J].云南冶金,2016,45(1):56-63.

    [4]盧曉華.基于數(shù)值模擬的企業(yè)地下水重金屬污染的環(huán)境影響預(yù)測(cè)評(píng)價(jià)[J].安全與環(huán)境工程,2014,21(1):93-97.

    [5]馮書順,王國(guó)瑞,馬自強(qiáng),等.基于Visual Modflow的礦井涌水量預(yù)測(cè)模擬研究[J].煤炭技術(shù),2016,35(2):239-242.

    [6]鄧鼎興.某礦山地下水污染分區(qū)治理效果的數(shù)值模擬[J].安全與環(huán)境工程,2013,20(4):27-31.

    [7]張曉葉.基于可滲透反應(yīng)格柵的地下水中2,4-二氯酚溶質(zhì)運(yùn)移模擬研究[D].北京工業(yè)大學(xué),2015,1-2.

    [8]焦友軍,施小清,吳吉春,等.基于PHT3D的地下水中六價(jià)鈾吸附反應(yīng)運(yùn)移數(shù)值模擬[J].地下水,2015(2):8-10.

    [9]Salmon S U,Prommer H,Park J,et al.A general reactive transportmodelingframeworkforsimulatingandinterpreting groundwater14Cageandδ13C[J].WaterResourcesResearch,2015,51(1):359-376.

    [10]孫驍磊.銀北井渠結(jié)合灌區(qū)地表水地下水耦合模擬及優(yōu)化配置研究[D].寧夏大學(xué),2016,2-3.

    [11]Ahad J M,Ganeshram C L,Bryant L M.,et al.Sources of n-alkanes in an urbanized estuary:Insights from molecular distributions and compound-specific stabl e a nd radiocarbonisotopes[J].MarineChemistry,2011,126(1),239-249.

    [12]張人權(quán),等.水文地質(zhì)學(xué)基礎(chǔ)(第六版)[M].北京:地質(zhì)出版社,2011:1-184.

    [13]薛禹群,吳吉春.地下水動(dòng)力學(xué)[M].北京:中國(guó)水利水電出版社,2009:1.

    [14]仵彥卿.多孔介質(zhì)滲流與污染物遷徙數(shù)學(xué)模型[M].北京:科學(xué)出版社,2012.

    [15]薛禹群,謝春紅.地下水?dāng)?shù)值模擬[M].北京:科學(xué)出版社,2007.

    [16]GelharLW,Welty C,Rehfeldt K R.Reply to Comments on“A critical review of data on field-scale dispersion in aquifers”by S.P.Neuman[J].Water Resources Research,1993,29(6):1867-1869.

    猜你喜歡
    居民點(diǎn)運(yùn)移河流
    曲流河復(fù)合點(diǎn)壩砂體構(gòu)型表征及流體運(yùn)移機(jī)理
    東營(yíng)凹陷北帶中淺層油氣運(yùn)移通道組合類型及成藏作用
    河流
    流放自己的河流
    散文詩(2018年20期)2018-05-06 08:03:44
    基于GIS的宜興市農(nóng)村居民點(diǎn)景觀格局變化研究
    開采過程中上覆急傾斜巖層運(yùn)移規(guī)律模擬與研究
    當(dāng)河流遇見海
    川西坳陷孝泉-新場(chǎng)地區(qū)陸相天然氣地球化學(xué)及運(yùn)移特征
    濟(jì)南市農(nóng)村居民點(diǎn)用地整理潛力
    基于3S的巴彥縣農(nóng)村居民點(diǎn)內(nèi)農(nóng)地測(cè)算研究
    51午夜福利影视在线观看| 久久国产亚洲av麻豆专区| 国产国语露脸激情在线看| av免费观看日本| 中文欧美无线码| 亚洲熟女精品中文字幕| 亚洲美女黄色视频免费看| 一区二区三区精品91| 只有这里有精品99| 国产麻豆69| 日本猛色少妇xxxxx猛交久久| 激情五月婷婷亚洲| 丁香六月欧美| 在线 av 中文字幕| 女的被弄到高潮叫床怎么办| 一区二区三区四区激情视频| 亚洲国产最新在线播放| 青春草国产在线视频| 51午夜福利影视在线观看| 99久久综合免费| 男人爽女人下面视频在线观看| 欧美日本中文国产一区发布| 久久久久精品国产欧美久久久 | 日本一区二区免费在线视频| 97精品久久久久久久久久精品| 最新在线观看一区二区三区 | 一级a爱视频在线免费观看| 亚洲av综合色区一区| 中文字幕制服av| videos熟女内射| 男女无遮挡免费网站观看| 亚洲欧美成人精品一区二区| 欧美日韩亚洲综合一区二区三区_| 久久99一区二区三区| 久久99精品国语久久久| 一区二区av电影网| 国产片内射在线| 啦啦啦中文免费视频观看日本| 侵犯人妻中文字幕一二三四区| 少妇精品久久久久久久| 久久精品国产亚洲av涩爱| 90打野战视频偷拍视频| 这个男人来自地球电影免费观看 | 亚洲婷婷狠狠爱综合网| 又黄又粗又硬又大视频| 美女视频免费永久观看网站| 国产xxxxx性猛交| 免费黄频网站在线观看国产| 欧美黄色片欧美黄色片| 最近最新中文字幕大全免费视频 | 美女大奶头黄色视频| 一区二区三区精品91| 少妇人妻久久综合中文| 9色porny在线观看| 亚洲欧美中文字幕日韩二区| av线在线观看网站| 久久亚洲国产成人精品v| 成人免费观看视频高清| 亚洲精品久久久久久婷婷小说| 久久精品久久久久久噜噜老黄| 一本色道久久久久久精品综合| 伦理电影免费视频| 色婷婷久久久亚洲欧美| 日韩制服丝袜自拍偷拍| 免费看不卡的av| 国产精品.久久久| 伊人亚洲综合成人网| 成人黄色视频免费在线看| 成人国产av品久久久| 精品国产一区二区三区久久久樱花| 51午夜福利影视在线观看| 观看av在线不卡| 国产精品人妻久久久影院| 精品第一国产精品| 国产片特级美女逼逼视频| 国产无遮挡羞羞视频在线观看| 精品亚洲成a人片在线观看| 欧美激情 高清一区二区三区| 亚洲成色77777| 十八禁人妻一区二区| 激情视频va一区二区三区| 精品亚洲乱码少妇综合久久| 夫妻性生交免费视频一级片| 久久久久视频综合| 精品久久蜜臀av无| 亚洲精品成人av观看孕妇| 日韩欧美一区视频在线观看| 国产黄频视频在线观看| 国产日韩一区二区三区精品不卡| 一级a爱视频在线免费观看| 国产有黄有色有爽视频| 亚洲av成人精品一二三区| 国产成人欧美在线观看 | 国产熟女午夜一区二区三区| 欧美精品一区二区免费开放| 国产深夜福利视频在线观看| 91精品伊人久久大香线蕉| 亚洲av日韩精品久久久久久密 | 亚洲精品国产区一区二| 国产日韩一区二区三区精品不卡| 视频区图区小说| 国产av一区二区精品久久| 国产片特级美女逼逼视频| 精品少妇久久久久久888优播| 久久久久国产一级毛片高清牌| 香蕉国产在线看| 亚洲av电影在线进入| 亚洲,欧美精品.| 久久久久久久大尺度免费视频| 赤兔流量卡办理| 亚洲欧美清纯卡通| 99久久精品国产亚洲精品| 国产成人午夜福利电影在线观看| tube8黄色片| 国产爽快片一区二区三区| 一级爰片在线观看| 纯流量卡能插随身wifi吗| 超碰97精品在线观看| 又黄又粗又硬又大视频| 久久97久久精品| 18在线观看网站| 精品人妻熟女毛片av久久网站| 高清视频免费观看一区二区| 精品国产乱码久久久久久男人| 久久这里只有精品19| 亚洲av电影在线观看一区二区三区| 女性被躁到高潮视频| 久久久久久人妻| 免费在线观看完整版高清| 亚洲欧美一区二区三区黑人| 亚洲精品乱久久久久久| 一本大道久久a久久精品| 国产欧美日韩一区二区三区在线| 久久婷婷青草| 交换朋友夫妻互换小说| 熟女少妇亚洲综合色aaa.| 亚洲精品国产色婷婷电影| 天天躁日日躁夜夜躁夜夜| 午夜激情av网站| www日本在线高清视频| 秋霞在线观看毛片| 精品视频人人做人人爽| videosex国产| 国产麻豆69| 久久热在线av| 精品国产超薄肉色丝袜足j| 久久婷婷青草| 国产精品 国内视频| 爱豆传媒免费全集在线观看| 国产精品女同一区二区软件| 水蜜桃什么品种好| 日韩精品免费视频一区二区三区| 亚洲一码二码三码区别大吗| 又粗又硬又长又爽又黄的视频| av卡一久久| 麻豆乱淫一区二区| h视频一区二区三区| 国产精品蜜桃在线观看| 国产av码专区亚洲av| 啦啦啦视频在线资源免费观看| 亚洲精华国产精华液的使用体验| 老司机靠b影院| 在线 av 中文字幕| av.在线天堂| 两个人免费观看高清视频| 999久久久国产精品视频| 少妇人妻精品综合一区二区| 婷婷色综合www| 亚洲情色 制服丝袜| 亚洲伊人久久精品综合| 街头女战士在线观看网站| 天天添夜夜摸| av天堂久久9| 男女国产视频网站| 91精品国产国语对白视频| 免费久久久久久久精品成人欧美视频| 亚洲美女黄色视频免费看| 一级片免费观看大全| 国产精品一区二区在线观看99| 啦啦啦 在线观看视频| 国产精品久久久av美女十八| 无遮挡黄片免费观看| 三上悠亚av全集在线观看| 侵犯人妻中文字幕一二三四区| 欧美日韩国产mv在线观看视频| 丰满少妇做爰视频| 黄色视频在线播放观看不卡| 久久人人97超碰香蕉20202| 涩涩av久久男人的天堂| 欧美xxⅹ黑人| 国产精品亚洲av一区麻豆 | 国产精品一二三区在线看| h视频一区二区三区| 视频在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲一区中文字幕在线| 熟女av电影| 满18在线观看网站| 性少妇av在线| 久久精品久久精品一区二区三区| 国产精品蜜桃在线观看| 成年人免费黄色播放视频| 国产精品一国产av| 成人漫画全彩无遮挡| 亚洲欧美日韩另类电影网站| 一区二区日韩欧美中文字幕| 大香蕉久久成人网| 成人漫画全彩无遮挡| 国产精品免费视频内射| kizo精华| 捣出白浆h1v1| 久久久久精品久久久久真实原创| 久久久久久人妻| 午夜福利乱码中文字幕| 亚洲成av片中文字幕在线观看| 三上悠亚av全集在线观看| 国产老妇伦熟女老妇高清| 日日摸夜夜添夜夜爱| 一区二区三区四区激情视频| 亚洲综合精品二区| 久久久久视频综合| 可以免费在线观看a视频的电影网站 | 国产亚洲欧美精品永久| 亚洲综合色网址| 亚洲欧洲日产国产| 亚洲综合色网址| 久久热在线av| 国产精品av久久久久免费| 大片免费播放器 马上看| 精品午夜福利在线看| 国产精品av久久久久免费| 久久女婷五月综合色啪小说| 巨乳人妻的诱惑在线观看| 久久精品aⅴ一区二区三区四区| 男人舔女人的私密视频| 午夜福利乱码中文字幕| 啦啦啦视频在线资源免费观看| 亚洲色图综合在线观看| 久热这里只有精品99| 亚洲欧美成人综合另类久久久| 日本91视频免费播放| 免费久久久久久久精品成人欧美视频| 人妻 亚洲 视频| 亚洲 欧美一区二区三区| 两个人免费观看高清视频| 人人澡人人妻人| 超碰成人久久| 亚洲少妇的诱惑av| 久久天躁狠狠躁夜夜2o2o | 午夜福利一区二区在线看| 亚洲色图综合在线观看| 男的添女的下面高潮视频| 亚洲图色成人| 精品午夜福利在线看| 国产不卡av网站在线观看| 在线观看国产h片| 国产黄色免费在线视频| 久久久久久免费高清国产稀缺| 国产麻豆69| 美女高潮到喷水免费观看| 一边摸一边抽搐一进一出视频| bbb黄色大片| 欧美黑人欧美精品刺激| 男人爽女人下面视频在线观看| 国产极品粉嫩免费观看在线| 国产又色又爽无遮挡免| 丁香六月欧美| 日本午夜av视频| 国产一卡二卡三卡精品 | 国产成人精品在线电影| 成人影院久久| 国产av一区二区精品久久| 亚洲一区二区三区欧美精品| 国产精品熟女久久久久浪| 男人操女人黄网站| 国产爽快片一区二区三区| 久久久国产一区二区| 国产黄色免费在线视频| av免费观看日本| 亚洲免费av在线视频| 在线观看国产h片| 亚洲色图综合在线观看| 9色porny在线观看| 亚洲国产欧美在线一区| 女人精品久久久久毛片| 国产激情久久老熟女| 黑人猛操日本美女一级片| 大码成人一级视频| 免费不卡黄色视频| 欧美国产精品一级二级三级| 亚洲情色 制服丝袜| 欧美日韩视频高清一区二区三区二| 亚洲,一卡二卡三卡| netflix在线观看网站| av国产久精品久网站免费入址| 咕卡用的链子| 欧美中文综合在线视频| 亚洲av电影在线观看一区二区三区| 久久热在线av| 啦啦啦视频在线资源免费观看| 久久综合国产亚洲精品| 卡戴珊不雅视频在线播放| 色婷婷久久久亚洲欧美| 国产亚洲av高清不卡| 在线亚洲精品国产二区图片欧美| 国产又爽黄色视频| 人人妻人人澡人人爽人人夜夜| 性高湖久久久久久久久免费观看| 国产成人精品在线电影| 欧美日韩综合久久久久久| 国产精品成人在线| 中文字幕人妻熟女乱码| 操出白浆在线播放| 亚洲美女视频黄频| 久久久久精品国产欧美久久久 | 欧美少妇被猛烈插入视频| 美女福利国产在线| 97精品久久久久久久久久精品| 精品亚洲乱码少妇综合久久| 久久久精品免费免费高清| 久久av网站| 五月开心婷婷网| 亚洲精品aⅴ在线观看| 欧美日韩福利视频一区二区| 最新在线观看一区二区三区 | 国产免费又黄又爽又色| 操美女的视频在线观看| 亚洲欧洲日产国产| 好男人视频免费观看在线| 性色av一级| 欧美av亚洲av综合av国产av | 少妇人妻久久综合中文| 欧美国产精品一级二级三级| 亚洲成人手机| 国语对白做爰xxxⅹ性视频网站| 晚上一个人看的免费电影| 国产亚洲av片在线观看秒播厂| av视频免费观看在线观看| 国产一区二区激情短视频 | 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产色婷婷电影| 嫩草影院入口| 亚洲欧美日韩另类电影网站| 亚洲av成人不卡在线观看播放网 | 日本一区二区免费在线视频| xxx大片免费视频| 少妇猛男粗大的猛烈进出视频| 一边摸一边做爽爽视频免费| 美女午夜性视频免费| 大陆偷拍与自拍| 两性夫妻黄色片| 国产激情久久老熟女| 欧美 亚洲 国产 日韩一| 九九爱精品视频在线观看| 精品少妇黑人巨大在线播放| 伊人亚洲综合成人网| 亚洲成国产人片在线观看| 啦啦啦中文免费视频观看日本| 色婷婷久久久亚洲欧美| 秋霞在线观看毛片| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日本中文国产一区发布| 国产黄频视频在线观看| 国产一区二区在线观看av| 夫妻午夜视频| 亚洲欧美激情在线| 亚洲第一青青草原| 亚洲精品视频女| 久久ye,这里只有精品| 久久精品久久久久久久性| 亚洲,欧美精品.| 久久久亚洲精品成人影院| 一边亲一边摸免费视频| 免费在线观看视频国产中文字幕亚洲 | svipshipincom国产片| 女性被躁到高潮视频| 午夜影院在线不卡| 欧美精品人与动牲交sv欧美| 国产一区亚洲一区在线观看| 五月天丁香电影| 青春草视频在线免费观看| 日韩一本色道免费dvd| 国产精品欧美亚洲77777| 最近2019中文字幕mv第一页| 久久毛片免费看一区二区三区| 51午夜福利影视在线观看| 亚洲欧美成人综合另类久久久| 午夜91福利影院| 精品人妻一区二区三区麻豆| 精品国产露脸久久av麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 一本大道久久a久久精品| 搡老岳熟女国产| 国产成人精品无人区| 大片电影免费在线观看免费| av女优亚洲男人天堂| 精品人妻在线不人妻| 18在线观看网站| 成人三级做爰电影| 老司机影院成人| 欧美在线黄色| 51午夜福利影视在线观看| 国产精品 欧美亚洲| 欧美日韩精品网址| 不卡视频在线观看欧美| 69精品国产乱码久久久| 午夜日韩欧美国产| 日韩电影二区| 一区二区三区精品91| 亚洲国产欧美网| 欧美精品高潮呻吟av久久| 亚洲精品久久久久久婷婷小说| 丁香六月欧美| 精品视频人人做人人爽| 久久精品亚洲av国产电影网| 成人影院久久| 久久精品国产a三级三级三级| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦在线观看免费高清www| 777久久人妻少妇嫩草av网站| h视频一区二区三区| 观看av在线不卡| 五月天丁香电影| 精品少妇黑人巨大在线播放| 各种免费的搞黄视频| 亚洲精品第二区| 国产在线视频一区二区| 亚洲成国产人片在线观看| 免费看av在线观看网站| 香蕉丝袜av| 久久影院123| 人人妻,人人澡人人爽秒播 | 人体艺术视频欧美日本| 亚洲精品国产一区二区精华液| 建设人人有责人人尽责人人享有的| 女人久久www免费人成看片| 久久久国产一区二区| 亚洲精品在线美女| 人体艺术视频欧美日本| 国产一区二区 视频在线| 天天躁夜夜躁狠狠躁躁| 欧美精品高潮呻吟av久久| 亚洲伊人色综图| 母亲3免费完整高清在线观看| 久久97久久精品| 亚洲国产看品久久| 欧美日韩一级在线毛片| 国语对白做爰xxxⅹ性视频网站| 亚洲国产精品999| 亚洲免费av在线视频| 老汉色∧v一级毛片| 国产福利在线免费观看视频| 啦啦啦在线免费观看视频4| 精品一区二区三卡| 国产精品偷伦视频观看了| 国产片特级美女逼逼视频| 老司机影院毛片| 校园人妻丝袜中文字幕| 欧美黑人精品巨大| 少妇人妻 视频| 天堂俺去俺来也www色官网| 纵有疾风起免费观看全集完整版| 天美传媒精品一区二区| 国产又爽黄色视频| 美女主播在线视频| 又黄又粗又硬又大视频| 新久久久久国产一级毛片| 午夜av观看不卡| 丝袜美足系列| 亚洲成人国产一区在线观看 | 美女脱内裤让男人舔精品视频| 亚洲综合色网址| 欧美国产精品va在线观看不卡| 亚洲第一青青草原| 777米奇影视久久| 国产精品久久久久久精品古装| 亚洲一码二码三码区别大吗| 色播在线永久视频| 亚洲在久久综合| 欧美国产精品一级二级三级| 天天躁狠狠躁夜夜躁狠狠躁| 国产黄频视频在线观看| 国产老妇伦熟女老妇高清| 国产免费现黄频在线看| 观看美女的网站| 精品人妻在线不人妻| 久久久久视频综合| 如何舔出高潮| 亚洲av国产av综合av卡| 亚洲综合精品二区| 亚洲在久久综合| 精品一区在线观看国产| 涩涩av久久男人的天堂| 日韩伦理黄色片| 国产片特级美女逼逼视频| 岛国毛片在线播放| 在线观看国产h片| 一级爰片在线观看| 日本午夜av视频| 中文乱码字字幕精品一区二区三区| 国产伦理片在线播放av一区| 久久久久国产一级毛片高清牌| 秋霞在线观看毛片| 大陆偷拍与自拍| 9色porny在线观看| 精品人妻在线不人妻| 十分钟在线观看高清视频www| 韩国高清视频一区二区三区| 久久99一区二区三区| 99久久99久久久精品蜜桃| 亚洲成av片中文字幕在线观看| 亚洲综合色网址| 在线观看人妻少妇| 最黄视频免费看| 国产无遮挡羞羞视频在线观看| 男女之事视频高清在线观看 | 在线看a的网站| 国产99久久九九免费精品| 黄片小视频在线播放| 中国三级夫妇交换| 亚洲国产欧美日韩在线播放| 高清黄色对白视频在线免费看| 色播在线永久视频| 欧美激情 高清一区二区三区| 久久久久国产一级毛片高清牌| 精品亚洲乱码少妇综合久久| 这个男人来自地球电影免费观看 | 黄频高清免费视频| 丝袜美足系列| 久久鲁丝午夜福利片| 两个人看的免费小视频| 丝袜喷水一区| 国产成人免费观看mmmm| 纯流量卡能插随身wifi吗| 午夜日本视频在线| 一边摸一边抽搐一进一出视频| 啦啦啦 在线观看视频| 丰满乱子伦码专区| 日韩电影二区| 在线观看三级黄色| 国产野战对白在线观看| 在线观看免费日韩欧美大片| a级毛片黄视频| 19禁男女啪啪无遮挡网站| 精品第一国产精品| 97在线人人人人妻| 日韩不卡一区二区三区视频在线| 99热全是精品| 亚洲av日韩在线播放| 热99国产精品久久久久久7| 国产精品一区二区在线观看99| 成人影院久久| 久久狼人影院| 最近最新中文字幕免费大全7| 伊人久久大香线蕉亚洲五| 看十八女毛片水多多多| 人人澡人人妻人| 欧美黑人精品巨大| 日韩大片免费观看网站| 免费黄网站久久成人精品| 免费观看a级毛片全部| 久久久精品区二区三区| 久久久亚洲精品成人影院| 一本久久精品| 国产高清不卡午夜福利| 成年动漫av网址| 精品久久久久久电影网| 极品少妇高潮喷水抽搐| 七月丁香在线播放| 亚洲国产看品久久| svipshipincom国产片| 国产一卡二卡三卡精品 | 水蜜桃什么品种好| 99精国产麻豆久久婷婷| 日韩欧美精品免费久久| 女人久久www免费人成看片| 久久97久久精品| 成年人免费黄色播放视频| 黑人猛操日本美女一级片| 国产又爽黄色视频| 色吧在线观看| 国产xxxxx性猛交| 亚洲国产欧美在线一区| 99久久人妻综合| 精品国产乱码久久久久久小说| 一级黄片播放器| 色精品久久人妻99蜜桃| 欧美精品亚洲一区二区| 满18在线观看网站| 免费日韩欧美在线观看| 亚洲欧美精品自产自拍| 一本久久精品| 最近最新中文字幕大全免费视频 | 毛片一级片免费看久久久久| 大香蕉久久网| 国产成人a∨麻豆精品| 国产欧美亚洲国产| 丝袜美足系列| 欧美激情高清一区二区三区 | 岛国毛片在线播放| 国产av国产精品国产| 九草在线视频观看| 久久国产亚洲av麻豆专区| 国产免费现黄频在线看| av电影中文网址| 免费女性裸体啪啪无遮挡网站| 在线 av 中文字幕| 美女大奶头黄色视频| 中文乱码字字幕精品一区二区三区| 久久久久久久久久久久大奶| 国产乱来视频区| 如日韩欧美国产精品一区二区三区| 如何舔出高潮| 日日撸夜夜添| 久久精品国产a三级三级三级|