MAPGIS工作平台在煤层气地质评价中的应用

黄辉强，吴顺彬

（四川省地矿局113地质队，四川 泸州 646000）

石宝-邱家祠井田位于四川省古蔺县石宝镇、长坪乡。大地构造位置归属杨子地台内的V级构造单元，川黔娄山关断褶带北缘，古蔺复式背斜南翼的次级褶皱-石宝向斜东段。区内构造不发育，地层产状较缓，煤层分布稳定，煤系地层较厚，总厚为76.25～130.11m，煤炭资源丰富，具有较大的煤层气资源潜力。

1 研究目的及方法

煤层埋深增大，压力和温度均增加。储层压力是煤层含气量预测的重要参数。因此，精确地获取煤层埋深，能够较准确地估算煤层储层压力，合理地预测煤层气资源量。

传统获取煤层埋深的方法是通过统计区内钻孔中煤层的埋深，再通过软件制作埋深等值线图。由于钻孔一般间距较大，钻孔数量有限，数据量少。因此，通过软件拟合的等值线图不能真实反应我们整个井田煤层埋深情况。鉴此情况，笔者利用煤层埋深与地形等高线和煤层底板等高线的关系，通过MAPGIS工作平台中DTM分析的功能，制作煤层埋深等值线图。其能够较真实地反映煤层埋深状况，为石宝-邱家祠井田煤层气地质评价提供可靠依据。

图1 石宝-邱家祠井田C25煤层埋深等值线图

2 DTM分析获取煤层埋深流程

DTM分析系统中GRD功能是以数字形式按一定结构组织在一起表示地形特征空间分布的数字模型，是一系列地面X、Y坐标与相应高程Z所组成。笔者通过处理地形等高线和底板等高线，分别获取规则网格数据文件，再通过简单的减运算，得到新的网格数据文件，作出煤层埋深等值线图。

其具体步骤为：第一步，在MAPGIS的图形处理功能中分别对两类等高线赋高程值；第二步，在DTM分析的处理点线中进行线高程提取，提取线高程值（Z）；第三步，在DTM分析的TIN模型中进行离散点三角网格化，设置相同的数据范围和网格间距（确保读取的点位一致），分别生成 GRD文件；第四步，在DTM分析的GRD模型中把生成的GRD文件进行DEM差值运算，生成新的GRD文件；最后，在Mapgis工作平台的DTM功能中进行平面等值线绘制处理，保存生成点、线、区文件，对图像进行修饰校正，生成埋深等值线图（图1）。

3 区块埋深数据读取流程

通过MAPGIS辅助软件sction的属性功能中高程属性点操作，就能够把区块内任意点高程值赋给区块内对应点，导出点属性，就可以得到我们每个区块的埋深。

其具体步骤为：第一步，在矢量化的煤层气资源量估算图中，新建点文件，每个区块均输入点图元，编辑其属性结构，使其属性结构中某一字段与埋深等值线图中线属性结构一致，保存点文件；第二步，将生成的点文件添加到煤层埋深等值线图中，在sction的属性功能中高程属性点处理，此时，点被赋埋深数据属性；最后，在sction的辅助工具中导出点属性数据，处理数据表格，获取各区块埋深。

4 结论及建议

通过该方法获取的各区块中的煤层埋深数据，我们解决了传统方法获取煤层埋深不准确的问题。石宝-邱家祠井田煤层气地质评价，通过上述方法获取了煤层埋深数据，使计算的煤层气储层压力，估算煤层气含气量更为合理，煤层气地质评价更为可靠。

上述方法步骤，充分利用了Mapgis工作平台中高程自动赋值、数据网格化、属性导出等功能。但由于制图的不规范以及Mapgis系统的空间分析对数据量处理的局限性，使生成图件不能够反映埋深实际情况。在此，作者通过大量操作，总结了两点建议。第一点，对等高线赋高程值时，应对线文件中微短线进行删除处理，排除微短线未赋值给数据读取带来的错误；第二点，在离散点三角网格化处理时，对要进行运算的文件，设置适当的数据间距，以免数据量过多导致软件停止工作。上述问题是作者在长期煤层气地质评价工作中总结得出。由于作者水平有限，文章中错误难免，敬请读者批评指正！

[1] 四川省地矿局一一三地质队. 四川省古蔺县石宝矿段石宝-邱家祠井田煤层气地质评价报告[R]. 2012：4～5

[2] 四川省地矿局一一三地质队. .四川省古蔺县川南煤田古叙矿区石宝矿段煤炭资源详查报告[R]. 2008：67～69.

[3] 傅雪海，秦勇, 等. 煤层气地质学[M]. 徐州：中国矿业出版社，2007：69～70