岩石中孔隙和裂隙的三维cfd模拟.docxVIP

岩石中孔隙和裂隙的三维cfd模拟 众所周知，介质的孔隙度和体积指数是影响地震反射波振幅和位移之间变化的重要因素。长期以来，通过振幅随位移变化及其avo属性的相交图进行储层预测和油藏规划，取得了良好的效果。尤其是对于具有较差固结不良祖母岩储量的破碎屑岩储量，有许多成功的例子。然而,随着油气勘探向复杂地区延伸,勘探目标越来越复杂,特别是在裂缝发育的地区,AVO应用受到了挑战,主要困难是难以确定岩石物性与AVO响应的关系。笔者通过裂缝介质的AVO正演模拟,研究裂缝对AVO的影响,并进一步分析对AVO交会图识别气层的影响。 AVO的地质基础在于不同岩石以及含有不同流体的同类岩石之间泊松比的差异,故在做AVO正演模拟时,泊松比与孔隙度、流体饱和度、孔隙纵横比变化关系的求取变得十分重要,利用Berryman散射模型方法模拟孔隙介质,通过设定不同的孔隙几何形状(孔隙纵横比)及其分布,来模拟裂隙分布。通过对一个两层砂泥岩介质模型的AVO正演模拟(保持孔隙度和流体饱和度不变,只改变孔隙纵横比),研究反射系数与孔隙纵横比的关系。根据Gassmann 方程,对孔隙流体进行替换,研究孔隙纵横比对AVO交会图识别气层的影响。 1 方法原理 1.1 表征及基本参数 把岩石当作是由不同形状的孔隙、各种矿物成分及孔隙流体组成的复合介质,矿物成分的百分含量可任意设定(所有矿物成分百分含量的总和为1),矿物颗粒形状及孔隙形状是可变的。Berryman根据长波长的一次散射理论,即弹性波波长远大于组成复合介质矿物颗粒和孔隙尺寸,此时可忽略高次散射,推导出了复合介质的等效密度ρm、等效体积弹性模量Κm和等效剪切模量μm与各组分的密度ρi、体积弹性模量Κi和剪切模量μi的关系 n∑i=1ci(ρi-ρm)=0,n∑i=1ci(Κi-Κm)Ρi=0,n∑i=1ci(μi-μm)Qi=0。(1)∑i=1nci(ρi?ρm)=0,∑i=1nci(Ki?Km)Pi=0,∑i=1nci(μi?μm)Qi=0。(1) 式中,n为组分数,ci为第i种介质的体积与复合介质的体积百分比,Pi和Qi是与组分和孔隙形状有关的参数,Berryman给出了球状、针状、碟状和椭球状4种形状的Pi和Qi表达式。笔者用椭球状孔隙模拟裂隙,对于椭球状孔隙,其表达式如下 Ρi=Κm+43μiΚi+43μi+παβm?(2)Pi=Km+43μiKi+43μi+παβm?(2) Qi=15[1+8μm4μi+πα(μm+2βm)+2Κi+23μi+23μmΚi+34μi+παβm]?(3) 1.2 含气砂岩的avo反演 AVO分析技术在实际应用中乃是利用地震反射的共成像点(CIP)道集资料分析反射振幅随炮检距的变化规律,或通过计算纵波反射振幅随入射角的变化参数,估算界面上的AVO属性参数和泊松比差,进一步推断储层岩性和含油气性质。AVO正演模拟是采用AVO方法进行烃类检测的基础,模拟的结果可以指导利用实际的地震道集的AVO反演结果进行可靠的含气砂岩解释,AVO主要的理论基础是Zoeppritz(1919)方程,对于各向同性介质,平面纵波以某一角度入射到某一平界面上,发生的反射和透射及波的转换可由Zoeppritz方程计算得到。由于Zoeppritz方程的复杂性,从应用的角度考虑,不少研究者从不同方面对其进行了简化(例如Aki等,1980;Bortfeld ,1961;Shuey,1985)。Shuey在基于一定假设的基础上,推导出了PP(纵波入射,纵波反射)反射系数随入射角变化的关系 R(x)=R0+[A0R0+Δσ(1-σ)2]sin2α+Δvp2vp(tan2α-sin2α)?A0=B-2(1-B)(1-2σ)1-σ?B=Δvpρ(Δvpρ+Δρvp)。(4) 式中,vp为相邻2层纵波速度的平均, Δvp为相邻2层的纵波速度差,R0为垂直入射时的反射系数,σ为相邻2层的泊松比的平均,Δσ为相邻2层的泊松比差,α为入射角。该公式是目前应用最多的Zoeppritz近似方程,笔者采用该公式计算AVO响应。 2 计算与分析 2.1 孔隙结构及厚度 依据上述理论,对一个上层为泥岩,下层为砂岩的二层介质模型,探讨孔隙纵横比对界面反射系数和振幅的影响。泥岩的纵波速度vp为3 400 m/s,横波速度vs为1 450 m/s,泊松比σ为0.38,密度ρ为2.5 g/cm3。下层的砂岩由固体骨架和孔隙水组成,孔隙度为0.15,骨架的主要矿物成分是石英,含少量长石和方解石,参数见表1。计算过程中,保持骨架和孔隙流体及孔隙度不变,只改变孔隙纵横比,根据Berryman散射理论模型计算复合介质的等效弹性参数和速度(图1a),图1b是饱水砂岩孔隙纵横比与速度的关系,这是根据Gassmann方程做流体替换得到的饱气砂岩孔隙纵横比与