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