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

    徑流-滲流作用下松散土坡非均勻顆粒臨界失穩(wěn)分析

    2021-11-29 06:01:56袁星宇符文熹李龍國
    中國農村水利水電 2021年11期
    關鍵詞:土坡坡面安全系數(shù)

    劉 悅,袁星宇,張 帥,符文熹,李龍國

    (1.四川大學水利水電學院,成都610065;2.四川大學水力學與山區(qū)河流保護國家重點實驗室,成都610065)

    0 引 言

    松散土坡淺表物質具有結構松散、膠結程度低、抗干擾能力弱、孔隙度大、透水性強等特性,斜坡顆粒極易在降雨、地震、人類活動等因素影響下發(fā)生失穩(wěn)破壞[1,2]。我國西南地區(qū)松散土坡分布廣泛,降水豐富,發(fā)生滑坡和碎屑流災害嚴重。因此,有必要開展坡面徑流與坡體滲流聯(lián)動作用下松散土坡顆粒沖刷啟動機理研究。

    近年來,松散土坡穩(wěn)定性問題受到國內外學者廣泛關注。Rijn 等[3]基于泥沙顆粒的受力特征,通過綜合考慮橫縱向底坡對泥沙起動的疊加效應,建立了任意坡面上泥沙顆粒處于臨界起動狀態(tài)時的拖曳力公式;謝立全等[4]在進行坡面泥沙受力分析時考慮了滲透力作用,發(fā)現(xiàn)向坡外的滲流將大幅降低岸坡泥沙起動條件;Kumar等[5]考慮到岸坡滲流的影響,提出無量綱的水流功率概念,并將其運用到滲流對坡面泥沙起動過程的描述中;Zamani 等[6]通過顆粒離散元法建立了三維土-基礎-結構體系,對松散體在不同多向振動作用下的響應進行模擬;田海等[7]進行離心模型試驗,對降雨條件下邊坡土體強度的變化進行了分析,并得出結論降雨會降低松散堆積體邊坡穩(wěn)定性;毛雪松等[8]基于Midas/GTS 建立松散堆積體路基上邊坡穩(wěn)定數(shù)值模型,分析了各坡形因素對邊坡穩(wěn)定性的影響;張巖巖[9]采用一維地表徑流運動波方程和二維地下水滲流Richards 方程,將地表徑流和地下水滲流耦合,并結合非飽和土Bishop 有效應力及極限平衡理論,研究了降雨作用下滑坡穩(wěn)定性問題;羅浩等[10]考慮棄渣場坡體結構特征,通過建立邊坡顆粒流模型對降雨工況下邊坡變形失穩(wěn)過程進行了模擬;雷小芹等[11]構建描述非飽和堆積土中細顆粒侵蝕-運移-沉積全過程的滲流潛蝕模型,定量研究了細顆粒運移對邊坡降雨穩(wěn)定性的影響。

    綜合查閱的文獻分析表明,有關坡面徑流和坡體滲流對松散土坡失穩(wěn)的聯(lián)合作用鮮有報道。本文從單顆粒的滑動和滾動角度進行分析,分別采用Navier-Stokes方程[12]和擴展Brinkman-Darcy方程[13,14]描述坡面徑流和坡體滲流,聯(lián)合連續(xù)性方程,構建徑流-滲流作用下的分析模型;將徑流-滲流界面視為應力跳躍邊界,推求出流固界面流速表達式,并結合顆粒的滑動和滾動受力分析,得到考慮徑流-滲流聯(lián)動作用下單顆粒穩(wěn)定安全系數(shù);最后討論了單顆粒滑動、滾動安全系數(shù)對各影響因子的敏感強弱,確定安全系數(shù)敏感參數(shù)。

    1 徑流-滲流作用下非均勻顆?;瑒雍蜐L動受力分析

    1.1 顆粒滑動分析

    松散土坡受外界因素(如降雨)影響,表層顆粒會發(fā)生運動。本文假定土坡顆粒為不可壓縮非均勻顆粒(等效粒徑為di),顆粒只沿土坡傾向方向運動(圖1)。

    如圖1 所示,徑流-滲流作用下單個顆粒受到拖曳力FD、上舉力FL、水下重力W,以及由于顆粒相對暴露度產生的附加質量力FM;FD、FL、W和FM的計算分別見式(1)~式(4)[15,16]。

    圖1 松散土坡淺層顆粒受力情況Fig.1 Force of shallow particles on loose soil slope

    式中:CD為拖曳力阻力系數(shù);ρ為水的密度,kg/m3;αD為與拖曳力相應的面積系數(shù),取值范圍0~1;ub為流固界面流速,m/s。

    式中:CL為上舉力系數(shù);αL為與上舉力相應的面積系數(shù),取值范圍0~1;其余符號同前。

    式中:αW為與水下重力相應的面積系數(shù),取值范圍0~1;γs為泥沙顆粒的容重,N/m3;γw為水的容重,N/m3;其余符號同前。

    式中:dmax、dmin、dm分別為顆粒級配中最大、最小和平均粒徑,m;αM為與附加質量力相應的面積系數(shù),取值范圍為0~1;di為研究顆粒等效粒徑,m;dj為研究對象周圍顆粒的粒徑,m;f(dj)為河床表面混合物以重量計算dj顆粒分布函數(shù);ζ為與顆粒相對暴露度(床面粗化度)有關的系數(shù),顆粒完全暴露時ζ=0,顆粒完全隱蔽時ζ=1。

    當土坡顆粒發(fā)生滑動時,顆粒還會受到周圍其他顆粒的摩阻力Ff,F(xiàn)f計算見式(5)。

    式中:φ為土體內摩擦角;其余符號同前。

    將式(3)、式(4)代入式(5)中,F(xiàn)f可改寫為式(6)。

    在此,定義顆??够€(wěn)定安全系數(shù)SFs為:

    式中:Frf為抗滑力,N;Fsf為滑動力,N。

    從圖1 可見,抗滑力僅為摩擦阻力Ff,即式(6)?;瑒恿Π轮亓ρ豿方向的分量Wx=Wsinθ、附加質量力沿x方向的分量FMx=FMsinθ和拖曳力FD。于是,式(7)可具體寫出如下:

    分析式(8)可知,顆??够€(wěn)定安全系數(shù)主要受到流固界面流速ub、土坡坡度θ、土體內摩擦角φ、顆粒等效粒徑di、顆粒容重γs的影響。

    1.2 顆粒滾動分析

    在此,定義顆???jié)L動穩(wěn)定安全系數(shù)SFr為:

    式中:Frm為抗?jié)L動力矩,N?m;Fsm為滾動力矩,N?m。

    從圖2 可見,當土坡顆粒發(fā)生滾動時,支持力FN=0,設顆粒繞O點發(fā)生滾動,且W、FD、FL、FM的力臂分別為lW、lD、lL、lM。于是,F(xiàn)rm=W lw+FMlM,F(xiàn)sm=FDlD+FLlL,相應地,式(9)可具體寫出如下。

    圖2 顆粒滾動趨勢下受力及力臂情況Fig.2 Force and force arm of rolling particles

    分析式(10)可知,顆???jié)L動穩(wěn)定安全系數(shù)主要受到流固界面流速ub、顆粒等效粒徑di、顆粒容重γs的影響。

    2 徑流-滲流作用土坡流場分析

    2.1 模型描述

    根據2.1 節(jié)、2.2 節(jié)顆粒安全系數(shù)計算公式可知,徑流-滲流聯(lián)合作用下松散土坡非均勻顆粒是否穩(wěn)定與流固界面流速ub密切相關;安全系數(shù)隨著ub的增加而減小。為推求松散土坡流速分布、求解流固界面流速ub,本文將坡面徑流及坡體滲流進行耦合分析,如圖3 所示,土坡坡長為L,松散層厚度為b,孔隙率為n,滲透率為K,下層基巖假設為不透水層;松散土坡坡面徑流、坡體滲流均為二維平面流,且沿y方向的流速為0;坡面徑流水深為h,沿x方向的水流流速為ux;坡體滲流沿x方向的水流流速為vx。

    圖3 松散土坡徑流-滲流耦合分析模型Fig.3 Runoff-Seepage coupling analysis model on loose soil slope

    假設圖3 中水流是不可壓縮的Newton 流體,為充分發(fā)展層流。采用連續(xù)性方程[式(11)]和Navier-Stokes 方程[式(12)]描述松散土坡坡面徑流運動[12,17]。

    式中:ux、uy分別為松散土坡坡面徑流沿x、y方向的流速,m/s;fx為沿x方向的質量力,N/kg;P為沿x方向的壓強,Pa;υ為水的運動黏滯系數(shù),m2/s;t為時間,s;其余符號同前。

    已知水流沿y方向的流速為0,故有?uy/?y=0,將其代入式(11)中可得?ux/?x=0;x方向流速ux在y方向不發(fā)生變化,即?ux/?y=0;水流為恒定流,則?ux/?t=0;沿x方向fx=gsinθ,?P/?x=-ΔP/L。因此,式(12)可寫為:

    式中:μ為水的動力黏滯系數(shù),Pa?m,其中μ=υ ρ;L為土坡長度,m。

    采用連續(xù)性方程[式(14)]和擴展Brinkman-Darcy 方程[式(15)]描述松散土坡坡體滲流運動[13,14,18]。

    式中:vx、vy分別為松散土坡坡體滲流沿x、y方向的流速,m/s;n為土顆??紫堵?;K為土坡滲透率,m2,k=K g/υ,其中k為滲透系數(shù),m/s,g為重力加速度,m/s2。

    同理,式(15)可簡化為:

    2.2 界面流速

    為便于計算,現(xiàn)對流速計算公式進行無量綱化處理。引入無量綱參數(shù)c0、c1、c2、c3、c4、V,令c0=b/D、c1=y/D、c2=1/n、c3=K/D2、c4=1/(c2c3)0.5、V=μ v/(ρ g D2tanθ),其中D=b+h,v為水流流速。則式(13)、式(16)可進一步簡化為:

    式中:Vr為松散土坡坡面徑流的無量綱流速;Vs為松散土坡坡體滲流的無量綱流速。

    求解式(17)、式(18),得無量綱流速表達式為:

    式中:a1、a2、b1、b2為待求系數(shù)。

    將無量綱參數(shù)V=μ v/(ρ g D2tanθ)代入式(19)和式(20)中,即可反求出有量綱流速:

    松散土坡坡面徑流ux和坡體滲流vx滿足以下邊界條件[19,20]:

    (1)在松散土坡坡面徑流上表面(y=b+h)處,流速ux達到最大,即滿足?ux/?y=0;

    (2)在土坡表面(y=b)處,滿足速度連續(xù)且剪應力跳躍的邊界條件,即ux=vx、c2?vx/?c1-?ux/?c1=β vx/c30.5,其中β為應力跳躍系數(shù),與多孔介質特性、界面位置密切相關,其取值范圍為-1~1.5;

    (3)在土坡底部(y=0)處,固體顆粒與表面流體不存在相對位移,即滿足ux=0。

    將上述邊界條件代入式(21)、式(22)中,解得待定系數(shù)a1、a2、b1、b2,見式(23)~(26)。

    式(24)~(26)中d1、d2、d3、d4、d5、d6分別見式(27)~(32)。

    因此,流固界面流速ub(y=b)為:

    分析式(33)可知,流固界面流速主要受到松散層厚度b、坡面徑流水深h、土坡坡度θ、應力跳躍系數(shù)β、松散土坡孔隙率n及滲透率K的影響。

    3 顆粒穩(wěn)定性影響因子分析

    根據式(8)、式(10)和式(33),松散土坡非均勻顆粒抗滑穩(wěn)定安全系數(shù)SFs的影響因素有松散層厚度b、坡面徑流水深h、土坡坡度θ、應力跳躍系數(shù)β、土體內摩擦角φ、顆粒等效粒徑di、顆粒容重γs、松散土坡孔隙率n、滲透率K;抗?jié)L動穩(wěn)定安全系數(shù)SFr的影響因素有b、h、θ、β、n、K、di及γs。

    為定量了解上述因素對松散土坡非均勻顆粒穩(wěn)定性的影響程度,對各因子進行敏感性分析。敏感系數(shù)Si計算見式(34),敏感系數(shù)越大,安全系數(shù)對該影響因子越敏感[21,22];敏感分析相關參數(shù)基準值如表1所示,其中拖曳力阻力系數(shù)CD、上舉力系數(shù)CL、顆粒所受各力的面積系數(shù)α及力臂l取值參考文獻[23]、[24]。

    表1 相關計算參數(shù)基準值[23,24]Tab.1 The base value of related calculation parameters

    式中:Si為安全系數(shù)影響因子xi的敏感系數(shù);|ΔSFi/SF0|為安全系數(shù)SFi的相對變化率;SF0為各參數(shù)基準值計算出的安全系數(shù)基準值;|Δxi/x0|為影響因子xi的相對變化率;x0為各參數(shù)基準值。

    3.1 坡面徑流水深h敏感性

    令坡面徑流水深h從0.6 mm 增加到1.8 mm,其余參數(shù)取值如表1 所示。根據不同的h值,求得與h對應的SFs、SFr值,徑流水深h與安全系數(shù)SF關系如圖4所示。

    圖4 安全系數(shù)與敏感系數(shù)隨徑流深變化Fig.4 Variations of safety factors and sensitivity factors with runoff depth of slope surface

    從圖4 可以看出,當h從0.6 mm 增加到1.8 mm 時,安全系數(shù)SFs從1.36 下降到0.47,SFr從4.60 下降到0.88,SFs、SFr隨著h的增加而降低;敏感系數(shù)Shs從1.13 降低到0.46,Shr從3.75 降低到0.39,Shs、Shr隨著h的增加而降低。當h一定時,顆粒更易發(fā)生滑動失穩(wěn)。

    取h=1.0 mm 時的敏感系數(shù)為影響因子h的敏感系數(shù),即Shs=0.90、Shr=1.51,明顯SFr相較于SFs對h更敏感。

    3.2 坡度θ敏感性

    令土坡坡度θ從10°增加到20°,其余參數(shù)取值如表1 所示,計算得與θ對應的SFs、SFr值。土坡坡度θ與安全系數(shù)SF關系如圖5所示。

    從圖5 可以看出,當θ從10°增加到20°時,安全系數(shù)SFs從1.85下降到0.56,SFr從5.11下降到1.26,SFs、SFr隨著θ的增加而降低;敏感系數(shù)Sθs從3.94降低到0.96,Sθr從5.93降低到0.92,Sθs、Sθr隨著θ的增加而降低。當θ一定時,顆粒更易發(fā)生滑動失穩(wěn)。

    取θ=15°時的敏感系數(shù)為影響因子θ的敏感系數(shù),即Sθs=1.77、Sθr=2.04。

    3.3 土體內摩擦角φ敏感性

    令顆粒內摩擦角φ從25°增加到40°,其余參數(shù)取值如表1所示,計算得與φ對應的SFs、SFr值。內摩擦角φ與安全系數(shù)SF關系如圖6所示。

    圖6 安全系數(shù)與敏感系數(shù)隨土體內摩擦角變化Fig.6 Variations of safety factors and sensitivity factors with soil internal friction angle

    從圖6 可以看出,當φ從25°增加到40°時,安全系數(shù)SFs從0.78上升到1.40,SFr恒為2.26,SFs隨著φ的增加而增大;敏感系數(shù)Sφs從1.11升高到1.52,Sφs隨著φ的增加而增大。當φ一定時,顆粒更易發(fā)生滑動失穩(wěn)。

    取φ=30°時的敏感系數(shù)為影響因子φ的敏感系數(shù),即Sφs=1.21。

    3.4 顆粒容重γs敏感性

    令顆粒容重γs從16 kN/m3增加到30 kN/m3,其余參數(shù)取值如表1 所示,計算得與γs對應的SFs、SFr值。顆粒容重γs與安全系數(shù)SF關系如圖7所示。

    圖7 安全系數(shù)與敏感系數(shù)隨顆粒容重變化Fig.7 Variations of safety factors and sensitivity factors with volumetric weight of particles

    從圖7可以看出,安全系數(shù)與容重基本呈線性關系,當γs從16 kN/m3增加到30 kN/m3時,安全系數(shù)SFs從0.48 上升到1.13,SFr從0.90 上升到3.01,SFs、SFr隨著γs的增加而增大;敏感系數(shù)Sγs從1.88降低到0.76,隨著γs的增加而降低,Sγr恒為1.67。當γs一定時,顆粒更易發(fā)生滑動失穩(wěn)。

    土壤顆粒容重的變化,會同時引起抗滑力和滑動力的變化。根據文獻[21]、[25]可知,邊坡抗滑穩(wěn)定安全系數(shù)一般隨著巖土體容重的增加而減小,與上述得出的結論相悖。造成這一結果的主要原因是,松散狀態(tài)下顆粒間的黏聚力c可忽略,即c=0,重力增加抗滑力增加的比例大于滑動力增加的比例,所以得出了抗滑安全系數(shù)隨著容重的增加而增大的結論。

    取γs=25 kN/m3時的敏感系數(shù)為影響因子γs的敏感系數(shù),即Sγs=1.00、Sγr=1.67,明顯SFr相較于SFs對γs更敏感。

    3.5 顆粒等效粒徑di敏感性

    令顆粒等效粒徑di從1 mm 增加到25 mm,其余參數(shù)取值如表1所示。根據不同的di值,求得與di對應的SFs、SFr值,等效粒徑di與安全系數(shù)SF關系如圖8所示。

    圖8 安全系數(shù)與敏感系數(shù)隨顆粒等效粒徑變化Fig.8 Variations of safety factors and sensitivity factors with equivalent diameter of particles

    從圖8可以看出,安全系數(shù)與等效粒徑基本呈線性關系,當di從1 mm增加到25 mm時,安全系數(shù)SFs從0.48上升到1.04,SFr從0.92 上升到2.61,SFs、SFr隨著di的增加而增大;敏感系數(shù)Sds從0.72 降低到0.33,隨著di的增加而降低,Sdr恒為0.62。當di一定時,顆粒更易發(fā)生滑動失穩(wěn)。

    取di=20 mm 時的敏感系數(shù)為影響因子di的敏感系數(shù),即Sds=0.37、Sdr=0.62。

    3.6 土坡孔隙率n敏感性

    令土坡孔隙率n從0.25 增加到0.4,其余參數(shù)取值如表1 所示,計算得與n對應的SFs、SFr值。土坡孔隙率n與安全系數(shù)SF關系如圖9所示。

    從圖9 可以看出,當n從0.25 增加到0.40 時,安全系數(shù)SFs從1.11 下降到0.90,SFr從2.94 下降到2.02,SFs、SFr隨著n的增加而降低;敏感系數(shù)Sns從0.63 降低到0.44,Snr從1.36 降低到0.66,Sns、Snr隨著n的增加而降低。當n一定時,顆粒更易發(fā)生滑動失穩(wěn)。

    圖9 安全系數(shù)與敏感系數(shù)隨土坡孔隙率變化Fig.9 Variations of safety factors and sensitivity factors with porosity of loose soil slope

    取n=0.35 時的敏感系數(shù)為影響因子n的敏感系數(shù),即Sns=0.49、Snr=0.81。

    3.7 土坡滲透率K敏感性

    令土坡滲透率K從0.03 mm2增加到0.08 mm2,其余參數(shù)取值如表1 所示,計算得與K對應的SFs、SFr值。土坡滲透率K與安全系數(shù)SF關系如圖10所示。

    圖10 安全系數(shù)與敏感系數(shù)隨土坡滲透率變化Fig.10 Variations of safety factors and sensitivity factors with permeability of loose soil slope

    從圖10 可以看出,當K從0.03 mm2增加到0.08 mm2時,安全系數(shù)SFs從1.15下降到0.48,SFr從3.15下降到0.90,SFs、SFr隨著K的增加而降低;敏感系數(shù)SKs從0.89 降低到0.32,SKr從1.90降低到0.28,SKs、SKr隨著K的增加而降低。當K一定時,顆粒更易發(fā)生滑動失穩(wěn)。

    取K=3.93×10-3mm2時的敏感系數(shù)為影響因子K的敏感系數(shù),即SKs=0.78、SKr=1.39,明顯SFr相較于SFs對K更敏感。

    3.8 應力跳躍系數(shù)β敏感性

    令應力跳躍系數(shù)β從-0.3 增加到0.7,其余參數(shù)取值如表1所示,計算得與β對應的SFs、SFr值。應力跳躍系數(shù)β與安全系數(shù)SF關系如圖11所示。

    圖11 安全系數(shù)與敏感系數(shù)隨應力跳躍系數(shù)變化Fig.11 Variations of safety factors and sensitivity factors with stress jump coefficient

    從圖11 可以看出,當β從-0.3 增加到0.7 時,安全系數(shù)SFs從1.22 下降到0.46,SFr從3.53 下降到0.87,SFs、SFr隨著β的增加而降低;敏感系數(shù)Sβs從0.06 升高到0.09,Sβr從0.15 降低到0.08,Sβs隨著β的增加而增大,Sβr隨著β的增加而降低。當β一定時,顆粒更易發(fā)生滑動失穩(wěn)。

    取β=0.1 時的敏感系數(shù)為影響因子β的敏感系數(shù),即Sβs=0.08、Sβr=0.13。

    經計算,松散層厚度b的敏感系數(shù)小于10-10,松散層厚度b變化對安全系數(shù)SF的影響過小,幾乎可忽略不計,故不再對其進行詳細的敏感性分析;安全系數(shù)影響因子的敏感系數(shù)值如表2所示。

    由表2 可知,Sθs>Sφs>Sγs≥1>Shs>SKs≥0.5>Sns>Sds>Sβs>Sbs、Sθr>Sγr>Shr>SKr≥1>Snr>Sdr≥0.5>Sβr>Sbr。土坡坡度θ、土體內摩擦角φ、顆粒容重γs對顆粒抗滑安全系數(shù)的影響最敏感,為抗滑安全系數(shù)敏感參數(shù);坡面徑流水深h、滲透率K影響次之,為次敏感參數(shù);顆粒等效粒徑di、松散土坡孔隙率n、應力跳躍系數(shù)β、松散層厚度b為不敏感參數(shù)。土坡坡度θ、顆粒容重γs、坡面徑流水深h、滲透率K對顆???jié)L動安全系數(shù)的影響最敏感,為抗?jié)L動安全系數(shù)敏感參數(shù);松散土坡孔隙率n、顆粒等效粒徑di、影響次之,為次敏感參數(shù);應力跳躍系數(shù)β、松散層厚度b為不敏感參數(shù)。

    表2 安全系數(shù)影響因子敏感系數(shù)Tab.2 The sensitivity coefficient of safety factor’s impact factors

    4 結 論

    (1)對松散土坡淺層非均勻顆粒進行受力分析,分別采用Navier-Stokes方程和擴展Brinkman-Darcy方程描述坡面徑流和坡體滲流運動,并結合連續(xù)性方程構建徑流-滲流聯(lián)合作用數(shù)學模型;視徑流-滲流界面為應力跳躍邊界,推求出流固界面流速解析解和考慮徑流-滲流聯(lián)動作用的單顆粒穩(wěn)定安全系數(shù)??芍?,顆粒穩(wěn)定安全系數(shù)主要受到坡面徑流水深、松散層厚度、土坡坡度、應力跳躍系數(shù)、土體內摩擦角、顆粒等效粒徑及容重、松散土坡孔隙率及滲透率的影響。

    (2)通過敏感性分析,討論了單顆粒滑動、滾動安全系數(shù)對各影響因子的敏感強弱。結果表明,土坡坡度對顆粒穩(wěn)定性的影響相對最大,敏感系數(shù)分別為1.77 和2.04;抗滑安全系數(shù)的敏感參數(shù)為土坡坡度、土體內摩擦角和顆粒容重;抗滑安全系數(shù)的敏感參數(shù)為土坡坡度、顆粒容重、坡面徑流水深及土坡滲透率?!?/p>

    猜你喜歡
    土坡坡面安全系數(shù)
    考慮材料性能分散性的航空發(fā)動機結構安全系數(shù)確定方法
    沖積扇油氣管道坡面侵蝕災害因子分析
    超音速流越過彎曲坡面的反問題
    重力式擋土墻抗滑穩(wěn)定性安全系數(shù)的異性分析及經驗安全系數(shù)方法
    閘室樁基處理后水平抗滑穩(wěn)定安全系數(shù)提高值的估算范圍研究
    面板堆石壩墊層施工及坡面防護
    上海SMP公園土坡場
    戶外探險(2015年2期)2015-02-26 23:30:21
    Overview of Urban PM 2.5 Numerical Forecast Models in China
    接近物體感測庫顯著提升安全系數(shù)
    汽車零部件(2014年6期)2014-09-20 06:29:36
    午夜成年电影在线免费观看| 亚洲天堂国产精品一区在线| 大香蕉久久成人网| 日韩av在线大香蕉| 成人亚洲精品一区在线观看| videosex国产| 91av网站免费观看| 久久久久久久午夜电影| 少妇裸体淫交视频免费看高清 | 日韩精品中文字幕看吧| 婷婷精品国产亚洲av| 日本熟妇午夜| 亚洲一区中文字幕在线| 一二三四在线观看免费中文在| 亚洲一区中文字幕在线| 亚洲国产高清在线一区二区三 | 午夜福利免费观看在线| 欧美在线黄色| 亚洲成国产人片在线观看| 18美女黄网站色大片免费观看| 男女视频在线观看网站免费 | 1024手机看黄色片| 好男人在线观看高清免费视频 | 18禁观看日本| 欧美又色又爽又黄视频| 久久国产精品男人的天堂亚洲| 18禁黄网站禁片免费观看直播| 精品国产美女av久久久久小说| 亚洲熟妇熟女久久| 嫩草影院精品99| 免费高清在线观看日韩| 免费看美女性在线毛片视频| 美女高潮到喷水免费观看| 叶爱在线成人免费视频播放| 夜夜爽天天搞| 欧美zozozo另类| 丝袜美腿诱惑在线| 亚洲精华国产精华精| 一个人观看的视频www高清免费观看 | 一本大道久久a久久精品| 50天的宝宝边吃奶边哭怎么回事| 欧美丝袜亚洲另类 | 午夜久久久久精精品| 怎么达到女性高潮| 成年版毛片免费区| 亚洲国产欧美网| 欧美中文综合在线视频| 女警被强在线播放| 亚洲国产中文字幕在线视频| 一级毛片高清免费大全| 激情在线观看视频在线高清| 亚洲av成人不卡在线观看播放网| 在线永久观看黄色视频| 久久人妻福利社区极品人妻图片| 日韩大码丰满熟妇| a级毛片在线看网站| 69av精品久久久久久| 一个人观看的视频www高清免费观看 | 真人一进一出gif抽搐免费| 国产色视频综合| 99热6这里只有精品| 日本一区二区免费在线视频| www国产在线视频色| 人妻久久中文字幕网| 亚洲一区高清亚洲精品| 一级毛片女人18水好多| netflix在线观看网站| 免费无遮挡裸体视频| 在线观看www视频免费| 国内少妇人妻偷人精品xxx网站 | 亚洲人成伊人成综合网2020| 18禁黄网站禁片午夜丰满| 老司机午夜十八禁免费视频| 成年女人毛片免费观看观看9| 日韩欧美国产在线观看| 黄片播放在线免费| 我的亚洲天堂| 欧美日本视频| av电影中文网址| 国产亚洲精品久久久久久毛片| 天天躁夜夜躁狠狠躁躁| 国产不卡一卡二| 亚洲美女黄片视频| 搡老岳熟女国产| 久久香蕉国产精品| 69av精品久久久久久| 日本一区二区免费在线视频| 久久精品国产99精品国产亚洲性色| www.精华液| 91老司机精品| 超碰成人久久| 亚洲成人免费电影在线观看| 免费看a级黄色片| 精品国产美女av久久久久小说| 琪琪午夜伦伦电影理论片6080| 国产精品 国内视频| 三级毛片av免费| 精品熟女少妇八av免费久了| 亚洲精品国产一区二区精华液| 免费观看精品视频网站| 天天躁夜夜躁狠狠躁躁| 国产精品av久久久久免费| 少妇的丰满在线观看| 欧美精品啪啪一区二区三区| 日韩精品青青久久久久久| 亚洲精品色激情综合| 一二三四在线观看免费中文在| 日本三级黄在线观看| 两人在一起打扑克的视频| 国产日本99.免费观看| 黄色片一级片一级黄色片| 美女 人体艺术 gogo| 欧美国产日韩亚洲一区| 欧美久久黑人一区二区| 亚洲av成人一区二区三| 亚洲 欧美一区二区三区| 九色国产91popny在线| 亚洲最大成人中文| 在线观看www视频免费| www日本黄色视频网| 欧美在线一区亚洲| 人妻丰满熟妇av一区二区三区| ponron亚洲| 亚洲精品粉嫩美女一区| 亚洲欧美精品综合久久99| 午夜激情av网站| 国产精品,欧美在线| 好男人在线观看高清免费视频 | 一a级毛片在线观看| 亚洲精品久久国产高清桃花| 91在线观看av| 少妇粗大呻吟视频| 久久精品aⅴ一区二区三区四区| 免费高清在线观看日韩| 在线观看午夜福利视频| 亚洲第一欧美日韩一区二区三区| 成人精品一区二区免费| 日韩精品中文字幕看吧| 老汉色∧v一级毛片| 欧美av亚洲av综合av国产av| 老汉色av国产亚洲站长工具| 日本免费a在线| 黄色视频,在线免费观看| cao死你这个sao货| 国产精华一区二区三区| 十分钟在线观看高清视频www| 一卡2卡三卡四卡精品乱码亚洲| 久久热在线av| 岛国在线观看网站| 国产精品免费视频内射| 悠悠久久av| 一区福利在线观看| 性色av乱码一区二区三区2| 91麻豆av在线| 亚洲成人久久性| 欧美激情高清一区二区三区| 亚洲欧美日韩高清在线视频| 欧美av亚洲av综合av国产av| 欧美成狂野欧美在线观看| 亚洲av成人一区二区三| 在线免费观看的www视频| 人人妻人人澡人人看| 亚洲国产精品合色在线| 国产精品二区激情视频| 精品午夜福利视频在线观看一区| 亚洲成人久久性| 性色av乱码一区二区三区2| 一本精品99久久精品77| av中文乱码字幕在线| cao死你这个sao货| 色综合婷婷激情| 国产高清激情床上av| 制服诱惑二区| 国产乱人伦免费视频| 在线天堂中文资源库| 久久久久久亚洲精品国产蜜桃av| 欧美三级亚洲精品| 国产成人精品无人区| 人妻丰满熟妇av一区二区三区| 无遮挡黄片免费观看| 又大又爽又粗| 亚洲自偷自拍图片 自拍| 国产私拍福利视频在线观看| 亚洲成国产人片在线观看| 亚洲人成网站在线播放欧美日韩| 欧美性猛交╳xxx乱大交人| 99热这里只有精品一区 | 亚洲第一电影网av| 99国产极品粉嫩在线观看| 一个人观看的视频www高清免费观看 | 久久人妻av系列| 三级毛片av免费| 欧洲精品卡2卡3卡4卡5卡区| 韩国av一区二区三区四区| 最好的美女福利视频网| 亚洲无线在线观看| x7x7x7水蜜桃| 日韩大尺度精品在线看网址| 国产视频一区二区在线看| 精品久久久久久,| 女性生殖器流出的白浆| 男人的好看免费观看在线视频 | 免费观看人在逋| 欧美日韩一级在线毛片| 国产aⅴ精品一区二区三区波| 国产黄片美女视频| 成人三级黄色视频| 国内毛片毛片毛片毛片毛片| 午夜老司机福利片| 最近最新中文字幕大全电影3 | 欧美激情 高清一区二区三区| 99久久无色码亚洲精品果冻| 色av中文字幕| 搞女人的毛片| 男人舔奶头视频| 日韩欧美一区视频在线观看| 51午夜福利影视在线观看| 97超级碰碰碰精品色视频在线观看| 淫妇啪啪啪对白视频| 国产精品亚洲一级av第二区| 一进一出好大好爽视频| 男女视频在线观看网站免费 | 国产又黄又爽又无遮挡在线| 精品福利观看| 国产蜜桃级精品一区二区三区| 18禁裸乳无遮挡免费网站照片 | 侵犯人妻中文字幕一二三四区| 精品卡一卡二卡四卡免费| 欧美最黄视频在线播放免费| av免费在线观看网站| 99re在线观看精品视频| 琪琪午夜伦伦电影理论片6080| 这个男人来自地球电影免费观看| 91麻豆av在线| 亚洲专区字幕在线| 老汉色∧v一级毛片| 亚洲国产精品999在线| 精品欧美国产一区二区三| 欧美成人免费av一区二区三区| 欧美乱码精品一区二区三区| av中文乱码字幕在线| av在线天堂中文字幕| 久久久久国内视频| 亚洲国产中文字幕在线视频| 欧美日韩精品网址| 欧美黑人精品巨大| 一个人观看的视频www高清免费观看 | 中文在线观看免费www的网站 | 一夜夜www| 一级毛片精品| 国产精品久久久人人做人人爽| 国产麻豆成人av免费视频| 日韩大尺度精品在线看网址| 成年版毛片免费区| 色在线成人网| 色综合站精品国产| 欧美日韩福利视频一区二区| 精品国产乱子伦一区二区三区| 黄色成人免费大全| 久久精品aⅴ一区二区三区四区| 亚洲av日韩精品久久久久久密| 久久性视频一级片| 成人18禁高潮啪啪吃奶动态图| 亚洲中文字幕日韩| 啦啦啦 在线观看视频| 国产真实乱freesex| 黄色毛片三级朝国网站| 男人舔女人的私密视频| 免费女性裸体啪啪无遮挡网站| 成人国语在线视频| av电影中文网址| 国产精品久久久av美女十八| 欧美成人性av电影在线观看| 亚洲国产精品久久男人天堂| 变态另类成人亚洲欧美熟女| av欧美777| 国产伦人伦偷精品视频| 久久久精品国产亚洲av高清涩受| 91在线观看av| 少妇粗大呻吟视频| 免费搜索国产男女视频| 亚洲精品av麻豆狂野| 亚洲精品国产区一区二| 在线天堂中文资源库| 久99久视频精品免费| 欧美乱码精品一区二区三区| 波多野结衣av一区二区av| 国产高清有码在线观看视频 | 久久精品91无色码中文字幕| 两个人看的免费小视频| 亚洲专区国产一区二区| 免费电影在线观看免费观看| 欧美成狂野欧美在线观看| 免费在线观看日本一区| 国产一区二区三区视频了| 亚洲av成人av| 美女 人体艺术 gogo| 99国产精品一区二区三区| 婷婷六月久久综合丁香| 亚洲国产欧美一区二区综合| 成人三级做爰电影| 亚洲欧美精品综合一区二区三区| 久久婷婷人人爽人人干人人爱| 在线观看午夜福利视频| 欧美久久黑人一区二区| 国产精品久久久久久精品电影 | 国产成年人精品一区二区| 亚洲精品国产一区二区精华液| 久久人人精品亚洲av| 午夜老司机福利片| 波多野结衣高清无吗| 国产真实乱freesex| 欧美三级亚洲精品| 这个男人来自地球电影免费观看| 777久久人妻少妇嫩草av网站| 黄频高清免费视频| 亚洲va日本ⅴa欧美va伊人久久| 黑丝袜美女国产一区| 真人做人爱边吃奶动态| 色综合站精品国产| 久久久久久久久中文| 国产黄a三级三级三级人| 99久久久亚洲精品蜜臀av| 99久久国产精品久久久| 亚洲成人久久性| 国产不卡一卡二| 久久精品aⅴ一区二区三区四区| 18禁国产床啪视频网站| 欧美丝袜亚洲另类 | 亚洲一区二区三区色噜噜| 久久精品国产清高在天天线| 精品电影一区二区在线| 美女高潮喷水抽搐中文字幕| 啦啦啦 在线观看视频| a在线观看视频网站| 国产aⅴ精品一区二区三区波| 男人舔奶头视频| 制服诱惑二区| 国产99白浆流出| 久久天躁狠狠躁夜夜2o2o| 欧美精品啪啪一区二区三区| 日韩欧美免费精品| 亚洲 欧美一区二区三区| 国产91精品成人一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 91字幕亚洲| 高清毛片免费观看视频网站| 成年版毛片免费区| 两个人视频免费观看高清| 国产aⅴ精品一区二区三区波| 国产av一区二区精品久久| 欧美丝袜亚洲另类 | 亚洲成人免费电影在线观看| 99国产精品一区二区蜜桃av| 欧美乱色亚洲激情| 久久中文字幕一级| 中文字幕高清在线视频| 久久久久久大精品| 中文字幕久久专区| 精品乱码久久久久久99久播| 观看免费一级毛片| 免费在线观看完整版高清| 精品久久蜜臀av无| 亚洲第一av免费看| 91麻豆av在线| 美女午夜性视频免费| 最新在线观看一区二区三区| 人人妻,人人澡人人爽秒播| 18禁黄网站禁片午夜丰满| 日本黄色视频三级网站网址| 香蕉国产在线看| 日本成人三级电影网站| 免费看十八禁软件| 天天一区二区日本电影三级| 亚洲av第一区精品v没综合| 香蕉久久夜色| 国产精品av久久久久免费| 在线国产一区二区在线| 淫秽高清视频在线观看| 精品国产国语对白av| 欧美日韩中文字幕国产精品一区二区三区| 脱女人内裤的视频| 美女高潮到喷水免费观看| 免费观看精品视频网站| 亚洲精品国产精品久久久不卡| 免费在线观看日本一区| 国产精品九九99| 婷婷丁香在线五月| 日本黄色视频三级网站网址| 久久中文看片网| 又黄又粗又硬又大视频| 国产亚洲精品第一综合不卡| 国产精品精品国产色婷婷| 最近在线观看免费完整版| 国产麻豆成人av免费视频| 亚洲专区字幕在线| av在线播放免费不卡| 女生性感内裤真人,穿戴方法视频| 后天国语完整版免费观看| 啦啦啦观看免费观看视频高清| 久久久久久久久中文| 亚洲精品在线美女| 欧美日韩亚洲综合一区二区三区_| 人人澡人人妻人| 亚洲真实伦在线观看| 欧美成人一区二区免费高清观看 | www国产在线视频色| 色在线成人网| 老司机午夜福利在线观看视频| 欧美av亚洲av综合av国产av| 国产亚洲欧美精品永久| 久久精品国产亚洲av香蕉五月| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久亚洲av鲁大| 久久久久久亚洲精品国产蜜桃av| 亚洲第一av免费看| 国产av一区在线观看免费| 亚洲,欧美精品.| 黑人欧美特级aaaaaa片| 亚洲中文av在线| 一级毛片高清免费大全| 99国产综合亚洲精品| 中文字幕精品免费在线观看视频| 亚洲七黄色美女视频| 国产免费男女视频| or卡值多少钱| 国产一区二区在线av高清观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲片人在线观看| 成人18禁高潮啪啪吃奶动态图| 99热只有精品国产| 成人一区二区视频在线观看| or卡值多少钱| 日韩精品中文字幕看吧| 1024视频免费在线观看| 999久久久精品免费观看国产| 国产熟女午夜一区二区三区| 亚洲精品一区av在线观看| 欧美黄色淫秽网站| 久久精品aⅴ一区二区三区四区| 一本精品99久久精品77| 长腿黑丝高跟| 9191精品国产免费久久| 国产精品免费视频内射| www日本在线高清视频| 久久香蕉国产精品| 级片在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 久久久精品国产亚洲av高清涩受| 妹子高潮喷水视频| 日日干狠狠操夜夜爽| 麻豆一二三区av精品| 免费电影在线观看免费观看| 精品久久蜜臀av无| 丁香欧美五月| 国产国语露脸激情在线看| 悠悠久久av| 午夜老司机福利片| 免费高清在线观看日韩| 亚洲va日本ⅴa欧美va伊人久久| 一个人观看的视频www高清免费观看 | 熟妇人妻久久中文字幕3abv| 啪啪无遮挡十八禁网站| 国产午夜福利久久久久久| 日韩有码中文字幕| 亚洲自拍偷在线| 一卡2卡三卡四卡精品乱码亚洲| 久久中文看片网| 亚洲av中文字字幕乱码综合 | 日本熟妇午夜| 国产亚洲欧美在线一区二区| 一级毛片女人18水好多| 亚洲精品美女久久久久99蜜臀| 久久久久久久久免费视频了| 香蕉av资源在线| 国产成人啪精品午夜网站| 精品久久久久久,| 亚洲欧美激情综合另类| 午夜福利高清视频| 久久性视频一级片| 日本成人三级电影网站| 最近最新免费中文字幕在线| 欧美性长视频在线观看| 欧美性猛交╳xxx乱大交人| 看片在线看免费视频| 午夜a级毛片| 校园春色视频在线观看| 一区福利在线观看| 亚洲av电影不卡..在线观看| 嫩草影视91久久| 久9热在线精品视频| 欧美成人午夜精品| 中亚洲国语对白在线视频| 美女国产高潮福利片在线看| 一夜夜www| 无限看片的www在线观看| 俄罗斯特黄特色一大片| 搡老熟女国产l中国老女人| 视频区欧美日本亚洲| 亚洲在线自拍视频| 精品乱码久久久久久99久播| 九色国产91popny在线| 香蕉av资源在线| 国产aⅴ精品一区二区三区波| 国产私拍福利视频在线观看| 高清在线国产一区| 国产免费av片在线观看野外av| av片东京热男人的天堂| 韩国精品一区二区三区| 女警被强在线播放| 成人av一区二区三区在线看| 国产精华一区二区三区| 午夜两性在线视频| 18禁裸乳无遮挡免费网站照片 | 他把我摸到了高潮在线观看| 美女扒开内裤让男人捅视频| 亚洲第一av免费看| 99久久综合精品五月天人人| 91成人精品电影| 桃色一区二区三区在线观看| 99re在线观看精品视频| 淫妇啪啪啪对白视频| 听说在线观看完整版免费高清| 色婷婷久久久亚洲欧美| 亚洲精品中文字幕一二三四区| 欧美另类亚洲清纯唯美| 99久久99久久久精品蜜桃| 老汉色av国产亚洲站长工具| 久久热在线av| 久久天堂一区二区三区四区| 国产片内射在线| 在线永久观看黄色视频| 性欧美人与动物交配| 亚洲国产毛片av蜜桃av| 国产亚洲欧美98| 午夜福利18| 很黄的视频免费| 国产真人三级小视频在线观看| 国产精品一区二区免费欧美| 精品国产乱码久久久久久男人| 动漫黄色视频在线观看| 欧美丝袜亚洲另类 | 亚洲中文字幕一区二区三区有码在线看 | 亚洲国产精品999在线| 成年免费大片在线观看| 国产三级在线视频| 国产成人欧美在线观看| 美女午夜性视频免费| 久久青草综合色| 欧美zozozo另类| 每晚都被弄得嗷嗷叫到高潮| 91麻豆av在线| 女生性感内裤真人,穿戴方法视频| 老司机在亚洲福利影院| 精品一区二区三区四区五区乱码| 91麻豆精品激情在线观看国产| 搡老熟女国产l中国老女人| 禁无遮挡网站| 久久国产精品男人的天堂亚洲| 夜夜躁狠狠躁天天躁| 精品免费久久久久久久清纯| 少妇粗大呻吟视频| 欧美另类亚洲清纯唯美| 午夜福利免费观看在线| 麻豆成人av在线观看| 亚洲人成伊人成综合网2020| 欧美午夜高清在线| 韩国av一区二区三区四区| 日韩大码丰满熟妇| 国产三级在线视频| 高潮久久久久久久久久久不卡| 精品免费久久久久久久清纯| 老熟妇仑乱视频hdxx| 国产精品永久免费网站| 国产一区二区三区在线臀色熟女| 精品不卡国产一区二区三区| 国产成人欧美在线观看| 黄色视频不卡| 欧美精品亚洲一区二区| 成年免费大片在线观看| 国产av不卡久久| 国产精品1区2区在线观看.| 99热6这里只有精品| 久久天堂一区二区三区四区| 欧美绝顶高潮抽搐喷水| 1024手机看黄色片| 亚洲国产看品久久| 久久国产亚洲av麻豆专区| 久久天躁狠狠躁夜夜2o2o| 午夜影院日韩av| 精品国内亚洲2022精品成人| 亚洲黑人精品在线| 男人舔奶头视频| 国产精品一区二区免费欧美| 日韩av在线大香蕉| 韩国精品一区二区三区| 两个人看的免费小视频| 男女之事视频高清在线观看| 天天躁夜夜躁狠狠躁躁| 正在播放国产对白刺激| 日韩欧美免费精品| 国产成人欧美| 欧洲精品卡2卡3卡4卡5卡区| 丝袜美腿诱惑在线| 欧美一区二区精品小视频在线| 欧洲精品卡2卡3卡4卡5卡区| 国产精品av久久久久免费| 亚洲精品国产区一区二| 亚洲男人的天堂狠狠| a级毛片a级免费在线| 亚洲 欧美 日韩 在线 免费| 婷婷丁香在线五月| 亚洲狠狠婷婷综合久久图片| 999精品在线视频| 在线播放国产精品三级|