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

    基于高密度電阻率法的近地表工程巖體速度結(jié)構(gòu)成像

    2018-01-03 01:04:55孟慶生劉圣彪肖志廣
    物探化探計算技術(shù) 2017年6期
    關(guān)鍵詞:縱波波速高密度

    趙 釗, 孟慶生, 韓 凱, 劉圣彪, 肖志廣

    (山東省海洋環(huán)境地質(zhì)工程重點實驗室(中國海洋大學(xué)),青島 266100)

    基于高密度電阻率法的近地表工程巖體速度結(jié)構(gòu)成像

    趙 釗, 孟慶生*, 韓 凱, 劉圣彪, 肖志廣

    (山東省海洋環(huán)境地質(zhì)工程重點實驗室(中國海洋大學(xué)),青島 266100)

    縱波速度劃分巖體等級對地下空間工程具有重要的意義?;贔aust公式,建立了表征電阻率和縱波速度的超定方程組,并研究了β系數(shù)的確定方法。通過數(shù)值模擬手段,驗證了近地表巖體縱波速度成像方法的有效性。結(jié)果表明,該方法對層狀地層、含斷層、高低阻地質(zhì)異常等復(fù)雜地質(zhì)結(jié)構(gòu)均具有良好的反映能力。應(yīng)用實例表明,計算值與實測值之間的誤差小于8.7%,可為工程勘察提供可信的縱波速度。

    高密度電阻率; Faust公式; 速度成像

    0 引言

    隨著我國城市化的高速發(fā)展,人均可利用的地面資源日益匱乏。為解決這一矛盾,修建城市地鐵,地下制作車間、電站、人防等地下空間工程已成為發(fā)達國家城市建設(shè)的必然選擇[1]。但是由于地下地質(zhì)條件的復(fù)雜性,給地下空間工程建設(shè)帶來了嚴(yán)峻挑戰(zhàn)[2]。為合理建設(shè)地下工程,需要在勘察階段詳細了解巖體的工程特性。其中巖體等級的準(zhǔn)確劃分對選擇合理的設(shè)計方案、施工方法,以及確定工程結(jié)構(gòu)荷載等具有重要的意義。

    在巖土工程領(lǐng)域,巖體的等級主要以圍巖中的縱波速度作為劃分標(biāo)準(zhǔn)[3-4]。傳統(tǒng)獲取巖石縱波速度的方法是通過鉆孔取芯,然后對巖芯進行室內(nèi)測試獲取。雖然該方法可以直接獲取巖體的縱波速度,但受場地條件和投入的限制不能大規(guī)模開展。近年來,一些以物探方法為主的巖體縱波速度測試手段得到了廣泛應(yīng)用,常用的有折射波法[5]、瑞利面波法[6],這類方法具有無損、快速和經(jīng)濟的優(yōu)點[7]。然而,在城市環(huán)境中開展上述以彈性振動為基礎(chǔ)的方法時,由于外界震動干擾源多,常常無法獲得高分辨率的地震數(shù)據(jù),導(dǎo)致無法獲得準(zhǔn)確的縱波速度。

    高密度電阻率法是以地下介質(zhì)的電性差異為基礎(chǔ)的一類方法,可依據(jù)電阻率的變化定性確定巖體的分界面,在地質(zhì)勘察和地質(zhì)災(zāi)害探測中發(fā)揮了重要的作用[8]。此外,通過對數(shù)百萬英尺測井資料的分析證實,電阻率與地下介質(zhì)的縱波速度具有良好的關(guān)系[9]。將兩種手段結(jié)合,可獲得工程巖體的縱波速度,進而為巖體等級劃分提供依據(jù)。筆者基于Faust公式,結(jié)合高密度電阻率法,并利用少量鉆孔標(biāo)定的方法,探討一種工程巖體速度建模新方法,為地下空間工程中巖體等級劃分提供有利的工具。

    1 近地表速度建模理論

    1.1 Faust 公式[10]

    上世紀(jì)40年代后期和50年代中期,F(xiàn)aust分析了數(shù)百萬尺的測井資料,將電阻率與速度進行對比,證實了巖石的縱波速度與埋深和巖石的電性(電阻率)有關(guān),即Faust方程,如式(1)所示。

    V=AZTR1/6

    (1)

    式中:A為常數(shù);Z為深度(m);T為巖性因子;R為電阻率(Ω·m)。

    結(jié)合阿爾奇(Archie)經(jīng)驗公式[11]和威利(Wyllie)時間平均方程[12],可將Faust公式改寫為式(2)。

    v=kραHβ

    (2)

    式中:v為縱波波速(m/s);k為與巖石性質(zhì)有關(guān)的經(jīng)驗常數(shù);H為深度(m);α、β為與地層有關(guān)的經(jīng)驗常數(shù);ρ為電阻率(Ω·m)。在沒有聲波測井資料但有電阻率數(shù)據(jù)的地區(qū),式(2)可以用來獲取巖體的縱波速度。

    1.2 基于廣義逆法的速度反演

    將Faust公式的兩邊取對數(shù)可得式(3)。

    lgv=lgk+αlgρ+βlgH

    (3)

    式(3)反映出的lgv、lgρ和lgk的關(guān)系只是近似的線性關(guān)系,不能精確反映速度與電阻率和深度之間復(fù)雜的非線性關(guān)系。為此,筆者以lgρ和lgH的高階多項式來擬合地層速度,在每一深度如式(4)所示。

    v=a1lgρHβ+a2(lgρHβ)2+a3

    (4)

    寫成矩陣形式為

    V=XA

    (5)

    式(5)中:

    V=|v1,v2,v3,…,vn|T|

    A=|a1,a2,a3|T

    對于方程中的系數(shù)β,可以通過實測場地鉆孔位置的波速、電阻率和深度資料,直接擬合分析v與lgρHβ之間的關(guān)系求得。式(4)為超定方程組,采用廣義逆法可以求出系數(shù)α1、α2、α3。在無鉆孔的位置,利用測得的電阻率剖面數(shù)據(jù)插值每一深度段的電阻率,再利用求得的系數(shù)即可求出縱波速度。

    對于求解出的各深度的速度值來說,可能會出現(xiàn)比實際波速大或小的情況。我們將作如下處理,首先對鉆孔波速測試數(shù)據(jù)在深度方向進行與電阻率數(shù)據(jù)等密度的插值,得到一組波速數(shù)據(jù)

    V1=|v1,v2,v3,v4|T

    (6)

    設(shè)利用電阻率成像算法計算得到的波速數(shù)據(jù)為v2,兩者之差為

    ▽V=V2-V1

    (7)

    令m=▽V/V1,則權(quán)系數(shù)定義為

    W=1-m

    (8)

    修正后的波速為式(9)。

    V加權(quán)=W·V

    (9)

    2 有效性驗證

    2.1 均勻地層

    設(shè)計層狀均勻地層,地電模型參數(shù)如表1所示。

    表1 均勻?qū)訝畹仉娔P?/p>

    圖1 均勻地層反演模擬結(jié)果Fig.1 The inverse modeling result of layered formation

    利用Wenner裝置對上述模型正、反演后的結(jié)果見圖1,裝置參數(shù),電極為101個,電極間距為1 m。

    利用公式(5)求得的方程系數(shù)為β=0.308 5,a1=-235.821 1,a2=52.435 0,a3=867.970 7,進而求得剖面中其他各點隨深度變化的縱波速度,其反演速度剖面見圖2。

    由圖2可見,電阻率波速也呈現(xiàn)良好的正相關(guān)關(guān)系,各層波速與模型基本一致并隨深度增加而增大。

    2.2 含高低阻帶地層

    圖2 均勻地層縱波速度成像剖面Fig.2 The velocity imaging of layered formation

    圖3 含高低阻地層地電模型Fig.3 The abnormal formation with low or high resistance model

    均勻地層是一種十分理想的地質(zhì)模型,實際的地層是非常復(fù)雜的,常常含有斷層、溶洞等不良地質(zhì)體,總體上可分為高阻體和低阻體兩類。為驗證速度建模方法在復(fù)雜地質(zhì)條件下的適用性,分別在均勻地層模型中設(shè)立高阻區(qū)和低阻區(qū)(圖3),其電阻率分別為3 000 Ω·m和300 Ω·m。

    利用上述方法,獲得相應(yīng)的電阻率剖面(圖4)和速度模型(圖5)。

    圖4 高阻、低阻模型電阻率剖面Fig.4 The resistivity model profile of abnormal formation whth high and low resistance(a)高阻模型;(b)低阻模型

    圖5 高阻、低阻反演速度成像剖面Fig.5 The inverse velocity image of high and low resistance(a)高阻;(b)低阻

    由圖5可知,在模型的高、低阻異常處,波速值也有明顯地變化,分別表現(xiàn)為高波速異常和低波速異常,與模型中設(shè)計的異常位置非常吻合。證實了該方法對復(fù)雜的地層具有很好的表現(xiàn)能力。在實際的工程應(yīng)用中,可以根據(jù)波速的變化推斷異常體的存在。

    3 應(yīng)用實例

    3.1 試驗場地概況

    試驗選取某工程場地進行。由淺到深地層依次為人工填土(0.5 m ~5.0 m),安山巖強風(fēng)化帶、中等風(fēng)化帶和微風(fēng)化帶。為獲得方程系數(shù),在場地中的一條高密度電阻率探測測線上分別布設(shè)了Z38、Z39和Z40三個鉆孔,并在孔內(nèi)測量各層的聲波速度(表2)。本次試驗中用Z38號孔作為波速反演的參照孔,Z39和Z40孔作為驗證孔。

    高密度電阻率法試驗采用吉林大學(xué)工程技術(shù)研究所生產(chǎn)的E60DN高密度電法儀,測量裝置為Wenner裝置,電極數(shù) 為64個,電極距為2.0 m,采集層數(shù)為14層。試驗數(shù)據(jù)經(jīng)反演后的剖面見圖6。

    表2 試驗場地孔內(nèi)縱波速度

    圖6 實際場地高密度電阻率法探測剖面Fig.6 The ground profile by the high density resistivity method

    圖7 深度H與v/lg ρ擬合曲線Fig.7 The line of H and v/lg ρ

    3.2 巖體速度結(jié)構(gòu)成像

    本次試驗分別獲得了0.50 m、1.50 m、2.55 m、3.70 m、4.98 m、6.37 m、7.91 m、9.60 m及11.46 m深度處的電阻率值。而波速測試的深度點分別為1.0 m、2.0 m、3.0 m、4.0 m、5.0 m、6.0 m、7.0 m、8.0 m、9.0 m、10.0 m、11.0 m,二者不一致。為此,需將電阻率數(shù)據(jù)進行插值,得到與波速數(shù)據(jù)點處相對應(yīng)的電阻率值。

    在得到Z38鉆孔處的波速數(shù)據(jù)及電阻率數(shù)據(jù)后,使用origin軟件中Allometric型冪指數(shù)擬合方式求取系數(shù),擬合曲線見圖7,其中橫坐標(biāo)為H,縱坐標(biāo)為v/lgρ,求得的β系數(shù)為0.33。

    利用Z38孔處波速、電阻率和β系數(shù),將數(shù)據(jù)應(yīng)用于超定方程組中,求得系數(shù)a1=-785.526 8,a2=-76.570 6,a3=-591.128 8。進而獲得各測點不同深度的縱波速度剖面(圖8)。

    圖8 場地速度結(jié)構(gòu)剖面圖Fig.8 The velocity structure section in the ground

    3.3 誤差分析

    為驗證基于Faust公式速度建模方法的準(zhǔn)確性,將獲得的波速與Z39、Z40鉆孔處測量的波速進行比較(表3)。通過對反演電阻率數(shù)據(jù)與Z39、Z40孔實測波速的對比發(fā)現(xiàn),反演成像后的波速數(shù)據(jù)與計算數(shù)據(jù)具有較高的一致性,最大誤差為8.7%。個別深度點波速出現(xiàn)過大或過小的原因是由于本試驗中擬合β系數(shù)的數(shù)據(jù)點較少,從而導(dǎo)致的擬合精度不高。

    4 結(jié)論

    表3 反演成像速度與Z39、Z40孔測量波速對比

    1)理論與實際結(jié)果證實,基于Faust公式和高密度電阻率法的縱波速度建模技術(shù),對均勻地層、含低阻異常體地層和含高阻異常體地層均具有良好的適用性。

    2)實際場地應(yīng)用結(jié)果表明:在該案例中,相對誤差不超過8.7%。利用電阻率數(shù)據(jù)反演出的縱波速度剖面與驗證鉆孔揭示的縱波速度具有良好的一致性。

    3)筆者對夾層、互層等電阻率突變地質(zhì)條件下未進行驗證,在今后的工作中需進一步驗證。

    [1] 趙晉友,黃松.城市地下空間開發(fā)面臨的機遇與挑戰(zhàn)[J].地質(zhì)與勘探,2013,49(5):0964-0969.

    ZHAO J Y,HUANG S.Opportunities and challenges in the development of urban underground space[J].Geology and Exploration,2013,49( 5) : 0964-0969.(In Chinese)

    [2] 萬漢斌.城市高密度地區(qū)地下空間開發(fā)策略研究[D].天津:天津大學(xué),2013.

    WAN H B. The strategy research on urban underground space in high-density area [D]. Tianjin: Tianjin University,2013. (In Chinese)

    [3] 王洪勇. 關(guān)于應(yīng)用彈性波速度劃分隧道圍巖類別的認(rèn)識[J]. 隧道建設(shè), 2000, 20(3): 19-25.

    WANG H Y. Some knowledge about the separation of tunnel surrounding rock via the elastic wave velocity [J].Tunnel Construction,2000,20(3):19-25. (In Chinese)

    [4] 許廣春. 大地電磁測深進行隧道圍巖分級的可行性研究[J].鐵道工程學(xué)報,2010,27(8):83-86.

    XU G C. Research on feasibility of sorting surrounding rock of tunnel with telluric electromagnetic sounding [J]. Journal of Railway Engineering Society, 2010,27(8):83-86. (In Chinese)

    [5] 陳愛萍,鄒文,何光明,等.初至波地震層析成像中自動生成初始速度模型的方法研究[J]. 物探化探計算技術(shù),2014,157(5):583-586.

    CHEN A P, ZOU W, HE G M, et al. Automatically building the initial velocity model in first-arrival seismic tomography [J]. Computing Techniques for Geophysical and Geochemical Exploration, 2014,157(5):583-586. (In Chinese)

    [6] 周竹生,劉喜亮,熊孝雨.彈性介質(zhì)中瑞雷面波有限差分法正演模擬[J].地球物理學(xué)報,2007,50(2):567-573.

    ZHOU Z S, LIU X L, XIONG XY. Finite-difference modelling of Rayleigh surface wave in elastic media [J]. Chinese Journal of Geophysics, 2007, 50(2):567-573.(In Chinese)

    [7] 程紅杰,胡祥云,張榮,等.近地表地球物理發(fā)展?fàn)顩r綜述[J].工程地球物理學(xué)報,2005,2(1):73-76.

    CHENG H J, HU X Y, ZHANG R, et al.The development of the near surface geophysics [J].Chinese Journal of Engineering Geophysics, 2005, 2(1):73-76. (In Chinese)

    [8] 嚴(yán)加永,孟貴祥,呂慶田,等.高密度電法的進展與展望[J].物探與化探,2012,36(4):576-584.

    YAN J Y, MENG G X, LV Q T, et al. The progress and prospect of the electrical resistivity imaging survey [J].Geophysical and Geochemical Exploration, 2012,36(4):576-584. (In Chinese)

    [9] FAUST L Y. A velocity function including lithologic variation [J]. Geophysics, 1953, 18(2): 271-288.

    [10] FAUST L Y. Seismic velocity as a function of depth and geologic time [J]. Geophysics, 1951, 16(2): 192-206.

    [11] 孫建國. 阿爾奇公式:提出背景與早期爭論[J].應(yīng)用物理學(xué),2007,22(2):472-486.

    SUN J G. Archie’s formula: historical background and earlier debates[J].Applied Physics,2007,22(2):472-486. (In Chinese)

    [12] WOOD W T, STOFFA P L, SHIPLEY T H. Quantitative detection of methane hydrate through high-resolution seismic velocity analysis [J]. J Geophys Res,1994, 99: 9681-9695.

    Rockmassvelocityimagingbasedonhighdensityresistivitymethodinnearsurfaceengineering

    ZHAO Zhao, MENG Qingsheng*, HAN Kai, LIU Shengbiao, XIAO Zhiguang

    (Shandong Provincial Key Laboratory of Marine Environment and Geological Engineering (Ocean University of China), Qingdao 266100, China)

    It is meaningful to the underground engineer that the P-wave velocity is divided by the rock grade. According to the Faust formula, the overdetermined equation of resistivity data and P-wave velocity is built, and the way of β determination is researched. Validity of this rock mass velocity map imaging way is proved according to the numerical simulation. The result shows this way is useful to the complex geological structure like layered formation, fault, and abnormal formations with high or low resistance. The application example shows the error of the calculate and the measure results is less than 8.7%, the P-wave velocity calculated by this way is credible to the geotechnical investigation and surveying.

    high density resistivity; Faust formula; velocity imaging

    2016-11-13 改回日期: 2017-02-21

    國家重大科研儀器研制項目(41427803)

    趙釗(1994-),男,碩士,主要從事環(huán)境工程地球物理探測方面的研究工作,E-mail:705586405@qq.com。

    孟慶生(1972-),男,副教授,主要從事環(huán)境工程地球物理探測方面的教學(xué)和研究工作,E-mail:sdmengs@163.com。

    1001-1749(2017)06-0742-06

    P 631.3

    A

    10.3969/j.issn.1001-1749.2017.06.05

    猜你喜歡
    縱波波速高密度
    高密度電法在斷裂構(gòu)造探測中的應(yīng)用
    基于實測波速探討地震反射波法超前預(yù)報解譯標(biāo)志
    高密度電法在尋找地下水中的應(yīng)用
    黃257井區(qū)疊前縱波方位各向異性裂縫分布預(yù)測
    城市高密度環(huán)境下的建筑學(xué)探討
    吉林地區(qū)波速比分布特征及構(gòu)造意義
    變截面階梯桿中的縱波傳播特性實驗
    高密度互連技術(shù)強勁發(fā)展
    基于分位數(shù)回歸的剪切波速變化規(guī)律
    德布羅意關(guān)系式的相對論協(xié)變形式及物質(zhì)波波速
    物理與工程(2011年2期)2011-03-25 10:02:59
    亚洲精品粉嫩美女一区| 国产精品国产三级国产av玫瑰| 国产精品久久久久久精品电影| 一a级毛片在线观看| 国产精品自产拍在线观看55亚洲| 男女视频在线观看网站免费| 亚洲国产精品久久男人天堂| 午夜免费男女啪啪视频观看 | 日本黄色视频三级网站网址| 成人毛片a级毛片在线播放| 欧美一区二区精品小视频在线| 国产精品一及| 午夜免费激情av| 久99久视频精品免费| 十八禁国产超污无遮挡网站| 99视频精品全部免费 在线| 午夜福利欧美成人| 日本 欧美在线| 欧美成人免费av一区二区三区| 在线国产一区二区在线| 精品人妻1区二区| 亚洲自拍偷在线| 综合色av麻豆| 精品无人区乱码1区二区| 久久精品国产亚洲av天美| 国产精品乱码一区二三区的特点| 2021天堂中文幕一二区在线观| 欧美日韩精品成人综合77777| a级毛片a级免费在线| 男人和女人高潮做爰伦理| 日本在线视频免费播放| 久久久久免费精品人妻一区二区| 亚洲精品一区av在线观看| 麻豆国产97在线/欧美| 国产亚洲精品久久久com| 亚洲精品乱码久久久v下载方式| 黄色配什么色好看| 精品久久久久久久久亚洲 | 在线观看66精品国产| 色噜噜av男人的天堂激情| 一个人免费在线观看电影| 午夜激情欧美在线| 亚洲av不卡在线观看| 免费看av在线观看网站| 香蕉av资源在线| 亚洲精品日韩av片在线观看| 欧美日韩精品成人综合77777| 国产精品国产三级国产av玫瑰| 久久久精品欧美日韩精品| 免费人成视频x8x8入口观看| 亚洲四区av| 国产日本99.免费观看| 欧美高清性xxxxhd video| 成人美女网站在线观看视频| 国产成人一区二区在线| 国产精品亚洲一级av第二区| 久久欧美精品欧美久久欧美| 精品一区二区免费观看| 亚洲无线在线观看| 性色avwww在线观看| 天堂影院成人在线观看| 亚洲中文字幕一区二区三区有码在线看| 99热只有精品国产| 国产精品免费一区二区三区在线| 久久精品91蜜桃| 久久亚洲精品不卡| 国产一区二区激情短视频| 嫩草影院入口| 国产成人av教育| 日本一本二区三区精品| 欧美日韩国产亚洲二区| 日韩中文字幕欧美一区二区| 国产三级中文精品| 2021天堂中文幕一二区在线观| 99久久中文字幕三级久久日本| 91久久精品国产一区二区三区| 岛国在线免费视频观看| 午夜精品在线福利| 又爽又黄无遮挡网站| 久久久久久久亚洲中文字幕| 熟妇人妻久久中文字幕3abv| 午夜福利在线观看免费完整高清在 | 国产精品,欧美在线| 制服丝袜大香蕉在线| 国产精品一区www在线观看 | 色精品久久人妻99蜜桃| 淫妇啪啪啪对白视频| 国产精品亚洲一级av第二区| 亚洲精品粉嫩美女一区| 国产单亲对白刺激| 一个人看视频在线观看www免费| 五月伊人婷婷丁香| 亚洲四区av| 99热这里只有是精品50| 国产高潮美女av| 亚洲精品色激情综合| 欧美xxxx黑人xx丫x性爽| 免费黄网站久久成人精品| 伊人久久精品亚洲午夜| 色噜噜av男人的天堂激情| 一区二区三区四区激情视频 | 亚洲熟妇熟女久久| 狂野欧美白嫩少妇大欣赏| 精品久久久噜噜| 不卡一级毛片| 国产av不卡久久| av女优亚洲男人天堂| 啦啦啦观看免费观看视频高清| 午夜福利在线观看免费完整高清在 | 一级av片app| 亚洲成av人片在线播放无| 久久久午夜欧美精品| 欧美一区二区国产精品久久精品| 国产视频内射| 久久久久久久久中文| 男人舔奶头视频| 久久久久久大精品| 三级毛片av免费| 中文字幕高清在线视频| 成人国产一区最新在线观看| 黄色视频,在线免费观看| 国产aⅴ精品一区二区三区波| 两性午夜刺激爽爽歪歪视频在线观看| 精品人妻视频免费看| 超碰av人人做人人爽久久| 亚洲一区高清亚洲精品| 变态另类成人亚洲欧美熟女| 波多野结衣高清作品| 久久午夜亚洲精品久久| 色综合婷婷激情| 婷婷精品国产亚洲av在线| 亚洲无线在线观看| 国产单亲对白刺激| av福利片在线观看| 亚洲熟妇中文字幕五十中出| 少妇人妻精品综合一区二区 | 神马国产精品三级电影在线观看| 丰满乱子伦码专区| 欧美色欧美亚洲另类二区| 亚洲国产日韩欧美精品在线观看| 精品久久久久久久久亚洲 | 熟妇人妻久久中文字幕3abv| 国产免费av片在线观看野外av| 成人无遮挡网站| 大又大粗又爽又黄少妇毛片口| 日韩欧美国产一区二区入口| 久久午夜福利片| 国产aⅴ精品一区二区三区波| 国产高清视频在线播放一区| 日韩av在线大香蕉| 免费人成在线观看视频色| 亚洲精品亚洲一区二区| 一边摸一边抽搐一进一小说| 色5月婷婷丁香| 国产一区二区在线观看日韩| 熟女人妻精品中文字幕| 黄色日韩在线| 舔av片在线| 女的被弄到高潮叫床怎么办 | 免费在线观看成人毛片| 欧美激情久久久久久爽电影| 国产av不卡久久| 国产亚洲欧美98| 国产精品一区二区三区四区免费观看 | 国产一区二区在线av高清观看| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区三区四区免费观看 | 男女下面进入的视频免费午夜| 蜜桃亚洲精品一区二区三区| 欧美成人性av电影在线观看| 毛片女人毛片| 啦啦啦韩国在线观看视频| 国产高清不卡午夜福利| 一a级毛片在线观看| 欧美日韩乱码在线| 非洲黑人性xxxx精品又粗又长| 午夜激情欧美在线| 精品国内亚洲2022精品成人| av在线蜜桃| 国产精品一区二区三区四区久久| 可以在线观看的亚洲视频| 国产精品一区二区免费欧美| 日韩中文字幕欧美一区二区| 国产男人的电影天堂91| 男人和女人高潮做爰伦理| 亚洲在线观看片| 一级黄色大片毛片| 男人舔女人下体高潮全视频| 丰满乱子伦码专区| 亚洲国产色片| 精品一区二区三区av网在线观看| av天堂中文字幕网| 男女下面进入的视频免费午夜| av视频在线观看入口| 婷婷丁香在线五月| 少妇熟女aⅴ在线视频| 亚洲精品成人久久久久久| 男女之事视频高清在线观看| 麻豆精品久久久久久蜜桃| 日韩在线高清观看一区二区三区 | 免费av毛片视频| 最后的刺客免费高清国语| 久久人妻av系列| 亚洲成人免费电影在线观看| 69人妻影院| 国产色婷婷99| www.色视频.com| 国内精品久久久久久久电影| 人妻夜夜爽99麻豆av| 搡老妇女老女人老熟妇| 欧美潮喷喷水| 欧美又色又爽又黄视频| 国产伦人伦偷精品视频| 看片在线看免费视频| 成人av一区二区三区在线看| 制服丝袜大香蕉在线| 淫妇啪啪啪对白视频| 成人三级黄色视频| bbb黄色大片| 少妇高潮的动态图| 琪琪午夜伦伦电影理论片6080| 精品久久久久久成人av| 成人毛片a级毛片在线播放| 欧美zozozo另类| 亚洲在线自拍视频| 国产一级毛片七仙女欲春2| 亚洲av免费在线观看| 欧美在线一区亚洲| 麻豆av噜噜一区二区三区| netflix在线观看网站| 亚洲成人精品中文字幕电影| 免费观看人在逋| 婷婷丁香在线五月| 白带黄色成豆腐渣| 他把我摸到了高潮在线观看| 高清毛片免费观看视频网站| 亚洲人成网站高清观看| 国产女主播在线喷水免费视频网站 | 天堂av国产一区二区熟女人妻| 男女啪啪激烈高潮av片| 日本免费a在线| 婷婷丁香在线五月| 波野结衣二区三区在线| 国产成年人精品一区二区| 日本欧美国产在线视频| 精品一区二区三区视频在线| 亚洲国产欧洲综合997久久,| 18禁黄网站禁片午夜丰满| 国产欧美日韩精品一区二区| 国产高潮美女av| 国产黄色小视频在线观看| 欧美高清性xxxxhd video| 在线天堂最新版资源| 亚洲18禁久久av| 一级黄片播放器| 国产伦在线观看视频一区| 国产精品福利在线免费观看| 国产精品女同一区二区软件 | 在线a可以看的网站| 九九久久精品国产亚洲av麻豆| 18禁裸乳无遮挡免费网站照片| 久久久久久久亚洲中文字幕| 亚洲国产日韩欧美精品在线观看| 此物有八面人人有两片| 俄罗斯特黄特色一大片| 91精品国产九色| 黄色一级大片看看| 99热网站在线观看| 在现免费观看毛片| 免费观看精品视频网站| 久久热精品热| 午夜福利在线观看吧| 韩国av一区二区三区四区| 天堂√8在线中文| 欧美高清性xxxxhd video| 人人妻人人澡欧美一区二区| 一级黄片播放器| 国产午夜福利久久久久久| 少妇猛男粗大的猛烈进出视频 | 天堂网av新在线| 一a级毛片在线观看| 麻豆精品久久久久久蜜桃| 亚洲一区二区三区色噜噜| 亚洲精华国产精华精| 成人永久免费在线观看视频| 日韩强制内射视频| 看黄色毛片网站| 99热只有精品国产| 老女人水多毛片| 国内揄拍国产精品人妻在线| 蜜桃亚洲精品一区二区三区| 午夜福利在线在线| 简卡轻食公司| 老熟妇仑乱视频hdxx| 国产精品久久视频播放| 日韩精品中文字幕看吧| 欧美成人性av电影在线观看| 亚洲欧美日韩东京热| 亚洲精品国产成人久久av| 又黄又爽又刺激的免费视频.| 成人国产麻豆网| 成年免费大片在线观看| 自拍偷自拍亚洲精品老妇| 在现免费观看毛片| 国内精品美女久久久久久| videossex国产| 日韩精品中文字幕看吧| av视频在线观看入口| 国产精品亚洲美女久久久| 成人特级av手机在线观看| 亚洲在线自拍视频| 88av欧美| 精品久久久噜噜| 国产91精品成人一区二区三区| 国产精品久久电影中文字幕| 婷婷色综合大香蕉| 如何舔出高潮| 少妇被粗大猛烈的视频| 老司机福利观看| av中文乱码字幕在线| 亚洲人与动物交配视频| 男人舔奶头视频| 99精品久久久久人妻精品| 国产午夜福利久久久久久| 亚洲精品久久国产高清桃花| 99视频精品全部免费 在线| 天天躁日日操中文字幕| 网址你懂的国产日韩在线| 午夜亚洲福利在线播放| 国产精品99久久久久久久久| 一夜夜www| 欧美人与善性xxx| 亚洲无线在线观看| 91在线观看av| 日韩欧美免费精品| 日韩国内少妇激情av| 美女xxoo啪啪120秒动态图| 亚洲av中文字字幕乱码综合| 男女边吃奶边做爰视频| 波多野结衣巨乳人妻| 亚洲av成人av| 日韩大尺度精品在线看网址| 精品久久久久久久久av| a级一级毛片免费在线观看| 精品午夜福利在线看| 中文在线观看免费www的网站| 久久99热这里只有精品18| 99久国产av精品| 国产 一区精品| 桃色一区二区三区在线观看| 亚洲国产精品sss在线观看| 国产亚洲欧美98| 午夜免费成人在线视频| 麻豆精品久久久久久蜜桃| 国产精品久久久久久亚洲av鲁大| 国产日本99.免费观看| 午夜福利在线在线| 又爽又黄无遮挡网站| 国产午夜精品久久久久久一区二区三区 | 亚洲精华国产精华精| 亚洲av五月六月丁香网| 成人无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 日本 欧美在线| 国产精品98久久久久久宅男小说| 久久天躁狠狠躁夜夜2o2o| 波野结衣二区三区在线| 俺也久久电影网| 国产精品久久久久久久电影| 给我免费播放毛片高清在线观看| 亚洲中文字幕一区二区三区有码在线看| 嫁个100分男人电影在线观看| 国产视频一区二区在线看| 99riav亚洲国产免费| 亚洲成人免费电影在线观看| 1024手机看黄色片| 97碰自拍视频| 欧美中文日本在线观看视频| avwww免费| 午夜福利成人在线免费观看| 日本 av在线| 中文字幕免费在线视频6| 日韩精品有码人妻一区| 中文字幕人妻熟人妻熟丝袜美| 国产91精品成人一区二区三区| 国产大屁股一区二区在线视频| 国产真实伦视频高清在线观看 | 一区二区三区免费毛片| 欧美成人免费av一区二区三区| 99视频精品全部免费 在线| 91久久精品国产一区二区三区| 国产大屁股一区二区在线视频| 搡老岳熟女国产| 成人av在线播放网站| 国产乱人伦免费视频| 亚洲狠狠婷婷综合久久图片| av在线老鸭窝| eeuss影院久久| 无人区码免费观看不卡| 熟女电影av网| а√天堂www在线а√下载| 美女被艹到高潮喷水动态| 久久婷婷人人爽人人干人人爱| 99热这里只有是精品在线观看| 国产一区二区三区视频了| 日韩av在线大香蕉| 久久久久免费精品人妻一区二区| 精品久久久久久久久av| 久久人妻av系列| 日本 av在线| 国产单亲对白刺激| 99国产极品粉嫩在线观看| 美女高潮的动态| 97人妻精品一区二区三区麻豆| 欧美极品一区二区三区四区| 高清日韩中文字幕在线| 国产精品电影一区二区三区| 乱人视频在线观看| 国产精品一区二区三区四区免费观看 | 毛片一级片免费看久久久久 | 91久久精品国产一区二区成人| 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久久末码| 身体一侧抽搐| 无人区码免费观看不卡| 18禁黄网站禁片免费观看直播| 欧洲精品卡2卡3卡4卡5卡区| 深夜精品福利| 人人妻人人看人人澡| 日本免费a在线| 看十八女毛片水多多多| 成年人黄色毛片网站| 亚洲成人久久性| 热99re8久久精品国产| 欧美日本亚洲视频在线播放| 国产精品久久电影中文字幕| 亚洲自偷自拍三级| 国产精品人妻久久久久久| 制服丝袜大香蕉在线| 免费无遮挡裸体视频| 午夜免费激情av| 国产人妻一区二区三区在| 国产一级毛片七仙女欲春2| 国产人妻一区二区三区在| 亚洲真实伦在线观看| 亚洲男人的天堂狠狠| 日本成人三级电影网站| 国内精品美女久久久久久| or卡值多少钱| 成人欧美大片| 国产私拍福利视频在线观看| 欧美另类亚洲清纯唯美| 精品人妻视频免费看| 中文字幕精品亚洲无线码一区| 麻豆成人av在线观看| 深夜精品福利| 精品一区二区三区av网在线观看| 99久久中文字幕三级久久日本| 午夜精品久久久久久毛片777| 亚洲中文字幕一区二区三区有码在线看| 一个人观看的视频www高清免费观看| 伦精品一区二区三区| 九九热线精品视视频播放| 婷婷丁香在线五月| 精品人妻1区二区| 亚洲欧美激情综合另类| 免费av不卡在线播放| 一本久久中文字幕| 免费观看在线日韩| 嫁个100分男人电影在线观看| 久久久久国产精品人妻aⅴ院| 天美传媒精品一区二区| 成人永久免费在线观看视频| 听说在线观看完整版免费高清| 可以在线观看毛片的网站| 国产一区二区激情短视频| 久久久久久久久中文| 中文字幕熟女人妻在线| 在线播放无遮挡| 亚洲性夜色夜夜综合| 国产精品乱码一区二三区的特点| 啪啪无遮挡十八禁网站| 日本熟妇午夜| 久久久久久久午夜电影| 亚洲三级黄色毛片| 欧美bdsm另类| 成熟少妇高潮喷水视频| 91av网一区二区| xxxwww97欧美| 午夜福利在线在线| 18禁黄网站禁片午夜丰满| 成人av一区二区三区在线看| 国产人妻一区二区三区在| 12—13女人毛片做爰片一| 国产麻豆成人av免费视频| 久久久午夜欧美精品| 成年免费大片在线观看| 国产三级中文精品| 两人在一起打扑克的视频| av在线观看视频网站免费| 日本撒尿小便嘘嘘汇集6| ponron亚洲| 无人区码免费观看不卡| 精品人妻偷拍中文字幕| 一进一出抽搐动态| 此物有八面人人有两片| 国产成人av教育| 日韩一本色道免费dvd| 精品人妻熟女av久视频| 精品人妻1区二区| 亚洲av电影不卡..在线观看| 国产三级在线视频| 欧美色欧美亚洲另类二区| 亚洲av中文av极速乱 | 特级一级黄色大片| 综合色av麻豆| av专区在线播放| 久久精品影院6| 99精品在免费线老司机午夜| 日本免费a在线| 久久精品91蜜桃| 麻豆av噜噜一区二区三区| 亚洲av.av天堂| 亚洲午夜理论影院| 久久九九热精品免费| 中文字幕av成人在线电影| 亚洲avbb在线观看| 国产欧美日韩精品亚洲av| 亚洲va日本ⅴa欧美va伊人久久| 久久九九热精品免费| 精品人妻视频免费看| 亚洲人成网站高清观看| 天美传媒精品一区二区| 精品一区二区三区人妻视频| 免费看av在线观看网站| 精品午夜福利在线看| 两性午夜刺激爽爽歪歪视频在线观看| 一本一本综合久久| 欧美日韩精品成人综合77777| 午夜精品一区二区三区免费看| 校园人妻丝袜中文字幕| 色哟哟哟哟哟哟| 特大巨黑吊av在线直播| 色哟哟·www| h日本视频在线播放| 少妇人妻精品综合一区二区 | 在线观看66精品国产| 非洲黑人性xxxx精品又粗又长| 色综合亚洲欧美另类图片| 中文字幕高清在线视频| 亚洲黑人精品在线| 午夜福利高清视频| 日本 av在线| 国产精品人妻久久久影院| 国产色婷婷99| 久久久久久国产a免费观看| 国产亚洲精品av在线| а√天堂www在线а√下载| 白带黄色成豆腐渣| 国产成人影院久久av| 色哟哟哟哟哟哟| 中文字幕精品亚洲无线码一区| 色在线成人网| 九九爱精品视频在线观看| 亚洲一级一片aⅴ在线观看| 免费在线观看日本一区| 嫩草影院精品99| 又黄又爽又免费观看的视频| 一边摸一边抽搐一进一小说| 麻豆久久精品国产亚洲av| 国产成人aa在线观看| 国产精品无大码| 国产精品久久电影中文字幕| 国产精品久久久久久精品电影| 国产在线男女| 一区二区三区免费毛片| 午夜福利18| 免费在线观看日本一区| 性插视频无遮挡在线免费观看| 精品人妻1区二区| 国产精品一区二区免费欧美| av视频在线观看入口| 国产乱人视频| 中文在线观看免费www的网站| 欧美日韩乱码在线| 国内少妇人妻偷人精品xxx网站| 亚洲第一电影网av| 欧美日韩乱码在线| 久久热精品热| 国产又黄又爽又无遮挡在线| 亚洲七黄色美女视频| 亚洲 国产 在线| 观看免费一级毛片| 亚洲欧美日韩卡通动漫| 神马国产精品三级电影在线观看| 99热这里只有是精品50| 色播亚洲综合网| 97碰自拍视频| 男女那种视频在线观看| 99国产精品一区二区蜜桃av| 日日啪夜夜撸| 99热这里只有是精品50| 啪啪无遮挡十八禁网站| 国产精品福利在线免费观看| 啦啦啦韩国在线观看视频| 一边摸一边抽搐一进一小说| 亚洲精品日韩av片在线观看| 国产在线男女| 久久九九热精品免费| 热99在线观看视频| 亚洲,欧美,日韩| 乱系列少妇在线播放| 免费电影在线观看免费观看| 又黄又爽又免费观看的视频| 97超级碰碰碰精品色视频在线观看|