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

    顧及抗差性的GNSS多尺度地殼運動場構(gòu)建方法

    2022-08-31 13:06:54蘇小寧朱慶孟國杰孫海麗閆浩文
    地球物理學報 2022年9期
    關(guān)鍵詞:模型

    蘇小寧,朱慶,孟國杰,孫海麗,閆浩文

    1 蘭州交通大學測繪與地理信息學院,蘭州 730070 2 西南交通大學地球科學與環(huán)境工程學院,成都 611756 3 中國地震局地震預(yù)測研究所,北京 100036 4 首都師范大學資源環(huán)境與旅游學院,北京 100048

    0 引言

    GNSS(Global Navigation Satellite System)觀測技術(shù)在研究全球地殼運動和構(gòu)造活動等方面發(fā)揮了重要的作用,“中國地殼運動觀測網(wǎng)絡(luò)”、“中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)”、首都圈GNSS觀測網(wǎng)和省屬地震氣象測繪部門的CORS(Continuously Operating Reference Stations)網(wǎng)等覆蓋全國和區(qū)域的多源觀測數(shù)據(jù),為獲取中國大陸較高空間分辨率的變形場提供了豐富的資料(Wang and Shen,2020).由于GNSS速度場相對于某個參考基準,因此選擇不同參考基準導(dǎo)致速度場的表現(xiàn)形式差異巨大.而應(yīng)變率場與參考基準選擇無關(guān),通過對比應(yīng)變積累和地震矩釋放的關(guān)系,在地震危險性分析和長期預(yù)報方面發(fā)揮著重要作用(Ojo et al., 2021),這其中的關(guān)鍵是解算出能反映測站空間不規(guī)則分布特征且精準的應(yīng)變率場.

    目前解算應(yīng)變場的方法可總結(jié)為兩類:一類是建立塊體-斷層模型,計算塊體內(nèi)部均勻應(yīng)變,稱為物理模型(申重陽等,2002;Tape et al., 2009);另一種是利用GNSS速度場,構(gòu)建空間連續(xù)的速度場模型,進而解算應(yīng)變率場,稱為數(shù)學模型(Tape et al., 2009).可見,在數(shù)學模型中關(guān)鍵是如何精準構(gòu)建空間連續(xù)的速度場模型,常用的方法有距離加權(quán)(Shen et al,1996)、多面函數(shù)(武艷強等,2009)、Kriging插值(劉曉霞等,2014)、最小二乘配置(江在森和劉經(jīng)南,2010)和多尺度球面小波(Tape et al., 2009;程鵬飛等,2015;蘇小寧等,2016)等,這些方法在確定觀測值權(quán)重時或?qū)NSS速度場看作為等精度的觀測值,或利用擬合的速度場形式誤差表征其不等精度的特征.但實際上經(jīng)常存在由于觀測數(shù)據(jù)質(zhì)量較差導(dǎo)致的誤差較大或者土層觀測點標墩不穩(wěn)定產(chǎn)生的局部速度差異,也就是在觀測數(shù)據(jù)中存在少量粗差.此時利用現(xiàn)有方法構(gòu)建速度場模型時,少量粗差將會對解算的應(yīng)變場產(chǎn)生重要的影響,即上述方法均不具備抗差功能.

    本文基于能夠體現(xiàn)GNSS測站不規(guī)則分布特征的多尺度球面小波方法構(gòu)建GNSS速度場,獲取不同空間尺度的地殼運動場,在此基礎(chǔ)上通過抗差迭代算法顧及速度場中存在的少量粗差對建立速度場模型的影響,通過在實測數(shù)據(jù)中人為加入粗差的模擬計算驗證方法的抗差性,結(jié)合川滇地區(qū)最新的GNSS速度場數(shù)據(jù)給出了區(qū)域的速度場和應(yīng)變率場模型.

    1 基于球面小波算法的多尺度地殼運動場建模

    1.1 速度場函數(shù)模型

    球面上觀測點的三維速度可以以南北向、東西向和垂直向3個分量來表示,本文只考慮水平向速度場,因此只涉及到南北向和東西向2個分量,考慮垂直向時的方法與本文所述相同.地面觀測站(λ,φ)的地殼運動速度ν(λ,φ)可表示為

    (1)

    基于球面小波算法構(gòu)建多尺度地殼運動場模型的本質(zhì)是采用小波基函數(shù)作為估計速度場的框架,利用對小波母函數(shù)進行平移和伸縮等仿射變換后得到的小波基函數(shù)來表示不同分量的速度場,因此公式(1)可以變換為

    (2)

    (3)

    1.2 空間球面小波基函數(shù)

    以高斯函數(shù)差建立的球面小波基函數(shù)(Difference of Gaussians,DOG)在構(gòu)建多尺度地殼運動場模型時得到了廣泛應(yīng)用(蘇小寧等,2016),DOG小波的函數(shù)模型為(包伯成等,2009)

    (4)

    1.3 空間球面小波分解

    空間球面小波分解的原則是根據(jù)觀測站的分布建立一系列球面小波基函數(shù).通過對球面進行三角形分解,可以獲得在球面上均勻分布的球面三角形網(wǎng)格,分解的層次與球面小波基函數(shù)的尺度相對應(yīng),也對應(yīng)于不同邊長的球面三角形,具體形式可參考 Wang和Dahlen(1995).首先根據(jù)分解尺度建立空間規(guī)則分布的小波基函數(shù),再根據(jù)地表觀測點的空間密度對規(guī)則分布的小波基函數(shù)進行篩選,使得最終的小波基函數(shù)集在空間分布上與地表觀測站的空間密度密切相關(guān),進而構(gòu)建的模型能夠凸顯測站不規(guī)則分布的特征.分解的具體步驟如下:

    (1)選取研究區(qū)經(jīng)緯度的覆蓋范圍形成球面矩形,計算球面矩形對角線的長度,搜索小于該長度且距離最近的三角形邊長,確定最小分解尺度;

    (2)根據(jù)不同分解尺度,建立多層次且規(guī)則分布的空間小波基函數(shù),標注對應(yīng)的球面三角形邊長;

    (3)根據(jù)地表觀測點分布,對所有小波基函數(shù)逐個搜索,以小波基函數(shù)所在位置為圓心和分解尺度對應(yīng)的球面三角形邊長為半徑,確定小波基函數(shù)所對應(yīng)的“空間支撐”范圍,統(tǒng)計“空間支撐”內(nèi)地面觀測站的個數(shù),篩選出“空間支撐”內(nèi)部多于3個站點的小波基函數(shù),獲得模型可用的最大分解尺度、空間不均勻分布的小波基函數(shù)及其數(shù)量;

    (4)統(tǒng)計分析不同分解尺度小波基函數(shù)的有效覆蓋范圍,選擇有效覆蓋范圍超過70%的最大分解尺度為最終的球面小波分解尺度,形成在空間不規(guī)則分布且不同分解尺度的小波基函數(shù)集.

    1.4 觀測方程病態(tài)問題及其最優(yōu)化解

    根據(jù)地表觀測站的空間分布建立的球面小波基函數(shù)個數(shù)為N個,因此與其對應(yīng)的系數(shù),也就是待求參數(shù)的個數(shù)也是N個,而觀測點的個數(shù)為M個,也就是觀測量個數(shù)為M.在實際建立速度場模型時,觀測點個數(shù)往往小于待求參數(shù)個數(shù),導(dǎo)致觀測方程的參數(shù)求解成為對欠定方程的求解問題.同時,由母函數(shù)經(jīng)過仿射變換得到的DOG小波基函數(shù)為非正交基,導(dǎo)致法方程秩虧,涉及到不適定方程的求解問題.

    在建立速度場模型時,可以通過最優(yōu)化算法中的正則化方法,求解出待求參數(shù)向量,目標函數(shù)和參數(shù)的最優(yōu)解分別為

    (5)

    m=(GTPDG+λ2I)-1GTPDd,

    (6)

    式中S(m)為目標函數(shù),通過以目標函數(shù)為最小準則獲取參數(shù)的最優(yōu)解.G為由小波基函數(shù)組成的系數(shù)矩陣,m為待求參數(shù)向量,d為地表GNSS觀測點觀測的速度場向量,CD為觀測數(shù)據(jù)的方差-協(xié)方差矩陣,其逆矩陣為權(quán)矩陣PD,考慮到不同位置測站計算的地殼運動速度為獨立觀測量,則CD和PD均為對角矩陣.CM為待求參數(shù)的方差-協(xié)方差矩陣,本文采用Tikhonov正則化方法(徐天河和楊元喜,2003),CM的具體形式為λ2I,I為單位對角陣,λ為正則化參數(shù).利用公式(6),通過多次迭代之后,則可以求解出參數(shù)m的平差值.可以看出,由于對法方程矩陣的對角元素進行了修正,顯然該平差值為有偏估計值.

    從公式(6)可以看出,在參數(shù)最優(yōu)解求解的過程中,如果已知基函數(shù)分布和最優(yōu)化解算方法,則參數(shù)的解只依賴于如何確定觀測數(shù)據(jù)的權(quán)矩陣.以往普遍采用擬合的地表測站GNSS速度場的形式誤差給定,但如果在觀測數(shù)據(jù)中存在異常值或者粗差值,利用形式誤差確定權(quán)矩陣的方法并不能精確地表達其隨機特征,此時公式(6)不具備抗差功能,則無法獲得精確的參數(shù)解,最終建立的速度場模型值得商榷,因此,有必要對公式(6)進行修正,使其具備抗差功能.

    2 顧及抗差迭代算法的多尺度地殼運動場模型構(gòu)建

    一個好的估計方法不僅要具有“穩(wěn)健性”,還得具有“抗干擾性”.“穩(wěn)健性”是指當假定模型與實際模型有微小差異時,其估計方法的性能只受到微弱的影響;“抗干擾性”是指當觀測樣本中混入少量粗差時,估計值受到的影響不大(楊元喜等,2002).地表不規(guī)則分布的地殼運動速度場屬于獨立觀測,且不等精度,適合采用按等價權(quán)原理作抗差估計,目的是解決不等精度獨立觀測值權(quán)函數(shù)的確定問題(周江文,1989;楊元喜,1994).

    2.1 顧及抗差性的最優(yōu)化算法原理

    在公式(5)所示最優(yōu)化算法的目標函數(shù)中,第一項為觀測數(shù)據(jù)殘差的加權(quán)平方和,顯然個別異常大殘差的出現(xiàn)會導(dǎo)致平方和顯著增大.為了達到平方和最小的目標函數(shù),平差值必然會遷就那些異常觀測值或者粗差值,使得平差值偏離最優(yōu)值,因此個別異常值或粗差會對整個估值產(chǎn)生較大的影響.抗差估計的本質(zhì)是利用增長較慢的函數(shù)代替平方和函數(shù),得到具有較好抗粗差性的估值.

    (7)

    2.2 顧及抗差性的迭代最優(yōu)化算法實現(xiàn)步驟

    (1)確定待估參數(shù)的初值,確定權(quán)矩陣的初值,給出觀測數(shù)據(jù)的方差-協(xié)方差矩陣;

    (2)按照公式(5)和公式(6)所示的傳統(tǒng)最優(yōu)化算法估計參數(shù)的平差值,計算觀測數(shù)據(jù)與模型值的殘差;

    (5)對參數(shù)平差值與初值的差值進行評價,如果其符合限差要求,則結(jié)束求解;否則重復(fù)步驟(2)—(4)進行迭代求解,直至待求參數(shù)平差值的差值符合要求為止.其中第k+1次的迭代解為

    (8)

    3 川滇地區(qū)地殼運動場模型

    3.1 多源數(shù)據(jù)融合的川滇地區(qū)GNSS震間速度場

    中國大陸觀測地殼運動的GNSS測站主要來自于中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò),該網(wǎng)絡(luò)包括260個連續(xù)觀測站和2000個流動觀測站,其中32個連續(xù)觀測站和1056個流動觀測站觀測起始于1998年,其余測站觀測起始于2009年,所有的數(shù)據(jù)均已經(jīng)獲得了超過10年的觀測資料.除此之外,在一些地震活動性較強的區(qū)域,例如首都圈地區(qū)、山西地塹帶和川滇地區(qū)等,也進行了GNSS流動觀測網(wǎng)的加密觀測;省屬的地震、氣象和測繪部門為了滿足區(qū)域性CORS網(wǎng)的需求,在所屬范圍也布設(shè)了連續(xù)觀測站,這些測站中的一部分能夠滿足地殼運動觀測的需求.Wang和Shen(2020)對以上數(shù)據(jù)進行了收集和整理,同時也收集了中國大陸周邊的GNSS觀測資料,通過高精度數(shù)據(jù)處理與震間形變場提取,給出了中國大陸高空間分辨率的GNSS震間速度場.數(shù)據(jù)處理采用GAMIT/GLOBK軟件(Herring et al., 2015, 2018),所有數(shù)據(jù)解算采用統(tǒng)一的策略,確保解算的坐標參數(shù)的高精度和強自洽性.在GNSS原始觀測數(shù)據(jù)解算階段,將地面測站坐標、地球定向參數(shù)(極移和協(xié)調(diào)時間時)和衛(wèi)星軌道一并估算,有效地消除了由于地球物理模型或者軌道參數(shù)變化所帶來的系統(tǒng)性偏差(Wang and Shen, 2020).在時間序列建模階段,同時考慮了大地震同震和震后位移的影響;考慮到不同天線類型和可能的天線設(shè)置誤差,對GAMIT給出的測站坐標誤差進行了重新分配,并利用有色噪聲隨機模型構(gòu)造觀測數(shù)據(jù)的方差-協(xié)方差矩陣,體現(xiàn)出時間序列中的時間相關(guān)性特征(Tao et al., 2021).

    在Wang和Shen(2020)給定的GNSS速度場中,中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)的數(shù)據(jù)截止于2015年,本文采用與其一致的觀測數(shù)據(jù)解算策略、時間序列建模策略和地殼運動解算策略,補充處理了該觀測網(wǎng)2017、2018、2019和2020年的觀測數(shù)據(jù),增加了觀測數(shù)據(jù)的長度,也提高了速度場的穩(wěn)定性,更新的川滇地區(qū)GNSS速度場如圖1所示.

    圖1 川滇地區(qū)相對于歐亞參考框架的GNSS速度場,誤差橢圓置信區(qū)間為95%Fig.1 GNSS velocity field with respect to Eurasian reference frame in the Sichuan and Yunnan regions. The confidence interval of error ellipse is 95%

    圖2 川滇地區(qū)GNSS測站分布的球面小波分解尺度Fig.2 Decomposition scale of spherical wavelet based on the distribution of GNSS stations in the Sichuan and Yunnan regions

    3.2 基于球面小波算法的多尺度地殼運動場模型

    根據(jù)GNSS測站分布,依據(jù)空間球面小波分解原則,獲得川滇地區(qū)的球面小波分解尺度,如圖2所示.研究區(qū)內(nèi)最大分解尺度的最小值為6,最大分解尺度的最大值為9.除西南角和東南角的邊緣區(qū)域,其他區(qū)域的最大分解尺度均在7以上.最大分解尺度能達到9的區(qū)域主要沿著主要構(gòu)造帶分布,例如小江斷裂帶南端、紅河斷裂與麗江—小金河斷裂交匯的區(qū)域和鮮水河斷裂與龍門山斷裂南端交匯的三岔口區(qū)域.川滇菱形塊體覆蓋的主要區(qū)域最大分解尺度都能達到8,因此本文最終選擇的最大分解尺度為8,分解尺度的選擇范圍為3~8.

    正則化參數(shù)的求解采用留一交叉驗證法(Leave-one-out Cross Validation,LOOCV),首先去除一個觀測點數(shù)據(jù),計算待求參數(shù),然后計算該點觀測值與模型估值之間的差值,通過全局搜索,選擇出最優(yōu)的正則化參數(shù),進而求解出待求參數(shù)的最優(yōu)解,正則化參數(shù)的搜索過程如圖3所示.選擇1~1010作為搜索范圍,可以看出該過程搜索到了函數(shù)的全局最小值,進而提取出正則化參數(shù)λ的取值,圖3中橫坐標取正則化參數(shù)λ的常用對數(shù),其在東西分量和南北分量分別為533和849.

    圖3 基于LOOCV算法的最優(yōu)化算法正則化參數(shù)Fig.3 Selection of regulation parameters based on optimized LOOCV method

    確定球面小波分解尺度和正則化參數(shù)之后,通過多次迭代即可求解出待估參數(shù)的最優(yōu)化值,進而正演出研究范圍內(nèi)空間任意一點的地殼運動速度,實現(xiàn)地殼形變場建模,給出空間連續(xù)的地殼運動圖像,如圖4所示,圖中冷色(藍色)代表小速率值,暖色(紅色)代表大速率值.可以看出,安寧河—則木河—小江斷裂組成的構(gòu)造帶是川滇菱形塊體地殼運動速率的分界帶.

    圖4 川滇地區(qū)地殼運動場模型(a) 東西分量; (b) 南北分量.Fig.4 Crustal deformation model in the Sichuan and Yunnan regions(a) The EW component; (b) The NS component.

    3.3 顧及抗差迭代算法的多尺度地殼運動場模型

    抗差迭代算法可有效削弱異常觀測值和粗差對地殼形變場建模的影響,進而獲得準確的應(yīng)變率場.通過多次抗差迭代計算,在所有參與建模的367個GNSS測站中,檢測出74個降權(quán)觀測值,占比為20%,檢測出7個除權(quán)觀測值,占比為2%,圖5中藍色三角形站點為正常觀測站點,黑色三角形站點為降權(quán)觀測站點,紅色三角形站點為除權(quán)觀測站點.

    圖5 抗差迭代算法給出的川滇地區(qū)三類GNSS測站分布圖Fig.5 The distribution of three types of GNSS stations in Sichuan and Yunnan regions based on robust iterative algorithm

    圖6給出了利用抗差迭代算法對部分測站進行降權(quán)和除權(quán)之后,基于多尺度球面小波算法構(gòu)建的川滇地區(qū)地殼運動場模型,可以看出考慮抗差性前后計算的速度場空間分布形態(tài)基本一致,安寧河—則木河—小江斷裂是川滇菱形塊體的東邊界這一基本特征未發(fā)生明顯改變.為了進一步分析抗差迭代算法在構(gòu)建地殼形變場時的作用,對本文顧及抗差性的方法與傳統(tǒng)多尺度球面小波方法獲得的地殼運動場模型進行對比,結(jié)果如圖7所示,圖中也給出了降權(quán)站點和除權(quán)站點的分布.可以看出,二者差異的量值在1.0 mm·a-1以內(nèi),說明雖然在構(gòu)建速度場模型時考慮了抗差性,但由于原有速度場本身數(shù)據(jù)質(zhì)量較好,因此得到的結(jié)果是穩(wěn)定的,二者的差異較小,而差異性較大的區(qū)域總是覆蓋著被降權(quán)或者除權(quán)的測站.

    4 討論

    4.1 GNSS速度場中存在粗差測站的原因

    利用目視法對基于GNSS觀測資料解算的中國大陸地殼運動場進行分析,發(fā)現(xiàn)不可避免地存在一些與周圍測站的運動特征明顯不協(xié)調(diào)的測站.其產(chǎn)生的原因主要包括:(1)因觀測時間跨度太短而無法擬合出真正的震間運動;(2)觀測標墩建于土層之上,由于季節(jié)變化等因素導(dǎo)致標墩本身存在非線性的變化;(3)測站觀測環(huán)境的改變導(dǎo)致觀測噪聲增加或運動特征偏離真值.

    本文所采用的觀測數(shù)據(jù)均超過了10年,超過解算震間變形時對時間跨度4.5年的要求(Blewitt et al., 2013).觀測環(huán)境的改變在測站的日志文件中都會有明確記載,在數(shù)據(jù)預(yù)處理階段已經(jīng)被考慮到.因此,為了分析測站被降權(quán)處理和被除權(quán)處理的原因,本文重點對這些測站的觀測標墩和數(shù)據(jù)質(zhì)量進行了進一步分析.其中測站標墩分為基巖點和土層點,通常認為基巖點能夠精準代表地殼的運動特征,而土層點會受到土壤隨著季節(jié)變化的熱脹冷縮和土層本身不穩(wěn)定等因素的影響(譚偉杰等,2017),產(chǎn)生局部效應(yīng)而可能無法精準代表地殼的運動特征.被降權(quán)和除權(quán)處理的81個站點全部為土層站點,可見觀測結(jié)果中存在明顯的局部影響.通過觀測數(shù)據(jù)形式誤差的大小體現(xiàn)出數(shù)據(jù)質(zhì)量的差異性,分析權(quán)不變、降權(quán)和除權(quán)三類測站形式誤差的差異性.計算三類數(shù)據(jù)東西向和南北向形式誤差的均值,權(quán)不變、降權(quán)和除權(quán)三類數(shù)據(jù)東西向形式誤差的均值分別為0.61、0.70和0.93 mm·a-1,南北向形式誤差的均值分別為0.52、0.60和0.83 mm·a-1,可以看出三類數(shù)據(jù)的形式誤差存在差異,表明其時間序列的離散度逐漸增大,數(shù)據(jù)質(zhì)量逐漸降低.

    圖6 顧及抗差迭代算法的川滇地區(qū)地殼運動場模型(a) 東西分量; (b) 南北分量.Fig.6 Crustal deformation model in Sichuan and Yunnan regions considering robust iterative algorithm(a) The EW component; (b) The NS component.

    圖7 顧及抗差迭代算法與傳統(tǒng)方法構(gòu)建的地殼運動場在東西分量的差值Fig.7 The difference of crustal deformation model in the EW component between our and traditional method

    4.2 抗差迭代算法的抗差性分析

    為了進一步檢驗抗差迭代算法的抗差性能,在速度場空間變化較大的川滇地區(qū)人為地引入粗差,分析粗差引入前后速度場模型的差異性特征.引入粗差的方式有在單個測站引入粗差和在相鄰的幾個測站引入粗差兩種,具體如圖8所示,其中測站H335、H387和H380為3個在區(qū)域上獨立的含有粗差的測站,測站H346、H111、H112、SCYY和H103為在區(qū)域上具有相關(guān)性的5個含有粗差的測站.用含有粗差的測站構(gòu)建運動場模型,并與不含粗差的數(shù)據(jù)構(gòu)建的結(jié)果進行對比,分析本文算法的抗差性.

    圖8 速度場中單個測站和區(qū)域測站加入人為誤差的兩種方案Fig.8 Two schemes of adding gross errors artificially of individual station and regional stations in velocity field

    圖9 加入人為粗差前后東西向速度場模型的差異分布Fig.9 The difference of crustal deformation model in EW component before and after adding gross errors artificially

    圖9給出了對二者速度場求差的結(jié)果,圖中-0.5~0.5 mm·a-1以內(nèi)的值用同一個色標表示,可以看出,其中測站H335和H387兩個單獨的存在粗差的測站經(jīng)過抗差迭代算法之后,對建模的結(jié)果并未產(chǎn)生影響,而另一個獨立的H380測站的粗差卻對建模結(jié)果產(chǎn)生了明顯的影響,進一步對測站分布分析可以看出,在H335和H387兩個測站10 km以內(nèi)的范圍還有其他測站,而H380測站在10 km以內(nèi)的范圍內(nèi)沒有其他測站,也就是說抗差迭代算法從測站H335和H387檢測出了粗差,而從H380測站沒有檢測到粗差,把含有粗差的觀測值當作好的觀測值對待,這是由于該測站周邊沒有其他觀測站作為對比,即便其本身含有粗差也被當做了有效值.由測站SCYY、H346、H103、H111和H112等5個測站組成的含有粗差的區(qū)域?qū)=Y(jié)果也產(chǎn)生了明顯的影響,表明區(qū)域性存在的粗差會對建模結(jié)果產(chǎn)生影響.

    4.3 抗差迭代算法對解算應(yīng)變率場的貢獻

    從圖10可以看出,川滇地區(qū)重要應(yīng)變積累的區(qū)域沿著主要的活動構(gòu)造帶分布,其中主要的一個條帶是安寧河—則木河—小江斷裂帶,這也是川滇菱形塊體的東邊界,紅河斷裂與麗江—小金河斷裂交匯的區(qū)域是另外一個重要的應(yīng)變積累區(qū)域,該分布特征與已有的研究成果保持一致(Qu et al., 2018),但在應(yīng)變集中區(qū)的量值上普遍高于已發(fā)表的結(jié)果(Pan and Shen, 2017;李長軍等,2019),這主要是由于本文采用的多尺度球面小波的應(yīng)變計算方法能夠凸顯出測站不均勻分布的特征,獲取不同空間分辨率的應(yīng)變率場,進而提取出應(yīng)變積累的集中區(qū)域.圖10也給出本文方法與傳統(tǒng)多尺度球面小波方法解算的應(yīng)變率場的差異,從第二不變量的量值來看,應(yīng)變差異值的量值遠低于應(yīng)變的量值,表明兩個結(jié)果均是穩(wěn)健的,對比圖10和圖7可以看出,應(yīng)變率差異較大的區(qū)域也是速度場差異較大的區(qū)域,同時也是存在被降權(quán)或除權(quán)測站的區(qū)域.

    4.4 川滇地區(qū)多尺度應(yīng)變積累特征

    基于多尺度球面小波算法解算應(yīng)變率場的優(yōu)點之一是能夠獲取不同空間尺度的應(yīng)變率場,從而體現(xiàn)出不同規(guī)模構(gòu)造運動的應(yīng)變積累特征.小的分解尺度能夠獲取大規(guī)模構(gòu)造運動的應(yīng)變積累,例如塊體邊界的應(yīng)變積累;大的分解尺度能夠獲取小規(guī)模構(gòu)造運動的應(yīng)變積累,例如塊體邊界應(yīng)變積累的差異性分布或塊體內(nèi)部次級斷層的應(yīng)變積累.圖11分別給出了分解尺度為3~7和8的應(yīng)變率第二不變量,二者之和為區(qū)域的總應(yīng)變,總體特征與已有結(jié)果保持一致(Jin et al, 2019; Zhang et al., 2019).從分解尺度為3~7的應(yīng)變率場可以看出,應(yīng)變積累的主要區(qū)域為塊體運動的邊界區(qū)域,作為川滇菱形塊體東邊界的安寧河—則木河—小江斷裂帶,其是川滇地區(qū)主要的應(yīng)變積累區(qū)域,其中以安寧河斷裂和則木河斷裂為東邊界的大涼山次級塊體內(nèi)部存在應(yīng)變積累,沿著小江斷裂的應(yīng)變積累穿過紅河斷裂之后一直向西南方向延伸.其他主要的應(yīng)變積累區(qū)域是位于紅河斷裂與麗江—小金河交匯的區(qū)域和中緬邊界的幾條西南向斷裂分布區(qū)域,主要包括畹町斷裂和南汀河斷裂等.分解尺度為8的應(yīng)變積累分布顯示區(qū)域內(nèi)除了沿著主要斷裂的應(yīng)變積累外,還存在局部小區(qū)域的應(yīng)變積累,這些應(yīng)變積累主要體現(xiàn)了沿著斷裂應(yīng)變積累的差異性特征.

    圖10 川滇地區(qū)應(yīng)變率第二不變量(a) 本文方法解算的應(yīng)變率場; (b) 本文方法與傳統(tǒng)多尺度球面小波方法解算結(jié)果的差異.Fig.10 The second invariant of strain rates in Sichuan and Yunnan regions(a) The strain rate field calculated by our method; (b) The difference between our and traditional method.

    圖11 川滇地區(qū)應(yīng)變率第二不變量(a) 分解尺度為3~7; (b) 分解尺度為8.Fig.11 The second invariant of strain rates in Sichuan and Yunnan regions(a) The decomposition scale is from 3 to 7; (b) The decomposition scale is 8.

    5 結(jié)論

    利用長期積累的GNSS觀測資料經(jīng)過高精度統(tǒng)一數(shù)據(jù)處理提取地殼運動圖像,然后解算GNSS應(yīng)變率場,進而對應(yīng)變積累與地震矩釋放進行對比分析,這些是開展區(qū)域地震危險性分析的主要工作,并在實際地震潛勢研究中得到了廣泛應(yīng)用和驗證.由于GNSS時間序列中存在與長期構(gòu)造運動無關(guān)的信息致使擬合的速度場無法真實表達實際的震間地殼運動特征,進而影響到應(yīng)變率場的結(jié)果.

    本文提出了顧及抗差性的速度場模型構(gòu)建方法,考慮到了測站因局部運動和數(shù)據(jù)質(zhì)量較差在速度場中產(chǎn)生的粗差.以解算的川滇地區(qū)最新GNSS觀測資料為例,在觀測的速度場中檢測出了20%的降權(quán)測站和2%的除權(quán)測站,對這些測站進一步分析發(fā)現(xiàn),粗差產(chǎn)生的原因是由于觀測數(shù)據(jù)質(zhì)量較差和測站為土層觀測點而存在局部非構(gòu)造變形,進而獲得穩(wěn)健性較好的速度場模型.通過在觀測數(shù)據(jù)人為加入粗差檢驗本文方法的抗差性,結(jié)果表明單個存在的粗差對速度場模型不會產(chǎn)生影響,而區(qū)域性存在的粗差會對結(jié)果產(chǎn)生影響,表明本文的方法具有好的抗差性.本文進一步分析了顧及抗差性算法對解算應(yīng)變率場的貢獻,進而給出了川滇地區(qū)精準的應(yīng)變率場.

    致謝感謝中國地震臺網(wǎng)中心提供的中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)GNSS觀測數(shù)據(jù).

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    交换朋友夫妻互换小说| 日韩视频在线欧美| 99久久综合免费| 国产视频首页在线观看| 国产一区二区在线观看日韩| 午夜日本视频在线| 久久久久久九九精品二区国产| 在线天堂最新版资源| 精品少妇黑人巨大在线播放| 成人一区二区视频在线观看| 老熟女久久久| 国产无遮挡羞羞视频在线观看| 青春草国产在线视频| 一级黄片播放器| 丝袜脚勾引网站| 丝袜脚勾引网站| 色婷婷久久久亚洲欧美| 精品少妇久久久久久888优播| 国产精品三级大全| 久久久久视频综合| 久久久久久伊人网av| 一级毛片 在线播放| 交换朋友夫妻互换小说| av在线蜜桃| 欧美日韩一区二区视频在线观看视频在线| 久久久久久久久久久免费av| 女性生殖器流出的白浆| 精品亚洲成国产av| 老女人水多毛片| 亚洲av成人精品一区久久| 国产成人精品一,二区| 97在线人人人人妻| 中文字幕久久专区| 久久这里有精品视频免费| 国产精品国产av在线观看| 国产伦精品一区二区三区四那| av专区在线播放| 日韩一区二区视频免费看| 成年女人在线观看亚洲视频| 日日啪夜夜爽| 高清在线视频一区二区三区| 国产日韩欧美在线精品| 亚洲三级黄色毛片| 国产精品人妻久久久久久| 欧美激情国产日韩精品一区| 伦理电影大哥的女人| 美女脱内裤让男人舔精品视频| 看十八女毛片水多多多| 色综合色国产| 亚洲av二区三区四区| 国产精品一区二区在线观看99| 免费观看在线日韩| 国产高潮美女av| 国产亚洲5aaaaa淫片| 国产午夜精品久久久久久一区二区三区| 亚洲精品,欧美精品| 精华霜和精华液先用哪个| 亚洲国产精品专区欧美| 色吧在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 99国产精品免费福利视频| 国产在线免费精品| 色网站视频免费| 两个人的视频大全免费| 婷婷色麻豆天堂久久| 日本色播在线视频| 美女国产视频在线观看| 九色成人免费人妻av| 中国国产av一级| 嫩草影院新地址| 欧美xxxx黑人xx丫x性爽| 春色校园在线视频观看| 色婷婷久久久亚洲欧美| 国产精品伦人一区二区| 国产日韩欧美在线精品| 国产伦理片在线播放av一区| av天堂中文字幕网| 成年美女黄网站色视频大全免费 | 男男h啪啪无遮挡| 青青草视频在线视频观看| 91在线精品国自产拍蜜月| 春色校园在线视频观看| 精品人妻熟女av久视频| 日本黄色片子视频| 欧美成人精品欧美一级黄| 国产久久久一区二区三区| 国产高清三级在线| 国产色婷婷99| 精品一品国产午夜福利视频| 纯流量卡能插随身wifi吗| 亚洲av不卡在线观看| 纯流量卡能插随身wifi吗| av黄色大香蕉| 久热久热在线精品观看| 伊人久久国产一区二区| 观看免费一级毛片| 久久久欧美国产精品| 亚洲精品一二三| 九草在线视频观看| 肉色欧美久久久久久久蜜桃| 国产免费福利视频在线观看| 久久久久久久久久人人人人人人| 久久亚洲国产成人精品v| 亚洲最大成人中文| 亚洲av二区三区四区| 制服丝袜香蕉在线| 久久久午夜欧美精品| 精品人妻视频免费看| 看非洲黑人一级黄片| 亚洲av男天堂| 国产亚洲最大av| 亚洲aⅴ乱码一区二区在线播放| 国内精品宾馆在线| 国产视频内射| 国产在线一区二区三区精| 久久精品熟女亚洲av麻豆精品| 蜜桃在线观看..| 久久久久久久久久成人| 国产黄片美女视频| 欧美精品一区二区大全| 日韩一本色道免费dvd| 3wmmmm亚洲av在线观看| av又黄又爽大尺度在线免费看| a级毛片免费高清观看在线播放| 国产精品精品国产色婷婷| 国产精品久久久久久av不卡| 丝袜脚勾引网站| 国产精品久久久久久久久免| 成年免费大片在线观看| 纯流量卡能插随身wifi吗| 免费黄频网站在线观看国产| 少妇的逼水好多| 国产欧美亚洲国产| 最近中文字幕2019免费版| 亚洲av欧美aⅴ国产| 亚洲精品日韩av片在线观看| 日韩三级伦理在线观看| 国产国拍精品亚洲av在线观看| 免费观看性生交大片5| 五月天丁香电影| 男人狂女人下面高潮的视频| 尤物成人国产欧美一区二区三区| 久久毛片免费看一区二区三区| 国产日韩欧美在线精品| 国产欧美日韩一区二区三区在线 | 亚洲怡红院男人天堂| 国产精品免费大片| 免费观看无遮挡的男女| 男人和女人高潮做爰伦理| a级毛色黄片| 丝袜喷水一区| 人人妻人人爽人人添夜夜欢视频 | 国产色婷婷99| 十八禁网站网址无遮挡 | 国产爱豆传媒在线观看| 日韩中字成人| 精品午夜福利在线看| 校园人妻丝袜中文字幕| 国产精品久久久久久精品电影小说 | 亚洲欧美清纯卡通| 国产精品一二三区在线看| 中文字幕人妻熟人妻熟丝袜美| 国产毛片在线视频| 黑丝袜美女国产一区| 亚洲av日韩在线播放| 黄片wwwwww| 成人影院久久| 性色avwww在线观看| a级毛色黄片| 大又大粗又爽又黄少妇毛片口| 免费大片黄手机在线观看| 交换朋友夫妻互换小说| 免费播放大片免费观看视频在线观看| 五月玫瑰六月丁香| 亚洲精品久久午夜乱码| 黄色视频在线播放观看不卡| 久久久久久久国产电影| 成人午夜精彩视频在线观看| 日本-黄色视频高清免费观看| 亚洲美女黄色视频免费看| 有码 亚洲区| 男男h啪啪无遮挡| 三级国产精品片| 最近最新中文字幕大全电影3| www.av在线官网国产| 亚洲av欧美aⅴ国产| 97超视频在线观看视频| 亚洲欧美日韩卡通动漫| 久久精品夜色国产| 老师上课跳d突然被开到最大视频| 激情 狠狠 欧美| 一区二区三区免费毛片| 妹子高潮喷水视频| 精品一品国产午夜福利视频| 欧美日韩一区二区视频在线观看视频在线| 日韩伦理黄色片| 成人影院久久| 色婷婷久久久亚洲欧美| 亚洲图色成人| 尾随美女入室| 亚洲av日韩在线播放| 日韩伦理黄色片| 久久久久久久精品精品| 亚洲人与动物交配视频| 久久精品国产a三级三级三级| 欧美激情国产日韩精品一区| 亚洲av欧美aⅴ国产| 亚洲欧洲国产日韩| 亚洲精品国产成人久久av| 国产 精品1| 夜夜爽夜夜爽视频| 国产精品一区二区三区四区免费观看| 只有这里有精品99| 国产欧美另类精品又又久久亚洲欧美| 91久久精品国产一区二区成人| 99热6这里只有精品| 亚洲av成人精品一区久久| 网址你懂的国产日韩在线| 亚洲美女黄色视频免费看| 国产精品欧美亚洲77777| av不卡在线播放| 国产 一区 欧美 日韩| 看免费成人av毛片| 插逼视频在线观看| 婷婷色综合大香蕉| 国产高清国产精品国产三级 | 亚洲久久久国产精品| 日产精品乱码卡一卡2卡三| 各种免费的搞黄视频| 美女cb高潮喷水在线观看| 成年av动漫网址| 少妇高潮的动态图| 22中文网久久字幕| 18禁裸乳无遮挡动漫免费视频| 91久久精品国产一区二区三区| 男人爽女人下面视频在线观看| 一个人免费看片子| av福利片在线观看| 国产一区二区在线观看日韩| 偷拍熟女少妇极品色| 国产美女午夜福利| 国产精品人妻久久久久久| 日本免费在线观看一区| 久久国产乱子免费精品| 国产精品女同一区二区软件| 亚洲国产色片| 小蜜桃在线观看免费完整版高清| 22中文网久久字幕| 亚洲精品色激情综合| 超碰av人人做人人爽久久| 国产精品一二三区在线看| kizo精华| 在线看a的网站| 午夜激情久久久久久久| 97热精品久久久久久| 国产久久久一区二区三区| 亚洲最大成人中文| 亚洲成人中文字幕在线播放| 最近中文字幕高清免费大全6| 午夜激情福利司机影院| 人人妻人人看人人澡| 精品人妻偷拍中文字幕| 男女边吃奶边做爰视频| 精品一区二区三区视频在线| 亚洲人成网站高清观看| 人体艺术视频欧美日本| av免费在线看不卡| 99久国产av精品国产电影| 99久久人妻综合| 久久久久久久久久久免费av| 日本wwww免费看| 免费观看无遮挡的男女| 韩国高清视频一区二区三区| 国产免费又黄又爽又色| 寂寞人妻少妇视频99o| 亚州av有码| 最近手机中文字幕大全| 又黄又爽又刺激的免费视频.| 国产成人精品久久久久久| 在线观看三级黄色| 亚洲欧美清纯卡通| 欧美一级a爱片免费观看看| 身体一侧抽搐| 国产亚洲午夜精品一区二区久久| 日韩 亚洲 欧美在线| 亚洲精华国产精华液的使用体验| 国产精品成人在线| 国产 一区 欧美 日韩| 国产综合精华液| 国产淫片久久久久久久久| 国产免费又黄又爽又色| 男人舔奶头视频| 搡老乐熟女国产| 成年免费大片在线观看| 国产成人freesex在线| 深夜a级毛片| 国产黄片视频在线免费观看| 欧美成人午夜免费资源| 精品国产一区二区三区久久久樱花 | h日本视频在线播放| 高清日韩中文字幕在线| 亚洲高清免费不卡视频| 直男gayav资源| 色综合色国产| 美女国产视频在线观看| 夜夜骑夜夜射夜夜干| 国产免费又黄又爽又色| 亚洲精品aⅴ在线观看| 夫妻午夜视频| 美女高潮的动态| 女人久久www免费人成看片| 国产精品国产三级国产av玫瑰| 毛片女人毛片| videos熟女内射| 亚洲一区二区三区欧美精品| 久久久国产一区二区| 乱系列少妇在线播放| 在线精品无人区一区二区三 | 国产视频内射| 大香蕉久久网| 黑人高潮一二区| 99久国产av精品国产电影| 国产在线男女| 观看av在线不卡| 在线观看免费高清a一片| 丝袜喷水一区| 亚洲av中文av极速乱| 青青草视频在线视频观看| 国产在线一区二区三区精| 成人毛片a级毛片在线播放| 深爱激情五月婷婷| 在线亚洲精品国产二区图片欧美 | 国产精品一二三区在线看| 网址你懂的国产日韩在线| 日日撸夜夜添| 乱系列少妇在线播放| 国产探花极品一区二区| 久久久久久久久大av| 国产亚洲精品久久久com| 秋霞伦理黄片| 精华霜和精华液先用哪个| 大片免费播放器 马上看| 国产日韩欧美在线精品| 亚洲精品日韩在线中文字幕| 日日摸夜夜添夜夜爱| 成年女人在线观看亚洲视频| 3wmmmm亚洲av在线观看| 一本久久精品| 欧美老熟妇乱子伦牲交| 99热国产这里只有精品6| 一区二区三区免费毛片| av专区在线播放| 深夜a级毛片| 日韩欧美精品免费久久| a级一级毛片免费在线观看| av播播在线观看一区| 日本爱情动作片www.在线观看| 高清在线视频一区二区三区| 天天躁夜夜躁狠狠久久av| 女人十人毛片免费观看3o分钟| 一级爰片在线观看| 国产黄频视频在线观看| 欧美成人一区二区免费高清观看| 免费人成在线观看视频色| 精品人妻偷拍中文字幕| 久久99热这里只频精品6学生| 国产亚洲av片在线观看秒播厂| 亚洲精品亚洲一区二区| 天天躁日日操中文字幕| 国产日韩欧美在线精品| av在线观看视频网站免费| 国产精品久久久久久av不卡| 精品一品国产午夜福利视频| av国产免费在线观看| 热99国产精品久久久久久7| 在线播放无遮挡| 日韩大片免费观看网站| 国产中年淑女户外野战色| 国产精品99久久久久久久久| 在线看a的网站| 国产有黄有色有爽视频| 日韩成人伦理影院| 伊人久久国产一区二区| 亚洲人成网站在线播| 一区二区三区乱码不卡18| 免费不卡的大黄色大毛片视频在线观看| 亚洲高清免费不卡视频| 精品一区二区三卡| 国产 一区精品| 一本一本综合久久| 免费看光身美女| 一本—道久久a久久精品蜜桃钙片| 午夜福利视频精品| 日韩电影二区| 日本午夜av视频| 亚洲四区av| 国产伦在线观看视频一区| 十分钟在线观看高清视频www | 久久久久视频综合| 三级国产精品片| 欧美区成人在线视频| 日韩强制内射视频| 久久久久久久大尺度免费视频| 天堂中文最新版在线下载| 超碰av人人做人人爽久久| 国产日韩欧美亚洲二区| 一区二区三区乱码不卡18| 在线观看免费视频网站a站| 亚洲精品,欧美精品| 国模一区二区三区四区视频| h视频一区二区三区| 这个男人来自地球电影免费观看 | 欧美 日韩 精品 国产| 高清午夜精品一区二区三区| 国产精品精品国产色婷婷| 极品少妇高潮喷水抽搐| 自拍偷自拍亚洲精品老妇| 大香蕉97超碰在线| 亚洲国产日韩一区二区| 一级a做视频免费观看| 久久精品国产亚洲av天美| 日本黄大片高清| 亚洲国产精品一区三区| 亚洲第一av免费看| 精品久久久久久久久亚洲| 亚洲av免费高清在线观看| 国产成人freesex在线| 丰满乱子伦码专区| 精品亚洲成a人片在线观看 | 深夜a级毛片| 综合色丁香网| 久久婷婷青草| 少妇人妻久久综合中文| 韩国av在线不卡| 欧美三级亚洲精品| 久久青草综合色| 国产午夜精品一二区理论片| 亚洲欧洲国产日韩| 免费黄色在线免费观看| 国产高清不卡午夜福利| 久久影院123| 国产色婷婷99| 少妇被粗大猛烈的视频| 久久久久精品性色| 国产在线免费精品| 久久久久人妻精品一区果冻| 大又大粗又爽又黄少妇毛片口| 国产亚洲5aaaaa淫片| videossex国产| 久久精品国产亚洲av涩爱| 2018国产大陆天天弄谢| 99久久人妻综合| 中文欧美无线码| 黄色视频在线播放观看不卡| 黑丝袜美女国产一区| av一本久久久久| 久久精品国产亚洲av天美| 亚洲av国产av综合av卡| 最后的刺客免费高清国语| 色5月婷婷丁香| 一边亲一边摸免费视频| 亚洲精品一区蜜桃| 精品人妻视频免费看| 国产成人aa在线观看| 午夜福利在线观看免费完整高清在| 久久精品久久久久久噜噜老黄| 久久久久久伊人网av| 性色av一级| 亚洲精品亚洲一区二区| 免费观看的影片在线观看| 国产黄片美女视频| 午夜免费鲁丝| 青春草视频在线免费观看| 久久精品国产自在天天线| 国产中年淑女户外野战色| 欧美日韩一区二区视频在线观看视频在线| 在线观看美女被高潮喷水网站| 久久国产精品男人的天堂亚洲 | 一区二区av电影网| 99热这里只有是精品在线观看| 日韩欧美 国产精品| 亚洲va在线va天堂va国产| 尾随美女入室| 一边亲一边摸免费视频| av在线老鸭窝| 黄片wwwwww| 三级国产精品片| 汤姆久久久久久久影院中文字幕| 国产精品久久久久久精品电影小说 | 在线看a的网站| 青青草视频在线视频观看| 午夜免费观看性视频| 国产精品秋霞免费鲁丝片| 干丝袜人妻中文字幕| 午夜福利高清视频| av不卡在线播放| 一个人看的www免费观看视频| 综合色丁香网| 尾随美女入室| 久久久久人妻精品一区果冻| 国产精品国产av在线观看| 久久久国产一区二区| 极品少妇高潮喷水抽搐| 3wmmmm亚洲av在线观看| 韩国高清视频一区二区三区| 日韩av在线免费看完整版不卡| 国产老妇伦熟女老妇高清| 99视频精品全部免费 在线| 啦啦啦视频在线资源免费观看| 99re6热这里在线精品视频| av一本久久久久| 午夜免费观看性视频| 中文精品一卡2卡3卡4更新| 国产一级毛片在线| 国产午夜精品久久久久久一区二区三区| 麻豆乱淫一区二区| 日韩一本色道免费dvd| av女优亚洲男人天堂| 亚洲电影在线观看av| 美女福利国产在线 | 亚洲天堂av无毛| 久久99热6这里只有精品| 91午夜精品亚洲一区二区三区| 中文乱码字字幕精品一区二区三区| 国产有黄有色有爽视频| 日韩免费高清中文字幕av| 人体艺术视频欧美日本| 好男人视频免费观看在线| 波野结衣二区三区在线| 美女内射精品一级片tv| 性高湖久久久久久久久免费观看| 久久99热这里只频精品6学生| 欧美人与善性xxx| 人人妻人人添人人爽欧美一区卜 | av国产免费在线观看| 一级av片app| 91aial.com中文字幕在线观看| 秋霞在线观看毛片| 国产极品天堂在线| 少妇 在线观看| av在线蜜桃| 国产精品久久久久久精品电影小说 | 五月伊人婷婷丁香| 香蕉精品网在线| 五月伊人婷婷丁香| 99热这里只有精品一区| 男女下面进入的视频免费午夜| 狠狠精品人妻久久久久久综合| 日韩伦理黄色片| 久久99热6这里只有精品| av播播在线观看一区| 亚洲欧美一区二区三区国产| 午夜福利高清视频| 国产黄色视频一区二区在线观看| 久久精品人妻少妇| 免费在线观看成人毛片| 只有这里有精品99| 国产视频内射| 久久热精品热| 久久精品夜色国产| 日产精品乱码卡一卡2卡三| 欧美丝袜亚洲另类| 日韩欧美 国产精品| 丝袜脚勾引网站| 国产极品天堂在线| 韩国高清视频一区二区三区| 深夜a级毛片| av黄色大香蕉| 久久久国产一区二区| 80岁老熟妇乱子伦牲交| 在线观看美女被高潮喷水网站| 伦精品一区二区三区| 久久6这里有精品| 国产欧美亚洲国产| 卡戴珊不雅视频在线播放| 亚洲中文av在线| 亚洲婷婷狠狠爱综合网| 精品国产一区二区三区久久久樱花 | 免费黄网站久久成人精品| 精品人妻偷拍中文字幕| 国产黄色视频一区二区在线观看| 22中文网久久字幕| 建设人人有责人人尽责人人享有的 | 这个男人来自地球电影免费观看 | 国产欧美亚洲国产| 91精品国产国语对白视频| 欧美区成人在线视频| 我要看黄色一级片免费的| 大话2 男鬼变身卡| 国产精品久久久久久精品电影小说 | 日韩不卡一区二区三区视频在线| 国产在线一区二区三区精| 我的女老师完整版在线观看| 亚洲av成人精品一区久久| 国产亚洲欧美精品永久| 中文字幕久久专区| 午夜激情久久久久久久| 日韩一区二区视频免费看| 美女xxoo啪啪120秒动态图| 亚洲性久久影院| 观看美女的网站| 99re6热这里在线精品视频| 国精品久久久久久国模美| 性高湖久久久久久久久免费观看| 日韩欧美精品免费久久| 我的老师免费观看完整版| 乱系列少妇在线播放| 午夜福利在线在线| 极品教师在线视频| 久久久久久久久大av| 精品视频人人做人人爽| 成人影院久久| www.av在线官网国产| 一本久久精品| 男的添女的下面高潮视频|