摘要:使用Geotomographie公司IPG跨孔地震成像系統(tǒng)獲取井間地震波速分布情況,使用統(tǒng)計學方法將地震數(shù)據(jù)與鉆孔滲透系數(shù)相關聯(lián),通過調(diào)整地下水模型的水動力邊界條件,使該模型與多個鉆孔觀測到的水頭相符,使用粒子追蹤來模擬示蹤劑到達時間,模型計算所得的示蹤劑到達時間與現(xiàn)場觀測結果基本一致?
地下水
是指地面以下賦存于土壤和巖石空隙中的水
往往具有
水質(zhì)好、分布廣、便于開采等特征,
是生活飲用水、工農(nóng)業(yè)生產(chǎn)用水的重要水源
同時,地下水具有一定的緩慢流動性
在巖土孔隙中做復雜運動
地下水流的運移特性
也使得地下水成了污染物的搬運工
且具有隱蔽性、長期性和難恢復性等特點
圖片來源:南水北調(diào)與水利科技
地下水怎么動?往哪動?
就成了評估污染風險的重要問題。
識別地下含水層的空間水動力學特性,如孔隙度和滲透系數(shù)及其相關結構等,在污染物風險評估中起著重要的作用。
德國Geotomographie公司的Thomas Fechner等在論文《A hydro-geophysical tomographic approach to model groundwater flow and transport》中,基于地震-水動力模型,成功地模擬了在德國北部的Belau試驗場進行的野外示蹤實驗。
水文地球物理層析方法在地下水流和運移建模中的應用
在小尺度上,平均沉積物性質(zhì)可能不適合用于描述污染物遷移的地下構造特性,而是需要詳細的含水層特征來描繪優(yōu)先水流路徑。此外,河流含水層沉積物的水動力學特性可能在幾十年內(nèi)在小距離上發(fā)生變化,只有少數(shù)高度連通且低阻的區(qū)域可能主導整個區(qū)域地下水流動狀態(tài)。
為彌補疏距鉆孔橫向分辨率差的問題,可以采用地球物理跨孔方法以較高的分辨率來描述地下特征。將鉆孔樣品資料作為直接測量資料,將具有優(yōu)異的垂直分辨率和“軟”高分辨率的地球物理資料,如層析成像數(shù)據(jù),作為間接測量資料,兩者相結合,可以克服在獲取地下環(huán)境特征以描述地下水流動與遷移溶質(zhì)運移過程中所遇到的一些主要問題。
在此案例中,作者將井間地震層析成像數(shù)據(jù)與水動力鉆孔數(shù)據(jù)相結合并建模,模擬在德國北部Belau試驗場進行的野外示蹤實驗。
示蹤劑的到達記錄由不同的多級鉆孔提供。鉆孔地震波速、滲透系數(shù)、地震波速與孔隙度的巖石物理學關系等三者之間的基礎關系被建立起來,并相結合建立了水動力地下水模型。
試驗場地情況
Belau環(huán)境研究實驗場成立于1990,其研究項目主要為淺層冰川沉積物中的小尺度流動和輸運過程以及地球化學參數(shù)的波動。研究開始以來約完成了13個鉆孔的鉆探。
圖1 鉆孔位置和場地簡圖
鉆孔被設置為多級測井(M)、全屏蔽測井(B)和測壓計(D)。5個多級測井可通過永久取樣裝置在含水層不同深度進行取樣。為了確定含水層的水動力學特性,鉆探過程中每米提取600個土壤樣品,并識別地質(zhì)序列變化。取約300個樣品測定粒度分布并計算水動力傳導系數(shù)。
圖2 地質(zhì)構造簡圖(根據(jù)Scheytt 1994)
根據(jù)Piotrowski (1994),該場地的沉積層是European Elster, Saale 和 Weichsel冰期的冰川沉積物。作為一個主要的地下不整合構造,一個Elster冰期的冰川通道由北向南貫通。Bornhoeved 湖區(qū)的湖泊—Belau湖是5個湖泊之一—與該通道對齊。第四系底下沉積層可根據(jù)水動力學性質(zhì)分為兩個主要的含水層系統(tǒng)-I和II。
位于Belau東側的Belau試驗場為調(diào)查上層含水層系統(tǒng)I而建立,該區(qū)域由三個總厚度約40米的融水層組成。該融水層由Saale-Sander, Wankendorfer-Sander和 Kaluebber-Sander沙質(zhì)沉積層構成,其中Kaluebber-Sander在該位置為不飽和層。
場地示蹤劑實驗
為描述實驗場小尺度水流和輸運過程,Rumohr (1995)進行了一個自然示蹤劑梯度實驗。示蹤劑實驗自1994年5月2日(零日)開始,至1995年10月共持續(xù)約200天。1小時內(nèi)約1500g伊紅注入至上游鉆孔B9,注入位置自海拔11m至24m。起始示蹤劑濃度約1g/L,下游所有多級井進行采樣對示蹤劑注射進行跟蹤。
使用一維方法計算流速Va和縱向彌散系數(shù)D1,為說明示蹤劑穿透曲線與深度的相關性,繪制與平流運輸相關的示蹤劑到達時間曲線。
圖3 與平流輸運相關的示蹤劑到達時間(三角形為主峰,圓形為次峰)
地震-水動力聯(lián)合模型
首先,采用Geotomographie的IPG地震層析成像系統(tǒng),在鉆孔B9和M7之間,沿平行于地下水流的方向進行了一次地震層析成像,獲得了鉆孔之間的地震波速度結構分布。
IPG跨孔地震成像系統(tǒng)是一款應用廣泛,以高精度低故障率而業(yè)內(nèi)口碑好的高精度的跨孔P波地震層析成像設備,適用范圍探測包括但不限于:
圖4 孔間波速結構
將高分辨率層析成像數(shù)據(jù)和水動力學鉆孔數(shù)據(jù)相結合,是將地下構造的非均勻性納入局部水文地質(zhì)場地模型的基礎。
接著,根據(jù)對比分析鉆孔地震波速與滲透系數(shù)數(shù)據(jù),將上一步所獲得的鉆空間地震波速轉換成滲透系數(shù),使用篩分析方法所獲得的滲透系數(shù)“硬”數(shù)據(jù)與地震資料“軟”數(shù)據(jù)相結合,并應用于層析成像,其成果即為B9和M7之間的高分辨率滲透系數(shù)結構。
圖5 地震波速與Sanders-沉積層相互關系
建立了包含高分辨率數(shù)據(jù)的水動力概念模型。地下水的流向向西,往湖泊方向。如前所述,湖區(qū)附近的水動力學條件是復雜的,在對進出關注區(qū)域的水流的不同組分進行量化方面有許多不確定之處。
圖6 地下水概念模型
通過地震層析成像進行研究的區(qū)域大致位于計算的地下水流動網(wǎng)格的中間。為將水動力學數(shù)值分配給層析成像區(qū)域以外的邊緣區(qū)域,對地震波速等級進行了外推。雖然這些外部區(qū)域影響整個流態(tài),但關注的焦點僅限于由層析成像圖像覆蓋的區(qū)域。根據(jù)地震波速等級分布,使用SPP方法,生成一個孔隙度分布成果和多個模式下的滲透系數(shù)場成果
圖7 根據(jù)SPP模式2獲得的滲透系數(shù)分布
試驗成果
在參考邊界條件下,分別采用滲透系數(shù)場的十個不同模式,所得到的鉆孔M10、M9和M7的示蹤劑顆粒到達時間與深度的相關曲線。可以看出,計算所得的到達時間與現(xiàn)場觀測結果符合的很好。
圖8 邊界條件下(10種模式),M10,M9和M7鉆孔計算所得的示蹤劑粒子到達時間。主峰為三角形,次峰為圓形,計算所得到達時間為實線
15m以上的Wankendorfer砂層運輸速度快,在這里模型再現(xiàn)了隨海拔高度的增加到達時間越快的趨勢。10到15m的Saale砂層顆粒到達時間相對很慢,特別是觀測鉆孔M9附近。對該區(qū)域10-15m的現(xiàn)場測試結果的對比表明,該區(qū)域地下水速度被高估了,導致計算所得的顆粒到達時間要遠遠早于現(xiàn)場數(shù)據(jù)。該區(qū)域對應于有較高滲透系數(shù)變化性的細沙質(zhì)物質(zhì),其滲透系數(shù)潛在的小尺度波動無法用地震層析成像來解決。因此稍低的平均滲透系數(shù)可以補償較快的模型到達時間。
上述案例基于地震-水動力模型,成功地模擬了一次場地實驗。使用高分辨率地震層析成像獲取井間地震波速分布情況,使用統(tǒng)計學方法將地震數(shù)據(jù)與鉆孔滲透系數(shù)相關聯(lián),通過調(diào)整地下水模型的水動力邊界條件,使該模型與多個鉆孔觀測到的水頭相符,粒子追蹤被用于模擬示蹤劑到達時間,模型計算所得的示蹤劑到達時間與現(xiàn)場觀測結果基本一致。
此項目的研究成果,為在小尺度上使用地球物理數(shù)據(jù)來模擬和預測地下水流量和輸運提供了一個研究思路和方向,具有一定的借鑒意義。
翻譯:王雷
參考文獻: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.
主站蜘蛛池模板: 韶山市| 金平| 深州市| 永城市| 瑞安市| 武宣县| 潮安县| 瓦房店市| 忻城县| 常德市| 翼城县| 白河县| 垫江县| 武夷山市| 香格里拉县| 开远市| 罗山县| 博罗县| 宁河县| 湖南省| 上杭县| 大田县| 天峨县| 涡阳县| 江华| 喀喇沁旗| 文水县| 广宗县| 奉新县| 乌恰县| 稻城县| 鄢陵县| 讷河市| 灯塔市| 石嘴山市| 淮阳县| 屏南县| 苍梧县| 镇远县| 兴和县| 祁门县|