Page 97 - 《应用声学》2022年第4期
P. 97

第 41 卷 第 4 期              石志奇等: Hudson 模型适用性及储层微观建模应用                                       595


             λ, µ 为各向同性背景基质的弹性模量。U kl 为散射                      场快照如图 1 所示。对于干燥裂缝,随着裂缝密度
             算子,与裂缝本身的连通状态和充填状态有关。对                            增大,准纵波和准 SV波的各向异性逐渐增强,当裂
             于孤立的干燥裂缝,                                         缝密度大于 0.1 时,一阶模型和二阶模型的波场快
                                   16(λ + 2µ)                 照差异变得明显,主要体现在一阶模型波场的各向
                       ¯     ¯               ,
                      U 11 = U 22 =
                                    3(3λ + 4µ)                 异性明显强于二阶模型,表明裂缝之间的相互作用
                      
                                                        (4)
                              4(λ + 2µ)
                                                              有降低准纵波、准SV波各向异性强度的作用。当裂
                       ¯              .
                      U 33 =
                               3(λ + µ)                        缝密度为 0.2时,一阶模型波场快照出现畸变。同时
             1.2 各向异性参数                                        随着裂缝密度增大,准SV波能量增强。
                 为了解裂缝存在而对各向异性特征的影响,本
                                                                          表 1  裂缝介质模型参数取值
             文利用 Thomsen    [17]  提出的 3 个无量纲的各向异性
                                                                  Table 1    Parameters for the fracture
             参数对裂缝诱导各向异性进行表征,Thomsen 系数
                                                                  medium model
             将物理意义不直观的弹性常数进行组合,简化了相
             速度公式的表达式,参数表达式为                                            组成         背景基质       水     天然气
                   
                         c 11 − c 33                                         −1
                   ε =          ,                                 纵波速度/(m·s   )     5800    1500    447
                   
                   
                          2c 33
                   
                                                                  横波速度/(m·s −1 )    3349     0       0
                   
                        c 66 − c 44
                     γ =         ,                      (5)         密度/(kg·m −3 )    2600    1000    0.65
                           2c 44
                   
                   
                                  2            2
                                                                   体积模量/GPa         87.46   2.25  0.00013
                        (c 13 + c 44 ) − (c 33 − c 44 )
                   
                   δ =                          ,
                   
                              2c 33 (c 33 − c 44 )                  剪切模量/GPa         29.16    0       0
             其中,ε 表征纵波各向异性,γ 表征横波各向异性,δ
                                                                   当裂缝饱含水时,Hudson 一阶模型和二阶模
             为纵波与横波之间的过渡参数。在弱各向异性介质
                                                   √           型在裂缝密度分别为0.05、0.1、0.15、0.2时的波场快
             中,采用3个各向异性参数和纵波速度α =                    c 33 /ρ、
                          √                                    照如图 2 所示。对于饱含水裂缝,随着裂缝密度增
             横波速度 β =       c 44 /ρ 便可将 3 种波的相速度表达
                                                               大,准纵波各向异性无明显变化,准 SV波的能量和
             式简化为
                                                               各向异性逐渐增强。对比一阶和二阶模型的波场快
               
                          (        2             4  )
                                        2
               V P (θ) ≈ α 1 + δ sin θ cos θ + ε sin θ ,      照,含水裂缝的一阶模型和二阶模型的波场快照无
               
               
               
               
                                                              明显区别。
                          (      2                 )
                                α           2    2
                 V SV (θ) ≈ β 1 +  (ε − δ) sin θ cos θ ,           对比图 1 和图 2 可以分析充填流体对波场传播
                               β 2
               
               
               
                                                              的影响。与干燥裂缝相比,裂缝含水时准纵波、准
                          (        2  )
                 V SH (θ) ≈ β 1 + γ sin θ .
               
                                                               SV波的能量比裂缝为干燥裂缝时弱,同时准纵波各
                                                        (6)
                                                               向异性减弱。
             1.3 数值模拟分析                                            为进一步分析 Hudson 一阶模型和二阶模型波
                 为了解在干燥裂缝和含水裂缝两种情况下裂                           场快照出现差异的原因,分别考虑裂缝密度和裂缝
             缝密度和裂缝充填状况对 Hudson 一阶模型和二                         倾角对弹性常数和各向异性参数的影响,为后期裂
             阶模型的地震波传播特征影响,建立了一个裂缝                             缝储层微观建模明确应用条件。
             介质模型,x、z 方向网格数为 500×500,网格间距                          由图 3(a) 可以看出,对于干燥裂缝,当裂缝密
             ∆x = ∆z = 2 m,震源采用集中力源,中心频率                       度小于 0.1 时,一阶与二阶模型基本表现良好。当
             60 Hz,位于模型中心位置。裂缝在模型中水平排列                         裂缝密度大于 0.1 时,一阶模型和二阶模型的差异
             分布,纵横比为 0.001,裂缝介质模型的背景基质和                        逐渐增大,一阶模型在裂缝密度接近 0.2 时c 33 甚至
             流体相关参数引自文献[18],如表1 所示。采用有限                        会出现负值。需要得到特别注意的是,虽然二阶模
             差分数值模拟方法,模拟地震波在裂缝介质中的传                            型在一定程度上扩展了裂缝密度的适用范围,但是
             播,本文以0.1 s时的波场快照为例进行对比分析。                         在裂缝密度较大时弹性系数随裂缝密度增大不减
                 当裂缝为干燥裂缝时,Hudson 一阶模型和二                       反增,即二阶展开并不收敛,违背了物理规律,在
             阶模型在裂缝密度分别为0.05、0.1、0.15、0.2时的波                   这种情况下只使用一阶更正比不恰当地使用二阶
   92   93   94   95   96   97   98   99   100   101   102