盛書中，萬永革，黃驥超，卜玉菲，李祥

1中國地震局地球物理研究所，北京 100081

2防災(zāi)科技學院，河北三河 065201

3中國地震局蘭州地震研究所，蘭州 730000

4江蘇省地震局南京基準地震臺，南京 210014

1 引言

鄂爾多斯塊體是一個穩(wěn)定的塊體，位于南北地震帶北段東側(cè)、又處于華北、華南和青藏高原東北緣三個陸塊的匯聚部位，塊體內(nèi)部地震活動頻度和強度均低，但其周邊地震活動強烈.鄂爾多斯周緣的活斷層在我國大陸內(nèi)部自成一個體系，是一組十分具有活動特色的正斷層系，沿此帶歷史上發(fā)生過5次8級及8級以上強震，約占我國歷史上有記載的8級及8級以上強震的四分之一（鄧起東和尤惠川，1985；國家地震局《鄂爾多斯周緣活動斷裂系》課題組，1988）.由于鄂爾多斯塊體地理位置及地震活動的特殊性，有關(guān)于鄂爾多斯形成及其演化和形變動力學研究，一直以來都是地學研究的熱點問題，有些學者強調(diào)華南和華北大型走滑斷裂拉分作用是形成山西斷陷帶的重要原因（Tapponnier et al.，1982，1986；Peltzer et al.，1985），但該模型很難解釋近東西向展布的渭河和河套盆地的形成.有的學者認為鄂爾多斯周緣新生代斷陷盆地與華北東部地區(qū)盆地一起，其形成與太平洋板塊向亞洲大陸俯沖產(chǎn)生的弧后擴張有關(guān)，屬于弧后盆地（Uyeda and Kanamori，1979；Tian et al.，1992；Northrup et al.，1995）.有的學者認為鄂爾多斯塊體具有逆時針旋轉(zhuǎn)運動（蘇剛，1984；徐錫偉等，1994；邵輝成和蘇剛，1999；張岳橋等，2006），也有少數(shù)學者認為其存在順時針旋轉(zhuǎn)（張希等，2002）.“擠壓剪切模式”認為地塊西南側(cè)青藏塊體北東向的擠壓和鄂爾多斯塊體本身的隆升是決定鄂爾多斯地塊新生代構(gòu)造活動的兩大動力，并認為鄂爾多斯塊體本身不存在旋轉(zhuǎn)問題（鄧起東和尤惠川，1985；鄧起東等，1999）.陳小斌等（2005）提出了鄂爾多斯塊體與其周緣塊體相互作用的“擠推阻礙模式”.本研究將借助于地震波資料揭示鄂爾多斯塊體周緣現(xiàn)今地殼應(yīng)力場的空間分布，為進一步研究其周緣現(xiàn)今應(yīng)力應(yīng)變狀態(tài)、地震發(fā)震機理等研究提供基礎(chǔ)資料.

有關(guān)鄂爾多斯塊體周緣應(yīng)力場研究，已有許多研究成果.有的研究工作直接針對鄂爾多斯塊體周緣的地殼應(yīng)力場（薛宏運和鄢家全，1984；鄧起東和尤惠川，1985；鄧起東等，1999；范俊喜等，2003），并獲得了鄂爾多斯塊體周緣地殼應(yīng)力場的總體特征.

有的研究工作在研究較大區(qū)域地殼應(yīng)力場問題時，其中包含鄂爾多斯塊體周緣地區(qū)地殼應(yīng)力場的研究，如關(guān)于華北地殼應(yīng)力場的研究（高名修，1979；李欽祖，1980；許忠淮等，1983）、關(guān)于東亞地區(qū)現(xiàn)今構(gòu)造應(yīng)力場研究（許忠淮，2001）、關(guān)于中國大陸及鄰區(qū)應(yīng)力場的研究（鄢家全等，1979；張東寧和高龍生，1989；汪素云等，1996；杜興信和邵輝成，1999；謝富仁等，2004；徐紀人和趙志新，2006；朱守彪和石耀霖，2006；徐紀人等，2008；Wan，2010；范桃園等，2012）等均含有鄂爾多斯地區(qū)地殼應(yīng)力場的研究.

鄂爾多斯塊體周緣地質(zhì)構(gòu)造復雜，地震活躍，故也有大量關(guān)于鄂爾多斯周緣某個局部地區(qū)應(yīng)力場研究工作.有許多學者對寧夏地區(qū)應(yīng)力場進行了研究（陳愛玲等，1981；李孟鑾和趙知軍，1981；李孟鑾和趙知軍，1986；趙知軍和劉秀景，1990；趙知軍等，2002），獲得了寧夏地區(qū)較詳細的應(yīng)力場特征.許多研究者利用地球物理和地震觀測等手段，對山西地區(qū)應(yīng)力場進行了較深入的研究（梁海慶等，1986；胡惠民，1987；成爾林和張美芳，1990；劉巍等，1993，1994，1995，1996；安美建和李方全，1998）.姜家蘭（1991）利用地震震源機制解資料（其中包括小地震綜合斷層面解）給出了陜西及鄰近地區(qū)現(xiàn)代構(gòu)造應(yīng)力場，研究結(jié)果表明該區(qū)域應(yīng)力場從北向南存在差異性.謝富仁等（2000）通過斷層滑動資料反演構(gòu)造應(yīng)力場，確定了海原、六盤山斷裂帶至銀川斷陷的應(yīng)力場.

關(guān)于鄂爾多斯塊體周緣的地殼應(yīng)力場研究，顯示出該區(qū)域應(yīng)力場的復雜性，總體上獲得了較為一致的認識，即本區(qū)域應(yīng)力場為北東—北東東向主壓應(yīng)力和北西—北北西向拉張應(yīng)力；主壓應(yīng)力方位，自北向南，自西向東呈現(xiàn)出規(guī)律性變化，這種應(yīng)力場格局與鄂爾多斯塊體同時兼受印度板塊和太平洋板塊擠壓并以前者為主有關(guān)（薛宏運和鄢家全，1984；李孟鑾，1985；胡惠民，1989；鄧起東等，1999；范俊喜等，2003；徐黎明等，2006）.但上述研究中，由于受到資料有限性的限制，給出的應(yīng)力場結(jié)果不夠精細，不能詳細地反映地質(zhì)構(gòu)造相對應(yīng)的構(gòu)造應(yīng)力場的變化；隨著地震觀測臺網(wǎng)的加密和地震資料的豐富，為鄂爾多斯周緣地殼精細應(yīng)力場的深入研究提供了機會.本研究采用了近期發(fā)生的大量地震的P波初動資料，以期得到更為精細的地殼應(yīng)力場，為鄂爾多斯地塊周緣地殼形變的動力學機制、地震孕育機理等研究提供基礎(chǔ)資料.

2 方法、資料和地殼速度模型

Aki（1966）首先提出采用大量小震P波初動資料推測應(yīng)力場方向的方法，該方法被稱為求解應(yīng)力場的綜合震源機制解法.該方法的優(yōu)點是可以利用大量不能單獨確定震源機制解的小地震P波初動極性數(shù)據(jù)推斷應(yīng)力場方向，因此，得到了較為廣泛的應(yīng)用（李欽祖等，1973；許忠淮等，1979，1983，1987，1989；薛宏運和鄢家全，1984；成爾林和張美芳，1990；姜家蘭，1991；徐鳴潔等，1996；Wan，2010；卜玉菲等，2013；盛書中等，2013）.該方法還被用于分析區(qū)域應(yīng)力場的變化過程，已有不少學者在這方面做了許多研究工作（華祥文，1980；成爾林等，1982；闞榮舉等，1983）.李欽祖等（1973）根據(jù)該方法求得了紅山臺和沙城臺兩個區(qū)域的應(yīng)力場.許忠淮等（1983）進一步提出了采用研究區(qū)域的多個地震和多個臺站求解區(qū)域應(yīng)力場的方法，并采用合成數(shù)據(jù)對此方法進行了嚴格測試，驗證了用該方法可以獲得地殼應(yīng)力場.薛宏運和鄢家全（1984）將該方法用于鄂爾多斯地塊周圍現(xiàn)代地殼應(yīng)力場的確定.汪素云和許忠淮（1985）用該方法給出了中國東部大陸的構(gòu)造應(yīng)力場特征.萬永革等（2011）利用臺灣自1991—2007年間觀測的P波初動資料給出了臺灣地區(qū)地殼三維應(yīng)力結(jié)構(gòu)的初步結(jié)果.隨著地震觀測臺網(wǎng)的密集，地震監(jiān)測能力的增強，求解應(yīng)力場的綜合震源機制解法將有著良好的應(yīng)用前景.

本研究中使用了研究區(qū)內(nèi)2007年8月1日—2013年7月21日鄂爾多斯塊體周緣發(fā)生的8499個地震的49844個P波初動符號，其中2007年8月1日—2011年8月31日的P波初動資料為項目組由中國地震局地球物理研究所數(shù)據(jù)備份中心（鄭秀芬等，2009）申請數(shù)據(jù)后自行拾取的，2011年8月31日之后的P波初動極性資料由中國地震臺網(wǎng)中心震相報告中讀取，研究中使用的地震震中分布情況見圖1.

我們對所用的地震震級進行統(tǒng)計，不同震級檔的地震個數(shù)及其比例見圖2和表1.可見該段時間內(nèi)，該區(qū)域2級以下地震約占66%，b值一般情況下在一個地區(qū)是較為穩(wěn)定的，也就是說我們所統(tǒng)計的震級分布，大體上可以視為該區(qū)域長期的震級分布情況.由該區(qū)域地震震級分布可見，在未考慮完整性震級的情況下，M2.0級以下小震數(shù)目占了六成以上.一般情況下小地震難以給出準確的震源機制解，為了充分利用該區(qū)域發(fā)生的大量小震資料來約束應(yīng)力場方向，綜合震源機制解法無疑是最好的選擇，故本研究應(yīng)用P波初動綜合震源機制解法研究鄂爾多斯周緣地區(qū)地殼應(yīng)力場方向.

表1 鄂爾多斯周緣地震震級分布統(tǒng)計表（2007－8－1—2013－7－21，M＞0）Table 1 The earthquake magnitude statistic table of Ordos and its adjacent area

本研究中應(yīng)力場計算仿照萬永革等（2011）的應(yīng)力場計算方法，將研究區(qū)域（經(jīng)緯度范圍分別為104°E—116°E，33°N—42°N）劃分為0.25°×0.25°的二維網(wǎng)格，根據(jù)網(wǎng)格點選取周圍地震臺站上觀測到的P波初動符號計算每個網(wǎng)格點的應(yīng)力場.由于每個地震距離網(wǎng)格點的距離不同，對所計算的網(wǎng)格點應(yīng)力場方向提供的約束也不一樣，因此我們根據(jù)其與網(wǎng)格點的距離給予不同地震P波初動符號不同權(quán)重.每個地震的P波初動符號權(quán)重W仿照Shen等（1996）的大地測量數(shù)據(jù)計算應(yīng)變的方式求?。?/p>

其中，D為距離衰減常數(shù)；r為折合距離，其值按式（2）計算.距離衰減常數(shù)D的選取直接影響到距離計算網(wǎng)格點不同距離地震P波初動符號的權(quán)重，從而影響到各網(wǎng)格點應(yīng)力場計算結(jié)果，同時，D的選取對應(yīng)力場計算結(jié)果也起到一定的平滑作用.網(wǎng)格點應(yīng)力場計算時，如果D取值較小，網(wǎng)格較大，意味著僅考慮距離網(wǎng)格點較近地震的P波初動符號，網(wǎng)格點附近網(wǎng)格內(nèi)地震P波初動符號沒有考慮或權(quán)重較小，所得各個網(wǎng)格點間應(yīng)力場計算結(jié)果相對較為獨立；反之，如果D取值較大，網(wǎng)格較小，意味著網(wǎng)格點附近網(wǎng)格中地震的P波初動符號也有一定的權(quán)重，那么所得各個網(wǎng)格點間應(yīng)力場則相對較為平滑，即周圍網(wǎng)格點的地震P波初動對其應(yīng)力場計算存在約束.在本研究中考慮以下兩個方面：一、對各網(wǎng)格點應(yīng)力場計算結(jié)果給予一定的平滑作用；二、部分網(wǎng)格點內(nèi)地震資料較少，充分利用周邊網(wǎng)格點內(nèi)地震資料給予約束，以便給出整個環(huán)鄂爾多斯塊體周緣的地殼應(yīng)力場，將D取為50km.

圖1 地震震中及斷層分布圖圓表示震中位置，黑線表示斷層.F1牛首山北麓斷裂；F2黃河斷裂；F3銀川—平羅斷裂；F4賀蘭山東麓斷裂帶；F5正誼關(guān)斷裂；F6狼山—巴音烏拉山東麓斷裂；F7色爾騰山山前斷裂；F8鄂爾多斯北緣斷裂；F9烏拉山山前斷裂；F10大青山山前斷裂；F11和林格爾斷裂；F12蠻漢山山前斷裂帶；F13岱?！S旗海南緣斷裂帶；F14口泉斷裂；F15桑干河斷裂；F16恒山北麓斷裂；F17五臺山北麓斷裂；F18系舟山山前斷裂；F19交城斷裂；F20太谷斷裂；F21霍山斷裂；F22羅云山山前斷裂；F23峨眉臺地北緣斷裂；F24臨猗斷裂；F25中條山北麓斷裂；F26渭河斷裂；F27秦嶺北麓斷裂；F28華山山前斷裂；F29六盤山東麓斷裂；F30南、西華山北麓斷裂；F31天景山北麓斷裂；F32煙筒山東麓斷裂；F33會寧—文崗斷裂；F34馬銜山斷裂；F35臨潭—宕昌斷裂；F36禮縣—羅家堡斷裂；F37麻河沿斷裂；F38成縣盆地南緣斷裂；F39成縣盆地北緣斷裂；F40光蓋山—迭山北麓斷裂；F41光蓋山—迭山南麓斷裂；F42迭部—白龍江斷裂；F43文縣斷裂；F44平武—青川斷裂.Fig.1 The distribution of earthquake epicenter and seismic faultCircles represent the epicenter，black lines represent fault.F1Niushoushan northern foot fault zone，F(xiàn)2Huanghe fault，F(xiàn)3Yinchuan－Pingluo fault，F(xiàn)4Eastern foot fault zone of Helanshan，F(xiàn)5Zhengyiguan fault，F(xiàn)6Eastern foot fault of Langshan－Bayinwulashan，F(xiàn)7 Piedmont fault of Seerteng mountains，F(xiàn)8Northen marginal fault of Ordos，F(xiàn)9Piedmont fault of Wulashan，F(xiàn)10Piedmont fault of Daqingshan，F(xiàn)11Horinger fault，F(xiàn)12Piedmont fault of Manhanshan，F(xiàn)13Southern marginal fault zone of Daihai－Huangqihai，F(xiàn)14Kouquan fault，F(xiàn)15Sanggan River fault，F(xiàn)16Northen foot fault of Hengshan，F(xiàn)17Northen foot fault of Wutaishan，F(xiàn)18Piedmont fault of Xizhoushan，F(xiàn)19Jiaocheng fault，F(xiàn)20Taigu fault，F(xiàn)21Huoshan fault，F(xiàn)22Piedmont fault of Luoyunshan，F(xiàn)23Northern marginal fault of Emei platform，F(xiàn)24Linyi fault，F(xiàn)25Northen foot fault of Zhongtiaoshan，F(xiàn)26Weihe fault，F(xiàn)27Northen foot fault of Qinling，F(xiàn)28Piedmont fault of Huashan，F(xiàn)29Eastern foot fault of Liupanshan，F(xiàn)30Northen foot fault of Nanhuashan and Xihuashan，F(xiàn)31Northen foot fault of Tianjingshan，F(xiàn)32Eastern foot fault of Yantongshan，F(xiàn)33Huining－Wengang fault，F(xiàn)34Maxianshan fault，F(xiàn)35Lintan—Dangchang fault，F(xiàn)36 Lixian—Luojiabu fault， F37 Maheyan fault， F38 Southern boundary fault of Chengxian Basin，F(xiàn)39Northern boundary fault of Chengxian Basin，F(xiàn)40Northen foot fault of Guanggaishan－Dieshan，F(xiàn)41Southen foot fault of Guanggaishan－Dieshan，F(xiàn)42Diebu－Bailongjiang fault，F(xiàn)43Wenxian fault，F(xiàn)44Pingwu－Qingchuan fault.

圖2 地震震級（a）和深度（b）分布柱狀圖Fig.2 The histogram of earthquake magnitude（a）and depth（b）

其中，φ、θ、z為地震的經(jīng)度、緯度和深度，φ0、θ0、z0為網(wǎng)格點的經(jīng)度、緯度和深度.我們根據(jù)所用地震事件的深度分布情況（見圖2b），將網(wǎng)格點深度z0取為10km，即應(yīng)力場計算結(jié)果為地殼內(nèi)10km處應(yīng)力場情況.由P波初動符號權(quán)重隨距離變化關(guān)系可知，距網(wǎng)格點折合距離為41.6km的地震的P波初動數(shù)據(jù)只有50%的權(quán)重，距網(wǎng)格點折合距離為75.9km的地震的P波初動數(shù)據(jù)只有10%的權(quán)重，而折合距離大于100km地震的P波初動數(shù)據(jù)的權(quán)重為2%.在本研究中，當距離超過100km時其權(quán)重設(shè)置為0（即未考慮折合距離大于100km的地震）.為了保證應(yīng)力場的反演質(zhì)量，我們只選取P波初動符號個數(shù)大于等于100個的網(wǎng)格點進行反演.

為了確定區(qū)域內(nèi)與P波初動數(shù)據(jù)擬合最好的綜合震源機制解，進而分析其平均P、T軸方向，我們采用1°×1°×1°的網(wǎng)格搜索P、B、T軸方位.每一次搜索給出與綜合震源機制解模型不符合的P波初動符號數(shù)（考慮權(quán)重）與總P波初動符號個數(shù)（考慮權(quán)重）之比，我們稱之為矛盾比.選擇最小矛盾比所對應(yīng)的P、B、T軸方位作為該區(qū)域的主應(yīng)力方位.由于該方法綜合了多個地震的P波初動數(shù)據(jù)，得到的綜合震源機制解的P和T軸方位大體可以認為是該研究區(qū)域的最大和最小主壓應(yīng)力方位（李欽祖等，1973；許忠淮等，1979，1983，1989）.當P波初動數(shù)據(jù)足夠多，且有很好的方位覆蓋，則可以獲得約束較好的綜合震源機制解，即給出網(wǎng)格點處最優(yōu)且可靠的應(yīng)力場方向.

本研究中用于計算P波射線離源角和方位角所使用的速度模型參數(shù)見表2，該模型參數(shù)依據(jù)由Crust2.0中獲得.

表2 鄂爾多斯塊體地殼速度模型Table 2 The crustal velocity model of the Ordos block

3 應(yīng)力場計算結(jié)果及分析

圖3 P軸方位和網(wǎng)格點使用P波初動符號數(shù)分布圖Fig.3 The distribution of the azimuth of P axis and the number of P wave first motion data used in each grid

應(yīng)用上文所述資料和方法，我們對研究區(qū)域進行了經(jīng)緯度網(wǎng)格為0.25°×0.25°的空間掃描，獲得了每個網(wǎng)格點內(nèi)的最優(yōu)綜合震源機制解，所得的綜合震源機制解的P、T軸方位空間分布、各個網(wǎng)格點所用數(shù)據(jù)個數(shù)和矛盾比情況見圖3和圖4.以下的圖件中，為了避免初動符號數(shù)過少對結(jié)果可靠性的影響，我們舍棄了網(wǎng)格內(nèi)數(shù)據(jù)加權(quán)后初動符號數(shù)小于50的數(shù)據(jù).為了去除地震發(fā)生較少的網(wǎng)格點，以及避免結(jié)果過分依賴于周邊區(qū)域地震數(shù)據(jù)，舍棄了網(wǎng)格點內(nèi)數(shù)據(jù)平均權(quán)重小于0.15網(wǎng)格點的反演結(jié)果.由圖3可見，各個網(wǎng)格點內(nèi)P波初動符號數(shù)據(jù)量較為充足，為應(yīng)力場反演結(jié)果的可靠性奠定了基礎(chǔ)，總體來說受地震活動性影響，東緣各個網(wǎng)格點使用的P波初動數(shù)據(jù)量較大，西緣次之，北緣地區(qū)P波初動數(shù)據(jù)量最少.由圖4可見，在整個鄂爾多斯周緣地區(qū)，各個網(wǎng)格內(nèi)的矛盾比較高，特別是在其西南緣地區(qū)矛盾比最高，一定程度上反映了該區(qū)域應(yīng)力場的復雜性，具體原因還有待于進一步分析.

參照世界應(yīng)力圖的劃分原則（Zoback，1992），根據(jù)震源機制解3個應(yīng)力軸傾角的大小，將震源機制解類型（即其反映的應(yīng)力狀態(tài)）分為6種：正斷型（NF）、正走滑型 （NS）、走滑型 （SS）、逆走滑型（TS）、逆斷型（TF）和不確定型（U）.具體分類標準見表3.根據(jù)上述分類標準，用不同顏色表示不同類型的震源機制解，不同類型震源機制解的分布情況見圖5.

以下我們根據(jù)P軸方位分布圖、T軸方位分布圖和綜合震源機制解分布圖，由銀川—吉蘭泰斷陷帶開始，沿順時針方向，對鄂爾多斯周緣地殼應(yīng)力場計算結(jié)果進行分析.

表3 震源機制解類型劃分表Table 3 Categories of tectonic stress regime for focal mechanism

圖4 T軸方位和矛盾比分布圖Fig.4 The distribution of the azimuth of Taxis and the contradictory ratio

3.1 銀川—吉蘭泰斷陷帶內(nèi)應(yīng)力場結(jié)果及分析

銀川—吉蘭泰斷陷帶位于鄂爾多斯塊體西北側(cè)，南起牛首山，北達狼山以南，由銀川盆地、賀蘭山地壘和吉蘭泰盆地作右行雁裂式排列組成，斷層十分發(fā)育，在第四紀大多有不同程度的活動，它們控制著區(qū)域地貌及沉積的形成和演化.按走向與性質(zhì)，可將斷裂分為三組，即北北東至北東向斷裂、北西向斷裂和東西向斷裂.本斷陷帶內(nèi)北北東至北東向斷裂是第四紀以來活動性最強的一組斷裂，它不僅控制地塹的形成，而且與地震活動關(guān)系也最緊密，其性質(zhì)為張性右旋；北西西向斷裂性質(zhì)為壓性左旋（鄧起東和尤惠川，1985；國家地震局《鄂爾多斯周緣活動斷裂系》課題組，1988）.由圖3可見，在該區(qū)域內(nèi)主壓應(yīng)力軸方位從南向北，總體上由北東向逐漸變?yōu)楸北睎|向，與先前的應(yīng)力場研究結(jié)果較為一致（李孟鑾和趙知軍，1981，1986；趙知軍和劉秀景，1990；趙知軍等，2002；許忠淮等，1989；范俊喜等，2003）；同時也存在一定差異，具體表現(xiàn)為在銀川以北至烏海以西地區(qū)，前人研究給出的P軸方位為北北東向（薛宏運和鄢家全，1984；趙知軍和劉秀景，1990；范俊喜等，2003），我們給出的P軸方位為近南北向，且傾角較大，差異的原因可能和應(yīng)力場計算分區(qū)大小關(guān)系較大，先前的結(jié)果均為較大區(qū)域的平均結(jié)果，本研究結(jié)果在分區(qū)上更為精細.由圖4可見，P軸存在差異地區(qū)的T軸方位為北西西至北西向，傾角近水平，在銀川盆地和吉蘭泰盆地內(nèi)T軸走向大體上與盆地走向相垂直，反映了盆地的拉張狀態(tài)，故上述差異可能是我們的分區(qū)更小，體現(xiàn)了該區(qū)域的局部應(yīng)力場特征.在圖5中，綜合震源機制解在銀川盆地和吉蘭泰盆地東南為正斷層，且走向大體為北北東向，和張性盆地及主要斷裂的展布相一致，進一步說明了拉張正斷層的性質(zhì).在石嘴山至烏海附近區(qū)域，震源機制解類型以不確定型為主，目前推測其可能和該處有近東西向的正宜關(guān)斷裂有關(guān)，正宜關(guān)斷裂現(xiàn)代活動為擠壓兼左旋走滑為主（國家地震局《鄂爾多斯周緣活動斷裂系》課題組，1988），該區(qū)域的應(yīng)力場主要為拉張走滑型，而應(yīng)力場反演的前提條件是被反演應(yīng)力場的研究區(qū)域具有均勻的應(yīng)力場，因此可能是此處應(yīng)力場的不均勻性所致，具體原因還有待進一步分析.在磴口以西地區(qū)，綜合震源機制解類型表現(xiàn)為走滑型為主，該走滑型機制解也表現(xiàn)出一定的逆沖分量，且有兩個為近純逆沖型，綜合震源機制解類型和該區(qū)域的狼山—巴音烏拉山東麓斷裂的正斷層性質(zhì)不符.但從P、T軸方位來看，T軸走向與狼山—巴音烏拉山斷裂相垂直（圖4），P軸的走向與主要斷裂相平行（圖3），與先前研究結(jié)果較為一致（趙知軍和劉秀景，1990；趙知軍等，2002）.

圖5 綜合震源機制解分布圖Fig.5 The distribution of composite focal mechanism

3.2 河套斷陷帶內(nèi)應(yīng)力場結(jié)果及分析

河套斷陷帶位于陰山隆起和鄂爾多斯隆起之間，西界為狼山山前斷裂，東界是和林格爾斷裂，北界為陰山山前斷裂，南界為鄂爾多斯北緣斷裂.廣闊的河套盆地東西長約440km，南北寬約40～80km，總體走向近東西.河套斷陷帶可劃分為三個凹陷，即臨河凹陷、白彥花凹陷和呼和凹陷.河套斷陷帶周邊主要發(fā)育正斷層，局部見張扭性正斷層，主要顯示左旋扭動特征，其新構(gòu)造運動以垂直差異運動為主.由圖4可見，T軸方位為北北西向，與先前研究結(jié)果較為一致（許忠淮等，1989；范俊喜等，2003；徐紀人和趙志新，2006），P軸方位為北北東向，但在包頭附近P軸方位出現(xiàn)了北北西向，這一現(xiàn)象在先前分區(qū)較大的研究中沒有出現(xiàn)（許忠淮等，1989；范俊喜等，2003），但在中國大陸現(xiàn)代應(yīng)力場圖中也有類似情況（謝富仁等；2004），我們認為先前研究僅給出了總體趨勢，本次研究對研究區(qū)域的細分，顯示出局部應(yīng)力場的特性和應(yīng)力場的不均勻性.圖5中，在整個河套斷陷帶內(nèi)，除臨河盆地西部綜合震源機制解結(jié)果為不確定型，其余區(qū)域的綜合震源機制解結(jié)果均為正斷層型，且綜合震源機制解的節(jié)面走向大體與邊界主要斷層走向相一致，節(jié)面由西部的近東西向，向東逐漸變?yōu)楸睎|向，與河套斷陷帶的東界和林格爾斷裂走向相一致，綜合震源機制解結(jié)果反映了河套斷陷帶現(xiàn)今應(yīng)力狀態(tài)為拉張.

3.3 岱海斷陷帶內(nèi)應(yīng)力場結(jié)果及分析

在岱海斷陷帶內(nèi)，控制岱海斷陷帶的主要斷裂為蠻漢山山前斷裂帶、鞍子山山前斷裂帶和岱?！S旗海南緣斷裂帶，它們沿北東東向展布，均為張性或張扭性正斷層（國家地震局《鄂爾多斯周緣活動斷裂系》課題組，1988）.由圖4可見，T軸的方位和控制岱海斷陷帶的主要斷裂近垂直，且傾角近水平，而此處P軸（見圖3）的傾角較大，說明岱海斷陷帶的應(yīng)力類型為拉張型應(yīng)力場，P、T軸方位與先前的研究結(jié)果相一致（薛宏運和鄢家全，1984；許忠淮等，1989；范俊喜等，2003；謝富仁等；2004）.圖5中綜合震源機制解結(jié)果為正斷層，且綜合震源機制解的節(jié)面大體上呈北東東向，該結(jié)果和控制岱海斷陷帶的主要斷裂在性質(zhì)及走向上均相一致，反映了岱海斷陷帶現(xiàn)今應(yīng)力狀態(tài)為拉張型.

3.4 山西斷陷帶內(nèi)應(yīng)力場結(jié)果及分析

山西斷陷帶位于鄂爾多斯塊體的東側(cè)，是由十來個大小不等的北東、北東東走向的地塹或半地塹式盆地作右行斜列組成，總體走向北北東，平面上呈“S”形.該帶與陰山隆起相鄰，南與渭河斷陷帶東端的運城盆地以隆起相隔，西為鄂爾多斯塊體東部的呂梁山隆起，東為太行山隆起.山西斷陷帶盆地的排列形式、盆地邊界的對稱分布特征、盆地邊緣山地與平原地貌的強烈對照性，均是斷陷帶整體作右旋剪切拉張的產(chǎn)物，山西斷陷帶是一條右旋剪切拉張帶（鄧起東和尤惠川，1985；胡惠民，1987；國家地震局《鄂爾多斯周緣活動斷裂系》課題組，1988）.由圖3和4可見，在山西斷陷帶內(nèi)P、T軸方位的空間分布總體趨勢和先前研究結(jié)果較為一致（梁海慶等，1986；許忠淮等，1989；成爾林和張美芳，1990；劉巍等，1993；鄧起東等，1999；范俊喜等，2003；謝富仁等；2004；徐黎明等，2006）.在圖5中，山西斷陷帶綜合震源機制解總體上為近純正斷層型地震，綜合震源機制解的節(jié)面大體呈北東東向展布，在各個盆地內(nèi)部節(jié)面隨盆地邊界的主要控制斷裂的走向變化.在臨汾盆地北端和太原盆地以南一段，震源機制解較為復雜.其中臨汾盆地北端的綜合震源機制解類型為正斷走滑型以及個別震源機制解表現(xiàn)為走滑型，節(jié)面走向也有別于臨汾盆地的純正斷層機制解，該處綜合震源機制解類型變化可能與霍山山前斷裂以及靈石隆起相關(guān).霍山山前斷裂發(fā)育在臨汾盆地東緣北段，是山西斷陷帶中結(jié)構(gòu)最復雜的張扭性斷裂，總體走向15°～20°，斷面傾向北西，傾角65°～75°.可見此處綜合震源機制解和霍山山前斷裂性質(zhì)相似，我們認為此處應(yīng)力場也主要受到霍山山前斷裂的影響而出現(xiàn)了變化.在霍州附近綜合震源機制解結(jié)構(gòu)為不確定型，呈現(xiàn)出不同于山西地塹帶的局部應(yīng)力場特征，由圖中構(gòu)造可見該處的地質(zhì)構(gòu)造復雜，斷裂走向與周邊的主要斷裂走向存在差異，而近東西向的橫向斷裂的第四紀活動以左旋正傾滑為主，個別有逆傾滑（國家地震局《鄂爾多斯周緣活動斷裂系》課題組，1988），此處也恰好為太原盆地和臨汾盆地間的靈石隆起，而盆地間的橫向隆起則為剪切帶次級破裂之間的巖橋區(qū)內(nèi)的一種推擠型隆起（鄧起東等，1973；鄧起東和尤惠川，1985），故我們推測此處綜合震源機制解的復雜性由此處的構(gòu)造復雜性所引起，或者說是此處地震類型復雜，有左旋正傾滑和逆傾滑地震，而應(yīng)力場計算是我們假定其具有均勻應(yīng)力場，故所得綜合震源機制解為不確定型，即為復雜應(yīng)力場的一個平均，霍州附近區(qū)域的應(yīng)力場還有待進一步深入研究.

3.5 渭河斷陷帶內(nèi)應(yīng)力場結(jié)果及分析

渭河斷陷帶位于鄂爾多斯塊體南緣、秦嶺斷塊山地以北，由渭河、運城、靈寶三個斷陷盆地及其他構(gòu)造單元如凸起和小型斷塊山地組成，總體呈近東西向分布，向東逐漸轉(zhuǎn)為北東東向，東北部與山西斷陷帶為鄰，西端與鄂爾多斯西南邊界弧形斷裂束相接.盆地兩側(cè)均為正斷層所控制，盆地為左旋剪切拉張帶（鄧起東和尤惠川，1985；胡惠民，1987；國家地震局《鄂爾多斯周緣活動斷裂系》課題組，1988）.由圖3和圖4可見，P、T軸總體方位分布和先前研究結(jié)果較為一致（薛宏運和鄢家全，1984；許忠淮等，1989；范俊喜等，2003；謝富仁等；2004；徐黎明等，2006）.T軸方位大體上呈北北西，大體與渭河盆地內(nèi)斷陷盆地走向近垂直，傾角近水平；P軸傾角較大，反映了渭河斷陷帶現(xiàn)今處于拉張型應(yīng)力狀態(tài)，與盆地構(gòu)造相吻合.在圖5中，綜合震源機制解結(jié)果在渭河斷陷帶內(nèi)均為正斷層型為主，部分綜合震源機制解有少量的走滑分量，節(jié)面總體分布與盆地走向相吻合.從綜合震源機制解類型來看，本區(qū)域內(nèi)現(xiàn)今應(yīng)力狀態(tài)以拉張為主.

3.6 鄂爾多斯塊體西南邊緣應(yīng)力場結(jié)果及分析

鄂爾多斯塊體的西南邊緣（從寧夏南部中寧、中衛(wèi)一帶至陜西寶雞）是一條特殊性質(zhì)的邊界，它沒有發(fā)育拉張性質(zhì)為主的地塹式盆地，而是以一組向北突出的弧形斷裂為特色.這些弧形斷裂第四紀以來以強烈的擠壓性質(zhì)及左旋走滑運動區(qū)別與鄂爾多斯邊緣的其他活動斷裂（鄧起東和尤惠川，1985；胡惠民，1987；國家地震局《鄂爾多斯周緣活動斷裂系》課題組，1988）.由圖3和圖4可見，該區(qū)域的P軸方位大體上以北東至北北東向為主，傾角較為水平；T軸呈北西至北西西向，傾角較小，但稍大于P軸傾角，P、T軸方位大體上和先前結(jié)果較為一致（薛宏運和鄢家全，1984；許忠淮等，1989；杜興信和邵輝成，1999；范俊喜等，2003；謝富仁等；2004；張輝等，2012）.從P、T軸傾角均較小，T軸傾角略大于P軸來看，該區(qū)域的應(yīng)力類型應(yīng)為剪切兼擠壓型.在圖5中，該區(qū)域的綜合震源機制解類型較為復雜，用上述震源機制解的分類方法所分的幾種類型在該區(qū)域均有，在海原—景泰斷裂至牛首山斷裂間，綜合震源機制解類型大多為正斷層型，與先前認為該區(qū)域的逆沖兼有左旋走滑性質(zhì)不一致，具體原因還有待于進一步研究.從海原向南綜合震源機制解類型變?yōu)槟鏇_兼走滑型，逐漸過渡到走滑型，但細看震源機制解，我們不難發(fā)現(xiàn)被分為不確定型的和走滑型的綜合震源機制解中都含有一定的逆沖分量，在一定程度上反映了該區(qū)域的應(yīng)力場中均存在擠壓分量.至于用上述震源機制解分類方法，分得的結(jié)果致使部分綜合震源機制解沒有顯示出其所含的逆沖性，可能是所用分類方法的區(qū)分度有限，只能歸結(jié)為上述的6種類型.由海原至隴縣的綜合震源機制解均具有逆沖和走滑分量，可見該區(qū)域的應(yīng)力場應(yīng)為擠壓走滑型.本區(qū)域的研究結(jié)果與先前的認識大體上一致，但逆沖強度較先前認識的要弱，具體原因仍需進一步研究.

3.7 鄂爾多斯塊體鄰近地區(qū)應(yīng)力場結(jié)果及分析

由圖3和圖4可見，在甘東南地區(qū)，P軸走向從蘭州向南至文縣，逐漸由北北東向轉(zhuǎn)為近東西向，T軸走向由北西西轉(zhuǎn)為近南北向，和卜玉菲等（2013）利用甘東南寬頻帶流動臺陣以及周邊固定臺站P波初動資料所得結(jié)果以及張輝等（2007）利用震源機制解所得應(yīng)力場結(jié)果大體上一致.從圖5看，該區(qū)域的綜合震源機制解結(jié)果較為復雜，在漳縣和岷縣一帶綜合震源機制解為正斷走滑型，P、T軸在方位上與卜玉菲等（2013）接近，本研究中給出的P軸傾角大于T軸傾角，與卜玉菲等（2013）的關(guān)于P、T軸傾角結(jié)果剛好相反，兩個研究所用數(shù)據(jù)量大小相近，但本研究結(jié)果此處矛盾比相對較高，此處應(yīng)力場結(jié)果還有待于進一步探討.

從河北石家莊向南至河南許昌地區(qū)，P軸方位由北向南，從西向東，逐漸由北東東向轉(zhuǎn)為近東西向，P軸的傾角從西向東，傾角逐漸變得更加水平，表明由鄂爾多斯周緣的拉張型應(yīng)力場向外逐漸過渡.T軸方位呈北北西至南北向，隨P軸一同呈現(xiàn)出規(guī)律性變化.上述P、T軸方位分布結(jié)果和已有研究結(jié)果較為一致（許忠淮等，1989；杜興信和邵輝成，1999；徐紀人和趙志新，2006）.

4 鄂爾多斯塊體周緣斷陷盆地成因推測

先前許多學者認為鄂爾多斯塊體存在旋轉(zhuǎn)（蘇剛，1984；徐錫偉等，1994；邵輝成和蘇剛，1999；張希等，2002；陳小斌等，2005）；也有學著認為鄂爾多斯塊體不存在旋轉(zhuǎn)（鄧起東和尤惠川，1985；鄧起東等，1999；許忠淮等，2000），因為在其東北緣地區(qū)和外圍是相連的，沒有環(huán)向的走滑斷層活動.且不論是順時針旋轉(zhuǎn)還是逆時針旋轉(zhuǎn)，和其東西兩側(cè)的右旋剪切拉張帶以及南北兩側(cè)的左旋剪切拉張帶性質(zhì)相矛盾；若鄂爾多斯塊體是相對堅硬且旋轉(zhuǎn)的塊體，則其周緣的運動速度或是地震活動性應(yīng)大體一致，但其西北角地震活動性明顯低于其他區(qū)域.我們的研究結(jié)果表明，在環(huán)繞鄂爾多斯周緣的銀川—吉蘭泰斷陷帶、河套斷陷帶、岱海斷陷帶、山西斷陷帶和渭河斷陷帶內(nèi)，綜合震源機制解結(jié)果均以正斷層型為主，且綜合震源機制解節(jié)面走向的變化與控制斷陷帶邊界的主要斷裂走向相一致，拉張軸近似垂直于斷層走向，與鄂爾多斯周緣斷陷帶現(xiàn)今的拉張狀態(tài)相一致.為了合理解釋這一現(xiàn)象，我們由鄧起東和尤惠川（1985）的研究入手，他們利用帶有剛性橢圓核的無限平板在遠場壓應(yīng)力作用下的平面應(yīng)力問題分析了鄂爾多斯周緣的構(gòu)造活動及應(yīng)力場，得到鄂爾多斯周緣地區(qū)可以分為四個區(qū)域，其中西北和東南兩區(qū)為凹陷區(qū)，東北和西南兩個區(qū)為擠壓隆起區(qū)，上述研究結(jié)果大體上和我們所獲得的應(yīng)力場結(jié)果相一致，但在東北區(qū)我們的應(yīng)力場結(jié)果為拉張型，與之相矛盾.鄂爾多斯東北區(qū)也為張性的盆地為主，與鄧起東和尤惠川（1985）的擠壓隆起區(qū)結(jié)果相矛盾，而和本研究結(jié)果較為符合.鄂爾多斯東北緣地區(qū)存在多個相互平行的拉張盆地，我國著名的大同第四紀火山群分布在這一區(qū)域，火山巖包裹體研究揭示，在第四紀時期，該地區(qū)上地幔仍然繼續(xù)其底辟上涌過程（陳孝德等，1997；陳孝德等，2001）.地震成像研究揭示，該地區(qū)在上地幔范圍內(nèi)存在明顯的低速異常（何正勤等，2009）.綜合上述研究結(jié)果，我們推測在鄂爾多斯東北緣地區(qū)，盆地下方深部物質(zhì)的上涌形成的拉張作用大于鄧起東和尤惠川（1985）推導出的擠壓作用，故表現(xiàn)為拉張性.鄂爾多斯塊體在其西南角受到來自青藏高原北東向的擠壓力作用，其余區(qū)域均受到北西—南東向的水平拉張力作用，由于其自身是一個較為穩(wěn)定的塊體，故而在其周緣形成了一系列的斷陷盆地.在鄂爾多斯塊體周緣及附近區(qū)域，主應(yīng)變場為NEE－SWW方向壓縮與NNW－SSE方向擴張；主張應(yīng)變大于主壓應(yīng)變，并推測其原因很可能是與華北地殼的底部受到地幔物質(zhì)上隆作用而產(chǎn)生的NW－SE方向的擴張力有關(guān)（張靜華等，2004）.沈正康等（2003）給出的鄂爾多斯周緣地區(qū)擴容率相對較大，鄂爾多斯周緣的應(yīng)變場以及擴容率與鄂爾多斯塊體處于北北西—南南東向的水平拉張力作用為主的應(yīng)力場相符.至此，鄂爾多斯塊體在西南角受到來自青藏高原的擠壓和其東北角深部物質(zhì)上涌形成的拉張力的聯(lián)合作用，上述作用使得鄂爾多斯周緣地區(qū)除西南區(qū)為擠壓區(qū)外，其余區(qū)域均為剪切拉張區(qū)得到了合理的解釋，并與先前研究認為鄂爾多斯周緣地區(qū)處于引張應(yīng)力場作用相符合（高名修，1979；成爾林和張美芳，1990；劉巍等，1995），較好地解釋了環(huán)鄂爾多斯周緣的斷陷盆地構(gòu)造，亦符合鄂爾多斯塊體東西兩側(cè)的右旋剪切拉張帶以及南北兩側(cè)的左旋剪切拉張帶.

5 結(jié)論

我們將綜合震源機制解法應(yīng)用于鄂爾多斯周緣地殼應(yīng)力場的計算，使用了2007年8月1日至2013年7月21日鄂爾多斯塊體周緣發(fā)生的8499個地震的49844個P波初動符號資料，所得的0.25°×0.25°應(yīng)力場結(jié)果基本上覆蓋了鄂爾多斯周緣地區(qū).

根據(jù)鄂爾多斯塊體周緣地區(qū)應(yīng)力場的初步結(jié)果，我們得到鄂爾多斯塊體周緣地殼應(yīng)力場具有以下特征：

（1）在環(huán)繞鄂爾多斯周緣的銀川—吉蘭泰斷陷帶、河套斷陷帶、岱海斷陷帶、山西斷陷帶和渭河斷陷帶內(nèi)，綜合震源機制解結(jié)果以正斷層型為主，且綜合震源機制解節(jié)面走向大體與控制斷陷帶邊界的主要斷裂走向相一致，與鄂爾多斯周緣斷陷帶現(xiàn)今的拉張狀態(tài)相一致.

（2）在鄂爾多斯西南緣，綜合震源機制解類型主要為逆沖、逆沖走滑和走滑型，反映了鄂爾多斯塊體在西南緣受到青藏高原北東向的擠壓作用.鄂爾多斯西南緣的應(yīng)力場的主壓應(yīng)力方向在遠處為東向，源自于青藏高原向東北擠壓的作用，靠近鄂爾多斯塊體表現(xiàn)為北東—南西向.

（3）P軸方位在局部地區(qū)變化較大，但總體呈現(xiàn)規(guī)律性變化.P軸方位在鄂爾多斯塊體西緣，從南向北，主壓應(yīng)力軸方位更加偏北；在其北緣，由西向東，主壓應(yīng)力軸方位更加偏東.在其南緣和東緣，主壓應(yīng)力軸方位變化不大，大體上平行于控制各斷陷帶主要斷裂的走向.P軸傾角在西南緣為近水平，在其周緣各盆地內(nèi)P軸傾角近直立.

（4）T軸方位總體表現(xiàn)為北西—南東向；在鄂爾多斯周緣的斷陷帶內(nèi)，T軸走向大體與控制斷陷帶主要斷裂走向以及斷陷盆地走向相垂直.

（5）鄂爾多斯塊體在其西南角受到來自青藏高原的北東向擠壓和其東北角深部物質(zhì)上涌形成的北西—南東向拉張力聯(lián)合作用，上述作用使得鄂爾多斯周緣地區(qū)除西南區(qū)為擠壓區(qū)外，其余區(qū)域均為剪切拉張區(qū)，與先前研究認為鄂爾多斯周緣地區(qū)處于引張應(yīng)力場作用相符合，較好地解釋了環(huán)鄂爾多斯周緣的斷陷盆地構(gòu)造，與鄂爾多斯塊體東西兩側(cè)的右旋剪切拉張帶以及南北兩側(cè)的左旋剪切拉張帶相符合.

6 地震目錄誤差對結(jié)果的影響與結(jié)果討論

圖6 加入隨機誤差后應(yīng)力場計算結(jié)果施密特圖實心圓、三角形和正方形分別表示未加入誤差時綜合震源機制解的P、B和T軸方向；空心圖形表示加入隨機誤差后的計算結(jié)果.圖（a）中震中位置和震源深度誤差均為［－5km，5km］之間，圖（b）中震中位置誤差同圖（a），震源深度誤差為［－10km，10km］.Fig.6 Schmidt projection of the stress field calculation results after adding random errorCircle，triangle and square indicate the direction of P，Band Taxis of the composite focal mechanism by using the data without random error，respectively.The corresponding hollow symbols indicate the calculation results of the data with random error.The error of hypocenter and focal depth is between－5km and 5km in figure（a），figure（b）with the same hypocenter error as figure（a），the error of focal depth is between－10km and 10km in figure（b）.

地震目錄中事件參數(shù)的測定誤差總是不可避免的，震中位置及震源深度測定誤差會影響到P波初動符號在震源球上的位置，這將會對計算結(jié)果產(chǎn)生何種影響及其影響程度如何？對此我們采用類似蔣長勝等（2013）的方法，對地震目錄中的地震位置參數(shù)加入隨機誤差的方式來考察研究結(jié)果的穩(wěn)健性.我們選擇了研究區(qū)域中經(jīng)、緯度分別為112.25°、40.25°的網(wǎng)格點進行測試，由圖3和圖4可見，該網(wǎng)格點的P波初動符號數(shù)和矛盾比均處于中等水平，故我們選擇該網(wǎng)格點進行測試.我們對該網(wǎng)格點計算應(yīng)力場時所用的地震震中位置和震源深度均加入隨機誤差，再根據(jù)加入隨機誤差后的數(shù)據(jù)計算應(yīng)力場，反復進行1000次試算以考察事件參數(shù)誤差對應(yīng)力場計算結(jié)果的影響.第一次試驗震中位置和震源深度加入的隨機誤差均為［－5km，5km］之間，所得結(jié)果見圖6a；第二次試驗震中位置誤差同第一次試驗，震源深度誤差為［－10km，10km］之間，所得結(jié)果見圖6b.由圖6可見，加入隨機誤差后應(yīng)力場計算結(jié)果受到了一定程度的影響，且深度誤差增大后，計算結(jié)果離散程度增大，但整體來看計算結(jié)果還是穩(wěn)健的，反映了雖然地震目錄中事件參數(shù)的不確定性會影響到應(yīng)力場計算結(jié)果，但當數(shù)據(jù)量較大時，應(yīng)力場計算結(jié)果受到地震事件參數(shù)測定誤差的影響較小.這一現(xiàn)象的原因可能是由于綜合震源機制解法應(yīng)用了大量的P波初動數(shù)據(jù)，大量數(shù)據(jù)間的隨機誤差相互抑制，因而所得應(yīng)力場計算結(jié)果相對較穩(wěn)定.

在本研究中地震參數(shù)我們直接使用了中國地震臺網(wǎng)給出的地震目錄，沒有對其進行進一步的精確定位，一方面因為在對地震進行精定位時，對地震觀測要求較高，如地震事件須被4個以上臺站記錄到，這樣會導致部分小地震事件不滿足精定位要求而被舍棄，其初動符號也就不能夠被用于應(yīng)力場計算，這將會使可用初動符號數(shù)大大減少；另一方面由上述試驗結(jié)果可見，在P波初動數(shù)據(jù)量較大時，地震事件的參數(shù)對應(yīng)力場計算結(jié)果的影響較小，故沒有對地震事件進行精定位.地震參數(shù)對于應(yīng)力場計算結(jié)果終究有一定程度的影響，且當參與計算的數(shù)據(jù)量較少時，影響可能會增大，因而在后續(xù)研究中應(yīng)予以重視.鄂爾多斯塊體周緣主應(yīng)力方向分布情況，和先前研究結(jié)果總體上較為一致（薛宏運和鄢家全，1984；許忠淮等，1989；范俊喜等，2003；謝富仁等；2004；徐黎明等，2006）.同時，所給出的精細應(yīng)力場中也出現(xiàn)了和先前研究不一致的地方，如鄂爾多斯塊體西南緣的逆沖程度不如先前認為的強烈等，有個別地區(qū)出現(xiàn)了應(yīng)力場方向的突然變化的現(xiàn)象，這種現(xiàn)象體現(xiàn)了應(yīng)力場的不均勻性還是由于其他原因所致，若是應(yīng)力場不均勻性的體現(xiàn)，其和局部地質(zhì)構(gòu)造的關(guān)系又是什么？上述現(xiàn)象的解釋和探究將是我們下一步的研究工作.

致謝中國地震局地球物理研究所“國家數(shù)字測震臺網(wǎng)數(shù)據(jù)備份中心”為本研究提供地震波形數(shù)據(jù)；中國地震局地球物理研究所蔣長勝研究員給予的幫助和指導；文中部分圖件利用Generic Mapping Tools（GMT）（Wessel ＆Smith，1 9 9 5）繪制的，審稿人的建議對本文質(zhì)量的提高有很大幫助，在此一并致謝.

Aki K.1966.Earthquake generating stress in Japan for the years 1961to 1963obtained by smoothing the first motion radiation patterns.Bull.Earthq.Res.Inst.，44（2）：447－471.

An M J，Li F Q.1998.Current tectonic stress field of Shanxi garben system.ActaSeismologicaSinica（in Chinese），20（5）：461－465.

Bu Y F，Zhang Y S，Wan Y G，et al.2013.The tectonic stress field in Southeastern area of Gansu province deduced from P wave polarity data.ChinaEarthquakeEngineeringJournal（in Chinese），35（1）：160－165.

Chen A L，Zhao S L，Lu D H，et al.1981.The characteristics of tectonic stress field and its relations to seismicity in Xiji－Haiyuan－Guyuan region.NorthwesternSeismologicalJournal（in Chinese），3（2）：31－38.

Chen X B，Zang S X，Liu Y G，et al.Horizontal movement of Ordos block and the interaction of Ordos block and adjacent blocks.JournaloftheGraduateSchooloftheChineseAcademyof Sciences（in Chinese），22（3）：309－314.

Chen X D，Lin C Y，Zhang X O，et al.1997.Deformation features of mantle xenoliths from Quaternary basalts in Datong，Shanxi province and their rheological implications.Seismologyand Geology（in Chinese），19（4）：313－320.

Chen X D，Shi L B，Lin C Y.2001.A research on Quaternary volcanism in north China.SeismologyandGeology（in Chinese），23（4）：564－573.

Cheng E L，Li G F，Chen H C.1982.On the direction of the maximum compressive principal stress before and after the 1976 Songpan－Pingwu earthquake（M＝7.2）of the Sichuan province.ActaSeismologicaSinica（in Chinese），4（2）：136－148.

Cheng E L，Zhang M F.1990.Modern NW－NNW tension stress field of rift valley series in Shanxi province.Northwestern SeismologicalJournal（in Chinese），12（4）：20－25.

Deng Q D，Wang K L，Wang Y P，et al.1973.On the tendency of seismicity and their geological set－up of the seismic belt of Shanxi garben.ScientiaGaologicaSinica（in Chinese），（1）：37－47.

Deng Q D，You H C.1985.Mechanism and structural activity of Ordos basin.／／Study on the Modern Crustal Movement（1）（in Chinese）.Beijing：Seismological Press，58－78.

Deng Q D，Cheng S P，Min W，et al.1999.Discussion on Cenozoic tectonics and dynamics of Ordos block.JournalofGeomechanics（in Chinese），5（3）：13－21.

Du X X，Shao H C.1999.Modern tectonic stress field in the Chinese mainland inverted from focal mechanism solutions.Acta SeismologicaSinica（in Chinese），21（4）：354－360.

Fan J X，Ma J，Diao G L.2003.Contemporary tectonic stress field around the Ordos fault block inferred from earthquake focal mechanisms.SeismologyandGeology（in Chinese），25（1）：88－99.

Fan T Y，Long C X，Yang Z Y，et al.2012.Comprehensive modeling on the present crustal stress of China mainland with the viscoelastic spherical shell.ChineseJ.Geophys.（in Chinese），55（4）：1249－1260，doi：10.6038／j.issn.00015733.2012.04.020.

Gao M X.1979.Current tensile stress field in the block－faulting tectonic areas of North China.SeismologyandGeology（in Chinese），1（2）：1－12.

He Z Q，Ye T L，Ding Z F.2009.Surface wave tomography for the phase velocity in the northeastern part of North China.Chinese J.Geophys.（in Chinese），52（5）：1233－1242，doi：10.3969／j.issn.0001－5733.2009.05.012.

Hu H M.1987.The characteristic of recent tectonic activity in Shanxi graben system.EarthquakeResearchinChina（in Chinese），3（4）：67－73.

Hu H M.1989.The modern tectonic activity and regional stress field of the Ordos block and the graben system around it.CollectedWorksofTectonicandCrustalStress（in Chinese），3：60－68.

Hua X W.1980.Stress variations in the region around Beijing and Tianjin before and after the 1976Tangshan earthquake.Acta SeismologicaSinica（in Chinese），2（2）：130－146.

Jiang C S，Wu Z L，Zhuang J C.2013.ETAS model applied to the earthquake－sequence association（ESA）problem：the Tangshan sequence.ChineseJ.Geophys.（in Chinese），56（9）：2971－2981，doi：10.6038／cj920130911.

Jiang J L.1991.Recent tectonic stress field in Shanxi province and its adjacent area.NorthwesternSeismologicalJournal（in Chinese），13（4）：85－88.

Kan R J，Liu Z Y，Wang S J，et al.1983.The temporal variation of the orientation of the mean stress axis both before and after the strong Longling earthquake.ChineseJ.Geophys.（in Chinese），26（6）：564－573.

Li M L.1985.The seismic activity and tectonic around Ordos block.NorthwesternSeismologicalJournal（in Chinese），7（Suppl.）：1－9.

Li M L，Zhao Z J.1981.Focal mechanism and stress field of Ningxia region.NorthwesternSeismologicalJournal（in Chinese），3（4）：49－55.

Li M L，Zhao Z J.1986.Modern tectonic stress field and its relation with seismicities in Ningxia and its neighborhood.Journalof SeismologicalResearch（in Chinese），9（3）：299－314.

Li Q Z.1980.General features of the stress field in the crust of North China.ChineseJ.Geophys.（in Chinese），23（3）：376－388.

Li Q Z，Wang Z G，Jia Y N，et al.1973.Stress field obtained for two regions from weak earthquake data recorded at a single seismic station.ChineseJ.Geophys.（in Chinese），16（1）：49－61.

Liang H Q，Yang Z X，Guo Q L.1986.Characteristics of present tectonic stress field of Shanxi graben system.Earthquake（in Chinese），（1）：47－52.

Liu W，Zhao X P，An W P，et al.1993.The crustal stress field in Shanxi province.EarthquakeResearchinShanxi（in Chinese），（3）：3－11.

Liu W，Zhao X P，An W P.1994.Characteristics of stress field in Taiyuan basin.EarthquakeResearchinShanxi（in Chinese），（1）：18－24.

Liu W，Meng Y Y，Wang Z L.1995.Present tectogenesis and characteristics of crustal stress field in Datong basin.Earthquake ResearchinShanxi（in Chinese），（1）：7－13.

Liu W，An W P，Zhao X P.1996.Recent tectonic movement and recent crust stress field in Yuncheng basin.EarthquakeResearchin Shanxi（in Chinese），（2）：9－13.

Northrup C J，Royden L H，Burchfiel B C.1995.Motion of the Pacific plate relative to Eurasia and its potential relation to Cenozoic extension along the eastern margin of Eurasia.Geology，23（8）：719－722.

Peltzer G，Tapponnier P，Zhang Z T，et al.1985.Neogene and Quaternary faulting in and along the Qinling Shan.Nature，317（6037）：500－505.

Shao H C，Su G.1999.Analysis on the seismicity trend around the Ordos block.NorthwesternSeismologicalJournal（in Chinese），21（4）：395－398.

Shen Z K，Jackson D D，Ge B X.1996.Crustal deformation across and beyond the Los Angeles basin from geodetic measurements.J.Geophys.Res.，101（B12）：27957－27980.

Shen Z K，Wang M，Gan W J，et al.2003.Contemporary tectonic strain rate field of Chinese continent and its geodynamic implications.EarthScienceFrontiers（ChinaUniversityof Geosciences，Beijing）（in Chinese），10（Suppl.）：93－100.

Sheng S Z，Wan Y G，Li H G，et al.2013.Direction of tectonic stress in Nantong city deduced from composite focal mechanism method.ProgressinGeophysics（in Chinese），28（3）：1297－1303，doi：10.6038／pg20130322.

Su G.1984.The research of regional seismicity using the movingblock as an unit （The Orduoce block and its seismicity）.NorthwesternSeismologicalJournal（in Chinese），6（2）：1－9.Tapponnier P，Peltzer G，Le Dain A Y，et al.1982.Propagating extrusion tectonics in Asia：New insights from simple experiments with plasticine.Geology，10（12）：611－616.

Tapponnier P，Peltzer G，Armijo R.1986.On the mechanics of the collision between India and Asia.GeologicalSociety，London，SpecialPublications，19（1）：113－157.

The Research Group on Active Fault System around Ordos Massif，SSB.1988.Active Fault System Around Ordos Massif （in Chinese）.Beijing：Seismological Press.

Tian Z Y，Han P，Xu K D.1992.The Mesozoic－Cenozoic East China rift system.Tectonophysics，208（1－3）：341－363.

Uyeda S，Kanamori H.1979.Back－arc opening and the mode of subduction.J.Geophys.Res.，84（B3）：1049－1061.

Wan Y G.2010.Contemporary tectonic stress field in China.Earthq.Sci.，23（4）：377－386，doi：10.1007／s11589－010－0735－5.

Wan Y G，Wu Y M，Sheng S Z，et al.2011.Priliminary result of Taiwan 3－D stress field from P wave polarity data.ChineseJ.Geophys.（in Chinese），54（11）：2809－2818，doi：10.3969／j.issn.0001－5733.2011.11.011.

Wang S Y，Xu Z H.1985.Seismo－tectonic stress field in east China.ActaSeismologicaSinica（in Chinese），7（1）：17－32.

Wang S Y，Xu Z H，Yu Y X，et al.1996.Inversion for the plate driving forces acting at the boundaries of China and its surroundings.ChineseJ.Geophys.（in Chinese），39（6）：764－771.

Wessel P，Smith W H F.1995.New version of the generic mapping tools.Eos，TransactionsAmericanGeophysicalUnion，76（33）：329－329.

Xie F R，Cui X F，Zhao J T，et al.2004.Regional division of the recent tectonic stress field in China and adjacent areas.Chinese J.Geophys.（in Chinese），47（4）：654－662，doi：10.3321／j.issn：0001－5733.2004.04.016.

Xie F R，Shu S B，Dou S Q，et al.2000.Quaternary tectonic stress field in the region of Haiyuan－Liupanshan fault zone to Yinchuan fault－depression.SeismologyandGeology（in Chinese），22（2）：139－146.

Xu J R，Zhao Z X.2006.Regional characteristics of the lithospheric stress field and tectonic motions in China and its adjacent areas.GeologyinChina（in Chinese），33（4）：782－792.

Xu J R，Zhao Z X，Ishikawa Y.2008.Regional characteristics of crustal stress field and tectonic motions in and around Chinese mainland.ChineseJ.Geophys.（in Chinese），51（3）：770－781，doi：10.3321／j.issn：0001－5733.2008.03.018.

Xu L M，Zhou L F，Zhang Y K，et al.2006.Characteristics and tectonic setting of tectono－stress field of Ordos basin.GeotectonicaetMetallogenia（in Chinese），30（4）：455－462.

Xu M J，Jiang Y J，Zhou C Y.1996.Preliminary analysis for the space distribution of present stress field in Jiangsu and its adjacent area.EarthquakeResearchinChina（in Chinese），12（4）：383－388.

Xu X W，Cheng G L，Ma X Y，et al.1994.Rotation model and dynamics of blocks in North China and its adjacent areas.Earth Science－JournalofChinaUniversityofGeosciences（in Chinese），19（2）：129－138.

Xu Z H，Liu Y F，Zhang Y Z，et al.1979.On the characteristic of direction of the earthquake stress field around the Beijing area.ActaSeismologicaSinica（in Chinese），1（2）：121－132.

Xu Z H，Yan M，Zhao Z H.1983.Evaluation of the direction of tectonic stress in North China from recorded data of a large number of small earthquakes.ActaSeismologicaSinica（in Chinese），5（3）：268－279.

Xu Z H，Wang S Y，Huang Y R，et al.1987.Direction of mean stress axes in southwestern China deduced from microearthquake data.ChineseJ.Geophys.（in Chinese），30（5）：476－486.

Xu Z H，Wang S Y，Huang Y R，et al.1989.The tectonic stress field of Chinese continent deduced from a great number of earthquakes.ChineseJ.Geophys.（in Chinese），32（6）：636－647.

Xu Z H，Wang S Y，Gao A J.2000.Present－day tectonic movement in the northeastern margin of the Qinghai－Xizang （Tibetan）plateau as revealed by earthquake activity.ActaSeismologica Sinica（in Chinese），22（5）：472－481.

Xu Z H.2001.A present－day tectonic stress map for eastern Asia region.ActaSeismologicaSinica（in Chinese），23（5）：492－501.

Xue H Y，Yan J Q.1984.The contemporary crustal stress field around the Ordos block.ChineseJ.Geophys.（in Chinese），27（2）：144－152.

Yan J Q，Shi Z L，Wang S Y，et al.1979.Some features of the recent tectonic stress field of China and environs.Acta SeismologicaSinica（in Chinese），1（1）：9－24.

Zhang D N，Gao L S.1989.Three dimensional numerical simulation of eastern Asian stress field.EarthqukeResearchinChina（in Chinese），5（4）：24－33.

Zhang H，Liu X F，Zhao L Y.2007.Characteristics analysis of the stress field in Gansu and its adjacent region with various focal mechanism solutions.SouthChinaJournalofSeismology（in Chinese），27（2）：33－39.

Zhang H，Gao Y，Shi Y T，et al.2012.Tectonic stress analysis based on the crustal seismic anisotropy in the northeastern margin of Tibetan plateau.ChineseJ.Geophys.（in Chinese），55（1）：95－104，doi：10.6038／j.issn.0001－5733.2012.01.009.Zhang J H，Li Y X，Guo L Q，et al.2004.Study on present－day deformation and strain field in North China by use of GPS data.JournalofGeodesyandGeodynamics（in Chinese），24（3）：40－46.

Zhang X，Jiang Z S，Zhang X L，et al.2002.A seismic negative dislocation inversion of recent horizontal crust movement in North China.JournalofGeodesyandGeodynamics（in Chinese），22（3）：40－45.

Zhang Y Q，Liao C Z，Shi W，et al.2006.Neotectonic evolution of the peripheral zones of the Ordos basin and geodynamic setting.GeologicalJournalofChinaUniversities（in Chinese），12（3）：285－297.

Zhao Z J，Liu X J.1990.Seismic activity and local tectonic stress field in Ningxia and nearby regions.SeismologyandGeology（in Chinese），12（1）：31－46.

Zhao Z J，Liu X J，Kang L Y.2002.Features of focal－mechanism solutions in Ningxia and its neighboring region.Northwestern SeismologicalJournal（in Chinese），24（2）：162－166，173.

Zheng X F，Ouyang B，Zhang D N，et al.2009.Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake.ChineseJ.Geophys.（in Chinese），52（5）：1412－1417，doi：10.3969／j.issn.0001－5733.2009.05.031.

Zhu S B，Shi Y L.2006.Study on the cause of tectonic stress field of Chinese mainland and its adjacent area.ScienceinChina SeriesD（EarthScience）（in Chinese），36（12）：1077－1083.

Zoback M L.1992.First－and second－order patterns of stress in the lithosphere：the world stress map project.J.Geophys.Res.，97（B8）：11703－11728.

附中文參考文獻

安美建，李方全.1998.山西地塹系現(xiàn)今構(gòu)造應(yīng)力場.地震學報，20（5）：461－465.

卜玉菲，張元生，萬永革等.2013.P波極性揭示的甘東南地區(qū)構(gòu)造應(yīng)力場特征.地震工程學報，35（1）：160－165.

陳愛玲，趙淑蘭，呂德徽等.1981.西海固地區(qū)構(gòu)造應(yīng)力場特征及其與地震活動性的關(guān)系.西北地震學報，3（2）：31－38.

陳小斌，臧紹先，劉永崗等.2005.鄂爾多斯地塊的現(xiàn)今水平運動狀態(tài)及其與周緣地塊的相互作用.中國科學院研究生院學報，22（3）：309－314.

陳孝德，林傳勇，張小鷗等.1997.山西大同第四紀火山巖中幔源包體的變形特征及其上地幔流變學意義.地震地質(zhì)，19（4）：313－320.

陳孝德，史蘭斌，林傳勇.2001.華北第四紀火山作用研究.地震地質(zhì)，23（4）：564－573.

成爾林，李桂芳，陳和川.1982.1976年四川省松潘—平武7.2地震前后主壓應(yīng)力軸的方向特征.地震學報，4（2）：136－148.

成爾林，張美芳.1990.山西地塹系現(xiàn)代北西—北北西向拉張應(yīng)力場.西北地震學報，12（4）：20－25.

鄧起東，王克魯，汪一鵬等.1973.山西隆起區(qū)斷陷地震帶地震地質(zhì)條件及地震發(fā)展趨勢概述.地質(zhì)科學，（1）：37－47.

鄧起東，尤惠川.1985.鄂爾多斯周緣斷陷盆地帶的構(gòu)造活動特征及其形成機制.／／現(xiàn)代地殼運動研究（1）.北京：地震出版社，58－78.

鄧起東，程紹平，閔偉等.1999.鄂爾多斯塊體新生代構(gòu)造活動和動力學的討論.地質(zhì)力學學報，5（3）：13－21.

杜興信，邵輝成.1999.由震源機制解反演中國大陸現(xiàn)代構(gòu)造應(yīng)力場.地震學報，21（4）：354－360.

范俊喜，馬瑾，刁桂苓.2003.由小震震源機制解得到的鄂爾多斯周邊構(gòu)造應(yīng)力場.地震地質(zhì)，25（1）：88－99.

范桃園，龍長興，楊振宇等.2012.中國大陸現(xiàn)今地應(yīng)力場黏彈性球殼數(shù)值模擬綜合研究.地球物理學報，55（4）：1249－1260，doi：10.6038／j.issn.00015733.2012.04.020.

高名修.1979.華北斷塊構(gòu)造區(qū)的現(xiàn)代引張應(yīng)力場.地震地質(zhì)，1（2）：1－12.

國家地震局《鄂爾多斯周緣活動斷裂系》課題組.1988.鄂爾多斯周緣活動斷裂系.北京：地震出版社.

何正勤，葉太蘭，丁志峰.2009.華北東北部的面波相速度層析成像研究.地球物理學報，52（5）：1233－1242，doi：10.3969／j.issn.0001－5733.2009.05.012.

胡惠民.1987.山西地塹系的現(xiàn)代構(gòu)造活動特征.中國地震，3（4）：67－73.

胡惠民.1989.鄂爾多斯及其周緣斷陷盆地現(xiàn)代構(gòu)造活動與區(qū)域應(yīng)力場.地殼構(gòu)造與地殼應(yīng)力文集，3：60－68.

華祥文.1980.唐山強震前后北京、天津周圍地區(qū)應(yīng)力的變化過程.地震學報，2（2）：130－146.

蔣長勝，吳忠良，莊建倉.2013.地震的“序列歸屬”問題與ETAS模型——以唐山序列為例.地球物理學報，56（9）：2971－2981，doi：10.6038／cj920130911.

姜家蘭.1991.陜西及其鄰近地區(qū)現(xiàn)代構(gòu)造應(yīng)力場.西北地震學報，13（4）：85－88.

闞榮舉，劉祖蔭，王紹晉等.1983.龍陵地區(qū)大地震前后平均應(yīng)力軸取向的時間變化.地球物理學報，26（6）：564－573.

李孟鑾.1985.鄂爾多斯塊體周圍的地震活動與構(gòu)造.西北地震學報，7（增刊）：1－9.

李孟鑾，趙知軍.1981.寧夏地區(qū)的震源機制與應(yīng)力場.西北地震學報，3（4）：49－55.

李孟鑾，趙知軍.1986.寧夏地區(qū)現(xiàn)代構(gòu)造應(yīng)力場及其與地震活動的關(guān)系.地震研究，9（3）：299－314.

李欽祖.1980.華北地殼應(yīng)力場的基本特征.地球物理學報，23（3）：376－388.

李欽祖，王澤皋，賈云年等.1973.由單臺小地震資料所得兩個區(qū)域的應(yīng)力場.地球物理學報，16（1）：49－61.

梁海慶，楊增學，郭啟良.1986.山西地塹系現(xiàn)今構(gòu)造應(yīng)力場特征.地震，（1）：47－52.

劉巍，趙新平，安衛(wèi)平等.1993.山西地區(qū)的地殼應(yīng)力場.山西地震，（3）：3－11.

劉巍，趙新平，安衛(wèi)平.1994.太原盆地的應(yīng)力場特征.山西地震，（1）：18－24.

劉巍，孟雁英，王趙麗.1995.大同盆地現(xiàn)今構(gòu)造活動及地殼應(yīng)力場特征.山西地震，（1）：7－13.

劉巍，安衛(wèi)平，趙新平.1996.運城盆地的現(xiàn)今構(gòu)造活動及現(xiàn)代地殼應(yīng)力場的基本特征.山西地震，（2）：9－13.

邵輝成，蘇剛.1999.鄂爾多斯周緣近期地震活動趨勢分析.西北地震學報，21（4）：395－398.

沈正康，王敏，甘衛(wèi)軍等.2003.中國大陸現(xiàn)今構(gòu)造應(yīng)變率場及其動力學成因研究.地學前緣 （中國地質(zhì)大學，北京），10（特刊）：93－100.

盛書中，萬永革，李紅光等.2013.應(yīng)用綜合震源機制解法推斷南通市地殼應(yīng)力場方向.地球物理學進展，28（3）：1297－1303，doi：10.6038／pg20130322.

蘇剛.1984.以運動地塊為單元的區(qū)域地震活動研究——鄂爾多斯地塊和它的地震活動.西北地震學報，6（2）：1－9.

萬永革，吳逸民，盛書中等.2011.P波極性數(shù)據(jù)所揭示的臺灣地區(qū)三維應(yīng)力結(jié)構(gòu)的初步結(jié)果.地球物理學報，54（11）：2809－2818，doi：10.3969／j.issn.0001－5733.2011.11.011.

汪素云，許忠淮.1985.中國東部大陸的地震構(gòu)造應(yīng)力場.地震學報，7（1）：17－32.

汪素云，許忠淮，俞言祥等.1996.中國及其鄰區(qū)周圍板塊作用力的研究.地球物理學報，39（6）：764－771.

謝富仁，崔效鋒，趙建濤等.2004.中國大陸及鄰區(qū)現(xiàn)代構(gòu)造應(yīng)力場分區(qū).地球物理學報，47（4）：654－662，doi：10.3321／j.issn：0001－5733.2004.04.016.

謝富仁，舒塞兵，竇素芹等.2000.海原，六盤山斷裂帶至銀川斷陷第四紀構(gòu)造應(yīng)力場分析.地震地質(zhì)，22（2）：139－146.

徐紀人，趙志新.2006.中國巖石圈應(yīng)力場與構(gòu)造運動區(qū)域特征.中國地質(zhì)，33（4）：782－792.

徐紀人，趙志新，石川有三.2008.中國大陸地殼應(yīng)力場與構(gòu)造運動區(qū)域特征研究.地球物理學報，51（3）：770－781，doi：10.3321／j.issn：0001－5733.2008.03.018.

徐黎明，周立發(fā)，張義楷等.2006.鄂爾多斯盆地構(gòu)造應(yīng)力場特征及其構(gòu)造背景.大地構(gòu)造與成礦學，30（4）：455－462.

徐鳴潔，姜永基，周翠英.1996.江蘇及鄰區(qū)現(xiàn)代應(yīng)力場空間分布特征分析.中國地震，12（4）：383－388.

徐錫偉，程國良，馬杏垣等.1994.華北及其鄰區(qū)塊體轉(zhuǎn)動模式和動力來源.地球科學——中國地質(zhì)大學學報，19（2）：129－138.

許忠淮，劉玉芬，張郢珍.1979.京、津、唐、張地區(qū)地震應(yīng)力場的方向特征.地震學報，1（2）：121－132.

許忠淮，閻明，趙仲和.1983.由多個小地震推斷的華北地區(qū)構(gòu)造應(yīng)力場的方向.地震學報，5（3）：268－279.

許忠淮，汪素云，黃雨蕊等.1987.由多個小震推斷的青甘和川滇地區(qū)地殼應(yīng)力場的方向特征.地球物理學報，30（5）：476－486.

許忠淮，汪素云，黃雨蕊等.1989.由大量的地震的資料推斷的我國大陸構(gòu)造應(yīng)力場.地球物理學報，32（6）：636－647.

許忠淮，汪素云，高阿甲.2000.地震活動反映的青藏高原東北地區(qū)現(xiàn)代構(gòu)造運動特征.地震學報，22（5）：472－481.

許忠淮.2001.東亞地區(qū)現(xiàn)今構(gòu)造應(yīng)力圖的編制.地震學報，23（5）：492－501.

薛宏運，鄢家全.1984.鄂爾多斯地塊周圍的現(xiàn)代地殼應(yīng)力場.地球物理學報，27（2）：144－152.

鄢家全，時振梁，汪素云等.1979.中國及鄰區(qū)現(xiàn)代構(gòu)造應(yīng)力場的區(qū)域特征.地震學報，1（1）：9－24.

張東寧，高龍生.1989.東亞地區(qū)應(yīng)力場的三維數(shù)值模擬.中國地震，5（4）：24－33.

張輝，劉小鳳，趙凌云.2007.由多個震源機制解分析甘肅及邊鄰地區(qū)應(yīng)力場特征.華南地震，27（2）：33－39.

張輝，高原，石玉濤等.2012.基于地殼介質(zhì)各向異性分析青藏高原東北緣構(gòu)造應(yīng)力場特征.地球物理學報，55（1）：95－104，doi：10.6038／j.issn.0001－5733.2012.01.009.

張靜華，李延興，郭良遷等.2004.用GPS測量結(jié)果研究華北現(xiàn)今構(gòu)造形變場.大地測量與地球動力學，24（3）：40－46.

張希，江在森，張曉亮等.2002.華北地區(qū)近期地殼水平運動的非震負位錯反演.大地測量與地球動力學，22（3）：40－45.

張岳橋，廖昌珍，施煒等.2006.鄂爾多斯盆地周邊地帶新構(gòu)造演化及其區(qū)域動力學背景.高校地質(zhì)學報，12（3）：285－297.

趙知軍，劉秀景.1990.寧夏及其鄰區(qū)地震活動帶與小區(qū)域構(gòu)造應(yīng)力場.地震地質(zhì)，12（1）：31－46.

趙知軍，劉秀景，康凌燕.2002.寧夏及鄰近地區(qū)震源機制解特征.西北地震學報，24（2）：162－166，173.

鄭秀芬，歐陽飚，張東寧等.2009.“國家數(shù)字測震臺網(wǎng)數(shù)據(jù)備份中心”技術(shù)系統(tǒng)建設(shè)及其對汶川大地震研究的數(shù)據(jù)支撐.地球物理學報，52（5）：1412－1417，doi：10.3969／j.issn.0001－5733.2009.05.031.

朱守彪，石耀霖.2006.中國大陸及鄰區(qū)構(gòu)造應(yīng)力場成因的研究.中國科學D輯 （地球科學），36（12）：1077－1083.