临汾地裂缝灾害与地下水开采相关关系.pdf

中国矿业大学学报990 12 2 中国矿业大学学报 JO U RNA L O F CH I NA U NI VERSI T Y O F M I NI NG T ECH NO LO G Y 1999年 第1期 第2 8 卷 Vo l . 2 8 No . 1 1999 临汾地裂缝灾害与地下水开采相关关系 董东林 武 强 孙桂敏 姜振泉 隋旺华 田宝霖 摘要 山西临汾市地裂缝比较集中在5个区域发育，其灾害造成了巨大的经济损失，严重影响着临汾市的经济发 展. 阐述了地下水开采对地裂缝灾害的影响，并通过地面沉降的计算机的模拟，对地下水超采引起的地面沉降及 临汾市地裂缝灾害的发生与发展的时空关系进行了详尽的探讨，为该市的城市规划部门提供了科学的依据. 关键词 地裂缝灾害，地下水超采，沉降计算，相关关系，计算模型 中图分类号 P 6 4 Co r r e l a t i v e Re l a t i o n Be t w e e n Ea r t h -Fi s s u r e H a z a r d s a n d O v e r -Ex p l o i t a t i o n o f G r o u n d w a t e r D o n g D o n g l i n W u Q i a n g Su n G u i m i n ‧ Ji a n g Zh e n q u a n Su i W a n g h u a T i a n Ba o l i n ‧ D e p a r t m e n t o f Re s o u r c e Ex p l o i t a t i o n En g i n e e r i n g , CU M T , Be i j i n g 10 0 0 8 3 Co l l e g e o f M i n e r a l Re s o u r c e a n d En v i r o n m e n t Sc i e n c e s , CU M T , Xu z h o u 2 2 10 0 8 A b s t r a c t T h e r e a r e f i v e d e v e l o p m e n t z o n e s o f e a r t h -f i s s u r e s i n Li n f e n Ci t y i n Sh a n x i Pr o v i n c e . T h e h a z a r d r e s u l t s i n e n o r m o u s e c o n o m i c l o s s e s a n d h i n d e r s e c o n o m i c d e v e l o p m e n t o f t h e c i t y s e r i o u s l y . I n t h i s p a p e r , t h e e f f e c t o f o v e r - d e v e l o p m e n t o f g r o u n d w a t e r o n e a r t h -f i s s u r e i s s t u d i e d . T h r o u g h t h e s i m u l a t e d r e s u l t s o f s u r f a c e s u b s i d e n c e , t h e s p a c e -t i m e r e l a t i o n b e t w e e n t h e s u b s i d e n c e c a u s e d b y o v e r -d e v e l o p m e n t g r o u n d w a t e r a n d e a r t h -f i s s u r e h a z a r d s i s r e s e a r c h e d c o m p l e t e l y . I t i s a s c i e n t i f i c b a s e f o r u r b a n p l a n n i n g . K e y w o r d s e a r t h -f i s s u r e h a z a r d , o v e r -d e v e l o p m e n t o f g r o u n d w a t e r , s u b s i d e n c e c a l c u l a t i o n , c o r r e l a t e d r e l a t i o n , c a l c u l a t e d m o d e l 临汾市位于山西省南部临汾盆地内，是晋南经济、文化的中心城市. 随着临汾市经济的发展，城市规模不断 增大，地下水的需求量也逐年增加. 由于多年来该市地下水处于无序的开采状态，地下水位也急剧下降，从而导 致地面的不均匀下沉. 再加之该地区频繁的、高强度的地震活动及新构造运动的叠加效应，促使了该区地裂缝的 形成与发展. 临汾市地裂缝比较集中在5个区域发育，分别见于高堆、贾得、周家庄、梁村和临汾师大. 地裂缝灾害已经造 成市中心建筑物楼房开裂破坏，数十户民房因房屋遭破坏而被迫搬迁，大面积农田灌溉用水漏失，多处公路被切 断裂开，铁路也存在着潜在的危机，严重影响着临汾市的经济发展. 1 地下水的开采状况及时空特征 根据含水层的埋藏条件及地下水补给、排泄、径流特征，临汾市含水层可划分为3层 表1 . 临汾市自7 0 年代中期，地下水开采强度逐年增加，从7 0 年代后期，开始处于超采状态. 由于持续超采，使得地 下水位持续大幅度下降，至1992 年，年开采量已达51. 45h m3. 表1 临汾市含水层划分 T a b l e 1 A q u i f e r c l a s s i f i c a t i o n i n Li n f e n Ci t y 含水层地层时代岩性特征埋深富水性 上部 Qp 3, Qh 上为粘性土及粉细砂混层，下为砂 砾石层 10 ～2 5 m上部较差，下部较好 中部 Qp 2 粗砂层 西部小于8 0 m ，东部50 ～8 0 m ， 中部50 ～8 7 m 西部与城区富水性较好，其它地区 则较差 f i l e / / / E| / q k / z g k y d x x b / z g k y 99/ z g k y 990 1/ 990 12 2 . h t m （第 1／5 页）2 0 10 -3-2 3 15 56 55 中国矿业大学学报990 12 2 下部 Qp 1 细砂、粗砂 混层 15～17 m较差 临汾市有4个地下水集中开采区坟上水源地、翟村水源地、城区及城北（屯里、梁村一带）水源地，主要 用于农灌、工业用水及生活用水. 在时间上，中层水的强化开采主要起自197 6 年（中深水井主要集中成井时期）. 在此之前，地下水补给、排泄 基本处于均衡状态. 197 8 年已在城区范围形成降落漏斗，中心水位下降10 m 左右. 在此之后，由于在城区大量厂矿企 事业单位自备中深水井的成井及各水源地中深水井的增加，使得开采强度进一步加大，地下水位持续大幅度下 降. 据山西地矿局提供的资料，至198 6 年，地下水降落漏斗分布区域较197 8 年大幅度扩展，形成了一个NESW 向 展布、波及面积超过50 k m2的椭圆形降落漏斗，其中心位于下康、屯里一带，中心水位较197 8 年下降了30 m ，年 降幅接近4m . 另据本次勘查所进行的中深水井水位统测结果，至1995年1月，降落漏斗中心区域的地下水水位降幅比197 8 ～198 6 年期间趋缓，年平均降幅约在3m 左右，但外围降幅加大与198 6 年相比，降落漏斗向东西两侧大幅度扩展， 并在城区南部形成了一个新的降落中心 图1 . 图1 临汾市1995年地下水位降落漏斗分布 Fi g . 1 T h e d i s t r i b u t i o n o f w a t e r c o n e i n Li n f e n Ci t y i n 1995 2 地下水超采与地面沉降 地下水开采导致地面沉降的主要原因是随地下水位的下降，孔隙水压力逐步消散，有效应力逐步增加，土层 逐渐被压密所致. 根据太沙基（T e r z a g h l ）有效应力原理 σ＝σ′＋μ ， 式中σ为总应力；σ′为有效应力；μ为孔隙水压力. 对于所研究的区域，不同深度土层所承受的荷载主要由上覆土层的自重及地表建筑物荷重构成，其中上覆土 层自重为主体部分. 在土层固结压缩过程中，不同土层的压缩程度和固结特点具有很大差别. 一般来讲，含水层 （砂性土）在地下水位下降的过程中，孔隙水压力消散的速度较快，在较短的时间内即可完成其固结过程. 而粘 性土层中孔隙水压力消散极其缓慢，因此，其固结压缩将持续很长时间才能够趋于稳定. 2 . 1 压缩土层分布及沉降计算 在研究区范围，上部含水层多年水位观测结果表明，尽管水位总体变化趋势呈下降状态，但变化幅度小，而 深部承压含水层开采只限于市区及邻近范围，水位变化所波及的范围有限，所以，导致土层沉降的主要部分是中 部承压含水层的水位降低. 其水位降低幅度大，波及范围广，基本涉及整个研究区，所以，工作的重点是分析中 部含水层水位下降所导致的地面沉降变形问题. 基于中部承压含水层的埋藏条件及动态变化特征，分析认为其主要影响的深度大体限于10 0 m 埋深范围. 在这 个范围，含水层组主要是细砂－砂砾层，而压缩性土层则以粉质粘土为主，粘土、粉土层分布厚度相对较小，且 不同位置差异较大. 鉴于本次工作的基本目的及计算条件，选用南京水科院开发的“软土地基弹塑性固结变形计算程序 f i l e / / / E| / q k / z g k y d x x b / z g k y 99/ z g k y 990 1/ 990 12 2 . h t m （第 2 ／5 页）2 0 10 -3-2 3 15 56 55 中国矿业大学学报990 12 2 CO ND EP”进行土层沉降计算. 这种计算方法的主要特点是能够根据实际水位降，模拟分析孔隙水压的变化，并 获得孔隙水压力变化过程中土层中的应力、应变分布及变化趋势. 此模型已应用于北伦港码头、天津新港油罐地 基、宝钢引水工程大堤坝基等多项大型工程的地基沉降问题的分析，其计算结果均与实测结果比较接近. 2 . 2 计算模型的计算条件及参数的选取 沉降分析建立6 条计算剖面，其主要依据为山西地矿局2 13地质队提供的水位地质勘探成果，其中有3条纵向 剖面和3条横向剖面. 纵向剖面沿地下水位降落漏斗长轴方向，其中，汾河东、西岸及中部顺河各取一条，西岸剖 面（金殿伍默泊庄吴村）主要穿过二级阶地，东岸剖面（西秦鹅舍靳家庄阎马庄）主要穿过三级阶 地，中间剖面（尧庙郭家庄东芦韩村）主要穿过汾河一级阶地. 横向剖面垂直于地下水位降落漏斗的长轴 方向，在研究区北、中、南各有一条（分别为南麻南羊韩村孙曲李仵、良种场九洲堡郭家庄刘 村、鹅舍小贾尧庙金殿苏村），均横穿研究区内的所有地貌单元. 2 ．2 ．1 计算模型 1 模型主要分析中层承压水超采引起的土层沉降问题，将计算深度（埋深2 0 0 m ）范围的土层分为压缩层 （粘性土层）和含水层（砂性土）2 个层组，并概化为含水层与压缩层相间排列，以此为基础进行剖分，模型如 图2 所示. 图2 沉降计算模型 Fi g . 2 Co m p u t a t i o n m o d e l o f s u b s i d e n c e 2 模型采用平面应变模型，边界条件处理为地面为自由位移边界、自由排水边界，底部为固定位移边界、 不排水边界，两侧垂直方向为自由位移边界，水平方向为固定位移边界. 模型的横向尺寸大于降落漏斗的范围， 纵向尺寸取至主要抽水含水层以下. 3 计算以197 3年为起始时间, 即假定197 3年的水位为初始水位, 在此基础上根据197 8 ，198 6 ，1995及2 0 0 5年的水 位降分别计算相应时间的地面沉降值. 计算所依据的地下水水位值分别为197 8 ，198 6 年的数据依据山西地矿局提 供的资料，1995年的数据依据本次勘查地下水水位统测绘制的研究区地下水位等值线确定，2 0 0 5年的数据则是根 据198 6 ～1995年期间的水位平均年降幅以每年1 m 的降深推算所得. 4 计算采用4节点等参单元，沉降盆地边缘单元划分尺寸较小，盆地中间和降落漏斗外缘划分适当放宽. 计算 假定0 ～10 0 m 和10 0 ～2 0 0 m 两个深度范围的含水层段具有相同的承压水位，并在计算中对各自的降深作为独立的 因素分别考虑. 5 材料模式采用剑桥弹塑性模型. 计算参数有水下浮重度、静止侧压系数、渗透系数、剪切模量、压缩系 数、回弹指数等，参数的选取主要根据两个勘探孔的试验数据及前人工程地质勘察和抽水试验成果. 2 ．2 ．2 计算参数的选取 计算参数在综合分析本次勘查的试验成果及所收集的区内有关土层试验资料基础上确定，主要计算参数如表 2 所列. 表2 沉降计算模型参数 T a b l e 2 T h e p a r a m e t e r s o f c a l c u l a t i n g s u r f a c e s u b s i d e n c e 层组 浮重度/ k N.m-3 静止侧压力 系数K 0 剪切模量 G 0/ Pa 渗透系数 K y / 10 -7 K y ∶K x压缩系数Cc回弹指数Ce试验参数M C1 9. 90 . 300 . 462 4. 02 . 00 . 1590 . 0 150 . 145 S1 9. 10 . 2 00 . 6 02 49 2 0 01. 00 . 18 30 . 0 190 . 2 7 7 C2 10 . 20 . 300 . 462 3. 02 . 00 . 17 40 . 0 170 . 57 3 S2 8 . 90 . 2 00 . 522 49 2 0 02 . 00 . 18 30 . 0 130 . 57 3 f i l e / / / E| / q k / z g k y d x x b / z g k y 99/ z g k y 990 1/ 990 12 2 . h t m （第 3／5 页）2 0 10 -3-2 3 15 56 55 中国矿业大学学报990 12 2 C3 9. 90 . 330 . 462 0 . 02 . 00 . 0 8 80 . 0 2 20 . 397 S3 9. 00 . 2 30 . 522 49 2 0 01. 00 . 18 30 . 0 130 . 38 0 C4 10 . 50 . 330 . 4618 . 02 . 00 . 17 10 . 0 330 . 30 0 S4 10 . 00 . 2 30 . 6 02 49 2 0 01. 00 . 18 30 . 0 130 . 47 7 C5 10 . 70 . 350 . 332 4. 02 . 00 . 12 30 . 0 2 40 . 350 S5 9. 00 . 2 50 . 6 02 49 2 0 01. 00 . 18 30 . 0 130 . 397 C6 10 . 30 . 350 . 4614. 02 . 00 . 1590 . 0 150 . 397 2 . 3 计算结果分析 通过计算模型的运算，得出该区的沉降结果［1～3］. 其特征如下 1 不同阶段地面沉降区分布与同期水位降深漏斗分布区在形状上相似，沉降中心位于沿汾河的漫滩、一级阶 地范围，但明显向西偏移. 2 与地下水位降深漏斗的扩展相比，沉降盆地的扩展极不明显，主要表现为盆地边缘沉降梯度的变化. 另 外，对比197 8 ，198 6 ，1995年的地面沉降情况，197 8 年地面沉降极不明显，沉降速率为6 m m / a ，在197 8 ～198 6 年期 间呈明显加速趋势，沉降盆地中心的沉降速度为30 m m / a ，而在198 6 ～1995年期间，地下水位下降速度明显趋于平 缓，而沉降速率仍保持10 m m / a 以上的较高水平，这表明地面沉降比地下水位下降明显滞后. 3 沉降中心两侧的垂直形变梯度与地貌单元有明显对比关系. 在沉降中心两侧二级阶地、一级阶地范围的形 变梯度与三级阶地、洪积扇范围的形变梯度呈明显差异. 汾河西岸二级阶地后缘与洪积扇、三级阶地交界的单元 界限方向存在一形变梯度陡变带，在东侧的形变梯度陡变带也基本沿单元界限方向展布. 4 沉降引起的水平位移方向指向盆地中央，水平位移最大值出现在沉降盆地外部边缘的拐点处，模型对应于 198 6 ，1995及2 0 0 5年3个计算时段沉降盆地的水平位移如图3所示. 图3 临汾沉降的水平位移 Fi g . 3 H o r i z o n t a l d i s p l a c e m e n t o f s u b s i d e n c e i n Li n f e n Ci t y 3 地面沉降与地裂缝灾害的时空关系 随着中层含水层的持续超采，地面沉降进一步加大. 当中层含水层降不均匀时，引起上部土层（包括浅部地 层）骨架颗粒不均匀沉降从而产生应力集中. 当应力集中达到一定程度时，则会产生地裂缝. 3. 1 空间关系特征 对比地裂缝发育位置与沉降区分布 以1995年为例，见图4所示 ，可以发现，鹅舍、周家庄、师大、梁村等处 地裂缝与高堆地裂缝分别分布于沉降盆地的东、西两侧，且在展布方向上二者也表现有明显的对应关系. f i l e / / / E| / q k / z g k y d x x b / z g k y 99/ z g k y 990 1/ 990 12 2 . h t m （第 4／5 页）2 0 10 -3-2 3 15 56 55 中国矿业大学学报990 12 2 图4 地裂缝发育位置与沉降形态 Fi g . 4 Ea r t h -f i s s u r e d e v e l o p m e n t s i t e s a n d s u b s i d e n c e 3. 2 时间关系特征 研究区范围地面沉降主要发生于197 8 年以后，从沉降形变特征看，高堆、周家庄、梁村及师大四处地裂缝均 处于197 8 年以后的沉降盆地边缘区，而鹅舍、龙祠地裂缝发育位置在198 6 年以前的地面沉降还很不明显，在198 6 ～1995年期间所产生的沉降也保持较低水平. 因此，除这两处地裂缝与地面沉降的关系不太明朗外，其他地裂缝与 沉降的关系还是比较明显的. 其中，周家庄、山西师大地裂缝在产生时间上与沉降表现有一致性，而高堆、梁村 地裂缝尽管在发生的时间上早于地面沉降期，但其明显活动时期与地面沉降在时间上的对应性也比较明显. 4 结 论 临汾市地裂缝的形成与发展与该地区地下水超采在时空上有明显的相关关系. 因此对地下水资源的开发应加 强管理，有效地限制临汾市地下水过量开采是当务之急. 作者简介 董东林，男，197 0 年生，助教，工学硕士 作者单位董东林 武 强 孙桂敏中国矿业大学资源开发工程系 北京 10 0 0 8 3 姜振泉 隋旺华 田宝霖中国矿业大学资源与环境科学学院 徐州 2 2 10 0 8 参考文献 1 姜振泉，王晓波，张 京. 山西断陷带地裂缝的成因及其发育条件. 中国矿业大学学报，1997 , 2 6 3 7 4～7 8 2 董东林，武 强，孙桂敏. 山西临汾地裂缝G I S预测的初步研究. 中国地质灾害与防治学报，1996 , 7 4 16 ～2 0 3 董东林. G I S支持下的山西临汾市地裂缝灾害敏感性评价与发展趋势预测 ［硕士学位论文］. 中国矿业大学北京 研究生部，1997 收稿日期 1998 -0 4-30 f i l e / / / E| / q k / z g k y d x x b / z g k y 99/ z g k y 990 1/ 990 12 2 . h t m （第 5／5 页）2 0 10 -3-2 3 15 56 55