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

    坡面薄層水流優(yōu)勢(shì)流速研究

    2016-11-12 06:40:59陳麗燕雷廷武啜瑞媛
    中國(guó)水土保持科學(xué) 2016年5期
    關(guān)鍵詞:薄層坡面坡度

    陳麗燕,雷廷武?,啜瑞媛

    (1.中國(guó)農(nóng)業(yè)大學(xué),水利與土木工程學(xué)院,100083,北京;2.天津市龍網(wǎng)科技發(fā)展有限公司,300181,天津)

    坡面薄層水流優(yōu)勢(shì)流速研究

    陳麗燕1,雷廷武1?,啜瑞媛2

    (1.中國(guó)農(nóng)業(yè)大學(xué),水利與土木工程學(xué)院,100083,北京;2.天津市龍網(wǎng)科技發(fā)展有限公司,300181,天津)

    降雨形成的徑流是產(chǎn)生坡面土壤侵蝕的主要?jiǎng)恿?lái)源,徑流流速是土壤侵蝕模型的重要參數(shù)之一。為研究電解質(zhì)示蹤法測(cè)量坡面水流流速過(guò)程中電解質(zhì)優(yōu)勢(shì)流速和水流流速的關(guān)系,本研究利用實(shí)驗(yàn)水槽,在坡度4°、8°、12°,流量12、24、48 L/min條件下,于距離電解質(zhì)注入位置0.3、0.6、0.9、1.2、1.5m處放置探針測(cè)量電解質(zhì)傳遞過(guò)程,計(jì)算不同工況下各測(cè)量斷面的電解質(zhì)優(yōu)勢(shì)流速。結(jié)果表明:流量對(duì)電解質(zhì)優(yōu)勢(shì)流速的影響大于坡度對(duì)其影響,電解質(zhì)優(yōu)勢(shì)流速隨距離增加而增大,采用指數(shù)函數(shù)擬合計(jì)算得到的電解質(zhì)優(yōu)勢(shì)流速隨距離的變化過(guò)程,得到穩(wěn)定的電解質(zhì)優(yōu)勢(shì)流速,即水流優(yōu)勢(shì)流速,其范圍在0.241~0.568m/s之間。隨坡度和流量的增大,水流優(yōu)勢(shì)流速均增大。流量對(duì)水流優(yōu)勢(shì)流速增長(zhǎng)的影響大于坡度對(duì)其的影響。不同坡度和流量條件下,水流優(yōu)勢(shì)流速與平均流速基本一致,二者的比值為1.007,水流優(yōu)勢(shì)流速與最大流速的比值為0.774,平均流速與最大流速的比值為0.776,符合坡面薄層水流的流態(tài)。結(jié)果可為研究坡面薄層水流動(dòng)力過(guò)程提供新的計(jì)算方法和參考數(shù)據(jù)。

    薄層水流;電解質(zhì)示蹤法;流量法;優(yōu)勢(shì)流速;平均流速

    降雨形成的坡面水流是地表土壤侵蝕的主要?jiǎng)恿χ?,地表徑流的水力學(xué)特性是研究土壤侵蝕過(guò)程物理模型的基礎(chǔ)[1]。坡面水流不同于明渠水流,其水深極淺,一般為幾厘米甚至幾毫米,稱(chēng)為薄層水流。薄層水流運(yùn)動(dòng)受降雨及地表下墊面狀況影響較大,如降雨強(qiáng)度、降雨歷時(shí)、土壤種類(lèi)或質(zhì)地、前期水分條件、植被密度和類(lèi)型、坡度、坡長(zhǎng)[2]和土壤表面石塊的密度以及埋沒(méi)深度[3]等。坡面薄層水流流速是土壤侵蝕模型中的重要參數(shù)之一,是計(jì)算其他水力學(xué)要素如佛羅德數(shù)、雷諾數(shù)等的重要參數(shù)。

    薄層水流的流態(tài)不同于明渠水流,因此其流速的計(jì)算方法也不能沿用明渠水流的計(jì)算方法。目前用于測(cè)量薄層水流流速的便攜式儀器比較少,一些儀器如聲學(xué)多普勒測(cè)速儀(ADV)[4]、粒子圖像測(cè)速儀[5]、熱膜流速計(jì)[6-7]、光電或電導(dǎo)傳感器[8-10]等,由于價(jià)格昂貴、日常維護(hù)費(fèi)用高等各種條件限制不能廣泛應(yīng)用。流量法只適用于斷面規(guī)則的水流,如實(shí)驗(yàn)室人工模擬水槽中的水流流速。常用的方法主要是示蹤法,示蹤劑通常為染色劑(一般為KMnO4溶液)或電解質(zhì)(如NaCl或KCl)[11-12],示蹤粒子通常為聚苯乙烯泡沫粒子等。染色劑示蹤法測(cè)量誤差較大,因此提出電解質(zhì)示蹤法。其基本原理是通過(guò)儀器測(cè)量水流斷面的電導(dǎo)率,來(lái)確定水流中電解質(zhì)到達(dá)被測(cè)斷面所需要的時(shí)間,從而計(jì)算水流流速,避免染色劑示蹤法目測(cè)帶來(lái)的人為操作誤差。受泥沙含量以及水流流態(tài)的影響,由最大流速計(jì)算平均流速的校正系數(shù)一直不能準(zhǔn)確確定[13-15]。

    在測(cè)得的電導(dǎo)率曲線(xiàn)中,由于水流的導(dǎo)電性以及水流的紊動(dòng),判斷電解質(zhì)溶液最早到達(dá)測(cè)點(diǎn)的時(shí)間有難度,而測(cè)量電導(dǎo)率達(dá)到最大值的時(shí)間比較容易確定;因此室內(nèi)實(shí)驗(yàn)更多采用測(cè)量得到的電導(dǎo)率計(jì)算水流的優(yōu)勢(shì)流速[16-17]。目前對(duì)于經(jīng)驗(yàn)系數(shù)α的研究較多[18-21],不同水流流態(tài)α值不同。國(guó)內(nèi)一般認(rèn)為水流為層流α=0.67,混合水流或過(guò)度流α=0.7,紊流α=0.8[22]。理論上水流速度增大使校正系數(shù)增大,但在夏衛(wèi)生等的實(shí)驗(yàn)中測(cè)量得到不同泥沙含量下水流速度對(duì)校正系數(shù)的影響并不顯著[23-24]。

    本研究采用電解質(zhì)示蹤法測(cè)量電解質(zhì)隨薄層水流運(yùn)動(dòng)不同測(cè)量斷面電導(dǎo)率隨時(shí)間的變化過(guò)程,計(jì)算不同位置的電解質(zhì)優(yōu)勢(shì)流速,用指數(shù)函數(shù)擬合不同工況下電解質(zhì)優(yōu)勢(shì)流速隨距離的變化得到穩(wěn)定的電解質(zhì)優(yōu)勢(shì)流速,即水流的優(yōu)勢(shì)流速。采用流量法計(jì)算坡面薄層水流的平均流速,采用漂浮物法計(jì)算水流最大流速,并分析三種流速之間的相關(guān)關(guān)系。

    1 研究方法

    1.1 坡面薄層水流溶質(zhì)遷移模型

    在較短距離內(nèi),忽略土壤入滲作用和降雨產(chǎn)生的影響,可假定水流流速變化不大[25-26],根據(jù)菲克定律和質(zhì)量守恒定律,當(dāng)上邊界條件假設(shè)為脈沖函數(shù)時(shí),描述溶質(zhì)在一維穩(wěn)態(tài)流中的對(duì)流彌散方程的解析解即為電解質(zhì)在水流中的遷移過(guò)程[23]。

    1.2 優(yōu)勢(shì)流速、最大流速以及平均流速的計(jì)算方法

    1.2.1 電解質(zhì)示蹤法 采用薄層水流流速測(cè)量?jī)x測(cè)量不同工況下水槽中各斷面的電導(dǎo)率隨時(shí)間變化,根據(jù)測(cè)量斷面電導(dǎo)率變化對(duì)應(yīng)的時(shí)間計(jì)算該斷面的流速。測(cè)量斷面電導(dǎo)率隨時(shí)間變化如圖1所示。

    圖1 測(cè)量斷面電導(dǎo)率隨時(shí)間的變化Fig.1 Measured electronic conductivity as a function of time

    圖1中,電解質(zhì)脈沖從加入水流到運(yùn)動(dòng)至測(cè)量斷面所用的時(shí)間為電解質(zhì)以最大流速傳遞通過(guò)給定距離所用的時(shí)間Te,測(cè)量斷面電解質(zhì)濃度達(dá)到最大值所用的時(shí)間為T(mén)p,電解質(zhì)質(zhì)心到達(dá)測(cè)量斷面所用的時(shí)間為T(mén)c。由此可以分別根據(jù)相應(yīng)的距離計(jì)算得到電解質(zhì)的最大流速,優(yōu)勢(shì)流速和質(zhì)心運(yùn)動(dòng)速度,其中質(zhì)心運(yùn)動(dòng)速度經(jīng)公式推導(dǎo)得出與平均流速相等。

    本實(shí)驗(yàn)中最大流速同時(shí)采用漂浮物法測(cè)量,平均流速采用流量法測(cè)量。電解質(zhì)優(yōu)勢(shì)流速的計(jì)算公式如下:

    式中:up為電解質(zhì)優(yōu)勢(shì)流速,m/s;Lp為電解質(zhì)脈沖注入點(diǎn)到測(cè)量斷面的距離,m;Tp為從電解質(zhì)脈沖注入水流到測(cè)量斷面電導(dǎo)率達(dá)到峰值所用的時(shí)間,s。

    1.2.2 最大流速的計(jì)算

    試驗(yàn)采用漂浮物法測(cè)量水流的最大流速。漂浮物法是將可視性及跟隨性較好的粒子加入到流動(dòng)的水體中,讓它跟隨流體流動(dòng),從而測(cè)定流體流動(dòng)情況。示蹤物有液體、懸浮物、漂浮物或沉淀物[27],本試驗(yàn)選擇直徑5 mm的泡沫粒子作為示蹤物,試驗(yàn)中將若干粒子同時(shí)快速加入水流中,記錄泡沫粒子從第2組探針到最后一組探針的時(shí)間,從而計(jì)算得到水流的最大流速u(mài)e。計(jì)算公式為

    式中:ue為最大流速,m/s;Le為第2組探針到最后一組探針的距離,m;Te為泡沫粒子從第2組探針運(yùn)動(dòng)到最后一組探針的時(shí)間,s。

    1.2.3 平均流速的計(jì)算

    土壤侵蝕研究中的水流流速通常用流量法標(biāo)定[27]。

    對(duì)于規(guī)則河床內(nèi)水流,流量適當(dāng)時(shí),可采用流量法測(cè)量計(jì)算水流流速,測(cè)量原理簡(jiǎn)單,計(jì)算公式為

    由于水流深度在厘米甚至毫米級(jí),并且極其不均勻,水流深度的測(cè)量會(huì)產(chǎn)生一定誤差。在實(shí)際坡面徑流測(cè)量中,斷面形狀的不規(guī)則也會(huì)影響最終結(jié)果的準(zhǔn)確性。

    2 試驗(yàn)材料、儀器與方法

    試驗(yàn)主要由薄層水流流速測(cè)量系統(tǒng)和水槽組成,薄層水流流速測(cè)量系統(tǒng)主要由電解質(zhì)脈沖發(fā)生器、感應(yīng)探針、數(shù)據(jù)采集管理器、操作控制計(jì)算機(jī)系統(tǒng)組成[28],如圖2所示。

    試驗(yàn)采用有機(jī)玻璃制作的長(zhǎng)4 m,寬14.7 cm,高50 cm的水槽,底部粘接粒砂紙,模擬土壤下墊面。水槽前端設(shè)有水流注入系統(tǒng),模擬徑流。依據(jù)馬氏瓶原理設(shè)計(jì)而成的恒壓水箱控制水槽上方來(lái)水流量,流量通過(guò)板閥控制。

    在水槽前端固定電解質(zhì)脈沖發(fā)生器,在距離脈沖發(fā)生器下游0.3、0.6、0.9、1.2、1.5m處各固定一組探針,將電解質(zhì)脈沖發(fā)生器和探針?lè)謩e與數(shù)據(jù)采集管理器的相應(yīng)端口連接。

    試驗(yàn)坡度為4°、8°、12°,流量設(shè)置為12、24、48 L/min,試驗(yàn)重復(fù)3次。

    3 結(jié)果與分析

    3.1 電解質(zhì)優(yōu)勢(shì)流速

    在給定的坡度和設(shè)定的流量條件下,利用電解質(zhì)示蹤法測(cè)量得到水流流經(jīng)各斷面的電導(dǎo)率隨時(shí)間變化數(shù)據(jù),根據(jù)式(1)計(jì)算得到各測(cè)量斷面的電解質(zhì)優(yōu)勢(shì)流速,同時(shí)采用流量法測(cè)量并計(jì)算水流的平均流速,如表1所示。

    圖2 薄層水流流速測(cè)量系統(tǒng)示意圖Fig.2 Schematic diagram for themeasurement system of shallow water flow velocity

    表1 平均流速和各斷面的電解質(zhì)優(yōu)勢(shì)流速Tab.1 Mean velocity and electrolyte peak velocity at different section m/s

    為研究坡度、流量和距離對(duì)電解質(zhì)優(yōu)勢(shì)流速的影響,作圖如圖3??梢钥闯?,不同坡度和流量條件下,隨著坡度、流量的增大,電解質(zhì)的優(yōu)勢(shì)流速增大,流量對(duì)流速的影響較流速更為明顯,流速隨距離呈增大趨勢(shì),增大程度逐漸減小,最終趨于穩(wěn)定,即與水流流速相同。

    3.2 穩(wěn)定電解質(zhì)優(yōu)勢(shì)流速及水流優(yōu)勢(shì)流速

    董月群等[28]建議的流速隨距離變化的模型計(jì)算公式為

    式中:u為計(jì)算得到的電解質(zhì)優(yōu)勢(shì)流速,m/s;v為達(dá)到穩(wěn)定的電解質(zhì)優(yōu)勢(shì)流速,即水流優(yōu)勢(shì)流速,m/s;x為探針距離電解質(zhì)脈沖發(fā)生器的距離,m;β為測(cè)量流速隨距離增加的速度,1/m。

    根據(jù)公式(4)將計(jì)算得到的各斷面水流中電解質(zhì)的優(yōu)勢(shì)流速利用Grapher進(jìn)行擬合,各工況下擬合結(jié)果如圖4所示。

    不同工況下,電解質(zhì)的優(yōu)勢(shì)流速隨距離增大而增大,變化規(guī)律符合指數(shù)函數(shù),擬合系數(shù)均大于0.9,擬合效果較好。

    電解質(zhì)注入水流后并不能立刻與水流流速相同,而是存在一個(gè)呈指數(shù)函數(shù)的加速過(guò)程,最終達(dá)到穩(wěn)定后與水流流速相同,因此擬合得到的穩(wěn)定電解質(zhì)優(yōu)勢(shì)流速即可作為水流優(yōu)勢(shì)流速。

    3.3 坡度和流量對(duì)水流優(yōu)勢(shì)流速變化的影響

    各工況下的水流優(yōu)勢(shì)流速以及坡度和流量對(duì)水流優(yōu)勢(shì)流速的影響見(jiàn)表2。

    不同工況下,水流優(yōu)勢(shì)流速隨流量、坡度的增大而增大,變化范圍在0.241~0.568 m/s之間。坡度固定時(shí),流量從12 L/min增大到24 L/min時(shí),水流優(yōu)勢(shì)流速的增大程度較流量從24 L/min增大到48 L/min時(shí)大,最大可達(dá)54.2%;但隨著坡度增大,流量從12 L/min增大到24 L/min時(shí)增大程度遞減,從54.2%減小到29.4%,流量從24 L/min到48 L/min時(shí)增大程度遞增,從16.2%增大到30%。流量固定時(shí),坡度從4°增加到8°時(shí)水流優(yōu)勢(shì)流速增大較明顯,在16.2~29.2%之間,坡度從8°增加到12°時(shí),水流優(yōu)勢(shì)流速的增大程度明顯減小,在2.3%~3.6%之間。流量對(duì)水流優(yōu)勢(shì)流速的影響在16.4%~54.2%之間,坡度對(duì)優(yōu)勢(shì)流速的影響在2.3%~29.4%之間;因此流量對(duì)水流流速的影響大于坡度的影響。

    圖3 不同試驗(yàn)條件下各測(cè)量斷面的電解質(zhì)優(yōu)勢(shì)流速Fig.3 Electrolyte peak velocities at different distance under different experimental conditions

    表2 流量和坡度對(duì)水流優(yōu)勢(shì)流速的影響Tab.2 The impact of flow rate and slope on flow peak velocities

    3.4 優(yōu)勢(shì)流速與平均流速的關(guān)系

    將擬合得到的水流優(yōu)勢(shì)流速與流量法測(cè)量得到的平均流速進(jìn)行比較,通過(guò)常數(shù)項(xiàng)為零的線(xiàn)性擬合,結(jié)果如圖5所示。可以看出,優(yōu)勢(shì)流速與平均流速的比值為1.007,近似相等;因此,在理論上可用此方法計(jì)算的優(yōu)勢(shì)流速作為水流的平均流速。

    3.5 優(yōu)勢(shì)流速與最大流速的關(guān)系

    將流量法測(cè)量得到的平均流速與漂浮物法測(cè)量得到的最大流速進(jìn)行常數(shù)項(xiàng)為零的線(xiàn)性擬合,結(jié)果如圖6所示。可以看出,平均流速為最大流速的0.774倍,平均流速與最大流速的比值即為校正系數(shù)α=0.774,介于前面提到的過(guò)度流和紊流的校正系數(shù)區(qū)間范圍內(nèi),與坡面薄層水流的流態(tài)相吻合。

    圖4 不同坡度與流量條件下擬合電解質(zhì)優(yōu)勢(shì)流速Fig.4 Fitting electrolyte peak velocities under different slopes and different flow rates

    圖5 優(yōu)勢(shì)流速與平均流速的關(guān)系Fig.5 Relationship between peak velocities and mean velocities

    圖6 平均流速與最大流速的關(guān)系Fig.6 Relationship between mean velocities and leading edge velocities

    將水流優(yōu)勢(shì)流速與漂浮物法測(cè)量得到的最大流速,進(jìn)行常數(shù)項(xiàng)為零的線(xiàn)性擬合,結(jié)果如圖7所示??梢钥闯?,水流的優(yōu)勢(shì)流速與最大流速的比值為0.776,與平均流速和最大流速的比值0.774相近;因此水流優(yōu)勢(shì)流速在實(shí)際上可代替平均流速計(jì)算校正系數(shù)α。電解質(zhì)示蹤法計(jì)算水流優(yōu)勢(shì)流速可用于不規(guī)則斷面的水流平均流速測(cè)量,并減小了測(cè)量斷面面積和水深過(guò)程中產(chǎn)生的人為誤差,在理論上優(yōu)于流量法。

    圖7 優(yōu)勢(shì)流速與最大流速的關(guān)系Fig.7 Relationship between peak velocities and leading edge velocities

    4 結(jié)論

    電解質(zhì)示蹤法測(cè)量得到的電解質(zhì)的優(yōu)勢(shì)流速隨距離的增加而增大,最終達(dá)到穩(wěn)定后與水流流速相同,隨距離變化規(guī)律符合指數(shù)函數(shù),采用Grapher進(jìn)行擬合,得到穩(wěn)定的電解質(zhì)優(yōu)勢(shì)流速,即水流優(yōu)勢(shì)流速。各工況下水流優(yōu)勢(shì)流速變化范圍在0.241~0.568m/s之間,流量對(duì)水流優(yōu)勢(shì)流速的影響較坡度對(duì)其的影響更為明顯。優(yōu)勢(shì)流速與平均流速的比值為1.007,在理論上可用優(yōu)勢(shì)流速代替平均流速,優(yōu)勢(shì)流速與最大流速的比值為0.776,與平均流速與最大流速的比值0.774相近,符合坡面薄層水流的流態(tài)特征。試驗(yàn)區(qū)分了電解質(zhì)優(yōu)勢(shì)流速和水流優(yōu)勢(shì)流速,分析了電解質(zhì)注入水流之后隨水流的運(yùn)動(dòng)過(guò)程,結(jié)果可為研究坡面薄層水流動(dòng)力過(guò)程提供新的計(jì)算方法和數(shù)據(jù),為土壤侵蝕研究提供參考。

    [1] 啜瑞媛,雷廷武,史曉楠,等.測(cè)量坡面薄層水流流速的電解質(zhì)示蹤真實(shí)邊界條件法與系統(tǒng)[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(2):77. Chuo Ruiyuan,Lei Tingwu,Shi Xiaonan,et al.Method and system for measuring hill-slope water flow velocity with realistic boundary condition ofelectrolyte tracer[J]. Transactions of the Chinese Society of Agricultural Engineering,2012,28(2):77.(in Chinese)

    [2] 張永東,吳淑芳,馮浩,等.土壤侵蝕過(guò)程中坡面流水力學(xué)特性及侵蝕動(dòng)力研究[J].土壤,2013,45(1):26. Zhang Yongdong,Wu Shufang,F(xiàn)eng Hao,et al.Review on hydraulic characteristics and erosion dynamics of overland flow in soil erosion process[J].Soils,2013,45(1):26.(in Chinese)

    [3] 陳麗燕,雷廷武,董月群,等.不同示蹤法測(cè)量礫石層中水流流速研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(4):141. Chen Liyan,Lei Tingwu,Dong Yuequn,et al.Estimating velocity of water flow within gravel layer by different tracer methods[J].Transactions of the Chinese Society for Agricultural Machinery,2015,46(4):141.(in Chinese)

    [4] Gimenez R,Planchon O,Silvera N,et al.Longitudinal velocity patterns and bed morphology interaction in a rill[J].Earth Surface Process and Landforms,2004,29(1):105.

    [5] Hyun B S,Balachandar R,Yu K,et al.Assessment of PIV to measure mean velocity and turbulence in openchannel flow[J].Experiments in Fluids,2003,35(3):262.

    [6] Robinson K M,Cook K R.Stress measurement upstream of an overfall[J].Transactions of the American Society of Agricultural Engineers,1998,41(4):1019.

    [7] 李永祥,苑明順,李春華.熱膜技術(shù)在水流測(cè)速中的應(yīng)用研究[J].流體力學(xué)實(shí)驗(yàn)與測(cè)量,1997,11(4):45. Li Yongxiang,Yuan Mingshun,Li Chunhua.Hot film technique in the application of flow velocity measurement[J].Experiments and Measurements in Fluid Mechanics,1997,11(4):45.(in Chinese)

    [8] 劉鵬,李小昱,王為,等.基于相關(guān)法的坡面徑流流速測(cè)量系統(tǒng)[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(3):48. Liu Peng,Li Xiaoyu,Wang Wei,et al.Runoff flow velocity measurement system based on correlation method[J].Transactions of the Chinese Society of Agricultural Engineering,2008,24(3):48.(in Chinese)

    [9] 劉鵬,李小昱,王為.基于光電傳感器和示蹤法的徑流流速測(cè)量系統(tǒng)的研究[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(5):116. Liu Peng,Li Xiaoyu,Wang Wei.Runoff flow velocity measurement system using photoelectric sensor and tracing method[J].Transactions of the Chinese Society of Agricultural Engineering,2007,23(5):116.(in Chinese)

    [10]王為,李小昱,張軍,等.基于電導(dǎo)式傳感器徑流流速測(cè)量系統(tǒng)的試驗(yàn)研究[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(2):1. WangWei,Li Xiaoyu,Zhang Jun,et al.Experimental study on measurement system of runoff velocity based on conductance sensor[J].Transactions of the Chinese Society of Agricultural Engineering,2007,23(2):1.(inChinese)

    [11]Luk S H,Merz W.Use of the salt tracing technique to determine the velocity ofoverland flow[J].Soil Technology,1992,5(4):289.

    [12]高素娟,王占禮,黃明斌,等.黃河中游多沙粗沙區(qū)坡面薄層水流水動(dòng)力學(xué)特性[J].水土保持通報(bào),2010,30(4):11. Gao Sujuan,Wang Zhanli,Huang Mingbin,et al.Hydraulic Properties of shallow flow in coarse sediment region of Yellow River middle reaches[J].Bulletin of Soil and Water Conservation,2010,30(4):11.(in Chinese)

    [13]Dunkerley D L.Estimating the mean speed of laminar overland flow using injection uncertainty on rough surface[J].Earth Surface Process Landforms,2001,26(4):363.

    [14]Oliver P,Norbert S,Raphael G,et al.An automated salt tracing gauge for flow velocity measurement[J].Earth Surface Process Landforms,2005,30(7):833.

    [15]Myers TG.Modelling laminar sheet flow over rough surfaces[J].Water Resource Research,2002,38(11):1230.

    [16]Elder JW.The dispersion of marked fluid in turbulent shear flow[J].Journal of Fluid Mechanics,1959(5):544.

    [17]Gilley R E,Kottwitz ER,Simanton JR.Hydraulic characteristic of rill[J].Transactions of the American Society of Agriculture Engineering,1970,33:1900.

    [18]Luk SH,Merz W.Use of the salt tracing technique to determine the velocity ofoverland flow[J].Soil Technology,1992,5:289.

    [19]Horton R E,Leach H R,Vilet V R.Laminar sheet flow[J].Eos Transactions American Geophysical Union,1934,15(2):393.

    [20]Taylor G.The dispersion of matter in turbulent flow through a pipe[J].Proceedings of the Royal Society A,1954,223(1155):446.

    [21]Phelps H O.Shallow laminar flows over rough granular surfaces[J].Journal of the Hydraulics Division,1975,101:367.

    [22]Abrahams A D,Parsons A J,Luk SH.Field measurement of the velocity of overland flow using dye tracing[J]. Earth Surface Processes and Landforms,1986,11(6):653.

    [23]夏衛(wèi)生,劉春平,雷廷武,等.示蹤法測(cè)量坡面水流速度理論缺陷分析[J].自然災(zāi)害學(xué)報(bào),2007,16(1):50. XiaWeisheng,Liu Chunping,Lei Tingwu,et al.Analysis of faults in theory for measuring flow velocity on slope with tracer method[J].Journal of Natural Disasters,2007,16(1):50.(in Chinese)

    [24]夏衛(wèi)生,雷廷武,趙軍.泥沙含量對(duì)鹽液示蹤法經(jīng)驗(yàn)系數(shù)影響的研究[J].農(nóng)業(yè)工程學(xué)報(bào),2003(4):97. XiaWeisheng,Lei Tingwu,Zhao Jun.Research of empirical coefficient of salt tracer method affected by sediment[J].Transactions of the Chinese Society of Agricultural Engineering,2003(4):97.(in Chinese)

    [25]Lei Tingwu,Xia Weisheng,Zhao Jun.Method for measuring velocity of shallow water flow for soil erosion with an electrolyte tracer[J].Journal of Hydrology,2005,301(1/2/3/4):139.

    [26]夏衛(wèi)生,雷廷武,趙軍,等.薄層水流流速測(cè)量系統(tǒng)的研究[J].水科學(xué)進(jìn)展,2003,14(6):781. Xia Weisheng,Lei Tingwu,Zhao Jun,et al.Velocity measuring system for sheet flow[J].Advances in Water Science,2003,14(6):781.(in Chinese)

    [27]禹明忠,詹秀玲,龐東明,等.流場(chǎng)實(shí)時(shí)測(cè)量中示蹤球性能的研究[J].水利水電技術(shù),2002,33(2):43. Yu Mingzhong,Zhan Xiuling,Pang Dongming,et al. Research on the properties of particles in real-time velocity field measurement[J].Water Resources and Hydropower Engineering,2002,33(2):43.(in Chinese)

    [28]董月群,莊曉輝,雷廷武,等.脈沖邊界模型測(cè)量?jī)鐾疗旅鎻搅髁魉倥c距離優(yōu)選[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(2):146. Dong Yuequn,Zhuang Xiaohui,Lei Tingwu,et al.Optimum measurement distance of water flow velocity over frozen slopes with pulse boundary model method[J]. Transactions of the Chinese Society for Agricultural Machinery,2015,46(2):146.(in Chinese)

    Studies on peak velocity of shallow water flow on slopes

    Chen Liyan1,Lei Tingwu1,Chuo Ruiyuan2

    (1.College of Water Resources and Civil Engineering,China Agricultural University,100083,Beijing,China;2.Tianjin Lon win Technology,300181,Tianjin,China)

    [Background]Rainfall runoff is the major dynamic source for hill-slope erosion,runoff velocity is one of the important parameters in soil erosion model.The accurate measurement of shallow water flow velocity is critical in hydrological process.Lei et al.proposed an electrolyte pulse method for measuring the velocity by fitting the solute transport process with time using the least square method and improved the system with a Normal Model and a Sine Model.This laboratory experiments were conducted to determine the relationship between electrolyte peak velocity and water flow peak velocity during the electrolyte transport process and verify the new computational method.[Methods]The experimental devices include a flume,4m long and 15 cm wide,a solute injector,a data logger for control and data acquisition and a computer with specially designed software for data measurement and storage.The experiments involved three flow rates(12,24 and 48 L/min)and three slope gradients(4°,8°and 12°).Five sensors were used to measure the electrolyte transport processes at 0.3 m,0.6 m,0.9 m,1.2m,and 1.5m from the location where the salt solute was injected into the water flow.During each experiment,five complete curves of electrolyte changes with time were recorded,which can be used tocalculate electrolyte peak velocity with the distance from the injection to the measuring sensor and the time used for the peak of the electrolyte to travel through the distance.The leading edge velocities were measured by floating objects method and mean velocities by volumetric method.[Results]The electrolyte peak velocity was between 0.15 to 0.54 m/s,increased with distance and tended to stable,and the flow rates caused greater effect on electrolyte peak velocity than slope gradient under different conditions.The steady electrolyte peak velocity,regarded as the water flow peak velocity,were computed through fitting the electrolyte peak velocity at different distances from the salt solution injector with exponential function,ranging from 0.241 to 0.568 m/s.The exponential function fitted the electrolyte peak velocities very well for all the experimental conditions.The flow rate had greater effect on flow peak velocity growth rate than that of slope gradient.The water flow peak velocity were 1.007 times of mean velocity,0.774 times of leading edge velocity,and mean velocity was 0.776 times of leading edge velocity.[Conclusions]The flow peak velocity agreed well with leading edge velocity and mean velocity.These demonstrated that the new computational method for measuring shallow water flow velocity was reasonable and valid.The results provide a new method for computing the mean velocity of sheet flow and relevant data for the dynamic process of sheet flow,which will be useful for the investigation of soil erosion.

    shallow water flow;electrolyte tracer method;flow method;peak velocity;mean velocity

    S157.9

    A

    1672-3007(2016)05-0130-08

    10.16843/j.sswc.2016.05.017

    2015- 10- 28

    2016- 09- 21

    項(xiàng)目名稱(chēng):國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目“高海拔寒區(qū)融水土壤侵蝕機(jī)理與過(guò)程模擬研究”(41230746)

    陳麗燕(1980—),女,博士研究生。主要研究方向:土壤侵蝕。E-mail:cly508@163.com

    ?通信作者簡(jiǎn)介:雷廷武(1958—),男,博士,博士生導(dǎo)師。主要研究方向:土壤侵蝕與旱地農(nóng)業(yè)。E-mail:leitingwu@cau.edu. cn

    猜你喜歡
    薄層坡面坡度
    關(guān)于公路超高漸變段合成坡度解析與應(yīng)用
    沖積扇油氣管道坡面侵蝕災(zāi)害因子分析
    超音速流越過(guò)彎曲坡面的反問(wèn)題
    維藥芹菜根的薄層鑒別
    SiN_x:H膜沉積壓強(qiáng)與擴(kuò)散薄層電阻的匹配性研究
    面板堆石壩墊層施工及坡面防護(hù)
    基于圖像處理的定位器坡度計(jì)算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    參芪苓口服液的薄層色譜鑒別
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    CT和MR對(duì)人上脛腓關(guān)節(jié)面坡度的比較研究
    亚洲专区中文字幕在线| 成人18禁在线播放| 国产成人av教育| 嫁个100分男人电影在线观看| 午夜福利免费观看在线| 国产精品久久久人人做人人爽| 午夜福利高清视频| 亚洲电影在线观看av| 成人精品一区二区免费| 一本精品99久久精品77| 精品日产1卡2卡| 人人妻人人澡人人看| 久久九九热精品免费| 一个人观看的视频www高清免费观看 | 激情在线观看视频在线高清| 国内久久婷婷六月综合欲色啪| 免费在线观看黄色视频的| 在线观看免费日韩欧美大片| 18禁美女被吸乳视频| 丝袜美腿诱惑在线| 久久香蕉国产精品| 波多野结衣高清作品| 日韩欧美一区视频在线观看| 男人舔奶头视频| 麻豆久久精品国产亚洲av| 精品不卡国产一区二区三区| 757午夜福利合集在线观看| 久久青草综合色| 国产高清videossex| 成人三级做爰电影| 成熟少妇高潮喷水视频| 欧美一级毛片孕妇| 亚洲人成伊人成综合网2020| 99精品久久久久人妻精品| tocl精华| 大香蕉久久成人网| 久久精品91无色码中文字幕| 91av网站免费观看| 久久香蕉国产精品| 亚洲国产欧洲综合997久久, | 免费在线观看亚洲国产| 欧美色视频一区免费| 桃红色精品国产亚洲av| 国产精品久久久久久人妻精品电影| 亚洲欧美一区二区三区黑人| 宅男免费午夜| 久久天堂一区二区三区四区| 在线观看一区二区三区| 亚洲第一av免费看| 久久久国产欧美日韩av| 久久久久久人人人人人| 一区二区三区高清视频在线| 每晚都被弄得嗷嗷叫到高潮| 99久久久亚洲精品蜜臀av| 身体一侧抽搐| 亚洲精品在线美女| 夜夜看夜夜爽夜夜摸| 色哟哟哟哟哟哟| 久久久久久免费高清国产稀缺| 国产亚洲精品久久久久久毛片| 精品不卡国产一区二区三区| 精品少妇一区二区三区视频日本电影| 亚洲av五月六月丁香网| 白带黄色成豆腐渣| 亚洲欧美精品综合一区二区三区| 国产乱人伦免费视频| 亚洲国产精品合色在线| 99国产精品一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | e午夜精品久久久久久久| 午夜福利在线观看吧| 99久久无色码亚洲精品果冻| 亚洲aⅴ乱码一区二区在线播放 | videosex国产| 变态另类丝袜制服| АⅤ资源中文在线天堂| 老汉色av国产亚洲站长工具| 国产午夜精品久久久久久| 日韩三级视频一区二区三区| 波多野结衣高清作品| 国产黄a三级三级三级人| 亚洲精品国产精品久久久不卡| АⅤ资源中文在线天堂| 19禁男女啪啪无遮挡网站| 啦啦啦 在线观看视频| 黄色视频不卡| 亚洲人成77777在线视频| 男女之事视频高清在线观看| 色播亚洲综合网| 亚洲五月天丁香| 又大又爽又粗| 制服诱惑二区| 一进一出抽搐动态| 在线av久久热| 校园春色视频在线观看| 欧美另类亚洲清纯唯美| 在线播放国产精品三级| 国产又黄又爽又无遮挡在线| 91老司机精品| 日日夜夜操网爽| 黑人操中国人逼视频| 欧美日本视频| 中文字幕精品免费在线观看视频| 天天一区二区日本电影三级| 欧美精品亚洲一区二区| 久久国产精品影院| 亚洲中文字幕一区二区三区有码在线看 | 国产成人一区二区三区免费视频网站| 国产日本99.免费观看| 成人国语在线视频| 国产精品 欧美亚洲| 亚洲天堂国产精品一区在线| 国内精品久久久久精免费| 亚洲国产欧美网| 免费在线观看视频国产中文字幕亚洲| 免费在线观看视频国产中文字幕亚洲| 老熟妇乱子伦视频在线观看| 亚洲 欧美 日韩 在线 免费| 黄色 视频免费看| 亚洲av成人不卡在线观看播放网| 久久久久亚洲av毛片大全| 午夜福利18| 久久婷婷人人爽人人干人人爱| 国产精品二区激情视频| 国产亚洲欧美在线一区二区| 国产一区二区三区在线臀色熟女| 亚洲九九香蕉| 麻豆av在线久日| 久久香蕉激情| 女人被狂操c到高潮| 丰满人妻熟妇乱又伦精品不卡| 窝窝影院91人妻| 国产黄a三级三级三级人| 韩国精品一区二区三区| 男女午夜视频在线观看| 国产主播在线观看一区二区| 亚洲,欧美精品.| 成人亚洲精品一区在线观看| 久久久久久久精品吃奶| 国产亚洲av嫩草精品影院| 久久午夜亚洲精品久久| 人人妻人人澡人人看| 高清在线国产一区| 日日摸夜夜添夜夜添小说| 搞女人的毛片| 亚洲 欧美 日韩 在线 免费| 久久精品人妻少妇| 日韩 欧美 亚洲 中文字幕| 波多野结衣高清作品| 亚洲人成77777在线视频| 夜夜爽天天搞| 成年人黄色毛片网站| 国产精品亚洲av一区麻豆| √禁漫天堂资源中文www| 精品第一国产精品| 黄色毛片三级朝国网站| 日韩国内少妇激情av| 成人特级黄色片久久久久久久| 少妇粗大呻吟视频| 婷婷精品国产亚洲av在线| 一区二区日韩欧美中文字幕| 亚洲精品美女久久久久99蜜臀| 国产亚洲欧美在线一区二区| 两个人免费观看高清视频| 国产精品亚洲一级av第二区| 午夜视频精品福利| 一区福利在线观看| 亚洲人成77777在线视频| 精品久久久久久久末码| 国产成人av教育| xxx96com| 成人永久免费在线观看视频| 午夜两性在线视频| 欧美又色又爽又黄视频| 国产免费av片在线观看野外av| 女性生殖器流出的白浆| av超薄肉色丝袜交足视频| 日本一区二区免费在线视频| 婷婷丁香在线五月| 淫妇啪啪啪对白视频| 在线十欧美十亚洲十日本专区| 999精品在线视频| 色在线成人网| 免费高清视频大片| 满18在线观看网站| 18禁黄网站禁片午夜丰满| 熟女少妇亚洲综合色aaa.| 免费观看精品视频网站| 两性午夜刺激爽爽歪歪视频在线观看 | 精品欧美国产一区二区三| 精品久久蜜臀av无| 欧美久久黑人一区二区| 欧美乱色亚洲激情| 亚洲一区二区三区色噜噜| 村上凉子中文字幕在线| 十八禁网站免费在线| 国产成人精品无人区| 99在线视频只有这里精品首页| 久久草成人影院| 我的亚洲天堂| 亚洲最大成人中文| 久久久久久九九精品二区国产 | 少妇粗大呻吟视频| 亚洲国产看品久久| 一个人免费在线观看的高清视频| 色播亚洲综合网| 国内精品久久久久久久电影| 成人18禁高潮啪啪吃奶动态图| 99国产精品一区二区三区| 国产爱豆传媒在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 老汉色av国产亚洲站长工具| 国产成人精品无人区| 午夜激情福利司机影院| 国产精品 欧美亚洲| 美国免费a级毛片| 身体一侧抽搐| 国产主播在线观看一区二区| 亚洲性夜色夜夜综合| 青草久久国产| 老司机福利观看| 热re99久久国产66热| 日本黄色视频三级网站网址| 精品国产一区二区三区四区第35| 国产精品九九99| 国产精品久久视频播放| 性欧美人与动物交配| 中文字幕高清在线视频| 国产精品1区2区在线观看.| 国产精品久久久av美女十八| 成人av一区二区三区在线看| 一区二区日韩欧美中文字幕| 国产一区二区三区在线臀色熟女| 午夜福利视频1000在线观看| 免费在线观看成人毛片| 91老司机精品| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av成人不卡在线观看播放网| 国产91精品成人一区二区三区| 日本免费a在线| 久久久久久亚洲精品国产蜜桃av| 日韩欧美国产在线观看| 婷婷亚洲欧美| 日韩欧美国产一区二区入口| 天堂动漫精品| 午夜福利视频1000在线观看| 一卡2卡三卡四卡精品乱码亚洲| 一级毛片精品| 夜夜看夜夜爽夜夜摸| 国内揄拍国产精品人妻在线 | 一级作爱视频免费观看| 午夜成年电影在线免费观看| 久久国产亚洲av麻豆专区| 天天一区二区日本电影三级| 99国产极品粉嫩在线观看| 黄片小视频在线播放| 亚洲五月天丁香| 人妻久久中文字幕网| 国产激情偷乱视频一区二区| 一区二区三区激情视频| 国产主播在线观看一区二区| 香蕉丝袜av| 性欧美人与动物交配| 91九色精品人成在线观看| 12—13女人毛片做爰片一| 黄片小视频在线播放| 亚洲五月天丁香| 校园春色视频在线观看| 国产伦一二天堂av在线观看| 亚洲av成人不卡在线观看播放网| 天堂动漫精品| 亚洲五月婷婷丁香| 国产高清videossex| 757午夜福利合集在线观看| 午夜影院日韩av| 国产精品自产拍在线观看55亚洲| 好男人在线观看高清免费视频 | 日本 av在线| 满18在线观看网站| 久久久国产成人精品二区| 国产亚洲精品一区二区www| 99久久久亚洲精品蜜臀av| 88av欧美| 在线观看午夜福利视频| 成人18禁高潮啪啪吃奶动态图| 色婷婷久久久亚洲欧美| 美女免费视频网站| 亚洲av第一区精品v没综合| 国产精品亚洲av一区麻豆| 欧美绝顶高潮抽搐喷水| 人人妻人人澡欧美一区二区| 欧美性猛交╳xxx乱大交人| 无人区码免费观看不卡| 免费在线观看日本一区| 国产成人欧美| 岛国视频午夜一区免费看| 18美女黄网站色大片免费观看| 欧洲精品卡2卡3卡4卡5卡区| 免费无遮挡裸体视频| 久久狼人影院| 亚洲七黄色美女视频| x7x7x7水蜜桃| 热re99久久国产66热| 伊人久久大香线蕉亚洲五| 91成年电影在线观看| 操出白浆在线播放| 成人一区二区视频在线观看| 欧美一区二区精品小视频在线| 久久久水蜜桃国产精品网| 麻豆国产av国片精品| 国产精品电影一区二区三区| 亚洲 欧美 日韩 在线 免费| 日本精品一区二区三区蜜桃| 男人的好看免费观看在线视频| 一本精品99久久精品77| 精品不卡国产一区二区三区| 日日撸夜夜添| 亚洲四区av| 日韩一区二区视频免费看| 22中文网久久字幕| a级毛片a级免费在线| 国产视频内射| 欧美日韩国产亚洲二区| 搡老岳熟女国产| 丝袜美腿在线中文| 国产精品永久免费网站| 一个人看的www免费观看视频| 久久这里只有精品中国| av卡一久久| 亚洲图色成人| 精品不卡国产一区二区三区| 在线免费观看的www视频| 久久久色成人| 国产伦精品一区二区三区四那| 久久欧美精品欧美久久欧美| av在线天堂中文字幕| 欧美丝袜亚洲另类| 特级一级黄色大片| www日本黄色视频网| 亚洲精品影视一区二区三区av| 级片在线观看| 中文亚洲av片在线观看爽| av免费在线看不卡| 国产精华一区二区三区| 精品久久久久久久久久免费视频| 日本黄大片高清| 日韩三级伦理在线观看| 精品一区二区三区视频在线| 中文字幕av在线有码专区| 欧美日本亚洲视频在线播放| 亚洲国产精品成人综合色| 久久久久久久久久久丰满| 亚洲中文字幕一区二区三区有码在线看| 久久精品影院6| 热99re8久久精品国产| 一个人看视频在线观看www免费| 午夜a级毛片| 一本久久中文字幕| 国产av麻豆久久久久久久| 午夜福利在线在线| 久久婷婷人人爽人人干人人爱| 热99re8久久精品国产| 伊人久久精品亚洲午夜| www日本黄色视频网| 欧美成人a在线观看| 国产综合懂色| 日本-黄色视频高清免费观看| 亚洲精品粉嫩美女一区| 桃色一区二区三区在线观看| 99久久九九国产精品国产免费| 久久这里只有精品中国| 给我免费播放毛片高清在线观看| 精品久久久久久久人妻蜜臀av| 午夜福利在线在线| 亚洲aⅴ乱码一区二区在线播放| av黄色大香蕉| 国产高清视频在线播放一区| 久久久精品欧美日韩精品| 免费看日本二区| 99久久九九国产精品国产免费| 国产精品一及| 美女 人体艺术 gogo| 欧美激情国产日韩精品一区| 成人国产麻豆网| 久久久色成人| 国产精品不卡视频一区二区| 国产不卡一卡二| 毛片女人毛片| 久久亚洲国产成人精品v| 五月伊人婷婷丁香| 亚洲婷婷狠狠爱综合网| 日韩大尺度精品在线看网址| 日韩精品青青久久久久久| 国产一区二区三区在线臀色熟女| 91精品国产九色| av.在线天堂| 深爱激情五月婷婷| 亚洲图色成人| 91久久精品电影网| 亚洲欧美日韩卡通动漫| 国产大屁股一区二区在线视频| 12—13女人毛片做爰片一| 亚洲av一区综合| 精品久久久久久久久久久久久| 国产极品精品免费视频能看的| 丰满乱子伦码专区| 精品午夜福利视频在线观看一区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av成人av| 免费看日本二区| 日本-黄色视频高清免费观看| 丰满乱子伦码专区| 亚洲人成网站在线播| 国产免费一级a男人的天堂| 搡老妇女老女人老熟妇| 精品久久久久久久久久久久久| 悠悠久久av| 一个人免费在线观看电影| 国内久久婷婷六月综合欲色啪| 成人av一区二区三区在线看| av在线播放精品| 丝袜喷水一区| 亚洲av中文av极速乱| 亚洲婷婷狠狠爱综合网| 亚洲人成网站在线观看播放| 国产不卡一卡二| 给我免费播放毛片高清在线观看| 国产伦一二天堂av在线观看| 日韩欧美三级三区| 女人被狂操c到高潮| 免费观看在线日韩| 舔av片在线| 国产又黄又爽又无遮挡在线| 国产精品福利在线免费观看| 国产精品久久视频播放| 亚洲中文日韩欧美视频| 久久久久久久久大av| 国模一区二区三区四区视频| 亚洲精品乱码久久久v下载方式| 精品少妇黑人巨大在线播放 | 91麻豆精品激情在线观看国产| 国产黄a三级三级三级人| 久久午夜福利片| 中文字幕免费在线视频6| 国产精品亚洲美女久久久| 12—13女人毛片做爰片一| 97超视频在线观看视频| 国产成人a∨麻豆精品| 韩国av在线不卡| 小蜜桃在线观看免费完整版高清| 免费在线观看影片大全网站| 日本 av在线| 91在线精品国自产拍蜜月| 在线看三级毛片| 麻豆一二三区av精品| 国内精品美女久久久久久| 亚州av有码| 97超级碰碰碰精品色视频在线观看| 婷婷亚洲欧美| 国产精品日韩av在线免费观看| 国产v大片淫在线免费观看| 国内精品美女久久久久久| 两个人视频免费观看高清| 精品久久久久久久久久免费视频| 无遮挡黄片免费观看| 男人和女人高潮做爰伦理| 欧美绝顶高潮抽搐喷水| АⅤ资源中文在线天堂| 亚洲av电影不卡..在线观看| 91久久精品电影网| 日韩av在线大香蕉| 国产av麻豆久久久久久久| 国内揄拍国产精品人妻在线| 非洲黑人性xxxx精品又粗又长| .国产精品久久| 欧美性感艳星| 又爽又黄a免费视频| 搡老岳熟女国产| 国产精品嫩草影院av在线观看| 网址你懂的国产日韩在线| 免费高清视频大片| 别揉我奶头 嗯啊视频| 亚洲激情五月婷婷啪啪| 18禁裸乳无遮挡免费网站照片| 亚洲成av人片在线播放无| 国产精品,欧美在线| 日本撒尿小便嘘嘘汇集6| 永久网站在线| 亚洲七黄色美女视频| 午夜日韩欧美国产| 日韩欧美三级三区| 91精品国产九色| 中文字幕av成人在线电影| 一区二区三区免费毛片| 天天一区二区日本电影三级| av.在线天堂| 久久久精品欧美日韩精品| 老师上课跳d突然被开到最大视频| 最好的美女福利视频网| 1024手机看黄色片| 白带黄色成豆腐渣| 中国国产av一级| 啦啦啦韩国在线观看视频| 热99在线观看视频| 色视频www国产| 免费在线观看影片大全网站| av卡一久久| 免费看光身美女| 国产精品日韩av在线免费观看| 天堂影院成人在线观看| 少妇丰满av| 夜夜夜夜夜久久久久| 日本黄大片高清| 久久天躁狠狠躁夜夜2o2o| 伦精品一区二区三区| 精品久久久久久久久av| 午夜久久久久精精品| 一a级毛片在线观看| 久久精品夜色国产| 国产一级毛片七仙女欲春2| 天天躁日日操中文字幕| 国产av一区在线观看免费| 国产免费男女视频| 欧美性猛交黑人性爽| 亚洲在线自拍视频| 欧美bdsm另类| 成人二区视频| 亚洲精品日韩在线中文字幕 | 欧美日韩乱码在线| 中文字幕av在线有码专区| 国产av不卡久久| 男女那种视频在线观看| 九九久久精品国产亚洲av麻豆| 女人被狂操c到高潮| 精品熟女少妇av免费看| 长腿黑丝高跟| 观看免费一级毛片| 久久精品夜夜夜夜夜久久蜜豆| 欧美成人精品欧美一级黄| 人人妻人人澡欧美一区二区| 亚洲国产精品成人久久小说 | 人妻制服诱惑在线中文字幕| 色综合站精品国产| 欧美最黄视频在线播放免费| 深夜精品福利| 乱码一卡2卡4卡精品| 精品国产三级普通话版| 美女高潮的动态| 亚洲美女黄片视频| 一级黄片播放器| 日韩精品青青久久久久久| 国产人妻一区二区三区在| 国内精品美女久久久久久| 国产精品久久视频播放| 婷婷精品国产亚洲av| 亚洲精品影视一区二区三区av| 精品福利观看| 男女视频在线观看网站免费| 亚洲国产精品合色在线| 日本爱情动作片www.在线观看 | 成人特级av手机在线观看| 中文在线观看免费www的网站| 成人国产麻豆网| 欧美成人精品欧美一级黄| 国产 一区精品| 男人狂女人下面高潮的视频| 亚洲欧美日韩无卡精品| 最新在线观看一区二区三区| 自拍偷自拍亚洲精品老妇| 一级黄片播放器| 国产三级在线视频| 成人鲁丝片一二三区免费| 国产真实伦视频高清在线观看| 一级av片app| 国产一区二区三区av在线 | 秋霞在线观看毛片| 波多野结衣高清无吗| 国产成人aa在线观看| 国产高清视频在线播放一区| 亚洲国产色片| 久久久久久九九精品二区国产| 能在线免费观看的黄片| 精品人妻偷拍中文字幕| 国产麻豆成人av免费视频| 欧美+亚洲+日韩+国产| 成年免费大片在线观看| av中文乱码字幕在线| 女人被狂操c到高潮| 3wmmmm亚洲av在线观看| 国产91av在线免费观看| 亚洲av.av天堂| 乱码一卡2卡4卡精品| 日本精品一区二区三区蜜桃| 99精品在免费线老司机午夜| 国产精品日韩av在线免费观看| 精品久久久久久久人妻蜜臀av| 亚洲精品国产成人久久av| 亚洲欧美成人精品一区二区| 18禁在线播放成人免费| 国产伦一二天堂av在线观看| 日韩成人av中文字幕在线观看 | 国产一区二区亚洲精品在线观看| 国产伦精品一区二区三区四那| 国国产精品蜜臀av免费| 少妇的逼水好多| 国产一区二区亚洲精品在线观看| 看黄色毛片网站| 欧美一区二区亚洲| 成人高潮视频无遮挡免费网站| 免费搜索国产男女视频| 小说图片视频综合网站| 亚洲最大成人手机在线| 国产高清视频在线播放一区| 少妇人妻精品综合一区二区 | 国产精品永久免费网站|