李路偉,余中元,陳柏旭,鄭榮熒,王世元,趙 錢
(1.防災(zāi)科技學(xué)院,河北 三河 065201;2.北京防災(zāi)科技有限公司,北京 101100;3.四川省地震局,四川 成都 610041)
2022年9月5日,四川省甘孜藏族自治州瀘定縣發(fā)生了強烈地震。經(jīng)中國地震臺網(wǎng)測定,本次 地 震 震 級 MS6.8,微 觀 震 中(29.59°N,102.08°E)位于磨西鎮(zhèn)西南約12km,震源深度約16km[1]。截至2022年9月11日17時,地震共造成93人遇難,25人失蹤,11余萬人受災(zāi),5萬余間房屋受損。地震引發(fā)的滑坡等地質(zhì)災(zāi)害極為嚴(yán)重,造成大量人員傷亡、房屋建筑和基礎(chǔ)設(shè)施受損。地震發(fā)生后,應(yīng)急管理部中國地震局組織專家對災(zāi)區(qū)200個調(diào)查點展開了實地調(diào)查,通過參考震區(qū)斷裂構(gòu)造背景、儀器烈度、余震分布、震源機(jī)制和無人機(jī)遙感解譯等方法,結(jié)合強震動觀測記錄,確定了此次地震的烈度分布:最高烈度為Ⅸ度,等震線長軸呈北西走向,長軸195km,短軸112km。其中,Ⅵ度區(qū)及以上面積19 089km2,Ⅸ度區(qū)面積約280km2,主要影響區(qū)域涉及甘孜藏族自治州瀘定縣磨西鎮(zhèn)、得妥鎮(zhèn)、燕子溝鎮(zhèn)、德威鎮(zhèn),以及雅安市石棉縣王崗坪彝族藏族鄉(xiāng)、草科藏族鄉(xiāng)、新民藏族彝族鄉(xiāng),共計7個鄉(xiāng)鎮(zhèn)[2]。根據(jù)余震分布刻畫的同震斷層三維模型,判斷發(fā)震斷層為鮮水河斷裂帶南端的磨西段[3]。此次地震初始破裂點深度約12.8km,是一次高傾角左旋走滑型地震,發(fā)震斷層西傾,走向160°,傾角80°,最大滑動量位于震中附近約184cm,破裂持續(xù)時間約18s[4]。震源輻射地震能量的效率偏高,造成地面震動和建筑破壞的能力較強[4]。從震源區(qū)地殼速度結(jié)構(gòu)來看,此次地震發(fā)生在地震波速度的高、低速過渡帶上,體現(xiàn)了青藏高原東部與四川盆地之間的地殼結(jié)構(gòu)差異對區(qū)域應(yīng)力積累的控制作用[4]。
鮮水河斷裂位于青藏高原東南緣,是一條分割巴顏喀拉塊體和川滇菱形塊體的重要邊界構(gòu)造,為左旋走滑運動的全新世活動斷裂,與NE走向的龍門山斷裂帶、NW 走向的安寧河斷裂帶交匯構(gòu)成了川西地區(qū)著名的“Y”字形斷裂帶[5-7](圖1a),在青藏高原東部及鄰區(qū)的構(gòu)造變形過程中起到了重要的控制作用[8]。鮮水河斷裂北西起于甘孜西北至治多一帶,向南東經(jīng)爐霍、道孚、乾寧、康定、瀘定、磨西,至石棉新民以南活動形跡逐漸減弱,最終消失于石棉公益海附近,向東南繼續(xù)分為多個分支斷層。斷裂橫切了松潘-甘孜造山帶的主體,系造山運動后期陸內(nèi)變形的產(chǎn)物[9]。該斷裂帶具有規(guī)模大、活動強、地震頻度高等特點[10-14]。自1725年有地震記錄以來,共發(fā)生6級以上地震23次,中強震地表破裂在空間分布上幾乎覆蓋了整條斷裂[15]。
圖1 研究區(qū)地震地質(zhì)構(gòu)造背景Fig.1 Seismic geological tectonic setting in the study area
研究區(qū)(圖1c)地處以鮮水河-安寧河-小江斷裂為界的川滇斷塊與涼山斷塊結(jié)合部位SN向構(gòu)造帶北段,為多組構(gòu)造帶的交匯復(fù)合部位,位于晉寧-澄江期花崗巖體之上。主要出露地層為前中生界、古近系、新近系和第四系,少量印支-燕山期侵入巖(圖1b)。前中生界可見鉀長花崗巖、二長花崗巖、混合質(zhì)花崗巖等元古代花崗巖;二疊系包含石英巖、角閃片巖;三疊系包含砂巖、粉砂巖等。第四系多為沖洪積、坡積黏土、砂土和礫石,少量分布在河谷兩岸階地上[16]。
本次地震的發(fā)震斷層為鮮水河斷裂南端的磨西段,位于鮮水河、安寧河和龍門山斷裂的交匯區(qū),為全新世強烈活動的斷裂帶。磨西段位于大渡河西側(cè)、貢嘎山東側(cè),北起雅家梗,經(jīng)二臺子、灣東、猛虎崗、田灣,南至安順場一帶。斷裂北端海拔3900m,向南至安順場降低至1000m左右[16-19]。前人研究結(jié)果表明,磨西段平均水平滑動速率為6.0~9.9mm/a[17]。聞學(xué)澤等[20]在躍近坪附近測得同震位錯為(4.5±0.5)m,結(jié)合磨西斷裂(8±2)mm/a的平均滑動速率[16,17],估計的強震復(fù)發(fā)間隔為400~500a;周榮軍等[17]二臺子探槽研究結(jié)果表明,鮮水河斷裂磨西段的強震復(fù)發(fā)間隔大致在 300a左右。綜合前人資料[8,18-22],磨西段的強震復(fù)發(fā)間隔約束為300~500a[16]。
本文大比例尺野外調(diào)查通過遙感解譯、探槽古地震和無人機(jī)航空攝影測量等方法,確定了地表破裂、滑坡等同震地表變形特征和和滑坡分布,對地表同震變形的幾何分布、破裂長度、最大位錯量等地表變形特征和地震的影響范圍等進(jìn)行調(diào)查。該研究結(jié)果有望為甘孜-玉樹-鮮水河斷裂帶的地震危險性研究提供新的數(shù)據(jù)支撐,有助于深刻認(rèn)識川滇地塊和巴顏喀拉地塊向東擠出的過程和亞歐板塊與印度板塊的擠壓碰撞效應(yīng)。
為了查明鮮水河斷裂磨西段同震破裂的空間展布特征,根據(jù)前人研究資料[18-20,23],利用奧維互動地圖進(jìn)行遙感影像解譯,確定了磨西段的幾何位置。在此基礎(chǔ)上進(jìn)行野外地質(zhì)調(diào)查,調(diào)查范圍北起老榆林村,經(jīng)磨西鎮(zhèn)、猛虎崗,灣東村、幸福村、團(tuán)結(jié)村、安順場鎮(zhèn),南至紫馬跨,沿線全長約120km。野外地質(zhì)調(diào)查對象主要為本次地震的同震破裂和滑坡災(zāi)害等,調(diào)查內(nèi)容包括同震地表破裂識別、典型破裂點定量測量、滑坡災(zāi)害分布范圍和規(guī)模等,結(jié)合探槽解譯和無人機(jī)航空攝影測量,旨在確定破裂長度和大型連續(xù)滑坡的分布范圍等同震地表變形特征。
通過遙感影像解譯,在調(diào)查區(qū)北、南兩端的老榆林村、紫馬跨分別發(fā)現(xiàn)了大規(guī)模的斷層槽谷,垂直于斷層走向開挖2個大型探槽。老榆林村調(diào)查點探槽剖面長約8m,高約2m;紫馬跨調(diào)查點探槽剖面長約13m,高約3.5m,并進(jìn)行掛線處理,網(wǎng)格面積1m2。
利用無人機(jī)對構(gòu)造地貌進(jìn)行精細(xì)測量。無人機(jī)采用大疆“精靈4RTK”,通過“攝影測量2D”規(guī)劃航線,所拍攝的高清圖片數(shù)據(jù)通過Agisoft PhotoScan進(jìn)行后處理,包括照片對齊、建立密集點云、生成網(wǎng)格、生成紋理等操作步驟,最終生成高精度DEM數(shù)據(jù)。老榆林村調(diào)查點無人機(jī)共飛行2架次,拍攝照片共112張,飛行高度500m,航向重疊度80%、旁向重疊度80%,作業(yè)面積663194.6m2,生成高精度DEM數(shù)據(jù)分辨率為0.91m;紫馬跨調(diào)查點無人機(jī)共飛行2架次,拍攝照片共326張,飛行高度230m,航向重疊度70%、旁向重疊度80%,作業(yè)面積878618.7m2,生成高精度DEM 數(shù)據(jù)分辨率為0.59m。利用Arc GIS10.7軟件對DEM 數(shù)據(jù)進(jìn)行山體陰影、生成等值線等后處理,并對數(shù)據(jù)結(jié)果進(jìn)行遙感解譯,以研究本次地震的同震變形。
2.1.1 典型破裂調(diào)查點及特征
燕子溝調(diào)查點位于燕子溝北側(cè),裂隙沿NWW 向延伸(圖2a),長度約300m,可見破裂垂直位移量(10±0.5)cm,線性特征明顯。因其穿過大面積的耕地,且受到地形制約,可觀察到的破裂分布零散,連續(xù)性差。
磨西鎮(zhèn)北側(cè)調(diào)查點發(fā)育線性裂縫(圖2b),沿NW 向延伸約130m,裂隙規(guī)模較小,受到植物覆蓋和地形起伏的影響,未向兩側(cè)追索,推測仍有裂縫呈線性發(fā)育。
圖2 瀘定M S6.8地震典型的同震地表變形特征Fig.2 Photos of the typical co-seismic surface deformation induced by the Luding M S6.8 Earthquake
猛虎崗調(diào)查點破裂表現(xiàn)為張性裂縫(圖2c~e),沿NNW 向線性展布,垂直位移量(15±1)cm,寬度為(7±0.5)cm(圖2c)。由此向東北50m,長約30cm的礫石被錯斷,表面發(fā)育擦痕,斷面平整且平行于破裂的延伸方向(圖2d)。由此向南西50m,巖石斷面發(fā)育擦痕(圖2e),擦痕傾向SE,傾角60°,斷面上表現(xiàn)為一組均勻的平行細(xì)紋,細(xì)紋上端粗且深,下端細(xì)而淺。
幸福村調(diào)查點破裂表現(xiàn)為右階雁裂式裂縫(圖2f),主破裂沿著NW 向呈線性延伸,伴隨著多個小型次級破裂,與主破裂構(gòu)成約45°的夾角。主破裂垂直位移量為(4±1)cm,形成小型陡坎。次級破裂可見2~3組,沿NE走向的平行展布,每組間隔2~3m。
2.1.2 典型滑坡調(diào)查點及特征
燕子溝調(diào)查點(圖3a)位于燕子溝景區(qū)內(nèi),大小滑坡呈連片分布,多為土質(zhì)滑坡?;碌钠矫嫘螒B(tài)主要以舌形為主,滑坡前緣位于河溝內(nèi),滑坡后壁的平面形態(tài)多呈現(xiàn)為弧形,開口方向與滑動方向基本一致,坡度在45°左右。坡腳前緣僅有少量堆積物,結(jié)合道路旁護(hù)坡攔截的落石的規(guī)模和數(shù)量,推測此次地震在燕子溝內(nèi)造成的滑坡規(guī)模中等。野外調(diào)查情況表明,自此向西滑坡現(xiàn)象不明顯。
磨西鎮(zhèn)調(diào)查點(圖3b)位于磨西鎮(zhèn)北西方向約3km處的河谷中,滑坡在河谷兩側(cè)斜坡上呈點狀分布,滑坡類型為淺層滑坡。該滑坡體主要由砂礫、碎石混粉黏粒組成,其結(jié)構(gòu)松散,力學(xué)強度低且遇水易軟化,在地震作用擾動下易形成滑坡?;旅婵傮w呈半圓形,滑坡后壁的平面形態(tài)為圈椅狀,滑坡前緣除了堆積物之外,在坡腳前不遠(yuǎn)處還堆積有礫石和碎石。
圖3 瀘定M S6.8地震典型滑坡Fig.3 Photos of typical landslides induced by the Luding M S6.8 Ear thquake
得妥鎮(zhèn)調(diào)查點(圖3c)位于得妥鎮(zhèn)北西約5km處,此處滑連片分布,滑坡類型為淺層-中淺層滑坡。該處滑坡平面形態(tài)特征呈不規(guī)則的長條狀,滑坡后壁均沿著山脊分布,開口方向與滑動方向基本相同,坡度較陡,運動路徑狹長,滑坡源通常位于坡頂。由于此處坡度較大,通道狹窄,所以堆積物在到達(dá)坡腳后會繼續(xù)向遠(yuǎn)方運動。由此沿著大渡河向北東方向前進(jìn),很長一段距離內(nèi)并未出現(xiàn)滑坡,而后大渡河兩岸斜坡上又有滑坡發(fā)育。
二臺子調(diào)查點(圖3d)位于二臺子南東約400m處,滑坡沿斷層槽谷連片發(fā)育,滑坡平面特征以舌形為主,滑坡后壁的平面形態(tài)多呈現(xiàn)弧形,開口方向與滑坡方向基本一致,坡度較大,多在45°~80°,滑坡前緣一般位于坡腳和公路內(nèi)側(cè)。在坡腳的堆積物中可以看到斜坡上原本生長的植被,滑坡后壁附近可以看到傾倒的林木,多呈“醉漢林”狀發(fā)育。
幸福村調(diào)查點(圖3e)位于幸福村東側(cè)約400m處,與得妥鎮(zhèn)滑坡體有相似之處,運動路徑狹窄,且平面形態(tài)特征均為不規(guī)則的長條狀,地表形貌以滑坡群為主要特征。此外,也發(fā)育著一些半圓形的滑坡面和弧形滑坡后壁,呈圈椅狀,共同構(gòu)成小型滑坡。
愛國村調(diào)查點(圖3f)位于愛國村入口處附近,可見大渡河西岸發(fā)育滑坡群,多為淺層-滾動型滑坡。在劇烈的震動下,斜坡滑坡體沿著坡體軟弱結(jié)構(gòu)面以滾動運動的方式堆積到坡腳部位,堆積物破碎嚴(yán)重,礫石粒徑大小不一[26]。
2.2.1 同震地表變形特征
本次地震的地表變形主要分布在貢嘎山東麓,海拔1300~2700m。同震地表變形以地裂縫和滑坡為主要特征。調(diào)查結(jié)果顯示,破裂整體沿鮮水河斷裂磨西段呈NNW 方向分布,形態(tài)單一,規(guī)模小且分布范圍窄,破裂寬度僅2~8cm,最大垂直位移僅(15±1)cm。破裂帶的組合特征較簡單,主要由一系列不連續(xù)的地裂縫為主,伴隨著擦痕、巖石位錯等構(gòu)造現(xiàn)象。破裂規(guī)模自猛虎崗向南北兩側(cè)遞減,主要表現(xiàn)在裂縫垂直位移量、寬度和構(gòu)造變形強度等方面。猛虎崗附近,破裂的同震變形特征明顯,規(guī)模最大,可見擦痕、巖石裂縫等同震變形現(xiàn)象,并發(fā)現(xiàn)了最大垂直位移。需要注意的是,調(diào)查點沿線地勢陡峭,并發(fā)育大量滑坡,因此裂縫的垂直位移量不能完全排除是構(gòu)造變形與滑坡共同作用的結(jié)果。
滑坡分布范圍北起燕子溝,經(jīng)磨西鎮(zhèn)、猛虎崗、灣東村、幸福村、團(tuán)結(jié)村,終止于安順場鎮(zhèn),總計全長約60km。其中,以磨西鎮(zhèn)為中心,分別向北西15km和南東25km最為嚴(yán)重。根據(jù)滑坡的幾何分布范圍,可分為東、西兩部分,整體分布受斷裂帶和大渡河水系的空間展布影響較大。西部滑坡以燕子溝至愛國村一帶最為嚴(yán)重,滑坡體規(guī)模大,范圍廣,整體呈NNW 向,平行于鮮水河斷裂磨西段延伸約40km,可能代表著鮮水河斷裂磨西段對于滑坡的控制作用;東部滑坡區(qū)主要分布于得妥鎮(zhèn)以北約5km處至瀘定縣城西南端,滑坡體沿大渡河兩岸呈NE向展布,延伸約10km,與NNW 向沿磨西段分布的滑坡共同構(gòu)成“Y”字型分布,長約40km。造成這種現(xiàn)象的原因,筆者分析認(rèn)為,可能主要受控于斷裂帶和大渡河水系的影響,同時還可能受到滑坡體的巖土破碎程度、力學(xué)性質(zhì)等多種物理因素的共同影響[24-25]。
2.2.2 破裂長度的討論
綜合野外調(diào)查結(jié)果,認(rèn)為此次地震的同震破裂向北終止于燕子溝,向南終止于幸福村,全長約35km。此外,我們在老榆林村東側(cè)山坡(康定市以南約5km)和紫馬跨分別進(jìn)行了遙感影像解譯和野外地質(zhì)調(diào)查,可見斷層槽谷延伸約2~3km。垂直于斷層走向清理了2個大型探槽,同時進(jìn)行無人機(jī)地貌精細(xì)測量以研究微地貌的變形特征,旨在揭示最新一次的地震事件,為確定本次地震同震破裂的終點位置提供依據(jù)。
老榆林村東側(cè)山坡的斷層槽谷沿NNW 向呈線性延伸約3km,寬度約210m,兩條斷層呈右階雁裂狀展布(圖4a),斷層斷錯沖溝,沖溝左旋位錯量分別為(172±3)m和(326±5)m(圖4a),反映了斷層的多期次左旋剪切運動。探槽剖面揭示了4條高角度逆斷層和8套地層單元(圖4b~e、表1)。斷層產(chǎn)狀分別為84°∠83°、79°∠80°、82°∠89°和85°∠86°。最新一次地震事件發(fā)生于U2沉積之后,U1沉積之前,發(fā)震斷層為F1。U1保存完好,未發(fā)現(xiàn)瀘定MS6.8地震同震地表破裂的證據(jù)。
紫馬跨村北側(cè)斷層槽谷沿NNE向延伸2~3km,寬度約260m。探槽剖面揭示了3條正斷層、1條逆斷層、1條裂縫和10套地層單元(圖5b~h、表2),斷層產(chǎn)狀分別為88°∠83°、87°∠43°、84°∠35°、89°∠64°和83°∠85°。最新一次地震事件發(fā)生于U2、U8沉積之后,U1沉積之前。U1保存完好,未發(fā)現(xiàn)本次瀘定MS6.8地震同震地表破裂的證據(jù)。
表2 紫馬跨探槽剖面地層描述Tab.2 Description of character istics of stratums of the trench profile at Zimakua
圖5 紫馬跨槽谷地貌及探槽剖面Fig.5 Valley landfor m and trench profile at the Zimakua
因此,綜合野外資料,瀘定地震的同震破裂以小型裂縫為主,北起燕子溝,經(jīng)磨西鎮(zhèn)、二臺子、猛虎崗,南至幸福村,全長約35km。
2.2.3 對甘孜-玉樹-鮮水河斷裂帶的地震危險性啟示
現(xiàn)有資料表明,鮮水河斷裂強震頻繁,中強震事件在空間分布上幾乎覆蓋了整條斷裂[8,18-22],7級以上地震的強震危險性較高。因此,瀘定MS6.8地震之后,大家普遍關(guān)心下一次的強震可能發(fā)生在何地,震級可能為多大。
陳桂華等[10]通過大比例尺野外填圖、探槽開挖等地質(zhì)調(diào)查方法,認(rèn)為1786年磨西M強震的地表破裂帶分布于康定至田灣,這與鮮水河斷裂磨西段高度吻合。因此,鮮水河斷裂磨西段的離逝時間為236a。前人對鮮水河斷裂開展了大量的滑動速率研究,如徐錫偉等[27]通過位錯地貌的定量研究和14C樣品結(jié)果,得到了鮮水河斷裂南東段平均水平滑動速率為(9.6±1.7)mm/a;Chen等[28]對不同時間尺度的位錯地貌進(jìn)行定量研究,利用最小二乘法對滑動速率進(jìn)行擬合,得到鮮水河斷裂東南段左旋走滑速率為9.3 mm/a;Zhang[29]利用上階地廢棄年代限定錯離河道一側(cè)陡坎位移的起始年齡,提出了鮮水河斷裂主斷層的平均滑動速率為8~11 mm/a。綜合現(xiàn)有資料分析認(rèn)為,自1786年磨西地震之后,產(chǎn)生的水平位移量的理論值應(yīng)為2.19~2.48m(表3)。依據(jù)鄧起東等[13]建立的青藏高原走滑斷層的水平位錯量(D)與震級(M)的經(jīng)驗關(guān)系式:M=7.13+0.68log(D),估算震級為7.4級(表3);如果利用Wells等[9]建立的走滑斷層的平均位移量(D)與震級(M)的經(jīng)驗關(guān)系式:M=7.04+0.89log(D),估算震級為7.3~7.4級(表3)。經(jīng)驗公式推算出的理論震級顯然高于公布的6.8級,這可能暗示本次地震未完全釋放所積累的應(yīng)力。
表3 震級與水平位錯量的經(jīng)驗關(guān)系Tab.3 Relationship between magnitude and horizontal dislocation
甘孜-玉樹斷裂帶和鮮水河斷裂帶均受到青藏高原東部及鄰區(qū)的構(gòu)造地塊所控制,具備相似的構(gòu) 造 背 景。綜 合 前 人 資 料 可 知[17,18,27-31],甘孜-玉樹斷裂帶的甘孜段、馬尼干戈段、鄧柯段和玉樹段,最新一次歷史地震分別發(fā)生于1866年、1854年、1896年和2010年(圖6),最長的離逝時間為168a,最短僅為12a。位于甘孜-玉樹斷裂帶最西段的當(dāng)江段作為地震空區(qū),最新一次地震事件為1738年當(dāng)江M7?2 地震[30],離逝時間為284a,接近斷裂的強震復(fù)發(fā)周期。聞學(xué)澤等[22]對斷錯地貌進(jìn)行地形測圖、斷距測量等定量研究,結(jié)合熱釋光樣品結(jié)果,得到了斷裂5萬年以來的左旋滑動速率為(12±2)mm/a。同時,Yu等[30]對甘孜-玉樹斷裂當(dāng)江段進(jìn)行構(gòu)造地貌調(diào)查、探槽古地震和年代學(xué)樣品測試等方法,得到當(dāng)江段滑動速率為(6.3±1.9)mm/a。分析認(rèn)為,自1738年當(dāng)江M7?2 地震之后,當(dāng)江段產(chǎn)生的水平位移量理論值為1.79~3.41m(表3)。分別利用震級-水平位移量的經(jīng)驗關(guān)系式:M=7.13+0.68log(D)、M=7.04+0.89log(D),估算震級為M7.2~7.5級(表3)。綜合分析并構(gòu)造類比認(rèn)為,當(dāng)江段可能是鮮水河-甘孜-玉樹斷裂帶的強震空區(qū),其大震離逝時間已接近或甚至超過強震復(fù)發(fā)間隔,其短期內(nèi)發(fā)生7級以上強震的危險性不容低估。
圖6 甘孜-玉樹斷裂帶歷史地震地表破裂分布圖[30-31]Fig.6 Sur face rupture distribution map of historical earthquakes in Ganzi-Yushu Fault Zone[30-31]
(1)2022年瀘定MS6.8地震地表的同震地表變形特征主要表現(xiàn)為裂縫及滑坡。裂縫北起燕子溝,經(jīng)磨西鎮(zhèn)、二臺子、猛虎崗,南至幸福村,整體呈NW 向展布,全長約35km?;率軘嗔褞Э刂坪痛蠖珊铀档目臻g展布影響較大。西部滑坡區(qū)以燕子溝至愛國村一帶最為嚴(yán)重,滑坡體規(guī)模大,范圍廣,整體呈NNW 向,平行于鮮水河斷裂磨西段延伸約40km;東部滑坡區(qū)主要分布于得妥鎮(zhèn)以北約5km 處至瀘定縣城西南端,滑坡體沿大渡河兩岸呈NE 向展布,延伸約10km。
(2)巴顏喀拉塊體及鄰近地區(qū)短期內(nèi)仍可能是強震活動的高風(fēng)險地區(qū),塊體南邊界帶鮮水河斷裂西北端當(dāng)江段的強震危險性尤需高度重視。
致謝:四川省地震局、甘孜藏族自治州應(yīng)急管理局、雅安市應(yīng)急管理局、磨西鎮(zhèn)燕子溝地震臺、防災(zāi)科技學(xué)院、北京防災(zāi)科技有限公司、相關(guān)防疫部門、搶險救災(zāi)部門給予了本次野外科考與調(diào)查工作極大的支持與幫助。文中部分基礎(chǔ)數(shù)據(jù)參考了四川省地震局、中國地震局地質(zhì)研究所、中國地震局地球物理研究所、中國地震局工程力學(xué)研究所、中國地震局地震預(yù)測所、中國地震臺網(wǎng)中心以及中國地震局震后科考的相關(guān)研究成果。在此謹(jǐn)對上述給予幫助和支持的單位與個人表示衷心的感謝并致以崇高的敬意!