• <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国片精品| 婷婷六月久久综合丁香| 正在播放国产对白刺激| av有码第一页| 国产精品亚洲一级av第二区| 午夜免费激情av| 欧美黄色片欧美黄色片| 亚洲专区中文字幕在线| 91成年电影在线观看| 最近最新中文字幕大全免费视频| 波多野结衣高清无吗| 精品国产国语对白av| 国产成人欧美在线观看| 真人一进一出gif抽搐免费| 亚洲九九香蕉| 国产不卡一卡二| 老汉色∧v一级毛片| 人人妻人人看人人澡| 亚洲狠狠婷婷综合久久图片| 久久久久久九九精品二区国产 | 国产精品国产高清国产av| 成人欧美大片| 日韩欧美免费精品| 女人高潮潮喷娇喘18禁视频| 亚洲第一电影网av| 母亲3免费完整高清在线观看| 丝袜在线中文字幕| 99久久国产精品久久久| 天堂动漫精品| 在线视频色国产色| 亚洲熟女毛片儿| 亚洲av中文字字幕乱码综合 | 欧美性长视频在线观看| 级片在线观看| 久久精品影院6| 人人妻人人澡欧美一区二区| 欧美亚洲日本最大视频资源| 又黄又爽又免费观看的视频| 久久 成人 亚洲| 在线观看日韩欧美| 日本 欧美在线| 久久久久国内视频| 欧美日韩一级在线毛片| 操出白浆在线播放| 男人操女人黄网站| 麻豆成人午夜福利视频| 宅男免费午夜| 日韩精品青青久久久久久| 国产精品免费视频内射| 91在线观看av| 久久 成人 亚洲| 国产人伦9x9x在线观看| 可以免费在线观看a视频的电影网站| 99久久综合精品五月天人人| 中亚洲国语对白在线视频| 欧美性猛交╳xxx乱大交人| 久久狼人影院| 99国产极品粉嫩在线观看| 成人亚洲精品一区在线观看| 欧美丝袜亚洲另类 | 久久久久久大精品| 国产精品一区二区三区四区久久 | 大型av网站在线播放| 大型av网站在线播放| av视频在线观看入口| 狠狠狠狠99中文字幕| 免费观看人在逋| 久久久国产精品麻豆| 极品教师在线免费播放| 亚洲国产欧洲综合997久久, | 精品久久蜜臀av无| 夜夜夜夜夜久久久久| 欧美日韩福利视频一区二区| 黄频高清免费视频| 黄色视频,在线免费观看| 久久久精品欧美日韩精品| 美女高潮喷水抽搐中文字幕| 亚洲熟妇中文字幕五十中出| 日韩欧美一区二区三区在线观看| 一区福利在线观看| 男女午夜视频在线观看| 人成视频在线观看免费观看| 国产亚洲av嫩草精品影院| 一级片免费观看大全| 两个人视频免费观看高清| 男女视频在线观看网站免费 | 亚洲片人在线观看| 久久久国产成人免费| 校园春色视频在线观看| av片东京热男人的天堂| 最近在线观看免费完整版| 日韩高清综合在线| 在线观看午夜福利视频| 青草久久国产| 热99re8久久精品国产| 久久午夜亚洲精品久久| 亚洲精华国产精华精| 久久久久国内视频| 99riav亚洲国产免费| 久久性视频一级片| 桃红色精品国产亚洲av| 俺也久久电影网| 制服诱惑二区| 十八禁人妻一区二区| 久9热在线精品视频| aaaaa片日本免费| cao死你这个sao货| 91成年电影在线观看| 亚洲国产精品sss在线观看| 亚洲国产精品999在线| 巨乳人妻的诱惑在线观看| 国产在线观看jvid| 国产真人三级小视频在线观看| 老汉色av国产亚洲站长工具| 亚洲avbb在线观看| 国产单亲对白刺激| 午夜久久久久精精品| 国产精品永久免费网站| 久久国产精品人妻蜜桃| 国产亚洲欧美精品永久| 怎么达到女性高潮| 精品久久久久久久久久久久久 | 成人精品一区二区免费| 成人国产综合亚洲| 老熟妇仑乱视频hdxx| 亚洲精品中文字幕在线视频| 91在线观看av| 在线十欧美十亚洲十日本专区| 久久精品国产综合久久久| 中文字幕久久专区| 两个人看的免费小视频| 我的亚洲天堂| 国产精品电影一区二区三区| 精品久久久久久久久久免费视频| 国产主播在线观看一区二区| 成人三级黄色视频| 中文字幕人成人乱码亚洲影| 国产在线观看jvid| 巨乳人妻的诱惑在线观看| 国产一区在线观看成人免费| 国产精品久久久av美女十八| 一区二区日韩欧美中文字幕| 一二三四社区在线视频社区8| 免费电影在线观看免费观看| 国产私拍福利视频在线观看| 国产三级黄色录像| 一个人免费在线观看的高清视频| 午夜影院日韩av| 亚洲第一欧美日韩一区二区三区| 可以免费在线观看a视频的电影网站| 久久久久久久久免费视频了| 中文字幕最新亚洲高清| 午夜免费鲁丝| 亚洲av成人一区二区三| 超碰成人久久| 久热爱精品视频在线9| 麻豆成人av在线观看| 亚洲中文字幕日韩| 国产精品久久久久久亚洲av鲁大| 国产又爽黄色视频| 91麻豆精品激情在线观看国产| 老汉色av国产亚洲站长工具| 久久香蕉国产精品| 午夜亚洲福利在线播放| 国产精品久久久久久人妻精品电影| 一级作爱视频免费观看| 99久久国产精品久久久| 国产高清videossex| 亚洲av中文字字幕乱码综合 | 国产亚洲av高清不卡| 波多野结衣高清无吗| 一级a爱片免费观看的视频| 在线永久观看黄色视频| 久久久久久大精品| 日本a在线网址| 久久久久国内视频| 久久久国产成人免费| 在线av久久热| 啦啦啦观看免费观看视频高清| 91老司机精品| 亚洲午夜理论影院| 黄色成人免费大全| 欧美亚洲日本最大视频资源| 婷婷精品国产亚洲av在线| 欧美午夜高清在线| 国产视频一区二区在线看| 国产亚洲av嫩草精品影院| 一进一出抽搐动态| 高清在线国产一区| 免费在线观看视频国产中文字幕亚洲| 日韩中文字幕欧美一区二区| 午夜免费成人在线视频| 国产精品免费一区二区三区在线| 久久精品影院6| 国内毛片毛片毛片毛片毛片| 1024视频免费在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| aaaaa片日本免费| 中文字幕最新亚洲高清| 高潮久久久久久久久久久不卡| 热re99久久国产66热| 曰老女人黄片| 国产一级毛片七仙女欲春2 | 国产三级黄色录像| 国产黄色小视频在线观看| 一二三四社区在线视频社区8| 久久精品影院6| 日本 av在线| 亚洲中文av在线| 久久久久国内视频| 99久久综合精品五月天人人| 成人精品一区二区免费| 亚洲一区二区三区色噜噜| 草草在线视频免费看| 女性生殖器流出的白浆| 亚洲精品久久国产高清桃花| 久久久国产成人精品二区| 久久青草综合色| 日本熟妇午夜| 999久久久国产精品视频| avwww免费| 久久中文字幕人妻熟女| 色播亚洲综合网| 久久人妻福利社区极品人妻图片| 亚洲黑人精品在线| 后天国语完整版免费观看| 久久久久久亚洲精品国产蜜桃av| 男女那种视频在线观看| 搞女人的毛片| 日本黄色视频三级网站网址| 国产av在哪里看| 亚洲精品美女久久久久99蜜臀| 老汉色∧v一级毛片| 亚洲午夜理论影院| av免费在线观看网站| 亚洲成人精品中文字幕电影| 免费av毛片视频| 少妇裸体淫交视频免费看高清 | 黑人巨大精品欧美一区二区mp4| 国产男靠女视频免费网站| 久久精品影院6| 精品国产乱码久久久久久男人| 99精品久久久久人妻精品| 国产高清videossex| 一本综合久久免费| 亚洲中文字幕日韩| 亚洲全国av大片| 他把我摸到了高潮在线观看| 极品教师在线免费播放| 香蕉国产在线看| 国产激情久久老熟女| 在线观看舔阴道视频| 999久久久精品免费观看国产| 国产一区在线观看成人免费| av有码第一页| 精品久久久久久久末码| 三级毛片av免费| 久久热在线av| 久久久久国产精品人妻aⅴ院| 成人国产一区最新在线观看| 亚洲,欧美精品.| 亚洲av第一区精品v没综合| 99精品久久久久人妻精品| 久久久久久久久免费视频了| 国产私拍福利视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 可以在线观看的亚洲视频| 精品久久久久久成人av| 欧美乱码精品一区二区三区| 国产午夜福利久久久久久| 香蕉av资源在线| 久久亚洲真实| 精品日产1卡2卡| 人人妻人人看人人澡| 久久婷婷成人综合色麻豆| 久久婷婷人人爽人人干人人爱| 黑丝袜美女国产一区| 亚洲人成电影免费在线| 欧美国产日韩亚洲一区| 美女国产高潮福利片在线看| 无遮挡黄片免费观看| 欧美日韩乱码在线| 99久久国产精品久久久| 国产人伦9x9x在线观看| 无限看片的www在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲精品美女久久久久99蜜臀| 精品欧美国产一区二区三| 成人18禁高潮啪啪吃奶动态图| 精品不卡国产一区二区三区| 久久久久久久久久黄片| 制服诱惑二区| 日韩欧美国产在线观看| 欧美亚洲日本最大视频资源| 亚洲av电影在线进入| 97超级碰碰碰精品色视频在线观看| 国产精品乱码一区二三区的特点| 久久亚洲精品不卡| 十八禁人妻一区二区| 国产1区2区3区精品| 美女午夜性视频免费| 久9热在线精品视频| av电影中文网址| 97人妻精品一区二区三区麻豆 | 欧美在线一区亚洲| 美女扒开内裤让男人捅视频| 亚洲午夜精品一区,二区,三区| 人人妻,人人澡人人爽秒播| 天天添夜夜摸| xxxwww97欧美| 国产伦在线观看视频一区| 99re在线观看精品视频| 欧美日韩福利视频一区二区| 宅男免费午夜| 两个人视频免费观看高清| 久久久久国产一级毛片高清牌| 看黄色毛片网站| 悠悠久久av| 国内揄拍国产精品人妻在线 | 一边摸一边抽搐一进一小说| 欧美另类亚洲清纯唯美| 欧美精品啪啪一区二区三区| 午夜福利18| 国内毛片毛片毛片毛片毛片| 午夜精品久久久久久毛片777| 久久国产精品人妻蜜桃| 女性被躁到高潮视频| 一a级毛片在线观看| 色播亚洲综合网| 国产v大片淫在线免费观看| 岛国在线观看网站| 国产三级在线视频| 午夜福利免费观看在线| 十分钟在线观看高清视频www| 成人永久免费在线观看视频| 成人三级黄色视频| 国产成人精品久久二区二区91| 亚洲欧美精品综合久久99| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美性猛交╳xxx乱大交人| 国内少妇人妻偷人精品xxx网站 | 中文字幕精品亚洲无线码一区 | 女人爽到高潮嗷嗷叫在线视频| 国产爱豆传媒在线观看 | 国产单亲对白刺激| 亚洲中文日韩欧美视频| 大型黄色视频在线免费观看| 国产三级黄色录像| 后天国语完整版免费观看| 淫妇啪啪啪对白视频| 国产午夜福利久久久久久| 亚洲国产欧洲综合997久久, | 亚洲精品一卡2卡三卡4卡5卡| 欧美激情 高清一区二区三区| 国产精华一区二区三区| 淫妇啪啪啪对白视频| 欧美+亚洲+日韩+国产| 亚洲狠狠婷婷综合久久图片| 精品久久蜜臀av无| 在线十欧美十亚洲十日本专区| 国内揄拍国产精品人妻在线 | 两人在一起打扑克的视频| 久久久久免费精品人妻一区二区 | videosex国产| 亚洲欧美激情综合另类| 中文字幕av电影在线播放| 啦啦啦免费观看视频1| 禁无遮挡网站| 看黄色毛片网站| 午夜免费成人在线视频| 国产高清有码在线观看视频 | 一二三四社区在线视频社区8| 中亚洲国语对白在线视频| 黄片播放在线免费| 精品少妇一区二区三区视频日本电影| 亚洲激情在线av| 精品国产乱码久久久久久男人| 老司机午夜福利在线观看视频| 美国免费a级毛片| 啦啦啦免费观看视频1| 亚洲五月天丁香| 国产精品永久免费网站| 欧美精品亚洲一区二区| 女同久久另类99精品国产91| 色综合婷婷激情| 搡老岳熟女国产| 禁无遮挡网站| 啦啦啦韩国在线观看视频| 久久国产乱子伦精品免费另类| 亚洲av熟女| 九色国产91popny在线| 午夜视频精品福利| 久久中文字幕一级| 精品久久久久久成人av| av免费在线观看网站| www.999成人在线观看| 国内精品久久久久久久电影| 曰老女人黄片| 日韩欧美三级三区| 99在线视频只有这里精品首页| 在线十欧美十亚洲十日本专区| 亚洲精品美女久久久久99蜜臀| 欧美一级毛片孕妇| 久久香蕉精品热| 亚洲无线在线观看| 可以免费在线观看a视频的电影网站| 欧美不卡视频在线免费观看 | 又黄又爽又免费观看的视频| 一区福利在线观看| 69av精品久久久久久| 看黄色毛片网站| 亚洲国产看品久久| 给我免费播放毛片高清在线观看| 成人精品一区二区免费| 一区福利在线观看| 黄片大片在线免费观看| 怎么达到女性高潮| 欧美又色又爽又黄视频| 国产一区二区三区在线臀色熟女| 精品午夜福利视频在线观看一区| 国产视频一区二区在线看| 90打野战视频偷拍视频| 国产av一区在线观看免费| 国产精品一区二区三区四区久久 | 99久久综合精品五月天人人| 又紧又爽又黄一区二区| 国产精品亚洲av一区麻豆| 午夜福利一区二区在线看| 亚洲人成电影免费在线| 亚洲国产日韩欧美精品在线观看 | 亚洲自拍偷在线| 国产成人精品久久二区二区91| 黑人欧美特级aaaaaa片| 日本五十路高清| 欧美黑人巨大hd| 怎么达到女性高潮| 丁香欧美五月| 国产熟女午夜一区二区三区| 免费在线观看视频国产中文字幕亚洲| 午夜福利在线在线| 国产精品综合久久久久久久免费| 香蕉丝袜av| 久久人妻av系列| 在线视频色国产色| 18禁黄网站禁片午夜丰满| 一本大道久久a久久精品| 久久狼人影院| 制服诱惑二区| 欧美大码av| 欧美三级亚洲精品| av有码第一页| 91麻豆av在线| 身体一侧抽搐| 国产亚洲精品久久久久5区| 亚洲国产中文字幕在线视频| 成人av一区二区三区在线看| 精品国产乱码久久久久久男人| 中文在线观看免费www的网站 | 亚洲午夜理论影院| 午夜影院日韩av| 在线观看午夜福利视频| 男男h啪啪无遮挡| www日本黄色视频网| 亚洲成av人片免费观看| netflix在线观看网站| 一进一出好大好爽视频| 国产成人一区二区三区免费视频网站| 岛国在线观看网站| 长腿黑丝高跟| 精品日产1卡2卡| 不卡一级毛片| 一夜夜www| 午夜免费鲁丝| 精华霜和精华液先用哪个| 操出白浆在线播放| 亚洲专区国产一区二区| 亚洲人成77777在线视频| 亚洲精品国产一区二区精华液| 欧美三级亚洲精品| 人成视频在线观看免费观看| 波多野结衣巨乳人妻| 香蕉av资源在线| 成人国语在线视频| 亚洲精品中文字幕在线视频| a级毛片在线看网站| 在线观看舔阴道视频| 两个人免费观看高清视频| 无人区码免费观看不卡| 国产精品亚洲一级av第二区| 国产一区二区三区在线臀色熟女| 在线播放国产精品三级| 国产av在哪里看| 18禁裸乳无遮挡免费网站照片 | 欧美三级亚洲精品| 91国产中文字幕| e午夜精品久久久久久久| 亚洲人成网站在线播放欧美日韩| 丝袜美腿诱惑在线| 人妻丰满熟妇av一区二区三区| 18禁黄网站禁片免费观看直播| 美女 人体艺术 gogo| 精品不卡国产一区二区三区| 久久精品国产亚洲av高清一级| 啦啦啦韩国在线观看视频| 久久中文看片网| 国产成人啪精品午夜网站| 国产男靠女视频免费网站| 在线观看免费午夜福利视频| 亚洲熟女毛片儿| 麻豆av在线久日| 又黄又粗又硬又大视频| 在线免费观看的www视频| 欧美一级毛片孕妇| 日本五十路高清| 好男人在线观看高清免费视频 | 国产av一区二区精品久久| 欧美成人免费av一区二区三区| 午夜影院日韩av| 三级毛片av免费| 亚洲专区中文字幕在线| 成人手机av| 亚洲人成网站在线播放欧美日韩| 制服丝袜大香蕉在线| 99国产精品99久久久久| 视频区欧美日本亚洲| 免费在线观看日本一区| 国产免费av片在线观看野外av| 亚洲第一欧美日韩一区二区三区| 悠悠久久av| 国产精品久久视频播放| 欧美性长视频在线观看| 美女高潮到喷水免费观看| 成人18禁在线播放| 丁香欧美五月| 国产精品久久电影中文字幕| 国产不卡一卡二| 日韩一卡2卡3卡4卡2021年| 身体一侧抽搐| 禁无遮挡网站| 日本黄色视频三级网站网址| cao死你这个sao货| 亚洲专区字幕在线| 露出奶头的视频| 黑人欧美特级aaaaaa片| 一区二区日韩欧美中文字幕| 国产精品美女特级片免费视频播放器 | 天堂影院成人在线观看| 久久人人精品亚洲av| 午夜老司机福利片| 午夜视频精品福利| 午夜精品在线福利| 亚洲五月色婷婷综合| 免费搜索国产男女视频| 亚洲专区字幕在线| 成人亚洲精品一区在线观看| 怎么达到女性高潮| 精品乱码久久久久久99久播| 国产成人啪精品午夜网站| 亚洲片人在线观看| 国产精品亚洲av一区麻豆| 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美98| 国产精品影院久久| 日本一区二区免费在线视频| 国产一区二区三区视频了| 久久亚洲真实| 757午夜福利合集在线观看| 欧美日韩亚洲综合一区二区三区_| 午夜免费鲁丝| 午夜福利一区二区在线看| 亚洲国产精品久久男人天堂| 久久精品aⅴ一区二区三区四区| 99久久综合精品五月天人人| 色在线成人网| av片东京热男人的天堂| 亚洲av成人不卡在线观看播放网| 午夜免费激情av| 麻豆成人av在线观看| 桃红色精品国产亚洲av| 中文亚洲av片在线观看爽| 视频在线观看一区二区三区| www日本在线高清视频| 91av网站免费观看| 天天一区二区日本电影三级| 欧美成人性av电影在线观看| 免费在线观看视频国产中文字幕亚洲| 久久中文字幕人妻熟女| 色综合婷婷激情| 天堂√8在线中文| 国产99久久九九免费精品| 国产欧美日韩精品亚洲av| 最新美女视频免费是黄的| 久久性视频一级片| 亚洲三区欧美一区| 欧美zozozo另类| 国产精品免费视频内射| 18美女黄网站色大片免费观看| 久久国产精品影院| 久久性视频一级片| 丝袜在线中文字幕| 国产麻豆成人av免费视频| 国产精品久久久久久亚洲av鲁大| 亚洲欧美激情综合另类| 91成人精品电影| 女人被狂操c到高潮| 欧美av亚洲av综合av国产av| 亚洲av电影不卡..在线观看| 国产高清激情床上av| 欧美中文综合在线视频| 国产精品香港三级国产av潘金莲| 亚洲精品国产一区二区精华液|