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

        ?

        基于地震沉積學(xué)的沉積相精細(xì)刻畫(huà):以山西寺家莊礦15號(hào)煤層聚煤前后為例

        2021-11-02 07:40:14常鎖亮王瑞瑞劉最亮楊智華
        中國(guó)煤炭地質(zhì) 2021年10期
        關(guān)鍵詞:沉積相振幅切片

        張 剛,常鎖亮,王瑞瑞,劉 波,曾 博,劉最亮,楊智華

        (1.中國(guó)礦業(yè)大學(xué)煤層氣資源與成藏過(guò)程教育部重點(diǎn)實(shí)驗(yàn)室,江蘇徐州 221008;2.中國(guó)礦業(yè)大學(xué) 資源與地球科學(xué)學(xué)院,江蘇徐州 221008;3.太原理工大學(xué) 地球科學(xué)與工程系,太原 030024;4.煤與煤系氣地質(zhì)山西省重點(diǎn)實(shí)驗(yàn)室,太原 030024;5.陽(yáng)泉煤業(yè)(集團(tuán))有限責(zé)任公司,山西陽(yáng)泉 045000)

        在聚煤盆地的演化過(guò)程中,不同沉積相的展布特征對(duì)聚煤盆地的海陸變遷、構(gòu)造演化及周?chē)h(huán)境的相互作用起到了很好的記錄作用,對(duì)聚煤盆地沉積相的研究成為研究聚煤盆地演化過(guò)程的核心內(nèi)容[1]。在海陸交互相的古地理背景下,晚古生代沁水盆地的沉積環(huán)境成為影響聚煤特征以及煤層氣賦存的主要因素。而聚煤前后的沉積環(huán)境又影響著煤層厚度及其空間展布特征和煤層氣的保存[2]及煤體結(jié)構(gòu)的分布。對(duì)聚煤前后沉積相的刻畫(huà),可以了解儲(chǔ)層物性結(jié)構(gòu)及砂體空間展布,進(jìn)而研究煤系氣的分布規(guī)律及形成條件,對(duì)后續(xù)煤系氣的勘探開(kāi)發(fā)起到指導(dǎo)作用,因此對(duì)含煤巖系聚煤前后沉積相的研究是十分有必要的[3]。

        前人對(duì)沁水盆地太原組地層的沉積環(huán)境已有很多研究,然而以往大量的研究多集中于沁水盆地西北部及東南部,對(duì)沁水盆地東北部沉積環(huán)境的研究相對(duì)缺少,盡管已有研究表明沁水盆地東北部太原組地層下段沉積環(huán)境為障壁島—潮坪—潟湖—潮坪相沉積[4],但前人的研究仍存在以下兩個(gè)方面的不足,一方面以往研究較籠統(tǒng),缺少對(duì)太原組15號(hào)煤層聚煤前后沉積相的精細(xì)刻畫(huà);另一方面,沁水盆地太原組地層15號(hào)煤層位于太原組地層下段,屬海陸交互相沉積。在海陸交互相的古地理背景下,研究區(qū)地層的巖性、巖相在橫向上變化較快[5]。而傳統(tǒng)的基于地質(zhì)和測(cè)井的研究方法利用的是測(cè)井在垂向上的高分辨率對(duì)地層中的沉積序列進(jìn)行研究,在橫向分辨率上難以滿(mǎn)足研究區(qū)沉積相的精細(xì)刻畫(huà)[6]。

        本文基于地震沉積學(xué)的理論方法,選擇寺家莊礦區(qū)太原組地層部分鉆孔及地震資料進(jìn)行研究,充分利用測(cè)井在垂向上的高分辨率和地震勘探在橫向上的高分辨率的優(yōu)勢(shì),將二者的優(yōu)勢(shì)相結(jié)合[7],采用合成記錄、井震聯(lián)合標(biāo)定、層位追蹤、地層切片、屬性提取等關(guān)鍵技術(shù),提取層間屬性,并通過(guò)融合屬性來(lái)預(yù)測(cè)目標(biāo)層段的砂地比[8],結(jié)合研究區(qū)的古地理背景和沉積相模式,對(duì)15號(hào)煤層聚煤前后沉積相進(jìn)行精細(xì)刻畫(huà),旨在為進(jìn)一步開(kāi)展寺家莊礦區(qū)煤及煤層氣資源開(kāi)發(fā)地質(zhì)條件研究提供重要的沉積因素支撐。

        1 地質(zhì)背景

        1.1 地層特征

        研究區(qū)位于沁水盆地陽(yáng)泉礦區(qū)(圖1),以砂巖,砂質(zhì)泥巖,泥巖為主,還發(fā)育三層較厚的灰?guī)r,可作為K2,K3,K4的標(biāo)志層。除此以外,煤層出現(xiàn)有十層以上,但大多為薄煤層,在太原組底部和中部都發(fā)育較厚的煤層,如15號(hào)煤層和9號(hào)煤層[9]。依據(jù)巖性特征可以將太原組分為三段:K1砂巖底到K2灰?guī)r底為太原組地層下段,根據(jù)測(cè)井資料可知,該段最小厚度為 34.95 m,最大厚度為57.9 m,主要為泥巖、砂質(zhì)泥巖、粉砂巖以及煤層,在底部發(fā)育15號(hào)煤層,厚度可達(dá)4~7 m,且屬于可開(kāi)采煤層,同時(shí)下段中也間夾1~2層不穩(wěn)定的薄煤層[10]。K1砂巖底到15號(hào)煤層頂屬于聚煤前沉積地層,15號(hào)煤層頂?shù)終2灰?guī)r底屬于聚煤后沉積地層[11]。K2灰?guī)r底到K4灰?guī)r底為太原組中段,厚度在34.25~51.23 m,含有較厚的石灰?guī)r和粉砂巖,還含有幾層薄煤層,其中包括11號(hào)、12號(hào)、13號(hào)煤層。K4灰?guī)r底部到K7砂巖底部為太原組上段地層,厚度在38.55~69.32 m,砂質(zhì)泥巖,砂巖,泥巖為主,夾有較厚煤層,其中9號(hào)煤層屬于厚煤層,厚度可達(dá)0.89~3.35m,同時(shí)太原組上段地層中也含有1~3層不穩(wěn)定的薄煤層。

        1.2 構(gòu)造特征

        陽(yáng)泉寺家莊礦區(qū)位于沁水盆地拗陷的東邊緣,太行山的隆起西翼,其基本構(gòu)造形態(tài)呈單斜,走向?yàn)楸北蔽?,傾向?yàn)槟衔魑鱗12]。在該單斜上又發(fā)育有次級(jí)褶曲,表現(xiàn)為“S”型[13]。地層傾角比較平穩(wěn),普遍為10°,部分受到褶曲的影響達(dá)到20°。寺家莊礦區(qū)斷層構(gòu)造較其他井田多,主要分布于東部地區(qū)[14]。

        (a)寺家莊礦區(qū)構(gòu)造分布 (b)寺家莊礦區(qū)井位分布圖1 寺家莊礦區(qū)圖Figure 1 Map of Sijiazhuang mining area

        2 研究區(qū)15號(hào)煤層聚煤前后的砂地比預(yù)測(cè)及沉積相的刻畫(huà)

        2.1 測(cè)井資料解釋

        本文選取研究區(qū)內(nèi)8口井進(jìn)行單井相解釋?zhuān)趯?duì)測(cè)井資料進(jìn)行標(biāo)準(zhǔn)化處理以后,便可根據(jù)測(cè)井曲線(xiàn)特征標(biāo)定巖性、標(biāo)志層、沉積相以及沉積亞相。本文主要根據(jù)自然伽馬和深側(cè)向電阻率曲線(xiàn)特征對(duì)巖性進(jìn)行標(biāo)定,在對(duì)標(biāo)志層進(jìn)行標(biāo)定時(shí)主要根據(jù)測(cè)井曲線(xiàn)特征相似的、相對(duì)穩(wěn)定的層段進(jìn)行劃分,此外,還可根據(jù)標(biāo)志層的特殊巖性對(duì)標(biāo)志層進(jìn)行劃分與標(biāo)定,例如K2、K3均為較厚的灰?guī)r,K1為較厚的砂巖[15]。依據(jù)研究區(qū)地質(zhì)背景以及沉積環(huán)境對(duì)沉積相進(jìn)行劃分,沁水盆地東北部太原組地層下段自下而上的沉積環(huán)境依次為障壁島—潮坪—潟湖—潮坪沉積。其中障壁島相沉積主要由厚層細(xì)粒和中粗粒砂巖組成,向上變?yōu)槟嗵科撼练e,由泥質(zhì)粉砂巖、泥巖和16號(hào)薄煤層組成。隨后海平面上升,形成了由砂質(zhì)泥巖和灰?guī)r所組成的潟湖相沉積。隨后水體變淺,形成較厚的混合坪和泥炭坪,在此基礎(chǔ)上形成了較厚的15號(hào)煤層,基于此對(duì)研究區(qū)8口井分別進(jìn)行單井相解釋?zhuān)瑸楹罄m(xù)地震資料的解釋工作打下基礎(chǔ)(圖2)。

        圖2 研究區(qū)2號(hào)井巖性及沉積相劃分Figure 2 Study area well No.2 lithology and sedimentary facies partitioning

        2.2 井震結(jié)合標(biāo)定層序界面

        本次研究在對(duì)地震資料進(jìn)行解釋之前,首先要確定反射波所對(duì)應(yīng)的標(biāo)志層,其中合成記錄是聯(lián)系地震資料與測(cè)井資料之間的橋梁,通過(guò)制作合成記錄并與時(shí)間剖面相比對(duì),確定其對(duì)應(yīng)關(guān)系,從而賦予地震反射波相應(yīng)的地質(zhì)屬性,便可在時(shí)間剖面上對(duì)地震反射波進(jìn)行層位標(biāo)定與追蹤[16],具體操作流程如下:

        1)確定含煤巖系沉積等時(shí)界面,其方法主要有定性法和定量法,本文采用定量法進(jìn)行確定,即如果某一同相軸不隨地震頻率變化而變化,那么該同相軸就可反映等時(shí)信息。通過(guò)對(duì)該地區(qū)地震反射波頻譜特征進(jìn)行分析,得出該地區(qū)地震反射波的主頻率約為50 Hz,頻帶寬度約為10 Hz~70 Hz。因此,選擇30 Hz和70 Hz進(jìn)行分頻標(biāo)定。

        2)經(jīng)過(guò)井震分頻標(biāo)定的地震標(biāo)志層K1(K1砂巖)、T15(太原組15號(hào)煤)以及K2(K2灰?guī)r)在低頻、中頻和高頻剖面上均穩(wěn)定存在。這些層位具有一定的等時(shí)地層意義,可作為建立等時(shí)地層格架的參考標(biāo)志層[17]。

        3)結(jié)合研究區(qū)制作的合成記錄與測(cè)井曲線(xiàn)特征,在地震時(shí)間剖面中追蹤了3個(gè)連續(xù)性較好的地震層位(圖3),自下而上分別是K1、T15、K2。

        圖3 地震地質(zhì)等時(shí)標(biāo)志層Figure 3 Seismic-geologic isochronous marker beds

        2.3 地震子波相位調(diào)整及去強(qiáng)反射

        我們利用90°相位轉(zhuǎn)換技術(shù)處理研究區(qū)地震數(shù)據(jù)體,使得反射波同相軸與薄層的相關(guān)關(guān)系得到建立,即薄層的中心部位可由同相軸指示。這時(shí)巖性地層的意義便被賦予到反射波同相軸之上,降低了在薄層解釋上的困難[18]。15號(hào)煤層層位在零相位地震剖面中為過(guò)零點(diǎn)追蹤界面,剖面形態(tài)不能清晰反映煤層巖性特征,而在90°相位地震剖面中T15層位被波谷直接表示,更加清楚地反映了煤的地震剖面特征。

        煤層與砂巖層呈整合接觸狀態(tài)且二者波阻抗差異較大,因此砂巖的地震反射波特征會(huì)被淹沒(méi)在煤層的強(qiáng)反射波中,這就加大了對(duì)研究區(qū)砂地比預(yù)測(cè)的難度[19]。所以本文對(duì)地震反射波進(jìn)行去強(qiáng)反射調(diào)整,實(shí)現(xiàn)在煤層強(qiáng)反射能量背景中獲取砂巖的弱響應(yīng)特征信息。

        去強(qiáng)反射前,太原組15號(hào)煤層對(duì)應(yīng)的是波谷位置且振幅較強(qiáng)。經(jīng)去強(qiáng)反射后,太原組15號(hào)煤層所在位置的振幅值明顯降低。隨著匹配次數(shù)的增加,所選時(shí)窗內(nèi)強(qiáng)反射振幅值逐步降低,形成新的波組特征(圖4)。

        2.4 基于地層切片技術(shù)的屬性提取與優(yōu)選

        地層切片技術(shù)具有等比例切片的特點(diǎn),即通過(guò)等比例內(nèi)插的手段,在選定的兩個(gè)層面內(nèi)插入切片,從而保證了切片的實(shí)用性,所以選取地層切片技術(shù)進(jìn)行切片的制作[20]。地層切片對(duì)于傳統(tǒng)的切片技術(shù)(時(shí)間切片和層位切片)而言,最大的特點(diǎn)就是可以在兩個(gè)選定界面之間等比例內(nèi)插出一系列切片,并根據(jù)內(nèi)插切片生成不同的地震屬性切片。

        (a)去強(qiáng)反射前地震時(shí)間剖面 (b)去強(qiáng)反射后地震時(shí)間剖面圖圖4 去強(qiáng)反射前后的地震時(shí)間剖面對(duì)比Figure 4 Comparison of seismic time sections before and after strong reflections removal

        地層切片技術(shù)可分三步實(shí)現(xiàn):第一,選擇等時(shí)地震的參考層,即井震聯(lián)合追蹤出的標(biāo)志層K2、T15和K1[21];第二,建立地層時(shí)間模型。由于切片是等比例內(nèi)插生成,故所有切片所包含的時(shí)間間隔一般是相同的,通過(guò)內(nèi)插線(xiàn)性函數(shù),可以創(chuàng)建出接近實(shí)際沉積地層的切片模型。第三,通過(guò)所建立的切片模型來(lái)提取相應(yīng)的地震屬性,從而形成對(duì)應(yīng)地震屬性切片。

        前人實(shí)踐表明:和振幅有關(guān)的地震屬性可以反映具波阻抗差異的沉積扇體;弧長(zhǎng)是反映反射關(guān)系的橫向變化,可以反映砂巖或泥巖地層;能量半衰時(shí)則可以反映地層中巖性巖相在橫向上的變化[22]。因此對(duì)地層切片進(jìn)行地震敏感屬性提取時(shí),根據(jù)前人總結(jié)的經(jīng)驗(yàn),分別對(duì)15號(hào)煤層頂、底板中的地層切片提取均方根振幅,能量半衰時(shí)、弧長(zhǎng)、高階累量方差、平均振幅、最大振幅6種地震屬性,表1、表2分別為15號(hào)煤層底板和頂板所對(duì)應(yīng)的地震屬性值。

        表1 15號(hào)煤層底板所對(duì)應(yīng)的屬性值Table 1 Attribute value corresponding to the bottom of No.15 coal seam

        表2 15號(hào)煤層頂板所對(duì)應(yīng)的屬性值Table 2 Coal No.15 roof corresponding attribute values

        在對(duì)單一地震屬性和各層序砂地比進(jìn)行交匯分析之前,要對(duì)每一個(gè)地震屬性先進(jìn)行標(biāo)準(zhǔn)化處理[23]。本次采用極差標(biāo)準(zhǔn)化的計(jì)算方法,設(shè)本次獲取的屬性值為xi:

        (1)

        表3 15號(hào)煤層底板對(duì)應(yīng)標(biāo)準(zhǔn)化后的屬性值Table 3 Coal No.15 floor corresponding normalized attribute values

        表4 15號(hào)煤層頂板對(duì)應(yīng)標(biāo)準(zhǔn)化后的屬性值Table 4 Coal No.15 roof corresponding normalized attribute values

        本文在對(duì)地震屬性進(jìn)行極差標(biāo)準(zhǔn)化后,建立每口井中的煤層底板砂地比與單一地震屬性的擬合關(guān)系圖(圖5)以及每口井中的煤層頂板砂地比與單一地震屬性的擬合關(guān)系圖(圖6)。

        各井中煤層頂?shù)装迳暗乇葦?shù)據(jù)與單一地震屬性經(jīng)交匯分析后發(fā)現(xiàn)相關(guān)性系數(shù)(R2)絕大多數(shù)在0.7以下,為減少單一屬性進(jìn)行預(yù)測(cè)時(shí)所帶來(lái)的誤差,需要選取兩個(gè)相關(guān)性系數(shù)較高的地震屬性進(jìn)行屬性融合。最終優(yōu)選出高階累量方差(R2=0.639)和平均振幅(R2=0.637)作為預(yù)測(cè)底板砂地比的屬性組合,高階累量方差(R2=0.750)和平均振幅(R2=0.642)作為預(yù)測(cè)頂板砂地比的屬性組合。

        2.5 砂地比預(yù)測(cè)

        經(jīng)優(yōu)選得出的屬性組合還需進(jìn)行回歸分析驗(yàn)證其是否能有效地運(yùn)用于研究區(qū)砂地比預(yù)測(cè),進(jìn)而得出屬性融合后的擬合公式[24]。首先對(duì)15號(hào)煤層底板優(yōu)選敏感屬性(平均振幅和高階累量方差)進(jìn)行回歸分析(表5):

        表5 15號(hào)煤層底板屬性組合回歸分析Table 5 Coal No.15 floor attribute combination regression analysis

        本次回歸分析得出的擬合公式:

        Y=0.084+0.525X1+0.476X2

        (2)

        式中:Y代表15號(hào)煤層底板砂地比;X1代表經(jīng)標(biāo)準(zhǔn)化后的高階累量方差值;X2代表經(jīng)標(biāo)準(zhǔn)化后的平均振幅值。

        (a)高階累量方差與砂地比擬合圖;(b)均方根振幅與砂地比擬合圖;(c)能量半衰時(shí)與砂地比擬合圖;(d)弧長(zhǎng)與砂地比擬合圖;(e)最大振幅與砂地比擬合圖;(f)平均振幅與砂地比擬合圖圖5 15號(hào)煤層底板砂地比與地震屬性交匯分析Figure 5 Coal No.15 floor total sandstone thickness to strata thickness ratio and seismic attribute cross analysis

        (a)高階累量方差與砂地比擬合圖;(b)能量半衰時(shí)與砂地比擬合圖;(c)均方根振幅與砂地比擬合圖;(d)弧長(zhǎng)與砂地比擬合圖;(e)最大振幅與砂地比擬合圖;(f)平均振幅與砂地比擬合圖圖6 15號(hào)煤層頂板砂地比與地震屬性交匯分析Figure 6 Coal No.15 rooftotal sandstone thickness to strata thickness ratio and seismic attribute cross analysis

        聚煤前地震屬性組合經(jīng)回歸分析之后得出線(xiàn)性回歸系數(shù)為0.800 339,大于0.8,說(shuō)明該屬性組合可用來(lái)預(yù)測(cè)煤層底板砂地比。因此可根據(jù)回歸分析所得出的擬合公式結(jié)合克里金插值方法對(duì)15號(hào)煤層底板砂地比進(jìn)行預(yù)測(cè),編制砂地比等值線(xiàn)圖[25](圖7)。

        圖7 基于融合屬性預(yù)測(cè)15號(hào)煤層底板砂地比分布Figure 7 Coal no.15 floortotal sandstone thickness to strata thickness ratio distributions based on attribute synthesis prediction

        其次對(duì)15號(hào)煤層頂板優(yōu)選敏感屬性(平均振幅和高階累量方差)進(jìn)行回歸分析(表6):

        表6 15號(hào)煤層頂板屬性組合回歸分析Table 6 Coal No.15 roof attribute combination regression analysis

        本次回歸分析得出的擬合公式:

        Y=1.325-2.461X1+0.046X2

        (3)

        式中:Y代表15號(hào)煤層頂板砂地比;X1代表經(jīng)標(biāo)準(zhǔn)化后的高階累量方差值;X2代表經(jīng)標(biāo)準(zhǔn)化后的平均振幅值。

        聚煤后地震屬性組合經(jīng)回歸分析之后得出其線(xiàn)性回歸系數(shù)為0.866 094,大于0.8,說(shuō)明該屬性組合可用來(lái)預(yù)測(cè)煤層頂板砂地比。因此可根據(jù)回歸分析所得出的擬合公式結(jié)合克里金插值方法對(duì)15號(hào)煤層頂板砂地比進(jìn)行預(yù)測(cè),編制砂地比等值線(xiàn)圖(圖8)。

        本文在屬性融合得出研究區(qū)聚煤前后砂地比等值線(xiàn)圖后,結(jié)合區(qū)域地質(zhì)資料、沉積環(huán)境、單井相以及地震資料對(duì)聚煤前后沉積相進(jìn)行刻畫(huà)。在晚古生代時(shí)期,沁水盆地主要為海陸交互相沉積,在太原組15號(hào)煤層聚煤前期,經(jīng)歷了障壁島—潮坪—潟湖—潮坪環(huán)境,在潮坪沉積相中識(shí)別出了混合坪沉積亞相,在15號(hào)煤層聚煤后期,沉積環(huán)境為潮坪相沉積,在潮坪相沉積中識(shí)別出了泥炭坪和混合坪兩種沉積亞相(圖9)。

        砂地比圖8 基于融合屬性預(yù)測(cè)15號(hào)煤層頂板砂地比分布Figure 8 Coal no.15 roof total sandstone thickness to strata thickness ratio distributions based on attribute synthesis prediction

        (a)聚煤前沉積相 (b)聚煤后沉積相圖9 15號(hào)煤層聚煤前后沉積相分布Figure 9 Sedimentary facies distribution before and after coal No.15 accumulation

        3 結(jié)論

        本文充分結(jié)合地球物理測(cè)井與三維地震勘探優(yōu)勢(shì),以陽(yáng)泉礦區(qū)寺家莊煤礦為例,基于地震沉積學(xué)的理論方法,利用測(cè)井資料與地震資料,結(jié)合研究區(qū)地質(zhì)資料,對(duì)研究區(qū)聚煤前后砂地比分布特征進(jìn)行了預(yù)測(cè),并對(duì)其沉積相進(jìn)行了刻畫(huà)。主要得出以下結(jié)論:

        1)高階累量方差和平均振幅可作為預(yù)測(cè)研究區(qū)15號(hào)煤層頂?shù)装迳暗乇葏?shù)的敏感屬性組合。

        2)結(jié)合單井相,地震資料以及地質(zhì)背景對(duì)研究區(qū)沉積相進(jìn)行了識(shí)別,聚煤前為障壁島—潮坪—潟湖—潮坪相沉積,聚煤后為潮坪相沉積。在潮坪相的基礎(chǔ)上又識(shí)別出了泥炭坪,混合坪兩種沉積亞相。

        猜你喜歡
        沉積相振幅切片
        四川盆地寒武系底部黑色巖段對(duì)比與沉積相初探
        十大漲跌幅、換手、振幅、資金流向
        十大漲跌幅、換手、振幅、資金流向
        基于SDN與NFV的網(wǎng)絡(luò)切片架構(gòu)
        十大漲跌幅、換手、振幅、資金流向
        滬市十大振幅
        腎穿刺組織冷凍切片技術(shù)的改進(jìn)方法
        冰凍切片、快速石蠟切片在中樞神經(jīng)系統(tǒng)腫瘤診斷中的應(yīng)用價(jià)值比較
        中東卡塔爾地區(qū)二疊系Unayzah組沉積相分析
        大安北油田泉頭組四段沉積相研究
        日本一道高清在线一区二区| 国产成人亚洲综合| www国产亚洲精品久久麻豆| 国产成人亚洲精品青草天美| 色偷偷久久一区二区三区| 国产精品午睡沙发系列| 精品综合久久久久久99| 成人影院羞羞的视频免费观看| 久久久久av综合网成人| 免费特级毛片| 亚洲高潮喷水中文字幕| 麻豆av一区二区天堂| 少妇高潮精品在线观看| 亚洲国产精品无码中文字| 国产午夜精品理论片| 亚洲另类自拍丝袜第五页| 动漫在线无码一区| 杨幂国产精品一区二区| 国产成人亚洲精品91专区高清 | 亚洲成a人无码| 五月天丁香久久| 久久精品国产只有精品96| 国产目拍亚洲精品二区| 久草中文在线这里只有精品| 少妇被粗大的猛进出69影院| 日本老熟妇50岁丰满| 精品视频在线观看免费无码| 亚洲小少妇一区二区三区| 最新欧美精品一区二区三区| 一本色道无码道在线观看| 日本亚洲欧美高清专区| 少妇人妻偷人中文字幕| 成人久久黑人中出内射青草| 十八禁视频网站在线观看| 久久精品国产亚洲av忘忧草18| 久久国产乱子伦精品免费强| 在线观看国产自拍视频| 人妖av手机在线观看| 一本色综合久久| 国产片三级视频播放| 亚洲日本中文字幕高清在线|