会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

   首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-泥石流泥沙体积浓度之研究.doc

  • 资源星级:
  • 资源大小:100.56KB   全文页数:12页
  • 资源格式: DOC        下载权限:注册会员/VIP会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-泥石流泥沙体积浓度之研究.doc

水利工程论文泥石流泥沙体积浓度之研究摘要本研究系探讨泥石流先端部和整体平均包括先端部及后续水流泥沙体积浓度随着溪床坡度之变化规律,以及其纵剖面泥沙体积浓度沿着水深的分布特性。经采用其他研究者渠槽试验资料进行检定及验证后发现,当系数ξ0.1时,本研究所提出泥石流平衡泥沙体积浓度之半经验公式,不仅可模拟泥石流先端部泥沙体积浓度随着溪床坡度之变化趋势,亦可适于推估不同溪床坡度下之泥石流整体平均泥沙体积浓度。同时,本研究也应用最大熵原理及平衡泥沙体积浓度半经验公式,导出了纵剖面泥沙体积浓度分布之理论公式,经Tsubaki实测资料验证结果显示,本文模式可以掌握纵剖面泥沙体积浓度之分布趋势而优于Tsubaki及Chen两研究者之理论公式。关键词泥石流体积浓度模式分析1前言泥石流是指发生于陡峻的溪谷或斜坡面上未固结之沉积物或经风化之砾石、岩屑等堆积土体,受豪雨形成的地表迳流或土层地下水位上升等水文过程之影响,而失去原有安定之状态,土砂砾石伴随着洪流在重力作用下沿着自然坡面形成的一种高浓度集体流动现象,对区域环境安全和生产的影响甚巨,而逐渐引起人们的重视。在过去的二十余年间,已有不少学者专家投入泥石流的形成、运动、停积及治理措施等各方面课题之研究,并获得若干成果。但受限于对泥石流形成过程及其流动特性变化规律的了解程度,致使在探讨其行为机制的各项相关因子时,仍不免存在着一些盲点有待解决,例如泥石流泥沙体积浓度或泥沙含量的变化规律,就是一个相当显著的例子。泥石流泥沙体积浓度系指单位体积中固体泥沙颗粒所占有的体积,属于泥石流基本流体特性之一,其变化规律密切影响着泥石流的行为机制,而与泥石流运动、流出土方量、泛滥范围及其成灾规模等问题关系密切,为泥石流防灾相当重要的指标参数之一,惟目前除了日本高桥Takahashi的理论模式外,相关研究成果却是寥寥可数。根据高桥建构泥石流平衡泥沙体积浓度理论公式所引用的理论基础,可以有两种不同的表达形式,一是以巴格诺Bagnold高泥沙含量固液两相流的本构关系为基础另一则是运用了水力学原理。兹分述如下巴格诺1从实验中发现,当固液两相流中泥沙含量很高时,由于固体颗粒间可以透过不断地碰撞来进行动量转换而形成两种不同的力量,其中沿着水流方向的动量交换会形成粒间离散剪应力τ,interparticledispersiveshearstress,而与水流垂直方向的动量交换则形成粒间离散应力Pd,interparticledispersivestress。同时,他也提出了粒间离散剪应力与粒间离散应力间的本构关系constitutiverelationship可表为τPdtanα1式中,tanα动摩擦系数,当水流处于完全惯性区fullyinertialrange时,tanα≈0.321。高桥2即利用此一关系式,经由泥石流固液两相的动量守恒方程导出了泥石流平衡泥沙体积浓度公式Cdρtanθ/ρsρtanαtanθ2式中,Cd平衡泥沙体积浓度θ溪床倾斜角度ρs、ρ固体泥沙颗粒及流体之密度。由上式得知,因不能使用于tanθ≥tanα≈0.32的溪床条件,加上泥石流常发生在溪床倾角θ≥15°的陡峻地形上,使得应用式2在计算泥沙体积浓度时,因合理值范围相当有限而难以符合实况。例如,取tanα0.32、tanθ0.268θ15°、ρs2.6g/cm3及ρ1.0g/cm3代入式2,则平衡泥沙体积浓度Cd3.2,该值已超出泥沙体积浓度之合理值,实际上是不存在的。因此,高桥运用水力学原理推导平衡泥沙体积浓度公式3Cdρtanθ/ρsρtanφtanθ3式中,tanθ静摩擦系数。与式2比较得知,上式分母中系采用静摩擦系数tanθ,而非动摩擦系数,是两方程式间的主要不同点。惟高桥经实验验证显示,式3的形式在模拟泥石流泥沙体积浓度变化规律上优于式2,故乃建议采用式3作为预测稳定均匀条件下的泥石流先端部泥沙体积浓度。不过,Chen4曾对此提出看法认为,泥石流既已达稳定均匀流况,其流动型态自应处于完全惯性状态,固体泥沙颗粒间的摩擦作用理当采用动摩擦系数作为参数,较为合理,而非静摩擦系数。笔者5则自非牛顿流体流变方程式出发,认为泥石流在高浓度、高剪应变率时,由于颗粒间始终保持直接接触,不同粒径的石砾相互叠置碰撞,使颗粒间除承受来自于颗粒间弹性碰撞所产生的离散剪应力作用外,还必须考虑颗粒间之摩擦阻力。这样,泥石流先端部泥沙体积浓度可表为Cd32βρtanθ/ρsρ2β3tanθ3tanα4式中β摩擦阻力系数resistancecoefficient,且0≤β≤1。根据上式,若不考虑粒间摩擦阻力系数,取β0时,可获得和式3度表达式。欧国强6于流槽通过实验方式研究泥石流先端部及其后续水流两部分之整体平均泥沙体积浓度,并建立了以溪流坡度为主的泥沙体积浓度经验公式,即Cdt4.3Cmtanθ1.5/14.3Cmtanθ1.55式中Cdt泥石流整体包括先端部及后续流平均泥沙体积浓度Cm溪床面静止泥沙之最大泥沙体积浓度thegrainconcentrationinvolumeinthestaticdebrisbed。综合以上各泥沙体积浓度公式绘出泥沙体积浓度与溪床坡度之关系曲线,如图1所示。图中显示图1泥石流泥沙体积浓度与溪床坡度之关系曲线Relationofsedimentconcentration,Cd,andchannelbedslope,tanθ,inthedebrisflow1.高桥与笔者之泥石流体积浓度曲线均呈向上扬升趋势,当溪床坡度略大于某数值时,其泥沙体积浓度即因快速增加而使之大于1.0,这与泥沙体积浓度的基本定义相互矛盾。2.欧国强泥沙体积浓度曲线系随着溪床坡度增加而渐趋平衡,即使在较大的溪床坡度下,泥沙体积浓度亦未超出合理的范围,惟因属一经验公式,其通用性仍有待检验。3.不论是泥石流先端部泥沙体积浓度或其整体平均泥沙体积浓度,从单一泥石流段波的整体流动特性来看,它们之间理应具有某种程度的共通特性,并可依相同或近似的机制建立其理论模式。惟目前系分别采用两种截然不同的机制及公式予以计算和模拟,使得先端部和整体平均泥沙体积浓度间的变化过程无法连贯,其机制自是无法统一地表达。因此,虞及泥沙体积浓度对泥石流数值模拟和实际运用上的重要性,本文拟通过一定的理论过程和渠槽试验资料,建立泥石流泥沙体积浓度之通用公式,不仅可以修正高桥公式在模拟泥石流先端部泥沙体积浓度上不尽合理的变化趋势,亦可运用于推估泥石流的整体平均泥沙体积浓度。2理论分析考虑一无限长且厚度均匀的无粘聚性饱和堆积土层,厚度为D,溪床面倾角为θ。当土层表面形成水深h0及泥沙体积浓度C0的水砂混合水流时,受重力作用之影响,水流不断地冲刷床面泥沙颗粒,使泥沙含量沿程增加,从而提高了流和力量。由于这种水流的行为机制系受沿程泥沙含量的增减所左右,不仅具有碰撞过程的应力转换效应,亦与溪床表面泥沙颗粒的冲刷互为相关。因此,同时考虑巴格诺流体内部固体泥沙颗粒间的本构关系和水力学溪床表面的泥沙冲刷机制,将水流中泥沙颗粒间的离散剪应力和剪阻力转化为净剪应力netshearstress,并与溪床表面泥沙颗粒的起动应力相结合,以决定水流对溪床表面泥沙的冲淤作用。这样,水流沿着流动方向上的净剪应力,可表为粒间离散剪应力与剪阻力之差,即τrτPdtanαρsρC0ρgh0sinθρsρgC0h0cosθtanαρgh0sinθ1ρsρ/ρ6C01tanα/tanθ图2无因次净剪应力与泥沙体积浓度之关系Relationofdimensionlessnetshearstressandsedimentconcentration图3表面流作用床面堆积土层模型图Thesimplifieddiagramofflowdepositlayersonthebedsurface由上式得知,水流净剪应力系随着泥沙体积浓度C0的变化而改变,且与溪床坡度和动摩擦系数相关。为此,将上式写为无因次形式,即Tτr/ρghosinθ1ρsρ/ρ1tanα/tanθ7如图2为在给定的动摩擦系数下,不同溪床坡度时无因次净剪应力与泥沙体积浓度之关系图。图中,无因次净剪应力T随着泥沙体积浓度C0之增加,分别呈现上升或下降,其趋势与溪床坡度和动摩擦系数间的大小互为相关当tanθtanα时,无因次净剪应力T呈缓和的上升趋势,而在tanθtanα的情形下,无因次净剪应力T呈快速的下降趋势。但是,不论水流净剪应力随着泥沙体积浓度呈上升或下降,来自于溪床表面泥沙颗粒之阻滞效应,水流净剪应力和床面泥沙颗粒间的互动机转中,总会朝着某种动态的平衡关系发展,而使床面泥沙颗粒不再被大量的冲刷或者淤积换句话说,在适当的床面条件配合下,水流泥沙含量终将趋于稳定,不再作激烈的增减,此时水流净剪应力和床面泥沙颗粒间会存在一种等势的关系,表为方程式,则可写为τrfcr8当水流和床面泥沙颗粒具有以上的关系时,水流与溪床间的泥沙交换将达到平衡,水流泥沙体积浓度趋于饱和,溪床泥沙既不冲刷也不淤积。式中fcr为溪床面泥沙颗粒起动之临界剪应力criticalshearstressforincipientmovementofsedimentparticlesonthechannelbed,可表为单位面积溪床面上泥沙颗粒水中自重沿着流动方向上的抗冲力,即fcrηρsρgdscosθ9式中ds床面泥沙特征粒径ξ综合因子,它是静摩擦系数、泥沙颗粒形状、级配、排列及其粘聚特性等因素综合影响的修正系数,需由实测资料决定之。根据式8得知,当水流净剪应力与床面泥沙颗粒起动之临界剪应力相等τrfcr时,因床沙与水流间的泥沙交换已趋平衡,水流中必然存在一饱和的泥沙体积浓度Cd,使其得以稳定均匀流条件向下游行进,故由式6和式9分别代入8式中,经整理可得,恒定均匀水流条件下之平衡泥沙体积浓度Cdρtanθ/ρsρtanαtanθη/tanαtanθds/h010上式适用于一般挟沙水流及泥石流之流况。式中,泥沙粒径与水深比值ds/h0为与水流之泥沙体积浓度相关,其关系式可透过一物理图形予以求解,如图3所示,为表面流作用床面堆积层之简单模型图。假设表面流从断面00至断面11的有限距离Δx中,可移动床面泥沙厚度平均达Δh,则由砂砾连续方程式可得C0U0h0CmUΔxCdU1h111式中床面泥沙冲刷速度UΔh/Δt,且ΔxU0U1Δt/2。为简化问题起见,假设入流速度U0与出流速度U1相接近时,即U0≈U1,又h1h0Δh,则上式经整理可得Δh/h0CdC0/CmCd12设表面流于断面00的泥沙体积浓度C0≈0同时,因床面泥沙移动厚度Δh至少必须大于或等于单一泥沙颗粒之特征粒径ds,即Δhnds,且n≥1,故上式可改写为ds/h01/nCd/CmCd13或者,选取图3断面11处单位底面积之柱体分析,当断面00泥沙体积浓度等于零时,则柱体全部泥沙量h1Cd应来自于床面堆积泥沙,故有h1CdΔhh0CdΔhCm14上式经化简后与式13完全一致,这显示式13在推导过程的部分假设尚称合理。因此,将式13代入式10中,经整理可得恒定均匀水流平衡泥沙体积浓度之理论表达式为15式中χ1ρtanθ/ρsρtanαtanθ16χ2η/n/tanαtanθξ/tanαtanθ17且tanθ≠tanα18同时,由于没有任何水流之泥沙含量,可以大于静止溪床面堆积土沙之最大泥沙体积浓度Cm值,故式15平衡泥沙体积浓度两个解的合理范围应介于之间。如图4所示,为于不同的系数ξ下,泥沙体积浓度随溪床坡度的变化状况。图中显示,当系数ξ越大时,因床面泥沙必须要在更强的水流条件下始能起动,故在相同的溪床坡度下,水流的泥沙体积浓度会较小。0≤CdCm19

注意事项

本文(水利工程论文-泥石流泥沙体积浓度之研究.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5