馬海勇　羅安湘　王朝陽　歐陽征健　馮娟萍

摘要：為研究弧形構造帶形成演化特征，以六盤山弧形構造帶早古生代晚期形成的弧形構造為例，采用野外面地質(zhì)調(diào)查、數(shù)值模擬等研究手段，對弧形構造帶形成規(guī)律進行區(qū)域塊體應力定量化的數(shù)值模擬分析。結果表明：六盤山弧形構造帶受到青藏地塊向北東向的推擠，北祁連地塊也向北東向位移并推擠北祁連弧后盆地，北祁連弧后盆地被相對堅硬的阿拉善地塊與鄂爾多斯地塊阻擋。早古生代晚期六盤山區(qū)域北東向擠壓的構造應力是早期弧形構造形成的主要原因，利用數(shù)值模擬方法構建的古生代晚期的數(shù)值模型，進行構造應力與地質(zhì)塊體形變規(guī)律數(shù)值模擬計算，使地質(zhì)模型處于平衡狀態(tài)后獲得x方向、y方向位移，在x方向的位移，計算結果顯示六盤山區(qū)域范圍內(nèi)古生代地質(zhì)塊體向東最大位移約為向北最大位移的25%，位移向量的位移方向總體表現(xiàn)為向北東弧形突起的構造形態(tài)。研究結果可為區(qū)域弧形構造演化研究提供定量化分析的參考和理論依據(jù)。

關鍵詞：弧形構造;六盤山;構造演化;數(shù)值模擬;位移向量

中圖分類號：TD 121文獻標志碼：A

文章編號：1672-9315（2022）04-0752-08

DOI：10.13800/j.cnki.xakjdxxb.2022.0415

Formation and evolution characteristics of

arcuate tectonic belt MA Haiyong LUO Anxiang WANG Chaoyang OUYANG Zhengjian FENG Juanping

（1.College of Geology and Environment，Xian University of Science and Technology，Xian 710054，China;

2.Research Institute of Petroleum Exploration and Development，Changqing Oilfield Company，Xian 710018，China;

3.Bailie School of Petroleum Engineering，Lanzhou City University，Lanzhou 730070，China）Abstract：In order to study the formation and evolution characteristics of arcuate tectonic belt，the arcuate structure formed in the Late Early Paleozoic in liupanshan arcuate tectonic belt as an example.By means of field profile geological survey and numerical simulation，an analysis was made of the formation rule of regional block stress in the initial stage of arc tectonic belt formation.The results show that the tectonic belt was pushed NE by the Qinghai-Tibet block，and the North Qilian blocks also moved NE and pushed the North Qilian back-arc basin，which was blocked by the relatively hard Alxa block and Ordos block.The tectonic stress action of the NE compression in the study area in the late Early Paleozoic is the main reason for the formation of the early arc-shaped structure.The numerical model of the tectonic belt in the late Paleozoic was constructed via the numerical simulation method to examine the tectonic stress and the deformation law of the geological block.After numerical simulation，the model was in equilibrium with the displacement in x direction and y direction? obtained.The displacement calculation results in x direction show that the maximum eastward displacement of the Paleozoic geological block in Liupanshan region is about 25% of the maximum northward displacement，and the displacement direction of the displacement vector is generally presented as the structural shape of arc-shaped protrusion in northeast direction.The results can provide a? theoretical guidance for quantitative analysis of regional arc tectonic evolution.

Key words：arcuate structure;Liupanshan;tectonic evolution;numerical simulation;displacement component

0引言

六盤山弧形構造帶位于青藏高原東北緣，處于阿拉善地塊、青藏地塊和鄂爾多斯地塊的交接部位（圖1），特殊的地理位置使其經(jīng)歷了多期次構造變形和復雜的構造演化[1-4]。張培震認為六盤山弧形構造帶受青藏地塊向北東向的推擠導致中、北祁連地塊也向北東向位移并推擠北祁連弧后盆地，北祁連弧后盆地被相對堅硬的阿拉善地塊與鄂爾多斯地塊阻擋的構造轉換帶[5]。王東等通過幾何學、運動學特征研究，認為六盤山弧形構造帶目前有7條近南北向的斷層，均呈現(xiàn)為北東向“弓”形態(tài)[6-8]。徐化超等通過動力學研究發(fā)現(xiàn)六盤山弧形構造帶僅新生代就經(jīng)歷了上新世晚期到第四紀早期、中更新世和晚更新世3個期次的構造變形和隆升，構造帶相對基底垂向抬升1 500 m[9-10]。近年來許多學者利用幾何學、地殼電性、磁性地層年代、GPS測量等方法對六盤山弧形構造帶弧形的形成時間和機制進行了大量的研究，并在中、新生代的弧形斷層走向、地殼結構及應力形變等方面取得了諸多成果[11-14]，但對于弧形構造形成初期古生代的演化過程研究相對較少[15-16]。向宏發(fā)等通過古地震方法分析六盤山區(qū)斷裂特征[17]，李天斌等通過運動學、動力學研究弧形推覆構造帶形成機理[18-19]，張進等通過寧南弧形構造動力分析認為早古生代晚期北東向擠壓是早期弧形構造形成的主要原因[20-21]，但是都缺乏地質(zhì)模型的定量數(shù)值分析的證據(jù)。張東寧等利用大陸巖石塊體動力學數(shù)值模型分析應力場、位移場變化[22-25]，因此，在對六盤山古生代弧形構造演化特征系統(tǒng)分析的基礎上，利用數(shù)值模擬方法構建了古生代晚期構造帶的數(shù)值模型，進行構造應力與地質(zhì)塊體形變規(guī)律研究，為弧形構造帶的大陸動力學定量解釋提供理論和實踐依據(jù)。

1區(qū)域地質(zhì)特征

六盤山弧形構造帶位于鄂爾多斯盆地西緣、青藏高原東北緣，夾持在鄂爾多斯地塊、阿拉善地塊、秦祁昆碰撞帶之間。六盤山弧形構造帶、其南西部位的秦祁昆碰撞帶、其北東部位的賀蘭山褶皺與其東部的鄂爾多斯盆地西緣沖斷帶都受到青藏地塊隆升的影響，地質(zhì)現(xiàn)象，構造復雜。六盤山弧形構造帶發(fā)育的斷層從北向南看，走向依次由北西逐漸變?yōu)槟媳保傮w呈現(xiàn)一個鏡像的“S”形狀，北部的斷層形態(tài)展布較為舒展，南部的斷層形態(tài)展布變得收束，弧形展布是該區(qū)的顯著特點。研究區(qū)弧形構造的基底是形成于中、晚元古代的秦—祁—賀三叉形裂谷的發(fā)育時期到志留紀、泥盆紀前陸盆地發(fā)育時期，蓋層經(jīng)歷了石炭紀-三疊紀、侏羅紀-白堊紀、新近紀等不同時代的發(fā)育，形成現(xiàn)今的構造形態(tài)。

六盤山弧形構造帶在早古生代晚期受到青藏地塊推擠而產(chǎn)生北祁連造山運動的力學作用。由于六盤山區(qū)域被阿拉善地塊、青藏地塊和鄂爾多斯地塊包圍，在研究區(qū)東北部形成銀川地塹;北祁連弧后盆地在中、北祁連地塊的推擠下逐漸閉合。東部過渡帶在逆沖推覆作用下形成許多逆沖推覆體，推覆體主要由下古生界的地層組成。隨著中、北祁連地塊繼續(xù)向北推進，六盤山弧形構造帶繼續(xù)向北東方向推移，并且其一部分緩慢匯入銀川地塹。對于六盤山弧形構造帶的地質(zhì)邊界條件是北部堅硬的阿拉善地塊和東部堅硬的鄂爾多斯地塊。西部同心等地區(qū)的下古生界地層厚度很大，東部中衛(wèi)、中寧、固原等地區(qū)下古生界地層相對較薄，區(qū)域整體東西方向的地層厚度存在較大變化。六盤山弧形構造帶的邊界形狀逐漸發(fā)生了變動，諸多原因導致整個六盤山弧形構造帶地層從上往下彎曲。這些現(xiàn)象正是早古生代晚期處于北祁連弧后盆地東部的六盤山弧形構造帶發(fā)育的幾個重要條件。六盤山弧形構造帶古生代形成了初期的弧形構造之后，經(jīng)歷了晚三疊世印支運動、早侏羅世燕山運動Ⅰ、中侏羅世燕山運動Ⅱ、晚侏羅世燕山運動Ⅲ、早白堊世燕山運動Ⅳ、晚白堊世燕山運動Ⅴ、古近紀早喜山運動與新近紀晚喜山運動等8次構造運動，構造帶主體向北推移[22-25]。

六盤山弧形構造帶主要發(fā)育3個斷裂帶，分別為海原-六盤山斷裂、中衛(wèi)-同心斷裂帶、牛首山-羅山斷裂帶。海原-六盤山斷裂帶在靖遠、固原和南華山向西南方向傾斜，傾角為45°～80°。斷層面上大量的近乎水平的條紋和臺階，這表明了斷層的反向正弦走滑運動。根據(jù)海原-六盤山斷裂帶發(fā)育的一系列拉脫盆地，可以確定該斷裂帶在中新世晚期表現(xiàn)為左旋走滑運動。中衛(wèi)-同心斷裂帶在金山地區(qū)向西南方向傾斜，傾角在63°～70°之間。中衛(wèi)-同心斷裂屬于左旋走滑斷層。牛首山-羅山斷裂帶位置在青銅峽和牛首山這2個地區(qū)的西南方向，向西南傾斜，傾角在62°～85°之間。在該區(qū)域巖石形變大，根據(jù)橫切關系顯示出2個明顯的構造變形期。牛首山-羅山斷裂帶是在漸新世之前形成的，為反向的左旋走滑斷層。

六盤山西部有許多斷層和褶皺，大多數(shù)是逆走滑構造。褶皺的軸面向西傾斜，說明它們受到了東西向的壓應力作用。六盤山地區(qū)構造應力整體為北東東向壓應力。研究區(qū)出露地層為白堊紀和新近地層，晚白堊世地層逆沖于新近地層之上，相對年齡確定主要構造活動出現(xiàn)在新近紀以后（圖2）。對比分析和應力分析發(fā)現(xiàn)，六盤山西部構造活動強烈，東部構造變形遠弱于西部，地層平緩?？梢酝茢?，青藏高原隆升的遠程效應空間僅限于六盤山西部。

2弧形構造形成機理

六盤山弧形構造帶的野外勘查發(fā)現(xiàn)早古生代地層已有大范圍的弧形斷層與褶皺發(fā)育。下古生界地層存在不整合，在二維平面上也能觀察到地層形態(tài)有一定的弧形展布特征。六盤山及其鄰區(qū)早古生代晚期的弧型構造，西部可寬達200 km左右，而東部可窄至40 km左右，弧型構造的西部是散布的，而東部向南則逐漸收緊聚集，弧形構造的規(guī)模比較大，且不屬于對稱的弧形構造。六盤山古生代弧型構造帶的構造應力主要從南東指向北西方向，即來自規(guī)模和強度都較大的北祁連加里東造山帶。

從平面分布上看，六盤山及鄰區(qū)弧形構造西翼較寬，南起香山、天井山地區(qū)，北至香根-達賴地區(qū)北部，南北方向約200 km。弧形構造東翼較窄，牛首山-達賴地區(qū)，小羅山-青龍山南北向，東西方向約40～50 km。

六盤山弧形構造帶古弧構造的形成是早古生代晚期中、北祁連地塊北東向擠壓對鄂爾多斯地塊西緣和阿拉善地塊南部的影響。鄂爾多斯地塊與阿拉善地塊的邊界分布非常特殊，在北東方向上形成了賀蘭拉張區(qū);北祁連弧后盆地在中、北祁連地塊的擠壓下逐漸關閉，形成了一系列以下古生界地層為主的逆沖推覆構造。中、北祁連地塊的逆沖楔繼續(xù)向北東方向擠壓并進入賀蘭拉張區(qū)，由于在逆沖推覆構造帶的走向方向上推進時遇到的地層厚度不同，導致逆沖推覆構造帶地體變形，在平面上形成了弧形的構造形態(tài)。

中生帶在古弧形構造的基礎上形成了現(xiàn)今六盤山弧形構造帶的主體部分。古太平洋板塊俯沖作用下，使得六盤山地區(qū)總體處于NW-NE向左旋擠壓的構造應力場環(huán)境，三疊紀的印支運動使該地區(qū)隆起，侏羅紀受到有限的拉張，早白堊世早期拉張斷陷，早白堊世晚期擠壓反轉，陸內(nèi)變形巨大，多個時期多次發(fā)生逆沖推覆現(xiàn)象。早白堊世時期弧形構造處于早期走滑拉張沉降階段，構造帶在強烈的擠壓作用下大幅抬升，處于短暫的應力松弛階段，在近東西方向的微弱的應力松弛產(chǎn)生的張力作用下，導致侏羅紀己經(jīng)形成的正斷層再次活動，具有走滑的性質(zhì)。東西向拉張作用繼續(xù)加大，六盤山弧形構造帶內(nèi)的正斷層開始變的活躍起來，六盤山盆地西側的邊界斷層由早期的左行走滑轉為右行走滑，六盤山盆地東側存在斷裂帶，在六盤山盆地南部整體遭受北東向壓應力后，六盤山弧形構造帶隨六盤山盆地發(fā)生了一定程度的逆時針旋轉。早白堊世燕山運動晚期弧形構造帶處于晚期湖盆萎縮消亡階段，六盤山盆地繼續(xù)受到北東向的擠壓應力的作用，六盤山盆地地層開始整體發(fā)生緩慢抬升。

3弧形構造數(shù)值模擬

3.1平面模型的建立

構建六盤山弧形構造帶古生代弧形構造形成初期的地質(zhì)模型。地質(zhì)模型忽略了弧形構造垂向隆升的變化及其對弧形構造形成的次要影響因素（圖4）。地質(zhì)模型使用二維有限差分程序FLAC軟件建立平面數(shù)值模型，將數(shù)值模擬研究的重點主要集中在確認早古生代晚期形成的六盤山弧形構造帶弧形的形態(tài)形成的力學機制。

在平面分布中明確構造帶弧形的彎折的力學條件是中、北祁連地塊向北東擠壓。利用FLAC2D構建地質(zhì)模型，考慮到數(shù)值模擬軟件基本計算單元網(wǎng)格劃分，按照地質(zhì)模型將區(qū)域地塊劃分為阿拉善地塊、鄂爾多斯地塊、中北祁連地塊、北祁連弧后盆地等地質(zhì)塊體，得到了六盤山弧形構造帶及其鄰區(qū)的平面模型圖。平面數(shù)值模型按1∶5 000的比例建立模型，最終模型分別代表東西向1 000 km，南北向760 km的研究范圍?；⌒螛嬙斓男纬勺鳛樵撃Ｐ偷哪M分析重點。數(shù)值模擬的主要目的是解釋其彎折的力學條件。研究區(qū)位于數(shù)值模型中心偏北東部，與模型4個邊界的距離都較大，主要是為了抵消邊界效應對研究區(qū)的影響。此外，賀蘭拉張區(qū)的砥柱作用沒有阿拉善地塊和鄂爾多斯地塊強烈，故將其并入北祁連弧后前陸盆地的一部分進行模擬（圖5）。

3.2數(shù)值模型初始條件和邊界

構建的數(shù)值模型北部為阿拉善地塊，東部為鄂爾多斯地塊，南部為中、北祁連地塊，中部為北祁連弧后前陸盆地，由于六盤山弧形構造帶形成于北祁連弧后盆地區(qū)域內(nèi)部，因此可以利用代表北祁連弧后前陸盆地的網(wǎng)格形變來反映六盤山弧形構造的形成過程。根據(jù)區(qū)域內(nèi)地塊的形狀特征，利用FLAC2D建立可計算的數(shù)值模型并生成計算單元網(wǎng)格。為了在構建的模型的北東向壓力作用下網(wǎng)格能產(chǎn)生充分形變且劃分盡量少的區(qū)塊以減少計算負荷，地質(zhì)塊體劃分左右邊界以南北向為主，整個模型劃分成子區(qū)塊，通過子區(qū)塊建立完整網(wǎng)格。模型中阿拉善地塊、鄂爾多斯地塊、北祁連弧后盆地、中北祁連地塊分別由不同的網(wǎng)格組成。

由于六盤山弧形構造帶北部受阿拉善地塊的阻擋，東部受鄂爾多斯地塊的阻擋，將模型邊界設置為南部邊界施加了均等應力;模型的東部和西部為滾動邊界，即在東西邊界處物質(zhì)不發(fā)生向東和向西運動。模型北部邊界為固定邊界，南部邊界的物質(zhì)運移是自由的。南部邊界的壓力被設置為2.0×107 Pa。由于各向同性彈性介質(zhì)的特點是受到應力作用時應力與應變呈現(xiàn)正相關的線性關系，應力均勻增大時應變也按一定的比率線性增大。本次數(shù)值模型將地質(zhì)塊體設定為各向同性彈性介質(zhì)來建立模型（表1），明確六盤山弧形構造帶形成的主要動力學過程。

3.3數(shù)值模型計算結果

六盤山弧形構造帶的數(shù)值模型進行模擬計算后模型處于平衡狀態(tài)，并得到了位移圖。導出x方向、y方向位移等值線圖，位移均為正值。對于y方向位移等值線圖，數(shù)值模型均為正值，即相對于初始階段的模型各網(wǎng)格節(jié)點向北位移，并且南部位移大于北部位移（圖6）。數(shù)值模型在x方向的位移計算結果表明六盤山區(qū)域范圍內(nèi)古生代地質(zhì)塊體位移的主要趨勢是向東，向東最大位移約為向北最大位移的25%（圖7）。將數(shù)值模擬計算結果中x方向位移與y方向位移合并后得到位移場圖，從位移向量圖指示的位移方向總體表現(xiàn)為向北東弧形突出（圖8）。

4結論

1）六盤山弧形構造帶經(jīng)歷了長期復雜的構造演化過程，其主體形成于印度-歐亞板塊碰撞所導致的青藏高原隆升對其東北緣產(chǎn)生的北東向的擠壓應力。六盤山及其鄰區(qū)早古生代晚期北東向擠壓變形的相對柔軟的地塊遇到堅硬的阿拉善地塊、鄂爾多斯地塊阻擋的共同作用下經(jīng)歷了早期的弧形構造。數(shù)值模擬研究古生代弧形構造形成初期的地球動力學條件對弧形構造的作用機制。

2）依據(jù)前人對六盤山及其鄰區(qū)早古生代晚期弧形構造研究成果構建了數(shù)值模型。六盤山弧形構造帶南部邊界不同時期受力大小與方向是存在變化的，其壓力大小雖然變化較大，但壓力方向為小幅度變化。數(shù)值模型將地質(zhì)塊體做為各向同性彈性介質(zhì)，并對不同地質(zhì)塊體賦予巖石力學參數(shù)，數(shù)值模型南部邊界施加的北東向應力為均等定值定向壓力。

3）數(shù)值模型模擬計算的位移結果與弧形的總趨勢相一致。數(shù)值模擬過程有效驗證了北東向擠壓作用是六盤山弧形構造帶古生代弧形構造形成主要因素。

參考文獻（References）：

[1]萬天豐.論六盤山——賀蘭山碰撞帶[J].大地構造與成礦學，2020，44（5）：845-851.WAN Tianfeng.Discussion on the Liupanshan-Helanshan collision zone[J].Geotectonica et Metallogenia，2020，44（5）：845-851.

[2]周鼎武，趙重遠，李銀德，等.鄂爾多斯西南緣地質(zhì)特征及其與秦嶺造山帶的關系[M].北京：地質(zhì)出版社，1994.

[3]張培震.六盤山地區(qū)的構造特征及地殼縮短量[M].北京：地震出版社，1990.

[4]劉鳳山，吳中海，張岳橋，等.青藏高原東緣新構造與活動構造研究新進展及展望[J].地質(zhì)通報，2014，33（4）：403-418.LIU Fengshan，WU Zhonghai，ZHANG Yueqiao，et al.New progress and prospects of neotectonics and active tectonics synthetically study on eastern edge Qinghai-Xizang plateau[J].Geological Bulletin of China，2014，33（4）：403-418.

[5]張培震，鄧起東，張國民，等.中國大陸的強震活動與活動地塊[J].中國科學D輯（地球科學），2003，33（增刊）：12-20.ZHANG Peizhen，DENG Qidong，ZHANG Guomin，et al.Active tectonic blocks and strong earth quakes in continent of China[J].Science in China（Series D），2003，33（S）：12-20.

[6]王東，張快樂，李偉.六盤山弧形構造帶幾何學、運動學及年代學研究[J].西部資源，2013（6）：84-86.WANG Dong，ZHANG Kuaile，LI Wei.Geometry，kinematics and chronology of Liupanshan arc-shaped structural belt[J].Western Resources，2013（6）：84-86.

[7]張韜，林暢松，杜玉良，等.六盤山弧形構造帶銅金屬礦床的成礦環(huán)境[J].地球學報，2002，23（6）：553-558.ZHANG Tiao，LIN Changsong，DU Yuliang，et al.Mineralization environments of the Cu polymetallic deposits in the Liupanshan arc structure belt[J].Acta Geoscientia Sinica，2002，23（6）：553-558.

[8]李強，江在森，武艷強，等.海原-六盤山斷裂帶現(xiàn)今構造變形特征[J].大地測量與地球動力學，2013，33（2）：18-22.LI Qiang，JIANG Zaisen，WU Yanqiang，et al.Present-day tectonic deformation characteristics of HaiYuan-Liupanshan fault zone[J].Journal of Geodesy and Geodynamics，2013，33（2）：18-22.

[9]王雙緒，蔣鋒云，郝明，等.青藏高原東緣現(xiàn)今三維地殼運動特征研究[J].地球物理學報，2013，56（10）：3334-3345.WANG Shuangxu，JIANG Fengyun，HAO Ming，et al.Investigation of features of present 3D crustal movement in eastern edge of Tibet plateau[J].Chinese Journal of Geophysics，2013，56（10）：3334-3345.

[10]徐化超，王輝，曹建玲.青藏高原東北緣主要斷裂滑動速率及其動力學意義[J].地震，2018，38（3）：13-23.XU Huachao，WANG Hui，CAO Jianling.Slip rates of the major faults in the northeastern Tibetan plateau and their geodynamic implications[J].Earthquake，2018，38（3）：13-23.

[11]李文輝，高銳，王海燕，等.六盤山斷裂帶及其鄰區(qū)地殼結構[J].地球物理學報，2017，60（6）：2265-2278.LI Wenhui，GAO Rui，WANG Haiyan，et al.Crustal structure beneath the Liupanshan fault zone and adjacent regions[J].Chinese Journal of Geophysics，2017，60（6）：2265-2278.

[12]張曉亮，師昭夢，蔣鋒云，等.海原-六盤山弧型斷裂及其附近最新構造變形演化分析[J].大地測量與地球動力學，2011，31（3）：20-24.ZHANG Xiaoliang，SHI Zhaomeng，JIANG Fengyun，et al.Research on late tectonic deformation evolvement of Huaiyuan-Liupanshan arc fault and its surrounding area[J].Journal of Geodesy and Geodynamics，2011，31（3）：20-24.

[13]LI Wei，DONG Yunpeng，GUO Anlin，et al.Chronology and tectonic significance of Cenozoic faults in the Liupanshan arcuate tectonic belt at the northeastern margin of the Qinghai-Tibet plateau[J].Journal of Asian Earth Sciences，2013，73：103-113.

[14]戴洪寶，許繼影.海原斷裂東段至六盤山斷裂西段GPS剖面地殼變形與應變積累分析[J].大地測量與地球動力學，2016，36（4）：343-345.DAI Hongbao，XU Jiying.Analysis of the GPS section deformation and strain accumulation of the eastern Haiyuan and Western Liupanshan faults[J].Journal of Geodesy and Geodynamics，2016，36（4）：343-345.

[15]陳虹，胡健民，渠洪杰，等.川滇南北向構造帶早中生代構造變形研究[J].中國科學D輯（地球科學），2011，41（9）：1281-1294.CHEN Hong，HU Jianmin，QU Hongjie，et al.Early Mesozoic structural deformation in the Chuandian N-S tectonic belt，China[J].Science China Earth Sciences，2011，41（9）：1281-1294.

[16]張岳橋，廖昌珍，施瑋，等.鄂爾多斯盆地周邊斷裂運動學分析與晚中生代構造應力體制轉換[J].地質(zhì)學報，2006，80（5）：639-647.ZHANG Yueqiao，LIAO Changzhen，SHI Wei，et al.Fault kinematic analysis and change in late Mesozoic tectonic stress regimes in the peripheral zones of the Ordos Basin，North China[J].Acta Geologica Sinica，2006，80（5）：639-647.

[17]向宏發(fā)，池田安隆，張晚霞，等.六盤山東麓的古地震研究[J].中國地震，1999，15（1）：74-81.XIANG Hongfa，YASUTAKA Ikeda，ZHANG Wanxia，et al.Study on paleoearthquakes of the eastern Liupanshan piedmont fault zone[J].Earthquake Research in China，1999，15（1）：74-81.

[18]李天斌.寧夏南部弧形推覆構造帶特征及演化[J].地質(zhì)力學學報，1999，5（3）：22-27.LI Tianbin.The feature and formation mechanism of arciform thrust-nappe structure zone of the southern NingXia[J].Journal of Geomechanics，1999，5（3）：22-27.

[19]李天斌.鄂爾多斯盆地西緣逆沖推覆構造特征及演化[D].北京：中國地質(zhì)大學（北京），2006.LI Tianbin.The characteristics and evolution of thrust nappe structure in the west margin of Ordos basin[D].Beijing：China University of Geoseienees（Beijing），2006.

[20]張進，馬宗晉，任文軍，等.寧夏中南部古生代弧型構造[J].大地構造與成礦學，2004（1）：29-37.ZHAGN Jin，MA Zongjin，REN Wenjun，et al.Paleozoic arcuate thrust belt in central and southern NingXia[J].Geotectonica et Metallogenia，2004，28（1）：29-37.

[21]ZHANG Jin，ZHANG Yinping，XIAO Wenxia，et al.Linking the Alxa terrane to the eastern gondwana during the early Paleozoic：constraints from detrital zircon U-Pb ages and cambrian sedimentary records[J].Gondwana Research，2015，28（3）：1168-1182.

[22]陳連旺，楊樹新，謝富仁，等.中國大陸構造應力應變場現(xiàn)今年變化特征的數(shù)值模擬[J].中國地震，2005，21（3）：341-349.CHEN Lianwang，YANG Shuxin，XIE Furen，et al.Numerical simulation on annual change patterns of present-day tectonic stress-strain field on the Chinese mainland[J].Earthqueake Research in China，2005，21（3）：341-349.

[23]劉鎏，魏東平.中國大陸及鄰區(qū)板內(nèi)應力場的數(shù)值模擬及動力機制探討[J].地震學報，2012，34（6）：727-740.LIU Liu，WEI Dongping.Numerical simulation of the intraplate stress field in China mainland and adjacent areas and its dynamic implication[J].Acta Seismologica Sinica，2012，34（6）：727-740.

[24]張東寧，許忠淮.中國大陸巖石層動力學數(shù)值模型的邊界條件[J].地震學報，1999，21（2）：21-27.ZHANG Dongning，XU Zhonghuai.Boundary conditions of the numerical model of lithosphere dynamics in mainland China[J].Acta Seismologica Sinica.1999，21（2）：21-27.

[25]馬德云，高振敏，楊世瑜，等.北衙金礦區(qū)構造應力場數(shù)值模擬[J].大地構造與成礦學，2003，27（2）：160-166.MA Deyun，GAO Zhenmin，YANG Shiyu，et al.Numerical simulation for tectonic stress field of the Beiya gold deposit[J].Geotectonica et Metallogenia，2003，27（2）：160-166.