疏松砂岩干岩石模量研究
魏冰涛, 杨国权, 曹文俊
中国石油大学(华东),山东 青岛 266580

作者简介: 魏冰涛(1990-),男,硕士在读研究生,主要从事岩石物理研究工作。E-mail:13156050650@163.com

摘要

目前对疏松砂岩的研究主要是基于一些模型的思想,这些模型都是以Hertz-Mindlin和Walton的接触理论作为基础,从而很好的预测疏松砂岩的干岩石模量。但是,这传统的疏松砂岩模型并不能很好地表征疏松砂岩的弹性性质,经过研究发现,HMR、HMS预测的干岩石的体积模量比实际测量值略低,如果用Hashin-Shtrikman下限方程内插,会产生更大的误差;HMR、WR预测的剪切模量都会比实际值高出很多,主要是没有考虑颗粒边界的相对运动;WS的适用性不是很好,有些地方预测效果很好,有些地区效果很差;HMS的预测效果很好,适用性很强,但是预测的结果都比实际值略低,如果用Hashin-Shtrikman上限进行内插,会取得很好的结果,并将该结果与传统的疏松砂岩模型对比分析,发现该方法具有很高的精度。

关键词: 疏松砂岩模型; Hertz-Mindlin; Walton; Hashin-Shtrikman; HMS; HMR; WR; WS
中图分类号:P631 文献标志码:A 文章编号:1000-8918(2016)05-0986-05 doi: 10.11720/wtyht.2016.5.23
Research on the dry modulus of the unconsolidated sandstone
WEI Bing-Tao, YANG Guo-Quan, CAO Wen-Jun
China University of Petroleum (East China),Qingdao 266580,China
Abstract

At present,the study of loose sand mainly focuses on some commonly used models which are based on the Hertz-Mindlin and Walton contact theories.In this paper,the authors tested the capability of various rock physics models to predict dry modulus of the sand samples by comparing the estimates with results measured in laboratory.The study shows that the traditional model of unconsolidated sand is not good at characterizing its elastic properties.Using either the HMR or HMS theories to define the high porosity endpoint gives a lower-prediction of the dry bulk modulus which will have a much bigger error when combined with HS lower bounding equation.Both HMR and WR result in an over-prediction of the dry shear modulus because they do not consider the deformation across the contact.As WS cannot give a good fit to all the areas,it is rarely used.Despite the fact that the results are somewhat lower than the real data,the dry modulus is best described by HMS.With the help of the HS higher bounding equation,the predicted dry modulus agrees fairly well with the measured data.So the new constructed method has better application effect and deserves much attention.

Keyword: unconsolidated sand model; Hertz-Mindlin; Walton; Hashin-Shtrikman; HMS; HMR; WR; WS

干岩石模量的求取, 一直以来都是一个未完全解决的问题[1]。干岩石模量本质上反映的是在孔隙压力不变的条件下围压变化所引起的岩石体积改变[2]。首先可以通过实验室直接测定干岩石的纵横波和密度, 进而计算干岩石模量, 其次借助经验公式以及等效介质模型[3]。对于地下深层胶结岩石, 可以通过实验室测量干岩石模量, 但是对于浅层的疏松砂岩, 由于没有足够多的岩芯资料, 因此对该方面的研究比较少, 所以不能很好地研究其干岩石模量。

Dvormkin和Nur引入高孔隙度砂岩理论模型, 结合Hertz-Mindlin理论和Hashin-Shtriman理论, 提出了一种计算疏松砂岩干岩石模量的有效模型, 即疏松砂岩模型。该模型对于富含石英的砂岩具有很好的应用效果, 但是黏土含量相对较大时(> 20%), 应用效果不是很好[4, 5, 6, 7], 同时, Bachrach、Dvormkin和Nur指出, Hertz-Mindlin理论没有考虑岩石颗粒间的相对运动, 比如滑动、旋转等, 所以预测的剪切模量高于实际测量的剪切模量[8]

当等效介质在流体静压力下, 介质是统计上的各向同性时, Walton考虑了球体是无限粗糙和无限光滑时的两种情况, 推导出了对应的弹性模量公式。但是地壳压力既不是静水压力也不是单轴压力, 所以在定性评价油气储层中的非固结砂岩的等效模量时, 该模型应谨慎使用[10]。Zimmer发现, Hertz-Mindlin理论中岩石对压力的敏感性比实际值小, 当假设岩石颗粒间的剪切应力很大时, 用该理论估算的干岩石速度远远大于实际测量值[11]。Andersen和Johansen通过实验发现, 当岩石颗粒间的剪切应力很大时, 前文提到的两种理论都会过高的估算岩石的速度, 纵波速度比实际值高1.5%, 横波速度比实际值高1.8%, 这种差异会随着孔隙度或者含水饱和度的增加而减小, 当假设岩石颗粒间的剪切应力很小时, 岩石颗粒会发生相对运动, 在该理论(Hertz-Mindlin-Smooth或Walton-Smooth)下预测的速度值更加准确[12]

上面提到的几种模型可以用来计算干岩石模量, 同时结合实际岩芯数据, 来测试模型的准确度。实际结果发现Hertz-Mindlin-Smooth的效果比较好, 但是预测值仍然比实际弹性模量小, 可以采用类似于疏松砂岩模型的方法[9], 结合Hashin-Shtriman上边界理论估算岩样的弹性模量。经过实验发现, 该估算值与实际值具有很好的一致性, 具有很好的研究意义。

1 方法原理
1.1 Walton模型

Walton假设球体和等效介质都是弹性均匀的, 球体任意充填且统计上是各向同性的, 当压力为0时, 邻近球体之间的接触是点接触, 从而推导出任意充填的弹性球体的等效弹性模量的一组方程[10]

当球体是无限粗糙时(Walton-Rough):

KWR=163(1-φ)2C2Peπ4B2, μWR=35KWR5B+A2B+A, A=14π1μs-1μs-λs, B=14π1μs+1μs-λs, λs=Ks-23μs,

其中:KWRμ WR分别是球体无限粗糙时的等效介质的体积模量和剪切模量, Pe是等效压力, φ 是孔隙度, λ sμ s是基质的拉梅常数和剪切模量, C是配位数, 代表岩石颗粒间的接触关系。

当球体是无限光滑时(Walton-Smooth):

KWS=KWR,   μWS=35KWR,

其中:KWSμ WS分别是球体无限光滑时的等效介质的体积模量和剪切模量。

1.2 Hertz-Mindlin模型

基于Hertz-Mindlin模型的等效模量理论常被用于定量描述疏松砂岩在一定储层条件下的弹性模量的特征, H-M模型假设组成粒状介质的颗粒接触边界是紧密连接的, 在应力的作用下不产生相对滑动, 该理论模型通常计算出的体积模量与实际测量较为吻合, 但剪切模量往往高于实际测量值。造成上述差异的主要原因是H-M模型是基于均匀应变近似的, 即在一个无限小的应力扰动下, 颗粒边界的力学特征与颗粒集合体的宏观力学特征一致, 符合规则排列颗粒的应力应变特征, 但是对于不规则排列的颗粒集合体会在局部产生不均匀应变, 从而在颗粒边界造成扭转和旋转运动[13]。类似于Walton模型一样, 该理论也被发展成两种颗粒接触理论: HMR(H-M-Rough)和HMS(H-M-Smooth)。

HMR理论:

KHMR=n2(1-φ)2μS218π2(1-σ)2Pe13, μHMR=5-4σ5(2-σ)3n2(1-φ)2μS22π2(1-σ)2Pe13,

其中:KHMRμ HMR分别是等效介质的体积模量和剪切模量; Pe是有效压力(即围压和空隙压力之差); μ σ 是固体相的剪切模量和泊松比; n是配位数, 代表岩石颗粒间的接触关系。

HMS理论:

KHMS=KHMR,   μHMS=35KHMR,

其中:KHMSμ HMS分别是等效介质的体积模量和剪切模量。

KHMμ HM分别是临界孔隙度φ c(即沉积孔隙度)时干岩石的体积模量和剪切模量; Pe是有效压力(即围压和空隙压力之差); μ σ 是固体相的剪切模量和泊松比; n是配位数, 代表岩石颗粒间的接触关系。

1.3 疏松砂岩模型

Dvorkin和Nur引入了两个高孔隙度砂岩理论模型, 疏松砂岩模型或“ 未固结线” , 来描述当分选变差时, 速度— 孔隙度关系如何变化。“ 分选很好” 的端元用分选很好的相似颗粒填充表示, 这些颗粒的弹性性质由颗粒接触点的弹性确定。“ 分选很好” 的端元具有大约40%左右典型的临界孔隙度。疏松砂岩模型把分选不好的砂岩表示为修改的“ 分选很好” 的端元, 修改时通过附加沉积在空隙空间的小颗粒完成的, 这些附加的颗粒使分选变差, 降低了孔隙度, 并且轻微的增加了岩石硬度[14]

在临界孔隙度, 干的分选很好的端元的弹性模量被模拟成一个易于受围压影响的弹性球充填。这些模量由HMR理论给出, 疏松砂岩模型的另一个端元在零孔隙度, 具有矿物体模量Ks和剪切模量μ s, 孔隙度在0和临界孔隙度之间分选不好的砂岩模量用Hashin-Shtrikmanxia下边界在矿物点和分选很好的端元之间内插。这里把分选不好的砂岩想象为一些由细小颗粒砂子软“ 壳” 包裹的大颗粒, 这是Hashin-Shtrikman界限的体现。

Kdry=φ/φcKHMR+4μHMR/3+1-φ/φcKS+4μHMR/3-1-4μHMR/3, μdry=φ/φcμHMR+Z+1-φ/φcμS+Z-1-Z, Z=μHMR69KHMR+8μHMRKHMR+2μHMR,

其中:φ c是临界孔隙度, 其他参数意义同上。

疏松砂岩模型代表了在砂岩单元内速度— 孔隙度— 分选变化。对于富食石英的砂岩, 分选变化是由小的石英颗粒填充进大石英颗粒之间的孔隙空间中。然而, 分选变差通常与黏土含量增加有关, 如果黏土含量相对比较大(> 20%), 该模型效果不是很好。

1.4 改进的疏松砂岩模型

上述提到的疏松砂岩模型, 是HMR模型与HS模型的结合, 但是基于HMR模型基本假设的缺陷, 剪切模量往往高于实际测量值。经过实验发现, HMS模型预测的弹性模量与实际值较为接近, 比实际值略低, 类似于构建疏松砂岩模型的方法, 在孔隙度在0和临界孔隙度之间分选不好的砂岩模量用Hashin-Shtrikmanxia上边界内插[18]:

Kdry=φ/φcKHMS+4μS/3+1-φ/φcKS+4μS/3-1-4μS/3, μdry=φ/φcμHMS+Z+1-φ/φcμS+Z-1-Z, Z=μS69KS+8μSKS+2μS

2 实际数据测试

应用Galveston Beach, Gulf of Mexico, Pomponio Beach, Merritt 4个地区的岩样数据, 对上述方法进行测试。图1是用4种不同的方法预测的4个地区的疏松砂岩的体积模量, 4种方法分别是HMR(HM-Rough), WR(W-Rough), HMS(HM-Smooth), WS(W-Smooth), 图1a、b、c、d中的蓝点是预测的岩样的体积模量, 即图中纵坐标、横坐标是实测的岩样的体积模量, 理论上预测值和实测值应该相等, 对应图中红线(y=x), 但是, 从预测的结果可以看出, WR和WS的适用性较差, 在Gulf of Mexico和Pomponio Beach地区的效果较好, 但是在Galveston Beach 和Merritt地区的应用效果较差; HMR和HMS的结果较好, 与实际值相差相差较小, 但是都小于实际值。

图1 预测的干岩石体积模量对比
a— Galveston Beach sands; b— Gulf of Mexico sands; c— Pomponio Beach sands; d— Merritt sands

图2是用4种不同的方法预测的4个地区的疏松砂岩的剪切模量, 4种方法分别是HMR(HM-Rough), WR(W-Rough), HMS(HM-Smooth), WS(W-Smooth), 从预测的结果可以看出, HMR和WR的应用效果很差, 预测值与实际值相差很大; WS在Galveston Beach地区预测的疏松砂岩的的剪切模量比实际值大很多, 但是在其他地区效果较好; 相比之下, HMS在4个地区取得了很好的效果。

图2 预测的干岩石剪切模量对比
a— Galveston Beach sands; b— Gulf of Mexico sands; c— Pomponio Beach sands; d— Merritt sands

图3是用两种不同的方法预测的4个地区的疏松砂岩的弹性模量, 两种方法分别是HMS-HS+(HMS-HS上限), HMR-HS-(HMR-HS下限), 从预测的结果可以看出, HMR-HS-应用效果很差, 预测值都很小, 预测值与实际值相差很大; HMS-HS+的预测结果很好, 具有很好的研究意义。

图3 预测的干岩石弹性模量对比
a— Galveston Beach sands; b— Gulf of Mexico sands; c— Pomponio Beach sands; d— Merritt sands

3 结论和认识

笔者主要是针对一些常用疏松砂岩模型展开研究, 并通过实际岩芯数据测试模型的准确性和适用性。研究发现, HMR和WR模型预测的剪切模量比实际值大很多, 主要是没有考虑岩石颗粒间的相对运动, WS模型较好, 但是适用性较差, HMS的模型的应用效果较好, 具有很好的研究意义。

同时发现HMR模型预测的疏松砂岩的体积模量比实际值小, 如果采用传统的疏松砂岩模型, 即用HS下限进行内插, 会使预测结果产生更大的误差。HMS预测的弹性模量比实际值都小, 但是和实际值具有很好的一致性, 采用HS上限进行内插后发现, 预测的弹性模量与实际值的的误差很小, 证明了该方法的有效性。

The authors have declared that no competing interests exist.

参考文献
[1] 张佳佳, 李宏兵, 刘怀山, . 几种岩石骨架模型的适用性研究[J]. 地球物理学进展, 2010, 25(5): 1697-1702. [本文引用:1]
[2] 李宏兵, 张佳佳, 姚逢昌. 岩石的等效孔隙纵横比反演及其应用[J]. 地球物理学报, 2013, 56(2): 608-615. [本文引用:1]
[3] 伍向阳, 陈祖安. 干岩石体积模量估计[C]//中国科学院地质与地球物理研究所2001学术论文汇编, 2001. [本文引用:1]
[4] 刘浩杰. 地震岩石物理研究综述[J]. 油气地球物理, 2009, 7(3): 1-8. [本文引用:1]
[5] Avseth P, Bachrach R, Fristad T, et al. Application of contact mechanics to hydrocarbon prediction in shallow sediments[C]//The 69th EAGE meeting, 2007. [本文引用:1]
[6] Avseth P, Mukerji T, Mavko G. Quantitative Seismic Interpretation: Applying Rock Physics Tools to Reduce Interpretation risk[M]. London: Cambridge University Press, 2005. [本文引用:1]
[7] Bachrach R, Avseth P. Rock physics modeling of unconsolidated sand s: Accounting for nonuniform contacts and heterogeneous stress fields in the effective media approximation with applications to hydrocarbon exploration[J]. Geophysics, 2008, 73(6): 197-209. [本文引用:1]
[8] Bachrach R, Dvorkin J, Nur A. Seismic velocities and Poisson's ratio of shallow unconsolidated sand s[J]. Geophysics, 2000, 65(2): 559-564. [本文引用:1]
[9] Dvorkin J, Nur A. Elasticity of high porosity sand stones: Theory for two North Sea datasets[J]. Geophysics, 1996, 61(5): 1363-1370. [本文引用:1]
[10] Hill R. Elastic properties of reinforced solids: Some theoretical principles[J]. Journal of the Mechanics and Physics of Solids, 1963, 11(5): 357-372. [本文引用:2]
[11] Mavko G, Mukerji T, Dvorkin J. The Rock Physics Hand book[M]. London: Cambridge University Press, 1998. [本文引用:1]
[12] Mindlin R D. Compliance of elastic bodies in contact[J]. Journal of Applied Mechanics, 1949, 16: 259-268. [本文引用:1]
[13] Walton K. The effective elastic moduli of a rand om packing of spheres[J]. Journal of the Mechanics and Physics of Solids, 1987, 35(2): 213-226. [本文引用:1]
[14] Zimmer M, Prasad M, Mavko G, et al. Seismic velocities of unconsolidated sand s: Part 1-Pressure trends from 0. 1 to 20 MPa[J]. Geophysics, 2007, 72(1): 1-13. [本文引用:1]
[15] Zimmer M, Prasad M, Mavko G, et al. Seismic velocities of unconsolidated sand s: Part 2-influence of sorting-and compaction-induced porosity variation[J]. Geophysics, 2007, 72(1): 15-25. [本文引用:1]
[16] Digby P J. The effective elastic moduli of porous granular rocks[J]. Journal of Applied Mechanics, 1981, 48(4): 803-808. [本文引用:1]
[17] Domenico S N. Elastic properties of unconsolidated porous sand reservoirs[J]. Geophysics, 1977, 42(5): 1339-1368. [本文引用:1]
[18] Hashin Z, Shtrikman S. A variational approach to the elastic behavior of multiphase materials[J]. Mech. Phys. Solids. , 1963, 11(2): 127-140. [本文引用:1]