張園園 何欣 杜紅旺

摘 要：在對(duì)工程地質(zhì)鉆孔、地層分層、剖面等數(shù)據(jù)分析研究的基礎(chǔ)上，進(jìn)行數(shù)據(jù)量化、篩選、標(biāo)準(zhǔn)化后，采用Dream Rocks軟件創(chuàng)建完成了通州區(qū)（906 km2）地下60 m工程建設(shè)層的三維地質(zhì)結(jié)構(gòu)與屬性模型。模型可以達(dá)到的效果：1）通過三維方式展示砂土液化在空間上的分布;2）通過層間拉伸的方式，展示上下層或相鄰各層之間的三維空間關(guān)系;3）模型成果形式能轉(zhuǎn)換成標(biāo)準(zhǔn)的obj格式，能較好地融入其他相關(guān)三維系統(tǒng)中，開展深層次的計(jì)算研究。

關(guān)鍵詞：Dream Rocks;三維地質(zhì)建模;工程鉆孔;砂土液化

Abstract： Based on the quantification， screening and standardization of the data obtained from engineering geological drilling， stratification and profiles， the 3D geological structure and attribute model of the 60 m underground engineering construction layer in Tongzhou District （906 km2） was established by using Dream Rocks software. This model can display the 3-D spatial distribution of sand liquefaction， as well as the 3D spatial relationship between upper and lower layers or adjacent layers by stretching between layers. The model results can be converted into standard obj format to be better integrated into other relevant 3D systems and carry out in-depth computational research.

Keywords： Dream Rocks; 3D geological modeling; boreholes; sand liquefaction

地質(zhì)學(xué)經(jīng)常利用地質(zhì)平面圖、剖面圖、等值線圖等二維圖件，輔助展示并應(yīng)用于解決實(shí)際的地質(zhì)問題。在城市大型項(xiàng)目建設(shè)工作中，僅靠缺乏整體性、連續(xù)性的二維圖件，很難直觀且準(zhǔn)確地展示真實(shí)的地質(zhì)環(huán)境條件， 給城市地下空間的整體規(guī)劃、設(shè)計(jì)、建設(shè)以及后期維護(hù)工作帶來相當(dāng)大的難度。

加拿大工程地質(zhì)學(xué)家Houlding 在1994年首次提出地質(zhì)理論與計(jì)算機(jī)三維可視化技術(shù)有機(jī)融合的三維地質(zhì)建模概念。三維地質(zhì)模型將各地層的空間分布及物理力學(xué)參數(shù)性質(zhì)以三維形式直觀展示，對(duì)合理利用地下空間資源以及有效保障地下空間地質(zhì)安全建設(shè)具有重要意義（何靜等，2019）?；诳辈煦@孔數(shù)據(jù)分析構(gòu)建三維地質(zhì)實(shí)體模型一直是國(guó)內(nèi)外研究的熱點(diǎn)，經(jīng)過多年的研究發(fā)展，我國(guó)在理論、程序開發(fā)以及實(shí)例應(yīng)用等方面均取得了大量的科研成果。理論方面編纂并修訂發(fā)行CJJ/T 157-2010 《城市三維建模技術(shù)規(guī)范》、CH/T 9015-2012 《三維地理信息模型數(shù)據(jù)產(chǎn)品規(guī)范》、CH/T 9017-2012 《三維地理信息模型數(shù)據(jù)庫(kù)規(guī)范》等相關(guān)理論規(guī)范。實(shí)例應(yīng)用有二維剖面圖的三維實(shí)體重建（李獻(xiàn)峰等，2008年），精細(xì)化三維地質(zhì)建模在城市地下空間開發(fā)中的應(yīng)用（古銳開等，2019），也有以某區(qū)（通州、雄安、柳州等）為例探討城市三維地質(zhì)模型構(gòu)建（李靜等，2016;李敏等，2018;馬震等，2019;楊向敏等，2019）。

2016—2019年北京市地質(zhì)工程勘察院先后承擔(dān)了通州區(qū)重大地質(zhì)問題調(diào)查與評(píng)價(jià)子課題，對(duì)該區(qū)建筑抗震設(shè)計(jì)條件、古河道分布、砂土液化等工程地質(zhì)條件進(jìn)行小范圍的專門研究，并基于Dream Rocks軟件，利用地質(zhì)鉆孔、地層分層、剖面等數(shù)據(jù)，建立了局部工程層三維地質(zhì)結(jié)構(gòu)模型，在模型中展示的砂土液化判別為后續(xù)很多地質(zhì)工程研究提供了理論依據(jù)，后續(xù)的地質(zhì)勘察研究成果讓模型精確度進(jìn)一步得到驗(yàn)證。隨著建設(shè)項(xiàng)目的推進(jìn)，巖土工程勘察取得了大量建模區(qū)域巖土體的物理力學(xué)指標(biāo)，對(duì)建模區(qū)域內(nèi)細(xì)微和局部的工程地質(zhì)條件有了進(jìn)一步的深入了解。為更好地服務(wù)于研究區(qū)大范圍的城市建設(shè)，本文在現(xiàn)有資料基礎(chǔ)上，通過搜集全區(qū)的工程地質(zhì)、水文地質(zhì)和區(qū)域構(gòu)造等地質(zhì)資料，建立了工程建設(shè)層三維地質(zhì)模型，直觀反映建模區(qū)域內(nèi)工程地質(zhì)層的地質(zhì)概況，為開展具體的工程建設(shè)提供指導(dǎo)和參考。

1 工程建設(shè)層概況

研究區(qū)地表下60 m深度范圍分為人工堆積層、新近沉積層及一般第四紀(jì)沉積層3大類。其中人工堆積層主要為人工填土①層（圖1） ，填土厚度由東向西遞減。新近沉積②大層主要為粉砂—細(xì)砂②層、粉土②1層、粉質(zhì)黏土②2層（圖2），廣泛分布在通州皮村、徐辛莊—張家灣一線的溫榆河、小中河、潮白河、北運(yùn)河、通惠河等現(xiàn)今或近代河道、河漫灘及其兩岸一級(jí)階地;通州北部的宋莊雙埠頭—尹各莊—葛渠村一帶、西部的永順紀(jì)莊—梨園崗墳—馬駒橋北門口一帶及南部永樂店北辛店—大務(wù)一帶缺失;一般第四紀(jì)沖洪積層主要為細(xì)砂與粉質(zhì)黏土互層，其中③大層主要為粉砂—細(xì)砂層③層、粉土③1層、粉質(zhì)黏土③2層（圖3），主要分布在通州北部的宋莊雙埠頭—尹各莊—葛渠村一帶、西部的永順紀(jì)莊—梨園崗墳—臺(tái)湖次渠—馬駒橋北門口一帶及南部永樂店北辛店—大務(wù)一帶，其余地區(qū)缺失。④大層中細(xì)砂④層從西到東厚度逐漸變厚（在北馬莊附近厚度約1.6 m，而在東部縣附近厚度可達(dá)16 m），粉質(zhì)黏土④1層在南部地區(qū)分布較連續(xù);粉土④2層則局部分布，多以透鏡體形式存在;而圓礫、卵石④3僅局部有揭露（圖4）。⑤大層以粉質(zhì)黏土⑤層、粉土⑤1層為主，廣泛分布，形成區(qū)域性的隔水層;而細(xì)砂—中砂⑤2則局部分布（圖5），并多以透鏡體形式存在于粉土及黏質(zhì)粉土層之間。在潞城胡各莊—張家灣大甘棠一帶因缺失隔水層，局部形成了上、下層地下水的聯(lián)系天窗。⑥大層中細(xì)砂—中砂⑥層分布連續(xù)穩(wěn)定，粉質(zhì)黏土⑥1層、粉土⑥2層分布不連續(xù)，多呈透鏡體存在;而圓礫—卵石⑥3僅在東北處極少鉆孔內(nèi)揭露（圖6）。受鉆孔深度限制，僅在通州南部60 m深度鉆孔內(nèi)揭露本層，且未揭穿本大層（北京市地質(zhì)工程勘察院，2007）。

2 三維地質(zhì)模型構(gòu)建

2.1 三維建模軟件介紹

本文采用DreamRocks工程地質(zhì)三維自動(dòng)建模系統(tǒng)軟件進(jìn)行模型構(gòu)建，建模采用基于鉆孔數(shù)據(jù)的不規(guī)則格網(wǎng)層序地層構(gòu)建方法。該軟件具有以下特點(diǎn)：1）能夠較好的處理地層尖滅、正確地反映地層時(shí)序，保證不同地質(zhì)年代的地層在三維空間中覆蓋與被覆蓋的關(guān)系與解釋成果一致;2）能在已經(jīng)建立的工程層三維地質(zhì)模型的基礎(chǔ)上，能以三維方式展示砂土液化在空間上的分布;3）能通過層間拉伸的方式，反映上下層或相鄰各層之間在三維空間中的關(guān)系;4）能通過垂直剖切或水平剖切方式，反映上下層或相鄰各層之間在剖切面上的關(guān)系;5）將三維模型成果轉(zhuǎn)換成標(biāo)準(zhǔn)的obj格式，使三維模型成果融入相關(guān)三維系統(tǒng)中。

2.2 三維建模方法與路線

（1）三維建模方法及特點(diǎn)

Dream Rocks建模軟件支持三維地質(zhì)空間信息獲取、管理、分析和可視化。該系統(tǒng)可以使用多種類型的數(shù)據(jù)源生成三維數(shù)據(jù)，并提供多種三維數(shù)據(jù)的生成方法。建模方法具有以下特點(diǎn)：1）提供基于地質(zhì)要素的建模方法，并合理地將地質(zhì)概念模型和實(shí)際的工程數(shù)據(jù)相結(jié)合，便于地質(zhì)人員使用;2）建模方法在地質(zhì)人員認(rèn)識(shí)工作的基礎(chǔ)上，盡量做到自動(dòng)化，減輕建模的工作量;3）針對(duì)不同類型的地質(zhì)建模工作能夠提供專門的建模方法支持，提高建模效率;4）充分考慮構(gòu)造對(duì)于地層模型的影響，地層模型中需要體現(xiàn)斷層控制;5）提供合理的地質(zhì)模型更新方法，便于隨地質(zhì)工作的不斷深入逐步細(xì)化地質(zhì)模型，并為后續(xù)地質(zhì)工作的持續(xù)進(jìn)行打下基礎(chǔ);6）提供模型的質(zhì)量控制功能，對(duì)于模型的面網(wǎng)質(zhì)量及厚度等參數(shù)進(jìn)行控制，便于后續(xù)的模型應(yīng)用。

（2）三維建模技術(shù)路線

項(xiàng)目工程層三維地質(zhì)結(jié)構(gòu)建模采用基于多源數(shù)據(jù)和垂向多精度模型構(gòu)建策略。充分收集建模區(qū)域的地質(zhì)資料，包括：鉆孔、各類剖面、等值線、地質(zhì)圖、地質(zhì)報(bào)告等;根據(jù)收集的地質(zhì)資料以及各種建模數(shù)據(jù)進(jìn)行仔細(xì)研究整合確定建模縱向標(biāo)準(zhǔn)分層，對(duì)所有數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化統(tǒng)一，保障多個(gè)的統(tǒng)一性，同時(shí)形成地質(zhì)建模概念模型，以便開展三維建模工作;根據(jù)地質(zhì)概念模型重新審視各類數(shù)據(jù)，進(jìn)行校正和取舍;根據(jù)建模區(qū)域的地質(zhì)單元及地質(zhì)要素的地質(zhì)特征，選擇性的劃分子建模區(qū)域，同時(shí)選擇相應(yīng)的建模方法。建模方法包括層狀地層鉆孔建模、多源復(fù)雜構(gòu)造建模（圖7）。

2.3 數(shù)據(jù)收集與處理

（1）數(shù)據(jù)來源

三維建模以研究區(qū)內(nèi)多方位的數(shù)據(jù)為基礎(chǔ)，包括225個(gè)孔深大于60 m鉆孔基本信息、分層信息、巖性描述信息、鉆孔土工試驗(yàn)信息、61個(gè)剖面、工作區(qū)基礎(chǔ)地理信息以及砂土液化判別結(jié)果信息等。

（2）工程鉆孔數(shù)據(jù)

對(duì)收集的各種鉆孔數(shù)據(jù)進(jìn)行整理、儲(chǔ)存、加工和統(tǒng)計(jì)，按照建模軟件對(duì)數(shù)據(jù)格式的要求進(jìn)行整理，建立建模所需的鉆孔數(shù)據(jù)源。

（3）鉆孔數(shù)據(jù)格式標(biāo)準(zhǔn)化

收集到的鉆孔坐標(biāo)系統(tǒng)是從不同數(shù)據(jù)來源的工程坐標(biāo)系轉(zhuǎn)換到北京54投影坐標(biāo)系。用于三維建模的鉆孔數(shù)據(jù)根據(jù)建模軟件的數(shù)據(jù)需求，分別對(duì)鉆孔數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)分層，建立含7個(gè)地層大層、20個(gè)地層亞層的地層結(jié)構(gòu)模型（表1）：

三維地質(zhì)屬性模型：根據(jù)獲取到得實(shí)驗(yàn)數(shù)據(jù)，完成了含水率、天然密度、干密度、比重、孔隙比、飽和度、液限、塑限、塑性指數(shù)、液性指數(shù)、自重壓力、壓縮模量、承載力等屬性模型得構(gòu)建（表2）。

3 三維模型成果可視化空間展示

（1）巖土體結(jié)構(gòu)三維空間展示

根據(jù)收集獲取的225個(gè)工程地質(zhì)鉆孔基本信息、分層信息、巖性描述信息，本次工作構(gòu)建了通州區(qū)地基土三維結(jié)構(gòu)模型（圖8、圖9），覆蓋面積約906 km2。模型能直觀反映工作區(qū)內(nèi)工程地質(zhì)層巖土體的空間展布情況以及軟弱土、液化砂土、密實(shí)砂土等地質(zhì)體的分布情況，為工作區(qū)地下空間開發(fā)、規(guī)劃與建設(shè)的地質(zhì)環(huán)境條件評(píng)價(jià)提供了依據(jù)。

（2）三維地質(zhì)體與地表影響信息疊加

通過工程層三維地質(zhì)模型疊加地表遙感影像和DEM，更加直觀和完整地反映出了該地區(qū)地表河流水系、地形地貌、道路居民地與地下地層分布、走向間的空間關(guān)系（圖10）。

（3）人工填土層綜合展示

三維地質(zhì)模型能較好地反映出人工填土的空間分布特征（圖11），為深入開展工程地質(zhì)研究提供了數(shù)據(jù)支撐。工作區(qū)表層廣泛分布人工填土，總體來說人工填土厚度未超過6 m，鉆孔揭露填土厚度0～5.8 m。人工填土由于其成因、成分、粒度變化很大，因而其密度、濕度、塑性等不穩(wěn)定，導(dǎo)致土體力學(xué)狀態(tài)極不平衡，其工程性質(zhì)特差，一般不能直接作為建（構(gòu)）筑物的天然地基持力層。

4 三維地質(zhì)模型應(yīng)用

可以對(duì)該三維地質(zhì)體進(jìn)行多種三維空間應(yīng)用分析，對(duì)關(guān)注區(qū)域重點(diǎn)剖析、深入了解模型內(nèi)部的地質(zhì)狀況。如剖面分割，柵欄圖制作分析，屬性模型切割分析、地質(zhì)體透視顯示、生成某一深度的地質(zhì)圖、海量剖面動(dòng)態(tài)顯示等。

（1）虛擬鉆孔分析評(píng)價(jià)

在該模型中，可以在建模范圍的任意處生成虛擬鉆孔及相應(yīng)的鉆孔分層信息（圖12）。本次研究中選取11個(gè)其他地勘項(xiàng)目中深60 m鉆孔柱狀圖與相同位置虛擬鉆孔作對(duì)比，結(jié)果顯示，該模型在地層信息準(zhǔn)確度上極高，除個(gè)別鉆孔薄層不一致外所有地層均相互對(duì)應(yīng)，判斷本模型的虛擬鉆孔可作為地質(zhì)研究參考。

（2）砂土液化評(píng)價(jià)

砂土地震液化是由多種原因綜合作用的結(jié)果，如土的巖性、顆粒組成、密度、埋藏條件、地下水位、沉積環(huán)境、地質(zhì)歷史、地震動(dòng)強(qiáng)度、頻譜特征和持續(xù)時(shí)間等。目前各種判別液化的方法都是經(jīng)驗(yàn)方法，都有一定的局限性和模糊性。在項(xiàng)目中參照GB 50011-2010《建筑抗震設(shè)計(jì)規(guī)范》（2016年版）、GB 50021-2001《巖土工程勘察規(guī)范》（2009年版）相關(guān)規(guī)定，首先進(jìn)行初判，然后對(duì)可能液化的地層采用標(biāo)準(zhǔn)貫入試驗(yàn)法進(jìn)行液化復(fù)判，并采用靜力觸探試驗(yàn)法和剪切波速法進(jìn)行驗(yàn)判（北京市地質(zhì)工程勘察院，2017）。在此基礎(chǔ)上，以標(biāo)準(zhǔn)貫入試驗(yàn)法判別結(jié)果為主，以其他兩種方法判別結(jié)果為輔助參考，將224個(gè)鉆孔的砂土液化信息導(dǎo)入到該三維模型中，獲取研究區(qū)砂土液化空間分布特征（圖13）并進(jìn)一步進(jìn)行分析評(píng)價(jià)，最后從模型中提取砂土液化參數(shù)，得出該區(qū)砂土液化等級(jí)可分為不液化、輕微液化、中等液化和嚴(yán)重液化4級(jí)。其中輕微液化區(qū)主要分布在該區(qū)北部的宋莊翟里—潞城郝家府一帶、東部的西集武辛莊—曹莊一帶及中部十里莊、棗林莊等地;中等液化區(qū)主要分布在北部潮白河西岸的港北村—小楊各莊村、中部張家灣梁各莊附近地區(qū)、縣鎮(zhèn)蘇莊—馬蹄—凌莊—南小莊—前東儀村狹窄U型地帶、南部的尚武集村和胡家村—臨溝屯村地區(qū);嚴(yán)重液化區(qū)主要分布東部潮白河西岸的Ⅰ級(jí)階上的廣大地區(qū)，如白廟—東堡—大沙務(wù)一線、崔家樓—耿樓村一帶;不液化區(qū)主要分布北部潮白河Ⅲ階地以及該區(qū)中西部地區(qū)（北京市地質(zhì)工程勘察院，2019）。因此在該區(qū)進(jìn)行工程建設(shè)需要考慮到砂土液化的影響。

5 結(jié)論

本次建模主要是通過工程地質(zhì)三維解釋好的鉆孔分層數(shù)據(jù)，建立工作區(qū)鉆孔地層控制三角格網(wǎng)，再通過地層控制三角格網(wǎng)，調(diào)整虛擬鉆孔插值個(gè)數(shù)，建立虛擬鉆孔，之后利用原始鉆孔和虛擬鉆孔進(jìn)行三角格網(wǎng)剖分，根據(jù)地層層序，依次構(gòu)建出各個(gè)地層的底板，以保證鉆孔上下的地層接觸關(guān)系一致，最后，圈定地層邊界，利用邊界分塊，根據(jù)抽取的各層底板規(guī)劃到各個(gè)地層塊中，通過融合建立了工作區(qū)的三維地質(zhì)結(jié)構(gòu)模型。

該三維地質(zhì)模型解決了以往二維空間的限制，具體形象地展示了研究區(qū)地下土層的空間三維分布，及評(píng)判結(jié)果的二維和三維分布，既揭露了研究區(qū)的地下土層結(jié)構(gòu)，又展示了易于發(fā)生砂土液化災(zāi)害的空間范圍，為砂土液化災(zāi)害的防治工作提供決策依據(jù)。

模型可服務(wù)于研究區(qū)地基穩(wěn)定性，工程建設(shè)適宜性，工程地質(zhì)分區(qū)和評(píng)價(jià)工作，其三維空間分析能夠快速有效避免因工程建設(shè)而導(dǎo)致的地質(zhì)災(zāi)害的發(fā)生，為城市建設(shè)開展具體的工程提供建設(shè)性的科學(xué)依據(jù)。

參考文獻(xiàn)：

北京市地質(zhì)工程勘察院，2007. 北京通州規(guī)劃新城前期區(qū)域工程地質(zhì)勘察報(bào)告[R].

北京市地質(zhì)工程勘察院，2017. 通州某區(qū)重大地質(zhì)問題調(diào)查與評(píng)價(jià)砂土液化專題成果報(bào)告[R].

北京市地質(zhì)工程勘察院，2019. 通州區(qū)南部（751 km2）重大地質(zhì)問題調(diào)查與評(píng)價(jià)工程地質(zhì)調(diào)查與評(píng)價(jià)專題成果報(bào)告[R].

古銳開，陳偉，歐陽(yáng)春飛，2019. 精細(xì)化三維地質(zhì)建模在城市地下空間開發(fā)中的應(yīng)用[J]. 中國(guó)勘察設(shè)計(jì)（4）：76-79.

何靜，何晗晗，鄭桂森，等，2019. 北京五環(huán)城區(qū)淺部沉積層的三維地質(zhì)結(jié)構(gòu)建模[J]. 中國(guó)地質(zhì)，46（2）：244-254.

李獻(xiàn)峰，薛惠鋒，蘇錦旗，2008. 基于鉆孔信息的二維剖面圖的三維實(shí)體重建[J]. 計(jì)算機(jī)應(yīng)用研究（7）：2107-2109.

李靜，趙帥，2016. 城市三維地質(zhì)建模在砂土液化分析中的應(yīng)用：以通州為例[J]. 中國(guó)礦業(yè)，25（5）：164-168+174.

李敏，劉釗，韓征，等，2018. 城市區(qū)域三維地質(zhì)結(jié)構(gòu)模型建設(shè)與集成方法[J]. 城市地質(zhì)，13（2）：93-97.

馬震，夏雨波，王小丹，等，2019. 雄安新區(qū)工程地質(zhì)勘查數(shù)據(jù)集成與三維地質(zhì)結(jié)構(gòu)模型構(gòu)建[J]. 中國(guó)地質(zhì)，46（S2）：123-138.

楊向敏，吳福，盧志文，等，2019. 淺層地下空間三維地質(zhì)結(jié)構(gòu)建模關(guān)鍵技術(shù)探討：以柳州市為例[J]. 南方國(guó)土資源（7）：81-84+87.

HOULDING S W， 1994. 3D Geoscience modeling：computer techniques for geological characterization[M]. Berlin： Springer-verlag.