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

    基于ANSYS的CSAMT法三維正演仿真模擬基礎研究

    2016-06-14 01:09鄧偉王少華張露
    電腦知識與技術 2016年12期

    鄧偉+王少華+張露

    摘要:介紹了基于ANSYS軟件的CSAMT法三維正演的建模及網(wǎng)格劃分方法,研究模型參數(shù)設置與卡尼亞視電阻率之間的關系,提出了CSAMT三維正演時模型參數(shù)設置建議及規(guī)律。研究結果表明:合理設置過渡區(qū)可以提高網(wǎng)格劃分精度,非均勻劃分網(wǎng)格可以有效減少網(wǎng)格數(shù)量,提高計算速度;線圈高度的變化會影響磁場垂直入射效果;線圈高度一定時,模型空氣層厚度會對視電阻率值產(chǎn)生影響;模型地層厚度會對視電阻率值產(chǎn)生很大影響,測量有效深度與模型地層厚度間存在一定比例關系,從而保證視電阻率值的準確性和規(guī)律性。

    關鍵詞:CSAMT法;三維正演;ANSYS仿真;臨界比例;視電阻率

    中圖分類號:TP391 文獻標識碼:A 文章編號:1009-3044(2016)12-0228-04

    Abstract: This paper makes an introduction about modeling and meshing of 3D CSAMT forward modeling by ANSYS, researching the relation between model parameters and Cagniard resistivity, coming up with a reference about setting of 3D CSAMT modeling parameters. The results indicate: meshing accuracy is improved by use of transition regionm, inhomogeneous meshing can effictively decrease quantity and enhance simulation rate. Coil height can influence the incoming of magnetic field which should be vertical. When coil height is certain, air thickness of the model influence apparent resistivity. Depth of the ground has a strong influence on apparent resistivity, and there is a proportion relation between surveying depth and ground depth. An appropriate ground depth can keep the simulated result accurate and regular.

    Key words:CSAMT; 3D forward modeling; Ansys simulation; critical proportion; apparent resistivity

    1 概述

    CSAMT法是一種易于操作且成本較低的地質(zhì)探測手段,它采用有限長線電流源作為激發(fā)源,在一定距離外布設測線接收電場和磁場,從而獲得卡尼亞視電阻率,并據(jù)此判斷地表下各地層的地質(zhì)構造[1]。目前,一維、二維以及2.5維CSAMT法正、反演研究已經(jīng)較為全面[2-5],理論、模擬都很成熟。但是從實際問題考慮,特別是對于復雜局部地質(zhì)體的勘探來說,三維CSAMT方法的正、反演才能更客觀、準確地反映地質(zhì)特征[6-9]。

    通過理論方法求解三維正反演問題是相當困難,而有限元仿真在解決場問題上有著極大優(yōu)勢,采用有限元仿真軟件ANSYS對三維地質(zhì)模型進行正演模擬[10-11],研究模型建立、網(wǎng)格劃分、參數(shù)設置、仿真計算、數(shù)據(jù)提取的方法,為三維CSAMT法正演提供一種參考。

    2 ANSYS仿真模擬原理

    2.1 CSAMT法視電阻率計算

    3 仿真模型建立及視電阻率計算

    考慮到工程實際情況,假設模型需要實現(xiàn)探測100m以內(nèi)的地層視電阻率,即模擬時的有效穿透深度[≥100m]。根據(jù)收發(fā)距要求,[r=(3~5)h有效≥300~500m],選擇收發(fā)距[r=500m],模型長、寬為[1km×1.3km];線圈內(nèi)、外徑[8m×12m];初始空氣層厚度設為300m,地層設為200m,空氣層幾乎不導電,電阻率設為[106Ω·m],地層電阻率[ρ0=100Ω·m],線圈施加環(huán)向電流,電流密度[Js=10Am2]。

    3.1 線圈高度變化與磁場入射角

    在線圈內(nèi)、外徑一定時,線圈高度的變化會影響磁場線入射地面時的角度,而卡尼亞視電阻率的計算要求是相互垂直的磁場和電場分量,所以在測量區(qū)應該讓磁場線盡量垂直入射,以保證計算的準確性。研究線圈磁場時采用二維模型,改變線圈高度,查看測量區(qū)[r≥500m]附近的磁場線入射角度。模擬結果如圖2所示。

    通過改變線圈高度可以看到,隨著線圈高度減小,磁場線垂直入射效果越好,當線圈內(nèi)、外徑一定時,高度決定線圈截面積,但場強變化并不會改變磁場方向。進一步研究發(fā)現(xiàn),當線圈高度為200m、100m、60m、30m時,入射角度近似為62°、85°、90°、90°,線圈高度下降到一定值以后就能具有很好的垂直入射效果。后續(xù)仿真模擬選取線圈高度為30m。

    3.2 網(wǎng)格

    三維模型示意如圖1所示,整個模型外形為立方體,但線圈處于立方體內(nèi)部,劃分網(wǎng)格時必須進行有效切分且進行有效的網(wǎng)格過渡處理,以保證線圈處以及測量區(qū)域網(wǎng)格細化。線圈與邊界、測量區(qū)與邊界之間需有一定距離的過渡區(qū),因為模型并不是無限大,所以模型邊界處數(shù)據(jù)必定存在誤差。三維模型網(wǎng)格劃分結果如圖3所示,在線圈網(wǎng)格外部設置了過渡網(wǎng)格,測量區(qū)與邊界之間網(wǎng)格由細到粗,自然過渡,空氣層網(wǎng)格等比例劃分。

    對于地層網(wǎng)格,有效深度以內(nèi)部分采用等距離細網(wǎng)格劃分,有效深度以外等比例網(wǎng)格劃分。這樣處理的好處是既保證了網(wǎng)格劃分精度,又可以有效地控制網(wǎng)格數(shù)量,提高運算速度。

    3.3 不同頻率下視電阻率值

    根據(jù)有效深度計算公式(3),不同激發(fā)頻率的線圈產(chǎn)生的磁場能夠穿透的地層深度不同。模型地層厚度為200m,由公式可推算出穿透10m~200m范圍大致的線圈激發(fā)頻率,首先在200~100000Hz激發(fā)頻率范圍內(nèi)運行仿真模型,提取測量帶中心處的電場、磁場分量,其值如表1所示。再利用公式計算出視電阻率,不同頻率下視電阻率計算值如圖4所示。

    從圖中可以看出,在高頻區(qū),視電阻率值與設定的電阻率值保持了較好地吻合,但是當場源頻率下降到某一值(1000Hz附近)之后,視電阻率值快速增大,偏離真實值。這個結果是與采用有限長導線作為激發(fā)源的模擬效果相符的[12]。

    3.4 空氣層厚度對視電阻率值影響

    空氣層雖然認為近似不導電,但作為磁場傳播介質(zhì),其厚度對穿過空氣層進入地層的磁場有影響??刂破渌麉?shù),僅改變空氣層厚度,計算有效深度為20m、30m、40m時的視電阻率值。不同空氣層厚度是測量區(qū)中心點數(shù)據(jù)對應視電阻率值如圖5所示。

    從圖中可以看出,在200m空氣層厚度時整體偏離地層電阻率值較大,300m~600m范圍內(nèi)視電阻率值趨于穩(wěn)定。這說明空氣層厚度確實會對模擬結果產(chǎn)生影響,可以推測在其他參數(shù)一定時,只有當模型空氣層厚度超過某一臨界值之后,仿真結果才趨于穩(wěn)定。具體是與線圈高度還是地層深度有關有待進一步探究。

    3.5 地層深度對視電阻率值影響

    保持模型既定參數(shù)不變,空氣層厚度設為300m,模型有效穿透深度為100m,根據(jù)有效深度計算公式,在10m~100m范圍內(nèi)按照每10m一個測深間隔的順序,計算出對應諧波電流頻率。在ANSYS中通過建立模型多載荷步分析,可一次完成10個頻段計算,地層深度200m時視電阻率計算值如表2所示,表中[σ]表示視電阻率值與真實電阻率的偏差。

    由表中數(shù)據(jù)可以看出,此時視電阻率的模擬效果并不理想,視電阻率值與真實值偏差范圍達到了-26%~66%,測深100m時的視電阻率值更是遠遠偏離實際值。通過查看淺層地層和深層地層的感應電場線可以發(fā)現(xiàn),淺層時的感應電場線是較為規(guī)則的,而一定深度后電場線將逐漸發(fā)生變形,越靠近邊界這種現(xiàn)象越明顯,這必將影響視電阻率計算值的準確性。

    當探測深度為一定值時,應該對模型的地層深度適當加深,以保證提取數(shù)據(jù)的深度范圍電場線變形較小。因此,改變模型地層深度,查看在不同深度下10m~100m范圍內(nèi)視電阻率值如何變化。

    在300m~500m范圍內(nèi)改變地層深度,計算200m、300m、400m、500m深度下各頻段的視電阻值,結果如圖6所示。

    可以看到,隨著深度的增加,視電阻率值計算值逐漸趨于穩(wěn)定,偏離程度減小。在400m時,10m~100m有效深度范圍內(nèi),如表3所示,視電阻率值偏離真實電阻率值范圍在-5%~15%,這已經(jīng)能夠很好的反映地層電阻率特征。

    同時,由圖6還可以發(fā)現(xiàn),測深范圍在10m~40m范圍內(nèi)時,地層深度200m~500m范圍視電阻率值幾乎相同,即測深[≤40m]時,地層深度設為200m即可確保計算準確性;而當有效深度為50m~60m時,只有在300m~500m范圍視電阻率值接近,即地層深度[Hd≤60m]時,地層深度設為300m可保證計算準確性;最后,當有效深度范圍為70m~100m時,地層深度取為400m~500m均可得到較為理想的視電阻率值。由上述分析可以發(fā)現(xiàn),為了確保模型計算出的視電阻率值準確性,有效測量深度和地層深度之間存在著一個比例關系,在本文研究的參數(shù)下,比例系數(shù)大致為1:5,即地層深度為有效深度5倍時,可保證較好仿真效果。

    地層深度為400m時,誤差上下限已經(jīng)比較理想,此時,再提取測量帶中心截面上以測量帶中心為原點,-80m~80m范圍內(nèi)9個點的磁場、電場值,計算截面視電阻率。結果如圖7所示。

    通過截面視電阻率圖可以看出,在-80m~80m范圍內(nèi),視電阻值保持了很好的準確性。說明采用上述的模型處理方法進行三維正演的仿真,可以得到較好模擬效果。

    4 結論

    通過建立三維有限元仿真模型并對模型進行有效地網(wǎng)格處理、參數(shù)設置,得出了進行CSAMT三維正演模擬一些有參考價值的結論,具體如下:

    ①采用線圈作為激發(fā)源進行模擬仿真,有別于傳統(tǒng)的偶極源激發(fā)源,卡尼亞視電阻率值能夠很好反映測量區(qū)電阻率。線圈高度會影響磁場入射地面的角度,劃分模型網(wǎng)格時,在線圈、測量區(qū)處要進行網(wǎng)格細化,模型各部分之間采用過渡網(wǎng)格、等比例劃分網(wǎng)格可以保證運算精度并提高運算速度。

    ②空氣層厚度對仿真結果有影響,當模型要求測深一定時,空氣層厚度需要超過一臨界值才能保證視電阻率計算值的準確性。

    ③地層厚度直接影響視電阻率計算值,只有合理調(diào)整縱向地層厚度才能保證一定范圍內(nèi)視電阻率值的準確性,總的規(guī)律是隨著電磁波測量有效深度的增加,地層厚度需要增加。且在探測范圍[≤]100m時,有效深度與地層厚度之間存在近似1:5的臨界比例關系,模型地層厚度設為有效探測深度的5倍可以保證很好的仿真效果。

    文中僅探究了以線圈為激發(fā)源的三維CSAMT法正演的有限元仿真實現(xiàn)方法以及部分參數(shù)的變化規(guī)律,主要是模型高度方向變化規(guī)律,對CSAMT三維正演研究提供一種參考,為下一步引入吸收邊界的視電阻率研究提供基礎支撐。模型長度與寬度對仿真應該也有一定影響,有待進一步研究并總結出更全面的規(guī)律。

    參考文獻:

    [1] 底青云,王若.可控源音頻大地電磁數(shù)據(jù)正反演及方法應用[M].北京:科學出版社,2008.

    [2] 王若,王妙月.一維全資料CSAMT反演[J].石油地球物理勘探,2007,42(1):107-110.

    [3] 龔強.CSAMT二維正演研究[D].北京:中國地質(zhì)大學,2007.

    [4] 鄭仁淑,鐘宏偉.地表上方電偶極源響應特征研究[J].工程地球物理學報,2005,2(5):349-350.

    [5] 許廣春.CSAMT法在實際應用中的若干問題[J].物探與化探,2011,35(6):810-812.

    [6] 王若.CSAMT三維單分量有限元正演[J].地球物理學進展,2014,29(2):840-843.

    [7] 陳銳.CSAMT三維交錯采樣有限元差分數(shù)值模擬并行算法研究[D].北京:中國地質(zhì)大學,2012.

    [8] 黃臨平,戴世坤. 復雜條件下三維電磁場有限元計算法[J].地球科學—中國地質(zhì)大學學報,2002,27(6):775-780.

    [9] 閻述,陳明生. 電偶極源頻率電磁測深三維地電模型有限元正演[J]. 煤田地質(zhì)與勘探,2000(3):50-56.

    [10] 張倩.ANSYS 12.0電磁學有限元分析從入門到精通[M].北京:機械工業(yè)出版社,2010.

    [11] 謝龍輝,耿煜,邱婉. ANSYS電磁場分析[M]. 北京:電子工業(yè)出版社,2012.

    [12] 王若,王妙云,盧元林. 三維三分量CSAMT法有限元正演模擬研究初探[J]. 地球物理學進展,2007,22(2):584.

    亚洲精品在线美女| 精品欧美国产一区二区三| 久久精品综合一区二区三区| 男人舔奶头视频| 欧美日本亚洲视频在线播放| 国产激情欧美一区二区| 特大巨黑吊av在线直播| 黄频高清免费视频| 亚洲国产欧美网| 人妻久久中文字幕网| 88av欧美| 一进一出抽搐gif免费好疼| 中文资源天堂在线| 蜜桃久久精品国产亚洲av| 91久久精品国产一区二区成人 | 国产精品影院久久| 99re在线观看精品视频| 国产av不卡久久| 午夜亚洲福利在线播放| 欧美乱妇无乱码| 美女 人体艺术 gogo| www国产在线视频色| 成人18禁在线播放| 免费观看精品视频网站| 超碰成人久久| 麻豆久久精品国产亚洲av| 国产三级在线视频| 99精品欧美一区二区三区四区| АⅤ资源中文在线天堂| 波多野结衣高清作品| 国产精品av久久久久免费| 亚洲黑人精品在线| 日本撒尿小便嘘嘘汇集6| 久久草成人影院| 级片在线观看| 亚洲无线观看免费| 黄色 视频免费看| 女人高潮潮喷娇喘18禁视频| 丰满人妻一区二区三区视频av | 国产亚洲精品一区二区www| 国产午夜福利久久久久久| 两个人的视频大全免费| 欧美日韩黄片免| 最新美女视频免费是黄的| 久久精品亚洲精品国产色婷小说| www国产在线视频色| 日韩三级视频一区二区三区| 97人妻精品一区二区三区麻豆| 色哟哟哟哟哟哟| 日韩国内少妇激情av| 99国产精品99久久久久| 少妇裸体淫交视频免费看高清| 一个人免费在线观看电影 | 中文字幕高清在线视频| 麻豆国产av国片精品| 亚洲五月天丁香| 日本 av在线| 欧美日韩国产亚洲二区| 国产单亲对白刺激| 欧美乱妇无乱码| 噜噜噜噜噜久久久久久91| 国产黄片美女视频| 少妇的逼水好多| 中文资源天堂在线| 亚洲 欧美一区二区三区| 国产在线精品亚洲第一网站| 最新美女视频免费是黄的| 免费在线观看视频国产中文字幕亚洲| 欧美大码av| 精品国内亚洲2022精品成人| 国产精品99久久久久久久久| 国产精品综合久久久久久久免费| 欧美午夜高清在线| 久久久久久大精品| 久久久久九九精品影院| 国产精品久久视频播放| 中国美女看黄片| 成人性生交大片免费视频hd| 午夜成年电影在线免费观看| 高潮久久久久久久久久久不卡| 91九色精品人成在线观看| 日本黄大片高清| 一a级毛片在线观看| 欧美性猛交╳xxx乱大交人| 一夜夜www| 成在线人永久免费视频| 搡老妇女老女人老熟妇| 毛片女人毛片| 熟女少妇亚洲综合色aaa.| 中文资源天堂在线| 我的老师免费观看完整版| 国产伦精品一区二区三区视频9 | 他把我摸到了高潮在线观看| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 亚洲av电影不卡..在线观看| 99久久99久久久精品蜜桃| 超碰成人久久| 亚洲电影在线观看av| 亚洲乱码一区二区免费版| 宅男免费午夜| 国产成人精品无人区| 最新在线观看一区二区三区| 一级毛片精品| 中文亚洲av片在线观看爽| 久久亚洲精品不卡| 日本免费a在线| 亚洲成人久久爱视频| 亚洲欧洲精品一区二区精品久久久| 亚洲专区字幕在线| 97碰自拍视频| 久久久水蜜桃国产精品网| 高清毛片免费观看视频网站| 成熟少妇高潮喷水视频| 亚洲午夜理论影院| 一进一出抽搐gif免费好疼| 久久精品影院6| 岛国在线免费视频观看| 在线看三级毛片| 国产成人欧美在线观看| 国产高清三级在线| 日日夜夜操网爽| 黑人操中国人逼视频| 国产乱人视频| 亚洲av免费在线观看| 十八禁人妻一区二区| 欧美日韩亚洲国产一区二区在线观看| 一本综合久久免费| av天堂在线播放| 亚洲国产欧美人成| 久久九九热精品免费| 亚洲国产精品久久男人天堂| 欧美不卡视频在线免费观看| 一级a爱片免费观看的视频| 久久精品综合一区二区三区| 好男人在线观看高清免费视频| 欧美大码av| 熟妇人妻久久中文字幕3abv| 日韩人妻高清精品专区| 亚洲18禁久久av| 禁无遮挡网站| 久久国产精品人妻蜜桃| 网址你懂的国产日韩在线| 精品日产1卡2卡| 中文字幕最新亚洲高清| 日本熟妇午夜| 国产av在哪里看| or卡值多少钱| 女警被强在线播放| 欧美最黄视频在线播放免费| 国产亚洲精品久久久久久毛片| avwww免费| а√天堂www在线а√下载| 国产成人欧美在线观看| 欧美绝顶高潮抽搐喷水| 欧美+亚洲+日韩+国产| 久久久国产精品麻豆| 99在线视频只有这里精品首页| 亚洲天堂国产精品一区在线| 天天一区二区日本电影三级| 美女大奶头视频| 国产精品精品国产色婷婷| 国产日本99.免费观看| 久久久久久久久中文| 大型黄色视频在线免费观看| 亚洲18禁久久av| 操出白浆在线播放| 成人欧美大片| 国产爱豆传媒在线观看| 长腿黑丝高跟| 亚洲人成电影免费在线| 国产伦精品一区二区三区四那| 日韩欧美 国产精品| 亚洲精品国产精品久久久不卡| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区免费欧美| 精品国产超薄肉色丝袜足j| 国产精品免费一区二区三区在线| 一边摸一边抽搐一进一小说| 欧美丝袜亚洲另类 | 婷婷六月久久综合丁香| 男女视频在线观看网站免费| 婷婷丁香在线五月| 色老头精品视频在线观看| 精品国产三级普通话版| 欧美zozozo另类| 日韩欧美在线二视频| 国语自产精品视频在线第100页| 日本 av在线| 丝袜人妻中文字幕| 国产aⅴ精品一区二区三区波| 国产亚洲精品久久久久久毛片| 亚洲精品色激情综合| 18禁裸乳无遮挡免费网站照片| 久久午夜综合久久蜜桃| 国产亚洲精品久久久久久毛片| 欧美午夜高清在线| 国产精华一区二区三区| 高潮久久久久久久久久久不卡| 午夜视频精品福利| 综合色av麻豆| 国产亚洲欧美98| 国产成+人综合+亚洲专区| 国产成人系列免费观看| 欧美zozozo另类| 99久久成人亚洲精品观看| 国产乱人伦免费视频| 一a级毛片在线观看| 性色avwww在线观看| 国产精品 欧美亚洲| 网址你懂的国产日韩在线| 国产私拍福利视频在线观看| 成人无遮挡网站| 一本精品99久久精品77| 熟女电影av网| 哪里可以看免费的av片| 免费观看人在逋| 成年免费大片在线观看| 一个人看视频在线观看www免费 | 午夜两性在线视频| 色吧在线观看| 国产97色在线日韩免费| 国产精品亚洲一级av第二区| 欧美性猛交黑人性爽| 亚洲男人的天堂狠狠| 悠悠久久av| 狠狠狠狠99中文字幕| 国产高清videossex| 成人国产一区最新在线观看| 国产精华一区二区三区| 不卡av一区二区三区| 长腿黑丝高跟| 国产综合懂色| 精品国产三级普通话版| 少妇的丰满在线观看| 午夜福利在线观看免费完整高清在 | 欧美3d第一页| ponron亚洲| 亚洲国产精品成人综合色| 成人欧美大片| 久久婷婷人人爽人人干人人爱| 国产麻豆成人av免费视频| 手机成人av网站| 国产一区二区在线av高清观看| 成人亚洲精品av一区二区| 麻豆久久精品国产亚洲av| 日韩av在线大香蕉| 亚洲av成人一区二区三| 久久久久久国产a免费观看| 国产av一区在线观看免费| 99国产极品粉嫩在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人一区二区三| 99在线人妻在线中文字幕| 成人鲁丝片一二三区免费| 亚洲成人久久性| 亚洲精华国产精华精| 国产午夜精品论理片| 国产成年人精品一区二区| 欧美一级毛片孕妇| 人人妻,人人澡人人爽秒播| 俺也久久电影网| 成人高潮视频无遮挡免费网站| 国产成人福利小说| 99热这里只有精品一区 | 高清毛片免费观看视频网站| 黄色日韩在线| 日韩av在线大香蕉| 91九色精品人成在线观看| 一二三四社区在线视频社区8| 国产激情久久老熟女| 日韩欧美国产一区二区入口| 国产精品美女特级片免费视频播放器 | 嫩草影院精品99| 高潮久久久久久久久久久不卡| 久久精品aⅴ一区二区三区四区| 免费在线观看亚洲国产| 午夜福利高清视频| 精品欧美国产一区二区三| www.999成人在线观看| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看| 一进一出抽搐动态| 久久久久久九九精品二区国产| 亚洲精华国产精华精| 欧美绝顶高潮抽搐喷水| 免费观看人在逋| 小蜜桃在线观看免费完整版高清| 18禁国产床啪视频网站| 国内久久婷婷六月综合欲色啪| 国产野战对白在线观看| av片东京热男人的天堂| 久久久久久人人人人人| 男人和女人高潮做爰伦理| 白带黄色成豆腐渣| 久久久久久久精品吃奶| 色老头精品视频在线观看| 国产精品av久久久久免费| 免费在线观看日本一区| 亚洲av免费在线观看| 午夜成年电影在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 国产69精品久久久久777片 | 欧美日韩亚洲国产一区二区在线观看| 五月伊人婷婷丁香| 美女高潮的动态| 国产三级中文精品| 久久久国产成人免费| 一级黄色大片毛片| 这个男人来自地球电影免费观看| 一级黄色大片毛片| 久久久久国内视频| 国产亚洲精品久久久久久毛片| 啪啪无遮挡十八禁网站| 欧美激情久久久久久爽电影| cao死你这个sao货| 手机成人av网站| 亚洲精品中文字幕一二三四区| 中文字幕av在线有码专区| 国产成年人精品一区二区| 久久久久精品国产欧美久久久| 国产成人系列免费观看| 欧美中文日本在线观看视频| 最好的美女福利视频网| 在线看三级毛片| 91字幕亚洲| 高清毛片免费观看视频网站| 久久精品91蜜桃| 亚洲精品久久国产高清桃花| 免费av不卡在线播放| cao死你这个sao货| 国语自产精品视频在线第100页| 亚洲七黄色美女视频| 久久精品国产综合久久久| 床上黄色一级片| 国产 一区 欧美 日韩| 亚洲片人在线观看| 亚洲国产中文字幕在线视频| 最新中文字幕久久久久 | 久久久久亚洲av毛片大全| 久9热在线精品视频| 午夜福利在线观看免费完整高清在 | 欧美午夜高清在线| 最新美女视频免费是黄的| 国产精品久久久人人做人人爽| 成熟少妇高潮喷水视频| 亚洲美女视频黄频| 国产高清videossex| 久久久成人免费电影| 校园春色视频在线观看| 国产真人三级小视频在线观看| 夜夜夜夜夜久久久久| 欧美+亚洲+日韩+国产| 18禁黄网站禁片午夜丰满| 黄色 视频免费看| 后天国语完整版免费观看| 中文亚洲av片在线观看爽| 成人无遮挡网站| 在线a可以看的网站| 99久久无色码亚洲精品果冻| 国产乱人伦免费视频| 久久亚洲真实| 黄色成人免费大全| 久久久久国内视频| 亚洲中文字幕一区二区三区有码在线看 | 看黄色毛片网站| 亚洲精品国产精品久久久不卡| 国产精品电影一区二区三区| 90打野战视频偷拍视频| 一本综合久久免费| 成人午夜高清在线视频| 亚洲av片天天在线观看| 亚洲片人在线观看| 岛国在线免费视频观看| 国产主播在线观看一区二区| 国产成人aa在线观看| 少妇熟女aⅴ在线视频| www.999成人在线观看| 久久伊人香网站| 国产又色又爽无遮挡免费看| 成人性生交大片免费视频hd| 国产精品久久久久久久电影 | 一二三四社区在线视频社区8| 国产激情久久老熟女| 国产精品久久久久久亚洲av鲁大| 免费观看的影片在线观看| 亚洲最大成人中文| 午夜福利免费观看在线| 天堂网av新在线| av国产免费在线观看| 免费观看的影片在线观看| 好看av亚洲va欧美ⅴa在| 亚洲av日韩精品久久久久久密| 午夜精品一区二区三区免费看| 一a级毛片在线观看| 成年免费大片在线观看| 欧美三级亚洲精品| 制服丝袜大香蕉在线| 亚洲在线自拍视频| 波多野结衣高清无吗| 日本免费a在线| 亚洲欧美日韩高清在线视频| 人人妻人人澡欧美一区二区| 黄频高清免费视频| 精品一区二区三区av网在线观看| 999精品在线视频| 两个人视频免费观看高清| 级片在线观看| 久久久久久九九精品二区国产| 国产成人精品久久二区二区91| 精品久久久久久久人妻蜜臀av| 久久婷婷人人爽人人干人人爱| 久9热在线精品视频| 免费看光身美女| 九色成人免费人妻av| 伦理电影免费视频| 久久精品影院6| 一级毛片高清免费大全| av在线天堂中文字幕| 黑人操中国人逼视频| 一夜夜www| 无限看片的www在线观看| 欧美黑人欧美精品刺激| 久久久久精品国产欧美久久久| 美女高潮喷水抽搐中文字幕| 亚洲在线自拍视频| 免费电影在线观看免费观看| 国产精品永久免费网站| 免费高清视频大片| 亚洲av片天天在线观看| 岛国在线免费视频观看| 亚洲熟女毛片儿| 久久九九热精品免费| 黄色 视频免费看| 美女 人体艺术 gogo| 最好的美女福利视频网| 欧美中文综合在线视频| 亚洲国产欧洲综合997久久,| 久久精品国产综合久久久| 亚洲欧美一区二区三区黑人| 九九在线视频观看精品| 国产成人aa在线观看| 亚洲美女黄片视频| 高清在线国产一区| 久久久水蜜桃国产精品网| 亚洲午夜理论影院| 免费搜索国产男女视频| 国产一区二区三区视频了| 三级毛片av免费| 亚洲国产精品999在线| svipshipincom国产片| 国产成年人精品一区二区| 波多野结衣巨乳人妻| 国产精品av久久久久免费| 欧美不卡视频在线免费观看| tocl精华| 色在线成人网| 天堂影院成人在线观看| 久久久久久国产a免费观看| 国产蜜桃级精品一区二区三区| 国产av一区在线观看免费| 黄色成人免费大全| 别揉我奶头~嗯~啊~动态视频| 免费看a级黄色片| 99热这里只有精品一区 | 国产单亲对白刺激| 久久婷婷人人爽人人干人人爱| 天堂网av新在线| 亚洲狠狠婷婷综合久久图片| 网址你懂的国产日韩在线| 午夜精品在线福利| 午夜免费成人在线视频| 国产久久久一区二区三区| 欧美黄色片欧美黄色片| 91在线观看av| 久久精品91无色码中文字幕| 国产极品精品免费视频能看的| 精品一区二区三区四区五区乱码| avwww免费| 精品国产乱子伦一区二区三区| 久久久久久大精品| 欧美日韩精品网址| 后天国语完整版免费观看| 男女午夜视频在线观看| 色尼玛亚洲综合影院| 亚洲熟妇熟女久久| 国产成+人综合+亚洲专区| 18美女黄网站色大片免费观看| 亚洲电影在线观看av| 午夜激情欧美在线| 国产不卡一卡二| 九色国产91popny在线| 精品久久久久久久毛片微露脸| 少妇熟女aⅴ在线视频| 禁无遮挡网站| 午夜福利18| 日韩欧美国产在线观看| 日韩三级视频一区二区三区| 91九色精品人成在线观看| 欧美黑人巨大hd| 色哟哟哟哟哟哟| 女人被狂操c到高潮| 最新中文字幕久久久久 | 亚洲国产精品合色在线| 免费人成视频x8x8入口观看| 成人性生交大片免费视频hd| 免费一级毛片在线播放高清视频| 宅男免费午夜| 久久人妻av系列| 在线免费观看的www视频| 精品午夜福利视频在线观看一区| 久久这里只有精品中国| 最近最新中文字幕大全电影3| 久久久久久久久久黄片| 少妇的逼水好多| 亚洲欧美一区二区三区黑人| 亚洲av五月六月丁香网| 18美女黄网站色大片免费观看| 一级黄色大片毛片| 真人一进一出gif抽搐免费| 亚洲中文字幕一区二区三区有码在线看 | 免费在线观看亚洲国产| 欧美一区二区精品小视频在线| 女警被强在线播放| 一进一出抽搐gif免费好疼| 精品欧美国产一区二区三| 国产三级中文精品| 国产成人欧美在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲av电影不卡..在线观看| 国产精品久久久av美女十八| 一卡2卡三卡四卡精品乱码亚洲| 国产高清激情床上av| 免费看日本二区| 日韩 欧美 亚洲 中文字幕| 日本在线视频免费播放| 欧美黄色片欧美黄色片| 久久久精品大字幕| 九色成人免费人妻av| 看黄色毛片网站| 午夜亚洲福利在线播放| 精品福利观看| 在线观看午夜福利视频| 欧美日韩黄片免| 热99re8久久精品国产| 国产av在哪里看| 亚洲欧美精品综合久久99| 国产免费av片在线观看野外av| 亚洲人成网站高清观看| 久久精品影院6| 黑人操中国人逼视频| а√天堂www在线а√下载| 国产成人精品久久二区二区免费| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久人妻蜜臀av| 国内毛片毛片毛片毛片毛片| 国产精品一区二区三区四区免费观看 | 日韩中文字幕欧美一区二区| 国产真实乱freesex| 国产av在哪里看| 国产成人啪精品午夜网站| 久久久久久久久久黄片| 一级毛片高清免费大全| 亚洲欧美一区二区三区黑人| 国产亚洲av高清不卡| 一级作爱视频免费观看| 男女午夜视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲,欧美精品.| 窝窝影院91人妻| 一个人免费在线观看电影 | 午夜福利视频1000在线观看| 三级国产精品欧美在线观看 | 999精品在线视频| 黄片大片在线免费观看| 国产亚洲精品久久久久久毛片| 在线观看免费午夜福利视频| 亚洲精品一区av在线观看| 精品午夜福利视频在线观看一区| 午夜福利欧美成人| 欧美大码av| 亚洲中文av在线| 不卡av一区二区三区| 欧美日韩福利视频一区二区| 国产主播在线观看一区二区| 韩国av一区二区三区四区| 午夜精品在线福利| 一区福利在线观看| 一区二区三区国产精品乱码| 黄片大片在线免费观看| 美女午夜性视频免费| 久久久久久久久中文| 性色avwww在线观看| 久久人人精品亚洲av| 夜夜躁狠狠躁天天躁| 母亲3免费完整高清在线观看| 日本免费一区二区三区高清不卡| 国产av一区在线观看免费| 免费观看人在逋| 午夜福利高清视频| 国产综合懂色| 国产成人精品无人区| 午夜福利在线在线| 国产高清三级在线| 亚洲国产精品成人综合色| 亚洲真实伦在线观看| 精品一区二区三区av网在线观看| 最近最新中文字幕大全免费视频| 97碰自拍视频| 国产成人av激情在线播放| 久久亚洲真实| 色在线成人网| 国产激情偷乱视频一区二区| 亚洲精品久久国产高清桃花| 精华霜和精华液先用哪个|