李雅靜 郭迎濤 劉立才 張訓(xùn)玉

摘 要：地下水封石油洞庫是通過人工開挖建于地下，利用穩(wěn)定地下水的水封作用密封儲(chǔ)存在洞室內(nèi)的石油。具有占地少、投資少、損耗少等優(yōu)點(diǎn)，但對(duì)地下水環(huán)境影響相對(duì)較大。為評(píng)價(jià)地下水封石油洞庫施工期和運(yùn)營(yíng)期對(duì)地下水環(huán)境的影響，以某大型地下水封石油儲(chǔ)庫為研究對(duì)象，結(jié)合現(xiàn)場(chǎng)試驗(yàn)數(shù)據(jù)，利用GMS軟件，建立了非均質(zhì)、各向異性、三維非穩(wěn)定地下水流模型和溶質(zhì)運(yùn)移模型。研究結(jié)果表明：有水幕系統(tǒng)時(shí)預(yù)測(cè)水位高于設(shè)計(jì)水位3 m以上，可以保證洞庫的水封效果，且石油局限于洞室周圍不外泄，為后期制定地下水污染防治措施及長(zhǎng)期監(jiān)測(cè)計(jì)劃提供了依據(jù)。

關(guān)鍵詞：地下水封石油洞庫；滲流場(chǎng)；涌水量；模擬計(jì)算

Groundwater seepage field and water inflow simulation calculation for an underground water sealing petroleum storage cavern

LI Yajing1， GUO Yingtao1， LIU Licai1， ZHANG Xunyu2

（1.Beijing Z.D.H.K Environmental Science & Technology Co.， Ltd.， Beijing 100028， China；

2.Beijing Institute of Ecological Geology， Beijing100011， China）

Abstract： The underground water sealing petroleum storage cavern is built underground by manual excavation， which stores oil by using stable groundwater water sealing effect. It has the advantages of less land occupation， less investment and less loss， but has relatively greater influence on the groundwater environment. In order to assess the groundwater environmental impact of an underground water sealing petroleum storage cavern during the construction period and the operation period， the heterogeneity， anisotropy and three-dimensional unsteady groundwater flow model is established by using GMS software and in-situ test data in the engineering background of a large water sealing petroleum storage caverns in China. The results show that the predicted water level is over 3 meters above the design water level with the water curtain system， which ensures the water sealing effect of the cavern， and oil is limited to the surrounding area of the cave and does not leak out. The study provides a basis for the development of groundwater pollution prevention measures and long-term monitoring plans in the later period.

Keywords： an underground water sealing petroleum storage cavern; seepage field; pouring quantity; groundwater numerical simulation

石油是國(guó)家的經(jīng)濟(jì)命脈，現(xiàn)代工業(yè)的“血液”，是國(guó)家生存和發(fā)展的重要戰(zhàn)略資源。我國(guó)從1993年開始成為石油凈進(jìn)口國(guó)，截至2020年，我國(guó)石油對(duì)外依存度高達(dá)72%，投入運(yùn)營(yíng)的大型地下水封洞庫僅3座（李印等，2022）。建設(shè)國(guó)家石油儲(chǔ)備是保證我國(guó)能源安全的重要舉措，在石油儲(chǔ)備體系建設(shè)中，地下洞庫儲(chǔ)備與地上油罐儲(chǔ)備相比，具有占地少、投資少、損耗少、污染小、運(yùn)營(yíng)管理費(fèi)用低、安全性能較高、裝卸速度較快等優(yōu)點(diǎn)，更有利于儲(chǔ)備。世界上地質(zhì)條件適宜的國(guó)家均趨向于建造地下水封石油洞庫（王婷婷，2019；王者超等，2022）。

在地下水封洞庫建設(shè)方面，中國(guó)尚處于起步階段。地下水封石油洞庫的建設(shè)對(duì)洞庫區(qū)域的工程地質(zhì)條件和水文地質(zhì)條件要求較高（郭書太等，2006；彭振華等，2011；蔡紅飚等，2000；時(shí)洪斌，2010），尤其是掌握地下水滲流場(chǎng)，對(duì)洞庫地下水環(huán)境影響評(píng)價(jià)和后期運(yùn)營(yíng)管理起到至關(guān)重要的作用。地下水環(huán)境影響評(píng)價(jià)是環(huán)境影響評(píng)價(jià)中的難點(diǎn)與重點(diǎn)，近年來眾多學(xué)者采用地下水?dāng)?shù)值模擬技術(shù)分析建設(shè)項(xiàng)目對(duì)地下水環(huán)境的影響（梁斌等，2017；洪淑娜，2023；胡成等，2022；韓曼，2007；李術(shù)才等，2013）。

目前，在世界范圍內(nèi)被廣泛應(yīng)用的地下水?dāng)?shù)值模擬軟件有Visual MODFLOW、Groundwater Modeling System（簡(jiǎn)稱GMS）、FEFLOW等。GMS是目前國(guó)際上先進(jìn)的、綜合性的地下水模擬軟件包，由MODFLOW、MODPATH、MT3D、FEMWATER、PEST、MAP、SUBSUR-FACE CHARACTERIZATION、Borehole Data、Triangulated Irregular Nets、Solid、GEO-STATISTICS等模塊組成的可視化三維地下水模擬軟件包，可進(jìn)行水流模擬、溶質(zhì)運(yùn)移模擬、反應(yīng)運(yùn)移模擬，建立三維地層實(shí)體，進(jìn)行鉆孔數(shù)據(jù)管理、二維（三維）地質(zhì)統(tǒng)計(jì)，可視化和打印二維（三維）模擬結(jié)果。它是唯一支持TIN、立體圖、鉆孔數(shù)據(jù)、2D和3D地質(zhì)統(tǒng)計(jì)、2D和3D有限差的集成系統(tǒng)。

本文以某地下水封石油洞庫為研究對(duì)象，利用GMS軟件建立了地下水模型，模擬了地下水封石油洞庫的滲流場(chǎng)、涌水量和溶質(zhì)運(yùn)移變化情況，驗(yàn)證了地下水封石油洞庫的儲(chǔ)油原理，為地下水封石油洞庫的建設(shè)提供了理論保證。

1 ?地下水封石油洞庫儲(chǔ)油原理

地下水封石油洞庫是通過人工在地下巖石中開挖形成的。通過水幕巷道向巖體內(nèi)注水使水壓高于內(nèi)部油氣壓力，從而形成密封的存儲(chǔ)設(shè)施。通過氣密性試驗(yàn)成果設(shè)定安全運(yùn)營(yíng)水位，當(dāng)監(jiān)測(cè)到水幕系統(tǒng)水位低于此水位時(shí)，從注水井向水幕中補(bǔ)水直至水位到達(dá)要求（夏喜林等，2004）。

運(yùn)營(yíng)期圍巖裂隙水不斷涌入洞內(nèi)，在主洞室底層形成水墊層。為調(diào)節(jié)水墊層的高度，防止儲(chǔ)油空間被地下水壓縮，每個(gè)洞罐設(shè)置一個(gè)泵坑定期抽水到地面，處理達(dá)標(biāo)后排放。

地下水封石油洞庫儲(chǔ)油原理見圖1。

2 ?模擬區(qū)水文地質(zhì)條件

根據(jù)項(xiàng)目前期水文地質(zhì)勘察成果，模擬區(qū)地下水類型從上到下分別為第四系松散巖類孔隙水、基巖裂隙水，基巖裂隙水又可分為淺層風(fēng)化裂隙水及深層脈狀裂隙水。地下工程所在層位的地下水類型為基巖裂隙水（圖2）。

第四系松散巖類孔隙水主要賦存于第四系松散覆蓋層中，巖性主要為沖積物夾雜淤泥質(zhì)、粉砂質(zhì)、黏土混合物，富水性較強(qiáng)，含水層單井涌水量小于50 m3·d-1，主要分布于庫址區(qū)的東部和東南部，厚度0.5～7 m不等。

基巖裂隙水主要賦存在中粒二長(zhǎng)花崗巖強(qiáng)風(fēng)化帶及中風(fēng)化帶淺部。強(qiáng)風(fēng)化帶中的基巖裂隙水與上部松散層孔隙水聯(lián)系較為密切，主要接受大氣降水的補(bǔ)給，富水性等級(jí)為中—弱；中風(fēng)化帶中深部及其以下巖體中，含水層富水性弱。壓水試驗(yàn)滲透系數(shù)為10-3～10-5 m·d-1，滲透系數(shù)總體隨深度有減小的趨勢(shì)，且洞庫深部滲透系數(shù)變化范圍較小，無大的導(dǎo)水裂隙。

第四系松散巖類孔隙水接受大氣降水入滲補(bǔ)給和淺層風(fēng)化裂隙水側(cè)向補(bǔ)給；受地形控制，地下水向溝谷方向徑流，或垂向補(bǔ)給淺層風(fēng)化裂隙水；以側(cè)向補(bǔ)給河水或出露成泉的形式排泄。

淺層風(fēng)化裂隙水的補(bǔ)給來源主要為大氣降水入滲，其次為上層松散巖類孔隙水垂向補(bǔ)給；其排泄主要為蒸發(fā)、側(cè)向排入溪流；其徑流方向是沿著地勢(shì)流向較低處的含水層或地表水系。深層脈狀裂隙水埋度變化不一，主要受構(gòu)造運(yùn)動(dòng)和結(jié)構(gòu)面發(fā)育的影響，與地表水的水力聯(lián)系相對(duì)較小。

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

3.1 ?水文地質(zhì)概念模型

洞庫地下工程主要包括地下水封洞庫及水幕系統(tǒng)2部分工程，水幕系統(tǒng)主要按平行主洞室軸線水平水幕孔的方式布置（圖3）。

結(jié)合水文地質(zhì)條件及工程垂向上的分布特征，將第四系松散巖類孔隙水及基巖裂隙水作為模擬含水層。考慮到工程垂向上的分布深度、巖性風(fēng)化及裂隙發(fā)育程度、現(xiàn)有壓水試驗(yàn)數(shù)據(jù)的支撐性，將模擬的含水層系統(tǒng)概化為非均質(zhì)、各向異性、三維非穩(wěn)定地下水流系統(tǒng)，模擬區(qū)面積18.32 km2。

側(cè)向邊界：模擬范圍以擬建庫區(qū)為中心，西側(cè)、東側(cè)以地表分水嶺為界，北側(cè)以北部山前花崗巖體出露線為界，南側(cè)以河流為界；東、西兩邊界垂直于地下水等水位線，均設(shè)定為零流量邊界，南、北兩側(cè)設(shè)定為流量邊界（圖4）。

垂向邊界：上邊界為潛水面，存在大氣降水入滲、地表水入滲補(bǔ)給及潛水蒸發(fā)排泄等垂向水量交換；模擬區(qū)底部均為新鮮基巖，概化為隔水底板。

3.2 ?地下水水流模型

3.2.1 ?數(shù)學(xué)模型

根據(jù)水文地質(zhì)概念模型，建立模擬區(qū)的數(shù)學(xué)模型。

（x，y，z）∈S_2，t>0））┤

式中：Ω表示地下水滲流區(qū)域；S1為模型的第一類邊界；S2為模型的第二類邊界；k_xx，k_yy，k_zz分別表示x，y，z主方向的滲透系數(shù)（m·d-1）；w表示源匯項(xiàng)，包括降水入滲補(bǔ)給、蒸發(fā)、井抽水量和泉排泄量（1/d）；μ_s第一層表示給水度與含水層厚度的比值，第二層至第十層表示彈性貯水率；H_0 （x，y，z）表示初始含水層水位標(biāo)高（m）；H_1 （x，y，z）為第一類邊界已知的含水層水位標(biāo)高（m）；q（x，y，z，t）為第二類邊界單位寬度流量（m2·d-1）。

3.2.2 ?模型設(shè)計(jì)

1）網(wǎng)格剖分

根據(jù)洞室、水幕系統(tǒng)結(jié)構(gòu)特點(diǎn)，本次模擬將洞庫所在區(qū)加密為10 m × 10 m矩形網(wǎng)格，外部區(qū)域剖分為30 m × 30 m矩形網(wǎng)格。網(wǎng)格剖分見表1和圖5。

2）源匯項(xiàng)處理

評(píng)價(jià)區(qū)內(nèi)地下水補(bǔ)給項(xiàng)主要為降水入滲及地下水側(cè)向流入，地下水排泄項(xiàng)主要為地下水側(cè)向流出、蒸發(fā)排泄及溝渠排泄。

降雨入滲補(bǔ)給：根據(jù)模擬區(qū)各月份平均降雨量數(shù)據(jù)，結(jié)合降雨入滲系數(shù)，換算成降水入滲補(bǔ)給強(qiáng)度給定到第一層的每個(gè)單元格當(dāng)中。降水入滲分區(qū)見圖6，其中：1、3、4、5區(qū)降水入滲系數(shù)為0.2，2區(qū)為0.1，6、7區(qū)為0.25。

潛水蒸發(fā)：根據(jù)模擬區(qū)內(nèi)多年平均蒸發(fā)量，利用軟件自帶蒸發(fā)程序包，按照水面蒸發(fā)強(qiáng)度、極限蒸發(fā)深度、地面標(biāo)高三者之間的關(guān)系，自動(dòng)進(jìn)行計(jì)算。極限蒸發(fā)深度設(shè)定為4 m，水面蒸發(fā)強(qiáng)度為6.62×10-4 m·d-1。

溝渠：對(duì)于模擬區(qū)內(nèi)的小型支溝、溪流，本次模擬采用DRAIN程序包進(jìn)行給定，各溝渠底部標(biāo)高參照溝渠與地下水的排泄關(guān)系，并最終通過調(diào)參進(jìn)行給定。

側(cè)向邊界：流量邊界根據(jù)該位置處的水力梯度、含水層厚度、滲透系數(shù)，利用達(dá)西定律估算進(jìn)行給定；水頭邊界根據(jù)邊界處的水位標(biāo)高進(jìn)行給定。

3）初始流場(chǎng)

本次模型的初始流場(chǎng)是依據(jù)2019年11月調(diào)查工作實(shí)測(cè)水位標(biāo)高及工程勘查鉆孔水位標(biāo)高進(jìn)行刻畫，最終給定的地下水水位初始流場(chǎng)見圖4。

4）含水層參數(shù)

本次含水層參數(shù)分區(qū)，以區(qū)域地形地貌、地質(zhì)構(gòu)造為基礎(chǔ)，結(jié)合庫區(qū)鉆探工作中的節(jié)理裂隙發(fā)育、巖心、提水試驗(yàn)、壓水試驗(yàn)相關(guān)成果，保守給定各層位含水層橫向滲透系數(shù)（Kxx）、垂向滲透系數(shù)（Kyy）、縱向滲透系數(shù)（Kzz）、給水度、彈性貯水率參數(shù)初始值。第一層為潛水含水層，給定了給水度；其他層位為微承壓水，給定了彈性貯水率。各層取值見表2。

3.2.3 ?模型的校正與檢驗(yàn)

本次模擬計(jì)算步長(zhǎng)為1 d，1年內(nèi)應(yīng)力期為90 d，1 年后間隔為1 a。計(jì)算水位數(shù)據(jù)與2020年4月（識(shí)別驗(yàn)證期）實(shí)測(cè)的水位數(shù)據(jù)，從區(qū)域流場(chǎng)形態(tài)、鉆孔水位擬合情況、模型地下水均衡項(xiàng)3個(gè)方面進(jìn)行檢驗(yàn)。從擬合結(jié)果來看，地下水流場(chǎng)基本符合區(qū)域地下水滲流趨勢(shì)，流場(chǎng)擬合結(jié)果較好。同時(shí)根據(jù)模型地下水均衡結(jié)果（表3）可以看出，模擬區(qū)識(shí)別期內(nèi)地下水總體呈負(fù)均衡狀態(tài)，與區(qū)域地下水水位變動(dòng)情況基本相符，因此，本模型的可靠性是有保障的，可進(jìn)行下一步流場(chǎng)及溶質(zhì)預(yù)測(cè)。區(qū)域地下水水位擬合結(jié)果及鉆孔水位擬合結(jié)果見圖7～9。

3.3 ?地下水溶質(zhì)運(yùn)移模型

本次污染物模擬預(yù)測(cè)過程不考慮污染物在含水層中的吸附、揮發(fā)、生物化學(xué)反應(yīng)，模型中各項(xiàng)參數(shù)予以保守性考慮。這樣選擇的理由有3點(diǎn)：1）有機(jī)污染物在地下水中的運(yùn)移非常復(fù)雜，影響因素除對(duì)流、彌散作用以外，還存在物理、化學(xué)、微生物等作用，這些作用常常會(huì)使污染物總量減少，運(yùn)移擴(kuò)散速度減慢。目前國(guó)際上對(duì)這些作用參數(shù)的準(zhǔn)確獲取還存在一定困難；2）從保守性角度考慮，假設(shè)污染物在運(yùn)移中不與含水層介質(zhì)發(fā)生反應(yīng)，可以被認(rèn)為是保守型污染物，只按保守型污染物計(jì)算，即只考慮運(yùn)移過程中的對(duì)流、彌散作用。在國(guó)際上有很多用保守型污染物作為模擬因子的環(huán)境影響評(píng)價(jià)成功實(shí)例；3）保守型考慮符合環(huán)境影響評(píng)價(jià)風(fēng)險(xiǎn)最大的原則。

3.3.1 ?數(shù)學(xué)模型

地下水中溶質(zhì)運(yùn)移的數(shù)學(xué)模型可表示為：

n_e ??C/?t=?/（?X_i ） （nD_ij ??C/（?X_j ））-?/（?x_i ） （nCV_i ）±C^' ?W

C（x，y，0）=C_0 （x，y）（x，y）∈Ω，t=0

（C ν ?-D_gradc）?n ?|┤_（Γ_2 ）=φ（x，y，t） （x，y）∈Γ_2，t≥0

式中：D_ij=α_ijmn ?（V_m V_n）/|V| ；αijmn為含水層的彌散度；Vm，Vn分別為m和n方向上的速度分量；|V|為速度模；C為模擬污染質(zhì)的濃度（mg·L-1）；t為時(shí)間（d）；ne為有效孔隙度；n為介質(zhì)孔隙度；W為源匯單位面積上的通量；Vi為滲流速度（m·d-1）；C＇為源匯的污染質(zhì)濃度（mg·L-1）；C0（x，y，z）為已知濃度分布，Ω為模型模擬區(qū)；Γ2為通量邊界，Dgradc為濃度梯度。

3.3.2 ?源匯項(xiàng)和邊界條件

模擬區(qū)內(nèi)的自然條件相對(duì)穩(wěn)定，主要表現(xiàn)在降雨量、蒸發(fā)量等氣象要素年際變化不大，模擬區(qū)內(nèi)地下水未來開采量變化不大，可近似等于現(xiàn)狀開采量。由于洞庫布置在水幕系統(tǒng)之下，因此按照水幕位置及底板高程確定水位高程。

本次模型將污染源以面源形式設(shè)定濃度邊界，污染源位置按實(shí)際設(shè)計(jì)概化。在模擬污染物擴(kuò)散時(shí)，不考慮吸附作用、化學(xué)反應(yīng)等因素，重點(diǎn)考慮對(duì)流、彌散作用。

為了分析運(yùn)營(yíng)期洞庫內(nèi)石油類隨地下水運(yùn)移對(duì)周邊地下水環(huán)境造成的影響，利用經(jīng)過校正檢驗(yàn)的水流模型進(jìn)行溶質(zhì)運(yùn)移預(yù)測(cè)。

3.3.3 ?彌散度的確定

由于水動(dòng)力彌散尺度效應(yīng)的存在，難以通過野外或室內(nèi)彌散試驗(yàn)獲得真實(shí)的彌散度。因此，本次參考蘇玉娟等（2019）裂隙承壓含水層彌散試驗(yàn)的研究成果，取縱向彌散度值為4.8×10-3 m，橫向彌散度值為4.4×10-4 m。

4 ?預(yù)測(cè)情景確定與結(jié)果分析

4.1 ?建設(shè)期地下水環(huán)境影響預(yù)測(cè)與評(píng)價(jià)

1）建設(shè)期預(yù)測(cè)情景確定及賦值

根據(jù)工程內(nèi)容分析，地下水封洞庫區(qū)地下工程施工建設(shè)主要是對(duì)地下水水位產(chǎn)生影響，對(duì)地下水水質(zhì)的影響較小。洞室開挖時(shí)水平水幕系統(tǒng)已經(jīng)投入使用，地下水模型將水平水幕概化為定水頭邊界，其位置、厚度和范圍按水幕系統(tǒng)實(shí)際設(shè)計(jì)概化，洞室層按頂板、底板平均標(biāo)高給定，水幕系統(tǒng)水頭值取設(shè)計(jì)水頭標(biāo)高。

2）建設(shè)期預(yù)測(cè)結(jié)果分析

施工期洞室涌水量逐年減小，第3年（1 095 d）涌水量為789 m3·d-1，第4年（1 460 d）為673 m3·d-1，趨于穩(wěn)定。施工期洞庫涌水量變化見圖10。

洞室開挖4年后，洞室紅線范圍內(nèi)潛水水位最大降深值為2.8 m，洞室紅線外潛水水位降深最大值為0.1 m，最大影響半徑為396 m。洞室開挖4年后洞室層地下水水位最大降深值為85 m，位于洞室正上方，最大影響半徑為210 m。洞室開挖4年后地下水水位下降值見圖11、圖12，與洞室開挖1年后相比，洞室開挖4 年后地下水漏斗有所擴(kuò)展。

4.2 ?運(yùn)營(yíng)期地下水環(huán)境影響預(yù)測(cè)與評(píng)價(jià)

1）運(yùn)營(yíng)期預(yù)測(cè)情景確定及賦值

在運(yùn)營(yíng)期，考慮地下洞室和水幕系統(tǒng)條件，其形狀和深度等按實(shí)際設(shè)計(jì)概化，水幕設(shè)置為定水頭邊界，覆蓋整個(gè)洞庫上方范圍。洞室層按運(yùn)行壓力給定水位標(biāo)高值，水幕系統(tǒng)頂板水頭值取設(shè)計(jì)標(biāo)高值。洞室設(shè)為定濃度邊界，預(yù)測(cè)石油類持續(xù)擴(kuò)散對(duì)地下水水質(zhì)的影響。

2）運(yùn)營(yíng)期預(yù)測(cè)結(jié)果分析

運(yùn)營(yíng)期前3年涌水量變化較大，從第4年至運(yùn)營(yíng)期結(jié)束涌水量較穩(wěn)定，約為420 m3·d-1。運(yùn)營(yíng)期洞庫涌水量變化見圖13。

洞室運(yùn)行30年后影響范圍趨于穩(wěn)定，達(dá)到最大范圍。洞室紅線范圍內(nèi)潛水水位最大降深值為20 m，降深值大于20 m的范圍占整個(gè)洞庫面積的0.1%；洞室紅線外潛水水位降深最大值為5 m，最大影響半徑為883 m。洞室層地下水水位最大降深值100 m，位于洞室正上方，下降后預(yù)測(cè)水位比水幕系統(tǒng)設(shè)計(jì)水位高3 m，最大影響半徑600 m。水位影響范圍見圖14。

在運(yùn)營(yíng)期內(nèi)，石油類污染范圍局限于洞室周圍19 m內(nèi)，其影響范圍很小。由于洞庫為排水邊界，洞內(nèi)石油不會(huì)向外泄露，而且?guī)r層滲透系數(shù)和彌散系數(shù)都很小，洞庫對(duì)地下水水質(zhì)影響局限在建筑界線附近。水質(zhì)影響范圍見圖15。

5 ?結(jié)論

本次利用地下水?dāng)?shù)值模擬有效刻畫出了水文地質(zhì)模型，反映了洞庫施工期和運(yùn)營(yíng)期地下水滲流場(chǎng)變化情況，提前預(yù)判了洞庫涌水量和對(duì)周圍地下水環(huán)境的影響，了解了洞庫的水封效果，為洞庫的設(shè)計(jì)和建設(shè)提供必要的技術(shù)支持。

1）在有水幕條件下：洞庫上方潛水水位最大降深值為20 m，洞庫四周最大降深值為5 m；洞室層地下水水位最大降深值為100 m，下降后水位比水幕系統(tǒng)的設(shè)計(jì)水位高3 m，能夠保證洞庫的水封條件，減小了原油泄漏的風(fēng)險(xiǎn)，對(duì)地下水環(huán)境影響較小。

2）洞庫建設(shè)引起地下水位的下降主要影響山體中上部基巖裂隙水，可能會(huì)導(dǎo)致部分泉流量減少甚至干枯，但影響范圍有限，對(duì)山腰以下松散巖的孔隙潛水基本不產(chǎn)生影響。

3）由于水幕系統(tǒng)補(bǔ)水，洞內(nèi)石油不會(huì)向外泄露，而且?guī)r層滲透系數(shù)和彌散系數(shù)都很小，洞庫對(duì)地下水水質(zhì)影響局限在建筑界線附近。

參考文獻(xiàn)

蔡紅飚，陳情來，周亮臣，2000.LPG地下水封洞庫工程地質(zhì)條件分析［C］//第六屆全國(guó)工程地質(zhì)大會(huì)論文集：504-506.

郭書太，高劍鋒，陳雪見，倪亮，代云清，2006.錦州某地下水封洞庫工程地質(zhì)條件適宜性分析［C］//第二屆全國(guó)巖土與工程學(xué)術(shù)大會(huì)論文集：277-282.

韓曼，2007.地下水封石油洞庫滲流場(chǎng)及溶質(zhì)運(yùn)移模擬研究［D］.青島：中國(guó)海洋大學(xué).

洪淑娜，2023.水封洞庫石油類污染物地下水時(shí)空分布模擬［J］.廣東化工，50（4）： 163- 165.

胡成，陳剛，曹孟雄，唐連松，鄭可，王繼剛，2022.基于離散裂隙網(wǎng)絡(luò)法和水流數(shù)值模擬技術(shù)的地下水封洞庫水封性研究［J］．地質(zhì)科技通報(bào)，41（1）：119-126.

李術(shù)才，張立，馬秀媛，薛翊國(guó)，王者超，李逸凡，平洋，姜彥彥，2013.大型地下水封石油洞庫滲流場(chǎng)時(shí)空演化特征研究［J］.巖土力學(xué)，34（7）：1 979-1 986.

李印，陳雪見，姜俊浩， 2022.大型地下水封洞庫運(yùn)營(yíng)期涌水量計(jì)算與注漿設(shè)計(jì)［J］.油氣儲(chǔ)運(yùn)，41（7）：847-851.

梁斌，陳剛，胡成，2017.某地下水封油庫地下水環(huán)境影響評(píng)價(jià)［C］//2017年中國(guó)環(huán)境科學(xué)學(xué)會(huì)科學(xué)與技術(shù)年會(huì)論文集（第二卷）：2 137-2 144.

彭振華，李俊彥，2011.某地下水封洞庫工程地質(zhì)特性研究［J］.勘察科學(xué)技術(shù)，29（1）：25-27.

時(shí)洪斌，2010.黃島地下水封洞庫水封條件和圍巖穩(wěn)定性分析與評(píng)價(jià)［D］.北京交通大學(xué).

蘇玉娟，陶勇，2019.花崗巖風(fēng)化殼網(wǎng)狀裂隙承壓含水層彌散試驗(yàn)研究［J］.勘察科學(xué)技術(shù)，37（6）：19-23.

王婷婷，2019.LPG地下水封洞庫的建設(shè)［J］.化工設(shè)計(jì)通訊，45（1）：236+250.

王者超，張彬，喬麗蘋，楊森，2022.中國(guó)地下水封儲(chǔ)存理論與關(guān)鍵技術(shù)研究進(jìn)展［J］.油氣儲(chǔ)運(yùn)，41（9）：995- 1 003.

夏喜林，劉燁， 2004.淺談我國(guó)地下油庫的建設(shè)［J］. 石油規(guī)劃設(shè)計(jì)， 15（4） ：26-27.

收稿日期：2023-10-23；修回日期：2023-12-22

第一作者簡(jiǎn)介：李雅靜（1987- ），女，碩士，工程師，主要從事水文地質(zhì)與環(huán)境地質(zhì)工作。E-mail：616225682@qq.com

引用格式：李雅靜，郭迎濤，劉立才，張訓(xùn)玉，2024.某地下水封石油洞庫地下水滲流場(chǎng)及涌水量模擬［J］.城市地質(zhì)，19（1）：70-79