• <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)化點源敏感性的鏡面視寧度評價
    手机成人av网站| 菩萨蛮人人尽说江南好唐韦庄| 亚洲,欧美精品.| 一级a爱视频在线免费观看| 婷婷丁香在线五月| 国产黄频视频在线观看| 蜜桃在线观看..| 黄色成人免费大全| 可以免费在线观看a视频的电影网站| 满18在线观看网站| 国产真人三级小视频在线观看| 欧美日韩av久久| 最新美女视频免费是黄的| 欧美日韩福利视频一区二区| 18禁美女被吸乳视频| 久久精品国产99精品国产亚洲性色 | 成人18禁在线播放| 电影成人av| 久久精品熟女亚洲av麻豆精品| 久久久久久免费高清国产稀缺| 夜夜骑夜夜射夜夜干| 久久ye,这里只有精品| 最新的欧美精品一区二区| 久久久水蜜桃国产精品网| 成人国产一区最新在线观看| 9热在线视频观看99| 日日摸夜夜添夜夜添小说| 一级片'在线观看视频| 国产欧美日韩一区二区精品| 色94色欧美一区二区| 久久国产精品男人的天堂亚洲| 精品国产超薄肉色丝袜足j| 精品国产乱子伦一区二区三区| 少妇的丰满在线观看| 视频在线观看一区二区三区| 国产在线视频一区二区| 99re在线观看精品视频| 精品国产超薄肉色丝袜足j| 国产精品98久久久久久宅男小说| 亚洲av片天天在线观看| 丁香六月欧美| 国产亚洲欧美精品永久| 久久九九热精品免费| 肉色欧美久久久久久久蜜桃| 午夜福利乱码中文字幕| 国产高清视频在线播放一区| 黄色成人免费大全| 国产在视频线精品| 国产高清videossex| 日韩一卡2卡3卡4卡2021年| 麻豆av在线久日| 久久中文字幕一级| 亚洲精品一二三| 免费看a级黄色片| 亚洲熟女毛片儿| 无限看片的www在线观看| 成人18禁高潮啪啪吃奶动态图| 久久久久视频综合| 老司机午夜十八禁免费视频| 91av网站免费观看| 亚洲精品美女久久久久99蜜臀| 精品亚洲成国产av| a级毛片黄视频| 女人爽到高潮嗷嗷叫在线视频| 人人澡人人妻人| 国产亚洲午夜精品一区二区久久| 操出白浆在线播放| 男女边摸边吃奶| 俄罗斯特黄特色一大片| 国产一区二区三区综合在线观看| 一进一出好大好爽视频| www.自偷自拍.com| 亚洲专区国产一区二区| 搡老熟女国产l中国老女人| 热re99久久国产66热| 国产aⅴ精品一区二区三区波| 久久影院123| 久久国产精品人妻蜜桃| 国产成人免费观看mmmm| 一区二区三区精品91| 十八禁人妻一区二区| 免费观看av网站的网址| 亚洲精品美女久久久久99蜜臀| 1024香蕉在线观看| 日本黄色日本黄色录像| 欧美黑人欧美精品刺激| 久久久精品免费免费高清| 波多野结衣av一区二区av| 每晚都被弄得嗷嗷叫到高潮| 黑人巨大精品欧美一区二区mp4| 看免费av毛片| 最近最新中文字幕大全电影3 | 咕卡用的链子| 国产精品电影一区二区三区 | 国产精品免费大片| 熟女少妇亚洲综合色aaa.| 国产欧美日韩综合在线一区二区| 欧美在线一区亚洲| 国产伦人伦偷精品视频| 丁香六月天网| 欧美+亚洲+日韩+国产| 欧美日韩亚洲高清精品| 欧美日韩黄片免| 妹子高潮喷水视频| 精品福利观看| 久热爱精品视频在线9| 欧美国产精品va在线观看不卡| 在线观看免费高清a一片| 黄色毛片三级朝国网站| 精品少妇内射三级| 国产av国产精品国产| 91老司机精品| 午夜久久久在线观看| 桃花免费在线播放| 日本一区二区免费在线视频| 久久性视频一级片| 日本wwww免费看| tocl精华| 久久精品91无色码中文字幕| 极品教师在线免费播放| 国产一区二区三区综合在线观看| 亚洲成av片中文字幕在线观看| 一进一出好大好爽视频| 91麻豆精品激情在线观看国产 | 国产成人欧美在线观看 | 久久香蕉激情| 男人操女人黄网站| 久久免费观看电影| 久久免费观看电影| 亚洲国产精品一区二区三区在线| 国产福利在线免费观看视频| 国产精品.久久久| 国产精品香港三级国产av潘金莲| 日本五十路高清| 日韩成人在线观看一区二区三区| 极品少妇高潮喷水抽搐| 大型黄色视频在线免费观看| 久久久久久人人人人人| 日韩免费av在线播放| 十八禁高潮呻吟视频| 精品一区二区三卡| 蜜桃在线观看..| 中文字幕人妻丝袜一区二区| 精品乱码久久久久久99久播| 亚洲av美国av| 国产欧美日韩精品亚洲av| 91成年电影在线观看| 久久精品人人爽人人爽视色| 丝袜美腿诱惑在线| 满18在线观看网站| 国产福利在线免费观看视频| 欧美日韩视频精品一区| 亚洲人成电影观看| 欧美日韩成人在线一区二区| 色精品久久人妻99蜜桃| 最新的欧美精品一区二区| 两性夫妻黄色片| 亚洲天堂av无毛| 99久久国产精品久久久| 免费在线观看黄色视频的| 50天的宝宝边吃奶边哭怎么回事| 黑人欧美特级aaaaaa片| 久久精品国产综合久久久| 下体分泌物呈黄色| 人人妻,人人澡人人爽秒播| 成年人免费黄色播放视频| 国产精品偷伦视频观看了| 伦理电影免费视频| 亚洲精品一卡2卡三卡4卡5卡| 色在线成人网| 久久久国产欧美日韩av| 成年动漫av网址| 他把我摸到了高潮在线观看 | 久久久国产一区二区| av一本久久久久| 在线观看免费高清a一片| 国产老妇伦熟女老妇高清| 手机成人av网站| 欧美成人午夜精品| 免费看a级黄色片| 99精品欧美一区二区三区四区| 最近最新免费中文字幕在线| 中亚洲国语对白在线视频| 又大又爽又粗| 日本一区二区免费在线视频| 成人国产av品久久久| 亚洲五月色婷婷综合| 亚洲 国产 在线| 新久久久久国产一级毛片| 国产在线观看jvid| 日韩欧美一区二区三区在线观看 | 国产亚洲av高清不卡| 日韩免费av在线播放| 国产一区有黄有色的免费视频| 日本欧美视频一区| 黄色视频,在线免费观看| 亚洲人成伊人成综合网2020| 窝窝影院91人妻| 麻豆乱淫一区二区| 欧美久久黑人一区二区| 国产亚洲欧美精品永久| 久久热在线av| 色综合婷婷激情| 欧美精品啪啪一区二区三区| 亚洲欧美激情在线| 老汉色∧v一级毛片| 久久精品国产a三级三级三级| 日韩免费av在线播放| 两个人免费观看高清视频| 日韩熟女老妇一区二区性免费视频| 纵有疾风起免费观看全集完整版| 亚洲天堂av无毛| 国产成人欧美在线观看 | 大片电影免费在线观看免费| 高潮久久久久久久久久久不卡| 国产在视频线精品| 亚洲精品国产区一区二| 日本欧美视频一区| 国产伦理片在线播放av一区| 色94色欧美一区二区| 亚洲av第一区精品v没综合| 欧美成狂野欧美在线观看| 国精品久久久久久国模美| 999精品在线视频| 黄色视频不卡| 国产色视频综合| av一本久久久久| 一边摸一边抽搐一进一出视频| tube8黄色片| 午夜福利在线免费观看网站| 精品国产乱码久久久久久小说| 国产成人啪精品午夜网站| 热99国产精品久久久久久7| 久久精品国产a三级三级三级| 亚洲成av片中文字幕在线观看| 视频在线观看一区二区三区| 一级片'在线观看视频| 国产成+人综合+亚洲专区| av线在线观看网站| 男女床上黄色一级片免费看| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲av一区麻豆| 99九九在线精品视频| 狠狠精品人妻久久久久久综合| 亚洲av美国av| 男女下面插进去视频免费观看| 最黄视频免费看| 亚洲精品国产色婷婷电影| 美女国产高潮福利片在线看| 丝袜在线中文字幕| 18禁国产床啪视频网站| 国产有黄有色有爽视频| 两性夫妻黄色片| 国产精品 欧美亚洲| 亚洲欧美一区二区三区久久| 午夜视频精品福利| 亚洲久久久国产精品| 天堂中文最新版在线下载| 亚洲一卡2卡3卡4卡5卡精品中文| 国产一区二区 视频在线| 亚洲一区中文字幕在线| 亚洲 国产 在线| 国产欧美日韩精品亚洲av| 久久免费观看电影| 首页视频小说图片口味搜索| 看免费av毛片| 丰满饥渴人妻一区二区三| 精品乱码久久久久久99久播| 这个男人来自地球电影免费观看| 国产精品久久电影中文字幕 | 国产区一区二久久| 亚洲国产欧美网| 波多野结衣av一区二区av| 国产男女超爽视频在线观看| 亚洲人成电影免费在线| 久久99热这里只频精品6学生| 美女视频免费永久观看网站| 免费在线观看完整版高清| 欧美精品一区二区免费开放| 9色porny在线观看| 少妇 在线观看| 日韩精品免费视频一区二区三区| 大型av网站在线播放| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕另类日韩欧美亚洲嫩草| av片东京热男人的天堂| 一级毛片精品| 亚洲精品中文字幕在线视频| 深夜精品福利| 极品人妻少妇av视频| 亚洲精品久久午夜乱码| 国产亚洲一区二区精品| 日韩视频在线欧美| 国产亚洲精品久久久久5区| 亚洲国产av新网站| 精品卡一卡二卡四卡免费| 亚洲专区中文字幕在线| 亚洲精品美女久久久久99蜜臀| 久久午夜亚洲精品久久| 免费av中文字幕在线| 国产老妇伦熟女老妇高清| 国产又色又爽无遮挡免费看| 黄频高清免费视频| 99久久精品国产亚洲精品| 欧美日韩国产mv在线观看视频| 午夜精品国产一区二区电影| 中文字幕人妻丝袜一区二区| 欧美人与性动交α欧美精品济南到| 欧美精品av麻豆av| 超碰97精品在线观看| 亚洲成av片中文字幕在线观看| 久久午夜亚洲精品久久| 9191精品国产免费久久| tube8黄色片| 国产亚洲精品久久久久5区| 乱人伦中国视频| 欧美黄色片欧美黄色片| 啪啪无遮挡十八禁网站| 日本撒尿小便嘘嘘汇集6| 黄色a级毛片大全视频| 在线观看66精品国产| 日韩 欧美 亚洲 中文字幕| 成在线人永久免费视频| 香蕉国产在线看| 久久久精品94久久精品| 极品少妇高潮喷水抽搐| 性高湖久久久久久久久免费观看| 久久精品亚洲av国产电影网| 伦理电影免费视频| 女人被躁到高潮嗷嗷叫费观| 天天操日日干夜夜撸| 久久人人爽av亚洲精品天堂| 日韩欧美国产一区二区入口| 欧美精品一区二区大全| 精品少妇久久久久久888优播| 91成人精品电影| 亚洲一区二区三区欧美精品| 国产成人精品久久二区二区免费| 国产精品久久久人人做人人爽| 国产国语露脸激情在线看| 亚洲欧美精品综合一区二区三区| 岛国毛片在线播放| 高清在线国产一区| 午夜免费成人在线视频| 日韩欧美免费精品| 国产精品一区二区免费欧美| 国产精品影院久久| 久久免费观看电影| xxxhd国产人妻xxx| 97在线人人人人妻| 交换朋友夫妻互换小说| 久久热在线av| 久久精品国产a三级三级三级| 国产精品 国内视频| 精品少妇久久久久久888优播| 91成人精品电影| 丝袜人妻中文字幕| 国产亚洲精品一区二区www | 日韩三级视频一区二区三区| 亚洲精品乱久久久久久| 激情视频va一区二区三区| h视频一区二区三区| 老鸭窝网址在线观看| 国产片内射在线| 亚洲精品中文字幕在线视频| 亚洲精品美女久久久久99蜜臀| 超色免费av| 一级毛片女人18水好多| 老司机亚洲免费影院| 国产深夜福利视频在线观看| 亚洲久久久国产精品| 成人手机av| 国产在线精品亚洲第一网站| 成人特级黄色片久久久久久久 | 黄频高清免费视频| 亚洲精品国产一区二区精华液| 久久精品亚洲熟妇少妇任你| 美女高潮喷水抽搐中文字幕| 狂野欧美激情性xxxx| 悠悠久久av| 亚洲欧美一区二区三区黑人| 亚洲精品av麻豆狂野| 美女午夜性视频免费| 亚洲人成电影免费在线| 日韩免费高清中文字幕av| 国产不卡av网站在线观看| 操美女的视频在线观看| 国产精品久久久久久精品电影小说| 日本黄色日本黄色录像| 久久久久国产一级毛片高清牌| 一边摸一边做爽爽视频免费| 成在线人永久免费视频| 国产日韩一区二区三区精品不卡| av天堂在线播放| 一边摸一边抽搐一进一小说 | 夜夜骑夜夜射夜夜干| 18禁国产床啪视频网站| aaaaa片日本免费| 天堂俺去俺来也www色官网| 亚洲av日韩在线播放| 一级毛片精品| 80岁老熟妇乱子伦牲交| 色综合欧美亚洲国产小说| 极品少妇高潮喷水抽搐| 亚洲熟女精品中文字幕| 亚洲自偷自拍图片 自拍| 亚洲,欧美精品.| 高清av免费在线| 伊人久久大香线蕉亚洲五| 超色免费av| 香蕉久久夜色| 在线天堂中文资源库| 亚洲人成伊人成综合网2020| 亚洲色图av天堂| 久久国产精品大桥未久av| 国产亚洲午夜精品一区二区久久| 色老头精品视频在线观看| 精品一品国产午夜福利视频| 18禁黄网站禁片午夜丰满| 免费在线观看日本一区| 女人高潮潮喷娇喘18禁视频| 在线观看免费日韩欧美大片| 久久精品亚洲熟妇少妇任你| 欧美激情极品国产一区二区三区| 国产精品一区二区在线不卡| 新久久久久国产一级毛片| 超色免费av| a级毛片黄视频| 国产一区二区三区在线臀色熟女 | 国产91精品成人一区二区三区 | 免费看a级黄色片| 一区二区三区精品91| www日本在线高清视频| 19禁男女啪啪无遮挡网站| 一本—道久久a久久精品蜜桃钙片| 国产视频一区二区在线看| 欧美激情高清一区二区三区| 国产片内射在线| 欧美av亚洲av综合av国产av| 一进一出抽搐动态| 少妇 在线观看| 免费女性裸体啪啪无遮挡网站| 国产成人av激情在线播放| 日韩视频在线欧美| 日本撒尿小便嘘嘘汇集6| 妹子高潮喷水视频| 99国产精品一区二区蜜桃av | 中文字幕av电影在线播放| 亚洲国产中文字幕在线视频| 日本av手机在线免费观看| 欧美成人午夜精品| 亚洲久久久国产精品| 可以免费在线观看a视频的电影网站| 91成人精品电影| 男女免费视频国产| 男女高潮啪啪啪动态图| 亚洲专区中文字幕在线| 国产黄频视频在线观看| 欧美激情久久久久久爽电影 | 国产欧美日韩一区二区三区在线| 国产精品av久久久久免费| 国产高清国产精品国产三级| 精品一品国产午夜福利视频| av片东京热男人的天堂| 久久久久久久大尺度免费视频| 亚洲av第一区精品v没综合| 亚洲情色 制服丝袜| 波多野结衣av一区二区av| 亚洲人成电影免费在线| 免费一级毛片在线播放高清视频 | 黄色a级毛片大全视频| 午夜福利视频在线观看免费| 亚洲七黄色美女视频| 午夜免费鲁丝| 美女午夜性视频免费| 午夜精品久久久久久毛片777| 色婷婷久久久亚洲欧美| 天堂动漫精品| 婷婷丁香在线五月| 丝袜在线中文字幕| 69av精品久久久久久 | 久久中文看片网| 99热网站在线观看| 在线亚洲精品国产二区图片欧美| 最近最新中文字幕大全免费视频| 国产成人av激情在线播放| 国产精品久久久av美女十八| 女警被强在线播放| 超碰97精品在线观看| 夜夜爽天天搞| 色综合婷婷激情| 欧美老熟妇乱子伦牲交| 老司机福利观看| videos熟女内射| 菩萨蛮人人尽说江南好唐韦庄| 国产精品影院久久| aaaaa片日本免费| 男女午夜视频在线观看| 热99re8久久精品国产| 老汉色∧v一级毛片| 欧美日韩视频精品一区| 黑人巨大精品欧美一区二区mp4| 夜夜爽天天搞| 色老头精品视频在线观看| 天天添夜夜摸| 99久久99久久久精品蜜桃| 女人久久www免费人成看片| 久久国产亚洲av麻豆专区| av有码第一页| 亚洲av片天天在线观看| 国产精品av久久久久免费| 午夜福利视频精品| 黑丝袜美女国产一区| 人人妻人人添人人爽欧美一区卜| 99re在线观看精品视频| 国产成人免费观看mmmm| 日韩免费av在线播放| 久久九九热精品免费| 欧美成人午夜精品| 亚洲精品自拍成人| 一边摸一边做爽爽视频免费| 免费人妻精品一区二区三区视频| 在线 av 中文字幕| 啪啪无遮挡十八禁网站| 97人妻天天添夜夜摸| 18禁美女被吸乳视频| 五月天丁香电影| 久久99一区二区三区| 757午夜福利合集在线观看| 最近最新中文字幕大全免费视频| 高清欧美精品videossex| 大陆偷拍与自拍| netflix在线观看网站| 91成年电影在线观看| 亚洲中文日韩欧美视频| 一个人免费在线观看的高清视频| 中文字幕高清在线视频| 国产欧美亚洲国产| 亚洲 欧美一区二区三区| 黄色视频不卡| 精品少妇久久久久久888优播| 国产在线视频一区二区| 制服人妻中文乱码| 另类亚洲欧美激情| 男女高潮啪啪啪动态图| 热99久久久久精品小说推荐| 多毛熟女@视频| 国产精品久久久人人做人人爽| 亚洲一区中文字幕在线| 亚洲精品美女久久av网站| 亚洲欧美日韩高清在线视频 | 国产成人欧美在线观看 | 欧美av亚洲av综合av国产av| 精品久久久久久电影网| 91字幕亚洲| 99精品久久久久人妻精品| 亚洲精品在线美女| 18在线观看网站| 国产精品免费视频内射| 69av精品久久久久久 | 性高湖久久久久久久久免费观看| 50天的宝宝边吃奶边哭怎么回事| 久久毛片免费看一区二区三区| 高清av免费在线| 欧美在线黄色| 老司机午夜福利在线观看视频 | 亚洲熟妇熟女久久| 高潮久久久久久久久久久不卡| 欧美老熟妇乱子伦牲交| 一级毛片电影观看| 国产主播在线观看一区二区| 色综合欧美亚洲国产小说| 国产精品麻豆人妻色哟哟久久| 国产成人免费观看mmmm| 午夜免费成人在线视频| 菩萨蛮人人尽说江南好唐韦庄| 好男人电影高清在线观看| 99热国产这里只有精品6| 女人被躁到高潮嗷嗷叫费观| 久久午夜亚洲精品久久| 精品国产超薄肉色丝袜足j| 香蕉国产在线看| 欧美日韩黄片免| 亚洲天堂av无毛| 亚洲专区国产一区二区| 久久精品成人免费网站| a级毛片黄视频| 久久久水蜜桃国产精品网| 多毛熟女@视频| 国产伦理片在线播放av一区| 中文字幕色久视频| 韩国精品一区二区三区| 性少妇av在线| 欧美精品高潮呻吟av久久| www.精华液| 黄片小视频在线播放| 久久青草综合色| 久久国产精品大桥未久av| 成年人午夜在线观看视频| 啪啪无遮挡十八禁网站| 亚洲欧美激情在线| 中文欧美无线码| 久久久久久久大尺度免费视频| 国产国语露脸激情在线看| 18禁观看日本| 亚洲人成电影观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美一区二区三区黑人| 天堂俺去俺来也www色官网| 99久久国产精品久久久| 亚洲精品久久成人aⅴ小说| 黄色a级毛片大全视频| 久久人妻福利社区极品人妻图片| 久久国产精品影院| 免费观看a级毛片全部|