考慮土性空間變異特征的智能地層劃分方法

摘" 要：同一土層中的巖土參數(shù)會(huì)隨著位置改變而變化，在一定距離范圍內(nèi)存在一定程度的相關(guān)關(guān)系，該相關(guān)性隨著距離的增大而衰減，直至不相關(guān)。地層劃分是巖土工程勘察設(shè)計(jì)的重要內(nèi)容，在地層劃分中應(yīng)考慮土性的空間相關(guān)性，且有必要提高地層劃分的智能化程度。該研究發(fā)展一套能“自然”地考慮土性空間變異特征的智能地層劃分方法。首先分析某場(chǎng)地巖土工程勘察數(shù)據(jù)的空間變異性特征，并據(jù)此生成一個(gè)隨機(jī)場(chǎng)實(shí)例，構(gòu)建“完備地層”模型；對(duì)構(gòu)建的地層模型進(jìn)行6個(gè)孔靜力觸探模擬，基于靜力觸探結(jié)果展示基于貝葉斯壓縮感知的數(shù)據(jù)擴(kuò)展和基于支持向量機(jī)、高斯模型和隱馬爾科夫隨機(jī)場(chǎng)的地層智能劃分方法，將劃分結(jié)果與“完備地層”模型對(duì)比，論證地層劃分方法的準(zhǔn)確性。

關(guān)鍵詞：土性變異性；地層劃分；相關(guān)距離；隨機(jī)場(chǎng)理論；靜力觸探

中圖分類(lèi)號(hào)：TU43" " " 文獻(xiàn)標(biāo)志碼：A" " " " " 文章編號(hào)：2095-2945（2024）36-0139-06

Abstract： Geotechnical parameters in the same soil layer will change with the change of location. There is a certain degree of correlation within a certain distance. This correlation will attenuate with the increase of distance until it is irrelevant. Stratum division is an important part of geotechnical engineering survey and design. The spatial correlation of soil properties should be considered in stratum division， and it is necessary to improve the intelligence level of stratum division. This study developed a set of intelligent stratigraphic classification methods that can \"naturally\" consider the spatial variation characteristics of soil properties. First， the spatial variability characteristics of geotechnical engineering survey data of a certain site are analyzed， and a random field example is generated based on this to build a \"complete stratum\" model; Six hole static penetration simulations were carried out on the constructed stratum model. Based on the static penetration results， data expansion based on Bayesian compressive sensing and intelligent stratum classification based on support vector machines， Gaussian models and hidden Markov random fields were demonstrated. The classification results were compared with the \"complete stratum\" model to demonstrate the accuracy of the stratum classification method.

Keywords： soil variability; stratum division; correlation distance; random field theory; static sounding

在巖土工程活動(dòng)中，準(zhǔn)確獲取地層分布是優(yōu)化工程設(shè)計(jì)、保障工程安全、提高工程質(zhì)量的關(guān)鍵一環(huán)。在實(shí)際巖土工程中，地層劃分多依賴(lài)于鉆探取樣和靜力觸探試驗(yàn)。鉆探取樣可以直觀地獲取該鉆孔的土層分布情況；靜力觸探試驗(yàn)獲取的數(shù)據(jù)，雖不能直觀提供土樣分層，但可以提供在深度方向近乎連續(xù)的數(shù)據(jù)點(diǎn)，參數(shù)沿深度變化剖面也有利于描述地下土體分層情況，并且根據(jù)所獲取的參數(shù)還能反映不同土層的力學(xué)特性。以上2種方法都只針對(duì)單一探測(cè)鉆孔獲取沿深度方向的地層數(shù)據(jù)，而真實(shí)的地層結(jié)構(gòu)往往在水平方向同樣呈現(xiàn)出變異性。由于時(shí)間和經(jīng)費(fèi)限制，不論是鉆探取樣還是靜力探測(cè)的數(shù)量都很有限且分布稀疏，很難獲取未知區(qū)域地層的完整分布情況，這是一些嚴(yán)重工程事故的原因之一。而巖土材料的空間相關(guān)性、土性分類(lèi)的不確定性和土層邊界復(fù)雜性又是地層劃分中不可忽視的關(guān)鍵因素，它們導(dǎo)致了地層劃分剖面不可避免的不確定性[1]。有必要在充分考慮實(shí)際巖土參數(shù)的隨機(jī)性及空間相關(guān)性的基礎(chǔ)上對(duì)地層進(jìn)行劃分。

Hegazy等[2]將聚類(lèi)分析應(yīng)用在地層界面分層的問(wèn)題中，該方法能夠利用聚類(lèi)分析測(cè)試出數(shù)據(jù)點(diǎn)之間的細(xì)微差別，從而發(fā)現(xiàn)地層類(lèi)型或者性質(zhì)的劇烈變化這一特性，客觀地劃分土壤剖面中相似的數(shù)據(jù)組。夏艷華等[3]引入地層層序序列和水平集理論來(lái)解決地層界面稀疏特性，實(shí)現(xiàn)了地質(zhì)解釋和空間分割等技術(shù)。Ching等[4]利用小波分析法可以精確地識(shí)別剖面中變化點(diǎn)的位置，可以生成可視化和直觀方法的圖形和圖表，并指出依據(jù)貝葉斯方法劃分土層更加嚴(yán)謹(jǐn)，但由于貝葉斯方法的局限性，其基礎(chǔ)理論和貝葉斯方程計(jì)算復(fù)雜，導(dǎo)致其很難在實(shí)際工程中應(yīng)用。Wang等[5-6]采用期望最大化算法來(lái)修正已經(jīng)建立的隱馬爾可夫模型，此方法提供了相關(guān)突出特性的合理估計(jì)，能夠更加精準(zhǔn)地識(shí)別土層分布情況，更加精準(zhǔn)地確定不同土層之間的邊界劃分情況。曹子君等[7]提出了一種基于靜力觸探試驗(yàn)數(shù)據(jù)自動(dòng)劃分土層的方法，此方法客觀反映了不同土層的統(tǒng)計(jì)特性，能夠合理地反映識(shí)別出的土層邊界的可靠性，不僅能劃分土層剖面，而且定量地表示了土層剖面的不確定性。許程[8]基于概率論統(tǒng)計(jì)分析同一土層力學(xué)參數(shù)之間的線性相關(guān)關(guān)系，并基于隨機(jī)場(chǎng)理論計(jì)算不同土層的相關(guān)距離，根據(jù)各類(lèi)資料總結(jié)鎮(zhèn)江地區(qū)不同地層年代所形成的土層分布規(guī)律。Molina-Gómez等[9]提出了一種基于聚類(lèi)分析的多元統(tǒng)計(jì)方法，能夠客觀地對(duì)土壤剖面中相似的數(shù)據(jù)進(jìn)行分類(lèi)。Suryasentana等[10]研究了一種概率方法，通過(guò)使用貝葉斯檢測(cè)方法處理動(dòng)力觸探數(shù)據(jù)，檢測(cè)不同土層之間的土壤數(shù)據(jù)趨勢(shì)的突變來(lái)圈定不同土層，實(shí)現(xiàn)了更可靠的土壤地層識(shí)別方案。

上述研究雖提出了眾多地層劃分方法，但其普適性驗(yàn)證不充分，且劃分結(jié)果無(wú)“標(biāo)準(zhǔn)答案”可對(duì)照，對(duì)地層本身存在的空間變異性特征也未充分考慮。為此，本研究首先基于某場(chǎng)地巖土工程勘察數(shù)據(jù)分析了各層土的土性空間變異性量化指標(biāo)，并據(jù)此生成該場(chǎng)地的一個(gè)隨機(jī)場(chǎng)實(shí)例，作為“完備地層”模型。對(duì)該“完備地層”，模擬其靜力觸探過(guò)程，基于靜力觸探結(jié)果進(jìn)行地層智能劃分方法探索，將劃分結(jié)果與“完備地層”模型對(duì)比，探討了地層劃分的準(zhǔn)確性。

1" 考慮空間變異性的“完備地層”模型構(gòu)建

1.1" 土性參數(shù)空間變異性分析

對(duì)某場(chǎng)地進(jìn)行的巖土工程勘查，共實(shí)施了14個(gè)鉆孔取樣和6個(gè)靜力觸探，識(shí)別到六層土（自上而下分別為填土1、黏土1、粉質(zhì)黏土1、粉土粉砂、粉質(zhì)黏土2和黏土2），其彈性模量、黏聚力、內(nèi)摩擦角均值和標(biāo)準(zhǔn)差見(jiàn)表1。

土性相關(guān)距離是土性剖面隨機(jī)場(chǎng)模型的重要度量。在地質(zhì)變化與土體沉積過(guò)程中，土層因地質(zhì)、環(huán)境、物理、化學(xué)等因素形成了土體天然的變異性。對(duì)于土體內(nèi)的任意兩點(diǎn)，在一定距離內(nèi)存在一定的相關(guān)性，大于此距離則可認(rèn)為這兩點(diǎn)之間不存在相關(guān)性，這個(gè)距離就稱(chēng)之為相關(guān)距離。因測(cè)試上述力學(xué)指標(biāo)的試樣的取樣深度已知，此處采用相關(guān)函數(shù)法、基于線指數(shù)型相關(guān)函數(shù)擬合得到6種土體的垂直相關(guān)距離，各個(gè)參數(shù)的相關(guān)距離非常接近（后文取其均值），擬合結(jié)果見(jiàn)表1。由于各孔間距較大，無(wú)法直接獲得水平方向的相關(guān)距離，在不失一般性的情況下，依據(jù)前人工作[11]，取10倍的垂直相關(guān)距離作為水平相關(guān)距離。

1.2" “完備地層”模型構(gòu)建

隨機(jī)變量的不可數(shù)無(wú)窮集構(gòu)成了隨機(jī)場(chǎng)，可用一系列離散的隨機(jī)數(shù)近似描述。隨機(jī)場(chǎng)處理離散問(wèn)題的方法有：局部平均法、中心點(diǎn)法、Karhunen-Loeve（K-L）展開(kāi)法等。其中，K-L展開(kāi)法由于優(yōu)越的解析性而在隨機(jī)場(chǎng)分析中廣泛應(yīng)用，其實(shí)質(zhì)是將隨機(jī)場(chǎng)分解為一系列不相關(guān)的隨機(jī)變量和確定系數(shù)（特征函數(shù)、特征值）；與其他隨機(jī)場(chǎng)離散方法相比，其對(duì)于任意類(lèi)型的隨機(jī)場(chǎng)均收斂，展開(kāi)階數(shù)相同時(shí)具有最小的均方誤差。既有文獻(xiàn)表明，K-L展開(kāi)對(duì)規(guī)則的矩形區(qū)域有較好的收斂性能，對(duì)不規(guī)則區(qū)域的收斂性不佳。參考前人經(jīng)驗(yàn)，對(duì)典型地質(zhì)剖面內(nèi)的各層土體，先單獨(dú)在整個(gè)剖面計(jì)算其K-L展開(kāi)結(jié)果，再根據(jù)指定的地層分界面對(duì)各矩形區(qū)域進(jìn)行剖切合并，由此構(gòu)成二維的地層剖面。二維剖面數(shù)據(jù)點(diǎn)維度為50×110，即由5 500個(gè)像素點(diǎn)組成矩陣，結(jié)果如圖1所示。圖中橫向黑虛線表示地層分界線，縱向黑虛線表示后續(xù)模擬靜力觸探孔位置。

圖1是基于實(shí)際地層有限的地勘數(shù)據(jù)、指定的地層界面，利用隨機(jī)場(chǎng)理論“擴(kuò)展”出的具有完備信息（包括彈性模量、黏聚力和內(nèi)摩擦角）的“完備地層”模型。這里的“完備”是指地層模型具有后續(xù)數(shù)值模擬所需的所有幾何與力學(xué)參數(shù)。其理論價(jià)值在于地層界面位置完全已知，力學(xué)參數(shù)完全已知；可基于該地層模型模擬地勘過(guò)程并獲取地勘參數(shù)（如靜力觸探曲線），然后基于地勘數(shù)據(jù)進(jìn)行地層劃分，將該劃分結(jié)果與完全已知的“正確”地層界面位置對(duì)比，實(shí)現(xiàn)從理論上探討基于靜力觸探曲線的智能地層劃分方法的可靠性。

2" 靜力觸探過(guò)程的數(shù)值模擬

本節(jié)基于“完備地層”模型，提取各個(gè)靜力觸探孔附近的地層剖面數(shù)據(jù)，鏈接到ABAQUS軟件中，將隨機(jī)場(chǎng)理論與有限元分析相結(jié)合，實(shí)現(xiàn)靜力觸探的有限元模擬。模擬結(jié)果將用作后續(xù)地層智能劃分的依據(jù)。

本節(jié)共模擬6個(gè)孔（孔7—孔12）的靜力觸探，孔位如圖1所示。此處采用耦合歐拉-拉格朗日法模擬，將土體視作歐拉體，網(wǎng)格保持不變，材料可在網(wǎng)格內(nèi)發(fā)生任意流動(dòng)；將觸探桿、錐設(shè)為拉格朗日體，材料變形與網(wǎng)格變形保持一致。該方法的優(yōu)點(diǎn)是網(wǎng)格不發(fā)生扭曲，并能準(zhǔn)確地計(jì)算材料邊界，有效解決了有關(guān)大變形、材料破壞和流體材料等諸多問(wèn)題。

土體模擬范圍是以探頭貫入路徑為中心、邊長(zhǎng)為0.5 m和深度為25 m的長(zhǎng)方體，模型底部為固定端，四周為法向約束，頂端為自由端。探頭采用變形體模擬，錐角為60°、錐底直徑為3.57 cm、探桿長(zhǎng)度30 m，貫入速度為20 mm/s。貫入前，探頭錐尖位于土體頂面處。土體采用摩爾-庫(kù)倫本構(gòu)模型描述參數(shù)按照?qǐng)D1賦值。采用ABAQUS內(nèi)置的通用接觸模擬土體與觸探錐和桿之間的相互作用，接觸算法為歐拉-拉格朗日接觸。切向接觸摩擦系數(shù)為0.2，接觸面法向行為定義為“硬”接觸。另對(duì)探頭變形體施加剛體約束，限制探頭本身在貫入的變形。

探頭、探桿采用八結(jié)點(diǎn)線性六面體的單元（C3D8R），減縮積分，沙漏控制。土體采用八節(jié)點(diǎn)線性歐拉六面體單元（EC3D8R）劃分網(wǎng)格。為減少計(jì)算時(shí)間和提高計(jì)算精度，探桿貫入范圍附近采用較密的網(wǎng)格劃分方式，外周采用相對(duì)稀疏的劃分方式。網(wǎng)格劃分如圖2所示。試算表明，圖中的網(wǎng)格劃分密度能得出合理的模擬結(jié)果。模型采用顯示動(dòng)力學(xué)求解器進(jìn)行求解計(jì)算，時(shí)間步長(zhǎng)0.000 1 s，并開(kāi)啟大變形選項(xiàng)。

從模擬結(jié)果中提取比貫入阻力沿深度的變化，并與地勘報(bào)告靜力觸探曲線對(duì)比，如圖3所示。由圖3可知，基于隨機(jī)場(chǎng)理論考慮彈性模量、黏聚力和內(nèi)摩擦角3種巖土參數(shù)的空間變異性，靜力觸探貫入模型能夠較好地模擬錐尖與土體之間的相互作用，比貫入阻力隨深度的變化趨勢(shì)與原位該處?kù)o力貫入試驗(yàn)結(jié)果基本一致，表明此處的數(shù)值模擬方法合理、可靠。

3" 地層智能劃分方法

本研究提出的智能地層劃分方法包括2個(gè)主要步驟：①基于貝葉斯壓縮感知（Bayesian Compressive Sensing，BCS）方法，對(duì)有限的地勘數(shù)據(jù)進(jìn)行空間擴(kuò)充，使得場(chǎng)地剖面各處均有擴(kuò)充的勘察數(shù)據(jù)；②基于數(shù)據(jù)聚類(lèi)識(shí)別各土類(lèi)的判別邊界，據(jù)此對(duì)地層剖面各點(diǎn)判別其土類(lèi)，形成地層劃分結(jié)果。以下對(duì)上述步驟逐一展開(kāi)介紹。

3.1" 基于貝葉斯壓縮感知的地層數(shù)據(jù)擴(kuò)充

BCS是一種廣泛應(yīng)用于雷達(dá)信號(hào)成像、人臉識(shí)別、無(wú)線通信、機(jī)器學(xué)習(xí)等領(lǐng)域的信號(hào)處理技術(shù)，其最大的優(yōu)點(diǎn)是以采集到的稀疏信號(hào)進(jìn)行插值和展開(kāi)，實(shí)現(xiàn)對(duì)原始數(shù)據(jù)的重構(gòu)。在巖土工程勘察中，靜力觸探數(shù)據(jù)可認(rèn)為是對(duì)地層的一次稀疏采樣，本研究采用Fang等[12]發(fā)展的BCS算法，基于少量觸探數(shù)據(jù)（此處為第2節(jié)模擬的靜力觸探結(jié)果）完成對(duì)地層剖面數(shù)據(jù)擴(kuò)充。

勘察剖面長(zhǎng)度為110 m、深度為25，靜力觸探孔共6個(gè)，從左到右各孔水平間隔分別為23、21、24、23、19 m。垂直分辨率為0.5 m，水平分辨率為1.0 m，則待擴(kuò)充的場(chǎng)地?cái)?shù)據(jù)維度為110×50，而已知稀疏數(shù)據(jù)為6個(gè)鉆孔模擬得到的比貫入阻力，維度為6×50，經(jīng)BCS算法擴(kuò)充后的比貫入阻力分布如圖4所示。由圖4可知，擴(kuò)充的結(jié)果水平向呈現(xiàn)帶狀分布特征，比貫入阻力總體沿深度方向的變化趨勢(shì)和對(duì)應(yīng)的模擬數(shù)據(jù)變化趨勢(shì)基本一致。

3.2" 智能地層劃分方法

綜合鉆孔取樣目測(cè)、物理力學(xué)特性測(cè)試等結(jié)果，可將土類(lèi)與靜力觸探比貫入阻力進(jìn)行關(guān)聯(lián)，如圖5所示，該“深度-比貫入阻力”二維平面內(nèi)以點(diǎn)的顏色區(qū)分不同土類(lèi)。為更加準(zhǔn)確地在圖5中劃分不同土類(lèi)的分界面，實(shí)現(xiàn)數(shù)據(jù)聚類(lèi)，首先采用基于線性核函數(shù)的支持向量機(jī)（Support Vector Machine，SVM）算法，獲得圖5中各土類(lèi)的線性邊界。SVM能尋找出邊界線使各土類(lèi)數(shù)據(jù)點(diǎn)盡可能地分布在邊界的兩側(cè)，結(jié)果如圖5所示（各聚類(lèi)分區(qū)以色塊顏色區(qū)分）。在上述擴(kuò)充的地層靜力觸探比貫入阻力（圖4）基礎(chǔ)上，利用SVM聚類(lèi)結(jié)果可判斷勘察斷面上各點(diǎn)所屬的土類(lèi)，結(jié)果如圖6所示?？梢钥闯?，本方法推測(cè)的二維剖面的土層分布基本符合原始的土層分布（黑虛線為原始邊界）。第2層土（黏土1）和第3層土（粉質(zhì)黏土1）相對(duì)于其他土層的估計(jì)準(zhǔn)確度略低。是由于該土層相對(duì)其他土層面積偏小，對(duì)插值結(jié)果最為敏感。且在進(jìn)行數(shù)據(jù)統(tǒng)計(jì)時(shí)，該土層數(shù)量也相對(duì)較少，準(zhǔn)確度也受到統(tǒng)計(jì)不確定的影響。綜上所述，總體地層的分布較為清晰，分類(lèi)較為準(zhǔn)確，但土層突變處劃分粗糙，與原始地層分布還存在一定的誤差。經(jīng)計(jì)算整個(gè)研究場(chǎng)地的二維剖面的5 500個(gè)像素點(diǎn)分類(lèi)準(zhǔn)確的個(gè)數(shù)為5 203個(gè)，準(zhǔn)確率為94.6%。

在SVM聚類(lèi)基礎(chǔ)上，進(jìn)一步采用高斯模型（Gaussian Model， GM）算法對(duì)上述聚類(lèi)進(jìn)行概率運(yùn)算，其實(shí)質(zhì)是給出圖5中各數(shù)據(jù)點(diǎn)屬于某類(lèi)土的概率。樣本數(shù)據(jù)x（深度，比貫入阻力）是多維數(shù)據(jù)，高斯分布用下列概率密度函數(shù)描述

式中：u為期望；？撞為協(xié)方差；D為數(shù)據(jù)維度。高斯模型參數(shù)？茲=（？滋，？撞）由期望和協(xié)方差組成。對(duì)于高斯模型中參數(shù)？茲的值可以通過(guò)極大似然函數(shù)估計(jì)

假設(shè)每個(gè)數(shù)據(jù)點(diǎn)是獨(dú)立的，N表示樣本點(diǎn)的總數(shù)。似然函數(shù)L（？茲）可以由概率密度函數(shù)給出

最終結(jié)果如圖5中各等概率橢圓所示：在每個(gè)聚類(lèi)分區(qū)中，距離同族橢圓中心越近則概率越大。將上述擴(kuò)充的地層中各點(diǎn)靜力觸探數(shù)據(jù)帶入即得各點(diǎn)屬于各土類(lèi)的概率，取概率最大值對(duì)應(yīng)的土類(lèi)為改點(diǎn)所屬土類(lèi)，結(jié)果如圖7所示?？梢钥吹?，第4層土的上、下界面的識(shí)別準(zhǔn)確性比圖6有明顯改善，土層劃分準(zhǔn)確性提高到95.6%。

SVM+GM算法對(duì)于地層分界處的細(xì)節(jié)劃分仍有所欠缺，這是由于各個(gè)點(diǎn)之間僅根據(jù)圖5的聚類(lèi)進(jìn)行土類(lèi)確定，相鄰點(diǎn)之間的聯(lián)系并未考慮。為了考慮相鄰點(diǎn)在土類(lèi)屬性上的潛在聯(lián)系以更加符合實(shí)際情況，進(jìn)一步采用隱馬爾可夫隨機(jī)場(chǎng)（Hidden Markov Random Field，HMRF）對(duì)分類(lèi)模型進(jìn)行進(jìn)一步細(xì)化，具體算法詳見(jiàn)Zhao等[13]，細(xì)化結(jié)果如圖8所示。結(jié)果表明，經(jīng)過(guò)HMRF優(yōu)化后，識(shí)別準(zhǔn)確率提高為97.35%，對(duì)土層與土層之間的邊界處像素點(diǎn)劃分更加精細(xì)準(zhǔn)確。第1層土（填土1）和第6層土（黏土2）與原始地層分布基本一致，第2層土（黏土1）上邊界處誤判為填土1的點(diǎn)更少。第4層土（粉土粉砂）的上下分界處有少量位置不準(zhǔn)確。

綜上，SVM+GM+HMRF算法具有最好的地層劃分精度。其中，SVM主要實(shí)現(xiàn)地層初步劃分，而GM、HMRF則逐步提高劃分的精細(xì)程度。BCS數(shù)據(jù)擴(kuò)充和地層劃分方法中涉及的參數(shù)都是基于地勘數(shù)據(jù)自動(dòng)識(shí)別獲取，無(wú)需人為干涉，因而可實(shí)現(xiàn)整個(gè)地層劃分的智能化。此外，土性參數(shù)的空間變異性特征“自然地”體現(xiàn)在BCS數(shù)據(jù)擴(kuò)充和HMRF中。

4" 結(jié)論

本文發(fā)展了基于靜力觸探結(jié)果對(duì)地層進(jìn)行智能劃分的方法，主要結(jié)論如下。

1）基于線指數(shù)型相關(guān)函數(shù)可擬合得到土性參數(shù)空間變異性的相關(guān)距離，運(yùn)用K-L展開(kāi)法能得到合理的隨機(jī)場(chǎng)實(shí)例，考慮參數(shù)空間變異性的數(shù)值模擬能夠較好地模擬錐尖與地基土體之間的相互作用，模擬得到的比貫入阻力分布與原位靜力貫入試驗(yàn)的結(jié)果基本一致。

2）基于貝葉斯壓縮感知理論對(duì)模擬數(shù)據(jù)進(jìn)行擴(kuò)充，可可靠地獲取整個(gè)研究剖面的比貫入阻力分布，其數(shù)值在地勘剖面上呈現(xiàn)帶狀分布，與原模擬數(shù)據(jù)趨勢(shì)吻合，數(shù)據(jù)擴(kuò)充方法可行。

3）采用支持向量機(jī)（SVM）+高斯模型（GM）+隱馬爾可夫隨機(jī)場(chǎng)（HMRF）算法相結(jié)合對(duì)擴(kuò)充后的地層數(shù)據(jù)進(jìn)行地層劃分，模型準(zhǔn)確率可達(dá)97.35%，各地層位置整體準(zhǔn)確性高，部分層位附近少量點(diǎn)位劃分準(zhǔn)確性稍低。驗(yàn)證了本文地層劃分方法的可行性。該方法只需輸入擴(kuò)充的地層數(shù)據(jù)，其余參數(shù)可自動(dòng)從輸入數(shù)據(jù)中學(xué)習(xí)得到，從而支持智能化地層劃分的實(shí)現(xiàn)。

參考文獻(xiàn)：

[1] 鄭碩.基于靜力觸探的土層自動(dòng)劃分方法及不確定性表征[D].武漢：武漢大學(xué)，2017.

[2] HEGAZY Y A， MAYNE P W. Objective site characterization using clustering of piezocone data[J]. Journal of Geotechnical and Geoenvironmental Engineering， 2002，128（12）：986-996.

[3] 夏艷華，白世偉.復(fù)雜三維地層模型地層劃分研究[J].安徽理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，31（3）：35-39.

[4] CHING J， WANG J S， JUANG C H， et al. Cone penetration test （CPT）-based stratigraphic profiling using the wavelet transform modulus maxima method[J]. Canadian Geotechnical Journal， 2015，52（12）：1993-2007.

[5] WANG H， WELLMANN J F， LI Z， et al. A segmentation approach for stochastic geological modeling using hidden Markov random fields[J]. Mathematical Geosciences， 2017，49（2）：145-177.

[6] WANG X R， WANG H， LIANG R Y， et al. A hidden Markov random field model based approach for probabilistic site characterization using multiple cone penetration test data[J]. Structural Safety， 2018，70：128-138.

[7] 曹子君，鄭碩，李典慶，等.基于靜力觸探的土層自動(dòng)劃分方法與不確定性表征[J].巖土工程學(xué)報(bào)，2018，40（2）：336-345.

[8] 許程.基于不確定性理論的鎮(zhèn)江第四紀(jì)土層層序劃分研究[D].南京：東南大學(xué)，2021.

[9] MOLINA-G？魷MEZ F， VIANA D F A， FERREIRA C， et al. Defining the soil stratigraphy from seismic piezocone data： A clustering approach[J]. Engineering Geology， 2021，287：106-111.

[10] SURYASENTANA S K， LAWLER M， SHEIL B B， et al. Probabilistic soil strata delineation using DPT data and Bayesian changepoint detection[J]. Journal of Geotechnical and Geoenvironmental Engineering， 2023，149（4）：6-23.

[11] 吳振君.土體參數(shù)空間變異性模擬和土坡可靠度分析方法應(yīng)用研究[D].武漢：中國(guó)科學(xué)院武漢巖土力學(xué)研究所，2009.

[12] FANG Y， WU J J， HUANG B M. 2D sparse signal recovery via 2D orthogonal matching pursuit[J]. Science China Information Sciences， 2012，55（4）：889-897.

[13] ZHAO T Y， WANG Y. Non-parametric simulation of non-stationary non-gaussian 3d random field samples directly from sparse measurements using signal decomposition and markov chain monte carlo[J]. Reliability Engineering and System Safety，2020，203（1）：107-108.