孫 倩

(塔里木河流域工程建設(shè)處，新疆 庫(kù)爾勒 841000)

了解不同水情條件下河流水量的去向?qū)φ麄€(gè)流域水資源的配置管理具有重大意義。河流水量的消耗包括自然耗散和人為利用，其中自然耗散指蒸發(fā)、滲漏和漫溢河道損失。塔里木河是典型的干旱內(nèi)陸河[1- 3]，塔里木河流域的地下水資源是由地表水轉(zhuǎn)化而成的，河道沉積物的滲透性將影響著地下水與地表水的補(bǔ)排關(guān)系以及補(bǔ)給量的大小。只有清楚掌握地表水與地下水的補(bǔ)排關(guān)系和補(bǔ)給量的大小，才能更合理地利用水資源，為此，需對(duì)各河段滲漏量進(jìn)行分析。

1 滲透系數(shù)計(jì)算方法

河水的下滲對(duì)于維持地下水與地表水之間的平衡起著重要作用，其最基本的影響參數(shù)即為河床沉積物的滲透系數(shù)。河道滲透系數(shù)K又稱為水力傳導(dǎo)系數(shù)[4]，是一個(gè)重要的水文地質(zhì)參數(shù)，表征流體通過(guò)多孔介質(zhì)的難易程度，直接影響到河流對(duì)含水層的補(bǔ)給量以及含水層對(duì)河流的排泄量，是對(duì)水資源開展評(píng)價(jià)的重要參數(shù)，現(xiàn)采用顆粒分析法來(lái)進(jìn)行滲透系數(shù)的測(cè)定。

1.1 顆粒分析原理

顆粒分析法是通過(guò)對(duì)沉積物的顆粒級(jí)配進(jìn)行分析，得到沉積物各個(gè)特征的粒徑(如有效粒徑d5，中值粒徑d50等)，建立滲透系數(shù)與特征粒徑之間的經(jīng)驗(yàn)關(guān)系，從而有效地得出測(cè)試點(diǎn)的K值。塔里木河干流長(zhǎng)1321km[5]，干流洪水由阿克蘇河、和田河和葉爾羌河山區(qū)暴雨及冰雪融水下泄塔河干流形成[6- 7]，上游阿拉爾站的多年平均流量為46.5×108m3，水量主要集中在汛期。在塔里木河流域，應(yīng)基于河床沉積物的顆粒級(jí)配曲線來(lái)推求河床滲透參數(shù)的大小。

1.2 數(shù)據(jù)收集處理

工作人員先后兩次沿塔里木河防洪堤壩進(jìn)行土樣采集，主要集中在上游阿拉爾至英巴扎、中游英巴扎至恰拉河段，從不同河段的起點(diǎn)開始，每隔10km取1次土樣，在距地表30、60、100、150cm處分別取土樣50～100g，全程共26個(gè)采樣點(diǎn)，土樣共計(jì)104個(gè)。然后在實(shí)驗(yàn)室對(duì)所采集的土樣顆粒進(jìn)行分析，在顆粒級(jí)配曲線的基礎(chǔ)上采用經(jīng)驗(yàn)公式計(jì)算土樣的滲透系數(shù)。通過(guò)對(duì)沉積物的顆粒級(jí)配進(jìn)行分析，得到沉積物的各特征粒徑。

以粒徑為橫坐標(biāo)，小于某粒徑的沉積物質(zhì)量百分比為縱坐標(biāo)，在半自然對(duì)數(shù)坐標(biāo)中繪制顆粒級(jí)配曲線，然后在顆粒級(jí)配曲線上查讀出所需的特征粒徑值。根據(jù)實(shí)驗(yàn)及分析結(jié)果，每個(gè)取樣點(diǎn)不同土層有效粒徑有d5、d25、d75、d95的值。

1.3 滲透系數(shù)的計(jì)算方法

基于顆粒分析法獲得河道土壤顆粒粒徑級(jí)配，根據(jù)河道本身固有特點(diǎn)選取合適合理的經(jīng)驗(yàn)公式進(jìn)行滲透系數(shù)的計(jì)算。此次研究，滲透系數(shù)的計(jì)算采用的統(tǒng)一公式如下：

(1)

式中，K—顆粒分析中計(jì)算的滲透系數(shù)(無(wú)方向)，m/d(1m/d=3.271ft/d)，C—滲透系數(shù)單位的換算系數(shù)；g—重力加速度，9.81m/s2(32ft/s2)；v—水的運(yùn)動(dòng)黏性系數(shù)，m2/s，v值取決于水溫；n—孔隙度，通過(guò)均勻度系數(shù)來(lái)近似確定；φ(n)—無(wú)量綱的孔隙度函數(shù)；de—累計(jì)重量占總重e%時(shí)的顆粒直徑，mm。系數(shù)C和φ(n)因選取方法不同而取值不同。

根據(jù)粒徑分析結(jié)果，選用Hazen計(jì)算公式來(lái)計(jì)算河床的滲透系數(shù)。在利用Hazen經(jīng)驗(yàn)公式過(guò)程中，其孔隙率函數(shù)φ(n)=1+10(n-0.26)，所對(duì)應(yīng)的有效粒徑為d5，所對(duì)應(yīng)的C值為6×10-4。則不同采樣點(diǎn)不同深度的土壤滲透系數(shù)K見表1。

表1 滲透系數(shù)計(jì)算結(jié)果

(續(xù)表)

Q=(b+2a)·L·K·t

2 滲漏量的計(jì)算方法

塔里木河干流局部河勢(shì)不穩(wěn)，河道彎曲過(guò)度，抗沖能力差[8]?；谒锬竞雍哟驳奶攸c(diǎn)，將河床近似地看做倒梯形，滲漏量的計(jì)算公式為：

(2)

式中，Q—滲漏量，m3/s;b—河床寬度，m；a—坡面長(zhǎng)度即水體與坡面接觸長(zhǎng)度，m；L—河段長(zhǎng)度，m；K—平均滲漏系數(shù);t—河水滲漏時(shí)間，d。根據(jù)宋郁東、樊自立等研究成果[3]，各監(jiān)測(cè)站水位Z與水面寬度B的計(jì)算公式為：

阿拉爾：Z1=0.409lnQ1+7.03；

(3)

B1=4×10-6exp(1.954Z1)

(4)

新其滿：Z2=0.502lnQ2+4.80；

(5)

B2=6×10-8exp(2.677Z2)

(6)

英巴扎：B3=15.22lnB2+67.68；

(7)

烏斯?jié)M：B4=9.87lnB3+25.67；

(8)

恰拉：B5=3.23lnB4+8.32

(9)

式中，Q—滲漏量，m3/s;Z—監(jiān)測(cè)站水位，m;B—水面寬度,m。

3 河道滲漏量的結(jié)果分析

基于滲透系數(shù)計(jì)算結(jié)果，利用滲漏量計(jì)算公式計(jì)算出不同水情下不同河段年內(nèi)滲漏量的結(jié)果。在合理的水情預(yù)報(bào)條件下，明確滲漏量等水量自然耗散，結(jié)合生態(tài)用水及農(nóng)業(yè)用水等，可以為水資源的合理分配提供科學(xué)依據(jù)。

3.1 阿拉爾—新其滿河段滲漏量

阿拉爾到新其滿河段為塔里木河上游[9]，在來(lái)水頻率為10%的情況下，該河段的平均滲透系數(shù)為1.07，在1—6月平均河寬為45.34m，平均水深為3m。將以上數(shù)值代入公式計(jì)算得，在來(lái)水頻率10%的情況下，1—6月阿拉爾-新期滿河段的滲漏量為1.438×108m3。7—9月，水量大，發(fā)生漫溢現(xiàn)象。計(jì)算時(shí)的河寬即為河道的最大寬度，為80m左右，水深在6～7m之間。經(jīng)計(jì)算而得，每個(gè)月的河道滲漏量均在0.39×108m3左右。10—12月，水量逐漸減少，河水的寬度也在減少。利用流量計(jì)算而得，在10—12月，河道的平均寬度為68.7m，河床平均深度為2.5m左右，則10—12月河道的總滲漏量為1.194×108m3。

在來(lái)水頻率為25%、50%、75%、90%的情況下，利用河道的流量以及河道的估算深度，可以計(jì)算出1年之內(nèi)該河段的總滲漏量，見表2。根據(jù)表2.，塔里木河阿拉爾—新其滿河段在5個(gè)來(lái)水頻率下的徑流滲漏量分別為3.805×108、3.843×108、3.163×108、3.302×108、3.59×108m3。

表2 阿拉爾—新其滿河段滲漏量(108m3)

3.2 新其滿—英巴扎河段滲漏量

該河段的平均滲透系數(shù)為0.89，由水的流量可得，在1—6月，河道的平均寬度為28m左右，水深為3.1m左右，該河段的長(zhǎng)度為258km，設(shè)河床的平均坡度為30°。當(dāng)來(lái)水頻率為10%時(shí)，在1—6月，河道的總滲漏量為0.86×108m3。在7、8、9三個(gè)月，來(lái)水過(guò)大，河水漫過(guò)河堤向外漫溢，河床的總寬度約為42m，經(jīng)計(jì)算得這三個(gè)月每個(gè)月的滲透量均為0.249×108m3。而到10—12月，河道的來(lái)水量又有所減少，經(jīng)計(jì)算而得，該段時(shí)間內(nèi)河道的總滲漏量為0.67×108m3。

利用相同的方法，可以得出在來(lái)水頻率為25%、50%、75%、90%的情況下該河段的滲透量(見表3)。根據(jù)表3，塔里木河新其滿—英巴扎河段在5個(gè)來(lái)水頻率下的徑流滲漏量分別為2.277×108、2.467×108、2.247×108、2.118×108、1.909×108m3。

表3 新其滿—英巴扎河段滲漏量(108m3)

3.3 英巴扎—烏斯?jié)M河段滲漏量

該河段河長(zhǎng)179km。平均滲透系數(shù)為1.97m/d，平均水深為2.6m，設(shè)河床的平均坡度為30°。當(dāng)河道的來(lái)水頻率為10%時(shí)，在1—6月，利用流量可得，該河段河道寬度平均為22m，河道的總滲漏量為0.994×108m3。在7—9月，河道的來(lái)水量較大，河水充滿整個(gè)河道并發(fā)生漫溢現(xiàn)象，在此時(shí)期，計(jì)算應(yīng)用的水面寬度為大約為35m，所以該段時(shí)間內(nèi)，每個(gè)月的滲漏量一致，均為0.319×108m3。10—12月，水量減少，河水的深度變淺，利用來(lái)水量與水深計(jì)算而得，該時(shí)段內(nèi)河道的滲透量為0.78×108m3。

利用來(lái)水頻率為25%、50%、75%、90%時(shí)的來(lái)水量與估測(cè)水深，便可得到這些來(lái)水頻率下的河道滲透量(見表4)。根據(jù)表4，塔里木河英巴扎—烏斯?jié)M河段在5個(gè)來(lái)水頻率下的徑流滲漏量分別為2.727×108、2.52×108、2.489×108、2.282×108、1.784×108m3。

表4 英巴扎—烏斯?jié)M河段滲漏量(108m3)

3.4 烏斯?jié)M—恰拉河段滲漏量

該河段河長(zhǎng)219km。利用顆粒分析法得該河段的平均滲透系數(shù)為2.35。當(dāng)來(lái)水頻率為10%時(shí)，在1—6月，利用流量可得該河段水面的平均寬度為18.3m，平均水深為2.1m，經(jīng)計(jì)算得河道的總滲漏量為1.22×108m3，在7—9月，河水發(fā)生漫溢現(xiàn)象，河床的總寬度為29m左右，水深大約為2.5m左右，所以每個(gè)月的滲漏量大約為0.39×108m3。到10—12月，水量減少，水面寬度下降，河水深度降低，經(jīng)計(jì)算得該時(shí)間段內(nèi)河道的總滲漏量為0.955×108m3。運(yùn)用相同的方法，計(jì)算來(lái)水頻率為25%、50%、75%、90%時(shí)的滲漏量(見表5)。根據(jù)表5，塔里木河英巴扎—烏斯?jié)M河段在5個(gè)來(lái)水頻率下的徑流滲漏量分別為3.345×108m3、3.082×108m3、3.043×108m3、2.788×108m3和2.181×108m3。

表5 烏斯?jié)M—恰拉河段滲漏量(108m3)

由以上結(jié)果可知，在以后的研究中，應(yīng)注意使用現(xiàn)場(chǎng)測(cè)定的方法，再結(jié)合其他方法進(jìn)行綜合確定，以便更為準(zhǔn)確地估算河床的滲漏量。同時(shí)，塔里木河中游還是治理的關(guān)鍵[10]。

4 結(jié)語(yǔ)

由于多種原因，塔里木河干流河道水資源滲漏量較大，直接影響沿河水資源的使用效益和成本。通過(guò)分析計(jì)算和探討，得出阿拉爾-新其滿、新其滿-英巴扎、英巴扎-烏斯?jié)M、烏斯?jié)M-恰拉各河段的平均滲漏系數(shù)分別為1.07、0.89、1.97、2.35m/d，以及各河道不同頻率下的水資源滲漏量。經(jīng)過(guò)數(shù)據(jù)對(duì)比可知，在不同來(lái)水頻率下，水資源滲漏量呈增大趨勢(shì)，即來(lái)水越多滲漏越多，越往下游，滲漏量越大，急需采取措施，加強(qiáng)河道滲漏量處置。本文的分析為工程建設(shè)提供了參考。