用于識(shí)別斷層的相干增強(qiáng)處理方法
【專利摘要】一種用于識(shí)別斷層的相干增強(qiáng)處理方法，通過(guò)確定目的層段在地震數(shù)據(jù)體上的雙程反射時(shí)間范圍；對(duì)目的層段對(duì)應(yīng)的原始地震資料進(jìn)行構(gòu)造導(dǎo)向?yàn)V波，再進(jìn)行相干運(yùn)算，得到構(gòu)造導(dǎo)向?yàn)V波本征值相干數(shù)據(jù)體；對(duì)數(shù)據(jù)體進(jìn)行反射強(qiáng)度交流分量的運(yùn)算，得到交流分量數(shù)據(jù)體，將數(shù)據(jù)體正負(fù)對(duì)調(diào)，得到正負(fù)值對(duì)調(diào)的相干反射強(qiáng)度交流分量數(shù)據(jù)體；將數(shù)據(jù)體進(jìn)行數(shù)據(jù)體充零運(yùn)算，將表現(xiàn)為斷層的低相干的正值全部保留，得到正值數(shù)據(jù)體；將相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體進(jìn)行螞蟻?zhàn)粉欉\(yùn)算，得到螞蟻?zhàn)粉檾?shù)據(jù)體；將螞蟻?zhàn)粉檾?shù)據(jù)體與原始地震進(jìn)行融合，生成融合數(shù)據(jù)體；利用在原始地震數(shù)據(jù)體上解釋的層位數(shù)據(jù)在融合數(shù)據(jù)體上提取沿層切片，得到斷層與原始地震平面圖。
【專利說(shuō)明】
用于識(shí)別斷層的相干增強(qiáng)處理方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地球物理勘探方法，是一種用于識(shí)別斷層的相干增強(qiáng)的處理方法。

【背景技術(shù)】
[0002]相干是識(shí)別斷層有效的一種地球物理方法，目前的相干有道相關(guān)相干、方差、相似性相干、本征值相干等相干方法。但無(wú)論哪種相干方法，在地震資料品質(zhì)差的情況下，都不能獲得理想的效果。為了突出相干效果，目前流行一些增強(qiáng)相干的方法，包括在相干運(yùn)算之前對(duì)地震資料進(jìn)行疊后處理，包括優(yōu)勢(shì)頻帶濾波、構(gòu)造導(dǎo)向?yàn)V波等，雖然相干得到了增強(qiáng)，但有時(shí)效果仍不理想。目前還有的方法是在相干之后進(jìn)行數(shù)據(jù)處理，包括對(duì)相干數(shù)據(jù)體進(jìn)行螞蟻?zhàn)粉櫟?，但螞蟻?zhàn)粉櫴苤朴谙喔尚Ч挠绊?，地層因素影響較重，所以螞蟻?zhàn)粉櫟慕Y(jié)果看起來(lái)過(guò)于破碎，除了斷層因素外地層因素也包含在內(nèi)，很難區(qū)分哪些是斷層，哪些是地層。


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

[0003]為解決上述問(wèn)題，本發(fā)明目的在于提供一種用于識(shí)別斷層的相干增強(qiáng)的處理方法，以此，有效解決上述地震資料品質(zhì)較差的問(wèn)題。
[0004]本發(fā)明所提供的用于識(shí)別斷層的相干增強(qiáng)處理方法，具體包含如下步驟:步驟一:通過(guò)鉆井分層、聲波時(shí)差測(cè)井曲線、密度測(cè)井曲線制作合成地震記錄，并進(jìn)行井震層位標(biāo)定，通過(guò)研究區(qū)內(nèi)已知井資料建立井震關(guān)系，確定目的層段在地震數(shù)據(jù)體上的雙程反射時(shí)間范圍；步驟二:將目的層段對(duì)應(yīng)的原始地震資料置于所述雙程反射時(shí)間范圍內(nèi)，以所述雙程反射時(shí)間范圍為時(shí)窗對(duì)目的層段對(duì)應(yīng)的原始地震資料進(jìn)行構(gòu)造導(dǎo)向?yàn)V波，得到導(dǎo)向?yàn)V波地震資料數(shù)據(jù)體一；步驟三:對(duì)導(dǎo)向?yàn)V波地震資料數(shù)據(jù)體一進(jìn)行相干運(yùn)算，得到構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體二 ；步驟四:對(duì)構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體二進(jìn)行反射強(qiáng)度交流分量的運(yùn)算，得到相干體反射強(qiáng)度交流分量數(shù)據(jù)體三，將相干體反射強(qiáng)度交流分量數(shù)據(jù)體三正負(fù)對(duì)調(diào)，得到正負(fù)值對(duì)調(diào)的相干反射強(qiáng)度交流分量數(shù)據(jù)體四；步驟五:將正負(fù)值對(duì)調(diào)的相干反射強(qiáng)度交流分量數(shù)據(jù)體四進(jìn)行數(shù)據(jù)體賦值運(yùn)算，將表現(xiàn)為斷層的低相干的正值全部保留，得到相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體五；步驟六:將相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體五進(jìn)行螞蟻?zhàn)粉欉\(yùn)算，得到相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六；步驟七:將相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六與原始地震進(jìn)行融合，生成融合數(shù)據(jù)體七；步驟八:利用在原始地震數(shù)據(jù)體上解釋的層位數(shù)據(jù)在融合數(shù)據(jù)體七上提取沿層切片，得到斷層與原始地震平面圖。
[0005]在上述實(shí)施例中，所述導(dǎo)向?yàn)V波的算法優(yōu)選的還包含:中值濾波或均值濾波或主分量濾波。
[0006]在上述實(shí)施例中，所述步驟三中的相干算法優(yōu)選的還包含:本征值相干或相關(guān)或相似性相干。
[0007]在上述實(shí)施例中，優(yōu)選的所述相干算法還包含:通過(guò)5點(diǎn)相干或7點(diǎn)相干或9點(diǎn)相干，得到構(gòu)造導(dǎo)向?yàn)V波本征值相干數(shù)據(jù)體二。
[0008]在上述實(shí)施例中，所述步驟四中優(yōu)選的還包含:數(shù)據(jù)體三在數(shù)值上是正值和負(fù)值相間的數(shù)據(jù)體，低相干的斷層表現(xiàn)為負(fù)值，高相干的地層信息表現(xiàn)為正值。
[0009]在上述實(shí)施例中，所述數(shù)據(jù)體四中優(yōu)選的還包含:斷層表現(xiàn)為正值，地層信息表現(xiàn)為負(fù)值。
[0010]在上述實(shí)施例中，所述步驟六中的蟻群算法優(yōu)選的還包含:非侵略性蟻群算法或侵略性算法，得到相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六。
[0011]在上述實(shí)施例中，所述步驟七中優(yōu)選的還包含:融合時(shí)螞蟻?zhàn)粉檾?shù)據(jù)體范圍選擇為 _0.7-1 ο
[0012]在上述實(shí)施例中，優(yōu)選的還包含:所述原始地震數(shù)據(jù)體選擇為8位數(shù)據(jù)體，生成融合數(shù)據(jù)體七。
[0013]本發(fā)明的有益技術(shù)效果在于:在地震資料品質(zhì)較差的情況下，通過(guò)本發(fā)明能夠得到更好的識(shí)別斷層的數(shù)據(jù)體，與實(shí)際鉆井情況吻合較好。

【專利附圖】

【附圖說(shuō)明】
[0014]此處所說(shuō)明的附圖用來(lái)提供對(duì)本發(fā)明的進(jìn)一步理解，構(gòu)成本申請(qǐng)的一部分，并不構(gòu)成對(duì)本發(fā)明的限定。在附圖中:
[0015]圖1為利用鉆井分層、聲波測(cè)井曲線、密度測(cè)井曲線、制作的合成地震記錄和層位標(biāo)定；
[0016]圖2為構(gòu)造導(dǎo)向?yàn)V波前后地震剖面的對(duì)比圖，左圖為原始地震剖面，右圖為構(gòu)造導(dǎo)向?yàn)V波剖面；
[0017]圖3為構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體的時(shí)間切片；
[0018]圖4為構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體的反射強(qiáng)度交流分量數(shù)據(jù)體的時(shí)間切片；
[0019]圖5為將正負(fù)對(duì)調(diào)的構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體的反射強(qiáng)度交流分量數(shù)據(jù)體負(fù)數(shù)賦值為O的數(shù)據(jù)體的時(shí)間切片；
[0020]圖6為將圖4負(fù)數(shù)賦值為O的數(shù)據(jù)體進(jìn)行螞蟻?zhàn)粉櫤蟮臄?shù)據(jù)體的時(shí)間切片；
[0021]圖7螞蟻?zhàn)粉檾?shù)據(jù)體與原始數(shù)據(jù)體融合的數(shù)據(jù)體的時(shí)間切片。

【具體實(shí)施方式】
[0022]為使本發(fā)明實(shí)施例的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下面結(jié)合實(shí)施例和附圖，對(duì)本發(fā)明做進(jìn)一步詳細(xì)說(shuō)明。在此，本發(fā)明的示意性實(shí)施例及其說(shuō)明用于解釋本發(fā)明，但并不作為對(duì)本發(fā)明的限定。
[0023]請(qǐng)參考圖1至圖7所示，本發(fā)明所提供的用于識(shí)別斷層的相干增強(qiáng)的處理方法具體如下:
[0024]步驟一:通過(guò)鉆井分層、聲波時(shí)差測(cè)井曲線、密度測(cè)井曲線制作合成地震記錄，并進(jìn)行井震層位標(biāo)定，通過(guò)研究區(qū)內(nèi)已知井資料建立井震關(guān)系，確定目的層段在地震數(shù)據(jù)體上的雙程反射時(shí)間范圍；例如:通過(guò)鉆井分層、聲波時(shí)差測(cè)井曲線、密度測(cè)井曲線制作合成地震記錄，并進(jìn)行井震層位標(biāo)定，通過(guò)研究區(qū)內(nèi)已知井資料建立井震關(guān)系，確定目的層段在地震數(shù)據(jù)體上的雙程反射時(shí)間范圍，根據(jù)圖1所示的標(biāo)定結(jié)果，目的層段的雙程反射時(shí)間確定為2400-2900毫秒。
[0025]步驟二:將目的層段對(duì)應(yīng)的原始地震資料置于所述雙程反射時(shí)間范圍內(nèi)，以所述雙程反射時(shí)間范圍為時(shí)窗對(duì)目的層段對(duì)應(yīng)的原始地震資料進(jìn)行構(gòu)造導(dǎo)向?yàn)V波，得到導(dǎo)向?yàn)V波地震資料數(shù)據(jù)體一；例如:將目的層段對(duì)應(yīng)的原始地震資料置于2400-2900ms內(nèi)，以該時(shí)間范圍為時(shí)窗對(duì)原始地震資料進(jìn)行構(gòu)造導(dǎo)向?yàn)V波，為了取得最佳效果，導(dǎo)向?yàn)V波的算法可采用主分量濾波，得到導(dǎo)向?yàn)V波后地震資料數(shù)據(jù)體一，圖2所示為構(gòu)造導(dǎo)向?yàn)V波后的地震數(shù)據(jù)體一的地震剖面。
[0026]在上述實(shí)施例中，所述導(dǎo)向?yàn)V波的算法采用中值濾波或均值濾波或主分量濾波均可得到數(shù)據(jù)體一；其中主分量濾波效果較好，但本發(fā)明并不限制該導(dǎo)向?yàn)V波僅僅為主分量濾波。
[0027]步驟三:對(duì)導(dǎo)向?yàn)V波地震資料數(shù)據(jù)體一進(jìn)行相干運(yùn)算，得到構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體二 ；例如:對(duì)導(dǎo)向?yàn)V波地震資料數(shù)據(jù)體一進(jìn)行相干運(yùn)算，為了得到最佳效果，相干算法可采用本征值相干，為了獲得較快的運(yùn)算速度，本征值相干選擇5點(diǎn)相干，得到構(gòu)造導(dǎo)向?yàn)V波本征值相干數(shù)據(jù)體二，圖3所示為得到的數(shù)據(jù)體二的2800ms時(shí)間切片，由于地震資料品質(zhì)差的原因，此時(shí)效果比較雜亂，斷層不清晰。
[0028]在上述實(shí)施例中，所述步驟三中的相干算法采用本征值相干、相關(guān)、相似性相干都可得到需要的數(shù)據(jù)體二 ；其中本征值相干算法效果較好，但本發(fā)明并不限制該相干算法僅僅為本征值相干。
[0029]該本征值相關(guān)算法采用5點(diǎn)相干或7點(diǎn)相干或9點(diǎn)相干，得到構(gòu)造導(dǎo)向?yàn)V波本征值相干數(shù)據(jù)體二 ；其中5點(diǎn)想干算法獲得的數(shù)據(jù)體二效果較好。
[0030]步驟四:對(duì)構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體二進(jìn)行反射強(qiáng)度交流分量的運(yùn)算，得到相干體反射強(qiáng)度交流分量數(shù)據(jù)體三，將相干體反射強(qiáng)度交流分量數(shù)據(jù)體三正負(fù)對(duì)調(diào)，得到正負(fù)值對(duì)調(diào)的相干反射強(qiáng)度交流分量數(shù)據(jù)體四；例如:對(duì)構(gòu)造導(dǎo)向?yàn)V波本征值相干數(shù)據(jù)體二進(jìn)行反射強(qiáng)度交流分量的運(yùn)算，得到相干體反射強(qiáng)度交流分量數(shù)據(jù)體三，數(shù)據(jù)體三在數(shù)值上是正值和負(fù)值相間的數(shù)據(jù)體，低相干的斷層表現(xiàn)為負(fù)值，高相干的地層信息表現(xiàn)為正值；將數(shù)據(jù)體三乘以-1得到正負(fù)符號(hào)對(duì)調(diào)的相干體反射強(qiáng)度交流分量數(shù)據(jù)體四，數(shù)據(jù)體中低相干的斷層表現(xiàn)為正值，高相干的地層信息表現(xiàn)為負(fù)值。圖4為數(shù)據(jù)體四的2800ms時(shí)間切片。
[0031]步驟五:將正負(fù)值對(duì)調(diào)的相干反射強(qiáng)度交流分量數(shù)據(jù)體四進(jìn)行數(shù)據(jù)體賦值運(yùn)算，將表現(xiàn)為斷層的低相干的正值全部保留，得到相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體五；例如:將相干體反射強(qiáng)度交流分量數(shù)據(jù)體四進(jìn)行數(shù)據(jù)體賦值運(yùn)算，將表現(xiàn)為地層因素的高相干的負(fù)數(shù)值全部賦值為0，將表現(xiàn)為斷層的低相干的正值全部保留，得到相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體五，圖5為數(shù)據(jù)體五的2800ms時(shí)間切片，高相干的因素全部被充零。
[0032]步驟六:將相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體五進(jìn)行螞蟻?zhàn)粉欉\(yùn)算，得到相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六；例如:將相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體五進(jìn)行螞蟻?zhàn)粉欉\(yùn)算，為了獲得最佳效果，運(yùn)算方法采用非侵略性蟻群算法，得到相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六，圖6為數(shù)據(jù)體六的2800ms時(shí)間切片，此時(shí)斷層非常清晰。
[0033]在上述實(shí)施例中，所述步驟六中的蟻群算法采用非侵略性蟻群算法和侵略性算法可得到需要的螞蟻?zhàn)粉檾?shù)據(jù)體六；其中非侵略性蟻群算法效果較好，但本發(fā)明并不限制該蟻群算法僅僅為非侵略性蟻群算法。
[0034]步驟七:將相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六與原始地震進(jìn)行融合，生成融合數(shù)據(jù)體七；例如:將相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六與原始地震進(jìn)行融合，為了獲得較好的效果，融合時(shí)螞蟻?zhàn)粉檾?shù)據(jù)體范圍選擇為-0.7-1，為了獲得較快的運(yùn)算速度，原始地震數(shù)據(jù)體選擇為8位數(shù)據(jù)體，生成融合數(shù)據(jù)體七。
[0035]最后，利用在原始地震數(shù)據(jù)體上解釋的層位數(shù)據(jù)在融合數(shù)據(jù)體七上提取沿層切片，得到斷層與原始地震平面圖，圖7為平面圖結(jié)果，可以用于斷層解釋。
[0036]通過(guò)上述實(shí)施例，在地震資料品質(zhì)較差的情況下，通過(guò)本發(fā)明能夠得到更好的識(shí)別斷層的數(shù)據(jù)體，與實(shí)際鉆井情況吻合較好。
[0037]以上所述的具體實(shí)施例，對(duì)本發(fā)明的目的、技術(shù)方案和有益效果進(jìn)行了進(jìn)一步詳細(xì)說(shuō)明，所應(yīng)理解的是，以上所述僅為本發(fā)明的具體實(shí)施例而已，并不用于限定本發(fā)明的保護(hù)范圍，凡在本發(fā)明的精神和原則之內(nèi)，所做的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。
【權(quán)利要求】
1.一種用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述處理方法包含: 步驟一:通過(guò)鉆井分層、聲波時(shí)差測(cè)井曲線、密度測(cè)井曲線制作合成地震記錄，并進(jìn)行井震層位標(biāo)定，通過(guò)研究區(qū)內(nèi)已知井資料建立井震關(guān)系，確定目的層段在地震數(shù)據(jù)體上的雙程反射時(shí)間范圍； 步驟二:將目的層段對(duì)應(yīng)的原始地震資料置于所述雙程反射時(shí)間范圍內(nèi)，以所述雙程反射時(shí)間范圍為時(shí)窗對(duì)目的層段對(duì)應(yīng)的原始地震資料進(jìn)行構(gòu)造導(dǎo)向?yàn)V波，得到導(dǎo)向?yàn)V波地震資料數(shù)據(jù)體一； 步驟三:對(duì)導(dǎo)向?yàn)V波地震資料數(shù)據(jù)體一進(jìn)行相干運(yùn)算，得到構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體-* ; 步驟四:對(duì)構(gòu)造導(dǎo)向?yàn)V波相干數(shù)據(jù)體二進(jìn)行反射強(qiáng)度交流分量的運(yùn)算，得到相干體反射強(qiáng)度交流分量數(shù)據(jù)體三，將相干體反射強(qiáng)度交流分量數(shù)據(jù)體三正負(fù)對(duì)調(diào)，得到正負(fù)值對(duì)調(diào)的相干反射強(qiáng)度交流分量數(shù)據(jù)體四； 步驟五:將正負(fù)值對(duì)調(diào)的相干反射強(qiáng)度交流分量數(shù)據(jù)體四進(jìn)行數(shù)據(jù)體賦值運(yùn)算，將表現(xiàn)為地層因素的高相干的負(fù)值全部賦值為零，將表現(xiàn)為斷層的低相干的正值全部保留，得到相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體五； 步驟六:將相干反射強(qiáng)度交流分量的正值數(shù)據(jù)體五進(jìn)行螞蟻?zhàn)粉欉\(yùn)算，得到相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六； 步驟七:將相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六與原始地震進(jìn)行融合，生成融合數(shù)據(jù)體七； 步驟八:利用在原始地震數(shù)據(jù)體上解釋的層位數(shù)據(jù)在融合數(shù)據(jù)體七上提取沿層切片，得到斷層與原始地震平面圖。
2.根據(jù)權(quán)利要求1所述的用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述導(dǎo)向?yàn)V波的算法包含:中值濾波或均值濾波或主分量濾波。
3.根據(jù)權(quán)利要求1所述的用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述步驟三中的相干算法包含:本征值相干或相關(guān)或相似性相干。
4.根據(jù)權(quán)利要求3所述的用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述相干算法包含:通過(guò)5點(diǎn)相干或7點(diǎn)相干或9點(diǎn)相干,得到構(gòu)造導(dǎo)向?yàn)V波本征值相干數(shù)據(jù)體二。
5.根據(jù)權(quán)利要求1所述的用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述步驟四中還包含:數(shù)據(jù)體三在數(shù)值上是正值和負(fù)值相間的數(shù)據(jù)體，低相干的斷層表現(xiàn)為負(fù)值，高相干的地層信息表現(xiàn)為正值。
6.根據(jù)權(quán)利要求5所述的用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述數(shù)據(jù)體四中，斷層表現(xiàn)為正值，地層信息表現(xiàn)為負(fù)值。
7.根據(jù)權(quán)利要求1所述的用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述步驟六中的蟻群算法包含:非侵略性蟻群算法或侵略性算法，得到相干反射強(qiáng)度交流分量正值體的螞蟻?zhàn)粉檾?shù)據(jù)體六。
8.根據(jù)權(quán)利要求1所述的用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述步驟七中還包含:融合時(shí)螞蟻?zhàn)粉檾?shù)據(jù)體范圍選擇為-0.7-1 ο
9.根據(jù)權(quán)利要求8所述的用于識(shí)別斷層的相干增強(qiáng)處理方法，其特征在于，所述原始地震數(shù)據(jù)體選擇為8位數(shù)據(jù)體，生成融合數(shù)據(jù)體七。
【文檔編號(hào)】G06F19/00GK104459801SQ201410756927
【公開日】2015年3月25日 申請(qǐng)日期:2014年12月10日 優(yōu)先權(quán)日:2014年12月10日
【發(fā)明者】陳志剛, 馬輝, 唐建超, 王亞玲, 潘良云 申請(qǐng)人:中國(guó)石油天然氣集團(tuán)公司, 中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司