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

    基于前后向擬合的色譜重疊峰分峰方法

    2018-10-11 01:35:36高文清張俊良HuJunJack俞建成Tangkeqi
    質(zhì)譜學(xué)報(bào) 2018年5期
    關(guān)鍵詞:峰高單峰極大值

    高文清,張俊良,王 艷,劉 蔓,Hu JunJack, 俞建成,Tang keqi

    (1.寧波大學(xué)信息科學(xué)與工程學(xué)院,浙江 寧波 315211;2.中國人民解放軍陸軍總醫(yī)院,北京 100700; 3.中國科學(xué)院寧波材料技術(shù)與工程研究所,浙江 寧波 315201)

    色譜技術(shù)作為復(fù)雜組分分離和定量分析的手段之一,具有分析速度快、分離效率高、樣品用量少等優(yōu)勢(shì)[1]。應(yīng)用色譜技術(shù)分離復(fù)雜成分,在理想實(shí)驗(yàn)條件下,每個(gè)單峰對(duì)應(yīng)一個(gè)組分,但實(shí)際上,當(dāng)兩個(gè)或多種成分含有相似的結(jié)構(gòu)或性質(zhì)時(shí),很容易形成重疊峰,這將為后續(xù)的定性和定量分析帶來困難[2]。因此,如何有效分離重疊色譜峰,是目前亟需解決的重要問題之一[3-4]。實(shí)際應(yīng)用中,大多數(shù)重疊峰是雙峰重疊,所以本工作只研究2個(gè)峰重疊的分離方法。

    針對(duì)色譜重疊峰分離問題,主要的分離方法有基于幾何的、基于代數(shù)的和基于模式識(shí)別的分峰方法[5]。其中,基于幾何的分峰方法主要有垂線法和切線法,計(jì)算速度快、直觀,適用于實(shí)時(shí)在線分析,已被商業(yè)集成所應(yīng)用,但其分離精度會(huì)隨著色譜峰重疊程度的增加而降低,在嚴(yán)重重疊的情況下,這兩種方法的定量誤差都很大[6]。針對(duì)這一問題,葉國陽等[7]提出了交點(diǎn)垂線法和比例分配法,使定量精度得到了大幅提高,但該方法依賴于重疊峰的兩個(gè)極大值點(diǎn),而且在只有一個(gè)極大值點(diǎn),即肩峰的情況下,還需要通過做重疊峰拐點(diǎn)切線,求切線交點(diǎn)的方法求得另一個(gè)極大值點(diǎn)。另外,拖尾峰和前延峰在實(shí)際應(yīng)用中也很常見,該研究并未探討這兩種峰的適用度。嚴(yán)彩娟[8]基于重疊峰相似性原理提出了縱向迭代法,該方法從遠(yuǎn)離重疊區(qū)域的前后沿開始擬合,對(duì)色譜峰重疊區(qū)域進(jìn)行修正,并研究了拖尾重疊峰、前延重疊峰,但該方法要求重疊峰具有谷點(diǎn)?;诖鷶?shù)的分峰方法主要有基于高斯函數(shù)的譜峰擬合算法[9]、基于最小二乘法的譜峰擬合算法[10]、小波變換算法[11]等。這些算法能夠得到較好的計(jì)算結(jié)果和分峰效果,然而在分峰過程中,它們需要一定的參數(shù)估計(jì)與優(yōu)化、模型選擇等步驟,計(jì)算量大,消耗時(shí)間長,不適用于實(shí)時(shí)在線處理[12]?;谀J阶R(shí)別的分峰方法中,免疫算法分離效果較好,但它僅適用于已知成分的重疊峰分離,應(yīng)用受限[13]。

    1 前后向擬合分離重疊峰理論

    1.1 重疊峰數(shù)學(xué)模型

    色譜分離過程中,由于性質(zhì)或結(jié)構(gòu)上的微小差異,各組分與流動(dòng)相和固定相的相對(duì)運(yùn)動(dòng)產(chǎn)生的作用力以及兩相的分配系數(shù)不同,使得各組分在色譜柱上的保留時(shí)間不同,當(dāng)兩種物質(zhì)具有相似的性質(zhì)或結(jié)構(gòu)時(shí),從色譜柱流出的時(shí)間非常接近,產(chǎn)生重疊峰[7]。根據(jù)重疊峰產(chǎn)生的原理,假設(shè)重疊峰的單組分峰具有相似性,具體特征如下:1) 所有單峰的峰形相似;2) 單峰之間可以通過一個(gè)系數(shù)相互轉(zhuǎn)換;3) 每個(gè)單峰的半峰寬相同。

    為了便于理解,構(gòu)建了重疊峰的數(shù)學(xué)模型,示于圖1,其方程如下:

    (1)

    式中:OC代表重疊色譜峰;OA、OB分別代表構(gòu)成重疊峰的兩個(gè)單峰;+代表兩個(gè)單峰相加;×代表兩個(gè)單峰之間的尺度變換;K代表尺度變換因子,即兩個(gè)單峰的峰高比;Δt為平移跨度(Δt>0)。

    1.2 前后向擬合算法

    從圖1可以看出,單峰OA的前沿、單峰OB的后沿分別與重疊峰OC的前沿、后沿相差很小,即它們受到重疊的影響很小,因此,本工作從此處著手對(duì)重疊峰進(jìn)行分離。以tM、tN、tS、tT分別代表M、N、S、T四個(gè)點(diǎn)的時(shí)間坐標(biāo),hM、hN、hS、hT分別代表其在相應(yīng)時(shí)間上的峰強(qiáng)。

    注:M、N分別是重疊峰的兩個(gè)極大值點(diǎn); S、T分別是后向擬合、前向擬合的起始點(diǎn) 圖1 重疊峰的數(shù)學(xué)模型Fig.1 Mathematical model of overlapping peaks

    前后向擬合算法主要分為前向擬合和后向擬合兩步。

    定義S、T分別為后向擬合和前向擬合的起始點(diǎn),對(duì)于第一次前向擬合和后向擬合,具體步驟如下:

    泄露500 d后,上覆第四系松散孔隙含水層進(jìn)一步向南側(cè)遷移,污染羽狀物擴(kuò)散范圍有逐漸增加,而下伏巖溶含水層仍然未見受到明顯污染,未出現(xiàn)明顯污染羽狀物聚集及滲漏。

    1) 從T點(diǎn)出發(fā)開始前向擬合。將重疊峰的極大值M點(diǎn)當(dāng)作單峰OA暫時(shí)的極大值點(diǎn)(以m點(diǎn)表示),T點(diǎn)當(dāng)作單峰OB暫時(shí)的極大值點(diǎn)(以t點(diǎn)表示),以此修正單峰OB的峰高點(diǎn)。此時(shí)峰高比可以定義為:

    kft=hm/ht

    (2)

    2) 定義τ(τ>0)為時(shí)間跨度,其大小可以根據(jù)需求調(diào)節(jié)。沿著時(shí)間軸將t點(diǎn)向重疊區(qū)域推進(jìn)τ,現(xiàn)在時(shí)間表示為tt-τ。

    3) 根據(jù)式(1),擬合單峰OB在時(shí)間tt-τ到tt部分的峰形,把擬合出來的峰對(duì)應(yīng)于tt-τ點(diǎn)記作t1。此過程可由式(3)表示:

    (3)

    式中,OBft1(t)為前向擬合第一次迭代之后得到的單峰OB的后沿,此時(shí)OBft1(t)暫時(shí)的峰高點(diǎn)被記作ht1,同時(shí)峰高比可以表示為:

    kft1=hm/ht1

    (4)

    4) 依此類推,前向擬合經(jīng)過第(n+1)(n=0,1…)次迭代之后,tt-(n+1)τ到tt之間的峰形可以表示為:

    (5)

    kftn+1=hm/htn+1

    (6)

    5) 比較前向擬合過程前后兩次迭代得到的htn與htn+1值。如果htn≥htn+1,并且兩者之間的差值處于設(shè)定的范圍內(nèi),就完成了第一次前向擬合,可得到單峰OB暫時(shí)的后沿OBftn(t)和極大值htn。

    6) 結(jié)束第一次前向擬合后,從S點(diǎn)開始類似的后向擬合過程。將S點(diǎn)當(dāng)作單峰OA暫時(shí)的極大值點(diǎn)(以s點(diǎn)表示),保持單峰OB暫時(shí)的極大值htn不變,重復(fù)上述5個(gè)步驟,可以得到ts到ts+(n+1)τ之間的峰形,記作OAbsn+1(t):

    (7)

    kbsn+1=htn/hsn+1

    (8)

    如果hsn≥hsn+1,而且它們的差值落在一定的范圍內(nèi),第一次后向擬合過程即停止。此時(shí),可得到單峰OA暫時(shí)的前沿OAbsn(t)和極大值hsn。

    第一次前向擬合和第一次后向擬合的結(jié)果示于圖2。

    經(jīng)過第一次完整的前后向擬合過程后,根據(jù)上述6個(gè)步驟,進(jìn)行第R(R=2,3,…)次前后向擬合。此時(shí),可以得到單峰OA新的前沿OAbRsn和極大值hRsn,單峰OB新的后沿OBfRtn和極大值點(diǎn)hRtn。最后,通過色譜重疊峰相似性原理,單峰OA新的后沿OAback可以由單峰OB新的后沿OBfRtn得到,單峰OB新的前沿OBfront可以由單峰OA新的前沿OAbRsn得到。因此,完整的單峰可以表示為:

    OA(t)=OAbRsn(t)+OAback(t)

    OAback(t)=(hRsn/hRtn)×OBfRtn(t-ttn+tsn)

    (9)

    OB(t)=OBfRtn(t)+OBfront(t)

    OBfront(t)=(hRtn/hRsn)×OAbRsn(t+ttn-tsn)

    (10)

    為了更好地理解該算法,將經(jīng)過3次前后向擬合得到的分峰效果圖示于圖3。

    圖2 第一次前向擬合得到的單峰OB的 后沿OBftn和第一次后向擬合得到的 單峰OA的前沿OAbsnFig.2 Back edge OBftn of single peak OB obtained from the first forward fitting and the front edge OAbsn of single peak OA obtained from the first backward fitting

    圖3 經(jīng)過3次完整的前后向擬合得到的 重疊峰分峰效果圖Fig.3 Resolution result of overlapping peaks after three complete forward and backward fittings

    2 實(shí)驗(yàn)部分

    2.1 儀器及試劑

    GC-MS 7890B-5977A氣相色譜-質(zhì)譜聯(lián)用儀:美國Agilent公司產(chǎn)品;間二甲苯甲醇溶液、對(duì)二甲苯甲醇溶液(1 000 mg/L):北京壇墨質(zhì)檢科技有限公司產(chǎn)品;實(shí)驗(yàn)數(shù)據(jù):使用Matlab R2014a編寫的程序進(jìn)行計(jì)算。

    2.2 實(shí)驗(yàn)條件

    2.2.1色譜條件 色譜柱:DB-WAX柱(30 m×0.25 mm×0.25 μm);升溫程序:80 ℃保持4 min,以20 ℃/min升至200 ℃,保持1 min;載氣(He)流速1 mL/min;進(jìn)樣量0.5 μL;分流比20∶1。

    2.2.2質(zhì)譜條件 電子轟擊(EI)離子源;電子能量70 eV;傳輸線溫度240 ℃;離子源溫度230 ℃;溶劑延遲2 min。

    2.3 仿真實(shí)驗(yàn)設(shè)計(jì)

    根據(jù)塔板理論,指數(shù)修正高斯(EMG)模型可以很好地描述色譜流出曲線,因此,本工作采用EMG模型來構(gòu)造色譜峰[14]。影響色譜重疊峰峰形的3個(gè)主要因素為分離度、兩個(gè)單峰的峰高比、拖尾因子。分別針對(duì)這3個(gè)影響因素對(duì)重疊峰峰形進(jìn)行劃分,然后對(duì)比垂線法、交點(diǎn)垂線法、比例分配法和前后向擬合法在不同影響因素情況下的分離性能。

    1) 有谷點(diǎn)重疊峰:設(shè)定分離度為0.797,峰高比分別為1∶4、2∶3、1∶1、 3∶2、4∶1。

    2) 無谷點(diǎn)重疊峰,即肩峰:設(shè)定峰高比分別為4∶1和1∶4,分離度分別為0.478、0.345。

    3) 前延重疊峰和拖尾重疊峰:設(shè)定峰高比為4∶1,分離度為0.797,拖尾因子分別為0.766 2、0.853 8、0.975 1、1.031 6、1.172 5、1.417 2。

    以上實(shí)驗(yàn)設(shè)計(jì)中,常見重疊峰類型及重疊程度高低均包含在內(nèi),選取的參數(shù)比較具有代表性。

    2.4 真實(shí)實(shí)驗(yàn)設(shè)計(jì)

    為了更好地驗(yàn)證前后向擬合算法的可行性,選擇對(duì)二甲苯和間二甲苯進(jìn)行分離實(shí)驗(yàn)[12]。配制一系列不同濃度比的間二甲苯和對(duì)二甲苯的混合甲醇溶液,具體情況列于表1。為了評(píng)估算法的有效性,同時(shí)分別配制了一系列與各混合溶液中濃度相同的對(duì)二甲苯、間二甲苯標(biāo)準(zhǔn)溶液,以獲取兩者的單峰。

    3 結(jié)果與討論

    通常情況下,峰面積是定量分析的一種衡量標(biāo)準(zhǔn),因此,本工作通過計(jì)算真實(shí)峰面積與垂線法、交點(diǎn)垂線法、比例分配法、前后向擬合法所得到的單峰面積之間的誤差,評(píng)估各個(gè)方法的性能。誤差可由式(11)計(jì)算:

    (11)

    式中,Amethod是由垂線法、交點(diǎn)垂線法、比例分配法、前后向擬合法計(jì)算得到的單峰面積;AR是構(gòu)成重疊峰的每個(gè)單峰的面積。

    表1 對(duì)二甲苯和間二甲苯混合溶液配比Table 1 Mixture ratios of p-xylene and m-xylene

    3.1 仿真結(jié)果分析

    有谷點(diǎn)重疊峰和無谷點(diǎn)重疊峰的計(jì)算誤差列于表2。從表2可以看出:對(duì)于有谷點(diǎn)重疊峰(No.1~5),分離誤差會(huì)隨著峰高比的增大而增大,垂線法的分離誤差最大,交點(diǎn)垂線法和比例分配法次之,而前后向擬合方法的誤差處于±0.01%范圍內(nèi);對(duì)于無谷點(diǎn)重疊峰(No.6~9),垂線法是無法使用的,通過做重疊峰拐點(diǎn)處的切線交點(diǎn)求得重疊峰的兩個(gè)極大值點(diǎn)后,使用交點(diǎn)垂線法和比例分配法進(jìn)行分峰,得到的最高誤差分別可達(dá)18.93%和26.24%,但使用前后向擬合方法進(jìn)行計(jì)算,其最大誤差僅為0.51%;從No.1、No.8、No.9和No.5、No.6、No.7這兩組數(shù)據(jù)來看,在同一峰高比的情況下,隨著分離度的降低,重疊峰的重疊程度增大,垂線法、交點(diǎn)垂線法、比例分配法的分離誤差變化幅度較大,而前后向擬合方法的分峰誤差則相對(duì)平穩(wěn)地處于較小范圍之內(nèi)。同時(shí),對(duì)于No.7無谷點(diǎn)重疊峰,雖然其分離度很小,但前后向擬合方法的分峰輪廓與真實(shí)單峰幾乎重合,實(shí)驗(yàn)效果圖示于圖4。

    在實(shí)際色譜分離實(shí)驗(yàn)中,由于進(jìn)樣量大、色譜柱污染等原因,可能導(dǎo)致色譜峰拖尾,因此,本實(shí)驗(yàn)對(duì)4種方法分離拖尾峰的性能進(jìn)行了研究,結(jié)果列于表3。在同一峰高比和分離度的情況下,隨著重疊峰拖尾程度的增加,4種方法的分離誤差均有隨之增大的趨勢(shì),與表2中No.5無拖尾實(shí)驗(yàn)數(shù)據(jù)相比,這種現(xiàn)象也得到了驗(yàn)證。隨著拖尾程度的增大,前后向擬合方法的分峰誤差明顯小于其他三種方法,整體誤差穩(wěn)定在±0.11%以內(nèi)。前后向擬合方法分離No.10重疊峰實(shí)驗(yàn)效果圖示于圖5。由圖5可見,對(duì)于存在嚴(yán)重拖尾的重疊峰,前后向擬合方法不僅可以將誤差控制在較小范圍內(nèi),其輪廓也與真實(shí)單峰較一致。

    表2 有谷點(diǎn)重疊峰(分離度為0.797)和無谷點(diǎn)重疊峰(峰高比分別為4∶1、1∶4)的計(jì)算誤差Table 2 Resolution results of overlapping peaks with valley point under resolution of 0.797 and those of overlapping paks without valley point under peak height ratios 4∶1 and 1∶4, respectively

    3.2 真實(shí)結(jié)果分析

    為了進(jìn)一步驗(yàn)證前后向擬合方法的性能,使用氣相色譜-質(zhì)譜聯(lián)用儀開展了一組對(duì)二間苯和間二甲苯同分異構(gòu)體的重疊峰分離實(shí)驗(yàn),并分別采用4種方法對(duì)重疊峰進(jìn)行分離,結(jié)果列于表4。

    圖4 前后向擬合方法分離No.7重疊峰

    序號(hào)No.拖尾因子Trailing factor單峰Single peak誤差Error/%前后向擬合法垂線法交點(diǎn)垂線法比例分配法100.7662單峰1-0. 013.25-2.70-3.26單峰20. 00-13.0114.8317.07110.8538單峰10. 002.58-0.56-1.24單峰20. 00-10.312.254.97120.9751單峰10. 001.630.57-0.03單峰2-0. 01-6.54-2.260.13131.0316單峰10. 011.230.500.13單峰2-0. 02-4.923.304.78141.1725單峰10. 03-0. 271.550.47單峰2-0. 061.07-4.180.13151.4172單峰10.04-3.561.860.26單峰2-0. 1114.25-7.45-1.04

    圖5 前后向擬合方法分離No.10重疊峰實(shí)驗(yàn)效果圖Fig.5 Resolved profiles of No.10 leading overlapping peak by forward-backward fitting method

    從表4可以看出,對(duì)于Ms3溶液,垂線法的分離誤差比前后向擬合法的誤差小,這是因?yàn)榇咕€法更適用于分離兩個(gè)含量相同的單組分,此時(shí)兩個(gè)單峰受到重疊影響的程度相似,因此垂線法的分峰誤差最小,但前后向擬合法仍能將Ms3溶液的分離效果控制在0.8%以內(nèi)。對(duì)于其他4種比例的混合溶液,前后向擬合方法的整體誤差均小于其他3種方法。隨著兩種成分比例的增加,其他3種方法的分離誤差變化較大,而前后向擬合方法的誤差變化則比較平穩(wěn)。前后向擬合方法分離Ms1重疊峰的結(jié)果示于圖6。

    表4 對(duì)二甲苯和間二甲苯的重疊峰分離結(jié)果Table 4 Resolution results of overlapping peaks of p-xylene and m-xylene

    圖6 前后向擬合方法分離Ms1重疊峰效果圖Fig.6 Resolved profiles of Ms1 overlapping peak by forward-backward fitting method

    4 結(jié)論

    根據(jù)重疊峰相似性原理,提出了前后向擬合分峰方法,通過仿真實(shí)驗(yàn)和真實(shí)實(shí)驗(yàn),與垂線法、交點(diǎn)垂線法、比例分配法的分離精度進(jìn)行對(duì)比。結(jié)果表明,前后向擬合方法的分離誤差可穩(wěn)定在1.8%以內(nèi),分離精度高于常規(guī)的3種分離方法,而且無需依賴兩個(gè)極大值點(diǎn),且能適應(yīng)多種常見的重疊峰類型,即使對(duì)于重疊程度較高或具有拖尾的重疊峰,前后向擬合方法也具有較好的分峰效果。

    猜你喜歡
    峰高單峰極大值
    金字塔
    圖的完全積的獨(dú)立數(shù)和獨(dú)立多項(xiàng)式
    老君山
    許慶賢
    書香兩岸(2020年3期)2020-06-29 12:33:45
    陳初良
    書香兩岸(2020年3期)2020-06-29 12:33:45
    Kirchhoff方程單峰解的局部唯一性
    基于小波模極大值理論的勵(lì)磁涌流新判據(jù)研究
    基于經(jīng)驗(yàn)?zāi)B(tài)分解的自適應(yīng)模極大值去噪方法
    行人檢測中非極大值抑制算法的改進(jìn)
    關(guān)于單峰偏好的注記
    99久久99久久久精品蜜桃| 欧美黄色片欧美黄色片| 欧美乱码精品一区二区三区| 在线观看免费午夜福利视频| 丝袜美腿诱惑在线| 一级a爱片免费观看的视频| 国产乱人伦免费视频| 国产亚洲精品av在线| svipshipincom国产片| 欧美人与性动交α欧美精品济南到| 久久久久久久精品吃奶| 国产又色又爽无遮挡免费看| 成人特级黄色片久久久久久久| 韩国av一区二区三区四区| 怎么达到女性高潮| 国产激情久久老熟女| 免费在线观看黄色视频的| 国产精品久久视频播放| 亚洲久久久国产精品| www日本黄色视频网| 又紧又爽又黄一区二区| 人妻丰满熟妇av一区二区三区| 精品日产1卡2卡| 一个人观看的视频www高清免费观看 | 精品久久久久久成人av| 欧美乱妇无乱码| 亚洲av第一区精品v没综合| 在线永久观看黄色视频| 一区二区三区高清视频在线| 99国产精品一区二区三区| 特大巨黑吊av在线直播 | 午夜免费鲁丝| 久久久久精品国产欧美久久久| 日韩高清综合在线| 久久久久久久精品吃奶| 亚洲精品美女久久av网站| 嫩草影视91久久| 国产精品98久久久久久宅男小说| 国产又黄又爽又无遮挡在线| 桃色一区二区三区在线观看| 男女视频在线观看网站免费 | 午夜视频精品福利| 精品国产国语对白av| 天堂影院成人在线观看| www国产在线视频色| 欧美丝袜亚洲另类 | 午夜a级毛片| 国产亚洲精品综合一区在线观看 | 夜夜夜夜夜久久久久| 一级片免费观看大全| 久久中文看片网| 女人高潮潮喷娇喘18禁视频| 亚洲国产高清在线一区二区三 | 男女下面进入的视频免费午夜 | 国产又色又爽无遮挡免费看| 12—13女人毛片做爰片一| 757午夜福利合集在线观看| 日韩大尺度精品在线看网址| 成年女人毛片免费观看观看9| 国产熟女午夜一区二区三区| 午夜亚洲福利在线播放| 亚洲成av片中文字幕在线观看| 亚洲国产精品999在线| 丝袜在线中文字幕| 亚洲专区国产一区二区| 国产亚洲精品av在线| 国产成人av激情在线播放| 国产亚洲欧美精品永久| 在线av久久热| 成人手机av| a级毛片在线看网站| 国内毛片毛片毛片毛片毛片| 丰满的人妻完整版| ponron亚洲| 熟女电影av网| 99国产精品一区二区三区| 亚洲中文日韩欧美视频| 日日摸夜夜添夜夜添小说| 精品久久久久久久久久久久久 | 中文字幕高清在线视频| 一本精品99久久精品77| 精品午夜福利视频在线观看一区| 香蕉丝袜av| av在线天堂中文字幕| 国产一区二区激情短视频| 久久亚洲精品不卡| 国产真人三级小视频在线观看| 亚洲专区中文字幕在线| 久久精品人妻少妇| 一级a爱视频在线免费观看| 不卡av一区二区三区| 黑人操中国人逼视频| 制服人妻中文乱码| 欧美国产日韩亚洲一区| 国产精品久久电影中文字幕| 女人被狂操c到高潮| 老司机午夜十八禁免费视频| 给我免费播放毛片高清在线观看| 啪啪无遮挡十八禁网站| 亚洲一码二码三码区别大吗| 露出奶头的视频| 欧美激情极品国产一区二区三区| 久久精品影院6| 亚洲av成人av| 精品乱码久久久久久99久播| 久久精品国产清高在天天线| 免费高清视频大片| 久久久国产欧美日韩av| 久久国产精品影院| 悠悠久久av| 在线视频色国产色| 国产精品久久久久久精品电影 | 夜夜夜夜夜久久久久| 大香蕉久久成人网| 一边摸一边抽搐一进一小说| 三级毛片av免费| 亚洲国产精品久久男人天堂| 激情在线观看视频在线高清| www日本黄色视频网| 香蕉丝袜av| 给我免费播放毛片高清在线观看| 高清在线国产一区| 国产免费男女视频| 精品人妻1区二区| 国产私拍福利视频在线观看| 国产精品日韩av在线免费观看| 听说在线观看完整版免费高清| 2021天堂中文幕一二区在线观 | 丝袜美腿诱惑在线| 在线天堂中文资源库| 在线免费观看的www视频| 国产蜜桃级精品一区二区三区| av欧美777| 精品国产亚洲在线| 韩国av一区二区三区四区| 特大巨黑吊av在线直播 | 精品一区二区三区视频在线观看免费| 欧美日韩乱码在线| 国产av一区二区精品久久| 精品国产美女av久久久久小说| 精品国产一区二区三区四区第35| 精品一区二区三区视频在线观看免费| 亚洲人成网站在线播放欧美日韩| 婷婷精品国产亚洲av在线| 国产蜜桃级精品一区二区三区| 听说在线观看完整版免费高清| 久久久久九九精品影院| 熟女少妇亚洲综合色aaa.| 精品久久久久久久久久免费视频| 精品久久久久久久久久免费视频| 一卡2卡三卡四卡精品乱码亚洲| 可以在线观看的亚洲视频| 露出奶头的视频| 欧美激情极品国产一区二区三区| avwww免费| 亚洲av电影在线进入| 俄罗斯特黄特色一大片| 1024视频免费在线观看| 国产成人欧美在线观看| 国产区一区二久久| www国产在线视频色| 免费看日本二区| 久久久久久久久免费视频了| 国产1区2区3区精品| 99久久99久久久精品蜜桃| 中文资源天堂在线| 搡老熟女国产l中国老女人| 亚洲一区二区三区不卡视频| 又黄又粗又硬又大视频| 欧美 亚洲 国产 日韩一| 久久精品亚洲精品国产色婷小说| 欧美精品啪啪一区二区三区| 久久久久免费精品人妻一区二区 | 亚洲在线自拍视频| 国产片内射在线| 亚洲第一电影网av| 亚洲精品国产一区二区精华液| 午夜成年电影在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 国产单亲对白刺激| 两人在一起打扑克的视频| 久久精品成人免费网站| 亚洲精品久久国产高清桃花| 久久天躁狠狠躁夜夜2o2o| 老鸭窝网址在线观看| 亚洲精品久久国产高清桃花| 日韩大码丰满熟妇| 免费看十八禁软件| 久99久视频精品免费| 久9热在线精品视频| 久久这里只有精品19| 久久草成人影院| 亚洲色图 男人天堂 中文字幕| 18美女黄网站色大片免费观看| 久久久久九九精品影院| 91国产中文字幕| 国产欧美日韩一区二区精品| 精品一区二区三区四区五区乱码| 国产v大片淫在线免费观看| 欧美中文日本在线观看视频| 亚洲第一电影网av| 国产av一区在线观看免费| www.熟女人妻精品国产| 久久国产精品人妻蜜桃| 男女那种视频在线观看| 黑人操中国人逼视频| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久久久精品电影 | 国产激情欧美一区二区| 亚洲精品中文字幕在线视频| 午夜福利一区二区在线看| 欧美激情高清一区二区三区| 18禁黄网站禁片午夜丰满| 亚洲欧美日韩高清在线视频| 少妇粗大呻吟视频| 国产亚洲欧美在线一区二区| 国产精品久久视频播放| √禁漫天堂资源中文www| 最近最新中文字幕大全电影3 | www.自偷自拍.com| 美女高潮到喷水免费观看| 欧美日韩福利视频一区二区| 级片在线观看| 美女国产高潮福利片在线看| 国产精品影院久久| 国产一区二区在线av高清观看| 真人一进一出gif抽搐免费| 国产在线精品亚洲第一网站| 成人av一区二区三区在线看| 国产精品九九99| 制服诱惑二区| 夜夜躁狠狠躁天天躁| 免费在线观看影片大全网站| 91麻豆精品激情在线观看国产| 欧美zozozo另类| 精品日产1卡2卡| 老熟妇乱子伦视频在线观看| av欧美777| 性色av乱码一区二区三区2| 热99re8久久精品国产| 高潮久久久久久久久久久不卡| 欧美日韩黄片免| 国产午夜精品久久久久久| 中文字幕人妻熟女乱码| 淫妇啪啪啪对白视频| 中文亚洲av片在线观看爽| 亚洲激情在线av| 亚洲成av人片免费观看| 精品国产乱子伦一区二区三区| 国产亚洲精品久久久久久毛片| 午夜福利在线在线| 国产男靠女视频免费网站| 免费高清在线观看日韩| 亚洲精华国产精华精| 国产一级毛片七仙女欲春2 | av欧美777| 国产又爽黄色视频| 亚洲第一av免费看| 欧美三级亚洲精品| 看黄色毛片网站| www.精华液| 成熟少妇高潮喷水视频| 精品久久久久久成人av| 50天的宝宝边吃奶边哭怎么回事| 精品久久久久久久久久免费视频| 日韩大尺度精品在线看网址| 老司机午夜十八禁免费视频| 老司机午夜福利在线观看视频| 校园春色视频在线观看| 一级毛片精品| 一本大道久久a久久精品| 妹子高潮喷水视频| 亚洲成人精品中文字幕电影| 1024香蕉在线观看| 啦啦啦 在线观看视频| 日韩视频一区二区在线观看| 亚洲av电影在线进入| 一夜夜www| 麻豆国产av国片精品| 成年人黄色毛片网站| 深夜精品福利| 亚洲专区中文字幕在线| 亚洲av成人不卡在线观看播放网| 久久精品国产清高在天天线| 午夜福利欧美成人| 久热这里只有精品99| 久久中文看片网| 一个人免费在线观看的高清视频| 99riav亚洲国产免费| 国产成人一区二区三区免费视频网站| 国产亚洲欧美98| 九色国产91popny在线| 看免费av毛片| 国产精品野战在线观看| 精品久久久久久久久久久久久 | 亚洲成国产人片在线观看| 91大片在线观看| 天天一区二区日本电影三级| 欧美 亚洲 国产 日韩一| 超碰成人久久| 日韩中文字幕欧美一区二区| 99久久国产精品久久久| 欧美色欧美亚洲另类二区| 一本精品99久久精品77| 中文资源天堂在线| 不卡一级毛片| 国产免费av片在线观看野外av| 国产精品免费一区二区三区在线| 51午夜福利影视在线观看| 免费高清在线观看日韩| 国产成人精品无人区| 这个男人来自地球电影免费观看| 一级片免费观看大全| 老鸭窝网址在线观看| 人人妻人人澡人人看| 成人三级做爰电影| 老司机福利观看| 最新美女视频免费是黄的| 热re99久久国产66热| 日本黄色视频三级网站网址| 久久婷婷人人爽人人干人人爱| 久久精品影院6| 啦啦啦观看免费观看视频高清| 99国产精品一区二区三区| 亚洲五月色婷婷综合| 男人舔奶头视频| 午夜免费观看网址| 日日干狠狠操夜夜爽| 男女之事视频高清在线观看| 一级片免费观看大全| 日韩有码中文字幕| 亚洲五月婷婷丁香| 国产97色在线日韩免费| 成人一区二区视频在线观看| 中文字幕av电影在线播放| 91九色精品人成在线观看| 国产精品1区2区在线观看.| 婷婷精品国产亚洲av在线| 久久中文看片网| www.www免费av| 好男人在线观看高清免费视频 | 国产亚洲精品一区二区www| 亚洲国产日韩欧美精品在线观看 | 女同久久另类99精品国产91| 国产亚洲精品第一综合不卡| 国产三级黄色录像| 国产精品久久视频播放| x7x7x7水蜜桃| 色婷婷久久久亚洲欧美| 久久欧美精品欧美久久欧美| 色综合站精品国产| 99riav亚洲国产免费| 欧美黄色淫秽网站| 一级毛片女人18水好多| 老司机深夜福利视频在线观看| 大型黄色视频在线免费观看| 岛国视频午夜一区免费看| av超薄肉色丝袜交足视频| 成人亚洲精品av一区二区| 欧美日韩亚洲国产一区二区在线观看| 精品卡一卡二卡四卡免费| 久久久久久人人人人人| 欧美av亚洲av综合av国产av| 亚洲激情在线av| 两个人视频免费观看高清| av福利片在线| 观看免费一级毛片| 国产一级毛片七仙女欲春2 | 怎么达到女性高潮| 免费看a级黄色片| 99riav亚洲国产免费| 亚洲av成人一区二区三| 久热这里只有精品99| 色综合欧美亚洲国产小说| 国产亚洲精品一区二区www| 美女高潮喷水抽搐中文字幕| 欧美色视频一区免费| 黑人欧美特级aaaaaa片| 亚洲国产欧美网| 一级a爱片免费观看的视频| 日韩av在线大香蕉| 亚洲av电影不卡..在线观看| 久久 成人 亚洲| 大型黄色视频在线免费观看| 久久香蕉国产精品| 麻豆av在线久日| 91大片在线观看| 脱女人内裤的视频| www.自偷自拍.com| 久久亚洲精品不卡| 久久久久久久精品吃奶| 精品一区二区三区四区五区乱码| 精品国内亚洲2022精品成人| 无限看片的www在线观看| 国产乱人伦免费视频| 久久99热这里只有精品18| 成人午夜高清在线视频 | 欧美亚洲日本最大视频资源| 一二三四社区在线视频社区8| 国产色视频综合| 波多野结衣巨乳人妻| 日韩成人在线观看一区二区三区| 少妇裸体淫交视频免费看高清 | 亚洲成人国产一区在线观看| 满18在线观看网站| 一级作爱视频免费观看| 精品第一国产精品| 国产高清视频在线播放一区| 国产亚洲精品久久久久5区| 免费在线观看成人毛片| 色综合婷婷激情| 色在线成人网| 国产黄a三级三级三级人| 成年人黄色毛片网站| 9191精品国产免费久久| 成人三级黄色视频| 久久久久国内视频| 可以在线观看毛片的网站| 午夜免费成人在线视频| 性色av乱码一区二区三区2| 国产精品永久免费网站| 婷婷亚洲欧美| 精品久久久久久,| 90打野战视频偷拍视频| 国产精品永久免费网站| 脱女人内裤的视频| 精品人妻1区二区| 最近最新中文字幕大全免费视频| 两个人看的免费小视频| 老司机靠b影院| 90打野战视频偷拍视频| 亚洲欧美精品综合久久99| 精品福利观看| 免费电影在线观看免费观看| 嫁个100分男人电影在线观看| 久久青草综合色| 自线自在国产av| 亚洲九九香蕉| 熟女少妇亚洲综合色aaa.| 淫妇啪啪啪对白视频| 欧美日本视频| 久久狼人影院| 1024手机看黄色片| 黑人操中国人逼视频| 99re在线观看精品视频| 精品欧美国产一区二区三| 无限看片的www在线观看| 18禁国产床啪视频网站| 亚洲欧美激情综合另类| 久久国产精品影院| 91成年电影在线观看| 国产精品 欧美亚洲| 嫩草影视91久久| a级毛片在线看网站| 日日爽夜夜爽网站| 在线观看一区二区三区| 大香蕉久久成人网| 91九色精品人成在线观看| 欧美午夜高清在线| 日本一区二区免费在线视频| 亚洲一区中文字幕在线| av有码第一页| 欧美成人免费av一区二区三区| 国产一级毛片七仙女欲春2 | 成人欧美大片| 午夜a级毛片| 一本久久中文字幕| 757午夜福利合集在线观看| 88av欧美| 性欧美人与动物交配| 日本黄色视频三级网站网址| 久久中文字幕人妻熟女| 熟女少妇亚洲综合色aaa.| 999久久久国产精品视频| 欧美在线一区亚洲| 国产1区2区3区精品| 欧美黑人精品巨大| 欧美日韩一级在线毛片| 欧美日韩精品网址| 少妇裸体淫交视频免费看高清 | 亚洲五月婷婷丁香| 国产黄片美女视频| av有码第一页| 丰满人妻熟妇乱又伦精品不卡| 999久久久国产精品视频| 俺也久久电影网| 一本大道久久a久久精品| 亚洲精品色激情综合| 性欧美人与动物交配| 老司机靠b影院| 夜夜爽天天搞| 久久久久免费精品人妻一区二区 | 一级片免费观看大全| 国产激情偷乱视频一区二区| 国产欧美日韩一区二区三| 中文字幕精品亚洲无线码一区 | www.999成人在线观看| 国产区一区二久久| 亚洲精品在线美女| 两个人视频免费观看高清| 看黄色毛片网站| 美女高潮喷水抽搐中文字幕| 欧美日韩亚洲国产一区二区在线观看| 自线自在国产av| 丝袜人妻中文字幕| 国产亚洲精品av在线| 日韩三级视频一区二区三区| 亚洲国产精品久久男人天堂| 女人高潮潮喷娇喘18禁视频| 国产欧美日韩一区二区精品| 好看av亚洲va欧美ⅴa在| 在线天堂中文资源库| 成人亚洲精品av一区二区| 午夜福利在线在线| 制服诱惑二区| 成年版毛片免费区| 91av网站免费观看| 国产99久久九九免费精品| 国产蜜桃级精品一区二区三区| 日韩一卡2卡3卡4卡2021年| 丝袜在线中文字幕| 日韩高清综合在线| 在线十欧美十亚洲十日本专区| 成人亚洲精品一区在线观看| 欧美大码av| 日韩精品免费视频一区二区三区| 免费在线观看成人毛片| 丝袜在线中文字幕| 久久久久九九精品影院| 亚洲午夜精品一区,二区,三区| 欧美性猛交黑人性爽| 久久人妻av系列| 久久狼人影院| 首页视频小说图片口味搜索| 每晚都被弄得嗷嗷叫到高潮| 一区二区三区精品91| 亚洲专区国产一区二区| 国产精品av久久久久免费| 精品久久久久久久人妻蜜臀av| 中文字幕人妻丝袜一区二区| 国产v大片淫在线免费观看| 人人妻,人人澡人人爽秒播| 在线观看免费午夜福利视频| 国产一区二区三区视频了| 在线观看免费日韩欧美大片| 黄色成人免费大全| 久久精品国产亚洲av高清一级| av电影中文网址| 久久久久久国产a免费观看| 欧美一级a爱片免费观看看 | 亚洲片人在线观看| 欧美一级毛片孕妇| 婷婷六月久久综合丁香| 久久人妻av系列| 天堂动漫精品| 免费一级毛片在线播放高清视频| 欧美性猛交╳xxx乱大交人| 亚洲va日本ⅴa欧美va伊人久久| 69av精品久久久久久| 精品少妇一区二区三区视频日本电影| 91国产中文字幕| 久久久久久九九精品二区国产 | 亚洲av电影不卡..在线观看| 午夜福利一区二区在线看| 久久国产精品人妻蜜桃| 啦啦啦韩国在线观看视频| 19禁男女啪啪无遮挡网站| 精品久久久久久久人妻蜜臀av| 欧美性猛交╳xxx乱大交人| 欧美日韩瑟瑟在线播放| 亚洲国产精品成人综合色| 侵犯人妻中文字幕一二三四区| 18禁裸乳无遮挡免费网站照片 | 免费av毛片视频| 亚洲欧洲精品一区二区精品久久久| 久久精品国产亚洲av香蕉五月| www.熟女人妻精品国产| 宅男免费午夜| 一本精品99久久精品77| 免费电影在线观看免费观看| 黄色a级毛片大全视频| 19禁男女啪啪无遮挡网站| 亚洲午夜理论影院| 性色av乱码一区二区三区2| 757午夜福利合集在线观看| 母亲3免费完整高清在线观看| 亚洲专区国产一区二区| 亚洲五月婷婷丁香| 老熟妇仑乱视频hdxx| x7x7x7水蜜桃| 欧美av亚洲av综合av国产av| 成人亚洲精品av一区二区| 宅男免费午夜| 免费观看精品视频网站| 久久久久久大精品| 美女国产高潮福利片在线看| tocl精华| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲天堂国产精品一区在线| 女生性感内裤真人,穿戴方法视频| 不卡av一区二区三区| 老鸭窝网址在线观看| 中文亚洲av片在线观看爽| 青草久久国产| 999精品在线视频| 亚洲成人久久爱视频| 神马国产精品三级电影在线观看 | 欧美日韩福利视频一区二区| 亚洲 国产 在线| 午夜激情福利司机影院| 波多野结衣av一区二区av| 99久久综合精品五月天人人| 久久精品国产亚洲av香蕉五月| 久久久久久久久免费视频了| 亚洲精品美女久久av网站|