呂 明 王 穎 徐 微

(中海油研究總院)

沉積模擬方法在Bonaparte盆地的應用

呂 明 王 穎 徐 微

(中海油研究總院)

在有一定的鉆井資料但地震資料分辨率不高的情況下,嘗試借助沉積模擬方法對澳大利亞Bonaparte盆地某研究區(qū)進行了平面沉積相編圖和儲層快速評價工作。根據(jù)鉆井巖心、巖性、古生物年齡、古環(huán)境等資料,在研究區(qū)Elang組—Frigate組識別出了3個層序界面,確定相對海平面呈脈動上升變化;通過類比分析,認為研究區(qū)沉積為水進型三角洲河口灣體系,主物源在盆地東南方向;輸入這些參數(shù)進行沉積模擬,得到了研究區(qū)巖性及沉積結構的時空展布特征,并以此指導平面沉積相圖的編制。研究區(qū)具有灣內潮控、灣外浪控,高位潮控、低位浪控,近源扇三角洲與相對遠源三角洲共存的復合型沉積特點。將沉積模擬方法用于儲層研究,是在特定條件下進行儲層快速評價的有效方法。

沉積模擬 沉積相編圖 儲層快速評價 Bonaparte盆地

Bonaparte盆地位于澳大利亞西北陸架區(qū)中段,面積約27萬km2(圖1)。研究區(qū)位于盆地西北部靠近陸架邊緣處,主體位于Sahu l向斜,東部為Sahul臺地,東南部為Petrel次盆,南部為Londonderry高地[1],面積約4 200 km2,水深大多在200 m以內;研究區(qū)內及周邊有不少鉆井,并有二維及部分三維地震資料覆蓋,現(xiàn)已發(fā)現(xiàn)油田和氣田各一個。

Bonaparte盆地是典型的古生界和中生界疊合裂谷盆地,具有由若干凹陷和隆起相間分布的構造格局,并且早期構造線呈NW向,晚古生代后轉為NE向[2]。盆地自下而上依次發(fā)育石炭系、二疊系、三疊系、侏羅系、白堊系和第三系[3]。前人的沉積研究多體現(xiàn)在區(qū)域上,認為早期為淺?！侵蕹练e,中、晚期為河流三角洲—邊緣海沉積[3-4],研究精度尚不能滿足目標評價的需要,且在物源方向、沉積類型上還存在分歧。

圖1 Bonaparte盆地研究區(qū)位置圖

雖然研究區(qū)鉆井資料較多,但現(xiàn)有的地震資料品質較差,且研究層段的厚度小,內部反射結構不清,難以用地震相進行平面沉積相編圖。因此,筆者探索用沉積模擬方法輔助進行研究區(qū)中生界上侏羅至下白堊統(tǒng)目的層段(Plover組—Flam ingo組)平面沉積相編圖及沉積模式、沉積演化分析,以期在有限的條件下為海外探區(qū)快速評價提供更可信的信息。本次研究中,利用巖心(層理、構造、遺跡化石)、古生物、巖石、測井、地震反射、構造與地形演化、水系及物源等資料,編制了22口井的層序及沉積相分析圖、連井對比圖以及區(qū)域地震層序及沉積相解釋剖面等大量基礎圖件,保證了沉積模擬中有關沉積模式、沉積類型、海平面變化、物源方向、沉積物供應量等參數(shù)的合理選取。

1 沉積模擬基本理論及工作方法

除了描述以外,現(xiàn)代地層和沉積學研究的另一主要任務是利用少量不規(guī)則分布的樣本來預測巖石類型、結構、分布、巖石物理性質等。然而,如果沒有客觀標準和方法來評估多種解釋的可能性或不確定性,那么這種預測的價值就不大。但是,對于是否存在這種客觀標準和方法,能否通過模擬的方法來進行評估,長期以來一直存在爭議。Bu rton等認為地層數(shù)值模擬在理論上是不可能的[5],因為地層資料的信息量不足和存在非唯一性,不同作用參數(shù)的組合會產生同樣的結果,即使根據(jù)傳統(tǒng)層序地層學理論的假設和原則進行這種反演也是不可能的[6]。

但是,Cross認為只有在下列5種情況下進行反演才是不可能的[6]:①作用參數(shù)不是獨立的(一個作用能替代另一個作用而產生同樣結果);②作用及響應體系是隨機的、不可預測的,或缺少基本的約束條件;③不同位置的地層信息沒有聯(lián)系;④正演模型很差不能模擬真正的作用及響應體系;⑤未正確測量或解釋觀察數(shù)據(jù)?；趯Φ貙有再|和沉積作用及響應體系的不同理解,Cross認為地層含有足夠信息量,用定量方法重建盆地歷史可行。這些不同的理解在于:①地層不僅是“大事件”的記錄,而且在多數(shù)情況下某些低幅、單向變化的作用要比高幅、間歇性、低頻、偶發(fā)、“災變性”的事件產生更顯著、更廣泛和更有意義的標志。②從四維地層的視點,時間是連續(xù)和完整的,在地層記錄中沒有缺失的時間,完整的時間以地層和地層不連續(xù)面或無沉積面為代表,巖石與面有相等的意義,無論何種成因的沉積間斷(侵蝕或無沉積)都攜帶有該時間段其它位置上的巖石實體的屬性信息。③地球上的各種作用不會突然全部停止。④沉積過程遵守物質守恒原則。⑤地層的作用及響應體系并非過于復雜而不能模擬,盡管有眾多的作用變量、復雜的反饋循環(huán)和非線性的相互關系,但地層卻并不復雜。復雜的地層實際上是高度有組織的,僅含有相對少的、規(guī)則的、重復的地層形態(tài)。地層的規(guī)則性是有組織的沉積形態(tài)與地層控制下的改造作用相結合的產物,這種改造減少了原始沉積的變化性。⑥地層作用參數(shù)的組合不能互相替代而產生完全相同的地層響應。這些不同的理解為使用模擬的方法重建沉積史,進行層序及沉積相的連井對比和平面分布特征的預測提供了科學基礎。

為此,Cross領導研究小組自上世紀80年代以來一直致力于層序及沉積相的計算機模擬工作,至90年代已成功地根據(jù)野外露頭獲得的海平面升降、構造沉降、沉積物供應等參數(shù)正演模擬出形態(tài)、厚度、相序等都與之相似的二維地質剖面。這些工作不僅驗證了高分辨率層序地層學的許多原理,而且顯示了計算機模擬方法強大的預測功能,展示了良好的應用前景,代表著層序地層學和沉積學研究的發(fā)展方向。近年來許多軟件公司也在著手開發(fā)這方面的軟件,現(xiàn)在的軟件可以根據(jù)鉆井、地震中獲得的海平面升降、構造沉降、沉積物性質和數(shù)量、物源方向、構造和地貌等參數(shù)正演或重建研究區(qū)三維的地層、巖性(砂巖或泥巖)分布及變化情況。

沉積模擬方法應用的關鍵主要體現(xiàn)在輸入參數(shù)選取和結果解釋及應用兩方面。筆者根據(jù)研究區(qū)鉆井資料多、地震資料信噪比低的現(xiàn)狀,提出了沉積模擬方法應用的技術路線(圖2):分析地質背景資料及前人研究成果,確定基本沉積類型和物源方向;結合古生物年齡資料進行單井基準面旋回和層序劃分、連井層序對比,結合區(qū)域地震剖面建立層序格架,確定沉積速率;根據(jù)古生物和區(qū)域地震剖面沉積中心變化情況分析相對海平面變化曲線;根據(jù)這些參數(shù)及相關構造圖進行沉積數(shù)值模擬;根據(jù)沉積模擬結果,結合等厚圖和地震屬性圖建立沉積模式,進行沉積相對比及平面沉積相編圖,完成儲層預測。

圖2 沉積模擬方法應用技術路線圖

2 沉積模擬參數(shù)分析

2.1 層序劃分與對比

參照Cross提出的基準面旋回識別模式[7],結合研究區(qū)地震剖面上的削蝕、上超等不整合層序界面的識別特征和鉆井在盆地中的相對位置信息,利用測井曲線(主要是自然伽馬)反映的地層疊加方式和巖性組合及突變,以及古生物缺帶及年齡等資料進行基準面旋回和層序劃分。在研究區(qū)22口井的Plover組—Flam ingo組中識別出了6個層序(表1),每個層序的延續(xù)時間在3～7M a,屬于中期基準面旋回,相當于Vail的三級層序[8]。

SB 1界面屬于中侏羅統(tǒng)Plover組頂界,為全區(qū)可追蹤的較大層序界面,區(qū)域地震剖面上可見對下伏層的明顯削蝕,具有約15 M a的沉積間斷(如Tam ar-1井)。SB4界面通常表現(xiàn)為由砂巖到泥巖的突變面,并普遍伴有約10 M a的沉積間斷(圖3)。SB5界面多與一套砂巖的底界面有關(如Iris-1井),在沒有砂巖的井中也可表現(xiàn)為由前積到退積的轉換,或泥巖顏色由深到淺的轉變(圖3)。

表1 Bonapar ter盆地層序關系表

基準面旋回的表現(xiàn)形式因鉆井位置而異。S3層序多表現(xiàn)為1個完整的中期旋回,且在凹陷部位或凸起帶上相對較低的鞍部,上升、下降半旋回轉換面處的深水泥巖較發(fā)育,反映了三角洲前緣或濱海灘壩與淺海為主的沉積特征(如圖3中Bu ller-1、Capung-1A井);而在向斜邊緣或凸起帶相對較高的部位,則泥巖多位于層序上部,使下降半旋回部分不發(fā)育或不明顯,反映了分流河道或潮道與濱海為主的沉積特征(如圖3中Buffalo-1、Krill-1井)。

古生物(主要用溝鞭藻)年齡資料為研究區(qū)層序的劃分與對比提供了較為可靠的依據(jù)。層序界面常表現(xiàn)為古生物年齡的間斷面,有時在大不整合面處由于再沉積的存在還會出現(xiàn)年齡反轉的現(xiàn)象,如Capung-1A、Buffalo-1、Krill等井在SB1附近上部測得年齡為169 M a,而下部測得為163或167 M a (圖3)。

總的來看,Plover組的年代跨度有40 M a,屬于Vail的二級層序[8]。Tam ar-1井揭露自下而上2個由深到淺的大旋回,包含4個三級層序,井中揭露了其中的3個。Elang組是一套砂巖為主的巖石地層單元,厚度100 m左右,年代跨度約5M a(160～165 M a)。該組內部具明顯旋回性,可分為S2、S3、S4等3個層序,其中S4層序的上升半旋回屬于Elang組,下降半旋回屬于以泥巖為主的Frigate組。S2與S3層序大致以163 M a為界,表現(xiàn)為2個較對稱的完整旋回,而且層序厚度薄,側向變化不大,反映當時盆地內地勢起伏不大,屬于構造運動相對和緩期。S4層序厚度側向變化增大,反映當時盆地內局部構造活動增強,凹陷處沉降加速,隆起處抬升強烈,造成地層減薄甚至缺失;該層序在盆內凸起及周緣以上升半旋回為主,在低洼處則保存有較多下降半旋回。Flam ingo組(分為S5、S6兩個層序)各井雖都以大套灰色泥巖為主,但地層厚度側向變化大,同樣反映當時盆地內有強烈局部構造活動。SB5界面上的濁積扇和SB6界面上15M a的沉積間斷都與較大的區(qū)域海平面下降有關(圖3)。

圖3 研究區(qū)連井層序對比剖面(剖面位置見圖11)

2.2 相對海平面分析

古生代至中生代晚期全球絕對海平面呈長期脈動上升趨勢,最高點出現(xiàn)在白堊系中部,此后轉為長期脈動下降[9](圖4左)。研究區(qū)的相對海平面變化與此不盡一致,可通過區(qū)域地震剖面和鉆井巖性及古生物信息得到相對海平面變化。由于受到信噪比的限制,在研究區(qū)難以用地震資料進行精細的層序分析,但對于大的層序界面還是可以通過削蝕、上超等特征幫助識別。另外,利用沿盆地傾向的長地震剖面進行區(qū)域層序解釋,通過沉降中心、沉積中心位置的變化也可分析研究區(qū)中—長期相對海平面變化(圖4)。

圖4示出的地震剖面中明顯可見3個大的構造旋回,大致對應于古生代克拉通內斷陷、中生代裂陷和新生代被動陸緣等3套沉積地層。從沉積地層厚度變化看,古生界二疊系的沉積中心在南部Petrel次盆;中生界三疊系厚度較小且未見明顯沉積中心,推測沉積中心北遷至研究區(qū)北部更遠處;侏羅系中部的沉積中心南移至Sahul向斜附近,上部地層厚度又變小,推測沉積中心曾再次北遷至現(xiàn)代陸架坡折之外;白堊系盆地邊緣的沉積厚度明顯加大,沉積中心又迅速南遷至Sahu l向斜以南;此后沉積中心隨陸坡的推進而逐漸向北遷移,但新生界中新統(tǒng)以上又有所南移。根據(jù)沉積物體積分配原理,推斷研究區(qū)的相對海平面變化趨勢是:二疊紀最高,三疊紀呈下降趨勢;侏羅紀早、中期再次上升,晚期先下降然后又迅速脈動上升,至白堊紀中期達到最高;此后隨陸架坡折的不斷推進海平面逐漸下降,僅新生代中新世后由于盆地中心沉積速率的增大而發(fā)生過局部性的海平面上升。從鉆井巖性特征看,Plover組向上泥質含量或泥巖夾層逐漸增多,至Flam ingo組則迅速變?yōu)榇筇啄鄮r,反映了侏羅紀末到白堊紀初沉積水深由緩慢增大到迅速增大的過程。

2.3 物源及沉積模式分析

據(jù)前人研究,Bonaparte盆地在中生代是一個由岡瓦納古陸邊緣破裂而形成的裂陷盆地,沉積中心在三疊紀偏北靠近蒂汶海,在侏羅紀南移至Sahul向斜、M alita地塹區(qū)[4]。由于古生代Petrel次盆的影響,盆地的南界向陸凹入約250 km,形成倒三角形或喇叭口狀的海灣地貌形態(tài)。此地形與現(xiàn)今的Bonaparte灣在位置和形態(tài)上基本一致,而且從區(qū)域地震剖面看,侏羅系與新生界都是一套由陸向海前積推進的地層(圖4),沉積格局相似,故適合用將今論古的方法研究目的層段的沉積特征。

圖4 Bonapar te盆地區(qū)域層序演化特征及相對海平面變化圖(區(qū)域地震剖面位置見圖11)

研究區(qū)衛(wèi)星圖片(圖5)顯示,現(xiàn)今該區(qū)的陸架較寬(約200 km),在陸架邊緣坡折以上有沉積物堆積,可見大致呈扇形或叢狀的水道形態(tài),反映出三角洲的特征。從水道的走向看,主要物源是來自有更大匯水范圍的河流體系,這些河流可沿低洼處匯集并經Londonderry高地兩側海灣處入海。而近岸局部山區(qū)的近源河水系也可作為局部物源形成小型三角洲。鑒于研究區(qū)現(xiàn)今為海灣且無長年性的遠源大河流注入,認為這些陸架上的三角洲是一般海退期的低位三角洲,由于物源區(qū)地勢低平,剝蝕量小,沉積物供應不充分,所以這種三角洲主體在陸架上發(fā)育,未形成典型的盆底扇,但大海退期的水道能以下切谷的形式切過這些三角洲及陸架坡折,如圖5中Bonaparte盆地區(qū)就見2條這種下切谷,分別向西北和東北方向延伸進入深海區(qū)。至于在下切谷的終端未見盆底扇發(fā)育,分析其原因是這2條下切谷都開口于蒂汶海溝,從海溝底沉積物的裸露和侵蝕情況判斷,這里有較強的沿岸海流或洋流的影響,會將沉積物向西搬運至海溝的開口之外沉積,如在B row se盆地深水盆底中的一套盆底扇堆積體與此有關。

圖5 Bonaparte區(qū)衛(wèi)星圖片

Bonaparte灣喇叭口狀的海岸形態(tài)會造成潮水的聚能作用,使潮差增大、潮流加強,同時較為寬闊的陸架緩坡會使波浪作用相對減弱,因此在正常情況下容易發(fā)育潮控型三角洲(理論上其特點是有港灣形的地貌,灣內三角洲在潮流的作用下被沖刷成若干垂直岸線分布的長條形潮流砂壩或砂脊),巴布亞灣三角洲就是此類潮控三角洲的典型代表,此外恒河-布拉馬普特拉河、伊洛瓦底江、湄公河以及珠江、長江、鴨綠江、遼河等三角洲都在一定程度上屬于潮控三角洲。研究區(qū)Elang組鉆井巖心中大量的潮汐層理、泄水構造和Skolithos、Ophiom orpha等生物遺跡化石也證實了潮汐作用的存在,如K ril-1井巖心中可見具槽狀及板狀交錯層理的三角洲分流河道、受潮汐改造的反粒序及脈狀層理的三角洲前緣河口壩、具羽狀交錯層理的潮汐砂壩、具波狀及透鏡狀層理的潮控河口灣和受到強烈生物擾動的潮坪等(圖6)。但是,從圖5示出的衛(wèi)星圖片中可以看到澳大利亞陸架邊緣的低位三角洲前緣是沿平行岸線方向呈條帶狀分布,這與陸架坡折附近較強的波浪作用和沿岸流改造有關。分析認為,在低海平面期,海灣的聚能作用減弱,而原來的陸架坡折區(qū)則成為三角洲前緣區(qū)使岸坡變陡,加之被狹窄的蒂汶海溝加強了的海流及沿岸流對岸坡的沖刷又進一步加劇了岸坡的變陡,三角洲的類型就變成浪控型三角洲了。浪控三角洲以平行岸線的沿岸砂壩為特點,砂巖分選好、側向連通性好,是理想的儲集體。

圖6 K r ill-1井巖心的潮控三角洲特征

綜合以上分析建立了研究區(qū)的沉積模式圖(圖7)。在低海平面期,發(fā)育浪控三角洲,波浪和沿岸流的改造作用使河口壩沿岸線呈條帶狀分布,形成厚度大且具有一定前積特征的三角洲前緣砂體;同時在風暴等突發(fā)事件的觸發(fā)下間歇性地發(fā)育以濁流為主的重力流沉積,在前三角洲或凹陷部位形成低位盆底扇或斜坡扇。在高海平面期,灣口處沉積物供應量減少,同時海水侵入地形平坦的原三角洲平原或海岸平原區(qū),水體的加深和岸坡的變緩使潮汐作用增強,形成了潮控三角洲體系,灣內的潮汐砂壩多垂直岸線分布,其間可為海灣泥巖所分隔。

圖7 研究區(qū)沉積模式圖

3 沉積模擬結果分析

以鉆井巖性、年代、古生物環(huán)境和主要地界面構造圖等資料為基礎,利用Dionisos軟件對澳大利亞Bonaparte盆地E lang組—Frigate組進行了沉積模擬。海平面變化曲線取向上脈動加深變化形式,給定東南方為主要物源方向,暫未進行壓實及剝蝕量等古構造、古地貌的恢復。

模擬結果可通過一系列平面圖和剖面圖以三維形式展示,這些圖件為連井對比和平面沉積相編圖提供了直觀的參考資料。從研究區(qū)沉積模擬平面圖(圖8)可見,在低海平面期(S2—S3層序,圖8中的a、b、c),研究區(qū)內砂質沉積物較為豐富,可覆蓋南部和中部的大部分地區(qū),這是由于此時受波浪影響較大,三角洲前緣砂體呈沿岸線的席狀展布。在高海平面期(S3—S4層序,圖8中的d、e、f),三角洲沉積中心向南遷移,研究區(qū)內砂質供應減少,Malita地塹部位出現(xiàn)泥質沉積,Sahul向斜部位仍有部分砂質充填,但多呈垂直岸線的條帶狀或網狀分布。

從研究區(qū)沉積模擬剖面圖(圖9)可見,前期(S2—S3層序)三角洲砂體主要由南向北沿Sahul向斜注入(圖9 f),后期(S3—S4層序)Londonderry、Tamar、Flamingo等高地對砂質沉積物的貢獻增大(圖9d),可在高地的周緣形成一些被潮汐改造、再搬運的砂體,成為Sahul向斜東翼很多構造目標的主要儲層類型。

4 平面沉積相分析

上述沉積模擬結果為研究區(qū)平面沉積相分析提供了新的信息,彌補了因地震相分辨率問題帶來的缺憾,這樣再結合鉆井沉積相解釋成果和參考地震屬性分析圖、構造圖、地層厚度圖等信息,可以編制平面沉積相圖,并進行研究區(qū)平面沉積相分析。

圖10 研究區(qū)S2層序底部均方根振幅圖

以研究區(qū)S2層序為例,其底部均方根振幅圖顯示由南向北振幅逐漸減弱(圖10),巖石地球物理相關性分析認為這反映了砂巖由南向北減少的變化趨勢,這種趨勢與鉆井揭示的砂巖厚度、粒度以及沉積模擬預測的砂巖百分比變化的結果(圖8)相一致。此外,從鉆井資料看,S2層序多為一套砂巖,僅粒度及砂巖含量略有變化(圖3)。然而,沉積模擬預測的砂巖分布及演化特征顯示,該層序除來自澳洲大陸遠源的三角洲外,還有來自工區(qū)內隆起或凸起帶近源扇三角洲的存在(圖8、9)。研究區(qū)鉆井資料與沉積模擬結果相結合而確定出的三角洲、扇三角洲及河口灣等沉積相帶的位置及三角洲分流河道的推進情況如圖11所示,可以看出,由于該層序沉積時的相對海平面較低,研究區(qū)內除Flamingo之外的幾個高地都出露水面而遭受剝蝕,在其周緣發(fā)育近源小型扇三角洲。同時,在研究區(qū)東南方向有一相對遠源的三角洲注入,受古地形的控制可形成幾個分流河道向不同方向散開,成為主要儲層發(fā)育的區(qū)域:第一支流向Malita地塹,由剖面圖可見砂巖主要來自南側(圖9 f),北面的Flamingo高地貢獻不大;第二支沿Flamingo高地西側的低地向北流,剖面上見向高地的上超(圖9e);第三支沿Sahul向斜向北流,為Laminaria高地南側的低凸起帶提供物源(圖9d);第四支沿Sahul向斜底部向西北流(圖9c),剖面顯示向北砂巖逐漸增多(圖9b、a),可作為生烴凹陷烴源巖的運移通道,其中水道可延伸至Fulica-1井附近,再向北則轉變?yōu)楹涌趬巍?/p>

圖11 研究區(qū)S2層序沉積相圖

總的來看,S2層序沉積時期研究區(qū)總體為三角洲-河口灣沉積體系,三角洲平原為含砂質的分流河道和偏泥質的分流間沼澤與泛濫平原,主要在Tamar高地東側分布;三角洲前緣以砂質為主,早期在波浪改造下可側向延伸成為河口壩復合體,晚期隨相對海平面上升潮汐作用增強,在向斜中多有垂直岸線的潮汐砂壩發(fā)育;前三角洲與盆地內淺海泥呈過渡關系。這種三角洲在演化上呈逐步向陸退縮的水進型特征,有利于蓋層的形成。

5 結束語

根據(jù)從鉆井或其它渠道提取的相關參數(shù)進行層序、沉積相及巖性體的三維數(shù)值模擬,進而確定沉積模式,指導平面沉積相編圖和儲層預測,這種方法有理論依據(jù),在技術上是可行的,特別是在資料有限而又需要及時得到盡可能準確、合理的評價結果的情況下進行沉積儲層快速評價更為有效。由于受條件的限制,本次研究中沒有進行研究區(qū)的古地貌恢復,同時軟件本身也存在一些不足,如物源方向只能選擇南-北、東-西方向,只能輸入統(tǒng)一的地層砂泥比,鉆井獲得的巖性資料不能在模擬時對其所在位置進行約束等,因此,沉積模擬的精度和可靠性均不夠,還不能作為沉積相和巖性解釋的直接依據(jù),但這畢竟是一次有益的嘗試。我們相信,隨著沉積研究向數(shù)字化、計算機化的進展,這種借助沉積模擬進行平面沉積相編圖和沉積儲層快速評價的工作方法會發(fā)揮越來越大的作用。

[1] BARRETTA G,H INDE A L,KENNARD JM.Undiscovered resource assessment methodologies and application to the Bonaparte Basin[G]//ELL ISG K,BA ILL IE PW,MUNSON T J.Timor Sea Petroleum Geoscience,Proceedings of the Timor Sea Symposium, Northern Territory,Australia,19-20 June,2003,Special Publication 1.Darwin:Northern Territory Geological Survey,2004:354-372.

[2] BA ILL IE P,FRASER T,HALL R,etal.Geological development of Eastern Indonesia and the northern Australia collision zone:a review[G]//ELL ISG K,BA ILL IE PW,MUNSON T J.Timor Sea Petroleum Geoscience,Proceedings of the Timor Sea Symposium, Northern Territory,Australia,19-20 June,2003,Special Publication 1.Darwin:Northern Territory Geo logical Survey,2004:539-550.

[3] AMBROSEG J.Jurassic sedimentation in the Bonaparte and northern Browse Basin:new models for reservoir-source rock development,hydrocarbon charge and entrapment[G]//ELL IS G K, BA ILL IE PW,MUNSON T J.Symposium on Timor Sea Petroleum Geoscience.Darw in:Bonaparte Basin and Surrounds,2004:125-142.

[4] LONGLEY IM,BUSSENSNSCHUETT C,CLYDSDALE L,et al. The North W est Shelf of Australia-a Woodside perspective[G]// KEEPM,MOSS S.The Sedimentary Basins of Western Australia v.3,Perth:Petroleum Exploration Society of Australia,2002:27-88.

[5]BURTON R,KENDALL C G C,LERCHE I.Out of our depth:on the impossibility of fathoming eustasy from the stratigraphic record [J].Earth-Science Review s,1987(24):237-277.

[6]CROSS T A.The scientific basis for stratigraphic inversion[G]. Ciencia Tecnologia Y Futuro,2000(1):824-845.

[7] CROSS T A.Sediment volume partitioning:rationale for stratigraphic model evaluation and high-resolution stratigraphic correlation[G]//SANDV IK K O,GRADSTE IN F,MILTON N.Predictive high resolution sequence stratigraphy.Norwegian Petroleum Society Special Publication,1998:171-196.

[8] VA IL PR,AUDEMARD E,BOWMAN SA,et al.The stratigraphic signatures of tectonics,eustasy,and sedimentology—an overview [G]//EINSELE G,R ICKENW,SEILACHERA.Cycles and events in stratigraphy.Berlin,Sp ringer-Verlag,1991:617-659.

[9] HAQ B U,HARDENBOL J,VA IL PR.The chrono logy of fluctuating sea level since the Triassic[J].Science,1987,235:1156-1167.

(編輯:崔護社)

An application of sedimentation simulation in Bonaparte basin

There are some d rilling data and low-resolution seismic data in an area of Bonaparte basin,Australia,where it is attempted to map sedimentary facies and quickly assess reservoir by means of the sedimentary simulation method.Based on the data such as d rilling core,lithology,palaeontology,and palaeoenvironment,three sequence boundaries in Elang to Frigate Formation were identifed,and a pulsed rise of relative sea level was determined in the area.An analogy has indicated that there is a transgrassive deposition system of deltaic estuary in the area,with its main provenance in the southeastern basin.These parameters were inputted to make a sedimentation simulation, resulting in the time-space distribution pattern of litho logy and sedimentary texture in the area,which candirectus to mapping sedimentary facies.As a composite deposition system,the area is characterized by tide dominance within the estuary,w ave dominance without the estuary,highst and tide dominance,lowstand wave dominance and coexistence of proximal fan delta and relative distal fan delta.Applying the sedimentary simulation method to reservoir rocks may be an efficient too l to quickly assess reservoirs under given conditions.

sedimentary simulation;mapping sedimentary facies;fast reservoir assessment;Bonaparte basin

呂明,男,高級工程師,1986年畢業(yè)于美國科羅拉多礦業(yè)大學地質與地質工程系,獲理科(地質)碩士學位,現(xiàn)任中海油研究總院勘探研究院沉積室主任,主要從事層序、沉積、儲層研究。地址:北京市東城區(qū)東直門外小街6號海油大廈409室(郵編:100027)。

2009-05-06 改回日期:2009-09-27