张 昭,张大明,殷全增,王进朝,王焕允,王春博,王 鹏,张龙飞
河北省煤田地质局物测地质队,河北 邢台 054000
2019年8月28日,以“太行名郡·园林生活”为主题的河北省第三届园林博览会在邢台开幕.园博园选址在邢东矿采煤塌陷区,总占地面积3 km2,作为邢台市中央生态公园中一颗璀璨明珠的园博园,实现了从煤矿开采塌陷区向城市核心绿肺的蜕变.
邢台作为资源城市,煤矿经过长期开采,形成了大量采空区,存在地面沉降、水土资源损毁和区域环境污染等问题.进行矿山生态修复工程需探查地下采空区范围及地面沉降情况.鉴于此,本研究通过收集分析以往的地震、电法等物探资料,确定研究区煤层底板赋存形态、构造发育情况及采空区范围,指导公园地面选址.利用遥感技术,开展无人机动态监测,获取矿山生态修复的基础数据,为治理过程提供影像资料,实现“物探+遥感”技术有效服务于采煤坍陷区矿山生态修复工作.
地震勘探是通过天然场或人工场激发的地震波在地下介质中的传播改变,确定地下地质体分布变化及岩石学弹性参数的物探方法[1].由于人工场地震勘探激发地点可控、激发时间可控、能量可控,地震勘探多采用人工场激发方式.其技术特点为:激发场产生地震波,预设地震检波器采集地震波信号,分析采集的地震波动力学和运动学特征(时差、振幅、波形、相位、频率等),搭建地震波地下传播路径,获取地下岩石构造及弹性参数信息.地震波分为纵波(P波)、横波(S波)和面波,岩石弹性参数不同波速不同,依据此可反映岩石性质.
地震勘探是勘测石油、天然气和固体矿产的重要手段,在煤田勘查、工程勘查、区域地质研究等方面广泛应用[2].目前,地震勘探方法在寻找地下构造、民用工程建设、地质灾害防治中发挥着重要作用,为城市建设服务.
电磁法勘探是利用天然交变电磁场来研究地球电性结构的地球物理勘探方法,通过研究介质电磁性差异,观测天然或人工的交变电磁场随空间或时间的变化规律[3-5].大地电磁法技术成熟,施工方便,成本低,受高阻层屏蔽影响小,低阻层分辨率高,在矿产勘查、地热调查、深部地质构造探究等方面也有很重要的应用.
可控源音频大地电磁法(CSAMT法)是一种人工源频率域电磁法,其方法原理与大地电磁测深相近[6],针对天然电磁场信号弱的特性,采用人工源发射方式,利用发射电偶极AB(1~3 km)向地下发送频率不同的交变电流,产生交变电磁场,在场源距离足够远的区域测量相互垂直的电场信号(Ex)和磁场信号(Hy),依据公式(Cagniard法)求取地层视电阻率和阻抗相位.CSAMT法采用人工源激发,野外采集数据质量优质,操作性强,利于强干扰区开展勘测工作[6-8].该方法处理解释简单,成果剖面横向分辨率高,且不受高阻屏蔽,在采空区调查等方面得到广泛应用[9].
遥感方法采集装置与测量目的有一定的距离[10].其采集探头安装在空中飞行器、飞机或卫星上,受地形影响小.数据采集阶段,采用不同的探头可以得到多种波段的电磁波信息,经软件处理得到地质信息.卫星遥感技术采集资料具有区域广、绿色、高效等特性[11].
地质构造在遥感图像上往往以图形、色调特征、地貌等表现,呈平直、微弯形的线性条带或者圆形、半圆形、椭圆形环状条带展示.电磁波辐射的差异导致不同的影像色调与形态,地表地貌或地质现象的异常可作为断裂构造的依据[12-13].
高分影像无人机(UAV)是利用遥控装置和控制系统进行飞行和工作的不载人式飞行器.通过搭载各种传感器,如可见光相机、热红外相机等,可以获取高分辨率遥感数据[14-16].无人机可以低空飞行,具有高机动性,能够多角度观测目标,快速获取目标地物的连续、高分辨率影像数据以及其他自然资源和环境条件等信息.随着无人机遥感技术的发展,其在生态监测、灾害勘察、精准农业等领域应用范围迅速扩大,对于矿山生态修复工程领域的应用也越来越重要[17].
冀中能源股份有限公司邢东矿位于邢台市东郊,地理位置优越,交通便利(图1).矿井设计年生产能力60×104t[18].2001年11月矿井建成投产,开采深度为+59.5~-1200 m,开采的3个煤层分别为山西组的2煤、2下煤和太原组的6煤.经过20多年的开采,形成了采空区,长时间的沉降作用,在采空区的上方产生了6.7 km2的沉陷区,最大地面沉陷值为2.9 m.目前沉陷区均在缓慢地沉陷,没有达到稳定状态.随着矿业活动的延续,矿山将形成新的采空区,采矿引起的采空沉陷及地裂缝将对地面村庄和城市建设产生严重危害.
图1 物探工作概况图Fig.1 Overview of geophysical prospecting work
邢台市政府经研究,决定在邢东矿采煤塌陷区建立中央生态公园——邢台园博园,有必要对采空区范围、地面破坏情况调查清楚.
多年来,为保障矿井安全开采和规避地质灾害,在该地区进行了大量的地面物探工作,其中河北省煤田地质局物测地质队分别于1999年、2007年、2010年、2014年在矿区进行了多次三维地震勘探工作,2018年开展了中央生态公园采空区电法勘查工作,2016—2019年,在邢台园博园开展了连续遥感监测工作,详见物探工作概况图(图1).多次工作查明了邢东矿煤层埋藏深度、地层起伏形态、构造发育、采空区范围等情况,为地质灾害防治、生态公园建设提供了科学数据支撑.
2.2.1 地震勘探结果
多年来,邢东矿多次开展三维地震勘探工作.以2014年邢东矿后备区为例,激发方式采用井中高爆速成型炸药激发,60 Hz高频检波器接收,24次覆盖单边下倾激发的观测系统,428XL数字地震仪0.5 ms采样录制,记录长度1.6 s.
勘探工作查明了区内新生界厚度变化、可采煤层的构造形态及赋存深度、断距大于5 m的断层、异常区和其他地质构造,能够准确划分采空区分布范围.从邢东矿后备区2煤层沿层切片图(图2)中可以看出,区内被多条近东西向张性断层切割,小型断层和裂缝较发育,构造面貌复杂.
图2 邢东矿后备区2煤层沿层切片Fig.2 Plan view of No.2 coal seam in the reserve area of Xingdong Mine
从后备区AA′剖面线图(图3)中可以看出,勘探区东北部2煤层已经开采,采空区域反射波频率低,能量弱,波形紊乱不规则,影响下部煤层的反射波.煤层开采后形成的采空区是造成地面塌陷的重要原因.
2.2.2 电法勘探结果
在邢东新区中央生态公园采用可控源音频大地电磁法开展物探勘查工作.中央生态公园占地面积17.478 km2,位于邢州大道以南、中兴大街以北、襄都路以东、东华路以西.邢东煤矿处于生态公园范围内,面积14.5149 km2,两者范围基本重合.
电法工作查明了区内的采空区位置、分布范围及积水情况.从D18线反演电阻率断面示意图(图3)中可以看出:上部电阻率较低,为第四系、新近系、古近系、三叠系、煤系地层的电性显示;下部电阻率较高,电性分布相对均匀,为奥陶系地层灰岩电性显示.图中,反演电阻率自上而下逐渐增大,反映出第四系、新近系、古近系、三叠系、煤系地层、奥陶系基底反演电阻率从浅到深逐渐增大的基本电性规律[19-20].
图3 AA′剖面线图Fig.3 AA′section line
图4中,-700~-800 m煤层赋存深度中,在960~1560 m点距之间出现相对高阻特征,与围岩电阻率差异明显.分析认为,该高阻异常区由煤系地层采空后引起,采空区富水性相对较弱.
图4 D18线反演电阻率断面示意图Fig.4 Resistivity inversion section along exploratory line D18
2.2.3 遥感监测结果
为在受采空区沉降影响的区域建立中央生态公园,需对涉及的村庄、厂房、土地进行重新设计和规划,达到矿山修复、绿色发展的目的.
场地调查内容包括生态修复区域内地表植被情况、土壤污染情况、地形地貌及地物等.无人机低空摄影获取高质量的遥感影像数据,获取地表基本信息,减少人工调查,提高工作效率和质量,对场地改造进行宏观把控.
高分遥感和无人机航摄可获取较高分辨率的正射影像图,可将生态修复区域的实际情况直观地展现出来.该技术可实施监测施工进度,加强生态修复施工管理,有助于把控施工进度,为了解生态恢复工程的完成情况提供依据.
生态修复工作是一个长期的动态过程.完成生态修复工程后,按照监测计划形成正射影像,对工程治理效果、地貌地物变化进行长期监测,能够帮助评估生态修复区的治理现状,推动生态修复工作进一步开展[21].
采用高分遥感和无人机航摄,得出2016年9月—2019年10月间邢台园博园建设情况变化(图5).2016年9月间,对区域内植被、建筑物等地貌地物等进行了观测;2017年—2019年8月间,进行了园博园建设全过程监测和工程实施进度监督;2019年10月,对园博园开放两个月以来地表植被变化进行了动态监测,可见绿色覆盖率明显提高,生态修复工程初见成效.
图5 邢台园博园遥感监测连续变化图Fig.5 Continuous changes in Xingtai Garden Expo Park by remote sensing monitoring
遥感技术在生态修复区获取基础资料、施工监管和长期监测方面均具有积极作用,作为新技术在采煤塌陷区改造修复进程中,为政府制定决策和生态修复结果评价提供基础数据.
在邢东矿开展地面物探和遥感工作,成果应用于矿山修复即邢台园博园建设工程.结果表明:
1)三维地震可准确查明地层起伏形态、构造发育情况、采空区分布范围;电法勘探能够圈定采空区范围、地层富水性状态;遥感技术可连续监测地表地物变化情况,为生态修复区域提供基础数据.
2)在进行矿山生态修复中,收集前期已开展的物探工作资料,详细分析数据,能够指导地质灾害防治、生态修复改造工作.