摘" 要：常規(guī)三維地質(zhì)建模方法處理復(fù)雜地質(zhì)結(jié)構(gòu)和空間屬性插值時，受限于地層連續(xù)性的處理和插值精度的不足，導致模型精度不夠高。為此，提出基于Kriging插值[ 2] 的鉆孔數(shù)據(jù)三維地質(zhì)可視化建模方法研究。通過標準化處理鉆孔數(shù)據(jù)并實現(xiàn)地層統(tǒng)一編碼，為建模提供數(shù)據(jù)基礎(chǔ)；利用Kriging插值技術(shù)，結(jié)合地質(zhì)屬性的空間變異性，進行空間屬性插值；最后通過三維網(wǎng)格劃分技術(shù)，將插值結(jié)果映射到網(wǎng)格單元中，構(gòu)建三維地質(zhì)可視化模型。實驗結(jié)果表明，該方法能夠提高三維地質(zhì)模型的精度、更準確地反映地下空間的地質(zhì)特征。

關(guān)鍵詞：kriging插值" 鉆孔數(shù)據(jù)" 三維地質(zhì)" 可視化建模" 空間屬性插值

中圖分類號：P624

Research on 3D Geological Visualization Modeling Method of Drilling Data Based on Kriging Interpolation

SHI De

Anhui Construction Engineering Water Resources Development Investment Group Co.， Ltd.， Bengbu， Anhui Province， 233010 China

Abstract： When conventional 3D geological modeling methods handle complex geological structures and spatial attribute interpolation， they are limited by the processing of stratigraphic continuity and insufficient interpolation accuracy， resulting in insufficient model accuracy. Therefore， research on three-dimensional geological visualization modeling method for drilling data based on Kriging interpolation is proposed. By standardizing drilling data and achieving unified coding of geological layers， a data foundation is provided for modeling; It uses Kriging interpolation technique， combined with the spatial variability of geological attributes， perform spatial attribute interpolation; Finally， using 3D mesh partitioning techniques， the interpolation results are mapped to grid cells to construct a 3D geological visualization model. The experimental results show that this method can improve the accuracy of three-dimensional geological models and more accurately reflect the geological characteristics of underground spaces.

Key Words： Kriging interpolation; Drilling data; 3D geology; Visualization modeling; Spatial attribute interpolation

三維地質(zhì)建模不僅能夠直觀地展示地下空間的地質(zhì)結(jié)構(gòu)，還為礦產(chǎn)資源評估、地質(zhì)災(zāi)害預(yù)測、地下工程建設(shè)等提供了重要的決策支持。然而，地質(zhì)數(shù)據(jù)的復(fù)雜性給三維地質(zhì)建模帶來挑戰(zhàn)。因此，研究一種高效、準確的三維地質(zhì)可視化建模方法，對于推動地質(zhì)科學的發(fā)展和應(yīng)用有重要意義。

針對三維地質(zhì)建模的需求，研究者們提出一系列方法。例如，基于構(gòu)造恢復(fù)理論的三維地質(zhì)建模方法，通過斷裂矢量場建模，實現(xiàn)含復(fù)雜斷裂網(wǎng)絡(luò)的地質(zhì)模型構(gòu)建[1]。但在處理地層連續(xù)性和屬性插值方面存在局限性?；贕emPy的三維地質(zhì)建模方法利用開源軟件GemPy創(chuàng)建三維地質(zhì)模型，通過提取剖面離散點坐標并進行曲線擬合，再轉(zhuǎn)換為方向數(shù)據(jù)，實現(xiàn)對地質(zhì)結(jié)構(gòu)的隱式表達[2]。然而，這種方法在地質(zhì)屬性插值方面可能存在精度不足的問題?；贐IM技術(shù)的邊坡地質(zhì)建模方法，通過運用空間插值插件實現(xiàn)地質(zhì)模型的建模[3]。但建模過程相對煩瑣，對于大規(guī)模地質(zhì)數(shù)據(jù)的處理效率有待提升。

因此，本文提出基于克里金插值（Kriging Interpolation，Kriging）鉆孔數(shù)據(jù)三維地質(zhì)可視化建模方法研究。通過標準化處理鉆孔數(shù)據(jù)、實現(xiàn)地層統(tǒng)一編碼，利用Kriging插值進行空間屬性插值，再通過三維網(wǎng)格劃分構(gòu)建三維地質(zhì)模型，實現(xiàn)地質(zhì)結(jié)構(gòu)的可視化表達，為地質(zhì)勘探、資源開發(fā)和地下工程建設(shè)提供有力的技術(shù)支持。

1" 基于kriging插值的鉆孔數(shù)據(jù)三維地質(zhì)可視化建模方法設(shè)計

1.1" 標準化鉆孔數(shù)據(jù)并實現(xiàn)地層統(tǒng)一編碼

鉆孔數(shù)據(jù)是三維地質(zhì)建模的基礎(chǔ)，但原始鉆孔柱狀圖數(shù)據(jù)難直接用于建模。需處理這些數(shù)據(jù)，提取位置、深度、地層類型及厚度等關(guān)鍵信息，并整理成統(tǒng)一格式。同時，統(tǒng)計地層類型，建立關(guān)系矩陣[ 3] ，如公式（1）所示。

關(guān)系矩陣的維度為地層類型數(shù)量×地層類型數(shù)量，矩陣中的元素表示兩個地層類型的關(guān)系。鉆孔根據(jù)地層在縱向位置上的出現(xiàn)順序，計算地層間的關(guān)系值。關(guān)系值的計算基于地層間的三種位置關(guān)系，如公式（2）所示。

統(tǒng)計各個鉆孔內(nèi)地層關(guān)系值的總和，得到關(guān)系矩陣[ 5] 的綜合值。根據(jù)關(guān)系矩陣的綜合值，對地層關(guān)系進行判斷，確定地層之間的相對位置關(guān)系[5]

1.2" 利用kriging插值進行空間屬性插值

在完成鉆孔數(shù)據(jù)的標準化處理及地層統(tǒng)一編碼后，利用Kriging插值方法對地質(zhì)屬性的空間分布進行估計。選擇地層厚度作為插值變量，根據(jù)鉆孔數(shù)據(jù)的空間分布和地質(zhì)特征，選擇合適的變異函數(shù)來描述地質(zhì)屬性的空間變異性，具體公式如下。

式[ 7] （3）中，表示位置處的地址屬性值；表示空間距離（或稱為滯后距）；表示在距離內(nèi)的樣本點對數(shù)量。選擇理論變異函數(shù)模型擬合觀測數(shù)據(jù)的半變異函數(shù)，并確定模型的參數(shù)。將標準化處理后的鉆孔數(shù)據(jù)作為樣本點集，每個樣本點包含位置信息和對應(yīng)的地質(zhì)屬性信息。Kriging插值的核心是權(quán)重系數(shù)的計算，系數(shù)用于將樣本點的屬性值加權(quán)求和以估計未知點的屬性值。權(quán)重系數(shù)的計算依賴變異函數(shù)和樣本點的空間配置。對于未知點處的地質(zhì)屬性值進行估計，插值公式如下。

式[ 9] （4）中，表示樣本點的數(shù)量；表示第個樣本點的屬性值；表示對應(yīng)的權(quán)重系數(shù)，在獲得權(quán)重系數(shù)后，對未知點的地質(zhì)屬性值進行估計。

1.3" 通過三維網(wǎng)格劃分構(gòu)建三維地質(zhì)模型

利用三維網(wǎng)格劃分技術(shù)來構(gòu)建三維地質(zhì)模型，以實現(xiàn)地質(zhì)結(jié)構(gòu)的可視化表達。選擇規(guī)則的立方體網(wǎng)格作為網(wǎng)格劃分的基礎(chǔ)方式?？紤]地質(zhì)屬性的空間變異性、計算資源的限制以及建模精度的要求確定網(wǎng)格尺寸。確定網(wǎng)格尺寸后，在研究區(qū)內(nèi)生成由個立方體單元組成的三維網(wǎng)格，其中、和分別代表網(wǎng)格在、和方向上的劃分數(shù)量。對每個網(wǎng)格單元，其中心點坐標為，而該單元的地質(zhì)屬性信息則被標記為。將插值得到的地質(zhì)屬性信息映射到網(wǎng)格單元中，映射公式如下。

網(wǎng)格單元根據(jù)位置信息獲取最近樣本點的地質(zhì)屬性，若受多個樣本點影響，則綜合信息。為精確表達地質(zhì)結(jié)構(gòu)，引入地層界面作為約束。利用鉆井數(shù)據(jù)確定地層界面位置形態(tài)，結(jié)合三維網(wǎng)格劃分和地質(zhì)屬性映射，最終構(gòu)建出包含地層分布、厚度、巖性的三維地質(zhì)可視化模型。

2" 實驗分析

2.1" 實驗對象

選取某典型地質(zhì)研究區(qū)作為實驗對象。該區(qū)域具有豐富的鉆孔數(shù)據(jù)，能夠反映該地區(qū)的地質(zhì)結(jié)構(gòu)和地層特征。鉆孔數(shù)據(jù)包括鉆孔位置、深度、地層類型、地層厚度等關(guān)鍵信息。實驗將鉆孔數(shù)據(jù)，利用Kriging插值和三維網(wǎng)格劃分技術(shù)來構(gòu)建三維地質(zhì)模型，并進行可視化表達。

2.2" 建模剖面分析

使用本文方法建立好的三維地質(zhì)模型進行剖面可視化，在三維地質(zhì)模型中，選擇特定的剖面進行切割，以展示模型內(nèi)部的地質(zhì)結(jié)構(gòu)和地層分布，如圖[ 11] 1所示。

2.3" 不同建模方法均方根誤差對比

為評估三維地質(zhì)模型的精度，選擇均方根誤差（Root Mean Square Error， RMSE）作為檢測模型精度的指標。實驗利用鉆孔數(shù)據(jù)中的地層厚度作為實際值，將模型預(yù)測的地層厚度作為預(yù)測值，計算均方根誤差來評估模型的精度。為進行對比分析，與3種建模方法進行對比，分別為方法1基于構(gòu)造恢復(fù)理論的三維地質(zhì)建模方法、方法2基于GemPy的三維地質(zhì)建模方法、方法3基于建筑信息模型[ 13] （uilding Information Modeling，BIM）技術(shù)的地質(zhì)建模方法。隨機選取[ 14] 5個鉆孔進行驗證，4種建模方法的均方根誤差結(jié)果如下表。

從表1數(shù)據(jù)看出，本文方法在[ 15] 5個鉆孔上的均方根誤差均小于其他3種方法。這表明該方法在預(yù)測地層厚度方面具有更高的精度和準確性。這一結(jié)果得益于本文方法利用Kriging插值進行空間屬性插值，選擇合適的變異函數(shù)來描述地質(zhì)屬性的空間變異性，并擬合觀測數(shù)據(jù)的半變異函數(shù)，從而更準確地估計未知點的地質(zhì)屬性值，提高三維地質(zhì)模型的精度。

3" 結(jié)語

本文提出的基于Kriging插值的鉆孔數(shù)據(jù)三維地質(zhì)可視化建模方法研究，通過標準化處理鉆孔數(shù)據(jù)、利用Kriging插值進行空間屬性估計與三維網(wǎng)格劃分技術(shù)，成功構(gòu)建了包含地層分布、厚度、巖性等信息的三維地質(zhì)可視化模型。實驗結(jié)果表明，該方法能夠準確地反映地質(zhì)結(jié)構(gòu)特征，提高三維地質(zhì)模型的精度。通過對比不同建模方法的均方根誤差，驗證本文方法的優(yōu)越性。未來，將繼續(xù)深入研究地質(zhì)建模技術(shù)，優(yōu)化算法流程，提高建模效率，為地質(zhì)學研究和工程應(yīng)用領(lǐng)域提供解決方案。

參考文獻

[1]花衛(wèi)華，曾新靈，郭丹陽，等.基于構(gòu)造恢復(fù)理論的含復(fù)雜斷層三維地質(zhì)建模方法[J].地球科學，2024，49（4）：1411-1420.

[2]趙勇，許國，盧鵬，等.基于GemPy的隱式三維地質(zhì)建模方法[J].人民長江，2023，54（10）：98-104.

[3]鄧云月，王章瓊，焦漫，等.基于BIM技術(shù)的邊坡地質(zhì)建模與有限元分析[J].土木建筑工程信息技術(shù)，2024，16（6）：122-126.

[4]于廣婷，劉同文，王琦，等.基于空間變異修正普通克里金的礦區(qū)復(fù)雜地形高精度建模[J].測繪通報，2022（11）：106-111.

[5]張琦，馬艷寧，王肖駿，等.基于優(yōu)化克里金插值方法的隧道主結(jié)構(gòu)面建模[J].土木工程學報，2022，55（S2）：74-82.