亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于BOTDR的煤層底板突水溫度場監(jiān)測模擬研究

        2016-05-25 00:37:01張平松孫斌楊許時昂
        關(guān)鍵詞:突水傳感溫度場

        張平松,孫斌楊,許時昂

        (1.安徽理工大學 地球與環(huán)境學院,安徽 淮南 232001;2.全國煤炭行業(yè)礦井水害綜合防治工程研究中心,安徽 淮南 232001)

        基于BOTDR的煤層底板突水溫度場監(jiān)測模擬研究

        張平松1, 2,孫斌楊1,許時昂1

        (1.安徽理工大學 地球與環(huán)境學院,安徽 淮南 232001;2.全國煤炭行業(yè)礦井水害綜合防治工程研究中心,安徽 淮南 232001)

        針對采煤工作面底板突水過程中常會引起突水區(qū)域溫度場的異常響應,利用分布式光纖測試技術(shù)對溫度場變化捕捉敏感的特性進行測試。通過設計光纖底板溫度場測試傳感裝置,構(gòu)建測試地質(zhì)地球物理模型,對底板突水所引起的巖層溫度場的變化進行測試模擬。結(jié)果表明:分布式光纖測試技術(shù)能很好分辨底板水導升過程中水溫改變所引起的溫度場變化特征,且對突水區(qū)域范圍判斷快速、有效。

        隧道工程;巖土工程;底板突水;溫度場;分布式光纖傳感;相似模擬

        0 引 言

        基于我國能源結(jié)構(gòu)的獨特性,新形勢下煤炭資源在整體能源消費中依然發(fā)揮主導作用。隨著礦井機械化生產(chǎn)的推進,煤炭資源的大規(guī)模開采與深度化發(fā)展使得礦井地質(zhì)環(huán)境發(fā)生較大的變化。此時受地質(zhì)條件的改變,礦井災害也呈現(xiàn)出相應的復雜性與多樣性,如礦井水害、煤與瓦斯突出、巖爆、高地溫等。這給礦井的深部化開采帶來較大的威脅,其中采煤工作面底板突水危害最為嚴重,使得華北巖溶地區(qū)大部分礦井的深部煤層難以正常開采[1],造成巨大的經(jīng)濟損失和自然資源的浪費。

        如何實現(xiàn)對底板突水的快速、有效探測,一直是眾多學者探究的重點和難點。目前在底板突水防治工作中主要采用理論分析和原位測試兩種方法。其中理論方法有“底板相對隔水層”、“突水系數(shù)法”[2]、“下三帶”[3]以及 “關(guān)鍵層”等理論[4-6]。這些理論研究對于底板突水的防治具有一定的指導作用,但是由于理論多基于較為理想化的狀態(tài),其與實際工程會有較大的差異。因此目前對于底板水的防治還多依據(jù)原位測試技術(shù)獲得相應的技術(shù)參數(shù)。其中常用的測試方法包括高密度電阻率法、三維并行電法、雙巷音頻電透視及礦井瞬變電磁法等[7-8]。這幾種方法對于礦井水害地質(zhì)條件及異常區(qū)的判斷具有良好效果,但在工作面的開采過程控制與監(jiān)測方面仍存在不足,其觀測系統(tǒng)的布置及現(xiàn)場操作較為繁瑣難以靈活應用。分布式光纖傳感技術(shù)對光纖各點的溫度、應力變化可以進行良好探測,已廣泛應用于井下管道、隧道、大壩水位以及倉庫火情的監(jiān)測[9-11],并取得良好成效。綜合考慮光纖本身具有安全性能好、抗電磁強、測量范圍廣、重量輕等優(yōu)點,同時遠距離傳輸其光損較小可以保證數(shù)據(jù)的可靠性,便于礦井實時和長期監(jiān)測。目前,分布式光纖測溫技術(shù)對于底板突水的溫度場測試研究還比較薄弱,因此其具有較大的應用實踐空間。

        采煤工作面底板突水的主要因素包括采動礦壓和水壓兩方面,大致可分為3個階段。為有效預防突水需進行實時監(jiān)測,主要對突水前期即導水裂隙初步形成和斷裂帶塑性擴展階段進行監(jiān)測和預警[12]。筆者以淮南某礦區(qū)工作面底板原始地質(zhì)資料為基礎,構(gòu)建一套底板突水相似模擬裝置,結(jié)合BOTDR分布式光纖傳感技術(shù)對底板突水進行溫度場監(jiān)測,研究了光纖在煤層底板巖層中的布設及安裝工藝,對突水期間煤層底板溫度場的分布及變化規(guī)律進行探究。

        1 布里淵光時域反射技術(shù)

        布里淵光時域反射技術(shù)(BOTDR)是分布式光纖傳感技術(shù)的一種[13]。通過檢測光纖中布里淵散射頻移的變化,根據(jù)變化量隨溫度和應變呈線性變化的關(guān)系來分析探測體物理性質(zhì)的改變。其工作原理如圖1[14-17]。

        圖1 BOTDR工作原理Fig.1 Working principle of BOTDR

        圖1中,利用測試儀在光纖的一端輸入脈沖光后,其內(nèi)部介質(zhì)之間的相互作用會發(fā)生布里淵背向散射,當光纖的某一點或某一區(qū)域的原始狀態(tài)發(fā)生改變時,原始光路也會隨之改變即布里淵散射的頻率發(fā)生變化,且變化量只與溫度和軸向應變有關(guān)。

        利用布里淵頻移(υB),光纖溫度(T)和軸向應變(ε)的關(guān)系,通過測量ΔυB進而得到光纖沿線各點T和ε的分布。T與υB之間的數(shù)值關(guān)系可用式(1)表達為:

        T=CT(υB-υB0)+T0

        (1)

        式中:T為光纖溫度;CT為常數(shù),即頻移與溫度的相關(guān)系數(shù);υB為布里淵頻移;υB0和T0分別為自然狀態(tài)下初始布里淵頻移量和溫度值。

        實驗所采用測試儀器為AV6419光纖應變分布測試儀,其量程范圍0.5~80 km,空間分辨率1 m,采樣間隔最小可達5 cm,對于溫度檢測具有0.5 ℃的分辨測試能力。底板突水相似模擬實驗采用基于BOTDR分布式光纖傳感技術(shù),可以對煤層底板溫度場進行實時監(jiān)測,當發(fā)生突水時測試儀可有效探測破壞區(qū)域和溫度場的變化,以及水位上升溫度場分布的演化過程,滿足工程實踐判斷的基本要求。

        2 底板突水溫度場監(jiān)測相似模擬

        我國華北石炭系煤層多形成于濱海、湖泊等海相沉積環(huán)境,含煤巖系中存在大量的灰?guī)r層,致使其基底巖溶發(fā)育,賦存大量的承壓水,給深部石炭系煤層的開采帶來極大的危害。筆者主要以兩淮煤田淮南某礦工作面底板為研究對象,設計和構(gòu)建煤層底板突水相似模擬實驗[18]。

        2.1 實驗模型設計

        底板突水相似模擬實驗根據(jù)礦井相關(guān)資料,經(jīng)過合理的簡化,以1∶100比例搭建長、寬、高均為0.4 m的平面應力模型,如圖2。構(gòu)建模型的圍壓較低,但是由于測試中水頭相對較高,因此并不影響模擬突水事件的發(fā)生及后期底板溫度場的監(jiān)測。

        圖2 相似模擬裝置Fig.2 Equipment of similar simulation

        結(jié)合模擬材料配比實驗研究成果,確定模型材料的配比方案,并利用式(2)對各分層材料用量進行計算(表1)。模型中各巖性分層采用砂子、石膏、石灰、水泥等分別以不同比例加入適量的水分層鋪設并壓實固結(jié)而成,搭載整個實驗模型共需要砂子71.25 kg,石灰4.61 kg,水泥2 kg,石膏3.77 kg。

        G=lbhγm

        (2)

        式中:G為模型分層材料總重(不包括含水量);l為模型長度;b為模型寬度;h為模擬巖層各分層高度;γm為模擬巖層各分層容重。

        表1 相似模擬材料配比方案

        2.2 光纖的布設及安裝

        模擬實驗傳感光纖為NZS-DSS-C07聚氨酯緊套光纜,其為單模光纖,光纜截面尺寸為2.0 mm,重量為2 kg/km。為提高AV6419光纖應變分布測試儀的空間分辨率,將傳感光纜以螺旋式的形式纏繞在PVC管外壁,共纏繞135圈加上光纖跳線總長為15 m(圖3),期間采用黏結(jié)劑使之與PVC管外壁更好地耦合。

        圖3 傳感光纜布設、安裝示意Fig.3 Layout and installation of sensing optical fiber cable

        在鋪設模擬材料之前,將傳感光纜連同PVC管傾斜放置于平面模型中,然后分層鋪設模擬材料。材料填充完畢后,在模擬煤層的上表面均勻施加外力使其搗固成型,介質(zhì)、光纜、PVC管這3者之間的耦合性將大大提高。本次模擬實驗光纜安裝工藝與現(xiàn)場施工有所不同,礦井實測時會在底板施工鉆孔,然后利用鉆孔安裝技術(shù)將光纜植入到孔中,后期進行注漿耦合。但是由于模擬試驗尺度較小無法達到施工鉆孔要求,只能前期預埋設光纜再鋪設巖層使之其與實際工況基本相符。

        2.3 模擬實驗過程

        實驗采用分層加水的方式,水溫總體控制在55 ℃左右,每次注水高度為2.5 cm。期間為了使得水流能夠從底板巖層均勻地上升,人工施加一定的水壓,將輸水軟管抬至2.5 m的高度,相當于0.025 MPa的水頭。模型搭載完畢,根據(jù)現(xiàn)場實驗條件設置儀器參數(shù)。首先在未注水條件下采集巖體空間背景溫度值,然后按注水間隔2.5 cm高度采集一次數(shù)據(jù),共擬采集16組數(shù)據(jù)。由于突水的突發(fā)性及難以控制性,在第7次注水時水流已完全覆蓋整個煤層工作面,實驗共采集10組數(shù)據(jù)。

        3 實驗數(shù)據(jù)分析

        在BOTDR分布式光纖傳感技術(shù)進行煤層工作面底板突水相似模擬實驗過程中,根據(jù)光纖的布里淵頻移與溫度之間存在相關(guān)關(guān)系[19-20],由式(1)可計算出不同布里淵頻移所對應的巖體溫度值,從而可獲取整個突水過程中底板巖層內(nèi)部溫度場的變化及分布情況。

        3.1 突水過程分析

        根據(jù)監(jiān)測數(shù)據(jù)獲得突水過程中的水位上升曲線,如圖4;突水時光纖各段溫度分布,如圖5。

        圖4 水位上升曲線Fig.4 Curve of water level rises

        圖5 突水時光纖各段溫度分布Fig.5 Temperature distribution in the fiber of the water inrush

        圖4中,突水界面將煤層底板巖石分為3個層次。其一為上部未受突水影響原始狀態(tài)的巖層(圖4中Ⅰ區(qū)域);其二為每次注水結(jié)束后穩(wěn)定狀態(tài)下的水源界面,判別原因是灰?guī)r承壓水沿著裂隙帶導升時,周圍巖石介質(zhì)會由于水溫高于自身的溫度而發(fā)生溫度突變,使得高于上部原始狀態(tài)的巖石溫度(圖4中Ⅱ區(qū)域);其三為已被承壓水貫穿的下部巖層,二者充分接觸使得巖體溫度明顯高于上部各個巖層的值(圖4中Ⅲ區(qū)域)。

        圖4中的Ⅱ區(qū)域水位上升曲線主要是由7次注水后的突水界面拼接而成,較為真實地反映出整個動態(tài)突水過程。同時,圖4和圖5中第7次突水都出現(xiàn)了陡增的現(xiàn)象,原因在于此時底板巖層中已富含大量的承壓水,并且?guī)r性為砂、泥巖,其孔隙度較大于底部的灰?guī)r層,因此以相同的水量和水壓注入底板時出現(xiàn)上述情形。

        3.2 溫度場變化及定位分析

        在實驗過程中,連同PVC管一起植入巖體內(nèi)部的光纜會由于溫度場的變化受影響,而裸露在外部的則不會發(fā)生布里淵頻移。實驗結(jié)果表明:BOTDR對溫度場的變化有良好的探測效果(圖5),溫度最終測量值基本都處于50~60 ℃之間,與實際注入的55 ℃水源在相對誤差允許范圍內(nèi)(55±5.5)℃基本吻合。

        分布式光纖傳感技術(shù)不但能夠靈敏地探測溫度場變化,而且可以有效定位破壞區(qū)域。由圖5可知,每次突水曲線的首個抬高點即突水區(qū)域零界點都不一樣,正符合于分層注水情況。選取第4次突水曲線做定量分析可得,突水后水源穩(wěn)定界面大致位于-32 cm,與實際注水高度10 cm(注水界面達到-30 cm)基本相符。

        綜上,基于BOTDR的分布式光纖傳感技術(shù)應用于采煤工作面底板突水的監(jiān)測,可以直觀地反映底板溫度場狀態(tài)并準確地定位突水水源位置,可對煤層工作面底板突水的預防提供技術(shù)參數(shù)。

        4 結(jié) 語

        鑒于BOTDR的先進性,筆者研究了用于監(jiān)測煤層工作面底板突水的傳感光纜布設方式及安裝工藝,并開展相關(guān)模擬實驗探究BOTDR分布式光纖傳感技術(shù)應用于底板突水監(jiān)測的有效性。

        研究結(jié)果表明:分布式光纖測試方法能夠有效分辨模型突水過程中溫度場的相關(guān)參數(shù)變化,通過監(jiān)測數(shù)據(jù)可有效圈定溫度場變化區(qū)域及溫度變化范圍。通過對比測試結(jié)果,光纖測試數(shù)據(jù)與實驗實際結(jié)果在誤差允許范圍內(nèi)吻合度較高,同時底板光纖測試施工工藝簡單、易操作,其測試方法和結(jié)果能為底板突水防治提供有效的技術(shù)指導。與傳統(tǒng)測試方法相比,基于BOTDR的分布式光纖傳感具有一定的優(yōu)越性,其應用前景十分廣闊,但具體的實踐應用需進一步深入研究。

        [1] 胡耀青.我國帶壓開采研究現(xiàn)狀及其展望[J].太原理工大學學報,2002,33(6):594-600. HU Yaoqing. The development and application situation on aquifer in China[J].JournalofTaiyuanUniversityofTechnology, 2002, 33(6): 594-600.

        [2] 施龍青.底板突水機理研究綜述[J].山東科技大學學報(自然科學版),2009,29(3):17-22. SHI Longqing. Summary of research on mechanism of water-inrush from seam floor[J].JournalofShandongUniversityofScienceandTechnology(NaturalScience), 2009, 29(3): 17-22.

        [3] 荊自剛,李白英.煤層底板突水機理的初步探討[J].煤田地質(zhì)與勘探,1980(2):27-29. JING Zigang, LI Baiying. Preliminary study about water-inrush from floor[J].CoalfieldGeologyandExploration, 1980(2): 27-29.

        [4] 錢鳴高,繆協(xié)興,許家林.巖層控制中的關(guān)鍵層理論研究[J].煤炭學報,1996,21(3):44-47. QIAN Minggao, MIAO Xiexing, XU Jialin. Theoretical study of key stratum in ground control[J].JournalofChinaCoalSociety, 1996, 21(3): 44-47.

        [5] 王作宇,張建華,劉鴻泉,等.承壓水上近距煤層重復采動的底板巖體移動規(guī)律[J].煤炭科學技術(shù),1995,23(2):9-12. WANG Zuoyu, ZHANG Jianhua, LIU Hongquan, et al. Strata movement law of repeat mining close coal layer[J].CoalScienceandTechnology, 1995, 23(2): 9-12.

        [6] 施龍青,韓進.開采煤層底板“四帶”劃分理論與實踐[J].中國礦業(yè)大學學報,2005,34(1):16-23. SHI Longqing, HAN Jin. Theory and practice of dividing coal mining area floor into four-zone[J].JournalofChinaUniversityofMining&Technology, 2005, 34(1): 16-23.

        [7] 付茂如,張平松,王大設,等.礦井工作面底板水害綜合探查技術(shù)研究[J].礦業(yè)安全與環(huán)保,2013,40(2):92-95. FU Maoru, ZHANG Pinsong, WANG Dashe, et al. Study on integrated detection for face floor water hazards[J].MiningSafety&EnvironmentalProtection, 2013, 40(2): 92-95.

        [8] 吳榮新,劉盛東,張平松.雙巷并行三維電法探測煤層工作面底板富水區(qū)[J].煤炭學報,2010,35(3):454-457. WU Rongxin, LIU Shengdong, ZHANG Pingsong. The exploration of two-gateways parallel 3-D electrical technology for water-rich area within coal face floor[J].JournalofChinaCoalSociety, 2010, 35(3): 454-457.

        [9] 蔣奇.分布式光纖溫度傳感技術(shù)在隧道監(jiān)測中的應用[J].應用光學,2005,26(3):20-22. JIANG Qi. Application of fiber-optic distributed temperature sensor to tunnel monitoring system[J].JournalofAppliedOptics, 2005, 26(3): 20-22.

        [10] 何寧,丁勇,吳玉龍,等.基于分布式光纖測溫技術(shù)的堤壩滲漏監(jiān)測[J].水利水運工程學報,2015(1):20-27. HE Ning, DING Yong, WU Yulong, et al. Experimental study of distributed optical fiber temperature measurement technology for measuring leakage in embankment dam[J].Hydro-ScienceandEngineering, 2015(1): 20-27.

        [11] 程宗澤,李奔鋒.分布式光纖感溫報警系統(tǒng)在礦井火災的監(jiān)測研究[J].煤礦機電,2013(5):45-48. CHENG Zongze, LI Benfeng. Research on distributed fiber temperature sensing alarm system for monitoring of ming fire[J].CollieryMechanical&ElectricalTechnology, 2013(5): 45-48.

        [12] 張長文,付斌,徐毅.礦井突水問題的研究[J].煤炭技術(shù),2004,23(5):65-66. ZHANG Changwen, FU Bin, XU Yi. Research on gushing water in coal-mine[J].CoalTechnology, 2004, 23(5): 65-66.

        [13] 周柏兵,萬永波,徐國龍,等.分布式光纖監(jiān)測系統(tǒng)在混凝土壩的研究與應用[J].水利信息化,2010(6):32-36. ZHOU Baibing, WAN Yongbo, XU Guolong, et al. Research in and application of distributed optical fiber monitoring system at concrete dam[J].WaterResourcesInformatization, 2010(6): 32-36.

        [14] 劉杰,施斌,張丹,等.基于BOTDR的基坑變形分布式監(jiān)測實驗研究[J].巖土力學,2006,27(7):1224-1228. LIU Jie, SHI Bin, ZHANG Dan, et al. Experimental study of foundation pit monitoring using BOTDR-based on distributed optical fiber sensor[J].RockandSoilMechanics, 2006,27(7):1224-1228.

        [15] 施斌,徐洪鐘,張丹,等.BOTDR應變監(jiān)測技術(shù)應用在大型基礎工程健康診斷中的可行性研究[J].巖石力學與工程學報,2004,23(3):493-499. SHI Bin, XU Hongzhong, ZHANG Dan, et al. Feasibility study on application of BOTDR to health monitoring for large infrastructure engineering[J].ChineseJournalofRockMechanicsandEngineering, 2004, 23(3): 493-499.

        [16] 隋海波,施斌,張丹,等.基于BOTDR的錨桿拉拔試驗研究[J].巖土工程學報,2008,30(5):755-759. SUI Haibo, SHI Bin, ZHANG Dan, et al. BOTDR-based pull-out tests on anchor bolts[J].ChineseJournalofGeotechnicalEngineering, 2008, 30(5): 755-759.

        [17] 黃民雙,陳偉民,黃尚廉,等.基于Brillouin散射的分布式光纖拉伸應變傳感器的理論分析[J].光電工程,1995,22(4):11-16. HUANG Minshuang, CHEN Weimin, HUANG Shanglian, et al. The theoretical analysis of distributed optic fiber strain sensor based on Brillouin scatter[J].Opto-ElectronicEngineering, 1995, 22(4): 11-16.

        [18] 高峰,胡蓉,譚緒凱.隧道注漿加固模型試驗研究[J].重慶交通大學學報(自然科學版),2014,33(4):44-46. GAO Feng, HU Rong, TAN Xukai. Model test of tunnel grouting reinforcement[J].JournalofChongqingJiaotongUniversity(NaturalScience), 2014, 33(4): 44-46.

        [19] 張丹,張平松,施斌,等.采場覆巖變形與破壞的分布式光纖監(jiān)測與分析[J].巖土工程學報,2015,37(5):952-957. ZHANG Dan, ZHANG Pingsong, SHI Bin, et al. Monitoring and analysis of overburden deformation and failure using distributed fiber optic sensing[J].ChineseJournalofGeotechnicalEngineering, 2015, 37(5): 952-957.

        [20] XU Shiang, ZHNAG Pingsong, ZHANG Dan, et al. Simulation study of fiber optic monitoring technology of surrounding rock deformation under deep mining conditions[J].JournalofCivilStructuralHealthMonitoring, 2015, 5(5): 563-571.

        Simulation Research on BOTDR-based Monitoring over Temperature Field of Water Inrushing from Coal Floor

        ZHANG Pingsong1, 2, SUN Binyang1, XU Shi’ang1

        (1. School of Earth and Environment, Anhui University of Science and Technology, Huainan 232001, Anhui, P.R.China; 2. Engineering Research Center of Mine Water Disaster Comprehensive Control of Nationwide Mining Industry, Huai’nan 232001, Anhui, P.R.China)

        During the process of floor water inrushing on coal face, abnormal responses of corresponding temperature field occured. Therefore, the sensitiveness of distributed fiber testing technology was used to test the changed temperature field. Through designing the sensor for detecting optical fiber floor temperature field, a geo-physical model was developed to test the geology and the earth a simulation test was conducted to simulate and analyze the rock stratum temprature field variation induced by inrush on coal face. The result indicates that the distributed fiber testing technology can well identify the temperature variations induced by water temperature change in process of floor water rising and can determine water bursting area quickly and effectively.

        tunnel engineering; geotechnical engineering; floor water bursting; temperature filed; distributed optical fiber sensing; similar simulation

        10.3969/j.issn.1674-0696.2016.05.07

        2015-09-15;

        2015-11-13

        中國煤炭工業(yè)協(xié)會科學技術(shù)研究指導性計劃項目(MTKJ2014-223);安徽省大學生創(chuàng)新創(chuàng)業(yè)訓練計劃項目(201510361126)

        張平松(1971—),男,安徽六安人,教授,博士,主要從事地球物理探查技術(shù)開發(fā)與應用方面的研究。E-mail:pszhang@sohu.com。

        孫斌楊(1992—),男,安徽淮南人,碩士研究生,主要從地球物理勘探方面的研究。E-mail:binyangsun1993@163.com。

        U453.6+1;TD15

        A

        1674-0696(2016)05-028-04

        猜你喜歡
        突水傳感溫度場
        《傳感技術(shù)學報》期刊征訂
        新型無酶便攜式傳感平臺 兩秒內(nèi)測出果蔬農(nóng)藥殘留
        礦井突水水源的判別方法
        礦井滯后突水機理探討
        鋁合金加筋板焊接溫度場和殘余應力數(shù)值模擬
        IPv6與ZigBee無線傳感網(wǎng)互聯(lián)網(wǎng)關(guān)的研究
        電子制作(2018年23期)2018-12-26 01:01:26
        基于紋影法的溫度場分布測量方法
        MJS工法與凍結(jié)法結(jié)合加固區(qū)溫度場研究
        建筑科技(2018年6期)2018-08-30 03:41:08
        李雅莊礦滯后突水主控因素分析
        大相嶺隧道高壓突水機理與預測分析
        亚洲桃色蜜桃av影院| 无码精品人妻一区二区三区人妻斩| 青青草原综合久久大伊人精品| 免费国产自拍在线观看| 国内精品久久久久伊人av| 亚洲日韩精品一区二区三区| 人妻仑乱a级毛片免费看| 中文字幕一区二区三区人妻少妇| 久久婷婷香蕉热狠狠综合| 国产无遮挡无码视频免费软件 | 免费在线亚洲视频观看| 久久综网色亚洲美女亚洲av| 久久精品亚洲精品国产色婷| 亚洲精品乱码久久久久久蜜桃不卡| 免费无码av片在线观看| 中文字幕无码免费久久| 日本亚洲欧美高清专区| 大屁股少妇一区二区无码| 日韩精品一区二区三区视频| 久草福利国产精品资源| 国产欧美精品一区二区三区四区 | 国产一区二区三区在线大屁股| 白浆国产精品一区二区| 曰韩无码无遮挡a级毛片| 久久人人爽人人爽人人av| 欧美伊人网| 大红酸枝极品老料颜色| av男人的天堂亚洲综合网| 人成综合视频在线播放| 国产精品成人3p一区二区三区| 亚洲av无码一区二区二三区| 婷婷亚洲综合五月天小说| 伊人网视频在线观看| 国产精品国产三级在线高清观看| 精品国产午夜久久久久九九| 日本高清一区在线你懂得 | 欧美一级欧美一级在线播放| 加勒比精品一区二区三区| 加勒比一本大道大香蕉| 在线免费观看蜜桃视频| 日韩av无码社区一区二区三区|