關(guān)鍵詞：近地表，品質(zhì)因子Q，Q值估算，表層Q場(chǎng)，四川盆地

0 引言

地震波在地下介質(zhì)中傳播，能量會(huì)被地層吸收，造成地震波的振幅衰減、相位畸變，導(dǎo)致地震資料分辨率降低。研究[1-2]表明：與深部地層介質(zhì)的吸收衰減效應(yīng)相比，近地表介質(zhì)對(duì)地震波高頻成分的吸收約占地層吸收總量的80%[3]，由于表層巖性、速度、厚度的橫向變化，吸收衰減程度會(huì)有所不同，還會(huì)造成地震波道間能量、頻率和相位的不一致，影響同相疊加結(jié)果的橫向分辨率。因此，研究地下介質(zhì)對(duì)地震波的吸收特性，對(duì)于改善地震資料分辨率、油氣地球物理屬性反演和儲(chǔ)層描述等有著重要的意義[4-5]。

地層品質(zhì)因子Q 值可定量刻畫地震波在黏彈性介質(zhì)中的能量吸收和速度頻散[6]，準(zhǔn)確求取近地表Q值再進(jìn)行反Q 濾波補(bǔ)償可有效消除表層地震波吸收衰減的影響，實(shí)現(xiàn)波形保真保幅并提高分辨率的效果。其中的關(guān)鍵是如何獲得準(zhǔn)確的近地表Q 值和Q 模型。

對(duì)于表層Q 值求取及建模方法，很多學(xué)者進(jìn)行了研究。于承業(yè)等[7]利用雙井微測(cè)井資料，通過建立和求解地面和井下檢波器的峰值頻率方程組，求得近地表Q 值的解析解；裴江云等[8]利用地震數(shù)據(jù)中不同類型的波場(chǎng)進(jìn)行近地表Q 反演；李國(guó)發(fā)等[9]采用井—地聯(lián)合觀測(cè)方式就近地表吸收結(jié)構(gòu)進(jìn)行了實(shí)驗(yàn)分析，提出了不受激發(fā)因素影響的吸收參數(shù)層析反演方法；肖永新等[10]提出了利用雙井微測(cè)井高速層中多個(gè)炮點(diǎn)的數(shù)據(jù)，通過擬合方法得到精確的近地表Q 值；李偉娜等[11]基于微測(cè)井資料發(fā)展了雙線性回歸穩(wěn)定Q估計(jì)方法，取得了較好的效果；李擁軍等[12]提出了一種復(fù)數(shù)域快速匹配追蹤分解并結(jié)合對(duì)數(shù)譜比法估算微測(cè)井?dāng)?shù)據(jù)近地表Q 值的方法；許李囡等[13]提出了一種基于S 變換和變分法提取品質(zhì)因子的方法，克服了Q值計(jì)算對(duì)地震子波類型及窗函數(shù)長(zhǎng)度的依賴。以上均是基于Q 值求取進(jìn)行的研究。宋吉杰等[14]提出基于信息融合的近地表介質(zhì)Q 估計(jì)方法，利用多井微測(cè)井和地面聯(lián)合的觀測(cè)系統(tǒng)，通過統(tǒng)計(jì)反演估算全區(qū)Q值[15]；蘇勤等[16]引入表層相對(duì)衰減系數(shù)的概念，在共炮檢域迭代計(jì)算相對(duì)衰減系數(shù)，進(jìn)而求取表層Q 模型；馬凱等[17]采用雙井微測(cè)井獲取的絕對(duì)Q 值和頻移法計(jì)算的相對(duì)Q 值，建立了準(zhǔn)確的Q 場(chǎng)并進(jìn)行Q 補(bǔ)償，獲得了良好的效果。以上方法都利用了微測(cè)井和地面地震數(shù)據(jù)，但都是先得到相對(duì)Q 值，通過標(biāo)定和插值得到最終表層Q 模型，該表層Q 模型具有橫向空間變化，但缺少縱向時(shí)間變化。

為提高近地表Q 建模精度，本文針對(duì)四川盆地近地表厚度較薄的特點(diǎn)，提出了一種結(jié)合單、雙微測(cè)井及地面地震數(shù)據(jù)信息進(jìn)行近地表Q 建模的方法。首先，對(duì)于單井微測(cè)井?dāng)?shù)據(jù)，本文提出了一種改進(jìn)的譜比法進(jìn)行Q 值計(jì)算，通過構(gòu)造統(tǒng)計(jì)子波譜，消除子波效應(yīng)，提高了計(jì)算結(jié)果的穩(wěn)定性；然后擬合Q 與速度v 的關(guān)系，將層析反演得到的近地表速度場(chǎng)轉(zhuǎn)換成Q 場(chǎng)，從而得到一個(gè)時(shí)空變的Q 場(chǎng)。由于層析反演的近地表速度場(chǎng)深度采樣間隔為10 m，難以反映近地表幾米厚的表層變化。為此結(jié)合雙井微測(cè)井的多炮擬合算法，計(jì)算表層絕對(duì)Q 值，再結(jié)合表層速度模型和地面地震振幅及走時(shí)信息計(jì)算相對(duì)Q 場(chǎng)，用絕對(duì)Q值對(duì)其標(biāo)定，得到表層等效的平面Q 場(chǎng)。將該表層平面Q 場(chǎng)鑲嵌在用Q?V 擬合關(guān)系轉(zhuǎn)成的近地表Q 場(chǎng)的淺層，得到該工區(qū)最終的近地表Q 場(chǎng)。最后對(duì)工區(qū)內(nèi)的地震數(shù)據(jù)進(jìn)行表層Q 補(bǔ)償，消除了近地表對(duì)地震資料的影響，提高了地震資料的分辨率。

1 近地表Q 場(chǎng)計(jì)算方法及實(shí)現(xiàn)步驟

1. 1 改進(jìn)的譜比法求Q 值

首先處理單井微測(cè)井于井下的12 級(jí)檢波器采集的數(shù)據(jù)。對(duì)該數(shù)據(jù)進(jìn)行初至拾取，通過波場(chǎng)分離得到下行波，選取子波時(shí)窗范圍內(nèi)的數(shù)據(jù)，對(duì)每道數(shù)據(jù)進(jìn)行傅里葉變換，求取每道頻譜對(duì)數(shù)譜的截距和斜率，根據(jù)截距和斜率將每道的對(duì)數(shù)譜拉平，將下行波對(duì)數(shù)譜疊加求取未衰減的擬下行子波譜，利用各道對(duì)數(shù)譜減去擬下行子波譜，得到初步消除子波影響的各道對(duì)數(shù)譜。在此基礎(chǔ)上求取每道對(duì)數(shù)譜的斜率，對(duì)其進(jìn)行平滑擬合，保證斜率值隨深度的增加而增大，再利用平滑后的斜率和初至求取等效Q，再根據(jù)等效Q 與層Q 的關(guān)系，求取層Q 值。通過以上步驟求取的Q 可用于后續(xù)的反Q 濾波，以獲得高分辨率的剖面成果。

依據(jù)該方法，首先進(jìn)行模型試算。設(shè)計(jì)一個(gè)5 層的水平層狀模型（表1），激發(fā)子波為Ricker 子波，主頻為50 Hz，檢波器的深度范圍為120～1660 m，檢波點(diǎn)間隔為20 m，井源距為0，垂直入射，正演含Q 值衰減的下行波記錄。從正演的下行波記錄（圖1）可以看出，隨著深度的增加，振幅能量發(fā)生衰減，相位改變，子波波形變胖，從記錄的頻譜上也可以看出，振幅能量減弱，主頻降低，頻帶變窄。

圖2 展示了改進(jìn)譜比法Q 值求取的過程。圖2a是圖1 中下行子波記錄的振幅譜對(duì)數(shù)譜，選擇兩個(gè)頻率點(diǎn)，一般選擇主頻和高頻點(diǎn)的位置，將每道振幅譜對(duì)數(shù)按照該斜率拉平，然后多道數(shù)據(jù)進(jìn)行疊加平均，得到擬下行子波的頻譜（圖2b）。然后將每道的振幅譜與擬下行子波頻譜進(jìn)行譜比，得到譜比曲線（圖2c）。在單調(diào)遞減的頻率段上，計(jì)算每道數(shù)據(jù)的斜率，繪出斜率隨深度變化的關(guān)系（圖2d），斜率值應(yīng)該隨深度的增加而增加。其中斜率G = πt/Qeff，t 為初至?xí)r間，Qeff 為每道的等效Q 值，所以Qeff = πt/G。再根據(jù)Qeff 轉(zhuǎn)層Q 的關(guān)系式（1），計(jì)算第n 層的Q 值（Qint）即層Q 值等于相鄰道間的時(shí)差除以相鄰道間的斜率差，初至?xí)r間和斜率必須是隨深度遞增的，才能保證計(jì)算的Q 值沒有負(fù)值。圖2e 是模型數(shù)據(jù)相鄰道間的斜率差，可以看出，同一地層的斜率差相等，圖2f 是根據(jù)式（2）計(jì)算的層Q 值。

模型試算的數(shù)據(jù)考慮的是垂直入射的情況，當(dāng)炮檢距為非零時(shí)，需要考慮對(duì)射線路徑進(jìn)行校正，將初至?xí)r間和斜率校正到垂直入射的情況，才能得到正確的解。如圖3 所示，可以看出不考慮射線路徑校正時(shí)，特別是在分界面位置存在較明顯的抖動(dòng)，淺層計(jì)算結(jié)果與真實(shí)值差距較大，但對(duì)深層結(jié)果影響較小。當(dāng)考慮射線路徑進(jìn)行校正后，可以得到正確Q 值解。

對(duì)實(shí)際單井微測(cè)井?dāng)?shù)據(jù)進(jìn)行試算。此次單井微測(cè)井?dāng)?shù)據(jù)，選用12 級(jí)檢波器同時(shí)接收的記錄，保證了激發(fā)子波的一致性，同時(shí)12 級(jí)檢波器已覆蓋表層，計(jì)算結(jié)果能夠反應(yīng)表層的Q 值。具體做法如下：

（1）井中檢波器接收，地面利用爆炸震源或重錘在近井口激發(fā)，多級(jí)檢波器同時(shí)采集數(shù)據(jù)，保證了激發(fā)子波的一致性。

（2）對(duì)接收數(shù)據(jù)進(jìn)行預(yù)處理，剔除壞道，經(jīng)過拾取縱波初至和上下行波分離處理，得到預(yù)處理后下行縱波數(shù)據(jù)和縱波初至?xí)r間。

（3）利用步驟（2）中所得到的數(shù)據(jù)，截取出下行子波數(shù)據(jù)，對(duì)每道數(shù)據(jù)進(jìn)行FFT 變換，得到每道數(shù)據(jù)的對(duì)數(shù)頻譜，即對(duì)FFT 變換后的振幅譜取對(duì)數(shù)。

（4）取頻譜中的峰值頻率和高截頻，計(jì)算每道對(duì)數(shù)頻譜的斜率和截距。

（5）對(duì)每道的對(duì)數(shù)頻譜按照步驟（4）中的斜率和截距進(jìn)行拉平處理，并進(jìn)行疊加，得到初始下行子波的對(duì)數(shù)頻譜。

（6）每道下行波對(duì)數(shù)頻譜減去步驟（5）的下行子波對(duì)數(shù)譜，得到去除下行子波影響的下行子波對(duì)數(shù)頻譜。

（7）在峰值頻率和高截頻范圍內(nèi)，做線性擬合，擬合的斜率可以計(jì)算該道數(shù)據(jù)對(duì)應(yīng)的等效Q 值，Qeff =πt/G。這里G 應(yīng)為增函數(shù)，但由于數(shù)值誤差，估算的G 值可能產(chǎn)生反向的波動(dòng)，為保證求解的穩(wěn)定性，可以按不同地層，對(duì)G 進(jìn)行多點(diǎn)平滑或擬合處理，保證G 隨深度變化為增函數(shù)，這樣求取的層Q 值不會(huì)出現(xiàn)負(fù)值。

（8）利用式（2）求第n 層的層Q 參數(shù)，即可得到每個(gè)深度點(diǎn)的層Q 值。

下面進(jìn)行實(shí)際單井微測(cè)井?dāng)?shù)據(jù)的Q 值估算。圖4a 是實(shí)際單井微測(cè)井?dāng)?shù)據(jù)及其振幅譜，圖4b 是截取的微測(cè)井的初至子波及其頻譜，由圖可見，隨著檢波點(diǎn)深度的增加，頻帶逐漸變窄。利用截取的子波數(shù)據(jù)來進(jìn)行Q 值的估算。

圖５是實(shí)際的單井微測(cè)井?dāng)?shù)據(jù)估算Q 值的過程圖。圖5a 是每道下行子波的經(jīng)傅里葉變換的振幅譜取對(duì)數(shù)的曲線；圖5b 是統(tǒng)計(jì)平均得到的擬下行子波的振幅對(duì)數(shù)譜；圖5c 是每道數(shù)據(jù)的振幅譜與統(tǒng)計(jì)子波譜的譜比曲線，選擇頻段范圍，每一道的譜比曲線可以求得一個(gè)斜率；圖5d 是每道數(shù)據(jù)的初至和斜率交會(huì)圖。由于計(jì)算的斜率與初至的關(guān)系不一定是單調(diào)遞增的，可以按照微測(cè)井解釋成果的層位，進(jìn)行擬合平滑，保證層位內(nèi)的斜率是單調(diào)遞增，從而保證結(jié)果的穩(wěn)定性。圖5e 是按照微測(cè)井分層的第一層，得到的斜率與時(shí)間的擬合關(guān)系；圖5f 是按照微測(cè)井分層第二層，得到的的斜率與時(shí)間的擬合關(guān)系；圖5g 是擬合后的斜率與計(jì)算的斜率的對(duì)比圖，擬合后的斜率與圖5j 每道數(shù)據(jù)的初至?xí)r間存在對(duì)應(yīng)關(guān)系；圖5h 相鄰道間斜率差和圖5k 相鄰道間時(shí)差也存在對(duì)應(yīng)關(guān)系。由斜率差可得出圖5i 層Q 值；由初至?xí)r差可得出圖5l 層速度。因此可見，層Q 值與層速度間存在著一定的對(duì)應(yīng)關(guān)系。

1. 2 多炮擬合法求Qeff

采用肖永新等[10]提出的多炮擬合算法，利用雙井微測(cè)井資料計(jì)算表層Q 值。雙井微測(cè)井調(diào)查方法以兩口淺井開展，一口為激發(fā)井，一口為接收井。具體做法是：將地層看成表層和高速層兩層模型，選取高速層中某一炮的地面道和井底道數(shù)據(jù)，利用譜比法求取表層的等效Q 值（Q0）。該方法對(duì)于表層較薄的地層適用。

圖6 為兩層速度模型的雙井微測(cè)井射線傳播示意圖。假設(shè)表層Q 值為Q0，高速層Q 值為Q1，將高速層中某一炮點(diǎn)的地面道R1 與井底道R4 進(jìn)行譜比

式中：B1、 B4 分別為R1、R4 道的振幅譜；f 為頻率；C為兩道譜比取對(duì)數(shù)的常數(shù)項(xiàng)。對(duì)同一激發(fā)點(diǎn)，地面道與井底道振幅譜比取對(duì)數(shù)，選擇頻段范圍，線性擬合，得到斜率πt2 /Q0，又已知速度模型，利用射線追蹤計(jì)算t2，即可得到表層Q0。

實(shí)際野外施工過程中，很難提前知道滿足t1－t3=0 的激發(fā)點(diǎn)位置，所以不可能在該深度設(shè)置激發(fā)點(diǎn)，因此選取高速頂以下多個(gè)炮點(diǎn)進(jìn)行表層Q 值計(jì)算，然后再對(duì)多個(gè)不同深度炮點(diǎn)計(jì)算的表層Q 值進(jìn)行擬合，找到t1－t3=0 位置時(shí)的Q 值，即為該點(diǎn)的表層等效Q值，即為Q0。

由圖7 可知，隨著激發(fā)點(diǎn)深度的增加，Q 值逐漸增大，這里的Q 值指的是該激發(fā)點(diǎn)深度位置以上地層的等效Q 值（Q0）。令DT=t1－t3，隨著DT 的增加，表層Q 值（Q0）也逐漸增大。通過擬合DT-Q 曲線，在DT=0 處的表層Q 值（Q0）即為所求位置處的表層等效Q 值。

圖8a 為四川盆地蓬萊氣區(qū)內(nèi)的一口雙井微測(cè)井觀測(cè)系統(tǒng)示意圖及近地表結(jié)構(gòu)圖，圖中紅色點(diǎn)為炮點(diǎn)，綠色點(diǎn)為接收點(diǎn)。該雙井微測(cè)井點(diǎn)井深17 m 兩口井之間的距離5 m；圖8 右為該口雙井微測(cè)井?dāng)?shù)據(jù)的地面接收道記錄和井底接收道記錄及其頻譜。

首先拾取地震數(shù)據(jù)的初至?xí)r間，通過常規(guī)的微測(cè)井解釋方法獲得近地表速度模型。然后根據(jù)近地表速度模型信息，選擇高速層中的炮點(diǎn)，本例中，選擇激發(fā)深度5～16 m 的炮點(diǎn)（激發(fā)深度17 m 處的地面接收點(diǎn)記錄存在異常，已剔除）。選擇高速頂之下的每一炮的一個(gè)地面接收道和井底接收道，每一炮的地面道和井底道采用常規(guī)譜比法計(jì)算該炮點(diǎn)深度位置以上的等效Q 值。根據(jù)圖8 所示近地表速度模型進(jìn)行射線追蹤，得到每炮地面道和井底道的射線路徑，及每一段射線的傳播時(shí)間。計(jì)算每一炮的地面道和井底道在高速層中的傳播時(shí)間差，即t1-t3。Q 值與t1-t3的關(guān)系如圖9右所示，多項(xiàng)式擬合t1-t3 和Q 值，當(dāng)t1-t3=0 時(shí)的Q值即為所求Q0。本例中，該微測(cè)井點(diǎn)表層Q0 為3. 84。

1. 3 平面Q 建場(chǎng)

采用地面地震的地表一致性振幅補(bǔ)償相對(duì)系數(shù)和表層旅行時(shí)，通過譜比法求取表層相對(duì)Q 場(chǎng)。此時(shí)求得的相對(duì)Q 場(chǎng)也是將表層看作一個(gè)等效地層。再用雙井微測(cè)井的計(jì)算出的Q 值進(jìn)行標(biāo)定，得到最終的表層平面Q 場(chǎng)。將雙井微測(cè)井求得的Q 值稱為絕對(duì)Q 值，將地面地震振幅信息求得的Q 場(chǎng)稱為相對(duì)Q場(chǎng)，然后用絕對(duì)Q 值對(duì)相對(duì)Q 場(chǎng)進(jìn)行約束，采用反距離加權(quán)內(nèi)插法[18]對(duì)相對(duì)Q 場(chǎng)進(jìn)行標(biāo)定。基于譜比法求相對(duì)Q 值

式中：R 為地表一致性振幅相對(duì)系數(shù)；scale 為人為給定比例因子，目的是使R × scale 小于1，保證計(jì)算出Q 的值大于0。

計(jì)算步驟：①由地表一致性振幅補(bǔ)償模塊做振幅相對(duì)系數(shù)的計(jì)算，沿初至開小窗；②導(dǎo)出文本，利用相對(duì)Q 值的表達(dá)式，求解相對(duì)Q 場(chǎng)；③再用雙井微測(cè)井計(jì)算的Q 值進(jìn)行標(biāo)定，可以求出平面Q 場(chǎng)。平面Q場(chǎng)的建立流程如圖10。

圖11 展示了工區(qū)實(shí)際表層平面Q 場(chǎng)的建立過程。由圖11a 的表層旅行時(shí)和圖11b 的相對(duì)振幅系數(shù)，利用相對(duì)Q 值的式（6），求解相對(duì)Q 場(chǎng)（圖11c）；再用單、雙井微測(cè)井計(jì)算的表層Q 值插值成平面的Q場(chǎng)（圖11d），對(duì)計(jì)算的相對(duì)Q 場(chǎng)采用反距離加權(quán)內(nèi)插方法進(jìn)行約束，可以求出等效的平面Q 場(chǎng)（圖11e）。

1. 4 融合立體Q 建場(chǎng)

根據(jù)實(shí)際的單井微測(cè)井資料得到速度和Q 值的擬合關(guān)系，利用該擬合關(guān)系，將層析反演得到的近地表速度場(chǎng)轉(zhuǎn)換成Q 場(chǎng)（圖12）。由于層析反演的近地表速度場(chǎng)的深度采樣間隔為10 m，而該工區(qū)的微測(cè)井解釋成果中，表層厚度較小，一般為2～8 m，表層反演的速度較微測(cè)井實(shí)測(cè)速度要大，所以速度場(chǎng)轉(zhuǎn)換的Q 場(chǎng)不能反應(yīng)表層真正的Q 值，所以將雙井微測(cè)井結(jié)合地面地震信息求取的表層平面Q 場(chǎng)鑲嵌在由Q?V 關(guān)系轉(zhuǎn)換后的近地表Q 場(chǎng)的表層，得到最終的近地表Q 場(chǎng)（圖13）。

2 應(yīng)用效果

為了能恢復(fù)工區(qū)復(fù)雜地表情況所導(dǎo)致的振幅能量衰減和相位畸變等情況，將工區(qū)近地表綜合調(diào)查微測(cè)井資料和地面地震資料綜合運(yùn)用，進(jìn)行近地表Q 補(bǔ)償，具體步驟如下：

（1）利用工區(qū)近地表微測(cè)井資料求取表層的絕對(duì)Q 值，雙井微測(cè)井資料利用多炮擬合法可求取表層薄層的Q 值，單井微測(cè)井資料利用改進(jìn)的譜比法求Q值，進(jìn)而可得到表層速度和Q 值的關(guān)系；

（2）根據(jù)地震資料計(jì)算各個(gè)檢波點(diǎn)響應(yīng)的相對(duì)振幅系數(shù)和近地表層析反演的表層旅行時(shí)，通過譜比法計(jì)算出相對(duì)Q 值；

（3）用求取的絕對(duì)Q 值來標(biāo)定相對(duì)Q 值，進(jìn)而獲得優(yōu)化后的表層平面Q 場(chǎng)；

（4）對(duì)于表層以下到高速頂處的降速帶部分，采用步驟（1）得到的Q?V 擬合關(guān)系，將近地表速度場(chǎng)轉(zhuǎn)換成Q 場(chǎng)；

（5）將優(yōu)化后的表層平面Q 場(chǎng)鑲嵌在Q?V 關(guān)系轉(zhuǎn)換后的Q 體的表層；

（6）利用求取的優(yōu)化Q 場(chǎng)，使用穩(wěn)定因子的Q 補(bǔ)償算法對(duì)工區(qū)實(shí)際地震數(shù)據(jù)體進(jìn)行Q 補(bǔ)償處理，具體補(bǔ)償效果詳見后面內(nèi)容。

通過使用上文所提及的方法對(duì)工區(qū)進(jìn)行Q 補(bǔ)償處理，對(duì)比分析Q 補(bǔ)償前、后的地震數(shù)據(jù)的一致性、分辨率和頻譜特征。圖14 為近地表Q 補(bǔ)償前、后的單炮記錄對(duì)比圖，可見Q 補(bǔ)償后的單炮記錄同相軸更細(xì)，橫向連續(xù)性得以改善，子波一致性更好。

由于本文涉及的工區(qū)低降速帶厚度較小，表層起伏，使得不同位置處的子波波形存在差異。如圖15所示，經(jīng)過近地表Q 補(bǔ)償后，地震剖面的分辨率得到了提高，部分小細(xì)軸得到了恢復(fù)，從自相關(guān)記錄上可以看出，子波的橫向一致性得到了改善，從頻譜圖上可以看出，應(yīng)用近地表Q 補(bǔ)償后，頻帶拓寬6 Hz 左右（圖16）。

對(duì)近地表Q 補(bǔ)償前、后的地震剖面的第一零交叉時(shí)平面屬性進(jìn)行質(zhì)控，如圖17 所示，可以看出，應(yīng)用近地表Q 補(bǔ)償后，第一零交叉時(shí)的值變小，說明子波得到了壓縮，而且從整個(gè)工區(qū)的平面圖上可以看出，應(yīng)用近地表Q 補(bǔ)償后，第一零交叉時(shí)的變化趨于一致，說明子波一致性變好。

將近地表Q 補(bǔ)償前、后的地震記錄與合成記錄進(jìn)行井震標(biāo)定，可以看出，應(yīng)用近地表Q 補(bǔ)償后的地震記錄與合成記錄的相關(guān)系數(shù)是0. 72，高于沒有應(yīng)用近地表Q補(bǔ)償?shù)牡卣鹩涗浥c合成記錄的相關(guān)系數(shù)，（圖18）。

3 結(jié)論

利用單井微測(cè)井資料，采用改進(jìn)的譜比法，計(jì)算Q 值，然后擬合Q?V 關(guān)系，可以很好地將層析反演出的降速帶的速度場(chǎng)轉(zhuǎn)換成Q 場(chǎng)。利用雙井微測(cè)井資料，采用多炮擬合的算法，計(jì)算井點(diǎn)位置處的絕對(duì)Q值，再結(jié)合地面地震的振幅和走時(shí)信息計(jì)算相對(duì)Q場(chǎng)，再利用絕對(duì)Q 值對(duì)相對(duì)Q 場(chǎng)進(jìn)行標(biāo)定，可以得到表層低速帶的等效平面Q 場(chǎng)。由于該工區(qū)表層厚度較薄，一般為2～8 m，而層析反演的速度場(chǎng)的采樣間隔是10 m，很難將小于10 m 的低速帶反應(yīng)出來。所以，將標(biāo)定后的表層平面Q 場(chǎng)鑲嵌在由Q?V 關(guān)系轉(zhuǎn)換的近地表Q 場(chǎng)的淺表層，得到最終的近地表Q 場(chǎng)。然后應(yīng)用此Q 場(chǎng)，對(duì)地面地震資料進(jìn)行近地表Q 補(bǔ)償，可以有效消除近地表對(duì)地震資料的吸收衰減影響，改善子波的一致性，有效拓寬地震資料的頻帶，提高資料的分辨率。