朱成林 賈 媛 李建奎 李 杰

1）山東省地震局，濟南 250014

2）山東省地震工程研究院，濟南 250021

剖面水準(zhǔn)測量的測線長、觀測周期長、序列長度短、觀測誤差較大、外界因素多等不利因素，給分析單測點的變化曲線帶來了困難。相比于跨斷層短水準(zhǔn)測量的定點觀測，剖面水準(zhǔn)測量更能從整體上反映斷裂帶在空間上和時間上的垂直形變規(guī)律。剖面測線上的各個測點間以及不同觀測周期的剖面曲線間是具有某種關(guān)聯(lián)性的，因此，可以利用這種關(guān)聯(lián)性，從整體上來研究斷裂帶垂直剖面形變的空間特征和時間分布特征，避開單獨分析各測點序列的不利因素。并利用灰色關(guān)聯(lián)分析方法研究了位于沂沭斷裂帶北段壽沂測線垂直剖面的形變規(guī)律。

利用灰色關(guān)聯(lián)分析分析測點曲線間及年度剖面曲線間的灰色關(guān)聯(lián)性。灰色關(guān)聯(lián)分析基于行為因子序列的微觀或宏觀幾何接近，通過分析曲線間幾何形狀來比較因素間的關(guān)聯(lián)性，是分析和確定因子間的影響程度的一種分析方法。設(shè)x0={x0(k)|k=1,2,…n}為參考序列；xi={xi(k)|k=1,2,…n}(i=1,2,…,m)為比較序列，則xi(k)與x0(k)的關(guān)聯(lián)系數(shù)可以表示為：

式中，ρ 為分辨系數(shù)，通常取ρ=0.5。關(guān)聯(lián)度γi按大小排序即為關(guān)聯(lián)序。

壽沂測線位于沂沭帶北段，自南向北主跨雙山-李家莊斷裂、上五井?dāng)嗔选⒁娑?無棣斷裂和昌樂-廣饒斷裂，自南向北的測點分別為壽沂8、壽沂6、壽沂5、壽沂3 和壽沂2。

分別以各測點觀測序列為參考序列，以所有測點的觀測序列為比較序列，可以求出測線上各測點序列的關(guān)聯(lián)度矩陣和關(guān)聯(lián)序。從測點序列間的關(guān)聯(lián)度及關(guān)聯(lián)序可以將所有測點按照變化趨勢可以分為2 組：一組為壽沂8、壽沂6 和壽沂3，垂直形變趨勢相對一致；另一組為壽沂5 與壽沂2。根據(jù)這些測點的位置分布認為這種分組現(xiàn)象與斷裂分布有關(guān)：益都-無棣斷裂和昌樂-廣饒斷裂的南側(cè)和北側(cè)分別具有同趨勢性。

將各個測點的形變量進行距離改化，即將各個測點的形變量除以該測點距離基準(zhǔn)點的距離，然后進行灰色關(guān)聯(lián)分析可以研究測線剖面垂直形變與距離之間的關(guān)系。距離改化前后關(guān)聯(lián)度矩陣變化說明壽沂測線各個測點的形變量與距離具有正相關(guān)性，亦說明了壽沂測線所處斷裂帶剖面具有整體傾斜特征，垂直形變及應(yīng)變隨距離增加而增強。

以同樣方法求得了不同年度剖面曲線間的關(guān)聯(lián)度矩陣和關(guān)聯(lián)序，據(jù)此發(fā)現(xiàn)壽沂測線垂直剖面形變在時間上具有以下時段性分布特征：①2001－2005 年垂直剖面曲線之間關(guān)聯(lián)度較大，具有趨勢上一致穩(wěn)定性，可以將這簇剖面曲線看作是該剖面垂直形變的正常背景變化；②2007－2008 年間垂直剖面曲線關(guān)聯(lián)度較大，而且從關(guān)聯(lián)序上看出與其2001－2005 年間剖面曲線簇的關(guān)聯(lián)度很低。說明了從時間特征上看，壽沂測線所描述的沂沭斷裂北段的垂直形變，也在2007－2008 年間出現(xiàn)了較大的應(yīng)力調(diào)整；③2009－2010 年間剖面曲線關(guān)聯(lián)度較大，從關(guān)聯(lián)序上看出其相比于2007－2008 年，與2001－2005 年間的剖面曲線簇關(guān)聯(lián)度較高，可以認為這一時段中的應(yīng)力處于釋放回調(diào)階段；④2011 年剖面曲線與以往各年份的剖面曲線關(guān)聯(lián)序基本上處于末尾，與2001－2005 年間的剖面曲線的關(guān)聯(lián)度明顯低于其他年份曲線的關(guān)聯(lián)度，可以看出2011 年以來壽沂測線所描述的沂沭斷裂北段的應(yīng)力調(diào)整處于最高位，需要進一步密切跟蹤監(jiān)測其垂直變化趨勢。

灰色關(guān)聯(lián)分析應(yīng)用于斷裂帶垂直剖面形變分析時，通過計算關(guān)聯(lián)度矩陣及關(guān)聯(lián)序可以排除分析個別測點序列的若干不利因素，從整體上分析斷裂帶的形變規(guī)律：①可以得到垂直形變的空間分布特性；②通過測點序列距離改化前后的關(guān)聯(lián)度矩陣的對比分析，可以直觀發(fā)現(xiàn)各測點各期形變量與距離的關(guān)系是否顯著，進而發(fā)現(xiàn)剖面測點各期變化是否具有整體傾斜性；③通過對各年份垂直剖面曲線的灰色關(guān)聯(lián)分析，可以從整體上分析垂直剖面形變的時間分布特征。