石旭飞,江 山,郭晓东,赵海卿
中国地质调查局沈阳地质调查中心,辽宁 沈阳 110000
长(长春)吉(吉林)图(图们)经济区是我国东北地区重要的工业基地,是振兴东北老工业基地国家战略的重要组成部分[1].但是,在城市建设与发展过程中,各种自然与人为诱发的崩塌、滑坡、泥石流、地面塌陷等地质灾害[2-4]严重阻碍了长吉图经济区经济社会的发展[5-6].因此,如何保持人类经济活动在无害少灾的环境中生存发展并保持与地质环境较为协调的关系,是摆在区域规划者、城市建设者、环境地质研究者面前一个难题.
解决上述问题的一个重要方法是对研究区地质灾害易发性进行区划.目前国内外研究地质灾害易发性区划的方法主要有GIS法[7-10]、数理统计法[11-15]、数学模型法[16-19].由于地质灾害系统的复杂性和不确定性[20-21],上述方法均是基于地质灾害的影响因素和地质现象作出宏观评价,有关地质灾害易发性区划定量研究还较为缺乏.
综上,本研究在调查分析长吉图经济区内存在的崩塌、滑坡、泥石流、地面塌陷等环境地质问题的数量及其影响因素的基础上,利用地质灾害综合危险性指数法对上述地质灾害的易发性进行定量研究并分区,为研究区优化配置利用地质资源与国土空间规划提供科学依据.
长吉图经济区位于吉林省中东部,包括吉林省长春市城区、德惠市、九台市和农安县,吉林市城区、蛟河市和永吉县,及延边朝鲜族自治州全境,面积约72 000 km2.研究区属温带大陆性季风气候,降水多集中于6—9月份,占全年降水量的80%以上.雨季多暴雨,常造成江河泛滥,诱发地质灾害.区内河网密布,以长白山群峰为中心,向西北和东北分别形成松花江、图们江二大水系(图1).
图1 长吉图经济区地貌分区及多年平均降水量等值线图Fig.1 Geomorphic division and multiyear average precipitation isoline in the CJT Economic Zone
研究区地势东南高西北低,起伏较大.全区自东向西以北东向为主的山岭、盆地平行排列,相间分布.研究区基岩主要有以奥陶系和侏罗系泥岩、砂质泥岩、粉砂岩互层为主要特征的泥质岩类以及变质的中基性-酸性岩浆岩、碳酸盐岩和三叠系、白垩系夹煤层的沉积碎屑岩.研究区第四系分布广泛,以河流冲积、冲洪积和冰川堆积为主.中更新统黄土状亚砂土、亚黏土比较发育.
研究区地下水资源丰富,在山区地下水主要赋存于基岩裂隙中,而在平原区及山间盆地、河谷区地下水主要赋存于松散堆积物的孔隙中,富水性差别大.
地质灾害易发区指具备地质灾害发生的地质构造、地形地貌和气候条件,或可能发生地质灾害的区域.作者在分析研究区内主要环境地质问题的特征与影响因素的基础上,利用原国土资源部下发的《<县(市)地质灾害调查与区划基本要求>实施细则》[22]中的地质灾害综合危险性指数法对研究区地质灾害易发区进行划分,进而分析研究区环境地质问题的分布特征.
将研究区进行网格剖分,每个单元面积为1 km×l km~3 km×3 km.对于地质条件变化不大的地区,单元面积可取高限;对于地质条件复杂的地区,单元面积可取低限.计算每个单元格的地质灾害综合危险性指数Z:
式中:Zq为潜在地质灾害强度指数;r1为潜在地质灾害强度权值;Zx为现状地质灾害强度指数;r2为现状地质灾害强度权值.
潜在地质灾害强度指数按以下公式计算:
式中,Ti分别为控制评价单元地质灾害形成的地质条件(D)、地形地貌条件(X)、气候植被条件(Q)、人为条件(R)充分程度的表度分值,各评价指标的选取与评判标准依据具体情况而定;Ai分别为各条件的权值,根据实际情况分配.
现状地质灾害强度指数可以用单元格灾害点密度、灾害面积密度以及灾害体积密度求得:
式中,a为归一化处理后的灾害个数密度系数;b为归一化处理后的灾害面积密度系数;c为归一化处理后的灾害体积密度系数,地面塌陷灾害体积密度系数为0.
根据各评价单元的地质、地形地貌、气候以及人类工程活动等条件,利用MapGIS软件空间分析功能计算各评价单元的潜在地质灾害强度指数与现状地质灾害强度指数,赋权后计算得到各评价单元的地质灾害综合危险性指数,并根据计算结果合并相同单元格,划定地质灾害易发区.
崩塌是长吉图经济区突发性地质灾害中发生频度最大的灾害类型,区内一般包括岩崩和土崩两种类型.
研究区岩崩多分布在坡角受人为破坏、岩体裸露、坡度较陡的山坡处,尤其是江河沿岸,具有崩落滚动特点,水平位移较小.区内土崩坡角基本处于天然状态,坡度40~45°左右,多发生于表层松散物质发育、植被稀少的山体坡角处.崩塌时具有一定的水平位移,面积较大,厚度较小.单体崩塌具有多次重复发生性,表现为某一崩塌带(群)一年内或多年内多次重复发生,有时每年发生2~5次,如长白山天池地区.
研究区崩塌的区域空间分布规律主要表现为在某些地区成群、成带、连片地集中分布,绝大多数分布在中低山区,其次为低山丘陵区,尤其是集中分布在老岭中山区及图们江沿岸,均发育在地形高差大(一般高差大于300 m)、坡度陡(坡角40°以上)的地区.区内崩塌点岩土体类型以岩浆岩为主,占72.8%;碎屑岩占15.9%;变质岩占11.2%.此与区内岩体类型分布面积比例基本一致.从发生时间上看,发生于6—9月的崩塌占研究区全年崩塌的87.4%,其中7—8月占65.6%,说明崩塌与降水关系密切,具有明显的正相关性.
泥石流是研究区东部中低山区分布较普遍、活动较频繁的地质灾害,主要以河谷型水石流和沟谷型泥石流为主.区内引发泥石流地质灾害的降雨临界值约为:24 h降水量大于100 mm,或48 h连续降水量大于150 mm.研究区大型泥石流周期为10~20年发生一次,小型泥石流灾害发生周期为2~3年发生一次.
研究区泥石流绝大部分分布在张广才岭-龙岗山的中低山区.河谷型水石流分布在山间坡降较大、地形开阔,物源广泛的河谷中,其形成区与流通区无明显界线,搬运路程较长,具有流速快、能量大、破坏力强的特点.沟谷型泥石流分布于中低山山区易于汇水、地形坡度大、冲沟发育、无常年流水、碎屑物质丰富的沟谷中,物源范围小,具有发生频率高、规模小的特点.
研究区泥石流多以稀性泥石流为主,一般均发育多条泥石流支沟,坡降20°左右,常成群集中分布.从泥石流中的固体碎屑物质组成上看,河谷型由碎石、砾石、砂及黏土组成,粒径大小相差不太悬殊,有一定的磨圆度和分选性;沟谷型由块石、角砾、砂、黏土及耕植土组成,粒径大小相差十分悬殊,磨圆度差,且无分选性.泥石流流体运动特征也有差别.河谷型泥石流属稀性泥石流,固体碎屑物质在流体中具有明显的垂直分带性,土粒与砂粒呈悬浮状态,紊流运动,被搬运距离远;碎石、砾石在流体中为半悬浮状态,在流体下部以翻滚、跳跃和被推移的搬运,被搬运距离相对较近.沟谷型泥石流属黏性,固体碎屑物质含量高,其运动除靠流体搬运外,由于坡度大,自身还具有较大的下滑力.因此,固体碎屑物质在流体中没有明显的垂直分带性,块石、碎石、砂、土混杂.研究区泥石流的堆积特征主要有:河谷型泥石流堆积区比例大,自上而下略有粒径由大到小的总体趋势,多次活动后,堆积扇轴方向变化较小,堆积物叠加为主;沟谷型泥石流堆积区比例很小,粒径分布无明显规律,多次活动后,堆积扇扇轴摆动大,堆积物相接,部分叠加方式分布.
坡面型泥石流主要分布在老岭山脉,一般分布在山体坡角40~50°的坡面稍凹处,其特点是流通区极短,供给区与堆积区几乎相连.堆积物近锥体型,规模一般较小,泥石流沟总长度百余米左右,供给区长度50~80 m,堆积区长30~60 m,堆积物无分选,暴发时间短.
研究区内滑坡分布较少,但危害较大,多属岩体蠕动滑坡,具有间歇性活动特点,已发生的滑坡均为小型.研究区内滑坡按动力成因多属自然滑坡,按滑坡体物质组成划分多属土体滑坡,滑动面为碎石土滑坡体与基岩接触面.区内滑坡具有较明显的形态特征一致性,即滑坡要素基本相同.滑坡体上部由破碎松散的土石组成,表面起伏不平,周界不明显;滑坡壁高度1~5 m左右,陡度60~70°左右;滑动面(带)倾向与坡体总体坡度基本一致,与滑动方向基本一致的剪切裂缝成群出现,一般延伸长度20~30 m,个别达50 m左右,多数宽度为5~10 cm,最宽达30 cm,相对错距一般0.3~1.0 m,最大为1.5 m.
地面塌陷是研究区造成单位面积国民经济损失最大的地质灾害,多为采矿引起.多年以来,各矿山均有不同程度的地面塌陷发生,尤其是珲春市、九台市、蛟河市等煤矿分布较多的城市.研究区地面塌陷坑一般为近椭圆形或近圆形,比采空区范围稍大,塌陷深度从内到外,由大变小,有时数个塌陷坑相连或部分重叠.
通过利用地质灾害综合危险性指数法对研究区地质灾害易发性进行计算,将研究区环境地质问题分为4个高易发区、2个中易发区和3个低易发区(图2).
图2 长吉图经济区环境地质问题分区图Fig.2 Environmental geological zoning map of the CJT Economic Zone
主要分布于地质构造密集,岩土体类型复杂,年降雨量较大的长白中低山区,分布面积16 048 km2,占研究区总面积的22.31%.区内有崩塌369处,滑坡61处,泥石流285处,地裂缝及地面塌陷11处,共计726处,占研究区地质灾害点总数的46.75%,平均地质灾害点密度为4.52个/100 km2.
(1)和龙-龙井-延吉-珲春高易发区(I1)
位于长吉图经济区的东南部的和龙市、龙井市、延吉市及珲春市和西北部的的安图县部分地区,分布面积11 751 km2,占研究区总面积的16.34%.区内岩土体类型主要为坚硬的中厚层状以砂砾岩为主的岩组、软弱的中层—薄层状以砂砾岩为主的岩组及坚硬块状以花岗岩为主的侵入岩岩组,北西向及北东向断裂发育.区内地貌由长白火山与熔岩台地区、延边中山低山区及延吉、和龙、图们、珲春丘陵盆地等组成.区内有崩塌336处,泥石流152处,滑坡58处、地裂缝及地面塌陷7处,共计553处,占研究区地质灾害点总数的35.61%,平均密度为4.71个/100 km2,为崩塌、滑坡、泥石流高易发区域.
(2)嘎呀河上游高易发区(I2)
位于汪清县嘎呀河上游沿岸,分布面积2 772 km2,占研究区总面积的3.85%.区内岩土体类型主要为砂砾岩及软弱的中厚—薄层状黏土,为泥石流的形成提供了丰富的物源条件.区内地貌类型主要为中低山.区内有泥石流76处,崩塌14处,滑坡2处,共计92处,占研究区地质灾害点总数的5.92%,平均密度为3.32个/100 km2,为泥石流高易发区.
(3)安图县松江-永庆高易发区(I3)
位于研究区东南部安图县所属三道乡、松江镇、两江镇、小沙河乡及永庆乡,分布面积1 406 km2,占研究区总面积的1.95%.区内岩土体类型以较坚硬夹软弱的中厚—薄层状砂砾岩黏土岩互层岩组为主,地貌类型主要为中低山.区内有泥石流55处,崩塌14处,地裂缝及地面塌陷4处,共计73处,占研究区地质灾害点总数的4.70%,平均密度为5.19个/100 km2,为泥石流高易发区域.
(4)长白山天池高易发区(I4)
位于长白山天池周边,分布面积119 km2,占研究区总面积的0.17%.区内岩土体类型主要为块状坚硬玄武岩,地貌属长白山火山中山区,地形坡度大.崩塌、泥石流地质灾害发生频率高,现有崩塌5处,泥石流2处,滑坡1处,总计8处,占研究区地质灾害点总数的0.52%,平均密度为6.72个/100 km2,为崩塌、泥石流、滑坡高易发区域.
主要分布在岩土体类型较多或地质构造较密集的东部低山丘陵及部分中低山地区,分布面积34 599 km2,占研究区总面积48.10%.区内有崩塌364处,滑坡63处,泥石流179处,地面塌陷及地裂缝65处,共计671处,占研究区地质灾害点总数的43.21%,平均密度为1.94个/100 km2.
(1)吉林-蛟河-敦化-汪清中易发区(II1)
位于研究区中部及东部的蛟河市、永吉县、敦化市、汪清及安图县的部分地区,分布面积32 595 km2,占研究区总面积的45.31%.区内岩土体类型复杂,变质岩、花岗岩、碎屑岩、土体均有分布,地貌类型大部分为低山丘陵,部分为中低山.区内有崩塌313处,滑坡62处,泥石流173处,地面塌陷及地裂缝53处,共计601处,占研究区地质灾害点总数的38.70%,平均密度为1.84个/100 km2.区内受人类工程活动影响比较强,岩土体较为松散,在公路沿线及岩土体类型接触带附近发生崩塌灾害较多,加之年平均降雨量较大,在汛期常发生泥石流等突发性地质灾害,为崩塌、泥石流中易发区域.
(2)九台-上河湾中易发区(II2)
位于九台市,呈带状分布,面积2 004 km2,占研究区总面积的2.79%.区内岩土体类型为砂砾岩、火山岩及砂类土,地貌类型为低山丘陵及冲洪积台地,地形坡度10~40°.区内煤矿、采石场较多,人类工程活动比较剧烈,地面塌陷灾害较为发育.区内有崩塌51处,滑坡1处,泥石流6处,地面塌陷及地裂缝12处,共计70处,占研究区地质灾害点总数的4.51%,平均密度为2.79个/100 km2.由于区内九台煤矿开采并且修建铁路、公路较多,因此引发地质灾害发生较多,属崩塌、地面塌陷为主的中等易发区域.
主要分布在岩土体类型较少,地质构造密度较低的地区以及平原区,分布面积21 288 km2,占研究区总面积的29.59%.区内有崩塌123处,滑坡10处,泥石流23处,地面塌陷及地裂缝1处,共计157处,占研究区地质灾害点总数的10.11%,平均密度为0.74个/100 km2.该区域大部地区地形比较平缓,地质构造较不发育,丘陵区森林植被覆盖较好,不易发生危害性较大的地质灾害.
(1)双阳-万昌-孤店子低易发区(III1)
位于伊舒槽地、伊通县及吉林市,分布面积5 270 km2,占研究区总面积的7.33%.区内有崩塌28处,滑坡3处,泥石流11处,地面塌陷及地裂缝1处,共计43处,占研究区地质灾害点总数的2.77%,平均密度为0.82个/100 km2.区内岩土体类型较简单,以砂类土和黏性土为主.该区域地质灾害以土质崩塌、水土流失为主,主要分布在河流沿岸,其表现形式多为人工切坡或河流冲刷形成的陡坎而导致的崩塌.
(2)林海-海浪低易发区(III2)
位于蛟河市与敦化市北部蛟河、牡丹江的上游,分布面积4 672 km2,占研究区总面积的6.49%.区内有崩塌28处,滑坡6处,泥石流12处,共计46处,占研究区地质灾害点总数的2.96%,平均密度为0.98个/100 km2.区内岩土体类型主要有坚硬块状中厚花岗岩、薄层状软弱花岗岩风化壳、较坚硬块状-厚层状火山碎屑岩,为崩塌低易发区域.
(3)平原地区为主地质灾害低易发区(III3)
主要分布在长春-德惠-农安一带,地形比较平坦,分布面积11 346 km2,占研究区总面积的15.77%.区内共有67处小型崩塌点,1处滑坡点,共计68处,占研究区地质灾害点总数的4.38%,平均密度为0.60个/100 km2.区内岩土体类型较简单,主要为黄土状亚黏土、亚砂土、砂、泥岩等.区内地质灾害以土质崩塌为主,多分布在河流沿岸,其表现形式多为人工切坡或河流冲刷形成的陡坎.
(1)研究区内崩塌、滑坡、泥石流地质灾害发生的主要影响因素是降水、地形地貌、岩土体类型、地质构造、人类工程活动.区内地质灾害的发生多由以上数种因素共同作用产生,其中降水是诱发地质灾害的最主要的原因,其次为地形地貌.
(2)研究区内的和龙市、龙井市、延吉市、珲春市、安图县部分地区、汪清县嘎呀河上游沿岸以及长白山天池周边等地为崩塌、滑坡、泥石流等地质灾害高易发区,共有地质灾害点726处,占研究区地质灾害点总数的46.75%;蛟河市、永吉县、敦化市、汪清县等低山丘陵地区为崩塌、滑坡、泥石流地质灾害中易发区,共有地质灾害点671处,占研究区地质灾害点总数的43.21%;伊通县、吉林市和长春-德惠-农安一带为崩塌、滑坡、泥石流地质灾害低易发区,共有地质灾害点157处,占研究区地质灾害点总数的10.11%.
(3)研究区内崩塌、滑坡、泥石流等地质灾害高易发区分布面积约占研究区总面积的22.31%,平均密度为4.52个/100 km2;中易发区分布面积约占研究区总面积48.10%,平均密度为1.94个/100 km2;低易发区分布面积约占研究区总面积29.59%,平均密度为0.74个/100 km2.