張秉銘 ， 韓 磊， 劉俊州， 劉振峰， 王勁松

(1.中國(guó)石化油田勘探開(kāi)發(fā)事業(yè)部，北京 100728；2.中國(guó)石化石油勘探開(kāi)發(fā)研究院，北京 100083)

井震結(jié)合速度建模在水平井軌跡設(shè)計(jì)中的應(yīng)用

張秉銘1， 韓 磊2， 劉俊州2， 劉振峰2， 王勁松2

(1.中國(guó)石化油田勘探開(kāi)發(fā)事業(yè)部，北京 100728；2.中國(guó)石化石油勘探開(kāi)發(fā)研究院，北京 100083)

分析了合成記錄標(biāo)定時(shí)深關(guān)系、地震速度分析和測(cè)井速度插值等目前常用的三種速度建模方法。這些方法分別在鉆井?dāng)?shù)量要求、構(gòu)造復(fù)雜程度和縱、橫向分辨率兼顧等方面存在一定局限性，限制了速度建模的精度。結(jié)合地震和測(cè)井速度資料各自優(yōu)勢(shì)，提出了一種新的在三維區(qū)域使用測(cè)井速度對(duì)地震速度進(jìn)行校正的方法：①?gòu)牡卣鹚俣润w中沿井軌跡抽取偽井速度曲線，并與測(cè)井速度比較得到殘差速度曲線；②在地質(zhì)格架約束下插值得到一個(gè)殘差速度體；③使用殘差速度體對(duì)地震速度進(jìn)行校正。將該方法應(yīng)用到川西水平井軌跡設(shè)計(jì)中，取得了較好的實(shí)際應(yīng)用效果。

速度建模； 時(shí)深轉(zhuǎn)換； 水平井軌跡設(shè)計(jì)

0 引言

速度預(yù)測(cè)一直是地球物理領(lǐng)域的難題之一。在前人持續(xù)不斷地努力下，形成了一系列速度建模的方法，這些方法概括起來(lái)分為3大類：

1)使用測(cè)井信息建立速度場(chǎng)(如井震標(biāo)定時(shí)深關(guān)系、測(cè)井速度空間插值和VSP等)[1-3]。這類方法建立的速度場(chǎng)在井點(diǎn)是比較準(zhǔn)確的，但在構(gòu)造復(fù)雜區(qū)域井間的橫向變化難以把握，只適合于構(gòu)造簡(jiǎn)單且鉆井較多的成熟開(kāi)發(fā)區(qū)塊。

2)使用地震處理速度，將地震疊加速度通過(guò)迪克斯(Dix)公式[4-7]或射線參數(shù)反演[8-12]轉(zhuǎn)為層速度用于時(shí)深轉(zhuǎn)換。Dix公式有水平層狀介質(zhì)射線垂直入射的假設(shè)，對(duì)于構(gòu)造復(fù)雜、各向異性強(qiáng)的區(qū)域會(huì)有較大誤差；射線參數(shù)反演具有更廣泛的地層產(chǎn)狀適應(yīng)性，但其反演是一種“逐層剝離”的過(guò)程，存在誤差逐層傳遞和放大的缺點(diǎn)。

3)使用鉆井或地質(zhì)信息約束地震速度建立速度場(chǎng)[13-20]，此類方法使用測(cè)井速度、分層對(duì)地震速度進(jìn)行了一定程度的校正，所得到的速度場(chǎng)相比前兩種方法更為準(zhǔn)確。但是也有兩方面局限性：①由于需要使用地質(zhì)分層進(jìn)行校正，對(duì)分層的準(zhǔn)確性要求較高；②只限于在井點(diǎn)或者二維平面對(duì)速度進(jìn)行校正，速度的空間變化特征難以有效把握。

1 常規(guī)速度建模方法

1.1 合成記錄標(biāo)定時(shí)深關(guān)系

利用合成記錄標(biāo)定得到的時(shí)深關(guān)系進(jìn)行時(shí)深轉(zhuǎn)換，是在實(shí)際工作中最常用的一種方法。合成記錄的目的是把時(shí)間域的地震記錄和深度域的測(cè)井資料建立對(duì)應(yīng)關(guān)系，通過(guò)井震標(biāo)定，可以得到一系列的“時(shí)深對(duì)”。這些“時(shí)深對(duì)”是由一些離散點(diǎn)組成的，通過(guò)軟件擬合就可以得到時(shí)深轉(zhuǎn)換公式。此方法的優(yōu)點(diǎn)是操作比較方便簡(jiǎn)單，對(duì)構(gòu)造起伏不大且鉆井?dāng)?shù)量較多的地區(qū)準(zhǔn)確度還是比較高的，但是如果地質(zhì)情況復(fù)雜、速度橫向變化大，將從井上得到的時(shí)深關(guān)系推廣應(yīng)用到整個(gè)研究區(qū)會(huì)產(chǎn)生很大的誤差。

1.2 測(cè)井速度插值

測(cè)井得到的速度在井點(diǎn)位置是比較準(zhǔn)確的，縱向分辨率也很高。如果鉆井在研究區(qū)分布均勻且數(shù)量較多，可以使用插值的方法建立三維速度場(chǎng)。在井間以距離作為權(quán)系數(shù)確定各井的速度數(shù)據(jù)在速度場(chǎng)中的貢獻(xiàn)值，一般只考慮距離因素，不考慮地質(zhì)結(jié)構(gòu)的影響，是一種純粹數(shù)學(xué)運(yùn)算的插值，與真實(shí)地層速度會(huì)有較大誤差。此方法得到的速度場(chǎng)模型化的痕跡較為嚴(yán)重，在井點(diǎn)周圍可能會(huì)出現(xiàn)明顯的異常區(qū)域，產(chǎn)生“牛眼”現(xiàn)象。其結(jié)果通常在井點(diǎn)位置是準(zhǔn)確的，但井間外推復(fù)雜構(gòu)造區(qū)域位置通常精度較差，而且由于使用的是速度曲線插值，對(duì)測(cè)井的聲波曲線質(zhì)量也要求較高。此類方法適用于開(kāi)發(fā)區(qū)塊測(cè)井?dāng)?shù)量較多且分布比較均勻的情況，對(duì)鉆井?dāng)?shù)量較少的區(qū)域，得到的速度場(chǎng)精度不高。

1.3 地震處理速度建模

地震速度的特點(diǎn)是分布均勻，含有較為豐富的橫向信息，并可根據(jù)需要進(jìn)行加密，形成一個(gè)全區(qū)的三維速度場(chǎng)。適合于勘探初期沒(méi)有井或者井很少的情況，由于布置測(cè)線時(shí)主測(cè)線會(huì)垂直于構(gòu)造走向，因此，一般情況下疊加速度的橫向變化能較好地反映勘探區(qū)的構(gòu)造形態(tài)。在實(shí)際應(yīng)用過(guò)程中，由于地震速度在拾取時(shí)會(huì)存在一定的異常點(diǎn)，所得到的速度場(chǎng)也只是在部分CDP點(diǎn)上有分布，因此需要對(duì)地震速度進(jìn)行散點(diǎn)校正消除異常值，并在縱橫向上進(jìn)行插值平滑處理。使用地震速度進(jìn)行時(shí)深轉(zhuǎn)換，大的構(gòu)造格架不會(huì)有問(wèn)題，但是由于速度譜資料是每隔幾百米拾取一個(gè)CDP點(diǎn)，局部速度變化被忽略了，用這種方法建立的速度場(chǎng)不夠精細(xì)。此外，常規(guī)地震速度在大偏移距、地層界面傾斜、鹽下等情況下的精度都不高[21]，而且由于速度譜在拾取過(guò)程中很大程度依賴于處理人員對(duì)區(qū)塊的整體地質(zhì)認(rèn)識(shí)，人為主觀影響比較大，這在一定程度上也會(huì)影響到時(shí)深轉(zhuǎn)換的精度。

2 井震殘差校正速度建模方法

2.1 方法流程

將單模透鏡光纖分別固定于三維精密調(diào)節(jié)平臺(tái)上,在微環(huán)腔的輸入輸出端進(jìn)行近似垂直光柵耦合,在實(shí)驗(yàn)過(guò)程中通過(guò)實(shí)時(shí)觀察CCD成像,不斷調(diào)節(jié)光纖位置,以達(dá)到最優(yōu)耦合。耦合信號(hào)通過(guò)光纖導(dǎo)出被光電探測(cè)器接收后由示波器顯示,最終得到微環(huán)諧振系統(tǒng)的透射譜線。

上述的三種速度求取方法是目前較為常用的，但由于這些方法是單一地使用地震或測(cè)井速度進(jìn)行速度建模，不可避免地存在一定的局限性。為了提高速度建模的精度，并盡可能地將測(cè)井、地震兩種速度在縱橫向上的優(yōu)勢(shì)結(jié)合起來(lái)，提出了一種新的使用井震三維殘差速度校正地震速度來(lái)建立速度場(chǎng)的方法。圖1是該方法計(jì)算的流程圖，具體步驟如下：

1)數(shù)據(jù)準(zhǔn)備。需要2類基礎(chǔ)數(shù)據(jù)：①地震處理過(guò)程中生成的疊加速度譜；②一定數(shù)量的有聲波時(shí)差曲線的井。目前大部分研究區(qū)都能滿足這兩個(gè)條件，對(duì)數(shù)據(jù)的要求不高。

2)疊加速度譜處理。地震疊加速度譜是人工拾取得到的，在拾取過(guò)程中由于每個(gè)人對(duì)區(qū)塊地質(zhì)認(rèn)識(shí)及拾取習(xí)慣不同，會(huì)存在一些誤差，需要對(duì)其進(jìn)行散點(diǎn)校正，結(jié)合研究區(qū)地質(zhì)認(rèn)識(shí)去掉一些不合理的樣點(diǎn)。此外，由于疊加速度譜的樣點(diǎn)是比較稀疏的，每隔幾百米才有一個(gè)點(diǎn)，需要對(duì)其進(jìn)行插值，得到三維速度場(chǎng)。最后，使用中值濾波等方法進(jìn)行縱橫向的平滑處理，使其看起來(lái)更光滑連續(xù)。圖2是散點(diǎn)校正前后速度譜對(duì)比圖，經(jīng)過(guò)散點(diǎn)校正后，去掉了一些異常的速度點(diǎn)。圖3是平滑前后的速度對(duì)比圖，平滑后消除了一些不連續(xù)點(diǎn)，使速度場(chǎng)看起來(lái)更連續(xù)自然。

3)疊加速度轉(zhuǎn)為層速度。地震處理后得到的速度一般是疊加速度，而測(cè)井速度是層速度，為了使二者匹配起來(lái)，需將疊加速度轉(zhuǎn)為層速度。一般使用Dix公式進(jìn)行轉(zhuǎn)化，其公式如式(1)所示：

(1)

其中：Vn為層速度；t0,n和Vδ,n分別為第n層的傳播時(shí)間和疊加速度。

4)層速度處理。Dix公式有一個(gè)基于水平層狀介質(zhì)的假設(shè)，對(duì)于復(fù)雜介質(zhì)在轉(zhuǎn)化過(guò)程中會(huì)產(chǎn)生一定的誤差，需要對(duì)層速度進(jìn)行校正處理來(lái)減小誤差。此處層速度已經(jīng)是三維體了，需要先將其進(jìn)行抽稀，然后再次進(jìn)行散點(diǎn)校正和平滑處理，其方法和流程2)類似。

圖1 井震殘差校正建立速度場(chǎng)的技術(shù)流程圖Fig.1 The flowchart of velocity model building through well seismic combination

圖2 速度譜散點(diǎn)校正前、后之間的比較Fig.2 Comparisons before and after velocity spectrum correction(a)散點(diǎn)校正前；(b)散點(diǎn)校正后

圖3 速度場(chǎng)平滑前、后之間的比較Fig.3 Comparisons before and after velocity smoothing(a)速度場(chǎng)平滑前；(b)速度場(chǎng)平滑后

5)原始聲波曲線預(yù)處理。聲波曲線是測(cè)井速度的來(lái)源，且對(duì)井震標(biāo)定影響較大，有必要對(duì)其進(jìn)行處理以提高質(zhì)量。影響聲波曲線質(zhì)量的因素主要有2個(gè)：①由于不同系列的測(cè)井儀器之間存在的系統(tǒng)誤差，導(dǎo)致不同井之間相同測(cè)井曲線在標(biāo)準(zhǔn)層的測(cè)井響應(yīng)存在差異；②由于井眼擴(kuò)徑或垮塌導(dǎo)致的曲線測(cè)量不準(zhǔn)。這兩個(gè)問(wèn)題可以使用多井一致性處理和環(huán)境校正來(lái)解決[22]。

6)建立地層格架模型。常規(guī)使用測(cè)井速度進(jìn)行插值一般不考慮地質(zhì)構(gòu)造的影響，是一種純粹數(shù)學(xué)運(yùn)算的插值。為了克服這一缺點(diǎn)，使用地層模型約束進(jìn)行井速度插值。對(duì)聲波曲線進(jìn)行校正之后，在井震標(biāo)定的基礎(chǔ)上對(duì)研究區(qū)主要目的層進(jìn)行層位解釋，然后使用解釋層位建立地層格架模型。該模型既含有井震標(biāo)定后生成的時(shí)深關(guān)系信息，又在橫向上指示構(gòu)造走向信息，可以作為后續(xù)三維速度建模的約束信息。

7)將步驟5)里經(jīng)過(guò)預(yù)處理的測(cè)井聲波時(shí)差曲線轉(zhuǎn)為速度曲線(此速度曲線是一個(gè)速度隨著時(shí)間變化的函數(shù)V1=f1(t))；從步驟4)生成的層速度體里沿井軌跡抽取偽井速度曲線(此曲線不是真實(shí)存在的，而是沿著已知井的井軌跡將地震層速度的值賦予井速度的概念，這樣就生成了一個(gè)地震速度曲線V2=f2(t))；將二者相減得到殘差速度曲線ΔV=V1-V2。最后在步驟6)生成的地層格架約束下對(duì)殘差速度曲線插值得到一個(gè)殘差速度體。

8)將步驟4)得到的層速度體和步驟7)得到的殘差速度體相加，即可得到最終速度體。此速度體在井點(diǎn)和測(cè)井速度吻合，在橫向上保留了地震速度的相對(duì)變化關(guān)系，是一個(gè)更加準(zhǔn)確的速度場(chǎng)。

2.2 方法優(yōu)勢(shì)

筆者提出的井震殘差校正三維速度建模方法，相比常規(guī)方法具有一定創(chuàng)新性。常規(guī)方法大多單一地使用地震、鉆井速度建立速度場(chǎng)，或者只是在井點(diǎn)或二維平面使用鉆井速度或地質(zhì)分層進(jìn)行校正。方法是在地層格架約束下在整個(gè)三維空間對(duì)地震速度進(jìn)行校正。(目前在專利和文獻(xiàn)中尚沒(méi)有出現(xiàn)此種方法)，該方法不僅適用于勘探程度低的區(qū)域，同時(shí)也適用于勘探程度較高區(qū)域，而且隨著勘探程度的提高、井?dāng)?shù)增加，其精度也會(huì)逐漸提高。能夠克服研究區(qū)地下結(jié)構(gòu)復(fù)雜、地層速度橫向變化劇烈等問(wèn)題。

3 應(yīng)用效果

速度場(chǎng)可以用于水平井軌跡設(shè)計(jì)、構(gòu)造圖編制和道集轉(zhuǎn)換等眾多方面。新馬-什邡地區(qū)處于四川盆地川西坳陷中段，隨著該區(qū)勘探開(kāi)發(fā)的腳步加快，越來(lái)越多的水平井被投入開(kāi)發(fā)。由于研究區(qū)儲(chǔ)層非均質(zhì)性較強(qiáng)，相鄰井甚至同一個(gè)井區(qū)其儲(chǔ)層厚度、泥巖夾層厚度變化較大；此外，研究區(qū)北部和東部的構(gòu)造起伏也比較大。在實(shí)際鉆進(jìn)過(guò)程中，稍有不慎，就會(huì)鉆到泥巖夾層里，需要不斷調(diào)整鉆頭位置，施工的難度很大。如果能得到深度域的數(shù)據(jù)體，就可以直接預(yù)測(cè)儲(chǔ)層頂、底深度，這樣會(huì)極大地降低施工難度。將儲(chǔ)層預(yù)測(cè)數(shù)據(jù)體從時(shí)間域轉(zhuǎn)換到深度域，最關(guān)鍵的是得到準(zhǔn)確的速度場(chǎng)。

圖4為在研究區(qū)使用井震殘差校正速度建立速度場(chǎng)的過(guò)程圖，其中圖4(a)是經(jīng)過(guò)插值、平滑后的疊加速度場(chǎng)，圖4(b)是從疊加速度轉(zhuǎn)化而來(lái)的層速度場(chǎng)，轉(zhuǎn)化之后同樣需要進(jìn)行平滑處理，圖4(c)中的藍(lán)色曲線是經(jīng)過(guò)濾波的測(cè)井速度曲線、紅色曲線是沿著井軌跡從層速度體中抽取的偽井速度曲線，兩種曲線比較得到一條殘差速度曲線，然后在地層格架的約束下對(duì)殘差曲線插值得到一個(gè)殘差速度體(圖4(d))，最后將殘差速度體和層速度體相加得到最終速度體(圖4(e))。可以看到原始的地震速度體(圖4(b))在橫向上包含一定的地層變化信息，但在縱向上分辨率較低，且通過(guò)兩種速度曲線比較可知，地震速度在井點(diǎn)和測(cè)井速度存在一定差異。在經(jīng)過(guò)三維校正之后，顯著提高了速度在縱向上的分辨率，既保證了在縱向上和測(cè)井速度吻合，同時(shí)在橫向上保留了地震速度的相對(duì)變化關(guān)系，是一個(gè)更加準(zhǔn)確的速度場(chǎng)。

為了評(píng)價(jià)不同方法建立速度場(chǎng)的準(zhǔn)確性，使用這些速度將時(shí)間域波阻抗體轉(zhuǎn)到深度域，然后將研究區(qū)的水平井投影在波阻抗剖面圖上。圖5為一口水平井在不同速度轉(zhuǎn)換的深度域波阻抗剖面上的投影，其中圖5(a)為井震標(biāo)定時(shí)深關(guān)系轉(zhuǎn)深度域波阻抗的結(jié)果，可以看到其相對(duì)井軌跡偏上，誤差在15 m左右；圖5(b)和圖5(c)分別使用的是地震速度和井插值速度，二者都相對(duì)井軌跡偏下，誤差在10 m左右；圖5(d)是使用殘差速度建模得到的速度轉(zhuǎn)深，其和井軌跡基本吻合，誤差較小。對(duì)研究區(qū)10口水平井進(jìn)行了A靶點(diǎn)深度的誤差統(tǒng)計(jì)(表1)。由表1可以看到，常規(guī)的三種方法平均絕對(duì)值誤差都超過(guò)了10 m，而使用本文方法的A靶點(diǎn)平均絕對(duì)值誤差為5 m左右，顯著提高了預(yù)測(cè)精度。

圖4 井震結(jié)合建立速度場(chǎng)技術(shù)的關(guān)鍵流程示意圖Fig.4 The key flowchart of velocity model building through well-seismic combination(a)疊加速度場(chǎng)；(b)層速度場(chǎng)；(c)測(cè)井、地震兩種速度曲線；(d)殘差速度場(chǎng)；(e)最終速度場(chǎng)

圖5 不同方法建立速度場(chǎng)時(shí)深轉(zhuǎn)換效果對(duì)比示意圖Fig.5 Time-depth conversion result comparisons between different velocity modeling methods(a)井震標(biāo)定時(shí)深關(guān)系轉(zhuǎn)深；(b)地震速度轉(zhuǎn)深；(c)測(cè)井插值速度轉(zhuǎn)深；(d)井震結(jié)合速度轉(zhuǎn)深表1 不同方法預(yù)測(cè)水平井A靶點(diǎn)深度誤差統(tǒng)計(jì)表Tab. 1 Comparisons of depth prediction error from different velocity modeling methods

井名實(shí)測(cè)A靶點(diǎn)深度/m井震標(biāo)定時(shí)深關(guān)系地震速度井插值速度殘差速度校正預(yù)測(cè)A靶點(diǎn)深度/m誤差/m預(yù)測(cè)A靶點(diǎn)深度/m誤差/m預(yù)測(cè)A靶點(diǎn)深度/m誤差/m預(yù)測(cè)A靶點(diǎn)深度/m誤差/mA114381425131430814251314335A211571147101169-1211451211525A31373.41380-6.61368.84.61380.4-71375.9-2.5A41436.4142016.41424.4121448.4-121429.17.3A51468.4145213.41480-11.614644.41462.63.8A61352.5134012.51365-13.51360-7.513511.5A7143414142014231114301414286A810171006111001161005121023-6A914331415181455-2214201314276A1012721293-11125418126571280-8平均絕對(duì)值誤差13.212.8710.195.11

4 結(jié)論

相比常規(guī)速度建模方法，筆者提出的井震殘差校正速度建模方法具有三個(gè)優(yōu)點(diǎn)：①不限制鉆井?dāng)?shù)量，幾口即可，多多益善；②由于使用構(gòu)造格架約束建模，對(duì)復(fù)雜構(gòu)造區(qū)域也有較好的適應(yīng)性；③保證了在縱向上和測(cè)井速度的吻合，又在橫向上保留了地震速度的相對(duì)變化關(guān)系。通過(guò)和常規(guī)方法在水平井靶點(diǎn)預(yù)測(cè)的應(yīng)用效果對(duì)比可知，提出的方法能建立一個(gè)更加準(zhǔn)確的速度場(chǎng)。由于速度預(yù)測(cè)是地球物理的一個(gè)難題，在水平井設(shè)計(jì)中可以使用幾種方法分別預(yù)測(cè)，然后將不同結(jié)果作為儲(chǔ)層的上下參考邊界。

[1] 凌云，郭建明，郭向宇，等．油藏描述中的井震時(shí)深轉(zhuǎn)換技術(shù)研究[J]．石油物探，2011，50(1)：1-13. LING Y，GUO J M，GUO X Y，et al．Research on time-Depth c-onversion by well-to-seismic in reserv-eoir characterization [J]．Geophysical Prospecting for Petroleum，2011，50(1)：1-13 .(In Chinese)

[2] 王西文，劉全新，蘇明軍，等．多井約束下的速度建模方法和應(yīng)用[J]．石油地球物理勘探，2003，38(03)：263-267. WANG X W，LIU Q X，SU M J,et al．M-ulti well constrained velocity model b-uilding method and its applicatio[J]．OGP，2003，38(3)：263-267.(In Chinese)

[3] 斯興焱，李錄明，胥良君，等．三維VSP多波速度分析方法及應(yīng)用[J]．石油物探，2012，51(1)：56-64. SI X Y，LI L M，XU L J，et al．3D-VSP multiple velocity analy-sis method and its application[J]．Geo-physical Prospecting for Petroleum，2012，51(1)：56-64.(In Chinese)

[4] 王磊．三維地震勘探中疊加速度成圖[J]．中國(guó)煤田地質(zhì)，2000( 2)：57-58. WANG L．Stacking velocity mapping in 3-D Seismic Exploration[J]．Coal g-eology of China，2000(2)：57-58.(In Chinese)

[5] DUCHESNE MJ，CLAPROOD M，GLOAGUEN E．Improving seismic velocity estimat-ion for 2d poststack time migration ofregional seismic data using kriging w-ith an external drift[J]．The Leading Edge，2012，31(10) ：1156-1166

[6] 李明娟，李守濟(jì)，牛濱華，等．地震速度譜在精細(xì)深度圖制作中的應(yīng)用[J]．石油物探，2004，43(3)：272-274. LI M J，LI S J，NIU B H，et al．The application of seismic velocity spectrum in making depth map[J]．GPP，2004，43(3)：272-274.(In Chinese)

[7] MARSDEN D，BUSH M D，JOHNG D S． Analytic velocity function[J]．The leading Edge，1995，14(7) ：775-782.

[8] 馬濤，嚴(yán)又生，劉文利．塔里木盆地速度場(chǎng)的建立及應(yīng)用[J]．石油地球物理勘探，1996，31(3)：382-393. MA T，YAN Y S，LIU W L．Det-ermination of velocity field in the Ta-rimbasin and its application[J]．OGP，1996，31(3)：382-393.(In Chinese)

[9] WHITMORE N D，GARING J D．Interval velocity estimation using iterative pre-stack depth migration in the constant angle domain．The Leading Edge，1993，12(7) ：757-762.

[10]陳業(yè)全，劉春曉，王道義．塔里木盆地中原探區(qū)速度場(chǎng)研究方法[J]．地球物理學(xué)進(jìn)展，2004，19(3)：656-662. CHEN Y Q，LIU C X，WANG D Y．The research method of velocity field in central plains prospect area oftarim basin[J]．Process In Geophysics，2004，19(3)：656-662.(In Chinese)

[11]馬海珍，雍學(xué)善，楊午陽(yáng)，等．地震速度場(chǎng)建立與變速構(gòu)造成圖的一種方法[J]．石油地球物理勘探，2002，31(01)：53-59. MA H Z，YONG X S，YANG W Y．A method for building of seismic velocity field and structure maping by variation in velocity[J]．OGP，2002，37(1) ：53-59.(In Chinese)

[12]王樹(shù)華，劉懷山，張?jiān)沏y，等．變速成圖方法及應(yīng)用研究[J]．中國(guó)海洋大學(xué)學(xué)報(bào)，2004，34(1)：139-146. WANG S H，LIU H S，ZHANG Y Y,et al. Research on variable-velocitystructure Mapping and its Application[J]．Journal of ocean university of Ch-ina，2004，34(1)：139-146 .(In Chinese)

[13]王華忠，馮波，李輝，等．各種速度分析與反演方法的對(duì)比研究[J]．巖性油氣藏，2012，24(5)：1-11. WANG H Z，F(xiàn)ENG B，LI H， et al．Comparison among velocity an-alysis and inversion methods[J]．Litho-logic Reservoirs，2012，24(5)：1-11.(In Chinese)

[14]賈義蓉，賀振華，石蘭亭，等．變速三維地震速度場(chǎng)的構(gòu)建與應(yīng)用[J]．物探化探計(jì)算技術(shù)，2011，22(3)：243-247. JIA Y R，HE Z H，SHI L T，et al．Structure mapping with transver-sely variable 3D seismic velocity mo-del in Erlian basin[J]．Compating techniques for geohysical and geochemical exploration,2011，22(3)：243-247.(In Chinese)

[15]膠存禮，周杰，尚雅珍．變速成圖系統(tǒng)技術(shù)關(guān)鍵及其在塔中西部的應(yīng)用[J]．物探與化探，2003，27(3):189-194. JIAO C L，ZHOU J，SHANG YZ． The varying velocity mapping techniq-ue：Its key problems and its applicati-on to western Tazhong area[J]．Compating technipues for geo;hysical and geochemical exploration,2003，27(3)：189-194.(In Chinese)

[16]王香文．東嶺地區(qū)三維速度模型的建立和應(yīng)用[J].勘探地球物理進(jìn)展，2006，29(6)：413-418. WANG X W．3-D velocity model building and its application in Dongling area[J]．PEG，2006，29(6)：413-41.(In Chinese)

[17]汪功懷，秦廣勝，蔡其新，等．東濮凹陷地震速度場(chǎng)建立方法與應(yīng)用研究[J].勘探技術(shù)，2011(2)：58-66. WANG G H，QIN G S，CAI Q X，et al．A method for establishing seismic velocity field and Its application in Dongpu Sag[J].China petroleum exploration，2011(2)：58-66.(In Chinese)

[18]邊立恩，于茜，谷志猛，等．低速異常帶平均速度場(chǎng)建模方法研究[J]．巖性油氣藏，2015，27(3)：122-126. BIAN L E，YU Q，GU Z M，et al．Research on mean velocity field modeling method for low velocity anomaly zone：A case study from Q oilfield，Bohai[J]．Lithologic Reservoirs，2015，27(3)：122-126.(In Chinese)

[19]吳超，寧斌，王曉云．優(yōu)化三維層位控制法速度建場(chǎng)及其應(yīng)用[J]．新疆石油地質(zhì)，2003，24(5)：458-460. WU C ，NING B，WANG X Y． Developing seismic velocity field by optimized 3D horizon control method and its application[J]．XinJiang Petrol-eum Geology，2003，24(5)：458-460.(In Chinese)

[20]KEHO T，SAMSU D．Depth conversion of Tangguh gas field[J]．The Leading Edge，2002，21(10)：966-971.

[21]潘宏勛，方伍寶．地震速度分析方法綜述[J]．勘探地球物理進(jìn)展，2006，29(5)：305-311. PAN H X，F(xiàn)ANG W B．The su-mmarize of seismic velocity analysis[J]．PEG，2006，29(5)：305-311.(In Chinese)

[22]苗凡勝．測(cè)井曲線預(yù)處理技術(shù)在時(shí)深標(biāo)定中的應(yīng)用[J]．長(zhǎng)江大學(xué)學(xué)報(bào)(自科版)，2014，11(13)：57-60. MIAO F S．Application of Pretre-atment Technology of Logging Curves in Time-depth Delimitation[J]．Journal of Yangtze University(Natural Science Edition)，2014，11(13)：57-60.(In Chinese)

Well-seismic velocity model building method and its application to horizon well deployment

ZHANG Bingming1， HAN Lei2， LIU Junzhou2， LIU Zhenfeng2， WANG Jinsong2

(1.Oil Exploration and Development Division of SINOPEC,Beijing 100728,China；2.Oil Exploration and Development Research Institute of SINOPEC,Beijing 100083,China)

The traditional three velocity model construction methods, including 1) by estimating from calibrated Time-Depth relationships, 2) by extracting from the seismic velocity and 3) by interpolating the well logging data, are often limited by the drilling/logging insufficiency, the defective description of the complicated geologic structure, as well as the imperfect resolution in both lateral and vertical direction. In this paper, we proposed a velocity construction method through spatially calibrating seismic velocity using well velocity. Firstly, the pseudo velocity curves along each wellbore trajectory extracted from the seismic velocity are compared with the well logging velocity. Then the residual velocity can be achieved, and a residual velocity volume can be computed via the interpolation under the constraint of geological framework. Finally, this volume is further employed to correct the seismic velocity and to produce a new velocity model. This method has been successfully applied in the horizontal well deployment in the Chuanxi area.

velocity model； time-depth conversion； horizontal well deployment

2017-05-10 改回日期：2017-05-20

國(guó)家科技重大專項(xiàng)(2016ZX05002-005)

張秉銘(1963-)，男，博士，主要從事各向異性介質(zhì)地震波場(chǎng)理論、油氣藏地球物理應(yīng)用研究， E-mail：zhangbm@sinopec.com。

韓磊(1984-)，男，工程師，主要從事地震資料解釋與儲(chǔ)層預(yù)測(cè)， E-mail：hanlei.syky@sinopec.com。

1001-1749(2017)04-0551-07

P 631.4

A

10.3969/j.issn.1001-1749.2017.04.16