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

    電磁測深數據地形影響的快速校正

    2016-12-07 07:34:55薛國強閆述陳衛(wèi)營
    地球物理學報 2016年12期
    關鍵詞:比值電阻率剖面

    薛國強, 閆述, 陳衛(wèi)營

    1 中國科學院礦產資源研究重點實驗室, 中國科學院地質與地球物理研究所, 北京 100029 2 江蘇大學計算機科學與通信工程學院, 江蘇鎮(zhèn)江 212013

    ?

    電磁測深數據地形影響的快速校正

    薛國強1, 閆述2, 陳衛(wèi)營1

    1 中國科學院礦產資源研究重點實驗室, 中國科學院地質與地球物理研究所, 北京 100029 2 江蘇大學計算機科學與通信工程學院, 江蘇鎮(zhèn)江 212013

    地形起伏會對電磁法的數據產生一定影響,尤其會影響淺部地層電性結果的準確性.本文通過對地形影響基本規(guī)律、經典比值校正原理的分析, 認為可以用實測標準電阻率替代數值模擬中均勻半空間電阻率,提出一種新的地形影響快速校正方法,即采用小極矩直流電阻率法獲得無地形影響的表層電阻率值,作為地形校正的標準電阻率,以此構造一個新的校正公式.分別對CSAMT和TEM仿真和實際測量資料進行地形校正處理,效果較好.說明新的比值校正公式,是一種快速、有效和實用的校正算法.

    電磁法; 地形影響; 校正; CSAMT; TEM

    1 引言

    在山地電磁法勘探中,地形起伏不僅使觀測點偏離水平位置,還使大地中電磁場的分布發(fā)生變化,導致觀測數據中包含了目標體異常和地形異常,對推斷解釋造成干擾(Jiracek G R,1990;王緒本等, 1999;H?rdt and Scholl,2004;Tang et al.,2011).消弱地形影響,是提高山地勘探精度、取得良好地質效果的關鍵問題之一.地形影響解決方案中,經典的方法是比值校正法(李金銘,2005),該方法先通過一定手段獲得純地形異常,然后用地形造成的異常倍數乘或除以實測數據達到消弱地形影響的效果.在早期,獲得純地形異常的手段有水槽、土槽、導電紙等物理模擬方法,直流電法還可利用場的解析解通過保角變換獲得.隨著計算機和計算技術的發(fā)展,二維和三維純地形正演模擬成為獲得純地形異常的主要方式(Wannamaker et al., 1986; Leppin, 1992; 閆述等, 1996; Xu and Zhou, 1997; 孫懷鳳等, 2013).進一步地,將地形作為二維、三維地質構造的組成部分代入反演也取得了大量的研究成果(Tong and Yang,1990;Yi et al., 2001;吳小平, 2005;肖懷宇, 2006;Haber et al., 2007;雷達, 2010;劉云, 2012;Yang and Oldenburg, 2012;董浩等;2014).二維、三維純地形影響正演和帶地形的反演具有強大的仿真能力,但對計算機內存容量和速度要求極高.在人工源野外勘探中,除必要的點位測量外,發(fā)收之間及周邊的廣大區(qū)域還須有相當密度的測地工作.因此,在數值正反演向實用化發(fā)展的同時,從各個不同的方向消除地形影響(Yi et al., 2001),對電磁勘探數據處理解釋的多樣性是有意義的.本文通過對地形影響基本規(guī)律及經典比值校正原理的分析,用實測標準電阻率替代數值模擬中均勻半空間電阻率,利用地形影響隨頻率降低(或早到晚期)逐漸減弱的特征,以頻點間或時間道間變化速率作為校正系數(閆述, 2015),實現了依賴實測數據進行地形校正的方法.新的比值校正方法增加的工作量很小,為山地探測野外實時數據處理解釋提供了快速簡便的手段.仿真CSAMT和實測TEM校正實例表明,新的比值校正公式是有效的.

    2 地形影響分析

    造成地形影響的有幾何因素和物理因素,主要包括(H?rdt and Scholl., 2004):(1)接收裝置高于水平地表;(2)山體作為導電體置換了空氣,或空氣作為絕緣體置換了原屬水平大地的一部分;(3)接收裝置處于非水平狀態(tài),接收的不僅有磁(電)場的垂直(水平)分量,還有水平(垂直)分量的成分;(4)山體和地下構造的耦合.

    2.1 地形影響的基本規(guī)律

    對以電壓值為基本觀測量的電和電磁法勘探來說,實測電壓或感應電動勢與地形的變化反相,即地勢高處降低,地勢低處增高;地形對頻域(CSAMT或MT)高頻段數據影響較大,對時域(TEM)早期數據影響較大,隨著頻率降低或觀測時間進入晚期,地形影響逐漸減弱.還有,地形變化劇烈處場變化劇烈,地形變化較平緩處場變化平緩;高阻表層地形影響小,低阻表層影響大.這種現象可由電流線在地形作用下的聚集和發(fā)散機理解釋:在導電的大地中,隨地勢的起伏,山峰處電流線發(fā)散(密度小),山谷處電流線聚集(密度大).電場強度與電流密度成正比(如圖1所示):

    J=σE,

    (1)

    電壓與電場強度成正比:

    V=∫lE·dl.

    (2)

    因此實測電壓值幅度的變化與地形的高低變化反相.

    圖1 地形對電流的聚集與發(fā)散作用Fig.1 The gathering and diverging effect of topography to current

    地形影響的電流聚散機理無論對直流電法還是時域、頻域的電磁法勘探都是相同的.時域的早期,相當于頻域的高頻或直流的小極距情況.但是,直流和頻域中的視電阻率與電壓(或電場)成正比,因此視電阻率變化與地形起伏仍然是反相的.時域瞬變場在低阻體中激發(fā)起較大的渦旋場,觀測電壓,也就是感應電動勢EMF(Electromotive Force)隨之增大;在高阻體中激發(fā)的渦旋場較小,EMF隨之減小.這種特性在視電阻率上的反映就是其值與EMF值的大小變化相反,即視電阻率的變化與地形同相,地勢高處視電阻率值增高,地勢低處視電阻率值降低,表現出與頻域和直流電法勘探不一樣的特點.但這在本質上是相同的,都是地形對電流線的發(fā)散與匯集作用造成觀測電壓(電場)幅值的改變.地形影響的基本規(guī)律,是判斷觀測數據是否受到起伏地形影響、和影響程度的準則.

    圖2 邊界元CSAMT電場地形影響模擬(a) 斜坡地形高程; (b) 斜坡地形視電阻率剖面.Fig.2 Boundary element modeling results of CSAMT(a) Slope model; (b) Apparent resistivity profile.

    圖3 陜西黃陵14線中心回線TEM實測地形影響(600 m×600 m)(a) 地形高程; (b) 實測EMF測道圖.Fig.3 The topographic influence on central loop TEM data (600 m×600 m)(a) Topography; (b) Measured EMF curves.

    2.2 地形起伏情況下數值模擬與實測數據

    圖2是二維純地形CSAMT邊界元模擬結果.其中圖2a為地形,圖2b是代表近、中、遠區(qū)的3個頻率沿線剖面圖,源置于坐標原點.圖3a是陜西黃陵煤礦中心回線TEM勘探L14線的地形剖面,圖3b是代表早、中、晚期的實測EMF測道圖.對照兩圖中地形起伏和電場或EMF沿線的變化,大致規(guī)律是:地形下凹處,EMF上凸;地形上凸處,EMF下凹;地形變化劇烈處EMF變化劇烈,隨著頻率的降低,或者從早期向晚期的時間過渡,地形影響逐漸減弱.因此,山地電磁法勘探中,接收點應盡量布置在地形較平緩處,減小地形異常對有用異常的干擾.但山勢起伏總的影響仍然存在,校正是必要的.

    3 比值校正方法

    3.1 經典的比值校正公式

    直流、頻域或時域經典比值校正公式有相同的形式,均以視電阻率為校正參量(李金銘, 2005):

    (3)

    則實測值除以1/n(即乘以n)恢復正常值.

    為避免比值校正時二、三維數值正演的巨大計算量,不因地形校正增加野外勘探工程的測地工作量,根據前述地形影響的基本規(guī)律,依據經典公式(3)的校正原理,給出新的比值校正公式利用實測數據本身達到校正的目的.

    3.2 快速比值校正公式

    新的頻域CSAMT和時域TEM比值校正公式為

    (4)

    式中頻點或時間道序數i從高頻(早期)向低頻(晚期)增加.公式(4)與經典比值校正公式(3)的不同之處為:

    ① 在各測點采用小極矩直流電法或通過巖土實驗室測試,獲得表層標準電阻率ρstandard,并用其替代公式(2)中均勻半空間電阻率ρ1.當表層較厚,如為巨厚黃土層,地表電阻率相差不大時,也可以用測區(qū)內的電測井資料作為各測點的標準電阻率.

    小極矩觀測的目的是避免地形影響,提供校正的標準.可先選適當極矩觀測,當極矩縮小后測得的電阻率變化可忽略時,即認為地形影響可忽略.

    (5)

    其中f(j)為頻率,從高頻到低頻編號1~N.

    (6)

    其中t(i)為時間道,從早期到晚期編號1~N.

    為提高系數C(i)精度,可加密頻點或時間道.需要說明的是,校正系數C(i)也可以是非線性的,這就需要通過鉆孔、地質資料、地下挖掘等其他手段確定不同層位的標準電阻率,此時校正結果的精度更高.

    3.3 快速地形校正適用范圍

    如果起伏地表下,目標地層平緩穩(wěn)定,后續(xù)數據解釋中可作為一維情況處理,那么新的、快速簡便的地形校正是有意義的.如果探測深度范圍內地質構造不能作為一維地層近似,采用帶地形的三維反演就是必要的.巨大的計算工作量和測地工作量,是復雜地質構造勘探取得良好地質效果的必要代價.

    快速校正法假設低頻或晚期視電阻率已不受地形影響,并以此和標準電阻率一同作為校正的基準.若極端情況下,未能觀測到不隨地形變化的最低或最晚數據,或者低頻端、特別是晚時間道信噪比不夠,無法利用時,公式(4)的校正效果下降.如對探測有較高的精度要求,仍需采用公式(3)校正,數值模擬不可避免.

    4 校正實例

    4.1 CSAMT數據校正

    圖4a是較圖2斜坡地形變化更為劇烈的正弦地形,圖4b顯示了典型的地形影響特征.對CSAMT仿真數據校正時,公式(3)中的ρstandard即為均勻半空間電阻率ρ1=100 Ωm.表1是所用的全部頻率列表,校正后的結果示于圖4.眾所周知,CSAMT單分量遠區(qū)視電阻率是大地真電阻率,是近區(qū)視電阻率的2倍,中區(qū)視電阻率隨頻率降低從遠區(qū)過渡到近區(qū)(陳明生和閆述, 1995).將圖4與圖2較后可以看到,代表遠區(qū)的f=7355 Hz視電阻率剖面為100 Ωm的直線,這是大地真實電阻率;代表中區(qū)的f=81.45 Hz視電阻率剖面成一斜線,代表近區(qū)的f=1.274 Hz視電阻率剖面為50 Ωm的直線.這表明,地形響應被消除,快速校正公式(4)對仿真的地形影響是有效的.

    圖4 仿真CSAMT地形影響校正(a) 地形剖面; (b) 校正前的CSAMT視電阻率剖面; (c) 校正后的CSAMT視電阻率剖面.Fig.4 Topographic influence correction for synthetic data of CSAMT(a) Topography profile; (b) Apparent resistivity curves before correction; (c) Apparent resistivity curves after correction.

    頻點觀測頻率f/Hz頻點觀測頻率f/Hz頻點觀測頻率f/Hz頻點觀測頻率f/Hz1104209651.51740.73252.5452*7355*10459.81828.75261.7983521011325.81920.3727*1.274*4367812229.82014.38280.8995260513167.02110.19290.4506183914115.0227.19300.2257130315*81.45*235.0958919.51657.5243.595

    注:加*數字為實例采用的頻點和頻率.

    圖5 實測TEM地形影響校正(a) 沿TEM測線的地形剖面; (b) 校正前的TEM視電阻率剖面; (c) 校正后的TEM視電阻率剖面.Fig.5 Topographic influence correction of field TEM data(a) Topography profile; (b) Apparent resistivity curves before correction; (c) Apparent resistivity curves after correction.

    4.2 TEM實測數據校正

    陜西黃陵煤礦地表為黃土沖蝕地形,下伏地層平緩、穩(wěn)定,傾角在5°以內,符合快速地形校正的適用條件.所用中心回線TEM發(fā)射線框600 m×600 m、接收線圈等效面積100 m2,觀測時長30 ms、20道.圖 5a是沿線地形剖面,為表示清晰起見,圖5b的實測TEM視電阻率剖面用了其中的t=0.318 ms、t=0.504 ms、t=1.005 ms、t=15.985 ms等4個時間道.表3列出了通過小極矩直流電阻率法獲得的各測點地表電阻率值,作為公式(4)中的標準電阻率ρstandard.時域TEM視電阻率和地形(圖5a)呈同相變化,即地形下降時視電阻率下降、地形上升時視電阻率上升,地形變化劇烈時(如圖5b中的240 m、960 m處)視電阻率的變化也劇烈.根據地形的時域TEM響應特征,可以判斷t=15.985 ms時刻地形影響消失,故取t(1)=0.318 ms、t(N)=15.985 ms,連同實測數據和表3作為標準電阻率的地表電阻率值代入公式(4)和(6)后,得到校正后的TEM視電阻率剖面(圖5c).由圖5可知,實測TEM視電阻率除了純地形響應外,還耦合了地質結構本身的變化.但是地形影響校正僅去除地形影響,并保留地質結構響應,這種新的校正方法取得了較好的效果.

    表2 實測TEM山地勘探時間道

    注:加*數字為實例采用的時間道和觀測時間.

    表3 小極矩直流電阻率法獲得的各測點標準電阻率值

    5 結論

    本文通過對地形影響基本規(guī)律、經典比值校正原理的分析,用實測標準電阻率替代數值模擬中均勻半空間電阻率,利用地形影響隨頻率降低(或早到晚期)逐漸減弱的特征,以頻點間或時間道間變化速率作為校正系數,針對地形對電磁法數據的影響提出了一種新的快速校正公式.本文提出的校正方法在實施過程中,除了測點以外,無須額外的地形高程測量,實現了依賴實測數據進行地形校正.同經典比值校正法一樣,僅校正地形影響,保留了地質結構的響應.新的比值校正方法增加的工作量很小,為山地探測野外實時數據處理解釋提供了快速簡便的手段.仿真CSAMT和實測TEM校正實例表明,新的比值校正公式是有效的.

    Chen M S, Yan S. 1995. Problems in Frequency Electromagnetic Soundings (in Chinese). Beijing: Geology Press.

    Dong H, Wei W B, Ye G F, et al. 2014. Study of three-dimensional magnetotelluric inversion including surface topography based on finite-difference method.ChineseJ.Geophys. (in Chinese), 57(3): 939-952, doi: 10.6038/cjg20140323. Haber E, Oldenburg D W, Shekhtman R. 2007. Inversion of time domain three-dimensional electromagnetic data.GeophysicalJournalInternational, 171(2): 550-564.

    H?rdt A, Scholl C. 2004. The effect of local distortions on time-domain electromagnetic measurements.Geophysics, 69(1): 87-96.

    Jiracek G R. 1990. Near-surface and topographic distortions in electromagnetic induction.SurveysinGeophysics, 1990, 11(2-3): 163-203.

    Lei D. 2010. Studies and applications of 2-D CSAMT modeling and inversion with a dipole source and topography.ChineseJ.Geophys. (in Chinese),53(4): 982-993, doi: 10.3969/j.issn.0001-5733.2010.04.023.

    Leppin M. 1992. Electromagnetic modeling of 3-D sources over 2-D inhomogeneities in the time domain.Geophysics, 57(8): 994-1003.

    Li J M. 2005. Geo-electric Field and Electrical Exploration Methods (in Chinese). Beijing; : Geology Press.

    Liu Y. 2012. Two-dimensional numerical modeling for topography Magnetotelluric/time-domain transient electromagnetic and direct inverse method [Ph.D. thesis] (in Chinese)., Chengdu: Chengdu University of Technology. Sun H F, Li X, Li S C, et al. 2013. Three-dimensional FDTD modeling of TEM excited by a loop source considering ramp time.ChineseJ.Geophys. (in Chinese), 56(3): 1049-1064, doi: 10.6038/cjg20130333.

    Tang X G, Hu W B, Yan L J. 2011. Topographic effects on long offset transient electromagnetic response.AppliedGeophysics, 8(4): 277-284. Tong L T, Yang C H. 1990. Incorporation of topography into two-dimensional resistivity inversion.Geophysics, 55(3): 354-361. Wang X B, Li Y N, Gao Y C. 1999. Two dimensional topographic responses in magneto telluric sounding and its correction methods.ComputingTechniquesforGeophysicalandGeochemicalExploration(in Chinese), 21(4): 327-332. Wannamaker P E, Stodt J A, Rijo L. 1986. Two-dimensional topographic responses in magnetotellurics modeled using finite elements.Geophysics, 51(11): 2131-2144.

    Wu X P. 2005. 3-D resistivity inversion under the condition of uneven terrain.ChineseJ.Geophys. (in Chinese), 48(4): 932-936. Xiao H Y. 2006. Three-dimensional numerical modeling considering the topography of TEM [Master thesis] (in Chinese). Beijing: China Universiity of Geosciences (Beijing).

    Xu S Z, Zhou H. 1997. Modelling the 2D terrain effect on MT by the boundary-element method.GeophysicalPprospecting, 45(6): 931-943.

    Yan S. 2015. One ratio method for correcting the topographic effect in electromagnetric exploration (in Chinese). Chinese Patent, ZL201310039614.6.

    Yan S, Chen M S, Li Z M. 1996. Boundary element for method forward solution of two dimensional topography effection electromagnetic frequency sounding.ComputingTechniquesforGeophysicalandGepochemicalExploration(in Chinese), 22(4): 4533-4842, 47.

    Yang D K, Oldenburg D W. 2012. Three-dimensional inversion of airborne time-domain electromagnetic data with applications to a porphyry deposit.Geophysics, 77(2): B23-B34.

    Yi M J, Kim J H, Song Y, et al. 2001. Three-dimensional imaging of subsurface structures using resistivity data.GeophysicalProspecting, 49(4): 483-497.

    附中文參考文獻

    陳明生, 閆述. 1995. 論頻率測深應用中的幾個問題. 北京: 地質出版社, 1995.

    董浩, 魏文博, 葉高峰等. 2014. 基于有限差分正演的帶地形三維大地電磁反演方法. 地球物理學報, 57(3): 939-952, doi: 10.6038/cjg20140323.

    雷達. 2010. 起伏地形下CSAMT二維正反演研究與應用. 地球物理學報, 53(4): 982-993, doi: 10.3969/j.issn.0001-5733.2010.04.023.

    李金銘. 2005. 地電場與電法勘探. 北京: 地質出版社.

    劉云. 2012. 起伏地形大地電磁、時間域瞬變電磁二維數值模擬及直接反演法[博士論文]. 成都: 成都理工大學.

    孫懷鳳, 李貅, 李術才等. 2013. 考慮關斷時間的回線源激發(fā)TEM三維時域有限差分正演. 地球物理學報, 56(3): 1049-1064, doi: 10.6038/cjg20130333.

    王緒本, 李永年, 高永才. 1999. 大地電磁測深二維地形影響及其校正方法研究. 物探化探計算技術, 21(4).: 327-332.

    吳小平. 2005. 非平坦地形條件下電阻率三維反演. 地球物理學報, 48(4): 932-936.

    肖懷宇. 2006. 帶地形的瞬變電磁法三維數值模擬[碩士論文]. 北京: 中國地質大學(北京).

    閆述. 2015. 一種校正電磁勘探中地形影響的比值方法. ZL201310039614.6.

    閆述, 陳明生, 李志民. 1996. 頻率測深二維地形影響的邊界元素法正演模擬. 物探化探計算技術, 22(4): 4533-4842, 47.

    (本文編輯 汪海英)

    A fast topographic correction method for electromagnetic data

    XUE Guo-Qiang1, YAN Shu2, CHEN Wei-Ying1

    1KeyLaboratoryofMineralResources,InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China2SchoolofComputerScienceandTelecommunicationEngineering,JiangsuUniversity,JiangsuZhenjiang212013,China

    Topography generally affects the electromagnetic data and will mislead the results of shallow layer′s resistivity. In this paper, by analyzing the basic law of topography effect and classic ratio correction principle, it is feasible to replace the resistivity of uniform half space from numerical simulation by standard measured resistivity. We proposed a new fast method to correct the topographic effect in electromagnetic data. We obtain the surface resistivity without topography by using a small DC configuration, and then take it as standard resistivity to construct a new correction formula. Both synthetic and field data of CSAMT and TEM were processed using this method. Results show that the new method is a fast, effective and practical tool for topographic correction of electromagnetic data.

    Electromagnetic method; Topographic effect; Correction; CSAMT; TEM

    10.6038/cjg20161202.Xue G Q, Yan S, Chen W Y. 2016. A fast topographic correction method for electromagnetic data.ChineseJ.Geophys. (in Chinese),59(12):4408-4413,doi:10.6038/cjg20161202.

    國家自然科學基金(41474095),中國科學院礦產資源研究重點實驗室開放課題,國土資源部煤炭資源勘查與綜合利用重點實驗室開放課題資助.

    薛國強,男,1966年生,研究員,主要從事瞬變電磁法理論與應用研究. E-mail: ppxueguoqiang@163.com

    10.6038/cjg20161202

    P631

    2016-02-04,2016-10-09收修定稿

    薛國強, 閆述, 陳衛(wèi)營. 2016. 電磁測深數據地形影響的快速校正. 地球物理學報,59(12):4408-4413,

    猜你喜歡
    比值電阻率剖面
    三點法定交叉剖面方法
    ——工程地質勘察中,一種做交叉剖面的新方法
    基于曲線擬合的投棄式剖面儀電感量算法
    電子測試(2017年12期)2017-12-18 06:35:46
    比值遙感蝕變信息提取及閾值確定(插圖)
    河北遙感(2017年2期)2017-08-07 14:49:00
    復雜多約束條件通航飛行垂直剖面規(guī)劃方法
    三維電阻率成像與高聚物注漿在水閘加固中的應用
    不同應變率比值計算方法在甲狀腺惡性腫瘤診斷中的應用
    隨鉆電阻率測井的固定探測深度合成方法
    海洋可控源電磁場視電阻率計算方法
    粉煤灰摻量對水泥漿體電阻率與自收縮的影響
    雙電機比值聯動控制系統(tǒng)
    国产真实伦视频高清在线观看 | 国产69精品久久久久777片| 日韩有码中文字幕| 成人av在线播放网站| 午夜激情福利司机影院| 精品无人区乱码1区二区| aaaaa片日本免费| 国产v大片淫在线免费观看| 免费黄网站久久成人精品 | 日韩中字成人| 波多野结衣高清作品| 欧美性感艳星| 亚洲黑人精品在线| 日本一二三区视频观看| 色综合站精品国产| 国产三级黄色录像| 日韩欧美精品v在线| www日本黄色视频网| 无遮挡黄片免费观看| 欧美+亚洲+日韩+国产| 久久人人精品亚洲av| 国产真实伦视频高清在线观看 | 在线免费观看的www视频| 90打野战视频偷拍视频| 真人一进一出gif抽搐免费| 免费黄网站久久成人精品 | 一本一本综合久久| 国产又黄又爽又无遮挡在线| 欧美日韩黄片免| 91麻豆精品激情在线观看国产| 亚洲国产精品成人综合色| 中出人妻视频一区二区| 看黄色毛片网站| 午夜福利高清视频| 一个人免费在线观看的高清视频| 夜夜夜夜夜久久久久| 亚洲中文日韩欧美视频| 国产野战对白在线观看| 岛国在线免费视频观看| 中文字幕免费在线视频6| 伊人久久精品亚洲午夜| 最新中文字幕久久久久| 中国美女看黄片| 国产成人啪精品午夜网站| 黄色配什么色好看| 精品国产三级普通话版| 久久久久久大精品| 久久人人爽人人爽人人片va | 日韩 亚洲 欧美在线| 国产亚洲精品久久久com| 成人av在线播放网站| 亚洲黑人精品在线| 九色国产91popny在线| 免费看日本二区| 欧美成人免费av一区二区三区| 亚洲一区高清亚洲精品| av黄色大香蕉| 人人妻人人看人人澡| 国产乱人伦免费视频| 午夜日韩欧美国产| 日日摸夜夜添夜夜添av毛片 | 亚洲国产色片| 男女床上黄色一级片免费看| 给我免费播放毛片高清在线观看| 免费一级毛片在线播放高清视频| 欧美日本亚洲视频在线播放| 日韩中字成人| 99热只有精品国产| 亚洲18禁久久av| 久久久久久国产a免费观看| 色播亚洲综合网| 好男人在线观看高清免费视频| 免费在线观看亚洲国产| 亚洲久久久久久中文字幕| 欧美最黄视频在线播放免费| 熟妇人妻久久中文字幕3abv| 身体一侧抽搐| 国产野战对白在线观看| 欧美精品国产亚洲| 丝袜美腿在线中文| 欧美日韩福利视频一区二区| 国内久久婷婷六月综合欲色啪| 在线播放无遮挡| 日韩欧美 国产精品| 91麻豆精品激情在线观看国产| 俺也久久电影网| 男女那种视频在线观看| 欧美日韩黄片免| 九色成人免费人妻av| 一级作爱视频免费观看| 欧美+日韩+精品| 欧美在线黄色| 青草久久国产| 9191精品国产免费久久| avwww免费| 亚洲精品久久国产高清桃花| 极品教师在线视频| 精品久久久久久久久久免费视频| 熟妇人妻久久中文字幕3abv| 欧美三级亚洲精品| 美女大奶头视频| 国产精品一区二区三区四区久久| 天美传媒精品一区二区| 精品欧美国产一区二区三| 精品日产1卡2卡| 99久久精品一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 国产男靠女视频免费网站| 久久中文看片网| 国产黄色小视频在线观看| 国产淫片久久久久久久久 | 99久久九九国产精品国产免费| 九九久久精品国产亚洲av麻豆| 午夜老司机福利剧场| 午夜精品一区二区三区免费看| 成人毛片a级毛片在线播放| 国产精品日韩av在线免费观看| 免费在线观看成人毛片| 亚洲色图av天堂| 国产精品免费一区二区三区在线| 日韩有码中文字幕| 搡老熟女国产l中国老女人| 午夜福利欧美成人| 久久久久国内视频| 99精品在免费线老司机午夜| 午夜免费成人在线视频| 欧美在线一区亚洲| 可以在线观看毛片的网站| 亚洲精品在线观看二区| 成人无遮挡网站| 九色国产91popny在线| 精品人妻偷拍中文字幕| 国产av不卡久久| 国产精华一区二区三区| 狂野欧美白嫩少妇大欣赏| 欧美丝袜亚洲另类 | 国产 一区 欧美 日韩| 人妻制服诱惑在线中文字幕| 国产高潮美女av| 午夜福利在线在线| 亚洲美女视频黄频| 床上黄色一级片| 欧美xxxx黑人xx丫x性爽| 男人狂女人下面高潮的视频| 中文字幕av在线有码专区| 亚洲中文字幕一区二区三区有码在线看| eeuss影院久久| 一边摸一边抽搐一进一小说| 好看av亚洲va欧美ⅴa在| 男女下面进入的视频免费午夜| 久久性视频一级片| 欧美绝顶高潮抽搐喷水| avwww免费| 亚洲av免费高清在线观看| 波多野结衣巨乳人妻| 国产午夜精品久久久久久一区二区三区 | 夜夜躁狠狠躁天天躁| 日韩欧美国产一区二区入口| 日韩高清综合在线| 搡老熟女国产l中国老女人| 男女做爰动态图高潮gif福利片| 午夜精品在线福利| 日韩国内少妇激情av| 久久久久久久久大av| 99在线视频只有这里精品首页| 亚洲成av人片在线播放无| 国产高清激情床上av| 国产精品嫩草影院av在线观看 | 亚洲午夜理论影院| 亚洲中文字幕日韩| 黄色女人牲交| 岛国在线免费视频观看| 亚洲人与动物交配视频| 亚洲欧美日韩高清在线视频| 精品不卡国产一区二区三区| 国产不卡一卡二| 一个人免费在线观看电影| 国产精品嫩草影院av在线观看 | 观看美女的网站| 日本撒尿小便嘘嘘汇集6| 丁香六月欧美| av在线蜜桃| 久久精品影院6| 老司机午夜十八禁免费视频| 久久精品91蜜桃| 亚洲性夜色夜夜综合| 亚洲第一电影网av| 日韩欧美精品v在线| 色播亚洲综合网| 久久久久精品国产欧美久久久| 精品福利观看| 国产精品久久久久久亚洲av鲁大| 国产成人福利小说| 久久久国产成人免费| 十八禁国产超污无遮挡网站| 精品一区二区三区视频在线观看免费| 狠狠狠狠99中文字幕| 久久伊人香网站| 午夜免费男女啪啪视频观看 | 1024手机看黄色片| 免费黄网站久久成人精品 | 久久午夜福利片| 免费一级毛片在线播放高清视频| avwww免费| 1024手机看黄色片| 三级男女做爰猛烈吃奶摸视频| 国产精品三级大全| 九色国产91popny在线| 亚洲精品一区av在线观看| 性插视频无遮挡在线免费观看| 国产在视频线在精品| 国产亚洲精品久久久久久毛片| 最新中文字幕久久久久| 一个人免费在线观看电影| 美女高潮的动态| 91狼人影院| 3wmmmm亚洲av在线观看| 免费在线观看成人毛片| 悠悠久久av| 69av精品久久久久久| 免费人成在线观看视频色| 美女免费视频网站| 欧美xxxx性猛交bbbb| 亚洲经典国产精华液单 | 亚洲成av人片在线播放无| 国产伦精品一区二区三区视频9| 18禁黄网站禁片午夜丰满| 日韩精品中文字幕看吧| 搡老岳熟女国产| 日韩精品中文字幕看吧| 久久国产精品人妻蜜桃| 精品不卡国产一区二区三区| 丰满的人妻完整版| 男女视频在线观看网站免费| 亚洲久久久久久中文字幕| 看十八女毛片水多多多| 国内揄拍国产精品人妻在线| 欧美成人一区二区免费高清观看| 高清在线国产一区| ponron亚洲| 亚洲熟妇熟女久久| 午夜福利成人在线免费观看| 日韩免费av在线播放| 老鸭窝网址在线观看| 91在线观看av| 国产高清视频在线观看网站| 日韩欧美国产在线观看| 最新中文字幕久久久久| 久久久久性生活片| 日本免费一区二区三区高清不卡| 给我免费播放毛片高清在线观看| 又黄又爽又免费观看的视频| 少妇裸体淫交视频免费看高清| 亚洲人成网站高清观看| 亚洲精品成人久久久久久| 别揉我奶头~嗯~啊~动态视频| 亚洲av二区三区四区| 亚洲性夜色夜夜综合| 亚洲在线自拍视频| 国产精品,欧美在线| 日本黄色片子视频| 波多野结衣高清无吗| 男插女下体视频免费在线播放| 免费看光身美女| 丰满乱子伦码专区| 亚洲在线观看片| 在线观看一区二区三区| 蜜桃亚洲精品一区二区三区| av在线观看视频网站免费| 久久久色成人| 日韩中文字幕欧美一区二区| 18美女黄网站色大片免费观看| 国产精品亚洲美女久久久| 亚洲内射少妇av| 国产精品永久免费网站| 精华霜和精华液先用哪个| 中文字幕精品亚洲无线码一区| 欧美三级亚洲精品| 国产精品久久视频播放| 国产毛片a区久久久久| 国产一区二区三区在线臀色熟女| 国产在线男女| 成年女人永久免费观看视频| 99精品在免费线老司机午夜| 长腿黑丝高跟| 国产成人a区在线观看| 男女下面进入的视频免费午夜| 久久99热6这里只有精品| 身体一侧抽搐| 性色avwww在线观看| 国产精品一区二区三区四区免费观看 | 少妇熟女aⅴ在线视频| or卡值多少钱| 中文字幕高清在线视频| 精品久久国产蜜桃| 脱女人内裤的视频| 免费高清视频大片| 久久久精品大字幕| 亚洲av美国av| 九色国产91popny在线| 婷婷六月久久综合丁香| 亚洲不卡免费看| 搡老岳熟女国产| 久久久久性生活片| 国产爱豆传媒在线观看| 久久久国产成人免费| 国产精品久久久久久久电影| 国模一区二区三区四区视频| 亚洲专区国产一区二区| 18美女黄网站色大片免费观看| 欧美日韩瑟瑟在线播放| 一夜夜www| 欧美日本视频| 国产精品女同一区二区软件 | 亚洲无线观看免费| 亚洲综合色惰| av在线老鸭窝| 成人av一区二区三区在线看| 天美传媒精品一区二区| 九色国产91popny在线| 老熟妇仑乱视频hdxx| av视频在线观看入口| 国产精品一区二区免费欧美| 自拍偷自拍亚洲精品老妇| 18美女黄网站色大片免费观看| 免费观看精品视频网站| av在线观看视频网站免费| 亚洲精华国产精华精| 啦啦啦观看免费观看视频高清| 最好的美女福利视频网| 午夜老司机福利剧场| 在线观看舔阴道视频| 国产私拍福利视频在线观看| 桃色一区二区三区在线观看| 亚洲成人免费电影在线观看| 18禁裸乳无遮挡免费网站照片| 激情在线观看视频在线高清| 国产精品女同一区二区软件 | 免费黄网站久久成人精品 | 91在线精品国自产拍蜜月| 国内久久婷婷六月综合欲色啪| 国产视频一区二区在线看| 露出奶头的视频| 亚洲第一电影网av| 我的老师免费观看完整版| а√天堂www在线а√下载| 好看av亚洲va欧美ⅴa在| or卡值多少钱| 欧美日韩国产亚洲二区| 最近中文字幕高清免费大全6 | 午夜两性在线视频| 给我免费播放毛片高清在线观看| 麻豆久久精品国产亚洲av| 欧美成狂野欧美在线观看| 亚洲 欧美 日韩 在线 免费| 好男人在线观看高清免费视频| 国产极品精品免费视频能看的| 夜夜爽天天搞| 午夜日韩欧美国产| 久久欧美精品欧美久久欧美| 亚洲成人免费电影在线观看| 久久亚洲精品不卡| 国产亚洲精品久久久com| 亚洲欧美日韩高清在线视频| 久久精品夜夜夜夜夜久久蜜豆| 性插视频无遮挡在线免费观看| 国产精品一区二区免费欧美| 男女下面进入的视频免费午夜| 日韩欧美一区二区三区在线观看| 成年女人毛片免费观看观看9| 亚洲 国产 在线| 日本 欧美在线| 国产精品99久久久久久久久| 我的老师免费观看完整版| 日本 欧美在线| 亚洲国产精品合色在线| 国产精品日韩av在线免费观看| 高清在线国产一区| 午夜福利18| 一级作爱视频免费观看| 亚洲精品一卡2卡三卡4卡5卡| 中文字幕熟女人妻在线| 真实男女啪啪啪动态图| 久久国产乱子伦精品免费另类| 久久久久国内视频| 国产伦一二天堂av在线观看| 亚洲av成人av| 免费黄网站久久成人精品 | 欧美成人一区二区免费高清观看| 国产探花极品一区二区| 久久婷婷人人爽人人干人人爱| 亚洲人成伊人成综合网2020| 一级黄色大片毛片| 国产黄a三级三级三级人| 男人狂女人下面高潮的视频| 成人毛片a级毛片在线播放| 色精品久久人妻99蜜桃| 韩国av一区二区三区四区| 国产91精品成人一区二区三区| 老司机午夜福利在线观看视频| ponron亚洲| 国模一区二区三区四区视频| 少妇熟女aⅴ在线视频| 变态另类丝袜制服| 99热这里只有是精品50| 日本五十路高清| 日韩欧美免费精品| 中文字幕精品亚洲无线码一区| 99热这里只有精品一区| 亚洲美女黄片视频| 精品人妻熟女av久视频| 婷婷亚洲欧美| 国产一区二区三区在线臀色熟女| 午夜久久久久精精品| av国产免费在线观看| 麻豆国产av国片精品| 亚洲色图av天堂| 18禁在线播放成人免费| 日日摸夜夜添夜夜添小说| 亚洲美女视频黄频| 国产成人欧美在线观看| 国产高潮美女av| a在线观看视频网站| 国产伦一二天堂av在线观看| 免费观看人在逋| 757午夜福利合集在线观看| 免费人成视频x8x8入口观看| 欧美三级亚洲精品| bbb黄色大片| 亚洲美女黄片视频| 直男gayav资源| 亚洲自偷自拍三级| 久久久久久久久中文| 欧美日韩乱码在线| 麻豆成人av在线观看| 亚洲最大成人手机在线| 亚洲av成人不卡在线观看播放网| 成人鲁丝片一二三区免费| av欧美777| 国产色爽女视频免费观看| 亚洲内射少妇av| 中文资源天堂在线| 欧美在线一区亚洲| 999久久久精品免费观看国产| 无人区码免费观看不卡| 美女大奶头视频| 高清在线国产一区| 国产一区二区在线av高清观看| 在线国产一区二区在线| 亚洲专区国产一区二区| 青草久久国产| 97碰自拍视频| 亚洲人成网站在线播| 精品乱码久久久久久99久播| 看片在线看免费视频| 欧美日韩乱码在线| 国产亚洲精品久久久久久毛片| 简卡轻食公司| 18禁黄网站禁片免费观看直播| av黄色大香蕉| 欧美日韩中文字幕国产精品一区二区三区| 成人高潮视频无遮挡免费网站| www.www免费av| 悠悠久久av| 国产69精品久久久久777片| 国产精品av视频在线免费观看| 淫妇啪啪啪对白视频| a级毛片免费高清观看在线播放| 国产伦一二天堂av在线观看| 18美女黄网站色大片免费观看| 亚洲成人久久性| 欧美丝袜亚洲另类 | 亚洲精品粉嫩美女一区| 人人妻人人澡欧美一区二区| 亚洲精品成人久久久久久| 亚洲欧美激情综合另类| 又粗又爽又猛毛片免费看| 好男人在线观看高清免费视频| 18+在线观看网站| 长腿黑丝高跟| 午夜久久久久精精品| 国产欧美日韩一区二区三| 欧美日韩亚洲国产一区二区在线观看| 神马国产精品三级电影在线观看| 日韩 亚洲 欧美在线| 久久国产乱子免费精品| 乱码一卡2卡4卡精品| 色综合站精品国产| 国产毛片a区久久久久| 亚洲av免费高清在线观看| 国产av在哪里看| 亚洲久久久久久中文字幕| 搡女人真爽免费视频火全软件 | 18禁在线播放成人免费| 在线观看免费视频日本深夜| 国产中年淑女户外野战色| 欧美bdsm另类| 窝窝影院91人妻| 最近中文字幕高清免费大全6 | 村上凉子中文字幕在线| 天堂影院成人在线观看| 噜噜噜噜噜久久久久久91| 亚洲经典国产精华液单 | 亚洲精品成人久久久久久| 又黄又爽又刺激的免费视频.| 中文字幕熟女人妻在线| 精品不卡国产一区二区三区| 狂野欧美白嫩少妇大欣赏| 午夜福利欧美成人| 日日夜夜操网爽| 一夜夜www| 精品一区二区三区视频在线观看免费| 久久99热6这里只有精品| 国产精品久久久久久精品电影| 波多野结衣巨乳人妻| 男女视频在线观看网站免费| 91字幕亚洲| 一卡2卡三卡四卡精品乱码亚洲| 欧美精品国产亚洲| 成人高潮视频无遮挡免费网站| 久久午夜亚洲精品久久| 久久国产乱子伦精品免费另类| av在线老鸭窝| 国产黄色小视频在线观看| 精品一区二区三区视频在线观看免费| 久久久国产成人精品二区| 亚洲av免费在线观看| 中文资源天堂在线| 男人的好看免费观看在线视频| 2021天堂中文幕一二区在线观| 国产综合懂色| 制服丝袜大香蕉在线| 国产午夜福利久久久久久| 国产单亲对白刺激| 99热只有精品国产| 一个人看视频在线观看www免费| 一本一本综合久久| eeuss影院久久| 国产精品电影一区二区三区| 亚洲黑人精品在线| 亚洲精品亚洲一区二区| av黄色大香蕉| АⅤ资源中文在线天堂| av专区在线播放| 国产精品自产拍在线观看55亚洲| 色5月婷婷丁香| 国产一区二区三区在线臀色熟女| 亚洲成av人片在线播放无| 男女之事视频高清在线观看| 99国产精品一区二区三区| av视频在线观看入口| 成人国产综合亚洲| 久久性视频一级片| 美女免费视频网站| 国产欧美日韩精品一区二区| 夜夜爽天天搞| 婷婷精品国产亚洲av在线| 黄色丝袜av网址大全| 黄色视频,在线免费观看| 国产精品久久久久久久电影| 男人狂女人下面高潮的视频| 变态另类丝袜制服| 波野结衣二区三区在线| 99热这里只有是精品在线观看 | 在线观看舔阴道视频| 给我免费播放毛片高清在线观看| 美女免费视频网站| 国产伦精品一区二区三区四那| 色精品久久人妻99蜜桃| 自拍偷自拍亚洲精品老妇| 淫妇啪啪啪对白视频| 99国产精品一区二区蜜桃av| 97碰自拍视频| 国产色爽女视频免费观看| 51午夜福利影视在线观看| 老熟妇仑乱视频hdxx| 在线播放国产精品三级| 免费看美女性在线毛片视频| 两个人的视频大全免费| 久久精品国产亚洲av天美| 国产精品精品国产色婷婷| 熟女电影av网| 亚洲欧美精品综合久久99| 国产一区二区在线av高清观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成a人片在线一区二区| 亚洲18禁久久av| 国产白丝娇喘喷水9色精品| 国产中年淑女户外野战色| 97超视频在线观看视频| 黄色一级大片看看| 国内精品久久久久久久电影| 丰满人妻一区二区三区视频av| 少妇熟女aⅴ在线视频| 国产成人av教育| 欧美黄色片欧美黄色片| 婷婷亚洲欧美| 国产视频内射| www.999成人在线观看| 国语自产精品视频在线第100页| x7x7x7水蜜桃| 亚州av有码| 国产白丝娇喘喷水9色精品| 亚洲成a人片在线一区二区| 十八禁网站免费在线| 波野结衣二区三区在线| 欧美中文日本在线观看视频| 国产av不卡久久| 小说图片视频综合网站| 亚洲第一欧美日韩一区二区三区| 成人特级av手机在线观看| 欧美不卡视频在线免费观看| 1024手机看黄色片| 日韩欧美精品v在线| 窝窝影院91人妻|