陳 麒 玉，劉 剛,2*，吳 沖 龍,2，李 新 川，張 志 庭

(1.中國地質(zhì)大學(xué)(武漢)計算機(jī)學(xué)院，湖北 武漢430074；2.智能地學(xué)信息處理湖北省重點實驗室，湖北 武漢 430074)

城市地質(zhì)調(diào)查中知識驅(qū)動的多尺度三維地質(zhì)體模型構(gòu)建方法

陳 麒 玉1，劉 剛1,2*，吳 沖 龍1,2，李 新 川1，張 志 庭1

(1.中國地質(zhì)大學(xué)(武漢)計算機(jī)學(xué)院，湖北 武漢430074；2.智能地學(xué)信息處理湖北省重點實驗室，湖北 武漢 430074)

在城市地質(zhì)三維模型構(gòu)建中，數(shù)據(jù)來源龐雜，且不同應(yīng)用領(lǐng)域?qū)δＰ途鹊囊蟛煌．?dāng)前地質(zhì)體三維模型構(gòu)建技術(shù)還無法很好地支持大范圍復(fù)雜地質(zhì)結(jié)構(gòu)模型的快速構(gòu)建，建模過程繁雜，并影響后續(xù)地質(zhì)分析工作的開展。為了解決這一問題，該文面向三維城市地質(zhì)調(diào)查工作，提出了一種知識驅(qū)動的多尺度三維地質(zhì)體模型快速構(gòu)建方法。以專家知識為指導(dǎo)，通過對研究區(qū)域進(jìn)行盆地分析和古沉積環(huán)境分析，建立了研究區(qū)域內(nèi)符合城市地質(zhì)要求的沉積地層分層標(biāo)準(zhǔn)，摸清了時代內(nèi)沉積相和沉積微相的空間分布。在此基礎(chǔ)上，融合多源多類地質(zhì)數(shù)據(jù)，運用Kriging插值、多約束條件Delaunay三角剖分與局部拓?fù)渲貥?gòu)技術(shù)，詳細(xì)探討了從鉆孔與剖面數(shù)據(jù)→地層格架面→三維地質(zhì)體模型的快速構(gòu)建過程。該方法已在中國地質(zhì)調(diào)查局試點項目“閩江口地區(qū)地質(zhì)環(huán)境調(diào)查”項目中得到應(yīng)用，并集成在QuantyView3D系列地礦應(yīng)用軟件中。

三維地質(zhì)體建模；知識驅(qū)動；構(gòu)造-地層格架；拓?fù)渲貥?gòu)

0 引言

目前，三維城市地質(zhì)調(diào)查已成為中國地質(zhì)調(diào)查局新時期地質(zhì)工作的重點之一。隨著地質(zhì)信息科學(xué)理論和技術(shù)方法的形成和發(fā)展，三維地質(zhì)建模技術(shù)已在數(shù)字礦山、煤田及石油勘察、水利工程、災(zāi)害地質(zhì)和城市地質(zhì)等領(lǐng)域取得廣泛、深入的應(yīng)用[1-5]。

已有的三維建模方法依據(jù)數(shù)據(jù)來源可將其分為基于鉆孔、剖面及多源數(shù)據(jù)混合的建模方法[6-11]。鉆孔數(shù)據(jù)在一定程度上能夠直接反映研究區(qū)域的地下構(gòu)造情況，但也僅限于有限的區(qū)域，而地質(zhì)體的實際分布是延續(xù)的。鉆孔之間存在的盲區(qū)，需要通過相互比較和分析來形成剖面[12]，才能夠較好地反映鉆孔之間地層和構(gòu)造的變化。采用多源數(shù)據(jù)進(jìn)行模型的構(gòu)建，雖然可以提高模型的精度，但卻因為資料的獲取途徑、數(shù)據(jù)標(biāo)準(zhǔn)的不同而不可避免地會產(chǎn)生一定的信息沖突，這便需要有更好的機(jī)制來協(xié)調(diào)、融合多源異構(gòu)地質(zhì)數(shù)據(jù)，使其更好地為后期的三維建模服務(wù)。另一個突出的問題是這些方法在一定程度上都停留在基于地層對比的地質(zhì)結(jié)構(gòu)建模，所涉及的地質(zhì)內(nèi)容僅僅是地層巖性的劃分，不包含研究區(qū)的地質(zhì)背景知識和沉積物的沉積相和沉積體系，造成所構(gòu)建模型與實際的地質(zhì)內(nèi)容脫節(jié)，難以用來進(jìn)行古沉積體系、工程地質(zhì)條件、地質(zhì)災(zāi)害孕育條件、含水層和隔水層空間分布等地質(zhì)分析。已有學(xué)者提出了基于知識規(guī)則對剖面數(shù)據(jù)進(jìn)行處理，再進(jìn)行三維模型構(gòu)建的方法[13-18]，但也沒有完全解決以上這一問題。

開展城市地下三維地質(zhì)建模的目的，是為進(jìn)行城區(qū)工程地質(zhì)條件評價、地面建筑選址和地基穩(wěn)定性評價、地下空間利用和地下工程建設(shè)、地面沉降與塌陷等地質(zhì)災(zāi)害預(yù)測預(yù)警、地下水污染源追蹤與后備水源評價等，提供易于實現(xiàn)各種空間分析和空間數(shù)據(jù)挖掘的三維可視化精細(xì)地質(zhì)模型。因此，其建模對象主要是第四紀(jì)沉積體系的結(jié)構(gòu)模型及其巖土體地層格架模型。

針對存在的問題和建模的目標(biāo)，本文面向三維城市地質(zhì)調(diào)查工作，提出了一種知識驅(qū)動的多尺度三維地質(zhì)體模型快速構(gòu)建方法。以鉆孔及基于鉆孔所建立的剖面為基本數(shù)據(jù)源，建立多源多類數(shù)據(jù)的融合方法。通過對研究區(qū)域進(jìn)行盆地分析和古沉積環(huán)境分析，摸清了時代內(nèi)沉積相和沉積微相的空間分布，從而使生成的模型既基于實際的鉆孔資料同時又具備了真正意義上的帶有研究區(qū)域地質(zhì)涵義。在此工作基礎(chǔ)上，詳細(xì)探討了從鉆孔與剖面數(shù)據(jù)→地層格架面→三維地質(zhì)體模型的快速構(gòu)建過程。

1 知識驅(qū)動的多尺度三維地質(zhì)體建模方法

1.1 多源多類地質(zhì)數(shù)據(jù)的融合與預(yù)處理

多源多類地質(zhì)數(shù)據(jù)融合主要研究如何對不同來源和類型的數(shù)據(jù)進(jìn)行加工與轉(zhuǎn)換，使其相互協(xié)同補(bǔ)充，得到更加充分的應(yīng)用，從而實現(xiàn)對同一地質(zhì)實體更加客觀、精確的表達(dá)與認(rèn)識。數(shù)據(jù)的豐富與準(zhǔn)確性是決定所建三維地質(zhì)模型是否精準(zhǔn)的重要因素，但很多時候，并不是缺少數(shù)據(jù)，而是不能將已有的多源、多類、異構(gòu)數(shù)據(jù)進(jìn)行充分的協(xié)調(diào)融合，使其更好地為模型的構(gòu)建服務(wù)。

本文主要以鉆孔、剖面和沉積環(huán)境區(qū)劃圖等為數(shù)據(jù)來源。數(shù)據(jù)融合的主要難點在于數(shù)據(jù)屬性的多樣性和外在表現(xiàn)形式的多樣性。從數(shù)據(jù)內(nèi)容上看，涉及的數(shù)據(jù)有鉆孔、剖面、平面圖等；從數(shù)據(jù)表現(xiàn)形式上，有點、線、面等多種特征數(shù)據(jù)；從數(shù)據(jù)來源看，有實測數(shù)據(jù)和經(jīng)過分析推測、包含地質(zhì)專家認(rèn)識的成果數(shù)據(jù)。這就要求必須用統(tǒng)一的編碼規(guī)則對這些數(shù)據(jù)進(jìn)行融合處理。

對鉆孔分層信息、剖面地層屬性和沉積相屬性描述進(jìn)行統(tǒng)一編碼。其編碼規(guī)則如圖1a所示，前兩位為時代代號，中間兩位為巖土類別大類(相)代號，最后兩位為巖土類別小類(微相)代號。這樣的6位編碼，便可將鉆孔、剖面的時代、分層屬性統(tǒng)一起來，方便根據(jù)具體研究需求快速構(gòu)建多尺度三維地質(zhì)模型(圖1b、圖1c)。如果要研究該區(qū)域大的地質(zhì)背景及時代成因，則只需按照時代代號，提取時代的分層信息即可；如需繼續(xù)細(xì)化，可依次根據(jù)后面幾位提取更加豐富的分層信息，構(gòu)建更加精細(xì)的三維地質(zhì)模型。對于沉積相分區(qū)圖，主要關(guān)注的是某一時代內(nèi)沉積相的空間劃分，所以需要按照以上編碼規(guī)則將每個分區(qū)的所屬時代、沉積相類型和沉積微相類型進(jìn)行編碼(圖1d)，這樣就可以根據(jù)設(shè)定自動提取特定沉積環(huán)境的分區(qū)邊界，用于同一時代內(nèi)地層的進(jìn)一步劃分。

按照以上編碼規(guī)則，地質(zhì)專家對全區(qū)第四紀(jì)地層及基巖風(fēng)化層進(jìn)行統(tǒng)一編碼。但在實際地質(zhì)建模過程中，并不是建模區(qū)域包含全部地層分層信息(圖1e)；再者，應(yīng)用需求不同，也可對地層進(jìn)行相應(yīng)的簡化歸并處理。所以需要根據(jù)應(yīng)用需求和區(qū)域內(nèi)鉆孔和剖面的對比分析，制定某次建模的地層序列定義表(圖1f)。之后將按照此表中的編碼及層序，進(jìn)行鉆孔信息、剖面分層界線和沉積相分區(qū)邊界線提取。

在鉆孔和剖面數(shù)據(jù)的提取與處理中，需特別注意地層缺失的情況，即尖滅、透鏡體等復(fù)雜地質(zhì)結(jié)構(gòu)。圖2a示例性說明了在地層缺失的情況下鉆孔地層分隔點的提取。假設(shè)總共分為C1、C2、C3三層，那么對應(yīng)的地層分割點就有p0、p1、p2、p3，當(dāng)C2層缺失時，分割點p1和p2便處于同一位置，當(dāng)C1、C2層都缺失時，分割點p0、p1和p2處于同一位置，同理可處理所有類似缺失地層的情況。如圖2b所示為剖面地層分割線的處理方法，從上到下地層依次為Ci-1、Ci和Ci+1，當(dāng)Ci層兩邊缺失時，為了保證缺失處地層層面重合，必須讓每一條地層分割線都從頭到尾，并且保證缺失處點的一致性，即Ci與Ci+1層的上界線在缺失處的點p(i,1)、p(i,2)、p(i,3)、p(i,10)、p(i,11)、p(i,12)與p(i+1，1)、p(i+1，2)、p(i+1，3)、p(i+1，12)、p(i+1，13)、p(i+1，14)必須對應(yīng)重合。

圖1 多源多類地質(zhì)數(shù)據(jù)的融合處理

依據(jù)以上所述數(shù)據(jù)融合方法，進(jìn)行標(biāo)準(zhǔn)化編碼處理，按照試驗區(qū)三維地質(zhì)建模地層序列表對鉆孔、剖面分層信息及沉積環(huán)境區(qū)劃邊界信息進(jìn)行提取，最終以真實坐標(biāo)位置將其展布在三維空間中。與傳統(tǒng)將鉆孔投影在一個平面上的勘探剖面相比，其顯著特征是鉆孔及其連接而成的剖面形態(tài)嚴(yán)格遵循勘探工程在三維空間的展布形式，故而能夠更加真實地反映出地質(zhì)體在勘探位置處的特征。圖3是鉆孔、剖面地層界線及沉積環(huán)境區(qū)劃邊界在三維空間的真實展布情況。

圖2 地層缺失處理

圖3 鉆孔、剖面地層界線及沉積環(huán)境區(qū)劃邊界在三維空間的展布

1.2 帶約束條件的地層層面校正

在以上數(shù)據(jù)的基礎(chǔ)上，將鉆孔及剖面線離散成帶有高程屬性的空間離散點，利用Kriging空間插值算法進(jìn)行某一特定地層層面的擬合。最后將擬合點(插值得到的帶高程屬性點)和型值點(原始鉆孔及剖面線離散點)一起，運用Delaunay三角剖分算法生成地層TIN曲面。型值點也參與曲面剖分，可保證所建模型的地層信息在鉆孔及有勘察剖面經(jīng)過的位置與實際勘察數(shù)據(jù)完全一致(圖4a)。

由于輸入的型值點是有限的，在依次進(jìn)行上下層面構(gòu)建時，利用插值得到的擬合點，難免會出現(xiàn)上下交錯的情況，導(dǎo)致上下地層層面交錯，這在地層實體模型構(gòu)建時是不允許的，所以需要對交錯層面進(jìn)行校正。如圖4所示為層面校正過程示意圖，可分為從上到下校正和從下到上校正。層面的構(gòu)建順序是由具體數(shù)據(jù)特征決定的，一般會選擇從上到下來構(gòu)建，因為地表有高精度的等高線數(shù)據(jù)，這樣以高精度層面為約束面校正低精度數(shù)據(jù)，才能保證模型的準(zhǔn)確性。

圖4 層面校正過程示意

型值點是勘察實測數(shù)據(jù)，必須保證所建模型與實測數(shù)據(jù)的一致性，所以型值點不參與層面校正。擬合點是通過空間插值計算得到的估計點，其分布如圖4a參考面所示，每一個層面的擬合點在XOY平面的投影是完全一致的，但通過空間插值之后其高程屬性Z各不相同。所以校正過程就是要保證上下層面任意位置擬合點的高程值不交錯。

1.3 地質(zhì)體結(jié)構(gòu)模型拓?fù)渲貥?gòu)

對經(jīng)過層面校正后的地質(zhì)體上下格架面，自動計算獲取側(cè)面，并圍合成體，生成對應(yīng)地層的地質(zhì)體模型。但是，經(jīng)上述層面校正后，上下層面在地層缺失位置完全重合，構(gòu)建地質(zhì)體模型之后，必須對重合位置進(jìn)行體拓?fù)渲貥?gòu)，才能使地層形態(tài)符合實際情況。

本文選擇實體模型(Solid Model)來表達(dá)地質(zhì)實體。在此數(shù)據(jù)結(jié)構(gòu)的基礎(chǔ)上，通過對每個三角面片法線方向的判斷，將其分為頂?shù)變刹糠?。然后依次遍歷查找每個三角面片上3個頂點坐標(biāo)之間的距離，當(dāng)兩個三角面片對應(yīng)的3個頂點之間距離都為零時，則兩個面片完全重合，將其記錄從其數(shù)據(jù)結(jié)構(gòu)的拓?fù)溆蛑袆h除。遍歷完成后，調(diào)整拓?fù)溆蛑忻娴膫€數(shù)，用調(diào)整后的拓?fù)浣Y(jié)構(gòu)重構(gòu)該體。圖5為體拓?fù)渲貥?gòu)的過程示意，經(jīng)過處理，將圖5a中間區(qū)域的三角面拓?fù)鋸捏w拓?fù)浣Y(jié)構(gòu)中刪除，得到重構(gòu)后的體；圖5b上下層面直接成體后，經(jīng)拓?fù)渲貥?gòu)，得到該地層的真實形態(tài)。

圖5 體拓?fù)渲貥?gòu)

2 技術(shù)實現(xiàn)

該方法已在三維地學(xué)信息系統(tǒng)軟件平臺QuantyView[19]上得到實現(xiàn)，圖6是知識驅(qū)動的多尺度三維地質(zhì)體建模方法的整體實現(xiàn)流程，詳細(xì)描述了從數(shù)據(jù)整理入庫，到多源多類數(shù)據(jù)融合，再到模型動態(tài)構(gòu)建的全過程。首先將國家及行業(yè)標(biāo)準(zhǔn)、地質(zhì)專家解釋的成果數(shù)據(jù)以及包含專家知識的區(qū)域地質(zhì)認(rèn)識等存入專家知識庫，并以此為標(biāo)準(zhǔn)指導(dǎo)空間數(shù)據(jù)的規(guī)范化與標(biāo)準(zhǔn)化入庫。在專家?guī)斓幕A(chǔ)上，定制試驗區(qū)“三維地質(zhì)建模地層序列表”，并以此為約束，融合鉆孔、剖面以及帶有專家認(rèn)識的綜合評價圖(如沉積環(huán)境區(qū)劃圖)等多源多類地質(zhì)數(shù)據(jù)，使其按照實際空間位置展布在三維空間中。最后結(jié)合本文的關(guān)鍵技術(shù)，通過空間插值擬合地層格架面，并通過層面校正與體拓?fù)渲貥?gòu)等技術(shù)方法，動態(tài)構(gòu)建地層格架模型。根據(jù)后續(xù)需要，可依次構(gòu)建地質(zhì)時代模型→沉積相模型→精細(xì)巖土體模型，使其更好地支持后續(xù)的地質(zhì)分析與應(yīng)用。

圖6 知識驅(qū)動的多尺度三維地質(zhì)體建模方法的整體實現(xiàn)流程

3 應(yīng)用實例

本文依托于中國地質(zhì)調(diào)查局試點項目“閩江口地區(qū)地質(zhì)環(huán)境調(diào)查”，基本數(shù)據(jù)源為該項目開展過程中對研究區(qū)域內(nèi)進(jìn)行工程、水文、地?zé)岬刭|(zhì)調(diào)查所獲得的各類資料，具體包括鉆探數(shù)據(jù)、測井?dāng)?shù)據(jù)及剖面數(shù)據(jù)等。

3.1 研究區(qū)沉積環(huán)境分析

閩江口是福州市及其輻射經(jīng)濟(jì)區(qū)的總稱，處于閩江下游近入海口，中生代火山巖及燕山期花崗巖構(gòu)成了福州盆地的基底以及盆地周邊的山地。福州市坐落于閩江河口區(qū)內(nèi)側(cè)的福州斷陷盆地第四紀(jì)沖積平原之上。自從中生代花崗巖侵入后，該盆地基底地殼長期處于上升剝蝕的狀態(tài)。直到晚更新世時期，盆地開始發(fā)生沉降，才開始充填并堆積第四紀(jì)沉積物。晚更新世以來，由于氣候的變化，海平面出現(xiàn)3次升降。海水的每次進(jìn)退都對盆地的沉積環(huán)境產(chǎn)生過重要的影響，在不同時期形成了不同的沉積相和沉積體系。正是因為以上原因，造成福州盆地第四紀(jì)沉積比較復(fù)雜，陸相沉積和海相沉積交錯出現(xiàn)，第四紀(jì)蓋層很薄，平均在30～40 m左右，同一性質(zhì)沉積物零碎且分散，這些都為巖性對比及后續(xù)的三維建模帶來了一定的困難。

依次對福州盆地各個時代做沉積環(huán)境分析，圈畫出不同時代沉積相及沉積微相的空間分布，矢量化提取沉積相邊界，并通過數(shù)據(jù)融合處理，使其與鉆孔及剖面一起參與到后續(xù)的三維地質(zhì)建模中。

3.2 建模實例

整個研究區(qū)面積為1 587 km2，布設(shè)勘查剖面50條，其中南北向45條，東西向5條，用到鉆孔1 374個，平均鉆孔間距為1 km?？紤]到僅依靠剖面數(shù)據(jù)無法保證三維模型的精度，又規(guī)范化處理了1 149個鉆孔，和50條剖面一起參與到三維地質(zhì)建模中。為了保證地表地形的精度，使用區(qū)域內(nèi)5 m等高線擬合生成地表面。

結(jié)合地質(zhì)專家認(rèn)識，并根據(jù)城市地質(zhì)實際需求，著重對第四紀(jì)沉積層進(jìn)行精細(xì)建模，依次構(gòu)建了地質(zhì)時代模型、沉積相模型和精細(xì)巖土體模型(圖7,彩圖見封2)。其中，圖7a是該區(qū)域地質(zhì)時代模型，從上到下依次為全新世上段(Qh3)、全新世中下段(Qh1-2)、更新世上段(Qp3)、更新世中段(Qp2)、更新世下段(Qp)和基巖。在時代模型的基礎(chǔ)上，對每一個時代體模型進(jìn)行沉積相劃分，構(gòu)建沉積相模型(圖7b)。從模型中可清楚地看出，由于海侵，在Qp3時期，研究區(qū)東部區(qū)域為海相沉積，主要為砂質(zhì)潮坪相和泥質(zhì)潮坪相。該模型還對古河道及河道砂、洪泛平原等都有清晰地刻畫。圖7c為更小區(qū)域的精細(xì)巖土體模型，按時代依次細(xì)化，總共17層，精細(xì)刻畫并可視化表達(dá)了區(qū)域內(nèi)各時代巖土體空間分布狀況。

圖7 試驗區(qū)建模實例

對不同尺度的區(qū)域模型，可通過設(shè)置剖分網(wǎng)格間距來控制數(shù)據(jù)點密度，從而控制模型精度。本次建模中，對于較大范圍的時代模型和沉積相模型，網(wǎng)格間距為500 m，對于更小范圍的精細(xì)巖土體模型，網(wǎng)格間距20 m。所建模型從不同尺度、不同細(xì)節(jié)層次刻畫了區(qū)域內(nèi)地下地質(zhì)構(gòu)造、沉積環(huán)境分布情況。尤其是第四紀(jì)沉積地層的精細(xì)刻畫，分層數(shù)較多，層與層之間形態(tài)各異，且同一層地層存在大量缺失，這都在一定程度上反映了該方法良好的建模質(zhì)量和效果，驗證了該方法的準(zhǔn)確性和穩(wěn)健性。

4 結(jié)論

本文面向三維城市地質(zhì)調(diào)查工作，提出了一種知識驅(qū)動的多尺度三維地質(zhì)體模型快速構(gòu)建方法。所選研究區(qū)為中國福建省閩江口地區(qū)，結(jié)合地質(zhì)工作人員的認(rèn)識，繪制區(qū)域內(nèi)50條地質(zhì)剖面，和鉆孔數(shù)據(jù)一起參與到模型的構(gòu)建中。同時，為了解決城市三維地質(zhì)模型缺乏地質(zhì)含義的問題，通過對研究區(qū)域進(jìn)行盆地分析和古沉積環(huán)境分析，摸清了時代內(nèi)沉積相和沉積微相的空間分布，從而使生成的模型既基于實際的鉆孔資料同時又具備真正意義上的帶有研究區(qū)域地質(zhì)涵義。建立了多源多類數(shù)據(jù)的融合方法，使鉆孔、剖面、沉積相區(qū)劃邊界等不同來源、不同類型的數(shù)據(jù)有效地融合在一起。在此工作基礎(chǔ)上，運用Kriging插值、多約束條件Delaunay三角剖分算法，擬合生成地層格架面，并通過層面校正和體局部拓?fù)渲貥?gòu)技術(shù)，生成符合要求的三維地質(zhì)體模型。根據(jù)后續(xù)地質(zhì)分析需要，動態(tài)生成多尺度構(gòu)造-地層格架模型：時代模型→沉積相模型→精細(xì)巖土體模型。該方法已集成于地學(xué)信息系統(tǒng)QuantyView3D平臺中，所建模型滿足城市地質(zhì)分析需求，進(jìn)一步驗證了模型的有效性和該方法的健壯性。包含斷層、褶皺及倒轉(zhuǎn)等復(fù)雜地質(zhì)結(jié)構(gòu)模型的動態(tài)構(gòu)建，以及地質(zhì)知識在這些特殊情況下的表達(dá)與融合方法等，需要做進(jìn)一步研究和討論。

[1] JACKSON I.OneGeology：Making geological map data for the earth accessible[J].Episodes,2007,30(1):60-61.

[2] 吳沖龍,毛小平,田宜平,等.三維數(shù)字盆地構(gòu)造-地層格架模擬技術(shù)[J].地質(zhì)科技情報,2006,25(4):1-8.

[3] 何珍文.地質(zhì)空間三維動態(tài)建模關(guān)鍵技術(shù)研究[D].武漢：華中科技大學(xué),2008.

[4] ZANCHI A,FRANCESCA S,STEFANO Z,et al.3D reconstruction of complex geological bodies:Examples from the alps[J].Computers and Geosciences,2009,35(1):49-69.

[5] 李章林,吳沖龍,張夏林,等.東勝煤田北部地段3D煤層模型的動態(tài)構(gòu)建及不確定性評估[J].地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報),2014,39(10):1398-1405.

[6] WU Q,XU H,ZOU X K.An effective method for 3D geological modeling with multi-source data integration[J].Computers and Geosciences,2005,31(1):35-43.

[7] MING J,PAN M,QU H G,et al.GSIS:A 3D geological multi-body modeling system from netty cross-sections with topology[J].Computers and Geosciences,2010,36(6):756-767.

[8] 朱良峰,潘信,吳信才,等.地質(zhì)斷層三維可視化模型的構(gòu)建方法與實現(xiàn)技術(shù)[J].軟件學(xué)報,2008,19(8):2004-2017.

[9] LEMON A M,JONES N L.Building solid models from boreholes and user-defined cross-sections[J].Computers and Geosciences,2003,29(5):547-555.

[10] HE Z W,WU C L,TIAN Y P,et al.Three-dimensional reconstruction of geological solids based on section topology reasoning[J].Geo-spatial Information Science,2008,11(3):201-208.

[11] 屈紅剛,潘懋,明鏡,等.基于交叉折剖面的高精度三維地質(zhì)模型快速構(gòu)建方法研究[J].北京大學(xué)學(xué)報(自然科學(xué)版),2008,44(6):915-920.

[12] 毛善君.灰色地理信息系統(tǒng)——動態(tài)修正地質(zhì)空間數(shù)據(jù)的理論和技術(shù)[J].北京大學(xué)學(xué)報(自然科學(xué)版),2002,38(4):556-562.

[13] KESSLER H,MATHERS S,SOBISCH H G.The capture and dissemination of integrated 3D geospatial knowledge at the British Geological Survey using GSI3D software and methodology[J].Computers and Geosciences,2009,35(6):1311-1321.

[14] PERRIN M,ZHU B T,RAINAUD J F,et al.Knowledge-driven applications for geological modeling[J].Journal of Petroleum Science and Engineering,2005,47(1):89-104.

[15] ROYSE K R.Combining numerical and cognitive 3D modelling approaches in order to determine the structure of the Chalk in the London Basin[J].Computers and Geosciences,2010,36(4):500-511.

[16] SCHARLING P B,RASMUSSEN E S,SONNENBORG T O,et al.Three-dimensional regional-scale hydrostratigraphic modeling based on sequence stratigraphic methods:A case study of the Miocene succession in Denmark[J].Hydrogeology Journal,2009,17(8):1913-1933.

[17] WYCISK P,HUBERT T,GOSSEL W,et al.High-resolution 3D spatial modelling of complex geological structures for an environmental risk assessment of abundant mining and industrial megasites[J].Computers and Geosciences,2009,35(1):165-182.

[18] 吳立新,陳學(xué)習(xí),車德福,等.一種基于GIP的地下真3D集成表達(dá)的實體模型[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2007,32(4):331-335.

[19] WU C L,LIU G,TIAN Y P,et al.Geoview:A computer-aided system for informatization of geological and mineral resources survey and exploration works[C].Proceedings of IAMG 2005:GIS and Spatial Analysis.Toronto,Canada,2005.958-963.

Knowledge-Driven Multiple Scale 3D Geological Modeling Method in Urban Geological Survey

CHEN Qi-yu1,LIU Gang1,2,WU Chong-long1,2,LI Xin-chuan1,ZHANG Zhi-ting1

(1.SchoolofComputerScience,ChinaUniversityofGeosciences,Wuhan430074;2.HubeiKeyLaboratoryofIntelligentGeo-InformationProcessing,ChinaUniversityofGeosciences,Wuhan430074,China)

3D urban geological survey has been a major part of the work of China Geological Survey in recent years.In the process of 3D urban geological modeling,the data are enormous and complex while the different application area require different precision.At present，3D geological modeling technology cannot well support the rapid modeling of complex geological structure,leading to the intricate modeling process and influencing the further geological analysis.In order to solve this problem,this paper proposes a rapid knowledge-driven multiple scale 3D geological modeling method according to 3D urban geological survey.Instructed by expert knowledge,the analysis of basin and sedimentary environment in the study area helps to build up the division standard of sedimentary stratum conforming to the urban geological requirement.Thereupon,a good knowledge of spatial distribution of sedimentary facies and sedimentary micro facies is obtained in the different geologic age.Based on the work above,the paper uses Kriging interpolation,multi-constraints Delaunay triangulation and surface topological reconstruction technology,to discuss the rapid modeling process adequately.This method has been applied in a pilot project "Investigation of Urban Geological Environment in Minjiang Port Area"which is launched by Geological Survey of China and it has been integrated in the QuantyView3D series application software in the geological and mineral field.These also show the availability and robustness of the proposed method.

3D geological modeling;knowledge-driven;tectonic-stratigraphic framework;topology reconstruction

2015-12-15；

2016-01-25

國家863計劃重點項目課題(2012AA121401)；國家自然科學(xué)基金面上項目(41172300)

陳麒玉(1990-)，男，博士研究生，主要從事三維地質(zhì)建模、時空數(shù)據(jù)模型和GIS應(yīng)用方面的研究和開發(fā)。*通訊作者E-mail:liugang67@163.com

10.3969/j.issn.1672-0504.2016.04.003

P642；TP391；P208

A

1672-0504(2016)04-0011-06