劉漢龍 章潤紅 劉東升 仉文崗

摘要：為了推動智慧城市的規(guī)劃和管理，為減輕重慶市地質(zhì)災(zāi)害風(fēng)險等提供科學(xué)依據(jù)和技術(shù)支持，基于大數(shù)據(jù)分析管理和數(shù)據(jù)處理思想，采用大數(shù)據(jù)交匯融合技術(shù)，構(gòu)建包括地區(qū)地層構(gòu)造及巖土材料參數(shù)值域等基礎(chǔ)特征信息的區(qū)域三維地質(zhì)模型?；跀?shù)據(jù)統(tǒng)計、圖像識別方法以及空間統(tǒng)計分析理論和插值方法，綜合利用MATLAB、ArcGIS以及BIM軟件，實(shí)現(xiàn)對重慶市某區(qū)域地質(zhì)鉆孔數(shù)據(jù)的虛擬地層插入、地層曲面擬合、三維實(shí)體模型生成、局部地形深化處理；根據(jù)Bayes參數(shù)估計理論，采用基于歷史大數(shù)據(jù)的先驗(yàn)估計值對工程實(shí)時采集的巖土測試試驗(yàn)參數(shù)進(jìn)行修正，得到各類不同地質(zhì)物理力學(xué)參數(shù)的期望值。

關(guān)鍵詞：圖像識別；地質(zhì)鉆孔數(shù)據(jù)；Bayes參數(shù)估計理論；大數(shù)據(jù)分析；三維地質(zhì)建模

中圖分類號：P642.4 文獻(xiàn)標(biāo)志碼：A 文章編號：2096-6717（2022）06-0001-11

Study on characteristics of physical and mechanical parameters of engineering geology based on data fusion technique LIU Hanlong 1，ZHANG Runhong 2，LIU Dongsheng 3，ZHANG Wengang 1，2

（1.School of Civil Engineering，Chongqing University，Chongqing 400045，P.R.China；2.Institute for Smart City of Chongqing University in Liyang，Liyang 213300，Jiangsu，P.R.China；3.Chongqing Bureau of Geological Survey and Mineral Exploration，Chongqing 401120，P.R.China）

Abstract：In order to promote the planming and management of smart city，and provide scientific basis and technical surpport for reducing the risk of geological disasters in Changqing，based on big data analysis and management as well as data processing，big data convergence and fusion technique are used to construct athree-dimensional geological model of the region，including basic feature information such as regional strata structure and range of geotechrical material parameters.Based on data statistics，image recognition methods，spatial statistical analysis theories and interpolation methods，the MATLAB，ArcGIS and BIM software are used to realize virtual stratum insertion on borehole data，stratum surface fitting，three-dimensional solid model generation，and local terrain deepening processing in aregion of Chongqing.By using Bayes parameter estimation theory and priori estimated values based on the historical big data，the expected values of various geophysical and mechanical parameters are obtained by revising the geotechnical test parameters collected in real-time on site.

Keywords：image recognition;geological borehole data;Bayes parameter estimation theory;big data analysis;three-dimensional geological modeling

城市的快速和高質(zhì)量發(fā)展離不開地質(zhì)環(huán)境條件的支撐。全面掌握城市所處地質(zhì)條件、數(shù)據(jù)化地對城市地質(zhì)特征進(jìn)行描述對智慧城市的建設(shè)尤為重要。近年來，各種三維地質(zhì)建模軟件的迅速發(fā)展、功能的穩(wěn)步提高、利于數(shù)據(jù)信息大力發(fā)展的政策使得三維地質(zhì)建模工作更有意義，更能服務(wù)于國民經(jīng)濟(jì)[1]。目前，三維地質(zhì)建模能將地質(zhì)環(huán)境數(shù)據(jù)、地質(zhì)學(xué)理論與計算機(jī)三維可視化技術(shù)有機(jī)結(jié)合，從而準(zhǔn)確地反映地下空間內(nèi)地質(zhì)構(gòu)造展布、相互接觸關(guān)系以及各地質(zhì)體物理力學(xué)或化學(xué)屬性等特征[2]。

20世紀(jì)90年代以來，隨著計算機(jī)軟硬件技術(shù)的不斷發(fā)展和地球空間信息科學(xué)的蓬勃發(fā)展，能夠反映地質(zhì)構(gòu)造形態(tài)、構(gòu)造關(guān)系及土體參數(shù)屬性變化規(guī)律的三維地質(zhì)建模技術(shù)得到了顯著發(fā)展，取得了眾多有意義的成果[3-8]。建立多維地質(zhì)信息系統(tǒng)，有利于一體化存儲、管理和處理地下—地上、地質(zhì)—地理、空間—屬性數(shù)據(jù)，進(jìn)行三維可視化地質(zhì)建模，有利于對地質(zhì)對象的總體把握和整體分析?；诤Ａ慷嘣础悩?gòu)、異質(zhì)的地質(zhì)勘查數(shù)據(jù)，構(gòu)建其三維地質(zhì)模型，可直觀而形象地展現(xiàn)地質(zhì)體和地質(zhì)結(jié)構(gòu)，從而顯著提高對地質(zhì)現(xiàn)象、地質(zhì)資源和地質(zhì)環(huán)境的認(rèn)知能力。

目前，開展三維地質(zhì)信息系統(tǒng)建設(shè)和三維數(shù)字圖幅、數(shù)字礦山、數(shù)字油田、數(shù)字工程和數(shù)字地災(zāi)點(diǎn)地質(zhì)建模，已經(jīng)成為國際地質(zhì)部門和地質(zhì)科學(xué)界技術(shù)創(chuàng)新的重要領(lǐng)域[5-8]。

對于城市地質(zhì)而言，三維地質(zhì)建模對地質(zhì)災(zāi)害防御和地下設(shè)施綜合管理具有巨大的作用，近20年來，三維地質(zhì)建模相關(guān)研究取得了許多可喜的成果[9-18]，三維地質(zhì)建模與地質(zhì)大數(shù)據(jù)的融合研究越來越受到重視[19-20]。薛林福等[21]提出了分塊三維地質(zhì)建模方法，以斷裂、巖體邊界、不整合接觸面等為邊界，將較為復(fù)雜的三維地質(zhì)建模工作區(qū)分解為多個內(nèi)部構(gòu)造相對簡單的、建模難度較小的地質(zhì)單元，為開展復(fù)雜山區(qū)三維地質(zhì)建模提供了精確依據(jù)。

董梅等[22]利用GOCAD建立了北京市朝陽區(qū)某建筑場地的土層地質(zhì)模型，并提出了以Kriging插值和離散光滑插值相互結(jié)合來建立地層界面以及采用GOCAD中的GRID功能代替SOLID功能建立三維地層實(shí)體的方法。王駿等[23]以重慶歌樂山典型背斜構(gòu)造區(qū)為建模研究對象，采用搜集資料→野外踏勘→室內(nèi)建模的工作路線，運(yùn)用SKUA-GOCAD的workflow建模流程，建立了構(gòu)造模型，并運(yùn)用數(shù)據(jù)趨勢分析方法建立了巖性模型。李青元等[24]從多角度全面總結(jié)了三維地質(zhì)建模技術(shù)的用途、現(xiàn)狀、存在問題、發(fā)展趨勢，并結(jié)合三維建模的現(xiàn)狀提出了“相關(guān)部門應(yīng)對三維地質(zhì)建模示范工程制定總體風(fēng)格與基本內(nèi)容一致、適合本行業(yè)專業(yè)特點(diǎn)的專用標(biāo)準(zhǔn)”。易宗旺等[25]利用區(qū)域調(diào)查中的地形數(shù)據(jù)、遙感影像數(shù)據(jù)、鉆孔數(shù)據(jù)、剖面數(shù)據(jù)等各種資料，在三維地質(zhì)建模系統(tǒng)（GeoBIM）平臺中進(jìn)行三維地質(zhì)建模等。劉東升等[26]研發(fā)了巖土材料特征大數(shù)據(jù)管理及分析軟件系統(tǒng)，實(shí)現(xiàn)了對巖土材料大數(shù)據(jù)實(shí)施收集、傳輸、分類、篩選、管理和統(tǒng)計分析的功能，并結(jié)合重慶地區(qū)收集到的7萬余條典型巖土材料數(shù)據(jù)對該地區(qū)的典型巖土材料特征進(jìn)行分析。由于既往積累的地質(zhì)環(huán)境數(shù)據(jù)相對較分散，盡管總量巨大，但缺乏整合和必要的分析，使得其成為“死數(shù)據(jù)”，未能建立起城市地質(zhì)基礎(chǔ)特征系統(tǒng)，無法為智慧城市規(guī)劃和建設(shè)提供幫助。

筆者基于重慶市地勘單位已有的大量巖土勘察與地形地貌數(shù)據(jù)資源，基于大數(shù)據(jù)分析管理和數(shù)據(jù)處理思想，采用大數(shù)據(jù)交匯融合技術(shù)，構(gòu)建重慶市區(qū)域三維地質(zhì)模型，生成了包括地區(qū)地層構(gòu)造及巖土材料參數(shù)值域等基礎(chǔ)特征信息的重慶區(qū)域地質(zhì)條件圖譜?；跀?shù)據(jù)統(tǒng)計、圖像識別方法以及空間統(tǒng)計分析理論和插值方法，綜合利用MATLAB、ArcGIS以及BIM軟件，對鉆孔數(shù)據(jù)進(jìn)行虛擬地層插入、地層曲面擬合、三維實(shí)體模型生成、局部地形深化處理；運(yùn)用Bayes參數(shù)估計理論，采用基于歷史2第6期劉漢龍，等：基于數(shù)據(jù)融合的工程地質(zhì)物理力學(xué)參數(shù)特征研究大數(shù)據(jù)的先驗(yàn)估計值對工程實(shí)時采集的巖土測試試驗(yàn)參數(shù)進(jìn)行修正，得到各類不同地質(zhì)物理力學(xué)參數(shù)的期望值。

1鉆孔地質(zhì)數(shù)據(jù)

收集地質(zhì)鉆孔數(shù)量共計3456個，鉆孔分布如圖1所示，典型地質(zhì)鉆孔具體信息如圖2所示。由收集到的地質(zhì)鉆孔分布可知，地質(zhì)勘探點(diǎn)主要分布在朝天門區(qū)域。

1.1空間幾何參數(shù)

根據(jù)鉆孔資料和區(qū)域資料顯示，研究區(qū)域的巖層產(chǎn)狀110°∠6°，地層巖性為砂、泥巖互層，巖層層面裂隙張開，無填充，結(jié)合差，為軟弱結(jié)構(gòu)面。區(qū)域范圍內(nèi)主要有兩組裂隙，其產(chǎn)狀分別為：1）L1產(chǎn)狀170°∠78°，裂隙間距0.8～1.5 m，裂面平，略呈波狀，多黑灰色，無充填，開口1～3 mm，局部充填鈣質(zhì)，延伸長約1～3 m，貫通性較差。裂隙面均為軟弱結(jié)構(gòu)面，結(jié)構(gòu)面結(jié)合差；2）裂隙產(chǎn)狀235°∠81°，裂隙間距1.0～2.5 m，裂面粗糙，呈波狀，多呈閉合狀，局部開口2～5 mm，無充填，貫通性較差。裂隙面均為軟弱結(jié)構(gòu)面，結(jié)構(gòu)面結(jié)合差。據(jù)地面調(diào)查及鉆探揭露，部分區(qū)域內(nèi)出露地層巖性有第四系人工填土層、侏羅系中統(tǒng)沙溪廟組砂巖及泥巖。

1.2物理力學(xué)參數(shù)

為揭示擬建模區(qū)域的地質(zhì)構(gòu)造、元素組成及力學(xué)性能，采用傳統(tǒng)鉆孔巖樣試驗(yàn)和基于圖像識別技術(shù)的方法獲取地質(zhì)物理力學(xué)參數(shù)?；贑T、SEM掃描和三維重構(gòu)技術(shù)的巖樣孔隙率，確定運(yùn)用Avizo8.0軟件對巖樣CT掃描切片進(jìn)行中值過濾、二值化處理和精細(xì)化加工，之后再對每組切片進(jìn)行數(shù)字巖芯重構(gòu)、孔吼模型建立、孔裂隙提取及孔隙率計算，如圖3所示。

圖3真實(shí)砂巖試件及其3種模型Fig.3 Real sandstone specimen and its three models基于EDS能譜分析技術(shù)對巖樣元素組成規(guī)律進(jìn)行確定。采用Quanta 250掃描電子顯微鏡和配套的X射線能譜儀進(jìn)行巖樣微結(jié)構(gòu)中的元素測定。

在每個巖樣切片的掃描電鏡圖像中，至少選取3個微結(jié)構(gòu)EDS能譜分析區(qū)，如圖4所示。對每個微結(jié)構(gòu)區(qū)均進(jìn)行元素測定分析，將三者的平均值作為最終結(jié)果。

基于巖石細(xì)觀圖像深度學(xué)習(xí)的巖性智能識別與分類。采用SN0745-60U2K型2K測量電子顯微鏡進(jìn)行巖樣細(xì)觀圖像拍攝，電子放大倍數(shù)為15～180倍，拍攝圖像像素高達(dá)1600萬，并具備自動尋邊、HDR高動態(tài)、自動存儲等功能，操作簡便，對巖樣物理力學(xué)測試結(jié)果影響小。

巖樣端面細(xì)觀圖像經(jīng)細(xì)致篩分后，將其中的70%作為訓(xùn)練集、20%作為驗(yàn)證集、10%作為測試集，分別運(yùn)用Inception-v3、Xception、NASNetLarge和InceptionResNetV2這4種深度卷積神經(jīng)網(wǎng)絡(luò)模型及GoogleNet圖像分類模型，同時，結(jié)合Softmax回歸模型等共同構(gòu)建巖（土）性材料的識別模型，并對其進(jìn)行循環(huán)訓(xùn)練學(xué)習(xí)，深度學(xué)習(xí)流程如圖5所示。

圖6為Inception-v3深度神經(jīng)網(wǎng)絡(luò)模型的結(jié)構(gòu)展示，該模型首先是3個卷積層，然后連接一個最大池化層，再設(shè)兩個卷積層，再連接一個最大池化層，最后連接11個不同的Inception結(jié)構(gòu)混合層。其余深度卷積神經(jīng)網(wǎng)絡(luò)模型不再列出。

采用該方法對所采集的4種典型巖樣圖像進(jìn)行旋轉(zhuǎn)、偏移處理，生成深度學(xué)習(xí)樣本，隨后，基于Inception-v3、Xception、NASNetLarge和Inception ResNetV2這4種深度卷積神經(jīng)網(wǎng)絡(luò)模型，利用遷移學(xué)習(xí)的方式形成巖石細(xì)觀圖像集分類的深度學(xué)習(xí)模型，最后使用集成學(xué)習(xí)算法，通過加權(quán)評估得出最終的判別結(jié)果，從而實(shí)現(xiàn)巖（土）性材料圖像特征的提取與識別，并對其進(jìn)行自動、高效、可靠地分類。

在訓(xùn)練過程中，對訓(xùn)練集和驗(yàn)證集中的巖石圖像分別進(jìn)行檢驗(yàn)評估。訓(xùn)練結(jié)果顯示：在循環(huán)訓(xùn)練1000次后，Inception-v3深度卷積神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練集中的巖石細(xì)觀圖像分類準(zhǔn)確率達(dá)到92.77%，驗(yàn)證集中的巖石細(xì)觀圖像分類準(zhǔn)確率為76.31%，如圖7（a）所示；在循環(huán)訓(xùn)練300次后，Xception深度卷積神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練集中的巖石細(xì)觀圖像分類準(zhǔn)確率為93.30%，驗(yàn)證集中的巖石細(xì)觀圖像分類準(zhǔn)確率為78.27%，如圖7（b）所示；在循環(huán)訓(xùn)練300次后，NASNetLarge深度卷積神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練集中的巖石細(xì)觀圖像分類準(zhǔn)確率為90.52%，驗(yàn)證集中的巖石細(xì)觀圖像分類準(zhǔn)確率為74.74%，如圖7（c）所示；在循環(huán)訓(xùn)練300次后，Inception?ResNetV2深度卷積神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練集中的巖石細(xì)觀圖像分類準(zhǔn)確率為90.81%，驗(yàn)證集中的巖石細(xì)觀圖像分類準(zhǔn)確率為78.81%，如圖7（d）所示。

經(jīng)加權(quán)評估后，最終訓(xùn)練集中的巖石細(xì)觀圖像分類準(zhǔn)確率將達(dá)到95%以上，驗(yàn)證集中的巖石細(xì)觀圖像分類準(zhǔn)確率高于80%。

相比于傳統(tǒng)的機(jī)器學(xué)習(xí)方法，該深度學(xué)習(xí)方法通過圖像像素點(diǎn)組成深度學(xué)習(xí)的輸入層，利用多個深度可分離卷積網(wǎng)絡(luò)構(gòu)建復(fù)雜的BP神經(jīng)網(wǎng)絡(luò)，最終得到的模型在巖石巖性的識別與分類上具有更高的魯棒性和泛化能力。此外，集成學(xué)習(xí)算法的運(yùn)用顯著提高了巖石細(xì)觀圖像分類的準(zhǔn)確率，使巖性判別結(jié)果更可信。

2地質(zhì)數(shù)據(jù)分析與建模

2.1三維地質(zhì)空間幾何數(shù)據(jù)分析與建模

為了解決數(shù)據(jù)與成本的矛盾，基于空間統(tǒng)計分析理論選擇合適的插值方法，利用有限的采樣點(diǎn)數(shù)據(jù)擬合其他相對準(zhǔn)確、精度較高的高程數(shù)據(jù)。常用的空間插值方法有三角剖分法、距離倒數(shù)加權(quán)法、最小曲率法、最近鄰點(diǎn)法、Kriging方法，其中，Kriging方法通過引進(jìn)以距離為自變量的半變差函數(shù)來計算權(quán)值，由于半變差函數(shù)既可以反應(yīng)變量的空間結(jié)構(gòu)特性，又可以反映變量的隨機(jī)分布特性，利用Kriging方法進(jìn)行空間數(shù)據(jù)插值往往可以取得理想的效果。另外，通過設(shè)計變差函數(shù)，Kriging方圖74種深度神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練結(jié)果Fig.7 Four types of deep neural network model training results5第44卷土木與環(huán)境工程學(xué)報（中英文）法很容易實(shí)現(xiàn)局部加權(quán)插值，克服了一般距離加權(quán)插值結(jié)果的不穩(wěn)定性。因此，采用Kriging插值對地質(zhì)巖性進(jìn)行預(yù)測。

利用單一軟件進(jìn)行三維地質(zhì)建模存在一些弊端，導(dǎo)致一些缺陷，無法完全實(shí)現(xiàn)預(yù)期的功能。因此，綜合利用MATLAB、ArcGIS以及BIM軟件進(jìn)行精細(xì)化三維地質(zhì)建模，具體建模流程見圖8。

第1步：將地勘報告內(nèi)的鉆孔及其他有用數(shù)據(jù)進(jìn)行提取和篩選，形成地質(zhì)信息基礎(chǔ)數(shù)據(jù)庫，儲存在Excel表格內(nèi)；

第2步：基于虛擬地層插入原理，利用MATLAB對鉆孔數(shù)據(jù)進(jìn)行預(yù)處理，形成標(biāo)準(zhǔn)鉆孔地層分布；

第3步：將標(biāo)準(zhǔn)化后的鉆孔數(shù)據(jù)導(dǎo)入ArcGIS中，基于ArcGIS自帶工具以及Python編程，建立ModelBuilder建模流程，得到TIN地層模型；

第4步：將TIN地層模型轉(zhuǎn)換為Auto Desk平臺能夠識別的DWG文件；

第5步：將地層模型文件導(dǎo)入Civil 3D內(nèi)，然后結(jié)合其他地層信息和地層巖性經(jīng)驗(yàn)分析進(jìn)行手動概化處理，建立三維地質(zhì)實(shí)體模型；

第6步：結(jié)合地表的建筑信息、基坑信息以及開洞信息，利用Revit軟件對Civil 3D模型進(jìn)行深化處理，建立較為精細(xì)的三維地質(zhì)模型；第7步：將所建立的模型導(dǎo)入Navisworks內(nèi)進(jìn)行三維漫游展示；第8步：在Auto Desk平臺進(jìn)行二次開發(fā)，編寫可以對模型添加屬性的程序，然后對前面建立的三維地質(zhì)模型指定巖層屬性，形成動態(tài)的特征信息圖譜。

2.2三維地質(zhì)物理力學(xué)參數(shù)分析與估計

2.2.1巖土參數(shù)先驗(yàn)分布及超參數(shù)確定方法 先驗(yàn)分布的確定受控于先驗(yàn)信息的提取模式，而先驗(yàn)信息的提取模式必須與歷史數(shù)據(jù)的特點(diǎn)相符合。

提出一種新的先驗(yàn)信息提取模式——概率模式，根據(jù)歷史數(shù)據(jù)，有R個工程的數(shù)據(jù)，其中每個工程有n i個巖土強(qiáng)度參數(shù)x的樣本，按實(shí)際工程分組分別統(tǒng)計均值xˉn1與方差s2 n1，形成一組新的二元隨機(jī)樣本{（xˉn1，s2n1），（xˉn2，s2n2）…（xˉ2 nR，s2 nR）}，認(rèn)為這組樣本來自先驗(yàn)分布π（μ，σ2），通過樣本參數(shù)估計對先驗(yàn)分布進(jìn)行擬合，從而確定先驗(yàn)分布π（μ，σ2），如圖9所示。

與單一信息源模式及多信息源權(quán)重融合模式相比，提出的概率模式有3個特點(diǎn)：1）歷史數(shù)據(jù)來自不同的工程，因此，認(rèn)為R個工程樣本來自R個正態(tài)或?qū)?shù)正態(tài)總體，這符合巖土參數(shù)隨機(jī)分布、具有明顯區(qū)域性的屬性。而單一信息源模式認(rèn)為所有數(shù)據(jù)服從同一個分布，多信息源差異權(quán)重融合模式雖然提出了存在多個總體，但對歷史數(shù)據(jù)不是按區(qū)域或工程進(jìn)行分組，分組具有一定主觀性，分組的數(shù)量也較少，只有3組。

2）按照貝葉斯理論，先驗(yàn)分布是針對某一個參數(shù)而言，這里的參數(shù)是指巖土強(qiáng)度指標(biāo)x所服從的正態(tài)總體的兩個未知參數(shù)（μ，σ2），而原始數(shù)據(jù)是關(guān)于x的樣本。單信息源模式是在設(shè)定σ2已知的情況下通過歷史數(shù)據(jù)得到x～N（a，b），然后間接地推測μ～N（a，b/m 0）即先驗(yàn)分布π（μ）；多信息源差異權(quán)重融合模式是在單信息源模式的基礎(chǔ)上認(rèn)為存在多個先驗(yàn)分布π i（θ），θ表示μ或σ2，根據(jù)差異權(quán)重獲得綜合先驗(yàn)分布。而概率模式是在原始數(shù)據(jù)的基礎(chǔ)上結(jié)合樣本來自不同工程的特點(diǎn)，進(jìn)行了二次處理，分別統(tǒng)計了各個工程的均值和方差，形成了一組新的歷史數(shù)據(jù)樣本{（xˉn1，s2n1），（xˉn2，s2n2）…（xˉnR，s 2 nR）}，是直接針對參數(shù)（μ，σ2）的先驗(yàn)樣本，這樣就減少了間接推測過程中主觀性的影響。另外，提出的概率模式將（μ，σ2）按二元隨機(jī)變量考慮，而其他兩種模式是在假設(shè)其中一個參數(shù)已知的情況下的一元隨機(jī)變量的先驗(yàn)分布。

3）單一信息源模式實(shí)質(zhì)是一種工程類比的模式，即認(rèn)為歷史數(shù)據(jù)源于統(tǒng)一總體x～N（a，b），認(rèn)為這個總體與待估總體x～N（μ，σ2）存在工程類比性，因此可以推測μ～N（a，b/m 0）。多信息源權(quán)重融合模式與單一信息源模式的核心相同，只是認(rèn)為存在多個相似工程，按差異權(quán)重形成一個綜合的相似工程，然后按單一信息源模式推測先驗(yàn)分布。而概率模式則是將（μ，σ2）視為二元隨機(jī)變量，根據(jù)其歷史樣本{（xˉn1，s2n1），（xˉn2，s2n2）…（xˉ2nR，s2 nR）}得到其先驗(yàn)分布。

2.2.2典型區(qū)域物理力學(xué)參數(shù)估計 巖土強(qiáng)度參數(shù)的貝葉斯估計是綜合先驗(yàn)信息、總體信息和現(xiàn)場信息得到的，其中，先驗(yàn)分布是基于某一隨機(jī)特征參數(shù)而言。提出一種針對巖土強(qiáng)度參數(shù)總體分布特征（μ，σ2）的聯(lián)合先驗(yàn)分布，該分布基于大量歷史數(shù)據(jù)，并根據(jù)巖土參數(shù)的隨機(jī)分布特征，經(jīng)過一定的處理提煉而得。根據(jù)提出的基于歷史大數(shù)據(jù)的巖土強(qiáng)度參數(shù)貝葉斯估計方法，利用重慶市某工程所處區(qū)域的歷史數(shù)據(jù)，結(jié)合現(xiàn)場取樣測試的結(jié)果，確定巖土強(qiáng)度參數(shù)的貝葉斯估計值，如表1、表2所示。表1中，先驗(yàn)信息樣本數(shù)量為38，現(xiàn)場樣本數(shù)量為8；表2中，先驗(yàn)信息樣本數(shù)量為47，現(xiàn)場樣本數(shù)量為8。

2.3空間幾何模型與物理力學(xué)參數(shù)模型的融合 地質(zhì)數(shù)據(jù)多源性是三維地質(zhì)建模最大的特點(diǎn)，模型構(gòu)建的關(guān)鍵是將這些數(shù)據(jù)有效地融合，以提高模型的應(yīng)用效率，指導(dǎo)工程設(shè)計。因此，在分別建立了三維地質(zhì)的空間幾何模型和物理力學(xué)參數(shù)模型后進(jìn)行數(shù)據(jù)的交匯融合，形成既包含幾何信息又包含物理信息的三維地質(zhì)模型（見圖10）。在數(shù)據(jù)交匯融合的過程中，涉及多方面的技術(shù)問題，如統(tǒng)一所有數(shù)據(jù)的坐標(biāo)系和比例尺、地表和地下數(shù)據(jù)的有效融合及主要數(shù)據(jù)和次要數(shù)據(jù)的有效融合等，主要包括：1）統(tǒng)一所有具有多源、多尺度、多分辨率數(shù)據(jù)的坐標(biāo)系和比例尺，構(gòu)建可以作為后期模型構(gòu)建的數(shù)據(jù)源和約束數(shù)據(jù)的數(shù)據(jù)庫；2）基于地表填圖和遙感解譯的地質(zhì)數(shù)據(jù)及鉆孔、中段平面圖、勘探線剖面圖、物探解譯的深部地質(zhì)數(shù)據(jù)進(jìn)行有效融合和插值處理；3）覆蓋整個建模區(qū)并足以構(gòu)建模型的主要數(shù)據(jù)和次要數(shù)據(jù)的有效融合。

3典型區(qū)域三維地質(zhì)模型

3.1空間幾何模型建模及深化處理

基于空間統(tǒng)計分析理論，在ArcGIS軟件中通過插值分析得到地層模型。為便于在Civil 3D中建模，將建立的三維雜填土、砂巖以及泥巖地層模型轉(zhuǎn)換為等高線模型輸出，如圖11所示。

在Civil 3D中使用曲面功能導(dǎo)入雜填土層的等高線，生成曲面，并通過提取實(shí)體生成三維實(shí)體；按此方法依次生成砂巖及泥巖三維實(shí)體，如圖12所示。由于所建三維地質(zhì)實(shí)體模型只是地質(zhì)模型，未反映出地下建筑等信息，因此，需要對模型進(jìn)行深化處理，以反映真實(shí)地層與建筑物的關(guān)系。

擬建模區(qū)域存在大量的地下建筑，因此，需在三維實(shí)體模型的相應(yīng)位置下方開挖基坑，其基坑建模步驟與隧道建模類似，需使用Civil 3D中實(shí)體編輯模塊的實(shí)體差集功能使三維地質(zhì)實(shí)體模型減去相應(yīng)的實(shí)體生成基坑。為更真實(shí)地反映擬建模區(qū)域的情況，在開挖隧道和基坑的三維地質(zhì)模型中建立相應(yīng)的建筑模型，其方式為在Revit軟件中參考實(shí)際建筑圖紙建立建筑模型，具體步驟為：將Revit中內(nèi)建模型通過實(shí)體拉伸和空模型拉伸建立建筑模型，將建立的建筑模型從Revit中導(dǎo)入Civil 3D中，并裝配在相應(yīng)位置，建立的模型如圖13所示。

3.2融合的三維地質(zhì)BIM模型地層物性參數(shù)

通過Civil 3D建立區(qū)域三維地質(zhì)模型，利用Civil 3D功能和Revit建模功能建立隧道、基坑和建筑三維地質(zhì)模型并進(jìn)行深化處理。為給后續(xù)其他工程提供設(shè)計指導(dǎo)和建議，還需對三維地質(zhì)模型的各個地層賦予相應(yīng)的地質(zhì)物性參數(shù)，為此，在Civil3D建模平臺的基礎(chǔ)上開發(fā)了可以賦予地質(zhì)物性參數(shù)的插件（RWCD.bundle）。其使用方法為：

1）安裝插件（RWCD.bundle）；

2）打開Civil 3D建模平臺；

3）讀取或設(shè)置數(shù)據(jù)。

單個對象的數(shù)據(jù)讀寫如圖14所示，選取目標(biāo)三維地質(zhì)實(shí)體，選擇寫入數(shù)據(jù)或讀取數(shù)據(jù)。

地層物性參數(shù)賦值流程為：

1）人工雜填土層物性參數(shù)賦值。根據(jù)地質(zhì)物性研究報告獲得人工雜填土的基底摩擦系數(shù)為0.25，在Civil 3D中輸入三維雜填土地層參數(shù)。

2）砂巖層物性參數(shù)賦值。根據(jù)地質(zhì)物性研究報告獲得砂巖的天然重度為24.53 kN/m 3，天然抗壓強(qiáng)度標(biāo)準(zhǔn)值為39.6 MPa，飽和抗壓強(qiáng)度標(biāo)準(zhǔn)值為32 MPa，地基承載力特征值為14.36 MPa，地基極限承載力標(biāo)準(zhǔn)值為43.54 MPa，黏聚力為2.12 MPa，內(nèi)摩擦角為36.48°，變形模量為6600 MPa，彈性模量為7100 MPa，泊松比為0.16，抗拉強(qiáng)度為0.84 MPa，基底摩擦系數(shù)為0.65，在Civil 3D中輸入三維砂巖層參數(shù)。

3）泥巖層物性參數(shù)賦值。根據(jù)地質(zhì)物性研究報告獲得泥巖的天然重度為24.7 kN/m 3，天然抗壓強(qiáng)度標(biāo)準(zhǔn)值為9.8 MPa，飽和抗壓強(qiáng)度標(biāo)準(zhǔn)值為6.5 MPa，地基承載力特征值為6.8 MPa，地基極限承載力標(biāo)準(zhǔn)值為10.78 MPa，黏聚力為0.67 MPa，內(nèi)摩擦角為33.55°，變形模量為1460 MPa，彈性模量為1700 MPa，泊松比為0.33，抗拉強(qiáng)度為0.27 MPa，基底摩擦系數(shù)為0.50，在Civil 3D中輸入三維泥巖層參數(shù)。

對各地層賦予地質(zhì)物性參數(shù)后，若需要讀取目標(biāo)地層的地質(zhì)物性參數(shù)，則只需選取目標(biāo)地層并選取讀取參數(shù)，即可獲得目標(biāo)地層參數(shù)。

4結(jié)論

1）利用大量的巖土材料勘察試驗(yàn)數(shù)據(jù)資源，結(jié)合先進(jìn)的巖土參數(shù)試驗(yàn)方法、圖像識別及機(jī)器學(xué)習(xí)方法對巖土材料物理力學(xué)參數(shù)進(jìn)行分析研究。在此基礎(chǔ)上，利用大數(shù)據(jù)分析管理和數(shù)據(jù)處理思想，采用大數(shù)據(jù)交匯融合技術(shù)，將試驗(yàn)和收集所得的地質(zhì)鉆孔數(shù)據(jù)（幾何數(shù)據(jù)和物理數(shù)據(jù)）進(jìn)行分類篩選。

基于已有鉆孔數(shù)據(jù)，利用數(shù)學(xué)插值方法對空白區(qū)虛擬鉆孔的幾何特征數(shù)據(jù)進(jìn)行分析和推送，有效生成虛擬數(shù)字鉆孔。

2）基于已有的巖土材料試驗(yàn)測試試驗(yàn)參數(shù)，利用貝葉斯理論對空白區(qū)同類巖性材料的物理力學(xué)參數(shù)進(jìn)行估計，獲得相應(yīng)的巖土材料物理力學(xué)參數(shù)估計值，彌補(bǔ)了數(shù)據(jù)融合中巖土材料參數(shù)不足的問題。

3）利用實(shí)際獲取和計算生成的鉆孔地質(zhì)特征數(shù)據(jù)，從幾何特征和物理特征兩方面進(jìn)行大數(shù)據(jù)融合并進(jìn)行數(shù)據(jù)建模，得到一定區(qū)域范圍內(nèi)地質(zhì)巖層特征的三維空間分布及物理力學(xué)特征完整表達(dá)模型，實(shí)現(xiàn)區(qū)域地層結(jié)構(gòu)特征描述的數(shù)字化和可視化。

參考文獻(xiàn)

[1]童亮.GOCAD在某橋基三維地質(zhì)建模中的應(yīng)用研究[D].成都：西南交通大學(xué)，2009.TONG L.The study on application of three-dimensional modeling in abridges foundation by using the software named GOCAD[D].Chengdu：Southwest Jiaotong University，2009.（in Chinese）

[2]張洋洋，周萬蓬，吳志春，等.三維地質(zhì)建模技術(shù)發(fā)展現(xiàn)狀及建模實(shí)例[J].東華理工大學(xué)學(xué)報（社會科學(xué)版），2013，32（3）：403-409.ZHANG YY，ZHOU WP，WU ZC，et al.The developmentstatusof3Dgeologicalmodeling technology and modeling instances[J].Journal of East China Institute of Technology（Social Science），2013，32??? （3）：403-409.（in Chinese）

[3]AKPOKODJEEG.Theengineering-geological characteristics and classification of the major superficial soils of the Niger Delta[J].Engineering Geology，1987，23（3/4）：193-211.

[4]DASSARGUESA，BIVERP，MONJOIEA.Geotechnical properties of the Quaternary sediments in Shanghai[J].Engineering Geology，1991，31（1）：71-90.

[5]ZHU LF，ZHANG CJ，LI MJ，et al.Building 3D solid models of sedimentary stratigraphic systems from borehole data：An automatic method and case studies[J].Engineering Geology，2012，127：1-13.

[6]TONINI A，GUASTALDI E，MASSA G，et al.3D geo-mapping based on surface data for preliminary study of underground works：A case study in Val Topina （Central Italy）[J].Engineering Geology，2008，99（1/2）：61-69.

[7]TURNER AK.Challenges and trends for geological modelling and visualisation[J].Bulletin of Engineering Geology and the Environment，2006，65（2）：109-127.

[8]WELLMANN JF，HOROWITZ FG，SCHILL E，et al.Towards incorporating uncertainty of structural data in 3D geological inversion[J].Tectonophysics，2010，490（3/4）：141-151.

[9]林良俊，李亞民，葛偉亞，等.中國城市地質(zhì)調(diào)查總體構(gòu)想與關(guān)鍵理論技術(shù)[J].中國地質(zhì)，2017，44（6）：1086-1101.LIN LJ，LI YM，GE WY，et al.General ideas for urban geological survey in China and key theory and techniques[J].Geology in China，2017，44（6）：1086-1101.（in Chinese）

[10]朱良峰，吳信才，劉修國，等.基于鉆孔數(shù)據(jù)的三維地層模型的構(gòu)建[J].地理與地理信息科學(xué)，2004，20（3）：26-30.ZHU LF，WU XC，LIU XG，et al.Reconstruction of3D strata model based on borehole data[J].Geography and Geo-Information Science，2004，20（3）：26-30.（in Chinese）

[11]潘懋，方裕，屈紅剛.三維地質(zhì)建模若干基本問題探討[J].地理與地理信息科學(xué)，2007，23（3）：1-5.PAN M，F(xiàn)ANG Y，QU HG.Discussion on several foundationalissuesinthree-dimensionalgeological modeling[J].Geography and Geo-Information Science，2007，23（3）：1-5.（in Chinese）

[12]盧鵬飛，龍奎，楊春，等.重慶都市經(jīng)濟(jì)圈城市地質(zhì)三維地質(zhì)結(jié)構(gòu)模型構(gòu)建方法[J].地質(zhì)學(xué)刊，2017，41（1）：79-84.LU PF，LONG K，YANG C，et al.Study on model construction methods of three-dimensional geological structure in the urban economic circle of Chongqing[J].Journal of Geology，2017，41（1）：79-84.（in Chinese）

[13]陳國旭，吳沖龍，張夏林，等.三維地質(zhì)建模與地礦勘查圖件編制一體化方法研究[J].地質(zhì)與勘探，2010，46 （3）：542-546.CHEN GX，WU CL，ZHANG XL，et al.Study on integration of 3D geological modeling and mineral resourceexplorationmapping[J].Geologyand Exploration，2010，46（3）：542-546.（in Chinese）

[14]張?jiān)海瑒⒁?，許苗娟，等.北京三維地質(zhì)建模研究與實(shí)踐[J].工程勘察，2015，43（6）：60-65.ZHANG Y，LIU Y，XU MJ，et al.Research and practiceofBeijing3Dgeologicalmodeling[J].Geotechnical Investigation&Surveying，2015，43（6）：60-65.（in Chinese）

[15]李青元，張洛宜，曹代勇，等.三維地質(zhì)建模的用途、現(xiàn)狀、問題、趨勢與建議[J].地質(zhì)與勘探，2016，52（4）：759-767.LI QY，ZHANG LY，CAO DY，et al.Usage，status，problems，trends and suggestions of 3D geological modeling[J].Geology and Exploration，2016，52（4）：759-767.（in Chinese）

[16]明鏡.三維地質(zhì)建模技術(shù)研究[J].地理與地理信息科學(xué)，2011，27（4）：14-18，56.MING J.A study on three-dimensional geological modeling[J].Geography and Geo-Information Science，2011，27（4）：14-18，56.（in Chinese）

[17]何靜，何晗晗，鄭桂森，等.北京五環(huán)城區(qū)淺部沉積層的三維地質(zhì)結(jié)構(gòu)建模[J].中國地質(zhì)，2019，46（2）：244-254.HE J，HE HH，ZHENG GS，et al.3D geological modelling of superficial deposits in Beijing City[J].Geology in China，2019，46（2）：244-254.（in Chinese）

[18]施木俊，熊毅明，甄云鵬.基于工程勘察鉆孔數(shù)據(jù)的三維地層模型的自動構(gòu)建[J].城市勘測，2006（5）：62-65.SHI MJ，XIONG YM，ZHEN YP.Auto-reconstruction of 3D strata model based on borehole data ofprojectinvestigation[J].UrbanGeotechnical Investigation&Surveying，2006（5）：62-65.（in Chinese）

[19]朱月琴，譚永杰，張建通，等.基于Hadoop的地質(zhì)大數(shù)據(jù)融合與挖掘技術(shù)框架[J].測繪學(xué)報，2015，44（Sup1）：152-159.ZHU YQ，TAN YJ，ZHANG JT，et al.A framework of hadoop based geology big data fusion and miningtechnologies[J].ActaGeodaeticaet Cartographica Sinica，2015，44（Sup1）：152-159.（in Chinese）

[20]潘婷婷，陳建平，吳永亮，等.多源異構(gòu)地質(zhì)數(shù)據(jù)集成方法應(yīng)用研究[J].地質(zhì)學(xué)刊，2018，42（1）：122-126.PAN TT，CHEN JP，WU YL，et al.The applied research of integration method for multi-source and heterogeneous spatial geological data[J].Journal of Geology，2018，42（1）：122-126.（in Chinese）

[21]薛林福，李文慶，張偉，等.分塊區(qū)域三維地質(zhì)建模方法[J].吉林大學(xué)學(xué)報（地球科學(xué)版），2014，44（6）：2051-2058.XUE LF，LI WQ，ZHANG W，et al.A method of block-divided 3D geologic modeling in regional scale[J].Journal of Jilin University（Earth Science Edition），2014，44（6）：2051-2058.（in Chinese）

[22]董梅，慎乃齊，胡輝，等.基于GOCAD的三維地質(zhì)模型構(gòu)建方法[J].桂林工學(xué)院學(xué)報，2008，28（2）：188-192. DONG M，SHEN NQ，HU H，et al.3D geological modeling method based on GOCAD[J].Journal of Guilin University of Technology，2008，28（2）：188-192.（in Chinese）

[23]王駿，陳立川，黃勤，等.基于數(shù)據(jù)趨勢分析的三維地質(zhì)模型構(gòu)建：以歌樂山地區(qū)為例[J].中國金屬通報，2018（8）：214-215.WANG J，CHEN LC，HUANG Q，et al.Construction of 3D geological model based on data trend analysis：Taking Geleshan area as an example[J].China Metal Bulletin，2018（8）：214-215.（in Chinese）

[24]李青元，張麗云，魏占營，等.三維地質(zhì)建模軟件發(fā)展現(xiàn)狀及問題探討[J].地質(zhì)學(xué)刊，2013，37（4）：554-561.LI QY，ZHANG LY，WEI ZY，et al.On 3D geologicalmodelingsoftwaredevelopmentand discussions on several issues[J].Journal of Geology，2013，37（4）：554-561.（in Chinese）

[25]易宗旺，胡力心，陳彬彬，等.三維地質(zhì)建模方法在區(qū)域地質(zhì)調(diào)查中的應(yīng)用研究[J].三峽生態(tài)環(huán)境監(jiān)測，2019，4（3）：74-79.YI ZW，HU LX，CHEN BB，et al.Application of 3D geological modeling method in regional geological survey[J].Ecology and Environmental Monitoring of Three Gorges，2019，4（3）：74-79.（in Chinese）

[26]劉東升，劉漢龍，吳越，等.巖土材料的基因特征及其測試方法[J].土木與環(huán)境工程學(xué)報（中英文），2022，44（4）：1-9.LIU D，LIU HL，WU Y，et al.Genetic features of geo-materials and their testing method[J].Journal of Civil and Environmental Engineering，2022，44（4）：1-9.

（編輯王秀玲）