泥石流动力学模型及数值模拟研究难点问题.docxVIP

泥石流动力学模型及数值模拟研究难点问题 1 泥石流动力学研究的难点 作为一种典型的灾难，泥石流形成有许多影响因素和复杂的动力学过程。这是尚未通过的最后一个问题。同时由于全球气候变化、山区经济发展的加快,以及泥石流对生命财产造成的严重破坏经常发生在远离源区的运移路径中等因素,泥石流灾害带来的风险日益加剧。例如2004年全国因地灾害造成死亡的有743人,125人失踪,直接经济损失48.7亿元。而仅仅2004年7月的云南滇西地区怒江、宝山、德宏3个州市的两次洪涝、群发泥石流灾害,就造成48人死亡,85人失踪,伤病497人,直接经济损失10.3亿元。 目前对于改善人类生存条件的计划,不论是与资源环境有关的计划,还是与灾害防御有关的计划,很大程度上都要依靠旨在增进人类对地球作用过程认识的基础研究的成果。只有对泥石流等地质灾害的动力学物理过程有充分的了解,并在GIS等新技术的支持下进行客观、准确地模拟预测,才能更有效地制定防灾减灾策略,减少不必要的人员伤亡和财产损失。因此,泥石流等环境地质灾害的动力学研究日益受到国内外包括政府部门和学术机构的重视。 在过去的几十年的时间内,许多国内外学者进行了旨在更好理解泥石流物理机制的理论探索和实践应用,如美国地质调查局(USGS)成立了专门的泥石流研究机构,他们提出的计算统计模型已经在一些地区得到实际的应用。杜榕桓、吴积善、康志成等学者在中国泥石流形成机理方面进行了开创性研究工作,并在云南小江流域建立了东川泥石流观测研究站,积累了大量宝贵的观测资料。王光谦等人在泥石流的本构模型、模拟问题上进行了较深入的研究,促进了国内泥石流本构模型研究的发展。但是根据大量的文献资料和实例研究,目前泥石流动力学研究仍然存在一些未能很好解决的难点问题,主要包括: (1)已有泥石流动力学本构模型的不完善和适用性问题; (2)已有数值模拟方法之间存在分歧,难以选择应用; (3)泥石流动力学复杂环境(地形、水文等)效应问题远未得到很好的解决,如:复杂环境下的泥石流应力分布、动能传递的定量分析;复杂环境下的泥石流孔隙水压力重分布和传播的定量分析等。 (4)泥石流动力学数值模型耦合及与GIS集成等。 由于以上问题的存在,严重影响了泥石流动力学模拟结果的准确性和实用性。本文将对这几个难点问题进行深入剖析,系统分析其有效解决和改进的途径,实现泥石流动力学过程的准确模拟预测。以下从4个方面对泥石流本构模型和动力学模拟的研究现状和存在问题进行具体阐述。 2 深度均化理论—泥石流动力学本构模型及适用性 泥石流动力学过程分析包括泥石流影响范围、程度的评估及与灾害有关的各种参数分布规律的研究,包括运移速度、流体深度、流体分布、堆积等。泥石流动力学研究有利于进行泥石流潜在破坏性的估计、防护设施的设计、确定和预测灾害的后期效应,如滑坡泥石流灾害造成的涌浪或洪水现象等。 虽然泥石流体是以重力为主要驱动力的流体,但其动力学机制非常复杂。仅就固体颗粒之间的相互作用而言,就包括:碰撞、摩擦、滑动、转动等多种形式。固体颗粒之间的流体压力作用的考虑大大增加了本构模型研究的复杂性。许多泥石流体还往往携带大量的树木、有机质等非岩土类的介质,这更增加了泥石流研究的难度。 泥石流动力学本构模型的研究是动力学过程模拟的基础,它描述泥石流体介质的应力应变关系、流体流动行为及流变特征。 早期的本构模型比较有代表性的是Johnson的粘塑模型。粘塑模型简单地将泥石流体看作单相的宾汉(Bingham)或库仑(Coulomb)均质体。在低剪切应力作用下,泥石流体近似看作弹性体;在高剪切应力下作为粘性流体。有关的流变参数如:屈服强度、内聚力和宾汉粘性系数等通过现场勘测获取。但粘塑模型不能精确地描述影响泥石流体流动和堆积的某些重要动力学特征:如颗粒与颗粒,颗粒与水体之间的相互作用。 为了描述这种颗粒与颗粒之间的相互作用,Takahashi 在Bagnold工作的基础上提出了惯性粒子流模型,并在1992年将模型进行了一定的改进,增加了侵蚀相的估计。该模型根据流体的剪切速率计算悬浮在水体中固体颗粒的正应力与剪切应力相关关系的变化。虽然Takahashi的模型考虑了固体颗粒之间的相互作用,但它假设固体颗粒在水体中均匀分布,不允许孔隙水压力发生变化,这与野外和实验的结果有很大的分歧。 Savage和Hutter对颗粒流的动力学特征进行了较为深入的研究,主要包括:(1)推导了一维浅层流的动量和质量守恒方程,采用库仑准则解决了流体底部的剪切力和正应力分布的问题。(2)采用Lagrangian有限差分算法对方程进行了数值求解。(3)进行模型实验,验证模型的准确性。后来许多学者在Savage-Hutter模型的基础上将该模型推广到多维空间,并将孔隙间粘性流体的作用纳入到模型中,产生了一些具有代表性