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

    井-井電阻率2.5維近似成像研究

    2016-08-05 08:56:53呂玉增
    物探化探計(jì)算技術(shù) 2016年3期
    關(guān)鍵詞:井間電阻率反演

    熊 飛 , 呂玉增, 李 勇

    (1.華安奧特(北京)科技股份有限公司,北京 100085;2.桂林理工大學(xué) 地球科學(xué)學(xué)院,桂林 541004;3.廣西礦冶與環(huán)境科學(xué)實(shí)驗(yàn)中心,桂林 541004)

    ?

    井-井電阻率2.5維近似成像研究

    熊飛1, 呂玉增2,3, 李勇2

    (1.華安奧特(北京)科技股份有限公司,北京100085;2.桂林理工大學(xué)地球科學(xué)學(xué)院,桂林541004;3.廣西礦冶與環(huán)境科學(xué)實(shí)驗(yàn)中心,桂林541004)

    摘要:在使用有限單元法解決了井-井2.5維正演問(wèn)題后,采用矩形剖分方式將井間剖面離散成若干個(gè)點(diǎn),再將正演所得的視電阻率數(shù)據(jù)采用一種新穎的方式分配重組到井間每個(gè)離散點(diǎn)上,繪制出了井間剖面電阻率近似成像圖。該近似成像方法避開(kāi)了反演成像中龐大的內(nèi)存需求和計(jì)算量,計(jì)算效率高,并且成像結(jié)果表明該方法具有較好的效果,是一種行之有效的井間視電阻率資料處理手段。

    關(guān)鍵詞:井-井電阻率法; 2.5維正演; 有限元; 近似成像

    0引言

    井-井電阻率法將供電點(diǎn)和測(cè)量點(diǎn)都置于鉆孔中,更加接近目標(biāo)體,因而在深部盲礦體的勘探方面相對(duì)地面電法更為有效。當(dāng)下現(xiàn)有的淺部礦產(chǎn)資源日益枯竭,急需對(duì)深部隱伏礦床進(jìn)行勘探,使得井-井電阻率成像研究具有一定的現(xiàn)實(shí)意義。

    在上世紀(jì)90年代,國(guó)、內(nèi)外就開(kāi)始了井間電阻率成像的研究,Daily等[1]首先將有限單元法應(yīng)用于井間二維電阻率成像中;Shima[2]又將有限單元法與α中心法聯(lián)合用于井間電阻率圖像的重建中;董清華[3]系統(tǒng)介紹了井間電阻率層析成像的一些進(jìn)展;王興泰等[4]使用佐迪反演進(jìn)行了電阻率成像研究;底青云等[5]采用積分法實(shí)現(xiàn)了電阻率三維成像;Zhdanov等[6]研究了井間三維成像;Wilkinson等[7]和Chambers等[8]實(shí)現(xiàn)了礦井電阻率成像;A.Arato等[9]實(shí)現(xiàn)了基于交錯(cuò)網(wǎng)格的二維井間電阻率層析成像,提高了井間成像效果。

    反演成像過(guò)程往往要涉及到大量的計(jì)算,使得成像時(shí)間增加,制約了算法的實(shí)用性,在提高成像速度方面,人們做了大量研究,取得了諸多成果:董清華[10]、劉國(guó)強(qiáng)等[11]對(duì)井間電阻率反演成像中偏導(dǎo)數(shù)的計(jì)算方法進(jìn)行了對(duì)比研究;Knoll[12]、吳小平等[13]采用直接計(jì)算Jacobian矩陣向量積的方法,避開(kāi)了Jacobian矩陣的計(jì)算和存儲(chǔ);呂玉增[14]、沈平等[15]直接避開(kāi)反演方法探討了井間電阻率的近似成像。

    這里借鑒前人的思路,從點(diǎn)電流場(chǎng)在地下介質(zhì)的分布情況以及地下異常球體的電位解析式出發(fā),提出了一種新穎的井-井電阻率近似成像方案。算例表明,該方法只需2 s左右便可完成一次近似成像,并且具有較好的成像效果。

    1井-井2.5維有限元數(shù)值模擬

    對(duì)于復(fù)雜的地下介質(zhì)分布,有限單元法是一種強(qiáng)有力的數(shù)值模擬方法[16],因而這里的2.5維正演問(wèn)題采用有限單元法解決??紤]到電位在電源點(diǎn)附近具有奇異性,如果直接采用總電位法求解,包含電源項(xiàng),勢(shì)必會(huì)導(dǎo)致求解的電源附近節(jié)點(diǎn)處的電位值具有較大誤差,為了避免這一問(wèn)題,作者采用異常電位法求解。

    1.1二極觀測(cè)

    井-井觀測(cè)一般采用在兩口不同的井中布設(shè)供電電極和測(cè)量電極的工作方式,主要包括二極裝置、三極裝置和四極裝置等。這里近似成像研究基于二極裝置。它將供電電極A置于供電井中,供電電極B置于地面無(wú)窮遠(yuǎn)處;將測(cè)量電極M置于測(cè)量井中,測(cè)量電極N置于無(wú)窮遠(yuǎn)處,通過(guò)移動(dòng)供電井中的供電電極A實(shí)現(xiàn)全孔供電,并且在每一個(gè)供電深度處,都移動(dòng)測(cè)量井中的測(cè)量電極M進(jìn)行全孔測(cè)量(圖1)。與地面二極觀測(cè)不同,井-井二極觀測(cè)中視電阻率計(jì)算公式定義如下:

    (1)

    其中:UM為M處的測(cè)量電位;KAM為裝置系數(shù)。

    圖1 井-井二極觀測(cè)示意圖Fig.1 Sketch of borehole-borehole pole-pole survey

    1.2點(diǎn)源異常電位的變分問(wèn)題及有限元正演

    與地面電法不同,井-井觀測(cè)供電點(diǎn)置于地下半空間中,因而在正常電位的表達(dá)式以及無(wú)窮遠(yuǎn)邊界條件的推導(dǎo)中,都應(yīng)設(shè)置關(guān)于地表對(duì)稱(chēng)的電源鏡像點(diǎn)。經(jīng)過(guò)推導(dǎo),得到井-井2.5維點(diǎn)源異常電位變分問(wèn)題:

    (2)

    其中:▽為哈密頓算子;σ為地下介質(zhì)電導(dǎo)率;σ′為介質(zhì)異常電導(dǎo)率;k為傅里葉變換中的波數(shù);U0為正常傅氏電位;U為待求的異常傅氏電位;Ω為計(jì)算區(qū)域;Γ∞為無(wú)窮遠(yuǎn)邊界;r1、r2分別為點(diǎn)電源A與其鏡像點(diǎn)A*至無(wú)窮遠(yuǎn)邊界?!薜木嚯x;n為?!尥夥ň€方向;K0、K1分別為第二類(lèi)零階和一階修正貝塞爾函數(shù)。

    同樣根據(jù)鏡像法原理,正常傅氏電位U0的表達(dá)式可寫(xiě)為式(3)。

    (3)

    其中:I為點(diǎn)源電流強(qiáng)度;σ0為介質(zhì)正常電導(dǎo)率;r為待求點(diǎn)與電源點(diǎn)之間的距離;r′為待求點(diǎn)與電源鏡像點(diǎn)之間的距離。

    采用以井間區(qū)域?yàn)橹行闹鸩较蛲庀∈璧钠史址绞?,將地下通過(guò)兩井的剖面剖分成若干個(gè)矩形單元,在每個(gè)單元內(nèi)進(jìn)行雙線性插值,將式(2)中F(U)分解為各個(gè)矩形單元上的積分,再將各單元積分相加,可得式(4)。

    (4)

    根據(jù)F(U)取極值的條件有:

    KU=-K′U0

    (5)

    其中:U、U0分別為所有節(jié)點(diǎn)處U、U0組成的列向量;系數(shù)矩陣K、K′為大型稀疏對(duì)稱(chēng)矩陣。

    為了節(jié)省內(nèi)存,采用改進(jìn)的行壓縮存儲(chǔ)方法[17](Modified Sparse Row ,MSR)來(lái)存儲(chǔ)系數(shù)矩陣?yán)锏姆橇阍?,為了提高?jì)算精度和效率,采用對(duì)稱(chēng)逐步超松弛預(yù)條件共軛梯度法[18](SSORPCG)法求解線性方程組。最后將求得的異常傅氏電位通過(guò)傅里葉反變換(表1為反變換中采用的一組波數(shù)k和系數(shù)g[16]),可得到各節(jié)點(diǎn)處的異常電位,與正常電位相加便得到了各節(jié)點(diǎn)處的總電位,再將電位換算成視電阻率。

    表1 傅里葉反變換的k和g

    2井-井電阻率近似成像

    基于上述正演算法,可以模擬得到井間不同模型二極觀測(cè)下的視電阻率數(shù)據(jù),再以這些數(shù)據(jù)為基礎(chǔ),依照如下方式,實(shí)現(xiàn)井-井2.5維電阻率近似成像。

    2.1近似成像原理

    根據(jù)地下點(diǎn)電源電場(chǎng)的分布規(guī)律以及視電阻率的定義,在井-井二極觀測(cè)中,測(cè)量點(diǎn)M處所得到的視電阻率值是供電點(diǎn)A和測(cè)量點(diǎn)M之間介質(zhì)的綜合反映,并且主要取決于AM連線附近一小塊鄰域內(nèi)介質(zhì)電阻率。

    從異常球體的電位公式[19]可以發(fā)現(xiàn):COS∠AOM直接影響著解析解(A為供電點(diǎn)、O為球心、M為測(cè)量點(diǎn))。由于井間剖面可視為無(wú)數(shù)個(gè)離散的點(diǎn)P,每個(gè)點(diǎn)又可以近似當(dāng)做一個(gè)“微小球體”,那么可以用COS∠APM的大小來(lái)限定AM連線鄰域,并近似作為該A-M觀測(cè)下視電阻率反映區(qū)域。如圖2所示,對(duì)于任意一條連線AM,存在無(wú)數(shù)個(gè)近似橢圓的區(qū)域,記區(qū)域邊界上任一點(diǎn)(A、M除外)為P*、區(qū)域內(nèi)的任意一點(diǎn)為P,使得∠AP*M為一定值,且COS∠APM

    (6)

    其中:S=AM/d。

    圖2 AM鄰域限定示意圖Fig.2 Sketch of AM neighbourhood limited

    由式(6)可以看出,參數(shù)S很好地確定了鄰域范圍;并且當(dāng)選定一個(gè)固定S值后,AM增大時(shí),d也會(huì)隨之增大,這與電流實(shí)際分布情況較為符合。所以作者提出的近似思路在理論上具有一定的合理性。

    基于以上思路,給出具體的井間視電阻率近似成像操作方案:

    1)采用矩形剖分的方式將井間剖面近似離散成若干個(gè)點(diǎn),橫向剖分?jǐn)?shù)為NX,縱向剖分?jǐn)?shù)為NZ(圖3)。

    圖3 近似成像網(wǎng)格化示意圖Fig.3 Sketch of meshing for approximate image

    2)對(duì)于一個(gè)格點(diǎn)P,依次試探每一組A-M,求出相應(yīng)的COS∠APM,并與給定值MAXCOS進(jìn)行比較,若COS∠APM≤MAXCOS,那么記P點(diǎn)處電阻率為該A-M觀測(cè)下的視電阻率值,否則記為“0”,將所有的非零值的平均視為P點(diǎn)處電阻率值;如此得到所有點(diǎn)處的電阻率值,便可繪制近似成像圖。其中,式(6)為MAXCOS的計(jì)算式,COS∠APM的計(jì)算式如(7)所示。

    (7)

    2.2井間模型的近似成像

    通過(guò)對(duì)單一簡(jiǎn)單模型進(jìn)行試算表明,當(dāng)S取9時(shí),會(huì)得到相對(duì)較好的成像效果,NX、NZ對(duì)成像影響較小。采用上述近似成像方法,取S=9、NX=10、NZ=50,對(duì)四種有代表性的井間二維模型、模型組合進(jìn)行近似成像。其中,井深為H=100 m,井間距為D=20 m,背景電阻率為ρ0=100Ω·m。供電極A在供電井中移動(dòng)范圍為0 m~100 m,點(diǎn)距為2 m,共進(jìn)行51次供電,供電強(qiáng)度為1A;每次供電時(shí),測(cè)量極M在測(cè)量井中從0 m至100 m每隔2 m依次進(jìn)行51次測(cè)量。

    模型一為厚度為4m的傾斜高阻板狀體,ρ1=1 000 Ω·m,具體如圖4(a)所示。

    模型二為厚度為4 m的“L”形高阻板狀體,ρ1=1 000 Ω·m,具體如圖4(b)所示。

    圖4 井間模型示意圖Fig.4 Sketch of models between boreholes

    模型三為井間上方賦存一4 m×4 m的低阻異常體,ρ1=1 0 Ω·m,下方賦存12 m×4 m的高阻異常體,ρ2=1 000 Ω·m。具體如圖4(c)所示。

    模型四為在同一深度賦存兩個(gè)尺寸都為4 m×8 m的異常體,左邊的為低阻異常:ρ1=1 0 Ω·m右邊為高阻異常,ρ2=1 000 Ω·m,具體如圖4(d)所示。

    模型一至模型四的近似成像圖如圖5所示,圖5中正方形、長(zhǎng)方形等對(duì)應(yīng)異常體截面實(shí)際形態(tài)和位置。

    從圖5中可以看出,①當(dāng)異常體為傾斜板狀時(shí),異常區(qū)域顯示為一傾斜橢圓,橢圓位置、范圍以及傾斜方向和角度都與異常體實(shí)際對(duì)應(yīng)較好;②當(dāng)異常體為“L”形時(shí),異常區(qū)域近似為“芒果“形態(tài),大致上反映了異常體的位置以及在橫向和縱向上的展布;③當(dāng)兩個(gè)異常體一上一下,形態(tài)不同、電阻率也不同時(shí),近似成像圖無(wú)論在位置、形態(tài)還是電阻率上都能加以反映;④當(dāng)兩個(gè)形態(tài)相同、電阻率不同的異常體處于同一深度,并且一左一右相互干擾時(shí),近似成像圖仍能較明顯地辨別。

    在AMD Phenom(tm) II N660 Dual-Core ProcessorCPU,3 GHz,2 GB內(nèi)存,Windows 7平臺(tái)下,每次成像時(shí)間只需2 s左右。除此之外,還應(yīng)說(shuō)明以下兩點(diǎn):

    1)由于近似成像圖是根據(jù)視電阻率數(shù)據(jù)重組而來(lái),因而繼承著視電阻率數(shù)據(jù)本身的特點(diǎn)——高阻區(qū)域?qū)?yīng)低阻異常、低阻區(qū)域?qū)?yīng)高阻異常,這點(diǎn)與反演成像有所區(qū)別。

    2)近似成像圖上下邊界處存在明顯的形如“山峰”的假異常,“峰頂”指向異常體。這時(shí)由于成像過(guò)程中采用的近似方式引起的,在圖像解釋中應(yīng)注意區(qū)別。

    圖5 井間模型近似成像圖Fig.5 The approximate image of models between borehole

    3結(jié)語(yǔ)

    作者采用有限單元法實(shí)現(xiàn)了井-井電阻率2.5維正演模擬,其中采用MSR法存儲(chǔ)超大型稀疏系數(shù)矩陣,運(yùn)用SSORPCG法求解線性方程組,保證了正演計(jì)算的效率和精度。根據(jù)地下點(diǎn)源電場(chǎng)的分布規(guī)律,提出了一種新穎的井-井電阻率2.5維近似成像方法,據(jù)此進(jìn)行了四種典型模型的近似成像。結(jié)果顯示,這里的井間電阻率近似成像方法能定性判斷異常體電阻率、較準(zhǔn)確反映異常體的位置、形態(tài),并且具有極快的成像速度,因此可用于井-井觀測(cè)數(shù)據(jù)的快速預(yù)處理,得到井間剖面較為直觀的構(gòu)造。

    參考文獻(xiàn):

    [1]DAILY W, OWEN E. Cross-borehole resistivity tomography[J]. Geophysics, 1991, 56( 8): 1228-1235.

    [2]HIROMASA SHIMA. 2D and 3D resistivity image reconstruction using crosshole data[J]. Geophysics, 1992, 57( 10): 1270-1281.

    [3]董清華. 井間電阻率層析成像的某些進(jìn)展[J]. 地球物理學(xué)進(jìn)展, 1997, 12(3): 77-89.

    DONG Q H. The latest developments of crosshole resistivity tomography[J]. Progress in Geophysics, 1997, 12(3): 77-89. (In Chinese)

    [4]王興泰, 李曉芹. 電阻率圖像重建的佐迪反演及其應(yīng)用效果[J]. 物探與化探, 1996, 20(3): 228-233.

    WANG X T, LI X Q. The image reconstruction of electric resistivity by Zohdy inversion and it’s application effects[J]. Geophysical&Geochemical Exploration, 1996, 20(3): 228-233. (In Chinese)

    [5]底青云, 王妙月. 積分法三維電阻率成像[J]. 地球物理學(xué)報(bào), 2001, 44(6): 844-851.

    DI Q Y, WANG M Y. 3-D resistivity tomography by integral method[J]. Chinese J. Geophys,2001, 44(6): 844-851. (In Chinese)

    [6] ZHDANOW M S, YOSHIOKA K. Cross-well electromagnetic imaging in three dimesions[J]. Exploration Geophysics, 2003,34:34-40.

    [7]WILKINSON P B, CHANBERS J E,LELLIOTT M,et al. Extreme sensitivity of crosshole electrical resistivity tomography measurements to geometric errers[J]. Geophys, 2008,173:49-62.

    [8]CHAMBERS J E, WILKINSON P B,WELLER AL,et al. Mineshaft imaging using surface and crosshole 3D electrical resistivity tomography :A case history from the East Pennine Coalfield[J]. UK. Journal of Applied-Geophysics, 2007,62(4):324-337.

    [9]A. ARATO, A. GODIO,L.SAMBUELLI. Staggered grid inversion of cross hole 2-D resistivity tomography[J]. Journal of Applied-Geophysics, 2014,107:60-70.

    [10]董清華, 田憲謨. 井間電阻率成像中Frchet導(dǎo)數(shù)的算法比較[J]. 物探化探計(jì)算技術(shù), 1999, 19(1): 41-45.

    DONG Q H, TIAN X M. A comparison of methods for calculating Frchet derivatives in crosshole resistivity image[J]. Computing Techniques for Geophysical and Geochemical Exploration, 1999, 19(1): 41-45. (In Chinese)

    [11]劉國(guó)強(qiáng), 陶果, 柯式鎮(zhèn), 等. 攝動(dòng)求解Frchet導(dǎo)數(shù)的電阻率層析成像方法[J]. 地質(zhì)地震, 2001, 23(2): 314-320.

    LIU G Q, TAO G, KE S Z, et al. Improved reconstruction algorithms for resistivity tomography[J]. Seismology and Geology, 2001, 23(2): 314-320. (In Chinese)

    [12]KNOLL D A, KEYES D E. Jacobian-free Newton-Krylov methods: a survey of approaches and applications[J]. Journal of Computational Physics, 2004, 193(2): 357-397.

    [13]吳小平, 徐果明. 電阻率三維反演中偏導(dǎo)數(shù)矩陣的求取和分析[J]. 石油地球物理勘探, 1999, 34(4): 363-372.

    WU X P, XU G M. Derivation and analysis of partial derivative matrix in resistivity 3-D inversion[J]. Oil Geophysical Prospecting, 1999, 34(4): 363-372. (In Chinese)

    [14]呂玉增, 阮百堯, 黃俊革. 直流電井間三維直接成像[J]. 物探化探計(jì)算技術(shù), 2003, 25(1): 60-64.

    LV Y Z, RUAN B Y, HUANG J G. The 3-D immediate crosshole tomography with direct current[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2003, 25(1): 60-64. (In Chinese)

    [15]沈平, 強(qiáng)建科, 李永軍,等. 井間視電阻率的幾何成像方法[J]. 中南大學(xué)學(xué)報(bào):自然科學(xué)版, 2010, 41(3): 1079-1084.

    SHEN P. QIANG J K, LI Y J, et al. Geometry image method of crosshole apparent resistivity[J]. Journal of Central South University (Science and Technology), 2010, 41(3): 1079-1084. (In Chinese)

    [16]徐世浙. 地球物理中的有限單元法[M]. 北京: 科學(xué)出版社, 1994.

    XU S Z. The finite element method in Geophysics[M]. Beijing: Science Press, 1994.(In Chinese)

    [17]SAAD Y. Iterative methods for sparse linear systems[M]. Siam, 2003.

    [18]吳小平, 汪彤彤. 利用共軛梯度算法的電阻率三維有限元正演[J]. 地球物理學(xué)報(bào), 2004,46(3): 428-432.

    WU X P, WANG T T. A 3-D finite-element resistivity forward modeling using conjugate gradient algorithm[J]. Chinese J. Geophys, 2004, 46(3): 428-432. (In Chinese)

    [19]李金銘. 地電場(chǎng)與電法勘探[M]. 北京: 地質(zhì)出版社, 2005.

    LI J M. Geoelectric field and electrical exploration[M]. Beijing: Geological Publishing House, 2005. (In Chinese)

    收稿日期:2015-04-10改回日期:2015-05-21

    基金項(xiàng)目:國(guó)家“十二五”國(guó)家重大科學(xué)儀器設(shè)備開(kāi)發(fā)專(zhuān)項(xiàng)(2012YQ030126);廣西自然科學(xué)基金(2013GXNSFBA019212);廣西高??茖W(xué)技術(shù)研究重點(diǎn)項(xiàng)目(2013ZD029)

    作者簡(jiǎn)介:熊飛(1990-)男,碩士,主要從事電磁法數(shù)值模擬研究,E-mail:xiongfei357@163.com。

    文章編號(hào):1001-1749(2016)03-0308-06

    中圖分類(lèi)號(hào):P 631.3

    文獻(xiàn)標(biāo)志碼:A

    DOI:10.3969/j.issn.1001-1749.2016.03.03

    A study on approximate imaging of borehole-borehole 2.5D resistivity method

    XIONG Fei1, LV Yu-zeng2,3, LI Yong2

    (1.HUAAN AOTE(Beijing) Technology Co.,Ltd,Beijing100085,China;2.Guilin University of Technology,College of Earth Science,Guilin541004,China;3.Guangxi Scientific Experiment Center of Mining,Metallurgy and Environment,Guilin541004,China)

    Abstract:After solving the borehole-borehole 2.5 dimensional forward problem by using the finite element method, we can use the rectangular partition discretized the section between boreholes into a number of points. According to the distribution of the electric field of a point source in the underground space, the apparent resistivity data can be re-allocated on each discrete points between boreholes through a novel way. Then we can draw the approximate resistivity image across hole profile. The approximate imaging method has high computational efficiency as avoiding the large memory requirements and calculations. And the imaging results show that the method is a kind of effective processing method of apparent resistivity data between boreholes because of it’s good effect.

    Key words:borehole-borehole resistivity method; 2.5D forward; finite element method; approximate imaging

    猜你喜歡
    井間電阻率反演
    反演對(duì)稱(chēng)變換在解決平面幾何問(wèn)題中的應(yīng)用
    煤層氣井間抽機(jī)理及故障處理方法研究及應(yīng)用
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    井間示蹤劑監(jiān)測(cè)在復(fù)雜斷塊油藏描述中的應(yīng)用
    錄井工程(2017年1期)2017-07-31 17:44:42
    三維電阻率成像與高聚物注漿在水閘加固中的應(yīng)用
    隨鉆電阻率測(cè)井的固定探測(cè)深度合成方法
    海洋可控源電磁場(chǎng)視電阻率計(jì)算方法
    疊前同步反演在港中油田的應(yīng)用
    粉煤灰摻量對(duì)水泥漿體電阻率與自收縮的影響
    亚洲欧美一区二区三区黑人 | 老司机亚洲免费影院| 国精品久久久久久国模美| 97超视频在线观看视频| 爱豆传媒免费全集在线观看| 国产欧美日韩综合在线一区二区| 国产欧美日韩综合在线一区二区| 黑人高潮一二区| 中文字幕最新亚洲高清| 国产视频首页在线观看| 视频中文字幕在线观看| 91在线精品国自产拍蜜月| 一区二区三区四区激情视频| 国产精品久久久久久av不卡| 91久久精品国产一区二区成人| 91精品国产国语对白视频| 国产精品久久久久久av不卡| 老司机亚洲免费影院| 久久人人爽人人片av| 狂野欧美激情性bbbbbb| 国产 精品1| 蜜桃国产av成人99| 九九爱精品视频在线观看| av黄色大香蕉| 免费黄网站久久成人精品| 国产成人精品一,二区| 美女中出高潮动态图| 80岁老熟妇乱子伦牲交| 国产又色又爽无遮挡免| 成人午夜精彩视频在线观看| kizo精华| 日韩,欧美,国产一区二区三区| 看免费成人av毛片| 国产精品久久久久久久电影| 亚洲av日韩在线播放| 汤姆久久久久久久影院中文字幕| 免费高清在线观看视频在线观看| 免费播放大片免费观看视频在线观看| 最黄视频免费看| 少妇被粗大猛烈的视频| 99久久中文字幕三级久久日本| av免费在线看不卡| 亚洲内射少妇av| 亚洲欧美日韩卡通动漫| 最近手机中文字幕大全| 51国产日韩欧美| 少妇熟女欧美另类| 性色av一级| 老熟女久久久| 久久国产精品男人的天堂亚洲 | 亚洲人成网站在线观看播放| 少妇被粗大的猛进出69影院 | 中文字幕人妻熟人妻熟丝袜美| 中文字幕人妻熟人妻熟丝袜美| 高清在线视频一区二区三区| 国产在线一区二区三区精| 亚洲精品国产av蜜桃| 欧美精品亚洲一区二区| 亚洲国产色片| 亚洲欧美日韩另类电影网站| 99久国产av精品国产电影| 啦啦啦中文免费视频观看日本| 久久97久久精品| 激情五月婷婷亚洲| 多毛熟女@视频| 日本猛色少妇xxxxx猛交久久| 色哟哟·www| 亚洲欧洲精品一区二区精品久久久 | 三上悠亚av全集在线观看| 九九爱精品视频在线观看| 新久久久久国产一级毛片| 少妇精品久久久久久久| 精品一品国产午夜福利视频| 亚洲欧美色中文字幕在线| 久久久国产欧美日韩av| 久久精品国产鲁丝片午夜精品| 丁香六月天网| 亚洲av成人精品一区久久| 国产亚洲精品第一综合不卡 | 99九九在线精品视频| av又黄又爽大尺度在线免费看| 色哟哟·www| 日韩成人av中文字幕在线观看| 国产精品久久久久久av不卡| 老司机亚洲免费影院| 麻豆精品久久久久久蜜桃| 免费看av在线观看网站| av黄色大香蕉| av视频免费观看在线观看| 在线精品无人区一区二区三| 伦理电影大哥的女人| 免费观看av网站的网址| 国产精品 国内视频| 亚洲熟女精品中文字幕| 精品午夜福利在线看| 在线观看国产h片| 国产成人av激情在线播放 | 亚洲精品国产av蜜桃| av在线app专区| 精品视频人人做人人爽| 午夜免费男女啪啪视频观看| 日本黄色日本黄色录像| 夫妻性生交免费视频一级片| 免费观看在线日韩| 欧美精品人与动牲交sv欧美| 在线观看免费日韩欧美大片 | 中文字幕免费在线视频6| 日日啪夜夜爽| 日日啪夜夜爽| 97超碰精品成人国产| 美女脱内裤让男人舔精品视频| 国产免费一级a男人的天堂| av在线观看视频网站免费| 国产国语露脸激情在线看| 五月玫瑰六月丁香| 国产男女内射视频| 中国国产av一级| 大片免费播放器 马上看| 欧美日韩av久久| 久久99热6这里只有精品| 91成人精品电影| 国产成人午夜福利电影在线观看| 国产精品.久久久| av.在线天堂| 欧美精品人与动牲交sv欧美| 人成视频在线观看免费观看| 69精品国产乱码久久久| 精品久久久久久久久亚洲| 亚州av有码| 少妇猛男粗大的猛烈进出视频| 亚洲人与动物交配视频| 亚洲欧美色中文字幕在线| 秋霞伦理黄片| 制服丝袜香蕉在线| 99热网站在线观看| 天堂8中文在线网| 久久影院123| 成年人午夜在线观看视频| 亚洲国产成人一精品久久久| 两个人免费观看高清视频| 全区人妻精品视频| 最黄视频免费看| 久久午夜综合久久蜜桃| 日韩亚洲欧美综合| 亚洲av不卡在线观看| 午夜精品国产一区二区电影| 中文字幕久久专区| 女性生殖器流出的白浆| 2022亚洲国产成人精品| 18禁裸乳无遮挡动漫免费视频| 日日摸夜夜添夜夜爱| 另类精品久久| 精品视频人人做人人爽| 欧美激情极品国产一区二区三区 | 久久99精品国语久久久| 一级片'在线观看视频| 王馨瑶露胸无遮挡在线观看| 三级国产精品片| 亚洲丝袜综合中文字幕| 在线天堂最新版资源| 亚洲色图 男人天堂 中文字幕 | 欧美性感艳星| 性色avwww在线观看| 在线亚洲精品国产二区图片欧美 | 国产精品免费大片| 黄色怎么调成土黄色| 午夜激情av网站| 国产日韩欧美亚洲二区| 欧美日本中文国产一区发布| 在线观看免费日韩欧美大片 | 涩涩av久久男人的天堂| 少妇高潮的动态图| 欧美日韩成人在线一区二区| 亚洲国产欧美在线一区| 啦啦啦视频在线资源免费观看| 简卡轻食公司| 97在线人人人人妻| 久久av网站| 狠狠婷婷综合久久久久久88av| 日本色播在线视频| 欧美国产精品一级二级三级| 一区二区三区免费毛片| 另类精品久久| 国产高清三级在线| 欧美日本中文国产一区发布| 国产精品国产三级专区第一集| 亚洲精品色激情综合| 热re99久久国产66热| 国产一级毛片在线| 国产又色又爽无遮挡免| 国产免费一级a男人的天堂| 一区二区av电影网| 色视频在线一区二区三区| 九色亚洲精品在线播放| 在线观看人妻少妇| 色网站视频免费| 精品久久久久久久久亚洲| 18禁在线播放成人免费| 免费日韩欧美在线观看| 久久国产精品男人的天堂亚洲 | 国产精品99久久久久久久久| 欧美日韩一区二区视频在线观看视频在线| 久久久久精品久久久久真实原创| 高清不卡的av网站| 三级国产精品片| 欧美人与善性xxx| 日本与韩国留学比较| 激情五月婷婷亚洲| 久久久久国产精品人妻一区二区| 欧美变态另类bdsm刘玥| 日本黄大片高清| 777米奇影视久久| 曰老女人黄片| 国产日韩欧美在线精品| 九九在线视频观看精品| 视频区图区小说| 亚洲欧美日韩卡通动漫| 亚洲国产欧美日韩在线播放| 考比视频在线观看| 国产免费现黄频在线看| 少妇被粗大猛烈的视频| 啦啦啦中文免费视频观看日本| 飞空精品影院首页| 国产精品99久久久久久久久| 只有这里有精品99| 中文字幕最新亚洲高清| 午夜激情av网站| 亚洲丝袜综合中文字幕| 国产免费一级a男人的天堂| 777米奇影视久久| 中文乱码字字幕精品一区二区三区| 涩涩av久久男人的天堂| 国产高清不卡午夜福利| 又大又黄又爽视频免费| 国产精品女同一区二区软件| 日韩一区二区三区影片| 欧美97在线视频| 交换朋友夫妻互换小说| 欧美精品国产亚洲| 人妻人人澡人人爽人人| videos熟女内射| 九九在线视频观看精品| 69精品国产乱码久久久| 一区二区三区四区激情视频| 18禁观看日本| 女性生殖器流出的白浆| 日韩不卡一区二区三区视频在线| 久久国产亚洲av麻豆专区| 精品99又大又爽又粗少妇毛片| 精品少妇黑人巨大在线播放| 亚洲美女黄色视频免费看| 视频中文字幕在线观看| 另类精品久久| 日韩熟女老妇一区二区性免费视频| 如日韩欧美国产精品一区二区三区 | 成人国产av品久久久| 汤姆久久久久久久影院中文字幕| 男女边吃奶边做爰视频| a 毛片基地| 久久久欧美国产精品| 欧美三级亚洲精品| 各种免费的搞黄视频| 99久国产av精品国产电影| 欧美精品高潮呻吟av久久| 狠狠婷婷综合久久久久久88av| 我的女老师完整版在线观看| 王馨瑶露胸无遮挡在线观看| 伊人久久国产一区二区| 熟女电影av网| 少妇的逼好多水| 亚洲国产精品成人久久小说| 在线观看一区二区三区激情| 亚洲欧美一区二区三区黑人 | 男女无遮挡免费网站观看| 夫妻性生交免费视频一级片| 日本wwww免费看| 丰满少妇做爰视频| 在线看a的网站| 久久久久视频综合| 美女国产高潮福利片在线看| 99久久精品一区二区三区| 母亲3免费完整高清在线观看 | 99热这里只有精品一区| 少妇高潮的动态图| 97超视频在线观看视频| 五月天丁香电影| 成年女人在线观看亚洲视频| 男女啪啪激烈高潮av片| √禁漫天堂资源中文www| 欧美精品亚洲一区二区| 欧美一级a爱片免费观看看| 永久网站在线| 国产一区亚洲一区在线观看| 亚洲国产最新在线播放| 亚洲三级黄色毛片| 久久免费观看电影| 视频中文字幕在线观看| 欧美 日韩 精品 国产| 在线播放无遮挡| 成人毛片60女人毛片免费| av有码第一页| 亚洲人成网站在线观看播放| 久久久久久久大尺度免费视频| a 毛片基地| 五月天丁香电影| 免费观看无遮挡的男女| 777米奇影视久久| 亚洲av男天堂| 91精品伊人久久大香线蕉| 久久精品国产a三级三级三级| 国产av码专区亚洲av| 18禁裸乳无遮挡动漫免费视频| 午夜激情av网站| 亚洲久久久国产精品| 美女国产视频在线观看| 一本大道久久a久久精品| 午夜免费观看性视频| 大话2 男鬼变身卡| 欧美另类一区| 国产精品国产三级国产专区5o| 伦精品一区二区三区| 国产欧美日韩综合在线一区二区| 一本大道久久a久久精品| 久久久精品免费免费高清| 街头女战士在线观看网站| 久久97久久精品| 涩涩av久久男人的天堂| 午夜福利网站1000一区二区三区| 男女国产视频网站| 嘟嘟电影网在线观看| 夫妻性生交免费视频一级片| 五月玫瑰六月丁香| 水蜜桃什么品种好| 一个人看视频在线观看www免费| 26uuu在线亚洲综合色| 一本大道久久a久久精品| 午夜免费男女啪啪视频观看| 哪个播放器可以免费观看大片| 亚洲成人手机| 日本91视频免费播放| av网站免费在线观看视频| 综合色丁香网| 99久久人妻综合| 三级国产精品欧美在线观看| 欧美性感艳星| 中文字幕免费在线视频6| 国产毛片在线视频| 搡女人真爽免费视频火全软件| 老司机亚洲免费影院| 国产精品国产av在线观看| 69精品国产乱码久久久| 日韩免费高清中文字幕av| 日本欧美视频一区| 22中文网久久字幕| 亚洲,一卡二卡三卡| tube8黄色片| 18在线观看网站| 免费观看无遮挡的男女| 国产精品麻豆人妻色哟哟久久| 婷婷色麻豆天堂久久| 91久久精品国产一区二区成人| av在线老鸭窝| 大香蕉久久成人网| 性色avwww在线观看| 在线观看免费视频网站a站| 美女主播在线视频| 国产免费一区二区三区四区乱码| 亚洲第一av免费看| 欧美精品一区二区免费开放| 国产成人a∨麻豆精品| 男人添女人高潮全过程视频| 欧美精品一区二区免费开放| 亚洲成人手机| av视频免费观看在线观看| 日韩伦理黄色片| 老司机亚洲免费影院| 91精品一卡2卡3卡4卡| 国产成人免费观看mmmm| 秋霞在线观看毛片| 亚洲成色77777| 久久久精品94久久精品| 三上悠亚av全集在线观看| 欧美三级亚洲精品| 制服人妻中文乱码| 亚洲人成网站在线播| 亚洲第一区二区三区不卡| 国产亚洲最大av| 自线自在国产av| 汤姆久久久久久久影院中文字幕| 91久久精品国产一区二区成人| 久久久国产欧美日韩av| 久久久午夜欧美精品| 亚洲婷婷狠狠爱综合网| 秋霞伦理黄片| 午夜av观看不卡| 99久国产av精品国产电影| 国产毛片在线视频| 97超视频在线观看视频| 国产精品国产三级国产专区5o| 国产高清不卡午夜福利| 男女高潮啪啪啪动态图| 丝瓜视频免费看黄片| 国产永久视频网站| 午夜福利视频精品| 涩涩av久久男人的天堂| 汤姆久久久久久久影院中文字幕| 一二三四中文在线观看免费高清| 久久狼人影院| 又大又黄又爽视频免费| 精品人妻熟女av久视频| 国产精品一区二区三区四区免费观看| 日韩av免费高清视频| 久久久久久人妻| 久久99一区二区三区| 美女xxoo啪啪120秒动态图| 欧美一级a爱片免费观看看| 少妇丰满av| 在线 av 中文字幕| 精品人妻一区二区三区麻豆| 久久久亚洲精品成人影院| 少妇猛男粗大的猛烈进出视频| 精品一区二区三卡| av福利片在线| 亚洲欧美一区二区三区国产| .国产精品久久| 在线观看国产h片| 极品少妇高潮喷水抽搐| 亚洲精品亚洲一区二区| 亚洲成人手机| 中国美白少妇内射xxxbb| 99热网站在线观看| 亚洲精品第二区| 高清午夜精品一区二区三区| 国产精品女同一区二区软件| 久久久久久久亚洲中文字幕| 亚洲精品,欧美精品| 午夜老司机福利剧场| 高清视频免费观看一区二区| 大片电影免费在线观看免费| 国产欧美日韩一区二区三区在线 | 在线观看免费高清a一片| 日本欧美视频一区| 日韩中字成人| a级毛片黄视频| 欧美激情极品国产一区二区三区 | 一级毛片 在线播放| 嘟嘟电影网在线观看| 男女边摸边吃奶| 最近手机中文字幕大全| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久99热这里只频精品6学生| 国产精品麻豆人妻色哟哟久久| 人人妻人人澡人人看| a级毛片黄视频| 人体艺术视频欧美日本| 精品少妇久久久久久888优播| 免费黄频网站在线观看国产| 成人亚洲精品一区在线观看| 日本av免费视频播放| 男人爽女人下面视频在线观看| 伦精品一区二区三区| 夫妻午夜视频| 国产成人精品一,二区| 哪个播放器可以免费观看大片| 亚洲欧美日韩另类电影网站| 我的老师免费观看完整版| 欧美日韩国产mv在线观看视频| 人妻一区二区av| 亚洲内射少妇av| 欧美一级a爱片免费观看看| 久久久久久久久大av| 嘟嘟电影网在线观看| 国产精品国产三级专区第一集| 成年人午夜在线观看视频| 99国产综合亚洲精品| 91久久精品国产一区二区成人| 另类亚洲欧美激情| 午夜激情久久久久久久| 新久久久久国产一级毛片| 婷婷色综合大香蕉| 99热这里只有是精品在线观看| 欧美精品人与动牲交sv欧美| 好男人视频免费观看在线| 黄色配什么色好看| 久久ye,这里只有精品| 男人添女人高潮全过程视频| 国产精品国产三级国产av玫瑰| 特大巨黑吊av在线直播| 91国产中文字幕| 久久久久精品性色| 成人亚洲欧美一区二区av| 亚洲欧美日韩卡通动漫| 伊人久久精品亚洲午夜| 青春草亚洲视频在线观看| 99久久精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 黑人巨大精品欧美一区二区蜜桃 | 欧美日韩精品成人综合77777| 欧美3d第一页| 日韩av不卡免费在线播放| 国产精品久久久久久精品电影小说| 全区人妻精品视频| videossex国产| 女性生殖器流出的白浆| 狂野欧美激情性xxxx在线观看| 亚洲丝袜综合中文字幕| av又黄又爽大尺度在线免费看| 国产色爽女视频免费观看| 亚洲伊人久久精品综合| 晚上一个人看的免费电影| 欧美日韩成人在线一区二区| 亚洲精品第二区| 人妻夜夜爽99麻豆av| 大香蕉97超碰在线| 肉色欧美久久久久久久蜜桃| 性色avwww在线观看| 色网站视频免费| 亚洲成人手机| 亚洲成色77777| 亚洲欧美日韩卡通动漫| 国产成人freesex在线| 国产综合精华液| 久久av网站| 一级二级三级毛片免费看| 国产成人午夜福利电影在线观看| 午夜久久久在线观看| 久久久精品免费免费高清| 日本猛色少妇xxxxx猛交久久| 日韩电影二区| 亚洲国产欧美在线一区| 日本欧美视频一区| xxx大片免费视频| 国产成人精品一,二区| 亚洲av男天堂| 国产高清三级在线| 精品国产一区二区久久| 国产一级毛片在线| av免费观看日本| 国产一级毛片在线| 日本色播在线视频| 美女视频免费永久观看网站| 日本午夜av视频| 晚上一个人看的免费电影| 国产 一区精品| 中文字幕制服av| 日韩精品有码人妻一区| 丝瓜视频免费看黄片| 精品国产一区二区久久| 老熟女久久久| 国产极品粉嫩免费观看在线 | 十八禁网站网址无遮挡| 国产成人av激情在线播放 | 成人手机av| 久久久久久久大尺度免费视频| 国产极品粉嫩免费观看在线 | 久热久热在线精品观看| 日本色播在线视频| 欧美 日韩 精品 国产| 街头女战士在线观看网站| 国产成人freesex在线| 中国国产av一级| 久久久久久久久久成人| 免费观看性生交大片5| 看免费成人av毛片| 亚洲成人手机| 少妇高潮的动态图| 18禁在线无遮挡免费观看视频| 国产av码专区亚洲av| 成人无遮挡网站| 美女脱内裤让男人舔精品视频| 人妻人人澡人人爽人人| 欧美三级亚洲精品| av国产久精品久网站免费入址| .国产精品久久| 国模一区二区三区四区视频| 成人亚洲欧美一区二区av| 欧美日韩国产mv在线观看视频| 老司机影院毛片| 亚洲欧美清纯卡通| 我的老师免费观看完整版| 国产黄色视频一区二区在线观看| 天天操日日干夜夜撸| 国产精品久久久久久精品电影小说| 中文天堂在线官网| 少妇 在线观看| 91国产中文字幕| 中文字幕免费在线视频6| 亚洲国产av影院在线观看| 精品亚洲成国产av| 国产精品国产三级国产av玫瑰| 午夜影院在线不卡| 欧美国产精品一级二级三级| 欧美97在线视频| 国产亚洲午夜精品一区二区久久| 国产av一区二区精品久久| 国产精品欧美亚洲77777| 少妇的逼水好多| 久久人人爽av亚洲精品天堂| 日韩成人伦理影院| 午夜免费鲁丝| 不卡视频在线观看欧美| 亚洲成人手机| 免费观看av网站的网址| 在线 av 中文字幕| 边亲边吃奶的免费视频| 99热6这里只有精品| 黄片无遮挡物在线观看| 九九爱精品视频在线观看| 最近最新中文字幕免费大全7| 精品一区二区三卡| 在线观看国产h片| 国产成人精品福利久久| 人妻人人澡人人爽人人| 永久免费av网站大全| 免费黄色在线免费观看|