吉倩倩

【摘 要】針對黃土高原小流域降雨徑流模擬試驗中，污染物運移規(guī)律不明確，污染物運移比尺難以確定的問題，本文利用環(huán)境水力學(xué)中污染物擴散定理及相似論的基本原理，進行了一系列相似轉(zhuǎn)換和推導(dǎo)，最終得出了污染物運移濃度的比尺為。該結(jié)論從理論上給出了原型與模型在降雨徑流過程中，污染物運移濃度之間的關(guān)系，為進一步進行模型試驗進行污染物濃度率定和驗證提供參考。

【關(guān)鍵詞】降雨徑流；污染物擴散；相似比尺；相似率

中圖分類號： X53 文獻標(biāo)識碼： A 文章編號： 2095-2457（2018）12-0123-002

DOI：10.19694/j.cnki.issn2095-2457.2018.12.053

0 引言

降雨沖刷地表沉積物而產(chǎn)生含有大量污染物的徑流流入水體，對水環(huán)境造成了嚴重的危害[1]。20世紀60年代以來，國外首先注意到了降雨徑流污染導(dǎo)致的水質(zhì)惡化的問題，并且在降雨徑流水質(zhì)特性研究、水環(huán)境影響評價、模型模擬方面做出大量的工作[2-3]。近些年，我國也開展了大量的降雨徑流污染物運移規(guī)律的研究，在黃土坡面土壤氮磷等化學(xué)物質(zhì)隨徑流遷移的研究中，張亞麗等[4]運用了降雨-徑流-土壤化學(xué)物質(zhì)混合層深度（簡稱混合層深度）的方法，模擬坡面土壤化學(xué)物質(zhì)隨徑流遷移的過程，王敬美等[5-6]運用模型試驗結(jié)合相似理論的方法研究了水中污染物的擴散規(guī)律，高建恩等[7-8]運用相似性原理對黃土高原小流域水力侵蝕模擬試驗的比尺進行了設(shè)計和驗證，可作為模擬小流域水力侵蝕調(diào)控，優(yōu)化治理措施研究的有效措施。這些研究通過運用模擬實驗、相似理論推導(dǎo)等方法，從不同的側(cè)重點探究降雨徑流中水流運動和污染物運移的規(guī)律，但沒有從定量的角度推求出污染物在模型與原型中運移濃度的關(guān)系。本研究是建立在相似理論的基礎(chǔ)上，通過將污染物三元對流擴散方程與相似理論建立聯(lián)系，從而推求出原型與模型的污染物運移濃度比尺，為模型試驗的驗證提供了理論依據(jù)。

1 研究背景

1.1 費克定律

污染物擴散是環(huán)境水力學(xué)的核心部分，水流紊動與擴散有密切的聯(lián)系。研究污染物在水體中的擴散規(guī)律是推導(dǎo)出污染物運移比尺的基礎(chǔ)工作，因此具有重要意義。

1885年費克提出擴散定律：在各向同性的介質(zhì)中，在一定方向上通過單位面積擴散輸送的物質(zhì)，與該斷面的濃度梯度成比例[9]。公式如下：

F=-D （1）

式中：F——單位時間通過單位面積的擴散物質(zhì)，又稱為通量；C——擴散物質(zhì)的濃度； ——x方向的濃度梯度；D——擴散系數(shù)，在各向同性的液體中，在x、y、z三個方向上，擴散系數(shù)為常數(shù)。

上述公式只表示出污染物在一個方向上的擴散，實際當(dāng)中，污染物質(zhì)不是在單一方向而是在x、y、z三個方向上都有擴散，因此，將其寫成矢量形式應(yīng)為：

F=-D（ + + ）=-D？犖C（2）

1.2 單純擴散方程

在菲克定律的基礎(chǔ)上，通過質(zhì)量守恒原理建立單位時間擴散物質(zhì)的量的變化等于x，y，z三個方向的通量之和的關(guān)系，推出靜止流體中污染物單純擴散方程。其基本形式為：

=D（ + ）（3）

如考慮兩個方向上的擴散，則為： =D（ + ）（4）

如只考慮一個方向的擴散，則： =D （5）

1.3 對流擴散方程

以上所述的單純擴散，流體是靜止的，但是如果水體在流動，在x、y、z三個方向上的流速為u、v、w。由于流體的平均運動而引起的輸送稱為對流，對流輸送與擴散輸送是可以分開又可以迭加的過程。

將對流通量和擴散通量相加即為總輸送率（通量），代入單純擴散方程（3）中，并首先考慮x一個方向，通過化簡可得：

+ =D 或 +？犖（uC）=D？犖2C（6）

因為，流體的容積守恒，？犖·u=0，則上式可變?yōu)椋?+u·？犖C=D？犖2C

將其寫成三元形式，則有：

+u +v +ω =D（ + + ）（7）

這就是三元對流擴散方程，簡稱三元擴散方程。

通常只考慮x方向的縱向流速，y和z方向忽略，則可忽略v 和ω ，x方向的擴散項與對流輸送項相比很小，因此可以忽略 ，簡化后可得：

+u =D（ + ）（8）

如果在y方向（垂向）的濃度梯度較小，可忽略y方向的擴散項 ，則為：

+u =D （9）

以上公式中，（6）為一元擴散方程、（7）（8）為三元擴散方程、（9）為二元擴散方程。

2 分析推導(dǎo)

由于實際水流多數(shù)屬于三度紊流，所以采用三元擴散方程（7）進行比尺推導(dǎo)。公式形式為：

+u +v +ω =D（ + + ）

假定采用正態(tài)模型，設(shè)長度比尺為λl，則λx=λy=λz=λl

首先，水流運動滿足連續(xù)性方程 + + =0，很容易由比尺關(guān)系式推出λu=λv=λω

其次，水流運動滿足比尺關(guān)系式 =1以及 =1。容易得出λu=λ 。

現(xiàn)將公式（7）中原型的有關(guān)物理量用比尺轉(zhuǎn)化為相應(yīng)的模型物理量，可得：

用 去除每一項，可得：

描述統(tǒng)一物理現(xiàn)象的的原型和模型的方程應(yīng)該相等[10]，所以：

擴散系數(shù)D表示單位濃度梯度下的擴散通量，它表達某個組分在介質(zhì)中的擴散快慢，是物質(zhì)的一種傳遞性質(zhì)。相同的溶質(zhì)在相同溶劑中的擴散系數(shù)是相同的，若原型與模型采用相同的污染物，則擴散系數(shù)D相同，λD等于1。

則上述比尺關(guān)系式可變?yōu)椋?=1，所以：

可得污染物的濃度比尺為λC=λ 。

3 驗證方法

3.1 確定擴散方式及其解的形式

上述多種形式的基本方程都是二階線性偏微分方程，從物理學(xué)的角度來看，必須給定初始條件及邊界條件，才可求解，從而得出在空間上或時間上的擴散物質(zhì)的濃度分布情況。

總結(jié)起來污染物在水中的擴散方式有以下五種情況：（1）一元瞬時平面源；（2）二元及三元的瞬時點源擴散；（3）起始空間分布源；（4）時間連續(xù)源；（5）對流擴散。對應(yīng)于每一種擴散方式，其都有不同的解。

以一元瞬時平面污染源為例，其污染物添加方式為在水管的中部瞬時投放污染源，投放點坐標(biāo)設(shè)為0，采用的方程為一元擴散方程： =D

其解為：C= exp- ，在t=0時，在x=0處投放的質(zhì)量為M；

3.2 模型試驗的設(shè)計與布置

按照一定的幾何比尺，水流運動比尺和坡面侵蝕產(chǎn)沙及泥沙運動比尺布置原型和模型，采用相同的污染物添加方式，在對應(yīng)的投加點施加污染物，如氮肥或磷肥，原型模型施加量可按推導(dǎo)比尺λC=λ 施加。

采用放水沖刷和降雨徑流兩種外部自然環(huán)境條件的模擬，在不同斷面位置、不同時間，測量徑流中污染物的濃度，將所測數(shù)值代入擴散方程所確定的解中，計算出原型與模型的污染物濃度，驗證污染物運移濃度比尺的準(zhǔn)確性。

4 結(jié)論

（1）本研究將環(huán)境水力學(xué)中污染物擴散方程與相似理論中比尺模擬的方法相結(jié)合，依據(jù)原型與模型滿足同一物理方程的思想，運用水力學(xué)的基本定理，推導(dǎo)出污染物運移相似比尺，是一種大膽的創(chuàng)新和嘗試。

（2）為進一步進行坡面降雨徑流污染物運移的研究提供了新的思路和方法。

（3）在污染物施加量的確定和驗證效果的對比上，還需要通過一系列模型試驗驗證和修改，有待于進一步提升。

【參考文獻】

[1]張德剛，陳永川，湯利.滇池流域典型城郊降雨徑流污染特征與排放過程[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2010，29（4）：733-739.

[2]Zhiguo He，Haoxuan Weng，Hao-Che Ho.Soil Erosion and Pollutant Transport during Rainfall-Runoff Processes.Water Resources，2014，vol41，No5，604-611.

[3]Francirose Shigaki，Andrew Sharpley，Luis Ignacio Prochnow. Rainfall intensity and phosphorus source effects on phosphorus transport in surface runoff from soil trays. Science of the Total Environment， 2007，373：334-343.

[4]張亞麗.黃土坡面土壤氮磷等化學(xué)物質(zhì)徑流遷移與混合層深度研究[D].西安：西安理工大學(xué)，2007.

[5]王美敬，羅麟，盧紅偉，等.水中污染物擴散模型實驗中的相似理論[J].四川大學(xué)學(xué)報（工程科學(xué)版），2004，36（2）：25-28.

[6]王美敬，羅 麟，黃小雪. 污染物橫向擴散在變態(tài)模型與河流中的相似性[J].四川大學(xué)學(xué)報（工程科學(xué)版），2006，38（4）：25-28.

[7]高建恩，楊世偉，吳普特，等.水力侵蝕調(diào)控物理模擬試驗相似律的初步確定[J].農(nóng)業(yè)工程學(xué)報，2006，22（1）：27-31.

[8]高建恩，吳普特，牛文全，等.黃土高原小流域水力侵蝕模擬試驗設(shè)計與驗證[J].農(nóng)業(yè)工程學(xué)報，2005，21（10）：41-45.

[9]張書農(nóng).環(huán)境水力學(xué)[M].南京：河海大學(xué)出版社， 1988：52-87.

[10]河流泥沙工程學(xué)（下冊）[M].北京：水利出版社，1982：167-183.