一种考虑动态毛管力的数值模拟方法技术

本发明专利技术公开了一种考虑动态毛管力的数值模拟方法，包括：对岩心设置多个测点，通过动态毛管力测试方法获得各测点下不同驱替速度下的动态毛管力曲线；利用渗流力学方法推导得到动态毛管力与含水饱和度、流速之间的函数关系，使用此函数关系对实验测试得到的动态毛管力曲线进行拟合，得到函数关系中的待定系数；根据函数关系对基本油相方程和基本水相方程进行差分，得到油相方程的差分离散方程，并且得到水相方程的与动态毛管力相关的差分离散方程，对油相方程的差分离散方程和水相方程的与动态毛管力相关的差分离散方程进行全隐式迭代求解，得到线性化的系数矩阵。本发明专利技术解决了现有技术中无法获得线性化的系数矩阵的问题。

A numerical simulation method considering dynamic capillary force

The invention discloses a method of numerical simulation, considering the dynamic capillary force includes: setting a plurality of measuring points on the core, obtained by the testing methods of dynamic capillary force of each measuring point under different displacement dynamic capillary pressure curve speed; function relation between seepage mechanics are derived using dynamic capillary force and water saturation, velocity the dynamic Mao Guanli curve using this function to test the fitting by undetermined coefficient function relationship; according to the function of basic equations and basic equations of water oil phase difference, get the oil phase equation discrete difference equation, and the equations of water related dynamic capillary force differential separation the dispersion equation, the oil phase equation of discrete difference equation and the equations of water related dynamic capillary force difference discrete equation and iterative implicit solution The linearized coefficient matrix is obtained. The invention solves the problem that the coefficient matrix which can not be linearized is not available in the prior art.

【技术实现步骤摘要】
一种考虑动态毛管力的数值模拟方法
本专利技术涉及油藏数值模拟技术应用领域，尤其涉及一种考虑动态毛管力的数值模拟方法。
技术介绍
毛管力是毛细管压力的简称，是指毛细管中弯液面两侧两种流体(非湿相流体与湿相流体)的压力差。油气藏开发过程中，储层中的油、水、气是始终流动的，由于流体流动必然会造成动润湿滞后(流体流动速度大于三相周界移动速度引起的润湿角改变)，从而影响流体流动过程中毛管力的数值。渗流力学和油藏数值模拟过程中涉及到的毛管力，应当是与流体流动速度相关的动态毛管力。但是，目前毛管力的测试实验所应用的岩心夹持器都是一端进入高压流体，另一端封闭，测试过程中流体的流动方式与油藏条件下的流体流动方式截然不同，整个实验过程的压力变化体现的是孔喉半径的差异，其实质是“静态毛管力”，并非油藏条件下流体流动的毛管力。在水动力学领域，动态毛管力的研究成果较多，但都是测试的气水两相在填砂模型中的动态毛管力，从Camps-Roach测试得到的气水两相的实验结果可知静态毛管力与动态毛管力的存在明显差异。Kalaydjian利用经过改装的测试装置，测试得到了油水两相的动态毛管力的实验结果，从测试结果分析，静态毛管力与动态毛管力之间依然存在较大差异。很多学者利用物质平衡方法计算得到了动态毛管力数值，并对动态毛管力的影响因素开展了分析，认为动态毛管力不仅受控于孔隙结构、渗透率、流体性质等静态因素，而且受控于流体流动速度等动态因素。HassanizadehandGray提出动态毛管力与静态毛管力之间的差异可用如下的数学方程进行描述：式(1)中，Pdyn是动态毛管力；Pstat是静态毛管力(处于平衡状态的毛管力，即τ是动态系数，用于描述饱和度变化对毛管力的影响程度。已有大量学者研究了动态毛管力在不同条件下的变化规律和取值方法。Hassanizadeh利用室内实验结果，提出τ的取值范围是5×104～6×104Pa·s。Stauffer对大量室内实验的结果进行总结，提出了τ的具体计算公式：式中，α是无因次系数，通常取0.1；φ是孔隙度，单位为f；K是渗透率，单位为md；μ是湿相流体粘度，单位为mPa·s；g是重力加速度，单位为N/kg；ρ是湿相流体密度，单位为g/m3；Pd是孔隙进入压力，单位为Mpa；λ是粒径分布系数，单位为f。然而，目前应用较为广泛的油藏数值模拟方法在处理毛管力时，都认为毛管力仅仅是含水饱和度的函数，数值模拟所用的毛管力在取值时，参考的是压汞法、离心法或者半渗透隔板法所测试得到的毛管力，考虑到润湿滞后效应，将毛管力分为吸吮过程的毛管力和驱替过程的毛管力。基本渗流方程建立过程和基本渗流方程差分求解过程中，没有考虑流体流动速度和饱和度变化速度对毛管力的影响。Sandoval-Torres等利用有限元解法，将动态毛管力引入到计算木头真空干燥的数值模拟方程中，将数值模拟的计算结果与实验结果进行了对比，认为建立的数值模拟模型准确可靠。Shu等同样利用有限元解法，将动态毛管力引入到计算油水两相渗流的数值模拟方程中，计算得到了一维条件下考虑动态毛管力的数值模拟结果，分析了一维条件下动态毛管力对油水运移规律的影响。但是，油藏数值模拟器更多的选用有限差分方法，且油藏数值模拟针对的是宏观的三维空间的油水运移规律。上述两种模型并不能很好的用于油藏条件下的数值模拟计算。Tian等将动态毛管力引入到油藏数值模拟基本渗流方程中，并用有限差分方法进行了求解，但是求解方法选用的是IMPES(显式饱和度，隐式压力)方法，计算的稳定性和准确性不如全隐式求解，并且，计算结果仅与室内实验结果进行了对比，并未与实际生产数据进行对比。现有技术中在油藏数值模拟中不考虑动态毛管力的方式影响模拟结果的精确性。
技术实现思路
为了解决上述技术问题，本专利技术一种考虑动态毛管力的数值模拟方法。本专利技术提供的考虑动态毛管力的数值模拟方法，包括：步骤1，对岩心设置多个测点，通过动态毛管力测试方法获得各测点下不同驱替速度下的动态毛管力曲线；步骤2，利用渗流力学方法推导得到动态毛管力与含水饱和度、流速之间的函数关系，使用此函数关系对实验测试得到的动态毛管力曲线进行拟合，得到函数关系中的待定系数；步骤3，根据所述函数关系对基本油相方程和基本水相方程进行差分，得到油相方程的差分离散方程，并且得到水相方程的与动态毛管力相关的差分离散方程，对所述油相方程的差分离散方程和所述水相方程的与动态毛管力相关的差分离散方程进行全隐式迭代求解，得到线性化的系数矩阵。上述考虑动态毛管力的数值模拟方法还具有以下特点：所述步骤1中动态毛管力曲线的横坐标为含水饱和度Sw，纵坐标为动态毛管力数值vw表示水相驱替流速；所述步骤2中，所述函数关系为：所述函数关系中的待定系数为k1和k2；其中，指静态毛管力，αs是实验常数，μw是润湿相相粘度，是孔隙度，Ct是压缩系数，K是绝对渗透率，λ和pe是Brook–Corey模型中的毛管压力因子，ρw是水密度，g是重力加速度。上述考虑动态毛管力的数值模拟方法还具有以下特点：所述步骤3中水相方程的与动态毛管力相关的差分离散方程为:是二阶差商算子，是第n+1个时间步的传导率，是第n+1个时间步的势，是指第n+1时刻的水相注入速度，Vijk是指坐标为(i，j，k)网格的孔隙体积，Δtn+1第n+1个时间步的时间间隔，φ是孔隙度，ρw是指润湿相密度。上述考虑动态毛管力的数值模拟方法还具有以下特点：步骤3中对所述油相方程的差分离散方程和所述水相方程的与动态毛管力相关的差分离散方程进行全隐式迭代求解包括：设置参数势为是指第n时刻的非润湿相压力，是指第n个时间步的毛管力，x是指坐标x的方向，p是指压力，t是指时间，γw是指水的重度，D是指垂向高度；根据所述参数势对所述油相方程的差分离散方程和所述水相方程的与动态毛管力相关的差分离散方程进行全隐式迭代求解，得到线性化的系数矩阵。上述考虑动态毛管力的数值模拟方法还具有以下特点：所述动态毛管力测试方法包括：步骤一，在岩心不同测点上设置油相压力传感器和水相压力传感器，在岩心的油输入口连接油容器，在岩心的水输入口连接水容器，在实验用油容器和实验用水容器的输入端连接驱替泵；步骤二，确定油水比例，确定驱替泵的驱替速度；步骤三，分别在实验用油容器和实验用水容器放置符合所述油水比例的实验用油和实验用水；步骤四，将所述油容器的实验用油和所述水容器的实验用水驱替入所述岩心，在驱替稳定后，通过所述油相压力传感器和所述水相压力传感器实时检测并记录油相压力和水相压力；步骤五，改变油水比例和驱替速度，重复执行步骤二和步骤三，获得不同驱替速度条件下的油相压力、水相压力和饱和度变化，从而获得各测点下不同驱替速度下的动态毛管力曲线离散参数。本专利技术解决了现有技术中无法获得线性化的系数矩阵的问题。本专利技术将动态毛管力引入到油藏数值模拟基本渗流方程中，并选用全隐式的方法对饱和度和压力进行了差分求解，提高了模拟结果的精确性。申请人将计算结果分别与室内实验结果和实际油田生产数据进行了对比，验证了本方法的可靠性和准确性。附图说明图1是本专利技术中考虑动态毛管力的数值模拟方法的流程图。具体实施例为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的本文档来自技高网...

【技术保护点】
一种考虑动态毛管力的数值模拟方法，其特征在于，包括：步骤1，对岩心设置多个测点，通过动态毛管力测试方法获得各测点下不同驱替速度下的动态毛管力曲线；步骤2，利用渗流力学方法推导得到动态毛管力与含水饱和度、流速之间的函数关系，使用此函数关系对实验测试得到的动态毛管力曲线进行拟合，得到函数关系中的待定系数；步骤3，根据所述函数关系对基本油相方程和基本水相方程进行差分，得到油相方程的差分离散方程，并且得到水相方程的与动态毛管力相关的差分离散方程，对所述油相方程的差分离散方程和所述水相方程的与动态毛管力相关的差分离散方程进行全隐式迭代求解，得到线性化的系数矩阵。

【技术特征摘要】
1.一种考虑动态毛管力的数值模拟方法，其特征在于，包括：步骤1，对岩心设置多个测点，通过动态毛管力测试方法获得各测点下不同驱替速度下的动态毛管力曲线；步骤2，利用渗流力学方法推导得到动态毛管力与含水饱和度、流速之间的函数关系，使用此函数关系对实验测试得到的动态毛管力曲线进行拟合，得到函数关系中的待定系数；步骤3，根据所述函数关系对基本油相方程和基本水相方程进行差分，得到油相方程的差分离散方程，并且得到水相方程的与动态毛管力相关的差分离散方程，对所述油相方程的差分离散方程和所述水相方程的与动态毛管力相关的差分离散方程进行全隐式迭代求解，得到线性化的系数矩阵。2.如权利要求1所述的考虑动态毛管力的数值模拟方法，其特征在于，所述步骤1中动态毛管力曲线的横坐标为含水饱和度Sw，纵坐标为动态毛管力数值vw表示水相驱替流速；所述步骤2中，所述函数关系为：所述函数关系中的待定系数为k1和k2；其中，指静态毛管力，αs是实验常数，μw是润湿相相粘度，是孔隙度，Ct是压缩系数，K是绝对渗透率，λ和pe是Brook–Corey模型中的毛管压力因子，ρw是水密度，g是重力加速度。3.如权利要求1所述的考虑动态毛管力的数值模拟方法，其特征在于，所述步骤3中水相方程的与动态毛管力相关的差分离散方程为:

【专利技术属性】
技术研发人员：王硕亮，于春磊，魏勇，于希南，张媛，付强，赵辉，李俊键，陈文斌，傅栋，王明，
申请(专利权)人：中国地质大学北京，
类型：发明
国别省市：北京,11