摘 要：隨著計算機建模技術的發(fā)展以及硬件運算水平的提高，計算機建模技術有了顯著的提高，三維地質(zhì)建模技術越來越成熟，在建模過程中已經(jīng)引入了地質(zhì)學概念。本文以SKUA-GOCAD三維地質(zhì)建模軟件為研究平臺，介紹了數(shù)據(jù)整理、結(jié)構(gòu)模型建設、屬性模型建設的總體思路，重點對比了基于地層約束的三維屬性模型和區(qū)域整體三維屬性模型的關聯(lián)和區(qū)別，從插值參數(shù)角度分析了產(chǎn)生這一現(xiàn)象的原因。結(jié)果表明在已完成三維地層結(jié)構(gòu)模型的基礎上，基于地層約束建設屬性模型，是一種提高屬性模型精度的成熟有效的方法。

　　關鍵詞：三維建模;屬性模型;分層插值;地層;鉆孔數(shù)據(jù)

　　0 前言

　　隨著計算機建模技術的發(fā)展以及硬件運算水平的提高，三維地質(zhì)建模技術越來越成熟(王亞靜，2016)。依據(jù)建模數(shù)據(jù)源的不同，三維地質(zhì)建模方法主要有以下5種：①鉆孔數(shù)據(jù)建模：直接利用鉆孔分層資料建立基礎數(shù)據(jù)庫，由建模軟件主導完成模型建設;②剖面數(shù)據(jù)建模：對二維的剖面數(shù)據(jù)進行“立體化”，利用剖面格網(wǎng)通過人機交互或建模軟件自動完成三維地質(zhì)模型的建設;③三維地震數(shù)據(jù)建模：以地震資料解譯技術為支撐，采用井震結(jié)合及相控的方法建立三維地質(zhì)模型(張洋洋等，2013);④屬性數(shù)據(jù)建模：在三維空間中建立格網(wǎng)，利用離散的屬性點作為插值點，按照特定的空間插值算法進行空間插值，反映的是地質(zhì)體內(nèi)部的屬性狀態(tài)情況;⑤多源數(shù)據(jù)建模：利用平面圖、鉆孔數(shù)據(jù)、剖面數(shù)據(jù)、地球物理或地球化學數(shù)據(jù)、屬性數(shù)據(jù)等開展模型建設，數(shù)據(jù)互補又互為約束對象，共同協(xié)調(diào)統(tǒng)一，提高了模型的精度。

　　在三維地質(zhì)模型建設及應用方面，國內(nèi)油氣勘探領域的應用程度最高，有色貴金屬礦山次之，城市地質(zhì)也已經(jīng)起步(李青元等，2016)。在油氣勘探領域，以油藏數(shù)值模擬為例，采用相控條件下的序貫高斯(SGS)算法，建立孔隙度、滲透率、氣飽和度模型，定量描述儲層參數(shù)的空間非均質(zhì)性(王威，2013);在金屬礦山領域，以銅礦礦床分析為例，建立不同類型的結(jié)構(gòu)模型，并對銅品位實施克立格插值(荊永濱，2010);在城市地質(zhì)領域，以北京市地鐵9號線某段為例，將同一土層內(nèi)賦予參數(shù)實驗值并進行插值，形成參數(shù)模型(董慧超，2011);在其他領域，以土壤屬性三維模型建設為例，利用地統(tǒng)計學插值方法，以10cm 為采樣梯度獲取土壤樣本，實現(xiàn)了農(nóng)田尺度土壤有機質(zhì)、pH、鈣、鋅的三維模型構(gòu)建和空間分異規(guī)律分析(周四維，2018)。

　　本文在前人研究基礎上，以當前城市地下空間資源開發(fā)利用主要深度(地下50m)以內(nèi)的地層為研究對象，選擇北京市平原區(qū)東部約60km2范圍作為研究區(qū)，開展三維結(jié)構(gòu)模型建設，并以結(jié)構(gòu)模型各個地層為約束體，逐層進行屬性插值，建設三維地質(zhì)屬性模型，更為精細的刻畫了該區(qū)域地層屬性的分布形態(tài)。

　　1 數(shù)據(jù)標準化

　　本次研究工作采用SKUA-GOCAD建模軟件輔助進行三維地質(zhì)模型建設，需按其數(shù)據(jù)格式要求進行數(shù)據(jù)標準化處理。

　　1.1 地層數(shù)據(jù)標準化

　　研究區(qū)內(nèi)共有15個鉆孔，深度均為50m左右，結(jié)合以往該地區(qū)工程層剖面資料，對鉆孔數(shù)據(jù)進行標準化。

　　參考以往該地區(qū)工程層剖面資料，與本次工作采用的15個鉆孔進行地層比對校準，在以往大地層韻律背景下依據(jù)15個鉆孔的實際分層情況進行校準，形成本次研究工作的標準分層，標準地層自上而下依次為：填土層(10)、粘質(zhì)粉土層(22)、粉質(zhì)粘土層(23)、粉細砂層(24)、粘質(zhì)粉土層(32)、粉質(zhì)粘土層(33)、粉細砂層(34)、粉質(zhì)粘土層(43)以及粉細砂層(44)，標準地層柱狀。

　　按照建模軟件數(shù)據(jù)格式要求，將15個鉆孔的分層數(shù)據(jù)拆分，形成鉆孔分層數(shù)據(jù)表。每個鉆孔分層數(shù)據(jù)表內(nèi)容均有15條，由該層對應的鉆孔分層數(shù)據(jù)構(gòu)成，數(shù)據(jù)表內(nèi)容包括鉆孔編號、地層代號、層頂深度、坐標。

　　1.2 屬性數(shù)據(jù)標準化

　　對15個鉆孔的室內(nèi)土工實驗表進行整理，按照標準化后的地層深度重新分段。如存在土工實驗取樣深度與地層深度交叉的情況，按照如下方式進行處理：根據(jù)該地層的頂?shù)装迳疃葘耐凉嶒灁?shù)據(jù)進行分段提取，將提取后的土工實驗數(shù)據(jù)按照厚度計算中心點坐標Z值(圖2)。

　　完成屬性數(shù)據(jù)分段后，重新整理數(shù)據(jù)表，按照地層順序?qū)?5個鉆孔所包含的屬性數(shù)據(jù)整理入庫，數(shù)據(jù)表內(nèi)容為鉆孔編號、參數(shù)值、坐標。

　　至此，完成建模所用鉆孔分層數(shù)據(jù)及屬性數(shù)據(jù)的整理、錄入工作。

　　2 模型建設

　　2.1 三維地層結(jié)構(gòu)模型建設

　　SKUA-GOCAD建模軟件基于工作流(Workflows)思路進行三維模型建設，地層結(jié)構(gòu)模型建設過程與人的思維保持一致。總體過程大致如下：數(shù)據(jù)錄入、地層整理、模型生成、模型調(diào)整完善。

　　在數(shù)據(jù)錄入階段，主要解決的是建模工作中相關數(shù)據(jù)的“地質(zhì)定義”，例如定義模型邊界線、定義斷層性質(zhì)等，及所有導入到建模軟件中參與建模的數(shù)據(jù)，都需要定義其本身的地質(zhì)意義。

　　在地層整理階段，主要解決的問題是各個地層的沉積順序以及各個地層之間的接觸關系(圖3)。在該階段還可以在各地層數(shù)據(jù)內(nèi)增加或減少控制點，調(diào)整各個地層的顏色等。

　　在模型生成階段，可以調(diào)整模型的精細度以及平滑程度，還可以通過軟件自帶的地層檢查功能，檢查每個參與建模的鉆孔地層一致性(圖4)，保證模型的數(shù)據(jù)準確性。

　　在模型調(diào)整完善階段，可以引入控制鉆孔與模型地層比對，發(fā)現(xiàn)誤差較大的位置可以利用控制鉆孔進行修正。因該軟件基于工作流模式，可以便捷的退回地層整理步驟直接在誤差地層上以鉆孔地層為依據(jù)增加控制點，重新執(zhí)行一次建模工作流即可完成模型調(diào)整完善。最終形成的三維地層結(jié)構(gòu)模型見圖5。

　　2.2 三維地質(zhì)屬性模型建設

　　SKUA-GOCAD建模軟件可以在三維地層結(jié)構(gòu)模型的基礎上，以地層為約束單獨在目標層內(nèi)進行屬性插值，建設三維地質(zhì)屬性模型。總體步驟可概括為：設置結(jié)構(gòu)模型網(wǎng)格大小，數(shù)據(jù)分析，變異函數(shù)求取，以及屬性插值。

　　(1)網(wǎng)格設置

　　結(jié)構(gòu)模型創(chuàng)建完成后，根據(jù)不同層位建模精度的要求，選擇性的設置不同網(wǎng)格密度。

　　以23層粉質(zhì)粘土為例：平面上，綜合工區(qū)范圍，將網(wǎng)格設置為20m，與整體保持一致;垂向上，該層厚度平均約4.5m，將網(wǎng)格設置為0.5m(可據(jù)實際需求調(diào)整)。據(jù)此，該層所包含的全部網(wǎng)格數(shù)為430×491×9=1900170個(圖6)。

　　(2)數(shù)據(jù)分析

　　數(shù)據(jù)分析階段，主要對原始屬性數(shù)據(jù)進行整理，利用人-機交互，按照數(shù)學方法對屬性數(shù)據(jù)進行優(yōu)化，根據(jù)地層的地質(zhì)學特點剔除屬性數(shù)據(jù)中的噪點，基于地質(zhì)背景分析屬性數(shù)據(jù)在目標層位的總體分布規(guī)律。若數(shù)據(jù)量大且分區(qū)復雜，可通過skua gocad平臺的data and trend analysis板塊實現(xiàn)分析(圖7)，獲取相應區(qū)塊(region)的趨勢，作為后續(xù)插值的約束條件;若數(shù)據(jù)量較少，未劃分或僅有少數(shù)區(qū)塊，可通過軟件的計算器等工具欄完成數(shù)據(jù)分析(圖8)。

　　(3)變異函數(shù)求取

　　以每個地層(已完成網(wǎng)格化處理)為約束條件，創(chuàng)建各地質(zhì)體獨立塊體(region)，并將其作為變異函數(shù)求取的約束條件(圖9)。因此，對同一屬性，每一層位對應不同的變異函數(shù)。需要注意的是，變異函數(shù)的求取并非一蹴而就，需結(jié)合地質(zhì)背景以及模型創(chuàng)建經(jīng)驗，不斷試驗，直至獲取最符合地質(zhì)認知的結(jié)果，插值方法采用克里格法。圖10為23層粉質(zhì)粘土的變異函數(shù)獲取示例。

　　(4)屬性插值

　　通過Reservoir Properties模塊，選擇合適的插值方法與已完成的變異函數(shù)，逐步完成屬性賦值。對比不同變異函數(shù)的插值結(jié)果，基于地質(zhì)認知選擇最優(yōu)者。處理各屬性區(qū)域的顏色、數(shù)據(jù)顯示方式等后，最終完成三維地質(zhì)屬性模型。圖11為本次創(chuàng)建的含水量屬性模型。

　　3 對比分析及應用

　　在屬性建模方面，較為傳統(tǒng)的是區(qū)域整體屬性建模方法。該方法不考慮建模工區(qū)的地層分布情況，沒有在工區(qū)地層內(nèi)逐層約束插值建設屬性模型，而是利用建模工區(qū)的整體厚度進行一次整體插值形成三維屬性模型。這種模型通常是為了展示區(qū)域內(nèi)地質(zhì)屬性整體情況，從宏觀尺度進行地質(zhì)分析應用，反應的是區(qū)域內(nèi)地質(zhì)屬性的分布規(guī)律;而本文提到的基于地層約束的屬性建模，充分考慮地層內(nèi)部地質(zhì)條件特征，以地層為約束，更能反應研究區(qū)范圍內(nèi)不同地層的屬性差異。相比于區(qū)域整體屬性模型，基于地層約束的屬性模型更為精細、更能反應地層實際屬性狀態(tài)，在模型分析應用方面也更加靈活和便捷。

　　為了驗證上述觀點，同時研究地層約束前提下三維地質(zhì)屬性建模與傳統(tǒng)三維地質(zhì)屬性建模數(shù)據(jù)成果的區(qū)別和關聯(lián)，本次工作按照工區(qū)范圍和深度，采用區(qū)域整體屬性建模方法創(chuàng)建了一個區(qū)域整體三維地質(zhì)屬性模型用于比對分析(圖12)。

　　通過二者比對，可證實存在以下關聯(lián)和區(qū)別：二者在屬性空間分布趨勢上基本一致(例如各個屬性值段分布的空間位置基本一致)，但存在屬性空間區(qū)域的范圍不同。產(chǎn)生這一結(jié)果的原因，主要是由于分層屬性插值建模方法逐層設置插值參數(shù)，更符合每個地層的特點，而整體插值方法只針對了數(shù)據(jù)本身制定插值算法，未考慮每個地層的實際情況。以23層粉質(zhì)粘土為例，單地層約束的插值結(jié)果參見圖13，與整體插值在該地層截取的結(jié)果存在差異(圖14)，但在反映層內(nèi)屬性數(shù)據(jù)的分布情況方面，更具代表性。

　　由于存在上述區(qū)別，基于地層約束的屬性模型在應用層面不可用區(qū)域整體模型替代。基于地層約束的屬性模型已經(jīng)按照地層進行了劃分，可以針對某一類地層開展模型應用研究，而區(qū)域整體模型無法指定某一類地層為對象開展研究工作;另外，利用地層約束的屬性模型統(tǒng)計功能，可統(tǒng)計每個地層內(nèi)屬性的分布情況，輔助相關工作的開展。以篩選出地層含水率小于30%的地區(qū)作為初步選定研究區(qū)域為例：若利用區(qū)域整體屬性模型進行分析，只能將區(qū)域內(nèi)所有含水率小于30%的地區(qū)篩選出來(圖15);而利用地層約束的屬性模型能將每一層中含水率小于30%的區(qū)域篩選出來，同時可以分析每一層中篩選出的區(qū)域在整體篩選結(jié)果中的占比(圖16)，結(jié)合地層巖性數(shù)據(jù)綜合分析，有針對性提出的開發(fā)利用規(guī)劃建議。

　　4 結(jié)論

　　三維地質(zhì)建模技術發(fā)展迅速，已經(jīng)有不少建模軟件從單一數(shù)學算法建模提升到地質(zhì)學概念結(jié)合數(shù)學算法建模的綜合建模理念，且對建模流程進行了優(yōu)化，提高了模型修改、完善的效率，本文所采用的SKUA-GOCAD三維地質(zhì)建模軟件就具備了以上特點，便于地質(zhì)人員理解和使用。

　　基于地層約束下的三維地質(zhì)屬性模型建設，從屬性數(shù)據(jù)分析、屬性插值等方面，充分考慮了不同地層的地質(zhì)屬性數(shù)據(jù)特點，與傳統(tǒng)的整體屬性建模方法相比，在每個地層內(nèi)分別對屬性數(shù)據(jù)進行了插值參數(shù)定義，提高了三維地質(zhì)屬性模型的精度，擴展了模型的應用層次。

　　受研究深度所限，本次并未涉及模型分析相關研究，將會在后續(xù)工作中根據(jù)實際需要，逐步開展相關方面研究工作，提升模型實際應用效果。

　　參考文獻

　　北京市地質(zhì)調(diào)查研究院， 2017. 北京市地下空間資源調(diào)查評價及關鍵技術研究成果報告[R].

　　董慧超， 2011. 北京地鐵九號線三維地質(zhì)建模研究[D]. 中國地質(zhì)大學(北京).

　　荊永濱， 2010. 礦床三維地質(zhì)混合建模與屬性插值技術的研究及應用[D]. 中南大學.

　　推薦閱讀：地質(zhì)專業(yè)核心期刊有哪些