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

    荊江沙市段分流比計(jì)算公式的改進(jìn)及應(yīng)用*

    2020-05-08 03:06:00立,陳帆*,張為,袁
    湖泊科學(xué) 2020年3期
    關(guān)鍵詞:枯水糙率水深

    陳 立,陳 帆*,張 為,袁 晶

    (1:武漢大學(xué)水資源與水電工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,武漢 430072) (2:長(zhǎng)江水利委員會(huì)水文局,武漢 430010)

    分汊河段在長(zhǎng)江中下游大量存在,占總河長(zhǎng)的75%以上[1]. 分汊河段具有主流擺動(dòng)頻繁,洲灘消長(zhǎng),主、支汊興衰交替的演變特征[2-3]. 在以懸移質(zhì)運(yùn)動(dòng)為主的沖積河流中,汊道分流很大程度上決定了汊道分沙[4],分流比的變化會(huì)引發(fā)各汊輸沙動(dòng)力強(qiáng)弱的改變,造成汊道的消長(zhǎng). 中枯水分流比反映了基本河槽的分流狀況,因此被作為判斷汊道發(fā)展?fàn)顟B(tài)的重要指標(biāo)[5]. 顯然,研究汊道分流比在不同水沙、地形邊界條件下的變化規(guī)律,對(duì)于研究汊道沖淤演變機(jī)理和航道治理實(shí)踐具有重大意義[6].

    從1940s開始,陸續(xù)有學(xué)者通過水槽試驗(yàn)來研究明渠分汊的水流結(jié)構(gòu). Lakshmana等[7]和Ramamurthy等[8]先后建立了直角明渠分汊中支汊分流比的計(jì)算式,但這些公式僅適用于直角分汊,無法應(yīng)用于邊界條件更為復(fù)雜的分汊河流. 天然汊道中分流比可基于水面線的推求[9]或數(shù)學(xué)模型計(jì)算[10-11]獲得,但對(duì)計(jì)算結(jié)果影響巨大的河道阻力分布大多通過實(shí)測(cè)資料率定反求[12]或直接用經(jīng)驗(yàn)公式[13]計(jì)算,難以反映水沙、地形邊界條件改變后阻力發(fā)生的調(diào)整. 而由河床演變基本理論基于各種假設(shè)推求的分流比計(jì)算式[14-16],在已知單個(gè)汊道斷面內(nèi)各汊的過水面積、水深、糙率、比降等數(shù)值后即可計(jì)算得到分流比,避免了復(fù)雜的流場(chǎng)求解,但對(duì)糙率、比降不同的處理方式造成了各分流比計(jì)算式精度的差異. 可見,無論是通過求解汊道段流場(chǎng)分布還是利用分流比計(jì)算公式來求得分流比,汊道阻力分布的確定都是重難點(diǎn)所在.

    本文研究?jī)H限于雙分汊河段的分流比計(jì)算式.

    1 研究區(qū)域與研究方法

    1.1 研究河段

    沙市段位于上荊江河段,上起陳家灣、下至玉和坪,總長(zhǎng)約18 km(圖1). 河段進(jìn)口左岸有沮漳河入?yún)R,右岸有太平口分流,自上而下根據(jù)平面外形的不同可分為太平口順直分汊和三八灘微彎分汊兩個(gè)河段,其中太平口分汊段由太平口心灘分為北(左)槽和南(右)槽,三八灘分汊段由三八灘分為北(左)汊和南(右)汊. 各分汊段的形態(tài)特征指標(biāo)見表1.

    表1 研究河段汊道形態(tài)特征

    1.2 研究數(shù)據(jù)及方法

    研究數(shù)據(jù)包括2003-2019年沙市段兩分汊段典型斷面1#、2#(斷面位置見圖1)在不同流量下測(cè)得的流速、水位、比降、河寬等,其中,1#斷面有3500~27000 m3/s流量下共23個(gè)測(cè)次的數(shù)據(jù),2#斷面有6000~8000 m3/s流量下共12個(gè)測(cè)次的數(shù)據(jù). 上述35個(gè)測(cè)次數(shù)據(jù)中,24個(gè)包含有比降實(shí)測(cè)值. 研究還收集了2003-2018年沙市段地形圖.

    分析計(jì)算過程中,過水面積由水深沿河寬積分得到,過水面積除以水面寬得到平均水深,汊道流量則由單寬流量(垂線平均流速與水深之積)沿汊寬積分得到,分流比根據(jù)兩汊的流量求出. 沿兩汊深泓線從分流點(diǎn)到匯流點(diǎn)間的距離作為相應(yīng)的汊長(zhǎng),心灘灘頭脊線與分流點(diǎn)處兩汊水動(dòng)力軸線的夾角為相應(yīng)的偏轉(zhuǎn)角(圖1),1-1#~1#、1#~1-2#、1-1#~1-2#之間的比降分別為汊道進(jìn)口段、出口段和整個(gè)汊道段的平均比降.

    本文的研究方法為:采用實(shí)測(cè)數(shù)據(jù)檢驗(yàn)并對(duì)比現(xiàn)有分流比計(jì)算公式的精度,在此基礎(chǔ)上,針對(duì)基于曼寧公式建立的分流比計(jì)算式,分析糙率、比降不同表達(dá)方法對(duì)分流比計(jì)算精度的影響,提出改進(jìn)后的分流比計(jì)算公式型式,利用實(shí)測(cè)數(shù)據(jù)率定公式的參數(shù)后得到新的分流比計(jì)算經(jīng)驗(yàn)公式,并進(jìn)一步分析沙市段汊道的分流比變化特點(diǎn).

    圖1 長(zhǎng)江荊江河段沙市段河勢(shì)圖Fig.1 Sketch of the Jingjiang River and Shashi anabranching reach

    2 分流比計(jì)算公式的比較

    2.1 分流比計(jì)算公式

    2.1.1 基于曼寧公式的分流比計(jì)算式 基于曼寧公式較易得到如下的分流比計(jì)算式[14]:

    (1)

    式中,A、H、n和J分別表示汊道斷面過水面積、平均水深、綜合糙率系數(shù)和縱比降,下標(biāo)1、2分別表示兩汊. 兩汊綜合糙率系數(shù)表征了兩汊的綜合阻力,由于缺乏天然河道糙率實(shí)測(cè)數(shù)據(jù),因此采用以下的近似處理辦法:

    鍵槽凸模的間歇工作通過智能控制裝置來實(shí)現(xiàn)。智能控制裝置由氣動(dòng)執(zhí)行元件(氣缸)、磁性開關(guān)、氣源處理元件、氣動(dòng)控制元件(方向控制閥)、PLC、接近開關(guān)等元件組成, 實(shí)現(xiàn)可換位沖鍵槽定子圓片落料模在沖床沖次在25~30次/min情況下,完成4個(gè)鍵槽凸模循環(huán)連續(xù)沖制沖片,并實(shí)時(shí)監(jiān)測(cè)模具狀態(tài),出現(xiàn)故障及時(shí)停機(jī)。智能控制裝置的控制系統(tǒng)由兩個(gè)部分組成:(1)PLC程序,用來檢測(cè)1號(hào)到4號(hào)氣缸頂出和縮回時(shí)是否到位。(2)指示燈,用來循環(huán)顯示1號(hào)到4號(hào)氣缸的狀態(tài)。

    1)假定兩汊糙率、比降均相差不大,即n1/n2(J2/J1)1/2≈1. 汊道分流比可利用水力幾何變量(A和H)計(jì)算得到[14]:

    (2)

    2)假定從分流點(diǎn)至匯流點(diǎn)兩汊水頭差相等,即Δh1=Δh2,J2/J1=(Δh2/L2)/(Δh1/L1)=L1/L2,其中L表示汊長(zhǎng),h表示水位. 此法較為常用,也稱等水位差法[14],利用相對(duì)較易測(cè)量的汊長(zhǎng)代替比降,即可得到分流比計(jì)算式:

    (3)

    2.1.2 等動(dòng)量法分流比計(jì)算式 等動(dòng)量法也稱動(dòng)量平衡法,由童朝鋒[15]最先提出. 他認(rèn)為分流段兩分流水體在垂直于河道主軸方向上的動(dòng)量相等,即Q1·U1·sinθ1=Q2·U2·sinθ2,其中Q和U分別表示汊道流量和平均流速,θ表示偏轉(zhuǎn)角. 將Q=U·A代入上式,即可得等動(dòng)量法分流比計(jì)算式:

    (4)

    2.1.3 等含沙量法分流比計(jì)算式 文獻(xiàn)[16]指出,對(duì)于河床相對(duì)穩(wěn)定的分汊河段,各汊的含沙量是比較接近的. 而處于平衡狀態(tài)的河流,含沙量近似等于其挾沙力. 根據(jù)張瑞瑾挾沙力計(jì)算公式可得到:

    (5)

    式中,下標(biāo)0表示上游單一段,K和m分別為挾沙力系數(shù)和指數(shù),一般可認(rèn)為變化不大. 進(jìn)一步假設(shè)分汊前后泥沙粒徑不變,則各汊內(nèi)泥沙沉速相等. 將Q=U·A代入上式后,即可得到等含沙量法分流比計(jì)算式[15]:

    (6)

    圖2 分流比計(jì)算結(jié)果對(duì)比Fig.2 Measured versus calculated flow division ratio

    表2 分流比計(jì)算式計(jì)算誤差

    2.2 不同分流比計(jì)算式的精度

    用含有比降的實(shí)測(cè)數(shù)據(jù)檢驗(yàn)各公式的計(jì)算精度,分別計(jì)算北槽和北汊的分流比并與實(shí)測(cè)值作比較,結(jié)果見圖2.

    式(2)不考慮兩汊糙率和比降的差異,最大計(jì)算偏差可達(dá)15%,中、洪水流量下計(jì)算值與實(shí)測(cè)值之間的吻合程度高,計(jì)算偏差不超過3%,說明枯水流量下忽略糙率、比降差異會(huì)造成較大的偏差,中洪水流量下影響不大;式(3)中保留式(2)n1≈n2的假設(shè)用汊長(zhǎng)代替比降后,計(jì)算誤差略有減小,但變化不大;式(4)在4個(gè)公式中計(jì)算誤差最大,說明偏轉(zhuǎn)角θ不能全面反映兩汊糙率、比降對(duì)分流比的影響;式(6)與式(2)形式相同,僅水深之比的指數(shù)不同,其最大計(jì)算偏差約為11%.

    總體上,現(xiàn)有分流比計(jì)算公式考慮兩汊糙率和比降差異的程度不夠,雖然型式簡(jiǎn)單,但精度不高.

    3 分流比計(jì)算式的改進(jìn)

    準(zhǔn)確反映兩汊糙率和比降的差異是提高分流比計(jì)算公式精度的關(guān)鍵. 本節(jié)將主要探討分流比計(jì)算中如何選取代表性比降以及能夠綜合反映糙率、比降差異的參數(shù).

    圖3 考慮不同的比降對(duì)分流比計(jì)算結(jié)果的影響Fig.3 Impacts on flow division ratio calculation by considering different water surface slopes

    表3 分流比計(jì)算式計(jì)算誤差

    3.1 選取不同的比降對(duì)分流比計(jì)算精度的影響

    汊道內(nèi)比降沿程變化[5],進(jìn)口段、出口段和整個(gè)汊道段的平均比降不同,選擇不同的比降計(jì)算分流比結(jié)果不同. 將式(1)中取n1/n2=1,然后分別采用進(jìn)口段、出口段和汊道的平均比降計(jì)算分流比,即在式(2)的基礎(chǔ)上考慮比降差異的影響,太平口心灘分汊段計(jì)算結(jié)果見圖3和表3. 可以看出:和不考慮比降影響的式(2)相比,采用進(jìn)口段比降和汊道平均比降后,分流比的計(jì)算誤差均有所減小,其中采用進(jìn)口段比降計(jì)算,誤差相對(duì)更??;但只考慮比降的差異,計(jì)算結(jié)果仍存在著較大的偏差(最大偏差約為9%).

    3.2 反映兩汊糙率、比降差異的綜合指標(biāo)

    對(duì)兩汊分別運(yùn)用曼寧公式相除后可得到:

    (7)

    定義n/J1/2為糙率比降綜合影響系數(shù)C(后文簡(jiǎn)稱綜合系數(shù)),則上式表達(dá)的是綜合系數(shù)之比C1/C2,與水力半徑之比和流速之比相關(guān),也與兩汊的平均水深之比和Froude數(shù)之比相關(guān),而水深與Froude數(shù)之間存在著內(nèi)在聯(lián)系. 分別點(diǎn)繪綜合系數(shù)比與平均水深比和Froude數(shù)比之間的關(guān)系(圖4).

    可以看出,兩汊綜合系數(shù)比與Froude數(shù)之比高度相關(guān)、與平均水深比強(qiáng)相關(guān),雖然Froude數(shù)之比能更充分反映綜合系數(shù)差異對(duì)分流比的影響,但由于Froude數(shù)中顯含流速U,而分流比的計(jì)算實(shí)際上就是要確定U,因此本文選取平均水深之比來反映綜合系數(shù)比對(duì)分流比的影響. 對(duì)圖4中24組數(shù)據(jù)擬合后可得平均水深之比與綜合系數(shù)之比的經(jīng)驗(yàn)表達(dá)式,將其代入式(1)中可得到:

    (8)

    式(8)即為采用平均水深之比代替綜合系數(shù)之比的分流比計(jì)算式,其中,a、b可通過實(shí)測(cè)資料線性擬合得到. 采用實(shí)測(cè)數(shù)據(jù)對(duì)該公式進(jìn)行驗(yàn)證,結(jié)果見表3和圖5. 由圖表可見,該計(jì)算式的計(jì)算誤差顯著減小,各級(jí)流量下的分流比計(jì)算偏差均小于5%.

    圖4 Froude數(shù)之比、平均水深之比與綜合系數(shù)之比的相關(guān)關(guān)系Fig.4 Relationships between Fr1/Fr2 and C1/C2, H1/H2 and C1/C2

    順直分汊段和微彎分汊段的平面形態(tài)不同,a、b的取值會(huì)有明顯的差異. 以2010年4月-2018年3月沖淤變化引起的水深比、綜合系數(shù)比的改變?yōu)槔?,太平口心灘分汊段北槽沖刷(圖6),北槽與南槽水深之比(H1/H2)增加約15%,兩汊綜合系數(shù)比(C1/C2)增長(zhǎng)約28%,北槽枯水分流比增加約17%;與之相反,三八灘北汊大幅萎縮(圖6),北汊與南汊水深之比H1/H2減小約36%,C1/C2僅減小約21%,北汊枯水分流比減小達(dá)64%. 顯然,不同類型的汊道內(nèi),兩汊水深之比的變化對(duì)綜合系數(shù)之比和分流比的改變程度不同,主、支均勢(shì)的順直分汊段分流比對(duì)于汊道沖淤的敏感度較低[19],A1/A2、H1/H2小幅度的增加往往引起C1/C2顯著的增加,因而η1增幅較??;而主、支差異明顯的微彎分汊段,分流比受斷面沖淤變化的調(diào)節(jié)作用更為明顯,當(dāng)北汊沖刷發(fā)展(或淤積萎縮)時(shí),隨著A1/A2、H1/H2的增加(或減小),C1/C2的增幅(或減幅)相對(duì)較小,η1隨之發(fā)生顯著的增加(或減小).

    綜上所述,本文提出的分流比計(jì)算式可準(zhǔn)確地反映分汊河段汊道沖淤變化條件下的分流比變化,由于不同的分汊河段內(nèi),兩汊水深之比的變化對(duì)綜合系數(shù)之比和分流比的改變程度不同,因此應(yīng)用在具體的分汊河段時(shí),須采用實(shí)測(cè)資料率定公式中的系數(shù)a、b.

    圖5 改進(jìn)后分流比計(jì)算式計(jì)算結(jié)果對(duì)比Fig.5 Measured versus calculated flow division ratio by the improved formula and other ones

    圖6 分汊斷面2010年4月-2018年3月沖淤變化Fig.6 Variations in the branching cross-section profiles from April of 2010 to March of 2018

    4 三峽水庫(kù)蓄水后沙市段分流比變化特點(diǎn)分析

    4.1 汊道分流比隨水位的變化

    利用式(8),計(jì)算1#斷面處A1/A2、H1/H2和η1隨水位的變化過程. 結(jié)果顯示,當(dāng)計(jì)算所選取的地形邊界不同時(shí),η1隨水位的變化情況是不同的,有先增后減、先減后增、逐漸增加和逐漸減小等多種變化趨勢(shì)(圖7),與以往研究認(rèn)為支汊分流比隨流量的增加而增加[20]并不相同. 這是因?yàn)樵谧匀粭l件下枯水期主流流經(jīng)主汊,隨著流量的增加,主流所具有的“小水坐彎,大水趨直”的特性使得其在一定程度上偏離主汊、偏向支汊,因而支汊的分流比隨流量的增加而增加. 而在兩汊差別不大,主流擺動(dòng)頻繁的順直分汊段,上游來沙減少引起汊道劇烈的沖淤調(diào)整,汊道沖淤交替導(dǎo)致主汊位置的不固定,使得順直分汊段分流比隨水位的變化規(guī)律呈多樣性. 盡管不同地形邊界下,分流比隨水位的變化趨勢(shì)有所不同,但H1/H2與η1的變化趨勢(shì)總體較為一致.

    圖7 不同地形邊界下A1/A2、H1/H2和η1隨水位的變化Fig.7 Variation curves of A1/A2, H1/H2 and η1 versus water level under different channel boundaries

    4.2 蓄水后沙市段汊道分流比年際變化特征

    利用式(8)計(jì)算不同地形邊界下枯、中、洪水流量分流比的變化,分析蓄水后不同流量級(jí)分流比的年際變化特征. 方法為:針對(duì)不同的實(shí)測(cè)1#、2#斷面地形,利用相應(yīng)的沙市水文站水位流量關(guān)系確定枯、中、洪水流量(8000、12000和25000 m3/s)所對(duì)應(yīng)的水位,由此得到各地形邊界下的A1/A2和H1/H2. 代入式(8)計(jì)算η1,結(jié)果見圖8.

    圖8 沙市分汊段不同流量分流比的年際變化Fig.8 Temporal variations in flow division ratio of different flow discharges at the Shashi anabranching reach

    可以看出,三峽水庫(kù)蓄水后,太平口北槽洪、中、枯水分流比均呈總體減小的趨勢(shì). 蓄水初期,北槽枯水分流比經(jīng)歷了從低于50%增加到超過50%,而后又減小至低于50%的轉(zhuǎn)變,即主支汊發(fā)生了較頻繁的易位. 2008年三峽水庫(kù)試驗(yàn)性蓄水后,各流量下北槽分流比有所減小,2010年后有增有減,整體變化不大,現(xiàn)枯水分流比穩(wěn)定在40%~45%;蓄水初期,三八灘分汊段各流量下分流比變化較小,始終保持南汊占優(yōu)的局面. 2005-2008年和2010-2012年,北汊分流比經(jīng)歷了兩次明顯的增加,中枯水分流比超過50%. 2013年后,北汊分流比迅速減小,至2015年小于50%,也發(fā)生了主支汊易位,其后減小趨勢(shì)仍在繼續(xù),但速度放緩,現(xiàn)枯水分流比低于20%.

    可見,太平口北槽、三八灘北汊分流比在三峽水庫(kù)蓄水初期變化劇烈,后期由于連續(xù)實(shí)施的航道整治工程對(duì)河勢(shì)的控制作用逐漸顯現(xiàn),加上上游來水持續(xù)偏枯,因而分流比調(diào)整的劇烈程度降低. 隨著上游梯級(jí)水庫(kù)群的建成運(yùn)用,下游分汊河段趨勢(shì)性調(diào)整仍將延續(xù)這一特點(diǎn),特別是三八灘北汊繼續(xù)趨于萎縮. 但在出現(xiàn)大洪水時(shí),河段內(nèi)分流比的變幅會(huì)增大.

    比較不同流量下分流比變化過程可發(fā)現(xiàn),不管是順直還是微彎分汊段,流量越大,分流比變化越小,流量越小,分流比的調(diào)整則越劇烈. 這主要是因?yàn)樾钏筱獾罌_淤主要集中于枯水河槽[18],同時(shí)枯水位明顯降低[21],使得兩汊枯水流量下水深變幅較大,對(duì)分流的影響也更大.

    4.3 同一流量下汊道分流比隨汊道沖淤的變化

    對(duì)同流量下某一固定的汊道斷面,設(shè)t1時(shí)刻左、右汊過水面積為A1、A2,左汊分流比為η1,t2時(shí)刻分別為A1′、A2′和η1′,則該時(shí)段內(nèi)兩汊單位河長(zhǎng)發(fā)生的沖淤體積可分別表示為 (A1′-A1)×1和 (A2′-A2)×1. 令ER=A1′-A1-(A2′-A2),Δη1=η1′-η1,不同流量下Δη1與ER的關(guān)系如圖9所示.

    由圖9可見,ER與左汊分流比變化之間呈較好的線性關(guān)系,左汊相對(duì)于右汊的沖刷發(fā)展(或淤積萎縮)會(huì)改善(或惡化)其入流條件,分流比隨之增加(或減小). 擬合出Δη1與ER之間的線性關(guān)系表達(dá)為:

    Δη1′=k1·ER+k2

    (9)

    式中,k1和k2分別表示直線斜率和截距,線性擬合的相關(guān)系數(shù)R及k1的取值見表4.

    從表4可以看出,隨著流量的增加,k1逐漸減小,說明沖淤變化對(duì)分流比的影響趨于減小,枯水分流比對(duì)于河道沖淤變化的敏感程度要大于洪水分流比. 各級(jí)流量下,微彎分汊段k1取值均大于順直分汊段,進(jìn)一步反映了兩汊分流差異較大的微彎分汊段內(nèi),沖淤變化對(duì)于分流比的影響要大于順直分汊段[19].

    圖9 不同流量下Δη1隨ER的變化Fig.9 Variations of Δη1 with ER at different flow discharges

    表4 線性擬合參數(shù)

    5 結(jié)論

    1)汊道分流比計(jì)算公式中,兩汊糙率、比降差異的準(zhǔn)確反映與否是影響其計(jì)算精度的關(guān)鍵. 完全忽略糙率、比降差異的影響會(huì)給枯水分流比的計(jì)算結(jié)果帶來較大的偏差;只考慮兩汊比降差異而忽略糙率差異的影響,分流比計(jì)算精度仍然較低.

    2)兩汊平均水深之比可較全面反映河道阻力、比降的差異,與糙率、比降綜合影響系數(shù)之比高度相關(guān),將平均水深之比與糙率、比降綜合影響系數(shù)之比引入到分流比計(jì)算公式,得到了較高精度的分流比計(jì)算式,各流量下計(jì)算偏差不大于5%.

    3)三峽水庫(kù)蓄水后,沙市段汊道枯水分流比的變化較為劇烈,洪水分流比變化較小. 汊道沖淤變化對(duì)枯水分流比的影響強(qiáng)于洪水分流比,對(duì)主、支差異較大的微彎分汊段分流比的影響強(qiáng)于主、支均勢(shì)的順直分汊段.

    猜你喜歡
    枯水糙率水深
    梯形河道行洪能力與護(hù)岸糙率的關(guān)系研究
    陜西水利(2023年7期)2023-07-28 09:30:26
    書法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    哲理漫畫
    丹棱縣晚熟桔橙枯水現(xiàn)狀及防控技術(shù)分析
    耕作與栽培(2021年3期)2021-12-01 02:38:13
    柑橘枯水研究進(jìn)展
    新疆阿勒泰哈巴河縣養(yǎng)殖渠人工渠道糙率的試驗(yàn)分析
    陜西水利(2019年5期)2019-06-26 06:27:34
    復(fù)式河道整治設(shè)計(jì)中綜合糙率研究
    大口徑玻璃鋼管道糙率及過流能力分析
    GPS RTK技術(shù)在水深測(cè)量中的應(yīng)用
    欧美性长视频在线观看| 亚洲成人免费av在线播放| 18禁观看日本| 成人av一区二区三区在线看 | 波多野结衣一区麻豆| 免费久久久久久久精品成人欧美视频| 免费日韩欧美在线观看| 女人高潮潮喷娇喘18禁视频| 欧美日韩中文字幕国产精品一区二区三区 | 午夜福利,免费看| 天天操日日干夜夜撸| 18禁裸乳无遮挡动漫免费视频| xxxhd国产人妻xxx| 国产无遮挡羞羞视频在线观看| 国产精品国产av在线观看| 国产成人精品无人区| 亚洲av男天堂| 99国产精品99久久久久| 亚洲伊人色综图| 美女高潮到喷水免费观看| 国产成人影院久久av| 国产精品 欧美亚洲| 欧美人与性动交α欧美软件| 久久久国产欧美日韩av| 日本a在线网址| av片东京热男人的天堂| 欧美日韩国产mv在线观看视频| 91字幕亚洲| 捣出白浆h1v1| 最黄视频免费看| xxxhd国产人妻xxx| 亚洲精品一区蜜桃| 狠狠婷婷综合久久久久久88av| 国产麻豆69| 丝袜美足系列| 精品亚洲乱码少妇综合久久| 狠狠狠狠99中文字幕| 亚洲精品乱久久久久久| 色婷婷久久久亚洲欧美| 男人添女人高潮全过程视频| 高清黄色对白视频在线免费看| 18禁裸乳无遮挡动漫免费视频| 人妻人人澡人人爽人人| 99国产极品粉嫩在线观看| 视频区欧美日本亚洲| 日本一区二区免费在线视频| 91九色精品人成在线观看| 99久久人妻综合| 女性生殖器流出的白浆| 欧美日韩福利视频一区二区| 久久精品国产综合久久久| 老熟女久久久| 丝袜美腿诱惑在线| 高清av免费在线| 久久久久网色| 久久精品国产亚洲av香蕉五月 | 欧美激情高清一区二区三区| 精品欧美一区二区三区在线| 国产精品1区2区在线观看. | 悠悠久久av| www.999成人在线观看| 在线观看免费日韩欧美大片| 亚洲 欧美一区二区三区| 欧美 日韩 精品 国产| 一二三四社区在线视频社区8| 国产精品成人在线| 日日爽夜夜爽网站| 男人舔女人的私密视频| 91av网站免费观看| 大码成人一级视频| 久久国产亚洲av麻豆专区| 亚洲avbb在线观看| 精品欧美一区二区三区在线| 国产人伦9x9x在线观看| 在线十欧美十亚洲十日本专区| 伊人久久大香线蕉亚洲五| 深夜精品福利| 成年美女黄网站色视频大全免费| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽av亚洲精品天堂| 亚洲av成人一区二区三| 热re99久久精品国产66热6| 国产淫语在线视频| 久久久久精品人妻al黑| 最近中文字幕2019免费版| 中国美女看黄片| 国产成人一区二区三区免费视频网站| 中文精品一卡2卡3卡4更新| 久久综合国产亚洲精品| 9191精品国产免费久久| 老司机在亚洲福利影院| 老汉色∧v一级毛片| 久久国产精品男人的天堂亚洲| 一二三四社区在线视频社区8| 成人亚洲精品一区在线观看| av超薄肉色丝袜交足视频| 亚洲国产av新网站| 99国产精品一区二区三区| 巨乳人妻的诱惑在线观看| 亚洲国产精品一区二区三区在线| 天堂8中文在线网| 悠悠久久av| 久久女婷五月综合色啪小说| 多毛熟女@视频| 一边摸一边做爽爽视频免费| 日韩熟女老妇一区二区性免费视频| 久久久精品国产亚洲av高清涩受| 女人爽到高潮嗷嗷叫在线视频| 少妇精品久久久久久久| 中国美女看黄片| 人人妻人人添人人爽欧美一区卜| 国产野战对白在线观看| 久久精品成人免费网站| xxxhd国产人妻xxx| 亚洲av成人不卡在线观看播放网 | 淫妇啪啪啪对白视频 | 1024视频免费在线观看| 国产成人免费无遮挡视频| 嫩草影视91久久| 国产伦人伦偷精品视频| 午夜成年电影在线免费观看| 国产深夜福利视频在线观看| 韩国精品一区二区三区| 国产成人a∨麻豆精品| 美女国产高潮福利片在线看| 欧美 日韩 精品 国产| 精品亚洲成国产av| 亚洲av美国av| 最近最新中文字幕大全免费视频| 午夜精品久久久久久毛片777| 亚洲精品国产一区二区精华液| 亚洲色图综合在线观看| 亚洲精品第二区| 伊人久久大香线蕉亚洲五| 黄片小视频在线播放| 满18在线观看网站| 一区二区av电影网| 久久综合国产亚洲精品| 日韩制服丝袜自拍偷拍| 亚洲精华国产精华精| 人人澡人人妻人| 亚洲av成人不卡在线观看播放网 | 午夜福利影视在线免费观看| 宅男免费午夜| 黄色片一级片一级黄色片| 美国免费a级毛片| avwww免费| 久久精品亚洲av国产电影网| 无限看片的www在线观看| 欧美黑人精品巨大| 国产成人av教育| 久久久久网色| 成年人免费黄色播放视频| 日韩大片免费观看网站| 动漫黄色视频在线观看| 国产精品一区二区在线观看99| 亚洲国产欧美在线一区| 悠悠久久av| 中文字幕人妻丝袜制服| 免费在线观看日本一区| 亚洲国产毛片av蜜桃av| 日本av免费视频播放| 色播在线永久视频| 亚洲精品久久成人aⅴ小说| 91国产中文字幕| 精品卡一卡二卡四卡免费| 日韩电影二区| 久久精品久久久久久噜噜老黄| 超碰成人久久| 日韩欧美一区二区三区在线观看 | 啦啦啦免费观看视频1| 汤姆久久久久久久影院中文字幕| 美国免费a级毛片| av在线老鸭窝| 欧美黄色片欧美黄色片| 999精品在线视频| 日本91视频免费播放| 国产真人三级小视频在线观看| 国产三级黄色录像| 制服诱惑二区| 久久久久精品国产欧美久久久 | 18禁观看日本| 十八禁网站网址无遮挡| 狂野欧美激情性xxxx| 99久久99久久久精品蜜桃| 欧美 日韩 精品 国产| 99国产综合亚洲精品| 免费av中文字幕在线| 三上悠亚av全集在线观看| 久久人妻熟女aⅴ| 中国国产av一级| 日韩欧美国产一区二区入口| 黄片大片在线免费观看| 精品卡一卡二卡四卡免费| 久久精品aⅴ一区二区三区四区| 亚洲熟女精品中文字幕| 欧美日韩av久久| 国产一区二区三区综合在线观看| 国产精品 欧美亚洲| 久久人人97超碰香蕉20202| 纵有疾风起免费观看全集完整版| 亚洲精品久久午夜乱码| 日韩欧美一区二区三区在线观看 | 深夜精品福利| 一级片'在线观看视频| 老汉色∧v一级毛片| 丝袜美腿诱惑在线| 窝窝影院91人妻| tube8黄色片| 久久精品亚洲av国产电影网| 亚洲专区字幕在线| 19禁男女啪啪无遮挡网站| 91麻豆av在线| 纯流量卡能插随身wifi吗| 国产福利在线免费观看视频| 99久久综合免费| 考比视频在线观看| 久久久久国产一级毛片高清牌| 在线观看www视频免费| 久久人妻福利社区极品人妻图片| 国产精品一区二区免费欧美 | 亚洲精品国产色婷婷电影| 黄色 视频免费看| 亚洲七黄色美女视频| 看免费av毛片| 另类亚洲欧美激情| 亚洲第一欧美日韩一区二区三区 | 国产真人三级小视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 久9热在线精品视频| 女性被躁到高潮视频| 啦啦啦 在线观看视频| 国产亚洲欧美精品永久| 午夜福利视频在线观看免费| 成年人黄色毛片网站| 老司机影院成人| 日韩,欧美,国产一区二区三区| 亚洲国产欧美在线一区| 极品少妇高潮喷水抽搐| 99精品久久久久人妻精品| 下体分泌物呈黄色| 久久人妻熟女aⅴ| 啦啦啦在线免费观看视频4| 久久久国产欧美日韩av| 日韩制服骚丝袜av| 12—13女人毛片做爰片一| 电影成人av| 黄色毛片三级朝国网站| 国产精品国产av在线观看| 亚洲国产av影院在线观看| 丰满饥渴人妻一区二区三| 亚洲av日韩在线播放| 国产1区2区3区精品| 蜜桃在线观看..| 精品久久蜜臀av无| 久久久国产欧美日韩av| 精品第一国产精品| 国产精品久久久久久精品电影小说| 国产一区二区 视频在线| 久久综合国产亚洲精品| 午夜激情av网站| 80岁老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 91大片在线观看| 动漫黄色视频在线观看| 国产在线视频一区二区| 麻豆国产av国片精品| 丁香六月欧美| 国产在线一区二区三区精| 女性生殖器流出的白浆| 日韩一卡2卡3卡4卡2021年| 在线观看舔阴道视频| 国产一区有黄有色的免费视频| 男女下面插进去视频免费观看| 免费久久久久久久精品成人欧美视频| 久久久久国产精品人妻一区二区| 中文字幕人妻丝袜制服| 人妻人人澡人人爽人人| 性高湖久久久久久久久免费观看| 一区二区三区四区激情视频| 久久天堂一区二区三区四区| 国产激情久久老熟女| 国产亚洲精品第一综合不卡| 日本vs欧美在线观看视频| 欧美亚洲日本最大视频资源| 久久天躁狠狠躁夜夜2o2o| 高清欧美精品videossex| 性高湖久久久久久久久免费观看| 啦啦啦在线免费观看视频4| 亚洲精品一二三| 成人三级做爰电影| 国产精品久久久av美女十八| 爱豆传媒免费全集在线观看| 精品人妻1区二区| 亚洲色图 男人天堂 中文字幕| 一级a爱视频在线免费观看| 99国产综合亚洲精品| 免费一级毛片在线播放高清视频 | 中文字幕另类日韩欧美亚洲嫩草| 天天躁日日躁夜夜躁夜夜| 日韩三级视频一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品999| 美女扒开内裤让男人捅视频| 精品国产一区二区三区四区第35| av又黄又爽大尺度在线免费看| 成年av动漫网址| 一区在线观看完整版| 建设人人有责人人尽责人人享有的| 欧美精品人与动牲交sv欧美| 岛国在线观看网站| 叶爱在线成人免费视频播放| 久久国产亚洲av麻豆专区| 91精品三级在线观看| 伦理电影免费视频| 他把我摸到了高潮在线观看 | 9热在线视频观看99| 久久人妻熟女aⅴ| 咕卡用的链子| 美女脱内裤让男人舔精品视频| 美女福利国产在线| 老汉色av国产亚洲站长工具| 精品第一国产精品| 日本精品一区二区三区蜜桃| 久久综合国产亚洲精品| 中文字幕人妻丝袜一区二区| 看免费av毛片| 国产在线免费精品| 少妇猛男粗大的猛烈进出视频| 一二三四社区在线视频社区8| 国产一区二区 视频在线| 久久久欧美国产精品| 久久久国产一区二区| 午夜影院在线不卡| 欧美一级毛片孕妇| 午夜福利视频在线观看免费| 日韩中文字幕视频在线看片| 美女午夜性视频免费| 首页视频小说图片口味搜索| 国产精品久久久av美女十八| av网站在线播放免费| tube8黄色片| 最新在线观看一区二区三区| 极品少妇高潮喷水抽搐| 爱豆传媒免费全集在线观看| 美女午夜性视频免费| av又黄又爽大尺度在线免费看| 丝袜喷水一区| 18禁裸乳无遮挡动漫免费视频| 涩涩av久久男人的天堂| 爱豆传媒免费全集在线观看| 欧美97在线视频| 国产男人的电影天堂91| 午夜免费观看性视频| av网站在线播放免费| 日韩视频在线欧美| 亚洲欧洲日产国产| 亚洲精品一二三| 久久这里只有精品19| 1024视频免费在线观看| 精品国产乱子伦一区二区三区 | 伊人亚洲综合成人网| 50天的宝宝边吃奶边哭怎么回事| 日韩欧美国产一区二区入口| 久久精品国产亚洲av高清一级| 色视频在线一区二区三区| tube8黄色片| 国产一区二区三区av在线| 日韩免费高清中文字幕av| av视频免费观看在线观看| 精品国产一区二区久久| 国产精品一区二区精品视频观看| 国产一区二区激情短视频 | 精品一品国产午夜福利视频| 亚洲一卡2卡3卡4卡5卡精品中文| 美女国产高潮福利片在线看| 久久精品人人爽人人爽视色| 9191精品国产免费久久| 国产伦人伦偷精品视频| 手机成人av网站| av超薄肉色丝袜交足视频| 一级a爱视频在线免费观看| 1024视频免费在线观看| 亚洲精品在线美女| 欧美日韩一级在线毛片| 亚洲少妇的诱惑av| 亚洲综合色网址| 欧美97在线视频| 欧美中文综合在线视频| 亚洲精品一二三| 精品第一国产精品| 国产三级黄色录像| 日韩免费高清中文字幕av| 欧美xxⅹ黑人| 少妇 在线观看| 久久精品成人免费网站| 美女视频免费永久观看网站| 国产麻豆69| 国产野战对白在线观看| 2018国产大陆天天弄谢| 久久毛片免费看一区二区三区| av一本久久久久| 他把我摸到了高潮在线观看 | 欧美日韩黄片免| 欧美人与性动交α欧美精品济南到| 国产高清视频在线播放一区 | 亚洲国产看品久久| 国产精品一区二区在线不卡| www.999成人在线观看| 中文字幕av电影在线播放| 成人免费观看视频高清| 亚洲专区字幕在线| 欧美精品人与动牲交sv欧美| 激情视频va一区二区三区| 亚洲天堂av无毛| 老熟妇仑乱视频hdxx| 国产麻豆69| 少妇精品久久久久久久| www.熟女人妻精品国产| 91精品三级在线观看| 91精品国产国语对白视频| 国产成人精品在线电影| 老司机影院成人| 男女边摸边吃奶| 手机成人av网站| 叶爱在线成人免费视频播放| 黄色a级毛片大全视频| 亚洲伊人色综图| 国产野战对白在线观看| 日韩大码丰满熟妇| 两性午夜刺激爽爽歪歪视频在线观看 | 国产在线免费精品| 人妻人人澡人人爽人人| av网站免费在线观看视频| 黑人巨大精品欧美一区二区mp4| 国产精品成人在线| 老熟女久久久| 2018国产大陆天天弄谢| 欧美黑人欧美精品刺激| 午夜福利视频精品| 最近最新中文字幕大全免费视频| 午夜两性在线视频| 国产成人一区二区三区免费视频网站| 久久国产精品男人的天堂亚洲| 亚洲精品国产精品久久久不卡| 91精品伊人久久大香线蕉| 少妇裸体淫交视频免费看高清 | 黄频高清免费视频| 亚洲专区国产一区二区| 99精品久久久久人妻精品| 女人高潮潮喷娇喘18禁视频| √禁漫天堂资源中文www| 深夜精品福利| 丰满少妇做爰视频| av超薄肉色丝袜交足视频| 超碰97精品在线观看| 亚洲欧美成人综合另类久久久| 国产亚洲精品第一综合不卡| 久久久精品94久久精品| 日本黄色日本黄色录像| 又黄又粗又硬又大视频| a在线观看视频网站| 国产欧美日韩一区二区三区在线| 成年女人毛片免费观看观看9 | 亚洲成人免费电影在线观看| 欧美日韩一级在线毛片| 亚洲第一av免费看| 两个人看的免费小视频| 国产精品自产拍在线观看55亚洲 | 亚洲专区字幕在线| a在线观看视频网站| 三上悠亚av全集在线观看| 亚洲精品乱久久久久久| 久久天堂一区二区三区四区| 午夜福利视频精品| 交换朋友夫妻互换小说| 999久久久精品免费观看国产| 一区二区av电影网| 在线观看免费视频网站a站| 久久国产精品人妻蜜桃| 久久女婷五月综合色啪小说| 水蜜桃什么品种好| 99热国产这里只有精品6| 国产在线观看jvid| 黄片播放在线免费| 黄色视频,在线免费观看| 777久久人妻少妇嫩草av网站| 亚洲熟女精品中文字幕| 国产精品久久久av美女十八| 黄片播放在线免费| 欧美激情高清一区二区三区| 亚洲熟女毛片儿| 日韩制服骚丝袜av| 80岁老熟妇乱子伦牲交| 真人做人爱边吃奶动态| 欧美精品人与动牲交sv欧美| 老熟妇仑乱视频hdxx| 亚洲国产欧美在线一区| 99re6热这里在线精品视频| 一本久久精品| 久久精品亚洲av国产电影网| videos熟女内射| 黑人猛操日本美女一级片| 国产精品九九99| 欧美精品啪啪一区二区三区 | a 毛片基地| 99香蕉大伊视频| 国产xxxxx性猛交| 91老司机精品| 男女之事视频高清在线观看| 精品国产超薄肉色丝袜足j| 少妇人妻久久综合中文| 午夜福利一区二区在线看| 国产免费视频播放在线视频| 美女大奶头黄色视频| 中文字幕制服av| 久久久国产一区二区| 久久久水蜜桃国产精品网| 超碰成人久久| 亚洲天堂av无毛| 巨乳人妻的诱惑在线观看| 亚洲欧美日韩高清在线视频 | 2018国产大陆天天弄谢| 岛国毛片在线播放| 下体分泌物呈黄色| 国产成人av教育| 人妻 亚洲 视频| 亚洲精品av麻豆狂野| 久久精品国产亚洲av高清一级| 不卡一级毛片| 亚洲第一欧美日韩一区二区三区 | 久久人妻熟女aⅴ| 999久久久精品免费观看国产| 18禁裸乳无遮挡动漫免费视频| 久久av网站| 亚洲精品国产av成人精品| 最新的欧美精品一区二区| av在线app专区| 成人18禁高潮啪啪吃奶动态图| 99国产综合亚洲精品| 亚洲成人免费电影在线观看| 首页视频小说图片口味搜索| 人人妻人人爽人人添夜夜欢视频| 欧美性长视频在线观看| 久9热在线精品视频| 精品一区二区三卡| 国产1区2区3区精品| a级毛片在线看网站| 乱人伦中国视频| 欧美精品av麻豆av| 大陆偷拍与自拍| 丰满饥渴人妻一区二区三| 久久久久视频综合| 亚洲精品美女久久av网站| 在线十欧美十亚洲十日本专区| 亚洲一卡2卡3卡4卡5卡精品中文| 一区福利在线观看| 免费一级毛片在线播放高清视频 | 成人18禁高潮啪啪吃奶动态图| 欧美日韩亚洲综合一区二区三区_| 99热全是精品| 日韩大片免费观看网站| 一二三四社区在线视频社区8| 色老头精品视频在线观看| 18禁黄网站禁片午夜丰满| 一区福利在线观看| 亚洲av成人一区二区三| 国产黄频视频在线观看| 欧美在线黄色| 亚洲中文日韩欧美视频| 免费看十八禁软件| 一个人免费看片子| 窝窝影院91人妻| 色精品久久人妻99蜜桃| 女人高潮潮喷娇喘18禁视频| 淫妇啪啪啪对白视频 | av电影中文网址| 如日韩欧美国产精品一区二区三区| 亚洲欧美色中文字幕在线| 亚洲精品av麻豆狂野| 亚洲一区二区三区欧美精品| www.自偷自拍.com| 国产日韩一区二区三区精品不卡| 国产伦人伦偷精品视频| 久久久久久久国产电影| 国产精品一区二区在线不卡| 日本a在线网址| 男人舔女人的私密视频| 中国国产av一级| 最近最新中文字幕大全免费视频| 亚洲全国av大片| 成人18禁高潮啪啪吃奶动态图| 国产免费福利视频在线观看| 亚洲专区中文字幕在线| 中国国产av一级| 国产成人精品久久二区二区91| 丝瓜视频免费看黄片| 亚洲国产精品一区三区| 亚洲欧美精品自产自拍| 欧美+亚洲+日韩+国产| 国产日韩一区二区三区精品不卡| 亚洲av日韩在线播放| 国产精品影院久久| 91老司机精品| 亚洲国产欧美一区二区综合| 亚洲少妇的诱惑av| 脱女人内裤的视频| 国产黄色免费在线视频| 国产在线观看jvid| 精品高清国产在线一区| a在线观看视频网站| 美女中出高潮动态图| 丰满少妇做爰视频| 午夜精品久久久久久毛片777| 欧美午夜高清在线|