國道G318藏東段碎屑斜坡分類與穩(wěn)定性評(píng)判*

葉唐進(jìn) 謝 強(qiáng) 王 鷹

(①西南交通大學(xué)地球科學(xué)與環(huán)境工程學(xué)院 成都 610031)(②西藏大學(xué)工學(xué)院 拉薩 850000)

0 引 言

國道G318藏東段，由于強(qiáng)烈構(gòu)造活動(dòng)、陡峻地形以及特殊的氣象條件(中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所，西藏自治區(qū)交通科學(xué)研究所， 1999)，沿線的溜砂、滾石等碎屑斜坡地質(zhì)災(zāi)害數(shù)量多、分布廣、發(fā)生頻繁、治理難度較大，嚴(yán)重威脅過往行人及車輛的安全(楊志法等， 2002)。公路部門采用擋土墻、SNS主動(dòng)和被動(dòng)防護(hù)網(wǎng)等多種措施，效果一直不甚理想，其危害沒有得到根本性的解決。碎屑斜坡是指在特殊地質(zhì)、氣候等作用下形成的，以碎石、礫石及中粗砂為主，穩(wěn)定性極為敏感的松散堆積物斜坡。碎屑斜坡在藏東高原高山峽谷區(qū)公路邊坡段往往集中發(fā)育并導(dǎo)致災(zāi)害頻發(fā)。

針對(duì)碎屑斜坡的研究，最早可以追溯到1938年， Sharpe(1938)首次提出了碎屑流的概念，后來，Clark(1987)針對(duì)碎屑斜坡的形成因素、運(yùn)動(dòng)特征和堆積形態(tài)展開了相關(guān)研究。2000年后，部分學(xué)者在碎屑斜坡自組織臨界規(guī)律等方面展開了研究(Dhar，2006； Mohanty et al.，2007)。而國內(nèi)的研究，最早始于1995年，羅德富等(1995)最先對(duì)國道G318西藏境內(nèi)碎屑斜坡展開研究； 隨后，袁廣祥等(2012)對(duì)川藏公路第四紀(jì)堆積體的成因及其分布規(guī)律做了研究； 王成華等(2007)對(duì)粒狀碎屑溜砂坡的形成和基本特征做了較為系統(tǒng)的分析研究； 湯明高等(2012)對(duì)國道G318典型地質(zhì)災(zāi)害成因機(jī)制及防治對(duì)策作了分析； 張路青等(2004a)對(duì)川藏公路南線八宿—林芝段滾石災(zāi)害進(jìn)行了調(diào)查和評(píng)價(jià)，對(duì)滾石擊中行人和車輛進(jìn)行了風(fēng)險(xiǎn)評(píng)價(jià)分析，但研究過程中將不同成分、成因及危害程度的碎屑斜坡混為一體； 齊得旭(2015)對(duì)散粒體斜坡穩(wěn)定性的分析方法進(jìn)行了實(shí)驗(yàn)研究。綜上所述，國內(nèi)外的研究都主要集中在個(gè)別典型碎屑斜坡的研究上，目前對(duì)碎屑斜坡災(zāi)害的重視程度仍不夠、調(diào)查不全面、研究方法較簡(jiǎn)單、研究深度不足，導(dǎo)致研究成果少、后期治理效果較差，道路的安全暢通沒有得到根本性的改善。

在研究區(qū)域上，主要集中在阿爾卑斯山脈，天山山脈和橫斷山脈(王成華等， 2007； 張?jiān)诺龋?2008)，也有喜馬拉雅山脈(Kanungo et al.，2013)和川西高原等地(曾慶利等， 2018)，但是部分學(xué)者主要研究滑坡引起的碎屑流災(zāi)害，其對(duì)象與G318藏東段碎屑斜坡的概念及類型存在明顯的差異。從碎屑斜坡的研究方法來看，主要針對(duì)單個(gè)碎屑斜坡的野外調(diào)查、理論分析、數(shù)值模擬和模型試驗(yàn)等方法。近年來，針對(duì)斜坡穩(wěn)定性利用現(xiàn)代數(shù)學(xué)工具分析逐漸增多(詹威威等， 2017)，模糊數(shù)學(xué)、灰色理論、隨機(jī)森林和BP神經(jīng)網(wǎng)絡(luò)等數(shù)學(xué)方法得到了廣泛應(yīng)用。在地質(zhì)災(zāi)害影響因素分析中，丁繼新等(2005)利用信息熵分析了川藏公路然烏—魯朗段地質(zhì)災(zāi)害影響因素，葉唐進(jìn)(2017)利用信息熵分析了甲瑪?shù)V區(qū)碎屑斜坡的失穩(wěn)因素，其研究取得了較好的效果。在邊坡穩(wěn)定性評(píng)價(jià)中，支持向量機(jī)(SVM)能夠很好地解決小樣本、非線性、高維數(shù)的問題，趙洪波等(2003)利用支持向量機(jī)函數(shù)擬合對(duì)邊坡穩(wěn)定性估計(jì)進(jìn)行研究，取得了理想的效果。因此，針對(duì)碎屑斜坡影響因素分析和穩(wěn)定性判別，信息熵和支持向量機(jī)是一種不錯(cuò)的選擇。

國道G318藏東段是進(jìn)藏物資和人員往來的重要交通要道，對(duì)鞏固國防、發(fā)展經(jīng)濟(jì)、改善民生以及一帶一路等具有重要的戰(zhàn)略價(jià)值。因此，如何保證該路段全天候安全暢通顯得尤為重要。從2014年8月至2017年8月，作者等人歷時(shí)3年，通過野外調(diào)查獲取了豐富的第一手資料，在此基礎(chǔ)上，根據(jù)斜坡碎屑物質(zhì)成分、顆粒大小和運(yùn)動(dòng)特征等，對(duì)斜坡碎屑進(jìn)行了類型劃分，并闡述了各類碎屑斜坡的基本特征； 然后，再運(yùn)用信息熵法和FLAC計(jì)算分析了主要影響因素與碎屑斜坡穩(wěn)定性的相關(guān)關(guān)系，最后，結(jié)合信息熵和FLAC計(jì)算分析的主要影響因素，建立基于支持向量機(jī)對(duì)碎屑斜坡的穩(wěn)定性進(jìn)行了評(píng)判，以期為今后的碎屑斜坡穩(wěn)定性判別和維修整治參考。

圖 1 國道G318藏東段碎屑斜坡分布圖Fig. 1 The spatial distribution map of debris slopes in eastern Tibet section of G318ihighway

1 碎屑斜坡野外調(diào)查

1.1 調(diào)查背景

從2014年8月至2017年8月，作者等人歷時(shí)3年，對(duì)藏東段1277ikm的碎屑斜坡進(jìn)行了全面調(diào)查，采集110個(gè)碎屑斜坡的數(shù)據(jù)，基本涵蓋了所有的碎屑斜坡類型。從調(diào)查統(tǒng)計(jì)結(jié)果來看(表 1)，以崩坡積物為主的溜砂坡分布最為集中，主要分布于林芝至八宿的帕隆藏布河谷； 以沖洪積物為主的滾石坡數(shù)量最多，主要分布于拉薩至林芝的尼洋河谷； 以殘坡積物為主的碎石坡分布最廣，從八宿至巴塘的怒江、瀾滄江河谷均有分布，具體分布如圖 1所示。碎屑斜坡的危害嚴(yán)重、治理難度較大，治理效果較差。

表 1 國道G318藏東段碎屑斜坡統(tǒng)計(jì)表Table 1 Statistics of debris slopes from eastern Tibet section of G318 highway

1.2 碎屑斜坡調(diào)查內(nèi)容

碎屑斜坡的野外調(diào)查內(nèi)容包括：(1)斜坡的概況(類型、坐標(biāo)、高程、地名)； (2)地質(zhì)背景(微地貌、地層巖性、地質(zhì)構(gòu)造)； (3)地理環(huán)境(降水量、水文、土地利用類型)； (4)結(jié)構(gòu)特征(名稱、密實(shí)度、稠度、下伏基巖的產(chǎn)狀、埋深、地層巖性及年代)； (5)地下水(類型、埋深、出露點(diǎn)、補(bǔ)給)； (6)斜坡的變形特征(類型、部位、特征、產(chǎn)生時(shí)間)； (7)斜坡可能的失穩(wěn)因素、目前的穩(wěn)定狀態(tài)、已造成的危害、潛在的危害等； (8)監(jiān)測(cè)和防治建議。在野外調(diào)查基礎(chǔ)上繪制了相關(guān)的平面圖、剖面圖、碎屑斜坡與公路的位置關(guān)系等圖件。

2 碎屑斜坡分類、基本特征及失穩(wěn)模式

2.1 碎屑斜坡分類

根據(jù)野外調(diào)查資料，結(jié)合碎屑斜坡成分、顆粒大小、膠結(jié)狀態(tài)、失穩(wěn)運(yùn)動(dòng)方式，對(duì)公路的危害程度以及整治的難易程度，將碎屑斜坡細(xì)分為：溜砂坡、滾石坡以及碎石坡。具體的類型及劃分依據(jù)如表 2所示。

表 2 碎屑斜坡分類及劃分依據(jù)Table 2 Dividing basis and basic classification of debris slopes

圖 2 3類碎屑斜坡照片F(xiàn)ig. 2 The photos of three types of the debris slopes(a) sliding sand slope(b) rolling rock slope(c) and gravel slopea. 溜砂坡; b. 滾石坡; c. 碎石坡

2.2 碎屑斜坡基本特征

溜砂坡主要集中在宗壩至然烏段，因中粗砂和礫石為主而得名溜砂坡，為二長(zhǎng)花崗巖、片麻巖的風(fēng)化產(chǎn)物(圖 2a)。在地震、振動(dòng)等外力作用下，斜坡上的砂粒及零星碎石則向下滾動(dòng)，極易砸毀擋墻、路面以及護(hù)欄等，威脅過往的車輛及行人。當(dāng)降雨融雪的地表水與砂?；旌铣伤傲飨蛳铝鲃?dòng)(葉唐進(jìn)等， 2016)，極易沖毀防護(hù)網(wǎng)，并堆積于公路之上，阻礙公路暢通。此類碎屑斜坡雖然規(guī)模不大，但數(shù)量眾多，突發(fā)性較強(qiáng)，而且連續(xù)成群出現(xiàn)。

滾石坡主要分布于尼洋河、帕隆藏布、怒江以及瀾滄江河谷路段，其主要成分為礫石、卵石等沖積物，屬于河谷階地斜坡(圖 2b)。由于公路工程建設(shè)或河水侵蝕坡腳，導(dǎo)致巖土顆粒裸露，當(dāng)降雨融雪時(shí)，裸露的巖土顆粒不斷剝落以及局部滑坡解體形成滾石災(zāi)害，砸毀擋墻、路面及防護(hù)網(wǎng)，阻斷交通，嚴(yán)重威脅過往行人及車輛安全。

碎石坡主要分布于八宿至芒康橫斷山脈干旱、半干旱的高原高山峽谷路段，主要為板巖、千枚巖等風(fēng)化的殘坡積物，以礫石和碎石顆粒為主(圖 2c)。當(dāng)含水量變化或振動(dòng)時(shí)，處于極限穩(wěn)定的顆粒失穩(wěn)運(yùn)動(dòng)(葉唐進(jìn)， 2017)，同時(shí)觸發(fā)其下游碎石運(yùn)動(dòng)形成“碎石雨”，并伴隨陣陣異響和滾滾粉塵。從而沖毀公路防護(hù)網(wǎng)、擋墻等設(shè)施，砸毀路面及過往車輛，砸死砸傷過往行人，極易阻斷公路暢通。該類碎屑災(zāi)害突發(fā)性極強(qiáng)，且能量大、速度快、具有一定規(guī)模。

2.3 碎屑斜坡失穩(wěn)模式

通過野外調(diào)查和統(tǒng)計(jì)分析，碎屑坡的3個(gè)亞類都有自己的失穩(wěn)模式(圖 3)。溜砂坡和碎石坡的破壞模式基本相同：形成區(qū)和流通區(qū)在降雨、融雪、地震以及振動(dòng)等條件下，極易形成顆粒滾落，且它們的堆積區(qū)在地震及振動(dòng)下也會(huì)形成顆粒滾落； 不同的是溜砂坡堆積區(qū)在降雨融雪條件下形成水砂流，而碎石坡則易形成滑坡。滾石坡只有形成區(qū)和堆積區(qū)，形成區(qū)有凹腔的斜坡在降雨、融雪、地震以及振動(dòng)等條件下極易形成滑坡，形成區(qū)沒有凹腔和堆積區(qū)在降雨、融雪、地震以及振動(dòng)等條件下均易形成顆粒滾落。

圖 3 碎屑斜坡的失穩(wěn)模式Fig. 3 The failure mode of debris slopes

3 碎屑斜坡的穩(wěn)定性評(píng)判

3.1 主要影響因素分析

碎屑斜坡災(zāi)害形成因素有：碎屑物、構(gòu)造、坡度、坡高、坡長(zhǎng)、植被、風(fēng)化； 誘發(fā)因素有：工程活動(dòng)、降雨融雪、河水侵蝕、地震、振動(dòng)、凍融、風(fēng)吹。通過野外調(diào)查分析，得出10個(gè)主要影響因素分別為坡度、坡高、碎屑物、風(fēng)化、植被、振動(dòng)、地震、河水侵蝕、工程開挖、降雨融雪，為了定量計(jì)算主要影響因素大小順序，利用信息熵(羅進(jìn)， 2014)對(duì)其進(jìn)行計(jì)算分析。將10個(gè)主要影響因素的評(píng)價(jià)指標(biāo)分為5級(jí)并賦值(丁繼新等， 2005)，通過專家打分，其評(píng)價(jià)指標(biāo)的值都在0～1之間， 0表示作用最小， 1表示作用最大。碎屑斜坡穩(wěn)定性評(píng)價(jià)指標(biāo)分級(jí)及標(biāo)準(zhǔn)化詳見表 3。

表 3 碎屑斜坡穩(wěn)定性評(píng)價(jià)指標(biāo)分級(jí)及標(biāo)準(zhǔn)化表Table 3 Classification and standardization of evaluation factors of debris slope stability

通過打分其評(píng)價(jià)指標(biāo)已經(jīng)全部變成無量綱值。再利用式(1)對(duì)隨機(jī)抽取的20個(gè)碎屑斜坡進(jìn)行標(biāo)準(zhǔn)化處理。

(1)

式中，Xi， j表示第i(i=1， 2，…，a)個(gè)樣本的j(j=1， 2，…，b)項(xiàng)評(píng)價(jià)指標(biāo)的值，Xmax(i)、Xmin(i)分別為最大值和最小值，然后根據(jù)式(2)求第j項(xiàng)評(píng)價(jià)指標(biāo)的信息熵。

(2)

表 4 碎屑斜坡穩(wěn)定性影響因素的信息熵及權(quán)重Table 4 Entropy and weight of evaluation factors of debris slope stability

(3)

式中，Wj為第j項(xiàng)評(píng)價(jià)指標(biāo)的權(quán)重，0≤Wj≤1，評(píng)價(jià)指標(biāo)的權(quán)重越大，表示該指標(biāo)在碎屑斜坡穩(wěn)定性中的影響因素就越大。

由表 4可以得出，各評(píng)價(jià)指標(biāo)的權(quán)重從大到小依次為：降雨融雪、平均坡度、碎屑物、平均坡高、植被覆蓋率、工程開挖、地震、風(fēng)化等級(jí)、振動(dòng)和河水侵蝕。為了驗(yàn)證信息熵分析的準(zhǔn)確性，保證穩(wěn)定性判別可靠，下面將繼續(xù)進(jìn)行數(shù)值模擬分析。

3.2 影響因素與穩(wěn)定性的數(shù)值模擬

表 5 碎屑斜坡數(shù)值模擬的概化模型參數(shù)取值Table 5 Parameter value of generalized model of numerical simulation of debris slope

數(shù)值模擬的關(guān)鍵在于模型準(zhǔn)確性驗(yàn)證，現(xiàn)將數(shù)值模擬的云圖與野外斜坡進(jìn)行對(duì)比，從圖 4a和圖4b對(duì)比不難看出，碎屑斜坡的破壞方式，應(yīng)力集中狀態(tài)基本一致，說明了數(shù)值模擬的準(zhǔn)確性。

圖 4 數(shù)值模擬云圖與野外照片對(duì)比Fig. 4 Numerical simulation of cloud compared with wild photos(a) images of shear strain and (b) debris slope photographa. 剪應(yīng)變?cè)茍D； b. 碎屑斜坡照片

現(xiàn)對(duì)碎屑斜坡的降雨融雪、坡度、坡高、工程建設(shè)、植被5個(gè)因數(shù)取值進(jìn)行數(shù)值模擬，各因素的取值如表 6。數(shù)值模擬因數(shù)取值時(shí)，其中坡度、坡高和降雨量的變化直接影響到斜坡的安全系數(shù)； 工程建設(shè)主要涉及削坡，如果削坡坡度一致，工程建設(shè)則可以用削坡的百分比來表示； 植被主要影響碎屑斜坡地表水滲流和坡面巖土顆粒的連接，因此，按照植被的覆蓋率來表示。

表 6 碎屑斜坡數(shù)值模擬因素變化取值Table 6 Value of factors of numerical simulation of debris slope

根據(jù)數(shù)值模擬的結(jié)論，利用式(1)對(duì)橫坐標(biāo)的降水、坡度、坡高、工程建設(shè)(削坡)、植被5個(gè)因素進(jìn)行歸一化處理，然后再分別比較其安全系數(shù)的變化(圖 5)。

圖 5 各影響因素的安全系數(shù)曲線Fig. 5 Safety factor curve of various influencing factors

從圖 5不難得出，碎屑斜坡各因素的斜率大小分別為：降雨融雪(0.8)>平均坡度(0.66)>平均坡高(0.52)>植被(0.46)>工程開挖(0.43)。其斜率越大，對(duì)碎屑斜坡的穩(wěn)定性影響也越大，反之則越小，因此，其結(jié)論與信息熵分析結(jié)論基本一致。結(jié)合野外調(diào)查、信息熵分析和數(shù)值模擬計(jì)算結(jié)論，現(xiàn)將前6個(gè)主要因子利用支持向量機(jī)(SVM)建立模型對(duì)其穩(wěn)定性進(jìn)行判別。

3.3 穩(wěn)定性判別

國道G318藏東段調(diào)查收集了110個(gè)碎屑斜坡實(shí)例，首先，將無效數(shù)據(jù)篩除，即利用箱線圖把離群值、極端值、模糊值和缺失值篩除。然后，從剩余的106個(gè)邊坡中隨機(jī)選取了78個(gè)樣本作為模型的學(xué)習(xí)樣本，其中破壞樣本43個(gè)，穩(wěn)定樣本35個(gè)。最后，將其余28個(gè)樣本用于測(cè)試建立的模型。

由于支持向量機(jī)(SVM)在分類方法中可以很好地解決小樣本、非線性、高維數(shù)的問題(趙洪波等， 2003)，因此，選用該方法對(duì)碎屑斜坡的穩(wěn)定性進(jìn)行判別分類。根據(jù)支持向量機(jī)理論(楊帆等， 2017)，碎屑斜坡安全系數(shù)預(yù)測(cè)判別模型的建立則為尋求SVM模型函數(shù)的成立：

(4)

式中，αi為少數(shù)訓(xùn)練樣本對(duì)應(yīng)的非零拉格朗日因子，yi為碎屑斜坡穩(wěn)定性的指標(biāo)，n為支持向量，是學(xué)習(xí)樣本集中樣本容量的總數(shù)，K(x，xi)為滿足mercer條件的核函數(shù)。支持向量回歸的目標(biāo)函數(shù)為：

(5)

式中，ξi為松弛變量，ξi≥0，C為懲罰因子。由于核函數(shù)和懲罰因子對(duì)預(yù)測(cè)結(jié)果的影響很大，因此，通過對(duì)各種核函數(shù)的測(cè)試，最終得出碎屑斜坡安全系數(shù)預(yù)測(cè)判別模型中核函數(shù)為多項(xiàng)式核函數(shù)K(x，y)=(xy+1)r(r=1， 2，…，m)，且m=2，懲罰因子C=1.6。然后求解出αi，yi，b，即可得到支持向量機(jī)碎屑斜坡安全系數(shù)判別模型。

把剩余的28個(gè)碎屑斜坡對(duì)判別模型進(jìn)行回判，用“1”表示穩(wěn)定，“0”表示不穩(wěn)定，“*”表示結(jié)論錯(cuò)誤。根據(jù)回判結(jié)果與實(shí)際值和BP人工神經(jīng)網(wǎng)絡(luò)(BPNN)判別結(jié)論對(duì)比，詳見表 7，得出SVM模型的判別準(zhǔn)確率達(dá)到89.29%，而BPNN的準(zhǔn)確率為85.71%，而且BPNN把7#不穩(wěn)定斜坡判別為穩(wěn)定斜坡，說明了SVM模型判別比BPNN判別更準(zhǔn)確可靠。分析SVM對(duì)6#、18#和24#3個(gè)碎屑斜坡判定錯(cuò)誤原因，均因安全系數(shù)略偏低，而為了工程安全考慮，所建模型計(jì)算偏于保守，因此，將其判別為不穩(wěn)定斜坡，也表明了建立的模型基本合理可靠。

表 7 支持向量機(jī)穩(wěn)定性評(píng)判模型檢驗(yàn)Table 7 Model test of stability evaluation of support vector machine

為了進(jìn)一步驗(yàn)證模型的推廣性和適用性，利用2017年補(bǔ)充采集國道G318拉薩至樟木段的16個(gè)碎屑斜坡對(duì)模型進(jìn)行檢驗(yàn)，具體預(yù)測(cè)結(jié)果與實(shí)際結(jié)果和BPNN判別結(jié)論對(duì)比(表 8)。SVM判別準(zhǔn)確率為93.75%，而BPNN判別準(zhǔn)確率為87.5%，仍然出現(xiàn)將1#不穩(wěn)定斜坡判別為穩(wěn)定斜坡，推廣應(yīng)用也說明了SVM模型判別比BPNN判別更準(zhǔn)確可靠。分析SVM模型存在1個(gè)斜坡判斷錯(cuò)誤的原因，同樣是安全系數(shù)略偏低，模型計(jì)算偏于保守造成。因此，SVM模型的評(píng)判和適用性檢驗(yàn)準(zhǔn)確率均超過89%，完全符合工程安全要求，表明模型具有較好的適用性和推廣性。

表 8 支持向量機(jī)穩(wěn)定性評(píng)判模型推廣預(yù)測(cè)Table 8 Generalization prediction of support vector machine stability evaluation model

4 結(jié) 論

基于大量的野外調(diào)查成果，對(duì)碎屑斜坡分類和基本特征進(jìn)行了分析論述，運(yùn)用信息熵和FLAC數(shù)值模擬分析了碎屑斜坡的主要影響因素，采用支持向量機(jī)評(píng)判模型進(jìn)行了穩(wěn)定性評(píng)價(jià)。得出以下幾點(diǎn)結(jié)論：

(1)根據(jù)野外調(diào)查的110個(gè)碎屑斜坡數(shù)據(jù)，結(jié)合斜坡碎屑物質(zhì)成分、顆粒大小及運(yùn)動(dòng)特征等對(duì)碎屑斜坡再細(xì)分為溜砂坡、滾石坡和碎石坡3類。

(2)溜砂坡主要為崩坡積物，規(guī)模不大，但數(shù)量眾多，突發(fā)性較強(qiáng)，而且連續(xù)成群出現(xiàn)，危害大； 滾石坡主要為沖洪積物的剝落或局部滑落，危害較大； 碎石坡主要為殘坡積物，其突發(fā)性極強(qiáng)，能量大、速度快、具有一定規(guī)模，危害很大。

(3)利用信息熵計(jì)算得出碎屑斜坡影響因素從大到小依次為：降雨融雪>平均坡度>碎屑物>平均坡高>植被覆蓋率>工程開挖>地震>風(fēng)化等級(jí)>振動(dòng)>河水侵蝕。

(4)通過數(shù)值模擬，統(tǒng)計(jì)得出碎屑斜坡安全系數(shù)斜率大小為：降雨融雪(0.80)>平均坡度(0.66)>平均坡高(0.52)>植被(0.46)>工程開挖(0.43)，其斜率越大則該因素影響作用越顯著，結(jié)論與野外調(diào)查分析一致。

(5)利用信息熵和FLAC數(shù)值模擬分析的主要因素，建立了基于支持向量機(jī)的碎屑斜坡穩(wěn)定性評(píng)判模型，模型檢驗(yàn)和推廣預(yù)測(cè)的準(zhǔn)確率分別為89.29%和93.75%，模型具有較好的可靠性和適用性。

致 謝非常感謝審稿專家對(duì)文章初稿提出的建設(shè)性意見，感謝同濟(jì)大學(xué)趙程教授提供的修改意見，同時(shí)感謝西藏大學(xué)工學(xué)院2203大學(xué)生創(chuàng)新實(shí)驗(yàn)室同學(xué)們?cè)谝巴庹{(diào)查中提供的幫助。