李治斌，袁立月，丁 黔，張淼鑫，王紅梅

(1.黑龍江大學(xué) 建筑工程學(xué)院，黑龍江 哈爾濱 150080；2.黑龍江大學(xué) 水利電力學(xué)院，黑龍江 哈爾濱 150080)

鹽漬土廣泛分布在我國(guó)不同地區(qū)，總面積約占我國(guó)國(guó)土面積的2%[1]，而鹽漬土本身具有溶陷、鹽脹或腐蝕等性質(zhì)[2]，力學(xué)和工程性質(zhì)不佳，因此，十分有必要對(duì)改良鹽漬土展開(kāi)深入的研究。

近年來(lái)，不少學(xué)者對(duì)凍結(jié)土和改良土單軸試驗(yàn)的數(shù)學(xué)模型進(jìn)行了研究。陳鑫[3]研究了尺寸和加載速率對(duì)凍結(jié)水泥土單軸壓縮影響，發(fā)現(xiàn)凍結(jié)水泥土的破壞應(yīng)變隨溫度的降低和加載速率的增加而增大，在1.94%～6.94% 之間變化，不同加載速率條件下破壞應(yīng)變與溫度呈冪函數(shù)關(guān)系，并建立了溫度和加載速率對(duì)水泥土單軸抗壓強(qiáng)度的數(shù)學(xué)預(yù)測(cè)模型和支持向量機(jī)的預(yù)測(cè)模型。曹智國(guó)[4]根據(jù)水泥土無(wú)側(cè)限抗壓強(qiáng)度試驗(yàn)，提出一個(gè)綜合反映水泥摻入量、養(yǎng)護(hù)齡期以及孔隙率等因素對(duì)水泥土強(qiáng)度影響規(guī)律的綜合表征參數(shù)，水泥加固土和水泥穩(wěn)定土的強(qiáng)度均與綜合表征參數(shù)呈冪函數(shù)關(guān)系。

而影響鹽漬土強(qiáng)度的因素有很多，這些影響因素與強(qiáng)度之間的關(guān)系往往是非線性的，難以用準(zhǔn)確的數(shù)學(xué)模型來(lái)表達(dá)，所以越來(lái)越多的學(xué)者利用神經(jīng)網(wǎng)絡(luò)來(lái)建立所需的預(yù)測(cè)模型。如宋啟卓等[5]利用BP神經(jīng)網(wǎng)絡(luò)研究了鹽漬土的鹽脹特性。張強(qiáng)等[6]利用小波神經(jīng)網(wǎng)絡(luò)建立了秸稈灰混凝土強(qiáng)度預(yù)測(cè)模型。王彥虎[7]、霍珍生[8]等利用BP神經(jīng)網(wǎng)絡(luò)建立了碳酸鹽漬土的凍脹模型，并利用該模型對(duì)影響凍脹的因素進(jìn)行敏感性分析。劉宇峰[9]則分別利用BP神經(jīng)網(wǎng)絡(luò)、自適應(yīng)模糊神經(jīng)推理系統(tǒng)和支持向量機(jī)建立了鹽漬土未凍水含量的預(yù)測(cè)模型。展超[10]基于BP神經(jīng)網(wǎng)絡(luò)建立了適用于富水砂層地區(qū)渣土改良試驗(yàn)的效果預(yù)測(cè)模型，用于土的坍落度、滲透系數(shù)和內(nèi)摩擦角的預(yù)測(cè)，效果良好。這些已有研究表明這些神經(jīng)網(wǎng)絡(luò)模型均能較好地反映出影響因素與預(yù)測(cè)目標(biāo)之間復(fù)雜的關(guān)系，可靠性比較高。

目前，國(guó)內(nèi)學(xué)者對(duì)鹽漬土的研究主要是針對(duì)氯鹽和硫酸鹽，對(duì)于宏觀影響因素影響下改良碳酸鹽漬土單軸壓縮強(qiáng)度的神經(jīng)網(wǎng)絡(luò)模型研究尚比較缺乏，因此，本文采用實(shí)驗(yàn)室配制碳酸鹽漬土的方式，考慮了固化劑摻量、溫度和加載速率對(duì)改良鹽漬土的無(wú)側(cè)限抗壓強(qiáng)度的影響，通過(guò)室內(nèi)無(wú)側(cè)限抗壓試驗(yàn)獲得了土樣的強(qiáng)度值。最后，利用BP神經(jīng)網(wǎng)絡(luò)建立了固化碳酸鹽漬土無(wú)側(cè)限抗壓強(qiáng)度的預(yù)測(cè)模型。

1 研究方法

1.1 樣本數(shù)據(jù)收集和預(yù)處理

試驗(yàn)用土來(lái)自黑龍江省哈爾濱市機(jī)場(chǎng)路第2通道的粉質(zhì)黏土。室內(nèi)進(jìn)行試驗(yàn)前，在試驗(yàn)土樣中加入1.5%的碳酸氫鈉固體顆粒，配制成人工鹽漬土。改良鹽漬土的無(wú)側(cè)限壓縮試驗(yàn)考慮石灰和粉煤灰摻量、溫度和加載速率3個(gè)因素的影響。其中，石灰和粉煤灰總摻量保持15%不變，設(shè)計(jì)4組摻量組合：3%石灰+12%粉煤灰、5%石灰+10%粉煤灰、6%石灰+9%粉煤灰、7%石灰+8%粉煤灰；試驗(yàn)溫度設(shè)置4組：0 ℃、-5 ℃、-10 ℃、-15 ℃；加載速率設(shè)置3組：0.5 mm/min、1 mm/min和2 mm/min。

試驗(yàn)前，將試樣在標(biāo)準(zhǔn)養(yǎng)護(hù)室養(yǎng)護(hù)28 d，然后將4組溫度的樣品放入凍箱中冷凍24 h。試樣的無(wú)側(cè)限抗壓強(qiáng)度采用WDW-100微機(jī)控制電子萬(wàn)能試驗(yàn)機(jī)測(cè)定。

1.2 模型結(jié)構(gòu)及樣本數(shù)據(jù)預(yù)處理

各因素與固化土無(wú)側(cè)限抗壓強(qiáng)度之間的關(guān)系難以用精確的數(shù)學(xué)模型來(lái)表達(dá)，而B(niǎo)P人工神經(jīng)網(wǎng)絡(luò)具有很強(qiáng)的自學(xué)習(xí)和自適應(yīng)能力，能很好地反映出這種復(fù)雜的非線性關(guān)系。因此，本文用Python語(yǔ)言編寫(xiě)了BP人工神經(jīng)網(wǎng)絡(luò)，嘗試建立二灰固化土的無(wú)側(cè)限抗壓強(qiáng)度預(yù)測(cè)模型。BP神經(jīng)網(wǎng)絡(luò)一般包括輸入層、隱含層和輸出層，采用誤差反向傳播學(xué)習(xí)算法，其學(xué)習(xí)過(guò)程主要包括正向傳播和反向傳播兩部分，在正向傳播中，輸入信息從輸入層經(jīng)過(guò)權(quán)值的處理傳到隱含層，在隱含層經(jīng)過(guò)激活函數(shù)處理傳向輸出層，最終輸出預(yù)測(cè)值。

首先需要確定BP神經(jīng)網(wǎng)絡(luò)各層的層數(shù)。本文確定神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)為三層，包括一個(gè)輸入層、一個(gè)隱含層和一個(gè)輸出層。在本次改良土的無(wú)側(cè)限壓縮試驗(yàn)中，考慮了改良鹽漬土的固化劑摻量、溫度、加載速率作為影響土樣強(qiáng)度大小的主要因素，因此選擇這3個(gè)因素的數(shù)據(jù)導(dǎo)入模型中，對(duì)這三種因素的改良鹽漬土強(qiáng)度預(yù)測(cè)。神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)，如圖1。

圖1 BP神經(jīng)網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)

為提高和保證預(yù)測(cè)精度，使模型訓(xùn)練的樣本數(shù)據(jù)要有較好的代表性，在開(kāi)始訓(xùn)練預(yù)測(cè)模型之前，要對(duì)無(wú)側(cè)限抗壓強(qiáng)度值及各影響因素?cái)?shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。本文從48個(gè)樣本數(shù)據(jù)中隨機(jī)選取40個(gè)作為訓(xùn)練樣本，剩下的8個(gè)作為預(yù)測(cè)樣本。將所有的數(shù)據(jù)參數(shù)X按公式(1)進(jìn)行歸一化處理。

(1)

式中：Xmin是參數(shù)X中的最小值；Xmax是參數(shù)X中的最大值；Xn是參數(shù)X歸一化后的值。

2 基于BP神經(jīng)網(wǎng)絡(luò)強(qiáng)度預(yù)測(cè)模型的建立

2.1 模型結(jié)構(gòu)參數(shù)選取

神輸入層3個(gè)節(jié)點(diǎn)，分別為固化劑摻量(本文用石灰摻量代替固化劑摻量參與模型計(jì)算)，溫度和加載速率；隱含層節(jié)點(diǎn)數(shù)待定；輸出層1個(gè)節(jié)點(diǎn)為無(wú)側(cè)限抗壓強(qiáng)度。每層的輸入數(shù)據(jù)通過(guò)上一層的輸出數(shù)據(jù)經(jīng)sigmoid函數(shù)獲得。

根據(jù)式(2)可確定[11]隱含層節(jié)點(diǎn)的數(shù)量范圍為2～12。

(2)

式中：n是隱含層的神經(jīng)元個(gè)數(shù)；ni是輸入層的神經(jīng)元個(gè)數(shù)；n0是輸出層的神經(jīng)元個(gè)數(shù)；a是一個(gè)0～10之間的固定值。

用訓(xùn)練樣本對(duì)不同隱含層節(jié)點(diǎn)的模型進(jìn)行訓(xùn)練，可得到不同隱含層節(jié)點(diǎn)模型訓(xùn)練后結(jié)果的均方根誤差(RMSE)和確定系數(shù)(R2)，如圖2所示。從圖中可以看出，當(dāng)隱含層的神經(jīng)元個(gè)數(shù)n=11時(shí)，數(shù)據(jù)訓(xùn)練后RMSE達(dá)到最小值，且R2達(dá)到最大值，所以本次隱含層神經(jīng)元個(gè)數(shù)設(shè)置為11。

經(jīng)過(guò)計(jì)算統(tǒng)計(jì)，本次神經(jīng)網(wǎng)絡(luò)模型的參數(shù)值選擇如下：輸入層神經(jīng)元個(gè)數(shù)為3；隱含層=1；

圖2 不同隱含層節(jié)點(diǎn)模型訓(xùn)練后結(jié)果

隱含層神經(jīng)元個(gè)數(shù)=11；輸出層個(gè)數(shù)=1；模擬精度=0.001；學(xué)習(xí)步數(shù)=10 000；學(xué)習(xí)算法為梯度下降法。

2.2 訓(xùn)練模型誤差分析

當(dāng)隱含層節(jié)點(diǎn)數(shù)為11時(shí)，無(wú)側(cè)限抗壓強(qiáng)度訓(xùn)練值和實(shí)測(cè)值如圖3所示。由圖知，訓(xùn)練值和實(shí)測(cè)值趨勢(shì)基本相同，誤差最小為0.0009，可以反映三因素對(duì)改良土無(wú)側(cè)限抗壓強(qiáng)度的影響；其中訓(xùn)練樣本編號(hào)為3、28、35、39的誤差較大，而這些誤差均發(fā)生在試驗(yàn)設(shè)置的邊界條件下，即固化劑摻量為7%(石灰)、溫度為0 ℃或-15 ℃、加載速率為0.5 mm/min或2 mm/min的條件下，模型的訓(xùn)練誤差較大，主要原因是模型收斂過(guò)程中為保持連續(xù)性，變化幅度較大處誤差變大，這說(shuō)明還有其他因素影響改良土無(wú)側(cè)限抗壓強(qiáng)度。

圖3 模型訓(xùn)練值和實(shí)測(cè)值對(duì)比

2.3 模型預(yù)測(cè)結(jié)果分析

將選中的預(yù)測(cè)樣本歸一化后輸入到建立好的模型中進(jìn)行預(yù)測(cè)。無(wú)側(cè)限抗壓強(qiáng)度預(yù)測(cè)值和實(shí)測(cè)值的對(duì)比如圖4所示。由圖4可以看出，兩條曲線基本重合，說(shuō)明模型預(yù)測(cè)基本準(zhǔn)確。最小誤差為0.0003，最大誤差為0.0599，所有預(yù)測(cè)結(jié)果誤差均小于0.1；且均方根誤差為0.0270，確定系數(shù)為0.985，說(shuō)明模型基本能反映改良土的無(wú)側(cè)限抗壓強(qiáng)度變化趨勢(shì)，也能很好地預(yù)測(cè)無(wú)側(cè)限抗壓強(qiáng)度值。

圖4 模型預(yù)測(cè)值和實(shí)測(cè)值對(duì)比

3 結(jié) 論

本文通過(guò)BP神經(jīng)網(wǎng)絡(luò)建立了改良鹽漬土的無(wú)側(cè)限抗壓強(qiáng)度預(yù)測(cè)模型，研究結(jié)論如下：

(1)無(wú)側(cè)限抗壓試驗(yàn)選擇3個(gè)影響因素：固化劑摻量、溫度和加載速率。影響因素作為輸入層，改良土的無(wú)側(cè)限抗壓強(qiáng)度作為輸出層，建立了BP神經(jīng)網(wǎng)絡(luò)無(wú)側(cè)限抗壓強(qiáng)度預(yù)測(cè)模型。

(2)從無(wú)側(cè)限抗壓試驗(yàn)獲得的48個(gè)樣本數(shù)據(jù)中隨機(jī)選取40個(gè)作為訓(xùn)練樣本，剩下的8個(gè)作為預(yù)測(cè)樣本。當(dāng)隱含層的神經(jīng)元個(gè)數(shù)n=11時(shí)，數(shù)據(jù)訓(xùn)練后RMSE達(dá)到最小值，且R2達(dá)到最大值，說(shuō)明此時(shí)模型收斂精度最高，所以本次隱含層神經(jīng)元個(gè)數(shù)設(shè)置為11。

(3)訓(xùn)練模型在試驗(yàn)設(shè)置的邊界條件處為保持收斂，會(huì)有較大的誤差，但總體訓(xùn)練值與實(shí)測(cè)值誤差不大，可以反映出各影響因素對(duì)改良土無(wú)側(cè)限抗壓強(qiáng)度的影響趨勢(shì)。

(4)預(yù)測(cè)值和實(shí)驗(yàn)值的最小誤差為0.0003，最大誤差為0.0599，所有預(yù)測(cè)結(jié)果誤差均小于0.1；且均方根誤差為0.0270，確定系數(shù)為0.985，說(shuō)明該強(qiáng)度預(yù)測(cè)模型預(yù)測(cè)效果較好，能很好地反映改良土的無(wú)側(cè)限抗壓強(qiáng)度變化趨勢(shì)。