河流模拟课程设计—水库一维泥沙淤积计算(编辑修改稿)

河流模拟课程设计—水库一维泥沙淤积计算(编辑修改稿)内容摘要：

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 一维泥沙淤积计算 ~ 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 坝前断面变化曲线 一维泥沙淤积计算 ~ 12 ~ 六、 结果分析 剖面形态分析 A、 纵断面淤积情况分析： 由深泓变化曲线图可以看出，随着时间的增加，河道深泓点高程成 增高趋势。 该水库运行 10年后，在 距坝约 10000米处，深泓点高程增加 40多米， 计算结果 表明 ， 该 水库泥沙 以淤积为主，致使水库库容逐年损失，且情况较为严重。 B、坝前断面变化曲线分析： 由坝前断面变化曲线图可以看出，随着时间的增加，坝前断面各点的总体高程也成增高趋势。 该水库运行 10后，该断面 河床抬高 20多米 ， 表明在该处河道逐年以淤积为 主 ，且淤积情况较为严重。 库容损失合理性分析 由水位库容曲线可以看出，水库在运行 10年后，经计算可知，在 273m水位条件下，库容的损失率为 %，在 241m的水位条件下，库容的损失率高达 %，计算结果表明，经过十年的运行，水库库容基本成全部淤满状态，所以在 10 年后将不能再继续使用。 结果表明，该水库的调节能力较差。 上述计算结果 与现实不太 相符 ，其原因是在计算该水库库容的泥沙损失时，计算冲往下游的泥沙 量不合理。 在给定条件下，该水库并没有考虑水位的变化及坝前冲沙情况。 该水库的计算条件为正常蓄水位保持不变，在高水位运行条件下，必然会导致水库的严重淤积，而实际情况是，在洪水来临时，坝前水位要降至防洪限制水位，以增大防洪库容，该条件下，水流的速度及水流挟沙力较大，会造成库区泥沙处于冲刷状态，以增加水库库容。 一维泥沙淤积计算 ~ 13 ~ 七、 计算程序 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。 逐年累计淤积量 .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。 逐年水位库容关系变化 .txt39。 ,status=39。 unknown39。 ) !**************读入地形数据 **************** call Qbed(mm,npxt,x,dxa,alow0,npoint) !**************计算沉降速度 **************** W=FW(T,r,rs,d,0) write(*,*)W !**************计算断面间距 **************** do i=1,npxt1 dx(i)=dxa(i)dxa(i+1) enddo !****计算初始水面线、深泓、库容水位关系 **** dWQST= dWQS= dWGb= 一维泥沙淤积计算 ~ 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) !*计算各个时段水库的淤积情况 * DO 40,I=1,ND WRITE(*,*)39。 读入天数 :39。 ,I !***读入流量、沙量数据 *** 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) !**计算各断面的水流挟沙力和推移质输沙率 ** CALL SGB(nn,a,b,Q,Sx,Gb,rough,w,ndisp) !******计算各断面的含沙量和冲淤厚度 ****** CALL DEPO(nn,mm,w,x,S,Sx,Gb,dt,B,Q,dx,dy,alow,alow0,npoint,faidy,ndisp) 一维泥沙淤积计算 ~ 15 ~ !**************检验质量守恒 ************** CALL CHECK(dy,B,S,Gb,dt,npxt,Zlevel,x,npoint,Q,dx,ndisp) 41 CONTINUE !******计算悬移质、推移质年均流量 ******** WQS1=Q(1)*S(1) Gb(1)=WQS1* WQSN=Q(nn)*S(nn) dWQST=dWQST+WQS1*dt/10000000 !悬移质累计输沙量 dWQS=dWQS+(WQS1WQSN)*dt/!悬移质累计冲淤量 dWGb=dWGb+(Gb(1)Gb(nn))*dt/!推移质累计淤积量 AVQ=AVQ+Q(1)*dt/10000 !累计流量 !****输出相关文件 *** 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 一维泥沙淤积计算 ~ 16 ~ !**************读入地形数据 ****************** 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 !*******计算深泓 ******* 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 !*泥沙沉降速度的计算 ,采用张瑞瑾公式 * FUNCTION FW(。