渭河盆地地形變觀測結果分析①

高好林，丁曉光，趙 韜，蘇利娜，古云鶴，關 昕，張永奇

（陜西省地震局，陜西 西安 710068）

0 引言

在渭河盆地開展地震流動重力和跨斷層形變測量已有三十多年，區域GPS網絡也積累了十幾年觀測資料。這三個觀測項目同屬于地形變監測，其觀測結果是否能共同解釋區域應力場變化和相關構造活動，進而為地震監測預報提供科學依據，是我們開展這些觀測項目的最終目的。

對于任何一種地震監測項目在研究其變化量和變化趨勢時，不同時間長度的選取所得觀測結果會有很大差別。時間越長，短期波動被抑制，長期變化和趨勢性變化得到顯現。

為此本文選取有GPS觀測結果以來較長時間段內以上各手段觀測結果進行對比分析，并結合其他資料和研究成果，探討渭河盆地地形變變化特征及其相互關系。

1 地質構造背景

位于鄂爾多斯塊體與秦嶺褶皺帶之間的渭河盆地位于秦嶺地背斜與鄂爾多斯地向斜之間，盆地南北兩側均為高角度正斷層所圍限，隴西旋卷構造體系東南撒開部分從西北邊插入，新華夏系的次級構造在盆地東部疊加復合，其相互穿插、干擾、遷就和利用構成盆地現今活動斷裂分布格局（圖1）。

在隴寶地區，隴西旋卷構造從隴縣向東南撒開插入秦嶺，該旋褶帶南端自西向東由桃園—龜川寺、固關—八渡、千陽—彪角和隴縣—岐山—馬召等四條斷裂帶組成。

據布設在寶雞市東側NS向穿越渭河盆地的地震流動觀測剖面資料，其接收函數偏移成像結果顯示，在鄂爾多斯西南緣的寶雞附近莫霍面存在間斷而不連續，表明青藏塊體對整體相對穩定的鄂爾多斯地塊具有持續而巨大的推擠作用力①陜西省地震局.2011年度陜西省地震趨勢研究報告.2010.。

在眉縣至岐山縣布設的一條近NS向深地震反射測線探測發現了10條斷裂。其中除相對應的秦嶺北側大斷裂、渭河斷裂和隴縣—岐山—馬召斷裂外，其余6條斷裂是新發現的，規模較小。這些斷裂根據各斷裂所錯斷的上部地層TE的形態都是正斷層，但是根據Tg及其下面的波組形態來看，岐山－馬召斷裂在淺部為正斷層，在深部可能為逆斷層，這表明該地區地層在古生代時期受壓應力作用，形成了一系列逆斷層；到新生代時期應力方向發生了變化，地層又受張應力作用形成了一系列的正斷層，老斷層也在張應力作用下形成了下逆上正的斷層性質［1］。

在長安區—咸陽市—禮泉縣布設的另一條深地震反射測線，也探測到規模大小不等的11條斷層，有的斷層深達康拉德界面。整個剖面反映出盆地南緣即秦嶺褶皺帶莫霍面呈俯沖形態伸展至渭河盆地下方莫霍面以下，而盆地下方即鄂爾多斯塊體南緣的莫霍面呈近似水平形態。同樣顯示出渭河盆地莫霍面南北向的不連續性［2］。

圖1 地形變監測網與活動斷裂分布圖Fig.1 Distribution of crustal deformation monitoring network and faults

以上探測結果表明渭河盆地西部斷層十分發育，尤其是岐山—馬召斷裂下逆上正的斷層性質，表明該斷裂深部構造的復雜性。在隴西旋卷構造南端，由于發育了一組NW向或NNW向的斷層，且該區域莫霍面間斷而不連續。于是在隴縣—寶雞—眉縣—周至—戶縣一帶布格重力異常等值線密集分布，形成一條布格重力異常高梯度帶，布格重力等值線呈NW展布，其中分別在寶雞和周至附近形成兩個中心區量值約－205×10－8m·s－2和－180×10－8m·s－2的布格重力異常凹陷區（圖2）。

布格重力異常梯度帶的存在反映了地殼中質量分布的不均勻性，從而使巖石圈因不同部位上所承受的負載不均一而發生彈性彎曲，導致它以下的上地幔所受到的靜壓力的不均衡，結果產生了垂直方向的作用力，造成了上地幔的下沉和上隆這種重力的均衡調整過程，而這種調整過程又必然引起地殼中物質的重新分布，結果又在地殼中造成了水平方向的作用力，這種橫向上的作用力在布格重力異常梯級帶或莫霍面的斜坡地帶為最大，這有可能是中強地震發生的動力來源之一。因此，布格重力異常梯級帶應該是我們識別中強地震潛在震源區的一個重要標志［3］。

2 地形變監測網絡概況與重力網平差方法

本文涉及到的渭河盆地地形變監測網絡由三部分組成：即陜西流動GPS網、關中地區流動重力網和跨斷層形變場地（圖1）。

2.1 陜西流動GPS監測網

2001年投入觀測的陜西流動GPS網西起隴縣，東至韓城，北達榆林，南抵留壩，由38個站點組成，主要站點分布在渭河盆地。每年第三季度觀測一期。由于陜北3個GPS點距本研究區域較遠，為了突出渭河盆地點網分布及相關構造，在圖1中將其舍去。

2.2 關中地區流動重力網及平差方法

圖2 渭河盆腔布格重力異常圖（據陜西省地質局）Fig.2 Bouguer gravity anomaly map of Weihe basin，Shannxi province（from Geological butrau of Shannxi province）

1992年投入觀測的關中地區高精度流動重力網覆蓋范圍東西與上相同，南至太白，北到黃龍，由11個閉合環、10條支線和116個測點組成。每期使用兩臺拉科斯特重力儀觀測，2008年以前每年觀測一期，其后每年兩期，至今已積累了二十余期觀測資料。

根據近十余年渭河盆地重力場分布圖像，筆者采用中國地震局地震研究所最新升級的流動重力資料處理軟件對近二十年的觀測資料進行擬穩平差處理。擬穩點的選取原則是：其一，擬穩點所在區域地質構造相對穩定；其二，在研究時間段內，點位沒有變動，觀測資料連續；其三，點位穩定可靠，觀測期間內沒有出現段差變化超限和測點重力變化異常，點值變化曲線平穩；其四，根據本文研究時段重力場變化分布特征，西部是重點關注區域，因此擬穩點應在東部選??；其五，大致位于研究區的中部，以減小傳遞誤差的影響。經反復比較分析所選擬穩點的組合方式以及所得重力變化圖像達到的最佳顯示效果。最終渭南附近位于斷層同一盤的幾個測點符合作為擬穩點的基本條件。

2.3 跨斷層形變場地

為了監測渭河盆地及鄰區的地震活動，1982年以后陜西省地震局先后布設了6處形變場地9條測線（包括口鎮定點形變1處兩條測線）。其后跨斷裂形變場地增至14處18條測線，其中王莊、康村、石堤峪、馬峪和廣濟等5個場地2000年投入觀測；灃峪口場地有3條測線分別布設在3條斷裂上；口鎮和大原場地各有兩條相連的測線跨過同一斷裂，其余場地均為一條測線。

跨斷裂形變測量每年復測三期以上，所得跨斷裂垂直形變資料均經尺長改正。

3 GPS觀測資料處理方法與應變場分布特征

3.1 GPS數據處理

GPS數據處理采用GAMIT10.4軟件。處理時聯合周邊15個IGS國際站，并進行了相應的各項誤差改正，得到精度量級在10－8的基線解協因數陣，最后融合全球IGS站的協因數陣，平差得到各站點在ITRF全球參考框架下的坐標。

利用GLOBK的GLRED命令可以處理多期數據并得到區域網站點的速率，此時的速率仍為全球框架下的速度。為了得到相對于歐亞板塊的站點運動，需要扣除歐亞板塊的整體旋轉，即公式（1）。

式中Vn、Ve為歐亞板塊旋轉引起的北分量、東分量線速率；λ、φ 和＝（Ωx，Ωy，Ωz）分別為歐亞板塊旋轉中心位置和旋轉角速度，（x，y，z）為GPS站點地心直角坐標；Ω為測站所處板塊的歐拉矢量。

3.2 GPS應變計算

利用上一步得到的相對于歐亞板塊的速度，就可以利用位移偏導與應變的關系式（3）計算應變值［4］：

式中，Δui為兩點間位移之差；Δxj為相應點間的經緯度之差。利用該式可得到主應變，面應變，最大剪應變，第一剪應變和第二剪應變。

求解位移之差可以有多種方法，比如三角形法，二次曲面法，最小二乘配置等，本文采用的是最小二乘配置法求解［5］。

3.3 GPS應變場分布特征

渭河盆地位于鄂爾多斯塊體與秦嶺褶皺帶之間，盆地內構造活動不但受到南北兩大塊體的牽制作用，而且來自青藏塊體東北緣NE向擠壓是鄂爾多斯地塊運動的主要驅動力，從而使鄂爾多斯地塊逆時針旋轉［6］。

2001－2010年渭河盆地GPS速度場顯示，鄂爾多斯塊體南緣所有站點運動方向為SE或SEE向。利用陜西地區流動GPS站的速度（或位移），可以得到2001－2010和2001－2012年主應變率、最大剪應變率和面應變率。

3.3.1 最大剪應變率

最大剪應變率的大小反映了構造活動的強烈程度。2001－2010年最大剪應變率空間分布顯示：在眉縣以北形成一個高值區，中心區在麟游與長武之間，約為3.0×10－8／a；其次，在黃陵以北和藍田以東分別形成兩個量值約為1.5×10－8／a和1.2×10－8／a的次高值區。盆地中東部腹地變化量在0.3×10－8／a～0.9×10－8／a之間。顯然西部活動水平明顯高于東部，東部南北兩側高于中間地帶（圖3（a））。

圖3 渭河盆地最大剪應變率與面應變率等值線圖Fig.3 Isogram of maximum shear strain rate and plane strain rate in Weihe basin

2001－2012年變化圖像顯示：西部高值區位置不變，量值僅為2.4×10－8／a，東部分別在宜君、三原和商洛附近各有一個中心量值約為1.8×10－8／a～2.1×10－8／a次高值區，似有東部增強西北減弱的跡象（圖3（b））。

3.3.2 面應變率

面應變率的正負值反映了局部地區的擠壓與拉張作用的相對大小，若為負值則以壓性活動為主，若為正則以張性活動為主［7］。

從2001－2010年分布圖像可以看出：除藍田－三原－旬邑一線兩側和 藍田至韓城一線以東為正值區外，渭河盆地面應變率以負值變化為主，并分別在岐山和蒲城形成兩個量值約為－1.9×10－8／a和－1.3×10－8／a的負值凹陷區（圖（3（c））。

2001－2012年圖顯示：前述東西部兩個面應變壓縮區已經貫通，并分別在鳳翔、戶縣和蒲城形成三個中心量值約為－1.5×10－8／a、－2.1×10－8／a和－0.9×10－8／a的強壓縮區。顯然戶縣附近壓縮變形較顯著（圖3（d））。

4 地形變觀測結果分析

4.1 GPS應變場與重力場變化

在相同時間段內，渭河盆地重力場變化圖像分析如下。

2001年5月－2010年11月，除乾縣、武功和三原所圍區域為小幅負向變化外，其余大部分區域為正向變化，西部重力變化明顯強于東部。隴縣－岐山－周至一線重力變化顯著，最大可達55×10－8m·s－2，并形成重力變化梯度帶。其次，蒲城附近也有一小范圍的正向變化區（圖4（a））。

圖4 重力場累積變化圖Fig.4 Accumulation changes in the gravitational field

2001年5月－2012年4月，雖然前述負值區仍存在，正值區范圍也有所擴大，但西部正值中心區位于周至以南，最大值僅為40×10－8m·s－2。重力場變化有減弱之勢（圖4（b））。

位于面應變壓縮區的重力點值變化曲線顯示，2001年以來處于波動上升變化中，累積變化幅度都在40×10－8m·s－2以上。其中位于面應變率強壓縮中心區的橫水測點累積重力變化量高達近80×10－8m·s－2（圖5）。顯然在面應變壓縮區強弱不同的區域重力變化量也存在明顯差異。

對比分析渭河盆地GPS應變場與重力場變化特征可以看出：GPS應變場和重力場空間分布有很好地相關性，尤其是面應變率與重力場變化相關性更好。即在剪應變率高值區和面應變率壓縮區重力場呈高梯度帶的正向變化；反之重力場則呈負向變化或微量變化。

在2001－2010年，以岐山為中心的面應變率負值區呈NW向展布，涵蓋渭河盆地中西部以及商洛地區。這一NW向負值區中心地帶在鳳翔－扶風—周至—戶縣一線，與區域重力場變化梯度帶基本重合。在此區域內，布格重力異常等值線也呈NW展布，等值線密集收縮在一個狹長的區域內（圖2）。顯然，由于這一時期該區域受到以擠壓為主的應力場作用，地殼介質表現為以擠壓變形為主，介質密度增加。于是在面應變收縮區重力場表現為正向變化，并形成一定規模的重力變化梯度帶［8］。

同樣，以蒲城為中心的面應變率負值區附近重力場也出現相應地正向變化。

以乾縣、武功和三原所圍區域的重力負值區位于東西兩大面應變率負值區的中間地帶，主要在面應變率的正值區。在面應變率正值區地殼處于拉張狀態，介質密度減小，于是重力場呈負向變化。

盡管2001－2012年面應變壓縮區連貫東西，與重力變化正值區仍然有較好地相關性。而且NW向的面應變率壓縮條帶與重力變化梯度帶依然存在，面應變強壓縮區與重力變化極值區也十分吻合。但與上期相比，最大剪應變率呈西部減弱東部增強的態勢，面應變強壓縮區轉移到周至、戶縣一帶，重力場變化也出現一定程度的回落。另外盆地西部主應變率方向由NE向轉為NEE向擠壓變形，東部北部為近EW向擠壓變形，南部為近EW向拉張變形。這些現象可能表明渭河盆地構造應力場目前處于調整階段，應變能仍有可能在渭河盆地西部或鄰區逐漸積累。

圖5 壓縮區重力點值變化曲線Fig.5 The change curves of gravity value in the compressed areas

4.2 GPS應變場與跨斷層形變

在圖3中，渭河盆地西部鳳翔—岐山—周至—戶縣一線形成面應變壓縮區，該區域西端最大剪應變率較大，主應變率和面應變率都表現為擠壓性質。隴西旋卷構造主旋帶南端的隴縣—岐山—馬召斷裂就在這一區域內。在此斷裂上有大塬、南大同和廣濟等三處跨斷層形變場地，圖6是其中兩個場地3個形變測段變化曲線（測段增量方向：自斷層上盤至下盤）。

位于鳳翔縣境內大塬2－3測段，2號點深埋于耕地中，無干擾破壞；3號點位于山坡基巖上，點位穩定可靠。兩點之間無過渡測點。2001－2013年，曲線經歷兩升兩降，最大波動幅度在7.3～－7.89 mm之間；2004年7月－2006年3月下降7.12 mm；于2007年7月恢復到高點位置以后，在此高點附近窄幅波動長達3年；2010年7月以后下盤相對上盤再次快速持續下降：2010年7月－2013年3月下降15.19mm，2001年3月－2013年3月下降7.89mm。在近三年時間內，斷裂活動表現為擠壓變形特征，逆繼承性構造運動顯著。

大塬4－1是該場地另一條跨斷層測線上的一個測段，兩個測點都是基巖點，其變化形態與前者基本同步，說明該場地跨斷層形變確實監測到了斷層的垂直形變信息。只是4號點2011年11月以后被破壞，該測段無法繼續觀測。

位于扶風縣境內南大同3－5測段，變化幅度在0.69～－3.1mm之間，比前者小很多，2011年3月－2013年3月下降2.15mm。廣濟場地相關測段2010年7月以后也急轉直下?？梢姅嗔鸦顒有再|與前者相同。至于這兩個場地變幅明顯小于大塬，可能有兩個原因：其一，大塬2－3和大塬4－1測段有關測點基本都在基巖上，能直接測得斷層的垂直形變；而后者是數百米厚的黃土沉積層測點，斷層微小的變化量被黃土層“淹沒”了。其二，岐山—馬召斷裂在鳳翔與扶風之間，斷層傾向發生了轉變，其西北段為SW65°，東南段為NE70°。由于斷層傾向相反，變化幅度也會存在差異。尤其是在青藏塊體北東向推擠和鄂爾多斯地塊逆時針旋轉的共同作用下，斷層在鳳翔以東運動受阻，大塬場地垂直形變加劇，更易聚集應變能。這或許就是該區域附近出現地形變異常的根源。

與此形成鮮明對比的是，主應變率和面應變率都表現為壓縮拉張性質兼具的渭河盆地東部，其主要活動斷裂上的垂直形變則表現為繼承性構造運動，甚至有的斷裂近兩年出現加速變化現象②陜西省地震局.2013年度陜西省地震趨勢研究報告.2012.。

由此可見，穿越面應變率強壓縮區的活動斷裂表現為逆繼承性構造運動，擠壓變形特征顯著；而在面應變膨脹區斷裂活動表現為繼承性構造運動，拉張變形特征明顯。

4.3 岐山—馬召斷裂逆繼承性構造運動分析

前文敘及，跨隴縣—岐山—馬召斷裂的大塬2－3和大塬4－1測段垂直形變曲線在2001年以后經歷過兩次較大幅度的下降變化，即斷裂表現為逆繼承性構造運動特征（圖6）。

在2004年7月－2006年3月的下降變化過程中，陜西省境內及鄰近區發生的最大地震分別為2006年7月26日鎮安3.6級和2004年8月10日甘肅迭部4.9級地震，形變場地距震中分別為212和323km；在2010年7月－2013年3月下降變化前后，2009年11月5日渭河盆地發生了高陵4.4級地震，鄰省發生了2013年7月22日甘肅岷縣6.6級地震，形變場地距震中分別為150和300km。

據分析，大塬2－3和大塬4－1測段兩次下降變化可能是甘肅迭部4.9和岷縣6.6地震孕育發生過程中的遠場效應（其震中相距約44km）。其實，該場地相關測段在1998年1月5日涇陽4.8級和2002年2月5日寶雞3.5級地震前后也曾出現較大幅度的下降變化，說明該場地具有較好的映震能 力（圖6和表1）［9－10］。

圖6 岐山—馬召斷裂垂直形變與地震對應關系Fig.6 Correspondence between vertical deformation of Qishan－Mazhao faults and the earthquakes

表1 大塬場地兩個測段在地震前后的下降幅度（垂直形變單位：mm）Table 1 The decline of vertical deformation of the two observation lines of Dayuan area before and after the earthquake（The unit of deformation：mm）

現以岷縣6.6級地震為例說明其映震機理：其一，渭河盆地西部與青藏塊體東北緣接壤，在斷裂構造上，南有西秦嶺北緣斷裂相連，北有隴西旋卷構造主旋帶——烏鞘嶺—六盤山旋褶帶自天祝經海原過六盤山進入本區，綿延數百公里，而隴縣—岐山—馬召斷裂就是其南端的一條分支斷裂。因此渭河盆地西部區域應力場的變化必然受制于青藏塊體東北緣地球動力學運動特征的影響。其二，六盤山－西秦嶺地區地形變觀測結果表明②該區在岷縣6.6級地震前處于顯著受壓應變積累狀態，斷裂活動呈現一定的擠壓閉鎖和壓扭性特征。2011－2012年GPS連續站和區域站觀測結果顯示海原—六盤山—西秦嶺構造區當前處于應變積累加速。區域水準觀測資料顯示西秦嶺北緣地區垂直形變背景異常顯著。區域流動重力場在2011年前后出現南部山區與北部盆地反向變化，在蘭州—青川一線以西形成NW向重力變化梯度帶。2010年10月～2011年10月該梯度帶以西為正向變化，以東為負向變化，岷縣6.6級地震震中區位于正值中心區附近，重力變化等值線約為40×10－8m·s－2左右；2011年10月－2012年10月正負區域反向，震中區位于負值中心區附近，重力變化等值線約為－70×10－8m·s－2左右。顯然，在2011年前后震中區附近重力差異變化較大，說明該區域附近構造變形和深部介質密度變化十分顯著。震中區附近重力場變化經歷了正向增強－反向減弱－發震的演化過程。如果把渭河盆地地形變觀測結果放在六盤山—西秦嶺地區地形變空間變化大背景中分析，可能其區域應力場變化具有統一的動力學背景。因此盆地西部地形變異常變化應該是六盤山—西秦嶺地區地形變空間分布的一部分，其異常變化可能是2013年7月22日甘肅岷縣6.6地震孕育過程中區域構造應力場逐漸增強，并在特定構造部位以地形變變化信息表現出來。其三，岷縣6.6級地震發生在岐山－馬召斷裂逆繼承性構造運動末期，表明該時段震中區附近區域應力場增強已達到臨震狀態。

根據我省現有的地形變觀測結果和其他地震前兆信息分析，在渭河盆地乃至陜西境內發生5.0級以上地震的條件似乎并不充分。往往某一區域存在地形變異常并不意味著該地區一定會發生強震［11］。因此岷縣6.6級地震的發生可能在一定程度上降低了在渭河盆地發生較大地震的可能性，但也不排除未來在渭河盆地西部及其鄰區有新的地震活動。

5 結論

（1）隴縣—岐山—周至—戶縣一線既是面應變壓縮區，又是重力變化梯度帶，還是布格重力異常梯度帶，說明這一區域應力場異常變化可能與深部構造有關。

（2）在主應變率以NE向擠壓變形為主的盆地西部，最大剪應變率高值區和面應變率壓縮區量值大范圍廣，且其二者的中心區基本重合。在穿越中心區附近的隴縣—岐山—馬召斷裂上，大塬和南大同等跨斷層形變場地2010年以來斷層活動表現為擠壓變形特征，逆繼承性構造運動顯著。在斷裂所經之處形成重力變化梯度帶，表明該時段渭河盆地西部尤其是岐山北部附近區域應力場有所增強，構造活動明顯。

（3）在面應變率壓縮區地殼處于壓縮狀態，介質密度增加，重力場呈正向變化，斷裂活動表現為擠壓變形特征。反之，重力場則呈負向變化或微量變化，斷裂活動表現為繼承性構造運動。

（4）岐山—馬召斷裂2004年7月以后的兩次逆繼承性構造運動可能是甘肅迭部4.9級和岷縣6.6級地震孕育發生過程中的遠場效應。

（References）

［1］任雋，馮希杰，王夫運，等.深地震反射剖面揭示的渭河盆地西部寶雞凸起斷裂深部特征［J］.國際地震動態，2012，85（6）：86.Ren J，Feng X J，Wang F Y，et al.The Characteristics of Deep Faults under Baoji Bulge in West of Weihe Basin by Deep Seismic Reflection Profile［J］.Recent Developments in World Seismology，2012，85（6）：86.（in Chinese）

［2］任雋，馮希杰，王夫運，等.深地震反射剖面揭示的渭河盆地西安坳陷的地殼精細結構［J］.地球物理學報，2013，56（2）：513－521.Ren J，Feng X J，Wang F Y，et al.Fine Crust Structures of Xi＇an Sag in the Weihe Basin Revealed by a Deep Seismic Reflection Profile［J］.Chinese Journal of Geophysics，2013，56（2）：513－521.（in Chinese）

［3］韓竹軍，鄔倫，于貴華，等.江淮地區布格重力異常與中強地震發生的構造環境分析［J］.中國地震，2002，18（3）：230－238.Han Z J，Wu L，Yu G H，et al.The High Gradients of BouguerGravity Anomalies and Seismotectonicsof Moderate Earthquakes in Jianghuai Region，China［J］.Earthquake Research in China，2002，18（3）：230－238.（in Chinese）

［4］李延興，黃珹，胡新康，等.板內塊體的剛性彈塑性運動模型與中國大陸主要塊體的應變狀態［J］.地震學報，2001，23（6）：565－572.Li Y X，Wu L，Huang C，Hu X K，et al.The Rigid and Elasticplastic Model of the Blocks in Intro－plate and Strain Status of Principal Blocks in the Continent of China［J］.Acta Seismologica Sinica，2001，23（6）：565－572.（in Chinese）

［5］張希，江在森.用最小二乘配置獲得地形變應變場動態圖像的幾個問題［J］.地殼形變與地震，1999，19（3）：32－39.Zhang X，Jiang Z S.Study on Some Questions of Dynamic Pictures of Crustal Deformation and Strain Fields Obtained by the Least Square Collocation［J］.Crustal Deformation and Earthquake，1999，19（3）：32－39.（in Chinese）

［6］范俊喜，馬瑾，甘衛軍.鄂爾多斯地塊運動的整體性與不同方向邊界活動的交替性［J］.中國科學（D輯），2003，33（增）：119－123.Fan J X，Ma J，Gan W J.Movement Integrity of Erodes Block and Its Boundary Activity Alternating in Different Directions［J］.Chinese Science（Series D），2003，33 （Supp）：119－123.（in Chinese）

［7］郭良遷，薄萬舉，胡新康.天津地區的水平應變場特征［J］.大地測量與地球動力學，2010，30（5）：8－14.Guo L Q，Bo W J，Hu X K.Characteristics of Horizontal Stra in Field in Tianjin［J］.Jouenal of Geodesy and Geodynamics，2010，30（5）：8－14.（in Chinese）

［8］劉代芹，李杰，王曉強，等.利用重力資料與GPS數據研究伽師及鄰區近期地殼運動特征［J］.中國地震，2011，27（2）：182－191.Liu D Q，Li J，Wang X Q，et.al.Characteristics of Recent Crustal Movements in Jiashi and Its Neighboring Regions on the Basis of Gravity and GPS Data［J］.Earthquake Research in China，2011，27（2）：182－191.（in Chinese）

［9］高好林，戴王強，曹建平，等.陜西涇陽4.8級地震前后重力變化與斷裂活動［J］.災害學，2006，21（1）：73－77.Gao H l，Dai W Q，Cao J P，et al.Gravity Variation and Fault Activity before and after the Jingyang MS4.8Earthquake in Shaanxi Province［J］.Journal of Catastrophology，2006，21（1）：73－77.（in Chinese）

［10］高好林，祝意青.重力段差變化與斷裂垂直位移的對比分析［J］.大地測量與地球動力學，2006，26（增）：31－36.GAO Hao－lin，ZHU Yi－qing.Comparative Analysis on Change of Gravity Segment Difference and Vertical Movement ofFault［J］.Jouenal of Geodesy and Geodynamics，2006，26（Supp）：31－36.（in Chinese）

［11］張四新，王雙緒，張希，等.地形變監測地震能力的評價［J］.西北地震學報，2011，33（1）：15－19.Zhang S X，Wang S X，Zhang X，et.al.Evaluation for Seismic Monitoring Ability of Crustal Deformation［J］.Northwestern Seismological Journal，2011，33（1）：15－19.（in Chinese）