本发明专利技术公开了一种基于Xu‑White模型的三维流体特征的反演方法，属于气田高效开发领域，包括：S1：构建Xu‑White模型，预测纵波速度V<subgt;p</subgt;、横波速度V<subgt;s</subgt;及密度ρ；S2：通过步骤S1得到的纵波速度V<subgt;p</subgt;、横波速度V<subgt;s</subgt;、密度ρ与地震反演得到的纵波速度V<subgt;p</subgt;、横波速度V<subgt;s</subgt;、密度ρ进行构建目标函数；S3：输入压力、含水饱和度，利用最优化算法对所述目标函数进行迭代优化，得到满足精度的含水饱和度S<subgt;w</subgt;和压力P。利用本申请的方法，技术精度较高，计算效率较快，和其他方法相比适用条件更为广泛。
【技术实现步骤摘要】
本专利技术属于气田高效开发领域，具体涉及一种基于xu-white模型的三维流体特征的反演方法。
技术介绍
1、随着不断开发，各气田目前已经相继进入产量递减期，其中以我国的某气田为例，历经十九年生产已进入开发中后期，但依然存在较多的动静储量差异，存在有较大挖潜调整空间。
2、由于其部分气组为低渗储层，空间非均质性强，造成储层连通性差，开发低效，为了实现气田高效开发，需要对气藏进行地震监测，比如气田开发过程中流体饱和度变化，流体压力变化等，进而对已处于开发中后期甚至开发末期的气藏进行剩余气预测并调整挖潜方案。
技术实现思路
1、本专利技术所为了解决
技术介绍
中存在的技术问题，目的在于提供了一种基于xu-white模型的三维流体特征的反演方法。
2、为了解决技术问题，本专利技术的技术方案是：
3、一种基于xu-white模型的三维流体特征的反演方法，所述方法包括：
4、s1：构建xu-white模型，预测纵波速度vp、横波速度vs及密度ρ；
5、s2：通过步骤s1得到的纵波速度vp、横波速度vs、密度ρ与地震反演得到的纵波速度vp、横波速度vs、密度ρ进行构建目标函数；
6、s3：输入压力、含水饱和度，利用最优化算法对所述目标函数进行迭代优化，得到满足精度的含水饱和度sw和压力p。
7、进一步，所述步骤s1包括：
8、利用v-r-h平均得到砂岩与泥岩混合物即基质的弹性模量并作为输入；
>9、使用融合了dem理论的k-t模型来计算干岩石即骨架的弹性模量并作为输入；10、利用加入了压力和含水饱和度参数的流体一起代入gassmann模型估算饱和流体岩石的等效纵、横波速度和密度。
11、进一步，在所述步骤s1中构建基质包括：
12、
13、
14、
15、式中：vsh-表示泥质含量,通过储层建模得到，直接进行输入；
16、k0_s-表示砂岩基质矿物体积模量；
17、k0_n-表示泥岩基质矿物体积模量；
18、μ0_s-表示砂岩基质矿物剪切模量；
19、μ0_n-表示砂泥岩基质矿物剪切模量；
20、ρm_s-表示砂岩基质矿物密度；
21、ρm_n-表示泥岩基质矿物密度；
22、k0-表示基质的体积模量；
23、μ0-表示基质的剪切模量；
24、ρm-表示基质的密度。
25、进一步，在所述步骤s1中在基质中插入孔隙度来构建骨架包括：
26、
27、
28、
29、
30、
31、
32、
33、
34、
35、
36、
37、f6＝1+a[1+g-r(g+γ)]+b(1-γ)(3-4r)
38、
39、
40、f9＝a[g(r-1)-rγ]+bγ(3-4r)
41、a＝-1，b＝0
42、
43、
44、
45、式中：α_s-为砂岩孔隙纵横比；
46、α_n-为泥岩孔隙纵横比；
47、-为岩石孔隙度，通过储层建模得到，直接进行输入；
48、kd-当孔隙度为时，岩石骨架的体积模量；
49、μd-当孔隙度为时，岩石骨架的剪切模量；
50、p和q-代表跟孔隙纵横比α相关的系数；
51、tiijj(α)、fiijj(α)、f1、f2、f3、f4、f5、f6、f7、f8、f9、r、g、γ-都是代表跟孔隙纵横比α相关的函数。
52、进一步，通过孔隙纵横比来描述孔隙形状，将泥质砂岩中的孔隙分成孔隙纵横比大的砂岩孔隙和孔隙纵横比小的泥岩孔隙，使用椭球的短轴与长轴之比(如图所示)定义椭球形的孔隙纵横比：
53、α＝c/a。
54、进一步，在所述步骤s1中，获取流体参数包括：
55、ρf＝swρw+(1-sw)ρg
56、
57、
58、
59、
60、ρw＝ρw0+s×{0.668+0.44s+1×10-6[300p-2400ps+t(80+3t-3300s-13p+47ps)]}
61、其中：
62、kb＝ρwvb2
63、vb＝vw+s(1170-9.6t+0.055t2-8.5×10-5t3+2.6p-0.0029tp-0.0476p2)+s1.5(780-10p+0.16p2)-820s2
64、vw＝x1+p(x2+p(x3+px4))
65、x1＝1402.85+t(4.871+t(-0.04783+t(1.487×10-4-2.197×10-7t)))
66、x2＝1.524+t(-0.0111+t(2.747×10-4+t(-6.503×10-7+7.987×10-10t)))
67、x3＝3.437×10-3+t(1.739×10-4+t(-2.135×10-6+t(-1.455×10-8+5.23×10-11t)))
68、x4＝-1.197×10-5+t(-1.628×10-6+t(1.237×10-8+t(1.327×10-10-4.614×10-13t)))
69、
70、ρw0＝1+1×10-6(-80t-3.3t2+0.00175t3+489p-2tp+0.016t2p-1.3×10-5t3p-0.333p2-0.002tp2)
71、
72、
73、
74、
75、式中t—普通温度，给的标准值80；
76、ta——绝对温度[ta＝t(℃)+273.15]，给的标准值353.15；
77、sw——含水饱和度，自变量；
78、p——压力，自变量；
79、s——为盐浓度，给的标准值0.05；
80、g——气体混合物的比重，即15.6℃大气压下的气体密度与空气密度之比，给的标准值0.1；
81、ρf、ρw、ρg—分别为流体密度、水密度、气密度；
82、ρw0—为盐水密度，温度和压力的函数；
83、kf—为流体体积模量；
84、kw—为水体积模量；
85、kg—为气体积模量；
86、vw-纯水在100℃和大约100mpa时的速度；
87、vb-—为不含气体的盐水的速度；
88、kb—为不含气体的盐水的体积模量；
89、rg—为室温下的气水比；
90、ppr、tpr—分别为拟对比本文档来自技高网...
【技术保护点】
1.一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，所述方法包括：
2.根据权利要求1所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，所述步骤S1包括：
3.根据权利要求2所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，在所述步骤S1中构建基质包括：
4.根据权利要求2所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，在所述步骤S1中在基质中插入孔隙度来构建骨架包括：
5.根据权利要求4所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，通过孔隙纵横比来描述孔隙形状，将泥质砂岩中的孔隙分成孔隙纵横比大的砂岩孔隙和孔隙纵横比小的泥岩孔隙，使用椭球的短轴与长轴之比(如图所示)定义椭球形的孔隙纵横比：
6.根据权利要求2所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，在所述步骤S1中，获取流体参数包括：
7.根据权利要求2所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，在所述步骤S1中，构建含流体的岩石包括：
8.根据权利要求2所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，在所述步骤S1中，输出岩石速度密度包括：
9.根据权利要求1所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，所述目标函数含有两个自变量：含水饱和度Sw和压力P；三个因变量：步骤S1得到的Vp、Vs、ρ的一个函数，即：
10.根据权利要求1所述的一种基于Xu-White模型的三维流体特征的反演方法，其特征在于，利用粒子群算法对所述目标函数进行迭代优化，得到三维流体的压力和饱和度特征。
...
【技术特征摘要】
1.一种基于xu-white模型的三维流体特征的反演方法，其特征在于，所述方法包括：
2.根据权利要求1所述的一种基于xu-white模型的三维流体特征的反演方法，其特征在于，所述步骤s1包括：
3.根据权利要求2所述的一种基于xu-white模型的三维流体特征的反演方法，其特征在于，在所述步骤s1中构建基质包括：
4.根据权利要求2所述的一种基于xu-white模型的三维流体特征的反演方法，其特征在于，在所述步骤s1中在基质中插入孔隙度来构建骨架包括：
5.根据权利要求4所述的一种基于xu-white模型的三维流体特征的反演方法，其特征在于，通过孔隙纵横比来描述孔隙形状，将泥质砂岩中的孔隙分成孔隙纵横比大的砂岩孔隙和孔隙纵横比小的泥岩孔隙，使用椭球的短轴与长轴之比(如图所示)定义椭球形的孔隙纵横比：
6.根据权利要求2...
【专利技术属性】
技术研发人员：王明强，唐静，李浩源，黄旭日，徐云贵，胡叶正，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。