【文章內(nèi)容簡介】 6 8 10 6 267 8 267 10 267 6 8 267 10 267 6 8 267 10 267 6 8 267 10 267 6 8 267 10 267 6 8 10 267 6 8 10 一維泥沙淤積計(jì)算 ~ 11 ~ 6 8 10 6 8 10 6 8 10 6 8 10 6 8 10 6 8 10 6 8 10 6 8 10 6 8 10 6 8 10 6 8 10 6 8 267 10 267 6 8 10 6 8 10 6 8 10 6 8 10 6 8 10 圖 5 壩前斷面變化曲線 一維泥沙淤積計(jì)算 ~ 12 ~ 六、 結(jié)果分析 剖面形態(tài)分析 A、 縱斷面淤積情況分析： 由深泓變化曲線圖可以看出，隨著時間的增加，河道深泓點(diǎn)高程成 增高趨勢。該水庫運(yùn)行 10年后，在 距壩約 10000米處，深泓點(diǎn)高程增加 40多米， 計(jì)算結(jié)果 表明 ， 該 水庫泥沙 以淤積為主，致使水庫庫容逐年損失，且情況較為嚴(yán)重。 B、壩前斷面變化曲線分析： 由壩前斷面變化曲線圖可以看出，隨著時間的增加，壩前斷面各點(diǎn)的總體高程也成增高趨勢。該水庫運(yùn)行 10后，該斷面 河床抬高 20多米 ， 表明在該處河道逐年以淤積為 主 ，且淤積情況較為嚴(yán)重。 庫容損失合理性分析 由水位庫容曲線可以看出，水庫在運(yùn)行 10年后，經(jīng)計(jì)算可知，在 273m水位條件下，庫容的損失率為 %，在 241m的水位條件下，庫容的損失率高達(dá) %，計(jì)算結(jié)果表明，經(jīng)過十年的運(yùn)行，水庫庫容基本成全部淤滿狀態(tài)，所以在 10 年后將不能再繼續(xù)使用。 結(jié)果表明，該水庫的調(diào)節(jié)能力較差。 上述計(jì)算結(jié)果 與現(xiàn)實(shí)不太 相符 ，其原因是在計(jì)算該水庫庫容的泥沙損失時，計(jì)算沖往下游的泥沙 量不合理。 在給定條件下，該水庫并沒有考慮水位的變化及壩前沖沙情況。該水庫的計(jì)算條件為正常蓄水位保持不變，在高水位運(yùn)行條件下，必然會導(dǎo)致水庫的嚴(yán)重淤積，而實(shí)際情況是，在洪水來臨時，壩前水位要降至防洪限制水位，以增大防洪庫容，該條件下，水流的速度及水流挾沙力較大，會造成庫區(qū)泥沙處于沖刷狀態(tài)，以增加水庫庫容。 一維泥沙淤積計(jì)算 ~ 13 ~ 七、 計(jì)算程序 program main parameter(nn=16,mm=60,ndisp=0,npxt=16,ND=3653,NY=10,T=15) paramete (rs=,r=,d=,failev=,faidy=,dt=86400) dimension x(mm,2,nn),rough(nn),npoint(nn),B(nn),A(nn),Gb(nn),dy(nn),Nday(NY) dimension dx(nn),Xw(nn),dxa(nn),alow(nn),alow0(nn),Q(nn),Sx(nn),S(nn),Zlevel(nn) data(Nday(m),m=1,10)/366,731,1096,1461,1827,2192,2557,2922,3288,3653/ !**************打開輸出文件 **************** open(10,file=39。河床地形 .txt39。,status=39。old39。) open(11,file=39。流量沙量 .txt39。,status=39。old39。) open(12,file=39。逐年累計(jì)淤積量 .txt39。,status=39。unknown39。) open(13,file=39。逐年深泓變化 .txt39。,status=39。unknown39。) open(14,file=39。逐年水面線變化 .txt39。,status=39。unknown39。) open(15,file=39。逐年壩前斷面淤積變化 .txt39。,status=39。unknown39。) open(16,file=39。逐年水位庫容關(guān)系變化 .txt39。,status=39。unknown39。) !**************讀入地形數(shù)據(jù) **************** call Qbed(mm,npxt,x,dxa,alow0,npoint) !**************計(jì)算沉降速度 **************** W=FW(T,r,rs,d,0) write(*,*)W !**************計(jì)算斷面間距 **************** do i=1,npxt1 dx(i)=dxa(i)dxa(i+1) enddo !****計(jì)算初始水面線、深泓、庫容水位關(guān)系 **** dWQST= dWQS= dWGb= 一維泥沙淤積計(jì)算 ~ 14 ~ AVQ= do I=1,NPXT rough(I)= enddo Q(1)=295. do I=2,NPXT Q(I)=Q(I1) enddo Zlevel(npxt)= CALL level(x,rough,nn,zlevel,dx,q,npoint,b,a,xw,nn,mm,failev,NDISP) !初始水面線 !********輸出到文件 ********** CALL FILE(12,13,14,15,16,0,x,npoint,Zlevel,alow,dxa,dx,mm,nn,dWQST,AVQ,dWQS,dWGb,ndisp) !*計(jì)算各個時段水庫的淤積情況 * DO 40,I=1,ND WRITE(*,*)39。讀入天數(shù) :39。,I !***讀入流量、沙量數(shù)據(jù) *** READ(11,*)Q(1),S(1) IF(Q(1)==0)THEN GOTO 41 ENDIF DO J=2,NPXT Q(J)=Q(1) ENDDO !***************************************** CALL level(x,rough,nn,zlevel,dx,Q,npoint,b,a,xw,nn,mm,failev,0) !**計(jì)算各斷面的水流挾沙力和推移質(zhì)輸沙率 ** CALL SGB(nn,a,b,Q,Sx,Gb,rough,w,ndisp) !******計(jì)算各斷面的含沙量和沖淤厚度 ****** CALL DEPO(nn,mm,w,x,S,Sx,Gb,dt,B,Q,dx,dy,alow,alow0,npoint,faidy,ndisp) 一維泥沙淤積計(jì)算 ~ 15 ~ !**************檢驗(yàn)質(zhì)量守恒 ************** CALL CHECK(dy,B,S,Gb,dt,npxt,Zlevel,x,npoint,Q,dx,ndisp) 41 CONTINUE !******計(jì)算懸移質(zhì)、推移質(zhì)年均流量 ******** WQS1=Q(1)*S(1) Gb(1)=WQS1* WQSN=Q(nn)*S(nn) dWQST=dWQST+WQS1*dt/10000000 !懸移質(zhì)累計(jì)輸沙量 dWQS=dWQS+(WQS1WQSN)*dt/!懸移質(zhì)累計(jì)沖淤量 dWGb=dWGb+(Gb(1)Gb(nn))*dt/!推移質(zhì)累計(jì)淤積量 AVQ=AVQ+Q(1)*dt/10000 !累計(jì)流量 !****輸出相關(guān)文件 *** DO m=1,NY IF(I==Nday(m))then CALL FILE(12,13,14,15,16,m,x,npoint,Zlevel,alow,dxa,dx,mm,nn,dWQST,AVQ,dWQS,dWGb,ndisp) ENDIF ENDDO 40 CONTINUE close(10) close(11) close(12) close(13) close(14) close(15) close(16) end 一維泥沙淤積計(jì)算 ~ 16 ~ !**************讀入地形數(shù)據(jù) ****************** subroutine Qbed(mm,nn,x,dxa,alow0,npoint) dimension x(mm,2,nn),dxa(nn),alow0(nn),npoint(nn) do i=1,nn read(10,*) read(10,*) dxa(i),npoint(i) read(10,*) do j=1,npoint(i) read(10,*) x(j,1,i),x(j,2,i) enddo !*******計(jì)算深泓 ******* alow0(i)=x(1,2,i) do m=2,npoint(i) if(alow0(i)x(m,2,i))then alow0(i)=x(m,2,i) endif enddo enddo end !*泥沙沉降速度的計(jì)算 ,采用張瑞瑾公式 * FUNCTION FW(