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

    氣候變化和區(qū)外調(diào)水對米山水庫水質(zhì)的影響

    2022-09-05 08:08:16張明芳黃林顯劉建光宋協(xié)良李文虎
    濟南大學學報(自然科學版) 2022年5期
    關(guān)鍵詞:水質(zhì)模型

    張明芳,黃林顯,劉 揚,張 杰,劉建光,宋協(xié)良,李文虎

    (1. 威海市水文中心,山東 威海 264209;2. 濟南大學 a. 水利與環(huán)境學院,b. 山東省地下水數(shù)值模擬與污染控制工程技術(shù)研究中心,山東 濟南 250022;3. 山東省地質(zhì)測繪院,山東 濟南 250013;4. 威海市文登區(qū)水利事務(wù)服務(wù)中心,山東 威海 264400)

    米山水庫是山東省威海市境內(nèi)最大的水庫,也是威海市最重要的城市供水水源。米山水庫作為南水北調(diào)工程山東段末端的調(diào)蓄水庫,主要承擔著威海市主城區(qū)和文登區(qū)的供水任務(wù),其水質(zhì)對城市供水安全和區(qū)域水生態(tài)穩(wěn)定至關(guān)重要[1]。湖庫水質(zhì)主要受氣候變化和人類活動的影響,其動態(tài)變化表現(xiàn)出復雜的趨勢性、季節(jié)性及隨機性等特征。對于米山水庫來說,一方面隨著全球性氣候的變化,近年來膠東半島連續(xù)多年降雨量減少,氣候變暖,導致米山水庫入庫水量減少,對米山水庫的水質(zhì)產(chǎn)生影響;另一方面,米山水庫自2015年12月21日開始區(qū)外調(diào)水,調(diào)水水源包括黃河水、長江水和當?shù)厮庹{(diào)水同樣對米山水庫的水質(zhì)產(chǎn)生了顯著影響。文獻[2]中的研究已經(jīng)指出,米山水庫2018年水體中的電導率、鉀、鈉、鈣、氯化物、硫酸鹽、重碳酸鹽等指標與2015年均值都有較大幅度的增長[2]。通過定量評估氣候變化和區(qū)外調(diào)水對米山水庫水質(zhì)演變的貢獻程度并評價在不同時期對水質(zhì)變化的作用過程,可以為米山水庫建立更加全面的水質(zhì)安全風險管理體系,實施更加科學的水量-水質(zhì)聯(lián)合調(diào)控決策以及維護供水安全提供重要借鑒和參考。

    前人關(guān)于調(diào)水工程對湖庫水質(zhì)影響的研究大多集中在水質(zhì)評價[3]、風險評估[4-5]和調(diào)控政策[6]等方面,從定量角度評價氣候變化和區(qū)外調(diào)水對湖庫水質(zhì)演變貢獻度的研究尚未見報道。氣候變化和區(qū)外調(diào)水貢獻度的定量評價需要科學選擇水體水質(zhì)的度量指標。在水體的眾多指標中,水的電導率主要用來表征水體傳導電流能力的強弱,通過電導率值的大小可以推測水中帶電荷物質(zhì)的總濃度,因此可以作為反映水體污染情況的一個重要度量[7]。掌握電導率在水體中的時空分布對進行湖庫生態(tài)環(huán)境變化趨勢分析具有重要的參考價值。

    本文中首先利用局部加權(quán)回歸的時間序列分解(STL)法[8-10]對降雨及電導率時間序列進行時序分解,揭示它們的變化趨勢、周期長度、隨機波動范圍等,分別提取出電導率的自然狀態(tài)(受氣候變化影響)和混合狀態(tài)(同時受氣候變化和區(qū)外調(diào)水影響)時間序列。利用自然狀態(tài)時間序列通過灰色系統(tǒng)GM(1,1)分解模型建立電導率模擬模型,得到自然狀態(tài)下的電導率模擬值,通過與混合狀態(tài)下實測電導率的比較,定量評價氣候變化和區(qū)外調(diào)水對米山水庫水質(zhì)演變的貢獻程度。

    1 數(shù)據(jù)序列和研究方法

    米山水庫位于山東省威海市文登區(qū)米山鎮(zhèn),是以防洪、城鎮(zhèn)供水、發(fā)電等綜合利用的大(二)型水庫,控制流域面積為440 km2,總庫容為2.8×108m3,死庫容為5.07×106m3。2015年12月22日南水北調(diào)東線黃河水入威開始通水,2016年3月10日,長江水首次調(diào)入米山水庫,威海市實現(xiàn)長江水、黃河水和當?shù)厮嗨绰?lián)合調(diào)度。米山水庫流域及雨量站分布如圖1所示。

    圖1 米山水庫流域及雨量站分布

    1.1 數(shù)據(jù)序列

    1.1.1 數(shù)據(jù)來源

    受制于月數(shù)據(jù)收集的困難,本研究中STL時間序列分解采用1980—2019年月均降雨量數(shù)據(jù)和2012—2019年月均電導率數(shù)據(jù)(見圖2),通過對時間序列進行統(tǒng)計分析來提取歷史數(shù)據(jù)信息,以此找出其動態(tài)變化規(guī)律,揭示行為的未來演變趨勢;其中2013年電導率數(shù)據(jù)部分缺失,采用樣條插值進行插補(見圖2)。GM(1,1)分解模型采用1987—2019年年均電導率數(shù)據(jù),通過建立電導率模擬模型,對電導率進行模擬計算。降雨數(shù)據(jù)來自于中國氣象科學數(shù)據(jù)共享服務(wù)網(wǎng),電導率數(shù)據(jù)來自于威海市水文中心監(jiān)測數(shù)據(jù)。

    1.1.2 降雨及電導率統(tǒng)計特征

    對米山水庫1980—2019年年均降雨量和1987—2019年年均電導率長時間序列進行統(tǒng)計分析,結(jié)果如表1所示。從表中可以看出,研究區(qū)多年年均降雨量為811.67 mm,降雨較為豐富;年均降雨量的標準差為205.66 mm,說明研究區(qū)降雨的波動性相對較為平穩(wěn);年均電導率為374.98 μS/cm,說明水體中離子含量較多;年均電導率的標準差為215.02 μS/cm,說明電導率波動性較大。

    (a)1980—2019年月均降雨量

    (b)2012—2019年月均電導率圖2 米山水庫月均降雨量和月均電導率及樣條插值數(shù)據(jù)

    表1 米山水庫年均降雨量和電導率時間序列數(shù)據(jù)統(tǒng)計特征

    1.2 研究方法

    1.2.1 STL法

    STL法是一種非參數(shù)統(tǒng)計方法,它以具有魯棒性的局部加權(quán)回歸散點平滑(LOESS)作為平滑方法,能夠?qū)r間序列以加法模型為原理分解為趨勢變化、季節(jié)性變化及不規(guī)則的隨機波動項。主要計算步驟參見文獻[8-9],計算公式為

    Yt=Tt+St+Rt,

    (1)

    式中:Yt是t時刻的觀測值;Tt是t時刻的趨勢變化,反映時間序列的變化趨勢,如增加或者減少;St是t時刻的季節(jié)性變化,反映時間序列的波動周期,如12個月;Rt是t時刻的隨機波動項,反映受外部因素影響的強弱。

    1.2.2 GM(1,1)分解模型

    目前常用的水質(zhì)模擬預測方法主要有數(shù)理統(tǒng)計法[11-12]、機器學習模型法[13-16]、數(shù)值模型法[17-18]和灰色模型法[19-22]等,其中灰色模型法中的GM(1,1)模型因用法簡單且所需實測數(shù)據(jù)少而被廣泛應用,但是存在灰色偏差和抗干擾能力弱的局限性,只適合于對呈指數(shù)規(guī)律變化且增長速率較低的短時間序列進行預測[23]。水質(zhì)、水量等水文要素的增長速率有快有慢,并且數(shù)據(jù)本身存在周期變化和隨機擾動,造成傳統(tǒng)GM(1,1)模型對于水文要素的預測精度不高。為了解決此類問題,Hao等[24]提出了GM(1,1)分解模型,并將其應用于水文時間序列的模擬預測,取得了非常理想的預測效果。GM(1,1)分解模型將時間序列分解為趨勢項、周期波動項和隨機變動項,因此能較好地擬合水文要素的動態(tài)變化特征。李華敏等[25]將GM(1,1)分解模型應用于山西柳林泉水流量的預測中,并把預測結(jié)果與自回歸差分移動平均 (autoregressive integrated moving average, ARIMA)模型進行對比,結(jié)果顯示通過周期修正能顯著提高GM(1,1)分解模型預測精度,其預測結(jié)果要好于ARIMA模型的。

    GM(1,1)分解模型的工作原理分述如下。

    1)時間序列的趨勢項。利用GM(1,1)模型模擬計算時間序列的趨勢項。首先對原始時間序列x(0)(t)采用累加的方式生成一組趨勢明顯的新時間序列x(1)(t),然后建立微分方程

    (2)

    式中:a為系統(tǒng)的發(fā)展系數(shù);u為灰色作用量。

    (3)

    然后,進行累減操作恢復原始數(shù)據(jù)序列,得到預測結(jié)果如下:

    (4)

    2)時間序列的周期波動項。周期波動項通過利用正弦函數(shù)對殘差序列進行擬合獲得。殘差序列ε(t)計算公式如下:

    (5)

    文獻[24]中對殘差序列采用分段正弦曲線擬合的方法,本文中為了避免人為選擇分段周期和振幅所帶來的誤差,采用帶有截距的正弦曲線進行擬合,

    (6)

    (7)

    3)時間序列的隨機波動項。時間序列的隨機波動項R(t)通過以下公式計算獲得:

    (8)

    合并趨勢項、周期波動項和隨機波動項,最終得到對時間序列動態(tài)變化的模擬預測值Q(t),

    (9)

    2 實驗結(jié)果與分析

    2.1 研究方法

    利用STL法對1980—2019年月降雨量及2012—2019年月電導率值進行時序分解,揭示它們的變化趨勢、周期長度、隨機波動范圍等,并進一步提取出電導率的自然狀態(tài)(受氣候變化影響)和混合狀態(tài)(同時受氣候變化和區(qū)外調(diào)水影響)時間序列。

    對降雨時間序列進行STL分解,獲取降雨量趨勢變化、季節(jié)性變化及隨機波動項,結(jié)果見圖3。從圖中的趨勢變化可以看出,降雨量從2012年開始有減少的趨勢,整體呈現(xiàn)穩(wěn)定中略有波動的狀態(tài);降雨時間序列表現(xiàn)出顯著的周期性(周期為12個月),且從2009年開始周期內(nèi)的季節(jié)性變化均呈遞減的趨勢;從隨機波動項可以看出,殘差波動程度在1998—2013年的時間序列要大于其他時間序列。

    圖3 1980—2019年米山水庫月均降雨量時間序列分解結(jié)果

    對月均電導率時間序列進行STL分解,獲取其趨勢變化、季節(jié)性變化及隨機波動項,結(jié)果見圖4。從圖4(a)2012—2019年月均電導率STL分解結(jié)果可以看出,電導率值在2016年之前呈現(xiàn)穩(wěn)定中略有波動的態(tài)勢,電導率波動范圍為250~350 μS/cm,此時電導率主要受氣候變化影響,即處于自然狀態(tài)。在水庫2016年開始調(diào)水后電導率迅速增大并在2017年達到最大值,波動范圍為857~977 μS/cm,之后處于穩(wěn)定波動狀態(tài),此時電導率同時受氣候變化和區(qū)外調(diào)水的影響,即電導率處于混合狀態(tài)。電導率值的波動呈現(xiàn)出明顯的周期性,周期為12個月。從隨機波動項可以看出,殘差的波動程度在2016年后的時間序列要明顯大于之前的時間序列,說明米山水庫水質(zhì)在2016年實施區(qū)外調(diào)水后受人類活動等外部因素的干擾有所增加。為了更清楚地了解自然狀態(tài)下電導率的變化特征,對2012—2015年月電導率時間序列進行STL分解,如圖4(b)所示。由圖可以看出,自然狀態(tài)下電導率處于平穩(wěn)波動狀態(tài)數(shù)值的波動同樣呈現(xiàn)12個月的周期性,且隨機波動項的波動程度在2013—2014年之間的時間序列要大于其他的時間序列。

    (a)2012—2019年月均電導率

    (b)2012—2015年月均電導率圖4 米山水庫月均電導率的時間序列分解結(jié)果

    2.2 自然狀態(tài)下電導率與降雨關(guān)系分析

    對自然狀態(tài)下電導率與降雨的關(guān)系進行分析。1987—2015年年均電導率與年均降雨量變化關(guān)系圖如圖5所示。由圖可以看出,電導率值與降雨量呈現(xiàn)明顯的負相關(guān)關(guān)系,當年均降雨量較大時,因降雨對水體的帶電荷物質(zhì)的稀釋作用而導致電導率年均值相應的較小: 1999年研究區(qū)降雨量出現(xiàn)較小值(360.6 mm),米山水庫電導率達到347 μS/cm,為歷年最大值;1990年降雨量出現(xiàn)較大值,為1 180.9 mm,電導率僅為232 μS/cm,為歷年最小值。由此可見,降雨對米山水庫電導率的稀釋作用較為明顯,進一步說明自然狀態(tài)下電導率主要受氣候變化的影響。

    圖5 1987—2015年米山水庫年均電導率與降雨量變化關(guān)系

    2.3 米山水庫電導率模擬計算

    通過上述分析可以看出,米山水庫電導率在2016年前處于自然狀態(tài),因此可以利用1987—2015年電導率時間序列通過GM(1,1)分解模型建立自然狀態(tài)電導率模型;運用該模型外推,獲得2016—2019年自然狀態(tài)下的電導率,利用同時間序列實測電導率減去自然狀態(tài)電導率,獲得區(qū)外調(diào)水對電導率增加的貢獻率,從而分別獲得氣候變化和區(qū)外調(diào)水對電導率增加的影響。

    2.3.1 電導率的趨勢項

    利用1987—2015年自然狀態(tài)下電導率值作為原始數(shù)據(jù),建立GM(1,1)模型并求得模擬模型為

    (10)

    2.3.2 電導率的周期波動項

    GM(1,1)模型的殘差值可以由式(5)求得,如圖7所示。利用式(6)帶截距的正弦曲線對殘差值進行擬合,得到殘差周期擬合曲線的函數(shù)表達式如下:

    圖6 灰色系統(tǒng)GM(1,1)模型對米山水庫電導率模擬值與實測值

    (11)

    由圖可以看出,所求得的正弦殘差周期擬合曲線與殘差值整體擬合較好,擬合精度較高,能夠反映殘差的周期性變化。

    2.3.3 電導率的隨機波動項

    根據(jù)式(9),合并電導率趨勢項、周期波動項和隨機波動變化項,可以對自然狀態(tài)下米山水庫電導率動態(tài)變化進行模擬預測,結(jié)果見圖8。由圖可以看出:米山水庫自然狀態(tài)下電導率值從1987年開始不斷增大,到2002年達到最高值;2002—2015年電導率值呈現(xiàn)減小的趨勢,從2016年開始電導率值又有增大的趨勢。同時,電導率模擬值與實際觀測值整體擬合較好,絕大部分觀測值均在模擬值的上、下限范圍內(nèi)。此外,模擬值與實測值一樣,都體現(xiàn)了明顯的上升和回落的周期性變化,說明GM(1,1)分解模型的模擬精度較高。

    圖7 米山水庫電導率殘差周期擬合曲線、隨機波動曲線及殘差值

    2.3.4 氣候變化和區(qū)外調(diào)水貢獻度定量評價

    利用米山水庫2016—2019年電導率實測值減去自然狀態(tài)下的電導率值,即可得到區(qū)外調(diào)水造成的電導率變化值,結(jié)果見表2。從表中數(shù)據(jù)可以看出,2016、2017、2018、2019年自然狀態(tài)下米山水庫的電導率分別為269.3~317.26、270.94~318.9、273.4~321.36、276.63~324.59 μS/cm,處于正常的波動范圍內(nèi);區(qū)外調(diào)水造成的電導率的增加量分別為540.41~588.37、643.70~691.66、656.18~704.14、629.74~677.70 μS/cm,分別是自然狀態(tài)下電導率的1.85~2.01、2.17~2.38、2.19~2.40、2.09~2.28倍,說明區(qū)外調(diào)水對米山水庫電導率的增加起著主要控制作用。通過對造成米山水庫電導率變化的主控因素進行定量研究表明,影響米山水庫水質(zhì)變化的主要原因不僅僅是氣候變化,這一結(jié)論對于米山水庫建立更加全面的水質(zhì)安全風險管理體系,實施更加科學的水量-水質(zhì)聯(lián)合調(diào)控決策,維護供水安全具有重要意義。

    圖8 灰色系統(tǒng)GM(1,1)分解模型對米山水庫電導率的模擬值與實測值

    表2 氣候變化和區(qū)外調(diào)水對米山水庫電導率的影響

    2.3.5 分析討論

    本文中利用GM(1,1)分解模型建立了天然狀態(tài)下電導率的模擬預測模型,分別對氣候變化和區(qū)外調(diào)水對米山水庫電導率的影響進行了定量評價,研究結(jié)果具有重要的參考價值,但仍然存在著以下不足:

    1)米山水庫從2016年才開始實施區(qū)外調(diào)水,較短的區(qū)外調(diào)水時間序列不能充分反映區(qū)外調(diào)水影響下電導率的變化。

    2)GM(1,1)分解模型雖然具有用法簡單、所需實測數(shù)據(jù)少和比傳統(tǒng)的GM(1, 1)模型模擬精度更高的優(yōu)勢;但湖庫水體中的污染物濃度受到多種因素的影響(如降雨、氣溫和調(diào)水量等),因此水質(zhì)預測模型涉及多種因素的非線性關(guān)系。GM(1, 1)分解模型無法充分考慮多種因素的影響,造成一定誤差。

    3)在今后的研究中一方面要增加電導率時間序列的長度,另一方面要采用長短期記憶(LSTM)神經(jīng)網(wǎng)絡(luò)等多參數(shù)機器學習方法,通過輸入多個影響變量(如降雨量、氣溫和調(diào)水量等)進行水質(zhì)預測,由于輸入的不同變量之間可以相互驗證、相互作用,因此能夠有效提高預測精度。

    3 結(jié)論

    本文中首先利用STL法對米山水庫降雨量與電導率進行時間序列分解,從宏觀角度揭示它們的趨勢變化、周期波動長度和隨機波動范圍,并提取出電導率的自然狀態(tài)和混合狀態(tài)時間序列;在此基礎(chǔ)上,通過GM(1,1)分解模型建立自然狀態(tài)下電導率模擬模型,并將其外推,定量評價氣候變化和區(qū)外調(diào)水對米山水庫水質(zhì)演變的貢獻,得出以下結(jié)論:

    1)通過STL分解結(jié)果可以看出,米山水庫降雨時間序列雖然從2012年開始有減少的趨勢,但整體上仍然呈現(xiàn)穩(wěn)定中略有波動的狀態(tài);降雨表現(xiàn)出顯著的周期性。電導率在2016年之前同樣呈現(xiàn)穩(wěn)定中略有波動的狀態(tài),此時電導率主要受氣候變化影響;2016年區(qū)外調(diào)水后電導率快速增大并在2017年達到最大值,之后處于穩(wěn)定波動狀態(tài),此時電導率同時受氣候變化和區(qū)外調(diào)水的影響。電導率值的波動同樣呈現(xiàn)明顯的周期性。

    2)GM(1,1)分解模型的電導率模擬值與實測值整體擬合較好,絕大部分實測值均在模擬值的上、下限范圍內(nèi),且模擬值與實測值都能反映明顯的上升和回落的周期性變化,說明GM(1,1)分解模型的模擬精度較高。

    3)模擬分析結(jié)果顯示,2016—2019年區(qū)外調(diào)水造成米山水庫電導率的增加量分別為540.41~588.37、643.7~691.66、656.18~704.14、629.74~677.7 μS/cm,分別是自然狀態(tài)下電導率的1.85~2.01、2.17~2.38、2.19~2.40、2.09~2.28倍,說明區(qū)外調(diào)水對米山水庫電導率的增加起到主要控制作用。

    4)由于米山水庫水質(zhì)的演變主要受區(qū)外調(diào)水量、客水水質(zhì)影響,因此應建立更加全面的區(qū)外調(diào)水水質(zhì)安全風險管理體系,實施更加科學的水量-水質(zhì)聯(lián)合調(diào)控決策,以此維護米山水庫的供水安全和水生態(tài)穩(wěn)定。

    猜你喜歡
    水質(zhì)模型
    一半模型
    水質(zhì)抽檢豈容造假
    環(huán)境(2023年5期)2023-06-30 01:20:01
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    一月冬棚養(yǎng)蝦常見水質(zhì)渾濁,要如何解決?這9大原因及處理方法你要知曉
    這條魚供不應求!蝦蟹養(yǎng)殖戶、垂釣者的最愛,不用投喂,還能凈化水質(zhì)
    圖像識別在水質(zhì)檢測中的應用
    電子制作(2018年14期)2018-08-21 01:38:16
    3D打印中的模型分割與打包
    濟下水庫徑流水質(zhì)和垂向水質(zhì)分析及評價
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    国产毛片a区久久久久| 少妇的逼好多水| 亚洲最大成人av| 18禁裸乳无遮挡免费网站照片| 桃色一区二区三区在线观看| 久久亚洲国产成人精品v| 国产极品精品免费视频能看的| 久久亚洲精品不卡| av在线播放精品| 午夜激情欧美在线| 亚洲人成网站在线观看播放| 波多野结衣高清无吗| 亚洲熟妇中文字幕五十中出| 一级毛片我不卡| 中文字幕免费在线视频6| 国产伦一二天堂av在线观看| 美女大奶头视频| 色哟哟·www| 亚洲四区av| 国产久久久一区二区三区| 成人永久免费在线观看视频| 久久国产乱子免费精品| 两个人视频免费观看高清| 国产69精品久久久久777片| 午夜免费男女啪啪视频观看| 日韩精品有码人妻一区| 欧美在线一区亚洲| 成熟少妇高潮喷水视频| h日本视频在线播放| av国产免费在线观看| 超碰av人人做人人爽久久| 少妇熟女欧美另类| 小蜜桃在线观看免费完整版高清| 日本免费一区二区三区高清不卡| 亚洲精品粉嫩美女一区| 男人和女人高潮做爰伦理| 国产在视频线在精品| 亚洲色图av天堂| 此物有八面人人有两片| 99九九线精品视频在线观看视频| 免费看日本二区| 国产伦一二天堂av在线观看| 午夜a级毛片| 女人十人毛片免费观看3o分钟| 久久精品国产99精品国产亚洲性色| 亚洲国产日韩欧美精品在线观看| 深夜a级毛片| 亚洲无线在线观看| 婷婷色av中文字幕| 国产毛片a区久久久久| 国产一区二区在线观看日韩| 亚洲高清免费不卡视频| 直男gayav资源| 精品久久久久久久久久免费视频| 99热这里只有是精品50| 精品99又大又爽又粗少妇毛片| 一本精品99久久精品77| 美女被艹到高潮喷水动态| 亚洲精华国产精华液的使用体验 | 国产不卡一卡二| 色综合亚洲欧美另类图片| 亚洲国产欧洲综合997久久,| or卡值多少钱| 六月丁香七月| 亚洲成人久久爱视频| 亚洲在线观看片| 麻豆久久精品国产亚洲av| 国产高清视频在线观看网站| 久久精品国产亚洲网站| 亚洲无线观看免费| 男女下面进入的视频免费午夜| 丝袜美腿在线中文| 亚洲天堂国产精品一区在线| 国产高潮美女av| 变态另类丝袜制服| 乱系列少妇在线播放| 看片在线看免费视频| 男的添女的下面高潮视频| 久久久久久久午夜电影| 一级黄片播放器| 亚洲国产日韩欧美精品在线观看| 亚洲成人中文字幕在线播放| 免费人成视频x8x8入口观看| 老女人水多毛片| 高清午夜精品一区二区三区 | 国产激情偷乱视频一区二区| 村上凉子中文字幕在线| 毛片一级片免费看久久久久| 成人美女网站在线观看视频| 神马国产精品三级电影在线观看| 一区二区三区四区激情视频 | 一边亲一边摸免费视频| 午夜激情欧美在线| 亚洲在久久综合| www日本黄色视频网| 色播亚洲综合网| 亚洲国产精品久久男人天堂| 国产极品精品免费视频能看的| 变态另类丝袜制服| 成人鲁丝片一二三区免费| 国内精品一区二区在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产精品合色在线| 国产精品av视频在线免费观看| 亚洲婷婷狠狠爱综合网| 99热这里只有是精品50| 日产精品乱码卡一卡2卡三| 在线观看美女被高潮喷水网站| 久久精品国产99精品国产亚洲性色| 大香蕉久久网| 九九久久精品国产亚洲av麻豆| 高清毛片免费观看视频网站| 欧美高清成人免费视频www| 免费看a级黄色片| 一区二区三区四区激情视频 | 国产一区二区在线av高清观看| 国产黄色视频一区二区在线观看 | 欧美日本亚洲视频在线播放| 精品久久久久久久末码| 成人美女网站在线观看视频| 亚洲成人精品中文字幕电影| 久久久久九九精品影院| 亚洲欧美成人精品一区二区| 男女做爰动态图高潮gif福利片| 最后的刺客免费高清国语| 日韩大尺度精品在线看网址| 亚洲欧美日韩无卡精品| 午夜免费激情av| 精品国内亚洲2022精品成人| 黄片无遮挡物在线观看| 性色avwww在线观看| 99热精品在线国产| 久久热精品热| 免费人成视频x8x8入口观看| 看片在线看免费视频| 亚洲国产精品sss在线观看| 亚洲aⅴ乱码一区二区在线播放| 国内久久婷婷六月综合欲色啪| 欧美日韩一区二区视频在线观看视频在线 | 一本久久中文字幕| 亚洲成av人片在线播放无| 国产 一区精品| 天天躁日日操中文字幕| 成人午夜精彩视频在线观看| 一区福利在线观看| 熟女人妻精品中文字幕| 国产精品国产高清国产av| 日韩欧美一区二区三区在线观看| 国产中年淑女户外野战色| 国产免费男女视频| 免费不卡的大黄色大毛片视频在线观看 | 国产精品人妻久久久影院| 国产片特级美女逼逼视频| 精品无人区乱码1区二区| 国产成人午夜福利电影在线观看| 日本欧美国产在线视频| 久久6这里有精品| 精品久久久久久久久久久久久| 国产精品一区二区性色av| 亚洲精品国产成人久久av| 久久草成人影院| 亚洲人成网站在线播放欧美日韩| 综合色av麻豆| 亚洲久久久久久中文字幕| 精品熟女少妇av免费看| 一夜夜www| 久久久久久久久久久丰满| 不卡一级毛片| 亚洲18禁久久av| 99久久无色码亚洲精品果冻| 精品一区二区三区视频在线| 国产色婷婷99| 在线观看一区二区三区| 在线免费观看不下载黄p国产| 免费看光身美女| 免费电影在线观看免费观看| 秋霞在线观看毛片| 高清午夜精品一区二区三区 | 成人高潮视频无遮挡免费网站| 99热只有精品国产| 成人二区视频| 日韩人妻高清精品专区| 亚洲va在线va天堂va国产| 国产午夜精品一二区理论片| 久久久久久久久大av| 午夜激情欧美在线| 日本黄色视频三级网站网址| 18禁在线播放成人免费| 亚洲欧美日韩高清在线视频| 久久精品国产自在天天线| 国产亚洲精品久久久com| 亚洲成av人片在线播放无| 男女做爰动态图高潮gif福利片| 亚洲国产色片| av天堂中文字幕网| 国产黄a三级三级三级人| 国产一级毛片在线| 三级国产精品欧美在线观看| av卡一久久| 免费一级毛片在线播放高清视频| av国产免费在线观看| 国产老妇伦熟女老妇高清| 美女被艹到高潮喷水动态| 一级黄色大片毛片| 3wmmmm亚洲av在线观看| 免费看a级黄色片| 老师上课跳d突然被开到最大视频| 免费黄网站久久成人精品| 赤兔流量卡办理| 日韩精品青青久久久久久| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩高清专用| 十八禁国产超污无遮挡网站| 国产精品1区2区在线观看.| 我的女老师完整版在线观看| 日韩高清综合在线| 乱码一卡2卡4卡精品| 国产高清不卡午夜福利| 国产精品一及| 久久国内精品自在自线图片| 老司机福利观看| 人妻制服诱惑在线中文字幕| www日本黄色视频网| 亚洲av第一区精品v没综合| 免费看日本二区| 亚洲最大成人中文| 亚洲精品成人久久久久久| 欧美性猛交╳xxx乱大交人| 久久精品久久久久久久性| 国产成人一区二区在线| 亚洲精品日韩av片在线观看| 91久久精品国产一区二区成人| 九九爱精品视频在线观看| 国产精品.久久久| 国产伦精品一区二区三区视频9| 精品久久久久久久末码| 观看免费一级毛片| 麻豆久久精品国产亚洲av| 亚洲内射少妇av| 久久精品夜色国产| 日韩一区二区视频免费看| 午夜福利高清视频| 九九热线精品视视频播放| 国产美女午夜福利| 特级一级黄色大片| 一区福利在线观看| 一级毛片我不卡| 国产成人精品一,二区 | 少妇裸体淫交视频免费看高清| 国产精品一二三区在线看| 午夜免费男女啪啪视频观看| av天堂在线播放| a级毛片a级免费在线| 国产成人aa在线观看| 美女被艹到高潮喷水动态| 欧美+亚洲+日韩+国产| 乱人视频在线观看| 国产精品久久久久久久久免| 不卡视频在线观看欧美| 亚洲乱码一区二区免费版| 我要搜黄色片| 99久国产av精品| 婷婷精品国产亚洲av| 亚洲18禁久久av| 1024手机看黄色片| 国产精品野战在线观看| 男的添女的下面高潮视频| 亚洲精品国产成人久久av| 桃色一区二区三区在线观看| 精品少妇黑人巨大在线播放 | 国产伦精品一区二区三区视频9| 久久韩国三级中文字幕| 午夜久久久久精精品| 亚洲中文字幕日韩| 黑人高潮一二区| 非洲黑人性xxxx精品又粗又长| videossex国产| 国产色爽女视频免费观看| 亚洲电影在线观看av| 中文精品一卡2卡3卡4更新| 乱系列少妇在线播放| 国产精品不卡视频一区二区| av免费观看日本| 亚洲精品影视一区二区三区av| 欧洲精品卡2卡3卡4卡5卡区| 国产黄片美女视频| 高清日韩中文字幕在线| 日韩亚洲欧美综合| 在线国产一区二区在线| 亚洲精品成人久久久久久| videossex国产| 国产成人a区在线观看| 亚洲精品亚洲一区二区| 国产精品女同一区二区软件| 国产精品蜜桃在线观看 | 1000部很黄的大片| 亚洲18禁久久av| 日本欧美国产在线视频| 成人欧美大片| 日韩人妻高清精品专区| 国产黄片视频在线免费观看| 国产乱人偷精品视频| 麻豆久久精品国产亚洲av| 在线观看美女被高潮喷水网站| 国产又黄又爽又无遮挡在线| 成人毛片a级毛片在线播放| 一个人免费在线观看电影| 欧洲精品卡2卡3卡4卡5卡区| 精品不卡国产一区二区三区| 亚洲精品色激情综合| 国产麻豆成人av免费视频| 欧美高清成人免费视频www| 久久国产乱子免费精品| 久久韩国三级中文字幕| 国内精品久久久久精免费| 91久久精品电影网| 日日摸夜夜添夜夜爱| 亚洲丝袜综合中文字幕| 亚洲自拍偷在线| 日韩欧美精品免费久久| 变态另类丝袜制服| 成年女人永久免费观看视频| 一级av片app| 在线a可以看的网站| 欧美精品国产亚洲| 亚洲欧美精品综合久久99| h日本视频在线播放| 中文精品一卡2卡3卡4更新| 青春草国产在线视频 | 最新中文字幕久久久久| 男人的好看免费观看在线视频| 亚洲丝袜综合中文字幕| 波多野结衣高清无吗| 丝袜美腿在线中文| 日韩大尺度精品在线看网址| 人人妻人人澡欧美一区二区| 欧美激情在线99| 精品久久久久久久久久久久久| 悠悠久久av| 精品免费久久久久久久清纯| 亚洲国产精品国产精品| 久久热精品热| 日韩av在线大香蕉| 国产色婷婷99| 大型黄色视频在线免费观看| 两个人视频免费观看高清| 国产高清三级在线| 男人舔奶头视频| 亚洲精品影视一区二区三区av| 可以在线观看毛片的网站| eeuss影院久久| 亚洲四区av| 国产伦理片在线播放av一区 | 嫩草影院入口| 成人一区二区视频在线观看| 18禁黄网站禁片免费观看直播| 春色校园在线视频观看| 中国美白少妇内射xxxbb| www.色视频.com| 久久久国产成人免费| 成人高潮视频无遮挡免费网站| 久久久久久久久久黄片| 日韩一区二区三区影片| 免费av不卡在线播放| 成年女人永久免费观看视频| 欧美成人精品欧美一级黄| 在线天堂最新版资源| 欧美一区二区国产精品久久精品| 国产精品蜜桃在线观看 | 亚洲av免费在线观看| 国产视频内射| 精品熟女少妇av免费看| 18禁在线播放成人免费| 精品久久久久久久久久久久久| 亚洲av中文av极速乱| 夜夜爽天天搞| 国产精品电影一区二区三区| eeuss影院久久| 日本欧美国产在线视频| 熟妇人妻久久中文字幕3abv| 美女黄网站色视频| 麻豆精品久久久久久蜜桃| 国产精品1区2区在线观看.| 亚洲精品国产成人久久av| 在线播放国产精品三级| 男女啪啪激烈高潮av片| 岛国在线免费视频观看| 国产色婷婷99| 看非洲黑人一级黄片| 国产老妇伦熟女老妇高清| 伦精品一区二区三区| 在线免费观看的www视频| 少妇人妻精品综合一区二区 | 久久人人爽人人爽人人片va| 精品人妻熟女av久视频| 欧美激情久久久久久爽电影| 国产一区二区在线av高清观看| 晚上一个人看的免费电影| 亚洲国产高清在线一区二区三| 免费大片18禁| 久久久精品大字幕| 久久精品久久久久久噜噜老黄 | 午夜激情欧美在线| 在线观看66精品国产| 久久久久久国产a免费观看| 午夜精品国产一区二区电影 | 亚洲精品亚洲一区二区| 午夜老司机福利剧场| 神马国产精品三级电影在线观看| 性插视频无遮挡在线免费观看| 99久久精品一区二区三区| av国产免费在线观看| 国产真实乱freesex| 国产精品一区二区三区四区免费观看| 国产av不卡久久| 直男gayav资源| 日本五十路高清| 久久久色成人| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美三级三区| 夜夜夜夜夜久久久久| 身体一侧抽搐| 嫩草影院精品99| 观看免费一级毛片| 18禁黄网站禁片免费观看直播| 国产真实乱freesex| 国产精品.久久久| 久久人人爽人人爽人人片va| 精品99又大又爽又粗少妇毛片| 波多野结衣高清无吗| 禁无遮挡网站| 亚洲av二区三区四区| 自拍偷自拍亚洲精品老妇| 日本欧美国产在线视频| 成人高潮视频无遮挡免费网站| 成人特级av手机在线观看| 国产精品久久电影中文字幕| 国产精品福利在线免费观看| 老司机影院成人| 少妇高潮的动态图| 国产精品日韩av在线免费观看| 一级毛片aaaaaa免费看小| 变态另类成人亚洲欧美熟女| 狂野欧美激情性xxxx在线观看| 最近最新中文字幕大全电影3| 少妇人妻一区二区三区视频| 欧美色视频一区免费| 色综合色国产| 狂野欧美激情性xxxx在线观看| 最近最新中文字幕大全电影3| 久久精品久久久久久噜噜老黄 | 99热全是精品| 一进一出抽搐gif免费好疼| 舔av片在线| 啦啦啦观看免费观看视频高清| 日韩欧美一区二区三区在线观看| 国产v大片淫在线免费观看| 精品人妻偷拍中文字幕| 欧美zozozo另类| 亚洲,欧美,日韩| 精品久久国产蜜桃| 中文字幕人妻熟人妻熟丝袜美| 日韩中字成人| 能在线免费观看的黄片| 久久久久国产网址| 成人特级av手机在线观看| 韩国av在线不卡| 成人特级黄色片久久久久久久| 99久久无色码亚洲精品果冻| 欧美人与善性xxx| 噜噜噜噜噜久久久久久91| 99热精品在线国产| 国产精品一及| 久久精品国产亚洲av天美| 自拍偷自拍亚洲精品老妇| 午夜激情欧美在线| 欧美精品一区二区大全| 天天躁日日操中文字幕| 亚洲在久久综合| 我要搜黄色片| 97超碰精品成人国产| 国产精品三级大全| 丝袜喷水一区| 男女下面进入的视频免费午夜| 国产精品一二三区在线看| 神马国产精品三级电影在线观看| 一本久久精品| 国语自产精品视频在线第100页| 日日干狠狠操夜夜爽| 中文资源天堂在线| 噜噜噜噜噜久久久久久91| 精品国内亚洲2022精品成人| 五月伊人婷婷丁香| 桃色一区二区三区在线观看| 欧美又色又爽又黄视频| 狂野欧美白嫩少妇大欣赏| 一进一出抽搐gif免费好疼| 春色校园在线视频观看| 插逼视频在线观看| 国产精品久久视频播放| 亚洲经典国产精华液单| 女人十人毛片免费观看3o分钟| 内射极品少妇av片p| 嫩草影院精品99| 国产伦精品一区二区三区四那| 床上黄色一级片| 99精品在免费线老司机午夜| 亚洲国产欧美在线一区| 国产成人精品一,二区 | 久久精品国产亚洲av天美| 色噜噜av男人的天堂激情| 极品教师在线视频| 精品人妻视频免费看| 青春草国产在线视频 | 激情 狠狠 欧美| 国产精品人妻久久久久久| 淫秽高清视频在线观看| 精品久久久久久久人妻蜜臀av| 国产一区二区激情短视频| 国产欧美日韩精品一区二区| 国产黄a三级三级三级人| 久久这里有精品视频免费| 九色成人免费人妻av| 美女被艹到高潮喷水动态| АⅤ资源中文在线天堂| 亚洲丝袜综合中文字幕| 国产精品日韩av在线免费观看| 99久久精品一区二区三区| 精品久久国产蜜桃| 可以在线观看毛片的网站| 久久久久久国产a免费观看| 国产伦精品一区二区三区四那| 欧美一区二区国产精品久久精品| 亚洲人成网站在线观看播放| 欧美一区二区国产精品久久精品| 天堂av国产一区二区熟女人妻| 一级毛片电影观看 | 一本一本综合久久| 色5月婷婷丁香| 激情 狠狠 欧美| 日日摸夜夜添夜夜添av毛片| 99九九线精品视频在线观看视频| 尾随美女入室| 免费观看人在逋| 日本熟妇午夜| 日日摸夜夜添夜夜爱| 看十八女毛片水多多多| 夜夜爽天天搞| 国产高清不卡午夜福利| 高清毛片免费观看视频网站| 亚洲第一区二区三区不卡| 亚洲欧美日韩高清专用| 三级国产精品欧美在线观看| 麻豆国产av国片精品| 亚洲精品乱码久久久久久按摩| 精品人妻熟女av久视频| 亚洲精品乱码久久久久久按摩| 久久久久性生活片| 亚洲欧美精品专区久久| 一级av片app| 久久亚洲精品不卡| 老师上课跳d突然被开到最大视频| 亚洲精品色激情综合| 可以在线观看毛片的网站| 国产av不卡久久| 免费一级毛片在线播放高清视频| 嫩草影院精品99| 亚洲在久久综合| 久久6这里有精品| 免费av不卡在线播放| 噜噜噜噜噜久久久久久91| 搞女人的毛片| 亚洲aⅴ乱码一区二区在线播放| 免费av不卡在线播放| 亚洲va在线va天堂va国产| 身体一侧抽搐| 午夜激情欧美在线| 国产精品福利在线免费观看| 特大巨黑吊av在线直播| 在线观看一区二区三区| 久久精品国产亚洲av天美| 99久久中文字幕三级久久日本| 又黄又爽又刺激的免费视频.| 国产国拍精品亚洲av在线观看| 国产精品人妻久久久久久| 免费大片18禁| 看十八女毛片水多多多| 日本一二三区视频观看| 天天躁夜夜躁狠狠久久av| 欧美3d第一页| 亚洲欧美成人综合另类久久久 | 午夜精品国产一区二区电影 | 少妇丰满av| 午夜久久久久精精品| 国产国拍精品亚洲av在线观看| 18+在线观看网站| 久久精品国产亚洲av香蕉五月| 级片在线观看| 哪里可以看免费的av片| 中文资源天堂在线| 亚洲aⅴ乱码一区二区在线播放| 高清午夜精品一区二区三区 | 人妻夜夜爽99麻豆av| 51国产日韩欧美| 婷婷色综合大香蕉| 九九爱精品视频在线观看| 99热只有精品国产| 精品久久久久久久人妻蜜臀av| 午夜亚洲福利在线播放| 成年版毛片免费区| 毛片女人毛片| 国产av在哪里看| 欧洲精品卡2卡3卡4卡5卡区| 国产成人精品婷婷| 18+在线观看网站|