• <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
    久久精品aⅴ一区二区三区四区| 噜噜噜噜噜久久久久久91| 99精品在免费线老司机午夜| 国产精品久久电影中文字幕| 亚洲人成电影免费在线| www.www免费av| 真人一进一出gif抽搐免费| 色精品久久人妻99蜜桃| 叶爱在线成人免费视频播放| 亚洲精品美女久久av网站| 亚洲av日韩精品久久久久久密| 欧美性猛交黑人性爽| 男人舔女人下体高潮全视频| 国产午夜精品久久久久久| 久久久水蜜桃国产精品网| 国产精品九九99| 国产精品精品国产色婷婷| 亚洲av电影在线进入| www.999成人在线观看| 日韩欧美三级三区| 美女免费视频网站| 午夜影院日韩av| 最新中文字幕久久久久 | 午夜福利欧美成人| 久久久色成人| tocl精华| 色视频www国产| www.精华液| 亚洲 欧美一区二区三区| 熟妇人妻久久中文字幕3abv| 成人永久免费在线观看视频| 一级毛片女人18水好多| 高潮久久久久久久久久久不卡| 在线观看舔阴道视频| 麻豆成人av在线观看| 亚洲欧美精品综合一区二区三区| 亚洲色图av天堂| 两性午夜刺激爽爽歪歪视频在线观看| 国产黄片美女视频| 亚洲欧美日韩东京热| 日韩精品中文字幕看吧| 国产免费男女视频| 成年版毛片免费区| 精品福利观看| 国产精品1区2区在线观看.| 国产精品亚洲美女久久久| 欧美绝顶高潮抽搐喷水| 全区人妻精品视频| 身体一侧抽搐| 男女下面进入的视频免费午夜| 久久99热这里只有精品18| 99精品欧美一区二区三区四区| 亚洲无线在线观看| 我要搜黄色片| 亚洲性夜色夜夜综合| 男女视频在线观看网站免费| 国产成人影院久久av| 俺也久久电影网| 欧美日韩一级在线毛片| 精品日产1卡2卡| 看片在线看免费视频| 亚洲av成人一区二区三| 一级毛片高清免费大全| 精品一区二区三区av网在线观看| 亚洲在线观看片| 中文字幕高清在线视频| 岛国在线观看网站| 全区人妻精品视频| 日韩免费av在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 国产黄a三级三级三级人| 一个人观看的视频www高清免费观看 | 九九久久精品国产亚洲av麻豆 | 亚洲九九香蕉| 俺也久久电影网| 99热6这里只有精品| 99久国产av精品| 高清在线国产一区| 亚洲av中文字字幕乱码综合| 人人妻人人看人人澡| 国产精品99久久99久久久不卡| 国产精品久久久久久久电影 | 好男人电影高清在线观看| 久久国产乱子伦精品免费另类| 久久99热这里只有精品18| 日本免费a在线| 日本一二三区视频观看| 特大巨黑吊av在线直播| 欧美中文日本在线观看视频| 一本精品99久久精品77| 真实男女啪啪啪动态图| 国产精品一区二区三区四区免费观看 | 最近视频中文字幕2019在线8| 国产精品女同一区二区软件 | tocl精华| 欧美一级毛片孕妇| 欧美大码av| 少妇熟女aⅴ在线视频| 97碰自拍视频| 少妇丰满av| 观看美女的网站| 久久久色成人| 后天国语完整版免费观看| 国产精品永久免费网站| 黑人欧美特级aaaaaa片| 精品久久久久久久久久免费视频| 亚洲av成人精品一区久久| 一级毛片精品| 亚洲国产日韩欧美精品在线观看 | 精品国产亚洲在线| 欧美+亚洲+日韩+国产| 亚洲av免费在线观看| 欧美在线一区亚洲| 久久这里只有精品19| 九色国产91popny在线| 9191精品国产免费久久| 精品国内亚洲2022精品成人| 国产成人精品久久二区二区免费| 日本黄大片高清| 亚洲美女视频黄频| 极品教师在线免费播放| 非洲黑人性xxxx精品又粗又长| 麻豆成人午夜福利视频| 真人一进一出gif抽搐免费| 桃色一区二区三区在线观看| 亚洲成a人片在线一区二区| 两个人看的免费小视频| 一本综合久久免费| www.自偷自拍.com| 国产伦精品一区二区三区视频9 | 婷婷精品国产亚洲av| 亚洲精品一区av在线观看| 午夜免费观看网址| 久久久久久久精品吃奶| 狂野欧美激情性xxxx| 人妻丰满熟妇av一区二区三区| 9191精品国产免费久久| 熟女电影av网| 日韩成人在线观看一区二区三区| 亚洲精品美女久久久久99蜜臀| 精品免费久久久久久久清纯| 成人午夜高清在线视频| 操出白浆在线播放| 日本在线视频免费播放| 99久久99久久久精品蜜桃| 黄色成人免费大全| www.自偷自拍.com| 国产欧美日韩一区二区三| 久久中文字幕一级| 男插女下体视频免费在线播放| 免费看日本二区| 国产aⅴ精品一区二区三区波| 国产精品亚洲美女久久久| 久久香蕉精品热| 黄色 视频免费看| 国产主播在线观看一区二区| 桃红色精品国产亚洲av| 国产精品亚洲美女久久久| 国产精品久久久久久精品电影| 后天国语完整版免费观看| 91在线观看av| 在线免费观看的www视频| 99久久国产精品久久久| 欧美精品啪啪一区二区三区| 精品乱码久久久久久99久播| 成年女人毛片免费观看观看9| 变态另类成人亚洲欧美熟女| or卡值多少钱| 中文字幕最新亚洲高清| 成人三级黄色视频| 欧美日韩福利视频一区二区| 国产精品亚洲美女久久久| 日韩欧美一区二区三区在线观看| 桃色一区二区三区在线观看| 黑人巨大精品欧美一区二区mp4| 麻豆国产97在线/欧美| 女人被狂操c到高潮| 欧美日韩综合久久久久久 | 亚洲第一电影网av| 日本与韩国留学比较| 欧美黑人巨大hd| 啦啦啦韩国在线观看视频| 亚洲乱码一区二区免费版| 午夜福利在线观看吧| 精品久久久久久久人妻蜜臀av| 国产真人三级小视频在线观看| 午夜福利视频1000在线观看| 国产精品乱码一区二三区的特点| 麻豆国产97在线/欧美| 老司机福利观看| 欧美日韩黄片免| 国产精品精品国产色婷婷| 伦理电影免费视频| 丁香欧美五月| 在线观看免费午夜福利视频| 99视频精品全部免费 在线 | 国产精品久久久久久久电影 | 国产单亲对白刺激| 一区二区三区高清视频在线| 麻豆av在线久日| 在线国产一区二区在线| 精品国产美女av久久久久小说| 99久久国产精品久久久| 国产极品精品免费视频能看的| 亚洲狠狠婷婷综合久久图片| 久久久精品大字幕| 一区福利在线观看| 黄色片一级片一级黄色片| 老熟妇乱子伦视频在线观看| 国产一区二区三区在线臀色熟女| 99热这里只有是精品50| 日韩成人在线观看一区二区三区| 一进一出抽搐动态| 日韩中文字幕欧美一区二区| 淫秽高清视频在线观看| 久久久久久久久免费视频了| 亚洲乱码一区二区免费版| 日本一二三区视频观看| 又黄又粗又硬又大视频| 丰满人妻一区二区三区视频av | 久久伊人香网站| 伊人久久大香线蕉亚洲五| 最近最新免费中文字幕在线| 一区二区三区国产精品乱码| 在线观看一区二区三区| 99久久精品国产亚洲精品| 亚洲 欧美 日韩 在线 免费| 日本 欧美在线| 国产亚洲精品综合一区在线观看| 日本精品一区二区三区蜜桃| 国产精品久久电影中文字幕| www.自偷自拍.com| 欧美性猛交黑人性爽| 男人和女人高潮做爰伦理| 香蕉久久夜色| 欧美一区二区国产精品久久精品| 高潮久久久久久久久久久不卡| 中文字幕人成人乱码亚洲影| xxxwww97欧美| 悠悠久久av| 嫩草影院入口| 日本熟妇午夜| 欧美性猛交╳xxx乱大交人| 黑人欧美特级aaaaaa片| 久久久成人免费电影| 最新在线观看一区二区三区| 大型黄色视频在线免费观看| 亚洲成av人片免费观看| 三级男女做爰猛烈吃奶摸视频| 综合色av麻豆| 丁香六月欧美| 一二三四社区在线视频社区8| 熟女少妇亚洲综合色aaa.| 日韩欧美 国产精品| 国产精品 国内视频| 亚洲国产日韩欧美精品在线观看 | 亚洲色图 男人天堂 中文字幕| 欧美激情久久久久久爽电影| 天天添夜夜摸| 亚洲av中文字字幕乱码综合| 99在线人妻在线中文字幕| 国产久久久一区二区三区| 国产成人欧美在线观看| 日本五十路高清| 给我免费播放毛片高清在线观看| 日本熟妇午夜| av在线天堂中文字幕| 中文字幕久久专区| 久久中文看片网| 成人亚洲精品av一区二区| 97碰自拍视频| bbb黄色大片| 国产精品亚洲一级av第二区| 亚洲av第一区精品v没综合| 国产伦精品一区二区三区四那| 久久99热这里只有精品18| www国产在线视频色| 好男人电影高清在线观看| 黄片大片在线免费观看| 小说图片视频综合网站| 亚洲欧美日韩卡通动漫| 亚洲最大成人中文| 亚洲欧美一区二区三区黑人| 国产精品亚洲美女久久久| 国产精品九九99| 黑人欧美特级aaaaaa片| 在线免费观看的www视频| 蜜桃久久精品国产亚洲av| 久久久国产欧美日韩av| 久99久视频精品免费| 午夜两性在线视频| 日韩高清综合在线| 国产精品1区2区在线观看.| 性色av乱码一区二区三区2| 一区二区三区国产精品乱码| 亚洲国产欧美人成| 窝窝影院91人妻| 88av欧美| 天堂网av新在线| 午夜免费观看网址| 级片在线观看| 国产私拍福利视频在线观看| 中文字幕最新亚洲高清| 亚洲激情在线av| 亚洲av电影不卡..在线观看| 五月伊人婷婷丁香| 久久久久国产一级毛片高清牌| 人人妻,人人澡人人爽秒播| 午夜精品久久久久久毛片777| 国产亚洲av高清不卡| 欧美zozozo另类| 韩国av一区二区三区四区| 亚洲精品粉嫩美女一区| 啦啦啦韩国在线观看视频| 岛国视频午夜一区免费看| 高清在线国产一区| 香蕉av资源在线| 亚洲精品乱码久久久v下载方式 | 欧美成人一区二区免费高清观看 | 99在线视频只有这里精品首页| 久久香蕉精品热| 嫩草影院入口| 又紧又爽又黄一区二区| 久久久久国内视频| 夜夜看夜夜爽夜夜摸| 麻豆av在线久日| 欧美日韩精品网址| 久久精品91无色码中文字幕| 欧美不卡视频在线免费观看| 色播亚洲综合网| 天天一区二区日本电影三级| 制服丝袜大香蕉在线| 久久精品国产综合久久久| 久久精品综合一区二区三区| 看黄色毛片网站| 亚洲精品国产精品久久久不卡| 男女那种视频在线观看| 亚洲熟妇中文字幕五十中出| 亚洲 欧美一区二区三区| 国产精品一区二区三区四区久久| 淫妇啪啪啪对白视频| 欧美午夜高清在线| 日本一本二区三区精品| 十八禁网站免费在线| 黑人欧美特级aaaaaa片| 天堂√8在线中文| 高潮久久久久久久久久久不卡| 亚洲av美国av| 久久久精品大字幕| 国产黄a三级三级三级人| 免费在线观看日本一区| 女生性感内裤真人,穿戴方法视频| 亚洲在线自拍视频| www.999成人在线观看| 国产成人av激情在线播放| 精品乱码久久久久久99久播| 国内精品久久久久精免费| 亚洲av美国av| 亚洲 欧美一区二区三区| 国产av不卡久久| 脱女人内裤的视频| 久久精品综合一区二区三区| netflix在线观看网站| 午夜福利免费观看在线| 亚洲一区二区三区不卡视频| 狠狠狠狠99中文字幕| 熟妇人妻久久中文字幕3abv| av视频在线观看入口| 免费在线观看视频国产中文字幕亚洲| 精品国产乱码久久久久久男人| 国产成+人综合+亚洲专区| 黄色视频,在线免费观看| 中文资源天堂在线| а√天堂www在线а√下载| 超碰成人久久| 啪啪无遮挡十八禁网站| 欧美日韩综合久久久久久 | 少妇熟女aⅴ在线视频| 狠狠狠狠99中文字幕| 国内精品一区二区在线观看| 最近最新免费中文字幕在线| 亚洲电影在线观看av| 欧美成狂野欧美在线观看| 999精品在线视频| av女优亚洲男人天堂 | 欧美成人免费av一区二区三区| 美女午夜性视频免费| 欧美极品一区二区三区四区| 久久婷婷人人爽人人干人人爱| svipshipincom国产片| 两个人视频免费观看高清| 欧美日韩瑟瑟在线播放| 精品日产1卡2卡| 日韩欧美 国产精品| 18禁黄网站禁片午夜丰满| 久久久久亚洲av毛片大全| 免费在线观看亚洲国产| 欧美av亚洲av综合av国产av| 亚洲精品乱码久久久v下载方式 | 久久九九热精品免费| av在线蜜桃| 国产美女午夜福利| 日韩欧美精品v在线| 国产单亲对白刺激| 99久久精品热视频| 夜夜躁狠狠躁天天躁| 国产一区二区激情短视频| 国产成人精品久久二区二区免费| 男女之事视频高清在线观看| 黄色成人免费大全| 真人一进一出gif抽搐免费| 国产爱豆传媒在线观看| 夜夜夜夜夜久久久久| 午夜免费成人在线视频| 精品久久蜜臀av无| 国产精品av视频在线免费观看| 亚洲天堂国产精品一区在线| 亚洲精品美女久久久久99蜜臀| 窝窝影院91人妻| 一个人免费在线观看的高清视频| 亚洲va日本ⅴa欧美va伊人久久| 免费人成视频x8x8入口观看| 欧美三级亚洲精品| 99久久国产精品久久久| 丁香欧美五月| 亚洲欧洲精品一区二区精品久久久| 国产精品一及| 天堂av国产一区二区熟女人妻| 天天添夜夜摸| 老司机在亚洲福利影院| 成年人黄色毛片网站| 999久久久国产精品视频| xxxwww97欧美| 亚洲美女视频黄频| 啦啦啦观看免费观看视频高清| 视频区欧美日本亚洲| 一本综合久久免费| 色精品久久人妻99蜜桃| 两个人的视频大全免费| 日本撒尿小便嘘嘘汇集6| 在线视频色国产色| 级片在线观看| 琪琪午夜伦伦电影理论片6080| 搡老熟女国产l中国老女人| 国产成人精品无人区| 变态另类成人亚洲欧美熟女| 99精品久久久久人妻精品| 国产男靠女视频免费网站| 午夜福利成人在线免费观看| 88av欧美| 韩国av一区二区三区四区| 欧美zozozo另类| 国产视频一区二区在线看| 99精品久久久久人妻精品| 成人三级做爰电影| 国产一区二区三区视频了| 久久久久九九精品影院| 一级黄色大片毛片| 嫩草影院入口| 亚洲专区国产一区二区| cao死你这个sao货| 国产免费男女视频| 亚洲av电影在线进入| 精品久久久久久久久久免费视频| 色吧在线观看| 国产精品久久电影中文字幕| 亚洲精品在线观看二区| 别揉我奶头~嗯~啊~动态视频| 小说图片视频综合网站| 最近最新中文字幕大全免费视频| 99久久99久久久精品蜜桃| 国产午夜精品久久久久久| 亚洲在线观看片| 精品乱码久久久久久99久播| 国内久久婷婷六月综合欲色啪| 午夜福利欧美成人| 老司机午夜福利在线观看视频| 99国产极品粉嫩在线观看| 亚洲人与动物交配视频| 每晚都被弄得嗷嗷叫到高潮| 18禁裸乳无遮挡免费网站照片| 老司机在亚洲福利影院| 99国产精品一区二区蜜桃av| 欧美精品啪啪一区二区三区| 男女之事视频高清在线观看| 国产亚洲精品一区二区www| 视频区欧美日本亚洲| 成人精品一区二区免费| 国产一区二区激情短视频| 国产免费男女视频| 免费观看的影片在线观看| 伊人久久大香线蕉亚洲五| 成人亚洲精品av一区二区| 亚洲18禁久久av| 99久久99久久久精品蜜桃| 久久亚洲真实| 国产乱人视频| 亚洲国产精品久久男人天堂| 亚洲精品国产精品久久久不卡| 亚洲va日本ⅴa欧美va伊人久久| 男插女下体视频免费在线播放| 久久久久亚洲av毛片大全| 国产成人啪精品午夜网站| 久久久水蜜桃国产精品网| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看的影片在线观看| 男女午夜视频在线观看| 最新美女视频免费是黄的| 国产一区二区三区在线臀色熟女| 黄色成人免费大全| 非洲黑人性xxxx精品又粗又长| 国内揄拍国产精品人妻在线| 一级毛片精品| 男女视频在线观看网站免费| 999久久久国产精品视频| 校园春色视频在线观看| 亚洲18禁久久av| 免费在线观看视频国产中文字幕亚洲| 毛片女人毛片| 丁香欧美五月| 一级毛片女人18水好多| 国产成人精品久久二区二区91| 午夜福利在线观看免费完整高清在 | 国产高清视频在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看| 欧美成狂野欧美在线观看| 精品久久久久久,| 日本一本二区三区精品| 男女之事视频高清在线观看| 999精品在线视频| 久久久久国产一级毛片高清牌| 波多野结衣高清作品| 亚洲专区字幕在线| 久久午夜亚洲精品久久| 啦啦啦免费观看视频1| 免费在线观看日本一区| 日韩av在线大香蕉| 久久久久久九九精品二区国产| 成熟少妇高潮喷水视频| 国产午夜精品久久久久久| 老熟妇乱子伦视频在线观看| 99久久国产精品久久久| 一进一出好大好爽视频| 亚洲一区二区三区不卡视频| 亚洲成av人片免费观看| 欧美在线黄色| 天堂av国产一区二区熟女人妻| 亚洲国产欧美一区二区综合| 神马国产精品三级电影在线观看| 日韩欧美国产在线观看| 日韩欧美免费精品| 麻豆久久精品国产亚洲av| 欧美一级毛片孕妇| 久久香蕉国产精品| 欧美日韩精品网址| 亚洲人成伊人成综合网2020| 国产精品免费一区二区三区在线| 国产午夜精品论理片| 男人和女人高潮做爰伦理| 国产一区二区在线观看日韩 | 曰老女人黄片| 黄色片一级片一级黄色片| 丰满的人妻完整版| 国产成+人综合+亚洲专区| 成年人黄色毛片网站| 99re在线观看精品视频| 国产av麻豆久久久久久久| 精品久久久久久久末码| 人人妻,人人澡人人爽秒播| 制服人妻中文乱码| 18禁美女被吸乳视频| ponron亚洲| 99久久成人亚洲精品观看| 久久中文字幕一级| 免费在线观看影片大全网站| 亚洲国产精品久久男人天堂| 欧美性猛交╳xxx乱大交人| 亚洲中文字幕一区二区三区有码在线看 | 日韩欧美 国产精品| 国产久久久一区二区三区| 天天添夜夜摸| 欧美日韩黄片免| 国产激情久久老熟女| 亚洲无线在线观看| 99视频精品全部免费 在线 | 美女高潮喷水抽搐中文字幕| 成人高潮视频无遮挡免费网站| 在线观看日韩欧美| 久久天堂一区二区三区四区| 国产伦一二天堂av在线观看| 久久久色成人| 婷婷六月久久综合丁香| 国产伦在线观看视频一区| 一级作爱视频免费观看| 99久久精品热视频| 熟妇人妻久久中文字幕3abv| 精品熟女少妇八av免费久了| 欧美午夜高清在线| 亚洲欧美日韩东京热| 热99re8久久精品国产| 1024香蕉在线观看| 岛国在线观看网站| 亚洲欧美激情综合另类| 亚洲av成人av| 又紧又爽又黄一区二区| 亚洲精品456在线播放app | 亚洲熟女毛片儿| 99久久精品热视频| 窝窝影院91人妻| 欧美在线一区亚洲| 欧美又色又爽又黄视频| 国产在线精品亚洲第一网站| 成人午夜高清在线视频| 午夜免费成人在线视频| 免费观看精品视频网站|