董火祥, 劉景彥*, 李勇, 史超群
(1.中國(guó)地質(zhì)大學(xué)(北京)能源學(xué)院, 北京 100083; 2.中國(guó)石油塔里木油田分公司勘探開(kāi)發(fā)研究院, 庫(kù)爾勒 841000)
塔里木盆地是中國(guó)主要含油氣盆地之一,三疊系作為塔里木盆地主要勘探目的層系之一[1-2],在早期的勘探中陸續(xù)發(fā)現(xiàn)了多個(gè)油氣藏[3],表明三疊系陸相地層對(duì)塔里木盆地具有十分重要的地位。但是由于復(fù)雜的構(gòu)造環(huán)境導(dǎo)致對(duì)塔里木三疊系的沉積環(huán)境認(rèn)識(shí)并不清晰。前人在對(duì)三疊系進(jìn)行研究探討時(shí),從層序地層學(xué)的角度出發(fā),得到的結(jié)論存在較大差異。
傅恒等[4]將三疊系劃分為下三疊1個(gè)長(zhǎng)周期,中上統(tǒng)均為3個(gè)長(zhǎng)周期,共7個(gè)長(zhǎng)周期,對(duì)應(yīng)于三級(jí)層序;呂雪雁等[5]將塔中塔北區(qū)帶三疊系劃分下三疊3個(gè)三級(jí)層序,中上統(tǒng)各1個(gè)三級(jí)層序;郭建華等[6]的劃分方案為下三疊1個(gè),中三疊3個(gè),上三疊2個(gè)三級(jí)層序。受限于資料品質(zhì)及采用的層序理論不同,諸多學(xué)者對(duì)三疊系的層序劃分方案存在較大差異[7]。層序方案劃分存在爭(zhēng)議,導(dǎo)致沉積相帶發(fā)育同樣存在差異性,制約了研究區(qū)的勘探進(jìn)展。因而對(duì)層序發(fā)育情況進(jìn)行客觀(guān)深入研究,有助于相帶研究突破。當(dāng)?shù)卣鹳Y料品質(zhì)較差,按照傳統(tǒng)的地震追蹤不整合面來(lái)識(shí)別三級(jí)層序具有一定的困難,導(dǎo)致劃分層序不一致時(shí),如何客觀(guān)地對(duì)連續(xù)沉積的地層劃分層序是一個(gè)關(guān)鍵問(wèn)題。現(xiàn)用小波分析的方法來(lái)對(duì)缺乏地震資料進(jìn)行層序劃分的地層進(jìn)行多級(jí)層序劃分,詳細(xì)刻畫(huà)塔里木三疊系各構(gòu)造單元內(nèi)層序發(fā)育情況,為研究區(qū)沉積相帶的研究及后續(xù)的生儲(chǔ)組合研究奠定一定的基礎(chǔ)。
由于測(cè)井?dāng)?shù)據(jù)蘊(yùn)含著豐富的地質(zhì)信息,利用其進(jìn)行小波變換來(lái)得到多級(jí)次旋回曲線(xiàn)有助于對(duì)內(nèi)藏的信息進(jìn)行揭示[8]。此外通過(guò)測(cè)井?dāng)?shù)據(jù)獲取的高頻旋回還能夠作為繪制Fischer圖解的基礎(chǔ),用來(lái)建立研究區(qū)的相對(duì)湖平面變化。
湖平面升降變化控制著各類(lèi)沉積體的分布情況和層序地層界面的形成,對(duì)于沉積體系分布及生儲(chǔ)蓋組合研究有著重要意義。在陸相湖盆中,沉積基準(zhǔn)面可以近似視為相對(duì)湖平面變化。近年來(lái),前人通過(guò)Fischer圖解曲線(xiàn)來(lái)刻畫(huà)相對(duì)湖平面升降,取得了諸多成果[9-10]。現(xiàn)綜合多口井的Fischer圖解,描繪可容納空間的變化情況,重建塔里木盆地三疊系相對(duì)湖平面變化曲線(xiàn),并與沉積相帶的變遷進(jìn)行對(duì)應(yīng)分析。研究結(jié)果將對(duì)高頻層序地層的研究和沉積相帶發(fā)育和有利儲(chǔ)層的預(yù)測(cè)具有一定的意義。
塔里木盆地是中國(guó)內(nèi)陸最大的含油氣盆地,總面積約56萬(wàn)km2,近似菱形,東西向展布。盆地處于阿爾金山脈、西昆侖山和天山山脈之間。三疊系時(shí)期,古特提斯洋向北俯沖,造成塔里木盆地西部隆起,沉積僅分布在坳陷內(nèi),即現(xiàn)今盆地的中北部[11-13]。庫(kù)車(chē)地區(qū)受沖斷負(fù)荷的作用發(fā)生大幅撓曲沉降,形成庫(kù)車(chē)前陸盆地。而中部盆地具有復(fù)合成因機(jī)制,不僅受到熱冷卻沉降,還受到新和前緣隆起的沖斷作用的影響,發(fā)生了撓曲沉降[14-15]。因此,該時(shí)期的原型盆地具有前陸盆地特征,發(fā)育有3個(gè)一級(jí)盆地:庫(kù)車(chē)前陸盆地、中部克拉通坳陷盆地及塔西南前陸盆地。但現(xiàn)今殘存的三疊系主要集中在庫(kù)車(chē)盆地及塔北-塔中盆地[12]。塔里木盆地現(xiàn)今劃分為庫(kù)車(chē)坳陷、塔北隆起、北部坳陷、中央隆起等7個(gè)一級(jí)構(gòu)造單元帶[16-17],研究區(qū)主要位于盆地的中北部,三疊紀(jì)時(shí)期,尚未形成如今的隆坳地貌格局,表現(xiàn)為東、西隆,中間坳,如圖1所示。
塔里木盆地三疊系地層自下而上可以劃分為俄霍布拉克組、克拉瑪依組、黃山街組和塔里奇克組,但在臺(tái)盆區(qū)塔里奇克組普遍缺失,有學(xué)者指出可能是構(gòu)造運(yùn)動(dòng)導(dǎo)致,但仍然存在爭(zhēng)議,具體成因尚不明確[13]。
臺(tái)盆區(qū)三疊系俄霍布拉克組發(fā)育湖泊-辮狀河三角洲相帶[5],在近物源的盆地邊緣處,該組底部發(fā)育粗砂巖、含礫砂巖等,與下部二疊系多呈不整合接觸;克拉瑪依組時(shí)期辮狀河三角洲、扇三角洲廣泛發(fā)育[18-19],近源粗碎屑巖發(fā)育有兩到三期,呈向上變細(xì)的特征;黃山街組多發(fā)育灰色、褐色泥巖,局部發(fā)育辮狀河三角洲[20]。由于三疊系末構(gòu)造運(yùn)動(dòng)形成大范圍剝蝕導(dǎo)致黃山街組殘余面積較小,與上部侏羅系多呈不整合接觸。臺(tái)盆區(qū)三疊系綜合柱狀圖如圖2所示。
圖1 塔里木盆地構(gòu)造單元及研究區(qū)井位示意圖Fig.1 Schematic diagram of structural units and well locations in the study area of Tarim Basin
圖2 臺(tái)盆區(qū)三疊系綜合柱狀圖Fig.2 Comprehensive histogram of Triassic in Platform Basin area
測(cè)井?dāng)?shù)據(jù)通常有著分辨率較高、連續(xù)性較好的優(yōu)點(diǎn),還具有豐富的地質(zhì)信息。當(dāng)?shù)卣鹳Y料分辨率較低、品質(zhì)較差,取芯資料缺乏時(shí),進(jìn)行高精度層序單元?jiǎng)澐窒鄬?duì)受限,對(duì)測(cè)井信息進(jìn)行挖掘能起到幫助作用。而測(cè)井小波分析方法有“數(shù)學(xué)顯微鏡”之稱(chēng),對(duì)數(shù)據(jù)變化的響應(yīng)精度高,可以定量識(shí)別各種級(jí)別的層序界面,從而彌補(bǔ)了依照傳統(tǒng)高分辨率層序劃分的不足[21-23]。
小波分析的基本原理是利用一個(gè)小波函數(shù)作為基本單位進(jìn)行伸縮和平移變換來(lái)模擬信號(hào)值,將信號(hào)值分解為多個(gè)不同尺度、多級(jí)次的信號(hào)[24-25],而在測(cè)井小波分析中,分解出的小波信號(hào)呈現(xiàn)周期性變化的特點(diǎn),與各個(gè)級(jí)次的層序具有較好的對(duì)應(yīng)關(guān)系。低頻小波信號(hào)一般與長(zhǎng)期旋回具有對(duì)應(yīng)關(guān)系,而高頻信號(hào)通常與高頻短期旋回對(duì)應(yīng),據(jù)此可以對(duì)各個(gè)級(jí)別的層序旋回進(jìn)行劃分[26-27]。
隨著小波分析方法的推廣及發(fā)展,應(yīng)用的小波函數(shù)也出現(xiàn)了諸多類(lèi)型,如Daubechies小波、Morlet小波等?;谇叭说难芯砍晒鸞28-30]及研究區(qū)的具體情況,經(jīng)過(guò)實(shí)驗(yàn)分析,最終選擇了效果最好的Daubechies小波及Dmeyer小波進(jìn)行一維連續(xù)小波變換。而進(jìn)行變換的測(cè)井信號(hào)選擇了伽馬曲線(xiàn),因?yàn)槠鋵?duì)砂泥巖變化敏感,能很好地揭示巖性變化[31]。
共選取了12口單井作為研究對(duì)象,且均勻分布在4個(gè)一級(jí)構(gòu)造帶,每個(gè)構(gòu)造帶均有3口單井,分布位置如圖1所示。以A2井為例,對(duì)其進(jìn)行小波變換及高頻層序劃分。當(dāng)選用設(shè)定階數(shù)為10,最大級(jí)數(shù)為12的小波系數(shù)曲線(xiàn),能夠較好地識(shí)別各級(jí)層序,且相互對(duì)應(yīng)。如圖3所示,對(duì)A2井測(cè)井伽馬值進(jìn)行一維連續(xù)小波變換后,根據(jù)db小波系數(shù)曲線(xiàn)中的d11、d10、d7級(jí)小波曲線(xiàn)(結(jié)合Dmeyer11級(jí)小波曲線(xiàn)),在目的層段中識(shí)別出4個(gè)長(zhǎng)周期旋回、8個(gè)中周期旋回及57個(gè)短周期旋回。研究區(qū)內(nèi)12口單井經(jīng)過(guò)小波變換劃分旋回結(jié)果如圖4所示。其中,塔北隆起帶中的A1~A3井長(zhǎng)、中、短周期旋回發(fā)育平均值分別為4、8、58;北部坳陷帶的B1~B3井長(zhǎng)、中、短周期旋回發(fā)育平均值分別為3、5、40;中央隆起帶C1~C3井長(zhǎng)、中、短周期旋回發(fā)育平均值分別為3、4、40;西南坳陷帶的D1~D3井長(zhǎng)、中、短周期旋回發(fā)育平均值分別為2、4、33。長(zhǎng)、中周期旋回?cái)?shù)比值在1.3~2范圍內(nèi),長(zhǎng)、短周期旋回?cái)?shù)比值在13~16.5。
從數(shù)據(jù)上來(lái)看,4個(gè)區(qū)帶旋回?cái)?shù)量從多到少排列為:塔北隆起帶,北部坳陷帶,中央隆起帶,西南坳陷帶??梢钥闯?塔北隆起區(qū)帶的地層旋回個(gè)數(shù)保留最多,說(shuō)明該區(qū)帶內(nèi)原始沉積地層保留相對(duì)完整,其次是北部坳陷帶及中央隆起帶,最差的是西南坳陷帶。
從層序角度分析,長(zhǎng)周期旋回大致對(duì)應(yīng)于三級(jí)層序,中周期旋回發(fā)育于三級(jí)層序內(nèi)部。從A2-B2-C1-D1連井層序格架圖(圖5)上可見(jiàn),從南到北,三級(jí)層序個(gè)數(shù)漸增,從1個(gè)增加至4個(gè),A2井層序、旋回?cái)?shù)最多,D1井層序、旋回?cái)?shù)最少。從垂向上來(lái)看,底部的SQ1發(fā)育有1~3個(gè)中期旋回;SQ2發(fā)育有2個(gè)中期旋回;SQ3與SQ4均發(fā)育有1個(gè)中期旋回。
綜上分析,A1~A3井所位于的塔北隆起帶旋回結(jié)構(gòu)發(fā)育最為完整,從底到頂,僅黃山街組存在剝蝕,可能導(dǎo)致旋回缺失,而其他區(qū)帶內(nèi)井旋回相對(duì)缺失,且缺失程度從北向南加劇。因此以塔北隆起帶井旋回?cái)?shù)據(jù)來(lái)分析原始沉積旋回比較恰當(dāng)。
借助Morlet小波變換來(lái)進(jìn)行能譜分析能夠?qū)⒁痪S時(shí)間域測(cè)井信號(hào)變換到二維時(shí)頻域,而不同的沉積旋回具有不同的時(shí)頻變換形態(tài),當(dāng)測(cè)井曲線(xiàn)的形態(tài)變化不明顯時(shí),通過(guò)對(duì)時(shí)頻變換形態(tài)組合的分析,可以快速地從測(cè)井資料中提取時(shí)頻特征變化規(guī)律,輔助研究層序旋回的變化及沉積環(huán)境的改變[32]。但分析時(shí)需要注意排除首尾處異常高能量團(tuán)的影響,這是由于軟件本身所導(dǎo)致的干擾因素[33]。Morlet小波的能量團(tuán)變化形態(tài)與層序旋回的對(duì)應(yīng)關(guān)系如圖6所示。退積型、加積型與進(jìn)積型砂體形態(tài)對(duì)應(yīng)的Morlet小波的能量團(tuán)的分布與遷移有規(guī)律性變化[34]。沉積地層呈現(xiàn)進(jìn)積型沉積樣式時(shí),砂體厚度向上減薄,沉積時(shí)水體能量減弱,能譜分析圖中能量團(tuán)垂向上能量減弱,橫向尺度上由低頻向高頻遷移;呈現(xiàn)加積沉積樣式時(shí),砂體垂向上厚度均一,沉積時(shí)水體能量不變,能量團(tuán)垂向上能量不增不減,橫向尺度不變;呈現(xiàn)退積沉積樣式時(shí),砂體垂向上厚度增大,沉積時(shí)水體能量增強(qiáng),能量團(tuán)垂向上能量增強(qiáng),橫向尺度上由低頻向高頻遷移。
根據(jù)db小波綜合Dmeyer小波對(duì)A2井進(jìn)行層序劃分的結(jié)果,對(duì)4個(gè)三級(jí)旋回分別進(jìn)行Morlet小波變換能譜分析,來(lái)判斷沉積時(shí)期水體能量強(qiáng)弱的變化,結(jié)果如圖7所示。每個(gè)層序中都存在多個(gè)水體環(huán)境能量變化旋回,SQ1時(shí)期在中部位置,水體能量達(dá)到最弱,相對(duì)應(yīng)的湖平面最高;SQ2早期水體能量相對(duì)較弱,相對(duì)湖平面處于高位,中晚期水體能量較強(qiáng);SQ3早期、晚期水體能量均較強(qiáng),相對(duì)湖平面處于低位,中期水體能量最弱,此時(shí)湖平面最高;SQ4中期水體能量達(dá)到最弱,相對(duì)湖平面達(dá)到高值。
Fischer圖解最先是在研究奧地利環(huán)潮坪碳酸鹽沉積時(shí)提出[35],后經(jīng)過(guò)發(fā)展完善,逐步演化為現(xiàn)今常用的形式。以旋回?cái)?shù)量、平均厚度偏移累積量分別作為橫軸與縱軸,以此得到Fischer圖解曲線(xiàn)[36-37]。早期的Fischer圖解主要應(yīng)用于海相碳酸鹽領(lǐng)域[38-40],后經(jīng)過(guò)不斷的發(fā)展,逐漸被推廣至陸相湖盆碎屑巖領(lǐng)域[10,30]。
在較穩(wěn)定的陸相湖盆中,湖平面的升降主要受可容納空間的控制,兩者的變化趨勢(shì)相一致,而Fischer圖解中平均厚度累積偏移曲線(xiàn)代表可容納空間的變化趨勢(shì),因此Fischer圖解曲線(xiàn)可以近似地表示相對(duì)湖平面的變化趨勢(shì)[41]。
將Fischer圖解運(yùn)用到研究區(qū)內(nèi)三疊系相對(duì)湖平面變化分析中,當(dāng)厚度累積偏移值正向偏移的特征時(shí),表明沉積地層厚度增大,是對(duì)湖平面上升導(dǎo)致可容納空間增長(zhǎng)的響應(yīng);當(dāng)厚度累積偏移值出現(xiàn)負(fù)向偏移時(shí),表明沉積地層厚度減小,是對(duì)湖平面下降導(dǎo)致可容納空間減小的響應(yīng)。
由于已對(duì)各區(qū)帶長(zhǎng)、中、短周期旋回進(jìn)行了劃分,依據(jù)劃分結(jié)果顯示,A區(qū)帶的高頻旋回保存最為完整,因此選取A1、A2、A3測(cè)井伽馬曲線(xiàn)進(jìn)行測(cè)井小波定量層序劃分。對(duì)所得到的高頻旋回經(jīng)過(guò)統(tǒng)計(jì)得到了高頻層序累積偏移量,并繪制了Fischer圖解(圖8)。由于缺乏校正地層厚度所需要的砂泥巖孔隙度的數(shù)據(jù),并且前人[41-43]認(rèn)為缺乏壓實(shí)校正過(guò)程不會(huì)對(duì)Fischer圖解曲線(xiàn)的趨勢(shì)形態(tài)及旋回組合產(chǎn)生影響,因此未進(jìn)行壓實(shí)校正。
①~⑤表示旋回結(jié)構(gòu)圖8 相對(duì)湖平面變化曲線(xiàn)Fig.8 Relative lake level change curve
Fischer圖解揭示出研究區(qū)內(nèi)的相對(duì)湖平面的升降規(guī)律,從圖8中可見(jiàn),A1井主體表現(xiàn)為5個(gè)旋回結(jié)構(gòu),第一個(gè)旋回僅發(fā)育湖退過(guò)程,第二至第四個(gè)旋回在可容納空間低值附近波動(dòng),表明當(dāng)時(shí)可容納空間較小,湖平面位于低值,伴隨著快速湖侵,第五個(gè)旋回主體發(fā)育在黃山街組內(nèi),同時(shí)包含了可容納空間最大值,這表明黃山街組時(shí)期可容納空間保持在較高水平。A2井發(fā)育有4個(gè)旋回,T2k發(fā)育有2個(gè)三級(jí)旋回,T1e和T3h內(nèi)部分別存在一個(gè)三級(jí)旋回,但A2井的可容納空間最大值點(diǎn)不同于A(yíng)1井,而是位于T2k時(shí)期內(nèi)。說(shuō)明同一個(gè)湖盆不同區(qū)帶的可容納空間的變化受到空間位置的影響,可能是由于物源供給所引起的可容納空間的不一致。A3井與A2井相似度極高,三級(jí)層序旋回分布相一致,可容納空間最低值點(diǎn)位于T2k時(shí)期內(nèi),最高值點(diǎn)位于T3h內(nèi)部。比較這3口井的可容納空間曲線(xiàn),發(fā)現(xiàn)A1井俄霍布拉克組存在一個(gè)旋回是A2、A3井沒(méi)有的,推測(cè)成因?yàn)锳2、A3井在俄霍布拉克組早期處于構(gòu)造高位,缺少地層沉積,因此未能記錄下第一個(gè)旋回。
將Fischer圖解與Morlet小波分析結(jié)果進(jìn)行對(duì)比,具有很高的一致性。整體來(lái)看,臺(tái)盆區(qū)三疊系相對(duì)湖平面呈現(xiàn)一次湖退到湖侵的過(guò)程,俄霍布拉克組與克拉瑪依組發(fā)育湖退,直至克拉瑪依中晚期達(dá)到相對(duì)湖平面最低值點(diǎn),克拉瑪依組晚期與黃山街組內(nèi)呈現(xiàn)湖侵的特點(diǎn),相對(duì)湖平面的最大值基本位于黃山街組內(nèi)部。
在A(yíng)1~A3井所在的塔北隆起帶內(nèi),選擇一條順物源方向的連井剖面進(jìn)行層序沉積對(duì)比分析。如圖9所示,A區(qū)共發(fā)育有5個(gè)三級(jí)層序,其中SQ1層序在以A1井附近的區(qū)域發(fā)育半深湖-深湖相帶,與一個(gè)長(zhǎng)周期旋回相對(duì)應(yīng),往東北方向沉積地層受到X井所在的英買(mǎi)力低凸的古地貌控制,無(wú)沉積發(fā)育。SQ1時(shí)期沉積相帶垂向上深湖-半深湖泥質(zhì)沉積過(guò)渡為半深湖灘壩及泥質(zhì)沉積,湖平面變化曲線(xiàn)主要體現(xiàn)為一個(gè)湖退的過(guò)程,與沉積物縱向變化特征一致。
SQ2發(fā)育在俄霍布拉克組中上段,這一時(shí)期,南北兩側(cè)無(wú)大型物源匯入,廣泛發(fā)育濱淺湖相,湖平面變化曲線(xiàn)在中高水位動(dòng)蕩,濱淺湖泥質(zhì)沉積與灘壩構(gòu)成主要沉積微相,湖平面變化曲線(xiàn)顯示,此時(shí)可容納空間小于SQ1時(shí)期,在中低水位波動(dòng)。
SQ3發(fā)育于克拉瑪依組下段,沉積相帶為辮狀河三角洲前緣-濱淺湖,水下分流河道、河口壩厚層砂體及遠(yuǎn)端灘壩發(fā)育,頂部發(fā)育有少量濱淺湖泥質(zhì)沉積,可容納空間主體呈水退的特征,在較低水位處動(dòng)蕩。
SQ4位于克拉瑪依組中上段,這一時(shí)期西部物源充足發(fā)育辮狀河三角洲前緣近端,相帶向東過(guò)度為三角洲前緣遠(yuǎn)端-濱淺湖。西部在大量陸源碎屑的注入下,湖平面變化曲線(xiàn)呈現(xiàn)水進(jìn)的特征,與錄井巖性呈現(xiàn)向上變細(xì)的正旋回相一致。東部物源較少,垂向相帶展布呈現(xiàn)三角洲前緣近端-濱淺湖-三角洲前緣遠(yuǎn)端的變遷。與湖平面變化曲線(xiàn)呈現(xiàn)的水進(jìn)-高位-水退的變化形態(tài)相吻合。
SQ5時(shí)期,黃山街底部存在東部推進(jìn)的大套三角洲前緣近端,向上過(guò)渡為半深湖-深湖的暗色泥質(zhì)沉積,而湖平面變化曲線(xiàn)呈現(xiàn)水進(jìn)-高位震蕩的特征,且可容納空間最大值位于黃山街組內(nèi),表明黃山街組沉積時(shí),湖平面位于較高水位。
將湖平面升降曲線(xiàn)與層序-沉積連井剖面進(jìn)行對(duì)比分析,發(fā)現(xiàn)湖平面的升降變化與沉積相帶縱向演變、橫向遷移相互契合,因此認(rèn)為本次得到的湖平面升降曲線(xiàn)能夠代表臺(tái)盆區(qū)三疊系的湖平面變化曲線(xiàn)。
(1)在相對(duì)穩(wěn)定的陸相湖盆中,借助小波分析可以有效識(shí)別各級(jí)層序旋回的變化。Morlet小波變換能譜分析圖能夠揭示出沉積時(shí)期水體能量的變化,Fischer圖解揭示出可容納空間的變化,精細(xì)地反映相對(duì)湖平面變化的趨勢(shì)。
(2)通過(guò)小波變換定量層序劃分在塔里木盆地三疊系臺(tái)盆區(qū)不同區(qū)帶識(shí)別出多個(gè)不同級(jí)別的高頻旋回,最完整區(qū)帶為塔北隆起帶,共識(shí)別出5個(gè)長(zhǎng)期層序、9個(gè)中期層序和61個(gè)短期層序。旋回保存完整性從高到低依次為塔北隆起帶、北部坳陷帶、中央隆起帶、西南坳陷帶。
(3)根據(jù)能譜分析圖可以劃分出沉積時(shí)期水體能量變化趨勢(shì)及沉積模式,能量團(tuán)減弱對(duì)應(yīng)湖平面上升;能量團(tuán)增強(qiáng),對(duì)應(yīng)湖平面下降。
(4)臺(tái)盆區(qū)三疊系湖平面變化以塔北隆起帶最具有代表性,整體呈現(xiàn)一次大規(guī)模的湖退到湖侵。俄霍布拉克組先緩慢湖退,直至克拉瑪依組中期,可容納空間達(dá)到最低值,此后緩緩湖侵,黃山街組時(shí)期達(dá)到可容納空間最大值,并一直在較高水位動(dòng)蕩。