• <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ì)水泥漿體電阻率與自收縮的影響
    日韩国内少妇激情av| 肉色欧美久久久久久久蜜桃| 一区在线观看完整版| 亚洲一区二区三区欧美精品| 国产爱豆传媒在线观看| 黄片无遮挡物在线观看| 亚洲成人中文字幕在线播放| 欧美精品人与动牲交sv欧美| 亚洲精华国产精华液的使用体验| 老师上课跳d突然被开到最大视频| 免费看av在线观看网站| 啦啦啦视频在线资源免费观看| 爱豆传媒免费全集在线观看| 亚洲一级一片aⅴ在线观看| 亚洲欧美成人综合另类久久久| 精品国产三级普通话版| 日日摸夜夜添夜夜添av毛片| 国产精品蜜桃在线观看| 国产久久久一区二区三区| 日本欧美视频一区| 尤物成人国产欧美一区二区三区| 蜜桃亚洲精品一区二区三区| 久久久久性生活片| 韩国av在线不卡| 特大巨黑吊av在线直播| 能在线免费看毛片的网站| 十分钟在线观看高清视频www | 久久久久国产精品人妻一区二区| 国产片特级美女逼逼视频| 久久久久网色| 日本av手机在线免费观看| av国产久精品久网站免费入址| 午夜福利影视在线免费观看| 国产精品99久久99久久久不卡 | 成人漫画全彩无遮挡| 国产av国产精品国产| 久久人人爽av亚洲精品天堂 | 国产精品久久久久久久电影| 亚洲欧美日韩无卡精品| 日韩中文字幕视频在线看片 | 亚洲国产精品999| 久久av网站| 又粗又硬又长又爽又黄的视频| 中文在线观看免费www的网站| 日日啪夜夜爽| 久久久久久久久久久丰满| av视频免费观看在线观看| 亚洲av中文字字幕乱码综合| 久久这里有精品视频免费| 日韩欧美一区视频在线观看 | 久久久久国产精品人妻一区二区| 在线亚洲精品国产二区图片欧美 | 六月丁香七月| 在线 av 中文字幕| 男人狂女人下面高潮的视频| 精品人妻一区二区三区麻豆| 观看免费一级毛片| 欧美老熟妇乱子伦牲交| 黄色配什么色好看| 永久免费av网站大全| 亚洲精品视频女| 建设人人有责人人尽责人人享有的 | a 毛片基地| 老熟女久久久| 观看免费一级毛片| 黄色日韩在线| 一级毛片黄色毛片免费观看视频| 亚洲国产高清在线一区二区三| 国产精品久久久久久久电影| 1000部很黄的大片| 国产永久视频网站| 亚洲欧美精品自产自拍| 亚洲精品,欧美精品| 99热国产这里只有精品6| 视频中文字幕在线观看| 日韩欧美精品免费久久| 欧美成人精品欧美一级黄| 国产伦理片在线播放av一区| 91午夜精品亚洲一区二区三区| 国产永久视频网站| 我的女老师完整版在线观看| av卡一久久| 日韩一区二区三区影片| 免费黄频网站在线观看国产| 人妻制服诱惑在线中文字幕| 交换朋友夫妻互换小说| 少妇丰满av| 成人特级av手机在线观看| 麻豆成人午夜福利视频| 国产一区有黄有色的免费视频| 成人免费观看视频高清| 国产午夜精品久久久久久一区二区三区| 色哟哟·www| 日韩亚洲欧美综合| 伊人久久国产一区二区| 看免费成人av毛片| 亚洲精品乱码久久久久久按摩| av又黄又爽大尺度在线免费看| 好男人视频免费观看在线| 晚上一个人看的免费电影| 人人妻人人爽人人添夜夜欢视频 | 日韩欧美一区视频在线观看 | 午夜福利在线在线| 各种免费的搞黄视频| 少妇猛男粗大的猛烈进出视频| 晚上一个人看的免费电影| 亚洲欧美日韩东京热| 亚洲成人中文字幕在线播放| 亚洲内射少妇av| 亚洲,一卡二卡三卡| 国产精品国产三级国产专区5o| 久久久a久久爽久久v久久| 黄色日韩在线| 99热国产这里只有精品6| 日本av手机在线免费观看| 日本黄色片子视频| 秋霞伦理黄片| 最近中文字幕2019免费版| 亚洲欧美日韩无卡精品| 噜噜噜噜噜久久久久久91| 简卡轻食公司| 精品一区二区三卡| 最近中文字幕2019免费版| 婷婷色综合大香蕉| 观看美女的网站| 亚洲aⅴ乱码一区二区在线播放| 久久99热这里只频精品6学生| 久久久色成人| 观看免费一级毛片| 国产伦在线观看视频一区| 日本与韩国留学比较| 深爱激情五月婷婷| 麻豆精品久久久久久蜜桃| 少妇的逼水好多| 亚洲av国产av综合av卡| 九色成人免费人妻av| 毛片女人毛片| 久久精品国产a三级三级三级| 五月开心婷婷网| 看非洲黑人一级黄片| 免费观看性生交大片5| 国产中年淑女户外野战色| 久久久a久久爽久久v久久| 国产又色又爽无遮挡免| 99国产精品免费福利视频| 全区人妻精品视频| 汤姆久久久久久久影院中文字幕| 日韩欧美 国产精品| 国产成人a∨麻豆精品| 国产亚洲91精品色在线| 久久综合国产亚洲精品| 国产精品国产三级国产专区5o| 又黄又爽又刺激的免费视频.| 精品国产一区二区三区久久久樱花 | 少妇熟女欧美另类| 中文字幕亚洲精品专区| 亚洲在久久综合| 男女啪啪激烈高潮av片| 色吧在线观看| 有码 亚洲区| 国产成人免费无遮挡视频| 丝袜脚勾引网站| 国产在线一区二区三区精| 免费看光身美女| 国内少妇人妻偷人精品xxx网站| 高清日韩中文字幕在线| 久久97久久精品| 麻豆乱淫一区二区| 国产淫语在线视频| 欧美三级亚洲精品| 国产伦精品一区二区三区视频9| 成人漫画全彩无遮挡| 简卡轻食公司| 99热这里只有是精品50| 免费av不卡在线播放| 亚洲精品,欧美精品| 亚洲内射少妇av| 国产精品久久久久久精品古装| 亚洲av二区三区四区| 日本猛色少妇xxxxx猛交久久| 婷婷色av中文字幕| 欧美精品一区二区大全| 国产免费视频播放在线视频| 中文字幕久久专区| freevideosex欧美| 亚洲精品国产色婷婷电影| 久久精品久久久久久久性| 欧美最新免费一区二区三区| 免费久久久久久久精品成人欧美视频 | 在线观看免费视频网站a站| 国产 一区精品| 我要看黄色一级片免费的| 国产精品不卡视频一区二区| 国产成人精品婷婷| 亚洲av不卡在线观看| 成年美女黄网站色视频大全免费 | 国产亚洲午夜精品一区二区久久| 中文字幕亚洲精品专区| 国产极品天堂在线| 亚洲精品国产av成人精品| 丝袜喷水一区| 波野结衣二区三区在线| 男女啪啪激烈高潮av片| 26uuu在线亚洲综合色| 97超视频在线观看视频| 观看美女的网站| 51国产日韩欧美| 又爽又黄a免费视频| 久久久久久久久久人人人人人人| 亚洲欧美日韩无卡精品| 激情五月婷婷亚洲| 国产视频内射| 日韩一区二区视频免费看| 国产男女超爽视频在线观看| 男女边摸边吃奶| 看非洲黑人一级黄片| 美女内射精品一级片tv| 久久久久久久亚洲中文字幕| 精品熟女少妇av免费看| 精品一区在线观看国产| 欧美精品亚洲一区二区| 青春草国产在线视频| 日韩亚洲欧美综合| 热99国产精品久久久久久7| 成人无遮挡网站| 超碰97精品在线观看| 日日摸夜夜添夜夜添av毛片| 3wmmmm亚洲av在线观看| 天堂8中文在线网| 女的被弄到高潮叫床怎么办| 小蜜桃在线观看免费完整版高清| 亚洲成人一二三区av| 乱系列少妇在线播放| 久久久久久久亚洲中文字幕| 亚洲欧美日韩东京热| 黄片无遮挡物在线观看| 男女边吃奶边做爰视频| 人人妻人人添人人爽欧美一区卜 | 国产精品国产三级国产专区5o| 国产免费又黄又爽又色| 欧美最新免费一区二区三区| 亚洲欧美精品专区久久| 日韩av免费高清视频| 亚洲综合精品二区| 久久久久久人妻| 91精品一卡2卡3卡4卡| av网站免费在线观看视频| 草草在线视频免费看| 国产精品偷伦视频观看了| 久久久久精品久久久久真实原创| 看非洲黑人一级黄片| 国产亚洲最大av| 日韩欧美精品免费久久| 亚洲精品456在线播放app| 久久毛片免费看一区二区三区| 欧美区成人在线视频| 久久久久久伊人网av| 高清不卡的av网站| 国产精品蜜桃在线观看| 欧美+日韩+精品| 日本与韩国留学比较| 色视频在线一区二区三区| 美女高潮的动态| 天天躁夜夜躁狠狠久久av| 国产精品免费大片| 精品人妻熟女av久视频| 蜜桃在线观看..| 国内揄拍国产精品人妻在线| 亚洲成人一二三区av| 男的添女的下面高潮视频| videossex国产| 欧美人与善性xxx| 男女边吃奶边做爰视频| 亚洲人成网站高清观看| 亚洲国产av新网站| 91精品国产九色| av在线蜜桃| 国产成人a∨麻豆精品| 人妻一区二区av| 在线观看三级黄色| 国产亚洲av片在线观看秒播厂| 日本午夜av视频| 秋霞在线观看毛片| 久久精品人妻少妇| 简卡轻食公司| 亚洲av日韩在线播放| 国产精品久久久久久av不卡| 国产一区二区三区综合在线观看 | 久久婷婷青草| 99热网站在线观看| 亚洲精品乱码久久久v下载方式| 久久99蜜桃精品久久| a 毛片基地| 夜夜看夜夜爽夜夜摸| 校园人妻丝袜中文字幕| 内射极品少妇av片p| 91久久精品国产一区二区三区| 国产成人a区在线观看| 草草在线视频免费看| 蜜桃亚洲精品一区二区三区| 亚洲欧美精品专区久久| 久久精品久久久久久久性| 在线观看av片永久免费下载| 精品人妻熟女av久视频| 国产免费福利视频在线观看| 亚洲精品乱久久久久久| 国产黄色视频一区二区在线观看| xxx大片免费视频| 国内揄拍国产精品人妻在线| 在现免费观看毛片| 久久久久久久久久成人| 国产欧美另类精品又又久久亚洲欧美| 国产乱来视频区| 国产一区二区三区综合在线观看 | 国产毛片在线视频| a级一级毛片免费在线观看| 亚洲熟女精品中文字幕| 大香蕉97超碰在线| 91狼人影院| 99久久中文字幕三级久久日本| 在线免费观看不下载黄p国产| 亚洲色图av天堂| 国产亚洲午夜精品一区二区久久| 久久久久久久国产电影| 22中文网久久字幕| 国产高清有码在线观看视频| 乱码一卡2卡4卡精品| 精华霜和精华液先用哪个| 51国产日韩欧美| 国产男女超爽视频在线观看| 观看美女的网站| 国产免费又黄又爽又色| 欧美bdsm另类| 少妇熟女欧美另类| 亚洲成人av在线免费| 啦啦啦在线观看免费高清www| 韩国高清视频一区二区三区| av福利片在线观看| 多毛熟女@视频| 亚洲精品456在线播放app| 中文天堂在线官网| 国产精品久久久久久av不卡| 成人午夜精彩视频在线观看| 少妇被粗大猛烈的视频| 五月天丁香电影| 中文资源天堂在线| 搡女人真爽免费视频火全软件| 丰满迷人的少妇在线观看| 亚洲国产色片| 久久亚洲国产成人精品v| 久久久久久久久大av| 在线免费观看不下载黄p国产| 色网站视频免费| 久久久久视频综合| 狂野欧美激情性bbbbbb| 亚洲中文av在线| 深爱激情五月婷婷| 日韩电影二区| 五月伊人婷婷丁香| 男男h啪啪无遮挡| 一级毛片aaaaaa免费看小| 亚洲美女视频黄频| 国产黄色视频一区二区在线观看| 成人高潮视频无遮挡免费网站| 菩萨蛮人人尽说江南好唐韦庄| 久久久久网色| 免费大片黄手机在线观看| 日本猛色少妇xxxxx猛交久久| 内地一区二区视频在线| 国产 一区精品| 国产淫片久久久久久久久| av福利片在线观看| 国产av国产精品国产| av视频免费观看在线观看| 久久久久视频综合| av黄色大香蕉| 国产黄色免费在线视频| 亚洲四区av| h日本视频在线播放| 国产亚洲91精品色在线| 熟妇人妻不卡中文字幕| 亚洲美女黄色视频免费看| 一本一本综合久久| 日韩在线高清观看一区二区三区| 国内揄拍国产精品人妻在线| 99热全是精品| 欧美精品人与动牲交sv欧美| 18禁裸乳无遮挡免费网站照片| 人妻制服诱惑在线中文字幕| 综合色丁香网| 97在线人人人人妻| 视频区图区小说| 啦啦啦啦在线视频资源| 在线看a的网站| 亚洲精品国产av成人精品| 熟妇人妻不卡中文字幕| 夜夜看夜夜爽夜夜摸| 欧美日韩国产mv在线观看视频 | av网站免费在线观看视频| 一本久久精品| 最近中文字幕2019免费版| 国产高清有码在线观看视频| 91精品一卡2卡3卡4卡| 丰满乱子伦码专区| 免费少妇av软件| av免费在线看不卡| 蜜桃亚洲精品一区二区三区| 亚洲激情五月婷婷啪啪| 男人和女人高潮做爰伦理| 中国美白少妇内射xxxbb| 国产成人免费观看mmmm| 国产一级毛片在线| 亚洲第一区二区三区不卡| 日韩欧美精品免费久久| 大香蕉久久网| 天美传媒精品一区二区| 91狼人影院| 精品久久国产蜜桃| 亚洲欧美一区二区三区黑人 | 国产白丝娇喘喷水9色精品| 男女边吃奶边做爰视频| av在线app专区| 蜜臀久久99精品久久宅男| 午夜福利影视在线免费观看| 尾随美女入室| 国产精品免费大片| 国产 一区 欧美 日韩| 五月伊人婷婷丁香| 久久久久性生活片| 中文天堂在线官网| 欧美xxxx黑人xx丫x性爽| 日韩,欧美,国产一区二区三区| 欧美另类一区| 亚洲国产毛片av蜜桃av| 看非洲黑人一级黄片| 性色av一级| 91午夜精品亚洲一区二区三区| 亚洲精品日本国产第一区| 你懂的网址亚洲精品在线观看| 美女福利国产在线 | 少妇人妻 视频| 国产69精品久久久久777片| 中文字幕av成人在线电影| 国产精品熟女久久久久浪| 国产爱豆传媒在线观看| 国产精品伦人一区二区| 国产爱豆传媒在线观看| 王馨瑶露胸无遮挡在线观看| 一级毛片黄色毛片免费观看视频| 最近的中文字幕免费完整| freevideosex欧美| 国产91av在线免费观看| 亚洲自偷自拍三级| 尾随美女入室| 欧美xxxx黑人xx丫x性爽| 高清黄色对白视频在线免费看 | 欧美日本视频| 日日摸夜夜添夜夜爱| 亚洲精品一区蜜桃| 亚洲av.av天堂| 在线观看一区二区三区激情| 中国国产av一级| 成人免费观看视频高清| 国产片特级美女逼逼视频| 自拍偷自拍亚洲精品老妇| 精品99又大又爽又粗少妇毛片| 各种免费的搞黄视频| 午夜免费男女啪啪视频观看| 在线观看免费高清a一片| 永久网站在线| 国产69精品久久久久777片| 国产伦精品一区二区三区视频9| 国产色婷婷99| 超碰av人人做人人爽久久| 免费av不卡在线播放| 老司机影院成人| 亚洲精品,欧美精品| 亚洲不卡免费看| 久久人人爽人人爽人人片va| 一个人看视频在线观看www免费| 欧美xxxx黑人xx丫x性爽| 天堂中文最新版在线下载| 国产精品久久久久成人av| 99视频精品全部免费 在线| 精品久久久精品久久久| 老司机影院毛片| 毛片一级片免费看久久久久| 十八禁网站网址无遮挡 | 性色avwww在线观看| 欧美激情极品国产一区二区三区 | 秋霞在线观看毛片| 久久精品人妻少妇| 国产精品99久久99久久久不卡 | 大香蕉97超碰在线| 日本-黄色视频高清免费观看| 国语对白做爰xxxⅹ性视频网站| 国产av精品麻豆| 性色avwww在线观看| 你懂的网址亚洲精品在线观看| 亚洲欧美清纯卡通| 日本免费在线观看一区| 亚洲精品乱码久久久v下载方式| 亚洲精品日韩av片在线观看| 日韩一区二区三区影片| 在线精品无人区一区二区三 | 免费黄网站久久成人精品| 日本wwww免费看| 亚洲人成网站在线播| 亚洲欧美成人综合另类久久久| 尤物成人国产欧美一区二区三区| 亚洲国产欧美人成| 看十八女毛片水多多多| 亚洲欧美中文字幕日韩二区| 国产精品偷伦视频观看了| 国产久久久一区二区三区| 精品久久国产蜜桃| 女性生殖器流出的白浆| 欧美精品亚洲一区二区| 亚洲欧美成人综合另类久久久| 亚洲丝袜综合中文字幕| 免费观看性生交大片5| 欧美日韩一区二区视频在线观看视频在线| 韩国av在线不卡| www.色视频.com| 国产又色又爽无遮挡免| 亚洲国产最新在线播放| 2022亚洲国产成人精品| av黄色大香蕉| 美女xxoo啪啪120秒动态图| 网址你懂的国产日韩在线| 国产精品久久久久久精品电影小说 | 欧美少妇被猛烈插入视频| 亚洲色图av天堂| 亚洲欧美一区二区三区国产| 日本wwww免费看| 另类亚洲欧美激情| 亚洲高清免费不卡视频| 亚洲激情五月婷婷啪啪| 高清在线视频一区二区三区| 国产精品精品国产色婷婷| 国产av码专区亚洲av| 亚洲四区av| 免费观看a级毛片全部| 日韩视频在线欧美| 亚洲色图综合在线观看| 国产av码专区亚洲av| 亚洲国产精品成人久久小说| 中文字幕av成人在线电影| 久久久色成人| 在线观看人妻少妇| 国产欧美另类精品又又久久亚洲欧美| 成人美女网站在线观看视频| av一本久久久久| 不卡视频在线观看欧美| 国产免费一区二区三区四区乱码| 欧美日韩综合久久久久久| 又大又黄又爽视频免费| 观看免费一级毛片| 免费久久久久久久精品成人欧美视频 | 免费观看av网站的网址| 尤物成人国产欧美一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 网址你懂的国产日韩在线| 2021少妇久久久久久久久久久| 亚洲激情五月婷婷啪啪| 91精品国产九色| 视频区图区小说| 欧美一级a爱片免费观看看| 成人影院久久| 亚洲色图av天堂| 一本一本综合久久| av在线观看视频网站免费| 成年av动漫网址| 草草在线视频免费看| 欧美97在线视频| 黑人高潮一二区| 爱豆传媒免费全集在线观看| 日韩视频在线欧美| 成人18禁高潮啪啪吃奶动态图 | 日本av免费视频播放| 中文字幕亚洲精品专区| 99久久精品热视频| 人妻 亚洲 视频| 一级毛片久久久久久久久女| 亚洲av二区三区四区| 亚洲国产精品成人久久小说| 日本av免费视频播放| 熟妇人妻不卡中文字幕| 又粗又硬又长又爽又黄的视频| 国产成人一区二区在线| 一区二区三区乱码不卡18| av天堂中文字幕网| 久久久精品94久久精品| 国产免费又黄又爽又色| 舔av片在线| 99re6热这里在线精品视频| 大片免费播放器 马上看| 肉色欧美久久久久久久蜜桃| 18禁裸乳无遮挡免费网站照片| 91狼人影院| 人妻少妇偷人精品九色| 亚洲精品日韩av片在线观看| 久久久久久久久久久丰满| 九色成人免费人妻av| 高清视频免费观看一区二区| 综合色丁香网| 午夜福利在线观看免费完整高清在| 少妇 在线观看| 欧美bdsm另类| 国产大屁股一区二区在线视频| 噜噜噜噜噜久久久久久91| .国产精品久久| 少妇熟女欧美另类| 色婷婷av一区二区三区视频| 人人妻人人看人人澡|