刘诗鹏,邱 天,吴志超,杨 新
(1.山东黄金地质矿产勘查有限公司,山东 莱州 261400;2.中国地质大学(武汉)环境学院,湖北 武汉 430000)
矿产资源是珍贵的自然资源,对于现阶段我国社会发展的进步至关重要。近年来,随着相关行业对矿产资源的需求的增加,矿井逐年加深,导致矿井热害问题愈发突出。本文结合焦家金成矿带矿区实际情况,模拟矿区地温场,预测可能会产生热害的矿井开采深度,并且分析人为抽水对地温场的影响。
胶东地区地质历史较为复杂,受早白垩纪岩浆热隆影响,侏罗纪形成的玲珑岩基强烈隆升,其上覆的变质岩层形成拆离断层,胶东地区三山岛—仓上断裂、新城—焦家断裂和招远—平度断裂均是铲形下滑断裂,属拆离断层的组成部分[1]。本文研究对象是新城—焦家断裂带,是胶东地区三大金矿成矿带之一。本文的地温场数值模拟中包含达西渗流场,依据水文地质边界条件,在研究区划分出模拟区作为三维数值模拟对象(图1)。
图1 研究区范围及测温孔位置图
收集及汇总焦家金成矿带矿区地温资料,共整理出15个测温孔数据,来源于矿区各个时期的勘探资料,测温孔位置如图1所示。
地温梯度受地层结构及热物理性质的影响比较明显,所以它能较好地反映地层的温度变化情况。部分测温孔ZK509和ZK652地温梯度为1.85℃/hm和2.4℃/hm,统计其他测温孔地温梯度值,均在1.6~3.0℃/hm的正常分布范围内,因此模拟区中出现地热异常区概率较低。
为进一步研究矿区地温场,分析水平和垂直方向的地温分布特征。根据测温孔数据,分别绘制了-400m、-800m标高水平温度等值线图,两个标高水平地温分布特征相似,-800m标高温度等值线图如下所示。
图2 模拟区-800m标高温度等值线图
等值线图上最大和最小温度值之间都是相差10℃左右,南北方向上看,地温西低东高;东西方向上看,地温等值线北疏南密,总体上从西到东,温度逐渐升高。
垂直方向上,地温与深度之间呈一元线性相关。设一元线性方程拟合数据,钻孔ZK509拟合公式为y=-0.0183x+8,R²值为0.996。
统计其他测温孔数据拟合的R²值,均大于0.95,表明拟合出的一元回归方程可信度非常高,证明地温在垂向上是随深度均匀增大的。
岩石在地球与大气层的热交换中起着介质的作用,是影响地球内部热量向外传递或地球内部热交换的重要因素[2]。根据钻孔揭露,研究区地层岩性中花岗岩和变辉长岩所占范围及比例都较大,而这两种岩性热导率也较大,所以在研究区范围内,因岩性而导致的地温热异常概率较低。
岩浆活动对地温场会产生影响,评定主要从两方面,一个是时间,一个岩浆体规模[3]。研究区内岩浆岩广布,但形成时期都比较早,对现今研究区内地温几乎没影响不大。
地下水是最活跃的地质因素,在地壳浅部广泛分布,热容量大易于流动,是围岩温度的重要影响因素[4]。研究区补给和径流条件较差,主要判断人为抽水活动是否对地温场有影响。
本次模拟范围为图1中模拟区,模拟区三条边界定为西北的莱州海岸线,东北的朱桥河,东南的地下分水岭。模型共分为4个地层,上覆岩层、断层上盘碎裂岩、断层下盘碎裂岩和基底岩层,再设置薄壁断层建立三维地质模型。为模拟矿区抽水条件,采用大井法概化各个矿区抽水井为一个主抽水井,并将其设置在模拟区中间位置,通过控制该井水位,来模拟井群抽水产生的水位下降效果。
本文设置有三种模拟模型:①初始模型,模型中的地下水位为自然水位;②验证模型,用于验证模型模拟精度模型;③预测模型:模拟抽水井水位为-1400m时情况。
达西渗流场模拟中共有三类边界,模拟采用定水头边界和隔水边界。模拟区朱桥河和地下分水岭定为隔水边界,莱州海岸线为定水头边界,水头为0.05m。
地温场模拟中有三种类型的边界条件,模拟采用给定温度值和给定热流值两种边界。模拟区上边界为给定温度值边界,下边界为给定热流值边界。
表1 各岩层主要参数
验证模型主要是验证模拟结果精度,测温孔数据是矿区开采过程中测出的,因此,通过条件概化,将模型地下水位设为-400m。该模型模拟结果作为与实际温度数据的对比数据。
在标高-100m位置,孔JFZK1和孔FJZK1实测值为14.81℃和15.2℃,模拟值为15.38℃和15.38℃,分别相差0.57℃和0.18℃,相对误差在3.85%和1.18%。统计不同深度模拟值实际值误差得出,两个测温孔温度值的平均相对误差为4.92%,相关系数均在0.99以上。综上所述,本次模拟结果产生的数据是可靠的。
初始模型为参照模型,模拟结果如图3所示,模拟的地温场为自然条件下,未受人类扰动的地温场。模型自上而下,温度值由标高23.62m的12.9℃逐渐增加到标高-2500m的64.4℃。
图3 三维数值模型及初始模型模拟结果(网格线为xyz坐标轴;抽水井位置)
预测模型主要是预测抽水井水位到-1400m标高时地温场的变化,以及可能发生热害的标高位置。
根据模拟结果显示,对比不同水位下抽水井和钻孔JFZK1不同深度的温度值,-400m和-1400m标高水位的抽水井温度值均大于初始水位抽水井不同深度的温度值;钻孔JFZK1不同深度的温度值在初始水位、-400m和-1400m标高水位差别不大。
模拟区大概会在-720m的标高位置达到我国《矿山安全条例》中规定的井下工人工作最高温度28℃。
(1)抽水对于地温场会产生影响,但是影响非常小,基本可以忽略。并且,抽水对地温的影响与距抽水井距离有关,距离越远,地温场所受影响越小。
(2)结合实际测温孔数据及模拟结果,焦家金成矿带矿区开采标高从-625m到-1090m都有可能产生地温热害,标高从-720m到-1090m产生热害的概率更高,在实际开采过程中应密切注意,采取相应防护措施。