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

    崇明東灘鳥類棲息地優(yōu)化工程區(qū)海堤水閘外側(cè)引水渠沖淤數(shù)值模擬

    2020-07-04 02:08:14魯佩儀朱建榮錢偉偉袁琳
    關(guān)鍵詞:數(shù)值模擬

    魯佩儀 朱建榮 錢偉偉 袁琳

    摘要:在河口水動(dòng)力模式的基礎(chǔ)上,耦合泥沙模塊和底部沖淤方程,建立了崇明東灘鳥類保護(hù)區(qū)海堤外側(cè)引水渠沖淤三維數(shù)值模式,模式上游閘門水通量由堰流公式計(jì)算,內(nèi)側(cè)水位由隨塘河水動(dòng)力模式計(jì)算,模式下游水位、鹽度、泥沙濃度由長江河口大區(qū)域數(shù)值模式計(jì)算,利用2018年7月28-30日優(yōu)化工程區(qū)水閘開閘沖淤前后實(shí)測引水渠底部高程進(jìn)行驗(yàn)證,模式計(jì)算的高程與實(shí)測值吻合良好,表明模式能較好地模擬引水渠沖淤變化,沖淤計(jì)算結(jié)果表明,閘下引水渠在自然狀態(tài)下,涵閘關(guān)閉60d后基本淤平,閘口淤積量最大,淤積程度自閘口向外海逐漸減小,引水渠自然淤積2個(gè)月后,隨塘河蓄水至3.0m開閘放水,引水渠平均沖淤量為96mm:隨塘河蓄水至3.8m開閘放水,引水渠平均沖淤量為133mm;納潮沖淤一天,引水渠平均沖淤量為625mm,納潮沖淤是最為有效的沖淤方式,研究成果可為崇明東灘國家級(jí)自然保護(hù)區(qū)生態(tài)修復(fù)工程區(qū)的保護(hù)與管理、引水渠沖淤方案提供科技指導(dǎo)。

    關(guān)鍵詞:崇明東灘;候鳥保護(hù)區(qū);閘下淤積:引水渠沖淤:數(shù)值模擬

    中圖分類號(hào):P731.23文獻(xiàn)標(biāo)志碼:A DOI:10.3969/j,issn,1000-5641.201941013

    0引言

    在我國沿海大部分中小人海河口建閘以擋潮御鹵、排洪蓄淡,以及防臺(tái)抗暴,但河口建閘面臨的一個(gè)重要問題是閘下淤積,我國建閘河口閘下淤積非常嚴(yán)重,顯著影響上游河道的防洪排澇,閘下淤積預(yù)報(bào)和減淤措施研究一直是水利工程中急需解決的問題之一,在以海相來沙為主的潮汐河口建閘,閘下淤積問題相當(dāng)普遍,海河流域35個(gè)河口建閘后,22座發(fā)生嚴(yán)重淤積,江蘇沿海修建的58座大中型擋潮閘中,5座基本淤死,浙江省甬江支流姚江建閘后閘下河段及甬江河床發(fā)生嚴(yán)重淤積,金元?dú)g和沈煥庭指出,建閘河口淤積的主要原因是水動(dòng)力條件發(fā)生了改變,淤泥質(zhì)河口與砂質(zhì)河口沖淤特性差異很大,大部分閘門第一年淤積速度最大,在泥沙含量較高的水域,由于密度差的作用,渾水會(huì)流向清水形成異重流,李大山、任汝述基于對(duì)射陽河建閘前后徑流量、潮動(dòng)力條件等水力要素的比對(duì),對(duì)該地區(qū)泥沙來源和泥沙輸運(yùn)量進(jìn)行分析,竇國仁得出射陽河發(fā)生閘下淤積的條件和影響因素,伍冬領(lǐng)等對(duì)永寧江樞紐引航道和擋潮閘下淤積等問題進(jìn)行了深入研究,認(rèn)為該處閘港泥沙粒徑較細(xì),易發(fā)生絮凝,曲紅嶺等采用數(shù)學(xué)模型研究了短引河河道型、短引河灘槽型和長引河河道型三種類型建閘河口閘下淤積特征,結(jié)果表明閘下淤積與閘下河道特征密切相關(guān),河口建閘后,由于潮波變形及徑流量的匱乏,浮泥隨納潮水體進(jìn)入閘下河道形成泥沙淤積,成建、劉德榮結(jié)合江蘇沿海一些擋潮閘,全面分析了閘下港道淤積泥沙來源和泥沙淤積的原因,比較全面地介紹了一些可行的治理措施,徐雪松等總結(jié)影響閘下淤積形態(tài)的多個(gè)因素,包括潮流作用、閘下引河長度、河口泥沙的動(dòng)力特性、閘下港道類型以及建閘位置等,王文杰認(rèn)為閘下河道發(fā)生持續(xù)性淤積是由水流挾沙能力和含沙量的不對(duì)稱所導(dǎo)致,楊元平等研究得出曹娥江大閘閘下淤積現(xiàn)狀及減淤對(duì)策,認(rèn)為尖山河段主槽靠近曹娥江口,曹娥江大閘閘下淤積不明顯;一旦出現(xiàn)主槽偏離南岸,閘下灘面較寬時(shí),閘下淤積問題就會(huì)出現(xiàn)Ⅲ,朱明成分析了閘下淤積的影響因素,論證了閘址條件不同時(shí)閘下淤積情況,其中閘下河道類型、閘下河道長度、徑流條件、植被條件和閘下河道人海方向都起到一定作用,錢建平等分析了小洋口閘下港槽淤積成因,提出了減淤的措施。

    以往研究提出的閘下減淤的對(duì)策有:(1)水力沖淤,將勢(shì)能轉(zhuǎn)化為動(dòng)能,利用水流的攜沙能力達(dá)到清淤的效果,水位清淤適用于閘門上下河道比降大且水源較為充足的河流,開閘放水的時(shí)機(jī)也要考慮到潮汐能量的變化,納潮沖淤指在漲潮期利用河道納潮,待有一定水頭差時(shí)開閘沖淤,一般在上游水源不足或潮汐能量強(qiáng)的地區(qū)可考慮該方法,或者,大潮期間打開閘門,增強(qiáng)納潮量和落潮流速,達(dá)到清淤效果,(2)機(jī)械清淤,機(jī)船拖耙清淤,是以機(jī)動(dòng)船拖動(dòng)耙具耙起淤泥,利用推進(jìn)器,攪動(dòng)水體,使部分耙起的淤泥變成懸沙,隨水流流走,水力挖塘機(jī)組,抓斗式或鏈斗式挖泥船亦可用于引河清淤,也可采用絞吸式、抓斗式或鏈斗式挖泥船清淤,(3)工程措施,整治閘上河道和閘下引河,擴(kuò)大排泄量;閘下裁彎、切灘,增加泄流能力;以清頂渾,對(duì)于較短的港槽可以在漲潮前2-3h,在引河中預(yù)先開閘放水,減少進(jìn)入引河的渾水,選好閘址,應(yīng)考慮沖淤水源、港口的穩(wěn)定性、閘下引河的長短等條件;圍墾時(shí)保留沖港灘面水;植物攔沙減淤,阻止灘面產(chǎn)沙、進(jìn)入引河;修筑導(dǎo)堤,束水攻沙,切斷由風(fēng)浪掀起的海灘泥沙補(bǔ)給源,減少泥沙進(jìn)入引河。

    長江口挾沙量巨大,崇明東灘鳥類保護(hù)區(qū)新建大堤外區(qū)域的淤漲明顯,海堤水閘外側(cè)引水渠極易淤積,為了提高沖淤效果,降低開閘能耗,保障調(diào)控優(yōu)化工程區(qū)生態(tài)用水安全,就需要了解引水渠自然淤積情況,以及開閘放水沖淤的效果,本文對(duì)此開展研究,給出在水閘關(guān)閉下引水渠自然淤平所需時(shí)間、開閘沖淤和納潮沖淤效果,研究成果可為崇明東灘保護(hù)區(qū)大型濕地修復(fù)工程水閘運(yùn)行提供關(guān)鍵技術(shù)支撐,對(duì)崇明東灘鳥類保護(hù)區(qū)濕地資源保育具有重要的理論意義和應(yīng)用價(jià)值。

    1區(qū)域概況

    崇明東灘位于長江口崇明島東端,是東亞一澳大利亞候鳥遷飛路線的中段,這里植被分布廣,有大量的底棲動(dòng)物,適宜水鳥停歇、越冬,是候鳥遷徙的中轉(zhuǎn)站,1998年上海市政府成立崇明東灘鳥類保護(hù)區(qū),2005年成為國家級(jí)自然保護(hù)區(qū),上海市政府以遏制互花米草擴(kuò)張,改善崇明東灘生態(tài)系統(tǒng),優(yōu)化鳥類棲息地,穩(wěn)定鳥類食物來源為目標(biāo),在崇明東灘鳥類國家級(jí)自然保護(hù)區(qū)內(nèi)開展了規(guī)模巨大的互花米草生態(tài)控制及鳥類棲息地優(yōu)化工程(以下簡稱優(yōu)化工程),該優(yōu)化工程北邊在八激閘門,南邊接崇明東灘1998大堤中部,西邊以崇明東灘1998大堤為界,東邊界為新建的海塘大堤,優(yōu)化工程范圍達(dá)240km2.于2016年12月竣工,本文研究區(qū)域?yàn)閮?yōu)化工程區(qū)c區(qū)(見圖1)。

    該優(yōu)化工程包括內(nèi)部生境區(qū)、隨塘河、大堤和涵閘,生態(tài)修復(fù)工程區(qū)的南部c區(qū),內(nèi)部生境區(qū)分為C1N、C1S、C2W、C2E、C3、C4六個(gè)不同的生境區(qū)塊(見圖1),隨塘河圍繞內(nèi)部生境區(qū)域一圈,是內(nèi)部生境與外界進(jìn)行水交換的重要水域,由于蒸發(fā)、降水、植物蒸散和土壤下滲,生境區(qū)的水深隨時(shí)間變化,當(dāng)水深不在候鳥合適的水深范圍內(nèi),就需要保護(hù)區(qū)管理人員從隨塘河引水或排水,以維持不同內(nèi)部生境區(qū)適宜水深,在調(diào)控優(yōu)化工程區(qū)新建海堤上建有水閘,當(dāng)內(nèi)部生境區(qū)和隨塘河水位低、不能滿足生態(tài)用水時(shí),需要開閘引進(jìn)閘外崇明東灘水體;當(dāng)降水量大,如梅雨季節(jié),內(nèi)部生境區(qū)和隨塘河水位超過臨界值時(shí),需要開閘放水,在一年中的大部分時(shí)間,由于降水的補(bǔ)充和生態(tài)耗水,內(nèi)部生境區(qū)和隨塘河的水位一般能滿足要求,大部分時(shí)間水閘是關(guān)閉的,長時(shí)間的水閘關(guān)閉,以及崇明東灘高泥沙濃度環(huán)境,極易導(dǎo)致閘外引水渠淤積,引水渠淤積后就難以開閘取水,需要開閘沖淤,為了提高閘下沖淤效率,優(yōu)化工程區(qū)在閘口外建立了寬50m、長300m的引水渠(見圖2),引水渠兩岸是由擋水建筑物建成的導(dǎo)堤,在高潮位時(shí)也不會(huì)被淹沒,閘門開啟期間,當(dāng)隨塘河水位低于外海時(shí),水體由引水渠自外海引入;當(dāng)隨塘河水位高于外海時(shí),水體自引水渠向外海流出,水流受到引水渠導(dǎo)堤的約束,動(dòng)能相對(duì)集中,作用在底床泥沙的時(shí)間增長,底切應(yīng)力增大,可提高沖刷效率。

    2數(shù)值模式建立和驗(yàn)證

    2.1水動(dòng)力數(shù)值模式

    水動(dòng)力模式采用作者課題組長期改進(jìn)和應(yīng)用的三維數(shù)值模式ECOM-si(semi-implicit Estuarine,Costal and Ocean Model),該模式是基于海洋數(shù)值模式POM(Princeton Ocean Model)的基礎(chǔ)上研發(fā)的可設(shè)置包括潮流、徑流、風(fēng)、波浪等各種動(dòng)力因子,改進(jìn)的ECOM-si模式在長江口水動(dòng)力過程和鹽水入侵等方面研究中得到廣泛的應(yīng)用,取得了大量研究成果,為了節(jié)省篇幅,有關(guān)模式介紹可參考已發(fā)表的文獻(xiàn)。

    2.2泥沙數(shù)值模式

    三維泥沙數(shù)值模式是在上述水動(dòng)力模式基礎(chǔ)上耦合泥沙模塊建立起來的,本文中泥沙模塊運(yùn)用是基于劉高峰改進(jìn)后的泥沙模塊,該模式考慮了泥沙的沉降過程和底部臨界啟動(dòng)應(yīng)力,具體模塊介紹和參數(shù)設(shè)置參照文獻(xiàn),這里簡單介紹泥沙控制方程,懸沙輸運(yùn)模塊控制方程為(4)式和(5)式中,τ為模式計(jì)算所得的底部切應(yīng)力,Te為臨界起動(dòng)應(yīng)力,Td為臨界沉降應(yīng)力,M為底床沖刷系數(shù),表示了水底沖刷的難易程度,取值范圍為1×10-5-4×10-3kg/(m2.a),α為泥沙發(fā)生沉降的概率,根據(jù)經(jīng)驗(yàn)一般取值為0.67-0.84。

    波浪與潮流共同作用下的底切應(yīng)力參考Liang and Li中給出的計(jì)算公式,泥沙沉速使用的是吳宇帆…經(jīng)過多次沉降實(shí)驗(yàn)建立的絮凝沉降公式,懸沙粒徑的設(shè)定決定了靜水沉速的計(jì)算結(jié)果,對(duì)于三維水動(dòng)力泥沙數(shù)值模式的模擬至關(guān)重要,長江口及鄰近海域懸沙組分多樣,在平面上和垂向上懸沙粒徑變化明顯,需依靠實(shí)測數(shù)據(jù)確定該海域懸沙粒徑,劉紅等在對(duì)崇明東灘懸浮泥沙進(jìn)行采集分析后指出,崇明東灘護(hù)花米草和泥灘表層沉積物較細(xì),懸浮泥沙中值粒徑為8-16um,自岸向海逐漸變粗,模式中泥沙中值粒徑統(tǒng)一取15um,詳細(xì)的長江口泥沙絮凝沉降公式可參見文獻(xiàn)。

    計(jì)算底部泥沙通量的關(guān)鍵參數(shù)是臨界起動(dòng)應(yīng)力,臨界起動(dòng)應(yīng)力受到多種因素影響,難以使用具體公式計(jì)算,通過實(shí)測含沙量進(jìn)行調(diào)試率定,設(shè)置經(jīng)驗(yàn)值是目前最有效的計(jì)算方法,本文根據(jù)臨界起動(dòng)應(yīng)力,利用經(jīng)驗(yàn)公式來換算臨界沉降應(yīng)力,即Td=4/9Te。

    河床底部高程的變化直接受到底部泥沙凈通量的影響,泥沙凈通量為負(fù)時(shí),河床發(fā)生沖刷,河床底高程減小;泥沙凈通量為正時(shí),發(fā)生淤積,河床底高程增大,引水渠內(nèi)局地高程變化Ah的計(jì)算公式為

    2.3計(jì)算區(qū)域和邊界條件

    2.3.1引水渠

    新建大堤水閘外側(cè)引水渠淤積和沖刷是本文研究的要點(diǎn),引水渠西邊接涵閘閘口,東邊直通崇明東灘,導(dǎo)堤兩邊為優(yōu)化工程建設(shè)擋水的堤壩,最高點(diǎn)為5.0m,外海潮水基本不會(huì)漫過導(dǎo)堤,因此可以將其作為固邊界設(shè)定(見圖2),使用Delft3D網(wǎng)格生成軟件設(shè)計(jì)的閘下引水渠正交化網(wǎng)格精度為2m(見圖3)。

    水表面風(fēng)應(yīng)力由長江口崇明東灘實(shí)測月平均風(fēng)速和風(fēng)向給出_28J,劉紅在對(duì)崇明東灘表層沉積物進(jìn)行采集分析后指出,崇明東灘護(hù)花米草和泥灘表層沉積物粒徑由岸向海逐漸增大,靠近外海處的中值粒徑為33~84um,本文取80um,本文依據(jù)淤泥質(zhì)海灘,泥沙臨界起動(dòng)應(yīng)力取值為0.05N/m2。

    閘下引水渠側(cè)向開邊界為西側(cè)的閘門和東側(cè)的崇明東灘,水閘閉合時(shí)為固邊界,打開時(shí)為開邊界,因閘門內(nèi)外水位不連續(xù)和網(wǎng)格不同,過閘水通量由堰流公式計(jì)算,引水渠東側(cè)開邊界由長江口大區(qū)域模式計(jì)算所得的水位、鹽度和泥沙濃度提供。

    2.3.2隨塘河

    建立隨塘河模式是為了計(jì)算閘內(nèi)隨時(shí)間變化的水位,采用堰流公式計(jì)算可靠的過閘水通量,隨塘河的網(wǎng)格見圖4.網(wǎng)格分辨率為10m,關(guān)閘時(shí),隨塘河的流動(dòng)主要由風(fēng)應(yīng)力驅(qū)動(dòng);開閘時(shí),受過閘流量和風(fēng)應(yīng)力的共同作用。

    由閘孔出流流量計(jì)算公式可知,閘門水通量的計(jì)算同時(shí)受到閘內(nèi)水位與閘外水位的制約,閘內(nèi)水位由隨塘河模式計(jì)算,閘外水位由長江口大區(qū)域模式計(jì)算。

    2.3.3長江口

    長江口大區(qū)域數(shù)值模式為引水渠模式提供水位、鹽度和泥沙濃度開邊界條件,模式的網(wǎng)格設(shè)置見圖5.模式中外海開邊界條件、風(fēng)場、泥沙沉降速度、底部泥沙臨界起動(dòng)應(yīng)力的設(shè)置和給出,以及波浪場的計(jì)算等,可參見羅志發(fā)博士的學(xué)位論文,模式的計(jì)算區(qū)域和網(wǎng)格見圖5.包括整個(gè)長江河口、杭州灣和鄰近海域。

    隨塘河閘外海域7月閘外水位、鹽度模擬計(jì)算結(jié)果過程見圖6.由于崇明東灘地勢(shì)較高,高潮位時(shí)淹沒,低潮位時(shí)露出,此時(shí)水位即為灘地高程,隨時(shí)間不變,鹽度不連續(xù),隨塘河三維泥沙模式采用長江口大區(qū)域三維泥沙模式計(jì)算的懸浮泥沙濃度作為泥沙邊界條件,7月隨塘河閘外泥沙濃度見圖6。

    2.4沖淤驗(yàn)證

    本文三維水動(dòng)力和泥沙數(shù)值模式已做了大量驗(yàn)證,結(jié)果表明此模式能很好地模擬水文、流速、流向和鹽度的變化,也能較好地模擬泥沙濃度的變化,為了節(jié)省篇幅,有關(guān)模式驗(yàn)證等可參考文獻(xiàn)。

    本文采用2018年7月引水渠開閘前后底高程觀測資料驗(yàn)證模式,崇明東灘候鳥保護(hù)區(qū)優(yōu)化工程區(qū)的涵閘開啟調(diào)度情況為:2018年7月27日21時(shí)開啟涵閘,從外海引水,7月28日1時(shí)關(guān)閉涵閘;7月28日21時(shí)30分開啟涵閘,從外海引水,7月29日2時(shí)關(guān)閉涵閘;7月29日22時(shí)開啟涵閘,從外海引水,7月30日1時(shí)20分關(guān)閉涵閘,閘門開啟時(shí),過閘水流不斷沖刷閘下引水渠,采用三維紅外激光掃描儀于7月27日開閘前測量閘下引水渠斷面底高程,于7月30日關(guān)閘后測量閘下引水渠斷面底高程,通過反復(fù)調(diào)整模型中底部泥沙臨界啟動(dòng)應(yīng)力,確定最佳值,實(shí)測高程與模型計(jì)算數(shù)據(jù)的比較見圖7.由實(shí)測數(shù)據(jù)可得引水渠內(nèi)平均沖淤122mm,模型計(jì)算值為123mm,RMSE(Root Mean squareError)、CC(相關(guān)系數(shù)Correlation Coefficent)、ss值(skill Score)分別為0.0625m、0.74、0.8511.驗(yàn)證結(jié)果良好。

    3模擬結(jié)果和分析

    3.1閘下自然淤積

    在自然狀態(tài)下,引水渠西側(cè)閘門關(guān)閉,只有東側(cè)與崇明東灘水域相連,隨漲落潮交換水體、鹽分和泥沙,模式設(shè)置計(jì)算時(shí)段為1月1日至3月1日,引水渠東側(cè)開邊界由水位、鹽度和泥沙濃度給出,資料來自大區(qū)域模式計(jì)算的結(jié)果,引水渠初始設(shè)置為無泥沙淤積,底高程為1.0m。

    模式計(jì)算的閘下淤積厚度的分布見圖8.涵閘關(guān)閉兩個(gè)月后,引水渠靠近閘門處基本淤平,引水渠內(nèi)淤積嚴(yán)重,涵閘關(guān)閉時(shí),引水渠內(nèi)受到外海潮漲潮落的影響,高潮位時(shí)長江口攜帶大量懸浮泥沙的鹽水進(jìn)入閘下引水渠,進(jìn)入引水渠后水體動(dòng)能逐漸減小,水流攜沙能力逐漸減小,到達(dá)閘口時(shí)水流流速最小,底切應(yīng)力減小,泥沙發(fā)生沉降,難以再次起懸,此時(shí)最易發(fā)生淤積,閘下引水渠靠近外海側(cè)受到外海潮動(dòng)力的影響,水動(dòng)力泥沙易懸浮,較難沉降,此處淤積程度不高。

    由關(guān)閘后引水渠中的底高程變化可以看出,引水渠閘口淤積量最大,靠近外海淤積量較小,引水渠內(nèi)淤積程度自閘口向外海逐漸減小,同時(shí)靠近西北側(cè)一岸淤積量較大,閘前水流流速很小,幾乎為零,在異重流的作用下,閘門附近的底部泥沙向閘下輸運(yùn),閘下底部切應(yīng)力很小,泥沙無法被落潮的水流帶走,基本全部沉降下來,因此閘前淤積程度最高,在科里奧利慣性力(科氏力)的作用下,漲潮時(shí),進(jìn)入引水渠內(nèi)的水體向北偏移,北部水體挾沙能力較大;落潮時(shí),離開引水渠的水流向南偏移,南部水體挾沙能力較大,泥沙沖刷量較大,引水渠南北兩岸挾沙能力的不平衡導(dǎo)致了引水渠內(nèi)淤積厚度不同。

    3.2閘下沖淤

    閘門淤積厚度較低時(shí),沖刷效果不明顯,頻繁開啟閘門進(jìn)行沖淤耗時(shí)耗力,不利于優(yōu)化工程區(qū)的管理;閘門淤積程度高時(shí),長時(shí)間不開啟閘門則很有可能淤死,本節(jié)分別選擇自然淤積兩個(gè)月后的淤積程度進(jìn)行沖淤效果模擬,分為3種工況進(jìn)行:(1)隨塘河蓄水至3.0m,開閘放水至2.0m;(2)隨塘河蓄水至3.8m,開閘放水至2.0m;(3)隨塘河蓄水至3.0m,選擇大潮時(shí)期開閘一天,本節(jié)分別模擬和分析自然淤積兩個(gè)月情況下3種不同工況的沖淤效果。

    模式初始底高程為淤積兩個(gè)月后的引水渠內(nèi)底高程,沖淤前引水渠底高程圖見圖8.隨塘河蓄水至3.0m,放水至2.0m,開閘放水持續(xù)時(shí)間長達(dá)20h后的沖淤效果見圖9.引水渠內(nèi)平均總沖淤厚度為96mm,閘前沖刷效果最好,這是因?yàn)榭拷l口處流態(tài)復(fù)雜,流線急劇彎曲,湍流動(dòng)能大,底床沖刷效果明顯,由于淤積時(shí)間長、引水渠內(nèi)淤積程度高,在相同的水力條件下,水深較小,底部切應(yīng)力較大,沖淤效果明顯。

    隨塘河蓄水至3.8m,放水至2.0m,開閘放水時(shí)長達(dá)22h后的沖淤效果見圖10.引水渠閘前沖刷效果較好,渠內(nèi)平均總沖淤厚度為133mm。

    隨塘河蓄水至3.0m,外海高水位進(jìn)水,開閘一天后的納潮沖淤效果見圖11.納潮沖淤時(shí),引水渠內(nèi)受到外海潮動(dòng)力的作用,引水渠內(nèi)整體沖淤效果較好,渠內(nèi)平均總沖淤厚度為625mm,閘前沖淤效果比較好,在沖淤的模式驗(yàn)證中,3d的平均沖淤幅度為122mm,遠(yuǎn)比這里的平均沖淤量值小,原因是:(1)模式驗(yàn)證時(shí)涵閘只在夜間漲潮時(shí)開啟,3d的開閘時(shí)間總和約為12h,而這里開閘時(shí)長為24h,少了一半,(2)驗(yàn)證時(shí)段前期,臺(tái)風(fēng)安比2018年7月22日12時(shí)30分在上海崇明登陸,大量降雨使得保護(hù)區(qū)水閘打開排澇,沖刷了引水渠,模式驗(yàn)證時(shí)引水渠初始淤積高程較低,導(dǎo)致模式驗(yàn)證時(shí)總的沖淤幅度較小,(3)這里納潮沖淤選擇的外邊界水位是7月最高潮位的一天,比驗(yàn)證計(jì)算中的外邊界水位更高。

    4結(jié)論

    本文在河口水動(dòng)力模式的基礎(chǔ)上,耦合泥沙模塊和底部沖淤方程,建立了崇明東灘鳥類保護(hù)區(qū)海堤外側(cè)引水渠沖淤三維數(shù)值模式,模式上游閘門水通量由堰流公式計(jì)算,內(nèi)側(cè)水位由隨塘河水動(dòng)力模式計(jì)算,模式下游水位、鹽度、泥沙濃度由長江河口大區(qū)域數(shù)值模式計(jì)算,利用2018年7月2830日優(yōu)化工程區(qū)水閘開閘沖淤前后實(shí)測引水渠底部高程,模式計(jì)算的高程與實(shí)測值吻合良好,表明模式能較好地模擬引水渠沖淤變化。

    模式計(jì)算結(jié)果表明,引水渠在自然狀態(tài)下關(guān)閉60d后閘口處基本淤平,無法使用,涵閘閘前淤積厚度約達(dá)2.0m,引水渠靠近外海處淤積厚度約為1.0m,需要開閘沖淤,淤積厚度較高時(shí),沖淤效果明顯,僅在隨塘河水力沖淤、無外海潮動(dòng)力作用時(shí),閘前沖淤效果明顯,靠近外海處沖淤厚度幾乎為0.

    引水渠自然淤積兩個(gè)月后,隨塘河蓄水至3.0m、放水至2.0m,引水渠內(nèi)平均沖淤厚度為96mm;隨塘河蓄水至3.8m、放水至2.0m,引水渠內(nèi)平均沖淤厚度為133mm,引水渠內(nèi)淤積程度高,開閘沖淤時(shí),流量相同的情況下,水深較淺,流速較大,沖淤量也因此增加,納潮沖淤后,引水渠內(nèi)整體沖淤效果好,平均沖淤厚度為625mm。

    分析比較不同工況下引水渠內(nèi)的沖淤厚度,淤積程度越高,沖淤效果越明顯,利用隨塘河蓄水沖淤時(shí),閘前沖淤效果好,總體沖淤厚度較小,水位下降快,水量損失較大,不利于維護(hù)優(yōu)化工程區(qū)閘內(nèi)生態(tài)系統(tǒng)的穩(wěn)定,在高潮位納潮沖淤時(shí),水動(dòng)力強(qiáng),閘下引水渠整體沖淤效果明顯,引水渠內(nèi)淤積的泥沙基本清除,此時(shí)隨塘河內(nèi)水量損失較少,有利于優(yōu)化工程區(qū)內(nèi)鳥類棲息地的建設(shè)與修復(fù),本文研究成果可為崇明東灘國家級(jí)自然保護(hù)區(qū)生態(tài)修復(fù)工程區(qū)的保護(hù)與管理、引水渠沖淤方案提供科技指導(dǎo),

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    亚洲成人免费电影在线观看| 两个人视频免费观看高清| 69av精品久久久久久| 99在线视频只有这里精品首页| 亚洲人成伊人成综合网2020| 久久久久久亚洲精品国产蜜桃av| 可以免费在线观看a视频的电影网站| www国产在线视频色| 国产精品 欧美亚洲| 日本一本二区三区精品| 高潮久久久久久久久久久不卡| 亚洲国产日韩欧美精品在线观看 | 欧美乱色亚洲激情| 99国产精品一区二区三区| 亚洲欧美精品综合一区二区三区| 成人特级黄色片久久久久久久| 91字幕亚洲| av在线播放免费不卡| 国产精品二区激情视频| 视频在线观看一区二区三区| 日韩欧美在线二视频| 少妇熟女aⅴ在线视频| 久久久国产欧美日韩av| 国产黄a三级三级三级人| 国产亚洲av嫩草精品影院| 精品久久蜜臀av无| 国产区一区二久久| 色哟哟哟哟哟哟| 婷婷精品国产亚洲av在线| 国产精品二区激情视频| 国产高清有码在线观看视频 | 亚洲精品在线观看二区| 国产又爽黄色视频| 88av欧美| 国产成人精品无人区| tocl精华| 欧美一级毛片孕妇| 国产成人一区二区三区免费视频网站| 法律面前人人平等表现在哪些方面| 99久久精品国产亚洲精品| 日本三级黄在线观看| 成熟少妇高潮喷水视频| 日韩欧美在线二视频| 免费一级毛片在线播放高清视频| 又黄又粗又硬又大视频| 日本成人三级电影网站| 天天躁夜夜躁狠狠躁躁| 成人三级黄色视频| 欧美黄色片欧美黄色片| 黑丝袜美女国产一区| 精品久久久久久久末码| 日本一区二区免费在线视频| 亚洲精品国产一区二区精华液| 亚洲精品av麻豆狂野| 在线观看舔阴道视频| 身体一侧抽搐| 日韩视频一区二区在线观看| aaaaa片日本免费| 在线观看舔阴道视频| 国产亚洲精品第一综合不卡| 叶爱在线成人免费视频播放| 熟女少妇亚洲综合色aaa.| 少妇熟女aⅴ在线视频| 日韩 欧美 亚洲 中文字幕| 搡老妇女老女人老熟妇| 亚洲av电影不卡..在线观看| 亚洲熟女毛片儿| 国产精品免费一区二区三区在线| 亚洲一卡2卡3卡4卡5卡精品中文| netflix在线观看网站| 亚洲成av片中文字幕在线观看| 一个人观看的视频www高清免费观看 | 欧美亚洲日本最大视频资源| 亚洲av电影在线进入| 非洲黑人性xxxx精品又粗又长| 波多野结衣巨乳人妻| 精品人妻1区二区| 黄色视频,在线免费观看| 国产日本99.免费观看| 久久精品国产综合久久久| 2021天堂中文幕一二区在线观 | 白带黄色成豆腐渣| 久久天堂一区二区三区四区| 亚洲人成电影免费在线| 欧美 亚洲 国产 日韩一| 成人三级黄色视频| 十八禁人妻一区二区| 亚洲欧美日韩高清在线视频| 波多野结衣高清无吗| 国内毛片毛片毛片毛片毛片| 国产久久久一区二区三区| 日本成人三级电影网站| 视频区欧美日本亚洲| 国内精品久久久久久久电影| 曰老女人黄片| 日韩欧美三级三区| 免费在线观看成人毛片| 欧美不卡视频在线免费观看 | 午夜福利视频1000在线观看| 中文字幕高清在线视频| 国产又色又爽无遮挡免费看| 一本一本综合久久| 成人18禁高潮啪啪吃奶动态图| 国产激情偷乱视频一区二区| 亚洲av第一区精品v没综合| 久久精品影院6| 精品乱码久久久久久99久播| 国产av一区二区精品久久| 国产私拍福利视频在线观看| 亚洲av电影在线进入| 亚洲五月婷婷丁香| 久久草成人影院| 我的亚洲天堂| 欧美不卡视频在线免费观看 | 国产成人欧美| 久久婷婷人人爽人人干人人爱| 免费看美女性在线毛片视频| 国产麻豆成人av免费视频| 一级毛片女人18水好多| 国产精品综合久久久久久久免费| 免费av毛片视频| 99久久无色码亚洲精品果冻| 伦理电影免费视频| 深夜精品福利| 操出白浆在线播放| 手机成人av网站| 国产91精品成人一区二区三区| 欧美人与性动交α欧美精品济南到| 精品高清国产在线一区| 日本 av在线| 一本久久中文字幕| 91在线观看av| 日韩欧美三级三区| 精品一区二区三区四区五区乱码| 国产精品1区2区在线观看.| 在线免费观看的www视频| 久久香蕉激情| 久久久久国产精品人妻aⅴ院| 精品国产美女av久久久久小说| 香蕉国产在线看| 国产成人精品久久二区二区免费| 97人妻精品一区二区三区麻豆 | 成人永久免费在线观看视频| 久久精品aⅴ一区二区三区四区| 午夜免费激情av| 啦啦啦韩国在线观看视频| 人人妻,人人澡人人爽秒播| 亚洲七黄色美女视频| 高清在线国产一区| 国产av在哪里看| 最近最新中文字幕大全免费视频| 一级毛片女人18水好多| 亚洲国产毛片av蜜桃av| 最近最新免费中文字幕在线| 国产成人精品久久二区二区免费| 日韩欧美一区二区三区在线观看| 午夜老司机福利片| 变态另类成人亚洲欧美熟女| 90打野战视频偷拍视频| 欧美另类亚洲清纯唯美| 成人国产一区最新在线观看| 欧美黑人巨大hd| 黄片播放在线免费| 一级a爱片免费观看的视频| 又大又爽又粗| 国产熟女xx| 久久精品aⅴ一区二区三区四区| 好男人在线观看高清免费视频 | 国产精品 国内视频| 淫妇啪啪啪对白视频| 欧美日本亚洲视频在线播放| 禁无遮挡网站| 成人国产综合亚洲| 天堂√8在线中文| 国产97色在线日韩免费| 亚洲电影在线观看av| 国产精品爽爽va在线观看网站 | 亚洲 欧美 日韩 在线 免费| 国产精品,欧美在线| 禁无遮挡网站| 欧美性长视频在线观看| 国产一区二区在线av高清观看| x7x7x7水蜜桃| 亚洲 欧美一区二区三区| www.自偷自拍.com| 亚洲精华国产精华精| 成人一区二区视频在线观看| 国产精品久久久人人做人人爽| 女警被强在线播放| 亚洲av五月六月丁香网| 极品教师在线免费播放| 正在播放国产对白刺激| 村上凉子中文字幕在线| 午夜精品在线福利| 久久久久国产精品人妻aⅴ院| 十八禁网站免费在线| 国产区一区二久久| 国产精品综合久久久久久久免费| 日韩av在线大香蕉| 国产1区2区3区精品| 啦啦啦免费观看视频1| 制服诱惑二区| 色老头精品视频在线观看| 国产精品一区二区三区四区久久 | 给我免费播放毛片高清在线观看| 久久久久久久久中文| 久久精品国产亚洲av香蕉五月| 99在线视频只有这里精品首页| 色婷婷久久久亚洲欧美| 视频区欧美日本亚洲| 久久伊人香网站| 午夜视频精品福利| 久久久久国内视频| 日本熟妇午夜| 亚洲人成77777在线视频| 日本三级黄在线观看| 嫩草影视91久久| 国产精品 国内视频| 国产成人精品久久二区二区免费| 99精品欧美一区二区三区四区| www国产在线视频色| 999精品在线视频| 99国产精品一区二区蜜桃av| netflix在线观看网站| 18禁国产床啪视频网站| 人妻久久中文字幕网| 午夜福利高清视频| 嫁个100分男人电影在线观看| 美国免费a级毛片| www日本黄色视频网| 中文字幕人妻熟女乱码| 手机成人av网站| 成人亚洲精品一区在线观看| 黑人欧美特级aaaaaa片| avwww免费| 久久久精品国产亚洲av高清涩受| 妹子高潮喷水视频| 中亚洲国语对白在线视频| 亚洲av成人一区二区三| av欧美777| 别揉我奶头~嗯~啊~动态视频| 午夜激情av网站| 午夜福利免费观看在线| 夜夜躁狠狠躁天天躁| 中文字幕人成人乱码亚洲影| 日本撒尿小便嘘嘘汇集6| 久久国产精品影院| 最近最新中文字幕大全免费视频| 日日摸夜夜添夜夜添小说| 亚洲精品久久成人aⅴ小说| 精品一区二区三区四区五区乱码| 12—13女人毛片做爰片一| 亚洲九九香蕉| 看片在线看免费视频| 久久精品国产清高在天天线| 欧美一级毛片孕妇| 亚洲专区中文字幕在线| 久9热在线精品视频| 国产黄片美女视频| 国产一区二区在线av高清观看| 日本撒尿小便嘘嘘汇集6| 久久这里只有精品19| 国产不卡一卡二| 欧美午夜高清在线| 久久久久久国产a免费观看| 久久这里只有精品19| 中文资源天堂在线| 日韩 欧美 亚洲 中文字幕| 亚洲精华国产精华精| √禁漫天堂资源中文www| 一级a爱视频在线免费观看| 国产成人欧美| 久久中文字幕一级| 成人三级做爰电影| 日韩精品中文字幕看吧| 男女做爰动态图高潮gif福利片| 免费人成视频x8x8入口观看| 国内久久婷婷六月综合欲色啪| ponron亚洲| 天堂√8在线中文| a级毛片a级免费在线| 最新在线观看一区二区三区| 一本大道久久a久久精品| 国产精品二区激情视频| 性欧美人与动物交配| 中文字幕精品免费在线观看视频| 在线视频色国产色| 日本免费a在线| 法律面前人人平等表现在哪些方面| 手机成人av网站| 欧美黄色片欧美黄色片| 午夜激情福利司机影院| 欧美日本视频| 久久婷婷人人爽人人干人人爱| 精品欧美国产一区二区三| 免费一级毛片在线播放高清视频| 午夜影院日韩av| cao死你这个sao货| 精品第一国产精品| 日韩免费av在线播放| 最近最新中文字幕大全免费视频| 级片在线观看| 国产爱豆传媒在线观看 | 国产亚洲欧美在线一区二区| 国产精品美女特级片免费视频播放器 | 久久久精品欧美日韩精品| 人人妻人人澡人人看| 免费在线观看视频国产中文字幕亚洲| 成人三级黄色视频| 久久久久久久精品吃奶| 亚洲精品国产一区二区精华液| 亚洲一区中文字幕在线| 1024香蕉在线观看| 欧美人与性动交α欧美精品济南到| 99久久99久久久精品蜜桃| www日本黄色视频网| 18美女黄网站色大片免费观看| 久久青草综合色| 久久久久亚洲av毛片大全| 大型av网站在线播放| 国产成人av激情在线播放| 人人妻,人人澡人人爽秒播| 老鸭窝网址在线观看| 欧美最黄视频在线播放免费| 亚洲成人精品中文字幕电影| 国产av又大| 好看av亚洲va欧美ⅴa在| 国产精品久久久人人做人人爽| 成熟少妇高潮喷水视频| 日韩 欧美 亚洲 中文字幕| 成人国语在线视频| 亚洲第一av免费看| 欧美成人免费av一区二区三区| 亚洲成人精品中文字幕电影| 欧美色欧美亚洲另类二区| www.自偷自拍.com| 三级毛片av免费| 成人av一区二区三区在线看| 99在线视频只有这里精品首页| 精品国产超薄肉色丝袜足j| 高清毛片免费观看视频网站| 午夜两性在线视频| 两性夫妻黄色片| 日本黄色视频三级网站网址| 桃红色精品国产亚洲av| 亚洲国产高清在线一区二区三 | 好看av亚洲va欧美ⅴa在| 久久久精品欧美日韩精品| 男男h啪啪无遮挡| 91国产中文字幕| 精品福利观看| 一边摸一边抽搐一进一小说| 中文亚洲av片在线观看爽| 又大又爽又粗| 999久久久精品免费观看国产| 午夜福利视频1000在线观看| 老汉色∧v一级毛片| 午夜福利视频1000在线观看| 欧美亚洲日本最大视频资源| 日本一区二区免费在线视频| 精品人妻1区二区| 精品一区二区三区四区五区乱码| 最好的美女福利视频网| 可以免费在线观看a视频的电影网站| 91字幕亚洲| 男人操女人黄网站| 亚洲精品中文字幕在线视频| 欧美人与性动交α欧美精品济南到| 黄色 视频免费看| 午夜精品在线福利| 日韩免费av在线播放| 国产极品粉嫩免费观看在线| 此物有八面人人有两片| 欧美黄色淫秽网站| 男人舔女人下体高潮全视频| 国产高清videossex| 国产高清激情床上av| 欧美日韩亚洲国产一区二区在线观看| 一a级毛片在线观看| 亚洲avbb在线观看| 国产精品久久久人人做人人爽| 黄色丝袜av网址大全| 91麻豆精品激情在线观看国产| 国产久久久一区二区三区| 久久99热这里只有精品18| 国产精品永久免费网站| 99热只有精品国产| 成人一区二区视频在线观看| 国产精品爽爽va在线观看网站 | 在线天堂中文资源库| 日日爽夜夜爽网站| 99国产综合亚洲精品| 一本综合久久免费| 免费人成视频x8x8入口观看| 亚洲人成网站在线播放欧美日韩| 好男人在线观看高清免费视频 | 精品日产1卡2卡| 精品国产一区二区三区四区第35| 日韩欧美在线二视频| 久久99热这里只有精品18| 午夜久久久久精精品| 18禁国产床啪视频网站| 日日夜夜操网爽| 欧美中文日本在线观看视频| 国产在线观看jvid| svipshipincom国产片| 久久精品国产清高在天天线| 免费在线观看日本一区| 90打野战视频偷拍视频| 欧美性长视频在线观看| 免费看十八禁软件| 老司机靠b影院| 欧美日韩福利视频一区二区| 久久人人精品亚洲av| 国产亚洲欧美98| 嫩草影视91久久| 亚洲精品久久国产高清桃花| 国内毛片毛片毛片毛片毛片| 亚洲av熟女| 日本在线视频免费播放| 成年版毛片免费区| 亚洲成人国产一区在线观看| 中文字幕人妻熟女乱码| 又紧又爽又黄一区二区| 国产黄色小视频在线观看| 久久 成人 亚洲| 国产91精品成人一区二区三区| 特大巨黑吊av在线直播 | 狠狠狠狠99中文字幕| 国产成人精品久久二区二区免费| 国产精品一区二区三区四区久久 | 动漫黄色视频在线观看| 啪啪无遮挡十八禁网站| 国产伦一二天堂av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 欧美av亚洲av综合av国产av| xxxwww97欧美| 国内揄拍国产精品人妻在线 | 亚洲自拍偷在线| 亚洲精品国产一区二区精华液| 亚洲激情在线av| 欧美久久黑人一区二区| 国产亚洲欧美在线一区二区| 日本熟妇午夜| 桃红色精品国产亚洲av| 日韩精品青青久久久久久| 女人被狂操c到高潮| 亚洲片人在线观看| 精品第一国产精品| 亚洲av电影在线进入| 色婷婷久久久亚洲欧美| 在线观看免费日韩欧美大片| 久久欧美精品欧美久久欧美| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产清高在天天线| 99国产精品99久久久久| 满18在线观看网站| 精品一区二区三区视频在线观看免费| 99热6这里只有精品| 亚洲欧美日韩无卡精品| 天堂影院成人在线观看| 日韩三级视频一区二区三区| 91在线观看av| 最近最新免费中文字幕在线| 国产精品综合久久久久久久免费| 可以免费在线观看a视频的电影网站| 久久精品国产亚洲av香蕉五月| 50天的宝宝边吃奶边哭怎么回事| or卡值多少钱| 欧美不卡视频在线免费观看 | 99精品在免费线老司机午夜| 日韩视频一区二区在线观看| 亚洲欧美精品综合一区二区三区| 午夜精品久久久久久毛片777| 国产亚洲精品第一综合不卡| 精品久久久久久久久久久久久 | 制服人妻中文乱码| 色综合站精品国产| 日本熟妇午夜| 香蕉国产在线看| 亚洲精品一区av在线观看| 午夜福利成人在线免费观看| 亚洲熟妇中文字幕五十中出| 好看av亚洲va欧美ⅴa在| 亚洲五月天丁香| 免费在线观看影片大全网站| 免费看日本二区| 宅男免费午夜| 国产亚洲精品av在线| 欧美黑人巨大hd| 少妇粗大呻吟视频| 色综合欧美亚洲国产小说| 亚洲av美国av| 成人永久免费在线观看视频| 桃红色精品国产亚洲av| 欧美又色又爽又黄视频| 国产精品永久免费网站| 欧美性猛交黑人性爽| 老熟妇乱子伦视频在线观看| 亚洲色图 男人天堂 中文字幕| 久久久久亚洲av毛片大全| 欧洲精品卡2卡3卡4卡5卡区| 欧美激情极品国产一区二区三区| 少妇熟女aⅴ在线视频| 18美女黄网站色大片免费观看| 国产99白浆流出| 亚洲成人久久爱视频| 亚洲成人精品中文字幕电影| 日韩欧美国产在线观看| 欧美激情极品国产一区二区三区| 成人永久免费在线观看视频| 色精品久久人妻99蜜桃| 久久婷婷成人综合色麻豆| 国产麻豆成人av免费视频| 国产又爽黄色视频| 欧美久久黑人一区二区| 91成年电影在线观看| 欧美日韩瑟瑟在线播放| 久久久久国内视频| 此物有八面人人有两片| 成熟少妇高潮喷水视频| 久久久久久久久久黄片| 精品久久久久久久人妻蜜臀av| 国产亚洲精品久久久久5区| 婷婷精品国产亚洲av| 亚洲色图 男人天堂 中文字幕| 久久久久亚洲av毛片大全| 国产亚洲欧美精品永久| 天堂动漫精品| 亚洲一区二区三区不卡视频| 国产成人影院久久av| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看亚洲国产| 日韩精品免费视频一区二区三区| 禁无遮挡网站| 精品人妻1区二区| 国产人伦9x9x在线观看| or卡值多少钱| 亚洲一区中文字幕在线| 欧美成狂野欧美在线观看| 波多野结衣高清无吗| 校园春色视频在线观看| 天堂动漫精品| 欧美性长视频在线观看| 亚洲五月婷婷丁香| 午夜免费成人在线视频| 国产精品永久免费网站| 美女高潮到喷水免费观看| 国产精品国产高清国产av| 最近最新免费中文字幕在线| 午夜福利高清视频| 亚洲男人的天堂狠狠| 国产精品二区激情视频| av电影中文网址| 国产成人精品久久二区二区91| 日本黄色视频三级网站网址| 99久久精品国产亚洲精品| 1024香蕉在线观看| 亚洲av片天天在线观看| 国产精品 国内视频| 老熟妇仑乱视频hdxx| 精品午夜福利视频在线观看一区| 亚洲,欧美精品.| netflix在线观看网站| 日本 欧美在线| 此物有八面人人有两片| 久久伊人香网站| 女生性感内裤真人,穿戴方法视频| 99国产精品一区二区三区| 国产精华一区二区三区| 欧美日本亚洲视频在线播放| 日韩欧美国产在线观看| 国产一区二区三区在线臀色熟女| 手机成人av网站| 女警被强在线播放| 国产精品香港三级国产av潘金莲| 精品久久久久久成人av| 丰满的人妻完整版| 国产单亲对白刺激| 日日夜夜操网爽| 亚洲精品美女久久av网站| 国产在线精品亚洲第一网站| 亚洲男人天堂网一区| 亚洲欧美日韩高清在线视频| 亚洲,欧美精品.| 两个人看的免费小视频| 啦啦啦 在线观看视频| 精品高清国产在线一区| 丝袜美腿诱惑在线| 欧美不卡视频在线免费观看 | 亚洲av熟女| 精品不卡国产一区二区三区| 国产色视频综合| 国产高清视频在线播放一区| 日本撒尿小便嘘嘘汇集6| 久久国产精品男人的天堂亚洲| 国内久久婷婷六月综合欲色啪| 天堂√8在线中文| 一本综合久久免费| 亚洲精华国产精华精| 少妇粗大呻吟视频| 国产97色在线日韩免费| 午夜福利一区二区在线看| 亚洲 欧美一区二区三区| 91av网站免费观看| 久热爱精品视频在线9| 男人操女人黄网站| 久久久国产欧美日韩av| 久久久久精品国产欧美久久久| 长腿黑丝高跟| 成人永久免费在线观看视频| 1024视频免费在线观看| 日韩免费av在线播放| 欧美成人免费av一区二区三区| 亚洲男人天堂网一区|