• <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)絡的天然氣水合物飽和度軟測量模型

    2020-06-29 12:13:40曹勝昌邢蘭昌偉2韓維峰2牛佳樂
    計算機測量與控制 2020年6期
    關(guān)鍵詞:模值水合物飽和度

    曹勝昌,邢蘭昌,魏 偉2,韓維峰2,牛佳樂

    (1.中國石油大學(華東)控制科學與工程學院,山東 青島 266580)(2.中國石油勘探開發(fā)研究院新能源研究所,河北 廊坊 065007)

    0 引言

    天然氣水合物又稱為“可燃冰”,是一種新型優(yōu)質(zhì)的潛在能源,是在低溫、高壓等條件下形成的一種類冰狀、具有籠形結(jié)構(gòu)的固態(tài)化合物[1-2],通常分布在大陸永久凍土帶和海洋沉積物中,具有非常巨大的能源和經(jīng)濟價值[3-4]。天然氣水合物飽和度是水合物儲層評價的主要參數(shù)之一,為認識水合物富集成藏規(guī)律、估算水合物資源量提供了必要的基礎(chǔ)參數(shù)[5-7]。

    水合物飽和度是含水合物沉積物電學特性的直接影響因素之一。周錫唐等[8]研究了甲烷水合物生成與分解過程中電導性的變化,實驗結(jié)果表明電阻隨石英砂中水合物生成量的增加而增大,隨水合物的分解而急速減小,通過參考溫度和壓力的同步變化特征證明了使用電導性來監(jiān)測甲烷水合物生成與分解過程的可靠性。陳玉鳳等[9-10]研究了南海沉積物中甲烷水合物飽和度與電阻率的關(guān)系,發(fā)現(xiàn)在水合物生成過程中,其電阻率會隨著水合物飽和度的增加而增大,而且含水合物沉積物存在非阿爾奇現(xiàn)象,其飽和度指數(shù)隨水合物飽和度的增加而增加。

    基于電學特性參數(shù)計算含水合物飽和度的模型主要有阿爾奇公式、含泥質(zhì)修正的阿爾奇公式、印度尼西亞公式和雙水模型等,其中阿爾奇公式和含泥質(zhì)修正的阿爾奇公式被廣泛應用[11-13]。Riedel等[14]通過分析不同地區(qū)的測井數(shù)據(jù),得到了不同的阿爾奇公式經(jīng)驗參數(shù)并估算了水合物的飽和度,指出不同地區(qū)阿爾奇公式經(jīng)驗參數(shù)的差異會導致水合物飽和度預測結(jié)果的不同。Spangenberg等[15-17]首先從理論上研究了水合物含量對體系電阻率性質(zhì)的影響,之后采用溶解氣在均勻玻璃珠內(nèi)開展了甲烷水合物生成實驗,測量結(jié)果顯示電阻率增大系數(shù)與含水飽和度之間的關(guān)系不能用恒定飽和度指數(shù)的阿爾奇公式來描述。含水合物沉積物的電學性質(zhì)受到沉積物孔隙結(jié)構(gòu)、孔隙水化學成分及濃度、孔隙度、飽和度、溫度等多種因素的影響,水合物飽和度與電學性質(zhì)參數(shù)呈現(xiàn)出顯著的復雜非線性、多值對應關(guān)系,以阿爾奇公式為代表的經(jīng)驗公式的適用性受到極大局限。

    神經(jīng)網(wǎng)絡具有強大的處理復雜非線性問題的能力以及高度自學習和自適應的能力,運用神經(jīng)網(wǎng)絡來建立預測模型成為一種有潛力的水合物飽和度計算方法。呂琳[18]等人在使用地球物理測井數(shù)據(jù)預測儲量參數(shù)方法的基礎(chǔ)上,選用中子測井、密度測井等8種測井數(shù)據(jù)作為BP神經(jīng)網(wǎng)絡的輸入量,預測了天然氣水合物的儲量參數(shù)(孔隙度、飽和度),結(jié)果表明該模型計算結(jié)果的精度高于傳統(tǒng)的經(jīng)驗公式。陸敬安[19]等利用神狐海域某井的電阻率測井數(shù)據(jù)作為樣本訓練并建立了徑向基神經(jīng)網(wǎng)絡模型,利用該模型計算出另一口井的水合物飽和度,其計算結(jié)果與現(xiàn)場孔隙水分析得到的飽和度值相吻合。現(xiàn)有的研究存在以下問題:(1)僅僅利用單個神經(jīng)網(wǎng)絡建立天然氣水合物飽和度的計算模型,模型的可靠性有待進一步驗證與提高;(2)神經(jīng)網(wǎng)絡的輸入?yún)?shù)缺少必要的預處理,導致神經(jīng)網(wǎng)絡模型的準確度受到限制。

    本文提出基于集成神經(jīng)網(wǎng)絡方法建立一種以電阻抗譜特性參數(shù)為特征向量的水合物飽和度軟測量模型,即將不同頻率的阻抗模值作為單個BP神經(jīng)網(wǎng)絡的輸入,通過對不同位置傳感器所對應的子網(wǎng)絡輸出進行集成,建立基于集成BP神經(jīng)網(wǎng)絡的水合物飽和度計算模型。

    1 天然氣水合物模擬實驗

    采用了自主設(shè)計的天然氣水合物電學探測實驗裝置開展實驗工作。實驗裝置的總體結(jié)構(gòu)如圖1所示,該裝置分為環(huán)境模擬部分和測試部分。環(huán)境模擬部分包括耐高壓反應釜、低溫冷卻系統(tǒng)、天然氣氣瓶、中間儲氣罐、軸壓系統(tǒng)等,可以模擬出天然氣水合物的生成分解所需的溫度和壓力條件、提供充足的氣源以及必要的存儲空間等;測試部分主要包括電學傳感器、溫度傳感器、壓力傳感器、電學測試單元、溫度測試單元、壓力測試單元、工控機及測試軟件,能夠?qū)崟r獲取反映水合物生成分解過程動態(tài)變化的實驗數(shù)據(jù)。

    圖1 實驗裝置總體結(jié)構(gòu)

    天然氣水合物在反應釜內(nèi)呈現(xiàn)出空間不均勻分布的特性,為了獲取水合物的空間分布信息,將電學傳感器分為上下兩層均勻安置在反應釜壁,每一層安裝8個傳感器(如圖2所示),采用分時輪流工作模式控制各傳感器的工作狀態(tài)。

    圖2 單層傳感器排列方式圖

    利用上述裝置開展了甲烷水合物生成分解模擬實驗,分析了多孔介質(zhì)(粒徑50~180 μm的天然海沙)中甲烷水合物生成和分解過程中電學特性的變化規(guī)律,以及水合物飽和度對含水合物的多孔介質(zhì)的電學特性影響,在此基礎(chǔ)上利用阻抗模值等特征參數(shù)建立起特征參數(shù)與飽和度之間的軟測量模型。

    2 神經(jīng)網(wǎng)絡樣本集的確定

    2.1 數(shù)據(jù)預處理

    實驗獲得的阻抗譜包括了從1 Hz~10 MHz范圍內(nèi)的22個頻率點,但由于系統(tǒng)誤差、環(huán)境干擾等因素影響,其中一些頻率點的數(shù)據(jù)會出現(xiàn)較大的無規(guī)律波動和明顯偏差,經(jīng)過篩選、剔除等步驟對數(shù)據(jù)進行選擇。最終選擇出5 Hz、10 Hz、20 Hz、100 Hz、200 Hz、500 Hz、1 kHz、2 kHz、4 kHz、10 kHz、20 kHz、40 kHz等12個頻率點作為特征頻率點,以這些測試頻率點的阻抗模值作為備選的特征參數(shù)。

    上述特征頻率點的原始數(shù)據(jù)中往往存在噪聲值、缺失值以及離群點等情況,不能直接用于軟測量模型的建立與測試,因此需要對原始數(shù)據(jù)進行預處理。數(shù)據(jù)預處理主要針對數(shù)據(jù)在數(shù)值上的各種異常情況進行處理。本文中數(shù)據(jù)預處理所采用的步驟為:數(shù)據(jù)離群點修正處理、數(shù)據(jù)去噪處理和數(shù)據(jù)歸一化處理。

    數(shù)據(jù)離群點修正處理針對原始數(shù)據(jù)中遠離一般水平的極端大值和極端小值,采用均值插補修正法,根據(jù)數(shù)據(jù)曲線的整體趨勢對數(shù)據(jù)進行估計值的修改,即用異常值數(shù)據(jù)所處位置左右相鄰兩個數(shù)據(jù)的平均值來替代異常值。本文采用小波去噪和滑動平均濾波去噪相結(jié)合的方式對數(shù)據(jù)進行去噪處理,以便獲得更好的去噪效果。數(shù)據(jù)歸一化處理將原始數(shù)據(jù)映射到0~1范圍之內(nèi),如式(1):

    (1)

    2.2 特征提取

    確定特征頻率點后,需要針對每組傳感器對在多頻率下的阻抗模值進行特征提取,即降維處理。通常神經(jīng)網(wǎng)絡輸入為多維數(shù)據(jù),計算復雜度隨著數(shù)據(jù)維數(shù)的增大成指數(shù)增加[20],數(shù)據(jù)維數(shù)過大會導致網(wǎng)絡訓練時間大大增加,也會導致訓練后的神經(jīng)網(wǎng)絡性能變差。

    本文采用主成分分析法(PCA)進行降維處理,消除不同頻率間阻抗模值之間的相關(guān)性,降低網(wǎng)絡輸入的維數(shù)、精簡網(wǎng)絡的結(jié)構(gòu),采用累計貢獻度達到99%的主成分作為神經(jīng)網(wǎng)絡的輸入。

    2.3 形成樣本集

    2.3.1 水合物飽和度計算

    在甲烷水合物模擬實驗中,基于耗氣量原理計算水合物飽和度,即利用實時采集的反應釜內(nèi)溫度和壓力數(shù)據(jù)并結(jié)合式(2)進行計算:

    (2)

    式中,Sh為水合物飽和度;P1為反應釜內(nèi)初始壓力;P2為水合物生成分解過程中反應釜內(nèi)壓力;T1為反應釜內(nèi)初始溫度,T2為水合物生成分解過程中反應釜內(nèi)溫度;Z1和Z2為對應壓力溫度條件下的氣體壓縮因子;Vg為反應釜氣相體積;Mh為水合物摩爾質(zhì)量,122.02 g/mol;R為氣體常數(shù),8.314 J/(mol·K);ρh為水合物的密度,0.91 g/mL;Vp為多孔介質(zhì)孔隙的體積。圖3為一輪典型實驗中水合物生成分解過程中溫度、壓力以及水合物飽和度的變化曲線。

    在水合物生成分解過程中,上層1-5傳感器對在部分特征頻率點(20 Hz、200 Hz、2 kHz、20 kHz)下的阻抗模值隨水合物飽和度變化的關(guān)系如圖4所示。通過分析阻抗模值參數(shù)與飽和度參數(shù)之間的關(guān)系可以發(fā)現(xiàn),阻抗模值隨著飽和度的增加逐漸增大,呈現(xiàn)出顯著的非線性。

    圖4 特征頻率點下阻抗模值隨水合物飽和度的變化圖

    2.3.2 劃分樣本集

    為了得到高效的神經(jīng)網(wǎng)絡模型,使之能夠有效地逼近訓練樣本所蘊含的內(nèi)在規(guī)律,需要將數(shù)據(jù)集劃分為訓練集和測試集。本文采用2/3的樣本作為訓練樣本來訓練網(wǎng)絡,1/3的樣本作為測試樣本來測試網(wǎng)絡,采用等間隔的方式從總樣本集中取出。

    3 水合物飽和度軟測量模型

    3.1 軟測量模型設(shè)計

    本文基于BP神經(jīng)網(wǎng)絡建立水合物飽和度軟測量模型。設(shè)計步驟如下:1)確定BP子神經(jīng)網(wǎng)絡最優(yōu)隱含層節(jié)點數(shù),進而確定最優(yōu)網(wǎng)絡結(jié)構(gòu)和訓練參數(shù);2)訓練并測試子神經(jīng)網(wǎng)絡,根據(jù)多種誤差指標評價網(wǎng)絡性能;3)采用平均法集成策略,建立集成BP神經(jīng)網(wǎng)絡的水合物飽和度軟測量模型。

    3.2 BP網(wǎng)絡結(jié)構(gòu)和訓練參數(shù)

    神經(jīng)網(wǎng)絡理論證明,一個三層結(jié)構(gòu)的BP網(wǎng)絡可實現(xiàn)任意m維到n維的映射,并且通過調(diào)整隱含層神經(jīng)元的個數(shù)來調(diào)節(jié)網(wǎng)絡的精度[21]。因此,本文選用三層BP神經(jīng)網(wǎng)絡。網(wǎng)絡的輸入為一組傳感器對在不同測試頻率下的阻抗模值,輸出為水合物飽和度。隱含層節(jié)點數(shù)的確定過程中參考式(3):

    (3)

    式中,L為隱含層節(jié)點數(shù),n為輸入層節(jié)點數(shù),m為輸出層節(jié)點數(shù),a為1~10內(nèi)的常數(shù)。

    在分析評價網(wǎng)絡的性能時采用以下4種誤差指標[22]:相對誤差(RE)反映測量結(jié)果偏離真值的實際大??;平均相對誤差(MRE)反映誤差的相對值,效果直觀;平均絕對誤差(MAE)反映預測值誤差的實際情況;均方根誤差(RMSE)反映網(wǎng)絡預測值相對于實際值的平均離散程度。通過分析比較以上誤差,可以確定更加合理地隱含層節(jié)點數(shù)。

    在網(wǎng)絡訓練時,本文中隱含層的激活函數(shù)采用正切Sigmoid函數(shù)(tansig),其輸入值可取任意值,輸出值在-1~+1之間;輸出層的激活函數(shù)采用線性函數(shù)(purelin),其輸入值與輸出值可取任意值;網(wǎng)絡學習方法采用動量梯度下降法(traingdm),該方法不但具有較高收斂速度,而且通過引入一個動量項可以避免局部極小問題。

    以上層1-5傳感器對為例(參見圖2),將隱含層節(jié)點數(shù)L取值由4增加到13,分別對不同隱含層節(jié)點數(shù)的網(wǎng)絡進行500次訓練和測試,通過分析比較網(wǎng)絡誤差來獲得最優(yōu)隱含層節(jié)點數(shù)。表1所示為500次測試誤差的平均結(jié)果,各種誤差指標的整體趨勢為先減小后增大,普遍存在一個最小值。當L=9時,平均相對誤差和均方根誤差兩項指標達到最優(yōu)值,雖然沒有實現(xiàn)全部指標都達到最優(yōu)值,但此時指標取值與最優(yōu)值相差較小。綜合考慮,隱含層節(jié)點數(shù)設(shè)為9,從而獲得針對上層1-5傳感器對的水合物飽和度計算模型。采用同樣的方法,根據(jù)表2、表3和表4的網(wǎng)絡測試誤差來確定針對2-6、3-7、4-8等傳感器對的網(wǎng)絡模型最優(yōu)隱含層節(jié)點數(shù),分別為9、10和10。

    表1 上層1-5傳感器對不同隱含層節(jié)點數(shù)網(wǎng)絡測試誤差表

    表2 上層2-6傳感器對不同隱含層節(jié)點數(shù)網(wǎng)絡測試誤差表

    表3 上層3-7傳感器對不同隱含層節(jié)點數(shù)網(wǎng)絡測試誤差表

    表4 上層4-8傳感器對不同隱含層節(jié)點數(shù)網(wǎng)絡測試誤差表

    3.3 測試和評價網(wǎng)絡

    將訓練樣本和測試樣本分別輸入上文所建立的具有最優(yōu)結(jié)構(gòu)和最佳性能的4個網(wǎng)絡模型,進行訓練和測試后得到如5所示的測試結(jié)果,誤差統(tǒng)計列于表5。4個網(wǎng)絡模型的平均相對誤差|RE|<6%,平均絕對誤差MAE<0.004,均方根誤差RMSE<13%。

    圖5 4個子網(wǎng)絡模型的輸出值與實際值比較

    表5 各個子網(wǎng)絡測試誤差統(tǒng)計

    3.4 集成網(wǎng)絡模型設(shè)計

    集成學習是指使用多個學習器對同一個問題進行學習,再使用某種策略將各個學習器進行集合,從而獲得更優(yōu)學習效果的機器學習方法。本文將以上所建立的4個子網(wǎng)絡進行集成,以期獲得性能更優(yōu)的軟測量模型。

    當集成神經(jīng)網(wǎng)絡用于回歸估計時,集成的輸出通常由各網(wǎng)絡的輸出通過簡單平均或者加權(quán)平均產(chǎn)生[25],即對各個學習器的輸出結(jié)果求平均值,進而得到最終的預測結(jié)果,但是由于采用加權(quán)平均可能會導致過度擬合,從而使集成網(wǎng)絡模型的泛化能力降低,因此本文選用簡單平均法作為集成策略。圖6中,ZTi表示通過不同傳感器對測量得到的阻抗模值,BP-Ti表示對應不同傳感器對的子網(wǎng)絡。

    圖6 集成網(wǎng)絡模型結(jié)構(gòu)示意圖

    圖7和圖8所示為集成網(wǎng)絡模型的測試結(jié)果。與獨立的子網(wǎng)絡相比(表5所示),利用集成網(wǎng)絡模型對飽和度進行預測時誤差更小。集成神經(jīng)網(wǎng)絡的平均相對誤差MRE為3.33%、平均絕對誤差MAE為0.001 4、均方根誤差RMSE為6.56%,三項誤差指標均低于各個子網(wǎng)絡的預測誤差。集成神經(jīng)網(wǎng)絡的輸入包含了通過4對傳感器(位于4個測量方位)所獲得的阻抗值,更多的信息量以及有效的集成策略對于提升預測模型的預測準確度提供了前提。

    圖7 集成網(wǎng)絡模型預測輸出

    圖8 集成網(wǎng)絡預測值相對誤差

    4 結(jié)束語

    針對沉積物中含天然氣水合物飽和度的計算問題,提出了一種基于電阻抗特性參數(shù)和集成神經(jīng)網(wǎng)絡的軟測量模型建立方法。開展了甲烷水合物模擬實驗并獲取了電阻抗譜、溫度和壓力數(shù)據(jù),在對電阻抗譜數(shù)據(jù)進行預處理、特征參數(shù)提取和選擇之后形成了樣本集,針對4對傳感器分別設(shè)計了BP神經(jīng)網(wǎng)絡,繼而將4個BP網(wǎng)絡作為子網(wǎng)絡進行集成得到集成網(wǎng)絡模型。模型測試結(jié)果表明:通過集成網(wǎng)絡模型計算得到的含水合物飽和度值平均相對誤差3.33%、平均絕對誤差0.001 4、均方根誤差為6.56%,三項誤差指標均低于各個子網(wǎng)絡的計算誤差。

    通過在寬頻范圍內(nèi)對含水合物沉積物進行電阻抗譜測試,能夠獲得沉積物的頻率響應特性以及特性描述參數(shù)(如不同頻率條件下的阻抗模值),從而能夠為神經(jīng)網(wǎng)絡模型提供大量的輸入信息。利用集成神經(jīng)網(wǎng)絡能夠綜合應用位于不同測量方位的多個傳感器的測量數(shù)據(jù),通過采用適合的集成策略能夠克服水合物空間分布不均勻?qū)︼柡投扔嬎銣蚀_度的不利影響。

    基于集成神經(jīng)網(wǎng)絡的軟測量模型能夠接受多維度、豐富的輸入信息,非常適合于解決復雜沉積物中含水合物飽和度的評價問題,下一步將在擴大飽和度范圍的基礎(chǔ)上,聯(lián)合應用聲學特征參數(shù)和電學特征參數(shù)作為模型的輸入,以期獲得更加可靠、精確的飽和度評價結(jié)果。

    猜你喜歡
    模值水合物飽和度
    糖臬之吻
    氣井用水合物自生熱解堵劑解堵效果數(shù)值模擬
    電纜參數(shù)對ZPW-2000A 軌道電路接收側(cè)傳輸性能的影響研究
    熱水吞吐開采水合物藏數(shù)值模擬研究
    天然氣水合物保壓轉(zhuǎn)移的壓力特性
    我國海域天然氣水合物試采成功
    制作一個泥土飽和度測試儀
    一種小波變換的數(shù)字信號調(diào)制方式識別方法研究
    巧用有機物的不飽和度
    柔情粉色
    優(yōu)雅(2014年4期)2014-04-18 10:18:25
    av一本久久久久| 麻豆国产av国片精品| 天天操日日干夜夜撸| 久久久久久久久久久久大奶| 少妇的丰满在线观看| 欧美另类一区| 啦啦啦视频在线资源免费观看| 91精品国产国语对白视频| 9191精品国产免费久久| 亚洲综合色网址| 校园人妻丝袜中文字幕| 一本综合久久免费| 亚洲三区欧美一区| 在线观看人妻少妇| 男女无遮挡免费网站观看| 91精品三级在线观看| 90打野战视频偷拍视频| 亚洲精品日韩在线中文字幕| 成年女人毛片免费观看观看9 | 大码成人一级视频| 如日韩欧美国产精品一区二区三区| 免费看十八禁软件| 大香蕉久久成人网| 日本午夜av视频| 欧美日韩福利视频一区二区| 亚洲国产欧美一区二区综合| www.av在线官网国产| 午夜av观看不卡| 80岁老熟妇乱子伦牲交| 啦啦啦啦在线视频资源| 久久国产亚洲av麻豆专区| 亚洲色图综合在线观看| 黄色a级毛片大全视频| 9色porny在线观看| 国产男人的电影天堂91| e午夜精品久久久久久久| 一本色道久久久久久精品综合| 欧美日韩福利视频一区二区| 亚洲国产av新网站| 少妇人妻久久综合中文| 国产免费视频播放在线视频| 日韩一区二区三区影片| 人人妻人人爽人人添夜夜欢视频| 久久99一区二区三区| 日韩人妻精品一区2区三区| av又黄又爽大尺度在线免费看| 大片免费播放器 马上看| 手机成人av网站| 国产日韩欧美在线精品| 看十八女毛片水多多多| 亚洲,欧美精品.| 久久天躁狠狠躁夜夜2o2o | 精品福利观看| 在线观看人妻少妇| 飞空精品影院首页| 精品少妇一区二区三区视频日本电影| 一级毛片我不卡| 久久久久久久精品精品| 欧美少妇被猛烈插入视频| 极品人妻少妇av视频| 男人操女人黄网站| 只有这里有精品99| 黄频高清免费视频| 大陆偷拍与自拍| 美女高潮到喷水免费观看| 另类亚洲欧美激情| 国产在线观看jvid| 久久天堂一区二区三区四区| 久久久久久久久久久久大奶| 黄色片一级片一级黄色片| 午夜激情久久久久久久| 美女福利国产在线| 日韩中文字幕欧美一区二区 | 国产成人精品久久久久久| 亚洲美女黄色视频免费看| 一级,二级,三级黄色视频| 国产一区二区三区av在线| 免费看十八禁软件| 国产99久久九九免费精品| 久久精品国产综合久久久| 丝瓜视频免费看黄片| 国产不卡av网站在线观看| a级毛片在线看网站| 人人妻人人澡人人看| 日本猛色少妇xxxxx猛交久久| 久久综合国产亚洲精品| 在线 av 中文字幕| 亚洲少妇的诱惑av| 满18在线观看网站| 女人精品久久久久毛片| 黄色毛片三级朝国网站| 日韩一区二区三区影片| 成人午夜精彩视频在线观看| 51午夜福利影视在线观看| 爱豆传媒免费全集在线观看| 久久精品国产综合久久久| 97人妻天天添夜夜摸| 久久精品aⅴ一区二区三区四区| 波野结衣二区三区在线| 观看av在线不卡| 久久99热这里只频精品6学生| 色婷婷久久久亚洲欧美| 99久久人妻综合| kizo精华| 波多野结衣一区麻豆| 日韩av不卡免费在线播放| 国产精品久久久人人做人人爽| 人体艺术视频欧美日本| 99久久人妻综合| 日韩中文字幕欧美一区二区 | 亚洲欧美一区二区三区久久| 日本vs欧美在线观看视频| 狠狠婷婷综合久久久久久88av| 国产精品九九99| 成人午夜精彩视频在线观看| 国产成人系列免费观看| 亚洲专区中文字幕在线| 一边摸一边做爽爽视频免费| 一本大道久久a久久精品| 亚洲国产欧美网| 50天的宝宝边吃奶边哭怎么回事| 日本av免费视频播放| 久久久久久久久久久久大奶| 久久久精品国产亚洲av高清涩受| 国产高清videossex| 秋霞在线观看毛片| 久热这里只有精品99| 最新的欧美精品一区二区| 精品欧美一区二区三区在线| 啦啦啦在线免费观看视频4| 国产97色在线日韩免费| 色婷婷久久久亚洲欧美| 一级毛片电影观看| av福利片在线| 肉色欧美久久久久久久蜜桃| 亚洲精品一卡2卡三卡4卡5卡 | 午夜福利乱码中文字幕| 精品人妻一区二区三区麻豆| 国产成人欧美| 一区二区三区激情视频| 老汉色∧v一级毛片| 丰满饥渴人妻一区二区三| 亚洲精品成人av观看孕妇| 美女国产高潮福利片在线看| 国产精品久久久久成人av| 日韩免费高清中文字幕av| 又黄又粗又硬又大视频| 999精品在线视频| 国产av一区二区精品久久| 国产免费福利视频在线观看| 亚洲成国产人片在线观看| 高清不卡的av网站| 自线自在国产av| 日本黄色日本黄色录像| 一区在线观看完整版| 激情五月婷婷亚洲| 侵犯人妻中文字幕一二三四区| 亚洲图色成人| 99国产综合亚洲精品| 免费人妻精品一区二区三区视频| 欧美国产精品一级二级三级| 国产成人免费无遮挡视频| 妹子高潮喷水视频| 亚洲成av片中文字幕在线观看| 国产黄色免费在线视频| 中文字幕最新亚洲高清| 国产日韩欧美视频二区| 亚洲久久久国产精品| 国产成人精品久久二区二区91| 少妇人妻 视频| 我要看黄色一级片免费的| 色94色欧美一区二区| 免费观看a级毛片全部| 国产成人a∨麻豆精品| 伊人久久大香线蕉亚洲五| 一二三四社区在线视频社区8| 国产成人av教育| 亚洲欧美日韩高清在线视频 | 波多野结衣一区麻豆| 国产人伦9x9x在线观看| 9热在线视频观看99| 久久国产精品男人的天堂亚洲| 99精国产麻豆久久婷婷| 亚洲美女黄色视频免费看| 丝袜人妻中文字幕| 亚洲美女黄色视频免费看| 99久久精品国产亚洲精品| e午夜精品久久久久久久| 欧美日韩亚洲国产一区二区在线观看 | 看免费av毛片| 国产97色在线日韩免费| 大码成人一级视频| 国产人伦9x9x在线观看| 水蜜桃什么品种好| 久久国产精品男人的天堂亚洲| 亚洲精品久久成人aⅴ小说| 欧美黄色片欧美黄色片| 国产91精品成人一区二区三区 | 男人操女人黄网站| 久久综合国产亚洲精品| 免费看av在线观看网站| 国产一区二区激情短视频 | 嫁个100分男人电影在线观看 | 新久久久久国产一级毛片| 黄色一级大片看看| 99久久精品国产亚洲精品| 久久精品国产综合久久久| 精品人妻1区二区| 久久国产精品大桥未久av| 18禁黄网站禁片午夜丰满| 啦啦啦在线免费观看视频4| 国产一级毛片在线| 亚洲欧美色中文字幕在线| 精品福利永久在线观看| 国产精品 欧美亚洲| 国产精品一区二区精品视频观看| 国产精品国产av在线观看| 国产女主播在线喷水免费视频网站| 国产精品人妻久久久影院| 我要看黄色一级片免费的| 国产91精品成人一区二区三区 | 在现免费观看毛片| 美女中出高潮动态图| 少妇裸体淫交视频免费看高清 | av线在线观看网站| 人成视频在线观看免费观看| 制服人妻中文乱码| 精品人妻一区二区三区麻豆| 免费观看人在逋| 亚洲精品乱久久久久久| 国产亚洲一区二区精品| 人人妻人人澡人人爽人人夜夜| 18在线观看网站| 黑人欧美特级aaaaaa片| 男人舔女人的私密视频| 男的添女的下面高潮视频| 国产精品免费大片| 精品亚洲乱码少妇综合久久| 国产免费福利视频在线观看| 一区二区三区四区激情视频| 国产亚洲精品第一综合不卡| 日本av免费视频播放| 美女福利国产在线| 精品人妻熟女毛片av久久网站| 操出白浆在线播放| 99九九在线精品视频| 国产91精品成人一区二区三区 | 91精品三级在线观看| 国产一区二区三区av在线| 亚洲成人手机| 国产一区有黄有色的免费视频| 大型av网站在线播放| 男女国产视频网站| 中文乱码字字幕精品一区二区三区| 久久久久精品人妻al黑| 精品久久久久久电影网| a 毛片基地| 亚洲国产精品999| 欧美日韩福利视频一区二区| 99久久综合免费| 老熟女久久久| 狠狠婷婷综合久久久久久88av| 色视频在线一区二区三区| 精品亚洲成a人片在线观看| 免费女性裸体啪啪无遮挡网站| 777米奇影视久久| 777久久人妻少妇嫩草av网站| 国产又色又爽无遮挡免| 成人黄色视频免费在线看| av在线播放精品| 欧美+亚洲+日韩+国产| 免费观看人在逋| 日韩中文字幕视频在线看片| 国产视频一区二区在线看| 亚洲精品久久久久久婷婷小说| 精品卡一卡二卡四卡免费| 男女床上黄色一级片免费看| 久久影院123| 国产日韩欧美在线精品| 18在线观看网站| 啦啦啦在线免费观看视频4| 乱人伦中国视频| 男女之事视频高清在线观看 | 人人澡人人妻人| 在线亚洲精品国产二区图片欧美| 精品免费久久久久久久清纯 | 夫妻午夜视频| 91成人精品电影| a级毛片在线看网站| 亚洲情色 制服丝袜| 黄色 视频免费看| 亚洲九九香蕉| www.熟女人妻精品国产| 中文字幕亚洲精品专区| 亚洲七黄色美女视频| 欧美日韩福利视频一区二区| 人人澡人人妻人| 国产精品麻豆人妻色哟哟久久| 国产福利在线免费观看视频| 欧美精品av麻豆av| 亚洲第一av免费看| 久久性视频一级片| 中文乱码字字幕精品一区二区三区| 久久亚洲国产成人精品v| 精品国产一区二区久久| 国产1区2区3区精品| 黄色a级毛片大全视频| 老司机在亚洲福利影院| 国产日韩欧美视频二区| 日韩电影二区| 波野结衣二区三区在线| 婷婷色综合大香蕉| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美黄色淫秽网站| 久久青草综合色| 九色亚洲精品在线播放| 亚洲 国产 在线| 欧美黄色片欧美黄色片| 国产精品久久久久久人妻精品电影 | 亚洲一码二码三码区别大吗| 国产日韩欧美亚洲二区| 国产精品亚洲av一区麻豆| 高清黄色对白视频在线免费看| 纯流量卡能插随身wifi吗| 超碰97精品在线观看| 午夜福利影视在线免费观看| 精品久久久久久电影网| 青春草视频在线免费观看| 别揉我奶头~嗯~啊~动态视频 | 成在线人永久免费视频| av一本久久久久| 考比视频在线观看| av在线app专区| 成年动漫av网址| 成人亚洲精品一区在线观看| 亚洲天堂av无毛| 国产色视频综合| 国产精品亚洲av一区麻豆| 亚洲精品av麻豆狂野| 日韩中文字幕欧美一区二区 | 国产欧美日韩一区二区三 | 女性生殖器流出的白浆| 免费观看a级毛片全部| 久久久久久久久免费视频了| 国产日韩欧美视频二区| 亚洲人成网站在线观看播放| 精品熟女少妇八av免费久了| 人人妻人人添人人爽欧美一区卜| 欧美精品亚洲一区二区| 精品亚洲乱码少妇综合久久| 大型av网站在线播放| 免费一级毛片在线播放高清视频 | 成人国产一区最新在线观看 | 大香蕉久久网| 亚洲天堂av无毛| 校园人妻丝袜中文字幕| 一边摸一边做爽爽视频免费| 赤兔流量卡办理| 日韩一本色道免费dvd| 美女扒开内裤让男人捅视频| 妹子高潮喷水视频| 久久av网站| 18在线观看网站| 中国国产av一级| 一级黄色大片毛片| 考比视频在线观看| 这个男人来自地球电影免费观看| 国产精品99久久99久久久不卡| 飞空精品影院首页| 国产精品久久久人人做人人爽| 99精国产麻豆久久婷婷| 亚洲欧洲国产日韩| 看免费av毛片| 国产无遮挡羞羞视频在线观看| av视频免费观看在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久精品成人免费网站| 婷婷丁香在线五月| 夜夜骑夜夜射夜夜干| 国产精品麻豆人妻色哟哟久久| 你懂的网址亚洲精品在线观看| 亚洲国产欧美一区二区综合| 美女扒开内裤让男人捅视频| 97人妻天天添夜夜摸| 赤兔流量卡办理| 成人亚洲欧美一区二区av| 日韩人妻精品一区2区三区| 国产成人91sexporn| 亚洲欧洲日产国产| 黑人欧美特级aaaaaa片| 黄色毛片三级朝国网站| 国产片特级美女逼逼视频| 国产高清视频在线播放一区 | 国产成人精品在线电影| 亚洲av欧美aⅴ国产| 各种免费的搞黄视频| 一区二区三区乱码不卡18| 日韩中文字幕欧美一区二区 | 免费观看a级毛片全部| 9色porny在线观看| 丝袜在线中文字幕| 久久精品人人爽人人爽视色| 制服人妻中文乱码| 亚洲欧美色中文字幕在线| 欧美精品一区二区免费开放| 欧美乱码精品一区二区三区| 欧美人与善性xxx| 国产成人欧美在线观看 | 亚洲人成网站在线观看播放| 欧美日本中文国产一区发布| 汤姆久久久久久久影院中文字幕| 欧美大码av| 九色亚洲精品在线播放| 少妇被粗大的猛进出69影院| 久久免费观看电影| 国产高清国产精品国产三级| 精品少妇久久久久久888优播| 国产精品一区二区在线不卡| 国产成人欧美在线观看 | 在线 av 中文字幕| 啦啦啦 在线观看视频| 成人国产av品久久久| 夜夜骑夜夜射夜夜干| 欧美日韩亚洲高清精品| 高清黄色对白视频在线免费看| 丁香六月欧美| 国产精品国产av在线观看| 欧美日韩视频精品一区| 成年女人毛片免费观看观看9 | 久久久久久人人人人人| 成人国产一区最新在线观看 | 日韩一卡2卡3卡4卡2021年| 妹子高潮喷水视频| 久久精品国产亚洲av涩爱| 欧美日韩福利视频一区二区| 亚洲成人手机| 亚洲欧美色中文字幕在线| 女警被强在线播放| 麻豆乱淫一区二区| 午夜福利视频在线观看免费| 国产福利在线免费观看视频| 青春草亚洲视频在线观看| 国产精品久久久久久精品古装| 国产高清不卡午夜福利| av网站在线播放免费| 又黄又粗又硬又大视频| tube8黄色片| 亚洲国产精品一区二区三区在线| 国产精品欧美亚洲77777| 秋霞在线观看毛片| 亚洲男人天堂网一区| 视频区图区小说| 亚洲精品日本国产第一区| 免费黄频网站在线观看国产| 欧美精品一区二区大全| 制服人妻中文乱码| 亚洲专区中文字幕在线| 亚洲欧美精品自产自拍| svipshipincom国产片| 在线 av 中文字幕| 婷婷色av中文字幕| 成年人免费黄色播放视频| 久9热在线精品视频| 免费在线观看日本一区| 纯流量卡能插随身wifi吗| 精品国产乱码久久久久久男人| 欧美日韩亚洲国产一区二区在线观看 | 精品国产乱码久久久久久小说| 黄色毛片三级朝国网站| 国产亚洲精品第一综合不卡| 成年美女黄网站色视频大全免费| 夫妻午夜视频| 啦啦啦啦在线视频资源| 在线观看免费视频网站a站| 亚洲人成电影观看| 我的亚洲天堂| 一级a爱视频在线免费观看| 美女午夜性视频免费| 久久久久久人人人人人| 精品免费久久久久久久清纯 | 久久热在线av| 性高湖久久久久久久久免费观看| 丰满迷人的少妇在线观看| 国产精品九九99| 久久青草综合色| 欧美人与善性xxx| 男女国产视频网站| avwww免费| 免费在线观看日本一区| 男女床上黄色一级片免费看| 黑人欧美特级aaaaaa片| 免费在线观看完整版高清| svipshipincom国产片| 国产在线一区二区三区精| 欧美日韩一级在线毛片| 免费看av在线观看网站| 制服诱惑二区| 亚洲精品国产区一区二| 欧美日韩亚洲综合一区二区三区_| 天堂8中文在线网| 在线精品无人区一区二区三| 久久久久久久久免费视频了| 午夜福利,免费看| 国产日韩一区二区三区精品不卡| 99久久人妻综合| 国产精品欧美亚洲77777| 成年动漫av网址| 国产男女内射视频| 美女午夜性视频免费| 青春草亚洲视频在线观看| 亚洲,欧美精品.| 亚洲av电影在线观看一区二区三区| 国产主播在线观看一区二区 | 国产日韩欧美在线精品| 成人国产av品久久久| 亚洲精品乱久久久久久| 女人高潮潮喷娇喘18禁视频| 日本wwww免费看| av在线app专区| 久久久精品国产亚洲av高清涩受| 男女无遮挡免费网站观看| 国产亚洲午夜精品一区二区久久| 日韩精品免费视频一区二区三区| 91麻豆av在线| 亚洲成人国产一区在线观看 | 男女无遮挡免费网站观看| 精品久久蜜臀av无| 欧美97在线视频| 2018国产大陆天天弄谢| 午夜免费观看性视频| 欧美国产精品一级二级三级| 久久久欧美国产精品| 成人影院久久| 欧美黄色片欧美黄色片| 黄色怎么调成土黄色| 中文字幕人妻丝袜一区二区| 啦啦啦中文免费视频观看日本| 啦啦啦 在线观看视频| 久久 成人 亚洲| 色94色欧美一区二区| 99九九在线精品视频| 久久中文字幕一级| 日本欧美国产在线视频| 永久免费av网站大全| 中文乱码字字幕精品一区二区三区| 成在线人永久免费视频| 国产精品一区二区在线不卡| 啦啦啦视频在线资源免费观看| 国产精品免费视频内射| 日韩欧美一区视频在线观看| 人人妻人人澡人人爽人人夜夜| 男女下面插进去视频免费观看| 美女主播在线视频| 一本久久精品| 亚洲,欧美,日韩| 丝袜美腿诱惑在线| 自拍欧美九色日韩亚洲蝌蚪91| 美女视频免费永久观看网站| 满18在线观看网站| 亚洲 欧美一区二区三区| 男女无遮挡免费网站观看| 久久99一区二区三区| 国产精品熟女久久久久浪| 欧美另类一区| 亚洲精品国产av蜜桃| 少妇人妻久久综合中文| 男女午夜视频在线观看| 可以免费在线观看a视频的电影网站| 巨乳人妻的诱惑在线观看| 制服人妻中文乱码| 成年美女黄网站色视频大全免费| 欧美日韩亚洲高清精品| 黑人猛操日本美女一级片| 成人亚洲欧美一区二区av| 性色av一级| av又黄又爽大尺度在线免费看| 黄色视频不卡| 97精品久久久久久久久久精品| 亚洲av男天堂| 久久av网站| 中文乱码字字幕精品一区二区三区| 国产极品粉嫩免费观看在线| 日韩欧美一区视频在线观看| 咕卡用的链子| 国产色视频综合| 精品视频人人做人人爽| 少妇猛男粗大的猛烈进出视频| 18在线观看网站| 亚洲欧美色中文字幕在线| 日日摸夜夜添夜夜爱| 日韩人妻精品一区2区三区| 97在线人人人人妻| 高潮久久久久久久久久久不卡| 久久99一区二区三区| 韩国精品一区二区三区| 丰满迷人的少妇在线观看| 久久99一区二区三区| 亚洲av成人不卡在线观看播放网 | 日日摸夜夜添夜夜爱| 婷婷色麻豆天堂久久| 搡老岳熟女国产| 男男h啪啪无遮挡| 国产精品成人在线| 精品一品国产午夜福利视频| 欧美日韩国产mv在线观看视频| 亚洲熟女精品中文字幕| 80岁老熟妇乱子伦牲交| 又黄又粗又硬又大视频| 久久精品久久久久久噜噜老黄| 1024香蕉在线观看| 女性被躁到高潮视频| 国产成人av教育| 人人妻人人爽人人添夜夜欢视频| a 毛片基地| 亚洲成国产人片在线观看|