本发明专利技术公开了一种水驱气藏相对渗透率曲线高效构建与分析方法，方法包括建立毛细管压力函数关系式；气水相对渗透率函数关系式；获得归一化毛细管压力曲线；获得毛细管压力函数关系式中的待定系数以及气水相对渗透率曲线；计算对比分析不同毛细管压力函数关系式的拟合误差以及不同方法计算相渗曲线差异；将由压汞实验测试成果间接计算获取的相渗曲线与直接实验测试获得的相渗曲线进行对比检验。本发明专利技术在实验室相渗曲线测试难度大、对设备要求高，以及当前相渗曲线测试成果不足的情况下，具有显著优势与应用价值。可为水驱气藏工程计算、产能预测和数值模拟等方面工作提供可靠基础相渗资料，同时为同类水驱气藏的相渗曲线的间接计算提供理论指导。
【技术实现步骤摘要】
本专利技术涉及水驱气藏气水相对渗透率曲线间接获取方法的应用领域，具体涉及水驱气藏相对渗透率曲线高效构建与分析方法。
技术介绍
1、水驱气藏作为一种重要的油气资源类型，在全球广泛分布。在我国，水驱气藏主要分布在鄂尔多斯、四川、塔里木等盆地，占据了我国气藏数量的较高比例。高效开发这类气藏对缓解我国对清洁能源的迫切需求具有重要意义。深入了解水驱气藏储层的孔隙结构特征和流体渗流规律是开发该类油气资源的前提与基础。
2、气水相对渗透率曲线(以下简称“相渗曲线”)是气藏开发方案与调整方案设计的必备基础资料。通常，室内驱替实验是获取相渗曲线最直接的手段。对于深层水驱气藏，储层温度、压力普遍较高。在高温和高压条件下进行相渗曲线测试实验，对实验设备要求高，实验测试成本高，测试费时费力。相渗曲线测试实验组数量过少或相渗测试结果代表性不足，测试成果将难以真实反映气藏(或开发层段)气水两相流体渗流特性。将不具代表性、不准确的相渗曲线用于气藏动态分析、数值模拟及产能预测等工作会导致计算结果出现较大偏差，进而导致气藏开发人员对气藏开发状况形成错误的认识。
3、本专利技术提供水驱气藏相对渗透率曲线高效构建与分析方法，尤其是在相渗曲线测试难度大、测试成果不足的情况下，该方法为水驱气藏代表性相渗曲线的获取提供一种可靠的技术思路与解决方案。
技术实现思路
1、鉴于上述问题，本专利技术提供水驱气藏相对渗透率曲线高效构建与分析方法，旨在解决水驱气藏气水相对渗透率曲线的间接高效获取难题。
>2、为解决上述至少一种技术问题，本专利技术提供了水驱气藏相对渗透率曲线高效构建与分析方法，该方法包括：3、(1)建立毛细管压力函数关系式。具体来说，收集实验室压汞测试成果，分析各岩心压汞曲线，根据压汞曲线形态特征，建立毛细管压力函数关系式；
4、(2)建立气水相对渗透率函数关系式。具体来说，结合多孔介质流体渗流理论与毛细管压力函数，采用purcell方法与burdine修正方法，建立用于表征多孔介质驱替与渗吸过程中气水相对渗透率函数关系式。
5、(3)获得归一化毛细管压力曲线。具体来说，根据储层、气井、气藏性质认识，明确岩心孔渗参数，对岩心进行分类，并采用毛细管压力曲线归一化方法，获得不同类型储层、气井及整个气藏的归一化毛细管压力曲线。所获归一化毛细管压力曲线能够反映不同类型储层、气井及气藏的整体孔隙结构与渗流特性。
6、(4)获得毛细管压力函数关系式中待定系数以及气水相对渗透率曲线。具体来说，采用非线性回归方法，获得毛细管压力函数关系式中的待定系数，将待定系数代入由毛细管压力函数关系式推导获得的气水相对渗透率函数关系式中，获得能够反映不同类型储层、气井及目标气藏整体孔隙结构与渗流特性的气水相对渗透率曲线。
7、(5)计算对比分析不同毛细管压力函数关系式的拟合误差以及不同方法计算得到的气水相对渗透率曲线差异。具体来说，分别采用本专利技术毛细管压力函数关系式与气水相对渗透率函数关系式、行业应用较为广泛的corey毛细管压力函数关系式与brooks&corey相对渗透率函数关系式，以及采用基于burdine修正方法的数值积分手段，结合实验室压汞测试成果，计算对比分析不同毛细管压力函数关系式的拟合误差，进而分析不同函数关系式在毛细管压力曲线定量表征方面的差异，以及不同方法计算的相对渗透率曲线差异，并给出使用建议。
8、(6)将实验室由压汞测试成果间接计算获取的气水相对渗透率曲线与直接实验测试获得的气水相对渗透率曲线进行对比检验，讨论本专利技术提出方法的适用性，并给出使用建议。
9、本申请实施例提供的方法可以帮助技术人员获得能够准确可靠反映目标气藏、不同类型储层及气井真实渗流特性的气水相对渗透率曲线，为水驱气藏的工程计算、数值模拟、动态分析等工作提供可靠的基础相渗资料。
10、在一种实施例中，毛细管压力函数关系式，包括：
11、
12、或为如下所述毛细管压力函数关系式：
13、
14、其中：pc为毛细管压力，mpa；pc,max为最大毛细管压力，mpa；pc,min为最小毛细管压力，mpa；swn为归一化润湿相饱和度；a、b、c、α、β、γ为拟合函数的待定系数或实验拟合系数，α＝(pc,max/pc,min)a。
15、在一种实施例中，气水相对渗透率函数关系式，包括：
16、基于所述purcell方法的气水相对渗透率函数关系式，如下式：
17、
18、
19、基于所述burdine方法修正上述气水相对渗透率函数关系式，如下式：
20、
21、
22、其中，krw、krg分别为液相与气相的相对渗透率，无因次；b、c为拟合函数待定系数或实验拟合系数；swn为归一化润湿相饱和度。
23、在一种实施例中，所述毛细管压力曲线归一化方法，通过采用leverett无因次j(swn)函数计算各岩心压汞曲线中不同归一化润湿相饱和度swn下的j(swn)函数值，然后利用如下方程对j(swn)函数值非线性拟合，获得j(swn)函数非线性拟合曲线，表达式如下：
24、
25、式中：j(swn)为leverett无因次“j(swn)函数”；ε1、ε2为拟合方程待定系数或实验拟合系数；swn为归一化润湿相饱和度。
26、在一种实施例中，基于所述j(swn)函数值非线性拟合曲线，是通过采用非线性拟合方法获得j(swn)函数拟合方程的待定系数，然后代入所述j(swn)拟合函数方程，绘制j(swn)函数非线性拟合曲线，然后根据所述j(swn)函数非线性拟合曲线，将归一化润湿相饱和度swn等分为n份，每个等分点采用算数平均方法，得到不同归一化润湿相饱和度swn下的平均j(swn)函数值，计算式如下：
27、
28、基于上式所述平均j(swn)函数值，利用leverett提出的j(swn)函数方程，计算毛细管压力，并对归一化润湿相饱和度swn反归一化，得到润湿相饱和度sw，反归一化公式如下：
29、进汞：
30、
31、退汞：
32、
33、式中：sw为润湿相饱和度，％；swn为归一化润湿相饱和度，％；i为j(swn)函数拟合曲线上归一化润湿相饱和度swn的第i个等分点，个；k为岩心样品数，个；为归一化润湿相饱和度swn的第i个等分点对应的j(swn)函数平均值；sw,i为归一化润湿相饱和度swn第i个等分点值对应的sw值，％；swn,i为归一化润湿相饱和度swn第i个等分点值，％；为所有岩心样品残余润湿相饱和度(或最小润湿相饱和度)的平均值，％；
34、为所有岩心样品最大润湿相饱和度的平均值，％。
35、在一种实施例中，所述获取毛细管压力与相渗函数方程与曲线，包括：1)依据定量指标(主要孔隙度和渗透本文档来自技高网...
【技术保护点】
1.一种水驱气藏相对渗透率曲线高效构建与分析方法，其特征在于，包括：
2.根据权利要求1所述的方法，其特征在于，所述毛细管压力函数关系式，包括：
3.根据权利要求1所述的方法，其特征在于，所述气水相对渗透率函数关系式，包括：
4.根据权利要求1所述的方法，其特征在于，所述毛细管压力曲线归一化方法，通过采用Leverett无因次J(Swn)函数计算各岩心压汞曲线中不同归一化润湿相饱和度Swn下的J(Swn)函数值，然后利用如下方程对J(Swn)函数值非线性拟合，获得J(Swn)函数非线性拟合曲线，表达式如下：
5.根据权利要求4所述的方法，其特征在于，基于所述J(Swn)函数非线性拟合曲线，是通过采用非线性拟合方法获得J(Swn)函数拟合方程的待定系数，然后代入所述J(Swn)拟合函数方程，绘制J(Swn)函数非线性拟合曲线，然后根据所述J(Swn)函数非线性拟合曲线，将归一化润湿相饱和度Swn等分为n份，每个等分点采用算数平均方法，得到不同归一化润湿相饱和度Swn下的平均J(Swn)函数值，计算式如下：
6.根据权利要求1所述的方法，其特征在于，还包括：
...
【技术特征摘要】
1.一种水驱气藏相对渗透率曲线高效构建与分析方法，其特征在于，包括：
2.根据权利要求1所述的方法，其特征在于，所述毛细管压力函数关系式，包括：
3.根据权利要求1所述的方法，其特征在于，所述气水相对渗透率函数关系式，包括：
4.根据权利要求1所述的方法，其特征在于，所述毛细管压力曲线归一化方法，通过采用leverett无因次j(swn)函数计算各岩心压汞曲线中不同归一化润湿相饱和度swn下的j(swn)函数值，然后利用如下方程对j(swn)函数值非线性拟合，获得j(swn...
【专利技术属性】
技术研发人员：秦正山，何勇明，陈志豪，丁洋洋，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。