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

    地面沉降-回彈及地下水位波動(dòng)的InSAR長(zhǎng)時(shí)序監(jiān)測(cè)——以德州市為例

    2014-09-26 02:23:26葛大慶殷躍平郭小方
    自然資源遙感 2014年1期
    關(guān)鍵詞:德州深層基線

    葛大慶,殷躍平,王 艷,張 玲,郭小方,王 毅

    (1.中國(guó)地質(zhì)大學(xué)(北京)水資源與環(huán)境學(xué)院,北京 100083;2.中國(guó)國(guó)土資源航空物探遙感中心,北京 100083;3.中國(guó)地質(zhì)環(huán)境監(jiān)測(cè)院,北京 100081)

    0 引言

    長(zhǎng)期過(guò)量開(kāi)采地下水,特別是深層地下水,是目前我國(guó)一些平原和三角洲地區(qū)產(chǎn)生地面沉降的主要原因[1]。地下水位的持續(xù)下降改變了原有的補(bǔ)徑排特征,形成了區(qū)域性地下水位降落漏斗,地面沉降隨之產(chǎn)生并擴(kuò)大。我國(guó)已形成了以華北平原、長(zhǎng)江三角洲及汾渭斷陷盆地為主的三大地面沉降嚴(yán)重區(qū),華北平原已出現(xiàn)近20個(gè)大型復(fù)合水位降落漏斗,地面沉降呈現(xiàn)出區(qū)域性且連片分布的特征[2-3]。

    雷達(dá)干涉測(cè)量(InSAR)技術(shù)為區(qū)域性地面沉降的快速準(zhǔn)確和連續(xù)監(jiān)測(cè)提供了重要手段。以永久性散射體干涉測(cè)量(permanent scatterer InSAR,PSIn-SAR)[4]為代表的相干目標(biāo)(coherent target,CT)時(shí)間序列分析技術(shù)克服了差分InSAR技術(shù)所面臨的失相干和大氣相位延遲等問(wèn)題,通過(guò)對(duì)相干目標(biāo)的差分干涉相位序列進(jìn)行時(shí)序分析,逐個(gè)估計(jì)和分離地形相位、形變相位、大氣相位和噪聲等信息,進(jìn)而解算相干目標(biāo)的形變速率和累積形變量。與之類似的還有短基線集(small baseline subset,SBAS)[5]以及干涉測(cè)量點(diǎn)目標(biāo)分析(interferometric point target analysis,IPTA)[6]等技術(shù),這些技術(shù)已在地面沉降、滑坡、開(kāi)采沉陷、地震和冰川滑移等多尺度緩慢地表形變研究中得以發(fā)展和應(yīng)用[7-9]。

    德州市位于山東省北部,是華北平原地面沉降和深層地下水位降落漏斗極為典型的地區(qū)之一。其地下水位降落漏斗在1978年前已經(jīng)形成,是目前華北平原水位降落深度最大的漏斗,與衡水、滄州等漏斗連成一片。該區(qū)的地面沉降在1989年以前就已經(jīng)出現(xiàn),至2010年累積沉降量超過(guò)1 m。與地面沉降對(duì)應(yīng)的是,德州市地下水位呈現(xiàn)出多年連續(xù)下降和年內(nèi)季節(jié)性波動(dòng)的特征[10]。由于地面沉降與地下水位下降相關(guān),因此在德州市深、淺層地下水位均呈現(xiàn)出季節(jié)性波動(dòng)的情況下,地面沉降場(chǎng)的時(shí)空演變特征值得研究,這包括:沉降場(chǎng)的空間變化特征,即沉降中心的移動(dòng)及其變化幅度;沉降場(chǎng)的時(shí)間變化特征,即季節(jié)性波動(dòng),包括是否會(huì)出現(xiàn)地面回彈,回彈幅度多大,地面沉降與地下水位波動(dòng)是否同步等方面。針對(duì)上述問(wèn)題,本文以相干目標(biāo)InSAR時(shí)序分析方法為主,利用2004年1月—2010年10月近7 a的ENVISAT衛(wèi)星數(shù)據(jù)研究分析德州地區(qū)地面沉降-回彈的動(dòng)態(tài)變化特征及其與地下水位波動(dòng)的對(duì)應(yīng)關(guān)系,以揭示該地區(qū)地面沉降季節(jié)性變化的產(chǎn)生原因和誘發(fā)機(jī)制。

    1 相干目標(biāo)InSAR時(shí)序分析方法

    InSAR時(shí)序分析主要有PSInSAR和SBAS兩類方法。根據(jù)干涉像對(duì)的組合模式,前者以單一主影像構(gòu)成差分干涉圖序列,后者則利用短基線原則以多個(gè)主影像構(gòu)成若干個(gè)干涉紋圖子集。單一主影像的干涉紋圖序列處理主要顧及大氣相位和非線性形變的統(tǒng)計(jì)模型,有利于對(duì)大氣相位的準(zhǔn)確估計(jì),但對(duì)數(shù)據(jù)量要求過(guò)高,且受相干性影響而對(duì)大變形的監(jiān)測(cè)能力不足;多個(gè)主影像則降低了對(duì)數(shù)據(jù)量的要求,將時(shí)空基線小于一定閾值的干涉像對(duì)組合,生成干涉紋圖序列,增加對(duì)同一信號(hào)的采樣密度,有利于準(zhǔn)確求解形變速率,在數(shù)據(jù)量較小的情況下也可以實(shí)現(xiàn)。在選擇相干目標(biāo)的方式上,前者以點(diǎn)目標(biāo)為主,后者兼顧點(diǎn)目標(biāo)和分布式目標(biāo)[5]。本文集成了2種方法的優(yōu)點(diǎn),以短基線為準(zhǔn)則構(gòu)成差分干涉相位圖,利用點(diǎn)目標(biāo)識(shí)別算法提取相干目標(biāo)。以Delanay三角網(wǎng)連接相鄰的相干像元,用二維周期圖估計(jì)點(diǎn)間形變速率和高程誤差改正。以此為基礎(chǔ),根據(jù)大氣相位特征,利用奇異值分解算法解算空間濾波后的殘余相位,并對(duì)解算結(jié)果進(jìn)行時(shí)域?yàn)V波,以求解非線性形變,得到不同時(shí)刻的累積變形量,再利用線性回歸法求解每個(gè)相干目標(biāo)的變形速率。

    1.1 干涉相位時(shí)序分析的基本模型

    對(duì)于給定的M景雷達(dá)數(shù)據(jù),在短基線條件下生成N個(gè)差分干涉圖(一般情況下,N>M),對(duì)于其中的任一差分干涉圖k(k=1,2,…,N),其相位模型為

    顯然,待求解量d與InSAR參數(shù)不直接相關(guān),但受其他分量的影響。對(duì)于大氣分量,其本質(zhì)上是2次大氣延遲影響的疊加,在空間上具有一定的相關(guān)性[4]。大氣波動(dòng)在時(shí)間上的隨機(jī)特性決定了2次互差也為隨機(jī)量。需要重點(diǎn)區(qū)分的是大氣相位與軌道誤差,雖然二者都具有空間低頻的特性,但各自的影響范圍有所不同。基線誤差多是全局性的,而大氣影響半徑有限,對(duì)于大范圍處理,二者應(yīng)區(qū)別對(duì)待。高程誤差與垂直基線相關(guān),垂直基線是變量。變形量與干涉圖的時(shí)間基線相關(guān),時(shí)間間隔是變量,相對(duì)變形的速率與時(shí)間相關(guān)。

    1.2 地表形變參數(shù)的提取算法

    1.2.1 線性分量的估計(jì)

    考慮到大氣的空間相關(guān)性,如果對(duì)相鄰2點(diǎn)(小于大氣相關(guān)距離)求差,則可削弱大氣影響。相鄰相干目標(biāo)i和j差分干涉相位的互差△為

    由于相鄰目標(biāo)間的相對(duì)變形速率與時(shí)間間隔相關(guān),干涉圖的時(shí)間基線為自變量,因而可將形變量d分解為線性形變速率和非線性形變量2部分。如此,時(shí)序分析過(guò)程則轉(zhuǎn)為參數(shù)估計(jì),即

    式中:CB為與垂直基線相關(guān)的系數(shù);T為時(shí)間基線;△ε為相對(duì)高程誤差;△υ為相對(duì)形變速率;μN(yùn)L為非線性形變量;α為大氣相位;n為噪聲。

    對(duì)于纏繞相位,其周期數(shù)由相對(duì)高程誤差和線性形變速率共同決定,采用二維參數(shù)估計(jì)方法可實(shí)現(xiàn)對(duì)相位周期的估算。由此,構(gòu)建目標(biāo)函數(shù)

    并將其從相位互差中減去,得到殘余相位,即

    顯然,在準(zhǔn)確估計(jì)目標(biāo)函數(shù)的參數(shù)△ε和△υ時(shí),殘余相位將最小化。

    由于相鄰2點(diǎn)存在多個(gè)干涉圖,即存在N維B和T,因而,對(duì)△ε和△v的估計(jì)使模型相關(guān)系數(shù)最大化,即

    在纏繞相位條件下,△ε和△v為周期函數(shù)的二維頻率,可利用二維周期圖估計(jì)使模型相關(guān)系數(shù)最大化;在解纏相位條件下,則轉(zhuǎn)換為二維線性函數(shù),利用參數(shù)估計(jì)方法可實(shí)現(xiàn)對(duì)△ε和△v的求解。

    在相鄰2點(diǎn)解算時(shí)需對(duì)所有相干目標(biāo)建立連接關(guān)系,整體求解形變場(chǎng)的速率和高程誤差改正量。這一目標(biāo)的實(shí)現(xiàn)可通過(guò)建立Delanay三角網(wǎng)的方法完成,也可以利用冗余網(wǎng)構(gòu)建更為復(fù)雜的連接關(guān)系,強(qiáng)化對(duì)待解算方程組的約束。利用鄰近法則將所有距離滿足大氣相關(guān)距離的相干目標(biāo)連接起來(lái),在求解完成相鄰點(diǎn)間的互差后,通過(guò)最小二乘或加權(quán)平均的方法求解每個(gè)目標(biāo)相對(duì)于參考點(diǎn)的變形速率,即總體形變速率場(chǎng)。

    1.2.2 非線性分量的估計(jì)

    對(duì)于具有顯著非線性形變過(guò)程,仍需對(duì)殘余相位進(jìn)行更為復(fù)雜的處理,以提取非線性形變量。從差分相位中去除高程誤差估計(jì)值和線性速率后,殘余相位為

    由于已解算出整周相位,因而殘余相位為相位主值,其大小

    這里,首先要實(shí)現(xiàn)對(duì)單個(gè)差分干涉圖中殘余相位的濾波處理。由于大氣表現(xiàn)出空間低頻特性,而非線性變形的空間變化范圍較小,相對(duì)大氣而言表現(xiàn)為高通特性。因而,對(duì)殘余相位圖進(jìn)行空間低通濾波可以進(jìn)一步弱化大氣影響。在短基線集條件下求解的非線性形變量,不同于PSInSAR經(jīng)典處理策略,不直接對(duì)大氣進(jìn)行估計(jì),而是通過(guò)濾波減弱大氣的影響。殘余相位經(jīng)空間高通濾波后,大氣分量已經(jīng)減弱。此時(shí),進(jìn)一步的處理是求解與雷達(dá)數(shù)據(jù)對(duì)應(yīng)的不同時(shí)刻的非線性變形量,包含噪聲等影響。根據(jù)主輔影像的關(guān)系,可將解纏后的殘余相位分解為

    式中:p和q分別表示生成第k(k=1,…,N)景差分干涉圖的主輔影像的獲取時(shí)間。由于采用短基線原則,在求解每個(gè)時(shí)刻對(duì)應(yīng)的殘余相位時(shí)出現(xiàn)秩虧方程組,解決這一問(wèn)題的辦法就是奇異值分解法(singular value decomposition,SVD)[3]。在此基礎(chǔ)上,對(duì)M景雷達(dá)數(shù)據(jù)的殘余相位進(jìn)行時(shí)域低通濾波處理,以提取最終的非線性形變量,可表示為

    在求解出非線性形變量后,每個(gè)相干目標(biāo)對(duì)應(yīng)的相變序列為

    式中vest為線性形變速率。

    1.2.3 平均形變速率的重新估計(jì)

    對(duì)于顯著的非線性形變過(guò)程,利用前述二維線性模型估計(jì)出來(lái)的線性速率往往偏小。在恢復(fù)整個(gè)形變序列的基礎(chǔ)上,可根據(jù)對(duì)應(yīng)的每個(gè)目標(biāo)的形變序列,利用最小二乘法重新估計(jì)形變速率,即

    式中d和T分別為相對(duì)起始時(shí)刻的形變量和時(shí)間間隔。

    2 雷達(dá)數(shù)據(jù)及InSAR時(shí)序分析處理

    研究中共獲取2004年1月—2010年10月近7 a的ENVISAT衛(wèi)星ASAR數(shù)據(jù)。該數(shù)據(jù)覆蓋了德州地區(qū)主要沉降區(qū),總數(shù)據(jù)量為46景(Track-447,F(xiàn)rame-2853),其基線分布如圖1所示。

    圖1 差分干涉圖序列的短基線關(guān)系Fig.1 Small baseline of the differential interferograms series

    需要說(shuō)明的是,利用幅度離散指數(shù)[1]、相干系數(shù)[5]以及子視相關(guān)[6]均可以很好地識(shí)別出相干目標(biāo)候選點(diǎn)。對(duì)于候選點(diǎn),只有在滿足模型相關(guān)函數(shù)最大化約束條件下才能稱為時(shí)序分析中的相干目標(biāo)。實(shí)際上,一些目標(biāo)后向散射特性穩(wěn)定,但非線性變化強(qiáng)烈,不滿足該條件。因而,可通過(guò)降低相關(guān)系數(shù)的閾值,將其納入到形變解算網(wǎng)絡(luò)中,以增加監(jiān)測(cè)點(diǎn)的密度。

    時(shí)空基線閾值的設(shè)置取決于SAR系統(tǒng)參數(shù)。空間基線主要考慮雷達(dá)系統(tǒng)臨界基線,而時(shí)間基線則應(yīng)根據(jù)變形幅度大小,短時(shí)間基線相當(dāng)于降低了形變場(chǎng)的梯度,便于相位解纏和變形參數(shù)估計(jì)。因而,本研究中設(shè)定的短基線條件為:時(shí)間基線<450 d,空間基線<300 m。

    圍繞德州地區(qū)地面沉降和地下水位漏斗的基本分布特征,分別針對(duì)區(qū)域性大范圍沉降和重點(diǎn)沉降中心進(jìn)行了2類處理:利用2007—2010年4 a間的數(shù)據(jù)提取了德州地區(qū)100 km×100 km平均地面沉降速率;分析全區(qū)地面沉降和地下水位變化關(guān)系,并確定時(shí)序分析過(guò)程中穩(wěn)定的地面參考點(diǎn)。在此基礎(chǔ)上,對(duì)德州城區(qū)內(nèi)20 km×20 km范圍的主要沉降中心的7 a數(shù)據(jù)進(jìn)行了時(shí)序分析處理,提取每個(gè)相干目標(biāo)的沉降率和累積沉降量序列,用以詳細(xì)分析多年和年度內(nèi)沉降中心的時(shí)空變化特征。

    3 區(qū)域性地面沉降與地下水位變化

    3.1 區(qū)域地面沉降分布

    圖2示出德州地區(qū)約50 km×60 km范圍內(nèi)2007年初—2010年底4 a間平均地面沉降速率。

    圖2 德州地區(qū)2007—2010年CTInSAR監(jiān)測(cè)平均地面沉降率Fig.2 Average subsidence velocity of Dezhou area from 2007 to 2010 derived by CTInSAR data

    從區(qū)域分布上來(lái)看,德州市區(qū)沉降較為顯著,沉降區(qū)向西北方向延伸,與河北省景縣和故城相連,構(gòu)成連片沉降區(qū)。區(qū)內(nèi)景縣沉降向西北方向擴(kuò)展,是目前衡水—德州沉降區(qū)的重要組成部分,沉降速率均達(dá)到40~50 mm/a。德州沉降區(qū)向東南方向,分布著3個(gè)明顯的沉降中心,分別為陵縣、臨邑縣和平原縣,其中以陵縣最為突出,中心位于縣城西北菜園村,最大沉降速率達(dá)75 mm/a。臨邑縣沉降中心位于臨盤(pán)采油廠,中心處沉降速率達(dá)50~60 mm/a。

    3.2 與深層地下水位漏斗的比較分析

    InSAR監(jiān)測(cè)結(jié)果準(zhǔn)確揭示了德州地區(qū)區(qū)域性沉降和單個(gè)漏斗的分布。連續(xù)多年的地下水位變化監(jiān)測(cè)結(jié)果表明,德州市深層地下水位降落漏斗是目前魯北地區(qū)最大的漏斗,與周邊的小漏斗相連接,引起區(qū)域性地面沉降[8-9]。深層地下水開(kāi)采(第Ⅳ,Ⅴ含水層組)是該地區(qū)地面沉降的主要原因。以臨邑為例,該沉降中心深層地下水降落漏斗的-30 m等水壓線與德城區(qū)深層漏斗相連,已成為德州漏斗的一部分;-40 m等水壓線包圍臨邑縣城和臨盤(pán)鎮(zhèn),圈閉面積達(dá)192 km2。漏斗中心在臨盤(pán)采油廠一帶,水位埋深超過(guò)59.70 m,水位標(biāo)高-41.82 m左右。圖2中的沉降中心準(zhǔn)確地反映了目前該漏斗的分布位置和影響范圍。由于本區(qū)地下水水位動(dòng)態(tài)屬于較穩(wěn)定的連續(xù)開(kāi)采消耗型,年水位變化趨勢(shì)與開(kāi)采量大小密切相關(guān)。臨邑深層地下水開(kāi)采量達(dá)3.56×104m3/a。其中,第Ⅳ和第Ⅴ含水層組地下水開(kāi)采量約占90%,地下水位多年平均降速2.0 m/a,從而使臨邑深層地下水降落漏斗迅速擴(kuò)展,成為較大的區(qū)域性降落漏斗。

    4 德城區(qū)沉降中心的沉降與回彈

    4.1 德城區(qū)主要沉降中心

    2010年之前連續(xù)7 a的InSAR時(shí)序分析結(jié)果(圖3)表明,德城區(qū)有3個(gè)典型沉降中心,分別位于德城區(qū)陳莊(德棉一廠)、長(zhǎng)莊(華源紡織廠)和肖何莊。其中陳莊沉降中心呈南北向分布,長(zhǎng)約2.2 km,寬1.5 km,中心最大沉降速率達(dá)到55 mm/a;長(zhǎng)莊沉降中心呈北東向,長(zhǎng)為2 km,寬為1.5 km,中心最大速率也超過(guò)50 mm/a。連片分布沉降出現(xiàn)在德城區(qū)東部的付莊和宋官屯鎮(zhèn),這一地區(qū)為近年來(lái)快速發(fā)展的新區(qū)。根據(jù)1991—2010年水準(zhǔn)測(cè)量結(jié)果[11],德城區(qū)沉降中心累計(jì)沉降量為-1 186.9~-636.9 mm,位于德棉一廠(即陳莊漏斗)位置,中心多年平均速率為59.35 mm/a。

    圖3 德城區(qū)2004—2010年CTInSAR監(jiān)測(cè)平均地面沉降率Fig.3 Average subsidence velocity of Decheng district from 2004 to 2010 derived by CTInSAR data

    其中,2005—2006年沉降量為-56 mm;2006—2007年沉降量為-89.0 mm;2007—2010年總沉降量為93.2 mm,年平均沉降量為31.1 mm,與 InSAR監(jiān)測(cè)反映的德城區(qū)陳莊沉降中心的變化完全一致。

    4.2 地面沉降與回彈的時(shí)間效應(yīng)

    由2004—2010年InSAR監(jiān)測(cè)結(jié)果發(fā)現(xiàn):德城區(qū)多個(gè)沉降中心年內(nèi)均表現(xiàn)出季節(jié)性沉降與回彈特征(圖4)。陳莊、長(zhǎng)莊、肖何莊以及德棉二廠4個(gè)典型沉降中心均出現(xiàn)回彈,回彈最顯著的位于陳莊、長(zhǎng)莊和德棉二廠沉降中心。其中陳莊的沉降和回彈幅度均為最大,對(duì)應(yīng)的年份為2007年,2—8月間沉降幅度超過(guò)80 mm,至翌年2月回彈15~20 mm。除德棉二廠外,其他3個(gè)中心2008年沉降幅度均有所減緩,范圍縮小,2009年繼續(xù)減緩,而市區(qū)東部在8月份之后沉降加快,2010年繼續(xù)增加。

    圖4 德城區(qū)典型沉降中心的季節(jié)性沉降與回彈Fig.4 Seasonal subsidence and rebound of the major center in Decheng district

    圖4 所示5個(gè)年度的變化特征表明,以陳莊為代表的沉降中心具有明顯的季節(jié)性變化特征,表現(xiàn)為每年的2—3月開(kāi)始快速沉降,至8—9月達(dá)到最大,此后逐步回彈,至翌年2—3月重新開(kāi)始沉降,沉降與回彈之比約為4∶1。

    5 地面沉降及回彈與地下水位變化的關(guān)系

    5.1 持續(xù)沉降與深層地下水位的關(guān)系

    德州地區(qū)地下水系統(tǒng)分為4個(gè)孔隙含水巖組:第Ⅰ含水巖組潛水層,是主要開(kāi)采層;第Ⅱ含水巖組屬于咸水層,基本不開(kāi)采;第Ⅲ含水巖組,底板埋深450~500 m;第Ⅳ含水巖組,底板埋深800~950 m。第Ⅲ,Ⅳ含水巖組屬于承壓含水層,水頭埋深70~100 m,水質(zhì)普遍較好,是城鎮(zhèn)居民生活和工業(yè)用水的主要開(kāi)采層,德州地下水位降落漏斗屬于該層位,其中心水位(位于陳莊國(guó)棉一廠)埋深超過(guò)120 m[10]。根據(jù)德州深層地下開(kāi)采引發(fā)地面沉降變化閾值的研究[11]:持續(xù)開(kāi)采深層地下水是地面沉降的主要原因,地面沉降與深層水位降落呈現(xiàn)出顯著的線性關(guān)系,深層地下水位每下降1 m,產(chǎn)生的地面沉降約為17.5 mm。該經(jīng)驗(yàn)值與近4 a來(lái)地面測(cè)量結(jié)果一致,速率為30~40 mm/a。

    圖5所示為德城區(qū)202長(zhǎng)觀孔(位于陳莊漏斗區(qū))1996—2010年的深層水位變化情況,表現(xiàn)出持續(xù)下降和年內(nèi)波動(dòng)的特點(diǎn)。德城區(qū)淺層地下水位變化曲線如圖6所示。

    圖5 德城區(qū)202長(zhǎng)觀孔深層地下水位變化曲線Fig.5 Dynamic curve of deep groundwater level of No.202 hole in Decheng district

    圖6 德城區(qū)淺層地下水位變化曲線Fig.6 Dynamic curve of shallow groundwaterlevel in Decheng district

    5.2 季節(jié)性沉降-回彈與地下水位波動(dòng)的關(guān)系

    以陳莊漏斗為例,其動(dòng)態(tài)變化(圖3和圖4)與深、淺層地下水位的季節(jié)性波動(dòng)密切相關(guān)。每年1—5月降水少,農(nóng)業(yè)用水量大,淺層地下水位急劇下降;7—8月雨量增多、河道水量補(bǔ)給和開(kāi)采量相對(duì)減少,地下水埋深逐漸回升,至年底達(dá)到次年最高點(diǎn)。深層承壓水位變化主要取決于開(kāi)采量,年水位變化與開(kāi)采量關(guān)系密切。每年3—7月,深層地下水連續(xù)大量開(kāi)采,漏斗內(nèi)壓力水頭急劇下降;8月—翌年2月由于開(kāi)采量減小,加之大量引入黃河水,水位有所回升,回升幅度一般小于2.0 m。

    開(kāi)采量決定了地面沉降的范圍和強(qiáng)度。目前,在深層水位降落漏斗中心區(qū)(陳莊),開(kāi)采井密度為334眼/km2。2007年以前,由于深層地下水開(kāi)采量大,漏斗內(nèi)壓力水逐年下降,年下降幅度較大,一般2~4 m;2008年以后,由于采取封井措施,深層地下水開(kāi)采量減少,水位呈現(xiàn)回升趨勢(shì)[11]。

    陳莊、長(zhǎng)莊等中心的沉降與回彈的波動(dòng)與上述變化十分吻合。圖7分別給出了陳莊(國(guó)棉一廠)、長(zhǎng)莊(華源紡織廠)、國(guó)棉二廠、肖何莊和宋官屯5個(gè)沉降中心連續(xù)7 a的沉降過(guò)程(監(jiān)測(cè)日期為每年的1月30日),均表現(xiàn)出季節(jié)性的非線性波動(dòng)(紅色線表示初次估計(jì)沉降速率)。以陳莊為例,其年度內(nèi)沉降與回彈與近年深層水位波動(dòng)相對(duì)應(yīng),沉降速率在2008年后出現(xiàn)回落,小于2007年之前。

    圖7 德城區(qū)典型沉降中心2004—2010年地面沉降及深層水位變化序列Fig.7 Subsidence history of coherent target in the typical subsidence center of Decheng district

    從地下水天然補(bǔ)給而言,春旱夏澇的氣候特征決定了德州地區(qū)春季地下水補(bǔ)給量的不足[11],而同時(shí)用水量又進(jìn)一步增大,特別是陳莊、長(zhǎng)莊和國(guó)棉二廠3個(gè)沉降中心,均為棉紡、化工等耗水企業(yè),需要大量采水,而汛期澇季時(shí),地表徑流補(bǔ)給及時(shí),地下水開(kāi)采總量減小,水位出現(xiàn)回升,地面沉降也隨之減緩,并在部分時(shí)段內(nèi)出現(xiàn)回彈的現(xiàn)象。

    6 結(jié)論

    1)2004—2010年間的連續(xù)InSAR監(jiān)測(cè)有效揭示了德州地區(qū)區(qū)域性地面沉降和典型漏斗季節(jié)性地面沉降-回彈的發(fā)展與變化特征。德州市及其周邊主要地下水降落漏斗與地面沉降中心相互一致,證明了深層地下水開(kāi)采是本地區(qū)地面沉降的主要原因。

    2)德城區(qū)陳莊、長(zhǎng)莊等典型漏斗的季節(jié)性沉降-回彈與該地區(qū)深淺層地下水位波動(dòng)相關(guān),地面沉降與地下水位變化表現(xiàn)出良好的相關(guān)關(guān)系,變化時(shí)間基本同步,滯后期約1個(gè)月。地面沉降漏斗呈現(xiàn)出3—8月快速沉降,沉降幅度達(dá)50~80 mm,8月—翌年2月地面回彈,回彈幅度達(dá)10~20 mm的變化特點(diǎn)。

    志謝:感謝山東省魯北工程地質(zhì)勘察院鄒祖光、陳松、楊亞賓以及河北水文地質(zhì)工程地質(zhì)四隊(duì)杜興明、田小偉等專家的指導(dǎo)和深入討論。

    [1]殷躍平,張作辰,張開(kāi)軍.我國(guó)地面沉降現(xiàn)狀及防治對(duì)策研究[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2005,16(2):1-8.Yin Y P,Zhang Z C,Zhang K J.Land subsidence and countermeasures for its prevention in China[J].The Chinese Jounal of Geological Hazard and Control,2005,16(2):1- 8.

    [2]葛大慶,郭小方,張 玲,等.華北平原地面沉降區(qū)InSAR監(jiān)測(cè)[R].北京:中國(guó)國(guó)土資源航空物探遙感中心,2010.Ge D Q,Guo X F,Zhang L,et al.Subsidence monitoring with SAR interferometry in the North China Plain[R].Beijing:China Aero Geophysical Surveying and Remote Sensing Center for Land and Resources,2010.

    [3]王 艷,葛大慶,郭小方,等.長(zhǎng)江三角洲地區(qū)地面沉降InSAR監(jiān)測(cè)[R].北京:中國(guó)國(guó)土資源航空物探遙感中心,2010.Wang Y,Ge D Q,Guo X F,et al.Subsidence monitoring with SAR interferometry in the Yangtze River Delta area[R].Beijing:China Aero Geophysical Surveying and Remote Sensing Center for Land and Resources,2010.

    [4]Ferretti A,Prati C,Rocca F.Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry[J].IEEE Trans Geosci Remote Sensing,2001,38(5):2202-2212.

    [5]Mora O,Mallorqui J J,Broquetas A.Linear and nonlinear terrain deformation maps from a reduced set of interferometric SAR images[J].IEEE Trans Geosci Remote Sensing,2003,41(10):2243-2253.

    [6]Charles W,Wegmuller U,Andreas W,et al.Interferometric point target analysis with JERS-1 L-band SAR data[C]//Proceedings of the IEEE Transactfions on Geoscience and Remote Sensing.Toulouse,F(xiàn)rance,2003:4359-4361.

    [7]Hoffmann J,Zebker H A,Galloway D L,et al.Seasonal subsidence and rebound in Las Vegas Valley,Nevada,observed by synethetic aperture Radar interferometry[J].Water Resources Research,2001,37(6):1551-1566.

    [8]王 艷,廖明生,李德仁,等.利用長(zhǎng)時(shí)間序列相干目標(biāo)獲取地面沉降場(chǎng)[J].地球物理學(xué)報(bào),2007,50(2):598-604.Wang Y,Liao M S,Li D R,et al.Subsidence velocity retrieveal from long term coherent targets in Radar interferometric stacks[J].Chinese Jounal Geophys,2007,50(2):598-604.

    [9]http://www.terrafirma.eu.com/documents.htm

    [10]王小剛,陳 松,王秀芹,等.德州市地面沉降成因及防治對(duì)策淺析[J].地質(zhì)災(zāi)害與環(huán)境保護(hù),2006,17(3):62-66.Wang X G,Chen S,Wang X Q,et al.A Discussion on the causes of the ground subsidence and its countermeasures in Dezhou City[J].Journal of Geological Hazards and Environment Preservation,2006,17(3):62-66.

    [11]楊麗芝.德州深層地下水位降落漏斗演變機(jī)制與可調(diào)控性研究[D].北京:中國(guó)地質(zhì)科學(xué)院,2009.Yang L Z.Form principle and controlling- adjusting research about deep ground water depression cone in Dezhou[D].Beijing:Chinese Academy of Geologieal Seienee,2009.

    [12]王勝嶺,宋 波,王德生,等.德州市臨盤(pán)采油區(qū)地面沉降監(jiān)測(cè)[J].山東國(guó)土資源,2009,25(1):25-28,32.Wang S L,Song B,Wang D S,et al.Land subsidence monitoring of Linpan oil mining area in Dezhou City[J].Land and Resources in Shandong Province,2009,25(1):25-28,32.

    [13]何國(guó)榮,馮在敏,馮克印,等.水準(zhǔn)測(cè)量網(wǎng)監(jiān)測(cè)在德州市地面沉降中的應(yīng)用[J].山東國(guó)土資源,2012,28(6):28-30.He G R,F(xiàn)eng Z M,F(xiàn)eng K Y,et al.Application of leveling surveying for subsidence monitoring in Dezhou City[J].Land and Resources in Shandong Province,2012,28(6):28-30.

    猜你喜歡
    德州深層基線
    德州大陸架石油工程技術(shù)有限公司
    鉆采工藝(2022年4期)2022-10-22 10:24:20
    適用于MAUV的變基線定位系統(tǒng)
    航天技術(shù)與甚長(zhǎng)基線陣的結(jié)合探索
    科學(xué)(2020年5期)2020-11-26 08:19:14
    德州魯源貨場(chǎng)信號(hào)聯(lián)鎖設(shè)備關(guān)鍵技術(shù)的應(yīng)用
    考慮各向異性滲流的重力壩深層抗滑穩(wěn)定分析
    SAM系統(tǒng)對(duì)TDCS數(shù)據(jù)的優(yōu)化處理與深層應(yīng)用
    一種改進(jìn)的干涉儀測(cè)向基線設(shè)計(jì)方法
    在德州,電力經(jīng)紀(jì)人幫你選電!
    能源(2016年3期)2016-12-01 05:10:48
    德州地區(qū)懸鈴木方翅網(wǎng)蝽的綜合防治措施
    對(duì)“醫(yī)患失去信任”的深層憂慮
    男人和女人高潮做爰伦理| 一边亲一边摸免费视频| 三级国产精品欧美在线观看| 男人添女人高潮全过程视频| 免费黄色在线免费观看| a级毛片免费高清观看在线播放| 国产精品国产三级专区第一集| 黄色配什么色好看| 国产亚洲av片在线观看秒播厂| 国产精品99久久99久久久不卡 | 97在线人人人人妻| 热99国产精品久久久久久7| 一边亲一边摸免费视频| 亚洲精品自拍成人| 搞女人的毛片| 乱码一卡2卡4卡精品| 免费黄色在线免费观看| 一二三四中文在线观看免费高清| 又粗又硬又长又爽又黄的视频| 五月伊人婷婷丁香| 欧美激情久久久久久爽电影| 男女那种视频在线观看| 成人鲁丝片一二三区免费| 免费观看性生交大片5| 久久人人爽人人爽人人片va| 免费看a级黄色片| 黄色日韩在线| 国产免费又黄又爽又色| 尤物成人国产欧美一区二区三区| 亚洲国产最新在线播放| av又黄又爽大尺度在线免费看| 久久99蜜桃精品久久| 亚洲成人av在线免费| 亚洲人成网站在线播| 黄色视频在线播放观看不卡| 女人十人毛片免费观看3o分钟| 中文精品一卡2卡3卡4更新| 黄片无遮挡物在线观看| 免费观看av网站的网址| 久久久精品免费免费高清| 极品少妇高潮喷水抽搐| 综合色丁香网| 精品国产乱码久久久久久小说| 日韩人妻高清精品专区| 精品酒店卫生间| 色5月婷婷丁香| 国产成人a区在线观看| 久久久久久久久久久丰满| 99热这里只有精品一区| 久久精品夜色国产| 午夜福利视频1000在线观看| 国产高清有码在线观看视频| 亚洲国产欧美人成| 九九久久精品国产亚洲av麻豆| 午夜福利在线在线| 在线天堂最新版资源| 久久人人爽av亚洲精品天堂 | 一级毛片aaaaaa免费看小| 久久精品国产鲁丝片午夜精品| 三级国产精品片| 一个人观看的视频www高清免费观看| 久久久成人免费电影| 91aial.com中文字幕在线观看| 亚洲欧美精品自产自拍| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久精品精品| 嫩草影院精品99| 狂野欧美激情性xxxx在线观看| av免费在线看不卡| 亚洲一区二区三区欧美精品 | 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品乱码久久久v下载方式| 成人鲁丝片一二三区免费| 美女高潮的动态| 午夜爱爱视频在线播放| 国产片特级美女逼逼视频| 大香蕉97超碰在线| 国模一区二区三区四区视频| 亚洲国产精品国产精品| 国产精品99久久99久久久不卡 | 欧美日本视频| 青春草视频在线免费观看| 国产大屁股一区二区在线视频| 亚洲精品456在线播放app| 黄色视频在线播放观看不卡| 久久国内精品自在自线图片| 国精品久久久久久国模美| 欧美性感艳星| 国国产精品蜜臀av免费| 美女国产视频在线观看| 黄片wwwwww| 欧美日韩在线观看h| 国产爱豆传媒在线观看| 欧美最新免费一区二区三区| 在线精品无人区一区二区三 | 亚洲av电影在线观看一区二区三区 | 高清av免费在线| 久久久成人免费电影| 少妇被粗大猛烈的视频| 久久99蜜桃精品久久| 婷婷色麻豆天堂久久| 又大又黄又爽视频免费| 国内精品美女久久久久久| 在线天堂最新版资源| 欧美激情久久久久久爽电影| 麻豆成人av视频| 精品久久久久久久人妻蜜臀av| 一级爰片在线观看| av线在线观看网站| 午夜精品国产一区二区电影 | 大又大粗又爽又黄少妇毛片口| 又爽又黄a免费视频| 中文在线观看免费www的网站| av在线app专区| 老女人水多毛片| a级一级毛片免费在线观看| 亚洲三级黄色毛片| 久久久久久久国产电影| 嘟嘟电影网在线观看| 一级av片app| 国产高清三级在线| 中国美白少妇内射xxxbb| 国产欧美日韩精品一区二区| 美女被艹到高潮喷水动态| 欧美潮喷喷水| 亚洲丝袜综合中文字幕| 国产精品国产av在线观看| kizo精华| 国产毛片在线视频| 国产成人一区二区在线| 18+在线观看网站| 91精品一卡2卡3卡4卡| 秋霞伦理黄片| 青春草视频在线免费观看| 亚洲在线观看片| 又黄又爽又刺激的免费视频.| 99热这里只有是精品50| 国产日韩欧美在线精品| 国产日韩欧美在线精品| av专区在线播放| 中国美白少妇内射xxxbb| 欧美人与善性xxx| 成人高潮视频无遮挡免费网站| 亚洲av二区三区四区| 国产高清不卡午夜福利| 成人亚洲欧美一区二区av| 国产乱人视频| 中文字幕av成人在线电影| 亚洲一区二区三区欧美精品 | 亚洲国产精品成人综合色| 久久99热这里只频精品6学生| 精品少妇久久久久久888优播| 久久综合国产亚洲精品| 久久人人爽av亚洲精品天堂 | 各种免费的搞黄视频| 日韩三级伦理在线观看| 欧美精品国产亚洲| 一级毛片久久久久久久久女| 丝袜美腿在线中文| 久久99热这里只频精品6学生| 观看美女的网站| 久久6这里有精品| 国模一区二区三区四区视频| 免费看av在线观看网站| 别揉我奶头 嗯啊视频| 新久久久久国产一级毛片| 久久精品国产a三级三级三级| 亚洲人成网站在线播| 欧美高清性xxxxhd video| 成人综合一区亚洲| 久久久久久久久久成人| 色网站视频免费| 日韩欧美 国产精品| 婷婷色综合大香蕉| 国产成人免费观看mmmm| 我的女老师完整版在线观看| av又黄又爽大尺度在线免费看| 国产亚洲91精品色在线| 国产一区亚洲一区在线观看| 五月伊人婷婷丁香| av在线观看视频网站免费| 日本wwww免费看| 欧美另类一区| 亚洲精品影视一区二区三区av| 午夜老司机福利剧场| 久久久久国产精品人妻一区二区| 欧美日韩亚洲高清精品| 午夜精品一区二区三区免费看| 国产毛片在线视频| 亚洲电影在线观看av| 久久久久久九九精品二区国产| 哪个播放器可以免费观看大片| 国产一区二区在线观看日韩| 美女国产视频在线观看| 大片电影免费在线观看免费| 久久久精品94久久精品| 久久精品久久精品一区二区三区| 一区二区三区免费毛片| 性色avwww在线观看| 国产美女午夜福利| 美女高潮的动态| 国产久久久一区二区三区| 乱码一卡2卡4卡精品| a级一级毛片免费在线观看| 国产毛片在线视频| 国产白丝娇喘喷水9色精品| 亚洲欧美一区二区三区黑人 | 国产探花极品一区二区| av福利片在线观看| 国产女主播在线喷水免费视频网站| 亚洲av中文av极速乱| 亚洲欧美日韩无卡精品| 精品一区二区免费观看| 午夜激情久久久久久久| av一本久久久久| 日本爱情动作片www.在线观看| 丰满少妇做爰视频| 精品人妻熟女av久视频| 一级av片app| 国产综合懂色| 少妇高潮的动态图| 在线观看一区二区三区| 亚洲av一区综合| 一区二区三区免费毛片| 国产人妻一区二区三区在| 超碰97精品在线观看| 尤物成人国产欧美一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲,欧美,日韩| 国产高清不卡午夜福利| 亚洲精品自拍成人| 亚洲怡红院男人天堂| 日本猛色少妇xxxxx猛交久久| 嫩草影院新地址| 亚洲四区av| 卡戴珊不雅视频在线播放| 最后的刺客免费高清国语| 深爱激情五月婷婷| av在线播放精品| 国产淫片久久久久久久久| 亚洲精品国产av蜜桃| 男人添女人高潮全过程视频| 久久久久国产网址| 少妇人妻一区二区三区视频| 国产黄a三级三级三级人| 亚洲,一卡二卡三卡| 丝袜脚勾引网站| 亚洲人成网站在线播| 97超碰精品成人国产| av网站免费在线观看视频| 精品人妻视频免费看| 国产欧美另类精品又又久久亚洲欧美| 综合色av麻豆| 一区二区三区乱码不卡18| 久久久久性生活片| 欧美高清性xxxxhd video| 99久久精品一区二区三区| av福利片在线观看| 交换朋友夫妻互换小说| 精品国产三级普通话版| 成人漫画全彩无遮挡| 亚洲四区av| 国产久久久一区二区三区| 秋霞伦理黄片| 波野结衣二区三区在线| 亚洲精品,欧美精品| xxx大片免费视频| 精品久久久久久久末码| 一级黄片播放器| 男女那种视频在线观看| 国产精品久久久久久久电影| 看非洲黑人一级黄片| 男女国产视频网站| 久久国内精品自在自线图片| 久久精品国产鲁丝片午夜精品| 不卡视频在线观看欧美| 97热精品久久久久久| 国产精品一区www在线观看| 亚洲av在线观看美女高潮| 欧美+日韩+精品| 成人毛片60女人毛片免费| 亚洲国产av新网站| 日韩强制内射视频| 高清欧美精品videossex| 中文字幕亚洲精品专区| 91狼人影院| 国产免费一级a男人的天堂| 26uuu在线亚洲综合色| 听说在线观看完整版免费高清| 国产老妇伦熟女老妇高清| 成人亚洲精品av一区二区| 91精品一卡2卡3卡4卡| 日韩欧美 国产精品| 国内揄拍国产精品人妻在线| 亚洲最大成人av| 久久久久精品性色| 别揉我奶头 嗯啊视频| 日韩国内少妇激情av| 97人妻精品一区二区三区麻豆| 一边亲一边摸免费视频| 亚洲国产高清在线一区二区三| av在线蜜桃| 亚洲av男天堂| 在线观看美女被高潮喷水网站| 国产精品99久久99久久久不卡 | av国产免费在线观看| 激情 狠狠 欧美| 大香蕉97超碰在线| 国产在线一区二区三区精| 精品少妇黑人巨大在线播放| 亚洲精品aⅴ在线观看| 18禁裸乳无遮挡动漫免费视频 | 小蜜桃在线观看免费完整版高清| 国产高潮美女av| 黑人高潮一二区| 国模一区二区三区四区视频| 午夜精品国产一区二区电影 | 国产免费福利视频在线观看| 精品一区二区免费观看| 91精品伊人久久大香线蕉| 91精品国产九色| 欧美高清成人免费视频www| 欧美最新免费一区二区三区| 久久精品久久久久久噜噜老黄| 成人特级av手机在线观看| 日本欧美国产在线视频| 亚洲在线观看片| 亚洲欧美精品自产自拍| 亚洲精品乱久久久久久| 国产成人a∨麻豆精品| 中文字幕人妻熟人妻熟丝袜美| 欧美精品一区二区大全| 激情五月婷婷亚洲| 欧美极品一区二区三区四区| 日韩国内少妇激情av| 亚洲自拍偷在线| 深爱激情五月婷婷| 性色av一级| 欧美人与善性xxx| 日韩制服骚丝袜av| 性色avwww在线观看| 日韩电影二区| 婷婷色麻豆天堂久久| 美女xxoo啪啪120秒动态图| 亚洲欧美一区二区三区国产| 91狼人影院| 精品久久久久久久久亚洲| 日日摸夜夜添夜夜爱| 观看免费一级毛片| 一级毛片黄色毛片免费观看视频| 丰满人妻一区二区三区视频av| 大香蕉97超碰在线| 国产免费又黄又爽又色| 纵有疾风起免费观看全集完整版| 欧美少妇被猛烈插入视频| 麻豆乱淫一区二区| 搞女人的毛片| 神马国产精品三级电影在线观看| 一边亲一边摸免费视频| 亚洲四区av| 精品一区二区三卡| 九九在线视频观看精品| 美女主播在线视频| 亚洲av一区综合| 午夜精品国产一区二区电影 | 欧美精品人与动牲交sv欧美| av在线天堂中文字幕| 秋霞在线观看毛片| 成人高潮视频无遮挡免费网站| 亚洲av成人精品一二三区| 最后的刺客免费高清国语| 亚洲在线观看片| 亚洲精品久久久久久婷婷小说| 欧美zozozo另类| 男女那种视频在线观看| 亚洲一区二区三区欧美精品 | 欧美一级a爱片免费观看看| 亚洲电影在线观看av| 精品午夜福利在线看| 日韩亚洲欧美综合| 男人添女人高潮全过程视频| 欧美精品一区二区大全| 日韩强制内射视频| 少妇熟女欧美另类| 精品少妇黑人巨大在线播放| 亚洲自拍偷在线| 亚洲欧美日韩卡通动漫| 精品99又大又爽又粗少妇毛片| 国产精品.久久久| 久久久成人免费电影| 婷婷色综合www| 精品久久久噜噜| 欧美精品人与动牲交sv欧美| 成人特级av手机在线观看| 日韩一区二区三区影片| 熟女人妻精品中文字幕| 可以在线观看毛片的网站| 午夜日本视频在线| 天美传媒精品一区二区| 亚洲久久久久久中文字幕| 在线a可以看的网站| 亚州av有码| 国产精品99久久99久久久不卡 | 狂野欧美激情性xxxx在线观看| 国产视频内射| 岛国毛片在线播放| 少妇人妻久久综合中文| 国产亚洲av嫩草精品影院| 秋霞伦理黄片| 一区二区三区精品91| 亚洲色图综合在线观看| 亚洲美女搞黄在线观看| 国产av码专区亚洲av| 亚洲一级一片aⅴ在线观看| 日韩亚洲欧美综合| 日韩在线高清观看一区二区三区| 只有这里有精品99| 国产成人免费无遮挡视频| 极品少妇高潮喷水抽搐| 日本黄大片高清| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 国产毛片a区久久久久| av播播在线观看一区| 久久韩国三级中文字幕| 最近最新中文字幕大全电影3| 91久久精品电影网| 激情 狠狠 欧美| 精品99又大又爽又粗少妇毛片| 在线天堂最新版资源| 亚洲国产欧美人成| 在线免费观看不下载黄p国产| 亚洲国产精品999| 婷婷色综合大香蕉| 亚洲综合精品二区| 小蜜桃在线观看免费完整版高清| 成人黄色视频免费在线看| 一级毛片我不卡| 日韩在线高清观看一区二区三区| 97超视频在线观看视频| 成年免费大片在线观看| 全区人妻精品视频| 欧美一区二区亚洲| 真实男女啪啪啪动态图| 身体一侧抽搐| 69av精品久久久久久| 欧美成人午夜免费资源| 成人毛片60女人毛片免费| 午夜免费鲁丝| 日本一本二区三区精品| 街头女战士在线观看网站| 亚洲精品亚洲一区二区| 两个人的视频大全免费| 成人漫画全彩无遮挡| 中文在线观看免费www的网站| 我的老师免费观看完整版| 欧美性猛交╳xxx乱大交人| 极品教师在线视频| 高清午夜精品一区二区三区| 日韩av免费高清视频| 久久精品久久久久久久性| 久久久久久久久久久免费av| 免费电影在线观看免费观看| 啦啦啦啦在线视频资源| 亚洲电影在线观看av| 听说在线观看完整版免费高清| 观看免费一级毛片| 欧美日韩视频精品一区| 日本欧美国产在线视频| 日本熟妇午夜| 中文精品一卡2卡3卡4更新| 蜜桃亚洲精品一区二区三区| 纵有疾风起免费观看全集完整版| 好男人视频免费观看在线| 另类亚洲欧美激情| 我的女老师完整版在线观看| 三级国产精品片| 国产黄片美女视频| 91久久精品国产一区二区三区| 日日啪夜夜爽| 国产精品成人在线| 美女cb高潮喷水在线观看| 久久这里有精品视频免费| www.色视频.com| 免费看av在线观看网站| 你懂的网址亚洲精品在线观看| 欧美少妇被猛烈插入视频| 特大巨黑吊av在线直播| 欧美精品一区二区大全| 欧美成人a在线观看| 久久久久久久久久成人| 大码成人一级视频| 亚洲美女搞黄在线观看| 亚洲成人精品中文字幕电影| 国产美女午夜福利| 少妇高潮的动态图| 国产成人福利小说| 亚洲精品乱码久久久v下载方式| 国产极品天堂在线| 91久久精品国产一区二区成人| 青春草视频在线免费观看| 亚洲国产精品专区欧美| 中文字幕制服av| 神马国产精品三级电影在线观看| av女优亚洲男人天堂| 国产老妇女一区| 久久97久久精品| 欧美潮喷喷水| 中文精品一卡2卡3卡4更新| 亚洲精品乱码久久久v下载方式| 亚洲av中文av极速乱| 丰满少妇做爰视频| 人妻制服诱惑在线中文字幕| 丰满乱子伦码专区| 久久久a久久爽久久v久久| 精品人妻视频免费看| 日本一二三区视频观看| 亚洲天堂av无毛| 国产免费视频播放在线视频| 久久久精品免费免费高清| 最近最新中文字幕免费大全7| 亚洲精品国产av蜜桃| av免费观看日本| 少妇 在线观看| 制服丝袜香蕉在线| av卡一久久| 亚洲性久久影院| 18禁裸乳无遮挡动漫免费视频 | 久久久久久久久久成人| av一本久久久久| 国产探花在线观看一区二区| 菩萨蛮人人尽说江南好唐韦庄| 青春草国产在线视频| 永久免费av网站大全| av卡一久久| 国产成人freesex在线| 亚洲精品国产成人久久av| 男人添女人高潮全过程视频| 99久久中文字幕三级久久日本| 精品99又大又爽又粗少妇毛片| 亚洲自拍偷在线| 久久久久国产精品人妻一区二区| 国产有黄有色有爽视频| 69av精品久久久久久| 欧美日本视频| 国产亚洲午夜精品一区二区久久 | 亚洲激情五月婷婷啪啪| 男插女下体视频免费在线播放| 人人妻人人澡人人爽人人夜夜| 一区二区三区四区激情视频| 特级一级黄色大片| 嫩草影院精品99| 国产又色又爽无遮挡免| 久久久久久国产a免费观看| 一级爰片在线观看| 最近最新中文字幕免费大全7| 亚洲图色成人| 免费播放大片免费观看视频在线观看| 男人狂女人下面高潮的视频| 国产黄频视频在线观看| 最近2019中文字幕mv第一页| 在线天堂最新版资源| 91精品伊人久久大香线蕉| 免费看光身美女| 日本猛色少妇xxxxx猛交久久| 色视频www国产| 永久网站在线| 一区二区三区乱码不卡18| 国产精品人妻久久久久久| 亚州av有码| 国产色婷婷99| xxx大片免费视频| 日本一本二区三区精品| 亚洲精品乱久久久久久| av黄色大香蕉| 久久99热这里只有精品18| 日本三级黄在线观看| 国内揄拍国产精品人妻在线| 中文字幕亚洲精品专区| 国产一区二区在线观看日韩| 免费看日本二区| 黄片无遮挡物在线观看| 日韩一区二区三区影片| 亚洲av国产av综合av卡| 国产成人aa在线观看| 国内少妇人妻偷人精品xxx网站| 热99国产精品久久久久久7| 亚洲欧美一区二区三区黑人 | 麻豆成人午夜福利视频| 国产av码专区亚洲av| 免费观看的影片在线观看| 午夜福利在线观看免费完整高清在| 99久久精品国产国产毛片| 免费观看av网站的网址| 日本一二三区视频观看| 国产午夜福利久久久久久| 亚洲av免费高清在线观看| 精品国产露脸久久av麻豆| 欧美成人一区二区免费高清观看| 黄色日韩在线| 成人无遮挡网站| 国产有黄有色有爽视频| 婷婷色综合www| 三级国产精品欧美在线观看| 成年人午夜在线观看视频| 边亲边吃奶的免费视频| 久久精品夜色国产| 日本色播在线视频| 97在线人人人人妻| 亚洲伊人久久精品综合| 尤物成人国产欧美一区二区三区| av免费在线看不卡| 涩涩av久久男人的天堂| 日韩三级伦理在线观看| 成人一区二区视频在线观看|