• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    鐵路路塹區(qū)域風(fēng)吹雪防雪柵效果研究

    2022-03-17 08:03:14李鵬翔白明洲邱樹茂李景賢
    關(guān)鍵詞:雪粒雪深吹雪

    李鵬翔,白明洲,2,邱樹茂,李景賢

    (1.北京交通大學(xué) 土木與建筑工程學(xué)院,北京 100044;2.北京市軌道工程重點實驗室(北京交通大學(xué)),北京 100044;3.新疆鐵道勘察設(shè)計院 地質(zhì)路基設(shè)計分處,烏魯木齊 830011)

    風(fēng)吹雪是雪粒在氣流的作用下發(fā)生蠕移、躍移及懸移等運動的一種自然現(xiàn)象[1]。隨著寒區(qū)鐵路網(wǎng)的日益完善,風(fēng)吹雪現(xiàn)象帶來的積雪災(zāi)害對鐵路運營安全的影響成為亟待解決的問題?,F(xiàn)場監(jiān)測是獲取風(fēng)吹雪特征最直接有效的方法,但交通線路延伸長,風(fēng)吹雪災(zāi)害往往出現(xiàn)在高寒邊疆地區(qū),難以長時間全面進行實地的風(fēng)吹雪觀測,且風(fēng)吹雪災(zāi)害具有明顯的區(qū)域特征,而現(xiàn)場監(jiān)測難以控制風(fēng)雪場條件,所得結(jié)果也具有隨機性[2]。

    通過數(shù)值計算是獲取積雪分布信息的有效方法之一,根據(jù)天氣條件、地理特征等多因素總結(jié)各種不易積雪的公路路基斷面型式和不同防護措施,提取不同斷面流場的風(fēng)速變化情況[3],從流場變化的角度分析積雪的變化特征[4];陳鳳箴[5]考慮防雪柵對公路風(fēng)吹雪災(zāi)害的防治作用,總結(jié)了公路風(fēng)吹雪災(zāi)害防治的防雪柵類型,吳鵬等[6]通過改變路堤高度及邊坡坡率等參數(shù),對不同路基結(jié)構(gòu)形式下的防雪柵作用進行了研究,并總結(jié)了防雪柵欄、擋雪墻、導(dǎo)風(fēng)板和防雪林等4種目前用于新疆公路風(fēng)吹雪的防治措施的相關(guān)設(shè)置參數(shù)[7]。

    目前風(fēng)吹雪災(zāi)害的防治一般通過改變風(fēng)雪場流向以達到目的,其中主要包括應(yīng)用防雪柵(墻)、防雪林及防雪走廊等,實踐表明防雪林在風(fēng)吹雪易發(fā)區(qū)域種植條件要求較高且難以維護[8],防雪走廊入口處易產(chǎn)生積雪[9],采用防雪柵來減輕積雪對道路的影響是相對簡便有效的方法[10]。懷俄明公路自設(shè)置防雪柵以來, 柵欄防護區(qū)域事故率減少了75%,每年減少近9 d因積雪封閉道路的時間[11],國內(nèi)近年來已將防雪柵較多地應(yīng)用在濱綏線、圖佳線、兩伊線等風(fēng)雪災(zāi)害易發(fā)區(qū)域,但由于早期工程建設(shè)的忽視,往往是出現(xiàn)災(zāi)害后才考慮采取防治措施,所采用的方法與布置形式也較多基于經(jīng)驗選擇。

    精伊霍鐵路是中國首次大規(guī)模在鐵路工程設(shè)計時考慮了風(fēng)吹雪災(zāi)害的防護,并選用了包括防雪墻和防雪柵在內(nèi)的組合設(shè)施,十余年來發(fā)揮了顯著的防雪作用,但由于兩種防雪設(shè)施距離較近,防雪墻作用效果沒有充分發(fā)揮,同時因為防雪柵欄距離路基結(jié)構(gòu)過近,也未能充分發(fā)揮防雪柵效果[12];在新建克塔鐵路則主要從鐵路選線設(shè)計原則和路基結(jié)構(gòu)形式(包括設(shè)計積雪平臺和對比不同邊坡坡率)方面對風(fēng)吹雪災(zāi)害進行了重點防護,對所設(shè)置的防雪柵效果則需要進一步進行觀察[13]。

    綜上所述,目前交通領(lǐng)域風(fēng)吹雪災(zāi)害防治成果主要以公路工程為代表,但鐵路路基形式與公路路基有較大不同,而且鐵路線路難以像公路可以靈活繞避雪害區(qū)域,鐵路沿線風(fēng)吹雪災(zāi)害呈現(xiàn)“點多面廣危害大”的特點,但目前關(guān)于鐵路路塹結(jié)構(gòu)形式的積雪量以及防雪柵作用下兩側(cè)堆積雪量和對路塹內(nèi)外雪量的影響尚無統(tǒng)一定論。因此,本文以新建鐵路工程背景為建模基礎(chǔ),采用實際所測風(fēng)雪場參數(shù)為計算條件,定量分析風(fēng)吹雪作用下鐵路路塹內(nèi)外和不同形式防雪柵兩側(cè)的積雪分布。

    1 數(shù)值模型的建立

    基于ANSYS中Fluent計算模塊,并采用user-defined-function自定義加載運動模型和邊界條件,建立三維多相流風(fēng)吹雪分析模型。

    1.1 風(fēng)雪場運動模型

    采用歐拉多相流模型來描述雪粒在空氣中的運動,其中空氣的運動受雷諾時均方程(RANS)支配,并采用RNGk-ε湍流模型使RANS方程封閉后求解,其中k為湍流動能,ε為湍流耗散率[14]。在多相流計算中認(rèn)為雪粒為連續(xù)流體,運動方程[15]為

    (1)

    式中:ρs為雪的密度;f為雪的體積分?jǐn)?shù)(ρs與f的乘積表示雪粒體積濃度);Sc為施密特數(shù),取值為1[16];w為雪粒降落速度。

    驅(qū)使雪粒運動的最小摩擦速度稱為“閾值速度u*t”,當(dāng)雪粒表面受到的摩擦速度u*大于自身運動的閾值速度u*t時便產(chǎn)生運動,雪粒表面的摩擦速度u*由其剪切力計算得出,如式(2)所示。雪表面的堆積和侵蝕運動取決于摩擦速度u*[17],當(dāng)u*大于u*t時,雪粒開始運動并發(fā)生雪量減少的侵蝕作用,侵蝕減少雪量為qero;u*小于u*t時雪粒無法運動產(chǎn)生堆積,堆積增加雪量為qero,如式(3)~(4)所示,并根據(jù)式(5)計算某一點處的總雪量變化。

    (2)

    (3)

    (4)

    qtotal=qdep+qero

    (5)

    式中:τ為雪粒受到的剪切力;qdep為雪粒堆積量;qero為雪粒侵蝕量;qtotal表示雪??偭髁?;Aero為雪粒間黏結(jié)強度,取值為7×10-4。

    在上述運動方程和堆積-侵蝕模型中,雪粒密度139 kg/m3來源于線路區(qū)域降雪后連續(xù)3 d實測的平均值;風(fēng)吹雪災(zāi)害現(xiàn)場測得3 m高處的雪粒宏觀起動風(fēng)速為6 m/s,根據(jù)近地層風(fēng)場對數(shù)率剖面模型算得閾值速度u*t為0.15 m/s,雪粒降落速度w為0.1 m/s,雪粒粒徑為200 μm[18]。

    1.2 計算域參數(shù)

    基于現(xiàn)場風(fēng)吹雪災(zāi)害防治試驗段路基形式,并根據(jù)TB 10098—2017《鐵路線路設(shè)計規(guī)范》要求[19],建立路基面寬度為8 m、路基邊坡坡度為1∶1.75、路基兩側(cè)均有寬度為5 m積雪平臺的單線鐵路路基,計算中來流風(fēng)向垂直于線路方向。

    計算域如圖1所示,其中為保證來流風(fēng)速與出流風(fēng)速的充分發(fā)展,并且減小計算域壁面對積雪分布的影響,流場方向長度為所有計算工況中防雪柵到鐵路路塹上風(fēng)側(cè)邊坡頂部的最大距離的10倍,高度為所有計算工況中最高模型高度的10倍,線路方向則盡可能地延長,最終計算域為800 m(流場方向)×600 m(線路方向)×50 m(高度)。

    圖1 計算域整體示意圖

    計算域邊界條件中入流面設(shè)置為速度入口(velocity-inlet),以監(jiān)測數(shù)據(jù)建立近地層風(fēng)場對數(shù)率剖面模型[20],其中3 m高處來流風(fēng)速為6 m/s,地面粗糙度為2×10-4m;出流面設(shè)置為自由流出口(outflow),該邊界不定義任何參數(shù);計算域頂面與兩側(cè)面設(shè)置為對稱邊界,該邊界上的剪應(yīng)力為0,各物理量的梯度為0,不影響區(qū)域內(nèi)計算;地面、防雪柵和路基模型表面設(shè)置為流體不可穿過的wall邊界。

    1.3 求解器參數(shù)

    計算采用非定常流,初始設(shè)置雪顆粒體積分?jǐn)?shù)vs為0.01,與雪量變化相關(guān)的參數(shù)α和輸出結(jié)果r分別用式(6)和式(7)表示。計算步長設(shè)置為0.01 s,總計算時長為1 800 s(30 min)。

    (6)

    (7)

    式中:r0為無防雪柵時計算前后的雪粒體積分?jǐn)?shù)比值,r為有雪防雪柵時計算前后的雪粒體積分?jǐn)?shù)比值;r0-1為無防雪柵作用下雪粒體積分?jǐn)?shù)的變化量;α為雪量變化率,表示防雪柵作用下雪量的變化程度,正值表示雪量減少,負值表示雪量增加;vf為計算結(jié)束時的測點雪粒體積分?jǐn)?shù),數(shù)值越大,表示雪量越多。

    2 模型驗證

    為驗證數(shù)值模擬結(jié)果的合理性,在新建鐵路試驗段設(shè)置氣象站和防雪柵,在降雪后對試驗段進行雪深測量,并與數(shù)值模擬計算結(jié)果進行比較。

    2.1 工程背景

    本段以新建阿勒泰—富蘊鐵路為研究對象,線路全長154.5 km,所經(jīng)區(qū)域地勢平坦開闊,沿線植被主要為低矮雜草,山坡相對平緩,有發(fā)達的丘陵平原。該地區(qū)氣候干燥,區(qū)域內(nèi)歷史最大瞬時風(fēng)速達到22.1 m/s, 最大雪深940 mm,是新疆雪災(zāi)較嚴(yán)重的地區(qū)之一。根據(jù)全線調(diào)查和前期自建氣象站的風(fēng)場監(jiān)測數(shù)據(jù),選擇某來流風(fēng)向垂直于線路,并出現(xiàn)風(fēng)吹雪現(xiàn)象的區(qū)域作為試驗段,試驗段全長400 m,分成兩段長度相等的區(qū)域Ⅰ和區(qū)域Ⅱ,其中區(qū)域Ⅰ線路外預(yù)設(shè)防雪柵,防雪柵孔隙率20%,高度3 m,與路塹距離40 m;區(qū)域Ⅱ無防雪柵。試驗段沿線路塹深度在0.8~1.2 m之間,路基兩側(cè)積雪平臺寬度為5 m。試驗段防雪柵與雪深測量點位如圖2所示。

    2.2 現(xiàn)場雪深分布

    2018年12月16日至17日自然降雪后,2018年12月18日使用精度為1 mm的鋼尺對試驗段內(nèi)積雪分布進行測量。圖3為試驗段內(nèi)雪深測量結(jié)果等高線圖。由圖3可以看出區(qū)域Ⅰ雪粒在柵兩側(cè)形成堆積,柵兩側(cè)雪深最大值為430 mm,在遠離防雪柵區(qū)域雪深逐漸減小;區(qū)域Ⅱ路塹外雪深波動較小,整體分布較為平均。試驗段區(qū)域自然降雪深度為135 mm,區(qū)域Ⅰ路塹外平均雪深為215 mm,區(qū)域Ⅱ路塹外平均雪深為193 mm,即防雪柵作用下路塹外平均雪深較無防雪柵區(qū)域增加了37.93%;路塹內(nèi)雪深與路塹外的雪深呈相反的趨勢。區(qū)域Ⅰ路基頂面平均雪深為218 mm, 區(qū)域Ⅱ平均雪深為231 mm,即防雪柵作用下路基頂面平均雪深較無防雪柵區(qū)域減少約13.61%。

    (a)試驗段防雪柵現(xiàn)場照片

    (b)測量點位布置

    圖3 試驗段實際測量雪深等高線圖

    2.3 數(shù)值模型驗證

    本文以無人機航測獲得的地形圖數(shù)據(jù)為基礎(chǔ),建立地形分辨率為5 m的模型。計算采用的多相流模型、堆積-侵蝕條件與邊界條件均與第1節(jié)保持一致。

    圖4為對試驗段鐵路路塹內(nèi)外進行風(fēng)吹雪數(shù)值計算得到的雪粒體積分?jǐn)?shù)等高線圖。由圖4可以看出,區(qū)域Ⅰ雪粒體積分?jǐn)?shù)平均值為1.032,比區(qū)域Ⅱ高69.96%。路基頂面雪粒體積分?jǐn)?shù)為1.109,防雪柵對路基頂面的雪量減少率為25.93%,計算得到的雪量變化趨勢與現(xiàn)場實測結(jié)果有較好的一致性。

    圖4 試驗段雪粒體積分?jǐn)?shù)等高線圖

    3 不同防雪柵結(jié)構(gòu)與布置形式

    通常高度過低(小于2 m)的防雪柵難以發(fā)揮阻雪作用,高度過大(大于6 m)的防雪柵在大風(fēng)條件下穩(wěn)定性存在問題并且造價較高[10]。結(jié)合孔隙率為0的防雪墻和存在孔隙結(jié)構(gòu)的防雪柵,并考慮防雪柵與路塹間的不同布置距離,本文對防雪柵不同結(jié)構(gòu)形式和布置距離對路基內(nèi)外積雪分布變化的影響進行了分析,計算變量和工況見表1。

    表1 計算變量和工況

    對各工況的不變量,防雪柵的高度和距離分別選擇5 m和80 m,使柵兩側(cè)盡可能多地積雪,孔隙率選擇中間值50%。試驗同時進行了無防雪柵工況的計算。

    3.1 孔隙率

    圖5為不同孔隙率作用下路塹內(nèi)外沉積雪量。由圖5可見,不同孔隙率防雪柵的柵前與柵后積雪分布形態(tài)明顯不同,也導(dǎo)致了路塹內(nèi)的雪量差異??紫堵蕿?時(無孔隙擋雪墻),柵兩側(cè)雪粒體積分?jǐn)?shù)比值最大分別為1.208,出現(xiàn)在柵前;孔隙率為25%、50%、75%時,柵兩側(cè)雪粒體積分?jǐn)?shù)比值最大分別為1.173、1.168、1.126,且75%孔隙率時該值出現(xiàn)在柵后,即孔隙率從0增至75%的過程中柵兩側(cè)積雪峰值逐漸減小,積雪堆積趨勢從柵前向柵后轉(zhuǎn)移。

    圖5 不同孔隙率作用下路塹內(nèi)外沉積雪量

    利用曲線擬合多項式方程,得到不同工況下路塹外柵前、柵后和柵兩側(cè)影響長度內(nèi)的雪粒體積分?jǐn)?shù)比值的變化情況。防雪柵孔隙率為50%時的計算公式如式(8)~(12)所示,其中防雪柵柵前影響長度按20 m計算,多項式方程的曲線擬合度均在0.9以上。

    rf50=28.316+7.65×10-2x+7.15×10-3x2+

    2.228×10-5x3

    (8)

    rb50=0.99-5.28×10-3x-1.32×10-4x2-

    9.935×10-7x3

    (9)

    (10)

    (11)

    (12)

    圖6為不同防雪柵孔隙率作用下不同位置的沉積雪量,其中折線圖對應(yīng)左縱軸,柱狀圖對應(yīng)右縱軸。從圖6中可以看出,柵前雪粒體積分?jǐn)?shù)比值的平均值與孔隙率成反比,說明柵前雪量隨著孔隙率的增加而減少,柵后雪粒體積分?jǐn)?shù)比值的平均值先增大后減小,即隨著孔隙率的增加?xùn)藕笱┝肯仍龆嗪鬁p小。50%孔隙率防雪柵作用下路塹外雪粒體積分?jǐn)?shù)比值的平均值最大,為1.049。與無防雪柵工況相比,孔隙率為0、25%、50%、75%時,路塹外雪粒體積分?jǐn)?shù)比值的平均值分別增加了51.45%、74.55%、92.95%、45.50%,即隨著孔隙率從0增加到75%,路塹外雪量先增大后減小。

    圖6 不同孔隙率作用下不同位置沉積雪量對比

    鐵路路塹內(nèi)積雪變化趨勢與路塹外相反。無防雪柵時路基頂面測點的雪粒體積分?jǐn)?shù)比值為1.096,當(dāng)孔隙率為0、25%、50%、75%時,路基頂面雪粒體積分?jǐn)?shù)比值分別為1.077、1.082、1.073、1.085,雪量減少率分別為19.44%、14.42%、23.52%、10.98%,即防雪柵孔隙率為50%時路基頂面雪量最少,防雪柵對路基頂面的雪量減少作用最明顯。

    3.2 作用距離

    防雪柵與路基距離除了影響防雪柵兩側(cè)流場和積雪分布外,也涉及占地規(guī)劃問題,使得該因素成為防雪柵設(shè)計乃至鐵路工程初期設(shè)計重要的考慮因素之一。

    圖7為防雪柵與路基不同作用距離時沉積雪量。由圖7可以看出,防雪柵與路基距離變化時柵后的積雪形態(tài)和路塹內(nèi)雪量明顯不同。隨著防雪柵與路基距離的減小,靠近路塹區(qū)域的雪粒體積分?jǐn)?shù)比值越來越大,柵兩側(cè)積雪峰值也逐漸增加。當(dāng)防雪柵和路基相距20 m時,路基頂面雪粒體積分?jǐn)?shù)比值沿風(fēng)向逐漸減小,說明此時路基頂面處于柵后的積雪沉積區(qū),路基頂面雪量甚至大于無防雪柵時的雪量。

    圖7 不同防雪柵作用距離時路塹內(nèi)外沉積雪量

    圖8為不同防雪柵與路基間距離作用下不同位置的雪量對比,其中折線圖對應(yīng)左縱軸,柱狀圖對應(yīng)右縱軸。從圖8可以看出,作用距離60 m時,路塹外雪粒體積分?jǐn)?shù)比值的平均值達到最大,同時路基頂面測點的雪粒體積分?jǐn)?shù)比值最小,這說明使用5 m高、孔隙率為50%的防雪柵對路基頂面的最優(yōu)防護距離約在60 m左右;在作用距離從80 m減小到20 m過程中,路基頂面雪量減少率先從23.52%增大到24.69%,后減小到-9.57%,負值說明雪量相對無防雪柵時增加。

    圖8 不同防雪柵作用距離時不同位置雪量對比

    防雪柵與路基距離較遠時防雪柵下風(fēng)側(cè)積雪沉積區(qū)得到充分發(fā)展,路塹外雪量增加,但此時防雪柵與路塹距離過遠,防雪柵對路塹內(nèi)雪量影響較弱;隨著作用距離的減小,防雪柵對路基頂面雪量的影響逐漸增加,而在過近時路基頂面會處于柵后積雪沉積區(qū),使得路基頂面雪量相對無防雪柵時增大,此時雪量減少率為負值。

    3.3 高度

    選定防雪柵距離為80 m,柵欄孔隙率為50%,圖9為不同高度防雪柵作用下的沉積雪量對比。由圖9可以看出,防雪柵高度為2、3、4、5 m時柵兩側(cè)的雪粒體積分?jǐn)?shù)比值最大分別為1.106、1.112、1.120、1.168,隨柵高增大而增大,同時防雪柵兩側(cè)積雪沉積距離也增大。防雪柵高度為2 m時,柵后雪粒沉積長度為15~20 m;防雪柵高度為3、4、5 m時,柵后雪粒沉積長度分別為20~25 m、30~35 m、40~50 m,即柵后雪粒沉積長度約為防雪柵高度的10~12倍。

    圖9 不同高度時路塹內(nèi)外沉積雪量

    表2為用相同的多項式擬合方法得到的不同柵高下雪量對比結(jié)果。由表2可以看出5 m高時鐵路路塹外雪粒體積分?jǐn)?shù)比值最大達到1.047,與無防雪柵工況相比,2、3、4、5 m高度下鐵路路塹外雪粒體積分?jǐn)?shù)比值分別增加了49.92%、70.90%、87.53%和92.95%。路塹內(nèi)積雪變化趨勢與路塹外積雪變化趨勢相反。柵高為2、3、4、5 m時,路基頂面的雪量減少率分別為10.42%、13.54%、19.79%和23.52%,即防雪柵高度為5 m時,路基頂面雪量最少,防雪柵對路基頂面的雪量減少作用最明顯。

    表2 不同高度時各位置處雪量對比

    4 雪量影響因素正交試驗

    4.1 正交試驗設(shè)計

    根據(jù)上節(jié)計算采用的防雪柵孔隙率、高度和布置距離,建立L16(45)的正交試驗。計算采用的多相流模型、堆積-侵蝕條件與邊界均與上節(jié)保持一致。計算工況與結(jié)果見表3,路塹外柵欄的影響距離內(nèi)為區(qū)域1,其中柵前取值為20 m,柵后取值為防雪柵到路塹的距離;路基頂面為區(qū)域2。

    (13)

    4.2 路塹外柵欄影響長度內(nèi)雪量

    表4為路塹外柵欄影響長度內(nèi)(區(qū)域1)雪粒體積分?jǐn)?shù)比值的極差分析,由表4算得防雪柵孔隙率、高度和作用距離的極差R值分別為3.81×10-3、1.743×10-2和3.047×10-3,即防雪柵高度是影響路塹外積雪沉降量的主要因素,防雪柵的孔隙率和作用距離主要影響防雪柵兩側(cè)積雪的分布形態(tài)和堆積位置,所以對路塹外整體的雪量影響較小。

    4.3 路基頂面雪量

    表5為路基頂面(區(qū)域2)雪粒體積分?jǐn)?shù)比值的極差分析,由表5算得防雪柵孔隙率、高度和作用距離的極差R值分別為7.55×10-3、7.36×10-3和1.36×10-3,即防雪柵的作用距離是影響路基頂面沉積雪量主要因素。當(dāng)防雪柵與路基距離過近時,防雪柵對路基頂面的雪量減少率為負值,隨著作用距離的增加,雪量減少率逐漸增加,路基頂面雪量減少,因此作用距離會對路基頂面的雪量產(chǎn)生較明顯的影響;柵高和孔隙率主要影響路塹外積雪分布和沉積雪量,并間接改變了路基頂面的雪量。

    5 結(jié) 論

    1)鐵路路塹外設(shè)置的防雪柵可以改變流場分布并影響雪粒運動,使雪粒沉積在防雪柵兩側(cè);當(dāng)路塹位于防雪柵下風(fēng)側(cè)時,路塹外沉積雪量增多,路塹內(nèi)雪量較無防雪柵時減少,可達到風(fēng)吹雪災(zāi)害防治的目的。

    2)防雪柵孔隙率主要影響柵兩側(cè)積雪分布形態(tài),隨著孔隙率的增加,柵兩側(cè)的沉積雪量從小范圍的劇烈變化向大范圍的平緩分布轉(zhuǎn)變,孔隙率為50%時柵兩側(cè)的沉積雪量最大。

    3)柵欄高度增加時,防雪柵兩側(cè)最大沉積雪量和防雪柵影響范圍均增加;當(dāng)孔隙率為50%時,防雪柵下風(fēng)側(cè)積雪沉積區(qū)長度約為防雪柵高度的10~12倍。

    4)影響路塹外沉積雪量的主要因素是防雪柵結(jié)構(gòu)形式,影響路塹內(nèi)沉積雪量的主要因素是防雪柵作用距離。16組正交試驗中,防雪柵所致的路塹外雪量增加17.14% ~ 79.80%,對路基頂面的雪量最大減少率為26.17%。

    猜你喜歡
    雪粒雪深吹雪
    我愛下雪
    天山北坡公路風(fēng)吹雪災(zāi)害分布特征及形成機理分析
    一種基于頻率與相位算法的雷達雪深探測方法
    大雪滿弓刀
    南風(fēng)(2018年5期)2018-07-14 02:23:40
    高原冬季雪深與重慶夏季降水的年際關(guān)系研究
    室內(nèi)聽雪〔組詩選二〕
    中國詩歌(2017年5期)2017-11-25 02:52:17
    風(fēng)吹雪的日子
    中國公路(2017年6期)2017-07-25 09:13:56
    室內(nèi)聽雪
    詩潮(2017年2期)2017-03-16 12:28:56
    鐵路防災(zāi)雪深圖像采集的設(shè)計和實現(xiàn)
    偏斜風(fēng)對地物周邊空間吹雪及積雪形成的影響
    国产欧美日韩一区二区三区在线 | 免费av不卡在线播放| 在线观看免费高清a一片| 日韩不卡一区二区三区视频在线| av.在线天堂| 久久av网站| 日韩在线高清观看一区二区三区| 少妇精品久久久久久久| 国产黄色视频一区二区在线观看| 久久久久久久精品精品| av免费观看日本| 欧美日韩视频精品一区| 日本av免费视频播放| 制服丝袜香蕉在线| 免费看光身美女| 成年女人在线观看亚洲视频| 2021少妇久久久久久久久久久| 中国三级夫妇交换| 最黄视频免费看| 午夜日本视频在线| 欧美人与善性xxx| 国产精品99久久久久久久久| 国产成人午夜福利电影在线观看| 91精品国产九色| 日韩一区二区三区影片| 免费在线观看成人毛片| 亚洲va在线va天堂va国产| 国产精品免费大片| 搡女人真爽免费视频火全软件| 日本-黄色视频高清免费观看| 国产成人精品福利久久| 亚洲成人手机| 最新中文字幕久久久久| 自拍偷自拍亚洲精品老妇| 伊人亚洲综合成人网| 午夜福利网站1000一区二区三区| 国产有黄有色有爽视频| 日日啪夜夜撸| 男的添女的下面高潮视频| 日韩视频在线欧美| 精品久久久噜噜| 在线观看www视频免费| av女优亚洲男人天堂| 女性被躁到高潮视频| 日韩强制内射视频| 在线天堂最新版资源| av在线老鸭窝| 麻豆乱淫一区二区| 国产美女午夜福利| 久久久久久久精品精品| av国产精品久久久久影院| 狂野欧美激情性xxxx在线观看| 日韩伦理黄色片| 一级毛片黄色毛片免费观看视频| 高清视频免费观看一区二区| 午夜视频国产福利| 久久这里有精品视频免费| 免费黄色在线免费观看| 日韩成人av中文字幕在线观看| 哪个播放器可以免费观看大片| 亚洲精品一二三| 亚洲精品456在线播放app| 99九九在线精品视频 | 国产一区二区三区av在线| 欧美精品人与动牲交sv欧美| 伊人久久精品亚洲午夜| 精品酒店卫生间| 成年美女黄网站色视频大全免费 | 91久久精品国产一区二区成人| 观看美女的网站| 日日撸夜夜添| 午夜视频国产福利| 在线观看免费视频网站a站| 亚洲图色成人| 国产一区有黄有色的免费视频| 久久99热6这里只有精品| 亚洲美女搞黄在线观看| 一级爰片在线观看| 久久久久久久久久久丰满| 精品亚洲乱码少妇综合久久| 精品人妻一区二区三区麻豆| 久久久久久人妻| 在线亚洲精品国产二区图片欧美 | 中国美白少妇内射xxxbb| av在线app专区| 乱人伦中国视频| 中文字幕免费在线视频6| 欧美日韩av久久| 99久久精品国产国产毛片| 亚洲国产精品专区欧美| 成人亚洲欧美一区二区av| 少妇精品久久久久久久| 97在线人人人人妻| 亚洲欧美日韩另类电影网站| 国产美女午夜福利| 丰满少妇做爰视频| 黑人猛操日本美女一级片| 97在线视频观看| 亚洲美女黄色视频免费看| 菩萨蛮人人尽说江南好唐韦庄| av免费观看日本| 久久99蜜桃精品久久| 街头女战士在线观看网站| 最新的欧美精品一区二区| av在线观看视频网站免费| 日本午夜av视频| 日韩伦理黄色片| 日本黄大片高清| 日韩av不卡免费在线播放| 22中文网久久字幕| 亚洲av日韩在线播放| 91午夜精品亚洲一区二区三区| 天天躁夜夜躁狠狠久久av| 亚洲国产精品一区二区三区在线| 国产av码专区亚洲av| 边亲边吃奶的免费视频| av网站免费在线观看视频| 成人亚洲欧美一区二区av| 精品酒店卫生间| 久热这里只有精品99| 蜜臀久久99精品久久宅男| 两个人的视频大全免费| 欧美区成人在线视频| 国产精品不卡视频一区二区| 免费黄频网站在线观看国产| 久久精品国产鲁丝片午夜精品| 国产一区有黄有色的免费视频| 国产伦理片在线播放av一区| 精品酒店卫生间| 婷婷色麻豆天堂久久| 亚洲av男天堂| 亚州av有码| 狂野欧美激情性xxxx在线观看| 美女福利国产在线| 三级国产精品片| 春色校园在线视频观看| 人妻少妇偷人精品九色| 亚洲国产色片| 国产一区二区三区av在线| 国产永久视频网站| 亚洲精品自拍成人| 草草在线视频免费看| 国产成人精品婷婷| 亚洲内射少妇av| 亚洲精品国产成人久久av| 69精品国产乱码久久久| 亚洲成色77777| 国产成人一区二区在线| 中文字幕人妻丝袜制服| 一区二区三区精品91| 国产精品偷伦视频观看了| 女性生殖器流出的白浆| 能在线免费看毛片的网站| 日韩成人av中文字幕在线观看| 777米奇影视久久| 亚洲久久久国产精品| av在线app专区| av不卡在线播放| 我要看黄色一级片免费的| 亚洲欧美日韩卡通动漫| 五月天丁香电影| 在线观看免费高清a一片| 一二三四中文在线观看免费高清| 一级毛片电影观看| 亚洲欧美成人综合另类久久久| 色吧在线观看| 国产视频首页在线观看| 亚洲av电影在线观看一区二区三区| 国产成人91sexporn| 黑人高潮一二区| 草草在线视频免费看| 少妇丰满av| 亚洲美女黄色视频免费看| 成人国产麻豆网| 精品国产一区二区久久| 亚洲伊人久久精品综合| √禁漫天堂资源中文www| 亚洲精品国产成人久久av| 日韩中文字幕视频在线看片| 国产男女内射视频| 日本猛色少妇xxxxx猛交久久| 国产男人的电影天堂91| 少妇裸体淫交视频免费看高清| 综合色丁香网| 亚洲,一卡二卡三卡| a级毛片免费高清观看在线播放| 国产视频首页在线观看| 王馨瑶露胸无遮挡在线观看| 少妇人妻久久综合中文| 美女国产视频在线观看| 国产黄色免费在线视频| 亚洲av欧美aⅴ国产| 久久这里有精品视频免费| 亚洲精品456在线播放app| 青春草国产在线视频| av国产久精品久网站免费入址| 精品午夜福利在线看| 99热国产这里只有精品6| 成人毛片60女人毛片免费| 你懂的网址亚洲精品在线观看| 国产永久视频网站| 免费播放大片免费观看视频在线观看| 亚洲精品乱久久久久久| 成人18禁高潮啪啪吃奶动态图 | 22中文网久久字幕| 日韩,欧美,国产一区二区三区| 国产欧美亚洲国产| av在线app专区| 中文字幕亚洲精品专区| 欧美丝袜亚洲另类| 女性被躁到高潮视频| 欧美 亚洲 国产 日韩一| 精品国产一区二区久久| 黄片无遮挡物在线观看| 中文在线观看免费www的网站| 日韩亚洲欧美综合| 免费大片18禁| 日日撸夜夜添| 成人黄色视频免费在线看| 视频区图区小说| 日韩不卡一区二区三区视频在线| 免费观看a级毛片全部| 日韩人妻高清精品专区| 国产精品成人在线| 美女大奶头黄色视频| 三级经典国产精品| 国产熟女欧美一区二区| 永久网站在线| 色5月婷婷丁香| 日韩精品免费视频一区二区三区 | 欧美日本中文国产一区发布| 中文字幕制服av| 中文字幕免费在线视频6| a级毛色黄片| 能在线免费看毛片的网站| a 毛片基地| 国产视频内射| 国产在线一区二区三区精| 亚洲人成网站在线观看播放| 看十八女毛片水多多多| 日韩av不卡免费在线播放| 久久久久久久久大av| 丰满迷人的少妇在线观看| 国产精品女同一区二区软件| 卡戴珊不雅视频在线播放| 老司机亚洲免费影院| 亚洲人成网站在线播| 一个人免费看片子| 日韩欧美一区视频在线观看 | 99热这里只有是精品50| 亚洲精品成人av观看孕妇| 51国产日韩欧美| 国产乱来视频区| 亚洲国产最新在线播放| 99久久综合免费| 日韩成人av中文字幕在线观看| 日韩成人av中文字幕在线观看| 观看美女的网站| 观看美女的网站| 激情五月婷婷亚洲| 日韩三级伦理在线观看| 亚洲成人手机| 成人漫画全彩无遮挡| 亚洲精品色激情综合| 国产黄频视频在线观看| 高清黄色对白视频在线免费看 | 少妇 在线观看| 男男h啪啪无遮挡| 久久97久久精品| 亚洲欧美日韩卡通动漫| 搡老乐熟女国产| 国产高清有码在线观看视频| 久久久久人妻精品一区果冻| 有码 亚洲区| 亚洲一级一片aⅴ在线观看| 一区在线观看完整版| 日本午夜av视频| 人人妻人人澡人人看| 三级国产精品欧美在线观看| 亚洲中文av在线| 你懂的网址亚洲精品在线观看| 麻豆成人av视频| 中文资源天堂在线| 中文资源天堂在线| 色网站视频免费| 日本免费在线观看一区| 国产精品.久久久| 国产精品久久久久久久电影| 日韩熟女老妇一区二区性免费视频| 色网站视频免费| 亚洲激情五月婷婷啪啪| 久久国产精品大桥未久av | 午夜影院在线不卡| 国产中年淑女户外野战色| 熟妇人妻不卡中文字幕| 一级爰片在线观看| 欧美老熟妇乱子伦牲交| av有码第一页| 十八禁网站网址无遮挡 | 亚洲电影在线观看av| 一区在线观看完整版| 亚洲欧洲国产日韩| 亚洲怡红院男人天堂| 18禁裸乳无遮挡动漫免费视频| 日本av免费视频播放| 又大又黄又爽视频免费| 嘟嘟电影网在线观看| 插阴视频在线观看视频| 好男人视频免费观看在线| 嫩草影院入口| 男人和女人高潮做爰伦理| 天美传媒精品一区二区| 六月丁香七月| 在线亚洲精品国产二区图片欧美 | 国产成人一区二区在线| 国产成人91sexporn| 97精品久久久久久久久久精品| 久久久精品免费免费高清| 久久精品熟女亚洲av麻豆精品| 少妇高潮的动态图| 91久久精品电影网| 最黄视频免费看| 欧美日韩视频高清一区二区三区二| 一级黄片播放器| 少妇的逼水好多| av福利片在线观看| 欧美日韩一区二区视频在线观看视频在线| 久久久久久久久久人人人人人人| 亚洲av中文av极速乱| 亚洲精品日韩av片在线观看| 九色成人免费人妻av| 97在线视频观看| 午夜老司机福利剧场| 久久99热6这里只有精品| a级毛色黄片| 婷婷色av中文字幕| 在线亚洲精品国产二区图片欧美 | 国产av一区二区精品久久| 中文字幕免费在线视频6| 老熟女久久久| 蜜桃在线观看..| 中国美白少妇内射xxxbb| 十八禁网站网址无遮挡 | 晚上一个人看的免费电影| 久久精品国产a三级三级三级| 久久午夜福利片| 热re99久久国产66热| 精品人妻熟女av久视频| 日日爽夜夜爽网站| 久久人人爽人人爽人人片va| 最近中文字幕高清免费大全6| 欧美精品亚洲一区二区| 在线观看www视频免费| av在线app专区| 日日摸夜夜添夜夜添av毛片| 中文字幕人妻丝袜制服| 亚洲综合色惰| 91久久精品国产一区二区三区| 一本一本综合久久| 韩国高清视频一区二区三区| av线在线观看网站| 久久亚洲国产成人精品v| 亚洲电影在线观看av| 国产成人免费观看mmmm| 99久久综合免费| 一级av片app| 人妻一区二区av| 亚洲精品日韩av片在线观看| 久久ye,这里只有精品| 国内精品宾馆在线| 青春草视频在线免费观看| 少妇高潮的动态图| 亚洲欧美精品专区久久| 一级a做视频免费观看| 日韩亚洲欧美综合| 欧美97在线视频| 熟女电影av网| 一级毛片 在线播放| 高清欧美精品videossex| 丝袜喷水一区| 国产熟女午夜一区二区三区 | a级一级毛片免费在线观看| 日韩制服骚丝袜av| 国产美女午夜福利| 另类亚洲欧美激情| 国产男人的电影天堂91| 国产无遮挡羞羞视频在线观看| 日韩,欧美,国产一区二区三区| xxx大片免费视频| 妹子高潮喷水视频| 日本-黄色视频高清免费观看| 久久人人爽av亚洲精品天堂| 免费人成在线观看视频色| 精品国产乱码久久久久久小说| 韩国av在线不卡| 熟女电影av网| 色视频在线一区二区三区| 大又大粗又爽又黄少妇毛片口| 国产精品成人在线| 一区二区三区免费毛片| 男女无遮挡免费网站观看| 搡女人真爽免费视频火全软件| 人妻系列 视频| 老熟女久久久| 欧美精品一区二区免费开放| 韩国av在线不卡| 永久免费av网站大全| 久久韩国三级中文字幕| 欧美日本中文国产一区发布| 久久久精品94久久精品| 丰满饥渴人妻一区二区三| 97超碰精品成人国产| 日韩中文字幕视频在线看片| 下体分泌物呈黄色| 亚洲国产精品一区二区三区在线| 国产成人aa在线观看| 蜜臀久久99精品久久宅男| 欧美日本中文国产一区发布| 国产黄片视频在线免费观看| 赤兔流量卡办理| 偷拍熟女少妇极品色| 免费大片黄手机在线观看| 如何舔出高潮| 亚洲欧美成人综合另类久久久| 亚洲精品乱码久久久v下载方式| 日日撸夜夜添| 午夜免费观看性视频| 狂野欧美激情性xxxx在线观看| 欧美人与善性xxx| 久久国产亚洲av麻豆专区| 欧美区成人在线视频| 日韩亚洲欧美综合| 欧美亚洲 丝袜 人妻 在线| 如何舔出高潮| 丰满人妻一区二区三区视频av| 久久狼人影院| 在线观看国产h片| 久久久久久久久久成人| 美女内射精品一级片tv| 寂寞人妻少妇视频99o| 久久精品夜色国产| 欧美日韩精品成人综合77777| 久久精品国产鲁丝片午夜精品| 国产一区亚洲一区在线观看| 97超碰精品成人国产| av在线播放精品| 欧美区成人在线视频| 麻豆成人午夜福利视频| 婷婷色麻豆天堂久久| 国产美女午夜福利| 婷婷色综合大香蕉| 99久久人妻综合| 一级毛片 在线播放| 精品一品国产午夜福利视频| 美女xxoo啪啪120秒动态图| 国产一区二区三区综合在线观看 | 日本猛色少妇xxxxx猛交久久| a级片在线免费高清观看视频| 天美传媒精品一区二区| 国产av码专区亚洲av| 中文乱码字字幕精品一区二区三区| 最新的欧美精品一区二区| 特大巨黑吊av在线直播| 又大又黄又爽视频免费| 免费少妇av软件| 中国国产av一级| 日韩熟女老妇一区二区性免费视频| 久久国内精品自在自线图片| 日韩av在线免费看完整版不卡| 香蕉精品网在线| 一级毛片 在线播放| 国产伦精品一区二区三区四那| 久久影院123| 日韩熟女老妇一区二区性免费视频| 2021少妇久久久久久久久久久| freevideosex欧美| 高清毛片免费看| 久久人人爽人人爽人人片va| 在线精品无人区一区二区三| 一本久久精品| 蜜臀久久99精品久久宅男| 免费黄网站久久成人精品| 国产男人的电影天堂91| 亚洲久久久国产精品| 日韩欧美精品免费久久| 男人舔奶头视频| 晚上一个人看的免费电影| 一本—道久久a久久精品蜜桃钙片| 少妇被粗大的猛进出69影院 | 在线观看一区二区三区激情| 午夜久久久在线观看| 亚洲无线观看免费| 亚洲婷婷狠狠爱综合网| 男女啪啪激烈高潮av片| 内地一区二区视频在线| 免费久久久久久久精品成人欧美视频 | 久久99精品国语久久久| 亚洲精品一区蜜桃| 国产亚洲午夜精品一区二区久久| 岛国毛片在线播放| 麻豆成人午夜福利视频| 精品久久久久久久久亚洲| 观看免费一级毛片| 一级二级三级毛片免费看| 三级经典国产精品| 欧美xxxx性猛交bbbb| 久久久久精品久久久久真实原创| 国产淫语在线视频| 日韩成人av中文字幕在线观看| 熟女电影av网| 性色av一级| 亚洲伊人久久精品综合| 成人美女网站在线观看视频| 天美传媒精品一区二区| 国产免费福利视频在线观看| 久久精品久久久久久久性| 精品一区二区免费观看| kizo精华| 插阴视频在线观看视频| 国产精品福利在线免费观看| 丝袜喷水一区| 午夜免费鲁丝| 国产精品一区二区性色av| 国产精品久久久久久久久免| 国产男人的电影天堂91| 久久久久国产网址| a 毛片基地| 国内揄拍国产精品人妻在线| 国产av一区二区精品久久| 精品少妇内射三级| 精品人妻偷拍中文字幕| av线在线观看网站| 久久久久人妻精品一区果冻| 熟女电影av网| 国产精品久久久久久精品古装| av免费在线看不卡| 高清视频免费观看一区二区| 日本欧美国产在线视频| 五月伊人婷婷丁香| 久久婷婷青草| 免费观看无遮挡的男女| 国产女主播在线喷水免费视频网站| 久久久久久久大尺度免费视频| 少妇的逼好多水| 五月开心婷婷网| 下体分泌物呈黄色| 日本wwww免费看| kizo精华| 高清欧美精品videossex| 如日韩欧美国产精品一区二区三区 | 亚洲欧美日韩另类电影网站| 成人黄色视频免费在线看| 国产乱来视频区| av网站免费在线观看视频| 日本vs欧美在线观看视频 | 精华霜和精华液先用哪个| 国产乱人偷精品视频| 2018国产大陆天天弄谢| 国产亚洲av片在线观看秒播厂| 日韩大片免费观看网站| 五月天丁香电影| 只有这里有精品99| 91精品国产九色| 久久久久精品久久久久真实原创| 极品少妇高潮喷水抽搐| 成人国产av品久久久| 亚洲av国产av综合av卡| 韩国av在线不卡| 久久久午夜欧美精品| 另类精品久久| 国产欧美日韩综合在线一区二区 | 人体艺术视频欧美日本| 中文字幕精品免费在线观看视频 | av黄色大香蕉| 精品午夜福利在线看| 欧美变态另类bdsm刘玥| 亚洲国产精品一区三区| 九草在线视频观看| 精品久久久精品久久久| 天天操日日干夜夜撸| 国产日韩欧美视频二区| 亚洲熟女精品中文字幕| 99久久人妻综合| 久久人人爽人人爽人人片va| 亚洲精品一二三| 一二三四中文在线观看免费高清| 国产av码专区亚洲av| 黄色视频在线播放观看不卡| 日本av免费视频播放| 99久久人妻综合| 一级片'在线观看视频| 亚洲欧美成人精品一区二区| 一级毛片我不卡| 日本午夜av视频| 日韩电影二区| 嫩草影院入口| 欧美日韩在线观看h| 精品一区二区三卡| 好男人视频免费观看在线| 91aial.com中文字幕在线观看| 精品久久久噜噜| 大片免费播放器 马上看| 中文欧美无线码| 精品一区二区三卡| 国产免费一区二区三区四区乱码| 91aial.com中文字幕在线观看| 久久婷婷青草| 国语对白做爰xxxⅹ性视频网站| 免费观看a级毛片全部| 少妇被粗大的猛进出69影院 | 精品久久久久久久久av| 人人妻人人澡人人爽人人夜夜| 青春草国产在线视频| 2018国产大陆天天弄谢| 欧美区成人在线视频|