• <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
    日本免费a在线| 亚洲人与动物交配视频| 在线播放国产精品三级| 九草在线视频观看| 免费看a级黄色片| 亚洲精品粉嫩美女一区| 能在线免费观看的黄片| 少妇的逼水好多| 欧美区成人在线视频| 日本色播在线视频| 亚洲av一区综合| 午夜激情福利司机影院| 国产亚洲av片在线观看秒播厂 | 精品久久久久久久久亚洲| 91aial.com中文字幕在线观看| 色噜噜av男人的天堂激情| 免费看a级黄色片| av天堂在线播放| av在线观看视频网站免费| 精品人妻视频免费看| 精品国产三级普通话版| 日韩欧美三级三区| 国产一区二区在线av高清观看| 在线观看美女被高潮喷水网站| 好男人视频免费观看在线| 成人性生交大片免费视频hd| 国产麻豆成人av免费视频| 日本黄色片子视频| 日本三级黄在线观看| 国产伦精品一区二区三区四那| 91精品一卡2卡3卡4卡| 伦理电影大哥的女人| 中文字幕av成人在线电影| 天堂√8在线中文| 亚洲av免费高清在线观看| 日本三级黄在线观看| 黄色日韩在线| 日韩欧美三级三区| 国产亚洲精品久久久com| 中文字幕人妻熟人妻熟丝袜美| 老熟妇乱子伦视频在线观看| 国产精品,欧美在线| 午夜福利在线观看吧| 亚洲美女搞黄在线观看| 国产熟女欧美一区二区| 亚州av有码| 男插女下体视频免费在线播放| 国产av一区在线观看免费| 变态另类成人亚洲欧美熟女| 久久精品国产清高在天天线| 亚洲欧美清纯卡通| 久久久久久久亚洲中文字幕| 天堂中文最新版在线下载 | 男人舔女人下体高潮全视频| 神马国产精品三级电影在线观看| av免费在线看不卡| 免费人成视频x8x8入口观看| 一本一本综合久久| 久久久国产成人精品二区| 欧美最黄视频在线播放免费| .国产精品久久| www日本黄色视频网| 噜噜噜噜噜久久久久久91| 禁无遮挡网站| 久久精品夜色国产| 欧美成人精品欧美一级黄| 晚上一个人看的免费电影| 校园春色视频在线观看| 精品人妻熟女av久视频| 亚洲精品国产av成人精品| 少妇人妻精品综合一区二区 | 晚上一个人看的免费电影| 国产精品永久免费网站| 国产探花在线观看一区二区| 夫妻性生交免费视频一级片| 自拍偷自拍亚洲精品老妇| 丰满人妻一区二区三区视频av| 免费人成视频x8x8入口观看| 最近2019中文字幕mv第一页| 99热这里只有精品一区| 99久久成人亚洲精品观看| 精品久久久久久久久久免费视频| 国产探花在线观看一区二区| 在线观看66精品国产| 免费电影在线观看免费观看| 99久国产av精品国产电影| 国产色婷婷99| 村上凉子中文字幕在线| 性插视频无遮挡在线免费观看| 国产久久久一区二区三区| av福利片在线观看| 有码 亚洲区| 在线免费观看不下载黄p国产| 18+在线观看网站| 久久精品人妻少妇| 青春草亚洲视频在线观看| 精品无人区乱码1区二区| 中国美白少妇内射xxxbb| 亚洲成人精品中文字幕电影| 午夜福利在线观看吧| 99久久成人亚洲精品观看| 成人特级黄色片久久久久久久| 波多野结衣高清作品| 中文字幕人妻熟人妻熟丝袜美| h日本视频在线播放| 91av网一区二区| av在线播放精品| 欧美高清性xxxxhd video| 伦理电影大哥的女人| 国产蜜桃级精品一区二区三区| 午夜免费激情av| 久久久精品大字幕| 欧美日本亚洲视频在线播放| 欧美xxxx黑人xx丫x性爽| 成熟少妇高潮喷水视频| 午夜福利在线观看免费完整高清在 | 欧美zozozo另类| 日本三级黄在线观看| 国产极品精品免费视频能看的| 少妇的逼好多水| 热99re8久久精品国产| 欧美人与善性xxx| 国内揄拍国产精品人妻在线| 久99久视频精品免费| 亚洲国产日韩欧美精品在线观看| 亚洲欧洲国产日韩| 夜夜爽天天搞| 可以在线观看的亚洲视频| 成人亚洲精品av一区二区| 久久这里有精品视频免费| 免费无遮挡裸体视频| 国产女主播在线喷水免费视频网站 | 人人妻人人看人人澡| 中文欧美无线码| 日韩成人av中文字幕在线观看| 成人亚洲精品av一区二区| 欧美区成人在线视频| 亚洲熟妇中文字幕五十中出| 成人鲁丝片一二三区免费| 老熟妇乱子伦视频在线观看| 少妇的逼好多水| 久久精品夜色国产| 国产精品一区二区在线观看99 | 身体一侧抽搐| 舔av片在线| 欧美极品一区二区三区四区| 一进一出抽搐gif免费好疼| 舔av片在线| 欧美另类亚洲清纯唯美| 国产精品综合久久久久久久免费| 免费人成视频x8x8入口观看| 亚洲人成网站在线观看播放| 亚洲国产精品sss在线观看| 夜夜爽天天搞| 国产精品不卡视频一区二区| 91麻豆精品激情在线观看国产| 久99久视频精品免费| 特级一级黄色大片| 国产精品,欧美在线| 日韩 亚洲 欧美在线| 日韩国内少妇激情av| 色综合亚洲欧美另类图片| 在现免费观看毛片| 岛国毛片在线播放| 久久久精品94久久精品| 日韩欧美国产在线观看| 2021天堂中文幕一二区在线观| 69av精品久久久久久| 亚洲在线自拍视频| 日韩精品有码人妻一区| 国产麻豆成人av免费视频| 亚洲国产精品合色在线| 亚洲国产欧美人成| 亚洲精品乱码久久久v下载方式| 亚洲第一区二区三区不卡| 欧美一区二区国产精品久久精品| 级片在线观看| 亚洲乱码一区二区免费版| 嫩草影院入口| 国产成人午夜福利电影在线观看| 99热网站在线观看| 国产精品一二三区在线看| 夫妻性生交免费视频一级片| 丰满乱子伦码专区| 亚洲天堂国产精品一区在线| 欧美日韩在线观看h| 最后的刺客免费高清国语| 国产av麻豆久久久久久久| 看非洲黑人一级黄片| 男人和女人高潮做爰伦理| 91精品一卡2卡3卡4卡| 久久久成人免费电影| 午夜爱爱视频在线播放| 亚洲欧美精品自产自拍| 久久久久久久久久黄片| 卡戴珊不雅视频在线播放| 干丝袜人妻中文字幕| 午夜免费男女啪啪视频观看| 99视频精品全部免费 在线| 国产三级中文精品| 久久精品国产清高在天天线| 国模一区二区三区四区视频| 九九爱精品视频在线观看| 九草在线视频观看| 一区福利在线观看| 99久久无色码亚洲精品果冻| 成年av动漫网址| 黄色日韩在线| 久久久久久伊人网av| 午夜老司机福利剧场| 老女人水多毛片| 午夜福利在线观看免费完整高清在 | 亚洲欧美日韩卡通动漫| 亚洲美女搞黄在线观看| 久久精品国产亚洲av涩爱 | 国产一级毛片在线| 亚洲国产精品合色在线| 久久99热这里只有精品18| 亚洲五月天丁香| 久久这里有精品视频免费| 成人国产麻豆网| 久久久国产成人免费| 看免费成人av毛片| 亚洲最大成人av| 天堂影院成人在线观看| 久久国产乱子免费精品| 亚洲人与动物交配视频| 永久网站在线| 激情 狠狠 欧美| 99久国产av精品国产电影| 国产伦精品一区二区三区视频9| 亚洲欧洲日产国产| 欧美成人a在线观看| 女的被弄到高潮叫床怎么办| av福利片在线观看| 午夜a级毛片| 亚洲国产精品成人综合色| 亚洲最大成人中文| 国产精品福利在线免费观看| 免费大片18禁| 极品教师在线视频| 91狼人影院| 日产精品乱码卡一卡2卡三| 综合色av麻豆| 日本黄色片子视频| 看免费成人av毛片| 亚洲欧美成人综合另类久久久 | 免费看美女性在线毛片视频| 国产精品久久久久久久久免| 久久久久国产网址| 精品久久久久久久久亚洲| 嫩草影院入口| 老师上课跳d突然被开到最大视频| av福利片在线观看| av免费在线看不卡| 好男人在线观看高清免费视频| av专区在线播放| 日韩一本色道免费dvd| or卡值多少钱| 熟女人妻精品中文字幕| 成人综合一区亚洲| 国产亚洲精品久久久久久毛片| 日韩制服骚丝袜av| 丝袜美腿在线中文| 级片在线观看| 亚洲一区二区三区色噜噜| 久久久精品94久久精品| 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 免费看美女性在线毛片视频| 亚洲精华国产精华液的使用体验 | 久久久久久久久中文| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 男女啪啪激烈高潮av片| av女优亚洲男人天堂| 小说图片视频综合网站| 精品免费久久久久久久清纯| 亚洲人成网站在线观看播放| 色哟哟·www| 国产大屁股一区二区在线视频| 国产 一区精品| 可以在线观看的亚洲视频| 亚洲国产高清在线一区二区三| 91久久精品国产一区二区成人| 久久久久久久久中文| 国产高清视频在线观看网站| 国产在线精品亚洲第一网站| 国产精品,欧美在线| 亚洲人成网站在线播| 亚洲欧美精品专区久久| 久久精品国产亚洲网站| 久久久久久九九精品二区国产| 日韩高清综合在线| 在线a可以看的网站| 日本免费a在线| 免费人成在线观看视频色| 亚洲av.av天堂| 老司机福利观看| 亚洲乱码一区二区免费版| 男人和女人高潮做爰伦理| 一个人看的www免费观看视频| 深爱激情五月婷婷| 国产午夜精品论理片| 久久久久久国产a免费观看| 久久久精品大字幕| 国产亚洲91精品色在线| 亚洲人成网站高清观看| 亚洲自偷自拍三级| 久久久精品大字幕| 又粗又硬又长又爽又黄的视频 | 国产极品精品免费视频能看的| 18禁在线无遮挡免费观看视频| 亚洲最大成人中文| 亚洲av熟女| 久久久久免费精品人妻一区二区| 精品久久久久久久久av| 99精品在免费线老司机午夜| 免费不卡的大黄色大毛片视频在线观看 | 超碰av人人做人人爽久久| 国产av麻豆久久久久久久| 国产成年人精品一区二区| 国产在视频线在精品| 少妇丰满av| 国产v大片淫在线免费观看| 2022亚洲国产成人精品| 国产日本99.免费观看| 欧美一区二区精品小视频在线| 亚洲av免费高清在线观看| 99视频精品全部免费 在线| 性插视频无遮挡在线免费观看| 国内精品一区二区在线观看| 搡老妇女老女人老熟妇| 国产伦在线观看视频一区| 一级毛片电影观看 | 亚洲欧美精品综合久久99| 黄色配什么色好看| 在线免费观看的www视频| 久久久久九九精品影院| 亚洲欧美日韩高清在线视频| 欧美又色又爽又黄视频| 久久午夜福利片| 精品久久久久久久人妻蜜臀av| 国产一区二区在线观看日韩| 麻豆成人午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看| 91麻豆精品激情在线观看国产| 久久人人爽人人片av| 在线观看66精品国产| 九九热线精品视视频播放| 日本三级黄在线观看| 99久久久亚洲精品蜜臀av| 色播亚洲综合网| 亚洲人成网站在线观看播放| 色噜噜av男人的天堂激情| 成人性生交大片免费视频hd| 少妇熟女aⅴ在线视频| 国产三级在线视频| 中国美女看黄片| 免费人成在线观看视频色| 国产高潮美女av| 99riav亚洲国产免费| av天堂中文字幕网| 午夜老司机福利剧场| 美女脱内裤让男人舔精品视频 | 国产探花极品一区二区| 丝袜美腿在线中文| 小说图片视频综合网站| 免费观看的影片在线观看| 小说图片视频综合网站| 给我免费播放毛片高清在线观看| 亚洲成人久久性| 成人性生交大片免费视频hd| 神马国产精品三级电影在线观看| 欧美丝袜亚洲另类| 亚洲真实伦在线观看| 精品人妻熟女av久视频| 亚洲真实伦在线观看| 午夜激情福利司机影院| 日韩欧美精品v在线| 国产单亲对白刺激| 午夜a级毛片| 别揉我奶头 嗯啊视频| 少妇丰满av| 精品久久久久久成人av| 色噜噜av男人的天堂激情| 女的被弄到高潮叫床怎么办| 搡女人真爽免费视频火全软件| 国产91av在线免费观看| av免费在线看不卡| 成人一区二区视频在线观看| 国产一区二区激情短视频| 国产欧美日韩精品一区二区| 国产精品三级大全| 免费av观看视频| 欧美一区二区亚洲| 欧美激情国产日韩精品一区| 国内精品美女久久久久久| 国产成人影院久久av| 成年女人看的毛片在线观看| 亚洲国产精品成人久久小说 | 国产一级毛片在线| 草草在线视频免费看| 级片在线观看| 少妇的逼好多水| a级一级毛片免费在线观看| 综合色av麻豆| 日本熟妇午夜| 一边摸一边抽搐一进一小说| 人妻夜夜爽99麻豆av| 国产美女午夜福利| 黄片无遮挡物在线观看| 亚洲精品乱码久久久v下载方式| 99热网站在线观看| 99热精品在线国产| 18禁在线播放成人免费| 亚洲色图av天堂| 一区二区三区免费毛片| 寂寞人妻少妇视频99o| 99在线视频只有这里精品首页| 欧美一区二区精品小视频在线| 久久人人爽人人片av| 久久欧美精品欧美久久欧美| 亚洲四区av| 亚洲av二区三区四区| 久久久久久国产a免费观看| av在线老鸭窝| av国产免费在线观看| 国产片特级美女逼逼视频| 大香蕉久久网| or卡值多少钱| 亚洲av成人精品一区久久| 成人国产麻豆网| 国产精品伦人一区二区| 国产精品永久免费网站| 成年女人看的毛片在线观看| 99久久无色码亚洲精品果冻| 色5月婷婷丁香| 变态另类丝袜制服| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产自在天天线| 日本撒尿小便嘘嘘汇集6| 中文字幕av成人在线电影| 99在线视频只有这里精品首页| 又粗又爽又猛毛片免费看| 女人被狂操c到高潮| 成人国产麻豆网| 国产高清三级在线| 色哟哟·www| 欧美+日韩+精品| 亚洲aⅴ乱码一区二区在线播放| 九九在线视频观看精品| 成人av在线播放网站| 亚洲精品久久久久久婷婷小说 | 国产亚洲av片在线观看秒播厂 | 99国产极品粉嫩在线观看| 美女高潮的动态| 草草在线视频免费看| 亚洲在线观看片| or卡值多少钱| 国产亚洲精品av在线| 亚洲最大成人av| 黄色配什么色好看| 亚洲欧美日韩卡通动漫| 三级经典国产精品| 国产精品久久久久久精品电影小说 | 免费观看a级毛片全部| 别揉我奶头 嗯啊视频| 亚洲av成人av| 久久久久久久久久成人| 久久午夜福利片| 91久久精品电影网| 极品教师在线视频| 色吧在线观看| 亚洲人成网站在线观看播放| 丝袜美腿在线中文| 成人毛片60女人毛片免费| 久久99热6这里只有精品| 久久亚洲精品不卡| 美女被艹到高潮喷水动态| 免费观看人在逋| 欧美日韩综合久久久久久| 麻豆国产av国片精品| 自拍偷自拍亚洲精品老妇| 别揉我奶头 嗯啊视频| 精品午夜福利在线看| 国产成人午夜福利电影在线观看| 一个人观看的视频www高清免费观看| 国产午夜福利久久久久久| 国产黄色视频一区二区在线观看 | 一个人观看的视频www高清免费观看| 天堂中文最新版在线下载 | 亚洲成人久久性| 婷婷色av中文字幕| 97超视频在线观看视频| 可以在线观看毛片的网站| 99在线人妻在线中文字幕| 久久久精品大字幕| 免费电影在线观看免费观看| 日本三级黄在线观看| 少妇丰满av| 麻豆av噜噜一区二区三区| 又粗又爽又猛毛片免费看| 国产黄色视频一区二区在线观看 | 免费观看人在逋| 亚洲人成网站在线播放欧美日韩| 日韩一区二区三区影片| 亚洲综合色惰| 99久久精品一区二区三区| 一夜夜www| 日本-黄色视频高清免费观看| 99热全是精品| 国产蜜桃级精品一区二区三区| 天堂av国产一区二区熟女人妻| 我要看日韩黄色一级片| 夫妻性生交免费视频一级片| 美女 人体艺术 gogo| 色噜噜av男人的天堂激情| 免费观看精品视频网站| 午夜亚洲福利在线播放| 欧美日韩一区二区视频在线观看视频在线 | 中文字幕av在线有码专区| 99国产精品一区二区蜜桃av| 久久精品91蜜桃| 欧美色欧美亚洲另类二区| 热99在线观看视频| 久久久久久久久中文| 亚洲精品影视一区二区三区av| 不卡一级毛片| 99久国产av精品| 久久人人爽人人爽人人片va| 成人高潮视频无遮挡免费网站| 久久国内精品自在自线图片| 久久人人精品亚洲av| 国产在视频线在精品| 中国国产av一级| 国产熟女欧美一区二区| 女人被狂操c到高潮| 深夜精品福利| 久久鲁丝午夜福利片| av女优亚洲男人天堂| av国产免费在线观看| 99热精品在线国产| 免费搜索国产男女视频| 91午夜精品亚洲一区二区三区| 日本在线视频免费播放| 国产精品爽爽va在线观看网站| 一级黄色大片毛片| 成年免费大片在线观看| 国产一区亚洲一区在线观看| 亚洲欧美日韩东京热| 老司机影院成人| 26uuu在线亚洲综合色| 国产伦精品一区二区三区四那| 国产成人a∨麻豆精品| 大型黄色视频在线免费观看| 久久精品国产鲁丝片午夜精品| 中国国产av一级| 欧美三级亚洲精品| 亚洲精品456在线播放app| 欧美精品国产亚洲| 伦理电影大哥的女人| 美女 人体艺术 gogo| av免费在线看不卡| 不卡视频在线观看欧美| 国产人妻一区二区三区在| 久久这里只有精品中国| 91精品一卡2卡3卡4卡| 国产精华一区二区三区| www日本黄色视频网| 免费人成在线观看视频色| 男女视频在线观看网站免费| 亚洲精品乱码久久久久久按摩| 国产乱人偷精品视频| 18禁在线播放成人免费| 国产精品一区二区性色av| 99热这里只有精品一区| 舔av片在线| 你懂的网址亚洲精品在线观看 | 成人永久免费在线观看视频| 精品一区二区三区视频在线| 搡女人真爽免费视频火全软件| 91精品一卡2卡3卡4卡| 久久九九热精品免费| 国产一区二区在线观看日韩| 国产午夜精品一二区理论片| 日韩欧美精品免费久久| 日韩av不卡免费在线播放| 日本三级黄在线观看| 偷拍熟女少妇极品色| 国产成人精品久久久久久| 成年av动漫网址| 成人漫画全彩无遮挡| 寂寞人妻少妇视频99o| 欧美高清性xxxxhd video| 中国国产av一级| 色哟哟哟哟哟哟| 日本三级黄在线观看| 国产精品女同一区二区软件| 国产成人精品久久久久久| 九九爱精品视频在线观看| 偷拍熟女少妇极品色| 国产高潮美女av| 午夜久久久久精精品| 两个人视频免费观看高清| 波多野结衣高清无吗| 亚洲人成网站高清观看| 99久久无色码亚洲精品果冻| 免费看a级黄色片| 一级毛片aaaaaa免费看小| 一边摸一边抽搐一进一小说| 哪个播放器可以免费观看大片| 成人午夜精彩视频在线观看| 99riav亚洲国产免费| 看非洲黑人一级黄片| 亚洲成av人片在线播放无|