渗透率各向异性条件下裂缝介质地震频散、衰减及各向异性分析.pdf

2023 年 8 月第 58 卷 第 4 期渗透率各向异性条件下裂缝介质地震频散、衰减及各向异性分析姜晓敏1，2，陈学华*1，2，吴昊杰2，罗鑫1，刘俊杰2（1.成都理工大学油气藏地质及开发工程国家重点实验室，四川成都 610059；2.成都理工大学地球勘探与信息技术教育部重点实验室，四川成都 610059）摘要：渗透率各向异性与地震各向异性是含流体裂缝孔隙岩石的关键特性，文中利用裂缝参数（倾角、方位角、半径等）表征裂缝介质的渗透率并分析渗透率各向异性变化规律，基于 Chapman多组裂缝等效介质模型研究了渗透率各向异性条件下正交介质的地震各向异性及地震波频散、衰减特征。数值模拟结果表明：渗透率各向异性与裂缝参数有关，裂缝方位角的变化对渗透率各向异性的影响具有对称性；裂缝半径和倾角的变化影响渗透率各向异性的强弱，其中裂缝半径的影响尤为显著。渗透率各向异性情况下，地震波的频散、衰减出现两个不同的弛豫时间，随着饱含流体裂缝介质中含气量的增加，特征频率呈现先减后增的趋势。渗透率差异引起的两个不同的特征频率与地震各向异性系数极值对应的频率具有相关性，表明渗透率各向异性会引起地震频散、衰减的各向异性。数值分析结果可为裂缝储层渗透率预测及流体识别提供有力的理论依据。关键词：裂缝参数，渗透率各向异性，地震各向异性，频散，衰减中图分类号：P631 文献标志码：A doi：10.13810/ki.issn.1000-7210.2023.04.011Analysis of seismic wave dispersion，attenuation and anisotropyof fractured media under anisotropic permeabilityJIANG Xiaomin1，2，CHEN Xuehua1，2，WU Haojie2，LUO Xin1，LIU Junjie2（1.State Key Laboratory of Oil&Gas Reservoir Geology and Exploitation，Chengdu University of Technology，Chengdu，Sichuan 610059，China；2.Key Lab of Earth Exploration&Information Techniques of Ministryof Education，Chengdu University of Technology，Chengdu，Sichuan 610059，China）Abstract：Permeability anisotropy and seismic anisotropy are the key characteristics of fluidbearing fracturepore rock.In this paper，fracture parameters（dip angle，azimuth angle，radius，etc.）are used to characterize the permeability of fractured media and analyze the anisotropy variation of permeability.The seismic anisotropy and dispersionattenuation characteristics of seismic waves of orthogonal media under permeability anisotropy are studied based on Chapman s multigroup fracture equivalent medium model.The numerical simulation results show that：The permeability anisotropy changes with fracture parameters and is symmetrically influenced by fracture azimuth angle；the variation of fracture radius and dip angle affects the permeability anisotropy，especially the fracture radius.In the case of permeability anisotropy，there are two different relaxation time for the dispersion and attenuation of seismic waves，and the characteristic frequency decreases first and then increases with the increase in gas content in the fluidbearing fractured medium.The two different characteristic frequencies produced by different permeability are related to the frequency corresponding to the extreme value of the seismic anisotropy coefficient.This indicates that permeability anisotropy can cause anisotropy of seismic dispersion and attenuation.The numerical analysis results in this paper can provide a powerful theoretical basis for the prediction of reservoir permeability and fluid identification in fractured reservoirs.Keywords：fracture parameters，permeability anisotropy，seismic anisotropy，dispersion，attenuation姜晓敏，陈学华，吴昊杰，等.渗透率各向异性条件下裂缝介质地震频散、衰减及各向异性分析 J.石油地球物 岩石物理 文章编号：1000-7210（2023）04-0902-11*四川省成都市成华区二仙桥东三路 1号成都理工大学地球物理学院 5301室，610059。Email：chen_ 本文于 2022年 8月 22日收到，最终修改稿于 2023年 4月 15日收到。本项研究受国家自然科学基金项目“致密储层裂缝系统诱发地震异常的机理及其与储层产能的关系”（41874143）和四川省自然科学基金重 点项目“深层储层裂缝渗流结构的地震岩石物理响应机理与检测方法”（2023NSFSC0019）联合资助。第 58 卷 第 4 期姜晓敏，等：渗透率各向异性条件下裂缝介质地震频散、衰减及各向异性分析理勘探，2023，58（4）：902912.JIANG Xiaomin，CHEN Xuehua，WU Haojie，et al.Analysis of seismic wave dispersion，attenuation and aniso tropy of fractured media under anisotropic permeability J.Oil Geophysical Prospecting，2023，58（4）：902912.0引言含流体裂缝孔隙岩石的地震各向异性及地震波频散、衰减特性对储层渗流特征提取及流体识别具有指导意义。当地震波穿过饱含流体的裂缝介质时，流体流动会影响地震各向异性特征，同时伴随着地震波的能量耗散。在裂缝介质中，孔隙单元之间的流体流动是引起地震频带内地震波频散、衰减的主要原因1，而流体的流动依赖于介质渗透率。因此，研究裂缝孔隙岩石渗透率的各向异性与地震各向异性间的影响机理具有重要意义。为了研究中尺度裂缝引起的地震频带内的强衰减特性，Chapman2 建立了包含一组定向排列的中尺度裂缝的局部喷射流模型，综合考虑流体性质、裂缝尺度等因素，分析并解释了裂缝介质中地震波的频 散、衰减现象，结果表明裂缝尺度是影响地震各向异性的关键因素；Brajanovski等3 类比线性滑动理论4，通过周期性排列的高孔隙度薄层等效裂缝介质，分析裂缝与孔隙之间流体交换引起的弹性波频散和衰减；2009年，Chapman5 将一组定向排列裂缝扩展到两组裂缝，通过控制裂缝闭合情况，分析了裂缝方位角对地震波频散、衰减和地震各向异性的影响。在上述研究中，岩石物理模型中的渗透率均为标量（各向同性），储层流体饱和仅考虑单一流体和均匀饱和的情况，而实际储层中多相流体及非均饱和情况普遍存在，裂缝介质的渗透率多表现为各向异性，并且可通过等效渗透率进行描述68。为了进一步了解裂缝储层中地震波的频散、衰减特性，学者们对多相流体等效黏滞度9及流体饱和情况1013开展了广泛研究。Snow1415 采用平行平板对裂缝进行等效建模，提出了包含任意数量、方向及孔径的裂缝介质等效模型，将裂缝与基质的渗透率线性叠加获得裂缝介质的等效渗透率张量；Long 等1617 提出的离散裂缝网（DFN）模型，利用统计数据模拟裂缝系统，但该模型忽略了基质固有渗透率，无法描述裂缝与基质间流体交换过程中的地震波能量耗散；Vu等18 假设裂缝受到无限远处的均匀压力梯度，基于 Darcy定律和 Eshelby等效夹杂理论1921 推导包含多组裂缝的裂缝介质等效渗透率表达式；Xiong等22 建立三维孔隙网络模型，探讨了流体性质和孔隙连通性对岩石渗透率的影响；熊繁升等23 提出三维裂缝网络模型，定量分析了岩石渗透率随裂缝参数（发育密度、纵横比）、孔隙流体类型和围压等因素的变化规律。上述研究表明，裂缝的半径、相对位置及发育密度等因素会使裂缝介质的渗透率表现出各向异性。诸多学者还研究了不同渗流条件下地震波频 散、衰减特性。Pride等24 研究表明地震数据中的振幅可能包含与渗透率相关的信息；Kozlov25 通过地震数据频谱分析，明确了地震频谱中包含渗透率信息；Rubino等26 通过改变储层渗透率的大小，分析渗透率对气水饱和储层中地震波频散、衰减的影响，结果表明随着储层渗透率增加地震波速度从高频极限不断下降到低频极限；Ren等27 建立周期性裂缝储层模型，模型的渗透率表现为纵向各向异性，分析了渗透率对含气储层地震反射振幅的影响；Mller等28 通过推导动态等效渗透率模型分析渗透率的空间变化规律，并且探讨了渗透率对地震波频散、衰减的影响，发现渗透率的变化会引起弛豫峰值沿频率轴移动；Rubino等2930 通过分析裂缝的连通性及其控制裂缝介质的各向异性渗流情况，研究了渗透率对地震波频散、衰减的影响，结果表明裂缝介质渗透率的变化会引起弛豫时间的差异。王丁等31 在非均匀裂隙孔隙介质中建立各向异性渗流模型，研究各向异性渗流条件下弹性波传播规律，结果表明渗透率的变化使地震波衰减曲线产生第二个弛豫峰。上述研究缺乏合适的介观尺度裂缝模型探讨裂缝参数变化导致的渗透率各向异性，同时，缺乏对渗透率各向异性条件下地震各向异性变化特征的研究。为此，本文首先通过裂缝参数表征裂缝介质的等效渗透率，分析裂缝参数变化引起的裂缝介质中渗透率的各向异性；然后基于 Chapman多组裂缝岩石物理模型5，研究各向异性渗流条件下，不同类型流体饱和的正交介质中地震波频散、衰减与地震各向异性的变化规律，探讨渗透率各向异性与地震各向异性间的响应机理，为由地震各向异性提取实际储层渗流特征及流体识别提供理论依据。903石 油 地 球 物 理 勘 探2023 年1方法原理1.1裂缝介质渗透率表征基于Vu等18 提出的裂缝介质等效渗透率半解析表征方法，综合考虑背景介质的渗透率、裂缝方向及半径等裂缝参数，对裂缝介质渗透率进行表征。假设含多组椭圆形裂缝的裂缝介质受到远场均匀压力，则第m个裂缝组的等效渗透率fm表征为fm=m16cmR3m3cm+163Rm(E-nm nTm)A（1）式中：表示各向同性介质的渗透率；cm为导水率；m为裂缝发育密度（单位体积中裂缝中心的数量）；Rm为裂缝半径（椭圆形裂缝的半长轴，见图 1）；nm为裂缝的法向向量；A表示远场压力梯度；E为单位矩阵；表示两个向量的张量积；“f”代表与裂缝相关的变量。如图1所示，假设第m个裂缝组的法向量nm与观测坐标系(x1，x2，x3)的夹角分别为m、m，则nm=(cos msin m，sin msin m，cos m)（2）由此，包含N个裂缝组的裂缝介质的等效渗透率可表示为eff=E+m=1Nfm（3）假 设 裂 缝 介 质 的 主 渗 透 方 向 与 观 测 坐 标 系(x1，x2，x3)重合，含多组裂缝的裂缝介质等效渗透率张量可写成由不同平面的渗透率张量元素组成的矩阵形式eff=111213122223132333（4）由此，裂缝介质不同平面内的渗透率各向异性可通过三个主渗透方向的渗透率比值来评估，即ij=|jjii-1|（5）式中：i=1，2；j=2，3，且i j；ij表示平面xi，xj内的渗透率各向异性。1.2正交介质中频率依赖的速度及各向异性系数Chapman5 将含有一组对齐裂缝的等效介质模型拓展到两组裂缝的情况，正交介质频率依赖的等效刚度矩阵由背景介质刚度矩阵、球形孔隙校正量和两个裂缝组校正量的总和给出，即Cij()=C0ij-Cpij()-m=12Cmij()（6）式中：C0ij表示由拉梅系数计算得到的各向同性刚度矩阵；Cpij代表球形孔隙的刚度矩阵校正量；Cmij代表第m个裂缝组的刚度矩阵校正量；为角频率。考虑裂缝介质渗透率各向异性情况下，第m个裂缝组与孔隙间流体交换的弛豫时间m为m=8l()1-()1+Kc3mRm（7）其中m=211+222+2333式中：l为颗粒尺度；为剪切模量；为泊松比；为流体黏滞度；Kc为流体体积模量相关系数。通过弛豫时间可计算出不同裂缝组中的流体压力，考虑裂缝孔隙间流体质量守恒并运用 Eshelby等效夹杂理论21，获得了各向异性渗流条件下正交介质频率依赖的等效刚度矩阵，即C11C12C13000C12C22C23000C13C23C33000000C44000000C55000000C66（8）根据上述等效刚度矩阵，垂直入射的地震波纵波复速度计算公式32 为v()=C33()（9）式中为裂缝介质密度。由此，根据复速度 v 计算频率依赖的纵波速度VP()及衰减系数1 QP()表达式26 为图 1椭圆形裂缝平面法向量投影示意图904第 58 卷 第 4 期姜晓敏，等：渗透率各向异性条件下裂缝介质地震频散、衰减及各向异性分析VP()=1Re 1v()1QP()=Imv2()Rev2()（10）式中Re 和Im 分别表示对复数求实部和虚部。Thomsen32 通过等效刚度矩阵计算各向异性系数以刻画各向异性介质的地震各向异性特征，Tsvankin33 在此基础上推导了正交介质的各向异性系数。据此，得到由等效刚度矩阵元素表示的正交介质各个平面内描述流体流动的地震各向异性系数。（1）在 x1，x3 平面内1=C11-C332C331=()C13+C552-()C33-C5522C33()C33-C55（11）（2）在 x2，x3 平面内2=C22-C332C332=()C23+C442-()C33-C4422C33()C33-C44（12）（3）在 x1，x2 平面内3=()C12+C662-()C11-C6622C11()C11-C66（13）式中、为各向异性参数，下标代表不同平面。2数值模拟分析2.1裂缝介质渗透率各向异性以下基于不同裂缝参数进行裂缝介质的渗透率表征及渗透率各向异性特征分析。假设远场均匀压力梯度与裂缝介质渗透率的主渗透方向重合，各向同性介质渗透率为20 mD，两组裂缝的基本参数如表1所示。2.1.1渗透率表征为了研究裂缝参数变化对裂缝介质渗透率张量的影响，参照表1调整裂缝组1的参数变化范围，并展示了当裂缝半径 R1=R2=0.50 m 时，不同倾角和方位角对应的等效渗透率张量（图 2），边框数值代表等效渗透率的元素下标。图2a、图2b分别为1=45时，不同裂缝倾角对应的等效渗透率张量eff和去除元素11、22、33后的渗透率张量eff-d；图 2c、图 2d 分别为 1=30时，不同裂缝方位角对应的eff和eff-d。由图 2a、图 2c可见，不同裂缝倾角和方位角对应的等效渗透率张量几乎一致；由图 2b、图 2d可见，不同裂缝倾角和方位角对应的eff-d 各不相同。这表明当裂缝介质受到同一均匀压力时，不同倾向和方位的裂缝对裂缝介质的渗透率贡献存在差异。图 3显示了主渗透方向渗透率随 1和1变化曲线，可更清晰地展示裂缝参数对裂缝介质等效渗透率张量的影响。由图3a可见，随着1增大，11、22逐渐减小，且11、22间的差异逐渐增大；不同方位角33对应曲线完全重合且逐渐增大；1=45时11、22的对应曲线重合。由图 3b可见，随着1增大，11逐渐增大，22逐渐减小，11、22间的差异先减小后增大；33不随方位角变化；1=45时 11、22相等。此外，图 3b还显示随裂缝方位角的增大，11、22、33的变化趋势具有对称性且以1=45为对称轴对称变化，这表明在一定范围内 （090），裂缝倾角的增加会减弱x1、x2方向的渗流能力，同时增强 x3方向的渗流能力；裂缝方位角的增加会增强x1方向的渗流能力，同时减弱x2方向的渗流能力，而 x3方向的渗流能力不改变，且 x1、x2方向渗流能力增强或减弱趋势互相对称。为了探究其他裂缝参数对裂缝介质渗透率的影响，绘制了1=45时，主渗透方向渗透率随R1及1变化曲线（图4）。由图可见，随着R1增加，11、22、33非线性增加，且 11、22的变化较 33更剧烈，变化范围从 20 mD到160 mD，这种变化跨越一个数量级，表明裂缝半径对裂缝介质渗透率的影响是极为显著的。图4a中33随R1变化的三条曲线完全重合；在1=45时，11、22随 R1变化的曲线完全重合；此外，图 4b 和 图3b中11、22、33随1的变化趋势一致，表明11、22、33随裂缝方位角的变化具有对称性。裂缝参数的变化引起11、22、33间的变化特征存在差异，表明裂缝参数的变化会导致裂缝介质渗透率表现为各向异性。表 1裂缝参数参数类型纵横比发育密度半径 R/m倾角/()方位角/()裂缝组 11 10-40.050.011.50090090裂缝组 21 10-40.050.5000905石 油 地 球 物 理 勘 探2023 年图 2不同裂缝方位角、裂缝倾角对应的不同渗透率张量（a）1=45，eff；（b）1=45，eff-d （c）1=30，eff；（b）1=30，eff-d 图 3裂缝介质主渗透方向渗透率随1（a）和1（b）的变化曲线906第 58 卷 第 4 期姜晓敏，等：渗透率各向异性条件下裂缝介质地震频散、衰减及各向异性分析2.1.2渗透率各向异性图5显示由于1和1变化引起的不同平面内裂缝介质渗透率各向异性平面图。由图可见，13、23随 1增大而减弱，13随1增大而增强，23随1的变化与13相反；12随 1增大而增强，12随1增大而先减弱后增强，并以1=45为对称轴对称分布。这表明裂缝参数变化造成裂缝介质在不同平面内渗透率各向异性的强弱变化各有特点。在 x1，x3 平面和 x2，x3 平面内，裂缝倾角的变化引起渗透率各向异性的强弱变化相对剧烈；而在 x1，x2 平面内，裂缝方位角的变化引起渗透率各向异性的强弱变化具有对称性。图6展示了1=45时，1及R1的变化引起裂缝介质不同平面内渗透率各向异性平面图。从图中可见，1较小时（135），13随R1增大而增强，23的变化特征与13相反。12随R1增大而增强，不同的1增强幅度不同，且12的强弱变化仍是对称的。结合图 5、图 6 可知，裂缝介质渗透率的各向异性是裂缝参数共同作用的结果，裂缝方位角的变化主要影响渗透率各向异性变化的对称性；裂缝半径、倾角的变化主要影响渗透率各向异性的强弱变化，其中裂缝半径的影响尤为显著。图 41=45裂缝介质主渗透方向渗透率随R1（a）和1（b）的变化曲线图 5裂缝倾角、方位角变化引起不同平面内渗透率各向异性平面图（a）x1，x3 平面；（b）x2，x3 平面；（c）x1，x2 平面图 6裂缝半径、方位角变化引起不同平面内渗透率各向异性平面图（a）x1，x3 平面；（b）x2，x3 平面；（c）x1，x2 平面907石 油 地 球 物 理 勘 探2023 年2.2各向异性渗流时地震频散、衰减与各向异性为了研究渗透率各向异性情况下，裂缝介质的地震各向异性变化规律和地震波的频散衰减特征，本文基于 Chapman 多组裂缝等效介质模型5设计了双相流体饱和的正交介质模型进行数值模拟。设置裂缝组 1的=0、=0；裂缝组 2的=90、=0，其他裂缝参数见表 1，背景介质参数见表 2。对于裂缝介质的流体饱和情况，本文基于 Jin等12 提出的各向异性介质的双相流体饱和模型，分别考虑裂缝介质为气水饱和及油水饱和两种情况（流体物理参数见表3）。此外，岩石中填充流体的基本参数由Batzle等34 的方程在压力为25MPa、温度为50条件下计算得到。2.2.1频率依赖的速度及衰减图7为岩石在不混溶双相流体饱和情况下，地震波速度及衰减系数随频率、含气（油）饱和度变化平面图。由图可见，无论气水饱和岩石还是油水饱和岩石，地震波衰减平面图中（图 7b）均出现了两个弛豫峰，但因饱含流体类型的不同，地震波速度频散、衰减的特征频率（地震波速度变化速率最高点或衰减系数的极值对应的频率）变化趋势不同，衰减系数的强弱也有区别。图 8更清晰地分析了地震波速度和衰减系数随饱和度和频率的变化。由图 8可见：地震波速度出现两个频段内速度剧烈变化的情况，且衰减系数曲线相应出现两个弛豫峰，对应的特征频率为 1和 2，这是由于岩石渗透率的各向异性而引起的。地震波经过岩石会引起孔隙流体诱导压力，岩石中流体运移能力强的位置平衡该诱导压力所需时间短、特征频率高；反之，则时间长、特征频率低，故地震波衰减曲线有低频段特征频率 1和高频段特征频率 2。由于气体与油的黏滞度存在差异，随着饱和度增加，气水饱和岩石对应衰减系数曲线的弛豫峰值增加且向高频方向移动；油水饱和岩石对应衰减系数曲线弛豫峰值增加且向低频方向移动。此外，由于水相与气相间湿润性的差异远大于水相与油相间的差异，在含气饱和度与含油饱和度相等情况下，气水饱和岩石中的地震波频散范围更宽、衰减更强，这是因为地震波经过气水饱和岩石时会产生更强的流体诱导压力，引起更多的流体流动，故气水饱和岩石中地震波衰减更强（图 8a、图 8b）。不同频率的地震波速度对含气（油）饱和度的敏感性不同。随着饱和度的增加，不同频率的地震波速度均呈非线性降低趋势，降低的速率与饱和度相关；气水饱和岩石的衰减系数曲线呈现先增加后减小的趋势，油水饱和岩表 2背景介质参数物理参数岩石拉梅系数/Pa岩石拉梅系数/Pa密度/(kgm3)孔隙度/(%)颗粒尺度/m基质渗透率/mD数值2.58 10109.9 1092750100.00120表 3流体参数流体参数密度/(kgm 3)体积模量/MPa黏滞度/cP油749.91248.81.2687天然气201.660.60.0254盐水1068.23024.40.7550图 7双相流体饱和岩石随含气（油）饱和度 和频率变化的速度（a）、衰减（b）平面图908第 58 卷 第 4 期姜晓敏，等：渗透率各向异性条件下裂缝介质地震频散、衰减及各向异性分析石的衰减系数曲线呈现单调增加趋势，不同频率的衰减系数随饱和度的变化速率各不相同，这种差异是岩石渗透率各向异性与流体饱和度共同作用的结果 （图8c、图8d）。2.2.2频率依赖的地震各向异性通过各向异性系数随频率和含气（油）饱和度的变化曲线（图9、图10）分析各向异性渗流条件下，正交介质中不同平面内的地震各向异性变化规律。从 图 9、图 10可见：受渗透率各向异性与流体饱和度的共同影响，裂缝介质不同平面内地震各向异性变化特征存在差异。随着频率增加，各向异性系数 非线性单调减小，不同频带内减小速率不同；各向异性系数1呈先增加后减小的趋势，2呈非线性减小的趋势。随着饱和度增加，各向异性系数和非线性增加，且不同平面内增加幅度、速率各不相同。这表明岩石中流体渗流的差异会导致岩石的地震各向异性特征发生变化。另外，由于两组裂缝的裂缝参数不同，岩石不同平面内的渗透率各向异性表现为132312，使得岩石各向异性系数表现为321，各向异性系数表现为12。由此可知，渗透率各向异性的强弱与岩石各向异性系数相反、与一致，表明裂缝介质渗透率的各向异性与地震各向异性具有一定相关性。由于气相与水相间的湿润性差异较大而产生相对强的诱导压力，进而引起大量流体定向流动。因此，同一个平面内气水饱和岩石的各向异性系数大于油水饱和岩石。这与页岩中裂隙与粘土矿物的定向排列导致页岩强烈的各向异性相类似。结合图8图10中3、2随频率变化的曲线可知，3、2在两个频段内剧烈减小，并且这两个频带将特征频率1和2包含在内。此外，还可见1随频率变化曲线的极值对应频率（3）大致位于3和2剧烈变化的两个频带的分界处。通过图 11显示双相流体饱和岩石中1、2及3随饱和度的变化曲线，以明确渗透率各向异性情况下特征频率与地震各向异性系数间的关系。由图可见：由于渗透率的各向异性，地震波速度频散衰减产生两个特征频率1和2，其中1集中在低频段（0100 Hz），2集中在高频段（3001 kHz）。随着含气饱和度增加，1、2呈现先减小后增加的趋势；随着含油饱和度增加，1和2呈现先增加后减小的趋势，这种差异与岩石中饱和流体的湿润性相关。含气饱和度越高，1和 2越高且跨越频带越宽，表明特征频率对高含气储层具有高敏感性。在同一流体饱图 8双相流体饱和岩石中纵波速度和衰减系数随频率、饱和度的变化曲线（a）纵波速度与频率；（b）衰减系数与频率；（c）纵波速度与饱和度；（d）衰减系数与饱和度909石 油 地 球 物 理 勘 探2023 年和岩石中，3随饱和度的变化趋势与1、2一致，且介于1和2之间，甚至是平均值。这表明含流体裂缝孔隙岩石的地震各向异性变化规律与岩石的渗透率各向异性情况相关，为实际裂缝储层中通过地震各向异性特征与地震波频散衰减特性预测储层渗流特征及流体识别提供了理论基础。图 9气水饱和岩石中 和随频率、饱和度的变化曲线（a）与频率；（b）与含气饱和度；（c）与频率；（d）与含气饱和度图 10油水饱和岩石中 和随频率、饱和度的变化曲线（a）与频率；（b）与含油饱和度；（c）与频率；（d）与含油饱和度910第 58 卷 第 4 期姜晓敏，等：渗透率各向异性条件下裂缝介质地震频散、衰减及各向异性分析3结论本文分析了裂缝参数变化引起的裂缝介质渗透率各向异性变化规律，探讨了各向异性渗流条件下，双相不混溶流体饱和的裂缝介质地震波频散、衰减和地震各向异性特征以及渗透率各向异性与地震各向异性的相关性，得到如下结论。（1）渗透率随裂缝参数（半径、倾向及方位）的变化规律各不相同，进而导致裂缝介质渗透率表现出各向异性。裂缝方位角引起的渗透率各向异性变化具有对称性；裂缝半径和倾角主要引起渗透率各向异性的强弱变化，其中裂缝半径的影响尤为显著。（2）渗透率各向异性条件下，地震波速度频散、衰减曲线呈现两个弛豫峰。同时，随饱和度增加，地震波速度的衰减增强，其中衰减系数对含气饱和度的变化极为敏感。（3）裂缝介质渗透率各向异性的强弱与地震各向异性系数相关，同时地震波速度频散、衰减的特征频率与各向异性系数极值或各向异性系数变化率的极值所对应的频率密切相关，这种相关性依赖于渗透率的各向异性情况、流体类型及流体饱和度。参 考 文 献1 MLLER T M,GUREVICH B,LEBEDEV M.Seismic wave attenuation and dispersion resulting from waveinduced flow in porous rocks:A reviewJ.Geophysics,2010,75(5):75A14775A164.2 CHAPMAN M.Frequencydependent anisotropy due to mesoscale fractures in the presence of equant poro sityJ.Geophysical Prospecting,2003,51(5):369379.3 BRAJANOVSKI M,GUREVICH B,SCHOEN BERG M.A model for Pwave attenuation and dispersion in a porous medium permeated by aligned fracturesJ.Geophysical Journal International,2005,163(1):372384.4 SCHOENBERG M.Elastic wave behavior across linear slip interfacesJ.The Journal of the Acoustical Society of America,1980,68(5):15161521.5 CHAPMAN M.Modeling the effect of multiple sets of mesoscale fractures in porous rock on frequencydependent anisotropyJ.Geophysics,2009,74(6):D97D103.6 ODA M.Permeability tensor for discontinuous rock massesJ.Gotechnique,1985,35(4):483495.7 KACHANOV M.Continuum model of medium with cracksJ.Journal of the Engineering Mechanics Divi sion,1980,106(5):10391051.8 JAKOBSEN M.Effective hydraulic properties of frac tured reservoirs and composite porous mediaJ.Journal of Seismic Exploration,2007,16(2):199224.9 CHEN X H,ZHONG W L,HE Z H,et al.Frequency dependent attenuation of compressional wave and seismic effects in porous reservoirs saturated with multiphase fluidsJ.Journal of Petroleum Science and Engineering,2016,147:371380.10 PAPAGEORGIOU G,AMALOKWU K,CHAP MAN M.Theoretical derivation of a Brielike fluid mixing lawJ.Geophysical Prospecting,2016,64(4):10481053.11 PAPAGEORGIOU G,CHAPMAN M.Wavepropa gation in rocks saturated by two immiscible fluidsJ.Geophysical Journal International,2017,209(3):17611767.12 JIN Z Y,CHAPMAN M,PAPAGEORGIOU G.Frequencydependent anisotropy in a partially saturated fractured rockJ.Geophysical Journal International,2018,215(3):19851998.13 刘芸菲,陈学华,罗鑫,等.双相不混溶流体饱和裂缝孔隙岩石依赖频率的地震响应数值分析J.石油地球物理勘探,2020,55(4):821830.LIU Yunfei,CHEN Xuehua,LUO Xin,et al.Numerical analysis of frequencydependent seismic responses from fracturedporous rock saturated with twophase immiscible fluidsJ.Oil Geophysical Prospecting,2020,55(4):821830.14 SNOW D T.Anisotropie permeability of fractured mediaJ.Water Resources Research,1969,5(6):1273 1289.15 SNOW D T.The frequency and apertures of fractures in rockJ.International Journal of Rock Mechanics and Mining Sciences&Geomechanics Abstracts,1970,7(1):2340.16 LONG J C,REMER J S,WILSON C R,et al.Porous media equivalents for networks of discontinuous fracturesJ.Water Resources Research,1982,18(3):645 658.17 LONG J C,GILMOUR P,WITHERSPOON P A.A model for steady fluid flow in random threedimen sional networks of discshaped fracturesJ.Water Resources Research,1985,21(8):11051115.18 VU M N,POUYA A,SEYEDI D M.Modelling of 图 11双相流体饱和岩石随饱和度变化的特征频率曲线红色曲线代表气水饱和岩石，绿色曲线代表油水饱和岩石；实线代表1，虚线表示2，点线代表3。911石 油 地 球 物 理 勘 探2023 年steadystate fluid flow in 3D fractured isotropic porous media:application to effective permeability calculationJ.International Journal for Numerical and Analytical Methods in Geomechanics,2013,37(14):22572277.19 SHAFIRO B,KACHANOV M.Anisotropic effective conductivity of materials with nonrandomly oriented inclusions of diverse ellipsoidal shapesJ.Journal of Applied Physics,2000,87(12):85618569.20 BARTHLMY J F.Effective permeability of media with a dense network