Page 71 - 《应用声学》2020年第1期

第 39 卷 第 1 期 印兴耀等： 五维地震油气识别方法 67

( )
(
)
2
2
2
2
d = sin θ cos φ 1 + tan θ sin φ , ρ (V )
+ c(θ) ln + d(θ, φ)δ
2
4
2
2
2
2
e = cos φ sin θ tan θ, f = −8k sin θ cos φ, ρ 0
+ e(θ, φ)ε (V ) + f(θ, φ)γ (V ) , (15)
M 0 、µ 0 和ρ 0 分别为纵横波模量和密度的均值，标准
化的方位弹性阻抗方程克服了不同角度弹性阻抗 1
其中，EI 0 = (M 0 ρ 0 ) 。从方程 (15) 可以看出利用
2
量纲不统一的不足，提高了算法的稳定可靠性。
线性化后的方位弹性阻抗方程至少需要 6 个不同方
为了能够从弹性阻抗数据体中反演出稳定的
位角、不同入射角的弹性阻抗数据体以提取储层的
弹性参数和各向异性参数数据体，需要对方程 (14)
弹性参数(M 0 、µ 0 和ρ 0 )和各向异性参数 (δ (V ) 、ε (V )
进行线性化处理，即：
( ) ( ) ( ) 和 γ (V ) )。将 6 个不同方位、不同入射角的弹性阻抗
EI M µ
ln = a(θ) ln + b(θ) ln 方程写成矩阵形式为
EI 0 M 0 µ 0
 ( EI (θ 1 , φ 1 ) )
ln
 ( )
M  EI 0 
 ln  ( ) 
  
   EI (θ 2 , φ 1 ) 
a (θ 1 ) b (θ 1 ) c (θ 1 ) d (θ 1 , φ 1 ) e (θ 1 , φ 1 ) f (θ 1 , φ 1 )  M 0  ln
( ) 
 
   µ   EI 0 

 a (θ 2 ) b (θ 2 ) c (θ 2 ) d (θ 2 , φ 1 ) e (θ 2 , φ 1 ) f (θ 2 , φ 1 )  ln   ( )
   EI (θ 3 , φ 1 ) 
 µ 0 
   ln 
 ( ) 
   
a (θ 3 ) b (θ 3 ) c (θ 3 ) d (θ 3 , φ 1 ) e (θ 3 , φ 1 ) f (θ 3 , φ 1 )  ρ  EI 0
   
 ln  ) . (16)
  =  (
 
  ρ 0  EI (θ 1 , φ 2 ) 
a (θ 1 ) b (θ 1 ) c (θ 1 ) d (θ 1 , φ 2 ) e (θ 1 , φ 2 ) f (θ 1 , φ 2 )   ln
  
 
  (V )  
 δ  EI 0
   
a (θ 2 ) b (θ 2 ) c (θ 2 ) d (θ 2 , φ 2 ) e (θ 2 , φ 2 ) f (θ 2 , φ 2 )   ( )
   
 (V )  EI (θ 2 , φ 2 ) 
ε ln

 
a (θ 3 ) b (θ 3 ) c (θ 3 ) d (θ 3 , φ 2 ) e (θ 3 , φ 2 ) f (θ 3 , φ 2 )   
 EI 0 
γ (V )  ( )
 EI (θ 3 , φ 2 ) 
ln
EI 0
利用已知的测井数据以及井旁道方位弹性阻 射振幅产生影响，因此在相同入射角的情况下，不同
抗值计算公式 (16) 中各参数对应的常系数。在储 方位的地震数据之间存在明显的差异。
层弹性性质横向变化较小的情况下，可将常系数看 根据图 2 五维地震数据，结合该区域井资料以
作是恒定不变的值。将叠前方位弹性阻抗反演得 及岩石物理信息作为约束，反演得到不同方位的大
到的 6 个不同方位、不同入射角的弹性阻抗代入公 中小角度弹性阻抗如图 3 所示。通过基于方位弹性
式 (16) 中，结合各参数对应的常系数，通过求解可 阻抗的各向异性参数反演方法从方位弹性阻抗数
得到各采样点处的弹性参 ρ 和各向异性参数。由于 据体中提取地层的纵波模量、横波模量、密度和HTI
该方法在进行反演时使用了测井资料进行约束，并 介质的各向异性参数，进而求取剖面含流体区域的
结合地震资料的先验信息，因此得到的弹性参数和 各向异性因子，如图4 所示，A 井和 B 井的有色部分
各向异性参数较为准确可靠。 代表流体区，白色部分是干岩区。图 4 中的橘色箭
头所指部分为测井解释的流体部分，发现该工区储
3 实际应用 层中，流体分布区域的各向异性因子在反演数值呈
现低值特征，而且各向异性因子数值在岩石发育位
五维地震资料来自中国东部某宽方位地震勘 置展现出高值特征，该特征与岩石物理分析规律相
探工区，该区域各向异性特征明显。对五维地震数 吻合。同时，含流体裂缝性储层各向异性流体的估
据首先进行预处理得到 6 个方位且每个方位包含 3 算剖面与 A 井和井 B 的流体解释结果具有良好的
个入射角的共18个分方位部分角度叠加道集，然后 匹配性。因此，通过各向异性反演各向异性因子可
选取其中6 个叠前方位地震数据进行叠前方位弹性 展示储层流体分布特征，为五维地震定量解释创建
阻抗反演，叠前方位地震数据如图2所示。从图2中 新的理论方法，该方法为复杂裂缝储层地震流体识
可以看出由于储层中存在的裂缝对地震资料的反 别提供了一种新的途径。