作者简介:张大林(1987-),男,山东济南人,博士研究生,主要从事地质灾害评估和预测预报研究.E-mail:zdl87@aliyun.com
崩岗泥砂流是降雨过程中在崩岗流域内形成的一种固液两相流,是崩岗侵蚀泥沙向外输移的主要方式,泥砂流流体类型的判别也是崩岗治理的依据之一。通过野外考察与现场采样,对崩岗泥砂流的粒度组成特性进行了分析。结果表明,泥砂流容重介于1.16~1.60 t/m3之间,含沙量为257.03~960.55 kg/m3,且均自沟道上部至下部逐渐降低。泥砂流浆体以粉砂和黏粒物质为主。随着容重的增加,粒度曲线由单峰型转变为与风化壳土体类似的双峰型,呈现无分选搬运的特点,且流体粒度逐渐粗化。泥砂流固体物质中值粒径与流体容重有较好的线性正相关关系。通过对泥砂流与泥石流和黄土丘陵沟壑区高含沙水流粒度特性的对比后发现,泥砂流属于高含沙水流向泥石流过渡的中间类型,但与泥石流具有更为密切的关系,可以认为泥砂流是广义泥石流的一个亚类,即崩岗型泥石流。
Mud sand flow is a kind of solid-liquid two-phase flow formed in collapsing hill and gully basin during rainfall. It is the main way to export erosion product. The discrimination of its fluid type is one of the collapsing hill and gully control theoretical basis. This paper analyzed the basic characteristics of mud sand flow like grain size and so on through fieldwork and sampling. The results show that the density of mud sand flow is between 1.16~1.60 t/m3 and the solids content is between 257.03~960.55 kg/m3, both of which decrease from the upper to the lower channel. The slurry of mud sand flow is composed mainly of silt and clay. As the density increases, the particle size distribution curve transforms from a single peak to the bimodal distribution similar to the weathering crust with no sorting, and the grain size of mud sand flow becomes coarser which shows a well positive linear correlation between the sediment median particle diameter and density of mud sand flow. The comparison during mud sand flow, hyperconcentrated flow and debris flow shows that collapsing hill and gully mud sand flow, which belongs to an intermediate class between hyperconcentrated flow and debris flow, has a closer link with debris flow. Therefore, mud sand flow can be considered as a sub-class of generalized debris flow that may be called as collapsing hill and gully type debris flow.
崩岗即崩塌的山岗,在地貌分类中属于花岗岩丘陵地貌,如同岩溶地貌、黄土地貌和丹霞地貌一样,同属于岩石地貌大类中的一种[ 1]。由于崩岗的发生发展与流水作用形成的沟谷地貌密切相关,崩岗亦曾归属为暴流地貌[ 2]。虽然暴流地貌这一名称目前已很少使用,但崩岗这一地貌学命名却被保留下来。崩岗这一称谓源自于广东梅州五华县,当地客家人将“丘陵山地冲沟源头汇水区围椅状崩塌崖壁地貌”称为崩岗。崩岗主要分布在我国南方花岗岩丘陵山区,全国崩岗总侵蚀面积约5万km2,面积≥60m2的崩岗总数达24万处,其中广东约11万处,江西近5万处,福建、湖南、广西均超过2万处[ 3]。20世纪20年代广东德庆县才最早发现有崩岗出现,至今也不过百年历史[ 1]。作为地表过程,崩岗研究具有综合性、动态性以及方法的多样性等特点[ 4],我国自20世纪60年代开始,已积累了许多研究成果[ 5, 6, 7],主要涉及崩岗的分布规律、形态特征、形成发展、开发治理等。近些年随着研究的深入和新技术方法的应用,崩岗研究领域不断拓宽,朝着多学科融合的趋势发展,如数值模拟、遥感监测、水动力学的应用等[ 8, 9, 10],为崩岗研究提供了新的思路。崩岗侵蚀模数巨大,危害严重,降雨过程常激发固体物质含量高的沟道流体,是崩岗侵蚀产物向外输移的主要途径及破坏方式。许炯心[ 11]从高含沙水流角度关注了此类现象,并初步研究了其发生条件、形成过程和泥沙输移特征。沈灿燊[ 12]运用地理学与水力学有关理论,研究了花岗岩地区崩岗地表片流和径流作用以及泥沙的冲刷和搬运作用,建立了华南崩岗地区水土流失产流产沙预报概念模型。林金石等[ 13]运用组合指纹法研究崩岗泥沙迁移的来源,发现崩岗冲积扇的组成物质主要来源于风化壳的砂土层及碎屑层。
刘希林等[ 14]从崩岗土体物理性质入手,通过与泥石流源地土体以及沟道堆积物的对比研究认为,崩岗具备泥石流产生的物质和地形条件,在降雨作用下形成的沟道流体,是泥石流大类中的一个新种——泥砂流。泥石流是介于土体运动与水流运动之间的一种土、水混合流,以此与一般挟沙水流相区分[ 15]。在实际工作中,根据不同的需求,对泥石流可有多种分类方式。例如,按泥石流动力条件,可将其划分为土力类和水力类泥石流。按地貌形态,可分为河谷型、沟谷型和山坡型泥石流。此外还可根据启动模式、暴发规模、发育阶段等多种方式进行泥石流的类型划分[ 16, 17, 18]。根据泥石流流体性质进行分类,能够体现泥石流的本质特性,通常采用泥石流的粒度特征及其容重来进行划分,亦有狭义和广义之分。狭义的泥石流包括黏性、亚黏性(过渡性)和稀性泥石流,广义的泥石流还包括泥流和水石流[ 19]。本文所指的崩岗泥砂流即是介于高含沙水流与狭义泥石流之间的过渡流体。崩岗泥砂流过程复杂,对其流体类型的判别是其内部属性和外部特征的综合体现,是在原有研究基础上,通过泥砂流实体样本的分析,对这种新的崩岗型泥石流——泥砂流的进一步剖析。
莲塘岗崩岗位于广东五华县华城镇NW方向约4 km处、省道S228桩号K67附近,地理坐标N24°06′10.5″,E15°34′57″,属亚热带季风气候,多年平均降水量1 547.5 mm,76%集中在4~9月,多暴雨。莲塘岗丘陵全部由风化物质组成,母岩为黑云母粗粒花岗岩,全区未见基岩出露,风化壳厚度达50m以上,机械组成主要为黏粒含量较高的含砾粗砂土,干容重在1.40g/cm3左右。三维激光扫描仪定位监测显示,莲塘岗崩岗为典型的瓢形崩岗,主沟长104 m,自沟口至崩壁顶部相对高差55 m,侵蚀面积5 125 m2,内部地形陡峭,≥25°的陡坡面积占总侵蚀面积的80.73%,其中≥60°近乎直立的崩壁面积占16.69%[ 9],崩壁、崩积锥、沟道、洪积扇(沙渍地)等崩岗地貌类型齐全,目前正处于发育旺盛的壮年期阶段[ 20](图1)。2011年6月至2012年5月间,侵蚀量高达1 007 m3/a,侵蚀模数达269 268 t/km2·a。每年雨季,有规模大小不一、历时长短不等的泥砂流发生(图2),是研究泥砂流比较理想的场所。
2013年5月21日和7月14日,在莲塘岗崩岗现场,遭遇了突然的降雨事件,目睹了崩岗泥砂流的形成、发展和堆积过程。2013年5月21日上午8时,五华县华城镇地区突降中到大雨,莲塘岗崩岗流域内产生泥砂流,规模较大,持续1个多小时,降雨停止后不久,泥砂流即行结束。2013年7月14日下午2时,莲塘岗崩岗及其周边地区开始降雨,雨量为小到中雨,约20 min后,泥砂流开始形成,规模不大,降雨持续约20 min,泥砂流全过程约30 min。2次泥砂流与降雨过程相关,表现出与降雨共生共灭的特点。
2013年5月21日采集了4个泥砂流样品,编号1~4号,其中1~3号样品为泥砂流产生后同一时间沟道内不同部位的样品。1号样品采集于崩岗崩积锥底部支沟沟头部位。2号样品采集于崩岗支沟汇合部位,即崩口位置。3号样品采集于崩岗主沟出口处。4号样品为泥砂流前锋样品,在崩口附近采集。2013年7月14日采集的4个泥砂流样品编号5~8号,分别为崩岗主沟道内同一时间自上游至下游采集的实时流体样品。
2次实地共取得8个泥砂流样品,带回实验室后,测量体积与质量并计算浆体容重,之后予以沉淀过夜,抽取表层清水后,在105℃烘箱温度下进行烘干称重,计算样品的固体物质含量。烘干后的固体颗粒进行粒径分析,粒径>2 mm砾石用筛析法测定颗粒级配,≤2 mm的泥沙颗粒用Mastersizer 2000激光粒度仪进行测定。
莲塘岗崩岗泥砂流容重及含沙量见表1。由表1可知,泥砂流容重介于1.16~1.60 t/m³之间,平均为1.29 t/m³,前锋部分的容重达到了稀性泥石流的状态(容重1.60 t/m³),含沙量介于257.03~960.55 kg/m³之间,平均为465.52 kg/m³。如此看来,泥砂流固体物质含量虽然部分达到或接近于稀性泥石流,但整体上尚未完全达到稀性泥石流的容重要求;同样,泥砂流类似于高含沙水流,但又不完全等同于高含沙水流(文中讨论部分将详述)。因此,泥砂流是一种介于稀性泥石流与高含沙水流之间的新型固液两相流。
同一时间支沟沟头、支沟汇合处、主沟出口处取得的1,2,3号样品的容重和含沙量逐渐减小,中值粒径亦由0.128 mm下降至0.012 mm,固体物质颗粒逐渐细化。这种现象可解释为:一是降雨产生的地表径流并未足够大,水流只能搬运2 mm以下的砂粒物质,因而>2 mm的砾粒物质在泥砂流流动过程中沿途沉积下来,从而导致泥砂中值粒径自沟道上游向下游逐渐变小,颗粒越来越细,容重也相应变小;二是降雨期间并非每条支沟都形成泥砂流,没有形成泥砂流的支沟以清水汇入主沟,导致泥砂流在流动过程中被稀释,容重逐渐变小、含沙量不断降低(图3)。
4号样品为泥砂流前锋样品,容重和含沙量分别为1.60 t/m³和960.55 kg/m³。泥砂流前锋即泥砂流形成后的第一波流体,能量集中,砾石聚集,冲刷和侵蚀能力强,泥砂流到来之前沟床粗糙,摩擦阻力大,泥砂流流动时拖曳和挟裹沟床物质,使其固体物质含量不断增大,容重逐渐接近稀性泥石流,从而表现出与稀性泥石流体相似的特性。泥砂流前锋浆体黏稠,砾石含量高(图4),冲击力较大、破坏力也相对较大,向下游推移过程中不断扩散并最终停积。
莲塘岗崩岗泥砂流粒度成分见表2,并与莲塘岗崩岗风化壳土体作对比分析[ 14]。除2,3号样品外,其它泥砂流样品的砾粒含量均较高,平均为32.3%,主要为2~5 mm的细砾部分,砾粒含量比崩岗风化壳土体样品的平均值26.87%高出5.43%。莲塘岗崩岗风化壳土体以含砾粗砂土为特征,砂粒(包括粉砂)平均含量为60.97%,其中粗砂平均含量为23.60%,粉砂为24.15%。泥砂流各样品砂粒平均含量为51.69%,其中粗砂仅为6.60%,粉砂高达35.31%。两者中砂和细沙含量差别较小。表明粗砂是泥砂流水力分选的主要沉积区段,粉砂则是泥砂流主要的组成物质。泥砂流黏粒含量平均高达16.02%,高于莲塘岗崩岗风化壳土体样品的12.17%。黏粒含量高,表明泥砂流浆体黏稠,但由于缺乏粗大砾石,泥砂流的性质与泥石流并不完全相同。样品中值粒径差异较大,4号样品的中值粒径为2 mm,3号样品的中值粒径仅为0.012mm。全部样品的平均中值粒径为0.419 mm,低于莲塘岗崩岗各类土体样品的平均中值粒径0.525 mm,表明泥砂流经过水流的短距离搬运,确实存在一定的分选作用。
泥石流中一般以2 mm粒径为分界,小于2 mm的颗粒作为浆体物质,是搬运介质,大于2 mm的颗粒则是被搬运的固体物质[ 18]。激光粒度仪测定的≤2 mm泥砂流粒度分布曲线表明(图5),泥砂流固体物质颗粒级配主要以双峰为主,但不同样品峰值粒径略有不同,高峰值显著集中在粒径0.01 mm附近,印证了泥砂流以粉砂(粒径0.005~0.075 mm)为基本组成的特点。除2,3号样品为单峰以外,其余6个样品在粒径1 mm附件的粗砂段(粒径0.5~2 mm)都具有第二个峰值,构成双峰型粒度曲线,其中崩岗风化壳土体以及容重较高的4号(泥砂流前锋)和5号样品(主沟道上部)的粗砂段峰值波形陡峭,明显高于粉砂段。泥砂流固体物质的组成既受崩岗原始土体成分的影响,也随泥砂流的输移过程而改变,随着泥砂流含沙量和容重的增加,流体的分选作用下降,重力流作用增强,其物质组成接近于崩岗风化壳土体,泥砂流粒度分布曲线由单峰型转变为双峰型,更多地表现出与泥石流相似的特性。
泥砂流中值粒径与容重具有较好的正相关关系(图6),泥砂流中值粒径越大,其容重越高;反之亦然。当无法实测泥砂流流体容重时,可参考图6中的回归方程式,通过泥砂流暴发后实地取样分析,以其平均中值粒径反算,也不失为得到泥砂流流体容重估算值的一种经验方法。
突发性短历时降水或堵塞—溃决作用可使径流量激增,形成高浓度泥砂流,如4号和5号样品,其固体物质粒度曲线与崩岗原始土体较为接近,分选相对较差。由于高容重的流体和黏稠的浆体对沟道底部具有较大的冲击力和拖曳力,产生的“揭底”作用吸纳了沟底的砾粒物质,4号样品粒径>2 mm的砾粒含量高达50.9%,明显高于崩岗原始土体。总的来说,崩岗泥砂流的侵蚀和堆积表现为“小水淤,大水冲”的特点。泥砂流以冲淤交替方式将崩岗土体物质逐渐搬运出沟口,在地貌发育“削高填低”的过程中,虽然崩岗绝对高度随着崩岗的发育进程而降低,但由于泥砂流的存在,在“削高”的同时又能够侵蚀沟底,这样维持崩岗的相对高差基本不变,从而促使崩岗持续发展并不断向丘陵内部掘进(图7)。
莲塘岗崩岗可分为崩壁集水区、崩积锥分布区和沟道流通区(图8)。流通区以外则是丘陵谷地堆积区,表现为洪积扇或沙渍地。崩壁是指丘陵坡地遭受崩岗侵蚀后上部出露的风化壳剖面,坡度较陡近乎直立,以重力侵蚀为主。崩壁后退发生在雨季,与降雨过程密切相关,以崩塌方式进行,因此具有一定的突发性和随机性。崩壁物质结构紧密,坡度较大,降水往往来不及下渗而快速流走,因此崩壁部位在泥砂流形成过程中主要起上游集水的作用。崩积锥是崩壁重力崩塌的产物,物质组成松散,坡度可达30°以上,上部有纹沟和细沟发育,深度从十几厘米到几十厘米不等,长度可达10 m以上。崩积锥地带的侵蚀作用以雨滴溅蚀和片流冲刷为主,崩积锥之间的地带则发育股状水流,细沟进一步发展成切沟,沟道横断面呈“V”字形,深度可达数米,强烈的下切及侧蚀造成沟壁垮塌,使得沟间地带呈现陡脊形态(图9),沟道水流在崩积锥分布区获得大量泥砂补给,是泥砂流主要固体物质来源地。崩岗支沟在崩口处汇合成主沟,地形变缓,泥砂流在此出现冲淤交替现象,随后沿沟道进入流通区,地形平缓,坡度在10°左右,泥砂开始沉积。流通区以外为开阔的堆积区,由于缺少沟道束缚,泥砂流迅速扩散沉积,形成洪积扇或沙渍地。莲塘岗崩岗主沟道较短,约47 m,沟口与相邻崩岗主沟道交汇(图3c),因而洪积扇不明显,表现为沿谷底分布的沙渍地。
崩岗泥砂流与泥石流均为降雨过程中暴发的阶段性洪流。典型的泥石流沟谷主要包括清水汇流区、形成区、流通段和堆积区。崩岗地形与泥石流沟谷在形态上基本相似,但崩岗是水力—重力复合侵蚀,两种侵蚀作用按照不同时空组合交替进行[ 20]。从时间顺序而言,水力是崩岗发育的第一驱动力,重力则是在水力侵蚀造成临空面后发挥作用,然后两者交替进行。从空间顺序来看,上部崩壁集水区以重力作用为主,中部崩积锥分布区重力和水力作用同时进行,下部沟道流通区则以水力作用为主。同一崩岗内两种外营力的对比关系可发生转换,很大程度上影响了泥砂流的形成,导致其性质和物质组成随空间和时间产生较大变化。泥石流作用营力相对单一,例如土力类泥石流主要依靠势能重力沿坡面的剪切分力发生运动,水力类泥石流则是主要依靠水力动能和冲刷作用维持运动。泥砂流与泥石流的差别还在于固体物质含量及其组成结构不同。研究表明[ 14],崩岗沟道中泥砂流堆积土体的平均中值粒径仅为泥石流沟道堆积土体的1/12,砾石含量仅为泥石流的1/4。受制于化学风化充分、基岩埋藏较深、以及发育规模相对较小等因素,崩岗泥砂流搬运物质中大于5mm的粗砾物质很少,特别是较大的石块和漂砾更是缺乏。泥石流中的粗大颗粒很多,较大的石块和漂砾一般直径在1~2 m,大者可达5~7m,因此具有很强的冲击力和破坏力[ 18]。
高含沙水流是指挟带一定黏性细颗粒、具有高浓度固体颗粒的一种复杂挟沙水流,也是一种典型而特殊的固液两相流[ 21],主要分布于我国西北的广大黄土高原地区。钱宁[ 22]曾提出高含沙水流按照物质组成的不同具有两种极端模式,一种是以细粉砂及黏土等细粒物质为主的高含沙水流,另一种是以细沙以上颗粒为主的高含沙水流。第二种模式或者两者间的过渡形态多为两相高含沙紊流,常出现于黄河干支流,具有常年流水的属性,因此两相高含沙紊流与崩岗泥砂流的区别显而易见。接近第一种模式的高含沙水流具有非牛顿体的性质,多发生在黄土丘陵沟壑区,亦即通常表达的黄土 “泥流”[ 23]。综合相关参考文献[ 18, 19, 23]对此类高含沙水流以及泥石流基本特征的描述,与泥砂流进行比较,三者的主要特性见表3。
从流体容重角来看,高含沙水流与崩岗泥砂流相差不大,一般均低于亚黏性泥石流1.85 t/m3的容重界限值,与稀性泥石流容重存在重叠。由于物质基础的不同,三者在流体物质组成方面具有明显差别。泥石流可发育于比降较大的各类岩性地区,主要由黏粒、砂粒、砾石及漂砾等多种物质构成,具有极宽的颗粒级配。高含沙水流受黄土粒径级配限制,主要由黏粒和粉砂组成,中值粒径小于0.025mm,大于2mm颗粒百分含量小于2%。崩岗泥砂流的物质组成介于两者之间,其发育的花岗岩风化壳地区以含砾粗砂土为主,砂粒是其主要组成部分,虽含有大量2mm以上的未风化石英颗粒,但相对于泥石流仍缺乏大的砾石和漂砾。高含沙水流呈单峰曲线,峰值在0.035mm左右,崩岗泥砂流与泥石流都具有单峰、双峰两种形式,当容重较高时,流体处于无分选状态下搬运,流体中固体物质组成具有双峰特征;当容重下降,流体以细粒物质组成的浆液为主,浆体在有分选状态下形成粒径的单峰分布。总体来看,崩岗泥砂流在地貌形态、形成过程及性质等方面与泥石流具有较多相似之处,但与高含沙水流和黏性泥石流又存在明显区别,从物质组成分析来看,应属于两者之间的过渡形态。因此,从广义泥石流角度,可将崩岗泥砂流作为泥石流的一个亚类,即崩岗型泥石流。
莲塘岗崩岗泥砂流容重介于1.16~1.60 t/m³之间,平均为1.29 t/m³,含沙量为257.03~960.55 kg/m³,平均为465.52 kg/m³。突发性短历时降水或堵塞—溃决作用可形成高浓度前锋阵性流,拖曳和“揭底”作用造成砾石聚集,可使泥砂流容重和含沙量达到极大值,超过稀性泥石流的容重1.46 g/cm3的标准。
泥砂流粒度分布曲线兼具单峰和双峰形态,峰值分别位于0.01mm的粉砂段和1mm的粗砂段。随着容重增加,粒度分布曲线逐渐由单峰型转变为双峰型,水力分选作用下降,物质组成向崩岗风化壳土体接近。粗砂是泥砂流水力分选的主要沉积区段,粉砂是泥砂流浆体的主要组成物质。泥砂流容重越高,物质组成越粗,其中值粒径与容重具有很好的线性关系,符合经验公式r =0.1946d50 +1.2073 (式中r为泥砂流容重,单位t/m³;d50为泥砂流沟道堆积物中值粒径,单位mm)。
崩岗流域可分为崩壁集水区、崩积锥分布区、沟道流通区和丘陵谷地堆积区。崩积锥分布区地形坡度较大,松散固体物质丰富,是泥砂流形成的主要固体物质来源地,区内细沟、切沟十分发育,降雨期间的片流侵蚀和股流侵蚀,极易启动泥砂流。在水力和重力共同作用下,泥砂流不断发展壮大,最终表现出冲淤交替以及侵蚀、搬运和堆积的完整地貌过程。
泥砂流与泥石流及高含沙水流既有联系又有区别,其物质组成是三者的主要区别所在。泥砂流属于高含沙水流向泥石流过渡的一种中间类型,但其地貌形态、形成过程和流体粒度特性与泥石流更为类似。泥石流的分类有多种,根据组成物质来分类,曾有过泥流、泥石流、水石流的分类方案,因此,我们认为可以将崩岗泥砂流作为广义泥石流的一个亚类,即崩岗型泥石流。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|