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

    基于人工神經(jīng)網(wǎng)絡(luò)模型的麗江黑龍?zhí)度簲嗔黝A(yù)測(cè)

    2020-11-06 05:53李恒麗李保珠
    中國(guó)水運(yùn) 2020年7期
    關(guān)鍵詞:人工神經(jīng)網(wǎng)絡(luò)

    李恒麗 李保珠

    摘 要:位于云南省麗江市古城區(qū)的黑龍?zhí)度?,在降雨偏少時(shí)會(huì)發(fā)生季節(jié)性斷流,特別是近年來(lái)由于降雨減少,人為活動(dòng)增多,泉群斷流情況越來(lái)越嚴(yán)重,斷流次數(shù)逐漸增多、時(shí)長(zhǎng)也增長(zhǎng),以致麗江市從“高原姑蘇”變成嚴(yán)重缺水的城市,直接影響了麗江市的形象,制約其社會(huì)經(jīng)濟(jì)的發(fā)展。因此對(duì)黑龍?zhí)度哼M(jìn)行動(dòng)態(tài)研究與斷流分析預(yù)測(cè),對(duì)于麗江市的發(fā)展、地下水資源的利用具有重要意義。本文采用人工神經(jīng)網(wǎng)絡(luò)方法分析黑龍?zhí)度航涤炅颗c流量的變化,應(yīng)用python語(yǔ)言及相關(guān)庫(kù)分析,計(jì)算黑龍?zhí)督涤炅颗c流量的相關(guān)性系數(shù),利用matlab人工神經(jīng)網(wǎng)絡(luò)方法預(yù)測(cè)黑龍?zhí)度旱某鏊考芭袛嗍欠駭嗔?。模型模擬及預(yù)測(cè)結(jié)果較好,與實(shí)際值相差無(wú)幾,能較好地預(yù)測(cè)流量的變化趨勢(shì)及斷流結(jié)果。

    關(guān)鍵詞:人工神經(jīng)網(wǎng)絡(luò);流量預(yù)測(cè);斷流預(yù)測(cè)

    中圖分類號(hào):TP183 ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼:A ? ? ? ? ? ?文章編號(hào):1006—7973(2020)07-0149-04

    1 引言

    黑龍?zhí)度何挥邴惤璧貣|北部的象山山麓,九子海復(fù)向斜南部翹起端,泉群出露于三疊系北衙組中上段裂隙發(fā)育的灰?guī)r、白云質(zhì)灰?guī)r中。泉群地下水類型為炭酸鹽巖裂隙水,補(bǔ)給來(lái)源主要是大氣降水,其次還有冰雪融水、基巖地下水、河流渠道水等,地下水補(bǔ)給來(lái)源十分豐富,使黑龍?zhí)度簱碛休^大的出水量[1~4]。黑龍?zhí)度簽辂惤懦翘峁┥?、景觀及服務(wù)業(yè)用水,其流量的變化直接關(guān)系到麗江市民的生活狀況。但上世紀(jì)八十年代以來(lái),其發(fā)生過多次斷流,給麗江的生產(chǎn)、生活各方面造成了巨大的影響,所以其流量的穩(wěn)定性及提前預(yù)測(cè)具有重要意義[5,6]。

    目前研究發(fā)現(xiàn),對(duì)黑龍?zhí)度毫髁康闹饕绊懸蛩貫椋核囱a(bǔ)給區(qū)域的降雨量、城市用水量、植被退化與地下水開采。其中降雨量為最主要因素,降雨量減少導(dǎo)致滲入地下的水流減少,最終造成黑龍?zhí)稊嗔鱗7]。目前對(duì)黑龍?zhí)读髁康难芯空撐囊灿泻芏啵饕酝A粼诙ㄐ苑治錾?,只是研究了大體的影響因素及如何作用還有一些統(tǒng)計(jì)學(xué)數(shù)據(jù);對(duì)流量的定量研究仍處于初步階段,很多模型都沒有嘗試過,當(dāng)然,對(duì)其進(jìn)行定量研究時(shí),整個(gè)系統(tǒng)的周期性可能不是按整月的,而數(shù)據(jù)統(tǒng)計(jì)時(shí)都按整月統(tǒng)計(jì),造成了變化不一致等;還有整個(gè)系統(tǒng)的響應(yīng)并不是線性的,自變量較多及存在很多潛在因素的影響等造成系統(tǒng)建模難度增大。而利用傳統(tǒng)的GMS進(jìn)行建模分析時(shí),需要鉆孔鉆探試驗(yàn)確認(rèn)地下含水層的厚度、頂?shù)装?,各處的滲透系數(shù)等,而這些數(shù)據(jù)因?yàn)楦鞣矫娴脑蚨紱]進(jìn)行實(shí)驗(yàn)采集,所以該方法也很難進(jìn)行。

    在這種模糊條件過多,響應(yīng)為非線性,各因素的影響系數(shù)未知的系統(tǒng)中,最合適的研究方法為人工神經(jīng)網(wǎng)絡(luò)。該方法在確定主要影響因素的情況下,可以將所有影響因素作為自變量導(dǎo)入,將要研究的對(duì)象作為因變量,讓計(jì)算機(jī)自動(dòng)尋找規(guī)律,而且如果對(duì)結(jié)果不滿意可以多次訓(xùn)練模型或者調(diào)整參數(shù)、自變量等,使結(jié)果最終符合預(yù)期目標(biāo),達(dá)到良好的預(yù)測(cè)效果。利用python先分析出各月降雨量、井水位與流量的相關(guān)性,確定主要影響因素,然后用matlab的神經(jīng)網(wǎng)絡(luò)工具箱,代入所有影響因素跟是否斷流后,訓(xùn)練模型,再用新數(shù)據(jù)檢驗(yàn)?zāi)P停粩嗾{(diào)整參數(shù),當(dāng)效果達(dá)到預(yù)期后最終導(dǎo)出模型。

    2 主要影響因素的確定

    當(dāng)前能為建模提供建?;A(chǔ)的數(shù)據(jù)主要參考為麗江站1988~2018年每月降雨量,九子海2000~2017年每月降雨量,麗江站1992~2018年每月井水位,黑龍?zhí)?990~2012年每月流量。由于九子海降雨量太小,為間接影響黑龍?zhí)读髁恳蛩厍遗c黑龍?zhí)读髁繑?shù)據(jù)的時(shí)間重合度太低,暫時(shí)不作為因素考慮。為了宏觀地觀察麗江站降雨量、麗江井水位與黑龍?zhí)读髁康年P(guān)系,先繪制出重合時(shí)間內(nèi)三者的趨勢(shì)圖。由于井水位的數(shù)據(jù)級(jí)別太小,為方便觀察,將其擴(kuò)大400倍與降雨量跟流量比較。

    圖1明顯表明:流量與井水位保持同步變化,當(dāng)井水位上升時(shí)流量隨之上升,且時(shí)間刻度上保持一致,未見滯后性。而降雨量與流量整體趨勢(shì)也相同,但圖上也很明顯有時(shí)間差,當(dāng)降雨量到達(dá)最大值時(shí),往往過二至三個(gè)月,流量才達(dá)到最大值,由此猜測(cè)本月麗江降雨量并不會(huì)對(duì)本月流量產(chǎn)生影響,而是影響兩到三月后黑龍?zhí)兜牧髁俊?/p>

    為了驗(yàn)證以上猜想,需要計(jì)算降雨量與流量的相關(guān)性系數(shù)。

    經(jīng)計(jì)算,本月流量與本月降雨量相關(guān)性系數(shù)為-0.166067,與往前一個(gè)月為0.110125,與往前兩個(gè)月為0.333704,與往前三個(gè)月為0.449788,與往前四個(gè)月為0.456936,本月流量與本月井水位的相關(guān)性系數(shù)為0.938674。很明顯,相關(guān)性系數(shù)的計(jì)算結(jié)果與假設(shè)相符。

    由于表中所有數(shù)據(jù)都是基于時(shí)間序列的離散型數(shù)據(jù),且流量與本月降雨,前月降雨相關(guān)性系數(shù)太低,由此利用傳統(tǒng)的分析方法,假設(shè)流量的變化量與前2、3、4月降雨量存在相關(guān)性。將本月流量與前月流量作差,得到流量的變化量。由于多元自變量不方便分析,同時(shí)為了均衡多個(gè)月份降雨量對(duì)流量的影響,現(xiàn)計(jì)算多個(gè)月流量的平均值與流量的相關(guān)性,然后選取相關(guān)性最大的作為單一自變量。

    經(jīng)計(jì)算,本月流量與前2、3月降雨量平均值的相關(guān)性系數(shù)為0.435440,與前3、4月份為0.503288,與前2、3、4月為0.503063。因此此處選擇相關(guān)性系數(shù)最大值作為自變量,即前3、4月份的平均降水量。由此,將前3、4月份的平均降雨量作為自變量,將流量的變化量作為因變量進(jìn)行線性分析并繪制出分布散點(diǎn)圖。

    圖2中縱坐標(biāo)為本月相對(duì)于上月的流量變化量,橫坐標(biāo)為本月前3、4月的平均降雨量,圖中直線為基于最小二乘法的線性回歸線。很明顯,分布呈現(xiàn)很大的離散性,即表明相同或相近的降雨量所導(dǎo)致的流量變化差異很大,如前3、4月平均降雨量在1500mm左右時(shí),本月流量的變化可能幾乎不變,也可能上升800mm。另外,點(diǎn)的分布也并未與任何規(guī)律曲線(二次或三次曲線)相符,由此,并不能基于此建模。

    由以上推論,降雨量不能作為影響流量的單一變量,而井水位與流量的相關(guān)性系數(shù)比較大,所以可以將上個(gè)月的井水位作為一個(gè)影響因素。另外,從原理上分析,也是在降雨后,雨水滲入地下,三個(gè)月以后才影響黑龍?zhí)度毫髁浚晁德浜鬂B入地下,優(yōu)先會(huì)作用表現(xiàn)在井水位上,然后井水位上升后,再體現(xiàn)在流量的變化上。且井水位可以看成是降雨量與流量雙重作用下長(zhǎng)期累積的結(jié)果,其承擔(dān)著降雨量轉(zhuǎn)化為黑龍?zhí)读髁康臉屑~關(guān)系,所以應(yīng)將其看成是影響流量的重要因素。另外,前一個(gè)月流量也必然影響著本月的流量,因?yàn)榱髁慷际腔谇耙粋€(gè)月的流量在變化著的,例如,上月流量很大,但前幾個(gè)月降雨量都很小,流量必然會(huì)降低,這種降低是基于前個(gè)月的流量下的,卻不會(huì)因?yàn)榻涤炅繙p少而突然斷流。上個(gè)月流量較大時(shí),即使降雨量很小,斷流的幾率也幾乎為零。最后,由于每年降雨量,蒸發(fā)量,農(nóng)業(yè)灌溉等也是周期性的,所以將月份納入到影響因素也很有必要。由此,初步確認(rèn)了所有的自變量即影響流量的所有因素:前月流量,前月井水位,月份,往前三個(gè)月,往前四個(gè)月與往前五個(gè)月的降雨量。

    3 ?MATLAB神經(jīng)網(wǎng)絡(luò)分析

    人工神經(jīng)網(wǎng)絡(luò)是模仿人腦工作原理設(shè)計(jì)的計(jì)算模型,它由多個(gè)節(jié)點(diǎn)互相連接組成,能對(duì)數(shù)據(jù)的復(fù)雜關(guān)系進(jìn)行建模。節(jié)點(diǎn)間的連接都被賦予不同的權(quán)重,權(quán)重的含義是一個(gè)節(jié)點(diǎn)對(duì)另一個(gè)節(jié)點(diǎn)的影響大小。每個(gè)節(jié)點(diǎn)表示一個(gè)特定的函數(shù),其他節(jié)點(diǎn)傳來(lái)的數(shù)據(jù)到達(dá)該節(jié)點(diǎn)后會(huì)進(jìn)行相應(yīng)的權(quán)重計(jì)算,輸入到一個(gè)激活函數(shù)中并計(jì)算新的值。從系統(tǒng)的角度上看,人工神經(jīng)網(wǎng)絡(luò)是由大量的神經(jīng)元經(jīng)過復(fù)雜的連接構(gòu)成的自適應(yīng)非線性動(dòng)態(tài)系統(tǒng)。一個(gè)兩層的神經(jīng)網(wǎng)絡(luò)就可以逼近任意函數(shù),所以,只要有足夠的神經(jīng)元和訓(xùn)練數(shù)據(jù),人工神經(jīng)網(wǎng)絡(luò)可以學(xué)到很多復(fù)雜的函數(shù)。

    流量的預(yù)測(cè)是一個(gè)復(fù)雜的非線性系統(tǒng),這個(gè)系統(tǒng)的真正狀態(tài)方程必然是由很多微分方程組成,同時(shí)系統(tǒng)又包含了很多的未知情況,有很多隱藏的因素作用在系統(tǒng)上,而且支持整個(gè)系統(tǒng)的理論也有很多需要深究之處。綜上,要想用數(shù)學(xué)方法準(zhǔn)確建模是十分困難的,因此,可以嘗試用神經(jīng)網(wǎng)絡(luò)來(lái)表達(dá)這個(gè)復(fù)雜的非線性系統(tǒng)。

    將整個(gè)系統(tǒng)看成是一個(gè)黑箱,所有影響因素作為系統(tǒng)輸入,把是否斷流作為輸出,最終整個(gè)系統(tǒng)也就是成為了一個(gè)二分類問題。通過所有的統(tǒng)計(jì)數(shù)據(jù)神經(jīng)網(wǎng)絡(luò),不斷地調(diào)整優(yōu)化后,使得該網(wǎng)絡(luò)有了能表達(dá)該未知系統(tǒng)的能力,這個(gè)神經(jīng)網(wǎng)絡(luò)也就能預(yù)測(cè)系統(tǒng)輸出了。

    此處使用MATLAB的神經(jīng)網(wǎng)絡(luò)工具箱對(duì)其進(jìn)行分析,該神經(jīng)網(wǎng)絡(luò)是一種多層前饋神經(jīng)網(wǎng)絡(luò),該網(wǎng)絡(luò)最大的特點(diǎn)就是信號(hào)前向傳播,誤差反向傳播,每層的神經(jīng)元狀態(tài)只會(huì)影響下一層的,如果輸出層達(dá)不到期望的輸出,就會(huì)把誤差進(jìn)行反向傳播,調(diào)整神經(jīng)網(wǎng)絡(luò)的權(quán)值與閾值,最終使神經(jīng)網(wǎng)絡(luò)的輸出不斷逼近預(yù)期值。神經(jīng)網(wǎng)絡(luò)的訓(xùn)練過程包含以下幾個(gè)步驟:

    (1)初始化神經(jīng)網(wǎng)絡(luò):根據(jù)系統(tǒng)的輸入輸出值(X,Y)大致推出網(wǎng)絡(luò)輸入層節(jié)點(diǎn)數(shù)n、隱含層節(jié)點(diǎn)數(shù) ? ,輸出層節(jié)點(diǎn)數(shù)m,初始化輸入層、隱含層和輸出層神經(jīng)元之間的連接權(quán)值ωij,ωjk,初始化隱含層閾值α,輸出層閾值b,給定學(xué)習(xí)速率和神經(jīng)元激勵(lì)函數(shù)。

    對(duì)于預(yù)測(cè)斷流,將是否斷流看成是一個(gè)二分類問題,將斷流定義為0,未斷流定義為1,而神經(jīng)網(wǎng)絡(luò)所有的預(yù)測(cè)并不會(huì)完全確定各月是否斷流,而是計(jì)算出一個(gè)概率,該概率值介于0與1之間,預(yù)測(cè)結(jié)果的數(shù)據(jù)越接近1就表示越不會(huì)斷流,越接近0,即表示斷流的概率越大。

    將麗江站降雨量,月份,黑龍?zhí)渡显铝髁?,上月井水位作為自變量,設(shè)置神經(jīng)網(wǎng)絡(luò)的隱藏層深度為8層,設(shè)置訓(xùn)練集為90%的總體樣本數(shù),驗(yàn)證集為5%的總體樣本數(shù),測(cè)試集為5%的總體樣本數(shù),對(duì)本月是否斷流進(jìn)行分析。在對(duì)神經(jīng)網(wǎng)絡(luò)進(jìn)行多次調(diào)整并重新學(xué)習(xí)后,訓(xùn)練出了預(yù)測(cè)準(zhǔn)確率較高的模型。用該模型重新預(yù)測(cè)全部月份的斷流情況,總數(shù)據(jù)量為237個(gè),預(yù)測(cè)斷流的準(zhǔn)確次數(shù)為227次,準(zhǔn)確率達(dá)到了95.78%,而且在預(yù)測(cè)會(huì)斷流而實(shí)際未發(fā)生斷流時(shí),該月實(shí)際流量也已經(jīng)很低了,接近于斷流狀態(tài)。實(shí)際發(fā)生了斷流卻未預(yù)測(cè)成功的只有一次。截取所有數(shù)據(jù)中1995年的情況,表1中列出了該年每月的實(shí)際流量,實(shí)際是否斷流,預(yù)測(cè)的斷流情況。

    表1對(duì)1995年6、7兩月的斷流已經(jīng)準(zhǔn)確地預(yù)測(cè)出來(lái)了,但對(duì)5月的預(yù)測(cè)有誤差,預(yù)測(cè)該月斷流的概率較大,但實(shí)際還有流量。從實(shí)際情況來(lái)看,卻也不能認(rèn)為預(yù)測(cè)錯(cuò)誤,因?yàn)閺牧髁縼?lái)說(shuō),該月的流量已經(jīng)很低了,可能該月月底也已經(jīng)斷流了,實(shí)際也需要采取補(bǔ)水措施。綜上,該模型已經(jīng)能對(duì)斷流提供良好的預(yù)測(cè)能力,能夠指導(dǎo)進(jìn)行提前預(yù)防控制了。

    4 結(jié)論

    (1)經(jīng)計(jì)算得出,黑龍?zhí)度毫髁颗c降雨量的相關(guān)性系數(shù)最大的月份為本月往前三、四月:與往前三個(gè)月降雨量為0.449788,與往前四個(gè)月為0.456936,另外本月流量與本月井水位的相關(guān)性系數(shù)為0.938674。由此得出:①流量與井水位基本保持同步變化,當(dāng)井水位上升時(shí)流量隨之上升,未見滯后性;②本月麗江降雨量并不會(huì)對(duì)本月流量產(chǎn)生影響,而是影響三到四月后黑龍?zhí)兜牧髁俊?/p>

    (2)前幾月平均降雨量與本月流量相關(guān)性系數(shù)計(jì)算中,得出前3、4月份平均降雨量相關(guān)性最大,為0.503288。而將前3、4月份的平均降雨量與流量的變化量進(jìn)行線性分析,發(fā)現(xiàn)準(zhǔn)確性太差,從而得出降雨量不能作為影響流量的單一影響因素,并且綜合分析初步確認(rèn)了所有的影響流量的所有因素:前月流量,前月井水位,月份,往前三個(gè)月,往前四個(gè)月與往前五個(gè)月的降雨量。

    (3)通過人工神經(jīng)網(wǎng)絡(luò)建立的模型已經(jīng)能較準(zhǔn)確地預(yù)測(cè)是否斷流,能使用MATLAB中建好的模型對(duì)是否斷流進(jìn)行指導(dǎo)性的預(yù)測(cè);本文通過選取黑龍?zhí)度?992年到2011年降雨量與出水量數(shù)據(jù)的90%搭建人工神經(jīng)網(wǎng)絡(luò),用該模型預(yù)測(cè)剩余10%的情況,結(jié)果與預(yù)期目標(biāo)相符。

    (4)通過人工神經(jīng)網(wǎng)絡(luò)模型如果有更多的數(shù)據(jù)進(jìn)行訓(xùn)練,預(yù)測(cè)的準(zhǔn)確性和通用性將會(huì)大大提高,所以可以將以后統(tǒng)計(jì)的實(shí)際數(shù)據(jù)納入到訓(xùn)練樣本中,從而提高生產(chǎn)力。對(duì)于預(yù)測(cè)到斷流,該給黑龍?zhí)堆a(bǔ)多少水的問題,由于沒有補(bǔ)水對(duì)流量影響的試驗(yàn)數(shù)據(jù),所以很難給出參考,需要進(jìn)行相關(guān)實(shí)驗(yàn)及相關(guān)的地質(zhì)考察后,才有條件模擬出。

    參考文獻(xiàn):

    [1] 康曉波,張華,王波. 麗江黑龍?zhí)度核牡刭|(zhì)特征及斷流的影響因素分析(云南省地質(zhì)環(huán)境監(jiān)測(cè)院)[C].巖溶水文地質(zhì)研討會(huì)暨第十八屆全國(guó)洞穴學(xué)術(shù)會(huì)議,2012.

    [2] 周友妹.麗江古城黑龍?zhí)度當(dāng)嗔髟蚍治鯷J].人民珠江.2015(01):54-58.

    [3] 和菊芳,方金鑫.麗江黑龍?zhí)度畡?dòng)態(tài)變化與降水關(guān)系初探[J].中國(guó)農(nóng)村水利水電. 2019(08): 26-27.

    [4] 譚乃元,張麗芳.麗江黑龍?zhí)度當(dāng)嗔鞒跆絒J].城市建設(shè), 2012(17).

    [5] 李豫馨,許模,高偉,漆繼紅. 基于時(shí)間序列方法預(yù)測(cè)麗江黑龍?zhí)度蛄髁縖J].人民珠江, 2016(3):6-9.

    [6] 李豫馨.基于時(shí)間序列分析的麗江黑龍?zhí)度騽?dòng)態(tài)研究[J].成都理工大學(xué),2016.

    [7] 曾成,楊睿,楊明明,胡君春,武貴華,樊宇紅 .麗江市黑龍?zhí)度簲嗔鞯娜斯ど窠?jīng)網(wǎng)絡(luò)模擬[J].中國(guó)巖溶, 2013(04) P391-397.

    [8]王小川,史峰,郁磊,李洋編著. MATLAP神經(jīng)網(wǎng)絡(luò)43個(gè)案例分析[Z]. MATLAP 技術(shù)論壇.

    [9] 陳明編著.MATLAP神經(jīng)網(wǎng)絡(luò)原理與實(shí)例精講[Z].

    猜你喜歡
    人工神經(jīng)網(wǎng)絡(luò)
    基于人工神經(jīng)網(wǎng)絡(luò)的故障診斷專利淺析
    人工神經(jīng)網(wǎng)絡(luò)發(fā)展歷史與訓(xùn)練算法概述
    基于人工神經(jīng)網(wǎng)絡(luò)的優(yōu)化配置研究
    基于人工神經(jīng)網(wǎng)絡(luò)的地震電場(chǎng)衛(wèi)星數(shù)據(jù)異常識(shí)別
    ANN、ANFIS和AR模型在日徑流時(shí)間序列預(yù)測(cè)中的應(yīng)用比較
    人工神經(jīng)網(wǎng)絡(luò)在旅游人數(shù)預(yù)測(cè)中的應(yīng)用研究
    中國(guó)房地產(chǎn)上市公司治理狀況預(yù)測(cè)
    淺析人工神經(jīng)網(wǎng)絡(luò)在汽車輪軸松曠間隙檢測(cè)中的應(yīng)用
    基于改進(jìn)隱馬爾科夫模型的畜禽全基因組關(guān)聯(lián)分析中的多重檢驗(yàn)方法
    基于人工神經(jīng)網(wǎng)絡(luò)的經(jīng)濟(jì)預(yù)測(cè)模型
    久久香蕉精品热| 精品午夜福利视频在线观看一区| 欧美日韩福利视频一区二区| 禁无遮挡网站| 日韩欧美国产一区二区入口| 日日爽夜夜爽网站| 首页视频小说图片口味搜索| 在线观看午夜福利视频| 制服丝袜大香蕉在线| 国产精品综合久久久久久久免费 | 久久久久精品国产欧美久久久| 丝袜美腿诱惑在线| 欧美激情 高清一区二区三区| 国产精品一区二区在线不卡| 99国产精品一区二区蜜桃av| 亚洲专区中文字幕在线| 国内精品久久久久久久电影| 91av网站免费观看| 国产一区二区激情短视频| 免费在线观看视频国产中文字幕亚洲| 后天国语完整版免费观看| 国产精品日韩av在线免费观看 | 夜夜夜夜夜久久久久| 18美女黄网站色大片免费观看| 亚洲av成人不卡在线观看播放网| 999久久久精品免费观看国产| 亚洲专区国产一区二区| 亚洲自偷自拍图片 自拍| 老鸭窝网址在线观看| 国产野战对白在线观看| 亚洲va日本ⅴa欧美va伊人久久| 精品卡一卡二卡四卡免费| 久久 成人 亚洲| 国产精华一区二区三区| 国产亚洲精品久久久久5区| 女生性感内裤真人,穿戴方法视频| 777久久人妻少妇嫩草av网站| 精品国产乱子伦一区二区三区| 熟女少妇亚洲综合色aaa.| www.熟女人妻精品国产| 精品欧美一区二区三区在线| 大码成人一级视频| 亚洲成国产人片在线观看| 宅男免费午夜| 好男人电影高清在线观看| 少妇 在线观看| 精品国产超薄肉色丝袜足j| 精品福利观看| 高潮久久久久久久久久久不卡| 国产成人av教育| 中亚洲国语对白在线视频| 视频区欧美日本亚洲| 久久欧美精品欧美久久欧美| 亚洲成人久久性| 欧美精品啪啪一区二区三区| www.精华液| 欧美乱妇无乱码| 国产亚洲精品久久久久5区| 亚洲av美国av| 亚洲成人精品中文字幕电影| 精品国产超薄肉色丝袜足j| 国产亚洲精品一区二区www| 中文字幕久久专区| 在线观看午夜福利视频| 久久人妻福利社区极品人妻图片| 一个人免费在线观看的高清视频| av欧美777| 一级黄色大片毛片| 免费在线观看日本一区| 人人澡人人妻人| 啦啦啦韩国在线观看视频| 久久 成人 亚洲| 亚洲天堂国产精品一区在线| 亚洲欧美激情综合另类| x7x7x7水蜜桃| 十分钟在线观看高清视频www| 午夜福利视频1000在线观看 | 夜夜爽天天搞| 99精品欧美一区二区三区四区| 88av欧美| 免费高清视频大片| 女生性感内裤真人,穿戴方法视频| 国产午夜福利久久久久久| 9191精品国产免费久久| 天堂影院成人在线观看| a级毛片在线看网站| 视频区欧美日本亚洲| 亚洲成av片中文字幕在线观看| 亚洲少妇的诱惑av| www.自偷自拍.com| 日本五十路高清| 搡老岳熟女国产| 久热爱精品视频在线9| 首页视频小说图片口味搜索| tocl精华| 日韩大码丰满熟妇| 丝袜人妻中文字幕| 制服诱惑二区| 大陆偷拍与自拍| 日韩精品免费视频一区二区三区| 女人被躁到高潮嗷嗷叫费观| 侵犯人妻中文字幕一二三四区| 亚洲最大成人中文| 亚洲男人的天堂狠狠| 国产精品秋霞免费鲁丝片| 久久中文字幕人妻熟女| 欧美午夜高清在线| 久久国产精品男人的天堂亚洲| 好男人电影高清在线观看| 久久国产亚洲av麻豆专区| 我的亚洲天堂| 啦啦啦韩国在线观看视频| 欧美乱码精品一区二区三区| svipshipincom国产片| 在线观看日韩欧美| 在线观看免费日韩欧美大片| 在线观看免费午夜福利视频| 热99re8久久精品国产| 国产97色在线日韩免费| 免费看十八禁软件| 男女下面插进去视频免费观看| 18禁裸乳无遮挡免费网站照片 | 国产亚洲av嫩草精品影院| 岛国视频午夜一区免费看| 欧美一区二区精品小视频在线| 精品一区二区三区视频在线观看免费| 久热爱精品视频在线9| 亚洲色图av天堂| 波多野结衣巨乳人妻| x7x7x7水蜜桃| 男人的好看免费观看在线视频 | 欧美黄色片欧美黄色片| 97人妻精品一区二区三区麻豆 | av片东京热男人的天堂| 久久午夜亚洲精品久久| 变态另类丝袜制服| 亚洲伊人色综图| 成人特级黄色片久久久久久久| www日本在线高清视频| 亚洲av电影不卡..在线观看| 午夜免费激情av| 国产成人精品无人区| 久久精品国产亚洲av高清一级| 免费女性裸体啪啪无遮挡网站| 极品教师在线免费播放| 欧美日韩亚洲综合一区二区三区_| 男人舔女人下体高潮全视频| 午夜福利一区二区在线看| 婷婷精品国产亚洲av在线| 在线观看免费视频日本深夜| 一本综合久久免费| 中文字幕av电影在线播放| 黄色a级毛片大全视频| 精品久久久久久久久久免费视频| 99久久精品国产亚洲精品| 国产精品美女特级片免费视频播放器 | 日韩一卡2卡3卡4卡2021年| 免费看十八禁软件| 黄片大片在线免费观看| 欧美色欧美亚洲另类二区 | avwww免费| 午夜福利欧美成人| 国产91精品成人一区二区三区| 久久久精品欧美日韩精品| 久久久久九九精品影院| 首页视频小说图片口味搜索| 日本在线视频免费播放| 欧美人与性动交α欧美精品济南到| 亚洲人成电影观看| 在线观看免费日韩欧美大片| 黄色片一级片一级黄色片| 黄色丝袜av网址大全| 国产精品亚洲美女久久久| 国产精品亚洲一级av第二区| 免费在线观看视频国产中文字幕亚洲| 国产精品99久久99久久久不卡| 亚洲全国av大片| 欧美 亚洲 国产 日韩一| 日韩国内少妇激情av| 亚洲色图av天堂| 精品一区二区三区视频在线观看免费| 一区二区三区高清视频在线| 亚洲av成人不卡在线观看播放网| 欧美久久黑人一区二区| 国产片内射在线| 国内精品久久久久精免费| 亚洲第一电影网av| 午夜精品久久久久久毛片777| 国产精品精品国产色婷婷| 亚洲成人精品中文字幕电影| 乱人伦中国视频| 国产伦一二天堂av在线观看| 亚洲欧美精品综合久久99| 精品熟女少妇八av免费久了| 天天添夜夜摸| 19禁男女啪啪无遮挡网站| 亚洲成人国产一区在线观看| 亚洲国产精品sss在线观看| 国产亚洲av高清不卡| 久久香蕉精品热| 在线国产一区二区在线| 久久久国产成人精品二区| 中文字幕人妻熟女乱码| 亚洲精品久久国产高清桃花| 中文字幕另类日韩欧美亚洲嫩草| 久久久水蜜桃国产精品网| 黑丝袜美女国产一区| 免费在线观看完整版高清| 亚洲片人在线观看| 久久精品人人爽人人爽视色| 欧美激情久久久久久爽电影 | 国产99久久九九免费精品| 亚洲精品国产色婷婷电影| 精品第一国产精品| 精品第一国产精品| 国产97色在线日韩免费| 女人被狂操c到高潮| 又大又爽又粗| 国产亚洲欧美在线一区二区| 亚洲成人精品中文字幕电影| 波多野结衣高清无吗| 999久久久精品免费观看国产| 黄网站色视频无遮挡免费观看| 午夜福利高清视频| 大型黄色视频在线免费观看| 久热爱精品视频在线9| 精品国产乱子伦一区二区三区| 国产高清videossex| 91大片在线观看| 制服人妻中文乱码| 国产成人一区二区三区免费视频网站| 99在线视频只有这里精品首页| 亚洲精品久久国产高清桃花| av在线天堂中文字幕| 大陆偷拍与自拍| 大陆偷拍与自拍| 久久久久久大精品| 国产一区二区在线av高清观看| 久久久久久大精品| 很黄的视频免费| 后天国语完整版免费观看| 黄色视频,在线免费观看| 欧美黑人欧美精品刺激| 国产精品久久久久久精品电影 | 母亲3免费完整高清在线观看| 99国产综合亚洲精品| 欧美日本中文国产一区发布| 国产精品,欧美在线| 91麻豆精品激情在线观看国产| 成人亚洲精品一区在线观看| 久久中文字幕一级| 亚洲黑人精品在线| 中文字幕另类日韩欧美亚洲嫩草| 精品一品国产午夜福利视频| 男人舔女人的私密视频| 在线观看免费视频日本深夜| 亚洲专区字幕在线| 欧美日本视频| tocl精华| 99精品久久久久人妻精品| 亚洲欧美精品综合久久99| 国产人伦9x9x在线观看| av天堂久久9| 日本三级黄在线观看| 午夜精品久久久久久毛片777| 欧美日韩乱码在线| 中文字幕另类日韩欧美亚洲嫩草| 免费搜索国产男女视频| 国产极品粉嫩免费观看在线| 国内精品久久久久久久电影| 很黄的视频免费| 熟女少妇亚洲综合色aaa.| 黄色a级毛片大全视频| 最新美女视频免费是黄的| 国产精品永久免费网站| 亚洲avbb在线观看| 久久精品国产99精品国产亚洲性色 | 午夜福利免费观看在线| 大型av网站在线播放| 国产欧美日韩综合在线一区二区| 校园春色视频在线观看| 国产aⅴ精品一区二区三区波| 欧美一级毛片孕妇| 91麻豆av在线| 视频在线观看一区二区三区| 丝袜人妻中文字幕| 88av欧美| 日韩欧美一区视频在线观看| 波多野结衣高清无吗| 国产成人系列免费观看| 午夜福利欧美成人| 国产成人系列免费观看| 国产欧美日韩一区二区三| 女生性感内裤真人,穿戴方法视频| 国产一区二区三区在线臀色熟女| 男女床上黄色一级片免费看| 变态另类丝袜制服| 国产av一区二区精品久久| 色综合亚洲欧美另类图片| 欧美午夜高清在线| 国产在线观看jvid| 久久精品国产亚洲av香蕉五月| 宅男免费午夜| 亚洲熟妇熟女久久| 国产一区在线观看成人免费| av视频免费观看在线观看| 美女高潮喷水抽搐中文字幕| 99热只有精品国产| 欧美激情极品国产一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 色播在线永久视频| 欧美日韩精品网址| 日韩精品中文字幕看吧| 午夜激情av网站| 欧美+亚洲+日韩+国产| 久久精品国产清高在天天线| www.999成人在线观看| 国产成年人精品一区二区| 国产av一区二区精品久久| 1024香蕉在线观看| 国产成年人精品一区二区| 狠狠狠狠99中文字幕| 久久久精品国产亚洲av高清涩受| 中出人妻视频一区二区| 91老司机精品| 国产伦一二天堂av在线观看| 999精品在线视频| 久久亚洲精品不卡| 黑丝袜美女国产一区| 男女之事视频高清在线观看| 日本免费a在线| 国产成人精品久久二区二区91| 男人的好看免费观看在线视频 | 精品欧美国产一区二区三| 国产精品久久视频播放| 1024视频免费在线观看| 日韩一卡2卡3卡4卡2021年| 国产国语露脸激情在线看| aaaaa片日本免费| 久久久久精品国产欧美久久久| 中文字幕av电影在线播放| 激情在线观看视频在线高清| 日本免费一区二区三区高清不卡 | 欧美在线一区亚洲| 人人妻人人澡人人看| x7x7x7水蜜桃| xxx96com| 老熟妇乱子伦视频在线观看| 国产亚洲av高清不卡| 亚洲精品一卡2卡三卡4卡5卡| 老司机靠b影院| 国产在线观看jvid| 国产又爽黄色视频| 精品人妻1区二区| 黑人巨大精品欧美一区二区mp4| 久久香蕉精品热| 国产精品乱码一区二三区的特点 | 国产在线观看jvid| 欧美午夜高清在线| 亚洲无线在线观看| 国产精品自产拍在线观看55亚洲| 高清毛片免费观看视频网站| 一本久久中文字幕| 性少妇av在线| 多毛熟女@视频| 国产成人啪精品午夜网站| 欧美最黄视频在线播放免费| 激情在线观看视频在线高清| 亚洲精品国产区一区二| 乱人伦中国视频| av在线播放免费不卡| 长腿黑丝高跟| 91av网站免费观看| 十八禁人妻一区二区| 最新在线观看一区二区三区| 欧美国产日韩亚洲一区| 国产1区2区3区精品| 亚洲av五月六月丁香网| 亚洲三区欧美一区| 久久人人精品亚洲av| 国产真人三级小视频在线观看| 精品人妻在线不人妻| 无遮挡黄片免费观看| 国产av一区二区精品久久| 美女高潮到喷水免费观看| 男女之事视频高清在线观看| 精品国产美女av久久久久小说| 精品国产美女av久久久久小说| 免费高清视频大片| 极品教师在线免费播放| 国产成人一区二区三区免费视频网站| 国产成人啪精品午夜网站| 久久精品亚洲精品国产色婷小说| 精品福利观看| 久久久久亚洲av毛片大全| 日韩欧美一区二区三区在线观看| videosex国产| 欧美日韩一级在线毛片| 精品久久蜜臀av无| 韩国av一区二区三区四区| 12—13女人毛片做爰片一| 一区二区日韩欧美中文字幕| 国产精品九九99| 成人手机av| xxx96com| 亚洲电影在线观看av| 亚洲第一av免费看| 男女床上黄色一级片免费看| 少妇裸体淫交视频免费看高清 | 国产精品乱码一区二三区的特点 | 一二三四社区在线视频社区8| 久久久国产精品麻豆| 国产熟女午夜一区二区三区| 夜夜躁狠狠躁天天躁| 亚洲欧美精品综合一区二区三区| 亚洲欧美日韩另类电影网站| 午夜影院日韩av| 国产成人精品久久二区二区91| 97超级碰碰碰精品色视频在线观看| 色av中文字幕| 久久精品亚洲熟妇少妇任你| 亚洲美女黄片视频| 久久精品人人爽人人爽视色| bbb黄色大片| 国产精品久久电影中文字幕| 久久久久久久久中文| 精品国产亚洲在线| 色在线成人网| 亚洲狠狠婷婷综合久久图片| 18禁黄网站禁片午夜丰满| 欧美丝袜亚洲另类 | 亚洲五月色婷婷综合| 乱人伦中国视频| 久久精品影院6| 日韩有码中文字幕| 午夜久久久久精精品| 欧美 亚洲 国产 日韩一| 久久久久亚洲av毛片大全| 亚洲熟妇中文字幕五十中出| 亚洲aⅴ乱码一区二区在线播放 | 99精品欧美一区二区三区四区| 叶爱在线成人免费视频播放| 欧美日韩一级在线毛片| 国产片内射在线| 超碰成人久久| 人人妻人人澡欧美一区二区 | 欧美在线一区亚洲| 极品人妻少妇av视频| 久久久久精品国产欧美久久久| av免费在线观看网站| 精品久久久久久,| а√天堂www在线а√下载| 精品无人区乱码1区二区| 黄片播放在线免费| 亚洲av成人一区二区三| 欧美精品啪啪一区二区三区| 侵犯人妻中文字幕一二三四区| 国产主播在线观看一区二区| 夜夜爽天天搞| 日韩欧美国产在线观看| 手机成人av网站| 亚洲国产精品999在线| 久久人妻熟女aⅴ| 在线av久久热| 国产精品亚洲美女久久久| 老汉色∧v一级毛片| 日韩欧美三级三区| 一边摸一边抽搐一进一小说| 美女 人体艺术 gogo| 国产99久久九九免费精品| av免费在线观看网站| 一区福利在线观看| 中亚洲国语对白在线视频| 久久午夜综合久久蜜桃| 成人av一区二区三区在线看| 19禁男女啪啪无遮挡网站| 三级毛片av免费| 日日爽夜夜爽网站| 久久精品91蜜桃| 在线永久观看黄色视频| 69av精品久久久久久| 亚洲欧美精品综合久久99| 午夜成年电影在线免费观看| 国产精品国产高清国产av| 无遮挡黄片免费观看| 露出奶头的视频| 老司机午夜福利在线观看视频| 一本大道久久a久久精品| 在线观看免费视频网站a站| 一本久久中文字幕| 九色国产91popny在线| 亚洲国产精品合色在线| 亚洲全国av大片| 琪琪午夜伦伦电影理论片6080| 波多野结衣一区麻豆| 国产精品国产高清国产av| 亚洲国产精品999在线| 国产亚洲精品久久久久5区| 免费在线观看黄色视频的| 一级毛片精品| 一区二区三区高清视频在线| 女生性感内裤真人,穿戴方法视频| 69精品国产乱码久久久| 最近最新中文字幕大全电影3 | 婷婷丁香在线五月| 99国产精品99久久久久| 国产精品久久久久久亚洲av鲁大| 美女高潮喷水抽搐中文字幕| 国产色视频综合| 久久国产精品人妻蜜桃| 动漫黄色视频在线观看| 少妇裸体淫交视频免费看高清 | 日日夜夜操网爽| 男人舔女人下体高潮全视频| 在线观看日韩欧美| 老熟妇乱子伦视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 此物有八面人人有两片| av在线天堂中文字幕| 一进一出抽搐动态| 丝袜美足系列| 老司机靠b影院| www.999成人在线观看| 一区二区三区精品91| 亚洲精品av麻豆狂野| 黑人巨大精品欧美一区二区蜜桃| 国产成人系列免费观看| 一级片免费观看大全| 成人18禁高潮啪啪吃奶动态图| 18禁观看日本| 久久精品国产亚洲av高清一级| 可以免费在线观看a视频的电影网站| 可以在线观看毛片的网站| 亚洲成人免费电影在线观看| 男女床上黄色一级片免费看| aaaaa片日本免费| 日韩国内少妇激情av| 99精品欧美一区二区三区四区| 久久久精品国产亚洲av高清涩受| 久久精品亚洲精品国产色婷小说| 欧美亚洲日本最大视频资源| 极品教师在线免费播放| 丰满人妻熟妇乱又伦精品不卡| 麻豆成人av在线观看| 人人妻人人爽人人添夜夜欢视频| 中文字幕人妻熟女乱码| 午夜福利影视在线免费观看| 男女下面进入的视频免费午夜 | 中文字幕精品免费在线观看视频| 国产不卡一卡二| 久久性视频一级片| 国产精品亚洲美女久久久| 精品一品国产午夜福利视频| 18美女黄网站色大片免费观看| 一区二区三区国产精品乱码| 欧美黑人精品巨大| 久久性视频一级片| 色播在线永久视频| 一区福利在线观看| 亚洲五月婷婷丁香| 成人欧美大片| 天天躁狠狠躁夜夜躁狠狠躁| 久久亚洲真实| 一级毛片精品| 国产成人系列免费观看| 大型av网站在线播放| 国内精品久久久久精免费| 99国产极品粉嫩在线观看| 夜夜爽天天搞| 丁香六月欧美| 精品国产亚洲在线| 国产黄a三级三级三级人| 午夜精品久久久久久毛片777| 最新美女视频免费是黄的| 久久人人97超碰香蕉20202| 免费观看人在逋| 人人澡人人妻人| 国产成人系列免费观看| 中国美女看黄片| 欧美久久黑人一区二区| 一区二区三区高清视频在线| 狠狠狠狠99中文字幕| 淫妇啪啪啪对白视频| 精品久久久久久,| 在线观看午夜福利视频| 黄色毛片三级朝国网站| 亚洲第一av免费看| 久久香蕉激情| 日韩国内少妇激情av| 色综合亚洲欧美另类图片| 极品人妻少妇av视频| 久热爱精品视频在线9| 天堂√8在线中文| 91麻豆av在线| 亚洲国产高清在线一区二区三 | 午夜久久久在线观看| 婷婷精品国产亚洲av在线| 亚洲人成电影免费在线| 免费无遮挡裸体视频| 黄色视频,在线免费观看| 99久久99久久久精品蜜桃| 亚洲,欧美精品.| 熟女少妇亚洲综合色aaa.| 国产精品香港三级国产av潘金莲| 久9热在线精品视频| 成人永久免费在线观看视频| 久久精品国产清高在天天线| 亚洲一区高清亚洲精品| 男女之事视频高清在线观看| 国产亚洲精品一区二区www| 国内毛片毛片毛片毛片毛片| 国产成人av激情在线播放| 老汉色∧v一级毛片| 满18在线观看网站| 亚洲中文日韩欧美视频|