矿井涌水量预测方法的改进及结果准确性判定.pdf

返回 相似 举报
矿井涌水量预测方法的改进及结果准确性判定.pdf_第1页
第1页 / 共8页
矿井涌水量预测方法的改进及结果准确性判定.pdf_第2页
第2页 / 共8页
矿井涌水量预测方法的改进及结果准确性判定.pdf_第3页
第3页 / 共8页
矿井涌水量预测方法的改进及结果准确性判定.pdf_第4页
第4页 / 共8页
矿井涌水量预测方法的改进及结果准确性判定.pdf_第5页
第5页 / 共8页
点击查看更多>>
资源描述:
第 48 卷第 8 期煤 炭 科 学 技 术Vol􀆱 48 No􀆱 8 2020 年8 月Coal Science and Technology Aug.2020 青年博士学术专栏 移动扫码阅读 贺晓浪ꎬ蒲治国ꎬ丁 湘ꎬ等.矿井涌水量预测方法的改进及结果准确性判定[J].煤炭科学技术ꎬ2020ꎬ488 229-236􀆱 doi10􀆱 13199/ j􀆱 cnki􀆱 cst􀆱 2020􀆱 08􀆱 029 HE XiaolangꎬPU ZhiguoꎬDING Xiangꎬet al.Improved methods for prediction of mine water inflow and determination of accuracyofresults [ J ]. CoalScienceandTechnologyꎬ 2020ꎬ 48 8 229 - 236􀆱 doi 10􀆱 13199/ j􀆱 cnki􀆱 cst􀆱 2020􀆱 08􀆱 029 矿井涌水量预测方法的改进及结果准确性判定 贺晓浪1ꎬ2ꎬ蒲治国1ꎬ丁 湘1ꎬ段东伟1 1.中煤能源研究院有限责任公司ꎬ陕西 西安 710054ꎻ2.西安科技大学 地质与环境学院ꎬ陕西 西安 710054 摘 要针对当前矿井涌水量包括正常涌水量、最大涌水量、疏放水量概念认识模糊、界定不清ꎬ涌 水量预测结果与实际水量差距较大、难以精细化以及预测结果准确性判定存在歧义的问题ꎬ综合当前 各类规范以及相关研究成果ꎬ提出了矿井涌水量的 3 重属性特征统计属性、水文属性和生产属性ꎮ 在此基础之上ꎬ探讨了矿井涌水量的基本定义ꎬ将矿井涌水量概化为充水水源与充水通道综合作用的 结果ꎬ重点结合充水通道的动态变化特征以及含水层与采动空间之间的补给关系ꎬ阐述了“大井法” 预测过程中的精细水量构成以及计算方法ꎬ提出了利用“同频率放大法”改进水文地质比拟法ꎬ得到 涌水量预测动态曲线的过程ꎻ将时间序列分析和水文地质比拟相结合提高预测步长即预测周期ꎬ 增强时间序列分析法的适用性ꎮ 比较了矿井正常涌水量、最大涌水量与预测正常涌水量、最大涌水量 在各类规范、文献中的定义与煤矿防治水实践过程中的统计差异ꎬ总结了对涌水量预测结果认识存在 分歧、忽视了涌水量预测的技术条件是当前矿井涌水量预测结果准确性判定存在的主要问题ꎻ论述了 矿井涌水量预测结果准确性判定的绝对差值法和相对差值法 2 类基本判定方法ꎬ并以涌水量预测的 主要功能为导向ꎬ提出了涌水量准确性判定的经济型、安全型、偏离型的分级标准ꎬ通过对比分析说明 了 2 类方法的优缺点ꎻ提出在矿井涌水量准确性判断的实际应用过程中应当综结合 2 类方法进行综 合判定的技术思路ꎮ 通过统一技术背景ꎬ可实现矿井涌水量的预测、实际涌水量的统计、预测结果准 确性判定在技术上的全过程闭环ꎬ进一步完善矿井水害防治的技术体系ꎮ 关键词矿井涌水量ꎻ涌水量预测ꎻ准确性判定ꎻ煤矿防治水 中图分类号TD742 文献标志码A 文章编号0253-2336202008-0229-08 Improved methods for prediction of mine water inflow and determination of accuracy of results HE Xiaolang1ꎬ2ꎬPU Zhiguo1ꎬDING Xiang1ꎬDUAN Dongwei1 1.China Coal Energy Research Institute Co.ꎬLtd.ꎬ Xi’an 710054ꎬChinaꎻ2. School of Geology and Environmentꎬ Xi’an University of Science and TechnologyꎬXi’an 710054ꎬChina 收稿日期2020-03-25ꎻ责任编辑曾康生 基金项目国家自然科学基金资助项目41272388ꎻ陕西省社会发展科技攻关资助项目2016SF-449 作者简介贺晓浪1989ꎬ男ꎬ陕西绥德人ꎬ工程师ꎬ博士研究生ꎮ E-mailhexiaolang2010@ 163.com AbstractBeing mindful of issues such as ambiguity in how to define and determine mine water inflow including normal water inflow and maximum water inflow and releasing waterꎬ the large gap between the predicted results and measured ones for mine water inflow and diffi ̄ culty for refining results plus lack of consensus on the accuracy of predicted valuesꎬ after reviewing currently available specifications and research resultsꎬ we present triple attribute characteristics for mine water inflow statistical attributeꎬ hydrological attribute and production attributeꎬ on the basis of which the basic definition of mine water inflow was discussedꎻ The mine water inflow is generalized as the result of the joint action of the water filling source and the water filling channelꎬ with special attention on the dynamic characteristics of the water filling channel and the recharge relationship between the aquifer and the mining space. This work analyzes the fine composition of water in ̄ 922 2020 年第 8 期煤 炭 科 学 技 术第 48 卷 flow and relevant calculation methods under “large well method”ꎬ and propose to employ “same frequency magnification method” to im ̄ prove hydrogeological analogy method and obtain the predicted dynamic curve of water inflowꎻ By combining time series analysis with hydrogeological analogyꎬ the prediction step size prediction periodic increased and the applicability of time series analysis was en ̄ hanced. By reviewing the statistical differences existing in all specificationsꎬ literature definitions and the flood control practice in actual mines for the normal water inflowꎬ the maximum water inflow and the predicted values for both normal water inflowꎬ the maximum water inflowꎬ it is concluded that there is lack of consensus on understanding the prediction results of water inflow and the present main issue posed in determining the accuracy in prediction of water inflow is neglecting the technical conditions .In this paperꎬ we discussed two main methods for determining the accuracy of mine water inflow predictionꎬ the “absolute difference method” and the “relative difference meth ̄ od”ꎬ proposed standards for how to classify water inflow prediction results into economic typeꎬ safety type and deviation type depending on main functions of mine water inflow predictionꎬ the advantages and disadvantages of the two methods are discussed by comparisonꎬ and we concluded that the two methods should be combined in practical application for a accurate prediction of mine water inflow. By unifying the technical backgroundsꎬ we could achieve the technical mutual verification on mine water inflow predictionꎬ real water inflow statistics and determining the accuracy of predicted resultsꎬ hence improve the technical system for mine water disaster prevention and control. Key wordsmine water inflowꎻ prediction of water inflowꎻ accuracy determinationꎻ mine water disaster prevention 0 引 言 矿井涌水量是矿井水害防治工作的一个十分关 键基础数据ꎬ伴随矿井整个生产周期ꎬ是矿井生产系 统特别是防排水系统设计及其水害防治技术方案与 技术路线选择的基础[1]ꎮ 在矿井水害防治技术进 程中ꎬ形成了多种关于矿井涌水量预测计算的理论 和方法[2-3]ꎬ随着绿色、高效矿井建设的持续推进ꎬ 煤矿对于涌水预测的结果精细度、预测方法的适用 性等要求也在逐步提高ꎮ 而相关理论研究和防治水 实践中关于矿井涌水量的定义、预测方法有待改进ꎬ 缺乏完善的矿井涌水量预测结果的准确性判定方 法ꎬ难以实现矿井涌水量的预测、实际涌水量的统 计、结果准确性判定等涌水量计算全过程的闭环ꎬ进 而提升涌水量预测的技术水平ꎮ 笔者针对以上问题ꎬ综合当前各类规范以及 科研成果和前人的论述ꎬ探讨矿井涌水量的定 义ꎬ阐述矿井涌水量的主要计算方法在实际应用 中的改进方案ꎬ分析矿井正常涌水量、最大涌水 量在各类定义和煤矿实践过程中的统计差异ꎬ论 述矿井涌水量预测结果准确性判定的基本方法ꎮ 从而促进矿井涌水量预测方法、理论的不断完 善ꎬ预测精度的持续提高ꎬ为煤矿的安全、绿色开 采提供支撑ꎮ 1 矿井预测涌水量与实际涌水量差异性 矿井涌水量一般分为正常涌水量和最大涌水 量ꎬ煤矿防治水细则 [4] 中以下简称细则对 矿井正常涌水量的定义为指矿井开采期间ꎬ单位时 间内流入矿井的平均水量ꎮ 一般以年度作为统计区 间[4]ꎮ 该定义强调了涌水量的时间段特征ꎬ为平均 水量ꎮ 矿井最大涌水量ꎬ是指矿井开采期间ꎬ正常情 况下矿井涌水量的高峰值[5]ꎮ 最大涌水量主要与 采动影响或降水量有关ꎬ不包括矿井灾害水量ꎬ一般 也以年度作为统计区间取最大值ꎮ 而对于矿井预测 涌水量也是如此ꎬ将正常涌水量和最大涌水量分别 预测ꎮ 不同之处在于预测正常涌水量和预测最大 涌水量针对的是某一生产阶段或时间节点ꎬ而且为 了保障所预测的水量包含了整个生产过程ꎬ各项参 数会取在阶段性生产过程完成时ꎬ采用各类计算方 法得到的涌水量值ꎬ一般多属于预测生产阶段内正 常涌水量的峰值ꎬ尤其是解析法计算的涌水量ꎮ 但 该预测值属于瞬时量ꎬ而非细则中计算的平均量 如矿井实际涌水量处在稳定增长阶段ꎬ则平均量 接近生产阶段中期的水量ꎬ强调了生产的空间过 程ꎮ 由于细则定义的矿井正常涌水量强调了时 间过程1 年内ꎬ而预测的正常涌水量则强调了空 间过程某个水平或者盘区、工作面ꎬ两者一般不 能完全对应起来图 1ꎬ表 1ꎮ 此外ꎬ预测最大涌水量一般在预测正常涌水量 的基础上乘以一个最大系数所得一般取 1.2 3.0ꎮ 其意义在于计算出该生产阶段涌水量的上限 值ꎬ目的在于保障计算误差处在矿井可接受的范围 之内ꎮ 预测最大涌水量也属于瞬时涌水量ꎬ而与 细则中所定义的某一个时间段内的峰值量的意 义不同ꎬ则有 矿井涌水量一直处在稳定增长阶段ꎬ 则 2 个数值才有可能重合ꎻ矿井实测涌水量的依据 为井下涌水量台账ꎬ每个生产节点都只会有一个涌 水量数据ꎬ而该数据是唯一的ꎬ该节点不可能存在有 最大涌水量ꎮ 但从整个年度来看ꎬ会有最大涌水量ꎬ 该涌水量数值与细则中定义的最大涌水量值相 对应ꎮ 032 贺晓浪等矿井涌水量预测方法的改进及结果准确性判定2020 年第 8 期 图 1 蒙陕地区某矿井涌水量历时曲线 Fig.1 Duration curve of water inflow from a mine in Inner Mongolia-Shaan Area 表 1 蒙陕地区某煤矿不同定义下矿井涌水量对比 Table 1 Comparison of mine water inflow under different definitions of a coal mine in Inner Mongolia-Shaan Area 年 份2016 年2017 年2018 年 细则 定义矿井正常涌水 量/ 最大涌水量/ m3􀅱d -1 750/805 1 130/1 480 1 394/1 500 预测的矿井正常涌水量/ 最大 涌水量12 月底 / m3􀅱d -1 610/860 1 400/1 530 1 325/1 620 实际矿井涌水量/ 最大涌水量 12 月底 / m3􀅱d -1 650/805 1 480/1 480 1 373/1 500 综合来看ꎬ矿井涌水量的预测与实践ꎬ必须以涌 水量台账为根本ꎮ ①涌水量台账的记录应尽量详 细ꎬ除了记录观测时间外ꎬ还应详细记录各观测点的 位置以及采掘进尺变化情况ꎻ②涌水量预测前ꎬ引用 涌水量台账过程中ꎬ在利用“大井法”和“廊道法”等 各类确定性分析方法时ꎬ应以不同回采距离对应的 瞬时水量台账值为准ꎬ而非细则中所定义的矿井 涌水量值统计意义上的数值ꎬ忽略了涌水量数据 随回采距离变化的规律性ꎮ ③对于煤矿来说ꎬ从 年度角度出发ꎬ可以参考细则进行正常涌水量和 最大涌水量的统计ꎬ但在具体实践过程中ꎬ应当以台 账为基础ꎬ结合当前矿井所处的开采阶段ꎬ来阐述矿 井的涌水情况ꎮ 2 矿井涌水量的概念与计算方法 2.1 矿井涌水量概念 在煤矿防治水细则中将矿井涌水量定义为 矿井开采期间ꎬ单位时间流入矿井的水量ꎮ 虎维岳 等[1]认为以上定义对矿井涌水量概念定义不清ꎬ使 得预测结果对生产实际的指导性不足ꎬ并重新提出 了矿井涌水量的定义“在矿井建设或生产过程中ꎬ 矿井采掘工程直接或间接揭露或导通某种充水水源 的水ꎬ使其在特定时段和区段单位时间内流入矿井 的总水量ꎬ特定的时间段是指矿井所处的特定生产 期间如井筒建设期、首采面回采期、一水平开采期 等ꎬ特定的区段是指矿井某个生产区域如水平、 采区、工作面等”ꎮ 主要特点为①强调了矿井涌 水量是因矿井采掘活动诱发的ꎬ充水水源的水自然 地、无法控制地流入矿井ꎬ没有把人为的、通过专门 的疏放水工程疏放出来的水量归结为矿井涌水量ꎮ ②强调了流入矿井的水总量是特定生产阶段的对应 水量ꎬ因此必须是个一个动态变化值[6]ꎮ 以上论述 值得商榷在生产实际的涌水量统计、预测过程中ꎬ 一般矿井涌水量主要为排水系统设计及水害防治提 供直接依据ꎬ对于疏放水量和井下生产用水等人为 流入水量ꎬ如果不计入矿井涌水量ꎬ则可能对矿井的 系统设计造成严重的误判比如ꎬ部分煤矿的井下 疏放水水量超过其他所有水量之和ꎬ进而威胁矿 井安全ꎮ 此外ꎬ矿井疏放水的目的、方法以及结果评 价均属于矿井水害防治的范畴ꎬ不将疏放水量纳入 到矿井涌水量之中ꎬ难以真正使得涌水量统计、预测 工作有效服务于矿井安全生产ꎮ 笔者认为ꎬ矿井涌水量的定义ꎬ应当从服务生产 实际出发ꎬ不应仅考虑矿井的水文地质特征以及采 掘实际ꎬ而忽略了疏放水工程、井下生产用水等产生 的水量ꎮ 与此同时ꎬ矿井涌水量的概念应该是从发 生事实角度的定义ꎬ其还有统计意义分析涌水量的 总体水平及矿井水文地质条件ꎬ总结水量预测结果 的准确性ꎬ优化排水系统布置ꎬ与矿井涌水量预测有 着十分明显的区别由于矿井的生产是个长期、缓慢 的发展过程ꎬ因此每个矿井涌水量必然对应着特定 的时间和生产阶段ꎬ也导致了矿井涌水量本身具备 了特定时间、特定生产阶段的属性ꎮ 相应地ꎬ矿井涌 水量预测则面向的是未来的长期过程ꎬ在进行预测 时尤其需要考虑或强调采掘阶段性特征来满足生产 实践ꎮ 通过这两者的计算可做明显对比矿井涌水 量计算当前主要通过井下水泵房的排水记录来推断 或者井下各涌水点的水量实测值相加ꎬ其目的还在 于进行水量的总体统计分析ꎻ而矿井涌水量预测的 计算方法则包含了解析法、数值法、水文地质比拟法 等各类方法ꎮ 但综合来看ꎬ虎维岳等对于矿井涌水 量的定义更多是结合了矿井涌水量预测的要求ꎬ而 非矿井涌水量本身意义ꎮ 进而从广义属性特征来 说ꎬ矿井涌水量应该同时具备统计属性、水文属性和 生产属性ꎮ 基于以上分析ꎬ笔者对矿井涌水量进行 了重新定义矿井涌水量是指矿井在某一生产阶段ꎬ 单位时间内井下各出水点涌水量总和ꎮ 该定义方式 可以体现涌水量的统计属性之外ꎬ也说明了矿井的 132 2020 年第 8 期煤 炭 科 学 技 术第 48 卷 涌水量是在不同生产阶段或者不同生产区域的水量 叠加的总和ꎬ兼有了上述 3 重属性ꎮ 2.2 矿井涌水量预测计算方法改进 如前文所述ꎬ矿井涌水量是井下各个涌水点的 水量之和ꎬ矿井涌水量预测实际上是结合生产规划ꎬ 在分析矿井涌水量各组成的前提下ꎬ运用水文地质、 数学地质方法ꎬ针对井下一个或多个涌水点ꎬ开展水 量预测ꎬ并将水量结果进行叠加分析ꎬ进而得到矿井 涌水量的过程ꎮ 因此ꎬ下文所述的矿井涌水量预测 方法皆是针对涌水点重点是回采工作面的计算 方法ꎬ不再赘述ꎮ 2.2.1 涌水量预测的基本要求 为了使预测涌水量也能够具备统计属性、水文 属性和生产属性ꎮ 提出涌水量的预测结果应该包含 了 3 个数值①矿井预测正常涌水量ꎬ取该生产阶段 结束时的预测矿井涌水量ꎻ②矿井预测最大涌水量ꎬ 一般取矿井预测正常涌水量的 1.23.0 倍ꎬ具体倍 数值可参考以往最大涌水量与实际正常涌水量的比 值确定ꎻ③矿井预测平均涌水量ꎬ取该生产阶段内典 型节点年初、年终、年底处进行正常涌水量预测ꎬ 结合生产进度、充水通道的变化特征以及含水层与 采动空间之间的补给关系进行涌水量的精细预测ꎬ 绘制涌水量预测曲线ꎬ利用该曲线求取矿井的平均 涌水量ꎮ 2.2.2 矿井涌水量预测改进方法 通过统计多个煤矿的工作面和开采年代久远的 矿井涌水量曲线发现ꎬ工作面的最大涌水量多不会 出现在工作面回采结束或者闭坑时ꎬ而这与传统的 水文地质理论计算结果相矛盾无限边界条件下ꎬ采 动影响的范围越大ꎬ则涌水量越大[7-8]ꎮ 为此ꎬ提出 以下涌水量预测的改进方法以使得涌水量预测结果 满足实际需求ꎮ 2.2.2.1 解析法 将矿井涌水量概化为充水水源与充水通道综合 作用的结果ꎬ重点结合充水通道的动态变化特征以 及含水层与充水通道间的补给关系ꎬ说明“解析法” 的改进方法ꎮ 1“大井法”的理论基础是井壁四周进水模型ꎬ 概化为抽水井ꎬ水井内部不含水[9]ꎮ 但作为井下主 要涌水点的回采工作面ꎬ当前现代化高效回采方式 下ꎬ单个工作面的走向长度一般为 2 0006 000 mꎬ 倾向长度 120400 mꎬ因此工作面内部的静储量 Qj 无法忽略ꎮ 而大井法实际上计算的是工作面回采后 的动态补给量 Qd1ꎮ 此外ꎬ对于导水裂隙带未直接 导通但其与直接充水含水层间无稳定隔水层的巨厚 含水层ꎬ其垂向动态补给量 Qd也不可忽略ꎮ 这部 分水量可以按照“顶板+井壁”的出水模式计算[10]ꎬ 顶板出水面概化为圆形ꎬ通过圆形补给周长与垂向 渗透系数、水位降深之积得到 Qd2ꎮ 则工作面的实 测涌水量应该是Q=Qj +Q d1 +Q d2ꎬ其中静储量 Qj 以 重力释水为主ꎬ可以用重力给水度、开采冒落面积和 充水含水层厚度的乘积求得ꎻ导水裂隙带直接导通 范围内的含水层厚度按照实际厚度计算ꎬ而对于未 直接导通但中间无稳定隔水层的其他部分含水层厚 度ꎬ也可引入抽水“有效带”的概念[11-12]ꎬ将非完整 井群涌水量计算转化为等效完整井群涌水量计算ꎮ 2对于“大井法”而言ꎬ关键在于确定参数ꎬ尤 其是渗透系数 K 和含水层厚度 Mꎬ这 2 个对涌水量 计算结果敏感度最高的参数[13]ꎮ 一般渗透系数是 通过抽水试验或者井下放水试验得到ꎬ但地面抽水 试验得到的参数结果普遍偏小ꎬ而井下放水试验的 结果却普遍偏大ꎮ 为此ꎬ为了能够得到较为准确的 Kꎬ应该以实际回采过的工作面涌水量进行 K 的反 算ꎬ并结合地面和井下的水文孔进行综合计算ꎬ而该 方法得到的数值较为准确ꎮ 对于未进行过工作面回 采的矿井ꎬ可以参考周边矿井的反算渗透系数和本 矿井的水文试验数据综合得到ꎮ 如前文所述ꎬ不同 的导水裂隙带发育高度下ꎬ对应着不同的涌水模式ꎮ 但对于一个工作面来说ꎬ导水裂隙带是逐步发育的ꎬ 因此含水层的厚度值也应当是动态变化的ꎮ 如一个 倾向长度 300 m 的工作面ꎬ其导水裂隙带一般在一 次见方附近发育到顶界面ꎬ因此对于前 300 m 的涌 水量预测ꎬ应该通过预测导水裂隙带的动态发育特 征来 确 定 不 同 的 充 水 含 水 层 厚 度 值 来 进 行 计算[14-15]ꎮ 2.2.2.2 水文地质比拟法 水文地质比拟法属于当前矿井涌水量预测中最 为简单、直接的一类方法ꎬ可以规避大量的参数计 算ꎬ默认部分参数相似或者相等来进行计算[16]ꎮ 对 于开采程度较高的矿井ꎬ预测的准确性也较高ꎮ 除 了常用的“富水系数比拟法”、“单宽流量比拟法”等 比拟方法外[17]ꎬ还可以用以下方法进行计算 1相邻次采工作面“同频率放大法”ꎮ 对于次 采的相邻工作面涌水量精细预测ꎬ由于难以直接判 定首采工作面回采对次采工作面的影响ꎬ因此预测 的精准程度较低ꎮ 可以参考工程水文学中对于“洪 水过程线”计算时采用的“同频率放大法” [18] 进行 计算依据首采工作面的单宽流量双侧进水比拟 计算出未采面的总体水量单侧进水ꎮ 将 2 个水 量的比值记作 kꎬ即 232 贺晓浪等矿井涌水量预测方法的改进及结果准确性判定2020 年第 8 期 图 2 巨厚含水层“顶板+井壁”涌水模式[10] Fig.1 Water irushing model diagram of “roof + sidewall”in thick aquifer[10] k = Q2/ Q11 式中Q1为首采工作面采空涌水量ꎬm3/ hꎻQ2为次采 工作面采空区涌水量ꎬm3/ hꎮ 在此基础之上ꎬ将首采工作面的涌水量曲线放 大缩小k 倍ꎬ即可得到次采工作面涌水量与进尺 的曲线图ꎮ 2远距离次采工作面多因素综合比拟法ꎮ 比 拟计算的前提是 2 个工作面参数相近ꎬ但受到回采 以及工作面距离远的影响ꎬ工作面的水文地质参数 发生变化ꎬ尤其是含水层的水位、厚度与已采工作面 有较大差异ꎬ因此还需要考虑多因素进行比拟[19]ꎮ 多因素比拟所采取的公式为 Q = Q0 F F0 X S S0 Y 2 式中Q 为待采工作面的预测涌水量ꎻQ0、F0、S0为已 知的工作面涌水量、开采面积和相应的水位降深ꎻ F、S 为预测工作面的开采面积和水位降深ꎻX、Y 为 待定系数ꎬ据已采工作面涌水量进行拟合确定ꎮ 2.2.2.3 时间序列分析法 其基本思想是利用分段线性化的方法来处理非 线性动态数据[20]ꎮ 由于门限值的控制作用ꎬ门限自 回归模型对于模拟类似跳跃、共振等非线性特征有 较强大的功能ꎬ适用于开采年代久远的矿井涌水量 预测[21]ꎮ 涌水量的时间序列{Qt}自相关可用如下 数学模型表示为 Qt =a j 0 +∑ kj i=1ai jQ t-i +ε j t rj-1 Qt-d rj j = 1ꎬ2 ꎬ􀆻ꎬl3 式中rj为门限值ꎻd 为延迟步长ꎻkj为第 j 个分段的 阶ꎻa j 0 、a j i 均为模型的系数ꎻε j t 为第 j 个门限自 回归模型的白噪声ꎻl 为模型的个数ꎮ 但受限于ꎮ 因上述方法受限于模型的预测步长 SETAR 模型的预测步长不能超过延迟步长 dꎬ通 过自激励门限回归建立的模型并不能够预测较长时 间段内的矿井涌水量ꎮ 而是对于矿井的生产建设来 说ꎬ必须要将短期及时性预测与长期总体性预测相 结合ꎬ才能更好的指导矿井的安全生产ꎮ 通过有效利用矿井多年来的涌水量资料ꎬ可采 用水文地质比拟法进行矿井涌水量的长期总体性预 测ꎻ应用自激励门限自回归进行矿井即时预测ꎬ两者 配合使用ꎬ实现对矿井涌水量预测周期的延长和预 测精度的提高[22-23]ꎻ该预测方法可以有效减少一般 回归分析对于其他相关因素数据的依赖ꎬ更适应于 水文地质条件复杂或者水文资料相对匮乏矿井的涌 水量预测ꎮ 3 矿井涌水量预测准确性判定 在进行涌水量预测过程中难以避免地要出现模 型误差、观测误差、截断误差、舍入误差[24]ꎮ 导致预 测结果与实际涌水量存在一定的偏差ꎬ对这种偏差 或者涌水量准确性的判定ꎬ不仅可以评价相关预测 成果的科学性和可靠程度ꎬ更是煤矿及相关科研单 位不断优化科研成果ꎬ开展后续防治水工作的基本ꎮ 3.1 当前涌水量预测准确性判定存在的问题 如前文所述ꎬ当前对涌水量的认识存在分歧以 及忽视技术背景条件使得对矿井预测涌水量准确性 判定中存在明显不足ꎮ 3.1.1 涌水量预测结果认识存在分歧 煤矿要判定某个回采工作面的涌水量预测结果 是否准确ꎬ当前判定方法多如下将该工作面预测最 大涌水量解析计算以及比拟法计算的最大涌水量 多出现在回采结束时与整个工作面回采过程中的 最大涌水量进行对比该最大涌水量一般多不出现 在回采结束时ꎬ图 3、图 4ꎬ如果两者相差较小ꎬ则认 为计算较为准确ꎬ反之则认为计算结果差距较大ꎮ 这样的比较方法ꎬ没有将涌水量预测与实际水 量的预测固定在同一参数背景下同一个生产节 点ꎬ也使得矿井涌水量的预测结果准确性判定存 在有明显的错误ꎬ对于矿井防治水工作存在严重的 误导ꎮ 3.1.2 忽视了涌水量预测的技术条件 经过多年的技术发展ꎬ矿井涌水量预测已经朝 着精细化方向发展ꎬ也使得矿井涌水量预测多以煤 332 2020 年第 8 期煤 炭 科 学 技 术第 48 卷 图 3 某矿 3101 工作面涌水量历时曲线 Fig.3 Duration curve of water inflow in No.3101 working face 图 4 某矿 201 工作面涌水量历时曲线 Fig.4 Duration curve of water inflow in No.201 working face 矿的采掘规划为依据ꎬ这个采掘规划包含了回采区 的范围、面积、接续方法、回采方式等技术参数ꎮ 而 由于井下施工条件的复杂性ꎬ回采范围以及接续方 式等涌水量计算的基础参数可能发生变化ꎬ但煤矿 易忽略技术条件的变化ꎬ此时依然用改变前的预测 结果去与现有技术参数下的实际涌水量进行对比ꎬ 进而使得预测结果多与实际差距较大ꎮ 3.2 涌水量预测结果准确性判定方法 文献[25-26]依据 GB1521894地下水资源 分类分级标准与 GB 1271991矿区水文地质工 程地质勘探规范ꎬ根据不同的勘查研究程度和计 算方法求得的矿井涌水量ꎬ可以认定为不同的精度 级别ꎮ 但其针对的是勘探阶段的矿井涌水量预测ꎬ 预测的范围多为整个井田或者盘区ꎬ难以真正有效 指导生产ꎮ 而当前的生产实践要求涌水量预测须结 合不同阶段的采掘规划进行精细计算ꎬ并在同一生 产节点下ꎬ与实际水量值进行对比分析ꎬ从而确定同 一生产阶段涌水量预测值的准确性ꎮ 当前的准确性 对比判定主要有如下 2 种方法ꎮ 1涌水量绝对差值判定法按照预测的正常涌 水量的与实际涌水量差值的绝对值进行确定ꎮ 如某 矿在某一具体的生产阶段涌水量为Qsꎬ而对应生产 阶段的预测涌水量结果为 Qyꎬ则 Qc =Q y -Q s为预测 涌水量 Qy的绝对差值ꎮ 当 Qc0 时ꎬQc为盈预测量或者 强预测涌水量ꎮ 而为了保障矿井的排水安全ꎬ多出 现了的是盈预测量或强预测量ꎮ 而对于具体绝对差 值为多少才算是偏差大ꎬ难以指导实践ꎬ尚没有可以 确定的标准ꎮ 主要在于用绝对差值来判定预测涌水 量的准确程度存在较大的局限性ꎬ在不同水文条件 的矿井ꎬ其难以具备普遍适用的特征ꎮ 如实际涌水 量为 1 000 和 100 m3/ h 的 2 个矿井ꎬ若其预测涌水 量和实际涌水量绝对差值均为 10 m3/ hꎬ但显然前 者的准确程度要高于后者ꎮ 由此可知ꎬ决定涌水量 预测准的确程度ꎬ除了要计算绝对差值ꎬ还有考虑该 矿井实际涌水量本身的大小ꎬ进而提出涌水量的相 对差值判定法ꎮ 2涌水量相对差值判定法该方法按照预测的 矿井涌水量和实际涌水量差值与实际值的比值确 定ꎮ 如某矿在确定的生产阶段的涌水量为Qsꎬ而对 应生产阶段的预测涌水量结果为 Qyꎬ则 Qr= Qy- Qs / Qs为预测涌水量 Qy的相对差值ꎮ 根据上述定 义可知ꎬ实际涌水量为 1 000 和 100 m3/ h 的 2 个矿 井ꎬ预测涌水量和实际涌水量绝对差值都为 10 m3/ h 的预测结果ꎬ前者的相对误差为 1%ꎬ而后者的相 对误差为 10%ꎬ可见前者的准确程度要远高于后 者ꎮ 类似与绝对差值限ꎬ为了使得涌水量预测结果 可以指导实践ꎬ须确定相对差值限或者相对差值上 限 δꎬ使得 | Qr| =| Qy - Q s / Qs| ≤δꎮ δ 越小ꎬ则 预测的精度越高ꎮ 在实践中如何确定 δ 至关重要ꎮ 笔者提出以矿井排水安全为根本ꎬδ 判定方法为考 虑到煤矿设置排水系统时要求能在 20 h 内排出 24 h 的水ꎬ即排水系统设置仅高于预测水量值的 20%ꎮ 因此涌水量的预测相对差值最大不超过- 20%ꎬ否则威胁矿井安全ꎮ 考虑 0.5 的保险系数ꎬ则 认为该涌水量预测值弱相对差值最大不超过 - 10%ꎬ而强相对差值则主要考虑涌水量预测过大带 来的防治水工程建设投入过大ꎬ将强相对差值限设 置为 20%ꎮ 此外ꎬ按照涌水量预测的相对差值ꎬ进 一步将涌水量预测的准确程度按照准确至不准确或 偏离ꎬ以涌水量预测的主要功能为导向ꎬ可以细分经 济型Ⅰ、安全型Ⅱ和偏离型Ⅲ3 个等级ꎬ具 体见表 2ꎮ 事实上ꎬ相对差值可以转化为绝对差值的表达 方式ꎮ 参考上述标准ꎬ可以根据矿井水文地质条件 432 贺晓浪等矿井涌水量预测方法的改进及结果准确性判定2020 年第 8 期 的不同ꎬ尤其是涌水量的不同ꎬ采用涌水量绝对差值 进行准确等级的判定ꎬ如水文地质条件中等ꎬ涌水量 小于且接近 600 m3/ h 的矿井ꎬ其涌水量预测结果的 准确程度判定标准可如表 3ꎬ其他水文条件的矿井ꎬ 也可以参考表 2 进行预测水量准确性的判定ꎮ 表 2 以相对差值计算的涌水量预测准确性等级判定标准 Table 2 Criteria for determining the accuracy of water inflow prediction calculated by relative difference 序号相对差值/ %预测准确性等级 1-1010经济型Ⅰ 21020安全型Ⅱ 320偏离型Ⅲ 表 3 以绝对差值计算的涌水量预测准确性等级判定标准 Table 3 Criteria for determining accuracy of water in flow prediction calculated by absolute difference 序号绝对差值/ %预测准确性等级 1-6060经济型Ⅰ 260120安全型Ⅱ 3120偏离型Ⅲ 总体而言ꎬ采用绝对差值和相对差值判定这两 种方法各有优缺点表 4ꎬ可以将两者结合使用ꎮ 但需要注意的是ꎬ对于预测涌水量准确性的判 定ꎬ其主要针对的是涌水量中等及以上的情况ꎬ而对 于正常涌水量小于 50 m3/ hꎬ最大涌水量小于 100 m3/ h 的矿井ꎬ由于难以避免的要出现模型误差、观 测误差、截断误差、舍入误差等因素存在ꎬ因而实际 值与预测值必然有一定差值ꎮ 而当实际值较小时ꎬ 绝对差值小ꎬ但相对差值大ꎬ进而使得两者判定方法 得出完全不同的结论ꎬ此时必须结合矿井实际的排 水系统设置进行准确性判定ꎬ而非仅依靠某一种判 定方法ꎮ 表 4 涌水量准确性判定方法对比 Table 4 Comparison of methods for determining the accuracy of water inflow 判定方法缺点优点 绝对差值判定法无法展示涌水特征数据对比直观ꎬ清晰 相对差值判定法数据对比不够直观可以展示本身涌水特点 4 结 论 1提出了矿井涌水量的属性特征包含统计属 性、水文属性和生产属性ꎬ矿井涌水量可以定义为 矿井在某一生产阶段时ꎬ单位时间内井下各涌水点 涌水量之和ꎮ 2阐述了“大井法”计算矿井涌水量应当包含 静储量 Qj和导水裂隙带直接导通含水层的动态补 给量 Qd1和巨厚含水层未直接导通部分的动态补给 量 Qd2三部分ꎬ提出了利用“同频率放大法”改进水 文地质比拟法ꎬ得到预测涌水量的动态曲线ꎻ将时间 序列分析和水文地质比拟相结合提高预测步长周 期ꎬ增强时间序列分析的适用性ꎬ从而实现涌水量 的精细、动态预测ꎮ 3提出采用“绝对差值法”和“相对差值法”进
展开阅读全文

资源标签

最新标签

长按识别或保存二维码,关注学链未来公众号

copyright@ 2019-2020“矿业文库”网

矿业文库合伙人QQ群 30735420