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

    明渠植物對(duì)沙波形態(tài)變化的影響研究

    2023-08-27 15:41:53趙杰孫雪嵐郝瑞霞姚歌張朝瑜黃典晨
    人民黃河 2023年7期
    關(guān)鍵詞:流場(chǎng)植物

    趙杰 孫雪嵐 郝瑞霞 姚歌 張朝瑜 黃典晨

    摘 要:沙波在自然河道中廣泛存在,其對(duì)河床形態(tài)的變化起著不可忽視的作用。通過(guò)在玻璃水槽中采用有機(jī)玻璃棒模擬剛性挺水植物對(duì)沙波各種運(yùn)動(dòng)形態(tài)變化的影響,并引入MIKE21 水動(dòng)力數(shù)值模型,從宏觀與微觀兩個(gè)方面,研究沙波的波長(zhǎng)、波寬、波高、波長(zhǎng)波寬比、迎流面與背流面波長(zhǎng)比等在水中植物影響下的變化現(xiàn)象,探究河道泥沙的群體運(yùn)動(dòng)機(jī)制,為自然河道中泥沙的治理提供一定的借鑒。結(jié)果表明,在挺水植株對(duì)水流流態(tài)變化的影響下,沙波的波長(zhǎng)、波寬、波高沿程形態(tài)變化顯著,尤其在植株繞流影響下,泥沙顆粒在植株兩側(cè)運(yùn)動(dòng)速度增大,呈現(xiàn)沖刷的態(tài)勢(shì),沙波也隨之出現(xiàn)復(fù)雜的運(yùn)動(dòng)形態(tài),植樹(shù)前后上游沙波形態(tài)亦有顯著差異。

    關(guān)鍵詞:沙波;植物;流場(chǎng);MIKE21

    中圖分類號(hào):TV83 文獻(xiàn)標(biāo)志碼:A doi:10.3969/ j.issn.1000-1379.2023.07.008

    引用格式:趙杰,孫雪嵐,郝瑞霞,等.明渠植物對(duì)沙波形態(tài)變化的影響研究[J].人民黃河,2023,45(7):42-47,67.

    隨著我國(guó)國(guó)民經(jīng)濟(jì)的不斷發(fā)展,河流在人民日常生活和國(guó)家經(jīng)濟(jì)建設(shè)中扮演著越來(lái)越重要的角色,尤其是在居民生活供水與農(nóng)業(yè)灌溉用水方面,更是起著不可忽視的作用。系統(tǒng)研究自然河道中推移質(zhì)泥沙的運(yùn)動(dòng)狀況,準(zhǔn)確掌握河床沖淤演變規(guī)律是充分發(fā)揮河流綜合效益的前提條件。沙波形態(tài)變化是水動(dòng)力、泥沙性質(zhì)、各種邊界條件共同作用的結(jié)果[1] ,是河床形態(tài)變化的重要表現(xiàn)形式[2] 。沙波的動(dòng)態(tài)變化經(jīng)常影響著人類的各類活動(dòng),例如航運(yùn)和近岸工程建設(shè)等,尤其是對(duì)淺水區(qū)的航運(yùn)影響顯著;同時(shí)對(duì)河底設(shè)施構(gòu)成安全隱患,如取水工程等。因此,研究沙波形態(tài)變化對(duì)水上航運(yùn)、設(shè)施建設(shè)等均有重要的意義。

    隨著水流的運(yùn)動(dòng),推移質(zhì)顆粒群體在宏觀上做各種不同形式的運(yùn)動(dòng),床面形態(tài)也隨著輸沙強(qiáng)度的不同而實(shí)時(shí)改變[3] 。沙波的形成發(fā)展影響著河床阻力及推移質(zhì)輸沙率等,對(duì)河道的自然演變有著重要的影響[4] 。關(guān)于沙波的形成機(jī)制,Cornish[5] 認(rèn)為初始沙波的形成是泥沙沖淤的結(jié)果,其中沙波橫向發(fā)展速度主要受近床水流影響。Damgaard 等[6] 分析了泥沙非均勻度對(duì)沙波形態(tài)的影響,認(rèn)為沙波波長(zhǎng)、波高隨水流強(qiáng)度的變化不受泥沙非均勻度的影響。Werner 等[7] 歸納了沙波波長(zhǎng)隨時(shí)間的變化規(guī)律。沙波運(yùn)動(dòng)規(guī)律與多種因素存在著密切關(guān)系。Coleman 等[8] 基于水槽試驗(yàn)數(shù)據(jù),建立了沙波波速、波高與推移質(zhì)輸沙率之間的函數(shù)關(guān)系。Raudkivi[9] 通過(guò)數(shù)值試驗(yàn)得出沙波是泥沙顆粒在各種因素共同作用下的宏觀表現(xiàn),并進(jìn)一步發(fā)現(xiàn)沙質(zhì)推移質(zhì)輸沙率與沙波的波高成負(fù)相關(guān)關(guān)系,河道中植物的阻力作用直接、間接影響著沙波的形態(tài)變化,尤其是對(duì)水流流態(tài)的影響,間接影響了沙波的運(yùn)動(dòng)變化規(guī)律。Dunn 等[10] 在矩形渠道中對(duì)淹沒(méi)的剛性與柔性圓棒進(jìn)行試驗(yàn)研究,發(fā)現(xiàn)水流阻力隨植被密度的增大而增大。張?jiān)h等[11] 利用床面形態(tài)控制數(shù)及床面形態(tài)判別方法,研究發(fā)現(xiàn)當(dāng)流量逐漸增大時(shí),床面形態(tài)逐漸由沙壟向過(guò)渡、動(dòng)平整方向發(fā)展。

    在天然河道中,外界環(huán)境條件復(fù)雜,不能保證沙波形態(tài)是在觀測(cè)到的水流條件下形成的,故研究采用室內(nèi)玻璃水槽試驗(yàn)方法,以基本確保沙波形態(tài)與當(dāng)前水流條件的一致性。在植物對(duì)沙波運(yùn)動(dòng)形態(tài)的影響領(lǐng)域,前人的研究甚少,本文基于推移質(zhì)模型沙運(yùn)動(dòng)的試驗(yàn)數(shù)據(jù),分析了植物種植前后,水槽沿程的沙波波高、波長(zhǎng)、波寬等的變化情況,并應(yīng)用MIKE21 軟件對(duì)植被影響下的水流流態(tài)進(jìn)行了數(shù)值模擬計(jì)算,分別從微觀與宏觀兩方面定性、定量地研究剛性植物對(duì)沙波形態(tài)變化的影響及其變化規(guī)律,進(jìn)而比較分析植物種植前后沙波形態(tài)的差異,以期深入理解河床推移質(zhì)運(yùn)動(dòng)的特點(diǎn)與規(guī)律,完善動(dòng)床河道的泥沙治理理論與方法。

    1 試驗(yàn)簡(jiǎn)介

    1.1 試驗(yàn)設(shè)備

    本試驗(yàn)是在太原理工大學(xué)水利科學(xué)與工程學(xué)院水工實(shí)驗(yàn)室的玻璃水槽內(nèi)進(jìn)行的(設(shè)備布置見(jiàn)圖1),水槽長(zhǎng)12 m、寬0.4 m、高0.5 m,玻璃水槽底部平整度小于0.5 mm。水槽模型試驗(yàn)用水由循環(huán)供水系統(tǒng)供給,能保證水流狀態(tài)為恒定流。整個(gè)循環(huán)供水系統(tǒng)由變頻控制水泵提供所需流量,電磁流量計(jì)觀測(cè)實(shí)時(shí)流量。槽首設(shè)有消能柵板,用于穩(wěn)定水流;水槽末端設(shè)有可調(diào)節(jié)的活頁(yè)式尾門(mén),垂直于出水口處設(shè)有一攔沙篩,水流經(jīng)尾門(mén)后穿過(guò)攔沙篩跌入水道。

    水槽上設(shè)有可活動(dòng)測(cè)針架,水槽頂部有帶刻度的軌道,可供測(cè)針架移動(dòng)與讀取所測(cè)斷面的位置。在玻璃水槽底部設(shè)有自動(dòng)升降螺桿,可用于調(diào)節(jié)水槽的坡度。沿水槽縱向自進(jìn)口起設(shè)置了10 個(gè)量測(cè)斷面,斷面間距為1.0 m,水槽軌道刻度零點(diǎn)所在的斷面記為0 斷面,其余斷面依次編號(hào)(見(jiàn)圖1)。

    1.2 試驗(yàn)方案

    試驗(yàn)采用中值粒徑d50 =0.34 mm 的非均勻塑料沙作為推移質(zhì)模型沙,塑料沙在目前河工模型試驗(yàn)中運(yùn)用較為廣泛,密度約為1 050 kg/ m3,呈白色半透明狀。塑料沙為憎水型輕質(zhì)沙,在水中的物理、力學(xué)性質(zhì)都較為穩(wěn)定,且顆粒間一般不存在黏結(jié)現(xiàn)象。水槽底部鋪約5 cm 厚的模型沙,試驗(yàn)過(guò)程中根據(jù)輸沙率在水槽首部加沙。

    為了較真實(shí)地模擬剛性挺水植物對(duì)水沙運(yùn)動(dòng)的影響[12] ,同時(shí)考慮到樹(shù)木的阻力相似,本試驗(yàn)采用有機(jī)玻璃棒模擬自然莖稈狀的挺水植物(見(jiàn)圖2),有機(jī)玻璃棒高26 cm、直徑4 mm。試驗(yàn)采用網(wǎng)格結(jié)構(gòu)來(lái)固定有機(jī)玻璃(模型樹(shù)),試驗(yàn)時(shí)將模型樹(shù)固定在網(wǎng)格的圓形孔內(nèi)。模型樹(shù)布置采用行列交錯(cuò)方式,其中行距為9 cm、株距為8 cm,共7 排,植物群長(zhǎng)60 cm、寬24cm,共25 棵,種植密度為174 棵/ m2,模型樹(shù)布置示意如圖3 所示。為避免平整床面因初始流量過(guò)大沖刷泥沙而破壞底床,試驗(yàn)中水槽流量逐步調(diào)節(jié),最終穩(wěn)定流量為7.69 L/ s,采用沙莫夫公式計(jì)算得泥沙顆粒起動(dòng)流速為0.041 m/ s。

    1.2.2 試驗(yàn)測(cè)試手段

    本試驗(yàn)通過(guò)可活動(dòng)式測(cè)針測(cè)讀水位、沙波各部位的水深等,測(cè)量精度達(dá)0.1 mm。采用多普勒流速儀測(cè)量流速,基于聲學(xué)多普勒原理測(cè)量水體中粒子的運(yùn)動(dòng)速度來(lái)代表水體運(yùn)動(dòng)速度,測(cè)量精度達(dá)0.001 m/ s。水槽流量通過(guò)電磁流量計(jì)測(cè)量,測(cè)量原理是基于法拉第電磁感應(yīng)定律,測(cè)量精度可達(dá)0.01 m3 / h。

    2 試驗(yàn)現(xiàn)象及分析

    2.1 植樹(shù)前沙波沿程變化

    大部分試驗(yàn)數(shù)據(jù)表明,當(dāng)河槽底部比較平坦時(shí),沙波的波高、波長(zhǎng)等都比較均勻。本試驗(yàn)選取了多組形態(tài)較為明顯的沙波進(jìn)行量測(cè),如圖5 所示。受玻璃水槽邊界的影響,部分沙波失真,所以數(shù)據(jù)較為離散,本文選取數(shù)目較多的沙波來(lái)分析其變化的總趨勢(shì)。植樹(shù)前沙波各種形態(tài)沿程變化規(guī)律見(jiàn)圖6,其中紅線為擬合的整體趨勢(shì)線。

    從圖6 中趨勢(shì)線可以看出沙波的波長(zhǎng)、波寬均沿程逐漸增大,主要原因是玻璃水槽中水流流速沿程逐漸加快。相比之下,波長(zhǎng)波寬比卻沿程逐漸減小,說(shuō)明波長(zhǎng)沿程變化幅度大于波寬的變化幅度。迎流面與背流面波長(zhǎng)比也在逐漸增大,說(shuō)明迎流面波長(zhǎng)增長(zhǎng)幅度大于背流面波長(zhǎng)的增長(zhǎng)幅度。試驗(yàn)中波高沿程下降,其趨勢(shì)線斜率比其他指標(biāo)趨勢(shì)線斜率大,直到水槽2.35 m刻度處,河床仍無(wú)明顯起伏狀態(tài)。大部分試驗(yàn)表明,沙波的波高隨水流流速的加快而增大,當(dāng)超過(guò)某一臨界值后,由于水流流速過(guò)快,沖刷加劇,沙波波高反而逐漸減小,直至消亡,說(shuō)明試驗(yàn)中雖然水流流速在水槽中沿程加快,但是水流初始速度已經(jīng)超過(guò)波高增加的臨界值,故研究中植樹(shù)前的沙波波高呈現(xiàn)逐漸減小的趨勢(shì)。從圖6 中可以看出沙波的各種形態(tài)均隨水槽長(zhǎng)度的增加而逐漸變化,直到約2.5 m 處,沙波形態(tài)逐漸消亡。從上述分析可以看出,在玻璃水槽中未植樹(shù)前,沙波形態(tài)沿程變化的趨勢(shì)均具有一定的規(guī)律性,且基本為較單一的線性變化。

    2.2 植樹(shù)后沙波形態(tài)沿程變化

    水槽中植樹(shù)后沙波形態(tài)的沿程變化有較明顯的規(guī)律性,如圖7 所示虛線為擬合的整體變化趨勢(shì)線。圖7(a)為沙波波高沿程變化趨勢(shì)圖,沙波波高在植株群上游呈現(xiàn)先增大后逐漸減小的趨勢(shì),而流速是沿程逐漸遞減的。這是因?yàn)樯巢úǜ唠S流速的增大而增大,但是當(dāng)流速超過(guò)一定值時(shí),水流的沖刷作用又會(huì)使沙波波高逐漸變小,所以在流速逐漸變小過(guò)程中,靠近進(jìn)口段水流流速較大,對(duì)波高有沖刷作用,而流速沿程減小后沖刷作用減弱,推移質(zhì)泥沙顆粒被水流帶到波峰較多,波高高度增大,之后水流流速繼續(xù)變小,泥沙顆粒運(yùn)動(dòng)變緩,波峰高度逐漸減小。圖7(b)為沙波波寬沿程變化趨勢(shì)圖,從趨勢(shì)線可以看出沙波波寬沿程變化幅度較小,但仍然呈現(xiàn)逐漸減小的趨勢(shì),水流流動(dòng)變緩,水流能量減弱,沙波波寬變小。圖7(c)為沙波波長(zhǎng)沿程變化趨勢(shì)圖,從圖中的變化趨勢(shì)線可以看出在植株群上游波長(zhǎng)整體趨勢(shì)變化不大,但有輕微減小,這是因?yàn)橹仓耆簩?duì)水流有一定的阻滯作用,在植株群上游出現(xiàn)壅水,水流流速沿程逐漸減小,沙波波長(zhǎng)隨流速的減小亦呈現(xiàn)逐漸減小的趨勢(shì),由于波峰處水流流速減小,波峰泥沙運(yùn)動(dòng)速度變慢,部分泥沙沉降,宏觀上表現(xiàn)為波峰運(yùn)動(dòng)被阻滯,速度變慢,沙波背流面向前推進(jìn)亦變緩,沙波迎流面流速大于波峰處流速,泥沙運(yùn)動(dòng)速度也同樣大于波峰處泥沙的運(yùn)動(dòng)速度,故波長(zhǎng)在宏觀上表現(xiàn)為變短。相反,波峰處流速大于迎流面的流速時(shí),波長(zhǎng)在宏觀上表現(xiàn)為變長(zhǎng)。圖7(d)為迎流面與背流面波長(zhǎng)之比沿程變化情況,從趨勢(shì)線可以看出,在植株群上游波長(zhǎng)比沿程呈減小趨勢(shì),至植株群前變?yōu)樽钚?,即背流面波長(zhǎng)變化幅度小于迎流面,說(shuō)明背流面淤積程度高于迎流面的沖刷程度。圖7(e)為波長(zhǎng)波寬比沿程變化趨勢(shì)圖,趨勢(shì)線表明波長(zhǎng)波寬比在植株群上游先變大后逐漸減小,與沙波波高沿程變化趨勢(shì)類似。波長(zhǎng)波寬比沿程變化亦有一個(gè)臨界值。波長(zhǎng)、波寬隨流速增大逐漸變大,說(shuō)明同一沙波的波長(zhǎng)沿程變化幅度先大于波寬的變化幅度,之后又逐漸小于波寬的變化幅度,除去一些特殊點(diǎn),波長(zhǎng)波寬比始終大于1,即波長(zhǎng)始終大于等于波寬。從圖7(e)中可以看出在植株群下游的沙波數(shù)量突然變少,這是因?yàn)橹仓耆合掠嗡畲蠓葴p小,流速大幅度增大,沙波形態(tài)逐漸消亡,一段距離后水深逐漸恢復(fù)至植樹(shù)前的水深,流速亦逐漸減小,河床底部再次產(chǎn)生沙波形態(tài)。

    2.3 基于MIKE21 的有機(jī)玻璃棒繞流運(yùn)動(dòng)分析

    物理試驗(yàn)對(duì)水槽流場(chǎng)的觀測(cè)難以保證連續(xù)性,且觀測(cè)結(jié)果受限于各種外界條件的影響,無(wú)法得到一個(gè)較準(zhǔn)確完整的水槽流速分布及變化,故此處基于MIKE21 軟件對(duì)水槽水流的流場(chǎng)變化進(jìn)行水動(dòng)力數(shù)值模擬,可以得到更加直觀的流場(chǎng)分布,更全面地分析模型樹(shù)周圍的流場(chǎng)變化,為進(jìn)一步研究沙波的模擬計(jì)算做基礎(chǔ)。MIKE21 是一個(gè)專業(yè)的工程軟件包,在平面二維自由表面流數(shù)值模擬方面具有強(qiáng)大的功能。為了更全面直觀地研究水槽中的流場(chǎng)變化, 本文基于MIKE21 水動(dòng)力模塊(HD)對(duì)有機(jī)玻璃棒的繞流運(yùn)動(dòng)形態(tài)進(jìn)行仿真模擬,從微觀角度更加細(xì)致地去觀察水流流態(tài)在圓棒阻力作用下的完整變化。本模擬試驗(yàn)采用圓柱體模擬剛性挺水植株,并分別設(shè)計(jì)了單圓柱體、并排雙圓柱體、3 個(gè)圓柱體3 種水流繞流運(yùn)動(dòng)的數(shù)值模擬方案(模擬結(jié)果如圖8 所示,不同顏色及其顏色的深淺程度代表不同的流速范圍)。

    結(jié)合現(xiàn)有的研究成果,從流體力學(xué)角度分析,每一棵模型樹(shù)實(shí)際上就是一個(gè)邊界條件極為復(fù)雜的柱體。圓柱體形體阻力、阻水范圍與拖曳力的存在,使有機(jī)玻璃棒迎流面一段范圍的流速降低,從圖8(a)中可以看出,有機(jī)玻璃棒周圍流速分布為層式分布,類似蠟燭火焰的分布形式。最內(nèi)側(cè)即緊貼圓柱體一層為藍(lán)色層,距離圓柱體越近、顏色越深、速度越小,其中顏色最深處流速小于0.096 m/ s。整個(gè)圓柱體周圍的流速分布呈拖曳狀。在順?biāo)鞣较?,圓柱體兩側(cè)各有一半圓形紅色流速區(qū),流速大于等于0.176 m/ s。結(jié)合已有試驗(yàn)結(jié)果可知,水流在圓柱體的前駐點(diǎn)速度為0,水流沿圓柱體兩側(cè)流動(dòng),水流壓力在圓柱體前半部分慢慢減小,流速逐漸增大到最大值,故圓柱體兩側(cè)有紅色流速區(qū),并大于來(lái)水流速。隨著水流向圓柱體后半部分流動(dòng),水流壓力慢慢增大,流速逐漸減小,當(dāng)水流到達(dá)圓柱體背流面駐點(diǎn)時(shí),流速重新變?yōu)椋啊?/p>

    圓柱體周圍泥沙顆粒運(yùn)動(dòng)隨流態(tài)變化而變化,泥沙顆粒在圓柱體周圍有部分沉積,但由于流速減少范圍較小,且考慮到有機(jī)玻璃棒在水中來(lái)回振動(dòng)作用,實(shí)際中植株群前流速并不為0,植株群前、后流速差異并不十分顯著,因此并未觀察到明顯的泥沙沉積現(xiàn)象。由于模型樹(shù)兩側(cè)的流速增大,因此迎流面波長(zhǎng)逐漸增大。沙波運(yùn)動(dòng)到模型樹(shù)前時(shí),受樹(shù)形阻力作用與水流作用,波寬增長(zhǎng)趨勢(shì)明顯。圖7(e)中顯示波長(zhǎng)波寬比在植株群內(nèi)逐漸減小,這是因?yàn)椴▽捵兓却笥诓ㄩL(zhǎng)的變化幅度,但植株群后的沖刷使波長(zhǎng)波寬比又突然增大。圓柱體兩側(cè)流速加快,泥沙顆粒起動(dòng)數(shù)量變多,運(yùn)動(dòng)速度加快,有沖刷趨勢(shì),故沙波波高有降低趨勢(shì)。

    當(dāng)2 個(gè)圓柱體并排同時(shí)存在于水流中時(shí),如圖8(b)所示,有2 個(gè)拖曳狀流速分布區(qū)并行存在,2 個(gè)圓柱體中間流速增大的范圍變大,但流速增大幅度變小,表明由于兩邊水流的對(duì)沖,消耗了一部分動(dòng)能,流速增大的速度變緩,但速度增大范圍較單圓柱體整體變大,對(duì)沙波影響范圍亦增大。

    當(dāng)3 個(gè)圓柱體同時(shí)存在于水流中時(shí),如圖8(c)所示,由于第二排圓柱體的阻力作用,第一排平行放置的2 個(gè)圓柱體間流速增大范圍較圖8(b)變小。同時(shí),受第一排兩個(gè)圓柱體背流面流速變化的影響,第二排圓柱體兩側(cè)流速增大范圍亦變小。此時(shí)沙波變化較復(fù)雜,整體波高呈減小趨勢(shì),波寬迅速增大。

    從圖7 中可看出,植株群范圍內(nèi)波長(zhǎng)逐漸增加,迎流面與背流面波長(zhǎng)比逐漸增大,相比較波長(zhǎng)波寬比的變化不是很顯著,并有減小的趨勢(shì)。結(jié)合圖8(c)流速分布圖,從流速解釋該變化現(xiàn)象,前后排圓柱體周圍流速整體偏大,故水流可帶起更多泥沙顆粒,整體偏沖刷態(tài)勢(shì)。模型樹(shù)形體阻力等的存在,使波寬、波長(zhǎng)增大的同時(shí)波高逐漸減小,但由于本試驗(yàn)中模型樹(shù)密度較小,因此沙波并未消失。

    考慮到有機(jī)玻璃棒在水中來(lái)回振動(dòng),植株群前流速并不為0,植株群前、后流速差異并不顯著,因此圓柱繞流形成的復(fù)雜條件導(dǎo)致沙波波峰高度并不十分顯著。但是模型樹(shù)周圍的流速增大,導(dǎo)致迎流面波長(zhǎng)逐漸增大。沙波運(yùn)動(dòng)到模型樹(shù)前時(shí),受模型樹(shù)的樹(shù)形阻力作用與水流作用,波寬增長(zhǎng)速度變快。圖7(e)中顯示波長(zhǎng)波寬比在植株群內(nèi)逐漸減小,這是因?yàn)椴▽捵兓却笥诓ㄩL(zhǎng)變化幅度,但植株群后的沖刷使波長(zhǎng)波寬比又突然增大。

    2.4 植樹(shù)前后對(duì)比分析

    采用spss 統(tǒng)計(jì)軟件對(duì)玻璃水槽中植樹(shù)前后0 ~3 m之間的沙波形態(tài)數(shù)據(jù)進(jìn)行獨(dú)立樣本t 檢驗(yàn)分析,分別選取植樹(shù)前45 組沙波與植樹(shù)后56 組沙波,判斷植樹(shù)條件變化對(duì)上游沙波形態(tài)變化的影響程度,結(jié)果見(jiàn)表1、表2。

    表1 是植樹(shù)前后沙波波高、波寬、波長(zhǎng)、迎流面與背流面波長(zhǎng)比、波長(zhǎng)波寬比試驗(yàn)結(jié)果統(tǒng)計(jì)??梢钥闯觯矘?shù)前后沙波數(shù)目相差不多,即植樹(shù)對(duì)上游沙波數(shù)目影響并不顯著。波寬、波長(zhǎng)的平均值植樹(shù)前均略大于植樹(shù)后;與之相反,波高與波長(zhǎng)波寬比植樹(shù)后均顯著大于植樹(shù)前,迎流面與背流面波長(zhǎng)比則基本相同。

    從標(biāo)準(zhǔn)偏差來(lái)看,波高、波寬、迎流面與背流面波長(zhǎng)比、波長(zhǎng)波寬比均是植樹(shù)后顯著大于植樹(shù)前,說(shuō)明植樹(shù)后數(shù)據(jù)變得更加紊亂,即模型樹(shù)的存在使沙波多種形態(tài)數(shù)據(jù)的離散程度變大。但是波長(zhǎng)的標(biāo)準(zhǔn)偏差在植樹(shù)后略大于植樹(shù)前,說(shuō)明植樹(shù)后波長(zhǎng)數(shù)據(jù)離散程度變化不大。

    波高、波寬、波長(zhǎng)的萊文方差等同性檢驗(yàn)的顯著性均大于0.05,即判斷植樹(shù)前后兩總體為方差齊性,選擇假定等方差t 檢驗(yàn),迎流面與背流面波長(zhǎng)比、波長(zhǎng)波寬比的萊文方差等同性檢驗(yàn)的顯著性均小于0.05,即判斷植樹(shù)前后兩總體為方差非齊性,選擇服從修正自由度的不假定等方差t 檢驗(yàn)。

    波高、波長(zhǎng)波寬比的Sig.(雙尾)均小于0.01,說(shuō)明這兩種沙波形態(tài)在植樹(shù)前后有顯著差異。波寬、波長(zhǎng)、迎流面與背流面波長(zhǎng)比的Sig.(雙尾)均大于0.01,說(shuō)明植樹(shù)前后無(wú)顯著差異。

    從圖6 可以看出,植樹(shù)前沙波集中出現(xiàn)在0 ~3.0 m范圍內(nèi),而植樹(shù)后則沿程均出現(xiàn)沙波形態(tài),說(shuō)明植樹(shù)對(duì)河床沙波影響顯著。在自然河流內(nèi)種植植物后,應(yīng)該重視河床沙波形態(tài)的變化,進(jìn)而對(duì)河道自然條件下河床沖淤演變規(guī)律進(jìn)行趨勢(shì)預(yù)測(cè),對(duì)于合理安排河流兩岸經(jīng)濟(jì)技術(shù)活動(dòng),保證擬建設(shè)工程安全并發(fā)揮應(yīng)有的效益至關(guān)重要。

    結(jié)合表1 的數(shù)據(jù),可判斷出植樹(shù)對(duì)沙波形態(tài)的影響。植樹(shù)條件的變化對(duì)沙波波寬、波長(zhǎng)、迎流面與背流面波長(zhǎng)比無(wú)顯著影響。但是,上游波高與波長(zhǎng)波寬比受植樹(shù)條件影響較大,植樹(shù)后波高與波長(zhǎng)波寬比顯著增加,沙波輪廓更加明顯,河床形態(tài)更加復(fù)雜。

    3 結(jié)論

    1)沙波的各項(xiàng)要素受植樹(shù)情況影響較大,沿程變化具有一定的規(guī)律性,從變化趨勢(shì)上看各有區(qū)別。

    2)MIKE21 水動(dòng)力模型的模擬結(jié)果表明,在植株群范圍內(nèi),水流條件復(fù)雜,沙波受植株形體阻力與阻力范圍影響顯著。

    3)獨(dú)立樣本t 檢驗(yàn)表明,植樹(shù)后波高、波長(zhǎng)波寬比顯著增大,沙波輪廓更加明顯。

    4)泥沙群體運(yùn)動(dòng)是一個(gè)非常復(fù)雜的問(wèn)題,目前研究還很不充分,尤其是放入植株群后使問(wèn)題更加復(fù)雜。進(jìn)行更深入的研究,除需要用先進(jìn)的測(cè)量手段、儀器和大量的模型試驗(yàn)外,還需要合理選擇參數(shù)并積累豐富可靠的實(shí)測(cè)資料。

    5)本文未對(duì)沙波形態(tài)變化進(jìn)行數(shù)值模擬分析,今后將進(jìn)一步研究,結(jié)合更多的試驗(yàn)數(shù)據(jù),總結(jié)沙波變化規(guī)律,得出相應(yīng)的計(jì)算公式和適用的計(jì)算方法,結(jié)合計(jì)算機(jī)技術(shù)解決實(shí)際工程中的泥沙問(wèn)題。

    參考文獻(xiàn):

    [1] 錢(qián)寧,萬(wàn)兆惠.泥沙運(yùn)動(dòng)力學(xué)[>M].北京:科學(xué)出版社,1983:140-182.

    [2] 段文忠.沙波尺度和運(yùn)動(dòng)速度與水力泥沙因素的關(guān)系[J].武漢水利電力學(xué)院學(xué)報(bào),1988,21(4):66-72.

    [3] 惠二青,江春波,潘應(yīng)旺.植被覆蓋的河道水流縱向流速垂向分布[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,49(6):834-837.

    [4] 高昂.沙波面水沙運(yùn)動(dòng)特性及沙波演化規(guī)律試驗(yàn)研究[D].鄭州:華北水利水電大學(xué),2016:23-48.

    [5] CORNISH V.On Sand?Waves in Tidal Currents[J].The Geo?graphical Joural,1991,18(2):170-200.

    [6] DAMGAARD J S,SOULSBY R,PEET A,et al.Sand Transporton Steeply Sloping Plane and Rippled Bed[J]. Journal ofHydrologic Engineering,2003,129(9):706-719.

    [7] WERNER B,KOCUREK G.Bedform Spacing from Defect Dy?namics[J].Geology,1999,27(8):727-730.

    [8] COLEMAN S E,MELVILLE B W.Bed?Form Development[J].Journal of Hydraulic Engineering,1994,120(4):544-560.

    [9] RAUDKIVI A J.Rippleson Stream Bed [J].Journal of Hy?Draulic Engineering,1997,123(1):58-64.

    [10] DUNN C,LOPEZ F,GARCIA M.Mean Flow and Turbulencein a Laboratory Channel with Simulated Vegetation[R].Chi?cago:Dept of Civil Engrg,Univ of Illinois at Urbana?Cham?paign,1996:19-105.

    [11] 張?jiān)h,王平.黃河下游游蕩型河段床面形態(tài)變化特征[J].人民黃河,2018,40(8):8-11.

    [12] 程年生.植被對(duì)明渠水流及泥沙運(yùn)動(dòng)的影響[J].水力發(fā)電學(xué)報(bào),2018,37(9):1-7.

    【責(zé)任編輯 簡(jiǎn) 群】

    猜你喜歡
    流場(chǎng)植物
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場(chǎng)計(jì)算
    轉(zhuǎn)杯紡排雜區(qū)流場(chǎng)與排雜性能
    植物的防身術(shù)
    把植物做成藥
    哦,不怕,不怕
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場(chǎng)統(tǒng)計(jì)分析
    將植物穿身上
    基于瞬態(tài)流場(chǎng)計(jì)算的滑動(dòng)軸承靜平衡位置求解
    植物罷工啦?
    植物也瘋狂
    99re6热这里在线精品视频| 在线看a的网站| 各种免费的搞黄视频| 叶爱在线成人免费视频播放| 午夜免费观看性视频| 在线亚洲精品国产二区图片欧美| 国产一级毛片在线| 妹子高潮喷水视频| 国产黄频视频在线观看| 日韩,欧美,国产一区二区三区| 两个人免费观看高清视频| 久久人人97超碰香蕉20202| 韩国av在线不卡| 18禁动态无遮挡网站| 777久久人妻少妇嫩草av网站| 99久国产av精品国产电影| 国产免费视频播放在线视频| 色视频在线一区二区三区| 欧美少妇被猛烈插入视频| 成年女人毛片免费观看观看9 | 男女边吃奶边做爰视频| 亚洲精品国产一区二区精华液| 午夜老司机福利片| 日本一区二区免费在线视频| 国产成人欧美在线观看 | 日本一区二区免费在线视频| 亚洲欧美成人综合另类久久久| 精品一区二区免费观看| 精品午夜福利在线看| 亚洲免费av在线视频| 一本色道久久久久久精品综合| 极品人妻少妇av视频| 如日韩欧美国产精品一区二区三区| 哪个播放器可以免费观看大片| 桃花免费在线播放| 亚洲国产av新网站| 别揉我奶头~嗯~啊~动态视频 | 国产爽快片一区二区三区| 自线自在国产av| 亚洲国产看品久久| 国产av一区二区精品久久| h视频一区二区三区| tube8黄色片| 少妇精品久久久久久久| 一区二区日韩欧美中文字幕| 精品视频人人做人人爽| 国产乱人偷精品视频| 久久精品国产a三级三级三级| 超碰97精品在线观看| 超碰成人久久| 我要看黄色一级片免费的| 一区二区三区乱码不卡18| 热99国产精品久久久久久7| 女性被躁到高潮视频| 午夜福利影视在线免费观看| 国产免费视频播放在线视频| 色播在线永久视频| 欧美乱码精品一区二区三区| 精品卡一卡二卡四卡免费| 免费观看性生交大片5| 日韩成人av中文字幕在线观看| 人人妻人人添人人爽欧美一区卜| av国产精品久久久久影院| www.熟女人妻精品国产| 激情五月婷婷亚洲| 国产精品女同一区二区软件| 一级毛片电影观看| 国产精品一区二区在线不卡| 综合色丁香网| 好男人视频免费观看在线| 丝袜美足系列| 精品卡一卡二卡四卡免费| 中文字幕精品免费在线观看视频| 下体分泌物呈黄色| 美女国产高潮福利片在线看| 欧美精品亚洲一区二区| 九九爱精品视频在线观看| 国产精品麻豆人妻色哟哟久久| 精品国产一区二区三区四区第35| 国产不卡av网站在线观看| 高清视频免费观看一区二区| 新久久久久国产一级毛片| 人人妻,人人澡人人爽秒播 | 免费观看性生交大片5| 热re99久久国产66热| 久久韩国三级中文字幕| 国产片特级美女逼逼视频| 丝袜脚勾引网站| 丰满饥渴人妻一区二区三| 老熟女久久久| 亚洲欧美成人精品一区二区| av.在线天堂| 精品福利永久在线观看| 亚洲精品久久久久久婷婷小说| 亚洲美女视频黄频| xxx大片免费视频| 一本色道久久久久久精品综合| 亚洲国产中文字幕在线视频| 卡戴珊不雅视频在线播放| av卡一久久| www日本在线高清视频| 曰老女人黄片| 亚洲欧美激情在线| 女的被弄到高潮叫床怎么办| 亚洲三区欧美一区| 国产无遮挡羞羞视频在线观看| 九九爱精品视频在线观看| 人人妻人人爽人人添夜夜欢视频| 最近中文字幕2019免费版| 亚洲精品aⅴ在线观看| 久久性视频一级片| 少妇人妻久久综合中文| 国产精品人妻久久久影院| 桃花免费在线播放| 看非洲黑人一级黄片| 午夜激情久久久久久久| 丰满迷人的少妇在线观看| 午夜91福利影院| 国产精品免费视频内射| 高清在线视频一区二区三区| 国产不卡av网站在线观看| 国产有黄有色有爽视频| 亚洲av日韩在线播放| 伊人亚洲综合成人网| 久久99精品国语久久久| 操美女的视频在线观看| av天堂久久9| 亚洲精品国产av成人精品| 精品少妇一区二区三区视频日本电影 | 欧美精品av麻豆av| 中文字幕高清在线视频| 丰满饥渴人妻一区二区三| 欧美av亚洲av综合av国产av | 9色porny在线观看| 老司机影院成人| 女的被弄到高潮叫床怎么办| 秋霞伦理黄片| 亚洲成色77777| 精品免费久久久久久久清纯 | 国产精品香港三级国产av潘金莲 | a级片在线免费高清观看视频| 欧美日韩综合久久久久久| 男女边摸边吃奶| 成年女人毛片免费观看观看9 | videosex国产| 久久免费观看电影| 最近最新中文字幕免费大全7| 天堂8中文在线网| 黄色视频在线播放观看不卡| 一区福利在线观看| 又大又黄又爽视频免费| 国产老妇伦熟女老妇高清| 亚洲欧美成人精品一区二区| 一本色道久久久久久精品综合| 九色亚洲精品在线播放| 久久女婷五月综合色啪小说| 久久久久国产一级毛片高清牌| 亚洲欧美精品综合一区二区三区| 日韩,欧美,国产一区二区三区| 七月丁香在线播放| 国产精品国产三级国产专区5o| 女性生殖器流出的白浆| 一边亲一边摸免费视频| 男人舔女人的私密视频| 亚洲av日韩在线播放| 青春草亚洲视频在线观看| 久久久久久久国产电影| a级毛片在线看网站| 日本av手机在线免费观看| 国产日韩欧美亚洲二区| 欧美亚洲 丝袜 人妻 在线| 久久韩国三级中文字幕| 亚洲欧美成人精品一区二区| 各种免费的搞黄视频| 麻豆乱淫一区二区| 51午夜福利影视在线观看| 伦理电影免费视频| tube8黄色片| 超色免费av| 国产成人a∨麻豆精品| 精品国产一区二区三区四区第35| 少妇人妻 视频| 9热在线视频观看99| 久久精品国产亚洲av涩爱| 七月丁香在线播放| 欧美亚洲日本最大视频资源| 欧美在线一区亚洲| 观看美女的网站| 女人高潮潮喷娇喘18禁视频| 如日韩欧美国产精品一区二区三区| 人妻人人澡人人爽人人| 久久精品国产亚洲av涩爱| 亚洲精品乱久久久久久| 日韩中文字幕欧美一区二区 | 一级毛片 在线播放| 亚洲五月色婷婷综合| 日韩熟女老妇一区二区性免费视频| 午夜激情久久久久久久| 欧美日韩亚洲高清精品| 久久天躁狠狠躁夜夜2o2o | 久久影院123| 丰满少妇做爰视频| 老司机深夜福利视频在线观看 | 久久99一区二区三区| 亚洲伊人久久精品综合| 亚洲情色 制服丝袜| 好男人视频免费观看在线| 女的被弄到高潮叫床怎么办| 午夜福利一区二区在线看| 久久人妻熟女aⅴ| 亚洲精品久久久久久婷婷小说| 少妇人妻 视频| svipshipincom国产片| 日韩大码丰满熟妇| 黄片无遮挡物在线观看| 国产亚洲午夜精品一区二区久久| 日本av手机在线免费观看| 自线自在国产av| 久久久亚洲精品成人影院| 亚洲 欧美一区二区三区| 一级爰片在线观看| 婷婷色综合大香蕉| 亚洲欧美中文字幕日韩二区| 久久久久久久久久久免费av| 久久久久国产一级毛片高清牌| 亚洲精品中文字幕在线视频| 久久国产精品大桥未久av| www.熟女人妻精品国产| 国产精品久久久久久久久免| 精品福利永久在线观看| 一级毛片我不卡| 欧美精品高潮呻吟av久久| 欧美日韩视频精品一区| 99九九在线精品视频| 国产又色又爽无遮挡免| 午夜日本视频在线| 欧美 亚洲 国产 日韩一| 亚洲美女搞黄在线观看| 51午夜福利影视在线观看| 大片电影免费在线观看免费| 国产日韩欧美亚洲二区| 男的添女的下面高潮视频| 欧美日韩视频高清一区二区三区二| 成人三级做爰电影| 一边摸一边做爽爽视频免费| 亚洲欧美一区二区三区黑人| www.精华液| 伊人久久大香线蕉亚洲五| 色网站视频免费| 精品一区二区三区av网在线观看 | 亚洲国产毛片av蜜桃av| 男女边摸边吃奶| 午夜福利影视在线免费观看| 女人久久www免费人成看片| 最新在线观看一区二区三区 | 久久精品国产综合久久久| 丝袜在线中文字幕| 亚洲av国产av综合av卡| 亚洲精品国产av成人精品| avwww免费| 欧美人与性动交α欧美精品济南到| 在线观看免费高清a一片| 国产一卡二卡三卡精品 | 老司机深夜福利视频在线观看 | 最近手机中文字幕大全| 日本黄色日本黄色录像| 桃花免费在线播放| 国产在线一区二区三区精| 大香蕉久久网| 精品一区在线观看国产| 欧美黑人精品巨大| 国产亚洲欧美精品永久| 国产精品麻豆人妻色哟哟久久| 免费在线观看黄色视频的| 99热全是精品| 亚洲国产欧美一区二区综合| 免费日韩欧美在线观看| 久久久久久久久久久免费av| 久久久久国产一级毛片高清牌| 亚洲专区中文字幕在线 | 国产精品嫩草影院av在线观看| 欧美激情极品国产一区二区三区| 久久久精品国产亚洲av高清涩受| 另类精品久久| 丝袜美腿诱惑在线| 精品少妇内射三级| 国产淫语在线视频| 精品久久蜜臀av无| 777久久人妻少妇嫩草av网站| 国产av码专区亚洲av| 久久久精品国产亚洲av高清涩受| 男人舔女人的私密视频| 国产高清不卡午夜福利| 久久 成人 亚洲| 久久天堂一区二区三区四区| 十分钟在线观看高清视频www| 老司机靠b影院| 欧美最新免费一区二区三区| 一区二区三区四区激情视频| 青春草视频在线免费观看| 精品一区二区三区四区五区乱码 | 亚洲成人av在线免费| 啦啦啦中文免费视频观看日本| 国产在线一区二区三区精| 69精品国产乱码久久久| 2021少妇久久久久久久久久久| 国产精品.久久久| 在线观看三级黄色| 久久久国产欧美日韩av| 亚洲人成电影观看| 免费黄色在线免费观看| 久久精品人人爽人人爽视色| av卡一久久| 午夜影院在线不卡| 欧美精品一区二区免费开放| 一本久久精品| 男人添女人高潮全过程视频| 少妇人妻 视频| 如何舔出高潮| 国产熟女欧美一区二区| 丝袜美足系列| 午夜福利一区二区在线看| 亚洲综合色网址| 青春草国产在线视频| 精品酒店卫生间| 一本色道久久久久久精品综合| 国产黄色视频一区二区在线观看| 色精品久久人妻99蜜桃| 侵犯人妻中文字幕一二三四区| 大话2 男鬼变身卡| 啦啦啦在线免费观看视频4| kizo精华| 一本色道久久久久久精品综合| 亚洲伊人色综图| 欧美人与性动交α欧美精品济南到| 亚洲欧美清纯卡通| 伦理电影大哥的女人| 精品一区二区三区四区五区乱码 | 精品第一国产精品| 久久精品国产a三级三级三级| 操出白浆在线播放| 我要看黄色一级片免费的| 大陆偷拍与自拍| 国产免费视频播放在线视频| 精品人妻在线不人妻| 日韩av免费高清视频| 亚洲成人一二三区av| 久久性视频一级片| 精品一区二区三区av网在线观看 | 大香蕉久久网| 日韩大片免费观看网站| 日日啪夜夜爽| 亚洲精品久久久久久婷婷小说| 少妇人妻精品综合一区二区| 欧美激情高清一区二区三区 | 久久99一区二区三区| 制服人妻中文乱码| av视频免费观看在线观看| 看非洲黑人一级黄片| 欧美日韩亚洲高清精品| 51午夜福利影视在线观看| 夜夜骑夜夜射夜夜干| 国产一区二区激情短视频 | 国语对白做爰xxxⅹ性视频网站| 久久天堂一区二区三区四区| 亚洲欧洲精品一区二区精品久久久 | 亚洲男人天堂网一区| 亚洲国产精品一区三区| 人妻一区二区av| 欧美另类一区| 免费看不卡的av| 免费不卡黄色视频| 国产日韩欧美在线精品| 日韩大片免费观看网站| 嫩草影视91久久| 99热国产这里只有精品6| 色婷婷久久久亚洲欧美| 伊人亚洲综合成人网| 亚洲,一卡二卡三卡| 日韩,欧美,国产一区二区三区| 波多野结衣一区麻豆| 国产精品一区二区精品视频观看| 久久精品人人爽人人爽视色| 一本大道久久a久久精品| 国产不卡av网站在线观看| 久久久久人妻精品一区果冻| 中文字幕亚洲精品专区| 精品久久久久久电影网| 人妻人人澡人人爽人人| 永久免费av网站大全| 五月天丁香电影| 夫妻性生交免费视频一级片| netflix在线观看网站| 尾随美女入室| 欧美日韩视频高清一区二区三区二| 欧美97在线视频| 成人三级做爰电影| 九九爱精品视频在线观看| 亚洲伊人久久精品综合| 免费看av在线观看网站| 又大又爽又粗| 男人爽女人下面视频在线观看| 亚洲一码二码三码区别大吗| 精品少妇内射三级| 精品免费久久久久久久清纯 | 黄色怎么调成土黄色| 妹子高潮喷水视频| 亚洲欧洲国产日韩| 九九爱精品视频在线观看| 国产极品天堂在线| 国产精品久久久久久精品古装| 男女免费视频国产| 一边摸一边抽搐一进一出视频| 亚洲在久久综合| 看免费av毛片| 两个人免费观看高清视频| 搡老乐熟女国产| 99re6热这里在线精品视频| 亚洲国产中文字幕在线视频| 亚洲男人天堂网一区| 成年女人毛片免费观看观看9 | 精品一区在线观看国产| 国产精品久久久人人做人人爽| 久久这里只有精品19| 午夜福利视频在线观看免费| av网站在线播放免费| 熟妇人妻不卡中文字幕| 精品久久久久久电影网| 久久久久网色| 国产成人精品久久二区二区91 | 色综合欧美亚洲国产小说| 久久精品熟女亚洲av麻豆精品| 日韩av免费高清视频| 少妇精品久久久久久久| 国产亚洲午夜精品一区二区久久| 99热全是精品| 亚洲精品aⅴ在线观看| 国产亚洲av高清不卡| 黄频高清免费视频| 国产精品一区二区精品视频观看| 婷婷色麻豆天堂久久| 9热在线视频观看99| av在线观看视频网站免费| 9热在线视频观看99| av在线观看视频网站免费| 国精品久久久久久国模美| 亚洲成色77777| 一区二区三区乱码不卡18| 久久久国产精品麻豆| 国语对白做爰xxxⅹ性视频网站| 精品人妻熟女毛片av久久网站| 大片免费播放器 马上看| 考比视频在线观看| 人人妻人人添人人爽欧美一区卜| 岛国毛片在线播放| 国产一区二区 视频在线| 五月开心婷婷网| 青春草亚洲视频在线观看| 高清视频免费观看一区二区| 一二三四中文在线观看免费高清| 亚洲精华国产精华液的使用体验| 国产成人免费无遮挡视频| 日韩电影二区| 久久人人爽av亚洲精品天堂| 亚洲精品美女久久av网站| 成人影院久久| 亚洲精品美女久久av网站| 亚洲欧洲精品一区二区精品久久久 | 久久久精品区二区三区| 亚洲欧洲国产日韩| 亚洲三区欧美一区| 菩萨蛮人人尽说江南好唐韦庄| 国产色婷婷99| 在线观看三级黄色| 高清视频免费观看一区二区| 国产免费一区二区三区四区乱码| 看十八女毛片水多多多| 2021少妇久久久久久久久久久| 国产免费现黄频在线看| 亚洲欧美中文字幕日韩二区| 国产 一区精品| 天堂俺去俺来也www色官网| 午夜老司机福利片| 国产精品久久久久久精品古装| 街头女战士在线观看网站| 亚洲av在线观看美女高潮| 69精品国产乱码久久久| 搡老岳熟女国产| 国产片内射在线| 午夜福利乱码中文字幕| www.精华液| 久久精品aⅴ一区二区三区四区| 看非洲黑人一级黄片| 日本欧美视频一区| 天美传媒精品一区二区| av视频免费观看在线观看| 国产精品久久久久久人妻精品电影 | 欧美精品人与动牲交sv欧美| 999精品在线视频| 最近最新中文字幕免费大全7| 97在线人人人人妻| 亚洲欧洲精品一区二区精品久久久 | 日日爽夜夜爽网站| 久久久久精品人妻al黑| 中文字幕色久视频| 欧美日韩精品网址| 天天操日日干夜夜撸| 啦啦啦啦在线视频资源| 中文天堂在线官网| 久久久久久久国产电影| 国产男女内射视频| 黄片小视频在线播放| 国产精品一区二区在线不卡| 大陆偷拍与自拍| 国产一区二区 视频在线| 国产成人精品在线电影| 久久久久国产精品人妻一区二区| 亚洲欧洲国产日韩| 一区福利在线观看| 亚洲精品国产色婷婷电影| 久久精品人人爽人人爽视色| 欧美久久黑人一区二区| 中文字幕另类日韩欧美亚洲嫩草| 建设人人有责人人尽责人人享有的| 国产成人精品无人区| 飞空精品影院首页| 国产精品一国产av| 巨乳人妻的诱惑在线观看| 日韩成人av中文字幕在线观看| 99re6热这里在线精品视频| 我要看黄色一级片免费的| 汤姆久久久久久久影院中文字幕| 激情视频va一区二区三区| 国产极品粉嫩免费观看在线| 九九爱精品视频在线观看| 香蕉国产在线看| 青春草国产在线视频| 亚洲精品美女久久av网站| 在线观看三级黄色| 国产又爽黄色视频| bbb黄色大片| 精品人妻一区二区三区麻豆| 一区在线观看完整版| 国产欧美日韩一区二区三区在线| 天天躁夜夜躁狠狠久久av| 久久人人爽人人片av| 肉色欧美久久久久久久蜜桃| 欧美人与性动交α欧美软件| 欧美黑人精品巨大| 视频区图区小说| 一区二区av电影网| 最近最新中文字幕大全免费视频 | 久久精品久久久久久久性| 国产成人欧美在线观看 | 日韩精品有码人妻一区| avwww免费| 国精品久久久久久国模美| 美女大奶头黄色视频| 99久久人妻综合| 在线观看国产h片| 黑丝袜美女国产一区| 一区二区三区激情视频| 捣出白浆h1v1| 免费观看a级毛片全部| 啦啦啦在线免费观看视频4| 丝袜美腿诱惑在线| 欧美日韩国产mv在线观看视频| 久久国产精品大桥未久av| 在现免费观看毛片| 日韩大码丰满熟妇| 丝袜在线中文字幕| 免费不卡黄色视频| 免费少妇av软件| 欧美少妇被猛烈插入视频| videos熟女内射| 国产av码专区亚洲av| 水蜜桃什么品种好| 中文欧美无线码| 如何舔出高潮| 日日撸夜夜添| 日韩欧美一区视频在线观看| 最近中文字幕高清免费大全6| 国产探花极品一区二区| 午夜福利视频精品| 欧美日韩一区二区视频在线观看视频在线| 亚洲综合精品二区| 纵有疾风起免费观看全集完整版| 黄色视频在线播放观看不卡| 国产精品国产三级专区第一集| 18禁动态无遮挡网站| 亚洲免费av在线视频| 日韩一区二区视频免费看| 欧美久久黑人一区二区| 午夜福利免费观看在线| 欧美国产精品va在线观看不卡| 18在线观看网站| 国产精品久久久久久久久免| 色婷婷av一区二区三区视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品乱久久久久久| 波多野结衣av一区二区av| 色综合欧美亚洲国产小说| 黄片小视频在线播放| 午夜日韩欧美国产| 亚洲七黄色美女视频| 麻豆精品久久久久久蜜桃| 美女主播在线视频| 丝袜喷水一区| 午夜免费鲁丝| 男女下面插进去视频免费观看| 亚洲美女搞黄在线观看| 天堂8中文在线网| 天堂中文最新版在线下载| 九九爱精品视频在线观看| 最新在线观看一区二区三区 |