• <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)軸承靜平衡位置求解
    植物罷工啦?
    植物也瘋狂
    久久精品久久久久久噜噜老黄 | 国产一区二区三区在线臀色熟女| 夜夜看夜夜爽夜夜摸| 亚洲va日本ⅴa欧美va伊人久久| 大又大粗又爽又黄少妇毛片口| 色尼玛亚洲综合影院| 免费看av在线观看网站| 麻豆精品久久久久久蜜桃| 成人性生交大片免费视频hd| 亚洲欧美清纯卡通| 中亚洲国语对白在线视频| 欧美激情在线99| 日本a在线网址| 99国产精品一区二区蜜桃av| 一个人免费在线观看电影| 国产精品av视频在线免费观看| a级一级毛片免费在线观看| 97超级碰碰碰精品色视频在线观看| 日本一二三区视频观看| 欧美成人一区二区免费高清观看| 69人妻影院| 嫩草影院新地址| 国内毛片毛片毛片毛片毛片| 国产一区二区亚洲精品在线观看| 免费不卡的大黄色大毛片视频在线观看 | 禁无遮挡网站| 国产精品三级大全| 国产精品免费一区二区三区在线| 久久久久久久久久黄片| 日韩欧美三级三区| 国产伦精品一区二区三区视频9| 偷拍熟女少妇极品色| 久久久久久久久久黄片| 亚洲av美国av| 深夜精品福利| 亚洲av美国av| 亚洲国产欧洲综合997久久,| 亚洲av免费在线观看| 五月玫瑰六月丁香| 亚洲自偷自拍三级| 天堂√8在线中文| 国产精品一区二区三区四区久久| 亚洲自偷自拍三级| 亚洲一级一片aⅴ在线观看| 亚洲av美国av| 国产一区二区亚洲精品在线观看| 欧美精品国产亚洲| 国产一区二区在线av高清观看| 免费高清视频大片| 丰满人妻一区二区三区视频av| 日日撸夜夜添| 国产精品,欧美在线| 国产高清视频在线观看网站| 97超视频在线观看视频| 少妇的逼好多水| 国国产精品蜜臀av免费| 黄色视频,在线免费观看| 国内毛片毛片毛片毛片毛片| 国内毛片毛片毛片毛片毛片| 狠狠狠狠99中文字幕| 男女视频在线观看网站免费| 久久久久久久精品吃奶| 91久久精品国产一区二区成人| 日本一二三区视频观看| 午夜老司机福利剧场| 一进一出抽搐动态| 欧美丝袜亚洲另类 | 夜夜看夜夜爽夜夜摸| 动漫黄色视频在线观看| 国产免费一级a男人的天堂| 国内少妇人妻偷人精品xxx网站| 亚洲 国产 在线| 亚洲美女黄片视频| 国产蜜桃级精品一区二区三区| 色综合亚洲欧美另类图片| 日韩欧美一区二区三区在线观看| 亚洲av.av天堂| 亚洲自偷自拍三级| 狂野欧美激情性xxxx在线观看| 中文亚洲av片在线观看爽| 免费观看在线日韩| 婷婷色综合大香蕉| 久久精品影院6| 免费看美女性在线毛片视频| 午夜老司机福利剧场| 久久6这里有精品| 久久国内精品自在自线图片| 最近在线观看免费完整版| 精品人妻熟女av久视频| 久久99热6这里只有精品| 性欧美人与动物交配| 如何舔出高潮| 久久精品国产亚洲av天美| 久久久精品大字幕| 亚洲av中文字字幕乱码综合| av在线亚洲专区| 欧美精品啪啪一区二区三区| 亚洲经典国产精华液单| 夜夜看夜夜爽夜夜摸| 波多野结衣高清无吗| www日本黄色视频网| 国产精品av视频在线免费观看| 亚洲一级一片aⅴ在线观看| 国产男人的电影天堂91| 久久久久久国产a免费观看| 深爱激情五月婷婷| 一区二区三区四区激情视频 | 又爽又黄a免费视频| 免费看av在线观看网站| 亚洲男人的天堂狠狠| 国产精品永久免费网站| 国产综合懂色| 国产精品人妻久久久久久| 亚洲欧美日韩东京热| 女同久久另类99精品国产91| 亚洲中文字幕日韩| 男人舔奶头视频| 美女免费视频网站| 精品乱码久久久久久99久播| 长腿黑丝高跟| 国产精品爽爽va在线观看网站| 国产黄色小视频在线观看| 欧美zozozo另类| 日韩欧美精品免费久久| 能在线免费观看的黄片| 午夜福利成人在线免费观看| 久久久久久久亚洲中文字幕| 免费一级毛片在线播放高清视频| 91狼人影院| 成人特级黄色片久久久久久久| 成人国产一区最新在线观看| 国产精品av视频在线免费观看| 久久久久久久久久黄片| 久久久久久国产a免费观看| 亚洲欧美日韩卡通动漫| 国产淫片久久久久久久久| 亚洲人与动物交配视频| 亚洲av电影不卡..在线观看| 国产在线精品亚洲第一网站| 亚洲精华国产精华液的使用体验 | 丰满的人妻完整版| 97超视频在线观看视频| 欧美性猛交╳xxx乱大交人| 日韩av在线大香蕉| 精品一区二区三区视频在线| 日本撒尿小便嘘嘘汇集6| 国产又黄又爽又无遮挡在线| 91在线观看av| 波野结衣二区三区在线| 此物有八面人人有两片| 草草在线视频免费看| 成人特级av手机在线观看| 国产精品久久视频播放| 狂野欧美白嫩少妇大欣赏| 成人鲁丝片一二三区免费| 级片在线观看| av在线亚洲专区| 国产午夜福利久久久久久| 欧美日本亚洲视频在线播放| 久久久久久久久大av| 综合色av麻豆| 中亚洲国语对白在线视频| 三级毛片av免费| 99久久成人亚洲精品观看| 99久久精品一区二区三区| 国产白丝娇喘喷水9色精品| 久久精品人妻少妇| 男女下面进入的视频免费午夜| 美女高潮喷水抽搐中文字幕| 亚洲av一区综合| 久久午夜亚洲精品久久| 欧美中文日本在线观看视频| 国产在线精品亚洲第一网站| 两性午夜刺激爽爽歪歪视频在线观看| 看免费成人av毛片| av在线蜜桃| 亚洲精品日韩av片在线观看| 亚洲综合色惰| 久久6这里有精品| 国产中年淑女户外野战色| 国产在视频线在精品| 老司机深夜福利视频在线观看| 成人特级黄色片久久久久久久| 日本色播在线视频| 国产一区二区在线av高清观看| 别揉我奶头 嗯啊视频| 亚洲av不卡在线观看| 欧美潮喷喷水| 国产精品国产高清国产av| 亚洲精品一区av在线观看| 亚洲自偷自拍三级| 中文字幕熟女人妻在线| 啪啪无遮挡十八禁网站| 国产精品美女特级片免费视频播放器| 直男gayav资源| 日本一本二区三区精品| 欧美激情国产日韩精品一区| 不卡一级毛片| 成人国产麻豆网| 亚洲av日韩精品久久久久久密| 亚洲va在线va天堂va国产| 两人在一起打扑克的视频| 黄色欧美视频在线观看| 欧美不卡视频在线免费观看| 国产精品福利在线免费观看| 久久人人爽人人爽人人片va| 亚洲最大成人av| 亚洲精品影视一区二区三区av| 精品无人区乱码1区二区| 亚洲国产精品久久男人天堂| 中文资源天堂在线| 三级国产精品欧美在线观看| avwww免费| 欧美日韩乱码在线| 麻豆国产av国片精品| 亚洲avbb在线观看| 国产亚洲精品av在线| 一区二区三区高清视频在线| 99久久精品热视频| 久久久久久久久久黄片| 观看美女的网站| 久久久国产成人免费| 色尼玛亚洲综合影院| 成人高潮视频无遮挡免费网站| 久久久午夜欧美精品| 黄色女人牲交| 偷拍熟女少妇极品色| 国内精品宾馆在线| 日本撒尿小便嘘嘘汇集6| 精品国内亚洲2022精品成人| 国产精品国产高清国产av| 午夜福利在线在线| 国产亚洲91精品色在线| 一个人观看的视频www高清免费观看| 精品人妻熟女av久视频| 天堂影院成人在线观看| 日本熟妇午夜| 色在线成人网| 久久久久久大精品| 亚洲欧美日韩高清专用| 久久精品国产亚洲av天美| 国产伦在线观看视频一区| 国产精品综合久久久久久久免费| 精品乱码久久久久久99久播| 性插视频无遮挡在线免费观看| 久久精品国产99精品国产亚洲性色| 联通29元200g的流量卡| 久久香蕉精品热| 国产伦精品一区二区三区视频9| 免费看美女性在线毛片视频| 国产美女午夜福利| 国产精品一区二区免费欧美| 免费无遮挡裸体视频| 少妇人妻一区二区三区视频| 亚洲精品国产成人久久av| 日本一二三区视频观看| 国产精品永久免费网站| 免费看光身美女| 国产真实伦视频高清在线观看 | 非洲黑人性xxxx精品又粗又长| 九色国产91popny在线| 变态另类成人亚洲欧美熟女| 五月伊人婷婷丁香| 欧美性猛交╳xxx乱大交人| 久99久视频精品免费| 天堂网av新在线| 熟女电影av网| 国产免费av片在线观看野外av| 搡老妇女老女人老熟妇| 国内精品久久久久久久电影| 日日摸夜夜添夜夜添小说| 亚洲一区二区三区色噜噜| 免费电影在线观看免费观看| 最近中文字幕高清免费大全6 | 国产色爽女视频免费观看| 欧美不卡视频在线免费观看| 日本与韩国留学比较| 国产精品爽爽va在线观看网站| 麻豆国产av国片精品| 国产在视频线在精品| 成年版毛片免费区| 国产 一区精品| АⅤ资源中文在线天堂| 午夜福利视频1000在线观看| 欧美极品一区二区三区四区| 免费av观看视频| 成人国产一区最新在线观看| 亚洲av美国av| 桃红色精品国产亚洲av| 校园春色视频在线观看| 在线国产一区二区在线| 国产乱人视频| 女生性感内裤真人,穿戴方法视频| 桃红色精品国产亚洲av| 国产精品永久免费网站| 内射极品少妇av片p| 免费看a级黄色片| 五月伊人婷婷丁香| 女生性感内裤真人,穿戴方法视频| 乱人视频在线观看| 日韩欧美精品免费久久| 天堂影院成人在线观看| 亚洲精品一卡2卡三卡4卡5卡| 精品国产三级普通话版| 岛国在线免费视频观看| 久久久午夜欧美精品| 亚洲成人中文字幕在线播放| 欧美性猛交黑人性爽| 午夜激情欧美在线| 国产在线男女| 欧美+亚洲+日韩+国产| 国产精品三级大全| 日韩 亚洲 欧美在线| 欧美高清成人免费视频www| 亚洲一区二区三区色噜噜| 亚洲av成人av| 高清毛片免费观看视频网站| 精品人妻一区二区三区麻豆 | 九九爱精品视频在线观看| 看免费成人av毛片| 男女下面进入的视频免费午夜| 最近中文字幕高清免费大全6 | 琪琪午夜伦伦电影理论片6080| 99热这里只有是精品50| 久久久成人免费电影| 三级毛片av免费| 久久精品国产清高在天天线| 国产高清视频在线观看网站| 春色校园在线视频观看| 校园春色视频在线观看| 欧美精品啪啪一区二区三区| 美女免费视频网站| 一区二区三区四区激情视频 | 国产高潮美女av| 精品一区二区免费观看| 欧美zozozo另类| 国产精品久久电影中文字幕| 在线观看66精品国产| 免费不卡的大黄色大毛片视频在线观看 | 两个人视频免费观看高清| 欧美丝袜亚洲另类 | 国产男人的电影天堂91| 男人和女人高潮做爰伦理| 国产色婷婷99| 日本黄大片高清| 色哟哟·www| 人妻少妇偷人精品九色| 成人欧美大片| 国产黄片美女视频| 男人舔奶头视频| 亚洲自偷自拍三级| 一级毛片久久久久久久久女| 成年女人毛片免费观看观看9| 特级一级黄色大片| 精品久久久久久久久亚洲 | 在线播放无遮挡| 国产男靠女视频免费网站| av在线天堂中文字幕| 99热这里只有是精品在线观看| 欧美日韩国产亚洲二区| 少妇被粗大猛烈的视频| 性插视频无遮挡在线免费观看| 日本三级黄在线观看| 亚洲专区国产一区二区| 97超视频在线观看视频| 欧美色视频一区免费| 国产成人a区在线观看| 欧美+亚洲+日韩+国产| 男女之事视频高清在线观看| 女的被弄到高潮叫床怎么办 | 国产私拍福利视频在线观看| 久久久久国内视频| 精品久久久久久久久久久久久| 日本成人三级电影网站| 久久99热这里只有精品18| 男插女下体视频免费在线播放| 大又大粗又爽又黄少妇毛片口| 波多野结衣高清无吗| 一级毛片久久久久久久久女| 欧美3d第一页| 熟妇人妻久久中文字幕3abv| 午夜福利在线在线| 日日干狠狠操夜夜爽| 国产色爽女视频免费观看| 亚洲一区二区三区色噜噜| 亚洲无线在线观看| 亚洲精品456在线播放app | 国产一区二区三区在线臀色熟女| 久久精品国产亚洲网站| av在线亚洲专区| 色综合亚洲欧美另类图片| 真人一进一出gif抽搐免费| 欧美日本视频| 成人国产麻豆网| 日韩亚洲欧美综合| 免费av观看视频| 3wmmmm亚洲av在线观看| 男女之事视频高清在线观看| 级片在线观看| 国产一区二区三区视频了| 精品久久久久久久久av| 韩国av在线不卡| 99久久九九国产精品国产免费| 国产精品av视频在线免费观看| 在现免费观看毛片| 中文资源天堂在线| 欧美bdsm另类| 久久午夜亚洲精品久久| 韩国av一区二区三区四区| 精品久久久久久久久亚洲 | 婷婷精品国产亚洲av| 日韩欧美三级三区| 韩国av在线不卡| 少妇高潮的动态图| 97超视频在线观看视频| 久久午夜亚洲精品久久| 少妇的逼水好多| 美女大奶头视频| 久久国内精品自在自线图片| 91在线观看av| 午夜a级毛片| 春色校园在线视频观看| 国产精品不卡视频一区二区| 男女做爰动态图高潮gif福利片| 国产淫片久久久久久久久| 国产精品99久久久久久久久| 久久午夜福利片| 老女人水多毛片| 神马国产精品三级电影在线观看| 亚洲av熟女| 精品一区二区三区av网在线观看| 自拍偷自拍亚洲精品老妇| 亚洲五月天丁香| 日韩大尺度精品在线看网址| av在线观看视频网站免费| 99久久无色码亚洲精品果冻| 99精品在免费线老司机午夜| av在线观看视频网站免费| 中亚洲国语对白在线视频| 国产探花在线观看一区二区| 18禁在线播放成人免费| 一本一本综合久久| 日韩国内少妇激情av| 免费不卡的大黄色大毛片视频在线观看 | 国产成人aa在线观看| 午夜日韩欧美国产| 精品国产三级普通话版| 一本一本综合久久| 欧美人与善性xxx| 国产伦精品一区二区三区视频9| 中文字幕精品亚洲无线码一区| 一级黄色大片毛片| 九色成人免费人妻av| 校园春色视频在线观看| 国产免费一级a男人的天堂| 午夜福利在线观看吧| 成人一区二区视频在线观看| 久久久久久九九精品二区国产| 欧美性猛交╳xxx乱大交人| 久久人妻av系列| 少妇人妻一区二区三区视频| 亚洲av免费高清在线观看| 欧美绝顶高潮抽搐喷水| 国产男人的电影天堂91| 精品日产1卡2卡| ponron亚洲| 欧美xxxx黑人xx丫x性爽| 国产亚洲精品综合一区在线观看| 国产蜜桃级精品一区二区三区| 成人性生交大片免费视频hd| 看十八女毛片水多多多| 综合色av麻豆| 尤物成人国产欧美一区二区三区| 日韩精品青青久久久久久| 亚洲精品影视一区二区三区av| 亚洲av中文字字幕乱码综合| 国产91精品成人一区二区三区| xxxwww97欧美| 日韩中文字幕欧美一区二区| 搡老妇女老女人老熟妇| 国产成人a区在线观看| 少妇丰满av| 欧美丝袜亚洲另类 | 国产高清有码在线观看视频| 窝窝影院91人妻| 九九久久精品国产亚洲av麻豆| 久久久久性生活片| av.在线天堂| 极品教师在线免费播放| 亚洲在线观看片| 久久精品国产自在天天线| 精品无人区乱码1区二区| 一级黄色大片毛片| 国内毛片毛片毛片毛片毛片| av在线老鸭窝| 精品久久久久久久久亚洲 | 1000部很黄的大片| 色精品久久人妻99蜜桃| 欧美绝顶高潮抽搐喷水| 丰满乱子伦码专区| 午夜久久久久精精品| 美女大奶头视频| 在线免费观看不下载黄p国产 | 天美传媒精品一区二区| 国产精品电影一区二区三区| 亚洲综合色惰| 国产v大片淫在线免费观看| 男女啪啪激烈高潮av片| 禁无遮挡网站| 亚洲国产精品合色在线| 亚洲av日韩精品久久久久久密| 真实男女啪啪啪动态图| 欧美日韩精品成人综合77777| 一进一出抽搐动态| 国产精品综合久久久久久久免费| videossex国产| 看免费成人av毛片| 中文字幕熟女人妻在线| 大又大粗又爽又黄少妇毛片口| 久久这里只有精品中国| 夜夜夜夜夜久久久久| 精品人妻偷拍中文字幕| 午夜福利视频1000在线观看| 亚洲经典国产精华液单| 嫩草影院精品99| 午夜精品一区二区三区免费看| 久久久成人免费电影| 人妻久久中文字幕网| 亚洲精品日韩av片在线观看| 少妇人妻精品综合一区二区 | 日本精品一区二区三区蜜桃| 亚洲精品乱码久久久v下载方式| av在线观看视频网站免费| 久9热在线精品视频| 亚洲成人免费电影在线观看| 久久亚洲真实| 亚洲精品一卡2卡三卡4卡5卡| 色播亚洲综合网| 亚洲最大成人中文| 国产精品自产拍在线观看55亚洲| 尤物成人国产欧美一区二区三区| avwww免费| 成熟少妇高潮喷水视频| 亚洲国产精品sss在线观看| 一级黄片播放器| 精品人妻偷拍中文字幕| 熟女人妻精品中文字幕| 亚洲国产高清在线一区二区三| 日本 欧美在线| 婷婷精品国产亚洲av| av视频在线观看入口| 欧美日韩瑟瑟在线播放| 丝袜美腿在线中文| 日韩精品中文字幕看吧| 男女之事视频高清在线观看| 欧美+日韩+精品| 俄罗斯特黄特色一大片| 亚洲精华国产精华液的使用体验 | 成人一区二区视频在线观看| 免费电影在线观看免费观看| 九色国产91popny在线| 我的女老师完整版在线观看| 欧美国产日韩亚洲一区| 丰满的人妻完整版| 在线看三级毛片| 亚洲精品一区av在线观看| 国产视频内射| 婷婷精品国产亚洲av在线| 韩国av在线不卡| 亚洲乱码一区二区免费版| 桃色一区二区三区在线观看| 久久欧美精品欧美久久欧美| 亚洲成a人片在线一区二区| 日日啪夜夜撸| 韩国av在线不卡| 国产成人aa在线观看| 色哟哟哟哟哟哟| 在现免费观看毛片| 国产精品综合久久久久久久免费| 成年人黄色毛片网站| 身体一侧抽搐| 亚洲成人久久爱视频| 亚洲精品亚洲一区二区| 亚洲欧美激情综合另类| 此物有八面人人有两片| 热99re8久久精品国产| 午夜免费成人在线视频| 国产69精品久久久久777片| 亚洲乱码一区二区免费版| 精品午夜福利视频在线观看一区| 午夜福利18| 久久久国产成人精品二区| 成人三级黄色视频| 国产v大片淫在线免费观看| 精品国产三级普通话版| 欧洲精品卡2卡3卡4卡5卡区| 少妇丰满av| 天堂影院成人在线观看| 97超视频在线观看视频| 欧美最黄视频在线播放免费| 哪里可以看免费的av片| 免费一级毛片在线播放高清视频| 欧美性猛交黑人性爽| av在线亚洲专区| 狠狠狠狠99中文字幕| 又紧又爽又黄一区二区| 男人的好看免费观看在线视频| 三级国产精品欧美在线观看| 国产精品一区www在线观看 | 在线免费十八禁| 亚洲 国产 在线| 欧美丝袜亚洲另类 | 美女 人体艺术 gogo| 亚洲美女黄片视频| 亚洲不卡免费看|