徐群洲

(新疆油田分公司勘探開發(fā)研究院 地球物理研究所，烏魯木齊 830013)

?

準噶爾盆地地震速度場的建立

徐群洲

(新疆油田分公司勘探開發(fā)研究院 地球物理研究所，烏魯木齊 830013)

分析盆地內(nèi)各層巖性速度特征，在準噶爾盆地里建立一個盆地級巨型速度場(面積約13×104km2，深度達8 km)。面對海量的地震資料，首先進行盆地地震剖面的解釋，建立盆地17層構(gòu)造解釋時間模型，利用疊加速度譜分析、VSP井速度標定、聲波測井的VSP轉(zhuǎn)換技術(shù)、測線交點疊加速度自動校正閉合差、野值自動剔除及井約束層速度分析，來研究各沉積層速度特征，模型模擬層速度建場、最終建立盆地平均速度場。用于區(qū)域?qū)铀俣确治?，?gòu)造轉(zhuǎn)深和地震資料處理的速度參考，用兩口VSP井實測速度資料進行驗證，盆地速度場平均速度絕對誤差平均為35 m/s, 相對誤差平均為1.8%,兩者大小和變化特征一致，表明盆地速度場精度較高。

地震剖面解釋； VSP井測試； 層速度； 三維速度場

0 前言

準噶爾盆地屬于大致呈三角形的內(nèi)陸盆地，面積約13×104km2。盆地構(gòu)造樣式多，巖性變化大。隨著勘探程度的提高對構(gòu)造圖的精度要求赿來赿高，單塊三維和小區(qū)域二維的地震速度研究由于資料少，邊界效應等因素影響，層速度變化趨勢很難準確分析，這就需要建立一個大范圍的地震速度場，研究盆地各層巖性速度的橫向變化特征，指導小區(qū)域的精細地震速度研究，建立高精度的速度場，提高變速成圖的精度。

一個盆地的形成，它的地震速度也就是唯一的，通過不斷的地震勘探和鉆井資料的積累，就能得到逐漸精確的盆地速度場。建立巨型盆地級速度場存在大量的問題，如盆地地質(zhì)層位的時間剖面的解釋，基準面、測線交點的速度誤差，疊加速度野值的剔除，VSP井分布不均、密度小、鉆井井淺、測線疊加速度的標定，以及海量數(shù)據(jù)的存儲、統(tǒng)一調(diào)用計算；深層鉆井少，深層層速度標定困難等。作者采用統(tǒng)一固定基準面、二維速度閉合差校正技術(shù)，自動去除野值、VSP井速度約束、聲波測井資料彌補VSP井資料的不足，三維空間網(wǎng)格化技術(shù)來校正鉆井和疊加速度的基準面，分析各層層速度資料，建立了盆地級17層層速度場，解決了上述的問題[1]。

1 原始資料分析

建立速度場所需的原始資料包括盆地時間層位解釋成果、疊加速度資料、VSP測井和聲波測井資料。VSP測井的頻帶與地震反射波法頻帶相同，因此由VSP測井得到的層速度和平均速度被認為是地層速度的真值。到2013年底新疆油田公司共作了365口井的VSP測井，平均350 km2才有一口VSP測井。

本次研究主要采用盆地內(nèi)二維測線的疊加速度資料，約94 865個。由疊加速度譜轉(zhuǎn)換得到的層平均速度是構(gòu)建速度場的基礎(chǔ)。用疊加速度譜計算的層速度精度比VSP測井與聲波測井計算的層速度精度要低，所以用VSP測井與聲波測井資料來約束校正速度譜計算的層速度，以提高速度場精度。另外聲波測井計算的層速度精度比VSP測井計算的層速度精度低，VSP測井計算的層速度精度最高，因此只用聲波測井計算的層速度彌補VSP測井密度不夠的地區(qū)的層速度校正計算。

2 綜合標定和地震解釋

盆地時間層位解釋主要采用盆地內(nèi)二維測線4 922條[2]。歷史上準噶爾盆地二維測線大多在浮動基準面上，首先數(shù)字化每條剖面的浮動基準面數(shù)據(jù)，將二維測線基準面統(tǒng)一調(diào)整到水平基準面1 500 m。通過地面露頭標定、標準井合成記錄標定引層區(qū)域統(tǒng)層大剖面9條，在標定好的基礎(chǔ)上進行層位和斷裂的解釋。再由這9條區(qū)域統(tǒng)層大剖面向其他剖面引層標定層位，進行層位和斷裂的精細解釋，根據(jù)鉆井資料對地震剖面時間控制層位進行分析解釋，分別解釋了N1s(沙灣組)、E1-2z(紫泥泉子組)、K2d(東溝組)、K1h(呼圖壁河組)、K1q(清水河組)、J3q(齊古組)、J2t(頭屯河組)、J2x(西山窯組)、J1s(三工河組)、J1b(八道灣組)、T3b(白堿灘組)、T2k(克拉瑪依組)、T1b(百口泉組)、P3w(上烏爾禾組)、P2p(平地泉組)、P1j(金溝組)共17層，獲得盆地各層的橫向變化規(guī)律及空間展布特征[1]。通過地震剖面時間控制層位解釋成果，建立三維t0構(gòu)造模型[3]。

3 地震速度場建立方法

3.1 模型迭代法求取平均速度

由疊加速度譜轉(zhuǎn)換得到的平均速度是構(gòu)建速度場的基礎(chǔ)，它的精度直接關(guān)系到速度場的精度。模型迭代法首先給出層速度模型，用射線追蹤法求出反射波時距曲線，然后模擬做速度譜的方法，用理論曲線對它進行擬合，求出疊加速度，把它和速度譜上的疊加速度進行比較，當絕對誤差小于5%時，認為所給的模型是正確的，輸出層速度值，轉(zhuǎn)換成平均速度。

準噶爾盆地的滴西12井，井深4 260 m，井位上有一個疊加速度譜(L200603測線，樁號2050)。使用模型迭代法計算出的平均速度與滴西12井的VSP測井的平均速度進行比較(圖1)，絕對誤差為21 m/s～190 m/s, 相對誤差為0.1%～9.5%；準噶爾盆地南緣山前高陡構(gòu)造上的高泉1井上正好有一個疊加速度譜(AG201113測線，樁號750)，井深為4 660 m。使用模型迭代法計算出的平均速度與高泉1井的VSP測井的平均速度進行比較(圖2)，絕對誤差為200 m/s～250 m/s, 相對誤差為20%～25%，全井段誤差大小基本一致[4]。

圖1 疊加速度譜模型法轉(zhuǎn)換平均速度與VSP井實測平均速度對比

圖2 盆地速度場平均速度與VSP井實測平均速度對比

對盆地內(nèi)VSP井點處的地震速度進行比較研究，結(jié)果表明，在盆地腹部地層比較平緩的地區(qū)用模型迭代法計算出疊加速度譜的平均速度，誤差最大在10%。在地層傾角比較大的區(qū)域或上下層地層傾角差異比較大的區(qū)域?qū)ΟB加速度譜的誤差更大，需要進行較正[5]。

3.2 速度資料分析

3.2.1 疊加速度譜去野值

疊加速度譜在處理點譜時、人工錄入數(shù)據(jù)庫的過程中常存在一些野值，使得層速度網(wǎng)格化時出現(xiàn)異常。但盆地內(nèi)有9.5×104個疊加速度譜，人工手動剔除難以實現(xiàn)，因此采用編程自動去野值。自動去野值程序采用圓半徑內(nèi)數(shù)據(jù)大于或小于圓半徑內(nèi)平均速度值的多少判定為野值，進行取舍，這個值一般定為10%。

3.2.2 測線交點疊加速度自動校正閉合差

由于野外地震數(shù)據(jù)采集參數(shù)不同、資料處理參數(shù)不同、采集時間跨度大，以及測量、地層傾角、施工誤差和靜校正的影響等，準噶爾盆地二維測線交點上不同測線的速度不一致，交點速度閉合差一般為20 m/s,最高可達200 m/s。在測線交點速度閉合運算時，先計算兩條測線交點處的平均速度—深度曲線，對兩條速度曲線的同一深度段，計算其平均時差，由其可靠性高低確定加權(quán)系數(shù)，作為該層閉合差[6]。

3.2.3 井中速度數(shù)據(jù)的處理

VSP速度和聲波測井速度數(shù)據(jù)都是井中速度資料，其中VSP速度資料最為重要，是層速度較正的標準[7]。

1)統(tǒng)一速度分析的基準面校正，就是將三種速度資料的時深零都統(tǒng)一到一個基準面上，定為準噶爾盆地浮動基準面。歷史上準噶爾盆地二維測線都在浮動基準面上，首先數(shù)字化每條剖面的浮動基準面數(shù)據(jù)，再對數(shù)據(jù)進行平面網(wǎng)格化，形成浮動基準面[8]。將位于地表高程的VSP資料校正到浮動基準面上。

2)VSP層速度提取。對365口井VSP的速度資料進行分析，計算出平均速度隨深度變化特征曲線，利用鉆井分層數(shù)據(jù)，按每一層位的頂界深度和底界深度計算出VSP井對應地層的層平均速度，用于層速度校正[9]。

3)用聲波測井速度資料層速度提取。在VSP井密度小對疊加速度譜約束力度不夠的地區(qū)，增選了400口井的聲波測井資料，做為層速度補充校正。但聲波測井很少測試到井口資料，缺少地面和淺層旅行時，就不能得到平均速度。如把聲波測井速度校正轉(zhuǎn)換成VSP速度，就可當成偽VSP井數(shù)據(jù)應用[10]。在計算時在盆地不同構(gòu)造單元選取VSP井的聲波數(shù)據(jù)，用同一口井的VSP數(shù)據(jù)和聲波數(shù)據(jù)進行統(tǒng)計分析，淺層旅行時用線性預測法進行計算，聲波漂移量用分層對比差值法計算，這樣就得到了一個小構(gòu)造單元的校正經(jīng)驗值。用在本區(qū)域其他井聲波測井數(shù)據(jù)的VSP校正轉(zhuǎn)換，就將聲波速度數(shù)據(jù)轉(zhuǎn)制成偽VSP數(shù)據(jù)。用校正后的聲波數(shù)據(jù)得到的速度，與VSP速度一同來標定疊加速度[11]。

3.3 井約束層速度校正

層速度數(shù)據(jù)分析時，VSP、聲波數(shù)據(jù)、疊加速度的時深零線都校正到了浮動基準面上，由疊加速度得到的層速度要用VSP和聲波測井的偽VSP速度進行標定[12]。

對每一沿層的層速度進行VSP和聲波層速度誤差校正(井約束層速度分析)，這個誤差是坐標函數(shù)Δv(x,y,z),即不同深度、不同坐標的校正量不同。在每口井點上疊加速度的層速度vsp(x,y,z)與VSP層速度vvsp(x,y,z)之差，即井位上的層速度校正量，寫成

Δvi(x,y)=vsp(x,y,z)-vvsp(x,y,z)

根據(jù)該層上各井的層速度誤差校正量Δvi(x,y)作曲面網(wǎng)格擬合，得到誤差網(wǎng)格曲面，便可以得到該層上無井處的速度校正量。

圖3 速度場應用效果實例

3.4 模擬迭代層速度建場

通過井約束層速度分析，確定各層層速度變化特征和對應17層時間層位的地質(zhì)模型。其中白堊系清水河組層速度，在昌吉凹陷處層速度為一個高速區(qū)，最大為5 200 m/s，向四周斜坡區(qū)層速度逐漸減小，最小值為3 200 m/s，烏倫古凹陷處有一個高速平臺，橫向上層速度受埋深控制，層速度變化趨勢和時間層位變化趨勢一致[13]。有了每一地質(zhì)層位的層速度數(shù)據(jù)，用盆地17層時間解釋模型構(gòu)建對應的三維空間層速度模型。模型模擬方法由用理論道集和實際地質(zhì)模型得到理論時距曲線，再根據(jù)理論時距曲線求取疊加效果最好的理論疊加速度，直到兩者相吻合，分析共反射道集中每一道的疊加情況反復迭代求出每一層的層速度。這樣地震速度模型和時間層位模型對應一致，附合地震波實際的傳播路徑。

4 速度場的精度和應用效果

盆地的速度場完成后，用鉆井比較深的高泉1井和最新的沙排1井VSP井資料進行驗證，用速度場中井位上的平均速度與各井的VSP實測的平均速度進行對比(圖1)[14]。在高泉1井處絕對誤差最大為85 m/s，最小為6 m/s，平均為36.5 m/s；相對誤差最大為5.0%，最小為0.2%，平均值為1.4%，兩者的平均速度變化趨勢基本一致。在沙排1井井位上絕對誤差最大為100 m/s，最小為16 m/s，平均為36.5 m/s；相對誤差最大為4.9%，最小為0.7%，平均值為1.7%，兩者的平均速度基本相同。對比結(jié)果顯示，淺層的誤差較大，主要是由于建場用的是充填速度2 000 m/s，而VSP速度由于淺層頻譜發(fā)散速度偏低且不準確，速度場的淺層速度更合理。中深部對比說明,速度場的精度很高。

地震速度研究的主要目的是變速成圖，由于地層巖層速度的變化能使地下構(gòu)造在地震時間剖面和等t0構(gòu)造圖上發(fā)生位移、畸變。速度場的重要作用就是在等深構(gòu)造圖上使這些被位移、畸變或消失的構(gòu)造真實地顯現(xiàn)出來。

陸南地區(qū)的A構(gòu)造是經(jīng)鉆井證實的一個低幅度背斜，在等t0構(gòu)造圖上不存在(圖3(a))，用本速度場提供的平均速度(圖3(b)),進行變速成圖，構(gòu)造轉(zhuǎn)深的等深度構(gòu)造圖上這個構(gòu)造顯現(xiàn)出來(圖3(c))。進一步證明建立高精度地震速度場的重要性，在勘探復雜構(gòu)造、大深度和低幅度構(gòu)造上的價值。

5 結(jié)論與建議

地震速度是地震勘探研究的靈魂，通過地質(zhì)資料、地質(zhì)構(gòu)造是無法和處理過程溝通的，要真正實現(xiàn)處理解釋一體化，只能以速度作為一個橋梁。建立盆地整體的地震速度場是基礎(chǔ)，特別是有了大范圍的層系的速度變化趨勢特征，不僅能為處理提供參考速度，也能為構(gòu)造成圖提供轉(zhuǎn)深速度。

[1] 雷德文，劉繼山，張越遷，等.構(gòu)造建模及三維可視化技術(shù)在準噶爾盆地的應用[J].新疆石油地質(zhì)，2012，33(5):589-591. LEI D W, LIU J SH, ZHANG Y Q, et al. Structural modeling and application of visualization technology in the junggar basin . Xinjiang Petroleum Geology, 2012,33 (5): 589-591.(In Chinese)

[2] 李航.南方陸相盆地復雜斷裂帶綜合解釋技術(shù)[J].大慶石油地質(zhì)與開發(fā)，2012，31(3):156-159. LI H. comprehensive interpretation technology of complicated fault zone in continental basins in the South . geology and development of Daqing petroleum, 2012,31 (3): 156-159.(In Chinese)

[3] 唐建人，李清仁，陳炳峰，等.變速校t0曲射線空間歸位成圖方法[J].大慶石油地質(zhì)與開發(fā)，1994，13(3):67-70. TANG J R, LI Q R, CHEN B F,et al. t0 speed school song-ray back into . Geology and development of Daqing petroleum, 1994,13 (3): 67-70.(In Chinese)

[4] 汪功杯，劉忠亮，李勤英，等.文明寨油田構(gòu)造精細解釋與儲層預測[J].石油物探，2006，45(1):94-95. WANG G B, LIU ZH L, LI Q Y,et al. Fine structural interpretation and reservoir prediction in wenmingzhai oilfield . Geophysical, 2006,45 (1): 94-95.(In Chinese)

[5] 朱廣生，楊金華，段天友，等.塔里木盆地巨型速度場的建立[J].新疆石油地質(zhì)，2002，23(6):507-510. ZHU G SH,YANG J H,DUAN T Y,et al. Established mega-velocity field in the Tarim basin . Xinjiang Petroleum Geology, 2002,23 (6): 507-510.(In Chinese)

[6] 李明娟，李守濟，牛濱華，等.地震速度譜在精細深度圖制作中的應用[J].石油物探，2004，43(3):272-274. LI M J,LI SH J,NIU B H, et al.Application of seismic velocity in the making of fine depth map . Geophysical, 2004,43 (3): 272-274.(In Chinese)

[7] 馬海珍，雍學善，楊午陽，等.地震速度場建立與變速構(gòu)造成圖的一種方法[J].石油地球物理勘探，2002，37(1)：53-59. MA H ZH,YONG X SH，YANG W Y,et al.Establishment of seismic velocity field and a method of variable structure mapping . Geophysical exploration for oil, 2002,37 (1): 53-59.(In Chinese)

[8] 范貴良.山前高陡構(gòu)造區(qū)地震成圖方法研究[J].石油物探，2000，39(3):16-25. FAN G L. Research on mountain seismic mapping method in high-steep structural area . Geophysical, 2000,39 (3): 16-25.(In Chinese)

[9] 李學義，李天明，王兵，等.準噶爾盆地低幅度背斜圈閉識別方法[J].新疆石油地質(zhì)，2005，26(2):199-201. LI X Y,LI T M,WANG B,et al. Identification method of low-amplitude anticline traps in Junggar basin . Xinjiang Petroleum Geology, 2005,26 (2): 199-201.(In Chinese)

[10]王興軍，彭更新，滿益志，等.復雜高陡構(gòu)造圈閉落實方法[J].新疆石油地質(zhì)，2009，30(2):258-260. WANG X J,PENG G X,MAN Y ZH,et al. Implementation of high steep complex structure trap. Xinjiang Petroleum Geology, 2009,30 (2): 258-260.(In Chinese)

[11]甘永忠，朱廣生.由疊加速度反演層速度[C].中國地球物理學會年刊，1992. GAN Y ZH,ZHU G SH.From stack velocity inversion layer speed . the Chinese Geophysical Society annual report, 1992.(In Chinese)

[12] 陳遵德，楊金華，楊舒，等.速度閉合差校正方法[J].石油地球物理勘探，2000，35(3):381-385. CHEN Z D,YANG J H,YANG SH,et al.Speed loop error correction methods . Geophysical Exploration for oil, 2000,35 (3): 381-385.(In Chinese)

[13]趙雨晴，張永華，羅家群，等.泌陽凹陷復雜斷裂帶三維地震資料精細解釋方法和效果分析[J].石油物探，2009，48(5):524-530. ZHAO Y Q,ZHANG Y H,LUO J Q,et al.Complex in Biyang fault zone fine interpretation method and effect analysis of three-dimensional seismic data . Geophysical, 2009,48 (5): 524-530.(In Chinese)

[14]汪功杯，劉忠亮，李勤英，等.文明寨油田構(gòu)造精細解釋與儲層預測[J].石油物探，2006，45(1):94-95. WANG G B,LIU ZH L,LI Q Y,et al.Fine structural interpretation and reservoir prediction in wenmingzhai oilfield . Geophysical, 2006,45 (1): 94-95.(In Chinese)

[15]楊勤林，王彥春，張靜，等.濱里海盆地M塊鹽下構(gòu)造變速成圖[J].石油物探，2012，51(4):377-381. YANG Q L,WANG Y CH,ZHANG J,et al.M block in the Caspian basin subsalt varying velocity mapping . Geophysical, 2012,51 (4): 377-381.(In Chinese)

The establishment of seismic velocity field in Junggar basin

XU Qun-zhou

(Institute of exploration and Development Research Institute of Xinjiang Oilfield Branch Ball Physics,Urumqi 830013,China)

By analysis of velocity characteristics of rock within the basin, it is established that a basin-level giant velocity field in the Junggar basin (an area of about 13x104 km2with depth of 8 km). With a great deal of seismic data, the interpretation was firstly carried out for the seismic profiles, and a structure of 17 layers was established in time domain. The velocity characterization of various sedimentary layers has been analyzed and the interval velocity fields have been setup by using stack velocity spectrum analysis, VSP well velocity calibration, sonic logging with VSP conversion technology, automatic mistie correction of line intersection and the wild value automatically removed and interval velocity analysis with well velocity constrains. The final average velocity field of whole basin has been built up for regional velocity analysis, time to depth conversion of structures and the reference velocity data for seismic data processing. The verification of the accuracy of the velocity field has been done with 2 VSP wells velocity data, the absolute error is averagely 35 m/s, the relative errors is averagely 1.8 %, and the both amount and characteristics of computing with measuring is coincide, which indicated the high precision of the basin velocity field.

seismic section interpretation; VSP well data; interval velocity; three-dimensional velocity field

2014-09-15 改回日期：2014-10-08

徐群洲(1967-)，男，高級工程師，長期從事油氣勘探石油地質(zhì)解釋工作，E-mail：xuqz@petrochina.com.cn。

1001-1749(2015)04-0502-06

P 631.4

A

10.3969/j.issn.1001-1749.2015.04.14