水利工程論文-沖積河流懸移質泥沙與床沙交換機理及計算方法研究.doc

水利工程論文-沖積河流懸移質泥沙與床沙交換機理及計算方法研究摘要：本文在分析前人研究成果的基礎上，從非恒定挾沙水流沙量平衡入手，建立了沖積河流懸移質泥沙與床沙交換關系表達式，并引入平衡沖淤物粒徑的概念，對沖淤物粒徑的計算方法進行了系統(tǒng)分析，提出了一整套計算非恒定挾沙水流懸移質泥沙與床沙交換的模式與方法，同時將該方法引用到數(shù)學模型中，采用黃河下游典型洪水進行了驗證。關鍵詞：懸移質泥沙床沙交換機理平衡沖淤物粒徑1引言懸移質泥沙與床沙交換機理及計算方法，作為河流泥沙動力學的重要問題之一，長期以來，被許多學者所關注，并取得了卓有成效的進展。該問題的研究，對河道演變分析、河工模型設計、河流泥沙數(shù)學模型計算以及河流泥沙動力學學科的發(fā)展都頗有意義。現(xiàn)有方法一般為將懸移質泥沙和床沙分成若干組，通過計算分組挾沙力實現(xiàn)懸沙與床沙的交換，從而求出懸沙與床沙級配的變化規(guī)律。其中韓其為是由懸移質級配求分組挾沙力，該方法假定含沙量級配與挾沙力級配相同，且各粒徑組泥沙獨自保持平衡而不相互影響，在此基礎上建立了淤積和沖刷兩種情況的級配計算公式1，但在計算過程中出現(xiàn)與實際河流沖淤不相符合的情況，如河床發(fā)生沖刷時，各粒徑組的泥沙均發(fā)生沖刷，而淤積時各粒徑組泥沙均發(fā)生淤積。HEC-6模型2是以床沙級配求分組挾沙力，其基本思路為：先求每一粒徑組泥沙的可能挾沙力，即這一粒徑組含沙量占100%的水流挾沙力，而實際挾沙力則是此粒徑組泥沙在床沙中所占百分比與上述可能水流挾沙力的乘積，這種方法的前提假定與實際挾沙圖形有較大出入。李義天3通過理論分析建立了輸沙平衡狀態(tài)下床沙質級配與床沙級配之間的關系為(1)(2)式中Pi,Pbi分別為床沙質與床沙級配；i為第i粒徑組泥沙沉速；u*為摩阻流速；為卡門常數(shù)；m為粒徑分組數(shù)；v為垂向紊動強度，原作者取v=u*。在已知床沙級配的條件下，根據(jù)式(1)計算出懸沙中床沙質級配后，再由全沙挾沙力乘以挾沙力級配Pi即可得出分組挾沙力，進一步采用非飽和輸沙模型計算河段內的沖淤，爾后根據(jù)沖淤結果調整床沙級配作為下一時段的計算依據(jù)。錢意穎、曲少軍、張啟衛(wèi)、韋直林4提出的分組挾沙力計算方法在數(shù)學模型計算中也頗有應用價值。但應該承認，水沙運行過程中，各粒徑組泥沙是相互影響的，人為地分組后分別計算顯然與實際出入較大。特別象黃河這樣的多沙細沙河流，細顆粒泥沙對水流流變特性及挾沙特性影響尤為顯著。此外，現(xiàn)有通過計算分組挾沙力來確定泥沙級配的方法多應用于恒定流模型，而天然河流一般為非恒定流，泥沙級配不僅沿程變化，而且還隨時間變化，故此進一步研究懸沙與床沙交換機理有重要意義。2沖積河流床沙與懸沙交換模式及計算方法初步研究2.1懸沙與床沙交換過程中懸沙粒徑變化規(guī)律研究如圖1所示，取1-2之間的x微段，Q為流量；S為含沙量；A為過水面積；(d)為懸移質泥沙粒徑分布函數(shù)；t為時間微段；為沖淤物干容重；A0為橫斷面淤積面積；c(d)為沖淤物粒徑分布函數(shù)；qL、SL、L(d)分別為側向入流單寬流量、含沙量、入流懸沙粒徑分布函數(shù)。則t時段內，x河段上di的凈流入量P1為(3)微段x的水體沙量增加中粒徑為di的泥沙所占的量P2為(4)側向入流與沖淤使di的增加量P3為(5)由質量守恒原理，則P1+P3=P2，將式(3)(5)代入，兩邊除以xt，并忽略微小量，整理得(6)式(6)兩邊乘以di，并對d積分后化簡，同時將懸沙平均粒徑、沖淤物平均粒徑、側向入流泥沙平均粒徑代入可得(7)根據(jù)多元函數(shù)求導法則將上式展開，并代入泥沙連續(xù)方程，進一步化簡得(8)采用式(7)或式(8)即可計算任一時刻任一斷面的懸沙平均粒徑大小。對式(6)兩邊乘以$d2i$,積分后并取懸沙粒徑二階圓心距；沖淤物粒徑二階圓心距；側向入流泥沙粒徑二階圓心距，則化簡可得(9)或(10)求得進口d過程后,即可依據(jù)式(9)或(10)計算出任一時刻任一段面的懸沙粒徑二階圓心距。2.2床沙與懸沙交換過程中床沙粒徑變化規(guī)律研究取單位寬度的斷面，設t時段床沙交換層厚度為Hc；河床高程為Z；床沙粒徑分布函數(shù)為(Di)；床沙平均粒徑為Dcp。則對于Di的泥沙含量,瞬初x1=Hc(Di)；瞬末；沖淤部分。由質量守恒知x2-x1=x3,同時等號兩邊同除以t，可化簡為(11)式兩邊乘以Di并對D積分得上(12)根據(jù)多元函數(shù)求導法則把式(12)展開，并將代入，可得(13)采用式(13)即可計算出床沙粒徑變化過程。同時，對式(11)兩邊乘以，并對D積分可導得(14)式中D為床沙粒徑二階圓心距。采用式(14)即可計算出床沙粒徑二階圓心距的變化過程。2.3沖淤物粒徑dc的確定天然沖積性河流的懸移質泥沙與床沙交換的過程中，當河床沖刷時，懸沙中較粗顆粒仍會落淤下來，而河道淤積時，床面上也有部分細顆粒由床面被水流沖起，當沖淤平衡時，也存在部分懸沙與床沙不停地進行著交換，只是落淤的懸沙與沖起的床沙的質量與粒徑大小相等。由式(8)、(13)可知，推求懸沙與床沙的變化規(guī)律集中反映在沖淤物平均粒徑dc的確定上。天然河流的懸沙一般較床沙為細，dc應介于Dcp與dcp之間。當沖淤平衡時，沖淤物平均粒徑可定義為平衡沖淤物粒徑，采用dph表示。2.3.1淤積時dc的計算方法進一步分析可知，河床發(fā)生淤積時，dc應介于dcp與dph之間。當淤積強度較小時，dc接近于dph，當淤積強度較大時，dc接近于dcp，于是可初步確定dc的如下關系式dc=mdcp+(1-m)dph(15)式中m為介于01之間的權重系數(shù)，其取值大小應與相對淤積強度有關。若令初始單位面積柱體中水體的含沙量為S，水深為h，則水體中含有的泥沙總量為hS；河床淤積后，落淤在河床上泥沙的厚度為z，則淤積泥沙質量為z。于是可將m表示為如下關系式m=(z/hS)k(16)式中k為指數(shù)。進一步可求出淤積后水體中剩余部分泥沙的平均粒徑為d=(hSdcp-zdc)/(hS-z)(17)在含沙量不太高的情況下，懸沙淤積時分選作用較明顯，當zhS時，應有d0，于是將式(15)、(16)代入式(17)后，并求極限，即可導出k與dcp及dph的關系式k=dcp/(dph-dcp)(18)將式(16)、(18)代入式(15)，即可得出淤積時dc的計算公式(19)2.3.2沖刷時dc的計算方法對于河床沖刷情況而言,dc應介于dph與床沙平均粒徑Dcp之間。沖刷強度較小時,dc接近于dph,當河床發(fā)生強烈沖刷時,dc接近于Dcp。采用上述類似的方法，可確定dc的計算式為dc=mDcp+(1-m)dph(20)式中m為介于01之間的系數(shù)，其取值大小與相對沖刷強度有關。若將床沙混合層分為兩層，上層為床沙與水流直接接觸分選層，定義為直接交換層，下層為分選后的床沙通過沙波運動進行床沙調整層。則對于單位面積的河床而言，令沖刷時床面的直接交換層厚度為Hc，沖刷厚度為z，可把m表示為m=(z/Hc)k(21)式中k為指數(shù)。于是可計算出沖刷后直接交換層內剩余部分床沙的粒徑D=(HcDcp-zdc)/(Hc-z)(22)當河床發(fā)生強烈沖刷時，zHc，則D應接近于床沙組成中的極粗部分床沙粒徑D。于是將式(20)、(21)代入式(22)并對其求極限可求得k的表達式為k=(D-Dcp)/(Dcp-dph)(23)將式(21)、(23)代入式(20)得(24)若采用此時水流條件求出的泥沙揚動粒徑為df，則強烈沖刷剩余部分床沙組成中泥沙粒徑應均大于df，這樣D可由初始床沙級配曲線反求得出。2.3.3dph的確定目前尚無dph的實測資料，且其影響因素極為復雜，故dph的確定十分困難。由上述分析可知，dph應介于dcp與Dcp之間，為簡單起見，不仿采用類似dc的方法來提出dph的如下表達式dph=(1-k1)dcp+k1Dcp(25)上式中的k1為介于01之間的權重系數(shù)。經分析發(fā)現(xiàn)，k1的取值大小和dcp與Dcp的比值有關。當Dcpdcp時,dph的大蝎只取決于dcp,而與Dcp無關。為此,可初步假定k1與dcp/Dcp成正比，即k1=k2(dcp/Dcp)(26)系數(shù)k2的取值區(qū)域為0,1，經數(shù)學模型初步