地表下沉系数与采厚采深关系的讨论_毕业设计

1、摘要随着经济的发展，对能源的消耗需求日益增加。由于“三下”（尤其建筑物）压煤严重，伴随着煤炭生产规模的扩大和开采技术的提高，“三下”开采已迫在眉睫，这对生产安全提出了很大的挑战；同时，对防范灾害的意识的增强和对防护措施的加大，使人们开始重视研究采动地表动态沉陷的预计理论及方法。越来越多的新理论被应用于开采沉陷预测中来。在众多的数值模拟软件中，FLAC数值模拟软件在模拟非线性、大变形或动态问题时更加有效。本文对开采沉陷动态过程机理、预计方法等方面的研究现状进行了论述，并使用FLAC2D软件对地表下沉数值进行模拟同时与实测数据进行比较，探讨地表下沉系数与采厚、采深的关系。结论是下沉系数随采深增大而

2、减小，呈现线性表现；采厚与下沉系数呈非线性关系，随着采厚的增加，下沉系数逐渐减小。关键词：下沉系数，采厚，采深，数值模拟，开采沉陷。AbstractWiththedevelopmentofeconomy,theenergyconsumptiondemandisincreasing.Duetotheminingunder"threebodies"(especiallybuildings)Serious,withtheexpansionofthescaleandimproveoftechnology,theminingunder"threebodies"is

3、imminent,thisisahugechallengeforproductionsafety;while,withtheenhancementofdisasterawarenessandincreaseofpreventionmeasures,thepredictiontheoriesandmethodsofsurfacedynamicssubsidencebyminingbegindrawpeople'attentionincreasingly.Moreandmorenewtheorywasappliedtopredictminingsubsidence.Inmanyofthen

4、umericalsimulationsoftware,FLACismosteffectivewhenappliedtononlinearorlarge-strainproblems,ortosituationsinwhichphysicalinstabilitymayoccur.Thestudystatusofthedynamicprocessmechanismofminingsubsidenceandpredictionmethodswerediscussed,anduseflca2dsoftwaretosimulationthesurfacesubsidencedata,thencompa

5、rewiththemeasureddata,therelationshipofsurfacesubsidencecoefficientandminingheightorminingdepthwasdiscussed.TheConclusionisthesubsidencecoefficientdecreaseswiththeincreaseofminingdepth,nearlylinearchange;withtheincreaseofminingheight,thesubsidencecoefficientdecreaseswithanonlinearrelationship.Keywor

6、ds:subsidencecoefficientminingheightminingdepthnumericalsimulationsubsidencefactor.目录1绪论11.1课题研究目的和意义11.2国内外开采沉陷数值模拟研究历史和现状11.3论文研究内容32模拟矿区地质条件52.1模拟矿区地质条件52.2可采煤层顶底板条件102.3实验模拟矿区综述1.13FLAC数值模拟实验133.1FLAC软件介绍和下沉系数简介133.2数值模拟实验153.3成图结果及分析213.4模拟方案举例334总结34参考文献35致谢错误！未定义书签。附录137附录2441绪论1.1课题研究目的和意义矿山

7、资源的大规模开发和利用给人类带来了巨大的经济效益和社会效益，同时也对人类生存环境产生了一系列影响，其中，矿山开采沉陷是其中最重要的问题。矿山开采引起了地表移动和变形，直接的影响了下沉地区范围内的建筑物、河流、铁路、管道及其它构筑物，改变甚至破坏它们原有的状态。地表沉陷也会破坏地形、地貌产生，影响农业耕种、景观和矿区地面生态环境。在中国，开采沉陷使东部平原矿区的耕地大量破坏，西部矿区的水土流失和沙漠化不断加剧，南部和西南部矿区的山体滑坡现象越来越严重。这所有的一系列问题不仅给国家和企业带来了巨大损失，同时还对矿区及周边生态环境产生较大的负面影响。因此，开采沉陷的防治是一项极其重要的工作。开采沉陷

8、的防治问题不仅关系到煤炭企业能否持续发展，而且也直接影响到国民经济是否健康的进行可持续发展。因此，进行煤炭开采沉陷的研究，具有重要的理论意义和现实意义。在如何减少地表沉陷造成的损失的工作中，对地表下沉的预测是一项卓有成效的工作,其中下沉系数是表征开采地表移动规律的重要参数。然而各矿区在地质条件和采掘条件上多有差异，故下沉系数的确定也不尽相同，这就为预测地表变形增加了难度，因而确定下沉系数与煤矿开采条件的关系成为了一项重要的工作。确定更精确的下沉系数就意味着能更好的预测下沉值，就能在经济上更加有效益，在安全方面更加有保障。同时，这项工作也切实可行，对矿区安全工作影响很大，具有很重要的意义。1.2

9、国内外开采沉陷数值模拟研究历史和现状1.2.1开采沉陷的概念煤炭从地下被开采出来后，原始应力状态遭到破坏，导致应力重新分布，从而达到新的平衡状态。在此过程中，岩体产生下沉、变形、移动，word文档可自由复制编辑向上波及地表，形成地表移动盆地，从而对地物造成破坏，这就是开采沉陷。开采沉陷的分布规律取决于地质和采矿因素的综合影响。自然地质因素包括煤层厚度、倾角、埋藏厚度、松散层厚度、构造因素、物理力学和化学性质等；采矿技术因素包括采出煤量、采空区的位置、开采厚度、巷道尺寸、开采方法、顶板管理方法、开采系统、工作面推进速度、选择开采等。1.2.2国外开采沉陷研究历史与现状开采沉陷对人类的生产、生活影

10、响较突出和较早被人们认识的应该算是地表移动和地下水的流失。由于工业和经济的发展导致对能源需求的大幅度增加，“三下”开采逐渐被提上日程，在如何最大限度地开采地下矿藏，同时乂能有效地保护目标，成为人们日益关注的课题。纵观开采沉陷的形成和发展过程，可分成三个阶段：第一阶段:开采沉陷在人类开始开采地下资源时就被观察到了。早在15世纪，关于允许开采深度的问题就已被比利时人写进了法律。18世纪下半叶，人们已经能够对移动范围进行估计，这就是“法线理论”和“自然斜面理论”。第二阶段:第二次世界大战后，工业革命的迅速发展增加了对煤炭的需求，采矿业成为工业国家重要的基础产业。由于受开采技术和装备的限制，多在浅部开

11、采，因而，开采沉陷问题非常突出。这一问题受到了政府和国民的重视，许多学者开始对开采沉陷计算进行了深入探索。第三阶段：随着计算机技术的发展，岩层和地表移动计算的数学模型和技术发展很快。从上世纪70年代至今，人们将经典理论（第二阶段发展的理论）的算法编成计算机程序，使得过去难于计算的问题成为可能。有限元、边界元、离散元等三种数值方法均在开采沉陷中得到广泛的应用，开采沉陷计算正在向着自动化、智能化、可视化的方向发展，计算结果能以各种图形方式给出、直观可视。1.2.3国内开采沉陷研究历史与现状我国的开采沉陷的研究工作是在新中国成立以后开始的。在50年代初即有矿区建立地表移动监测站，在此后的时间更多的矿

12、区建立观测站，通过对这几十年的实测资料进行综合分析，初步提出了移动与变形的计算公式，以及选定有关参数的方法。借鉴外国先进技术、经验基础之上，再加上工程技术人员的创新改进，我国也逐渐形成了具有中国特色的且完全符合我国实际的一套开采沉陷理论。近年来，由于数学、力学、建筑、计算机等学科新成果不断地引入到开采沉陷学科中来，从而形成了许多新的方法和理论。在数值分析方面，应用离散元、非线性三维有限元来模拟开采沉陷的岩层移动动态过程。例如：国内于20世纪90年代初引进FLAC数值模拟软件，主要应用于工程地质和岩土力学分析，如矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等。1.2.4

13、开采沉陷学的发展方向岩层与地表移动的科学研究正处于向纵深方向发展的阶段，从研究水平和缓倾斜煤层到研究倾斜和急倾斜煤层的地表沉陷规律；从研究简单地质条件到研究复杂地质条件的移动和变形规律；从研究最终稳定的静态移动和变形到研究开采过程中动态移动和变形分布规律；从研究二维平面模型到研究三维立体模型。研究方面更加复杂化、全面化和实际化，计算机模拟也更加可视化和智能化，不仅预测结果，也能在动态过程中进行预测。1.3内文研究内容1、模拟矿区地质条件简介，通过对模拟地区的蒋庄煤矿进行介绍，进而模拟整个滕南地区。该阶段详细介绍该地区的各地层岩石性质，为后面的数值模拟提供准确的参数。2、FLAC数值模拟实验。对

14、软件和下沉系数进行原理上的介绍，并编写程序进行实验，根据地层资料，随着采深的变化改变各地层的厚度。在固定单一采深的情况下模拟不同采厚的下沉情况；在固定单一采厚的情况下模拟不同采深的下沉情况。3、分析模拟图（主要为下沉曲线图）,对下沉系数与采深、采厚的关系进行说明。根据得出的关系进行数学建模，得出数学关系式。再结合数学模型对煤矿开采边界条件的开采工作提出几点建议。2模拟矿区地质条件2.1模拟矿区地质条件枣庄矿业集团公司蒋庄煤矿位于滕县煤田的滕南矿区，行政区划届于滕州市西岗镇和微山县欢城镇。井田东西宽约3.5km,南北长约10km,面积36.6394km2。该矿设计生产能力为150万t/a,200

15、4年核定生产能力240万t/a,现年产量达270万to2.1.1井田位置枣庄矿业集团公司蒋庄煤矿位于滕县煤田的滕南矿区，位于山东省西南部，京沪铁路线南沙河站至枣庄西车站的西侧，地跨滕州、微山两县（市），西频昭阳、微山等湖，与7工苏省丰沛煤田隔湖相望。行政区划届于滕州市西岗镇和微山县欢城镇。井田东西宽约3.5km，南北长约10km,面积236.6394km2。蒋庄井田边界及四邻关系如下：1、北部西段以程楼断层，中段以蒋庄井田3870000线、杜庙断层、刘仙庄断层为界与郭庄煤矿相邻。东段以田陈井田第七勘探线为界与田陈煤矿相邻。2、南部以李集断层为界，东段与付村煤矿接壤，西段与高庄煤矿为邻。3、西部

16、北段以高庙断层、徐庄断层与柴里煤矿相邻。南段以高庙断层与崔庄煤矿接壤。4、东部以尹家洼断层为界东南段与岱庄煤矿相邻。东北段与田陈煤矿相邻。蒋庄井田内地面标局48.3733.37m,东北局西南低，井田内没有河流，仅西界的南端有旧运河流过，并与邵阳湖相接。2.1.2交通条件概况铁路以井筒位置为中心，东至京沪铁路线官桥车站13km,由官桥车站南至徐州车站82km,北至滕州车站19km。从官桥车站接轨至柴里煤矿的专用铁路线全长17km,丁1972年建成通车。官桥一柴里铁路专用线穿越蒋庄井田中部，距井筒仅110m。公路在井田中部有滕州至微山的公路穿过，井筒在该公路西侧2.1km处。井筒东至104国道15

17、km。以井筒为中心，北距西岗镇5km,经过西岗至滕州25km；南距欢城5km，经欢城东至官桥车站16km,至枣庄49km,西过微山湖二级坝到江苏省沛县34km,南过微山县城至枣庄矿业集团公司约25kmo南四湖和京杭大运河山东以南段已经疏浚完工，可以通行汽轮，从井口到微山县欢城的微山湖港大约10kmo矿井铁路、公路、水陆交通十分方便（图2.1）。图2.1交通位置示意图岩性描述柱状21.5573.40m0.602.85m。上组：顶部罗229m最大残厚度以垒色砂岩、耕砂届粘好、耕砂3上、CordaitesSchwagerinaSp3层（第四、五、六左右,是本区主40m11-12第九层石灰岩,左右,为

18、第2m左右:7m1m5m3416、:十四灰）FusilinaSp岩组柱状图青灰色石灰岩、豹皮灰岩岩溶发育,该层为井田煤层开采的间接充水含水层上部由灰、灰绿色细砂岩、粉砂岩互层、由色细砂岩、青灰色田砂岩与灰色粉砂岩互层、杂色粘土岩、红及紫灰色粉砂岩夹紫红色细砂岩组成。下部紫红、紫灰、灰绿色砾岩间夹数层紫红、紫灰色细砂岩、粉砂岩组成。图图2.尹沟组7-8-9-层兴1'刃更色、更伺色袖工、吵，吵安婚质粘土及34层砂砾层组成。下组：由粘土、砂质粘土、含砂姜结核粘土、粘土质砂及砂砾层组成。最大残余厚度为353m'根据岩性可分为上下两部。地（m）七上10-岩层厚度（m）层（16层煤全区普遍

19、发育,为稳定可17层顶板,含泥岩、常集IIII'*hTTTO憎带jj!l_IhiJ/（十二灰）岩为主'间夹灰绿、灰色中粒砂岩及紫灰色中细粒砂岩'属纯陆相沉积。由于受上侏罗统红色砂岩、砾岩的剥蚀地层保留极不完整。底部以一层分选差、矽质胶结的、浅含中粗粒砂岩或含睇岩与山西组分界。本组与太原组为连续沉积,以第一层灰岩上部之海相泥岩顶面或其相当层位为界,其上过渡相的砂质沉积物则划归山西组,本组下部为过度相沉积上部为纯陆相沉积'是本区主要含煤地层。本组底部以一层具浑浊层理和底栖动物通道的细砂岩与太原组分界。中下部以浅灰、灰白色厚层状中粒砂岩为主,夹粉砂岩、泥岩、砂质泥岩、

20、粘土岩,共含煤三层（第）。其中.41"、31岩内富含PrincipalisKoraiense_2下犀层为主要可采煤层,在泥岩、粉砂SphenopteriscftenuisSchenk、Sphenophllumthonii、Callipteridium、TaeniopterisSp、Empiectopteridiumalatum-Pecopterisarcuata等植物化石。上部为杂色泥岩、粉砂质粘土岩,粉砂岩夹灰色、灰绿色中粒砂岩。3上、3下煤层相的泥岩、砂岩、薄层石灰岩和少量粘土岩组成。含化足类、纺缝虫等动物化石,共含薄煤定可采层（16煤）,局部可采者两层（'在11层石灰岩

21、中以第三,第十层石灰岩层位稳定,厚度亦较大'为良好的标志层。本组富含下列动物化石：Dictyoclostustaiyuanfuensis-FusilinaSp、QuasifusilinaSp等。1. 第三层石灰岩至本组顶界海相泥岩：该段由泥岩,砂质泥岩,粉砂岩和薄层灰岩两层（第一、第二灰岩）组成,夹不可采:煤,偶尔出现可采点该段岩性特点是质细均匀、黑灰色、富含黄铁矿细晶和菱铁矿结核'在泥岩中常含有腕足类、海百合茎化石。2. 第三层石灰岩：一月罅度为8m左右,为探灰色,顶部色较浅,致密坚硬,含绿垂虫、珊瑚、海百合茎化石,层位与厚度均稳定'距3下煤层一般为要的标志层之一。3

22、. 第三层石灰岩至第12下煤层段：由以灰至深灰色粉砂岩为主,其间夹膨岩、溉、粘土岩和45层石灰岩（第四、五、六、-七下灰）及67层（第12下煤）组成。该段显著特点是石灰岩和煤层的层数多而厚度薄,说明是下幅度频繁箧荡环境下的沉积物一因而不利于可采煤层的堆积,除在该段有一层局部可祢层外（12下）其余煤层均不可采,且变层、尖灭现象显著。在石灰岩中以第五层石灰岩层位稳定'一般厚度为3m左右,上距第三层灰岩一般为1420m,可作为该段地层的对比标准层。4. 第层石灰岩至第十层石灰岩层段：该段以灰色至中粒砂岩为主,砂岩含较多的绿泥石，可与山西组的砂岩相区别,其次为粉砂岩、细岩岩与粉砂岩互层,夹灰岩

23、两层（第八、第九）'为太原组内的粗碎屑岩含量最高的层段含不稳定的薄层煤两层（第14、15煤）均不可采。层位稳定'其岩性为灰及棕灰色'低部质不纯,常富集纺缝虫化石碎屑,一般厚度为层的顶板,可作为该段地层对比标准层。5. 第十层石灰岩：一般厚度为由两个小层组成十上石灰岩厚一般为'质较粗糙'含动物碎屑化石,在局部地区相变为泥岩'十下石灰岩含有大量纺缝虫化石'常堆集成群'层层致密坚硬,层位稳定,厚度一般促持在'是第16层层的顶板,为本区主要标志层。6. 第16煤18煤层段：其岩性以细纱岩为主和粉沙岩、粘土岩组成,为太原组主要含煤段

24、该段含煤上、18下,其中第层煤层,第十一层灰岩为第层垂虫化石,有时相变为泥岩。为覆盖土16层,其中稳12下煤）上、15煤左右左右17、1824层石灰岩和厚3842m,由杂色泥岩、粘土岩、少量砂岩组成,顶部偶含蒲煤一层。底部为灰绿色铝土泥岩,下部为青灰、灰绿、紫红等杂色泥岩,美极不稳定的第十五层灰岩质不纯,含黄铁矿结核。中部为石灰岩（十四灰）,厚8.2513.30m'乳白灰色'致密质纯'常含灰绿色粘土薄层及粒状黄铁矿。上部为杂色泥岩'夹土质粉砂岩及不稳定的第十三层石灰岩,顶部以第十二层石灰与与太原组分界,第十二层灰岩为灰色至深灰色,顶部含粘土'呈跚犬

25、9;含有FusilinellaspPseudostaffellaSp,海百合茎化石。2.1.3矿井开采地质条件2.1.3.1地层结构特点本区地层届华北型沉积，以石炭二迭系为煤系地层，地层系统自老向新有震旦、寒武、奥陶、石炭、二迭、侏罗、第四系等。除奥陶系以前地层在煤田外围有零星出露外，余皆隐伏于第四系和侏罗系地层之下。如图2.2所示。下面对3下煤所在的二迭系及其覆盖情况叙述如下：1、二迭系（P）二迭系地层自下而上有山西组和石盒子组。（1）山西组（Pi1）该组厚度95140m,一股为110m左右。本组与太原群为连续沉积，以第一层灰岩上部之海相泥岩顶面或其相当层位为界，其上过渡相的砂质沉积物则划归

26、山西组，本组下部为过渡相沉积，上部为存陆相沉积，是本区重要的含煤地层。本组底部以一层具混浊层理和底栖动物通道的细砂岩与太原群分界，中下部为浅灰、灰白色厚层状中粒砂岩为主，加粉砂岩、泥岩、砂质泥岩、粘土岩共含煤三层（第2、3上、3下），其中3上、3下煤层为主要可采煤层。在泥岩、粉砂岩内富含楔叶羊齿、可达木、带羊齿等植物化石。上部为杂色泥岩、粉砂质粘土岩以及粉砂岩夹灰色、灰绿色中粒砂岩。（2）石盒子组（P12+P21）本组最大残厚229m。岩性以杂色砂岩、粉砂质粘土岩、粉砂岩为主，间夹灰绿、灰色中粒砂岩及紫灰色中细粒砂岩，届纯陆相沉积。底部以一层分选差、矽质胶结的、浅灰至灰绿色中粗粒砂岩或含砾砂岩

27、与山西组分界。距底部含砾砂岩2672m,一股是60m左右处含薄煤14层（柴煤），煤厚一般为1m左右，位居柴煤附近的粉砂岩和粘土岩颜色常变深，在粉砂岩内含大羽羊齿、翼羽木等植物化石，柴煤层位较稳定，可作为煤系地层上部的对比标志层。2、上侏罗统（J3）本区最大残厚为353m，根据岩性可分为上下两部。上部有灰、灰绿色细砂岩、粉砂岩互层、灰绿色细砂岩、宵灰色细砂岩与灰色粉砂岩互层、杂色粘土岩、红及紫灰色粉砂岩夹紫红色细砂岩组成。下部一般为厚80190m的砾岩，间夹数层紫红、紫灰色细砂岩、粉砂岩组成。砾岩为紫红、紫灰、灰绿色，砾岩成分以石英砾、灰岩砾为主，砾石滚圆度较好，分选较差，砾径一般1520mm，

28、胶结物一般为钙质，部分为铁质和泥质胶结，一般是以砾岩覆盖于石盒子组之上，但在部分块段亦有以砂岩与石盒子组地层相接触的。3、第四系（Q）本区第四系厚140.55181.40m,平均160.54m,区内分布东北薄、西南厚，分上下两组。上组：顶部0.602.85m为覆盖土层，其下为黄色、黄褐色粘土、砂质粘土、砂姜结核砂质粘土及3-4层砂砾层组成。下组：由粘土、含砂姜结构粘土、粘土质砂及砂砾层组成。2.1.3.2 煤层本井田主要煤系为下二叠统山西组和上石炭统太原组。含煤地层平均总厚为280m,共含煤4层，煤层平均总厚7.77m,含煤系数3.20%。山西组含煤2层，其中全区可采为3下煤层，较大面积可采为

29、3上煤层，可采总厚度为6.14m。太原组含煤2层（自上而下编号为12下、16）,其中全区可采为16层煤，局部可采为12下层煤，其余煤层均不可采，可采总厚度为1.63m。2.1.3.3 煤质井田内各煤层的成因类型届陆植煤类，变质程度为低中变质阶段的烟煤，工业牌号3上、3下煤层为气煤，12下煤层为气煤，16煤层为肥煤；3下、16煤层平均为中灰分，3上、12下煤层平均为富灰分；3上、3下煤层为低硫煤，12下煤层为中硫煤、16煤层为富硫煤；3上煤挥发分（精煤）为34.9143.29%,3下煤挥发分（精煤）为33.5845.02%,12下煤挥发分（精煤）为37.0648.36%,16煤挥发分（精煤）为4

30、1.2246.37%。2.1.3.4地质构造蒋庄井田受高庙断层、刘仙庄断层、尹家洼断层等控制，总体构造线展布为NNE向，形成一断裂构造发育，以地堑、地垒为主要特点的宽缓褶皱区。井田内断层很多，含煤地层沿走向、倾向产状均有变化，个别地段有岩浆侵入。届构造中等区。2.2可采煤层顶底板条件2.2.1 3上煤层：顶板岩性变化较大，顶板为泥岩、泥质粉砂岩，厚0.002.46m。老顶常为中、细粒砂岩，厚1425m。部分地段顶板为中细粒砂岩。砂岩单向抗压强度十燥状态下平均为60.198MPa,饱和吸水状态下平均为41.010MPa；泥岩单向抗压强度十燥状态下平均为45.449MPa。届中等冒落难冒落顶板。底

31、板多为泥岩、砂质泥岩、粉砂岩细砂岩互层，厚度一股为24m。2.2.2 3下煤层3上、3下煤层合并区内顶板为砂质泥岩、泥岩，厚0.065.25m。分义区内大部分地段顶板为浅灰、灰白色细、中砂岩，有时夹砂、泥岩互层。老顶常为中、细粒砂岩，厚度一般为2030m,当3上煤层尖灭时，可增至50m左右或更大。砂岩没有强度试验资料，据肉眼鉴定，其坚固性优于3上煤层顶板砂岩。届中等冒落极难冒落的顶板。14勘探线以北部分地段及以南的大部分地段有泥岩伪顶，其厚度多在0.50m以下，有时增厚至23m或更大些，成直接顶，容易冒落，应予以重视。底板为砂质泥岩，向下渐变为粉砂岩，砂泥岩互层，一般为35m。未发现底鼓现象。

32、2.2.3 12下煤层顶板为泥岩、砂质泥岩，向上渐变为粉砂岩、细砂岩，厚57m。稳定性较差。届冒落中等冒落顶板。底板为泥岩、砂质泥岩及粘土岩，向下为石灰岩（八灰）或砂岩，有石灰岩（八灰）存在时，泥岩厚度不超过0.5m。2.2.4 16煤层顶板为第十下层石灰岩，厚3.389.40m,一般5m左右。层位稳定，岩性坚固。据柴里井田资料，石灰岩单向抗压强度干燥状态下平均为166.557MPa,饱和吸水状态下平均为146.355MP&届极难冒落顶板。有时具泥岩伪顶，厚度0.10m左右，易随煤脱落。底板为泥岩，厚度多在1m以内，向下变为砂岩。表2.1可采煤层顶底板岩石强度一览表岩性抗压强度（MPa

33、）”强度（MPa）抗男强度（MPa）中、细砂岩581365.514.312.637.2泥岩19402.85.51623.8石灰岩541677.914.12.3实验模拟矿区综述上述为蒋庄煤矿地质报告，在此次实验中，模拟地区为整个滕南矿区，以济宁地区为中心。蒋庄煤矿地质代表着滕南矿区一般矿区的地质条件。在整个滕南矿区地形地貌特征为平原，采煤主要为3下、3下煤，还有其余的煤层例如12下、16煤层，在整个滕南矿区各部分的煤层埋藏深度和厚度各不相同。地质条件大体相近，但由丁地理中有断层，各矿区同时还具有自己的特点。在相同中还具有自己的不同。考虑到此次模拟矿区地质条件相似，但采矿条件和埋藏条件稍有不同的情

34、况，因此，在模拟的过程中采用对比方法。具体的做法是方案一：选取固定的采深，分别模拟不同采厚条件下的下沉值并根据下沉值确定数学关系式；方案二：选取固定的采厚，分别模拟不同采深条件下的下沉值并确定数学关系式。考虑到滕南矿区各矿开挖时间较长，有些矿服务时间已经很长，实测资料比较全，根据实际经验进行模拟，并与实测数据进行比对，确定模拟数据的可用性和准确性，并以此为根据进行沉陷预报模拟。根据最终的结果讨论如何在保证矿区各种公共设施和建筑物的安全的前提下创造更大的经济效益。为矿区沉陷预防方法提供一个佐证。有的模拟条件并不是真实情况,故在建模选取地质条件参数时并不是地质报告上的参数，主要是为了考虑模拟各个条

35、件的全面性，使结论更具全面性，能够适应较大的方面。但由于地质条件选取为滕南矿区，故结论的使用方面还具有局限性，只是适用于滕南矿区。3FLAC数值模拟实验3.1FLAC软件介绍和下沉系数简介3.1.1 FLAC使用原理和方法FLAC软件由美国Itascs司于20世纪80年代开发并于20世纪90年代得以广泛利用。目前该软件在国外已被广泛应用于工程地质、岩土力学以及构造地质学和成矿学等研究领域。它是基于拉格朗日差分法的一种显式有限差分程序，它不仅适宜于处理大尺度、大变形工程和地质问题，而且可以在初始模型中加入诸如断裂、节理构造等地质因素。国内于20世纪90年代初引进该软件，主要应用于工程地质和岩土力

36、学分析，如矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等。FLAC软件设计的基本算法是拉格朗日差分法，是一种利用拖带坐标系分析大变形问题的数值方法，并利用差分格式按时步积分求解。3.1.2 FLAC软件的特点FLAC应用范围广泛，拥有11种材料本构模型和5种计算模式，同时还可以模拟多种结构形式和确定任意的边界条件。下面介绍一下其余的模拟方法如有限元法、概率积分法等，1、有限元法是一种高效能、常用的计算方法.有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。自从1969年以

37、来，某些学者在流体力学中应用加权余数法中的迦辽金法或最小二乘法等同样获得了有限元方程，从而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值问题有所联系。2、概率积分法数学模型，当前，我国各矿区在进行地表移动预计时普遍采用概率积分法，其主要原因为：概率积分法采用随机介质力学理论，数学模型可靠；预计公式推导过程严密，便于进行理论分析和编程解算；对下沉盆地的描述较为精确、合理；我国学者对概率积分法进行了行之有效的发展、修改和补充工作，并且大力进行推广和应用。上述方法代表着现有的开采沉陷理论，传统的定量理论分析方法假设较多，由于地下工程的结构和岩石破坏过程十分复杂，

38、相似材料模拟成本高且误差大，这样就会影响研究结果的真实性。而数值计算由于其可以综合考虑多种因素，且计算结果直观，在地下工程稳定性研究中得到了日益广泛的应用，已经成为岩石力学的一种重要研究手段，也是分析水平承载层（如顶板）稳定性的有力工具。而且大部分现有的开采沉陷理论基本上都是以均匀连续介质假设作为理论研究前提，不能考虑岩层中存在的不连续面，如节理、裂隙以及断层的影响。而实际上岩层中存在大量的节理裂隙和规模不等的断层，这些不连续面的存在影响了开采沉陷规律，尤其是当开采区域断层比较发育时，断层对开采沉陷规律的影响十分明显，在这种情况下就不得不考虑断层等不连续面的影响；另一方面，由于目前的开采沉陷预

39、计理论（如概率积分法）的关键参数必须经过现场观测才能确定，因而给实际观测工作造成了一定的难度，因此，此次实验使用FLAC模拟软件对煤矿开采沉陷进行预计.该方法以岩石力学理论为基础，以煤岩物理力学参数和地层构造特性为计算依据,无需作任何假设或确定某些关键参数，从而克服了经典预计方法的不足。此外，它能较好地模拟地质材料在达到强度极限或屈服极限时发生的破坏或塑性流动的力学行为，分析渐进破坏和失稳，特别适用于模拟大变形。FLAC程序具有很强的前后处理功能。只要设定某些控制点的坐标，程序就可以自动生成网络，界面美观。用户或根据实际情况通过某些命令修改网络，还可以直接在屏幕上绘制或以文件形式创建和输出打印

40、多种形式的图形。还可根据需要，将若十个变量合并在同一副图形中进行研究分析。3.1.3下沉系数下沉系数（subsidencefactor是在充分采动条件下，开采近水平煤层时地表最大下沉值与开采厚度之比。其大小主要与上覆岩层的性质、地层结构、采煤方法、顶板管理方法和采动性质等地质采矿因素有关。地表下沉系数是预计地表沉陷的关键性参数，其取值的准确性决定了确定地表移动变形的精度，对地表沉陷移动预计的影响非常大。即使是相同的或相似的地质采矿条件，地表下沉机理的不同会导致下沉系数存在着很大的差异。例如：在同一矿区，利用常规方法计算出下沉系数来进行地表沉陷预计时，一些预计结果与实测结果相互之间存在着相当大的

41、误差，甚至出现错误。目前，大量的实测资料已证实下沉系数与覆岩岩性有关，岩体越坚硬，其值越小。然而，下沉系数与覆岩岩性之间的定量关系，以及是否与采深和其它因素有关等，处理不是很多，本文在构造模型时选取不同采厚和不同采深的情况，探讨下沉系数与采厚采深之间的关系，根据得出的关系进行数学建模，得出数学关系式。再结合数学模型对煤矿开采边界条件的开采工作提出几点建议。3.2数值模拟实验3.2.1数值模型建立与参数选取（1）基本思路实验考虑到研究区基岩及第四系厚松散层的变形破坏的特征和低抗拉强度性状，采用岩体力学弹塑性模型离散元分析和低拉力分析原理相藕合的方法，借助计算机模拟开采后地表及覆岩应力、应变和移动

42、规律。其基本思路为：在综合考虑影响矿区开采覆岩变形破坏特征因素的基础上，首先对工程地质体和地质环境进行深入的调查，特别是对煤（岩）层的赋存状态、物理力学性质、原岩应力场、地质边界条件的组合特征等方面的研究，概化出全面的、符合实际的地质模型。地质模型与地质原型的吻合程度是数值模拟的前提和基础。在地质模型的基础上，通过合理的抽象、简化、建立数值分析的物理模型，并进行模型的受力机制的分析研究。在此基础上，提取控制性工程地质问题的主导因素，建立数学模型，准确反映地质体的客观实际，同时乂具有力学分析的可能性和计算机条件下模拟的可行性。最后考虑岩土体的力学性质，选择合理的塑性准则、流动法则和本构关系，结合

43、工程实际实施计算。（2）物理模型此次研究计算模型主要依据所选矿区工程地质岩组柱状图（见图2.2）,在所选例子（即采深450m采厚3m）中煤层平均埋深420m,表土层厚160m,煤层厚3m。将模拟的长度和模型高度增加，以便消除边界对模型的影响程度。矿山工程岩体的应力总是处于三维空间应力状态，一般来说最好采用三维软件计算，才能较准确地确定岩体的空间应力状态，然而三维软件计算及成图复杂，对所用的软件要求高，再加上本次研究的范围很大，划分的单元多，给计算带来难度。为此，决定采用取一个具有代表性的垂直剖面，来研究煤系岩层及地表的变形破坏情况。（3）边界条件模型的左右边界上，x方向的位移较y方向的位移为小

44、，因此将水平方向的位移置为零，即给水平方向的约束，此边界条件定义为单约束边界：下部边界上，x-y两个方向上的位移都很微小，均可作为零位移边界，此边界条件定义为全约束边界；上部边界上，水平和垂直两个方向的位移都比较大，不予约束，此边界可以定义为自由边界（如图3.1）。即模型两侧表示x约束边界，即只有y方向的位移；模型下部边界表示全约束边界，即没有自由度；上边界为自由边界。H表土层匚.§i-hh各岩层h底板图3.1地层岩石模拟图此次实验分为两部分：1、下沉系数与采深关系的数值模拟；2、下沉系数与采厚关系的数值模拟。在1中分别在采厚为1m、3m、5m的情况下，选取其中一个固定的采厚，分别模

45、拟在300m、450m、600m、750m、900m采深条件下的下沉情况，观察下沉曲线图，探讨其中的规律。最终构造数学模型，并制成图表。在2中分别在采深为300m、450m、600m的情况下，选取其中一个采深固定，分别模拟在1m、2m、4m、6m、8m采厚条件下的下沉情况。同样，观察下沉曲线图，探讨其中的规律，形成公式并制成图表。在考虑到模拟实验的边界处是否受采深影响的情况下，考虑到只要在模拟的试验中地区中间部分达到最大下沉即可。因此计算模型举例如下：采深450m采厚3m的各个条件：计算条件及计算范围计算模型剖面长1000m,采长800m,两边各留长度为100m的边界。深度为450m,煤层平均

46、埋深420m，表土层厚160m，煤层厚3m。(1) 单元划分采用四边形平面单元，单元长10m,高110m。模型共划分成6000个四边形平面单元。(2) 边界条件模型两侧的边界条件为限定水平移动的滑动支撑，模型底部的边界条件为限定垂直位移的滑动支撑，模型上部为自由边界。(3) 岩石力学参数的选取岩石是一种脆性材料，当荷载达到屈服强度后将发生破坏、弱化，应届丁弹塑性体，本次试验选择摩尔一库仑准则。计算中采用的岩石力学参数与工程地质力学模型模拟试验采用的数据相同。模型中采用的岩层力学参数包括弹性模量、泊松比、粘聚力、内摩擦角、抗拉强度和密度等参数见表3.1。表3.1工作面计算模型的岩土层力学参数表序

47、号工程地质类型密度(kg/m)E/MPac/MPa抗拉强度/MPa4/(°)1土层r2400i9.6570.30010.120.00221.22粉岩250010080.2722.682.12632.53泥炭250053960.2791.411.14331.24细砂岩r260078360.23410.152.34034.55煤层14009890.2402.020.17930.26粉砂240010260.2796.121.91831.4基本假设 模拟区段内各分层内部为均匀介质。 开采模型视为平面应变模型。3.2.2模拟程序的代码编写在程序设计初，首先确定合适的网格，此次实验网格数为100

48、X60,共6000个网格。选取走向中心剖面为研究对象，岩石性质如上述，各岩层厚度依据经验设定，在程序编写的前期，首先将各岩层与网格联系起来，即建立坐标x、y与网格的关系，根据网格的疏密程度的不同，实际控制网格的坐标。其次是设定岩石性质、设定外力与边界条件，求取起始平衡并保存文件。再求工程影响的时候，首先调出应力平衡状态，在设定新的开挖word文档可自由复制编辑条件和外力、边界条件，再求取最终的应力平衡。检查最终的结果是否符合要求。不符合再调试参数继续求解。以下写出其中之一的程序代码：newgrid10060modelmgen0,00,301000,301000,0i=1101j=1,5rat=

49、1,1gen0,300,331000,331000,30i=1,101j=6,7rat=1,1gen0,330,1001000,1001000,33i=1,101j=8,14rat=1,1gen0,1000,2201000,2201000,100i=1,101j=15,28rat=1,1gen0,2200,2901000,901000,220i=1,101j=29,39rat=1,1gen0,2900,4501000,4501000,290i=1,101j=40,61rat=1,1propd=2400b=1.4e9s=5.05e8f=32.5t=2.126e6c=2.68e6j=1,盼砂pro

50、pd=1400b=9.77e8s=3.966e8f=30t=0.18e6c=2.00e6j=6,7®propd=2600b=1.29e10s=5.77e9f=34t=2.616e6c=9.04e6j=8,14砂propd=2500b=3.2e8s=1.1e8f=31.1t=1.142e6c=1.41e6j=15,2酬炭propd=2500b=3.39e9s=1.28e9f=32.7t=1.936e6c=5.41e6j=29,3盼岩propd=2000b=1.207e7s=3.71e6f=20t=0.002e6c=0.12e6j=40,6俵土层mnj=5mnj=7mnj=14mnj=2

51、8mnj=39interface1asidefrom1,5to101,5bsidefrom1,6to101,6interface1ks=14.7e6kn=588e6c=0.1e6fri=25interface2asidefrom1,7to101,7bsidefrom1,8to101,8interface2ks=14.7e6kn=588e6c=0.1e6fri=25interface3asidefrom1,14to101,14bsidefrom1,15to101,15interface3ks=14.7e6kn=588e6c=0.1e6fri=25interface4asidefrom1,28to

52、101,28bsidefrom1,29to101,29interface4ks=14.7e6kn=588e6c=0.1e6fri=25interface5asidefrom1,39to101,39bsidefrom1,40to10140interface5ks=14.7e6kn=588e6c=0.1e6fri=25fixxi=1fixxi=101fixyj=1setg=9.81appsyy=0j=61hisnstep=10HISUNBALsolvesave3m4501.savinixdisp=0ydisp=0inixvel=0yvel=0modelnulli=20,80j=6,7step820

53、save3m4502.sav以上为450m采深3m采厚的程序源代码。采用的是摩尔-库伦模型,在地层模型分为6层，分别是粉砂层、煤层、细砂层、泥炭层、粉岩层、表土层。各个地层的参数在上面程序中可以看出。各地层交界面的参数为ks=14.7e6kn=588e6c=0.1e6fri=253.3成图结果及分析3.3.1下沉系数随采深关系的模拟3.3.1.11m采厚，不同采深的下沉曲线图（图3.2-图3.6）蠕51E*皿畋；1HE避LhhM.Y-ixeYdspiMra脚iRf,蠕加图3.2300m采深下沉曲线图图3.3450m采深下沉曲线图诵IH抑-典而EJOMS¥-ue44'roira

54、EHUSHEtfT|junr正LKWl御?l>巧强E«Uw1H6EHX1-1瞬噌职LnwATlfe"Mt&uaMkeFisrnnE4D；ruE*<?图3.4600m采深下沉曲线图图3.5750m采深下沉曲线图$410H1DDOE+CQ*i«11DQEtfQ-15UE2四1EEMflCtakfwFhiUE%建皿图3.6900m采深下沉曲线图从图3.2图至图3.6中可以看出，在采厚1m的情况下，模拟结果是采深300m下沉值为0.58m,采深450m下沉值为0.51m,采深600m下沉值为0.45m,采深750m下沉值为0.38m。采深900m下沉值

55、为0.29m。可以看出，下沉系数随采深增大而减小，趋势接近线性变化。3.3.1.23m采厚，不同采深的下沉曲线图（图3.7-图3.11）510观际！MN3HE1ME部4S5芷个E$Hf«®¥«：.Y-dtvKMMnr：QQEMQfflHKi囱iat部煜茏xe四JUealA*Ul-lj3TJLM<牌Y«b.r-±¥SnWHO加也项"芸图3.9600m采深下沉曲线图图3.10750m采深下沉曲线图月单附*土峋lastXt5ap5S-IMMUtilimMi-ismwiCEMi时4on4d4M&_?acTHI

56、JDDE图3.11900m采深下沉曲线图从图3.7至图3.11中可以看出：在采厚为3m的情况下，模拟的结果是采深300m下沉值为1.62m,采深450m下沉值为1.41m,采深600m下沉值为1.19m,采深750m下沉值为1.03m。采深900m下沉值为0.85m。可以看出，下沉系数随采深增大而减小，总的趋势是线性变化。与1m采厚的模拟结果相近。3.3.1.35m采厚，不同采深的下沉曲线图（图3.12-图3.16）图3.12300m采深下沉曲线图图3.13450m采深下沉曲线图JU£nwlN|iirtfrPie*-图3.14600m采深下沉曲线图图3.15750m采深下沉曲线图Zd

57、rtt»1虹变，H诞式图3.16900m采深下沉曲线图从图3.12至图3.16中，可以得出在采厚5m的情况下，模拟的结果是采深300m下沉值为2.45m,采深450m下沉值为2.16m,采深600m下沉值为1.85m,采深750m下沉值为1.57m。采深900m下沉值为1.24m。可以看出，下沉系数随采深增大而减小，总的趋势是线性关系。与1m、3m采厚的模拟结果相近。根据以上模拟下沉曲线图可以看出下沉系数随采深增加而线性关系的减小。3.3.1.4下沉系数与采深的数据模拟结果分析表3.2下沉系数与采深的数据模拟结果下采沉,深300m450m600m750m900m1m0.58m0.51m0.45m0.38m0.29m3m1.62m1.41m1.19m1.03m0.85m5m2.45m2.16m1.85m1.57m1.24m下沉量随采深的变化大小由表3.2可见，从中可以看出下沉系数