馬爾可夫鏈在小南海壩基軟弱夾層分析中的應(yīng)用

聶瓊，項(xiàng)偉, ，杜水祥，顧晶

(1. 中國地質(zhì)大學(xué)(武漢) 工程學(xué)院，湖北 武漢，430074；2. 中國地質(zhì)大學(xué)(武漢) 教育部長江三峽庫區(qū)地質(zhì)災(zāi)害研究中心，湖北 武漢，430074)

重慶長江小南海水電站壩址位于重慶市江津區(qū)珞璜鎮(zhèn)下游約1.4 km，是國務(wù)院批準(zhǔn)的《長江流域綜合利用規(guī)劃簡要報(bào)告》推薦梯級開發(fā)方案的重要樞紐。建壩地層為侏羅系上統(tǒng)遂寧組(J3s)含礫砂巖、砂巖、粉砂巖、黏土質(zhì)粉砂巖與黏土巖互層等組成的一套河流相地層。該地層巖性、巖相變化較大，結(jié)構(gòu)上為軟硬相間的層狀巖體[1]，層間分布多條黏土巖類軟弱夾層，其性狀及分布特征決定了大壩及廠房的安全穩(wěn)定及經(jīng)濟(jì)和理性，是工程勘察過程中亟待研究和解決的工程地質(zhì)問題[2]。軟弱夾層的分布規(guī)律取決于整套巖層的沉積相[3]，由于沉積過程、沉積作用的不斷續(xù)性和事件沉積、背景沉積的交替進(jìn)行，地層巖相序列被認(rèn)為是在過程和作用上相互獨(dú)立的離散空間單元和離散時(shí)間單元序列，導(dǎo)致了巖相剖面的馬爾科夫性[4]。因此，許多學(xué)者采用馬爾可夫鏈數(shù)學(xué)模擬方法，分析巖相在垂向上的繼承關(guān)系、韻律結(jié)構(gòu)模式[5-7]，客觀地研究隨機(jī)轉(zhuǎn)化的相序列并建立相模型，恢復(fù)古環(huán)境[8-16]。前人的研究主要從巖性、粒度和沉積構(gòu)造方面將地層劃分為不同類型，進(jìn)行馬爾可夫鏈模型計(jì)算，但缺少對不同巖性狀態(tài)沉積厚度的考慮，且將此分析方法應(yīng)用于四川盆地紅層地區(qū)軟弱夾層分布規(guī)律的研究很少見。為此，本文作者應(yīng)用馬爾可夫鏈方法對小海南主要建壩地層J3S4-1剖面中軟弱夾層的沉積環(huán)境進(jìn)行概率動態(tài)模擬，充分考慮不同巖性層厚特征，結(jié)合實(shí)際地質(zhì)情況分析巖層的垂向序列分布規(guī)律及組合關(guān)系；綜合考慮軟弱夾層巖性、粒度成分、泥化程度、構(gòu)造特征及層序分析得到的軟弱夾層與上下巖層的組合關(guān)系，最終將壩址區(qū)軟弱夾層劃分為不同類型。

1 地質(zhì)概況

研究區(qū)域地貌屬長江流域第二階梯中部的四川盆地區(qū)[1]，處于龍泉山脈北東側(cè)與大婁山脈北西側(cè)之間，是最具代表性的近水平紅層分布區(qū)。小南海水電站位于觀音峽背斜與南溫泉背斜間的金鰲寺向斜內(nèi)，中壩址兩岸地貌為構(gòu)造剝蝕低山丘陵，壩址區(qū)相伴的地質(zhì)構(gòu)造主要有褶曲、斷層破碎帶、剪切帶及裂隙等。

出露的地層除第四系外，基巖主要為侏羅系上統(tǒng)遂寧組(J3s)和中統(tǒng)沙溪廟組(J2s)。侏羅世早期四川盆地為相對穩(wěn)定的陸相大型淡水湖相沉積，到侏羅世中期開始向河流相沉積轉(zhuǎn)變，侏羅中晚期湖相沉積范圍明顯縮小，河流相沉積范圍大為擴(kuò)展，表現(xiàn)為沙溪廟組總體上一套黃綠色灰紫色厚層砂巖與泥巖的組合，到侏羅晚期，為泛濫平原間夾湖相(時(shí)令湖)沉積，形成了侏羅世早期遂寧組厚300～500 m 鮮紅色細(xì)粒沉積。壩基巖體層狀構(gòu)造發(fā)育，尤其在J3s4-1段地層中多見粉細(xì)砂巖與黏土巖不等厚狀互層，由于物理力學(xué)性質(zhì)差異明顯，在構(gòu)造作用下，巖層間相互錯動容易形成巖性、結(jié)構(gòu)、厚度、延伸長度等性狀不同的軟弱夾層，是工程勘察中的重點(diǎn)研究對象。

2 J3s4-1 地層馬爾可夫鏈過程分析

小南海地層剖面在巖性及成因標(biāo)志上出現(xiàn)頻繁交互變化，閉合的半韻律旋回特征具有馬爾可夫鏈性質(zhì)。馬爾可夫鏈分析過程可以理解為將復(fù)雜的剖面地層按照不同的巖性狀態(tài)進(jìn)行分類，通過研究具有特定含義的巖性狀態(tài)的轉(zhuǎn)移屬性來把握剖面旋回特征、定量分析相變過程中的時(shí)空關(guān)系。具體方法包括極限概率分析、差值分析、置換分析和熵分析，各個(gè)環(huán)節(jié)的數(shù)學(xué)計(jì)算結(jié)果均賦以明確的地質(zhì)含義。

2.1 巖性狀態(tài)n 的選擇

研究地層劃分出的巖性狀態(tài)稱為系統(tǒng)的狀態(tài)n，狀態(tài)選擇是應(yīng)用馬爾可夫鏈分析地層的關(guān)鍵。通過對剖面巖性特征分析，本文選擇5 種巖性狀態(tài)，即含礫砂巖(Ssg)、砂巖-細(xì)砂巖(Ss)、粉砂巖(Fs)、黏土質(zhì)粉砂巖(Zf)、粉砂質(zhì)黏土巖-黏土巖(Cr)分別表示為狀態(tài)A，B，C，D 和E。單層層厚是巖相序列中一個(gè)重要的參數(shù)，主要反映了沉積規(guī)模、沉積速率及物源供給情況等，同時(shí)也是影響軟弱夾層形成的重要因素之一，因此，本文將剖面地層按單層層厚分為2 類：① 單層層厚≥0.5 m；② 單層層厚＜0.5 m。各個(gè)狀態(tài)地質(zhì)含義分述如下。

(1) A 雜色含礫砂巖：一般在砂巖底部夾有一層厚層灰綠色含礫砂巖條帶，礫石為紫紅色、灰黃色砂礫，大小混雜，粒徑為3～5 cm。通常河道流量已逐漸減小的情況下，沉積了含礫砂巖。由于大量粗碎屑快速分離出，使得流體密度大為降低，形成了較為正常的牽引流，這時(shí)含礫砂巖向上過渡為發(fā)育不好的邊灘沉積。

(2) B 灰白色、紫灰色砂巖-細(xì)砂巖：一般沉積厚度較大，夾少量黏土巖，槽狀交錯層理發(fā)育，是典型的邊灘沉積特點(diǎn)。也有較薄層的平行層理砂巖，出現(xiàn)在含礫砂巖之上，見零星分布的水流線理。

(3) C 紫紅色粉砂巖：C1 厚層粉砂巖多發(fā)育在洪泛平原沉積中，是洪水所攜帶的細(xì)、粉砂級物質(zhì)越過堤岸淹沒整個(gè)河間地形成的，分布范圍廣，厚度大，發(fā)育小型槽狀交錯層理及水平層理；C2 薄—中厚層粉砂巖在天然堤沉積及洪泛平原沉積中均有，薄層粉砂巖常與黏土巖形成砂、泥薄互層。

(4) D 紫紅色黏土質(zhì)粉砂巖：D 類巖石多發(fā)育在洪泛平原沉積中，D1 厚層黏土質(zhì)粉砂巖黏粒質(zhì)量分?jǐn)?shù)(粒徑＜5 μm)為24.39%～31.85%，性狀較C1 差；D2薄—中厚層黏土質(zhì)粉砂巖常夾在粉砂巖與粉砂質(zhì)黏土巖之間。

(5) E 紫紅色、磚紅色粉砂質(zhì)黏土巖、黏土巖：E1厚層黏土巖類一般在遠(yuǎn)離河床的積水低洼區(qū)沉積，稱洪泛盆地沉積，分布較穩(wěn)定，黏粒質(zhì)量分?jǐn)?shù)(粒徑＜5 μm)為36.31%～49.80%；E2 類薄—中厚層黏土巖類通常發(fā)育在堤岸及洪泛平原沉積中，黏粒質(zhì)量分?jǐn)?shù)(粒徑＜5 μm)為60.58%～63.91%，由于上、下所夾砂巖為良好的透水巖體，因此，其間的黏土巖類容易泥化，形成泥化夾層。

經(jīng)統(tǒng)計(jì)，J3s4-1地層總厚度45.76 m，共劃分出60層，其中A 類共1 層、B 類共1 層、C 類共17 層、D類共25 層、E 類共16 層，如表1 所示。

表1 J3s4-1 地層剖面巖性Table 1 Lithologies of J3s4-1

2.2 計(jì)算轉(zhuǎn)移頻數(shù)矩陣M

將各種巖性狀態(tài)按在巖層剖面中(按地層由老到新的順序)的彼此實(shí)際上覆接觸關(guān)系的次數(shù)統(tǒng)計(jì)出來，分別填入矩陣，稱轉(zhuǎn)移頻數(shù)矩陣M。計(jì)算出巖性狀態(tài)的轉(zhuǎn)移頻數(shù)矩陣M=[Mij]，Mij為狀態(tài)i 轉(zhuǎn)移到狀態(tài)j的頻數(shù)，且Mii=0。而

計(jì)算出J3S4-1剖面轉(zhuǎn)移頻數(shù)矩陣M 見表2。

表2 轉(zhuǎn)移頻數(shù)矩陣MTable 2 Transfer frequency matrix M

2.3 計(jì)算轉(zhuǎn)移概率矩陣P 和Q

轉(zhuǎn)移頻數(shù)矩陣每行總數(shù)除該行各元素，可得向上轉(zhuǎn)移概率矩陣P，表示某種巖性狀態(tài)被其他各種巖性狀態(tài)上覆的概率。而每列總數(shù)除該列各元素，可得向下轉(zhuǎn)移概率矩陣Q，表示各種巖性狀態(tài)被其中某一種巖性狀態(tài)上覆的概率。計(jì)算出的向上轉(zhuǎn)移概率矩陣和向下轉(zhuǎn)移概率矩陣分別為 P=[Pij]，Q=[Qij]，且Pii=Qii=0。

式中：P 為向上轉(zhuǎn)移概率矩陣；Q 為向下轉(zhuǎn)移概率矩陣。

計(jì)算出J3S4-1剖面轉(zhuǎn)移概率矩陣P 和Q 分別見表3和表4。

表3 向上轉(zhuǎn)移概率矩陣PTable 3 Upward transition probability matrix P

表4 向下轉(zhuǎn)移概率矩陣QTable 4 Downward transition probability matrix Q

2.4 χ2 檢驗(yàn)

經(jīng)計(jì)算求得自由度為49，χ2=73.3，查χ2分布表，

2.5 計(jì)算極限轉(zhuǎn)移概率

利用向上轉(zhuǎn)移概率矩陣自乘19 次后計(jì)算J3S4-1剖面的極限概率見表5。該剖面中占比例最大的是D2(薄—中厚層黏土質(zhì)粉砂巖)，其極限概率為0.221，其次為E2(薄—中厚層黏土巖類)，最小的是A(含礫砂巖)和B(砂巖類)，僅為0.016。由此可見剖面中主要巖性以薄—中厚層黏土質(zhì)粉砂巖及黏土巖類的細(xì)粒沉積為主，反映了大部分地層屬于水動力條件弱的洪泛平原沉積環(huán)境。通常取極限矩陣中的最大或最小概率的巖性作為旋回的起點(diǎn)，本文結(jié)合實(shí)際地層分析，將A 或B 作為旋回分析的起點(diǎn)。

表5 極限概率Table 5 Limit probability

從巖性狀態(tài)振蕩曲線(圖1)可以看出：J3S4-1剖面地層中以D2(薄—中厚層黏土質(zhì)粉砂巖)的振蕩幅度最大，且向著固定向量收斂速度較慢，說明D2 在地層序列中表現(xiàn)出明顯的韻律性，在旋回中具有較固定的位置；C1(厚層粉砂巖)、D1(厚層黏土質(zhì)粉砂巖)、E2(薄—中厚層黏土巖類)振幅較大，說明其巖性狀態(tài)旋回性清楚，在旋回中穩(wěn)定性較高；A(含礫砂巖)和B(砂巖類)的振蕩幅度最小且收斂快，在旋回中穩(wěn)定性低；而E1(厚層黏土巖類)為單向逼近，不具有旋回性。

圖1 J3S4-1 剖面巖性狀態(tài)振蕩曲線Fig.1 Lithological status vibration curve

2.6 計(jì)算插值矩陣D

利用插值矩陣較好地作出旋回模式，其方法是先由轉(zhuǎn)移頻數(shù)矩陣M=[Mij]構(gòu)造1 個(gè)E=[Eij]矩陣，其中

表示由狀態(tài)i 轉(zhuǎn)移到狀態(tài)j 的概率，然后用矩陣D 來表示轉(zhuǎn)移的傾向，即

計(jì)算出J3S4-1剖面插值矩陣D(表6)，根據(jù)表6 將差值矩陣中所有非零元素從大到小排成向量，見表7，將并取差值0.1 作為閥值，舍掉小于0.1 的實(shí)際觀察相變，作出比隨機(jī)更常出現(xiàn)的簡化的相變關(guān)系圖(圖2)。

從圖2 可見：(1) A(含礫砂巖)和B(砂巖類)作為顆粒最粗的巖性處在旋回的最底部，由于J3S4-1剖面地層以洪泛平原沉積為主，A 和B 僅出現(xiàn)1 次，但可以看出它們的上覆地層均為D2(薄—中厚層黏土質(zhì)粉砂巖)；(2) D2 與E2(薄—中厚層黏土巖類)常表現(xiàn)為粉砂質(zhì)黏土巖與黏土巖類薄互層；(3) C 類粉砂巖與D1(厚層黏土質(zhì)粉砂巖)常交替出現(xiàn)，代表了典型的洪泛平原細(xì)粒沉積；(4) E1(厚層黏土巖類)多形成于洪泛盆地沉積中，其上覆地層通常為D 類(黏土質(zhì)粉砂巖)。

表6 差值矩陣DTable 6 Difference matrix D

表7 差值矩陣D 的各正值元素Table 7 Positive value of difference matrix D

圖2 相變關(guān)系圖Fig.2 Relation of facies change

2.7 計(jì)算置換矩陣L，R 和C

置換分析是以轉(zhuǎn)移矩陣P 為基礎(chǔ)，用矩陣?yán)碚搧硌芯康貙悠拭嫘蛄械哪撤N共性的方法，其著眼點(diǎn)是反映狀態(tài)間共性的P 矩陣中，行與行或列與列的元素的相似性概念，其計(jì)算公式如下。

(1) 左置換矩陣。在地層觀測序列中，不同狀態(tài)后面常出現(xiàn)同一種后繼狀態(tài)時(shí)，彼此可以置換，由于都位于同一后繼狀態(tài)的左邊，故稱為左置換。

(2) 右置換矩陣。同樣，在不同狀態(tài)的前面常出現(xiàn)同一種前行狀態(tài)時(shí)，彼此間可以置換，因都位于同一前行狀態(tài)的右邊而稱之為右置換。

(3) 互置換矩陣。在地層觀測序列中，不同狀態(tài)前后常出現(xiàn)相同的狀態(tài)組合時(shí)，彼此間可以置換，因同時(shí)考慮了左、右2 種可置換性，故稱為互置換。若Cij很高，說明2 種狀態(tài)i 和j 可以互置換。在地質(zhì)意義上，表明i 和j 兩狀態(tài)之間發(fā)生相變的可能性大，因此，可以用一種狀態(tài)代表，即可簡化巖性分類，以突出沉積序列的旋回性。若Cij很低，則表明該狀態(tài)很穩(wěn)定，在地層序列研究中可作為標(biāo)志層。

Lij和Rij均介于0～1之間，越接近于1，表示可置換性越大，且Lii=Rii=1。計(jì)算出J3S4-1剖面置換矩陣L，R 和C，分別見表8～10。

表8 左置換矩陣LTable 8 Left replacement matrix L

表9 右置換矩陣RTable 9 Right replacement matrix R

表10 互置換矩陣CTable 10 Mutual replacement matrix C

左置換矩陣L(表8)表明：狀態(tài)C2(薄—中厚層粉砂巖)和E1(厚層黏土巖類)在概率為0.95 的條件下有共同的頂層巖性D1(厚層黏土質(zhì)粉砂巖)。右置換矩陣R(表9)表明：狀態(tài)B(砂巖類)和E2(薄層黏土巖類)在概率為0.85 的條件下有共同的底層巖性D2(薄—中厚層黏土質(zhì)粉砂巖)。互置換矩陣C(表10)表明：狀態(tài)B(砂巖類)和E2(薄—中厚層黏土巖類)概率值最大，為0.72，即B 和E2 常常緊密共生，有共同的頂、底層巖性D2(薄—中厚層黏土質(zhì)粉砂巖)。巖性類別相近的狀態(tài)互相代替，反映了相變的存在；而巖性類別差別較大的狀態(tài)相互代替，反映出巖相橫向變化的不穩(wěn)定性。C1(厚層粉砂巖)和D1(厚層黏土粉砂巖)的概率值最小，表明C1 和D1 相似性較小，即它們有不同的頂?shù)谞顟B(tài)，該層位穩(wěn)定性好。

2.8 熵分析

在隨即實(shí)驗(yàn)范圍內(nèi)，熵是一種不確定性程度的度量。為了確定沉積類型，Hattori[18]于1976 年提出2種熵類型：一種為沉積后熵(Epost)，另一種為沉積前熵(Epre)。其表達(dá)式分別為

式中：Eipost和 Eipre分別表示狀態(tài)i 出現(xiàn)之后和出現(xiàn)之前的巖性轉(zhuǎn)移的變化規(guī)律，若 Eipost＞ Eipre，說明狀態(tài)i 后繼狀態(tài)隨機(jī)性較大；若 Eipost＜ Eipre，則說明狀態(tài)i 的后繼狀態(tài)比前行狀態(tài)有較大的確定性。J3S4-1剖面熵分析結(jié)果如表11 所示。

表11 熵分析結(jié)果Table 11 Result of entropy analysis

圖3 熵相關(guān)圖Fig.3 Entropy correlation

3 軟弱夾層分布特征

重慶長江小南海水電站軟弱夾層從成因來考慮屬于原生構(gòu)造型。它是巖體在構(gòu)造應(yīng)力作用下，沿軟、硬巖層接觸帶或軟巖內(nèi)部發(fā)生層間剪切錯動，剪切錯動帶受到多期構(gòu)造變動而發(fā)生剪切破壞，形成大量細(xì)顆粒物質(zhì)和裂隙，經(jīng)地下水的滲透和物理化學(xué)作用而使原生夾層軟化、泥化形成的一種復(fù)雜的軟弱夾層。由于野外J3S4-1段地層軟弱夾層樣品多取自壩址區(qū)左岸淺表層，受構(gòu)造應(yīng)力作用影響不大，但受地表風(fēng)化及水巖相互作用強(qiáng)烈，物理化學(xué)性質(zhì)較壩基軟弱夾層變化較大。因此，本文在考慮軟弱夾層巖性、粒度成分、泥化程度、構(gòu)造特征等因素的基礎(chǔ)上，結(jié)合地層層序分析得到的軟弱夾層與上下巖層的組合關(guān)系將壩址區(qū)軟弱夾層劃分為三大類。

(1) Ⅰ類破碎夾層(巖塊巖屑型)。Ⅰ類夾層多為黏土質(zhì)粉砂巖在構(gòu)造應(yīng)力作用下，沿軟、硬巖層接觸帶發(fā)生層間剪切錯動而形成的破碎帶，多同時(shí)具有節(jié)理帶、劈理帶和泥化光面等分帶現(xiàn)象。物質(zhì)組成以粗顆粒為主，細(xì)顆粒成分較少，黏粒極少。

壩址區(qū)巖塊巖屑型軟弱夾層中粉粒質(zhì)量分?jǐn)?shù)多大于黏粒質(zhì)量分?jǐn)?shù)，多發(fā)育于D2(薄—中厚層黏土質(zhì)粉砂巖)和E1(厚層黏土巖類)中。從巖性狀態(tài)振動曲線(圖1)可以看出：E1 由于處在洪泛盆地沉積中，不具有韻律性，因此，E1 中沉積的Ⅰ類夾層分布不連續(xù)。由相變關(guān)系圖(圖2)可看出，D2 為E1 的上覆地層時(shí)，由于處在與粉砂巖的軟硬過渡層之間，也易形成Ⅰ類夾層，因此，此類夾層只可能對壩基局部穩(wěn)定性帶來影響。

從形態(tài)特征上來看(圖4(a))，Ⅰ類夾層野外露頭宏觀上均可見其裂隙發(fā)育，往往形成一條基本連續(xù)或斷續(xù)的主剪裂面，主剪裂面一般彎曲起伏，起伏差幾毫米到幾厘米不等，面上一般斷續(xù)附泥膜。主剪裂面一側(cè)或兩側(cè)裂隙發(fā)育，據(jù)裂隙發(fā)育程度可劃分為碎裂帶和節(jié)理帶。近主剪裂面?zhèn)攘严睹芗?，形成網(wǎng)紋格架，切割形成菱形或透鏡體狀碎塊或碎屑，并大致定向呈帶狀與主剪裂面平行展布，裂隙間距幾厘米不等，碎屑、碎塊間多填泥，形成手感軟化的碎裂帶。透鏡體狀或菱形狀碎塊往往可見沿主剪裂面方向偏轉(zhuǎn)的跡象，表明后期剪切位移較大。碎裂帶外為節(jié)理稀疏的節(jié)理帶，裂隙間距幾十厘米不等。

(2) Ⅱ類破碎夾泥層(巖屑夾泥型)。Ⅱ類夾層實(shí)為Ⅰ類夾層進(jìn)一步錯動形成。這一類型軟弱夾層一般在軟巖內(nèi)部發(fā)育，頂?shù)撞客鶠閹r質(zhì)相對較硬的C(粉砂巖)或B(砂巖類)。破碎夾層在層間剪切作用下進(jìn)一步破碎產(chǎn)生粉粒，同時(shí)產(chǎn)生許多裂隙，地下水活動加強(qiáng)，在地下水的作用下部分產(chǎn)生泥化。組成物質(zhì)以細(xì)碎屑為主，細(xì)顆粒中黏粒質(zhì)量分?jǐn)?shù)和粉粒質(zhì)量分?jǐn)?shù)多小于10%。如圖4(b)所示，Ⅱ類夾層界面實(shí)為密集裂隙錯斷巖層形成，界面由裂隙面和裂隙間巖塊錯斷面組成，起伏較大，起伏差幾毫米至幾厘米。頂?shù)捉缑嫱鶖嗬m(xù)附泥膜，在起伏的拐點(diǎn)段往往充填泥夾碎屑，泥膜厚度幾毫米不等。碎屑排列具定向性，泥膜面可見擦痕。

(3) Ⅲ類泥化夾層。Ⅲ類泥化夾層構(gòu)造特征與巖屑夾泥型相似，一般發(fā)育在E(黏土巖類)中。極限概率分析可知E1(厚層黏土巖類)不具有旋回性，而E2(薄—中厚層黏土巖類)旋回性清楚，且在旋回中穩(wěn)定性較高；差值矩陣和熵分析均表明E2 與D2 常表現(xiàn)為薄互層的半韻律，由于E2 上、下所夾砂巖為良好的透水巖體，因此，其間的黏土巖類物質(zhì)容易泥化，形成Ⅲ類泥化夾層。韻律層厚度最大為1.48 m，最小為0.05 m，且韻律層越厚此類夾層出現(xiàn)頻率越大。

圖4 軟弱夾層分類Fig.4 Classification of weak intercalation

根據(jù)Ⅲ類泥化夾層黏粒質(zhì)量分?jǐn)?shù)將其分為Ⅲ1泥夾巖屑型和Ⅲ2全泥型。Ⅲ1類泥化夾層如圖4(c)所示。碎塊巖屑間填充泥質(zhì)物較多，呈泥包碎塊狀，黏粒質(zhì)量分?jǐn)?shù)多在10%～30%之間，上下層面多見薄層泥化帶；Ⅲ2全泥型軟弱夾層是薄層黏土巖、粉砂質(zhì)黏土巖或破碎夾泥層后期經(jīng)地下水及風(fēng)化作用而發(fā)生軟化和泥化，形成一定厚度、斷續(xù)或連續(xù)的泥層，如圖4(d)所示，由粗碎屑、細(xì)碎屑、粉粒、黏粒組成，黏粒質(zhì)量分?jǐn)?shù)大于30%，黏土巖底部有1～5 mm 厚的泥膜。

4 結(jié)論

(1) 小南海壩基J3S4-1段地層沉積物較細(xì)，一般為泥、砂沉積，主要形成于河流相中的洪泛平原亞相，沉積厚度大，分布廣泛，為軟弱夾層的形成提供了物質(zhì)基礎(chǔ)。

(2) J3S4-1段地層的韻律結(jié)構(gòu)模式主要分2 種：一種為粉砂巖-厚層黏土質(zhì)粉砂巖互層半韻律沉積，韻律層厚度大，代表了典型的洪泛平原細(xì)粒沉積；另一種為黏土質(zhì)粉砂巖-黏土巖薄互層半韻律沉積，韻律層厚度小，為泥化夾層的形成提供了良好的條件。

(3) 軟弱夾層形成于特定的巖相沉積環(huán)境下，在考慮軟弱夾層巖性、粒度成分、泥化程度、構(gòu)造特征等因素的基礎(chǔ)上，結(jié)合地層層序分析得到的軟弱夾層與上下巖層的組合關(guān)系把J3S4-1段地層中的軟弱夾層劃分為Ⅰ破碎夾層(巖塊巖屑型)、Ⅱ破碎夾泥層(巖屑夾泥型)和Ⅲ泥化夾層(Ⅲ1泥夾巖屑型和Ⅲ2全泥型)三大類，其中Ⅲ1類泥化夾層發(fā)育最為廣泛。

(4) Ⅰ類夾層多形成于洪泛盆地沉積的厚層黏土巖中或與厚層粉砂巖交界處，僅對工程的局部穩(wěn)定性產(chǎn)生影響；Ⅱ類夾層實(shí)為Ⅰ類夾層進(jìn)一步錯動形成，一般在軟巖內(nèi)部發(fā)育，頂?shù)撞客鶠閹r質(zhì)相對較硬的粉砂巖或砂巖；Ⅲ類軟弱夾層主要形成于黏土質(zhì)粉砂巖與黏土巖薄互層的半韻律沉積中，分布較穩(wěn)定，對工程穩(wěn)定性十分不利。

[1] 王子忠. 四川盆地紅層巖體主要水利水電工程地質(zhì)問題系統(tǒng)研究[D]. 成都: 成都理工大學(xué)地質(zhì)工程系, 2011: 6-10.WANG Zizhong. Systematic researches on red bed rock mass engineering geological problems of water resources and hydropower projects in Sichuan Basin, China[D]. Chengdu:Chengdu University of Technology. Department of Geological Engineering, 2011: 6-10.

[2] 柴波, 殷坤龍. 三峽庫區(qū)巴東新城區(qū)庫岸三疊系巴東組層間軟弱帶[J]. 工程地質(zhì)學(xué)報(bào), 2009, 17(6): 809-816.CHAI Bo, YIN Kunlong. Interlayer weakness zones in Badong formation of middle Triassic forming bank slopes of Three Gorges Reservoir in new Badong Country[J]. Journal of Engineering Geology, 2009, 17(6): 809-816.

[3] 宋自成, 鄭德順, 李明龍, 等. 濟(jì)源盆地中侏羅統(tǒng)馬凹組沉積環(huán)境演化及其控制因素分析[J]. 河南理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 31(5): 543-550.SONG Zicheng, ZHENG Deshun, LI Minglong, et al. Analysis on the evolution of sedimentary environments and its controlling factors of the Middle Jurassic Maao Formation in Jiyuan Basin[J]. Journal of Henan Polytechnic University (Natural Science), 2012, 31(5): 543-550.

[4] 龔一鳴. 巖相剖面的定量歷史沉積學(xué)分析[J]. 地球科學(xué)——武漢地質(zhì)學(xué)院學(xué)報(bào), 1987, 12(6): 613-620.GONG Yiming. The quantitative analysis of historical sedimentology for lithofacies successions[J]. Earth Sciences:Journal of Wuhan College of Geology, 1987, 12(6): 613-620.

[5] YE Ming, Khaleel R. A Markov chain model for characterizing medium heterogeneity and sediment layering structure[J]. Water Resources Research, 2008, 44(9): 1-15.

[6] MEI Mingxiang, MA Yongsheng, GUO Qingyin, et al. Basic lithofacies-succession model for the Wumishan cyclothems:Their Markov chain analysis and regularly vertical stacking patterns in the third-order sequences[J]. Acta Geologica Sinica(English Edition), 2001, 75(4): 421-431.

[7] Hein F J, Robb G A, Wolberg A C, et al. Facies descriptions and associations in ancient reworked transgressive shelf sandstones Cambrian and Cretaceous examples[J]. Sedimentology, 1991,38(3): 405-431.

[8] Summer E J, Talling P J, Amy L A, et al. Facies architecture of individual basin-plain turbidites: Comparison with existing models and implications for flow processes[J]. Sedimentology,2012, 59(6): 1850-1887.

[9] Charvin K, Gallagher K, Hampson G L, et al. A Bayesian approach to inverse modelling of stratigraphy. Part 1: Method[J].Basin Research, 2009, 21(1): 5-25.

[10] Verwer K, Porta G D, Oscar M T, et al. Controls and predictability of carbonate facies architecture in a Lower Jurassic three-dimensional barrier-shoal complex (Djebel Bou Dahar,High Atlas, Morocco)[J]. Sedimentology, 2009, 56(6):1801-1831.

[11] Chutko K J, Lamoureux S F. Biolaminated sedimentation in a high arctic freshwater lake[J]. Sedimentology, 2009, 56(6):1642-1654.

[12] Riegl B M, Purkis S J. Markov models for linking environments and facies in space and time (Recent Arabian Gulf, Miocene Paratethys)[J]. Int Assoc Sedimentol Spec Publ, 2009, 41:337-360.

[13] 張成鳳, 呂洪波, 夏邦棟, 等. 南盤江盆地中三疊統(tǒng)復(fù)理石韻律的馬爾科夫鏈特征及其地質(zhì)意義[J]. 地質(zhì)評論, 2011, 57(5):632-640.ZHANG Chengfeng, Lü Hongbo, XIA Bangdong, et al. Markov chain simulation on the rhythmic sequences of middle Triassic flysch in Nanpanjing basin, SW China[J]. Geological Review,2011, 57(5): 632-640.

[14] 王秋良, 李長安. 馬爾柯夫鏈在江漢平原第四紀(jì)沉積環(huán)境分析中的應(yīng)用[J]. 地質(zhì)科技情報(bào), 2008, 27(1): 38-42.WANG Qiuliang, LI Changan. Markov chain in quaternary sedimentary environment of Jianghan plain[J]. Geological Review, 2008, 27(1): 38-42.

[15] 郭飛, 葛成, 韓宇. 嵌入式馬爾科夫鏈的地質(zhì)屬性建模與應(yīng)用[J]. 地理與地理信息科學(xué), 2012, 28(1): 47-50.GUO Fei, GE Cheng, HAN Yu. Markov chain-based geological properties modeling and its application[J]. Geography and Geo-Information Science, 2012, 28(1): 47-50.

[16] 張成鳳, 呂洪波. 馬爾科夫鏈模擬在濟(jì)陽坳陷中的應(yīng)用[D].青島: 中國石油大學(xué)地質(zhì)系, 2010: 5-9.ZANG Chengfeng, Lü Hongbo. Application of Markov Chain simulation in the Jiyang depression[D]. Qingdao: China University of Petroleum. Department of Geology. 2010: 5-9.

[17] Anderson T W, Goodman L A. Statistical inference about Markov chinas[J]. Ann Math Statist, 1957, 28(1): 89-110.

[18] Hattori I. Entropy in Markov chains and discrimination of cyclic patterns in lithologic succession[J]. Math Geol, 1976, 8(4):477-497.