一種大傾角地區(qū)的垂直地震數(shù)據(jù)橋式標(biāo)定方法
【專利摘要】本發(fā)明是油氣地震勘探的垂直地震橋式標(biāo)定方法，將NMO數(shù)據(jù)分離成兩部分，其一為數(shù)據(jù)格式SEGY的卷頭、道頭等說明部分，其二為純數(shù)值部分，對NMO數(shù)據(jù)純數(shù)值部分進(jìn)行NMO數(shù)據(jù)轉(zhuǎn)換，卷頭、道頭等說明信息，與轉(zhuǎn)換后矩陣數(shù)據(jù)拼合，得到SEGY格式的新NMO地震數(shù)據(jù)，對VSP橋式標(biāo)定。本發(fā)明與目前的VSP橋式標(biāo)定方法不同的是，標(biāo)定結(jié)果中不再需要VSP走廊疊加剖面，這是因?yàn)橥ㄟ^上述系列轉(zhuǎn)換，已將井口附近的地震反射軸轉(zhuǎn)移至NMO剖面一側(cè)，與地面地震剖面可直接進(jìn)行對比，在回避大傾角地區(qū)走廊不能同相疊加問題的同時，得到更直觀、更可靠的橋式標(biāo)定結(jié)果。
【專利說明】一種大傾角地區(qū)的垂直地震數(shù)據(jù)橋式標(biāo)定方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于非常規(guī)油氣地震勘探領(lǐng)域，是一種利用轉(zhuǎn)換正常動校正（NMO)數(shù)據(jù)進(jìn) 行大傾角地區(qū)的垂直地震橋式標(biāo)定方法。

【背景技術(shù)】
[0002] 橋式標(biāo)定是垂直地震（VSP)資料的重要用途之一，是利用VSP數(shù)據(jù)建立地面地震 反射與井中地層之間的關(guān)系，通常的做法是基于VSP時深關(guān)系將VSP的NMO剖面和走廊疊 加剖面、過井地震剖面、測井曲線、巖合柱狀圖組合起來，得到VSP橋式標(biāo)定成果。
[0003] 隨著近些年來VSP應(yīng)用范圍的不斷擴(kuò)大，在許多大傾角地區(qū)進(jìn)行了 VSP采集，現(xiàn)有 的VSP走廊標(biāo)定流程為：首先處理得到的NMO數(shù)據(jù)，然后在近井口位置沿某一時間（或深度） 進(jìn)行窄走廊切除，進(jìn)而利用數(shù)據(jù)疊加的方法得到走廊疊加剖面。由于地層傾角的存在，數(shù)據(jù) NMO波場中上行波反射同相軸與時間軸(或深度軸）存在夾角，因此利用現(xiàn)有的方法不能同 相疊加得到可靠的VSP走廊剖面，并且現(xiàn)有的橋式標(biāo)定方法主要是基于VSP走廊和過井地 震剖面的波阻特征進(jìn)行對比標(biāo)定，走廊不可靠就大大的影響了 VSP橋式標(biāo)定的準(zhǔn)確性。


【發(fā)明內(nèi)容】

[0004] 本發(fā)明的目的是提供一種標(biāo)定準(zhǔn)確、可靠性高的利用轉(zhuǎn)換NMO數(shù)據(jù)進(jìn)行大傾角地 區(qū)的垂直地震數(shù)據(jù)橋式標(biāo)定方法。
[0005] 本發(fā)明通過以下步驟實(shí)現(xiàn)：
[0006] 1)采集VSP數(shù)據(jù)，在VSP處理中，進(jìn)行TAR補(bǔ)償(真振幅恢復(fù)）和道均衡處理，得到 保真的轉(zhuǎn)換前NMO數(shù)據(jù)；
[0007] 將NMO數(shù)據(jù)分離成兩部分，其一為數(shù)據(jù)格式SEGY的卷頭、道頭等說明部分，其二為 純數(shù)值部分；
[0008] 所述的TAR補(bǔ)償和道均衡處理是保證剖面中不同位置的能量相對穩(wěn)定，不因轉(zhuǎn)換 處理出現(xiàn)能量不均衡現(xiàn)象。
[0009] 所述的處理在時間域和深度域中方法完全相同。
[0010] 所述的數(shù)據(jù)格式為SEGY格式。
[0011] 2)按照以下方式對NMO數(shù)據(jù)純數(shù)值部分進(jìn)行NMO數(shù)據(jù)轉(zhuǎn)換：
[0012] (1)將NMO純數(shù)值部分以矩陣的方式排列，矩陣的行數(shù)為最大采樣點(diǎn)數(shù)，列數(shù)為 NMO數(shù)據(jù)所包含的總道數(shù)；
[0013] (2)利用地震處理中初至拾取的方法，拾取NMO純數(shù)值部分?jǐn)?shù)據(jù)的原始初至?xí)r間 序列，每一道對應(yīng)一個初至?xí)r間；
[0014] (3)對初至?xí)r間序列進(jìn)行整數(shù)采樣計算，使各點(diǎn)時間值均為采樣率的整數(shù)倍，得到 新初至?xí)r間序列T' :
[0015] 所述的計算依下式進(jìn)行：
[0016] r, = int(會)* ddt U )
[0017] 其中T'為計算得到的新初至?xí)r間序列，T為步驟（2)拾取的原始初至?xí)r間序列， A t為數(shù)據(jù)的采樣率，int為取整函數(shù)；
[0018] 所述的采樣率為數(shù)據(jù)記錄時間除以最大采樣點(diǎn)數(shù)。
[0019] (4)依據(jù)新初至?xí)r間序列，將步驟（1)矩陣的每一道進(jìn)行分段，得到包含系列子矩 陣的新矩陣；
[0020] (5)將得到新矩陣進(jìn)行數(shù)據(jù)變換運(yùn)算；
[0021] 所述的數(shù)據(jù)變換運(yùn)算依下述方法進(jìn)行：
[0022] 首先取出矩陣的右下三角陣，然后對三角陣進(jìn)行左右翻轉(zhuǎn)，并將翻轉(zhuǎn)后三角陣中 各子矩陣順序不變的重新排布回原三角陣的子矩陣位置，進(jìn)而將新三角陣替換步驟（5)中 新矩陣的右下三角部分，最后展開新矩陣中各個子矩陣，得到與步驟（1)大小完全相同的轉(zhuǎn) 換后矩陣數(shù)據(jù)。
[0023] (6)將步驟1)中分離得到的卷頭、道頭等說明信息，與轉(zhuǎn)換后矩陣數(shù)據(jù)拼合，得到 SEGY格式的新NMO地震數(shù)據(jù)；
[0024] 3)按照以下方式進(jìn)行VSP橋式標(biāo)定：
[0025] 將過井地震剖面、測井?dāng)?shù)據(jù)、錄井剖面、以及步驟2)得到的轉(zhuǎn)換后矩陣數(shù)據(jù)以圖像 方式輸出，然后拼合，再以實(shí)鉆分層數(shù)據(jù)進(jìn)行VSP標(biāo)定。
[0026] 本發(fā)明與目前的VSP橋式標(biāo)定方法不同的是，標(biāo)定結(jié)果中不再需要VSP走廊疊加 剖面，這是因?yàn)橥ㄟ^上述系列轉(zhuǎn)換，已將井口附近的地震反射軸轉(zhuǎn)移至NMO剖面一側(cè)，與地 面地震剖面可直接進(jìn)行對比，在回避大傾角地區(qū)走廊不能同相疊加問題的同時，得到更直 觀、更可靠的橋式標(biāo)定結(jié)果。

【專利附圖】

【附圖說明】
[0027] 圖1原始NMO數(shù)據(jù)(左）和轉(zhuǎn)換后NMO數(shù)據(jù)(右）示意圖；
[0028] 圖2實(shí)際VSP資料的原始的NMO波場(左）和轉(zhuǎn)換后的NMO波場(右）；
[0029] 圖3原始的NMO波場的近井口切除剖面(左)、VSP走廊疊加剖面（中）、VSP走廊鑲 嵌到過井地震剖面(右）；
[0030] 圖4轉(zhuǎn)換后NMO波場與過井地震剖面拼接對比；
[0031] 圖5實(shí)驗(yàn)地區(qū)某井轉(zhuǎn)換NMO波場的VSP標(biāo)定實(shí)例。

【具體實(shí)施方式】
[0032] 本發(fā)明【具體實(shí)施方式】通過以下步驟實(shí)現(xiàn)：
[0033] 1)利用采集得到的VSP數(shù)據(jù)，進(jìn)行常規(guī)地震資料處理，其中TAR補(bǔ)償和道均衡處理 兩步流程影響到剖面中的能量分布，要通過對比觀察的方式確定最佳參數(shù)，進(jìn)而處理得到 保真的轉(zhuǎn)換前NMO數(shù)據(jù)；
[0034] 將NMO數(shù)據(jù)分離成兩部分，其一為數(shù)據(jù)格式SEGY的卷頭、道頭等說明部分，其二為 純數(shù)值部分。
[0035] 2)按照以下方式對NMO數(shù)據(jù)純數(shù)值部分進(jìn)行數(shù)據(jù)轉(zhuǎn)換：
[0036] (1)將上步得到的NMO純數(shù)值部分記作X，它由mXn個元素x構(gòu)成，其中x為數(shù)據(jù) 樣點(diǎn)值，雙下標(biāo)指示元素在矩陣中的位置，如X 11在矩陣中的位置是第1行第1列，其物理意 義是NMO數(shù)據(jù)中第1道的第1個樣點(diǎn)值，以此類推，n為NMO數(shù)據(jù)所包含的總道數(shù)，m為數(shù)據(jù) 最大采樣點(diǎn)數(shù)，如下式所示：
[0037]

【權(quán)利要求】
1. 一種大傾角地區(qū)的垂直地震數(shù)據(jù)橋式標(biāo)定方法，特點(diǎn)是通過以下步驟實(shí)現(xiàn)： 1) 采集VSP數(shù)據(jù)，在VSP處理中，進(jìn)行TAR補(bǔ)償和道均衡處理，得到保真的轉(zhuǎn)換前NMO 數(shù)據(jù)； 將NMO數(shù)據(jù)分離成兩部分，其一為數(shù)據(jù)格式SEGY的卷頭、道頭等說明部分，其二為純數(shù) 值部分； 2) 按照以下方式對NMO數(shù)據(jù)純數(shù)值部分進(jìn)行NMO數(shù)據(jù)轉(zhuǎn)換： (1) 將NMO純數(shù)值部分以矩陣的方式排列，矩陣的行數(shù)為最大采樣點(diǎn)數(shù)，列數(shù)為NMO數(shù) 據(jù)所包含的總道數(shù)； (2) 利用地震處理中初至拾取的方法，拾取NMO純數(shù)值部分?jǐn)?shù)據(jù)的原始初至?xí)r間序列， 每一道對應(yīng)一個初至?xí)r間； (3) 對初至?xí)r間序列進(jìn)行整數(shù)采樣計算，使各點(diǎn)時間值均為采樣率的整數(shù)倍，得到新初 至?xí)r間序列T': (4) 依據(jù)新初至?xí)r間序列T'，將步驟（1)矩陣的每一道進(jìn)行分段，得到包含系列子矩 陣的新矩陣； (5) 將得到新矩陣進(jìn)行數(shù)據(jù)變換運(yùn)算； (6) 將步驟1)中分離得到的卷頭、道頭等說明信息，與轉(zhuǎn)換后矩陣數(shù)據(jù)拼合，得到SEGY 格式的新NMO地震數(shù)據(jù)； 3) 按照以下方式進(jìn)行VSP橋式標(biāo)定： 將過井地震剖面、測井?dāng)?shù)據(jù)、錄井剖面、以及步驟2)得到的轉(zhuǎn)換后矩陣數(shù)據(jù)以圖像方式 輸出，然后拼合，再以實(shí)鉆分層數(shù)據(jù)進(jìn)行VSP標(biāo)定。
2. 根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟1)所述的TAR補(bǔ)償和道均衡處理是保證剖面中 不同位置的能量相對穩(wěn)定，不因轉(zhuǎn)換處理出現(xiàn)能量不均衡現(xiàn)象。
3. 根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟1)所述的道均衡處理在時間域和深度域中方法 完全相同。
4. 根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟1)所述的NMO數(shù)據(jù)格式為SEGY格式。
5. 根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟2)所述的整數(shù)采樣計算依下式進(jìn)行：
其中：T'為計算得到的新初至?xí)r間序列，T為步驟（2)拾取的原始初至?xí)r間序列，At 為數(shù)據(jù)的采樣率，int為取整函數(shù)； 所述的采樣率為數(shù)據(jù)記錄時間除以最大采樣點(diǎn)數(shù)。
6. 根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟2)所述的數(shù)據(jù)變換運(yùn)算依下述方法進(jìn)行： 首先取出矩陣的右下三角陣，然后對三角陣進(jìn)行左右翻轉(zhuǎn)，并將翻轉(zhuǎn)后三角陣中各子 矩陣順序不變的重新排布回原三角陣的子矩陣位置，進(jìn)而將新三角陣替換步驟（5)中新矩 陣的右下三角部分，最后展開新矩陣中各個子矩陣，得到與步驟（1)大小完全相同的轉(zhuǎn)換后 矩陣數(shù)據(jù)。
【文檔編號】G01V1/32GK104375174SQ201310356189
【公開日】2015年2月25日 申請日期:2013年8月15日 優(yōu)先權(quán)日:2013年8月15日
【發(fā)明者】蔡志東, 張慶紅, 彭繼新, 王沖 申請人:中國石油天然氣集團(tuán)公司, 中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司