盛書中, 萬永革, 胡曉輝, 宮猛, 張?zhí)K祥

1 東華理工大學地球物理與測控技術學院， 南昌 330013 2 防災科技學院， 河北三河 065201 3 中國科學技術大學， 合肥 230026

0 引言

地殼應力場的積累與增強會引起斷層錯動，斷層錯動又會引起地殼應力場的釋放，斷層運動與地殼應力場間不斷地調(diào)整，因此，詳細研究震源區(qū)應力場變化是認識地震孕育規(guī)律的一個關鍵因素.Yamakawa(1971)討論了地震會對震源區(qū)應力場產(chǎn)生擾動，從而引起震源區(qū)應力偏轉.許忠淮(1985)研究唐山地震震源區(qū)應力場時，發(fā)現(xiàn)震后最大主應力軸偏轉了約30°，后續(xù)大量研究均表明震后震源區(qū)存在應力偏轉現(xiàn)象(Michael，1987；Hauksson，1994；刁桂苓等，1994，2005；Zhao et al.，1997；Hardebeck and Hauksson，2001；Ickrath et al.，2014；Sheng and Meng，2020).Ickrath等(2014)研究了1999年土耳其MW7.4伊茲米特地震震源區(qū)應力場，研究結果表明主震改變了震源區(qū)應力場性質(zhì)，震前震源區(qū)應力場為走滑型應力場，主震后震源區(qū)應力場變?yōu)檎龜嘈蛻?，主震后?jīng)過2個月的調(diào)整，震源區(qū)應力場又恢復到震前走滑型應力場.Hardebeck和Okada(2018)較系統(tǒng)地總結了國外地震引起的應力場偏轉現(xiàn)象方面的研究，發(fā)現(xiàn)主震震級小于～6.2級時引起的應力偏轉現(xiàn)象報道較少，認為引起顯著應力偏轉現(xiàn)象的臨界震級約為6.2級，但也有震級較小地震引起顯著應力場變化的震例，如美國加利福尼亞州圣哈辛托斷裂上的一次5.5級地震則引起主應力軸傾角的變化(Martínez-Garzón et al., 2016)，他們估計該臨界震級受無偏轉震例研究報道較少的影響，因此，仍需要開展大量的研究工作去進一步確定應力偏轉的震級閾值.震源區(qū)震后應力偏轉現(xiàn)象為我們揭示了地震對應力場的影響，同時也提供了估計震源處應力量值的新途徑(Yin and Rogers，1995；Hardebeck and Hauksson，2001；萬永革等，2006；Wesson and Boyd，2007；盛書中和萬永革，2011；萬永革等，2012).可見，同震或震后應力場變化可以為我們揭示地震過程中震源區(qū)應力演化特征以及估計震源區(qū)應力量值，所以，震源區(qū)應力場變化具有重要的研究意義.

云南盈江地區(qū)中強震頻發(fā)，備受研究者關注，研究者圍繞歷次中強震序列的發(fā)震構造、震源機制解和地震重定位等開展了一系列研究工作.徐彥等(2012)基于2008年盈江地震序列中強震震源機制解測定與地震序列重定位結果，對該序列的發(fā)震構造及地殼流變特征進行了分析.2011年3月10日盈江發(fā)生MS5.9地震，震后研究人員就該地震序列震源機制解、重定位和應力降等開展了一系列研究工作(房立華等，2011；趙小艷等，2012；鄧菲等，2016).2014年5月盈江陸續(xù)發(fā)生2次中強震，震后研究人員針對該雙震開展了詳細的研究工作，如許力生等(2014)利用云南省地震臺網(wǎng)記錄的波形資料研究了該雙震的破裂歷史；為了研究其發(fā)震構造，研究人員測定了震源機制解，并對地震序列開展了重定位工作(趙旭等，2014；楊婷等，2016)；此外，對本次地震序列的剪切波分裂也開展了詳細的研究(李金等，2015；孫楠等，2017；吳朋等，2018)，研究結果顯示主震前快剪切波優(yōu)勢偏振方向與區(qū)域主壓應力方向一致，快剪切波偏振方向在地震序列過程中存在變化.Xu等(2015)對2008—2011年間盈江地區(qū)發(fā)生的5個地震序列進行了重定位研究，結果表明2011年3月10日的地震序列發(fā)生在以71°向南西傾斜的雙層斷層上，其余震序列分布和區(qū)域已知主要斷裂走向相共軛.Xu等(2020)對云南地區(qū)震源機制解開展了詳細的研究工作，并給出云南地區(qū)自2000年1月至2014年12月的統(tǒng)一震源機制解目錄，為深入研究該地區(qū)應力場提供了基礎.

有關云南地區(qū)地殼應力場已經(jīng)開展了大量研究工作，大多數(shù)研究工作均涉及盈江地區(qū)(錢曉東等，2011；Zhao et al., 2013；Xu et al., 2016；孫業(yè)君等，2017)，但缺少針對盈江地區(qū)中強震對該地區(qū)地殼應力場影響的研究.因此，本研究目標為利用2008年以來盈江地區(qū)地震震源機制解，詳細研究該地區(qū)發(fā)生的中強震是否對區(qū)域地殼應力場產(chǎn)生擾動，即中強震對應力場的影響，確定盈江地區(qū)地殼應力場及其是否存在時空變化，為盈江地區(qū)地震孕育環(huán)境、地質(zhì)構造活動分析和中強震對應力場影響等研究提供參考，以及為應力場時空變化研究中是否需要考慮中強震對應力場影響提供參考.

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

本研究的研究區(qū)域為北緯24.5°—25.4°、東經(jīng)97.6°—98.2°，時間范圍為2008年3月21日—2014年11月19日，從Xu等(2020)的研究中收集到震源機制解262個，震源機制解空間分布情況見圖1，所用地震的矩震級范圍為3.03～5.82，震源深度主要分布在20 km以內(nèi).震源機制解類型為15個正斷層型，10個正斷走滑型，189個走滑型，6個逆沖兼走滑型，6個逆沖型，36個過渡型(Zoback，1992)，可見，盈江地區(qū)地震以走滑型地震為主.

圖1 各時段震源機制解空間分布圖(a) 、(b)、 (c)、 (d)和(e)分別對應表2中的5個時段；(f)為全部震源機制解空間分布圖；紅點示意中強震震中，圖中數(shù)字對應表1中的中強震序號.Fig.1 The focal mechanism solutions spatial distribution of each time period(a) to (e) are consistent with the 5 periods in Table 2； (f) shows the spatial distribution of all focal mechanisms solutions, red dots represent the moderate earthquakes epicenter, Numbers in graphs correspond to the serial numbers of moderate earthquakes in Table 1.

本研究中應力場反演我們采用MSATSI軟件(Martínez-Garzón et al., 2014)，該軟件包是Hardebeck和Michael(2006)給出的區(qū)域應力場阻尼反演程序SATSI(Spatial And Temporal Stress Inversion)的Matlab版，增加了相應的繪圖功能.MSATSI軟件在應力場反演和結果成圖方面較為便捷，因此，被廣泛應用于區(qū)域應力場反演工作(Ickrath et al.，2014；王曉山等，2015；Xu et al.，2016；Sheng and Meng，2020).考慮應力場的時空連續(xù)性，Hardebeck和Michael(2006)基于Michael(1987)線性應力場反演方法，增加相鄰網(wǎng)格點間應力場連續(xù)性約束，在綜合考慮應力場反演殘差和相鄰網(wǎng)格點間差異的基礎上，給出了區(qū)域應力場阻尼反演方法.該方法考慮了應力場的時空連續(xù)性，因此，可以有效地抑制由于單獨反演各個網(wǎng)格點而產(chǎn)生相鄰網(wǎng)格點間應力場出現(xiàn)突然變化的假象.

2 應力場反演結果及分析

使用上述數(shù)據(jù)和方法，我們研究中強震對盈江地區(qū)應力場的影響.首先，為了研究5次中強震對研究區(qū)應力場的擾動，中強震參數(shù)見表1，根據(jù)發(fā)震時間將數(shù)據(jù)分為5段，反演各時段的應力場，研究結果見表2和圖3.其次，根據(jù)震源機制解空間分布的叢集性，我們將研究區(qū)分為南北兩部分，分別反演應力場，以考察南北部應力場是否存在差異.再次，為了考察研究時段內(nèi)，盈江地區(qū)地殼應力場是否存在時空演化，我們在全部數(shù)據(jù)中隨機抽取震源機制解反演應力場.最后，鑒于前面的研究表明盈江地區(qū)在研究時段內(nèi)存在穩(wěn)定的應力場，我們將整個研究時段的全部震源機制解作為整體進行了應力場反演，獲得盈江地區(qū)該時段內(nèi)的總體應力場，研究結果見表3和圖4.

表1 中強震參數(shù)表Table 1 Parameters of the moderate earthquakes

表2 各時段應力場反演結果表Table 2 The stress field inversion result of each time period

2.1 各時段應力場反演結果及分析

先前研究認為引起應力場偏轉的地震震級下限為～6.2級(Hardebeck and Okada，2018)，但也有較小地震引起應力場偏轉的震例(Martínez-Garzón et al., 2016)，本次研究的中強震震級小于該震級下限，這意味著這些地震的發(fā)生可能對應力場沒有產(chǎn)生顯著的擾動，即盈江地區(qū)的應力場在時間上仍是連續(xù)的.因此，考慮到時間上應力場可能是連續(xù)的，我們使用阻尼應力場反演了該地區(qū)應力場；同時，為了捕獲中強震可能引起的應力場變化，我們分時間段單獨反演各時段應力場，即進行了阻尼和無阻尼的分時段應力場反演.阻尼應力場反演時，根據(jù)數(shù)據(jù)誤差和模型長度的折中曲線確定阻尼系數(shù)大小為1.6(見圖2).由表2可見，前3個時段內(nèi)所用震源機制解的旋轉角均方根值基本滿足要求，后兩個時段的旋轉角均方根值相對較小，總體滿足震源機制解多樣性要求.除第三個時段的單獨反演殘差平均值為36.1°外，其余各時段應力場反演結果的殘差平均值均小于35°，反映各時段內(nèi)所用的震源機制解滿足應力場均勻性假設.應力場反演的兩個基本假設均得到滿足，反映了應力場反演結果是可靠的.

圖2 模型長度與數(shù)據(jù)擬合誤差之間的折中曲線Fig.2 The trade-off curve of model length and data misfit

盡管阻尼應力場反演給出的置信區(qū)間相對較小，但總體上兩種應力場反演結果是一致的(見表2和圖3)，反映了研究結果的穩(wěn)定性.

各時段應力場反演結果均表明：盈江地區(qū)的應力場為走滑型，最大主應力軸呈NNE-SSW向，最小主應力軸呈SEE-NWW向，最大和最小主應力軸傾角近水平，中間主應力軸傾角近直立，該結果與先前的研究結果一致(錢曉東等，2011；Zhao et al.，2013；Xu et al.，2016；孫業(yè)君等，2017；Jin et al.，2019).

總體上兩種反演獲得的應力場結果均表明：中強震對應力場產(chǎn)生了擾動，但主應力軸方位變化不顯著(見表2和圖3).阻尼應力場反演結果表明：各時段主應力軸走向和傾角最優(yōu)值變化較小，最大主應力軸方位角的最大和最小值分別為215°和206°(26°的反方向).主應力軸的95%置信區(qū)間均非常集中(見圖3)，各時段主應力軸走向和傾角的最優(yōu)值差異較小且它們的95%置信區(qū)間有重疊，可見，阻尼應力場反演結果中最大主應力軸存在順時針偏轉現(xiàn)象但并不顯著.無阻尼應力場反演結果表明：各時段應力場發(fā)生了相對較大的變化，最大主應力軸走向的最大旋轉角達到19°，發(fā)生在第3個和第4個地震之后的應力場間，同時，這兩個時段應力場結果的置信區(qū)間較大且有相互重疊部分，因此，雖然這兩個時段應力場存在較大差異，但考慮置信區(qū)間后，則它們的差異不顯著.綜上可見，兩種反演所得應力場結果均顯示各時段應力場存在差異，但考慮到各時段主應力軸的95%置信區(qū)間存在重疊區(qū)域 (見圖3)，則中強震引起的應力場變化并不顯著.因此，這5次中強震對應力場存在影響但并不顯著.

圖3 各時段應力場反演結果圖圖中黑色加號示意應力場反演的最優(yōu)值，紅色、綠色和藍色點示意2000次抽樣給出的最大、中間和最小主應力軸的95%置信區(qū)間； (a)為阻尼應力場反演結果， (b) 為無阻尼應力場反演結果；數(shù)字表示各時段與表2中序號一致.Fig.3 The stress field inversion result of each time segmentThe black ‘+’ symbols denote the best fit principle stresses, while the red, green and blue dots denote the 95% confidence areas of maximum, intermediate and minimum principal stresses by 2000 times bootstrap resampling, respectively. (a) and (b) represent the damping and no damping inversion result respectively. Numbers represent each time period and consistent with the serial number in Table 2.

2.2 盈江地區(qū)應力場時空變化研究

為了進一步確定盈江地區(qū)地殼應力場是否存在時空變化，我們進行了如下兩方面研究.首先，為了確定盈江地區(qū)應力場在空間上是否存在變化，根據(jù)地震主要發(fā)生在大盈江斷裂和卡場—大竹寨附近(見圖1f)，大體上可以分為南北兩部分，因此，以北緯24.8°為界，將研究區(qū)劃分為南北兩個區(qū)域進行應力場反演，所得盈江地區(qū)南北部應力場結果見圖4a、4b和表3.反演得到的應力場最優(yōu)解相差較小且其95%置信區(qū)間相互重疊，可見南北部應力場可能存在較小的差異.其次，為了進一步確定南北部應力場差異的存在，我們將整體數(shù)據(jù)隨機劃分為與北部和南部區(qū)域震源機制解數(shù)量相同的數(shù)據(jù)集，分別進行應力場反演，并重復上述過程1000次，獲得這1000次應力場反演的平均反演殘差分布(見圖4d).從圖4d可見，總體作為一個區(qū)域反演應力場的平均殘差大于南北部應力場反演平均殘差以及隨機劃分數(shù)據(jù)為南北兩部分的應力場反演平均殘差，一定程度上是由于應力場擬合參數(shù)的增加導致反演殘差降低，同時也反映了應力場可能存在空間差異.南北地區(qū)震源機制解單獨反演應力場時的平均殘差小于隨機選取與南北部數(shù)量相等的震源機制解的應力場反演平均殘差，這個結果反映南北地區(qū)的應力場確實存在差異，即盈江地區(qū)南北部應力場存在空間差異；但從平均殘差的差異大小來看，這個差異較小，反映其南北部應力場接近.最后，為了進一步確定盈江地區(qū)應力場是否存在時空變化，我們從全部震源機制解中進行1000次隨機抽取26個(5個時間段震源機制解數(shù)的中位數(shù))震源機制解反演應力場，所得結果見圖5.可見，隨機抽取26個震源機制解反演應力場，反演的平均殘差總體上均滿足應力場均勻性要求(見圖5c)，反映任意抽取的地震集均滿足應力場均勻性假設，且總體上應力場反演結果較為穩(wěn)定(見圖5a).圖5a顯示有部分結果較為離散，為了定量地衡量隨機抽取26個地震所得應力場結果與區(qū)域應力場結果間的差異大小，我們用Kagan(1991)提出的方法計算他們與區(qū)域應力場間的3D旋轉角大小，所得結果見圖5b.由圖5b可見，隨機抽取26個地震所得應力場結果與區(qū)域應力場結果間差異較小，有78.3%的3D旋轉角小于等于近距離相似地震對間的3D旋轉角15°(Kagan，2000)，可見盈江地區(qū)在研究時段內(nèi)應力場是穩(wěn)定的.綜上可見，盈江地區(qū)的地殼應力場在我們研究的時段內(nèi)是穩(wěn)定的，且可以作為一個均勻應力場來看.

圖4 北部(a)、南部(b)、區(qū)域總體(c)應力場反演結果和平均殘差柱狀圖(d)圖(d)中紅線表示南北地區(qū)單獨反演應力場時的平均殘差，藍線表示全部數(shù)據(jù)反演應力場的平均殘差，灰色柱表示1000次隨機選取與南北地震數(shù)量一致震源機制解反演應力場時的平均殘差.圖例同圖3.Fig.4 The stress field inversion result of the Northern region (a), Southern region (b), the whole regional (c) of the study area and figure (d) is the histogram of The red line in Figure (d) represents the when all the focal mechanism solutions inverted together, and the histogram show 1000 tests in which the focal mechanism solutions were reshuffled and divided into two sets the same sizes as the northern and southern data sets. The legend is the same as Fig.3.

圖5 1000次隨機抽取26個震源機制解反演應力場結果圖(a)、與總體應力場間3D旋轉角(b)和平均殘差柱狀圖(c)圖(b)中紅線對應15°旋轉角.圖例同圖3.Fig.5 The 1000 times stress field inversion results of 26 focal mechanism solutions selected randomly from the whole data set (a), the 3-D rotation angle between the inversion result and the regional stress field (b) and the histogram of (c)The red line in figure (b) corresponds to the rotation angle of 15°. The legend is the same as Fig.3.

2.3 盈江地區(qū)總體應力場研究結果

鑒于前面研究表明盈江地區(qū)在研究時段內(nèi)存在時空均勻應力場，因此，我們使用本研究中全部震源機制解反演盈江地區(qū)的總體應力場.由表3可見，全部震源機制解旋轉角均方根值為38.5°，反映所用震源機制解基本上符合多樣性假設；總體數(shù)據(jù)反演的殘差平均值為23.2°，反映全部數(shù)據(jù)作為整體反演應力場是滿足應力場均勻性假設，說明用全部數(shù)據(jù)反演應力場時滿足應力場反演的兩個基本假設，所得應力場結果是可靠的.整個區(qū)域的應力場為走滑型(見圖4c)，可見總體應力場反演結果與各時段應力場反演結果是一致的.前面各時段反演結果表明應力場存在較小變化，且各時段應力場的95%置信區(qū)間存在重疊區(qū)域；整個研究區(qū)作總體應力場反演時，所用震源機制解依然滿足應力場均勻性要求，且反演結果和各時段結果相近.綜上可見，盈江地區(qū)在研究時段內(nèi)的應力場為走滑型，最大主應力軸呈NNE-SSW向，最小主應力軸呈SEE-NWW向，最大和最小主應力軸傾角近水平，中間主應力軸傾角近直立，該結果與先前的研究結果一致(錢曉東等，2011；Zhao et al.，2013；Xu et al.，2016；孫業(yè)君等，2017；Jin et al.，2019).

表3 研究區(qū)總體和南北部應力場反演結果表Table 3 The stress field inversion results of the whole study area and the Northern and Southern region

3 討論與結論

距離盈江縣～70 km的龍陵縣，在1976年5月29日曾發(fā)生7.3級和7.4級雙震，在本研究中，我們未考慮該雙震對研究區(qū)應力場的影響.一方面，先前的觀測條件有限，可以收集到的震源機制解資料較少，且不同方法給出的震源機制解誤差水平不一樣，因此，難以為應力場反演提供充足的約束以及給出可靠的誤差范圍估計.另一方面，先前研究表明大地震會引起震源區(qū)應力場旋轉，但在震后幾年內(nèi)主應力軸方位會恢復到震前方位.如1994年1月17日加州北嶺MW6.4地震引起震源區(qū)最大主應力軸發(fā)生17°的逆時針旋轉，在1995年底主應力軸方位恢復到震前方位(Zhao et al.，1997).1999年土耳其伊茲米特MW7.4地震發(fā)生后，震源區(qū)應力場由走滑型變?yōu)檎龜嘈?，?jīng)過兩個月的調(diào)整后，震源區(qū)應力場又恢復到震前的走滑型(Ickrath et al., 2014).2004年蘇門答臘—安達曼9.2級地震、2010年智利莫爾8.8級地震和2011年日本9.0級地震均引起了震源區(qū)主應力軸方位的劇烈變化，使得應力場為擠壓型的震源區(qū)在震后短期內(nèi)出現(xiàn)大量正斷層型地震，日本9.0級地震震源區(qū)、莫爾8.8級地震震源區(qū)和蘇門答臘—安達曼9.2級地震震源區(qū)南部主應力軸方位在震后幾個月內(nèi)又恢復到震前狀態(tài)(Hardebeck, 2012).基于上述兩方面原因，本研究未考慮1976年龍陵雙震對本研究區(qū)所研究時段的影響.

盈江地區(qū)各時段應力場反演，均滿足應力場均勻性假設和震源機制解多樣性假設，雖然阻尼和無阻尼應力場反演的最優(yōu)結果均顯示各時段應力場間存在變化，但是各時段應力場的95%置信區(qū)間相互重疊，說明這5次中強震對應力場產(chǎn)生的擾動較小.盈江地區(qū)作為整體反演應力場時，所用震源機制解能夠滿足多樣性假設，反演殘差也滿足應力場均勻性假設，因此，整體應力場反演結果是合理和可靠的，說明盈江地區(qū)在研究時段內(nèi)存在相對均勻的應力場.對盈江地區(qū)隨機抽取地震數(shù)據(jù)集反演應力場，反演平均殘差滿足應力場均勻性假設，最優(yōu)主應力軸分布較為集中，且與區(qū)域應力場結果間的3D旋轉角較小，反映盈江地殼在研究時段內(nèi)，該地區(qū)應力場的時空變化較小.綜上可見，5次中強震對該地區(qū)應力場產(chǎn)生的擾動較小，擾動后整個區(qū)域應力場仍可以視為均勻應力場，足見中強震對應力場擾動非常小.5次中強震中，震級最大的為MW5.82，根據(jù)震級和破裂尺度間經(jīng)驗關系，可知其破裂尺度為～10 km(Wells and Coppersmith，1994)；震級最小的地震MW～5.0，則破裂尺度～3 km，可見中強震的破裂尺度較小，則其對應力場的影響也較小，所以盈江地區(qū)雖然發(fā)生了5次中強震事件，但該地區(qū)的地殼應力場依然保持時空均勻性.這一結果表明，我們可以利用中強震后余震震源機制解獲取研究區(qū)的背景應力場.當主震對應力場產(chǎn)生較大擾動時，用余震獲得的應力場只是震后應力場，與背景應力場間存在差異(Hardebeck and Okada，2018；Sheng and Meng, 2020)，因此，利用余震反演應力場，所得應力場結果是反映區(qū)域應力場還是震源區(qū)震后應力場，取決于主震對應力場的影響.

在本研究中，我們使用了阻尼和無阻尼應力場反演方法反演各時段應力場，研究結果表明兩者所得的應力場存在一定的差異，因此，研究者在反演應力場時應該視具體問題的需要決定是否添加阻尼.當使用阻尼應力場反演時，所得應力場結果非常穩(wěn)健；雖然各時段震源機制解資料數(shù)量差異較大，但它們的置信區(qū)間均非常集中，反映使用阻尼應力場反演時，基于時空上應力場變化是連續(xù)和平滑的假設，在資料較少的網(wǎng)格點依然可以給出穩(wěn)健的結果.使用無阻尼應力場反演時，第1、3和4三個時段反演結果的置信區(qū)間非常離散，說明參與反演的震源機制解雖然大體上滿足震源機制解多樣性要求，但其數(shù)量相對較少，如Martínez-Garzón等(2016)的研究認為要獲得可靠的反演結果至少需要約30個震源機制解，Yoshida等反演應力場時要求網(wǎng)格點至少有10個震源機制解(Yoshida et al., 2014).特別是第3個時段兩種反演方式獲得的最優(yōu)解存在較大的差異，阻尼和無阻尼反演結果給出的最大主應力軸走向分別為208°和198°，相對與第2個時段最大主應力軸走向的旋轉角分別為4°和15°，可見，阻尼應力場反演方法會減小相鄰網(wǎng)格間應力場的差異(Hardebeck and Michael, 2006).Hardebeck和Michael(2006)用阻尼應力場反演方法研究1987年科林加(Coalinga)地震引起的最大水平主應力方向逆時針旋轉為～10°，無阻尼應力場反演結果得到的旋轉角為15°且與Michael(1984)的研究結果一致，上述結果表明阻尼應力場反演結果會減小時空相鄰網(wǎng)格點應力場的差異，這與本研究所得的結果是一致的.在研究地震是否引起應力場變化時，應謹慎使用阻尼應力場反演方法，因為，地震事件會引起震源區(qū)應力場的突然釋放，必然會對應力場時空連續(xù)性產(chǎn)生一定的影響，阻尼應力場反演會在一定程度上減小該影響.

當前應力場研究中，趨于將研究區(qū)細化為更小的網(wǎng)格，如研究云南地區(qū)地殼應力場時，Xu等(2016)使用0.7°的網(wǎng)格、孫業(yè)君等(2017)使用0.5°的網(wǎng)格以及Tian等(2019)使用0.3°的網(wǎng)格，以便獲得更為精細的應力場結果.但從本研究來看，總體可以作為一個區(qū)域來研究時，研究區(qū)進一步細化的必要性大大降低，或者說當前的資料與方法能觀測到的應力場下限是一定的.若要研究更為精細的應力場，首先務必要確保研究資料充足且滿足基本假設；再者，研究區(qū)或網(wǎng)格應力場反演殘差較大，進一步細化后應力場反演殘差明顯降低，或是當研究區(qū)存在明顯的構造差異時，則有必要進一步分區(qū)，進行更為精細的應力場研究，否則細化研究區(qū)不一定能獲得更為精確的應力場.

綜合上述研究結果，本研究獲得以下結論：(1) 5次中強震對盈江地區(qū)應力場產(chǎn)生的擾動較小，沒有引起盈江地區(qū)應力場在時間上的顯著變化；(2) 從隨機抽取地震反演應力場以及南北部地殼應力場研究結果可見，在研究時段內(nèi)，盈江地區(qū)地殼應力場在時空上可以視為均勻應力場；(3) 盈江地區(qū)地殼應力場為走滑型，最大主壓應力軸走向為NNE-SSW，最小主壓應力軸走向為SEE-NWW，最大和最小主應力軸傾角近水平，中間主應力軸傾角近直立.

致謝云南大學徐彥教授為本研究提供了震源機制解資料，本文部分圖件使用GMT繪圖軟件進行繪制(Wessel and Smith, 1991)，兩位匿名專家提出了寶貴的修改意見，特此致謝.