一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法及裝置制造方法
【專利摘要】本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法及裝置,其中,所述評(píng)價(jià)方法包括:根據(jù)勘探地區(qū)的地質(zhì)任務(wù)設(shè)計(jì)觀測系統(tǒng),并模擬放炮,得到觀測系統(tǒng)的每個(gè)面元的炮檢距分布;選取代表勘探地區(qū)特征的炮集記錄;在所述炮集記錄上選定出地震反射波時(shí)窗;計(jì)算地震反射波時(shí)窗內(nèi)炮集記錄中的每一道地震數(shù)據(jù)的地震屬性值;對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式;根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值;根據(jù)所述每個(gè)面元的地震屬性值獲取觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差;根據(jù)所述觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差對(duì)所述觀測系統(tǒng)進(jìn)行評(píng)價(jià);其中,所述標(biāo)準(zhǔn)方差越小,表明所述觀測系統(tǒng)的地震面元屬性越好。
【專利說明】一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法及裝置【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及石油勘探領(lǐng)域,特別涉及一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法及裝置。
【背景技術(shù)】
[0002]地震勘探是石油勘探的重要手段。地震數(shù)據(jù)勘探采集技術(shù)設(shè)計(jì)是為保證探區(qū)地質(zhì)勘探設(shè)計(jì)提出的地質(zhì)任務(wù)、技術(shù)要求、勘探時(shí)間要求而進(jìn)行的方法論證、技術(shù)分析等工作。其主要內(nèi)容包括勘探任務(wù)來源、收集探區(qū)以往資料、分析探區(qū)情況、勘探難點(diǎn)與技術(shù)對(duì)策、采集方法設(shè)計(jì)與觀測系統(tǒng)方案等內(nèi)容。觀測系統(tǒng)設(shè)計(jì)是否合理是影響地震勘探效果的關(guān)鍵因素,因此如何評(píng)價(jià)觀測系統(tǒng)設(shè)計(jì)方案的優(yōu)劣非常重要。
[0003]目前國內(nèi)外觀測系統(tǒng)的優(yōu)劣評(píng)價(jià)方法主要有三種:
[0004]1、基于目的層CMP面元信息分析的觀測系統(tǒng)評(píng)價(jià)方法。該方法采取計(jì)算地表和目的層的水平層狀假設(shè),通過計(jì)算規(guī)則共中心點(diǎn)CMP面元的覆蓋次數(shù)分布、炮檢距分布、方位角分布來評(píng)價(jià)觀測系統(tǒng) 。基于水平層狀介質(zhì)假設(shè)的常規(guī)三維地震勘探觀測系統(tǒng)設(shè)計(jì)方法難以適應(yīng)劇烈起伏地表和復(fù)雜構(gòu)造區(qū)。實(shí)際上當(dāng)?shù)乇砬闆r和地下構(gòu)造比較復(fù)雜時(shí),CMP共中心點(diǎn)面元與目的層共反射點(diǎn)(CRP)面元有很大的差別。
[0005]2、基于目的層CRP面元信息分析的觀測系統(tǒng)評(píng)價(jià)方法。通過建立工區(qū)的二維或者三維地質(zhì)模型,利用射線追蹤技術(shù)對(duì)地下目的層進(jìn)行射線追蹤,然后根據(jù)追蹤結(jié)果計(jì)算目的層的CRP面元的屬性,通過分析CRP面元屬性的優(yōu)劣來判斷是否可以得到目的層的反射資料,從而判斷設(shè)計(jì)方案是否合理。利用射線追蹤方法進(jìn)行CRP地震屬性分析,雖然能較準(zhǔn)確地提供地下各界面射線參數(shù),但要求速度模型橫向變化不能太大,其適用性和有效性由于存在高頻近似和射線理論本身的缺陷在復(fù)雜構(gòu)造區(qū)受到挑戰(zhàn)。
[0006]3、基于目的層照明信息的觀測系統(tǒng)評(píng)價(jià)。通過建立二維或者三維地質(zhì)模型,利用波動(dòng)方程照明分析技術(shù)計(jì)算觀測系統(tǒng)對(duì)地下目的層的照明強(qiáng)度,通過分析目的層的照明情況也評(píng)價(jià)觀測系統(tǒng)的優(yōu)劣。該方法的計(jì)算量比較大,實(shí)用性不是很強(qiáng),并且其結(jié)果依賴模型的精確性。
[0007]以上對(duì)三種觀測系統(tǒng)的評(píng)價(jià)分析技術(shù)基本上從理論參數(shù)和地質(zhì)模型出發(fā),通過抽取共中心點(diǎn)、射線追蹤或波場數(shù)值模擬,以CMP、CRP面元屬性和照明為分析對(duì)象來優(yōu)化觀測系統(tǒng),缺少與以往實(shí)際地震數(shù)據(jù)的融合論證分析。
[0008]目前各油田多為勘探開發(fā)的高成熟區(qū),基本上為三維所覆蓋,具有豐富的2D、3D地震數(shù)據(jù),但這些資料在實(shí)際的采集設(shè)計(jì)過程中往往很少被采用。常規(guī)地震資料包含眾多的野外采集信息,利用這些信息評(píng)價(jià)基于地震屬性的觀測系統(tǒng)是石油勘探開發(fā)中急需重點(diǎn)解決的關(guān)鍵性難題。
【發(fā)明內(nèi)容】
[0009]為解決上述問題,本發(fā)明提出一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法及裝置,通過探區(qū)已有的實(shí)際地震數(shù)據(jù)計(jì)算描述每個(gè)面元的地震屬性值,根據(jù)地震屬性值分布評(píng)價(jià)觀測系統(tǒng),使得評(píng)價(jià)更加科學(xué)、合理。
[0010]為實(shí)現(xiàn)上述目的,本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法,所述評(píng)價(jià)方法包括:
[0011]根據(jù)勘探地區(qū)的地質(zhì)任務(wù)設(shè)計(jì)觀測系統(tǒng),并模擬放炮,得到觀測系統(tǒng)的每個(gè)面元的炮檢距分布;
[0012]選取代表勘探地區(qū)特征的炮集記錄;
[0013]在所述炮集記錄上選定出地震反射波時(shí)窗;
[0014]計(jì)算地震反射波時(shí)窗內(nèi)炮集記錄中的每一道地震數(shù)據(jù)的地震屬性值;
[0015]對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式;
[0016]根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值;
[0017]根據(jù)所述每個(gè)面元的地震屬性值獲取觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差;
[0018]根據(jù)所述觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差對(duì)所述觀測系統(tǒng)進(jìn)行評(píng)價(jià);其中,所述標(biāo)準(zhǔn)方差越小,表明所述觀測系統(tǒng)的地震面元屬性越好。
[0019]可選的,在本發(fā)明一實(shí)施例中,所述勘探地區(qū)特征為山地、丘陵、河流、湖泊、平原、海洋、濕地、沙漠或城鎮(zhèn)。
[0020]可選的,在本發(fā)明一實(shí)施例中,所述地震屬性值為所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的能量、所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的主頻或所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的信噪比。
[0021]可選的,在本發(fā)明一實(shí)施例中,所述對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式的步驟包括:
[0022]根據(jù)中值濾波法對(duì)所述地震屬性值進(jìn)行濾波;
[0023]對(duì)濾波后的地震屬性值根據(jù)炮檢距擬合獲取地震屬性擬合表達(dá)式。
[0024]可選的,在本發(fā)明一實(shí)施例中,所述根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值的步驟包括:
[0025]根據(jù)所述地震屬性擬合表達(dá)式,從觀測系統(tǒng)的每個(gè)面元所有的炮檢距中抽取對(duì)應(yīng)的地震屬性值;
[0026]將獲取的地震屬性值相加并除以每個(gè)面元內(nèi)炮檢距個(gè)數(shù),得到觀測系統(tǒng)中每個(gè)面元的地震屬性值。
[0027]為實(shí)現(xiàn)上述目的,本發(fā)明還提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)裝置,所述評(píng)價(jià)裝置包括:
[0028]面元的炮檢距分布獲取單元,用于根據(jù)勘探地區(qū)的地質(zhì)任務(wù)設(shè)計(jì)觀測系統(tǒng),并模擬放炮,得到觀測系統(tǒng)的每個(gè)面元的炮檢距分布;
[0029]炮集記錄獲取單元,用于選取代表勘探地區(qū)特征的炮集記錄;
[0030]地震反射波時(shí)窗獲取單元,用于在所述炮集記錄上選定出地震反射波時(shí)窗;
[0031]時(shí)窗內(nèi)地震屬性值獲取單元,用于計(jì)算地震反射波時(shí)窗內(nèi)炮集記錄中的每一道地震數(shù)據(jù)的地震屬性值;
[0032]地震屬性擬合表達(dá)式獲取單元,用于對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式;[0033]面元的地震屬性值獲取單元,用于根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值;
[0034]標(biāo)準(zhǔn)方差獲取單元,用于根據(jù)所述每個(gè)面元的地震屬性值獲取觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差;
[0035]評(píng)價(jià)單元,用于根據(jù)所述觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差對(duì)所述觀測系統(tǒng)進(jìn)行評(píng)價(jià);其中,所述標(biāo)準(zhǔn)方差越小,表明所述觀測系統(tǒng)的地震面元屬性越好。
[0036]可選的,在本發(fā)明一實(shí)施例中,所述炮集記錄獲取單元選取勘探地區(qū)特征為山地、丘陵、河流、湖泊、平原、海洋、濕地、沙漠或城鎮(zhèn)。
[0037]可選的,在本發(fā)明一實(shí)施例中,所述時(shí)窗內(nèi)地震屬性值獲取單元獲取的地震屬性值為所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的能量、所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的主頻或所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的信噪比。
[0038]可選的,在本發(fā)明一實(shí)施例中,所述地震屬性擬合表達(dá)式獲取單元包括:
[0039]濾波模塊,用于根據(jù)中值濾波法對(duì)所述地震屬性值進(jìn)行濾波;
[0040]擬合模塊,用于對(duì)濾波后的地震屬性值根據(jù)炮檢距擬合獲取地震屬性擬合表達(dá)式。
[0041]可選的,在本發(fā)明一實(shí)施例中,所述面元的地震屬性值獲取單元包括:
[0042]抽取模塊,用于根據(jù)所述地震屬性擬合表達(dá)式,從觀測系統(tǒng)的每個(gè)面元所有的炮檢距中抽取對(duì)應(yīng)的地震屬性值;
[0043]計(jì)算模塊,用于將獲取的地震屬性值相加并除以每個(gè)面元內(nèi)炮檢距個(gè)數(shù),得到觀測系統(tǒng)中每個(gè)面元的地震屬性值。
[0044]上述技術(shù)方案具有如下有益效果:
[0045]本發(fā)明的技術(shù)方案為觀測系統(tǒng)提供了一種更科學(xué)、定量的評(píng)價(jià)方法,減少由于觀測系統(tǒng)不合理產(chǎn)生的“采集腳印”,提高復(fù)雜地區(qū)地震資料的采集效果。
【專利附圖】
【附圖說明】
[0046]為了更清楚地說明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案,下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡單地介紹,顯而易見地,下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例,對(duì)于本領(lǐng)域普通技術(shù)人員來講,在不付出創(chuàng)造性勞動(dòng)的前提下,還可以根據(jù)這些附圖獲得其他的附圖。
[0047]圖1為本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法流程圖;
[0048]圖2為本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)裝置框圖;
[0049]圖3為本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)裝置中地震屬性擬合表達(dá)式獲取單元框圖;
[0050]圖4為本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)裝置中面元的地震屬性值獲取單元框圖;
[0051]圖5為本實(shí)施例的12線4炮120道正交觀測系統(tǒng)I模版示意圖;
[0052]圖6為本實(shí)施例的12線24炮120道正交觀測系統(tǒng)2模版示意圖;
[0053]圖7為本實(shí)施例的從原始資料中選取代表工區(qū)特征的炮集記錄示意圖;
[0054]圖8為本實(shí)施例的基于圖7所示炮集記錄上用時(shí)窗選定所要研究的目的層反射波示意圖;
[0055]圖9為本實(shí)施例中時(shí)窗內(nèi)炮點(diǎn)右邊接收道均方根能量擬合曲線示意圖;
[0056]圖10為本實(shí)施例中時(shí)窗內(nèi)炮點(diǎn)右邊接收道主頻擬合曲線示意圖;
[0057]圖11為本實(shí)施例中觀測系統(tǒng)I面元內(nèi)均方根能量疊加平面示意圖;
[0058]圖12為本實(shí)施例中基于圖11的一部分放大面元內(nèi)均方根能量疊加平面圖;
[0059]圖13為本實(shí)施例中觀測系統(tǒng)I面元內(nèi)主頻疊加平面示意圖;
[0060]圖14為本實(shí)施例中基于圖13的一部分放大面元內(nèi)主頻疊加平面圖;
[0061]圖15為本實(shí)施例中觀測系統(tǒng)2面元內(nèi)均方根能量疊加平面示意圖;
[0062]圖16為本實(shí)施例中基于圖15的一部分放大面兀內(nèi)均方根能量疊加平面圖;
[0063]圖17為本實(shí)施例中觀測系統(tǒng)2面元內(nèi)主頻疊加平面示意圖;
[0064]圖18為本實(shí)施例中基于圖17的一部分放大面元內(nèi)主頻疊加平面圖;
【具體實(shí)施方式】
[0065]下面將結(jié)合本發(fā)明實(shí)施例中的附圖,對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述。顯然,所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例,而不是全部的實(shí)施例。基于本發(fā)明中的實(shí)施例,本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例,都屬于本發(fā)明保護(hù)的范圍。
[0066]本發(fā)明是石油地球物理勘探中一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)的技術(shù)方案。首先獲取設(shè)計(jì)的觀測系統(tǒng)每一個(gè)面元的炮檢距分布,其次從工區(qū)采集的實(shí)際地震資料中選取能代表工區(qū)每個(gè)區(qū)塊特征的炮集或CMP道集,然后在選取的炮集記錄上針對(duì)需要分析的地震波類型選擇時(shí)窗,計(jì)算每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的地震屬性值(包括能量、頻率、信噪比等)并進(jìn)行濾波,對(duì)濾波后的地震屬性值按炮檢距進(jìn)行曲線擬合,求取地震屬性擬合表達(dá)式,利用地震屬性擬合表達(dá)式對(duì)設(shè)計(jì)的觀測系統(tǒng)每個(gè)面元所有的炮檢距抽取對(duì)應(yīng)的地震屬性值,將地震屬性值相加并除以面元內(nèi)炮檢距個(gè)數(shù),得到每一個(gè)面元的地震屬性值,根據(jù)觀測系統(tǒng)中面元的地震屬性值獲取地震屬性值分布的標(biāo)準(zhǔn)方差,通過地震屬性值分布的標(biāo)準(zhǔn)方差確定合理的觀測系統(tǒng)。
[0067]如圖1所示,為本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法流程圖。所述評(píng)價(jià)方法包括:
[0068]步驟101):根據(jù)勘探地區(qū)的地質(zhì)任務(wù)設(shè)計(jì)觀測系統(tǒng),并模擬放炮,得到觀測系統(tǒng)的每個(gè)面元的炮檢距分布;
[0069]步驟102):選取代表勘探地區(qū)特征的炮集記錄;
[0070]步驟103):在所述炮集記錄上選定出地震反射波時(shí)窗;
[0071]步驟104):計(jì)算地震反射波時(shí)窗內(nèi)炮集記錄中的每一道地震數(shù)據(jù)的地震屬性值;
[0072]步驟105):對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式;
[0073]步驟106):根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值;
[0074]步驟107):根據(jù)所述每個(gè)面元的地震屬性值獲取觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差;
[0075]步驟108):根據(jù)所述觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差對(duì)所述觀測系統(tǒng)進(jìn)行評(píng)價(jià);其中,所述標(biāo)準(zhǔn)方差越小,表明所述觀測系統(tǒng)的地震面元屬性越好。[0076]可選的,在本發(fā)明一實(shí)施例中,所述勘探地區(qū)特征為山地、丘陵、河流、湖泊、平原、海洋、濕地、沙漠或城鎮(zhèn)。
[0077]可選的,在本發(fā)明一實(shí)施例中,所述地震屬性值為所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的能量、所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的主頻或所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的信噪比。
[0078]可選的,在本發(fā)明一實(shí)施例中,所述對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式的步驟包括:
[0079]根據(jù)中值濾波法對(duì)所述地震屬性值進(jìn)行濾波;
[0080]對(duì)濾波后的地震屬性值根據(jù)炮檢距擬合獲取地震屬性擬合表達(dá)式。
[0081]可選的,在本發(fā)明一實(shí)施例中,所述根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值的步驟包括:
[0082]根據(jù)所述地震屬性擬合表達(dá)式,從觀測系統(tǒng)的每個(gè)面元所有的炮檢距中抽取對(duì)應(yīng)的地震屬性值;
[0083]將獲取的地震屬性值相加并除以每個(gè)面元內(nèi)炮檢距個(gè)數(shù),得到觀測系統(tǒng)中每個(gè)面元的地震屬性值。
[0084]如圖2所示,為本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)裝置框圖。所述評(píng)價(jià)裝置包括:
[0085]面元的炮檢距分布獲取單元201,用于根據(jù)勘探地區(qū)的地質(zhì)任務(wù)設(shè)計(jì)觀測系統(tǒng),并模擬放炮,得到觀測系統(tǒng)的每個(gè)面元的炮檢距分布;
[0086]炮集記錄獲取單元202,用于選取代表勘探地區(qū)特征的炮集記錄;
[0087]地震反射波時(shí)窗獲取單元203,用于在所述炮集記錄上選定出地震反射波時(shí)窗;
[0088]時(shí)窗內(nèi)地震屬性值獲取單元204,用于計(jì)算地震反射波時(shí)窗內(nèi)炮集記錄中的每一道地震數(shù)據(jù)的地震屬性值;
[0089]地震屬性擬合表達(dá)式獲取單元205,用于對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式;
[0090]面元的地震屬性值獲取單元206,用于根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值;
[0091]標(biāo)準(zhǔn)方差獲取單元207,用于根據(jù)所述每個(gè)面元的地震屬性值獲取觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差;
[0092]評(píng)價(jià)單元208,用于根據(jù)所述觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差對(duì)所述觀測系統(tǒng)進(jìn)行評(píng)價(jià);其中,所述標(biāo)準(zhǔn)方差越小,表明所述觀測系統(tǒng)的地震面元屬性越好。
[0093]可選的,在本發(fā)明一實(shí)施例中,所述炮集記錄獲取單元選取勘探地區(qū)特征為山地、丘陵、河流、湖泊、平原、海洋、濕地、沙漠或城鎮(zhèn)。
[0094]可選的,在本發(fā)明一實(shí)施例中,所述時(shí)窗內(nèi)地震屬性值獲取單元獲取的地震屬性值為所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的能量、所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的主頻或所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的信噪比。
[0095]如圖3所示,為本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)裝置中地震屬性擬合表達(dá)式獲取單元框圖。所述地震屬性擬合表達(dá)式獲取單元205包括:
[0096]濾波模塊2051,用于根據(jù)中值濾波法對(duì)所述地震屬性值進(jìn)行濾波;[0097]擬合模塊2052,用于對(duì)濾波后的地震屬性值根據(jù)炮檢距擬合獲取地震屬性擬合表達(dá)式。
[0098]如圖4所示,為本發(fā)明提供了一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)裝置中面元的地震屬性值獲取單元框圖。所述面元的地震屬性值獲取單元206包括:
[0099]抽取模塊2061,用于根據(jù)所述地震屬性擬合表達(dá)式,從觀測系統(tǒng)的每個(gè)面元所有的炮檢距中抽取對(duì)應(yīng)的地震屬性值;
[0100]計(jì)算模塊2062,用于將獲取的地震屬性值相加并除以每個(gè)面元內(nèi)炮檢距個(gè)數(shù),得到觀測系統(tǒng)中每個(gè)面元的地震屬性值。
[0101]實(shí)施例:
[0102]本實(shí)施例的步驟包括:
[0103]步驟I):收集勘探地區(qū)已有的原始資料;
[0104]步驟2):根據(jù)勘探地區(qū)的地質(zhì)任務(wù)設(shè)計(jì)了兩個(gè)觀測系統(tǒng),如圖5和圖6,圖5為12線4炮120道正交觀測系統(tǒng)I,覆蓋次數(shù)為90次,最大炮檢距3199m,圖6為12線24炮120道正交觀測系統(tǒng)2,覆蓋次數(shù)為90次,最大炮檢距3414m,對(duì)設(shè)計(jì)的觀測系統(tǒng)模擬放炮,得到觀測系統(tǒng)中每一個(gè)面元的炮檢距分布;
[0105]步驟3):探區(qū)為平原,從原始資料中選取能代表整個(gè)工區(qū)特征的炮集。如圖7所示,為本實(shí)施例的從原始資料中選取代表工區(qū)特征的炮集記錄示意圖。中間放炮,接收線上220道,道距25m,縱坐標(biāo)為時(shí)間,單位為ms,橫坐標(biāo)為道號(hào)。
[0106]步驟4):在圖7所示炮集記錄上用時(shí)窗選定所要研究的目的層反射波,如圖8所示,為本實(shí)施例的基于圖7所示炮集記錄上用時(shí)窗選定所要研究的目的層反射波示意圖。兩條雙曲線中間的位置為選定的目的層位置,縱坐標(biāo)為時(shí)間,單位為ms,橫坐標(biāo)為道號(hào)。在本實(shí)施例中,計(jì)算炮點(diǎn)右邊接收道每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的均方根能量和主頻;
[0107]步驟5):利用三點(diǎn)中值濾波對(duì)步驟4)中求出的均方根能量和主頻進(jìn)行濾波,濾波后均方根能量為:2.27,4.10,……,0.23,0.18,如圖9中“.”所示,主頻為:19.74,17.02,……,8.83,7.06,如圖10中“.”所示,對(duì)濾波后的數(shù)據(jù)擬合,最大炮檢距 2737.5m。
[0108]如圖9所示,為本實(shí)施例中時(shí)窗內(nèi)炮點(diǎn)右邊接收道均方根能量擬合曲線示意圖。圖9中均方根能量的擬合表達(dá)式為:
[0109]y = 9.10-η.χ6_3.10-8.χ5+4.10-6.χ4_0.0OO2.χ3+0.0054.χ2_0.0237.χ+3.3648
[0110]上式中,X為炮檢距,y為均方根能量。
[0111]如圖10所示,為本實(shí)施例中時(shí)窗內(nèi)炮點(diǎn)右邊接收道主頻擬合曲線示意圖。圖10中主頻的擬合表達(dá)式為:
[0112]y = -10-9.χ6+3.10-7.χ5_4.10-5.χ4+0.0016.χ3_0.0228.χ2_0.2124.χ+18.33
[0113]上式中,X為炮檢距,y為主頻。
[0114]步驟6):對(duì)步驟2)設(shè)計(jì)的觀測系統(tǒng)模擬放炮得到的每個(gè)面元所有的炮檢距按步驟5)中的擬合表達(dá)式抽取對(duì)應(yīng)的均方根能量或主頻,每個(gè)面元所有的均方根能量或主頻相加并除以面元內(nèi)炮檢距個(gè)數(shù),得到每個(gè)面元的均方根能量或主頻,如圖11所示,為本實(shí)施例中觀測系統(tǒng)I面元內(nèi)均方根能量疊加平面示意圖;如圖12所示,為本實(shí)施例中基于圖11的一部分放大面元內(nèi)均方根能量疊加平面圖;如圖13所示,為本實(shí)施例中觀測系統(tǒng)I面元內(nèi)主頻疊加平面示意圖;如圖14所示,為本實(shí)施例中基于圖13的一部分放大面元內(nèi)主頻疊加平面圖;如圖15所示,為本實(shí)施例中觀測系統(tǒng)2面元內(nèi)均方根能量疊加平面示意圖;如圖16所不,為本實(shí)施例中基于圖15的一部分放大面兀內(nèi)均方根能量疊加平面圖;如圖17所示,為本實(shí)施例中觀測系統(tǒng)2面元內(nèi)主頻疊加平面示意圖;如圖18所示,為本實(shí)施例中基于圖17的一部分放大面元內(nèi)主頻疊加平面圖。
[0115]比較面元均方根能量和主頻在兩個(gè)觀測系統(tǒng)中的分布,觀測系統(tǒng)2中的均方根能量分布圖15和主頻分布圖17中出現(xiàn)了較強(qiáng)的“采集腳印”。
[0116]步驟7):計(jì)算兩種觀測系統(tǒng)地震面元屬性值分布的標(biāo)準(zhǔn)方差,觀測系統(tǒng)I對(duì)應(yīng)的面元均方根能量標(biāo)準(zhǔn)方差為0.1965,主頻的標(biāo)準(zhǔn)方差為0.2831 ;觀測系統(tǒng)2對(duì)應(yīng)的面元均方根能量標(biāo)準(zhǔn)方差為0.2174,主頻的標(biāo)準(zhǔn)方差為0.3017,觀測系統(tǒng)I對(duì)應(yīng)的面元均方根能量和主頻的標(biāo)準(zhǔn)方差分別比觀測系統(tǒng)2對(duì)應(yīng)的面元均方根能量和主頻的標(biāo)準(zhǔn)方差小。觀測系統(tǒng)地震面元屬性值分布的標(biāo)準(zhǔn)方差越小,表明觀測系統(tǒng)面元的地震屬性值波動(dòng)范圍越小,觀測系統(tǒng)地震面元屬性值分布越均勻,標(biāo)準(zhǔn)方差小的觀測系統(tǒng)具有較好的地震面元屬性值分布,因此,觀測系統(tǒng)I采集更合理。
[0117]最后應(yīng)說明的是:上述僅用以說明本發(fā)明而并非限制本發(fā)明所描述的技術(shù)方案;盡管本說明書對(duì)本發(fā)明已進(jìn)行了詳細(xì)的說明,但是,本領(lǐng)域的技術(shù)人員仍然可以對(duì)本發(fā)明進(jìn)行修改或等同替換,一切不脫離本發(fā)明的精神和范圍的技術(shù)方案及其改進(jìn),其均應(yīng)涵蓋在本發(fā)明的權(quán)利要求范圍中。
【權(quán)利要求】
1.一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)方法,其特征在于,所述評(píng)價(jià)方法包括: 根據(jù)勘探地區(qū)的地質(zhì)任務(wù)設(shè)計(jì)觀測系統(tǒng),并模擬放炮,得到觀測系統(tǒng)的每個(gè)面元的炮檢距分布; 選取代表勘探地區(qū)特征的炮集記錄; 在所述炮集記錄上選定出地震反射波時(shí)窗; 計(jì)算地震反射波時(shí)窗內(nèi)炮集記錄中的每一道地震數(shù)據(jù)的地震屬性值; 對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式; 根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值; 根據(jù)所述每個(gè)面元的地震屬性值獲取觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差; 根據(jù)所述觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差對(duì)所述觀測系統(tǒng)進(jìn)行評(píng)價(jià);其中,所述標(biāo)準(zhǔn)方差越小,表明所述觀測系統(tǒng)的地震面元屬性越好。
2.如權(quán)利要求1所述的方法,其特征在于,所述勘探地區(qū)特征為山地、丘陵、河流、湖泊、平原、海洋、濕地、沙漠或城鎮(zhèn)。
3.如權(quán)利要求1所述的方法,其特征在于,所述地震屬性值為所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的能量、所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的主頻或所選炮集記錄中的每一 道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的信噪比。
4.如權(quán)利要求1所述的方法,其特征在于,所述對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式的步驟包括: 根據(jù)中值濾波法對(duì)所述地震屬性值進(jìn)行濾波; 對(duì)濾波后的地震屬性值根據(jù)炮檢距擬合獲取地震屬性擬合表達(dá)式。
5.如權(quán)利要求1所述的方法,其特征在于,所述根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值的步驟包括: 根據(jù)所述地震屬性擬合表達(dá)式,從觀測系統(tǒng)的每個(gè)面元所有的炮檢距中抽取對(duì)應(yīng)的地震屬性值; 將獲取的地震屬性值相加并除以每個(gè)面元內(nèi)炮檢距個(gè)數(shù),得到觀測系統(tǒng)中每個(gè)面元的地震屬性值。
6.一種基于地震屬性的觀測系統(tǒng)評(píng)價(jià)裝置,其特征在于,所述評(píng)價(jià)裝置包括: 面元的炮檢距分布獲取單元,用于根據(jù)勘探地區(qū)的地質(zhì)任務(wù)設(shè)計(jì)觀測系統(tǒng),并模擬放炮,得到觀測系統(tǒng)的每個(gè)面元的炮檢距分布; 炮集記錄獲取單元,用于選取代表勘探地區(qū)特征的炮集記錄; 地震反射波時(shí)窗獲取單元,用于在所述炮集記錄上選定出地震反射波時(shí)窗; 時(shí)窗內(nèi)地震屬性值獲取單元,用于計(jì)算地震反射波時(shí)窗內(nèi)炮集記錄中的每一道地震數(shù)據(jù)的地震屬性值; 地震屬性擬合表達(dá)式獲取單元,用于對(duì)所述地震屬性值進(jìn)行炮檢距擬合求取地震屬性擬合表達(dá)式; 面元的地震屬性值獲取單元,用于根據(jù)所述地震屬性擬合表達(dá)式獲取每個(gè)面元的地震屬性值; 標(biāo)準(zhǔn)方差獲取單元,用于根據(jù)所述每個(gè)面元的地震屬性值獲取觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差;評(píng)價(jià)單元,用于根據(jù)所述觀測系統(tǒng)的地震屬性值分布的標(biāo)準(zhǔn)方差對(duì)所述觀測系統(tǒng)進(jìn)行評(píng)價(jià);其中,所述標(biāo)準(zhǔn)方差越小,表明所述觀測系統(tǒng)的地震面元屬性越好。
7.如權(quán)利要求6所述的裝置,其特征在于,所述炮集記錄獲取單元選取勘探地區(qū)特征為山地、丘陵、河流、湖泊、平原、海洋、濕地、沙漠或城鎮(zhèn)。
8.如權(quán)利要求6所述的裝置,其特征在于,所述時(shí)窗內(nèi)地震屬性值獲取單元獲取的地震屬性值為所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的能量、所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的主頻或所選炮集記錄中的每一道地震數(shù)據(jù)在選定時(shí)窗內(nèi)的信噪比。
9.如權(quán)利要求6所述的裝置,其特征在于,所述地震屬性擬合表達(dá)式獲取單元包括: 濾波模塊,用于根據(jù)中值濾波法對(duì)所述地震屬性值進(jìn)行濾波; 擬合模塊,用于對(duì)濾波后的地震屬性值根據(jù)炮檢距擬合獲取地震屬性擬合表達(dá)式。
10.如權(quán)利要求6所述的裝置,其特征在于,所述面元的地震屬性值獲取單元包括: 抽取模塊,用于根據(jù)所述地震屬性擬合表達(dá)式,從觀測系統(tǒng)的每個(gè)面元所有的炮檢距中抽取對(duì)應(yīng)的地震屬性值; 計(jì)算模塊,用于將獲取的地震屬性值相加并除以每個(gè)面元內(nèi)炮檢距個(gè)數(shù),得到觀測系統(tǒng)中每個(gè)面元的地震屬性值?!?br>
【文檔編號(hào)】G01V13/00GK103592698SQ201310531176
【公開日】2014年2月19日 申請(qǐng)日期:2013年10月31日 優(yōu)先權(quán)日:2013年10月31日
【發(fā)明者】許銀坡, 鄒雪峰, 詹仕凡, 李培明, 蔣先藝, 楊劍, 宋衛(wèi)鋒, 于香妮 申請(qǐng)人:中國石油天然氣集團(tuán)公司, 中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司