楊偉利，王毅，張眾琛，梁晨，李京昌

中國石化石油勘探開發(fā)研究院，北京 100083

塔里木盆地順南地區(qū)中—下奧陶統(tǒng)正演沉積數(shù)值模擬

楊偉利，王毅，張眾琛，梁晨，李京昌

中國石化石油勘探開發(fā)研究院，北京 100083

利用正演沉積數(shù)值模擬，恢復(fù)了塔里木盆地順南地區(qū)中—下奧陶統(tǒng)沉積演化過程。綜合分析基礎(chǔ)地質(zhì)條件，以Dionisos軟件為平臺，建立了塔里木順南地區(qū)中—下奧陶統(tǒng)正演沉積數(shù)值模擬的初始模型。通過調(diào)整構(gòu)造沉降、海平面變化、沉積物產(chǎn)率、波浪強度和方向等參數(shù)及模型校驗，以正演方法再現(xiàn)了該地區(qū)三維沉積演化過程，并分析了沉積旋回演化、關(guān)鍵界面特征和內(nèi)部沉積格局。早中奧陶世，塔里木盆地順南地區(qū)發(fā)育四個沉積旋回，四次相對海平面下降形成了中—下奧陶統(tǒng)蓬萊壩組下段、蓬萊壩組上段、鷹山組下段、鷹山組上段等頂面多個局部或區(qū)域短期暴露面，形成了非斷裂控制的短期暴露巖溶。

正演沉積數(shù)值模擬；碳酸鹽巖；中—下奧陶統(tǒng)；順南地區(qū)；塔里木盆地

0 引言

順南地區(qū)位于塔里木盆地塔中北坡，近幾年在中—下奧陶統(tǒng)碳酸鹽巖層系見到良好油氣顯示，儲層是控制油氣富集的最主要因素。目前認(rèn)為油氣發(fā)現(xiàn)的主要儲集空間為受北東向斷裂控制的縫洞型儲層，因此油氣主要沿北東向斷裂帶展布。碳同位素地層學(xué)研究認(rèn)為塔中北坡中—下奧陶統(tǒng)不存在地層缺失[1]，地震剖面上也難以發(fā)現(xiàn)明顯的不整合接觸特征，因此認(rèn)為暴露條件下形成的巖溶儲層欠發(fā)育。這些認(rèn)識把制約油氣富集的儲層空間局限在斷裂帶附近，也制約了勘探空間的拓展。因此，通過沉積正演數(shù)值模擬恢復(fù)順南地區(qū)中—下奧陶統(tǒng)沉積過程，為探索該地區(qū)是否發(fā)育形成巖溶儲層的暴露背景具有重要意義，為下步油氣勘探的拓展提供依據(jù)。

正演沉積數(shù)值模擬選用法國石油研究院的Dionisos軟件實現(xiàn)。該軟件是一個為沉積過程建模的工業(yè)應(yīng)用軟件，在碳酸鹽巖建造和碎屑巖沉積過程中都得到了成功的應(yīng)用。

1 Dionisos沉積模擬基本原理及方法

沉積模擬是以沉積學(xué)理論研究為基礎(chǔ)的一種重要的實驗手段和技術(shù)方法，是對自然界真實過程的一種還原[2]。通過沉積實驗，還原沉積過程和結(jié)果。對沉積過程進行抽象，提取其中的數(shù)學(xué)、物理規(guī)律以及經(jīng)驗公式，然后通過計算機將這種規(guī)律加以實現(xiàn)，直觀、形象地再現(xiàn)各種不同的沉積過程。

數(shù)值模擬根據(jù)模擬過程中的數(shù)據(jù)流向，可分為正演模擬和反演模擬。正演數(shù)值模擬輸入的是控制參數(shù)，而反演數(shù)值模擬輸入的是觀察數(shù)據(jù)。正演模擬建立在假定過程參數(shù)和地層響應(yīng)之間是相互依存的，通過調(diào)試過程參數(shù)的相互作用所產(chǎn)生的地層響應(yīng)來模擬真實地質(zhì)現(xiàn)象。

三維沉積數(shù)值模擬軟件Dionisos是第一個再現(xiàn)沉積演化歷史的模擬工業(yè)應(yīng)用軟件?？梢猿浞掷脦r芯、測井、地震等數(shù)據(jù)，應(yīng)用層序地層學(xué)原理，通過模擬沉積物的沉積過程，定量評估構(gòu)造沉降、海平面變化、沉積物供給、沉積搬運和可容空間之間復(fù)雜的相互作用關(guān)系，為沉積盆地的演化建立一個理想的地質(zhì)歷史模型，如地層、巖相、沉積建造、古水深在時間和空間上的演化。物質(zhì)平衡原理的應(yīng)用使得Dionisos 可以模擬盆地每個時間和空間點的沉積或者剝蝕作用，因此可以進一步模擬盆地的變形和沉積充填，三維定量描述沉積體的形態(tài)，提供沉積盆地中沉積建造、地層巖相、古水深在空間和時間上的演化。

2 工區(qū)概況

模擬工區(qū)包括塔里木塔中北坡的順南地區(qū)和古城墟隆起西部地區(qū)(圖1)。寒武紀(jì)—早奧陶世，塔里木周緣處于拉張背景，盆地內(nèi)部構(gòu)造穩(wěn)定，塔中隆起尚未形成，該地區(qū)為東盆西臺的結(jié)構(gòu)。中晚奧陶世，隨著塔里木周緣區(qū)域動力學(xué)背景的轉(zhuǎn)換及盆地內(nèi)部響應(yīng)，塔中隆起逐漸形成，加里東中晚期斷裂體系開始發(fā)育。順南斜坡位于塔中Ⅰ號帶下盤，現(xiàn)今整體表現(xiàn)為南東高、北西低的斜坡，并在西北部和東南部發(fā)育順托果勒低隆和古城墟隆起， “東西分區(qū)”的構(gòu)造格局明顯。多條北東向、北東東向走滑構(gòu)造帶將塔中北坡與卡塔克隆起分割成東西不同的塊體。

圖1 模擬工區(qū)位置圖Fig.1 The simulated area in Tarim Basin

研究區(qū)發(fā)育以碳酸鹽巖為主的寒武系、中—下奧陶統(tǒng)和上奧陶統(tǒng)良里塔格組，以碎屑巖為主的上奧陶統(tǒng)桑塔木組—新生界。本次研究層位為中、下奧陶統(tǒng)蓬萊壩組、鷹山組、一間房組。塔里木主體區(qū)域在早奧陶世繼承了晚寒武世的沉積面貌[3]，主要發(fā)育開闊臺地、臺地邊緣、臺地斜坡、深水盆地相。蓬萊壩組巖性為淺灰、褐灰色薄—中厚層狀泥—粉晶灰?guī)r、砂屑灰?guī)r、云質(zhì)灰?guī)r與粉—細晶白云巖互層；鷹山組巖性為灰色中厚層狀球粒泥晶灰?guī)r、藻凝塊泥晶灰?guī)r夾砂屑灰?guī)r、云質(zhì)灰?guī)r、灰質(zhì)云巖和白云巖；一間房組巖性為灰、深灰色厚層狀泥晶生物碎屑灰?guī)r及淺灰、灰白色厚層狀礁灰?guī)r。

3 模擬初始模型

初始模型需要確立工區(qū)的范圍、初始基底和模擬時間。本次模擬工區(qū)為長136 km、寬86 km的矩形區(qū)域，模擬時間為中—下奧陶統(tǒng)蓬萊壩組—一間房組沉積時期，即488.3～460.9 Ma，初始模型采用蓬萊壩組沉積初期即488.3 Ma時的古水深。

古水深的確定有多種方法[4-6]，根據(jù)塔里木盆地實際情況，主要采用以沉積相為主、輔助古生物特征確定古水深。一般認(rèn)為克拉通盆地水深在5～10 m[7]，對塔里木早古生代古水深而言，處于斜坡以上的臺地沉積區(qū)域其海水深度小于10 m。金之鈞等[8]曾依據(jù)沉積相帶和古生物資料對塔里木盆地的古水深進行了系統(tǒng)分析，提出塔東地區(qū)早寒武世古水深200～1 000 m， 中寒武世古水深200～2 000 m，晚寒武世古水深200～700 m。早奧陶世繼承了晚寒武世的古地理格局，在斜坡區(qū)以下的區(qū)域，水深應(yīng)該不小于700 m，但具體多少尚無方法來確定。

綜合沉積相、沉積巖性、古生物特征及前人研究成果，半定量地確定出早古生代塔里木盆地的古水深，臺地相區(qū)水深10 m，盆地相區(qū)700 m，斜坡帶處于過度水深。在區(qū)域相帶轉(zhuǎn)換處的蓬萊壩組沉積初期水深圖上，截取研究區(qū)古水深，由散點的形式導(dǎo)入Dionisos軟件，并經(jīng)過變差函數(shù)運算和平滑運算，最終得到初始基底的模型(圖2)。

圖2 順南地區(qū)正演數(shù)值模擬初始模型Fig.2 The forward numerical sedimentary simulation model of Shunnan area

4 模擬參數(shù)

碳酸鹽巖沉積數(shù)值正演模擬中需要三類參數(shù)：可容空間(構(gòu)造沉降、海平面變化和沉積物供給)、物質(zhì)(沉積物類型、產(chǎn)率)、能量(波浪方向、波浪能量)。

4.1 構(gòu)造沉降

構(gòu)造沉降采用艾里均衡公式計算的結(jié)果*云金表，楊偉利，李京昌，等.塔里木下古生界碳酸鹽巖層系古構(gòu)造演化研究，2013.?；驹硎歉鶕?jù)質(zhì)量守恒法則，假設(shè)地層的骨架厚度不隨埋深的增加而發(fā)生改變。構(gòu)造作用引起的沉降(STT)和負(fù)荷均衡作用引起的沉降(SL)組成了盆地在某一時刻的基底總沉降量(ST)。構(gòu)造作用引起的沉降運動是被構(gòu)造沉降所反映的，表示為：

STT(t)=ST(t)-SL(t)

(1)

式中STT(t)為時刻t的構(gòu)造沉降量；ST(t)為時刻t的基底總沉降量，可以用時刻t盆地的古地層厚度Hs(t)+古水深Hw(t)表示；SL(t)是指時刻t由盆地地層負(fù)荷引起的沉降量，可由艾里均衡公式求得。即：

SL(t)=Hs(t)(ρs-ρw)/ (ρm-ρw)

(2)

式中古地層厚度Hs(t)是用回剝方法確定的，在這里可用現(xiàn)今測得的沉積厚度進行計算，因為白云巖和灰?guī)r是寒武系的主要巖性，可以認(rèn)為其不存在沉積壓實作用。ρs是盆地沉積層的平均密度，ρm和ρw是地幔和水在0℃時的密度。由于白云巖為寒武系的主要巖性，故ρs取值2 870 kg/m3，ρm取值3 330 kg/m3，ρw取值1 030 kg/m3。

由上面的(1)式和(2)式可得出構(gòu)造沉降量的計算公式：

STT(t)=Hs(t)+Hw-Hs(t)(ρs-ρw)/ (ρm-ρw)

(3)

由上式得到的即為總沉積量去掉沉積負(fù)載作用影響之后的構(gòu)造沉降量(圖3)。從構(gòu)造沉降看，工區(qū)南部沉降幅度小，北部幅度大，總體上分布在900～1 200 m。

圖3 順南模擬工區(qū)早中奧陶紀(jì)構(gòu)造沉降量Fig.3 Early-Middle Ordovician tectonic subsidence quantity in Shunnan area

4.2 海平面升降

海平面的變化對于沉積相的遷移、碳酸鹽緩坡—臺地的形成和演化、沉積類型都有影響。海平面分為絕對海平面和相對海平面兩種，Dionisos軟件中使用的是絕對海平面。目前有多種方法確定海平面變化[9-10]，包括上超曲線法、測井曲線法、層序體系域法、海岸上超法、沉積相序法、群落生態(tài)法、穩(wěn)定同位素法、數(shù)學(xué)模型法、Fishier 圖解法等，這些方法反應(yīng)的是海平面的變化趨勢或者相對海平面變化，目前很難得到絕對海平面變化。通過對全球海平面變化曲線、塔里木地區(qū)不同方法得到的海平面變化曲線分析，本次選用高志前等[11]根據(jù)塔里木盆地巴楚、塔中地區(qū)4口井和柯坪露頭剖面的微量元素數(shù)據(jù)、鍶同位素、碳氧同位素數(shù)據(jù)，利用塔中1井的碳氧同位素資料編制的海平面升降變化曲線，雖然也代表了相對海平面變化，但是就碳同位素背景演化曲線而言，它具有全球的一致性，可能反映了全球大范圍的海平面變化規(guī)律[12]。

4.3 沉積物產(chǎn)率

碳酸鹽巖產(chǎn)率主控因素包括確定各種巖性產(chǎn)率與水深、時間的定量關(guān)系。通過對鉆井和地震反射結(jié)構(gòu)分析，確定了研究區(qū)三種主要沉積物類型：灰?guī)r、泥巖和礁灘灰?guī)r。礁灘灰?guī)r代表生物礁障壁體系(珊瑚、藻類)及高能沉積物的所有沉積，灰?guī)r代表正常水動力條件下的碳酸鹽巖沉積，泥代表開放海洋環(huán)境下的沉積。在Dionisos中，沉積物產(chǎn)率與水深的關(guān)系用0～1的范圍表示，1代表該類型沉積物達到最大產(chǎn)率，0代表沉積物停止沉積或碳酸鹽巖停止生長。我國現(xiàn)代珊瑚礁的生長帶水深范圍在3～18 m[13]，Doddetal.[14]認(rèn)為，現(xiàn)代珊瑚礁在15 m以淺水深范圍內(nèi)生長最豐富，吳亞生等人[6]認(rèn)為珊瑚礁亮晶骨架相水深范圍為2～10 m。大部分灰?guī)r形成在水體清淺暖環(huán)境下，水深以下10～20 m是光照強度飽和區(qū)[15]。考慮到塔里木盆地順南地區(qū)礁灘發(fā)育的高能環(huán)境，定義礁灘、灰?guī)r沉積物最大產(chǎn)率的水深是0～20 m，水深大于該值則逐漸降低，礁灘沉積產(chǎn)率隨著水深的增加呈指數(shù)遞減，灰?guī)r產(chǎn)率則隨著水深的增加逐漸遞減。由于沒有外部物源，泥的沉積主要是在斜坡及盆地區(qū)，在一定深度內(nèi)隨水深增加產(chǎn)率增加(圖4)。產(chǎn)率和時間的關(guān)系可以根據(jù)鉆井地層厚度計算，該地區(qū)選用古隆1井計算灰?guī)r產(chǎn)率，結(jié)果如圖5所示，礁灘產(chǎn)率很高，泥產(chǎn)率很低。沉積物產(chǎn)率，特別是水深和產(chǎn)率的關(guān)系是一個可控參數(shù)，根據(jù)實際情況進行反復(fù)嘗試與調(diào)整，最終采用圖4、5中所示結(jié)果。

圖4 沉積物產(chǎn)率與水深的關(guān)系Fig.4 The relation between deposition rate and water depth

圖5 沉積物產(chǎn)率與時間的關(guān)系Fig.5 The relation between deposition rate and geological time

4.4 能量參數(shù)

能量參數(shù)需要定義波浪大小和方向。軟件中可以定義兩套波浪：正常海浪和風(fēng)暴浪，用基準(zhǔn)深度、波浪傳播角度和頻率表達?；鶞?zhǔn)深度是波浪影響開始減少的參考深度，它等于浪基面除以2.5。浪基面深度認(rèn)識不一，一般認(rèn)為在20多米，此處選取25 m。風(fēng)暴浪浪基面深度認(rèn)識差異也比較大，一般風(fēng)暴浪波及深度在40～50 m[16-17]。因此，在能量參數(shù)的確定中，正常海浪基準(zhǔn)深度設(shè)定在5～15 m范圍，風(fēng)暴浪基準(zhǔn)深度范圍設(shè)定為10～25 m。波浪傳播的角度以斜交和垂直岸線為主，變化范圍15°～150°。波浪的頻率指波浪類型在一年中出現(xiàn)的頻率，全年中都是好天氣則波浪定義為100%，風(fēng)暴天氣僅僅在一年中很小一部分時間發(fā)生。設(shè)置正常波浪作用頻率是100%，風(fēng)暴浪的作用頻率為25%。

在構(gòu)造沉降、海平面變化、沉積物產(chǎn)率和能量參數(shù)中，構(gòu)造沉降和海平面變化是不變的量，沉積物產(chǎn)率和能量參數(shù)是可控的。沉積物產(chǎn)率調(diào)整的是隨水深的關(guān)系，最大產(chǎn)率的水深范圍參數(shù)在0～10 m、0～20 m兩組中調(diào)整。能量參數(shù)中，基準(zhǔn)深度中正常和風(fēng)暴浪基準(zhǔn)深度在5 m和10 m、10 m和15 m、15 m和20 m三組中調(diào)整；波浪傳播角度則在15°、45°、75°、105°、135°中調(diào)整。通過對各種條件下沉積演化過程模型的分析，最終確定最大產(chǎn)率為0～10 m、基準(zhǔn)深度10 m和15 m、波浪傳播角度75°條件下的模型和實際最接近。

5 沉積模擬結(jié)果正確性分析

通過對模擬結(jié)果的模擬單井沉積旋回、連井剖面及臺地結(jié)構(gòu)等方面與實際資料進行對比、局部地區(qū)和三維地震資料恢復(fù)的古地貌做了對比驗證了模擬結(jié)果的正確性。

5.1 單井對比

以沉積正演數(shù)值模擬的結(jié)果生成模擬井剖面，通過與實際井剖面對比驗證模擬結(jié)果正確性。根據(jù)該地區(qū)鉆井完鉆實際層位大部分在鷹山組，對比分析了順南1、順南5、順南2、古隆2和古隆1等單井及連井鷹山組—一間房組沉積旋回。以位于臺地邊緣的古隆1井鷹山組為例(圖6)，古隆1井鷹山組實鉆厚度450 m，模擬井沉積厚度430 m；實際鉆井中高能沉積環(huán)境中沉積的亮晶顆?；?guī)r與模擬井中高能沉積物進行對比基本吻合；實鉆井小波變換得到的準(zhǔn)層序與模擬井巖性變化劃分的準(zhǔn)層序?qū)Ρ蕊@示，大部分準(zhǔn)層序是一致的，只是中—下部略微有差異，實際井小波變換得到的準(zhǔn)層序頻率比模擬井略多，符合率在80%以上，說明正演模擬結(jié)果與實際相符。

圖6 模擬井與實鉆井剖面對比(古隆1井)Fig.6 The comparison of the modelling section and the drilling section in Well GL1

5.2 地層結(jié)構(gòu)對比

圖7 模擬地層結(jié)構(gòu)與地震反射結(jié)構(gòu)對比Fig.7 The geometry comparison of the seismic reflection and the simulation model

6 討論

數(shù)值沉積模擬結(jié)果顯示，該地區(qū)從奧陶紀(jì)蓬萊壩組沉積開始，經(jīng)歷了四期沉積旋回，在臺地區(qū)形成了多期暴露或局部暴露，即形成了暴露不整合(圖8)。蓬萊壩組沉積初期(488.3 Ma)，臺地區(qū)全部被水體覆蓋；此后，海平面上升，水體加深，到487.05 Ma，工區(qū)內(nèi)部分化，水體分布不均，形成臺地內(nèi)部洼地，代表一次短暫海侵；隨后海平面下降，到485.55 Ma，洼地周緣部分露出水面，接受短時間的暴露侵蝕，代表了一次水退結(jié)果。到482.05 Ma，相對海平面上升，研究區(qū)被淹沒，早期暴露地區(qū)形成淹沒不整合；到蓬萊壩組沉積結(jié)束(479.3 Ma)，該地區(qū)水體完全退出到斜坡帶以下，臺地區(qū)完全暴露，是一次比較大的海退。鷹山組沉積時期，隨著可容空間變化，該地區(qū)也經(jīng)歷了幾次相對海平面變化。在這種短期暴露之下，局部地區(qū)形成了中—下奧陶統(tǒng)內(nèi)部小規(guī)模的平行不整合、微角度不整合、暴露侵蝕面等，在該界面之下發(fā)育一系列具有暴露或侵蝕標(biāo)識的溶洞特征或在地震上界面之下串珠集中分布。

從模擬結(jié)果看，順南地區(qū)中—下奧陶統(tǒng)內(nèi)部并不是目前所認(rèn)識的不存在地層缺失，而是存在多個內(nèi)部暴露界面。碳同位素地層學(xué)揭示的地層間斷是一個大尺度間斷界面，相對小尺度間斷則無法識別。模擬工區(qū)顯示即為一個短時間內(nèi)的沉積間斷，沉積間斷時間在1 Ma左右，剝蝕作用比較弱，在地震剖面上也沒有形成明顯的削截、上超等特征的反射結(jié)構(gòu)。區(qū)域動力學(xué)環(huán)境分析認(rèn)為該時期是一個區(qū)域應(yīng)力場由拉張轉(zhuǎn)換為擠壓的過程，相對海平面下降應(yīng)該是應(yīng)力轉(zhuǎn)換在盆地內(nèi)部的響應(yīng)和區(qū)域海平面變化綜合作用的結(jié)果。這種綜合結(jié)果導(dǎo)致了中—下奧陶統(tǒng)內(nèi)部發(fā)育多次短時間暴露，形成了短期暴露巖溶背景。

References)

[1] 張智禮，李慧莉，譚廣輝，等. 塔里木中央隆起區(qū)奧陶紀(jì)碳同位素特征及其地層意義[J]. 地層學(xué)雜志，2014，38(2)：181-189. [Zhang Zhili, Li Huili, Tan Guanghui, et al. Carbon isotope chemostratigraphy of the Ordovician system in Central Uplift of Tarim Basin[J]. Journal of Stratigraphy, 2014, 38(2): 181-189.]

[2] 李翔. 沉積盆地充填過程的計算機模擬研究[J]. 現(xiàn)代地質(zhì)，1994，8(4)：393-402. [Li Xiang. Computer simulation of sedimentary basin fill[J]. Geoscience, 1994, 8(4): 393-402.]

[3] 劉迪，張哨楠，謝世文，等. 塔里木盆地塔中地區(qū)奧陶系鷹山組層序地層特征[J]. 礦物巖石地球化學(xué)通報，2014，33(1)：55-64. [Liu Di, Zhang Shaonan, Xie Shiwen, et al. Characteristics of sequence stratigraphy of the Ordovician Yingshan Formation in the central part of Tarim Basin[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2014, 33(1): 55-64.]

[4] 萬錦峰，鮮本忠，佘源琦，等. 基于伽馬能譜測井信息的古水深恢復(fù)方法——以塔河油田4區(qū)巴楚組為例[J]. 石油天然氣學(xué)報，2011，33(6)：98-103. [Wan Jinfeng, Xian Benzhong, She Yuanqi, et al. Palaeobathymetric reconstruction based on natural gamma ray spectrometry logging data-by taking Bachu Formation in region 4 of Tahe oilfield for example[J]. Journal of Oil and Gas Technology, 2011, 33(6): 98-103.]

[5] 王敏芳，焦養(yǎng)泉，任建業(yè)，等. 沉積盆地中古地貌恢復(fù)的方法與思路——以準(zhǔn)噶爾盆地西山窯組沉積期為例[J]. 新疆地質(zhì)，2006，24(3)：326-330. [Wang Minfang, Jiao Yangquan, Ren Jianye, et al. Method and thinking of Palaeogeomorphologic reconstruction in sedimentary basin-Example from depositional stage of Xishanyao Formation in Junggar Basin[J]. Xinjiang Geology, 2006, 24(3): 326-330.]

[6] 吳亞生，范嘉松. 根據(jù)生物礁定量計算茅口期全球海平面變化幅度[J]. 中國科學(xué)(D輯)：地球科學(xué)，2001，31(3)：233-242. [Wu Yasheng, Fan Jiasong. Calculating eustatic amplitude of Middle Permian from reefs[J]. Science China (Seri. D): Earth Sciences, 2001, 31(3): 233-242.]

[7] 賈振遠，李之琪. 碳酸鹽巖沉積相和沉積環(huán)境[M]. 武漢：中國地質(zhì)大學(xué)出版社，1989：1-186. [Jia Zhenyuan, Li Zhiqi. Carbonate Sedimentary Facies and Sedimentary Environment[M]. Wuhan: China University of Geosciences Press, 1989: 1-186.]

[8] 金之鈞，呂修祥，王毅，等. 塔里木盆地波動過程及其控油規(guī)律[M]. 北京：石油工業(yè)出版社，2003：1-134. [Jin Zhijun, LüXiuxiang, Wang Yi, et al. Wave Process ant Its Oil-Controlling Law of Tarim Basin[M]. Beijing: Petroleum Industry Press, 2003: 1-134.]

[9] 李祥輝，張潔. 海平面及海平面變化綜述[J]. 巖相古地理，1999，19(4)：61-72，41. [Li Xianghui, Zhang Jie. Problems and suggestions on sea-level and its changes[J]. Sedimentary Facies and Paleogeography, 1999, 19(4): 61-72, 41.]

[10] 趙玉光，劉寶珺，許校松. 地質(zhì)歷史時期海平面變化曲線的數(shù)學(xué)模型[J]. 地球科學(xué)，2000，25(2)：117-121. [Zhao Yuguang, Liu Baojun, Xu Xiaosong. Mathematical model of sea-level change curves in geological history[J]. Earth Science, 2000, 25(2): 117-121.]

[11] 高志前，樊太亮，李巖，等. 塔里木盆地寒武奧陶紀(jì)海平面升降變化規(guī)律研究[J]. 吉林大學(xué)學(xué)報：地球科學(xué)版，2006，36(4)：549-556. [Gao Zhiqian, Fan Tailiang, Li Yan, et al. Study on eustatic sea-level change rule in Cambrian-Ordovician in Tarim Basin[J]. Journal of Jilin University: Earth Science Edition, 2006, 36(4): 549-556.]

[12] 李玉成. 華南晚二疊世碳酸鹽巖碳同位素旋回對海平面變化的響應(yīng)[J]. 沉積學(xué)報，1998，16(3)：52-53，55-57. [Li Yucheng. The carbon isotope cyclostratigraphic responses to sea level change in Upper Permian limestone from South China[J]. Acta Sedimentologica Sinica, 1998, 16(3): 52-53, 55-57.]

[13] 呂炳全，王國忠，全松青. 海南島珊瑚岸礁的特征[J]. 地理研究，1984，3(3)：1-16. [Lü Bingquan, Wang Guozhong, Quan Songqing. The characteristics of fringing reefs of Hainan Island[J]. Geographical Research, 1984, 3(3): 1-16.]

[14] Dodd J R, Stanton R J. Paleoecology, Concepts and Applications[M]. New York: John Wiley and Sons, 1981: 1-559.

[15] 鄭艷剛，袁美云，徐維勝. 臺地邊緣生物礁生長規(guī)律和控制因素綜合研究[J]. 內(nèi)蒙古石油化工，2009，34(6)：111-112. [Zheng Yangang, Yuan Meiyun, Xu Weisheng. Study on growth rule and control factors of platform edge reef[J]. Inner Mongolia Petrochemical Industry, 2009, 34(6): 111-112.]

[16] 金振奎，石良，高白水，等. 碳酸鹽巖沉積相及相模式[J]. 沉積學(xué)報，2013，31(6)：965-979. [Jin Zhenkui, Shi Liang, Gao Baishui, et al. Carbonate facies and facies models[J]. Acta Sedimentologica Sinica, 2013, 31(6): 965-979.]

[17] Brett C E, Boucot A J, Jones B. Absolute depths of Silurian benthic assemblages[J]. Lethaia, 1993, 26(1): 25-40.

A Three-dimensional Numerical Sedimentary Simulation Model for the Lower-Middle Ordovician of Shunnan Area, Tarim Basin

YANG WeiLi，WANG Yi，ZHANG ZhongChen，LIANG Chen，LI JingChang

Petroleum Exploration and Production Research Institute, SINOPEC, Beijing 100083, China

The numerical sedimentary forward simulation of Lower-Middle Ordovician of Shunnan area in Tarim Basin was presented by means of software Dionisos, which are controlled by primary model, tectonic subsidence, fluctuation of sea level, sedimentation rate and wave energy. Through the validation, the appropriate model has been chosen and then the sedimentary process is clear. The evolution of this area has four deposititional cycles. The distribution of carbonate reef, beach and short-time exposed surfaces has been analyzed. During Early-Middle Ordovician, there developed four frequent short-time exposure surfaces which occurred at the top of the depositional cycle in Shunnan area and resulted in the formation of vugular-pore in carbonate deposits subjected to leaching and dissolution of meteoric fresh water,leading to which makes a change that only the fracture dominates the development of the carbonate reservoir.

sedimentary simulation; the lower-middle Ordovician; carbonate; Shunnan area; Tarim Basin

1000-0550(2017)01-0117-07

10.14027/j.cnki.cjxb.2017.01.012

2015-12-07；收修改稿日期： 2016-05-01

國家科技重大專項(2016ZX05005-002)[Foundation: National Science and Technology Major Project, No. 2016ZX05005-002]

楊偉利，男，1973年生，高級工程師，沉積、石油地質(zhì)，E-mail:yangwl.syky@sinopec.com

P618.13

A