農(nóng)田土壤表面干縮裂隙形態(tài)分析與數(shù)值模擬.pdf

1、目前，土壤干縮裂隙的發(fā)展規(guī)律以及形成機理尚不明確，為揭示農(nóng)田土壤表面干縮裂隙發(fā)展規(guī)律以及形成機理，本文對農(nóng)田土壤裂隙進行了試驗研究。應(yīng)用數(shù)字圖像處理技術(shù)對農(nóng)田土壤裂隙試驗圖像進行了二值化、去除噪點、提取裂隙邊緣等處理。閔可夫斯基面積、長度、歐拉數(shù)3個可以描述裂隙空間分布以及連通性的形態(tài)參數(shù)作為裂隙定量分析指標。對裂隙試驗圖像進行形態(tài)學膨脹與腐蝕運算，分析了運算后裂隙試驗圖像的面積、長度、歐拉數(shù)密度的變化規(guī)律。同時構(gòu)建了基于胡克定律的二維

2、土壤開裂模型，模擬土壤收縮動態(tài)開裂的過程。模型假定水分蒸發(fā)是裂隙開裂的主要原因，生成裂隙的土層厚度足夠小。設(shè)置均值為ε，方差為δ2的正態(tài)分布N(ε，δ2）的土壤臨界應(yīng)變隨機場，并引入摩擦力f作為土壤靜態(tài)粘合力，以及表示土壤彈塑性性能的參數(shù)網(wǎng)絡(luò)結(jié)構(gòu)連接強度n(u)，在由結(jié)點組成的六邊形網(wǎng)絡(luò)結(jié)構(gòu)中持續(xù)減小結(jié)點間的距離模擬水分蒸發(fā)，由此引起網(wǎng)絡(luò)結(jié)點受力改變發(fā)生斷裂的過程。本文比較了模擬裂隙圖像與試驗圖像面積、長度、歐拉數(shù)密度。結(jié)果表明:

3、>　　(1)應(yīng)用閔可夫斯基面積、長度、歐拉數(shù)密度函數(shù)可以有效描述裂隙形態(tài);形態(tài)學分析結(jié)果顯示裂隙的閔可夫斯基面積、長度、歐拉數(shù)密度函數(shù)具有不同的變化規(guī)律，其中面積密度函數(shù)表現(xiàn)為0～1之間的增函數(shù);長度密度函數(shù)表現(xiàn)為先增大后減小的變化規(guī)律，歐拉數(shù)密度函數(shù)表現(xiàn)為先增大，后減小再增大的變化規(guī)律。
　　(2)運用基于胡克定律構(gòu)建的二維土壤開裂模型，模擬由于水分蒸發(fā)而引起土壤收縮開裂過程，計算采用蒙特卡洛方法模擬獲得裂隙面積、長度、歐拉數(shù)密

4、度平均值μ以及標準差δ，試驗平均密度值基本分布在模擬密度值μ±δ之間，表明模擬結(jié)果較好。
　　(3)試驗和模擬面積密度、歐拉數(shù)密度平均值相對誤差較小，面積密度相對誤差基本小于10％;歐拉數(shù)密度相對誤差在15％左右;長度密度平均值相對誤差較大，除過38-48 h之間小于10％，其余時段內(nèi)均大于25％，最大相對誤差達到37.38％。應(yīng)用數(shù)值方法模擬土壤表面裂隙，有助于研究農(nóng)田土壤裂隙的形成機理以及土壤裂隙隨時間變化過程中動態(tài)演化規(guī)律。