一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法及裝置的制造方法
【專利摘要】本發(fā)明提供了一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法及裝置,涉及地震沉積學(xué)技術(shù)領(lǐng)域。方法包括:獲取第一疊后地震數(shù)據(jù)體,并對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后地震數(shù)據(jù)體;對所述第一零相位疊后地震數(shù)據(jù)體以一預(yù)設(shè)采樣率進(jìn)行重采樣,形成第二零相位疊后地震數(shù)據(jù)體;將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,形成寬頻帶地震數(shù)據(jù)體;對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第三零相位疊后地震數(shù)據(jù)體;通過有色反演法對所述第三零相位疊后地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù)據(jù)體。本發(fā)明可以解決現(xiàn)有技術(shù)中簡單地通過90°相位轉(zhuǎn)換技術(shù)不能達(dá)到區(qū)分巖性的功能,難以把地震數(shù)據(jù)體轉(zhuǎn)換為巖性數(shù)據(jù)體的問題。
【專利說明】
-種巖性數(shù)據(jù)體的轉(zhuǎn)換方法及裝置
技術(shù)領(lǐng)域
[0001] 本發(fā)明設(shè)及地震沉積學(xué)技術(shù)領(lǐng)域,特別設(shè)及地震沉積學(xué)中的一種巖性數(shù)據(jù)體的轉(zhuǎn) 換方法及裝置。
【背景技術(shù)】
[0002] 目前,地震沉積學(xué)是近年來發(fā)展起來的油氣精細(xì)勘探的新方法,它是在高精度= 維地震資料的基礎(chǔ)上,運(yùn)用精細(xì)沉積學(xué)模式,結(jié)合地球物理技術(shù)方法,對沉積體系分布特征 及其演化進(jìn)行精細(xì)研究的一口新興學(xué)科;它是地震地層學(xué)、層序地層學(xué)和高分辨率層序地 層學(xué)的繼承和發(fā)展。當(dāng)然,由于地震分辨率和研究手段的限制,地震沉積學(xué)還沒有形成一套 系統(tǒng)的理論體系,國內(nèi)近幾年雖然廣泛開展利用地震進(jìn)行沉積相、地層巖性識別的研究,但 還沒有出現(xiàn)有關(guān)地震沉積學(xué)的系統(tǒng)研究。因此,地震沉積學(xué)研究一方面顯示出了它的優(yōu)勢 和廣闊前景,另一方面相應(yīng)的技術(shù)手段又尚未成熟。
[0003] 地震巖性學(xué)和地震地貌學(xué)是地震沉積學(xué)的兩個重要分支學(xué)科。在現(xiàn)行的地震巖性 學(xué)中,巖性數(shù)據(jù)體的轉(zhuǎn)換方法采用90°相位轉(zhuǎn)換技術(shù),即將地震數(shù)據(jù)體轉(zhuǎn)化成巖性數(shù)據(jù)體, 從而使地震反射同相軸像測井曲線可W直接進(jìn)行巖性解釋。但是,由于陸相沉積盆地具有 沉積體系類型多、沉積相變快、砂體厚度變化大、成巖作用和成巖序列復(fù)雜等地質(zhì)難題,砂 泥巖阻抗疊置明顯,運(yùn)給利用地震速度資料預(yù)測沉積體系巖性、識別砂體特別是薄層砂體 帶來了極大困難。而簡單地通過90°相位轉(zhuǎn)換技術(shù)并不能達(dá)到區(qū)分巖性的功能,存在難W把 地震數(shù)據(jù)體轉(zhuǎn)換為巖性數(shù)據(jù)體的問題。
【發(fā)明內(nèi)容】
[0004] 本發(fā)明的實施例提供一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法及裝置,W解決現(xiàn)有技術(shù)中簡單 地通過90°相位轉(zhuǎn)換技術(shù)并不能達(dá)到區(qū)分巖性的功能,難W把地震數(shù)據(jù)體轉(zhuǎn)換為巖性數(shù)據(jù) 體的問題。
[0005] 為達(dá)到上述目的,本發(fā)明采用如下技術(shù)方案:
[0006] -種巖性數(shù)據(jù)體的轉(zhuǎn)換方法,包括:
[0007] 獲取第一疊后地震數(shù)據(jù)體,并對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正, 生成第一零相位疊后地震數(shù)據(jù)體;
[000引對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重采樣,形成第二零相位 疊后地震數(shù)據(jù)體;
[0009] 將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,形成寬頻帶地震數(shù)據(jù)體;
[0010] 對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第=零相位疊后地震數(shù)據(jù) 體;
[0011] 通過有色反演法對所述第=零相位疊后地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖 性數(shù)據(jù)體。
[0012] 具體的,對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后 地震數(shù)據(jù)體,包括:
[0013] 根據(jù)所述第一疊后地震數(shù)據(jù)體確定一第一預(yù)設(shè)時窗內(nèi)的第一原始地震道和第i個 相移后的第一地震道;
[0014] 報據(jù)公式:
[0015]
[0016] 確定第一原始地震道和第一地震道的第一相似系數(shù)Ri;其中,a(t)為第一原始地 震道;ai(t)為第i個相移后的第一地震道;A為所述第一預(yù)設(shè)時窗;
[0017] 根據(jù)所述第一相似系數(shù)Ri確定第一預(yù)設(shè)時窗內(nèi)第一疊后地震數(shù)據(jù)體相對于零相 位的左右波形的相位,并進(jìn)行相位旋轉(zhuǎn);
[0018] 當(dāng)所述第一相似系數(shù)Ri取最大值時,確定第一疊后地震數(shù)據(jù)體對應(yīng)的第一地震剩 余相位量;
[0019] 將所述第一地震剩余相位量進(jìn)行校正消除,生成第一零相位疊后地震數(shù)據(jù)體。
[0020] 具體的,對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第=零相位疊后地 震數(shù)據(jù)體,包括:
[0021] 根據(jù)所述寬頻帶地震數(shù)據(jù)體確定一第二預(yù)設(shè)時窗內(nèi)的第二原始地震道和第i個相 移后的第二地震道;
[0022] 巧據(jù)公式:
[0023]
[0024] 確定第二原始地震道和第二地震道的第二相似系數(shù)Ri';其中,a'(t)為第一原始 地震道;ai'(t)為第i個相移后的第二地震道;A'為所述第二預(yù)設(shè)時窗;
[0025] 根據(jù)所述第二相似系數(shù)Ri '確定第二預(yù)設(shè)時窗內(nèi)寬頻帶地震數(shù)據(jù)體相對于零相位 的左右波形的相位,并進(jìn)行相位旋轉(zhuǎn);
[0026] 當(dāng)所述第二相似系數(shù)Ri '取最大值時,確定寬頻帶地震數(shù)據(jù)體對應(yīng)的第二地震剩 余相位量;
[0027] 將所述第二地震剩余相位量進(jìn)行校正消除,生成第=零相位疊后地震數(shù)據(jù)體。 [00%]此外,所述預(yù)設(shè)采樣率為Ims。
[0029] 此外,該巖性數(shù)據(jù)體的轉(zhuǎn)換方法,還包括:
[0030] 對所述巖性數(shù)據(jù)體進(jìn)行Wheeler變換;
[0031] 或者,對所述巖性數(shù)據(jù)體進(jìn)行地層切片。
[0032] -種巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,包括:
[0033] 第一零相位疊后地震數(shù)據(jù)體生成單元,用于獲取第一疊后地震數(shù)據(jù)體,并對所述 第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后地震數(shù)據(jù)體;
[0034] 重采樣單元,用于對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重采 樣,形成第二零相位疊后地震數(shù)據(jù)體;
[0035] 反射系數(shù)反演單元,用于將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演, 形成寬頻帶地震數(shù)據(jù)體;
[0036] 第=零相位疊后地震數(shù)據(jù)體生成單元,用于對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相位估 算及校正,生成第=零相位疊后地震數(shù)據(jù)體;
[0037] 巖性數(shù)據(jù)體轉(zhuǎn)換單元,用于通過有色反演法對所述第=零相位疊后地震數(shù)據(jù)體進(jìn) 行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù)據(jù)體。
[0038] 此外,所述第一零相位疊后地震數(shù)據(jù)體生成單元,具體用于:
[0039] 根據(jù)所述第一疊后地震數(shù)據(jù)體確定一第一預(yù)設(shè)時窗內(nèi)的第一原始地震道和第i個 相移后的第一地震道;
[0040] 根據(jù)公式:
[0041]
[0042] 確定第一原始地震道和第一地震道的第一相似系數(shù)Ri;其中,a(t)為第一原始地 震道;ai(t)為第i個相移后的第一地震道;A為所述第一預(yù)設(shè)時窗;
[0043] 根據(jù)所述第一相似系數(shù)Ri確定第一預(yù)設(shè)時窗內(nèi)第一疊后地震數(shù)據(jù)體相對于零相 位的左右波形的相位,并進(jìn)行相位旋轉(zhuǎn);
[0044] 當(dāng)所述第一相似系數(shù)Ri取最大值時,確定第一疊后地震數(shù)據(jù)體對應(yīng)的第一地震剩 余相位量;
[0045] 將所述第一地震剩余相位量進(jìn)行校正消除,生成第一零相位疊后地震數(shù)據(jù)體。
[0046] 此外,所述第=零相位疊后地震數(shù)據(jù)體生成單元,具體用于:
[0047] 根據(jù)所述寬頻帶地震數(shù)據(jù)體確定一第二預(yù)設(shè)時窗內(nèi)的第二原始地震道和第i個相 移后的第二地震道;
[004引 巧據(jù)公古,
[0049
[0050] 確定第二原始地震道和第二地震道的第二相似系數(shù)Ri';其中,a'(t)為第一原始 地震道;ai'(t)為第i個相移后的第二地震道;A'為所述第二預(yù)設(shè)時窗;
[0051] 根據(jù)所述第二相似系數(shù)Ri '確定第二預(yù)設(shè)時窗內(nèi)寬頻帶地震數(shù)據(jù)體相對于零相位 的左右波形的相位,并進(jìn)行相位旋轉(zhuǎn);
[0052] 當(dāng)所述第二相似系數(shù)Ri '取最大值時,確定寬頻帶地震數(shù)據(jù)體對應(yīng)的第二地震剩 余相位量;
[0053] 將所述第二地震剩余相位量進(jìn)行校正消除,生成第=零相位疊后地震數(shù)據(jù)體。
[0054] 另外,所述重采樣單元中的預(yù)設(shè)采樣率為1ms。
[0055] 進(jìn)一步的,該巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,還包括:
[0056] Wheeler變換單元,用于對所述巖性數(shù)據(jù)體進(jìn)行Wheeler變換;
[0057] 地層切片單元,用于對所述巖性數(shù)據(jù)體進(jìn)行地層切片。
[0058] 本發(fā)明實施例提供的一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法及裝置,首先獲取第一疊后地震 數(shù)據(jù)體,并對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后地震數(shù) 據(jù)體;之后,對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重采樣,形成第二零相 位疊后地震數(shù)據(jù)體;進(jìn)而,將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,提高了第 二零相位疊后地震數(shù)據(jù)體的分辨率,形成寬頻帶地震數(shù)據(jù)體;對所述寬頻帶地震數(shù)據(jù)體進(jìn) 行相位估算及校正,生成第=零相位疊后地震數(shù)據(jù)體;通過有色反演法對所述第=零相位 疊后地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù)據(jù)體。本發(fā)明實施例能夠解決現(xiàn)有技術(shù) 中簡單地通過90°相位轉(zhuǎn)換技術(shù)并不能達(dá)到區(qū)分巖性的功能,難W把地震數(shù)據(jù)體轉(zhuǎn)換為巖 性數(shù)據(jù)體的問題。
【附圖說明】
[0059] 為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案,下面將對實施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡單地介紹,顯而易見地,下面描述中的附圖僅僅是本 發(fā)明的一些實施例,對于本領(lǐng)域普通技術(shù)人員來講,在不付出創(chuàng)造性勞動性的前提下,還可 W根據(jù)運(yùn)些附圖獲得其他的附圖。
[0060] 圖1為本發(fā)明實施例提供的一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法的流程圖一;
[0061] 圖2為本發(fā)明實施例提供的一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法的流程圖二;
[0062] 圖3為本發(fā)明實施例中的反射系數(shù)反演前后效果對比示意圖;
[0063] 圖4為本發(fā)明實施例中的地層切片與現(xiàn)有技術(shù)的90°相位轉(zhuǎn)換體切片的對比示意 圖;
[0064] 圖5為本發(fā)明實施例提供的一種巖性數(shù)據(jù)體的轉(zhuǎn)換裝置的結(jié)構(gòu)示意圖。
【具體實施方式】
[0065] 下面將結(jié)合本發(fā)明實施例中的附圖,對本發(fā)明實施例中的技術(shù)方案進(jìn)行清楚、完 整地描述,顯然,所描述的實施例僅僅是本發(fā)明一部分實施例,而不是全部的實施例?;?本發(fā)明中的實施例,本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他 實施例,都屬于本發(fā)明保護(hù)的范圍。
[0066] 如圖1所示,本發(fā)明實施例提供一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法,包括:
[0067] 步驟101、獲取第一疊后地震數(shù)據(jù)體,并對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算 及校正,生成第一零相位疊后地震數(shù)據(jù)體。
[0068] 步驟102、對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重采樣,形成第 二零相位疊后地震數(shù)據(jù)體。
[0069] 步驟103、將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,形成寬頻帶地震 數(shù)據(jù)體。
[0070] 步驟104、對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第=零相位疊后地 震數(shù)據(jù)體。
[0071] 步驟105、通過有色反演法對所述第=零相位疊后地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn) 換,生成巖性數(shù)據(jù)體。
[0072] 本發(fā)明實施例提供的一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法,首先獲取第一疊后地震數(shù)據(jù) 體,并對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后地震數(shù)據(jù)體; 之后,對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重采樣,形成第二零相位疊 后地震數(shù)據(jù)體;進(jìn)而,將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,提高了第二零 相位疊后地震數(shù)據(jù)體的分辨率,形成寬頻帶地震數(shù)據(jù)體;對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相 位估算及校正,生成第=零相位疊后地震數(shù)據(jù)體;通過有色反演法對所述第=零相位疊后 地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù)據(jù)體。本發(fā)明實施例能夠解決現(xiàn)有技術(shù)中簡 單地通過90°相位轉(zhuǎn)換技術(shù)并不能達(dá)到區(qū)分巖性的功能,難W把地震數(shù)據(jù)體轉(zhuǎn)換為巖性數(shù) 據(jù)體的問題。
[0073] 為了使本領(lǐng)域的技術(shù)人員更好的了解本發(fā)明,下面列舉一個更為詳細(xì)的實施例, 如圖2所示,本發(fā)明實施例提供一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法,包括:
[0074] 步驟201、獲取第一疊后地震數(shù)據(jù)體,根據(jù)所述第一疊后地震數(shù)據(jù)體確定一第一預(yù) 設(shè)時窗內(nèi)的第一原始地震道和第i個相移后的第一地震道。
[00巧]步驟202、根據(jù)公式:
[0076]
[0077] 確定第一原始地震道和第一地震道的第一相似系數(shù)Ri;其中,a(t)為第一原始地 震道;ai(t)為第i個相移后的第一地震道;A為所述第一預(yù)設(shè)時窗。
[0078] 步驟203、根據(jù)所述第一相似系數(shù)Ri確定第一預(yù)設(shè)時窗內(nèi)第一疊后地震數(shù)據(jù)體相 對于零相位的左右波形的相位,并進(jìn)行相位旋轉(zhuǎn)。
[00巧]此處,該相位旋轉(zhuǎn)的范圍為-90。至90。。
[0080]步驟204、當(dāng)所述第一相似系數(shù)Ri取最大值時,確定第一疊后地震數(shù)據(jù)體對應(yīng)的第 一地震剩余相位量。
[0081 ]步驟205、將所述第一地震剩余相位量進(jìn)行校正消除,生成第一零相位疊后地震數(shù) 據(jù)體。
[0082] 步驟206、對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重采樣,形成第 二零相位疊后地震數(shù)據(jù)體。
[0083] 此處,該預(yù)設(shè)采樣率可W為Ims。
[0084] 步驟207、將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,形成寬頻帶地震 數(shù)據(jù)體。
[0085] 此處,反射系數(shù)反演的方式可W參見現(xiàn)有技術(shù),如[UPorniaguine 0,Cas化gna J P.Spectral inversion:Lessons from modeling and Boonesville case study.SEG Technical Program Expanded Abstracts,2005,24:1638-1641. D或者,[2]張營,楊立英. 反射系數(shù)反演方法研究及其在薄層識別中的應(yīng)用.山東化工,2015,44(4) :95-97。此處步驟 207的目的為提高分辨率,恢復(fù)反射系數(shù),輸出高分辨率的寬頻帶地震數(shù)據(jù)體。
[0086] 此處,為了強(qiáng)調(diào)反射系數(shù)反演前后的效果,如圖3所示,其中圖3中的(a)為一原始 的地震數(shù)據(jù)體,(b)為反射系數(shù)反演之后的地震數(shù)據(jù)體,通過對比,通過反射系數(shù)反演獲得 的寬頻地震數(shù)據(jù)體明顯提高了分辨率,(b)中目的層內(nèi)圓圈處同相軸增多,而且整體面貌與 (a)差異不大。同時,高頻剖面上地層反射特征更加清楚,有利于層序的精細(xì)解釋。
[0087] 步驟208、根據(jù)所述寬頻帶地震數(shù)據(jù)體確定一第二預(yù)設(shè)時窗內(nèi)的第二原始地震道 和第i個相移后的第二地震道。
[0088] 步驟209、根據(jù)公式:
[0089]
[0090] 確定第二原始地震道和第二地震道的第二相似系數(shù)Ri';其中,a'(t)為第一原始 地震道;ai'(t)為第i個相移后的第二地震道;A'為所述第二預(yù)設(shè)時窗。
[0091] 步驟210、根據(jù)所述第二相似系數(shù)Ri'確定第二預(yù)設(shè)時窗內(nèi)寬頻帶地震數(shù)據(jù)體相對 于零相位的左右波形的相位,并進(jìn)行相位旋轉(zhuǎn)。
[0092] 此處,該相位旋轉(zhuǎn)的范圍為-90°至90°。
[0093] 步驟211、當(dāng)所述第二相似系數(shù)Ri '取最大值時,確定寬頻帶地震數(shù)據(jù)體對應(yīng)的第 二地震剩余相位量。
[0094] 步驟212、將所述第二地震剩余相位量進(jìn)行校正消除,生成第=零相位疊后地震數(shù) 據(jù)體。
[00M]步驟213、通過有色反演法對所述第=零相位疊后地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn) 換,生成巖性數(shù)據(jù)體。
[0096] 此處,具體的有色反演法可W參見文獻(xiàn):畢俊鳳,楊培杰.有色反演技術(shù)在少井區(qū) 巖性體預(yù)測中的應(yīng)用.物探與化探,2014,38(3): 558-565。此處不再寶述。通過有色反演法, 可W最大限度地保留原始的地震分辨率,不存在基于模型方法的多解性問題,能夠明顯地 反映巖相、巖性的空間變化。
[0097] 步驟214、對所述巖性數(shù)據(jù)體進(jìn)行Wheeler變換,或者對所述巖性數(shù)據(jù)體進(jìn)行地層 切片。
[0098] 例如,W通過巖性數(shù)據(jù)體進(jìn)行地層切片為例,如圖4所示,其中在圖4中的(a)為現(xiàn) 有技術(shù)的90°相位轉(zhuǎn)換體切片,而圖4中的(b)為本發(fā)明通過有色反演法得到的巖性數(shù)據(jù)體 進(jìn)行的地層切片。經(jīng)過對比,可W發(fā)現(xiàn),本發(fā)明的切片相對于現(xiàn)有技術(shù)的切片不損失地震橫 向分辨率,刻畫河道邊界更清晰,更適合用地層切片來表征地震相。
[0099] 本發(fā)明實施例提供的一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法,首先獲取第一疊后地震數(shù)據(jù) 體,并對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后地震數(shù)據(jù)體; 之后,對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重采樣,形成第二零相位疊 后地震數(shù)據(jù)體;進(jìn)而,將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,提高了第二零 相位疊后地震數(shù)據(jù)體的分辨率,形成寬頻帶地震數(shù)據(jù)體;對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相 位估算及校正,生成第=零相位疊后地震數(shù)據(jù)體;通過有色反演法對所述第=零相位疊后 地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù)據(jù)體。本發(fā)明實施例能夠解決現(xiàn)有技術(shù)中簡 單地通過90°相位轉(zhuǎn)換技術(shù)并不能達(dá)到區(qū)分巖性的功能,難W把地震數(shù)據(jù)體轉(zhuǎn)換為巖性數(shù) 據(jù)體的問題。
[0100] 對應(yīng)于上述圖1、圖2所述的方法實施例,如圖5所示,本發(fā)明實施例提供一種巖性 數(shù)據(jù)體的轉(zhuǎn)換裝置,包括:
[0101] 第一零相位疊后地震數(shù)據(jù)體生成單元301,用于獲取第一疊后地震數(shù)據(jù)體,并對所 述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后地震數(shù)據(jù)體。
[0102] 重采樣單元302,用于對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重 采樣,形成第二零相位疊后地震數(shù)據(jù)體。
[0103] 反射系數(shù)反演單元303,用于將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反 演,形成寬頻帶地震數(shù)據(jù)體。
[0104] 第=零相位疊后地震數(shù)據(jù)體生成單元304,用于對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相 位估算及校正,生成第=零相位疊后地震數(shù)據(jù)體。
[0105] 巖性數(shù)據(jù)體轉(zhuǎn)換單元305,用于通過有色反演法對所述第=零相位疊后地震數(shù)據(jù) 體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù)據(jù)體。
[0106] 此外,所述第一零相位疊后地震數(shù)據(jù)體生成單元301,具體可W :
[0107] 根據(jù)所述第一疊后地震數(shù)據(jù)體確定一第一預(yù)設(shè)時窗內(nèi)的第一原始地震道和第i個 相移后的第一地震道。
[010引 巧巧從才.
[0109]
[0110] 確定第一原始地震道和第一地震道的第一相似系數(shù)Ri;其中,a(t)為第一原始地 震道;ai(t)為第i個相移后的第一地震道;A為所述第一預(yù)設(shè)時窗。
[0111] 根據(jù)所述第一相似系數(shù)Ri確定第一預(yù)設(shè)時窗內(nèi)第一疊后地震數(shù)據(jù)體相對于零相 位的左右波形的相位,并進(jìn)行相位旋轉(zhuǎn)。
[0112] 當(dāng)所述第一相似系數(shù)Ri取最大值時,確定第一疊后地震數(shù)據(jù)體對應(yīng)的第一地震剩 余相位量。
[0113] 將所述第一地震剩余相位量進(jìn)行校正消除,生成第一零相位疊后地震數(shù)據(jù)體。
[0114] 此外,該第=零相位疊后地震數(shù)據(jù)體生成單元304,具體可W :
[0115] 根據(jù)所述寬頻帶地震數(shù)據(jù)體確定一第二預(yù)設(shè)時窗內(nèi)的第二原始地震道和第i個相 移后的第二地震道。
[0116] 根據(jù)公式:
[0117
[0118」佛足弟二原貓化晨迫和弟二化晨迫的弟二相似糸數(shù)Ri';其中,a'(t)為第一原始 地震道;ai'(t)為第i個相移后的第二地震道;A'為所述第二預(yù)設(shè)時窗。
[0119] 根據(jù)所述第二相似系數(shù)Ri '確定第二預(yù)設(shè)時窗內(nèi)寬頻帶地震數(shù)據(jù)體相對于零相位 的左右波形的相位,并進(jìn)行相位旋轉(zhuǎn)。
[0120] 當(dāng)所述第二相似系數(shù)Ri '取最大值時,確定寬頻帶地震數(shù)據(jù)體對應(yīng)的第二地震剩 余相位量。
[0121] 將所述第二地震剩余相位量進(jìn)行校正消除,生成第=零相位疊后地震數(shù)據(jù)體。
[0122] 另外,所述重采樣單元302,中的預(yù)設(shè)采樣率可W為1ms。
[0123] 進(jìn)一步的,如圖5所示,該巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,還可W包括:
[0124] Whee Ier變換單元306,可W對所述巖性數(shù)據(jù)體進(jìn)行Whee Ier變換。
[0125] 地層切片單元307,可W對所述巖性數(shù)據(jù)體進(jìn)行地層切片。
[0126] 值得說明的是,本發(fā)明實施例提供的一種巖性數(shù)據(jù)體的轉(zhuǎn)換裝置的具體實現(xiàn)方式 可W參見上述圖1和圖2對應(yīng)的方法實施例,此處不再寶述。
[0127] 本發(fā)明實施例提供的一種巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,通過該裝置,首先可W獲取第 一疊后地震數(shù)據(jù)體,并對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位 疊后地震數(shù)據(jù)體;之后,對所述第一零相位疊后地震數(shù)據(jù)體W-預(yù)設(shè)采樣率進(jìn)行重采樣,形 成第二零相位疊后地震數(shù)據(jù)體;進(jìn)而,將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反 演,提高了第二零相位疊后地震數(shù)據(jù)體的分辨率,形成寬頻帶地震數(shù)據(jù)體;對所述寬頻帶地 震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第=零相位疊后地震數(shù)據(jù)體;通過有色反演法對所述 第=零相位疊后地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù)據(jù)體。本發(fā)明實施例能夠解 決現(xiàn)有技術(shù)中簡單地通過90°相位轉(zhuǎn)換技術(shù)并不能達(dá)到區(qū)分巖性的功能,難W把地震數(shù)據(jù) 體轉(zhuǎn)換為巖性數(shù)據(jù)體的問題。
[0128] 本領(lǐng)域內(nèi)的技術(shù)人員應(yīng)明白,本發(fā)明的實施例可提供為方法、系統(tǒng)、或計算機(jī)程序 產(chǎn)品。因此,本發(fā)明可采用完全硬件實施例、完全軟件實施例、或結(jié)合軟件和硬件方面的實 施例的形式。而且,本發(fā)明可采用在一個或多個其中包含有計算機(jī)可用程序代碼的計算機(jī) 可用存儲介質(zhì)(包括但不限于磁盤存儲器、CD-ROM、光學(xué)存儲器等)上實施的計算機(jī)程序產(chǎn) 品的形式。
[0129] 本發(fā)明是參照根據(jù)本發(fā)明實施例的方法、設(shè)備(系統(tǒng))、和計算機(jī)程序產(chǎn)品的流程 圖和/或方框圖來描述的。應(yīng)理解可由計算機(jī)程序指令實現(xiàn)流程圖和/或方框圖中的每一流 程和/或方框、W及流程圖和/或方框圖中的流程和/或方框的結(jié)合??商峁┻\(yùn)些計算機(jī)程序 指令到通用計算機(jī)、專用計算機(jī)、嵌入式處理機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器W產(chǎn) 生一個機(jī)器,使得通過計算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實 現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的裝置。
[0130] 運(yùn)些計算機(jī)程序指令也可存儲在能引導(dǎo)計算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備W特 定方式工作的計算機(jī)可讀存儲器中,使得存儲在該計算機(jī)可讀存儲器中的指令產(chǎn)生包括指 令裝置的制造品,該指令裝置實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或 多個方框中指定的功能。
[0131] 運(yùn)些計算機(jī)程序指令也可裝載到計算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備上,使得在計 算機(jī)或其他可編程設(shè)備上執(zhí)行一系列操作步驟W產(chǎn)生計算機(jī)實現(xiàn)的處理,從而在計算機(jī)或 其他可編程設(shè)備上執(zhí)行的指令提供用于實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一 個方框或多個方框中指定的功能的步驟。
[0132] 本發(fā)明中應(yīng)用了具體實施例對本發(fā)明的原理及實施方式進(jìn)行了闡述,W上實施例 的說明只是用于幫助理解本發(fā)明的方法及其核屯、思想;同時,對于本領(lǐng)域的一般技術(shù)人員, 依據(jù)本發(fā)明的思想,在【具體實施方式】及應(yīng)用范圍上均會有改變之處,綜上所述,本說明書內(nèi) 容不應(yīng)理解為對本發(fā)明的限制。
【主權(quán)項】
1. 一種巖性數(shù)據(jù)體的轉(zhuǎn)換方法,其特征在于,包括: 獲取第一疊后地震數(shù)據(jù)體,并對所述第一疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成 第一零相位疊后地震數(shù)據(jù)體; 對所述第一零相位疊后地震數(shù)據(jù)體以一預(yù)設(shè)采樣率進(jìn)行重采樣,形成第二零相位疊后 地震數(shù)據(jù)體; 將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,形成寬頻帶地震數(shù)據(jù)體; 對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第三零相位疊后地震數(shù)據(jù)體; 通過有色反演法對所述第三零相位疊后地震數(shù)據(jù)體進(jìn)行巖性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù) 據(jù)體。2. 根據(jù)權(quán)利要求1所述的巖性數(shù)據(jù)體的轉(zhuǎn)換方法,其特征在于,對所述第一疊后地震數(shù) 據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后地震數(shù)據(jù)體,包括: 根據(jù)所述第一疊后地震數(shù)據(jù)體確定一第一預(yù)設(shè)時窗內(nèi)的第一原始地震道和第i個相移 后的第一地震道; 根據(jù)公式:確定第一原始地震道和第一地震道的第一相似系數(shù)R1;其中,a(t)為第一原始地震道; aKt)為第i個相移后的第一地震道;A為所述第一預(yù)設(shè)時窗; 根據(jù)所述第一相似系數(shù)心確定第一預(yù)設(shè)時窗內(nèi)第一疊后地震數(shù)據(jù)體相對于零相位的左 右波形的相位,并進(jìn)行相位旋轉(zhuǎn); 當(dāng)所述第一相似系數(shù)心取最大值時,確定第一疊后地震數(shù)據(jù)體對應(yīng)的第一地震剩余相 位量; 將所述第一地震剩余相位量進(jìn)行校正消除,生成第一零相位疊后地震數(shù)據(jù)體。3. 根據(jù)權(quán)利要求1所述的巖性數(shù)據(jù)體的轉(zhuǎn)換方法,其特征在于,對所述寬頻帶地震數(shù)據(jù) 體進(jìn)行相位估算及校正,生成第三零相位疊后地震數(shù)據(jù)體,包括: 根據(jù)所述寬頻帶地震數(shù)據(jù)體確定一第二預(yù)設(shè)時窗內(nèi)的第二原始地震道和第i個相移后 的第二地震道; 根據(jù)公式:確定第二原始地震道和第二地震道的第二相似系數(shù)心';其中,a'(t)為第一原始地震 道為第i個相移后的第二地震道;A'為所述第二預(yù)設(shè)時窗; 根據(jù)所述第二相似系數(shù)心'確定第二預(yù)設(shè)時窗內(nèi)寬頻帶地震數(shù)據(jù)體相對于零相位的左 右波形的相位,并進(jìn)行相位旋轉(zhuǎn); 當(dāng)所述第二相似系數(shù)& '取最大值時,確定寬頻帶地震數(shù)據(jù)體對應(yīng)的第二地震剩余相位 量; 將所述第二地震剩余相位量進(jìn)行校正消除,生成第三零相位疊后地震數(shù)據(jù)體。4. 根據(jù)權(quán)利要求1所述的巖性數(shù)據(jù)體的轉(zhuǎn)換方法,其特征在于,所述預(yù)設(shè)采樣率為lms。5. 根據(jù)權(quán)利要求1所述的巖性數(shù)據(jù)體的轉(zhuǎn)換方法,其特征在于,還包括: 對所述巖性數(shù)據(jù)體進(jìn)行Wheeler變換; 或者,對所述巖性數(shù)據(jù)體進(jìn)行地層切片。6. -種巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,其特征在于,包括: 第一零相位疊后地震數(shù)據(jù)體生成單元,用于獲取第一疊后地震數(shù)據(jù)體,并對所述第一 疊后地震數(shù)據(jù)體進(jìn)行相位估算及校正,生成第一零相位疊后地震數(shù)據(jù)體; 重采樣單元,用于對所述第一零相位疊后地震數(shù)據(jù)體以一預(yù)設(shè)采樣率進(jìn)行重采樣,形 成第二零相位疊后地震數(shù)據(jù)體; 反射系數(shù)反演單元,用于將所述第二零相位疊后地震數(shù)據(jù)體進(jìn)行反射系數(shù)反演,形成 寬頻帶地震數(shù)據(jù)體; 第三零相位疊后地震數(shù)據(jù)體生成單元,用于對所述寬頻帶地震數(shù)據(jù)體進(jìn)行相位估算及 校正,生成第三零相位疊后地震數(shù)據(jù)體; 巖性數(shù)據(jù)體轉(zhuǎn)換單元,用于通過有色反演法對所述第三零相位疊后地震數(shù)據(jù)體進(jìn)行巖 性數(shù)據(jù)體轉(zhuǎn)換,生成巖性數(shù)據(jù)體。7. 根據(jù)權(quán)利要求6所述的巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,其特征在于,所述第一零相位疊后地 震數(shù)據(jù)體生成單元,具體用于: 根據(jù)所述第一疊后地震數(shù)據(jù)體確定一第一預(yù)設(shè)時窗內(nèi)的第一原始地震道和第i個相移 后的第一地震道; 根據(jù)公式:確定第一原始地震道和第一地震道的第一相似系數(shù)R1;其中,a(t)為第一原始地震道; aKt)為第i個相移后的第一地震道;A為所述第一預(yù)設(shè)時窗; 根據(jù)所述第一相似系數(shù)心確定第一預(yù)設(shè)時窗內(nèi)第一疊后地震數(shù)據(jù)體相對于零相位的左 右波形的相位,并進(jìn)行相位旋轉(zhuǎn); 當(dāng)所述第一相似系數(shù)心取最大值時,確定第一疊后地震數(shù)據(jù)體對應(yīng)的第一地震剩余相 位量; 將所述第一地震剩余相位量進(jìn)行校正消除,生成第一零相位疊后地震數(shù)據(jù)體。8. 根據(jù)權(quán)利要求6所述的巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,其特征在于,所述第三零相位疊后地 震數(shù)據(jù)體生成單元,具體用于: 根據(jù)所述寬頻帶地震數(shù)據(jù)體確定一第二預(yù)設(shè)時窗內(nèi)的第二原始地震道和第i個相移后 的第二地震道; 根據(jù)公式:確定第二原始地震道和第二地震道的第二相似系數(shù)心';其中,a'(t)為第一原始地震 道為第i個相移后的第二地震道;A'為所述第二預(yù)設(shè)時窗; 根據(jù)所述第二相似系數(shù)心'確定第二預(yù)設(shè)時窗內(nèi)寬頻帶地震數(shù)據(jù)體相對于零相位的左 右波形的相位,并進(jìn)行相位旋轉(zhuǎn); 當(dāng)所述第二相似系數(shù)& '取最大值時,確定寬頻帶地震數(shù)據(jù)體對應(yīng)的第二地震剩余相位 量; 將所述第二地震剩余相位量進(jìn)行校正消除,生成第三零相位疊后地震數(shù)據(jù)體。9. 根據(jù)權(quán)利要求6所述的巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,其特征在于,所述重采樣單元中的預(yù) 設(shè)采樣率為lms。10. 根據(jù)權(quán)利要求6所述的巖性數(shù)據(jù)體的轉(zhuǎn)換裝置,其特征在于,還包括: Whee 1 er變換單元,用于對所述巖性數(shù)據(jù)體進(jìn)行Whee 1 er變換; 地層切片單元,用于對所述巖性數(shù)據(jù)體進(jìn)行地層切片。
【文檔編號】G01V1/32GK105954802SQ201610444652
【公開日】2016年9月21日
【申請日】2016年6月20日
【發(fā)明人】竇玉壇, 蘇明軍, 劉化清
【申請人】中國石油天然氣股份有限公司