一种建立饱和重油岩石波速随温度变化的理论模型的方法

1.本发明属于非常规油气储层岩性定量预测领域，具体涉及一种建立饱和重油岩石波速随温度变化的理论模型的方法。


背景技术：

2.储层岩石可以被认为是非均质多孔粘弹性介质，流体完全或部分饱和，且流体在深度、温度和压力不同条件作用下，表现出不同的性质。因此，研究饱和流体对这些性质和流体-骨架相互作用的影响是很重要的。
3.然而，由于重油的高黏度，很难对重油的物理性质进行建模。许多作者通过实验室测量手段研究了重油饱和岩石的弹性特性(nur等.1984；eastwood 1993；schmitt 1999；yuan等.2017)。重油在室温和压力下是固体，但加热后会融化，从而降低了剪切波速度(behura等.2007)。han(2007)等人的研究表明，gassmann-biot理论得到的预测值与实测实验数据不符(gassmann 1951；biot 1956)。
4.由于重油黏度高，常规的流体替换方法在这里不适用。一些学者采用cpa方法对该问题进行了研究，该方法考虑了温度对流体黏度的影响，并将流体视为粘弹性介质(berryman 1980；gurevich等.2008；makarynska等.2010)。然而，所提出的模型仍存在不确定性，限制了实际应用。
5.cpa方法允许考虑由maxwell粘弹性模型描述的高黏度孔隙流体。carcione(2020)等人已经成功地使用这个模型来描述地球上的熔融物质，但并未涉及温度和压力对微裂隙的影响。


技术实现要素：

6.针对现有技术的不足，本发明提供一种建立饱和重油岩石波速随温度变化的理论模型的方法，基于cpa方法以及david和zimmerman(2012)提出的微观孔隙结构理论，建立了一个温度依赖的岩石物理模型。在这里考虑了压力和温度对微裂隙的影响，然后，将随温度变化的波速与实验数据进行了比较，这是之前的方法所没有的。本发明的目的在于替代传统重油岩石温度和波速关系的模型，同时这种方法也可以推广到其他岩性。
7.本发明是通过以下技术方案实现的：
8.一种建立饱和重油岩石波速随温度变化的理论模型的方法，包括以下步骤：
9.步骤1)获取样本的超声实验数据；
10.步骤2)建立具有温度和频率叠加效应的流体模型，计算流体的热学性质；
11.步骤3)计算裂隙相关参数进而计算裂隙孔隙度，建立裂隙孔隙度随温度的变化关系；
12.步骤4)基于重油性质，应用单孔cpa模型；
13.步骤5)结合单孔cpa模型，构建与温度有依赖性的双重孔隙cpa模型；
14.步骤6)利用迭代方法，计算饱和岩石的等效体积模量和剪切模量，进而建立温度
与速度的关系，分析温度对样本波速的影响。
15.优选地，其特征在于，步骤1)中所述样本为油砂，主要矿物成分包括石英和粘土，孔隙流体为重油，api密度为6.6
°
，孔隙度为40.96％，颗粒密度为2490kg/m3，在频率为1mhz的4个温度点进行超声实验，在实验过程中，将孔隙压力设置为0psi，获得重油砂的岩石物理参数，并利用voigt-reuss-hill平均值计算得到样品矿物的体积模量和剪切模量。
16.优选地，所述步骤2)中采用maxwell模型得到孔隙流体剪切模量sf随温度和频率变化的函数，和batzle-wang经验公式得出流体的热物理性质，建立流体模型，如公式(1)：
[0017][0018]
公式(1)中：s
∞
为高频下流体的剪切模量，τ为弛豫时间，ω为角频率，η为随温度变化的流体黏度，
[0019]
所述sf还可以表示成如公式(2)：
[0020]
sf＝re(sf) iim(sf)
[0021]
(2)；
[0022]
公式(2)中：re(sf)和iim(sf)分别为存储模量和损失模量；
[0023]
所述流体剪切模量是指流体的储存模量。
[0024]
优选地，所述步骤3)中根据dz模型计算裂隙孔隙度φ2(p)，如公式(3)：
[0025][0026]
公式(3)中：β为裂隙密度，α
p
为孔隙纵横比；
[0027]
基于在不同温度下随有效压力变化的超声实验数据，利用微分等效介质理论计算相同温度下裂隙密度随有效压力的变化，建立孔隙纵横比α
p
与有效压力p的定量关系，如公式(4)：
[0028][0029]
公式(4)中：es为高有效压力下的杨氏模量，定义为es＝3ts[1-2vs]，其中vs为泊松比。
[0030]
优选地，所述步骤4)中应用单孔cpa模型，计算方程如公式(5)和(6)：
[0031]
φ(t
f-t)pf (1-φ)(t
s-t)ps＝0
[0032]
(5)；
[0033]
φ(s
f-s)qf (1-φ)(s
s-s)qs＝0
[0034]
(6)；
[0035]
公式(5)和(6)中：φ为孔隙度，t和s分别为岩石的体积模量和剪切模量，tf和sf分别为孔隙流体的体积模量和剪切模量，ts和ss分别为岩石基质的体积模量和剪切模量，pf和qf均为孔隙流体的形状因子；ps和qs均为岩石基质形状因子。
[0036]
优选地，所述步骤5)中基于岩石样本，利用voigt-reuss-hill平均值计算得到样品矿物的体积模量和剪切模量，并根据batzle-wang经验公式得出流体的热物理性质，利用maxwell模型得到流体的复剪切模量，然后基于dz模型计算求得裂隙孔隙度，最后将求得的参数全部带入步骤4)提到的单孔cpa模型中，得到温度相关的双重孔隙cpa模型；
[0037]
具体为：将公式(1)至(4)代入公式(5)和(6)，得到温度依赖的双重孔隙cpa模型，如公式(7)和(8)：
[0038]
φ1(t
f-t)p
f1
φ2(t
f-t)p
f2
(1-φ)(t
s-t)ps＝0
[0039]
(7)；
[0040]
φ1(s
f-s)q
f1
φ2(s
f-s)q
f2
(1-φ)(s
s-s)qs＝0
[0041]
(8)；
[0042]
公式(7)和(8)中：φ为孔隙度，φ1和φ2分别为孔1和孔2的孔隙度，且φ＝φ1 φ2，p
f1
、q
f1
和p
f2
、q
f2
分别为孔1和孔2孔隙流体的形状因子。
[0043]
优选地，所述步骤6)中结合温度相关的双重cpa模型，采用迭代法计算饱和岩石的等效体积模量和剪切模量，求出相应饱和岩石的纵、横波速度，从而分析温度对油砂波速的影响；
[0044]
具体为：结合公式(7)和(8)，采用迭代法计算饱和岩石的n 1次迭代时的体积模量t
n 1
和n 1次迭代时的剪切模量s
n 1
，如公式(9)和(10)：
[0045][0046][0047]
公式(9)和(10)中：和分别为n次迭代时孔1和孔2孔隙流体的形状因子，分别为n次迭代时岩石基质的形状因子；
[0048]
根据公式(9)和(10)的结果，求出饱和岩石的纵波速度v
p
、横波速度vs和总密度ρ，如公式(11)、(12)和(13)：
[0049][0050][0051]
ρ＝(1-φ)ρs φρf[0052]
(13)；
[0053]
公式(13)中：ρs为岩石基质密度，ρf为流体密度。
[0054]
优选地，步骤1)所述样本为油砂或碳酸盐岩。
[0055]
本发明的有益效果如下：
[0056]
本发明考虑了压力和温度对微裂隙的影响，基于相干势近似法(cpa)，结合温度和频率相关的孔隙流体经验方程和david-zimmerman(dz)模型，提出了一个温度依赖的双重孔隙模型，并采用maxwell模型得到流体复剪切模量随温度和频率变化的函数，来描述饱和重油岩石波速随温度变化的情况。为了验证模型的有效性，对油砂样品进行测试，通过预测的波速与实测的进行对比，表明该模型能定量地描述波速随温度的变化规律。并且该方法可以推广到其他岩性，如碳酸盐岩。
附图说明
[0057]
图1为一种建立饱和重油岩石波速随温度变化的理论模型示意图；
[0058]
图2为饱和水砂岩和饱和重油砂岩(api＝-5
°
)的体积模量(a)和剪切模量(b)与频率的关系图；
[0059]
图3为重油(api＝6.6
°
)黏度随温度的变化曲线图；
[0060]
图4为重油剪切模量(api＝6.6
°
)随温度和不同频率变化的关系图；
[0061]
图5为1100psi有效压力下，重油砂的p波和s波速度随温度变化的关系图。
具体实施方式
[0062]
以下结合附图附表及具体实施例对本发明作进一步的详细说明。需要说明的是，以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
[0063]
实施例1
[0064]
一种建立饱和重油岩石波速随温度变化的理论模型的方法，如图1所示，具体步骤如下：
[0065]
1、获取样本的超声实验数据。
[0066]
采用li等.(2016)测量的油砂样品，样品主要矿物成分包括石英和粘土，孔隙流体为重油，api密度为6.6
°
，孔隙度为40.96％，颗粒密度为2490kg/m3，在频率为1mhz的4个温度点进行超声实验，在实验过程中，将孔隙压力设置为0psi，获得重油砂的岩石物理参数，并利用voigt-reuss-hill(v-r-h)平均值计算得到样品矿物体积模量和剪切模量。
[0067]
2、建立具有温度和频率叠加效应的流体模型，计算流体的热学性质。
[0068]
采用maxwell模型得到孔隙流体剪切模量sf随温度和频率变化的函数，和batzle-wang经验公式得出流体的热物理性质，建立流体模型，如公式(1)：
[0069][0070]
公式(1)中：s
∞
为高频下流体的剪切模量，τ为弛豫时间，ω为角频率，η为随温度变化的流体黏度，
[0071]
所述sf还可以表示成如公式(2)：
[0072]
sf＝re(sf) iim(sf)
[0073]
(2)
[0074]
公式(2)中：re(sf)和iim(sf)分别为存储模量和损失模量；
[0075]
本实施例所考虑的流体剪切模量是指流体的储存模量。
[0076]
3、计算裂隙相关参数进而计算裂隙孔隙度，建立裂隙孔隙度随温度的变化关系。
[0077]
根据dz模型计算裂隙孔隙度φ2(p)，如公式(3)：
[0078][0079]
公式(3)中：β为裂隙密度，α
p
为孔隙纵横比。
[0080]
基于在不同温度下随有效压力变化的超声实验数据，利用微分等效介质理论计算相同温度下裂隙密度随有效压力的变化，建立孔隙纵横比α
p
与有效压力p的定量关系，如公式(4)：
[0081][0082]
公式(4)中：es为高有效压力下的杨氏模量，定义为es＝3ts[1-2vs]，其中vs为泊松比。
[0083]
4、基于重油性质，应用单孔cpa模型。
[0084]
基于传统的建模方法不适用饱和重油岩石，gassmann-biot理论的假设对其无效，为了分析高黏度流体的影响，本实施例基于单孔cpa模型进行实现。计算方程如公式(5)和(6)：
[0085]
φ(t
f-t)pf (1-φ)(t
s-t)ps＝0
[0086]
(5)
[0087]
φ(s
f-s)qf (1-φ)(s
s-s)qs＝0
[0088]
(6)
[0089]
公式(5)和(6)中：φ为孔隙度，t和s分别为岩石的体积模量和剪切模量，tf和sf分别为孔隙流体的体积模量和剪切模量，ts和ss分别为岩石基质的体积模量和剪切模量，pf和qf均为孔隙流体的形状因子；ps和qs均为岩石基质形状因子。
[0090]
5、结合单孔cpa模型，构建与温度有依赖性的双重孔隙cpa模型。
[0091]
基于岩石样本，利用v-r-h平均值计算得到样品矿物的体积模量和剪切模量，并根据batzle-wang经验公式得出流体的热物理性质，利用maxwell模型得到流体的复剪切模量，然后基于dz模型计算求得裂隙孔隙度，最后将求得的参数全部带入步骤4提到的单孔cpa模型中，得到温度相关的双重孔隙cpa模型。
[0092]
将公式(1)至(4)代入公式(5)和(6)，得到温度依赖的双重孔隙cpa模型，如公式(7)和(8)：
[0093]
φ1(t
f-t)p
f1
φ2(t
f-t)p
f2
(1-φ)(t
s-t)ps＝0
[0094]
(7)
[0095]
φ1(s
f-s)q
f1
φ2(s
f-s)q
f2
(1-φ)(s
s-s)qs＝0
[0096]
(8)
[0097]
公式(7)和(8)中：φ为孔隙度，φ1和φ2分别为孔1和孔2的孔隙度，且φ＝φ1 φ2，p
f1
、q
f1
和p
f2
、q
f2
分别为孔1和孔2孔隙流体的形状因子。
[0098]
6、利用迭代方法，计算饱和岩石的等效体积模量和剪切模量，进而建立温度与速度的关系，分析温度对样本波速的影响。
[0099]
结合温度相关的双重cpa模型，采用迭代法计算饱和岩石的等效体积模量和剪切模量，求出相应饱和岩石的纵、横波速度，从而分析温度对油砂波速的影响。
[0100]
结合公式(7)和(8)，采用迭代法计算饱和岩石的n 1次迭代时的体积模量t
n 1
和n 1次迭代时的剪切模量s
n 1
，如公式(9)和(10)：
[0101][0102][0103]
公式(9)和(10)中：和分别为n次迭代时孔1和孔2孔隙流体的形状因子，分别为n次迭代时岩石基质的形状因子。
[0104]
根据公式(9)和(10)的结果，求出饱和岩石的纵波速度v
p
、横波速度vs和总密度ρ，如公式(11)、(12)和(13)：
[0105][0106][0107]
ρ＝(1-φ)ρs φρf[0108]
(13)；
[0109]
公式(13)中：ρs为岩石基质密度，ρf为流体密度。
[0110]
7、分析温度对油砂波速的影响。
[0111]
以下为本实施例的一个实际应用，说明基于一种建立饱和重油岩石波速随温度变化的理论模型的实现过程。
[0112]
本实施例先从流体性质角度出发，考虑了低黏度和高黏度流体体积模量和剪切模量随孔隙度及温度的变化情况，即水和重油(api＝-5
°
)，假设了硬孔隙的纵横比为0.5。图2
显示了体积模量(a)和剪切模量(b)随频率和温度变化的关系，从图2明显可以看到，当温度从40℃升高到60℃时，模量频散急剧降低。
[0113]
紧接着，本实施例给出了温度对重油性质的影响，如图3、图4所示。图3显示了黏度与温度的关系，可以看到黏度随温度的升高而降低。图4为不同频率下重油剪切模量与温度的关系，由式(1)和(2)可知。在这种情况下，高频下重油的剪切模量(s
∞
)可通过对超声波剪切模量的线性回归拟合得到(li等.2016)。在超声频率下，剪切模量在48.7℃左右为零。结果表明，随着频率的增加，流体点向右侧移动。在1mhz时，重油流体点黏度约为1000cp。
[0114]
图5所示为基于温度相关的双重孔隙模型得到的预测结果，与油砂超声实验测得的实验数据基本一致，说明该模型能够有效地描述了油砂波速随温度的变化规律。当温度小于流体点温度时，重油处于黏度迅速降低的准固相，速度随温度升高而急剧下降。当温度超过流体点时，重油为流体，速度呈线性下降。
[0115]
综上所述，本发明是基于相干势近似法(cpa)，结合温度和频率相关的孔隙流体经验方程和david-zimmerman(dz)模型，提出了一个温度依赖的双重孔隙模型，并采用maxwell模型得到流体复剪切模量随温度和频率变化的函数，来描述饱和重油岩石波速随温度变化的情况。结果表明：当温度低于流体点温度时，重油为高黏度准固相非牛顿流体，剪切模量不可忽略，此时波速随温度升高而快速降低；随着温度升高超过流体点温度时，重油接近流体相，此时剪切模量作用可忽略不计，波速呈线性下降趋势。在实施例1中考虑了油砂中超声波速的频率依赖特征，但本发明的方法可以推广到其他岩性，如碳酸盐岩。
[0116]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。