范紅晨，孫曉梁，杜 堯*，鄧婭敏，劉廣寧

(1.中國地質(zhì)大學(xué)(武漢)環(huán)境學(xué)院，湖北 武漢 430078；2.中國地質(zhì)大學(xué)(武漢)生物地質(zhì)與環(huán)境地質(zhì)國家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430078；3.中國地質(zhì)大學(xué)(武漢)地質(zhì)調(diào)查研究院，湖北 武漢 430074；4.中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心，湖北 武漢 430205)

地下水是影響湖泊水量和水質(zhì)的重要端元，地下水向湖泊排泄(即湖底地下水排泄)的量化對(duì)于評(píng)估地下水在湖泊水量和物質(zhì)均衡中的貢獻(xiàn)具有重要作用。不同的技術(shù)方法被用于湖底地下水排泄的量化，包括滲流儀測(cè)量、水量均衡模型、惰性示蹤劑質(zhì)量平衡模型、放射性示蹤劑質(zhì)量平衡模型、溫度示蹤等，其中放射性示蹤劑氡(Rn)質(zhì)量平衡模型是量化湖底地下水排泄最為普遍而有效的方法之一。

在Rn質(zhì)量平衡模型中，地下水中Rn活度是最重要的端元之一，也是基于Rn質(zhì)量平衡模型估算湖底地下水排泄通量時(shí)最主要的不確定性來源，這是因?yàn)閷⑹聚檮┩哭D(zhuǎn)換為水文通量時(shí)需要除以地下水端元的示蹤劑活度。一般來說，地下水端元的Rn活度有3種選擇，分別是湖區(qū)周邊民井或監(jiān)測(cè)井地下水中Rn活度、湖岸區(qū)孔隙水中Rn活度和沉積物平衡試驗(yàn)所得孔隙水中Rn活度。不同地下水端元的Rn活度可能存在較大的差異，因此地下水端元的選取對(duì)于Rn質(zhì)量平衡模型結(jié)果的輸出以及湖底地下水排泄通量的估算至關(guān)重要。

Burnett等在研究柬埔寨洞里薩湖的湖底地下水排泄通量時(shí)選取湖區(qū)周邊井水中平均Rn活度作為地下水端元值；Dabrowski等采用井水中平均Rn活度作為地下水端元值，估算了阿拉斯加育空-庫斯科維姆三角洲一個(gè)淺水苔原湖的湖底地下水排泄通量；Dimova等利用Rn質(zhì)量平衡模型評(píng)價(jià)佛羅里達(dá)州中部和北部湖泊群的湖底地下水排泄時(shí)，分別采用湖岸孔隙水中Rn活度和沉積物培養(yǎng)試驗(yàn)所得孔隙水中Rn活度作為地下水端元值，估算了Haines和Josephine湖泊的湖底地下水排泄通量，結(jié)果發(fā)現(xiàn)估算值之間存在顯著的差別。然而，在基于質(zhì)量平衡模型評(píng)價(jià)地下水排泄的研究領(lǐng)域中，對(duì)于地下水端元選取的適宜性一直沒有定論，且對(duì)于不同地下水端元選取對(duì)Rn質(zhì)量平衡模型量化湖底地下水排泄的影響研究十分薄弱。

針對(duì)這一問題，本文以長江中游典型牛軛湖——天鵝洲濕地為研究區(qū)，以Rn質(zhì)量平衡模型為主要研究方法，分別選取湖區(qū)周邊井水、湖岸孔隙水、沉積物平衡試驗(yàn)所得孔隙水這3個(gè)端元或端元組合的Rn活度作為地下水端元值，對(duì)地下水向天鵝洲濕地的排泄通量進(jìn)行了估算，評(píng)估了不同地下水端元或端元組合的選取對(duì)Rn質(zhì)量平衡模型量化湖底地下水排泄的影響。本研究可為更精確地量化地下水向湖泊的排泄通量提供理論與技術(shù)支撐。

1 研究區(qū)概況

天鵝洲濕地位于長江中游，是長江主河道于1972年截彎取直后形成的牛軛湖，見圖1。區(qū)內(nèi)有麋鹿自然保護(hù)區(qū)和白鰭豚自然保護(hù)區(qū)兩個(gè)國家級(jí)自然保護(hù)區(qū)，同時(shí)棲息著大量水生生物和陸生生物，其中包括各種珍稀水禽。天鵝洲濕地長約21 km，豐水期總水面面積約20 km，蓄水量為1×10m，最大深度為20 m，平均寬度為1 200 m。為了控制洪水，1998年在長江主河道和天鵝洲之間修建了沙灘子大堤，使天鵝洲濕地與長江主河道之間的水力聯(lián)系被切斷。目前僅在天鵝洲濕地東北和東南建有兩處水閘(即馮家潭閘和天鵝洲閘)與外界連通，其中在枯水期時(shí)閘門關(guān)閉。

天鵝洲濕地區(qū)周邊地層為第四系松散沉積物，其中表層普遍分布有厚度1～15 m不等的黏性土層，巖性主要為黏土和亞黏土，局部含粉砂，富水性較弱。而從黏性土層往下，巖性主要為粉砂和細(xì)砂，富水性較強(qiáng)，為區(qū)內(nèi)主要的含水層位(即孔隙承壓含水層)，也是湖區(qū)周邊民井主要的取水層位。天鵝洲濕地與這一主要含水層位直接連通，并與其具有良好的水力聯(lián)系(見圖1)。

圖1 研究區(qū)位置(a)、典型水文地質(zhì)剖面(b)和采樣點(diǎn) 分布圖(c)Fig.1 Location(a), typical hydrogeological section(b) and sampling point distribution(c) of the study area

2 野外采樣及研究方法

2.1 野外采樣工作

于枯水期(2019年12月)在天鵝洲濕地區(qū)進(jìn)行了野外采樣工作，共采集49個(gè)樣品，其中湖水樣26個(gè)，井水樣11個(gè)，湖岸孔隙水樣7個(gè)，沉積物樣5個(gè)(見圖1)。利用差分GPS技術(shù)對(duì)湖水(地表水)和井水(地下水)的水位進(jìn)行了測(cè)量，發(fā)現(xiàn)井水水位顯著高于湖水水位(見圖2)，說明在枯水期天鵝洲濕地與地下水相互作用模式為地下水向湖泊排泄。

圖2 湖水(地表水)水位和井水(地下水)水位對(duì)比圖Fig.2 Comparison between the lake water surface water level and well water(groundwater) level

湖水樣是從湖水面以下0.5 m處直接采集；井水樣是通過手壓井或蠕動(dòng)泵采集，在抽水至少10 min后開始采集井水樣，以保證水樣的新鮮度；采集湖岸孔隙水樣時(shí)，將測(cè)壓管插入地表以下1 m左右，測(cè)壓管末端0～20 cm有孔徑3 mm左右的小孔，以保證孔隙水可以進(jìn)入測(cè)壓管，然后用手持式蠕動(dòng)泵抽取測(cè)壓管中的孔隙水；采取沉積物樣時(shí)，將PVC管插入湖底，然后將管中沉積物取出，刮去外層與湖水混合的部分，保留內(nèi)部沉積物。在野外用RAD7 HO儀測(cè)定了湖水、井水、湖岸孔隙水以及大氣中Rn活度。通過中國氣象網(wǎng)(http://data.cma.cn)獲取了野外工作期間逐小時(shí)的風(fēng)速數(shù)據(jù)。

2.2 222Rn質(zhì)量平衡模型

Rn是Ra的子體，其單質(zhì)形態(tài)是氡氣，半衰期為3.823 d。對(duì)湖水中氡的源匯關(guān)系建立Rn質(zhì)量平衡模型，將地下水排泄通量作為唯一的未知項(xiàng)進(jìn)行求取(見圖3)。

圖3 222Rn質(zhì)量平衡模型的概念圖解Fig.3 Conceptual model of the 222Rn mass balance

對(duì)天鵝洲濕地而言，由于沒有地表水的流入和流出，湖水中Rn通量的源項(xiàng)主要包括地下水排泄、沉積物擴(kuò)散，而湖水中Rn通量的匯項(xiàng)包括大氣擴(kuò)散、自身衰變，故Rn質(zhì)量平衡模型可表示如下:

(1)

式中：

F

、

F

、

F

分別表示地下水排泄、沉積物擴(kuò)散、大氣擴(kuò)散的Rn通量[Bq/(m·d)]；

I

Rn表示湖水中Rn的儲(chǔ)量(Bq/m)，其值等于湖水中Rn的活度(Bq/m)乘以湖水深度(m)；

λ

Rn表示Rn的衰變常數(shù)(d)，其值為0.186 d；?

I

Rn/?

t

表示湖水中Rn儲(chǔ)量隨時(shí)間的變化，其值為零。

地下水排泄速率通過如下公式計(jì)算：

(2)

式中：

V

表示地下水排泄速率(mm/d)；

F

表示地下水排泄的Rn通量[Bq/(m·d)]；

C

表示地下水端元的Rn活度(Bq/m)。

2.2.1Rn的大氣擴(kuò)散

一般來說，湖水中氡活度高于大氣中氡活度。在本研究中，湖水和大氣的Rn平均活度分別為333.69 Bq/m和16.2 Bq/m。因此，由于濃度梯度大，Rn會(huì)從水中逃逸到大氣中，相應(yīng)的Rn大氣損失(

F

)可用如下公式計(jì)算：

F

=

k

×

(C

-

αC

)

(3)

式中

:C

和

C

分別為湖水和上覆大氣中Rn活度(Bq/m)；

α

為溶解分配系數(shù)，由公式(4)計(jì)算；

k

為氡氣的氣體傳遞系數(shù)(m/d)，是確定

F

的關(guān)鍵因子。

α

=0.105+0.405e-0.050 2

(4)

式中：

T

為湖水的溫度(℃)。

(5)

式中：

μ

表示風(fēng)速(m/s)；

Sc

表示特定水溫下氡氣的施密特?cái)?shù)，其定義為運(yùn)動(dòng)黏滯系數(shù)(

v)

與分子擴(kuò)散系數(shù)

(D

)的比值。Pilson在1998年提出了根據(jù)湖水溫度

T

計(jì)算施密特?cái)?shù)的方法，其計(jì)算公式為

Sc

=3 417.6e-0.063 4×

(6)

式中的

Sc

除以600即為標(biāo)準(zhǔn)化到20℃時(shí)的施密特?cái)?shù)。

2.2.2Rn的沉積物擴(kuò)散

在湖水-沉積物界面，沉積物孔隙水中的氡濃度遠(yuǎn)高于湖水中的氡濃度，從而使Rn從湖底沉積物向湖水?dāng)U散。Rn的沉積物擴(kuò)散通量可表示為

在國內(nèi)外出現(xiàn)了一些新的項(xiàng)目合作態(tài)勢(shì)。即業(yè)主與承包商從斗智斗勇的競(jìng)爭(zhēng)者轉(zhuǎn)變?yōu)閼?zhàn)略合作者，雙方達(dá)成戰(zhàn)略聯(lián)盟，建立一種合作伙伴的關(guān)系。強(qiáng)調(diào)了企業(yè)之間的合作共贏、資源共享、風(fēng)險(xiǎn)共擔(dān)、利益均攤的關(guān)系。這種新的管理理念尤其強(qiáng)調(diào)了信息的交流與共享對(duì)于增強(qiáng)企業(yè)合作競(jìng)爭(zhēng)力的作用，有效的降低了信息不對(duì)稱的程度。

F

=(

λ

Rn×

D

)

(C

-

C

)

(7)

式中：

C

表示平衡時(shí)沉積物孔隙水中氡濃度(Bq/m)；

D

表示Rn從沉積物孔隙水向上覆水體中擴(kuò)散時(shí)的擴(kuò)散系數(shù)(m/d) ；

C

表示實(shí)際測(cè)得的上覆水體中Rn活度(Bq/m) 。

D

值大致等于沉積物孔隙度(

n

)乘以Rn的分子擴(kuò)散系數(shù)(

D

)，而在數(shù)值上分子擴(kuò)散系數(shù)(

D

)是溫度(

T

)的函數(shù)，其表達(dá)式如下：

D

=nD

(8)

(9)

C

值需要通過沉積物平衡培養(yǎng)試驗(yàn)得到，具體方法為將沉積物表層樣品取回實(shí)驗(yàn)室后，取約150 g沉積物和500 mL原位湖水置于錐形瓶內(nèi)并密封，放入搖床培養(yǎng)30 d直至沉積物孔隙水中Rn活度和其上覆湖水中Rn活度達(dá)到平衡，再將平衡后的湖水用溢流法轉(zhuǎn)移至250 mL采樣瓶內(nèi)，并用RAD7及RAD HO水中氡配件進(jìn)行測(cè)量。

C

值可以根據(jù)沉積物培養(yǎng)試驗(yàn)的結(jié)果按下式計(jì)算26：

(10)

式中：

C

為根據(jù)試驗(yàn)得到的沉積物濕樣中Rn濃度(Bq/kg)；

ρ

為沉積物的濕密度(kg/m)；

n

為沉積物的孔隙度。

3 結(jié)果與討論

3.1 湖水和地下水中222Rn活度的空間分布特征

天鵝洲濕地水體中Rn活度的空間分布，見圖4。

圖4 天鵝洲濕地水體中222Rn活度的空間分布Fig.4 Spatial distribution of 222Rn activity in lake water

由圖4可見，在空間分布上，天鵝洲濕地湖岸地帶及湖中心湖水中Rn活度范圍為147.00～1 245.80 Bq/m，平均值為427.37 Bq/m。其中，具有較大Rn活度的湖水樣點(diǎn)均分布在湖岸地帶，湖岸地帶湖水中Rn活度范圍為178.36～1 245.80 Bq/m，平均值為542.80 Bq/m；而湖中心湖水中Rn活度范圍為147.00～364.10 Bq/m，平均值為242.69 Bq/m。由此可見，天鵝洲濕地湖岸地帶湖水中Rn活度顯著高于湖中心湖水中Rn活度，這說明地下水向天鵝洲的排泄主要集中在湖岸地帶，攜帶大量Rn的地下水更多地排泄到湖岸地帶使湖水中具有更高的Rn活度。

湖區(qū)周邊井水中Rn活度范圍為1 814.61～5 030.35 Bq/m，平均值為3 742.00 Bq/m；而湖岸孔隙水中Rn活度范圍為1 988.48～11 583.54 Bq/m，平均值為5 207.30 Bq/m，其具有比井水更大的Rn活度變化區(qū)間和總體上更高的Rn活度。整體上看，天鵝洲濕地井水或湖岸孔隙水中Rn活度僅比湖水高大約一個(gè)數(shù)量級(jí)，與國內(nèi)外已有的大多數(shù)研究結(jié)果相比，其地下水與湖水中Rn活度的差異更小，暗示了研究區(qū)較為強(qiáng)烈的地下水向湖泊排泄的過程。

3.2 222Rn質(zhì)量平衡模型的源匯項(xiàng)

在湖泊Rn通量的匯項(xiàng)中，通過野外測(cè)量得到的湖泊上方大氣中Rn活度為16.2 Bq/m。湖水與大氣之間的Rn活度梯度會(huì)導(dǎo)致Rn向大氣中逸出。一般來說，湖泊表層風(fēng)速和溫度對(duì)湖水中Rn的損失通量影響很大。在野外工作期間，湖面上方風(fēng)速范圍為0～7.60 m/s，平均值為2.23 m/s；水溫范圍為8.40～14.60℃，平均值為10.46℃。根據(jù)公式(3)～(6)，可計(jì)算得到大氣擴(kuò)散的Rn通量為334.10 Bq/(m·d)。湖水中Rn的自身衰變通量為357.71 Bq/(m·d)。

在湖泊Rn通量的源項(xiàng)中，通過沉積物培養(yǎng)試驗(yàn)得到的沉積物孔隙水中Rn活度范圍為2 140.0～ 21 280.0 Bq/m，平均值為7 772.0 Bq/m。5個(gè)沉積物采樣點(diǎn)均勻分布于湖泊的不同區(qū)域，但它們之間Rn活度的差異較大，為了使得到的沉積物擴(kuò)散通量更加精確，分別用5個(gè)沉積物培養(yǎng)試驗(yàn)所得沉積物孔隙水中的Rn活度求得沉積物擴(kuò)散通量并取平均值。根據(jù)公式(7)～(10)，可計(jì)算得到Rn從沉積物表層擴(kuò)散到湖泊水體的通量為1.90～18.90 Bq/(m·d)，平均值為6.88 Bq/(m·d)(見表1)?；诠?1)，可計(jì)算得到地下水排泄的Rn通量為684.92 Bq/(m·d)。

表1 天鵝洲濕地水體222Rn質(zhì)量平衡模型參數(shù)Table 1 Used parameters in 222Rn mass balance model ofTian-E-Zhou wetland water body

天鵝洲濕地水體中Rn通量的源匯項(xiàng)，見圖5。

圖5 天鵝洲濕地水體中222Rn通量的源匯項(xiàng)Fig.5 Percentage of sources and sinks of 222Rn fluxes of Tian-E-Zhou wetland water body

由圖5可見，由于沒有外來地表水的輸入或輸出，在湖泊Rn通量的源項(xiàng)中只有地下水排泄和湖底沉積物擴(kuò)散，而地下水排泄的Rn通量是Rn在湖水中的絕對(duì)主導(dǎo)來源，其通量占源項(xiàng)總通量的比例為99.00%；在湖泊Rn通量的匯項(xiàng)中有大氣擴(kuò)散和自身衰變，自身衰變的Rn通量略大于大氣擴(kuò)散的Rn通量，其通量占匯項(xiàng)總通量的比例分別為51.71%、48.29%。

3.3 地下水端元選取對(duì)湖底地下水排泄通量估算的影響

已有研究發(fā)現(xiàn)，地下水端元的Rn活度值在很大程度上決定著地下水向湖泊排泄通量的大小，有時(shí)會(huì)導(dǎo)致結(jié)果相差2～3倍，因此在估算湖底地下水排泄通量時(shí)選取合適的地下水端元是非常重要的。

野外工作期間，天鵝洲濕地水面的平均面積為10.92 km，平均湖深為4.5 m。以Rn質(zhì)量平衡模型為基礎(chǔ)，計(jì)算得到地下水排泄輸入湖水中的Rn通量為684.92 Bq/(m·d)，其除以地下水端元的Rn活度，即可計(jì)算得到地下水向湖泊的排泄速率[見公式(2)]?，F(xiàn)分別以湖區(qū)周邊井水中Rn活度、湖岸孔隙水中Rn活度和沉積物培養(yǎng)試驗(yàn)所得沉積物孔隙水中Rn活度作為地下水端元，計(jì)算地下水向天鵝洲濕地的排泄通量，其計(jì)算結(jié)果見表2。

表2 選取不同地下水端元值所得的地下水排泄速率和排泄通量結(jié)果Table 2 Results of the rates and fluxes of lacustrinegroundwater discharge based on different endmember values of 222Rn activity in groundwater

由表2可知，以井水中Rn活度作為地下水端元值得到的地下水排泄速率最大，為183.04 mm/d，地下水排泄通量為2.00×10m/d，大約是以沉積物培養(yǎng)試驗(yàn)所得沉積物孔隙水中Rn活度作為地下水端元值得到的地下水排泄通量的2倍左右；其次是以井水和野外采集的湖岸孔隙水中Rn活度共同作為地下水端元值，得到的地下水排泄速率為158.85 mm/d，地下水排泄通量為1.73×10m/d，稍高于僅以野外采集的湖岸孔隙水中Rn活度作為地下水端元值得到的地下水排泄速率和地下水排泄通量。

已有研究對(duì)于Rn質(zhì)量平衡模型中地下水端元的選取一直沒有統(tǒng)一的定論，這很可能是由于湖泊類型、湖泊周邊含水系統(tǒng)結(jié)構(gòu)、地下水與湖泊的相互作用模式等的差異所引起。因此，在利用Rn質(zhì)量平衡模型量化湖底地下水排泄時(shí)，應(yīng)提前對(duì)湖泊類型、湖泊周邊含水系統(tǒng)結(jié)構(gòu)、地下水與湖泊相互作用模式等進(jìn)行詳細(xì)研究，從而得到最為可靠的地下水端元值。在本研究所在的天鵝洲濕地區(qū)，湖泊底部被區(qū)內(nèi)主要的含水層(孔隙承壓含水層)直接切割，孔隙承壓含水層中的地下水與湖水具有良好的水力聯(lián)系，而湖區(qū)周邊井水即主要取自孔隙承壓含水層中的地下水，因此在地下水端元選取時(shí)，從大尺度上看，將井水中Rn活度囊括進(jìn)Rn質(zhì)量平衡模型是十分必要的；同時(shí)，在天鵝洲濕地內(nèi)部小尺度上，通過野外測(cè)量可觀察到明顯的淺層地下徑流排泄進(jìn)入湖泊，淺層地下徑流主要發(fā)生于上覆弱透水層中的局部優(yōu)先通道(以粉砂質(zhì)為主)，因此從小尺度上看，也需要將所采集的湖岸孔隙水中Rn活度考慮進(jìn)Rn質(zhì)量平衡模型，因?yàn)槠浯砹撕^(qū)內(nèi)部淺層地下徑流。此外，由Rn的沉積物擴(kuò)散通量可知，盡管沉積物培養(yǎng)試驗(yàn)所得沉積物孔隙水中Rn活度高于井水和湖岸孔隙水中Rn活度，但其擴(kuò)散通量僅占Rn通量源項(xiàng)的1%，說明擴(kuò)散滲透的強(qiáng)度及其貢獻(xiàn)十分有限，因此地下水端元選取時(shí)可不考慮沉積物孔隙水中Rn活度。

綜上分析可知，對(duì)于天鵝洲濕地區(qū)，使用湖區(qū)周邊井水和湖岸孔隙水中Rn平均活度作為地下水端元值來進(jìn)行湖底地下水排泄通量的估算最為合適。

4 結(jié) 論

本文針對(duì)利用Rn質(zhì)量平衡模型估算湖底地下水排泄通量時(shí)地下水端元選取的不確定性問題，以長江中游典型牛軛湖——天鵝洲濕地為研究區(qū)，評(píng)估了不同地下水端元選取對(duì)Rn質(zhì)量平衡模型量化湖底地下水排泄的影響。研究發(fā)現(xiàn)，湖水中Rn活度僅比不同類型地下水低1個(gè)數(shù)量級(jí)，指示了地下水強(qiáng)烈地向湖泊排泄的過程；在湖泊Rn通量的源項(xiàng)中，地下水排泄的Rn通量占99%，而沉積物擴(kuò)散的Rn通量僅占1%；選取不同地下水端元所得的湖底地下水排泄通量從大到小依次為湖區(qū)周邊井水、湖岸孔隙水、沉積物孔隙水。根據(jù)天鵝洲濕地區(qū)孔隙承壓含水層與湖泊在大尺度上較為強(qiáng)烈的水力聯(lián)系和湖泊內(nèi)部小尺度上淺層地下徑流向湖泊的排泄，以及十分有限的沉積物擴(kuò)散滲透，確定了采用湖區(qū)周邊井水和湖岸孔隙水中Rn平均活度作為地下水端元值，所得地下水排泄速率為158.85 mm/d，地下水排泄通量為1.73×10m/d。本研究結(jié)果表明：在Rn質(zhì)量平衡模型中，地下水端元應(yīng)基于湖泊類型、湖泊周邊含水系統(tǒng)結(jié)構(gòu)、地下水與湖泊的相互作用模式等來進(jìn)行合理選定。