基于等效造槽流量的三维水沙数值模拟.docxVIP

基于等效造槽流量的三维水沙数值模拟 0 模拟方法及模拟条件选择 由于能够精细模拟河流的复杂演变过程，三维水沙模型得到了迅速发展和应用，但与此同时，计算量也在增加。因此探讨采用有效模拟方法,减少计算用时,提高其计算效率,对于三维水沙数值模拟具有重要实践意义。 解决上述困难的途径主要有两种:其一、开展基于并行环境的模拟研究;其二、采用概化恒定水沙过程来近似天然来水来沙过程,模拟过程中可加大计算时间步长,提高计算效率。对于前者,属于另一问题,本文不作讨论。对于后者,现有的长历时三维数值模拟多采用分级概化来代替天然水沙过程。对于河岸崩塌剧烈的冲积河流,如何采用一恒定水沙过程(等效造槽流量)作为计算条件来近似天然水沙过程,以模拟其显著平面变形特征则还有待进一步研究。 为此,本文针对河岸崩塌剧烈、河道平面变形特征显著的冲积河流,采用流量的大小来直接反映水流的强弱,并从河岸变形角度出发,针对大小流量下冲积河道平面变形存在的显著差异,提出了基于河道平面变形的等效造槽流量的计算方法,并以等效造槽流量作为计算条件,成功应用于下荆江石首河段平面变形的三维水沙数值模拟中。 1 等式制造流量 1.1 造床流量和侵蚀速率 造床流量是其造床强度和实际流量过程的综合造床作用相等的某—流量,目前造床流量的计算方法有多种,代表性的计算方法有以下三种:(1)利用实测水位流量关系,根据平滩水位确定平滩流量,代表造床流量;(2)由Wolman等提出绘制各级流量的频率和输沙率乘积的地貌功曲线,取与曲线峰值相应的流量作为造床流量;(3)韩其为提出的计算第一造床流量和第二造床流量计算方法,可根据实测水沙资料计算。 对于河岸崩塌剧烈、河道平面变形特征显著的冲积河流,本文从河岸变形角度出发,提出等效造槽流量的计算方法,其依据为在某时段内该流量作用下的河道平面变形与天然水沙过程条件下的造槽作用相当。其中河岸侵蚀速度按下式计算: 式中ζ为河岸侵蚀速度;Ee为岸坡侵蚀速率系数;Qcb为临界造槽流量;m为造槽指数,与河岸组成相关,当河段大水造槽效果越显著时,该值越大。 由时段内河岸变形幅度相当可知,等效造槽流量Qe与流量过程有如下关系: 则ΔT时段内河段的等效造槽流量Qe表示为: 1.2 荆江河段崩岸与流量的关系 由式(3)可知,确定ΔT时段内河段等效造槽流量Qe的两关键参数为:临界造槽流量Qcb及造槽指数m。现选取具有显著平面变形的下荆江石首河段作为计算实例加以说明。 临界造槽流量Qcb可根据河段发生崩岸的临界时段来反推。该河段坍岸的年内变化,一般为5～7月坍塌强度较小,7～10月最大,10～12月又逐渐减弱。因此可以坍塌强度较小月份的平均流量来估算其值,根据1996-1998年河段实测流量资料,临界流量Qcb取值8500m3/s。 在大江大河中,观测河道崩岸及其与流量的关系是很困难的工作,20世纪60年代初,长江科学院与荆江河床试验站合作,对荆江河段崩岸进行了专门研究,唐日长等人通过分析,发现河岸崩塌总量与时段内平均流量的平方(Q2)及持续时间(Δt)存在较好关系,河岸崩塌总量随Q2Δt的增大而增大,如图1所示。为此,本实例取式(3)中的造槽指数m=2.0。从等效造槽流量的计算式可知,该值越大,表示高水期的造槽效果越明显,可解释为高水期崩岸强度较大的同时,崩塌下来的泥沙被带走的比例亦较大,而对近岸的掩护作用则较弱;与此相反,在流量较小时,崩岸强度较小,崩塌下来的泥沙被带走的比例亦较小,而坍塌物对近岸的掩护作用则较强。采用等效造槽流量可较好地反映大小流量下河道崩岸过程的差异,充分体现高水期显著的造槽效果。对于长江中下游来说,河岸坍塌主要发生在汛期及汛后退水期,其中7～8月发生崩岸的次数几乎占全年的50%。因此对于石首河段来说,其造槽系数m较大,为2.0左右。 根据式(3)可计算出该河段1996-1998年期间的等效造槽流量Qe约为19100m3/s,与实际工程应用较多的平滩流量法所确定的数值(19000m3/s左右)较为接近,大于该时段平均流量13000m3/s。计算表明,临界流量Qcb的取值对等效造槽流量Qe值的影响并不十分显著,Qcb取值7500 m3/s～9500 m3/s时,该河段1996-1998年期间的等效造槽流量Qe的取值范围约为18900 m3/s～19700m3/s。 2 下荆江河岸崩岸建模 水流模型求解及边界条件与文献[8、9]一致,泥沙模型及崩岸模拟简介如下。 1.悬沙输移方程 三维悬移质输沙方程在笛卡尔坐标系下可表示为: 式中u、v、w分别为沿x、y、z方向流速;sk为第k组含沙量;ωsk为相应泥沙沉速;εs为泥沙扩散系数。 2. 推移质不平衡输沙方程 推移质不平衡输沙方程: 式中qb及qb*分别为实际推移质输沙率和平衡推移质输沙率;Db、Eb分别表示悬移质