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

    從同震和震后形變分析1668年M8.5級郯城地震對周邊地震活動性的影響

    2016-12-05 02:28:57孫玉軍吳中海
    地質力學學報 2016年3期
    關鍵詞:郯城同震庫侖

    徐 昊,孫玉軍,吳中海

    (中國地質科學院地質力學研究所,北京 100081)

    ?

    從同震和震后形變分析1668年M8.5級郯城地震對周邊地震活動性的影響

    徐 昊,孫玉軍,吳中海

    (中國地質科學院地質力學研究所,北京 100081)

    以山東郯城1668年大地震為例,以前人地表地質調查結果為約束,利用彈性位錯理論初步獲取了該地震的同震破裂模型;在此基礎上,基于粘彈性分層模型分析了該地震的同震和震后形變,同時以主震斷層為接收斷層計算了庫侖應力分布,進一步討論了地幔不同粘滯性系數(shù)對地表形變和庫侖應力變化的影響。計算結果顯示,該地震是一個右旋走滑為主兼有一定逆沖性質的地震,其同震位移巨大,能量釋放較徹底;同震破裂造成震中郯城縣西北、東北和南部部分斷層庫侖應力增加,而震后形變使得這些斷層庫侖應力進一步增加,在單縣、宿遷和日照等地,地震后350 a庫侖應力變化量達到+1bar—+1MPa量級;地幔粘滯性系數(shù)不同,形變量和庫侖應力變化達到穩(wěn)定的時間不同,但最終趨于穩(wěn)定的數(shù)值基本一致。

    郯城大地震;同震和震后形變;庫侖應力變化;粘滯性

    0 引言

    郯廬斷裂帶是中國東部一條重要的活動斷裂帶,歷史上沿該斷裂帶發(fā)生過多次大地震。1668年7月25日晚,在山東南部發(fā)生了一次曠古未有的特大地震,根據歷史資料得到該地震的震級約為Ms8.5級,極震區(qū)位于山東省郯城、臨沭、臨沂交界處,震中位置(34.8°N,118.5°E),極震區(qū)烈度達Ⅻ級[1]。由于極震區(qū)大部分位于郯城縣境內,故稱為郯城地震。這是我國大陸東部板塊內部最強烈的地震,造成巨大財產損失和人員傷亡。

    關于1668郯城地震,我國學者已經對其進行了相關研究,分析了這次地震發(fā)震斷層的構造特征和破裂機制[2~3],采用不同手段得到了此次地震震源參數(shù)[4~6],估算了郯城8.5級地震的復發(fā)周期[7~8]。但關于此次地震造成的應力場變化研究很少。近年來,越來越多的研究表明,地震發(fā)生會隨之釋放斷層上積聚的應力,但應力并不會憑空消失,部分應力轉移到其他地區(qū),導致應力積聚,進而誘發(fā)后續(xù)地震,此即地震的應力觸發(fā)理論[9]。根據應力觸發(fā)理論,庫侖應力增強相當于斷層額外負荷的增加,容易誘發(fā)后續(xù)地震;反之,庫侖應力減弱的區(qū)域——應力影區(qū),斷層負荷部分卸載,地震發(fā)生的概率降低。1668年郯城8.5級特大地震的同震及震后應力特征如何?其造成的應力場變化對周圍地區(qū)的地震活動性有何影響?本文圍繞這些問題,利用彈性位錯理論和分層巖石圈模型,計算了郯城1668年地震同震和震后形變以及庫侖應力的變化,探討本次地震對郯廬斷裂帶及附近區(qū)域地震活動性的影響,以期為進一步研究該區(qū)域的地震活動性提供參考。

    1 區(qū)域構造背景

    地震的發(fā)生與活動斷層聯(lián)系密切,斷層活動是地震產生的原因。郯廬斷裂帶是中國大陸東部一條北北東向的深大斷裂,也是中國東部有巨大影響的強震活動帶,史料記載曾發(fā)生過公元前70年安丘7級地震、1597年渤海7級地震、1668年郯城8.5級地震、1888年渤海7.5級地震、1969年渤海7.4級地震、1975年海城7.3級地震等[10]。其中1668年郯城8.5級地震的地震斷層為一高角度右旋走滑逆斷層,呈北北東向展布(見圖1),斷層面向東傾[6]。

    圖1 研究區(qū)構造背景圖(橢圓內加粗斷裂為此次地震主震斷層)Fig.1 Tectonic background map of Tanlu fault and its surrounding regions

    2 計算方法

    2.1 分層模型

    用于地震研究的彈性位錯模型主要可分為3類,第一類主要基于均勻彈性半空間介質模型[11],第二類主要基于分層介質模型[12],第三類基于成層球對稱地球模型的球形位錯模型[13]。本文基于地殼的分層模型采用格林函數(shù)解求解這次地震引起的同震位移以及引起的應力場變化。震后很長時間,有許多影響震后應力變化的因素,如震后余滑、孔隙效應和粘性松弛等,其中震后粘性松弛效應對震后形變和應力起到非常重要的作用[14]。Wang[15]提出了利用正交歸一法計算地震應力場Green函數(shù)方法,并在此基礎上建立了粘彈松弛分層模型下的地震同震及震后形變模型,編寫了數(shù)值計算程序——PSGRN/PSCMP軟件,用該程序基于分層、重力作用下的彈性模型來研究由于斷層錯動引起的同震和震后形變。假設郯城及周邊地區(qū)震后彈性松弛的過程主要發(fā)生在30 km上地殼以下。震后粘彈松弛應力變化是指一定深度下,受高溫高壓環(huán)境的影響,巖石層的力學性質逐漸由脆性向粘塑性轉化。由于地震發(fā)生時破裂的速度很快,在短時間內粘彈性的中下地殼以及上地幔力學性質表現(xiàn)為彈性體,然后震后隨時間推移,中下地殼和上地幔所積累的應力和應變逐漸釋放,并向上傳輸?shù)綇椥缘纳系貧ぶ?,導致上地殼中應力狀態(tài)發(fā)生改變[16]。對于地殼和上地幔來說,粘滯性系數(shù)不同,對震后形變和庫侖應力變化影響不同。粘滯性系數(shù)越低,松弛時間越短,震后短期內形變和庫侖應力增加的都比較快,很快能達到穩(wěn)定狀態(tài);相反,粘滯性系數(shù)越高,震后短期形變和庫侖應力增加的都比較慢,但最終的穩(wěn)定值相對較高[17]。模型假定地殼厚度為30 km,視為完全彈性體,粘滯性系數(shù)取無窮大;研究區(qū)缺乏長期的形變觀測資料和相關的震后形變研究成果,模型復雜化需要更多的參數(shù)來約束,因此采用簡單的Maxwell體分析地幔粘彈性變形。由于未獲取研究區(qū)地幔準確粘滯性系數(shù),本文模型中計算所需的粘滯性系數(shù)取1.0×1018Pa·s,最后會分析不同地幔粘滯性系數(shù)條件下,計算得到的形變量和庫侖應力的差異,并結合李勇紅等[5]給出的郯城8.5級地震震源區(qū)附近地殼上地幔速度結構,給出最終的分層模型及參數(shù)(見表1)。

    表1 分層模型及參數(shù)

    2.2 庫侖應力變化的計算

    隨著斷層上應力不斷累積,最終剪應力能夠克服斷裂上的靜摩擦力,鎖閉的斷層會發(fā)生破裂,斷層的破裂符合庫倫破裂準則[18]。利用彈性位錯理論,可以計算地震同震破裂所引起的應力場變化[19]。所研究的接收斷層庫侖應力變化為:

    (1)

    式中:Δσf為計算斷層面上庫侖應力的變化,Pa;Δτs是斷層面上剪應力的變化,Pa;μ′為有效摩擦系數(shù);Δσn是計算斷層面上正應力的變化,Pa。關于μ′的取值,不同學者有所差別。King等[20]認為有效摩擦系數(shù)對庫侖應力變化的結果將產生中等大小的影響,本文在計算過程中取中間值0.4。

    通過主震斷層模型計算靜態(tài)庫侖應力不但能夠有效反映余震的發(fā)展態(tài)勢,同時也能夠獲取對周圍活動斷層的影響[17]。由于此次地震發(fā)生時間距今已有350 a,本文僅選取一條主震斷層計算庫侖應力。采用之前學者研究所獲得的斷層相關參數(shù),走向21.57°、傾角89.51°、滑移角148.9°;起始破裂坐標為34.5°N、118.35°E,斷層面起始深度及其實際寬度取0 km和50 km[21],斷層破裂長度為150 km[6]。

    3 計算結果與分析

    3.1 同震和震后形變

    1668年郯城8.5級主震斷層為一高角度右旋走滑兼有逆沖性質的斷層,呈北北東向展布,斷層面東傾,震后觀察到水平最大形變量和垂直最大形變量分別為10 m和3 m[6]。盡管計算過程是在半無限空間假設模型下得到的,未考慮相關地形的影響,但通過圖像反映出的形變,依然可以判斷主震斷層為一右旋走滑逆斷層。計算得到同震時刻地表水平變形量和350 a后地表水平變形量(見圖2a,2c),同震時刻地表垂直變形量和350 a后地表垂直變形量(見圖2b,2d),其中350 a以后水平方向和垂直方向變形量分別為8.33 m和2.73 m,與王華林等[6]觀測數(shù)據存在一定誤差,由于在模型構建中缺少斷層面上精準的滑移數(shù)據,因此認為斷層模型依然能夠在一定程度上反映當時的地震破裂變形。

    a—同震時刻地表水平方向形變量;b—同震時刻地表垂直方向形變量(向上為正);c—350 a后地表水平方向形變量;d—350 a后地表垂直方向形變量(向上為正)圖2 同震及震后0 km處地表形變量Fig.2 The coseismic and postseismic deformation of the surface

    對比地表處水平方向的形變可以發(fā)現(xiàn),同震時刻和350 a以后斷層右旋的性質沒有改變,但350 a后水平形變量稍大于同震時刻水平形變量,距離主震斷層越遠增量越大。對比地表處垂直方向的形變可以發(fā)現(xiàn),同震時刻和350 a以后斷層逆沖的性質也未改變,但局部影區(qū)形變方向發(fā)生了倒轉,以右下角區(qū)域為例,同震時刻位移方向垂直向上,350 a后則向下變形。推測以上兩種變化原因與地殼巖石粘滯性有關:同震時刻斷層兩盤發(fā)生水平方向錯動,從而使距離破裂帶較近的區(qū)域瞬間獲得較大的形變量,而震后距破裂帶較遠的區(qū)域水平形變在巖石粘滯性作用下緩慢增加;已知主震斷層滑移角為149.5°,同震時刻斷層上下盤沿著走向錯動,必然會引起所在盤沿滑移方向一側的抬升,因此同震時刻右下和左上形變均為向上方向;隨著時間的推移,由于巖石粘滯性影響,變形逐漸恢復甚至引起反向形變,從而使圖中原先向上方向的變形轉變?yōu)橄蛳伦冃?。關于這種變化的量化分析,會在后來的部分進一步論述。

    3.2 庫侖應力變化

    每一次大地震發(fā)生都會引起周圍區(qū)域的應力增加或降低。一般來說,如果一個斷層的庫侖應力增加,則該斷層發(fā)生地震的危險性增強,反之則降低。同震破裂造成震中郯城縣西北、東北和南部部分地區(qū)庫侖應力增加,而震后形變使得這些地區(qū)庫侖應力進一步增加,在單縣、宿遷和日照等地,地震350 a后庫侖應力變化量達到+1 bar—+1 MPa量級(見圖3)。觀察研究區(qū)內震后350 a間Ms1—Ms7地震事件發(fā)現(xiàn),地震主要發(fā)生在這些庫侖應力增加的區(qū)域內。雖然地震的發(fā)生與構造運動直接相關,但這樣的特大型地震發(fā)生后對周圍斷層的應力積累影響十分明顯,而且計算結果顯示震后350 a該影響依然存在,只是目前已經趨于穩(wěn)定。

    同震庫侖應力增加最明顯區(qū)域主要位于斷層破裂帶的兩端及同震破裂不完全的區(qū)域。斷層破裂帶應力釋放明顯,而破裂帶的兩端則為應力集中的區(qū)域。本次地震的震源深度為20 km左右,因此,20 km深度同震應力釋放較地表同震應力釋放更為明顯(見圖3a,3b),沿著破裂帶主要為應力影區(qū);而在地表有部分滑動位移小,破裂不完全,應力增加的區(qū)域(見圖3a,3c)。對比同震時刻和350 a的同震庫侖應力分布,可以看出震后350 a在地表和20 km處應力增加的區(qū)域在面積和數(shù)值上都明顯增大。將本次地震后該區(qū)域的所有地震投影在平面后發(fā)現(xiàn),該區(qū)域的大部分地震事件也主要位于應力增加的區(qū)域。

    為了更清楚地觀察各應力影區(qū)的形變和庫侖應力變化趨勢,筆者選取南京、宿遷、單縣和日照4座城市計算并繪制其在深度20 km處350 a間位移和庫侖應力變化曲線。南京和日照在地震發(fā)生后庫侖應力都減少,不同的是南京在地震發(fā)生約20 a后庫侖應力逐漸增加,至今增加到-25 kPa左右(見圖4a);而日照在地震反生后庫侖應力稍微減少,隨后一直增加,到后來變化值穩(wěn)定在+1000 kPa左右(見圖4b)。宿遷和單縣則在地震發(fā)生后庫侖應力一直增加,最終宿遷庫侖應力增加穩(wěn)定在+1300 kPa左右(見圖4c),而單縣僅增加了+150 kPa左右(見圖4d)。庫侖應力增量表明宿遷和日照等地與郯廬斷裂帶主震斷層相同性質的斷層地震危險性增強。

    考慮同震和震后粘彈松弛效應對應力場變化的影響時,地幔的粘滯性系數(shù)不同取值可能會改變研究區(qū)的應力場分布。由于U(x)反映南北方向形變量,斷層走向為21.57°且為右旋斷層,U(x)的值能清楚反映斷層的運動性質。本文選取不同地幔粘滯性系數(shù)計算了郯城地震發(fā)生后350 a內,南京市深度20 km處的U(x)和庫侖應力變化值(見圖5),結果表明,震后形變量和庫侖應力變化量初始值相同,受粘滯性系數(shù)影響較小。隨著時間的推移,不管粘滯性大小,最終形變量穩(wěn)定值都近似相同,均為-1 m,但粘滯性系數(shù)越小,U(x)越快達到穩(wěn)定值,當粘滯性系數(shù)為1×1017Pa·s時,10 a左右就可達到穩(wěn)定值,當粘滯性系數(shù)為7×1018Pa·s時,則需要300 a左右才能達到穩(wěn)定值(見圖5a)。對于庫侖應力變化來說,震后初值變化也都相同,均為-17.5 kPa,曲線變化模式類似,都是庫侖應力先減少后增加最后趨于穩(wěn)定;當粘滯性系數(shù)為1×1017Pa·s時,50 a左右即可達到穩(wěn)定值,約為-24 kPa;當粘滯性系數(shù)為1×1018Pa·s時,150 a左右達到穩(wěn)定,約為-22.5 kPa;當粘滯性系數(shù)為7×1018Pa·s時,350 a時仍未達到穩(wěn)定狀態(tài),但可以估計其達到穩(wěn)定狀態(tài)的值基本相同(見圖5b)。

    a—地表處同震庫侖應力變化分布;b—20 km深處同震庫侖應力變化分布;c—地表處震后350 a庫侖應力變化分布;d—20 km深處震后350 a庫侖應力變化分布圖3 同震庫侖應力變化分布圖(灰色圓圈為1668年郯城地震后至今研究區(qū)內發(fā)生的M1—M7地震事件)Fig.3 The distribution map of coseismic Coulomb stress change

    a—南京市;b—日照市;c—宿遷市;d—單縣市△CFS—庫侖應力變化量;U(x)—南北向位移(北方向為正);U(y)—東西向位移(東方向為正);U(z)—垂向位移(向上為正)圖4 研究區(qū)內不同城市震后350 a內形變量和庫侖應力變化Fig.4 Deformation and Coulomb stress change in different cities in study regions 350 years after the earthquake

    圖5 不同地幔粘滯性系數(shù)對應的震后形變和庫侖應力變化(南京市)Fig.5 The postseismic deformation and Coulomb stress change of different mantle viscosities(Nanjing)

    4 結論

    本文基于粘彈松弛分層模型下的地震同震及震后形變模型,采用數(shù)值計算程序——PSGRN/PSCMP軟件,利用前人地表觀測和反演的地震數(shù)據計算了主震斷層同震以及震后的地表形變和庫侖應力變化量。

    通過對比同震和350 a后震中郯城縣東南部分位移形變量,從垂直向上的位移量轉變?yōu)橄蛳挛灰屏?,可知震后最終位移方向在受地殼粘彈性效應的影響下會發(fā)生倒轉。

    以主震斷層為接收斷層計算了同震和震后庫侖應力。同震庫侖應力增加的區(qū)域主要位于斷層破裂帶的兩端及同震破裂不完全的區(qū)域,具體分布于震中郯城縣西北、東北和東南的部分區(qū)域,震后形變使得這些地區(qū)庫侖應力進一步增加,在單縣、宿遷和日照等地,地震350 a后庫侖應力變化量達到+1 bar—+1 MPa量級。其他地區(qū)庫侖應力降低,南京地震350 a后庫侖應力變化量為-25 kPa。震后該區(qū)域發(fā)生的地震主要分布于庫侖應力增加的區(qū)域。

    通過對比地幔不同粘滯性系數(shù)對震后形變和庫侖應力的影響,發(fā)現(xiàn)對于同一區(qū)域,粘滯性系數(shù)越小,震后形變量和庫侖應力變化達到穩(wěn)定的時間越短,反之則時間越長;不同的地幔粘滯性系數(shù)對最終穩(wěn)定的形變量和庫侖應力影響不大。

    在討論此次地震與區(qū)域內周邊地震活動性關系時,主要依據庫侖應力變化量,雖然地震的發(fā)生與構造運動直接相關,但這樣的特大型地震發(fā)生后對周圍斷層的應力積累影響十分明顯,而且計算結果顯示震后350 a后該影響依然存在,但目前已經趨于穩(wěn)定,實際地震的發(fā)生還要綜合考慮構造活動等其他因素。此外由于1668—1970年之間區(qū)域地震事件缺少深度數(shù)據,所以在考慮后續(xù)地震在0 km和20 km處分布與庫侖應力變化的關系時誤差較大。

    [1] 馬玉香,鐘普裕.1668年山東郯城8.5級地震綜述[J].國際地震動態(tài),2009,(2):9~18.

    MA Yu-xiang,ZHONG Pu-yu.Summarization of the Tancheng Earthquake with 8.5 in 1668[J]. Recent developments in World Seismology, 2009, (2): 9~18.

    [2] 高維明,鄭朗蓀,李家靈,等.1668年郯城8.5級地震的發(fā)震構造[J].中國地震,1988,4(3):9~15.

    GAO Wei-ming,ZHENG Lang-sun,LI Jia-ling.The triggering structure of theMs8.5 earthquake of Tan Cheng in 1668[J].Earthquake Research in China,1988, 4(3): 9~15.

    [3] 鞠林雪.1668年郯城8.5級地震發(fā)震斷層構造特征與應力場分析[D].合肥:合肥工業(yè)大學,2012.

    JU Lin-xue.The research on the structure of triggring faults and stress field in 1668: The Tancheng earthquake[D]. Hefei: Hefei University of Technology, 2012.

    [4] 李清河,張元生,鮑海英,等.據地殼速度結構推測1668年山東郯城8.5級大地震震中[J].中國地震,2014,30(1):30~42.

    LI Qing-he, ZHANG Yuan-sheng, BAO Hai-ying, et al. Crust veolcity structure of the epicenter ofM8.5 earthquake at Tancheng, Shandong, China, in 1668[J]. Earthquake Research in China, 2014, 30 (1): 30~42.

    [5] 李永紅,胡新亮,許萍,等.以雙差定位方法對郯城8.5級地震震中附近現(xiàn)代小震重新定位[J].西北地震學報,2011,33(1):71~75.

    LI Yong-hong, HU Xin-liang, XU Ping, et al. Relocation of modern small earthquakes near the epicenter of Tancheng 8.5 earthquake using the Double Difference Earthquaake Location Algorithm[J]. Northwestern Seismological Journal, 2011, 33(1): 71~75.

    [6] 王華林,耿杰.1668年郯城8.5級地震震源參數(shù)及其討論[J].地震學刊,1996,(4):27~33.

    WANG Hua-lin, GENG Jie. Discussion about focus parameters of Tancheng earthquake ofMs8.5 in 1668[J]. Journal of Seismology, 1996, (4): 27~33.

    [7] 王華林.1668年郯城8.5級地震斷裂的全新世滑動速率,古地震和強震復發(fā)周期[J].西北地震學報,1995,17(4):1~12,22.

    WANG Hua-lin. Holocene displacement rate,paleoearthquakes and recurrence intervals of strong earthquakes along the 1668 Tancheng Earthquake (Ms=8.5) Fault[J]. Northwestern Seismological Journal, 1995, 17(4): 1~12.

    [8] 王啟鳴.三維有限元方法估算郯城8.5級大震重復周期[J].中國地震,1988,4(3):54~58.

    WANG Qi-ming. The calculation of recurrence intervals of theMs8.5 earthquake in Tancheng by the three dimension finite element method[J]. Earthquake Research in China, 1988, 4(3): 54~58.

    [9] 萬永革,吳忠良,周公威,等.地震應力觸發(fā)研究[J].地震學報,2002,24(5):533~551.

    WAN Yong-ge, WU Zhong-liang, Zhou Gong-wei, et al. Research on seismic stress triggering[J]. Acta Seismologica Sinica, 2002, 24(5): 533~551.

    [10] 魏光興.郯廬帶地震活動性研究[M].北京:地震出版社,1993.

    WEI Guang-xing. Research on the seismic activity of Tanlu earthquake belt[M]. Beijing: Seismological Press, 1993.

    [11] Okada Y. Surface deformation to shear and tensile faults in a halfspace[J]. Bulletin of the Seismological Society of America, 1985, 75 (4): 1135~1154.

    [12] Pollitz F F. Postseismic relaxation theory on a spherical Earth[J]. Bulletin of the Seismological Society of America, 1992, 82(1): 422~453.

    [13] 付廣裕,孫文科.球體位錯理論計算程序的總體設計與具體實現(xiàn)[J].地震,2012,32(2):73~87.

    FU Guang-yu, SUN Wen-ke. Overall design and special strctures of the computing codes for coseismic deformations on a layered sperical earth[J]. Earthquake, 2012, 32(2): 73~87.

    [14] 石耀霖.關于應力觸發(fā)和應力影概念在地震預報中應用的一些思考[J].地震,2001, 21(3):1~7.

    SHI Yao-lin. Sress triggers and stress shadows How to apply these concepts to earthquake prediction[J]. Earthquake, 2001, 21(3): 1~7.

    [15] Wang R, Lorenzo-Martín F, Roth F. PSGRN/PSCMP: A new code for calculating co-and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory[J]. Computers & Geosciences, 2006, 32: 527~541.

    [16] Freed A M. Earthquake triggering by static, dynamic, and postseismic stress transfer[J]. Annual Review of Earth & Planetary Sciences, 2005, 33: 335~367.

    [17] 孫玉軍,董樹文,范桃園,等.從同震和震后形變分析日本東北Mw9.0級大地震對近場地震活動性的影響[J].地球物理學進展,2013,28(3):1131~1139.

    SUN Yu-jun, DONG Shu-wen, FAN Tao-yuan, et al. The effect of TohotuMw9.0 earthquake on the near field seismic activity from coseismic and postseismic deformation[J]. Progress in Geophys, 2013, 28(3): 1131~1139.

    [18] 江國焰.利用庫侖破裂準則評估活動斷層地震危險性[D].武漢:武漢大學,2013.

    JIANG Guo-yan. Assessing the seismic hazards of active faults using Coulomb Failure Criterion[D]. Wuhan: Wuhan University, 2013.

    [19] Okada Y. Internal deformation due to shear and tensile fault in a half space[J]. Bulletin of the Seismological Society of America, 1992, 92: 1018~1040.

    [20] King G C P, Stein R S, Lin J. Static stress changes and the triggering of earthquakes[J]. Bull.seismol.soc.amer, 1994, 84: 935~953.

    [21] 周翠英,刁桂苓,耿杰,等.1668年郯城大地震震源斷層三維特征反演[J].地球物理學進展,2013,28(6):2814~2824.

    ZHOU Cui-ying, DIAO Gui-ling, GENG jie, et al. 3-D characteristics inversion of hypocenter fault-plane of the 1668 Tancheng great earthquake[J]. Progress in Geophys, 2013, 28(6): 2814~2824.

    THE EFFECT OF 1668 TANCHENG M8.5 EARTHQUAKE ON THE SEISMIC ACTIVITY OF THE VICINITY FROM COSEISMIC AND POSTSEISMIC DEFORMATION

    XU Hao, SUN Yu-jun, WU Zhong-hai

    (Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100081, China)

    The Tanlu fault zone is an important active fault zone in east China. Many great earthquakes have occurred along it. The 1668M8.5 Tancheng earthquake is the largest recorded event. With the constraint of surface rupture, we get a coseismic rupture model based on the elastic dislocation theory. In addition, we calculate the coseismic, postseismic deformation and Coulomb stress change of this earthquake with a viscoelastic multilayered model. The results show that the rupture fault of the earthquake is adextral strike-slip fault, which has slight thrust. The coseismic displacement of this event is very large. The accumulated energy is released thoroughly. The coseismic rupture of the 1668 Tancheng earthquake increased the Coulomb stress on the northwest, northeast and southeast faults of Tanlu fault zone. The Coulomb stress increases further due to the postseismic deformation. After 350 years of the Tancheng earthquake, the Coulomb stress changes in Shanxian, Suqian and Rizhaoget to +1bar-+1MPa. When the mantle is imposed different viscosities, the duration time which the deformation and Coulomb changes become stable will be different. However, the stable values of them are generally same.

    1668 TanchengM8.5 earthquake; coseismic and postseismic deformation; Coulomb stress change; viscosity

    1006-6616(2016)03-0568-09

    2016-04-24

    中國地質調查局地質調查項目“長江經濟帶活動構造與區(qū)域地殼穩(wěn)定性調查”(DD20160268)

    徐昊(1993-),男,碩士研究生,研究方向地球動力學。E-mail:xhdky2011@163.com

    孫玉軍(1983-),男,博士,副研究員,主要從事地球動力學方面的研究。E-mail:sunyujunabc@163.com

    P315.5

    A

    猜你喜歡
    郯城同震庫侖
    1976年唐山強震群震后庫侖應力演化及其與2020年古冶5.1級地震的關系
    地震研究(2021年1期)2021-04-13 01:04:46
    My Mother Taught Me About……
    山東郯城:多肉種植成花卉苗木產業(yè)“新寵”
    云南思茅大寨井水位地震同震響應特征分析*
    地震研究(2018年4期)2018-11-23 02:29:36
    2種模式下2001—2016年郯城主要糧食作物氣象適宜度對比分析
    基于粘彈庫侖應力變化的后續(xù)最大地震震級估計及2008、2014年于田2次7.3級地震之間關系的討論
    中國地震(2015年1期)2015-11-08 11:11:18
    蘆山地震前后介質波速變化與GPS應變場相關性研究?
    蘆山Ms7.0地震引起的水位同震響應特征分析
    四川地震(2014年2期)2014-12-02 04:16:30
    一種周期庫侖作用勢優(yōu)化法的改進
    計算物理(2014年1期)2014-03-11 17:01:03
    長程庫侖勢對高溫超導渦旋電荷的影響
    97人妻精品一区二区三区麻豆| 女的被弄到高潮叫床怎么办| 97超碰精品成人国产| 黄色日韩在线| 免费观看精品视频网站| 成人三级黄色视频| 男女下面进入的视频免费午夜| 美女黄网站色视频| 国国产精品蜜臀av免费| 十八禁网站免费在线| 国产男人的电影天堂91| 国产私拍福利视频在线观看| 亚洲av.av天堂| 五月玫瑰六月丁香| 久久久精品94久久精品| 欧美日韩乱码在线| 国产精品嫩草影院av在线观看| 乱系列少妇在线播放| 美女高潮的动态| 午夜福利成人在线免费观看| 久久精品国产亚洲av天美| 亚洲人成网站高清观看| 日日摸夜夜添夜夜添av毛片| 欧美激情久久久久久爽电影| 精品一区二区三区视频在线观看免费| 亚洲欧美日韩无卡精品| 搡老妇女老女人老熟妇| 中文字幕免费在线视频6| 一个人看视频在线观看www免费| 久久久午夜欧美精品| 九九久久精品国产亚洲av麻豆| 最新中文字幕久久久久| av天堂在线播放| 在线观看av片永久免费下载| 不卡一级毛片| 国内精品宾馆在线| 久久精品综合一区二区三区| 亚洲七黄色美女视频| 中文在线观看免费www的网站| 国产片特级美女逼逼视频| 亚洲av.av天堂| 国产男靠女视频免费网站| 国产 一区精品| 久久久久久久久大av| 色综合亚洲欧美另类图片| 亚洲最大成人av| 久久久色成人| 高清毛片免费看| 日韩欧美在线乱码| 国产高潮美女av| 精品乱码久久久久久99久播| 免费看美女性在线毛片视频| 一边摸一边抽搐一进一小说| 午夜福利视频1000在线观看| 麻豆成人午夜福利视频| 国产精品一区www在线观看| 国产成人aa在线观看| 99riav亚洲国产免费| 最近的中文字幕免费完整| 校园人妻丝袜中文字幕| 尾随美女入室| 大香蕉久久网| 免费无遮挡裸体视频| 搡老熟女国产l中国老女人| 黄色欧美视频在线观看| 亚洲丝袜综合中文字幕| 欧美+日韩+精品| 全区人妻精品视频| av在线观看视频网站免费| 99久久精品热视频| a级毛片免费高清观看在线播放| 国产精品三级大全| 男女做爰动态图高潮gif福利片| 激情 狠狠 欧美| 国产精品三级大全| 狂野欧美白嫩少妇大欣赏| 在线国产一区二区在线| 亚洲欧美日韩无卡精品| 精品久久久久久久末码| 精品一区二区三区av网在线观看| 寂寞人妻少妇视频99o| 天堂动漫精品| 国产美女午夜福利| 搞女人的毛片| 在线观看一区二区三区| 最近手机中文字幕大全| 亚洲欧美日韩高清在线视频| 1000部很黄的大片| 深爱激情五月婷婷| 国产久久久一区二区三区| 国产精品1区2区在线观看.| 婷婷六月久久综合丁香| 麻豆国产av国片精品| av专区在线播放| 日韩av不卡免费在线播放| 成人特级av手机在线观看| 三级毛片av免费| 97超碰精品成人国产| 亚洲成人精品中文字幕电影| 亚洲精品影视一区二区三区av| 欧美国产日韩亚洲一区| 成年女人看的毛片在线观看| 露出奶头的视频| a级毛色黄片| 搞女人的毛片| 亚洲av免费在线观看| 亚洲欧美日韩高清专用| 91狼人影院| 免费无遮挡裸体视频| 国产成人freesex在线 | 大香蕉久久网| 欧美一区二区精品小视频在线| 精品熟女少妇av免费看| 午夜福利在线观看免费完整高清在 | 能在线免费观看的黄片| 欧美另类亚洲清纯唯美| 赤兔流量卡办理| 亚洲人成网站在线播放欧美日韩| 国产精品无大码| 91久久精品国产一区二区三区| 搡老妇女老女人老熟妇| 午夜免费激情av| 99在线视频只有这里精品首页| 午夜精品国产一区二区电影 | 成人av一区二区三区在线看| 少妇被粗大猛烈的视频| 国产精品三级大全| 一边摸一边抽搐一进一小说| 国产单亲对白刺激| 午夜视频国产福利| 小说图片视频综合网站| 国产 一区精品| 国产精品爽爽va在线观看网站| 亚洲欧美清纯卡通| 亚洲av电影不卡..在线观看| 看黄色毛片网站| 99久久久亚洲精品蜜臀av| 1000部很黄的大片| av国产免费在线观看| 22中文网久久字幕| 91午夜精品亚洲一区二区三区| 国产精品久久视频播放| 久久午夜亚洲精品久久| 久久午夜亚洲精品久久| 成人二区视频| 99久久精品一区二区三区| av福利片在线观看| 欧美日韩在线观看h| 国产精品三级大全| 免费观看的影片在线观看| 亚洲久久久久久中文字幕| 天堂av国产一区二区熟女人妻| 夜夜看夜夜爽夜夜摸| 久久精品夜夜夜夜夜久久蜜豆| 成人午夜高清在线视频| 久久99热6这里只有精品| 在线播放无遮挡| 久久热精品热| 久久久久久久午夜电影| 精品久久久久久久末码| 看黄色毛片网站| 精品久久久久久久久亚洲| 一本久久中文字幕| 日韩欧美在线乱码| 国产成人91sexporn| 国产精品一区www在线观看| 欧美激情国产日韩精品一区| 天堂动漫精品| 成年av动漫网址| 国产午夜福利久久久久久| 中文字幕熟女人妻在线| 级片在线观看| 欧美日韩在线观看h| 最后的刺客免费高清国语| 精品久久久久久久人妻蜜臀av| 日本色播在线视频| 别揉我奶头~嗯~啊~动态视频| 欧美另类亚洲清纯唯美| 国产伦在线观看视频一区| 亚洲中文字幕日韩| 国产真实伦视频高清在线观看| 哪里可以看免费的av片| 99久国产av精品国产电影| 最新在线观看一区二区三区| 又爽又黄a免费视频| 日日摸夜夜添夜夜添av毛片| 精品人妻偷拍中文字幕| 久久久久久久久久黄片| 啦啦啦啦在线视频资源| 国产高潮美女av| 亚洲av第一区精品v没综合| 国产综合懂色| 婷婷色综合大香蕉| 日本色播在线视频| 狂野欧美白嫩少妇大欣赏| 大香蕉久久网| 变态另类丝袜制服| 欧美国产日韩亚洲一区| 国产 一区精品| 美女cb高潮喷水在线观看| 人人妻人人看人人澡| 久久精品国产自在天天线| 麻豆国产av国片精品| 伊人久久精品亚洲午夜| 亚洲无线在线观看| 天堂网av新在线| 午夜福利18| 男人的好看免费观看在线视频| 中国美白少妇内射xxxbb| 色综合亚洲欧美另类图片| 中文字幕精品亚洲无线码一区| 亚洲经典国产精华液单| 国产成人福利小说| 蜜桃久久精品国产亚洲av| 国产女主播在线喷水免费视频网站 | 在现免费观看毛片| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 国内久久婷婷六月综合欲色啪| 又粗又爽又猛毛片免费看| 亚洲国产精品成人综合色| 成人欧美大片| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久久久av| 日本在线视频免费播放| 亚洲成人中文字幕在线播放| 欧美xxxx性猛交bbbb| 国产黄片美女视频| 在线观看美女被高潮喷水网站| 亚洲av成人精品一区久久| 国产精品无大码| 日韩亚洲欧美综合| 国产精品1区2区在线观看.| 噜噜噜噜噜久久久久久91| 亚洲国产精品国产精品| 日韩中字成人| 国产亚洲欧美98| 亚洲人成网站在线播放欧美日韩| 我要搜黄色片| 少妇的逼水好多| 久久久久久久久大av| 欧美丝袜亚洲另类| 99久久久亚洲精品蜜臀av| 夜夜爽天天搞| 午夜精品一区二区三区免费看| 精品熟女少妇av免费看| 免费电影在线观看免费观看| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲网站| 亚洲人成网站高清观看| 亚洲国产精品合色在线| 日韩中字成人| 六月丁香七月| 久久久久精品国产欧美久久久| 九九在线视频观看精品| 久久久久久久久大av| 一本久久中文字幕| 啦啦啦观看免费观看视频高清| 国产爱豆传媒在线观看| a级毛片a级免费在线| 成人鲁丝片一二三区免费| 亚洲欧美日韩东京热| 久久精品国产亚洲av香蕉五月| 俺也久久电影网| 国产色婷婷99| 男人舔奶头视频| 国产精品精品国产色婷婷| 国产精品伦人一区二区| 亚洲专区国产一区二区| 天堂动漫精品| 成年免费大片在线观看| 色5月婷婷丁香| 国产精品福利在线免费观看| 精品久久国产蜜桃| 99热6这里只有精品| 免费av观看视频| 无遮挡黄片免费观看| 日日干狠狠操夜夜爽| 在线播放国产精品三级| 三级经典国产精品| 国产私拍福利视频在线观看| av女优亚洲男人天堂| 麻豆一二三区av精品| av天堂在线播放| 两个人视频免费观看高清| 99久久成人亚洲精品观看| 国内精品美女久久久久久| 噜噜噜噜噜久久久久久91| 欧美另类亚洲清纯唯美| 午夜福利在线观看免费完整高清在 | 大型黄色视频在线免费观看| 最近中文字幕高清免费大全6| 色尼玛亚洲综合影院| 人人妻人人澡欧美一区二区| aaaaa片日本免费| 99久久久亚洲精品蜜臀av| 亚洲va在线va天堂va国产| 伊人久久精品亚洲午夜| 一个人免费在线观看电影| 午夜视频国产福利| 久久精品综合一区二区三区| 日日摸夜夜添夜夜添小说| 久久精品国产清高在天天线| 九九爱精品视频在线观看| 国产一区二区三区av在线 | 人人妻人人看人人澡| 国产av不卡久久| 99久久成人亚洲精品观看| 少妇被粗大猛烈的视频| 亚洲av电影不卡..在线观看| 黄片wwwwww| 干丝袜人妻中文字幕| 亚洲精品日韩在线中文字幕 | 国产精品电影一区二区三区| 亚洲18禁久久av| 天堂影院成人在线观看| 国产私拍福利视频在线观看| 乱码一卡2卡4卡精品| 热99在线观看视频| 别揉我奶头 嗯啊视频| 18禁黄网站禁片免费观看直播| 国产日本99.免费观看| 久久久久久九九精品二区国产| 毛片女人毛片| 国产一级毛片七仙女欲春2| 性色avwww在线观看| 日韩精品中文字幕看吧| 免费搜索国产男女视频| 99久久九九国产精品国产免费| 国产黄色小视频在线观看| 人人妻人人看人人澡| 有码 亚洲区| 日本一本二区三区精品| 精品午夜福利在线看| 成人亚洲精品av一区二区| 欧美三级亚洲精品| 国产精品三级大全| 国产精品美女特级片免费视频播放器| 亚洲丝袜综合中文字幕| 少妇被粗大猛烈的视频| 亚洲精品久久国产高清桃花| 秋霞在线观看毛片| 伊人久久精品亚洲午夜| 九九久久精品国产亚洲av麻豆| 国内精品宾馆在线| 91在线观看av| 国产成人91sexporn| 99久国产av精品| 午夜影院日韩av| 非洲黑人性xxxx精品又粗又长| 国产爱豆传媒在线观看| 午夜激情福利司机影院| 日本五十路高清| 日韩 亚洲 欧美在线| 日本黄大片高清| 国内精品宾馆在线| 日韩欧美 国产精品| 好男人在线观看高清免费视频| 国产熟女欧美一区二区| 春色校园在线视频观看| 免费观看在线日韩| 午夜精品在线福利| 国产爱豆传媒在线观看| 秋霞在线观看毛片| 亚洲自拍偷在线| 欧美丝袜亚洲另类| 成人三级黄色视频| 亚洲人成网站在线播| 97人妻精品一区二区三区麻豆| 成人二区视频| 午夜a级毛片| 亚洲美女搞黄在线观看 | 又黄又爽又刺激的免费视频.| 一进一出好大好爽视频| 久久久久久久亚洲中文字幕| 国产男人的电影天堂91| 看片在线看免费视频| 成年女人永久免费观看视频| 国产精品不卡视频一区二区| 欧美一区二区亚洲| 又黄又爽又免费观看的视频| 久久韩国三级中文字幕| 国产视频一区二区在线看| 欧美xxxx黑人xx丫x性爽| 国语自产精品视频在线第100页| 六月丁香七月| 精品熟女少妇av免费看| 国产精品一区www在线观看| 日韩人妻高清精品专区| 桃色一区二区三区在线观看| 国产三级在线视频| 成熟少妇高潮喷水视频| 99热网站在线观看| 国产精品一区二区免费欧美| 精品久久久久久久久av| 中国美白少妇内射xxxbb| 我要搜黄色片| 三级经典国产精品| 一边摸一边抽搐一进一小说| 身体一侧抽搐| 亚洲成人久久爱视频| 久久欧美精品欧美久久欧美| 久久久国产成人精品二区| 无遮挡黄片免费观看| 精品午夜福利视频在线观看一区| 少妇丰满av| 国产精品永久免费网站| 97超视频在线观看视频| 啦啦啦韩国在线观看视频| 欧美一区二区国产精品久久精品| 国产伦在线观看视频一区| а√天堂www在线а√下载| 成年版毛片免费区| 久久国内精品自在自线图片| 国内精品美女久久久久久| 国产真实乱freesex| 日韩精品青青久久久久久| 色视频www国产| 深夜a级毛片| 精品久久久久久成人av| 久久久久久伊人网av| 精品久久久久久久久久免费视频| 国产黄a三级三级三级人| 麻豆一二三区av精品| 18禁在线播放成人免费| 亚洲色图av天堂| 亚洲va在线va天堂va国产| 三级国产精品欧美在线观看| 亚洲av熟女| 精品久久久噜噜| 亚洲精品日韩av片在线观看| 最近2019中文字幕mv第一页| 国产高清不卡午夜福利| 嫩草影院精品99| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕 | 一本久久中文字幕| 一个人看视频在线观看www免费| 日本a在线网址| 久久久久国产精品人妻aⅴ院| 国模一区二区三区四区视频| 成人av一区二区三区在线看| 最新在线观看一区二区三区| 婷婷六月久久综合丁香| 亚洲三级黄色毛片| 美女xxoo啪啪120秒动态图| 中国美白少妇内射xxxbb| 精品久久久久久久久av| 午夜久久久久精精品| 3wmmmm亚洲av在线观看| 国产人妻一区二区三区在| 亚洲成人av在线免费| 亚洲av第一区精品v没综合| a级一级毛片免费在线观看| 91av网一区二区| 成年av动漫网址| 99热只有精品国产| 免费av观看视频| av视频在线观看入口| 变态另类成人亚洲欧美熟女| 亚洲国产精品成人久久小说 | 18禁在线播放成人免费| 老熟妇乱子伦视频在线观看| 日韩亚洲欧美综合| 色在线成人网| 国产乱人偷精品视频| 少妇的逼好多水| 亚洲性久久影院| 18禁在线无遮挡免费观看视频 | 非洲黑人性xxxx精品又粗又长| 91av网一区二区| 级片在线观看| 丰满乱子伦码专区| 国产伦一二天堂av在线观看| 久久韩国三级中文字幕| 免费在线观看影片大全网站| 国产视频内射| 一级a爱片免费观看的视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲内射少妇av| 特级一级黄色大片| 日韩欧美精品免费久久| 欧美+亚洲+日韩+国产| 亚洲av不卡在线观看| 97超视频在线观看视频| 日韩成人伦理影院| 欧美成人一区二区免费高清观看| 日韩欧美 国产精品| 在线看三级毛片| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av成人精品一区久久| 嫩草影院精品99| 三级经典国产精品| 能在线免费观看的黄片| 国产一区二区在线观看日韩| 国产色婷婷99| 日本a在线网址| 亚洲国产高清在线一区二区三| 精品无人区乱码1区二区| 久久精品夜夜夜夜夜久久蜜豆| 人妻丰满熟妇av一区二区三区| 国产精品爽爽va在线观看网站| 欧美高清成人免费视频www| 一夜夜www| 国产免费男女视频| 亚洲精品国产av成人精品 | 亚洲国产精品久久男人天堂| 亚洲精品粉嫩美女一区| 欧美激情久久久久久爽电影| 国产淫片久久久久久久久| 成人三级黄色视频| 成人美女网站在线观看视频| 99九九线精品视频在线观看视频| 在线观看免费视频日本深夜| aaaaa片日本免费| 日韩欧美一区二区三区在线观看| 中文字幕av在线有码专区| 夜夜爽天天搞| 日韩一本色道免费dvd| 欧美成人免费av一区二区三区| 亚洲国产欧洲综合997久久,| 天天一区二区日本电影三级| 好男人在线观看高清免费视频| 日韩欧美精品v在线| 国产色爽女视频免费观看| 18禁在线播放成人免费| 欧美日韩一区二区视频在线观看视频在线 | 性插视频无遮挡在线免费观看| 熟女电影av网| 亚洲一区高清亚洲精品| 国产一区亚洲一区在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成av人片在线播放无| 有码 亚洲区| 国产精品电影一区二区三区| 床上黄色一级片| 联通29元200g的流量卡| 国产亚洲91精品色在线| 成人美女网站在线观看视频| 女生性感内裤真人,穿戴方法视频| 国产av一区在线观看免费| 无遮挡黄片免费观看| 免费看av在线观看网站| 久久热精品热| 最近的中文字幕免费完整| 国产三级中文精品| 一进一出抽搐动态| 天天躁夜夜躁狠狠久久av| 如何舔出高潮| 精品久久国产蜜桃| 波野结衣二区三区在线| 午夜激情欧美在线| 久久精品国产清高在天天线| 国产伦一二天堂av在线观看| 一本一本综合久久| 久久精品人妻少妇| 尤物成人国产欧美一区二区三区| 听说在线观看完整版免费高清| 熟女电影av网| 男人狂女人下面高潮的视频| 日本精品一区二区三区蜜桃| 看免费成人av毛片| 麻豆国产av国片精品| 久久精品久久久久久噜噜老黄 | 日本在线视频免费播放| 黄色视频,在线免费观看| 亚洲中文字幕日韩| 极品教师在线视频| 在线免费观看的www视频| 成人永久免费在线观看视频| 国产欧美日韩一区二区精品| 久久精品夜夜夜夜夜久久蜜豆| 久久久欧美国产精品| 久久久久九九精品影院| 亚洲精品日韩av片在线观看| 亚洲丝袜综合中文字幕| 久久午夜福利片| 午夜影院日韩av| 黄色欧美视频在线观看| 色播亚洲综合网| 亚洲av中文av极速乱| 久久人妻av系列| 国产高清有码在线观看视频| 日韩av不卡免费在线播放| 中文字幕人妻熟人妻熟丝袜美| 在线观看美女被高潮喷水网站| 亚洲成人久久性| 直男gayav资源| 欧美一区二区亚洲| 免费在线观看影片大全网站| 国产久久久一区二区三区| 黄色视频,在线免费观看| 亚洲一区二区三区色噜噜| 内地一区二区视频在线| 国产精品精品国产色婷婷| av在线蜜桃| 精品99又大又爽又粗少妇毛片| a级毛片a级免费在线| 免费观看在线日韩| 中文字幕免费在线视频6| 国产一区二区三区在线臀色熟女| 午夜福利在线观看免费完整高清在 | 99久久成人亚洲精品观看| 亚洲av一区综合| ponron亚洲| 国产aⅴ精品一区二区三区波| 一夜夜www| 日韩精品青青久久久久久| 亚洲人成网站高清观看| 卡戴珊不雅视频在线播放| av中文乱码字幕在线| 联通29元200g的流量卡| 亚洲精品久久国产高清桃花| av国产免费在线观看|