基于深度學(xué)習(xí)的黃土干濕循環(huán)損傷分析

(西安科技大學(xué) 建筑與土木工程學(xué)院，西安 710054)

1 研究背景

黃土作為特殊土的一種，廣泛分布于我國(guó)西北、華北和東北地區(qū)。黃土的組成成分多樣，含有多種氧化物，以SiO2為主，然后是Al2O3、CaO，其次為Fe2O3、MgO和K2O等。黃土顆粒之間較為疏松，含有許多可溶性物質(zhì)，使其具有多孔性[1]、垂直節(jié)理[2]和濕陷性[3]等特殊性質(zhì)。在工程建設(shè)中黃土往往容易被流水侵蝕，破壞其原有的結(jié)構(gòu)，導(dǎo)致塌方[4]、滑坡[5]和不均勻沉降[6]等各種問(wèn)題。當(dāng)氣候處于旱季時(shí)，黃土較為干燥，而氣候處于雨季時(shí)黃土較為濕潤(rùn)，旱-雨季類型的干濕循環(huán)對(duì)黃土造成循環(huán)性風(fēng)化損傷破壞。在這個(gè)干濕循環(huán)過(guò)程中，黃土原有顆粒的尺寸、位置和組成成分均發(fā)生改變，宏觀上導(dǎo)致黃土的抗剪強(qiáng)度、滲透性和壓縮模量產(chǎn)生改性[7]。

干濕循環(huán)條件下的土壤微觀結(jié)構(gòu)變化規(guī)律，可通過(guò)掃描電鏡、X射線衍射和核磁共振技術(shù)等技術(shù)從物理層面直觀描述。Liu等[8]試驗(yàn)研究了干濕循環(huán)對(duì)壓實(shí)低塑性黏土基質(zhì)吸力和體積特性的影響，發(fā)現(xiàn)隨著循環(huán)次數(shù)的增加，干濕循環(huán)下土壤保水曲線(SWRC)的位置向下移動(dòng)。Chen等[9]進(jìn)行直接剪切試驗(yàn)以研究多個(gè)干濕循環(huán)對(duì)不飽和壓實(shí)黏土剪切行為的影響。Chu等[10]通過(guò)掃描電鏡和壓汞法研究了干濕循環(huán)對(duì)改性膨脹土的膨脹行為和可壓縮性的影響，膨脹土微觀結(jié)構(gòu)表現(xiàn)出退化趨勢(shì)。Liu等[11]通過(guò)電鏡掃描發(fā)現(xiàn)在粉質(zhì)黏土干燥/濕潤(rùn)循環(huán)過(guò)程中，土壤顆粒被拉近，導(dǎo)致微觀結(jié)構(gòu)的變化，從而導(dǎo)致力學(xué)行為的改變。

葉萬(wàn)軍等[12]通過(guò)干濕循環(huán)試驗(yàn)揭示細(xì)微觀結(jié)構(gòu)變化與宏觀力學(xué)性質(zhì)的相關(guān)關(guān)系，并區(qū)分古土壤與黃土間存在的差異性。趙貴濤等[13]通過(guò)壓汞試驗(yàn)研究?jī)鋈谘h(huán)、干濕循環(huán)和交替的干濕-凍融循環(huán)對(duì)壓實(shí)膨脹土的微觀結(jié)構(gòu)、土-水特征曲線和收縮特征曲線的影響。胡長(zhǎng)明等[14]通過(guò)電鏡掃描從微觀角度研究干濕循環(huán)引起的壓實(shí)黃土強(qiáng)度劣化，同時(shí)采用ABAQUS數(shù)值模擬分析干濕循環(huán)作用下壓實(shí)黃土填方邊坡穩(wěn)定性。萬(wàn)勇等[15]從微觀層次揭示了壓實(shí)黏土在干濕循環(huán)作用下的變形特性和強(qiáng)度衰減內(nèi)在本質(zhì)。袁志輝等[16]采用單軸拉伸法對(duì)干濕循環(huán)條件下黃土抗拉強(qiáng)度進(jìn)行試驗(yàn)研究，多次的干濕循環(huán)作用打破了原狀黃土的原有結(jié)構(gòu)，使得其抗拉結(jié)構(gòu)強(qiáng)度消失。目前土壤干濕循環(huán)損傷處于試驗(yàn)探索階段，其理論尚不完善，基于土壤微觀紋理建立土壤干濕循環(huán)損傷理論迫在眉睫。

本文采用灰度共生矩陣法[17]，利用陜西西安地區(qū)的黃土電鏡掃描圖片提取黃土微觀紋理特征，通過(guò)建立基于長(zhǎng)短期記憶網(wǎng)絡(luò)(LSTM)[18]的深度學(xué)習(xí)時(shí)序回歸預(yù)測(cè)模型描述黃土紋理特征和土壤微觀裂紋與孔隙占比之間的聯(lián)系，依據(jù)干濕循環(huán)損傷模型計(jì)算土壤干濕循環(huán)損傷因子，揭示土壤微觀結(jié)構(gòu)在干濕循環(huán)過(guò)程的變化規(guī)律。

2 試驗(yàn)過(guò)程

2.1 試驗(yàn)材料

試驗(yàn)用土為取自陜西西安地區(qū)的黃土，氣候?yàn)榕瘻貛О霛駶?rùn)大陸性季風(fēng)氣候，在施工季土壤干濕循環(huán)較為頻繁，其基本性質(zhì)如表1所示。

表1 西安試驗(yàn)黃土基本性質(zhì)Table 1 Basic properties of Xi’an loess for experiment

2.2 試驗(yàn)裝置

試樣微觀結(jié)構(gòu)試驗(yàn)裝置采用荷蘭飛納Phenon/pro臺(tái)式掃描電子顯微鏡，如圖1所示，其SEM分辨率優(yōu)于14 nm，放大倍數(shù)區(qū)間為80～130 000，采用背散射電子探測(cè)器。

圖1 掃描電子顯微鏡Fig.1 Scanning electron microscope(SEM)

2.3 試驗(yàn)方案

采用膠頭滴管將水滴滴在試樣的濾紙上以均勻濕潤(rùn)土壤，再將土壤試樣密封靜置48 h，此時(shí)土壤試樣含水率達(dá)到設(shè)定狀態(tài)。將濕潤(rùn)土樣置于恒溫烘箱進(jìn)行干燥，當(dāng)含水率達(dá)到16%可認(rèn)為土壤試樣處于干燥狀態(tài)，視為土壤進(jìn)行了一次干濕循環(huán)。通過(guò)滴入濾紙的水滴數(shù)量控制土樣的含水率，試驗(yàn)土壤含水率控制在16%～24%，將土壤試樣分別編號(hào)為1組、2組，分別進(jìn)行0次、1次、2次、4次和6次干濕循環(huán)試驗(yàn)，用美工刀將試樣修剪成1 cm3的試塊在噴金處理后進(jìn)行電鏡掃描。將試樣分別放大500、2 000、5 000倍，得到3組電鏡掃描照片。

3 試驗(yàn)結(jié)果及分析

3.1 灰度共生矩陣

電鏡掃描圖像灰度在空間位置上的分布反映了土壤試樣的微觀紋理結(jié)構(gòu)，掃描圖像的灰度空間相關(guān)特性可以用來(lái)描述土壤試樣微觀紋理結(jié)構(gòu)在干濕循環(huán)條件下的變化規(guī)律?；叶裙采仃囀菍?duì)土壤電鏡掃描圖像中固定距離的兩像素具有某種灰度進(jìn)行的統(tǒng)計(jì)，可用來(lái)描述土壤試樣圖像灰度的空間相關(guān)特性。

可通過(guò)編程讀取電鏡掃描灰度圖像，可表示為f(x，y)，則反映該圖像灰度一定空間關(guān)系的灰度共生矩陣P為

P(g1，g2)=

(1)

式中：#為求集合中的元素個(gè)數(shù)；S為土壤圖像中一定空間關(guān)系的像素對(duì)集合；等號(hào)右邊分子為滿足一定空間關(guān)系的像素對(duì)灰度分別為g1和g2的個(gè)數(shù)；分母為電鏡掃描圖像的像素對(duì)總個(gè)數(shù)。在求解土壤試樣圖像的灰度共生矩陣時(shí)，設(shè)置兩像素之間的距離為1保持不變，兩像素之間的角度選取0°、45°、90°、135°。

在具體進(jìn)行土壤試樣微觀紋理結(jié)構(gòu)的評(píng)價(jià)時(shí)，選取5個(gè)指標(biāo)作為土壤試樣灰度圖像的紋理特征統(tǒng)計(jì)值，分別為對(duì)比度、相關(guān)性、熵、平穩(wěn)度和二階矩。

其中對(duì)比度CON計(jì)算公式為

(2)

式中：i與j分別為相鄰像素的灰度；P(i，j)為相鄰像素灰度差為|i-j|的像素分布概率?；叶葓D像對(duì)比度反映了土壤試樣灰度圖像的微觀紋理結(jié)構(gòu)深淺的程度，是對(duì)黃土試樣內(nèi)部微觀裂紋和顆?？紫渡顪\程度的描述，灰度圖像對(duì)比度越大，則說(shuō)明微觀裂紋和孔隙深度越深。

相關(guān)性COR的計(jì)算公式為

(3)

其中：

熵ENT的計(jì)算公式為

(4)

土壤試樣電鏡掃描灰度圖像熵描述了土壤試樣微觀結(jié)構(gòu)紋理在二維灰度圖像中的空間分布狀況，當(dāng)土壤試樣微觀裂紋和孔隙數(shù)量增多時(shí)，土壤試樣微觀紋理結(jié)構(gòu)趨于破壞和混亂化，土壤試樣灰度圖像熵變大。

平穩(wěn)度IDM的計(jì)算公式為

(5)

土壤試樣灰度圖像平穩(wěn)度反映了土壤試樣微觀紋理結(jié)構(gòu)的粗糙程度，土壤試樣微觀裂紋與孔隙在二位灰度圖像中截面積越大，灰度連續(xù)性越強(qiáng)，則平穩(wěn)度越大。

二階矩ASM的計(jì)算公式為

(6)

土壤試樣灰度圖像二階矩描述了土壤試樣微觀紋理均勻程度狀態(tài)，二階矩越小則說(shuō)明土壤試樣灰度圖像灰度共生矩陣灰度越接近，結(jié)構(gòu)紋理比較細(xì)致。灰度圖像二階矩隨土壤干濕循環(huán)的變化反映了土壤微觀裂隙和孔隙的變化狀態(tài)，隨著裂隙和孔隙的增大，二階矩也逐漸變大。

3.2 選擇合適的電鏡掃描圖像放大倍數(shù)

電鏡掃描圖像可以將土壤試樣放大不同的倍數(shù)，選取合適的放大倍數(shù)對(duì)基于灰度共生矩陣的土壤試樣微觀結(jié)構(gòu)分析至關(guān)重要。以1組未進(jìn)行干濕循環(huán)的土壤試樣為例，電鏡掃描圖像如圖2所示，土壤呈現(xiàn)團(tuán)粒結(jié)構(gòu)，多為針孔狀孔隙，含鈣質(zhì)條紋及少量鈣質(zhì)結(jié)核。通過(guò)編程求取土壤電鏡掃描圖像的灰度共生矩陣，計(jì)算灰度圖像的對(duì)比度、相關(guān)性、熵、平穩(wěn)度和二階矩，如圖3所示。

圖2 不同放大倍數(shù)的電鏡掃描圖像Fig.2 SEM images with different magnifications

圖3 灰度共生矩陣各指標(biāo)量與放大倍數(shù)的關(guān)系Fig.3 Relations between indices of gray level co-occurrence matrix and magnification

一般來(lái)說(shuō)，灰度圖像指標(biāo)隨圖像放大倍數(shù)的變化規(guī)律呈現(xiàn)一致性，而圖3中指標(biāo)量熵和二階矩出現(xiàn)突變現(xiàn)象，在放大倍數(shù)為5 000倍時(shí)突然增大，是由于圖像放大倍數(shù)過(guò)大，所計(jì)算指標(biāo)過(guò)于局部，使計(jì)算結(jié)果過(guò)小。但若只放大500倍，又顯得放大倍數(shù)過(guò)小，難以將干濕循環(huán)過(guò)程中土壤試樣顆粒和孔隙發(fā)展規(guī)律表現(xiàn)出來(lái)，因此對(duì)干濕循環(huán)條件下土壤微觀結(jié)構(gòu)機(jī)理的研究采用放大倍數(shù)為2 000倍的電鏡圖像。

從圖3可知熵和二階矩的發(fā)展趨勢(shì)較為相似，對(duì)5個(gè)指標(biāo)量相關(guān)性進(jìn)行分析，結(jié)果如表2所示。熵和二階矩相關(guān)性顯著，在進(jìn)行土壤試樣干濕循環(huán)微觀結(jié)構(gòu)機(jī)理研究時(shí)，本文只選取指標(biāo)量熵進(jìn)行分析。

表2 灰度共生矩陣指標(biāo)量之間的相關(guān)關(guān)系

3.3 干濕循環(huán)條件下灰度指標(biāo)變量變化

將1組土壤試樣電鏡掃描圖像放大倍數(shù)設(shè)定為2 000倍，在0、1、2、4、6次干濕循環(huán)條件下，土壤試樣灰度圖像如圖4所示。

圖4 不同干濕循環(huán)次數(shù)下土壤試樣灰度圖像Fig.4 Gray-scale images of loess samples undergone different dry-wet cycles

根據(jù)試樣灰度圖計(jì)算土壤試樣圖像灰度共生矩陣指標(biāo)量，如圖5所示。圖5中土壤試樣指標(biāo)量對(duì)比度隨著干濕循環(huán)次數(shù)的增加從74.186 5減小到37.684 5，呈現(xiàn)反“S”型曲線，減小速度逐漸增大后又逐漸減小。土壤試樣結(jié)構(gòu)紋理逐漸變淺，是由于干濕循環(huán)過(guò)程中土壤試樣水分會(huì)產(chǎn)生遷移和蒸發(fā)現(xiàn)象，不斷沖刷土壤試樣團(tuán)粒致使土壤試樣整體紋理變淺。

圖5 指標(biāo)量隨干濕循環(huán)次數(shù)變化Fig.5 Changes of gray index quantities with the number of dry-wet cycles

灰度共生矩陣指標(biāo)量相關(guān)性從3.56增加到3.81，增加速度逐漸減小，表明土壤試樣的紋理走向隨干濕循環(huán)次數(shù)的增加逐漸明顯，說(shuō)明在干濕循環(huán)過(guò)程中土壤試樣原有微觀裂紋和孔隙在水分遷移的不斷沖刷下不斷擴(kuò)大，與此同時(shí)在水氣的轉(zhuǎn)化作用下土壤試樣又不斷產(chǎn)生新的微觀裂紋和孔隙，使土壤試樣變得更加蓬松，在水分遷移沖刷作用下土壤試樣中形成了沿原有水分遷移通道的紋理走向。

試樣灰度共生矩陣指標(biāo)量熵隨干濕循環(huán)次數(shù)的增加呈現(xiàn)“S”走向，從0.011增加到0.014，表明土壤試樣微觀結(jié)構(gòu)紋理復(fù)雜度增加。在干濕循環(huán)過(guò)程中，在土壤水分的遷移沖刷下試樣微觀裂紋和孔隙數(shù)量和大小增加，土壤試樣團(tuán)粒結(jié)構(gòu)邊緣不斷圓潤(rùn)，顆粒豐度不斷提高，紋理結(jié)構(gòu)復(fù)雜度不斷增大。

隨著干濕循環(huán)次數(shù)的增加，土壤試樣灰度共生矩陣指標(biāo)量平穩(wěn)度不斷增大，這表明隨著干濕循環(huán)次數(shù)的增加，土壤試樣灰度圖像中灰度較為連續(xù)，在圖像不同區(qū)域之間微觀裂紋和孔隙變化逐漸趨于均勻化，這是由于干濕循環(huán)過(guò)程中水分的遷移和蒸發(fā)使得在原試樣中未產(chǎn)生微觀裂紋和孔隙的區(qū)域不斷產(chǎn)生新的微觀裂紋和孔隙，使得土壤試樣微觀裂紋和孔隙在二維灰度圖像的分布中趨于均勻化，土壤試樣不同區(qū)域間都出現(xiàn)或多或少的水分遷移通道，土壤試樣中的針狀孔隙之間彼此連通，出現(xiàn)條形孔隙，如圖4所示。相關(guān)性和平穩(wěn)度沿凸曲線發(fā)展，在5次干濕循環(huán)后達(dá)到極值。

3.4 干濕循環(huán)條件下土壤孔隙面積測(cè)量

對(duì)放大2 000倍的土壤電鏡掃描圖像進(jìn)行偽彩色處理，增強(qiáng)圖像中土壤微觀裂紋和孔隙與土壤團(tuán)粒的對(duì)比顯示效果，如圖6所示。隨著干濕循環(huán)次數(shù)的增加，土壤中原有的微觀裂紋和針狀孔隙不斷受到水分遷移的沖刷，彼此間逐漸連通形成較大的土壤微觀裂紋，土壤紋理更加淺顯。土壤中水分的遷移將土壤團(tuán)粒周圍的分散顆粒沖刷到一起形成新的土壤顆粒。土壤中原有結(jié)構(gòu)的點(diǎn)-點(diǎn)接觸、點(diǎn)-邊接觸、邊-邊接觸、邊-面接觸和面-面接觸在干濕循環(huán)過(guò)程中持續(xù)發(fā)育，土壤團(tuán)粒邊緣由尖銳變得圓潤(rùn)。

圖6 不同干濕循環(huán)次數(shù)下土壤試樣偽彩色圖像Fig.6 Pseudo-color images of loess samples undergone different dry-wet cycles

將偽彩色圖像變?yōu)閳D像分割模式，將灰度分布概率密度最大處定為土壤微觀裂紋和孔隙面積閾值，提取土壤微觀裂紋和孔隙面積，將已提取土壤微觀裂紋和孔隙面積的土壤圖像轉(zhuǎn)換為二值化圖像，如圖7所示，紅色部分為土壤微觀裂紋和孔隙區(qū)域，灰色部分為土壤顆粒。由圖7中紅色部分的發(fā)展趨勢(shì)可以看出，隨著干濕循環(huán)次數(shù)的增加，土壤原有結(jié)構(gòu)在水分遷移作用下被破壞并且形成新的結(jié)構(gòu)，原始細(xì)小深刻的微觀裂紋和孔隙變得更加寬大淺顯，散亂分布的土壤細(xì)小顆粒在水分遷移的沖刷下集中到一起形成新的土壤顆粒，使得土壤中原有顆粒變得圓潤(rùn)，土壤紋理發(fā)育得成熟。

提取圖7中前4次干濕循環(huán)的微觀裂紋和孔隙像素面積，計(jì)算土微觀裂紋和孔隙部分在土壤掃描部分的占比，其隨干濕循環(huán)次數(shù)的關(guān)系曲線如圖8所示。隨著干濕循環(huán)次數(shù)的增加，土壤試樣微觀裂紋與孔隙的占比逐漸增加，從最初67.12%經(jīng)歷4次干濕循環(huán)后增加到74.83%，在第6次干濕循環(huán)時(shí)土壤微觀裂紋與孔隙占比達(dá)到75.45%，增加速度不斷減小。在5次干濕循環(huán)后，其紋理性基本不再增強(qiáng)，說(shuō)明土體紋理化已經(jīng)完成。

圖8 土壤微觀裂紋和孔隙占比與干濕循環(huán)次數(shù)關(guān)系Fig.8 Relationship between the proportion of soil cracks and pores and the number of dry-wet cycles

4 干濕循環(huán)條件下黃土損傷分析

4.1 干濕循環(huán)條件下黃土損傷模型

隨著干濕循環(huán)次數(shù)的增加，土壤原有結(jié)構(gòu)被破壞，不斷有新的微觀裂紋和孔隙產(chǎn)生，舊有微觀裂紋和孔隙繼續(xù)發(fā)育，土壤結(jié)構(gòu)紋理變得寬大淺顯。干濕循環(huán)對(duì)土壤結(jié)構(gòu)的損傷微細(xì)觀方面主要表現(xiàn)在土壤微觀裂紋和孔隙面積增加，宏觀方面則主要表現(xiàn)在土壤強(qiáng)度的減小。土壤微觀裂紋和孔隙面積占比在一定程度上可以反映干濕循環(huán)對(duì)土壤原有結(jié)構(gòu)的損傷，為此建立黃土干濕循環(huán)損傷模型。對(duì)各灰度圖像指標(biāo)量進(jìn)行擬合，得到各灰度指標(biāo)量隨干濕循環(huán)次數(shù)增加的變化規(guī)律f(n)(n為干濕循環(huán)次數(shù))，對(duì)土壤微觀裂紋和孔隙面積占比隨干濕循環(huán)次數(shù)的變化曲線同樣進(jìn)行擬合，得到土壤微觀裂紋和孔隙面積占比隨干濕循環(huán)次數(shù)增加的變化規(guī)律y(n)。

以干濕循環(huán)次數(shù)為時(shí)間軸，采用LSTM深度學(xué)習(xí)算法建立土壤試樣灰度圖像指標(biāo)量和土壤微觀裂紋和孔隙面積占比之間的時(shí)序回歸預(yù)測(cè)模型，得到土壤微觀裂紋和孔隙面積占比隨灰度圖像指標(biāo)量的變化規(guī)律，如式(7)所示。

y=φ(f1(n)，f2(n)，f3(n)，f4(n)) 。

(7)

式中fi(n)為對(duì)應(yīng)的4個(gè)灰度圖像指標(biāo)變量關(guān)于干濕循環(huán)次數(shù)n的擬合函數(shù)。

土壤微觀裂紋和孔隙面積占比擬合曲線表明土壤在未進(jìn)行干濕循環(huán)時(shí)本身已存在微觀裂紋和孔隙，說(shuō)明土壤存在初始損傷。實(shí)際上，未進(jìn)行干濕循環(huán)的土壤的壓力強(qiáng)度是基于初始損傷條件下得出的。土壤微觀裂紋和孔隙面積占比變化導(dǎo)致了土體應(yīng)力強(qiáng)度的變化。土體有效應(yīng)力與應(yīng)力之比為

(8)

式中：土體應(yīng)力σ=P/S；土體有效應(yīng)力σ′=P/S′；P為土體微元受到的壓力；S為無(wú)損土體截面積；S′為有損土體截面積。

土體干濕循環(huán)有效應(yīng)力與土體未進(jìn)行干濕循環(huán)有效應(yīng)力之比為

(9)

根據(jù)式(10)定義土壤干濕循環(huán)損傷因子D為

(11)

4.2 干濕循環(huán)損傷時(shí)序回歸預(yù)測(cè)

以1組黃土試樣為例，前4次干濕循環(huán)灰度指標(biāo)量擬合曲線如圖9所示。在前2次干濕循環(huán)中，試樣對(duì)比度和熵發(fā)生劇烈變化，說(shuō)明土體中微觀裂紋和孔隙迅速增加。土壤干濕循環(huán)過(guò)程可以分為兩階段：微觀裂紋和孔隙增加、土壤結(jié)構(gòu)紋理化。第1個(gè)階段主要在前2次干濕循環(huán)內(nèi)發(fā)生，第2個(gè)階段主要發(fā)生在前5次干濕循環(huán)內(nèi)。

圖9 灰度圖像指標(biāo)量擬合曲線Fig.9 Fitted curves of gray image index quantities

按0～12次干濕循環(huán)分成1 201個(gè)點(diǎn)，按照擬合曲線可得到灰度指標(biāo)量和土壤微觀裂紋和孔隙面積占比的1 201個(gè)數(shù)據(jù)。同時(shí)建立深度學(xué)習(xí)時(shí)序回歸預(yù)測(cè)模型，包含LSTM層、全連接層和丟棄層。其中LSTM層具有200個(gè)隱含單元，全連接層大小為50，丟棄層丟棄概率為0.5。以0.01的學(xué)習(xí)率，小批量為220進(jìn)行300輪訓(xùn)練，梯度閾值為1，以1組土壤試樣灰度圖像指標(biāo)變量和土壤微觀裂紋與孔隙占比作為訓(xùn)練集合，輸入2組土壤試樣灰度圖像指標(biāo)量，對(duì)土壤微觀裂紋和孔隙占比進(jìn)行時(shí)序回歸預(yù)測(cè)。

第2組土壤試樣損傷因子深度學(xué)習(xí)時(shí)序回歸預(yù)測(cè)曲線與實(shí)際土壤損傷因子變化曲線如圖10所示。

圖10 土壤損傷因子變化曲線Fig.10 Change curve of loess damage factor

深度學(xué)習(xí)時(shí)序回歸預(yù)測(cè)曲線與實(shí)際土壤損傷因子變化曲線吻合度較高，干濕循環(huán)深度學(xué)習(xí)模型具有較好的可靠性。第1組第6次干濕循環(huán)損傷因子為0.249，深度學(xué)習(xí)預(yù)測(cè)曲線第6次損傷因子為0.246，可見(jiàn)，深度學(xué)習(xí)結(jié)果與第1組實(shí)測(cè)數(shù)據(jù)相比也較為貼合。隨著土壤干濕循環(huán)次數(shù)的增加，土壤損傷因子逐漸增大。土壤試樣經(jīng)歷4次干濕循環(huán)后干濕循環(huán)損傷因子損傷比例達(dá)到93.10%。在6次干濕循環(huán)后，土壤損傷因子增加速度幾乎為0，趨近于極限值0.253。在干濕循環(huán)過(guò)程，土壤試樣內(nèi)部水分出現(xiàn)遷移現(xiàn)象，會(huì)產(chǎn)生水分遷移通道，在土壤試樣中形成新的微觀裂紋和孔隙，土壤中的微觀裂紋和孔隙占比不斷增大，在6次干濕循環(huán)后土壤中微觀裂紋和孔隙占比增加到一個(gè)極限值。此時(shí)，在干濕循環(huán)損傷作用下土壤顆粒只進(jìn)行重組，土壤試樣不再產(chǎn)生新的微觀裂紋和孔隙，干濕循環(huán)不再對(duì)土壤試樣產(chǎn)生損傷作用。

5 結(jié) 論

本文對(duì)西安地區(qū)的黃土進(jìn)行干濕循環(huán)處理，基于灰度共生矩陣?yán)命S土試樣電鏡掃描灰度圖像提取土壤紋理特征數(shù)據(jù)，分析土壤微觀結(jié)構(gòu)變化。以灰度圖像最大概率密度分布的灰度為閾值計(jì)算土壤微觀裂紋和孔隙占比，通過(guò)深度學(xué)習(xí)時(shí)序回歸預(yù)測(cè)模型建立土壤紋理特征數(shù)據(jù)和土壤微觀裂紋和孔隙占比之間的聯(lián)系，根據(jù)建立的數(shù)學(xué)模型計(jì)算土壤干濕循環(huán)損傷程度，得到如下結(jié)論：

(1)黃土試樣灰度圖像在放大2 000倍時(shí)能夠反映黃土干濕循環(huán)后微觀結(jié)構(gòu)變化。經(jīng)過(guò)對(duì)比，灰度指標(biāo)量熵和二階矩呈現(xiàn)顯著相關(guān)性，對(duì)比度隨干濕循環(huán)次數(shù)成反“S”型曲線，熵呈“S”型曲線發(fā)展，經(jīng)過(guò)2次干濕循環(huán)后達(dá)到極限值，相關(guān)性和平穩(wěn)度沿凸曲線發(fā)展，在5次干濕循環(huán)后達(dá)到極值。

(2)在干濕循環(huán)后，對(duì)比度減小，相關(guān)性、熵和平穩(wěn)度不斷增大，說(shuō)明土壤試樣微觀裂紋和孔隙不斷增大，紋理性不斷增強(qiáng)。在2次干濕循環(huán)后土壤試樣中基本不再產(chǎn)生新的微觀裂紋和孔隙，土壤顆粒之間只是重復(fù)進(jìn)行重組，形成干濕循環(huán)中與水分遷移方向趨向平行的紋理走向，在5次干濕循環(huán)后紋理性基本不再增強(qiáng)。

(3)建立灰度指標(biāo)變量和土壤試樣微觀裂紋和孔隙占比深度學(xué)習(xí)時(shí)序回歸預(yù)測(cè)模型，可以較好地反映土壤的干濕循環(huán)損傷，經(jīng)6次干濕循環(huán)后，干濕循環(huán)的損傷趨近損傷極限值0.253，說(shuō)明土壤的干濕循環(huán)損傷是有限的，在6次干濕循環(huán)后損傷不再增大。

(4)土壤干濕循環(huán)損傷分為兩個(gè)階段：微觀裂紋和孔隙增加、土壤結(jié)構(gòu)紋理化。第1個(gè)階段主要在前2次干濕循環(huán)內(nèi)發(fā)生，第2個(gè)階段主要發(fā)生在前5次干濕循環(huán)內(nèi)。