劉躍蔭

摘要：Visual MODFLOW4.1是由加拿大Watedoo水文地質(zhì)公司在MODFLOW的基礎(chǔ)上開(kāi)發(fā)研制而成，曾是國(guó)際上最為流行且被各國(guó)同行一致認(rèn)可的三維地下水流和溶質(zhì)運(yùn)移模擬評(píng)價(jià)的標(biāo)準(zhǔn)可視化專(zhuān)業(yè)軟件系統(tǒng)。該軟件提供了求解地下水流有限差分公式的很多種方法，如強(qiáng)隱式迭代法SIP、逐次超松弛迭代法SOR、預(yù)調(diào)共扼梯度迭代法PCGZ、SSOR等，使模擬計(jì)算更加貼近實(shí)際，計(jì)算精度進(jìn)一步提高。

關(guān)鍵詞：數(shù)值法；巖溶地區(qū)；地下水；可開(kāi)采資源量

數(shù)學(xué)模型是建立在對(duì)天然地質(zhì)體概念模型高度認(rèn)識(shí)的基礎(chǔ)上進(jìn)行的，要確定一個(gè)地下水流問(wèn)題的數(shù)學(xué)模型，就必須進(jìn)行水文地質(zhì)測(cè)繪、鉆探、試驗(yàn)，查明地質(zhì)、水文地質(zhì)條件。為了更精確地評(píng)價(jià)研究區(qū)的地下水資源量，預(yù)測(cè)地下水系統(tǒng)狀態(tài)及其變化趨勢(shì)。本次根據(jù)柳東區(qū)的水文地質(zhì)條件、鉆孔、地下水水位動(dòng)態(tài)，建立地下水系統(tǒng)模擬模型；在地層巖性、地質(zhì)構(gòu)造和地形地貌條件的控制下，柳東區(qū)的水文地質(zhì)邊界清晰，屬于一個(gè)較完整的水文地質(zhì)單元，模擬范圍以河間地塊地下水分水嶺與柳江為界。

1.地下水賦存介質(zhì)概化

柳東區(qū)地形為峰林谷地一洼地地貌，地面高程一般80～380m，北面及南面有碳酸鹽巖山體出露，其余皆為第四系所覆蓋。模型三維立體展示如圖1所示。

本次地下水資源數(shù)值法評(píng)價(jià)采用地下水模擬軟件VisualMODFLOW4.1進(jìn)行模擬。

Visual MODFLOW4.1是由加拿大Waterloo水文地質(zhì)公司在MODFLOW的基礎(chǔ)上開(kāi)發(fā)研制而成，曾是國(guó)際上最為流行且被各國(guó)同行一致認(rèn)可的三維地下水流和溶質(zhì)運(yùn)移模擬評(píng)價(jià)的標(biāo)準(zhǔn)可視化專(zhuān)業(yè)軟件系統(tǒng)。該軟件提供了求解地下水流有限差分公式的很多種方法，如強(qiáng)隱式迭代法SIP、逐次超松弛迭代法SOR、預(yù)調(diào)共扼梯度迭代法PCGZ、SSOR等，使模擬計(jì)算更加貼近實(shí)際，計(jì)算精度進(jìn)一步提高。

如圖，第四系覆蓋層之下地層有孤峰組下段（P2g1）、棲霞組（P2g）、南丹組（C2Pn）、馬平組（C2Pm）和黃龍組下段（c2h2）灰?guī)r以及黃龍組上段（Ceh1）、大埔組（C2d）白云巖，受巖性及構(gòu)造影響，該區(qū)巖溶裂隙、溶洞較發(fā)育，為地下水賦存提供有利條件。根據(jù)模擬區(qū)地下水系統(tǒng)含水介質(zhì)的特點(diǎn)及水文地質(zhì)條件，將含水層概化為兩層。

2含水層水力特征的概化

數(shù)學(xué)模型是建立在對(duì)天然地質(zhì)體概念模型高度認(rèn)識(shí)的基礎(chǔ)上進(jìn)行的，要確定一個(gè)地下水流問(wèn)題的數(shù)學(xué)模型，就必須進(jìn)行水文地質(zhì)測(cè)繪、鉆探、試驗(yàn)，查明地質(zhì)、水文地質(zhì)條件。為了更精確地評(píng)價(jià)研究區(qū)的地下水資源量，預(yù)測(cè)地下水系統(tǒng)狀態(tài)及其變化趨勢(shì)。本次根據(jù)柳東區(qū)的水文地質(zhì)條件、鉆孔、地下水水位動(dòng)態(tài)，建立地下水系統(tǒng)模擬模型；在地層巖性、地質(zhì)構(gòu)造和地形地貌條件的控制下，柳東區(qū)的水文地質(zhì)邊界清晰，屬于一個(gè)較完整的水文地質(zhì)單元，模擬范圍以河間地塊地下水分水嶺與柳江為界。從空間上看，地下水流整體上以水平運(yùn)動(dòng)為主、垂向運(yùn)動(dòng)為輔，地下水系統(tǒng)的流入、輸出隨著時(shí)間、空間變化而變化，故該系統(tǒng)地下水為非穩(wěn)定流；參數(shù)隨時(shí)間和空間變化，體現(xiàn)了系統(tǒng)的非均質(zhì)性，在水平與垂向上有明顯的差異性，所以參數(shù)概化成各向異性。故將模擬區(qū)的地下水流作為非均質(zhì)各向異性三維非穩(wěn)定流處理。

3.柳東區(qū)邊界條件的概化

柳東區(qū)地形為峰林谷地一洼地地貌，地面高程一般80～380m，北面及南面有碳酸鹽巖山體出露，其余皆為第四系所覆蓋。上覆第四系粘土可視為相對(duì)隔水層，下伏碳酸鹽巖溶水含水層，稍具承壓性，以下對(duì)各個(gè)邊界概化作具體闡述。

1）側(cè)向邊界概化

①東西兩側(cè)為柳江河，河床切割較深，均切割基巖，河水流量較穩(wěn)定，水流平緩，水位變幅小，可作為碳酸鹽巖溶水含水層第一類(lèi)邊界，即河流定水頭邊界。

②北面牛車(chē)坪村一馬鹿村一帶，山體高程達(dá)180～380m，為地下水分水嶺，主要為碎屑巖分布區(qū)，巖性以砂、泥及頁(yè)巖為主，透性及富水性差，可作為隔水邊界，即為第二類(lèi)邊界。

③南面以架鶴村一市福利院一帶的局部地下水分水嶺作為隔水邊界，即為第二類(lèi)邊界。

4.水文地質(zhì)參數(shù)分區(qū)和初值的確定

模型主要的水文地質(zhì)參數(shù)為滲透系數(shù)、給水度、貯水率。在不影響計(jì)算精度的情況下，為了使參數(shù)分區(qū)簡(jiǎn)化，將每一個(gè)滲透系數(shù)分區(qū)同時(shí)視為給水度、貯水率分區(qū)。將第四系土層分區(qū)劃分為一層，根據(jù)模擬區(qū)的水文地質(zhì)條件，將碳酸鹽巖含水層在平面上進(jìn)行參數(shù)初步分區(qū)，其參數(shù)大小的初始值按照本次各地段施工鉆孔及收集以往資料的鉆孔揭露的巖溶發(fā)育特征和相應(yīng)有鉆孔抽水實(shí)驗(yàn)資料，初步確定每個(gè)分區(qū)相對(duì)應(yīng)的參數(shù)初值。

5.滲流場(chǎng)數(shù)值模型的識(shí)別和驗(yàn)證

取SZ10、SZ20、SZ23鉆孔觀測(cè)水位數(shù)據(jù)與模擬水位數(shù)據(jù)進(jìn)行對(duì)比，作為模型的識(shí)別和檢驗(yàn)，根據(jù)收集到的地下水觀測(cè)資料，取模型識(shí)別的時(shí)段為2015年5月20日到2016年4月20日，以實(shí)際水位和計(jì)算水位的擬合情況作為調(diào)整參數(shù)的依據(jù)。通過(guò)多次修正滲透系數(shù)、給水度、貯水率等，得出了含水層滲透系數(shù)和貯水率的大小及分區(qū)與觀測(cè)孔的計(jì)算與實(shí)測(cè)水位的對(duì)比。

計(jì)算水位較實(shí)測(cè)水位滯后，模擬水位比較接近實(shí)測(cè)水位，水位線(xiàn)趨勢(shì)相同，水位擬合誤差小于0.5m。水位變化值較少的情況下誤差在容許范圍內(nèi)。從圖中可看出，各觀測(cè)孔模擬值與實(shí)際值動(dòng)態(tài)擬合效果較好，說(shuō)明模擬所用的參數(shù)是合理的。

6.模擬區(qū)地下水流場(chǎng)分析

通過(guò)數(shù)學(xué)模型識(shí)別和驗(yàn)證，表明含水層結(jié)構(gòu)、邊界條件概化、水文地質(zhì)參數(shù)的選取較合理的，數(shù)學(xué)模型能較真實(shí)地反映研究區(qū)地下水系統(tǒng)特征，可用模型數(shù)據(jù)進(jìn)行地下水流場(chǎng)分析。天然條件下河間地塊地下水均向柳江排泄，在中部形成了一條近南北向的地下水分水嶺，將模擬區(qū)分成兩個(gè)次級(jí)水文地質(zhì)單元。東面水力坡度比西面的小。隨著地下水水位的下降，水力坡度逐漸減小。整體地下水位略高于柳江水位，水流較平緩。

7.模擬區(qū)地下水允許開(kāi)采量的預(yù)測(cè)

1）預(yù)測(cè)方案

在降水保證率95%條件下，利用建立的地下水模型對(duì)設(shè)計(jì)方案進(jìn)行模擬預(yù)測(cè)，本次預(yù)測(cè)模型的邊界條件和水文地質(zhì)參數(shù)保持不變，補(bǔ)給項(xiàng)和排泄項(xiàng)按95%保證率重新計(jì)算，開(kāi)采量擬采用東環(huán)水源地開(kāi)采井來(lái)設(shè)計(jì)開(kāi)采量，通過(guò)多次調(diào)試來(lái)確定。預(yù)測(cè)期為2016～2026年，以月、年為時(shí)間段，共21個(gè)應(yīng)力期，第一年以—個(gè)月為—個(gè)應(yīng)力期，第一年以后以一年為—個(gè)應(yīng)力期，每個(gè)時(shí)間段內(nèi)包括不同的時(shí)間步長(zhǎng)。初始流場(chǎng)選擇2015年12月末的流場(chǎng)。

2）預(yù)測(cè)結(jié)果

通過(guò)調(diào)試最終以開(kāi)采量48000 m3/d進(jìn)行預(yù)測(cè)，按設(shè)計(jì)方案開(kāi)采3個(gè)月后，水源地中心地下水位大幅下降形成降落漏斗，開(kāi)采1年、3年后，地下水位整體持續(xù)下降，漏斗面積不斷增大，直至6年后，地下水流場(chǎng)達(dá)到穩(wěn)定狀態(tài)，水位基本保持不變，與10年后相比降幅不超過(guò)O.02m，漏斗面積不再增加。因此，預(yù)測(cè)該模擬區(qū)允許開(kāi)采量為48000 m3/d（17.52×106m3/a）。

結(jié)語(yǔ)：在礦坑開(kāi)采過(guò)程中，開(kāi)采工作需要符合國(guó)家的技術(shù)規(guī)范以及要求，勘測(cè)礦山的地址水文條件以及礦山地址水文的類(lèi)型復(fù)雜程度，結(jié)合實(shí)際環(huán)境設(shè)計(jì)出符合條件的工作任務(wù)，并確定大致的工期范圍。在施工開(kāi)采的過(guò)程中，專(zhuān)業(yè)人員應(yīng)時(shí)時(shí)檢測(cè)地下水的流量變化，以防止意外事件發(fā)生，及時(shí)完成工作任務(wù)，提升安全防范意識(shí)。endprint