基于滾動(dòng)時(shí)間序列SVR的地鐵咽喉區(qū)小凈距隧道圍巖位移預(yù)測

陳濤

基于滾動(dòng)時(shí)間序列SVR的地鐵咽喉區(qū)小凈距隧道圍巖位移預(yù)測

陳濤

(中鐵第四勘察設(shè)計(jì)院集團(tuán)有限公司，湖北 武漢 430063)

為更有效地預(yù)測車輛段咽喉區(qū)小凈距隧道的圍巖位移，采用時(shí)間序列支持向量機(jī)(SVM)方法，基于實(shí)際監(jiān)測數(shù)據(jù)分析小凈距隧道的圍巖沉降，并比較不同SVM參數(shù)下的預(yù)測精度評價(jià)指標(biāo)。研究結(jié)果表明：基于滾動(dòng)時(shí)間序列SVR方法在預(yù)測25 d以內(nèi)的數(shù)據(jù)時(shí)，可以保持足夠的精度，其均方根誤差(RMSE)和平均絕對相對誤差(MAPE)分別能保持在0.3和0.5%以下。SVM輸出的預(yù)測數(shù)據(jù)不僅可以較好地?cái)M合原有數(shù)據(jù)的趨勢，還可以有效地模擬數(shù)據(jù)在小范圍內(nèi)的波動(dòng)，并可較為準(zhǔn)確地判斷襯砌結(jié)構(gòu)沉降是否收斂。因此，使用基于滾動(dòng)時(shí)間序列SVR的圍巖位移預(yù)測方法對隧道和圍巖的位移進(jìn)行預(yù)測可為隧道施工方法的優(yōu)化決策提供參考。

小凈距隧道；支持向量機(jī)；圍巖位移；位移預(yù)測

隧道襯砌結(jié)構(gòu)穩(wěn)定性分析方法有經(jīng)驗(yàn)類比法、應(yīng)力解析法、數(shù)值模擬法和物理試驗(yàn)法四類，對于地鐵車輛段咽喉區(qū)小凈距隧道可以從不同的角度嘗試新的方法進(jìn)行分析。隧道位移是典型的非線性變化的數(shù)據(jù)，具有很強(qiáng)的不確定性[1]，而隧道的施工需要對位移進(jìn)行嚴(yán)密的關(guān)注。從效果上看，數(shù)值模擬通過系統(tǒng)內(nèi)部的力學(xué)演化在宏觀上把握工程的安全性，而各部位實(shí)測的數(shù)據(jù)變化則能反應(yīng)細(xì)節(jié)上的變形特征，如何把握這樣的趨勢用來預(yù)測未來的變形趨勢并將其反饋到施工決策中，是一個(gè)復(fù)雜的非線性回歸問題[2]，借助于機(jī)器學(xué)習(xí)的方法便可以將這樣的問題化解。機(jī)器學(xué)習(xí)是大樣本學(xué)習(xí)機(jī)器，它利用已有數(shù)據(jù)對未來數(shù)據(jù)做出預(yù)測和響應(yīng)，將機(jī)器學(xué)習(xí)應(yīng)用于地鐵車輛段咽喉區(qū)小間距隧道工程中，會(huì)給工程監(jiān)測預(yù)警系統(tǒng)提供另一個(gè)維度的可靠指標(biāo)。本文以實(shí)際工程作為案例，探索基于機(jī)器學(xué)習(xí)對隧道工程的應(yīng)用，利用時(shí)間序列支持向量機(jī)，對隧道節(jié)點(diǎn)位移的變化進(jìn)行預(yù)測。機(jī)器學(xué)習(xí)的核心思想是利用算法從大量已有的資料或者數(shù)據(jù)中找到其中潛在的規(guī)律，并進(jìn)行預(yù)測或者分類。其中人工神經(jīng)網(wǎng)絡(luò)(Artificial Neural Network，即ANN)應(yīng)用到土木工程領(lǐng)域[3?6]是20世紀(jì)80年代末才逐漸浮現(xiàn)。而支持向量機(jī)(Support Vector Machine)[7?8]屬于有導(dǎo)師的學(xué)習(xí)網(wǎng)絡(luò)，是20世紀(jì)90年代中期發(fā)展起來的機(jī)器學(xué)習(xí)技術(shù)，與人工神經(jīng)網(wǎng)絡(luò)不同，前者基于結(jié)構(gòu)風(fēng)險(xiǎn)最小化原理，后者基于經(jīng)驗(yàn)風(fēng)險(xiǎn)最小化原理，它是針對二分類問題提出的，其機(jī)器學(xué)習(xí)思想是同時(shí)最小化經(jīng)驗(yàn)風(fēng)險(xiǎn)和置信范圍，最終解決的是一個(gè)凸二次規(guī)劃問題，從理論上說，得到的將是全局最優(yōu)解，解決了在神經(jīng)網(wǎng)絡(luò)方法中無法避免的局部極值問題。趙洪波等[2]采用基于統(tǒng)計(jì)學(xué)習(xí)理論的支持向量機(jī)結(jié)合時(shí)間序列預(yù)測滑坡的位移，這也是國內(nèi)學(xué)者初次將支持向量機(jī)進(jìn)行巖土領(lǐng)域的應(yīng)用。高瑋等[9]運(yùn)用人工神經(jīng)網(wǎng)絡(luò)的方法對邊坡穩(wěn)定性和滑坡變形進(jìn)行了預(yù)測分析，針對滑坡位移序列預(yù)測將人工神經(jīng)網(wǎng)絡(luò)改進(jìn)為灰色系統(tǒng)—進(jìn)化神經(jīng)網(wǎng)絡(luò)，認(rèn)為此種方法對于滑坡監(jiān)測具有較好的推廣性。此外, 王衛(wèi)東等[10]利用支持向量機(jī)模型對四川省的滑坡災(zāi)害易發(fā)性進(jìn)行區(qū)劃。陳昌富等[11]利用支持向量機(jī)對復(fù)合地基的沉降設(shè)計(jì)進(jìn)行優(yōu)化。周奇才等[12]采用改進(jìn)的支持向量機(jī)對隧道變形進(jìn)行了預(yù)測，6個(gè)測試樣本的預(yù)測均方根誤差最小達(dá)到0.070 5。汪華斌等[13]利用支持向量機(jī)的分類性能對巖爆進(jìn)行了預(yù)測，預(yù)測準(zhǔn)確率達(dá)到100%，即都不會(huì)發(fā)生巖爆。雖然支持向量機(jī)方法已被廣泛應(yīng)用在隧道、地基、邊坡等工程的位移和穩(wěn)定性評價(jià)中，但是使用該方法對小凈距隧道圍巖的位移進(jìn)行預(yù)測的研究報(bào)道較為少見。小凈距隧道具有隧道間距極小、圍巖位移易受施工影響等特點(diǎn)，圍巖位移的影響因素復(fù)雜，而支持向量機(jī)方法在解決類似問題上具有其獨(dú)特的優(yōu)勢，因此，本文使用支持向量機(jī)方法對小凈距隧道的圍巖位移進(jìn)行預(yù)測并對該方法的適用性進(jìn)行研究。

1 支持向量機(jī)的回歸與分類

支持向量機(jī)所解決的分類問題是通過尋找最優(yōu)超平面的方法，條件是使得間隔最大，反之，如果使得間隔最小就可以尋找到一個(gè)回歸平面，因此，支持向量機(jī)也被推廣到解決回歸問題上，即支持向量回歸(Support Vector Regression)，其數(shù)據(jù)集形式為(1,1),…,(x,y)，x可以是多維向量，即預(yù)測因子，y是回歸中的預(yù)測對象。SVR中最為普遍和有效的算法是基于-insensitive損失函數(shù)的-SVR算法。對于-SVR的應(yīng)用包括線性回歸和非線性回歸。

對于線性問題，約束條件為：

對于非線性回歸問題，首先也是使用一個(gè)非線性映射把樣本點(diǎn)映射到高維特征空間，再在高維特征空間進(jìn)行回歸，關(guān)鍵也是核函數(shù)的采用，用(,)來表示核函數(shù)，非線性優(yōu)化方程的目標(biāo)函數(shù)為：

其約束條件同樣為式(2)。

此處引入KKT(Karush–Kuhn–Tucker)條件，即對一個(gè)非線性問題有最優(yōu)解解法的充要條件，具體分為3種情況無約束條件，等式約束條件和不等式約束條件3種情況。對于等式約束條件，定義拉格朗日函數(shù)

其中：是各個(gè)約束條件的待定系數(shù)，h()為約束條件函數(shù)。再求解各變量的偏導(dǎo)方程：

如果有個(gè)約束條件，則得到+1個(gè)方程，求出的方程組的解就可能是最優(yōu)化值，將結(jié)果代入原方程驗(yàn)證就可得到解。因此對于式(5)~(7)得到：

上式即可解得的計(jì)算式。令

當(dāng)為非零向量時(shí)，其對應(yīng)的訓(xùn)練樣本就是支持向量。

SVR的創(chuàng)建過程需要對參數(shù)進(jìn)行優(yōu)化，參數(shù)的好壞直接決定了SVR性能的高低。

2 滾動(dòng)時(shí)間序列SVR創(chuàng)建及參數(shù)優(yōu)化

現(xiàn)實(shí)中許多情況下的數(shù)據(jù)通常就是以時(shí)間序列的形式給出，以時(shí)間序列法進(jìn)行預(yù)測，就是要尋找前個(gè)樣本與第+1個(gè)樣本的關(guān)系，即如下式的關(guān)系，輸出數(shù)據(jù)隨著時(shí)間向前推動(dòng)，訓(xùn)練樣本不斷地變化，但是樣本內(nèi)元素量不會(huì)變動(dòng)。

(8)

對向量機(jī)預(yù)測效果的評價(jià)由均方根誤差(Mean Absolute Error)、平均絕對相對誤差(Mean Absolute Percentage Error)、3項(xiàng)指差(Root Mean Square Error)、平均絕對誤標(biāo)來評判，計(jì)算公式如下：

均方根誤差(RMSE)是觀測值與真值偏差的平方和與觀測次數(shù)比值的平方根，用來衡量觀測值同真值之間的偏差。平均絕對誤差(MAE)是絕對誤差的平均值，能更好地反映預(yù)測值誤差的實(shí)際情況。平均絕對相對誤差(MAPE)是絕對誤差與真實(shí)值的比值的平均值，反映了誤差與真實(shí)值的比例關(guān)系。這3項(xiàng)評價(jià)指標(biāo)越小，說明預(yù)測精度越高。

3 基于滾動(dòng)時(shí)間序列SVR的圍巖位移預(yù)測

3.1 工程概況

深圳市城市軌道交通8號線一期工程8133標(biāo)段—望基湖停車場及出入線綜合工程，位于鹽田區(qū)鹽排高速以西及深圳外國語學(xué)校西南側(cè)的梧桐山山區(qū)內(nèi)。望基湖停車場出入線分別由鹽田港站站后及深外站站前接出，鹽田港站站后接出線下穿林長隧道，左右線(即入場線和出場線)交匯后以2.59 m的極小間距并行為咽喉區(qū)的單洞雙線隧道，以=300 m半徑自南向北轉(zhuǎn)入望基湖停車場。咽喉區(qū)，顧名思義即為地鐵線路的咽喉部位，對于地鐵運(yùn)營的意義重大，本工程中出現(xiàn)的連接咽喉區(qū)部位的小間距隧道無論是從施工難度上，還是從設(shè)施的重要性上都是值得注意和研究的。圖1為單洞雙線隧道與左右線隧道結(jié)構(gòu)設(shè)計(jì)圖。

單位：mm

3.2 位移預(yù)測結(jié)果分析

咽喉區(qū)小凈距隧道在施工過程中有嚴(yán)密的監(jiān)測預(yù)警系統(tǒng)，摘取WRDK2+094斷面從2018年7月10日起的拱頂沉降變化監(jiān)測數(shù)據(jù)，繪成折線圖如圖2所示，共150 d的累計(jì)沉降數(shù)據(jù)，該斷面位于入場線隧道與單洞雙線隧道的交界處，具有代表意義，可以看出拱頂?shù)某两道塾?jì)值一直隨著時(shí)間在增大，在小范圍內(nèi)存在一定的波動(dòng)。

圖2 150 d內(nèi)WRDK2+094斷面沉降值變化曲線

通過采用時(shí)間序列的方式構(gòu)建預(yù)測模型，即輸入訓(xùn)練樣本為滾動(dòng)的時(shí)間序列{X}={1,2,…,X}，樣本的大小由主觀確定，同樣，最終測試樣本數(shù)也可以自定義為。由于和值的不確定性，預(yù)測模型的性能也會(huì)隨之波動(dòng)，將數(shù)種情況下計(jì)算結(jié)果的評價(jià)指標(biāo)列入表1中，可以得出結(jié)論：當(dāng)預(yù)測范圍為25 d以內(nèi)的數(shù)據(jù)時(shí)，通過調(diào)整時(shí)間序列長度便可以得到非常高的預(yù)測精度，通常均方根誤差(RMSE)維持在0.5以內(nèi)，平均絕對相對誤差(MAPE)在3%以內(nèi)，但是對于不同訓(xùn)練集SVR會(huì)有不同的性能，不能保證這樣的預(yù)測精度適用于任意場景。當(dāng)預(yù)測天數(shù)達(dá)到30 d時(shí)，預(yù)測精度明顯下降，多數(shù)情況下均方根誤差(RMSE)超過了1，平均絕對相對誤差(MAPE)全部超過5%，最高為20.99%。

通過對比輸出樣本的折線圖將更加直觀，如圖3和圖4分別為=21時(shí)預(yù)測25 d和30 d的輸出樣本。另外，觀察表中數(shù)據(jù)可以看出，對于預(yù)測天數(shù)小的情況，值偏小時(shí)預(yù)測結(jié)果會(huì)相對更加準(zhǔn)確，預(yù)測天數(shù)大時(shí)，值也應(yīng)盡量選擇更大。由于時(shí)間序列本身的復(fù)雜性，不同的時(shí)間序列攜帶的信息差異巨大，所以值的選取應(yīng)該根據(jù)具體的時(shí)間序列及其測試結(jié)果進(jìn)行調(diào)整。

表1 沉降預(yù)測計(jì)算結(jié)果評價(jià)指標(biāo)

圖3 m=25，k=21時(shí)計(jì)算結(jié)果對比

圖4 m=30，k=21時(shí)計(jì)算結(jié)果對比

一般來說，精度預(yù)測精度并不一定隨著預(yù)測天數(shù)的減少而提高，而是與測試樣本的變異系數(shù)(C.V)有關(guān)，計(jì)算各個(gè)區(qū)間測試樣本的變異系數(shù)(C.V)如表2所示，取分別為13，17和21時(shí)得到的平均絕對相對誤差(MAPE)，計(jì)算得到C.V與MAPE的相關(guān)系數(shù)分別為0.943 0，0.855 6和0.946 8，均處于強(qiáng)相關(guān)的范圍，說明測試樣本的變化越大，變異系數(shù)越大，則預(yù)測難度增大從而精度降低，極端情況下，為1，預(yù)測天數(shù)為1時(shí)，此時(shí)測試集變異系數(shù)(C.V)為0，計(jì)算結(jié)果的3項(xiàng)評價(jià)指標(biāo)分別為0.105 8%，0.105 8%和0.89%，預(yù)測精度極高。

同樣方法對150 d內(nèi)的WRDK2+094斷面拱肩部位的水平收斂(絕對值)監(jiān)測數(shù)據(jù)進(jìn)行處理，如圖5所示，為收斂值的折線圖。

表2 各預(yù)測區(qū)間變異系數(shù)

圖5 150 d內(nèi)WRDK2+094斷面拱肩水平收斂值

(絕對值)變化曲線

Fig. 5 Variation curve of horizontal convergence value (absolute value) of WRDK2+094 cross section spandrel within 150 days

根據(jù)表3中的結(jié)果可以看出，對于同樣的測試樣本，不同的值決定了SVR的性能，當(dāng)=1時(shí)，MAPE均小于5%，當(dāng)=17時(shí)，MAPE最小為1.42%，最大達(dá)38.28%。故對于不同的預(yù)測天數(shù)時(shí)，應(yīng)尋找最優(yōu)的值以使SVR的性能最優(yōu)。圖6和圖7分別為=25，=9和=30，=9時(shí)的輸出樣本的曲線圖，可以看出訓(xùn)練輸出值仍可以很好的擬合訓(xùn)練樣本值，在測試樣本中，同樣可以擬合出大致的趨勢，但是在峰值處還不能很好的還原。

預(yù)測結(jié)果對于值選取的敏感表明對于SVR性能至關(guān)重要，由于SVR的表現(xiàn)依賴于樣本所給出的信息，而值則決定了能夠讀取的信息量大小，對于越長的預(yù)測時(shí)間，SVR訓(xùn)練所應(yīng)采用的值應(yīng)更大，從而使得更多的信息被SVR所讀取，以保證預(yù)測結(jié)果的精度。

如圖6為當(dāng)=21時(shí)計(jì)算得到的25 d的沉降預(yù)測值與監(jiān)測值的對比，可以看出預(yù)測沉降值基本收斂于12 mm，且此后的151~167 d的監(jiān)測數(shù)據(jù)也驗(yàn)證了這一結(jié)論。因此沉降的收斂值遠(yuǎn)低于現(xiàn)場的報(bào)警值24 mm，再將=9時(shí)計(jì)算得到的30 d拱肩水平收斂預(yù)測值與監(jiān)測值進(jìn)行對比，如圖7，可以看出，預(yù)測值有一定的波動(dòng)，但也逐漸的有收斂趨勢，且未大于10 mm，151~167 d的實(shí)測值也收斂于10 mm，同樣遠(yuǎn)低于報(bào)警值24 mm。

表3 水平收斂值預(yù)測結(jié)果評價(jià)指標(biāo)

圖6 沉降預(yù)測值與實(shí)測值的對比

圖7 拱肩水平收斂預(yù)測值與實(shí)測值的對比

4 結(jié)論

1) 采用時(shí)間序列法創(chuàng)建SVR預(yù)測圍巖位移，其輸出結(jié)果不僅較好地?cái)M合原數(shù)據(jù)的趨勢，并可還原實(shí)測數(shù)據(jù)小范圍內(nèi)的波動(dòng)，其預(yù)測精度也保持在較高的水平。

2) 采用時(shí)間序列法SVR對拱肩的水平收斂值(絕對值)進(jìn)行預(yù)測時(shí)，對于不同的預(yù)測天數(shù)需要對應(yīng)地調(diào)整值的大小以獲得較高的預(yù)測精度和較好的擬合效果，另外由于時(shí)間序列包含信息的不確定性，SVR所能提供的預(yù)測精度會(huì)有所不同。

3) 在已有一定數(shù)據(jù)積累的情況下，通過調(diào)整值可以獲得好的預(yù)測結(jié)果，預(yù)測值與后期實(shí)測數(shù)據(jù)的對比結(jié)果表明，本文所用方法可以較好地預(yù)測小凈距隧道圍巖位移的發(fā)展趨勢和收斂值。

[1] 梁慶國, 樊純壇, 趙志強(qiáng). 大斷面鐵路隧道室內(nèi)模型試驗(yàn)與分析研究[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2017, 14(9): 1915?1924. LIANG Qingguo, FAN Chuntan, ZHAO Zhiqiang. Analysis and research on laboratory model test for large section railway tunnel[J]. Journal of Railway Science and Engineering, 2017, 14(9): 1915?1924.

[2] 趙洪波, 馮夏庭. 非線性位移時(shí)間序列預(yù)測的進(jìn)化-支持向量機(jī)方法及應(yīng)用[J]. 巖土工程學(xué)報(bào), 2003(4): 468?471. ZHAO Hongbo, FENG Xiating. Study and application of genetic-support vector machine for nonlinear displacement time series forecasting[J]. Chinese Journal of Geotechnical Engineering, 2003(4): 468?471.

[3] Flood I, Kartam N. Neural networks in civil engineering. I: Principles and understanding[J]. Journal of Computing in Civil Engineering, 1994, 8(2): 131?148.

[4] 高浪, 謝康和. 人工神經(jīng)網(wǎng)絡(luò)在巖土工程中的應(yīng)用[J].土木工程學(xué)報(bào), 2002(4): 77?81. GAO Lang, XIE Kanghe. Application of artificial neural networks to geotechnical engineering[J]. China Civil Engineering Journal, 2002(4): 77?81.

[5] Dede T, Kankal M, Vosoughi A R, et al. Artificial intelligence applications in civil engineering[J]. Advances in Civil Engineering, 2019: 1?3.

[6] Knezevic M, Cvetkovska M, Hanák T, et al. Artificial neural networks and fuzzy neural networks for solving civil engineering problems[J]. Complexity, 2018: 1?2.

[7] CHEN B, FU X, GUO X, et al. Zoning elastic modulus inversion for high arch dams based on the PSOGSA- SVM method[J]. Advances in Civil Engineering, 2019(9): 1?13.

[8] Kirts S, Nam B H, Panagopoulos O P, et al. Settlement prediction using support vector machine (SVM)-based compressibility models: A case study[J]. International Journal of Civil Engineering, 2019: 1?11.

[9] 高瑋, 馮夏庭. 基于灰色—進(jìn)化神經(jīng)網(wǎng)絡(luò)的滑坡變形預(yù)測研究[J]. 巖土力學(xué), 2004(4): 514?517. GAO Wei, FENG Xiating. Study on displacement predication of landslide based on grey system and evolutionary neural network[J]. Rock and Soil Mechanics, 2004(4): 514?517.

[10] 王衛(wèi)東, 劉攀, 龔陸. 基于支持向量機(jī)模型的四川省滑坡災(zāi)害易發(fā)性區(qū)劃[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2019, 16(5): 1194?1200. WANG Weidong, LIU Pan, GONG Lu. Landslide susceptibility mapping of Sichuan province based on support vector machine[J]. Journal of Railway Science and Engineering, 2019, 16(5): 1194?1200.

[11] 陳昌富, 李欣, 龔曉南, 等. 基于支持向量機(jī)沉降代理模型復(fù)合地基優(yōu)化設(shè)計(jì)方法[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2018, 15(6): 1424?1429. CHEN Changfu, LI Xin, GONG Xiaonan, et al. Composite foundation’s optimization design of settlement’s surrogate models based on support vector machines[J]. Journal of Railway Science and Engineering, 2018, 15(6): 1424?1429.

[12] 周奇才, 范思遐, 趙炯, 等. 基于改進(jìn)的支持向量機(jī)隧道變形預(yù)測模型[J]. 鐵道工程學(xué)報(bào), 2015, 32(3): 67? 72. ZHOU Qicai, FAN Sixia, ZHAO Jiong, et al. Tunnel deformation prediction based on modified support vector machine[J]. Journal of Railway Engineering Society, 2015, 32(3): 67?72.

[13] 汪華斌, 盧自立, 邱杰漢, 等. 基于粒子群算法優(yōu)化支持向量機(jī)的巖爆預(yù)測研究[J]. 地下空間與工程學(xué)報(bào), 2017, 13(2): 364?369. WANG Huabin, LU Zili, QIU Jiehan, et al. Prediction of rock burst by improved particle swam optimization based support vector machine[J]. Chinese Journal of Underground Space and Engineering, 2017, 13(2): 364? 369.

Prediction of small net distance tunnel surrounding rock displacement in subway throat area based on rolling time series SVR

CHEN Tao

(China Railway Siyuan Survey and Design Group Co., Ltd., Wuhan 430063, China)

In order to more effectively predict the displacement of the surrounding rock of the small clear-distance tunnel in the throat area of the vehicle section, this paper used the time series support vector machine (SVM) method to analyze the surrounding rock settlement of the small clear-distance tunnel based on the actual monitoring data and compared the different SVM parameters Index of prediction accuracy. The results show that the rolling time series SVR method can maintain sufficient accuracy when predicting data within 25 days, and its root mean square error (RMSE) and average absolute relative error (MAPE) can be kept below 0.3 and 0.5%, respectively. The prediction data output by SVM can not only well fit the trend of the original data, but also effectively simulate the fluctuation of the data in a small range, and can more accurately judge whether the settlement of the lining structure has converged. Therefore, using the prediction method of surrounding rock displacement based on rolling time series SVR to predict the displacement of the tunnel and surrounding rock can provide a reference for the optimization decision of the tunnel construction method.

small net distance tunnel; support vector machine; surrounding rock displacement; displacement prediction

U455

A

1672 ? 7029(2020)09 ? 2338 ? 08

10.19713/j.cnki.43?1423/u. T20191066

2019?12?02

國家自然科學(xué)基金資助項(xiàng)目(51708564)

陳濤(1982?)，男，湖北荊州人，高級工程師，從事地下工程勘察設(shè)計(jì)與研究；E?mail：tsyct@126.com

(編輯 陽麗霞)