本站小编为你精心准备了测氡法探测采空区的应用参考范文,愿这些范文能点燃您思维的火花,激发您的写作灵感。欢迎深入阅读并收藏。
《煤矿开采杂志》2014年第三期
1奇异性分析
在低放射区域,场的局部变化性对于识别异常场具有重要意义。对于元素含量进行差值时应尽量保持局部奇异性。常规方法并不能有效识别少数奇异数据所反映的局部异常。奇异性分析法以单个测点作为计算对象,其结果具有明显的局部灵敏性。作为对含量-面积法使用的补充,文中引入奇异性分析法[10]。奇异性分析法的提出同样基于分形理论。具有自相似性的地质过程,在改变度量尺度时,场依然保持相似形。其可以由以下指数函数来表达C(r)∝ra(2)式中,C(r)表示基于尺度r的邻域内量或场;a为奇异性指数。此处,C(r)可以表示元素的平均面含量。如以元素平均密度ρ(r)代替C(r),则式(2)改写为ρ(r)∝ra-2(3)不同的地质作用具有不同的自相似性,对应不同的奇异性指数。当a=2时,ρ(r)为常数,此时表示元素平均密度为常数,与度量尺度无关,对应没有奇异性的区域,为背景场。当a≠2时,为异常场。数据处理中,应用窗口法[11]。从测区原始氡值分布图上任意点开始,选定窗口形态(圆形或正方形),半径r依次递增,对外进行窗口滑动,对每个窗口中氡值计算平均值ρ(r)。然后对ρ(r)-r取双对数,用最小二乘法对lgr-lgρ(r)进行线性关系拟合,所得直线斜率即为对(a-2)的估计。从而由其结果是否为零得出该点所在场。
2实例分析
以山西某公司场地采空区地面氡值为例进行数据处理。图1为探测区实测氡值分布直方图及正态分布图。由图可得出,测区氡值数据属于偏态分布,数据重心明显左偏。传统统计方法使用受到局限,应用含量-面积法对其进行数据处理。
2.1含量-面积法采用该方法对测区实测氡值进行数据分析,如表1所示。应用最小二乘法对lgr-lgN(r)进行线性关系最优化拟合,见图2。图2中,最终拟合得①,②,③3条直线,3条直线斜率不同,即分维数不同,反映了具有不同自相似性的地质作用场。拟合度R2≥0.8754,表示lgr与lgN(r)具有良好的相关性,表明测氡法探测采空区所得氡值分形特征明显。本次应用分形理论进行数据处理,最终拟合时得到2个交点。直线①与②的交点所对应的r值即为圈定该区采空区的氡值异常下限,而直线②与③不同的分维数反映了采空区中伴随有另一不同的地质作用,即煤矿采空的基础上伴随有陷落柱。对于陷落柱,相较柱外采空区,其柱体内较大程度的岩石碎裂,使碎岩间具有良好的连通性,更有利于氡气的释放和向上迁移,因此其地表附近氡浓度较其他区域高。同时,对于柱体未塌至基岩顶部的陷落柱(即封闭型柱体),地表附近氡值异常表现为幅值较高的单峰状异常,氡值等值线图上大体表现为封闭的等值圆圈簇,且极大值对应柱体中部。计算其交点对应r值,可得采空区氡值异常下限为502个计数/(3min),确定采空区中陷落柱的氡值异常下限为813个计数/(3min),由此确定采空区范围,见图3(探测区为18-26测线,1-63测点区域)。与该区域利用地震勘探、瞬变电磁法等综合物探方法所得采空区范围(图中线条圈定区域)相比,基本吻合(采空区域之间的连通性将在下节进行论述)。从而验证了含量-面积法运用于低放射区氡值数据处理的可行性。而基于传统统计方法确定的异常下限为814~918个计数/(3min),很明显,其划定的采空区范围只是探测区的相对高放射区域。在低放射区域,煤矿采空区氡值的高异常区域一般只发生于局部,如伴随有陷落柱的区域。而这一结果又与由含量-面积法所确定的地质构造或陷落柱区域相吻合。
2.2奇异性分析对比其他物探方法所划定的采空区范围,为了对测区内应用含量-面积法所圈定的采空区域之间的连通性做进一步检验,利用奇异性分析法的局部灵敏性,对该采空区遗漏区测点做奇异性分析。为了提高所测数据的使用率,从而保证计算结果的准确性,选取测线22为检验线(检测点依次选取测点5号、8号、11号…,45号测点附近检验点加密)。在此,ρ(r)表示基于尺度r的邻域内氡含量的平均值,计算窗口采用正方形,以检测点为中心,半边长r分别为15m,25m,35m,45m。(图3中测点点号、线号间距均为10m)。其计算结果见表2。当奇异性指数a=2时,ρ(r)变化不大,约等于常数,对应区域为背景场。当a≠2时,则对应异常场。与含量-面积法所得结果进行对比,结果见表3。
分析对照结果,检测结果基本吻合。除个别试验点之外,拟合相关系数R≥0.8631,可见,lgr与lgρ(r)指数关系成立,奇异性指数估计可靠。试验点22-23,22-38,ρ(r)≠常数,首先排除其存在于背景场的可能。另外,氡的运移受裂隙控制,碎岩的裂隙程度直接影响了地表积聚氡气浓度的大小。对于采空区内个别岩石碎裂程度相对较低的地方,反映在氡值等值线图上时则会表现为局部低值。图3中,试验点22-23,22-38均位于采空区内背景场与异常场交界处,岩石碎裂程度较其他地方复杂,氡向上迁移受到影响,地表所测氡值较真实值误差较大,而奇异性分析拟合过程对测量数据准确度要求较高,故拟合相关指数较低。特别是试验点22-44,22-45经奇异性分析均为异常场,从而表明,图3中由含量-面积法所确定的采空区域是连通的。该结论与探测区其他综合物探所测的采空区范围吻合,从而证明了奇异性分析法的局部灵敏性在处理氡值数据时具有明显优势。
3结论
(1)传统统计方法在早期处理高放射性区域异常数据时是一种较有效的方法,而在低放射性区域,其应用则具有一定的局限性。(2)应用含量-面积法圈定氡值异常区域,相较传统统计方法,能有效强化弱异常信息,简单行,但仍存在局部遗漏。而奇异性分析法对异常数值反映敏感,但数据处理过程较为繁琐。结合二者优势,首先运用含量-面积法对整体实验区数据进行分析、采空区划定,然后对边界部分进行奇异性分析。结果表明,试验结果较为理想。(3)在低放射区探测地下采空区的过程中,采用测氡法是可行的,其所得采空区范围与对该区进行的其他物探方法综合解译成果基本相符。
作者:徐萍张新军刘鸿福李峰梁桂玲单位:太原理工大学矿业工程学院新奥气化采煤有限公司