• <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é)面坡度的比較研究
    国产精品香港三级国产av潘金莲| 在线观看免费视频日本深夜| 成人国产综合亚洲| 国产精品久久久av美女十八| 国内精品久久久久久久电影| 巨乳人妻的诱惑在线观看| 日本三级黄在线观看| 婷婷精品国产亚洲av| 中亚洲国语对白在线视频| 亚洲 国产 在线| 91在线观看av| 精品日产1卡2卡| a级毛片a级免费在线| 久久人人精品亚洲av| 色婷婷久久久亚洲欧美| 日韩有码中文字幕| 少妇被粗大的猛进出69影院| 老熟妇乱子伦视频在线观看| 欧美一级a爱片免费观看看 | 日本黄色视频三级网站网址| 精品一区二区三区av网在线观看| 777久久人妻少妇嫩草av网站| 夜夜爽天天搞| 婷婷精品国产亚洲av在线| 制服人妻中文乱码| 天堂动漫精品| 午夜日韩欧美国产| 久久人妻福利社区极品人妻图片| 欧美久久黑人一区二区| 久久狼人影院| 在线国产一区二区在线| 国产精品美女特级片免费视频播放器 | 性色av乱码一区二区三区2| 最好的美女福利视频网| 男人的好看免费观看在线视频 | 黄片大片在线免费观看| 国产三级在线视频| 国产成人精品久久二区二区免费| 制服诱惑二区| 国内毛片毛片毛片毛片毛片| av在线播放免费不卡| 熟女电影av网| 日日爽夜夜爽网站| 亚洲免费av在线视频| 最近最新中文字幕大全免费视频| 日韩精品免费视频一区二区三区| 不卡一级毛片| 午夜日韩欧美国产| 国产久久久一区二区三区| 亚洲成人精品中文字幕电影| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片女人18水好多| 久久精品夜夜夜夜夜久久蜜豆 | 精品久久蜜臀av无| 欧美精品亚洲一区二区| 国产不卡一卡二| 色老头精品视频在线观看| 日本三级黄在线观看| 国产成人一区二区三区免费视频网站| 一夜夜www| 亚洲av美国av| 老熟妇仑乱视频hdxx| 脱女人内裤的视频| 久久久久国内视频| 亚洲av成人不卡在线观看播放网| 中文字幕久久专区| 免费高清视频大片| 精品久久久久久久久久免费视频| 日韩欧美三级三区| 国产熟女xx| 久久香蕉激情| 日韩欧美一区视频在线观看| 啪啪无遮挡十八禁网站| 国产精品乱码一区二三区的特点| 欧美+亚洲+日韩+国产| 精品一区二区三区四区五区乱码| 亚洲aⅴ乱码一区二区在线播放 | 国内毛片毛片毛片毛片毛片| 黄色女人牲交| 国产男靠女视频免费网站| 成人特级黄色片久久久久久久| 又黄又粗又硬又大视频| 免费在线观看日本一区| 亚洲一区高清亚洲精品| 青草久久国产| 黄片小视频在线播放| 啦啦啦观看免费观看视频高清| 一本一本综合久久| 黄色视频不卡| 欧美av亚洲av综合av国产av| 两个人看的免费小视频| 无限看片的www在线观看| 亚洲av美国av| 国产蜜桃级精品一区二区三区| 动漫黄色视频在线观看| 又黄又爽又免费观看的视频| 99精品在免费线老司机午夜| 高潮久久久久久久久久久不卡| 两个人看的免费小视频| 最新在线观看一区二区三区| 午夜福利视频1000在线观看| 18美女黄网站色大片免费观看| 1024视频免费在线观看| 听说在线观看完整版免费高清| 久久久久久大精品| 露出奶头的视频| 国产午夜精品久久久久久| 国产一区二区激情短视频| 欧美性猛交黑人性爽| 欧美黑人精品巨大| 少妇熟女aⅴ在线视频| 99热6这里只有精品| or卡值多少钱| 黄片小视频在线播放| 精品不卡国产一区二区三区| 午夜久久久在线观看| 很黄的视频免费| 国产亚洲欧美在线一区二区| 亚洲国产欧洲综合997久久, | 久久婷婷人人爽人人干人人爱| 美女午夜性视频免费| 欧美av亚洲av综合av国产av| 国内揄拍国产精品人妻在线 | 男女之事视频高清在线观看| 狠狠狠狠99中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 男女之事视频高清在线观看| 琪琪午夜伦伦电影理论片6080| 国产免费av片在线观看野外av| 亚洲国产高清在线一区二区三 | 久久亚洲真实| 母亲3免费完整高清在线观看| 国产一区二区三区在线臀色熟女| 成人免费观看视频高清| 国产精品免费视频内射| 久久精品国产清高在天天线| 欧美 亚洲 国产 日韩一| 给我免费播放毛片高清在线观看| 每晚都被弄得嗷嗷叫到高潮| 美国免费a级毛片| 欧美日本视频| 熟妇人妻久久中文字幕3abv| 日韩欧美 国产精品| 日韩 欧美 亚洲 中文字幕| 欧美日韩中文字幕国产精品一区二区三区| 一本久久中文字幕| aaaaa片日本免费| 操出白浆在线播放| 久久久精品国产亚洲av高清涩受| 大型av网站在线播放| 女性生殖器流出的白浆| 在线永久观看黄色视频| 人人妻人人澡欧美一区二区| 一级黄色大片毛片| 免费在线观看视频国产中文字幕亚洲| 亚洲精品av麻豆狂野| 精品国产乱子伦一区二区三区| 欧美国产日韩亚洲一区| xxxwww97欧美| 国产成人欧美在线观看| 婷婷精品国产亚洲av在线| 看黄色毛片网站| 丰满人妻熟妇乱又伦精品不卡| 熟女少妇亚洲综合色aaa.| 成人国语在线视频| 91麻豆av在线| 美女扒开内裤让男人捅视频| 人人澡人人妻人| 午夜日韩欧美国产| 欧美性长视频在线观看| 熟妇人妻久久中文字幕3abv| 高潮久久久久久久久久久不卡| 国产精品久久久人人做人人爽| 久久中文看片网| 一级a爱片免费观看的视频| 亚洲片人在线观看| 午夜亚洲福利在线播放| 国产成人精品久久二区二区91| 亚洲色图 男人天堂 中文字幕| 人妻久久中文字幕网| 91九色精品人成在线观看| 一级a爱视频在线免费观看| 99国产精品99久久久久| 91老司机精品| 欧美中文综合在线视频| 欧美黑人巨大hd| 一区二区日韩欧美中文字幕| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲精品一区二区www| aaaaa片日本免费| 夜夜爽天天搞| www.999成人在线观看| 日本免费一区二区三区高清不卡| 香蕉久久夜色| 久久久久精品国产欧美久久久| 中出人妻视频一区二区| 校园春色视频在线观看| 制服诱惑二区| 国产成人精品久久二区二区91| 美女大奶头视频| 免费无遮挡裸体视频| 欧美黑人欧美精品刺激| 日本精品一区二区三区蜜桃| 免费在线观看完整版高清| 亚洲精品国产精品久久久不卡| 国产精品爽爽va在线观看网站 | 亚洲三区欧美一区| 成人永久免费在线观看视频| 看免费av毛片| 日本在线视频免费播放| 日韩欧美国产在线观看| 最近最新中文字幕大全电影3 | 国产区一区二久久| 免费在线观看黄色视频的| 国产国语露脸激情在线看| 女人被狂操c到高潮| 日本三级黄在线观看| 怎么达到女性高潮| www日本在线高清视频| 一级黄色大片毛片| 18禁观看日本| 国产av在哪里看| 18禁黄网站禁片免费观看直播| 草草在线视频免费看| 天堂动漫精品| 黑丝袜美女国产一区| 听说在线观看完整版免费高清| а√天堂www在线а√下载| 夜夜爽天天搞| 两个人免费观看高清视频| 男人舔女人下体高潮全视频| 亚洲精品在线美女| 中文字幕高清在线视频| 国产精品久久电影中文字幕| 嫁个100分男人电影在线观看| 极品教师在线免费播放| 国产成年人精品一区二区| 黄色 视频免费看| 神马国产精品三级电影在线观看 | 欧美另类亚洲清纯唯美| 两个人看的免费小视频| 色播在线永久视频| 亚洲av五月六月丁香网| 动漫黄色视频在线观看| 日韩欧美国产一区二区入口| 午夜免费成人在线视频| 久久亚洲精品不卡| 亚洲国产欧洲综合997久久, | x7x7x7水蜜桃| 午夜久久久在线观看| 亚洲人成电影免费在线| 国产精品免费一区二区三区在线| 国产真人三级小视频在线观看| 国产av一区二区精品久久| av在线播放免费不卡| 黄色女人牲交| 亚洲国产欧洲综合997久久, | 精品午夜福利视频在线观看一区| 校园春色视频在线观看| 日本一本二区三区精品| 看黄色毛片网站| 亚洲一区二区三区不卡视频| 亚洲午夜精品一区,二区,三区| 少妇粗大呻吟视频| 美女免费视频网站| 中文字幕人成人乱码亚洲影| 日韩三级视频一区二区三区| 中文字幕精品免费在线观看视频| 久久中文字幕人妻熟女| 亚洲精品色激情综合| 久久久久久久久久黄片| 日韩三级视频一区二区三区| 最新美女视频免费是黄的| 欧美黄色淫秽网站| 高清在线国产一区| 日日爽夜夜爽网站| 日本成人三级电影网站| 人人妻人人澡欧美一区二区| 黄片小视频在线播放| 亚洲一码二码三码区别大吗| 国产欧美日韩一区二区精品| 欧美乱妇无乱码| 最新在线观看一区二区三区| 亚洲人成网站高清观看| 99国产精品一区二区三区| 久久久久久人人人人人| 久久精品aⅴ一区二区三区四区| 日日干狠狠操夜夜爽| 国产成人av教育| 亚洲精品国产一区二区精华液| 久久久久久久久中文| 熟女电影av网| 国产主播在线观看一区二区| av视频在线观看入口| 欧美精品亚洲一区二区| 午夜福利在线在线| 女性生殖器流出的白浆| www.www免费av| 精品久久蜜臀av无| 欧美最黄视频在线播放免费| 一卡2卡三卡四卡精品乱码亚洲| 99热6这里只有精品| 免费看日本二区| 性欧美人与动物交配| 久久精品国产99精品国产亚洲性色| 午夜精品在线福利| 亚洲激情在线av| 国产精品久久久久久亚洲av鲁大| 久久 成人 亚洲| av片东京热男人的天堂| 国产亚洲精品久久久久5区| 亚洲国产精品sss在线观看| 好男人电影高清在线观看| 色综合婷婷激情| 亚洲av成人一区二区三| www国产在线视频色| 欧美大码av| 久久精品亚洲精品国产色婷小说| 日韩精品青青久久久久久| 长腿黑丝高跟| 久久性视频一级片| 正在播放国产对白刺激| 久久久久国产精品人妻aⅴ院| 搡老岳熟女国产| 亚洲欧美一区二区三区黑人| 色av中文字幕| 国产精品99久久99久久久不卡| 国产一区在线观看成人免费| 香蕉丝袜av| 国产97色在线日韩免费| 一级毛片高清免费大全| 欧美+亚洲+日韩+国产| 啦啦啦免费观看视频1| 黄色视频不卡| 欧美日韩一级在线毛片| 中出人妻视频一区二区| 黑人欧美特级aaaaaa片| 欧美激情久久久久久爽电影| 国产精品99久久99久久久不卡| 黄色成人免费大全| 久久久国产欧美日韩av| 日韩三级视频一区二区三区| 这个男人来自地球电影免费观看| 老司机福利观看| 熟妇人妻久久中文字幕3abv| 国产精品影院久久| 十八禁网站免费在线| 日韩大码丰满熟妇| 日本三级黄在线观看| 搞女人的毛片| 午夜福利在线在线| 天堂√8在线中文| 高清毛片免费观看视频网站| 色综合欧美亚洲国产小说| 亚洲国产欧美网| 日本免费一区二区三区高清不卡| 成人三级黄色视频| 淫秽高清视频在线观看| 国产av又大| 精品一区二区三区av网在线观看| 12—13女人毛片做爰片一| 波多野结衣高清作品| 精品熟女少妇八av免费久了| 中文在线观看免费www的网站 | 精品久久久久久久久久久久久 | 波多野结衣av一区二区av| 99热只有精品国产| 国内毛片毛片毛片毛片毛片| 久久精品成人免费网站| 精品久久久久久久久久久久久 | 亚洲国产精品999在线| 男女之事视频高清在线观看| 正在播放国产对白刺激| 亚洲电影在线观看av| 国产精品久久久人人做人人爽| 亚洲七黄色美女视频| 97超级碰碰碰精品色视频在线观看| 午夜日韩欧美国产| 国产av在哪里看| 国产高清有码在线观看视频 | 亚洲成av片中文字幕在线观看| 黄片大片在线免费观看| 欧美久久黑人一区二区| 悠悠久久av| 国产一区二区三区在线臀色熟女| 国产一区二区激情短视频| 欧美 亚洲 国产 日韩一| 人人妻人人澡欧美一区二区| 久久中文字幕一级| 日本在线视频免费播放| 老司机在亚洲福利影院| 日本精品一区二区三区蜜桃| 18禁黄网站禁片午夜丰满| 欧美乱码精品一区二区三区| 久久久国产欧美日韩av| 欧美黑人巨大hd| 国产精品一区二区免费欧美| 大型黄色视频在线免费观看| 日韩视频一区二区在线观看| 中出人妻视频一区二区| 久热爱精品视频在线9| 国产av不卡久久| 在线观看免费午夜福利视频| 久久精品亚洲精品国产色婷小说| 欧美一区二区精品小视频在线| 久久精品影院6| 国产真实乱freesex| 黄片播放在线免费| 国产v大片淫在线免费观看| 美女扒开内裤让男人捅视频| 黄色片一级片一级黄色片| 搞女人的毛片| 啦啦啦 在线观看视频| 日本三级黄在线观看| 人人妻,人人澡人人爽秒播| 男女床上黄色一级片免费看| 首页视频小说图片口味搜索| 成人特级黄色片久久久久久久| 一个人免费在线观看的高清视频| 校园春色视频在线观看| 久久久久久免费高清国产稀缺| 成人免费观看视频高清| 十分钟在线观看高清视频www| 黄片播放在线免费| 欧美性猛交╳xxx乱大交人| 亚洲国产中文字幕在线视频| 一区二区三区高清视频在线| 亚洲av日韩精品久久久久久密| 熟女少妇亚洲综合色aaa.| 大型黄色视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放 | 91国产中文字幕| 9191精品国产免费久久| 一a级毛片在线观看| 精品无人区乱码1区二区| 999久久久精品免费观看国产| 亚洲人成77777在线视频| 怎么达到女性高潮| 久热这里只有精品99| 国产区一区二久久| 动漫黄色视频在线观看| 亚洲七黄色美女视频| 最新美女视频免费是黄的| 久久精品国产清高在天天线| 国产成人av教育| 黄色成人免费大全| 99国产综合亚洲精品| www.熟女人妻精品国产| 久久精品国产亚洲av香蕉五月| 亚洲国产中文字幕在线视频| 欧美av亚洲av综合av国产av| 夜夜夜夜夜久久久久| 韩国精品一区二区三区| 国产精品亚洲av一区麻豆| 欧美国产精品va在线观看不卡| 女同久久另类99精品国产91| 日本精品一区二区三区蜜桃| 99久久无色码亚洲精品果冻| 宅男免费午夜| 麻豆国产av国片精品| 久久精品亚洲精品国产色婷小说| 哪里可以看免费的av片| 51午夜福利影视在线观看| 一级片免费观看大全| 一本精品99久久精品77| 久久婷婷成人综合色麻豆| 国产亚洲av嫩草精品影院| 国产精品永久免费网站| av天堂在线播放| 99久久无色码亚洲精品果冻| 久久精品国产综合久久久| 亚洲精品色激情综合| 精品无人区乱码1区二区| 亚洲午夜精品一区,二区,三区| 搡老岳熟女国产| 一区二区三区高清视频在线| 999久久久国产精品视频| 欧美zozozo另类| 国产麻豆成人av免费视频| 亚洲av第一区精品v没综合| 无人区码免费观看不卡| 中文字幕久久专区| 欧美日韩黄片免| 国产一区二区激情短视频| xxxwww97欧美| 怎么达到女性高潮| 久久这里只有精品19| 精品久久久久久久久久久久久 | 日韩高清综合在线| 日本三级黄在线观看| 欧美三级亚洲精品| 国内揄拍国产精品人妻在线 | 欧美丝袜亚洲另类 | 女性被躁到高潮视频| 美女大奶头视频| 69av精品久久久久久| 午夜久久久在线观看| 侵犯人妻中文字幕一二三四区| 99精品欧美一区二区三区四区| www国产在线视频色| 侵犯人妻中文字幕一二三四区| 啦啦啦 在线观看视频| 99热这里只有精品一区 | 国产高清激情床上av| 国产极品粉嫩免费观看在线| 久久久久久久久中文| 欧美激情 高清一区二区三区| 日本撒尿小便嘘嘘汇集6| 精品不卡国产一区二区三区| 91成年电影在线观看| 男女那种视频在线观看| 在线观看一区二区三区| 免费在线观看成人毛片| 夜夜看夜夜爽夜夜摸| 久久久久久久精品吃奶| netflix在线观看网站| 在线观看一区二区三区| 一本一本综合久久| 亚洲国产欧美网| 最好的美女福利视频网| 午夜亚洲福利在线播放| 制服诱惑二区| 国产精品野战在线观看| 国产精品亚洲一级av第二区| 十八禁网站免费在线| 中亚洲国语对白在线视频| 非洲黑人性xxxx精品又粗又长| 中文在线观看免费www的网站 | 亚洲自拍偷在线| 免费高清在线观看日韩| 国产熟女午夜一区二区三区| 精品日产1卡2卡| 可以在线观看毛片的网站| 国产亚洲精品久久久久久毛片| 久久精品国产清高在天天线| 99国产综合亚洲精品| 欧美乱妇无乱码| 在线播放国产精品三级| 中文字幕精品免费在线观看视频| 少妇粗大呻吟视频| 国产成年人精品一区二区| 两个人视频免费观看高清| 欧美色视频一区免费| 国产成人av激情在线播放| 国语自产精品视频在线第100页| 99久久久亚洲精品蜜臀av| 亚洲第一青青草原| 欧美乱妇无乱码| av电影中文网址| 国产免费av片在线观看野外av| 欧美成人一区二区免费高清观看 | 免费在线观看日本一区| 国产一区二区三区视频了| 深夜精品福利| 很黄的视频免费| 脱女人内裤的视频| 神马国产精品三级电影在线观看 | 国产激情欧美一区二区| 搡老岳熟女国产| 母亲3免费完整高清在线观看| 亚洲avbb在线观看| 欧美乱色亚洲激情| 国产亚洲欧美在线一区二区| 欧美日韩一级在线毛片| 9191精品国产免费久久| 欧美激情 高清一区二区三区| 老司机午夜福利在线观看视频| 国语自产精品视频在线第100页| 亚洲国产欧洲综合997久久, | 国产精品 欧美亚洲| 激情在线观看视频在线高清| 国产亚洲欧美在线一区二区| 黄片大片在线免费观看| 国产成+人综合+亚洲专区| 国产精品久久久人人做人人爽| 可以免费在线观看a视频的电影网站| 亚洲成人久久性| 成人午夜高清在线视频 | 亚洲五月色婷婷综合| 亚洲五月天丁香| 午夜老司机福利片| 亚洲色图av天堂| 免费看美女性在线毛片视频| 91麻豆av在线| 777久久人妻少妇嫩草av网站| 欧美激情极品国产一区二区三区| 婷婷丁香在线五月| 欧美日韩黄片免| 久久九九热精品免费| 中文字幕高清在线视频| 国产精品二区激情视频| 在线av久久热| 亚洲成人免费电影在线观看| 免费看美女性在线毛片视频| 好看av亚洲va欧美ⅴa在| 91麻豆精品激情在线观看国产| 久久久国产成人免费| 美国免费a级毛片| 欧美性猛交黑人性爽| 婷婷丁香在线五月| 男人舔女人的私密视频| 久久香蕉激情| 国产亚洲精品av在线| 一级a爱片免费观看的视频| 午夜福利在线观看吧| 中出人妻视频一区二区| 又紧又爽又黄一区二区| 嫩草影院精品99| 91大片在线观看| 精品电影一区二区在线| 欧美成狂野欧美在线观看| 不卡av一区二区三区| 久久人妻福利社区极品人妻图片| 一进一出抽搐动态| 精品高清国产在线一区|