矿井物探Word文档格式.docx
- 文档编号:21050376
- 上传时间:2023-01-27
- 格式:DOCX
- 页数:24
- 大小:617.70KB
矿井物探Word文档格式.docx
《矿井物探Word文档格式.docx》由会员分享,可在线阅读,更多相关《矿井物探Word文档格式.docx(24页珍藏版)》请在冰豆网上搜索。
(AB为供电电极,MN为测量电极,I测量供电电流,△U测量电位差)
具体测量方法为:
首先以固定点距沿巷道测线布置一系列电极,相邻电极间距为
,取装置电极距
,将相距为
的一组电极排列经转换开关接到仪器上,通过转换开关改变装置类型,一次完成该测点各种装置形式的视电阻率
观测(电极排列中点为测点或记录点,记录深度为
,图2所示为
时
,对于6号8号10号12号电极组成的排列,9号点是该排列中点,即为记录点,记录深度为
);
一个测点观测完后,通过开关转换到下一相邻测点对应的电极,以同样方法进行观测,直到电极距为
的整条剖面观测完为止。
改变电极距,重复以上观测,直到
测量结束。
图2.2高密度电阻率成像法测量断面示意图
点距
的选择,主要依据探测精度要求,精度要求越高,
应越小。
当
确定后,最大电极距
决定于预期探测深度,深度越深,
要求越大,但一般隔离系数
最大值不超过15为好。
煤矿井下工业电干扰极其严重,在测量方法上要采取一些技术措施,保证采集的数据准确可靠。
另外,在高密度电阻率成像中,建议选择温纳对称四极的供电和测量方式,因为温纳四极排列对于旁侧不均匀影响不像偶极或微分排列那样敏感,由其测量结果经过统计处理所形成的地电断面已经是非常光滑的地下电阻率的近似分布。
2.2数据处理与资料解释
高密度电阻率成像法数据处理主要是电阻率成像,一般利用专用处理软件包。
根据测区岩层的正常电性特征,把电阻率成像断面色谱图分成正常区和异常区。
对于异常区,按异常幅值的大小和高低划分高阻异常和低阻异常,高阻异常一般是由局部岩性变化或不含水的断裂破碎带引起,这种异常体若在采矿活动中形成不了突水通道,则对回采无不良影响。
低阻异常一般是由含水体(如溶洞)、含水的断裂破碎带以及局部岩性变化引起,对于异常体为含水体时,回采工作面中一定要提前作好防治水工作,避免突水事故的发生。
2.3应用实例[3]
某煤矿3501工作面煤层赋存条件好,煤厚5.0~6.2m,为气煤,是矿井下一步生产主要工作面,由于该面邻近受奥灰承压水威胁的采区,开拓过程中发现断层发育,且部分地段底板非常潮湿并向外渗水。
为了确保该面安全开采,需要查清工作面底板岩层结构弱面及富水性,为此,在已开拓的回风巷、运输巷分别进行高密度电阻率成像法探测。
图2.33501工作面回风巷和运输巷底板高密度电阻率成像断面图
(a)回风巷底板电阻率断面图(b)运输巷底板电阻率断面图
现场数据采集采用温纳三电位电极装置,测点距6m,测线总长1Km,最大隔离系数n=15,最大探测深度90m,经孔旁测深对比校正系数取0.7,有效深度大于60m。
数据采集过程中,工业电干扰严重,采取一系列技术措施,确保原始数据质量可靠。
图2.3(a,b)分别为3501工作面回风巷和运输巷高密度电阻率成像断面色谱图,两个断面的电阻率分布形态相似,色谱图电性分布全面反映出巷道底板以下岩层的岩性分布,浅部电阻率偏高,主要为3下煤层的反映;
中部以下基本上反映了煤系砂岩电性特征。
从图2.3(a)可以看出比较明显的两处低阻异常:
测点80m~105m区段,中部有一较低阻区;
测点160m~260m区段,断面深部有一明显的喇叭形低阻区,其视电阻率小于35
,低于正常岩层的电阻率,解释为存在隐伏含水体,异常区域较大,通道比较明显,与浅部未连通,但不能排除采矿活动的后期破坏会造成该隐伏含水体的再导升。
图2.3(b)中没有比较明显的低阻异常存在,这一趋势与工作面底板岩层倾伏方向一致。
3.1层析成像计算原理
70年代以来,由于CT技术(计算机辅助层析成像技术)在医学领域得到广泛应用,发展迅速,在技术和计算方法方面都已比较成熟。
80年代开始,CT技术在地球物理领域中应用逐渐开展起来,取得了许多令人瞩目的成果[4]。
CT技术的图象重建计算方法从数学上有变换法和代数迭代法两大类[5]。
目前在地学中主要应用代数迭代法进行图像重建。
1970年Gordon等人提出代数重建法,其基本思想是依据射线原理。
首先对成像条件提出一个初始模型,然后把模型网格化,计算出投影函数的观测值与理论值的残差量,将每条射线的残差量以它穿经每一网格的路径长度为权分摊到该网格中去,修正模型,反复迭代直到满足收敛条件为止。
工作面电磁波透视法采用偶极子天线发射,在介质中任意点的磁场表达式可表示为:
(3.1)
式中H为实测场强值,H0为理论初始场强值,
为介质吸收系数,r为发射点到接收点距离,sin
为方向性因子,一般可认为等于1。
图3.1工作面电磁波CT成像单元离散示意图
如图3.1所示,把工作面划分成有不同吸收系数的若干小单元(像元),每一小单元内可视为介质均匀的。
假设电磁波的第i个传播路径为
,它可以表示为若干小单元的距离之和:
(3.2)
没有射线穿过的小单元,可视dij=0,于是公式(3.1)变成
(3.3)
对(3.3)式两端取对数有:
(3.4)
若在多个发射点上对场强分别进行多重观测,便可形成矩阵方程:
[X][D]=[Y](3.5)
式中:
[X]是
未知数矩阵;
[D]是
系数矩阵,[Y]是已知数矩阵,即实测值。
利用SIRT算法(SimultaneousIterativeReconstructionTechniques同时迭代重构技术),对方程(3.5)进行计算,可以反演各像元吸收系数值,从而实现了工作面成像区内吸收系数反演成像。
利用反演计算结果可以绘制成像区吸收系数等值线图和色谱图。
为了确保电磁波衰减系数的高精度成像,必须保证像元尽可能小,由于数据不完全,当像元小到一定程度,势必造成像元内无射线通过,这样,该像元就失去意义。
为克服这一不足,可以先用大像元成像,然后分割成更小像元再继续成像,确保成像精度。
3.2工作面电磁波CT工作方法
工作面电磁波CT工作方法与坑透常规工作方法基本上一致,主要采用定点法观测方式,即发射机相对固定于某巷道事先确定好的发射点位置上,接收机在另一巷道一定范围内逐点沿巷道观测场强值,观测射线呈扇形分布。
CT工作方法要求在数据采集中测点应适当加密,观测密度应加大。
为了保证层析反演成像的高精度,坑透CT工作方法要求在坑透施测范围内保证足够的射线密度的同时,射线分布基本均匀,避免出现“盲区”,因此,在工作之前,要根据工作面的实际情况预先画出观测系统,看观测射线分布是否均匀,并适当调整发射点和接收点位置,从关键性环节确保成像质量。
构造复杂区应适当加密射线数量。
不同发射点的发射参数要基本上保持不变,以减少吸收系数反演的误差,保证收敛速度。
3.3解释方法
工作面电磁波CT解释主要依据吸收衰减系数成像结果(包括吸收衰减系数色谱图和等值线图)来进行解释。
不同的地质条件和不同的地质现象所引起的电磁波衰减特征不同,因此,根据不同的异常特征即可以进行地质解释。
一般说来,正常煤层对电磁波的吸收较围岩对电磁波的吸收小,而当煤层破碎、煤层中裂隙发育以及裂隙中含水时,就会造成吸收衰减增大,甚至其引起的衰减较围岩大得多。
工作面中隐伏断层视其落差大小引起的衰减也不一样,落差大于煤厚的断层较落差小于煤厚的断层引起的衰减大,落差小于煤厚的断层引起的衰减随落差的减小而减小。
当煤层厚度变薄以及煤层中存在夹矸时,电磁波的衰减随煤层变薄和夹矸增厚而增大。
当煤层有火成岩侵入体存在时,其衰减依据火成岩的电性特征不同而不同,出现低值异常和高值异常都是可能的,但对同一采面内的侵入体具有单一性。
工作面中陷落柱同样引起高吸收衰减。
从吸收系数异常范围来说,断层引起异常范围一般呈条带状,落差大的断层较落差小的断层的异常带宽。
煤层厚度变化区、夹矸增厚区、煤层破碎区的异常范围一般都较大。
工作面中侵入体及陷落柱引起的异常范围与其形状基本上一致。
根据上述吸收系数的异常特征,就可以进行不同地质异常的地质解释。
图3.27513工作面电磁波CT成像吸收系数色谱图及解释成果
(a)吸收系数色谱图(b)解释成果
3.4应用实例
3.4.1复杂断裂构造的探测[6]
某矿7513工作面走向长650m,倾倾斜长145m,煤层厚度5.1~6.3m,煤层赋存稳定。
该工作面南北两端构造简单,但中部构造复杂,巷道中已揭露断层断点10个,断层之间的切割关系难以推断。
为了查明这些断层向采面内的延伸情况、断层之间的切割关系、断层最大落差位置以及是否存在隐伏断层,采用工作面电磁波CT成像方法进行探测,共发现异常14个,解释断层12条(如图3.2),工作面电磁波衰减系数CT成像色谱图清晰地显示出中间区段断层的切割关系以及最大落差位置,校正了原推断结论,电磁波CT工作成果为该面的正常生产提供了准确可靠的地质依据,主要解释成果为回采证实。
图3.33392工作面电磁波CT成像吸收系数色谱图及解释成果
(a)吸收系数色谱图(b)解释成果(——实际揭露,----CT解释结果与实际偏离位置)
3.4.2煤层结构变化及煤厚变化的预测[7]
某矿3392工作面巷道揭露局部地段煤层夹矸增厚,煤层变薄,为了圈定煤厚变薄带,指导下一步生产,采用工作面电磁波方法进行探测。
图3.3为该面电磁波CT成像吸收系数色谱图和解释成果,从图3.3(a)中可以清楚地看出工作面中部靠近运输巷一侧存在一明显的较高吸收区段(4#异常),解释为煤层变薄带,煤厚相对变薄0.4~0.7m,工作面推进中实际揭露煤厚变薄区位置与CT解释一致,最大相对变薄0.6m;
1#、2#异常区为煤层变薄区,实际揭露最大相对变薄0.4m,为底鼓引起。
另外,CT解释的七条断层(见图3.3(b))与实际揭露位置基本一致,准确率达86%。
该面夹矸厚度变化一般小于0.4m,因而成像结果反映不太明显。
4.1弹性波CT探测原理
如图4.1所示,在钻孔1放置发射换能器T1,T2,T3,…,在钻孔2放置接收换能器R1,R2,R3,…。
根据接收到的弹性波初至走时数据来反演两孔间的慢度分布
或速度
。
假设弹性波的第i个传播路径为
,其旅行时为
,则
(4.1)
上式是一曲线积分,ds是弧长微元。
式
(1)中
和
都是未知的,
是已知的。
这是一个非线性问题,在速度场变化不大的情况下,可以近似地把射线路径看作是直线,即假设
为直线,而只反演
把反演区域离散化,离散成
个小单元(像元),在每个单元内可视为介质均匀的,速度是常数,按一定顺序记为
其相应的倒数即慢度记为
,这样,第i个射线的旅行时表示为:
(4.2)
图4.1CT成像数据采集示意图
其中,
是第i条射线穿过第j个像元的长度(或距离)。
当有大量射线(如m条射线)穿过反演区域时,根据(4.2)式就可以得到关于未知量
的m个方程,形成一矩阵方程
(4.3)
其中:
[A]=(aij)m×
n称为距离矩阵,[T]=(ti)m×
1为走时向量,即测量值,[S]=(Si)n×
1为慢度列向量。
上式是大型稀疏矩阵方程,求解方法较多[2],采用收敛性好的SIRT法(同时迭代重构技术SimultaneousIterativeReconstructionTechniques)。
慢度的初值取射线的慢度均值
(4.4)
式中
为总射线数目,
分别为任意条射线的走时和射线长度,计算过程中第
次迭代像元
上的取值
为
(4.5)
其中
是
中第
列非零元素的个数,
为根据慢度近似值
计算的走时,
,
为驰豫参数,
,相当于对慢度修正量给定权系数,加快迭代初期收敛速度,增加计算的稳定性[9]。
为了保证反演结果的精度和分辨率,采用弯曲射线追踪方法,直接应用费马原理和Dijkstra最佳路径算法[10]。
利用SIRT法求解(4.3)就可以得到各单元离散慢度分布,相应地离散速度分布也已知,从而实现了钻孔1与钻孔2之间的速度场反演成像。
利用计算值便可绘出成像区域速度分布图。
下面结合煤层开采引起的底板破坏深度探测的例子[11],探讨弹性波CT技术的探测方法和探测效果。
4.2探测方法
测试方法采用一孔发射,另一孔多点接收的弹性波CT透射方式,其观测系统如图4.2所示。
图4.22#-4#孔观测系统图
仪器采用SD-1声波仪、ATF-1A型大功率发射机和换能器(频率5KHz)。
共进行了三对孔穿透测试:
2#-4#、5#-8#、4#-8#,发射点间距4m,接收点间距1m。
2#-4#孔位于A钻窝,两孔形成的断面垂直煤层走向;
5#-8#孔位于B钻窝,两孔形成的断面垂直煤层走向;
4#和8#孔分别位于A和B钻窝,间距7.2m,两孔形成的切平面与煤层斜交,交角约
上述三对钻孔形成立体空间控制煤层底板,如图4.3所示,因而利用弹性波CT速度场成像技术可以从三维空间中去观测底板在采动矿压作用下的应力变化和破坏规律。
图4.3测试钻孔工程布置平面图(a)和剖面图(b)
4.3测试结果分析
按采矿活动的不同阶段对三个断(切)面分别进行了10次弹性波CT透视动态探测,速度成像结果不仅较好地反映了断面的岩性变化及结构面发育状况,而且较全面地反映了底板岩体的应力变化和破坏程度及深度。
下面以2#-4#孔间断面为主进行分析。
图4.4为2#-4#孔间断面不同时段((a)采前工作面距2#-4#孔断面20m、(b)工作面推过2#-4#孔断面3m、(c)工作面推过2#-4#孔断面32m)三次CT测试速度场分布图。
从图中可以看出断面岩体的岩性分布和结构面状况,断面浅部速度比较均匀,一般4400m/s~4600m/s,岩性为细砂岩;
断面中部存在一明显低速层,速度为3300~3400m/s,层厚约0.7m,岩性为泥岩;
断面深部速度变化较大,一般3800m/s~5000m/s,下部靠近4号孔速度较低,上部靠近2号孔速度较高,从钻探资料知岩性为细砂岩,且4号钻孔岩芯有破碎迹象,说明速度场分布准确地反映出上部岩石比较致密,而下部岩体裂隙发育,原生结构面较多。
当工作面刚推过,观测断面正处于卸压区,此时的底板岩体速度分布发生了很大变化,如图4.4(b)所示,左上部岩体速度普遍降低,特别是在浅部细砂岩速度降低幅度较大,达1100m/s,说明这部分岩体裂隙发育,且以纵向裂隙为主。
中间泥岩层及其下层细砂岩上部速度也明显降低,再往下部岩体速度变化较小,这就反映出在此阶段底板岩层应力变化较大,主要体现在上部岩层中应力释放,结构面增多,裂隙发育,岩体遭受破坏。
图4.42#-4#孔间断面不同时段声波CT测试速度场分布图
(a)采前(b)采中(c)采后
随着采面的继续推进,顶板岩体冒落,底板岩体经历采后再压缩过程,图4.4(c)为采后底板岩体基本稳定后速度场分布,它表征了岩体最终的裂隙发育状态,可以看出此阶段的速度场分布比图4.4(b)中速度场分布在浅部岩层中略有增大,说明岩体裂隙受压后孔隙度有所降低,而在中下部略有减小,按速度减小20%(相对于采前正常场)且相对稳定的区域划分破碎后导水裂隙发育区,可以得到2号孔孔深30m以浅位置和4号孔孔深27m以浅位置划定为底板破坏导水裂隙发育范围,相对煤层底板垂直深度11.6~12.3m。
利用弹性波CT探测煤层开采引起的底板破坏深度直观准确,效果可靠,实现了面积测量,特别是探测结果反映出底板破坏的动态变化规律,对于承压水上采煤具有重要的指导意义。
弹性波CT探测方法对于采矿活动引起的其它岩体破坏,如覆岩变形破坏、巷道开挖围岩破坏等,探测方法同样适用。
地质雷达(又称探地雷达)是利用超高频(106~109HZ)脉冲电磁波探测地下介质分布的一种地球物理探测方法。
实践表明,地质雷达是一种高分辨率探测技术,它可以分辨地下10-1m尺度的介质分布,可以对浅层地质问题进行详细分层调查,也可以对地下浅部掩埋目的体进行无损探测[12]。
随着计算机技术及数字处理技术的发展,地质雷达在工程地质勘察,灾害地质调查,公路工程质量的无损检测,考古调查以及工程施工质量监测等诸多领域中得到广泛应用。
近年来,在国内外全面开展了地质雷达技术用于矿区井下探测顶、底板及回采工作面前方小断层、老窑、岩溶分布及探测煤厚、充水小构造、陷落柱和巷道围岩松动等地质问题的研究工作,并取得了较好的地质效果[13]。
5.1探测原理
探地雷达是利用高频电磁波(主频为数十兆赫至数百兆赫)以宽频带短脉冲形式,由地面通过天线T送入地下,经地下地层或目的体反射后返回地面,为另一天线R所接收,当地下介质中的波速v为已知时,可以根据测到的精确时间t值,求出反射体的深度(m),如图5.1所示。
电磁波在介质中传播时,其电磁波强度与波形将随所通过介质的电性质及几何形态而变化。
因此,根据接收到波的旅行时间(亦称双程走时)、幅度与波形资料,可推断介质的结构或地质体的空间位置。
雷达图形常以脉冲反射波的波形形式记录,波形的正负峰分别以黑、白表示,或者以灰阶或彩色表示。
这样同相轴以等灰度或等色线即可形象地表征出地下反射面。
图5.1探地雷达原理示意图
5.2工作方法
利用地质雷达进行探测时,根据探测目标和位置的不同,通常采用以下几种方法。
5.2.1同位发射——接收法
在同一点位设置发射天线和接收天线,同时发射脉冲信号和接收目标反射信号,根据回波信号走时来计算目标距离和位置。
在井下掘进工作面超前探测时,发射点(T)和接收点(R)相距在0.5~1m的位置上架设天线。
架设天线时,要先铲平工作面岩层,将天线紧贴煤壁或岩壁,并用金属网将四周围好,避免漏场。
5.2.2剖面法与多次覆盖法[12]
(1)剖面法
发射天线(T)和接收天线(R)以固定间距沿测线同步移动(图5.1)。
发射天线和接收天线同时移动一次便获得一个记录。
当发一一收天线同步沿测线移动时,就可以得到由多个记录组成的地质雷达时间剖面图像。
横坐标为天线在地表测线上的位置,纵坐标为雷达脉冲从发射天线出发经地下界面反射回到接受天线的双程走时。
这种记录能准确地反映测线下方各个反射面的起伏变化。
这种方法能在地面施工、井下巷道底板探测和侧壁探测中使用。
在使用同位发射——接收法和剖面法施工时,反射电磁波的旅行时间由下式计算:
(5.1)
为回波双程走时,ns;
VG为电磁波在地层中的传播速度,m/ns;
S为发射天线与接收天线的距离,m;
d为反射界面的法线深度,m。
(2)多次覆盖(共深点)方法
地质雷达探测来自深部界面的反射波时,由于信噪比过低,不易识别回波。
这时可采用类似于地震的多次覆盖技术[14],应用不同天线距的发射——接收天线对同一测线进行重复测量,然后把所得的测量记录中测点位置(共深点)相同记录进行叠加处理,能增加所得记录对地下介质的分辨率。
5.3仪器设备
地质雷达仪多采用脉冲发射,而脉冲波有调制方式和直接脉冲方式两种。
目前使用的多为时域地质雷达仪,主要有加拿大EKKO系列、美国SIR系列和瑞典RAMAC系列,这些系列为调制方式脉冲波;
国产仪器主要是煤科院重庆分院KDL系列,为直接脉冲方式。
5.4应用实例[14]
某煤矿地表三采十一层煤露头区及七采二、四,六层煤露头区位于小汶河北岸,由于盗采盗挖及滥取沙土,造成大量盗采坑和采沙坑,形成地表水泻入井下的良好通道。
而且六层煤底板一灰与十一层煤底板三灰在岩溶裂隙发育地段也会形成地表水导入井下的主要途径。
特别是在雨季小汶河水倒灌极易造成溃水淹井事故,严重威胁矿井安全生产,增加了矿井的排水费用。
因此,拟在此区进行灰岩浅截帷幕注浆工程进行治理。
但由于灰岩岩溶裂隙发育地段分布不清,给注浆布孔带来困难。
通过采用探地雷达探测,查明了煤柱破坏情况和灰岩岩溶裂隙发育地段分布范围,成功地解决了这一问题。
5.4.1工作方法及工程布置
采用瑞典RAMAC/GPR型探地雷达,通过条件实验,选用25MHz探测天线。
发射天线与接收天线间距为4m。
采样时窗为2000ns,可以保证探测深度达50m,点距为0.5m,重叠次数128次,采用反射剖面法进行工作。
根据探测目的和任务,共设计布置2条探测线。
5.4.2标志层的雷达波反射特征
通过对全区雷达时间剖面与已知地质资料的对比分析,本区主要地层介质的雷达反射波特征如下:
①第四系底界面反射波第四系底部大部分为流砂层,与下伏的不同风化程度的砂岩电性差异较大,能形成较强的反射波,该反射波的能量与下伏岩石岩性及其风化程度关系较大,基岩风化越严重,则两者电性差异越小,该反射波的能量则越弱。
如图5.2所示
②灰岩反射波测区内浅部地层中需探测的灰岩主要有三灰、一灰,三灰为泥质灰岩,厚度0.7~1.2m,是十一层煤的底板,其下为砂质页岩。
三灰与煤11电性差异不大,二者能形成一组复合波,在煤层厚度及结构变化不大的情况下,其能量强弱主要是三灰的岩性结
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 矿井 物探