塔里木盆地塘西北反沖斷層及其形成機(jī)制

陳書平 云金表 劉志娜 劉士林 李偉 季弘瑩

摘 要：前陸褶皺帶反向沖斷作用的形成機(jī)理一直沒有得到很好的認(rèn)識。塔里木盆地塘西北構(gòu)造帶發(fā)育很好的反沖斷層，給這種構(gòu)造形成機(jī)理的探討提供了很好的實(shí)例。地震資料顯示，塘西北構(gòu)造帶發(fā)育多排逆沖斷層及相關(guān)構(gòu)造，逆沖斷層傾向北西，向造山帶方向即南東向逆沖，斷層形成時間為奧陶紀(jì)晚期，受阿爾金造山帶或西昆侖造山帶的影響，是一組典型的反沖斷層。地質(zhì)分析及離散元顆粒流模型計算證明，反沖斷層形成的可能原因包括兩個方面，一是寒武系底部膏鹽層滑脫面的存在；二是與阿爾金斷裂帶相關(guān)，塘南古隆起對塘北構(gòu)造區(qū)的不均勻擠壓。

關(guān)鍵詞：反沖斷層 滑脫面 差異擠壓 機(jī)制 離散元模擬

中圖分類號：P542.3 文獻(xiàn)標(biāo)識碼：A 文章編號：1674-098X（2018）01（c）-0127-05

前沖斷層是沖斷方向與擠壓力方向相同的逆沖斷層，其傾向造山帶方向，向前陸方向沖斷。反沖斷層則相反，其傾向前陸向造山帶方向逆沖。在前陸盆地的盆山耦合帶，多數(shù)情況下發(fā)育前沖斷層，但反沖斷層也時有發(fā)生，如塔里木盆地玉東-瑪東斷裂帶發(fā)育向腹陸阿爾金山方向逆沖的沖斷層[1]、揚(yáng)子北緣弧形構(gòu)造發(fā)育的隔檔式褶皺也伴生反沖斷層[2]、大巴山前陸弧形構(gòu)造帶同樣也伴生反沖斷層[3]。國外同樣不乏這樣的實(shí)例[4，5]。但是，這些反沖斷層形成的原因仍然沒有得到很好的解釋。本文基于地震資料和已有的研究成果，并應(yīng)用數(shù)字砂箱離散元計算方法[6-8]，分析研究塘西北構(gòu)造帶變形時間、變形巖石組成以及變形邊界條件（見圖1），以期對該區(qū)反沖斷層的形成機(jī)理提出新認(rèn)識，也給其他類似構(gòu)造發(fā)育地區(qū)的構(gòu)造成因解釋提供參考。這里的塘西北構(gòu)造帶包括通常意義的玉東斷裂帶、瑪東斷裂帶、塘北斷裂帶、中3-塔中3井?dāng)嗔褞А?/p>

1 反沖構(gòu)造特點(diǎn)

塘西北構(gòu)造帶由多條逆沖斷層組成，構(gòu)造帶走向由北東向轉(zhuǎn)變?yōu)楸睎|東向，整體表現(xiàn)為向北西方向突出的弧形，被走向與構(gòu)造帶近垂直的平移斷層切割。剖面上，逆沖斷層傾向北西或北北西，向南或南南東方向逆沖（見圖2）。由于其變形力源來自南東方向的阿爾金褶皺帶或西昆侖褶皺帶，因此這些逆沖斷層是典型的反沖斷層。

不整合顯示，構(gòu)造主形成期為奧陶紀(jì)晚期的加里東運(yùn)動（見圖2）。變形巖石組成為寒武系和奧陶系，由砂巖、泥巖和碳酸鹽巖組成。中下寒武統(tǒng)（T6-T9上部），發(fā)育多層膏鹽，累積厚度在400～700m，且在背斜下方有加厚現(xiàn)象。這些膏鹽層起滑脫面的作用，分隔上下層不同變形樣式?；搶拥拇嬖谝彩欠礇_斷層形成的必要條件[4，5，9]。

塘西北構(gòu)造帶變形動力條件與阿爾金構(gòu)造帶的形成息息相關(guān)（見圖3），是與阿爾金斷裂帶形成相關(guān)的動力條件遠(yuǎn)程效應(yīng)的結(jié)果，而阿爾金斷裂帶的形成則與阿爾金地體、昆侖地體（西昆侖、東昆侖和柴達(dá)木盆地）與塔里木地體的碰撞作用相關(guān)。中奧陶世末期，在祁連阿爾金北昆侖洋盆俯沖消減作用下，阿爾金地體與塔里木地體拼貼，塔里木盆地由被動陸緣的伸展體制轉(zhuǎn)變?yōu)閿D壓的克拉通沉積盆地[10，11]，盆地南部發(fā)生臺盆分異[12-14]，形成隆升幅度較大的和田古隆起和孤立型塘南古隆起（臺地）。阿爾金山鋯石U-Pb年代學(xué)研究，也證明了450～500Ma的熱事件[15]。塔里木盆地東南緣區(qū)晚奧陶世，西昆侖東段發(fā)生的板塊碰撞作用，得到了碰撞型花崗巖類發(fā)育的證明，如中粒含斑黑云母二長花崗巖、中粒似斑狀黑云母二長花崗巖、中細(xì)粒黑云母二長花崗巖等[16]。

2 反沖斷層形成機(jī)理討論

上地殼巖石變形符合庫倫摩擦理論，斷層的形成呈共軛狀，兩條斷層的交線對應(yīng)中間主應(yīng)力軸，銳角平分線對應(yīng)最大主應(yīng)力軸。另外，理論上認(rèn)為，重力勢能增加量一定的情況下，緩傾斜逆沖斷層比陡傾斜逆沖斷層更容易發(fā)育，因?yàn)榍罢吣軌蛭崭嗟乃娇s短量[4，5]。因此，巖性均質(zhì)、厚度均勻、水平擠壓下，形成兩組傾角相同、同等發(fā)育的斷裂（見圖4a）。在厚度不均勻，且考慮不同摩擦系數(shù)底面情況下，雖然區(qū)域?yàn)樗綌D壓力，但局部最大主壓應(yīng)力軸會發(fā)生傾斜，這時所伴生的兩組共軛斷層的活動性，會有所不同。當(dāng)最大主壓應(yīng)力向前陸方向抬起時，伴生前沖斷層（見圖4b）。當(dāng)最大主壓應(yīng)力向腹陸造山帶方向抬起時，伴生后沖斷層（見圖4c）。

在實(shí)際情況中，多種因素可能影響著最大主應(yīng)力軸向腹陸方向傾斜，如變形體形狀、底面摩擦、施力邊界形狀等。Mackay（1995）認(rèn)為，斜坡面上的沉積楔形體，當(dāng)?shù)酌嫔夏Σ亮苄r，可以使最大主應(yīng)力軸向腹陸方向傾斜[5]。Byrne等（1993）認(rèn)為，底面摩擦力很小、施力邊界為楔形時，同樣可以使最大主應(yīng)力向施力端（腹陸）傾斜，形成反沖斷層[4]。

對于塘西北構(gòu)造帶，作為滑脫面的寒武系中、下統(tǒng)，發(fā)育厚層膏鹽層。膏鹽層具有低的抗變形能力，在地質(zhì)時間尺度可以看作流體。從現(xiàn)今變形層厚度看（見圖2），構(gòu)造發(fā)育期不存在斜坡上的變形楔形體。相反，應(yīng)該是厚度均勻的近水平地質(zhì)體。另一方面，中奧陶世末期形成的塘南古隆起[1]，對于塘西北構(gòu)造帶來說，是楔形施力邊界，類似于Byrne等（1993）[4]模擬試驗(yàn)中的邊界條件。

3 反沖構(gòu)造形成條件離散元模擬

為了驗(yàn)證上述對反沖構(gòu)造形成條件的推斷，設(shè)計了離散元顆粒流模型，進(jìn)行了模擬計算。

3.1 原理簡介

在離散元顆粒流模型中，利用綁定顆粒來代替天然巖石，以研究宏觀變形的微觀機(jī)理和預(yù)測這些宏觀變形的行為[8]。該方法由Cundall創(chuàng)立[6，7]，作為離散元的一種，主要應(yīng)用于模擬巖石類材料基本特性、顆粒物質(zhì)動力響應(yīng)、巖石類介質(zhì)內(nèi)斷層、節(jié)理及其他地質(zhì)構(gòu)造的發(fā)展等基礎(chǔ)性問題。其原理是，巖石的力學(xué)行為表現(xiàn)為微觀裂紋的形成、生長和相互作用。這種行為可以用綁定顆粒模型來代替。在牛頓運(yùn)動定律（力-加速度關(guān)系）、力-位移（應(yīng)力-應(yīng)變）關(guān)系約束下進(jìn)行計算，可以得到速度場和應(yīng)變場，從而研究巖石的宏觀變形。該方法稱為顆粒流離散元模擬，由程序PFC2D和PFC3D實(shí)現(xiàn)二維和三維計算，已廣泛應(yīng)用于斜坡變形和巖石變形研究[18-22]。

3.2 模擬結(jié)果分析

離散元顆粒流模擬過程中，對比試驗(yàn)了斜邊推進(jìn)和旋轉(zhuǎn)斜邊推進(jìn)、底面存在塑性層和不存在塑性層等，發(fā)現(xiàn)只有底面存在塑性層和旋轉(zhuǎn)斜邊施力時，才形成向施力端（腹陸）方向沖斷的沖斷層，或陡翼在施力端（腹陸）一側(cè)的褶皺。下面主要介紹一下最后這個模型的結(jié)果。

模擬模型為二維模型，56m長，4.65m高（見圖6a）。顆粒半徑0.03～0.0498m，粒徑服從高斯分布?？紫抖?6%，脆性層彈性模量80GPa，泊松比0.3，摩擦系數(shù)0.57；塑性層彈性模量1GPa，泊松比0.3，摩擦系數(shù)0.34；組內(nèi)顆粒間粘結(jié)力3×104Pa。1、2層為塑性層，3～8為脆性層。1層與底部邊界摩擦力為0，推板整體向前（左）推移速度為2mm/s，同時擋板圍繞底部A點(diǎn)做順時針旋轉(zhuǎn)，角速度為6.7×10-5rad/s，相當(dāng)于上部比A點(diǎn)向前推進(jìn)的速度慢，模擬施力端施力不均勻。

早期變形都是集中在推板附近，形成前沖斷層（見圖6b）。進(jìn)一步擠壓，變形向前傳遞，形成褶皺。在圖6c中，后翼（靠近推板）傾角42°，前翼傾角22°，已顯示出反向褶皺（斷層）的特點(diǎn)。圖6d中，后翼（靠近推板）傾角58°，前翼傾角33°，反向褶皺（斷層）的特點(diǎn)已經(jīng)很清楚了。

離散元顆粒流模擬計算證明，兩個因素在構(gòu)造形成中起重要作用，一是高塑性層的存在；二是深淺層不均勻擠壓。

4 結(jié)語

塘南構(gòu)造帶反沖斷層及向腹陸倒向的褶皺，與兩個因素有關(guān)，一是寒武系底部發(fā)育的厚層膏鹽層；二是阿爾金山造山帶擠壓古塘南隆起，塘南古隆起給塘西北構(gòu)造帶施加由深到淺不均勻擠壓，下部擠壓強(qiáng)（快），上部擠壓弱（慢）。這種不均勻擠壓相當(dāng)于在原來水平擠壓的基礎(chǔ)上，附加了層間滑動引起的剪切力，引起最大主壓應(yīng)力方向發(fā)生偏轉(zhuǎn)，使向腹陸沖斷的沖斷層（反向斷層）或陡翼靠近腹陸的褶皺更容易發(fā)生。

參考文獻(xiàn)

[1] 郭穎，湯良杰，余騰孝.等.塔里木盆地塘古巴斯坳陷瑪東構(gòu)造帶斷裂特征及成因探討[J].大地構(gòu)造與成礦學(xué)，2016，40（4）：643-653.

[2] 王平，劉少峰，鄭洪波，等.揚(yáng)子北緣晚造山階段弧形構(gòu)造特征與盆地演化[J].古地理學(xué)報，2013，15（6）：819-838.

[3] 王瑞瑞，張岳橋，解國愛，等.大巴山前陸弧形構(gòu)造的成因：來自砂箱實(shí)驗(yàn)的認(rèn)識[J].地質(zhì)學(xué)報，2011，85（9）：1409-1419.

[4] Byrne，D.E.，W.Wang，D.M.Davis.Mechanical role of backstops in the growth of forearcs[J]. Tectonics，1993，12（1）：123-144.

[5] Mackay，M.E.Structural variation and landward vergence at the toe of the Oregon accretionary prism[J].Tectonics，1995，14（5）：1309-1320.

[6] Cundall PA.A computer model for simulating progressive large scale movements in blocky rock systems[J].In： Proceedings of the Symposium of International Society of Rock Mechanics，1971，11（ii-b）：8-11.

[7] Cundall PA，Strack ODL.A discrete numerical model for granular assemblies[J].Geotechnique，1979，29（1）：47-65.

[8] Potyondy，D.O.，Cundall，P.A.Abonded-particle model for rock[J].International Journal of Rock Mechanics & Mining Sciences，2004（41）：1329-1364.

[9] Seely，D.R.The significance of landward vergence and oblique structural trends on trench inner slopes[J].In Island Arcs，Deep Sea Trenches，and Back-Arc Basins，Maurice Ewing Ser.，1977（1）：187-198.

[10]Torsvik TH，Cocks LRM.Gondwana from top to base in space and time[J].Gondwana Research，2013，4（3-4）：999-1030.

[11]張光亞，劉偉，張磊，等.塔里木克拉通寒武紀(jì)奧陶紀(jì)原型盆地、巖性古地理與油氣[J].地學(xué)前緣，2015， 22（3）：269-276.

[12]羅金海，雷剛林，劉良，等.阿爾金構(gòu)造帶對塔東南油氣地質(zhì)條件的制約[J].大地構(gòu)造與成礦學(xué)，2009，33（1）：76-85.

[13]許志琴，李思田，張建新，等.塔里木地塊與古亞洲/特提斯構(gòu)造體系的對接[J].巖石學(xué)報，2011，27（1）：1-22.

[14]Lin C.S.，Yang H.J.，Liu J.Y.，et al. Distribution and erosion of the Paleozoic tectonicunconformities in the Tarim Basin， Northwest China：Significance for the evolution of paleo-uplifts and tectonic geography during deformation[J].Journal of Asian Earth Sciences， 2012（46）：1-19.

[15]張建新，李懷坤，孟繁聰，等.塔里木盆地東南緣（阿爾金山）“變質(zhì)基底”記錄的多期構(gòu)造熱事件：鋯石U-Pb年代學(xué)的制約[J].巖石學(xué)報，2011，27（1）：23-46.

[16]崔建堂，邊小衛(wèi)，王根寶.西昆侖地質(zhì)組成與演化[J].陜西地質(zhì)，2006，24（1）：1-11.

[17]王步清，潘正中，郭群英，等.塔東南構(gòu)造演化和上新世以來的構(gòu)造格架[J].新疆石油地質(zhì)，2007，28（6）：704-706.

[18]Strayer，L.M.，Suppe，J.Out-of-plane motion of a thrust sheet during along-strike propagation of a thrust ramp：a distinct-element approach[J].Journal of Structural Geology， 2002（24）：637-650.

[19]Hardy，S.，McClay，K.Munoz，J.A.Deformation and fault activity in space and time in high-resolution numerical models of doubly vergent thrust wedges[J].Marine and Petroleum Geology，2009（26）：232-248.

[20]Zhang，J.，Morgan，J.K.，Gray，G.G.，et al. Comparative FEM and DEM modeling of basement-involved thrust structures， with application to Sheep Mountain，Greybull area，Wyoming[J].Tectonophysics，2013（608）：408-417.

[21]Liu，Z.N.，Koyi，H.A.，Swantesson，J.O.H.，et al. Kinematics and 3-D internal deformation of granular slopes：analogue models and natural landslides[J].Journal of Structural geology， 2013（53）：27-42.

[22]Liu，Z.N.，Koyi，H.A.The impact of a weak horizon on kinematics and internal deformation of a failure mass using discrete element method[J].Tectonophysics，2013（586）：95-111.