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

    FPS異常的2種不同解釋方法對比

    2015-10-13 11:59:37肖波劉海飛戴前偉
    關(guān)鍵詞:柱體點源交匯

    肖波,劉海飛,戴前偉

    ?

    FPS異常的2種不同解釋方法對比

    肖波1,劉海飛2,戴前偉2

    (1.中國能源建設(shè)集團 廣東省電力設(shè)計研究院,廣東 廣州,510600;2.中南大學(xué) 地球科學(xué)與信息物理學(xué)院,湖南 長沙,410083)

    從有限單元法正演模擬入手,采用最小二乘法進行反演計算,通過VC++程序?qū)崿F(xiàn),對多個地電模型實例進行計算,并將最小二乘法反演結(jié)果與圓弧交匯?相對強度計算結(jié)果進行對比。研究結(jié)果表明:圓弧交匯?相對強度法對簡單模型適用性高,由其反演成果所確定異常體中心位置甚至?xí)?yōu)于最小二乘法最優(yōu)化反演結(jié)果;最小二乘法能適應(yīng)各種地電條件對多種模型進行反演計算,其計算結(jié)果優(yōu)于圓弧交匯?相對強度法計算結(jié)果。

    電阻率;極化率;數(shù)值模擬;有限元法;反演;定量解釋

    固定點源激電測深法(FPS)又稱固定點源雙邊三極測深法。該方法裝置在實際工作中采用多臺接受機同時測量,測量電極極距根據(jù)實際情況可靈活變化,導(dǎo)體極化場位置隨著供電點位置不同而變化,具有較高工作效率的同時可更好地確定異常體的位置。但該裝置在數(shù)據(jù)處理過程中,常規(guī)的相對強度?圓弧交匯解釋方法受場源的分布位置影響大,其解釋成果對模型形狀、埋深半徑比的依賴性較大;該方法對于單個模型的處理解釋通常能取得較好效果,但對組合模型的解釋通常會產(chǎn)生明顯的假異常,從而掩蓋了真實異常體的信息,為勘探解釋工作帶來了極大的困難。在此,本文作者從有限元正演模擬和最小二乘反演方法的角度入手,對多個地電模型進行實例計算,并將最小二乘反演結(jié)果與圓弧交匯?相對強度計算結(jié)果進行對比。

    1 正演模擬

    采用有限單元法進行正演模擬,將整個計算空間區(qū)域離散成許多相互連接的單元網(wǎng)格,通過對單元網(wǎng)格節(jié)點的電位進行求解得到計算域的電位分布,然后將電位分布值轉(zhuǎn)化為視電阻率,再通過“等效電阻率”求取視極化率。

    在二維地電條件下,點電流源場各節(jié)點電位的計算可歸結(jié)為對若干個給定波數(shù)求解電位的傅里葉變換(,,)所滿足的二維偏微分方程的邊值問題:

    式中:為電導(dǎo)率;和1為地面邊界;∞為其他邊界;I為第個供電電極的電流;為場源到測點的距離;邊界處法線向量;和分別為零階和一階修正貝塞爾函數(shù)[1?2]。

    與方程式(1)等價的變分問題為

    通過有限單元法數(shù)值模擬可得地表任意一點的電位,從而得到相應(yīng)的視電阻率[3]:

    其中:為裝置系數(shù);為供電電流;Δ為兩節(jié)點之間的電勢差。

    當(dāng)?shù)叵陆橘|(zhì)有極化特性時,其等效電阻率為

    式中:為介質(zhì)電阻率;為介質(zhì)極化率。利用上面所述的電阻率二維正演方法對等效電阻率進行正演,可以得到等效視電阻率:

    由式(5)可得到視極化率η

    2 解釋方法

    2.1 常規(guī)定量解釋方法

    電法勘探中,有2個或者2個以上的供電點,分別以供電點到視電阻率s曲線、視極化率s曲線極值點的距離為半徑畫圓弧,根據(jù)圓弧的交點確定異常體的中心位置,稱為圓弧交匯法[3?4],如圖1所示。

    圖1 圓弧交匯法示意圖

    相對強度法是將二維地電斷面網(wǎng)格化,以各網(wǎng)格節(jié)點到點源的距離為半徑作圓弧,取圓弧與水平地表交點位置對應(yīng)的視電阻率、視極化率為該網(wǎng)格節(jié)點的值,位于點源坐標(biāo)左邊的網(wǎng)格節(jié)點由視電阻率、視極化率左支實測值確定,位于點源坐標(biāo)右邊的網(wǎng)格節(jié)點由視電阻率、視極化率右支實測值確定。對于多個點源,先對各節(jié)點上求和,進而計算出他們的平均值找出斷面上最大的平均值,用各網(wǎng)格節(jié)點的平均值去除以,最后得到了相對強度;對于低阻體,提取斷面上最小平均值,然后用除以各網(wǎng)格節(jié)點的平均值,從而突出斷面異常體的中心位置。相對強度法的作圖示意圖如圖2所示。

    圖2 相對強度法的成圖示意圖

    2.2 最小二乘反演法

    采用基于圓滑約束最小二乘反演方法,利用正演模型和實測數(shù)據(jù)構(gòu)造一目標(biāo)函數(shù),并使其達到極小。圓滑約束最小二乘法是基于以下方程[5?8]:

    對方程組(7)進行求解,得到模型參數(shù)修正矢量Δ,將其代入式(8)便得到新的預(yù)測模型參數(shù)矢量(k):

    如此重復(fù),直至模擬數(shù)據(jù)和實測數(shù)據(jù)之間的平均均方誤差滿足要求為止[9?14]。平均均方誤差為

    在反演過程中,為提高求解大型、病態(tài)線性方程組的計算速度和反演效率,采用廣義共軛梯度算法(GCG)[15?16]對偏導(dǎo)數(shù)矩陣進行求解,其基本公式為

    其實現(xiàn)步驟如下。

    1) 初始向量(0)=0,(0)=,(0)=T(0)=(0);誤差限,阻尼因子;

    如果||(0)||≤,則終止,否則轉(zhuǎn)入步驟3);對于=2,3,…,max,計算到步驟10);

    8) 如果||(j+1)||≤終止, 否則轉(zhuǎn)入步驟9) ;

    3 模型計算

    采用有限單元法數(shù)值計算進行正演模擬,分別用圓弧交匯?相對強度法成圖反演和最小二乘最優(yōu)化反演進行反演計算。對多個模型實例的計算結(jié)果進行對比分析。

    綜合考慮數(shù)據(jù)量、反演效果及工作量之間的關(guān)系,擬設(shè)7個測深供電點:1=40 m,2=60 m,3=80 m,4=100 m,5=120 m,6=140 m,7=160 m。為達到有效勘探深度,勘探區(qū)域為?50~250 m,取0~200 m進行剖面成圖,反演結(jié)果剖面圖中方框內(nèi)為模型投影圖。

    3.1 正方柱體模型Ⅰ

    模型Ⅰ為長×寬為20 m×20 m低阻、高極化正方柱體(見圖3),中心位置為(100 m,?20 m),電阻率為1=50 Ω?m,極化率為1=5.0%;背景電阻率為0=100 Ω?m,極化率為0=1.0%。正演模擬各點源對應(yīng)的視極化率曲線如圖4所示。視極化率相對強度如圖5所示。視極化率最小二乘反演結(jié)果如圖6所示。

    圖3 正方柱體模型示意圖

    1—A1;2—A2;3—A3;4—A4;5—A5;6—A6;7—A7

    圖5 正方柱體模型視極化率相對強度

    圖6 正方柱體模型視極化率最小二乘反演結(jié)果

    對比圖5和6可知:最小二乘反演、圓弧交匯?相對強度法能通過視極化率反演成圖能較好的體現(xiàn)出正方柱體的位置。視極化率相對強度法計算結(jié)果能較準(zhǔn)確地反映出正方柱體的中心埋深位置;而最小二乘反演結(jié)果未能體現(xiàn)模型極化體中心,反演極化體中心位于實際模型中上部,但異常體處于模型投影范圍內(nèi)。

    3.2 正方柱體模型Ⅱ

    模型Ⅱ為長×寬為20 m×20 m正方低阻、高極化柱狀體(見圖7),中心位置為(40 m,?17 m),電阻率為1=50 Ω?m,極化率為1=5.0%;背景電阻率為0=100 Ω?m,極化率為0=1.0%。正演模擬各點源對應(yīng)的視極化率曲線如圖8所示。視極化率相對強度如圖9所示。視極化率最小二乘反演結(jié)果如圖10所示。

    圖7 正方柱體模型示意圖

    1—A1;2—A2;3—A3;4—A4;5—A5;6—A6;7—A7

    圖9 正方柱體模型正演視極化率相對強度圖

    圖10 正方柱體模型視極化率最小二乘反演結(jié)果

    由圖9可知:視極化率的相對強度中心與模型實際位置中心較吻合,但相對強度等值線不能確定模型的基本形態(tài)。相對強度中心周邊等值線向右成圓弧形傾斜,其原因是采用圓弧交匯法,與點源距離相等的各網(wǎng)格節(jié)點視電阻率、視極化率的取值相同,當(dāng)點源分布相對模型中心不對稱時,各網(wǎng)格節(jié)點視電阻率、視極化率難以均衡,從而導(dǎo)致相對強度中心周邊等值線呈圓弧形發(fā)生明顯傾斜。

    圖10表明:最小二乘反演能基本反映了模型的形態(tài)和位置;反演結(jié)果異常體頂面與實際模型頂面較吻合。但異常體中心位于實際模型中的中上部,模型的底界面無法通過異常體形態(tài)確定。圖10還表明:最小二乘反演結(jié)果受點源位置分布影響小。

    3.3 水平板狀體模型

    模型為長×寬為40 m×40 m水平低阻、高極化板狀體(見圖11),中心位置為(100 m,?10 m),電阻率1=50 Ω?m,極化率1=5.0%;背景電阻率0=100 Ω?m,極化率0=1.0%。正演模擬各點源對應(yīng)的視極化率曲線如圖12所示。視極化率相對強度如圖13所示。視極化率最小二乘反演結(jié)果如圖14所示。

    圖11 水平板狀體模型示意圖

    1—A1;2—A2;3—A3;4—A4;5—A5;6—A6;7—A7

    圖13 水平板狀體模型正演視極化率相對強度

    圖14 水平板狀體模型視極化率最小二乘反演結(jié)果

    對比圖13和圖14可知:最小二乘反演能很好的反映出水平低阻、極化板狀體的埋深、基本形態(tài),異常體位置與模型實際位置對應(yīng)關(guān)系較好;圓弧交匯?相對強度法能基本確定水平低阻、高極化板狀體位置,但視極化率相對強度中心位于低阻極化體模型位置正下部邊緣處,從視極化率相對強度等值線圖無法確定模型的基本形態(tài)。

    3.4 組合模型

    模型為2個長×寬為20 m×20 m正方柱形低阻、高極化模型(見圖15),中心位置分別為(80 m,?18 m)、(120 m,?18 m),電阻率1=50 Ω?m,極化率1=5.0%;背景電阻率0=100 Ω?m,極化率0=1.0%。正演模擬各點源對應(yīng)的視極化率曲線如圖16所示,視極化率相對強度如圖17所示,視極化率最小二乘反演結(jié)果如圖18所示。

    圖15 組合模型示意圖

    1—A1;2—A2;3—A3;4—A4;5—A5;6—A6;7—A7

    圖17 組合模型正演視極化率相對強度

    圖18 組合模型視極化率最小二乘反演結(jié)果

    圖17表明:組合模型視極化率相對強度中心與模型的實際位置中心相吻合;但模型的中垂線下方位置存在1個明顯的假異常,其原因是圓弧交匯法和相對強度法是建立在單個模型之上。對于組合模型,由于同一點源計算的視極化率曲線會產(chǎn)生2個極值點,從而進行圓弧交匯時就會產(chǎn)生其他相對強度,即為假異常,假異常通常隨著模型間的距離增大而變深。

    圖18表明:最小二乘視極化率反演較好地反映了組合模型的實際位置,反演結(jié)果與實際模型較吻合,無明顯假異常存在,但反演結(jié)果難以反映出模型的底界面。

    后期對高阻、高極化模型進行模型計算對比,其對比結(jié)論與低阻、高極化模型計算對比結(jié)果一致:最小二乘法反演結(jié)果優(yōu)于圓弧交匯?相對強度法。

    4 結(jié)論

    1) 固定點源激電測深具有簡單、高效的特點,其解釋方法圓弧交匯?相對強度法與最小二乘反演都能較準(zhǔn)確地反映出異常體中心位置,這2種解釋方法直觀且其解釋結(jié)果與實際結(jié)果較符合。

    2) 圓弧交匯?相對強度法對于單個模型能準(zhǔn)確地反映出模型的中心埋深位置;但該解釋方法受場源位置影響大;對于稍復(fù)雜模型未能精確的體現(xiàn)出異常體的位置及形狀,對于組合模型通常會產(chǎn)生較明顯的假異常。

    3) 最小二乘反演不受場源位置影響,能較好地反映出復(fù)雜模型的實際位置及輪廓,在對組合模型進行反演解釋時候,不會產(chǎn)生假異常;但該解釋方法難以體現(xiàn)模型的下底面位置。

    [1] 李小康. 時間域激電二維正演研究[D]. 北京: 中國地質(zhì)大學(xué)地球物理與信息技術(shù)學(xué)院, 2008: 27?47.
    LI Xiaokang. 2D forward numerical simulation of ip method in time domain[D]. Beijing: China University of Geosciences. School of Geophysics and Geoinformation Systems, 2008: 27?47.

    [2] 戴前偉, 肖波, 馮德山, 等. 基于二維高密度電阻率勘探數(shù)據(jù)的三維反演及應(yīng)用[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2012, 43(1): 293?310.
    DAI Qianwei, XIAO Bo, FENG Deshan, et al. 3-D inversion of the high density resistivity method based on 2-D exploration data and its application[J]. Journal of Central South University (Science and Technology), 2012, 43(1): 293?310.

    [3] 劉國興. 電法勘探原理與方法[M]. 北京: 地質(zhì)出版社, 2005: 111?128.
    LIU Guoxing, Electrical prospecting principle and method[M]. Beijing: Geological Publishing House, 2005: 111?128.

    [4] 李金銘. 地電場與電法勘探[M]. 北京: 地質(zhì)出版社, 2005: 258?260.
    LI Jinming. Electric field and electrical prospecting[M]. Beijing: Geological Publishing House, 2005: 258?260.

    [5] 李金銘, 魏文博, 陳本池, 等. 固定點源測深法定量解釋研究[J]. 物探與化探, 1997, 21(3): 187?196.
    LI Jinming, WEI Wenbo, CHEN Benchi, et al. The quantitative interpretation of the fixed point source sounding methods[J]. Geophysical & Geochemical Exploration, 1997, 21(3): 187?196.

    [6] 陳本池, 李金銘, 魏文博. 起伏地形條件下固定點源測深法定量解釋研究[J]. 現(xiàn)代地質(zhì), 1998, 12(1): 123?129.
    CHEN Benchi, LI Jinming, WEI Wenbo. A study of quantitative interpretation for fixed point source sounding method under undulatory terrai[J]. Geoscience, 1998, 12(1): 123?129.

    [7] 強建科, 阮百堯. 用FPS定量解釋法實現(xiàn)三維近似成像的探討[J]. CT理論與應(yīng)用研究, 2003, 12(1): 13?16.
    QIANG Jianke, RUAN Bairao. The study of 3-D IP anomaly by fixed point source sounding[J]. Study on the theory and Application of CT, 2003, 12(1): 13?16.

    [8] 強建科, 羅延鐘, 熊彬. 固定點源測深激電異常研究[J]. 地球物理學(xué)進展, 2005, 20(4): 1176?1183.
    QIANG Jianke, LUO Yanzhong, XIONG Bin. Study on the abnormal fixed point source sounding IP[J]. Progress in Geophysics, 2005, 20(4): 1176?1183.

    [9] 熊彬, 阮百堯, 黃俊革. 直流電阻率測深中二維反演程序?qū)θS數(shù)據(jù)的近似解釋[J]. 地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報), 2003, 28(1): 102?106.
    XIONG Bin, RUAN Baiyao, HUANG Junge. Approximate interpretation of 3-D data using 2-D inversion program in the DC resistivity sounding[J]. Earth Science (Journal of China University of Geosciences), 2003, 28(1): 102?106.

    [10] 黃俊革, 阮百堯. 直流電阻率測深中二維與三維反演結(jié)果的對比與分析[J]. 物探與化探, 2004, 28(5): 447?450.
    HUANG Junge, RUAN Baiyao. An analytical comparison between 2D and 3D inversions in dc resistivity sounding[J]. Geophysical & Geochemical Explopation, 2004, 28(5): 447?450.

    [11] 黃俊革, 王家林, 阮百堯. 三維高密度電阻率E-SCAN法有限元模擬異常特征研究[J]. 地球物理學(xué)報, 2006, 49(4): 1206?1214.
    HUANG Junge, WANG Jialin, RUAN Baiyao. A study on FEM modeling of anomalies of 3-D high-density E-SCAN resistivity survey[J]. Chinese Journal of Geophysics, 2006, 49(4): 1206?1214.

    [12] 劉海飛. 高密度數(shù)據(jù)處理方法研究[D]. 長沙: 中南大學(xué)信息物理工程學(xué)院, 2003: 30?60.
    LIU Haifei. The processing method of High-density data[D]. Changsha: Central South University. School of Info-physics and Geomatics Engineering, 2003: 30?60.

    [13] 黃俊革. 三維電阻率/極化率有限元正演模擬與反演成像[D]. 長沙: 中南大學(xué)信息物理工程學(xué)院, 2003: 20?60.
    HUANG Junge. 3-Dresistivity/IP modeling and inversion based on FEM[D]. Changsha: Central South University. School of Info-physics and Geomatics Engineering, 2003: 20?60.

    [14] 張大海, 王興泰. 二維視電阻率斷面的快速最小二乘反演[J]. 物探化探計算技術(shù), 1999, 21(1): 1?12.
    ZHANG Dahai, WANG Xingtai. Rapd least-squares inversion of 2-D apparent resistivity pseudosection[J]. Computing Techniques for Geophysical And Geochemical Exploration, 1999, 21(1): 1?12.

    [15] 吳小平, 徐果明. 利用共軛梯度法的電阻率三維反演研究[J]. 地球物理學(xué)報, 2000, 43(3): 420?427.
    WU Xiaoping, XU Guoming. Study 3-D resistivity inversion using conjugate gradient method[J]. Chinese Journal of Geophysics, 2000, 43(3): 420?427.

    [16] 周竹生. 廣義共軛梯度算法[J]. 物探與化探, 1996, 20(5): 351?358.
    ZHOU Zhusheng. The generalized conjugate gradient algorithm[J]. Geophysical and Geochemical Exploration, 1996, 20(5): 351?358.

    Two different interpretation methods about FPS anomaly

    XIAO Bo1, LIU Haifei2, DAI Qianwei2

    (1. Guangdong Electric Power Design Institute, China Energy Engineering Group, Guangzhou 510600, China;2. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China)

    The plurality of geoelectric model instance was calculated based on the forward modeling of the finite element method, using the least square method and the program with VC++. Then the least square inversion results and arc intersection?relative intensity calculation results were compared. The results show that arc intersection?relative intensity method has high applicability in simple model, the inverse mapping effect is even superior to the least square optimization inversion results in determining the center of abnormal body position. The least square inversion method can adapt various geoelectric conditions on a variety of complex model inversion, and the inversion effects of the complex model and combined model are superior to the calculation results of the arc intersection?relative intensity method.

    resistivity; polarizability; numerical simulation; finite element method; inversion; quantitative interpretation

    P631

    A

    1672?7207(2015)02?0595?08

    2014?02?13;

    2014?04?20

    國家自然科學(xué)基金資助項目(41374118,41174102,41074085);教育部博士點基金資助項目(20120162110015)(Project (41374118, 41174102, 41074085) supported by the National Natural Science Foundation of China; Project (20120162110015) supported by the Doctoral Fund of Ministry of Education of China)

    戴前偉,博士,教授,從事電磁法方法及理論、工程地球物理勘探等研究;E-mail:qwdai@mail.csu.edu.cn

    10.11817/j.issn.1672-7207.2015.02.030

    (編輯 趙俊)

    猜你喜歡
    柱體點源交匯
    不同倒角半徑四柱體繞流數(shù)值模擬及水動力特性分析
    海洋工程(2021年1期)2021-02-02 02:48:12
    兩水交匯,涇渭分明
    三角的交匯創(chuàng)新問題
    基于多介質(zhì)ALE算法的柱體高速垂直入水仿真
    關(guān)于脈沖積累對雙點源干擾影響研究
    聚焦離散型隨機變量的交匯創(chuàng)新
    靜止軌道閃電探測性能實驗室驗證技術(shù)研究
    談擬柱體的體積
    外注式單體液壓支柱頂蓋與活柱體連接結(jié)構(gòu)的改進
    基于標(biāo)準(zhǔn)化點源敏感性的鏡面視寧度評價
    热99re8久久精品国产| 国产淫片久久久久久久久| 久久综合国产亚洲精品| 九色成人免费人妻av| 久久久久久伊人网av| 又粗又硬又长又爽又黄的视频 | 日韩人妻高清精品专区| 午夜福利在线观看免费完整高清在 | 欧美日韩一区二区视频在线观看视频在线 | 麻豆精品久久久久久蜜桃| 精品人妻熟女av久视频| 欧美性猛交黑人性爽| 久久婷婷人人爽人人干人人爱| 久久精品人妻少妇| 在线观看av片永久免费下载| 2022亚洲国产成人精品| 神马国产精品三级电影在线观看| 亚洲美女视频黄频| 成人特级av手机在线观看| 超碰av人人做人人爽久久| 亚洲高清免费不卡视频| 嫩草影院新地址| 午夜福利在线观看吧| 亚洲最大成人手机在线| 午夜爱爱视频在线播放| 国产v大片淫在线免费观看| 高清在线视频一区二区三区 | 亚洲av一区综合| 91久久精品国产一区二区成人| 久久午夜福利片| 1024手机看黄色片| 亚洲精品国产av成人精品| 国产国拍精品亚洲av在线观看| 国产老妇女一区| 真实男女啪啪啪动态图| 99久久精品国产国产毛片| 99在线人妻在线中文字幕| avwww免费| 日韩精品有码人妻一区| 一进一出抽搐gif免费好疼| 久久亚洲国产成人精品v| 中文字幕久久专区| 日韩大尺度精品在线看网址| 麻豆一二三区av精品| 午夜精品一区二区三区免费看| 国产又黄又爽又无遮挡在线| 久久人妻av系列| 精品日产1卡2卡| 日本在线视频免费播放| 欧美极品一区二区三区四区| 看片在线看免费视频| 看免费成人av毛片| 神马国产精品三级电影在线观看| 男女做爰动态图高潮gif福利片| 男人舔奶头视频| 久久99蜜桃精品久久| 最新中文字幕久久久久| 亚洲熟妇中文字幕五十中出| 最近视频中文字幕2019在线8| 免费不卡的大黄色大毛片视频在线观看 | 九草在线视频观看| 一本一本综合久久| 久久综合国产亚洲精品| 丰满的人妻完整版| 久久精品国产亚洲av香蕉五月| 丰满的人妻完整版| 日韩欧美 国产精品| 欧美激情国产日韩精品一区| 亚洲精品亚洲一区二区| 久久99蜜桃精品久久| 国产成人一区二区在线| 亚洲国产精品国产精品| 你懂的网址亚洲精品在线观看 | 久久久久久久久久久免费av| 国产成人freesex在线| 国产色爽女视频免费观看| 少妇熟女欧美另类| 熟女人妻精品中文字幕| 97热精品久久久久久| 两个人视频免费观看高清| 色哟哟·www| 国产一区二区三区在线臀色熟女| 干丝袜人妻中文字幕| 国产成人午夜福利电影在线观看| 中国美白少妇内射xxxbb| 伦精品一区二区三区| 欧美最黄视频在线播放免费| 老司机福利观看| 精品人妻熟女av久视频| 麻豆精品久久久久久蜜桃| 国产中年淑女户外野战色| 国内久久婷婷六月综合欲色啪| 免费av不卡在线播放| 精品一区二区免费观看| 老司机影院成人| 99久久人妻综合| 在线播放无遮挡| 天堂av国产一区二区熟女人妻| 伦精品一区二区三区| 国产女主播在线喷水免费视频网站 | 婷婷色综合大香蕉| 国内揄拍国产精品人妻在线| 久久中文看片网| 精华霜和精华液先用哪个| 国产精品福利在线免费观看| 亚洲美女视频黄频| 99久国产av精品国产电影| 精品久久久久久久久久免费视频| 欧美+亚洲+日韩+国产| a级毛片a级免费在线| 中文字幕精品亚洲无线码一区| 日韩欧美在线乱码| 一区福利在线观看| 99久国产av精品国产电影| 亚洲丝袜综合中文字幕| 91av网一区二区| 夜夜爽天天搞| 日本一本二区三区精品| 亚洲精品亚洲一区二区| 免费搜索国产男女视频| 精品久久久噜噜| 一进一出抽搐gif免费好疼| 国产一区二区三区在线臀色熟女| 亚洲不卡免费看| 欧美日韩综合久久久久久| 国产一区亚洲一区在线观看| 国产人妻一区二区三区在| 麻豆一二三区av精品| 久久人人精品亚洲av| 亚洲av第一区精品v没综合| 免费黄网站久久成人精品| 变态另类成人亚洲欧美熟女| 美女黄网站色视频| 麻豆精品久久久久久蜜桃| 伊人久久精品亚洲午夜| 久久亚洲精品不卡| 久久九九热精品免费| 国产在线男女| 日韩精品有码人妻一区| 久久午夜亚洲精品久久| 中文资源天堂在线| 91麻豆精品激情在线观看国产| 给我免费播放毛片高清在线观看| 成人永久免费在线观看视频| 99热网站在线观看| 99热全是精品| 一级黄色大片毛片| 欧美性猛交黑人性爽| 中出人妻视频一区二区| 天堂中文最新版在线下载 | 国产v大片淫在线免费观看| 波野结衣二区三区在线| 久久精品综合一区二区三区| 日本在线视频免费播放| 欧美又色又爽又黄视频| 亚洲七黄色美女视频| 老师上课跳d突然被开到最大视频| 日韩欧美在线乱码| 国产 一区精品| 搡老妇女老女人老熟妇| 国产中年淑女户外野战色| 特级一级黄色大片| 亚州av有码| 日韩 亚洲 欧美在线| 欧美高清成人免费视频www| 亚洲精品日韩av片在线观看| 精品一区二区三区视频在线| 成人高潮视频无遮挡免费网站| 菩萨蛮人人尽说江南好唐韦庄 | 一进一出抽搐动态| 变态另类成人亚洲欧美熟女| 激情 狠狠 欧美| 欧美极品一区二区三区四区| 亚洲一区二区三区色噜噜| 久久久精品欧美日韩精品| 女的被弄到高潮叫床怎么办| 日韩一本色道免费dvd| 日韩成人av中文字幕在线观看| 精品久久久久久久久久久久久| 天堂中文最新版在线下载 | av黄色大香蕉| av在线观看视频网站免费| 欧美成人a在线观看| 国产成人精品一,二区 | 色尼玛亚洲综合影院| 噜噜噜噜噜久久久久久91| 国产精品一二三区在线看| 国产精品.久久久| 波多野结衣高清作品| 免费搜索国产男女视频| 乱系列少妇在线播放| 精品99又大又爽又粗少妇毛片| 国产伦理片在线播放av一区 | 我要搜黄色片| 精品久久久噜噜| 国产一区二区三区av在线 | 日本免费一区二区三区高清不卡| 99久久无色码亚洲精品果冻| 免费av观看视频| 99久久精品热视频| 99热全是精品| 成人亚洲欧美一区二区av| 精品久久国产蜜桃| 国产精品久久电影中文字幕| 自拍偷自拍亚洲精品老妇| 美女内射精品一级片tv| 国产伦理片在线播放av一区 | 干丝袜人妻中文字幕| 高清午夜精品一区二区三区 | 人人妻人人澡人人爽人人夜夜 | 成人毛片a级毛片在线播放| 观看美女的网站| 国产综合懂色| 国产精品无大码| 久久久久九九精品影院| av在线亚洲专区| 免费电影在线观看免费观看| 中出人妻视频一区二区| 男人舔奶头视频| 亚洲精品久久国产高清桃花| 国内精品一区二区在线观看| 综合色av麻豆| 色尼玛亚洲综合影院| 波多野结衣高清无吗| 一边摸一边抽搐一进一小说| 看免费成人av毛片| 日韩强制内射视频| 丝袜喷水一区| 男女做爰动态图高潮gif福利片| 狂野欧美白嫩少妇大欣赏| 欧美不卡视频在线免费观看| 亚洲最大成人av| 久久精品国产鲁丝片午夜精品| 波多野结衣高清无吗| 午夜老司机福利剧场| 亚洲成av人片在线播放无| 国产精品久久久久久精品电影| 久久精品夜色国产| 精品熟女少妇av免费看| 色播亚洲综合网| 日韩一区二区视频免费看| 在线观看一区二区三区| 非洲黑人性xxxx精品又粗又长| 日本成人三级电影网站| 丝袜美腿在线中文| 国产伦一二天堂av在线观看| 亚洲精品粉嫩美女一区| 国产日韩欧美在线精品| 九九久久精品国产亚洲av麻豆| 少妇人妻精品综合一区二区 | 91在线精品国自产拍蜜月| 联通29元200g的流量卡| 美女国产视频在线观看| 亚洲久久久久久中文字幕| 欧美又色又爽又黄视频| 亚洲欧美中文字幕日韩二区| 天堂中文最新版在线下载 | 精品久久久久久久人妻蜜臀av| 菩萨蛮人人尽说江南好唐韦庄 | 久久精品国产亚洲av天美| 男女啪啪激烈高潮av片| 久久国内精品自在自线图片| 女同久久另类99精品国产91| 丝袜美腿在线中文| 国产伦一二天堂av在线观看| 国产精品爽爽va在线观看网站| 级片在线观看| 国产伦精品一区二区三区视频9| 亚洲成a人片在线一区二区| 久久精品91蜜桃| 日本爱情动作片www.在线观看| 可以在线观看毛片的网站| 最新中文字幕久久久久| 亚洲欧美精品综合久久99| 亚洲精品日韩在线中文字幕 | 国产亚洲5aaaaa淫片| 午夜精品在线福利| 人人妻人人澡人人爽人人夜夜 | 国产亚洲av嫩草精品影院| 亚洲婷婷狠狠爱综合网| 免费一级毛片在线播放高清视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品国产亚洲av香蕉五月| 在现免费观看毛片| 日韩成人伦理影院| 午夜a级毛片| 99riav亚洲国产免费| 欧洲精品卡2卡3卡4卡5卡区| 成年免费大片在线观看| av福利片在线观看| 舔av片在线| 婷婷色av中文字幕| 国产成人freesex在线| av又黄又爽大尺度在线免费看 | 成人毛片60女人毛片免费| 国产一区二区在线av高清观看| 九九爱精品视频在线观看| 久久99精品国语久久久| 韩国av在线不卡| 色哟哟哟哟哟哟| 性插视频无遮挡在线免费观看| 欧美激情国产日韩精品一区| 我要搜黄色片| 国产伦在线观看视频一区| 变态另类丝袜制服| 欧美最黄视频在线播放免费| 又粗又硬又长又爽又黄的视频 | 欧美又色又爽又黄视频| 美女 人体艺术 gogo| av在线蜜桃| 97热精品久久久久久| 色综合色国产| av视频在线观看入口| 91狼人影院| 九九热线精品视视频播放| 国产成人精品一,二区 | 成人永久免费在线观看视频| av在线天堂中文字幕| 精品久久久久久久人妻蜜臀av| 亚洲国产欧美在线一区| 精品久久久噜噜| 免费电影在线观看免费观看| 观看免费一级毛片| 久久久久久久久久黄片| 亚洲在线自拍视频| 国产亚洲91精品色在线| 99精品在免费线老司机午夜| 亚洲av成人av| 国产探花在线观看一区二区| 一边亲一边摸免费视频| 大又大粗又爽又黄少妇毛片口| 国产精品,欧美在线| 亚洲久久久久久中文字幕| 久久这里有精品视频免费| 久久精品久久久久久噜噜老黄 | 国产午夜精品论理片| 亚洲欧洲日产国产| 老师上课跳d突然被开到最大视频| 午夜福利在线观看免费完整高清在 | 亚洲av中文字字幕乱码综合| 亚洲国产高清在线一区二区三| 日韩欧美三级三区| 国产精品美女特级片免费视频播放器| 可以在线观看毛片的网站| 99国产极品粉嫩在线观看| 天天躁日日操中文字幕| 久久久a久久爽久久v久久| 国产极品精品免费视频能看的| 亚洲无线在线观看| 特大巨黑吊av在线直播| 国产精品.久久久| 99在线人妻在线中文字幕| 狂野欧美白嫩少妇大欣赏| 国产午夜精品论理片| 丰满的人妻完整版| 午夜精品在线福利| 亚洲,欧美,日韩| 丰满的人妻完整版| 午夜激情福利司机影院| 婷婷亚洲欧美| 一级毛片久久久久久久久女| 国产精品乱码一区二三区的特点| 亚洲七黄色美女视频| 国产精品乱码一区二三区的特点| 69av精品久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区四区激情视频 | 亚洲高清免费不卡视频| 国产老妇女一区| 日韩欧美精品免费久久| 亚洲中文字幕一区二区三区有码在线看| 久久久久九九精品影院| 三级国产精品欧美在线观看| 亚洲国产色片| 中文字幕av成人在线电影| 亚洲成人久久爱视频| 一进一出抽搐gif免费好疼| 精品人妻熟女av久视频| 国产亚洲5aaaaa淫片| 麻豆精品久久久久久蜜桃| 人人妻人人澡欧美一区二区| 亚洲一区高清亚洲精品| 少妇被粗大猛烈的视频| 99久久成人亚洲精品观看| 国产精品av视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| av国产免费在线观看| 一级av片app| 精品人妻视频免费看| 国产精品久久久久久久久免| 亚洲美女视频黄频| 久久久a久久爽久久v久久| 亚洲av免费在线观看| av在线老鸭窝| 亚洲国产欧美人成| 亚洲无线观看免费| 两个人视频免费观看高清| 婷婷亚洲欧美| 亚洲va在线va天堂va国产| 日本-黄色视频高清免费观看| 久久午夜亚洲精品久久| 成人毛片a级毛片在线播放| 可以在线观看的亚洲视频| 国产爱豆传媒在线观看| 99热这里只有是精品在线观看| 黄片无遮挡物在线观看| 国产91av在线免费观看| 久久中文看片网| 国内精品久久久久精免费| 亚洲av成人av| 亚洲精华国产精华液的使用体验 | 国产高清不卡午夜福利| 午夜精品在线福利| 日韩精品青青久久久久久| 午夜激情欧美在线| 97人妻精品一区二区三区麻豆| 免费看日本二区| 久久精品久久久久久噜噜老黄 | 国产高清不卡午夜福利| 精品欧美国产一区二区三| 成人毛片60女人毛片免费| 国产麻豆成人av免费视频| 亚洲精品日韩在线中文字幕 | 美女国产视频在线观看| 青春草视频在线免费观看| 欧美一区二区国产精品久久精品| 欧美人与善性xxx| 亚洲aⅴ乱码一区二区在线播放| 少妇高潮的动态图| 国产精华一区二区三区| 亚洲精品亚洲一区二区| 如何舔出高潮| 国产精品免费一区二区三区在线| 女人被狂操c到高潮| 全区人妻精品视频| 午夜爱爱视频在线播放| 中文欧美无线码| 成人亚洲精品av一区二区| 亚洲美女视频黄频| 日韩视频在线欧美| 一级二级三级毛片免费看| 国内精品一区二区在线观看| 亚洲精品国产av成人精品| 悠悠久久av| 麻豆乱淫一区二区| 国产真实伦视频高清在线观看| 国产v大片淫在线免费观看| 国产伦精品一区二区三区视频9| 十八禁国产超污无遮挡网站| 男女边吃奶边做爰视频| 日韩国内少妇激情av| 亚洲婷婷狠狠爱综合网| 边亲边吃奶的免费视频| 99热这里只有是精品在线观看| 欧美日韩精品成人综合77777| 尤物成人国产欧美一区二区三区| 亚洲内射少妇av| 久久亚洲国产成人精品v| 久久6这里有精品| 国产高清不卡午夜福利| 能在线免费观看的黄片| 免费观看a级毛片全部| 日本三级黄在线观看| 国产精品,欧美在线| 久久这里有精品视频免费| 在线观看美女被高潮喷水网站| 国语自产精品视频在线第100页| 国产亚洲精品久久久com| 91久久精品电影网| 久久草成人影院| 日韩制服骚丝袜av| 中文字幕久久专区| 免费大片18禁| 日本成人三级电影网站| 亚洲人与动物交配视频| 久久精品国产自在天天线| 久久久a久久爽久久v久久| 久久久久久大精品| 国产白丝娇喘喷水9色精品| 毛片女人毛片| 国产精品福利在线免费观看| 尤物成人国产欧美一区二区三区| 色综合站精品国产| 黄色配什么色好看| 1024手机看黄色片| 成人三级黄色视频| 久久久a久久爽久久v久久| 一级毛片电影观看 | 人体艺术视频欧美日本| 美女 人体艺术 gogo| 色噜噜av男人的天堂激情| 亚洲av成人精品一区久久| 我要看日韩黄色一级片| 国产高清不卡午夜福利| 热99re8久久精品国产| 国产不卡一卡二| 亚洲人成网站高清观看| 一级毛片久久久久久久久女| 亚洲人成网站在线播放欧美日韩| 精品无人区乱码1区二区| 在线观看一区二区三区| 欧美最黄视频在线播放免费| 日日干狠狠操夜夜爽| 中文亚洲av片在线观看爽| 少妇人妻精品综合一区二区 | 国产视频内射| 美女国产视频在线观看| 日韩欧美在线乱码| 毛片一级片免费看久久久久| 悠悠久久av| 91精品国产九色| 亚洲av.av天堂| 久久久欧美国产精品| av卡一久久| av在线蜜桃| 天堂av国产一区二区熟女人妻| 69人妻影院| 国产极品精品免费视频能看的| 成年免费大片在线观看| 色吧在线观看| 免费av观看视频| 黄色视频,在线免费观看| 亚洲精华国产精华液的使用体验 | 99精品在免费线老司机午夜| 深夜精品福利| 中文在线观看免费www的网站| 国产精品av视频在线免费观看| 欧美另类亚洲清纯唯美| 日韩中字成人| 人妻久久中文字幕网| 午夜a级毛片| 国产精品永久免费网站| 久久99热这里只有精品18| 99久久成人亚洲精品观看| 久久精品国产清高在天天线| 又爽又黄a免费视频| 免费电影在线观看免费观看| 十八禁国产超污无遮挡网站| 中文精品一卡2卡3卡4更新| 亚洲在线观看片| 国产亚洲91精品色在线| 国产精品电影一区二区三区| 有码 亚洲区| 欧美高清成人免费视频www| 听说在线观看完整版免费高清| 国产伦在线观看视频一区| 亚洲七黄色美女视频| 国产精品.久久久| 国产成人影院久久av| 国产一级毛片七仙女欲春2| 国产片特级美女逼逼视频| 深夜a级毛片| 日韩视频在线欧美| 日韩av在线大香蕉| 91午夜精品亚洲一区二区三区| 别揉我奶头 嗯啊视频| 日韩强制内射视频| 99视频精品全部免费 在线| 久久精品国产亚洲网站| 国产成人91sexporn| 91精品一卡2卡3卡4卡| 一边摸一边抽搐一进一小说| 国产精品福利在线免费观看| 成人美女网站在线观看视频| 成人午夜高清在线视频| 日韩欧美 国产精品| 欧洲精品卡2卡3卡4卡5卡区| 噜噜噜噜噜久久久久久91| 熟女电影av网| 国产69精品久久久久777片| 丰满人妻一区二区三区视频av| 特级一级黄色大片| 亚洲国产色片| 日韩制服骚丝袜av| 色5月婷婷丁香| 一本一本综合久久| 日韩精品有码人妻一区| kizo精华| 99久国产av精品国产电影| 看非洲黑人一级黄片| 赤兔流量卡办理| 久久久久九九精品影院| av免费观看日本| a级一级毛片免费在线观看| 中国美女看黄片| 91aial.com中文字幕在线观看| 成人综合一区亚洲| 欧美成人a在线观看| 免费看a级黄色片| 真实男女啪啪啪动态图| 成人午夜高清在线视频| 男人舔奶头视频| 日本一本二区三区精品| 亚洲欧美日韩卡通动漫| 精品一区二区三区人妻视频| 亚洲精品乱码久久久v下载方式| 国产精品精品国产色婷婷| 亚洲国产高清在线一区二区三| 欧美成人免费av一区二区三区| 成人特级av手机在线观看| av.在线天堂| 午夜精品国产一区二区电影 | 老师上课跳d突然被开到最大视频| 少妇裸体淫交视频免费看高清| 91久久精品电影网| 尤物成人国产欧美一区二区三区| 免费观看人在逋| 久久热精品热| 久久久a久久爽久久v久久| 岛国在线免费视频观看| 久久久午夜欧美精品| 国产成人a区在线观看| 天天躁日日操中文字幕| 中文欧美无线码| 看片在线看免费视频|