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

    基于DS-InSAR的西部礦區(qū)地表時序沉降監(jiān)測與分析

    2022-05-30 01:28:30譚志祥鄧喀中
    金屬礦山 2022年5期
    關(guān)鍵詞:時序盆地礦區(qū)

    王 波 譚志祥,2 鄧喀中,2

    (1.中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院,江蘇 徐州 221116;2.自然資源部國土環(huán)境與災(zāi)害監(jiān)測重點實驗室,江蘇 徐州 221116)

    多年來對東部礦區(qū)煤炭資源的高強(qiáng)度開發(fā),使得 東部礦區(qū)煤炭產(chǎn)量不斷下降。為滿足國民經(jīng)濟(jì)發(fā)展對于煤炭資源的需求,我國加大了對西部礦區(qū)煤炭資源的開發(fā)利用力度。與東部礦區(qū)相比,西部礦區(qū)煤炭資源儲量更大,而且具有煤層厚度大、傾角小等優(yōu)點,因此有著更大的生產(chǎn)潛力[1]。然而由于歷史原因,西部礦區(qū)煤礦開采時間比東部礦區(qū)要短,因此針對西部礦區(qū)開采引起的地表形變的監(jiān)測工作開展較少,相應(yīng)的研究成果也比較匱乏。而且由于東西部礦區(qū)之間地質(zhì)采礦條件的巨大差異,導(dǎo)致以東部礦區(qū)為研究對象所得到的相關(guān)研究成果在西部礦區(qū)適用性不強(qiáng)[2],這給西部礦區(qū)地表形變預(yù)計、村莊保護(hù)煤柱留設(shè)以及各類建筑物保護(hù)帶來了困難,不利于煤礦生產(chǎn)工作安全高效進(jìn)行。因此有必要針對西部礦區(qū)特定的地質(zhì)采礦條件,開展礦區(qū)地表形變監(jiān)測工作,并以此為基礎(chǔ)進(jìn)行西部礦區(qū)地表形變規(guī)律和機(jī)理的研究,為后續(xù)煤礦生產(chǎn)工作提供技術(shù)依據(jù)。

    水準(zhǔn)測量等傳統(tǒng)的地表形變監(jiān)測方法在西部礦區(qū)的應(yīng)用具有較大的局限性,主要體現(xiàn)在受西部礦區(qū)高原堆積型沙丘地貌的影響,地面水準(zhǔn)點以及地表特征點的選取和觀測困難。而且傳統(tǒng)的測量方法往往費時費力,難以進(jìn)行大范圍面狀地表形變監(jiān)測[3]。合成孔徑雷達(dá)差分干涉測量技術(shù)(Differential InSAR,D-InSAR)能夠不依賴地面特征點來獲取地表形變信息,因此被廣泛應(yīng)用于地表形變的監(jiān)測與分析中[4]。然而該技術(shù)極易受到時空去相干因素的影響,難以獲取礦區(qū)開采引起的地表形變范圍。針對這一問題,科研工作者將注意力集中在了那些能夠在長時間序列合成孔徑雷達(dá)(Synthetic Aperture Radar,SAR)影像中保持高相干性的點目標(biāo)上,提出了永久散射體干涉測量技術(shù)(Permanent Scatterer InSAR,PS-InSAR)[5],在此基礎(chǔ)上進(jìn)行了相干點目標(biāo)識別[6]、三維相位解纏算法[7]等研究,將D-InSAR技術(shù)推進(jìn)到了時間序列InSAR的新階段。以PS-InSAR和小基線集方法(Small Baseline Subsets InSAR,SBAS-InSAR)[8]為代表的常規(guī)時序InSAR技術(shù)在干涉對優(yōu)化的基礎(chǔ)上,提取分析了影像中受時空去相干影響較小而表現(xiàn)出相對穩(wěn)定散射特性的高相干像元在整個監(jiān)測時段內(nèi)的變化,有效克服了D-InSAR技術(shù)在應(yīng)用過程中的失相干問題。然而礦區(qū)及周邊地表的覆蓋多為植被和沙地,導(dǎo)致常規(guī)的時序InSAR技術(shù)在礦區(qū)及其周圍難以獲得足夠數(shù)量、分布均勻的高相干點,地表形變監(jiān)測精度受到了很大限制。

    近些年發(fā)展起來的融合分布式目標(biāo)的時序InSAR技術(shù)[9]將研究的重點從散射特性相對穩(wěn)定的高相干點目標(biāo)擴(kuò)大到了散射特性中等的分布式目標(biāo),可同時利用這兩種目標(biāo)進(jìn)行地表形變監(jiān)測,從而極大地提高了監(jiān)測區(qū)域點目標(biāo)的密度,有效克服了常規(guī)時序InSAR技術(shù)在低相干區(qū)域由于無法獲取足夠數(shù)量的點目標(biāo)而無法對地表形變進(jìn)行準(zhǔn)確監(jiān)測的不足。該技術(shù)先是在地震形變監(jiān)測、山區(qū)滑坡災(zāi)害監(jiān)測、冰川運動監(jiān)測等方面進(jìn)行應(yīng)用并取得了良好的效果[10],證明了該方法與常規(guī)時序InSAR技術(shù)相比具有優(yōu)越性。之后部分學(xué)者開始嘗試將該技術(shù)應(yīng)用到礦區(qū)地表形變監(jiān)測中并取得了一定的成果[11-18]。但是上述研究的試驗區(qū)域大多位于東部礦區(qū),對于該技術(shù)在西部礦區(qū)地表形變監(jiān)測方面的應(yīng)用研究涉及較少,缺乏針對性的研究。為此,本研究在現(xiàn)有成果的基礎(chǔ)上,利用一種針對西部礦區(qū)高原堆積型沙丘地貌的基于DS-InSAR的時序地表形變監(jiān)測方法,獲取了西部礦區(qū)內(nèi)蒙古石拉烏素煤礦地表形變及其時空分布信息,并驗證了該方法的可靠性,豐富了DS-InSAR技術(shù)在西部礦區(qū)地表形變監(jiān)測應(yīng)用方面的研究。

    1 研究區(qū)概況與數(shù)據(jù)源

    1.1 研究區(qū)概況

    石拉烏素煤礦位于內(nèi)蒙古自治區(qū)鄂爾多斯市伊金霍洛旗東勝煤田內(nèi),礦區(qū)面積約70 km2,設(shè)計生產(chǎn)能力為1 000萬t/a。礦區(qū)位于鄂爾多斯高原的中南部,具有典型的高原堆積型沙丘地貌特征,地表大部分被第四系風(fēng)積沙所覆蓋,局部有白堊下統(tǒng)地層露出,植被稀疏,為沙漠—半沙漠地區(qū)。本文研究區(qū)域所在的221采區(qū)位于石拉烏素煤礦的中部(圖1),采區(qū)地形總體趨勢是東高西低、北高南低,地面標(biāo)高+1 332~+1 411 m,采區(qū)內(nèi)煤層埋深近700 m,上覆大厚度的中生代砂巖,上覆巖層性質(zhì)具體如表1所示,其地質(zhì)采礦條件與東部礦區(qū)相比存在明顯差異。在監(jiān)測時段內(nèi)221采區(qū)有17、18、06A 3個工作面開采,采區(qū)北部的17、18工作面采用綜合機(jī)械化一次采全高工藝開采,平均采厚約5.40 m,南部的06A工作面采用綜合機(jī)械化放頂煤工藝開采,平均采厚約9.78 m。各工作面的參數(shù)取值見表2。

    表1 煤層上覆巖層性質(zhì)Table 1 Properties of overlying strata of coal seam

    表2 221采區(qū)工作面參數(shù)Table 2 Parameters of working face in 221mining area

    圖1 研究區(qū)地理位置Fig.1 Geographical locationof the study area

    1.2 數(shù)據(jù)源

    Sentinel-1衛(wèi)星由歐空局發(fā)射,并免費向全球提供衛(wèi)星影像下載服務(wù),為SAR相關(guān)研究提供了極大的便利。本研究選用了2016年9月—2018年6月共52景覆蓋研究區(qū)域的IW模式下VV極化的Sentinel-1A升軌影像。Sentinel-1衛(wèi)星影像的分辨率為5 m×20 m,衛(wèi)星的重訪周期為12 d。此外本研究采用30m分辨率的航天飛機(jī)雷達(dá)測圖計劃(ShuttleRadar Topography Mission,SRTM)數(shù)字高程模型(Digital Elevation Model,DEM)從干涉圖中去除地形相位的影響和地理編碼。

    2 監(jiān)測方法選擇及其原理

    2.1 高相干點目標(biāo)選取

    從時序SAR影像中準(zhǔn)確地識別高相干點目標(biāo)是時間序列InSAR分析的重要環(huán)節(jié),也是實現(xiàn)聯(lián)合高相干點目標(biāo)與分布式目標(biāo)精確提取地表形變信息的前提條件。受西部礦區(qū)高原堆積型沙丘地貌的影響,研究區(qū)域內(nèi)能夠在長時間序列中保持相位穩(wěn)定的永久散射體目標(biāo)的數(shù)量稀少,但存在數(shù)量較多在短時間內(nèi)經(jīng)濾波后相位失相關(guān)性變小而保持高相干性的點目標(biāo)。本研究對所選SAR影像進(jìn)行處理得到濾波后的差分干涉圖集,通過幅度差分離差法[19]選取出了此類點目標(biāo)作為高相干點目標(biāo)。幅度差分離差法原理可用下式表示

    式中,DΔA為幅度差分離差值;σΔA為主、輔影像幅度差的標(biāo)準(zhǔn)差;μA為差分干涉圖幅度的均值。通過設(shè)定幅度差分離差閾值,便可實現(xiàn)高相干點目標(biāo)的選取。

    2.2 分布式目標(biāo)識別與相位優(yōu)化

    2.2.1 FaSHPS同質(zhì)像元識別

    受礦區(qū)特定地質(zhì)條件的影響,研究區(qū)域內(nèi)高相干點目標(biāo)的空間分布極不均勻,大部分高相干點目標(biāo)集中分布在地表裸露巖石較多的區(qū)域,而在地表大范圍的積沙區(qū)域內(nèi)此類點目標(biāo)數(shù)量十分稀少。為實現(xiàn)對研究區(qū)完整地表形變信息的提取,必須在高相干點目標(biāo)的基礎(chǔ)上充分利用廣大積沙區(qū)域內(nèi)的分布式目標(biāo)。參數(shù)假設(shè)檢驗方法中的快速同質(zhì)點選取方法(Fast Statistically Homogeneous Pixel Selection,FaSHPS)[20]具有同質(zhì)像元選取效率高且提取出的同質(zhì)像元之間的異質(zhì)性小的優(yōu)點,為此本研究選用該方法對研究區(qū)域內(nèi)的同質(zhì)像元進(jìn)行識別,并通過設(shè)置時間相干性閾值實現(xiàn)對分布式目標(biāo)的提取。FaSHPS方法的思路為,假設(shè)時序SAR影像數(shù)據(jù)集內(nèi)任一像元時間維上振幅的均值服從高斯分布,從而可以將傳統(tǒng)的假設(shè)檢驗方法轉(zhuǎn)化為置信區(qū)間的估計問題,只要利用簡單的邏輯計算來判斷兩個像元服從的分布函數(shù)是否相同,便可實現(xiàn)同質(zhì)像元的識別。

    設(shè)由N景時序SAR影像組成數(shù)據(jù)集,其中任一像元L在時間維上的振幅均值為(L),由中心極限定律可知,隨著樣本數(shù)目N加大,振幅均值(L)趨于高斯分布,此時其區(qū)間估計可表示為

    式中,P{·}表示概率;μ(L)為像元L的數(shù)學(xué)期望;z1-α/2為標(biāo)準(zhǔn)正態(tài)分布中置信度為α?xí)r對應(yīng)的分位點;Var(A(L))為像元L時間維上振幅的真實方差。

    當(dāng)樣本的數(shù)量足夠大時,單視復(fù)數(shù)SAR影像的時序振幅信息近似服從瑞利分布,此時變異系數(shù)CV為常量,其值約為0.52。若像元在時間維上的后向散射特性穩(wěn)定,則式(2)可表示為

    當(dāng)給定置信度α?xí)r,式(3)便給定了一個區(qū)間,若待估像元在時間維上的平均強(qiáng)度值在目標(biāo)像元所對應(yīng)的置信區(qū)間內(nèi),則待估像元與目標(biāo)像元為同質(zhì)像元。

    2.2.2 特征值分解相位優(yōu)化

    分布式目標(biāo)像元內(nèi)地物勻質(zhì)分布[21],其回波信號由像元內(nèi)各子散射體相干疊加而成,易受時空失相干影響,相位穩(wěn)定性較差。因此需要基于分布式目標(biāo)像元散射特性的整體統(tǒng)計特征來去除噪聲相位,實現(xiàn)對相位的優(yōu)化。特征值分解法在相位優(yōu)化的過程中無需迭代計算,因此其優(yōu)化效率較高,同時該方法的相位優(yōu)化質(zhì)量也較好。故本研究選用該方法對分布式目標(biāo)像元相位進(jìn)行優(yōu)化。特征值分解法的核心思想為從相干矩陣中分離,估計出多元散射機(jī)制并得到其對應(yīng)的相位分量,通過分離出主導(dǎo)散射機(jī)制的相位分量來達(dá)到優(yōu)化相位的目的[22-23]。

    式中,Λ為N階對角矩陣,對角線上的元素為相干矩陣的N個特征值,用λi表示;每個不同的特征值分別對應(yīng)著不同的特征向量,用μi表示;特征值不同,其對應(yīng)的散射機(jī)制也不同。特征值越大,在像元整體散射信號中其對應(yīng)的子散射體所占比重就越大。因此將最大特征值對應(yīng)的特征向量作為主散射體,其余的特征向量作為噪聲,就可以實現(xiàn)同質(zhì)像元相位的優(yōu)化,此時相干矩陣可以表示為

    在完成同質(zhì)像元識別以及相位優(yōu)化之后,通過計算相位優(yōu)化前后干涉相位之間的擬合度便可對相位優(yōu)化的質(zhì)量進(jìn)行評價。相位優(yōu)化質(zhì)量評估完成之后,設(shè)置時間相干性閾值,從分布式目標(biāo)中篩選出相干性較好、信噪比較高的部分。將其與篩選出的高相干點目標(biāo)一并融入時序InSAR處理流程中,便可對整個研究區(qū)域地表形變進(jìn)行高精度監(jiān)測。針對西部礦區(qū)高原堆積型沙丘地貌的基于DS-InSAR的時序地表形變監(jiān)測方法的主要流程如2所示。

    3 監(jiān)測結(jié)果與可靠性分析

    3.1 SBAS-InSAR與DS-InSAR監(jiān)測結(jié)果對比

    本研究分別利用傳統(tǒng)的SBAS-InSAR方法和圖2所示的基于DS-InSAR的時序地表形變監(jiān)測方法提取了2016年9月—2018年6月研究區(qū)地表沉降的時空分布信息。兩種方法獲取的研究區(qū)地表衛(wèi)星視線向的年平均沉降速率如圖3所示。傳統(tǒng)的SBASInSAR方法共提取到了57 803個監(jiān)測點,監(jiān)測點數(shù)目較少,且空間分布極不均勻。SBAS-InSAR方法監(jiān)測到研究區(qū)地表存在一處較為明顯的下沉盆地,最大年平均沉降速率出現(xiàn)在該下沉盆地的中心位置,為139 mm/a,在該下沉盆地的南部同樣監(jiān)測到一處較為明顯的沉降區(qū)域,但由于在該沉降區(qū)域內(nèi)獲取的監(jiān)測點數(shù)量稀少,只能大致反映該區(qū)域的外圍輪廓,無法得到其整體的沉降信息。DS-InSAR方法共提取到了239 750個監(jiān)測點,與傳統(tǒng)的SBAS-InSAR方法相比該方法提取到的監(jiān)測點數(shù)量顯著增加,約為前者的4倍,且監(jiān)測點空間分布均勻。DS-InSAR方法監(jiān)測到研究區(qū)地表存在南北方向相鄰的2個明顯的下沉盆地,且在2個下沉盆地中北部下沉盆地的沉降更為明顯,沉降范圍也更大。在北部下沉盆地的中心位置監(jiān)測到了研究區(qū)域內(nèi)的最大年平均沉降速率,其值為179 mm/a。

    圖2 基于DS-InSAR的時序地表形變監(jiān)測方法數(shù)據(jù)處理流程Fig.2 Data processing flow of time series surface deformation monitoring method based on DS-InSAR

    圖3 研究區(qū)地表衛(wèi)星視線向年平均沉降速率(單位:mm/a)Fig.3 Annual average subsidence rate of LOS direction in study area

    結(jié)合衛(wèi)星影像底圖對比分析兩種方法的監(jiān)測結(jié)果可知,受西部礦區(qū)高原堆積型沙丘地貌的影響,傳統(tǒng)的SBAS-InSAR方法只能在礦區(qū)地表裸露巖石較多的區(qū)域獲取足夠數(shù)量的高相干點目標(biāo),在積沙區(qū)域由于地物的后向散射特性較弱獲取的觀測點數(shù)量稀少,監(jiān)測結(jié)果難以準(zhǔn)確反映礦區(qū)地表整體的沉降信息,不能對沉降區(qū)域范圍進(jìn)行準(zhǔn)確探測。本研究基于DS-InSAR的時序地表形變監(jiān)測方法有效克服了傳統(tǒng)SBAS-InSAR方法在西部礦區(qū)地表沉降監(jiān)測中的不足,在充分利用裸露巖石較多區(qū)域的高相干點目標(biāo)的基礎(chǔ)上,通過提取研究區(qū)域內(nèi)的分布式目標(biāo)顯著增加了積沙區(qū)域觀測點的數(shù)量,從而能夠?qū)ΦV區(qū)地表形變進(jìn)行更高精度地監(jiān)測,獲取礦區(qū)地表整體豐富的沉降細(xì)節(jié)信息。

    3.2 DS-InSAR時序監(jiān)測結(jié)果

    以SAR影像數(shù)據(jù)集中首期影像的日期為開始時間,計算其余影像對應(yīng)時間相對于開始時間的累計沉降值,將觀測結(jié)果從衛(wèi)星的視線方向轉(zhuǎn)換為垂直方向,即可得到礦區(qū)地表的時序累積沉降。為更加直觀地顯示整個監(jiān)測時段內(nèi)研究區(qū)地表沉降區(qū)域的形成與地下工作面開采之間的關(guān)系,將地下工作面的位置疊加到時序累積沉降圖上,并提取出研究區(qū)域內(nèi)受采動影響地表發(fā)生沉降的區(qū)域,即可得到沉降區(qū)域的累計沉降及沉降細(xì)節(jié)信息。研究區(qū)地表沉降區(qū)域范圍如圖4所示,沉降區(qū)域內(nèi)的累計沉降及沉降細(xì)節(jié)信息如圖5所示。

    圖4 研究區(qū)地表沉降區(qū)域范圍Fig.4 Scope of ground subsidence in the study area

    圖5 研究區(qū)地表累計沉降及沉降細(xì)節(jié)信息Fig.5 Surface cumulative subsidence and sunsidence details in the study area

    由圖5可知:地表發(fā)生沉降的區(qū)域與地下工作面的位置一致。17、18兩個工作面的開采導(dǎo)致研究區(qū)地表北部下沉盆地形成,本次試驗監(jiān)測到的最大累計沉降出現(xiàn)在該下沉盆地的中心位置,為402 mm。受06A工作面開采的影響,研究區(qū)地表形成了南部下沉盆地,該盆地內(nèi)監(jiān)測到的累計沉降較小,其中心位置的最大累計沉降約220 mm。監(jiān)測結(jié)果表明,針對西部礦區(qū)高原堆積型沙丘地貌的基于DS-InSAR的時序地表形變監(jiān)測方法在研究區(qū)域內(nèi)獲取了更多數(shù)量、空間分布均勻的監(jiān)測點目標(biāo),能夠?qū)ΦV區(qū)地表開采引起的沉降區(qū)域的形成和發(fā)展過程進(jìn)行直觀地描述,并準(zhǔn)確地顯示地表沉降區(qū)域的范圍。

    3.3 DS-InSAR監(jiān)測結(jié)果可靠性分析

    本研究選取了17、18工作面上方2條觀測線在4個不同時段內(nèi)的水準(zhǔn)實測數(shù)據(jù),通過將觀測線上DS-InSAR監(jiān)測結(jié)果與水準(zhǔn)實測數(shù)據(jù)進(jìn)行對比分析來驗證其可靠性。觀測線與工作面的位置關(guān)系如圖6所示。因水準(zhǔn)實測數(shù)據(jù)的日期與SAR影像日期并非完全一致,故將DS-InSAR監(jiān)測結(jié)果按時間進(jìn)行插值使其與水準(zhǔn)實測數(shù)據(jù)的日期保持一致。本研究在2條觀測線上共選取了93個監(jiān)測點目標(biāo),觀測線上DS-InSAR沉降監(jiān)測結(jié)果與水準(zhǔn)實測數(shù)據(jù)對比如圖7所示。

    圖6 觀測線與工作面位置關(guān)系Fig.6 Position relationship between observation lines and working faces

    由圖7可知:DS-InSAR監(jiān)測結(jié)果與水準(zhǔn)實測數(shù)據(jù)擬合較好,除極個別點該方法監(jiān)測結(jié)果與水準(zhǔn)實測數(shù)據(jù)偏差略大外,各觀測線上二者得到的累計沉降值及沿觀測線的沉降趨勢基本一致。為更加直觀地體現(xiàn)兩種方法監(jiān)測結(jié)果間的相關(guān)性與誤差分布,繪制了如圖8所示的相關(guān)性及誤差分布圖,并通過計算每條觀測線上的誤差均值、絕對誤差均值及標(biāo)準(zhǔn)差,對監(jiān)測結(jié)果的精度進(jìn)行評價,結(jié)果見表3。

    圖7 DS-InSAR監(jiān)測結(jié)果與水準(zhǔn)實測數(shù)據(jù)對比Fig.7 Comparison of DS-InSAR monitoring results and leveling measured data

    圖8 DS-InSAR監(jiān)測結(jié)果與水準(zhǔn)實測數(shù)據(jù)相關(guān)性及誤差分布Fig.8 Correlation and error distribution between DS-InSAR monitoring results and leveling measured data

    由圖8和表3可知:兩種方法監(jiān)測結(jié)果間的相關(guān)性較高,皮爾遜相關(guān)系數(shù)為0.97,二者間的誤差較小且整體分布合理,證明針對西部礦區(qū)高原堆積型沙丘地貌的基于DS-InSAR的時序地表沉降監(jiān)測方法在研究區(qū)地表非充分采動情況下的監(jiān)測結(jié)果具有較高的可靠性。

    表3 DS-InSAR監(jiān)測結(jié)果精度評價Table 3 Accuracy evaluation of DS-InSAR monitoring results

    4 沉降監(jiān)測結(jié)果分析

    4.1 地表沉降與地下開采關(guān)聯(lián)性分析

    結(jié)合采工圖和DS-InSAR時序監(jiān)測結(jié)果可知,2016年9月17工作面開始掘進(jìn),但此后較長一段時間內(nèi)地表并未監(jiān)測到沉降,這是由于工作面的采深大且上方覆有大厚度的中生代砂巖,工作面開采影響尚未傳播到地表。地表開始發(fā)生沉降時,工作面推進(jìn)距離已達(dá)采動距,之后隨著工作面繼續(xù)推進(jìn),地表沉降區(qū)域的范圍和累計沉降值不斷增大。DS-InSAR監(jiān)測結(jié)果顯示,2016年11月份左右研究區(qū)地表開始出現(xiàn)沉降,12月份地表已經(jīng)形成一處較為明顯的沉降區(qū)域,2017年5月份17工作面開采結(jié)束后,沉降區(qū)域發(fā)展成為一個明顯的下沉盆地,其形狀近似為圓形,盆地內(nèi)各點的下沉值隨著與盆地中心距離的增加而減小。2017年6月,17工作面西鄰的18工作面開始回采,由于兩個工作面之間相距很近,18工作面的開采并沒有在礦區(qū)地表形成新的下沉盆地,而是使原有下沉盆地的范圍和累計沉降值進(jìn)一步增大,但受兩個工作面之間留設(shè)的120 m寬保護(hù)煤柱以及上覆巖層應(yīng)力狀態(tài)的影響,地表仍未達(dá)到充分采動且下沉極不充分。2017年12月份18工作面回采結(jié)束,其平均采厚達(dá)5.47 m,但直至2018年6月監(jiān)測結(jié)束該下沉盆地內(nèi)監(jiān)測到的最大沉降也僅為402 mm。

    06A工作面從2017年9月開始掘進(jìn),因此根據(jù)2017年12月地表沉降監(jiān)測結(jié)果可知,除了在開采已經(jīng)結(jié)束的17、18工作面上方存在一處下沉盆地外,在其南部也出現(xiàn)了一處較為明顯的沉降區(qū)域。2018年4月06A工作面推進(jìn)結(jié)束,到2018年6月該工作面上方已經(jīng)形成了一處明顯的下沉盆地,且其沉降范圍已經(jīng)與北部下沉盆地相當(dāng)。通過比較南北兩個下沉盆地的形成和發(fā)展過程不難發(fā)現(xiàn),從地表開始發(fā)生沉降到形成明顯的下沉盆地過程中南部下沉盆地所用時間要明顯短于北部下沉盆地。這是因為雖然3個工作面的采深和上覆巖層的性質(zhì)十分接近,但06A工作面采用綜合機(jī)械化放頂煤工藝開采,其采厚達(dá)9.78 m,遠(yuǎn)大于17、18兩個工作面,因此其開采影響傳播到地表所用的時間更短,但受該礦區(qū)特殊地質(zhì)采礦條件的影響,直至監(jiān)測時段結(jié)束該工作面上方地表同樣未達(dá)到充分采動。

    4.2 南部沉降區(qū)地表動態(tài)變形分析

    為進(jìn)一步驗證基于DS-InSAR的時序地表沉降監(jiān)測方法的性能,本研究利用該方法獲取了圖9所示06A工作面開采期間南部沉降區(qū)地表沿該工作面走向和傾向主斷面的動態(tài)沉降曲線,對該區(qū)域缺少的水準(zhǔn)實測數(shù)據(jù)起到了很好的補充作用。為充分反映工作面推進(jìn)過程中地表的動態(tài)變形情況,將沉降監(jiān)測結(jié)果進(jìn)一步處理得到了沿06A工作面走向主斷面的動態(tài)傾斜與動態(tài)曲率曲線,分別如圖10和圖11所示。由圖10、圖11可知:在06A工作面推進(jìn)過程中,沿該工作面走向的傾斜和曲率變形不斷增大,至2018年4月底開始趨于穩(wěn)定,傾斜在20號點附近(采空區(qū)中心)近似為0,在8號點和26號點附近(采空區(qū)邊界)達(dá)到最大,為0.41 mm/m;地表曲率變形較小,均小于0.02 mm/m2。

    圖9 06A工作面走向與傾向主斷面動態(tài)沉降曲線Fig.9 Dynamic subsidence curves of main section of strike and incline of 06A working face

    圖10 06A工作面走向主斷面動態(tài)傾斜曲線Fig.10 Dynamic inclination curves of main section of strike of 06A working face

    圖11 06A工作面走向主斷面動態(tài)曲率曲線Fig.11 Dynamic curvature curves of main section of strike of 06A working face

    5 結(jié) 論

    本研究針對西部礦區(qū)高原堆積型沙丘地貌,以內(nèi)蒙古石拉烏素煤礦為例,采用基于DS-InSAR的時序地表形變監(jiān)測方法對2016—2018年間覆蓋研究區(qū)域的52景Sentinel-1A影像數(shù)據(jù)進(jìn)行處理,獲取了研究區(qū)域的地表時序沉降信息,并將監(jiān)測結(jié)果結(jié)合井下開采資料進(jìn)行了分析。所得結(jié)論如下:

    (1)與傳統(tǒng)的SBAS-InSAR方法相比,本研究方法在充分利用高相干點目標(biāo)的基礎(chǔ)上,依靠FaSHPS方法對研究區(qū)域內(nèi)的分布式目標(biāo)進(jìn)行提取,并通過特征值分解對其相位進(jìn)行優(yōu)化,之后將二者聯(lián)合進(jìn)行時序處理來獲取地表形變信息,可有效提升研究區(qū)地表監(jiān)測點密度和形變監(jiān)測精度。

    (2)在監(jiān)測時段內(nèi),該方法監(jiān)測到研究區(qū)地表形成了南北相鄰的2處明顯的沉降區(qū)域,最大沉降約400 mm。與井下開采資料對比分析表明,該方法監(jiān)測到的地表沉降區(qū)域的位置和范圍較為準(zhǔn)確,且沉降區(qū)域的形成及發(fā)展能與地下工作面開采進(jìn)程相對應(yīng),可為開采沉陷影響邊界確定及其機(jī)理研究提供參考。

    (3)驗證了基于DS-InSAR的時序地表形變監(jiān)測方法的可靠性,并利用該方法獲取了礦區(qū)06A工作面開采期間南部沉降區(qū)域沿06A工作面走向和傾向主斷面的動態(tài)沉降曲線,對該區(qū)域缺少的水準(zhǔn)實測數(shù)據(jù)起到了很好的補充作用,并將沉降監(jiān)測結(jié)果進(jìn)一步處理,得到了沿06A工作面走向主斷面的動態(tài)傾斜和曲率變形曲線,可為地表動態(tài)變形規(guī)律研究及形變控制提供依據(jù)。

    猜你喜歡
    時序盆地礦區(qū)
    時序坐標(biāo)
    基于Sentinel-2時序NDVI的麥冬識別研究
    盆地是怎樣形成的
    加納Amanforom礦區(qū)Ⅲ號隱伏金礦帶的發(fā)現(xiàn)與評價
    加納Amanforom礦區(qū)Ⅲ號隱伏金礦帶的發(fā)現(xiàn)與評價
    湖北省??悼h堰邊上礦區(qū)發(fā)現(xiàn)超大型磷礦
    廣東省蕉嶺縣作壁坑礦區(qū)探明超大型銣礦
    二疊盆地Wolfcamp統(tǒng)致密油成藏特征及主控因素
    一種毫米波放大器時序直流電源的設(shè)計
    電子制作(2016年15期)2017-01-15 13:39:08
    楚雄盆地扭動構(gòu)造及其演化
    国产黄色视频一区二区在线观看| 国产一级毛片七仙女欲春2| 亚洲国产精品sss在线观看| 国产成人精品婷婷| 亚洲av免费高清在线观看| 国产精品伦人一区二区| 国产色爽女视频免费观看| 伦理电影大哥的女人| 男人舔奶头视频| 天堂√8在线中文| 最后的刺客免费高清国语| 久久久久久久国产电影| 亚洲最大成人中文| 少妇丰满av| 欧美日韩视频高清一区二区三区二| 国产人妻一区二区三区在| 亚洲国产精品国产精品| 国产有黄有色有爽视频| 我的女老师完整版在线观看| 国产亚洲精品久久久com| 色尼玛亚洲综合影院| 日韩,欧美,国产一区二区三区| 日日啪夜夜撸| 超碰av人人做人人爽久久| 亚洲人与动物交配视频| 中国美白少妇内射xxxbb| 久久热精品热| 激情 狠狠 欧美| 亚洲欧美中文字幕日韩二区| a级一级毛片免费在线观看| 国产成人精品久久久久久| 七月丁香在线播放| 2021少妇久久久久久久久久久| 精品亚洲乱码少妇综合久久| .国产精品久久| 成人国产麻豆网| 久久6这里有精品| 女人被狂操c到高潮| 国产亚洲精品av在线| 男女边摸边吃奶| 精品久久久久久久久av| 人妻少妇偷人精品九色| 久久亚洲国产成人精品v| 日韩人妻高清精品专区| 91久久精品国产一区二区三区| 国产片特级美女逼逼视频| 男女那种视频在线观看| 中文天堂在线官网| 嫩草影院精品99| 国产中年淑女户外野战色| 久久这里有精品视频免费| 一区二区三区高清视频在线| 久久久精品欧美日韩精品| 欧美xxxx性猛交bbbb| 国产成人精品一,二区| 熟女电影av网| 九色成人免费人妻av| 好男人在线观看高清免费视频| 国产一区有黄有色的免费视频 | 国产一区二区三区av在线| 亚洲最大成人av| 99热6这里只有精品| 赤兔流量卡办理| 亚洲精品亚洲一区二区| 嫩草影院新地址| 成人二区视频| 三级男女做爰猛烈吃奶摸视频| 99久久精品热视频| 一本一本综合久久| 国产午夜精品一二区理论片| 久久99蜜桃精品久久| 偷拍熟女少妇极品色| 国产乱来视频区| 亚洲精品国产成人久久av| 一本久久精品| 成人午夜精彩视频在线观看| 国产精品国产三级专区第一集| 噜噜噜噜噜久久久久久91| 人妻系列 视频| 国产亚洲精品av在线| 中文欧美无线码| 亚洲av一区综合| 免费看av在线观看网站| 国产高清不卡午夜福利| 亚洲18禁久久av| 亚洲精品日本国产第一区| 亚洲三级黄色毛片| 国产爱豆传媒在线观看| 综合色丁香网| 午夜精品国产一区二区电影 | 亚洲欧美一区二区三区黑人 | 日韩,欧美,国产一区二区三区| 26uuu在线亚洲综合色| 国产精品久久久久久av不卡| 男女那种视频在线观看| 成人毛片a级毛片在线播放| 精品久久国产蜜桃| 青春草国产在线视频| 亚洲性久久影院| 黄片wwwwww| 在线观看免费高清a一片| 欧美日韩亚洲高清精品| 99视频精品全部免费 在线| 美女黄网站色视频| 51国产日韩欧美| 欧美日本视频| 美女主播在线视频| 亚洲综合精品二区| 少妇的逼水好多| 国产精品国产三级专区第一集| 日韩精品有码人妻一区| 超碰av人人做人人爽久久| 特大巨黑吊av在线直播| 日本-黄色视频高清免费观看| 中文精品一卡2卡3卡4更新| 成年人午夜在线观看视频 | 免费av毛片视频| 春色校园在线视频观看| 日韩成人伦理影院| 亚洲怡红院男人天堂| 亚洲精品国产av成人精品| 日韩视频在线欧美| 亚洲成人精品中文字幕电影| 一区二区三区乱码不卡18| 天堂网av新在线| 午夜福利高清视频| 乱人视频在线观看| 亚洲av在线观看美女高潮| 亚洲电影在线观看av| 午夜激情久久久久久久| 麻豆av噜噜一区二区三区| 人妻制服诱惑在线中文字幕| 欧美精品国产亚洲| 成人二区视频| 国产成人精品久久久久久| av一本久久久久| 干丝袜人妻中文字幕| 欧美激情久久久久久爽电影| 国产国拍精品亚洲av在线观看| 麻豆乱淫一区二区| 一个人看的www免费观看视频| 日日撸夜夜添| 久久久久性生活片| 成年版毛片免费区| 有码 亚洲区| 国产精品一二三区在线看| 日韩,欧美,国产一区二区三区| 99热网站在线观看| 水蜜桃什么品种好| 国产精品.久久久| 亚洲av男天堂| 日韩一区二区三区影片| 国产毛片a区久久久久| 久久久久久伊人网av| 久久久久网色| 亚洲性久久影院| 九九久久精品国产亚洲av麻豆| 国产亚洲精品久久久com| 国产午夜精品一二区理论片| 3wmmmm亚洲av在线观看| 在线 av 中文字幕| 久久鲁丝午夜福利片| 免费av毛片视频| 熟女人妻精品中文字幕| 久久久久久久久久久免费av| 激情五月婷婷亚洲| 精品人妻一区二区三区麻豆| 亚洲欧美精品专区久久| 亚洲精品影视一区二区三区av| 国产成人精品一,二区| 免费高清在线观看视频在线观看| 婷婷色综合www| 国产午夜精品一二区理论片| 国产成人免费观看mmmm| 99热6这里只有精品| 亚洲成色77777| 在线观看三级黄色| 久久久久久久国产电影| 大话2 男鬼变身卡| 一区二区三区四区激情视频| h视频一区二区三区| 观看av在线不卡| 夜夜骑夜夜射夜夜干| 久久久欧美国产精品| 看非洲黑人一级黄片| 免费日韩欧美在线观看| 日韩熟女老妇一区二区性免费视频| 日韩中文字幕视频在线看片| 中文字幕亚洲精品专区| 亚洲激情五月婷婷啪啪| 精品久久久精品久久久| 免费看不卡的av| 26uuu在线亚洲综合色| 中文乱码字字幕精品一区二区三区| 亚洲成av片中文字幕在线观看 | 999久久久国产精品视频| 免费观看av网站的网址| 午夜免费鲁丝| 精品人妻在线不人妻| 美女主播在线视频| 亚洲五月色婷婷综合| 90打野战视频偷拍视频| 欧美国产精品一级二级三级| av在线app专区| 亚洲av中文av极速乱| 午夜福利视频精品| 久久这里有精品视频免费| 看非洲黑人一级黄片| 曰老女人黄片| 国产av精品麻豆| 国产视频首页在线观看| 天天影视国产精品| 日韩制服骚丝袜av| 欧美亚洲日本最大视频资源| 视频在线观看一区二区三区| 老司机亚洲免费影院| 深夜精品福利| 色哟哟·www| 亚洲精品国产av蜜桃| 精品人妻偷拍中文字幕| 日日爽夜夜爽网站| 国产精品不卡视频一区二区| 久久久国产一区二区| 欧美日韩亚洲高清精品| 汤姆久久久久久久影院中文字幕| 色婷婷久久久亚洲欧美| 精品国产一区二区三区久久久樱花| av网站在线播放免费| 最近2019中文字幕mv第一页| 久久久欧美国产精品| 黄片无遮挡物在线观看| 日日爽夜夜爽网站| 亚洲伊人色综图| 国产成人精品在线电影| 99久久中文字幕三级久久日本| 女性被躁到高潮视频| 欧美97在线视频| 国产免费福利视频在线观看| 免费久久久久久久精品成人欧美视频| 9热在线视频观看99| 亚洲国产看品久久| 午夜福利网站1000一区二区三区| 多毛熟女@视频| 国产无遮挡羞羞视频在线观看| 日韩,欧美,国产一区二区三区| 精品国产一区二区久久| 日韩欧美精品免费久久| 精品亚洲成a人片在线观看| 永久免费av网站大全| 国产亚洲午夜精品一区二区久久| 在线免费观看不下载黄p国产| 国产精品不卡视频一区二区| 69精品国产乱码久久久| 亚洲欧美一区二区三区久久| 精品视频人人做人人爽| 久久免费观看电影| 人人妻人人添人人爽欧美一区卜| 久久久久久久国产电影| 欧美日韩综合久久久久久| 各种免费的搞黄视频| 乱人伦中国视频| 香蕉丝袜av| 一级片'在线观看视频| 一级黄片播放器| 精品酒店卫生间| 国产精品亚洲av一区麻豆 | 午夜91福利影院| 91成人精品电影| 国产精品av久久久久免费| av网站在线播放免费| 亚洲精品成人av观看孕妇| 99re6热这里在线精品视频| 日韩伦理黄色片| 午夜老司机福利剧场| 男女无遮挡免费网站观看| 视频在线观看一区二区三区| 午夜免费鲁丝| 欧美成人午夜免费资源| 肉色欧美久久久久久久蜜桃| 熟女少妇亚洲综合色aaa.| 麻豆精品久久久久久蜜桃| 看十八女毛片水多多多| 最近最新中文字幕免费大全7| 女性被躁到高潮视频| 久久亚洲国产成人精品v| 日韩伦理黄色片| 亚洲av免费高清在线观看| 国产精品不卡视频一区二区| 国产极品粉嫩免费观看在线| 一本久久精品| kizo精华| 精品一区二区三区四区五区乱码 | 美女福利国产在线| 999久久久国产精品视频| 美女高潮到喷水免费观看| 国产成人91sexporn| 久久久久网色| 午夜福利,免费看| 精品国产乱码久久久久久小说| 高清不卡的av网站| 女人高潮潮喷娇喘18禁视频| 国产成人精品无人区| 国产精品女同一区二区软件| 三上悠亚av全集在线观看| 制服丝袜香蕉在线| 国产日韩一区二区三区精品不卡| 国产欧美日韩一区二区三区在线| 国产免费福利视频在线观看| 国产片内射在线| 午夜精品国产一区二区电影| 久久人妻熟女aⅴ| 久久热在线av| 国产亚洲精品第一综合不卡| 欧美亚洲 丝袜 人妻 在线| 黄色毛片三级朝国网站| 国产日韩欧美在线精品| 亚洲熟女精品中文字幕| 一级毛片我不卡| 免费在线观看完整版高清| 免费播放大片免费观看视频在线观看| 亚洲四区av| a级片在线免费高清观看视频| 免费在线观看黄色视频的| 日韩免费高清中文字幕av| 国产精品 欧美亚洲| 中国国产av一级| 999久久久国产精品视频| 亚洲av日韩在线播放| 欧美精品国产亚洲| 精品国产超薄肉色丝袜足j| 欧美最新免费一区二区三区| √禁漫天堂资源中文www| 蜜桃国产av成人99| 欧美激情高清一区二区三区 | 中文字幕色久视频| 建设人人有责人人尽责人人享有的| 90打野战视频偷拍视频| 国产精品欧美亚洲77777| 18禁观看日本| 寂寞人妻少妇视频99o| 精品第一国产精品| 欧美人与善性xxx| 久久久久久伊人网av| 日韩一区二区三区影片| 中文字幕另类日韩欧美亚洲嫩草| 丰满迷人的少妇在线观看| 日日啪夜夜爽| 2022亚洲国产成人精品| 好男人视频免费观看在线| 一区二区日韩欧美中文字幕| 国产麻豆69| 亚洲五月色婷婷综合| 国产精品av久久久久免费| 9色porny在线观看| 国产一区有黄有色的免费视频| 国产黄频视频在线观看| 最近手机中文字幕大全| 国产成人一区二区在线| 亚洲av福利一区| 欧美 日韩 精品 国产| 天天操日日干夜夜撸| 一二三四中文在线观看免费高清| 久久狼人影院| 亚洲av成人精品一二三区| 在线免费观看不下载黄p国产| 人妻一区二区av| 看免费成人av毛片| 26uuu在线亚洲综合色| 在线观看www视频免费| 精品人妻偷拍中文字幕| 91精品国产国语对白视频| 日韩在线高清观看一区二区三区| 99热网站在线观看| 一级黄片播放器| 99精国产麻豆久久婷婷| 久久久久久久久免费视频了| 18禁动态无遮挡网站| 麻豆av在线久日| 国产精品熟女久久久久浪| 人人妻人人澡人人爽人人夜夜| 亚洲欧洲精品一区二区精品久久久 | 久久精品aⅴ一区二区三区四区 | 免费观看性生交大片5| 成年女人在线观看亚洲视频| 精品国产一区二区久久| 99热全是精品| 永久网站在线| 免费观看无遮挡的男女| 欧美变态另类bdsm刘玥| 最近最新中文字幕免费大全7| 蜜桃国产av成人99| 永久免费av网站大全| 精品福利永久在线观看| 精品国产国语对白av| 国产伦理片在线播放av一区| 免费不卡的大黄色大毛片视频在线观看| 婷婷成人精品国产| 精品人妻偷拍中文字幕| 一级毛片我不卡| 亚洲久久久国产精品| 日韩av免费高清视频| 亚洲伊人色综图| 26uuu在线亚洲综合色| 黑人巨大精品欧美一区二区蜜桃| 免费不卡的大黄色大毛片视频在线观看| 国产日韩欧美在线精品| 精品一区在线观看国产| 国产精品二区激情视频| 男男h啪啪无遮挡| 99国产综合亚洲精品| 日本欧美国产在线视频| 一级,二级,三级黄色视频| 久久久久国产精品人妻一区二区| 香蕉国产在线看| 纯流量卡能插随身wifi吗| 69精品国产乱码久久久| 国产亚洲一区二区精品| 欧美精品一区二区大全| 亚洲国产最新在线播放| 人妻一区二区av| 一边亲一边摸免费视频| 精品亚洲乱码少妇综合久久| 日本av免费视频播放| 免费黄色在线免费观看| 免费观看av网站的网址| av卡一久久| 岛国毛片在线播放| 黑人欧美特级aaaaaa片| 成年女人在线观看亚洲视频| 最近最新中文字幕免费大全7| 久久久久久人人人人人| 久久99精品国语久久久| 亚洲成国产人片在线观看| 夜夜骑夜夜射夜夜干| 欧美日韩综合久久久久久| av在线老鸭窝| 黄频高清免费视频| 久久韩国三级中文字幕| 久久av网站| 国产精品欧美亚洲77777| 一级爰片在线观看| 亚洲国产av新网站| 性色av一级| 一级片免费观看大全| 高清欧美精品videossex| 大香蕉久久网| 成人18禁高潮啪啪吃奶动态图| 日韩电影二区| 成人国语在线视频| 18禁动态无遮挡网站| 免费人妻精品一区二区三区视频| 99久久精品国产国产毛片| 五月开心婷婷网| 国产综合精华液| 少妇被粗大的猛进出69影院| 高清av免费在线| 国产成人精品一,二区| 国产成人91sexporn| 久久人人97超碰香蕉20202| 亚洲熟女精品中文字幕| 99九九在线精品视频| av女优亚洲男人天堂| 熟女电影av网| 在线观看免费高清a一片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男女啪啪激烈高潮av片| 免费高清在线观看视频在线观看| 日本欧美国产在线视频| 国产av码专区亚洲av| 亚洲成人手机| 亚洲精品美女久久av网站| 日本爱情动作片www.在线观看| 高清视频免费观看一区二区| 两个人看的免费小视频| 黄片小视频在线播放| 亚洲av电影在线观看一区二区三区| 一区二区日韩欧美中文字幕| 日韩中文字幕视频在线看片| 日日啪夜夜爽| 亚洲欧洲国产日韩| 天天躁夜夜躁狠狠躁躁| 久久99蜜桃精品久久| 一级黄片播放器| 午夜日本视频在线| av女优亚洲男人天堂| 高清黄色对白视频在线免费看| 2021少妇久久久久久久久久久| 在线天堂最新版资源| 一区二区三区乱码不卡18| 最近中文字幕2019免费版| 日本爱情动作片www.在线观看| 五月开心婷婷网| 在线 av 中文字幕| 免费观看a级毛片全部| 一个人免费看片子| 两个人看的免费小视频| 国产av码专区亚洲av| 亚洲欧美日韩另类电影网站| 日韩伦理黄色片| 美女午夜性视频免费| 天堂8中文在线网| 夫妻性生交免费视频一级片| 狂野欧美激情性bbbbbb| 国产极品粉嫩免费观看在线| 久久精品熟女亚洲av麻豆精品| 日本av手机在线免费观看| 你懂的网址亚洲精品在线观看| 麻豆乱淫一区二区| 99热网站在线观看| 国产精品不卡视频一区二区| 天天躁日日躁夜夜躁夜夜| 亚洲精品日韩在线中文字幕| 亚洲精品一区蜜桃| 国产视频首页在线观看| 久久久久久久久久人人人人人人| 欧美精品亚洲一区二区| 欧美bdsm另类| 久久这里有精品视频免费| 高清不卡的av网站| 欧美精品人与动牲交sv欧美| 久久久久久久国产电影| 国产精品三级大全| 久久久久久久久久久免费av| 最近2019中文字幕mv第一页| 精品一区二区三区四区五区乱码 | 亚洲国产欧美网| 日韩一本色道免费dvd| 天天躁夜夜躁狠狠躁躁| √禁漫天堂资源中文www| 精品酒店卫生间| 人妻少妇偷人精品九色| 日韩不卡一区二区三区视频在线| 日韩一区二区视频免费看| 高清不卡的av网站| 久久久久久久久久久久大奶| 欧美成人午夜免费资源| 亚洲欧洲精品一区二区精品久久久 | 久久久精品免费免费高清| 国产免费视频播放在线视频| 日韩一本色道免费dvd| 久久精品国产亚洲av天美| 精品一区在线观看国产| 大陆偷拍与自拍| 久久久久国产精品人妻一区二区| 国产精品成人在线| 如何舔出高潮| 如何舔出高潮| 叶爱在线成人免费视频播放| av电影中文网址| 一二三四在线观看免费中文在| 国产一级毛片在线| 欧美日韩一级在线毛片| 一区二区三区激情视频| 在线观看一区二区三区激情| 中文字幕色久视频| 欧美在线黄色| 国产激情久久老熟女| 亚洲欧美色中文字幕在线| 国产在线视频一区二区| xxxhd国产人妻xxx| 黄色一级大片看看| 亚洲成人手机| 亚洲成av片中文字幕在线观看 | 七月丁香在线播放| 国产精品久久久av美女十八| 中文精品一卡2卡3卡4更新| 一区二区日韩欧美中文字幕| 成人18禁高潮啪啪吃奶动态图| 大香蕉久久网| av不卡在线播放| 国产黄频视频在线观看| 亚洲精品国产av蜜桃| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品一二三| 男女边吃奶边做爰视频| 捣出白浆h1v1| 2021少妇久久久久久久久久久| 在线观看免费日韩欧美大片| 狠狠精品人妻久久久久久综合| 中文精品一卡2卡3卡4更新| 熟女av电影| 男男h啪啪无遮挡| 国产成人精品无人区| 久久精品久久久久久久性| 久久久久久伊人网av| 亚洲国产精品999| 大片免费播放器 马上看| 边亲边吃奶的免费视频| 国产无遮挡羞羞视频在线观看| 国产欧美日韩综合在线一区二区| 亚洲伊人色综图| 最近最新中文字幕免费大全7| 在线观看免费视频网站a站| 国产综合精华液| 精品国产一区二区三区久久久樱花| 亚洲欧美一区二区三区国产| 婷婷色综合大香蕉| 老女人水多毛片| 色哟哟·www| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品熟女久久久久浪| 亚洲精品自拍成人| 秋霞在线观看毛片| 最新中文字幕久久久久| av在线老鸭窝| av天堂久久9| 女性生殖器流出的白浆| 国产精品熟女久久久久浪| 免费高清在线观看视频在线观看| 亚洲国产成人一精品久久久| 777米奇影视久久| 国产免费现黄频在线看| 亚洲精品在线美女| 国产精品三级大全| 一级毛片黄色毛片免费观看视频|