木兰溪下游裁弯工程数学模型的研究.pdfVIP

木兰溪下游裁弯工程数学模型研究 赵穗生 沈福新 蒙j蛊 (水利部南京水文泰资源研究所 ．o享2』00l钔 蔡一美 罗恭和 (福建省莆田市水利局莆田jjJ』口D)(福建省不利规划院福州350DOJ) 摘要本文采用非恒定流泥沙数学模型模拟术兰溪：游Ⅱ方案裁弯工程引起的 河床冲淤变化，为木兰溪下游防洪工程可行性研究和初步设计提供科学依据． 关键词木兰溪裁弯非恒定流推移质 一，木兰溪下游基本情况 本兰溪位于阌中沿海，干流由西向东注入兴化湾．其下游自木兰破至三江 口河段长25．8公里，河床平均比降O．2％，属感潮罚段． 该河段大致可分为四段： ④驶头至郑坂段，河段弯曲不大，一般河宽100多米；②郑坂至港利之 间，河段长12公里左右，平面形态为单股蜿蜒型：有十处大弯，大多为反向 弯曲．河段行洪断面小，冲淤变化大，河道演变尉烈．f可床处于淤泥层中，河 岸土质松软，河底表面有中粒径黄沙覆盖．弯段盟孝冲刷，凸岸则相应前淤； ③港利至桥兜河段，河面逐渐宽到250m左右，、可深逐渐加深到8m以上；④ 桥兜至三江口，河面更加宽敞，约j00m左右，河霖达10m以上． 木兰溪年最大洪水多发生在5～9月，洪水过基一般为三天，踱头至三江 口河段多年平均最大洪峰流量1 兴化湾潮汐为半日潮，河口多蔗平均潮差5．32m，暴大潮差7．54m，年内大潮 汛通常在农历八月或九月，月大潮汛出现在农历韧三和十八日，通常潮流界在 港利至跛头之间变动。泥沙多分布在汛期洪水季节．拈水季节舍沙量很少。 二、数学模型建立 姬前所述：木兰溪下游Ⅱ方案裁湾河段处于_{f可道演变剧烈的感潮河段。因 一143— 此，本模型采用非恒定流泥沙数学模型来模拟裁湾工程引起舟勺珂床冲淤变亿 考虑到裁弯工崔实施后河道顺直，故拟采用一维水流泥沙模型进仃模致． (一)一维非恒定流水流模型 1、基本方程组 一维渐变不恒定流方程组以流量Qfx、t)和水位Z(x、t)为未知变量，并考 虑了漫滩和旁侧人流，其表达式为： 器+日，鲁～ (1) 署+嘉f譬h筹俐警+导删 (2) 式中Z：水位： Q：流量； x：沿水流方向距离 t：时间： K：流量模数；置：且c矿i；C：土R”6 玎 q：单位河长旁铡人流； A：主槽过水断面面积； g：重力加速度； 砩：调蓄宽度，指包括滩地在内的全部河宽。 2、建立线愁叱代数方程组 (】)基本公式 本模型应用jmssmann格式对圣维南方程组进行离散。首乇设定 厂。；(△，J+．一土，J)+i(厂+，+jj 等z掣 善ze掣+簪 式中：△t——埘问步长； △x——距离步长； 日——双重系数，一般目在O．5～I．O之间，本模型取口=0．6： f—一n三一变量，如Q，Z、A及BT等。 将上述设定代人圣维南方程组，则得如下线性代数方程组： (3j 一f，△Q．。口l，△z，+Cl，△Q，+2+D“△Z，+l=Gl， (4) 2『△ZJ+l=G斗 』玎△Q，}82，△z，+C野△QJ+1+D (，=1，2，……．Ⅳ一1)