云南漾濞6.4級地震震源區(qū)及周邊的重力均衡特征

石 磊 李永華 張瑞青

1)中國地震局地球物理研究所， 北京 100081 2)中國地震局震源物理重點(diǎn)實(shí)驗(yàn)室， 北京 100081

0 引言

2021年5月21日云南省大理市漾濞縣發(fā)生MS6.4 地震， 震中位置為(99.87°E， 25.67°N)。該地震造成35人傷亡， 受災(zāi)面積近6500km2。震源機(jī)制研究結(jié)果表明， 此次地震為右旋走滑型， 具有正斷分量。該地震的發(fā)震位置位于紅河斷裂帶北段延伸方向的維西-喬后-巍山斷裂附近(圖 1)， 但現(xiàn)有資料顯示， 此次地震未發(fā)生在已知的維西-喬后-巍山斷裂上， 至于是否發(fā)生在隱伏斷裂上， 還有待深入研究(雷興林等， 2021)。因此， 了解震源區(qū)及周邊地殼結(jié)構(gòu)與物性特征， 對進(jìn)一步認(rèn)識該地震的發(fā)震構(gòu)造具有重要意義。

圖 1 研究區(qū)主要構(gòu)造與地形Fig. 1 Tectonic features and topography in the study area， colors show the elevations.LCJF 瀾滄江斷裂帶； WX-QH-WS 維西-喬后-巍山斷裂帶； RRF 紅河斷裂帶； XJH-LJF 小金河-麗江斷裂帶； JHF 菁河斷裂帶； ANHF 安寧河斷裂帶； ZMH-XJF 則木河-小江斷裂帶。紅色五角星為本次漾濞 MS6.4 地震的震中位置， 藍(lán)色實(shí)線為斷裂(Deng et al.， 2002)， 暗紅色圓點(diǎn)為實(shí)際測量剖面的位置

青藏高原東南緣由于受到印度板塊對青藏高原推擠作用的影響， 陸內(nèi)變形強(qiáng)烈， 發(fā)育了眾多深大斷裂， 地震活動(dòng)頻繁。重力是研究殼內(nèi)結(jié)構(gòu)與深部孕震環(huán)境的學(xué)科之一， 前人已在研究區(qū)開展了大量研究， 對重力場特征(鐘鍇等， 2005； 姜文亮等， 2011)、 莫霍面深度(黃建平等， 2006； Shietal.， 2015； Zhaoetal.， 2018)、 密度結(jié)構(gòu)(Zhangetal.， 2013； 石磊等， 2015； 楊文采等， 2015)和均衡異常(王謙身等， 2009； 張永謙等， 2010)等都進(jìn)行了深入分析。重力均衡可較好地反映地殼結(jié)構(gòu)及與構(gòu)造應(yīng)力的關(guān)系， 與地震活動(dòng)密切相關(guān)(方盛明等， 1997)。但現(xiàn)有研究對于重力均衡異常與地震分布的關(guān)系還存在爭議。如王謙身等(2009)研究結(jié)果顯示， 汶川地震位于均衡異常極值區(qū)， 其認(rèn)為艾里均衡理論莫霍面深度與實(shí)際莫霍面深度相差越大的區(qū)域越有可能發(fā)生地震。張永謙等(2010)的研究結(jié)果表明， 川西藏東地區(qū)的強(qiáng)震多發(fā)生在均衡異常梯度帶上， 頻率明顯高于極值區(qū)。當(dāng)前計(jì)算艾里均衡地殼厚度時(shí)常將殼幔密度差取為常數(shù)， 前人根據(jù)不同資料和方法選取的殼幔密度差并不相同(Wangetal.， 2010； 張恩會(huì)， 2015)， 且實(shí)際地殼密度分布也存在差異(樓海等， 2005)， 考慮密度差橫向變化的重力均衡異常能更好地反映地殼結(jié)構(gòu)。

本文基于區(qū)域布格重力異常與前人采用接收函數(shù)H-κ疊加掃描法估計(jì)得到的莫霍面深度， 獲得漾濞6.4級地震震源區(qū)及周邊的橫向變化殼幔密度差。根據(jù)艾里均衡理論計(jì)算出由地形起伏對應(yīng)的莫霍面深度， 通過對比莫霍面的理論均衡深度與接收函數(shù)得到的實(shí)際深度的差異分析整個(gè)研究區(qū)的均衡異常， 并探討均衡異常與中強(qiáng)地震分布的關(guān)系， 進(jìn)一步厘定漾濞地震震源區(qū)的均衡狀態(tài)， 以期為后續(xù)研究提供一些重力學(xué)參考。

1 數(shù)據(jù)與方法

1.1 重力與地形數(shù)據(jù)

布格異常數(shù)據(jù)來自世界重力圖網(wǎng)格數(shù)據(jù)庫(WGM2012)， 網(wǎng)格間距為10km×10km(圖 2)。該數(shù)據(jù)庫是由國際重力測量局依據(jù)高分辨率的地球重力模型EGM2008(Pavlisetal.， 2012)和高程模型ETOPO1(Amanteetal.， 2009)采用球諧方法計(jì)算得到的(Balminoetal.， 2012)， 計(jì)算中考慮了實(shí)際的地球模型和大多數(shù)地表物質(zhì)的質(zhì)量分布(空氣、 陸地、 海洋、 內(nèi)陸海洋、 湖泊、 冰蓋和冰架)。大部分地區(qū)的分辨率為5′， 部分插值地區(qū)(如天山、 青藏高原東南緣等)的分辨率為15′(Steffenetal.， 2011)。Li等(2009)對青藏高原及周邊地區(qū)開展了研究， 發(fā)現(xiàn)由EGM2008自由空氣異常計(jì)算得到的布格異常與地面數(shù)據(jù)較接近， 標(biāo)準(zhǔn)差為18mGal。Kiamehr等(2009)分別使用EGM2008和實(shí)測重力數(shù)據(jù)進(jìn)行反演， 獲得的莫霍面形態(tài)和深度只有較小的差異。

圖 2 研究區(qū)的布格重力異常Fig. 2 The Bouguer gravity anomaly in the study area.

研究區(qū)的高程數(shù)據(jù)來自計(jì)算布格異常時(shí)采用的全球地形模型ETOPO1， 分辨率為1′(圖 1)。ETOPO1模型由美國國家地球物理數(shù)據(jù)中心發(fā)布， 基于全球和地區(qū)的不同數(shù)據(jù)集經(jīng)過系統(tǒng)處理和編輯生成(Amanteetal.， 2009)。中國地震局于2011年在研究區(qū)對2條剖面進(jìn)行了重力和地形探測工作(石磊等， 2015)， 實(shí)際剖面位置見圖 1。將實(shí)際測量的布格重力和地形分別與WGM2012和ETOPO1模型提取的相同點(diǎn)數(shù)據(jù)進(jìn)行對比(圖 3)， 剖面AA′的重力和地形標(biāo)準(zhǔn)差分別為4.89mGal和17.86m， 剖面BB′的重力和地形標(biāo)準(zhǔn)差分別為8.42mGal和28.74m。對比可見， 研究區(qū)內(nèi)WGM2012和ETOPO1模型數(shù)據(jù)都與實(shí)際測量結(jié)果有較好的對應(yīng)， 可用于后續(xù)的均衡研究。

圖 3 剖面實(shí)際測量結(jié)果與WGM2012和ETOPO1模型數(shù)據(jù)的對比Fig. 3 The comparison of measurement gravity anomaly(a， c)and topography(b， d) with the models of WGM2012 and ETOPO1.紫色實(shí)線為實(shí)際測量重力數(shù)據(jù)得到的布格異常， 藍(lán)色虛線為WGM2012模型提取的相同點(diǎn)的布格異常， 黃色實(shí)線為實(shí)際測量的測點(diǎn)地形， 紅色虛線為ETOPO1模型提取的相同點(diǎn)的高程。剖面位置見圖1

1.2 莫霍面深度

遠(yuǎn)震接收函數(shù)H-κ疊加掃描法(Zhuetal.， 2000)是獲得殼幔速度界面最有效的方法之一， 并已得到廣泛應(yīng)用(Tianetal.， 2013； Wangetal.， 2016)。本文收集了基于以往川西流動(dòng)臺陣和 “中國地震科學(xué)臺陣探測”項(xiàng)目(Liuetal.， 2014； Wangetal.， 2017)在南北地震帶南段布設(shè)的358個(gè)流動(dòng)臺站記錄的遠(yuǎn)震數(shù)據(jù)利用H-κ掃描法獲得的地殼厚度(Zhengetal.， 2019)。上述流動(dòng)臺陣的平均間距為30～40km。將接收函數(shù)估算的地殼厚度減去臺站的實(shí)際高程， 即可獲得研究區(qū)的莫霍面埋深分布(圖 4)。由圖 4 可知， 研究區(qū)莫霍面埋深的變化范圍為28.8～63.9km， 且總體上呈現(xiàn)出由NW向SE減薄的趨勢。

圖 4 基于接收函數(shù)估計(jì)得到的莫霍面深度(Zheng et al.， 2019)Fig. 4 The Moho depth estimated from receiver functions(Zheng et al.， 2019).黑色三角形為研究區(qū)流動(dòng)臺站分布位置

1.3 估計(jì)殼幔密度差

假設(shè)地殼和上地幔的平均密度分別為ρc和ρm， 則莫霍面起伏Δh在地表任一點(diǎn)(x，y)引起的重力異常為

(1)

其中，G為引力常數(shù)，ρm-c=ρm-ρc為殼幔密度差。其在頻率域的形式為(Parker， 1973)

(2)

式中，k為波數(shù)，h0為莫霍面參考深度。

若重力異常Δg、 莫霍面參考深度h0和實(shí)際深度均已知， 界面起伏Δh=H-h0， 則根據(jù)式(2)可反演得到殼幔密度差ρm-c(ξ，η)。

1.4 計(jì)算參考莫霍面深度

假設(shè)布格重力異常由不同波長的異常疊加組成， 則異常的功率譜密度函數(shù)可看作是不同深度場源的功率譜密度函數(shù)的總和(Guoetal.， 2013)。重力異常的功率譜密度函數(shù)可利用不同深度的格林等效層構(gòu)建(Pwalowski， 1994)， 任意深度h處的格林等效層功率譜密度函數(shù)為(Naidu， 1968； Dampney， 1969)

P=se-2kh

(3)

式中，s為與格林等效層強(qiáng)度有關(guān)的常數(shù)，k為波數(shù)。

若采用直線近似擬合重力異常徑向平均對數(shù)功率譜的某個(gè)分段， 則利用式(4)可估計(jì)等效層的近似深度(郭良輝等， 2012)：

(4)

我們對研究區(qū)的布格異常進(jìn)行傅里葉變換， 計(jì)算徑向平均對數(shù)功率譜， 根據(jù)功率譜的形狀特征和波數(shù)范圍確定參考莫霍面對應(yīng)的格林等效層數(shù)及近似深度。

2 艾里均衡理論莫霍面

重力均衡理論認(rèn)為， 在地下某個(gè)補(bǔ)償深度單位截面上承受的質(zhì)量應(yīng)該是相等的(Dutton， 1925)。如果地球表面存在過剩的負(fù)載， 那么在補(bǔ)償深度之上一定存在等效的質(zhì)量虧損， 反之則會(huì)存在質(zhì)量盈余。根據(jù)艾里均衡理論， 在陸地上， 山根厚度t=ρcT/(ρm-ρc)=ρcT/ρm-c， 理論均衡莫霍面深度D=h0+t。其中，t為山根厚度，T為地表高程，ρc、ρm和ρm-c分別為上節(jié)提到的地殼平均密度、 上地幔頂部密度和殼幔密度差，h0為參考莫霍面深度。

本文采用計(jì)算重力異常徑向平均對數(shù)功率譜的方法估計(jì)參考深度(Spectoretal.， 1970)。研究區(qū)的布格重力異常功率譜呈指數(shù)遞減(圖 5)。通過分析我們認(rèn)為， <0.01km-1波數(shù)值范圍的異常是由淺層場源引起的。在一定波數(shù)范圍內(nèi)， 徑向平均對數(shù)功率譜回歸線的斜率(直線近似擬合)代表了等效層的深度。在0.01～0.02km-1的波數(shù)范圍內(nèi)， 擬合的回歸線表示研究區(qū)的莫霍基準(zhǔn)面深度h0為38km， 與前人在青藏高原及鄰區(qū)的研究結(jié)果一致(Jimenez-Muntetal.， 2008； Wangetal.， 2010)。

圖 5 參考莫霍面深度估計(jì)Fig. 5 The estimation of reference Moho.紅色直線為功率譜的擬合回歸線

基于收集到的前人(Zhengetal.， 2019)使用接收函數(shù)H-κ疊加法得到的莫霍面深度H和本文采用重力異常徑向平均對數(shù)功率譜法估計(jì)的莫霍面參考深度h0(38km)， 根據(jù)Δh=H-h0可計(jì)算出研究區(qū)莫霍界面起伏Δh。已知布格重力異常、 波數(shù)、 莫霍界面參考深度和起伏， 我們根據(jù)式(2)反演得到了研究區(qū)橫向變化的殼幔密度差ρm-c(ξ，η)。上地幔頂部的密度與地殼平均密度相差0.3～0.55g/cm3， 總體上自北向南、 自西向東逐漸減小(圖 6)。

圖 6 研究區(qū)的殼幔密度差Fig. 6 The crust-mantle density contrast in the research area.

漾濞地震的震源區(qū)位于密度過渡帶， 過渡帶變化的方向與維西-喬后-巍山斷裂帶的走向近垂直。四川盆地表現(xiàn)為相對低的密度差， 與前人的研究結(jié)果基本一致(王謙身等， 2013； 張恩會(huì)等， 2018)， 反映四川盆地為穩(wěn)定的構(gòu)造塊體(Huangetal.， 2002； 李永華等， 2014)。青藏高原東南緣的布格重力異常為較低的負(fù)值， 使得反演結(jié)果為高密度差(0.5～0.55g/cm3)， 推測與前人層析成像研究揭示的中下地殼廣泛存在低速層有關(guān)(Wangetal.， 2003； Lietal.， 2016)。深地震測深剖面(Tsengetal.， 2009)采用較小的殼幔密度差(0.45g/cm3)， 研究發(fā)現(xiàn)松潘-甘孜地塊的艾里均衡理論莫霍面深度比實(shí)際的莫霍面更深， 與本文的結(jié)果相反。這是由于山根的厚度與殼幔密度差成反比， 當(dāng)密度差取值較小時(shí)山根厚度偏厚， 相反地， 密度差取值較大時(shí)山根厚度則偏淺。

基于莫霍基準(zhǔn)面深度h0和殼幔密度差ρm-c(ξ，η)， 進(jìn)一步得到研究區(qū)理論均衡莫霍面深度(圖 7)。均衡理論莫霍面在青藏高原東南緣的平均深度約為55km， 向SE逐漸減??； 滇中次級塊體的深度約45km； 四川盆地最淺， 為40km。攀西地區(qū)表現(xiàn)為局部高值， 是由于該區(qū)殼幔密度差較小， 而地形與鄰區(qū)相近， 理論莫霍面深度與殼幔密度差成反比。人工寬角地震與天然地震層析成像研究(吳建平等， 2009； Baoetal.， 2015； 徐濤等， 2015)結(jié)果都顯示攀西地區(qū)中下地殼存在明顯的高速異常， 高速度一般對應(yīng)高密度， 與該區(qū)殼幔密度差較小一致。

圖 7 研究區(qū)的艾里均衡理論莫霍面深度Fig. 7 The Airy Moho in the research area.

3 重力均衡分布特征

3.1 區(qū)域均衡與地震活動(dòng)

研究區(qū)的重力均衡異常I=D-H(圖 8)， 其中，D為艾里均衡理論莫霍面深度，H為接收函數(shù)H-κ疊加法得到的深度。異常較小或接近零時(shí)表明相對 “穩(wěn)定”， 異常絕對值較大則說明該區(qū)處于不均衡的情況， 構(gòu)造應(yīng)力將調(diào)整深部結(jié)構(gòu)使其達(dá)到穩(wěn)定狀態(tài)(王謙身等， 2009)。我們將研究區(qū)-2km≤均衡異?！?km的地區(qū)劃分為基本均衡的狀態(tài)。青藏高原東南緣為負(fù)值區(qū)， 即理論艾里莫霍面深度小于實(shí)際深度， 表明山根厚度對地表地形起伏的補(bǔ)償過剩， 該地區(qū)應(yīng)該 “上升”(劉云龍等， 1977)， 使理論深度D增加以達(dá)到均衡。紅河斷裂帶西南和東側(cè)表現(xiàn)為正異常， 即理論艾里莫霍面深度大于實(shí)際深度， 山根厚度對地表地形起伏的補(bǔ)償不足， 這些地區(qū)應(yīng)該 “下降”， 使理論深度D減小， 從而達(dá)到均衡條件。

為進(jìn)一步分析均衡異常與中強(qiáng)地震分布的關(guān)系， 將研究區(qū)1900年以來6.5級及以上地震和此次漾濞地震的震源投影到重力均衡異常圖中(圖 8)。由圖8可見， 崇明8.0級古地震發(fā)生在均衡異常為4km的梯度帶上。研究區(qū)7.0～7.9級地震共29個(gè)， 其中16個(gè)(55%)地震的震中位于均衡異常正、 負(fù)過渡帶， 5個(gè)(17%)地震震中位于異常高值區(qū)伴隨的梯度帶上。42個(gè)6.5～6.9級地震中， 處于均衡異常正、 負(fù)過渡帶的有21個(gè)(50%)， 處于高值區(qū)伴隨梯度帶的有13個(gè)(31%)， 分布在均衡異常極值區(qū)的有9個(gè)(21%)。前人在川滇地區(qū)的研究結(jié)果也顯示， 中強(qiáng)地震多發(fā)生在均衡異常梯度帶附近(張永謙等， 2010； 張恩會(huì)； 2015)。均衡異常梯度帶為異常變化較快的地區(qū)， 推測這些區(qū)域可能處于更不穩(wěn)定的狀態(tài)， 有利于應(yīng)變能的積累， 為使地殼達(dá)到均衡， 易于從這些地方開始調(diào)整， 從而導(dǎo)致地震的發(fā)生。

圖 8 研究區(qū)的重力均衡異常Fig. 8 The gravity isostatic anomaly in the research area.圖中不同大小的圓點(diǎn)代表6.5級及以上震級地震的震中， 紅色三角形為漾濞地震震中， 黑色圓點(diǎn)為余震震中， 棕色實(shí)線為穿過震中和余震的剖面位置

3.2 漾濞地震地殼深部結(jié)構(gòu)

研究中我們還分別對穿過漾濞地震主震和余震的2條剖面的高程、 布格重力異常和艾里理論與實(shí)際莫霍面深度差異進(jìn)行了分析(圖 9)。剖面AA′主震震中西側(cè)的地形起伏變化不大， 東側(cè)維西-喬后-巍山斷裂帶附近為整條剖面高程最大的地區(qū)， 與震中存在600m的高差。布格異常全部為負(fù)值， 變化范圍為-256～-180mGal， 震中位于局部低異常區(qū)。剖面BB′余震西側(cè)地形存在一定的起伏變化， 東側(cè)高程最大值同樣出現(xiàn)在已知的維西-喬后-巍山斷裂帶， 存在近1000m的變化。布格異常的變化范圍與AA′類似， 余震同樣分布在局部低值區(qū)。

圖 9 2條剖面的地形、 布格異常和D、 H差異曲線Fig. 9 Topography(a， b)， Bouguer gravity anomaly(c， d)and the differences between D and H(e， f)for two profiles.

由圖8 可以看出， 漾濞地震同樣位于重力均衡異常梯級帶上， 且梯級帶的方向垂直于維西-喬后-巍山斷裂帶的走向。剖面AA′上D與H的最大差異出現(xiàn)在瀾滄江斷裂帶以西(10km)， 跨過瀾滄江斷裂帶后差異逐漸減小至2km內(nèi)， 達(dá)到近似穩(wěn)定的狀態(tài)。漾濞地震主震震中D與H存在約4km的差異， 至維西-喬后-巍山斷裂帶差異增加到6km。剖面BB′上瀾滄江斷裂帶以西D與H的差異為11km， 瀾滄江斷裂帶至維西-喬后-巍山斷裂帶一段D與H的差異在4～6km范圍內(nèi)波動(dòng)變化， 但沒有出現(xiàn)快速增大或減小的現(xiàn)象。對比分析可知， 漾濞地震同樣發(fā)生在重力均衡異常變化較快、 地殼結(jié)構(gòu)由不均衡—均衡—不均衡狀態(tài)的過渡區(qū)。

4 結(jié)論

本文基于WGM2012布格重力異常、 ETOPO1高程和前人采用接收函數(shù)H-κ疊加法得到的莫霍面深度獲得了漾濞6.4級地震震源區(qū)和鄰區(qū)重力均衡異常。不同于以往計(jì)算山根厚度時(shí)將殼幔密度差取為常數(shù)， 我們考慮了密度差的橫向變化， 得到了與實(shí)際地殼密度分布更為相符的研究區(qū)重力均衡異常。

研究結(jié)果表明， 地殼的平均密度與上地幔頂部密度的差異為0.3～0.55g/cm3， 漾濞地震的震源區(qū)位于密度差過渡帶上。四川盆地表現(xiàn)出相對較小的密度差， 反映了四川盆地為穩(wěn)定的構(gòu)造塊體(Huangetal.， 2002； 李永華等， 2014)。青藏高原東南緣的密度差較大(0.5～0.55g/cm3)， 與中下地殼廣泛存在的低速層有關(guān)(Wangetal.， 2003； Lietal.， 2016； 鄭晨等， 2016)。

研究區(qū)的理論均衡莫霍面深度為35～60km， 總體上呈現(xiàn)由NW向SE減薄的趨勢。攀西地區(qū)表現(xiàn)為局部高值， 是由于該區(qū)殼幔密度差較小， 而地形與鄰區(qū)相近， 理論莫霍面深度與殼幔密度差成反比。人工寬角地震與天然地震層析成像研究(吳建平等， 2009； Baoetal.， 2015； 徐濤等， 2015)結(jié)果都顯示攀西地區(qū)中下地殼存在明顯的高速異常， 高速度一般對應(yīng)高密度， 與該區(qū)殼幔密度差較小一致。

我們將艾里均衡理論莫霍面深度與接收函數(shù)H-κ疊加掃描法得到莫霍面深度求差， 得到研究區(qū)的均衡異常。崇明8.0級古地震發(fā)生在均衡異常為4km的梯度帶上。研究區(qū)內(nèi)發(fā)生的7.0～7.9級地震共29個(gè)， 其中16個(gè)(55%)地震的震中位于均衡異常正、 負(fù)過渡帶上， 5個(gè)(17%)地震的震中位于異常高值區(qū)伴隨的梯度帶上。在42個(gè)6.5～6.9級地震中， 處于均衡異常正、 負(fù)過渡帶和高值區(qū)伴隨梯度帶的分別有21個(gè)(50%)和13個(gè)(31%)， 分布在均衡異常極值區(qū)的有9個(gè)(21%)。

漾濞地震震源區(qū)處于均衡異常為4km的梯度帶上， 且梯度帶與維西-喬后-巍山斷裂帶的走向基本一致， 地殼深部結(jié)構(gòu)處于不均衡—均衡—不均衡狀態(tài)的過渡區(qū)。我們推測重力均衡異常梯度帶可能處于不穩(wěn)定狀態(tài)， 從而有利于應(yīng)變能的積累， 為使地殼達(dá)到均衡， 多從這些地方開始調(diào)整， 利于地震的發(fā)生。

致謝中國地震臺網(wǎng)中心提供了漾濞地震余震序列資料； 中國地震局地球物理研究所樓海研究員和評審專家為本文提出了寶貴的意見與建議。 在此一并表示感謝!