馬 虎, 謝宗齊, 鄧 利, 薛賽男, 周長(zhǎng)省

(1. 南京理工大學(xué) 機(jī)械工程學(xué)院, 南京 210094; 2. 北京空天技術(shù)研究所, 北京 100074)

0 引 言

旋轉(zhuǎn)爆震發(fā)動(dòng)機(jī)(Rotating Detonation Engine，簡(jiǎn)稱RDE)是一種采用一道或多道爆震波在燃燒室內(nèi)連續(xù)周向傳播進(jìn)而產(chǎn)生推力的新型推力動(dòng)力裝置，它具有熱效率高、結(jié)構(gòu)緊湊以及連續(xù)工作的優(yōu)點(diǎn)[1]，近年來受到廣泛關(guān)注。

旋轉(zhuǎn)爆震發(fā)動(dòng)機(jī)研究在過去幾年取得了顯著的進(jìn)展，主要包括旋轉(zhuǎn)爆震波穩(wěn)定傳播時(shí)的質(zhì)量流率以及當(dāng)量比范圍[2- 5]，發(fā)動(dòng)機(jī)的單雙波工作模態(tài)及其影響因素[4- 11]，燃料與氧化劑噴注結(jié)構(gòu)對(duì)發(fā)動(dòng)機(jī)工作過程的影響[2,7- 9,11]，RDE的推力性能[12- 15]以及RDE與其他發(fā)動(dòng)機(jī)的組合應(yīng)用[16- 17]。相關(guān)實(shí)驗(yàn)結(jié)果表明，對(duì)于給定的燃燒室構(gòu)型，發(fā)動(dòng)機(jī)工作模態(tài)與質(zhì)量流率、反應(yīng)物活性以及出口阻塞比有關(guān)。出口阻塞比為零時(shí)，增加質(zhì)量流率使得發(fā)動(dòng)機(jī)從單波逐漸往雙波傳播

過渡[7- 8]；增加反應(yīng)物活性，減小了旋轉(zhuǎn)爆震波傳播所需的臨界高度，有利于切向多波的形成[6]；增加出口阻塞比，提高了燃燒室壓力，使得爆震波的特征化學(xué)反應(yīng)尺寸減小，也為切向多波的形成提供了條件[9]。

旋轉(zhuǎn)爆震發(fā)動(dòng)機(jī)的推力性能研究表明，為了實(shí)現(xiàn)RDE熱效率高以及自增壓的特性，需采用雍塞的噴管將高溫高壓爆震產(chǎn)物等熵膨脹到環(huán)境狀態(tài)[12]；并且，在尋求RDE與渦輪的組合應(yīng)用時(shí)，需采取措施抑制RDE出口流動(dòng)的不均勻性以及添加渦輪導(dǎo)向器。然而，在燃燒室出口集成雍塞型噴管或渦輪導(dǎo)向器時(shí)，燃燒室出口的流動(dòng)狀態(tài)將發(fā)生改變，進(jìn)而導(dǎo)致發(fā)動(dòng)機(jī)工作模態(tài)發(fā)生變化，在一定條件下燃燒室內(nèi)將產(chǎn)生軸向脈沖爆震現(xiàn)象[10]，即爆震波沿燃燒室軸向傳播。此外，采用氧化劑亞聲速噴注時(shí)，當(dāng)燃燒室壓力與反壓比值較小且爆震產(chǎn)物在燃燒室出口發(fā)生過膨脹時(shí)也將出現(xiàn)軸向脈沖爆震現(xiàn)象[18]。Bykovskii[18]、Anand[11]以及王超[19]對(duì)環(huán)形燃燒室內(nèi)的軸向脈沖爆震現(xiàn)象進(jìn)行了較為詳細(xì)的研究，認(rèn)為軸向脈沖爆震主要發(fā)生在氧化劑亞聲速噴注的條件下，且對(duì)應(yīng)的當(dāng)量比較低。然而，不同結(jié)構(gòu)的燃燒室產(chǎn)生軸向脈沖爆震的方式存在一定的差異，Anand[11]在燃燒室出口安裝收斂噴管，軸向傳播的激波在出口處反射，再通過SWACER機(jī)制發(fā)展成為爆震波；而Bykovskii[18]以及王超[19]的實(shí)驗(yàn)結(jié)果表明，由于燃燒室內(nèi)的爆震產(chǎn)物在出口處出現(xiàn)了過膨脹，進(jìn)而產(chǎn)生了激波，波后未完全燃燒的反應(yīng)物進(jìn)一步燃燒釋放能量，增加激波的強(qiáng)度并推動(dòng)其往燃燒室上游傳播。

軸向脈沖爆震作為旋轉(zhuǎn)爆震燃燒室內(nèi)的一種特殊工作模態(tài)，目前人們對(duì)其認(rèn)識(shí)仍然有限，其產(chǎn)生條件和機(jī)制尚不明確，需進(jìn)一步開展研究，以促進(jìn)對(duì)RDE燃燒室內(nèi)工作模態(tài)形成和轉(zhuǎn)變的理解，推動(dòng)RDE走向工程應(yīng)用。因此，本文在已有研究的基礎(chǔ)上，進(jìn)一步探索發(fā)動(dòng)機(jī)以軸向脈沖爆震模態(tài)工作時(shí)的特點(diǎn)，結(jié)合瞬態(tài)壓力與離子信號(hào)以及高速攝影分析軸向脈沖爆震的傳播過程，并采用線性聲學(xué)理論對(duì)該模態(tài)下的工作頻率進(jìn)行預(yù)測(cè)。

1 實(shí)驗(yàn)設(shè)備與方法介紹

實(shí)驗(yàn)系統(tǒng)示意圖如圖1所示，主要由模型發(fā)動(dòng)機(jī)、燃料和氧化劑供給系統(tǒng)、時(shí)序控制系統(tǒng)以及高頻數(shù)據(jù)采集系統(tǒng)構(gòu)成。RDE模型發(fā)動(dòng)機(jī)內(nèi)外徑分別為70和80mm，長(zhǎng)度為40mm。氧化劑與燃料分別采用環(huán)縫與小孔對(duì)撞的模式進(jìn)行混合，其中燃料H2通過均勻分布在發(fā)動(dòng)機(jī)內(nèi)環(huán)面的90個(gè)小孔噴射入燃燒室內(nèi)，氧化劑空氣則通過Laval類型的環(huán)縫進(jìn)入燃燒室與燃料進(jìn)行快速混合，如圖2所示。燃料和氧化劑的流量通過安裝在管路中的聲速噴嘴進(jìn)行計(jì)算，實(shí)驗(yàn)中的全局工作當(dāng)量比采用管路中的燃料與氧化劑質(zhì)量流率計(jì)算。發(fā)動(dòng)機(jī)利用切向安裝的預(yù)爆震管進(jìn)行起爆，預(yù)爆震管內(nèi)徑為6mm，長(zhǎng)度600mm，其內(nèi)填充滿H2和O2的混合氣體，點(diǎn)火時(shí)刻選擇為燃料和氧化劑積氣腔壓力處于穩(wěn)定的區(qū)間。為了減少燃燒室填充階段對(duì)預(yù)爆震管內(nèi)混合氣體的影響，預(yù)爆震管出口與燃燒室連接處采用薄膜隔離。此外，燃燒室出口阻塞比通過安裝收斂噴管獲得。

圖1 實(shí)驗(yàn)系統(tǒng)示意圖

圖2 噴注結(jié)構(gòu)示意圖

本文主要采用壓力、離子信號(hào)以及高速攝影的方式來對(duì)發(fā)動(dòng)機(jī)是否成功起爆及其工作模態(tài)進(jìn)行判斷。高頻壓力傳感器與離子探針位置如圖3所示，圖中PCB1、PCB2為高頻壓電式壓力傳感器，I1為離子探針。高頻壓電傳感器型號(hào)為PCB113B24，其諧振頻率≥500kHz，上升時(shí)間≤1μs。燃料與氧化劑積氣腔壓力分別采用擴(kuò)散硅式傳感器進(jìn)行測(cè)量，測(cè)量精度≤0.5%FS。高頻數(shù)據(jù)采集采用NI高頻數(shù)采系統(tǒng)，其中基于USB的NI X系列多功能DAQ憑借NI- STC3定時(shí)和同步技術(shù)具有8路同步模擬信號(hào)輸入，單通道的采樣頻率最高可達(dá)2MHz，整個(gè)工作過程中每通道的采樣頻率設(shè)置為500kHz，足以捕捉到燃燒室內(nèi)爆震波傳播的信息。為了得到發(fā)動(dòng)機(jī)工作模態(tài)的轉(zhuǎn)變過程，采用高速相機(jī)(Phantom V210)從發(fā)動(dòng)機(jī)尾部對(duì)燃燒室內(nèi)傳播的火焰進(jìn)行拍攝，相機(jī)拍攝頻率為47 000幀/s，每幅照片分辨率為208pixel×120pixel。由于平齊安裝的PCB傳感器受到的熱載荷嚴(yán)重，因此，為了保護(hù)傳感器，實(shí)驗(yàn)工作時(shí)間設(shè)置為0.2s。

圖3 傳感器與離子探針安裝位置

2 實(shí)驗(yàn)結(jié)果與討論

2.1 軸向脈沖爆震模態(tài)工作過程

在一定的出口阻塞比條件下，增加反應(yīng)物質(zhì)量流率，燃燒室內(nèi)將出現(xiàn)自持傳播的軸向脈沖爆震。圖4和5分別為燃燒室出口阻塞比BR為0.6和0.7、空氣質(zhì)量流率為196.7g/s、當(dāng)量比為1.63時(shí)，發(fā)動(dòng)機(jī)以軸向脈沖爆震模態(tài)工作的壓力與離子信號(hào)。

圖4(a)和5(a)為工作過程中的壓力與離子信號(hào)，可以發(fā)現(xiàn)在質(zhì)量流率以及當(dāng)量比相同時(shí)，出口阻塞比為0.6和0.7的工作過程基本相同。在初始階段，PCB測(cè)量得到的動(dòng)態(tài)壓力幅值較小，相應(yīng)的離子信號(hào)也比較微弱；經(jīng)過一段時(shí)間之后，發(fā)動(dòng)機(jī)以爆震模態(tài)工作，測(cè)量的動(dòng)態(tài)壓力幅值也隨之增加，相應(yīng)的離子信號(hào)大大增強(qiáng)。此外，積氣腔的平均壓力表明，在未建立爆震工作模態(tài)之前，空氣積氣腔受到燃燒室的影響較弱，而在建立爆震工作模態(tài)之后，空氣積氣腔平均壓力明顯增加。圖4(b)和5(b)為對(duì)瞬態(tài)壓力信號(hào)進(jìn)行STFT處理的結(jié)果，可以明顯看出，發(fā)動(dòng)機(jī)在整個(gè)工作過程中出現(xiàn)了2種不同的工作模態(tài)，在工作初始階段，發(fā)動(dòng)機(jī)工作主頻處于5000～5500Hz之間，而以爆震模態(tài)工作時(shí)，發(fā)動(dòng)機(jī)的工作主頻接近13 000Hz。

為進(jìn)一步了解發(fā)動(dòng)機(jī)在不同階段內(nèi)的工作細(xì)節(jié)，將其信號(hào)放大。其中，圖4(c)與5(c)分別表示阻塞比為0.6和0.7時(shí)，發(fā)動(dòng)機(jī)以非爆震模態(tài)工作時(shí)的壓力與離子信號(hào)。在此階段內(nèi)，PCB測(cè)量的動(dòng)態(tài)壓力幅值較低且呈正弦波形，離子信號(hào)在其幅值附近有較小的波動(dòng)；并且PCB1與PCB2的動(dòng)態(tài)壓力信號(hào)存在一定的相位差，表示該波動(dòng)是沿燃燒室切向傳播的。綜合以上特征，推斷此時(shí)燃燒室內(nèi)可能發(fā)生了一階切向的聲學(xué)耦合燃燒，在該工況條件下，采用等壓燃燒假設(shè)得到的燃燒產(chǎn)物聲速約為1056m/s，計(jì)算得到聲學(xué)耦合燃燒的一階振頻為4500Hz，與實(shí)際工作主頻存在一定的差距，這可能是由于實(shí)際工作過程中燃料與氧化劑的混合不均勻，導(dǎo)致局部區(qū)域的當(dāng)量比較大，

(a) 工作過程總體信號(hào)

(b) 短時(shí)傅里葉變換

(c) 切向聲學(xué)耦合燃燒

(d) 軸向脈沖爆震

圖4 質(zhì)量流率196.7g/s，出口阻塞比為0.6，當(dāng)量比為1.63的信號(hào)特征

Fig.4Signalcharacteristicswithairmassflowrate196.7g/s,blockageratio0.6,equivalenceratio1.63

(a) 工作過程總體信號(hào)

(b) 短時(shí)傅里葉變換

(c) 切向聲學(xué)耦合燃燒

(d) 軸向脈沖爆震

圖5 質(zhì)量流率196.7g/s，出口阻塞比為0.7，當(dāng)量比為1.63的信號(hào)特征

Fig.5Signalcharacteristicswithairmassflowrate196.7g/s,blockageratio0.7,equivalenceratio1.63

從而使得燃燒室內(nèi)的聲速在不同的區(qū)域存在差別。

圖4(d)和5(d)為發(fā)動(dòng)機(jī)以爆震模態(tài)工作時(shí)的瞬態(tài)壓力與離子信號(hào)放大圖。在該工作階段內(nèi)，PCB1與PCB2的信號(hào)在同一時(shí)刻幾乎是重合的?？紤]PCB1與PCB2的位置分布，推斷此時(shí)燃燒室可能存在均勻分布的同向三波或軸向傳播的爆震波2種情況。為進(jìn)一步分辨發(fā)動(dòng)機(jī)的工作模態(tài)，結(jié)合該工況對(duì)應(yīng)的尾部高速攝影對(duì)爆震波的傳播方向進(jìn)行分析。圖6(a)～(h)為圖5對(duì)應(yīng)工況下的連續(xù)8張尾部高速攝影。首先，圖6顯示的熒光圖亮度變化表明，在每個(gè)時(shí)刻燃燒室周向位置都充滿火焰，這與燃燒室內(nèi)爆震波切向傳播的尾部高速攝影存在明顯區(qū)別[20]；其次，圖6中的熒光亮度顯示燃燒室內(nèi)的燃燒存在明顯的周期行為，其中圖6(a)～(d)為一個(gè)周期，而(e)～(h)為下一個(gè)周期。由于高速攝影拍攝頻率為47 000幀/s，由此計(jì)算的工作頻率約為11 750Hz，與壓力信號(hào)的主頻相近。因此，綜合壓力信號(hào)以及發(fā)動(dòng)機(jī)尾部高速攝影，判定燃燒室內(nèi)存在軸向傳播的爆震波。

圖6 質(zhì)量流率196.7g/s，出口阻塞比為0.7，當(dāng)量比為1.63的尾部高速攝影

Fig.6Highspeedimagesundertheconditionoftheairmassflowrate196.7g/s,blockageratio0.7,andequivalenceratio1.63

2.2 軸向脈沖爆震的傳播過程

環(huán)形燃燒室的軸向脈沖爆震傳播與脈沖爆震發(fā)動(dòng)機(jī)工作類似，但環(huán)形燃燒室的軸向脈沖爆震模態(tài)不需要額外的點(diǎn)火源，而是依靠激波與火焰的相互作用完成軸向脈沖爆震的引發(fā)。圖7為軸向脈沖爆震的傳播示意圖。從燃燒室入口反射的激波向下游傳播至PCB2處(如圖7(b)中的 Ⅰ過程)，PCB2檢測(cè)出壓力上升信號(hào)(如圖7(a)中的P2壓力尖峰)；向燃燒室出口傳播的激波在燃燒室出口處經(jīng)收斂段反射后向燃燒室入口傳播(對(duì)應(yīng)圖7(b)中的Ⅱ→Ⅲ過程)，在經(jīng)過PCB2時(shí)再次檢測(cè)出壓力上升信號(hào)(如圖7(a)中的P1壓力尖峰)；反傳激波繼續(xù)向燃燒室入口傳播，由于此時(shí)燃燒室入口已被填充一定量的新鮮反應(yīng)物，反傳激波與火焰以及壁面相互作用后，可觸發(fā)新的爆震波并向燃燒室入口傳播(對(duì)應(yīng)圖7(b)中的Ⅳ過程)。由于燃燒室下游沒有新鮮反應(yīng)物維持爆震波傳播，因此爆震波在燃燒室入口反射后解耦，以激波形態(tài)再次向燃燒室出口傳播，重復(fù)過程 Ⅰ和 Ⅱ。圖7(a)的壓力信號(hào)表明，從燃燒室尾部反射的激波強(qiáng)度較弱，如圖中的P1點(diǎn)；而從燃燒室頭部反射的激波強(qiáng)度較強(qiáng)，如圖中的P2點(diǎn)。

(a) 壓力與離子信號(hào)

(b) 軸向脈沖爆震傳播示意圖

考慮軸向爆震在環(huán)形燃燒室內(nèi)的傳播過程，并結(jié)合Anand[11]以及Bykovskii[18]提出的軸向脈沖爆震形成機(jī)制，從出口反傳回燃燒室內(nèi)的激波在靠近燃燒室頭部時(shí)發(fā)展成為爆震波，并伴隨劇烈的燃燒發(fā)光現(xiàn)象。Bykovskii[18]采用速度補(bǔ)償法得到了軸向脈沖爆震狀態(tài)時(shí)的火焰?zhèn)鞑ミ^程，如圖8所示。其中紅色虛線a、b為激波傳播過程，而紅色實(shí)線c、d為火焰?zhèn)鞑ミ^程，可以發(fā)現(xiàn)在激波回傳階段，火焰與激波逐漸加速并耦合，該階段與圖7(b)中Ⅳ階段對(duì)應(yīng)；Anand[11]也通過對(duì)軸向激波時(shí)間與位置的分析，提出反傳的激波在靠近燃燒室頭部時(shí)通過SWACER機(jī)制放大成為爆震波。結(jié)合高速攝影圖6，燃燒室內(nèi)出現(xiàn)周期性強(qiáng)弱交替的發(fā)光現(xiàn)象，對(duì)比圖8的火焰?zhèn)鞑ミ^程，可以推斷出圖6(d)和(h)對(duì)應(yīng)的高亮?xí)r刻處于圖7(b)中的Ⅳ階段，此時(shí)激波加速成爆震波。當(dāng)爆震波傳到燃燒室頭部并消耗完新鮮混合物后，在頭部反射的激波缺少能量供給而逐漸減速。從尾部高速攝影也可以看出，在劇烈發(fā)光之后燃燒室內(nèi)的發(fā)光強(qiáng)度急劇降低，如圖6(d)和(e)所示，表明激波發(fā)展為爆震波的過程中消耗完新鮮反應(yīng)物，并在燃燒室頭部反射后解耦。此外，對(duì)Ⅳ→Ⅰ過程和Ⅱ→Ⅲ過程中的平均波速進(jìn)行了計(jì)算，Ⅳ→Ⅰ過程的平均波速約為1400m/s(取對(duì)應(yīng)的傳播距離為PCB2傳感器距燃燒室入口截面距離的2倍)，而Ⅱ→Ⅲ過程平均波速約為800m/s(取對(duì)應(yīng)的傳播距離為PCB2傳感器距燃燒室出口截面距離的2倍)，這也表明激波在燃燒室頭部進(jìn)行了加速，而在往下游傳播時(shí)速度逐漸下降。

圖8 軸向脈沖爆震模態(tài)下火焰?zhèn)鞑ミ^程[18]

Fig.8Shockandflamepropagationunderlongitudinalpulsedetonationmode[18]

2.3 軸向脈沖爆震模態(tài)的工作頻率

當(dāng)燃燒室出口阻塞比以及質(zhì)量通量滿足一定條件時(shí)，環(huán)形燃燒室將會(huì)出現(xiàn)軸向傳播的爆震波，不同質(zhì)量通量以及出口阻塞比條件下發(fā)動(dòng)機(jī)的工作頻率如圖9所示。對(duì)于出口阻塞比為0.6和0.7的幾何約束條件，在質(zhì)量通量小于一定值時(shí)，發(fā)動(dòng)機(jī)以切向單波模態(tài)工作；而增加質(zhì)量通量，發(fā)動(dòng)機(jī)將從切向單波模態(tài)往軸向脈沖爆震模態(tài)轉(zhuǎn)變，對(duì)于本文工況，轉(zhuǎn)變的臨界質(zhì)量通量處于200kg/(m2·s)附近。Anand[11]對(duì)燃燒室出口添加收斂噴管導(dǎo)致軸向脈沖爆震的現(xiàn)象進(jìn)行了較為詳細(xì)的研究，認(rèn)為軸向脈沖爆震產(chǎn)生時(shí)空氣積氣腔與燃燒室的平均壓力比值需滿足一定條件。而在本文實(shí)驗(yàn)條件下，質(zhì)量通量高于200kg/(m2·s)時(shí)，燃燒室以軸向脈沖爆震模態(tài)工作，這是因?yàn)樵黾淤|(zhì)量通量將導(dǎo)致燃燒室平均壓力上升，燃燒室內(nèi)形成的新鮮混氣層的活性提高，有利于在圖7(b)的Ⅳ過程中觸發(fā)新的爆震波。此外，增加質(zhì)量通量將同時(shí)增加積氣腔與燃燒室內(nèi)的壓力，使得其壓力比值滿足一定條件，與Anand[11]的實(shí)驗(yàn)結(jié)果是相容的。此外，圖9也表明發(fā)動(dòng)機(jī)以軸向脈沖爆震模態(tài)工作時(shí)，其工作頻率隨質(zhì)量通量的增加而增加。

圖9 軸向脈沖工作模態(tài)頻率

爆震波在環(huán)形燃燒室切向傳播時(shí)，由于燃料與氧化劑混合不均勻、爆震波高度隨質(zhì)量流率的變化等因素，很難對(duì)發(fā)動(dòng)機(jī)工作頻率進(jìn)行預(yù)測(cè)。然而，RDE以軸向脈沖爆震模態(tài)工作時(shí)，其平均傳播速度與燃燒產(chǎn)物中的聲速接近，因此可利用線性聲學(xué)理論的假設(shè)對(duì)其傳播頻率進(jìn)行計(jì)算。燃燒室進(jìn)氣面積與出口流通面積相對(duì)于燃燒室的截面積較小，在利用線性聲學(xué)理論[21]時(shí)，可將燃燒室頭部與尾部當(dāng)做固體壁面來對(duì)待，則軸向模態(tài)不穩(wěn)定燃燒的頻率fn可通過式(1)進(jìn)行計(jì)算：

(1)

其中，n為整數(shù)1，代表軸向傳播模態(tài)階數(shù)；c0為燃燒室內(nèi)介質(zhì)的聲速；L為燃燒室長(zhǎng)度，其等效長(zhǎng)度為44mm。

由于實(shí)驗(yàn)采用的當(dāng)量比基本在1.4～1.6附近，采用CEA計(jì)算聲速c0為1056m/s，此時(shí)計(jì)算的一階振頻fn為12 000Hz。采用式(1)計(jì)算的頻率與出口阻塞比為0.6和0.7的實(shí)驗(yàn)結(jié)果對(duì)比如圖10所示。從圖中可以看出，出口阻塞比為0.6與0.7的實(shí)驗(yàn)頻率與采用線性聲學(xué)理論的計(jì)算結(jié)果最大誤差小于15%，這表明采用線性聲學(xué)理論可以對(duì)發(fā)動(dòng)機(jī)出現(xiàn)軸向脈沖爆震模態(tài)的工作頻率進(jìn)行較好的預(yù)測(cè)。

圖10 線性聲學(xué)理論計(jì)算偏差

Fig.10Deviationbetweenthelinearacoustictheoryandtheexperimentalresults

3 結(jié) 論

(1) 在本文研究的2種出口阻塞比(0.6、0.7)條件下，增加空氣混合物的質(zhì)量通量超過200kg/(m2·s)時(shí)，環(huán)形燃燒室內(nèi)將出現(xiàn)軸向傳播的爆震波，且發(fā)動(dòng)機(jī)工作頻率將隨反應(yīng)物質(zhì)量通量的增加而增加；

(2) 發(fā)動(dòng)機(jī)以軸向脈沖爆震模態(tài)工作時(shí)，將伴隨一段切向模態(tài)的聲學(xué)耦合燃燒階段，然后出現(xiàn)軸向傳播的爆震波；

(3) 經(jīng)燃燒室出口收斂型面反射的激波將在燃燒室頭部發(fā)展成為爆震波；

(4) 軸向爆震波在傳播過程中經(jīng)歷解耦與重新發(fā)展的過程，其總體傳播速度與燃燒產(chǎn)物的聲速相當(dāng)，采用線性聲學(xué)理論可對(duì)該模態(tài)下的發(fā)動(dòng)機(jī)工作頻率進(jìn)行較好的預(yù)測(cè)。