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

    基于GIS 的地震時空統(tǒng)計(jì)分析
    ——以長三角地區(qū)為例

    2022-06-16 07:48:26余其鵬董建平
    關(guān)鍵詞:震級長三角時空

    余其鵬,董建平

    (鎮(zhèn)江市地震臺,江蘇 鎮(zhèn)江 212000)

    0 引言

    長江三角洲地區(qū)(以下簡稱長三角地區(qū)),包括上海市、江蘇省、浙江省、安徽省,包含地級市和直轄市共41 個,是中國經(jīng)濟(jì)發(fā)展最活躍、開放程度最高、創(chuàng)新能力最強(qiáng)的區(qū)域之一,在國家現(xiàn)代化建設(shè)大局和全方位開放格局中具有舉足輕重的戰(zhàn)略地位。同時,歷史上長三角地區(qū)也是地震多發(fā)的地區(qū),多年來一直被列為全國地震重點(diǎn)監(jiān)視防御區(qū),僅次于首都圈。該地區(qū)歷史上發(fā)生過多次中強(qiáng)地震活動,其中5.0級以上地震38 次,6.0 級至7.0 級地震6 次,目前該地區(qū)沒有7.0 級以上地震的記錄[1]。因此針對該地區(qū)地震數(shù)據(jù)進(jìn)行時空分析,研究地震活動規(guī)律,對科學(xué)指導(dǎo)防震減災(zāi)工作以及城市群規(guī)劃和重大工程建設(shè)具有一定的指導(dǎo)意義和現(xiàn)實(shí)意義。

    關(guān)于對地震數(shù)據(jù)進(jìn)行時空統(tǒng)計(jì)分析的相關(guān)研究,國內(nèi)外許多學(xué)者都對其進(jìn)行過探索。如有學(xué)者在我國西南地區(qū),通過對1500—1992 年4.7 級以上的地震進(jìn)行統(tǒng)計(jì)分析,探索了地震活動的時間和空間變化情況[2]。對地震做時空統(tǒng)計(jì)分析,從而發(fā)現(xiàn)地震熱點(diǎn)區(qū)域,也是一種普遍的研究方法,有學(xué)者通過對2005—2007 年四川地區(qū)地震數(shù)據(jù)做空間統(tǒng)計(jì)分析,從而發(fā)現(xiàn)明顯的地震熱點(diǎn)區(qū)域[3]。地震空間統(tǒng)計(jì)方法很多,有學(xué)者就利用核密度估計(jì)法來在大時間尺度上研究涇洛河流域的地震時空格局,最終發(fā)現(xiàn)該流域空間分布呈W 型,高密度地震區(qū)域與斷裂區(qū)域有部分重疊的情況[4]。國外也有許多學(xué)者通過時空統(tǒng)計(jì)方法研究地震事件,著名的Gutenberg-Richter 定律,就是采用統(tǒng)計(jì)學(xué)方法對全球地震數(shù)據(jù)進(jìn)行研究,得出震級和頻率的關(guān)系[5]。有采用非廣延統(tǒng)計(jì)方法研究連續(xù)地震的時空分布,研究時間分布表現(xiàn)和冪律統(tǒng)計(jì)關(guān)系[6],也有采用GIS 中聚類分析方法研究地震的空間分布,通過Gutenberg-Richter 定律來分析地震發(fā)生的概率,探討地震在空間分布上的聚集特征[7]。

    這些研究都在地震時空統(tǒng)計(jì)分析方面作了深入的探索,為地震研究做出了重要的貢獻(xiàn)。研究中針對地震數(shù)據(jù)使用的統(tǒng)計(jì)方法大多是空間統(tǒng)計(jì)中的某一種,對地震數(shù)據(jù)開展較為全面的分析相對較少。本文通過GIS 空間分析方法,對長三角地區(qū)近14 年的ML≥0.1 的地震數(shù)據(jù)進(jìn)行時空統(tǒng)計(jì)分析,通過利用密度分析方法和空間統(tǒng)計(jì)分析方法,對研究區(qū)地震點(diǎn)的空間分布和地震事件震級、頻率、位置、時間等進(jìn)行全面研究,對區(qū)域內(nèi)地震的分布規(guī)律研究具有一定的借鑒意義。

    1 研究區(qū)概況與數(shù)據(jù)來源

    1.1 研究區(qū)概況

    長三角地區(qū)位于中國長江的下游地區(qū),瀕臨黃海與東海,地處江海交匯之地,沿江沿海港口眾多,是長江入海之前形成的沖積平原。區(qū)域面積35.8 萬平方千米,境內(nèi)有多條地震帶,如郯廬斷裂帶、揚(yáng)州—南京—江南斷裂帶、上?!贾輸嗔褞?、霍山斷裂帶等[8],是地震多發(fā)區(qū)域,2008 年12 月30 至2021 年7 月20 日近14 年來,長三角地區(qū)大約發(fā)生了ML≥0.1 的地震4835 次,如圖1 所示。

    圖1 長三角地區(qū)近14 年地震分布圖Fig.1 Seismic distribution map of Yangtze River Delta region in recent 14 years

    1.2 數(shù)據(jù)來源與處理

    關(guān)于地震研究數(shù)據(jù)的獲取,互聯(lián)網(wǎng)上有著豐富的資源,比如美國地質(zhì)調(diào)查局(USGS)、中國地震臺網(wǎng)中心等。本文所用的研究區(qū)地震點(diǎn)數(shù)據(jù)和斷層數(shù)據(jù)來源于中國地震臺網(wǎng)中心國家地震科學(xué)數(shù)據(jù)中心(http://data.earthquake.cn)。地震點(diǎn)數(shù)據(jù),詳細(xì)記錄了發(fā)震時間、經(jīng)緯度、震源深度、震級、參考位置等信息。斷層數(shù)據(jù)來源于中國大陸1/400 萬中國活動構(gòu)造空間數(shù)據(jù)庫。本文獲取的地震點(diǎn)數(shù)據(jù),存儲格式為.xls,由于數(shù)據(jù)多,為了方便管理,研究中將其點(diǎn)位信息(坐標(biāo)、震級、發(fā)震時間、地點(diǎn)等)以及斷層數(shù)據(jù)導(dǎo)入新建的Geodatabase 空間數(shù)據(jù)庫[9]。

    2 研究框架和分析方法

    本文依托GIS 強(qiáng)大的空間分析方法,借助ArcGIS 軟件,對長三角地區(qū)近14 年ML≥0.1 的地震數(shù)據(jù)進(jìn)行時空統(tǒng)計(jì)分析,通過核密度分析、震級頻率分析、區(qū)域-頻率分析、空間統(tǒng)計(jì)分析方法等多種方法進(jìn)行全面研究。

    2.1 地震事件點(diǎn)密度分析

    采用核密度估計(jì)法,在同一面積尺度上對所有地震點(diǎn)進(jìn)行密度分析,將分析結(jié)果用二維灰度進(jìn)行可視化表達(dá),這樣可以直觀簡明展現(xiàn)研究區(qū)域內(nèi)地震事件的聚集、離散情況。

    本文對長三角地區(qū)2008 年12 月30 至2021 年7 月20 日近14 年來ML≥0.1 的地震事件進(jìn)行統(tǒng)計(jì)分析,在ArcGIS 軟件中進(jìn)行具體操作,使用Spatial Analyst 工具的核密度分析,搜索半徑為默認(rèn),Population 為震級字段,賦予震級比其他要求更大的權(quán)重,這樣更加合理。將每一要素按照震級大小計(jì)算,最后對結(jié)果進(jìn)行重分類,將生成的密度圖分為6 個等級,最終生成近14 年來長三角地區(qū)地震熱力圖,如圖2 所示。由圖可以看出:長三角地區(qū)地震空間分布的核中心主要集中在六安、溫州、鹽城、揚(yáng)州、寧波、合肥、淮安幾個地方。

    圖2 ML≥0.1 的地震事件的空間分布Fig.2 Spatial distribution of seismic events

    另外,將長三角地區(qū)斷裂帶分布和核密度分布進(jìn)行疊加分析,如圖3 所示,霍山斷裂、余姚—麗水?dāng)嗔?、奉化—麗水?dāng)嗔?、?zhèn)海—溫州斷裂、洪澤—流均溝斷裂、涇口—流均斷裂、鹽城斷裂、穎上—定遠(yuǎn)斷裂、昌化—普陀斷裂為地震主要發(fā)生區(qū)域,重點(diǎn)區(qū)域均有地震斷裂帶穿過。

    圖3 斷層和ML≥0.1 的地震事件空間分布圖Fig.3 Spatial distribution of faults and seismic events

    2.2 地震震級與頻率分析

    大量研究者通過對地震頻率和震級(G-R)關(guān)系式的研究,發(fā)現(xiàn)地震震級和地震頻率的對數(shù)有線性關(guān)系,如圖4 所示。此關(guān)系式在地震活動性和危險性預(yù)測方面有比較好的效果[10],而且一般關(guān)系模型中震級和頻率的關(guān)系往往是冪律分布。根據(jù)Guterg 等[5]研究發(fā)現(xiàn)得出的一個經(jīng)驗(yàn)公式,如公式(1)所示:

    圖4 一般G-R 關(guān)系分布圖Fig.4 General G-R relationship distribution diagram

    其中,N表示地震頻率,M表示地震震級,a表示區(qū)域內(nèi)一定事件的地震活動強(qiáng)度,b表示地震的危險程度,它通常是地震研究的一個重要參數(shù),a、b兩個值都可以通過最小二乘法計(jì)算得到[11]。

    對長三角地區(qū)2008 年12 月30 日至2021年7 月20 日近14 年所有ML≥0.1 的震級和地震地點(diǎn)進(jìn)行統(tǒng)計(jì),得到震級和頻率對應(yīng)的表格,部分如表1 所示,由統(tǒng)計(jì)結(jié)果可以看出地震頻率隨震級增大而逐漸減小。

    表1 長三角地區(qū)部分震級-頻率表

    利用表格繪制長三角地區(qū)地震震級和地震頻率之間的關(guān)系圖,如圖5 所示,可以直觀的看到隨著震級增大,地震發(fā)生的頻率在逐漸減小,2008 年12 月30 至2021 年7 月20 以來發(fā)生ML≥4.0 地震只有7 次,而ML<4 的地震就有4828 次。利用最小二乘法計(jì)算得到a、b的值,進(jìn)一步得到震級-頻率的函數(shù)關(guān)系,得到其公式大致如公式(2) 所示,結(jié)果表明,研究區(qū)內(nèi)地震頻率和震級的關(guān)系符合Gutenberg-Richter定律,即G-R 關(guān)系式。

    圖5 研究區(qū)ML≥0.1 的地震震級-頻率關(guān)系圖Fig.5 Magnitude-frequency diagram of earthquakes in the study area

    雖然利用這種震級-頻率分布函數(shù)來準(zhǔn)確預(yù)測地震并不科學(xué),但這種研究方法對于研究地震數(shù)據(jù)的動力學(xué)和相關(guān)長期預(yù)測是有一定幫助的。

    2.3 地震區(qū)域與頻率分析

    研究區(qū)區(qū)域面積35.8 萬平方千米,包含40個地級市和1 個直轄市,共307 個縣區(qū),研究以41 個行政區(qū)為對象,通過臺網(wǎng)中心提供的震中參考位置,對2008 年12 月30 日至2021 年7 月20 日期間各個地區(qū)發(fā)生的ML≥0.1 地震的頻率進(jìn)行統(tǒng)計(jì),其部分結(jié)果如圖6 所示。

    圖6 長三角部分地區(qū)ML≥0.1 的地震發(fā)生的數(shù)量圖Fig.6 Numbes of earthquakes in parts of the Yangtze River Delta

    同時結(jié)合41 個行政區(qū)域內(nèi)發(fā)生的ML≥0.1地震的頻率,按照由高到低的頻率進(jìn)行分級,劃分為4 個等級,頻率為0 的為一級,共2 個地區(qū),占該區(qū)域的4.9%;頻率在1~100 的為一級,共32 個地區(qū),占總區(qū)域的78%;頻率在101~500 的為一級,共5 個地區(qū),占區(qū)域的12.2%;頻率在501~1640 的為一級,共2 個地區(qū),占區(qū)域的4.9%。長三角地區(qū)地震頻率的空間分布,如圖7 所示。

    圖7 長三角地區(qū)ML≥0.1 的地震頻率空間分布圖Fig.7 Spatial distribution of seismic frequency in Yangtze River Delta region

    由統(tǒng)計(jì)結(jié)果和空間分布圖可以看出,在統(tǒng)計(jì)時間內(nèi),研究區(qū)共7 個區(qū)域(占研究區(qū)17.1%)發(fā)生了100 次以上ML≥0.1 的地震,這7 個區(qū)域地震總體頻率達(dá)到3662 次。研究表明,這7個區(qū)域是地震發(fā)生的高頻地區(qū)。

    2.4 地震發(fā)生時間和頻率分析

    基于以上研究,本文將長三角區(qū)域在不同時間內(nèi)發(fā)生的ML≥0.1 的地震事件進(jìn)行統(tǒng)計(jì),通過統(tǒng)計(jì)長三角地區(qū)2009—2021 年每一年發(fā)生地震的次數(shù)(統(tǒng)計(jì)中X 軸表示時間,以年為單位,Y 軸表示一年當(dāng)中發(fā)生地震的次數(shù)),以研究長三角地區(qū)每年地震序列的變化情況,最終的統(tǒng)計(jì)結(jié)果如圖8 所示。由圖中可以看出:長三角地區(qū)在2009—2021 年每年發(fā)生地震的頻率相對比較一致,略有起伏。據(jù)中國地震局統(tǒng)計(jì)顯示,2014 年大陸地區(qū)發(fā)生地震災(zāi)害事件10次,其中特別重大地震災(zāi)害事件1 次,較大地震災(zāi)害事件4 次,一般地震災(zāi)害5 次。而其中10 月25 日在浙江省文成縣發(fā)生了4.2 級地震,在這個時間點(diǎn)前后,該地區(qū)發(fā)生了多次小震和余震,故2014 年長三角地區(qū)的地震頻率高達(dá)1592 次,遠(yuǎn)多于其他年份。由此可見,在某些較大地震前后,發(fā)生地震的頻率也會增高。

    圖8 長三角地區(qū)2009—2021 年每年ML≥0.1 的地震發(fā)生次數(shù)Fig.8 Annual number of earthquakes in the Yangtze River Delta region from 2009 to 2021

    另外,針對2014 年地震頻率較高的情況,我們對2014 年1 月1 日至12 月31 日長三角地區(qū)每日ML≥0.1 的地震發(fā)生頻率進(jìn)行統(tǒng)計(jì),其中Y 軸表示每天發(fā)生地震的頻率,X 軸表示時間,統(tǒng)計(jì)結(jié)果如圖9 所示,由統(tǒng)計(jì)圖可以看出:長三角地區(qū)2014 年8 月份之前發(fā)生地震的頻率都比較低,每天幾乎都是個位數(shù),但8 月以后地震次數(shù)明顯增加,最高一天多達(dá)65 次地震。在10 月25 日浙江省文成縣發(fā)生的4.2 級地震前后,地震頻率明顯變高。從地質(zhì)構(gòu)造角度來解釋,一般在強(qiáng)震發(fā)生后,由于受到地震斷裂帶影響,造成板塊之間相互碰撞、相互擠壓,這就造成了地殼中的應(yīng)力在一定程度上會產(chǎn)生集聚,隨著應(yīng)力不斷積累,達(dá)到一定臨界點(diǎn)就會釋放,原來地質(zhì)構(gòu)造中的系統(tǒng)平衡就會被打破。所以,板塊構(gòu)造中的應(yīng)力場也會隨之進(jìn)行新的調(diào)整和分配,因此導(dǎo)致震中區(qū)域的應(yīng)力場會隨著新的調(diào)整產(chǎn)生新的變化,這樣就會造成一系列小的地震。因此,從分析中我們可以得出,在2014 年該地區(qū)在較大地震發(fā)生之前,每天地震發(fā)生的頻率不高,相對比較穩(wěn)定,發(fā)生較大地震后,地震頻率會明顯發(fā)生變化,由統(tǒng)計(jì)結(jié)果可以看出每天能發(fā)生幾十次甚至上百次地震。

    圖9 長三角地區(qū)2014 年每天發(fā)生ML≥0.1 的地震的頻率圖Fig.9 Frequency chart of daily earthquakes in the Yangtze River Delta region in 2014

    3 結(jié)論與討論

    本文通過GIS 空間統(tǒng)計(jì)和分析方法,詳細(xì)介紹了利用ArcGIS 軟件對長三角地區(qū)近14 年ML≥0.1 的地震數(shù)據(jù)進(jìn)行時空統(tǒng)計(jì)分析方法,闡述了通過利用核密度分析方法和空間統(tǒng)計(jì)分析方法,對研究區(qū)地震點(diǎn)的空間分布和地震事件震級、頻率、位置等進(jìn)行全面研究。研究表明,長三角地區(qū)地震頻率的對數(shù)和地震震級之間存在線性關(guān)系,符合Gutenberg-Richter 定律。大部分地震事件主要發(fā)生在比較集中的幾個區(qū)域,且較大地震發(fā)生前后,存在區(qū)域內(nèi)的地震頻率有較明顯變化。研究對探尋一定時期內(nèi)長三角地區(qū)地震活動規(guī)律,以及對該地區(qū)城市群規(guī)劃及重大工程建設(shè)有一定指導(dǎo)意義,對科學(xué)指導(dǎo)防災(zāi)減災(zāi)工作具有一定現(xiàn)實(shí)意義。

    本文的研究還有不足和需要改進(jìn)的地方,比如受限于臺網(wǎng)中心提供的數(shù)據(jù),用到的地震數(shù)據(jù)事件時間的跨度還不大,地震斷裂數(shù)據(jù)庫比例尺較小、驅(qū)動機(jī)制分析不足,分析還存?zhèn)€例性和片面性等問題,如何根據(jù)各地區(qū)更加詳細(xì)的地震資料,結(jié)合地震造成的房屋破壞、人員傷亡、經(jīng)濟(jì)損失,依托規(guī)劃、交通、人口等數(shù)據(jù)進(jìn)行更加深入分析將是進(jìn)一步思考的方向。

    猜你喜歡
    震級長三角時空
    基于累積絕對位移值的震級估算方法
    跨越時空的相遇
    緊扣一體化 全面融入長三角
    地震后各國發(fā)布的震級可能不一樣?
    鏡中的時空穿梭
    新震級國家標(biāo)準(zhǔn)在大同臺的應(yīng)用與評估
    山西地震(2020年1期)2020-04-08 07:34:26
    “首屆長三角新青年改稿會”作品選
    詩歌月刊(2019年7期)2019-08-29 01:46:44
    玩一次時空大“穿越”
    時空之門
    長三角瞭望
    丝袜喷水一区| 国产精品一二三区在线看| 亚洲精品国产av成人精品 | 日韩中字成人| 国产av在哪里看| 国产精品久久久久久亚洲av鲁大| 男女做爰动态图高潮gif福利片| 99热网站在线观看| 天堂动漫精品| 亚洲国产欧美人成| 男女视频在线观看网站免费| 精品人妻熟女av久视频| 国产精品一区www在线观看| 麻豆久久精品国产亚洲av| 国内少妇人妻偷人精品xxx网站| 国产精品福利在线免费观看| 国产一区亚洲一区在线观看| 女人十人毛片免费观看3o分钟| 九色成人免费人妻av| 露出奶头的视频| 亚洲国产精品成人综合色| 欧美在线一区亚洲| 性插视频无遮挡在线免费观看| 校园人妻丝袜中文字幕| 老司机影院成人| 嫩草影院入口| 99热6这里只有精品| 国产黄a三级三级三级人| 一级毛片久久久久久久久女| 欧美一区二区亚洲| 我的女老师完整版在线观看| 久久精品人妻少妇| 你懂的网址亚洲精品在线观看 | 一夜夜www| av视频在线观看入口| 午夜精品在线福利| 亚洲精品456在线播放app| 中文字幕av在线有码专区| 国产精品国产高清国产av| 中文亚洲av片在线观看爽| 悠悠久久av| 日日啪夜夜撸| 在线天堂最新版资源| 成人毛片a级毛片在线播放| 色在线成人网| 1000部很黄的大片| 精品国产三级普通话版| 成人三级黄色视频| 国国产精品蜜臀av免费| 亚洲va在线va天堂va国产| 成人午夜高清在线视频| 女生性感内裤真人,穿戴方法视频| 日韩一本色道免费dvd| 久久久久国产精品人妻aⅴ院| av在线播放精品| 伦精品一区二区三区| 乱系列少妇在线播放| 久久人人精品亚洲av| 日韩大尺度精品在线看网址| 亚洲精品456在线播放app| 亚洲欧美中文字幕日韩二区| 少妇人妻精品综合一区二区 | 少妇高潮的动态图| 久久精品国产亚洲av涩爱 | 极品教师在线视频| 午夜福利在线观看吧| 日本三级黄在线观看| 一个人观看的视频www高清免费观看| 国产私拍福利视频在线观看| 免费av毛片视频| 亚洲成人中文字幕在线播放| 午夜福利18| 国产一区二区三区av在线 | 草草在线视频免费看| 国产黄片美女视频| 国产大屁股一区二区在线视频| 搡老妇女老女人老熟妇| 亚洲欧美日韩卡通动漫| 国产av不卡久久| 99热只有精品国产| 日日撸夜夜添| 香蕉av资源在线| 成人三级黄色视频| 午夜福利18| 成人毛片a级毛片在线播放| av.在线天堂| 高清午夜精品一区二区三区 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲av熟女| 少妇的逼好多水| 国产精品亚洲一级av第二区| 亚洲欧美精品自产自拍| 两个人的视频大全免费| 久久精品91蜜桃| 小说图片视频综合网站| 成人永久免费在线观看视频| 成人高潮视频无遮挡免费网站| 一个人看的www免费观看视频| 亚洲色图av天堂| 大又大粗又爽又黄少妇毛片口| 亚洲av美国av| 亚洲中文字幕日韩| 麻豆乱淫一区二区| 97在线视频观看| 亚洲国产欧洲综合997久久,| 五月玫瑰六月丁香| 亚洲成人中文字幕在线播放| 欧美日韩在线观看h| 日日摸夜夜添夜夜添av毛片| 亚洲av第一区精品v没综合| 黄色一级大片看看| 欧美潮喷喷水| 天堂影院成人在线观看| 插阴视频在线观看视频| 一a级毛片在线观看| 国产精品久久久久久久电影| 露出奶头的视频| 日本 av在线| 久久亚洲国产成人精品v| 成人二区视频| 岛国在线免费视频观看| 国产又黄又爽又无遮挡在线| 国产精品美女特级片免费视频播放器| 一区二区三区免费毛片| 亚洲最大成人av| 日本 av在线| 久久久精品94久久精品| 欧美一区二区亚洲| av在线天堂中文字幕| 床上黄色一级片| 狂野欧美激情性xxxx在线观看| 少妇猛男粗大的猛烈进出视频 | 99久久中文字幕三级久久日本| 免费人成视频x8x8入口观看| 亚洲最大成人av| 国产欧美日韩精品亚洲av| 三级国产精品欧美在线观看| 国产高潮美女av| 国语自产精品视频在线第100页| 欧美人与善性xxx| 国产精品一区二区三区四区免费观看 | 国产一区亚洲一区在线观看| 99热网站在线观看| 日韩三级伦理在线观看| 久久久久久久久大av| 日本免费一区二区三区高清不卡| 麻豆久久精品国产亚洲av| 欧美一区二区精品小视频在线| 亚洲精品日韩在线中文字幕 | 国产av一区在线观看免费| a级毛色黄片| 日本色播在线视频| 久久久午夜欧美精品| 精品一区二区三区人妻视频| 精品午夜福利视频在线观看一区| а√天堂www在线а√下载| 麻豆国产av国片精品| 久久久久性生活片| 久久久久性生活片| 国产精品永久免费网站| 午夜日韩欧美国产| 网址你懂的国产日韩在线| 两个人的视频大全免费| 波多野结衣高清无吗| 2021天堂中文幕一二区在线观| 97碰自拍视频| 亚洲国产精品sss在线观看| 久久婷婷人人爽人人干人人爱| 变态另类丝袜制服| 黄色视频,在线免费观看| av视频在线观看入口| 十八禁网站免费在线| 老熟妇乱子伦视频在线观看| 在线观看免费视频日本深夜| 99国产极品粉嫩在线观看| 男女做爰动态图高潮gif福利片| 男人舔女人下体高潮全视频| 久久人人精品亚洲av| 亚洲精品一区av在线观看| 搡老岳熟女国产| 国产精品,欧美在线| 在线观看一区二区三区| 国国产精品蜜臀av免费| 97热精品久久久久久| 欧美丝袜亚洲另类| 男人狂女人下面高潮的视频| 成人av在线播放网站| 午夜精品国产一区二区电影 | 国产成人a∨麻豆精品| 亚洲va在线va天堂va国产| 色播亚洲综合网| 18禁裸乳无遮挡免费网站照片| 精品欧美国产一区二区三| 国产真实乱freesex| 美女高潮的动态| 成人鲁丝片一二三区免费| 亚洲人与动物交配视频| 午夜爱爱视频在线播放| 精品日产1卡2卡| 亚洲天堂国产精品一区在线| 日日摸夜夜添夜夜添小说| 午夜激情欧美在线| 国产精品久久久久久亚洲av鲁大| 一本久久中文字幕| 啦啦啦韩国在线观看视频| 美女黄网站色视频| 亚洲熟妇中文字幕五十中出| 色噜噜av男人的天堂激情| 亚洲激情五月婷婷啪啪| 91在线观看av| 欧美xxxx性猛交bbbb| 精品人妻熟女av久视频| 一进一出抽搐动态| 一级av片app| 成年女人毛片免费观看观看9| 99久久成人亚洲精品观看| 少妇猛男粗大的猛烈进出视频 | 国产成人精品久久久久久| 亚洲国产欧洲综合997久久,| 国产一区亚洲一区在线观看| 欧美日本亚洲视频在线播放| 99在线视频只有这里精品首页| 简卡轻食公司| 啦啦啦啦在线视频资源| 非洲黑人性xxxx精品又粗又长| 亚洲国产精品国产精品| 真人做人爱边吃奶动态| 久久久久国产精品人妻aⅴ院| 国产免费一级a男人的天堂| 日日啪夜夜撸| 成年版毛片免费区| 国产精品久久电影中文字幕| 男人狂女人下面高潮的视频| 大又大粗又爽又黄少妇毛片口| 大香蕉久久网| 久久久久国产网址| 色在线成人网| 国产一区二区亚洲精品在线观看| 丝袜美腿在线中文| 美女 人体艺术 gogo| 直男gayav资源| 插阴视频在线观看视频| 春色校园在线视频观看| 亚洲高清免费不卡视频| 亚洲美女视频黄频| eeuss影院久久| 久久国产乱子免费精品| 少妇人妻一区二区三区视频| 天堂√8在线中文| 97热精品久久久久久| 亚洲内射少妇av| 天堂动漫精品| 夜夜爽天天搞| 国产成人aa在线观看| 嫩草影院入口| 91久久精品电影网| 伊人久久精品亚洲午夜| 伦理电影大哥的女人| 亚洲人成网站在线播放欧美日韩| 亚洲国产精品sss在线观看| 免费av毛片视频| 成人漫画全彩无遮挡| 精品欧美国产一区二区三| 国产91av在线免费观看| 中文字幕av在线有码专区| 日本黄色视频三级网站网址| 男人舔奶头视频| 高清日韩中文字幕在线| 91在线精品国自产拍蜜月| 一进一出抽搐动态| av在线老鸭窝| 免费av不卡在线播放| 在线观看一区二区三区| 噜噜噜噜噜久久久久久91| 日本爱情动作片www.在线观看 | 午夜精品国产一区二区电影 | 青春草视频在线免费观看| 免费看光身美女| 伊人久久精品亚洲午夜| 可以在线观看毛片的网站| 国产高清不卡午夜福利| 国产一区二区三区av在线 | 日韩欧美 国产精品| 又爽又黄无遮挡网站| 99在线视频只有这里精品首页| av在线播放精品| 99久久精品国产国产毛片| 秋霞在线观看毛片| 在线播放国产精品三级| 久久久久免费精品人妻一区二区| 在线a可以看的网站| 精华霜和精华液先用哪个| a级毛片免费高清观看在线播放| 女生性感内裤真人,穿戴方法视频| 97在线视频观看| 欧美一区二区亚洲| 一边摸一边抽搐一进一小说| 国产亚洲精品久久久com| 男插女下体视频免费在线播放| 成人亚洲精品av一区二区| 亚洲在线自拍视频| 亚洲国产精品国产精品| 日韩一区二区视频免费看| 真实男女啪啪啪动态图| 欧美潮喷喷水| 99视频精品全部免费 在线| 午夜影院日韩av| 国产女主播在线喷水免费视频网站 | 国产综合懂色| 国产精品一区二区性色av| 色综合色国产| 亚洲av二区三区四区| avwww免费| 毛片一级片免费看久久久久| 校园春色视频在线观看| 在线观看免费视频日本深夜| 国产欧美日韩精品亚洲av| 99热6这里只有精品| 欧美最黄视频在线播放免费| 国产乱人视频| 国产精品三级大全| 欧美丝袜亚洲另类| 久久久久国产网址| a级毛片a级免费在线| 久久久久精品国产欧美久久久| 成人av一区二区三区在线看| 嫩草影院精品99| 最近2019中文字幕mv第一页| 性插视频无遮挡在线免费观看| 国产精品一区二区三区四区久久| 日韩国内少妇激情av| 亚洲经典国产精华液单| 午夜福利在线在线| 久久久久久伊人网av| 国产欧美日韩精品亚洲av| 成人综合一区亚洲| 国产成人福利小说| 国产一区二区三区在线臀色熟女| 蜜臀久久99精品久久宅男| 国产亚洲精品综合一区在线观看| 精品人妻熟女av久视频| 女生性感内裤真人,穿戴方法视频| av免费在线看不卡| 91av网一区二区| 高清午夜精品一区二区三区 | 国产精品乱码一区二三区的特点| 菩萨蛮人人尽说江南好唐韦庄 | 女人被狂操c到高潮| 国产色爽女视频免费观看| 亚洲人成网站高清观看| 大香蕉久久网| 人妻夜夜爽99麻豆av| 精品人妻偷拍中文字幕| 欧美bdsm另类| 亚洲av成人av| 超碰av人人做人人爽久久| 校园人妻丝袜中文字幕| 我要搜黄色片| 麻豆成人午夜福利视频| 成人综合一区亚洲| 黄色一级大片看看| 久久亚洲国产成人精品v| 91午夜精品亚洲一区二区三区| 日韩欧美三级三区| 男插女下体视频免费在线播放| 少妇人妻精品综合一区二区 | 精品国产三级普通话版| 色综合亚洲欧美另类图片| 国产精品久久久久久精品电影| 久久99热这里只有精品18| 黑人高潮一二区| 日本免费a在线| 久久精品夜色国产| 亚洲av成人精品一区久久| 国产又黄又爽又无遮挡在线| 直男gayav资源| 欧美潮喷喷水| 国产探花极品一区二区| 久久久午夜欧美精品| 国产蜜桃级精品一区二区三区| АⅤ资源中文在线天堂| 九九热线精品视视频播放| 男人和女人高潮做爰伦理| 老司机福利观看| 美女 人体艺术 gogo| 欧美成人精品欧美一级黄| 国产精品一区二区免费欧美| 午夜日韩欧美国产| 嫩草影视91久久| 亚洲18禁久久av| 国产精品日韩av在线免费观看| 成人特级av手机在线观看| 床上黄色一级片| 给我免费播放毛片高清在线观看| 成人无遮挡网站| 亚洲成人久久爱视频| 国产又黄又爽又无遮挡在线| 变态另类成人亚洲欧美熟女| 亚洲综合色惰| 给我免费播放毛片高清在线观看| 国产免费一级a男人的天堂| 男女那种视频在线观看| 欧美日韩国产亚洲二区| av天堂在线播放| 69av精品久久久久久| 美女 人体艺术 gogo| 欧美潮喷喷水| 国产精品一区二区三区四区免费观看 | 可以在线观看毛片的网站| 久久精品国产亚洲av涩爱 | 国产白丝娇喘喷水9色精品| 天天躁夜夜躁狠狠久久av| 成人特级av手机在线观看| av视频在线观看入口| 老司机福利观看| 寂寞人妻少妇视频99o| 国产成年人精品一区二区| 美女 人体艺术 gogo| 久久久欧美国产精品| 午夜福利在线观看免费完整高清在 | 国产91av在线免费观看| 欧美激情国产日韩精品一区| 国产精品久久久久久亚洲av鲁大| 乱人视频在线观看| 极品教师在线视频| 欧美不卡视频在线免费观看| 少妇被粗大猛烈的视频| 少妇裸体淫交视频免费看高清| 床上黄色一级片| 国产真实乱freesex| 亚洲在线观看片| 欧美最新免费一区二区三区| 日本 av在线| 国产一区二区在线av高清观看| 欧美性猛交黑人性爽| 禁无遮挡网站| 日本撒尿小便嘘嘘汇集6| 日韩在线高清观看一区二区三区| 成人特级黄色片久久久久久久| 五月玫瑰六月丁香| 国产淫片久久久久久久久| 麻豆一二三区av精品| 最近视频中文字幕2019在线8| 欧美绝顶高潮抽搐喷水| 丝袜美腿在线中文| 插逼视频在线观看| 日韩高清综合在线| 日韩欧美 国产精品| 久久精品国产亚洲网站| 精华霜和精华液先用哪个| 欧美国产日韩亚洲一区| 国产精品人妻久久久影院| 国产激情偷乱视频一区二区| 噜噜噜噜噜久久久久久91| 欧美+日韩+精品| 毛片一级片免费看久久久久| av.在线天堂| 亚洲成人久久爱视频| 久久人人精品亚洲av| 女的被弄到高潮叫床怎么办| 国产精品久久久久久精品电影| 日韩欧美国产在线观看| 99在线人妻在线中文字幕| 欧美日韩精品成人综合77777| 99在线视频只有这里精品首页| 一级黄色大片毛片| 少妇猛男粗大的猛烈进出视频 | 一级av片app| 又粗又爽又猛毛片免费看| 99热这里只有是精品在线观看| 在线观看av片永久免费下载| 国产淫片久久久久久久久| 亚洲中文字幕一区二区三区有码在线看| 亚洲成人久久爱视频| 人人妻人人澡欧美一区二区| 欧美高清性xxxxhd video| 国产成年人精品一区二区| 亚洲人成网站在线观看播放| 亚洲欧美日韩卡通动漫| 亚洲第一区二区三区不卡| 成人性生交大片免费视频hd| 美女大奶头视频| 欧美高清性xxxxhd video| 国产精品人妻久久久影院| 国产精品人妻久久久久久| 日本免费a在线| 成人亚洲欧美一区二区av| 精品欧美国产一区二区三| 蜜臀久久99精品久久宅男| 午夜精品国产一区二区电影 | 在线免费观看的www视频| 精品国内亚洲2022精品成人| 日韩成人av中文字幕在线观看 | 久99久视频精品免费| 日韩欧美精品v在线| 精品久久久久久久久久久久久| 女生性感内裤真人,穿戴方法视频| 日韩国内少妇激情av| 一边摸一边抽搐一进一小说| 久久久久久九九精品二区国产| 国产v大片淫在线免费观看| 久久久久久久久大av| 成人欧美大片| 日韩av在线大香蕉| 老熟妇仑乱视频hdxx| 国产69精品久久久久777片| 99热网站在线观看| 真人做人爱边吃奶动态| 亚洲欧美日韩卡通动漫| 成人av一区二区三区在线看| 免费人成在线观看视频色| 国产黄a三级三级三级人| 色5月婷婷丁香| 六月丁香七月| 亚洲精品亚洲一区二区| 欧美另类亚洲清纯唯美| 国产成年人精品一区二区| 97超碰精品成人国产| 少妇猛男粗大的猛烈进出视频 | 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利视频1000在线观看| 国产精品福利在线免费观看| 丝袜喷水一区| 少妇裸体淫交视频免费看高清| 高清午夜精品一区二区三区 | 欧美人与善性xxx| 日本一二三区视频观看| 蜜臀久久99精品久久宅男| 在线播放国产精品三级| 日韩制服骚丝袜av| 欧美丝袜亚洲另类| 菩萨蛮人人尽说江南好唐韦庄 | 欧美激情在线99| 三级国产精品欧美在线观看| aaaaa片日本免费| 亚洲欧美成人综合另类久久久 | 国产 一区 欧美 日韩| 99九九线精品视频在线观看视频| 老熟妇乱子伦视频在线观看| 一进一出抽搐动态| 一区福利在线观看| 久久久久久久亚洲中文字幕| 日韩成人av中文字幕在线观看 | 亚洲一区高清亚洲精品| 日韩三级伦理在线观看| 不卡视频在线观看欧美| 久久久国产成人免费| 少妇丰满av| 精品国内亚洲2022精品成人| 少妇猛男粗大的猛烈进出视频 | 欧美精品国产亚洲| 久久精品国产亚洲av涩爱 | 啦啦啦韩国在线观看视频| 丝袜美腿在线中文| 男女啪啪激烈高潮av片| 精品久久久久久成人av| 综合色av麻豆| 一级黄色大片毛片| 最新在线观看一区二区三区| 国产精品国产三级国产av玫瑰| 中国美白少妇内射xxxbb| 美女cb高潮喷水在线观看| 欧美日本亚洲视频在线播放| 午夜精品国产一区二区电影 | 色综合站精品国产| 国产高清不卡午夜福利| 色视频www国产| 精品午夜福利在线看| 搞女人的毛片| 国产人妻一区二区三区在| 欧美bdsm另类| 中文字幕av在线有码专区| 99riav亚洲国产免费| 国产一区二区三区av在线 | 久久精品人妻少妇| 中文字幕人妻熟人妻熟丝袜美| 超碰av人人做人人爽久久| av在线蜜桃| 亚洲av电影不卡..在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲精品一卡2卡三卡4卡5卡| 午夜激情欧美在线| 亚洲美女视频黄频| 久久久久免费精品人妻一区二区| 男人的好看免费观看在线视频| 国产精品无大码| 一级毛片aaaaaa免费看小| 亚洲婷婷狠狠爱综合网| 99国产极品粉嫩在线观看| 亚洲在线自拍视频| av视频在线观看入口| 亚洲精品日韩在线中文字幕 | 精品熟女少妇av免费看| 亚洲精品456在线播放app| 男女边吃奶边做爰视频| 国产精品亚洲一级av第二区| 亚洲一区高清亚洲精品| 国产成人福利小说| 国产高清不卡午夜福利| 亚洲精品456在线播放app| 99热只有精品国产| 亚洲av五月六月丁香网| 久久精品国产亚洲网站| 波多野结衣高清作品| 黄色日韩在线| 成人二区视频| 亚洲内射少妇av| 国产视频一区二区在线看| 晚上一个人看的免费电影| a级毛片a级免费在线| 永久网站在线| 日本三级黄在线观看| 深爱激情五月婷婷| 成年女人看的毛片在线观看| 热99re8久久精品国产|