[发明专利]确定碳酸盐岩气水相对渗透率的方法及装置

摘要

本发明公开了一种确定碳酸盐岩气水相对渗透率的方法及装置，该方法包括确定碳酸盐岩中孔隙的分布；确定各孔隙的气相及水相流速；确定各孔隙的气相及水相流量；根据各孔隙的气相及水相流速、气相及水相流量，确定碳酸盐岩气水相对渗透率。本发明无需人工参与，也无需采用水或者汞等的注入过程，因此能够解决碳酸盐岩气水相对渗透率确定时出现的时间长、误差大的不足。

主权项

一种确定碳酸盐岩气水相对渗透率的方法，其特征在于，该方法包括：确定碳酸盐岩中孔隙的分布；确定各孔隙的气相及水相流速；确定各孔隙的气相及水相流量；根据各孔隙的气相及水相流速、气相及水相流量，确定碳酸盐岩气水相对渗透率；所述碳酸盐岩中孔隙的分布根据多孔介质分形几何原理确定，包括：取特征长度为L0的正方形岩样，则毛细管束直线长度为L0，单元截面面积为A，单元截面上孔隙半径在r和r+dr之间的孔隙个数为：N=(rmax/rmin)Df;]]>孔隙大小分布的概率密度函数为：f(r)=DfrminDfr-(Df+1);]]>多孔介质毛细管实际弯曲长度为：L=(2r)1-DTL0DT;]]>所述确定各孔隙的气相及水相流速，包括：根据以下方式确定各孔隙的气相及水相流速：vw=-(p1-p2)r24μwL+(p1-p2)(r0-δ)24μwL;]]>vg=-(p1-p2)r24μgL+(p1-p2)r124μgL+(p1-p2)[(r0-δ)2-r12]4μwL;]]>其中：Df为孔隙分形维数；DT为迂曲度分形维数；下标：max表示最大；min表示最小；r为毛细管中心到任意位置的径向距离，μm；dr为径向距离增量，μm；r1为气、水交界面相半径，μm；r0为毛细管半径，μm；μw为水相黏度，mPa·s；μg为气相黏度，mPa·s；p1为入口端压力，MPa；p2为出口端压力，MPa；δ为束缚水膜厚度，μm；下标：w为水相；g为气相；所述确定各孔隙的气相及水相流量，包括：根据达西公式按以下方式确定各孔隙的气相及水相流量为：Qw=AKKrw(p1-p2)μwL0;]]>Qg=AKKrg(p1-p2)μgL0;]]>其中，Qw为水相总流量，cm3/s；Qg为气相总流量，cm3/s；A为单元截面面积，cm2，K为岩石绝对渗透率，μm2；Krw为水的相对渗透率，百分数；Krg为气的相对渗透率，百分数。