[發(fā)明專利]基于散射波地震響應(yīng)特征的裂縫間距定量分析方法和設(shè)備

有效

• 申請(qǐng)?zhí)枺?/td>: 201910818522.5
• 申請(qǐng)日: 2019-08-30
• 公開（公告）號(hào): CN112444871B
• 公開（公告）日: 2023-07-04
• 發(fā)明（設(shè)計(jì)）人: 肖云飛
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油化工股份有限公司;中國(guó)石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號(hào): G01V1/30
• 分類號(hào): G01V1/30
• 代理公司: 北京聿宏知識(shí)產(chǎn)權(quán)代理有限公司 11372
• 代理人: 吳大建;張杰

1.一種基于散射波地震響應(yīng)特征的裂縫間距定量分析方法，其特征在于，包括以下步驟：

S100，對(duì)目標(biāo)儲(chǔ)層的地震響應(yīng)進(jìn)行分析，以獲取目標(biāo)儲(chǔ)層的包含裂縫地震響應(yīng)的散射波場(chǎng)記錄；

S200，從所述散射波場(chǎng)記錄中提取含有裂縫信息的裂縫散射波場(chǎng)；

S300，從所述裂縫散射波場(chǎng)中分解出后向散射波場(chǎng)；

S400，將所述后向散射波場(chǎng)變換到頻率－波數(shù)域中，并確定所述后向散射波場(chǎng)的頻率－波數(shù)域譜中最大后向散射能量值所對(duì)應(yīng)的波數(shù)值；

S500，根據(jù)所述波數(shù)值估算所述目標(biāo)儲(chǔ)層中裂縫的間距；

利用局部波場(chǎng)分解法，從所述裂縫散射波場(chǎng)中分解出后向散射波場(chǎng)；

所述利用局部波場(chǎng)分解法，從所述裂縫散射波場(chǎng)中分解出后向散射波場(chǎng)，具體為：

基于局部波場(chǎng)分解算子從所述裂縫散射波場(chǎng)中分解出后向散射波場(chǎng)，其中所述局部波場(chǎng)分解算子為：

式中，a為半偏移孔徑算子，x表示空間位置坐標(biāo)，h(x)為位于x位置處的空間錐形濾波器，用于把空間錐形區(qū)域里面的數(shù)據(jù)分解出來(lái)，為地震子波的傅里葉變化，φ(p,x)為積分路徑，ω為頻率，p為射線參數(shù)，i表示虛數(shù)單位；

所述基于局部波場(chǎng)分解算子利用局部波場(chǎng)分解法從所述裂縫散射波場(chǎng)中分解出后向散射波場(chǎng)，具體包括以下步驟：

基于局部波場(chǎng)分解算子，從所述裂縫散射波場(chǎng)中分解出不同射線參數(shù)p值所對(duì)應(yīng)的散射波場(chǎng)記錄；

對(duì)不同射線參數(shù)p值所對(duì)應(yīng)的散射波場(chǎng)記錄進(jìn)行能量歸一化分析，基于能量歸一化分析結(jié)果從所有散射波場(chǎng)記錄中找出能量大于預(yù)設(shè)百分比閾值且射線參數(shù)p值為負(fù)值的散射波場(chǎng)記錄，作為后向散射波場(chǎng)。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟S200包括以下步驟：

比較所述散射波場(chǎng)記錄中垂直裂縫方向和平行裂縫方向的波場(chǎng)散射特征；若垂直裂縫方向/平行裂縫方向的波場(chǎng)散射特征較強(qiáng)，則提取所述散射波場(chǎng)記錄中垂直裂縫方向/平行裂縫方向的波場(chǎng)分量，作為含有裂縫信息的裂縫散射波場(chǎng)。

3.根據(jù)權(quán)利要求1所述的方法，其特征在于，對(duì)不同射線參數(shù)值所對(duì)應(yīng)的散射波場(chǎng)記錄進(jìn)行能量歸一化分析，基于能量歸一化分析結(jié)果從所有散射波場(chǎng)記錄中找出能量大于預(yù)設(shè)百分比閾值且射線參數(shù)值為負(fù)值的散射波場(chǎng)記錄，作為后向散射波場(chǎng)，具體包括以下步驟：

統(tǒng)計(jì)每個(gè)p值對(duì)應(yīng)的散射波場(chǎng)記錄的能量，全部繪制在橫坐標(biāo)為p值、縱坐標(biāo)為能量的曲線圖中，然后進(jìn)行能量歸一化處理，取所有散射波場(chǎng)記錄中能量大于預(yù)設(shè)百分比閾值且射線參數(shù)值為負(fù)值的散射波場(chǎng)記錄，作為后向散射波場(chǎng)。

4.根據(jù)權(quán)利要求1或3所述的方法，其特征在于，所述百分比閾值為80％。

5.根據(jù)權(quán)利要求1所述的方法，其特征在于，在所述步驟S500中，按照下式根據(jù)所述波數(shù)值估算所述目標(biāo)儲(chǔ)層中裂縫的間距：

式中，L是裂縫的間距，k_a是最大后向散射能量值所對(duì)應(yīng)的波數(shù)值。

6.一種計(jì)算機(jī)存儲(chǔ)介質(zhì)，其特征在于，其中存儲(chǔ)有用于實(shí)現(xiàn)上述權(quán)利要求1至5中任意一項(xiàng)所述方法的計(jì)算機(jī)程序。

7.一種計(jì)算機(jī)設(shè)備，其特征在于，包括存儲(chǔ)器和處理器，所述處理器用于執(zhí)行所述存儲(chǔ)器中存儲(chǔ)的計(jì)算機(jī)程序，所述計(jì)算機(jī)程序用于實(shí)現(xiàn)上述權(quán)利要求1至5中任意一項(xiàng)所述方法。