数值法和解析法在矿井涌水量预测中的应用分析

第３９卷第４期　２０１７年　４月　能源与环保　Ｃｈｉｎａ　Ｅｎｅｒｇｙ　ａｎｄ　Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｐｒｏｔｅｃｔｉｏｎ　Ｖ０１．３９　Ｎｏ．４　Ａｐｒ．　２０１７　数值法和解析法在矿井涌水量预测中的应用分析　张莉丽　，张耀文　，武　强　，折书群　（１．防灾科技学院，河北三河０６５２０１；２．中国矿业大学（北京），北京１０００８３；　３．华北有色工程勘察院有限公司，河北石家庄０５００２４）　摘要：针对马城铁矿矿井水文地质条件、工程地质条件复杂及开采技术条件为复杂型复合问题矿床这　一特点，在分析矿井地层资料以及以往的水文地质资料、查明矿区水文地质条件基础上，采用数值法　和解析法对该矿井的涌水量进行了预测计算，所得结果可以作为矿井建井的设计依据。将２种方法　的预测结果进行分析比较，综合考虑这２种方法的优缺点和矿井实际情况，认为数值法计算结果比较　合理，研究结论可为其他类似矿井涌水量的预测提供借鉴。　关键词：水害；数值法；解析法；马城铁矿；矿井涌水量　中图分类号：ＴＤ７４２．１　文献标志码：Ａ　文章编号：１００３—０５０６（２０１７）０４—００５２—０５　Ａｐｐｌｉｃａｔｉｏｎ　ａｎａｌｙｓｉｓ　ｏｆ　ｎｕｍｅｒｉｃａｌ　ｍｅｔｈｏｄ　ａｎｄ　ａｎａｌｙｔｉｃａｌ　ｍｅｔｈｏｄ　ｉｎ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｍｉｎｅ　ｄｉｓｃｈａｒｇｅ　ｉｎ　ｍｉｎｅ　ｓｈａｆｔ　Ｚｈａｎｇ　Ｌｉｌｉ　，Ｚｈａｎｇ　Ｙａｏｗｅｎ　，Ｗｕ　Ｑｉａｎｇ　，Ｓｈｅ　Ｓｈｕｑｕｎ　（１．Ｉｎｓｔｉｔｕｔｅ　ｏｆＤｉｓａｓｔｅｒ　Ｐｒｅｖｅｎｔｉｏｎ，Ｓａｎｈｅ　０６５２０１，Ｃｈｉｎａ；　２．Ｃｈｉｎａ　Ｕｎｉｖｅｒｓｉｔｙ　ｏｆＭｉｎｉｎｇ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ，Ｂｅｉｊｉｎｇ　１０００８３，Ｃｈｉｎａ；　３．Ｎｏｒｔｈ　Ｃｈｉｎａ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｉｎｖｅｓｔｉｇａｔｉｏｎ　Ｉｎｓｔｉｔｕｔｅ　Ｃｏ．，Ｌｔｄ，Ｓｈｉｉｆａｚｈｕａｎｇ　０５００２４，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｔｈｅ　ｇｅｏｌｏｇｉｃａｌ　ｃｏｎｄｉｔｉｏｎｓ　ａｎｄ　ｅｎｇｉｎｅｅｒｉｎｇ　ｃｏｎｄｉｔｉｏｎｓ　ｏｆ　Ｍａｃｈｅｎｇ　ｉｒｏｎ　ｏｒｅ　ａｒｅ　ｃｏｍｐｌｅｘ，ａｎｄ　ｔｈｅ　ｍｉｎｉｎｇ　ｍｅｔｈｏｄ　ｕｓｅｄ　ｉｎ　ｔｈｅ　ｏｒｅ　ｉｓ　ｃｏｍｐｌｅｘ　ａｎｄ　ｃｏｍｂｉｎｅｄ．Ｂａｓｅｄ　ｏｎ　ｔｈｅ　ｗａｔｅｒ　ｉｎｆｌｏｗ　ｏｂｓｅｒｖａｔｉｏｎ　ｄａｔａ　ａｎｄ　ｔｈｅ　ｓｔｒａｔｕｍ　ｄａｔａ　ｏｆ　Ｍａｃｈｅｎｇ　ｉｒｏｎ　ｏｒｅ，ｔｈｅ　ｗａｔｅｒ　ｉｎｆｌｏｗ　ｗａｓ　ｃａｌ—　ｃｕｌａｔｅｄ　ｂｙ　ｔｈｅ　ｎｕｍｅｒｉｃａｌ　ｍｅｔｈｏｄ　ａｎｄ　ａｎａｌｙｔｉｃ　ｍｅｔｈｏｄ．Ｔｗｏ　ｐｒｅｄｉｃｔｅｄ　ｒｅｓｕｈｓ　ｗｅｒｅ　ｓｉｍｉｌａｒ　ａｎｄ　ｃｏｕｌｄ　ｂｅ　ｕｓｅｄ　ａｓ　ｔｈｅ　ｄｅｓｉｇｎ　ｂａｓｉｓ　ｆｏｒ　ｉｒｏｎ　ｏｒｅ　ｃｏｎｓｔｒｕｃｔｉｏｎ．Ｔｈｅ　ｃｏｒｒｅｓｐｏｎｄｉｎｇ　ａｐｐｌｉｃａｔｉｏｎ　ｃｏｎｄｉｔｉｏｎｓ　ｆｏｒ　ｔｈｅｓｅ　ｔｗｏ　ｍｅｔｈｏｄｓ　ｗｅｒｅ　ｇｉｖｅｎ　ｏｕｔ　ａｆｔｅｒ　ｔｈｅ　ａｎａｌｙｓｉｓ　ａｎｄ　ｔｈｅ　ｃｏｍｐａｒｉｓｏｎ，　ａｎｄ　ｉｔ　ｉｓ　ｃｏｎｓｉｄｅｒｅｄ　ｔｈａｔ　ｔｈｅ　ｒｅｓｕｈ　ｏｆ　ｎｕｍｅｒｉｃａｌ　ｃａｌｃｕｌａｔｉｏｎ　ｉｓ　ｒｅａｓｏｎａｂｌｅ．Ｉｔ　ｐｒｏｖｉｄｅｓ　ａ　ｒｅｆｅｒｅｎｃｅ　ｆｏｒ　ｏｔｈｅｒ　ｍｉｎｅｓ　ｔｏ　ｐｒｅｄｉｃｔ　ｔｈｅ　ｗａｔｅｒ　ｉｎ—　ｌｆｏｗ．　Ｋｅｙｗｏｒｄｓ：ｗａｔｅｒ　ｈａｚａｒｄ；ｌｍｍｅｒｉｃａｌ　ｍｅｔｈｏｄ；ａｎａｌｙｔｉｃａｌ　ｍｅｔｈｏｄ；Ｍａｃｈｅｎｇ　ｉｒｏｎ　ｏｒｅ；ｗａｔｅｒ　ｉｎｆｌｏｗ　０　引言　水害是影响和制约矿井安全生产的主要问题之　一型选择不合理等　。目前，预测矿井涌水量的方法　主要分为２类：①确定性的数学模型法，其中常用的　方法为数值法　、解析法　和水均衡法¨　；②统计　分析方法，如水文地质比拟法¨　、相关分析法¨　、　，如何准确地预测涌水量，一直都是矿井防治水工　作中一项复杂而困难的工作。尽管多年来，许多专　家学者们进行了大量的研究工作　。　，但预测涌水　量与实际开采涌水量相比，误差还是很大，有的误差　甚至达到数倍甚至数十倍。究其原因，主要是水文　涌水量降深曲线方程法　。。等。对于第①类方法来　说，如果水文地质条件没有查清楚，或者水文地质条　件比较复杂而概化不准确，都有可能导致计算结果　错误。而第②类方法则只能应用于地质条件相似的　矿井，而且无法考虑开采条件下影响矿井充水因素　地质参数不准确，或者水文地质条件不清楚，或者模　收稿日期：２０１６—１２—２１；责任编辑：秦爱新ＤＯＩ：１０．１９３８９／ｊ．ｃｎｋｉ．１００３—０５０６．２０１７．０４．Ｏｌ１　基金项目：中国地震局教师科研基金资助项目（２０１５０１０２）；国家自然科学基金资助项目（４１５７２２２２）　作者简介：张莉丽（１９７８一），女，内蒙古呼伦贝尔人，讲师，博士后，主要从事水文地质、工程地质和矿井防治水方面的研究。　引用格式：张莉丽，张耀文，武强，等．数值法和解析法在矿井涌水量预测中的应用分析［Ｊ］．能源与环保，２０１７，３９（４）：５２－５６．　Ｚｈａｎｇ　Ｌｉｌｉ，Ｚｈａｎｇ　Ｙａｏｗｅｎ，Ｗｕ　Ｑｉａｎｇ，ｅｔ　ａ１．Ａｐｐｌｉｃａｔｉｏｎ　ａｎａｌｙｓｉｓ　ｏｆ　ｎｕｍｅｒｉｃａｌ　ｍｅｔｈｏｄ　ａｎｄ　ａｎａｌｙｔｉｃａｌ　ｍｅｔｈｏｄ　ｉｎ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｍｉｎｅ　ｄｉｓｃｈａｒｇｅ　ｉｎ　ｍｉｎｅ　ｓｈａｆｔ［Ｊ＿．Ｃｈｉｎａ　Ｅｎｅｒｇｙ　ａｎｄ　Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｐｒｏｔｅｃｔｉｏｎ，２０１７，３９（４）：５２－５６．　・５２・　２０１７年第４期　张莉丽，等：数值法和解析法在矿井涌水量预测中的应用分析　第３９卷　的动态变化，不具有普遍性，大多情况下很难达到这　样的要求。马城铁矿周围没有与其地质条件相似的　矿井，不具备第②类方法预测矿坑涌水量的基础，因　此，根据矿区水文地质条件和矿山开采方式，采用数　刻的水头空间分布状况；Ｗ为源汇项，此模型中即　为矿坑涌水量；　为弹性给水度；ｑ为二类边界流　量；日为矿区四周无限边界；１２为模拟范围。　１．３数值模型　值法和解析法预测矿坑涌水量，两种方法相互印证，　力求获得精确的涌水量预测值　。　１数值法预测矿井涌水量　１．１水文地质模型概化　（１）含水系统。矿区含水层自上而下可大致划　分为第四系上部强孔隙含水层、下部中等孔隙含水　层、基岩风化裂隙弱含水层、基岩构造裂隙含水层　（带）。上述各含水层构成统一的地下含水系统，该　含水系统形似“蘑菇”。第四系含水层和基岩风化　裂隙含水层形似“蘑菇头”，是矿坑充水水源，基岩　构造裂隙含水层（带）形似“蘑菇茎”，是矿坑充水　通道。　（２）地下水流动系统。在天然状态下，地下水　大致由北向南运动，开采条件下，深部疏干排水，地　下水垂向越流补给，垂向上存在水头梯度，地下水呈　三维空间流场。　（３）边界条件。第四系含水层富水性极强，在　矿山开采过程中是矿坑充水水源，为定水头补给边　界，基岩风化带为无限补给边界。构造裂隙含水　（层）带两翼基岩裂隙不发育，透水性弱，为隔水边　界。由于９００　ｒｎ深以下基岩完整，透水性弱，视为计　算模型底边界，即隔水边界。为避免数值模型过小　造成较大计算误差，设计长１８　ｋｍ、宽１７　ｋｍ的计算　模型。故将含水系统概化为非均质各向异性含水系　统。　１．２数学模型　根据前述研究区的边界条件，含水介质特征以　及地下水的补给、径流、排泄条件，可以将研究区的　水文地质模型概化为三维非稳定流，其数学模型为：　Ｋｘ　（　）＋　０　（Ｏ　Ｈ）＋　去（　）＋Ｗ＝　ＯＨ（　，，　ｙ，ｚ）∈　∈　ｑ（　，Ｙ，　，ｔ）＝０（　，ｙ，　）∈Ｂ　（２）　Ｈ（ｐｃ，ｙ，ｚ，ｔ）＝Ｈｏ（ＰＣ，ｙ，　，０）　（　，ｙ，　）∈１２（３）　式中，　，Ｋ，　分别为　、Ｙ、　方向含水层的渗透系　数；日为任意空间点（　，ｙ，　）的水头随时间ｔ函数；　为任意空间点（　，），，ｚ）的初始水头函数，即ｔ。时　（１）渗流区剖分及初始渗流场确定。模型采用　有限差分法对物理模型进行矩形剖分，对地层分界　线、抽水孔、观测孔附近进行加密处理，共剖分成　１２０行，１０４列，２８层，共计３４９　４４０个单元格。根据　抽水试验前的统测水位资料，绘制地下水流场图，作　为模型识别的初始流场，初始渗流场如图１所示。　图１初始渗流场　Ｆｉｇ．１　Ｉｎｉｔｉａｌ　ｆｌｏｗ　ｆｉｅｌｄ　ｍａｐ　（２）水文地质参数确定。第四系含水层参数主　要依据马城水文地质勘探成果，第四系上部强含水　层渗透系数取值２３０　ｍ／ｄ，中部极弱透水层取值　０．００４　ｍ／ｄ，下部中等含水层取值３２　ｍ／ｄ。基岩含　水层参数主要依据抽水试验确定，并充分结合岩心　的ＲＱＤ值、线裂隙率、简易水文等资料给出不同标　高含水层水文地质参数初值；模型弹性给水度初值　的选取依据ＤＣＫ０２多孔抽水试验配线法确定，给出　初始值范围１０～～１０～。通过拟合调整参数，最终　确定基岩含水层在不同标高下的参数分区及参数值　见表１、表２。　（３）矿井涌水量预测。依据矿山开采设计方　案，向上每隔１２０　ｍ设计一个开采水平，即一９００　ｍ　与一５４０　ｍ进行联合开采、一７８０　ｍ与一４２０　ｍ进行　联合开采、一６６０　ｍ与一３００　１ＴＩ进行联合开采，一５７０　ｍ与一２４０　ｍ进行联合开采，２个开采高程采取上行　分层充填法，最上面一个开采水平至一２４０　ｍ标高，　一２４０　ｍ标高以上矿体及围岩留设矿坑安全顶板，　不再开采。结合开采设计方案，对一２４０　ｍ开采水　平进行涌水量预测，涌水量预测结果：正常涌水量为　６．１１万ｍ　／ｄ，最大涌水量为７．３３万ｍ　／ｄ。　．５３・　２０１７年第４期　能源表１水文地质模型渗透系数　Ｔａｂ．１　Ｐｅｒｍｅａｂｉｌｉｔｙ　ｃｏｅｆｉｆｃｉｅｎｔ　ｏｆ　ｔｈｅ　ｈｙｄｒｏｇｅｏｌｏｇｉｃａｌ　ｍｏｄｅｌ　参数分渗透系数（ｍ・ｄＩ１）　参数分渗透系数（ｍ・ｄ　）　区编号　区编号　Ｋｙ　基岩强　１　０．２５　０．２５　０．０２５　基岩弱　１　０．２５１　０．２５１　０．０２５　风化带　２　５ｅ一　５ｅ一３　５ｅ一４　风化带　２　６ｅ一　６ｅ一３　６ｅ一４　１　６．２５　６．２５　６．２ｌ　ｌ　１．１０　１．１０　ｌ　１０　２　２．３０　２．３Ｏ　２．３０　２　２．３０　２．３Ｏ　２．３０　３　０．３０　０．３０　０．３０　—５００　３　１．８０　１．８０　１．８０　—２５０　４　１．９５　１．９５　１．９３　４　１．９５　１．９５　ｌ　９３　ｉｎ　５　１．１０　１．１Ｏ　１．１０　５　３ｅ一４　３ｅ一４　３ｅ一４　６　１．８０　１．８Ｏ　１．８０　ｌ　１．８０　１．８０　１　８Ｏ　７　３ｅ一４　３ｅ一４　３ｅ一４　一６００　２　２．３０　２．３０　２．３Ｏ　ｌ　６．２５　６．２５　６．２１　ｎｌ　３　１．１Ｏ　１．１Ｏ　１．１Ｏ　２　２．３０　２．３０　２．３Ｏ　４　３ｅ一４　３ｅ一４　３ｅ一４　３００　３　１．９５　１．９５　１．９３　１　１．１Ｏ　１　１Ｏ　１．１Ｏ　ｍ　４　１．１０　１．１０　１．１Ｏ　一７００　２　２．３０　２．３０　２．３０　５　１．８Ｏ　１．８０　１．８０　３　３ｅ一４　３ｅ一４　３ｅ一４　６　３ｅ一　３ｅ－４　３ｅ　４　１　１．１Ｏ　１．１Ｏ　１．１０　ｌ　６．２５　６．２５　６．２１　—８００　２　２．３０　２．３０　２．３０　ｍ　２　２．３０　２．３０　２．３０　３　３ｅ一４　３ｅ一４　３ｅ一４　４００　３　１．９５　１．９５　１．９３　１　１．１０　１．１０　１．１０　ｉｎ　４　１．１Ｏ　１．１Ｏ　１．１Ｏ　９００　２　２．３０　２．３０　２．３０　５　１．８Ｏ　１．８０　１．８０　Ｉｎ　６　３ｅ一４　３ｅ一４　３ｅ一４　３　３ｅ一　３ｅ一４　３ｅ～４　表２水文地质模型给水度　Ｔａｂ．２　Ｓｐｅｃ￣ｃ　ｙｉｅｌｄ　ｏｆ　ｔｈｅ　ｈｙｄｒｏｇｅｏｌｏｇｉｃａｌ　ｍｏｄｅｌ　含水参数分　给水度　含水参数分　给水度　层区编号弹性给水度重力给水度　层区编号弹性给水度重力给水度　Ａ　６．９×１０－。１．０×１０—３　Ａ　６．９　ｘ１０＿。１．０×１０　３　一２５０　Ｂ　１４×１０—６　２０×１０—３　一６００　Ｂ　１ｍ．．．４×１０—６　２．０×１０—３　ｍＣ　３．１×１０—６　２　４×１０—３　Ｃ　３　ｌ×１０—６　２．４×１０—３　Ａ　６．９×１０－。１　０×１０—３　Ａ　６９×１０～１．０×１０—３　．一　Ｂ　１．４×１０～２．０×１０＿３　一　Ｂ　１．４×１０　２．０×１０＿３　Ｃ　３．１　ｘ１０—６　２４×１０—３　Ｃ　３１×１０—６　２　４×１０—３　．．Ａ　６．９×１０＿　１．０×１０—３　８０ｏ　Ａ　６．９×１０～１　０×１０—３　一　Ｂ　１．４×１０　２．０×１０＿］　Ｃ　３　１×１０—６　２　４×１０—３　Ｂ　１．４×１０—６　２．０×１０—３　Ａ　６　９×１０～１．０×１０—３　一５０Ｏ　Ｂ　ｌ　４９ｏｏ　Ａ　６．９×１０～１－０×１０—３　×１０—６　２．０×１０—３　ｍＣ　３．１×１０—６　２４×１０—３　Ｂ　１．４×１０—６　２．０×１０—３　．２解析法预测矿井涌水量　矿区第四系水量丰富，且接受滦河、新滦河渗漏　补给和大气降水的入渗补给，可视为定水头补给边　界。第四系底部强风化带分布广泛，厚度稳定，可视　为弱透水边界，控制着第四系水进入矿坑的“阀　门”。弱风化带为风化裂隙水的主要赋存部位，可　视为均质无限承压含水层。３线以南和２４线以北　弱风化带（Ｂ区）以下基岩完整，含水微弱，可视为　隔水边界。３＿２４线基岩构造破碎为矿坑充水的主　要通道（Ａ区）。在开采条件下，３＿２４线构造破碎　含水带压力释放，引起第四系水的垂向越流补给以　－５４・　与环保　第３９卷　及弱风化带裂隙水的侧向补给，其中，基岩弱风化带　的侧向补给量由两部分组成：①弱风化带裂隙水接　受北部裸露山区地下水侧向补给量；②第四系水通　过强风化带垂向越流补给弱风化带水量（图２）。　降　第四系定水头补给边界　、＼　灞　线　２４线以北＼　３—２４线　，　／３线以南　Ｉ　』＼　Ｉ＋　　＋Ｉ　　＋　＋ｌ　／　、　、　、　强风化带弱透水层　弱ｊ　相对隔水层　Ｑｚ　Ｗ　Ｑ－弱　相瓦对　隔ｌ水　层含　在　图２—２４０　ｉｌｌ标高矿坑涌水示意　Ｆｉｇ．２　Ｓｋｅｔｃｈ　ｍａｐ　ｏｆ　ｔｈｅ　ｗａｔｅｒ　ｇｕｓｈｉｎｇ　ａｔ一２４Ｏ　ｍ　由于矿床开采自下而上进行，开采至最上一个　水平时，矿坑顶板最薄，第四系水和风化裂隙水流入　矿坑所通过的“瓶颈”最大，地下水垂向运动过程中　受的阻力最小，矿坑涌水量最大，因此，本次解析法　对一２４０　ｍ标高矿坑涌水量进行估算。开采至标高　一２４０　ｉｎ时，基岩地下水位降至强风化带底板以下，　Ａ区构造破碎带成为无压区，接受第四系越流量不　再随降落漏斗向外扩展而增加，Ｑ　为定值。Ｂ区弱　风化带裂隙水Ｑ　随着漏斗向外扩展而增大（主要　第四系水越流量增大），对Ｑ　计算可采用存在垂向　越流补给的承压转无压“大井法”进行计算。各部　分水量计算见式（４）：　Ｑ总＝Ｑ１＋Ｑ２　Ｑ　＝Ｋ　１　１　Ｆ　㈡　（４）　式中，Ｑ　为２４０　ｍ标高矿坑涌水量；Ｐ，为Ａ区第四　系越流量；Ｑ　为Ｂ区弱风化带侧向补给量（内含第　四系越流量）；　。为强风化带垂向渗透系数；　为曰　区弱风化带渗透系数；Ｍ　为强风化带平均厚度；Ｍ：　为弱风化带平均厚度；ｎｏ为第四系水位；Ｈ　为　区　强风化带底板水位；Ｓ为大井法预测降深值；Ｆ为Ａ　区接收第四系水越流补给面积；ｒｎ为大井引用半径；　／ｒ　、　Ｂ为越流因素；Ｋｏｆ　ｌ为零阶第二类虚宗量Ｂｅｓｓｅｌ　２０１７年第４期　张莉丽，等：数值法和解析法在矿井涌水量预测中的应用分析　第３９卷　函数；　ｆ罟Ｉ＼Ｄ／　为一阶第二类虚宗量Ｂｅｓｓｅｌ函数。　经计算，一２４０　ｍ开采中段矿坑正常涌水量为　６．３３万ｍ　／ｄ，最大涌水量为７．６０万ｍ　／ｄ。对于正　常涌水量来说，　区第四系越流量占５．２４万ｍ。／ｄ，　占８３％，Ｂ区弱风化带侧向补给量占１．０９万ｍ　／ｄ，　占１７％，可见第四系越流量为矿坑涌水的主要来　源。　３结果对比　数值法模型介质场构建严格按照矿区钻孔水文　地质编录进行，重点剖分了构造破碎带、水文地质勘　探区域。在模型的识别过程中，采取了群孔抽水试　验中不同深度、不同地层的１３个观测孔３０９　ｈ，共计　１　５４７个观测数据作为模型识别验证目标，最终获得　了能够代表马城铁矿水文地质特点的水文地质数值　模型，据此预测了矿井最大涌水量和正常涌水量。　预测基岩含水层降落漏斗扩展范围小于４．５　ｋｍ　Ｘ　５．０　ｋｍ，小于模型范围，该计算结果可靠，能够反映　马城铁矿坑下开采阶段涌水量特点。　解析法估算了一２４０　ｍ标高的矿坑涌水量，与　数值法预测结果相近。由于解析法各参数的取值代　表性较差，假设条件与实际有一定偏差，矿坑涌水量　预测以数值法为准，以解析法为参考。数值法预测　矿坑涌水量精度Ｂ级，可信度０．５；解析法估算矿坑　涌水量精度ｃ级，可信度０．４ｌ２¨。因此，数值法相　对于解析法的结果准确可靠（表３）。　表３预测结果对比　Ｔａｂ・３　Ｃｏｍｐａｒｉｓｏｎ　ｏｆ　ｐｒｅｄｉｃｔｉｏｎ　ｒｅｓｕｌｔｓ　算法　正常涌水量／（万ｍ　・ｄ　）　最大涌水量／（万ｍ’・ｄ　）　数值法　６．１１　７．３３　壁堑鎏　：　！：　数值法预测矿坑涌水量是建立在矿山上向分层　充填开采法的基础上，即不得破坏矿床上覆第四系、　风化裂隙含水层组介质场结构。该区新构造运动频　繁，如果由于新构造运动或采矿活动造成矿床上覆　顶板坍塌，第四系上部强含水层组中的孔隙水将大　量涌入矿坑，形成矿山突涌水及矿坑泥石流灾害，后　果不堪设想。　４　结论　（１）基于矿井实际情况，建立了符合矿区水文　地质特点的水文地质数值模型来预测矿井涌水量，　可以处理复杂边界条件，结果更切合实际，所以结果　准确率较高。　（２）解析法对水文地质条件要求比较苛刻，适　用于均质、形状规则的含水层，而实际含水层不会这　样理想，这就使解析法的实际应用受到了一定的限　制，预测结果可信度不是很高，只能作为辅助的预测　方法。　参考文献（Ｒｅｆｅｒｅｎｃｅｓ）：　［１］　杜敏铭，邓英尔，许模．矿井涌水量预测方法综述［Ｊ］．地质学　报，２００９，２９（１）：７０－７３．　Ｄｕ　Ｍｉｎｍｉｎｇ，Ｄｅｎｇ　Ｙｉｎｇ　ｅｒ，Ｘｕ　Ｍｕｏ．Ｒｅｖｉｅｗ　ｏｆ　ｍｅｔｈｏｄｏｌｏｇｙ　ｆｏｒ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｗａｔｅｒ　ｗｉｅｌｄ　ｏｆ　ｍｉｎｅ［Ｊ］．Ａｃｔａ　Ｇｅｏｌｏｇｉｃａ　Ｓｉｎｉｃａ，２００９，　２９（１）：７０－７３．　…　武强，崔芳鹏，赵苏启，等．矿井水害类型划分及主要特征分析　［Ｊ］．煤炭学报，２０１３，３８（４）：５６１—５６５．　Ｗｕ　Ｑｉａｎｇ，Ｃｕｉ　Ｆａｎｇｐｅｎｇ，Ｚｈａｏ　Ｓｕｑｉ，ｅｔ　ａ１．Ａｎａｌｙｓｉｓ　ｏｎ　ｃｌａｓｓｉｆｉｃａ－　ｔｉｏｎ　ａｎｄ　ｍａｉｎ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｍｉｎｅ　ｗａｔｅｒ　ｈａｚａｒｄ　ｆ　Ｊ］．Ｊｏｕｍａｌ　ｏｆ　．　Ｃｈｉｎａ　Ｃｏａｌ　Ｓｏｃｉｅｔｙ，２０１３，３８（４）：５６１－５６５．　：　武强．煤矿防治水手册［Ｍ］．北京：煤炭工业出版社，２０１３．　武强．我国矿井水防控与资源化利用的研究进展、问题和展望　［Ｊ］．煤炭学报，２０１４，３９（５）：７９５－８０５．　Ｗｕ　Ｑｉａｎｇ．Ｐｒｏｇｒｅｓｓ，ｐｒｏｂｌｅｍｓ　ａｎｄ　ｐｒｏｓｐｅｃｔｓ　ｏｆ　ｐｒｅｖｅｎｔｉｏｎ　ａｎｄ　ｃｏｎ—　ｔｒｏｌ　ｔｅｃｈｎｏｌｏｇｙ　ｏｆ　ｍｉｎｅ　ｗａｔｅｒ　ｎａｄ　ｒｅｕｔｉｌｉｚａｔｉｏｎ　ｉｎ　Ｃｈｉｎａ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｃｈｉｎａ　Ｃｏａｌ　Ｓｏｃｉｅｔｙ，２０１４，３９（５）：７９５－８０５．　武强，赵苏启，孙文洁，等．中国煤矿水文地质类型划分与特征　分析［Ｊ］．煤炭学报，２０１３，３８（６）：９０１－９０５．　Ｗｕ　Ｑｉａｎｇ，Ｚｈａｏ　Ｓｕｑｉ，Ｓｕｎ　Ｗｅｎｊｉｅ，ｅｔ　ａ１．Ｃｌａｓｓｉｉｆｃａｔｉｏｎ　ｏｆ　ｔｈｅ　ｈｙｄｒｏｇｅｏｌｏｇｉｅａｌ　ｔｙｐｅ　ｏｆ　ｃｏａｌ　ｍｉｎｅ　ａｎｄ　ａｎａｌｙｓｉｓ　ｏｆ　ｉｔｓ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｉｎ　Ｃｈｉｎａ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｃｈｉｎａ　Ｃｏａｌ　Ｓｏｃｉｅｔｙ，２０１３，３８（６）：９０１—　９０５．　［６］　杨永国．灰色控制系统及其在矿井涌水量预测中的应用［Ｊ］．　中州煤炭，１９８９（５）：３０—３２．　Ｙａｎｇ　Ｙｏｎｇｇｕｏ．Ｇｒｅｙ　ｃｏｎｔｒｏｌ　ｓｙｓｔｅｍ　ａｎｄ　ｉｔｓ　ａｐｐｌｉｃａｔｉｏｎ　ｉｎ　ｐｒｅｄｉｃ—　ｔｉｏｎ　ｏｆ　ｍｉｎｅ　ｗａｔｅｒ　ｉｎｆｌｏｗ［Ｊ］．Ｚｈｏｎｇｚｈｏｕ　Ｃｏａｌ，１９８９（５）：３０－３２．　［７］　王义明，娄源改，杨新胜，等．注浆井筒的质量检验方法及其涌　水量计算［Ｊ］．中州煤炭，２００５（２）：３４．　Ｗａｎｇ　Ｙｉｍｉｎｇ，Ｌｏｕ　Ｙｕａｎｇａｉ，Ｙａｎｇ　Ｘｉｎｓｈｅｎｇ，ｅｔ　ａ１．Ｑｕａｌｉｔｙ　ｉｎｓｐｅｃ—　ｔｉｏｎ　ｍｅｔｈｏｄ　ｏｆ　ｇｒｏｕｔｉｎｇ　ｗｅｌｌｂｏｒｅ　ａｎｄ　ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ｗａｔｅｒ　ｉｎｆｌｏｗ　［Ｊ］．Ｚｈｏｎｇｚｈｏｕ　Ｃｏａｌ，２００５（２）：３－４．　『８］　祝金峰，宁超．正村井田地下水赋存特征分析及矿井涌水量计　算［Ｊ］．中州煤炭，２００８（１）：２３＿２４．　Ｚｈｕ　Ｊｉｎｆｅｎｇ，Ｎｉｎｇ　Ｃｈａｏ．Ｍｉｎｅ　ｈｙｄｒｏｇｅｏｌｏｇｙ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ａｎａｌｙｓｉｓ　ａｎｄ　ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ｗａｔｅｒ　ｙｉｅｌｄ　ｏｆ　ｍｉｎｅ　ｆｒｏｍ　Ｚｈｅｎｇｃｕｎ　ｍｉｎｅ　ｆｉｅｌｄ　［Ｊ］．Ｚｈｏｎｇｚｈｏｕ　Ｃｏａｌ，２００８（１）：２３－２４．　［９］　武强，许珂，张维．再论煤层顶板涌（突）水危险性预测评价的　“三图一双预测法”［Ｊ］．煤炭学报，２０１６，４１（６）：１—１３４７．　Ｗｕ　Ｑｉａｎｇ，Ｘｕ　Ｋｅ，Ｚｈａｎｇ　Ｗｅｉ．Ｆｕｒｔｈｅｒ　ｒｅｓｅａｒｃｈ　ｏｎ”ｔｈｒｅｅ　ｍａｐｓ—ｔｗｏ　・５５・　２０１７年第４期　能源与环保　第３９卷　ｐｒｅｄｉｃｔｉｏｎｓ”ｍｅｔｈｏｄ　ｆｏｒ　ｐｒｅｄｉｃｔｉｏｎ　ｏｎ　ｃｏａｌ　ｓｅａｍ　ｒｏｏｆ　ｗａｔｅｒ　ｂｕｒｓｔｉｎｇ　［１６］袁广祥，杨莉，李刚．解析法预测隧道可能最大涌水量的对比　研究［Ｊ］．华北水利水电学院学报，２０１３，３４（６）：１２．１５．　Ｙｕａｎ　Ｇｕａｎｇｘｉａｎｇ，Ｙａｎｇ　Ｌｉ，Ｌｉ　Ｇａｎｇ．Ｃｏｎｔｒａｓｔｉｖｅ　ｒｅｓｅａｒｃｈ　ｏｎ　ａｎａ－　ｌｙｒｉｃａｌ　ｍｅｔｈｏｄｓ　ｔｏ　ｐｒｅｄｉｃｔ　ｐｒｏｂａｂｌｅ　ｍａｘｉｍｕｍ　ｗａｔｅｒ　ｙｉｅｌｄ　ｉｎｔｏ　ｔｕｎ—　ｒｉｓｋ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｃｈｉｎａ　Ｃｏａｌ　Ｓｏｃｉｅｔｙ，２０１６，４１（６）：１３４１—１３４７．　［１０］武强，李博．煤层底板突水变权评价中变权区间及调权参数确　定方法［Ｊ］．煤炭学报，２０１６，４１（９）：２１４３－２１４９．　Ｗｕ　Ｑｉａｎｇ，Ｌｉ　Ｂｏ．Ｄｅｔｅｒｍｉｎａｔｉｏｎ　ｏｆ　ｖａｒｉａｂｌｅ　ｗｅｉｇｈｔ　ｉｎｔｅｒｖａｌ　ａｎｄ　ａｄ—　ｊｕｓｔ　ｗｅｉｇｈｔ　ｐａｒａｍｅｔｅｒｓ　ｉｎ　ｔｈｅ　ｖａｒｉａｂｌｅ　ｗｅｉｇｈｔ　ａｓｓｅｓｓｍｅｎｔ　ｍｏｄｅｌ　ｏｆ　ｎｅｌ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　ＮｏＡｈ　Ｃｈｉｎａ　Ｉｎｓｔｉｔｕｔｅ　ｏｆ　Ｗａｔｅｒ　Ｃｏｎｓｅｒｖａｎｃｙ　ａｎｄ　Ｈｙｄｒｏｅｌｅｃｔｉｒｃ　Ｐｏｗｅｒ，２０１３，３４（６）：１２—１５．　ｗａｔｅｒ　ｉｎｒｕｓｈ　ｆｒｏｍ　ｃｏａｔ　ｆｌｏｏｒ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｃｈｉｎａ　Ｃｏａｌ　Ｓｏｃｉｅｔｙ，　２０１６，４１（９）：２１４３・２１４９．　［１７］董兴文．矿井最大涌水量计算和参数确定——试论南方岩溶　矿井的水均衡预测方法［Ｊ］．四川地质学报，１９９５，１５（１）：７５　８０．　［１１］吴燕军，申建军，赵颖旺．工作面不同采宽与涌水量关系研究　［Ｊ］．中州煤炭，２０１６（６）：１８－２１．　Ｗｕ　Ｙａｈｊａｎ，Ｓｈｅｎ　Ｊｉａｎｊｕｎ，Ｚｈａｏ　Ｙｉｎｇｗａｎｇ．Ｒｅｓｅａｒｃｈ　ｏｎ　ｒｅｌａｔｉｏｎ－　Ｄｏｎｇ　Ｘｉｎｇｗｅｎ．Ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ｍａｘｉｍｕｍ　ｇｕｓｈｉｎｇ　ｗａｔｅｒ　ｉｎ　ｍｉｎｅｓ　ａｎｄ　ｉｔｓ　ｐａｒａｍｅｔｅｒ　ｄｅｔｅｒｍｉｎａｔｉｏｎ　－ｏｎ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｇｒｏｕｎｄｗａｔｅｒ　ｂａｌ－－　ｓｈｉｐ　ｂｅｔｗｅｅｎ　ｄｉｆｆｅｒｅｎｔ　ｗｉｄｔｈ　ｏｆ　ｗｏｒｋｉｎｇ　ｆａｃｅ　ａｎｄ　ｗａｔｅｒ　ｉｎｆｌｏｗ［Ｊ］．　Ｚｈｏｎｇｚｈｏｕ　Ｃｏａｌ，２０１６（６）：１８－２１．　ａｎｃｅ　ｏｆ　ｋａｒｓｔ　ｍｉｎｅｓ　ｉｎ　Ｃｈｉｎａ［Ｊ］．Ａｃｔａ　Ｇｅｏｌｏｇｉｃａ　Ｓｉｃｈｕａｎ，１９９５，　１５（１）：７５－８０．　［１２］翁荔玉，年宾．青东矿充水因素分析及涌水量预计［Ｊ］．中州煤　炭，２０１５（２）：１１３－１１５．　Ｗｅｎｇ　Ｌｉｙｕ，Ｎｉａｎ　Ｂｉｎ．Ａｎａｌｙｓｉｓ　ｏｆ　ｗａｔｅｒ　ｆｉｌｌｉｎｇ　ｆａｃｔｏｒｓ　ａｎｄ　ｆｏｒｅｃａｓｔ　［１８］李小勇，林坚，邱凤．水文地质比拟法在矿坑涌水量预测中的　应用［Ｊ］．资源环境与工程，２０１４，２８（１）：６６－６８．　Ｌｉ　Ｘｉａｏｙｏｎｇ，Ｌｉｎ　Ｊｉａｎ，Ｑｉｕ　Ｆｅｎｇ．Ｔｈｅ　ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ｈｙｄｒｏｌｏｇｉｃａｌ　ｏｆ　ｗａｔｅｒ　ｉｎｆｌｏｗ　ｉｎ　ｔｈｅ　Ｑｉｎｇｄｏｎｇ　Ｍｉｎｅ［Ｊ］．Ｚｈｏｎｇｚｈｏｕ　Ｃｏａｌ，２０１５　（２）：１１３－ｌ１５．　ｇｅｏｌｏｇｉｃａｌ　ａｎａｌｏｇｙ　ｍｅｔｈｏｄ　ｉｎ　ｔｈｅ　ｍｉｎｅ　ｗａｔｅｒ　ｉｎｆｌｏｗ　ｆｏｒｅｃａｓｔ［Ｊ］．　Ｒｅｓｏｕｒｃｅｓ　Ｅｎｖｉｒｏｎｍｅｎｔ　ａｎｄ　Ｅｎｇｉｎｅｅｒｉｎｇ，２０１４，２８（１）：６６－６８．　［１３］李双军，万小强，王照亮．新庄矿主井井筒涌水量预算［Ｊ］．中　州煤炭，２０１４（８）：３７－４０．　Ｌｉ　Ｓｈｕａｎｇｊｕｎ，Ｗａｎ　Ｘｉａｏｑｉａｎｇ，Ｗａｎｇ　Ｚｈａｏｌｉａｎｇ．Ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｍａｉｎ　［１９］段俭君，徐会军，王子河．相关分析法在矿井涌水量预测中的　应用［Ｊ］．煤炭科学技术，２０１３，４１（６）：１１３．１１６．　Ｄｕａｎ　Ｊｉａｎｊｕｎ，Ｘｕ　Ｈｕｉｊｕｎ，Ｗａｎｇ　Ｚｉｈｅ．Ｃｏｒｒｅｌａｔｉｏｎａｌ　ａｎａｌｙｓｉｓ　ｍｅｔｈ－　ｓｈａｆｔ　ｗａｔｅｒ　ｉｎｆｌｏｗ　ｉｎ　Ｘｉｎｚｈｕａｎｇ　Ｃｏａｌ　Ｍｉｎｅ［Ｊ］．Ｚｈｏｎｇｚｈｏｕ　Ｃｏａｌ，　２０１４（８）：３７－４０．　ｏｄ　ａｐｐｌｉｅｄ　ｔｏ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｍｉｎｅ　ｗａｔｅｒ　ｉｎｆｌｏｗ　ｑｕａｎｔｉｔｙ［Ｊ］．Ｃｏａｌ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ，２０１３，４１（６）：１１３—１１６．　［１４］陈酩知，刘树才，杨国勇．矿井涌水量预测方法的发展［Ｊ］．工　程地球物理学报，２００９，６（１）：６８－７２．　Ｃｈｅｎ　Ｍｉｎｇｚｈｉ，Ｌｉｕ　Ｓｈｕｃａｉ，Ｙａｎｇ　Ｇｕｏｙｏｎｇ．Ｔｈｅ　ｄｅｖｅｌｏｐｍｅｎｔ　ｏｆ　［２０］潘志．矿井涌水量的数学模型与预测［Ｊ］．宁夏大学学报（自然　科学版），１９９８，１９（３）：２２９－２３２．　Ｐａｎ　Ｚｈｉ．Ｍａｔｈｅｍａｔｉｃａ　ｍｏｄｅｌ　ｏｆ　ｍｉｎｅ　ｇｕｓｈｅｄ　ｗａｔｅｒ　ｑｕａｎｔｉｔｙ［Ｊ］．　Ｊｏｕｒｎａｌ　ｏｆ　Ｎｉｎｇｘｉａ　Ｕｎｉｖｅｒｓｉｔｙ（Ｎａｔｕｒａｌ　Ｓｃｉｅｎｃｅ　Ｅｄｉｔｉｏｎ），１９９８，１９　（３）：２２９－２３２．　ｍｉｎｉｎｇ　ｗａｔｅｒ　ｉｎｆｌｏｗ　ｐｒｅｄｉｃｔ　ｍｅｔｈｏｄ［Ｊ］．Ｃｈｉｎｅｓｅ　Ｊｏｕｒｎａｌ　ｏｆ　Ｅｎｇｉ—　ｎｅｅｒｉｎｇ　Ｇｅｏｐｈｙｓｉｃｓ，２００９，６（１）：６８－７２．　［１５］李平，郭会荣，吴孔军，等．王河煤矿矿井涌水量数值模拟及预　［２１］钱学溥，雷小乔，许超，等．矿坑涌水量的计算及其精度级别和　可信度［Ｊ］．中国煤炭地质，２０１６，２８（５）：３９－４１．　Ｑｉａｎ　Ｘｕｅｐｕ，Ｌｅｉ　Ｘｉａｏｑｉａｏ，Ｘｕ　Ｃｈａｏ，ｅｔ　ａ１．Ｍｉｎｅ　ｗａｔｅｒ　ｉｎｆｌｏｗ　ｅｓｔｉ・　测［Ｊ］．地球科学一中国地质大学学报，２０１１，３６（４）：７５５－７６０．　Ｌｉ　Ｐｉｎｇ，Ｇｕｏ　Ｈｕｉｒｏｎｇ，Ｗｕ　Ｋｏｎｇｊｕｎ，ｅｔ　ａ１．Ｎｕｍｅｒｉｃａｌ　ｓｉｍｕｌａｔｉｏｎ　ａｎｄ　ｆｏｒｅｃａｓｔ　ｏｆ　ｍｉｎｅ　ｄｉｓｃｈａｒｇｅ　ｉｎ　Ｗａｎｇｈｅ　Ｃｏａｌ　Ｍｉｎｅ［Ｊ］．Ｅａｒｔｈ　Ｓｃｉｅｎｃｅ－Ｊｏｕｒｎａｌ　ｏｆ　Ｃｈｉｎａ　Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　Ｇｅｏｓｃｉｅｎｃｅｓ，２０１１，３６（４）：　７５Ｓ．７６０　ｍａｔｉｏｎ，ａｃｃｕｒａｃｙ　ｌｅｖｅｌ　ａｎｄ　ｒｅｌｉａｂｉｌｉｔｙ［Ｊ］．Ｃｏａｌ　Ｇｅｏｌｏｇｙ　ｏｆ　Ｃｈｉｎａ，　２０１６，２８（５）：３９－４１．　（上接第５ｌ页）　［１６］陈美娅，谢宇琦，吴玉珍，等．牛粪、矿渣与污泥烧制高强陶粒　研究［Ｊ］．硅酸盐通报，２０１５，３４（１２）：３４４６－３４５０．　Ｃｈｅｎ　Ｍｅｉｙａ，Ｘｉｅ　Ｙｕｑｉ，Ｗｕ　Ｙｕｚｈｅｎ，ｅｔ　ａ１．Ｐｒｅｐａｒａｔｉｏｎ　ｏｆ　ｈｉｇｈ—　Ｗｅｎｇ　Ｈｕａｎｘｉｎ，Ｚｈａｎｇ　Ｊｉｎｊｕｎ，Ｃａｏ　Ｙａｎｓｈｅｎｇ，ｅｔ　ａ１．Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ａｎｄ　ｓｉｎｔｅｒｉｎｇ　ｔｅｃｈｎｏｌｏｇｙ　ｏｆ　ｈａｙｄｉｔｅ　ｍａｄｅ　ｏｆ　ｓｅｗａｇｅ　ｓｌｕｄｇｅ［Ｊ］．　Ｊｏｕｒｎａｌ　ｏｆ　Ｚｈｅｊｉａｎｇ　Ｕｎｉｖｅｒｓｉｔｙ（Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ），２０１　１，４５　（１０）：１８７７—１８８３．　ｓｔｒｅｎｇｔｈ　ｃｅｒａｍｓｉｔｅ　ｂｙ　ｆｉｉｎｇ　ｗｉｒｔｈ　ｃｏｗ　ｄｕｎｇ，ｓｌａｇ　ａｎｄ　ｓｌｕｄｇｅ　ａｓ　ｒａｗ　［１９ｌ　彭建军，刘恒波，宋美，等．赤泥高强陶粒的研制［Ｊ］．砖瓦，　２０１１（１０）：９－ｌ１．　ｍａｔｅｒｉａｌ［Ｊ］．Ｂｕｌｌｅｔｉｎ　ｏｆ　ｔｈｅ　Ｃｈｉｎｅｓｅ　Ｃｅｒａｍｉｃ　Ｓｏｃｉｅｔｙ，２０１５，３４　（１２）：３４４６＿３４５０．　Ｐｅｎｇ　Ｊｉｎｊｕｎ，Ｌｉａｕ　Ｈｅｎｇｂｏ，Ｓｏｎｇ　Ｍｅｉ，ｅｔ　ａ１．Ｒｅｓｅａｒｃｈ　ｏｎ　ｈｉｇｈ－　［１７］王佳福，吕剑明．利用城市污泥制备陶粒的研究［Ｊ］．硅酸盐通　报，２０１２，３１（３）：７０６－７１０．　Ｗａｎｇ　Ｊｉｆｕ．Ｌ（ａｉ　Ｊｉａｎｍｉｎｇ．Ｒｅｓｅａｒｃｈ　ｏｎ　ｐｒｅｐａｒａｔｉｏｎ　ｏｆ　ｃｅｒａｍｓｉｔｅ　ｎ—　ｓｔｒｅｎｇｔｈ　ｒｅｄ　ｍｕｄ　ｃｅｒａｍｉｓｉｔｅ［Ｊ］．Ｂｒｉｃｋ　ａｎｄ　Ｔｉｌｅ，２０１１（１０）：９－１１．　［２０ｌ　贺君，王启山，任爱玲．给水厂与污水厂污泥制陶粒技术研究　Ｊ［Ｊ］．环境工程学报，２００９，３（９）：１６５３—１６５７．　Ｈｅ　Ｊｕｎ，Ｗａｎｇ　Ｑｉｓｈａｎ，Ｒｅｎ　Ａｉｌｉｎｇ．Ｔｅｃｈｎｏｌｏｇｙ　ｒｅｓｅａｒｃｈ　ｏｎ　ｗａｔｅｒ　ｓｉｎｇ　ｃｉｔｙ　ｓｌｕｄｇｅ［Ｊ］．Ｂｕｌｌｅｔｉｎ　ｏｆ　ｔｈｅ　Ｃｈｉｎｅｓｅ　Ｃｅｒａｍｉｃ　Ｓｏｃｉｅｔｙ，　２０１２，３１（３）：７０６－７１０．　ｗｏｒｋｓ　ｓｌｕｄｇｅ　ａｎｄ　ｓｅｗａｇｅ　ｓｌｕｄｇｅ　ｆｏｒ　ｃｅｒａｍｓｉｔｅ［Ｊ］．Ｃｈｉｎｅｓｅ　Ｊｏｕｍａｌ　ｏｆ　Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，２００９，３（９）：１６５３—１６５７．　［１８］翁焕新，章金骏，曹彦圣，等．污泥陶粒的性能特征与烧制工艺　［Ｊ］．浙江大学学报（工学版），２０１１，４５（１０）：１８７７　１８８３．　・５６－