1. 引言 - beplay安卓登录
近年来为实现三维地质分析与可视化预测等,三维地质建模技术兴起并得到发展
[1]
-
[3]
。目前基于钻孔数据生成实体三维模型的算法不断被提出,如Kriging插值算法
[4]
、Shepard插值算法、最近临点(Nearest Neighbor)插值算法
[5]
、Fastrbf插值算法等,根据不同的地质分布特征可选取不同的插值算法。
目前国内学者针对三维地质建模已经做了许多研究,盛世强
[6]
等使用GOCAD建立三维地质模型,通过对模型特征的详细分析,识别出了构造带及其周缘地区褶皱变形样式,以及潜山内幕断层的交切关系。雷赟
[7]
等利用EVS建立的地质模型为解决污染物的控制以及合理有效开采地下水提供了模型支持。陈兵
[8]
等介绍了EVS的三种建模方法,及各种方法适用的地质条件。周念清
[9]
等提出了一种能同时表达复杂透镜体和清晰地质层面的EVS地层–岩性建模方法,能够对比较复杂的地形进行层序划分。孟美杉
[10]
等利用EVS建立出能准确刻画地下水污染状况的概念模型。
通过对上述学者的研究成果分析发现,目前的EVS地质建模是根据钻孔数据插值模拟出的地貌情况,不能真实地反映实际地貌与地层倾向,因此在复杂地貌情况下采用传统EVS建模法建立的模型无法满足实际需求。本文借鉴前人已有的研究成果,采用EVS软件地质建模与无人机正直摄影测量相结合的方式,提出一种EVS耦合地形–优化钻孔数据的建模方法,对贵州金象煤矿采煤区的地质情况进行三维地质建模,并采用该模型对采煤区的地层分布情况进行分析和统计采煤区的煤层储量情况。
3. 耦合地形–优化钻孔数据建模方法
耦合地形–优化钻孔数据的建模流程如
图3
所示。单独采用地层建模或者岩性建模无法反映出实际的地貌情况,为了取得更真实且能够满足实际需要的模型,可考虑在地层建模的基础上导入无人机正直摄影测量
[11]
[12]
的地貌数据,得到能真实反映地貌与地层倾向的模型。为此,提出了以下具体思路。
1) 根据无人机正直摄影的照片生成带坐标、高程的数字高程模型(DEM)
[13]
[14]
。
2) DEM文件转换为EVS软件能识别的gmf格式文件。
3) gmf格式文件导入优化钻孔数据的建模方法建立的模型中,生成能够满足实际需要的模型。
3.1. 优化钻孔数据的EVS地质建模方法
地层形成后往往会受到外力作用的影响,地形地貌会发生较大的改变。特别是在风化与流水作用下导致地层不断被侵蚀,形成低洼和高地地区,地层将出现断层如
图4
所示,地层在断层处会出现露头。采用传统地层建模在进行层序划分时无法划出地层的露头点,建立的模型在露头处会出现地层向下歼灭的情况,为了防止三维模型中出现地层向下歼灭,需用直线内插法在钻孔数据中补充出低洼地区被侵蚀的地层。如
图4
所示根据A、C钻孔中2号地层标高,利用直线内插法计算出B钻孔中未被侵蚀前2号地层的标高,以此补充出未被侵蚀前B号钻孔中的完整地层信息。
Figure 3
Figure 3. Coupled terrain-optimized borehole data modeling flow--图3. 耦合地形–优化钻孔数据建模流程--
Figure 4
Figure 4. Map of eroded landform--图4. 侵蚀地貌图--
“直线内插法”又称“数学内插法”,其原理如下:以2号地层为例,2号地层在C钻孔中坐标为C2 (x1 , y1 );在A钻孔中坐标为A2 (x2 , y2 );2号地层未被侵蚀前在B钻孔中坐标为B2 (x, y),B2 在A2 、C2 点确定的直线上。A2 、B2 、C2 三点共线,则(y − y1 )/(x − x1 ) = (y2 − y1 )/(x2 − x1 ) = 直线斜率,变换即得所求地层标高公式:
y
=
y
1
+
y
2
−
y
1
x
2
−
x
1
(
x
−
x
1
)
(1)
3.2. 耦合地形建模方法
优化钻孔数据的建模法建立的地质模型仍无法反映实际地貌情况,还需将无人机正直摄影测量的地貌数据耦合到三维地质模型中得到能真实反映地貌与地层倾向的模型。
随着无人机技术的快速发展,无人机低空摄影测量由于精度高、效率高、成本低廉的特点在勘察工程中得到了广泛应用。无人机正直摄影测量是一种通过真实的图像,丰富的色彩客观地反映真实的地貌状况,具有广阔应用前景的新型数字测绘技术。
1) 通过无人机正直摄影测量进行数据重建,生成数字高程模型,数字高程模型中包含有坐标、高程信息的DEM栅格文件。
2) DEM栅格文件导入Arcgis软件如
图5
所示,在Arcgis软件中可以将DEM栅格转为点,每个点都有地貌的坐标和高程信息,再将点文件导出为CAD文件格式如
图6
所示。
3) 采用CAD工具栏中数据提取工具提取CAD点文件的X、Y、Z信息如
图7
所示,提取结果保存为txt文本格式如
图8
所示。
4) 将txt文件后缀名改为gmf后缀名,并重新定义文件属性如
图9
所示,定义为units m (单位:米)、surface (面),使EVS软件能识别。
Figure 5
Figure 5. DEM raster image--图5. DEM栅格图--
Figure 6
Figure 6. CAD point file diagram--图6. CAD点文件图--
Figure 7
Figure 7. Data Extraction Steps--图7. 数据提取步骤--
Figure 8
Figure 8. Data extraction results--图8. 数据提取结果--
Figure 9
Figure 9. Attribute redefinition--图9. 属性重定义--
3.3. 耦合地形–优化钻孔数据建模的程序模块
地貌数据经过处理后成为能被EVS软件识别的gmf格式,将gmf格式文件导入建立好的模型中,使用模块如
图10
所示:gridding and horizons#1、gridding and horizons、combine horizons、horizons to 3d、explode and scale、file statistics、plume、viewer、post samples、legend。gmf地形文件导入gridding and horizons#1模块,gmf钻孔数据文件导入gridding and horizons模块,combine horizons模块用于耦合地形数据与钻孔数据,horizons to 3d、plume模块用于生成三维地质模型, viewer模块用于显示三维地质模型,post samples、legend模块用于显示钻孔数据信息,file statistics模块用于显示地形点信息。为了保证模型的高精度,gridding and horizons#1模块导入地形文件后需要与gridding and horizons模块坐标统一。
Figure 10
Figure 10. Coupled terrain-optimized borehole data modeling application module processes--图10. 耦合地形–优化钻孔数据建模应用模块流程--
4. 工程案列
4.1. 工程概况
金象煤矿隶属于织金县板桥乡管辖,2023年5月,位于织金县板桥镇、以那镇两镇交界处发生危岩垮塌,对矿区生成经营造成严重影响,需对矿区北侧危岩崩塌区进行应急抢险治理。经查询地质资料得知,危岩崩塌区下存在煤层,为科学化实施抢险治理,需对煤层分布情况进行探明。
矿区段内勘探钻孔揭露的地层主要有:三叠系下统飞仙关组1段、三叠系下统飞仙关组2段、三叠系下统飞仙关组3段、二叠系上统长兴+大隆组、二叠系上统龙潭组、6号煤层、7号煤层、14号煤层、16号煤层、27号煤层、30号煤层、二叠系上系统峨眉山玄武岩组。矿区内钻孔示意图如
图11
所示。
Figure 11
Figure 11. Drilling diagram--图11. 钻孔图--
4.2. EVS建模分析
为了能更好地探明煤矿区内的地质情况和煤层分布情况,分别采用地层建模与耦合地形–优化钻孔数据的建模方法进行建模,并将两种地质建模方法进行比较和阐明建模效果。
1) 地层建模方法。地层建模能够根据钻孔数据很好的进行层序划分,但在地层出现露头情况时,如
图12
所示,进行层序划分无法划出地层的露头点,建立的模型会出现地层向下歼灭的情况,如
图13
所示,这与实际地层倾向不符。且该方法建立的模型地形是根据孔口高程采用插值算法计算得到的,很难真实地反映矿区的地形地貌情况。
Figure 12
Figure 12. Ground investigation section--图12. 地勘剖面--
Figure 13
Figure 13. Stratigraphic model section--图13. 地层模型剖面--
Figure 14
Figure 14. Stratigraphic model--图14. 地层模型--
2) 耦合地形–优化钻孔数据的EVS地质建模方法。地层建模建立的剖面图与地勘剖面图进行比较,具有较大的差距。采用耦合地形–优化钻孔数据的建模方法建立的地质模型,与实际地层倾向一致且具有真实的地形地貌,如
图15
所示。将地勘剖面图和模型剖面图进行对比如
图12
、
图16
所示,可以看出两者的轮廓和地层倾向情况是较相近的。
Figure 15
Figure 15. Coupled terrain-optimized borehole data model--图15. 耦合地形–优化钻孔数据模型--
Figure 16
Figure 16. Coupled terrain-optimized borehole data model profiles--图16. 耦合地形–优化钻孔数据模型剖面--
4.3. 模型应用
EVS软件具有开挖功能,能根据模型计算出危岩崩塌区内各个地层的体积,使用模块如
图17
所示:import cad、create tin、distance to surface、select cell sets、intersection、volumetrics、format string、titles。
Figure 17
Figure 17. Excavation flow chart--图17. 开挖流程--
在煤矿危岩崩塌区内,三种建模方式计算的煤层储量与传统地勘方式勘测的煤层储量进行对比,详见
表1
。两种建模方式的相对误差值:地层建模6、7号煤层的误差分别为22.31%、6.89%;耦合地形–优化钻孔数据建模6、7号煤层的误差分别为11.19%、2.99%。耦合地形–优化钻孔数据的建模方法具有较低的误差值,证明该方法能更好地符合实际情况。
Table 1
Table 1. Coal seam reserve tableTable 1. Coal seam reserve table 表1. 煤层储量表
方法
6号煤层储量(万t)
7号煤层储量(万t)
传统勘探
30.533
69.570
地层建模
37.345
74.362
耦合地形–优化钻孔数据建模
33.950
67.486