• <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
    免费观看精品视频网站| 亚洲av中文字字幕乱码综合| 亚洲最大成人中文| 亚洲精品成人久久久久久| 国产视频首页在线观看| 热99re8久久精品国产| 男女那种视频在线观看| 国产老妇女一区| 亚洲成色77777| 波多野结衣高清无吗| 中文字幕av成人在线电影| 亚洲国产精品专区欧美| 亚洲国产最新在线播放| 久久人妻av系列| 亚洲成色77777| 黄片wwwwww| 久久综合国产亚洲精品| 久久人人爽人人片av| 非洲黑人性xxxx精品又粗又长| 在线天堂最新版资源| 99久久九九国产精品国产免费| 在线播放国产精品三级| 精品久久久久久成人av| 免费在线观看成人毛片| 蜜桃亚洲精品一区二区三区| 日韩,欧美,国产一区二区三区 | 超碰av人人做人人爽久久| 三级毛片av免费| 99热精品在线国产| 亚洲精品自拍成人| 少妇裸体淫交视频免费看高清| 国产精品美女特级片免费视频播放器| 国产伦精品一区二区三区四那| 人妻夜夜爽99麻豆av| 爱豆传媒免费全集在线观看| 国产亚洲精品久久久com| 尤物成人国产欧美一区二区三区| 国产精品久久久久久久电影| 国产av不卡久久| 禁无遮挡网站| 精品无人区乱码1区二区| 18禁裸乳无遮挡免费网站照片| 淫秽高清视频在线观看| 丰满人妻一区二区三区视频av| 日韩精品有码人妻一区| 久久亚洲精品不卡| 日韩三级伦理在线观看| 国产亚洲午夜精品一区二区久久 | 看免费成人av毛片| 三级男女做爰猛烈吃奶摸视频| 日韩欧美 国产精品| 男人舔奶头视频| 免费播放大片免费观看视频在线观看 | 国产精品综合久久久久久久免费| 波多野结衣高清无吗| 国产视频首页在线观看| 老司机福利观看| 91狼人影院| 只有这里有精品99| 69人妻影院| 18禁动态无遮挡网站| 黄色配什么色好看| 久久精品久久久久久噜噜老黄 | 三级男女做爰猛烈吃奶摸视频| 国产av在哪里看| 中文乱码字字幕精品一区二区三区 | 亚洲国产日韩欧美精品在线观看| 欧美高清成人免费视频www| 高清视频免费观看一区二区 | 91久久精品电影网| av福利片在线观看| 国内精品一区二区在线观看| 美女脱内裤让男人舔精品视频| 成人毛片a级毛片在线播放| 村上凉子中文字幕在线| 色综合亚洲欧美另类图片| 欧美精品国产亚洲| 女人十人毛片免费观看3o分钟| 国产在视频线精品| 亚洲av中文av极速乱| 国产探花极品一区二区| 99热精品在线国产| 亚洲人成网站在线观看播放| 一级毛片aaaaaa免费看小| 国产高清三级在线| 伦精品一区二区三区| 国产老妇女一区| 欧美三级亚洲精品| 2021天堂中文幕一二区在线观| 久久这里只有精品中国| 亚洲自偷自拍三级| 免费看美女性在线毛片视频| 一区二区三区乱码不卡18| 又爽又黄无遮挡网站| 大又大粗又爽又黄少妇毛片口| 精品酒店卫生间| 久久婷婷人人爽人人干人人爱| 插阴视频在线观看视频| 欧美变态另类bdsm刘玥| 成人三级黄色视频| 久久久久九九精品影院| 小蜜桃在线观看免费完整版高清| 久久精品国产自在天天线| 舔av片在线| 最后的刺客免费高清国语| 男人和女人高潮做爰伦理| 亚洲av熟女| 国产精品久久视频播放| 久久人妻av系列| 天堂中文最新版在线下载 | 免费人成在线观看视频色| 日韩av在线免费看完整版不卡| 日韩强制内射视频| 国产乱人偷精品视频| 你懂的网址亚洲精品在线观看 | 麻豆成人av视频| 精华霜和精华液先用哪个| 日韩av在线大香蕉| 国产色婷婷99| 亚洲精品日韩av片在线观看| 国产精华一区二区三区| 中国美白少妇内射xxxbb| 亚洲国产精品sss在线观看| 卡戴珊不雅视频在线播放| 午夜久久久久精精品| 亚洲中文字幕一区二区三区有码在线看| 色哟哟·www| 青春草国产在线视频| 国产成人freesex在线| 婷婷六月久久综合丁香| 欧美xxxx性猛交bbbb| 亚洲性久久影院| 婷婷色综合大香蕉| 国产亚洲av嫩草精品影院| 国产激情偷乱视频一区二区| 搡女人真爽免费视频火全软件| 大香蕉97超碰在线| 国产成人freesex在线| 国产成人91sexporn| 麻豆成人午夜福利视频| 尾随美女入室| av免费在线看不卡| 国产精品国产三级国产专区5o | 国产淫片久久久久久久久| 日韩欧美在线乱码| 国产精品久久久久久精品电影| 好男人在线观看高清免费视频| 小蜜桃在线观看免费完整版高清| av免费观看日本| 精品久久久久久久末码| 少妇的逼水好多| 日本午夜av视频| 亚洲国产最新在线播放| 国产成人一区二区在线| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 夫妻性生交免费视频一级片| 久久久精品94久久精品| 直男gayav资源| 一区二区三区高清视频在线| 免费大片18禁| 日本一本二区三区精品| 成人三级黄色视频| 午夜爱爱视频在线播放| 看十八女毛片水多多多| 国产亚洲91精品色在线| 亚洲国产精品久久男人天堂| 2022亚洲国产成人精品| 91精品伊人久久大香线蕉| 国产成人精品一,二区| 人人妻人人澡欧美一区二区| 久久国内精品自在自线图片| 大香蕉久久网| 男女国产视频网站| 真实男女啪啪啪动态图| 精华霜和精华液先用哪个| 日本一二三区视频观看| 久久久久久久久中文| 久久亚洲精品不卡| 色视频www国产| 欧美最新免费一区二区三区| 国产高清国产精品国产三级 | 国产一区二区亚洲精品在线观看| 成人亚洲精品av一区二区| 亚洲第一区二区三区不卡| 久久久精品欧美日韩精品| 欧美极品一区二区三区四区| 国产精品99久久久久久久久| 欧美三级亚洲精品| 在线观看av片永久免费下载| 精品人妻视频免费看| 国产国拍精品亚洲av在线观看| 校园人妻丝袜中文字幕| 成人性生交大片免费视频hd| 高清在线视频一区二区三区 | 国产不卡一卡二| 亚洲av一区综合| 国产片特级美女逼逼视频| 级片在线观看| 看十八女毛片水多多多| 美女大奶头视频| 国产视频首页在线观看| 国产美女午夜福利| 精品国内亚洲2022精品成人| 成人亚洲欧美一区二区av| 久久国产乱子免费精品| 久久综合国产亚洲精品| 又爽又黄无遮挡网站| 伦精品一区二区三区| 欧美一区二区精品小视频在线| 日韩制服骚丝袜av| 国产视频首页在线观看| 欧美xxxx性猛交bbbb| 中文资源天堂在线| 日本猛色少妇xxxxx猛交久久| 18禁动态无遮挡网站| 日韩欧美 国产精品| 日日摸夜夜添夜夜爱| 精品国内亚洲2022精品成人| 精品久久久久久电影网 | 小说图片视频综合网站| 亚洲成人精品中文字幕电影| 日本猛色少妇xxxxx猛交久久| 青青草视频在线视频观看| 美女内射精品一级片tv| 亚洲国产成人一精品久久久| 午夜日本视频在线| 午夜免费激情av| 国产午夜精品一二区理论片| 亚洲欧美精品专区久久| 精品久久久久久久久av| 国产黄色小视频在线观看| 国产亚洲av嫩草精品影院| 欧美成人精品欧美一级黄| 麻豆一二三区av精品| 亚洲av二区三区四区| 国产女主播在线喷水免费视频网站 | 人妻制服诱惑在线中文字幕| 女人久久www免费人成看片 | 亚洲av二区三区四区| 综合色av麻豆| 国产大屁股一区二区在线视频| 亚洲精品,欧美精品| 免费观看的影片在线观看| 欧美激情国产日韩精品一区| 国产高清国产精品国产三级 | 在线a可以看的网站| 亚洲av中文av极速乱| 国产男人的电影天堂91| 亚洲精品国产av成人精品| a级毛色黄片| kizo精华| 午夜福利在线在线| 亚洲最大成人av| 婷婷色av中文字幕| 天堂中文最新版在线下载 | 草草在线视频免费看| 国产麻豆成人av免费视频| 夫妻性生交免费视频一级片| av在线播放精品| 99久久精品一区二区三区| 美女国产视频在线观看| 国产精品美女特级片免费视频播放器| 中文字幕免费在线视频6| 色吧在线观看| 97在线视频观看| 成人av在线播放网站| a级一级毛片免费在线观看| eeuss影院久久| 久99久视频精品免费| 99久久九九国产精品国产免费| 中文字幕熟女人妻在线| 国产精品野战在线观看| 九色成人免费人妻av| av黄色大香蕉| 久久久成人免费电影| 国产精品伦人一区二区| 日韩av在线免费看完整版不卡| 午夜福利视频1000在线观看| 人人妻人人澡人人爽人人夜夜 | 欧美性猛交黑人性爽| 男女那种视频在线观看| 日韩大片免费观看网站 | 国产av在哪里看| 国产淫片久久久久久久久| 精品国产一区二区三区久久久樱花 | 三级国产精品片| 午夜福利网站1000一区二区三区| 免费大片18禁| 成人亚洲精品av一区二区| 高清视频免费观看一区二区 | 亚洲乱码一区二区免费版| 免费观看人在逋| 蜜桃久久精品国产亚洲av| 人妻夜夜爽99麻豆av| 成人漫画全彩无遮挡| 成人二区视频| 色综合色国产| 一边亲一边摸免费视频| 超碰97精品在线观看| 国语自产精品视频在线第100页| 国产免费一级a男人的天堂| av又黄又爽大尺度在线免费看 | 天堂网av新在线| 亚洲国产高清在线一区二区三| 久久久成人免费电影| 少妇猛男粗大的猛烈进出视频 | 两个人视频免费观看高清| 在线播放国产精品三级| 欧美日韩精品成人综合77777| 非洲黑人性xxxx精品又粗又长| 两个人的视频大全免费| 在现免费观看毛片| av在线蜜桃| 午夜老司机福利剧场| 欧美激情久久久久久爽电影| 久久6这里有精品| 欧美日本视频| 美女国产视频在线观看| 午夜激情欧美在线| 国产黄片视频在线免费观看| 成年免费大片在线观看| 国产成人午夜福利电影在线观看| 亚洲精品日韩在线中文字幕| 一级毛片久久久久久久久女| 别揉我奶头 嗯啊视频| 日本免费一区二区三区高清不卡| 亚洲av电影在线观看一区二区三区 | 日韩欧美国产在线观看| 亚洲欧美精品自产自拍| 国产国拍精品亚洲av在线观看| 99久久精品国产国产毛片| 亚洲丝袜综合中文字幕| 日本与韩国留学比较| 日韩强制内射视频| 亚洲美女搞黄在线观看| 1000部很黄的大片| 亚洲国产精品成人久久小说| 午夜福利高清视频| 在线免费观看不下载黄p国产| 久久久久久九九精品二区国产| 亚洲欧美日韩高清专用| 成人毛片a级毛片在线播放| 欧美色视频一区免费| 国产午夜福利久久久久久| 国模一区二区三区四区视频| 夜夜看夜夜爽夜夜摸| 能在线免费看毛片的网站| .国产精品久久| 亚洲五月天丁香| 色综合色国产| 一区二区三区免费毛片| 最近中文字幕2019免费版| 好男人视频免费观看在线| 色综合色国产| av卡一久久| 黄色日韩在线| 亚洲国产最新在线播放| 免费观看在线日韩| 中文乱码字字幕精品一区二区三区 | 精品人妻偷拍中文字幕| 国产成人aa在线观看| 边亲边吃奶的免费视频| 水蜜桃什么品种好| 日韩精品青青久久久久久| 热99re8久久精品国产| 亚洲激情五月婷婷啪啪| 麻豆乱淫一区二区| 汤姆久久久久久久影院中文字幕 | 亚洲不卡免费看| 亚洲欧洲国产日韩| 久久精品国产鲁丝片午夜精品| 久久婷婷人人爽人人干人人爱| av线在线观看网站| 日韩欧美三级三区| 国产淫语在线视频| 噜噜噜噜噜久久久久久91| 少妇的逼水好多| 国产中年淑女户外野战色| 成人毛片60女人毛片免费| 精品欧美国产一区二区三| 三级国产精品欧美在线观看| 成人av在线播放网站| 日韩人妻高清精品专区| 桃色一区二区三区在线观看| 只有这里有精品99| 女的被弄到高潮叫床怎么办| 亚洲怡红院男人天堂| 级片在线观看| 天堂√8在线中文| 日韩一区二区视频免费看| 国产黄a三级三级三级人| 麻豆一二三区av精品| 亚洲精品色激情综合| 极品教师在线视频| 91aial.com中文字幕在线观看| 韩国高清视频一区二区三区| 午夜福利网站1000一区二区三区| 亚洲国产最新在线播放| 亚洲中文字幕一区二区三区有码在线看| 久久久色成人| 久久欧美精品欧美久久欧美| 嘟嘟电影网在线观看| 在线观看一区二区三区| 又爽又黄a免费视频| 亚洲av.av天堂| 老师上课跳d突然被开到最大视频| 日韩欧美三级三区| 天堂av国产一区二区熟女人妻| 91久久精品国产一区二区三区| 尾随美女入室| 久久6这里有精品| 免费av不卡在线播放| 最近手机中文字幕大全| 久久热精品热| 欧美丝袜亚洲另类| 黄片wwwwww| 久久久久久久久久久丰满| 日韩大片免费观看网站 | 亚洲国产精品专区欧美| 日韩视频在线欧美| 舔av片在线| 亚洲欧美日韩东京热| 嫩草影院入口| av在线蜜桃| 99热6这里只有精品| 亚洲欧美日韩无卡精品| 在线播放国产精品三级| 国产av不卡久久| 久久久欧美国产精品| 日本免费a在线| 三级经典国产精品| 亚洲四区av| 我的老师免费观看完整版| 三级毛片av免费| 少妇人妻一区二区三区视频| 色尼玛亚洲综合影院| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美日韩高清专用| 蜜臀久久99精品久久宅男| 免费看av在线观看网站| 久久久成人免费电影| 亚洲在线观看片| 男女边吃奶边做爰视频| 一个人看视频在线观看www免费| 少妇猛男粗大的猛烈进出视频 | 国产精品1区2区在线观看.| 午夜福利视频1000在线观看| 日韩av在线免费看完整版不卡| 久久久国产成人精品二区| 国产精品av视频在线免费观看| 又黄又爽又刺激的免费视频.| 有码 亚洲区| 欧美日本视频| av卡一久久| 精品一区二区三区人妻视频| 建设人人有责人人尽责人人享有的 | 特大巨黑吊av在线直播| 日韩制服骚丝袜av| 国产色婷婷99| 一级毛片久久久久久久久女| 成人亚洲欧美一区二区av| 国产片特级美女逼逼视频| 亚洲真实伦在线观看| 精品99又大又爽又粗少妇毛片| 97超视频在线观看视频| 国产成人午夜福利电影在线观看| 亚洲国产精品成人综合色| 日本-黄色视频高清免费观看| 51国产日韩欧美| 婷婷色综合大香蕉| 99热这里只有是精品50| 成人漫画全彩无遮挡| 久久久久久久久久久免费av| 精品一区二区三区视频在线| 2021天堂中文幕一二区在线观| 九九爱精品视频在线观看| 亚洲国产欧美在线一区| 国产伦在线观看视频一区| 91精品伊人久久大香线蕉| 日本黄大片高清| 麻豆av噜噜一区二区三区| 国产伦精品一区二区三区视频9| 啦啦啦观看免费观看视频高清| 亚洲av电影在线观看一区二区三区 | 成年av动漫网址| 亚洲av福利一区| 91狼人影院| 国模一区二区三区四区视频| 亚洲精品亚洲一区二区| 欧美高清性xxxxhd video| 五月玫瑰六月丁香| 久久久久免费精品人妻一区二区| 午夜福利在线观看吧| 在线观看av片永久免费下载| 欧美另类亚洲清纯唯美| 国产av不卡久久| 搞女人的毛片| 男插女下体视频免费在线播放| 婷婷色麻豆天堂久久 | 亚洲av成人av| 国产免费男女视频| 午夜免费激情av| 国产极品天堂在线| 精品久久久久久久久久久久久| 国国产精品蜜臀av免费| 女人久久www免费人成看片 | 欧美bdsm另类| 久久精品熟女亚洲av麻豆精品 | 99久久精品国产国产毛片| 久久99精品国语久久久| 在线播放国产精品三级| 亚洲最大成人av| 免费av不卡在线播放| 最近2019中文字幕mv第一页| 亚洲美女视频黄频| 秋霞伦理黄片| 国产精品久久久久久av不卡| 日韩一本色道免费dvd| 蜜桃久久精品国产亚洲av| 精品久久久久久久久av| 伦理电影大哥的女人| 成人性生交大片免费视频hd| 久久久久精品久久久久真实原创| 国产欧美日韩精品一区二区| 亚洲国产精品成人综合色| 中文天堂在线官网| 内地一区二区视频在线| 亚洲欧洲日产国产| 国语自产精品视频在线第100页| 国产精品美女特级片免费视频播放器| 精品国产三级普通话版| 国产在线男女| 欧美成人一区二区免费高清观看| 午夜日本视频在线| 天堂影院成人在线观看| 偷拍熟女少妇极品色| av.在线天堂| 国产亚洲一区二区精品| 日韩三级伦理在线观看| 日本熟妇午夜| 国产真实乱freesex| 久久精品久久久久久噜噜老黄 | 日日摸夜夜添夜夜爱| 国产爱豆传媒在线观看| 亚洲不卡免费看| 九九热线精品视视频播放| 一级黄色大片毛片| 成年免费大片在线观看| 五月玫瑰六月丁香| 九草在线视频观看| 最近2019中文字幕mv第一页| 国产成人91sexporn| 中国美白少妇内射xxxbb| 男插女下体视频免费在线播放| 国产一区二区三区av在线| 国产色婷婷99| 亚洲aⅴ乱码一区二区在线播放| 久久久久久九九精品二区国产| 高清毛片免费看| 亚洲最大成人av| 亚洲中文字幕日韩| 国产v大片淫在线免费观看| 高清av免费在线| 久久久久久久久久久免费av| 精品无人区乱码1区二区| 免费电影在线观看免费观看| 亚洲,欧美,日韩| 国内精品宾馆在线| 26uuu在线亚洲综合色| 少妇的逼水好多| 亚洲综合精品二区| 亚洲自拍偷在线| 大又大粗又爽又黄少妇毛片口| 免费搜索国产男女视频| 亚洲精品日韩在线中文字幕| 国产 一区 欧美 日韩| 国产高清国产精品国产三级 | 国产69精品久久久久777片| 国产精品一区二区三区四区久久| 国产精品一区二区三区四区免费观看| 久久精品夜色国产| 亚洲一级一片aⅴ在线观看| 日本熟妇午夜| 久久精品夜色国产| 在线免费观看的www视频| 伊人久久精品亚洲午夜| 国产精品美女特级片免费视频播放器| 美女国产视频在线观看| 日韩国内少妇激情av| 国产精品久久久久久精品电影小说 | av国产免费在线观看| 国产激情偷乱视频一区二区| 一本久久精品| 3wmmmm亚洲av在线观看| 人人妻人人看人人澡| 2021少妇久久久久久久久久久| 色噜噜av男人的天堂激情| 久久久欧美国产精品| 毛片女人毛片| 亚洲av成人精品一二三区| 国产午夜精品论理片| 在线免费十八禁| 精品一区二区免费观看| 国产探花极品一区二区| 国产成人aa在线观看| 亚洲一级一片aⅴ在线观看| 欧美激情久久久久久爽电影| 久久精品国产亚洲网站| 中文精品一卡2卡3卡4更新| 久久久久久久久久久免费av| 日本黄色片子视频| 精品久久国产蜜桃| 夜夜看夜夜爽夜夜摸| 中文字幕免费在线视频6|