芮學來

(淮南礦業(yè)集團 地質(zhì)勘探工程處物探隊，安徽 淮南 232052）

0 引言

利用數(shù)字測井資料計算巖石抗壓強度是煤田測井所具有的一個基本功能。但現(xiàn)有測井處理軟件如CLGIS、LOGANALYST等均只能提供抗壓強度成果曲線，不能提供地質(zhì)部門需要的分層抗壓強度成果表。該表提供了各分層的巖石名稱、深度、厚度及抗壓強度的最小值、平均值和最大值。以往這項工作只能完全由人工來完成，其做法是：先將測井抗壓強度曲線打印出來，然后用鉛筆、直尺分別量出每層抗壓強度的最小值，最大值，并根據(jù)曲線中心線大致估算出該層的平均值。最后統(tǒng)計成表，用WORD或EXCELL打印出來。因此這是一項耗時耗力的工作，不僅效率極低，而且計算結果不準確，迫切需要一種高效的方法取代之。

為此，我們通過對測井抗壓強度數(shù)據(jù)文件和測井數(shù)據(jù)庫柱狀文件格式的解析，用VB6.0語言開發(fā)了這個抗壓強度成果格式化輸出程序。該程序在輸入抗壓強度文件名、柱狀文件名及需要輸出的深度范圍，即可自動計算并輸出各分層的巖石名稱、深度、厚度、抗壓強度最小值、平均值及最大值等。該程序不僅操作方便，容易使用，而且大大提高了工作效率。

1 原始數(shù)據(jù)的獲取

我們現(xiàn)在使用的測井處理軟件有LOGSYS、CLGIS和LOGANALYST。它們對單孔數(shù)字測井資料的處理流程大致相同。即：1）建立測井數(shù)據(jù)庫，根據(jù)設置不同，庫內(nèi)可存放10條、20條、30條或更多條測井曲線；2）將現(xiàn)場采集到的測井數(shù)據(jù)讀入數(shù)據(jù)庫，不同測井參數(shù)分別以相應測井曲線的形式存放；3）對測井數(shù)據(jù)進行預處理；4）選用至少三種測井參數(shù)對煤巖層進行分層定厚，并生成測井柱狀文件；5）根據(jù)設計要求計算相關物性參數(shù)如巖石抗壓強度等，計算結果也以相應參數(shù)曲線的形式存放于數(shù)據(jù)庫中；6）輸出相關成果圖、表、曲線等。

為自動計算并輸出巖石抗壓強度成果表，我們需要巖石抗壓強度原始數(shù)據(jù)和鉆孔柱狀文件。在上述測井資料處理過程中所得抗壓強度參數(shù)是以抗壓強度曲線的形式存放于測井數(shù)據(jù)庫中，該數(shù)據(jù)不能被直接引用。但可通過編程或使用LOGANALYST軟件的數(shù)據(jù)導出功能將其從數(shù)據(jù)庫中導出，形成一個文本文件。該文件由11行文件頭及隨后由深度和抗壓強度數(shù)據(jù)組成的文件體，抗壓強度數(shù)據(jù)以5cm為采樣間隔，1米20個點。如下所示：

井號,ZJFJ

起始深度,330.00

終止深度,965.10

采樣間隔,0.050

空數(shù)據(jù)表示,

曲線數(shù),1

數(shù)據(jù)點數(shù),12703

曲線1代號,KYQD

曲線1名稱,

曲線1單位,

深度,KYQD

330.00 ,171.326920

330.05 ,173.547119

330.10 ,173.522171

330.15 ,173.459915

330.20 ,173.385437

330.25 ,175.607376

……

測井柱狀文件存放測井分層定厚結果，它以數(shù)據(jù)庫名稱重復為名，以GOL為后綴，存放于數(shù)據(jù)庫目錄內(nèi)。該文件內(nèi)容可以被直接使用。文件內(nèi)容與格式如下所示：

抗壓強度數(shù)據(jù)文件是我們進行分層抗壓強度計算的源數(shù)據(jù)，而具體怎么分層，每層的名稱、每層有多厚，則要用到測井柱狀文件中的數(shù)據(jù)。二才缺一不可。

2 程序的編制

根據(jù)程序功能要求及輸入輸出關系，設計程序框圖。然后用Visual Basic 6.0強大的編程功能，對一些常用的功能如對話框、提示框等，進行程序的編寫，由于篇幅關系，程序代碼此略。當代碼調(diào)試完畢，用VB6.0“文件”菜單中的“生成程序”功能，生成一個可獨立執(zhí)行程序，取名“分層抗壓強度格式化輸出程序”。以后即可隨時運行。

3 程序的應用

本程序使用簡單方便，具體使用步驟如下：

1）對目的層段進行分層定厚，生成測井柱狀文件XXXXXXXX.GOL，XXXX為數(shù)據(jù)庫文件名；

2）利用測井處理軟件計算并生成抗壓強度源數(shù)據(jù)，然后使用LOGANALYST軟件的“文件導出”功能，將源數(shù)據(jù)以文本格式導出，取名kyqd.txt。

為方便操作，我們一般將上述兩個文件置于當前數(shù)據(jù)庫目錄。

3）運行程序，界面如圖1所示：

圖1

分別輸入鉆孔名稱、開始深度、結束深度、抗壓強度文件名及路徑、GOL文件名及路徑。程序默認輸出對象是打印機。也可輸出到文件，或同時輸出到兩個對象。當點選輸出到文件時，程序會提示輸入存放結果的文件名，路徑默認與抗壓強度數(shù)據(jù)源文件相同。

4）單擊“打印成果表”，程序自動運行。打印輸出如下圖。對于輸出到文件內(nèi)容，則可用Excell表格軟件進行隨意編輯和排版，更方便地質(zhì)部門使用。

表1 張集礦中央?yún)^(qū)新副檢孔分層抗壓強度計算成果表

4 結語

單孔抗壓強度數(shù)據(jù)以1000米深鉆孔為例，5厘米采樣間隔，1米20個點，1000米共有20000個數(shù)據(jù)。這么大的數(shù)據(jù)量用人工去分層計算抗壓強度平均值、最小值、最大值，幾乎是不可能。所以只有采取如圖1所示那樣大致估算的方法，即使如此，僅一個鉆孔也需要一個人2至3天的時間才能做完。而本程序可以根據(jù)測井柱狀圖中的分層結果，自動從抗壓強度源文件中取數(shù)、計算并輸出成果表，只需幾分鐘。若不打印只輸出到文件，則瞬間完成，極大提高了工作效率。

[1]尹乾.中文Visual Basic 6.0實用教程[M].北京希望電子出版社.