井控各向異性速度建模技術(shù)在YKL地區(qū)的應(yīng)用

裴云龍，王立歆，鄔達(dá)理，蔣 波，鄧 鋒

（1．中國(guó)石油化工股份有限公司石油物探技術(shù)研究院，江蘇南京211103；2．中國(guó)石油化工股份有限公司西北油田分公司，新疆烏魯木齊830011）

井控各向異性速度建模技術(shù)在YKL地區(qū)的應(yīng)用

裴云龍1，王立歆1，鄔達(dá)理1，蔣 波1，鄧 鋒2

（1．中國(guó)石油化工股份有限公司石油物探技術(shù)研究院，江蘇南京211103；2．中國(guó)石油化工股份有限公司西北油田分公司，新疆烏魯木齊830011）

YKL氣藏開發(fā)迫切需要高品質(zhì)地震成像資料的支撐，尤其是要求有高精度的構(gòu)造成像資料。由于該區(qū)塊三維地震勘探資料存在缺失近炮檢距、覆蓋次數(shù)低等不利因素，因此，常規(guī)單一的速度建模方法難以解決該區(qū)井震誤差大的問題。為此，針對(duì)性地提出并聯(lián)合應(yīng)用了井控各向異性網(wǎng)格層析和逐層速度建模技術(shù)。在選好區(qū)內(nèi)標(biāo)志層、建好初始速度模型的基礎(chǔ)上，利用區(qū)內(nèi)多口已知井信息求取各向異性參數(shù)，先消除淺部大套穩(wěn)定地層速度對(duì)氣藏深度的影響，再由淺到深逐層對(duì)標(biāo)志層反演，更新迭代各向異性參數(shù)、速度及深度成像，直至各標(biāo)志層的井震誤差降低到規(guī)定范圍內(nèi)。3口驗(yàn)證井的應(yīng)用結(jié)果顯示，實(shí)際資料最終成像成果的氣藏構(gòu)造深度預(yù)測(cè)值與實(shí)鉆井深度誤差均在5m內(nèi)，氣藏構(gòu)造成像的精度得到了顯著提高，斷裂成像也得到明顯改善，降低了開發(fā)鉆井的風(fēng)險(xiǎn)。

疊前深度偏移；井控；各向異性；網(wǎng)格層析；速度建模；井震誤差

隨著我國(guó)西部YKL氣藏開發(fā)的推進(jìn)，含水量上 升成為該氣藏開發(fā)的主要矛盾，氣藏開發(fā)難度越來越大。YKL氣藏歷年實(shí)鉆井?dāng)?shù)據(jù)顯示，氣藏構(gòu)造深度預(yù)測(cè)值與實(shí)鉆數(shù)據(jù)之間誤差較大，尤其在構(gòu)造邊部更是如此，地震勘探成果資料不能精確顯示構(gòu)造和斷裂位置是開發(fā)過程中面臨的主要問題。

用時(shí)間域偏移成像數(shù)據(jù)制作構(gòu)造圖，傳統(tǒng)的做法是變速時(shí)深轉(zhuǎn)換成圖［1－4］。利用DIX公式將均方根速度轉(zhuǎn)換成層速度，再將時(shí)間域?qū)游晦D(zhuǎn)換到深度域；或者利用區(qū)域綜合平均速度方法得到時(shí) 深關(guān)系，直接將T0圖轉(zhuǎn)換為深度域。受地層傾角和速度橫向變化等因素的影響，時(shí)間域偏移成像數(shù)據(jù)的成像精度不夠準(zhǔn)確。另外，因速度誤差以及變速成圖方法差異的影響，井震誤差往往也較大，準(zhǔn)確預(yù)測(cè)鉆探目標(biāo)深度難度很大［5－6］。

當(dāng)?shù)叵麓嬖谒俣犬惓ｓw時(shí)，射線經(jīng)過這些異常體將產(chǎn)生大的偏折，此時(shí)，若按疊前時(shí)間偏移（pre－stack time migration，PSTM）的共中心點(diǎn)道集進(jìn)行偏移疊加，顯然得不到真實(shí)的反射。疊前深度偏移的巨大優(yōu)勢(shì)在于，在共深度點(diǎn)道集中，如果速度模型正確，則對(duì)于某個(gè)深度點(diǎn)來說，其道集反射應(yīng)是水平的，由此我們可以判斷模型正確與否［7］。相對(duì)于PSTM，疊前深度偏移（pre－stack depth migration，PSDM）在成像精度等方面具有顯著的優(yōu)勢(shì)［8－12］，但因其是一個(gè)對(duì)速度很敏感的處理過程，往往會(huì)遇到建模難題［7］。

目前，三維疊前深度偏移速度模型建立的方法有多種，如：測(cè)井約束速度模型優(yōu)化速度建模方法［13］；初始速度模型的建立、模型的修改與驗(yàn)證及模型的修飾性處理和可視化顯示方法［14］；利用地震同相軸的多種信息進(jìn)行層析反演，并在反演中加入含有地質(zhì)含義的傾角信息以增加反演過程的穩(wěn)定性和收斂速度，利用類似于疊加速度譜的斜率速度譜來拾取地震同相軸斜率方法［15］；疊前逐層層速度反演速度建模方法［1617］；以疊前深度偏移為基礎(chǔ)進(jìn)行偏移速度和層速度分析，對(duì)速度模型的層位結(jié)構(gòu)及速度縱、橫向變化，采用三維可視化監(jiān)控和交互方式修改的方法［18］；共聚焦點(diǎn)層析速度建模方法［19］；多域雙重迭代法優(yōu)化更新層速度模型的方法［20］；聯(lián)合應(yīng)用基于層位的層析反演和基于網(wǎng)格的層析反演兩種速度建模方法［21］；VTI介質(zhì)多參數(shù)聯(lián)合走時(shí)層析成像的方法［22］等。目前應(yīng)用最廣、最為成熟的速度建模方法是網(wǎng)格層析速度建模方法［23］。各種速度建模方法都有其適應(yīng)性，由于YKL氣藏區(qū)三維地震勘探資料存在缺失近炮檢距資料等不利因素，單獨(dú)采用某一種速度建模方法難以解決該區(qū)塊三維資料的深度偏移準(zhǔn)確成像問題。

針對(duì)上述難點(diǎn)問題，綜合采用了井控［24］各向異性 網(wǎng)格層析 逐層反演多種速度建模技術(shù)，建立了合理的深度偏移速度模型，采用Kirchhoff方法進(jìn)行疊前深度偏移成像，最終在深度域?qū)⒛康膶臃謱訑?shù)據(jù)直接生成構(gòu)造圖。對(duì)3口井進(jìn)行驗(yàn)證，實(shí)鉆井深與地震預(yù)測(cè)深度誤差都在5m內(nèi)，較好地解決了YKL氣藏井震誤差大的問題。同時(shí)，通過重新處理，成像的品質(zhì)也得到明顯提高，為氣藏開發(fā)提供了可靠的基礎(chǔ)資料。

1 YKL氣藏井震誤差及問題分析

1．1 井震誤差分析

Kirchhoff疊前深度偏移層速度往往要高于測(cè)井的層速度，若將深度偏移成像深度直接進(jìn)行井深預(yù)測(cè)，絕對(duì)誤差會(huì)很大（見表1中的各向同性速度建模誤差）。在探區(qū)速度橫向變化很小的情況下，利用實(shí)際測(cè)井的層速度計(jì)算深度，將時(shí)間域成像數(shù)據(jù)通過井校變速成圖能夠得到較準(zhǔn)確的預(yù)測(cè)深度。YKL探區(qū)的地層產(chǎn)狀較為平坦，屬弱各向異性地區(qū)，但由于橫向速度存在一定的變化，時(shí)間域成像數(shù)據(jù)通過井校變速成圖存在較大誤差（見表1中的井校變速成圖誤差，19m以內(nèi)），尤其是井少的情況下，井間數(shù)據(jù)通過插值連接，難以得到精確的構(gòu)造形態(tài)。而采用井控各向異性網(wǎng)格層析逐層速度建模技術(shù)則不同，井間數(shù)據(jù)連接是靠數(shù)據(jù)本身驅(qū)動(dòng)的，在深度偏移成像條件約束下形成的構(gòu)造形態(tài)更為可靠，井震誤差?。ㄒ姳?中的井控VTI速度建模誤差，3米以內(nèi)）。

1．2 誤差原因分析

分析認(rèn)為，影響氣藏井震誤差的主要原因有以下幾個(gè)方面。

1）YKL氣藏探區(qū)三維地震資料采集參數(shù)為：8線8炮，接收線距200m，接收點(diǎn)距50m，炮線距300 m，炮點(diǎn)距100m，滿覆蓋次數(shù)24次，面元25m× 50m，最小偏移距600m，單邊接收。原始資料覆蓋次數(shù)太低，且缺失近偏移距數(shù)據(jù)，導(dǎo)致偏移后的道集在0～600ms基本沒有數(shù)據(jù)（圖1），不利于建立近地表速度模型；淺層0～2 500m的共成像點(diǎn)（CIP）道集中有效反射同相軸信息很少（圖2），淺層剩余延遲（RMO）無法準(zhǔn)確求取，嚴(yán)重影響速度建模精度。

2）淺層的測(cè)井曲線數(shù)據(jù)不可靠（圖3），標(biāo)定困難，無法準(zhǔn)確獲取淺層地質(zhì)分層，影響淺層速度估算的精度。

3）存在一定的速度各向異性問題，常規(guī)的各向同性偏移不能滿足該區(qū)的深度偏移需要，需采用各向異性參數(shù)建立速度模型。

表1 YKL氣藏實(shí)鉆井深及井震誤差統(tǒng)計(jì)結(jié)果 單位：m

圖1 CRP道集

圖2 CIP道集

圖3 測(cè)井曲線

2 井控各向異性網(wǎng)格層析速度建模

2．1 井控各向異性網(wǎng)格層析速度建模技術(shù)思路及流程

首先，利用疊前時(shí)間偏移得到的成像剖面及均方根速度，結(jié)合區(qū)內(nèi)地質(zhì)認(rèn)識(shí)，建立初始速度模型。在此基礎(chǔ)上進(jìn)行各向同性網(wǎng)格層析速度建模及疊前深度偏移，對(duì)偏移后CIP道集進(jìn)行剩余延遲分析，逐步拉平CIP道集，優(yōu)化速度模型，將每個(gè)深度偏移CIP道集的同相軸拉平。

其次，生成各向異性初始速度模型。選好區(qū)內(nèi)標(biāo)志層，利用區(qū)內(nèi)多口已知井信息求取各向異性參數(shù)，先消除大套穩(wěn)定地層速度對(duì)氣藏深度的影響，將第1個(gè)淺部穩(wěn)定標(biāo)志層的井震誤差消除。

然后，采用逐層建模技術(shù)，從淺到深在已知井的控制下，經(jīng)過多次更新各向異性參數(shù)和深度域?qū)铀俣饶Ｐ停饘幼龊弥猩畈繕?biāo)志層的偏移成像，直至各標(biāo)志層的井震誤差降低到規(guī)定范圍內(nèi)，確保最終YKL氣藏構(gòu)造成像準(zhǔn)確。

井控各向異性網(wǎng)格層析逐層速度建模技術(shù)流程如圖4所示。

2．2 初始速度模型建立

實(shí)際處理YKL地區(qū)資料，建立各向異性參數(shù)的初始模型包括以下3方面工作。

1）預(yù)處理，特別是做好靜校正處理。本區(qū)采用低速帶微測(cè)井成果約束下的折射層析靜校正及地表一致性剩余靜校正聯(lián)合應(yīng)用，避免靜校正不當(dāng)引起不利影響，生成可靠的疊前道集，再通過多次Kirchhoff疊前時(shí)間偏移均方根速度優(yōu)化迭代得到偏移成像數(shù)據(jù)體，然后對(duì)偏移體進(jìn)行井震標(biāo)定落實(shí)大套穩(wěn)定地層，并對(duì)其進(jìn)行構(gòu)造解釋。

2）利用解釋的大套層位按照空間插值方式約束均方根速度。通過DIX公式將均方根速度轉(zhuǎn)換為深度域?qū)铀俣?，再?duì)層速度進(jìn)行編輯、沿層外推和沿層平滑，使層速度與大套地層趨勢(shì)更吻合，得到初始深度域?qū)铀俣饶Ｐ?，進(jìn)行目標(biāo)線Kirchhoff疊前深度偏移。

3）從淺至深優(yōu)化初始速度模型。經(jīng)過多輪各向同性網(wǎng)格層析速度模型更新優(yōu)化及目標(biāo)線Kirchhoff疊前深度偏移迭代，確保主要地層CIP道集拉平以及偏移成像效果最佳，將最終迭代更新的各向同性網(wǎng)格層析速度模型V作為井控各向異性（VTI）網(wǎng)格層析的初始速度模型V0。另外，通過各向異性參數(shù)掃描，獲得ε與δ的比例關(guān)系，當(dāng)設(shè)定ε＝δ時(shí)，道集同相軸更平，模型速度與井速度擬合更好，層速度網(wǎng)格層析反演后，得到的構(gòu)造成像效果更好，且井震誤差小。各向異性參數(shù)的初始值設(shè)為：δ0＝0，ε0＝0。

2．3 井控各向異性網(wǎng)格層析速度建模概念

各向異性網(wǎng)格層析速度反演原理是從地下成像點(diǎn)向地表進(jìn)行射線追蹤，將之前在成像道集上拾取的剩余延遲沿著該射線路經(jīng)進(jìn)行反投影從而獲得各向異性參數(shù)和速度的更新量。需要注意的是，速度及各向異性參數(shù)是在規(guī)則的稀疏網(wǎng)格上更新的。各向異性網(wǎng)格層析方程為：

式中：Δtray表示深度偏移道集的延遲剩余時(shí)（RMOs）；flin表示各向異性參數(shù)的線性函數(shù)；Δv，Δδ，Δε表示解層析成像方程得到的各向異性中間參數(shù)剩余量［25］。

當(dāng)?shù)卣鸩ㄍㄟ^射線追蹤從地下構(gòu)造反射到地表時(shí)，時(shí)差是沿著追蹤射線的各向異性參數(shù)的一個(gè)線性函數(shù)，那么可以首先進(jìn)行射線追蹤，然后通過反演計(jì)算來更新各向異性參數(shù)，其中這些參數(shù)描述了以相速度波傳播的彈性參數(shù)，也就是湯姆森（Thomsen）各向異性介質(zhì)參數(shù)。在各向異性介質(zhì)中，Δtray／ΔWmisstie（旅行時(shí)／井震誤差）也可以認(rèn)為是各向異性參數(shù)誤差的線性函數(shù)［26］，因此，公式（1）也可表示為：

公式（2）定義為井控各向異性（VTI＋WellTie）網(wǎng)格層析速度迭代方程。井控各向異性參數(shù)計(jì)算公式為：

圖4 井控各向異性網(wǎng)格層析逐層速度建模技術(shù)流程

式中：ΔzI為由疊前深度偏移剖面獲得的各向同性層厚度；ΔzA為由測(cè)井時(shí)深關(guān)系獲得的各向異性層厚度；Va0為各向異性層速度；V0為各向異性初始層速度；δ為P波在垂直方向上的差異［27－30］。

2．4 井控各向異性參數(shù)模型修正

VTI各向異性介質(zhì)中的P波橫向速度比垂向速度快，測(cè)井速度是垂向分量，地震波速度包括垂向和橫向分量。在各向異性介質(zhì)中測(cè)井速度通常比地震各向同性的層速度小，如圖5所示。Thomsen參數(shù)ε代表P波水平方向和垂直方向的關(guān)系，其大小體現(xiàn)了由于各向異性而導(dǎo)致的同相軸彎曲程度。ε的求取通常依靠拾取剩余延遲譜，經(jīng)過井控層析校正后，剩余延遲譜基本歸于零值附近，中遠(yuǎn)偏移距處的同相軸得到拉平，從而進(jìn)一步提高偏移成像的效果。鑒于該區(qū)資料缺失近偏移距信息且覆蓋次數(shù)低，過少的CIP無法形成可靠的剩余延遲譜，為此，采取井控各向異性網(wǎng)格層析速度建模方法，借助已知井信息，在偏移成像數(shù)據(jù)體上逐層統(tǒng)計(jì)井震誤差，并根據(jù)公式（3）求取δ，且令ε＝δ，獲得初始ε，通過逐層多次迭代，直到誤差完全消除。

井控各向異性參數(shù)模型修正的步驟如下。

1）在時(shí)間域偏移成果數(shù)據(jù)上，通過井震標(biāo)定精確解釋主要標(biāo)志層（如：T42，T23，T43，其中T43是YKL氣藏目的層）。

2）將時(shí)間域T42轉(zhuǎn)到深度域，得到該層位的深度值，依據(jù)27口井信息求得的井震誤差，并計(jì)算Δδ0，然后將初始速度模型V0和Δδ0，Δε0共3個(gè)參數(shù)代入各向異性（VTI）網(wǎng)格層析迭代方程（2），求解該方程獲得速度更新量ΔV1，從而得到新的速度模型V1和新的偏移成像數(shù)據(jù)體。

4）按照步驟3）進(jìn)行網(wǎng)格層析迭代，經(jīng)過4～6次迭代，最后Δδn更新量趨于0，可認(rèn)為井震誤差基本消除，即消除了以上淺層速度對(duì)下覆地層的影響。同時(shí)和也隨著誤差的減小而減?。煌恚凑粘上襁^程地層的井震誤差經(jīng)過3～5次迭代可消除，最后對(duì)目的層進(jìn)行2～4次迭代消除井震誤差。

2．5 各向異性參數(shù)模型

經(jīng)過逐層井控，采用各向異性網(wǎng)格層析速度建模對(duì)速度進(jìn)行更新，最終使各向異性的層速度與實(shí)鉆井分層速度一致，求取的3層的各向異性參數(shù)δ分布在－0．12～0．07（圖6）。

圖5 YKL27口井實(shí)鉆層速度與各向同性偏移T42層速度對(duì)比

圖6 氣藏區(qū)求取的3層（）各向異性參數(shù)δ沿層平面分布

3 應(yīng)用效果分析

3．1 速度模型與氣藏分層的對(duì)比

YKL氣藏區(qū)不同速度建模方法得到的速度剖面結(jié)果如圖7所示。從圖7可以看出，不同速度建模方法得到的速度剖面存在明顯差異，井控各向異性網(wǎng)格層析速度建模方法得到的中間層速度變小了（圖7b）。

YKL氣藏區(qū)不同速度建模方法主要標(biāo)志層與層速度的匹配對(duì)比結(jié)果如圖8所示。從圖8可以看出，井控各向異性網(wǎng)格層析速度場(chǎng)分層界面與氣藏分層界面得到較好的匹配，說明速度得到了修正。

3．2 井震誤差情況

圖7 YKL氣藏區(qū)不同速度建模方法得到的速度剖面

圖8 YKL氣藏區(qū)不同速度建模方法主要標(biāo)志層與層速度的匹配對(duì)比

YKL氣藏區(qū)不同速度建模方法過部分井的疊前深度偏移剖面對(duì)比結(jié)果如圖9所示。從圖9可以看出，各向同性網(wǎng)格層析速度建模法存在明顯的深度誤差問題；應(yīng)用井控各向異性網(wǎng)格層析逐層速度建模后，成像精度得到明顯改善，氣藏構(gòu)造地震分層深度與測(cè)井實(shí)鉆分層深度基本閉合。

3．3 實(shí)際井的驗(yàn)證

圖9 YKL氣藏區(qū)不同速度建模方法過部分井的疊前深度偏移剖面

圖10 YKL氣藏區(qū)不同速度建模方法3口驗(yàn)證井的井震誤差情況對(duì)比

YKL氣藏區(qū)不同速度建模方法3口驗(yàn)證井的井震誤差情況對(duì)比結(jié)果如圖10所示。由圖10a可以看出，各向同性網(wǎng)格層析成像速度建模方法的井震誤差大；由圖10b可以看出，井控各向異性網(wǎng)格層析成像速度建模的井震誤差小，目的層氣藏深度預(yù)測(cè)值與實(shí)鉆井分層深度誤差在5m內(nèi)。

3．4 剖面品質(zhì)及構(gòu)造圖

YKL氣藏區(qū)新、老偏移剖面對(duì)比結(jié)果如圖11所示。從圖11可以看出，經(jīng)過井控各向異性網(wǎng)格層析速度建模和疊前深度偏移等技術(shù)的應(yīng)用，獲得的處理成像剖面品質(zhì)明顯提高，構(gòu)造成像精度更高，波組特征更清楚，斷裂清晰。

YKL氣藏區(qū)不同方法得到的構(gòu)造圖對(duì)比結(jié)果如圖12所示。從圖12可以看出，兩種方法得到的構(gòu)造形態(tài)有一定差異，經(jīng)已有鉆井資料證實(shí)白堊系構(gòu)造成像精度得到明顯提高，為YKL氣藏主體區(qū)井位加密調(diào)整以及外圍滾動(dòng)外擴(kuò)提供可靠數(shù)據(jù)支撐。

圖11 YKL氣藏區(qū)新、老偏移剖面

圖12 YKL氣藏區(qū)采用不同方法得到的構(gòu)造圖對(duì)比

4 結(jié)束語(yǔ)

井控各向異性網(wǎng)格層析速度建模和逐層成像速度建模方法技術(shù)的聯(lián)合應(yīng)用，是目前深度域提高速度建模精度和成像精度較好的方法，也是目前工業(yè)界較為普遍的做法。類似YKL氣藏區(qū)地震勘探覆蓋次數(shù)低，缺失近偏移距資料，近地表速度無法準(zhǔn)確獲取的難題，采用速度掃描、井控各向異性網(wǎng)格層析迭代及標(biāo)志反射層控制的手段解決淺層地層速度變化及各向異性等問題，是獲得高精度成像結(jié)果的關(guān)鍵，較好地解決了YKL氣藏井震誤差大的問題，得到的成果資料更好地滿足了氣藏開發(fā)對(duì)地震資料的高精度要求。在速度建模中重點(diǎn)考慮以下幾點(diǎn)：

1）處理解釋人員緊密結(jié)合，落實(shí)地震波組與測(cè)井分層的對(duì)應(yīng)關(guān)系，層位解釋、標(biāo)定必須準(zhǔn)確；

2）消除井震誤差首先需消除淺層大套穩(wěn)定地層速度對(duì)氣藏深度的影響，為后續(xù)深度成像迭代打好基礎(chǔ)；

3）對(duì)井震誤差要進(jìn)行合理性分析，做好井控約束外推的邊界處理，逐步降低井震誤差以提高成像精度；

4）盡可能充分利用多口已知井資料，提高約束的可靠性。

致謝：本研究得到了西北油田分公司勘探開發(fā)研究院及項(xiàng)目組等有關(guān)專家的大力支持，在此致以衷心的感謝。

［1］ 孫開峰，葉勇．泛克里格法在三維地震資料變速成圖中的應(yīng)用［J］．石油物探，1998，37（2）：112－117

SUN K F，YE Y．Applying the universal Kirging to the variable－velocity mapping of 3－D seismic data［J］．Geophysical Prospecting for Petroleum，1998，37（2）：112－117

［2］ 楊江峰，王咸彬，洪太遠(yuǎn)，等．準(zhǔn)噶爾盆地南緣山前帶構(gòu)造變速成圖［J］．石油物探，2008，47（2）：179－182

YANG J F，WANG X B，HONG T Y，et al．Application of varying velocity mapping technique for mountain front in sourthern Junggar basin［J］．Geophysical Prospecting for Petroleum，2008，47（2）：179－182

［3］ 葉勇，孫武亮，孫開峰．塔中圍斜區(qū)連片精細(xì)速度建模與變速成圖［J］．石油物探，2010，49（4）：401－406

YE Y，SUN W L，SUN K F．Detailed velocity modeling and velocity－variable structure mapping of eight－piecing processing in Tazhong area［J］．Geophysical Prospecting for Petroleum，2010，49（4）：401－406

［4］ 楊勤林，王彥春，張靜，等．濱里海盆地M區(qū)塊鹽下構(gòu)造變速成圖［J］．石油物探，2012，51（4）：377－382

YANG Q L，WANG Y C，ZHANG J，et al．Variablevelocity mapping for sub－salt structure in Block M，Pre－Caspian Basin［J］．Geophysical Prospecting for Petroleum，2012，51（4）：377－382

［5］ JONES I F．3－D prestack depth migration and velocity model building［J］．The Leading Edge，1998，19（7）：897－906

［6］ RESHEL M．The use of 3－D prestack depth imaging to estimate layer velocities and reflector position［J］．Geophysics，1997，62（1）：206－210

［7］ 李家康，樊佳芳，邊西燕，等．大規(guī)模并行處理機(jī)三維疊前深度偏移［J］．石油物探，1998，37（2）：77－87

LI J K，F(xiàn)AN J F，BIAN X Y，et al．3－D prestack depth migration on massive parallel processing system［J］．Geophysical Prospecting for Petroleum，1998，37（2）：77－87

［8］ 齊中山，王汝珍，敬朋貴，等．川東北地區(qū)膏鹽巖影響下的儲(chǔ)層成像技術(shù)與效果分析［J］．石油物探，2011，50（6）：607－611

QI Z S，WANG R Z，JING P G，et al．The imaging techniques of reservoirs influenced by saline deposit in the Northeastern Sichuan Basin［J］．Geophysical Prospecting for Petroleum，2011，50（6）：607－611

［9］ 劉西寧，劉司紅，馬秀國(guó)，等．疊前深度偏移技術(shù)在焉耆盆地寶南地區(qū)的應(yīng)用［J］．石油物探，2006，45（2）：197－201

LIU X N，LIU S H，MA X G，et al．The application and investigation of 3－D prestack depth migration in Baonan area of Yanqi basin［J］．Geophysical Prospecting for Petroleum，2006，45（2）：197－201

［10］ 郝守玲，趙群．橫向速度變化對(duì)構(gòu)造成像影響的物理模擬研究［J］．石油物探，2008，47（1）：49－54

HAO S L，ZHAO Q．Physical modeling research for the effect of lateral velocity on struc－tural imaging［J］．Geophysical Prospecting for Petroleum，2008，47（1）：49－54

［11］ 李鵬，劉志成，楊楠，等．玉北地區(qū)碳酸鹽巖縫洞型儲(chǔ)層精細(xì)成像方法應(yīng)用研究［J］．石油物探，2015，54（4）：443－451

LI P，LIU Z C，YANG N，et al．Application of fine target reservoir imaging in fractured－vuggy carbonate reservoir Yubei area［J］．Geophysical Prospecting for Petroleum，2015，54（4）：443－451

［12］ 王麗，屠世杰，趙傳雪，等．疊前深度偏移技術(shù)在江蘇火成巖地區(qū)的應(yīng)用［J］．石油物探，2006，45（1）：79－82

WANG L，TU S J，ZHAO C X，et al．The application of prestack depth migration in igneous rocks area of the Jiangsu Oilfield［J］．Geophysical Prospecting for Petroleum，2006，45（1）：79－82

［13］ 胡英，姚逢昌．用于疊前深度偏移速度建模的一種新方法［J］．石油勘探與開發(fā)，2000，27（2）：62－64

HU Y，YAO F C．A new approach to prestack depth migration model building［J］．Petroleum Exploration and Development，2000，27（2）：62－64

［14］ 方伍寶．三維疊前深度偏移的建模技術(shù)［J］．石油物探，2002，41（2）：132－135

FANG W B．The model building technology in 3－D prestack depth migration［J］．Geophysical Prospecting for Petroleum，2002，41（2）：132－135

［15］ 井西利，楊長(zhǎng)春，李幼銘．建立速度模型的層析成像方法研究［J］．石油物探，2002，41（1）：72－76

JING X L，YANG C C，LI Y M．Study of tomographic inversion in the construction of velocity model［J］．Geophysical Prospecting for Petroleum，2002，41（1）：72－76

［16］ 陳遵德，徐廣民，賀振華，等．三維逐層層速度反演方法及誤差分析［J］．石油物探，2000，39（4）：12－18

CHEN Z D，XU G M，HE Z H，et al．The 3Dinversion method of layer－by－layer interval velocity and erroranalysis［J］．Geophysical Prospecting for Petroleum，2000，39（4）：12－18

［17］ 王者順，樊佳芳，王衛(wèi)江．疊前逐層成像速度建模方法及其應(yīng)用［J］．石油物探，2004，43（4）：384－387

WANG Z S，F(xiàn)AN J F，WANG W J．Layer－by－layerprestack imaging velocity model building and application［J］．Geophysical Prospecting for Petroleum，2004，43（4）：384－387

［18］ 羅省賢，李錄明．三維疊前深度偏移速度模型建立方法［J］．石油物探，1999，38（4）：1－6

LUO X X，Li L M．Velocity model construction methods for 3－D prestack depth migration［J］．Geophysical Prospecting for Petroleum，1999，38（4）：1－6

［19］ 辛可鋒，王華忠，馬在田，等．共聚焦點(diǎn)層析速度建模方法［J］．石油物探，2005，44（4）：329－333

XIN K F，WANG H Z，MA Z T，et al．Tomographic velocity modeling based on CFP technology［J］．Geophysical Prospecting for Petroleum，2005，44（4）：329－333

［20］ 胡英，張研，陳立康，等．速度建模的影響因素與技術(shù)對(duì)策［J］．石油物探，2006，45（5）：503－507

HU Y，ZHANG Y，CHEN L K，et al．Influencing factors and technical strategies for velocity modeling［J］．Geophysical Prospecting for Petroleum，2006，45（5）：503－507

［21］ 馬彥彥，李國(guó)發(fā)，張星宇，等．疊前深度偏移速度建模方法分析［J］．石油地球物理勘探，2014，49（4）：687－693

MA Y Y，LI G F，ZHANG X Y，et al．Strategy of velocity model building in pre－stack depth migration［J］．Oil Geophysical Prospecting，2014，49（4）：687－693

［22］ 劉玉柱，王光銀，董良國(guó)，等．VTI介質(zhì)多參數(shù)聯(lián)合走時(shí)層析成像方法［J］．地球物理學(xué)報(bào)，2014，57（10）：3402－3410

LIU Y Z，WANG G Y，DONG L G，et al．Joint inversion of VTI parameters using nonlinear traveltime tomography［J］．Chinese Journal of Geophysics，2014，57（10）：3402－3410

［23］ 李慧，成德安，金婧．網(wǎng)格層析成像速度建模方法與應(yīng)用［J］．石油地球物理勘探，2013，48（增刊1）：12－16

LI H，CHENG D A，JIN J．Velocity model building based on grid tomography［J］．Oil Geophysical Prospecting，2013，48（S1）：12－16

［24］ 康有元，馬豐臣，滿紅霞，等．面向油田開發(fā)的井控地震資料處理技術(shù)及應(yīng)用［J］．石油物探，2013，52（4）：402－408

KANG Y Y，MA F C，MAN H X，et al．The well－controlling seismic data processing techniques faced oilfield development and their application［J］．Geophysical Prospecting for Petroleum，2013，52（4）：402－408

［25］ THOMSEN L．Weak elastic anisotropy［J］．Geophysics，1986，51（10）：1954－1966

［26］ 羅小明，凌云，牛濱華，等．井地聯(lián)合提取VTI介質(zhì)各向異性參數(shù)［J］．地學(xué)前緣，2005，12（4）：576－580

LUO X M，LING Y，NIU B H，et al．Estimation of anisotropy parameters in VTI media from surface P－wave seismic data and VSP［J］．Earth Science Frontiers，2005，12（4）：576－580

［27］ 鄧懷群，劉雯林，趙正茂．橫向各向同性介質(zhì)中縱波和轉(zhuǎn)換橫波的快速射線追蹤方法［J］．石油物探，2000，39（4）：1－11

DENG H Q，LIU W L，ZHAO Z M．Fast Ray－tracing Method for Compresional and Converted Waves in Transversely Isotropic Media［J］．Geophysical Prospecting for Petroleum，2000，39（4）：1－11

［28］ 徐亦鳴，黃中玉，劉路佳．各向異性介質(zhì)縱波速度分析［J］．石油物探，2004，43（5）：438－440

XU Y M，HUANG Z Y，LIU L J．P－wave velocity analysis in anistotropy medium［J］．Geophysical Prospecting for Petroleum，2004，43（5）：438－440

［29］ 周巍，郭全仕，劉旭躍，等．Thomsen參數(shù)對(duì)VTI介質(zhì)克?；舴虔B前時(shí)間深度偏移的影響［J］．地球物理學(xué)進(jìn)展，2014，29（6）：2866－2873

ZHOU W，GUO Q S，LIU X Y，et al．The influence of Thomsen parameters on Kirchhoff prestack depth migration in VTI media［J］．Progress in Geophysics，2014，29（6）：2866－2873

［30］ 李磊．橫向各向同性介質(zhì)Thomsen近似公式的適應(yīng)范圍［J］．石油物探，2008，47（2）：116－122

LI L．Applicability of Thomsen approximage formula in transversely isotropic media［J］．Geophysical Prospecting for Petroleum，2008，47（2）：116－122

（編輯：顧石慶）

The application of well－controlled anisotropy velocity modeling in YKL region

PEI Yunlong1，WANG Lixin1，WU Dali1，JIANG Bo1，DENG Feng2
（1．Sinopec Geophysical Research Institute，Nanjing211103，China；2．Sinopec Northwest Oilfield Branch Company，Urumqi 830011，China）

The urgent demand for high quality imaging data，especially high－precision geological structure，is a key problem in the development of YKL gas field．Because of missing near－offset data and low fold of seismic data，well－seismic misfit for structural depth prediction in this region cant meet the exploration requirement if we only use conventional velocity modeling．Therefore，well－controlled anisotropy grid－tomography technique was adopted to conduct velocity modeling layer－by－layer．Based on precise choosing the reasonable markers and building the initial velocity model in the selected area，we use known logging data to calculate anisotropic parameters within this region．Then we eliminate the influence of large set for stable formation velocity on gas reservoir depth，and from shallow to deep layer by layer inversion on marker layer，update iterative anisotropic parameters，velocity and depth migration imaging until the predicted depth errors from well log and seismic data for each marker layer to be reduced．In addition，it has been confirmed by 3drilling wells that the misfit between the predicted structural depths respectively from the final seismic migration results and logging depth measurements in this gas reservoir are all within 5m．So we can draw a conclusion that our technique can significantly improve the imaging accuracy for structure and fracture of gas reservoir so as to reduce the risk of production drilling．

PSDM，well－controlled，anisotropy，grid－tomo，velocity modeling，well－seismic depth error

P631

A

1000－1441（2017）03－0390－10

10．3969／j．issn．1000－1441．2017．03．009

裴云龍，王立歆，鄔達(dá)理，等．井控各向異性速度建模技術(shù)在YKL地區(qū)的應(yīng)用［J］．石油物探，2017，56（3）：390－399

PEI Yunlong，WANG Lixin，WU Dali，et al．The application of well－controlled anisotropy velocity modeling in YKL region［J］．Geophysical Prospecting for Petroleum，2017，56（3）：390－399

2016－02－18；改回日期：2016－10－26。

裴云龍（1984—），男，碩士，工程師，主要從事地震資料處理應(yīng)用研究。

國(guó)家科技重大專項(xiàng)“塔河油田6－7區(qū)室內(nèi)參數(shù)組合處理采集方案優(yōu)化研究”（2011ZX05049－001－002－004）資助。

This research is financially supported by the National Science and Technology Major Project of China（Grant No．2011ZX05049－001－002－004）．