黄土高原地震作用下黄土滑坡滑距预测方法

王鼐;王兰民;王谦;王峻;王平

【摘 要】为评估黄土高原地区地震滑坡致灾区域,提出一种基于模糊信息优化处理的地震滑坡滑距预测方法。在野外调查、室内试验和分析的基础上,对黄土高原的地震滑坡类型、性质和影响因素进行分析。以摩根斯坦-普莱斯法计算黄土边坡的稳定性,建立边坡最小安全系数与影响因素的关系,其影响因素主要包括地震烈度、比高、坡角、容重、内黏聚力、内摩擦角等,并在此基础上得到纯黄土地震滑坡滑距的影响关系。将黄土高原地区数次大震中采集到的93个样本数据进行归纳分析,建立模糊信息优化处理模型。所得的计算结果通过误差校验和与其他滑坡滑距预测方法进行对比,来证明本模型有较高的准确性。最后对天水市22个潜在滑坡的边坡进行滑距预测。

【期刊名称】《地震工程学报》

【年(卷),期】2016(038)004

【总页数】8页(P533-540)

【关键词】滑动距离;地震滑坡;预测;黄土高原;模糊信息模型

【作 者】王鼐;王兰民;王谦;王峻;王平

【作者单位】[1]中国石油勘探开发研究院,北京100080;[2]中国地震局兰州地震研究所,甘肃兰州730000;[3]甘肃省地震局（中国地震局黄土地震工程重点实验室）,甘肃兰州730000

【正文语种】中 文

【中图分类】TU43

黄土主要分布于我国华北和西北黄土高原地区。黄土高原面积辽阔,土层深厚,地貌复杂,长期受内外应力作用,形成千沟万壑,地形支离破碎。黄土高原沿断裂带地震活跃,以山西隆起、六盘山和贺兰山三个地震带的地震较频繁[1],作为地震多发区,历史上曾发生8级以上地震7次,7～7.9级地震22次,6～6.9级67次以及5～5.9级地震258次。地震造成了逾140万人次的伤亡,而地震引发的滑坡则是造成重大人员伤亡和经济损失的主要原因。根据地震资料研究,黄土高原触发滑坡产生的最小地震震级为5级。由于岩性松软,斜坡临空面大,受地震惯性力的作用易丧失抗剪强度,故每当发生较强地震时,斜坡上不稳定的土体会在瞬时抛出并埋压村庄、城市,毁坏道路、渠道,摧毁工厂、矿山,堵塞江河、湖泊,破坏农田、水库等。地震滑坡造成的损失远远大于地震本身造成的直接损失。历史上发生的众多地震引起了不计其数的黄土滑坡。据统计[2],1654年天水8.0级地震滑坡分布面积约为1 000 km2,滑距大于500 m的滑坡59个;1718年通渭7.5级地震滑坡分布面积约为3 500 km2,滑距大于500 m的滑坡337个;1920年宁夏海原地震滑坡分布面积约为50 000 km2,形成滑坡657个。本文通过收集整理甘肃省境内历史上较大的黄土地震滑坡资料,提出基于模糊信息优化处理方法的黄土地震滑坡滑距的预测方法。

滑坡滑距的预测研究目前主要有三个方向:一是基于质点运动学研究滑坡运动过程。通过上世纪90年代对随机地震荷载的深入研究,建立随机地震荷载下非饱和黄土的动本构关系,结果表明,其本构关系依然服从双曲线模型,但模型参数随所施加的随机地震荷载时程的不同而不同。并初步建立了动本构模型参数与所施加的随机地震荷载卓越周期间的定量关系[3-6]。另外坡体波动震荡效应被用来解释地震滑坡的形成,地震动使坡体波动震荡效应产生的启程剧发速度直接影响滑行速度和滑移距离[7]。Kyoji Sassa滑坡运动模型假设运动的滑坡符合牛顿流体运动方程,在此基础上可以推导出二维滑坡运动模拟微分方程。这是一种简单适用的方法[8];另一种是通过设计的振动台试验。森宽协等[9-10]提出以能量守恒法则为基础的能量法,结合土体物性参数,建立滑坡滑距预测模型。章健[11]通过分析三类滑坡的运动轨迹方程,在此基础上建立运动学模型,并根据实验及实测参数计算每类黄土滑坡的滑距,最后分析利用运动学模型计算黄土滑坡滑距方法的准确性;再一种是利用统计学原理。刘悦等[12]采用模糊信息原理用少量的滑坡数据建立模型,认为其具有普适性。随后的研究中,加入了更多的影响因素,而统计模型越来越具有针对性[13]。针对滑坡滑距的预测方法计算公式见表1。

由于已知的滑坡分类众多,需对影响滑坡的众多因素进行分析和对比,确定一种有利于对黄土地震滑坡的特征及稳定性研究工作的分类方法,因此提出黄土地区滑坡类别。由于影响黄土地震滑坡的因素众多,在对各个因素进行分析后,通过对比各因素对滑坡影响的程度,并结合前人的研究[14-17],本文以斜坡的结构组成特点作为分类依据来对黄土地震滑坡进行研究分析。通过对147个黄土地震滑坡进行分析统计,将其主要分为四类:纯黄土滑坡、黄土-基岩接触面滑坡、黄土-基岩滑坡、黄土-砾石-基岩滑坡(表2)。

本文主要以纯黄土滑坡为研究对象。据统计,纯黄土地震滑坡在地震滑坡中占据了极

大部分。在147个滑坡中,黄土滑坡共计93个,占总滑坡数的63.3%,;黄土-基岩接触面滑坡9个,占总滑坡数的6%;黄土-基岩滑坡33个,占22.5%;黄土-砾石-基岩滑坡12个,占8.2%。

黄土斜坡的地震稳定性用安全系数来定量判定[18]。本文基于极限平衡分析原理中的摩根斯坦-普莱斯方法进行计算,包括所有的条件(法向力和剪切力)同时满足力矩平衡方法和所有静力平衡方程。

地震滑坡的发生、形成以及滑距受到多种影响因素的控制。这些因素包括烈度、比高、坡角、土层岩性等,是滑坡发生以及产生一定规模的充分条件。只有处于一定条件下,在地震发生时才有可能失稳产生滑坡。分析不同因素对于滑坡的影响规律,对河谷地区黄土地震滑坡的特征类型研究、数值模型建立及其稳定性分析有着重要的意义[19-20]。

(1)烈度对地震滑坡的影响

烈度是衡量一个地区受地震破坏程度的量。当地震发生时,不同烈度区内地震滑坡规模及滑坡单位面积产生数目不同。通过黄土区多次地震震害结果统计,当震级大于7级时地震影响的高烈度区内则会出现成规模的滑坡和滑坡密集区(图1)。

(2)坡角对地震滑坡的影响

坡角是控制滑坡发生的一个主要内在因素,是岩土稳定性的衡量指标之一,其大小反映了土体产生滑动的可能性大小。坡角决定了斜坡的形态,虽然理论上认为坡角越陡越容易发生滑坡,然而由于其他因素制衡,地震滑坡坡角主要分布于10°～25°(图2)。

(3)比高对地震滑坡的影响

滑坡的比高体现了地表的起伏变化,反应了岩土的势能差。理论上讲,比高越大越容易发生滑坡,但统计发现地震滑坡存在一个发生频率最高的范围区间。滑坡比高主要分布于30～400 m,当比高在50～250 m时为滑坡易发区域(图3)。

(4)地震滑坡受土层岩性影响规律

对斜坡稳定性产生影响的黄土物理力学特性参数主要包括内黏聚力、内摩擦角与容重。选取几个典型黄土滑坡土样做数值模型,在确定两个变量的同时对第三个变量进行规律性分析,可以看出,斜坡安全系数在黄土内黏聚力与内摩擦角变化范围内呈线性上升趋势,且随黄土容重的增加而降低(图4)。

根据图4中曲线呈现的趋势可以看出:斜坡安全系数随内摩擦角增加的上升趋势比随内黏聚力增加的趋势略微明显,即从单一因素角度分析,内摩擦角对安全系数的影响较内黏聚力影响更大。分析单一土的物性参数并不能充分研究其对斜坡安全系数的影响。对多因素进行拟合,将多个物性参数纳入计算范围才能更加准确地计算和分析河谷地区地震作用下黄土斜坡的安全系数。

不同学者根据自身研究的特点提出了各自对滑坡滑距的定义。针对滑坡的致灾距离进行研究[1],将纯黄土地震滑坡滑距定义为从滑体与滑面接触的后缘至滑体滑动最远地点之间的距离。如图5所示,其中L为滑坡滑距。

模糊信息优化处理模式针对地震滑坡取样的局限性以及数据的不完整性进行优化处理,

可以用来分析信息的模糊性,以便更好地帮助决策。因此,对地震作用下黄土滑坡致灾区域的研究采用模糊信息优化处理模式可以得到很好的解决。

本文拟采用的模糊信息处理方法原理如下所示:

设X是自变量论域或原因论域,Y是因变量论域或结果论域,在实用的近似推理模型中,一个论域总是由有限个离散点来表示。

令A0是一观测事件,X0是原因,Y0是结果,则A0=A0(x0,y0}。设R描述观测事件A,为刻画原因X和结果Y的因果规律,则有:

因此,当已知原因X0,可以通过合适的近似推理规则“|”推出Y0。

为了验证模糊信息优化处理模型的误差,计算一个简单的样本来与其他滑坡滑距预测方法进行对比。1718年通渭级地震造成了逾40 000人员伤亡。基于实地测量的数据,选取15个样本进行验证,其中10个作为计算数据,5个作为结果误差。具体内容如表4所列。

按照Y0=X0|R可以得到可能性分布Y0,“|”采用普通乘运算。表5为YOH利用信息集中处理求取最终结果。

同理可以算得烈度YOM和坡角YOα的滑距结果:

通过信息集中求取最终结果如下:

为验证方法的精度,将求得的最终预测滑距与实际滑距相比,求标准差:

为进一步验本模型的可靠性,分别采用经验统计公式、舒德格尔公式和森宽协公式(表1)进行计算,并将所得结果相比较。

国内经验公式:

舒德格尔公式:

森宽协公式:

将模糊数据优化处理模型公式预测得到的滑距与其他三种方法的结果进行对比,如图6所示。

由图6可知,模糊信息优化处理方法对于黄土地震滑坡的滑距估算要优于其他几种方法。这是因为其他方法针对滑坡滑距的估算并没有考虑到地震的作用,而在模糊信息优化处理模型中由于加入了烈度的参数,使模型对黄土地震滑坡的预测更加准确。各方法计算结果与实际滑距的误差验算也证明了这一点。

采用模糊信息优化处理模型进行地震滑坡致灾区域预测。

通过数据分析,95%滑距L变化范围为40～1 500 m,比高H为20～240 m,平均坡角α为7°～70°,烈度M为6～9万m3。故按照模糊信息处理有关理论,依次取论域如下:

VL={L1,L2,L3,…,L8}=

{40,140,240,340,440,540,640,740,840,940,1040,1140,1240,1340,1440,1540}

μα={α1,α2,α3,…,α7}=

{7,14,21,28,35,42,49,56,63,70}

其中:VL为地震滑坡滑距论域;μH、μα、μM分别为比高论域、平均坡角论域以及烈度论域。按照二维正态信息扩散公式分别对VL和μH、VL和μα、VL和μM建立模糊关系矩阵。其中m为样本数,m=93。通过二维正态信息扩散公式分别得到原始信息分布矩阵QHL、QαL和QML,如下可见:

将原始信息分布矩阵正规化处理可得各元素与滑距之间的模糊信息矩阵RHL、RαL和RML,由于篇幅有限,本文仅以烈度M为例(表6)。

将原始信息分布矩阵正规化处理可得烈度与滑距之间的模糊信息矩阵RML(表7)。

根据模糊近似推理公式Y=X|R,得到天水市22个滑坡滑距烈度的可能性分布:

{0.03,1.00,0.08,0.30,0.42,0.38,0.67,0.07,0.01,0.00,0.00,0.00,0.00,0.00,0.00,0.00}

根据信息集中公式求取最终结果:

可以求得滑坡滑距如表8所列。

分析结果可知,天水市区所测斜坡在地震作用下,当地震烈度大于Ⅷ度强时会发生滑坡,滑动距离从20～150 m不等。值得注意的是低速缓动的斜坡,滑距计算结果反而相对其他滑坡滑距较大,这也符合了低角度滑坡滑速缓、滑距远、致灾区域大的特点。

本文观测的147个黄土滑坡中纯黄土层滑坡共计93个,占总滑坡数的63.3%。采用模糊信息优化处理方法对这类滑坡进行地震诱发致灾距离预测研究,得出以下结论:

(1)通过对滑坡的特性研究,将滑坡类型分为纯黄土层滑坡、黄土-基岩接触面滑坡、黄土-基岩滑坡和黄土-砾石-基岩滑坡。

(2)采用摩根斯坦-普莱斯法进行斜坡的地震稳定性计算。通过对纯黄土地震滑坡滑距的不同影响因素进行分析,发现影响斜坡稳定性和滑坡滑距的主要因素为烈度、坡角、比高等。

(3)通过对比模糊信息优化处理方法计算结果和其他方法结果,以计算误差为量化数据,证明该方法具有较高的准确性。

随着对地震滑坡的深入研究,数字化和定量化在研究过程中占据了越来越重要的地位。因此,建立地震信息数据库和滑坡信息数据库,将二者科学地联合在一起并及时地共享数据,对于黄土地区地震滑坡的预防和监测具有十分重要的意义。