【技术实现步骤摘要】
一种考虑毛细管力的双相饱和页岩速度频散和预测方法
[0001]本专利技术涉及预测页岩储层速度频散和衰减
,尤其是涉及一种考虑毛细管力的双相饱和页岩速度频散和预测方法。
技术介绍
[0002]White等人(1983)观察到Pierre页岩中声波测井的P波速度测量值始终高于来自VSP的速度(平均6%)。Homby(1995)发现,10KHz测量的页岩速度比1MHz频率下测量的速度低15%。Suarez
‑
Rivera等人(2001)针对Pierre页岩进行了频率依赖性研究,发现在10Hz和500kHz之间存在显著的纵波速度频散,高达45%。Szewczyk等人(2018)研究了地震频带内(1
‑
100Hz)流体饱和度对页岩频散和衰减的影响,认为地震波诱导的流体流动是主要原因,即波致流。已有技术从宏观、中观、微观三个尺度研究了页岩中的波致流现象,分析了流体分布特征、裂缝等诱发的地震波频散衰减和频变特性,并提出了相应的数值理论模型,结果表明多相流体混合和裂缝是导致介质弱非均质性及各向异性,并进一步导致频散和衰减的主要成因。然而,在页岩储层中,复杂的裂隙网络导致毛细管力发育,其对频散和衰减的影响未知。
技术实现思路
[0003]本专利技术的目的是提供一种考虑毛细管力的双相饱和页岩速度频散和预测方法,弥补现有技术方案对毛细管力影响估计不足,利用有限元方法来计算有效VT I介质和频率相关刚度的数值模拟方法,能够正确预测包含毛细管力的页岩储层地震波频散和衰减特征。
...
【技术保护点】
【技术特征摘要】
1.一种考虑毛细管力的双相饱和页岩速度频散和预测方法,其特征在于:包括以下步骤:S1:根据CT扫描结果建立页岩样品模型,获得岩石的物性参数;S2:基于步骤1在切面上进行三角网格划分,在轴向上进行扫描,完成三维空间的网格划分;S3:计算毛细关系和已知参数;S4:通过VTI介质的应力
‑
应变方程,得到复刚度,利用有限元计算出岩石样本的纵横波速度绘制曲线。2.根据权利要求1所述的一种考虑毛细管力的双相饱和页岩速度频散和预测方法,其特征在于:在步骤S3中,毛细关系为:其中,P
w
定义为润湿流体压力的无穷小变化的傅里叶变换,P
n
定义为非润湿流体压力的无穷小变化的傅里叶变换;表示非浸润流体饱和多孔弹性介质的饱和度,表示润湿流体饱和多孔弹性介质的饱和度;表示润湿流体压力,表示非润湿流体压力,表示孔隙率;润湿性流体和非润湿性流体界面接触位置的压力系数为:式中,P'
ca
为毛细管力的导数,C
n
表示非润湿性流体、C
w
表示润湿性流体的柔度;构建毛细管力变化系数,获得包含毛细管力的Biot
‑
willms系数α,计算混合流体的体积模量,计算孔隙空间模量,式中,K
m
表示骨架体积模量,K
s
矿物组分体积模量,K
f
是流体体积模量;计算饱和流体岩石的体积模量,K
u
=K
s
(K
m
+Ξ)/(K
s
+Ξ)其中,K
u
和μ分别表示湿岩体积和干岩剪切模量;流体相的相对位移u
θ
定义为:其中u
s
为固体位移;流体含量ξ
θ
变化定义为变化定义为是哈密尔顿算子;
建立本构关系:τ
ij
(u)=2με
ij
+δ
ij
(λ
u
e
s
‑
B1ξ
n
‑
B2ξ
w
)
ꢀꢀꢀꢀ
(1)(1)e
s
=ε
ii
(u
s
)其中,ε
ij
(u
s
)是固体应变张量的傅里叶变换,e
s...
【专利技术属性】
技术研发人员:孙超,姜志海,李良钰,高崇诚,刘晓阳,旦增措姆,
申请(专利权)人:中国矿业大学,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。