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

        ?

        基于反演的穩(wěn)定高效衰減補償方法

        2014-09-25 02:16:16王本鋒陳小宏李景葉陳增保劉國昌
        地球物理學報 2014年4期
        關(guān)鍵詞:反射系數(shù)振幅增益

        王本鋒,陳小宏,李景葉,陳增保,劉國昌

        1中國石油大學(北京)油氣資源與探測國家重點實驗室,北京 102249 2中國石油大學(北京)CNPC物探重點實驗室,北京 102249

        1 引言

        通過衰減補償獲得高分辨的地震資料,有利于地層或巖性解釋,提高油氣藏檢測或儲層反演的精度,是地質(zhì)學家和地球物理學家共同關(guān)注的熱點問題之一.地震波在地下介質(zhì)中傳播時,由于地下介質(zhì)的黏滯性以及非均勻性,經(jīng)歷了大地的Q值濾波過程,使得子波的高頻成分迅速衰減,垂向分辨率降低(Hargreaves and Calvert,1991;Varela etal.,1993).很多衰減模型都可以用來描述地震波的衰減過程,Toverud和Ursin(2005)利用零偏VSP數(shù)據(jù)比較了八種不同的衰減模型,效果相當,但是最常用的衰減頻散模型為Kolsky-Futterman模型.基于該模型,Wang(2002,2006)在波場延拓的基礎(chǔ)上推導了具有大地濾波效應(yīng)的衰減補償公式以及正Q濾波公式;Wang(2011)在爆炸反射面的基礎(chǔ)上推導了正Q濾波公式;Zhang和Ulrych(2007)利用貝葉斯理論和最小二乘反演策略,實現(xiàn)了衰減補償過程,得到了高分辨率的反射系數(shù)序列.

        反Q濾波過程可看成逆波場延拓或偏移的過程,消除大地Q值濾波的影響.相位校正是無條件穩(wěn)定的過程,但振幅補償具有不穩(wěn)定性,可以用振幅增益限制或其他穩(wěn)定性方法解決.Hale(1991)指出,在反Q濾波過程中,波場的直接顯式外推方法不穩(wěn)定,能量隨著深度的增加呈指數(shù)增加;隱式外推方法一定程度上是穩(wěn)定的,但是計算量大,效率低,難以向多維推廣,因此設(shè)計穩(wěn)定的顯式外推方法成為關(guān)鍵問題.Hargreaves(1992)分析比較了一些簡單的反Q濾波方法,為利用濾波器的相似性質(zhì)實現(xiàn)反Q濾波奠定基礎(chǔ).趙建勛和倪克森(1992)借鑒串聯(lián)偏移的思想,將常Q算法推廣到深變Q模型;將串聯(lián)頻散補償算法與限幅補償算法相結(jié)合,更好地改善地震資料的分辨率.Wang(2002,2006)在波場延拓的基礎(chǔ)上,推導了反Q濾波公式,由于振幅補償?shù)牟环€(wěn)定性,提出了增益限制以及穩(wěn)定性方法;Yan和Liu(2009)將該穩(wěn)定性方法推廣到轉(zhuǎn)換波中,實現(xiàn)了疊前反射縱波以及PS轉(zhuǎn)換波的反Q濾波,但是該穩(wěn)定性方法對高頻成分沒有起到補償作用,以致在Q值較小時,不能起到很好的補償作用.余振等(2009)對反Q濾波方法進行了詳細研究,根據(jù)其實現(xiàn)形式將其分為三類:用級數(shù)展開作近似高頻補償?shù)姆碤濾波方法、基于波場延拓的反Q濾波方法以及其他的反Q濾波方法,并對其穩(wěn)定性進行了分析.Wang(2011)在爆炸反射面的基礎(chǔ)上,利用正Q濾波公式,借鑒反演的思想實現(xiàn)了地震記錄的衰減補償;忽略子波的影響進行求解,即可實現(xiàn)補償?shù)倪^程,但物理解釋困難.Zhang和Ulrych(2007)利用貝葉斯理論和最小二乘反演策略,實現(xiàn)了補償?shù)倪^程,得到了高分辨率反射系數(shù)序列,但是它只是利用時變子波與反射系數(shù)進行褶積得到地震記錄,由于此方法涉及子波的提取過程,若子波提取不準,其誤差會影響反射系數(shù)反演的精度.姚振興等(2003)研究了深度域衰減補償方法,因為深度域地震剖面是石油地震勘探的最終成果.

        反Q濾波的過程需要知道Q的信息,常規(guī)的Q值提取方法有頻譜比法、上升時間法等(Zhang and Ulrych,2007);趙偉和葛艷(2008)在小波域利用零偏移距VSP資料對Q值提取進行了分析.由于Q值求取困難,但其可利用Gabor變換,通過非穩(wěn)態(tài)反褶積提高地震記錄的分辨率,避開求取Q值的過程,直接求取反射系數(shù),得到高分辨的地震資料(Margrave etal.,2011).但是它要求子波以及衰減函數(shù)都是最小相位的,而實際資料一般為混合相位,假設(shè)偏離實際,會導致其在實際資料處理中失效.

        本文基于波場延拓正Q濾波公式,借鑒反演的思想,實現(xiàn)地震記錄的衰減補償,最終得到高分辨的地震記錄.由于反問題求解是一個不穩(wěn)定的問題,這也詮釋了衰減補償是一個不穩(wěn)定的過程;反問題的求解可通過正則化方法進行,最終得到補償后的地震記錄.該過程避免了常規(guī)反Q濾波方法振幅補償不穩(wěn)定的缺陷,改善了振幅增益限制方法以及穩(wěn)定的反Q濾波方法對Q值較小時欠補償?shù)木窒?

        2 方法原理

        在波場延拓的基礎(chǔ)上,Wang(2002)推導了反Q濾波公式,

        式中的兩個指數(shù)項分別為對振幅以及相位的校正項.振幅校正項隨著時間和頻率的增加呈指數(shù)增加,這是導致振幅補償不穩(wěn)定的原因.鑒于其不穩(wěn)定性,Wang提出了增益限制方法和穩(wěn)定性方法,限制對高頻信號的補償.設(shè)振幅補償項為

        則增益限制方法對該補償項修改為

        整理公式(1),再利用其遞推關(guān)系式得

        比較穩(wěn)定反Q濾波方法與振幅增益限制方法,對無噪數(shù)據(jù),兩種方法都能對衰減的記錄起到穩(wěn)定化補償?shù)淖饔茫赒值較小的情況下,均表現(xiàn)為補償不足;對含噪數(shù)據(jù),增益限制會出現(xiàn)不穩(wěn)定情況,穩(wěn)定的反Q濾波方法的補償效果表現(xiàn)為欠補償.

        鑒于振幅增益限制方法以及穩(wěn)定的反Q濾波方法對Q較小情況下補償不足的缺陷,本文結(jié)合Wang(2011)以及Zhang和Ulrych(2007)的思想,基于正Q濾波公式,利用反演的策略,實現(xiàn)地震記

        錄的衰減補償.具有大地吸收效應(yīng)的正Q濾波公式為(Wang,2008)

        由于地震數(shù)據(jù)的帶限性,一般在有效頻帶內(nèi)求出有效頻率成分.公式(9)可離散化寫成方程組的形式,如公式(10)所示,

        其中,S為衰減的地震記錄,為實值序列;A 為包含衰減信息的復值矩陣;U為有效頻帶內(nèi)未衰減的地震記錄的頻譜,為復數(shù)序列.將方程(10)轉(zhuǎn)化為實值方程,如公式(11)所示,

        簡記為

        其中,d=S, L= [Re(A),-Im(A)], m=[Re(U);Im(U)].由于地震數(shù)據(jù)頻帶有限以及噪聲的污染,方程(12)的求解是不適定的.為此,考慮到地震數(shù)據(jù)頻譜的光滑性,采用正則化策略以及L2范數(shù)約束,構(gòu)建目標泛函如下:

        其中,μ為正則化因子,D可以取為單位算子、一階微分算子或二階微分算子,其對應(yīng)的解方程為

        通過求解得到未衰減的地震記錄的頻譜,再經(jīng)傅里葉逆變換,得到補償后的地震記錄.若已知子波的信息,可利用反褶積技術(shù),得到高分辨的反射系數(shù)序列.

        另外,也可直接反演得到補償后的地震記錄.由公式(9)得

        將公式(15)離散化,得矩陣的形式,

        其中,u1是衰減的地震記錄,為實值序列;A是與衰減有關(guān)的復值矩陣;B是與傅里葉變換有關(guān)的復值矩陣;u2是未衰減的原始地震記錄,為實值序列.將方程(16)改寫為實值方程,

        其中,L2=Re(AB).由于地震數(shù)據(jù)頻帶有限以及噪聲的污染,方程(17)的求解是不適定的.為此,考慮到地震數(shù)據(jù)的光滑性,采用正則化策略以及L2范數(shù)約束,構(gòu)建目標泛函如下:

        其中,μ為正則化因子,D可以取為單位算子、一階微分算子或二階微分算子.其對應(yīng)的解方程為

        通過求解方程(19),可直接得到補償后的地震記錄.

        對比公式(14)與公式(19),公式(14)中待求解的參數(shù)為地震記錄的有效頻率分量,而公式(19)中待求解的參數(shù)為地震記錄時間序列,相當于全頻段的信息,計算量較大;且高頻成分易被噪聲干擾,不易處理.因此本文僅采用公式(14)對衰減問題進行分析求解,得到有效頻帶內(nèi)的頻譜,再利用傅里葉逆變換得到未衰減的地震記錄.

        3 數(shù)值算例

        首先,對稀疏反射系數(shù)進行試驗.設(shè)地下有5個反射界面,分別位于200ms、600ms、1000ms、1400ms、1800ms處,利用雷克子波(主頻20Hz,采樣率為2ms),首先合成無衰減的地震記錄;其次,在Q值分別為400、200、100、50、25的情況下合成衰減的地震記錄,如圖1所示.

        圖1 (a)無衰減地震記錄;(b)不同Q值對應(yīng)的衰減記錄Fig.1 (a)Seismic record without attenuation;(b)Attenuated record with different Qvalues

        圖2 (a)振幅增益限制方法以及(b)穩(wěn)定化反Q濾波方法補償結(jié)果Fig.2 Compensated results by the(a)gain-limited method and(b)stabilized inverse Qmethod

        由圖1b可以看出,隨著Q值的減小,波衰減越快;傳播時間越長,波形畸變越厲害.首先利用振幅增益限制方法以及穩(wěn)定化的反Q濾波方法分別對衰減地震記錄進行補償,結(jié)果如圖2a、2b所示;再利用本文方法對衰減地震記錄進行補償,D取為單位矩陣,求得阻尼最小二乘解,得到補償后的地震記錄如圖3所示,另外D也可取為一階微分算子或二階微分算子.

        圖3 本文方法補償結(jié)果Fig.3 Compensated results by the proposed method

        由圖2a、2b可以看出,在無噪聲時,增益限制方法以及穩(wěn)定的反Q濾波方法均能對衰減地震記錄進行穩(wěn)定性補償,但在傳播時間較長以及Q值較小時,均不同程度地表現(xiàn)出補償不足;比較圖2a、2b可知,穩(wěn)定的反Q濾波方法補償結(jié)果相對較好.本文方法只利用正Q濾波公式,借鑒反演的策略實現(xiàn)衰減補償?shù)玫降挠涗?,對振幅和相位進行了很好的校正,彌補了振幅增益限制方法以及穩(wěn)定的反Q濾波方法對Q值較小時欠補償?shù)娜毕?該模型結(jié)果闡釋了該方法的穩(wěn)定性好、精度高的優(yōu)點.

        為了驗證本文方法的抗噪性,對含弱隨機噪聲的地震記錄進行處理.每道噪聲的最大振幅與信號的最大振幅比為2%,由于噪聲水平較小,在圖上幾乎觀察不到.含噪記錄、振幅增益限制方法補償結(jié)果、穩(wěn)定的反Q濾波方法補償結(jié)果以及本文方法補償結(jié)果如圖4所示.

        由圖4可以看出,振幅增益控制方法在Q值較小時,振幅補償?shù)耐瑫r放大了噪聲,出現(xiàn)了數(shù)值不穩(wěn)定的情形;穩(wěn)定的反Q濾波方法對噪聲有一定的適應(yīng)性,相對振幅增益控制方法,穩(wěn)定的反Q濾波方法的實用性更強,但Q值較小時表現(xiàn)為欠補償.本文方法對含噪地震記錄能進行較好的補償,補償結(jié)果優(yōu)于振幅增益控制方法以及穩(wěn)定的反Q濾波方法.其關(guān)鍵在于正則化因子μ的選擇,噪聲水平較大時,選擇較大的μ,會提高穩(wěn)定性,但是會降低補償效果,反之亦然.實際資料處理中,可通過反復試驗,最終確定μ的取值.

        圖4 (a)含噪衰減記錄;(b)增益控制補償方法;(c)穩(wěn)定的反Q濾波方法;(d)本文方法補償后的記錄Fig.4 (a)Noisy seismic records;(b)Compensated results by the gain-limited method;(c)Compensated results by the stabilized inverse Qmethod;(d)Compensated results by the proposed method

        上述討論的模型,反射系數(shù)為稀疏的,偏離實際情況,下面對隨機產(chǎn)生的偽反射系數(shù)序列進行試驗,一定程度上更接近實際情況.圖5是反射系數(shù)序列以及利用雷克子波(主頻30Hz,時間采樣率為2ms)合成的無衰減的地震記錄.

        當Q值分別取400、200、100、50、25時,合成衰減地震記錄如圖6a所示,穩(wěn)定的反Q濾波方法以及本文方法對衰減地震記錄進行補償,得到補償后的記錄如圖6b—6c所示.分析圖6可知,隨著Q值的減少以及傳播時間的增大,衰減越嚴重,且本文方法的補償結(jié)果優(yōu)于穩(wěn)定的反Q濾波方法.對比無衰減的地震記錄圖5b以及穩(wěn)定的反Q濾波方法圖6b和本文方法圖6c的補償結(jié)果可知,本文方法的精度最高,基本消除了衰減的影響;克服了傳統(tǒng)反Q濾波方法不穩(wěn)定的缺陷,彌補了Q值較低時,穩(wěn)定的反Q濾波方法補償不足的局限性,其穩(wěn)定性好、精度高、補償效果好.

        圖5 (a)隨機反射系數(shù)序列;(b)無衰減的地震記錄Fig.5 (a)Random reflectivity series;(b)Seismic records without attenuation

        最后,將本文方法應(yīng)用于實際資料處理中,衰減補償前的地震記錄如圖7a所示,從圖上可以看出,地震記錄的振幅發(fā)生衰減,相位發(fā)生畸變.利用穩(wěn)定的反Q濾波方法以及本文方法對該記錄進行補償;基于常Q模型(Q=100)得到衰減補償后的地震記錄如圖7b—7c所示,分別為穩(wěn)定的反Q濾波方法以及本文方法的補償結(jié)果.通過比較圖7a—7c知,地震記錄的分辨率得到提高(同相軸變細且連續(xù)),振幅得到補償、相位得到校正;特別是黑色圓圈以及箭頭所指位置,同相軸更加清晰且連續(xù),分辨率得到提高.圖7a—7c對應(yīng)的平均振幅譜如圖8所示,通過分析振幅譜曲線可知,補償后的地震記錄的振幅譜得到展寬,本文方法以及穩(wěn)定的反Q濾波方法均能達到提高分辨率的目的,且本文方法的補償效果更好.但補償后的歸一化振幅譜曲線低頻相對變?nèi)?,分析原因如下?/p>

        (1)由公式(9)可知,高頻相對低頻振幅值衰減較快,衰減補償前,峰值頻率較低,低頻為優(yōu)勢頻率;衰減補償后,相位得到校正、振幅得到補償,峰值頻率向中高頻移動、中高頻相對低頻振幅補償較多,且高低頻放大不成比例;

        (2)實際Q值隨頻率變化,但本文采用常Q模型對衰減地震記錄進行補償,Q值不精確也是導致歸一化的振幅譜低頻相對變?nèi)醯脑蛑?相對直接求解補償后的地震記錄而言,該方法只需求解有效頻帶內(nèi)的頻率分量,求解方程組的規(guī)模相對較小,計算效率較高,具體分析詳見討論部分.

        4 討論

        本文方法具有較好的穩(wěn)定性以及較高的精度,而且只計算有效頻帶內(nèi)的頻率分量,與直接計算未衰減的地震記錄相比具有較高的效率.分析如下:設(shè)N為地震記錄的長度,直接求解地震記錄方法,要求解的方程組的規(guī)模為N×N,本文方法只需求解規(guī)模為2 M×2 M的方程組,其中M是有效頻帶內(nèi)的頻點數(shù),一般情況下M<N/4.例如,dt=0.002s時,則奈奎斯特頻率fn=250Hz,假設(shè)地震數(shù)據(jù)的最高頻率為fmax=120Hz,則.用高斯消去法解方程組時,其計算量為,其中n為方程組的規(guī)模,從該計算量上可以看出,本文方法具有較高的效率.

        Wang(2011)基于反演的衰減補償,其求解的參數(shù)為反射系數(shù)序列,由上面的計算量分析可知,其計算量較大;其忽略子波的影響,即能起到補償?shù)淖饔茫锢沓梢蚪忉岆y.Zhang和Ulrych(2007)利用最小二乘的觀點求解反射系數(shù)序列,計算量與Wang的方法相當,但是其正演記錄為衰減子波與反射系數(shù)的褶積,它基于褶積機制,而不是接近實際的波動方程機制;其次,子波估計不準會對反射系數(shù)反演帶來誤差,影響反演的精度.本文的方法不涉及子波的問題,只起到衰減補償?shù)淖饔?,來提高地震記錄的分辨率,穩(wěn)定性好、精度高,且具有較高的計算效率.

        圖6 (a)不同Q值對應(yīng)的衰減地震記錄;(b)穩(wěn)定的反Q濾波方法;(c)本文方法衰減補償結(jié)果Fig.6 (a)Attenuated seismic records with different Qvalues;(b)Compensated records by the stabilized inverse Qmethod;(c)Compensated records by the proposed method

        本文的模擬實驗以及實際數(shù)據(jù)處理都是在Q值已知的情況下進行的,且Q為等效Q值.但是其可推廣到層Q以及連續(xù)變化的Q情形,再利用本文的方法進行求解.但是,Q值的精確提取是一個比較困難的工作,Q值估計誤差會影響本文方法的處理效果,因此高精度的Q值提取方法,如基于全波形反演的Q值估計成為熱點;再者,不直接求取Q值而實現(xiàn)提高分辨率目的的非穩(wěn)態(tài)反褶積,也逐漸成為了研究趨勢.

        5 結(jié)論

        圖7 (a)補償前的地震記錄;(b)穩(wěn)定的反Q濾波方法;(c)本文方法補償后的地震記錄Fig.7 (a)Seismic records before compensation;(b)Compensated results by the stabilized inverse Qmethod;(c)Compensated records by the proposed method

        本文基于大地濾波的正Q濾波機制,借鑒反演的思想討論了一種新的衰減補償方法,可以得到補償后的地震數(shù)據(jù).該方法克服了傳統(tǒng)反Q濾波方法不穩(wěn)定的弊端,改善了振幅增益限制方法以及穩(wěn)定的反Q濾波方法補償不足的局限,具有穩(wěn)定性好、精度高的優(yōu)點;與其他直接反演地震記錄時間序列的方法相比,其只需求解有效頻帶內(nèi)的頻率成分,具有較高的效率;且該方法可以推廣到層Q以及Q值連續(xù)變化的情況,對振幅和相位同時進行校正.模擬實驗以及實際數(shù)據(jù)處理效果均驗證了本文方法的有效性.

        圖8 地震剖面的平均振幅譜對比黑線對應(yīng)圖7a,藍線對應(yīng)圖7b,紅線對應(yīng)圖7c.Fig.8 Comparison of average amplitude spectra of seismic data Black line corresponds to Fig.7a;Blue line to Fig.7b;Red line to Fig.7c.

        Hale D.1991.Stable explicit depth extrapolation of seismic wavefields.Geophysics,56(11):1770-1777,doi:10.1190/1.1442989.

        Hargreaves N,Calvert A.1991.Inverse Qfiltering by Fourier transform.Geophysics,56(4):519-527,doi:10.1190/1.1443067.

        Hargreaves N D.1992.Similarity and the inverse Qfilter:some simple algorithms for inverse Qfiltering.Geophysics,57(7):944-947,doi:10.1190/1.1443307.

        Margrave G F,Lamoureux M P,Henley D C.2011.Gabor deconvolution:Estimating reflectivity by nonstationary deconvolution of seismic data.Geophysics,76(3):W15-W30,doi:10.1190/1.3560167.

        Toverud T,Ursin B R.2005.Comparison of seismic attenuation models using zero-offset vertical seismic profiling(VSP)data.Geophysics,70(2):F17-F25,doi:10.1190/1.1884827.

        Varela C L,Rosa A L R,Ulrych T J.1993.Modeling of attenuation and dispersion.Geophysics,58(8):1167-1173,doi:10.1190/1.1443500.

        Wang S D.2011.Attenuation compensation method based on inversion.Applied Geophysics,8(2):150-157,doi:10.1007/s11770-011-0275-3.

        Wang Y H.2002.A stable and efficient approach of inverse Q filtering.Geophysics,67(2):657-663,doi:10.1190/1.1468627.

        Wang Y H.2006.Inverse Q-filter for seismic resolution enhancement.Geophysics,71(3):V51-V60,doi:10.1190/1.2192912.

        Wang Y H.2008.Seismic Inverse QFiltering.Wiley-Blackwell.

        Yan H Y,Liu Y.2009.Estimation of Qand inverse Qfiltering for prestack reflected PP-and converted PS-waves.Applied Geophysics,6(1):59-69,doi:10.1007/s11770-009-0009-y.

        Yao Z X,Gao X,Li W X.2003.The forward Q method for compensating attenuation and frequency dispersion used in the seismic profile of depth domain.Chinese J.Geophys.(in Chinese),46(2):229-233.

        Yu Z,Wang Y C,He J,etal.2009.A review of inverse Qfiltering methods.Progress in Exploration Geophysics(in Chinese),32(5):309-314,325.

        Zhang C J,Ulrych T J.2007.Seismic absorption compensation:A least squares inverse scheme.Geophysics,72(6):R109-R114,doi:10.1190/1.2766467.

        Zhao J X,Ni K S.1992.Cascaded inverse Qfiltering and the application.Oil Geophys.Prosp.(in Chinese),27(6):722-730.

        Zhao W,Ge Y.2008.Estimation of Qfrom VSP data with zero offset in wavelet domain.Chinese J.Geophys.(in Chinese),51(4):1202-1208.

        附中文參考文獻

        姚振興,高星,李維新.2003.用于深度域地震剖面衰減與頻散補償?shù)姆碤濾波方法.地球物理學報,46(2):229-233.

        余振,王彥春,何靜等.2009.反Q濾波方法研究綜述.勘探地球物理進展,32(5):309-314,325.

        趙建勛,倪克森.1992.串聯(lián)反Q濾波及其應(yīng)用.石油地球物理勘探,27(6):722-730.

        趙偉,葛艷.2008.利用零偏移距VSP資料在小波域計算介質(zhì)Q值.地球物理學報,51(4):1202-1208.

        猜你喜歡
        反射系數(shù)振幅增益
        基于增益調(diào)度與光滑切換的傾轉(zhuǎn)旋翼機最優(yōu)控制
        多道隨機稀疏反射系數(shù)反演
        石油物探(2020年6期)2020-11-25 02:38:46
        基于單片機的程控增益放大器設(shè)計
        電子制作(2019年19期)2019-11-23 08:41:36
        基于Multisim10和AD603的程控增益放大器仿真研究
        電子制作(2018年19期)2018-11-14 02:37:02
        十大漲跌幅、換手、振幅、資金流向
        十大漲跌幅、換手、振幅、資金流向
        球面波PP反射系數(shù)的頻變特征研究
        十大漲跌幅、換手、振幅、資金流向
        滬市十大振幅
        沙質(zhì)沉積物反射系數(shù)的寬帶測量方法
        久久精品国产亚洲av热九九热| 天堂aⅴ无码一区二区三区| 人人妻人人澡av天堂香蕉| 亚洲va欧美va人人爽夜夜嗨| 色偷偷亚洲精品一区二区| 爽爽影院免费观看| 国产成人精品一区二区不卡| 无码一级视频在线| 国产喷白浆精品一区二区豆腐 | 国产精品久久久久久2021| 久久婷婷夜色精品国产| 久久久国产精品123| 亚洲精品乱码久久久久久蜜桃不卡 | 欧美粗大无套gay| 久久久久人妻精品一区5555| 综合激情五月三开心五月| 婷婷五月深深久久精品| 性一乱一搞一交一伦一性| 国产成人无精品久久久| 亚洲成人av一区免费看| 久久久久人妻精品一区三寸| 国产午夜福利短视频| 久久久婷婷综合五月天| 日本五十路人妻在线一区二区| 无码成人一区二区| 色丁香在线观看| 亚洲综合色视频在线免费观看| 完整版免费av片| 亚洲啪啪综合av一区| 日本精品视频一视频高清| 亚洲精品岛国av一区二区| 国内成+人 亚洲+欧美+综合在线| 国产人成精品免费视频| 美女草逼视频免费播放| 国产精品久久久久久久久久红粉| 国产成人精品日本亚洲| 国产粉嫩高清| 一区二区三区国产黄色| 男女爽爽无遮挡午夜视频| 婷婷激情六月| 亚洲国产日韩一区二区三区四区 |