模型正演技術在巴中地區須四段主河道砂體識別中的應用

肖 偉， 王明飛， 何志勇

(中石化勘探南方分公司，成都 610041)

?

模型正演技術在巴中地區須四段主河道砂體識別中的應用

肖 偉， 王明飛， 何志勇

(中石化勘探南方分公司，成都 610041)

巴中地區須四段主要為陸相辮狀河三角洲前緣沉積，前人的研究表明，主河道沉積微相是該區須四段高產富集的主控因素之一，但是由于須四段分流河道側向遷移頻繁，縱向上相互疊置，導致地震資料橫向變化較大，利用這些資料研究主河道砂體在剖面上展布其可靠性較差。運用模型正演技術，借助該區域三維剖面反射特征，利用已鉆井資料提供巖石物理參數，設計地質模型，進而研究和分析模型正演結果，明確了巴中地區須四段主河道砂體的地震響應特征，建立了主河道地震識別模式，在此基礎上追蹤出主河道砂體的分布范圍，取得了較好的應用效果。

模型正演； 主河道； 砂體識別； 巴中地區

0 引言

模型正演技術是解決地震多解性的重要手段,也是勘探復雜目標的有效技術。在辮狀河三角洲前緣巖性油氣藏勘探開發中，主河道砂體的識別是關鍵。由于河道砂體尺度小、散度大，橫向非均質性強，再加上地震資料分辨率有限，砂體識別存在不確定性。前人在研究時，一般都是根據地質特征設計不同的地質模型，進行正演模擬，得出河道砂體的地震識別模式，對提高地震資料巖性解釋的可靠性及提高油氣藏勘探成功率方面具有重要的參考價值[1-4]。借助前人的研究方法，針對優質河流相砂體地震識別難度大、預測存在多解性的特點，這里運用模型正演技術，首先在已知井點建立地質模型與地震剖面對比分析，在此基礎上設計出不同主頻、主河道砂體厚度、地球物理參數、巖性組合的多種地質模型，正演模擬各種模型的地震響應，綜合得出巴中地區須四段主河道砂體的地震識別模式，并在平面上追蹤出主河道砂體的分布。

1 研究區概況

巴中地區區域構造上位于通南巴構造帶的西南，九龍山構造東南側，西側緊鄰梓潼凹陷，南部為川中古隆起北斜坡，東南部為儀隴-平昌平緩構造帶；區域構造主要屬于通南巴構造隆起帶西南端和川北坳陷西部。上三疊統須四段砂巖儲層是巴中地區重要目的層，已有兩口井在須四段試獲工業氣流，其中YL17井須四段測試獲得22.635 9×104m3/d 工業氣流。

前期研究成果表明，巴中地區須四段為多期辮狀河三角洲前緣沉積，水下分流河道相互疊置，縱向上可以分為兩套砂體，①上部以細-中粒長石巖屑砂巖沉積為主；②下部以細粒、中粒巖屑砂巖沉積為主。其中上部長石巖屑砂巖為須四段主要儲層發育段，主河道沉積微相和裂縫是儲層富集高產的主控因素，作者研究范圍在須四段上部長石巖屑砂巖巖性圈閉內。

2 正演模擬

地震正演模擬就是利用已有資料建立地下地質模型，根據地震波在地下介質中的傳播原理，通過射線追蹤或波動方程偏移等方法，正演模擬計算出對應于建立地質模型的地震記錄。其主要目的是通過模擬記錄與實際地震記錄的對比分析，校正初始模型，使之更接近地下真實的地質情況。模型正演包括模型建立、正演模擬、記錄分析等過程。通過建立已知的地震地質模型進行地震正演，能夠幫助人們直觀地認識地震波在地層中的傳播規律，識別地質構造及油氣藏的地震響應，從而指導地震數據的采集、處理和解釋，模型正演研究是地震勘探方法的基礎。這一過程涉及到巖石物理、地質構造以及數學方法等方面[5-7]。對于須四段主河道的識別，由于研究區鉆井資料缺乏，因此需要根據實際地震剖面，不斷調整模型設計方案，對主河道厚度、物性及巖性變化進行正演模擬，進而建立主河道砂體的地震識別模式。

2.1 地質模型的建立及巖石物理參數分析

依據過YL17井實際的地震剖面(圖1(a))及研究區地質情況，建立地質模型(圖1(b))，該井段在須四段上部發育主河道砂體，從地震剖面上可以明顯看出，須四段上砂體底界波峰下拉，振幅減弱，頂界出現復波，根據地震反射外部形態來看，是較為典型的河道充填相特征。

圖1 巴中地區地震剖面及主河道地質模型

巖石物理參數是模型正演的重要參數之一，研究區內共有兩口已鉆至須四段底界的鉆井，通過單井沉積微相分析，只有YL17井在須四段鉆遇主河道砂體，因此利用YL17測井資料來獲取須四段基礎的巖石物理參數。首先根據YL17井須四段的錄井巖性(圖2)，將須四段自上而下分為鈣屑砂巖、泥巖、主河道砂巖、長石巖屑砂巖、巖屑砂巖五種巖性，再加上須四段頂、底兩套泥巖段共七個不同巖性段，利用聲波、密度測井資料，分別統計出各個巖性段的速度和密度參數(圖1(b))。

圖2 YL17井綜合柱狀圖

2.2 模型正演模擬分析

地震子波是正演模擬的另一個重要參數，為了便于研究，這里選用Ricker子波作為正演模擬子波，而Richer子波主頻決定了地震垂向分辨率[8]，直接影響模型正演的結果。因為巴中地區須四段地層橫向上巖性及厚度變化較快，不同的子波主頻對正演結果影響較大，因此首先討論研究區須四段Ricker子波主頻的選取。

選擇不同主頻的Ricker子波進行正演模擬，從圖3可以看出，隨著子波主頻增加，地震分辨率逐漸提高，在主頻大于30 Hz時可以很好地將主河道砂體底界與須四段上砂底界分開，但由于主河道砂體頂部夾有兩套10 m左右的薄層，即使子波主頻達到45 Hz也難以區分開主河道砂體頂界。

圖3 不同主頻的Ricker子波正演模擬結果

為了更好地將正演模擬結果與實際地震剖面進行對比，作者分析了巴中地區須四段頻譜(圖4)，從頻譜圖上可以看出，巴中地區須四段子波主頻在30 Hz左右，因此正演模擬子波都選用主頻為30 Hz的Ricker子波。

圖4 巴中地區須四段地震資料頻譜圖

由于在主河道遷移過程中，隨著水動力、物源供給量等因素的影響，主河道砂體在砂體厚度、物性及巖性組合在空間上會發生變化，為了研究主河道砂體在各種條件下的地震響應特征，在過YL17井主河道砂體模型的基礎上，分別設計了主河道砂體厚度變化、河道砂體物性變化以及須四段頂部巖性組合變化的模型進行正演模擬。

首先研究主河道砂體厚度變化對正演結果的影響，分別設計了主河道砂體厚度由20 m變化到70 m的地質模型，其他巖性段的厚度保持不變，進行正演模擬(圖5)。從正演模擬結果來看，隨著主河道砂體厚度增大：①主河道頂界波峰變化不大；②主河道底部波峰振幅逐漸減弱，當主河道砂體厚度達到一定程度(即主河道砂體底界與須四段上砂底界之間厚度小于四分之一波長，由井上統計可知巴中地區須四段平均速度約為 5 100 m/s，主頻約 30 Hz，因此四分之一波長約為43 m)，主河道砂體底界與須四段上砂底界相互干涉形成復波，隨著主河道砂體厚度繼續增大，復波振幅又逐漸增強；③須四段上砂底界波峰下拉，振幅逐漸減弱，當主河道砂體厚度增加到一定程度時，變化趨勢與主河道砂體底界相同；④須四上砂內部同相軸增多。

圖5 主河道砂體厚度變化模型及正演模擬結果

其次研究主河道砂體物性的變化對正演結果的影響。由于巴中地區須四段為致密砂巖儲層，儲層物性變化范圍不大，因此作者分別設計了主河道砂體速度為4 700 m/s、4 600 m/s兩個地質模型，物性比原始模型的物性稍好，然后進行正演模擬，從正演結果來看(圖6)，隨著河道砂體物性逐漸變好：①主河道砂體頂界波峰變化不大；②主河道砂體底界波峰振幅逐漸增強；③須四段上砂底界波峰下拉，振幅減弱。

圖6 主河道砂體物性變化模型及正演模擬結果

最后研究須四段上砂體內部巖性組合變化對正演結果的影響，由于巴中地區須四段整體為辮狀河三角洲前緣沉積，從測井、錄井資料上分析，巴中地區須四段為大套的砂巖夾薄層泥巖的“砂包泥”結構，因此在原始地質模型的基礎上，分別設計了無主河道砂體模型、無頂部泥巖夾層模型和無主河道且內部夾泥巖薄層模型進行正演模擬(圖7)，從正演模擬結果來看：①不發育主河道砂體時，須四段上砂底界波峰無下拉現象，內部為寬波谷反射；②主河道砂體頂部無泥巖夾層時，主河道砂巖頂界波峰上提現象與原始模型相比不太明顯，主河道砂體底界與須四段上砂底界反射變化不大；③主河道砂體不發育且須四上砂體內部夾薄層泥巖(泥巖夾層厚度為2 m～5 m，速度為4 000 m/s)時，須四段上砂底界波峰無下拉現象，振幅在夾泥巖段有所增強，須四上砂內部為寬波谷反射。

圖7 須四段砂體巖性組合變化模型及正演模擬結果

綜合以上正演模擬結果可以發現：主河道砂體發育時，主河道砂體頂界反射變化不大，須四段上砂底界波峰下拉，振幅減弱，須四段上砂內部同相軸增多，主河道砂體厚度增大時，須四段上砂底界出現復波，且復波振幅隨著主河道砂體的厚度增大而增強，主河道砂體物性變好時，須四段上砂內部同相軸振幅增強。因此在正演模擬結果的基礎上，可以建立巴中地區須四段主河道砂體的地震識別模式：①須四段上砂底界波峰下拉、波峰振幅減弱或出現復波；②須四段上砂體內部同相軸增多，出現不連續的短波峰反射。

3 主河道砂體的地震識別

3.1 須四段層序地層特征

區域構造研究表明，須三段與須四段沉積不整合界面為“安縣運動”的產物，通常將此界面作為劃分須家河早期和須家河晚期兩個成盆期的重要依據[9]。在巴中通江鐵溪剖面，須四段發育厚層礫巖沉積，可見底部對須三段炭質泥巖沖刷面，可作為層序界面識別標準之一(圖8)。

圖8 須四段底部礫巖沖刷接觸面

圖9 YL3-YL17井須家河組SQ4連井層序格架

巴中地區須四段為辮狀河三角洲前緣沉積，三角洲前緣水下分流河道側向遷移頻繁，結合露頭資料，YL3-YL17井層序-沉積特征(圖9)，以及地震層序研究表明，主河道沉積特征清晰，在地震剖面上，地震反射特征明顯，須三段地層被須四段削截，綜合反映了須四底界面為具有區域沖刷侵蝕性質的不整合面，為須四段底界面；須四段上覆須五段地層，須五段以泥巖沉積為主，地震上為連續中低頻強反射特征，區域上可連續追蹤，作為須四段頂界面(圖10)，總的來說，須四段層序具有頂部整合，底部沖刷侵蝕不整合的層序地層特點。

圖10 YL3-YL17井須四段地震層序特征

圖11 巴中地區須四段斷裂分布平面圖

3.2 層序地層格架內主河道砂體精細刻畫

由YL17井須四段速度、密度統計可以看出，主河道砂體縱波波阻抗平均為12 411 m/s*g/cm3，分支河道砂體縱波阻抗平均為13 535 m/s*g/cm3，從單井波阻抗上分析，利用波阻抗反演能夠較好地將主河道砂體和分支河道砂體分開，但是川東北地區須四段整體表現為低孔低滲的致密砂巖儲層。另外巴中地區須四段平均埋深大于4 000 m，因此主河道砂體與分支河道砂體在物性上不會相差太大，主河道砂體波阻抗與分支河道砂體波阻抗會有較大部分重疊，利用波阻抗反演來識別主河道砂體較困難，且結果存在多解性；由主河道砂體模型正演結果可以看出，主河道砂體地震反射外部形態、內部形態、振幅、頻率等地震相特征與分支河道地震相特征有較大的差別，單就模型正演結果來看，利用波形分類、RMS振幅等地震屬性應該能很好地刻畫主河道砂體平面展布特征，但是由于巴中地區須四段斷裂特別發育(圖11)，地層破壞較嚴重，再加上主河道砂體橫向上變化較快，利用地震屬性識別主河道砂體也存在較大的不確定性。

須四段地震主頻不高，地震上識別單個薄砂體較困難，當多個單砂體縱向上相互疊置，形成了厚度較大、物性較好的主河道砂體，地震上就會產生響應，使利用地震剖面來識別主河道砂體成為可能。依據巴中地區須四段主河道砂體地震識別模式，在層序地層格架內采用“相面法”在須四段長石巖屑砂巖圈閉內精細追蹤出主河道砂體的展布(圖12、13)。由主河道砂體平面分布圖可以看出，巴中地區須四段主河道砂體呈近東西向展布，與YL17井古水流方向一致(圖14)，南北兩套主河道砂體被中間富泥巖夾層的分支河道隔開，東部靠近物源，河道寬度大，縱向切割較深，砂體厚度較大，向西主河道逐漸分流，河道變窄，切割較淺，砂體厚度較小。

圖12 巴中地區須四段主河道分布平面

圖13 巴中地區須四段主河道砂體典型地震剖面

4 結論

模型正演技術是解決未知問題的有效手段，針對研究區主河道砂體的各種變化情況，設計對應的地質模型，采用合適的正演模擬方法，選擇了合理的子波主頻，開展正演模擬研究，建立須四段主河道砂體地震識別模式，在此基礎上，在地震剖面上精細追蹤出主河道砂體的分布，為后續儲層預測工作提供了指導

圖14 YL17井須四段古水流分析圖

[1] 許輝群，桂志先，佘曉宇，等.模型正演技術在SZ地區砂體識別中的應用[J].工程地球物理學報，2010，7(4)：399-402. XU H Q,G ZH X,SHE X Y, et al. Application of forward modeling to sand identification in SZ area.Journal of Engineering Geophysics，2010，7(4):399-402.(In Chinese)

[2] 劉軍迎，雍學善，高建虎,等.模型正演技術在碳酸鹽巖油氣藏地震資料解釋中的應用[J].巖性油氣藏，2007，19(1)：109-112． LIU J Y,YONG X SH,GAO J H, et al. Application of farward modeling to seismic data interpretation in carbonate Reservoir.Lithologic Hydrocarbon Reservoir，2007，19(1):109-112.(In Chinese)

[3] 彭真明，陳艷玲.砂、泥巖互層透鏡體模型正演及分析[J].物化探計算機技術，2001，23(2)：101-104. PENG ZH M,CHEN Y L. Farward modeling and analysis to lens in interbeded sandstones and mudstones.Computing echniques for Geophysical and Geochemical Exploration，2001，23(2):101-104.(In Chinese)

[4] 周義軍，熊玉萍.模型正演在鄂托克前旗地區的應用[J].石油物探，2008，47(2)：161-166. ZHOU Y J,XIONG Y P. Application of forward modeling in Etuoke area.Geophysical Prospecting for Petroleum，2008，47(2):161-166.(In Chinese)

[5] 高秋菊，姜效典.地震模型正演技術在復雜油氣藏勘探中的應用[J].內蒙古石油化工，2008，16:118-119. GAO Q J,JIANG X D. Application of farward modeling in exploration of complex reservoir.Inner Mongolia Petrochemical industry，2008，16:118-119.(In Chinese)

[6] 許輝群，桂志先，孫贊東，等.基于模型正演的研究區砂體識別方法研究[J].石油天然氣學報，2011，33(1):118-119. XU H Q,GUI ZH X,SUN Z D, et al. The research of sandstones identification based on forward modeling.Journal of Oil and Gas Technology，2011，33(1):118-119.(In Chinese)

[7] 鞏建強. 正演模擬法在油氣勘探中的應用.油氣地球物理，2009，7(2):13-16. GONG J Q. The application of the forward modeling in exploration.Petroleum Geophysics，2009，7(2):13-16.(In Chinese)

[8] 李慶忠.走向精確勘探的道路—高分辨率地震勘探系統工程剖析[M].北京：石油工業出版社，1994. LI Q ZH.Engineering analysis of high resolution seismic exploration system to the way of accurate exploration.Beijing：Petroleum industry publishing house，1994.(In Chinese)

[9] 王金琪.龍門山印支運動主幕辨析—再論安縣構造運動[J].四川地質學報，2003，23(2)：65-69． WANG J Q. Longmen mountain indosinian orogeny of the Anxian movement.Acta geologica Sichuan，2003，23(2):65-69.(In Chinese)

Application of sand identification from the main channal of the 4thmember of Xujiahe formation in BZ area using forward modeling

XIAO Wei, WANG Ming-fei, HE Zhi-yong

(South Exploration Company,Sinopec,Chengdu 610041,China)

Forward modeling can be used to indicate the reflection characteristics of geological body, especially in the area where have little wells. The 4th member of Xujiahe formation of BZ area belongs to terrestrial braided river delta front deposits, including the channel migration frequent lateral and vertical stacked on each other. This leads to greater lateral variation of seismic data. In addition, there have little wells in BZ area. It is unreliable to research distribution of main channel sand using these above materials. In this paper, the geological model is designed firstly with the typical three-dimensional profile in BZ area. The petrophysical parameters are then obtained from wells which have been drilled. The results of forward modeling are finally studyed and analyzed to clear the seismic response characteristics of main channel sand. Based on the above processes, the distribution of the main channel sand is tracked in whole BZ area.

forward modeling; main channel; sand identification; BZ area

2014-08-12 改回日期：2015-03-26

肖偉(1981-)，男，工程師，現從事地震資料解釋和儲層預測研究，E-mail: diliaha@163.com。

1001-1749(2015)04-0488-06

P 631.4

A

10.3969/j.issn.1001-1749.2015.04.12