081-基于Matlab的开采沉陷预计及可视化研究

开采沉陷知识总结

开采沉陷知识总结 地表到充分采动时的采空区面积。半无限开采：工作面煤壁一侧的煤层未被采动，而另一侧的煤层全部采空的开采情况。填空（包含解释） 1、岩层移动的形式：弯曲、冒落、片帮、底板隆起、滚动、岩石沿层面滑移（1）弯曲岩体部分下落，向下移动的现象。特点：不在保持其原有的层状性，岩体积增大。（3）片帮当底板角软时，在矿体采出后，底板在垂直方向减压，水平方向受压，导致底板向采空区方向隆起的现象。（5）滚动在开采倾斜矿体时，岩石在自重力的作用下，除产生沿层面法线方向的弯曲外，还会产生沿层面方向的移动。特点：岩层沿层面滑动，层间有错动。 2、岩层三个移动区特征：充分采动区（减压区），最大弯曲区，岩石压缩区。（1）充分采动区移动特征是：1）下部岩体破碎成块状，上部岩体断裂、离层、裂隙；2）竖向方向受拉、横向方向受压，3）层内各点的移动向量与煤层法线方向一致，在同一层内的移动向量彼此相等；4）移动量最大。（2）岩石压缩区（支承压力区）特征：沿层面方向受拉、层面法线方向受压。（3）最大弯曲区（Ⅱ、Ⅱ’）特征：1）在此范围内岩层向下弯曲的程度（曲率）最大。2）在层内产生沿层面方向的拉伸变形和压缩变形。

3、移动稳定后岩层内的三带：冒落带，裂缝带，弯曲带。冒落带（垮落带）垮落带是指由采矿引起的上覆岩层破裂并向采空区垮落的岩层范围。、垮落带内岩层破坏特点是：（1）分带性（2）碎胀性，碎胀系数值恒大于1，一般在 1、5～ 1、80之间。（3）可压缩性（4）垮落带高度主要取决于采出厚度、上覆岩层的碎胀系数、岩性等，通常为采出厚度的3～5倍。薄煤层开采时垮落高度较小，一般为采出厚度的 1、7倍左右。顶板岩石坚硬时，垮落带高度为采出厚度的5～6倍；顶板软弱时，垮落带高度为采出厚度的2～4倍。实践中可用下式近似估算垮落带高度：h=m/(k-1)cosα式中 h采出矿层厚度；K矿层倾角。（5）形态，随倾角不同而变化裂缝带（断裂带）在采空区上覆岩层中产生裂缝、离层及断裂，但仍保持层状结构的那部分岩层称为裂缝带。特征：（1）岩体内裂隙多，但仍保持层状结构（2）裂隙带高度与岩性密切相关：硬，高，软，低冒落带+裂缝带=冒落裂缝带（导水裂隙带）软弱岩石：9～12m （m为开采厚度）；中硬岩石：12～18m；坚硬岩石：18～28m。弯曲带弯曲带指的是裂缝带之上直至地表的整个岩系。特点：1）保持整体性和层状结构，2）移动过程连续而有规律，不存在或极少存在离层裂缝。3）在竖直面内，各部分的移动值相差很小。4）H 弯＞＞H裂（当H大时）三者的关系：从采空区往上：冒落带→裂

开采沉陷

名词解释： 1.半无限开采：沿着工作面推进方向x 区间0到正无穷上被开采，而沿垂直工作面推进方向的开采尺寸足够大，使之达到充分采动。（1.2--1.4H ） 2.主要影响半径：半无限开采主要的地表移动和变形均发生在r x -=～r +的范围之内，称r 为主要影响半径。 2.下沉盆地的角度参数： 边界角：开采达到或接近充分采动是，移动盆地主断面上盆地边界点和采空区边界点连续与采空区外侧水平线的夹角 移动角：开采达到或接近充分采动时，移动盆地主断面上临界变形点和采空区边界点连线与采空区外侧水平线的夹角。下山、上山和走向方向的移动角分别用β、γ和δ来表示；i=± 3.0mm/m ；E=±2.0mm/m ；K=±0.2×10-3/m 。 裂缝角：开采达到或接近充分采动时，采空区上方地表最外侧位置裂缝和采空区边界点连线与采空区外侧水平线的夹角 最大下沉角：在移动盆地的倾斜主断面上，采空区的中点与地表下沉盆地中点的连线与矿层下山方向水平线的夹角 3.启动距：地表开始移动式工作面推进距离 地表开始移动：观测地表下沉值达到10mm 地表移动时间：从地表开始移动到地表停止移动的持续时间。分为启动。活跃。衰减阶段，1.67mm/d ，百分之85 4.减沉开采：是通过改变采场顶板管理方法控制顶板下沉量，达到减缓地表沉陷量 5.协调开采：根据开采引起地表移动变形分布规律，通过合理的开采布局，开采顺序，方向时间等方法，减缓开采地表变形值 6.变形缓冲沟：是在建筑物周围地表挖掘的一定深度的沟槽。沟深超过基础底面深200--300mm ，沟槽不小于600m ，沟外缘建筑物外侧1--2m 7.变形缝：是将建筑物从屋顶到地基分成若干长度较小，刚度较大，自成变形体系的独立单元 8安全开采上限：安全开采边界的标高 9.安全开采深度：地表至安全开采边界的距离，即地面标高与安全开采上限的标高之差 10.安全煤岩厚度：水体地面向下至安全开采上边界水平面之间的距离 11.“三下”采煤：是指在建筑物下、铁路和公路下、水体下进行开采。 12.观测站：是指在开采影响范围内的地表、岩层内部或其他研究对象上，按一定要求设立一系列互相联系的观测点。 13.土地复垦：是指对矿区土地的恢复和再利用。 14.在x 、y 两个方向或任意一个方向未达到该地质采矿条件下的充分采动尺寸时为有限开采； 16.水平变形：是指相邻两点的水平移动差值与两点之间水平距离的比值 填空 1.覆岩移动破坏形式：弯曲，跨落，片帮，岩石沿层面滑移，跨落岩石的下滑。地板岩层隆起。 2.稳态移动盆地划分为三个区域：中性区，压缩区和拉伸区 3.下沉盆地的角度参数： 边界角：开采达到或接近充分采动是，移动盆地主断面上盆地边界点和采空区边界点连续与采空区外侧水平线的夹角 移动角：开采达到或接近充分采动时，移动盆地主断面上临界变形点和采空区边界点连线与采空区外侧水平线的夹角 裂缝角：开采达到或接近充分采动时，采空区上方地表最外侧位置裂缝和采空区边界点连线与采空区外侧水平线的夹角 最大下沉角：在移动盆地的倾斜主断面上，采空区的中点与地表下沉盆地中点的连线与矿层下山方向水平线的夹角 4.启动距：地表开始移动式工作面推进距离 地表开始移动：观测地表下沉值达到10mm 地表移动时间：从地表开始移动到地表停止移动的持续时间。分为启动。活跃。衰减阶段，1.67mm/d ，百分之85 5.开采沉陷预计理论有：影响函数方法，理论模型方

开采沉陷预计方法概述

开采沉陷预计方法概述 摘要：本文主要介绍了当前使用的开采沉陷预计方法（基于实测资料的经验方法、影响函数法和理论模拟法）的原理、特点及应用情况，并简要介绍了开采沉陷预计的发展趋势，相信会对开采沉陷工作具有一定的帮助意义。 关键词：开采沉陷；预计方法；概率积分法；理论模拟法 1 引言 开采沉陷预计是矿山开采沉陷的核心内容之一，它对开采沉陷的理论研究和生产实践都有重要意义[1]。由于采矿引起的地面沉陷损坏地面建筑、公路、铁路等，不但给人民生活带来了威胁，而且破坏环境。开采沉陷的预计，对建筑物和生态环境的保护有重要意义。因此，有必要对开采沉陷预计方法进行探讨，以指导矿山的开采。开采沉陷预计方法很多，按建立预计方法的途径可分可分为三类：基于实测资料的经验方法、影响函数法和理论模拟法[2-4]。 2 开采沉陷方法简介 基于实测资料的经验方法是通过对大量的已知开采沉陷实测资料进行数据处理，确定开采沉陷中各种移动变形值的函数形式和计算预计参数的经验公式。这种方法在预测时，首先根据开采的地质条件，确定经验公式中的预计参数，再代入公式确定预计函数进而求出移动和变形值。这种方法是当前最为可靠的一种预测方法，常见的经验方法有：典型曲线法和剖面函数法等。 理论模拟法把岩体抽象为某个数学的、力学或数学-力学的理论模型，按照这个模型计算受开采影响岩体产生的移动、变形和应力的分布情况。如认为岩层和地表是一种连续的介质，则此模型属于连续介质模型；否则，就属于非连续介质模型。此法所用的函数一般均由理论研究得出，所用的参数常用实验室试验或理论推导求得，一般与现场实测资料没有直接关系，常用的理论模型法主要有连续介质力学法等。 影响函数法是介于经验方法和理论模型方法之间的一种方法，它的实质是根据理论研究或其他方法确定微小单元开采对岩层或地表的影响（以影响函数表示），把整个开采对岩层和地表的影响看作采区内所有微小单元开采影响的总和，并据此计算整个开采引起的岩层和地表的移动和变形，目前此方法中所用的参数根据实测资料获得。常用的影响函数方法有概率积分法等[5]。 下面分别对各个方法进行简单介绍。 2.1 典型曲线法 典型曲线法是用无因次的典型曲线表示移动盆地主断面上的移动和变形曲线的一种方法，它适用于矩形或者近似矩形的采区的地表移动变形预计。典型曲线法由于其分布和参数均是直接基于实测资料，因此其预计误差较小。但是建立典型曲线需要大量的观测数据，在实测数据不足的地区不能使用典型曲线法。；另外，此方法原则上只适用于矩形或近似矩形采区的地表移动和变形预计，在形状不规则的工作面开采时预计误差较大，这些限制了典型

实验五_MATLAB计算的可视化

实验五 MATLAB 计算的可视化（一） 实验目的 1. 熟练掌握MATLAB 二维曲线的绘制 2.掌握图形的修饰 3.掌握三维图的绘制 4.了解各种特殊图形的绘制 内容与步骤 1.在同一幅图形窗口中分别绘制y1=sin(t)和y2=cos(t)二条函数曲线，t 的取值范围为[0,10]。y1用红色虚线表示，y2用蓝色实线表示，横坐标轴名称为“时间t ”，纵坐标轴名称为“正弦、余弦”，整个图形的标题为“正弦和余弦曲线”。在坐标(1.7*pi ，-0.3)处添加文字“sin(t)”， 在坐标(1.6*pi ，0.8)处添加文字“cos(t)”，并在右上角添加图例，其运行界面图如下图所示。之后并尝试修改坐标轴刻度。 2.用subplot 命令在同一个窗口的不同子窗口绘制曲线y=sin(t)，y1=sin(t+0.25) y2=sin(t+0.5)，其中t=[0 10]。 3.绘制三维曲线：?? ? ??=≤≤==)cos()sin()200() cos()sin(t t t z t t y t x π （注意：用plot3命令） 4.三维网线图：绘制z=sin(y)cos(x) 三维网线图。 5. 三维曲面图 绘制22y x z +=的三维曲面图，x 在[-5,5]范围，y 在[-5,5]范围。将曲面图颜色用shading 命令连续变化，并用颜色标尺显示色图（使用函数colorbar 生成）。生成的图形如下图所示。

6.请绘制一个饼形图，数据如下表所示 7. 用semilogx命令绘制传递函数为1//(s+1)(0.5s+1)的对数幅频特性曲线，横坐标为w，纵坐标为Lw，w的范围为10-2-103，按对数分布。

矿山开采沉陷学复习题

矿山开采沉陷学》复习题 一、名词解释 1.地表移动盆地主断面：地表移动盆地内通过最大下沉点（或者说移动盆地的中心）所作的沿煤层走向的垂直断面。（P35） 2.临界开采面积：地表到充分采动时的采空区面积。 3.防砂安全煤岩柱：在松散弱含水层底界面至煤层开采上限之间为防止流砂溃入井下而保留的煤和岩层块段。（P321） 4.垮落带（冒落带）：用全部垮落法管理顶板时，回采工作面放顶后引起煤层直接顶板岩层产生破坏的范围。（P26） 5.移动角：在充分采动或接近充分采动的条件下，地表移动盆地主断面上三个临界变形值中最外边的一个临界变形值点至采空区边界的连线与水平线在煤柱一侧的夹角。（P43） 6.地表移动起动距：地表开始移动时工作面的推进距离。（P85） 7.半无限开采：工作面煤壁一侧的煤层未被采动，而另一侧的煤层全部采空的开采情况。 8. 超前影响角：将工作面前方地表开始移动（即下沉10mm）的点与当时工作面的连续，此连线与水平线在煤柱一侧的夹角。（P85） 二、填空题 三、简答题 1.什么是开采沉陷预计，其目的是什么？（P116） （1）对一个计划进行的开采，在开采进行以前，根据其地质采矿条件和选用的预计函数、参数，预先计算出受此开采影响的岩层和（或）地表的移动和变形工作，称为开采沉陷预计。（2）开采沉陷预计对开采沉陷理论的研究和生产实践都有重要意义。①在理论研究生的作用在于，利用预计的结果可以定量地研究受开采影响的岩层与地表移动在时间上和空间上的分布规律。②对指导建筑物下、铁路下、水体下的开采实践具有重要的作用。 2.岩层移动稳定后，覆岩采入影响分为哪几个带？各影响带的主要特征是什么？（P25—P27）（1）冒落带：①随着煤层的开采，其直接顶板在自重力的作用下，发生法向弯曲，当岩层内部的拉应力超过岩石的抗拉强度时，便产生断裂、破碎成块而垮落，冒落岩块大小不一，无规则地堆积在采空区内；②冒落岩石具有一定的碎胀性，冒落岩块间空隙较大，连通性好，有利于水、砂、泥土通过。冒落岩石的体积大于冒落前的原岩体积；③冒落岩石具有可压缩性；④冒落带的高度主要取决于采出厚度和上覆岩石的碎胀系数，通常为采出厚度的3~5倍。实践中可用下式近似估算冒落带高度：h=m/[(k-1)cosα]。（2）裂隙（或称断裂）带：该带内岩层产生较大的弯曲、变形及破坏，其破坏特征是，裂缝带内岩层不仅发生垂直于层理面的裂缝或断裂，而且产生顺层理面的离层裂缝。（3）弯曲（又称整体移动）带：①弯曲带内岩层在自重力的作用下产生层面法向弯曲，在水平方向处于双向受压缩状态，因而其压实程度较好，一般情况下具有隔水性，特别是当岩性较软时，隔水性能更好，成为水下开采时的良好保护层，当透水的松散层在弯曲带内就不能起到这种作用；②弯曲带内岩层的移动过程是连续而有规律的，并保持其整体性和层状结构，不存在或极少存在离层裂缝。在竖直面内，各部分的移动值相差很小；③弯曲带内的高度主要受开采深度的影响。 3.影响地表移动盆地分布规律的地质采矿因素有哪些？(P95-P104) ①覆岩力学性质、岩层层位的影响；②松散层对地表移动特征的影响；③煤层倾角的影响； ④开采厚度与开采深度的影响；⑤采区尺寸大小的影响；⑥重复采动的影响；⑦采煤方法及顶板管理方法的影响。 4.采矿引起开采空间周围岩层移动的主要形式有哪些？每种移动形式具备的条件是什么？（P23-P24） （1）弯曲：当地下煤层采出后，上覆岩层中的各个分层，从直接顶板开始沿层理面的法线

开采沉陷形成机理及其预测方法

开采沉陷形成机理及其预测方法 有用矿物被采出以后，开采区域周围的岩体的原始应力平衡状态受到破坏，应力重新分布，达到新的平衡。在此过程中，使岩层和地表产生连续的移动、变形和非连续的破坏(开裂、冒落等)，这种现象称为“开采沉陷”(Mining subsidence)。 有用矿物的开采可以是井工方法开采，也可以是露天方法开采；开采的有用矿物可以是层状的也可以是非层状的。本材料主要指的是层状有用矿物(特别是煤层)的井工开采，“开采沉陷”也是特指煤层地下开采后产生的开采沉陷。 岩体本身是一种非常复杂的介质，它不仅是出各种不同性质的岩层组成，而且还由于各种地质作用(如褶皱、断层、开裂、火成岩侵入、陷落柱等)而产生了大量的不连续面。岩体在受到各种不同开采方法的开采影响时，产生的开采沉陷是一个在时间和空间上都是非常复杂的过程。在时间上来说，在移动过程中，开采沉陷的形式和大小在不同的时间是不同的，也就是说，此时的开采沉陷是“动态的”；随着时间的推移，开采沉陷的形式和大小逐渐趋向于稳定，开采沉陷变成“静态的”或“最终的”。从空间上来说，若地下开采的范围较小、开采的矿物的埋藏深度较大，则开采沉陷波及的范围往往只局限于开采区域周围的岩体；若开采范围较大、开采矿物的埋藏深度较小，则开采沉陷波及的范围就会从岩体发展到地表，引起“地表移动”。由于人类的生产和生活活动大部分都是在地表进行，所以地表移动对人类的影响更为普遍。 第一节煤矿地下开采引起的地表移动与变形 一、地表移动的形式 所谓地表移动，是指采空区面积扩大到一定范围后，岩层移动发展到地表，使地表产生移动和变形，在地表沉陷的研究中称这一过程和现象为地表移动。开采引起的地表移动过程，受多种地质采矿因素的影响，因此，随开采深度、开采厚度、采煤方法及煤层产状等因素的不同，地表移动和破坏的形式也不完全相同。在采深和采厚的比值较大时，地表的移动和变形在空间和时间上是连续的、渐变的，具有明显的规律性。当采深和采厚的比值较小(一般小于30)或具有较大的地质构造时，地表的移动和变形在空间和时间上将是不连续的，移动和变形的分布没有严格的规律性，地表可能出现较大的裂缝或塌陷坑。地表移动和破坏的形式，

《矿山开采沉陷学》复习题

《矿山开采沉陷学》复习题 一、名词解释 1.地表移动盆地主断面：地表移动盆地内通过最大下沉点（或者说移动盆地的中心）所作的沿煤层走向的垂直断面。（P35） 2.临界开采面积：地表到充分采动时的采空区面积。 3.防砂安全煤岩柱：在松散弱含水层底界面至煤层开采上限之间为防止流砂溃入井下而保留的煤和岩层块段。（P321） 4.垮落带（冒落带）：用全部垮落法管理顶板时，回采工作面放顶后引起煤层直接顶板岩层产生破坏的范围。（P26） 5.移动角：在充分采动或接近充分采动的条件下，地表移动盆地主断面上三个临界变形值中最外边的一个临界变形值点至采空区边界的连线与水平线在煤柱一侧的夹角。（P43） 6.地表移动起动距：地表开始移动时工作面的推进距离。（P85） 7.半无限开采：工作面煤壁一侧的煤层未被采动，而另一侧的煤层全部采空的开采情况。 8. 超前影响角：将工作面前方地表开始移动（即下沉10mm）的点与当时工作面的连续，此连线与水平线在煤柱一侧的夹角。（P85） 二、填空题 三、简答题 1.什么是开采沉陷预计，其目的是什么？（P116） （1）对一个计划进行的开采，在开采进行以前，根据其地质采矿条件和选用的预计函数、参数，预先计算出受此开采影响的岩层和（或）地表的移动和变形工作，称为开采沉陷预计。（2）开采沉陷预计对开采沉陷理论的研究和生产实践都有重要意义。①在理论研究生的作用在于，利用预计的结果可以定量地研究受开采影响的岩层与地表移动在时间上和空间上的分布规律。②对指导建筑物下、铁路下、水体下的开采实践具有重要的作用。 2.岩层移动稳定后，覆岩采入影响分为哪几个带？各影响带的主要特征是什么？（P25—P27） （1）冒落带：①随着煤层的开采，其直接顶板在自重力的作用下，发生法向弯曲，当岩层内部的拉应力超过岩石的抗拉强度时，便产生断裂、破碎成块而垮落，冒落岩块大小不一，无规则地堆积在采空区内；②冒落岩石具有一定的碎胀性，冒落岩块间空隙较大，连通性好，有利于水、砂、泥土通过。冒落岩石的体积大于冒落前的原岩体积；③冒落岩石具有可压缩性；④冒落带的高度主要取决于采出厚度和上覆岩石的碎胀系数，通常为采出厚度的3~5倍。实践中可用下式近似估算冒落带高度：h=m/[(k-1)cosα]。 （2）裂隙（或称断裂）带：该带内岩层产生较大的弯曲、变形及破坏，其破坏特征是，裂缝带内岩层不仅发生垂直于层理面的裂缝或断裂，而且产生顺层理面的离层裂缝。 （3）弯曲（又称整体移动）带：①弯曲带内岩层在自重力的作用下产生层面法向弯曲，在水平方向处于双向受压缩状态，因而其压实程度较好，一般情况下具有隔水性，特别是当岩性较软时，隔水性能更好，成为水下开采时的良好保护层，当透水的松散层在弯曲带内就不能起到这种作用；②弯曲带内岩层的移动过程是连续而有规律的，并保持其整体性和层状结构，不存在或极少存在离层裂缝。在竖直面内，各部分的移动值相差很小；③弯曲带内的高度主要受开采深度的影响。 3.影响地表移动盆地分布规律的地质采矿因素有哪些？(P95-P104) ①覆岩力学性质、岩层层位的影响；②松散层对地表移动特征的影响；③煤层倾角的影响； ④开采厚度与开采深度的影响；⑤采区尺寸大小的影响；⑥重复采动的影响；⑦采煤方法及顶板管理方法的影响。 4.采矿引起开采空间周围岩层移动的主要形式有哪些？每种移动形式具备的条件是什么？

矿山开采沉陷学答案整理 2

1．“三带”的定义? 答：冒落带是指用全部垮落法管理顶板时，回采工作面放顶后引起煤层直接顶板岩层产生破坏的范围。 裂缝带：在采空区上覆岩层中产生裂缝、离层及断裂，但仍保持层状结构的那部分岩层。弯曲带：又称整体移动带，位于裂缝带之上直至地表。 2.地表移动盆地边界的确定（此题答案不确定） 一、地表移动盆地边界的划分 地表移动盆地划分成如下三个边界： (一)移动盆地的最外边界 移动盆地最外边界是以地表移动和变形都为零的盆地边界点所固定的边界。这个边界由仪器观测确定。考虑到观测误差一般取下沉为10mm的点为边界点。所以，最外边界实际上是下沉为10mm的点圈定的边界。（图中ABCD） (二)移动盆地的危险移动边界 危险移动边界是以盆地内的地表移动与变形对建筑物有无危害而划分的边界。（图中A’B’C’D’） 不同结构的建筑机能承受最大变形的能力不一样，所以各种类型的建筑物都应有对应的临界变形值。在确定移动盆地内危险移动边界时，用相应建筑物的临界变形值圈定，会更接近于实际。 (三)移动盆地的裂缝边界 裂缝边界是根据移动盆地内最外侧的裂缝圈定的边界。 3.地表移动观测站设计内容有哪些？ 答：观测站设计包括便携设计说明书和绘制设计图两部分工作。 设计说明书应包括下列内容： 1)建立观测站的目的和任务 2)设站地区的地形、地物及地质采矿条件 3)观测站设计时所用的开采沉陷参数 4)观测线的位置及长度的确定，测点及控制点的数目、位置及其编号 5)工作测点和控制点的构造及其埋设方法 6)观测内容及所用仪器，与矿区控制网的联测方法，精度要求，联测的起始数据，定期观 测时间、方法及精度要求，有关地表采动影响的测定，编录方法。 7)经费估算：包括观测站所需材料、购地、人工等费用的预算 8)观测成果的整理方法与分析步骤，所需获得的成果 4.水平煤层(或沿煤层定向主颁)非充分采动时主断面内下沉曲线特征? 答：判别：水平煤层开采时的采动程度可用走向充分采动角φ3来判别。当用φ3 角作的两直线交于岩层内部而未及地表时，此时地表为非充分采动。 （一）下沉曲线 下沉曲线表示地表移动盆地内下沉的分布规律。设沿主断面方向为x轴，下沉曲线为W(x)=F(x) 在讨论分布规律时，先要确定下沉曲线上的三个特征点： 1．最大下沉点o：下沉值最大。在水平煤层开采时，在采区中央正上方。 2．盆地边界点A、B：据走向边界角δ0作边界点A、B，此处下沉值为零。 3．拐点E：拐点是指下沉曲线凹凸的分界点。拐点从理论上讲应位于工作面开采边界的正上方，但由于工作面边界附近的顶板并不切煤壁冒落或呈阶状弯曲，存在悬顶距，因此在四周没采情况下，拐点E不在工作面开采边界的正上方而是略偏向采空区一侧。在地表达充

开采沉陷研究的意义

1绪论 1．1开采沉陷研究的意义 煤炭资源是我国的主要能源，已探明总储最在9000×10。t以上，含煤面积达 55万多平方公里，是世界上煤炭资源最丰富的国家之一。据预测．到2050年煤炭 资源仍占我国能源需求的70％眦上…，因此对地下各种赋存条件煤炭资源的开采 将是一项长期的任务。我国煤炭资源的分布也十分广泛，平原、丘陵、山区的地下 蕴藏着丰富的煤炭资源，在一些建筑物下、铁路下、水体下(简称“三下”)也压着大 量的煤炭资源．据原煤炭部1982年的不完全统计，我国仅统配煤矿的生产矿井“三 下”压煤就达137 9×10。t。我国煤炭的赋存条件也千差万别．其中急倾斜或大倾 角煤层的开采越来越受到人们关注，其原因有三个方面：一是在全国重点煤矿区有 20处100多个矿井是急倾斜(也称“大倾角”)煤层开采，急倾斜或大倾角煤炭储量 约占全国煤炭储量的15％～20％，特别是在我国的西部矿区50％阻上矿井开采的 是急倾斜煤层，如主要产煤省、市、自治区有四川、重庆、云南、贵州、新疆、甘肃、宁 夏等的急倾斜(大倾角)煤层是许多矿区或矿井的主呆煤层。随着我国西部大开发 战略的实施，我国矿产资源开采重点西移，煤炭资源开采已成为西部地区区域经济 发展的重要支柱，加强对急倾斜(大倾角)煤层开采地表沉陷机理的研究是西部煤 炭开采的重大课题。二是东部矿区赋存条件较好的煤层越来越少．加上多年来高 强度的开采，浅部赋存条件较好煤层的开采也日益枯竭，从而使急倾斜(大倾角)煤 层问题迅速进入了人们的视野，引起了人们的高度重视，如山东兖州矿区、河北邢 台和开滦矿区、安徽淮南与淮北矿区、江苏徐州矿区等的许多矿井不得不由条件相 对优越的煤层开采转向复杂的急倾斜(大倾角)煤层开采。很显然．要保持矿区生 产的高产、高效和可持续发展．也必须加强对急倾斜(大倾角)煤层开采岩层移动与 地表沉陷问题的研究。三是对急颊斜(大倾角)煤层开采岩层移动与地表沉陷问题

开采沉陷计算过程说明

开采沉陷计算过程说明 一、数据准备 1．地形图数据 （1）AutoCAD图数据 首先将AutoCAD文件备份，然后在AutoCAD中将地形图中的等高线、高程点分别设置到图层上，并逐条等高线赋高程值，高程点赋高程值，然后将其它内容全部删除，处理完后另存为AutoCAD2000格式的dxf类型文件。 注意：一定要将不是等高线、高程点的内容删除，特别是块信息须全部删除，否则会引起读入数据不正确。 （2）扫描图矢量化数据 将扫描的地形图应用矢量化软件对等高线、高程点进行矢量化，并对所有等高线、高程点赋高程值，处理完后转换成AutoCAD2000格式的dxf类型文件。（3）文本格式数据文件 地形图数据一般包括三个数据文件：线数据文件、点数据文件、文字注记数据文件。 线数据文件格式： 序号曲线开闭标志(开=1;闭=2) 等高值线型颜色值点数n X1 Y1 ……Xn Yn。数据之间用空格隔开，一条等高线数据为一行，第二条等高线数据另起一行。例： 1 2 500.0 255 6 100.0 100.0 120.0 120.0 140.0 140.0 140.0 160.0 120.0 180.0 100.0 100.0 2 1 505.0 0 8 128.284271 100.0 134.142136 105.857864 154.142136 125.857864 174.142136 145.857864 194.142136 165.857864 214.142136 185.857864 234.142136 205.857864 254.142136 225.857864 …………… 64 1 515.0 0 4 100.0 976.812409 101.593796 978.406204 121.593796 998.406204 123.187591 1000.0 点数据格式: X坐标 Y坐标 Z坐标（高程值）。例： 99973.751000 86063.929000 1022.200000 98981.372000 86069.009000 897.800000 98699.407000 86035.140000 935.900000 98412.786000 86079.170000 1000.400000 98039.797000 86055.038000 943.100000 98181.626000 86245.978000 947.000000 98436.918000 86150.296000 984.100000 98430.991000 86331.499000 935.200000 98653.683000 86258.679000 891.900000 …………… 文字注记数据格式：

实验4、matlab的计算可视化和GUI设计

p345 subplot(2,2,1) t1=0:0.1:2; y1=sin(2*pi*t1); plot(t1,y1); title('y=sin(2\pit)') 练习： subplot(2,2,2) t2=0:0.1:2; y2=[exp(-t2);exp(-2*t2);exp(-3*t2)]; plot(t2,y2) axis([0 2 -0.2 1.2]); title('y=e-t,y=e-2t,y=e-3t') 练习： subplot(2,2,3); t3=[0 1 1 2 2 3 4]; y3=[0 0 2 2 0 0 0]; plot(t3,y3); axis([0 4 -0.5 3]); title('脉冲信号') 练习： subplot(2,2,4); t4=0:0.1:2*pi; plot(sin(t4),cos(t4));

axis([-1.2 1.2 -1.2 1.2]); axis equal; title('圆') 练习： P346 x=0:0.1:20; zeta=0 y1=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y1) zeta=0.3; y2=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); hold on plot(x,y2,'r:') zeta=0.5; y3=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y3,'g*') zeta=0.707; y4=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y4,'m-') title('二阶系统曲线') legend('\zeta=0','\zeta=0.3','\zeta=0.5','\zeta=0. 707') grid on gtext('\zeta=0') gtext('\zeta=0.3') gtext('\zeta=0.5') gtext('\zeta=0.707') ginput(3) zeta = ans = 2.6037 0.9035 13.1106 2.0029 4.2166 1.0380 P347 h_fig=gcf h_axis=gca h_line1=gco h_title=get(gca,'title') h_text2=findobj(h_fig,'string','\zeta=0.3') h_fig = 1 h_axis = 151.0018 h_line1 = 1 h_title = 152.0018 h_text2 = Empty matrix: 0-by-1 set(h_line1,'linewidth',5)

matlab计算结果的可视化

第五讲计算结果的可视化 本节介绍MATLAB 的两种基本绘图功能：二维平面图形和三维立体图形。 5.1 二维平面图形 5.1.1 基本图形函数 plot 是绘制二维图形的最基本函数，它是针对向量或矩阵的列来绘制曲线的。也就是说，使用plot 函数之前，必须首先定义好曲线上每一点的x 及y 坐标，常用格式为：（1）plot(x) 当x 为一向量时，以x 元素的值为纵坐标，x 的序号为横坐标值绘制 曲线。当x 为一实矩阵时，则以其序号为横坐标，按列绘制每列元素值相对于其序号的曲线， 当x 为m× n 矩阵时，就由n 条曲线。 （2）plot(x,y) 以x 元素为横坐标值，y 元素为纵坐标值绘制曲线。 （3）plot(x,y1,x,y2,…) 以公共的x 元素为横坐标值，以y1,y2,… 元素为纵坐标值绘 制多条曲线。 例5.1.1 画出一条正弦曲线和一条余弦曲线。 >> x=0:pi/10:2*pi; >> y1=sin(x); >> y2=cos(x); >> plot(x,y1,x,y2) 图5.1.1 函数plot 绘制的正弦曲线 在绘制曲线图形时，常常采用多种颜色或线型来区分不同的数据组，MATLAB 软件专 门提供了这方面的参数选项（见表5.1.1），我们只要在每个坐标后加上相关字符串，就可实 现它们的功能。 - 2 - 表5.1.1 绘图参数表 色彩字符颜色线型字符线型格式标记符号数据点形式标记符号数据点形式 y 黄- 实线. 点<小于号 m 紫：点线o 圆s 正方形 c 青-. 点划线x 叉号 d 菱形 r 红- - 虚线+ 加号h 六角星 g 绿* 星号p 五角星 b 蓝v 向下三角形 w 白^ 向上三角形 k 黑>大于号 例如，在上例中输入 >> plot(x,y1,'r+-',x,y2,'k*:') 则得图5.1.2 图5.1.2 使用不同标记的plot 函数绘制的正弦曲线 5.1.2 图形修饰 MATLAB 软件为用户提供了一些特殊的图形函数，用于修饰已经绘制好的图形。 表5.1.2 图形修饰函数表

第七章 矿山开采沉陷预测

第七章矿山开采沉陷预测 开采沉陷预计：根据已知的地质采矿条件在开采之前预先算出地表可能产生的移动和变形叫开采沉陷预计 预计参数：指在预计函数中所用到的一系列数据 按预计方法的形式： ①剖面函数；②影响函数；③典型曲线 （一） （1）充分采动条件下地表最大下沉值 Wmax=qmcosα m——煤层法向开采厚度，mm；α——煤层的倾角；q}——充分采动条件下的下沉系数 影响最大下沉值的因素：1）采厚；2）岩性；3）倾角；4）开采方法及顶板管理方法；5）采空区尺寸大小；6）采动次数；7）采深等 （2）非充分采动条件下的最大下沉值 Wmax=qmcosα.k√n1n2 k——系数，取2~3 n1，n2——沿倾向和走向的充分采动程度系数 n1=D1/D01，n2=D2/D02 D1，D2分别是采空区沿倾向和走向的长度；D01，D02分别为地表达到充分采动时采空区的临界长度 当倾向和走向的充分采动程度系数n1，n2同时等于或大于1时，地表达到充分采动，计算时取n1=n21=1，否则为非充分采动 （二）最大水平移动值预测 在充分采动或接近充分采动条件下，最大水平移动： （1）走向方向 Umax =bWmax Umax——最大水平移动；b——水平移动系数，b=0.2~0.3 （2）倾向方向 Umax= bαWmax 表土层较薄时：bα=b+0.7cotθ 表土层较厚时：bα=b+0.7（tanɑ-h/（H0-h）），其中H0——平均开采深度；θ开采影响传播角；h——表土层厚度 ，P=0 概率积分法 作为开采沉陷的研究主体——岩层可以用两种完全不用的介质模型来模拟：一种是连续介质模型，另一种是非连续介质模型 基本假定： （1）岩体是各向同性的、均质的、不连续介质，开采引起的各方向移动与方向无关（等影响原理）：

中国矿业大学 实验六 MATLAB数据可视化

实验六MATLAB数据可视化 一、实验目的 掌握MATLAB 二维、三维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。 二、实验内容 (1) 二维图形绘制。 (2) 三维曲线和三维曲面绘制。 三、实验步骤 1．二维图形绘制 (1) 二维图形绘制主要使用函数plot。 >> clear all; >> x=linspace(0,2*pi,100); >> y1=sin(x); >> plot(x,y) >> hold on >> y2=cos(x) >> plot(x,y) >> hold off

注：hold on 用于保持图形窗口中原有的图形，hold off解除保持。 (2) 函数plot 的参数也可以是矩阵。 >> close all >> x=linspace(0,2*pi,100); >> y1=sin(x); >> y2=cos(x); >> A=[y1 ; y2]'; >> B=[x ; x]' >> plot(B,A)

(3) 选用绘图线形和颜色。>> close all >> plot(x,y1,'g+',x,y2, 'r:') >> grid on

(4) 添加文字标注。 >> title('正弦曲线和余弦曲线') >> ylabel('幅度') >> xlabel('时间') >> legend('sin(x)', 'cos(x)') >> gtext('\leftarrowsinx')

(5) 修改坐标轴范围。 >> axis equal >> axis normal >> axis([0 pi 0 1.5]) 程序如下： x=linspace(0,2*pi,100); y1=sin(x); y2=cos(x); A=[y1 ; y2]'; B=[x ; x]' plot(B,A) plot(x,y1,'g+',x,y2, 'r:') axis equal axis normal axis([0 pi 0 1.5])

矿山开采沉陷学知识点整理

矿山开采沉陷学 第一章： 1：在地下开采前，岩体在地应力场作用下处于相对平衡状态。局部矿体被采出后，在岩体内部形成一个采空区，导致周围岩体应力状态发生变化，引起应力重分布，从而使岩体产生移动变形和破坏，直至达到新的平衡。随着采矿工作的进行，这一过程不断重复。它是一个十分复杂的物理、力学变化过程，也是岩层产生移动和破坏过程，这一过程和现象称为岩层移动。 2： 充分采动区COD位于采空区中部上方，其移动特征是：煤层顶板在上覆岩体重力作用下，先向采空区方向弯曲，然后破碎成大小不一的岩块向下冒落而充填采空区。此后，岩层成层状向下弯曲，同时伴随有离层、裂隙、断裂等现象。成层状弯曲的岩层下沉，使冒落破碎的岩块逐渐被压实。移动结束后，此区内下沉的岩层仍平行于它的原始层位，层内各点的移动向量与煤层法线方向一致，在同一层内的移动向量彼此相等。 3：岩层移动形式 （一）弯曲，这岩层移动的主要形式。当地下开采后，从直接顶板开始沿层面法线方向弯曲，直到地表。 （二）岩层的垮落（或称冒落）。当煤层采出后，采空区附近上方岩层弯曲而产生拉伸变形。当拉伸变形超过岩层的允许抗拉强度时，岩层破碎成大小不一的岩块，冒落充填于采空区。此时，岩层不再保持其原有的层状结构。这是岩层移动过程中最剧烈的形式，通常只发生在采空区直接顶板岩层中。 （三）煤的挤出（又称片帮）。采空区边界煤层在支承压力作用下，一部分被压碎挤向采空区，这种现象称为片帮。由于增压区的存在，煤层顶底板岩层在支承压力作用下产生竖向压缩，从而使采空区边界以外的上覆岩岩层和地表产生移动。 （四）岩石沿层面的滑移。在开采倾斜煤层时，岩石在自重力的作用下，除产生沿层面法线方向的弯曲外，还会产生沿层面方向的移动。岩层倾角越大，岩层沿层面滑移越明显。沿层面滑移的结果，使采空区上山方向的部分岩层受拉伸，甚至剪断，而下山方向的部分岩层受压缩。 （五）垮落岩石的下滑（或滚动）。煤层采出后，采空区为冒落岩块所充填。当煤层倾角较大，而且开采自上而下顺序进行，下山部分煤层继续开采而形成新的采空区时，采空区上部垮落的岩石可能下滑而充填新采空区，从而使采空区上部的空间增大，下部空间减小，使位于采空区上山部分的岩层移动加剧，而下山部分的岩层移动减弱。 （六）底板岩层的隆起。当底板岩层较软时，在煤层采出后，底板在垂直方向减压，水平方向受压，导致底板向采空区方向隆起。

【原创】MATLAB实验报告-第二次-用MATLAB实现计算数据可视化-北京交通大学

MATLAB 上机实验报告( 2 ) 实验内容: 一、试用如下几种方法来建立向量，观察结果 ( 1) x=1:5, x=(1:5) ' 实验结果：x=1:5是行向量,x=(1:5)是列向量.且1为初始值，5为终止值，默认的步长为 1. >> x=1:5 1 2 3 4 5 >> x=(1:5)' x = 1 2

3 4 5 ( 2) x=0:pi/4:pi 实验结果：x=0:pi/4:pi 指的是x=(0,0.25*pi,0.50*pi,0.75*pi,pi). 其中pi为圆周率初始值为0,终止值为pi,步长为pi/4. >> x=0:pi/4:pi x = 0 0.7854 1.5708 2.3562 3.1416 (3)x=(0：0.2：3) ', y=e-x)p.(*sin(x) 实验结果：x的初始值为0,终止值为3,步长为0.2.而函数y表示将x向量中的每一个数代入函数y=e%x)*sin(x)得到的函数值组成的向量. >> x=(0:0.2:3)', y=exp(-x).*sin(x)

x = 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000 2.2000 2.4000 2.6000 2.8000 3.0000

0.1627 0.2610 0.3099 0.3223 0.3096 0.2807 0.2430 0.2018 0.1610 0.1231 0.0896 0.0613 0.0383 0.0204 0.0070 (4) k=linspace(-pi,pi,5), k=logspace(-3,-1,5) 实验结果：k=linspace(-pi,pi,5)产生的是初始值为-pi，终止值为 pi,元素总数为5的行向量，即k的步长为pi/2. k=logspace(-3,-1,5产生的是初始值为10八(-3),终止值为10八(-1),元素总数为5的列向量.

实验二MATLAB计算的可视化

课程实验报告 学年学期2011-2012学年第1学期 课程名称MATLAB与科学计算 实验名称实验二MATLAB计算的可视化实验室测量测绘实验中心计算机室专业年级热动113 学生姓名白治朋 学生学号2011012106 提交时间2013年10月23日 成绩 任课教师许景辉、牛亚斌 水利与建筑工程学院

实验二 MATLAB 计算的可视化 1、目的和要求 （1）熟练掌握MATLAB 二维曲线、三维图形的绘制。 （2）熟练掌握各种特殊图形的绘制。 （3）熟练掌握三维图形绘制命令。 （4）了解GUI 设计的一般过程和方法。 2、内容和步骤 参见教材实验四。 3、实验报告提交要求 （1） x=[1 2 3]，y=[1 2;2 3;5 8]，z=[2 6 9;3 8 8;1 5 7]，绘制plot （x ，y ）、plot （x ，z ），说明其各 自绘制的内容。 （2） 绘制如下图形，建立figure （2），绘图同样曲线，但标题为“你的姓名（黑体，16号字）”， 在x 坐标和y 坐标上分别标识学号和班级名称，并将网格线打开。 数组X 的列个数与矩阵y 的行个数相同， plot （ x ，y ）绘制的是x 为横坐标y 的每列为纵坐标的图像。如图1。 图1 数组X 的列个数与方阵z 的行列个数相同，plot （x ，z ）绘制的是x 为横坐标z 的每列为纵坐标的图像。如图2。 图2

（3）演示P133页，例题4.17 。

（4）完成课本P336图S 4.1实验，并用.m文件显示其程序内容。 （5）完成P302第四章例题4.

（6）通过绘制二阶系统阶跃响应，综合演示图形标识，请注释每条命令的含义。 clf; %清除图形窗口 t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t); % 数据准备 tt=t(find(abs(y-1)>0.05)); %找出符合条件(y-1)的绝对值>0.05的对应t，赋值给tt ts=max(tt); %ts为tt中最大值ts=9.6133 plot(t,y,'r-','LineWidth',3) %画曲线t-y,红色实线，线粗3磅 axis([-inf,6*pi,0.6,inf]) %设置坐标轴范围。x轴下限自动产生,上限为6*pi；y轴下限0.6,上限自动产生 set(gca,'Xtick',[2*pi,4*pi,6*pi],'Ytick',[0.95,1,1.05,max(y)]) %二维坐标刻度设置。x轴刻度线取2*pi,4*pi,6*pi,y轴取0.95,1,1.05,max(y) grid on %显示坐标刻度线 title('\it y = 1 - e^{ -\alphat}cos{\omegat}') %用斜体1书写图名 text(13.5,1.2,'\fontsize{12}{\alpha}=0.3') %图形标识，添加文字注释。在x=13.5,y=1.2处,字体大小12磅，标注α=0.3 text(13.5,1.1,'\fontsize{12}{\omega}=0.7') %图形标识，添加文字注释。在x=13.5,y=1.1处,字体大小12磅，标注ω=0.7 hold on; %保持原有图形 plot(ts,0.95,'bo','MarkerSize',10); %在x=ts，y=0.95处画蓝色的空心圆圈,大小为