刘川顺;陈曦濛;鲁晓义
【摘 要】River bank erosion induced by ship waves can be a serious threat to the land resources and the dykes. Dangers given by the yacht-generated waves usually are neglected because of yacht's small sizes, however, the high speeds of yachts can result in significant damages to the river banks. As there is generally little quantitative information available to assess the magnitude of erosion,the objective of this paper is to establish an available model for estimating the recession of the river bank caused by ship-generated waves and provide an application example of the model. The study includes three main components:field measurement and investigation, laboratory experiment, and numerical modeling. The field measurement consists of soil sampling from banks of the Changgang River in Ezhou city and measurement of the river characteristics; site investigation includes learning about the types of recreational yachts sailing on the Changgang River and their parameters, as well as daily navigation frequencies. The laboratory experiments of the soil
samples were carried out to determinate the soil particle size distribution curves and the critical shear stresses,so as to estimate the riverbank erodibility. Finally,the yacht-generated wave factors were calculated, and the numerical modeling of the bank erosion rates of the Changgang River was conducted by using a lateral erosion model, from which, an annual bank recession of the Changgang River was obtained as 54. 4 cm, approximately consistent with the field survey result.%船行波导致的河岸侵蚀对土地资源和堤防构成威胁.游艇体量小,其危害往往被忽视,然而快速航行的游艇会对河岸产生严重损坏.目前量化评估船行波侵蚀河岸的文献很少.建立船行波引起河岸侵蚀后退的估算模型,并结合鄂州市长港河侵蚀进行实例计算.研究过程由三部分组成:一是现场勘测和调研,包括采集河岸土样,测量河道特征参数,调查长港河船舶类型、船舶技术参数及各类船舶每天航行频次;二是开展室内土工试验,测定土样颗粒级配曲线和临界剪应力,以评估河岸的抗侵蚀性;三是开展河岸侵蚀的数值计算,包括计算游艇船行波要素,采用侧向侵蚀模型计算河岸侵蚀率.计算结果表明,在游艇船行波侵蚀作用下,长港河河岸年平均侵蚀后退54.4 cm,与现场调研结果基本相符. 【期刊名称】《水利水运工程学报》
【年(卷),期】2017(000)004
【总页数】6页(P8-13)
三体船【关键词】游艇;船行波;河岸侵蚀;波浪要素;侧向侵蚀模型
【作 者】刘川顺;陈曦濛;鲁晓义
【作者单位】武汉大学 水资源与水电工程科学国家重点实验室, 湖北 武汉 430072;武汉大学 水资源与水电工程科学国家重点实验室, 湖北 武汉 430072;武汉大学 水资源与水电工程科学国家重点实验室, 湖北 武汉 430072
【正文语种】中 文
【中图分类】TV861
波浪侵蚀率指单位时间内波浪造成的侵蚀发展速率，主要与波浪要素、土质、水域特征有关。在某一特定地区其土质与水域特征一定，则波浪侵蚀率主要由波浪要素决定。船行波的波浪要素主要包括波速、波长、最大波高以及波周期。
船行波的研究始于19世纪60年代，Rankine于1868年最早发表了关于船行波问题的论著，英国学者Froude则在1877年对船行波的形态做了一系列研究，并绘制出了船行波形态图[1]。船行波理论研究的主要奠基人造船学家Kelvin在1887年首次从理论上对船行波进行了分析研究[2]，他以一个在水面上作等速直线运动的压力扰动点代替船舶研究其对水体的作用，对于无限水深的情况，得出压力扰动点在水面移动将激起一系列的散波和横波；激起的散波和横波波峰线可以用以波浪传播角θ为参数的等相位线方程式表示；横波将在扰动点处以扰动速度的50%向后扩散，形成新波；散波系相对航线向外凸出，任意点波速等于压力扰动点在垂直于单个波峰线上的速度分量vcosφ(φ为航线与波传播方向之间的夹角)，散波波长等于航线上的波长乘以cos2φ。1949年，前苏联包瑞奇在莫斯科运河实船试验和模型试验的基础上提出了计算波高的包瑞奇公式，把船速作为主要影响因素。1958年，前苏联向金基对于包瑞奇的研究，在理论分析和数学推导的基础上，引入干扰系数K，进一步得到了船舶沿运河轴线行驶时，靠近岸坡的船行波波高[3]。1977年Dawson提出了Rankine源方法即直接对各片体物面压力积分求得三体船的兴波阻力[4]，被视为船行波波浪载荷问题计算的突破性进展。1990年荷兰Delft水工试验所提出计算波高的经验公式，第1次使用了船型修正系数[5]。1997年Sorensen发表了不同航速区除波高外的波浪要素计算公式[6]。
我国针对船行波的研究起步较晚，且理论方面的研究不多，系统的船模试验开始于20世纪60年代。1963年，余广明等对大运河船行波问题进行了室内模型研究，提出了紧靠船边和临近岸坡的最大波高影响因素(相对船速、断面系数、相对水深)及各因素分别与紧靠船边的相对最大波高的相关曲线[7]。1992年，南京水利科学研究院、江苏省交通规划设计院合作进行现场测试和室内船模试验，研究成果验证了Delft波高公式的合理性[8]。
在船行波侵蚀河岸研究方面，荷兰Delft水工研究所1981—1983年在Hartel航道系统开展了护岸工程结构型式和滤层材料的模型试验研究[8]。Osman等[9]提出了水流对黏性岸坡侧向侵蚀模型，同时给出了相关参数侵蚀速率经验常数Se和土壤临界切应力τc的相关关系。Trimbak等利用侧向侵蚀模型，在结合实测数据系列的前提下，对船行波作用下密西西比河流域的岸坡冲蚀进行了数值模拟[10]。
目前我国量化评估船行波侵蚀河岸的理论研究和工程应用比较少见，本文在总结国内外相关研究的基础上，试图系统地建立船行波河岸侵蚀的计算模型，并结合鄂州市长港河侵蚀进行实例计算。
河岸在船行波作用下侵蚀后退距离取决于单个船舶船行波对河岸的侵蚀速率和各船舶船行
波的作用时间。船舶类型及其数量、船行波作用时间可由现场观测和调查获得，而单个船舶船行波的侵蚀速率则需要根据船舶参数(吨位、尺寸、航速)通过以下步骤计算获得。
1.1 船行波波浪要素计算
船行波由初生波和次生波组成。Hemphill提出，当船行波相对断面系数k<0.05时，船舶航行时激起的波浪主要以次生波为主，初生波可忽略[11]。
式中：As为该工况下船舶吃水线以下的船体横截面积；Ac为静水位时的河道横断面面积。
船行波理论研究、原体观测及室内模型试验结果都表明，次生波的波态与船速及水深有密切关系。1997年Sorensen提出用水深弗汝德数Fd 来判断快艇航行时所处航速区[6]：
式中: Fd为水深弗汝德数，无量纲；V为船速；g为重力加速度；d为水深。
当Fd<0.75时，航速位于亚临界流速区，对应深水次生波波态。该波态由横波和散波两种类型的波组成。横波传播方向与船行驶方向相同，散波是接近平行的一组短波，沿航向对称分布，传播方向与航向的夹角约为55°。
当Fd>0.75时，对应浅水次生波波态。其中当0.75≤Fd<1.00时，航速位于跨临界流速区，次生波散波波峰线与航线夹角随着Fd的增大逐渐增大接近90°，形成首横波和尾横波，散波与横波合一，波高达到最大值，波峰线长度随速度的增大而增大。当Fd≥1.00时，航速位于超临界流速区，横波消失，只剩下散波。散波的波峰线较长，波峰线与航线夹角随航速增加而逐渐减小。
根据1997年Sorensen发表的不同航速区波浪要素(不含波高)计算公式[6]可知，亚临界流速区波浪要素计算式如下：
式中：θ为波浪传播角(取θ=35.3°)；c为波速；L为波长；T为周期。
跨临界流速区波速c的计算同式(3)，周期T同式(5)，其余波浪要素计算式如下：
式中：波长L通过试算确定；tanhx为双曲余切函数。
超临界流速区波速c的计算同式(3)，波长L同式(7)，周期T同式(5)，波浪传播角计算式如下：
南京水利科学研究院的研究成果说明我国的内河航道最大波高在进行计算时基本符合Delft公式[8]，因此本次计算采用该公式。
式中：Hm为最大波高；A为船型修正系数，无量纲，根据参考文献[8]，巡逻艇、满载内河马达船取1.0，空载货船取0.5，空载马达船取0.35，拖牵船队取0.42； S为船舶与岸坡的距离。
1.2 船行波侵蚀河岸的计算
Osman和Thorne研究了船行波对黏土岸坡的侧向侵蚀，建立了计算平均侵蚀速率的侧向侵蚀计算模型 [9]：
式中：b为岸坡侵蚀后退距离; db/dt为岸坡侵蚀速率；τ为土壤有效切应力；τc为土壤临界切应力；se为侵蚀速率经验常数，Osman和Thorne通过开展黏土河岸侵蚀试验，给出了侵蚀速率常数与土壤有效切应力之间的关系式[9]：
式中：ρ为水的密度。
Temple等[12]提出了确定各种土壤临界切应力的公式，其中黏性土临界切应力按下式计算：
式中：τcb为土壤基本允许切应力；Ce为孔隙比修正系数。根据参考文献[12]，黏土孔隙比修正系数和基本允许切应力由下式确定：
式中：e为孔隙比；Ip为塑性指数。
波浪对河岸土体产生的有效切应力τ用Parchure计算式进行计算[10]：
式中：fw为波浪摩阻系数；ρ为水的密度，取1 000 kg/m3；ub为近底波浪质点水平速度。
根据微幅波理论，fw及ub可采用以下计算式：
式中：ɑ为波浪振幅，即有效波高的一半，有效波高取最大波高的一半；k为波数，即在波浪传播方向上2π长度内出现的全波数目；h为波浪序列在岸坡影响区底部的水深，即波高的一半；Ks为相对粗糙度，无量纲，根据Soulsby[13]推荐的Ks与粒径间的关系为Ks=2.5d50。
综合以上数据可求得各船型船行波对土质岸坡的侵蚀速率，然后乘以各船型1年内通过次数及每次船行波的作用时间，累计便可得到该河岸在1年内的侵蚀后退距离。
2.1 长港河河岸侵蚀崩塌
长港河位于湖北鄂州市，河道全长46.5 km，流域总面积3 265 km2，是梁子湖湖水向长江排泄的唯一通道。长港河旅游通航河段为梁子湖湖口至磨刀矶节制闸间1.5 km长河段，主要通航小型快艇和小型游艇。该河段为冲积平原湖区地形地貌，河道纵坡1/20 000，常水位18.00 m，常水位时平均河面宽度约55 m，河底深泓高程约12.00 m，两侧地面高程一般在17～20 m(见图1)。根据现场调查，长港河通航河段一般年份河岸侵蚀后退0.5 m以上，局部甚至每年崩塌后退1.0～1.5 m，威胁着两岸房屋建筑及农田安全。
2.2 河岸侵蚀崩塌原因分析
(1)河岸抗侵蚀能力较弱。河岸土层分布从上至下依次为：第四系人工堆积物(Q4ml)，主要为黏土或粉质黏土夹少量碎石、砂砾，厚3～8 m；第四系全新统湖积物(Q4l)，属湖积成因的黏土，夹极薄粉细砂层，厚3～6 m；第四系上更新统残积物(Q3el)，黄，以黏土为主，
夹少量砂粒，厚0.3～2.0 m；侏罗系下统武昌(Jlwc)。岸坡表层主要是黏土，从试验获得：土壤粒度累积分布曲线见图2，河道边坡粉质黏土层重度18.7 kN/m3，孔隙比0.904，液限36.7%，塑限18.8%，塑性指数17.8，液性指数0.82，内摩擦角10°，凝聚力24 kPa，抗冲性能较弱。长港河常水位18 m，汛期平均水位20 m，水位变幅区除了零散生长的柳树、杨树，基本无其他植被，河岸抗波浪侵蚀的能力较差，在船行波的频繁淘刷作用下，树根间的土壤被淘空，河岸边的树最终倾倒，缺乏植被保护的土质河岸被船行波侵蚀成陡坎，逐渐崩塌后退。
(2)旅游快艇的频繁波浪冲刷。水流、风浪、船行波浪都可能侵蚀河岸，但是对长港河而言，前两者可以忽略不计：其一，崩塌河段的河道坡降仅1/20 000，河道内水流流速很小，而且一年之中绝大部分时间磨刀矶节制闸是关闭的，即河水不流动；其二，长港河水面狭窄，风浪很小。此外，长港河侵蚀崩塌只发生在通航河段的事实更是佐证了该河段的河岸侵蚀崩塌是船行波导致的。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议