桂西南地區三疊系沉積型錳礦床成因研究

尹青，伊海生，周凱

1.成都理工大學沉積地質研究院，成都 610059 2.中國石化中原油田分公司勘探開發研究院，河南濮陽 457001 3.成都理工大學油氣藏地質及開發工程國家重點實驗室，成都 610059

桂西南地區三疊系沉積型錳礦床成因研究

尹青1,2，伊海生1,3，周凱2

1.成都理工大學沉積地質研究院，成都 610059 2.中國石化中原油田分公司勘探開發研究院，河南濮陽 457001 3.成都理工大學油氣藏地質及開發工程國家重點實驗室，成都 610059

桂西南地區錳礦床資源富集并具有較好的找礦前景，然而目前有關該區三疊系主要沉積地層的地層格架劃分尚較薄弱且鮮有探討。利用桂西南地區東平錳礦區四口鉆孔的自然伽馬測井曲線進行去趨勢、高通濾波處理以及累計殘差計算等方法，對桂西南地區含錳巖系進行了旋回地層學研究。研究結果表明：含錳巖系下三疊統北泗組共可劃分為四個海進—海退旋回層序，碳酸錳礦層發育的旋回層序位置主要出現于各旋回層序內海退沉積序列中，并位于旋回層序內海退—海侵界面附近，錳礦沉積主要發生在相對海平面由緩慢下降至逐漸上升的時期的轉換時期。綜合研究區錳礦形成的海平面變化、古地理條件、錳質來源、氧化還原條件等特征的基礎上，筆者提出了研究區錳礦成因模式，即具有地理隔絕形式的、并由同生斷裂控制的“局限海”臺盆相碳酸錳礦沉積模式，碳酸錳礦的形成主要經歷以下3個階段：1)海平面下降—氧化錳礦物形成階段；2)海平面上升—氧化錳礦物轉化碳酸錳礦物階段；3)最高海泛面—少量碳酸錳礦物發育階段。該研究成果可為桂西南地區今后錳礦勘探工作中的地層分析、建立區域內高頻層序單元格架及錳礦成因等研究提供一種新的方法和手段。

桂西南地區；三疊系；含錳巖系；旋回地層；海平面變化

桂西南地區是我國重要的優質富錳礦產區之一，該區發育的錳礦床(礦點)多、儲量大，已探明的錳礦儲量近2億噸，錳礦資源量居全國之首。且錳礦賦存層位較多，根據產出特征錳礦類型總體可分為地下深部沉積型碳酸錳礦床[1-3]和地表淺部次生氧化錳礦床[4-6]兩大類。三疊系是桂西南地區重要的含錳礦地層之一，含礦層位為下三疊統北泗組，是近碳酸鹽臺地的淺海臺盆相碳酸鹽巖建造。其下部夾4～13層含錳硅質巖或含錳灰巖，淺部形成錳帽型氧化錳礦。隨著近年來桂西南地區地表次生氧化錳礦床的開采接近枯竭，以及鋼鐵、錳系合金等錳礦深加工產品的需求量持續加大，錳礦已被國家列為緊缺礦種之一，亟需后續大量的錳資源來支持。因此，桂西南地區下三疊統北泗組錳礦勘查工作也由地表氧化錳礦床逐漸轉入深部沉積碳酸錳礦床，其地質研究工作也日益受到人們的重視。

然而目前有關桂西南地區三疊系地層格架劃分尚較薄弱，沉積旋回和錳礦床的沉積成因鮮有探討。研究區下三疊統含錳巖系巖性較為單調，總體為灰質、泥質、硅質、碳質等混合沉積，錳礦層雖然數量較多但標志層不明顯。如何進行井間地層以及礦層對比一直是制約該地區錳礦勘探的難題。對于研究區沉積型碳酸錳礦床盆地，旋回地層研究可作為其重要的地層分析方法之一。由于旋回地層學用于識別和對比沉積旋回周期性變化具有較高的精確性和分辨率，其層序劃分方案具有較統一和科學的定量方法，可大大減少了人為因素在高分辨率地層劃分和對比研究中的影響[7]，近年來在前寒武紀至第四紀地層高精度劃分與對比的研究中得到廣泛的應用[8-11]。

2013年，由中國冶金地質總局承擔，成都理工大學與四川省核工業地質調查院共同負責完成，首次將油田、煤田及鈾礦等廣泛采用的測井技術應用到桂西南地區錳礦床試驗研究。本文應用研究區四口鉆孔的自然伽馬測井曲線，采取頻譜分析、濾波處理及累計殘差等方法，對含錳巖系進行了旋回地層學方面的研究，揭示了旋回層序在下三疊統北泗組含錳巖系的發育規律和特點、旋回疊加模式以及層序成因機制，并提出了研究區錳礦成因模式。旨為今后在探討含錳巖系成因機制，建立勘探區內中、下三疊統高分辨率地層格架，以及應用定量模擬方法解決地層和沉積方面的問題等方面研究提供新的技術和思路。

1 區域地質背景

桂西南錳礦區位于中國南方大陸構造域與古特提斯—喜馬拉雅構造域交匯部位附近，早古生代屬于揚子準地臺與華南地槽的過渡區[12]。該區一級大地構造單元處于南華準地臺西南段，二級構造單元為右江再生地槽區[13]。研究區東南緣界限為南寧斷裂，東北緣界限為南丹—昆侖關斷裂，西北緣界限為彌勒—師宗斷裂，南緣界限紅河斷裂，在平面上表現為由東向西呈菱形狀展布的構造區，也被稱作為右江褶皺構造帶[14](圖1)。該區經歷了加里東構造運動期地槽階段—海西構造運動期地臺階段—印支構造運動期再生地槽階段的發展演化歷程。

桂西南地區普遍缺失上三疊統地層，發育以陸源碎屑巖沉積為主的中三疊統百逢組及以碳酸鹽巖為主的下三疊統地層北泗組與馬腳嶺組。其中，下三疊統北泗組是本區重要的含礦地層之一，含錳巖系主要分布于天等、田東、德保三縣交界部位，為近碳酸鹽臺地的淺海臺盆相碳酸鹽巖沉積。東部天等一帶主要發育灰—深灰色薄層泥灰巖、硅質泥灰巖，并含多層含錳硅質泥灰巖及碳酸錳礦層，地層總厚度182 m，含錳礦層厚約15 m，向北向西層數減少厚度變薄[15]。西部大旺一帶，含錳地層頂部為淺色薄層頁巖，泥巖及粉砂巖，夾有含錳灰巖，上部為細砂巖，中下部為內碎屑灰巖及鮞粒灰巖。

2 旋回層序劃分

2.1 原理與方法

地層中符合采樣密度要求的地球化學數據(如穩定碳、氧同位素、87Sr/86Sr比值和CaCO3含量等)和地球物理數據(如自然伽馬、磁化率、巖石密度和色率等)均可用于旋回地層學的研究[16]。在地下鉆井剖面中，自然伽馬測井曲線(GR)具有等間距采樣、連續測量和縱向分辨率高的特點，自然伽馬強度可以敏感地反映巖性旋回和巖相韻律互層在深度域中的變化[17-19]，進而能夠反映古環境和古氣候的變化[20-21]，可作為檢測旋回層序及識別海平面或湖平面變化趨勢的資料[22-23]，因而在旋回地層研究中得到了廣泛的應用[24-26]。

圖1 桂西南構造綱要圖(據文獻[15]修改)Fig.1 Structure units of the study area

自然伽馬測井曲線可以敏感的反映陸源碎屑巖—碳酸鹽巖混合沉積巖中黏土礦物含量的變化。在海侵期，由于相對海平面的上升，海域擴張，泥巖含量增加，因此所沉積的巖石中放射性黏土礦物如伊利石、高嶺石等含量較多，致使自然伽馬測井數值較高；而在海退期，隨著相對海平面的下降和海岸線逐漸向海遷移，沉積物中粉砂巖與碳酸鹽巖的含量增加，黏土礦物的含量逐漸減少，導致自然伽馬測井數值較低[27]。因此，可根據對泥質含量反映較為敏感的自然伽馬測井數值的高低變化的規律來劃分沉積地層中海侵—海退沉積旋回序列，并且沉積序列轉換界面的位置可作為識別沉積層序界面的標志。

在旋回地層研究中，殘差分析是最簡單常用的方法之一[28]。其通過計算參數值與參數平均值之差的累計值所提供的信息，進行參數數據序列中所包含的周期性信號分析。本文應用自然伽馬濾波曲線數據進行殘差運算，其變量為自然伽馬幅度殘差，即自然伽馬濾波數值與其平均值之差。其計算過程與應用在潮坪碳酸鹽巖沉積層序研究中的圖解相類似，而差別是Fischer圖解是采用旋回厚度作為計算參數[29-31]。應用測井曲線數據進行旋回幅度殘差計算，再利用波形法解析向上變深或向上變淺的沉積旋回疊置形式可以較好的劃分旋回層序。其計算方法如下：

(3) 對各測井參數數據殘差Yn進行累計計算，得到新數列Zn(n=1,2,3…)。Zn=Y1+Y2+Y3+…+Yn-1。

繪制自然伽馬幅度累計殘差圖解時，其橫坐標為地層深度，縱坐標為自然伽馬濾波數值殘差的累積偏移量，殘差累計數值的平均值則作為偏移數據的基準線。

2.2 測井數據選擇

為了對含錳巖系下三疊統北泗組層序地層格架有一個全面的了解，本文選用東平礦區鉆遇中、下三疊統地層較完整的ZK0302井的自然伽馬測井曲線作為旋回地層分析的參數曲線，該井所鉆遇的地層從頂至底分別為中三疊統百逢組，下三疊統北泗組和馬腳嶺組，鉆井深度為310 m，測井采樣間隔為0.05 m。其中，中三疊統百逢組未見頂，下三疊統馬腳嶺組未見底，而含礦地層下三疊統北泗組頂底界齊全。根據鉆孔ZK0302編錄資料，中三疊統百逢組位于井深10～84.71 m處，地層厚度為74.71 m，自然伽馬測井數值范圍為12.13～73.36 API，巖性主要以泥巖與泥質粉砂巖為主，底部以出現凝灰巖與下伏北泗組分界。含錳礦地層下三疊統北泗組約在井深84.71～250.58 m處，地層厚度為165.87 m，自然伽馬測井數值范圍為5.62～92.31 API，其巖性主要以硅質泥灰巖、含錳硅質泥灰巖夾碳酸錳礦為主。ZK0302井鉆遇的錳礦層較多，共為五層。該井北泗組上部只出現一層碳酸錳礦層，其厚度較薄，位于井深107.63～108.23 m處，厚度約為0.6 m。而北泗組下部錳礦層出現較為密集，共為四層，主要集中于井深150～170 m處，其分別出現于156.11～157.21 m、158.56～159.49 m、161.43～163.76 m、166.07～167.4 m處，厚度分別為1.1 m、0.93 m、2.3 m和1.33 m。下三疊統馬腳嶺組以出現條帶狀泥質灰巖為標志，與上覆下三疊統北泗組分界。其井深位于250.58～310.8 m處，地層厚度約為60.22 m，自然伽馬測井數值范圍為2.1～23.49 API。就測井曲線整體而言，自然伽馬測井數值的高值與地層中含泥質含量較高的巖性相對應，低值與含泥質含量較低的巖性相對應，自然伽馬曲線本身顯示出較好的旋回性變化。

2.3 測井頻譜分析

在實際情況下，測井信號在采集和測量的過程中，會受到大量噪聲信號的干擾，導致所測曲線形成較多尖峰和突變部分，造成含有較多的高頻干擾信號。并且在干井段和鉆井液侵染井段以及井壁受泥粉混染、鉆井液密度變化等干擾因素的影響，會造成測井曲線發生基線偏移，形成各種低頻波信號(圖2A，B)。因此在對實際測井數據進行分析之前，為消除噪音和人為干擾因素的影響，首先要進行測井數據的預處理，消除測井曲線中所夾雜的干擾信號。

本文通過去趨勢方法對原始數據進行預處理，去除原始測井數據的低頻信號的干擾。首先利用Past軟件[32-33]對原始自然伽馬測井曲線進行去趨勢處理，調用Past軟件下Model窗口下Polynomial程序，選擇需要濾掉的趨勢線階數，本文所有曲線均選用去除5階趨勢線的方法處理。之后再利用Matlab軟件，采用高通濾波方法，提取波長小于100 m(即頻率為0.01)的濾波數據進行分析。從圖中可見，經過處理的自然伽馬測井曲線數值既有正值，亦有負值(圖2C)。

圖2 東平礦區ZK0302井自然伽馬原始曲線、濾波曲線與累計殘差曲線及頻譜分析圖Fig.2 Original GR data，filtering GR data，cumulative residual curves and their spectral analysis diagram from the Well ZK0302 of Dongping manganese mine-zone

原始自然伽馬測井曲線進行信號預處理之后，需要對所得的自然伽馬測井濾波曲線數據進行頻譜分析。利用Past軟件進行頻譜分析，調用Time窗口下Spectral analysis程序，可得到濾波曲線中各頻率與功率的數據結果。再根據波長與頻率之間的倒數關系進行換算，即可知自然伽馬測井濾波曲線中各波長功率譜關系圖。

通過對ZK0302鉆孔自然伽馬濾波曲線數據進行頻譜分析(如圖2D)，其結果表明，地層厚度為56 m、33.42 m和11.64 m的主旋回波長為優勢頻譜峰，但高通濾波曲線仍受一系列波長小于十米的高頻波信號干擾。進一步進行累計殘差計算(圖2E)，并根據頻譜分析結果可知(圖2F)，該曲線僅保留了地層厚度為56 m與33.42 m的主旋回波長，波長小于十米的高頻波均無功率顯示。

2.4 旋回層序劃分

根據自然伽馬幅值累積殘差曲線的判別原理，ZK0302井所鉆遇的中、下三疊統地層共可劃分為七個海侵—海退沉積序列(圖3，4)。而中、下三疊統三個地層界面在累積殘差曲線上出現于井深85 m與250 m處負半波與正半波曲線的轉換界面位置附近，并且這個界面也是向上變淺與向上變深的沉積層序序列過渡的位置：下三疊統北泗組與頂部中三疊統百逢組界線在深度域上由深至淺變現為海侵沉積序列向海退沉積序列轉換的位置，為一正半波曲線向負半波曲線過渡的形態；而下三疊統北泗組與底部下三疊統馬腳嶺組地層界線在深度域上由淺至深表現為海退沉積序列向海侵沉積序列轉換的位置，其地層界線為負半波曲線與正半波曲線的轉換界面。含錳巖系下三疊統北泗組共可劃分為四個海進—海退旋回層序，錳礦層發育在Cs4與Cs6中的海退沉積序列內。同時可以看出，在北泗組地層內部還存在至少4～5個高頻旋回層序。

圖4 ZK0302鉆孔旋回層序綜合柱狀圖Fig.4 Analysis of cyclostratigraphy of the Well ZK0302 in Dongping manganese mine-zone

2.5 旋回結構解析

在以泥巖和灰巖為主的碳酸鹽巖地層剖面中，按照巖層的組成結構和厚度變化劃分，旋回地層可以識別出兩類變化趨勢(圖5)。一類是以泥質含量的變化作為指標，地層旋回可劃分為向上變深的海侵沉積序列或向上變淺的海退沉積序列。如果泥質含量以自然伽馬數值的變化作為指標，則可以通過解析沉積地層記錄中向上變淺和向上變深的旋回疊加形式，并可以確定旋回地層層序的轉換界面。另一類則是以地層的厚度作為參數，表現為向上變薄或變厚的旋回疊置樣式。這類分析也可以根據自然伽馬累積殘差曲線中正、負半波信號轉換點的位置計算波長求取旋回厚度在深度域上的變化，進而能夠反映沉積地層記錄中旋回層序厚度向上變薄或向上變厚的旋回疊加樣式的變化趨勢。

這兩類變化趨勢可以是同步變化的，也可以是無明顯相關變化的。一般解釋認為，碳酸鹽巖地層的組成結構變化與沉積環境的水動力強弱或相對海平面的變化有關，而巖層厚度的變化則受沉積物供應量大小以及沉積速率變化所控制[34-38]。因此，在旋回地層學的分析過程中，即可以通過測井信號幅度的變化規律劃分海侵—海退旋回層序序列，也可以根據地層旋回厚度的變化分析海侵—海退旋回中，各沉積序列沉積物供應量與沉積速率的變化趨勢。

圖5 沉積記錄中沉積旋回信號幅度和波長變化形式示意圖(據文獻[28]修改)Fig.5 Amplitude and wavelength directional changes of sedimentary cycle signals in stratigraphic records

對ZK0302井自然伽馬測井曲線進行高通濾波處理，濾波頻率的上限為0.01(即波長為100 m)，繪制其累積幅度殘差曲線。根據濾波曲線數據編繪的累積幅度殘差曲線上正負半波信號轉換點的位置計算其波長，可以求取并繪制深度域上每一個海侵沉積序列或海退沉積序列的厚度(圖6)。

觀察圖6中旋回波長隨井深變化可以發現，含錳礦地層下三疊統北泗組沉積序列厚度普遍大于中三疊統百逢組和下三疊統馬腳嶺組。同時根據計算求得下三疊統北泗組沉積序列厚度最大值為31.125 m，出現在井深172.25 m處，而旋回厚度最小值為1.25 m，位于井深203.375 m處。北泗組沉積序列厚度普遍較大，但有兩個低值區，位于井深100.875～118 m和203.375～212.125 m處，這些旋回厚度減薄的低值區間是沉積物供應量低，相對沉積速率小的標志。

通過ZK0302井自然伽馬測井曲線進行濾波頻率的上限為0.01的高通濾波處理計算得到的累積幅度殘差曲線中可知，在下三疊統北泗組地層內部還存在至少4～5個次一級的高頻旋回層序。由于頻率與波長存在倒數關系，因此增加濾波頻率的上限，即減小濾波波長的長度，可更為細致的研究地層旋回中沉積序列厚度的變化規律。

由前文所述，ZK0302鉆孔自然伽馬濾波曲線數據頻譜分析結果中，還存在主旋回波長為33.42 m的優勢頻譜峰，因此，繼續對自然伽馬測井數據進行濾波處理，將濾波上限增加到0.3，用所得到的濾波數據再進行累積幅度殘差計算，可進一步研究旋回層序厚度在更高頻率段上向上變薄或變厚的旋回疊置關系的變化趨勢。從旋回波長隨井深變化的曲線中可以觀察到(圖7)，下三疊統北泗組與馬腳嶺組的沉積序列厚度總體大于中三疊統百逢組。同時，下三疊統北泗組的頂部約85 m處與上覆中三疊統百逢組和底部約250 m處與下伏下三疊統馬腳嶺組的地層界限附近，有兩個較為明顯的旋回地層厚度變薄的地層界面，而這兩個界面深度恰是海侵沉積序列與海退沉積序列過渡的部位，其可能反映了沉積物的供應量由減少到增加、沉積速率從下降到上升的沉積過程。

對比上述兩種不同頻率的高通濾波曲線所得到旋回波長變化圖的可以發現，濾波頻率較低的波長變化圖反映的是低頻長旋回地層厚度的變化規律，可用來解析沉積地層記錄中旋回厚度向上變厚和向上變薄變化的整體趨勢。而濾波頻率較高的旋回波長變化圖可以表現出高頻短旋回地層厚度的變化特征，這能夠對旋回地層沉積序列厚度的變化規律進行細節上的研究。同時可以看出，濾波頻率越高，則對旋回厚度向上變厚和向上變薄變化規律描繪的越為細致，對沉積物供應量與沉積速率變化特征的研究也越為準確，并且可更加精確的標定地層界面。

圖6 鉆孔ZK0302自然伽馬累計殘差曲線及波長變化圖(濾波上限0.01)Fig.6 Cumulative residual curves and cyclic wavelength variation diagram of GR filtering data from the Well ZK0302 of Dongping manganese mine-zone(filter upper limit 0.01)

圖7 鉆孔ZK0302自然伽馬累計殘差曲線及波長變化圖(濾波上限0.3)Fig.7 Cumulative residual curves and cyclic wavelength variation diagram of GR filtering data from the Well ZK0302 of Dongping manganese mine-zone(filter upper limit 0.3)

鉆孔ZK0302鉆遇的碳酸錳礦層主要位于下三疊統北泗組頂部井深約107.63～108.23 m處和底部井深約156～167 m處。頂部碳酸錳礦層厚度較薄，厚度約為0.6 m。底部錳礦層發育較為密集，共為四層，總厚度約5.66 m。同步分析自然伽馬測井曲線累計殘差曲線和旋回波長變化曲線可以發現，北泗組中發育的碳酸錳礦層在自然伽馬累計殘差曲線上出現于海退旋回層序內，而在旋回波長變化曲線上位于旋回波長較大、地層向上變厚的相對沉積速率高值區。這些旋回厚度增厚的高值區表明沉積物供應量高、沉積速率大有利于碳酸錳礦的發育。并且通過進一步研究可知，在下三疊統北泗組地層內部，碳酸錳礦出現的底部層位旋回地層厚度明顯較頂部層位大，這說明在碳酸錳礦沉積期，北泗組底部沉積物供應量較高，沉積速率也較大，因此北泗組底部所發育的碳酸錳礦層也較頂部礦層數量多，礦層厚度也較大。

2.6 聯井旋回層序劃分與對比

通過上述方法，對鉆孔ZK1102、ZK7601與ZK7213自然伽馬測井曲線進行處理，并建立起研究區旋回層序地層格架后發現(圖8)，旋回層序的組成特征以及層序格架內部含錳巖系的橫向分布特點都具有規律性，能夠合理的發現碳酸錳礦層發育的有序性：1)含錳巖系下三疊統北泗組與頂部中三疊統百逢組以及底部下三疊統馬腳嶺組地層界線在自然伽馬幅值累積殘差曲線深度域上由淺至深均表現為由負半波向正半波轉換的曲線形態，即為海退—海侵沉積序列轉換界面位置附近。2)含錳巖系下三疊統北泗組可劃分為四個海進—海退旋回層序，大多表現為對稱狀或者近似對稱狀旋回，這表明北泗組至少經歷了四次相對海平面變化，且整個研究區所處的臺盆相沉積環境相似且穩定。并且在某些旋回層序內部，由于高頻海平面變化的影響，鉆孔自然伽馬幅值累積殘差曲線波動較大，這說明在地層內部還存在多個次一級的高頻旋回層序。3)碳酸錳礦層發育的旋回層序位置不同，但主要出現于各旋回層序內海退沉積序列中，并位于旋回層序內海退—海侵沉積序列界面附近，這表明錳礦的沉積主要發生在相對海平面由緩慢下降至逐漸上升的時期的轉換時期。

3 錳礦沉積成因討論

晚二疊世至中三疊世，在印支構造運動拉張作用的持續影響下，使桂西南地區在晚古生代形成的準地臺逐漸消失，致使研究區進入了裂谷坳陷期而形成右江再生地槽。在此階段，晚古生代末已存在的裂陷帶開始急劇擴張，使碳酸鹽臺地發生裂解而破碎，形成了大面積的拗陷區[39-40]。在早三疊世，在構造拉張作用下因同沉積斷裂帶發育了眾多的地勢較低的溝槽或盆地，形成了較為典型的臺盆相間的沉積古地理格局，這對之后錳礦的錳質的來源、沉積環境以及成礦作用提供了有利條件。

圖8 研究區中、下三疊統旋回地層劃分與對比圖Fig.8 The Triassic sequence stratigraphic framework in the southwestern Guangxi area

研究區下三疊統北泗組含錳巖系的發育經歷了多次海平面變化，錳礦層主要發育于旋回層序格架的海退沉積序列中，并位于海退沉積序列與海侵沉積序列之間的轉換界面附近。據張舜新[41]對本區田東縣作登剖面與平果縣太平圩剖面的牙形石帶的時代鑒定，本區下三疊統北泗組時代屬晚三疊世奧倫尼克階(245～249.7 Ma)。通過ZK0302鉆孔累計殘差曲線可進一步作出反映其低頻長趨勢海平面變化曲線，與全球海平面對比后發現[42]，其與全球三級海平面變化具有有較好的同步性，這說明研究區區海平面變化主要受全球三級海平面的影響(圖4)。同時可以發現，早三疊世在全球三級海平面上升的大背景下，研究區錳礦層的發育受更高一級的海平面海退期所控制。

自然水體系統的氧化還原環境決定了錳元素的溶解度及賦存狀態(圖9)。錳元素的價態會隨著環境的不同而發生變化：在氧化環境下，錳離子會呈高價態的Mn4+賦存于錳氧化物中(如軟錳礦)；而在還原環境中，錳離子會呈低價態的Mn2+溶解于水中或賦存于錳碳酸鹽中(如菱錳礦)。在分層海水中，錳質在水體深部的還原環境中以Mn2+離子狀態穩定存在，使錳在海水中的保持較高溶解度。當水體中氧化還原界面降低時，使原先處于還原狀態的Mn2+離子進入富氧的狀態，而被氧化成高價態的Mn4+離子形成氧化物而發生沉淀；而當水體中氧化還原界面再次升高時，含錳的氧化物重新進入還原環境發生溶解或與碳酸氫鹽發生反應形成錳碳酸鹽的沉積[43-45]。因此，氧化還原界面的波動控制著錳質賦存形態的性質，并決定了錳礦的發育。

圖9 氧化錳與碳酸錳在自然水體中的穩定性 (據文獻[45]修改)Fig.9 Stability of Mn oxide and Mn carbonate deposits in natural water

海平面變化影響著氧化還原界面的變化，因此在沉積錳礦床的形成過程中起著重要的作用。含錳巖系中錳含量的高值與海平面上升期最低含氧帶擴展有關，錳質的富集出現在層序地層格架中海侵體系域與高水位體系域轉換期內，并且大多與黑色頁巖共生[46]。在海平面高水位期，海灣和河口彎是錳質沉積產生的重要場所，結核狀和團塊狀的錳礦受到潮汐、河流的搬運而富集或分散。在海侵期海平面上升產生缺氧環境、沉積速率降低有利于錳的沉淀聚集。如果氧化還原界面深度隨著海平面升降發生變化，當處于還原狀態的海水底層低價態游離Mn2+隨海平面的上升混入到處于表層氧化狀態的海水中時，溶解錳就會被氧化而發生沉淀[47-48]。由此可知，海平面變化是通過控制氧化還原界面的波動來間接的影響碳酸錳礦形成，即盆地內氧化錳與碳酸錳的轉化環境以及碳酸錳成礦后的賦存層位主要受由海平面變化引起的氧化還原界面的升降所控制。

目前，現代大洋沉積物還沒有錳碳酸鹽沉積的報道，現代大洋底部的錳結核和錳結殼一般為氧化錳礦物形式出現，但在地質記錄中錳碳酸鹽礦床分布較為廣泛。氧化錳和碳酸錳的沉淀出現在缺氧或貧氧盆地環境，碳酸錳形成在早期成巖階段，它是在有機質的氧化過程發生氧化錳還原反應的結果，并且錳的富集與細菌活動密切或與海底甲烷氣滲漏事件有成因聯系[49-50]。通過綜合考慮研究區錳礦形成的海平面變化、古地理條件、錳質來源、氧化還原條件等特征的基礎上，筆者提出了研究區錳礦成因模式，即具有地理隔絕形式的，并由同生斷裂控制的“局限海”臺盆相碳酸錳礦沉積模式，研究區碳酸錳礦的形成主要經歷以下三個階段(圖10)：

(1) 海平面下降—氧化錳礦物形成階段

在海平面水位較低的階段，來自海岸區域的陸源地表水大量地注入海水，此時錳質除了來源于臺盆內生的熱液噴口之外，還有大量溶解在河水以及有機質中的錳質膠體與離子提供。海平面下降導致氧化還原界面下降，使臺盆內之前處于還原環境的深部區域逐漸轉化為水體較淺的氧化環境，并且由于地表水的注入，使臺盆內水體循環增強，富氧的水層與淡化海水中的的錳質大量混合，導致處于海水表層溶解狀態的錳質以較穩定的高價態氧化錳礦物(如MnO2)的形式大量沉積下來形成氧化錳礦物(圖10A)。

(2) 海平面上升—氧化錳礦物轉化碳酸錳礦物階段

在經歷低海平面階段后，海平面逐漸上升，使海岸區域的陸源地表流入海水的注入量逐漸減少。在這個階段，受控于海平面變化影響的氧化還原界面逐漸上升，之前形成的氧化錳礦物帶由水體較淺的氧化環境也逐漸過渡于水體較深的低氧甚至缺氧的還原環境。并且由于地表水的流入，海水中有機質的含量增加，導致藻類以及各類浮游類生物大量繁殖，導致海水中二氧化碳含量逐漸增加，貧氧的還原帶控制區域也逐漸擴大。同時，在水體中賦存的大量碳酸根、碳酸氫根離子與高價態的氧化錳礦物發生反應，成巖的過程中在貧氧的還原環境中不斷轉化為錳碳酸鹽，造成在旋回層序轉換界面附近碳酸錳礦層的形成。若海平面上升速度較緩慢，則會使氧化錳礦物溶解深度大，并與水體的反應也較為充分，碳酸錳礦層發育的厚度也會較大(圖10B)。

圖10 研究區碳酸錳礦沉積模式圖Fig.10 Schematic representation of Mn deposition in the southwestern Guangxi area

(3) 最大海泛面—少量碳酸錳礦物發育階段

在經歷海平面上升、水體加深階段之后，便進入了最高海泛面階段。在此階段，注入海體內的陸源地表流水大量減少，海體中的錳質主要由臺盆內生的熱液噴口提供。由于水體較深，氧化還原界面較高，在沒有充足的氧氣提供下，錳質在還原環境下主要以低價態的Mn2+離子形式溶解在海水中。但由于高水位期浮游生物在富氧層大量繁殖，有機質產率高，導致貧氧層逐漸擴大，在臺盆邊緣處于富氧的氧化環境逐漸轉化為貧氧的還原環境，因此也會導致少量厚度較薄、錳質含量較低的錳礦層發育(圖10C)。

4 結論

(1) 研究區含錳巖系下三疊統北泗組與頂部中三疊統百逢組以及底部下三疊統馬腳嶺組地層界線在自然伽馬幅值累積殘差曲線深度域上由淺至深均表現為由負半波向正半波轉換的曲線形態，即為海退沉積序列與海侵沉積序列轉換界面位置附近。北泗組內部共可劃分為四個海進—海退旋回層序，大多表現為對稱狀或者近似對稱狀旋回，這表明北泗組至少經歷了四次相對海平面變化，且整個研究區所處的臺盆相沉積環境相似且穩定。

(2) 海平面變化主要是通過控制氧化還原界面的波動來間接的影響碳酸錳礦形成。研究區碳酸錳礦層發育的旋回層序位置不同，但主要出現于各旋回層序內海退沉積序列中，并位于旋回層序內海退—海侵沉積序列界面附近，這表明錳礦的沉積主要發生在相對海平面由緩慢下降至逐漸上升的時期的轉換時期。早三疊世在全球三級海平面上升的大背景下，研究區錳礦層的發育受更高一級的海平面海退期所控制。研究區碳酸錳礦的形成主要經歷以下3個階段：1)海平面下降—氧化錳礦物形成階段；2)海平面上升—氧化錳礦物轉化碳酸錳礦物階段；3)最高海泛面—少量碳酸錳礦物發育階段。

致謝 感謝貴刊編輯部和審稿專家對文章進行審閱并提出寶貴意見。

)

[1] 秦元奎，徐柏安，姚敬劬. 桂西南上泥盆統含錳巖系形成古海水深度條件分析[J]. 地質與勘探，2010，46(2)：277-284.[Qin Yuankui, Xu Bai’an, Yao Jingqu. Depth of ancient sea related with formation of the upper Devonian Mn-bearing rock series in southwestern Guangxi province[J]. Geology and Exploration, 2010, 46(2): 277-284.]

[2] 秦元奎，張華成，姚敬劬. 廣西大新縣下雷錳礦床的地球化學特征及其意義[J]. 地質論評，2010，56(5)：664-672. [Qin Yuankui, Zhang Huacheng, Yao Jingqu. Geochemical characteristics and geological implication of the Xialei manganese deposit, Daxin county, Guangxi[J]. Geological Review, 2010, 56(5): 664-672.]

[3] 王榮庚. 廣西大新下雷錳礦沉積環境特征及成礦物質來源分析[J]. 中國錳業，2013，30(4)：8-10，21. [Wang Ronggeng. An elementary analysis of Xialei Mn-ore with the sedimentary features and material sources of ore-forming[J]. China’s Manganese Industry, 2013, 30(4): 8-10, 21.]

[4] 祝壽泉. 東平紅土型錳礦的地質特征[J]. 中國錳業，1998，16(1)：9-13，22. [Zhu Shouquan. Geological features of lateritic type manganese deposit in Dongping[J]. China’s Manganese Industry, 1998, 16(1): 9-13, 22.]

[5] 祝壽泉. 廣西東平錳礦半氧化帶中的菱錳礦[J]. 地質與勘探，2001，37(2)：58-61. [Zhu Shouquan. Phodochrosite of hemioxidative zone in Dongping manganese ore[J]. Geology and Prospecting, 2001, 37(2): 58-61.]

[6] 駱華寶，我國優質錳礦的勘查方向[J]. 地質與勘探，2002，38(4)：8-11. [Luo Huabao. The direction of prospecting for high quality manganese ore resources in China[J]. Geology and Prospecting, 2002, 38(4): 8-11.]

[7] 張海峰，張林曄，姚益民，等. 東營凹陷沙河街組四段灘壩砂米氏旋回層序劃分與對比[J]. 地層學雜志，2008，32(3)：278-284. [Zhang Haifeng, Zhang Linye, Yao Yimin, et al. Milankovitch cycles and sequence division and correlation of the member 4 of the Shahejie Formation in the Dongying depression, Shandong[J]. Journal of Stratigraphy, 2008, 32(3): 278-284.]

[8] Osleger D, Read J F. Relation of eustasy to stacking patterns of meter-scale carbonate cycles, Late Cambrian, U.S.A.[J]. Journal of Sedimentary Research, 1991, 61(7): 1225-1252.

[9] 李慶謀. 測井曲線Milankovitch周期分析與應用[J]. 地球物理學報，1996，39(5)：699-704. [Li Qingmei. The analysis and application of Milankovitch cycles by logging data[J]. Acta Geophysica Sinica, 1996, 39(5): 699-704.]

[10] 黃春菊. 旋回地層學和天文年代學及其在中生代的研究現狀[J]. 地學前緣，2014，21(2)：48-66. [Huang Chunju. The current status of cyclostratigraphy and astrochronology in the Mesozoic[J]. Earth Science Frontiers, 2014, 21(2): 48-66.]

[11] 尹青，伊海生，夏國清，等. 基于測井曲線頻譜分析在倫坡拉盆地古近系米氏旋回層序及可容空間變化趨勢中的研究[J]. 地球物理學進展，2015，30(3)：1288-1297. [Yin Qing, Yi Haisheng, Xia Quoqing, et al. Accommodation space and Milankovitch orbit cycle sequence of the Paleocene stratigraphic frames in Lunpola Basin based on the spectrum analysis of the logging curve[J]. Progress in Geophysics, 2015, 30(3): 1288-1297.]

[12] 毛健全，顧尚義，張啟厚. 右江—南盤江裂谷構造格局[J]. 貴州地質，1999，16(3)：188-194. [Mao Jianquan, Gu Shangyi, Zhang Qihou. Framework of fault tectonics in Youjiang-Nanpanjiang Rift[J]. Guizhou Geology, 1999, 16(3): 188-194.]

[13] 歐莉華. 桂西南地區上泥盆統錳礦沉積特征與成礦機理研究[D]. 成都：成都理工大學，2013. [Ou Lihua. Research on depositional feature and mineralization mechanism of the upper Devonian manganese deposit in southwestern Guangxi area[D]. Chengdu: Chengdu University of Technology, 2013.]

[14] 侯明才，陳洪德，田景春，等. 泥盆紀右江盆地演化與層序充填響應[J]. 地層學雜志，2005，29(1)：62-70. [Hou Mingcai, Chen Hongde, Tian Jingchun. Devonian sequence filling response to the evolution of the Youjiang Basin[J]. Journal of Stratigraphy, 2005, 29(1): 62-70.]

[15] 李升福，王澤華，李朗田，等. 桂西南優質錳礦成礦機理分析[J]. 資源環境與工程，2009，23(4)：363-370. [Li Shengfu, Wang Zehua, Li Langtian, et al. Analysis of metallogenic mechanism of high-grade manganese ore in Southwest Guangxi[J]. Resources Environment & Engineering, 2009, 23(4): 363-370.]

[16] Abels H A, Aziz H A, Krijgsman W, et al. Long-period eccentricity control on sedimentary sequences in the continental Madrid Basin (middle Miocene, Spain)[J]. Earth and Planetary Science Letters, 2010, 289(1/2): 220-231.

[17] 陳茂山. 測井資料的兩種深度域頻譜分析方法及在層序地層學研究中的應用[J]. 石油地球物理勘探，1999，34(1)：57-64. [Chen Maoshan. Two novel depth-domain frequency spectrum analysis methods for logging data and their application to sequence stratigraphy research[J]. Oil Geophysical Prospecting, 1999, 34(1): 57-64.]

[18] 李鳳杰，王多云，程微. 應用自然伽馬曲線反演隴東地區延安組沉積旋回[J]. 成都理工大學學報(自然科學版)，2004，31(5)：473-477. [Li Fengjie, Wang Duoyun, Cheng Wei. Use of natural gamma-ray well to study sedimentary cycles of Yan’an Formation in eastern Gansu, Ordos Basin, China[J]. Journal of Chengdu University of Technology (Science & Technology Edition), 2004, 31(5): 473-477.]

[19] 劉紅歧，陳平，夏宏泉. 測井沉積微相自動識別與應用[J]. 測井技術，2006，30(3)：233-236. [Liu Hongqi, Chen Ping, Xia Hongquan. Automatic identification of sedimentary microfacies with log data and its application[J]. Well Logging Technology, 2006, 30(3): 233-236.]

[20] Wonik T. Gamma-Ray measurements in the Kirchrode Ⅰ and Ⅱ boreholes[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2001, 174(1/2/3): 97-105.

[21] 陳中紅，查明，金強. 自然伽瑪及自然伽瑪能譜測井在沉積盆地古環境反演中的應用[J]. 地球物理學報，2004，47(6)：1145-1150. [Chen Zhonghong, Zha Ming, Jin Qiang. Application of natural gamma ray logging and natural gamma spectrometry logging to recovering paleoenvironment of sedimentary basin[J]. Chinese Journal of Geophysics, 2004, 47(6): 1145-1150.]

[22] 鄭興平，周進高，吳興寧. 碳酸鹽巖高頻層序定量分析技術及其應用[J]. 中國石油勘探，2004，9(5)：26-30. [Zheng Xingping, Zhou Jingao, Wu Xingning. High-frequency sequence quantitative analysis technology of carbonate rock and its application[J]. China Petroleum Exploration, 2004, 9(5): 26-30.]

[23] 徐麗萍. 測井約束地震反演技術在塔河油田碎屑巖儲層預測中的應用[J]. 石油物探，2009，48(6)：591-595. [Xu Liping. Application of log constrained seismic inversion for clastic reservoir prediction in Tahe oilfield[J]. Geophysical Prospecting for Petroleum, 2009, 48(6): 591-595.]

[24] 張小會，趙重遠. 鄂爾多斯盆地上三疊統延長組米蘭科維奇旋回的確定[J]. 石油與天然氣地質，2002，23(4)：372-375. [Zhang Xiaohui, Zhao Zhongyuan. Definition of Milankovitch cycles for Yangchang Formation of the upper Triassic in Ordos Basin[J]. Oil & Gas Geology, 2002, 23(4): 372-375.]

[25] 鄭興平，羅平. 川東渝北飛仙關組的米蘭克維奇周期及其應用[J]. 天然氣勘探與開發，2004，27(1)：16-19. [Zheng Xingping, Luo Ping. Analysis and application of Milankovitch cycles on Feixianguan Formation, northeast Sichuan Basin, China[J]. Gas Exploration & Development, 2004, 27(1): 16-19.]

[26] 朱劍兵. 基于測井資料的高分辨率層序旋回劃分方法探討[J]. 油氣地球物理，2011，9(4)：6-11. [Zhu Jianbin. The discussion on the division methods of high resolution sequence cycle based on logging data[J]. Petroleum Geophysics, 2011, 9(4): 6-11.]

[27] 伊海生. 測井曲線旋回分析在碳酸鹽巖層序地層研究中的應用[J]. 古地理學報，2011，13(4)：456-466. [Yi Haisheng. Application of well log cycle analysis in studies of sequence stratigraphy of carbonate rocks[J]. Journal of Palaeogeography, 2011, 13(4): 456-466.]

[28] 伊海生. 沉積旋回疊置形式的波形分析及旋回層序劃分方法[J]. 沉積學報，2015，33(5)：855-864. [Yi Haisheng. The Waveform graphic analysis of cyclic stacking patterns in sedimentary successions and detection methods of cyclostratigraphic sequences[J]. Acta Sedimentologica Sinica, 2015, 33(5): 855-864.]

[29] Read J F, Goldhammer R K. Use of Fischer plots to define third-order sea-level curves in Ordovician peritidal cyclic carbonates, Appalachians[J]. Geology, 1988, 16(10): 895-899.

[30] Goldhammer R K, Dunn P A, Hardie L A. Depositional cycles, composite sea-level changes, cycle stacking patterns, and the hierarchy of stratigraphic forcing: examples from Alpine Triassic platform carbonates[J]. Geological Society of America Bulletin, 1990, 102(5): 535-562.

[31] Husinec A, Basch D, Rose B, et al. FISCHERPLOTS: an excel spreadsheet for computing Fischer plots of accommodation change in cyclic carbonate successions in both the time and depth domains[J]. Computers & Geosciences, 2008, 34(3): 269-277.

[32] Hammer ?, Harper D A T, Ryan P D. Past: paleontological statistics software package for education and data analysis[J]. Palaeontologia Electronica, 2001, 4(1): 1-9.

[33] Hammer ?, Harper D A T. Paleontological Data Analysis[M]. Hoboken, United States: Wiley, 2007.

[34] 伊海生. 地層記錄中旋回層序界面的識別方法及原理[J]. 沉積學報，2012，30(6)：991-998. [Yi Haisheng. Detection of cyclostratigraphic sequence surfaces in stratigraphic record: its principle and approach[J]. Acta Sedimentologica Sinica, 2012, 30(6): 991-998.]

[35] Groetsch J. Cycle stacking and long-term sea-level history in the lower Cretaceous (Gavrovo Platform, NW Greece)[J]. Journal of Sedimentary Research, 1996, 66(4): 723-736.

[36] Hinnov L A, Park J. Detection of astronomical cycles in the stratigraphic record by frequency modulation (FM) analysis[J]. Journal of Sedimentary Research, 1998, 68(4): 524-539.

[37] Cunha A A S, Koutsoukos E A M. Orbital cyclicity in a Turonian sequence of the Cotinguiba Formation, Sergipe Basin, NE Brazil[J]. Cretaceous Research, 2001, 22(5): 529-548.

[38] Rodíguez-Tovar F J, Pardo-Igúzquiza E. Strong evidence of high-frequency (sub-Milankovitch) orbital forcing by amplitude modulation of Milankovitch signals[J]. Earth and Planetary Science Letters, 2003, 210(1/2): 179-189.

[39] 曾允孚，劉文均，陳洪德，等. 華南右江復合盆地的沉積構造演化[J]. 地質學報，1995，69(2)：113-124. [Zeng Yunfu, Liu Wenjun, Chen Hongde, et al. Evolution of sedimentation and tectonics of the Youjiang composite basin, South China[J]. Acta Geologica Sinica, 1995, 69(2): 113-124.]

[40] 鐘怡江，陳洪德，侯明才，等. 右江盆地北部中晚泥盆世層序充填特征分析[J]. 中國地質，2011，38(2)：338-345. [Zhong Yijiang, Chen Hongde, Hou Mingcai, et al. Sequence filling features of middle-late Devonian period in northern Youjiang Basin[J]. Geology in China, 2011, 38(2): 338-345.]

[41] 張舜新. 桂西下三疊統牙形石序列的新認識[J]. 現代地質，1990，4(2)：1-15. [Zhang Shunxin. On the lower Triassic conodont sequence of western Guangxi[J]. Geoscience, 1990, 4(2): 1-15.]

[42] Haq B U, Hardenbol J, Vail P R. Chronology of fluctuating sea levels since the Triassic[J]. Science, 1987, 235(4793): 1156-1167.

[43] Hem J D. Redox processes at surfaces of manganese oxide and their effects on aqueous metal ions[J]. Chemical Geology, 1978, 21(3/4): 199-218.

[44] Hem J D. Rates of manganese oxidation in aqueous systems[J]. Geochimica et Cosmochimica Acta, 1981, 45(8): 1369-1374.

[45] Roy S. Sedimentary manganese metallogenesis in response to the evolution of the Earth system[J]. Earth-Science Reviews, 2006, 77(4): 273-305.

[46] Jarvis I, Murphy A M, Gale A S. Geochemistry of pelagic and hemipelagic sediments: criteria for identifying systems tracts and sea-level change[J]. Journal of the Geological Society, 2001, 158(4): 685-696.

[47] Frakes L A, Bolton B R. Origin of manganese giants: sea-level change and anoxic-oxic history[J]. Geology, 1984, 12(2): 83-86.

[48] 唐云鳳，伊海生. 滇東南地區斗南沉積型錳礦床礦物相變化及沉積模式[J]. 中國地質，2011，38(2)：451-461. [Tang Yunfeng, Yi Haisheng. Mineral phase changes and depositional model of sedimentary manganese deposits in Dounan area of southeastern Yunnan[J]. Geology in China, 2011, 38(2): 451-461.]

[49] Polgari M, Okita P M, Hein J R. Stable isotope evidence for the origin of therkút manganese ore deposit, Hungary[J]. Journal of Sedimentary Research, 1991, 61(3): 384-393.

[50] Okita P M, Shanks III W C. Origin of stratiform sediment-hosted manganese carbonate ore deposits: examples from Molango, Mexico, and TaoJiang, China[J]. Chemical Geology, 1992, 99(1/2/3): 139-163.

ResearchonCyclostratigraphicandDepositionalModelofSedimentaryManganeseDeposits(Triassic)inSouthwesternGuangxiArea

YIN Qing1,2, YI HaiSheng1,3, ZHOU Kai2

1.InstituteofSedimentaryGeology,ChengduUniversityofTechnology,Chengdu610059,China2.ResearchInstituteofExplorationandDevelopmengt,SINOPECZhongyuanOilfieldCompany,Puyang,Henan457001,China3.StateKeyLaboratoryofOil/GasReservoirGeologyandExploitation,Chengdu610059,China

The manganese deposits are mainly distributed in Southwestern Guangxi area, which was regarded as one of the most predominant ore types in this area. However, the Triassic stratigraphic frameworks of manganese deposit in the southwestern Guangxi area is rarely discussed. The spectral analysis and filter methods can be adopted to divide cyclic sequence hierarchy and discuss depositional mechanism, and the cumulative residual curves from average amplitude variations of good logs are introduced and demonstrated as a valid graphical method to define depositional regime transformation and check boundary surfaces of cyclic sequences. The cumulative residual curves results show the positive half waves and negative half waves of cumulative residual curves are, respectively, corresponds to upward-deepening and upward-shallowing sedimentary successions in low frequency cyclic sequences in a well interval. Transform points of positive half waves and negative half waves are calibrated to boundary surface depth of cyclic sequences. The GR log data from Well 0302 of Southwestern Guangxi Area as a case example, cumulative residual curves are calculated and constructed on well log data. The low and high frequency cyclic sequences are distinguished and its procedure and methods are illustrated in detail. The significance and interpretation of amplitudes and wavelength changes of sedimentary cycles are discussed in this paper. The Beisi Formation of the Lower Triassic can be divided into four forth-order sequences, and the manganese ore deposit developed in the regression cycle. The origin model of manganese ore can be divided in three stages: 1) Sea-level fall-the oxidizing manganese is formation; 2) Sea-level rise-the oxidizing manganese is transformed in manganese carbonate; 3) The maximum flooding surface-a little manganese carbonate is formed. The investigation provides a new effective approach to build the correlating framework of high-resolution sequence in Southwestern Guangxi Area.

southwestern Guangxi area; Triassic; manganese-bearing sequence; cyclostratigraphy; sea level change

1000-0550(2017)06-1144-14

10.14027/j.cnki.cjxb.2017.06.006

2016-01-05；收修改稿日期2016-11-17

“十二五”國家科技支撐計劃項目(2011BAB04B10)；國家自然科學基金項目(40972084)；高等學校博士點專項科研基金課題(20105122110012)[Foundation“Twelfth five-year” National Science and Technology Plan, No. 2011BAB04B10; The National Natural Science Foundation of China, No.40972084; Specialized Research Fund for the Doctoral Program of Higher Education, No.20105122110012]

尹青，男，1986年出生，博士后，助理研究員，沉積地質學研究，E-mail: 168347532@qq.com

伊海生，男，教授，E-mail: yhs@cdut.edu.cn

TE121.3

A