CAOPONR在线视频进入,亚洲国产不卡久久久久久,狠狠狼鲁亚洲综合网,欧美色偷偷

400-700-9998
歐美大地 工程物探
工程物探
科學(xué)建立地下水流和運(yùn)移模型
發(fā)布時(shí)間:2020-09-18 瀏覽次數(shù):48930 來源:歐美大地

 

摘要:使用Geotomographie公司IPG跨孔地震成像系統(tǒng)獲取井間地震波速分布情況,使用統(tǒng)計(jì)學(xué)方法將地震數(shù)據(jù)與鉆孔滲透系數(shù)相關(guān)聯(lián),通過調(diào)整地下水模型的水動(dòng)力邊界條件,使該模型與多個(gè)鉆孔觀測(cè)到的水頭相符,使用粒子追蹤來模擬示蹤劑到達(dá)時(shí)間,模型計(jì)算所得的示蹤劑到達(dá)時(shí)間與現(xiàn)場(chǎng)觀測(cè)結(jié)果基本一致?


地下水

是指地面以下賦存于土壤和巖石空隙中的水

往往具有

水質(zhì)好、分布廣、便于開采等特征,

是生活飲用水、工農(nóng)業(yè)生產(chǎn)用水的重要水源

同時(shí),地下水具有一定的緩慢流動(dòng)性

在巖土孔隙中做復(fù)雜運(yùn)動(dòng)

地下水流的運(yùn)移特性

也使得地下水成了污染物的搬運(yùn)工

且具有隱蔽性、長(zhǎng)期性和難恢復(fù)性等特點(diǎn)

圖片來源:南水北調(diào)與水利科技

地下水怎么動(dòng)?往哪動(dòng)?

就成了評(píng)估污染風(fēng)險(xiǎn)的重要問題。

 

識(shí)別地下含水層的空間水動(dòng)力學(xué)特性,如孔隙度和滲透系數(shù)及其相關(guān)結(jié)構(gòu)等,在污染物風(fēng)險(xiǎn)評(píng)估中起著重要的作用。

德國(guó)Geotomographie公司的Thomas Fechner等在論文《A hydro-geophysical tomographic approach to model groundwater flow and transport》中,基于地震-水動(dòng)力模型,成功地模擬了在德國(guó)北部的Belau試驗(yàn)場(chǎng)進(jìn)行的野外示蹤實(shí)驗(yàn)。
 

 
 

水文地球物理層析方法在地下水流和運(yùn)移建模中的應(yīng)用

 
 

 

在小尺度上,平均沉積物性質(zhì)可能不適合用于描述污染物遷移的地下構(gòu)造特性,而是需要詳細(xì)的含水層特征來描繪優(yōu)先水流路徑。此外,河流含水層沉積物的水動(dòng)力學(xué)特性可能在幾十年內(nèi)在小距離上發(fā)生變化,只有少數(shù)高度連通且低阻的區(qū)域可能主導(dǎo)整個(gè)區(qū)域地下水流動(dòng)狀態(tài)。

為彌補(bǔ)疏距鉆孔橫向分辨率差的問題,可以采用地球物理跨孔方法以較高的分辨率來描述地下特征。將鉆孔樣品資料作為直接測(cè)量資料,將具有優(yōu)異的垂直分辨率和“軟”高分辨率的地球物理資料,如層析成像數(shù)據(jù),作為間接測(cè)量資料,兩者相結(jié)合,可以克服在獲取地下環(huán)境特征以描述地下水流動(dòng)與遷移溶質(zhì)運(yùn)移過程中所遇到的一些主要問題。

 

在此案例中,作者將井間地震層析成像數(shù)據(jù)與水動(dòng)力鉆孔數(shù)據(jù)相結(jié)合并建模,模擬在德國(guó)北部Belau試驗(yàn)場(chǎng)進(jìn)行的野外示蹤實(shí)驗(yàn)。

示蹤劑的到達(dá)記錄由不同的多級(jí)鉆孔提供。鉆孔地震波速、滲透系數(shù)、地震波速與孔隙度的巖石物理學(xué)關(guān)系等三者之間的基礎(chǔ)關(guān)系被建立起來,并相結(jié)合建立了水動(dòng)力地下水模型。

 

試驗(yàn)場(chǎng)地情況

Belau環(huán)境研究實(shí)驗(yàn)場(chǎng)成立于1990,其研究項(xiàng)目主要為淺層冰川沉積物中的小尺度流動(dòng)和輸運(yùn)過程以及地球化學(xué)參數(shù)的波動(dòng)。研究開始以來約完成了13個(gè)鉆孔的鉆探。

圖1 鉆孔位置和場(chǎng)地簡(jiǎn)圖

鉆孔被設(shè)置為多級(jí)測(cè)井(M)、全屏蔽測(cè)井(B)和測(cè)壓計(jì)(D)。5個(gè)多級(jí)測(cè)井可通過永久取樣裝置在含水層不同深度進(jìn)行取樣。為了確定含水層的水動(dòng)力學(xué)特性,鉆探過程中每米提取600個(gè)土壤樣品,并識(shí)別地質(zhì)序列變化。取約300個(gè)樣品測(cè)定粒度分布并計(jì)算水動(dòng)力傳導(dǎo)系數(shù)。

圖2 地質(zhì)構(gòu)造簡(jiǎn)圖(根據(jù)Scheytt 1994)

根據(jù)Piotrowski (1994),該場(chǎng)地的沉積層是European Elster, Saale 和 Weichsel冰期的冰川沉積物。作為一個(gè)主要的地下不整合構(gòu)造,一個(gè)Elster冰期的冰川通道由北向南貫通。Bornhoeved 湖區(qū)的湖泊—Belau湖是5個(gè)湖泊之一—與該通道對(duì)齊。第四系底下沉積層可根據(jù)水動(dòng)力學(xué)性質(zhì)分為兩個(gè)主要的含水層系統(tǒng)-I和II。

位于Belau東側(cè)的Belau試驗(yàn)場(chǎng)為調(diào)查上層含水層系統(tǒng)I而建立,該區(qū)域由三個(gè)總厚度約40米的融水層組成。該融水層由Saale-Sander, Wankendorfer-Sander和 Kaluebber-Sander沙質(zhì)沉積層構(gòu)成,其中Kaluebber-Sander在該位置為不飽和層。

 

場(chǎng)地示蹤劑實(shí)驗(yàn)

為描述實(shí)驗(yàn)場(chǎng)小尺度水流和輸運(yùn)過程,Rumohr (1995)進(jìn)行了一個(gè)自然示蹤劑梯度實(shí)驗(yàn)。示蹤劑實(shí)驗(yàn)自1994年5月2日(零日)開始,至1995年10月共持續(xù)約200天。1小時(shí)內(nèi)約1500g伊紅注入至上游鉆孔B9,注入位置自海拔11m至24m。起始示蹤劑濃度約1g/L,下游所有多級(jí)井進(jìn)行采樣對(duì)示蹤劑注射進(jìn)行跟蹤。

使用一維方法計(jì)算流速Va和縱向彌散系數(shù)D1,為說明示蹤劑穿透曲線與深度的相關(guān)性,繪制與平流運(yùn)輸相關(guān)的示蹤劑到達(dá)時(shí)間曲線。

圖3 與平流輸運(yùn)相關(guān)的示蹤劑到達(dá)時(shí)間(三角形為主峰,圓形為次峰)

 

地震-水動(dòng)力聯(lián)合模型

首先,采用Geotomographie的IPG地震層析成像系統(tǒng),在鉆孔B9和M7之間,沿平行于地下水流的方向進(jìn)行了一次地震層析成像,獲得了鉆孔之間的地震波速度結(jié)構(gòu)分布。

 

IPG跨孔地震成像系統(tǒng)是一款應(yīng)用廣泛,以高精度低故障率而業(yè)內(nèi)口碑好的高精度的跨孔P波地震層析成像設(shè)備,適用范圍探測(cè)包括但不限于:

  • 地下巖溶、古洞、空洞、埋設(shè)物、礦區(qū)采空區(qū);
  • 查明地下構(gòu)造;
  • 滲透帶、水流通道和防衛(wèi),圈定破碎帶位置和范圍;
  • 建筑物地基、鐵路公路路基等不良地質(zhì)體監(jiān)測(cè)、水電站、核電站選址勘查;
  • 樁基質(zhì)量檢測(cè)、庫壩灌漿幄幕和高噴防滲板墻質(zhì)量檢測(cè),水庫、壩基檢漏等?

圖4 孔間波速結(jié)構(gòu)

將高分辨率層析成像數(shù)據(jù)和水動(dòng)力學(xué)鉆孔數(shù)據(jù)相結(jié)合,是將地下構(gòu)造的非均勻性納入局部水文地質(zhì)場(chǎng)地模型的基礎(chǔ)。

接著,根據(jù)對(duì)比分析鉆孔地震波速與滲透系數(shù)數(shù)據(jù),將上一步所獲得的鉆空間地震波速轉(zhuǎn)換成滲透系數(shù),使用篩分析方法所獲得的滲透系數(shù)“硬”數(shù)據(jù)與地震資料“軟”數(shù)據(jù)相結(jié)合,并應(yīng)用于層析成像,其成果即為B9和M7之間的高分辨率滲透系數(shù)結(jié)構(gòu)。

圖5 地震波速與Sanders-沉積層相互關(guān)系

建立了包含高分辨率數(shù)據(jù)的水動(dòng)力概念模型。地下水的流向向西,往湖泊方向。如前所述,湖區(qū)附近的水動(dòng)力學(xué)條件是復(fù)雜的,在對(duì)進(jìn)出關(guān)注區(qū)域的水流的不同組分進(jìn)行量化方面有許多不確定之處。

圖6 地下水概念模型

通過地震層析成像進(jìn)行研究的區(qū)域大致位于計(jì)算的地下水流動(dòng)網(wǎng)格的中間。為將水動(dòng)力學(xué)數(shù)值分配給層析成像區(qū)域以外的邊緣區(qū)域,對(duì)地震波速等級(jí)進(jìn)行了外推。雖然這些外部區(qū)域影響整個(gè)流態(tài),但關(guān)注的焦點(diǎn)僅限于由層析成像圖像覆蓋的區(qū)域。根據(jù)地震波速等級(jí)分布,使用SPP方法,生成一個(gè)孔隙度分布成果和多個(gè)模式下的滲透系數(shù)場(chǎng)成果

圖7 根據(jù)SPP模式2獲得的滲透系數(shù)分布

 

試驗(yàn)成果

在參考邊界條件下,分別采用滲透系數(shù)場(chǎng)的十個(gè)不同模式,所得到的鉆孔M10、M9和M7的示蹤劑顆粒到達(dá)時(shí)間與深度的相關(guān)曲線??梢钥闯觯?jì)算所得的到達(dá)時(shí)間與現(xiàn)場(chǎng)觀測(cè)結(jié)果符合的很好。

圖8 邊界條件下(10種模式),M10,M9和M7鉆孔計(jì)算所得的示蹤劑粒子到達(dá)時(shí)間。主峰為三角形,次峰為圓形,計(jì)算所得到達(dá)時(shí)間為實(shí)線

15m以上的Wankendorfer砂層運(yùn)輸速度快,在這里模型再現(xiàn)了隨海拔高度的增加到達(dá)時(shí)間越快的趨勢(shì)。10到15m的Saale砂層顆粒到達(dá)時(shí)間相對(duì)很慢,特別是觀測(cè)鉆孔M9附近。對(duì)該區(qū)域10-15m的現(xiàn)場(chǎng)測(cè)試結(jié)果的對(duì)比表明,該區(qū)域地下水速度被高估了,導(dǎo)致計(jì)算所得的顆粒到達(dá)時(shí)間要遠(yuǎn)遠(yuǎn)早于現(xiàn)場(chǎng)數(shù)據(jù)。該區(qū)域?qū)?yīng)于有較高滲透系數(shù)變化性的細(xì)沙質(zhì)物質(zhì),其滲透系數(shù)潛在的小尺度波動(dòng)無法用地震層析成像來解決。因此稍低的平均滲透系數(shù)可以補(bǔ)償較快的模型到達(dá)時(shí)間。

 

上述案例基于地震-水動(dòng)力模型,成功地模擬了一次場(chǎng)地實(shí)驗(yàn)。使用高分辨率地震層析成像獲取井間地震波速分布情況,使用統(tǒng)計(jì)學(xué)方法將地震數(shù)據(jù)與鉆孔滲透系數(shù)相關(guān)聯(lián),通過調(diào)整地下水模型的水動(dòng)力邊界條件,使該模型與多個(gè)鉆孔觀測(cè)到的水頭相符,粒子追蹤被用于模擬示蹤劑到達(dá)時(shí)間,模型計(jì)算所得的示蹤劑到達(dá)時(shí)間與現(xiàn)場(chǎng)觀測(cè)結(jié)果基本一致。

此項(xiàng)目的研究成果,為在小尺度上使用地球物理數(shù)據(jù)來模擬和預(yù)測(cè)地下水流量和輸運(yùn)提供了一個(gè)研究思路和方向,具有一定的借鑒意義。

 

翻譯:王雷

參考文獻(xiàn):Thomas Fechner.

A hydro-geophysical tomographic approach to model groundwater flow and transport[A]. International Conference on Engineering Geophysics[C]. Al Ain, United Arab Emirates: 2017.

 
關(guān)注歐美大地儀器公眾號(hào),獲取最新資訊
為您推薦