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

    三峽水庫泄水- 蓄水過程香溪河庫灣水流振蕩特性

    2022-05-10 13:21:34龍良紅紀(jì)道斌楊忠勇楊正健劉德富湯正陽
    水科學(xué)進(jìn)展 2022年2期
    關(guān)鍵詞:三峽水庫溪河支流

    徐 慧,龍良紅,2,紀(jì)道斌,2,楊忠勇,2,楊正健,2,劉德富,湯正陽

    (1. 三峽大學(xué)水利與環(huán)境學(xué)院,湖北 宜昌 443002;2. 三峽大學(xué)三峽水庫生態(tài)系統(tǒng)湖北省野外科學(xué)觀測(cè)研究站,湖北 宜昌 443002;3. 湖北工業(yè)大學(xué)河湖生態(tài)修復(fù)與藻類利用湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430068;4. 中國長(zhǎng)江三峽集團(tuán)公司,北京 100038)

    三峽水庫是三峽水電站建成后蓄水形成的人工湖泊,總面積約1 084 km2,是中國重要的戰(zhàn)略水資源庫,在長(zhǎng)江流域防洪、發(fā)電、航運(yùn)、灌溉等方面起著重要作用。目前三峽水庫大部分支流回水區(qū)水體呈富營養(yǎng)化狀態(tài),而且部分支流庫灣暴發(fā)了嚴(yán)重的藻類水華,這已成為國內(nèi)外廣泛關(guān)注的庫區(qū)水環(huán)境問題[1- 2]。針對(duì)三峽水庫的水環(huán)境問題,不少學(xué)者通過建立水流特性與水質(zhì)、水華、沉積物、溫室氣體動(dòng)態(tài)的內(nèi)在聯(lián)系[3- 6],認(rèn)為支流庫灣的水動(dòng)力過程是影響其水環(huán)境現(xiàn)狀和發(fā)展的重要物理驅(qū)動(dòng)。基于此認(rèn)識(shí),許多學(xué)者提出通過水庫生態(tài)調(diào)度來調(diào)控支流庫灣水動(dòng)力過程,進(jìn)而改善水質(zhì)、水華等問題,成為緩解水環(huán)境問題的有力措施和新途徑[7- 8]。

    水庫生態(tài)調(diào)度的本質(zhì)是通過短時(shí)間內(nèi)調(diào)控水位的波動(dòng)實(shí)現(xiàn)增大干支流間的水體交換、改變水流形態(tài)、破壞水體分層狀態(tài)等來改善庫灣水環(huán)境[9- 11]。前期關(guān)于三峽水庫流量調(diào)節(jié)的研究主要集中在三峽下游河段的生態(tài)影響,例如三峽—葛洲壩兩壩間的水位、流量波動(dòng)特征對(duì)航運(yùn)安全和發(fā)電效率等造成的影響[12- 13]。近年來,針對(duì)支流庫灣生態(tài)環(huán)境,卜英[9]研究指出日調(diào)峰運(yùn)行方式可以提高支流水體流速,加強(qiáng)干支流的水體交換和摻混,但對(duì)其影響機(jī)制缺乏合理的解釋;徐雅倩[14]進(jìn)一步研究表明水位的高頻波動(dòng)可以增加支流庫灣的水體摻混,改變其水溫分層/混合過程,進(jìn)而對(duì)水環(huán)境產(chǎn)生深遠(yuǎn)影響。但前期的研究多基于長(zhǎng)周期的流量調(diào)節(jié)尺度,針對(duì)日調(diào)節(jié)尺度下的水位高頻波動(dòng)、水動(dòng)力響應(yīng)機(jī)制研究較少。2020年,Long等[15]發(fā)現(xiàn)香溪河庫灣蓄水期存在高頻的水流振蕩現(xiàn)象,并且表明這一物理過程與水庫日調(diào)度過程密切相關(guān);2021年紀(jì)道斌等[16]進(jìn)一步報(bào)道了三峽出庫流量日調(diào)節(jié)所驅(qū)動(dòng)的干流重力波問題,解析了形成機(jī)制,并闡述了在庫區(qū)干流的傳播特征。這些發(fā)現(xiàn)不僅揭示了三峽水庫長(zhǎng)期存在但被忽視的高頻水動(dòng)力現(xiàn)象,同時(shí)也為水庫調(diào)度提供了新的研究視角。

    然而,目前的研究仍缺乏對(duì)支流庫灣振蕩物理特性的系統(tǒng)性報(bào)道。針對(duì)這一問題,本研究基于野外高頻的流速、水位監(jiān)測(cè),對(duì)泄水- 蓄水典型調(diào)度過程下支流庫灣水體振蕩特性進(jìn)行深入探討,主要包括香溪河庫灣水位- 水流特性、庫灣振蕩的物理特征參數(shù)、庫灣振蕩的生成機(jī)制及影響因素等方面的問題。通過該研究,進(jìn)一步拓寬和深入認(rèn)識(shí)水庫水動(dòng)力過程及其對(duì)三峽水庫調(diào)度過程的響應(yīng)規(guī)律,以期為后續(xù)研究和實(shí)施生態(tài)調(diào)度解決庫灣水體富營養(yǎng)化及水華問題提供理論參考。

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

    1.1 研究區(qū)域

    三峽水庫位于長(zhǎng)江流域腹心地帶,年均徑流量達(dá)4 510億m3,約占長(zhǎng)江年總徑流量的49%,壩址斷面多年平均流量14 300 m3/s。三峽工程位于長(zhǎng)江三峽西陵峽河段,正常蓄水位175 m,防洪庫容為221.5億m3,為季調(diào)節(jié)水庫。香溪河是三峽水庫庫首湖北境內(nèi)第一大支流,全長(zhǎng)94 km,發(fā)源于神農(nóng)架林區(qū),向南流經(jīng)興山縣和秭歸縣,于郭家壩鎮(zhèn)注入長(zhǎng)江。三峽水庫蓄水至145 m時(shí),距離河口28 km范圍形成回水區(qū);當(dāng)蓄水至175 m時(shí),回水區(qū)范圍長(zhǎng)達(dá)42 km(河道長(zhǎng)度估算來源于谷歌地圖)。本研究水位觀測(cè)點(diǎn)分別位于香溪河常年回水區(qū)的上游、中游、下游,流速監(jiān)測(cè)點(diǎn)分別在回水區(qū)中游峽口鎮(zhèn)和末端平邑口(見圖1),位于香溪河庫灣的水華敏感區(qū)。

    圖1 研究區(qū)域及監(jiān)測(cè)點(diǎn)位Fig.1 Map of study area and monitoring sites

    1.2 監(jiān)測(cè)指標(biāo)與方法

    本研究中所需數(shù)據(jù)主要包括香溪河庫灣水位、流速和三峽水庫壩前水位、出庫流量(三峽大壩)等。香溪河庫灣水位監(jiān)測(cè)點(diǎn)分布如圖1(b)所示,從河口至上游共3個(gè)測(cè)點(diǎn),各測(cè)點(diǎn)的位置及河道斷面等信息提取于谷歌地圖,如表1所示。流速剖面數(shù)據(jù)采用聲學(xué)多普勒流速剖面儀(Signature 1 000 ADCP,Nortek,挪威)收集,將流速剖面儀固定在河床上并采用上視測(cè)量,應(yīng)用平均流速監(jiān)測(cè)模式,監(jiān)測(cè)時(shí)間間隔為10 min,測(cè)量時(shí)段為2019年5月9—21日,監(jiān)測(cè)點(diǎn)位位于S2;9月15—27日監(jiān)測(cè)點(diǎn)位位于S3。同時(shí)段的水位采用高精度水位計(jì)(HOBO U20L- 01,HOBO,加拿大)進(jìn)行觀測(cè),將水位計(jì)固定在河床底部,進(jìn)行水位校正后得到長(zhǎng)期高頻的水位波動(dòng)過程,測(cè)量精度0.5 cm,分辨率0.21 cm,采樣間隔為5 min。各觀測(cè)期間垂向水溫剖面采用多參數(shù)水質(zhì)儀(YSI,美國)監(jiān)測(cè)。流量數(shù)據(jù)包括2019年全年三峽大壩的出庫流量,流量數(shù)據(jù)的時(shí)間分辨率為2 h,即每天12個(gè)數(shù)據(jù);同時(shí)段的壩前水位數(shù)據(jù)時(shí)間分辨率為1 h,即每天24個(gè)數(shù)據(jù)。壩前水位及出庫流量數(shù)據(jù)由三峽集團(tuán)提供。

    表1 各監(jiān)測(cè)點(diǎn)位的地理位置及河道剖面信息

    1.3 數(shù)據(jù)分析方法

    基于高頻連續(xù)的水位、流量時(shí)間序列數(shù)據(jù),本研究采用的主要研究方法包括濾波分析、頻譜分析等。濾波分析主要用于濾除支流庫灣水位變化中的低頻波動(dòng)過程,從而顯示出水位的日波動(dòng)或更高頻的波動(dòng)過程。文中采用Butterworth濾波器對(duì)水位進(jìn)行高通濾波,濾波器的主要參數(shù)包括:濾波器階數(shù)N=20,通過頻率7.7×10-6Hz(約36 h)和6.9×10-5Hz(約4 h),水位樣本頻率3.33×10-3Hz(5 min),流速樣本頻率1.67×10-3Hz(10 min)等;使用標(biāo)準(zhǔn)單位為5 min或10 min。頻譜分析主要用于將時(shí)域下的信號(hào)變換為頻域下的信號(hào),通過對(duì)信號(hào)進(jìn)行傅里葉變換找到信號(hào)在不同頻率下的信息,對(duì)流速數(shù)據(jù)進(jìn)行頻譜分析,得到流速波動(dòng)變化的主要周期。

    2 三峽支流庫灣的水位特性

    2.1 水位和流量調(diào)度過程

    三峽水庫全年一般采取枯期高水位發(fā)電、汛期低水位防洪的運(yùn)營策略。圖2(a)中顯示了2019年大壩出庫流量和茅坪壩前水位變化過程。在1月,水庫水位呈緩慢下降趨勢(shì)并于月底降至170 m,隨后至3月中旬庫區(qū)水位較為穩(wěn)定,在170 m上下小幅波動(dòng);從3月下旬到6月初逐漸騰出庫容防洪,水位持續(xù)下降至147 m 左右;在整個(gè)汛期,水位在150 m上下浮動(dòng),中小洪水調(diào)度過程中水位有一定的抬升;從9月開始,水庫逐漸蓄水,到10月底庫區(qū)水位上升至175 m,隨后保持高水位運(yùn)行。從圖2中還可以看出,除了汛期流量波動(dòng)較大,1—5月和9—12月枯水期日出庫流量并非恒定值,而是一直處于高頻波動(dòng)中。

    圖2(b)和圖2(c)所示分別為汛前泄水期和汛后蓄水期觀測(cè)期內(nèi)(對(duì)應(yīng)圖2(a)中灰色時(shí)段)的出庫流量和庫灣S3斷面日尺度水位波動(dòng)過程。2個(gè)觀測(cè)期內(nèi),出庫流量變化都是典型的雙峰調(diào)節(jié),即早上出庫流量增大,午后出庫流量有小幅降低,傍晚出庫流量再次增大,午夜再次降低。相應(yīng)的支流庫灣水位日波動(dòng)也表現(xiàn)出雙峰波動(dòng),尤其是在出庫流量較大的時(shí)間段,水位日變幅在泄水期和蓄水期分別達(dá)到0.3 m和0.4 m。除了典型的日尺度波動(dòng)外,支流庫灣水位還表現(xiàn)出更高頻的水位波動(dòng),由于日尺度的水位上漲和下落幅度遠(yuǎn)大于高頻的水位波動(dòng),水位上升過程中高頻的波動(dòng)信息并不顯著,而在日尺度波峰和波谷的時(shí)刻,高頻的水位波動(dòng)較為明顯。

    2.2 庫灣沿程的水位波動(dòng)

    為進(jìn)一步分析高頻的水位波動(dòng)過程,對(duì)水位數(shù)據(jù)進(jìn)行高頻濾波分析,結(jié)果如圖3所示。2個(gè)觀測(cè)期內(nèi),這種周期性的水位波動(dòng)都十分顯著。從河口至上游,水位波動(dòng)的振幅沿程增大,且各測(cè)點(diǎn)波動(dòng)基本同步,都沒有明顯的相位差。在S1觀測(cè)點(diǎn),靠近河口的位置,水位波動(dòng)振幅很小(<0.05 m)且波動(dòng)的噪音較大,這可能與該位置受到外界其他干擾有關(guān)(行船擾動(dòng)、風(fēng)浪作用等);而在S3觀測(cè)點(diǎn),水位波動(dòng)振幅可達(dá)到0.1~0.2 m,可見越靠近庫尾的位置,水位波動(dòng)越強(qiáng)。此外,在5月9—15日,振蕩過程表現(xiàn)出一定的日尺度內(nèi)的衰減,日內(nèi)振幅先增加再減??;而在5月15日之后以及9月觀測(cè)期間,所產(chǎn)生的振蕩過程更為持久,能量更強(qiáng),振幅衰減特征不顯著,這可能與該時(shí)段內(nèi)出庫流量日變化相對(duì)較大有關(guān)。

    圖3 不同測(cè)點(diǎn)的水位高頻波動(dòng)Fig.3 High frequency water level fluctuation at different sites

    3 三峽支流庫灣的水流特性

    3.1 分層異重流特征

    三峽水庫作為長(zhǎng)江流域最典型的河道型水庫,其特殊的水動(dòng)力特性一直備受關(guān)注。前期的研究表明[17]:普遍存在的分層異重流(即支流上游入流水體主要從底層流出庫灣,長(zhǎng)江干流水體倒灌潛入支流庫灣)是三峽水庫支流庫灣最特殊、最顯著的水動(dòng)力過程。圖4(a)、圖4(b)分別顯示了泄水期和蓄水期香溪河庫灣監(jiān)測(cè)斷面的水平流速,流速為正表示水流流向庫灣上游。如圖4(a)所示,在泄水觀測(cè)期內(nèi)庫灣水位下降約3 m(圖中紅色實(shí)線為測(cè)點(diǎn)水深),監(jiān)測(cè)斷面流速顯示出顯著的分層異向流特征。距底10 m內(nèi)的水柱流速為負(fù),為上游低溫來流,以順坡異重流的形式流向庫灣下游,最大流速出現(xiàn)在距底2~5 m的位置,流速在0.15 m/s左右。庫灣中上層水體則為下游倒灌潛入的水體,流向相反,最大流速主要集中于距水面10 m 以內(nèi)的位置且表現(xiàn)出一定的上下波動(dòng)特征。在蓄水觀測(cè)期內(nèi)(圖4(b)),庫灣水位上升近8 m,主要集中在前3 d(<1 m),這段時(shí)間干流倒灌水體潛入體量較大,底層上游順坡異重流很弱。后9 d水位上升緩慢,底部順坡異重流僅出現(xiàn)在距底3 m 以內(nèi)且流速較低。除了分層異重流特征,2個(gè)觀測(cè)期內(nèi)剖面流速都還表現(xiàn)出顯著的高頻波動(dòng),如圖中間歇式的紅藍(lán)條帶所示。在三峽水庫中這種高頻的水流波動(dòng)過程長(zhǎng)期以來并未被廣泛報(bào)道。

    圖4 實(shí)測(cè)水平剖面流速Fig.4 Monitored horizontal profile velocity

    3.2 水流振蕩特征

    為進(jìn)一步分析這種顯著的高頻水流波動(dòng)特性,通過高通濾波分析過濾掉低頻的分層異重流成分,得到振蕩流速的變化過程,舉例分析了24 h內(nèi)振蕩流速和水位波動(dòng)過程,如圖5所示(圖中顯示4 h高頻濾波分析)。在泄水期,振蕩流速相對(duì)較小(< 0.05 m/s),水位波動(dòng)振幅也較弱;在蓄水期,振蕩過程更加強(qiáng)烈,振蕩流速達(dá)到0.1 m/s,水位的波動(dòng)幅度可達(dá)0.2 m,水位振幅和流速大小具有很高的一致性。其次,通過對(duì)比實(shí)際流速和振蕩流速的大小,庫灣水流振蕩流速與庫灣的分層異重流處于同一量級(jí),可見振蕩過程對(duì)庫灣水流流速的貢獻(xiàn)不能忽略。觀測(cè)結(jié)果也表明:泄水- 蓄水調(diào)度過程中,流速的這種間歇式增強(qiáng)和減弱在整個(gè)水柱剖面都可見,且具有固定的波動(dòng)周期。在整個(gè)水柱流速保持同步的周期性變化,這一特性表明庫灣振蕩為正壓波[18]。圖5中也顯示了振蕩流速和水位波動(dòng)兩者之間存在π/2個(gè)相位差(T/4,T為庫灣振蕩周期),即水位的波峰和波谷對(duì)應(yīng)于振蕩流速的平衡位置(流速為0的位置),這一特性表明水流振蕩為水平方向傳輸?shù)鸟v波。

    圖5 24 h內(nèi)振蕩流速和水位波動(dòng)過程Fig.5 Water level fluctuation and velocity

    4 三峽支流庫灣振蕩的物理特性

    4.1 振蕩模式

    圖6 庫灣振蕩示意Fig.6 Schematic diagram of tributary oscillation

    根據(jù)香溪河庫灣的水流特性可知,其振蕩的垂向模式為最基本的一層結(jié)構(gòu),具有正壓波的特征。根據(jù)流速和水位波動(dòng)的相位差特性則表明振蕩為水平傳輸?shù)鸟v波,并結(jié)合庫灣水位沿程波動(dòng)的特點(diǎn)(圖3),可知庫灣振蕩為最基本的水平模態(tài)(n),即Helmholtz模式(n=0)[19]。如圖6所示,概括反映了支流庫灣振蕩的模式,圖中箱圖為3個(gè)測(cè)點(diǎn)的實(shí)測(cè)水位波動(dòng)振幅。由于長(zhǎng)江干流水域體量較大,在庫灣振蕩反射過程中,引起的香溪河河口高頻的水位波動(dòng)振幅很小,可忽略不計(jì),該區(qū)域?yàn)檎袷幍牟ü?jié)(水位固定);而在庫灣上游盡頭,振蕩引起的水位波動(dòng)最強(qiáng),為波腹。該模式唯一的波節(jié)線和波腹點(diǎn)分別處于庫灣入口和庫灣盡頭,且該模式的波長(zhǎng)等于庫灣長(zhǎng)度的4倍,以這種方式振蕩的水域被稱為1/4波振蕩器[20]。振蕩過程中,波節(jié)線處受振蕩影響最顯著,流向水平且流速最大[19],庫灣振蕩通過河口周期性的水流交換輸送灣內(nèi)外物質(zhì),對(duì)庫灣物質(zhì)縱向交換將起到重要作用。越靠近上游,在波腹處流速最小且流向垂直[19],庫灣振蕩引起的垂向混合過程更為強(qiáng)烈。

    4.2 振蕩周期

    在振蕩系統(tǒng)中,庫灣所具有的本征頻率(固有頻率)以及模態(tài)結(jié)構(gòu)是其基本屬性[19- 21],主要由水域的幾何形狀和深度決定,并不受外力的影響。三峽水庫支流庫灣在形態(tài)上不同于湖泊水域,相對(duì)于支流庫灣,長(zhǎng)江干流水域體量極大,其形態(tài)特征與海灣、港灣類似,支流庫灣為典型的半封閉系統(tǒng)。在振蕩模式和周期方面的研究中,針對(duì)半封閉的狹長(zhǎng)矩形水域,通過三角函數(shù)構(gòu)建了相應(yīng)的波動(dòng)方程,并推導(dǎo)了在半封閉等深矩形湖庫系統(tǒng)中振蕩周期的參數(shù)方程[22- 23]。根據(jù)香溪河庫灣的實(shí)際幾何特征(縱剖面為三角形)修正其振蕩周期公式,并參照Wilson[23]研究結(jié)果,采用修正系數(shù)α=2.618,得到香溪河庫灣振蕩的理論周期計(jì)算公式。

    (1)

    式中:n為水平模態(tài),取值為0,1,2,…;L為庫灣長(zhǎng)度;g為重力加速度;H為河口處水深。

    由圖4和圖5可知,泄水期和蓄水期庫灣水流振蕩在垂向結(jié)構(gòu)上并無差異,且各自具有固定的振蕩理論和觀測(cè)周期。進(jìn)一步對(duì)2個(gè)觀測(cè)期內(nèi)不同水深的水平流速進(jìn)行頻譜分析,結(jié)果如圖7所示。結(jié)果進(jìn)一步驗(yàn)證了在垂向上振蕩過程的一致性,且具有十分相近的振蕩能量。但泄水- 蓄水過程中的振蕩觀測(cè)周期存在一定的差別,泄水期的振蕩觀測(cè)周期為1.78 h,蓄水期的振蕩觀測(cè)周期為1.94 h,兩者相差約10 min。2次觀測(cè)期間的河口平均水深分別為75 m和80 m,對(duì)應(yīng)的平均庫灣長(zhǎng)度分別為33.6 km和36.5 km,通過式(1)計(jì)算發(fā)現(xiàn):泄水期和蓄水期理論振蕩周期分別為1.80 h和1.89 h,與實(shí)測(cè)值較為吻合;而泄水期和蓄水期振蕩觀測(cè)周期的差異(10 min)是由2個(gè)觀測(cè)期間庫灣的長(zhǎng)度和水深差異造成的。類似地,在三峽水庫全年近30 m的水位波動(dòng)下,由于香溪河庫長(zhǎng)(L=28~42 km)和河口水深(H=63~93 m)等參數(shù)的變化較大,庫灣振蕩的理論周期應(yīng)為T=1.6~2.0 h。

    圖7 不同距底深度處水平流速的能譜分析Fig.7 Energy spectrum analysis of horizontal velocity at different depth from the bottom

    5 支流庫灣振蕩生成機(jī)制及影響因素

    在自然湖泊中,受到風(fēng)或其他外界因素的驅(qū)動(dòng)形成內(nèi)波或湖面波動(dòng)[24],而在本系統(tǒng)中,庫灣振蕩的驅(qū)動(dòng)力來源于水庫的日流量調(diào)節(jié)。其次,香溪河庫灣有別于典型的自然湖泊,可視為半封閉系統(tǒng)(類似于港灣、海港等水體),干支流的交匯作用、相互影響也增加庫灣水動(dòng)力過程的復(fù)雜性。結(jié)合前期研究結(jié)果,本研究總結(jié)了三峽水庫支流庫灣振蕩的形成過程:在三峽水庫的日調(diào)節(jié)調(diào)度過程中,出庫流量的劇烈增加或降低,都會(huì)導(dǎo)致壩前水域水位的降低或升高,這種水位的波動(dòng)會(huì)以重力波的形式向干流上游傳播;在干支流河口,干流的水位擾動(dòng)驅(qū)使干支流的水體交換,并在支流形成駐波,當(dāng)驅(qū)動(dòng)頻率接近支流庫灣的固有頻率時(shí),支流庫灣水位振幅被加強(qiáng),形成高頻的庫灣振蕩。

    對(duì)比不同時(shí)期庫灣振蕩的特征參數(shù),雖然振蕩模式一致,但振蕩周期、振幅等物理參數(shù)存在一定的差異。周期雖然是庫灣系統(tǒng)的固有屬性,但由于三峽水庫不同運(yùn)行期水位差異較大,水位從145 m上升至175 m的過程中庫灣長(zhǎng)度和河口水深都增加了近50 %,庫長(zhǎng)和水深的變化導(dǎo)致了年內(nèi)振蕩周期在1.6~2.0 h之間變化。而振蕩過程中水位振幅大小很大程度上依賴于外界作用力的能量[19],在香溪河庫灣,其驅(qū)動(dòng)能量來源于出庫流量調(diào)節(jié)引起的干流水位波動(dòng),因此,支流庫灣水位波動(dòng)過程可能與出庫流量的變化過程(流量改變大小、變化速率、持續(xù)時(shí)間等)密切相關(guān)。其次,若考慮到振蕩過程的自衰減,河口水深和河道寬度變化也會(huì)影響振蕩的衰減系數(shù)[23];不同水位下河道束窄也會(huì)影響水位波動(dòng)和能量的衰減過程[24]。另一方面,考慮到三峽水庫全年具有多種日調(diào)節(jié)模式,如單峰、雙峰等,紀(jì)道斌等[16]研究也表明了單峰調(diào)節(jié)、雙峰調(diào)節(jié)引起的重力波存在差異,以及多次日調(diào)度過程下庫灣振蕩如何相互疊加,庫灣振蕩如何和分層異重流相互影響,其相互作用機(jī)制有待后續(xù)進(jìn)一步研究。

    前人的研究也報(bào)道過類似由水庫日調(diào)度誘發(fā)的高頻水動(dòng)力過程。謝奇珂等[25]研究發(fā)現(xiàn)水庫日調(diào)節(jié)調(diào)度能誘發(fā)內(nèi)波,造成溫躍層水溫的高頻波動(dòng)。余真真等[26]通過數(shù)值模擬也發(fā)現(xiàn)了同樣來源于水庫日調(diào)度的“假潮”運(yùn)動(dòng),這種類似于長(zhǎng)江河口潮波運(yùn)動(dòng)[27]的往復(fù)水流,可驅(qū)動(dòng)三峽水庫香溪河庫灣形成潮成內(nèi)波,但該研究缺乏實(shí)測(cè)數(shù)據(jù)的驗(yàn)證。本研究中,通過對(duì)泄水和蓄水期間的觀測(cè)發(fā)現(xiàn):支流庫灣振蕩在垂向上始終為正壓波,水平模態(tài)為最基本的模式(Helmholtz模式)。在9月蓄水過程中,庫灣垂向幾乎沒有水溫分層,不具備形成內(nèi)波的條件。在5月,香溪河庫灣已逐漸形成水溫分層,但其振蕩的垂向結(jié)構(gòu)仍為正壓波,在具備密度分層和外界驅(qū)動(dòng)干擾的必要條件下[24]下,并沒有形成內(nèi)波;主要原因可能在于受分層異重流的影響,溫躍層不夠穩(wěn)定。相比于自然湖泊,在頻繁變化的分層異重流、庫灣振蕩的耦合影響下,庫灣振蕩也可能與溫躍層厚度、深度的變化密切相關(guān),在后續(xù)的研究中,將進(jìn)一步論證。

    6 結(jié) 論

    基于三峽水庫香溪河支流庫灣高頻連續(xù)的水位、流場(chǎng)結(jié)構(gòu)等觀測(cè)數(shù)據(jù),系統(tǒng)地分析了三峽水庫泄水期(5月)和蓄水期(9月)支流庫灣的水位、水流波動(dòng)特性,以及庫灣振蕩的物理特征、生成機(jī)制及其主要影響因素,得到了如下結(jié)論:

    (1) 在三峽水庫香溪河庫灣,分層異重流和庫灣振蕩相互疊加作用是該支流庫灣普遍存在的兩大顯著水動(dòng)力過程,其中水流振蕩是長(zhǎng)期被忽視的另一水動(dòng)力過程。

    (2) 庫灣振蕩在垂向上為正壓波,具有典型的駐波特性,水平上表現(xiàn)為最基本的模態(tài),即Helmholtz模式(n=0)。庫灣振蕩引起的水位波動(dòng)在庫尾可達(dá)0.2 m,振蕩流速可達(dá)0.1 m/s,全年的振蕩理論周期約1.6~2.0 h。

    (3) 支流庫灣振蕩形成的外界驅(qū)動(dòng)力為三峽庫區(qū)的干流重力波,其產(chǎn)生于三峽大壩流量日調(diào)節(jié)過程,大壩流量日調(diào)節(jié)的形式及流量調(diào)節(jié)大小是影響支流庫灣振蕩強(qiáng)度的主要因素。

    猜你喜歡
    三峽水庫溪河支流
    世界上最大的升船機(jī)
    安徽裕溪河流域發(fā)現(xiàn)的石制品研究
    東方考古(2020年0期)2020-11-06 05:34:02
    三峽水庫迎戰(zhàn)長(zhǎng)江2020年第2號(hào)洪水
    中國水利(2020年14期)2020-08-02 16:19:30
    美麗的泗溪河
    美麗的泗溪河
    阿什河某支流河岸帶草本植物重要性初步研究
    金沙江支流東川玉碑地遺址
    大眾考古(2014年6期)2014-06-26 08:31:40
    三峽水庫澎溪河水-氣界面CO2、CH4擴(kuò)散通量晝夜動(dòng)態(tài)初探
    三峽水庫綜合管理創(chuàng)新研究
    宛溪河矮塔斜拉橋索塔鞍座區(qū)應(yīng)力分析
    一级片免费观看大全| 亚洲avbb在线观看| 国产在视频线精品| 老司机靠b影院| 国产真人三级小视频在线观看| 久久精品国产亚洲av香蕉五月 | 欧美性长视频在线观看| 久久久久久免费高清国产稀缺| 丝瓜视频免费看黄片| 日韩电影二区| 99精国产麻豆久久婷婷| 99香蕉大伊视频| 欧美国产精品一级二级三级| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美激情久久久久久爽电影 | 日韩熟女老妇一区二区性免费视频| 久久久久精品人妻al黑| 日本精品一区二区三区蜜桃| 丝袜美腿诱惑在线| 国产97色在线日韩免费| 中文字幕av电影在线播放| 色老头精品视频在线观看| 在线亚洲精品国产二区图片欧美| 深夜精品福利| 久久九九热精品免费| √禁漫天堂资源中文www| 一区二区三区乱码不卡18| 亚洲精品中文字幕在线视频| 亚洲成人国产一区在线观看| 久久精品成人免费网站| 夫妻午夜视频| 又紧又爽又黄一区二区| 欧美变态另类bdsm刘玥| 啪啪无遮挡十八禁网站| 2018国产大陆天天弄谢| 女性被躁到高潮视频| 女人被躁到高潮嗷嗷叫费观| 一二三四社区在线视频社区8| 电影成人av| 国产成人一区二区三区免费视频网站| 后天国语完整版免费观看| 91大片在线观看| 成人影院久久| 男男h啪啪无遮挡| 91精品伊人久久大香线蕉| 亚洲va日本ⅴa欧美va伊人久久 | 黄色 视频免费看| 日韩 欧美 亚洲 中文字幕| 巨乳人妻的诱惑在线观看| 午夜免费成人在线视频| 久久人人97超碰香蕉20202| 少妇粗大呻吟视频| 女人高潮潮喷娇喘18禁视频| 国产精品久久久av美女十八| h视频一区二区三区| 人妻 亚洲 视频| 18禁国产床啪视频网站| 国产精品 欧美亚洲| 欧美国产精品一级二级三级| 91九色精品人成在线观看| 一级a爱视频在线免费观看| 国产精品99久久99久久久不卡| 涩涩av久久男人的天堂| 性色av乱码一区二区三区2| 性少妇av在线| 亚洲午夜精品一区,二区,三区| 女人被躁到高潮嗷嗷叫费观| 亚洲色图 男人天堂 中文字幕| 欧美亚洲 丝袜 人妻 在线| 成人av一区二区三区在线看 | 99精品欧美一区二区三区四区| 午夜福利乱码中文字幕| 考比视频在线观看| 精品国产一区二区三区四区第35| 99热网站在线观看| 男人操女人黄网站| 日韩有码中文字幕| 午夜激情久久久久久久| 欧美国产精品一级二级三级| 日本av手机在线免费观看| 国产成人精品久久二区二区免费| 在线观看舔阴道视频| 欧美激情久久久久久爽电影 | 咕卡用的链子| 久久久精品区二区三区| 精品少妇黑人巨大在线播放| 久9热在线精品视频| 搡老乐熟女国产| 欧美乱码精品一区二区三区| 国产黄频视频在线观看| 天堂俺去俺来也www色官网| 亚洲欧美精品综合一区二区三区| 97人妻天天添夜夜摸| 超碰成人久久| 国产亚洲av高清不卡| 少妇 在线观看| 久久久久久久国产电影| 久久久久久久国产电影| 三上悠亚av全集在线观看| 精品亚洲乱码少妇综合久久| 国产日韩欧美亚洲二区| 黄片小视频在线播放| 男女下面插进去视频免费观看| 日本av免费视频播放| 精品欧美一区二区三区在线| 一区二区av电影网| 少妇人妻久久综合中文| 色婷婷av一区二区三区视频| 成人三级做爰电影| 中文字幕av电影在线播放| 岛国在线观看网站| 婷婷色av中文字幕| 亚洲精品中文字幕在线视频| 日韩 亚洲 欧美在线| 热re99久久国产66热| 亚洲av片天天在线观看| 国产亚洲精品一区二区www | 亚洲,欧美精品.| 黑丝袜美女国产一区| 纵有疾风起免费观看全集完整版| 国产精品成人在线| 爱豆传媒免费全集在线观看| 青春草亚洲视频在线观看| 免费一级毛片在线播放高清视频 | 国产激情久久老熟女| 脱女人内裤的视频| 蜜桃国产av成人99| 女人爽到高潮嗷嗷叫在线视频| 成年美女黄网站色视频大全免费| 黄频高清免费视频| 91国产中文字幕| 两人在一起打扑克的视频| 丰满迷人的少妇在线观看| av天堂久久9| 91麻豆av在线| 五月开心婷婷网| 亚洲国产日韩一区二区| 欧美老熟妇乱子伦牲交| 国产真人三级小视频在线观看| 国产日韩欧美在线精品| 欧美在线黄色| 国产精品一区二区在线观看99| 1024视频免费在线观看| 视频区图区小说| 精品第一国产精品| 在线十欧美十亚洲十日本专区| 亚洲精品一区蜜桃| 中文字幕制服av| 国产片内射在线| 丰满少妇做爰视频| 国产精品99久久99久久久不卡| 日韩一区二区三区影片| 精品国产超薄肉色丝袜足j| 男女无遮挡免费网站观看| 日韩,欧美,国产一区二区三区| 欧美日韩一级在线毛片| 人妻人人澡人人爽人人| 一区二区日韩欧美中文字幕| 99精品久久久久人妻精品| 91精品伊人久久大香线蕉| 日本五十路高清| 美女脱内裤让男人舔精品视频| 欧美精品人与动牲交sv欧美| 后天国语完整版免费观看| 免费在线观看完整版高清| 丰满饥渴人妻一区二区三| 成年动漫av网址| 91av网站免费观看| 一进一出抽搐动态| 黄片播放在线免费| 亚洲熟女毛片儿| 久久毛片免费看一区二区三区| 天天添夜夜摸| 国产成人一区二区三区免费视频网站| 亚洲国产精品一区二区三区在线| 黑丝袜美女国产一区| 伊人亚洲综合成人网| 在线观看免费午夜福利视频| 18禁黄网站禁片午夜丰满| 午夜成年电影在线免费观看| 正在播放国产对白刺激| 国产99久久九九免费精品| 久久久久网色| 老熟女久久久| 久久国产精品大桥未久av| 一区在线观看完整版| 亚洲情色 制服丝袜| 97人妻天天添夜夜摸| 黑人巨大精品欧美一区二区蜜桃| 亚洲视频免费观看视频| cao死你这个sao货| 亚洲avbb在线观看| 又紧又爽又黄一区二区| 久久精品aⅴ一区二区三区四区| 90打野战视频偷拍视频| 久久性视频一级片| 两个人看的免费小视频| 一二三四在线观看免费中文在| 青春草视频在线免费观看| av有码第一页| 国产黄色免费在线视频| 免费少妇av软件| 久久久国产精品麻豆| av福利片在线| 一级毛片女人18水好多| 一二三四在线观看免费中文在| 97在线人人人人妻| 久久久欧美国产精品| 狠狠婷婷综合久久久久久88av| tocl精华| 精品一区在线观看国产| 97人妻天天添夜夜摸| 久久天堂一区二区三区四区| 一区二区三区精品91| 91老司机精品| 99国产极品粉嫩在线观看| 黄频高清免费视频| 伦理电影免费视频| 亚洲精品国产区一区二| 成年人免费黄色播放视频| 欧美精品一区二区免费开放| 精品一区二区三区av网在线观看 | 婷婷丁香在线五月| 老熟妇乱子伦视频在线观看 | 大香蕉久久网| 99热国产这里只有精品6| 国产日韩一区二区三区精品不卡| 下体分泌物呈黄色| 国产日韩欧美在线精品| 女人久久www免费人成看片| 麻豆av在线久日| 中文欧美无线码| 人妻一区二区av| 久9热在线精品视频| 久久午夜综合久久蜜桃| 日韩中文字幕欧美一区二区| 亚洲,欧美精品.| 亚洲黑人精品在线| 国产男人的电影天堂91| 欧美日韩av久久| 欧美激情 高清一区二区三区| 一区在线观看完整版| 国产99久久九九免费精品| 自线自在国产av| 国内毛片毛片毛片毛片毛片| 十分钟在线观看高清视频www| 亚洲va日本ⅴa欧美va伊人久久 | 日韩三级视频一区二区三区| 99国产精品一区二区蜜桃av | 老汉色∧v一级毛片| 丝瓜视频免费看黄片| 在线观看免费视频网站a站| 免费高清在线观看日韩| 日韩大码丰满熟妇| 欧美一级毛片孕妇| 高清欧美精品videossex| 妹子高潮喷水视频| 欧美97在线视频| 满18在线观看网站| 久久久精品94久久精品| 一本大道久久a久久精品| 人人妻,人人澡人人爽秒播| 亚洲精品美女久久av网站| 久久精品成人免费网站| 国产亚洲午夜精品一区二区久久| 下体分泌物呈黄色| 欧美激情 高清一区二区三区| 久久免费观看电影| 2018国产大陆天天弄谢| 国产精品1区2区在线观看. | 日韩一卡2卡3卡4卡2021年| 国产人伦9x9x在线观看| 黄网站色视频无遮挡免费观看| 国产三级黄色录像| 十分钟在线观看高清视频www| 看免费av毛片| 精品视频人人做人人爽| 大香蕉久久成人网| 久久久久久久国产电影| 欧美午夜高清在线| 三级毛片av免费| 亚洲精品一卡2卡三卡4卡5卡 | 中国美女看黄片| 2018国产大陆天天弄谢| 欧美精品亚洲一区二区| 亚洲男人天堂网一区| 99精品久久久久人妻精品| 亚洲欧美一区二区三区久久| 男女无遮挡免费网站观看| 国产成人精品久久二区二区91| 亚洲avbb在线观看| 欧美精品亚洲一区二区| 亚洲av美国av| 十分钟在线观看高清视频www| 亚洲va日本ⅴa欧美va伊人久久 | 欧美在线一区亚洲| 亚洲avbb在线观看| 一区二区三区四区激情视频| 中文欧美无线码| av在线播放精品| 亚洲精品国产av蜜桃| 99国产精品免费福利视频| 亚洲美女黄色视频免费看| 黄色怎么调成土黄色| 建设人人有责人人尽责人人享有的| 亚洲国产看品久久| 久久性视频一级片| 女性被躁到高潮视频| 十八禁网站免费在线| 99热国产这里只有精品6| 欧美激情久久久久久爽电影 | 男人爽女人下面视频在线观看| av在线app专区| 久久精品亚洲av国产电影网| av网站在线播放免费| 国产一级毛片在线| 天天操日日干夜夜撸| 国产日韩一区二区三区精品不卡| 伊人久久大香线蕉亚洲五| 男女免费视频国产| 黄色怎么调成土黄色| 亚洲av成人不卡在线观看播放网 | 99国产精品99久久久久| av在线app专区| 久久天躁狠狠躁夜夜2o2o| 亚洲精华国产精华精| 国产欧美日韩精品亚洲av| tocl精华| 一本大道久久a久久精品| 亚洲精品一二三| 99精品久久久久人妻精品| 国产1区2区3区精品| 天天躁夜夜躁狠狠躁躁| 国产三级黄色录像| 亚洲免费av在线视频| 人妻人人澡人人爽人人| 午夜免费观看性视频| 美女主播在线视频| 无限看片的www在线观看| 久久精品国产亚洲av香蕉五月 | 美女国产高潮福利片在线看| 国产伦理片在线播放av一区| 国产精品秋霞免费鲁丝片| 欧美性长视频在线观看| 久热这里只有精品99| 69av精品久久久久久 | 国产精品99久久99久久久不卡| 在线观看一区二区三区激情| a在线观看视频网站| 国产片内射在线| 亚洲精品国产区一区二| av有码第一页| 一级片'在线观看视频| 亚洲专区国产一区二区| 国产一区二区三区综合在线观看| 美女主播在线视频| 国产成+人综合+亚洲专区| 欧美久久黑人一区二区| 日本撒尿小便嘘嘘汇集6| netflix在线观看网站| 涩涩av久久男人的天堂| 老鸭窝网址在线观看| 黄色怎么调成土黄色| 狠狠精品人妻久久久久久综合| 国产免费福利视频在线观看| av网站在线播放免费| 亚洲欧美激情在线| 视频区欧美日本亚洲| 菩萨蛮人人尽说江南好唐韦庄| 男人爽女人下面视频在线观看| 19禁男女啪啪无遮挡网站| 十分钟在线观看高清视频www| 成年人午夜在线观看视频| 热99久久久久精品小说推荐| 日本wwww免费看| 亚洲欧美精品综合一区二区三区| www.999成人在线观看| 亚洲国产av新网站| 久久久久国产精品人妻一区二区| 欧美激情极品国产一区二区三区| 亚洲精品乱久久久久久| 水蜜桃什么品种好| 国产伦理片在线播放av一区| 免费在线观看黄色视频的| 日韩免费高清中文字幕av| 国产又爽黄色视频| 免费看十八禁软件| 国产免费一区二区三区四区乱码| 国产亚洲av片在线观看秒播厂| 中国美女看黄片| 超碰成人久久| 午夜福利在线免费观看网站| 国产精品 欧美亚洲| 纯流量卡能插随身wifi吗| 91精品三级在线观看| 丝袜喷水一区| 国产精品九九99| 欧美另类亚洲清纯唯美| 欧美亚洲日本最大视频资源| 国产精品一区二区在线观看99| 男女国产视频网站| 亚洲全国av大片| 亚洲熟女精品中文字幕| www日本在线高清视频| 日本精品一区二区三区蜜桃| 久久久久久亚洲精品国产蜜桃av| 久久性视频一级片| 黑人欧美特级aaaaaa片| 老熟女久久久| 两性夫妻黄色片| 亚洲va日本ⅴa欧美va伊人久久 | 超色免费av| a级毛片黄视频| 亚洲精品粉嫩美女一区| 欧美另类亚洲清纯唯美| 亚洲精品国产av蜜桃| 母亲3免费完整高清在线观看| 一级毛片精品| 久久毛片免费看一区二区三区| 国产成人欧美在线观看 | 久久久久久久久免费视频了| av欧美777| 日韩大片免费观看网站| 亚洲精品中文字幕在线视频| av超薄肉色丝袜交足视频| 又黄又粗又硬又大视频| 美女高潮喷水抽搐中文字幕| 精品久久久久久电影网| 国产成+人综合+亚洲专区| 亚洲男人天堂网一区| 免费不卡黄色视频| 国产视频一区二区在线看| 丁香六月天网| 在线永久观看黄色视频| 叶爱在线成人免费视频播放| 无遮挡黄片免费观看| 免费观看av网站的网址| 天堂8中文在线网| 精品国产一区二区三区四区第35| 免费高清在线观看视频在线观看| 国产免费一区二区三区四区乱码| 国产av一区二区精品久久| 视频区图区小说| 另类精品久久| 视频在线观看一区二区三区| 日日夜夜操网爽| 另类亚洲欧美激情| 国产97色在线日韩免费| 一本大道久久a久久精品| 久久精品熟女亚洲av麻豆精品| 欧美精品一区二区大全| 欧美成狂野欧美在线观看| 久久久久精品人妻al黑| 可以免费在线观看a视频的电影网站| 国产精品一二三区在线看| 亚洲人成77777在线视频| 亚洲精品国产色婷婷电影| 亚洲少妇的诱惑av| 又大又爽又粗| 精品少妇一区二区三区视频日本电影| 国产一区二区三区在线臀色熟女 | 少妇人妻久久综合中文| 国产亚洲av高清不卡| 性色av乱码一区二区三区2| 精品国内亚洲2022精品成人 | e午夜精品久久久久久久| 蜜桃国产av成人99| 免费不卡黄色视频| 男女无遮挡免费网站观看| 国产精品 欧美亚洲| 亚洲视频免费观看视频| 精品一区二区三区av网在线观看 | 日韩熟女老妇一区二区性免费视频| 亚洲成人手机| 最近最新中文字幕大全免费视频| 九色亚洲精品在线播放| 国产色视频综合| 国产精品影院久久| 欧美另类亚洲清纯唯美| 久久中文看片网| 91老司机精品| 国产在线视频一区二区| 青春草亚洲视频在线观看| 狠狠精品人妻久久久久久综合| 秋霞在线观看毛片| 免费观看av网站的网址| 青青草视频在线视频观看| 国产日韩欧美亚洲二区| 国产精品av久久久久免费| 国产亚洲午夜精品一区二区久久| 国产亚洲av片在线观看秒播厂| 91国产中文字幕| 丝袜脚勾引网站| 69av精品久久久久久 | 国产1区2区3区精品| 国产伦理片在线播放av一区| 久久久久久久精品精品| 色94色欧美一区二区| 无限看片的www在线观看| 交换朋友夫妻互换小说| 在线观看www视频免费| 国产欧美日韩综合在线一区二区| 国产亚洲av片在线观看秒播厂| 日本wwww免费看| 韩国精品一区二区三区| 亚洲av欧美aⅴ国产| 欧美精品人与动牲交sv欧美| 一级毛片精品| 国产免费av片在线观看野外av| 精品一区二区三区四区五区乱码| 啦啦啦啦在线视频资源| 日韩人妻精品一区2区三区| 精品国产国语对白av| 欧美激情 高清一区二区三区| 制服人妻中文乱码| 国产淫语在线视频| 精品国产乱码久久久久久小说| 亚洲欧美激情在线| 一区二区三区精品91| 亚洲欧美日韩另类电影网站| 大型av网站在线播放| 黄片播放在线免费| 国产不卡av网站在线观看| 久久精品久久久久久噜噜老黄| 丝袜在线中文字幕| 一进一出抽搐动态| 一边摸一边做爽爽视频免费| 中文字幕人妻丝袜制服| 久久久国产一区二区| 久久亚洲精品不卡| 高清在线国产一区| 深夜精品福利| 窝窝影院91人妻| 99精品欧美一区二区三区四区| 亚洲精品国产色婷婷电影| 少妇精品久久久久久久| 女警被强在线播放| 亚洲伊人色综图| 亚洲av男天堂| 午夜影院在线不卡| 一级片'在线观看视频| 欧美精品高潮呻吟av久久| 国产97色在线日韩免费| 成年av动漫网址| 91麻豆av在线| 亚洲国产日韩一区二区| 国产成+人综合+亚洲专区| 麻豆av在线久日| 日本猛色少妇xxxxx猛交久久| 一本一本久久a久久精品综合妖精| 精品少妇内射三级| 中国美女看黄片| 精品人妻一区二区三区麻豆| 黑人猛操日本美女一级片| 捣出白浆h1v1| 精品福利观看| 菩萨蛮人人尽说江南好唐韦庄| avwww免费| 中文字幕av电影在线播放| 亚洲av美国av| 国精品久久久久久国模美| 妹子高潮喷水视频| 欧美变态另类bdsm刘玥| 国产精品亚洲av一区麻豆| 丝袜在线中文字幕| 亚洲精品日韩在线中文字幕| 正在播放国产对白刺激| 亚洲午夜精品一区,二区,三区| 叶爱在线成人免费视频播放| 精品卡一卡二卡四卡免费| 国产在线免费精品| 国产av国产精品国产| 啦啦啦 在线观看视频| av国产精品久久久久影院| 国产在线视频一区二区| 免费在线观看日本一区| 一二三四社区在线视频社区8| 成人影院久久| 国产欧美日韩精品亚洲av| 亚洲欧美清纯卡通| 后天国语完整版免费观看| 亚洲欧美色中文字幕在线| 人人澡人人妻人| 永久免费av网站大全| 老司机靠b影院| 自拍欧美九色日韩亚洲蝌蚪91| 91精品三级在线观看| 老司机深夜福利视频在线观看 | 午夜福利影视在线免费观看| 十八禁人妻一区二区| 精品人妻一区二区三区麻豆| 欧美精品av麻豆av| 性色av乱码一区二区三区2| 不卡一级毛片| 亚洲欧洲日产国产| 一个人免费在线观看的高清视频 | 肉色欧美久久久久久久蜜桃| 丰满饥渴人妻一区二区三| 成人国语在线视频| 欧美激情极品国产一区二区三区| 国产精品久久久久久精品电影小说| 午夜精品久久久久久毛片777| 欧美人与性动交α欧美软件| 中文字幕制服av| 夜夜夜夜夜久久久久| a级片在线免费高清观看视频| 高潮久久久久久久久久久不卡| 亚洲欧美一区二区三区黑人| 韩国高清视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 人妻一区二区av| av不卡在线播放| 欧美日韩福利视频一区二区| 久久亚洲精品不卡| 亚洲精品国产av成人精品|