999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于滲流場變化的基坑止水帷幕缺陷判別研究與實踐

2016-02-06 09:50:47龐振勇崔王洪
都市快軌交通 2016年6期

龐振勇 崔王洪

(北京城建設計發展集團股份有限公司 北京 100037)

?

基于滲流場變化的基坑止水帷幕缺陷判別研究與實踐

龐振勇 崔王洪

(北京城建設計發展集團股份有限公司 北京 100037)

止水帷幕存在缺陷一直是基坑工程中的重大風險,若能在基坑開挖前判斷出止水帷幕缺陷的位置,便能預先采取加固堵漏措施,可避免基坑開挖時的風險。局部滲漏必然會引起附近滲流場的變化,利用有限元軟件ABAQUS分別對缺陷的長度、深度以及土體滲透性與帷幕外側滲流場分布規律的關系進行研究,得出不同工況下缺陷對滲流場分布的橫向和豎向影響范圍變化規律,證明了根據流速變化判別缺陷存在及位置、規模的可行性。提出采用聲納法可對土體中滲流場分布進行檢測,進而可根據滲流場分布判別出缺陷的規模與位置。通過工程實踐,驗證了研究結論及聲納法檢測技術的準確性,為軌道交通深基坑工程止水帷幕缺陷判別提供指導和借鑒方法。 關鍵詞 城市軌道交通;聲納法;滲漏檢測;ABAQUS軟件;滲流場

地下水的處理一直是透水地層深基坑工程中的重點和難點,對于地質條件復雜、環境保護要求高的區域,止水帷幕若存在缺陷,將導致基坑內外水力聯系無法隔斷,基坑開挖時會引起坑外地面沉降,從而對周邊環境造成較大危害。

1 研究現狀

目前工程中多是采用觀測坑外水位,通過坑內預降水后坑外水位變化來判斷是否存在止水帷幕缺陷,但實際工程中土層條件非常復雜,用此方法并不能準確判別缺陷的位置及規模。為找到一 種 方 法,對 滲 漏進行檢測,并對滲漏點準確定位,從而在工程實施前采取預加固措施,國內外已經進行了大量的研究與實踐。

高密度電阻率法利用水的流動對地層電阻率分布的影響,在同一剖面上測量不同位置和深度的土體電阻率,分析電阻率分布規律的異常,即可確定地層中滲漏的位置與規模[1-2]。瞬變電磁法通過不接地回線向地下發射一次脈沖磁場,地下低阻介質產生感應渦流并在衰減中產生二次磁場傳至地面回線。通過對地面接收的二次磁場空間分布規律的研究,判別出滲漏的位置[3]。探地雷達法利用高頻電磁波以寬頻帶短脈沖的形式,通過天線發射入地下,經過地層差異性反射,再由地面的另一天線接收,對接收到的信號進行處理分析,可知地下介質的空間結構情況并判斷出滲漏的位置[4-5]。同位素示蹤法將利用放射性同位素制作的示蹤液投入待檢測區域水中,若附近存在滲漏點,示蹤液會隨著水體流動,在滲漏點附近集中并被附近土體吸附,通過核探測器檢查附近的核輻射量即可判斷出滲漏點位置及規模[6-8]。溫度示蹤法地層的溫度隨著深度的增加呈規律性變化,但滲漏引起的水流變化會引起溫度場分布異常,通過在一定位置鉆孔并埋設光纖測量地層溫度,即可根據溫度分布異常判別出滲漏點位置[9-12]。以上各種檢測方法,受檢測精度、污染或技術條件的制約,且其經濟性和便利性較差,因此在工程中的應用受到一定的限制。聲納法利用聲波在水體中傳播的方向特性(若與水體流動方向相同其傳播速度加快;若相反,則會減慢),通過采集地層中的聲波信號進行處理即可得到土體中的滲流場分布,并進一步分析出滲漏點的位置及規模[13]。此方法輕便、高效、精確,且經濟性較好,能夠很好地應用于各種規模工程的滲漏檢測。

目前聲納法主要用于水利大壩工程滲漏點檢測,在深基坑特別是地鐵基坑工程中應用較少,理論仍需進一步完善。本文旨在通過有限元模擬,研究地層中三維滲流場分布規律。通過對比不同滲漏參數對地層中滲流場分布規律的影響,根據滲漏點位置及規模判別聲納法檢測技術的可行性,為此提供理論依據,并為其工程應用提供指導。

2 基本理論

由于土體孔隙的形狀和大小很不規則,地下水沿孔隙流動的實際路徑十分復雜,通常研究時假設其符合達西定律。采用有限元法進行計算,需為所取模型設置合適的邊界條件。在滲流計算中,通常采用二類邊界[14]。

2.1 第一類邊界(定水頭邊界)

對模型邊界給定一固定水頭值,并認為其在計算過程中保持不變。

(1)

2.2 第二類邊界(定流量邊界)

認為模型邊界處單位時間內的流量為一定值,不隨時間變化。

(2)

式中,kn為邊界上外法線方向的滲透系數;n為邊界的外法線方向。

Γ1和Γ2構成了三向空間滲流場的全部邊界。

3 計算模型

3.1 基本假定

在數值模擬計算時,為了有針對性地分析所關注的對象,將忽略掉一些次要影響因素,在模擬時采取一些假定,具體如下[15]。

1) 假定土體為均質材料,力學參數各向同性,且在計算過程中土體的密度、彈性模量、滲透系數等參數保持不變。

2) 管井降水之前,土體已在自重作用下完成固結。

3) 設止水帷幕為均質、連續彈性體,假定其滲透系數為0。

4) 為使滲流規律更加明顯,假定止水帷幕為懸掛式止水帷幕。

3.2 模型尺寸

如圖1所示:取基坑大小為100 m×100 m,土體范圍400 m×400 m,止水帷幕按連續墻體考慮,厚度為1 m,長度為35 m。基坑深度為14 m,取含水層厚度為50 m,下部10 m為滲透系數較小的巖層,為相對隔水層,模型豎向共60 m。

圖1 模型尺寸

3.3 邊界條件

通常,第一類邊界條件在滲流剛開始時對滲流場的分布起到支配作用,所以在進行非穩定滲流計算時,必須確定滲流場的初始水頭分布條件。

假定影響范圍之外的土體和流體受基坑降水的影響可忽略不計,故在土體四周邊界約束水平方向位移,流量邊界取定水頭邊界,靜水壓力從地面開始呈線性分布;土體底部約束3個方向的位移??觾仍O置降水井,在降水井底部設置流量邊界,通過調整流量大小控制坑內水位降至坑底以下1 m深度。

3.4 材料參數

以基坑工程中所遇到的較不利土層(粉土、粉砂層)為研究對象,參照JGJ 120—2012《建筑基坑支護技術規程》所提供巖土層滲透系數k的經驗值,取土的滲透系數在0.5~2.5 m/d之間,土體黏聚力c=8 kPa、φ=30°,其他物理參數如表1所示。

表1 材料參數

4 計算結果及分析

通過建立有限元模型,分別研究止水帷幕存在缺陷時,其尺寸、深度及土體的滲透性對滲流場分布規律的影響。因實際工程條件十分復雜,為定性研究不同因素的影響,數值分析時采用單因素變化的方法對比研究。

考慮到實際施工的便利性,流速檢測孔宜設置在帷幕外側一定距離且對應樁(墻)接縫處,在帷幕外側1 m處取平行于帷幕方向的2-2剖面(見圖1),如圖2流速分布云所示,在距缺陷中心線(缺陷位置1-1與2-2剖面的交線)不同距離處分別提取沿土層豎向不同深度的流速數據,并對比分析,研究缺陷對周邊土體中三維流速分布的影響規律。

圖2 2-2剖面上數據提取位置

取標準工況下止水帷幕深度35 m,降水井深度25 m,土體滲透性按各向同性處理,水平、豎向滲透系數均為0.5 m/d。各單因素分析時,其他參數均與標準工況相同。

4.1 止水帷幕缺陷長度的影響

研究缺陷尺寸變化時,帷幕外側流速分布變化規律。考慮缺陷未完全貫通,缺陷體處的滲透系數小于周邊土體,取值0.25 m/d。考慮到實際工程中若止水帷幕出現缺陷,一般在墻(樁)接縫處,缺陷形狀多為長條形,取缺陷寬度為0.1 m(以下分析中,缺陷寬度均按帷幕接縫止水局部失效后因垂直度偏差產生縫寬的平均值0.1 m計取,厚度同止水帷幕),缺陷中心埋深為9 m(以下缺陷深度即指缺陷中心深度),豎向長度分別取0.2、0.5、1.0、2.0、3.0 m。

由圖3可以看出,止水帷幕缺陷體處流速值遠大于2-2剖面上缺陷深度處流速值且分布規律相反。在缺陷寬度相同的條件下,缺陷體處流速值變化與長度呈反比關系。當缺陷長度大于1 m后,最大流速值隨長度增大而減小的速率降低。原因是缺陷尺寸變小造成過流斷面也減小,引起局部流速急劇增大,但隨著缺陷長度的增加,其影響逐漸減小。帷幕外側1 m處與缺陷體處最大流速值分布規律相反的原因是:缺陷尺寸較小時,雖然缺陷體處因過流斷面小而流速增大,但其對帷幕外側土體中流速分布的影響范圍隨缺陷尺寸減小而減小。

圖3 缺陷體與2-2剖面缺陷深度處最大流速值變化

由圖4可以看出,存在缺陷時,缺陷深度附近流速值明顯大于帷幕完整時的數值。當缺陷長度較小時(小于1 m),隨著長度的增加,缺陷深度附近的流速值明顯增大,但當缺陷長度較大時(大于1 m),其影響范圍遠大于1 m,缺陷長度增加對在帷幕外側1 m剖面中心線上缺陷深度處的流速峰值影響較小,故缺陷深度附近的流速值變化不大,只是流速異常變化的深度范圍隨缺陷長度的增加而增大。因此,憑借流速值大小無法判定缺陷規模,但可根據流速異常變化的豎向范圍確定缺陷的長度。

圖4 缺陷中心線流速值分布(長度影響)

結合圖5可以看出,缺陷長度增加,其橫向影響范圍隨之增大。缺陷長度為0.2 m時,距離缺陷中心線6 m處流速值變化較小,因此為保證檢測的精確度,帷幕外側的測點需保證一定的密度,測點距離缺陷隱患處不宜大于6 m。

圖5 距缺陷中心線不同距離處流速值分布(長度影響)

4.2 止水帷幕缺陷深度的影響

取帷幕缺陷長度均為2 m,缺陷中心埋深(缺陷深度)分別為9、13、17、21、25 m,考慮缺陷未完全貫通,缺陷體處的滲透系數按周邊土體的一半取值為0.25 m/d,對比分析缺陷深度對滲流場分布的影響。

由圖6可以看出,隨著缺陷深度的增加,流速值的變化主要出現在缺陷所處深度附近,其他部位基本不受影響,當缺陷長度相同的條件下,缺陷對滲流場的豎向影響范圍隨缺陷所處深度增加整體向下移動,影響范圍大小并沒有明顯變化。

圖6 缺陷中心線流速值分布(深度影響)

圖7 距缺陷中心線不同距離處流速值分布(深度影響)

結合圖7可以發現,與帷幕完整時對比,當缺陷所處深度為9 m時,距離缺陷中心8.5 m的缺陷所處深度的流速變化異常仍非常明顯,當缺陷所處深度為25 m時,距離缺陷中心6.0 m的流速變化已趨于平緩。這說明,隨著缺陷所處深度的增加,其橫向影響范圍逐漸減小,原因在于隨著深度的增加,流速值越來越大,缺陷引起的流速值的增加與相應位置正常情況下的流速相比越來越不明顯。因此,應根據實際工程中止水帷幕深度的增加,適當加密橫向數據采集點。

4.3 土體滲透性的影響

取缺陷中心埋深為9 m,缺陷長度為2 m,土體的滲透系數取0.5、1.0、1.5、2.0、2.5 m/d 5種工況,考慮缺陷未完全貫通,缺陷體處的滲透系數取值為周邊土體的一半,對比分析土層滲透性對滲流場分布的影響。

由圖8可以看出,流速值隨著滲透系數的增大而明顯增加,但增加幅度逐漸降低,說明隨著滲透系數的增大,其對流速值的影響逐漸降低。結合圖9可以看出,隨著土體滲透系數的增加,缺陷的橫向和豎向影響范圍并沒有明顯變化,但滲透系數較小時,流速值整體較小,對流速測量的精度要求更高,故測點分布密度應結合實際地層參數確定。

圖8 缺陷中心線流速值分布(滲透性影響)

圖9 距缺陷中心線不同距離處流速值分布(滲透性影響)

5 實踐應用

5.1 測量原理

若被測水體存在滲流,則必然在測點附近產生滲流場,而水體的流動必然會產生聲波,聲波在水體中傳播的方向若與水體流動方向相同,其傳播速度加快;若相反,則會減慢。聲納法滲流檢測技術,正是利用聲波在水中的優異傳播特性,實現對水體流速場的測量,通過在檢測孔不同深度處采集滲流所產生的聲波信號,并對接收到的信號進行處理,可判斷滲流場分布是否存在異常,進而判斷出是否存在缺陷及其位置。根據《水工混凝土缺陷檢測技術規程》(SL713-2015)中的要求,聲納流速矢量測量儀測量流速精度應達到8.64×10-3m/d(1×10-5cm/s),可精確檢測到透水性高于弱透水地層的流速變化,目前市場上工程應用的聲納矢量測量儀均能滿足規范要求。

5.2 工程實踐

基于上述滲流場理論研究成果并結合聲納法檢測原理,在南京河西新鴻基CBD項目的基坑工程中進行了抽水試驗,并利用聲納法滲流檢測技術對地下連續墻的質量進行檢測,采用的DB-5型聲納測量儀精度為8.64×10-3m/d,判別出了地下連續墻結構存在缺陷的位置并進行修補,有效保證了基坑開挖的順利進行。

5.2.1 工程概況

本工程位于南京地鐵1、2號線元通站東南側,場地為長江漫灘地貌單元,基坑長250 m,寬105 m,開挖深度為20.3~22.8 m。東北側和西側平行基坑方向有運行中的地鐵1、2號線隧道,基坑北側距1號線隧道僅13.6 m,距車站風亭僅5.8 m,西側距離2號線元通站最近出入口約為11.5 m。

基坑支護結構采用1 m厚地下連續墻加4道鋼筋砼內支撐體系。為確?;娱_挖承壓水穩定并減小坑外水位變化對既有地鐵結構的影響,地下連續墻插入不透水風化巖3.5 m,隔斷承壓水層,形成封閉隔水體系。

根據地勘資料,基坑開挖范圍及坑底到基巖面之間存在較厚的粉細砂層,因此工程對止水帷幕的隔水效果要求非常高。為保證帷幕(地下連續墻)的質量,避免基坑開挖對周邊環境的影響,采用聲納法滲流檢測技術對工程止水帷幕質量進行檢測。

圖10 滲流檢測孔布置

考慮到地連墻帷幕缺陷多是出現在墻幅接縫處,故觀測孔結合地連墻分幅,在臨近地鐵側沿基坑邊對應地連墻接縫處設置,并根據施工日志的記錄選擇存在質量隱患的位置重點檢測。觀測孔布置如圖10所示。檢測孔統一布置在帷幕外側1.2 m的位置。

5.2.2 數據分析

根據各測孔所得數據,分析得出不同測孔的豎向流速分布曲線,發現有個別測孔流速分布有異常,如圖11所示黃色塊部位。

圖11 止水帷幕滲漏檢測數據三維分布及剖面

第4節研究結果表明,當止水帷幕存在缺陷時,其所在深度附近一定范圍內的流速值異常變化都會比較明顯,因此試驗數據中流速異常變化范圍較小時可以忽略。圖11中數據表明,有6個測孔在15~25 m深度范圍出現了較大的滲漏流速值,分別為S2、S4、S6、S12、S21、S24,尤其是S12孔和S24孔兩個測孔的流速值最大,它們的滲透流速最大值分別達到了15 m/d和4.6 m/d,推測是由于此處缺陷規模較大且缺陷所處地層滲透性較強。在35~55 m深度范圍也有4個測孔的數據明顯較大,分別為S2、S25、S26、S27,尤其是S25測孔的流速值最大,滲流流速最大值均超過1 m/d。根據對比可以初步判定這些測孔附近的止水帷幕在相應深度范圍內存在滲漏(見圖12)。

根據4.1節關于缺陷尺寸對滲流場的影響規律的討論可以判斷,測孔S12和S24附近的止水帷幕缺陷尺寸最大。對S24孔數據進行處理得到三維成像(見圖13),發現該測孔在深度17.25 m附近流速明顯大于其他區域,可判斷此測點附近止水帷幕存在滲漏,根據流向數據可判斷滲漏點與測孔在平面上的方向關系,進而可綜合判定帷幕缺陷的具體位置。根據此方法對其他測孔數據進行處理,可綜合判斷出不同區域滲漏點的具體位置。

本次檢測之后,根據檢測結果,對滲漏量較大的S2、S12、S21、S24、S25和S26附近的地下連續墻接縫處采用了高壓旋噴樁預加固,加固深度范圍為5.0~60.9 m。在基坑開挖過程中對存在滲漏的位置進行驗證,如圖14所示,S24孔附近地下連續墻在地下17.25 m附近接縫處確實存在施工缺陷,因在開挖之前采取了預加固措施,故基坑開挖期間并未出現明顯滲漏,開挖得以順利進行,驗證了上述理論分析及聲納法檢測技術的可行性與準確性。目前該工程已順利通過竣工驗收,南京河西長江漫灘后續的幾個深基坑工程也采用了此技術進行基坑止水帷幕滲流檢測,均取得良好效果。

6 結論

1) 缺陷的存在導致附近滲流場出現明顯變化,進一步發展會引起周邊土體的破壞。缺陷長度較大時,缺陷長度的增加對帷幕外側一定距離處缺陷深度附近的流速值影響減小,根據檢測流速值無法判定缺陷的規模,可根據流速異常變化的范圍確定缺陷長度。為保證檢測精度,測點距離缺陷隱患位置距離不宜大于6 m。

圖12 各測孔流速沿深度分布曲線

圖13 S24檢測數據三維成像圖

圖14 滲漏缺陷處開挖驗證

2) 隨著缺陷所處深度增加,缺陷對滲流場的豎向影響范圍沒有明顯變化,但橫向影響范圍有所減小,因此缺陷所處深度不同時仍可根據其豎向影響范圍確定缺陷尺寸,但應根據所要檢測的帷幕深度選擇適當的測孔間距。

3) 土體滲透系數增大時,滲流場的橫向和豎向影響范圍沒有明顯變化,但土體滲透系數較小時,因流速值較小,所以對流速測量的精度要求較高,故測點分布密度應根據實際地層參數確定。

4) 根據現場檢測數據,結合數值分析得出各種因素對滲流場分布的影響規律,推測出工程滲漏風險點并預先采取加固措施,保證了工程順利完工,驗證了研究結論的正確性,可為類似軌道交通深基坑工程利用聲納法滲流檢測技術判別帷幕缺陷位置及規模,提供指導和借鑒意義。

[1] 趙志偉,徐旺敏,傅瓊華.高密度電阻率法在土壩滲漏檢測中的應用[J].江西水利科技,2011, 37(4):266-268.

[2] 李富,劉樹才,曹軍,等.高密度電阻率法在工程勘察中的應用[J].工程地球物理學報,2006, 3(2):119-123.

[3] 肖衡林,周錦華.滲漏監測技術研究進展[J].中國水運,2007, 7(2):87-91.

[4] 陳義群,肖柏勛.論探地雷達現狀與發展[J].工程地球物理學報,2005, 2(2):149-155.

[5] 葛雙成,江影,顏學軍.綜合物探技術在堤壩隱患探測中的應用[J].地球物理學進展,2006, 21(1):263-272.

[6] YAJUN L.Application of chemical tracing experiment technique in leakage detection of hydraulic engineering[J].agricultural engineering and agricultural machinery,2011, 12(9):1385-1387.

[7] 于瑞蓮,胡恭任,袁星,等.同位素示蹤技術在沉積物重金屬污染溯源中的應用[J].地球與環境,2008, 36(3):245-250.

[8] BLAIS J M.Using isotopic tracers in lake sediments to assess

atmospheric transport of lead in Eastern Canada[J].Water, air and soil pollution,1996, 92(3-4): 329-342.

[9] 吳志偉,宋漢周.地下水溫度示蹤理論與方法研究進展[J].水科學進展,2011, 22(5):733-740.

[10] 董海洲,陳建生.利用溫度示蹤方法探測基坑滲漏[J].巖土力學與工程學報,2004, 23(12):2085-2090.

[11] 肖衡林,張晉鋒,何俊.基于分布式光纖傳感技術的流速測量方法研究[J].巖土工程,2009, 30(11):3543-3547.[12] BECKER M W, Georgian T, Ambrose H, et al.Estimating flow and flux of ground water discharge using water temperature and velocity[J].Journal of hydrology,2004(296): 221-233.[13] 劉迪,李雪嬌,于艷秋.聲納滲流檢測于橋水庫大壩滲漏點的應用研究[J].河海水利,2013(3):46-47.

[14] 張學飛.基坑工程非穩定滲流場的有限元分析[D].太原: 太原理工大學, 2010.

[15] 金小榮.基坑降水對周圍環境影響的數值模擬分析[D].杭州: 浙江大學, 2004.

(編輯:郝京紅)

Research and Practice on Identifying Defects of Waterproof Curtain for Excavation Engineering Based on Changes of Seepage Field

Pang Zhenyong Cui Wanghong

(Beijing Urban Construction Design & Development Group Co., Ltd., Beijing, 100037)

The existence of defects in the sealing curtain has been a major risk in a pit project. If the location of defects of a water curtain can be identified before a pit is excavated, measures of plugging and reinforcement can be taken in advance to avoid the risk of pit excavation. The local seepage will inevitably cause the change of the seepage field. The relationship between the length and depth of the defect and the distribution of the seepage field of the curtain is studied by the finite element software ABAQUS. The variation of horizontal and vertical influence scope of defects on the seepage field under different conditions are explored. It is proved that it is feasible to determine whether the defects exists, where they are located and how large they are. It is proposed that the distribution of seepage field in soil can be detected by sonar method, and the size and position of defects can be detected according to the distribution of the seepage field. Through the engineering practice, it verifies the conclusion of the research and the accuracy of the sonar method, and provides guidance and references to detecting defects of waterproof curtain for deep excavation engineering in urban rail transit projects.Key words: urban rail transit; sonar technology; leakage detection; ABAQUS; seepage field

10.3969/j.issn.1672-6073.2016.06.020

2016-08-04

2016-10-23

龐振勇,男,碩士,高級工程師,從事軌道交通結構設計研究工作,41676902@qq.com

U231

A

1672-6073(2016)06-0099-07

主站蜘蛛池模板: 六月婷婷综合| 婷婷色婷婷| 成人精品免费视频| 四虎永久免费网站| 中国一级特黄视频| 亚洲最大福利网站| 69视频国产| 91精品在线视频观看| 制服丝袜亚洲| 欧美在线国产| 亚洲精品国产精品乱码不卞 | 国产精品男人的天堂| 国产尤物在线播放| 欧美黄色a| 九九热精品视频在线| 久久人人97超碰人人澡爱香蕉 | 国产精品 欧美激情 在线播放| 高清无码手机在线观看| 国产成人精品男人的天堂下载| YW尤物AV无码国产在线观看| 国产欧美精品一区aⅴ影院| 一本无码在线观看| 久久综合干| 亚洲欧美日韩高清综合678| 国产99热| 欧美日韩福利| 亚洲欧美综合精品久久成人网| 亚洲欧美成aⅴ人在线观看| 国产免费精彩视频| 任我操在线视频| 熟妇丰满人妻| 91成人在线免费视频| 欧美成人看片一区二区三区| 精品视频福利| 日本精品影院| 亚洲三级影院| 日韩在线视频网站| 久久永久视频| 波多野结衣一级毛片| 另类欧美日韩| 欧美区一区| 国产一区免费在线观看| 国产精品成人啪精品视频| 久久综合激情网| 青青草欧美| 色综合久久88色综合天天提莫| 国产第一页免费浮力影院| 国产免费久久精品99re丫丫一 | 国产免费高清无需播放器| 狼友视频一区二区三区| 亚洲天堂精品视频| 欧美成人在线免费| 国产成人免费| 国产全黄a一级毛片| 精品国产成人av免费| 日韩一区二区三免费高清| 色国产视频| 成人在线综合| 国产中文一区二区苍井空| 国产午夜不卡| 国产无遮挡猛进猛出免费软件| 无遮挡国产高潮视频免费观看 | 国产日韩欧美在线播放| 久久影院一区二区h| 亚洲人成网址| 国产一区亚洲一区| 国产久草视频| 亚洲六月丁香六月婷婷蜜芽| 亚洲无线国产观看| 操操操综合网| 国产人人射| 国产精品制服| 热这里只有精品国产热门精品| 波多野结衣无码视频在线观看| 久久99这里精品8国产| 国产成人免费高清AⅤ| 无码内射中文字幕岛国片| 亚洲一区国色天香| 国产午夜福利在线小视频| 毛片免费高清免费| 久久一日本道色综合久久| 日本a级免费|