采煤塌陷区水体淹没范围预测方法研究.pdf

声明声明下面论文由免费论文教育网 http//www.PaperE 用 户转载自互联网，版权归原作者所有，本文档仅供参考，严禁抄袭 免费免费论文论文教育教育网网 -1- 中国中国科技论文在线 采煤塌陷区水体淹没范围预测方法研究∗ 科技论文在线 采煤塌陷区水体淹没范围预测方法研究∗ 柳聪亮 1，2，谭志祥1，2，李培现1，2，白力改1，2，邓喀中1，2 1、中国矿业大学江苏省资源环境信息工程重点实验室，江苏 徐州（221008） 2、中国矿业大学环境与测绘学院，江苏 徐州（221008） E-mailliucongliang66 摘摘 要要为了准确预测开采后塌陷区地表的淹没范围，确保煤矿生产安全，本文用概率积分 法进行地表下沉预计，综合利用 AutoCAD 和 Surfer 绘图软件实现了沉陷前后矿区地形图的 变更，对地表水体水位进行分析，提出了采煤塌陷区淹没范围的预测方法。以山西某矿开采 沉陷后塌陷区淹没范围的计算为例， 说明了进行塌陷区地表淹没范围预测的过程， 提出了相 应的安全措施，为确保该矿生产安全和塌陷淹没区域治理提供了科学依据。 关键词关键词淹没范围；开采沉陷；塌陷；等高线；预测 中图分类号中图分类号TD3 1.引言引言 我国幅员辽阔，地表水系丰富。江河、湖泊等地表水体下压煤量大[1]，如何安全采出地 表水体下的煤炭资源和实现塌陷积水区环境的有效治理是矿山企业亟待解决的难题。 地下煤 层开采面积达到一定范围后，移动与变形将波及到地表，使地表产生大面积塌陷，导致塌陷 区内水系遭受破坏，势必造成沉陷区积水深度和淹没范围扩大[2]。准确预测开采沉陷后塌陷 区淹没范围，对预测煤矿可能渗水点，预防淹井事故，评价开采沉陷对土地和环境的影响、 制定后期治理方案； 合理安排和调整矿区建筑物与生产设备及煤矿安全生产有重要意义。 传 统的预计方法，仅计算开采后地表移动和变形大小[1]，没有根据地形情况对塌陷区内的地形 图进行更新和计算，不能实现塌陷区淹没范围的预测。 AutoCAD 是由美国 Autodesk 公司开发的交互式计算机辅助绘图和设计软件包，功能强 大，实用。它具有强大的图形编辑分析功能和开放的开发环境,可以实现图形和数据之间的 相互转换，图形坐标还原，方位旋转等，因此倍受矿山测量工作人员的青睐，已成为大多数 矿山企业矿图资料的管理平台[34]。Surfer 软件是美国 Golden Software 公司研制的基于 Windows 操作系统的二维、三维插值绘图软件。 该软件提供了 12 种内插方法，不等间距 的数据也可以用它作图。它还提供了多种图形、图像、文件格式输入输出接口，可与各大 GIS 软件实现文件和数据交流[5]。 山西某矿进行河下采煤，矿区为丘陵地貌，地形较为复杂，用常规的方法难以预测开采 沉陷后塌陷区淹没范围。本文在开采沉陷预计的基础上，综合利用了 AutoCAD 的数据、图 形相互转换功能、 地形图接边功能及图形分析功能和 Surfer 强大的等值线绘制功能对该矿塌 陷区淹没范围进行预测， 并提出了相应的安全措施， 为该矿生产管理和塌陷区治理提供了科 学依据。建立了一种采煤塌陷区地表淹没范围预测的新方法。 2 淹没范围预测方法淹没范围预测方法 从现有地形图中提取采动影响区域的原始地形数据；采用概率积分法预计相应点下沉 值；根据预计结果，实现地形数据更新和变更后地形图绘制；然后结合地表水体水文特征和 地形变更情况进行水位分析，进而实现淹没范围预测；沉陷预计等值线图、淹没范围和地形 ∗国家自然科学基金资助项目（40772191）；“十一五”国家科技支撑计划重点项目（2006BAC09B01） -2- 中国中国科技论文在线科技论文在线 图叠加分析可进一步对矿区积水的次生效应（灾害）进行预测和研究。具体过程见图 1。 2.1、沉陷预计数据准备、沉陷预计数据准备 （1）确定采动影响区煤层开采对地表的影响范围与采空区大小、煤层埋藏深度、边 界角及松散层移动角有关[1]。采空区边界至影响区边界的水平偏距 走 L 、 下 L 、 上 L 分别满足 图 1 淹没范围预计流程 Fig1 the area of mining inundation prediction flow 关系式（1）-（3） ϕδctghctgHL 0平均走 （1） ϕβctghctgHL 0下下 （2） ϕγctghctgHL 0上上 （3） 式中 走 L走向方向偏距，m； 下 L下山一侧偏距，m； 上 L上山一侧偏距， m；H上覆岩层厚度，m； 0 δ移动影响边界角,度；h松散层厚度，m；ϕ松 散层移动角，（）。 为了降低地形图形接边误差对淹没范围预测的影响， 应适当扩大地形图更新范围。 可取 走 L 、 下 L 、 上 L 中的较大者再加 10m为最终的水平偏距 L，将采空区边界线向外平移 Lm， 矿图 预计系统 预计文件 下沉预计文件 预计点三维坐标文件 更新高程数据文件 Surfer AutoCAD 等值线图 开采沉陷后预计地形图 淹 没 面 积预计 土 地 损 失预算 道路桥梁 淹没及治 理分析 房屋淹没及 损害分析 各种治理 工程土方 量预算 淹 没 次 生 灾 害 分 析 （滑坡等） 采动影响区地形图 -3- 中国中国科技论文在线科技论文在线 绘出地表采动影响区边界，然后提取影响区内的地形信息。 （2）生成预计文件用 GIS 软件提供的“查询”功能得到高程注记点的三维坐标数据， 并按顺序编号制成预计点三维坐标文件。删除该文件的高程信息，平面坐标和点号不变，然 后添加相应的地质采矿条件及预计参数数据生成预计文件。 2.2、开采沉陷预计、开采沉陷预计 常用的开采沉陷预计方法有影响函数法、剖面函数法、数值分析法和力学分析法等。其 中基于随机介质理论的概率积分法在我国应用最为广泛。 随机介质理论认为岩层移动是一个 服从统计规律的随机过程， 从统计观点出发， 可以把整个开采区域分解为无限个微小单元的 开采，整个开采对岩层及地表的影响等于各单元开采对岩层及地表影响之和[1]。按随机介质 理论，单元开采引起的地表单元下沉盆地呈正态分布，且与概率密度的分布一致。因此，整 个开采引起的任意点下沉方程可以表示为二元概率密度函数的二重积分 dtdse r WyxW Z r tysx Q ∫∫ −− − 0 2 0 2 22 0 1 , π （4） ,yxW 任意点 ,yx 的下沉值； 0 W 最大下沉值；Q开采区域倾向方向 边长；Z开采区域走向方向边长。 式（4）是常规概率积分模型的数学表达式，在工程实践中，国内外学者对其不断研究 改进，改进的模型有扩展概率积分模型[6]、非充分采动概率积分模型[7]、倾角变化的概率 积分模型[8]和概率积分分层预计模型[9]等。经过不断的改进，概率积分法已能满足复杂地质 采矿条件下的地表下沉预计。并且该方法便于编程实现。 根据矿区具体的地质采矿条件选择合适的概率积分模型和预计参数对预计点进行下沉 预计，得到下沉预计文件。下沉预计文件包含预计点的平面坐标和下沉预计值，其点号顺序 与预计点三维坐标文件相同。 2.3、采动影响区地形图更新、采动影响区地形图更新 由于预计点三维坐标文件与下沉预计文件点号顺序完全相同， 将两个文件中相同点的高 程减去其下沉值得到该点的预计高程值。用预计高程取代预计点三维坐标文件中的原始高 程，点号和平面坐标不变，由此生成一个更新高程数据文件。用更新后的地形数据重新生成 等高线，将等高线拷贝到原地形图中进行遇地物断开处理。为了使直线恢复共线状态，对其 进行直线平差处理。对折线、曲线，采用加权平均值的方法进行修正，实现对接等高线平滑 处理，完成矿区地形图变更。 2.4、淹没范围预测及相关分析、淹没范围预测及相关分析 矿区地表水系类型不同，开采沉陷后水位变化的机理也不同[1012]。须结合矿区地形现 状、 地表水体现状以及预计的塌陷区地形特征进行水位变化分析。 从而确定开采沉陷后塌陷 区积水水位。 然后通过内插计算， 绘制出预测的水位线。 水位线包围的范围即积水淹没范围。 对淹没范围进行统计分析。可计算淹没面积，进行淹没区形状特征分析。通过沉陷预计 等值线图、淹没范围和地形图叠加分析可以预测将被淹的建筑物、生产设备、矿区居民地、 水利设施和电力设施等；预测矿区可能的渗水点；可预测房屋、道路、桥梁和堤坝等建筑物 的损害等级；山体滑坡等地质灾害发生的可能性，为生产管理提供可靠的信息。 -4- 中国中国科技论文在线科技论文在线 3 工程实例工程实例 山西某矿进行河下采煤，该矿 5602 和 5603 工作面开采二叠纪山西组 3煤层，走向东 西向。具体地质采矿条件及预计参数见表（1）。河水自西向东流经 5603 工作面的正上方， 开采前河岸标高由 874m 逐渐降至 870m，西部落差较大，东部地势较为平缓，河床最低标 高 865m，最大流量 6～7 m3/s，最小了流量 0.3 m3/s，为常年有水河流。该河南邻“一高地” 地势陡而高，最高点标高 915.3m河岸以北地势缓而低，最高标高 882m，属于山区地形。 具体情况见图 2。 5602、5603 工作面开采后，其上方将形成塌陷区，地形发生变化，河流不会被疏干， 河床下沉会导致河面加宽，河的积水面积和深度都将明显扩大。以 AutoCAD 和 Surfer 为平 台用上述方法对塌陷区进行淹没范围预测。 表 1 地质采矿条件及预计参数 Tab.1 geological and mining digital and prediction parameters 工作面 岩层厚，m 土层厚，m 采深，m煤层厚， m走向长，m倾向长，m 煤层倾角（） 5602 295 10 305 6 1800 240 6～8 5603 290 10 300 6 1800 240 6～8 工作面 下沉系数 q 主要影响角 正切 tgβ 水平移动 系数 b 最大下沉 角（） 拐点偏移距综合边界 角（） 5602 0.86 2.5 0.31 84.4 0.05～0.1H 59.8 5603 0.86 2.5 0.31 84.4 0.05～0.1H 59.8 注表中岩层厚、土层厚、采深、煤层厚均为平均值。 图 2 采动影响区域采前三维地形 Fig2 the three-dimensional relief map before mining 1生成预计点文件由于该矿区表土层厚度、煤层倾角较小，用平均深度和综合移动 边界角计算采动影响范围即可。采动影响区域现状地形图见图 3。 -5- 中国中国科技论文在线科技论文在线 图 3 采动影响区域现状地形 Fig3 the relief map before mining （2）下沉预计用表（1）中的地质采矿条件及预计参数数据和预计点文件制成预计文 件，用概率积分法对采动影响区域各地形点进行下沉预计。预计效果图见图 4。 410800 410800 410600 410600 410400 410400 410200 410200 410000 410000 409800 409800 409600 409600 409400 409400 409200 409200 409000 409000 408800 408800 408600 408600 4041400 4041200 4041000 4040800 40406004040600 40404004040400 40402004040200 4040000 4040000 4041400 4041200 4041000 4040800 Ⅸ Ⅷ Ⅶ Ⅵ Ⅴ Ⅳ 4 75 . 13 0 3 20 479. 2403 18 4 82 .4 5 0 31 6 4 8 6 .41 0 312 48 8. 51 0310485.04 0308 50 2.0 0 11 3 0 1 115 03 .1 0 105 5 0 3. 2 50 3 . 7 01 03 Ⅲ 预制厂 沥青厂 预制厂 Ⅱ Ⅰ Ⅺ 20 0 7 Ⅻ 0. 80 0 .2 5 1 .78 0.0 3 3. 10 大庄渔塘 大黄庄渔塘保安煤柱边界（未批） 410444,40408108 88 8 07 56 0 2 工作 面 56 水仓 560 3 回风 巷 506.60101 56 -2 495.6L031 505. 4L032 518. 5L 034 515.5L 035 475.3L038 484.7L037 491.3L036 502.5J30506.7J 28 512.6J 27 513.9J26 507.8J25 496.4J24 489. 7J23 485.6 J22 485. 1 J20 487.2 J21 4 82 .9 C01 483. 2 C0 2 48 78 -F6 H 1m ∠80 5.30 504.47 872.57 南-59 巷 车 猴 4 86 .7 0 30 6 0 3 02 5.80 588.75 883.13 南-3 5 .70 552.62 893.28 南-12 4.90 528.33 877.20 南-1 7 7.60 524.45 896.60 1 -4 6.6 0 578.48 908 .84 南-7 5.10 576.98 874.88 南-8 5.65 558.18 875.28 1 -2 3.25 52 0.63 880. 13 五-22 2 6.80 577.29 872.34 五-2 20 6.20 605.12 869.15 南-4 6.42 619.50 877.70 9 1 6.04 612.27 885.31 2 7 5.50 502. 80 874.84 南-1 8 5.80 南-31 5 1 0 55 0 5 60 570 5 80 5 90 600 4 90 500 510 5 20 530 540 550 560 5 7 0 580 480 490 500 5 1 0 520 490 500 520 4 60 470 480 4 90 50 0 5 10 520 530 5 40 5 50 5 60 5 70 580 590 60 0 61 0 H3∠6 5 S 16 断层崔家庄3 X5 560 2运输巷 5602 回风巷 56- 3 所 电 变 巷 带 皮 巷 道 轨 高速公路连接路 H3 ∠ 65 S 17 06 马庄 马村 提水站 断层防隔水 煤 柱 界 边 柱 煤 安 保 庄 马 保安煤柱边 界 保安煤柱边界 傅 村保安 煤 柱边 界未批 断层防隔水煤柱 潞矿地字 [2 0 0 5 ] 第 17 号 崔家庄 3 断 层防隔水煤 柱 潞矿地字 [20 0 5]第17号 崔家 庄3 4 9 7 .3 5 0 2. 8 5 1 0.0 515 .7 521 .6 5 2 7 . 3 5 3 1.9 5 4 4 .6 5 5 2 .9 55 7 .8 5 6 6 .7 5 7 0 .8 5 7 5.2 5 7 7. 9 58 2 .2 584 .1 5 8 5.4 5 9 9 .7 6 0 6.3 6 1 3 . 4 613 .0 604.5 582 .6 574. 8 5 6 8 . 5 5 5 7 .8 5 5 3 .0 5 4 8 . 7 5 4 7 .1 5 4 6 . 0 5 4 5 .3 5 4 5.8 543. 6 5 3 4.3 532. 1 528. 05 2 5 . 95 2 1.6 5 1 2 .9 5 0 6.3 4 9 0 .0 4 8 4 . 3 4 8 4 .2 485.2 482.9 484.6 490.1 488.7 490 .0 494.2 490.1 496.8 499.9501.2 508.0 513. 8 513.2 514.6 515 .7 514.8 517.2 509.8 512.4 513.3 504.9 507.4 495.0 493. 8 490.5 489.5 488.8 480.4 479.5 477.1 1 2 0 0 1 2 0 2 1 20 4 1 2 0 6 1 2 08 1 21 0 1 2 1 2 1 2 1 4 1 2 1 6 1 2 1 8 1220 1222 1224 1 2 2 6 1 22 8 12 3 0 1 2 3 2 1 23 6 1 23 8 1 2 4 0 1245 1243 1 2 4 1 1 2 3 9 1237 1 2 3 5 1 2 3 3 1 2 3 1 12 2 9 1 22 7 12 2 5 1 2 2 3 12 2 1 1 2 1 9 1217 1 2 1 5 1 2 13 1211 1209 12 0 7 12 0 5 1 20 3 1 20 1 B 18 A17 B1 7 A16 A15 G6 G4 B16 A14 B15B 14 B 13 B3 B2 B1 B4 新B 12 A12 A1 1 新B11 L033 A10 B11 A9 B10B09 L030 A07 L 029 B0 8 饭店 提水站 公 路 公 漳 河 漳 河 黄庄水渠 十字 道背 斜 十字道背斜 H2 m 7 802 -5 H1 .6m7 80 2-4 H 1.5 m 7 80 2-3 H 1.5 m 7 802 -2 H 1 .3 m 7 8 02 -1 H 3 . 5m∠ 4 0 H3 .5 m∠65断层56 -2 断层5 6-1 H 5 m∠ 70 H6 m∠ 6 0 78-F5 H2.2m∠40 巷 风 回 用 专 图 4 沉陷预计结果 Fig3 mining subsidence prediction result 前河床在采动影响区东侧标高变化不明显，所以取下游未受采动河段平均水位 868.5m 为采动后采空区上方的水位来预测沉陷后该河段的河水淹没范围。 预计淹没范围见图 5 蓝线 包围区域。 （5）相关分析根据上述对 5602、5603 两工作面开采地表下沉和河水淹没范围预（3） 矿区地形图变更 用预计得到的下沉预计数据文件和预计点三维坐标文件生成更新高程数据 -6- 中国中国科技论文在线科技论文在线 文件，在 Surfer 中绘制采动影响区域的变更地形图输出 AutoCAD 能辨认的（*.dxf）格式的 等高线图，然后在 AutoCAD 中实现地形图更新，更新效果图见图 5。 （4）淹没范围预测河上游的水库已被疏干，且河与水库连接处已垫高，地表塌陷后 河水不会流进水库中。 地面塌陷后河流上游未受采动影响的河段其水位不变； 受采动影响的 河段水位随地表的沉陷而略有下降但不会低于塌陷区下游未受采动影响河段水位； 又因为未 开采测结果分析，积水河段东西长 1100m，南北最宽处 430m，最窄处 130m，淹没面积约 249190m2。高速公路连接路、公路部分路段及过河拱桥等交通设施将会被淹没；饭店及部 分民房等建筑物将被淹；电力设施将有部分线路处于塌陷盆地积水区；河岸南边“高地”在汛 期有滑坡的可能。另外，上述淹没范围是根据最近三年水位的平均标高预测的，实际上，水 位随补给量和蒸发量变化而变化，在汛期淹没面积将会扩大，为了保证矿区交通、电力、水 利设施及民用建筑的安全，应做好各项防汛工作。 （6）安全措施及建议在上述研究的基础上，本文提出了塌陷区淹没后防治技术措施 a、采煤应避开汛期进行；b、连续开采避免形成永久性边界对河床的损害；c、及时封堵地 表缝隙减少渗流量；d、应加强对河南岸“山地”的变形监测，避免滑坡灾害发生；e、对电力 设施进行加固、道路改道、淹没区域的房屋要做好搬迁工作，淹没区外受采动影响的建筑物 进行加固；f、地表稳沉后对淹没区河道进行“挖深垫浅”治理，有效利用土地。 开采影响范围 淹没范围边界 淹没范围边界 410800 410800 410600 410600 410400 410400 410200 410200 410000 410000 409800 409800 409600 409600 409400 409400 409200 409200 409000 409000 408800 408800 408600 408600 4041400 4041200 4041000 4040800 40406004040600 40404004040400 40402004040200 4040000 4040000 4041400 4041200 4041000 4040800 Ⅸ Ⅷ Ⅶ Ⅵ Ⅴ Ⅳ 475. 13 0 320 479. 2 4 0 318 4 82.45 031 6 486.41 0312 48 8.5103 10 4 85 . 04030 8 5 02 .001 13 01 11503 .1 01 05 5 03.2 50 3.7 01 03 Ⅲ 预制厂 沥青厂 预制厂 ⅡⅠ Ⅺ 20 07 Ⅻ 0.80 0.25 1.78 0.033.10 大庄渔塘 大黄庄渔塘保安煤柱边界 （未批） 41 044 4,40 408 10 8 88 8 0 7 56 02 工 作面 56 水仓 5 603 回风巷 50 6.601 01 56-2 49 5.6L031505 .4L0 32 518. 5L034 515. 5L035 47 5.3L03848 4.7L0 37 491.3L0 36502 .5J30 5 06.7J2851 2.6 J27 51 3.9J26 5 07.8J25 4 96 . 4J2 4 489 . 7J23 4 85 . 6 J 22 485.1 J20 487.2J21 4 8 2 .9 C01 4 83 .2 C0 2 48 78-F 6 H1m ∠ 80 5.30 504.47 872.57 南-59 巷 车 猴 486 . 7 0 306 030 2 5.80 588.75 883.13 南-3 5.70 552.62 893.28 南-1 2 4.90 528.33 877.2 0 南-17 7.60 524.4 5 896.60 1 -4 6 . 60 578.48 908.84 南-7 5.10 576.98 874.88 南-8 5.65 558.18 875.28 1 -2 3.25 520.63 880.13 五-2 22 6.8 0 577.29 872.3 4 五-22 0 6.20 6 0 5.12 869.15 南-4 6.42 619.50 877.7 0 91 6.04 612.27 885.31 2 7 5.50 502.80 874.84 南-1 8 5 . 80 南-31 51 0 550 560 570 5 80 59 0 6 00 49 0 500 510 5 2 0 530 540 55 0 560 57 0 5 80 4 80 490 5 00 5 10 520 490 5 0 0 520 460 4 70 480 4 9 0 500 510 520 5 30 540 550 560 570 580 590 6 00 610 H3∠65 S16 断层 崔家庄 3 X 5 560 2运输巷 5602 回风巷 56 -3 所 电 变 巷 带 皮 巷 道 轨 高速公路连接路 H 3∠65 S1 7 0 6 马庄 马村 提水站 断层防 隔 水 煤柱 界 边 柱 煤 安 保 庄 马 保 安煤 柱边界 保安煤柱边界 傅村 保 安 煤柱边界 未 批 断层防隔水煤柱 潞矿地字 [2005] 第 17 号 崔家庄 3 断层 防 隔水煤 柱 潞矿地字 [2005] 第17号 崔家 庄3 49 7.3 5 02 . 8 510 .0 5 15.7 5 21 .6 52 7.3 53 1 .9 54 4. 6 55 2. 9 55 7 .8 5 66.7 570 .8 57 5 .2 577 .9 582 .2 584 .1 5 85 . 4 599 .7 60 6 .3 6 13.4 613. 0 6 04.5 582 .6 574 .8 568 .5 55 7. 8 553.0 548. 7 547 .1 546 .0 5 45 . 3 5 45 . 8 54 3. 6 534 .3 532 .1 52 8.0 52 5.9 5 21 . 6 512. 9 50 6.3 4 90 . 0 48 4 .3 484. 2 48 5.2 4 82.9 4 84 . 6 49 0.1 4 88.7 490. 0 4 94.2 49 0.1 496.8 4 99 . 9501 .2 508.0 513 .8 513.2 514 . 6 5 15 . 7 514.8 517.2 509 . 8 5 12 . 4 5 13.3 504.9 5 07.4 49 5.0 493 . 8 49 0.5 4 89.5 4 88 . 8 4 80.4 47 9.5 477.1 1200 1202 1204 120612 08 1210 12 12 12 14 12 16 1 21 8 1 22 0 12 22 12 24 1226 1 22 8 1230 1 23 2 1 23 6 1238 12 40 1245 1243 1241 1239 1237 1235 1233 1231 1229 1227 1225 1223 1221 1219 1217 1215 1213 1211 1209 1207 1205 1203 1201 B18 A17 B 17 A16A15 G6 G4 B16 A14 B 15B14 B13 B3 B2 B1 B4 新B1 2 A 12 A11 新B 11 L 033 A10 B 11 A9 B 10 B09 L030 A07 L029 B08 饭店 提水站 公 路公 漳 河 漳 河 黄庄水 渠 十 字 道 背 斜 十字道背 斜 H 2m 780 2-5 H 1.6m7802-4 H 1.5m7802-3 H1. 5m 7802-2 H1 .3m 7802 -1 H3 .5 m∠40 H3. 5m∠65断层 56-2 断层 56-1 H5m ∠ 70 H6m∠60 7 8-F5 H2.2m∠40 巷 风 回 用 专 图 5 沉陷后淹没范围 Fig3 the inundation after mining subsidence 4 结束语结束语 在地表沉陷预计和水位分析的基础上，AutoCAD 和 Surfer 相结合进行塌陷区地表淹没 范围预测，综合利用了 Surfer 强大的等值线绘制功能和 AutoCAD 的图形编辑分析功能，该 方法简单实用高效，易于掌握，是一种可以推广的预计方法。本文用该方法预测了山西某矿 开采后淹没范围，给采煤塌陷区积水区域的治理和合理利用提供了科学依据。 -7- 中国中国科技论文在线 参考文献 科技论文在线 参考文献 [1]何国清，杨伦，凌赓娣等． 矿山开采沉陷学[M]． 徐州中国矿业出版社．1991 [2]徐良骥，严加平.煤矿塌陷区地表水系综合治理[J].煤炭学报，2007，5（32）317～320 [3]毛胜明， 胡丙中．谈 AutoCAD 技术在地形图修测中的应用[J]．四川测绘，2002（2）22～25 [4]马文祥， 金芳丽．在 AutoCAD 中实现高程信息的快速提取[J]．测绘通报，2006（3）54～55 [5]马天勤．Surfer 软件在开采沉陷可视化中的应用[J]．矿业工程，2008（6）48～50 [6]李永树.任意分布形式煤层开采地表移动预计方法[J].煤炭学报，1995，20（6）619～624 [7]戴华阳，王金庄.非充分开采地表移动预计模型[J].煤炭学报，2003，28（6）584～587 [8]戴华阳.基于倾角变化的开采沉陷模型及其 GIS 可视化应用研究[J].岩石力学与工程学报，2002，21（1） 148～149 [9]汤伏全,西北后黄土层矿区地表移动预计方法[J].西安科技大学学报，2005，25（3）469～471 [10]韩丛丛，吴泉源，李轶平．基于 DEM 的龙口煤矿塌陷区水域淹没分析[J]．山东师范大学学报（自然科 学版），2008（23）63～65 [11]叶芳毅，张力．利用大比例尺地形图计算库区土地淹没面积．人民长江，2007（38）39～40 [12]刘小生 黄玉生 基于 Arc/Info 的洪水淹没面积的计算方法 测绘通报 2003（6）46～48 Study on Prediction to Inundation Area of Mining Subsidence’s Surface Liu Congliang1, 2, Tan Zhixiang1, 2 , Li Peixian1, 2 ,Bai Ligai1, 2 ,Deng Kazhong1, 2 1 China University of Mining Mining Subsidence;Subsidence ;Contour ; Prediction 作者简介作者简介柳聪亮1985 ，男，河北省保定人，硕士研究生，主要研究方向从事矿山 开采沉陷与矿山空间信息及其数字化方面的研究工作。