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

    TASCC方法在研究山區(qū)流域氨氮輸移轉(zhuǎn)化中的應(yīng)用

    2023-11-22 03:58:26潘佳佳余甜雪李紀(jì)宇
    水力發(fā)電 2023年11期

    潘佳佳,余甜雪,李紀(jì)宇

    (1.浙江大學(xué)水科學(xué)與工程研究所,浙江 杭州 310058;2.浙江省水利水電勘測設(shè)計院有限責(zé)任公司,浙江 杭州 310002)

    0 引 言

    天然河道特別是源頭溪流自身對氮、磷等營養(yǎng)物質(zhì)具有滯留和吸收功能,通過生物同化、硝化反硝化和沉積物吸附等作用,可顯著降低養(yǎng)分濃度[1]。較大的河床底棲生物表面積促進(jìn)了底棲生物群落對營養(yǎng)物質(zhì)的吸收截留[2],37%~76%的氮進(jìn)入水體會被1~4級小河攔截[3],近年來的養(yǎng)分滯留研究都集中在上游小溪流。為了量化河道的養(yǎng)分吸收,Covino等[4]提出了基于營養(yǎng)鹽瞬時添加示蹤的TASCC(Tracer additions for spiraling curve characterization)方法,優(yōu)化了Newbold等[5]使用的穩(wěn)定同位素示蹤測量法因?qū)嶋H操作困難且分析樣品的設(shè)備昂貴存在的問題。

    1 研究區(qū)域及試驗方法

    1.1 研究區(qū)域

    小港流域位于浙江省麗水市松陽縣西南部山區(qū),流域面積500.3 km2。該區(qū)水系眾多且大多為山溪性河流,河道蜿蜒曲折,上游溪流流速快,河道狹窄而坡降較大;部分下游溪流地形呈深潭或淺灘,流量受降水控制明顯,水位易漲易落,洪水、枯水間差異顯著。該地屬于亞熱帶季風(fēng)氣候,四季分明,溫暖潮濕,年平均氣溫17.7 ℃,1月最冷,7月最熱,年平均降水和蒸發(fā)量分別為1 658.6 mm和1 291 mm。雨季在3月~9月,3月~6月為多雨季節(jié),降水量占全年的53.1%,其中3月~4月屬于春雨期,5月~6月屬于梅雨期;7月~9月為臺風(fēng)雨季,7月、8月氣溫高而降水少,蒸發(fā)大于降水,易出現(xiàn)晴熱少雨的伏旱天氣,降水量占全年的25.6%。11月雨量最小,僅40~50 mm。

    1.2 TASCC營養(yǎng)鹽瞬時添加示蹤試驗

    經(jīng)過前期實地考察,確定了小港流域內(nèi)8處600 m長的試驗河段,基本情況如圖1所示。2022年3月,對小港流域內(nèi)8處試驗河段各開展了一場TASCC氨氮瞬時添加示蹤試驗,每場試驗現(xiàn)場測量,采集溪流特征數(shù)據(jù)。

    圖1 小港流域水系及試驗站點分布

    野外氨氮添加示蹤試驗分為3個步驟:準(zhǔn)備工作、現(xiàn)場試驗和水樣測定。試驗開始前,采集河段的背景濃度水樣,并測定水溫、溪流流速及流量,以此計算NaCl(保守示蹤劑)和(NH4)2SO4(生物活性示蹤劑)示蹤劑投放量。將便攜式電導(dǎo)率儀探頭插入溪流中固定,顯示水中離子濃度實時數(shù)值。在上游起點處將混合示蹤溶液瞬時傾倒至溪流中并開始計時。在600 m的下游采樣處,每分鐘記錄一次電導(dǎo)率,根據(jù)電導(dǎo)率數(shù)據(jù)畫出曲線的上升和下降的斜率決定采樣頻率。曲線上升或下降得越快(即濃度快速變化時),采樣越頻繁,直至電導(dǎo)率曲線回歸平穩(wěn),獲得完整的示蹤劑濃度穿透曲線。采集的水樣需要過濾、冷藏并在24 h內(nèi)運(yùn)往實驗室測定濃度。

    1.3 流域特征數(shù)據(jù)采集及測定

    在野外開展氨氮添加示蹤試驗并計算相關(guān)參數(shù)的基礎(chǔ)上,采集了流域內(nèi)溪流的水文、水質(zhì)和水生態(tài)數(shù)據(jù),包括:①2021年小港水文站的月均流量數(shù)據(jù),直接從小港水文站的水文年鑒中獲??;②小港流域內(nèi)干流下游的試驗點——大東壩溪流口月度水質(zhì)采樣數(shù)據(jù);③小港流域內(nèi)各試驗點溪流的水質(zhì)生態(tài)實地采樣數(shù)據(jù),現(xiàn)場采集并在實驗室測定。

    2 基于TASCC方法的養(yǎng)分吸收參數(shù)計算

    2.1 添加氨氮動態(tài)螺旋指標(biāo)計算

    先利用下式確定添加營養(yǎng)鹽的動態(tài)綜合衰減系數(shù)(kw-add-dyn,m-1),即

    (1)

    添加養(yǎng)分的三項動態(tài)螺旋指標(biāo)Sw-add-dyn、Uadd-dyn和Vf-add-dyn的計算公式為

    (2)

    (3)

    (4)

    (5)

    2.2 氨氮背景濃度螺旋指標(biāo)計算

    (6)

    (7)

    (8)

    2.3 氨氮總動態(tài)螺旋指標(biāo)計算

    在整個添加試驗過程中,背景吸收速率和添加氨氮的動態(tài)吸收速率之和代表了總動態(tài)吸收速率Utot-dyn,公式為

    Utot-dyn=Uamb+Uadd-dyn

    (9)

    總動態(tài)吸收速度Vf-tot-dyn可以計算為

    (10)

    2.4 基于Michaelis-Menten的吸收動力學(xué)模型

    (11)

    3 研究結(jié)果

    3.1 流域水文生態(tài)環(huán)境背景數(shù)據(jù)

    3.1.1 小港水文站月均流量數(shù)據(jù)

    2021年小港水文站逐月月均流量如圖2所示。由圖2可知,小港流域的月均流量呈顯著季節(jié)變化,春末至夏末4月~9月為豐水期,峰值在5月~6月;枯水期從12月~次年3月,流量最低值在1月,與流域雨季分布大體一致。

    圖2 2021年小港水文站逐月月均流量

    3.1.2 大東壩溪流逐月水生態(tài)環(huán)境數(shù)據(jù)

    2021年大東壩溪流逐月水文水質(zhì)數(shù)據(jù)如圖3所示,將大東壩溪流出口處的水文水質(zhì)數(shù)據(jù)按時間順序展示,空白表示實測值空缺。從圖2與圖3a對比可知,大東壩與小港流域的流量變化趨勢一致;圖3b中的水溫與小港流域所處的季節(jié)氣溫變化一致。

    圖3 2021年大東壩逐月水文水質(zhì)數(shù)據(jù)

    由圖3c可知,氨氮濃度的季節(jié)性變化規(guī)律不顯著。由圖3d可知,夏季溪流的pH值略高于冬季,主要原因是冬季CO2溶解度高。由圖3e可知,由于電導(dǎo)率與鹽分濃度呈正相關(guān),而春夏季溪流的電導(dǎo)率高于秋冬季,因此其鹽分濃度也較秋冬高。由圖3f可知,夏季溪流的水的溶解氧(DO)濃度最低,秋冬季節(jié)上升,1月達(dá)峰,春季遞減。這與溫度及水生生物數(shù)量有關(guān),冬季水溫低,DO飽和度高,且生物量少活性差耗氧少,因此溶解氧濃度高。高濃度DO會強(qiáng)化水體的硝化作用而加速消耗氨氮,促進(jìn)消耗大量硝酸鹽氮的藻類繁殖過程,因此,推斷大東壩的氨氮濃度一般應(yīng)是夏季低、冬季高。

    由圖3g可知,化學(xué)需氧量(CODMn)夏季高、冬季低,可能由于夏季陽光充足、氣溫高以及水體生物生長迅速,增加了泵吸作用和底泥中營養(yǎng)鹽的釋放,導(dǎo)致水體需氧量增加。由圖3h可知,生化需氧量(BOD)在春夏高而秋冬低,與水體中生物的生長周期一致。由圖3i可知,總氮(TN)濃度大致呈春夏高,秋冬低,推測由于春夏溫度增高會加強(qiáng)水體微生物的反硝化作用。由圖3j可知,總磷(TP)濃度在5月~9月顯著升高與雨季來臨有關(guān),流域土地內(nèi)積存的磷容易隨雨水沖刷進(jìn)入河流,且夏季微生物生長活躍,聚磷菌的聚磷速度加快。CODMn、BOD、TP的濃度基本都是在春夏高而秋冬低,特別是在豐水期5月、6月都很高。

    由圖3k可知,氟化物在春夏汛期濃度低,冬季枯水期濃度高,枯水期降水少流量小,由于蒸發(fā)和溶解氟化物達(dá)到了最高濃度。由圖3l可知,糞大腸菌群濃度在冬季最低,春夏季節(jié)濃度顯著升高,原因是冬季氣溫較低,不利于水中糞大腸菌群的繁殖,而自春始,水中溫度、光照等條件的上升,為菌群提供了適宜的生長環(huán)境,春夏季菌群數(shù)量有較強(qiáng)的上升趨勢,并在7月達(dá)到最大值。

    3.1.3 各試驗點水生態(tài)環(huán)境采樣數(shù)據(jù)

    圖4 小港流域各試驗點水質(zhì)采樣數(shù)據(jù)

    小港流域各試驗點水文水生態(tài)采樣數(shù)據(jù)如圖5所示。由圖5a、5c可知,從采集試驗點的空間分布來看,在流量和河寬的數(shù)值上,干流試驗點幾乎顯著大于支流。此外,葉綠素a和河道無灰干質(zhì)量(AFDM)作為衡量底棲藻類群落生物量的參數(shù),其中溪流葉綠素a的干流濃度整體都大于支流(見圖5e),底棲葉綠素a支流大于干流,上游大于下游(見圖5f);AFDM在作為干流最下游試驗點的大東壩,數(shù)值明顯高于上游其他位置(見圖5d)。

    圖5 小港流域各試驗點水文水生態(tài)采樣數(shù)據(jù)

    3.2 TASCC方法在小港流域的應(yīng)用分析

    圖6 下宅街試驗點的氯離子與氨氮濃度穿透曲線

    各試驗點穿透曲線的氨氮濃度峰值存在較大區(qū)別,可以較好地反映不同試驗點河道對氨氮吸收能力的差異性。從穿透曲線的形狀看,在小港流域的8場試驗中,大部分試驗都完整呈現(xiàn)了類似的過程,僅有大東壩試驗點的效果不理想,原因是此河段河寬坡緩流速小,河床與氨氮的接觸時間和面積增加。這反映了TASCC法在下游干流的試驗不如上游小溪流適用。此外,根據(jù)氨氮濃度在穿透曲線中到達(dá)峰值的時間分布,推測同時受到流量以及不同河段形態(tài)特征的影響,其中主要是受河段形態(tài)特征的影響,特別是坡度平緩、河道曲度大和存在深潭的情況。下宅街試驗點的穿透曲線中氨氮峰值濃度遠(yuǎn)高于其他試驗點,排居口次之,濃度峰值高是由于不同試驗地點溪流對氨氮的吸收能力有很大差異。

    養(yǎng)分螺旋指標(biāo)包括背景吸收長度Sw-amb、背景吸收速度Vf-amb和背景吸收率Uamb,Sw-amb為營養(yǎng)鹽在溪流水體中行進(jìn)的平均距離,Vf-amb為營養(yǎng)物溶質(zhì)通過沉積物和水交界面沉降遷移的垂直速度,Uamb為單位面積沉積物對水中營養(yǎng)鹽的吸收速率。小港流域八場空間分布試驗的氨氮背景螺旋指標(biāo)見表1。本研究中,Sw-amb整體變化了一個數(shù)量級,而Uamb和Vf-amb變化幅度卻接近三個數(shù)量級。這充分說明Uamb、Vf-amb在不同試驗地點之間存在很大差異,且Sw-amb沒有很好體現(xiàn)出這種差異??梢婋m然都是反映吸收能力大小,不同參數(shù)表示的吸收能力大小的側(cè)重點不同,會反映在氨氮吸收效率的實測結(jié)果中。擬合Michaelis-Menten吸收動力學(xué)模型的參數(shù)Umax和Km,結(jié)果發(fā)現(xiàn),參數(shù)大小會隨不同河道而發(fā)生變化。無論是Umax還是Km都存在很強(qiáng)的空間異質(zhì)性,但與干支流、上下游位置的關(guān)聯(lián)不大。Umax表示生態(tài)系統(tǒng)對營養(yǎng)鹽的滯留潛力,而Km則能夠間接反映生態(tài)系統(tǒng)中生物對營養(yǎng)鹽的親和力,值越低意味著生物對其親和力越高,也就是生物對營養(yǎng)鹽越發(fā)敏感,響應(yīng)速度也就越快。雖然朱岱的吸收速率在較低濃度下響應(yīng)更快,但在大潘坑、下宅街、斗潭的Umax在高濃度下更大。

    表1 小港流域八場空間分布試驗的氨氮背景螺旋指標(biāo)

    當(dāng)把所有試驗點的總氨氮濃度放在同一標(biāo)準(zhǔn)時,不同標(biāo)準(zhǔn)會導(dǎo)致吸收能力的大小及空間分布的排序變化。這表明,若僅使用一次或一個地點計算的吸收參數(shù)的模擬可能不能完全代表水生生態(tài)系統(tǒng)在不同環(huán)境下的養(yǎng)分保留能力[11]。分別用河道環(huán)境中測得的最低(0.05 mg/L)、最高(0.5 mg/L)和平均(0.15 mg/L)氨氮濃度并且額外模擬濃度達(dá)到地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)中的一定限值濃度,如0.15(Ⅰ類)、0.5(Ⅱ類)、1(Ⅲ類)、2(Ⅴ類)mg/L作為基準(zhǔn)氨氮總動態(tài)濃度,代入擬合的Michaelis-Menten吸收動力學(xué)模型,計算得到相對應(yīng)的總吸收速率Utot,用來比較各種濃度水平下時空吸收能力的變化。這組濃度值代表了研究地點測得氨氮背景濃度的范圍,因此稱為基準(zhǔn)濃度[12]。

    各試驗點在不同基準(zhǔn)濃度對應(yīng)的總吸收速率如圖7所示。由圖7可知,不同基準(zhǔn)濃度下的吸收速率在不同河道之間相對大小的關(guān)系有所變化,在基準(zhǔn)氨氮總濃度較低(小于0.25 mg/L)時,普遍呈現(xiàn)出上游的總吸收速率大于下游,干流的總吸收速率大于支流;而當(dāng)基準(zhǔn)濃度增大,甚至超過Ⅴ類水標(biāo)準(zhǔn)時,下游的下宅街河道總吸收速率就超過了上游的大潘坑和朱岱。這意味著在低氨氮濃度輸入的情況下,大潘坑和朱岱河道的氨氮吸收能力反而會比吸收潛力更大的下宅街強(qiáng),只有當(dāng)總動態(tài)濃度達(dá)到一定高度,大潘坑和朱岱才會因為吸收達(dá)到飽和而不及下宅街。這一結(jié)果反映了TASCC法的試驗結(jié)果對從背景濃度到飽和濃度范圍內(nèi)氨氮吸收呈現(xiàn)動態(tài)變化性的過程具有較好的模擬效果。

    圖7 各試驗點在不同基準(zhǔn)濃度對應(yīng)的總吸收速率

    總之,TASCC方法能夠較好地定量刻畫流域內(nèi)不同水文生態(tài)條件的河道從背景濃度至飽和濃度范圍的氨氮動態(tài)吸收過程,展示了養(yǎng)分滯留過程中的動態(tài)變化,為量化氨氮背景濃度下的螺旋指標(biāo)和吸收動力學(xué)模型參數(shù)化以及評估河流吸收的飽和程度提供了一種快速、高效且相對簡單的技術(shù)。

    3.3 各類吸收參數(shù)的適用性分析

    養(yǎng)分螺旋指標(biāo)雖然都代表在養(yǎng)分背景濃度下的吸收效率,但彼此有一定區(qū)別和聯(lián)系。Sw-amb代表吸收能力受到生物和非生物綜合影響的結(jié)果[13],這一指標(biāo)同時受到水文輸運(yùn)和底棲生物交換過程的影響;而Vf-amb更接近于描述非生物吸收的影響,可比較在不同水文條件下河段對養(yǎng)分的吸收能力,甚至可以由河流水力特性估計[14]。氨氮的Uamb、Vf-amb數(shù)值在各試驗點之間的相對大小關(guān)系都基本一致,大致由上游向下游遞減,由大流量向小流量的河道遞減,與Covino等[4]結(jié)果類似,因此認(rèn)為養(yǎng)分螺旋指標(biāo)在小港流域內(nèi)的應(yīng)用效果較好,且與流量存在相關(guān)性。比如,在所有試驗的結(jié)果中Uamb和Vf-amb最大的都是排居口,最小的都是大東壩。排居口的流量最大,其河道較為寬闊,流速很快、水深較深,河床底質(zhì)以卵石和基巖為主;大東壩雖然水面寬度是最大的,但是流量相對較小,流速很緩慢,河岸兩側(cè)碎石塊堆積,植被較少。基于小港流域8處試驗點的的結(jié)果,背景吸收長度Sw-amb數(shù)值越大則越不利于營養(yǎng)鹽在河段中被截留及凈化,除了石倉源、朱岱、排居口和大潘坑外,其余點的Sw-amb并不小于溪流本身長度。說明不同特征的溪流對氨氮的截流凈化能力不盡相同,且變化很大,受到了生物和非生物因素綜合影響,也體現(xiàn)了Sw-amb參數(shù)本身的特征。若將研究專注于Vf-amb和Uamb值上,可以最大程度地減少流量差異造成的影響。因此,建議在比較不同水文條件下不同河道的養(yǎng)分吸收能力差異時,不但要使用綜合性參數(shù)Sw-amb,還非常有必要引入Uamb和Vf-amb進(jìn)行補(bǔ)充對比。

    4 結(jié) 論

    (1)本研究將TASCC養(yǎng)分添加示蹤方法應(yīng)用到小港流域內(nèi)多個試驗河道,此方法具有高效、簡單且成本較低,還能準(zhǔn)確刻畫從氨氮背景濃度到吸收飽和的完整過程的優(yōu)勢,對大部分河道都有很好的適用性,特別是上游小型源頭溪流,但也存在對于下游流速較小的寬河道淺灘試驗耗時過久的缺點。該方法應(yīng)用的條件可能和河道流量、流速、曲度、坡度、斷面水力形態(tài)是否存在深潭的水文特征等因子有關(guān),未來還要進(jìn)一步考慮生態(tài)環(huán)境、基準(zhǔn)濃度的影響。

    (2)養(yǎng)分吸收參數(shù)可以檢驗河道是否存在截留凈化的能力,并比較河道之間吸收能力大小,但不同參數(shù)對于表現(xiàn)吸收能力的側(cè)重點不同。Uamb、Vf-amb主要體現(xiàn)吸收受到的非生物影響,在流域內(nèi)試驗地點之間存在很大差異;Km和Umax則側(cè)重于反映底棲生物吸收能力;而Sw-amb受生物和非生物因素綜合影響,反而沒有體現(xiàn)出這種差別。因此,在比較不同水文水生態(tài)條件下河道的養(yǎng)分吸收能力時,除Sw-amb外也有必要引入其他參數(shù)補(bǔ)充分析。

    (3)流域內(nèi)氨氮吸收能力的空間分布差異可能與河道斷面形態(tài)、流量流速、背景(或基準(zhǔn))濃度和水生態(tài)環(huán)境等因素有關(guān),而小港流域的流量流速、背景濃度和水生態(tài)環(huán)境因素在季節(jié)之間的變化也很大。因此,可以繼續(xù)探究流域氨氮吸收能力的季節(jié)性波動,并考慮溫度、光照、底棲生物等生物相關(guān)因素,還有背景濃度、河道水文特征、降水等非生物因素,以及在季節(jié)變化中這些因素的影響。

    亚洲全国av大片| 搡老妇女老女人老熟妇| 亚洲精品一卡2卡三卡4卡5卡| 一本大道久久a久久精品| 国产一区二区在线av高清观看| 嫩草影视91久久| 韩国精品一区二区三区| 黑人操中国人逼视频| svipshipincom国产片| 两人在一起打扑克的视频| 久久久精品国产亚洲av高清涩受| 成年女人毛片免费观看观看9| 日本三级黄在线观看| 免费观看精品视频网站| 精品久久久久久久人妻蜜臀av| 国产熟女xx| 老司机靠b影院| 一区二区日韩欧美中文字幕| 免费观看人在逋| 国产麻豆成人av免费视频| 久久人妻av系列| 中文字幕精品免费在线观看视频| 99国产精品一区二区三区| 国产一区二区三区视频了| 久久精品夜夜夜夜夜久久蜜豆 | 97碰自拍视频| 欧美激情 高清一区二区三区| 少妇裸体淫交视频免费看高清 | 日本精品一区二区三区蜜桃| 香蕉丝袜av| 男女视频在线观看网站免费 | 在线观看免费日韩欧美大片| 久久久久久亚洲精品国产蜜桃av| 亚洲欧洲精品一区二区精品久久久| 午夜日韩欧美国产| 国产黄片美女视频| 精品国产亚洲在线| 午夜视频精品福利| 欧美成人午夜精品| 久久香蕉精品热| 午夜福利高清视频| 国产精品久久久av美女十八| 久热这里只有精品99| 美女高潮喷水抽搐中文字幕| 亚洲av成人av| av在线天堂中文字幕| 久久久久久久精品吃奶| 国产亚洲av嫩草精品影院| 亚洲成人久久爱视频| 欧美成人免费av一区二区三区| 两个人看的免费小视频| 中文字幕精品免费在线观看视频| 女同久久另类99精品国产91| 1024香蕉在线观看| 母亲3免费完整高清在线观看| 91麻豆av在线| 精品第一国产精品| 一级黄色大片毛片| 欧美大码av| 午夜久久久久精精品| 中文字幕av电影在线播放| 香蕉久久夜色| 又紧又爽又黄一区二区| 久久久久亚洲av毛片大全| 久久天堂一区二区三区四区| or卡值多少钱| 真人一进一出gif抽搐免费| 国产97色在线日韩免费| 日本一本二区三区精品| 欧美国产精品va在线观看不卡| 久久久久久免费高清国产稀缺| 美女免费视频网站| av有码第一页| 天天添夜夜摸| 操出白浆在线播放| 视频在线观看一区二区三区| 少妇的丰满在线观看| 国产激情欧美一区二区| 亚洲 欧美 日韩 在线 免费| 亚洲精品久久成人aⅴ小说| 91成年电影在线观看| 18美女黄网站色大片免费观看| 免费在线观看黄色视频的| 成年免费大片在线观看| 丰满人妻熟妇乱又伦精品不卡| 男男h啪啪无遮挡| 少妇熟女aⅴ在线视频| 一级黄色大片毛片| 久久精品人妻少妇| 一进一出好大好爽视频| 天天添夜夜摸| 黄片大片在线免费观看| 欧美丝袜亚洲另类 | 亚洲国产欧美网| 少妇 在线观看| 人人妻,人人澡人人爽秒播| 国产蜜桃级精品一区二区三区| av有码第一页| 在线观看午夜福利视频| 最新在线观看一区二区三区| 啦啦啦观看免费观看视频高清| 欧美激情久久久久久爽电影| 亚洲专区中文字幕在线| 午夜福利高清视频| 午夜福利视频1000在线观看| svipshipincom国产片| 欧美日韩瑟瑟在线播放| 在线观看www视频免费| av中文乱码字幕在线| 午夜福利免费观看在线| 这个男人来自地球电影免费观看| 99热6这里只有精品| 精品欧美一区二区三区在线| 最近最新中文字幕大全电影3 | 国产成+人综合+亚洲专区| 国产99白浆流出| 美女免费视频网站| 成人免费观看视频高清| 在线观看午夜福利视频| 十分钟在线观看高清视频www| 久久人妻av系列| 美女国产高潮福利片在线看| 亚洲av成人不卡在线观看播放网| 久久精品国产清高在天天线| 日本免费a在线| 别揉我奶头~嗯~啊~动态视频| 日本三级黄在线观看| 国产免费男女视频| 成人av一区二区三区在线看| 男人舔女人下体高潮全视频| 在线观看免费日韩欧美大片| 国产激情欧美一区二区| 精品久久久久久,| 久久久国产精品麻豆| 午夜老司机福利片| 午夜激情av网站| 91麻豆精品激情在线观看国产| 一本一本综合久久| 欧美精品啪啪一区二区三区| 韩国精品一区二区三区| 18禁黄网站禁片免费观看直播| avwww免费| 十八禁网站免费在线| 日本 欧美在线| 在线av久久热| 中文字幕最新亚洲高清| 最新美女视频免费是黄的| 可以在线观看毛片的网站| 国产精品乱码一区二三区的特点| 12—13女人毛片做爰片一| 国产真人三级小视频在线观看| 亚洲狠狠婷婷综合久久图片| tocl精华| www.熟女人妻精品国产| 免费在线观看视频国产中文字幕亚洲| 少妇粗大呻吟视频| 国产精品,欧美在线| av视频在线观看入口| 亚洲成a人片在线一区二区| 日韩大尺度精品在线看网址| 一级a爱视频在线免费观看| 97人妻精品一区二区三区麻豆 | 最近在线观看免费完整版| 国产精品电影一区二区三区| 男人舔女人下体高潮全视频| 亚洲成av人片免费观看| 青草久久国产| 此物有八面人人有两片| 日本五十路高清| 国产精品香港三级国产av潘金莲| 黄频高清免费视频| 很黄的视频免费| 国产成人av激情在线播放| 一级a爱视频在线免费观看| 亚洲国产精品sss在线观看| 国产一区二区在线av高清观看| 禁无遮挡网站| 级片在线观看| 午夜激情av网站| 国产不卡一卡二| 亚洲av熟女| www国产在线视频色| 美国免费a级毛片| 91字幕亚洲| 亚洲精品色激情综合| 好男人在线观看高清免费视频 | 亚洲第一av免费看| 国产av又大| 久久久久免费精品人妻一区二区 | 97超级碰碰碰精品色视频在线观看| 成人国产综合亚洲| 久久久水蜜桃国产精品网| 97超级碰碰碰精品色视频在线观看| 精品福利观看| 天堂√8在线中文| 国产亚洲欧美98| 国产精品一区二区精品视频观看| 久久久国产精品麻豆| 欧美人与性动交α欧美精品济南到| 久久久久国产一级毛片高清牌| 国产激情偷乱视频一区二区| 国产黄片美女视频| 啦啦啦免费观看视频1| 国产黄色小视频在线观看| 国产97色在线日韩免费| 午夜久久久在线观看| 99精品在免费线老司机午夜| 女同久久另类99精品国产91| 草草在线视频免费看| 成在线人永久免费视频| 亚洲人成电影免费在线| 人妻丰满熟妇av一区二区三区| 国产午夜福利久久久久久| 国产精品,欧美在线| 桃红色精品国产亚洲av| 亚洲天堂国产精品一区在线| 国产精品亚洲一级av第二区| 久久精品aⅴ一区二区三区四区| 免费搜索国产男女视频| 国产久久久一区二区三区| 少妇 在线观看| 宅男免费午夜| 青草久久国产| 成人欧美大片| 神马国产精品三级电影在线观看 | 国产av一区二区精品久久| 搡老妇女老女人老熟妇| 99久久精品国产亚洲精品| 一级毛片高清免费大全| 国产99白浆流出| 99热只有精品国产| 大香蕉久久成人网| 久久国产精品人妻蜜桃| 女同久久另类99精品国产91| 宅男免费午夜| 成年人黄色毛片网站| 亚洲欧美日韩高清在线视频| 亚洲精品中文字幕在线视频| 国产av不卡久久| 美女 人体艺术 gogo| 一本综合久久免费| 亚洲国产看品久久| a级毛片在线看网站| 午夜久久久久精精品| 欧美日韩黄片免| 免费高清视频大片| 欧美日韩一级在线毛片| 日日爽夜夜爽网站| 黄片小视频在线播放| 性欧美人与动物交配| 老汉色∧v一级毛片| 大型av网站在线播放| 成人精品一区二区免费| 亚洲av成人av| 久久国产精品人妻蜜桃| 在线永久观看黄色视频| 精品欧美一区二区三区在线| 最好的美女福利视频网| 久久久精品国产亚洲av高清涩受| 12—13女人毛片做爰片一| 久久久久免费精品人妻一区二区 | 日韩有码中文字幕| 国语自产精品视频在线第100页| 露出奶头的视频| 制服人妻中文乱码| 欧美av亚洲av综合av国产av| 精品无人区乱码1区二区| 在线观看www视频免费| 成人午夜高清在线视频 | 亚洲精品中文字幕一二三四区| 脱女人内裤的视频| 天堂动漫精品| 最新美女视频免费是黄的| 久久久久亚洲av毛片大全| 久久性视频一级片| 精品不卡国产一区二区三区| 日日干狠狠操夜夜爽| 精品高清国产在线一区| 日韩三级视频一区二区三区| 夜夜躁狠狠躁天天躁| avwww免费| 两人在一起打扑克的视频| 国产一区二区三区视频了| 亚洲黑人精品在线| 不卡av一区二区三区| 亚洲男人天堂网一区| 久久香蕉激情| 国产在线精品亚洲第一网站| 亚洲精品中文字幕在线视频| 久久久久久久午夜电影| 婷婷丁香在线五月| 午夜激情av网站| 美女免费视频网站| 中文字幕另类日韩欧美亚洲嫩草| 少妇的丰满在线观看| 久久久久国产精品人妻aⅴ院| 日本 av在线| 国产亚洲av嫩草精品影院| 国产99久久九九免费精品| 精品人妻1区二区| cao死你这个sao货| 国产又黄又爽又无遮挡在线| 悠悠久久av| 国产精品av久久久久免费| 精品日产1卡2卡| 国产亚洲av高清不卡| 亚洲无线在线观看| 亚洲欧美精品综合久久99| 一边摸一边做爽爽视频免费| 国产一区二区激情短视频| 成年免费大片在线观看| 国产精品一区二区精品视频观看| 亚洲,欧美精品.| 国产精品美女特级片免费视频播放器 | 久久久久国产精品人妻aⅴ院| 欧美人与性动交α欧美精品济南到| 99久久99久久久精品蜜桃| 国产av在哪里看| cao死你这个sao货| 啦啦啦观看免费观看视频高清| 亚洲精品一区av在线观看| 色精品久久人妻99蜜桃| 啦啦啦观看免费观看视频高清| 色在线成人网| 久久天躁狠狠躁夜夜2o2o| 怎么达到女性高潮| 成人三级黄色视频| 在线观看午夜福利视频| 色播在线永久视频| 久久亚洲精品不卡| 国产熟女午夜一区二区三区| 国产伦在线观看视频一区| 欧美午夜高清在线| 久久久精品国产亚洲av高清涩受| 老司机在亚洲福利影院| 国产三级黄色录像| svipshipincom国产片| 欧美又色又爽又黄视频| 精品无人区乱码1区二区| 久久久久久九九精品二区国产 | 日韩欧美国产在线观看| 国产精品,欧美在线| 亚洲第一电影网av| 欧美午夜高清在线| 精品第一国产精品| 日本精品一区二区三区蜜桃| xxx96com| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩高清在线视频| 日韩三级视频一区二区三区| 一级毛片高清免费大全| 日韩视频一区二区在线观看| 中文资源天堂在线| 午夜亚洲福利在线播放| 久久久久国产精品人妻aⅴ院| 免费女性裸体啪啪无遮挡网站| 成人一区二区视频在线观看| 一级a爱片免费观看的视频| 99riav亚洲国产免费| 一区二区三区国产精品乱码| 国产成人系列免费观看| 亚洲成av人片免费观看| 欧美日韩瑟瑟在线播放| 宅男免费午夜| 国产精品美女特级片免费视频播放器 | 女人爽到高潮嗷嗷叫在线视频| 18禁黄网站禁片免费观看直播| 亚洲狠狠婷婷综合久久图片| 色av中文字幕| 老司机靠b影院| 亚洲av熟女| 久久国产亚洲av麻豆专区| 午夜a级毛片| 欧美日韩福利视频一区二区| 手机成人av网站| 91成年电影在线观看| 久久午夜综合久久蜜桃| 亚洲成人精品中文字幕电影| tocl精华| 亚洲午夜精品一区,二区,三区| av有码第一页| 级片在线观看| 美女免费视频网站| 午夜老司机福利片| 中文亚洲av片在线观看爽| 天天添夜夜摸| 人人澡人人妻人| 亚洲国产毛片av蜜桃av| 两性夫妻黄色片| 无限看片的www在线观看| 欧美黑人欧美精品刺激| 欧美丝袜亚洲另类 | 免费电影在线观看免费观看| 欧美色欧美亚洲另类二区| 欧美成人午夜精品| 亚洲九九香蕉| 精品国内亚洲2022精品成人| 黄色女人牲交| 国内精品久久久久久久电影| 亚洲国产高清在线一区二区三 | 成熟少妇高潮喷水视频| 日本免费一区二区三区高清不卡| 欧美黑人巨大hd| 91九色精品人成在线观看| 亚洲精品美女久久久久99蜜臀| 精品午夜福利视频在线观看一区| 级片在线观看| 久久国产精品男人的天堂亚洲| 三级毛片av免费| 神马国产精品三级电影在线观看 | 麻豆一二三区av精品| 美女高潮喷水抽搐中文字幕| 欧美成人性av电影在线观看| 18美女黄网站色大片免费观看| 日韩大尺度精品在线看网址| 国产又色又爽无遮挡免费看| 亚洲色图 男人天堂 中文字幕| 一二三四在线观看免费中文在| 欧美成人性av电影在线观看| 久久久久免费精品人妻一区二区 | 欧美激情 高清一区二区三区| 很黄的视频免费| 侵犯人妻中文字幕一二三四区| 亚洲国产欧洲综合997久久, | 日韩免费av在线播放| 这个男人来自地球电影免费观看| 岛国在线观看网站| 精品久久久久久久毛片微露脸| 亚洲av中文字字幕乱码综合 | 成人一区二区视频在线观看| 91成人精品电影| 国产精品精品国产色婷婷| 精品国内亚洲2022精品成人| 美女高潮喷水抽搐中文字幕| 国产av一区在线观看免费| 97人妻精品一区二区三区麻豆 | 欧美精品啪啪一区二区三区| 黑丝袜美女国产一区| 一本精品99久久精品77| 国产三级在线视频| 久久这里只有精品19| 亚洲 欧美 日韩 在线 免费| cao死你这个sao货| 99国产精品99久久久久| 亚洲熟妇熟女久久| 高清在线国产一区| 亚洲国产精品999在线| 三级毛片av免费| 国产一区二区三区视频了| 成人免费观看视频高清| 伦理电影免费视频| 一a级毛片在线观看| 国产日本99.免费观看| 91在线观看av| 亚洲成av片中文字幕在线观看| 在线av久久热| 久久婷婷成人综合色麻豆| 色精品久久人妻99蜜桃| 天天一区二区日本电影三级| 午夜福利免费观看在线| 国内精品久久久久精免费| 视频区欧美日本亚洲| 中文字幕精品免费在线观看视频| 中亚洲国语对白在线视频| 麻豆成人av在线观看| 俺也久久电影网| 国产精品综合久久久久久久免费| 首页视频小说图片口味搜索| 国产精品久久久人人做人人爽| 亚洲在线自拍视频| 听说在线观看完整版免费高清| 欧美日韩瑟瑟在线播放| 成人手机av| 亚洲国产高清在线一区二区三 | 国产高清视频在线播放一区| 亚洲中文av在线| 91老司机精品| 男女做爰动态图高潮gif福利片| bbb黄色大片| 欧美成人午夜精品| 亚洲一码二码三码区别大吗| 精品国产国语对白av| 19禁男女啪啪无遮挡网站| 精品第一国产精品| 18禁国产床啪视频网站| 国产精品永久免费网站| 亚洲国产精品sss在线观看| 午夜激情福利司机影院| 国产野战对白在线观看| svipshipincom国产片| 变态另类丝袜制服| 脱女人内裤的视频| 日本精品一区二区三区蜜桃| 制服人妻中文乱码| av在线播放免费不卡| 国产成人一区二区三区免费视频网站| 亚洲熟女毛片儿| 69av精品久久久久久| 美国免费a级毛片| 黄色视频不卡| 九色国产91popny在线| 婷婷六月久久综合丁香| www国产在线视频色| 国产男靠女视频免费网站| 久久久久久久午夜电影| 手机成人av网站| 真人一进一出gif抽搐免费| 免费看a级黄色片| 国产精品永久免费网站| 视频区欧美日本亚洲| 又黄又爽又免费观看的视频| 欧美乱妇无乱码| 2021天堂中文幕一二区在线观 | 精品国产美女av久久久久小说| 中文字幕av电影在线播放| 亚洲精品中文字幕一二三四区| 在线观看免费视频日本深夜| 99热6这里只有精品| 国产不卡一卡二| 免费观看精品视频网站| 狠狠狠狠99中文字幕| 精品国产亚洲在线| 国产精品亚洲av一区麻豆| 免费电影在线观看免费观看| 欧美日本视频| 99riav亚洲国产免费| 自线自在国产av| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区在线av高清观看| 岛国在线观看网站| 中亚洲国语对白在线视频| 久久香蕉精品热| 大香蕉久久成人网| 久久国产精品影院| 亚洲av熟女| 国产熟女xx| 黑丝袜美女国产一区| 日日夜夜操网爽| 十八禁网站免费在线| 观看免费一级毛片| 免费在线观看完整版高清| 黄片播放在线免费| 国产精品一区二区三区四区久久 | 午夜两性在线视频| 热99re8久久精品国产| 国产高清激情床上av| 成人精品一区二区免费| 1024视频免费在线观看| 亚洲av成人不卡在线观看播放网| 欧美色欧美亚洲另类二区| 久久人人精品亚洲av| 少妇的丰满在线观看| 日韩欧美国产一区二区入口| 免费在线观看亚洲国产| 国产高清视频在线播放一区| 亚洲成人久久爱视频| 国内少妇人妻偷人精品xxx网站 | 亚洲国产精品sss在线观看| 成人三级做爰电影| 久久精品aⅴ一区二区三区四区| 中文字幕最新亚洲高清| 十八禁网站免费在线| 黑人操中国人逼视频| 高清毛片免费观看视频网站| 亚洲成人精品中文字幕电影| 国内精品久久久久精免费| av视频在线观看入口| 国产爱豆传媒在线观看 | 手机成人av网站| 18禁黄网站禁片午夜丰满| 久久精品国产亚洲av香蕉五月| 在线国产一区二区在线| 一本久久中文字幕| 国产成人欧美| aaaaa片日本免费| 亚洲国产精品sss在线观看| www国产在线视频色| 男人舔女人下体高潮全视频| 日本在线视频免费播放| 一区二区日韩欧美中文字幕| 中文字幕最新亚洲高清| 欧美黑人欧美精品刺激| 真人做人爱边吃奶动态| 91av网站免费观看| 国产精品99久久99久久久不卡| 99riav亚洲国产免费| 国产精品综合久久久久久久免费| 无人区码免费观看不卡| 波多野结衣av一区二区av| 国产又黄又爽又无遮挡在线| 亚洲欧美一区二区三区黑人| 18禁观看日本| 制服人妻中文乱码| 亚洲av第一区精品v没综合| 999久久久国产精品视频| 在线av久久热| 久久精品国产亚洲av香蕉五月| 久久热在线av| 变态另类丝袜制服| 亚洲av成人一区二区三| 国产精品av久久久久免费| 欧美精品啪啪一区二区三区| 久久精品成人免费网站| 身体一侧抽搐| 老司机午夜十八禁免费视频| 午夜福利欧美成人| 91老司机精品| 黄色视频,在线免费观看| 欧美精品啪啪一区二区三区| 高清毛片免费观看视频网站| 国产精品综合久久久久久久免费| 在线av久久热| 18禁国产床啪视频网站| 欧美精品啪啪一区二区三区| 久久久精品国产亚洲av高清涩受|