999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

低頻可控震源“兩寬一高”地震勘探的應用

2017-12-18 10:28:09張麗艷于常青
石油地球物理勘探 2017年6期

張麗艷 李 昂 于常青

(①中國地質調查局沈陽地質調查中心,遼寧沈陽 110034; ②中國地質科學院地質研究所,北京100037)

·綜合研究·

低頻可控震源“兩寬一高”地震勘探的應用

張麗艷*①李 昂①于常青②

(①中國地質調查局沈陽地質調查中心,遼寧沈陽 110034; ②中國地質科學院地質研究所,北京100037)

“兩寬一高”是指寬方位、寬頻帶、高密度三維地震勘探,能夠記錄更完整的波場信息,假頻少、低頻信息豐富,有利于薄互層油氣藏的刻畫和地層各向異性研究。以實際的“兩寬一高”資料為基礎,從原始資料的波場特征、信噪比及頻帶等三方面對比分析了“兩寬一高”資料的優勢和特點,確定了針對“兩寬一高”資料的處理技術和流程。利用最終成像結果對斷層和砂體進行刻畫和描述,體現了“兩寬一高”地震勘探的優勢。

“兩寬一高” 可控震源 低頻 寬頻帶 高分辨率

1 引言

如今油氣勘探開發面臨的地質問題越來越復雜,難度越來越大[1,2]。在油田勘探開發的中后期對三維地震技術的要求越來越高,如何綜合利用地質和地震資料進行油氣勘探,并且使地震的作用最大化,是地球物理學家追求的目標之一,地震勘探也由原來的常規三維發展到“兩寬一高”三維。“兩寬一高”地震指寬方位、寬頻帶和高密度的地震采集。寬方位是指地震勘探三維排列的橫縱比大于0.5,目前寬方位勘探的橫縱比一般都在0.7以上[3]; 寬方位地震勘探可以獲得較完整的地震波場信息,提高斷層及巖性的識別能力。寬頻帶是指采集的原始資料相對頻寬——倍頻程增加,即在同一分貝下,低截頻值變小,高截頻值變大,低頻端和高頻端同時向兩端擴展。低頻有效成份在地震勘探采集時越來越重要,直接關系到后續的寬頻處理; 另外,低頻信息在全波形反演及油氣檢測中的作用也非常關鍵[4]。高密度是指小面元地震采集,同時覆蓋次數也大幅度增加,有利于地質體的橫向分辨率和縱向分辨率的提高。“兩寬一高”三維地震勘探技術逐漸成熟,在國內外地震勘探中得到了廣泛應用,如新疆塔里木油田“兩寬一高”地震采集的密度已達到1070萬道/km2; 在準噶爾盆地實施了多塊“兩寬一高”三維地震采集,大幅度提高了地震資料的品質; 遼河油田的“兩寬一高”地震勘探面元大小為10m×10m,覆蓋次數為256次,橫縱比為0.91,接近于全方位采集。在海外的AGM區塊及乍得Bonger盆地都采集了“兩寬一高”地震資料,取得了較好的地震勘探效果[5-7]。“兩寬一高”三維地震勘探激發方式大多采用低頻可控震源激發[8]。與井炮激發方式相比,低頻可控震源有其獨特的優勢,每日采集可達到2~4萬炮,在提高采集效率、降低成本的同時,降低了HSE風險。低頻可控震源的掃描頻帶為1.5~96Hz,可達6個倍頻程,對于薄儲層的研究具有重要意義。大慶油田已全面進入精細勘探階段,儲層高含水,單砂體薄,2~3m的薄儲層占80%,儲層橫向變化快,水平井井位部署難度大,如何利用高分辨率三維精細地震勘探提高2~3m薄儲層的預測精度是目前三維地震勘探面臨的主要難題。2015年大慶油田在三肇地區利用低頻可控震源采集了“兩寬一高”三維地震資料,本文以此為基礎,分析了可控震源資料的固有特征和資料處理難點,針對難點開展研究,確立了針對可控震源的“兩寬一高”的處理關鍵技術,重點突出寬頻帶、寬方位的高分辨率保幅處理,并取得了較好的成像效果,為后續的高精度薄儲層預測和反演提供了高質量的數據。

2 采集參數及原始資料特征分析

2015年大慶油田在三肇地區采集了“兩寬一高”地震資料,表1為新、老(常規三維地震)采集參數的對比。“兩寬一高”資料由低頻可控震源激發,掃描頻帶為1.5~96Hz,單點模擬檢波器接收,橫縱比為0.83,覆蓋次數為480次,面元尺寸為5m×10m;老資料采用炸藥震源激發,組合檢波器接收,面元尺寸為20m×20m,橫縱比為0.47,屬常規窄方位采集。“兩寬一高”觀測系統的炮密度是常規老資料的6.4倍,道密度是53倍,接收線距是0.6倍,是真正的寬方位、高密度、高覆蓋、寬頻帶三維地震采集。

圖1是低頻可控震源“兩寬一高”和常規三維地震的原始資料的對比,也是單點檢波器接收和組合檢波器接收的對比,可以看出: 由于“兩寬一高”是高密度地震,排列長、道距小、面元小,地震信號的空間采樣率得到了提高,在滿足空間采樣定理的前提下,有效波和干擾波都能得到較完整的刻畫(圖中紅色方框),更有利于處理過程中的噪聲去除; 從信噪比方面看,由于常規三維地震勘探采用的是組合檢波器接收,具有一定的去噪作用,信噪比相對較高,面波視速度小、頻率低,油井干擾較少; 單點檢波器接收,面波頻帶寬、視速度大,油井和外源干擾也非常嚴重(圖1a中藍色方框)。除此之外,還存在可控震源的特有噪聲干擾——諧波[9,10]。當可控震源進行滑動掃描時,由于前一次掃描還沒有結束,后一次掃描已經開始,對于升頻掃描方式,下一組震源的諧波與上一組震源的基波信號混疊在一起,在相關后的地震記錄上形成了諧波干擾,而且滑動時間越短,諧波干擾越強,當勘探目標為弱反射層或者薄層時,諧波的影響不可忽略。諧波干擾橫向上波及十幾道的范圍,隨著時間的持續,能量逐漸增強,頻率主要在40~60Hz之間。圖1c和圖1d為諧波干擾在單炮上的表現特征、頻譜及時頻譜。

表1 “兩寬一高”和常規三維地震觀測系統的對比

圖1 “兩寬一高”資料和常規三維地震資料原始單炮特征分析

圖2是圖1a、圖1b單炮記錄1.0~2.0s時窗的頻譜,可以看出,“兩寬一高”資料的頻寬為4~45Hz,常規三維地震的頻寬為9~50Hz,可控震源資料在低頻端信息更豐富,高頻端的信號能量比炸藥震源稍弱。

選取“兩寬一高”資料和常規三維資料原始單炮記錄的相同時窗進行六個頻段的頻率掃描:2Hz以下、4~8Hz、8~16Hz、16~32Hz、32~64Hz、70Hz以上。圖3和圖4分別是二者的單炮分頻掃面記錄。從低頻端的頻率掃描來看,“兩寬一高”資料在2Hz頻率以下頻段目的層1.4s處(圖3a中紅色方框內)還能見到較強能量,而常規三維地震資料在2Hz以下幾乎見不到有效信號;隨著頻率的增加,在優勢頻帶內,“兩寬一高”資料的品質非常高。“兩寬一高”資料在70Hz以上頻段0.8s以后能量非常弱(圖3f中紅色箭頭),常規三維地震資料70H以上頻段1.4s以后的能量相對于“兩寬一高”資料較強,這是由于“兩寬一高”資料采用的是低頻可控震源激發,掃描截止頻率到96Hz,常規三維地震采用炸藥震源激發。如何提高可控震源的高頻信號的能量,是今后需要加強研究的方向。

圖2 “兩寬一高”(a)和常規(b)三維地震單炮記錄頻譜對比

圖3 “兩寬一高”地震資料的不同頻段掃描結果

圖4 常規三維地震資料的不同頻段掃描結果

3 “兩寬一高”處理技術

由于“兩寬一高”資料與常規三維資料相比具有高密度、寬方位、寬頻帶、高覆蓋的特點和優勢,所以在處理時應采取與常規三維地震不同的技術,在處理技術上應體現出“兩寬一高”資料的特點。圖5是二者處理技術流程的對比。圖5a流程中綠色部分是“兩寬一高”資料處理的關鍵技術。從處理流程可以看出,“兩寬一高”的處理流程重點體現了在激發方式不同的情況下,針對可控震源的處理技術,主要包括波場的特征分析、信噪分離、如何增強低頻弱信號的能量及擴展高頻能量的寬頻處理、OVT域的寬方位處理等。具體表現為以下四方面。

(1)噪聲去除。由于“兩寬一高”資料能夠記錄更完整的波場信息,所以噪聲的特征也表現得相對完整,尤其是面波的刻畫。針對“兩寬一高”資料面波的去除,應該采取更適合的方法。本文主要采用自適應面波衰減技術對面波進行有效去除。該技術提取面波的速度和頻率參數,通過可逆的子波變換,將輸入數據分解為不同頻率—波數域子集,在每個子集中通過波動方程進行彈性建模并計算自適應參數,最終再根據每個子集中面波特征在f-x域進行扇形濾波,減去不同階數的面波,該方法更適應“兩寬一高”資料對面波刻畫更完整的情況。第二類噪聲是針對諧波干擾的去除,方法很多[9-14],本文主要采取基于模型的相關方法對諧波干擾進行有效的去除,并取得了較好的效果。

(2)反褶積方面。“兩寬一高”的可控震源資料是零相位的,由于可控震源掃描信號激發(圖6a)所采集的地震數據是由可控震源子波與地下地層介質的反射系數序列的褶積結果,而可控震源子波是可控震源激發采集的地震數據經過自相關后得到的,為零相位(圖6b),所以可控震源采集的“兩寬一高”資料是零相位的。但是,為了滿足反褶積的假設條件,在反褶積前應該做最小相位化處理,這與炸藥震源資料處理不同。

(3)寬方位處理。“兩寬一高”資料的橫縱比達到0.87,是真正的寬方位資料,能夠全方位記錄地下的波場信息和照明引起的不均現象,所以針對“兩寬一高”資料的寬方位處理應該采取針對性的、能體現寬方位特征的處理,既OVT域的處理,在OVT域進行去噪、規則化、偏移等處理[15-19]。因為方位各向異性特征在寬方位資料上表現得更完整,有利于方位各向異性研究,進行方位各向異性校正,消除各向異性的影響,提高地震資料的分辨率。圖7a、圖7b為“兩寬一高”資料OVT域偏移后的道集校正前、后的對比,該OVT域的方位道集是在360°范圍內,地震波場隨方位的變化,體現了地下方位各向異性的特征,其中紫色的線代表炮檢距的變化,綠色的線代表方位角0°~360°內的變化。從圖中可以看出,由于方位各向異性的影響,偏移后的方位道集存在波浪形曲線形狀,這種特征在圖中1.0~1.2s表現尤其明顯(圖中紅色方框標識處),波浪形曲線的頂點和底點分別代表各向異性的主方向和垂直方向(也就是裂縫的走向和垂直方向),通過方位各向異性校正后,消除了這種方位各向異性的影響,同相軸校平。圖7c和圖7d分別是方位各向異性校正前、后成像剖面的對比,可以看出,校正后的剖面在成像精度和信噪比方面都要優于校正前的剖面,而且經過方位各向異性校正后,層間的弱信號能量得到加強(圖中黑色箭頭指示處)。

圖5 “兩寬一高”資料(a)和常規三維(b)處理技術流程對比

圖6 可控震源掃描信號(上)及自相關子波(下)

(4)寬頻處理。可控震源掃描頻率為1.5~96Hz,低頻端具有較強的優勢,但是能量相對較弱,為了增強低頻弱信號的能量,同時拓寬高頻有效信號成分,在提高倍頻程的前提下,提高其絕對頻寬,

圖7 OVT域偏移后的方位各向異性校正效果

使得“兩寬一高”資料達到真正意義的寬頻,所以在“兩寬一高”的處理流程上體現了針對可控震源的近地表吸收衰減補償技術和低頻補償技術[20]。可控震源的近地表吸收衰減補償技術,主要是在進行相位調整的同時,進行高頻能量的擴展。由于可控震源激發是在地表進行的,近地表的吸收衰減既包含從炮點傳播經過非彈性地層的衰減,又包含檢波點在非彈性地層中傳播時的吸收衰減,所以針對可控震源資料的近地表吸收衰減補償需要對炮點和檢波點都進行補償。低頻補償技術主要是針對可控震源的低頻信號進行的,低頻震源可以輸出豐富的低頻信號,但是由于儀器響應以及地層的吸收衰減等因素的影響會造成低頻的損失,這部分損失的低頻能量需要進行合理的補償處理,才能充分發揮低頻信息的作用。基于地震數據子波的低頻補償技術,針對頻譜特點合理地進行低頻補償,拓寬有效信號的低頻成分。

圖8為采用上述兩流程處理得到的最終成像剖面及頻譜對比。從頻譜圖(圖8c和圖8d)可以看出:常規三維資料的有效頻帶為10~82Hz,約為3個倍頻程; “兩寬一高”有效頻帶為3~80Hz,約為4個倍頻程,高頻成份略低、但低頻信息相對較豐富。倍頻程的增加可以減少視分辨率的假象,低頻成份對提高地震資料分辨率、識別薄儲層也起著重要作用。從成像剖面可以看出“兩寬一高”資料能夠使小斷層成像更清晰,斷點收斂更干脆(圖8a、圖8b中黑色箭頭標注處),圖8a中左側兩處斷層,在1.0~1.4s之間整個斷面成像從淺到深成像都非常清楚(綠色箭頭標注處),優于常規三維資料。

圖8 “兩寬一高”資料和常規三維處理結果對比

4 “兩寬一高”資料識別小斷層和單砂體

圖9是“兩寬一高”三維地震資料和常規三維地震資料的相干體切片對比,可以看出,兩種資料都能反映出該地區斷層的空間展布、走向及發育情況,但是“兩寬一高”資料能更清楚地反映斷層平面組合關系、連續性和線性特征(圖中橢圓標注處)。由于“兩寬一高”資料的低頻較豐富,成像剖面的低頻可達3Hz。低頻成分具有相對穩定、傳播距離遠、穿透能力強,能更好地識別深層或屏蔽層下的地質目標,增強特殊巖性體的識別能力。低頻成分的增加可以減少視分辨率的假象,在反演時能夠降低地震反演對測井資料的依賴,有效提高波阻抗反演的精度。從反演剖面(圖10)可以看出,“兩寬一高”資料縱向分辨率比常規三維資料有較為明顯的提升,圖中黑色方框內標注的Fangfu124-53井在0.8~1.0s之間兩套砂體的疊置關系得以突出,反映出第一套砂體為指狀的三角洲前緣相沉積,縱向上砂層較薄、橫向上有一定的延伸長度,第二套為具有正韻律結構的河道砂體沉積,縱向砂體較厚,橫向上成“短軸狀”。而常規資料由于分辨率較低,上面的砂體完全與上覆地層混在一起,下面的砂體也沒有刻畫出來。整體上“兩寬一高”資料與井上砂巖有較好的對應關系,吻合程度也高于常規三維資料。

圖9 “兩寬一高”資料和常規三維對小斷層刻畫對比

圖10 “兩寬一高”資料和常規三維資料對砂體刻畫的對比

5 結論和認識

(1)“兩寬一高”是精細三維地震的發展趨勢,將成為巖性油藏地震勘探的主要技術手段;

(2)低頻可控震源掃描頻率為1.5~96Hz,能夠實現“兩寬一高”三維地震經濟、高效采集,增加低頻成分、提高分辨率,實現真正意義的寬頻采集;

(3)“兩寬一高”三維地震處理技術與常規三維地震處理技術存在很大的差異,主要體現在如何提高低頻和高頻有效信號能量的寬頻處理和更加注重了OVT域內的方位處理及方位各向異性的研究;

(4)“兩寬一高”地震資料面元小、覆蓋次數高,縱向分辨率和橫向分辨率都優于常規三維地震資料,無論在刻畫斷層平面展布還是砂體描述方面都優于常規三維地震資料,且與井的吻合度較高。

感謝東方地球物理公司研究院提供了技術支持。

[1] 張保慶,周輝,左黃金等.寬方位地震資料處理技術及應用效果.石油地球物理勘探,2011,46(3):396-400.

Zhang Baoqing,Zhou Hui,Zuo Huangjin et al.Wide azimuth data processing techniques and their application.OGP,2011,46(3):396-400.

[2] 張軍華,朱煥,鄭旭剛等.寬方位角地震勘探技術評述.石油地球物理勘探,2007,42(5):603-609.

Zhang Junhua,Zhu Huan,Zheng Xugang et al.Summary of wide-azimuth seismic exploration technique.OGP,2007,42(5):603-609.

[3] 王學軍,于寶利,趙小輝等.油氣勘探中“兩寬一高”技術問題的探討與應用.中國石油勘探,2015,20(5):41-53.

Wang Xuejun,Yu Baoli,Zhao Xiaohui et al.Development and application of “2W1H” technique in oil and gas exploration.China Petroleum Exploration,2015,20(5):41-53.

[4] 孔德政,于敏杰,劉新文等.兩寬一高地震采集技術在復雜山前帶的應用及效果分析.新疆石油天然氣,2016,12(1):33-38.

Kong Dezheng,Yu Minjie,Liu Xinwen et al.Application and effects of the “two-wide one-high” seismic acquisition technique in complex piedmont.Xinjiang Oil & Gas,2016,12(1):33-38.

[5] 王兆磊,公亭,李隆梅等.高密度寬方位地震資料處理技術研究進展.物探化探計算技術,2015,37(4):465-471.

Wang Zhaolei,Gong Ting,Li Longmei et al.The new progress of high-density and width azimuth seismic data processing.Computing Techniques for Geophysical and Geochemical Exploration,2015,37(4):465-471.

[6] 凌云研究小組,寬方位角地震勘探研究.石油地球物理勘探,2003,38(4): 350-357.

Ling Yun Study Group.Application and study on wide azimuth seismic exploration.OGP,2003,38(4):350-357.

[7] 裴家學.寬方位地震資料在陸西凹陷勘探中的應用.大慶石油地質與開發,2015,34(5):146-150.

Pei Jiaxue.Application of the wide azimuth seismic data in the exploration of Luxi Sag.Petroleum Geology and Oilfield Development in Daqing,2015,34(5):146-150.

[8] 汪恩華,趙邦六,王喜雙等.中國石油可控震源高效地震采集技術應用與展望.中國石油勘探,2013,18(5):24-34.

Wang Enhua,Zhao Bangliu,Wang Xishuang et al.Application and outlook of vibroseis acquisition techniques with high efficiency of CNPC.China Petroleum Exploration,2013,18(5):24-34.

[9] 伍建,王潤秋,魏加明等.可控震源地震數據諧波濾除方法.石油地球物理勘探,2014,49(1):47-52.

Wu Jian,Wang Runqiu,Wei Jiaming et al.A method for vibroseis data harmonic filtering.OGP,2014,49(1):47-52.

[10] 曲英銘,李振春,韓文功等.可控震源高效采集數據特征干擾壓制技術.石油物探,2016,55(3):395-407.

Qu Yingming,Li Zhenchun,Han Wengong et al.The elimination technology for special interference in vibroseis efficient acquisition data.GPP,2016,55(3):395-407.

[11] 趙殿棟,胡立新,宋桂橋等.可控震源時變滑動掃描及低頻能量補償技術.中國石油勘探,2016,21(3):116-123.

Zhao Diandong,Hu Lixin,Song Guiqiao et al.Tech-nology of variable slip time sweep and low-frequency energy compensation with vibroseis operation.China Petroleum Exploration,2016,21(3):116-123.

[12] 林娟,羅勇,劉宜文等.可控震源滑動掃描諧波干擾壓制方法.石油地球物理勘探,2014,49(5):852-856.

Lin Juan,Luo Yong,Liu Yiwen et al.Harmonic interference suppression on vibroseis slip sweep data.OGP,2014,49(5):852-856.

[13] 胡江濤,王華忠,王兆磊.利用可控震源諧波提高地震數據分辨率.石油地球物理勘探,2013,48(6):876-881.

Hu Jiangtao,Wang Huazhong,Wang Zhaolei.Broadband seismic data with vibroseis harmonics.OGP,2013,48(6):876-881.

[14] 劉玉海,尹成,潘樹林等.基于非線性掃描技術的可控震源地震勘探分辨率提高方法.東北石油大學學報,2013,37(2):55-62.

Liu Yuhai,Yin Cheng,Pan Shulin et al.A high resolution research with nonlinear signal to improve vibrator exploration based on nonlinear scanning.Journal of Northeast Petroleum University,2013,37(2):55-62.

[15] 段文勝,李飛,黃錄忠等.OVT域寬方位層析速度建模與深度域成像.石油地球物理勘探,2016,51(3):521-528.

Duan Wensheng,Li Fei,Huang Luzhong et al.Wide azimuth tomographic velocity model building and depth imaging in the OVT domain.OGP,2016,51(3):521-528.

[16] 詹仕凡,陳茂山,李磊等.OVT域寬方位疊前地震屬性分析方法.石油地球物理勘探,2015,50(5):956-966.

Zhan Shifan,Chen Maoshan,Li Lei et al.OVT-domain azimuth prestack seismic attribute analysis.OGP,2015,50(5):956-966.

[17] 段文勝,李飛,王彥春等.面向寬方位地震處理的炮檢距向量片技術.石油地球物理勘探,2013,48(2):206-213.

Duan Wensheng,Li Fei,Wang Yanchun et al.Offset vector tile for wide-azimuth seismic processing.OGP,2013,48(2):206-213.

[18] Schapper S,Jefferson R,Calvert A et al.Anisotropic velocities and offset vector tile prestack-migration processing of the Durham Ranch 3D,Northwest Colorado.The Leading Edge,2009,28(11):1352-1361.

[19] Wang Shirui,Ling Yun,Guo Xiangyu.Application of wide azimuth land 3D - a case history.SEG Technical Program Expanded Abstracts,2003,22:490-493

[20] 公亭,王兆磊,顧小弟等.寬頻地震資料處理配套技術.石油地球物理勘探,2016,51(3):457-466.

Gong Ting,Wang Zhaolei,Gu Xiaodi et al.Broadband seismic data matching processing.OGP,2016,51(3):457-466.

*遼寧省沈陽市黃河北大街280號中國地質調查局沈陽地質調查中心,110034。Email:moonliyan@sohu.com

本文于2017年1月4日收到,最終修改稿于同年9月8日收到。

本項研究受國家自然科學基金面上項目(41374128)資助。

1000-7210(2017)06-1236-10

張麗艷,李昂,于常青.低頻可控震源“兩寬一高”地震勘探的應用.石油地球物理勘探,2017,52(6):1236-1245.

P631

A

10.13810/j.cnki.issn.1000-7210.2017.06.014

(本文編輯:宜明理)

張麗艷 博士,工程師,1980年生;2003年獲中國石油大學(華東)勘查技術與工程專業學士學位;2008年獲中國石油大學(北京)地球探測與信息技術專業工學博士學位。現在中國地質調查局沈陽地質調查中心主要從事“兩寬一高”高分辨率地震資料處理及各向異性成像等方面的研究。

主站蜘蛛池模板: 在线亚洲小视频| 激情五月婷婷综合网| 97se亚洲综合| 欧美日韩国产在线播放| 伊人激情综合网| 国产成人一区二区| www.av男人.com| 国产精品尤物铁牛tv| 人妻中文久热无码丝袜| 国产欧美亚洲精品第3页在线| 精品亚洲国产成人AV| 免费毛片全部不收费的| 啪啪免费视频一区二区| 国产美女无遮挡免费视频网站 | 亚洲视频免费在线| 国产一区二区色淫影院| 免费国产高清精品一区在线| 久久99国产视频| 白浆视频在线观看| 久久超级碰| 欧美人在线一区二区三区| 欧美 亚洲 日韩 国产| 国产毛片基地| 成人福利视频网| 久久国产高清视频| 在线播放91| 欧美一级99在线观看国产| 毛片久久网站小视频| 婷婷综合在线观看丁香| 亚洲AV无码一区二区三区牲色| jizz在线免费播放| 亚洲视频在线青青| 国产在线精品人成导航| 青青青视频91在线 | 一级毛片a女人刺激视频免费| 欧美亚洲综合免费精品高清在线观看 | 欧洲亚洲一区| 久久久黄色片| 国产精品专区第一页在线观看| 天天躁狠狠躁| 欧美日韩导航| 男女男免费视频网站国产| 日韩国产黄色网站| 欧美激情视频一区| 精品欧美日韩国产日漫一区不卡| 夜夜高潮夜夜爽国产伦精品| 国产精品黑色丝袜的老师| аⅴ资源中文在线天堂| 青青草原国产免费av观看| 欧美成人免费一区在线播放| 成年av福利永久免费观看| 老司机久久99久久精品播放 | 日韩免费毛片视频| 国产精品亚洲精品爽爽 | 国产网友愉拍精品视频| 免费福利视频网站| WWW丫丫国产成人精品| 国产免费自拍视频| 夜夜操国产| 少妇精品网站| 亚洲日韩精品欧美中文字幕| 天天干伊人| 不卡的在线视频免费观看| 国产主播在线一区| 青青草综合网| 99在线视频精品| 欧美国产中文| 国产主播福利在线观看| 久久无码高潮喷水| 亚洲精品不卡午夜精品| 日韩在线欧美在线| 成人免费视频一区| 亚洲日韩在线满18点击进入| 在线日韩日本国产亚洲| 亚洲中文在线看视频一区| 麻豆精品视频在线原创| 天堂av综合网| 99热这里只有免费国产精品| 色悠久久综合| 亚洲精品你懂的| 日韩一区二区三免费高清| 99re热精品视频国产免费|