[發(fā)明專利]一種地表主動(dòng)源反射波干涉成像方法

有效

• 申請(qǐng)?zhí)枺?/td>: 202110631426.7
• 申請(qǐng)日: 2021-06-07
• 公開(kāi)（公告）號(hào): CN113406700B
• 公開(kāi)（公告）日: 2022-11-22
• 發(fā)明（設(shè)計(jì)）人: 喬寶平;曹成寅;潘自強(qiáng);趙丹;黃偉傳;黃昱丞
• 申請(qǐng)（專利權(quán)）人: 核工業(yè)北京地質(zhì)研究院
• 主分類號(hào): G01V1/36
• 分類號(hào): G01V1/36
• 代理公司: 核工業(yè)專利中心 11007
• 代理人: 陳麗麗

1.一種地表主動(dòng)源反射波干涉成像方法，其特征在于，所述方法包括以下步驟：

步驟(1)、獲取共成像點(diǎn)道集地震數(shù)據(jù)；

步驟(2)、獲取基于一致性譜相關(guān)干涉重構(gòu)的共成像點(diǎn)虛反射道集地震數(shù)據(jù)；

步驟(3)、獲取基于互相關(guān)干涉重構(gòu)的共成像點(diǎn)道集地震數(shù)據(jù)；

步驟(4)、獲取基于水平疊加的目標(biāo)層反射波成像結(jié)果；

所述步驟(1)包括：

步驟(1.1)、檢查地震野外班報(bào)、SPS文件和地震原始數(shù)據(jù)；

步驟(1.2)、將地震野外班報(bào)和SPS文件中的觀測(cè)系統(tǒng)信息加載到地震原始數(shù)據(jù)中，計(jì)算生成共中心點(diǎn)號(hào)、偏移距和方位角的觀測(cè)系統(tǒng)信息，獲取具有觀測(cè)系統(tǒng)信息的地震數(shù)據(jù)；

步驟(1.3)、利用加載完觀測(cè)系統(tǒng)信息的地震數(shù)據(jù)，采用自動(dòng)加人工修正的方法拾取初至波走時(shí)信息，并利用初至波走時(shí)信息開(kāi)展折射靜校正處理，計(jì)算折射靜校正后的地震數(shù)據(jù)；

步驟(1.4)、利用折射靜校正后的地震數(shù)據(jù)，采用自適應(yīng)噪聲壓制方法開(kāi)展隨機(jī)噪聲、異常極值、側(cè)面反射波、繞射波和地震面波的噪聲壓制處理，生成噪聲壓制處理后的地震數(shù)據(jù)；

步驟(1.5)、利用噪聲壓制處理后的地震數(shù)據(jù)，開(kāi)展地表一致性振幅補(bǔ)償處理，彌補(bǔ)震源點(diǎn)和接收點(diǎn)產(chǎn)生的振幅非均勻性，獲取振幅補(bǔ)償后的地震數(shù)據(jù)；

步驟(1.6)、利用振幅補(bǔ)償后的地震數(shù)據(jù)，開(kāi)展地表一致性反褶積和脈沖反褶積處理，生成反褶積處理后的地震數(shù)據(jù)；

步驟(1.7)、利用反褶積處理后的地震數(shù)據(jù)，開(kāi)展分頻去噪處理，壓制步驟(1.6)中產(chǎn)生的高頻噪聲干擾，獲取高頻噪聲壓制后的地震數(shù)據(jù)。

2.根據(jù)權(quán)利要求1所述的一種地表主動(dòng)源反射波干涉成像方法，其特征在于，所述步驟(2)包括：

步驟(2.1)、利用高頻噪聲壓制后的地震數(shù)據(jù)，針對(duì)目標(biāo)層開(kāi)展高精度的速度分析，獲取目標(biāo)層速度場(chǎng)數(shù)據(jù)；

步驟(2.2)、利用步驟(1.7)生成的高頻噪聲壓制后的地震數(shù)據(jù)和步驟(2.1)生成的目標(biāo)層速度場(chǎng)數(shù)據(jù)，開(kāi)展局部拉伸清零的動(dòng)校正處理，獲取局部拉伸清零動(dòng)校正處理后的地震數(shù)據(jù)；

步驟(2.3)、利用局部拉伸清零動(dòng)校正處理后的地震數(shù)據(jù)，通過(guò)水平疊加，獲取目標(biāo)層的反射波場(chǎng)數(shù)據(jù)；

步驟(2.4)、利用步驟(1.7)生成的地震數(shù)據(jù)與步驟(2.3)獲取的反射波場(chǎng)數(shù)據(jù)，開(kāi)展一致性譜相關(guān)地震干涉處理，以求取代表動(dòng)校正量的虛反射道集。

3.根據(jù)權(quán)利要求2所述的一種地表主動(dòng)源反射波干涉成像方法，其特征在于，所述步驟(2.4)中求取代表動(dòng)校正量的虛反射道集的公式為：

其中，*代表了復(fù)共軛，ω代表了角頻率，α代表了穩(wěn)定因子，R_ZERO(ω)代表了步驟(2.3)生成的目標(biāo)層位反射波場(chǎng)數(shù)據(jù)，G(ω)代表了步驟(1.7)生成的高頻噪聲壓制后的地震數(shù)據(jù)，V(ω)代表了步驟(2.4)重構(gòu)的虛反射道集。

4.根據(jù)權(quán)利要求3所述的一種地表主動(dòng)源反射波干涉成像方法，其特征在于，所述步驟(3)具體為：利用步驟(2.4)生成的虛反射道集與步驟(1.7)獲取的高頻噪聲壓制后的地震數(shù)據(jù)，開(kāi)展基于互相關(guān)的地震干涉處理，以重構(gòu)相位校正后的共成像點(diǎn)道集地震數(shù)據(jù)。

5.根據(jù)權(quán)利要求4所述的一種地表主動(dòng)源反射波干涉成像方法，其特征在于，所述步驟(3)中求取相位校正后的共成像點(diǎn)道集地震數(shù)據(jù)的公式為：

C(ω)＝G(ω)V(ω)^*

其中，C(ω)代表了相位校正后的共成像點(diǎn)道集地震數(shù)據(jù)。

6.根據(jù)權(quán)利要求4或5所述的一種地表主動(dòng)源反射波干涉成像方法，其特征在于，所述步驟(4)具體為：針對(duì)步驟(3)生成的相位校正后的共成像點(diǎn)道集地震數(shù)據(jù)，開(kāi)展水平疊加處理，計(jì)算生成針對(duì)目標(biāo)層的反射波成像結(jié)果。

7.根據(jù)權(quán)利要求6所述的一種地表主動(dòng)源反射波干涉成像方法，其特征在于，所述步驟(4)中求取針對(duì)目標(biāo)層的反射波成像結(jié)果的公式為：

I(ω)＝∑C(ω)

其中，∑代表了疊加求和，I(ω)代表了通過(guò)水平疊加生成的成像結(jié)果。