應(yīng)用邊界層方法確定溶質(zhì)遷移參數(shù)的實驗研究

1、應(yīng)用邊界層方法確定溶質(zhì)遷移參數(shù)的實驗研究摘要：本文以時域反射儀(TieDainRefletetry,TDR)作為溶質(zhì)鋒的探測手段，對利用邊界層方法確定的參數(shù)和利用穿透曲線擬合法確定的參數(shù)進展了比擬，結(jié)果說明邊界層方法確定的擴散彌散系數(shù)(D)與穿透曲線擬合法確定的擴散彌散系數(shù)根本相近。但，因為TDR靈敏度的限制，邊界層方法確定的延遲因子(R)大于穿透曲線擬合法確定的延遲因子。應(yīng)用邊界層方法確定的參數(shù)，比擬了用邊界層方法和準確方法預(yù)測的濃度剖面圖，說明邊界層方法在一定時間內(nèi)可以準確地預(yù)測污染物的動態(tài)濃度分布。實驗結(jié)果不僅說明了TDR作為溶質(zhì)鋒探測手段的可行性，同時說明了邊界層方法在確定溶質(zhì)遷移參數(shù)

2、方面有一定的可靠性。關(guān)鍵詞：多孔介質(zhì)溶質(zhì)遷移邊界層理論溶質(zhì)鋒DE模型是描繪化學(xué)物質(zhì)遷移規(guī)律的數(shù)學(xué)模型，該模型的應(yīng)用需要兩個根本的參數(shù)：水動力彌散系數(shù)D和延遲因子R，因此準確確定化學(xué)物質(zhì)在特定多孔介質(zhì)中遷移的模型參數(shù)D和R是模型得以應(yīng)用的前提。關(guān)于模型參數(shù)的估計，已提出了許多方法14，而這些方法存在著收斂性和參數(shù)唯一性的問題。邊界層方法是利用溶質(zhì)鋒遷移位置與時間的關(guān)系進展DE方程參數(shù)的估計，可以省去做穿透曲線的費事。但此方法還僅僅停留在理論層次上，尚未從理論的角度加以驗證，況且如何準確確定溶質(zhì)鋒的位置還是一個有待解決的問題。本文通過對邊界層方法和傳統(tǒng)的穿透曲線擬合方法進展比擬，初步確定利用TDR

3、探測確定溶質(zhì)鋒位置具有可行性。1理論根底邵明安等(1998)通過對DE方程的研究，在均質(zhì)穩(wěn)態(tài)條件下，結(jié)合邊界層理論，提出了確定DE方程中參數(shù)的邊界層理論模型5。該理論通過對農(nóng)業(yè)污染物在多孔介質(zhì)中遷移過程中濃度鋒面的描繪和概化，得到了下述表達式：(1)式中：d(t)表示t時間溶質(zhì)鋒的位置。D是水動力彌散系數(shù)，R是延遲因子，V是平均孔隙流速。土壤中的駐留濃度：r(x,t)=vd(t)0/vd(t)+3D(1-x/d(t)3(2)對于半無限土柱或田間土壤剖面，相應(yīng)于通量型入流邊界條件，土壤溶質(zhì)遷移過程可以利用邊界層理論進展描繪。由方程(1)可知，只要求得時間t時的溶質(zhì)鋒的深度d(t)，我們就可確定溶

4、質(zhì)遷移參數(shù)D、R。為了便于對數(shù)據(jù)t和d(t)進展處理，下面對方程(1)進展進一步簡化。兩邊同時平方，得d2(t)=4vt/Rd(t)+12Dt/R(2)方程兩邊同除以t2，得1/t=R/12Dd2(t)/t2-v/3Dd(t)/t(3)2材料與方法表1實驗設(shè)置及土柱物理參數(shù)土樣高度/s/(g/3)v/(g/3)Vp/(/h)探測點個數(shù)水頭/黃綿土462.741.302.1448.0沙壤土682.701.341.0565.5婁土482.651.271.2549.0供試土樣分別為榆林沙壤土耕層、安塞黃綿土耕層及楊凌婁土耕層，過2篩，實驗所用的有機玻璃管高為100，管內(nèi)徑為14.玻璃管壁上每隔10鉆

5、內(nèi)徑為2的小孔，以便插TDR(時域反射儀)探頭，共有6個。TDR為波蘭EASYTEST公司消費的F/ts型，可以同時測定土壤中的水分、鹽度和土壤溫度。探頭選擇10長的雙針式探頭?？刂埔欢ǖ娜葜兀捎梅謱犹钛b法填裝土柱，采用馬氏瓶控制水頭。根據(jù)邊界層方法，實驗在穩(wěn)態(tài)飽和條件下進展。由于邊界層方法是在均質(zhì)穩(wěn)態(tài)不考慮源匯項的條件下建立的，所以本實驗中選擇化學(xué)性質(zhì)不活潑的l-為遷移離子。微機自動控制TDR，連續(xù)測定探測點電解質(zhì)的濃度變化。下端用100l的容量瓶接取出流液，獲得該土柱l-的穿透曲線。AgN3滴定法測定出流液中l(wèi)-含量。轉(zhuǎn)貼于論文聯(lián)盟.ll.3結(jié)果與討論3.1溶質(zhì)鋒位置的判斷邊界層方法的應(yīng)

6、用首先需要確定溶質(zhì)鋒的位置，而溶質(zhì)鋒是一個遷移物質(zhì)在遷移過程中從無到有的界面，而所采用的設(shè)備手段總是有一定的靈敏度，所以準確探測溶質(zhì)鋒有一定的難度。常用的監(jiān)測溶質(zhì)遷移的技術(shù)設(shè)備一般有溶液采樣器和電阻感應(yīng)探頭，而TDR那么是另一種新的用于溶質(zhì)遷移研究的設(shè)備手段68。本實驗所采用的TDR的詳細性能參數(shù)見表2。實驗中采用長度為10的雙針式探頭，它的感測范圍是一個內(nèi)徑為5，長度為13圓柱體，與實驗所采用的土柱相適宜。表2F/ts型TDR性能參數(shù)測量范圍測量精度分辨率土壤含水量0%100%2%0.1%總鹽/(s/)0.000110%0.0001土壤溫度()-20+500.8%0.1表3邊界層測量數(shù)據(jù)沙壤

7、土黃綿土婁土探測點d(t)/t/hd(t)/t/hd(t)/t/h1st114.28112.63113.452nd219.62216.32218.783rd3113.48319.953113.584th4119.174113.654118.675th5126.226th6136.18由于實驗處理為飽和穩(wěn)態(tài)條件，在溶質(zhì)鋒未到達探測點前，TDR所探測到電解質(zhì)的濃度為一恒定值。當該值發(fā)生變化時，也就是電解質(zhì)濃度開場增大時，就斷定在這個時刻，溶質(zhì)遷移到了該點，這樣就得到了溶質(zhì)鋒遷移到第一個探測點時的時間。隨后將TDR移向下一個探測點。各探測點位置及溶質(zhì)鋒到達的時間見表3.圖1為沙壤土土柱中TDR在前5

8、個點監(jiān)測電解質(zhì)濃度變化的局部過程圖，圖2為TDR在第6個點的完好過程圖。從圖1可以看出，每個監(jiān)測點的電解質(zhì)濃度變化曲線從開場的一段時間內(nèi)是一條直線，濃度沒有變化，此時我們認為溶質(zhì)鋒還沒有到達該點。隨時間的推移，濃度開場逐步增加，我們認為濃度開場增加的那個點對應(yīng)的時間是溶質(zhì)鋒到達該點的時間。從圖1還可以看出，隨著探測點的下移，電解質(zhì)濃度變化的曲線逐漸變緩，也就是遷移物質(zhì)的擴散間隔 逐漸增大，濃度梯度逐漸變小，這與實際情況相符。從圖2可以看出，TDR記錄的完好的電解質(zhì)濃度變化過程圖的形狀跟穿透曲線的形狀極其相似，這也說明了用TDR探測溶質(zhì)鋒的可行性。圖1沙壤土前5個探測點濃度變化過程圖2沙壤土第6

9、探測點濃度變化過程3.2穿透曲線方法求D、R參數(shù)穿透曲線形狀不僅可以反映溶質(zhì)遷移的機理和溶質(zhì)與土壤之間的作用，而且可以把穿透曲線的數(shù)據(jù)借助XTFIT程序用最小二乘法進展擬合，得到溶質(zhì)遷移的兩個重要參數(shù)擴散彌散系數(shù)D和延遲因子R.因為l-的化學(xué)性質(zhì)不活潑且?guī)ж撾姾桑瑧?yīng)選用XTFIT程序中的第二個數(shù)學(xué)模型即線形平衡吸附模型進展擬合。其控制方程為：圖3邊界層數(shù)據(jù)擬合曲線(4)式中：f為通量型溶質(zhì)濃度。求得的3種土壤的l-遷移參數(shù)，見表4.表4兩種方法求得D、R參數(shù)對照沙壤土黃綿土婁土參數(shù)BTBLTBTBLTBTBLTD2.092.264.445.131.473.39R0.932.441.163.34

10、0.752.883.3邊界層方法求D、R參數(shù)對表3中的數(shù)據(jù)按式(3)的形式再進一步處理，進展二次多項式擬合，分別得到l-溶質(zhì)鋒在3種土壤中運動的數(shù)學(xué)表達式(5)(7)及圖3.沙壤土：y=0.0902x2-0.1548x(5)r2=0.8091黃綿土：y=0.0552x2-0.1391x(6)r2=0.9971婁土：y=0.0633x2-0.1103x(7)r2=0.9896擬合的相關(guān)系數(shù)較高，擬合的結(jié)果具有可信性。根據(jù)公式各項系數(shù)可求得溶質(zhì)l-遷移參數(shù)D、R，見表4.3.4兩種方法的比擬與分析從表4可以看出，兩種方法得到的水動力彌散系數(shù)D值除婁土的相差比擬大外，其他兩種土樣的D值比擬相近，而用

11、邊界層方法得到的3種土樣的延遲因子R值都比用穿透曲線方法獲得的R值大。其原因可能是由于TDR所斷定的溶質(zhì)鋒的位置比實際溶質(zhì)鋒的位置靠后所致。因為TDR有一個靈敏度，只有當電解質(zhì)的濃度變化量到達儀器的靈敏度時，濃度的變化才被指示出來。而理論上邊界層是遷移物質(zhì)從無到有的一個層面，層面上的濃度變化很小時，TDR探測不到，所以等TDR的探測值發(fā)生變化時，溶質(zhì)鋒已經(jīng)過了探測點的位置。溶質(zhì)鋒位置確定的延遲在表觀上就等同于遷移離子的延遲因子R增大，使得溶質(zhì)鋒向前推移的速度變慢。從式(2)中也可以看出，在其他參數(shù)不變的前提下，隨d(t)增大，R也隨之增大，這正與我們的上述分析相符。為檢驗邊界層方法獲得參數(shù)的可

12、應(yīng)用性，我們可以用邊界層方法擬合的駐留濃度剖面與用DE方程的準確解擬合的濃度剖面進展比擬。根據(jù)方程(2)r(x,t)=vd(t)0/vd(t)+3D(1-x/d(t)3(8)在d(t)、D的情況下可求得從溶質(zhì)鋒到土柱上邊界的濃度剖面圖。而DE方程的解析解的表達式為9：r(x,t)/0=1/2erfRx-vt/2(DRt)0.5+(v2t/DR)1/2exp-(Rx-vt)2/4DRt-f(x,t)(9)f(x,t)=1/2(1+vx/D+v2t/DR)exp(vx/D)erfRx+vt/2(DRt)0.5(10)圖4圖93種土樣以邊界層方法獲得D、R參數(shù)代入方程(8)和方程(9)計算的濃度剖面

13、比照圖。圖4沙壤土4.28h濃度剖面圖5沙壤土26.2h濃度剖面圖6黃綿土2.62h濃度剖面圖7黃綿土13.65h濃度剖面圖8婁土3.45h濃度剖面圖9婁土18.67h濃度剖面由圖4圖9可以看出，在一定時間內(nèi)，邊界層計算的結(jié)果同準確方法計算的結(jié)果根本一致。但隨時間的增加，邊界層方法計算的結(jié)果與準確方法計算的結(jié)果之間的誤差會逐漸增大，說明用邊界層方法在給定遷移參數(shù)情況下預(yù)測遷移物質(zhì)濃度分布時，給定的時間不能過長，即不能預(yù)測日期太久后的濃度分布。4結(jié)論研究結(jié)果說明，TDR可以用來監(jiān)測溶質(zhì)鋒的遷移，與擬合穿透曲線法相比，邊界層方法確定的參數(shù)具有一定可靠性。TDR因受其靈敏度的制約，所判斷的溶質(zhì)鋒比實

14、際的位置靠后，導(dǎo)致邊界層方法確定的延遲因子R大于擬合穿透曲線方法得到的R，但兩種方法得到的D值相差不大。從兩種方法預(yù)測的在不同時段內(nèi)的遷移物質(zhì)動態(tài)分布來看，在一定的時間內(nèi)，用邊界層方法預(yù)測遷移物質(zhì)的濃度分布剖面具有較高的精度。參考文獻：1RifaiNE,KaufanJ,TddDK.DispersinphenenainlainarflthrughprusediaJ。Inst.fEng.Res.1956,93(2).SanitaryEng.Lab.Univ.falif.,Berkeley.2ElprineA,DayPR.Fittingslutebreakthrughequatinstdatausi

15、ngtadjustableparaetersJ。SilSi.S.A.J.,1977,41:3941.3KlJB,ParkJ,VanGenuhtenTH.Paraeterestiatinfrunsaturatedflandtransprtdel-arevieJ。Hydrl.1987,91:255293.4ButhterB,Hiniz,Flury,FluhlerH.HetergeneusFlandslutetransprtinanunsaturatedstnysilnlithJ。SilSi.S.A.J.,1995,59:1421.5Shainganetal.Anapprxiateslutintthenvetin-dispersinequatinfslutetransprtinsilJ。SilSiene,1998,163:339345.6BearJ.SeexperientsinDispersinJ。GephysRes.1961,66:24552467.7allantsD,Vanlster,NTride,JVanderbrght,JFeyen.parisnfthreeethdstalibrateTDRfrnitringsluteveentinundisturbe