一种考虑粗糙表面特征的岩石毛管力确定方法
【技术领域】
[0001] 本发明涉及微观渗流力学研究领域，尤其涉及一种考虑粗糙表面特征的岩石毛管 力确定方法，适用于确定各种粗糙表面多孔介质中的毛管力大小。
【背景技术】
[0002] 在自然界和工业界中，许多多孔介质表面凹凸不平，明显地表现出粗糙特性，比如 被纳米流体吸附的岩石表面等。粗糙特性影响了多孔介质润湿性、润湿滞后及驱替过程中 由于表面接触造成的的做功损耗。因此多孔介质粗糙表面下的毛细行为与其在光滑表面下 的表现会有较大差别。
[0003] 人们通常从孔隙水平甚至更微观的水平入手研究岩石孔隙间的复杂渗流现象， 进一步不断升级到介观多孔介质及宏观多孔介质，因此孔隙尺度下毛管力的大小是计 算多孔介质毛管力曲线和相对渗透率曲线等宏观特性参数的先决条件。经典的毛管力 公式只适用于圆形截面且未考虑由粗糙表面引起的润湿状态变化。确定不规则截面多 孔介质毛管力的模型最常用的是孔隙网络模型模拟中的MS-P方法（Oren P E，Bakke S,Arntzen 0 J. Extending predictive capabilities to network models[J].SPE Journal, 1998, 3 (04) : 324-336.)，该方法以能量方程为基础，假设驱替前沿的曲率半径与 孔隙截面角隅的界面曲率半径相同，可以考虑不同润湿性及不同截面形状的情况。粗糙表 面润湿可分为Wenzel状态和C-B状态，不同状态下液滴在表面微结构上的分布状态及接触 角不同。建立一种考虑粗糙表面特征的岩石毛管力确定方法将有助于精确计算多孔介质渗 流参数。
【发明内容】
[0004] 本发明的目的是提高粗糙表面多孔介质中毛管力确定的精确性，从而快速、准确 地预测多孔介质渗流参数。为了达到以上目的，本发明提供了一种考虑粗糙表面特征的岩 石毛管力确定方法，该方法充分考虑了粗糙特性对岩石微观结构造成的变化及润湿状态的 改变。
[0005] 本发明技术方案具体步骤如下：
[0006] 1. -种考虑粗糙表面特征的岩石毛管力确定方法，其特征在于以下步骤：
[0007] (1)采用基于CT图像构建孔隙网络模型的方法对岩石孔隙截面进行形状等价，得 到等价形状截面的周长、面积及内角角度。
[0008] (2)利用三维激光表面测量仪扫描岩石表面形貌并成像。
[0009] (3)在三维形貌图中分别截取N个XZ平面和N个YZ平面的岩石轮廓并编号为1， 2, 3…2N，N -般取10,统计凸起高度h、凸起宽度a、凹槽宽度d等粗糙表面平均特性参数。 [0010] (4)利用凸起高度、凸起宽度及凹槽宽度进行粗糙表面等价，同时确定表面平均粗 糙因子rf、接触面中凸起固体所占面积分数f s。
[0011] (5)利用光学接触角测量仪测量非湿相流体液滴在湿相流体环境下光滑岩石表面 的本征接触角，判断液滴在粗糙岩石表面上的润湿状态。
[0012] (6)根据润湿状态建立相应的毛管力方程确定岩石毛管力。
[0013] 其中，所述基于CT图像构建孔隙网络模型的方法为现有方法，具体可参见参考 文献：Hou Jian, Li Zhenquan, Zhang Sunkang, Cao Xulong, Song Xinwang, Gao Debo. Experiment and simulation study on construction of three-dimensional network model. Science in China, Series G, 2008,51 (9):1761-1774. 〇
[0014] 其所述步骤（3)中的粗糙表面特性参数统计方法如下：
[0016] 其中Ral为第i个截面岩石轮廓粗糙度，ζ为第i个截面岩石轮廓最小二乘中线与 波峰交线平均长度，式为第i个截面岩石轮廓最小二乘中线与波谷交线的平均长度。
[0017] 其所述步骤（4)中的等价的粗糙表面特性参数确定公式如下：
[0020] 其所述步骤（6)中的毛管力方程如下：
[0021] 如果液滴润湿为Wenzel状态：
[0023] 如果液滴润湿为C-B状态：
[0025] 其中rd是驱替前沿半径，Arff是截面中非湿相所占面积，L nw是非湿相与湿相的界 面周长，Lns是非湿相与管壁接触周长，γ nw是非湿相与湿相间的界面张力，γ ns是湿相与管 壁间的界面张力，Θ是光滑表面下的本征接触角，0$是粗糙表面下的表观接触角。
[0026] 本发明具有以下有益效果及优点：
[0027] (1)考虑了粗糙表面对多孔介质微观结构的影响，并且适用于不规则截面的孔隙。
[0028] (2)考虑了粗糙表面导致的液滴润湿状态及接触角的变化，可以更加准确地确定 粗糙表面岩石毛管力。
【附图说明】
[0029] 图1为本发明的步骤流程图。
[0030] 图2为粗糙表面三维形貌图。
[0031 ] 图3为截面轮廓示意图。
[0032] 图4为粗糖表面等价不意图。
[0033] 图5为液滴在光滑表面上的接触角测试示意图。
[0034] 图6为粗糖表面上润湿状态不意图。
[0035] 图7为非湿相驱替湿相过程中三角形毛细管截面中的流体分布。
【具体实施方式】
[0036] 结合附图及实施例对本发明作进一步说明：
[0037] 如图1所示，利用一种考虑粗糙表面特征的不规则岩石毛管力确定方法，其步骤 如下：
[0038] (1)采用基于CT图像构建孔隙网络模型的方法对岩石孔隙截面进行形状等价，得 到等价形状截面的周长、面积及内角角度。
[0039] 其中，所述基于CT图像构建孔隙网络模型的方法为现有方法，具体可参见参考 文献：Hou Jian, Li Zhenquan, Zhang Sunkang, Cao Xulong, Song Xinwang, Gao Debo. Experiment and simulation study on construction of three-dimensional network model. Science in China，Series G，2008, 51 (9) : 1761-1774. 〇
[0040] (2)利用三维激光表面测量仪扫描岩石表面形貌并成像。图2为某岩石粗糙表面 三维激光表面测量仪成像后的形貌图。
[0041] (3)在三维形貌图中分别截取N个XZ平面和N个YZ平面的岩石轮廓并编号为1， 2, 3…2N，N -般取10,统计凸起高度h、凸起宽度a、凹槽宽度d等粗糙表面平均特性参数。
[0042] 图3为某XZ平面截面轮廓示意图，其编号记为i，在示意图中Ox轴表示取样轮廓 的最小二乘中线，保证在取样长度内轮廓偏离该线的平方和最小。该截面的凸起高度匕为 2倍的粗糙度R al，Ral表示阴影部分的面积与取样长度的比值。该截面的凸起宽度再为截面 岩石轮廓最小二乘中线与波峰交线平均长度等于
凹槽宽度式为截面岩石轮廓最 小二乘中线与波谷交线的平均长度等于
那么对于2N个截面有：
[0044] (4)利用凸起高度、凸起宽度及凹槽宽度进行粗糙表面等价，同时确定表面平均粗 糙因子rf、接触面中凸起固体所占面积分数fs。图4为粗糙表面等价示意图，表面计算参数 为：
[0047] (5)利用光学接触角测量仪测量非湿相流体液滴在湿相流体环境下光滑岩石表面 的本征接触角，判断液滴在粗糙岩石表面上的润湿状态。图5为液滴在光滑表面上的接触 角测试示意图，湿相本征接触角为Θ。图6为粗糙表面上润湿状态判断标准，在C-B状态下 液滴未侵入表面凹槽，而在Wenzel状态下液滴会侵入表面凹槽，此时的接触角称为表观接 触角θ'
[0048] (6)根据液滴润湿状态建立毛管力方程。
[0049] 在非湿相驱替湿相过程中，图7为岩石孔隙截面通过孔隙网络模型等价后的三角 形毛细管截面中的流体分布示意图，白色代表非湿相，深色代表湿相，灰色代表粗糙体。其 中图7(a)为光滑毛细管截面流体分布，图7(b)为粗糙表面Wenzel状态下的截面流体分 布，图7 (c)为粗糙表面C-B状态下的截面流体分布。MS-P方法认为驱替前沿半径和截面角 隅界面曲率半径相等，在Wenzel润湿状态下，由能量守恒方程及杨-拉普拉斯公式建立毛 管力Wenzel方程组为：
[0051] 如果液滴润湿为C-B状态，建立的毛管力C-B方程组为：
[0053] 其中rd是驱替前沿半径，Arff是截面中非湿相所占面积，L nw是非湿相与湿相的界 面周长，Lns是非湿相与管壁接触周长，γ nw是非湿相与湿相间的界面张力，γ ns是湿相与管 壁间的界面张力，Θ是光滑表面下的本征接触角，0$是粗糙表面下的表观接触角。
[0054] 求解出巧后，利用P C= γ nw/rd即得到Wenzel或C-B润湿状态下的岩石毛管力。
【主权项】
1. 一种考虑粗糖表面特征的岩石毛管力确定方法，其特征在于W下步骤： (1) 采用基于CT图像构建孔隙网络模型的方法对岩石孔隙截面进行形状等价，得到等 价形状截面的周长、面积及内角角度。 (2) 利用Ξ维激光表面测量仪扫描岩石表面形貌并成像。 (3) 在Ξ维形貌图中分别截取N个XZ平面和N个YZ平面的岩石轮廓并编号为1，2,3··· 2N，N-般取10,统计凸起高度h、凸起宽度a、凹槽宽度d等粗糖表面平均特性参数。 (4) 利用凸起高度、凸起宽度及凹槽宽度进行粗糖表面等价，同时确定表面平均粗糖因 子ff、接触面中凸起固体所占面积分数片。 (5) 利用光学接触角测量仪测量非湿相流体液滴在湿相流体环境下光滑岩石表面的本 征接触角，判断液滴在粗糖岩石表面上的润湿状态。 (6) 根据润湿状态建立相应的毛管力方程确定岩石毛管力。2. 如权利要求1中的一种考虑粗糖表面特征的岩石毛管力确定方法，其特征在于，其 所述步骤（3)中的粗糖表面特性参数统计方法如下：其中Ru为第i个截面岩石轮廓粗糖度，瑪为第i个截面岩石轮廓最小二乘中线与波峰 交线平均长度，巧为第i个截面岩石轮廓最小二乘中线与波谷交线的平均长度。3. 如权利要求1中的一种考虑粗糖表面特征的岩石毛管力确定方法，其特征在于，其 所述步骤（4)中的等价的粗糖表面特性参数确定公式如下：4. 如权利要求1中的一种考虑粗糖表面特征的岩石毛管力确定方法，其特征在于，其 所述步骤化）中的毛管力方程如下： 如果液滴润湿为Wenzel状态：其中rd是驱替前沿半径，Apff是截面中非湿相所占面积，L。点非湿相与湿相的界面周 长，Lm是非湿相与管壁接触周长，γ是非湿相与湿相间的界面张力，γΜ是湿相与管壁间 的界面张力，Θ是光滑表面下的本征接触角，是粗糖表面下的表观接触角。
【专利摘要】本发明公开了微观渗流力学研究领域中一种考虑粗糙表面特征的岩石毛管力确定方法，其主要步骤包括：采用基于CT图像构建孔隙网络模型的方法对岩石孔隙截面进行形状等价并统计截面几何参数；扫描岩石表面并进行三维形貌成像；截取岩石三维形貌图轮廓，统计粗糙表面平均特性参数；进行粗糙表面等价；测量岩石本征接触角，判断液滴在粗糙表面下的润湿状态；根据润湿状态建立毛管力方程确定毛管力。该方法充分考虑了粗糙表面对岩石微观结构及润湿性的影响，适用于不规则截面孔隙的岩石毛管力的确定，提高了毛管力的确定精度，从而能准确预测多孔介质渗流参数。
【IPC分类】G01N15/08
【公开号】CN105319155
【申请号】CN201510818606
【发明人】侯健, 韦贝, 姚军, 付红斐, 刘建林, 刘东滢, 杨勇, 曹绪龙, 郭兰磊, 于波, 杜庆军, 刘永革, 李淑霞, 姚传进, 陆努
【申请人】中国石油大学(华东)
【公开日】2016年2月10日
【申请日】2015年11月23日