井控能量補(bǔ)償文獻(xiàn)調(diào)研

1、疊前處理部分：目前，將井資料與地震數(shù)據(jù)結(jié)合，來提高地面地震數(shù)據(jù)的處理質(zhì)量的方法大都基于疊后地震數(shù)據(jù)來開展，具體包括井(vsp)控的q補(bǔ)償方法，井(聲波、密度)約束的脈沖反褶積方法，井(聲波)約束的小波變換提高分辨率方法，井約束的零相位化處理等。國內(nèi)文獻(xiàn)調(diào)研：1井約束的能量補(bǔ)償方法研究1.1 球面擴(kuò)散補(bǔ)償?shù)卣鸩ㄔ诘叵陆橘|(zhì)中傳播時(shí)，由于球面擴(kuò)散和地層吸收衰減，振幅減弱，在地震資料處理中，須對(duì)其進(jìn)行補(bǔ)償，這是地震數(shù)據(jù)處理中的必要工作。真振幅恢復(fù)是從地面檢波器記錄到的振幅中消除波前擴(kuò)散的影響，使其恢復(fù)到僅與地下反射界面的反射系數(shù)大小有關(guān)的真振幅值。真振幅的恢復(fù)是基于球面波傳播假設(shè)，利用時(shí)變?cè)鲆孢M(jìn)行恢復(fù)

2、， 即： （1）其中：a為球面擴(kuò)散因子，t為記錄時(shí)間，t0為每個(gè)記錄道的初至?xí)r間?？紤]到球面擴(kuò)散和透射損失，利用指數(shù)增益可以模擬球面波場。井控處理真振幅恢復(fù)的方法是：用零井源距vsp資料求取球面擴(kuò)散補(bǔ)償因子（tar因子），然后對(duì)地震資料進(jìn)行真振幅恢復(fù)。圖1 (a)球面擴(kuò)散補(bǔ)償前單炮記錄 (b)球面擴(kuò)散補(bǔ)償后單炮記錄如圖1所示，應(yīng)用該方法后淺、中、深層均得到較好的補(bǔ)償，低、中、高頻能量總體均衡。1.2 q補(bǔ)償方法地震波在地下介質(zhì)中傳播時(shí)，由于球面擴(kuò)散和地層吸收衰減，振幅減弱。球面擴(kuò)散補(bǔ)償是根據(jù)經(jīng)驗(yàn)和試驗(yàn)給定補(bǔ)償參數(shù)，地層衰減的補(bǔ)償則是利用疊加(或偏移)速度進(jìn)行轉(zhuǎn)換，得到q補(bǔ)償參數(shù)，這樣補(bǔ)償精度取

3、決于處理人員的經(jīng)驗(yàn)和疊加速度的精度。利用井(vsp)資料求取q補(bǔ)償因子，并在獲得的q值附近進(jìn)行掃描，然后將做過不同q補(bǔ)償?shù)钠破拭媾c井資料進(jìn)行匹配分析，最終獲得精確的q值。這樣的井控地震處理技術(shù)最大程度地利用已有井的測井或 vsp等資料，將井點(diǎn)數(shù)據(jù)和地面地震數(shù)據(jù)進(jìn)行匹配，從而生成一系列的匹配屬性(包括基于頻率的可信度、可預(yù)測度、測井資料的固有子波和地震子波的傳遞函數(shù)、零相位算子等等)，然后用這些屬性來檢驗(yàn)地震資料與井資料的匹配程度以達(dá)到最佳匹配，最終便可得到高保真、高分辨率的地震資料。vsp采集到的下行波場數(shù)據(jù)是進(jìn)行q分析的理想數(shù)據(jù)。一般情況下，我們選取一定長度的時(shí)窗 (一般包含直達(dá)縱波即可)

4、，利用速度濾波提取下行直達(dá)縱波，然后對(duì)時(shí)窗內(nèi)的直達(dá)縱波數(shù)據(jù)進(jìn)行q分析，利用譜比法來計(jì)算q值。具體來說是針對(duì)某頻率范圍內(nèi)的vsp數(shù)據(jù)，分析頻率隨深度的變化，從而得到一系列不同深度點(diǎn)的 q值，然后統(tǒng)計(jì)出一個(gè)全局q值作為該深度段的q值。假設(shè)vsp記錄中某深度層段，有兩個(gè)時(shí)間點(diǎn)t1和t2，這兩點(diǎn)位置分別做短時(shí)傅里葉變換，得到兩點(diǎn)的振幅譜，s1(w)，s2 (w)，有： （2）對(duì)q() 做直線擬合，求取斜率可得： （3）就是底層等效吸收系數(shù)，t2-t1為選取的兩個(gè)時(shí)段的中心時(shí)間差，因此進(jìn)一步可以求得品質(zhì)因子： （4）在分析q值的過程中，我們一方面利用 vsp資料來估計(jì)井眼周圍的q因子，同時(shí)也利用了地面地

5、震資料在vsp資料提取出的q值附近掃描，并將以不同q補(bǔ)償?shù)玫降钠破拭媾c井資料(vsp )進(jìn)行匹配分析，達(dá)到最佳匹配的就是精確的q值。關(guān)于q值的可預(yù)測度、可信度及匹配度分別定義如下:可預(yù)測度p定義： （5）其中為某一定義的時(shí)窗范圍內(nèi)的井?dāng)?shù)據(jù)的自相關(guān)，為某一定義的時(shí)窗范圍內(nèi)地震資料的自相關(guān)。為井?dāng)?shù)據(jù)和地震資料的互相關(guān)。可信度定義： （6）匹配度定義： （7）其中為地震數(shù)據(jù)和子波的褶積，為井?dāng)?shù)據(jù)和子波的褶積。對(duì)于以上所定義的可預(yù)測度、可信度及匹配度，其值越大，則表示相應(yīng)的q值越佳。q補(bǔ)償即反q濾波是基于波場延拓理論進(jìn)行的，表達(dá)式如下： （9） （10）式中：u（t，）為合成信號(hào)廣義s變換的時(shí)頻譜；

6、u（t+t）為廣義s變換的反變換。圖2 q補(bǔ)償前(左)與q補(bǔ)償后(右)的疊前時(shí)間偏移剖面對(duì)比從q補(bǔ)償前與q補(bǔ)償后的疊前時(shí)間偏移對(duì)比結(jié)果可以看出(如圖2)，補(bǔ)償后比補(bǔ)償前的疊前時(shí)間偏移剖面在信噪比和分辨率方面有了較大程度的提高，同相軸變得豐富。1.3井約束小波變換提高分辨率地震資料分辨率的降低是由于地層對(duì)其中傳播的地震波的吸收所致。地層的吸收作用一方面表現(xiàn)為能量隨傳播時(shí)間衰減，另一方面由于吸收頻率的選擇性，表現(xiàn)為由于波形的延長而導(dǎo)致記錄分辨率的降低。而且由于地層巖性的變化、吸收頻率的選擇性不一致，使記錄在各時(shí)間的分辨率也不同。在真振幅恢復(fù)中，雖然能量隨時(shí)間的衰減已得到了某種補(bǔ)償，但是依賴頻率的衰

7、減依然存在。小波變換是尺度因子a和時(shí)移因子b的二元函數(shù)。即a與頻率相對(duì)應(yīng),b與時(shí)間相對(duì)應(yīng),因此可以把小波變換看成頻率和時(shí)間的二元函數(shù),這就為我們?cè)谛〔ㄗ儞Q域內(nèi)按不同的時(shí)間對(duì)衰減了的頻率進(jìn)行補(bǔ)償提供了思路。根據(jù)測井?dāng)?shù)據(jù)綜合分析知：由于實(shí)際地層的反射系數(shù)隨頻率增高，因而振幅能量呈線性增強(qiáng)趨勢，在小波變換域按這一規(guī)律對(duì)實(shí)際記錄進(jìn)行處理，可以提高分辨率。圖3 聲波測井曲線（a），該測井曲線的反射系數(shù)（b），該測井曲線30hz合成記錄（c），該測井曲線50hz合成記錄（d），經(jīng)小波變換處理后提高分辨率的結(jié)果（e）圖3是理論聲波測井曲線經(jīng)小波變換處理后提高分辨率的結(jié)果與50hz的合成記錄與經(jīng)小波變換處理提

8、高分辨率后的記錄相比十分相近，只是在某些細(xì)節(jié)上略有差異,但從總體上看,已經(jīng)吻合得很好了。這樣就可以利用井合成記錄取得小波處理參數(shù),進(jìn)而將其應(yīng)用于地震剖面處理中。為了將主頻為40hz的地震剖面經(jīng)小波變換處理后，主頻提高60hz，要先對(duì)井聲波曲線進(jìn)行一系列處理；為了做到這一點(diǎn)，首先在井合成記錄上取得合理的處理參數(shù)，然后將該參數(shù)應(yīng)用于地震剖面中，便得到了提高分辨率處理后的剖面，如圖4所示。圖4（a）二維測線與40hz合成記錄的對(duì)比圖圖4（b）經(jīng)小波變換處理后分辨率增加的結(jié)果與60hz合成記錄的對(duì)比圖 從圖可以看出，地震剖面與井合成記錄吻合的相當(dāng)好，從而進(jìn)一步說明，在井聲波測井曲線的約束下，經(jīng)小波變換

9、處理后地震剖面的增幅可以得到合理的控制，使處理結(jié)果的可信度大大增加。1.4 基于emd分解的地震數(shù)據(jù)自動(dòng)增益方法希爾伯特-黃變換( hht) 是由huang 等于1998年提出的，其核心思想是將時(shí)間序列資料通過經(jīng)驗(yàn)?zāi)B(tài)分解( emd) 分解成數(shù)個(gè)固有模態(tài)函數(shù)( imf)，然后利用希爾伯特變換構(gòu)造解析信號(hào)，得出資料的瞬時(shí)頻率和振幅， 進(jìn)而得到希爾伯特譜。地震波在非完全彈性介質(zhì)( 巖性介質(zhì)) 內(nèi)從激發(fā)點(diǎn)向外傳播時(shí)， 將引起地震波的球面發(fā)散和吸收衰減。常規(guī)處理軟件是通過時(shí)間域振幅衰減分析和頻譜域振幅分析來近似描述實(shí)際球面發(fā)散與吸收衰減特性， 但該類分析方法只能提供隨時(shí)間或頻率變化的振幅衰減特性。這種

10、單一的時(shí)間域、頻率域分析難以給出實(shí)際大地球面發(fā)散與吸收衰減的整體直觀信息， 加之近地表干擾波的影響及單道、單時(shí)窗計(jì)算中的誤差， 使得難以達(dá)到分析球面發(fā)散與吸收衰減函數(shù)的目的。從已發(fā)表的文章和實(shí)際應(yīng)用分析軟件來看， 至今仍沒有見到十分有效的從地面采集地震數(shù)據(jù)求取球面發(fā)散與吸收衰減的分析方法。hht振幅頻率恢復(fù)處理技術(shù)從實(shí)際地震數(shù)據(jù)本身出發(fā)，從多尺度時(shí)頻域角度， 在多道統(tǒng)計(jì)分析的基礎(chǔ)上得到實(shí)際反射波的球面發(fā)散與吸收衰減的相對(duì)變化特性。以實(shí)際地震數(shù)據(jù)作為驅(qū)動(dòng)， 避免了以往常規(guī)思路中由于地下介質(zhì)特性未知， 振幅補(bǔ)償結(jié)果人為因素影響較大的問題， 從思路上來說是一個(gè)創(chuàng)新。hht能量補(bǔ)償?shù)幕舅枷胧窍葘?duì)信號(hào)

11、進(jìn)行emd分解，然后就可以對(duì)不同尺度的地震記錄進(jìn)行補(bǔ)償，然后再重構(gòu)， 從而實(shí)現(xiàn)對(duì)地震記錄振幅頻率的恢復(fù)補(bǔ)償。圖5 hht振幅頻率恢復(fù)處理前后結(jié)果對(duì)比圖5是hht振幅頻率恢復(fù)技術(shù)處理前后剖面的比較，處理后數(shù)據(jù)頻帶展寬( 有效頻帶范圍從處理前的8 25hz展寬到處理后的845hz) ，薄層反射突出，地層接觸關(guān)系清楚，處理前剖面上的一些低頻噪音經(jīng)hht振幅頻率恢復(fù)處理后，隨著地震信號(hào)高頻反射信號(hào)的恢復(fù)，低頻噪音相對(duì)得到了壓制，處理后剖面信噪比得到了提高。2.井約束的零相位化處理地震記錄可用地震子波和反射系數(shù)序列的褶積來描述，而反射系數(shù)是地層固有的特性，可以通過測井資料得到比較精確的表達(dá)；地震子波則是

12、提高地震資料分辨率需要考慮的主要因素。在相同帶寬條件下，相對(duì)最小相位子波而言，零相位子波的旁瓣最小，能量最集中，具有較高的分辨率。此外，其峰值時(shí)間就是反射波的到達(dá)時(shí)間，是地震資料解釋中不可缺少的參數(shù)。井約束的零相位化處理思路如下：首先用零相位雷克子波與聲波測井?dāng)?shù)據(jù)(地震反射系數(shù)序列)進(jìn)行褶積，產(chǎn)生零相位合成記錄，其中雷克子波主頻由井旁地震數(shù)據(jù)頻譜分析得到。然后在landmark解釋系統(tǒng)中將該合成記錄進(jìn)行相位旋轉(zhuǎn)掃描，與經(jīng)過+90°相位旋轉(zhuǎn)后的實(shí)際地震數(shù)據(jù)進(jìn)行匹配。最后，在最大相似度時(shí)，求取剩余的相位角。具體的處理流程見圖4所示。圖6零相位化處理流程圖這套流程還有很多需要改進(jìn)和完善的地

13、方：（1）在實(shí)際處理中，我們用零相位雷克子波取代實(shí)際的地震子波生成合成記錄。實(shí)際上，子波在各道之間是變化的，而且是旅行時(shí)間函數(shù)，用雷克子波來取代是近似的，可以考慮針對(duì)井旁數(shù)據(jù)的某個(gè)層位提取平均子波，保證子波具有較好的穩(wěn)定性。（2）在這套處理流程里，我們利用互相關(guān)求最大相似度的方法來確定剩余相位角，其精度不是很高。若應(yīng)用匹配濾波技術(shù)，可以更精確地提供相位差異信息。（3）現(xiàn)有的一些處理軟件在求取匹配算子時(shí)，同時(shí)可提供零相位合成記錄和實(shí)際地震記錄之間的時(shí)移、相移等信息，可將其用于資料的分析，并對(duì)最終的零相位化結(jié)果進(jìn)行驗(yàn)證。3. 井約束脈沖反褶積方法由于地震資料的高頻段能量弱并且受噪聲影響嚴(yán)重，單純的

14、依靠地震的高頻信息進(jìn)行反褶積容易在提高地震資料分辨率的同時(shí)降低資料的信噪比，不利于薄儲(chǔ)層的解釋和儲(chǔ)存預(yù)測。測井信息作為地震資料的重要補(bǔ)充部分，提供了真實(shí)的高頻反射系數(shù)。引入測井中的高頻信息，確保了高頻信息的真實(shí)性，改進(jìn)常規(guī)反褶積方法，具體原理如下。常規(guī)的雙邊最小平方反褶積方程， 根據(jù)誤差平方和最小的原則， 將子波 w( t) 壓縮為尖脈沖，求取反子波： （11）式中，m為采樣點(diǎn)數(shù)，子波和反子波。引入井中反射系數(shù)r(t)對(duì)方法進(jìn)行改進(jìn)。根據(jù)誤差平方和最小的原則， 將地震記錄x(t)壓縮為反射系數(shù)r(t)， 求取反子波a(t)。上式變?yōu)椋?（12）式中，為地震記錄，為反子波。反射系數(shù)為。由此可推導(dǎo)

15、出反子波的計(jì)算公式如下： （13）用反褶積算子與地震記錄褶積從而得到反褶積輸出結(jié)果。由圖7看出反褶積后過井剖面分辨率明顯提高，信息更加豐富、細(xì)節(jié)更加突出，從頻譜上看，主頻提高，頻帶變寬。 （a）反褶積前過井剖面及頻譜 （b）反褶積后過井剖面及頻譜圖7 反褶積前后過井剖面及頻譜國外文獻(xiàn)調(diào)研：near-surface high-frequency absorption compensation in alluvial plainqinghui cui*,yongjun rui ,and xinmin shang, shengli gri sinopec（2011 seg年會(huì)）由松散沉積的淤泥形成的

16、沖積平原，其近地表地震波能量衰減嚴(yán)重，特別是它的高頻部分，高頻能量補(bǔ)償成為提高地震資料分辨率的重要的方法。為了消除近地表地震波衰減，提出了一種粘彈性波動(dòng)方程延拓。首先，我們建立一個(gè)近地表數(shù)學(xué)模型，再加上近地表q值和速度模型，然后利用粘彈性波方程延拓，把接收的地震波信號(hào)延伸到高速層頂部。第一步是從微測井?dāng)?shù)據(jù)中計(jì)算近地表品質(zhì)因子q和p波速度v，這個(gè)決定了近地表數(shù)學(xué)模型以及q和v之間的關(guān)系，然后利用微測井約束層析方法計(jì)算近地表全局區(qū)域p波的速度，最后用相移加插值算法進(jìn)行波場延拓。流程圖如下：圖1 基于粘彈性波動(dòng)方程延拓的近地表高頻衰減補(bǔ)償方法流程圖kjartansson提出的粘彈性介質(zhì)模型指出，q因

17、子與頻率的關(guān)系，當(dāng)頻率小于某一值時(shí)，速度和q因子可以描述為： （1）目前，主要有兩種近地表q反演的方法，一個(gè)是利用微測井?dāng)?shù)據(jù)在井位置處進(jìn)行q反演，然后利用空間解釋求取整個(gè)區(qū)域的q值，另一種方法是利用地震反射記錄的初至波直接反演。前者不精確，后者不夠穩(wěn)定，兩種方法都不太能得到令人滿意的結(jié)果。品質(zhì)因子與速度之間有非常緊密的聯(lián)系，可以通過分析所有的微測井?dāng)?shù)據(jù)對(duì)整個(gè)區(qū)域求取一個(gè)統(tǒng)一的公式：， （2） 其中a,b,c,d是常數(shù)。頻移加插值算法用于波動(dòng)方程延拓，該算法是一個(gè)改進(jìn)的頻移方法，能克服水平層狀介質(zhì)的假設(shè)，適于處理橫向變速的介質(zhì)。這個(gè)算法是在頻率空間域?qū)崿F(xiàn)的，可以通過頻移方法，利用不同的速度，分別

18、計(jì)算一系列的延拓波場，公式如下： （3） 其中， （4） ， （5）通過橫向上的反傅里葉變換，波場可以被轉(zhuǎn)換到頻率空間域。利用拉格朗日插值方法計(jì)算，用一系列速度計(jì)算延拓波場。經(jīng)過這些步驟，就可以逐漸延伸到高速度層頂部。由于近地表強(qiáng)吸收性質(zhì)，經(jīng)過長時(shí)間的地震波衰減過程，高頻成分振幅遠(yuǎn)遠(yuǎn)低于周圍噪聲的水平，幾乎為零，這樣會(huì)導(dǎo)致高頻能量補(bǔ)償時(shí)計(jì)算不穩(wěn)定。為了克服這個(gè)問題，高頻能量補(bǔ)償不能通過實(shí)際衰減，而應(yīng)該利用實(shí)際記錄的背景噪聲能量水平，所以每個(gè)延拓層的高頻能量補(bǔ)償穩(wěn)定條件變成了下式： （6）建立近地表速度和品質(zhì)因子模型，通過速度聲波方程正演模型得到無衰減地震記錄，用粘彈性波動(dòng)方程正演模型得到衰減地

19、震記錄，然后進(jìn)行補(bǔ)償。遠(yuǎn)偏移距道集的補(bǔ)償后的結(jié)果非常接近聲波記錄，說明這種補(bǔ)償算法是可行的。某一區(qū)域主要是砂泥巖為主，地震波的能量吸收明顯，主頻低，頻帶窄，利用波場延拓方法能量補(bǔ)償之后結(jié)果得到了明顯的改進(jìn)，地震數(shù)據(jù)的分辨率得到了明顯的提高，特別是重疊的波形被分離開來，弱相位得到了加強(qiáng)，頻率的有效帶寬得到了擴(kuò)展。a:聲波模擬記錄 b:粘彈性波模擬記錄 c:補(bǔ)償后的記錄圖2 模型a:原始剖面及其頻譜 b:補(bǔ)償后的剖面及其頻譜圖3 實(shí)例辦公室做的工作：1 利用合成地震記錄校正井旁道及道集地震資料與測井?dāng)?shù)據(jù)之間的匹配問題已成為地震巖性解釋、油藏描述、儲(chǔ)層反演以及儲(chǔ)層監(jiān)測過程中的主要問題。目前，用于地震

20、與測井匹配校正的方法主要有經(jīng)典的閉合差校正法、基于楊氏諧振q模型的直接匹配法以及基于多分辨分析技術(shù)的匹配方法等。這種匹配過程可以通過兩種方式實(shí)現(xiàn)：以井旁實(shí)際地震記錄為主，使合成地震記錄逼近井旁實(shí)際地震道；以合成地震記錄為主，使井旁地震道逼近合成地震記錄。通常我們所說的匹配校正應(yīng)當(dāng)是第一種方式，但有時(shí)也采用第二種方式對(duì)地面地震數(shù)據(jù)進(jìn)行時(shí)移和相位校正。我們利用第二種方式進(jìn)行井控能量補(bǔ)償。利用合成記錄校正井旁道及道集的具體思路如下：1）利用井?dāng)?shù)據(jù)制作反射系數(shù)時(shí)間序列或者合成地震記錄；2）求反射系數(shù)序列或者合成記錄整體振幅包絡(luò)線a；3）對(duì)井旁自激自收地震記錄進(jìn)行能量補(bǔ)償，假設(shè)補(bǔ)償參數(shù)向量為eta(t)，然后求解補(bǔ)償后記錄整體振幅包絡(luò)線b(eta)；4）匹配a與b(eta)，目標(biāo)函數(shù)設(shè)為 （1）此過程可以通過反演的方法實(shí)現(xiàn)，得到最佳匹配參數(shù)eta(t)；5）將依據(jù)道間差異調(diào)整eta(t)并依次向鄰道應(yīng)用。利用合成地震記錄校正井旁道及道集存在著一定的問題。在水平層狀介質(zhì)假設(shè)條件下自激自收道與井合成記錄是時(shí)間對(duì)應(yīng)的，而實(shí)際情況，水平層狀介質(zhì)假設(shè)不成