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

    超聲波檢測套管井水泥環(huán)缺陷仿真模擬

    2017-09-18 08:54:51解廷菊龍士國
    聲學(xué)技術(shù) 2017年4期
    關(guān)鍵詞:首波聲場方位

    解廷菊,龍士國

    ?

    超聲波檢測套管井水泥環(huán)缺陷仿真模擬

    解廷菊,龍士國

    (湘潭大學(xué)材料科學(xué)與工程學(xué)院,湖南湘潭411105)

    針對固井過程中由于頂替效率過低、氣竄等原因造成水泥環(huán)局部缺失的現(xiàn)象,基于超聲脈沖反射原理,運用有限元方法建立不同大小、不同方位的水泥環(huán)缺失模型,數(shù)值模擬了超聲波在套管井多層介質(zhì)中的傳播規(guī)律。通過提取回波首波幅度、首波的接收時間,對回波進行頻譜分析,對比了膠結(jié)良好和水泥環(huán)缺失兩種情況下的時域特征和頻域特征。在研究不同方位缺失影響時,引入了缺陷敏感頻率,發(fā)現(xiàn)在一定激發(fā)頻率下,有缺陷的位置上敏感頻率會出現(xiàn)極大值,不同方位的缺陷有不同的敏感頻率,而且區(qū)別較為顯著,以此可以判斷缺失方位。

    超聲檢測;水泥環(huán)缺失;缺失方位;頻譜

    0 引言

    在國內(nèi)外的油田井建井過程中,固井是油田井建井的關(guān)鍵環(huán)節(jié),也是保護油氣層和保證油氣井生產(chǎn)壽命的關(guān)鍵。固井質(zhì)量的好壞關(guān)系到油田井勘探技術(shù)的開發(fā)效益和產(chǎn)能建設(shè),具有十分重要的意義[1-2]。然而固井注水泥過程中水泥竄槽、開發(fā)過程中所采取的增產(chǎn)措施以及稠油熱采井中出砂所造成的應(yīng)力集中、注汽過程中的高溫?zé)釕?yīng)力等都會造成水泥環(huán)局部缺失[3]的現(xiàn)象。水泥環(huán)缺失情況檢測變得越來越重要,引起了人們的廣泛重視。國內(nèi)許多研究人員在水泥膠結(jié)測井方面進行了大量理論分析和實驗[4],但是在檢測水泥環(huán)缺失方面研究較少。車小花[5]等人采用物理模擬的方法研究了不同膠結(jié)狀況下的聲全波波形,通過分析波形判斷第一界面和第二界面的膠結(jié)質(zhì)量。周吟秋[6]等人討論了套管波首播幅度關(guān)于激發(fā)聲源中心頻率變化的關(guān)系。許飛龍[7]等人利用超聲反射法通過濾波檢測水泥膠結(jié)質(zhì)量,并建立有限聲束模型計算了斜入射時的反射全波,探討了入射角對反射波的影響。謝丹艷[8]等人對水平井中可能存在的水泥環(huán)上層弓形缺失,利用有限差分法實現(xiàn)了水平井弓形缺失復(fù)雜井孔聲場的數(shù)值模擬。宋若龍[9]等人利用應(yīng)力速度有限差分法并行程序?qū)μ坠芫葏^(qū)水泥缺失的非對稱井孔聲場進行研究,為反演方法提供數(shù)據(jù),完善了分區(qū)水泥膠結(jié)測井解釋方法。目前,中海油田服務(wù)有限公司研制出多功能超聲成像測井儀器,采用超聲脈沖反射測量技術(shù),進行固井質(zhì)量評價。

    有限元法是解決諸多復(fù)雜物理問題的有效手段,ANSYS是目前世界頂端的有限元商業(yè)應(yīng)用程序,是融結(jié)構(gòu)、流體、電場、磁場、聲場分析于一體的大型通用有限元分析軟件[10]。劉長福[11]等人使用ANSYS對超聲波傳播進行了仿真計算,結(jié)果能有效表現(xiàn)超聲場的特性?;诔暶}沖反射原理,本文采用ANSYS有限元計算軟件對超聲波在水泥環(huán)缺陷套管井多層介質(zhì)中的聲場進行仿真模擬,研究了不同缺陷大小以及不同方位缺陷對超聲反射回波的影響。

    1 井孔聲場模擬方法

    1.1 有限元模型的建立

    根據(jù)實際油田井結(jié)構(gòu)的形貌,本文建立了兩類套管井膠結(jié)良好和水泥環(huán)缺失的二維有限元物理模型。如圖1所示,模型由內(nèi)而外介質(zhì)層依次為流體、套管、水泥層、地層,套管內(nèi)半徑為40 mm,各層介質(zhì)參數(shù)見表1。換能器收發(fā)位置如圖1(a)所示,圖1(b)是水泥環(huán)缺陷模型,圖1(c)是套管井30°方位水泥環(huán)缺失15°模型。其中,對水泥環(huán)缺陷模型設(shè)置了11組不同的尺寸,分別沿軸周向缺失30°~360°。選擇井孔中心作為聲源激發(fā)點,管壁四周作為接收點。

    圖1 有或無水泥環(huán)缺失的套管井結(jié)構(gòu)物理模型

    表1 模型材料參數(shù)

    在確定有限元單元類型、模型材料參數(shù)后,則根據(jù)如圖1中所示的簡化物理模型生成相應(yīng)的有限元網(wǎng)格劃分模型。有限元模型的網(wǎng)格劃分在有限元建模過程中極其重要,網(wǎng)格劃分的合理性與計算結(jié)果的可靠性和準確性在很大程度上呈直接相關(guān)。如果將網(wǎng)格劃分太粗,計算精度將會有一定的偏差,如果將網(wǎng)格劃分太細,雖然提高了有限元模擬計算的精度,但會造成計算時間過長或占用電腦的內(nèi)存過大,甚至出現(xiàn)電腦存儲不夠的問題。

    為了保證計算精度,聲波沿傳播方向每個波長內(nèi)至少有10~20個單元[12]。既要滿足最大頻率下波長的要求,又要綜合考慮計算量。最大模型單元尺寸的選擇見表2。此外,對于邊界的處理,取地層厚度為管壁厚的10倍,聲波到達地層后接收到的回波在150 μs左右,對前面波形的分析干擾不大,加之ANSYS某些單元的局限性不能施加吸收邊界,故添加了約束。

    表2 模型最大單元尺寸

    1.2 聲源信號的選取

    本文采用的激勵信號由5個單音頻信號疊加并經(jīng)HANNING函數(shù)調(diào)制而成,為了更好地模擬出套管井結(jié)構(gòu)中的超聲波傳播,選擇的脈沖信號頻率=100 kHz,脈沖寬度=50 μs,計算時間=500 μs。這是因為超聲波頻率越高,其波長就越短,在傳播過程中衰減就越快,尤其在缺陷模型中需要繞過缺陷傳播,易產(chǎn)生頻散現(xiàn)象、穿透能力差、超聲波能量利用率低。而超聲波頻率越低,聲波能量就低,在應(yīng)力云圖中觀察不到明顯的傳播路徑。經(jīng)過多次模擬比對,選擇頻率為100 kHz進行模擬時接收到的波形較好,易于提取分析。聲源信號的波形及頻譜曲線如圖2、3所示。相對于一般如半正弦信號或矩形波、尖脈沖等信號,從應(yīng)力云圖明顯看到在傳播過程中能量衰減較快,接收到的波形雜亂,不利于后處理數(shù)據(jù)提取,而HANNING窗脈沖信號主瓣高,旁瓣瞬間衰減很快,其中心頻率附近集中著絕大部分頻譜能量,在信號識別中頻率敏感度高,這種窄帶激勵既可以使應(yīng)力波的傳播距離增大,又可以增強自身信號強度。

    在本模擬過程中先編寫信號函數(shù)表達式,然后以瞬時載荷的形式加載到井孔中央位置所對應(yīng)的節(jié)點上(即有限元創(chuàng)建模型井孔中央為坐標原點,劃分網(wǎng)格后在List/Keypoint/Coords+Attributes可查看到它所對應(yīng)的節(jié)點為34)。所加載激勵信號[13]的數(shù)學(xué)表達式為:

    其中:為脈沖信號頻率;為脈沖信號中單音頻數(shù);為時間。

    圖2 聲源信號時域波形

    Fig.2 Time domain waveform of the sound source signal

    2 模擬結(jié)果分析

    2.1 聲波在介質(zhì)中的傳播聲場

    圖4是超聲波在膠結(jié)良好和水泥環(huán)缺失套管井中不同時刻傳播的聲場快照。由圖4可見,在套管井中央加載一個聲源信號,超聲波就會在井孔流體中傳播,由于在鋼管中傳播速度快,傳至管壁后迅速在套管四周傳播,并在水泥界面處發(fā)生反射和透射。因為存在阻抗差異,大部分能量反射回去,只有一小部分能量進入水泥環(huán)向地層傳播,在水泥環(huán)與地層界面處又會反射和透射。當聲波遇到缺陷時,聲波會以擾動的形式繞到缺陷外側(cè)進而在外層介質(zhì)中傳遞能量,并在缺陷界面發(fā)生反射,如圖4 (b)所示;當缺陷增大時,透射能量明顯減少,缺陷兩側(cè)的反射能量增強。最后管壁四周接收到的為各個界面的多次反射波。

    2.2 波形結(jié)果分析

    圖5是不同缺陷大小模型在0°方位上接收到的回波對比圖,回波在波形上呈現(xiàn)幾個波包。通過時域分析,在30 μs時首先接收到波形,與預(yù)期計算時間29.6 μs基本一致。48.5 μs前的第一個波包主要為套管波,即井內(nèi)流體-套管界面與套管-水泥界面之間的套管內(nèi)來回反射波。在水泥環(huán)缺失角度較小時,大部分能量可以通過水泥環(huán)向地層傳播,套管波的波包幅度較小。當缺失角度達到180°后,水泥環(huán)的缺失造成套管內(nèi)外聲阻抗差異較大,能穿透套管進入地層的波能量較低,大部分能量在套管內(nèi)相互反射,套管波波包的幅度隨著缺失角度的繼續(xù)擴大呈增加趨勢,直至自由套管達到最大,即波能量只能在套管內(nèi)傳播,不會傳遞到外層介質(zhì)中去。這一結(jié)果符合聲波套管井水泥環(huán)缺失的基本理論,驗證了該模擬的準確性。

    圖6為水泥膠結(jié)良好和不同水泥環(huán)缺失大小套管井在0°方位接收到的首波幅度對比圖,提取其所有情況的幅度信息如圖7所示(幅度進行了同比放大)。由圖7可以看到,在水泥環(huán)缺失范圍小于60°內(nèi),無缺陷模型首波幅度最小。但隨著缺失角度的增大,首波幅度隨之緩慢增大,而且從圖6可觀察到首波信號較完整模型出現(xiàn)得較晚,可能局部缺失使得信號強度較小,到180°后幅度急劇減小,并保持上下微弱波動。這是因為在缺失角度較小時,在缺陷界面處也會發(fā)生反射,這些反射疊加會使得首波幅度有一定增加。而在缺失180°之后,缺失范圍較大,向外傳播的能量較低,0°方位接收到的多是套管內(nèi)的來回反射波,幅度自然減小。

    2.3 頻譜分析

    快速傅里葉變換(Fast Fourier Transform,F(xiàn)FT)是一種快速有效的時頻轉(zhuǎn)換方法。通過計算,二次回波出現(xiàn)的時間為197.7 μs,設(shè)置采樣頻率為1 MHz,截取回波250個點作為采樣點,對所有模型在0°、30°、60°、90°方位接收到的波形進行FFT,圖8至圖11分別是不同缺陷在0°方位的回波頻譜曲線。

    圖12將所有方位FFT頻譜圖中的高主頻值對應(yīng)的振幅幅值提取出來,與缺失大小形成對應(yīng)關(guān)系(幅度進行了同比放大)。頻譜圖中所示的幅值是由于介質(zhì)界面之間的多次反射產(chǎn)生瞬態(tài)共振所致,對應(yīng)的振幅幅值越高,表明反射次數(shù)越多,累計的反射能量越高。分析結(jié)果顯示,圖8無缺陷模型主頻值為996 kHz,而圖9缺陷模型出現(xiàn)單次主頻左移現(xiàn)象,降至956 kHz,這是由缺陷使得超聲波在介質(zhì)中傳播的高頻部分衰減大引起的,而在不同方位上接收到的幅值都呈一致性變化,隨著缺陷角度的增大,幅值先緩慢增大后明顯增大,當缺陷角度達到180°之后乃至自由套管,幅值急劇減小,且變化趨勢不明顯,趨于平衡。這說明缺陷的存在使得界面之間的反射加強,幅值增大,但當缺陷增大到一定程度時,能量不會傳遞到外層介質(zhì),幅值會明顯減小。

    3 水泥環(huán)缺失方位的確定

    為了分析水泥環(huán)局部缺失對聲場的影響,根據(jù)聲場的對稱性,分別模擬了水泥環(huán)在0°、30°、60°、90°不同方位含有15°的周向缺陷,其有限元物理模型如圖1(c)所示。各層介質(zhì)參數(shù)和模型單元尺寸選擇同表1、表2。通過對井壁四周不同方位接收到的回波進行分析,二次回波=250 μs,選取=1 MHz,分別截取無缺陷和有缺陷時的回波進行FFT,由頻譜幅度譜計算缺陷對超聲信號影響的敏感頻率[14],引入敏感頻率主要是為了更直觀得出在某缺失方位下回波FFT的不同特征,定義為

    從不同方位缺陷敏感頻率圖13可知:0°方位缺陷在203.2、422.3、446.2 kHz有明顯變化,且在446.2 kHz附近出現(xiàn)了極大值;30°方位缺陷在139.4 kHz附近最為敏感,同時在163.3、330.7 kHz也比較敏感;60°方位缺陷在338.6 kHz出現(xiàn)極大值,此頻率為最敏感,此外在322.7 kHz也較為敏感;90°方位缺陷變化幅度相對平緩,幅值較低,在366.5 kHz比較敏感。當改變激勵脈沖頻率時,敏感頻率也會發(fā)生變化,這樣可以結(jié)合不同激勵頻率檢測缺陷,對于更多外界因素如水泥層厚度、井徑等,都可以通過控制單一變量來研究敏感頻率的變化。

    (a)0°方位

    (b)30°方位

    (c)60°方位

    (d)90°方位

    圖13 不同方位上含有15°缺陷的敏感頻率

    Fig.13 Sensitive frequencies in different directions for defect of 15°

    由以上分析可知,不同方位的缺陷有不同的敏感頻率,而且區(qū)別較為顯著。通過分析不同方位的敏感頻率變化,結(jié)合頻域信息可以為進一步判斷水泥環(huán)缺失方位提供依據(jù)。

    4 結(jié)論

    (1) 當水泥環(huán)缺失范圍小于60°內(nèi),無缺陷模型首波幅度最小,但隨著缺失角度的增大,首波信號出現(xiàn)較晚,首波幅度隨之增大,到180°后幅度急劇減小,并保持平穩(wěn)。

    (2) 對所有回波進行頻譜分析,發(fā)現(xiàn)在不同方位上幅值都呈一致性變化,即隨著缺陷角度的增大,幅值增大,當達到一定程度(缺陷角度≥180°)后,幅值急劇減小,變化趨勢不明顯。

    (3) 不同方位的缺陷有不同的敏感頻率,這對判斷水泥環(huán)缺失方位有一定幫助。

    [1] ZHANG H, XIE D, SHANG Z, et al. Simulated various characteristic waves in acoustic full waveform relating to cement bond on the secondary interface[J]. Journal of Applied Geophysics, 2011, 73(2): 139-154.

    [2] 龐秉謙, 楊松, 竇益華, 等. 套管磨損與水泥環(huán)缺陷位置對套管應(yīng)力的影響[J]. 石油機械, 2009, 37(10): 1-3. PANG Bingqian, YANG Song, DOU Yihua, et al. The influence of casing wear and cement ring defect position on the casing stress[J]. Journal of petroleum machinery, 2009, 37(10): 1-3.

    [3] 趙鵬, 趙寶樂, 劉樂華. 水泥環(huán)缺陷對套管損壞的影響分析[J]. 油氣田地面工程, 2009, 28(10): 40-41. ZHAO Peng, ZHAO Baole, LIU Lehua. Cement ring defect analysis of the influence of the casing damage[J]. Journal of Oil and Gas Field Surface Engineering, 2009, 28(10): 40-41.

    [4] TANG Jun, ZHANG Chengguang, ZHANG Bixing, et al. Cement bond quality evaluation based on acoustic variable density logging[J]. Petroleum Exploration and Development Online, 2016, 43(03): 514-521.

    [5] 車小花, 喬文孝, 閆向宏. 不同膠結(jié)狀況對固井質(zhì)量測井影響的模擬實驗研究[J]. 測井技術(shù), 2005, 29(03): 185-187. CHE Xiaohua, QIAO Wenxiao, YAN Xianghong. Simulated experimental investigation about effects of different bongdings on cement bond quality[J]. Journal of Well Logging Technology, 2005, 29(03): 185-187.

    [6] 周吟秋, 王秀明, 陳德華. 套管井聲波測井中首波幅度衰減研究[J]. 測井技術(shù), 2007, 31(04): 321-326. ZHOU Yinqiu, WANG Xiuming, CHEN Dehua. First wave amplitude attenuation in cased well logging research[J]. Journal of Well Logging Technology, 2007, 31(04): 321-326.

    [7] 許飛龍, 胡恒山, 尹誠剛. 通過濾波檢測水泥膠結(jié)質(zhì)量的超聲反射法[J]. 聲學(xué)學(xué)報, 2015, 40(03): 429-436. XU Feilong, HU Hengshan, YIN Chenggang. Detection of cement bond quality by filtering the ultrasonic reflection method[J]. Acta Actustica, 2015, 40(03): 429-436.

    [8] 謝丹艷, 張宏兵, 孫樹林. 弓形缺失水平井井孔聲場數(shù)值模擬[J]. 測井技術(shù), 2015, 39(02): 150-154. XIE Danyan, ZHANG Hongbing, SUN Shulin. Numerical simulation of borehole acoustic field in horizontal well with missing arch[J]. Journal of Well Logging Technology, 2015, 39(02): 150- 154.

    [9] 宋若龍, 劉金霞, 姚桂錦, 等. 非軸對稱套管井中聲場的并行有限差分模擬[J]. 地球物理學(xué)報, 2010, 53(11): 2767-2775. SONG Ruolong, LIU Jinxia, YAO Guijin, et al. The parallel finite difference simulation of sound field in axisymmetric casing well[J]. Journal of Geophysics, 2010, 53(11): 2767-2775.

    [10] 高興軍, 趙恒華. 大型通用有限元分析軟件ANSYS簡介[J]. 遼寧石油化工大學(xué)學(xué)報, 2004, 24(03): 94-98. GAO Xingjun, ZHAO Henghua. Abstract of large finite element analysis software ANSYS[J]. Journal of Liaoning University of Petroleum Chemical Industry, 2004, 24(03): 94-98.

    [11] 劉長福, 牛曉光, 李中偉. 基于ANSYS的超聲縱/橫波傳播仿真計算[J]. 無損檢測, 2011, 30(06): 15-18. LIU Changfu, NIU Xiaoguang, LI Zhongwei. Ultrasonic longitudinal/transverse wave propagation based on ANSYS simulation[J]. Journal of Nondestructive Testing, 2011, 30(06): 15-18.

    [12] Song K I, Cho G C. Numerical study on the evaluation of tunnel shotcrete using the impact-echo method coupled with Fourier transform and short-time Fourier transform[J]. International Journal of Rock Mechanics and Mining Sciences, 2010, 47(8): 1274-88.

    [13] Aggelis D G, Tsimpris N, Chai H K, et al. Numerical simulation of elastic waves for visualization of defects[J]. Construction and Building Materials, 2011, 25(4): 1503-1507.

    [14] 孫繼華, 李書光, 倪云鹿, 等. 固體中缺陷對聲波頻譜特性影響的有限元仿真[J]. 計算機仿真, 2008, 25(08): 308-311. SUN Jihua, LI Shuguang, NI Yunlu, et al. Finite element simulation of the effect of defects in solid on ultrasonic spectral feature[J]. Computer Simulation, 2008, 25(08): 308-311.

    Simulation of cement defect in casing well with ultrasonic

    XIE Ting-ju, LONG Shi-guo

    (College of Materials Science and Engineering, Xiangtan University, Xiangtan 411105,Hu’nan,China)

    For studying the patial defect phenomenon of cement sheath caused by low displacement efficiency and gas channeling in the cementing process, the finite element method is used to establish different size and different orientation models of cement sheath defect and to simulate the propagation rules of ultrasound in multilayer medium of casing well. By extracting the echo amplitude and the headwave time in different directions and by analyzing ccho frequency spectrum, the time domain and frequency domain characteristics of the models are compared effectively. For determining the cement sheath defect at different orientations, the defect’s sensitive frequency (Dev) is introduced, and it is found that a great value of Dev in the defective position appears under a certain excitation frequency and the defects in different directions show different Devs, which can determine the orientation of the defect.

    ultrasonic testing; cement sheath defect; defect orientation; frequency spectrum

    TB553

    A

    1000-3630(2017)-04-0346-07

    10.16300/j.cnki.1000-3630.2017.04.009

    2016-08-10;

    2016-10-24

    湖南省教育廳重點科學(xué)研究資助項目(16A209)

    解廷菊(1990-), 女, 山西朔州人, 碩士研究生, 研究方向為材料測控技術(shù)與裝備。

    解廷菊, E-mail: 1318346903@qq.com

    猜你喜歡
    首波聲場方位
    認方位
    幼兒園(2021年12期)2021-11-06 05:10:20
    超聲波測量鉆井卡點的方法研究
    基于BIM的鐵路車站聲場仿真分析研究
    探尋360°全聲場發(fā)聲門道
    基于首波寬度調(diào)整自動增益控制方法的研究
    借助方位法的拆字
    中國修辭(2016年0期)2016-03-20 05:54:32
    說方位
    幼兒100(2016年28期)2016-02-28 21:26:17
    基于TMS320C6678的SAR方位向預(yù)濾波器的并行實現(xiàn)
    板結(jié)構(gòu)-聲場耦合分析的FE-LSPIM/FE法
    一種脈沖貝塞爾波的構(gòu)造及其非線性聲場的仿真
    亚洲精品国产精品久久久不卡| 丁香欧美五月| 不卡一级毛片| 三上悠亚av全集在线观看| av线在线观看网站| 国产一区有黄有色的免费视频| 黄色女人牲交| 1024香蕉在线观看| 99re6热这里在线精品视频| 精品一品国产午夜福利视频| 高清视频免费观看一区二区| 天天添夜夜摸| 久久亚洲真实| 久久久精品国产亚洲av高清涩受| 熟女少妇亚洲综合色aaa.| 黑人欧美特级aaaaaa片| 精品国产国语对白av| 久久精品亚洲熟妇少妇任你| 久久精品亚洲精品国产色婷小说| 色老头精品视频在线观看| 超碰97精品在线观看| av天堂在线播放| 欧美黑人精品巨大| 国产一区二区三区综合在线观看| 久久中文字幕人妻熟女| 99国产精品一区二区蜜桃av | 中文字幕高清在线视频| 搡老熟女国产l中国老女人| 午夜91福利影院| 成人精品一区二区免费| 久久国产亚洲av麻豆专区| 国产成人精品久久二区二区免费| 国产成+人综合+亚洲专区| 色婷婷久久久亚洲欧美| 美国免费a级毛片| 一边摸一边抽搐一进一出视频| 麻豆成人av在线观看| 不卡一级毛片| 涩涩av久久男人的天堂| 国产精品国产av在线观看| 美女高潮喷水抽搐中文字幕| 日韩熟女老妇一区二区性免费视频| 性色av乱码一区二区三区2| 老司机午夜十八禁免费视频| 亚洲va日本ⅴa欧美va伊人久久| 免费看a级黄色片| 日韩熟女老妇一区二区性免费视频| 在线观看免费午夜福利视频| 精品国产国语对白av| 国产精品影院久久| 一夜夜www| 18禁裸乳无遮挡动漫免费视频| 搡老熟女国产l中国老女人| 亚洲少妇的诱惑av| 黄频高清免费视频| √禁漫天堂资源中文www| 欧美激情久久久久久爽电影 | 午夜视频精品福利| 日韩视频一区二区在线观看| 三级毛片av免费| 国产成人精品久久二区二区免费| 色在线成人网| 老鸭窝网址在线观看| www.精华液| 欧美不卡视频在线免费观看 | 在线看a的网站| 国产单亲对白刺激| 国产高清国产精品国产三级| 窝窝影院91人妻| 欧美日韩国产mv在线观看视频| 一进一出抽搐动态| 成熟少妇高潮喷水视频| 亚洲在线自拍视频| 久久ye,这里只有精品| svipshipincom国产片| www.精华液| 人妻丰满熟妇av一区二区三区 | 男人舔女人的私密视频| 婷婷精品国产亚洲av在线 | 亚洲精品国产一区二区精华液| 欧美日韩亚洲国产一区二区在线观看 | 日本vs欧美在线观看视频| netflix在线观看网站| av天堂在线播放| 91av网站免费观看| 国产高清国产精品国产三级| 国产精品久久久久久精品古装| 99re6热这里在线精品视频| 欧美人与性动交α欧美精品济南到| 亚洲精品久久午夜乱码| 每晚都被弄得嗷嗷叫到高潮| 亚洲性夜色夜夜综合| 亚洲免费av在线视频| 老司机深夜福利视频在线观看| 免费少妇av软件| 国产成人精品久久二区二区91| 大片电影免费在线观看免费| 成人特级黄色片久久久久久久| 人人妻人人爽人人添夜夜欢视频| 日本一区二区免费在线视频| 不卡av一区二区三区| 日韩 欧美 亚洲 中文字幕| 欧美在线一区亚洲| 精品久久久久久久久久免费视频 | 免费看a级黄色片| 天堂中文最新版在线下载| 777米奇影视久久| 久久精品亚洲精品国产色婷小说| 久久精品国产清高在天天线| 后天国语完整版免费观看| 啦啦啦在线免费观看视频4| 伊人久久大香线蕉亚洲五| 亚洲av电影在线进入| 国产97色在线日韩免费| 麻豆成人av在线观看| 一本一本久久a久久精品综合妖精| 欧美日韩视频精品一区| 一本综合久久免费| 午夜免费观看网址| 精品国产一区二区久久| 国产亚洲精品一区二区www | 91精品国产国语对白视频| 精品亚洲成a人片在线观看| 午夜亚洲福利在线播放| 91大片在线观看| 欧美亚洲日本最大视频资源| 中文字幕人妻丝袜制服| 欧美大码av| 欧美亚洲日本最大视频资源| 精品久久久久久久毛片微露脸| 国产精品亚洲av一区麻豆| 在线国产一区二区在线| a级毛片黄视频| 波多野结衣一区麻豆| 免费观看a级毛片全部| 久久精品亚洲av国产电影网| 国产成人av教育| 9色porny在线观看| 97人妻天天添夜夜摸| 在线观看www视频免费| 久久精品熟女亚洲av麻豆精品| 欧美激情极品国产一区二区三区| 在线观看免费日韩欧美大片| 在线观看日韩欧美| 国产高清videossex| 人妻久久中文字幕网| 搡老熟女国产l中国老女人| 欧美日韩亚洲高清精品| 不卡一级毛片| 欧美黑人欧美精品刺激| 亚洲精华国产精华精| 美女 人体艺术 gogo| 久久久久精品国产欧美久久久| 69av精品久久久久久| 国产精华一区二区三区| 91麻豆av在线| 久久久久久久久免费视频了| 欧美国产精品va在线观看不卡| 新久久久久国产一级毛片| 日韩精品免费视频一区二区三区| 亚洲人成77777在线视频| 日韩大码丰满熟妇| 久99久视频精品免费| xxx96com| 一级a爱视频在线免费观看| 欧美成人午夜精品| 亚洲中文字幕日韩| tube8黄色片| 黄片小视频在线播放| 另类亚洲欧美激情| 露出奶头的视频| 亚洲精品美女久久av网站| 欧美精品人与动牲交sv欧美| 日韩一卡2卡3卡4卡2021年| 黄色怎么调成土黄色| 欧美国产精品va在线观看不卡| 1024香蕉在线观看| 欧美精品一区二区免费开放| 亚洲人成伊人成综合网2020| 一区二区日韩欧美中文字幕| 91麻豆精品激情在线观看国产 | 在线观看免费视频网站a站| 国产精品免费一区二区三区在线 | 多毛熟女@视频| 一区二区三区国产精品乱码| 自线自在国产av| 国产高清激情床上av| 老司机影院毛片| av福利片在线| 免费观看精品视频网站| 不卡一级毛片| 精品国产美女av久久久久小说| 天堂中文最新版在线下载| 最近最新免费中文字幕在线| 久久亚洲精品不卡| 51午夜福利影视在线观看| 操出白浆在线播放| 国产成人av激情在线播放| 女人爽到高潮嗷嗷叫在线视频| 国产无遮挡羞羞视频在线观看| 国产av一区二区精品久久| 亚洲七黄色美女视频| 99热网站在线观看| 一级毛片高清免费大全| 久久久久久免费高清国产稀缺| 午夜福利一区二区在线看| 最近最新中文字幕大全免费视频| 国产区一区二久久| 久久精品国产综合久久久| 热re99久久精品国产66热6| 亚洲第一青青草原| 欧美亚洲 丝袜 人妻 在线| 午夜成年电影在线免费观看| 色婷婷久久久亚洲欧美| 国产精品久久久av美女十八| 欧美日韩亚洲综合一区二区三区_| 国产成人免费观看mmmm| 午夜两性在线视频| 亚洲aⅴ乱码一区二区在线播放 | av不卡在线播放| 精品乱码久久久久久99久播| 国产男女内射视频| 国产野战对白在线观看| 日韩大码丰满熟妇| 欧美日韩av久久| 久久香蕉国产精品| 黑人欧美特级aaaaaa片| 午夜福利一区二区在线看| 老司机午夜十八禁免费视频| 悠悠久久av| 免费人成视频x8x8入口观看| 午夜91福利影院| 黄片小视频在线播放| e午夜精品久久久久久久| netflix在线观看网站| 欧美在线黄色| 在线观看免费午夜福利视频| 建设人人有责人人尽责人人享有的| 亚洲人成77777在线视频| 夜夜夜夜夜久久久久| 精品卡一卡二卡四卡免费| 国产蜜桃级精品一区二区三区 | 黄色成人免费大全| 国内毛片毛片毛片毛片毛片| 一级毛片精品| 一本一本久久a久久精品综合妖精| 中文字幕色久视频| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 色综合婷婷激情| 久久精品国产清高在天天线| 校园春色视频在线观看| 女人高潮潮喷娇喘18禁视频| 一区二区三区精品91| av不卡在线播放| 久久香蕉国产精品| 50天的宝宝边吃奶边哭怎么回事| 人妻 亚洲 视频| 亚洲人成伊人成综合网2020| 少妇粗大呻吟视频| 黄片大片在线免费观看| 下体分泌物呈黄色| 中文字幕另类日韩欧美亚洲嫩草| 免费观看精品视频网站| 一本一本久久a久久精品综合妖精| 精品午夜福利视频在线观看一区| 极品少妇高潮喷水抽搐| 操出白浆在线播放| 久久婷婷成人综合色麻豆| 最近最新中文字幕大全免费视频| 国产在视频线精品| 久久久国产精品麻豆| 国产人伦9x9x在线观看| 久久久精品免费免费高清| 日本vs欧美在线观看视频| 999精品在线视频| av超薄肉色丝袜交足视频| 欧美乱码精品一区二区三区| 岛国毛片在线播放| 怎么达到女性高潮| 人人妻人人爽人人添夜夜欢视频| 中文字幕最新亚洲高清| 国产免费av片在线观看野外av| 女人被狂操c到高潮| 中文字幕制服av| 高清毛片免费观看视频网站 | 热re99久久精品国产66热6| 人成视频在线观看免费观看| 国产精品久久久久成人av| 女人久久www免费人成看片| 久久热在线av| 久久精品国产亚洲av香蕉五月 | 精品国产美女av久久久久小说| 久久久国产成人精品二区 | 色94色欧美一区二区| 成人18禁高潮啪啪吃奶动态图| 国产xxxxx性猛交| 国产麻豆69| 建设人人有责人人尽责人人享有的| 男女下面插进去视频免费观看| 91在线观看av| 亚洲av日韩精品久久久久久密| 久久久久久久久免费视频了| 精品乱码久久久久久99久播| 久久影院123| 日韩 欧美 亚洲 中文字幕| e午夜精品久久久久久久| 老熟妇仑乱视频hdxx| 亚洲一区二区三区欧美精品| 老司机在亚洲福利影院| 国产成人精品无人区| 精品第一国产精品| av免费在线观看网站| 人妻丰满熟妇av一区二区三区 | 久久久国产欧美日韩av| 久久久久国内视频| 国产精品电影一区二区三区 | 午夜精品久久久久久毛片777| 天天躁日日躁夜夜躁夜夜| 一进一出好大好爽视频| 99精品久久久久人妻精品| 久久精品亚洲熟妇少妇任你| 午夜福利视频在线观看免费| 1024视频免费在线观看| 窝窝影院91人妻| 丁香欧美五月| 亚洲av日韩精品久久久久久密| 高清av免费在线| 999精品在线视频| 国产成人精品久久二区二区91| a级毛片黄视频| 宅男免费午夜| 亚洲一区高清亚洲精品| 欧美亚洲日本最大视频资源| 久久99一区二区三区| 丰满迷人的少妇在线观看| 亚洲精品在线美女| 久久人妻福利社区极品人妻图片| 人人妻人人添人人爽欧美一区卜| 色在线成人网| 一级片免费观看大全| 脱女人内裤的视频| 十八禁高潮呻吟视频| 一级,二级,三级黄色视频| 国产在线精品亚洲第一网站| 人人妻,人人澡人人爽秒播| 天堂中文最新版在线下载| 久久人妻av系列| 久久人人97超碰香蕉20202| 高清黄色对白视频在线免费看| 久久久久久久久免费视频了| 黑人巨大精品欧美一区二区蜜桃| 久久亚洲真实| 国产成人免费无遮挡视频| 日本欧美视频一区| 高清黄色对白视频在线免费看| 操美女的视频在线观看| 在线观看66精品国产| 日韩三级视频一区二区三区| 欧美成人免费av一区二区三区 | 成人亚洲精品一区在线观看| 成熟少妇高潮喷水视频| 精品一区二区三区视频在线观看免费 | 亚洲专区国产一区二区| 老司机影院毛片| 精品久久久久久久久久免费视频 | 午夜免费成人在线视频| 老鸭窝网址在线观看| 亚洲av片天天在线观看| 亚洲久久久国产精品| 亚洲精品国产一区二区精华液| 精品国产超薄肉色丝袜足j| 久久中文字幕人妻熟女| 怎么达到女性高潮| 国产欧美日韩一区二区三| 在线观看日韩欧美| 欧美国产精品va在线观看不卡| 午夜福利乱码中文字幕| 久9热在线精品视频| 中文字幕高清在线视频| 少妇被粗大的猛进出69影院| 桃红色精品国产亚洲av| 国产精品一区二区免费欧美| 国产精品影院久久| 国产精品久久视频播放| 国精品久久久久久国模美| av线在线观看网站| 狠狠狠狠99中文字幕| 夫妻午夜视频| 老熟妇乱子伦视频在线观看| 日本黄色日本黄色录像| 久久精品国产99精品国产亚洲性色 | 美女视频免费永久观看网站| 亚洲人成77777在线视频| 最近最新中文字幕大全电影3 | 十八禁网站免费在线| 亚洲 欧美一区二区三区| 国产一区二区三区综合在线观看| 丝袜美足系列| 中文字幕人妻丝袜制服| 成年版毛片免费区| 一级片免费观看大全| 亚洲五月色婷婷综合| 高清毛片免费观看视频网站 | 国产亚洲精品第一综合不卡| 夜夜夜夜夜久久久久| 国产精品久久久av美女十八| 国产精品永久免费网站| 另类亚洲欧美激情| xxxhd国产人妻xxx| 90打野战视频偷拍视频| 超碰成人久久| 最近最新免费中文字幕在线| 不卡一级毛片| x7x7x7水蜜桃| tube8黄色片| 熟女少妇亚洲综合色aaa.| 久久久精品免费免费高清| 制服人妻中文乱码| 大香蕉久久网| 好看av亚洲va欧美ⅴa在| 久久人妻熟女aⅴ| 老司机在亚洲福利影院| 欧美精品一区二区免费开放| 久久人人97超碰香蕉20202| 两人在一起打扑克的视频| 午夜影院日韩av| e午夜精品久久久久久久| 巨乳人妻的诱惑在线观看| 国产精品欧美亚洲77777| 亚洲视频免费观看视频| 国产激情欧美一区二区| 成人免费观看视频高清| 久久精品亚洲精品国产色婷小说| 日韩人妻精品一区2区三区| 一区二区三区激情视频| 热99re8久久精品国产| 男女高潮啪啪啪动态图| 美女高潮喷水抽搐中文字幕| 亚洲成人手机| 亚洲一码二码三码区别大吗| 少妇 在线观看| 亚洲五月婷婷丁香| 亚洲专区中文字幕在线| 久久草成人影院| 国产蜜桃级精品一区二区三区 | 欧美日韩亚洲综合一区二区三区_| 欧美日韩瑟瑟在线播放| 国产精华一区二区三区| 91字幕亚洲| 极品人妻少妇av视频| 免费少妇av软件| 国产片内射在线| 高清视频免费观看一区二区| 亚洲av成人一区二区三| cao死你这个sao货| 天堂俺去俺来也www色官网| 一进一出抽搐动态| 国产真人三级小视频在线观看| 欧美国产精品一级二级三级| 欧美日本中文国产一区发布| 亚洲一区二区三区不卡视频| 老鸭窝网址在线观看| 亚洲精品国产区一区二| 国产亚洲欧美98| 午夜福利欧美成人| 我的亚洲天堂| 日韩视频一区二区在线观看| 国产精品国产高清国产av | www.自偷自拍.com| 欧美日韩亚洲国产一区二区在线观看 | 最近最新免费中文字幕在线| 12—13女人毛片做爰片一| 国产区一区二久久| 热99久久久久精品小说推荐| 日日摸夜夜添夜夜添小说| 狠狠婷婷综合久久久久久88av| 成年版毛片免费区| 久9热在线精品视频| 国产亚洲精品第一综合不卡| 精品熟女少妇八av免费久了| 欧美乱码精品一区二区三区| 久久精品人人爽人人爽视色| 久久 成人 亚洲| 成年人免费黄色播放视频| 高潮久久久久久久久久久不卡| 精品乱码久久久久久99久播| 最新美女视频免费是黄的| 在线免费观看的www视频| 国产高清国产精品国产三级| 丰满的人妻完整版| 国产高清国产精品国产三级| 好看av亚洲va欧美ⅴa在| 久久国产精品男人的天堂亚洲| 女人精品久久久久毛片| 亚洲专区中文字幕在线| 日韩有码中文字幕| 欧美+亚洲+日韩+国产| 一级作爱视频免费观看| 国产欧美日韩一区二区三区在线| 99久久国产精品久久久| 亚洲少妇的诱惑av| 免费观看a级毛片全部| 免费在线观看视频国产中文字幕亚洲| 极品教师在线免费播放| 久久 成人 亚洲| 国产一区有黄有色的免费视频| 精品熟女少妇八av免费久了| 久久久久久久国产电影| 午夜老司机福利片| 人人妻,人人澡人人爽秒播| 欧美人与性动交α欧美精品济南到| 999精品在线视频| 国产av又大| 丝袜在线中文字幕| 国产欧美日韩综合在线一区二区| 欧美亚洲日本最大视频资源| 久9热在线精品视频| 久久久久久久国产电影| 国产国语露脸激情在线看| 日韩制服丝袜自拍偷拍| 久久久精品国产亚洲av高清涩受| 国产精品影院久久| 女人精品久久久久毛片| 色在线成人网| 亚洲第一欧美日韩一区二区三区| 精品亚洲成a人片在线观看| 精品一区二区三区四区五区乱码| 91字幕亚洲| 美女 人体艺术 gogo| 老熟女久久久| 人人妻人人澡人人看| 国产在线一区二区三区精| 99热网站在线观看| 免费黄频网站在线观看国产| 久久久久国产一级毛片高清牌| 国产又色又爽无遮挡免费看| 午夜久久久在线观看| 少妇被粗大的猛进出69影院| 伦理电影免费视频| 日韩欧美三级三区| 美女福利国产在线| 欧美日韩中文字幕国产精品一区二区三区 | 淫妇啪啪啪对白视频| 久久精品国产清高在天天线| 国产成人av激情在线播放| 啦啦啦 在线观看视频| 少妇猛男粗大的猛烈进出视频| 亚洲成av片中文字幕在线观看| 亚洲精品国产色婷婷电影| 久久天堂一区二区三区四区| 国产aⅴ精品一区二区三区波| av视频免费观看在线观看| 黄色丝袜av网址大全| 高清视频免费观看一区二区| 国产激情欧美一区二区| 成人亚洲精品一区在线观看| 美国免费a级毛片| 满18在线观看网站| 久热爱精品视频在线9| 国产成人影院久久av| 999精品在线视频| 日本一区二区免费在线视频| 亚洲少妇的诱惑av| 精品午夜福利视频在线观看一区| 一区二区日韩欧美中文字幕| 久久热在线av| 国产一区二区激情短视频| 精品久久久精品久久久| 满18在线观看网站| 侵犯人妻中文字幕一二三四区| 两性夫妻黄色片| 久久人人97超碰香蕉20202| 欧美黄色淫秽网站| 久久精品熟女亚洲av麻豆精品| 国产高清国产精品国产三级| 亚洲精品美女久久久久99蜜臀| 国产成人免费无遮挡视频| 欧美激情高清一区二区三区| a级毛片黄视频| 一区二区三区精品91| 嫩草影视91久久| 高清在线国产一区| 久久人妻福利社区极品人妻图片| 老司机靠b影院| 欧美国产精品一级二级三级| 黄色视频,在线免费观看| 国产99白浆流出| 亚洲精品国产一区二区精华液| 一个人免费在线观看的高清视频| 天堂俺去俺来也www色官网| 精品午夜福利视频在线观看一区| 成人免费观看视频高清| 亚洲中文字幕日韩| 精品午夜福利视频在线观看一区| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲人成电影观看| 免费少妇av软件| 男男h啪啪无遮挡| 久久性视频一级片| 国产真人三级小视频在线观看| 一级毛片高清免费大全| 嫁个100分男人电影在线观看| 亚洲国产毛片av蜜桃av| 免费一级毛片在线播放高清视频 | 男女下面插进去视频免费观看| 老熟女久久久| 50天的宝宝边吃奶边哭怎么回事| 欧美精品人与动牲交sv欧美| videos熟女内射| 一级毛片精品| 嫁个100分男人电影在线观看| 午夜福利欧美成人| 欧美日韩国产mv在线观看视频| 久久久国产成人免费|