泰安棕壤土導(dǎo)氣率與飽和導(dǎo)水率關(guān)系初步研究

張 倩，張 瓊，呂 艷，陳官運(yùn)，胡心雨，曾 健，仇 強(qiáng)，張 鵬

(1.山東農(nóng)業(yè)大學(xué)水利土木工程學(xué)院，山東 泰安 271018；2. 山東省果樹(shù)研究所，山東 泰安 271000；3. 山東省海河淮河小清河流域水利管理服務(wù)中心，濟(jì)南 250000；4. 濰坊市農(nóng)業(yè)農(nóng)村局，山東 濰坊 261061)

0 引 言

土壤飽和導(dǎo)水率是土壤重要的物理性質(zhì)之一，它是計(jì)算土壤剖面中水的通量和設(shè)計(jì)灌溉，排水系統(tǒng)工程的一個(gè)重要土壤水力參數(shù)，也是水文模型中的重要參數(shù)[1]。飽和導(dǎo)水率的準(zhǔn)確性影響到模型的準(zhǔn)確性。而土壤飽和導(dǎo)水率的測(cè)定是眾多土壤試驗(yàn)的基礎(chǔ)部分，能為其他試驗(yàn)提供基礎(chǔ)參數(shù)，具有重要的參考價(jià)值。土壤飽和導(dǎo)水率測(cè)定過(guò)程中容易受測(cè)定時(shí)間，測(cè)定方法，取樣尺寸，所選材料等各種因素的影響，導(dǎo)致飽和導(dǎo)水率測(cè)定結(jié)果的有效性受到影響。關(guān)于土壤飽和導(dǎo)水率的測(cè)定，國(guó)內(nèi)外許多學(xué)者進(jìn)行了相關(guān)研究，國(guó)外學(xué)者M(jìn)ertens、Jacques等利用單環(huán)入滲儀測(cè)定了田間飽和導(dǎo)水率[2]；Elrick和Reynolds利用積水入滲法測(cè)定了滲流區(qū)的水力傳導(dǎo)相關(guān)參數(shù)[3]。國(guó)內(nèi)學(xué)者陽(yáng)勇，陳仁生在位于祁連山區(qū)的黑河上游平坦草地開(kāi)展多次野外土壤入滲試驗(yàn)，對(duì)比分析雙環(huán)入滲儀內(nèi)外環(huán)尺寸對(duì)土壤入滲率和飽和導(dǎo)水率的影響[4]，雍晨旭，樊軍分析了坡度對(duì)三種單環(huán)法測(cè)量坡地飽和導(dǎo)水率的影響[5]。當(dāng)前研究中關(guān)于單環(huán)法測(cè)土壤飽和含水率，入滲儀直徑對(duì)測(cè)量結(jié)果的影響研究比較少，本研究在前人的基礎(chǔ)上，通過(guò)三種不同直徑的單環(huán)入滲儀探究分析入滲儀直徑與飽和導(dǎo)水率之間的關(guān)系，以期為土壤飽和導(dǎo)水率的測(cè)量提供參考。

土壤導(dǎo)氣率反映了氣體通過(guò)土壤的能力，是土壤空氣更新快慢的綜合指標(biāo)。土壤導(dǎo)氣率通常是通過(guò)向土壤的一端施加一定的氣體壓力，然后測(cè)量通過(guò)土壤的恒定氣體流量而得出的，土壤導(dǎo)氣能力還可作為分析土壤孔隙幾何形狀、結(jié)構(gòu)以及土壤穩(wěn)定性等的指標(biāo)[6]。因此，對(duì)土壤導(dǎo)氣特性的研究受到很高的重視。國(guó)外學(xué)者M(jìn)assmann等在空氣不可壓縮的情況下室內(nèi)測(cè)定土壤導(dǎo)氣率，試驗(yàn)方法僅適用于測(cè)定氣壓變化范圍較小時(shí)的導(dǎo)氣率[7]。Evans和Kirkham在田間實(shí)地應(yīng)用土壤導(dǎo)氣率的瞬態(tài)測(cè)量方法進(jìn)行導(dǎo)氣率的測(cè)量[8]。Jalbert和Dane提出了侵入式和非侵入式測(cè)量田間土壤導(dǎo)氣率的裝置，并推導(dǎo)出了穩(wěn)態(tài)土壤導(dǎo)氣率測(cè)算模型[9]。國(guó)內(nèi)學(xué)者王德勝，張振華對(duì)河南省五個(gè)地區(qū)代表性土壤的導(dǎo)氣率進(jìn)行了研究分析[10]。

在理想狀況下，完全干燥的土壤樣品的導(dǎo)氣率和飽和導(dǎo)水率具有相同的值。然而，在實(shí)際中飽和導(dǎo)水率和導(dǎo)氣率沒(méi)有可比性[11]。國(guó)內(nèi)外一些學(xué)者研究發(fā)現(xiàn)并預(yù)測(cè)Ka與Ks的關(guān)系[11-15]，Loll等人建議在-10 kPa的壓強(qiáng)下測(cè)量原狀土樣本獲得導(dǎo)氣率并推斷飽和導(dǎo)水率的空間變異性[12]。王衛(wèi)華、王全九等使用土壤導(dǎo)氣率測(cè)量?jī)x對(duì)長(zhǎng)武小麥試驗(yàn)田不同含水率、取樣方向、土層、根系密度下的導(dǎo)氣率進(jìn)行了研究，同時(shí)對(duì)土壤飽和導(dǎo)水率與導(dǎo)氣率的關(guān)系進(jìn)行了研究[11]。雖然國(guó)內(nèi)外學(xué)者在導(dǎo)氣率和飽和導(dǎo)水率關(guān)系研究中做了一些工作，但仍需要更多的試驗(yàn)研究來(lái)驗(yàn)證這種方法的可行性。

本研究在前人研究的基礎(chǔ)上應(yīng)用單環(huán)入滲法和穩(wěn)態(tài)三維土壤導(dǎo)氣率測(cè)量法測(cè)定泰安棕壤土飽和導(dǎo)水率與導(dǎo)氣率，研究不同直徑入滲環(huán)對(duì)飽和導(dǎo)水率的影響并將實(shí)驗(yàn)中得到的導(dǎo)氣率與飽和導(dǎo)水率建立關(guān)系，進(jìn)一步明確土壤飽和導(dǎo)水率與導(dǎo)氣率之間存在的函數(shù)關(guān)系，以期為土壤飽和導(dǎo)水率的測(cè)定以及通過(guò)土壤導(dǎo)氣率預(yù)測(cè)土壤飽和導(dǎo)水率的研究提供參考。

1 試驗(yàn)方法

1.1 試驗(yàn)地概況

試驗(yàn)所選場(chǎng)地為泰安地區(qū)山東省農(nóng)業(yè)科學(xué)院核桃國(guó)家種植資源圃。泰安屬于溫帶大陸性半濕潤(rùn)季風(fēng)氣候區(qū)，平均氣溫12.9 ℃，年平均日照時(shí)間2 627 h，無(wú)霜期平均195 d，多年平均降雨697 mm。受季風(fēng)氣候的影響，年際降雨變幅比較大。

試驗(yàn)場(chǎng)地的土壤質(zhì)地均勻疏松，通透性好，土壤類型為棕壤土。在現(xiàn)場(chǎng)除去表層浮土后，取深度為0～10 cm的土壤。從試驗(yàn)地取樣回實(shí)驗(yàn)室后，將土放在通風(fēng)處，自然風(fēng)干，過(guò)2 mm孔徑的土篩后利用激光粒度儀法對(duì)土壤進(jìn)行顆粒分析。如表1所示。

表1 土壤的顆粒組成Tab.1 The granular composition of soil

土壤的容重采用烘干法測(cè)定。用100 cm3的環(huán)刀切割未攪動(dòng)的自然狀態(tài)土樣，使土樣充滿其中，烘干后稱量并計(jì)算單位容積的土重量。土壤容重的計(jì)算公式：

(1)

式中：γs為土壤容重，g/cm3；M烘干為烘干土的質(zhì)量，g；V為環(huán)刀的容積，cm3。

土壤孔隙度的計(jì)算公式：

(2)

式中：Ф為孔隙度，%；γs為土壤容重，g/cm3；γd為土粒密度，g/cm3；試驗(yàn)中γd取2.65 g/cm3。

試驗(yàn)場(chǎng)地任意選取5個(gè)測(cè)點(diǎn)，測(cè)點(diǎn)編號(hào)分別命名為G1、G2、G3、G4、G5，每個(gè)測(cè)點(diǎn)重復(fù)取3個(gè)100 cm3環(huán)刀土樣、3個(gè)100 cm3鋁盒土樣。各測(cè)點(diǎn)測(cè)得的土樣參數(shù)取平均數(shù)。

(3)

試驗(yàn)場(chǎng)地各測(cè)點(diǎn)的土壤容重，孔隙度，初始含水率，飽和含水率如表2所示。

表2 不同直徑入滲環(huán)所測(cè)值Tab.2 Values measured by infiltration rings with different diameters

1.2 試驗(yàn)材料

本研究采用較為常用單環(huán)入滲法測(cè)定土壤飽和導(dǎo)水率。試驗(yàn)裝置(圖1)包括量筒、鐵架臺(tái)、輸水管、閥門(mén)、直徑10 cm的入滲環(huán)等。量筒放置于水平鐵架臺(tái)上作為供水裝置；輸水管用于向入滲環(huán)中按要求輸水；閥門(mén)用于控制水流流量大小以及輸水管的開(kāi)關(guān)。

圖1 單環(huán)入滲土壤飽和導(dǎo)水率測(cè)量裝置Fig.1 The measurement device Single ring infiltration soil Saturated water conductivity

土壤導(dǎo)氣率的測(cè)量采用穩(wěn)態(tài)三維土壤導(dǎo)氣率測(cè)量裝置。試驗(yàn)裝置主要有：空氣壓縮機(jī)、閥門(mén)、氣體流量計(jì)、輸氣管、土柱管和U形管壓力計(jì)等，土柱管的直徑D為10 cm。試驗(yàn)裝置如圖2所示。

1-空氣壓縮機(jī)；2-充電開(kāi)關(guān)；3-壓力指示表；4-閥門(mén)；5-減壓閥；6-氣體流量計(jì)；7-導(dǎo)氣軟管；8-導(dǎo)氣硬管；9-土柱管；10-橡皮塞；11-U形管壓力計(jì)；12-地面圖2 穩(wěn)態(tài)三維土壤導(dǎo)氣率測(cè)量裝置Fig.2 Steady-state three-dimensional soil air permeability measuring device

1.3 數(shù)學(xué)模型

1.3.1 土壤飽和導(dǎo)水率模型

單環(huán)入滲試驗(yàn)是測(cè)量土壤水分入滲速率及土壤飽和導(dǎo)水率較為常用的方法?；谕寥罈l件相同情況下單環(huán)入滲的入滲速率是一維入滲時(shí)的f倍(f是與試驗(yàn)土壤初始條件、周邊條件及入滲環(huán)形狀有關(guān)的參數(shù))的假定以及入滲環(huán)中水分在土壤中屬于三維運(yùn)動(dòng)過(guò)程，Wu和Pan[16]在前人研究基礎(chǔ)上得出了單環(huán)入滲過(guò)程的概化解，概化解與測(cè)定入滲曲線結(jié)合，得累積入滲量的求解公式[17]：

I=afKst+2bfKs(tTc)0.5

(4)

I=At+Bt0.5

(5)

A和B可通過(guò)實(shí)測(cè)的累積入滲量與觀測(cè)時(shí)間的擬合關(guān)系求得。土壤飽和導(dǎo)水率的計(jì)算公式：

(6)

(7)

(8)

(9)

G*=d+r/2

(10)

Δθ=θ0-θi

(11)

式中：a，b均為常數(shù)，a= 0.908 4，b= 0.168 2；t為入滲時(shí)間，min；H為入滲環(huán)中水的深度，cm；d為入滲環(huán)插入土中的深度，cm；r為入滲環(huán)半徑，cm；Tc為時(shí)間尺度特征值；Ks為飽和導(dǎo)水率，cm/min；θ0和θi分別為土壤飽和含水率和初始含水率。

1.3.2 土壤導(dǎo)氣率模型

Jalbert等推導(dǎo)出的三維穩(wěn)態(tài)土壤導(dǎo)氣率關(guān)系式[9]：

(12)

(13)

式中：μ為空氣動(dòng)力黏滯系數(shù)，Pa·s；D為土柱管直徑，m；Q為氣體流量，m3/s；ΔP為土柱管內(nèi)密閉空間氣體的壓強(qiáng)值，Pa；H為土柱管插入土壤中的深度，m；G為土柱管的形狀系數(shù)；Ka為穩(wěn)態(tài)三維土壤導(dǎo)氣率，m2。

1.3.3 土壤導(dǎo)氣率與飽和導(dǎo)水率關(guān)系模型

經(jīng)過(guò)Loll等人的探索，發(fā)現(xiàn)了土壤導(dǎo)氣率和飽和導(dǎo)水率之間對(duì)數(shù)關(guān)系[11,12]：

Lg(Ks)=αLg(Ka)+β

(14)

式中：α和β為待求參數(shù)。

1.4 試驗(yàn)設(shè)計(jì)

1.4.1 不同入滲環(huán)直徑對(duì)飽和導(dǎo)水率Ks的影響試驗(yàn)

選取試驗(yàn)場(chǎng)內(nèi)五處測(cè)點(diǎn)G1、G2、G3、G4、G5，分別采用10、20、30 cm 3種直徑入滲環(huán)對(duì)土壤飽和導(dǎo)水率進(jìn)行測(cè)量，獲取入滲過(guò)程中各時(shí)刻所對(duì)應(yīng)的入滲速率和累積入滲量。試驗(yàn)過(guò)程中，清除土壤表面枝葉等雜物，分別將3種不同直徑的入滲環(huán)垂直打入土中，打入深度為10 cm，打入過(guò)程保證入滲環(huán)不發(fā)生傾斜，以保證土壤的原狀結(jié)構(gòu)，并使入滲環(huán)周圍與土壤緊貼，以防水流沿邊壁流走。入滲環(huán)安裝完成后，將環(huán)內(nèi)部橫鋪一層低阻力尼龍網(wǎng)，以防止水流對(duì)土體表面的沖刷，破壞土體表面的孔隙。測(cè)量直徑10 cm的入滲環(huán)，按照每入滲50 mL水記錄一次時(shí)間，測(cè)量直徑20和30 cm的入滲環(huán)，按照每入滲100 mL水記錄一次時(shí)間，每次記錄數(shù)據(jù)時(shí)保證入滲環(huán)內(nèi)未見(jiàn)明顯積水，每個(gè)測(cè)點(diǎn)測(cè)量時(shí)間為120 min。根據(jù)所記錄的時(shí)間和入滲水量，求出直徑為10、20、30 cm的入滲環(huán)所測(cè)的累積入滲量和入滲速率，并繪制入滲速率-時(shí)間和累計(jì)入滲量-時(shí)間關(guān)系曲線，將所得數(shù)據(jù)代入土壤飽和導(dǎo)水率數(shù)學(xué)模型，最終求得土壤飽和導(dǎo)水率。

1.4.2 三維穩(wěn)態(tài)土壤導(dǎo)氣率試驗(yàn)

選取試驗(yàn)場(chǎng)內(nèi)飽和導(dǎo)水率測(cè)量點(diǎn)G1、G2、G3、G4、G5附近的點(diǎn)，作為測(cè)量土壤導(dǎo)氣率與土壤飽和導(dǎo)水率關(guān)系的試驗(yàn)點(diǎn)。將直徑為10 cm的土柱管垂直打入土中，打入深度為10 cm，為保證原狀土不受擾動(dòng)，打入過(guò)程要保證土柱管不發(fā)生傾斜，并使土柱管周圍與土壤緊貼。連接好儀器設(shè)備，接通電源，打開(kāi)空氣壓縮機(jī)閥門(mén)，調(diào)節(jié)減壓閥使氣壓控制在一定合適的范圍，調(diào)節(jié)氣體流量計(jì)，使壓縮機(jī)輸出的氣體在一定流量下通過(guò)測(cè)量裝置(試驗(yàn)采用3、4、5、6、7、8、9、10 L/min氣體通過(guò)流量)，當(dāng)氣體傳導(dǎo)速率達(dá)到穩(wěn)定時(shí)，通過(guò)U形測(cè)壓計(jì)兩側(cè)水柱高度差讀取土柱管內(nèi)土體密封端壓強(qiáng)值，記錄氣體流量計(jì)讀數(shù)Q和U形測(cè)壓計(jì)兩側(cè)水柱高度差Δh。將高度差Δh轉(zhuǎn)化為封閉氣體壓強(qiáng)ΔP，即可得到通過(guò)的氣體流量Q與儀器內(nèi)封閉氣體壓強(qiáng)ΔP之間的關(guān)系。根據(jù)三維穩(wěn)態(tài)導(dǎo)氣率數(shù)學(xué)模型求出各氣體流量下土壤導(dǎo)氣率Ka，對(duì)各流量測(cè)得的導(dǎo)氣率結(jié)果取平均值作為本測(cè)點(diǎn)的土壤導(dǎo)氣率。

2 結(jié)果與分析

2.1 不同入滲環(huán)直徑對(duì)飽和導(dǎo)水率Ks的影響

2.1.1 不同入滲環(huán)直徑所測(cè)入滲速率與累計(jì)入滲量的研究結(jié)果

根據(jù)所記錄的入滲時(shí)間和水量，求出直徑為10、20、30 cm的入滲環(huán)所測(cè)的累積入滲量和入滲速率，并繪制入滲速率-時(shí)間和累計(jì)入滲量-時(shí)間曲線。圖3分別為測(cè)點(diǎn)G1、G2、G3、G4、G5對(duì)應(yīng)的不同直徑入滲環(huán)所測(cè)的土壤水分入滲速率變化曲線。圖4分別為測(cè)點(diǎn)G1、G2、G3、G4、G5對(duì)應(yīng)的不同直徑入滲環(huán)所測(cè)累計(jì)入滲量變化曲線。由圖3可知，初始入滲階段前10 min左右，由于表層土壤初始含水率不高，所以水分入滲較為迅速，隨著入滲時(shí)間的增加，入滲速率急劇下降，80 min后，土壤入滲速率逐漸趨于穩(wěn)定。從5個(gè)測(cè)點(diǎn)的入滲速率的整個(gè)變化過(guò)程來(lái)看，直徑為10 cm入滲環(huán)所測(cè)入滲速率最大，20 cm入滲環(huán)次之，30 cm入滲環(huán)最小。由圖4可知，不同直徑入滲環(huán)所測(cè)的累計(jì)入滲量隨著時(shí)間的增加不斷增大，其中直徑10 cm入滲環(huán)所測(cè)累計(jì)入滲量最大，20 cm入滲環(huán)次之，30 cm入滲環(huán)最小。隨著入滲時(shí)間的增加，累計(jì)入滲量的增加值基本趨于穩(wěn)定，這與入滲速率的變化規(guī)律吻合。

圖3 G1、G2、G3、G4、G5入滲速率Fig.3 Infiltration rates of G1, G2, G3, G4, G5

圖4 G1、G2、G3、G4、G5累計(jì)入滲量Fig.4 Cumulative infiltration amount of G1, G2, G3, G4, G5

2.1.2 不同入滲環(huán)直徑土壤飽和導(dǎo)水率Ks的求解

根據(jù)上述所測(cè)的累計(jì)入滲量和入滲時(shí)間，擬合公式(5)的參數(shù)A、B(各測(cè)點(diǎn)所得擬合曲線的R2均接近1)，并將所測(cè)數(shù)據(jù)代入式(6)，得出土壤飽和導(dǎo)水率。不同直徑入滲環(huán)所測(cè)Ks值見(jiàn)表3。

通過(guò)分析表3的數(shù)據(jù)可以發(fā)現(xiàn)，利用直徑為10、20、30 cm的入滲環(huán)獲得的數(shù)據(jù)測(cè)得的飽和導(dǎo)水率具有明顯的差異性，且不具有規(guī)律性，入滲環(huán)直徑的取值影響單環(huán)入滲法測(cè)量飽和導(dǎo)水率的準(zhǔn)確性。為了避免入滲環(huán)直徑選取不同在求解飽和導(dǎo)水率時(shí)產(chǎn)生的影響，本研究在前人研究的基礎(chǔ)上，通過(guò)建立泰安棕壤土導(dǎo)氣率與飽和導(dǎo)水率的關(guān)系，預(yù)測(cè)土壤飽和導(dǎo)水率。

2.2 土壤導(dǎo)氣率與飽和導(dǎo)水率的關(guān)系

根據(jù)氣體流量數(shù)值和氣體壓強(qiáng)值繪制出二者關(guān)系圖，如圖5所示，結(jié)果表明隨著氣體流量的增加，對(duì)應(yīng)的封閉氣體壓強(qiáng)值ΔP也相應(yīng)變大；相同氣體流量下G5點(diǎn)的氣體壓強(qiáng)最大，G2點(diǎn)最小。對(duì)氣體流量與氣體壓強(qiáng)值進(jìn)行線性擬合如表4所示，由表4可知?dú)怏w流量Q與儀器內(nèi)封閉氣體壓強(qiáng)ΔP二者存在明顯的線性關(guān)系。

對(duì)5個(gè)測(cè)點(diǎn)G1、G2、G3、G4、G5三維穩(wěn)態(tài)法測(cè)量的土壤導(dǎo)氣率指標(biāo)統(tǒng)計(jì)分析見(jiàn)表5。由表5可知，5個(gè)測(cè)點(diǎn)測(cè)量值在1.35×10-11與7.03×10-11m2之間，變異系數(shù)Cv在0.07與0.13之間，可見(jiàn)各測(cè)點(diǎn)數(shù)據(jù)變異程度較小，數(shù)據(jù)波動(dòng)較小。

擬合實(shí)驗(yàn)測(cè)得導(dǎo)氣率與飽和導(dǎo)水率的對(duì)數(shù)關(guān)系如圖6所示，可知擬合線性方程系數(shù)α和β分別為1.737 9、-4.230 1，擬合Ka與Ks取對(duì)數(shù)后關(guān)系直線的R2達(dá)到0.874 1，說(shuō)明田間試驗(yàn)測(cè)得的土壤導(dǎo)氣率與飽和導(dǎo)水率關(guān)系密切，導(dǎo)氣率與飽和導(dǎo)水率的對(duì)數(shù)關(guān)系曲線可用于泰安地區(qū)棕壤土土壤飽和導(dǎo)水率預(yù)測(cè)。

圖5 氣體流量Q與封閉氣體壓強(qiáng)值ΔP關(guān)系圖Fig.5 Relation diagram between gas flow rate Q and closed gas pressure value ΔP

表4 線性擬合的相關(guān)性分析Tab.4 Correlation analysis of linear fitting

3 結(jié) 論

本文利用快速而簡(jiǎn)單的土壤飽和導(dǎo)水率單環(huán)入滲法，以泰安地區(qū)棕壤土飽和導(dǎo)水率為研究對(duì)象，探究10、20和30 cm 3個(gè)不同入滲環(huán)直徑對(duì)土壤飽和導(dǎo)水率測(cè)量值的影響，結(jié)果表明：

表5 三維穩(wěn)態(tài)法導(dǎo)氣率特征值Tab.5 Characteristic values of air permeability in three-dimensional steady state method

圖6 泰安地區(qū)導(dǎo)氣率與導(dǎo)水率對(duì)數(shù)關(guān)系圖Fig.6 Log-log relationship between saturated water conductivity and air permeability in taian county

(1)整個(gè)入滲過(guò)程入滲速率、累計(jì)入滲量具有相同的變化趨勢(shì)，數(shù)值最大的是直徑10 cm入滲環(huán)，20 cm入滲環(huán)次之，最小的是30 cm入滲環(huán)。利用不同直徑的入滲環(huán)所測(cè)飽和導(dǎo)水率值具有顯著差異性，但不具有規(guī)律性，入滲環(huán)直徑的取值影響單環(huán)入滲法測(cè)量飽和導(dǎo)水率的準(zhǔn)確性。

(2)利用三維穩(wěn)態(tài)裝置測(cè)量5個(gè)測(cè)點(diǎn)的導(dǎo)氣率，結(jié)果表明隨著氣體流量的增加，對(duì)應(yīng)的封閉氣體壓強(qiáng)值ΔP也相應(yīng)變大；氣體流量Q與封閉氣體壓強(qiáng)ΔP之間具有明顯的線性關(guān)系。

(3)擬合導(dǎo)氣率與飽和導(dǎo)水率的對(duì)數(shù)關(guān)系曲線R2=0.874 1，導(dǎo)氣率與飽和導(dǎo)水率之間有較明顯的線性關(guān)系，導(dǎo)氣率與飽和導(dǎo)水率的對(duì)數(shù)關(guān)系曲線可用于泰安地區(qū)棕壤土飽和導(dǎo)水率預(yù)測(cè)。