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

    基于FDTD的空間步進頻雷達煤層異質(zhì)區(qū)探測正演分析

    2022-01-25 06:49羅彩禎楊志偉賈遠航
    現(xiàn)代電子技術 2022年1期
    關鍵詞:探地富水異質(zhì)

    賀 順,羅彩禎,楊志偉,賈遠航

    (1.西安科技大學 通信與信息工程學院,陜西 西安 710054;2.西安電子科技大學 雷達信號處理重點實驗室,陜西 西安 710071)

    0 引 言

    現(xiàn)如今由煤層異質(zhì)區(qū)如采空區(qū)、富水區(qū)、瓦斯突出區(qū)等引發(fā)的礦井事故多有發(fā)生,了解掘進面前方煤體的結構對煤礦安全生產(chǎn)尤為重要。目前,多種地球物理勘探技術已被廣泛應用于煤礦井下地質(zhì)預測中,如礦井直流電法、瞬變電磁法、地震勘探法等。礦井直流法受地形的影響較大,需要花費大量的時間進行設備的鋪設;瞬變電磁法存在淺部探測盲區(qū)且容易受到鐵磁介質(zhì)的影響;在多種勘探技術中,探地雷達有抗干擾能力強、受地形的影響較小、操作方便、不需要進行龐雜的設備布置等優(yōu)點。許多學者對探地雷達應用于煤礦井下探測進行了研究。文獻[2?3]使用地質(zhì)雷達探測井下瓦斯富集區(qū)并對其特征進行研究;文獻[4]使用GprMax和Matlab數(shù)值軟件對井下二維空間磁場模型進行正演模擬,研究不同激勵源與中心頻率對模擬結果的影響。文獻[5]使用超寬帶探地雷達開展了煤層厚度的探測試驗。井下探測使用的探地雷達多為單發(fā)單收雷達系統(tǒng),這種系統(tǒng)雖然操作簡單,探測方便,但需要多次移動進行多點測量,可獲得的信息量有限。

    多輸入多輸出探地雷達探測一次可獲得多組數(shù)據(jù),相比單站雷達能獲得更多的信息量,更加充分地利用了回波信號的空時信息,并提高了系統(tǒng)的空間分辨能力。文獻[6?7]先后對陣列式探地雷達目標成像問題以及陣列式探地雷達信號極化場特征進行了研究;文獻[8]對車載陣列探地雷達實時檢測地雷進行了研究。目前,多輸入多輸出探地雷達多應用于檢測地下藏雷、管線結構、路面地下空洞等方面,在煤礦井下探測的應用較少。

    探地雷達通過增大信號帶寬和天線合成孔徑來提高雷達探測的距離分辨力和方位分辨力。硬件發(fā)展的水平限制了超寬帶信號的產(chǎn)生。多輸入多輸出體制下的空間步進頻雷達通過在相鄰的陣元上依次發(fā)射頻率步進的信號,在接收端合成一個寬帶信號,解決了探地雷達系統(tǒng)中目標的距離分辨力與雷達信號帶寬之間的矛盾,同時相比單站雷達方位分辨力又有所提高。本文主要研究多輸入多輸出雷達體制下的空間步進頻信號在煤礦井下環(huán)境的傳播,建立二維煤礦井下異質(zhì)區(qū)環(huán)境模型,使用時域有限差分方法(Finite?difference Time?domain,F(xiàn)DTD)進行空間步進頻信號探測異質(zhì)區(qū)的數(shù)值模擬,分析不同異質(zhì)區(qū)模型下的回波結果,得到對應模型下的傳播規(guī)律。

    1 空間步進頻信號模型

    空間步進頻雷達對煤礦井下異質(zhì)區(qū)進行探測時,將雷達天線陣列放置在井下巷道中。天線陣列同時平行于巷道地面與煤礦開采掘進面,如圖1所示。發(fā)射天線向掘進面發(fā)射電磁波信號,信號在電性參數(shù)不同的物質(zhì)分界面上發(fā)生反射和透射現(xiàn)象,接收天線接收反射信號,通過分析接收信號的時間、相位等信息對異質(zhì)區(qū)的位置等特征進行解釋。

    圖1 空間步進頻雷達井下放置示意圖

    電磁波信號在介質(zhì)中傳播的快慢,影響回波響應出現(xiàn)的時間,根據(jù)相鄰回波響應之間的時間差,可計算電磁波傳播的距離。雷達發(fā)射信號在介質(zhì)中的傳播速度可根據(jù)式(1)求得:

    式中:為電磁波在自由空間傳播的速度;′為介質(zhì)的相對介電常數(shù)??芍殡姵?shù)越大,信號在介質(zhì)中的傳播速度越慢。

    空間步進頻雷達實現(xiàn)空間上的頻率步進信號,使用等間隔的均勻線陣組成線陣的每個陣元發(fā)射頻率步進的信號??臻g步進頻雷達發(fā)射陣列結構如圖2所示,發(fā)射天線為個陣元的均勻線陣,其陣元間距為。第個陣元的發(fā)射信號x ()為:

    圖2 空間步進頻雷達發(fā)射陣列結構

    (1)空間位置管理需要采用矢量數(shù)據(jù)結構。在地理位置表述時,矢量數(shù)據(jù)具有精確度高、數(shù)據(jù)量小的特性,以方便全面系統(tǒng)的闡述空間位置關系,對于數(shù)據(jù)圖形、屬性及數(shù)據(jù)恢復均能在短時間內(nèi)實現(xiàn),其具體格式如圖4。

    2 FDTD時域有限差分原理

    井下探地雷達探測異質(zhì)區(qū)屬于電磁類探測方法。FDTD十分適用于探地雷達建模,因為它能有效地模擬具有不同幾何形狀和介電特性的地下分層結構。FDTD方法將微分形式的時域麥克斯韋方程轉(zhuǎn)化為差分方程,使用迭代的方法求解有限空間、時間范圍內(nèi)的電場和磁場分量。一般情況下,對于非均勻、各項同性、線性的媒質(zhì)材料,其麥克斯韋電磁方程組微分形式如下所示:

    式中:表示電場強度矢量,單位為V/m;表示磁場強度矢量,單位為A/m;表示電位移矢量,單位為C/m;表示磁感應強度,單位為T;表示電流密度矢量,單位為A/m;J 表示磁通量密度矢量,單位為V/m;表示電荷密度,單位為C/m。

    同時有電磁場對材料介質(zhì)作用的關系式,即本構方程:

    式中:表示物質(zhì)材料的介電常數(shù),單位為F/m;表示物質(zhì)材料的電導率,單位為S/m;表示物質(zhì)材料的磁導率,單位為H/m;σ表示計算磁損耗的磁阻率,單位為Ω/m。

    FDTD方法對計算空間進行網(wǎng)格的剖分,圖3為著名的Yee元胞,元胞在坐標軸方向上的大小稱為該方向上的空間步長,規(guī)定每一個(或)場分量周圍有4個(或)場分量環(huán)繞,應用這種離散方式將含時間變量的麥克斯韋旋度方程轉(zhuǎn)化為差分方程。

    圖3 FDTD離散中的元胞

    對于二維問題,所有物理量與無關,??=0。以TM波為例,其麥克斯韋旋度方程為式(8),式(9)~式(11)為推導所得任意節(jié)點處電場和磁場的FDTD差分公式,式中分別為(,+1 2),(+1 2,),(,)。

    由于使用差分方法對麥克斯韋方程進行數(shù)值計算會產(chǎn)生色散現(xiàn)象,即波的傳播速度會隨波長而改變。這種色散會引起波形的畸變,且網(wǎng)格越粗色散越嚴重。為了使差分方程的計算結果能更近似于原方程的解,需要對時間步長和空間步長進行一定的限制,在二維問題中應該滿足式(12)。

    3 正演模擬

    3.1 實驗模型及參數(shù)設置

    在煤礦開采的過程中,地質(zhì)條件復雜多變,將存在于煤體中的非煤區(qū)域稱為異質(zhì)區(qū)。煤層中常見的異質(zhì)區(qū)有陷落柱、斷層、采空區(qū)、富水區(qū)、瓦斯富集區(qū)等。本文對采空區(qū)、富水區(qū)以及巖層三種主要煤礦異質(zhì)區(qū)進行模擬。未考慮巷道中各種儀器、線纜等產(chǎn)生的干擾,并假設介質(zhì)為均勻同性,建立掘進工作方向煤層內(nèi)部的物理模型如圖4所示。模型設置異質(zhì)區(qū)上方煤層厚度為0.5 m,采空區(qū)厚度為0.5 m,富水區(qū)厚度為0.5 m,薄泥巖層厚度為0.5 m。

    圖4 異質(zhì)區(qū)物理模型

    本文使用收發(fā)分置的天線陣,收發(fā)天線均為5陣元的均勻線陣,陣元間隔為0.1 m。收發(fā)天線陣水平放置,與掘進方向煤層表面平行。收發(fā)天線位于同一水平面,發(fā)射天線陣放置距煤層表面0.8 m,接收天線陣放置距煤層表面0.6 m。發(fā)射天線發(fā)射時長為1 ns的正弦信號,選定線陣一端的陣元作為第一陣元,其所發(fā)射信號的中心頻率為200 MHz,其余陣元的中心頻率依次以200 MHz的頻率增量進行步進,這樣第5個陣元發(fā)射信號的中心頻率為1 GHz。

    模擬中的最高頻率選擇為陣元中最高中心頻率的3倍,由最高頻率確定的最小波長為0.1 m,計算得到空間步長不應超過0.01 m,本文實驗中3個二維異質(zhì)區(qū)模型在方向與方向上的空間網(wǎng)格步長均設定為0.003 m。由式(12)可算得時間步長為7.076 0×10s。本文使用完美匹配層(Perfectly Matched Layer,PML)作為吸收邊界條件對空間進行截斷。異質(zhì)區(qū)模型中所用介質(zhì)的電性參數(shù)如表1所示。

    表1 模型電性參數(shù)

    3.2 傳播特性分析

    不改變發(fā)射天線以及接收天線的位置,對上述三種異質(zhì)區(qū)進行數(shù)值計算,得到接收天線5個陣元處隨時間迭代次數(shù)變化的電場信號數(shù)值。將5個接收點結果進行繪圖顯示,得到三種異質(zhì)區(qū)的仿真結果如圖5~圖7所示。

    圖5 采空區(qū)結果

    圖7 薄泥巖層結果

    圖5為采空區(qū)模擬結果,煤與空氣界面表現(xiàn)為一條明顯的軸線,出現(xiàn)在第1 157次時間迭代處,通過計算可知該點時間為8.187 ns。采空區(qū)上下表面回波對應的反射軸線分別出現(xiàn)在第2 337次與第2 789次時間迭代處,計算得到對應時間分別為16.537 ns和19.735 ns。某位置的時間等于該位置的時間迭代次數(shù)乘以時間步進間隔,時間步進間隔可由式(12)求得。

    圖6為富水區(qū)模型的模擬結果,可以清楚地看到在煤與空氣界面形成的響應后僅有一條清楚的響應軸線,該反射軸線表示煤層與富水區(qū)上表面的回波。由于水的介電常數(shù)大,電磁波在水中的衰減速度快。當信號穿過煤層進入富水區(qū)之后很快被吸收,因此僅能看到煤與水上表面形成的回波響應。

    圖6 富水區(qū)結果

    圖7為薄泥巖層模型的模擬結果,與采空區(qū)模型相似,同樣可以看到煤與泥巖上下交界面處的兩個回波,但與采空區(qū)相比較為微弱。

    實驗設置天線陣平行于掘進煤層表面,即發(fā)射信號為垂直入射,根據(jù)雙程走時′和速度可計算煤體到異質(zhì)區(qū)上表面厚度=′2,電磁波在煤中的傳播速度利用式(1)求得。表2為根據(jù)實驗數(shù)據(jù)對模型厚度進行計算的結果,表3為對模型內(nèi)介質(zhì)速度進行計算的結果。速度從大到小排序為空氣、煤、泥巖、水。根據(jù)仿真數(shù)據(jù)得到的計算結果與理論值存在一定的誤差,厚度的誤差不超過0.05 m,速度的誤差不超過8%。除仿真過程中產(chǎn)生的誤差外,還有結果讀取和計算時不精確造成的誤差。

    表2 模型厚度分析 m

    表3 模型介質(zhì)速度分析

    異質(zhì)區(qū)的物理模型不變,改變異質(zhì)區(qū)上方煤層厚度為5 m,采空區(qū)、富水區(qū)、泥巖區(qū)三種異質(zhì)區(qū)的厚度均為5 m,和方向的空間步長調(diào)整為0.01 m,其余參數(shù)不變,得到的結果如圖8~圖10所示。同樣為采空區(qū)與泥巖區(qū)表現(xiàn)為兩條響應軸線,泥巖區(qū)的響應較弱,富水區(qū)表現(xiàn)為一條響應軸線,與厚度為0.5 m模型下的模擬結果一致,僅到達時間不同。

    圖8 5 m厚采空區(qū)模擬結果

    圖9 5 m厚富水區(qū)模擬結果

    圖10 5 m厚泥巖區(qū)模擬結果

    為了更好地觀察空間步進頻信號在整個仿真區(qū)域的傳播情況,計算了三種模型空間中所有網(wǎng)格的電場強度,繪制頻率為1 GHz情況下二維仿真模型平面上的電場分布圖,如圖11所示。

    圖11 三種模型的電場分布圖

    電場分布圖中的橫縱坐標分別為模型掘進方向與垂直掘進方向上的距離坐標。設置=0的界面為煤礦開采工作面,橫坐標的正半軸表示掘進面上方的巷道;橫坐標的負半軸表示掘進面下方的煤層內(nèi)部。

    右側(cè)各發(fā)射陣元發(fā)射頻率步進的信號向左側(cè)煤層中進行傳播。圖11a)為采空區(qū)模型下的電場分布情況,電磁波信號由自由空間傳入煤層中,采空區(qū)部分波紋與煤層上表面自由空間處波紋一致;圖11b)為富水區(qū)的電場分布情況,由于能量衰減快在傳播中消耗殆盡,電磁波基本不能穿過水層;圖11c)為薄泥巖層的電場分布情況,電磁波從煤層到泥巖中波紋變密。電磁波信號傳播過程中波紋有密有疏,波紋的稀疏情況隨波長大小發(fā)生變化。在頻率不變的情況下,波速越大波長就越大,波紋越稀疏,根據(jù)電場分布圖得到模型中介質(zhì)速度由大到小排序分別為空氣、煤、泥巖、水,與理論一致。

    4 結 語

    本文基于時域有限差分方法對空間步進頻信號在三種不同井下異質(zhì)區(qū)模型中的傳播情況進行正演模擬,模擬結果表明:通過使用空間步進頻信號的時域波形進行數(shù)值計算,得到的煤層厚度與正演模型設置的參數(shù)基本吻合;在三種異質(zhì)區(qū)正演結果中,采空區(qū)模型與薄泥巖層模型下均可看到信號在異質(zhì)區(qū)上下表面產(chǎn)生的兩條反射軸線,但采空區(qū)的回波更為明顯,而富水區(qū)模型下僅可看到水煤上交界面的反射軸線;根據(jù)模擬的電場分布圖,通過波紋疏密程度可以對介質(zhì)的速度大小進行大致的排序;信號傳播時衰減的快慢不同,在水中的衰減最大,泥巖次之,空氣衰減最小。

    由以上仿真結果可知:進行實際探測時,若探測數(shù)據(jù)出現(xiàn)明顯的回波異常區(qū)域,則可以判斷該處存在與煤體不同的異質(zhì)區(qū),并可根據(jù)異?;夭ǔ霈F(xiàn)的時間點對異質(zhì)區(qū)域上方的煤層厚度進行估計;若信號在產(chǎn)生響應后急速衰減則可判斷該區(qū)域存在水,驗證了使用多發(fā)多收體制下的空間步進頻雷達對煤體異質(zhì)區(qū)探測的可行性。

    猜你喜歡
    探地富水異質(zhì)
    基于超表面的探地雷達增強探測研究
    全極化探地雷達系統(tǒng)
    富水砂卵石地層錨索施工工藝工法
    超級高密度電法在新疆某礦區(qū)富水性評價應用研究
    基于探地雷達法的地下管線探測頻譜分析
    隨機與異質(zhì)網(wǎng)絡共存的SIS傳染病模型的定性分析
    掘進巷道構造富水性電法探測綜合應用研究
    Ag2CO3/Ag2O異質(zhì)p-n結光催化劑的制備及其可見光光催化性能
    MoS2/ZnO異質(zhì)結的光電特性
    淺埋偏壓富水隧道掘進支護工藝分析
    国产精品国产av在线观看| av在线观看视频网站免费| 日本黄大片高清| 亚洲真实伦在线观看| 国产亚洲91精品色在线| 最后的刺客免费高清国语| 91精品国产九色| 小蜜桃在线观看免费完整版高清| av天堂中文字幕网| 国产真实伦视频高清在线观看| 九色成人免费人妻av| 国产免费视频播放在线视频| 久久精品久久久久久久性| 少妇精品久久久久久久| 美女脱内裤让男人舔精品视频| 亚洲精品,欧美精品| 久久人人爽人人爽人人片va| 女性被躁到高潮视频| 精品亚洲乱码少妇综合久久| .国产精品久久| 亚洲天堂av无毛| 激情 狠狠 欧美| 国产老妇伦熟女老妇高清| 尤物成人国产欧美一区二区三区| 国产精品不卡视频一区二区| 汤姆久久久久久久影院中文字幕| 91久久精品国产一区二区三区| 国内精品宾馆在线| 欧美97在线视频| 看十八女毛片水多多多| 亚洲美女视频黄频| 大片免费播放器 马上看| 嘟嘟电影网在线观看| 日韩视频在线欧美| 国产精品一区二区性色av| 国产精品一及| 好男人视频免费观看在线| 欧美日韩亚洲高清精品| 中文资源天堂在线| 亚洲精品日韩av片在线观看| 超碰av人人做人人爽久久| 99久国产av精品国产电影| 亚洲,一卡二卡三卡| 国产老妇伦熟女老妇高清| 国产 一区 欧美 日韩| 女性生殖器流出的白浆| 日韩精品有码人妻一区| 一二三四中文在线观看免费高清| 国产久久久一区二区三区| 亚洲欧美日韩东京热| 自拍偷自拍亚洲精品老妇| 伦理电影免费视频| 久久久久性生活片| 欧美+日韩+精品| 亚洲精品aⅴ在线观看| 18禁动态无遮挡网站| 99热国产这里只有精品6| 久久影院123| 在线免费观看不下载黄p国产| 久久精品熟女亚洲av麻豆精品| 男女免费视频国产| 亚洲精品乱久久久久久| 国产精品嫩草影院av在线观看| 精品亚洲成国产av| 两个人的视频大全免费| 国产亚洲最大av| 最近的中文字幕免费完整| 成人美女网站在线观看视频| 亚洲av日韩在线播放| xxx大片免费视频| 国产亚洲精品久久久com| 九九爱精品视频在线观看| 黄色配什么色好看| 成人国产麻豆网| 欧美成人a在线观看| 蜜桃久久精品国产亚洲av| 九草在线视频观看| 丰满人妻一区二区三区视频av| 午夜老司机福利剧场| 久久久久久久久久久丰满| av天堂中文字幕网| 一边亲一边摸免费视频| 人妻少妇偷人精品九色| 亚洲,欧美,日韩| 婷婷色麻豆天堂久久| 噜噜噜噜噜久久久久久91| 日韩一区二区三区影片| 久久久久国产网址| 中国三级夫妇交换| 精品久久久精品久久久| 久久这里有精品视频免费| 中文天堂在线官网| 久久人人爽av亚洲精品天堂 | 久久av网站| 只有这里有精品99| 亚洲婷婷狠狠爱综合网| 一个人免费看片子| 亚洲一区二区三区欧美精品| 91在线精品国自产拍蜜月| 欧美极品一区二区三区四区| 你懂的网址亚洲精品在线观看| 97超视频在线观看视频| 51国产日韩欧美| 国产高潮美女av| 欧美精品人与动牲交sv欧美| 2018国产大陆天天弄谢| 日本黄色片子视频| 夫妻性生交免费视频一级片| 高清日韩中文字幕在线| 在现免费观看毛片| 在线播放无遮挡| 黄片无遮挡物在线观看| 嫩草影院新地址| 国产欧美日韩一区二区三区在线 | 精品一区在线观看国产| 国产一级毛片在线| 免费av不卡在线播放| 中国三级夫妇交换| 麻豆精品久久久久久蜜桃| 亚洲久久久国产精品| 婷婷色综合大香蕉| 一本一本综合久久| 高清视频免费观看一区二区| 观看av在线不卡| 啦啦啦视频在线资源免费观看| 国产久久久一区二区三区| 亚洲国产高清在线一区二区三| xxx大片免费视频| 狂野欧美激情性bbbbbb| 五月开心婷婷网| 99热网站在线观看| 少妇 在线观看| 欧美+日韩+精品| 日韩精品有码人妻一区| 下体分泌物呈黄色| 少妇人妻 视频| 国语对白做爰xxxⅹ性视频网站| 在线观看一区二区三区激情| 欧美三级亚洲精品| 国产综合精华液| 亚洲不卡免费看| 全区人妻精品视频| 多毛熟女@视频| 啦啦啦啦在线视频资源| 一级毛片黄色毛片免费观看视频| 国产国拍精品亚洲av在线观看| 国产美女午夜福利| 免费大片黄手机在线观看| 精品久久久久久电影网| 少妇人妻精品综合一区二区| 高清视频免费观看一区二区| 各种免费的搞黄视频| 欧美精品亚洲一区二区| 国产探花极品一区二区| 黑人高潮一二区| 最近的中文字幕免费完整| 国产成人freesex在线| 精品久久久久久电影网| 99热这里只有是精品在线观看| 少妇高潮的动态图| 亚洲欧洲日产国产| 欧美成人午夜免费资源| freevideosex欧美| 最近手机中文字幕大全| 久久久久久久久久成人| 欧美精品国产亚洲| 中文字幕精品免费在线观看视频 | 成人黄色视频免费在线看| 18禁在线无遮挡免费观看视频| 亚洲第一av免费看| 99热国产这里只有精品6| 亚洲四区av| 欧美老熟妇乱子伦牲交| 成人国产麻豆网| 老司机影院成人| 亚洲无线观看免费| 美女国产视频在线观看| 国产精品99久久99久久久不卡 | 18+在线观看网站| 欧美bdsm另类| 狂野欧美激情性bbbbbb| 全区人妻精品视频| 一边亲一边摸免费视频| 老司机影院成人| 国产高清三级在线| 日本黄色日本黄色录像| 国产免费福利视频在线观看| 搡女人真爽免费视频火全软件| 亚洲av福利一区| 成人漫画全彩无遮挡| 男人和女人高潮做爰伦理| 国产爽快片一区二区三区| 中国美白少妇内射xxxbb| 亚洲av日韩在线播放| 亚洲欧洲国产日韩| 啦啦啦视频在线资源免费观看| 乱系列少妇在线播放| 乱码一卡2卡4卡精品| 亚洲欧美日韩东京热| 只有这里有精品99| 丰满乱子伦码专区| 精品一区二区三区视频在线| 高清在线视频一区二区三区| 女性生殖器流出的白浆| 国产精品av视频在线免费观看| 日日摸夜夜添夜夜爱| 夜夜爽夜夜爽视频| 亚洲美女视频黄频| 久久精品久久久久久久性| 九九在线视频观看精品| 亚洲av电影在线观看一区二区三区| 九九爱精品视频在线观看| 天天躁夜夜躁狠狠久久av| 午夜福利视频精品| 91狼人影院| 爱豆传媒免费全集在线观看| 久久精品久久久久久久性| 国产熟女欧美一区二区| 久久97久久精品| 欧美97在线视频| 美女主播在线视频| 少妇人妻精品综合一区二区| 久久久久久久久久人人人人人人| 超碰av人人做人人爽久久| 一级爰片在线观看| 精品一区二区免费观看| 国产精品久久久久久精品古装| 熟女电影av网| 成人影院久久| 大又大粗又爽又黄少妇毛片口| 哪个播放器可以免费观看大片| 国产精品熟女久久久久浪| 国产精品99久久久久久久久| 欧美另类一区| 国产成人a区在线观看| 日本色播在线视频| 乱系列少妇在线播放| 三级经典国产精品| 人人妻人人添人人爽欧美一区卜 | 日韩不卡一区二区三区视频在线| 干丝袜人妻中文字幕| 日本-黄色视频高清免费观看| 综合色丁香网| 美女中出高潮动态图| 久久精品久久久久久噜噜老黄| 国产美女午夜福利| 日本黄大片高清| 免费观看性生交大片5| 国产淫语在线视频| 国内揄拍国产精品人妻在线| 欧美xxxx性猛交bbbb| 欧美日韩综合久久久久久| 国产高清国产精品国产三级 | 国产视频首页在线观看| 最近最新中文字幕大全电影3| 少妇人妻一区二区三区视频| 国产精品熟女久久久久浪| 少妇高潮的动态图| 青青草视频在线视频观看| 青春草视频在线免费观看| 日韩不卡一区二区三区视频在线| 99热全是精品| 午夜福利网站1000一区二区三区| 天天躁日日操中文字幕| 不卡视频在线观看欧美| 亚洲av电影在线观看一区二区三区| 少妇精品久久久久久久| 免费播放大片免费观看视频在线观看| 亚洲成人一二三区av| 国产女主播在线喷水免费视频网站| 天天躁日日操中文字幕| 国产成人精品福利久久| 午夜免费鲁丝| 中国美白少妇内射xxxbb| 午夜激情久久久久久久| 国产精品麻豆人妻色哟哟久久| 国产 一区精品| 亚洲精品日本国产第一区| 免费av不卡在线播放| 视频中文字幕在线观看| 在线播放无遮挡| 啦啦啦在线观看免费高清www| 欧美极品一区二区三区四区| 亚洲真实伦在线观看| 久久6这里有精品| 秋霞伦理黄片| 久久人人爽人人爽人人片va| 尤物成人国产欧美一区二区三区| 日本午夜av视频| 在线免费观看不下载黄p国产| 亚洲欧美精品自产自拍| 国产在视频线精品| 久久国产精品大桥未久av | 国产一级毛片在线| 国产精品无大码| 国产美女午夜福利| 欧美xxⅹ黑人| 偷拍熟女少妇极品色| 纯流量卡能插随身wifi吗| 国产精品一区二区三区四区免费观看| 日韩伦理黄色片| 久久午夜福利片| 国产乱人偷精品视频| 久久人妻熟女aⅴ| 免费播放大片免费观看视频在线观看| 午夜福利视频精品| 一个人免费看片子| 亚洲精品,欧美精品| 亚洲aⅴ乱码一区二区在线播放| av黄色大香蕉| 日韩欧美 国产精品| 激情五月婷婷亚洲| 狠狠精品人妻久久久久久综合| 全区人妻精品视频| 成人特级av手机在线观看| 亚洲精品中文字幕在线视频 | 成人综合一区亚洲| 汤姆久久久久久久影院中文字幕| 高清不卡的av网站| 丝袜喷水一区| 亚洲精品视频女| 高清不卡的av网站| 成人综合一区亚洲| 国内揄拍国产精品人妻在线| 免费看av在线观看网站| 国产真实伦视频高清在线观看| 亚洲精品成人av观看孕妇| 1000部很黄的大片| 国产片特级美女逼逼视频| 你懂的网址亚洲精品在线观看| 18禁在线无遮挡免费观看视频| 看免费成人av毛片| 高清不卡的av网站| 国产白丝娇喘喷水9色精品| 黄色怎么调成土黄色| 尾随美女入室| 偷拍熟女少妇极品色| 国产精品人妻久久久久久| 国产亚洲午夜精品一区二区久久| 秋霞伦理黄片| 91久久精品国产一区二区三区| 美女中出高潮动态图| 亚洲自偷自拍三级| 男女无遮挡免费网站观看| 在线观看三级黄色| 大又大粗又爽又黄少妇毛片口| 久久久久视频综合| 国产日韩欧美在线精品| 自拍欧美九色日韩亚洲蝌蚪91 | 久久人妻熟女aⅴ| 精品久久久精品久久久| 免费观看av网站的网址| av不卡在线播放| 亚洲成人av在线免费| 视频中文字幕在线观看| 久久av网站| 欧美性感艳星| 欧美+日韩+精品| 18禁在线播放成人免费| 国产探花极品一区二区| 欧美zozozo另类| 岛国毛片在线播放| 国产精品久久久久久av不卡| 久久鲁丝午夜福利片| 人人妻人人添人人爽欧美一区卜 | 国产精品欧美亚洲77777| 男女无遮挡免费网站观看| 午夜激情久久久久久久| 久久久精品免费免费高清| 国产欧美另类精品又又久久亚洲欧美| 美女中出高潮动态图| 国产国拍精品亚洲av在线观看| 中文字幕久久专区| 2018国产大陆天天弄谢| 国产成人午夜福利电影在线观看| 偷拍熟女少妇极品色| 91狼人影院| 超碰av人人做人人爽久久| 欧美+日韩+精品| 日韩成人av中文字幕在线观看| 人人妻人人爽人人添夜夜欢视频 | 一级毛片我不卡| 黑丝袜美女国产一区| 能在线免费看毛片的网站| 熟女电影av网| 亚洲国产色片| 亚洲av国产av综合av卡| 免费看av在线观看网站| 最新中文字幕久久久久| av黄色大香蕉| 又粗又硬又长又爽又黄的视频| 亚洲人成网站高清观看| 日本猛色少妇xxxxx猛交久久| 精品一品国产午夜福利视频| 久久久久人妻精品一区果冻| 伦精品一区二区三区| 久久久久久久精品精品| 亚洲经典国产精华液单| 日韩免费高清中文字幕av| 一级毛片黄色毛片免费观看视频| 大香蕉97超碰在线| 肉色欧美久久久久久久蜜桃| 亚洲av免费高清在线观看| 男女免费视频国产| 少妇人妻精品综合一区二区| 精品人妻视频免费看| 免费看日本二区| 日韩伦理黄色片| 99久久人妻综合| 久久久久久伊人网av| 老女人水多毛片| 亚州av有码| 精品一区二区三卡| 日韩强制内射视频| 亚洲久久久国产精品| 亚洲,一卡二卡三卡| 亚洲成色77777| 亚洲va在线va天堂va国产| 国产精品久久久久久久电影| 国产精品嫩草影院av在线观看| 国产爱豆传媒在线观看| 建设人人有责人人尽责人人享有的 | 午夜日本视频在线| 国产一级毛片在线| 美女国产视频在线观看| 一级毛片aaaaaa免费看小| 久久久精品94久久精品| 日本黄色日本黄色录像| 一区二区三区乱码不卡18| 男女免费视频国产| 久久精品国产亚洲网站| 99热国产这里只有精品6| 精品人妻一区二区三区麻豆| 大码成人一级视频| 亚洲精品一二三| 校园人妻丝袜中文字幕| 夫妻性生交免费视频一级片| 在现免费观看毛片| 亚洲激情五月婷婷啪啪| 国产精品蜜桃在线观看| 日韩大片免费观看网站| 久久这里有精品视频免费| 日韩av免费高清视频| 97在线人人人人妻| 国产精品秋霞免费鲁丝片| 欧美成人精品欧美一级黄| 视频区图区小说| 久久久久久久亚洲中文字幕| 黑人猛操日本美女一级片| 欧美人与善性xxx| 91久久精品国产一区二区成人| 看免费成人av毛片| 我要看黄色一级片免费的| 欧美精品人与动牲交sv欧美| 国产无遮挡羞羞视频在线观看| 国内揄拍国产精品人妻在线| 国产色婷婷99| 亚洲av电影在线观看一区二区三区| 欧美激情极品国产一区二区三区 | 男人舔奶头视频| 极品教师在线视频| 各种免费的搞黄视频| 国产在线男女| av在线播放精品| 高清午夜精品一区二区三区| 男女国产视频网站| 99热这里只有是精品50| 在线观看国产h片| 人妻制服诱惑在线中文字幕| 18禁裸乳无遮挡动漫免费视频| 直男gayav资源| 搡女人真爽免费视频火全软件| 国产 精品1| av黄色大香蕉| 蜜臀久久99精品久久宅男| 免费黄色在线免费观看| 精品熟女少妇av免费看| 一区二区三区乱码不卡18| 国产av精品麻豆| 国产成人免费无遮挡视频| 精品酒店卫生间| 91精品国产国语对白视频| 久热这里只有精品99| 久久99热这里只有精品18| 美女cb高潮喷水在线观看| 网址你懂的国产日韩在线| 国产精品女同一区二区软件| 一级毛片aaaaaa免费看小| 七月丁香在线播放| 国产精品久久久久久精品古装| 2022亚洲国产成人精品| a 毛片基地| 身体一侧抽搐| 日本一二三区视频观看| 久久国产精品大桥未久av | 国产精品三级大全| 美女中出高潮动态图| 国产一区亚洲一区在线观看| 亚洲国产精品国产精品| 男女边摸边吃奶| 毛片一级片免费看久久久久| 国产淫语在线视频| 99热国产这里只有精品6| 亚洲国产精品成人久久小说| 国产在线免费精品| 麻豆国产97在线/欧美| 狠狠精品人妻久久久久久综合| 交换朋友夫妻互换小说| 啦啦啦中文免费视频观看日本| 国产成人a区在线观看| 亚洲欧美日韩无卡精品| 日本欧美国产在线视频| 尾随美女入室| 一级a做视频免费观看| 精品视频人人做人人爽| 97在线视频观看| 久热这里只有精品99| a级毛片免费高清观看在线播放| a级一级毛片免费在线观看| 欧美少妇被猛烈插入视频| 国产精品国产av在线观看| 91久久精品国产一区二区三区| 国产精品人妻久久久久久| 久久影院123| 亚洲性久久影院| 黄色怎么调成土黄色| 九草在线视频观看| 18禁在线无遮挡免费观看视频| 亚洲成人av在线免费| 最近手机中文字幕大全| 午夜福利网站1000一区二区三区| 国产精品99久久99久久久不卡 | 欧美日韩国产mv在线观看视频 | 日日撸夜夜添| 观看免费一级毛片| 亚洲国产欧美在线一区| 精品视频人人做人人爽| 亚洲精品国产色婷婷电影| 中国美白少妇内射xxxbb| 六月丁香七月| 亚洲精品自拍成人| 国产在线一区二区三区精| 国产视频首页在线观看| 99精国产麻豆久久婷婷| 国产视频内射| 各种免费的搞黄视频| 欧美人与善性xxx| av不卡在线播放| 日本色播在线视频| 国产淫语在线视频| 久久韩国三级中文字幕| 亚洲精品自拍成人| 国产日韩欧美亚洲二区| 超碰av人人做人人爽久久| 亚洲第一区二区三区不卡| 妹子高潮喷水视频| 国产亚洲91精品色在线| 大香蕉久久网| 日韩 亚洲 欧美在线| 国产精品av视频在线免费观看| 少妇的逼水好多| 丰满少妇做爰视频| 自拍偷自拍亚洲精品老妇| tube8黄色片| 嫩草影院新地址| 欧美精品一区二区大全| 亚洲精品乱码久久久v下载方式| 久久久精品94久久精品| 成人亚洲欧美一区二区av| 一级黄片播放器| 成人国产av品久久久| 青春草国产在线视频| 亚洲精品乱久久久久久| 国产在线视频一区二区| 乱系列少妇在线播放| 久久久成人免费电影| 日韩在线高清观看一区二区三区| 美女cb高潮喷水在线观看| 成人黄色视频免费在线看| 联通29元200g的流量卡| 六月丁香七月| 国产精品成人在线| 老师上课跳d突然被开到最大视频| 天堂8中文在线网| 亚洲精品aⅴ在线观看| 国产深夜福利视频在线观看| 亚洲av电影在线观看一区二区三区| 久久人妻熟女aⅴ| 国产欧美日韩一区二区三区在线 | 日韩中文字幕视频在线看片 | 老司机影院成人| 免费播放大片免费观看视频在线观看| 人人妻人人看人人澡| 夜夜骑夜夜射夜夜干| 日本与韩国留学比较| 精品一区二区三卡| 嫩草影院入口| 黄片无遮挡物在线观看| 免费人成在线观看视频色| 好男人视频免费观看在线| 街头女战士在线观看网站| 国产日韩欧美在线精品| 国产成人a∨麻豆精品| 一区二区三区精品91| 夜夜爽夜夜爽视频| 大片免费播放器 马上看| 国产精品成人在线| 久久婷婷青草| 插阴视频在线观看视频| 亚洲成人中文字幕在线播放| 赤兔流量卡办理| 婷婷色综合www| 99热国产这里只有精品6| 亚洲欧美成人精品一区二区| 日本欧美国产在线视频|