宋吉杰,禹金營(yíng),王 成,張 猛

(1.中國(guó)石油化工股份有限公司油田勘探開發(fā)事業(yè)部,北京100728；2.中國(guó)石油化工股份有限公司西北分公司研究院,新疆烏魯木齊830000；3.中國(guó)石油天然氣股份有限公司大慶勘探開發(fā)研究院,黑龍江大慶163712;4.中國(guó)石油化工股份有限公司勝利油田物探研究院,山東東營(yíng)257100)

決定地震資料分辨率的關(guān)鍵是反射信號(hào)的有效頻帶寬度,而制約有效頻帶寬度的兩個(gè)主要因素是激發(fā)信號(hào)的能量(頻帶寬度)和信號(hào)傳播過(guò)程中的能量衰減。在能量衰減中,近地表吸收衰減的影響尤為嚴(yán)重。地層品質(zhì)因子Q是表征介質(zhì)中地震波衰減強(qiáng)弱的重要參數(shù)。在實(shí)際生產(chǎn)過(guò)程中,一般需要對(duì)地震資料進(jìn)行有效的Q補(bǔ)償處理,通過(guò)補(bǔ)償?shù)卣鸩ㄔ诘叵陆橘|(zhì)傳播過(guò)程中的衰減,也是恢復(fù)地下地層反射系數(shù)的過(guò)程。高效Q補(bǔ)償處理離不開高精度的Q值估計(jì)。

為了描述大地吸收效應(yīng),1962年FUTTERMAN[1]詳細(xì)論述了巖石對(duì)地震波的吸收衰減是地層的基本特性,并基于Q與頻率無(wú)關(guān)的假設(shè),推導(dǎo)出振幅吸收與相速度的頻散公式。HALE[2]在FUTTERMAN模型的基礎(chǔ)上提出了一種反Q濾波方法,該濾波函數(shù)的相位函數(shù)由對(duì)數(shù)振幅譜的希爾伯特變換求得;KJARTANSSON[3]基于上升時(shí)間成像原理估算地層的Q值;SCHOENBERGER等[4]利用兩口井的VSP合成記錄計(jì)算了反射波的衰減;HARGREAVES等[5]提出的反Q濾波方法都只針對(duì)相位補(bǔ)償,不考慮振幅的損失;WANG[6]基于KOLSKY模型,從一維波動(dòng)方程出發(fā),推導(dǎo)了能量補(bǔ)償和相位校正的公式,提出了一種穩(wěn)定有效的頻率域反Q濾波方法;TANER等[7]對(duì)7種常見(jiàn)Q值估算方法進(jìn)行了對(duì)比分析。

隨著油氣勘探對(duì)地震分辨率需求的不斷提高,近地表吸收衰減補(bǔ)償問(wèn)題得到重視。裴江云等[8]利用面波衰減的Q模型估算近地表Q值,以消除近地表不一致性影響,恢復(fù)地震波相對(duì)振幅關(guān)系,但并未考慮相位校正的問(wèn)題。尹喜玲[9]等用線性過(guò)渡體理論初步討論了沙漠地區(qū)近地表的吸收規(guī)律,認(rèn)為近地表處處存在速度界面,多次反射和透射是造成衰減的主要原因。于承業(yè)等[10]、張文等[11]提出利用雙井微測(cè)井資料通過(guò)頻移法估算近地表Q值,強(qiáng)調(diào)了近地表衰減補(bǔ)償能夠在時(shí)間、空間和頻率三域內(nèi)有效地消除近地表影響。

因地震波衰減機(jī)理極為復(fù)雜,故衰減參數(shù)估算的精度有限,造成相關(guān)應(yīng)用的發(fā)展較為緩慢。常規(guī)近地表衰減補(bǔ)償主要通過(guò)地表一致性振幅補(bǔ)償聯(lián)合反褶積處理來(lái)實(shí)現(xiàn),難以恢復(fù)近地表空變吸收對(duì)子波能量和相位的改造,甚至可能會(huì)引起油氣儲(chǔ)層信息的畸變,給儲(chǔ)層解釋帶來(lái)假象或困難。本文在前人研究的基礎(chǔ)上,開展了近地表介質(zhì)Q估計(jì)方法、衰減模型構(gòu)建以及衰減補(bǔ)償處理技術(shù)的研究并將其應(yīng)用于塔河北部油田的勘探開發(fā)。

1 近地表介質(zhì)Q估計(jì)方法

在綜合考慮現(xiàn)有微測(cè)井Q估計(jì)及大炮初至波Q估計(jì)方法優(yōu)缺點(diǎn)的基礎(chǔ)上,提出了一種新的多井微測(cè)井與地面聯(lián)合觀測(cè)系統(tǒng),實(shí)現(xiàn)控制點(diǎn)的精細(xì)Q估計(jì)。同時(shí),通過(guò)統(tǒng)計(jì)炮檢點(diǎn)相對(duì)衰減系數(shù)來(lái)反演相對(duì)Q值,進(jìn)一步利用不同尺度的信息融合技術(shù)將兩種Q值相融合,構(gòu)建近地表吸收衰減模型,通過(guò)穩(wěn)定反Q濾波實(shí)現(xiàn)近地表衰減補(bǔ)償。具體技術(shù)流程見(jiàn)圖1。

1.1 控制點(diǎn)精細(xì)Q估計(jì)

1.1.1 折射波法調(diào)查和透射波法調(diào)查

折射波法調(diào)查采用的觀測(cè)系統(tǒng)如圖2所示,左側(cè)設(shè)置了6口激發(fā)井(炸藥激發(fā)),右側(cè)為等道間距的接收排列(單檢波器接收),中間道位置布設(shè)一口打穿高速層的井,并設(shè)置井底道和對(duì)應(yīng)的井口道接收,其基本原理是通過(guò)截取井底道直達(dá)波和井口道折射波的初至信號(hào),基于質(zhì)心頻移法求取出表層的品質(zhì)因子Q。該調(diào)查數(shù)據(jù)可選取多炮共道(井底道)觀測(cè)記錄，或共炮多道(井口道)觀測(cè)記錄,基于信號(hào)隨炮檢距的變化完成高速層品質(zhì)因子Q的吸收衰減調(diào)查研究。

圖1 近地表介質(zhì)吸收衰減補(bǔ)償技術(shù)實(shí)現(xiàn)流程

圖2 折射波法調(diào)查觀測(cè)系統(tǒng)示意

透射波法調(diào)查采用了多井微測(cè)井觀測(cè)系統(tǒng),如圖3 所示,與常規(guī)的雙井微測(cè)井[12-13]不同之處在于多出了4口低速層井中道接收井。激發(fā)方式采用雷管。工作原理是在激發(fā)井中高速層內(nèi)激發(fā),鄰近接收井中的高速層井中道和地面井口道分別獲得高速層直達(dá)地震波和透射波,兩者初至波形的差異實(shí)際上就是由近地表的吸收衰減引起的。與上述折射波法計(jì)算出的高速層品質(zhì)因子Q相結(jié)合,即可計(jì)算出低速層的等效品質(zhì)因子Q。對(duì)于低速層厚度不太大的地區(qū),該方法可同時(shí)完成表層結(jié)構(gòu)調(diào)查。井深以打穿高速層為宜(塔河地區(qū)采用20m,因?yàn)樵搮^(qū)低速層厚度小于10m),接收井距離激發(fā)井4m,其它4口井間距2m。在離激發(fā)井最近的一口接收井內(nèi),距離高速層頂面2m處設(shè)置一個(gè)檢波器。為了更精細(xì)地研究低速層吸收情況,圖3右側(cè)設(shè)置了不同深度的低速層井中道,并在井底各插入一個(gè)檢波器,其埋深可以參照表層結(jié)構(gòu)調(diào)查結(jié)果。

為了取得可靠的吸收衰減調(diào)查記錄,要確保井中道的檢波器不漏電,與井底地層耦合良好;放炮前要做背景檢測(cè),擇“靜”采集。記錄評(píng)價(jià)的準(zhǔn)則:井中道地震波振幅遠(yuǎn)強(qiáng)于地面井口道,主頻遠(yuǎn)高于地面井口道,初至區(qū)無(wú)干擾。

上述兩種方法所求取的Q值可以相互印證,主要差異在于激發(fā)子波的不同(炸藥和雷管)。

圖3 透射波法(多井微測(cè)井)調(diào)查觀測(cè)系統(tǒng)示意

1.1.2 質(zhì)心頻移法Q估計(jì)

質(zhì)心頻移方法[14-15]是基于地震子波傳播中出現(xiàn)的子波質(zhì)心頻率向低頻端移動(dòng)這一特征進(jìn)行Q值估算的方法。振幅比法及譜比法等常用方法主要利用地震波振幅信息對(duì)Q值進(jìn)行估算。振幅受幾何擴(kuò)散、反射、透射等效應(yīng)的影響較大,僅由振幅衰減難以準(zhǔn)確估算非彈性介質(zhì)對(duì)地震波的吸收情況。而子波振幅譜質(zhì)心頻率下降或波形加寬等衰減特征受上述因素影響較小[16],因此,質(zhì)心頻移法較振幅比法和譜比法可以獲得更加穩(wěn)定可靠的Q估計(jì)。

共炮多道(地面道)或多炮共道(井底道)記錄的初至波可用來(lái)估計(jì)高速層Q值,通常采用不同炮檢距的兩道初至波,基于質(zhì)心頻移法求取Q值。

高速層中接收的直達(dá)波子波作為參考子波,令參考子波振幅譜為S(f),則在低速層中接收子波的振幅譜R(f)表示為:

(1)

式中:G(f)為幾何擴(kuò)散、反射、透射傳播過(guò)程中相位累積及儀器響應(yīng)等效應(yīng);H(f)描述非彈性介質(zhì)對(duì)振幅的吸收效應(yīng)。令入射波振幅譜S(f)為Gaussian型展布,即:

(2)

沿傳播路徑對(duì)Q值逐層進(jìn)行估算[7]:

(5)

對(duì)于品質(zhì)因子為Q,旅行時(shí)為t的表層而言,(5)式可寫為:

(6)

1.2 相對(duì)Q值估計(jì)

疊前Q補(bǔ)償處理通常需要輸入各檢波點(diǎn)的近地表Q值,即需要建立表層Q模型,而利用微測(cè)井?dāng)?shù)據(jù)只能得到調(diào)查點(diǎn)的Q值,而工區(qū)內(nèi)各檢波點(diǎn)的Q值只能通過(guò)地面地震數(shù)據(jù)來(lái)求取。地震初至波未經(jīng)地下界面反射,由炮點(diǎn)激發(fā)經(jīng)表層到達(dá)檢波點(diǎn),其橫向振幅、頻率和波形差異主要由表層變化引起。我們提出了一種穩(wěn)定統(tǒng)計(jì)性反演估計(jì)方法,具體如下:先利用地表一致性振幅補(bǔ)償模塊,再通過(guò)高斯塞德?tīng)柕鷮⒊踔聊芰康牟町惙纸鉃榕邳c(diǎn)、檢波點(diǎn)和偏移距三項(xiàng)差異,計(jì)算炮點(diǎn)、檢波點(diǎn)相對(duì)衰減系數(shù)r,最后利用譜比法通過(guò)相對(duì)衰減系數(shù)r、參考頻率和表層旅行時(shí)t求得相對(duì)Q值。

1.3 信息融合與近地表Q模型構(gòu)建

計(jì)算得到了野外調(diào)查估計(jì)的Q值和大炮統(tǒng)計(jì)反演的Q值,前者是精細(xì)的、低波數(shù)的和絕對(duì)的;后者是統(tǒng)計(jì)的、高波數(shù)的和相對(duì)的。利用信息融合策略,對(duì)不同來(lái)源、不同尺度和不同精度的數(shù)據(jù)進(jìn)行融合,起到絕對(duì)Q約束相對(duì)Q、高波數(shù)Q補(bǔ)充低波數(shù)Q的作用,并提高最終結(jié)果的精度和可靠性,由此建立起一個(gè)精確的近地表Q模型。

本著先標(biāo)定、后插值、再計(jì)算的原則,該方法具體實(shí)現(xiàn)的3步解決方案如下:①基于吸收理論譜比法求取調(diào)查點(diǎn)上的標(biāo)定因子;②采用協(xié)克里金插值方法將標(biāo)定因子內(nèi)插到接收點(diǎn);③再由譜比法計(jì)算每個(gè)接收點(diǎn)Q值。上述3步法達(dá)到了信息融合的目的,也實(shí)現(xiàn)了近地表Q模型的構(gòu)建。

1.4 近地表Q補(bǔ)償處理方法

近地表Q模型建立之后,就可以對(duì)疊前數(shù)據(jù)在頻率域進(jìn)行近地表Q補(bǔ)償處理,通過(guò)振幅補(bǔ)償項(xiàng)和相位補(bǔ)償項(xiàng)構(gòu)成完整的反Q濾波,由波動(dòng)方程推導(dǎo)的補(bǔ)償算式如下:

(7)

式中:U(τ,ω)為未經(jīng)Q補(bǔ)償?shù)念l率域數(shù)據(jù),即對(duì)輸入地震道數(shù)據(jù)進(jìn)行傅氏變換的結(jié)果;U(τ+Δτ,ω)為Q補(bǔ)償后的頻率域數(shù)據(jù);τ為傳播時(shí)間;△τ為表層旅行時(shí);ω為角頻率,是地震頻帶內(nèi)與最高頻率有關(guān)的調(diào)諧頻率。

近地表Q補(bǔ)償輸入的是疊前數(shù)據(jù),實(shí)現(xiàn)上有著比深層Q補(bǔ)償更大的難度,需要在提高分辨率的同時(shí)兼顧信噪比原則。因此,基于穩(wěn)定的策略,本文引用了WANG[6]提出的穩(wěn)定反Q濾波方法。WANG對(duì)振幅項(xiàng)做了如下改進(jìn):

(8)

式中:Λ(ω)為穩(wěn)定的振幅補(bǔ)償量;σ2是穩(wěn)定因子;β為經(jīng)驗(yàn)穩(wěn)定公式。

(9)

式中:Glim為增益限制,單位dB,是一個(gè)可調(diào)整的參數(shù);ωh為中心頻率(角頻率),也是一個(gè)可調(diào)整的參數(shù),跟地震波頻帶的最高頻率有關(guān)系;γ=(2/π)·tan-1(1/2Q);t為表層旅行時(shí)。

上述反Q濾波算法穩(wěn)定性體現(xiàn)在可以通過(guò)增益的門檻值控制補(bǔ)償?shù)念l帶范圍,防止對(duì)高頻噪聲的過(guò)分補(bǔ)償。由于補(bǔ)償量是隨頻率和時(shí)間的增大逐漸增多,通常的補(bǔ)償算法都會(huì)對(duì)高頻端過(guò)分補(bǔ)償,造成噪聲過(guò)量、信號(hào)失真,而WANG[6]采用的策略是對(duì)吸收嚴(yán)重的頻率成分不再進(jìn)行補(bǔ)償,這樣可以避免對(duì)高頻噪聲的過(guò)量補(bǔ)償。該算法的另一優(yōu)點(diǎn)是補(bǔ)償振幅的同時(shí)調(diào)整相位,這解決了近地表吸收衰減造成的地震波能量損失和頻散問(wèn)題。

算法中可調(diào)整的參數(shù)如增益限制Glim、中心頻率ωh,對(duì)補(bǔ)償也起著重要作用。在實(shí)際補(bǔ)償工作中需要根據(jù)數(shù)據(jù)的具體情況及試驗(yàn)效果選擇合理的參數(shù)。增益限制Glim的選擇從兩方面入手,一方面依據(jù)期望補(bǔ)償后頻寬的指標(biāo)來(lái)估算,另一方面依據(jù)補(bǔ)償后數(shù)據(jù)的頻譜和信噪比來(lái)定,過(guò)小則補(bǔ)償效果不理想,過(guò)大則會(huì)導(dǎo)致高頻噪聲放大、信噪比降低;中心頻率ωh通常要稍高于實(shí)際地震資料有效頻帶的上限值(以dB譜-20Hz上限值為參考),保證吸收補(bǔ)償充分作用于有效頻帶內(nèi)。

2 實(shí)際應(yīng)用及效果分析

研究區(qū)位于塔河北部,面積約為300km2。區(qū)內(nèi)地表主要有浮土、農(nóng)田與河灘。近地表低速帶沉積松散,干燥疏松,厚度變化范圍6～10m,整體速度較低,平均約400m/s。表層介質(zhì)對(duì)地震波的高頻信號(hào)有強(qiáng)烈的垂向吸收衰減作用,導(dǎo)致地下反射信號(hào)的有效頻帶變窄,如圖4所示,地面道接收的地震波相較于井中道有強(qiáng)衰減現(xiàn)象;并且表層橫向吸收衰減的差異很容易造成地震道間反射信號(hào)的振幅差異和相位不一致。該區(qū)儲(chǔ)層是白堊系的碎屑巖和奧陶系的碳酸鹽巖,因高頻信息的保真成像關(guān)乎薄層砂體的刻畫和縫洞串珠成像的可靠性,故上述近地表吸收衰減問(wèn)題將直接影響有效儲(chǔ)層預(yù)測(cè)。也就是說(shuō),近地表吸收衰減是該研究區(qū)勘探開發(fā)的瓶頸。該區(qū)以往主要通過(guò)地表一致性振幅補(bǔ)償和反褶積方法進(jìn)行地表一致性處理,但這種處理對(duì)于吸收補(bǔ)償作用不明顯。采用本文提出的基于信息融合的近地表介質(zhì)Q估計(jì)方法及穩(wěn)定的反Q濾波處理技術(shù),利用調(diào)查點(diǎn)精細(xì)Q值約束相對(duì)Q值整合建立的近地表Q模型,通過(guò)反Q濾波實(shí)現(xiàn)近地表衰減補(bǔ)償,提高了地震資料的分辨率。

2.1 吸收衰減調(diào)查記錄

為了對(duì)近地表結(jié)構(gòu)及衰減特性進(jìn)行詳細(xì)調(diào)查,本區(qū)共設(shè)計(jì)完成近地表結(jié)構(gòu)調(diào)查點(diǎn)484個(gè),吸收衰減調(diào)查點(diǎn)6個(gè)。圖4對(duì)比了某衰減調(diào)查點(diǎn)折射波法采集的記錄(圖4a)和透射波法采集的記錄(圖4b),可以看出,高速層井中道(藍(lán)色波形)初至清晰干脆,振幅強(qiáng)、頻率高(該特征也作為評(píng)價(jià)吸收衰減調(diào)查資料是否合格的重要依據(jù)),高速層井中道和地面道初至信號(hào)之間、高速層井中道和其右側(cè)相鄰的低速層的井中道初至信號(hào)之間均構(gòu)成較大的反差。圖4中地面道之間隨炮檢距的增大初至波的衰減規(guī)律性越發(fā)明顯。

2.2 近地表Q模型

表1為基于6個(gè)吸收衰減調(diào)查點(diǎn)所得的井中與地面記錄,采用質(zhì)心頻率法估算得到的近地表低速層Q值?？梢钥闯?近地表Q值整體偏小,變化范圍為2.0～3.5。

圖4 折射波法采集的地震記錄(a)和透射波法采集的地震記錄(b)對(duì)比

表1 6個(gè)吸收衰減規(guī)律調(diào)查點(diǎn)的Q值

圖5是基于表層結(jié)構(gòu)調(diào)查數(shù)據(jù)建立表層模型求取的接收點(diǎn)表層旅行時(shí)t的平面分布圖(圖5a)和信息融合后接收點(diǎn)近地表Q值分布圖(圖5b)。可以看出,所求Q值與表層旅行時(shí)t有較好的反比對(duì)應(yīng)關(guān)系,即表層旅行時(shí)t小,相對(duì)吸收小,Q值大,這也表明該區(qū)所求Q值符合客觀衰減規(guī)律。

2.3 近地表Q補(bǔ)償效果

圖6對(duì)比了近地表吸收衰減補(bǔ)償前后的單炮記錄和相應(yīng)目的層段反射波的振幅譜?？梢钥闯?通過(guò)近地表吸收衰減補(bǔ)償,單炮記錄波組特征橫向連續(xù)性得到改善,目的層縱向分辨率得到較大提高,在不降低信噪比的前提下,有效頻帶得到了合理拓寬。進(jìn)一步對(duì)比反射波的振幅譜可見(jiàn),補(bǔ)償后的反射波頻帶展寬超過(guò)20Hz。

圖7對(duì)比了兩種處理方法得到地震疊加成像剖面,圖7a是采用地表一致性振幅補(bǔ)償和反褶積聯(lián)合處理得到的疊加剖面;圖7b是近地表Q補(bǔ)償處理得到的剖面?？梢钥闯?Q補(bǔ)償后地震剖面分辨率顯著提高,圖7b中目的層的弱反射連續(xù)性得到加強(qiáng)。

圖5 檢波點(diǎn)表層旅行時(shí)(a)與信息融合后檢波點(diǎn)表層Q值(b)的平面分布

圖6 Q值補(bǔ)償前(a)、后(b)的單炮記錄對(duì)比及Q值補(bǔ)償前(c)、后(d)的目的層反射波振幅譜對(duì)比

圖7 地表一致性振幅補(bǔ)償和反褶積聯(lián)合處理剖面(a)和近地表Q補(bǔ)償處理剖面(b)的對(duì)比

3 結(jié)論

本文提出一種基于信息融合的近地表介質(zhì)Q估計(jì)方法以及穩(wěn)定的反Q濾波處理技術(shù),在此過(guò)程中得到以下結(jié)論和認(rèn)識(shí):

1) 基于多井微測(cè)井與地面聯(lián)合觀測(cè)記錄,采用質(zhì)心頻率法可以實(shí)現(xiàn)調(diào)查點(diǎn)近地表Q值的精細(xì)估計(jì),不失為一種近地表衰減調(diào)查或Q估計(jì)的有效方法技術(shù)。

2) 基于地面地震記錄通過(guò)穩(wěn)定統(tǒng)計(jì)性反演可以獲得近地表相對(duì)Q值平面分布。

3) 利用信息融合策略,本著先標(biāo)定、后插值、再計(jì)算的3步法方案,實(shí)現(xiàn)調(diào)查點(diǎn)近地表Q值與相對(duì)Q值的有機(jī)融合,構(gòu)建精細(xì)的近地表衰減模型,為進(jìn)一步衰減補(bǔ)償處理奠定基礎(chǔ)。

4) 塔河北部油田實(shí)際應(yīng)用結(jié)果表明采用不同尺度Q值融合的近地表衰減模型構(gòu)建方法,并借助于改進(jìn)的穩(wěn)定反Q濾波方法可使目的層反射波頻帶展寬超過(guò)20Hz,極大地提高了地震資料分辨率。

目前該方法僅基于近地表為層狀結(jié)構(gòu)的假設(shè),對(duì)于非層狀結(jié)構(gòu)的情況還需要繼續(xù)探索,以增強(qiáng)方法的適應(yīng)性。

參 考 文 獻(xiàn)

[1] FUTTERMAN W I.Dispersive body waves[J].Journal of Geophysical Research,1962,67(13):5279-5291

[2] HALE D.An inverseQ-filter[J].SEP Report 26,1981:231-243

[3] KJARTANSSON E.ConstantQ-wave propagation and attenuation[J].Journal of Geophysical Research (Solid Earth),1979,84(B9):4737-4748

[4] SCHOENBERGER M,LEVIN FK.Apparent attenuation due to intrabed multiples[J].Geophysics,1974,39(3):278-291

[5] HARGREAVES N D,CALVERT A J.InverseQfiltering by Fourier transform[J].Geophysics,1991,56(4):519-527

[6] WANG Y H.A stable and efficient approach of inverseQfiltering[J].Geophysics,2002,67(2):657-663

[7] TANER M T,TREITEL S T.A Robust method forQestimation[J].Expanded Abstracts of 73rdAnnual Internat Mtg,2003:710-713

[8] 裴江云,陳樹民,劉振寬,等.近地表Q值求取及振幅補(bǔ)償[J].地球物理學(xué)進(jìn)展,2001,16(4):18-22

PEI J Y,CHEN S M,LIU Z K,et al.Near surfaceQvalue extraction and amplitude compensation[J].Progress in Geophysics,2001,16(4):18-22

[9] 尹喜玲,石戰(zhàn)結(jié),田鋼.近地表低降速帶地震波傳播規(guī)律初探[J].地球物理學(xué)進(jìn)展,2009,24(2):398-407

YIN X L,SHI Z G,TIAN G.Study on transmittable laws of seismic waves in unconsolidated layers[J].Progress in Geophysics,2009,24(2):398-407

[10] 于承業(yè),周志才.利用雙井微測(cè)井資料估算近地表Q值[J].石油地球物理勘探,2011,46(1):89-92

YU C Y,ZHOU Z C.Estimation of near surfaceQvalue based on the datasets of the uphole survey in double hole[J].Oil Geophysical Prospecting,2011,46(1):89-92

[11] 張文,周志才,于承業(yè).一種改進(jìn)的近地表強(qiáng)吸收補(bǔ)償方法研究[J].石油物探,2017,56(2):210-215

ZHANG W,ZHOU Z C,YU C Y.An improved compensation method for near surface strong absorption[J].Geophysical Prospecting for Petroleum,2017,56(2):210-215

[12] 李天樹,陳寶德,蘇德仁.雙井微測(cè)井技術(shù)在表層結(jié)構(gòu)調(diào)查中的應(yīng)用[J].石油物探,2004,43(5):471-474

LI T S,CHEN B D,SU D R.Application of twin-well microlog in near surface investigation[J].Geophysical Prospecting for Petroleum,2004,43(5):471-474

[13] 王永卓,杜桂鋒,傅朝奎,等.虛反射面對(duì)地震激發(fā)效果的影響分析及其應(yīng)用[J].地球物理學(xué)進(jìn)展,2009,24(4):1454-1460

WANG Y Z,DU G F,FU Z K,et al.Influence analysis of ghosting interface on seismic shooting effect and its application[J].Progress in Geophysics,2009,24(4):1454-1460

[14] QUAN Y L,HARRIS J M.Seismic attenuation tomography using the frequency shift method[J].Geophysics,1997,62(3):895-905

[15] 宮同舉,孫成禹,彭洪超,等.幾種提取品質(zhì)因子方法的對(duì)比分析[J].勘探地球物理進(jìn)展,2009,32(4):252-256

GONG T J,SUN C Y,PENG H C,et al.Comparison and analysis of several methods for extracting quality factors[J].Progress in Exploration Geophysics,2009,32(4):252-256

[16] YANG S L,Gao J H.Seismic quality factors estimation from spectral correlation[J].IEEE Geoscience and Remote Sensing Letters,2008,5(4):740-744

[17] 李偉娜,云美厚,黨鵬飛,等.基于微測(cè)井資料的雙線性回歸穩(wěn)定Q估計(jì)[J].石油物探,2017,56(4):483-490

LI W N,YUN M H,DANG P F,et al.StabilityQestimation by dual linear regression based on uphole survey data[J].Geophysical Prospecting for Petroleum,2017,56(4):483-490