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

    唐山老震區(qū)現(xiàn)今地震活動(dòng)特征

    2024-05-29 07:13:50周依曹筠馬旭東單連君
    地震研究 2024年3期
    關(guān)鍵詞:小震震區(qū)余震

    周依 曹筠 馬旭東 單連君

    摘要:為分析唐山老震區(qū)地震活動(dòng)背景,對(duì)該區(qū)域現(xiàn)今地震活動(dòng)特征進(jìn)行研究。老震區(qū)現(xiàn)今地震活動(dòng)與1976年唐山MS7.8地震早期余震密集區(qū)域存在一定差異,結(jié)合區(qū)域背景地震活動(dòng)統(tǒng)計(jì)結(jié)果,認(rèn)為該區(qū)域地震活動(dòng)受唐山地震序列影響逐漸減弱,逐步趨近于背景地震。地震分布呈現(xiàn)出明顯的3分區(qū)差異,主要集中于中區(qū),東區(qū)次之,西區(qū)最弱。中等地震的空間分布特征顯示,地震活動(dòng)存在東、中、西區(qū)交替發(fā)震的規(guī)律。目前,唐山老震區(qū)構(gòu)造應(yīng)力積累加速,存在地震活動(dòng)水平由平靜向活躍轉(zhuǎn)變的可能,需要重點(diǎn)關(guān)注唐山斷裂帶北段以及灤縣—樂(lè)亭斷裂和盧龍斷裂的交匯區(qū)域。

    關(guān)鍵詞:唐山地震;地震活動(dòng)特征;背景地震;遷移規(guī)律

    中圖分類(lèi)號(hào):P315.72文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):1000-0666(2024)03-0430-11

    doi:10.20015/j.cnki.ISSN1000-0666.2024.0046

    0引言

    1976年7月28日3時(shí)42分,唐山發(fā)生MS7.8地震,15個(gè)小時(shí)之后在其東部的灤縣(現(xiàn)灤州市)發(fā)生MS7.1地震,3個(gè)半月后在其西部的寧河再次發(fā)生MS6.9地震,主震及兩次強(qiáng)余震之后形成3個(gè)明顯的余震區(qū),震中最高烈度XI度,造成了嚴(yán)重的災(zāi)害損失(國(guó)家地震局《一九七六年唐山地震》編輯組,1982;孫麗娜等,2021)。唐山地震序列復(fù)雜且豐富,由7.8級(jí)主震、兩次7級(jí)左右強(qiáng)余震和一系列5級(jí)以上余震構(gòu)成。隨著時(shí)間推移,序列總體呈衰減特征,衰減過(guò)程中雖有階段性起伏,但是符合強(qiáng)震的衰減規(guī)律(張躍剛等,2013;郭蕾等,2021)。

    截至2022年,唐山地震已過(guò)去46年,該地區(qū)現(xiàn)今小震活動(dòng)仍十分頻繁。唐山老震區(qū)巖石圈黏彈性松弛模擬結(jié)果顯示,1976年唐山地震發(fā)生15年內(nèi)庫(kù)倫應(yīng)力變化顯著,而15年之后庫(kù)倫應(yīng)力呈現(xiàn)緩慢的穩(wěn)定增加狀態(tài),顯示區(qū)域應(yīng)力調(diào)整已基本穩(wěn)定,因此15年左右的時(shí)間尺度可能是唐山余震活動(dòng)的一個(gè)時(shí)間分界點(diǎn)(朱琳等,2021)。王想等(2012)認(rèn)為唐山地區(qū)1992年進(jìn)入余震活動(dòng)晚期階段之后,受唐山MS7.8主震和強(qiáng)余震的控制減弱,地震活動(dòng)趨于穩(wěn)定狀態(tài)。1996年以來(lái)唐山老震區(qū)中等地震活動(dòng)顯示,唐山MS7.8地震震中附近區(qū)域基本沒(méi)有4級(jí)以上地震,余震活動(dòng)可能已經(jīng)不受主震引起的局部應(yīng)力狀態(tài)變化所制約(張素欣等,2020)。鄭建常等(2014,2022)通過(guò)對(duì)唐山地震序列不同時(shí)段進(jìn)行叢集地震識(shí)別,認(rèn)為2011年之后唐山地震序列受主震和較大余震影響的事件越來(lái)越少,更多小震已表現(xiàn)為彼此獨(dú)立的隨機(jī)事件。以上研究均表明,唐山余震活動(dòng)進(jìn)入晚期階段后,地震活動(dòng)受主震影響逐漸減小,目前該區(qū)域發(fā)生的地震可看作背景地震。

    唐山MS7.8地震屬于大陸板塊內(nèi)部地震(Robinson,Zhou,2005)。[HJ2mm]研究顯示,板內(nèi)地震衰減較慢,余震持續(xù)時(shí)間可達(dá)數(shù)百年(Stein,Liu,2009)。仲秋和史保平(2012)根據(jù)對(duì)數(shù)線性回歸和余震觸發(fā)理論計(jì)算,得到唐山地震的余震序列還將持續(xù)近100年的結(jié)論。唐山老震區(qū)的GPS測(cè)量結(jié)果顯示,唐山地區(qū)的地殼應(yīng)變率高于華北地區(qū)的背景值,王輝等(2013)認(rèn)為這一現(xiàn)象表明該地區(qū)現(xiàn)今仍受1976年唐山地震序列影響。關(guān)于唐山老震區(qū)現(xiàn)今地震活動(dòng)是否為唐山MS7.8地震余震的討論始終沒(méi)有明確的結(jié)論。唐山地區(qū)如今經(jīng)過(guò)40多年的恢復(fù)重建,經(jīng)濟(jì)發(fā)達(dá)、人口稠密,其所處的首都圈地區(qū),一直是我國(guó)地震重點(diǎn)監(jiān)視防御區(qū)(孫啟凱等,2017)。因此,開(kāi)展唐山老震區(qū)現(xiàn)今地震活動(dòng)特征研究,分析其是否為該區(qū)域的背景地震活動(dòng),對(duì)于進(jìn)行唐山地區(qū)地震風(fēng)險(xiǎn)評(píng)估以及助力當(dāng)?shù)胤勒饻p災(zāi)建設(shè)具有一定的現(xiàn)實(shí)意義。

    1區(qū)域地質(zhì)構(gòu)造背景與資料選取

    1.1區(qū)域地質(zhì)構(gòu)造背景

    唐山地區(qū)位于張渤地震帶和華北平原地震帶交匯處,受NEE向豐臺(tái)—野雞坨斷裂、寧河—昌黎斷裂和NW向?yàn)纯h—樂(lè)亭斷裂、薊運(yùn)河斷裂的切割,形成了一個(gè)NEE向菱形塊體(虢順民等,1977)。菱形塊體內(nèi)部蓋層發(fā)育有唐山斷裂帶,西南界的薊運(yùn)河斷裂與塊體西南發(fā)育的NNE向滄東斷裂相交,東邊界的灤縣—樂(lè)亭斷裂與NNE向盧龍斷裂相交(圖1)。塊體內(nèi)部的唐山斷裂帶是一條走向NE、具右旋走滑特征的褶皺斷裂帶,由陡河斷裂、巍山—長(zhǎng)山南坡斷裂、唐山—古冶斷裂等數(shù)條次級(jí)斷裂組成,1976年唐山MS7.8地震序列主要沿唐山斷裂帶展布(徐杰等,2012)。本文研究區(qū)唐山老震區(qū)主要指唐山地震序列的小震密集區(qū)(張素欣等,2017),其深部地球物理探測(cè)顯示,在唐山斷裂帶的正下方存在一條傾角陡立,具有走滑性質(zhì)的深大斷裂,被稱(chēng)為“唐山地殼深斷裂”。該斷裂錯(cuò)斷了上、下地殼分界和莫霍面,為幔源物質(zhì)的侵入提供了通道(劉保金等,2011)。據(jù)此推測(cè),唐山地區(qū)深、淺部耦合的地殼構(gòu)造系統(tǒng)可能是該區(qū)域地震活動(dòng)頻繁的原因(李赫等,2020)。

    1.2資料選取

    1992年以來(lái),唐山余震活動(dòng)進(jìn)入晚期階段(王想等,2012)。據(jù)此,本文基于河北地震臺(tái)網(wǎng)

    資料,整理了1992—2021年唐山地區(qū)(38.8°~40.5°N,117°~120°E)地震目錄,對(duì)該區(qū)域現(xiàn)今地震活動(dòng)特征進(jìn)行分析。為保證分析結(jié)果的可靠性,采用震級(jí)-序號(hào)法和多統(tǒng)計(jì)方法(蔣長(zhǎng)勝等,2013)對(duì)最小完整性震級(jí)MC進(jìn)行評(píng)估。其中,多統(tǒng)計(jì)方法是使用最大曲率法(MAXC)和擬合優(yōu)度法(GFT-90%、GFT-95%)的一種綜合評(píng)估方法。如圖2所示,兩種方法計(jì)算得到的結(jié)果較為一致,1992年以來(lái),唐山老震區(qū)最小完整性震級(jí)MC為ML2.0左右;2008年以來(lái),隨著“十五”數(shù)字地震觀測(cè)網(wǎng)絡(luò)項(xiàng)目投入運(yùn)行,研究區(qū)的地震監(jiān)測(cè)能力顯著提升,MC已低于ML1.5;2015年以后,得益于首都圈測(cè)震臺(tái)站升級(jí)改造項(xiàng)目的完成,研究區(qū)地震監(jiān)測(cè)能力進(jìn)一步提升,MC達(dá)到ML1.0。因此,本文分析唐山老震區(qū)現(xiàn)今地震活動(dòng)特征時(shí),將地震起始震級(jí)設(shè)定為ML2.0。

    2小震密集分布特征分析

    根據(jù)1992年以來(lái)唐山老震區(qū)最小完整性震級(jí)MC,設(shè)置起始震級(jí)為ML2.0,繪制區(qū)域地震震中分布圖。如圖3a所示,地震表現(xiàn)出明顯的密集分布特征,主要沿唐山斷裂帶和盧龍斷裂呈NE向密集分布,沿著薊運(yùn)河斷裂也存在少量地震呈NW向帶狀分布。在盧龍斷裂和灤縣—樂(lè)亭斷裂交匯區(qū)域附近,地震密集活動(dòng)出現(xiàn)“丁”字形轉(zhuǎn)折。為進(jìn)一步分析地震活動(dòng)的時(shí)空分布特征,沿地震密集區(qū)繪制了AB軸,長(zhǎng)軸方向NE。在分析小震時(shí)空分布圖像時(shí),設(shè)定A點(diǎn)為研究區(qū)軸線原點(diǎn),根據(jù)距離軸線20 km范圍內(nèi)地震相對(duì)原點(diǎn)的距離,得到唐山地區(qū)小震時(shí)空演化圖像(圖3b)。地震密集區(qū)集中在距離A點(diǎn)50~110 km處,該區(qū)域主要為唐山斷裂帶北段和盧龍斷裂南段,大部分ML4.0以上地震都發(fā)生在該區(qū)域,其中包括所有ML5.0以上地震,其余的ML4.0以上地震主要分布在薊運(yùn)河斷裂與唐山斷裂帶交匯附近。2004年1月20日昌黎ML5.0地震后,雖然小震活動(dòng)仍然集中于該區(qū)域,但是地震密集范圍明顯收縮,更為集中在唐山斷裂帶北端和盧龍斷裂與灤縣—樂(lè)亭斷裂的交匯區(qū)域,使得在唐山斷裂帶和灤縣—樂(lè)亭斷裂交而不匯的位置出現(xiàn)了一條清晰的地震分界線。

    為了更直觀地反映唐山老震區(qū)小震分布情況,區(qū)分地震活動(dòng)相對(duì)密集地區(qū)與稀疏地區(qū),本文基于小震震中分布繪制唐山地區(qū)地震密集分布圖。對(duì)研究區(qū)域進(jìn)行0.1°×0.1°的網(wǎng)格劃分,統(tǒng)計(jì)網(wǎng)格節(jié)點(diǎn)附近10 km范圍內(nèi)的地震次數(shù)并對(duì)其進(jìn)行歸一化處理,由此得到的密集程度值能夠直觀描述研究區(qū)域的地震密集程度和密集區(qū)域。從小震密集分布特征來(lái)看,唐山老震區(qū)現(xiàn)今地震活動(dòng)存在明顯差異,呈現(xiàn)空間分布不均勻特征。結(jié)合區(qū)域構(gòu)造環(huán)境并參考前人的分區(qū)規(guī)則(楊理華,趙喜柱,1985;張宏志等,2008),本文將唐山老震區(qū)劃分為東區(qū)、中區(qū)和西區(qū),其中灤縣—樂(lè)亭斷裂及以東為東區(qū),唐山斷裂帶所在的菱形塊體內(nèi)部為中區(qū),薊運(yùn)河斷裂及以西為西區(qū)。如圖4a所示,1992年以來(lái),中區(qū)地震密集程度尤為突出,地震密集區(qū)域主要集中在唐山斷裂帶北段,其南段地震分布相對(duì)稀疏;東區(qū)地震密集程度次之,主要集中在灤縣—樂(lè)亭斷裂與盧龍斷裂交匯區(qū)域,盧龍斷裂北段有少量地震聚集;西區(qū)地震活動(dòng)沿著薊運(yùn)河斷裂呈NW向分布,地震密集程度較弱。將1976年唐山MS7.8主震及兩次強(qiáng)余震的震中位置投影到地震密集分布圖上發(fā)現(xiàn),唐山MS7.8和灤縣MS7.1地震震中均位于現(xiàn)今地震密集區(qū)域的邊緣附近,寧河MS6.9地震震中位于現(xiàn)今地震活動(dòng)較弱區(qū)域。而1992年以前的余震密集分布圖顯示(圖4b),唐山MS7.8和灤縣MS7.1地震震中均處于地震密集程度較高區(qū)域,且寧河MS6.9地震震中也位于地震密集區(qū)域的邊緣。與1992年以前相比,現(xiàn)今地震密集區(qū)域收縮,偏離3次強(qiáng)震震中,可能預(yù)示區(qū)域地震受唐山MS7.8強(qiáng)震序列影響減弱。

    3中等地震時(shí)空分布特征分析

    3.1時(shí)間叢集特征

    區(qū)域地震活動(dòng)時(shí)間分布特征研究是在一定統(tǒng)計(jì)區(qū)內(nèi)(地震帶或地震區(qū)內(nèi))進(jìn)行的,描述在時(shí)間進(jìn)程中不同時(shí)間段的地震活動(dòng)特征。1992年以來(lái)唐山老震區(qū)ML4.0以上地震累積應(yīng)變能量釋放曲線整體表現(xiàn)比較平緩,未出現(xiàn)較大波動(dòng),說(shuō)明整體地震活動(dòng)比較平穩(wěn)。從細(xì)節(jié)上看,累積應(yīng)變能量釋放曲線在2004年和2019年出現(xiàn)了小的臺(tái)階,參考1992年以來(lái)唐山老震區(qū)ML4.0以上地震的M-t圖,可以將該區(qū)域現(xiàn)今中等地震活動(dòng)大致分為3個(gè)階段(圖5)。第1階段1992—2004年,累積應(yīng)變能量釋放用J來(lái)表示,其速率為11.2 J/a,ML4.0以上地震平均年頻次為1.7次。期間發(fā)生了3次ML5.0以上地震,最大地震為1995年10月6日唐山ML5.4。該階段地震發(fā)震時(shí)間間隔多集中在0.5~1 a,最長(zhǎng)發(fā)震間隔為2 a,期間多次出現(xiàn)短期(3個(gè)月內(nèi))地震成對(duì)發(fā)生。第2階段2005—2018年,累積應(yīng)變能量釋放速率和地震平均年頻次較第1階段有所降低,分別為6.0 J/a和0.9次。期間發(fā)生ML5.0以上地震1次,為2012年5月28日唐山ML5.1。該階段地震存在2 a左右的準(zhǔn)周期特征,并且“雙震”特征更加顯著。但2015年9月14日昌黎ML4.8地震擾亂了2 a左右的準(zhǔn)周期特征,此后發(fā)震時(shí)間間隔縮短,此時(shí)可能處于這1階段末期的剩余應(yīng)變釋放調(diào)整階段。2016年9月10日唐山ML4.6地震后,唐山老震區(qū)進(jìn)入了長(zhǎng)達(dá)3.2 a的ML4地震平靜,直至2019年12月5日豐南ML4.9地震的發(fā)生,預(yù)示著進(jìn)入了地震活動(dòng)的第3階段。目前處于第3階段早期,共計(jì)發(fā)生3次ML4.0以上地震,最大地震為2020年7月12日古冶ML5.4地震,該階段地震發(fā)震時(shí)間間隔保持在0.5~1.0 a,目前累積應(yīng)變能量釋放速率為10.5 J/a。

    前2個(gè)階段的持續(xù)時(shí)間均在10 a左右,第2階段累積應(yīng)變能量釋放速率、地震平均年頻次和最大地震震級(jí)均低于第1階段,符合老震區(qū)地震衰減規(guī)律。但是從第3階段唐山老震區(qū)目前發(fā)生的中等地震活動(dòng)情況來(lái)看,累積應(yīng)變能量釋放速率和最大地震震級(jí)均有所回升,ML4.0以上地震的發(fā)震時(shí)間間隔縮短,存在地震活動(dòng)水平由平靜向活躍轉(zhuǎn)變的可能,與衰減趨勢(shì)相悖,可能意味著該區(qū)域地震的發(fā)生受唐山MS7.8強(qiáng)震序列約束進(jìn)一步減弱。

    3.2空間遷移特征

    1976年唐山MS7.8主震及兩次強(qiáng)余震依次發(fā)生于3分區(qū)中的中、東和西區(qū),存在空間遷移特征(圖4)。眾多專(zhuān)家學(xué)者已對(duì)其成因及動(dòng)力機(jī)制進(jìn)行了闡述(Butler et al,1979;國(guó)家地震局《一九七六年唐山地震》編輯組,1982;聞學(xué)澤,馬勝利,2006)。本文對(duì)唐山老震區(qū)現(xiàn)今中等地震空間分布特征進(jìn)行分析時(shí),發(fā)現(xiàn)這一特征同樣顯著。

    將1992年以來(lái)唐山老震區(qū)ML4.0以上地震分組討論,若以發(fā)震時(shí)間間隔超過(guò)2.0 a作為分組條件,則大致可以分為5組,3個(gè)階段分別存在2組、2組和1組地震活動(dòng)。從地震遷移曲線上看(圖6),存在東、中、西3個(gè)分區(qū)交替發(fā)震的特點(diǎn)。5組地震活動(dòng)中有3組的始發(fā)地震出現(xiàn)在中區(qū),即遷移過(guò)程大多從中區(qū)開(kāi)始,而結(jié)束地震多發(fā)生在東區(qū)。從遷移特征來(lái)看,雖未發(fā)現(xiàn)較明顯的優(yōu)勢(shì)遷移方向,但遷移過(guò)程中存在一定的有序性,當(dāng)?shù)卣鸹顒?dòng)按經(jīng)度從東向西或從西向東活動(dòng)時(shí),會(huì)在3個(gè)分區(qū)依次進(jìn)行遷移,每個(gè)遷移過(guò)程中區(qū)都會(huì)參加。在19次不同分區(qū)的地震遷移過(guò)程中,僅有3次沒(méi)有經(jīng)過(guò)中區(qū),其中有2次從東區(qū)直接跳躍到西區(qū)(第1、4組),還有1次則從西區(qū)跳躍到東區(qū)(第1組)。由此可以將中區(qū)作為唐山老震區(qū)地震活動(dòng)的窗口。研究表明,在東、中、西3個(gè)分區(qū)均存在應(yīng)力集中點(diǎn)(馬瑾等,1980;劉蒲雄,呂曉健,2011),當(dāng)區(qū)域應(yīng)力積累到一定程度后,在其中一個(gè)應(yīng)力集中點(diǎn)進(jìn)行應(yīng)力釋放。1976年唐山7.8級(jí)地震后,由于老震區(qū)地質(zhì)介質(zhì)相對(duì)破碎,介質(zhì)強(qiáng)度不易積累更大能量(Li et al,2007),較難阻止地震波能量的傳播,初始破裂的沖擊波誘發(fā)了其余兩個(gè)應(yīng)力集中點(diǎn)的應(yīng)力釋放,導(dǎo)致發(fā)震。3個(gè)分區(qū)中等地震空間遷移的現(xiàn)象,就是在3個(gè)應(yīng)力集中點(diǎn)依次應(yīng)力釋放的過(guò)程。

    若以遍歷東、中、西區(qū)為一重遷移,再一次遍歷為二重遷移,則前2組地震活動(dòng)存在二重遷移,后3組地震活動(dòng)僅存在一重遷移。5組地震活動(dòng)存在的7次遷移過(guò)程中,遷移路徑為“西→中→東”模式的情況最多(圖7a、b、e),“東→中→西”模式次之(圖7b、c)。還有2次遷移過(guò)程為跨區(qū)跳躍的模式:第1組地震活動(dòng)的第一重遷移過(guò)程中,天津發(fā)生ML4.9地震后,直接在東區(qū)發(fā)生了灤縣ML4.0地震(圖7a);第4組地震活動(dòng)遷移過(guò)程中,東區(qū)發(fā)生昌黎ML4.7地震后,接著西區(qū)發(fā)生豐南ML4.0地震(圖7d)。但是從統(tǒng)計(jì)來(lái)看,跨區(qū)跳躍的情況較少發(fā)生,遷移過(guò)程主要遵循3個(gè)分區(qū)依次遷移的規(guī)律。

    研究時(shí)段內(nèi)最后一次ML4.0以上地震為2021年4月16日灤州ML4.8,位于東區(qū)。自此,第5組地震活動(dòng)已經(jīng)完成了一重遷移,遷移模式為“西→中→東”。參考前幾次該模式下的遷移過(guò)程,在第1組地震活動(dòng)的第二重遷移過(guò)程中,中區(qū)依次發(fā)生了唐山ML4.0和ML5.4地震后,再往東區(qū)遷移;第2組地震活動(dòng)的第二重遷移過(guò)程,在東區(qū)接連發(fā)生了昌黎ML4.1和ML5.0兩次結(jié)束地震。考慮到“西→中→東”遷移模式下“雙震”現(xiàn)象的存在,第5組地震遷移過(guò)程有可能還沒(méi)有結(jié)束,即下一次ML4.0以上地震存在發(fā)生在東區(qū)的可能。

    4討論

    4.1現(xiàn)今地震活動(dòng)性質(zhì)判定

    1976年唐山MS7.8地震早期余震主要分布在主震及兩次強(qiáng)余震附近(Shedlock et al,1987;張素欣等,2017,2020),且1992年前、后的余震密集分布也表現(xiàn)出不同特征(圖4)。許忠淮等(1983,1985)利用唐山地震前10年的小震P波初動(dòng)數(shù)據(jù),得到震前區(qū)域主壓應(yīng)力方位約為N60°E,而震后早期表現(xiàn)為近EW方向,順時(shí)針旋轉(zhuǎn)了約30°。之后,唐山老震區(qū)的主壓應(yīng)力方向長(zhǎng)期處在近EW向(李瑞莎等,2008)。但是,2002年以來(lái)的地殼應(yīng)力場(chǎng)反演結(jié)果和GPS應(yīng)變場(chǎng)顯示,該區(qū)域目前主要受NEE向的主壓應(yīng)變控制(陳長(zhǎng)云,2016;王曉山等,2020)。表明唐山老震區(qū)現(xiàn)今主壓應(yīng)力方向較1976年唐山地震序列早期可能發(fā)生了調(diào)整,趨于恢復(fù)到震前水平。而唐山老震區(qū)現(xiàn)今地震活動(dòng)密集區(qū)域的變化,也是區(qū)域應(yīng)力場(chǎng)調(diào)整下的一種間接表現(xiàn)。

    從唐山MS7.8地震序列時(shí)序變化來(lái)看(圖8),地震頻度隨時(shí)間呈指數(shù)衰減,地震強(qiáng)度逐漸減弱但存在起伏變化。在序列早期,余震衰減較快,2年內(nèi)最大余震衰減至ML5.0左右,85%左右的余震活動(dòng)集中在該時(shí)間段內(nèi)。其后余震年頻次減少,地震強(qiáng)度較為穩(wěn)定,直至1991年5月30日發(fā)生唐山ML5.6地震,地震強(qiáng)度出現(xiàn)起伏。此后老震區(qū)ML4.0以上地震年頻次進(jìn)一步減少,地震活動(dòng)呈現(xiàn)階段性變化,與余震衰減趨勢(shì)相悖,可能意味著唐山老震區(qū)現(xiàn)今地震活動(dòng)受1976年唐山MS7.8地震序列影響較弱。通常定義余震活動(dòng)持續(xù)時(shí)間為余震活動(dòng)減弱到與背景地震活動(dòng)水平相當(dāng)?shù)臅r(shí)間(蔣海昆等,2006),1992年以來(lái)區(qū)域ML4.0以上地震平均年頻次為1.2次,與震前區(qū)域地震平均年頻次(1.3次)基本相當(dāng)(圖8b)。推測(cè)唐山老震區(qū)現(xiàn)今地震活動(dòng)趨近于震前活動(dòng)水平,發(fā)生的地震可能為背景地震。

    4.2地震活動(dòng)趨勢(shì)分析

    唐山老震區(qū)現(xiàn)今中等地震活動(dòng)存在3個(gè)階段,目前處于第3階段早期。從累積應(yīng)變能量釋放速率、地震平均年頻次和最大地震震級(jí)的計(jì)算及統(tǒng)計(jì)結(jié)果來(lái)看,第2階段均低于第1階段,地震活動(dòng)由活躍轉(zhuǎn)為平靜,符合老震區(qū)地震衰減規(guī)律。小震密集分布顯示,第2階段之后(即2004年之后),唐山老震區(qū)密集活動(dòng)范圍收縮。因此,大致可將2004年看作該區(qū)域地震活動(dòng)由活躍轉(zhuǎn)為平靜的一個(gè)時(shí)間節(jié)點(diǎn)。2004年以來(lái)唐山老震區(qū)中等地震活動(dòng)形態(tài)呈現(xiàn)出“平靜—叢集”的周期性變化。2019年以后,累積應(yīng)變能量釋放曲線顯示區(qū)域構(gòu)造應(yīng)力積累加速(圖5b),存在地震活動(dòng)水平由平靜向活躍轉(zhuǎn)變的可能。結(jié)合小震密集分布特征,唐山斷裂帶北段以及灤縣—樂(lè)亭斷裂和盧龍斷裂的交匯區(qū)域可能是未來(lái)中強(qiáng)地震活躍的場(chǎng)所。

    唐山老震區(qū)3個(gè)分區(qū)的地震活動(dòng)差異顯著,并且存在中等地震東、中、西區(qū)交替發(fā)震的規(guī)律。依據(jù)對(duì)中等地震遷移過(guò)程的統(tǒng)計(jì),遷移過(guò)程大多從中區(qū)開(kāi)始,即構(gòu)造應(yīng)力多數(shù)從中區(qū)積累釋放。并且中區(qū)也是小震活動(dòng)最為密集的區(qū)域,因此可以將中區(qū)作為區(qū)域地震活動(dòng)的窗口,中區(qū)開(kāi)始活動(dòng),意味著唐山老震區(qū)可能開(kāi)啟活動(dòng)。

    本文對(duì)唐山老震區(qū)現(xiàn)今中等地震遷移現(xiàn)象的分析,結(jié)果可以嘗試應(yīng)用于區(qū)域地震遷移趨勢(shì)分析,對(duì)未來(lái)中強(qiáng)地震發(fā)震位置的判定提供一定參考。根據(jù)遷移模式分析,認(rèn)為2021年4月16日灤州ML4.8地震是一重遷移過(guò)程的結(jié)束地震,但考慮到以往震例中雙震情況的出現(xiàn),認(rèn)為東區(qū)還有可能發(fā)生一次ML4.0以上地震。2022年5月21日盧龍ML4.1地震的發(fā)生也驗(yàn)證了這一預(yù)測(cè)結(jié)論。

    5結(jié)論

    本文根據(jù)唐山老震區(qū)現(xiàn)今地震的最小完整性震級(jí),利用1992—2021年的小震目錄,對(duì)該區(qū)域現(xiàn)今地震活動(dòng)特征進(jìn)行了分析,得到以下結(jié)論:

    (1)1992年以來(lái)唐山老震區(qū)小震密集分布于唐山斷裂帶北段以及灤縣—樂(lè)亭斷裂和盧龍斷裂的交匯區(qū)域,與1976年唐山MS7.8地震早期余震密集區(qū)域存在一定差異,可能是區(qū)域地震活動(dòng)受強(qiáng)震序列影響減弱的表現(xiàn)。唐山MS7.8地震前、后老震區(qū)中等地震頻次統(tǒng)計(jì)結(jié)果顯示,現(xiàn)今中等地震年頻次與震前基本相當(dāng),區(qū)域地震活動(dòng)趨近于震前活動(dòng)水平,由此認(rèn)為該區(qū)域所發(fā)生地震逐步趨近于背景地震。

    (2)唐山老震區(qū)現(xiàn)今地震分布3個(gè)分區(qū)的差異顯著,主要集中于中區(qū),東區(qū)次之,西區(qū)最弱。中等地震活動(dòng)存在較為明顯的遷移規(guī)律,即地震在東、中、西區(qū)交替發(fā)震。將這一規(guī)律應(yīng)用于區(qū)域地震遷移趨勢(shì)預(yù)測(cè),可以對(duì)未來(lái)中強(qiáng)地震發(fā)震位置的判定提供一定的參考。

    (3)唐山老震區(qū)現(xiàn)今地震活動(dòng)呈現(xiàn)階段性變化,2019年以來(lái),區(qū)域構(gòu)造應(yīng)力積累加速,地震活動(dòng)水平存在由平靜向活躍轉(zhuǎn)變的可能。現(xiàn)今小震活動(dòng)密集分布于唐山斷裂帶北段以及灤縣—樂(lè)亭斷裂和盧龍斷裂的交匯區(qū)域,上述地區(qū)可能成為未來(lái)中強(qiáng)地震活躍的場(chǎng)所。[HTK]

    文中圖件使用GMT軟件繪制,在此表示感謝。

    參考文獻(xiàn):

    陳長(zhǎng)云.2016.張家口—渤海斷裂帶分段運(yùn)動(dòng)變形特征分析[J].地震,36(1):1-11.Chen C Y.2016.Characteristics of segmentary motion and deformation along the Zhangjiakou-Bohai fault[J].Earthquake,36(1):1-11.(in Chinese)

    郭蕾,孫麗娜,周依,等.2021.唐山地區(qū)近十年地震的重定位研究[J].地震學(xué)報(bào),43(6):706-715.Guo L,Sun L N,Zhou Y,et al.2021.Relocation of earthquakes in Tangshan area in the recent decade[J].Acta Seismologica Sinica,43(6):706-715.(in Chinese)

    國(guó)家地震局《一九七六年唐山地震》編輯組.1982.一九七六年唐山地震[M].北京:地震出版社.Editorial Group of The Tangshan Earthquake of 1976,State Seismological Bureau.1982.Tangshan Earthquake of 1976[M].Beijing:Seismological Press.(in Chinese)

    虢順民,李志義,程紹平,等.1977.唐山地震區(qū)域構(gòu)造背景和發(fā)震模式的討論[J].地質(zhì)科學(xué),12(4):305-321.Guo S M,Li Z Y,Cheng S P,et al.1977.Discussion on the regional structural background and the seismogenic model of the Tangshan earthquake[J].Scientia Geologica Sinica,12(4):305-321.(in Chinese)

    蔣長(zhǎng)勝,吳忠良,莊建倉(cāng).2013.地震的“序列歸屬”問(wèn)題與ETAS模型——以唐山序列為例[J].地球物理學(xué)報(bào),56(9):2971-2981.Jiang C S,Wu Z L,Zhuang J C.2013.ETAS model applied to the Earthquake-Sequence Association(ESA)problem:the Tangshan sequence[J].Chinese Journal of Geophysics,56(9):2971-2981.(in Chinese)

    蔣海昆,曲延軍,李永莉,等.2006.中國(guó)大陸中強(qiáng)地震余震序列的部分統(tǒng)計(jì)特征[J].地球物理學(xué)報(bào),49(4):1110-1117.Jiang H K,Qu Y J,Li Y L,et al.2006.Some statistic features of aftershock sequences in Chinese mainland[J].Chinese Journal of Geophysics,49(4):1110-1117.(in Chinese)

    李赫,董一兵,王熠熙,等.2020.唐山地區(qū)地震成因初探——以豐南2010年M4.1和2019年M4.5地震序列為例[J].地球物理學(xué)進(jìn)展,35(4):1276-1283.Li H,Dong Y B,Wang Y X,et al.2020.Discussion on the seismogenesis of Tangshan area:take the sequence Fengnan M4.1 in 2010 and M4.5 in 2019 as an example[J].Progress in Geophysics,35(4):1276-1283.(in Chinese)

    李瑞莎,崔效鋒,刁桂苓,等.2008.華北北部地區(qū)現(xiàn)今應(yīng)力場(chǎng)時(shí)空變化特征研究[J].地震學(xué)報(bào),30(6):570-580.Li R S,Cui X F,Diao G L,et al.2008.Temporal and spatial variation of the present crustal stress in northern part of North China[J].Acta Seismologica Sinica,30(6):570-580.(in Chinese)

    劉保金,曲國(guó)勝,孫銘心,等.2011.唐山地震區(qū)地殼結(jié)構(gòu)和構(gòu)造:深地震反射剖面結(jié)果[J].地震地質(zhì),33(4):901-912.Liu B J,Qu G S,Sun M X,et al.2011.Crustal structures and tectonics of Tangshan earthquake area:results from deep seismic reflection profiling[J].Seismology and Geology,33(4):901-912.(in Chinese)

    劉蒲雄,呂曉健.2011.唐山地震序列空間分布復(fù)雜性原因研究[J].地震,3(4):1-14.Liu P,Lyu X J.2011.On the causes of the spatial distribution complexity of the 1976 Tangshan earthquake sequence[J].Earthquake,31(4):1-14.(in Chinese)

    馬瑾,張渤濤,袁淑榮.1980.唐山地震與地震危險(xiǎn)區(qū)[J].地震地質(zhì),2(2):43-54.Ma J,Zhang B T,Yuan S R.1980.Tangshan earthquake and the earthquake risk area[J].Seismology and Geology,2(2):43-54.(in Chinese)

    孫麗娜,齊玉妍,陳婷,等.2021.震源機(jī)制和場(chǎng)地條件在河北地震影響場(chǎng)判定中的應(yīng)用研究[J].地震學(xué)報(bào),43(4):508-520.Sun L N,Qi Y Y,Chen T,et al.2021.Research on application of focal mechanism and site conditions in judgment of Hebei earthquake influence field[J].Acta Seismologica Sinica,43(4):508-520.(in Chinese)

    孫啟凱,池國(guó)民,徐東卓.2017.首都圈地區(qū)地殼垂直形變特征及剖面分析[J].大地測(cè)量與地球動(dòng)力學(xué),37(5):497-501.Sun Q K,Chi G M,Xu D Z.2017.Characteristics of vertical deformation and profile analysis for the capital circle[J].Journal of Geodesy and Geodynamics,37(5):497-501.(in Chinese)

    王輝,劉勉,李振,等.2013.華北地區(qū)大地震的長(zhǎng)余震序列[J].地震,33(3):1-12.Wang H,Liu M,Li Z,et al.2013.Long aftershock sequences of major earthquakes in North China[J].Earthquake,33(3):1-12.(in Chinese)

    王想,馮向東,宮猛,等.2012.唐山老震區(qū)地震活動(dòng)特征[J].地震地磁觀測(cè)與研究,33(S2):44-48.Wang X,F(xiàn)eng X D,Gong M,et al.2012.Study on seismicity features of Tangshan seismic region[J].Seismological and Geomagnetic Observation and Research,33(S2):44-48.(in Chinese)

    王曉山,馮向東,趙英萍.2020.京津冀地區(qū)地殼應(yīng)力場(chǎng)特征[J].地震研究,43(4):610-619.Wang X S,F(xiàn)eng X D,Zhao Y P.2020.Characteristics of crustal stress field in Beijing-Tianjin-Hebei region[J].Journal of Seismological Research,43(4):610-619.(in Chinese)

    聞學(xué)澤,馬勝利.2006.唐山大地震對(duì)相鄰斷裂段地震復(fù)發(fā)的影響[J].自然科學(xué)進(jìn)展,16(10):1346-1350.Wen X Z,Ma S L.2006.Effect of strong Tangshan earthquake on seismic recurrence in adjacent faults[J].Progress in Natural Science,16(10):1346-1350.(in Chinese)

    徐杰,周本剛,計(jì)鳳桔,等.2012.華北渤海灣盆地區(qū)大震發(fā)震構(gòu)造的基本特征[J].地震地質(zhì),34(4):618-636.Xu J,Zhou B G,Ji F J,et al.2012.Features of seismogenic structures of great earthquakes in the Bohai Bay Basin area,North China[J].Seismology and Geology,34(4):618-636.(in Chinese)

    許忠淮,閻明,趙仲和,等.1983.由多個(gè)小地震推斷的華北地區(qū)構(gòu)造應(yīng)力場(chǎng)的方向[J].地震學(xué)報(bào),5(3):268-279.Xu Z H,Yan M,Zhao Z H,et al.1983.Evaluation of the direction of tectonic stress in North China from recorded data of a large number of small earthquake[J].Acta Seismologica Sinica,5(3):268-279.(in Chinese)

    許忠淮.1985.用滑動(dòng)方向擬合法反演唐山余震區(qū)的平均應(yīng)力場(chǎng)[J].地震學(xué)報(bào),7(4):349-362.Xu Z H.1985.Mean stress field in Tangshan aftershock area obtained from focal mechanism data by fitting slip directions[J].Acta Seismologica Sinica,7(4):349-362.(in Chinese)

    楊理華,趙喜柱.1984.1976年唐山地震為三主震型地震[J].地震研究,7(6):639-648.Yang L H,Zhao X Z.1984.The Tangshan earthquake of 1976-an event of triple main shock pattern[J].Journal of Seismological Research,7(6):639-648.(in Chinese)

    張宏志,刁桂苓,陳祺福,等.2008.1976年唐山7.8級(jí)地震震區(qū)現(xiàn)今地震震源機(jī)制分析[J].地震研究,31(1):1-6.Zhang H Z,Diao G L,Chen Q F,et al.2008.Focal mechanism analysis of the recent earthquakes in Tangshan seismic region of M7.8 in 1976[J].Journal of Seismological Research,31(1):1-6.(in Chinese)

    張素欣,孫麗娜,李永慶,等.2020.2019年12月5日河北唐山M4.5地震的發(fā)震背景[J].華北地震科學(xué),38(S1):76-79.Zhang S X,Sun L N,Li Y Q,et al.2020.The seismogenic background of the M4.5 earthquake in Hebei Tangshan on December 5,2019[J].North China Earthquake Sciences,38(S1):76-79.(in Chinese)

    張素欣,王曉山,陳婷,等.2017.唐山老震區(qū)40年地震時(shí)空演化特征分析[J].華北地震科學(xué),35(1):32-37.Zhang S X,Wang X S,Chen T,et al.2017.Seismic space-time evolution characteristics in Tangshan old seismic region during the past 40 years[J].North China Earthquake Sciences,35(1):32-37.(in Chinese)

    張躍剛,王玉珍,尹寶軍,等.2013.2012年5月28日河北省唐山4.8級(jí)地震[J].中國(guó)地震,29(2):219-229.Zhang Y G,Wang Y Z,Yin B J,et al.2013.The Tangshan MS4.8 earthquake on May 28,2012 in Hebei Province[J].Earthquake Research in China,29(2):219-229.(in Chinese)

    鄭建常,李東梅,王鵬,等.2014.基于最鄰近事件距離算法的叢集地震識(shí)別方法及應(yīng)用[J].地震,34(4):100-109.Zhang J C,Li D M,Wang P,et al.2014.Method and application of clustering seismicity identification based on nearest-neighbor distance[J].Earthquake,34(4):100-109.(in Chinese)

    鄭建常,李冬梅,冀東普,等.2022.日本與華北地區(qū)強(qiáng)震活動(dòng)相關(guān)性再討論[J].華南地震,42(4):7-19.Zheng J C,Li D M,Ji D P,et al.2022.Rediscussion on the correlation of strong earthquake activities between Japan and North China[J].South China Journal of Seismology,42(4):7-19.(in Chinese)仲秋,史保平.2012.1976年唐山MS7.8地震余震序列持續(xù)時(shí)間及對(duì)地震危險(xiǎn)性分析的意義[J].地震學(xué)報(bào),34(4):494-508.Zhong Q,Shi B P.2012.Aftershock duration of the 1976 MS7.8 Tangshan earthquake and implication for seismic hazard estimation[J].Acta Seismologica Sinica,34(4):494-508.(in Chinese)

    朱琳,李騰飛,石富強(qiáng),等.2021.1976年唐山強(qiáng)震群震后庫(kù)侖應(yīng)力演化及其與2020年古冶5.1級(jí)地震的關(guān)系[J].地震研究,44(1):1-8.Zhu L,Li T F,Shi F Q,et al.2021.The coulomb stress evolution after the 1976 Tangshan strong earthquake swarm and its relationship with the 2020 Guye M5.2 earthquake[J].Journal of Seismological Research,44(1):1-8.(in Chinese)

    Butler R,Stewart G S,Kanamori H.1979.The July 27,1976 Tangshan,China,earthquake:a complex sequence of intraplate events[J].Bulletin of the Seismological Society of America,69(1):207-220.

    Li L,Chen Q F,Cheng X,et al.2007.Spatial clustering and repeating of seismic events observed along the 1976 Tangshan fault,north China[J].Geophysical Research Letters,34(23):L23309.

    Robinson R,Zhou S.2005.Stress interactions within the Tangshan,China,earthquake sequence of 1976[J].Bulletin of the Seismological Society of America,95(6):2501-2505.

    Shedlock K M,Baranowski J,Xiao W W,et al.1987.The Tangshan aftershock sequence[J].Journal of Geophysical Research,92(B3):2791-2803.

    Stein S,Liu M.2009.Long aftershock sequences within continents and implications for earthquake hazard assessment[J].Nature,462(5):87-89.

    Characteristics of the Recent Seismicity in the Tangshan Old Seismic Region

    ZHOU Yi1,2,CAO Jun2,MA Xudong2,SHAN Lianjun2

    (1.Hebei Hongshan National Observatory on Thick Sediments and Seismic Hazards,Xingtai 054000,Hebei,China)?(2.Hebei Earthquake Agency,Shijiazhuang 050021,Hebei,China)

    Abstract

    To analyze the seismic background in the Tangshan old seismic region and provide earthquake trend predictions,we describe the characteristics of the recent seismicity in this region.Results show that the recent seismicity in the Tangshan old seismic region is different from the early aftershock intensive area of the 1976 Tangshan MS7.8 earthquake.Considering the statistical results of regional background seismicity,the influence of the Tangshan earthquake is gradual weakening,and the recent seismicity in the Tangshan old seismic region can be considered as background seismicity.The characteristics of seismicity are different in the three sub-areas of the Tangshan old seismic region,the earthquakes mainly concentrate in the central area,followed by the eastern area,and the western area was the weakest,and the moderate earthquake following a law of alternate occurrence in the east,central and west sub-areas.This seismic migration law constitutes a reference for the determination of the locations of future moderate-to-strong earthquakes.At present,the accumulation of tectonic stress is accelerating in the old seismic region of Tangshan.Seismicity changing from quiet to active is possible,and moderate-to-strong earthquakes may occur at the northern segment of the Tangshan Fault and the intersection of the Luanxian-Laoting Fault and the Lulong Fault in the future.

    Keywords:?the 1976 Tangshan MS7.8 earthquake;characteristics of seismicity;background seismicity;seismic migration law

    收稿日期:2022-10-27.

    基金項(xiàng)目:河北省重點(diǎn)研發(fā)項(xiàng)目(20375404D);河北省地震科技星火計(jì)劃項(xiàng)目(DZ2021121700014、DZ2024112100001);河北省地震局創(chuàng)新團(tuán)隊(duì)項(xiàng)目(DZ202109080114).

    第一作者簡(jiǎn)介:周依(1991-),工程師,主要從事數(shù)字地震學(xué)和地震活動(dòng)性研究.E-mail:zhouyimonday@163.com.

    周依,曹筠,馬旭東,等.2024.唐山老震區(qū)現(xiàn)今地震活動(dòng)特征[J].地震研究,47(3):430-440,doi:10.20015/j.cnki.ISSN1000-0666.2024.0046.

    Zhou Y, Cao J, Ma X D,et al.2024.Characteristics of the recent seismicity in the Tangshan old seismic region[J].Journal of Seismological Research,47(3):430-440,doi:10.20015/j.cnki.ISSN1000-0666.2024.0046.

    猜你喜歡
    小震震區(qū)余震
    云南小江地區(qū)小震震源機(jī)制及構(gòu)造應(yīng)力場(chǎng)研究
    “超長(zhǎng)待機(jī)”的余震
    流浪衛(wèi)星
    2017年9月2日固原4.6級(jí)地震前后小震調(diào)制比分析
    跟蹤導(dǎo)練(五)3
    生死之間的靈魂救贖——《余震》和《云中記》的倫理問(wèn)題
    蘆山震區(qū)大田壩崩塌發(fā)育特征及其防治措施
    三次8級(jí)以上大地震的余震活動(dòng)特征分析*
    地震研究(2015年4期)2015-12-25 05:33:44
    小震調(diào)制比在寧夏及鄰區(qū)映震能力的研究
    強(qiáng)震區(qū)軟弱地基上承式連拱橋設(shè)計(jì)總結(jié)
    最近最新中文字幕免费大全7| av国产精品久久久久影院| 性色avwww在线观看| 在线观看www视频免费| 亚洲欧美色中文字幕在线| 日韩熟女老妇一区二区性免费视频| 日本免费在线观看一区| 91aial.com中文字幕在线观看| 又粗又硬又长又爽又黄的视频| 丰满少妇做爰视频| av线在线观看网站| 国产精品.久久久| 欧美日韩视频精品一区| 男女边吃奶边做爰视频| 韩国av在线不卡| 久久婷婷青草| 日本黄大片高清| 桃花免费在线播放| 久久精品国产综合久久久 | 日本色播在线视频| 久久 成人 亚洲| 在线观看www视频免费| 国产精品人妻久久久久久| 边亲边吃奶的免费视频| 黑人猛操日本美女一级片| 成年人午夜在线观看视频| 老女人水多毛片| 捣出白浆h1v1| 少妇人妻 视频| 亚洲国产精品一区三区| 一级黄片播放器| av视频免费观看在线观看| 欧美精品一区二区大全| 久久人人爽人人爽人人片va| 99久国产av精品国产电影| 午夜福利视频在线观看免费| 国产成人精品久久久久久| 色吧在线观看| 亚洲av免费高清在线观看| 丝袜喷水一区| 永久免费av网站大全| 亚洲精品色激情综合| 最近手机中文字幕大全| 成年动漫av网址| 成人毛片60女人毛片免费| 国产欧美另类精品又又久久亚洲欧美| 欧美日韩av久久| 亚洲三级黄色毛片| 这个男人来自地球电影免费观看 | av免费在线看不卡| 中文乱码字字幕精品一区二区三区| 免费黄频网站在线观看国产| 边亲边吃奶的免费视频| 黄片播放在线免费| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产精品成人久久小说| 美女中出高潮动态图| 大话2 男鬼变身卡| 久久ye,这里只有精品| 侵犯人妻中文字幕一二三四区| 成年动漫av网址| 母亲3免费完整高清在线观看 | 免费观看性生交大片5| 婷婷色av中文字幕| 少妇高潮的动态图| 各种免费的搞黄视频| 一本大道久久a久久精品| 国产亚洲av片在线观看秒播厂| 欧美日韩综合久久久久久| 少妇人妻久久综合中文| 亚洲精品国产av成人精品| 久久久a久久爽久久v久久| 亚洲欧美精品自产自拍| 97在线视频观看| 日韩 亚洲 欧美在线| av视频免费观看在线观看| 18在线观看网站| 欧美精品av麻豆av| 纯流量卡能插随身wifi吗| 女人久久www免费人成看片| 久久久久精品久久久久真实原创| 热99国产精品久久久久久7| 国产又爽黄色视频| av在线观看视频网站免费| 国产综合精华液| 丝袜在线中文字幕| 亚洲成人av在线免费| 我的女老师完整版在线观看| 美女国产高潮福利片在线看| 高清视频免费观看一区二区| www.色视频.com| 久久精品国产综合久久久 | 两个人看的免费小视频| 高清毛片免费看| 国产精品不卡视频一区二区| 国产成人精品久久久久久| 国产欧美日韩综合在线一区二区| 激情五月婷婷亚洲| 欧美日韩视频高清一区二区三区二| 精品国产乱码久久久久久小说| 免费大片18禁| 美女内射精品一级片tv| 精品一区二区三区视频在线| 欧美 亚洲 国产 日韩一| 黑人巨大精品欧美一区二区蜜桃 | 亚洲激情五月婷婷啪啪| 高清黄色对白视频在线免费看| 日日撸夜夜添| 亚洲色图综合在线观看| 国产精品一区www在线观看| 狠狠精品人妻久久久久久综合| 婷婷色av中文字幕| 男人添女人高潮全过程视频| 亚洲人成77777在线视频| 看非洲黑人一级黄片| 校园人妻丝袜中文字幕| 亚洲欧美日韩卡通动漫| 午夜福利影视在线免费观看| 精品卡一卡二卡四卡免费| 2021少妇久久久久久久久久久| 欧美日韩亚洲高清精品| 少妇的逼水好多| av.在线天堂| 亚洲天堂av无毛| 有码 亚洲区| av在线播放精品| 大香蕉97超碰在线| 亚洲美女搞黄在线观看| 少妇被粗大的猛进出69影院 | 久久久久精品性色| 飞空精品影院首页| h视频一区二区三区| 狠狠婷婷综合久久久久久88av| 免费不卡的大黄色大毛片视频在线观看| www日本在线高清视频| 一区二区日韩欧美中文字幕 | 午夜激情av网站| a级片在线免费高清观看视频| 日本与韩国留学比较| 亚洲欧美日韩卡通动漫| 日韩欧美一区视频在线观看| 又大又黄又爽视频免费| 国产精品99久久99久久久不卡 | 一区二区日韩欧美中文字幕 | 侵犯人妻中文字幕一二三四区| 91精品三级在线观看| 亚洲国产精品专区欧美| 亚洲五月色婷婷综合| 91精品三级在线观看| 美女国产视频在线观看| 麻豆乱淫一区二区| 国产av码专区亚洲av| 亚洲内射少妇av| 欧美丝袜亚洲另类| 欧美最新免费一区二区三区| 久久久久人妻精品一区果冻| 久久久a久久爽久久v久久| 亚洲一区二区三区欧美精品| 日日爽夜夜爽网站| 又黄又爽又刺激的免费视频.| 欧美精品人与动牲交sv欧美| 免费观看性生交大片5| 精品亚洲成国产av| 一区二区三区精品91| 久久韩国三级中文字幕| 国产精品国产av在线观看| 哪个播放器可以免费观看大片| 人妻人人澡人人爽人人| 亚洲久久久国产精品| 日本vs欧美在线观看视频| a 毛片基地| 久久鲁丝午夜福利片| 国产欧美日韩综合在线一区二区| 国产精品一区二区在线不卡| 亚洲精品日本国产第一区| 中文字幕制服av| 午夜免费观看性视频| 全区人妻精品视频| 美女国产视频在线观看| 午夜福利视频精品| 久久久久国产精品人妻一区二区| 乱码一卡2卡4卡精品| 久久国产精品大桥未久av| 搡女人真爽免费视频火全软件| 国产1区2区3区精品| 中文字幕亚洲精品专区| www.av在线官网国产| 波多野结衣一区麻豆| 免费大片18禁| 婷婷色av中文字幕| 亚洲,欧美,日韩| 日本与韩国留学比较| 久久影院123| 麻豆乱淫一区二区| 午夜免费观看性视频| 中文字幕人妻熟女乱码| 久久久久精品性色| 女人久久www免费人成看片| 一边亲一边摸免费视频| 精品人妻一区二区三区麻豆| 日本午夜av视频| 又黄又粗又硬又大视频| 国产一区二区激情短视频 | 亚洲精品日本国产第一区| 免费看光身美女| 亚洲av福利一区| 亚洲国产欧美日韩在线播放| 亚洲精品自拍成人| 99久久综合免费| 午夜视频国产福利| xxx大片免费视频| 如何舔出高潮| 欧美亚洲 丝袜 人妻 在线| 久久精品国产综合久久久 | 好男人视频免费观看在线| 久久精品aⅴ一区二区三区四区 | 日韩三级伦理在线观看| 精品国产国语对白av| 男人舔女人的私密视频| 免费大片黄手机在线观看| 五月伊人婷婷丁香| 中文字幕av电影在线播放| 精品99又大又爽又粗少妇毛片| 精品酒店卫生间| 日本wwww免费看| 丰满少妇做爰视频| 天天影视国产精品| 精品熟女少妇av免费看| 一级片'在线观看视频| 国产 精品1| 亚洲丝袜综合中文字幕| 欧美精品人与动牲交sv欧美| 在线观看国产h片| 亚洲精品456在线播放app| 欧美人与善性xxx| 精品一区二区三区四区五区乱码 | 国产福利在线免费观看视频| 91精品国产国语对白视频| 午夜福利乱码中文字幕| 黑人高潮一二区| 999精品在线视频| 国产成人精品无人区| 国产免费视频播放在线视频| 久久久久久久久久成人| 日韩av免费高清视频| 18禁观看日本| 汤姆久久久久久久影院中文字幕| 热99国产精品久久久久久7| 岛国毛片在线播放| 亚洲婷婷狠狠爱综合网| 国产日韩一区二区三区精品不卡| videos熟女内射| 日本av免费视频播放| 曰老女人黄片| 大香蕉久久网| 日韩三级伦理在线观看| 成人18禁高潮啪啪吃奶动态图| 欧美亚洲日本最大视频资源| 男女下面插进去视频免费观看 | 国产熟女欧美一区二区| 看免费av毛片| 亚洲,欧美精品.| 黑人高潮一二区| 成人二区视频| 只有这里有精品99| 天堂中文最新版在线下载| 欧美精品高潮呻吟av久久| 天天操日日干夜夜撸| 成人影院久久| 久热久热在线精品观看| 国产白丝娇喘喷水9色精品| 久久这里只有精品19| 2021少妇久久久久久久久久久| 日韩欧美一区视频在线观看| 黄片无遮挡物在线观看| 久久精品久久精品一区二区三区| 亚洲伊人色综图| 亚洲欧美精品自产自拍| 久久午夜综合久久蜜桃| 十八禁高潮呻吟视频| 纵有疾风起免费观看全集完整版| 91久久精品国产一区二区三区| 国产精品久久久久久久久免| 免费观看av网站的网址| 狂野欧美激情性xxxx在线观看| 人妻系列 视频| 亚洲精品一二三| 最新的欧美精品一区二区| 日韩av免费高清视频| 国产精品一区www在线观看| 日韩在线高清观看一区二区三区| 日本色播在线视频| 精品久久久精品久久久| 久久这里有精品视频免费| 在线观看免费日韩欧美大片| 91在线精品国自产拍蜜月| 日韩制服丝袜自拍偷拍| 欧美xxⅹ黑人| 丰满乱子伦码专区| 国产精品蜜桃在线观看| 99国产综合亚洲精品| 大片电影免费在线观看免费| 欧美3d第一页| 久久精品人人爽人人爽视色| 久久人人爽人人爽人人片va| 免费看av在线观看网站| 国产av精品麻豆| 亚洲精品av麻豆狂野| 尾随美女入室| 狂野欧美激情性xxxx在线观看| 精品视频人人做人人爽| 国产精品免费大片| 精品一品国产午夜福利视频| 高清av免费在线| 国产亚洲精品第一综合不卡 | 国产精品一二三区在线看| 精品人妻偷拍中文字幕| 曰老女人黄片| 高清黄色对白视频在线免费看| 美国免费a级毛片| 色94色欧美一区二区| 日本爱情动作片www.在线观看| 熟女av电影| 午夜久久久在线观看| 久久国产精品男人的天堂亚洲 | 内地一区二区视频在线| 亚洲av欧美aⅴ国产| 国精品久久久久久国模美| 久久国产精品男人的天堂亚洲 | 国产精品久久久久久久久免| 国产高清不卡午夜福利| 91久久精品国产一区二区三区| 成年av动漫网址| 伊人亚洲综合成人网| 亚洲 欧美一区二区三区| 欧美精品高潮呻吟av久久| 国产黄频视频在线观看| 9191精品国产免费久久| 成人午夜精彩视频在线观看| 国产一区二区在线观看日韩| 五月天丁香电影| 免费看光身美女| 国产男人的电影天堂91| 啦啦啦啦在线视频资源| 免费观看在线日韩| 综合色丁香网| 亚洲国产毛片av蜜桃av| 一区二区三区精品91| 美国免费a级毛片| 99re6热这里在线精品视频| 大码成人一级视频| 久热久热在线精品观看| 国产女主播在线喷水免费视频网站| 看免费av毛片| 成人综合一区亚洲| 哪个播放器可以免费观看大片| 国产精品久久久久成人av| 宅男免费午夜| 色吧在线观看| 国产日韩欧美视频二区| 国产精品一区二区在线观看99| 午夜久久久在线观看| 在线亚洲精品国产二区图片欧美| 亚洲精品aⅴ在线观看| 亚洲欧美日韩卡通动漫| 精品一品国产午夜福利视频| 国产国语露脸激情在线看| 51国产日韩欧美| 国产精品一区二区在线观看99| av天堂久久9| 久久久国产欧美日韩av| 街头女战士在线观看网站| 国产精品熟女久久久久浪| 中文字幕人妻丝袜制服| 在现免费观看毛片| 亚洲精品aⅴ在线观看| 9热在线视频观看99| 亚洲一码二码三码区别大吗| 涩涩av久久男人的天堂| 在线观看三级黄色| 考比视频在线观看| 久久人人爽av亚洲精品天堂| 在线天堂中文资源库| 亚洲精品av麻豆狂野| 国产精品.久久久| 大片免费播放器 马上看| 国产精品一区二区在线观看99| 好男人视频免费观看在线| 久久久久久久久久人人人人人人| av电影中文网址| 肉色欧美久久久久久久蜜桃| 亚洲四区av| 国产女主播在线喷水免费视频网站| 亚洲国产精品成人久久小说| 2022亚洲国产成人精品| 亚洲精品中文字幕在线视频| 人人妻人人澡人人看| 久久精品国产鲁丝片午夜精品| 精品久久久久久电影网| 亚洲欧美色中文字幕在线| 高清av免费在线| 成人午夜精彩视频在线观看| 青青草视频在线视频观看| 久久精品国产亚洲av涩爱| 高清不卡的av网站| 国产有黄有色有爽视频| 欧美日韩av久久| 亚洲av国产av综合av卡| 妹子高潮喷水视频| 最黄视频免费看| 成年av动漫网址| 亚洲av男天堂| 成人综合一区亚洲| 久久午夜综合久久蜜桃| 九草在线视频观看| 性色avwww在线观看| 国产成人精品久久久久久| 9热在线视频观看99| 亚洲精品日韩在线中文字幕| 人妻系列 视频| 成人黄色视频免费在线看| 国产欧美日韩一区二区三区在线| 美女国产视频在线观看| 久久精品熟女亚洲av麻豆精品| 男女高潮啪啪啪动态图| 亚洲久久久国产精品| 国产成人免费无遮挡视频| 欧美日韩综合久久久久久| 国产精品熟女久久久久浪| 日韩人妻精品一区2区三区| 人妻系列 视频| 国产熟女午夜一区二区三区| 99久久人妻综合| 成人国产av品久久久| 下体分泌物呈黄色| 黄片无遮挡物在线观看| 亚洲伊人色综图| 777米奇影视久久| 少妇的逼水好多| 色视频在线一区二区三区| 国产精品.久久久| 亚洲精品国产av成人精品| 中文天堂在线官网| 婷婷成人精品国产| 精品一区二区免费观看| 久久久国产欧美日韩av| 亚洲av.av天堂| 九九爱精品视频在线观看| 国产极品天堂在线| 又黄又粗又硬又大视频| 日韩熟女老妇一区二区性免费视频| 黄片播放在线免费| 少妇被粗大的猛进出69影院 | 久久久久精品性色| 日本-黄色视频高清免费观看| 久久午夜福利片| 久久鲁丝午夜福利片| 黑人欧美特级aaaaaa片| 久久久精品免费免费高清| 国产男人的电影天堂91| 男女午夜视频在线观看 | 高清毛片免费看| 欧美日韩精品成人综合77777| 欧美3d第一页| 亚洲国产欧美日韩在线播放| 一区二区av电影网| 免费人成在线观看视频色| 成人18禁高潮啪啪吃奶动态图| 大香蕉久久网| 女性被躁到高潮视频| 汤姆久久久久久久影院中文字幕| 熟妇人妻不卡中文字幕| 亚洲国产欧美在线一区| 中国美白少妇内射xxxbb| 亚洲,一卡二卡三卡| 国产av精品麻豆| 熟女人妻精品中文字幕| 99香蕉大伊视频| 国产色爽女视频免费观看| 欧美少妇被猛烈插入视频| 亚洲国产日韩一区二区| 热99国产精品久久久久久7| 欧美亚洲 丝袜 人妻 在线| 熟妇人妻不卡中文字幕| 日韩av免费高清视频| 欧美激情国产日韩精品一区| 午夜老司机福利剧场| 国产永久视频网站| 久热久热在线精品观看| 日韩中文字幕视频在线看片| 国产成人一区二区在线| 亚洲国产最新在线播放| 色网站视频免费| 好男人视频免费观看在线| 成人免费观看视频高清| xxxhd国产人妻xxx| 欧美人与善性xxx| 欧美另类一区| 日韩精品免费视频一区二区三区 | 国产精品久久久久久精品电影小说| 国产有黄有色有爽视频| 好男人视频免费观看在线| 日韩熟女老妇一区二区性免费视频| 高清av免费在线| 精品一区二区三卡| 国产精品国产三级国产av玫瑰| 自线自在国产av| 99香蕉大伊视频| 欧美日韩亚洲高清精品| 欧美日韩视频精品一区| 啦啦啦啦在线视频资源| 超色免费av| a级毛色黄片| 菩萨蛮人人尽说江南好唐韦庄| 久久婷婷青草| 91精品国产国语对白视频| 边亲边吃奶的免费视频| 午夜精品国产一区二区电影| 日韩av免费高清视频| 成年人免费黄色播放视频| 色婷婷久久久亚洲欧美| 日本午夜av视频| 国产又爽黄色视频| 啦啦啦中文免费视频观看日本| 欧美xxⅹ黑人| 国产精品久久久久久久久免| 亚洲综合色网址| 欧美激情国产日韩精品一区| av不卡在线播放| 大香蕉久久成人网| 国产精品麻豆人妻色哟哟久久| 免费大片18禁| 精品熟女少妇av免费看| 亚洲成人一二三区av| 国产欧美另类精品又又久久亚洲欧美| 亚洲久久久国产精品| 国内精品宾馆在线| 免费观看a级毛片全部| 国产免费又黄又爽又色| 亚洲色图综合在线观看| 伊人亚洲综合成人网| 天天躁夜夜躁狠狠躁躁| 国国产精品蜜臀av免费| 亚洲成色77777| 久久久久国产精品人妻一区二区| 日韩制服骚丝袜av| 亚洲欧美成人精品一区二区| 午夜福利影视在线免费观看| 日本av免费视频播放| 日本欧美视频一区| 免费观看无遮挡的男女| 永久免费av网站大全| 国产爽快片一区二区三区| 夫妻午夜视频| 亚洲人成77777在线视频| 三级国产精品片| 99久久中文字幕三级久久日本| 深夜精品福利| 在现免费观看毛片| 久久精品久久久久久噜噜老黄| 精品人妻在线不人妻| 久久午夜福利片| 一区二区三区精品91| 亚洲欧美色中文字幕在线| 日本爱情动作片www.在线观看| 午夜福利网站1000一区二区三区| 纯流量卡能插随身wifi吗| 免费av中文字幕在线| 99精国产麻豆久久婷婷| 国产无遮挡羞羞视频在线观看| 黑人高潮一二区| 日韩,欧美,国产一区二区三区| 国产老妇伦熟女老妇高清| 久久久久精品久久久久真实原创| 午夜免费男女啪啪视频观看| av电影中文网址| 免费观看a级毛片全部| 午夜免费观看性视频| 午夜91福利影院| 成人手机av| 免费人成在线观看视频色| 男女国产视频网站| 香蕉丝袜av| 女人被躁到高潮嗷嗷叫费观| 久久久精品94久久精品| 色网站视频免费| 亚洲精品中文字幕在线视频| 免费久久久久久久精品成人欧美视频 | 99re6热这里在线精品视频| 王馨瑶露胸无遮挡在线观看| 一级毛片黄色毛片免费观看视频| 成年人午夜在线观看视频| 国产极品天堂在线| 午夜精品国产一区二区电影| 亚洲一级一片aⅴ在线观看| 日韩欧美一区视频在线观看| 国产成人av激情在线播放| 人人妻人人澡人人爽人人夜夜| 99久久精品国产国产毛片| 视频在线观看一区二区三区| 欧美成人精品欧美一级黄| 一级毛片 在线播放| 日韩欧美一区视频在线观看| 午夜免费鲁丝| 插逼视频在线观看| 最近的中文字幕免费完整| 婷婷色综合大香蕉| 少妇精品久久久久久久| 免费av中文字幕在线| 黄色配什么色好看| 热99国产精品久久久久久7| 久久99精品国语久久久| 黑人巨大精品欧美一区二区蜜桃 | 两个人免费观看高清视频| av电影中文网址|