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

    基于SBAS技術(shù)的金屬礦山沉陷規(guī)律研究

    2015-03-28 06:11:22雷廣淵
    測繪工程 2015年3期
    關(guān)鍵詞:中段采空區(qū)礦床

    雷廣淵,周 輝

    (1.廣州市城市規(guī)劃勘測設(shè)計研究院,廣東 廣州510060;2.廣東交通職業(yè)技術(shù)學(xué)院 公路學(xué)院,廣東 廣州510650)

    礦產(chǎn)資源是國民經(jīng)濟(jì)中重要組成部分,也是工業(yè)發(fā)展的重要物質(zhì)保障。礦產(chǎn)資源的開采造成了大量地表塌陷、建筑物變形、道路破壞等,嚴(yán)重危害人類生命財產(chǎn)安全。目前煤炭開采的沉陷規(guī)律的研究已有顯著成果,但是金屬礦山開采沉陷研究國內(nèi)外仍然處于探索階段,尚且沒有形成成熟的理論,主要原因在于金屬礦山地質(zhì)采礦條件與煤炭相比要復(fù)雜很多。研究金屬礦山開采沉陷規(guī)律可以為指導(dǎo)和管理礦山開采提供依據(jù)。

    地下有用礦物被開采之后,原有的應(yīng)力平衡被破壞,巖層移動波及到地表,損害礦區(qū)地表的建筑物、鐵路等。文獻(xiàn)[1]對構(gòu)造應(yīng)力下的地表移動規(guī)律做了研究,總結(jié)出了急傾斜崩落法的地表移動規(guī)律。文獻(xiàn)[2-3]對鐵礦的地表移動規(guī)律做了探索性研究,總結(jié)出崩落法開采鐵礦石可采用概率積分法實(shí)現(xiàn),并將其應(yīng)用在程潮、北洺河鐵礦中。但是目前對礦區(qū)開采沉陷的研究都是基于傳統(tǒng)測量方式為基礎(chǔ)進(jìn)行分析研究,本文將采用小基線集(Smallbaseline Subset,SBAS)技術(shù)對某金屬礦山的沉陷規(guī)律進(jìn)行探索研究。針對復(fù)雜的時空沉陷過程,研究其沉陷規(guī)律和下沉機(jī)理,SBAS技術(shù)避免了基線的限制,可以生成時序的地表形變。本文針對湘中地區(qū)某大型的銻金屬礦山,利用SBAS技術(shù)綜合分析該金屬礦山的巖層與地表移動規(guī)律,指導(dǎo)礦山開采并作為類似金屬礦山的工程類比依據(jù)。

    1 SBAS技術(shù)原理

    SBAS技術(shù)將影像通過基線的長短分成多個小集合,對于每個基線差異較小的集合,利用最小二乘方法計算,得到地表形變時間序列。在解算的過程中,由于影像之間存在相關(guān)性以及小基線集的不同組合會使方程產(chǎn)生秩虧,對于秩虧問題,一般采用奇異值分解附加最小范數(shù)的方法解算。

    假設(shè)針對N+1幅影像,時間向量t=[t0,…,tN]T,選取合適的基線進(jìn)行干涉,生成G個差分干涉圖,其中

    忽略大氣誤差、去相關(guān)以及高程誤差的影響,對于時間tA和tB模型簡化為

    其中:λ為雷達(dá)波長,d(t,i,j)為t時刻雷達(dá)視線方向相對參考時刻的形變。

    SBAS是按像元逐個計算時間序列的形變,所以不失一般性,以某一像元建立關(guān)系。設(shè)SAR影像中該 像 元 的 相 位 表 示 為 φT= [φ(t0),φ(t1),…,φ(tN)],為待求量,相干圖中的觀測相位表示為δφT=[δφ1,δφ2,…,δφM],為觀測量。其中φT表示像元在不同時刻的相位,δφT表示干涉圖中相對于主參考點(diǎn)解纏后的相位。若對主輔影像按時間進(jìn)行排序,則主影像的時間大于輔影像的時間,對于干涉圖中的相位可表示為

    δφi=φ(tMaster(i))-φ(tSlave(i)),i=1,2,…,M.(3)用矩陣表示為δφ=Aφ

    A矩陣每一行有兩個非零元素(1和-1),代表一個干涉圖,矩陣的列代表對應(yīng)時間的一景SAR影像。若M>N,可以利用最小二乘方法獲得唯一解,若M與N接近或者小于N,需采用奇異值分解以及最小范數(shù)求解。公式可表示為

    帶入公式δφ=Aφ,可得

    表示第j幅干涉圖的相位值可以用兩幅主輔影像的獲取時間區(qū)間內(nèi)各時間段的沉降速率的積分求得,寫成矩陣表示為

    其中B為M×N的矩陣,在主輔影像獲取時間之間,B的系數(shù)不全為0,否則B=0。與式(5)同樣需要奇異值分解,得到速度的最小范數(shù)解。時間段上的形變可以利用速度在時間內(nèi)的積分計算。

    2 工程概況

    研究礦山位于湘中地區(qū),該礦山屬于大型銻礦,銻儲量非常豐富。礦山分南北兩個礦區(qū),本文針對北礦區(qū)的某礦床進(jìn)行分析。該礦區(qū)地質(zhì)構(gòu)造復(fù)雜,處于多種構(gòu)造體系的交界部位,而且受到不同性質(zhì)的構(gòu)造應(yīng)力作用。礦體主要受背斜和F3斷層控制。礦田賦存Ⅰ、Ⅱ、Ⅲ號3個礦體,Ⅰ、Ⅱ號為主要開采礦體。礦體呈不規(guī)則的似層狀和透鏡狀,傾角15°左右,礦體厚度為0.56~8.81 m,平均2.43 m。頂板巖層為頁巖、硅化灰?guī)r、破碎的灰、頁巖。主要采礦方法有普通房柱法和膠結(jié)充填法。目前礦床共有8個中段,每個中段間隔距離30 m,一中段標(biāo)高為+567 m,八中段標(biāo)高為+357 m。四中段以上礦體已開采完畢,但是對于部分礦體殘留,仍進(jìn)行殘采。五、六中段為主要開采中段開采礦體為Ⅰ、Ⅱ號礦體,八中段目前仍在開拓階段。

    3 實(shí)驗(yàn)分析

    3.1 實(shí)驗(yàn)數(shù)據(jù)

    實(shí)驗(yàn)采用的數(shù)據(jù)包括15景日本空間局ALOS衛(wèi)星的PALSAR數(shù)據(jù)以及NASA的SRT M(Shuttle Radar Topography Mission,航天飛機(jī)雷達(dá)地形測繪任務(wù))數(shù)據(jù)。ALOS衛(wèi)星獲取SAR影像的時間分別 為:2007-01-25、2007-06-12、2007-07-28、2007-12-13、2008-08-14、2008-12-15、2009-01-30、2009-08-02、2009-09-17、2009-12-18、2010-02-02、2010-03-20、2010-06-20、2010-09-20、2011-02-05。SRT M數(shù)據(jù)在網(wǎng)上公開發(fā)布,其DEM水平分辨率為90 m,數(shù)據(jù)可以從JPL提供的地址:htt p://dds.cr.usgs.gov/srt m/下載。覆蓋實(shí)驗(yàn)區(qū)域以及周邊區(qū)域的DEM 有 9 塊,分 別 為 N26E110、N26E111、N26E112 N27E110 N27E111 N27E112 N28E110、N28E111、N28E112。

    3.2 SBAS監(jiān)測金屬礦山

    采用瑞士Gamma軟件作為實(shí)驗(yàn)平臺,對15景PALSAR數(shù)據(jù)使用SBAS技術(shù)進(jìn)行處理。假設(shè)2007-01-25作為參考影像,即形變?yōu)?,那么可以得到14個時間段的形變序列,地表形變序列如圖1所示。

    圖1 地理編碼后形變序列圖

    由圖1可知,影像的中間區(qū)域在監(jiān)測時間段內(nèi)發(fā)生較大的形變,根據(jù)資料分析得知,該區(qū)域有大量煤礦,煤炭的大量開采導(dǎo)致沉降范圍較大,形變也較為顯著,在監(jiān)測周期內(nèi)發(fā)生最大變形約為60 c m。研究區(qū)域位于該影像的左下角,該區(qū)域在監(jiān)測周期內(nèi)發(fā)生地表形變達(dá)16 c m。該區(qū)域?yàn)槟炒笮徒饘俚V山,本文主要針對這一局域進(jìn)行分析,從圖像中提取該礦地表沉降序列,見圖2。

    圖2 地理編碼后地表形變序列

    選取2011-02-05的形變進(jìn)行分析,如圖3所示,該區(qū)域?yàn)榈V床地表的沉降情況,該區(qū)域在2007-01-25至 2011-02-05 監(jiān) 測 期 間 地 表 最 大 沉 降 達(dá)16 c m,沉降走向?yàn)楸逼珫|50°左右。

    圖3 2011-02-05累積沉降等值線圖(單位:m)

    3.3 礦山地表移動規(guī)律探討

    為了深入了解該礦床地表的沉陷規(guī)律,本文對沉降序列進(jìn)行分析。該礦床的礦體走向?yàn)楸睎|東,與地表沉降的發(fā)展方向一致。分別選取礦體走向方向和傾向方向的兩個斷面進(jìn)行分析,圖4、圖5分別為礦體走向方向上和傾向方向上的斷面沉降序列。

    圖4 礦體走向方向上的地表下沉曲線

    圖5 礦體傾向方向上的地表下沉曲線

    由圖4、圖5所示,該礦床的地表塌陷在空間上與高斯曲線相似,這說明在一定范圍內(nèi),該金屬礦床的地表移動規(guī)律宏觀上符合概率積分法的原理,經(jīng)實(shí)地調(diào)查發(fā)現(xiàn),沉降較大的區(qū)域主要是受該礦床第三中段開采的影響,四中段和五中段在不同程度上也影響上覆地表。由于無法獲得地下完整的采空區(qū)數(shù)據(jù),本文采用高斯曲線擬合沉降值,從曲線的形態(tài)上分析下沉盆地的規(guī)律。擬合的結(jié)果見圖6、圖7,使用的高斯擬合函數(shù)為

    通過曲線擬合求得的傾向方向和走向方向上的高斯曲線結(jié)果可知,曲線在宏觀上與實(shí)際數(shù)據(jù)擬合效果良好,但在局部區(qū)域擬合結(jié)果較差,可能由于SBAS觀測精度的問題或者地表變化不連續(xù)。在曲線的邊緣擬合效果較差,這一問題可能的主要原因有:①概率積分法模型誤差,較多工程實(shí)例表明,概率積分法在邊緣擬合效果較差。②實(shí)驗(yàn)區(qū)域?yàn)榻饘俚V山,采礦方法也是采用房柱法等,上覆巖層并不是理想的隨機(jī)游動模型,部分區(qū)域可能存在斷層、構(gòu)造應(yīng)力等的多重影響。

    圖6 傾向方向的高斯曲線擬合以及殘差分析

    走向方向和傾向方向求得的高斯曲線的參數(shù)和殘差中誤差見表1。

    圖7 走向方向的高斯曲線擬合以及殘差分析

    表1 高斯擬合函數(shù)相關(guān)指標(biāo)

    表1中均方根誤差

    3.4 地質(zhì)采礦因素與地表下沉聯(lián)系

    礦山開采導(dǎo)致了地表的下沉,大范圍的地表下沉不僅對地表建筑物和道路造成損害,而且使得在開采過程中,由于上覆巖層的破裂,坍塌等,造成采場、運(yùn)輸巷道、通風(fēng)系統(tǒng)等受到影響,無法進(jìn)行安全生產(chǎn),從而危害人類的生命財產(chǎn)安全。為了深入了解礦床的地表移動規(guī)律,進(jìn)一步分析礦山開采與地表下沉之間的聯(lián)系。圖8反映了該礦床35號勘探線3個礦體空間分布特征以及相應(yīng)的地質(zhì)條件。由圖中可知,該礦床主要受背斜以及F3(右第3條斷層)斷層控制,即控礦因素為“背斜加一刀”,I和II號礦體為主要開采礦體。

    為了分析礦體開采與地表沉降的因果關(guān)系,將InSAR獲取的沉降場與地形圖配準(zhǔn),配準(zhǔn)結(jié)果見圖9。

    圖8 35勘探線剖面圖

    圖9 InSAR沉降值與礦區(qū)開采圖配準(zhǔn)結(jié)果

    分析圖中A和B兩個區(qū)域的變形特點(diǎn),圖中所示的時間是2011年2月(最后一次觀測)的累積形變,從圖中可知,在14次的監(jiān)測中,A區(qū)域沒有明顯的沉降,B區(qū)域有較大的變形。為分析A、B區(qū)域的形變原因,對以上區(qū)域進(jìn)行的采礦活動進(jìn)行統(tǒng)計。通過2012采選廠儲量計算圖、各中段采礦驗(yàn)收報告中得知:

    A區(qū)域是第五中段,采場編號為5100~5127以及采場5211~5217等,采深大約在150 m左右,在2008年至2009年期間,該區(qū)域進(jìn)行了采礦活動,主要對I號薄礦體和II號厚礦體進(jìn)行了開采。開采的方法為充填法開采,由于開采離地表較遠(yuǎn),而且對大部分采空區(qū)進(jìn)行了管理,地表在此期間沒有發(fā)生大范圍、劇烈的形變。

    B區(qū)域?yàn)?4線到80線之間的區(qū)域,該區(qū)域礦體深度主要分布在90~150 m,即礦體主要分布在三、四、五中段。由于礦區(qū)開采的年限較長,三中段和四中段采空區(qū)未進(jìn)行處理,留有較大的采空區(qū)。采空區(qū)由于仍有礦體遺留(礦柱、邊角等礦體),在2007年至2011年期間,該區(qū)域三中段和四中段一直存在殘采,礦柱和邊角礦體被開采之后,對采空區(qū)未進(jìn)行處理,采空區(qū)離地表又較為接近,所以直接導(dǎo)致了地表產(chǎn)生了下沉。五中段開采主要由于離地表較遠(yuǎn),部分采場也進(jìn)行了充填,對上覆巖體進(jìn)行了支撐,故可認(rèn)為對地表影響較小。所以該區(qū)域主要的地表下沉一方面是由于存在大量空區(qū)所致,另一方面,礦柱和邊角礦體的回收,導(dǎo)致上覆巖體失去支撐,重復(fù)采動也對地表影響加劇,最終導(dǎo)致了該區(qū)域地表下沉較大。表2為對該礦床空區(qū)的統(tǒng)計,從表中可知,三中段采空區(qū)體積最大,四中段次之,從圖9中可知,下沉較大的區(qū)域主要范圍是三中段的巷道之上,這與地表下沉較大的區(qū)域較為符合。

    表2 礦區(qū)空區(qū)統(tǒng)計

    為更加清楚地認(rèn)識礦山開采與地表下沉的聯(lián)系,本文對66號勘探線上的沉降值進(jìn)行分析。66號勘探線位于礦床的中上部,在勘探線上每10 m取一個下沉點(diǎn),下沉序列見圖10。64~70勘探線在2008年之后對礦床進(jìn)行殘采工作,15號點(diǎn)附近為下沉較為劇烈的區(qū)域,該區(qū)域?yàn)槿卸尾蓤?49~259殘采范圍。隨著點(diǎn)號的增加,地表下沉變緩,點(diǎn)號23附近,在2008-2010年期間五中段5217~5243采場進(jìn)行了采礦活動,由于離地表相對較遠(yuǎn),地表變化變緩。

    圖10 66號勘探線地表下沉序列

    傾斜曲線表示地表的傾斜以及地表的下沉速度大小,在獲得了66號勘探線上沉降曲線之后,利用下沉值計算傾斜曲線,對最后一期數(shù)據(jù)進(jìn)行計算獲得的傾斜曲線如圖11所示。由于影像受到大氣、失相關(guān)、影像分辨率等因素影響,雖然下沉曲線在趨勢上呈U型,但是在局部區(qū)域,相鄰地表沉降值出現(xiàn)相等或者相近,所以在傾斜計算中,使得曲線出現(xiàn)跳動,與實(shí)際的物理意義不符合。為了削弱這種由于上述因素影響下的相鄰沉降值出現(xiàn)相等或者相近而引起的傾斜曲線跳動,本文使用一種濾波去噪算法,既要去掉大部分干擾噪聲,又不失原始下沉曲線的沉降特點(diǎn),并且保證修正后的數(shù)據(jù)與原始數(shù)據(jù)差距在一定范圍內(nèi),擬合的數(shù)據(jù)能夠反應(yīng)地表下沉的主要特征。本文采用滑動平均對原始數(shù)據(jù)進(jìn)行平滑,選用的窗口為5,目的是消除相鄰兩個點(diǎn)沉降值相等或者相近以及噪聲等因素的影響,實(shí)驗(yàn)對原始數(shù)據(jù)分別進(jìn)行了一次、二次、三次平滑,結(jié)果見圖11、圖12。

    圖11 原始沉降進(jìn)行修正

    圖12 修正后的傾斜曲線

    平滑3次后,結(jié)果見圖13、圖14。

    圖13 利用滑動平均修正沉降數(shù)據(jù)圖

    圖14 利用滑動平均修正沉降數(shù)據(jù)后傾斜曲線

    由圖13、圖14知:地表最大下沉值在點(diǎn)號15附近,此處計算的傾斜接近0,傾斜在15號點(diǎn)之前,有一個增大和下降的過程,與開采沉陷的理論中傾斜曲線和下沉曲線的關(guān)系符合。點(diǎn)號15之后傾斜仍然在下降,下降達(dá)到最小(負(fù)極值)后開始遞增,這里傾斜并沒有繼續(xù)增大,形成一個對稱的曲線,原因是由于20點(diǎn)之后,礦山開采的深度和采礦方法以及頂板管理方法等因素的制約,導(dǎo)致曲線變緩。

    4 結(jié) 論

    1)該金屬礦山淺層(三、四中段)有較大采空區(qū),地表塌陷規(guī)律在宏觀上服從高斯曲線,即地表受重力影響較大,沖擊力破壞較少。

    2)由觀測數(shù)據(jù)可知,在部分區(qū)域出現(xiàn)跳動,可能是由于構(gòu)造應(yīng)力存在或者是由于小斷層、背斜構(gòu)造綜合影響。

    3)礦山地表下沉主要因素是由于三中段和四中段大量的空區(qū)以及礦柱等回采引起,五中段以上對采空區(qū)進(jìn)行了不同程度的支護(hù)、填充。礦山開采對地表沉降影響較小。

    總之,小基線應(yīng)用在金屬礦山開采規(guī)律分析還處于探索階段,從實(shí)驗(yàn)中得出該礦床地表塌陷規(guī)律在淺層由于空區(qū)范圍大,上覆巖層為頁巖等,受重力影響顯著,在宏觀上可以利用概率積分法進(jìn)行分析。局部區(qū)域由于受到復(fù)雜的地質(zhì)構(gòu)造影響,對沉陷規(guī)律分析需要綜合構(gòu)造應(yīng)力、斷層以及背斜等因素綜合分析。

    [1] 賀躍光,顏榮貴,曾卓喬.構(gòu)造應(yīng)力作用下的地表移動規(guī)律研究[J].礦冶工程,2000,20(3):12-14.

    [2] 白義如,谷志孟,白世偉.程潮鐵礦東區(qū)地下采礦引起地表沉降和巖層移動初探[J].巖石力學(xué)與工程學(xué)報,2002,21(3):340-342.

    [3] 李春雷,謝謨文,李曉璐.基于GIS和概率積分法的北洺河鐵礦開采沉陷預(yù)測及應(yīng)用[J].巖石力學(xué)與工程學(xué)報,2007,26(6):1243-1250.

    [4] 葉粵文.金屬礦山地壓災(zāi)害發(fā)生機(jī)理研究[J].有色金屬(礦山部分),2006,58(1):30-32.

    [5] 李永樹,王金莊.背斜構(gòu)造礦層開采地面沉陷規(guī)律研究[J].巖土工程學(xué)報,1998,20(6):67-69.

    [6] 袁仁茂,馬鳳山,鄧清海.急傾斜厚大金屬礦山地下開挖巖移發(fā)生機(jī)理[J]中國地質(zhì)災(zāi)害與防治學(xué)報,2008,19(1):62-67.

    [7] 尹宏杰,朱建軍,李志偉,等.基于SBAS的礦區(qū)形變研究[J].測繪學(xué)報,2011,40(1):52-58.

    [8] 唐佑輝,黃騰,趙兵帥.GOM和時間序列組合模型在地表沉降預(yù)測中的應(yīng)用[J].測繪工程,2014,23(2):70-72.

    [9] 王行風(fēng),汪云甲,杜培軍.利用差分干涉測量技術(shù)監(jiān)測煤礦區(qū)開采沉陷變形的初步研究[J].中國礦業(yè),2007,16(7):76-79.

    [10]朱建軍,邢學(xué)敏,胡俊,等.利用InSAR技術(shù)監(jiān)測礦區(qū)地表形變[J].中國有色金屬報,2011,21(10):2564-2576.

    猜你喜歡
    中段采空區(qū)礦床
    老采空區(qū)建設(shè)場地采空塌陷地質(zhì)災(zāi)害及防治
    尿檢時如何取中段尿
    瞬變電磁法在煤礦采空區(qū)探測中的應(yīng)用
    構(gòu)造疊加暈法在深部找礦中的應(yīng)用——以河南小秦嶺楊砦峪金礦床S60號礦脈為例
    東天山中段晚古生代剪切帶疊加特征及構(gòu)造控礦作用
    黑龍江省林口縣三合村探明超大型石墨礦床
    邁好從低段到中段的“坎”
    考試周刊(2018年15期)2018-01-21 17:24:26
    西昆侖新發(fā)現(xiàn)鹽湖型鹵水硼鋰礦床
    某礦山采空區(qū)處理方案
    新疆鋼鐵(2015年2期)2015-11-07 03:27:52
    回風(fēng)井底附近采空區(qū)防滅火技術(shù)探討
    同煤科技(2015年2期)2015-02-28 16:59:15
    av网站在线播放免费| 99精国产麻豆久久婷婷| 久久人人97超碰香蕉20202| 中文字幕精品免费在线观看视频| 午夜免费鲁丝| 黄片无遮挡物在线观看| 一区二区三区激情视频| 亚洲内射少妇av| 男人添女人高潮全过程视频| 欧美日韩亚洲高清精品| 亚洲三区欧美一区| 日韩一区二区三区影片| videossex国产| √禁漫天堂资源中文www| 哪个播放器可以免费观看大片| 亚洲精品视频女| 老司机影院毛片| 久久精品熟女亚洲av麻豆精品| 久久午夜福利片| 久久久精品国产亚洲av高清涩受| 亚洲av中文av极速乱| 女性生殖器流出的白浆| 国产伦理片在线播放av一区| 国产日韩欧美视频二区| 丰满迷人的少妇在线观看| 91国产中文字幕| 精品人妻一区二区三区麻豆| 久久99热这里只频精品6学生| 久久综合国产亚洲精品| 成人国产av品久久久| 国产精品久久久久成人av| 这个男人来自地球电影免费观看 | 波多野结衣av一区二区av| 热re99久久国产66热| 国产亚洲欧美精品永久| 极品人妻少妇av视频| 午夜激情久久久久久久| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品麻豆人妻色哟哟久久| 蜜桃国产av成人99| 国产在视频线精品| 国产高清国产精品国产三级| 久久综合国产亚洲精品| 色视频在线一区二区三区| 日本色播在线视频| 成人毛片a级毛片在线播放| 免费黄频网站在线观看国产| 在线 av 中文字幕| 国产成人免费无遮挡视频| 黄色怎么调成土黄色| 欧美日韩成人在线一区二区| 精品人妻熟女毛片av久久网站| 中文字幕人妻丝袜一区二区 | 狂野欧美激情性bbbbbb| 爱豆传媒免费全集在线观看| 黄频高清免费视频| 高清在线视频一区二区三区| 久久ye,这里只有精品| 熟女少妇亚洲综合色aaa.| 香蕉国产在线看| 一级毛片黄色毛片免费观看视频| 欧美日韩综合久久久久久| 亚洲欧美一区二区三区黑人 | 国产综合精华液| 久久精品亚洲av国产电影网| 日韩一本色道免费dvd| 中文乱码字字幕精品一区二区三区| 国产探花极品一区二区| 夜夜骑夜夜射夜夜干| 99久久中文字幕三级久久日本| 晚上一个人看的免费电影| 欧美精品人与动牲交sv欧美| 国产成人精品久久久久久| 国产精品久久久av美女十八| 亚洲精品一二三| 国产乱人偷精品视频| 三级国产精品片| 9热在线视频观看99| 国产精品嫩草影院av在线观看| 中国国产av一级| 边亲边吃奶的免费视频| 满18在线观看网站| 老司机影院毛片| 在线看a的网站| av网站免费在线观看视频| 满18在线观看网站| 国产成人精品久久二区二区91 | 国产一区二区三区综合在线观看| 国产精品国产av在线观看| 日本91视频免费播放| 日韩欧美一区视频在线观看| 久热久热在线精品观看| 晚上一个人看的免费电影| 涩涩av久久男人的天堂| 精品亚洲成a人片在线观看| 国产乱人偷精品视频| 久久久久网色| 中文字幕精品免费在线观看视频| 水蜜桃什么品种好| 欧美 日韩 精品 国产| 免费看av在线观看网站| 性色av一级| 又黄又粗又硬又大视频| 男人爽女人下面视频在线观看| 晚上一个人看的免费电影| 黄色 视频免费看| 自线自在国产av| 超碰成人久久| 中文字幕色久视频| 永久免费av网站大全| 国产欧美亚洲国产| 亚洲第一青青草原| 日韩一区二区视频免费看| 国产精品久久久av美女十八| 18禁动态无遮挡网站| 久久久久久人人人人人| 欧美成人精品欧美一级黄| 午夜激情av网站| 久久久久国产精品人妻一区二区| 国产精品国产三级国产专区5o| 建设人人有责人人尽责人人享有的| 侵犯人妻中文字幕一二三四区| 麻豆精品久久久久久蜜桃| 一级黄片播放器| 国产片特级美女逼逼视频| 亚洲国产精品成人久久小说| 少妇人妻精品综合一区二区| 99久久精品国产国产毛片| 国产片特级美女逼逼视频| 激情视频va一区二区三区| 日韩中文字幕欧美一区二区 | 久久久久久久久久久久大奶| 亚洲av日韩在线播放| 亚洲激情五月婷婷啪啪| 黄片无遮挡物在线观看| 国产精品无大码| 麻豆av在线久日| 亚洲美女黄色视频免费看| 欧美日韩精品成人综合77777| 少妇人妻久久综合中文| 自拍欧美九色日韩亚洲蝌蚪91| 久久ye,这里只有精品| 免费观看在线日韩| 交换朋友夫妻互换小说| 哪个播放器可以免费观看大片| 亚洲国产欧美在线一区| freevideosex欧美| 国产有黄有色有爽视频| 成人亚洲欧美一区二区av| 国产毛片在线视频| 欧美激情极品国产一区二区三区| 老司机影院毛片| 在线天堂最新版资源| 久久久亚洲精品成人影院| 国产av精品麻豆| 久久精品久久久久久久性| 一边摸一边做爽爽视频免费| 国产精品国产三级国产专区5o| 美女福利国产在线| 一二三四中文在线观看免费高清| videossex国产| 久久韩国三级中文字幕| 午夜91福利影院| 一二三四在线观看免费中文在| 高清av免费在线| 97在线视频观看| 欧美少妇被猛烈插入视频| 国产片内射在线| 高清不卡的av网站| 色吧在线观看| 少妇的丰满在线观看| 中文精品一卡2卡3卡4更新| 美女主播在线视频| 日韩不卡一区二区三区视频在线| 七月丁香在线播放| 看非洲黑人一级黄片| 十分钟在线观看高清视频www| 国产精品免费大片| 青青草视频在线视频观看| 啦啦啦在线免费观看视频4| 色播在线永久视频| 国产片特级美女逼逼视频| 9191精品国产免费久久| 色婷婷久久久亚洲欧美| 在线 av 中文字幕| 久久久精品区二区三区| 丰满饥渴人妻一区二区三| 美女视频免费永久观看网站| 最近的中文字幕免费完整| 九九爱精品视频在线观看| 精品亚洲成国产av| 18+在线观看网站| 欧美黄色片欧美黄色片| 国产成人精品久久二区二区91 | 三上悠亚av全集在线观看| 在线观看人妻少妇| 国产精品一二三区在线看| 亚洲国产精品成人久久小说| 日韩 亚洲 欧美在线| 精品国产超薄肉色丝袜足j| 精品人妻一区二区三区麻豆| a级片在线免费高清观看视频| 嫩草影院入口| 久久久精品94久久精品| 日本爱情动作片www.在线观看| 中文字幕av电影在线播放| 9191精品国产免费久久| 我要看黄色一级片免费的| 美国免费a级毛片| 久久久精品区二区三区| 亚洲伊人色综图| 考比视频在线观看| av线在线观看网站| 国产精品蜜桃在线观看| 巨乳人妻的诱惑在线观看| 欧美老熟妇乱子伦牲交| 女人久久www免费人成看片| 永久网站在线| 欧美日韩成人在线一区二区| 久久久久视频综合| 国精品久久久久久国模美| 精品少妇久久久久久888优播| 午夜日韩欧美国产| 亚洲欧美清纯卡通| 精品少妇一区二区三区视频日本电影 | 少妇猛男粗大的猛烈进出视频| 色视频在线一区二区三区| 九色亚洲精品在线播放| 女的被弄到高潮叫床怎么办| 国产又色又爽无遮挡免| 美女脱内裤让男人舔精品视频| 久久久久久久久久久久大奶| 欧美日韩精品成人综合77777| 久久人人爽av亚洲精品天堂| 9热在线视频观看99| 香蕉精品网在线| 男女无遮挡免费网站观看| 久久国产精品男人的天堂亚洲| 午夜福利视频在线观看免费| h视频一区二区三区| 日本黄色日本黄色录像| 永久网站在线| 成人黄色视频免费在线看| 国产亚洲欧美精品永久| 春色校园在线视频观看| 亚洲欧美日韩另类电影网站| 国产一区有黄有色的免费视频| 人妻一区二区av| 91aial.com中文字幕在线观看| 国产熟女午夜一区二区三区| 可以免费在线观看a视频的电影网站 | 少妇被粗大的猛进出69影院| 国产无遮挡羞羞视频在线观看| 国产爽快片一区二区三区| 少妇精品久久久久久久| 日韩,欧美,国产一区二区三区| 亚洲国产精品一区二区三区在线| 美女主播在线视频| 美女xxoo啪啪120秒动态图| 精品少妇内射三级| 欧美精品亚洲一区二区| 亚洲精品一区蜜桃| 欧美黄色片欧美黄色片| 在线天堂最新版资源| av在线观看视频网站免费| 大码成人一级视频| 精品少妇内射三级| 国产精品免费视频内射| 欧美成人午夜免费资源| 亚洲视频免费观看视频| 欧美亚洲 丝袜 人妻 在线| 日韩制服骚丝袜av| av在线播放精品| 亚洲欧美中文字幕日韩二区| 男人舔女人的私密视频| 波野结衣二区三区在线| 丝袜人妻中文字幕| 高清av免费在线| 中文字幕精品免费在线观看视频| 国产精品国产三级国产专区5o| 日本欧美视频一区| 久久女婷五月综合色啪小说| 一区二区三区精品91| av网站在线播放免费| 色播在线永久视频| 97在线人人人人妻| 九草在线视频观看| 91国产中文字幕| 色视频在线一区二区三区| 久久久久国产精品人妻一区二区| 国产精品久久久av美女十八| 国产成人一区二区在线| 美女脱内裤让男人舔精品视频| 99久久中文字幕三级久久日本| 黑人猛操日本美女一级片| 2021少妇久久久久久久久久久| 国产在线一区二区三区精| 亚洲男人天堂网一区| 五月伊人婷婷丁香| 两个人看的免费小视频| 欧美精品亚洲一区二区| 久久亚洲国产成人精品v| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩精品网址| 午夜免费男女啪啪视频观看| 亚洲av成人精品一二三区| 国产精品嫩草影院av在线观看| 国产不卡av网站在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲国产最新在线播放| 国产免费又黄又爽又色| 国产精品无大码| 人妻人人澡人人爽人人| 中文字幕最新亚洲高清| 麻豆乱淫一区二区| 最近中文字幕2019免费版| 欧美日韩亚洲国产一区二区在线观看 | 天堂8中文在线网| 亚洲国产精品成人久久小说| a级片在线免费高清观看视频| 欧美日韩一级在线毛片| www.熟女人妻精品国产| 国产片特级美女逼逼视频| 国产毛片在线视频| 亚洲欧洲国产日韩| 日本欧美国产在线视频| 久久久a久久爽久久v久久| 中文天堂在线官网| 一级片'在线观看视频| 亚洲五月色婷婷综合| 日本黄色日本黄色录像| 精品国产一区二区久久| 中文字幕精品免费在线观看视频| 欧美av亚洲av综合av国产av | 国产精品一国产av| 亚洲色图综合在线观看| 亚洲图色成人| 日本免费在线观看一区| 一级片'在线观看视频| 一二三四在线观看免费中文在| 女人久久www免费人成看片| 国产极品粉嫩免费观看在线| 国产老妇伦熟女老妇高清| 亚洲第一青青草原| 亚洲国产av新网站| a级毛片黄视频| 欧美激情 高清一区二区三区| 欧美精品一区二区免费开放| 日日摸夜夜添夜夜爱| 久久这里有精品视频免费| 日本午夜av视频| 一二三四中文在线观看免费高清| 久久人人97超碰香蕉20202| 人成视频在线观看免费观看| 国产亚洲欧美精品永久| 国产精品久久久久久久久免| 国产精品.久久久| 亚洲综合色惰| 狂野欧美激情性bbbbbb| 国产成人精品福利久久| 亚洲综合色惰| 最近的中文字幕免费完整| 亚洲中文av在线| 午夜激情久久久久久久| 亚洲三区欧美一区| 黄色配什么色好看| 亚洲av.av天堂| 精品少妇一区二区三区视频日本电影 | 18禁裸乳无遮挡动漫免费视频| www日本在线高清视频| 国产成人精品一,二区| 亚洲av国产av综合av卡| 国产免费视频播放在线视频| 日本欧美视频一区| 国产黄频视频在线观看| 99久久综合免费| 久久婷婷青草| 欧美精品人与动牲交sv欧美| 黄片播放在线免费| 熟妇人妻不卡中文字幕| 久久久久久久久久人人人人人人| 欧美人与性动交α欧美精品济南到 | 欧美日韩国产mv在线观看视频| 国产午夜精品一二区理论片| 久久这里只有精品19| 亚洲精品成人av观看孕妇| 久久久国产精品麻豆| 一级毛片电影观看| 国产精品久久久久久精品古装| 亚洲成人av在线免费| 精品久久蜜臀av无| 欧美成人午夜精品| 国产成人精品福利久久| 国产精品亚洲av一区麻豆 | 建设人人有责人人尽责人人享有的| 国产男女超爽视频在线观看| 色哟哟·www| 伊人久久国产一区二区| 美女大奶头黄色视频| 日日摸夜夜添夜夜爱| 一级黄片播放器| 中文字幕人妻丝袜制服| 在线观看美女被高潮喷水网站| 丁香六月天网| 99香蕉大伊视频| 欧美精品亚洲一区二区| 日韩熟女老妇一区二区性免费视频| 国产亚洲欧美精品永久| 999久久久国产精品视频| 一级片免费观看大全| av天堂久久9| 巨乳人妻的诱惑在线观看| 国产一区二区在线观看av| 日韩中文字幕欧美一区二区 | 国产成人aa在线观看| 免费高清在线观看日韩| 黄色怎么调成土黄色| 一区福利在线观看| 如日韩欧美国产精品一区二区三区| 亚洲欧美中文字幕日韩二区| 国产白丝娇喘喷水9色精品| 男女高潮啪啪啪动态图| 极品少妇高潮喷水抽搐| 香蕉国产在线看| 久久久久网色| 欧美激情极品国产一区二区三区| 在线观看三级黄色| 午夜精品国产一区二区电影| a级毛片在线看网站| 性色avwww在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人av激情在线播放| 国产精品麻豆人妻色哟哟久久| 男人操女人黄网站| 在现免费观看毛片| 国产白丝娇喘喷水9色精品| 国产精品 欧美亚洲| 精品亚洲乱码少妇综合久久| 这个男人来自地球电影免费观看 | 一本大道久久a久久精品| 精品酒店卫生间| 一级毛片电影观看| 在线观看人妻少妇| 亚洲美女视频黄频| 一区二区三区激情视频| 亚洲欧美中文字幕日韩二区| 哪个播放器可以免费观看大片| 超碰成人久久| 黄色一级大片看看| 又粗又硬又长又爽又黄的视频| 久久人人爽人人片av| 熟女电影av网| 一区二区三区四区激情视频| 国产成人午夜福利电影在线观看| 99热网站在线观看| 国产精品久久久久久精品古装| 一级,二级,三级黄色视频| 亚洲婷婷狠狠爱综合网| 亚洲第一青青草原| 亚洲欧洲日产国产| 国产日韩欧美亚洲二区| 大香蕉久久成人网| 97人妻天天添夜夜摸| 欧美黄色片欧美黄色片| 国产成人精品婷婷| 久热久热在线精品观看| 美女大奶头黄色视频| 午夜日韩欧美国产| 十八禁网站网址无遮挡| 欧美精品国产亚洲| av在线app专区| 欧美日韩视频精品一区| 婷婷色综合www| 九草在线视频观看| 美女xxoo啪啪120秒动态图| 欧美在线黄色| 久久精品国产自在天天线| 在线免费观看不下载黄p国产| 精品酒店卫生间| 2018国产大陆天天弄谢| 超碰97精品在线观看| 美女大奶头黄色视频| 欧美日韩精品网址| kizo精华| 啦啦啦在线观看免费高清www| 亚洲图色成人| 69精品国产乱码久久久| 一级a爱视频在线免费观看| 黄频高清免费视频| 一区二区三区乱码不卡18| 免费在线观看完整版高清| av网站在线播放免费| 亚洲美女搞黄在线观看| 日产精品乱码卡一卡2卡三| 超碰成人久久| 人人妻人人澡人人爽人人夜夜| 青春草亚洲视频在线观看| 在线观看一区二区三区激情| 青草久久国产| 亚洲五月色婷婷综合| 国产福利在线免费观看视频| 在线观看美女被高潮喷水网站| 亚洲精品一二三| 丰满乱子伦码专区| 国产日韩欧美亚洲二区| 熟女av电影| 亚洲综合色网址| 久久精品国产鲁丝片午夜精品| 国产精品无大码| 三级国产精品片| 边亲边吃奶的免费视频| 亚洲国产欧美网| 亚洲久久久国产精品| 成人影院久久| 成年动漫av网址| 亚洲伊人久久精品综合| 丝瓜视频免费看黄片| 国产精品偷伦视频观看了| 午夜老司机福利剧场| 亚洲久久久国产精品| 丰满乱子伦码专区| 欧美在线黄色| 国产一区亚洲一区在线观看| av在线app专区| 日韩不卡一区二区三区视频在线| 国产精品久久久久久精品古装| 91午夜精品亚洲一区二区三区| 亚洲一区中文字幕在线| 欧美日韩综合久久久久久| 久久青草综合色| 精品少妇黑人巨大在线播放| 欧美少妇被猛烈插入视频| 老女人水多毛片| 久久这里只有精品19| 下体分泌物呈黄色| 黄频高清免费视频| av国产久精品久网站免费入址| 亚洲中文av在线| 高清不卡的av网站| 亚洲久久久国产精品| 国产人伦9x9x在线观看 | 欧美+日韩+精品| 黑人欧美特级aaaaaa片| 深夜精品福利| 天堂8中文在线网| 精品少妇黑人巨大在线播放| 观看美女的网站| 久久久久精品人妻al黑| 久久青草综合色| 亚洲男人天堂网一区| 男女国产视频网站| 蜜桃在线观看..| 欧美精品av麻豆av| av视频免费观看在线观看| 高清av免费在线| 国产日韩欧美在线精品| 精品久久蜜臀av无| 精品国产一区二区三区四区第35| 国产麻豆69| 精品人妻一区二区三区麻豆| 国产欧美日韩综合在线一区二区| 国产一级毛片在线| 日本av免费视频播放| 国产日韩欧美视频二区| 国产成人精品久久二区二区91 | av网站在线播放免费| 欧美97在线视频| 国产极品粉嫩免费观看在线| 熟妇人妻不卡中文字幕| 久久精品国产自在天天线| 天天躁夜夜躁狠狠躁躁| 一区二区av电影网| 亚洲伊人色综图| av免费在线看不卡| videos熟女内射| 91在线精品国自产拍蜜月| 亚洲,欧美精品.| 久久久久久久精品精品| 女性生殖器流出的白浆| 王馨瑶露胸无遮挡在线观看| 香蕉丝袜av| 男人爽女人下面视频在线观看| 亚洲精品aⅴ在线观看| 国产成人aa在线观看| 最新的欧美精品一区二区| 日韩一区二区三区影片| 国产一区二区三区av在线| 纯流量卡能插随身wifi吗| 免费观看a级毛片全部| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩综合在线一区二区| 18禁动态无遮挡网站| 国产麻豆69| 欧美av亚洲av综合av国产av | 国产成人精品久久久久久| 一本大道久久a久久精品| 欧美日韩一区二区视频在线观看视频在线| 黄色怎么调成土黄色| 不卡视频在线观看欧美| 美国免费a级毛片| 精品酒店卫生间| 不卡av一区二区三区| 2022亚洲国产成人精品| 精品99又大又爽又粗少妇毛片| 成年女人毛片免费观看观看9 | 午夜福利在线观看免费完整高清在| 国产精品久久久久成人av| 亚洲精品自拍成人| 另类精品久久| 国产一区二区在线观看av| 欧美av亚洲av综合av国产av | 大码成人一级视频| 一区福利在线观看| 黄色怎么调成土黄色| 90打野战视频偷拍视频|