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

    基于優(yōu)化fastICA盲源分離算法的地震屬性融合方法研究

    2018-11-07 01:18:00孫成禹
    石油物探 2018年5期
    關(guān)鍵詞:屬性數(shù)據(jù)均值規(guī)則

    邱 玥,孫成禹,唐 杰

    (中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島266580)

    目前地震屬性融合方法大致分為人工疊合、基于屬性顏色融合和基于屬性數(shù)據(jù)融合3類(lèi)方法。人工疊合法操作簡(jiǎn)單但精度較低,對(duì)信息的綜合處理及分析能力較差;基于屬性顏色融合法以RGB技術(shù)[1-2]為代表,充分利用了地震屬性中地質(zhì)構(gòu)造和巖性變化信息,但無(wú)法消除屬性間的冗余信息;基于屬性數(shù)據(jù)融合法包括多元線性回歸、人工神經(jīng)網(wǎng)絡(luò)、獨(dú)立分量分析等。這些方法直接對(duì)由數(shù)學(xué)統(tǒng)計(jì)和人工智能等手段提取得到的最優(yōu)屬性數(shù)據(jù)進(jìn)行融合,能夠獲得更高精度和更可靠的融合效果。

    與主成分分析(principal component analysis,PCA)方法相比,獨(dú)立分量分析(independent component analysis,ICA)算法主要針對(duì)數(shù)據(jù)的高階統(tǒng)計(jì)特性,能夠有效地去除高階統(tǒng)計(jì)特性中的冗余信息,更好地揭示數(shù)據(jù)的本質(zhì)結(jié)構(gòu)[3-4]。ZARZOSO等[5]提出了以最小互信息(minimum mutual information)為目標(biāo)函數(shù)的方法。BELL等[6]提出了名為Informax的自適應(yīng)算法,該算法通過(guò)選擇逼近源信號(hào)概率分布的最佳非線性函數(shù)來(lái)有效地恢復(fù)源信號(hào)分量。YUAN等[7]提出了極大似然目標(biāo)函數(shù)與非線性PCA方法。HYVARINEN等[8]提出了以峰度為目標(biāo)函數(shù)的一種快速獨(dú)立分量分析(fastICA)算法。普艷香[9]將fastICA算法引用到地震多屬性分析中。李小霞[10]研究了ICA算法在寬頻帶數(shù)據(jù)譜融合中的應(yīng)用。上述方法都未針對(duì)含噪地震屬性融合進(jìn)行研究。

    本文采用具有魯棒性強(qiáng)、統(tǒng)計(jì)特性明確的負(fù)熵近似[11]作為目標(biāo)函數(shù)。避免了上述方法中目標(biāo)函數(shù)復(fù)雜、非線性函數(shù)選取困難的缺點(diǎn)。目標(biāo)函數(shù)的優(yōu)化方法選用基于固定點(diǎn)迭代理論的fastICA算法,通過(guò)快速尋優(yōu)迭代找出非高斯性最大值。避開(kāi)了常規(guī)ICA梯度算法[12]需要選取參數(shù)和步長(zhǎng)的步驟。同時(shí)對(duì)初值采取降敏感性處理,選取具有5階收斂速度的迭代公式提高迭代穩(wěn)定性,減少迭代發(fā)散的可能性。引入滿足非高斯分布且實(shí)用性較強(qiáng)的貝葉斯閾值函數(shù)對(duì)含噪屬性數(shù)據(jù)進(jìn)行隨機(jī)噪聲壓制,改善了fastICA模型對(duì)含噪屬性數(shù)據(jù)去噪效果不顯著的局限性,進(jìn)一步提高了融合資料的品質(zhì)。最后用理論模型數(shù)據(jù)及實(shí)際地震資料進(jìn)行了應(yīng)用測(cè)試。

    1 基本原理

    1.1 地震屬性融合的fastICA算法

    當(dāng)輸入信號(hào)處于未知狀態(tài),且信號(hào)傳輸所處環(huán)境亦未知的情況下,僅依據(jù)儀器檢測(cè)到的混合信號(hào)估計(jì)出源信號(hào),是盲源分離技術(shù)研究的內(nèi)容[13]。針對(duì)某一目標(biāo)儲(chǔ)層,假設(shè)最終觀測(cè)得到的地震屬性是表征目標(biāo)儲(chǔ)層特征的各種巖性、物性、含流體性質(zhì)參數(shù)的線性混合,將這些儲(chǔ)層特性參數(shù)等效為輸入的未知源信號(hào),借鑒盲源分離思想,可將這些混合特征參數(shù)恢復(fù)成相互獨(dú)立的成分。

    本文采用fastICA算法進(jìn)行盲源分離。經(jīng)過(guò)該算法得到的屬性數(shù)據(jù)分量間去掉了高階相關(guān)性,實(shí)現(xiàn)各分量的統(tǒng)計(jì)獨(dú)立[14]。利用fastICA算法進(jìn)行地震多屬性融合的思想來(lái)源于多源信息融合,將由不同儀器模式或采集技術(shù)得到的二維地震屬性等效成圖像。為了加快算法處理效率,首先將根據(jù)任務(wù)目標(biāo)選擇出的地震屬性數(shù)據(jù)進(jìn)行平滑分塊,并從中隨機(jī)挑選出不少于1000的屬性塊作為fastICA算法輸入[15];然后利用優(yōu)化的fastICA基函數(shù)完成屬性數(shù)據(jù)塊從空間域至ICA域的映射,并在ICA域中利用構(gòu)造的貝葉斯閾值函數(shù)和4種效果不同的圖像融合規(guī)則對(duì)含噪聲屬性數(shù)據(jù)進(jìn)行去噪和融合;最后將經(jīng)過(guò)隨機(jī)噪聲壓制的融合地震屬性數(shù)據(jù)映射回空間域。

    1.2 融合規(guī)則

    不同的融合規(guī)則對(duì)應(yīng)于不同的數(shù)據(jù)取舍,均值融合規(guī)則側(cè)重于保持?jǐn)?shù)據(jù)背景信息,最大絕對(duì)值融合規(guī)則主要保持?jǐn)?shù)據(jù)邊界特征,加權(quán)融合規(guī)則是最大絕對(duì)值融合規(guī)則的拓展,而基于區(qū)域的融合規(guī)則綜合了最大絕對(duì)值和均值兩種融合規(guī)則。針對(duì)具體目標(biāo)要求,可以選擇不同的融合規(guī)則。

    1) 均值融合規(guī)則。

    (1)

    2) 最大絕對(duì)值融合規(guī)則。

    (2)

    3) 加權(quán)融合規(guī)則。

    對(duì)第i個(gè)屬性塊數(shù)據(jù),以其在變化域中均值作為指示因子,計(jì)算其權(quán)值:

    (3)

    則融合表達(dá)式為:

    (4)

    某些情況下,權(quán)值wi的分母可能很小,會(huì)導(dǎo)致計(jì)算的數(shù)值不穩(wěn)定,此時(shí)可對(duì)這些屬性塊采用均值融合規(guī)則或最大絕對(duì)值融合規(guī)則。

    4) 基于區(qū)域的融合規(guī)則。

    (5)

    mi=1的區(qū)域采取保持邊界特征的最大絕對(duì)值融合規(guī)則或加權(quán)融合規(guī)則;mi=0的區(qū)域采取保持背景信息為主的均值融合規(guī)則[16]。上述公式求取的是第i個(gè)屬性塊的融合結(jié)果,按此方法類(lèi)推,便能求出所有屬性塊在ICA域的融合結(jié)果。

    2 fastICA算法優(yōu)化

    fastICA算法能夠從多維數(shù)據(jù)中尋找出暗含在混合信號(hào)中彼此獨(dú)立的成分并使此種統(tǒng)計(jì)意義上的獨(dú)立性最大化[17]。假設(shè)t時(shí)刻觀測(cè)得到的混合信號(hào)為:

    (6)

    式中:X代表在時(shí)刻t經(jīng)觀測(cè)獲取的線性混合信號(hào);S代表原始信號(hào);A代表混合矩陣。目標(biāo)是根據(jù)唯一已知的X分離或估算出S[18]。通常無(wú)法得知混合系統(tǒng)的特征,故允許比例不定性及順序不定性現(xiàn)象存在。當(dāng)計(jì)算出的矩陣A存在逆矩陣時(shí),只要源信號(hào)數(shù)據(jù)滿足非高斯性的前提條件成立,便可以估算A的逆矩陣W,經(jīng)公式:

    (7)

    得到彼此之間完全獨(dú)立的估計(jì)源信號(hào)數(shù)據(jù)Y,Y是源數(shù)據(jù)S的最優(yōu)逼近。首先,需要對(duì)觀測(cè)數(shù)據(jù)進(jìn)行預(yù)處理,包括零均值化及白化處理。其中白化處理通過(guò)PCA算法[19]實(shí)現(xiàn):

    (8)

    在快速尋優(yōu)的過(guò)程中,本文選用負(fù)熵的近似作為目標(biāo)函數(shù):

    (9)

    式中:J(·)表示求負(fù)熵,即目標(biāo)函數(shù);E[·]表示求變量均值;g(·)代表非二次型函數(shù),其表達(dá)形式有很多,本方法選用g(Y)=-(1/4)Y4;Yg是高斯隨機(jī)變量。目標(biāo)是求(9)式最大值,等效于找出(7)式中具有最大非高斯性的分離矩陣W。為了提高迭代計(jì)算速度,采用具有5階收斂速度的改進(jìn)迭代公式[20-21]:

    (10)

    式中:g′(·),g″(·)分別代表g(·)的一階和二階導(dǎo)數(shù)。

    (11)

    (12)

    利用fastICA基及其逆矩陣可實(shí)現(xiàn)地震屬性數(shù)據(jù)空間域與ICA域的相互映射。但是傳統(tǒng)的ICA算法對(duì)隨機(jī)噪聲壓制效果不強(qiáng)。因此,利用貝葉斯方法構(gòu)造出滿足非高斯分布的閾值函數(shù)在ICA域中對(duì)隨機(jī)噪聲進(jìn)行去噪處理。貝葉斯閾值函數(shù)相較于常規(guī)的硬閾值等函數(shù)來(lái)說(shuō)更符合地震數(shù)據(jù)特征的分布規(guī)律,對(duì)隨機(jī)噪聲具有更好的壓制效果。

    假設(shè)有效信號(hào)S中混入了均值為0,方差為σ的高斯白噪聲N,則接收到的觀測(cè)數(shù)據(jù)Z可以表示為:

    (13)

    h(zi)=

    (14)

    綜合上述算法,基于優(yōu)化fastICA的地震屬性融合流程如圖1所示。

    圖1 基于優(yōu)化fastICA算法的地震屬性融合流程

    3 模型測(cè)試與實(shí)際資料效果分析

    3.1 優(yōu)化基函數(shù)模型測(cè)試

    選取圖2所示的4個(gè)信號(hào)作為原始信號(hào),將其通過(guò)隨機(jī)產(chǎn)生的矩陣,經(jīng)線性變換得到如圖3所示的混合信號(hào)。分別采用傳統(tǒng)fastICA算法和對(duì)初值進(jìn)行降敏感性并具有5階收斂速度的優(yōu)化后的fastICA算法進(jìn)行分離,得到的信號(hào)如圖4和圖5所示。

    計(jì)算兩種方法分離出的有效信號(hào)的波形相似度和信噪比,結(jié)果如表1所示??梢钥闯?采用優(yōu)化后的fastICA算法分離得到的信號(hào),波形相似度和信噪比都有了一定的提高。

    同時(shí),比較運(yùn)行10次后得到的平均迭代次數(shù),fastICA平均迭代次數(shù)為19.20次,方差為5.156,而優(yōu)化后的算法平均迭代次數(shù)為16.80次,方差為2.56,說(shuō)明優(yōu)化后的算法收斂性更強(qiáng),穩(wěn)定性更高。

    3.2 實(shí)際屬性數(shù)據(jù)對(duì)比分析

    3.2.1 實(shí)例一

    圖6為提取的某河道砂體的平均峰值振幅、總能量、均方根振幅和相干體4種含噪屬性,其中圖6a,圖6b和圖6c的振幅類(lèi)屬性中橢圓位置表明地層中發(fā)育古河道,且含有一些砂體信息。圖6d的相干體包絡(luò)反映了河道的邊界信息,具有前3個(gè)屬性不具備的特有信息。圖7是經(jīng)fastICA處理后得到的4種地震屬性數(shù)據(jù)。圖8是利用圖像融合中最大絕對(duì)值融合、均值融合、加權(quán)融合和基于區(qū)域融合4種不同融合規(guī)則對(duì)圖6中4種地震屬性進(jìn)行融合后的數(shù)據(jù)。圖9是加入貝葉斯閾值函數(shù)對(duì)原始各單一屬性進(jìn)行去噪處理后再利用融合規(guī)則進(jìn)行融合的結(jié)果。圖8和圖9中色標(biāo)代表經(jīng)歸一化處理后各融合規(guī)則下融合數(shù)據(jù)的相對(duì)大小。

    圖2 原始信號(hào)

    圖3 混合信號(hào)

    圖4 采用原始fastICA算法分離出的信號(hào)

    圖5 采用優(yōu)化fastICA算法分離出的信號(hào)

    表1 算法性能比較

    圖9與圖8相比,隨機(jī)噪聲得到明顯壓制,剖面的清晰度變高,有效信息得到凸顯。將融合后的地震屬性(圖9)與原始單一屬性(圖6)進(jìn)行對(duì)比,發(fā)現(xiàn)圖9 中蘊(yùn)含更多的地質(zhì)信息,既能顯示出河道的內(nèi)部變化特征,又能突出其邊界特征,類(lèi)似于給河道進(jìn)行了鑲邊處理,更能多角度地反映出目標(biāo)地質(zhì)體的響應(yīng)特征。對(duì)比4種不同的融合結(jié)果可以發(fā)現(xiàn),這種鑲邊效果在取最大絕對(duì)值融合結(jié)果(圖9a)中尤為明顯,主要是因?yàn)槿∽畲蠼^對(duì)值融合規(guī)則側(cè)重于保持邊界信息,故突出保留了河道的邊界特征,但該規(guī)則忽略了常值背景中的信息,使得融合后的地震資料中河道內(nèi)部信息不明顯。取均值融合(圖9b)則保留了河道內(nèi)部的信息,融合后的結(jié)果較為平滑,由于其邊緣信息過(guò)度平滑,所以并沒(méi)有出現(xiàn)鑲邊效果。加權(quán)融合規(guī)則(圖9c)是對(duì)上述兩種融合規(guī)則的綜合,融合后的結(jié)果較為折中,但是該算法不能考慮數(shù)據(jù)中目標(biāo)的特性,而一般數(shù)據(jù)中的目標(biāo)輪廓是不規(guī)則形狀,所以有時(shí)會(huì)引起融合圖像中背景對(duì)比度下降,造成融合圖像對(duì)比度低于原始圖像數(shù)據(jù)對(duì)比度。而基于區(qū)域的融合規(guī)則(圖9d)在融合過(guò)程中對(duì)數(shù)據(jù)中的邊界和背景信息分別采用與之相適應(yīng)的融合規(guī)則,使最終的融合結(jié)果中綜合了背景及邊界的信息,更能體現(xiàn)數(shù)據(jù)中目標(biāo)的特征。對(duì)比圖9c和圖9d黑色虛線框中的部分也可看出針對(duì)該數(shù)據(jù)基于區(qū)域的融合規(guī)則效果更好。

    圖7 經(jīng)fastICA處理后的4種地震屬性a 平均峰值振幅;b 總能量屬性;c 均方根振幅屬性;d 相干體屬性

    圖8 不同融合規(guī)則下的融合效果(實(shí)例一)a 最大絕對(duì)值融合;b 均值融合;c 加權(quán)融合;d 基于區(qū)域融合

    3.2.2 實(shí)例二

    圖10為XQ工區(qū)內(nèi)選取的總能量和相干體屬性(含噪聲),其中總能量屬性展示了該工區(qū)內(nèi)砂體空間分布情況,相干體屬性展示了該工區(qū)斷層分布情況。僅通過(guò)fastICA處理圖10中的兩種屬性而不加入融合規(guī)則時(shí)得到圖11所示結(jié)果。在傳統(tǒng)fastICA算法和加入閾值去噪函數(shù)的優(yōu)化fastICA算法下采用不同融合規(guī)則對(duì)總能量屬性和相干體屬性進(jìn)行融合,分別得到圖12和圖13所示結(jié)果。

    由圖11可看出,經(jīng)fastICA處理,能將反映地質(zhì)體本質(zhì)特征的源參數(shù)特征從屬性空間中提取出來(lái)。但是僅由fastICA處理的數(shù)據(jù),圖11a中橢圓部分缺少斷層信息,圖11b中又不能準(zhǔn)確描繪砂體邊界,所以需要加入融合規(guī)則進(jìn)一步綜合地質(zhì)信息。由圖12和圖13可以看出,加入融合規(guī)則后,不僅能同時(shí)得到砂巖與斷層的地質(zhì)信息,較為清晰地反映研究區(qū)域中砂巖與斷層的相互切割關(guān)系,而且與圖11相比砂體的范圍描述更加精確。圖13d中綠色虛線代表砂體區(qū)域,黑色虛線代表斷層分布。對(duì)比圖12與圖13可以明顯看出,圖13中包含的噪聲更少。同時(shí),為了驗(yàn)證優(yōu)化fastICA算法對(duì)隨機(jī)噪聲的壓制效果,對(duì)基于區(qū)域融合規(guī)則下的數(shù)據(jù)選用硬閾值、軟閾值函數(shù)與該優(yōu)化方法進(jìn)行對(duì)比,所得峰值信噪比參數(shù)見(jiàn)表2。

    表2 含噪屬性融合性能對(duì)比

    表2中峰值信噪比參數(shù)表明,優(yōu)化fastICA算法融合后的結(jié)果中隨機(jī)噪聲較硬閾值和軟閾值去噪方法更少,該方法對(duì)隨機(jī)噪聲的壓制效果更好,能夠獲得品質(zhì)更高的地震屬性融合結(jié)果。對(duì)比圖13中4種融合規(guī)則的融合效果,可以看出,圖13b比圖13a更注重背景信息,缺少細(xì)節(jié)信息;而圖13c和圖13d效果類(lèi)似,都是圖13a和圖13b的折中表現(xiàn),只是圖13c在斷層分布和融合數(shù)據(jù)的上部常值背景中對(duì)比度較圖13d稍差,所得結(jié)論與實(shí)例一相似。

    圖9 去噪后不同融合規(guī)則下的融合效果(實(shí)例一)a 最大絕對(duì)值融合;b 均值融合;c 加權(quán)融合;d 基于區(qū)域融合

    圖10 XQ地區(qū)總能量屬性(a)和相干體屬性(b)

    圖11 采用fastICA對(duì)XQ地區(qū)總能量屬性(a)和相干體屬性(b)進(jìn)行處理的結(jié)果

    圖12 不同融合規(guī)則下的融合效果(實(shí)例二)a 最大絕對(duì)值融合;b 均值融合;c 加權(quán)融合;d 基于區(qū)域融合

    圖13 去噪后不同融合規(guī)則下的融合效果(實(shí)例二)a 最大絕對(duì)值融合;b 均值融合;c 加權(quán)融合;d 基于區(qū)域融合

    4 結(jié)論

    本文將基于優(yōu)化fastICA盲源分離的算法引入地震屬性融合技術(shù)中,得到了以下結(jié)論。

    1) 優(yōu)化后的fastICA算法在迭代次數(shù)、波形相似度和信噪比方面都有一定的提高,而且突出了地震屬性間的高階統(tǒng)計(jì)特征,是分解混合觀測(cè)數(shù)據(jù)中所蘊(yùn)含的獨(dú)立信息的一種有效方法。

    2) 地震屬性數(shù)據(jù)資料試驗(yàn)結(jié)果表明,經(jīng)優(yōu)化fastICA算法融合后,資料中基本包含了單一屬性所蘊(yùn)含的地質(zhì)信息,更有利于突出與研究目標(biāo)相關(guān)的地質(zhì)特征,從而使研究人員能多角度、全方位地觀察目標(biāo)體的響應(yīng)特征。

    3) 優(yōu)化fastICA算法在實(shí)現(xiàn)地震屬性融合的過(guò)程中對(duì)隨機(jī)噪聲具有較好的壓制效果,能得到較高信噪比的地震屬性融合數(shù)據(jù),而不必再采用其它去噪算法對(duì)融合后的地震數(shù)據(jù)進(jìn)行去噪處理。算法快速穩(wěn)健,運(yùn)用方便,具有十分良好的應(yīng)用前景。

    4) 不同的融合規(guī)則對(duì)應(yīng)于不同的融合效果,取最大絕對(duì)值融合注重邊界信息,取均值融合注重背景信息,加權(quán)融合和基于區(qū)域融合是取最大絕對(duì)值和均值融合的折中表現(xiàn),但鑒于加權(quán)融合有時(shí)會(huì)降低融合圖像中背景對(duì)比度,故建議使用基于區(qū)域的融合規(guī)則。

    猜你喜歡
    屬性數(shù)據(jù)均值規(guī)則
    撐竿跳規(guī)則的制定
    數(shù)獨(dú)的規(guī)則和演變
    基于GIS的房產(chǎn)測(cè)繪管理信息系統(tǒng)架構(gòu)研究
    科技資訊(2019年18期)2019-09-17 11:03:28
    無(wú)源多傳感器綜合數(shù)據(jù)關(guān)聯(lián)算法研究
    屬性數(shù)據(jù)分析教學(xué)改革初探
    讓規(guī)則不規(guī)則
    Coco薇(2017年11期)2018-01-03 20:59:57
    TPP反腐敗規(guī)則對(duì)我國(guó)的啟示
    均值不等式失效時(shí)的解決方法
    均值與方差在生活中的應(yīng)用
    關(guān)于均值有界變差函數(shù)的重要不等式
    精品午夜福利视频在线观看一区| 国产亚洲精品一区二区www| 国产成人系列免费观看| 精品午夜福利视频在线观看一区| 日韩人妻精品一区2区三区| 亚洲激情在线av| 在线观看一区二区三区| 亚洲欧美日韩另类电影网站| 亚洲男人的天堂狠狠| 无人区码免费观看不卡| 午夜福利在线观看吧| 国产精品自产拍在线观看55亚洲| 国产视频一区二区在线看| 日本wwww免费看| 亚洲欧美日韩另类电影网站| 亚洲av第一区精品v没综合| 日本免费一区二区三区高清不卡 | 国产精品亚洲av一区麻豆| www.www免费av| 国产伦人伦偷精品视频| 日韩三级视频一区二区三区| 精品人妻在线不人妻| 精品国产一区二区三区四区第35| 亚洲全国av大片| 日本黄色视频三级网站网址| 亚洲中文字幕日韩| 亚洲精品一卡2卡三卡4卡5卡| 一级毛片精品| 看黄色毛片网站| 真人做人爱边吃奶动态| 国产精品一区二区免费欧美| 人人妻人人爽人人添夜夜欢视频| 欧美+亚洲+日韩+国产| 久久国产精品影院| 搡老熟女国产l中国老女人| 久久青草综合色| 黄色丝袜av网址大全| √禁漫天堂资源中文www| 欧美另类亚洲清纯唯美| 国产亚洲欧美98| bbb黄色大片| 女警被强在线播放| 精品国产一区二区久久| tocl精华| 中文字幕另类日韩欧美亚洲嫩草| 亚洲男人天堂网一区| 在线观看一区二区三区| 久久国产乱子伦精品免费另类| 男人舔女人的私密视频| 亚洲精品在线美女| 9热在线视频观看99| 欧美国产精品va在线观看不卡| 精品久久久久久久久久免费视频 | 国产成人精品久久二区二区免费| 国产在线精品亚洲第一网站| 国内久久婷婷六月综合欲色啪| 欧美黄色淫秽网站| 国产精华一区二区三区| 国产精品偷伦视频观看了| 麻豆av在线久日| 天堂中文最新版在线下载| 天堂中文最新版在线下载| 亚洲欧美一区二区三区黑人| 国产精品爽爽va在线观看网站 | 久99久视频精品免费| 91国产中文字幕| 别揉我奶头~嗯~啊~动态视频| 亚洲av成人av| 精品人妻在线不人妻| 色在线成人网| 神马国产精品三级电影在线观看 | 又黄又爽又免费观看的视频| 无限看片的www在线观看| 精品国产国语对白av| 看黄色毛片网站| 欧美久久黑人一区二区| √禁漫天堂资源中文www| 欧美不卡视频在线免费观看 | 夜夜躁狠狠躁天天躁| 99精品久久久久人妻精品| 两性夫妻黄色片| 在线观看日韩欧美| 制服诱惑二区| 大香蕉久久成人网| 久久99一区二区三区| 黄色视频,在线免费观看| 欧美激情极品国产一区二区三区| 黄频高清免费视频| 又大又爽又粗| 国产色视频综合| 久久午夜亚洲精品久久| 美女高潮到喷水免费观看| 90打野战视频偷拍视频| 欧美日本中文国产一区发布| 国产一卡二卡三卡精品| 日韩欧美一区视频在线观看| 免费在线观看完整版高清| 精品人妻1区二区| 巨乳人妻的诱惑在线观看| 啦啦啦免费观看视频1| 久久国产精品男人的天堂亚洲| av在线播放免费不卡| 国产精品香港三级国产av潘金莲| 久久人人97超碰香蕉20202| 乱人伦中国视频| 三上悠亚av全集在线观看| 国产精品综合久久久久久久免费 | 亚洲成a人片在线一区二区| 老司机午夜十八禁免费视频| 在线观看日韩欧美| 日本wwww免费看| 欧美黄色淫秽网站| 少妇的丰满在线观看| 国产激情久久老熟女| 久久影院123| 久久 成人 亚洲| 亚洲一区中文字幕在线| 国产乱人伦免费视频| 动漫黄色视频在线观看| 在线观看66精品国产| 日本黄色视频三级网站网址| 老司机靠b影院| 黄网站色视频无遮挡免费观看| 国产乱人伦免费视频| 97碰自拍视频| 国产xxxxx性猛交| 91麻豆精品激情在线观看国产 | 国产激情欧美一区二区| 人人妻人人添人人爽欧美一区卜| 欧美乱色亚洲激情| 亚洲精品在线美女| 国产男靠女视频免费网站| 波多野结衣高清无吗| 亚洲av成人不卡在线观看播放网| 成人三级黄色视频| 亚洲精品国产精品久久久不卡| 嫩草影院精品99| 国产精品成人在线| 免费在线观看黄色视频的| 亚洲精品国产色婷婷电影| 一本综合久久免费| 少妇 在线观看| 国产男靠女视频免费网站| 久久精品国产清高在天天线| 亚洲成人国产一区在线观看| 亚洲av第一区精品v没综合| 男人操女人黄网站| 欧美中文日本在线观看视频| 成人三级做爰电影| 校园春色视频在线观看| 黄片播放在线免费| 亚洲午夜理论影院| 女同久久另类99精品国产91| 伊人久久大香线蕉亚洲五| 美女大奶头视频| 中文字幕人妻丝袜一区二区| 人人妻人人爽人人添夜夜欢视频| 俄罗斯特黄特色一大片| 久久人妻熟女aⅴ| 国产欧美日韩一区二区精品| 性欧美人与动物交配| 看黄色毛片网站| av电影中文网址| 国产亚洲精品一区二区www| 麻豆成人av在线观看| 女人被躁到高潮嗷嗷叫费观| 日韩 欧美 亚洲 中文字幕| 日韩欧美三级三区| 国产有黄有色有爽视频| 大香蕉久久成人网| 国产麻豆69| 熟女少妇亚洲综合色aaa.| 最近最新中文字幕大全电影3 | 激情视频va一区二区三区| 999久久久国产精品视频| 嫩草影院精品99| 19禁男女啪啪无遮挡网站| 亚洲性夜色夜夜综合| 高清黄色对白视频在线免费看| 交换朋友夫妻互换小说| 人人妻人人澡人人看| 亚洲精品国产一区二区精华液| xxxhd国产人妻xxx| 免费日韩欧美在线观看| 999久久久精品免费观看国产| 亚洲 国产 在线| 99久久人妻综合| 日本欧美视频一区| 国产一卡二卡三卡精品| 久久久久久大精品| 午夜福利在线观看吧| 99热国产这里只有精品6| 婷婷丁香在线五月| 制服人妻中文乱码| 看黄色毛片网站| 亚洲欧美日韩另类电影网站| 午夜免费成人在线视频| 久久精品国产清高在天天线| 国产精品自产拍在线观看55亚洲| 啪啪无遮挡十八禁网站| 成年人黄色毛片网站| 日本vs欧美在线观看视频| 性欧美人与动物交配| 日韩大尺度精品在线看网址 | 成人影院久久| xxxhd国产人妻xxx| www日本在线高清视频| 日韩中文字幕欧美一区二区| 亚洲熟女毛片儿| 久久久国产欧美日韩av| 亚洲精品美女久久久久99蜜臀| 免费观看人在逋| 国产精品美女特级片免费视频播放器 | 黄色视频,在线免费观看| 人成视频在线观看免费观看| 男女下面插进去视频免费观看| 在线观看免费午夜福利视频| 国产精品亚洲av一区麻豆| 国产激情久久老熟女| www.熟女人妻精品国产| 国产在线精品亚洲第一网站| 在线观看免费高清a一片| 国产熟女xx| 国产精品国产高清国产av| 黑人巨大精品欧美一区二区蜜桃| 88av欧美| a级毛片在线看网站| 电影成人av| 热99re8久久精品国产| 自拍欧美九色日韩亚洲蝌蚪91| www.熟女人妻精品国产| 亚洲在线自拍视频| 精品国内亚洲2022精品成人| 50天的宝宝边吃奶边哭怎么回事| 亚洲av片天天在线观看| 在线观看66精品国产| 伊人久久大香线蕉亚洲五| 在线观看免费视频网站a站| 亚洲av电影在线进入| 电影成人av| 国产精品美女特级片免费视频播放器 | e午夜精品久久久久久久| 国产精品秋霞免费鲁丝片| 午夜久久久在线观看| 国产精品乱码一区二三区的特点 | 在线观看免费高清a一片| 窝窝影院91人妻| 亚洲一区中文字幕在线| ponron亚洲| 国产真人三级小视频在线观看| 一区福利在线观看| 日本免费a在线| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久国产一区二区| 亚洲av成人一区二区三| 亚洲精品在线美女| 亚洲第一av免费看| 99精品欧美一区二区三区四区| 啪啪无遮挡十八禁网站| 视频区图区小说| 久久人人爽av亚洲精品天堂| 黑人欧美特级aaaaaa片| 中文字幕av电影在线播放| 久久久久国产精品人妻aⅴ院| 色哟哟哟哟哟哟| 搡老岳熟女国产| 五月开心婷婷网| 99riav亚洲国产免费| av片东京热男人的天堂| 久久国产精品人妻蜜桃| 9191精品国产免费久久| 麻豆久久精品国产亚洲av | 一区二区三区精品91| 久久精品国产99精品国产亚洲性色 | 日韩欧美免费精品| 亚洲av美国av| 亚洲九九香蕉| 国产熟女午夜一区二区三区| 老司机福利观看| 精品国产国语对白av| 少妇 在线观看| 亚洲精品一区av在线观看| 亚洲精品美女久久av网站| 嫁个100分男人电影在线观看| 久热这里只有精品99| 日韩三级视频一区二区三区| 中文欧美无线码| 日韩欧美在线二视频| 欧美日韩av久久| 国产真人三级小视频在线观看| 亚洲人成77777在线视频| 日韩有码中文字幕| 人人澡人人妻人| 色老头精品视频在线观看| 亚洲国产精品999在线| 精品国产乱码久久久久久男人| 男人操女人黄网站| 午夜成年电影在线免费观看| 欧美乱色亚洲激情| 国产深夜福利视频在线观看| 午夜福利,免费看| 国产精品偷伦视频观看了| 国产高清国产精品国产三级| 一二三四在线观看免费中文在| 国产91精品成人一区二区三区| 国产黄色免费在线视频| 日韩一卡2卡3卡4卡2021年| 国产成人精品久久二区二区免费| av天堂在线播放| 日本精品一区二区三区蜜桃| 桃红色精品国产亚洲av| 国产精品一区二区免费欧美| 精品少妇一区二区三区视频日本电影| 亚洲第一青青草原| 国产黄色免费在线视频| 9色porny在线观看| 日本黄色视频三级网站网址| 亚洲成av片中文字幕在线观看| 香蕉国产在线看| 久久久久国产一级毛片高清牌| 亚洲午夜理论影院| 最好的美女福利视频网| 免费高清在线观看日韩| 日本vs欧美在线观看视频| 黄色毛片三级朝国网站| 大型黄色视频在线免费观看| 91精品三级在线观看| 99香蕉大伊视频| xxx96com| 亚洲 欧美 日韩 在线 免费| 男女床上黄色一级片免费看| 制服诱惑二区| 国产精品电影一区二区三区| 大型黄色视频在线免费观看| 亚洲精品国产精品久久久不卡| 亚洲人成电影免费在线| 久久久久久久久中文| 国产野战对白在线观看| 亚洲av熟女| 午夜久久久在线观看| 在线观看66精品国产| 亚洲男人天堂网一区| 国产欧美日韩一区二区精品| 午夜日韩欧美国产| 99久久人妻综合| 久久中文字幕一级| 婷婷丁香在线五月| av免费在线观看网站| 久热这里只有精品99| 最近最新中文字幕大全免费视频| 成人三级做爰电影| 手机成人av网站| 久久精品aⅴ一区二区三区四区| 69av精品久久久久久| www.精华液| 免费在线观看影片大全网站| 国产一区二区三区综合在线观看| 视频区图区小说| 国产精品久久久久成人av| 国产亚洲精品第一综合不卡| 亚洲一码二码三码区别大吗| 亚洲avbb在线观看| 亚洲欧美日韩另类电影网站| 天堂动漫精品| 久久久精品国产亚洲av高清涩受| 午夜福利影视在线免费观看| 神马国产精品三级电影在线观看 | 亚洲熟妇熟女久久| 别揉我奶头~嗯~啊~动态视频| 欧美最黄视频在线播放免费 | 老汉色av国产亚洲站长工具| 操出白浆在线播放| 国产精品野战在线观看 | 亚洲一区二区三区欧美精品| 久久久久国产一级毛片高清牌| 不卡av一区二区三区| 欧美+亚洲+日韩+国产| av在线天堂中文字幕 | 麻豆一二三区av精品| 一级黄色大片毛片| 国产精品九九99| 亚洲欧美激情综合另类| 亚洲中文日韩欧美视频| 亚洲成人免费av在线播放| 18禁美女被吸乳视频| 国产精品久久电影中文字幕| 亚洲欧美激情在线| 国产蜜桃级精品一区二区三区| 18禁美女被吸乳视频| 国产av一区在线观看免费| 激情视频va一区二区三区| 久久精品人人爽人人爽视色| 国产亚洲欧美精品永久| 亚洲精品在线美女| 国内毛片毛片毛片毛片毛片| xxx96com| 亚洲自拍偷在线| 黄片大片在线免费观看| 韩国av一区二区三区四区| 在线观看免费视频网站a站| 99久久国产精品久久久| 嫁个100分男人电影在线观看| 夜夜夜夜夜久久久久| 一级毛片精品| 久久久国产欧美日韩av| а√天堂www在线а√下载| 黄色丝袜av网址大全| 日韩精品青青久久久久久| 日韩视频一区二区在线观看| 亚洲男人的天堂狠狠| 久久婷婷成人综合色麻豆| 狠狠狠狠99中文字幕| 高潮久久久久久久久久久不卡| 无遮挡黄片免费观看| 国产成人影院久久av| 黑人操中国人逼视频| 精品国产一区二区三区四区第35| 两性夫妻黄色片| 欧美乱色亚洲激情| 国产又色又爽无遮挡免费看| 桃红色精品国产亚洲av| 韩国av一区二区三区四区| 久久中文看片网| 欧美黄色淫秽网站| 身体一侧抽搐| 黄频高清免费视频| 男女高潮啪啪啪动态图| 嫩草影视91久久| 在线观看午夜福利视频| 啦啦啦在线免费观看视频4| 午夜福利欧美成人| 香蕉国产在线看| 国产成+人综合+亚洲专区| 亚洲美女黄片视频| 中文欧美无线码| 日韩国内少妇激情av| 国产成年人精品一区二区 | 亚洲精品在线美女| 俄罗斯特黄特色一大片| 国产高清videossex| 亚洲一区二区三区色噜噜 | 大码成人一级视频| 亚洲三区欧美一区| 亚洲av日韩精品久久久久久密| 欧美成人免费av一区二区三区| 午夜两性在线视频| 黄色视频,在线免费观看| 久久影院123| 亚洲成人精品中文字幕电影 | 99在线视频只有这里精品首页| 国产伦一二天堂av在线观看| 天堂动漫精品| 欧美一级毛片孕妇| 麻豆久久精品国产亚洲av | 亚洲一区二区三区不卡视频| 天堂动漫精品| 一级片'在线观看视频| 天天躁夜夜躁狠狠躁躁| 最近最新免费中文字幕在线| 亚洲精品av麻豆狂野| 久久亚洲精品不卡| 免费在线观看亚洲国产| 日韩中文字幕欧美一区二区| 天堂中文最新版在线下载| xxxhd国产人妻xxx| 男女下面进入的视频免费午夜 | 丁香六月欧美| 午夜a级毛片| 国产97色在线日韩免费| 不卡一级毛片| 母亲3免费完整高清在线观看| 啦啦啦在线免费观看视频4| 黄色女人牲交| 欧美中文日本在线观看视频| 自线自在国产av| av网站免费在线观看视频| 国产精品久久久久成人av| 久久久水蜜桃国产精品网| 天堂影院成人在线观看| 午夜福利欧美成人| 无人区码免费观看不卡| 手机成人av网站| 啦啦啦 在线观看视频| 午夜精品久久久久久毛片777| 黄色怎么调成土黄色| 超碰97精品在线观看| 两个人免费观看高清视频| 亚洲第一青青草原| 黄色视频,在线免费观看| 一区在线观看完整版| 午夜免费鲁丝| 老司机靠b影院| 久久香蕉国产精品| √禁漫天堂资源中文www| 最好的美女福利视频网| 国产国语露脸激情在线看| 水蜜桃什么品种好| 日本a在线网址| 欧美日韩国产mv在线观看视频| 国产黄a三级三级三级人| 成人影院久久| 一区二区日韩欧美中文字幕| 老司机午夜福利在线观看视频| 国产精品影院久久| 免费在线观看黄色视频的| 成人特级黄色片久久久久久久| 交换朋友夫妻互换小说| 欧美 亚洲 国产 日韩一| 天天添夜夜摸| 91精品国产国语对白视频| 亚洲熟妇中文字幕五十中出 | 涩涩av久久男人的天堂| 久久香蕉激情| 在线观看www视频免费| 色精品久久人妻99蜜桃| 热99国产精品久久久久久7| 视频区欧美日本亚洲| 黄色怎么调成土黄色| 在线永久观看黄色视频| 久久精品国产99精品国产亚洲性色 | 亚洲一码二码三码区别大吗| 成人永久免费在线观看视频| 午夜福利在线免费观看网站| 人妻丰满熟妇av一区二区三区| 国产成年人精品一区二区 | 男人舔女人下体高潮全视频| 十分钟在线观看高清视频www| 亚洲 国产 在线| 欧美日韩福利视频一区二区| 久久婷婷成人综合色麻豆| 久久精品亚洲av国产电影网| 国产亚洲av高清不卡| av网站免费在线观看视频| 欧美激情久久久久久爽电影 | 亚洲七黄色美女视频| 老司机深夜福利视频在线观看| 亚洲九九香蕉| 欧美日本亚洲视频在线播放| 麻豆国产av国片精品| 久久久久久免费高清国产稀缺| 人成视频在线观看免费观看| 亚洲午夜精品一区,二区,三区| 久久99一区二区三区| 日韩精品中文字幕看吧| 亚洲精品中文字幕一二三四区| 不卡av一区二区三区| 日本免费a在线| 在线十欧美十亚洲十日本专区| 亚洲av熟女| 亚洲色图av天堂| 免费搜索国产男女视频| 亚洲专区国产一区二区| 久久天躁狠狠躁夜夜2o2o| 9热在线视频观看99| 免费人成视频x8x8入口观看| 欧美日韩瑟瑟在线播放| 日日干狠狠操夜夜爽| av电影中文网址| 欧美日韩中文字幕国产精品一区二区三区 | 国产在线观看jvid| 亚洲成国产人片在线观看| 村上凉子中文字幕在线| 免费在线观看亚洲国产| 国产三级在线视频| 极品人妻少妇av视频| 乱人伦中国视频| 亚洲成a人片在线一区二区| 亚洲自拍偷在线| 久久草成人影院| 美女午夜性视频免费| 成熟少妇高潮喷水视频| 一区福利在线观看| 日本三级黄在线观看| 国产真人三级小视频在线观看| av网站在线播放免费| 日本欧美视频一区| 久久久国产一区二区| 日本 av在线| 久久午夜综合久久蜜桃| 亚洲成人国产一区在线观看| 国产高清videossex| 亚洲av美国av| 国产99久久九九免费精品| 99在线人妻在线中文字幕| 窝窝影院91人妻| 一个人免费在线观看的高清视频| 黑人操中国人逼视频| 欧美日韩亚洲综合一区二区三区_| 亚洲av片天天在线观看| 国产av精品麻豆| 亚洲专区字幕在线| 狠狠狠狠99中文字幕| 曰老女人黄片| 国产精品久久久久久人妻精品电影| 亚洲一区中文字幕在线| av网站免费在线观看视频| 在线观看免费视频日本深夜| 91国产中文字幕| 亚洲国产欧美日韩在线播放| 精品熟女少妇八av免费久了| 大型黄色视频在线免费观看| 日本vs欧美在线观看视频| 久久久国产欧美日韩av| 可以免费在线观看a视频的电影网站| 啦啦啦免费观看视频1| 国产精品永久免费网站| 亚洲精品一二三| 亚洲色图综合在线观看| 五月开心婷婷网| 一本综合久久免费| 99国产精品一区二区蜜桃av| 久热爱精品视频在线9| 欧美日韩视频精品一区| 波多野结衣高清无吗| 久久久国产成人精品二区 |