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

    十三陵地震臺(tái)gPhone重力儀的儀器性能與水文響應(yīng)分析

    2015-02-15 01:06:32佘雅文付廣裕
    關(guān)鍵詞:重力儀殘差重力

    佘雅文 付廣裕 韋 進(jìn)

    1 中國地震局地震預(yù)測(cè)重點(diǎn)實(shí)驗(yàn)室(中國地震局地震預(yù)測(cè)研究所),北京市復(fù)興路63號(hào),100036

    2 中國地震局地震研究所(地震大地測(cè)量實(shí)驗(yàn)室),武漢市洪山側(cè)路40號(hào),430071

    地下水變化與降雨是影響重力觀測(cè)的重要因素之一。近年來水文影響對(duì)連續(xù)重力觀測(cè)的研究已經(jīng)從經(jīng)驗(yàn)性轉(zhuǎn)變?yōu)槲锢砟P停?]。Tanaka使用兩臺(tái)gPhone獲取的重力殘差時(shí)間序列較好地反映了重力變化和降雨的響應(yīng)關(guān)系,并用Kazama建立的水文模型較好地模擬了重力變化[2-3]。利用降雨和靜水位的變化來確定長(zhǎng)期定點(diǎn)重力觀測(cè)的水文響應(yīng)模型,對(duì)去除干擾、獲取當(dāng)?shù)氐貧?nèi)部物質(zhì)遷移以及深部構(gòu)造活動(dòng)具有重要意義[4-5]。

    本文利用十三陵地震臺(tái)的兩臺(tái)gPhone重力儀(序號(hào)109、118)2013-04~08 連續(xù)觀測(cè)數(shù)據(jù)進(jìn)行潮汐分析和提取重力殘差處理,從而對(duì)儀器的性能進(jìn)行判定,并在此基礎(chǔ)上分析gPhone連續(xù)重力觀測(cè)的水文響應(yīng)特征。

    1 連續(xù)重力觀測(cè)

    十三陵地震臺(tái)位于山區(qū)向平原過渡的斜坡地帶。兩臺(tái)gPhone重力儀并排放置在近似水平的平整瓷磚上,儀器旁放置數(shù)據(jù)采集和網(wǎng)絡(luò)通信設(shè)備,儀器和地面的接觸較為平整。利用兩臺(tái)重力儀在同一觀測(cè)條件下得到的數(shù)據(jù)來驗(yàn)證數(shù)據(jù)的可靠性,排除儀器自身原因造成的數(shù)據(jù)失真。由于斷電等原因,2013-06-19~07-04數(shù)據(jù)缺失。

    2 數(shù)據(jù)處理

    儀器原始數(shù)據(jù)為秒采樣,利用加漢寧窗的數(shù)字低通濾波器將秒采樣數(shù)據(jù)濾波為分鐘采樣,使用Tsoft軟件以480為窗口數(shù)、12cpd為采樣率將分鐘采樣降為小時(shí)采樣。采用Tsoft預(yù)處理軟件[6]對(duì)觀測(cè)曲線中的尖峰、地震和掉格等干擾信號(hào)進(jìn)行手工修正,對(duì)較短的記錄中斷用三次多項(xiàng)式進(jìn)行插值填補(bǔ),保留長(zhǎng)時(shí)間中斷,并把秒采樣數(shù)據(jù)轉(zhuǎn)換成1min至1h采樣,對(duì)小時(shí)值數(shù)據(jù)進(jìn)行極移校正和水平校正。

    利用BAYTAP-G 軟件進(jìn)行數(shù)據(jù)處理分析。BAYTAP-G 軟件可以將重力數(shù)據(jù)分離為隨機(jī)噪聲、固體潮汐、響應(yīng)數(shù)據(jù)項(xiàng)和漂移項(xiàng)4個(gè)部分,同時(shí)給出對(duì)數(shù)據(jù)潮汐因子的估計(jì)值[7]。使用輔助氣壓、傳感器溫度、水平校正值和相對(duì)重力數(shù)據(jù)一起作為輸入,其中,輔助氣壓、傳感器溫度和水平校正值為輔助數(shù)據(jù)。氣壓和傳感器溫度是重力觀測(cè)的重要影響因素,氣壓的影響大于傳感器溫度的影響。傳感器溫度雖然是自相關(guān)的,但其突變會(huì)引起觀測(cè)數(shù)據(jù)μGal級(jí)的變化,即每0.001 ℃可能會(huì)引起2.6μGal的重力變化[3]。選取水平校正作為輔助項(xiàng)是因?yàn)?,雖然儀器本身的預(yù)處理中包括水平校正,但是儀器的傾斜角度是隨時(shí)間不斷變化的,由儀器記錄的水平校正數(shù)據(jù)來看,影響在幾十μGal,且變化無規(guī)律。BAYTAP-G 在分離數(shù)據(jù)的過程中并不會(huì)因?yàn)檩o助數(shù)據(jù)在原始數(shù)據(jù)中不存在而產(chǎn)生錯(cuò)誤結(jié)果,因而忽略輔助數(shù)據(jù)的影響[7]。因?yàn)間Phone的漂移不是線性的[3],所以選擇二階最小二乘法擬合觀測(cè)數(shù)據(jù)。對(duì)分離出的漂移信號(hào)作二階最小二乘法擬合,用漂移項(xiàng)數(shù)據(jù)減去擬合數(shù)據(jù)即可得到重力殘差[1]。

    3 數(shù)據(jù)分析

    3.1 兩臺(tái)gPhone重力儀觀測(cè)數(shù)據(jù)的對(duì)比分析

    圖1給出了兩臺(tái)gPhone重力儀在2013-04-20~06-19觀測(cè)到的預(yù)處理數(shù)據(jù)、輔助觀測(cè)值和殘差時(shí)間序列,其中1(a)為gPhone重力儀觀測(cè)到的原始數(shù)據(jù)下采樣后的小時(shí)值數(shù)據(jù),1(b)為儀器傳感器溫度,1(c)為gPhone重力儀水平校正值,1(d)為儀器圍壓值,1(e)為重力殘差值。水平校正值和傳感器溫度作為BAYTAP-G的輔助項(xiàng)得到的重力漂移再減去其二階最小二乘法擬合值,得到重力殘差時(shí)間序列(05-09 23:00~05-11 12:00之間gPhone109無測(cè)量數(shù)據(jù))。從圖1(b)和圖1(c)可以看出,gPhone109和gPhone118雖然放在同一個(gè)觀測(cè)室內(nèi),但兩臺(tái)儀器的傳感器溫度和水平校正項(xiàng)有較大的差異。前者在2013-04-20~06-20時(shí)段內(nèi)傳感器溫度僅在53.463 ℃~53.465 ℃變化,且大部分時(shí)間維持恒溫,而后者傳感器溫度在53.094 ℃~53.103 ℃之間呈階梯狀變化,可見118號(hào)儀器的傳感器溫度變化沒有109號(hào)儀器穩(wěn)定。109和118號(hào)儀器的水平校正項(xiàng)變化幅度分別為10μGal和35μGal左右,水平穩(wěn)定性上109 號(hào)儀器也強(qiáng)于118 號(hào)。使用BAYTAP-G 分離出的109號(hào)和118號(hào)儀器漂移項(xiàng)通過一階最小二乘法擬合后計(jì)算的漂移速率分別約為5.2μGal/d和5.9μGal/d,109號(hào)儀器的漂移速率比118號(hào)低,造成這種情況的原因可能是109號(hào)儀器的水平性能和傳感器溫度變化都比118號(hào)穩(wěn)定??傮w而言,109號(hào)gPhone重力儀的性能優(yōu)于118號(hào)重力儀。

    圖2為兩臺(tái)重力儀在2013-04-20~06-20的觀測(cè)數(shù)據(jù)處理得到的重力殘差變化曲線,其中圖2(a)為BAYTAP-G 處理兩臺(tái)重力儀數(shù)據(jù)得到漂移項(xiàng)減去其二階最小二乘法擬合得到的重力殘差,圖2(b)為兩臺(tái)儀器二階重力殘差的三次樣條插值。從圖2(b)可以看出,重力變化的趨勢(shì)總體是一致的,局部略有差異。由于兩臺(tái)儀器外部環(huán)境相同,所以差異來源于儀器本身。從圖2(a)中可以看到,109號(hào)重力儀比118 號(hào)重力儀的重力殘差變化要平緩,且基本看不到周期性變化,而118重力殘差中可看到明顯的周期變化成分。

    圖1 gPhone109和118號(hào)重力儀2013-04-20~06-20預(yù)處理數(shù)據(jù)、輔助觀測(cè)值和殘差時(shí)間序列Fig.1 Time series of preprocessing data,auxiliary data and gravity residual of gPhone meters(serial number109and 118)from 2013-04-20~06-20

    圖2 2013-04-20~06-20兩臺(tái)重力儀殘差時(shí)間序列比較Fig.2 Comparison of time series of the gravity residuals of two gPhone gravimeters from 2013-04-20~06-20

    3.2 兩臺(tái)重力儀記錄數(shù)據(jù)的調(diào)和分析

    選取109和118號(hào)gPhone重力儀2013-04-20~06-20的連續(xù)觀測(cè)記錄進(jìn)行處理和分析。表1是109和118號(hào)gPhone重力儀記錄的重力數(shù)據(jù)通過下采樣為小時(shí)值后,使用BAYTAP-G[8]以及ETERNA33[9]進(jìn)行潮汐分析的結(jié)果。

    表1 BAYTAP-G、ETERNA33和DEHANT理論潮汐因子對(duì)比Tab.1 Comparison of tidal factors of BAYTAP-G,ETERNA33and DEHANT

    從表1可見,兩臺(tái)儀器的數(shù)據(jù)使用不同處理軟件得到的潮汐振幅因子及其中誤差較為接近,證明了本文數(shù)據(jù)處理的正確性。為了提供對(duì)比參考,將BAYTAP-G 和ETERNA33的振幅因子與DEHANT 理論模型振幅因子[10]進(jìn)行比較??梢钥闯?,BAYTAP-G 處理得到的潮汐振幅因子整體較為接近理論值,其原因可能是BAYTAP-G可使用輔助數(shù)據(jù)抑制數(shù)據(jù)校正不充分在潮汐分析中的干擾。此外,109號(hào)儀器與118號(hào)儀器的潮汐振幅因子相比整體更接近于理論值。

    綜上,109號(hào)gPhone儀器整體性能優(yōu)于118號(hào)儀器。由于這個(gè)原因,之后有關(guān)水文活動(dòng)造成重力變化的討論都將圍繞109號(hào)gPhone重力儀的觀測(cè)數(shù)據(jù)進(jìn)行。

    3.3 降雨造成的重力變化分析

    圖3(a)和圖3(b)分別給出了2013-05-25~06-19和2013-07-05~08-04的109號(hào)gPhone重力儀重力殘差變化和降雨的時(shí)間序列圖。從圖中可以看到,大部分重力變化與降雨存在較好的對(duì)應(yīng)關(guān)系。大部分情況下,它們之間的對(duì)應(yīng)關(guān)系滿足降雨時(shí)重力變小、降雨后重力逐漸變大的規(guī)律。但是也有少數(shù)情況不符合這種模式,這可能是由于其他未知因素的重力變化引起的,其原因尚待進(jìn)一步研究。

    為解釋這種重力隨降雨變化的特征,建立簡(jiǎn)單的物理模型(圖4(a))。圖4(b)中箭頭所指位置為109和118號(hào)gPhone重力儀的觀測(cè)室,觀測(cè)室西側(cè)緊鄰海拔為161m 的龍山,觀測(cè)室的海拔為101m。圖4(a)中,降雨開始時(shí)水層覆蓋在地表,由于龍山海拔高于觀測(cè)室,在龍山上的雨水對(duì)重力儀產(chǎn)生引力傾斜向上,與重力儀處于同一高度的雨水對(duì)重力儀的引力相互抵消,這時(shí)重力值變小。降雨結(jié)束后,由于山體的坡度,雨水更多地流向地勢(shì)低的地方,滲透到低于重力儀海拔的地下,這時(shí)重力值增大。Naujoks使用超導(dǎo)重力儀也觀測(cè)到了類似現(xiàn)象[4]。

    圖3 gPhone109二階最小二乘法擬合殘差與降雨對(duì)應(yīng)的時(shí)間序列Fig.3 Correlation of time series of 2nd-order curve fitting residual gravity of gPhone meter(serial number 109)and rainfall

    圖4 重力隨降雨變化的物理模型示意圖和儀器位置與周邊環(huán)境Fig.4 Response model of gravity to rainfall and the position and surrounding environment of gPhone meter

    3.4 重力變化與地下水靜水位變化的對(duì)應(yīng)關(guān)系

    圖5是119號(hào)gPhone重力儀觀測(cè)結(jié)果的二階最小二乘法擬合殘差和沙河地震臺(tái)靜水位關(guān)系圖。沙河地震臺(tái)與十三陵地震臺(tái)的距離約為15 km。由圖5可知,連續(xù)重力觀測(cè)殘差曲線與降雨量之間存在較好的相關(guān)性,水位下降時(shí)重力減小,水位上升時(shí)重力增大。可見在對(duì)連續(xù)重力觀測(cè)進(jìn)行分析時(shí),地下水的影響不可忽視,因此建議在全國的連續(xù)重力觀測(cè)網(wǎng)中配備地下水位觀測(cè)。

    4月左右的沙河靜水位下降在重力觀測(cè)信號(hào)中無顯示,可能是由于gPhone重力儀在03-24才啟動(dòng)觀測(cè),尚未完全穩(wěn)定所致。由圖5 可知,gPhone重力觀測(cè)大約需要1~2個(gè)月時(shí)間來穩(wěn)定信號(hào),之后才能進(jìn)行較好的分析。109號(hào)gPhone重力儀觀測(cè)到的重力變化和沙河地震臺(tái)的靜水位存在4d左右的延遲對(duì)應(yīng)關(guān)系。延遲的原因可能是地下水流動(dòng)造成的,由于地下結(jié)構(gòu)比較復(fù)雜,地下水的活動(dòng)受很多地質(zhì)環(huán)境的影響,如地下水流動(dòng)、斷層的阻斷等,都會(huì)造成重力變化對(duì)水文活動(dòng)的延遲[4,11]。

    圖5 gPhone二階最小二乘法擬合殘差與沙河地震臺(tái)靜水位Fig.5 Time series of 2nd-order curve fitting residual gravity and static water level of Shahe seismic station

    4 結(jié) 語

    本文對(duì)十三陵地震臺(tái)的兩臺(tái)gPhone重力儀(109和118號(hào))在2013-04~08的連續(xù)觀測(cè)數(shù)據(jù)進(jìn)行潮汐分析和提取重力殘差處理,對(duì)儀器的性能進(jìn)行判定,并在此基礎(chǔ)上分析gPhone重力觀測(cè)的水文響應(yīng)。結(jié)果顯示:1)影響gPhone重力儀觀測(cè)精度的儀器自身因素主要來自傳感器溫度和儀器的水平穩(wěn)定性。2)利用BAYTAP-G 軟件和ETERNA33軟件獲得的潮汐分析結(jié)果大體一致,證實(shí)了兩套處理軟件的一致性。相對(duì)而言,BAYTAP-G 得到的主要振幅因子更接近DEHANT 給出的理論值。3)依據(jù)109號(hào)gPhone重力儀觀測(cè)結(jié)果給出的二階重力殘差變化和十三陵地震臺(tái)的降雨有μGal量級(jí)的對(duì)應(yīng)關(guān)系,并且降雨開始時(shí)重力值減小,降雨結(jié)束后重力值增大。該現(xiàn)象可根據(jù)觀測(cè)站周圍特定的地形條件定性解釋。4)109號(hào)gPhone重力儀獲取的二階重力殘差變化和距離十三陵地震臺(tái)15km 的沙河地震臺(tái)靜水位變化存在4d左右延遲的正相關(guān)對(duì)應(yīng)關(guān)系,該對(duì)應(yīng)關(guān)系顯示了地下水位觀測(cè)對(duì)重力觀測(cè)解釋的必要性,而4d左右的延遲則是由于重力觀測(cè)與地下水位觀測(cè)點(diǎn)存在15km 的差異所致。重力變化與地下水位變化之間的正相關(guān)對(duì)應(yīng)關(guān)系表明,在全國性連續(xù)重力觀測(cè)網(wǎng)中,有必要進(jìn)一步配備地下水位觀測(cè),以提高連續(xù)重力觀測(cè)的數(shù)據(jù)處理精度。

    [1]Tanaka T,Miyajima R,Asai H,et al.Hydrological Gravity Response Detection Using agPhone Below-and Aboveground[J].Earth,Planets and Space,2013,65(2):59-66

    [2]Kazama T,Okubo S.Hydrological Modeling of Groundwater Disturbances to Observed Gravity:Theory and Application to Asama Volcano,Central Japan[J].Journal of Geophysical Research:Solid Earth,2009,114(B8)

    [3]Tanaka,T.Introduction/Adjustment of gPhone and Its Present Condition of the Data[J].Rep Tono Res Inst Earthq Sci,2010,25:75-81

    [4]Naujoks M,Kroner C,Weise A,et al.Evaluating Local Hydrological Modelling by Temporal Gravity Observations and a Gravimetric Three-Dimensional Model[J].Geophysical Journal International,2010,182(1):233-249

    [5]韋進(jìn),李輝,劉子維,等.武漢九峰地震臺(tái)超導(dǎo)重力儀觀測(cè)分析研究[J].地球物理學(xué)報(bào),2012(6):1 894-1 902(Wei Jin,Li Hui,Liu Ziwei,et al.Observation of Superconducting Gravimeter at Jiufeng Seismic Station[J].Chinese J Geophys,2012(6):1 894-1 902)

    [6]Camp M V,Vauterin P.Tsoft:Graphical and Interactive Software for the Analysis of Time Series and Earth Tides[J].Computers &Geosciences,2005,31(5):631-640

    [7]Tamura Y,Sato T,Ooe M,et al.A Procedure for Tidal Analysis with a Bayesian Information Criterion[J].Geophysical Journal International,1991,104(3):507-516

    [8]Akaike H.Likelihood and the Bayes Procedure[J].Trabajos de Estadística Y de Investigación Operativa,1980,31(1):143-166

    [9]Wenzel H G.The Nanogal Software:Earth Tide Data Processing Package ETERNA 3.30[J].Bull Inf Marées Terrestres,1996,124:9 425-9 439

    [10]Dehant V.Tidal Parameters for an Inelastic Earth[J].Physics of the Earth and Planetary Interiors,1987,49(1),97-116

    [11]Pool D R,Schmidt W.Measurement of Ground-Water Storage Change and Specific Yield Using the Temporal-Gravity Method near Rillito Creek,Tucson,Arizona[R].Tucson:US Geological Survey,1997

    猜你喜歡
    重力儀殘差重力
    瘋狂過山車——重力是什么
    基于雙向GRU與殘差擬合的車輛跟馳建模
    工程化原子重力儀綜述
    gPhone重力儀的面波頻段響應(yīng)實(shí)測(cè)研究
    地震研究(2021年1期)2021-04-13 01:04:56
    基于殘差學(xué)習(xí)的自適應(yīng)無人機(jī)目標(biāo)跟蹤算法
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    基于組合滑??刂频慕^對(duì)重力儀兩級(jí)主動(dòng)減振設(shè)計(jì)
    仰斜式重力擋土墻穩(wěn)定計(jì)算復(fù)核
    一張紙的承重力有多大?
    平穩(wěn)自相關(guān)過程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    男人狂女人下面高潮的视频| 日日摸夜夜添夜夜爱| 日韩三级伦理在线观看| 亚洲精品乱久久久久久| av免费在线看不卡| 大码成人一级视频| 日韩av免费高清视频| 九色成人免费人妻av| 欧美国产精品一级二级三级 | 国产一区二区在线观看日韩| 黑丝袜美女国产一区| 伊人久久国产一区二区| 欧美 日韩 精品 国产| 日韩亚洲欧美综合| 精品久久久精品久久久| 高清不卡的av网站| 男女边吃奶边做爰视频| 色婷婷av一区二区三区视频| 亚洲美女黄色视频免费看| 免费人妻精品一区二区三区视频| 国产精品一区二区三区四区免费观看| 少妇猛男粗大的猛烈进出视频| 中国国产av一级| 偷拍熟女少妇极品色| 黄色毛片三级朝国网站 | 两个人免费观看高清视频 | 久久精品久久精品一区二区三区| 人妻人人澡人人爽人人| 久久 成人 亚洲| 国产男人的电影天堂91| 视频区图区小说| 精品久久久久久久久亚洲| 久久久午夜欧美精品| 国产精品一区二区性色av| 欧美最新免费一区二区三区| 日本黄大片高清| 日本黄大片高清| 汤姆久久久久久久影院中文字幕| 色婷婷av一区二区三区视频| 高清av免费在线| 我要看日韩黄色一级片| 天美传媒精品一区二区| 校园人妻丝袜中文字幕| 一边亲一边摸免费视频| 大片电影免费在线观看免费| 如日韩欧美国产精品一区二区三区 | 亚州av有码| xxx大片免费视频| 亚洲自偷自拍三级| 午夜激情久久久久久久| 桃花免费在线播放| 日本wwww免费看| 日日撸夜夜添| 日本猛色少妇xxxxx猛交久久| 欧美日韩国产mv在线观看视频| 欧美成人精品欧美一级黄| 国产精品国产av在线观看| 国产精品一区二区性色av| 亚洲一级一片aⅴ在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲av综合色区一区| 久久毛片免费看一区二区三区| 亚洲欧美精品自产自拍| 少妇的逼好多水| 色视频在线一区二区三区| 国产美女午夜福利| 在线观看免费视频网站a站| 精品熟女少妇av免费看| 在线观看三级黄色| 国产淫片久久久久久久久| 久久久a久久爽久久v久久| 人人妻人人澡人人爽人人夜夜| 五月伊人婷婷丁香| 在线观看美女被高潮喷水网站| 在线免费观看不下载黄p国产| 日韩视频在线欧美| 性色avwww在线观看| 午夜福利影视在线免费观看| 国产精品人妻久久久影院| 一区在线观看完整版| 欧美日韩精品成人综合77777| 亚洲av国产av综合av卡| 永久网站在线| 在线天堂最新版资源| 欧美xxⅹ黑人| 婷婷色麻豆天堂久久| 欧美bdsm另类| 亚洲av成人精品一二三区| 精品久久久久久久久亚洲| 精品一区二区免费观看| 少妇的逼好多水| 国产片特级美女逼逼视频| 欧美一级a爱片免费观看看| 日本-黄色视频高清免费观看| 中国三级夫妇交换| 精品人妻熟女av久视频| 少妇的逼水好多| 国产永久视频网站| 一级毛片久久久久久久久女| 欧美日韩av久久| 青春草国产在线视频| 久久毛片免费看一区二区三区| 国产女主播在线喷水免费视频网站| 一区在线观看完整版| 我的女老师完整版在线观看| av国产精品久久久久影院| 亚洲高清免费不卡视频| 午夜精品国产一区二区电影| 亚洲成色77777| 日本wwww免费看| 美女内射精品一级片tv| 91成人精品电影| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦在线观看免费高清www| 亚洲四区av| 亚洲av成人精品一区久久| 欧美日韩综合久久久久久| 少妇人妻久久综合中文| 插阴视频在线观看视频| 99久久精品一区二区三区| av网站免费在线观看视频| 九色成人免费人妻av| 国产黄色视频一区二区在线观看| 久久这里有精品视频免费| 在线观看av片永久免费下载| 视频中文字幕在线观看| 丝袜在线中文字幕| 人人妻人人澡人人爽人人夜夜| 亚洲国产精品专区欧美| 日韩免费高清中文字幕av| 国产日韩欧美亚洲二区| 高清黄色对白视频在线免费看 | 欧美3d第一页| 国产成人freesex在线| 五月伊人婷婷丁香| 最近手机中文字幕大全| 欧美变态另类bdsm刘玥| 欧美亚洲 丝袜 人妻 在线| 久久综合国产亚洲精品| 最近中文字幕高清免费大全6| 老司机亚洲免费影院| 日韩精品免费视频一区二区三区 | 国产精品麻豆人妻色哟哟久久| 肉色欧美久久久久久久蜜桃| 色吧在线观看| 熟女人妻精品中文字幕| 老女人水多毛片| 久久av网站| a级毛片在线看网站| 日韩人妻高清精品专区| 一级,二级,三级黄色视频| 免费观看a级毛片全部| 亚洲一级一片aⅴ在线观看| 欧美另类一区| 人人澡人人妻人| 国产亚洲91精品色在线| 精华霜和精华液先用哪个| 亚洲激情五月婷婷啪啪| 色网站视频免费| 国产毛片在线视频| 欧美精品国产亚洲| 亚洲精品日韩av片在线观看| av网站免费在线观看视频| 国产精品不卡视频一区二区| 国产爽快片一区二区三区| 久久午夜福利片| 一级a做视频免费观看| 男女边摸边吃奶| 熟女av电影| 乱人伦中国视频| 国产黄色视频一区二区在线观看| 久久久久久久精品精品| xxx大片免费视频| 国产 一区精品| 免费观看av网站的网址| 不卡视频在线观看欧美| 一个人免费看片子| 日韩视频在线欧美| 日韩人妻高清精品专区| 美女福利国产在线| 国产爽快片一区二区三区| 99国产精品免费福利视频| 欧美日韩在线观看h| 久久97久久精品| 欧美日韩亚洲高清精品| 国产精品久久久久久久电影| 丰满乱子伦码专区| 亚洲图色成人| 国产精品一区二区性色av| 欧美精品国产亚洲| www.av在线官网国产| 中文天堂在线官网| 日韩强制内射视频| 欧美成人午夜免费资源| 久久久亚洲精品成人影院| 成人18禁高潮啪啪吃奶动态图 | 男人爽女人下面视频在线观看| 国产av一区二区精品久久| 国产成人精品婷婷| 啦啦啦啦在线视频资源| 亚洲欧美精品自产自拍| 啦啦啦中文免费视频观看日本| 丝瓜视频免费看黄片| 中文字幕人妻丝袜制服| 亚洲精品久久午夜乱码| 日本欧美视频一区| 观看av在线不卡| 亚洲综合精品二区| 欧美变态另类bdsm刘玥| 七月丁香在线播放| 成人二区视频| 中文天堂在线官网| 18禁在线无遮挡免费观看视频| 欧美97在线视频| av在线观看视频网站免费| 国产av国产精品国产| 国产精品麻豆人妻色哟哟久久| 日韩一本色道免费dvd| av在线播放精品| 在线观看人妻少妇| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久久av不卡| 亚洲三级黄色毛片| 视频区图区小说| av女优亚洲男人天堂| 久久久久久久久久成人| 777米奇影视久久| 青春草国产在线视频| 国产精品一区二区在线观看99| 日韩中文字幕视频在线看片| 97精品久久久久久久久久精品| 亚洲国产精品999| 国产伦精品一区二区三区视频9| 精品一区二区三区视频在线| 久久久精品免费免费高清| 日本黄色日本黄色录像| 免费看不卡的av| 日本欧美视频一区| av女优亚洲男人天堂| 午夜老司机福利剧场| 视频中文字幕在线观看| 夫妻午夜视频| 久久久久视频综合| 亚洲美女黄色视频免费看| 国产熟女欧美一区二区| 成年人免费黄色播放视频 | 少妇被粗大的猛进出69影院 | 插逼视频在线观看| 亚洲美女视频黄频| 一级二级三级毛片免费看| 春色校园在线视频观看| 香蕉精品网在线| 精品一品国产午夜福利视频| 国产在线免费精品| 日本午夜av视频| 日韩免费高清中文字幕av| 亚洲精品一二三| 桃花免费在线播放| 久久久久久久大尺度免费视频| 人妻夜夜爽99麻豆av| 国产综合精华液| 免费少妇av软件| 免费观看a级毛片全部| 老司机亚洲免费影院| 欧美日本中文国产一区发布| 亚洲自偷自拍三级| 久久精品国产鲁丝片午夜精品| 国产精品国产三级专区第一集| 成人无遮挡网站| a级片在线免费高清观看视频| www.色视频.com| 国产精品国产三级国产专区5o| 一级毛片黄色毛片免费观看视频| 成人漫画全彩无遮挡| 亚洲国产毛片av蜜桃av| 亚洲国产av新网站| 国产精品久久久久久久久免| 寂寞人妻少妇视频99o| 免费观看性生交大片5| 国产淫语在线视频| 亚洲国产精品专区欧美| 亚洲国产成人一精品久久久| 日韩人妻高清精品专区| 国产精品嫩草影院av在线观看| www.av在线官网国产| 久久久精品94久久精品| 黄色欧美视频在线观看| 美女国产视频在线观看| 久久久久久久大尺度免费视频| 在线观看av片永久免费下载| 午夜老司机福利剧场| 欧美激情国产日韩精品一区| 国产国拍精品亚洲av在线观看| 亚洲av不卡在线观看| 在线观看av片永久免费下载| 婷婷色综合www| 看免费成人av毛片| 青青草视频在线视频观看| 日韩视频在线欧美| 18禁裸乳无遮挡动漫免费视频| 嘟嘟电影网在线观看| 纵有疾风起免费观看全集完整版| 99九九在线精品视频 | 亚洲成人手机| 免费大片黄手机在线观看| 亚洲欧洲日产国产| 大又大粗又爽又黄少妇毛片口| 人人妻人人澡人人看| 免费黄色在线免费观看| 久久99一区二区三区| 亚洲高清免费不卡视频| 女人久久www免费人成看片| 久久午夜综合久久蜜桃| 国产欧美日韩精品一区二区| 99久久精品热视频| 午夜精品国产一区二区电影| 亚洲精品亚洲一区二区| 水蜜桃什么品种好| 交换朋友夫妻互换小说| 老司机亚洲免费影院| av天堂中文字幕网| 国产成人免费无遮挡视频| 97超碰精品成人国产| 大码成人一级视频| 又爽又黄a免费视频| 哪个播放器可以免费观看大片| 国产成人aa在线观看| 夜夜看夜夜爽夜夜摸| 少妇猛男粗大的猛烈进出视频| 国产毛片在线视频| 欧美性感艳星| 亚洲精品一二三| 亚洲内射少妇av| 一个人免费看片子| 国产精品国产三级国产av玫瑰| 三上悠亚av全集在线观看 | 国产深夜福利视频在线观看| 国产 一区精品| 亚洲精品第二区| 永久网站在线| 热99国产精品久久久久久7| 亚洲成色77777| av在线播放精品| 亚洲性久久影院| 亚洲精品国产色婷婷电影| 欧美97在线视频| 99久久中文字幕三级久久日本| 午夜精品国产一区二区电影| 91aial.com中文字幕在线观看| 日韩 亚洲 欧美在线| 亚洲第一区二区三区不卡| videossex国产| 欧美激情极品国产一区二区三区 | 国产欧美另类精品又又久久亚洲欧美| 免费黄频网站在线观看国产| 国产一区二区在线观看av| 丝袜在线中文字幕| 婷婷色综合大香蕉| 香蕉精品网在线| 亚洲精品,欧美精品| 午夜激情福利司机影院| 日韩人妻高清精品专区| 欧美xxxx性猛交bbbb| 寂寞人妻少妇视频99o| 一本久久精品| 国产精品久久久久久久电影| 午夜91福利影院| 亚洲经典国产精华液单| 久久韩国三级中文字幕| 日韩熟女老妇一区二区性免费视频| 色视频在线一区二区三区| 国产精品人妻久久久久久| 精品久久久噜噜| 十分钟在线观看高清视频www | 汤姆久久久久久久影院中文字幕| 国产 精品1| 亚洲精品aⅴ在线观看| 亚洲av免费高清在线观看| 熟女电影av网| 国产亚洲精品久久久com| av天堂中文字幕网| 一本久久精品| 寂寞人妻少妇视频99o| 国产精品久久久久成人av| 久久国产精品大桥未久av | 永久网站在线| 日本欧美国产在线视频| 最近中文字幕2019免费版| 嫩草影院入口| 99精国产麻豆久久婷婷| 亚洲va在线va天堂va国产| av有码第一页| 精品国产国语对白av| 黄色日韩在线| 精品久久久精品久久久| 欧美日韩一区二区视频在线观看视频在线| 大香蕉97超碰在线| 日韩欧美 国产精品| 色婷婷久久久亚洲欧美| 22中文网久久字幕| 狂野欧美激情性bbbbbb| 亚洲精品亚洲一区二区| 国产精品一二三区在线看| 亚洲精品,欧美精品| 久久99热6这里只有精品| 国产精品偷伦视频观看了| 国产精品一二三区在线看| 自拍偷自拍亚洲精品老妇| 人妻系列 视频| 在线观看免费日韩欧美大片 | 色5月婷婷丁香| 久久综合国产亚洲精品| 亚洲欧美成人精品一区二区| 美女xxoo啪啪120秒动态图| 中文资源天堂在线| 久久精品久久精品一区二区三区| 看免费成人av毛片| 好男人视频免费观看在线| 老司机影院毛片| 99久久综合免费| 亚洲av二区三区四区| 国产精品久久久久久av不卡| 国产精品久久久久久久电影| 热99国产精品久久久久久7| 纵有疾风起免费观看全集完整版| 国产 精品1| 亚洲av国产av综合av卡| 国产无遮挡羞羞视频在线观看| 91aial.com中文字幕在线观看| 久热久热在线精品观看| 亚洲国产av新网站| 一级爰片在线观看| 日韩av免费高清视频| 赤兔流量卡办理| 成人影院久久| 亚洲国产日韩一区二区| av国产久精品久网站免费入址| 欧美区成人在线视频| 99久久精品国产国产毛片| 午夜福利在线观看免费完整高清在| 五月天丁香电影| 欧美精品一区二区大全| 99热这里只有是精品在线观看| 欧美亚洲 丝袜 人妻 在线| 超碰97精品在线观看| 欧美日韩精品成人综合77777| 欧美成人午夜免费资源| 日韩三级伦理在线观看| 99久久精品一区二区三区| 我要看日韩黄色一级片| 黄色一级大片看看| 免费久久久久久久精品成人欧美视频 | 久久精品国产自在天天线| 99热这里只有精品一区| 国产成人免费观看mmmm| 男女无遮挡免费网站观看| 成人特级av手机在线观看| 尾随美女入室| 少妇猛男粗大的猛烈进出视频| 色网站视频免费| 老司机亚洲免费影院| 亚洲精品亚洲一区二区| 如何舔出高潮| 丰满人妻一区二区三区视频av| 美女中出高潮动态图| 成人影院久久| 女人久久www免费人成看片| 插阴视频在线观看视频| 天天躁夜夜躁狠狠久久av| 韩国高清视频一区二区三区| 成人国产麻豆网| 国产一区有黄有色的免费视频| 男女边吃奶边做爰视频| 在线精品无人区一区二区三| 男男h啪啪无遮挡| 亚洲电影在线观看av| 国产 精品1| 黄片无遮挡物在线观看| 亚洲av在线观看美女高潮| 精品午夜福利在线看| 亚洲人成网站在线观看播放| 久久6这里有精品| 九九在线视频观看精品| 狂野欧美激情性bbbbbb| 精品久久久久久久久亚洲| 国产日韩欧美亚洲二区| 黑人猛操日本美女一级片| 99国产精品免费福利视频| 久久午夜福利片| 成人毛片60女人毛片免费| 国产中年淑女户外野战色| 美女主播在线视频| 一个人免费看片子| 国产高清国产精品国产三级| 看十八女毛片水多多多| 老司机影院成人| 只有这里有精品99| 精品99又大又爽又粗少妇毛片| 亚州av有码| 18+在线观看网站| 嫩草影院新地址| 国产av码专区亚洲av| 亚洲欧美一区二区三区国产| 最黄视频免费看| 日本91视频免费播放| av一本久久久久| 国产成人精品福利久久| 99久国产av精品国产电影| 亚洲av国产av综合av卡| 自线自在国产av| 一级片'在线观看视频| 亚洲国产精品999| 国产免费一区二区三区四区乱码| 国内少妇人妻偷人精品xxx网站| 一级av片app| av免费在线看不卡| 欧美高清成人免费视频www| 国产免费视频播放在线视频| 纵有疾风起免费观看全集完整版| 国产精品伦人一区二区| 一本—道久久a久久精品蜜桃钙片| av视频免费观看在线观看| 免费观看无遮挡的男女| 成年美女黄网站色视频大全免费 | 久久毛片免费看一区二区三区| 亚洲,一卡二卡三卡| 一级毛片黄色毛片免费观看视频| 久久精品国产亚洲网站| 亚洲怡红院男人天堂| 成人综合一区亚洲| 黄色怎么调成土黄色| 精品国产露脸久久av麻豆| 熟女人妻精品中文字幕| 国产极品粉嫩免费观看在线 | 欧美高清成人免费视频www| 少妇猛男粗大的猛烈进出视频| 成人毛片60女人毛片免费| 亚洲精品中文字幕在线视频 | 特大巨黑吊av在线直播| 中文精品一卡2卡3卡4更新| 熟女av电影| 午夜91福利影院| 亚洲精品日韩在线中文字幕| 国产深夜福利视频在线观看| 欧美少妇被猛烈插入视频| 五月开心婷婷网| 欧美日韩av久久| 国产一区二区在线观看日韩| 国产淫语在线视频| 免费观看无遮挡的男女| 边亲边吃奶的免费视频| 亚洲国产欧美日韩在线播放 | 人妻一区二区av| 午夜福利影视在线免费观看| 日本与韩国留学比较| 国产乱人偷精品视频| 男女国产视频网站| 久久国产精品大桥未久av | 视频区图区小说| h视频一区二区三区| 黄色毛片三级朝国网站 | 成人二区视频| 国产男人的电影天堂91| 人妻制服诱惑在线中文字幕| 亚洲国产精品成人久久小说| 少妇丰满av| 久久久精品免费免费高清| 精品少妇内射三级| 这个男人来自地球电影免费观看 | 国产白丝娇喘喷水9色精品| 麻豆成人午夜福利视频| 国产欧美日韩综合在线一区二区 | 国产精品秋霞免费鲁丝片| 赤兔流量卡办理| 99九九在线精品视频 | 乱人伦中国视频| av天堂久久9| 在线天堂最新版资源| 免费看av在线观看网站| 777米奇影视久久| 黄色配什么色好看| 男人爽女人下面视频在线观看| 久久精品国产自在天天线| 久久这里有精品视频免费| 热re99久久国产66热| 91午夜精品亚洲一区二区三区| 精品亚洲乱码少妇综合久久| 久久影院123| 欧美三级亚洲精品| 极品教师在线视频| 国产精品三级大全| 国产 精品1| 国产一区二区在线观看av| a级片在线免费高清观看视频| 午夜福利影视在线免费观看| 久久国产乱子免费精品| 日韩精品免费视频一区二区三区 | 精品国产乱码久久久久久小说| 日韩一区二区视频免费看| 午夜福利视频精品| 国产精品国产三级专区第一集| 又黄又爽又刺激的免费视频.| 久久青草综合色| 久久久久久久久久久久大奶| 久久精品国产鲁丝片午夜精品| 一级二级三级毛片免费看| 国产成人精品久久久久久| 午夜精品国产一区二区电影| 亚洲精品国产成人久久av| 日本免费在线观看一区| 女的被弄到高潮叫床怎么办| 一级毛片电影观看| 一级毛片黄色毛片免费观看视频| 极品教师在线视频| 一级毛片电影观看|