第五章 海岸工程建后冲淤预报模式教案.docVIP

第五章 海岸工程建后冲淤预报模式 §5.1 淤积率计算公式的推导 窦国仁的潮汐水流悬沙运动微分方程为 对式(1)在一个潮周期内积分,可近似得到时均方程(略去时均符号)为: 采用有限差分法,式(2)可写作: 则工程后单位面积上的淤积体积,即一个潮周期的淤积强度为: 假定建造工程后不改变边界的来沙量,即: 则一年中淤积强度为： (8) 式中潮周期T单位为秒，工程前后的挟沙能力S*1，S*2的单位为kg/m3,n是一年中潮数。 §5.2 计算参数的确定 式(8)中一年中潮数n=706，S*的计算，前人已作过很多研究，较常用的有武汉水利学院公式： （9） 此外，李昌华曾总结得出： （10） 如将上式化成武汉水利学院公式相同的形式，则得： （11） 对于潮汐河口，其挟沙能力可以写成： （12） 式中k和m由当地实测资料确定，文献[33]曾分析认为k和m是水流及泥沙条件的函数。不同河口水流泥沙条件均不一样，因此，各地的k和m也不一样。 采用何种形式的挟沙能力公式，最重要的是要符合当地的实际情况，通常采用当地的实测流场和含沙量资料，通过相关分析的方法来建立适合当地情况的挟沙能力公式。 在通常情况下,不论淤泥质泥沙在分散状态下的粒径多么小,其絮凝沉降的当量粒径均在0.015～0.03mm范围内,在淤泥质海岸的回淤问题计算中,絮凝沉速取常值,即ω=0.0004～0.0005m/s。 根据公式(8)以及选定的参数，结合潮流数学模型计算工程前后流场变化的数值，可以计算回淤强度。 §5.5 淤积预报模式的应用 考虑到南汇附近海域的实际情况，通过相关分析，得到该水域的水流挟沙能力公式： 式中u为全潮平均流速，g是重力加速度，H为全潮平均水深，A=290，m=1.28。将S*代入式（8），得到： （13） 泥沙的干容重和沉降机率的选取与前一章相同。 在本海域，由于一个潮周期的泥沙沉降距离大于海域的水深，即 ω*TH，所以实际计算中以各区域的全潮平均水深H来代替ω*T。H=h+ζ，其中h为理论基准面以下水深，ζ为平均潮位。这样式(13)则为： （14） 回淤强度的计算区域见图（27）。计算结果见表5-1。 图（27） 表5-1 工程实施后年回淤强度（94年～95年） 单位：米/年 区域编号 计算值 实测值 区域编号 计算值 实测值 1 0.20 0.14 11 0.35 0.26 2 0.30 0.20 12 0.36 0.30 3 0.29 0.15 13 0.37 0.31 4 0.34 0.22 14 0.25 0.20 5 0.30 0.22 15 0.27 0.22 6 0.38 0.28 16 0.26 0.21 7 0.40 0.34 17 0.22 0.15 8 0.31 0.25 18 0.21 0.18 9 0.35 0.30 19 0.20 0.14 10 0.39 0.32 从表中可以看出，用所建立的淤积强度预报模式计算出的数值与实测值是比较接近的，因此，可以用该淤积预报模式来预估后几年的淤积强度。 我们用得到的回淤强度，对计算区域的地形资料进行修正，再利用数模在新的地形条件下算出的流场，用该淤积预报模式来预估后一年的回淤强度。以此类推，可得到今后几年各年的淤积强度，计算结果见表5-2。 表5-2 工程实施后年回淤强度计算值 单位：米/年 区域编号 1996年 1997年 1998年 1999年 2000年 1 0.19 0.16 0.12 0.09 0.07 2 0.28 0.24 0.21 0.15 0.10 3 0.30 0.23 0.18 0.14 0.08 4 0.32 0.25 0.22 0.16 0.09 5 0.29 0.24 0.20 0.15 0.08 6 0.35 0.30 0.23 0.16 0.10 7 0.38 0.34 0.24 0.18 0.10 8 0.29 0.25 0.20 0.15 0.07 9 0.30 0.26 0.22 0.16 0.07 10 0.35 0.28 0.23 0.17 0.09 11 0.33 0.27 0.21 0.14 0.08 12 0.34 0.29 0.22 0.15 0.09 13 0.35 0.30 0.24 0.16 0.09 14 0.23 0.20 0.16 0.13 0.08 15 0.25 0.21 0.17 0.13 0.09 16 0.25 0.2