亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于時間序列的高爐水溫差多維度模糊綜合評判

        2019-11-18 05:44:54崔桂梅張勝男張勇馬祥
        中國測試 2019年9期
        關(guān)鍵詞:模糊綜合評判多維度

        崔桂梅 張勝男 張勇 馬祥

        摘要:為克服爐缸局部熱狀態(tài)評判的局限性,全面考慮高爐各冷卻壁段水溫差對熱狀態(tài)的影響,該文建立基于時間序列的高爐水溫差多維度模糊綜合評判模型。首先,選取爐身、爐腰、爐腹、爐缸各冷卻壁段水溫差時間序列,用相關(guān)系數(shù)法確定上部各段水溫差對下部水溫差影響的滯后時間,并篩選相關(guān)系數(shù)大于閾值的參數(shù)為主要指標參數(shù);其次,用概率統(tǒng)計法劃分指標參數(shù)最優(yōu)區(qū)間,計算其相對劣化度值,選擇分段嶺形分布函數(shù)計算各指標參數(shù)隸屬度值;最后,對各層指標參數(shù)自下而上進行模糊合成,結(jié)合指標參數(shù)權(quán)重和最大隸屬度原則綜合獲得最終評判等級。經(jīng)實際數(shù)據(jù)仿真分析,模型評判命中率為81.7%,能實時地評判高爐實際水溫差狀態(tài)。

        關(guān)鍵詞:多維度;時間序列;相關(guān)系數(shù)法;概率統(tǒng)計分析;相對劣化度;模糊綜合評判

        中圖分類號:TP391 文獻標志碼:A 文章編號:1674-5124(2019)09-0013-08

        收稿日期:2019-01-18;收到修改稿日期:2019-03-27

        基金項目:國家自然科學基金(61763039)

        作者簡介:崔桂梅(1963-),女,河北保定市人,教授,博士,研究方向為復(fù)雜過程系統(tǒng)的建模及運行優(yōu)化控制研究。

        0 引言

        高爐系統(tǒng)具有非線性、時滯性及強耦合性。冶煉過程實質(zhì)上是高爐煉鐵系統(tǒng)與外界環(huán)境之間頻繁的物質(zhì)、能量、信息交換和轉(zhuǎn)化的過程。系統(tǒng)中的能量以焦炭、煤粉等物質(zhì)為載體進入高爐,經(jīng)過爐內(nèi)的復(fù)雜反應(yīng)再以鐵水、高爐煤氣等物質(zhì)為載體從高爐輸出。圖1為某高爐的整體結(jié)構(gòu)圖,有效容積為2500m3,爐身下部是軟融帶的起始位置[1],該區(qū)域主要受上升煤氣流的沖刷、化學侵蝕和熱震剝落;爐腹是軟融帶的根部,主要靠渣皮工作,長期經(jīng)歷著爐渣、鐵水的沖刷和高溫煤氣流侵蝕[2]。高爐各冷卻部位的水溫差,從爐底至高爐上部,水溫差逐漸增大,從爐腹開始,水溫差增大幅度加劇,高爐冷卻壁段水溫差的跳變可反映爐內(nèi)爐溫的波動性[3]爐內(nèi)渣皮的穩(wěn)定性、局部煤氣流強度及軟熔帶的相對位置等。根據(jù)高爐冷卻強度和冶煉強度的適應(yīng)性,將高爐立體水溫差作為表征高爐立體熱狀態(tài)的重要指標,高爐立體水溫差即對空間立體高爐冷卻壁段水溫多維度分析[4]。高爐水溫差的多維度評判為預(yù)判爐況奠定基礎(chǔ),從而提高生鐵質(zhì)量和實現(xiàn)節(jié)能降耗。

        近年來,國內(nèi)外學者采用一系列智能算法對高爐熱狀態(tài)進行研究。王華秋等[5]通過自反饋RBF神經(jīng)網(wǎng)絡(luò)模型預(yù)測鐵水中硅的含量來預(yù)報高爐熱狀態(tài),但該模型忽略了水溫差對高爐熱狀態(tài)的影響。周朝營等[6]用改進的PSO算法訓練BP神經(jīng)網(wǎng)絡(luò),建立爐溫預(yù)測模型,結(jié)合爐缸局部水溫差判斷爐體熱狀態(tài)。宋小鵬等[7]提出用爐缸水溫差和流量監(jiān)測爐體熱流判斷爐體熱狀態(tài)。針對傳統(tǒng)高爐熱狀態(tài),尚無高爐立體水溫差綜合評價方面的研究,高爐水溫差多維度綜合評價模型克服爐缸局部水溫差考慮的觀念,依據(jù)高爐煉鐵過程中能量流的傳遞特性及高爐冷卻壁段水溫差間的滯后性影響,對前一時刻和當前時刻水溫差建立評判模型,采用下一時刻的采樣數(shù)據(jù)驗證模型的有效性。由于實際工藝限制,如圖1所示,對高爐本體的冷卻,選取爐缸(H1、H2)、爐腹(H3)、爐腰(H4)、爐身(H5、H6)等6段冷卻壁,每段冷卻壁周圍等間隔采取48個測量點,圖1右側(cè)圖為冷卻壁測量點橫向刨面簡圖。建模主要過程為:根據(jù)相關(guān)性分析法提取有效特征[8],選出的指標參數(shù)經(jīng)劣化度分析后進行歸一化;利用相關(guān)系數(shù)法和專家經(jīng)驗對各層指標分配權(quán)重;計算覆蓋相對劣化度區(qū)間的隸屬度值;最后通過模糊綜合評判模型自下而上對每個子系統(tǒng)的狀態(tài)進行多屬性綜合評價。

        1 數(shù)據(jù)預(yù)處理

        選取該高爐在2018年10月21日至11月21日期間不同工況下的工業(yè)現(xiàn)場生產(chǎn)數(shù)據(jù)585組,評價指標H1、H2、H3、H4、H5、H6當前時刻及歷史時刻水溫差。

        1.1 異常值檢測與修補

        高爐平穩(wěn)運行時,采集的過程參數(shù)應(yīng)在小范圍內(nèi)波動,不會出現(xiàn)突跳點[9],3σ原則是根據(jù)過程數(shù)據(jù)的波動情況選取分布在(μ-3σ,μ+3σ)中數(shù)值,在此原則下,突跳點數(shù)據(jù)與平均值的偏差大于3σ。其中,μ和σ分別為樣本數(shù)據(jù)的平均值和標準差。

        原始樣本集中檢測出的異常數(shù)據(jù)用牛頓插值法修補[10]。對修補后的整體數(shù)據(jù),計算48個點的平均水溫差得到每段冷卻壁的實際水溫差。

        1.2 相關(guān)系數(shù)分析法提取關(guān)鍵特征

        由于建立模型時采集數(shù)據(jù)較多,相關(guān)分析法是一種采取數(shù)學降維的思想,選擇相關(guān)性大的因素作為評價體系的參數(shù)。圖2為高爐立體水溫差指標參數(shù)相關(guān)系數(shù)矩陣可視化圖,橫、縱軸均為高爐冷卻壁段k-1~k時刻的水溫差,圖中不同顏色的方塊為相關(guān)系數(shù)的可視化,右側(cè)的條形顏色圖代表每個相關(guān)系數(shù)對應(yīng)的顏色,相關(guān)系數(shù)為正時,方塊兒顏色趨于紅色系,反之,方塊兒顏色趨于藍色系。相關(guān)系數(shù)的絕對值越大,方塊面積越大,反之越小。

        相關(guān)系數(shù)絕對值越接近于1,參數(shù)相關(guān)度越強,越接近于0,相關(guān)度越弱。結(jié)合專家經(jīng)驗和圖2分析結(jié)果,選取|g|≥0.2的因素作為相關(guān)因素[11],爐缸保留H1(k-1)、H1(k)、H2(k-1),H2(k);爐腹保留H1(k-1)、H1(k)、H2(k-1)、H2(k)、H3(k-1)、H3(k):爐腰保留H2(k-1)、H2(k)、H3(k-1)、H3(k)、H4(k-1)、H4(k);爐身保留H3(k-1)、H3(k)、H4(k-1)、H4(k)、H5(k-1)、H5(k)、H6(k-1)、H6(k)。

        1.3 參數(shù)等級區(qū)間確定

        概率統(tǒng)計法:高爐長期運行于穩(wěn)定狀態(tài),各過程參數(shù)近似服從以理想工作點為均值的正態(tài)分布,通過公式(2)可將非標準正態(tài)曲線的橫軸區(qū)間折合為標準正態(tài)分布曲線區(qū)間進行計算,由查表法選取正態(tài)分布曲線下,橫軸(μ±10.19σ)區(qū)間內(nèi)面積為15%的區(qū)間為最優(yōu)區(qū)間[12]。圖3為H1、H2、H3、H4、H5、H6的概率統(tǒng)計分析圖。

        z=X-μ/δ(2)其中μ和δ分別為正態(tài)分布曲線的均值和方差,X為樣本數(shù)據(jù),z為標準化后的正態(tài)分布數(shù)據(jù)。

        由于評價模型需在爐況穩(wěn)定的情況下建立,各過程參數(shù)應(yīng)在小范圍內(nèi)波動,根據(jù)圖中各個區(qū)間擬合程度得到各參數(shù)取值區(qū)間及最優(yōu)區(qū)間如表1所示。

        2 高爐水溫差多維度模糊綜合評判模型

        模糊綜合評判法是一種基于模糊數(shù)學的方法,對受多種因素影響的事件進行綜合評判。模糊綜合評判能更好地解決具有模糊、難以量化和不確定性對象的狀態(tài)評估問題,它具有清晰、系統(tǒng)性強的特點,主要涉及4個要素:評語集L、指標體系集X、隸屬度矩陣R和模糊權(quán)向量W,本文中將高爐冷卻壁水溫差多維度評價狀態(tài)分為4個等級[13],記為L=[理想,良好,一般,預(yù)警]。根據(jù)影響因素的層次性,模糊綜合評判體系可分為一級和多級,對高爐立體水溫差采用兩級模糊綜合評判體系[14]。高爐立體水溫差多維度綜合評判模型建立的流程見圖4。

        2.1 構(gòu)建高爐水溫差多維度綜合評價指標體系

        由圖2的相關(guān)分析結(jié)果構(gòu)建高爐冷卻壁水溫差多維度綜合評價體系,一級指標集為X=[X1,X2,X3,X4],二級指標集分別為:

        X1指標參數(shù):H1(1)(k-1)、H1(1)(k)、H2(1)(k-1)、H2(1)(k)。

        X2指標參數(shù):H1(2)(k-1)、H1(2)(k)、H2(2)(k-1)、H2(2)(k)、H3(2)(k-1)、H3(2)(k)。

        X3指標參數(shù):H2(2)(k-1)、H2(2)(k)、H3(2)(k-1)、H3(2)(k)、H4(2)(k-1)、H4(2)(k)。

        X4指標參數(shù):H3(2)(k-1)、H3(2)(k)、H4(2)(k-1)、H4(2)(k)、H5(2)(k-1)、H6(2)(k-1)、H6(2)(k)。

        2.2 確定模糊權(quán)向量

        由于各因素對評價對象作用程度不同,應(yīng)設(shè)法給予不同的權(quán)重,確定評價指標體系權(quán)重的方法可分為兩類:主觀賦權(quán)和客觀賦權(quán)。主觀賦權(quán)主要依據(jù)專家經(jīng)驗賦權(quán)值,首先聘請有關(guān)專家,由每位專家根據(jù)總分s對每個考核指標的重要性打分,然后對每個考核指標的分數(shù)取平均值,作為最終權(quán)重??陀^賦權(quán)依據(jù)數(shù)據(jù)計算賦權(quán)值,如相關(guān)系數(shù)比值法和熵權(quán)法[15]。高爐水溫差多維度模糊綜合評價體系利用專家經(jīng)驗對一級指標賦權(quán)值,設(shè)權(quán)重矩陣為W,相關(guān)系數(shù)法對二級指標進行客觀賦權(quán),權(quán)重矩陣為w。

        專家經(jīng)驗對一級指標X賦為:

        W=[W1,W2,W3,W4]

        式中:aj——第j個指標參數(shù)權(quán)重;

        n——專家總數(shù);

        m——指標參數(shù)個數(shù);

        spj——第p個專家對第j個指標打分。

        相關(guān)系數(shù)比重計算二級指標參數(shù)權(quán)重為

        w=[wi1,wi2,…,wij]

        式中:cij——第i個子系統(tǒng)第j個指標的相關(guān)系數(shù);

        wij——第i個子系統(tǒng)第j個指標的權(quán)值;

        m——子系統(tǒng)X1、X2、X3、X4參數(shù)指標的個數(shù)。

        2.3 評估指標相對劣化度分析

        相對劣化度是指與故障狀態(tài)下的劣化程度相比,其當前實際狀態(tài)的劣化程度,對比各指標因素,需對所有數(shù)據(jù)進行相對劣化度處理消除量綱的影響,相對劣化度的取值范圍為[0,1],取值越小表示指標的狀態(tài)越好,相對劣化度計算方法有2類,即越小越優(yōu)型和中間型[16]。高爐立體水溫差評價體系中,過程參數(shù)工作在區(qū)間內(nèi),故選擇中間型函數(shù)計算劣化度值。中間型劣化度函數(shù)的總體參數(shù)需要估計變量的最大值xmax、最小值xmin和最佳范圍[xa,xb],最佳取值范圍可由概率統(tǒng)計分析法確定,其劣化度函數(shù)為

        2.4 因子隸屬度矩陣計算

        隸屬度是將評估指標的劣化度轉(zhuǎn)化成評估等級隸屬度的重要工具,各指標參數(shù)經(jīng)相對劣化度分析得到歸一化后的值,模糊綜合評價的關(guān)鍵是合理選取覆蓋相對劣化度區(qū)間的隸屬度函數(shù)。將高爐冷卻壁水溫差多維度評價的4個評價狀態(tài)記為L=[理想,良好,一般,預(yù)警],計算各等級劣化范圍內(nèi)數(shù)據(jù)的隸屬度值,區(qū)間內(nèi)取值偏小的數(shù)據(jù),隸屬度函數(shù)采用降嶺形分布,中間取值采用中間型嶺形分布函數(shù),區(qū)間內(nèi)取值偏大的數(shù)據(jù)采用升嶺形分布函數(shù),隸屬度的交叉重疊率控制在0.2~0.6[17]。根據(jù)隸屬度隨相對劣化度變化的原則,得到各評價等級的隸屬度函數(shù)為

        dij為第i個子系統(tǒng)第j個指標參數(shù)的劣化度,rij1、rij2、rij3、rij4為第i個子系統(tǒng)第j個指標參數(shù)對4個評價等級的隸屬度函數(shù),構(gòu)建隸屬度矩陣

        rij=[rij1,rij2,rij3,rij4j

        4個評價等級在不同劣化度區(qū)間的隸屬度函數(shù)如圖5所示。

        2.5 等級評價原則

        對評價體系指標參數(shù)從下到上逐層進行模糊合成,計算隸屬度矩陣向量得到最終評價向量,模糊合成過程為其中,bi表示第i個子系統(tǒng)的評判等級,C是各類指標參數(shù)的綜合評判矩陣,rji為第i個子系統(tǒng)第j個指標參數(shù)的隸屬度值。由以上模糊合成過程可計算最終評判向量,根據(jù)最大隸屬度原則有

        bi0=max{bi},1≤i≤4(12)

        等級評判原則:

        經(jīng)過多層權(quán)值的傳遞,劣化現(xiàn)象在最終層可能被掩蓋,最終評價原則不采用“最大隸屬度原則”,針對最終評價向量分情況討論:

        1)當最大隸屬度遠大于其他等級項的隸屬度,且權(quán)重最大的等級項隸屬度極低(<0.3)時,采用“最大隸屬度原則”進行評判。

        2)當最大隸屬度遠大于其他等級項的隸屬度,權(quán)重最大的等級項隸屬度也較高(≥0.3)時,則表明權(quán)重最大的等級項為評價對象的狀態(tài)。

        3)當有2個以上的等級項隸屬度與最大隸屬度差別較小時,需對其權(quán)重分析,權(quán)重較大的這些等級項采取“取隸屬度>0的最低等級項”的原則。

        4)當所有的等級項隸屬度較相近時,評價原則以權(quán)重為主導,取權(quán)重最大的等級項。

        3 實例分析

        選取該高爐10月21日至11月21日期間不同爐況下的現(xiàn)場生產(chǎn)數(shù)585組,以其中兩組數(shù)據(jù)為例進行詳細分析,對高爐立體水溫差進行兩級模糊綜合評判,驗證本文所提模型的有效性,詳細分析過程如下:

        1)利用相關(guān)系數(shù)及專家經(jīng)驗確定模糊權(quán)向量,由相關(guān)系數(shù)計算得到二級指標權(quán)重:

        由專家經(jīng)驗分配一級指標權(quán)重:

        W=[0.37,0.23,0.21,0.19]

        2)根據(jù)圖2選取的最優(yōu)區(qū)間,采用中間型劣化度函數(shù)計算劣化度值如表2所示。

        3)根據(jù)隸屬度函數(shù)計算二級指標的隸屬度矩陣,設(shè)第一組數(shù)據(jù)隸屬度分別為R11、R12、R13、R14,第二組數(shù)據(jù)隸屬度分別為R21、R22、R23、R24,詳細結(jié)果如下所示。由公式(11)計算一級指標隸屬度矩陣

        由公式(10)計算最終評價向量:

        C1=W1R1=[0.715 0.264 0 0]

        C2=W2R2=[0.270 0.043 0.680 0]

        綜合權(quán)重和最終評價矩陣分析:

        第一組數(shù)據(jù)評判結(jié)果為理想狀態(tài),第二組數(shù)據(jù)評判結(jié)果為一般狀態(tài),觀察該時刻高爐其他數(shù)據(jù)如表3所示。

        綜合表格數(shù)據(jù)分析:

        第一組生產(chǎn)數(shù)據(jù):由十字測溫法計得到的z/w值位于最優(yōu)區(qū)間內(nèi),此時高爐邊緣煤氣流分布合理,渣皮產(chǎn)生穩(wěn)定,軟化區(qū)間和熔滴區(qū)間的溫度均位于最優(yōu)區(qū)間內(nèi),軟熔帶位置合適,初步判斷該時刻爐況穩(wěn)定。

        第二組數(shù)據(jù):z/w值減小,邊緣煤氣流旺盛,渣皮脫落,冷卻壁熱面溫度升高,軟化區(qū)間和熔滴區(qū)間的溫度均低于最優(yōu)區(qū)間內(nèi)取值,軟熔帶位置偏高,間接還原區(qū)縮減,縱觀該時刻附近數(shù)據(jù),爐溫波動造成礦石軟熔帶根部位置頻繁變化,爐況稍有波動。

        在實際建模過程中,對選取的585組數(shù)據(jù)均做了仿真分析,圖6為585組數(shù)據(jù)實際仿真結(jié)果,模型的命中率為81.7%(允許誤差為0)。該模型命中率高,能實時地評價高爐實際水溫差狀態(tài)。

        4 結(jié)束語

        為克服傳統(tǒng)熱狀態(tài)爐缸局部水溫差評價,考慮冷卻壁段水溫差間的滯后性影響,本文提出一種基于時間序列的高爐水溫差多維度綜合評判模型。首先針對數(shù)據(jù)預(yù)處理后的高爐生產(chǎn)過程數(shù)據(jù)建立單因素評判矩陣,同時引入相關(guān)性分析法計算各指標參數(shù)的權(quán)重;其次結(jié)合單因素隸屬度評判矩陣和專家經(jīng)驗確定的模糊權(quán)向量進行二級評判;最后依據(jù)二級隸屬度評判矩陣和評價等級原則確定最終等級。高爐立體水溫差可間接反映高爐爐況的穩(wěn)定性,在此基礎(chǔ)上做高爐立體熱狀態(tài)評判,更能準確地預(yù)判爐況的穩(wěn)定性,實驗結(jié)果表明該研究具有良好的延展性,可指導爐長提前、準確地判斷爐況的穩(wěn)定性,對企業(yè)降低成本、節(jié)焦降耗及高爐優(yōu)化運行具有現(xiàn)實意義。

        參考文獻

        [1]楊貴軍,蔣朝輝,桂衛(wèi)華,等.基于熵權(quán)-可拓理論的高爐軟熔帶位置狀態(tài)模糊綜合評判方法[J].自動化學報,2015,41(1):75-83.

        [2]馬洪斌,張賀順.首鋼2號高爐銅冷卻壁使用的體會[J].煉鐵,2008,27(5):9-12.

        [3]杜國萍.包鋼6#高爐冷卻系統(tǒng)運行實踐[J].包鋼科技,2011,37(3):6-8.

        [4]余斌,張少偉,劉洋,等.首鋼京唐5500m-3高爐冷卻壁水溫差無線監(jiān)測系統(tǒng)[J].山東冶金,2018,40(1):60-62.

        [5]王華秋,廖曉峰,鄒航,等.自反饋RBF網(wǎng)絡(luò)在高爐熱狀態(tài)模型預(yù)測中的應(yīng)用[J].系統(tǒng)工程與電子技術(shù),2008(5):929-934.

        [6]周朝萱.基于粒子群的BP神經(jīng)網(wǎng)絡(luò)在高爐爐熱狀態(tài)預(yù)報中的應(yīng)用[A].Intelligent Information Technology ApplicationAssociation.Agricultural and Natural ResourcesEngineering(ANRE 2011 ABE V3)[C].Intelligent InformationTechnology Application Association:智能信息技術(shù)應(yīng)用學會,2011.

        [7]宋小鵬.高爐爐體熱狀態(tài)在線監(jiān)測系統(tǒng)[A].中國金屬學會.2012年全國煉鐵生產(chǎn)技術(shù)會議暨煉鐵學術(shù)年會文集(下)[C].中國金屬學會:中國金屬學會,2012.

        [8]崔桂梅,李靜,張勇,等.高爐鐵水溫度的多元時間序列建模和預(yù)測[J].鋼鐵研究學報,2014,26(4):33-37.

        [9]張良均,楊坦,肖鋼,等.數(shù)據(jù)分析與挖掘?qū)崙?zhàn)[M].北京:機械工業(yè)出版社,2016:47-53.

        [10]吳禮斌,李柏年,張孔生,等.MATLAB數(shù)據(jù)分析方法[M].北京:機械工業(yè)出版社,2017:39-57.

        [11]周英,舊金武,卞月青.大數(shù)據(jù)挖掘系統(tǒng)方法與實例分析[M].北京:機械工業(yè)出版社,2016:106-107.

        [12]楊立才,張朋飛,王德偉.基于主成分分析的代謝綜合征模糊綜合評價方法[J].生物醫(yī)學工程學雜志,2013,30(1):67-70,79.

        [13]陳龍,李濤,徐嘯峰.基于概率密度函數(shù)的時滯依賴故障檢測與診斷[J].計算機工程與應(yīng)用,2018,54(6):257-263.

        [14]代雪靜,田衛(wèi).水質(zhì)模糊評價模型中賦權(quán)方法的選擇[J].中國科學院研究生院學報,2011,28(2):169-176.

        [15]李晨曦,孫哲,蔣景英,等.近紅外光譜主成分分析與模糊聚類的典型地面目標物識別[J].光譜學與光譜分析,2017,37(11):3386-3390.

        [16]黃必清,何焱,王婷艷.基于模糊綜合評價的海上直驅(qū)風電機組運行狀態(tài)評估[J].清華大學學報(自然科學版),2015,55(5):543-549.

        [17]徐銘銘,曹文思,姚森,等.基于模糊層次分析法的配電網(wǎng)重復(fù)多發(fā)性停電風險評估[J].電力自動化設(shè)備2018,38(10):19-25,31.

        (編輯:商丹丹)

        猜你喜歡
        模糊綜合評判多維度
        空間角與距離的多維度解法
        “多維度評改”方法初探
        多維度市南
        商周刊(2017年7期)2017-08-22 03:36:22
        模糊綜合評判在優(yōu)化電機冷卻系統(tǒng)中的應(yīng)用
        模糊綜合評判在優(yōu)化電機冷卻系統(tǒng)中的應(yīng)用
        公路橋梁臺背回填質(zhì)量評價方法研究
        四維主體視閾下地方應(yīng)用型本科院校學生綜合素質(zhì)的模糊評價與應(yīng)用
        基于模糊層次分析法的公路橋梁施工安全風險評價研究
        中國市場(2016年41期)2016-11-28 05:25:17
        校企合作模式下的卓越軟件工程師培養(yǎng)質(zhì)量的研究
        will與be going to的多維度意義對比
        国产精品一级黄色大片| 国产精品久久久久久妇女6080| 久久久亚洲欧洲日产国产成人无码| 亚洲av综合色区久久精品天堂| 久久这里都是精品99| 国产亚洲精品美女久久久| 窝窝影院午夜看片| 亚洲色欲色欲欲www在线| 女女同女同一区二区三区| 国产精久久一区二区三区| 国产精品久久久av久久久| 中文亚洲成a人片在线观看| 国产视频一区二区三区观看 | 久久久无码一区二区三区| 无码中文字幕专区一二三| 特级黄色大片性久久久| 无码人妻一区二区三区兔费| 久久久久久av无码免费看大片| 色综合久久久久综合一本到桃花网| 国产精品一区久久综合| 免费无码av一区二区| 国产精品国产三级农村妇女| 婷婷开心五月综合基地| 97无码人妻一区二区三区蜜臀| 亚洲成人激情深爱影院在线| 人妻精品久久无码区| 午夜家庭影院| 国产精品一区二区久久毛片| 丰满少妇按摩被扣逼高潮| 亚洲av永久无码天堂网毛片| 国产一区二区三区韩国| 日韩精品一区二区亚洲观看av| 亚洲中文字幕久久精品无码a| 国产女人成人精品视频| 丝袜美腿亚洲综合玉足| 无码熟妇人妻av在线网站| 欧美老熟妇欲乱高清视频| baoyu网址国产最新| 少妇免费av一区二区三区久久| 成人爽a毛片在线视频| 在线av野外国语对白|