大连市地质灾害评价研究

大连市地处辽东半岛最南端，属于中国环渤海经济圈北部前沿，大连市的经济在最近几年发展迅速，城市人口增多，大连市在发展经济的同时也没有忽视地质灾害的预防和治理，大连市存在的主要地质灾害类型有泥石流、地面崩塌、山体滑坡以及塌陷等。就现阶段对大连市地质进行勘查的情况来看，大连市存在着不少的安全隐患。基于此，对一个城市所可能遭遇的地质灾害进行风险评价，进而预测其可能发生的地点以及损害程度，可以帮助城市进一步的制定好应对地质灾害的策略，利于实现人与自然的和谐共处，同时帮助政府在预防自然灾害方面提供一定的科学依据。

1 研究区概况

大连的位置在 120°58′至 123°31′，北纬 38°43′至40°10′之间，大连市位于北半球的暖温带地区，是东北地区最温暖的地方。年平均气温10.5℃，极端气温夏季最高37.8℃，冬季最低-19.13℃。大连市的年降水量550-950毫米，60%-70%的降水集中于夏季，多暴雨，且夜雨多于日雨。大连市河流主要涉及两大水系，分别是黄海流域、渤海流域，水流流向黄海的河流包括英那河、碧流河、庄河、大沙河、赞子河、清水河、登沙河以及马栏河等，水流流向渤海的河流包括复州河、三十里堡河以及李官村河等。

大连地区位于千山山脉西南延伸部分，辽东半岛的最南端，濒临黄海、渤海之间，形成碧海环抱、低山丘陵起伏的地形。地形起伏较小，平原海拔在20-60米之间，地形坡度主要在0-10°之间，面积约占全市总面积的77%，剩下为低山丘陵区，坡度起伏较小在10-20°之间，坡度最大的是53°，分布在步云山区和老帽山区，低山丘陵占全市面积的16.69%。大连山脉是NNE向太古界鞍山群变质基底变质深成侵入岩为主体，青白口积累震旦系古生界寒武系和奥陶系，侏罗纪和白垩纪沉积中生代和古生代的第四世纪。太古宙鞍山群主要分布在锦州东部、福州西部、沙河、瓦房店也暴露出驼山。元古代辽河群主要由变质火山沉积岩经历了低绿片岩和大理岩角闪岩相变质的麻粒岩相变质，主要分布在瓦房店地区。新元古代地层是本区最发达的地层。他们分布广泛，充分暴露。该平台的特点是碎屑岩盆地和克拉通盆地碳酸盐岩储量大，典型的沉积特征。由于后期构造运动和局部变质地层，南部的片岩、千枚岩、板岩，而北福州是弱变质或不变质。层序地层学是连续的、完整的，形成良好的接触，清晰的接触，和厚度大，共有7501 9317m沉积厚度[1]。

大连市大地构造位置隶属于中朝准地台胶辽台隆的复州台陷，Ⅳ级构造单元包括城子坦断块和复州—大连凹陷，构造形态复杂，具有地台型二元结构特征。

太古代晶基底组成了区内东部城子坦断块，其褶皱形态单斜化，片麻理走向北东—东西向，向南倾斜。上元古界与古生界盖层，褶皱复杂，按褶皱轴迹展布可划分为近东西向、北东向和北西向及平卧褶皱。

综上所述，高龄剖宫产术后孕产妇产后出血发生率较高，而影响其产后出血的危险因素众多，临床工作中可依据其危险因素为高龄剖宫产术后孕产妇实施合理防治措施。

本区断裂构造发育，可分为脆性断裂和韧性断裂两种。韧性断裂（韧性剪切带）主要有大黑山韧性剪切带、尹家村韧性剪切带和信台子韧性剪切带。大连地区是地震的多发区，大连及邻区有9条比较大的与地震活动有关的活动断裂，即营维断裂带（郯庐断裂北中段）、金州断裂、普兰店断裂、大黑山断裂（大和尚山山前断裂）、庄河断裂、海城河断裂、碧流河断裂、皮口断裂和北海断裂带。

2 地质灾害特征

2.1 地质灾害类型

4.1 通过搜集的资料研究得知，大连市境内地质灾害共有772处，主要有四种类型，分别为崩塌474处，滑坡108处，泥石流154处，地面塌陷34处。

  

图1 大连市地质灾害类型饼状图

2.2 地质灾害分布规律

崩塌是大连的主要地质灾害，各个地区都有分布,由于修建道路、景区、居民楼、厂房等人为因素开挖削坡，形成坡度大于60°的陡崖，认为是形成临空高耸的地形条件。

4.3 确定了地质灾害危险性评价指标，包括灾点密度、坡度、地层、地质构造、河流、降水、地震、人类工程活动八个指标。运用层次分析法的各项指标的权重分别为0.0936、0.1521、0.0178、0.0341、0.0757、0.3979、0.0489、0.1800。

地质灾害时间上分布主要集中于夏季，地质灾害的诱发因素主要是降雨，大连市的降水主要集中在夏季的6-9月，由此所引发的地质灾害也集中在降雨多且集中的月份。

  

图2 大连市地质灾害空间分布

 

表1 地质灾害分布（按行政区）统计表

  

名称 面积（km2）人口（万人）人口密度（人 /km2）行政乡、镇、街道（个）崩塌（处）滑坡（处）泥石流（处）地面塌陷（处）小计(处)密度（处/km2）中山区西岗区沙河口区甘井子区旅顺口区金州区长海县瓦房店市普兰店市庄河市43.85 26.6 48.32 502 506.8 1390 119 3794 2896 4073 0.410 0.602 0.228 0.104 0.037 0.054 0.546 0.054 0.060 0.034合计 160 474 108 156 34 772 0.058 33.95 30.57 69.31 132.18 32.48 110.28 7.80 94.22 74.12 84.13 7742.30 11492.48 14343.96 2633.07 640.88 793.38 655.46 248.34 255.94 206.56 8791 5 8 10 13 24 5 30 23 26 12 7 47 15 53 57 96 95 84 42541 581 194 0 00200007 8 00000702 65 11 052 18 16 11 52 19 75 65 205 173 137

地理信息系统（GIS）能够存储、显示和分析地理数据。主要包括数据录入与编辑、数据操作、数据管理、数据显示和输出等。空间数据作为获取、处理、管理和分析空间数据的重要工具、技术和学科，得到了广泛的关注和迅速的发展。

层次分析法的基本原理是基于目标的层次结构、子目标、约束来评价方案、选取两两比较特征向量矩阵，然后对应于相应的最大权系数矩阵。

3 地质灾害的危险性评价

在地质灾害危险性评估中，影响评价结果准确性、科学性的关键是评价指标的选择以及权重确定的合理性，目前对权重的确定方法主要有：专家评分法、模糊综合评判法、AHP层次分析法、调查统计法、数理统计法、信息模型法、隶属函数法等[2]。在这些定权方法中，专家评分法、调查统计法、均属于定性分析方法，调查统计法是在专家打分法基础上发展起来的，其实质尚未脱离定性分析，数理统计需要更多涉及到的数据较多，要经过较为复杂的处理，涉及到大量的计算过程[3]。层次分析法主要将复杂的研究对象来当成一个整体，决策问题分解为不同的层次，然后构造判断矩阵最大值的相邻层之间，特征向量和相应的水平，总加权指数权重相对目标水平。我国是在1982年开始引入层次分析法，它已被广泛的应用在社会的各个方面。它是灵活的、简单的、需要较少的定量数据。这是一个系统的分析方法[4]。因此，本文采用层次分析法确定各指标的权重。

3.1 层次分析法的理论与方法

层次分析法（AHP）是一种定性和定量相结合的决策分析法。它是一个由决策者对复杂系统的决策过程进行建模和量化的过程。使用这种方法，决策者将一个复杂的问题分成若干层次和若干因素，与各种因素相比，操作简单，可以得到不同方案的权重，为程序优化提供依据。

2.2.4 混播比例对混播植物Pro含量的影响 不同混播比例对禾草Pro含量均无显著影响(表3)，而对苜蓿Pro含量有明显影响(表4)。B1处理Pro含量最高，比B2，B3和B4处理高24.48%，27.41%和49.50%，差异极显著(P<0.01)；单播苜蓿苜蓿Pro含量分别比B3和B4处理高17.77%(P<0.05)和38.19%(P<0.01)，B2处理Pro含量比B4处理高20.10%(P<0.05)，其余处理间苜蓿Pro含量无显著差异(P>0.05)。

云南少数民族地区农村劳动力外出务工意愿的Order Logisitic回归模型，通过了极大似然法的卡方显著性检验，说明因变量与解释变量的概率关系显著，具体概率关系见表2。

建立层次结构。在这一次步骤中，要求将问题所含的因素进行分组，把每一组作为一个层次，按照最高层（目标层）、若干中间层（准则层）以及最低层（方案层）的形式排列起来。如果某一个元素与下一层的所有元素均有联系，则称这个元素与下一层次存在有完全层次关系；如果某一元素与下一层的部分元素有联系，则称这一元素与下一层次存在有不完全层次关系。层次之间可以建立子层次，子层次从属于主层次中的一元素，它的元素与下一层次的元素有联系，但不能形成独立层次。

构造判断矩阵。这一个步骤是层次分析方法的一个关键步骤。判断矩阵表示针对上一个层次中的某元素而言，评定该层次中各有关元素相对重要性的状况。设有n个指标{A1,A2,...,An}，aij表示Ai相对于Aj的重要程度判断值。aij一般取 1，3，5，7，9等 5 个等级标度,对于判断矩阵算出来的各要素的相对权重进行一致性检验，当结果小于0.1时，表示判断矩阵具有一致性，当结果大于0.1时则需对判断矩阵实施调整，一直达到满意的一致性。

3.2 大连市地质灾害危险性分区评价

根据对大连市的地质灾害影响因素的分析，将地质灾害危险性分区的评价指标体系分为灾害历史、基础因子和诱发因子，取灾点密度、坡度、地层、地质构造、河流、地震、大气降水、人类工程活动八个评价因子，如图3。

云南省财政厅相关负责人表示，截至今年10月底，云南省纳入财政部PPP综合信息平台项目437个，总投资10856亿元，示范项目86个，投资额2983亿元，项目个数和投资额均位居全国第一位。下一步，省财政厅将加强PPP项目规划、方案设计、规范实施和督导检查，强化绩效管理，加快落地实施，更加规范、有序地开展PPP工作，发挥好PPP项目拉动投资、弥补短板、提质增效、防范风险的作用。

由上述模型可构造判断矩阵，判断矩阵中各要素的数值由多名专业的专家分别对各要素的相对重要性进行评估打分，然后综合各个专家的智慧，得到判断矩阵，并进行一致性检验。结果如表2-表6所示。

1996年，适逢团场机车改制，在很多机车驾驶员犹豫时，孙改会不管不顾，硬是靠东挪西凑的钱，在连队带头把他那辆开惯了的铁牛55拖拉机买了下来。同时，他和妻子承包的棉花地也增加到了100余亩。尽管孙改会领着妻子每天忙得跟个陀螺似地，但一年年不断递增的经济收入让他把想要创业、创新的视野和心胸放得更长、更远。

  

图3 地质灾害危险性分区的指标体系

研究区内地质灾害的成因与地形地貌、岩土结构、气象因素、人类工程活动有直接关系。其中人类活动是地质灾害的主要原因，而降雨是主要的诱发因素。

在上述指标权重计算和各要素数据归一化的基础上，通过ArcGIS 10.0的统计分析功能，统计每个要素的指标值，得出再利用叠加功能，将每一个归一化后的要素图层与层次分析法的权重相结合进行叠加，如图4所示。

通过各个要素图层的叠加，如表7，得到地质灾害危险性程度评价结果，如图5，图6所示。

4.2 研究区内地质灾害主要受新构造运动、地形坡度、降水、地质构造、人类工程活动等因素的影响。

4 结论

大连市位于辽东半岛南部，地质地形条件复杂，地质灾害总类多，影响范围大，通过对大连市地质灾害危险性评价的研究，得到以下几点结论。

普洱茶产业的壮大需要将普洱茶目前自然摆放的陈化过程发展成为一道专业的加工工艺，关于普洱茶专业仓储醇化工艺的相关研究对未来普洱茶产业的提升意义重大。从上述国内普洱茶仓储贮藏和普洱茶陈化工艺的研究现状来看，普洱茶在不同的环境条件下进行一定时间内的醇化贮藏，普洱茶的品质将产生很大的变化，所以寻找既能加快普洱茶醇化速度，又能提高醇化品质的方式或方法是普洱茶专业仓储醇化研究的根本目的。目前云南省已有多家分别位于昆明、临沧、元阳等地的专业普洱茶仓储企业开展了相关醇化研究。

大连市的地质灾害类型多种多样，主要包含有崩塌、滑坡、泥石流以及地面塌陷等，其中以崩塌为主。经排查，现有地质灾害隐患点772处，包括崩塌474处，滑坡108处，泥石流156处，地面塌陷34处，这几类地质灾害分别占地质灾害总数的61.4%、14%、19.9%、4.4%,如图 1所示。

  

图4 要素图层在ArcGIS平台下运算过程

 

表2 危险性区划判断矩阵标度层次排序表

  

一致性比例：0.0825，对评价指标体系的权重1.0000，λmax：3.0858

 

?

 

表3 灾害历史判断矩阵标度层次排序表

  

一致性比例：0.0000，对评价指标体系的权重0.0936，λmax：1.0000

 

?

 

表5 诱发因子判断矩阵标度层次排序表

  

一致性比例：0.0904，对评价指标体系的权重0.6267，λmax：3.0940

 

?

 

表4 主控因子判断矩阵标度层次排序表

  

一致性比例：0.0742，对评价指标体系的权重0.2797，λmax：4.1981

 

?

3.3 地质灾害危险区分区评价依据地质灾害危险程度叠加结果，地质灾害危险程度划分为高危险区、中危险区、低危险区三个级别，如图7所示。

大连的崩塌灾害主要分布在各学校、住宅小区、公矿企业、道路主干道、旅游景点等人工活动强烈的地区，开挖的坡脚坡度大，容易发生崩塌灾害。滑坡主要分布在甘井子、金州、长海县、瓦房店、普兰店庄河等地的低山丘陵地区。泥石流主要分布在瓦房店与普兰店交界的老帽山地区，庄河北部山区的步云山镇、蓉花山镇、仙人洞镇、塔岭镇等山区。地面塌陷主要分布在金普新区的石棉矿采空区、董家沟煤矿采空区、瓦房店李官镇华铜矿采空区、庄河市荷花山芙蓉铜矿采空区、普兰店双塔—星台一带的花岗岩矿采空区等各种矿区形成的采空塌陷，岩溶塌陷主要分布在瓦房店市三家子岩溶区，金州北山、大连湾十三里台、金石滩等地的十三里台组岩溶地区，灾害分布如图2和表1所示。

4.4 大连市地质灾害危险分区划分三个等级，即高危险区、中危险区、低危险区。其中高危险区面积2311.4209km2，占大连陆地面积的18.13%；中危险区面积7210.96km2，占大连陆地面积的56.57%；低危险区面积3224.45km2，占大连陆地面积的25.3%。

 

表6 判断矩阵层次总排序

  

评价指标体系 灾害历史 主控因子 诱发因子 总排序权值 排序中间层—方案层灾点密度坡度地质构造河流地层降水人类工程活动地震0.0936 1.0000 0.2797 0.5439 0.2706 0.1219 0.0636 0.6267 0.6348 0.2872 0.0780 0.0936 0.1521 0.0757 0.0341 0.0178 0.3979 0.1800 0.0489 43578126

 

表7 评价因子的选取及指标量化

  

序号 指标 权重 分层图像 说明1 灾点密度指标 0.0936 2 坡度指标 0.1521 122°0＇ 123°0＇E 40°0＇N 39°0＇122°0＇ 123°0＇E 40°0＇N 39°0＇已有地质灾害群体统计评价指标应包括滑坡、崩塌（以及不稳定斜坡）、泥石流自然地质现象的数量和规模。此次调查仅采用已有的崩塌、滑坡、泥石流、地面塌陷的数量指标计算单元格内的数量密度。利用GIS从DEM数据中分别提取调查区的坡度信息，然后进行归一化。将斜坡的危险程度按照不同坡度区间地质灾害、自然地质现象发生的概率，进行0-1之间的归一化，得到坡度归一化结果。3 地层指标 0.0178 122°0＇ 123°0＇E 40°0＇N 39°0＇此次调查，将地层岩石结构对滑坡、崩塌、泥石流、地面塌陷危险性程度的影响进行0-1之间归一化差值处理。122°0＇ 123°0＇E 40°0＇N 4 河流指标 0.0341 39°0＇河流上游沟谷是地质灾害易发区，危险性程度高，将全市河流进行0-1之间归一化后得到河流归一化结果。122°0＇ 123°0＇E 40°0＇N 5 地质构造 0.0757 39°0＇按照不同断裂发育对区间自然地质现象发生概率的影响，进行0-1之间的线性归一化，得到地质构造归一化结果。122°0＇ 123°0＇E 40°0＇N 6 地震指标 0.0489 39°0＇按照地峰值加速度对地质灾害的影响程度，进行0-1之间归一化处理，得到地震指标归一化结果。

 

续表7 评价因子的选取及指标量化

  

序号 指标 权重 分层图像 说明122°0＇ 123°0＇E 0 25 50 100千米 N 40°0＇N 7 降水指标 0.3979 39°0＇对本区多年平均降水量等值线进行0-1之间线性归一化，得到降水归一化结果。灾害指标高：1低：0 122°0＇ 123°0＇E 40°0＇N 0 25 50 100千米 N 8人类工程活动指标0.1800 39°0＇交通指标高：1低：0人类工程活动对地质灾害形成发育的影响是极为复杂的，如何定量化反映是个难题。考虑到公路铁路等交通建设是区内最具代表性的人类工程活动，对灾害影响最明显，且具有贯穿或覆盖全区的特点，本次人类工程活动的量化是以调查区的公路为基准线，经栅格化和归一化处理，参与评价。

  

图5 地质灾害危险程度信息叠加

  

图6 计算结果直方图

  

图7 大连市地质灾害危险性分区图

4.5 基于危险性评价结果，提出相应的防治对策。高危险区，应当采取搬迁避让的措施，在中危险区和低危险区采取工程措施进行治理，加强群策群防与政府宣传，并建立地质灾害的相关数据库，成立地质灾害信息管理系统，形成三维立体的监管体系，能够快速和有效的实现对灾害的预警和避让。

通过搜索旅游地信息可以较快地得到反馈结果，省事且高效。所以，在以往的发展中，人们将注意力都放在了网络搜索上，忽视了论坛、问答社区等渠道的发展。

参考文献：

[1]王国良.层次分析法在地质灾害危险性评估中的应用[J].西部探矿工程，2006（09）：286-288.

临上警车，刘雁衡察觉有人牵他衣角，回头一看，是黄莺。黄莺眉头紧锁，眼中有泪。刘雁衡低声说：“放心吧，不会有事的。”

[2]辽宁省地质矿产局.辽宁省区域地质志[M].地质出版社，1989.

[3]李晚良，王春耕，孙熙芝，等.大连泉水—南关岭岩溶地下水库可行性勘察报告 [R].大连：辽宁水文地质工程地质勘察院，1993.

纵观中华传统文化的内涵广泛，传承是个急需思考的问题。像单田芳先生的评书，其实就是传承传统文化的一种形式，可以从儿时听到老时，它是精彩的，它是丰富的，在我们的精神方向上，它犹如一盏明灯，指引我们前行。所以人们在单田芳先生离开时，才会有着那一份敬意，那一份感激。

[4]《中国水利百科全书》第二版编辑委员会，中国水利水电出社.《中国水利百科全书》第二版.北京：中国水利水电出版社，2006：576.