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

    高潛水位采煤沉陷區(qū)積水范圍動態(tài)演化規(guī)律

    2020-06-05 02:11:34陳曉謝張文濤朱曉峻孫愛國楊曉玉
    煤田地質(zhì)與勘探 2020年2期
    關(guān)鍵詞:增長量采動積水

    陳曉謝,張文濤,朱曉峻,劉 輝,4,張 敏,孫愛國,楊曉玉

    高潛水位采煤沉陷區(qū)積水范圍動態(tài)演化規(guī)律

    陳曉謝1,2,張文濤3,朱曉峻1,2,劉 輝1,2,4,張 敏1,2,孫愛國5,楊曉玉6

    (1. 安徽大學(xué) 資源與環(huán)境工程學(xué)院,安徽 合肥 230601;2. 安徽省礦山生態(tài)修復(fù)工程實驗室,安徽 合肥 230601;3. 安徽理工大學(xué) 地球與環(huán)境學(xué)院,安徽 淮南 232001;4. 河北工程大學(xué) 礦業(yè)與測繪學(xué)院,河北 邯鄲 056038;5. 安徽恒源煤電股份有限公司五溝煤礦,安徽 淮北 235131;6. 合肥財經(jīng)職業(yè)學(xué)院,安徽 合肥 230601)

    為了研究高潛水位采煤沉陷區(qū)積水范圍動態(tài)演化規(guī)律,以安徽淮北五溝煤礦1031工作面為研究區(qū)域,基于Landsat-8遙感數(shù)據(jù)解譯結(jié)果,總結(jié)了2013—2017年1031工作面上方地表積水范圍演化規(guī)律,分析凈降水量、地下水埋深、工作面推進(jìn)距離等因素對積水范圍演化的影響,并提出積水邊界角的概念,建立了積水邊界角隨采動時間變化的函數(shù)關(guān)系式。研究發(fā)現(xiàn):高潛水位采煤沉陷區(qū)地表積水演化分為4個階段:未形成期、同步增長期、殘余增長期、相對穩(wěn)定期;在同步增長期和殘余增長期,工作面推進(jìn)距離、凈降水量與積水面積日增長量呈正相關(guān),地下水埋深與積水面積日增長量呈負(fù)相關(guān);工作面推進(jìn)距離是積水范圍演化的關(guān)鍵影響因素,在工作面推進(jìn)距離為476 m左右時,地表產(chǎn)生積水,工作面停采后5個月左右,積水面積趨于穩(wěn)定;積水邊界角先呈減小趨勢,再趨于穩(wěn)定,在同步增長期和殘余增長期,積水邊界角總體呈減小趨勢,但受氣象、地下水埋深等因素影響而波動;相對穩(wěn)定期,積水邊界角整體趨于90°。本研究為高潛水位采煤沉陷積水區(qū)土地利用規(guī)劃、土地復(fù)墾、水陸復(fù)合生態(tài)系統(tǒng)建立等提供理論依據(jù)。

    高潛水位;采煤沉陷區(qū);遙感監(jiān)測;積水演化階段;積水邊界角;五溝煤礦

    煤炭資源在我國能源體系中占非常重要的地位,占國家能源消費總量的64%[1],煤炭的大規(guī)模開采造成了嚴(yán)重的地表塌陷、地裂縫等地質(zhì)災(zāi)害[2-3]。兩淮地區(qū)是我國重要的煤炭基地,該地區(qū)地下潛水位高、松散層厚、下沉系數(shù)大,煤炭資源的開采引起地表下沉形成盆地,加上大氣降水與地下水補給等因素的共同影響,沉陷盆地內(nèi)形成大面積的積水,導(dǎo)致沉陷區(qū)水土流失嚴(yán)重,村莊、耕地被積水淹沒,原有的陸地生態(tài)系統(tǒng)轉(zhuǎn)化為水陸復(fù)合生態(tài)系統(tǒng)[4-5],對當(dāng)?shù)厣鷳B(tài)環(huán)境、社會經(jīng)濟造成了嚴(yán)重的影響。

    關(guān)于沉陷區(qū)積水演化研究,目前主要集中在積水面積演化和水循環(huán)兩個方面。在積水演化方面,彭蘇萍等[6]以淮南市為例,利用多時相TM遙感影像監(jiān)測地表積水塌陷區(qū)動態(tài)變化信息;魏礦靈等[7]以淮南礦區(qū)為例,通過監(jiān)測不同時期SPOT-5遙感影像,確定開采沉陷的位置及積水變化范圍;馬建威等[8]利用Landsat系列衛(wèi)星數(shù)據(jù),研究了武漢地區(qū)1973—2015年的城區(qū)水面范圍變化,分析多種因素對水體面積的影響等。在水循環(huán)方面,陸垂裕等[9-10]提出了面向?qū)ο竽K化的分布式水文模型MODCYCLE,并借助相關(guān)水文推理和數(shù)值模擬分析,對淮南典型高潛水位采煤沉陷區(qū)的地下水機制進(jìn)行辨析,得到高潛水位采煤沉陷區(qū)積水的補給來源和排泄渠道;李慧等[11-12]利用分布式水文模型,深入研究東部高潛水位采煤沉陷積水區(qū)庫容變化的主要影響因素及其相互關(guān)系。

    針對高潛水位采煤沉陷區(qū)積水演化規(guī)律研究,多針對地表或淺層地表的積水范圍,結(jié)合工作面采動影響的積水動態(tài)演化規(guī)律方面研究較少。本文以五溝煤礦1031工作面為研究區(qū),基于Landsat-8衛(wèi)星遙感數(shù)據(jù),對1031工作面地表積水范圍進(jìn)行長時序的動態(tài)監(jiān)測,提取積水邊界矢量并統(tǒng)計積水面積,分析積水范圍動態(tài)演化規(guī)律及其影響因素,以期為高潛水位采煤沉陷區(qū)土地利用規(guī)劃、土地復(fù)墾、水陸復(fù)合生態(tài)系統(tǒng)建立提供理論依據(jù)。

    1 研究區(qū)概況

    1031工作面位于安徽省濉溪縣五溝煤礦,是典型的高潛水位采煤區(qū),該地區(qū)降水量豐富,年均降水量827 mm,地下水位埋深1.5~4.0 m,多年平均埋深2.48 m[13]。1031工作面采動時間為2013年10月至2015年3月,工作面所在位置及現(xiàn)狀如圖1所示,底圖為2019年11月11日研究區(qū)遙感圖像,其中1031工作面全長1 411 m,寬180 m,煤層傾角為3°~16°,平均傾角8°,煤層平均厚度3.8 m,平均埋深(0)約364 m。工作面上覆巖層以砂巖和泥巖為主,厚度約94 m,松散層厚度約270 m,鉆孔柱狀如圖2所示。工作面所在區(qū)域地表地勢平坦,標(biāo)高+26.50~+27.56 m,無大流量河流經(jīng)過。由于該地區(qū)地下水埋深淺,降水量大[14],且地表下沉系數(shù)大,工作面采動后地表很快將形成積水[15]。

    圖1 五溝煤礦研究區(qū)所在位置及現(xiàn)狀

    2 積水范圍遙感監(jiān)測

    礦區(qū)地表積水周期長、范圍大,且所處位置交通不便,常規(guī)的觀測手段難以獲取大范圍、長時序、完整統(tǒng)一的積水范圍演化數(shù)據(jù);遙感觀測手段具有觀測周期短、數(shù)據(jù)量大、精度高等特點,因此,選擇遙感資料作為數(shù)據(jù)源,提取采煤沉陷區(qū)地表積水范圍信息[16]。本次選用Landsat-8系列衛(wèi)星的遙感數(shù)據(jù),影像空間分辨率為30 m,衛(wèi)星回歸周期為16 d。

    2.1 遙感圖像獲取與預(yù)處理

    五溝煤礦1031工作面采動時間為2013年10月至2015年3月,地表積水盆地不受其他工作面采動影響的截止時間為2016年10月,故選取2013—2016年Landsat-8獲取的遙感影像。剔除云量較大、能見度低的遙感影像,最終選取16期滿足水體提取基本要求的遙感影像。對選取的每幅影像進(jìn)行預(yù)處理,步驟包括:輻射定標(biāo)、大氣校正和圖像裁剪,得到最終影像[17]。

    圖2 五溝煤礦1031工作面鉆孔柱狀

    2.2 積水范圍提取

    目前,較高效的水體提取方法為水體指數(shù)法,常用的水體指數(shù)包括歸一化差異水體指數(shù)NDWI(Normalized Difference Water Index)[18]、改進(jìn)后的歸一化差異水體指數(shù)MNDWI[19](Modified NDWI)及自動水體提取指數(shù)AWEI(Automated Water Extrac-tion Index)[20]。NDWI利用近紅外波段和綠波段的反射率差值提取水體,該方法提取時可能有部分不透水層像素會被誤判為水體;故提出了MNDWI指數(shù),將NDWI中的近紅外波段換成短波紅外波段,增加水體的提取精度[19];AWEI針對Landsat-5影像設(shè)計,通過最大化水體和非水體區(qū)分度而得到最優(yōu)的系數(shù)組合[21]。結(jié)合研究區(qū)條件,對比各種水體指數(shù)特點,選擇MNDWI指數(shù)用于該研究區(qū)的水體提取。

    選擇Landsat-8影像中的綠波段(Band 3)和短波紅外波段(Band 6)反射值進(jìn)行計算,公式如下:

    式中:Green和MIR分別為綠光波段和短波紅外波段的反射值。

    通過MNDWI指數(shù)對遙感圖像進(jìn)行處理,根據(jù)影像內(nèi)的水體特征,統(tǒng)計積水面積并將積水邊界轉(zhuǎn)成矢量文件輸出。如圖3所示,底圖為Landsat-8標(biāo)準(zhǔn)彩色合成影像,白框部分為同期研究區(qū)積水矢量邊界,積水矢量邊界與水體邊界基本吻合,可知,通過MNDWI提取的積水邊界與水體邊界基本吻合,具有較高的準(zhǔn)確性,能夠滿足研究需求。

    圖3 2016年9月2日積水矢量邊界提取效果

    3 積水動態(tài)演化規(guī)律

    通過MNDWI方法,得到16幅積水范圍矢量并統(tǒng)計積水面積,結(jié)合工作面推進(jìn)位置,監(jiān)測積水范圍動態(tài)變化,分析地表積水范圍演化規(guī)律。

    圖4為1031工作面推進(jìn)位置與地表積水范圍變化疊加圖像,由圖4可以看出,地表積水形成明顯滯后于地下工作面采動時間,當(dāng)工作面推進(jìn)476 m后,即推進(jìn)距離接近采深(0)時,地表積水才開始形成;工作面推進(jìn)過程中,積水范圍演化方向與工作面推進(jìn)方向一致,隨著工作面推進(jìn),積水范圍沿工作面推進(jìn)方向擴大,積水面積日增長量約440 m2/d;2015年3月3日(工作面停止采動)至2015年7月30日,積水范圍仍然沿工作面推進(jìn)方向擴大,積水面積日增長量約120 m2/d;工作面停采后5個月左右,積水范圍不再擴大,趨于穩(wěn)定。

    為更方便研究積水范圍演化過程中的影響因素,通過工作面采動狀態(tài)及積水范圍的演化規(guī)律分析,將高潛水位采煤沉陷區(qū)積水范圍演化劃分為4個時期,如圖5所示。

    a. 未形成期 該階段為工作面開始采動至地表積水開始形成期間,地表積水尚未形成;

    圖4 2013—2016年工作面推進(jìn)距離與地表積水范圍變化

    圖5 五溝煤礦1031工作面積水面積變化

    b.同步增長期 該階段為地表首次觀測到積水至工作面停止采動期間,工作面處于采動狀態(tài),地表積水范圍隨工作面采動同步增長;

    c. 殘余增長期 該階段為工作面停止采動至停采后5個月左右,地表受采動殘余沉降影響仍呈動態(tài)變化,積水范圍沿工作面推進(jìn)方向呈增長趨勢,積水面積呈緩慢增長趨勢;

    d. 相對穩(wěn)定期 該階段為工作面停采5個月后,地表基本趨于穩(wěn)定,積水邊界不再沿采動方向擴大,積水面積整體趨于穩(wěn)定。

    4 地表積水演化的關(guān)鍵影響因素

    沉陷區(qū)積水范圍演化中,首先受工作面采動影響產(chǎn)生地表沉陷;另外,大氣降水和地下水補給是沉陷區(qū)積水的主要來源,占積水區(qū)水量來源90%以上,水分蒸發(fā)是主要的排泄方式[10],其他影響因素對積水范圍影響較小,非主要影響因素。

    為研究主要影響因素對積水范圍的作用,首先,根據(jù)2013—2016年國家氣象觀測數(shù)據(jù)計算研究區(qū)的降水量、蒸發(fā)量[22]及兩者差值,得到凈降水量,分析凈降水量對積水范圍演化各階段的影響;然后,根據(jù)濉溪縣地下水位埋深觀測數(shù)據(jù),分析地下水位埋深對積水演化各時期的影響;再次,結(jié)合工作面推進(jìn)距離,分析工作面推進(jìn)距離對積水范圍演化各時期的影響;最后,根據(jù)上述分析得到積水范圍演化的關(guān)鍵影響因素。

    4.1 凈降水量

    大氣降水、蒸發(fā)作用同步進(jìn)行,均對積水范圍產(chǎn)生影響,因此,應(yīng)綜合分析兩者對沉陷區(qū)積水范圍演化的作用。通過計算研究時期內(nèi)日降水量、日蒸發(fā)量差值之和,得到研究區(qū)凈降水量累計值。圖6為自首次觀測到積水以后的凈降水量與積水面積變化圖。

    圖6 凈降水量與積水面積相關(guān)關(guān)系

    積水未形成期,該沉陷區(qū)未產(chǎn)生積水。同步增長期和殘余增長期,積水面積整體呈增長趨勢,凈降水量與積水面積日增長量呈正相關(guān),隨著凈降水量增大,積水面積日增長量越大;凈降水量減小,積水面積日增長量相對減小。據(jù)圖6可知,2014年5月8日至2014年12月2日,凈降水量由–57 mm左右先減小至–170 mm左右,再增長到200 mm左右,同期積水面積由1.8 hm2增長到14.0 hm2,平均日增長量約584 m2;2014年12月2日至2015年3月3日(工作面停止采動)期間,凈降水量由200 mm減小至105 mm左右,同期積水面積由14.0 hm2增長到15.3 hm2,平均日增長量約116 m2,積水面積日增長量較前者明顯減小。相對穩(wěn)定期,凈降水量與積水面積日增長量負(fù)相關(guān),如2015年7月30日至2016年3月10日,凈降水量由180 mm減小至75 mm,由于凈降水量整體為正,對積水仍呈補充狀態(tài),同期積水面積呈小幅度增長,平均日增長量約為72 m2;2016年6月4日至2016年8月1日,凈降水量整體由15 mm增長至39 mm,但過程中凈降水量出現(xiàn)負(fù)值,同期積水面積呈小幅度減小,平均日減小量約170 m2。

    4.2 地下水位埋深

    淺層地下水在沉陷區(qū)形成積水后,與地表水體產(chǎn)生水力聯(lián)系,對地表水有一定的補充作用,是沉陷區(qū)積水主要來源之一[23],圖7為研究時期地下水位埋深與積水面積變化圖。

    圖7 地下水埋深與積水面積相關(guān)關(guān)系

    積水未形成期,地表未觀測到積水,該階段地表下沉量較小,未與地下水產(chǎn)生水力聯(lián)系;同步增長期與殘余增長期,地表下沉量變大,沉陷區(qū)形成積水區(qū)域,該階段地下水位埋深與積水面積日增加量呈負(fù)相關(guān),地下水位埋深快速減小,積水面積日增長量增大;地下水位埋深增加,地下水對沉陷區(qū)補給減少,積水面積日增長量減小。據(jù)圖7可知,2014年3月21日至2014年12月2日,地下水位埋深由3.7 m 左右減小到2.4 m左右,同期積水面積平均日增長量約504 m2;2015年3月24日至2015年7月30日,地下水位埋深由2.9 m左右減小至1.8 m,積水面積平均日增長量為139 m2;2014年12月2日至2015年3月24日,地下水位埋深由2.4 m左右增大至2.9 m,同期凈降水量由約210 mm減小到115 mm左右,積水面積平均日增長量約116 m2,較前兩區(qū)間增長速度明顯減??;相對穩(wěn)定期,沉陷區(qū)基本穩(wěn)定,地下水與沉陷區(qū)積水已產(chǎn)生穩(wěn)定水力聯(lián)系,淺層地下水為積水區(qū)域提供涵養(yǎng)條件[12],維持地表水的相對穩(wěn)定,該階段地下水位埋深與積水范圍演化相關(guān)性不明顯。

    4.3 工作面推進(jìn)距離

    工作面采動后,采區(qū)周圍巖體原始應(yīng)力平衡被破壞,在應(yīng)力達(dá)到重新平衡的過程中,上覆巖層產(chǎn)生連續(xù)的移動、變形和非連續(xù)破壞,地表產(chǎn)生移動盆地。當(dāng)移動盆地達(dá)到一定深度且與地下水產(chǎn)生水力聯(lián)系時,地表將形成積水區(qū)域,隨著地下工作面采動,地表積水區(qū)域隨之同步擴大。由積水范圍演化規(guī)律可知,地表積水范圍沿著工作面采動方向演化[24],可見工作面推進(jìn)距離與積水范圍演化具有明顯的相關(guān)性。

    圖8為工作面推進(jìn)距離與積水面積變化關(guān)系圖。據(jù)圖8可知,地表積水形成明顯滯后于工作面采動,在積水未形成期,受工作面采動影響巖層發(fā)生破壞,并逐漸傳遞至地表,地表形成下沉盆地,但下沉量較小,地表尚未形成穩(wěn)定積水區(qū)域;同步增長期,隨著下沉量增大,地表移動盆地已形成積水,截至工作面停止采動,積水面積隨工作面采動增長至15.3 hm2左右;殘余增長期,工作面停采后,采動導(dǎo)致的巖體移動尚未穩(wěn)定,該階段地表積水面積仍有小范圍增加,由15.3 hm2增長至17.1 hm2;相對穩(wěn)定期,該階段巖體應(yīng)力平衡基本穩(wěn)定,沉陷區(qū)地表趨于穩(wěn)定,但受降水量、蒸發(fā)量、地下水等因素影響,積水面積仍然存在小幅度波動。

    圖8 工作面推進(jìn)距離與地表積水面積相關(guān)關(guān)系

    4.4 關(guān)鍵影響因素確定

    綜上可知,積水未形成期,地表主要受工作面采動影響形成下沉盆地,由于下沉較小,地表尚未形成積水;同步增長期,地表下沉量增大,地表形成積水,該階段凈降水量、采動距離與積水面積日增長量呈正相關(guān),地下水埋深與積水面積日增長量呈負(fù)相關(guān);殘余增長期,該階段工作面停止采動,受殘余沉降影響,地表仍然有少量下沉,該階段積水面積仍然呈上升趨勢,凈降水量與積水面積日增長量呈正相關(guān),地下水埋深與積水面積日增長量呈負(fù)相關(guān);相對穩(wěn)定期,該階段地表下沉盆地基本穩(wěn)定,受降水量、蒸發(fā)量和地下水埋深的波動,積水區(qū)面積整體穩(wěn)定呈小范圍波動。

    凈降水量、地下水埋深、工作面推進(jìn)距離均與地表積水范圍演化有明顯相關(guān)性,但凈降水量、地下水位埋深受氣候條件影響均呈周期性變化[25],根據(jù)積水演化規(guī)律分析可知,工作面推進(jìn)距離增加為積水范圍演化的關(guān)鍵影響因素。

    5 工作面推進(jìn)距離與積水邊界角的關(guān)系

    為量化工作面采動距離與積水范圍演化的關(guān)系,提出積水邊界角的概念。工作面推進(jìn)距離是積水范圍演化的關(guān)鍵影響因素,與積水面積變化呈正相關(guān)。根據(jù)積水邊界角隨工作面采動距離的變化,建立工作面推進(jìn)方向上推進(jìn)位置與地表積水范圍之間的動態(tài)關(guān)系。

    如圖9所示,工作面推進(jìn)位置,工作面推進(jìn)方向剖面正上方積水邊界,連接,與工作面推進(jìn)方向的夾角定義為積水邊界角。圖中,在水平面投影距離為,2為工作面到積水水面的高差,0為采深,1為積水水面到原地表的距離,由于2遠(yuǎn)大于1,故2約等于0。

    圖9 積水邊界角示意圖

    積水邊界角計算公式為:

    式中:為積水邊界角,(°);0為采深,m;為在水平面的投影距離,m。

    根據(jù)式(2),計算每個觀測時間段積水邊界角。2013年10月至2014年3月21日地表未形成積水,積水邊界角記為0,從2014年3月21日開始統(tǒng)計,并利用統(tǒng)計數(shù)據(jù)繪制積水邊界角變化曲線,如圖10所示。

    圖10 積水邊界角變化曲線

    為便于研究積水演化周期內(nèi)積水邊界角的變化特征,根據(jù)工作面推進(jìn)速度與推進(jìn)距離,建立推進(jìn)距離與時間的關(guān)系,最終建立整個演化過程中積水邊界角與時間的關(guān)系。

    a. 未形成期 地表未形成積水,積水邊界角記為0。

    c. 殘余增長期 建立積水邊界角與工作面停采時間的關(guān)系。

    d. 相對穩(wěn)定期 該階段積水范圍整體動態(tài)平衡,積水邊界角小范圍波動,趨于90°。

    積水演化各階段中,積水邊界角與工作面采動時間的關(guān)系式為:

    由圖10、式(3)可得積水演化周期內(nèi)積水邊界角與工作面采動時間的關(guān)系。積水未形成期,不觀測積水邊界角;同步增長期和殘余增長期,積水邊界角整體呈減小趨勢,其中,同步增長期積水邊界角與時間呈三次函數(shù)關(guān)系,積水邊界角整體隨時間變化呈減小趨勢,即積水邊界角隨著工作面采動距離的增大而減??;殘余增長期,積水邊界角與采動時間呈一次函數(shù)關(guān)系,積水邊界角隨工作面停采時長的增加而減小;相對穩(wěn)定期,積水邊界角趨于穩(wěn)定,接近90°。

    6 結(jié)論

    a.高潛水位采煤沉陷區(qū)地表積水的形成分為4個階段:未形成期,開始采動至采動距離476 m;同步增長期,采動距離476 m至工作面停采;殘余增長期,工作面停采至停采后5個月左右;相對穩(wěn)定期,停采5個月后。

    b.積水范圍演化受到大氣降水、蒸發(fā)、地下水位埋深、工作面采動等因素的影響,在同步增長和殘余增長期,凈降水量、工作面采動距離與積水面積日增長量呈正相關(guān),地下水位埋深與積水面積日增長量呈負(fù)相關(guān);工作面推進(jìn)距離為積水范圍演化的關(guān)鍵影響因素。

    c.積水邊界角在同步增長期和殘余增長期整體呈減小趨勢;同步增長期的積水邊界角隨工作面推進(jìn)距離的增加而減??;殘余增長期的積水邊界角隨工作面停采時長的增加而減??;積水穩(wěn)定期,積水邊界角趨于穩(wěn)定,接近90°。

    請聽作者語音介紹創(chuàng)新技術(shù)成果等信息,歡迎與作者進(jìn)行交流

    [1] 英國石油公司.《BP世界能源統(tǒng)計年鑒》[EB/OL]. (2019- 06-14)[2019-10-10] https://www.bp.com/zh_cn/china/home/news/ reports.html. British Petroleum Company.《BP Statistical Review of World Energy》[EB/OL]. (2019-06-14)[2019-10-10] https://www.bp. com/zh_cn/china/home/news/reports.html.

    [2] 劉輝,鄧喀中,雷少剛,等. 采動地裂縫動態(tài)發(fā)育規(guī)律及治理標(biāo)準(zhǔn)探討[J]. 采礦與安全工程學(xué)報,2017,34(5):884–890. LIU Hui,DENG Kazhong,LEI Shaogang,et al. Dynamic developing law and governance standard of ground fissures caused by underground mining[J]. Journal of Mining & Safety Engineering,2017,34(5):884–890.

    [3] LIU Hui,DENG Kazhong,ZHU Xiaojun,et al. Effects of mining speed on the developmental features of mining-induced ground fissures[J]. Bulletin of Engineering Geology and the Environment,2019,78(8):6297–6309.

    [4] ZHU Xiaojun,GUO Guangli,ZHA Jianfeng,et al. Surface dynamic subsidence prediction model of solid backfill mining[J]. Environmental Earth Sciences,2016,75(6):1–9.

    [5] 孔令健,姜春露,鄭劉根,等. 淮北臨渙礦采煤沉陷區(qū)不同水體水化學(xué)特征及其影響因素[J]. 湖泊科學(xué),2017,29(5):1158–1167. KONG Lingjian,JIANG Chunlu,ZHENG Liugen,et al. Characters of hydrochemistry and their influenced factors of different waters in the Linhuan coal mining subsidence area of Huaibei City[J]. Journal of Lake Sciences,2017,29(5):1158–1167.

    [6] 彭蘇萍,王磊,孟召平,等. 遙感技術(shù)在煤礦區(qū)積水塌陷動態(tài)監(jiān)測中的應(yīng)用:以淮南礦區(qū)為例[J]. 煤炭學(xué)報,2002,27(4):374–378. PENG Suping,WANG Lei,MENG Zhaoping,et al. Monitoring the seeper subside in coal district by the remote sensing:Examples from Huainan coal district[J]. Journal of China Coal Society,2002,27(4):374–378.

    [7] 魏礦靈,王啟春,郭廣禮,等. 利用遙感數(shù)據(jù)監(jiān)測礦區(qū)開采沉陷積水變化[J]. 煤礦安全,2014,45(1):13–16. WEI Kuangling,WANG Qichun,GUO Guangli,et al. Using remote sensing data to monitor the change of mining subsidence water accumulation[J]. Safety in Coal Mines,2014,45(1):13–16.

    [8] 馬建威,黃詩峰,許宗男. 基于遙感的1973—2015年武漢市湖泊水域面積動態(tài)監(jiān)測與分析研究[J]. 水利學(xué)報,2017,48(8):903–913. MA Jianwei,HUANG Shifeng,XU Zongnan. Satellite remote sensing of lake area in Wuhan from 1973 to 2015[J]. Journal of Hydraulic Engineering,2017,48(8):903–913.

    [9] 陸垂裕,秦大庸,張俊娥,等. 面向?qū)ο竽K化的分布式水文模型MODCYCLEⅠ:模型原理與開發(fā)篇[J]. 水利學(xué)報,2012,43(10):1135–1145. LU Chuiyu,QIN Dayong,ZHANG Jun’e,et al. MODCYCLE-an object oriented modularized hydrological model Ⅰ:Theory and development[J]. Journal of Hydraulic Engineering,2012,43(10):1135–1145.

    [10] 陸垂裕,陸春輝,李慧,等. 淮南采煤沉陷區(qū)積水過程地下水作用機制[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(10):122–131. LU Chuiyu,LU Chunhui,LI Hui,et al. Groundwater influencing mechanism in water-ponding process of Huainan coal mining subsidence area[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2015,31(10):122–131.

    [11] 李慧,陸垂裕,孫青言,等. 分布式“河道–沉陷區(qū)–地下水”水循環(huán)耦合模型—Ⅰ模型原理與開發(fā)[J]. 水科學(xué)進(jìn)展,2016,27(2):214–223.LI Hui,LU Chuiyu,SUN Qingyan,et al. Research on the distributed “river-subsidence area-groundwater” coupling simulation of coal mining subsidence areaⅠ:Theory and development[J]. Advances in Water Science,2016,27(2):214–223.

    [12] 李慧. 采煤沉陷區(qū)分布式水循環(huán)模型研究[D]. 北京:中國水利水電科學(xué)研究院,2016. LI Hui. Research on the distributed water cycle model of coal mining subsidence area[D]. Beijing:China Institute of Water Resources and Hydropower Research,2016.

    [13] 王發(fā)信,柏菊. 淮北平原淺層地下水埋深區(qū)域分布特點[J]. 地下水,2014,36(5):51–53. WANG Faxin,BO Ju. Distribution characteristics of shallow groundwater depth in Huaibei plain[J]. Ground Water,2014,36(5):51–53.

    [14] 喬叢林,史明禮,蘇婭,等. 淮北平原地區(qū)水文特征[J]. 水文,2000,20(3):55–58. QIAO Conglin,SHI Mingli,SU Ya,et al. Hydrological features of Huaibei plain[J]. Journal of China Hydrology,2000,20(3):55–58.

    [15] 徐良驥,嚴(yán)家平,高永梅. 安徽省兩淮礦區(qū)地質(zhì)環(huán)境治理技術(shù)[J]. 煤田地質(zhì)與勘探,2007,35(6):37–40. XU Liangji,YAN Jiaping,GAO Yongmei. Comprehensive treatment technique of geological environment in Huainan & Huaibei mining area[J]. Coal Geology & Exploration,2007,35(6):37–40.

    [16] 王海波,馬明國. 基于遙感的湖泊水域動態(tài)變化監(jiān)測研究進(jìn)展[J]. 遙感技術(shù)與應(yīng)用,2009,24(5):674–684. WANG Haibo,MA Mingguo. A review of monitoring change in lake water areas based on remote sensing[J]. Remote Sensing Technology and Application,2009,24(5):674–684.

    [17] 祝佳. Landsat8衛(wèi)星遙感數(shù)據(jù)預(yù)處理方法[J]. 國土資源遙感,2016,28(2):21–27. ZHU Jia. Analysis of Landsat8 satellite remote sensing data preprocessing[J]. Remote Sensing for Land & Resources,2016,28(2):21–27.

    [18] MCFEETERS S K. The use of the Normalized Difference Water Index(NDWI) in the delineation of open water features[J]. International Journal of Remote Sensing,1996,17(7):1425–1432.

    [19] XU Hanqiu. Modification of Normalised Difference Water Index(NDWI) to enhance open water features in remotely sensed imagery[J]. International Journal of Remote Sensing,2006,27(14):3025–3033.

    [20] FEYISA G L,MEILBY H,F(xiàn)ENSHOLT R,et al. Automated water extraction index:A new technique for surface water mapping using Landsat imagery[J]. Remote Sensing of Environment,2014,140:23–35.

    [21] 計璐艷,尹丹艷,宮鵬. Landsat長時間序列的陽澄湖湖面圍網(wǎng)時空變化[J]. 遙感學(xué)報,2019,23(4):717–729.JI Luyan,YIN Danyan,GONG Peng. Temporal-spatial study on enclosure culture area in Yangcheng lake with long-term Landsat time series[J]. Journal of Remote Sensing,2019,23(4):717–729.

    [22] 張晶,李妍清. 基于地統(tǒng)計學(xué)的柬埔寨年降水量數(shù)據(jù)空間插值[J]. 人民長江,2018,49(22):100–103. ZHANG Jing,LI Yanqing. Spatial interpolation of annual precipitation data in Cambodia based on geostatistics[J]. Yangtze River,2018,49(22):100–103.

    [23] 張磊,秦小光,劉嘉麒,等. 淮南采煤沉陷區(qū)積水來源的氫氧穩(wěn)定同位素證據(jù)[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版),2015,45(5):1502–1514. ZHANG Lei,QIN Xiaoguang,LIU Jiaqi,et al. Characters of hydrogen and oxygen stable isotope of different water bodies in Huainan coal mining area[J]. Journal of Jilin University(Earth Science Edition),2015,45(5):1502–1514.

    [24] 何國清,楊倫,凌賡娣,等. 礦山開采沉陷學(xué)[M]. 徐州:中國礦業(yè)大學(xué)出版社,1994:84–92. HE Guoqing,YANG Lun,LING Gengti,et al. Mining subsidence[M]. Xuzhou:China University of Mining and Technology Press,1994:84–92.

    [25] 陸小明,楊秒,王振龍,等. 淮北平原氣象因素對裸地潛水蒸發(fā)的影響[J]. 灌溉排水學(xué)報,2019,38(4):84–91. LU Xiaoming,YANG Miao,WANG Zhenlong,et al. Metrological factors affecting evaporation of shallow groundwater in the absence of plants in Huaibei plain[J]. Journal of Irrigation and Drainage,2019,38(4):84–91.

    Dynamic evolution law of water accumulation range in coal mining subsidence area with high-level groundwater

    CHEN Xiaoxie1,2, ZHANG Wentao3, ZHU Xiaojun1,2, LIU Hui1,2,4, ZHANG Min1,2, SUN Aiguo5, YANG Xiaoyu6

    (1. College of Resources and Environmental Engineering, Anhui University, Hefei 230601, China; 2. Anhui Province Engineering Laboratory for Mine Ecological Remediation, Hefei 230601, China; 3. School of Earth and Environment, Anhui University of Science & Technology, Huainan 232001, China; 4. School of Mining and Geomatics, Hebei University of Engineering, Handan 056038, China; 5. Wugou Coal Mine of Anhui Hengyuan Coal Power Co. Ltd.,Huaibei 235131, China; 6. Hefei College of Finance & Economics, Hefei 230601, China)

    In order to study the dynamic evolution law of the water accumulation range with high-level groundwater, the working face 1031 of Wugou coal mine was taken as the research area.Based on Landsat-8 remote sensing data,the evolution law of water accumulation range above working face 1031 during 2013-2017 was summarized and the influence of factors such as net rainfall, groundwater depth and working surface propulsion distance on the evolution of water accumulation range was analyzed.The concept of the boundary angle of the water accumulation range was proposed, and the functional relationship between the angle and the mining time was established.The study found that the evolution of surface water in the coal mining subsidence area was divided into four periods: unformed, synchronous growth, residual growth, stable. The groundwater depth is negatively correlated with the daily growth of the water accumulation range. The working face advancement distance is the driving factor for the change of the water accumulation range. When the working distance is about 476 m, water is accumulated on the surface, about 5 months after stopping mining, the area of water accumulation range tends to be stable, in the synchronous growth and residual growth period, the angle generally decreases, in the stable period the angle tends to be 90°as a whole, but it fluctuates due to other factors.This study provides a theoretical basis for land use planning, land reclamation, and establishment of a land-water complex ecosystem in the coal mining subsidence area with high-level groundwater.

    high-level groundwater; coal mining subsidence area; remote sensing monitoring; evolution of subsidence water accumulation range; the boundary angle of water accumulation range; Wugou coal mine

    TD803

    A

    10.3969/j.issn.1001-1986.2020.02.020

    1001-1986(2020)02-0126-08

    2019-10-18;

    2020-02-22

    國家自然科學(xué)基金項目(51804001);安徽省自然科學(xué)基金項目(1808085QE147);安徽省高校自然科學(xué)研究重點項目(KJ2018A0003)

    National Natural Science Foundation of China(51804001);Natural Science Foundation of Anhui Province, China(1808085QE147);Key Project of Natural Science Research of Universities in Anhui Province(KJ2018A0003)

    陳曉謝,1993年生,男,安徽合肥人,碩士研究生,從事礦山開采沉陷方向研究. E-mail:chenxiaoxie99@163.com

    劉輝,1982年生,男,山東肥城人,博士,副教授,從事礦山開采沉陷及生態(tài)環(huán)境治理等方面的教學(xué)與科研工作. E-mail:lhui99@aliyun.com

    陳曉謝,張文濤,朱曉峻,等. 高潛水位采煤沉陷區(qū)積水范圍動態(tài)演化規(guī)律[J]. 煤田地質(zhì)與勘探,2020,48(2):126–133.

    CHEN Xiaoxie,ZHANG Wentao,ZHU Xiaojun,et al. Dynamic evolution law of water accumulation range in coal mining subsidence area with high-level groundwater[J]. Coal Geology & Exploration,2020,48(2):126–133.

    (責(zé)任編輯 周建軍)

    猜你喜歡
    增長量采動積水
    原來是輸卵管積水惹的禍
    不同配比的基質(zhì)對不同多肉植物生長的影響
    小熊當(dāng)當(dāng)玩積水
    不同坡度和不同坡向?qū)鸹ú枭L量的影響
    采動影響下淺埋輸氣管道與土體耦合作用機理
    原來是輸卵管積水惹的禍
    金都1號火龍果果實發(fā)育過程中品質(zhì)的變化規(guī)律分析
    施用新型肥料對閩南桉樹生長的影響
    深部采動巷道頂板穩(wěn)定性分析與控制
    采動巖體滲流力學(xué)研究進(jìn)展
    亚洲男人天堂网一区| 搡老妇女老女人老熟妇| 亚洲 欧美 日韩 在线 免费| 国产午夜福利久久久久久| 国产爱豆传媒在线观看 | 国产视频一区二区在线看| 悠悠久久av| 亚洲真实伦在线观看| 天堂动漫精品| 一二三四社区在线视频社区8| 午夜福利免费观看在线| 18禁裸乳无遮挡免费网站照片 | 色综合欧美亚洲国产小说| 亚洲欧美激情综合另类| 精品久久久久久久久久久久久 | 身体一侧抽搐| 神马国产精品三级电影在线观看 | 欧美性猛交黑人性爽| 日韩av在线大香蕉| 搡老妇女老女人老熟妇| 一个人观看的视频www高清免费观看 | 午夜激情av网站| 一区二区三区激情视频| 日韩成人在线观看一区二区三区| 男女床上黄色一级片免费看| 午夜福利免费观看在线| 亚洲av熟女| 久久香蕉激情| 18禁黄网站禁片午夜丰满| 成人三级做爰电影| 亚洲成国产人片在线观看| 午夜精品在线福利| 国产成人影院久久av| 国产亚洲av嫩草精品影院| 久久久久九九精品影院| 最新在线观看一区二区三区| 99精品在免费线老司机午夜| 亚洲国产中文字幕在线视频| 校园春色视频在线观看| 99热这里只有精品一区 | 国产高清视频在线播放一区| 精品国产超薄肉色丝袜足j| 欧美三级亚洲精品| 99久久久亚洲精品蜜臀av| 精品国产美女av久久久久小说| 99久久99久久久精品蜜桃| 男人舔女人的私密视频| 免费高清在线观看日韩| 精品欧美一区二区三区在线| 亚洲欧美精品综合一区二区三区| 午夜福利一区二区在线看| 精品国产乱码久久久久久男人| 99久久精品国产亚洲精品| 精品卡一卡二卡四卡免费| 好男人电影高清在线观看| 欧美激情极品国产一区二区三区| 成年人黄色毛片网站| 成人18禁高潮啪啪吃奶动态图| 女同久久另类99精品国产91| 日韩欧美一区二区三区在线观看| 国产av又大| 50天的宝宝边吃奶边哭怎么回事| 999久久久国产精品视频| 国产午夜精品久久久久久| 国产熟女午夜一区二区三区| 精品久久久久久久人妻蜜臀av| 国产高清视频在线播放一区| 精华霜和精华液先用哪个| 天天躁夜夜躁狠狠躁躁| 可以在线观看的亚洲视频| 中国美女看黄片| 国内揄拍国产精品人妻在线 | 他把我摸到了高潮在线观看| 黄色视频不卡| 国产精品九九99| 一个人免费在线观看的高清视频| 国产精品av久久久久免费| 免费观看人在逋| 免费人成视频x8x8入口观看| 午夜免费激情av| av福利片在线| 校园春色视频在线观看| 国产成+人综合+亚洲专区| 亚洲欧美日韩无卡精品| 视频在线观看一区二区三区| 国产又黄又爽又无遮挡在线| 久久欧美精品欧美久久欧美| 18禁美女被吸乳视频| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品久久久久久毛片| 国产精品,欧美在线| 麻豆av在线久日| 高清在线国产一区| 亚洲 国产 在线| 久久久精品欧美日韩精品| 精品无人区乱码1区二区| 精品高清国产在线一区| 午夜影院日韩av| 18禁美女被吸乳视频| 一卡2卡三卡四卡精品乱码亚洲| 精品福利观看| 亚洲av美国av| 在线观看免费午夜福利视频| 成人永久免费在线观看视频| 少妇被粗大的猛进出69影院| 18禁黄网站禁片午夜丰满| 国产熟女xx| 在线视频色国产色| 高清毛片免费观看视频网站| 久久国产亚洲av麻豆专区| 亚洲真实伦在线观看| 日本a在线网址| 一卡2卡三卡四卡精品乱码亚洲| 在线观看66精品国产| 欧美黑人欧美精品刺激| 十分钟在线观看高清视频www| 欧美 亚洲 国产 日韩一| 国产高清激情床上av| 视频在线观看一区二区三区| 在线国产一区二区在线| 啦啦啦免费观看视频1| 免费人成视频x8x8入口观看| 一级毛片女人18水好多| 午夜福利18| 久久久久久久精品吃奶| 最近最新中文字幕大全电影3 | 麻豆av在线久日| 黄色a级毛片大全视频| 日本免费一区二区三区高清不卡| 叶爱在线成人免费视频播放| 最近最新中文字幕大全电影3 | 一本久久中文字幕| 精品卡一卡二卡四卡免费| 在线观看日韩欧美| 欧美精品啪啪一区二区三区| www.999成人在线观看| 国产国语露脸激情在线看| 中文字幕精品免费在线观看视频| 国产主播在线观看一区二区| 在线观看日韩欧美| 亚洲色图av天堂| 欧美激情 高清一区二区三区| 久9热在线精品视频| 亚洲国产欧美网| 97人妻精品一区二区三区麻豆 | 精品卡一卡二卡四卡免费| 成人三级黄色视频| 欧美成人一区二区免费高清观看 | 成人三级黄色视频| 丰满人妻熟妇乱又伦精品不卡| 成人欧美大片| 国产高清视频在线播放一区| 天堂动漫精品| 欧美黑人巨大hd| 欧美日韩亚洲国产一区二区在线观看| 老司机在亚洲福利影院| 国产精品,欧美在线| 给我免费播放毛片高清在线观看| 日日摸夜夜添夜夜添小说| 欧美又色又爽又黄视频| 久久人妻福利社区极品人妻图片| 特大巨黑吊av在线直播 | 十八禁网站免费在线| 18禁黄网站禁片午夜丰满| 久久精品国产综合久久久| 日韩精品青青久久久久久| 99在线人妻在线中文字幕| 欧美激情高清一区二区三区| 国产精华一区二区三区| 亚洲中文av在线| 精品久久久久久久末码| 免费在线观看亚洲国产| 久久久久免费精品人妻一区二区 | avwww免费| 老汉色av国产亚洲站长工具| 欧美性猛交╳xxx乱大交人| 黄频高清免费视频| 男女视频在线观看网站免费 | aaaaa片日本免费| 怎么达到女性高潮| 国产一区在线观看成人免费| 青草久久国产| 又黄又爽又免费观看的视频| 91成年电影在线观看| 成熟少妇高潮喷水视频| 久久久久久免费高清国产稀缺| 国产精品爽爽va在线观看网站 | 人人妻人人澡欧美一区二区| 久久久久国内视频| 成人永久免费在线观看视频| 在线天堂中文资源库| 国产单亲对白刺激| 日韩中文字幕欧美一区二区| 欧美成人性av电影在线观看| xxx96com| 亚洲精品在线美女| 精华霜和精华液先用哪个| 午夜成年电影在线免费观看| 亚洲av电影不卡..在线观看| 中文字幕人成人乱码亚洲影| 视频区欧美日本亚洲| videosex国产| 成人午夜高清在线视频 | 亚洲av电影不卡..在线观看| 亚洲一区二区三区不卡视频| 白带黄色成豆腐渣| 色哟哟哟哟哟哟| 啦啦啦免费观看视频1| 亚洲性夜色夜夜综合| 免费搜索国产男女视频| 色综合欧美亚洲国产小说| 9191精品国产免费久久| 成人特级黄色片久久久久久久| 国产免费av片在线观看野外av| 99在线人妻在线中文字幕| 男人的好看免费观看在线视频 | 激情在线观看视频在线高清| 淫秽高清视频在线观看| 一级毛片女人18水好多| 国产一卡二卡三卡精品| 少妇的丰满在线观看| 在线观看www视频免费| 亚洲性夜色夜夜综合| 成人18禁高潮啪啪吃奶动态图| 村上凉子中文字幕在线| 女警被强在线播放| 少妇 在线观看| 90打野战视频偷拍视频| 亚洲在线自拍视频| 少妇被粗大的猛进出69影院| 在线看三级毛片| 一边摸一边做爽爽视频免费| 亚洲黑人精品在线| 国产精品亚洲美女久久久| 亚洲av美国av| 特大巨黑吊av在线直播 | 国产亚洲欧美精品永久| 一区二区三区高清视频在线| 18美女黄网站色大片免费观看| 国产精品野战在线观看| 男女床上黄色一级片免费看| 国产97色在线日韩免费| 国产99白浆流出| 香蕉丝袜av| 一级毛片女人18水好多| 亚洲第一青青草原| www.精华液| 男人舔奶头视频| 黄片大片在线免费观看| 亚洲午夜精品一区,二区,三区| 欧美激情极品国产一区二区三区| 精品久久久久久成人av| 日韩成人在线观看一区二区三区| 精品一区二区三区四区五区乱码| 老熟妇仑乱视频hdxx| 色精品久久人妻99蜜桃| 成人国产综合亚洲| 久久人妻福利社区极品人妻图片| 色综合婷婷激情| 美女高潮喷水抽搐中文字幕| 欧美国产日韩亚洲一区| 久久狼人影院| 精品欧美一区二区三区在线| 视频在线观看一区二区三区| 精品福利观看| 成人午夜高清在线视频 | 夜夜爽天天搞| 两性午夜刺激爽爽歪歪视频在线观看 | 国产一区二区在线av高清观看| 国产激情久久老熟女| 男人舔奶头视频| 欧美黑人巨大hd| 日日摸夜夜添夜夜添小说| 1024香蕉在线观看| 亚洲av电影不卡..在线观看| 97人妻精品一区二区三区麻豆 | 午夜激情av网站| 极品教师在线免费播放| 天堂影院成人在线观看| 黄色女人牲交| 国产成人欧美在线观看| 丝袜在线中文字幕| tocl精华| 精品国产超薄肉色丝袜足j| 亚洲国产精品成人综合色| 人人妻人人看人人澡| 国产不卡一卡二| 亚洲精华国产精华精| 国产成人欧美在线观看| 在线观看午夜福利视频| 成人18禁在线播放| 男女午夜视频在线观看| av超薄肉色丝袜交足视频| 日日摸夜夜添夜夜添小说| 村上凉子中文字幕在线| 老鸭窝网址在线观看| 久久久久久久久中文| 老司机靠b影院| 一个人免费在线观看的高清视频| 国内久久婷婷六月综合欲色啪| 成人三级做爰电影| 亚洲中文字幕一区二区三区有码在线看 | 禁无遮挡网站| 99国产精品99久久久久| 亚洲精品国产区一区二| 亚洲精品中文字幕一二三四区| 久久人妻av系列| 国产成人欧美| 香蕉丝袜av| 在线播放国产精品三级| 久久久久久人人人人人| 国产黄色小视频在线观看| 成人18禁在线播放| 亚洲精品国产区一区二| 亚洲精品美女久久av网站| 精品无人区乱码1区二区| 激情在线观看视频在线高清| 悠悠久久av| 女生性感内裤真人,穿戴方法视频| 欧美色欧美亚洲另类二区| 白带黄色成豆腐渣| 老汉色av国产亚洲站长工具| 国产亚洲欧美98| 久久国产精品影院| 亚洲精品粉嫩美女一区| 国产单亲对白刺激| 亚洲国产高清在线一区二区三 | 桃色一区二区三区在线观看| 婷婷精品国产亚洲av| 国产精品一区二区三区四区久久 | www日本在线高清视频| 亚洲第一青青草原| 日韩一卡2卡3卡4卡2021年| 精品福利观看| 亚洲精品粉嫩美女一区| 一卡2卡三卡四卡精品乱码亚洲| 欧美一区二区精品小视频在线| 丝袜美腿诱惑在线| 日日夜夜操网爽| 男女下面进入的视频免费午夜 | 亚洲国产日韩欧美精品在线观看 | 一边摸一边抽搐一进一小说| 真人做人爱边吃奶动态| 午夜福利欧美成人| 亚洲一卡2卡3卡4卡5卡精品中文| av电影中文网址| 宅男免费午夜| 给我免费播放毛片高清在线观看| 久久久久久久久中文| 在线视频色国产色| 嫁个100分男人电影在线观看| 久久精品国产综合久久久| 精品国产美女av久久久久小说| 中亚洲国语对白在线视频| 制服人妻中文乱码| 欧美日本亚洲视频在线播放| 欧美日韩亚洲综合一区二区三区_| 欧洲精品卡2卡3卡4卡5卡区| 老汉色∧v一级毛片| 免费观看人在逋| 淫妇啪啪啪对白视频| 色在线成人网| 色精品久久人妻99蜜桃| 国产1区2区3区精品| 男女下面进入的视频免费午夜 | 真人做人爱边吃奶动态| 久久精品国产亚洲av香蕉五月| 国产av又大| 欧美乱色亚洲激情| 亚洲,欧美精品.| 桃色一区二区三区在线观看| 国产精品亚洲一级av第二区| 久久久久久国产a免费观看| 精品日产1卡2卡| 亚洲国产欧美日韩在线播放| 久久精品国产亚洲av高清一级| 精品高清国产在线一区| 1024香蕉在线观看| 91老司机精品| 女性生殖器流出的白浆| 国产亚洲精品一区二区www| 免费电影在线观看免费观看| 18禁黄网站禁片午夜丰满| 国产精品电影一区二区三区| 1024香蕉在线观看| 午夜福利成人在线免费观看| 妹子高潮喷水视频| 亚洲性夜色夜夜综合| 亚洲在线自拍视频| 国产真人三级小视频在线观看| 1024手机看黄色片| 黄色成人免费大全| 中文字幕另类日韩欧美亚洲嫩草| 黄色丝袜av网址大全| 此物有八面人人有两片| 俄罗斯特黄特色一大片| 在线播放国产精品三级| 男人舔女人的私密视频| 日本 欧美在线| 国产亚洲精品av在线| 黄色视频不卡| 一级毛片女人18水好多| 亚洲av美国av| 观看免费一级毛片| 丰满的人妻完整版| www.精华液| 1024香蕉在线观看| 99久久精品国产亚洲精品| 欧美三级亚洲精品| 黄频高清免费视频| 可以在线观看毛片的网站| 久久国产精品男人的天堂亚洲| 好男人在线观看高清免费视频 | 男女之事视频高清在线观看| 亚洲欧美日韩高清在线视频| 免费av毛片视频| 欧美黄色淫秽网站| 91字幕亚洲| 好男人在线观看高清免费视频 | 日韩视频一区二区在线观看| 一本大道久久a久久精品| 日日干狠狠操夜夜爽| 欧美不卡视频在线免费观看 | 最近最新中文字幕大全电影3 | √禁漫天堂资源中文www| 精品国产美女av久久久久小说| 亚洲人成网站高清观看| 久久久久久久午夜电影| 两性午夜刺激爽爽歪歪视频在线观看 | а√天堂www在线а√下载| 欧美不卡视频在线免费观看 | 国产片内射在线| 亚洲成人国产一区在线观看| 久久 成人 亚洲| 国产精品av久久久久免费| 久久性视频一级片| 亚洲精品中文字幕在线视频| 欧美日韩福利视频一区二区| 在线视频色国产色| 国产高清videossex| www.999成人在线观看| 男男h啪啪无遮挡| 国产亚洲精品久久久久5区| 岛国在线观看网站| 男人的好看免费观看在线视频 | 国产三级黄色录像| 超碰成人久久| a在线观看视频网站| 欧美亚洲日本最大视频资源| 97超级碰碰碰精品色视频在线观看| ponron亚洲| 禁无遮挡网站| 在线观看66精品国产| 欧美日韩乱码在线| 久久 成人 亚洲| 午夜福利一区二区在线看| 人妻久久中文字幕网| 亚洲精品一区av在线观看| 久久久国产成人免费| 国产亚洲av高清不卡| 国产高清有码在线观看视频 | 法律面前人人平等表现在哪些方面| 精品人妻1区二区| 午夜激情福利司机影院| 色哟哟哟哟哟哟| 国产精华一区二区三区| 日本成人三级电影网站| 亚洲黑人精品在线| 在线十欧美十亚洲十日本专区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av成人一区二区三| 国产三级黄色录像| 啦啦啦韩国在线观看视频| 亚洲第一青青草原| 成人18禁高潮啪啪吃奶动态图| 日韩欧美一区视频在线观看| 久久人妻av系列| 亚洲av日韩精品久久久久久密| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品久久国产高清桃花| 好男人电影高清在线观看| 可以免费在线观看a视频的电影网站| 国产免费av片在线观看野外av| 精品国内亚洲2022精品成人| 国产伦人伦偷精品视频| 亚洲三区欧美一区| 精品一区二区三区四区五区乱码| 欧美黑人精品巨大| 可以在线观看毛片的网站| 国产伦人伦偷精品视频| 亚洲国产高清在线一区二区三 | 身体一侧抽搐| 女人被狂操c到高潮| 国产欧美日韩一区二区三| 免费一级毛片在线播放高清视频| 狠狠狠狠99中文字幕| 一边摸一边做爽爽视频免费| 亚洲人成伊人成综合网2020| 草草在线视频免费看| 啦啦啦观看免费观看视频高清| 国内毛片毛片毛片毛片毛片| 国产亚洲av高清不卡| 国产成人一区二区三区免费视频网站| 日本 av在线| 久久青草综合色| 国产黄色小视频在线观看| 亚洲熟女毛片儿| 免费在线观看成人毛片| 国产精品久久久人人做人人爽| 久久 成人 亚洲| 国产成人精品久久二区二区91| 国产一区二区三区视频了| 欧美丝袜亚洲另类 | 亚洲精品一区av在线观看| 成人18禁在线播放| 欧美性猛交黑人性爽| 亚洲色图av天堂| 日韩有码中文字幕| 亚洲精品久久成人aⅴ小说| 久久青草综合色| 日本a在线网址| 欧美三级亚洲精品| 性欧美人与动物交配| 黄片小视频在线播放| 色播在线永久视频| 亚洲精品久久国产高清桃花| 亚洲精品在线观看二区| 国产高清videossex| 国产精品爽爽va在线观看网站 | 一夜夜www| 国产一区二区三区视频了| 久久精品国产99精品国产亚洲性色| 欧美黄色片欧美黄色片| 国产免费男女视频| tocl精华| 国产又色又爽无遮挡免费看| 可以在线观看毛片的网站| 欧美zozozo另类| 一区二区三区高清视频在线| 精品久久久久久成人av| av有码第一页| 一级a爱视频在线免费观看| 国产成人一区二区三区免费视频网站| 黄色女人牲交| 久久精品国产综合久久久| 丝袜人妻中文字幕| 久久人人精品亚洲av| 日本一本二区三区精品| 免费女性裸体啪啪无遮挡网站| 不卡一级毛片| 午夜福利成人在线免费观看| 国产一区二区激情短视频| 美女国产高潮福利片在线看| 精品欧美国产一区二区三| 97人妻精品一区二区三区麻豆 | 在线永久观看黄色视频| 国产亚洲精品av在线| 嫁个100分男人电影在线观看| 久久久国产欧美日韩av| 高潮久久久久久久久久久不卡| 黄频高清免费视频| 99久久久亚洲精品蜜臀av| 美女扒开内裤让男人捅视频| 深夜精品福利| 麻豆一二三区av精品| 久久中文看片网| 欧美zozozo另类| 日韩国内少妇激情av| 天天添夜夜摸| 国产精品亚洲美女久久久| 亚洲av第一区精品v没综合| 国产乱人伦免费视频| 一本久久中文字幕| 一区二区三区高清视频在线| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品sss在线观看| 日韩欧美在线二视频| 美女高潮喷水抽搐中文字幕| 日日摸夜夜添夜夜添小说| 制服丝袜大香蕉在线| 免费无遮挡裸体视频| 给我免费播放毛片高清在线观看| 成人国产一区最新在线观看| 久久久久久亚洲精品国产蜜桃av| 欧美黑人精品巨大| 黄色a级毛片大全视频| 欧美日韩中文字幕国产精品一区二区三区| 国产午夜精品久久久久久| 看黄色毛片网站| 两个人视频免费观看高清| 欧美不卡视频在线免费观看 | 视频在线观看一区二区三区| 好男人电影高清在线观看| 亚洲av电影不卡..在线观看| 久久中文看片网| 成人精品一区二区免费| 51午夜福利影视在线观看| 精品国产亚洲在线| 欧美另类亚洲清纯唯美| 嫩草影院精品99| 身体一侧抽搐| 日本在线视频免费播放| 中文字幕另类日韩欧美亚洲嫩草| ponron亚洲| 国产又黄又爽又无遮挡在线| 满18在线观看网站| 午夜福利在线在线| 人妻丰满熟妇av一区二区三区| av超薄肉色丝袜交足视频| 村上凉子中文字幕在线| 丁香六月欧美| 十八禁人妻一区二区| 搡老岳熟女国产| 看片在线看免费视频| 青草久久国产|