煤矿涌水量预测知识

煤矿涌水量预测煤矿涌水量预测目的：煤矿涌水量不仅是对煤田进行技术经济评价的重要指标，而且也是设计和生产部分制订采掘方案、确定排水能力和防治措施的重要依据。任务：煤田勘探（详勘）阶段要进行首采区及第一开采水平的正常和最大涌水量预计；煤矿建设和生产阶段，要在勘探阶段涌水量预计的基础上，结合煤矿建设和生产过程中获得的观测资料进行验证对比和加以确定，对下一开采水平及水文地质条件复杂的地段单独进行巷道、采区或采面涌水量预算。在有条件时，应对有可能突水的地段进行突水量预计。煤矿涌水量预测煤矿涌水量预测是建立在定性分析基础上的定量评价。所以，预测结果是否准确可靠，在很大程度上取决于定性分析成果（如充水因素、边界条件等）。涌水量预测一般有三个步骤：建立预测煤矿涌水量的水文地质模型建立相应的数学模型代入参数解算数学模型煤矿涌水量预测涌水量预测一般有三个步骤：建立预测煤矿涌水量的水文地质模型

包括概化自然条件下的煤矿水文地质条件，如矿区含水层及隔水层的基本结构特征，以及地下水的补给、径流、排泄条件等；

根据煤矿开采方案（如开采深度、范围、方式等），确定煤矿的内边界条件；

分析未来充水因素，预测不同时期开采时的外边界条件。煤矿涌水量预测涌水量预测一般有三个步骤：建立相应的数学模型

目前，常用的数学模型分确定性和非确定性两种，确定性模型包括解析法、数值法、均衡法等，不确定性模型如经验方程（比拟法）模型和回归方程（相关分析）模型。煤矿涌水量预测涌水量预测一般有三个步骤：代入参数解算数学模型

将参数代入数学模型求解，并对解算结果进行分析评价，如发现解算结果不合理（或用不同方法解算结果相关悬殊，或与类似水文地质模型的解算结果相差甚远），应作认真检查，或修改水文地质模型或修改参数进行重新解算，直至合理时结束。相关比拟法研究涌水量与影响因素之间的数学规律为基础，建立某种可以表达这种规律的函数关系，并依此来外推未来设计疏干条件下的煤矿涌涌水量。常用的有：水文地质比拟法Q－s曲线法相关比拟法水文地质比拟法

原理：用地质、水文地质条件相似、开采方法相同的生产矿井的水文地质资料，来预测新矿井的涌涌水量。

应用条件：新、老矿井的水文地质条件要基本相似老矿井有长期的涌水量观测资料

计算方法：相关比拟法

相关比拟法

相关比拟法水文地质比拟法

计算方法关键在于涌水量与开采面积和水位降深之间的关系是否成直线，应按实际情况修改公式水文地质比拟法计算公式相关比拟法Q－s曲线法

原理：根据抽（放）水试验所获得的资料建立起来的Q－s曲线方程，预测井筒、矿井开采水平或开采地段的涌水量，反之，也可以根据涌水量预测水位降深。

应用条件：预测地区与试验地区的水文地质条件基本相似，同时，要有三个或三个以上的稳定降深和阶梯流量抽水试验资料。

计算方法：相关比拟法Q－s曲线法

计算方法：（1）分析整理抽水试验资料，一般列出这样的表格：抽水试验资料表相关比拟法Q－s曲线法

计算方法：（2）判别曲线类型，选择计算公式①曲线Ⅰ,当含水层均质、等厚且抽水试验水位降深不大，水井附近地下水运动状态保持层流时，呈直线关系；②曲线Ⅱ，在富水性强的承压含水层中进行强烈抽水时，抽水井附近水流呈紊流状态，而在离抽水进较远的地方则水流仍保持层流状态；大裂隙中的水呈紊流状态，而小裂隙中的水仍呈层流状态时，Q－s曲线呈抛物线型；③曲线Ⅲ，在地下水以储存量为主且补给来源差和导水性强的地区，水位降深小时，涌水量随降深大幅度增加；当水位降深过一定深度后，涌水量随降深增加的幅度很小，曲线有明显的下垂现象，曲线呈幂函数型，④曲线Ⅳ，对数曲线型，在富水性弱或分布范围有限，以及地下水补给贫乏且储量不大的含水层中抽水时；⑤曲线Ⅴ通常表明试验有错误或资料不可靠相关比拟法Q－s曲线法

计算方法：（2）判别曲线类型，选择计算公式

相关比拟法Q－s曲线法

计算方法：（2）判别曲线类型，选择计算公式当Q－s曲线是直线时，可以直接用公式计算。当Q－s曲线不是直线时，需要判断是何种曲线类型。这时可以用伸直法、差分法、曲度法进行判断。相关比拟法Q－s曲线法

计算方法：（2）判别曲线类型，选择计算公式伸直法相关比拟法Q－s曲线法

计算方法：（2）判别曲线类型，选择计算公式差分法一阶差分误差的大小，可用曲线拟合误差c来表示：C越小，拟合的越好。相关比拟法Q－s曲线法

计算方法：（2）判别曲线类型，选择计算公式曲度法判别式：当n=1时，为直线；1＜n＜2时，为幂函数曲线；当n=2时，为抛物线；当n＞2时，为对数曲线；如果n＜1，则抽水资料可能有误。相关比拟法Q－s曲线法

计算方法：（3）确定涌水量方程参数a、b，计算预测涌水量可以使用图解法、均衡误差法、最小二乘法确定参数a、b。图解法相关比拟法

相关比拟法实例：东庄煤矿竖井设计井深118m，预计将揭穿煤系地层30m、岩溶灰岩86m，勘探阶段在建井地段布置一水文地质孔，并进行了分层抽水试验。抽水结果表明，煤系地层含水微弱，计算涌水量时可以忽略，故未来竖井的总涌水量即为其揭露灰岩含水层的涌水量。计算步骤：解析法解析法是根据地下水动力学原理，用数学分析的方法，对一定的边界和初始条件下的地下水运动，建立定解公式，应用这些公式来预测煤矿涌水量。解析法又称地下水动力学法，它沿用了地下水动力学中的基本公式，如达西公式、裘布衣公式、泰斯公式等。在实际应用时应结合矿区的边界条件、开采条件、含水层条件，选用合适的公式进行涌水量预测。这里仅讨论边界条件、参数的选用及特定条件下的煤矿涌水量计算方法。确定边界条件矿区构造条件复杂，水文网纵横切割，使矿区边界条件形态不规则，相当复杂，为了获得理想的解析解，需要对此进行概化。1、边界进水类型根据解析法的要求，边界进水类型可以分为以下两种：隔水边界：指含水层与弱透水层、隔水层或隔水断层间的界线。在自然界中，由于绝对的隔水层是不存在的，因此常用相对隔水层的概念，即将弱透水或局部透水局部富水的岩层（这些岩层与该区含水层相比都弱得多），均作为相对隔水层处理。确定边界条件1、边界进水类型供水边界：理论上的供水边界指具有无限补给的定水头边界轮廓线，如含水层与地表水体（具有强烈水力联系）的接触线。此外，一些强含水层也可成为弱含水层的定水头供水边界。上述边界性质不是一成不变的。确定边界条件2、边界形态解析法要求将不规则的边界形态，简化为一些理想化的几何图形，如半无限直线边界、直交边界、斜交边界、平等边界。确定计算参数参数的选用直接影响煤矿涌水量预测的精度。为此，必须根据公式要求，结合矿区的水文地质条件及未来的开拓方案，合理地确定各项参数。1、渗透系数值的确定渗透系数K由抽水试验获得。在实际应用中，因为含水层的非均质性和抽水试验人为的误差，往往使求得的K值在同一含水层的不同地段差异很大，同一抽水孔中用不同方法和不同深度所获得的K值也不相同。确定计算参数1、渗透系数值的确定一般地，在抽水试验的渗流场中，都可以找到一个裘布依公式的适用区。裘布衣公式的适用区：16M≤r≤0.178R式中：M为含水层厚度，mR为补给半径（影响半径），mr为抽水孔至观测孔的距离，m在该区域测得的水位降深值代入裘布衣公式，才能求得真实的K值。确定计算参数1、渗透系数值的确定

复杂多变的地质体是很难符合上述要求的，因此，在煤矿涌水量计算时，通常采用以下两种方法利用抽水试验获得的K值。①加权平均法分以下三种情况：a）当垂直方向渗透性有变化时，如彼此之间有水力联系的几个透水性不同的砂层、砾石层或坚硬裂隙地层等，采用以下加权平均法

确定计算参数1、渗透系数值的确定①加权平均法分以下三种情况：b）沿水平各向岩石透水性有变化时，渗透系数值可由下式求得：式中：Li——不同方向渗透段的长度，m确定计算参数1、渗透系数值的确定①加权平均法分以下三种情况：c）对平面非均质情况，即含水层在水平方向上渗透性有变化时，应作渗透系数分布图，采用下式计算渗透系数：式中：Ai——某块段的面积，m2确定计算参数1、渗透系数值的确定②流场分析法它是利用抽（放）水资料绘制等水位线图，然后根据流场特征，采用闭合等值线法及辐射流法计算渗透系数。A）闭合等值线法根据达西定律，式中L1,L2——任意两条闭合等值线的长度，m

Δr——两条闭合等值线的平均距离，m Δh——两条闭合等值线的水位差，m确定计算参数1、渗透系数值的确定②流场分析法b）辐射流法采用下式计算：式中b1，b2——分别为辐射状上下游断面上的宽度，m； h1，h2——分别为上下游断面隔水底板上的水头高度，m； L——b1，b2断面之间的距离，m确定计算参数2、含水层厚度值的确定

通常含水层的厚度是依据钻孔简易水文地质观测、物探油井和抽水试验的成果分析确定的。对于岩溶充水矿井，特别是巨厚层岩溶含水层充水的矿井，

确定其厚度时，必须根据岩溶的发育规律及其富水性变化，划分出强弱含水带，决不能把整个岩溶地层的厚度作为含水层的厚度。

计算时，分别计算各区含水层的平均厚度，最后加权平均，求得矿区含水层的总体平均厚度。式中：MCP(HCP)——承压水（或潜水）含水层总体平均厚度，mAi——某一计算区的面积，m2Mi(Hi)——承压水（或潜水）含水层的厚度，m确定计算参数3、引用半径和巷道系统面积的确定

在预测巷道系统和露天采矿场的涌水量时，常把矿井的形状复杂的巷道系统（或露天采矿场的轮廓）所包围的面积，看作以r0为半径的圆形大井的面积。R0称为引用半径。由于矿井四周边界所包围面积的形状均不相同，下表罗列了几种几何形状及其r0的表达式。

巷道系统面积的确定，是用巷道系统、采区所占的水平投影面积；用设计巷道所包围地段的面积；用靠排水巷道最近一条封闭等水位（水压）线所圈定的面积。确定计算参数4、影响半径和引用影响半径的确定 1）影响半径的确定影响半径的确定方法很多，可采用裘布依公式、经验公式、抽水试验资料来确定，也可采用下述两种方法来确定：

a）直接观测法：利用地下水的天然露头或观测孔为观测点，直接确定抽（排）水时各降深的影响半径。确定计算参数4、影响半径和引用影响半径的确定 1）影响半径的确定b）作图法：在直角坐标上，将抽水孔与分布在同一直线上的各观测孔所测得的水位连接起来，并沿曲线趋势延长，与抽水前的静止水位线相交，该交点至抽水孔的距离即为影响半径。

在观测孔较多时，用本方法确定影响半径较精确。作图法确定影响半径示意图1——静止水位2——动水位3——观测水位确定计算参数4、影响半径和引用影响半径的确定 2）引用影响半径的确定a）当巷道系统所穿过的含水层原始水位近于水平时，可用下式计算：

R0=R+r0式中：R0——巷道系统的引用影响半径，m；

r0——巷道系统的引用半径，m。b）当巷道靠近地表水体时，可利用弗尔赫格依米尔公式计算引用影响半径，即：R0＝2a+r0式中：a——巷道边界线至地表水体的距离，mc）对于外形较复杂的巷道系统，当巷道系统位于地表水体附近，或者巷道分布在透水性很强的含水层、富水断层带附近时，可以用巷道边界线与地表水体或其它天然水文地质边界之间距离的加权平均值来计算引用影响半径：式中：aCP——巷道边界线与地表水体之间的平均距离，m；l——相邻两条剖面线之间的距离，m。确定计算参数5、给水度、储水系数和导水系数的确定储水系数、导水系数利用非稳定流抽水试验，通过图解法就可以获得，这里只强调一下给水度。给水度的确定一般有以下3种方法：1）对于裂隙、岩溶化含水层，可以近似用裂隙率、岩溶率代替。2）根据抽（放）水试验资料获得稳定流抽水时：式中：

V——疏干漏斗体积，可以通过绘制等降深图求得，m3；Q——抽水量，m3；Qa——补给量消耗（在抽水开始的瞬间，Qa＝0；水位稳定后，Qa=Q），m3；Qc——储存量消耗（在抽水开始的瞬间，Qc=Q；水位稳定后，Qc=0），m3。其中ΣQc可根据Q－t和s-t曲线资料通过作图求得。确定计算参数5、给水度、储水系数和导水系数的确定2）根据抽（放）水试验资料获得非稳定流抽水时：

确定计算参数5、给水度、储水系数和导水系数的确定3）根据无补给季节的动态观测资料，利用有限差分方程计算给水度计算公式的选用利用稳定流公式计算矿井涌水量无限边界完整井涌水量计算竖井、水平巷道涌水量计算倾斜巷道涌水量计算巷道系统涌水量计算有限边界完整井涌水量计算利用非稳定流公式计算矿井涌水量无限边界完整井涌水量计算有限边界完整井涌水量计算其它预测方法简介其它预测方法简介目前，在矿井涌水量预测方面，使用较为广泛的有解析法、水文地质比拟法、相关分析法和应用Benoulli能量方程进行预测等方法。而发展较快的有人工神经网络法和有限元法等预测方法。数值法原理：通过把渗透区分布有限个单元，并将这些单元内的水头函数近似地看作线性函数，然后，利用有限差分法或变分有限元法，将地下水连续方程式转化成有限个线性差分方程组或泛函数极值线性方程组，用电子计算机解这些方程组，即可得出方程式的解。有限元法方法：先给定边界条件、初始条件和水文地质参数等一组初值，看计算的水位与抽放水试验的实际观测资料是否一致。不一致则需要对边界每件和含水层参数进行修改，反复调试，直到基本一致为止，最后利用调试满意的参数和边界条件进行涌水量预测。有限元法应用：由于数值法能使考虑的数学模型更接近于实际的水文地质条件，所以其价值不低于作了许多概化的精确解。因此，自数值解问世以来，已经解决了许多解析解难以解决的复杂问题”。数值法按基于有限元或有限差分的思想，分为有限元法和有限差分法，分别通过VisualModflow、Feflow两地下水数值模拟软件实现。水均衡法原理：是应用水均衡原理来预测矿井涌水量的一种方法。它是通过对地下水动态规律的研究，建立矿区在某一期间(均衡期)地下水各收支部分间的变化关系，以建立均衡方程来预测开采地段的涌水量。方法：查明矿区范围内地下水补给、排泄条件，研究矿区在疏干过程中将要产生的变化，合理确定均衡项目和取得各项目数据，方能进行水的均衡计算。水均衡法应用：适用于均衡项目容易确定，均衡项目数据容易取得的矿区，如具有独立水文地质单元的露天采矿场或开采浅部矿体的地下巷道系统的涌水量计算。特点：不能分水平预计矿井涌水量，计算结果作为全矿最大可能涌水量的依据。人工神经网络法原理：人工神经网络(ANN)是指用大量的神经元构成的非线性系统，在一定程度和层次上模仿人脑神经系统的信息处理、存储及检索功能，具有学习、记忆和计算等智能处理功能。由于水文地质条件的复杂性，很难用准确的数学模型对矿井涌水量进行有效预测，而比较可行的方法就是建立一个矿井充水因素与涌水量之间的映射关系，人工神经网络BP（backpropagation）模型恰好体现了这方面的优点，不需要建立精确的数学模型，就能实现从输入空间到输出空间的非线性映射，因此BP网络适用于矿井涌水量预测。人工神经网络法方法：1.水文地质调查试验，确定充水水源、充水通道等影响矿井涌水量的因素。控制矿井涌突水的主要因素有充水水源和充水通道。充水水源主要包括大气有效降水（年降水量大小及季节性变化、降水性质与矿区地形、煤层埋藏与上覆岩层的透水性）、含水层水（含水层岩性、空隙性、含水层分布、厚度与补给条件）、地表水（地表水体性质与规模、地表水体与充水含水层间的水力联系和地表水体与矿井开采深度的相对位置及二者间岩层的透水性关系）和老窑水等。人工神经网络法方法：1.水文地质调查试验，确定充水水源、充水通道等影响矿井涌水量的因素。充水通道主要包括构造断裂（或岩溶发育带）、开采冒落导水裂隙带、底板隔水层卸荷裂隙带和人工导水通道等。人工神经网络法方法：2.建立人工神经网络BP模型。根据问题的实际情况建立人工神经网络模型，确定隐层神经元数。隐层神经元数的确定有2种办法，一种是通过编程，一种是利用matlab提供的神经网络工具箱（NNTOOL）。并对模型进行训练。3.运用该模型预测矿井涌水量。人工神经网络法特点：在涌水量预测中，不同类型的矿井充水因素各不相同，因此，建立普遍适用的涌水量预测模型是很不现实的，必须针对具体类型的矿井进行分析，才能具有良好的适用性。选择合适的方法进行涌水量预测目前矿井涌水量预测方法多种多样，大体上可以将它们分为确定性分析方法和不确定性(随机)分析方法两类。以上介绍的各种方法原理不同，结果表示方式不一，各有其优缺点和最佳应用环境。因而，实际工程中，应根据矿井的具体特点，最好能同时应用多种方法进行综合分析、验证，力求得出一个更加客观、可靠、合理的评价结果。选择合适的方法进行涌水量预测《化妆品术语》起草情况汇报中国疾病预防控制中心环境与健康相关产品安全所一、标准的立项和下达时间2006年卫生部政法司要求各标委会都要建立自己的术语标准。1ONE二、标准经费标准研制经费：3.8万三、标准的立项意义术语标准有利于行业间技术交流、提高标准一致性、消除贸易误差，作为标准体系中的基础标准，术语标准在各个领域的标准体系中均起着重要的作用。随着我国化妆品卫生标准体系建设逐步加快，所涉及的术语和定义的数量也在迅速增长，在此情形下，化妆品术语标准的制定就显得尤为重要。四、标准的制订原则1.合法性遵守《化妆品卫生监督条例》、《化妆品卫生监督条例实施细则》中关于化妆品的定义。2.协调性直接引用或修改采用的方式，与相关标准中的术语和定义相协调。3.科学性对于没有国标或定义不统一的术语，在定义时体现科学性的原则。4.实用性在标准体系中出现频率较高，与行业联系较紧密的术语优先选用。五、标准的起草经过

第一阶段：资料搜集

搜集国内外相关法规、标准、文献并对国外文献如美国21CFR进行翻译。第二阶段：2007年末形成初稿

初稿内容包括一般术语、卫生化学术语、毒理学术语、微生物术语、产品术语、人体安全和功效评价术语，常用英文成份术语等７部分。第三阶段：专家统稿1.2007年12月第一次专家统稿会（修订情况：1.在结构上增加原料功能术语、相关国际组织和科研机构等内容；2.在内容上增加一般术语、产品术语的种类，将化妆品行业的新产品类别纳入本标准；3.对于毒理学、卫生化学、微生物学术语进行修改；4.删除与化妆品联系不紧密、无存在必要的常用英文成分术语。2.2009年1月第二次专家统稿会会议意见：1.修改能引用国家标准的尽量引用国家标准；对存在歧义的个别用词进行修改。2.删除由于本标准中的“产品术语”一章和香化协会所制定的某个标准存在重复，因此删除“产品术语”一章的内容；对“原料功能术语”的内容进行梳理，删除了20余条内容。3.增加专家建议增加“化妆品限用物质”等若干项术语。第四阶段：征求意见2009年2月面向全国公开征求意见。第五阶段：征求意见的处理与形成送审稿。在征求意见的处理阶段再次征求了相关专家的意见。六、标准的内容依据1.《化妆品卫生监督条例》、《化妆品卫生监督条例实施细则》;2.《化妆品卫生规范》;3.美国21CFR；4.相关领域国家标准如：GB5296.3-2008消费品使用说明化妆品通用标签，GB/T14666-2003分析化学术语等；5.国内外化妆品的相关文献，如《化妆品监督管