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

    古登堡-里希特定律中的b值統(tǒng)計樣本量研究1

    2018-11-13 02:09:46李世杰呂悅軍劉靜偉
    震災(zāi)防御技術(shù) 2018年3期
    關(guān)鍵詞:震級樣本量乘法

    李世杰 呂悅軍 劉靜偉

    ?

    古登堡-里希特定律中的值統(tǒng)計樣本量研究1

    李世杰 呂悅軍 劉靜偉

    (中國地震局地殼應(yīng)力研究所,北京 100085)

    值是研究地震活動的重要指標,其廣泛應(yīng)用于地震危險性分析和地震預(yù)測研究之中,與實際資料的完整性、樣本量的大小、計算方法等因素有著重要的關(guān)系。常見的值計算方法有最小二乘法和最大似然法,樣本量的大小對這2種方法影響很大。本文利用蒙特卡羅模擬地震目錄和汾渭地震帶實際目錄作為樣本,從中抽取不同大小的樣本量進行計算,研究不同樣本量下這2種方法計算得到的值與設(shè)定值或真實值之間的差別。結(jié)果表明,最小二乘法需要的最低樣本量為1000,最大似然法為200;當樣本量達不到要求時,計算出的值是不可靠的;由于對樣本量的要求不同,前者適用于計算區(qū)域的整體值,而后者在研究某區(qū)域值在時間軸上的變化方面更有優(yōu)勢。本研究為確定2種值計算方法對樣本量的最低要求提供了參考依據(jù)。

    值 最小二乘法 最大似然法 樣本量大小

    引言

    古登堡和里希特在1941年美國地質(zhì)學(xué)會會刊中提出全球的地震活動服從經(jīng)驗關(guān)系lg=-(式中表示震級,表示震級≥的地震次數(shù),、是常數(shù))。一般認為值代表地區(qū)的地震活動總體水平(Rundle,1989);值代表地震活動的大小地震數(shù)量的比例,是地震活動研究中的重要參量。研究表明,值具有明確的物理意義,與地殼的介質(zhì)特性、應(yīng)力狀態(tài)和不均勻性有關(guān),能反映所研究區(qū)域的地質(zhì)構(gòu)造特征(王熠熙等,2015;謝卓娟等,2015)以及地震的震源特征(Schorlemmer等,2005;Gulia等,2010;劉靜偉等,2016;張廣偉,2016)。因此,值廣泛應(yīng)用于地震危險性分析和地震預(yù)測研究之中。在地震危險性分析中,值和地震年發(fā)生率共同用于確定地震活動的水平(胡聿賢,1999),其取值對地震危險性分析結(jié)果的影響較大(鄢家全等,1996;黃瑋瓊等,1998;謝卓娟等,2013);而在地震預(yù)測研究中,值作為基本的地震活動性參數(shù),成為地震預(yù)測的常用指標參數(shù)(韓渭賓,2003;沈建文等,2007)。

    通常根據(jù)實際地震資料統(tǒng)計得到值。目前,常用的值估算方法是最大似然估計法(Aki,1965)和最小二乘法等,這些方法不僅對實際資料的完整性和精度有一定要求,同時也需要足夠的地震資料。在實際工作中,對于歷史地震資料短缺或地震活動水平低的地區(qū),在計算時常常將現(xiàn)代小震資料與歷史地震資料聯(lián)合使用,以彌補地震資料樣本量的不足(黃瑋瓊等,1989;鄢家全等,1996;胡聿賢,1999;潘華等,2006)。統(tǒng)計值對地震資料樣本量需求的定量研究,國內(nèi)尚無專門的研究報道。在國外,Nava等(2017)利用蒙特卡羅模擬地震目錄進行抽樣估計,研究最大似然法計算值時對樣本量的需求,得出計算值時樣本量和精度之間的相互關(guān)聯(lián)關(guān)系,但沒有應(yīng)用實際地震目錄進行分析。國內(nèi)僅有少部分學(xué)者的研究涉及相關(guān)內(nèi)容,如韓曉明等(2016)研究河套地震帶值時空變化特征的文章中,探討了地震前后值的變化規(guī)律,分別使用最小二乘法和最大似然法對值時間和空間進行了掃描計算,在最小二乘法計算中,設(shè)定每次計算的窗長內(nèi)包含的樣本數(shù)目不少于100,而最大似然法的掃描窗內(nèi)包含的樣本數(shù)目不少于20;劉方斌等(2017)在魯西南聊考斷裂帶地震危險性評價與活動性分布的研究中,利用最大似然法和最小二乘法估算值并進行了對比研究,但是沒有給出樣本量的具體數(shù)目,且沒有用分震級段的方法進行值估算。

    本文采用Utsu(1965)提出的最大似然法和最小二乘法,利用模擬地震目錄和實際地震目錄,定量分析最小二乘法和最大似然法計算值時分別對地震資料樣本量的需求。

    1 資料概況

    1.1 模擬地震目錄

    蒙特卡羅法是以抽樣和隨機數(shù)的產(chǎn)生為基礎(chǔ)的隨機性方法,也稱為隨機抽樣法、計算機隨機模擬法等。蒙特卡羅方法的基本原理是通過數(shù)字模擬試驗,得到所要求解的出現(xiàn)某種事件的概率作為問題的近似解。其基本思想是:為了求解數(shù)學(xué)、物理、工程技術(shù)以及管理等方面的問題,首先建立概率模型或隨機過程,使用相應(yīng)的參數(shù),得到某些問題(如概率分布或數(shù)學(xué)期望等問題)的解;然后通過對模型或過程的觀察或抽樣試驗來計算所求參數(shù)的統(tǒng)計特征,并用算術(shù)平均值作為所求解的近似值。因此,只要盡可能完整地表示地震震源模型,即獲得某一區(qū)域內(nèi)地震事件的發(fā)生時間和空間上的分布規(guī)律,那么就可以直接利用蒙特卡羅方法產(chǎn)生合成地震序列。

    蒙特卡羅方法的基本依據(jù)是隨機模擬次數(shù)足夠多的情況下,事件發(fā)生的頻率可以反映事件發(fā)生的概率。當模擬次數(shù)有限時,計算結(jié)果和真實值之間必然存在誤差,這種誤差隨著模擬次數(shù)的增加而減少。

    蒙特卡羅模擬中隨機變量的簡單子樣(1,2,3,……,X)是獨立分布的,即每次模擬中事件發(fā)生的次數(shù)與其它任意1次中事件發(fā)生的次數(shù)無關(guān),那么隨機變量就是服從泊松分布的,當隨機變量的期望值趨向于無窮時,泊松分布趨近于正態(tài)分布。

    本文利用蒙特卡羅方法模擬地震目錄(張建中,1974a,1974b;任雪梅等,2011),研究不同樣本量下計算方法對值的影響。

    首先,構(gòu)造概率分布模型:

    其中,()為震級的概率密度函數(shù),表示為:

    其次,將公式(2)帶入公式(1),得到關(guān)于模擬震級的函數(shù):

    設(shè)定=1、0=2.0、u=8.0,利用MATLAB隨機數(shù)程序模擬生成不同樣本量的地震目錄,樣本數(shù)分別為5000、4500、4000、3500、3000、2500、2000、1750、1500、1250、1000、750、500、250、100、80、60、40、20、10,每個目錄模擬10000組。由于實際大地震非常稀缺,所以震級上限的限制對計算結(jié)果沒有顯著影響。

    1.2 汾渭地震帶實際地震目錄

    本研究使用的實際地震資料來自中國地震臺網(wǎng)中心的地震數(shù)據(jù)庫。研究表明地震資料的完整性對地震活動分析非常重要(Wiemer等,2000;Rotondi等,2002)。經(jīng)統(tǒng)計分析(焦遠碧等,1990;黃偉瓊等,1994a,1994b;周公威等,2007;謝卓娟等,2012a;徐偉進等,2014),1970年以來汾渭地震帶2.0以上的地震資料基本完整,故本研究選取1970年1月—2010年12月,≥2.0的8184次地震資料進行研究分析,其中2.0—2.9地震7044次、3.0—3.9地震964次、4.0—4.9地震148次、5.0—5.9地震25次、6.0—6.9地震3次,所用的地震目錄均已刪除前震、余震。汾渭地震帶4.7級以上地震的震中分布見圖1。

    圖1 汾渭地震帶范圍及震中分布

    2 模擬地震目錄的結(jié)果分析

    利用最小二乘法和最大似然法對模擬的地震目錄進行定量研究分析,比較2種方法在不同的樣本量下計算得到的平均值以及樣本量的大小對計算結(jié)果的影響。

    文中使用的最大似然法為Utsu(1965)提出的公式:

    采用上述最大似然法和最小二乘法計算,分別得到平均值隨樣本量的變化,如圖2、3所示。

    圖2 不同樣本量下最大似然法計算得到的平均b值擬合圖

    圖3 不同樣本量下最小二乘法計算得到的平均b值擬合圖

    由圖2、3可知,當樣本量大于300時,最大似然法計算的平均值能取得符合預(yù)期的數(shù)值;樣本量大于1000時,最小二乘法計算的平均值趨于穩(wěn)定,向理論值收斂。

    表1 最大似然法模擬結(jié)果

    表2 最小二乘法模擬結(jié)果

    當設(shè)定=1時,不同樣本下得到的值分布直方圖見圖4。由圖可知,當樣本量≥200時,最大似然法計算的值很好地向設(shè)定值收斂;樣本量≥500時,最小二乘法計算的值有較好的收斂效果。從估值概率來看,取得的值一般都偏低。

    根據(jù)上述分析,2種方法雖然都能得到符合預(yù)期的數(shù)值,但最小二乘法計算的平均值通常比最大似然法的小,而標準差比最大似然法的大;樣本量小于1000時,利用最小二乘法得到的平均值精確度較低;樣本量大于1000時,2種計算方法得出的結(jié)果差別不大。進一步分析認為,雖然2種方法都需要一定的樣本量,且樣本量越大、得到的結(jié)果越準確,但最大似然法對樣本量的要求要比最小二乘法低,樣本量大于200時,計算得到的平均值與設(shè)定值一致性較好;樣本量大于1000時,基本等于設(shè)定值。

    3 實際地震目錄的計算結(jié)果分析

    根據(jù)汾渭地震帶地震資料的完整性和可靠性研究,自1970年以來,臺站記錄得到的≥2.0地震目錄基本完整,故本研究以汾渭地震帶1970年1月—2010年12月的地震記錄為計算樣本,研究實際地震目錄樣本量對值計算的影響。表3給出該地震帶不同震級檔的地震數(shù)目,并計算出不同震級檔的年平均發(fā)生率(黃瑋瓊等,1989;潘華等,2006;吳兆營等2005;謝卓娟等,2012b)。

    圖4 不同樣本量的b值分布直方圖

    表3 汾渭地震帶地震分檔統(tǒng)計和年平均發(fā)生率

    續(xù)表

    震級檔M地震個數(shù)年平均發(fā)生率 5.0—5.4150.36 5.5—5.9100.23 6.0—8.530.07

    注:年平均發(fā)生率指震級≥的年均地震數(shù),代表地震活動水平。

    根據(jù)表3中統(tǒng)計的各震級檔的地震個數(shù)和年平均發(fā)生率,利用最小二乘法和最大似然法計算得到的值分別為0.75和0.75335。

    為研究實際地震目錄下值計算對樣本量的要求,采用任雪梅(2011)的抽樣原則,對汾渭地震帶的地震目錄進行均勻抽樣分析,每次抽樣完將樣本放回進行下一次抽樣,得到樣本量分別為10、50、100、200、300、500、700、1000、5000的地震目錄,每個目錄重復(fù)抽樣10000次。同樣,采用2種方法分別計算不同樣本量下的平均值(表4),并分析不同樣本量下的值變化。

    表4 汾渭地震帶b值擬合情況(1500—2010年)

    由表4可以看出,2種計算方法得到值隨著樣本量的增加逐漸接近真實值;當樣本量大于200時,最大似然法能夠得出相對穩(wěn)定可靠的值;當樣本量大于500時,最小二乘法才能得到相對穩(wěn)定可靠的值。因此,對實際的地震目錄,最小二乘法對樣本量的要求也比較高,且受地震目錄完整性的影響較大。

    4 結(jié)論

    通過上述研究和對比分析,得到以下認識:

    (1)利用最大似然法計算值時,樣本量至少要在200以上;對于最小二乘法,樣本量要求不少于1000;一般來說,計算的值都低于設(shè)定值。

    (2)最小二乘法利用震級-頻度的線性關(guān)系來擬合計算值,樣本量不足會影響其線性關(guān)系。因此,最小二乘法受樣本量影響較大,樣本量小于1000時,值的計算值與設(shè)定值相差較大,數(shù)值也不穩(wěn)定。

    (3)最大似然法方便快捷,受樣本量影響小,計算出的值相對穩(wěn)定,但誤差估計值偏大,其利用平均震級計算,只與地震個數(shù)有關(guān),受數(shù)據(jù)質(zhì)量的影響較?。蛔钚《朔ㄊ軜颖玖坑绊戄^大,在樣本量充足的情況下計算的值比較準確,而在樣本量不足時值波動較大,誤差也隨著樣本量的減少而增加。在實際應(yīng)用中,以地震目錄充足為前提,研究不同區(qū)域的值優(yōu)先選用最小二乘法;而研究某區(qū)域值時間上的變化時,采用最大似然法估算的值相對穩(wěn)定,更能體現(xiàn)值長期變化的趨勢。

    (4)從誤差和計算量來看,最大似然法比最小二乘法要小,但隨著樣本量的增加,2種方法計算結(jié)果的差異越來越小。但由于半對數(shù)坐標下不同震級檔數(shù)據(jù)權(quán)重的不對等性,國外已極少使用最小二乘法。

    韓渭賓,2003.值在地震預(yù)測中的三類應(yīng)用及其物理基礎(chǔ)與須注意的問題.四川地震,(1):1—5.

    韓曉明,張文韜,王樹波等,2016.河套地震帶的值時空變化特征分析.中國地震,32(3):522—532.

    胡聿賢,1999.地震安全性評價技術(shù)教程.北京:地震出版社.

    黃瑋瓊,時振梁,曹學(xué)鋒,1989.值統(tǒng)計中的影響因素及危險性分析中值的選取.地震學(xué)報,11(4):351—361.

    黃瑋瓊,李文香,曹學(xué)鋒,1994a.中國大陸地震資料完整性研究之一——以華北地區(qū)為例.地震學(xué)報,16(3):273—280.

    黃瑋瓊,李文香,曹學(xué)鋒,1994b.中國大陸地震資料完整性研究之二——分區(qū)地震資料基本完整的起始年分布圖象.地震學(xué)報,16(4):423—432.

    黃瑋瓊,李文香,1998.地震區(qū)劃中值統(tǒng)計時空范圍的確定.地震學(xué)報,20(5):2—6.

    焦遠碧,吳開統(tǒng),楊滿棟,1990.我國地震臺網(wǎng)監(jiān)測能力及臺網(wǎng)觀測條件質(zhì)量評定.中國地震,6(4):3—9.

    劉方斌,曲均浩,田兆陽等,2017.魯西南聊考斷裂帶地震危險性評價與活動性分布.大地測量與地球動力學(xué),37(8):802—807.

    劉靜偉,呂悅軍,2016.川滇地區(qū)值空間分布特征及其與震源類型關(guān)系的初步探討.震災(zāi)防御技術(shù),11(3):561—572.

    潘華,李金臣,2006.地震統(tǒng)計區(qū)地震活動性參數(shù)值及4不確定性研究.震災(zāi)防御技術(shù),1(3):218—224.

    任雪梅,2011.地震區(qū)劃中值統(tǒng)計的若干問題研究.北京:中國地震局地球物理研究所.

    任雪梅,高孟潭,馮靜,2011.地震目錄的完整性對值計算的影響.震災(zāi)防御技術(shù),6(3):257—268.

    沈建文,余湛,邱瑛,2007.地震安評中地震活動性的統(tǒng)計區(qū)域與值.國際地震動態(tài),(3):1—6.

    王熠熙,張輝,劉雙慶等,2015.河北平原地震帶值時空變化特征.地震工程學(xué)報,37(01):188—195.

    吳兆營,薄景山,劉志平等,2005.東北地震區(qū)值和地震年平均發(fā)生率的統(tǒng)計分析.東北地震研究,21(3):27—32,72.

    謝卓娟,呂悅軍,張力方,等,2012a.基于現(xiàn)代地震資料確定汾渭地震帶分區(qū)及其地震活動性參數(shù).地球物理學(xué)進展,27(3):894—902.

    謝卓娟,呂悅軍,彭艷菊等,2012b.東北地震區(qū)小震資料完整性分析及其對地震活動性參數(shù)的影響研究.中國地震,28(3):256—265.

    謝卓娟,呂悅軍,蘭景巖等,2013.值和4的統(tǒng)計分析及其不確定性對地震危險性分析結(jié)果的影響研究.地震研究,36(1):86—92.

    謝卓娟,李山有,呂悅軍,2015.滇西南地區(qū)主要活動斷裂的值空間分布特征.地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報),40(10):1755—1766.

    徐偉進,高孟潭,2014.中國大陸及周緣地震目錄完整性統(tǒng)計分析.地球物理學(xué)報,57(9):2802—2812.

    鄢家全,韓煒,高孟潭,1996.地震活動性參數(shù)的不確定性及其對區(qū)劃結(jié)果的影響.中國地震,12(S1):71—77.

    張廣偉,2016.云南地區(qū)地震的重新定位及值研究.中國地震,32(1):54—62.

    張建中,1974a.蒙特卡洛方法(Ⅰ).數(shù)學(xué)的實踐與認識,(01):28—40.

    張建中,1974b.蒙特卡洛方法(Ⅱ).數(shù)學(xué)的實踐與認識,(02):43—56.

    周公威,張伯明,吳忠良等,2007.中國數(shù)字地震臺網(wǎng)(CDSN)的近期發(fā)展.地球物理學(xué)進展,22(4):1130—1134.

    Maximum likelihood estimate ofin the formula log=-and its confidence limits. Bulletin of the Earthquake Research Institute, 43(2): 237—239.

    Gulia L., Wiemer S., 2010. The influence of tectonic regimes on the earthquake size distribution: A case study for Italy. Geophysical Research Letters, 37(10): L10305.

    Nava F. A., Márquez-Ramírez V. H., Zú?iga F. R., et al., 2017. Gutenberg-Richter-value maximum likelihood estimation and sample size. Journal of Seismology, 21(1): 127—135.

    Rundle J. B., 1989. Derivation of the complete Gutenberg-Richter magnitude-frequency relation using the principle of scale invariance. Journal of Geophysical Research: Solid Earth, 94(B9): 12337—12342.

    Rotondi R., Garavaglia E., 2002. Statistical Analysis of the Completeness of a Seismic Catalogue. Natural Hazards, 25(3):245—258.

    Schorlemmer D., Wiemer S., Wyss M., 2005. Variations in earthquake-size distribution across different stress regimes. Nature, 437(7058): 538—542.

    Utsu T., 1965. A method for determining the value ofin a formula log=-showing the magnitude-frequency relation for earthquakes. Geophysical Bulletin of Hokkaido University, (13): 99—103.

    Wiemer S., Wyss M., 2000. Minimum magnitude of completeness in earthquake catalogs: examples from Alaska, the Western United States, and Japan. Bulletin of the Seismological Society of America, 90(4): 859—869.

    The Study of Sample Size on-value Statistics in the Gutenberg-Richter's Law

    Li Shijie, Lü Yuejun and Liu Jingwei

    (Institute of Crustal Dynamics, China Earthquake Administration, Beijing 100085, China)

    The-value is an important indicator for evaluating the level of seismicity, and is widely used in seismic hazard analysis and earthquake prediction research. The-value is basically effected by some factors, such as the actual data integrity, the quantity of earthquake samples, calculation methods and so on. The most widely used two methods for calculatingvalues include the least square method (LSM) and the maximum likelihood method (MLM), and the results from both methods can remarkably vary with different the quantity of earthquake samples. In this paper, based on the real catalog of the Fen-Wei seismic zone and the Monte-Carlo simulated earthquake catalog, we use different number of samples to calculate thevalues and try to find the differences between the calculated results and the set values/true values. It turns out that the threshold value of samples is 1000 and 200 for LSM and MLM, respectively. When the sample size does not meet the requirements, the calculated-value is not reliable. It suggets that LSM is suitable for calculating the-value of the whole region, while MLM is more advantageous in studying the-value of a region on the time axis. This study provides a reference for the threshold values of the samples for two methods, which is of great significance to the seismic hazard analysis and earthquake prediction.

    -value; The least square method; The maximum likelihood method; Sample size

    李世杰,呂悅軍,劉靜偉,2018.古登堡-里希特定律中的值統(tǒng)計樣本量研究.震災(zāi)防御技術(shù),13(3):636—645.

    10.11899/zzfy20180315

    中國地震局地殼應(yīng)力研究所中央級公益性科研院所基本科研業(yè)務(wù)專項(ZDJ2016-03)

    2018-02-28

    李世杰,男,生于1991年。碩士研究生。主要從事工程地震研究。E-mail:729883299@qq.com

    呂悅軍,男,生于1966年。研究員。研究領(lǐng)域:工程地震。E-mail:luyj1@263.net

    猜你喜歡
    震級樣本量乘法
    算乘法
    基于累積絕對位移值的震級估算方法
    醫(yī)學(xué)研究中樣本量的選擇
    我們一起來學(xué)習(xí)“乘法的初步認識”
    地震后各國發(fā)布的震級可能不一樣?
    《整式的乘法與因式分解》鞏固練習(xí)
    新震級國家標準在大同臺的應(yīng)用與評估
    山西地震(2020年1期)2020-04-08 07:34:26
    把加法變成乘法
    航空裝備測試性試驗樣本量確定方法
    Sample Size Calculations for Comparing Groups with Binary Outcomes
    一区二区三区精品91| 色婷婷av一区二区三区视频| 在线观看免费视频网站a站| 咕卡用的链子| 国产成人av教育| 99国产综合亚洲精品| 国产亚洲精品第一综合不卡| 国产伦一二天堂av在线观看| 一个人免费在线观看的高清视频| 国产精品亚洲一级av第二区| √禁漫天堂资源中文www| 午夜福利免费观看在线| 黄色毛片三级朝国网站| 久久午夜亚洲精品久久| 天堂动漫精品| 国产野战对白在线观看| 91麻豆精品激情在线观看国产 | 中文字幕精品免费在线观看视频| av片东京热男人的天堂| 精品福利永久在线观看| 久久精品国产清高在天天线| 久久精品亚洲精品国产色婷小说| 国产熟女午夜一区二区三区| 视频区欧美日本亚洲| 久久人人97超碰香蕉20202| 亚洲欧美一区二区三区久久| 亚洲成国产人片在线观看| 日本免费a在线| 亚洲成人精品中文字幕电影 | 在线观看免费午夜福利视频| 十八禁网站免费在线| 亚洲专区国产一区二区| 欧美日韩亚洲国产一区二区在线观看| 久久香蕉国产精品| 亚洲一区二区三区色噜噜 | 欧美日韩一级在线毛片| av欧美777| 一边摸一边抽搐一进一出视频| www.精华液| 大香蕉久久成人网| 美女高潮到喷水免费观看| 色在线成人网| 一本综合久久免费| 精品第一国产精品| 久久精品国产99精品国产亚洲性色 | 亚洲精品一二三| 在线观看www视频免费| 亚洲精品av麻豆狂野| 久久久精品国产亚洲av高清涩受| 香蕉丝袜av| 日韩欧美三级三区| 神马国产精品三级电影在线观看 | 97超级碰碰碰精品色视频在线观看| 久久久精品欧美日韩精品| av视频免费观看在线观看| 成人特级黄色片久久久久久久| 两个人看的免费小视频| 别揉我奶头~嗯~啊~动态视频| 亚洲性夜色夜夜综合| 嫩草影视91久久| 精品福利观看| 久久 成人 亚洲| 久久久国产成人免费| 精品少妇一区二区三区视频日本电影| 亚洲av电影在线进入| 美女国产高潮福利片在线看| 亚洲七黄色美女视频| 日本免费一区二区三区高清不卡 | 精品午夜福利视频在线观看一区| 国产精品 国内视频| 日韩三级视频一区二区三区| 国产三级在线视频| 水蜜桃什么品种好| xxxhd国产人妻xxx| 大码成人一级视频| 国产国语露脸激情在线看| 天堂中文最新版在线下载| 国产精品久久久人人做人人爽| 一进一出抽搐gif免费好疼 | 国产蜜桃级精品一区二区三区| 国产深夜福利视频在线观看| 啦啦啦在线免费观看视频4| 日本撒尿小便嘘嘘汇集6| 夜夜爽天天搞| 亚洲性夜色夜夜综合| 国产亚洲欧美在线一区二区| 精品国产超薄肉色丝袜足j| 国产精品国产高清国产av| 中文字幕av电影在线播放| 夫妻午夜视频| 欧美黄色淫秽网站| 久久性视频一级片| 啪啪无遮挡十八禁网站| 久久久精品欧美日韩精品| 天堂动漫精品| 亚洲国产看品久久| 纯流量卡能插随身wifi吗| 在线十欧美十亚洲十日本专区| 夜夜躁狠狠躁天天躁| 黑人猛操日本美女一级片| 日韩有码中文字幕| 精品福利永久在线观看| 黄频高清免费视频| 国产黄a三级三级三级人| 国产精品国产av在线观看| 黑人巨大精品欧美一区二区蜜桃| 又紧又爽又黄一区二区| 午夜91福利影院| xxxhd国产人妻xxx| 黑人欧美特级aaaaaa片| 高清av免费在线| 亚洲色图 男人天堂 中文字幕| 中国美女看黄片| 亚洲欧美一区二区三区黑人| 亚洲性夜色夜夜综合| 国产男靠女视频免费网站| 水蜜桃什么品种好| 亚洲在线自拍视频| 日韩一卡2卡3卡4卡2021年| 午夜福利在线观看吧| 99精品欧美一区二区三区四区| 老汉色∧v一级毛片| 欧美久久黑人一区二区| 国产熟女午夜一区二区三区| 满18在线观看网站| 国产精品九九99| 日韩 欧美 亚洲 中文字幕| 激情视频va一区二区三区| 好看av亚洲va欧美ⅴa在| 国产成人啪精品午夜网站| 后天国语完整版免费观看| 电影成人av| 黄色视频,在线免费观看| 91麻豆精品激情在线观看国产 | 精品国产乱码久久久久久男人| 国产欧美日韩一区二区三| 嫩草影视91久久| 国产精品99久久99久久久不卡| 老司机亚洲免费影院| 国产激情欧美一区二区| 无人区码免费观看不卡| 国产精品久久久久久人妻精品电影| 免费高清在线观看日韩| 中文字幕人妻丝袜一区二区| 国产精品久久久av美女十八| 91成人精品电影| 亚洲成人久久性| 丰满饥渴人妻一区二区三| 日韩中文字幕欧美一区二区| 欧美成人午夜精品| 欧美精品亚洲一区二区| 美女国产高潮福利片在线看| 不卡av一区二区三区| 午夜91福利影院| 久久青草综合色| 可以免费在线观看a视频的电影网站| 久久久精品国产亚洲av高清涩受| 18禁裸乳无遮挡免费网站照片 | 99在线人妻在线中文字幕| 一区二区日韩欧美中文字幕| 巨乳人妻的诱惑在线观看| 久久人人精品亚洲av| 村上凉子中文字幕在线| 国产精品影院久久| 亚洲国产欧美一区二区综合| www.自偷自拍.com| 午夜精品国产一区二区电影| 久热爱精品视频在线9| 美女国产高潮福利片在线看| 女性被躁到高潮视频| 亚洲七黄色美女视频| 国产有黄有色有爽视频| 国产精品日韩av在线免费观看 | 老汉色∧v一级毛片| 变态另类成人亚洲欧美熟女 | 操美女的视频在线观看| 国产熟女午夜一区二区三区| 啪啪无遮挡十八禁网站| 免费观看人在逋| 亚洲精品粉嫩美女一区| 啦啦啦 在线观看视频| 视频区欧美日本亚洲| 一个人免费在线观看的高清视频| 久久香蕉国产精品| 免费高清在线观看日韩| 在线看a的网站| 69av精品久久久久久| 国产一区二区三区在线臀色熟女 | 成年人黄色毛片网站| 1024视频免费在线观看| av超薄肉色丝袜交足视频| 免费观看人在逋| 免费高清视频大片| 久久天堂一区二区三区四区| 国产精品久久久人人做人人爽| 波多野结衣高清无吗| 国产一区二区三区视频了| 午夜福利影视在线免费观看| 久久久久久久久中文| 丰满迷人的少妇在线观看| 国产aⅴ精品一区二区三区波| √禁漫天堂资源中文www| 自拍欧美九色日韩亚洲蝌蚪91| 久久这里只有精品19| 99久久国产精品久久久| 欧美日韩视频精品一区| 这个男人来自地球电影免费观看| 久久九九热精品免费| 国产精品偷伦视频观看了| 在线观看日韩欧美| aaaaa片日本免费| 成年女人毛片免费观看观看9| 成年女人毛片免费观看观看9| 精品久久久久久成人av| 女人被狂操c到高潮| 亚洲精品国产色婷婷电影| 成人18禁高潮啪啪吃奶动态图| 级片在线观看| 精品久久久久久久久久免费视频 | 美女福利国产在线| 国产精品乱码一区二三区的特点 | 亚洲专区字幕在线| 亚洲专区字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 一边摸一边做爽爽视频免费| 亚洲成人久久性| 波多野结衣av一区二区av| 国产精品国产av在线观看| 亚洲久久久国产精品| 少妇裸体淫交视频免费看高清 | 欧美日韩精品网址| 麻豆成人av在线观看| 国产99白浆流出| 黑人巨大精品欧美一区二区蜜桃| 日韩成人在线观看一区二区三区| 1024香蕉在线观看| 国产精品日韩av在线免费观看 | 50天的宝宝边吃奶边哭怎么回事| 一区在线观看完整版| 中文欧美无线码| 久久精品成人免费网站| a级毛片在线看网站| 后天国语完整版免费观看| bbb黄色大片| x7x7x7水蜜桃| 中文字幕色久视频| 久久久国产一区二区| 色尼玛亚洲综合影院| 久久伊人香网站| 国产单亲对白刺激| 欧美一级毛片孕妇| 美女扒开内裤让男人捅视频| 久久国产精品人妻蜜桃| 一级黄色大片毛片| 黄色 视频免费看| 在线观看www视频免费| 久久婷婷成人综合色麻豆| 亚洲专区中文字幕在线| av中文乱码字幕在线| www.熟女人妻精品国产| 成熟少妇高潮喷水视频| √禁漫天堂资源中文www| 精品久久久久久电影网| 亚洲国产欧美一区二区综合| 母亲3免费完整高清在线观看| 在线天堂中文资源库| 国产亚洲欧美精品永久| 一本大道久久a久久精品| 大型黄色视频在线免费观看| 亚洲专区中文字幕在线| 欧美乱色亚洲激情| 老司机深夜福利视频在线观看| 在线观看66精品国产| 成人国产一区最新在线观看| 日韩 欧美 亚洲 中文字幕| 精品乱码久久久久久99久播| 国产精品爽爽va在线观看网站 | 午夜亚洲福利在线播放| 99国产精品一区二区蜜桃av| av福利片在线| www.熟女人妻精品国产| 国产亚洲精品一区二区www| 国产精品亚洲一级av第二区| 日韩欧美国产一区二区入口| 精品一品国产午夜福利视频| 天堂√8在线中文| 国内毛片毛片毛片毛片毛片| 久久久久九九精品影院| 国产区一区二久久| 最新美女视频免费是黄的| 国产成人精品久久二区二区免费| 久久久久国产精品人妻aⅴ院| 亚洲专区中文字幕在线| 日韩欧美三级三区| 精品免费久久久久久久清纯| 久9热在线精品视频| 美女扒开内裤让男人捅视频| 我的亚洲天堂| 亚洲欧美精品综合一区二区三区| 嫩草影视91久久| 免费av毛片视频| 一个人观看的视频www高清免费观看 | 欧美+亚洲+日韩+国产| 国产高清激情床上av| 黑丝袜美女国产一区| 男人的好看免费观看在线视频 | 欧美日韩乱码在线| 91在线观看av| 国产成人影院久久av| 黄色片一级片一级黄色片| 午夜成年电影在线免费观看| av国产精品久久久久影院| 色婷婷久久久亚洲欧美| 精品少妇一区二区三区视频日本电影| 日本五十路高清| 岛国视频午夜一区免费看| 精品福利永久在线观看| 精品一区二区三区四区五区乱码| 午夜成年电影在线免费观看| 好看av亚洲va欧美ⅴa在| 色哟哟哟哟哟哟| 一a级毛片在线观看| 久久精品人人爽人人爽视色| 亚洲色图 男人天堂 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 不卡av一区二区三区| 亚洲精品在线美女| 久久精品国产亚洲av香蕉五月| 女性生殖器流出的白浆| 高清毛片免费观看视频网站 | 久久伊人香网站| 黄色怎么调成土黄色| 日韩中文字幕欧美一区二区| 巨乳人妻的诱惑在线观看| 51午夜福利影视在线观看| 国产区一区二久久| 黑人操中国人逼视频| 悠悠久久av| 国产成年人精品一区二区 | 香蕉久久夜色| 嫩草影视91久久| 亚洲第一青青草原| 在线观看午夜福利视频| 欧美中文日本在线观看视频| 9191精品国产免费久久| 91字幕亚洲| 免费在线观看完整版高清| 大香蕉久久成人网| 国产高清videossex| 曰老女人黄片| 波多野结衣av一区二区av| 久久香蕉国产精品| 日韩有码中文字幕| 精品日产1卡2卡| 久久精品国产综合久久久| 亚洲av成人不卡在线观看播放网| 欧美激情极品国产一区二区三区| 成人黄色视频免费在线看| 人人澡人人妻人| 久久久久久亚洲精品国产蜜桃av| 99精品在免费线老司机午夜| 久久久久久人人人人人| 男女床上黄色一级片免费看| 丝袜人妻中文字幕| 老司机福利观看| 欧美性长视频在线观看| av中文乱码字幕在线| 深夜精品福利| 久久久久国内视频| 高清在线国产一区| 亚洲欧美日韩另类电影网站| 日本 av在线| 国产精品 国内视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品亚洲熟妇少妇任你| 亚洲国产精品sss在线观看 | 人成视频在线观看免费观看| 成人国语在线视频| 咕卡用的链子| 韩国精品一区二区三区| 啦啦啦 在线观看视频| tocl精华| 我的亚洲天堂| 国产亚洲精品久久久久久毛片| 日韩欧美一区二区三区在线观看| 18禁美女被吸乳视频| 首页视频小说图片口味搜索| 国产熟女xx| 91成人精品电影| 丝袜美腿诱惑在线| 亚洲成国产人片在线观看| 久久精品亚洲熟妇少妇任你| 国产成人系列免费观看| 黄频高清免费视频| 夜夜夜夜夜久久久久| 精品久久久精品久久久| 久久久久亚洲av毛片大全| 满18在线观看网站| 国产亚洲欧美98| 香蕉久久夜色| 欧美日本亚洲视频在线播放| 久久人妻福利社区极品人妻图片| 亚洲成国产人片在线观看| 亚洲国产精品合色在线| 精品一区二区三卡| 久久久久久久久久久久大奶| 多毛熟女@视频| 欧美成人午夜精品| 欧美老熟妇乱子伦牲交| 一区二区三区精品91| 免费在线观看日本一区| 久久欧美精品欧美久久欧美| 日本免费a在线| 高清在线国产一区| 亚洲成人免费电影在线观看| 黄片小视频在线播放| 日韩高清综合在线| 午夜免费鲁丝| 自线自在国产av| 精品国产乱码久久久久久男人| 美女扒开内裤让男人捅视频| 亚洲精品成人av观看孕妇| 亚洲av美国av| 精品一区二区三区四区五区乱码| 国产免费男女视频| 亚洲片人在线观看| 桃红色精品国产亚洲av| 丁香六月欧美| 波多野结衣av一区二区av| 色精品久久人妻99蜜桃| 十八禁国产超污无遮挡网站| 悠悠久久av| 日韩欧美免费精品| 99久国产av精品| 国产精品98久久久久久宅男小说| 色吧在线观看| 乱人视频在线观看| 免费在线观看成人毛片| av天堂中文字幕网| 亚洲内射少妇av| 亚洲不卡免费看| 亚洲国产高清在线一区二区三| 亚洲av成人不卡在线观看播放网| 久久国产精品人妻蜜桃| 一进一出抽搐gif免费好疼| 观看美女的网站| 日韩欧美国产在线观看| 日本免费一区二区三区高清不卡| 俄罗斯特黄特色一大片| 日本在线视频免费播放| 国产精品综合久久久久久久免费| 少妇裸体淫交视频免费看高清| 日韩欧美在线二视频| 免费在线观看影片大全网站| 乱人视频在线观看| 精品久久久久久成人av| 免费看光身美女| 亚州av有码| 97超级碰碰碰精品色视频在线观看| 欧美绝顶高潮抽搐喷水| 99久久精品热视频| 欧美色视频一区免费| 色5月婷婷丁香| 色吧在线观看| 免费在线观看亚洲国产| 欧美最新免费一区二区三区 | 嫩草影院新地址| 天堂影院成人在线观看| 欧美日本亚洲视频在线播放| 午夜老司机福利剧场| 国产免费男女视频| 国产精品三级大全| 校园春色视频在线观看| 国产高清视频在线播放一区| 日本a在线网址| 国产 一区 欧美 日韩| 好男人电影高清在线观看| 老鸭窝网址在线观看| 波多野结衣巨乳人妻| 国产主播在线观看一区二区| 亚洲一区二区三区色噜噜| 亚洲av第一区精品v没综合| 亚洲欧美日韩高清在线视频| 亚洲乱码一区二区免费版| 禁无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 免费高清视频大片| 欧美中文日本在线观看视频| 国产精品一区二区三区四区免费观看 | 小说图片视频综合网站| 人人妻人人看人人澡| 看十八女毛片水多多多| 色哟哟·www| 色av中文字幕| 两个人的视频大全免费| 国产伦人伦偷精品视频| 国产又黄又爽又无遮挡在线| 午夜福利18| 给我免费播放毛片高清在线观看| 国产亚洲精品av在线| 国产欧美日韩一区二区三| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 美女xxoo啪啪120秒动态图 | 精品人妻偷拍中文字幕| 亚洲av第一区精品v没综合| 脱女人内裤的视频| av欧美777| 18禁黄网站禁片午夜丰满| 在线观看66精品国产| 我要看日韩黄色一级片| 亚洲av成人精品一区久久| 18禁黄网站禁片免费观看直播| 国产精品嫩草影院av在线观看 | 国产欧美日韩一区二区精品| 757午夜福利合集在线观看| 亚洲精品色激情综合| 亚洲欧美精品综合久久99| 亚洲 国产 在线| 99riav亚洲国产免费| 好看av亚洲va欧美ⅴa在| 国产欧美日韩精品亚洲av| 久久精品国产亚洲av天美| 日本与韩国留学比较| 久久久国产成人免费| 一级黄色大片毛片| 亚洲七黄色美女视频| 可以在线观看毛片的网站| 丰满的人妻完整版| 国内精品久久久久久久电影| 我的老师免费观看完整版| 亚洲精华国产精华精| 老司机深夜福利视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 一级黄片播放器| 国产精品久久久久久久电影| 国产精品久久久久久人妻精品电影| 亚洲成av人片在线播放无| 中文字幕人成人乱码亚洲影| 一本一本综合久久| 老女人水多毛片| 国语自产精品视频在线第100页| 久久国产乱子免费精品| netflix在线观看网站| 国产精品永久免费网站| 亚洲人成伊人成综合网2020| 熟女人妻精品中文字幕| 色综合欧美亚洲国产小说| 中文字幕av成人在线电影| 成人无遮挡网站| 国产精品乱码一区二三区的特点| 国产欧美日韩精品亚洲av| 国内久久婷婷六月综合欲色啪| 亚洲自拍偷在线| 一级作爱视频免费观看| 国产老妇女一区| 日日干狠狠操夜夜爽| 91在线观看av| www.999成人在线观看| 校园春色视频在线观看| 欧美又色又爽又黄视频| 精品人妻熟女av久视频| 成人毛片a级毛片在线播放| 午夜福利视频1000在线观看| 欧美日韩瑟瑟在线播放| 日韩欧美免费精品| 看免费av毛片| 美女高潮喷水抽搐中文字幕| 精品99又大又爽又粗少妇毛片 | 国产极品精品免费视频能看的| 一个人观看的视频www高清免费观看| 精品人妻一区二区三区麻豆 | 制服丝袜大香蕉在线| 午夜视频国产福利| 国产亚洲欧美在线一区二区| 一本精品99久久精品77| 丰满人妻熟妇乱又伦精品不卡| 日韩中文字幕欧美一区二区| 午夜老司机福利剧场| 看十八女毛片水多多多| 超碰av人人做人人爽久久| 怎么达到女性高潮| 91麻豆精品激情在线观看国产| 亚洲avbb在线观看| 亚洲精品日韩av片在线观看| 国产精品伦人一区二区| 男人舔奶头视频| 国产成人av教育| 国产极品精品免费视频能看的| 久久久久国内视频| 99视频精品全部免费 在线| 成年人黄色毛片网站| 日本黄色视频三级网站网址| 国产精品爽爽va在线观看网站| 国产综合懂色| 国产91精品成人一区二区三区| www.熟女人妻精品国产| 国产精品久久久久久精品电影| 久久精品国产亚洲av天美| 99久久精品热视频| 国产av一区在线观看免费| 国产精品久久久久久精品电影| 校园春色视频在线观看| 少妇人妻一区二区三区视频| 91麻豆av在线| 搡老岳熟女国产| 婷婷丁香在线五月| av国产免费在线观看| x7x7x7水蜜桃| 亚洲乱码一区二区免费版| 麻豆av噜噜一区二区三区| 亚洲美女黄片视频| 国产69精品久久久久777片| 亚洲成人精品中文字幕电影| 又爽又黄无遮挡网站|