四川省劍門關(guān)區(qū)晚侏羅—早白堊世沉積體系的數(shù)學模擬

楊 帥，謝小平，李 冰，李佳麗，王 芳

(曲阜師范大學 地理與旅游學院， 山東 日照 276826)

?

四川省劍門關(guān)區(qū)晚侏羅—早白堊世沉積體系的數(shù)學模擬

楊 帥，謝小平，李 冰，李佳麗，王 芳

(曲阜師范大學 地理與旅游學院， 山東 日照 276826)

四川省劍門關(guān)區(qū)沉積體系受龍門山控制，區(qū)內(nèi)分布來自龍門山的碎屑沉積物，以礫巖和砂巖為主。由于晚侏羅—早白堊世的濕熱氣候條件，導致劍門關(guān)區(qū)內(nèi)形成獨特的紅色巖層，經(jīng)過內(nèi)外營力剝蝕，形成獨特的丹霞地貌。劍門關(guān)區(qū)巖層成層性良好，分界線明顯，自下而上巖層組依次為上侏羅統(tǒng)遂寧組、蓮花口組、下白堊統(tǒng)劍門關(guān)組。依據(jù)物質(zhì)恒定原理對區(qū)內(nèi)沉積物進行數(shù)學模擬，建立劍門關(guān)區(qū)物質(zhì)恒定模型，確定沉積物主要堆積于115～100 Ma年之間，以此確定地質(zhì)事件的發(fā)生時間和古環(huán)境變化；運用馬爾科夫原理對區(qū)內(nèi)垂向上巖性、構(gòu)造進行數(shù)學模擬，建立馬爾科夫模型，得到巖性、構(gòu)造垂向上的變化規(guī)律。

數(shù)學模擬；沉積體系；地質(zhì)事件；四川省劍門關(guān)區(qū)

四川省劍門關(guān)區(qū)位于龍門山構(gòu)造帶北東段前緣，四川盆地的西北緣，處在龍門山傳遞帶[1-2]前端低凹處，屬于橫向逆沖斷層傳遞帶的下盤(圖1)。自20世紀30年代地質(zhì)學家趙亞曾、黃汲清越嶺入川，對劍門關(guān)區(qū)進行地質(zhì)資源調(diào)查，到侯德封、王現(xiàn)珩[3-5]命名的“劍門關(guān)礫巖”，劍門關(guān)地質(zhì)方面的研究開始被諸多學者所關(guān)注。1982年劍門關(guān)地質(zhì)公園核心景區(qū)劍門蜀道風景區(qū)被列為國家重點風景名勝區(qū)，自此劍門關(guān)區(qū)地質(zhì)遺跡資源逐一被發(fā)掘，地質(zhì)景觀得到開發(fā)。21世紀初，劍門關(guān)區(qū)旅游業(yè)發(fā)展迅速，地質(zhì)旅游人數(shù)越來越多，地質(zhì)遺跡遭到相當?shù)钠茐?。研究機構(gòu)和學者紛紛對劍門關(guān)區(qū)地質(zhì)遺跡資源保護開發(fā)方面提出學術(shù)性的建議，期間學術(shù)性論文出現(xiàn)數(shù)10篇，但大多是對區(qū)內(nèi)地質(zhì)遺跡資源開發(fā)保護方面，純地質(zhì)的甚少，僅有徐世球、王紅梅等[6-9]分別在沉積、紅層分子化石和磁化率方面進行研究。

圖1 龍門山北段傳遞帶簡圖

數(shù)學地質(zhì)屬地質(zhì)學科新的分支，是地質(zhì)學向著定量化方向發(fā)展的產(chǎn)物，主要運用各種模型對地質(zhì)演化進行模擬，能夠較準確地預測演化過程和變化規(guī)律，對古環(huán)境分析和古地理恢復具有重要的意義。數(shù)學地質(zhì)的首次應用者是美國沉積學家W.C.Krumbein，被公認為“數(shù)學地質(zhì)之父”。隨后在國際地質(zhì)大會上成立國際數(shù)學地質(zhì)協(xié)會，自此數(shù)學地質(zhì)正式被國際上認可。我國數(shù)學地質(zhì)先驅(qū)者為劉承祚、趙鵬大等[10-17]，分別在數(shù)學模型在地質(zhì)中應用、礦產(chǎn)統(tǒng)計學等方面卓有建樹。我國數(shù)學地質(zhì)工作始于20世紀60年代初，70年代得到了大發(fā)展，八九十年代以來，進入了一個更高水平和更深入發(fā)展的新階段[18]。

1 劍門關(guān)區(qū)沉積背景

劍門關(guān)區(qū)位于龍門山前陸盆地的西北緣，沉積體系受龍門山構(gòu)造帶控制(圖2)。龍門山逆沖推覆構(gòu)造帶屬于上揚子古陸塊，揚子陸塊東南緣從寒武—早奧陶世末被動大陸邊緣轉(zhuǎn)化為晚奧陶世—志留紀的前陸盆地[19]。中生代早起的晚三疊世龍門山受到來自北西方向的擠壓作用，褶皺變形，此時龍門山南東方向的四川盆地還處于海洋時期。中三疊世龍門山進一步擠壓褶皺，形成山脈，使四川盆地西緣受到相應的擠壓并伴隨海拔升高，海水退出四川盆地，結(jié)束了碳酸鹽為主的臺地沉積階段，進入了陸相沉積盆地的發(fā)展階段。劍門關(guān)區(qū)的沉積背景也在此時正式形成。李勇等學者根據(jù)龍門山前陸盆地充填序列中不整合面的層位和性質(zhì)，確定龍門山?jīng)_斷帶自晚三疊世諾利克期以來，至少有9次逆沖推覆事件，晚侏羅世蓮花口沉積期的推覆事件是規(guī)模較大的事件之一[20-21]。晚侏羅世龍門山構(gòu)造帶處于活動時期，構(gòu)造運動使龍門山抬升，逆沖推覆運動使其在南東方向產(chǎn)生巨大擠壓力，造成劍門關(guān)區(qū)地質(zhì)事件頻繁發(fā)生，沉積環(huán)境頻繁變化，從而產(chǎn)生一系列沉積相及其微相。每一次地質(zhì)事件的發(fā)生都會產(chǎn)生地質(zhì)要素的變化，形成一套有規(guī)律的韻律層。劍門關(guān)區(qū)受龍門山逆沖推覆影響，巖層巖性、構(gòu)造發(fā)生變化，巖性在垂向上表現(xiàn)為粗—細—粗循環(huán)，構(gòu)造表現(xiàn)為板狀—槽狀—波狀交替進行。數(shù)學地質(zhì)中認為這些變化隨著地質(zhì)事件呈有規(guī)律的變化，運用數(shù)學方法可以模擬這些變化。

圖2 劍門關(guān)區(qū)域構(gòu)造簡圖

2 劍門關(guān)區(qū)沉積體系

根據(jù)地層之間接觸關(guān)系以及巖性特征，劍門關(guān)區(qū)侏羅統(tǒng)巖層不整合超覆在三疊統(tǒng)巖層之上，上侏羅統(tǒng)—早白堊統(tǒng)之間存在上超和削截兩種不整合接觸(圖3)[22]，每種不整合接觸都表明發(fā)生過一次大的沉積事件。削截的巖層被上覆的往前推移的巖層覆蓋，往上巖石粒徑會變大，意味在削截的層面往上發(fā)生了正旋回。上超的巖層被往后推移的巖層覆蓋，往上巖石粒徑會變小，意味在上超的層面往上發(fā)生了負旋回。劍門關(guān)區(qū)垂向上沉積體系由各種沉積相組成，包括洪積扇相、河流相、湖泊相和三角洲相。

圖3 上超和削截示意圖

2.1 洪積扇相

沖積扇相發(fā)育在上侏羅統(tǒng)蓮花口組的下部，是強降水沖開出山口帶走巖屑和礫石，以水、泥、石混合物形式流動，在流動過程中形成的沉積模式，屬于近源沉積。劍門關(guān)區(qū)洪積扇相巖性組合為礫巖、砂礫巖和粉砂巖，且分選性、磨圓度極差，存在板狀、波狀交錯層理，總厚度在800 m左右。根據(jù)沉積特征將洪積扇相分為扇根、扇中、扇端亞相，主體為扇中亞相。

2.2 河流相

河流相發(fā)育在上侏羅統(tǒng)蓮花口組中、上部，是在水流穩(wěn)定流動狀態(tài)下形成的沉積模式，屬于中源沉積。劍門關(guān)區(qū)河流相巖性組合為礫巖、砂礫巖、粗砂巖和粉砂巖，分選性和磨圓度隨離龍門山由近及遠分別表現(xiàn)為差—良好—較好，存在槽狀、板狀、波狀交錯、平行層理，厚度在30 m左右。根據(jù)沉積特征將河流相分為曲流河、辮狀河亞相，劍門關(guān)區(qū)以曲流河亞相為主。

2.3 湖泊相

湖泊相發(fā)育在上侏羅統(tǒng)的遂寧組上部，是河流攜帶泥沙進入湖泊而形成的沉積模式，屬于遠源淡水湖泊沉積模式。劍門關(guān)區(qū)湖泊相巖性組合為粉砂巖、泥巖，分選性、磨圓度極好，存在波狀、板狀交錯、平行層理，并有大量生物潛穴和擾動發(fā)育，厚度在80 m左右。根據(jù)沉積特征景湖泊相分為湖濱和淺湖亞相，劍門關(guān)區(qū)以淺湖亞相為主。

2.4 三角洲相

三角洲相發(fā)育在上侏羅統(tǒng)遂寧組下部，處在河流和湖泊交匯處，受到河流和湖流共同影響而形成的沉積模式。劍門關(guān)區(qū)三角洲相巖性組合為砂巖、粉砂巖和泥巖，分選性、磨圓度良好，存在波狀交錯和平行層理，厚度在30 m左右。根據(jù)沉積特征分析，劍門關(guān)區(qū)三角洲相分為前三角洲、三角洲前緣和三角洲平原亞相，主體為三角洲平原亞相。由于受到湖浪、潮汐作用造成湖進和湖退兩種模式，從而形成局部的正旋回和負旋回的交替存在。

3 沉積相的數(shù)學模擬

地質(zhì)中的數(shù)學模擬就是運用數(shù)學方法預測成礦巖層、演化地質(zhì)作用過程、模擬地層垂向規(guī)律以及沉積過程。本文運用兩種數(shù)學方法模擬劍門關(guān)區(qū)晚侏羅世—早白堊世沉積演化過程及沉積特點，進而推測這一時期的古環(huán)境變化。物質(zhì)恒定模型和馬爾科夫模型是在沉積學領(lǐng)域應用較廣的數(shù)學方法，前者模擬沉積物的時空展布，后者模擬沉積物垂向變化。

3.1 物質(zhì)恒定模型

3.1.1 基本原理

物質(zhì)恒定模型是定量描述沉積過程的基本數(shù)學方法，將沉積過程看作剝蝕與沉積的對立統(tǒng)一過程，剝蝕產(chǎn)生沉積物來源，沉積是剝蝕物的歸宿(圖4)[22]。物質(zhì)恒定模型認為一定時間內(nèi)發(fā)生沉積過程的范圍內(nèi)沉積物總量保持不變，在這段時間中沉積物和剝蝕物相等。如圖4所示，把沉積過程分為3個階段：階段1：總沉積物為剝蝕；階段2：t0時沉積物發(fā)生剝蝕，剝蝕物在t1時沉積，且t0時剝蝕物等于t1時沉積物；階段3：t0、t1時沉積物發(fā)生剝蝕，剝蝕物在t2時沉積，t0、t1時剝蝕物等于t2時沉積物。

3.1.2 模擬過程

根據(jù)物質(zhì)恒定模型，單位時間內(nèi)剝蝕減少的沉積巖數(shù)量為定值，沉積巖總數(shù)量不變，因此

(1)

式中：M為沉積巖質(zhì)量；t表示時間；α為沉積系數(shù)/剝蝕系數(shù)，負號表明巖石剝蝕過程。

對式(1)進行積分可得

(2)

式中：Mt為巖石剝蝕后在時間t剩余的質(zhì)量；M0為t=0時沉積巖數(shù)量。

圖4 物質(zhì)恒定模型示意圖

在物質(zhì)恒定模型中沉積物的半剝蝕期為一定的時間間隔，這一時間間隔中沉積物質(zhì)的一半由于侵蝕作用而消失，半剝蝕期與沉積系數(shù)之間存在以下關(guān)系：

(3)

式中t0.5為半剝蝕期。

將t0.5看作未知數(shù)求解得

根據(jù)物質(zhì)恒定模型原理：一定時間沉積下來的物質(zhì)數(shù)量等于剝蝕的物質(zhì)數(shù)量。因此，在每一個相同時間的增量沉積下來的物質(zhì)數(shù)量等于沉積巖總數(shù)量減去剝蝕后剩余的巖石數(shù)量，即

(4)

式中：Md為每個時間間隔沉積物數(shù)量； Δt表示時間間隔。

為研究地質(zhì)事件序列可將研究時間劃分為n個時間段，每一個時間增量為Δt，并將其編號為i(其中i=1,2,3，…，n)。當i=1時，沉積巖年代最老；當i=n時，沉積巖年代最年青；在第i個時間增量中沉積的巖石年代為(n-i)Δt，并且任一時間增量期間沉積下來的物質(zhì)數(shù)量為Md?？梢杂嬎阍诘趇時間增量期間沉積下來的物質(zhì)經(jīng)過剝蝕以后到目前剩余下來的物質(zhì)數(shù)量為Mi，即

(5)

將式(4)代入式(5)，得

(6)

式(6)為物質(zhì)恒定模型的基本方程。劍門關(guān)區(qū)沉積模式遵循物質(zhì)恒定模型原理，研究晚侏羅—早白堊世的沉積特征，地質(zhì)時間序列為年150 Ma前至100 Ma之間的50 Ma(Δt=5 Ma，n=10)。劍門關(guān)區(qū)上侏羅—早白堊統(tǒng)地層垂向上在2 800 m左右，本文用其中間值對應的時間為半剝蝕期，根據(jù)化石測年可得t0.5=20 Ma，則α=0.034 66可得劍門關(guān)區(qū)數(shù)據(jù)，見表1。

表1 劍門關(guān)物質(zhì)恒定模型數(shù)據(jù)

根據(jù)表1進行劍門關(guān)區(qū)物質(zhì)恒定模型的線性模擬(圖5)。由圖5可以看出：沉積物主要沉積于7～9年代之間，即110～100 Ma之間，對于研究劍門關(guān)區(qū)晚侏羅—早白堊世的沉積特征、地質(zhì)事件及構(gòu)造過程具有重要意義。

圖5 劍門關(guān)區(qū)物質(zhì)恒定模型線性模擬

3.2 馬爾科夫模型

3.2.1 基本原理

馬爾科夫模型是概率統(tǒng)計學中預測發(fā)展趨勢常用的模型，里丁(Reading)和沃克(Walker)較早運用統(tǒng)計學的方法來確定不同沉積相的組合關(guān)系。沉積相研究中的馬爾科夫模型的關(guān)鍵在于狀態(tài)的確定和概率矩陣的計算，然后通過門坎值作出沉積相序的馬爾科夫鏈，以此為根據(jù)來反映地層沉積序列的沉積規(guī)律。

3.2.2 模擬過程

巖性、構(gòu)造是研究沉積相中主要的相標志，根據(jù)劍門關(guān)區(qū)晚侏羅—早白堊世沉積巖的巖性、構(gòu)造分析，可將巖性狀態(tài)確定為疊瓦狀礫巖、槽狀交錯砂礫巖、板狀及槽狀交錯粗砂巖、板狀及槽狀交錯中砂巖、板狀及槽狀交錯細砂巖、水平交錯粉砂巖、水平及板狀交錯泥巖7種狀態(tài)，本文分別用A、B、C、D、E、F、G代表。根據(jù)劍門關(guān)剖面晚侏羅世—早白堊世巖心巖性分析，確定巖層自下而上7種狀態(tài)分布(表2)。

表2 劍門關(guān)區(qū)狀態(tài)分布

據(jù)劍門丹霞地貌晚侏羅世—早白堊世遂寧組、蓮花口組、劍門關(guān)組的巖性、構(gòu)造狀態(tài)分布以及其他相標志得到觀察到的相變矩陣(表3)，通過馬爾科夫鏈的相變概率公式求出觀察到的轉(zhuǎn)移概率矩陣(表4)，公式如下：

式中：N為相變總數(shù)，nij是i相轉(zhuǎn)變?yōu)閖相次數(shù)；M為i相轉(zhuǎn)變?yōu)閖，j+1，j+2，… 相的總次數(shù)。

根據(jù)所有相變均是隨機的假設(shè)求出隨機層序的相變概率(表5)，公式如下：

式中：N為相變總數(shù)；ni、nj分別為i、j相出現(xiàn)的次數(shù)。

表3 觀察到的相變矩陣

表4 觀察到的相變概率矩陣

表5 隨機層序的相變矩陣

求出觀察的和隨機的相變概率之差(表6)，此差值范圍在-1～1之間。差值為正值表示觀察到的相變比隨機相變常見，差值為負值表示隨機相變比觀察到的相變常見。

選取觀察到的相變概率和隨機相變概率差值為正值的某數(shù)作為門坎值，去掉小于門坎值的實際觀察相變，作出差值大于門坎值的相變關(guān)系圖。本文選取0.10作為門坎值，得到相變關(guān)系圖，如圖6所示。

表6 觀察到的和隨機的相變之差矩陣

圖6 劍門關(guān)區(qū)馬爾科夫鏈

根據(jù)劍門關(guān)區(qū)馬爾科夫鏈表示的相變關(guān)系可以看出：巖性及構(gòu)造的垂向變化呈現(xiàn)一定規(guī)律，最主要的規(guī)律是疊瓦狀礫巖—板狀及槽狀交錯粗砂巖—板狀及槽狀交錯中砂巖—板狀及槽狀交錯細砂巖—槽狀交錯砂礫巖—疊瓦狀礫巖，其次是板狀及槽狀交錯中砂巖—水平交錯粉砂巖—水平及板狀交錯泥巖—板狀及槽狀交錯中砂巖，疊瓦狀礫巖—槽狀交錯砂礫巖—疊瓦狀礫巖。

4 結(jié)論

1) 劍門關(guān)區(qū)的沉積物的沉積時間主要發(fā)生在115～100 Ma之間，約占總沉積物的50.1%；

2) 劍門關(guān)區(qū)垂向上巖性變化呈粗—中—細—粗的規(guī)律，粒徑劇變說明有大的地質(zhì)事件發(fā)生，由此可推論劍門關(guān)區(qū)在晚侏羅—早白堊世發(fā)生的地質(zhì)事件；

3) 沉積特征表現(xiàn)為自下而上的正旋回沉積，局部有負旋回沉積；

4) 構(gòu)造的變化意味劍門關(guān)區(qū)構(gòu)造活動運動頻繁，導致巖石層理呈現(xiàn)各種變化，以槽狀及板狀交錯層理為主要標志。

[1] 陳發(fā)景,賈慶素,張洪年.傳遞帶及其在砂體發(fā)育中的作用[J].石油與天然氣地質(zhì),2004(2):144-148.

[2] 陳昭年,陳發(fā)景,王琦.正斷層軟聯(lián)接及其傳遞帶類型[J].現(xiàn)代地質(zhì),2005(4):495-499.

[3] 侯德封.四川北部之三疊紀地層[J].地質(zhì)論評,1939(2):97-100.

[4] 侯德封,楊敬之.四川盆地中的幾種地形與其形成史[J].地質(zhì)論評,1939(5):315-321,395-397.

[5] 侯德封.四川北部古生代地層的兩個剖面[J].地質(zhì)論評,1939(6):423-434,515-516.

[6] 徐世球,佘繼完.四川省劍門關(guān)一帶晚侏羅世沉積體系與古環(huán)境演化[J].地球科學,2001(3):235-240.

[7] 楊家,徐世球.川黔湘交境寒武紀二級層序的劃分及海平面變化[J].地球科學,1997(5):23-27.

[8] 徐世球,李富強,肖詩宇.川黔湘交境早寒武世“變馬沖組”沉積層序與古環(huán)境[J].地球科學,1997(5):77-82.

[9] 王紅梅,劉育燕,王志遠.四川劍門關(guān)侏羅——白堊系紅層分子化石的古環(huán)境和古氣候意義[J].地球科學,2001(3):229-234.

[10]劉承祚.數(shù)學地質(zhì)在1997—1999年期間的主要進展[J].地質(zhì)論評,2000(S1):1-5.

[11]劉承祚,張彥波,薛恩.數(shù)學方法在礦產(chǎn)地質(zhì)工作中的應用[J].地質(zhì)科學,1979(2):167-176.

[12]劉承祚.數(shù)學地質(zhì)的回顧與展望——一門地質(zhì)學的、新興的邊緣學科在我國的產(chǎn)生與發(fā)展[J].地質(zhì)論評,1982(3):286-290.

[13]劉承祚.沉積模擬的數(shù)學模型及其在我國中部一個含油氣沉積盆地中的應用[J].河北地質(zhì)學院學報,1988(3):9-19.

[14]趙鵬大.大數(shù)據(jù)時代數(shù)字找礦與定量評價[J].地質(zhì)通報,2015(7):1255-1259.

[15]趙鵬大.成礦定量預測與深部找礦[J].地學前緣,2007(5):1-10.

[16]趙鵬大.數(shù)字地質(zhì)與礦產(chǎn)資源評價[J].地質(zhì)學刊,2012(3):225-228.

[17]趙鵬大,夏慶霖.中國學者在數(shù)學地質(zhì)學科發(fā)展中的成就與貢獻[J].地球科學,2009,34(2):225-231.

[18]劉承祚.數(shù)學地質(zhì)的主要進展和發(fā)展趨勢[J].地質(zhì)論評,1996(4):364-368.

[19]潘桂棠,肖慶輝,陸松年,等.中國大地構(gòu)造單元劃分[J].中國地質(zhì),2009,36(1):1-4.

[20]李勇,周榮軍,DENSMORE A L,等.青藏高原東緣龍門山晚新生代走滑-逆沖作用的地貌標志[J].第四紀研究,2006,26(1):40-51.

(責任編輯 陳 艷)

Mathematical Modeling of Sedimentary System of the Late Jurassic-Early Cretaceous in Jianmenguan Area of Sichuan

YANG Shuai, XIE Xiaoping, LI Bing, LI Jiali, WANG Fang

(School of Geography and Tourism, Qufu Normal University, Rizhao 276826, China)

Sichuan customs Jianmen sedimentary system was controlled by Longmen mountain, where distributed classic sediments from Longmen mountain, given priority to the conglomerate and sandstone. Because of the humid climate conditions in the late Jurassic to early Cretaceous, there are the unique red stratums in the Jianmenguan scenic spot, which formed a unique Danxia landform through the denudation of the internal and external force. The rock stratum group of the area of Jianmenguan, with excellent layered property and obvious boundary, are the upper Jurassic Suining Formation, Lianhuakou Formation and the lower Cretaceous Jianmenguan Formation successively from top to bottom. This article, based on the principle of material constant proceeding mathematical modeling for sediment zone, establishes Jianmen customs material constant model, and determines the main sediment accumulation in the 115 Ma to 100 Ma years, in order to determine the occurrence time and paleoenvironmental changes of geological events; using the Markoff principle in the area to simulate vertical lithology and structure, we establish Markoff model to get change law of lithology and structure vertically.

mathematical modeling; sedimentary system; geological events; Jianmenguan area of Sichuan

2016-12-25

國家自然科學基金資助項目(41072164)

楊帥(1991—)，男，碩士，主要從事地質(zhì)學研究，E-mail:653383422@qq.com。

楊帥，謝小平，李冰，等.四川省劍門關(guān)區(qū)晚侏羅—早白堊世沉積體系的數(shù)學模擬[J].重慶理工大學學報(自然科學)，2017(7):133-139.

format：YANG Shuai, XIE Xiaoping, LI Bing, et al.Mathematical Modeling of Sedimentary System of the Late Jurassic-Early Cretaceous in Jianmenguan Area of Sichuan[J].Journal of Chongqing University of Technology(Natural Science)，2017(7):133-139.

10.3969/j.issn.1674-8425(z).2017.07.021

P534.52

A

1674-8425(2017)07-0133-07