[发明专利]基于纵、横波测井资料的三种孔隙类型定量反演方法

说明书

技术领域

本发明涉及岩石物理及测井评价技术领域，尤其涉及一种基于纵、横波测井资料的三种孔隙类型定量反演方法。

背景技术

孔隙形状会明显地影响着在储层中传播的地震波的速度，关于这点已得到许多理论和实验研究结果的证实（如等，1976；Yan等，2002；Sun等，2004）。搞清楚储层的孔隙结构形状将是非常重要的，这即是孔隙形状的定量反演问题。这个问题的研究尤其对储集空间复杂的储层更为重要，例如我国西部塔里木盆地的以溶蚀孔、洞、缝等各种各样的孔隙类型共存的奥陶系碳酸盐岩储层。面对如此复杂的这类储层，如果能够从井上定量地刻画出每个储层段所发育的各种孔隙类型及其百分比，那将会为储层的评价和描述提供一个全新的参考信息，由于测井信息是硬数据，所以如果能找到一个合理的方法能较准确地定量确定出井点位置的各种孔隙结构类型，这将使我们更好地理解和刻画井眼附近的地震波的传播特征，同时储集空间类型的反演结果也为地震尺度下的缝洞雕刻研究提供一个参考信息。如果能精细刻画出储层的孔隙类型及其所占比例，将会对地区的储层描述给出更加完善的孔隙空间参数。对于这个问题，Anselmetti和Eberli(1999)在讨论碳酸盐岩速度控制因素时提出了一种利用速度偏离（即利用实测速度与利用Wyllie平均时间公式预测的速度的差值）为正时即代表对应储层为以铸模孔为主导孔隙的识别孔隙类型的思想，但这只是定性的，并不能定量地算出铸模孔所占的总孔隙的比例。Kumar和Han(2005)提出了一种在孔隙度和纵波速度已知情况下利用DEM和Gassmann方程估算孔隙表面比的方法，但是横波速度则没有考虑进去；Xu等(2009)将原始的砂泥岩的Xu-White模型推广到碳酸盐岩储层中，并利用该模型根据建立的孔隙度与纵波速度的交汇反演出了“类裂缝”孔隙和刚性孔隙所占的百分比。然后，Sain等(2008)提出了一种将该拓展的Xu-White模型反演的地震孔隙（类裂缝孔隙、刚性孔隙）进一步拆分为所谓的岩石物理孔隙（微孔隙、粒间孔隙、粒内孔隙、洞穴、铸模孔、晶间孔和裂缝）的构想，但是这种构想可行性及结果的可靠性尚有待验证。所述Gassmann方程中文译文为盖斯漫方程，DEM中文译文为微分等效介质。发明内容

本发明要解决的技术问题是提供一种基于纵、横波测井资料的三种孔隙类型定量反演方法，所求结果与井壁成像测井资料及岩心资料吻合性较高，可作为测井评价的一种新的手段。

为了解决上述技术问题，本发明的基于纵、横波测井资料的三种孔隙类型定量反演方法是，将井壁上的复杂的孔隙形状假设为球形孔隙，针形孔隙，硬币状的裂缝孔隙三种，利用测井纵波速度、横波速度，孔隙度，综合使用微分等效介质岩石物理模型和盖斯漫流体替代方程反算地层中球形孔隙、针形孔隙，硬币状的裂缝孔隙三种孔隙类型之间的体积比，并得出这三种孔隙类型的孔隙度。

包括如下步骤：

步骤1，使用实测的纵、横波速度及密度反算出饱和流体岩石的体积模量和剪切模量；然后，利用盖斯漫方程反算出干岩石的体积模量和剪切模量，将其作为反演孔隙类型的输入数据；根据测井密度数据求取孔隙度；

步骤2：根据步骤1所求干岩石的体积模量和剪切模量，利用任意孔隙纵横比的微分等效介质岩石物理模型通过不断迭代计算确定出孔隙纵横比；

步骤3：根据步骤2所计算的孔隙纵横比对孔隙形状进行划分，当孔隙纵横比小于0.02时，将岩石中的孔隙视为全部为硬币状的裂缝孔隙，当孔隙纵横比大于0.12时，将孔隙视为全部为刚性孔隙，如果孔隙在0.02与0.12之间则将孔隙划分为裂缝形孔隙与刚性孔隙的组合，划分其各自占总孔隙的比例；

步骤4：在步骤2划分的基础之上，保持硬币状的裂缝孔隙的体积比不变，利用三维孔隙下的微分等效介质岩石物理模型将步骤2中所得的刚性孔隙进行划分，通过迭代计算求算出球形洞和针形孔所占的比例，约束条件为裂缝形孔隙、针形孔、球形洞的体积大于或者等于0，且三者的和为1；

步骤5：将步骤3和步骤4中所求得的三种孔隙类型，即球形孔隙，针形孔隙，硬币状的裂缝形孔隙的相对体积分别与总孔隙度相乘得到三种孔隙类型的孔隙度。

在步骤3中将步骤2中划分出的孔隙划分为硬币状裂缝形孔隙与刚性孔隙的比值，进而在步骤4中在硬币状裂缝孔隙体积比保持不变的情况下，通过三维孔隙的微分等效介质模型，将刚性孔隙划分为球形孔隙和针形孔隙，通过两次迭代计算得到球形孔隙，针形孔隙，硬币状裂缝形孔隙的体积比。