裂缝性碳酸盐岩基质与裂缝耦合作用表征方法技术

本发明专利技术提供一种裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，该裂缝性碳酸盐岩基质与裂缝耦合作用表征方法包括：步骤1，构建不同裂缝开度的平板裂缝模型；步骤2，将不同裂缝开度的平板裂缝模型施加到CT扫描碳酸盐岩基质数字岩心中构建裂缝性碳酸盐岩三维数字岩心；步骤3，利用格子玻尔兹曼方法模拟三维数字岩心渗流场分布，并计算三维数字岩心模型的渗透率；步骤4，引入渗透率耦合系数，分析渗透率结果并确定渗透率耦合系数表达式。该裂缝性碳酸盐岩基质与裂缝耦合作用表征方法解决裂缝性碳酸盐岩渗流耦合作用对总渗透率造成的影响，更真实的反映碳酸盐岩裂缝性储层流体流动能力，为测井评价储层以及油田开发提供帮助。

【技术实现步骤摘要】
裂缝性碳酸盐岩基质与裂缝耦合作用表征方法
本专利技术涉及油田开发
，特别是涉及到一种裂缝性碳酸盐岩基质与裂缝耦合作用表征方法。
技术介绍
随着我国勘探开发程度的不断加深，近年来新探明的油藏类型越来越复杂，特别是裂缝性油藏的探明储量不断增加。裂缝性碳酸盐岩油藏由于天然裂缝的存在，构成了双重介质储层，在勘探开发中裂缝性碳酸盐岩储层渗透性参数对整个油藏的评价及解释具有重要的意义。裂缝性碳酸盐岩储层的孔隙空间由原生粒间孔隙和次生裂缝构成，裂缝性碳酸盐岩储层渗透率定义与单孔隙介质渗透率定义一致，但可以进一步细化为基质渗透率和裂缝渗透率。基质渗透率与单孔隙介质渗透率相同，裂缝渗透率受多种因素影响，主要包括裂缝孔隙度、裂缝张开度、裂缝密度、裂缝倾角和裂缝长度等。由于实际裂缝储集层中裂缝的分布极为复杂，要建立双重介质的渗透率模型，国内外常将裂缝系统进行简化，建立简化模型。裂缝储集层的简化模型主要有平行板模型，Kazemi模型和Warren-Root模型。平行板模型将实际模型简化为由水平基质层和裂缝相互交替组成，但不考虑基质渗透率；Kazemi模型将实际模型简化为由水平基质层和裂缝相互交替组成，其总渗透率为基质渗透率和裂缝渗透率之和；Warren-Root模型将双重孔隙介质油藏简化为正交裂缝切割基质岩块呈六面体的地质模型，裂缝方向与主渗透率方向一致，并假设裂缝宽度为一个常数，其总渗透率没有考虑基质渗透率。简化情况下推导得到的渗透率模型不能反映真实情况下的渗透率，只是一种理论性模型，有的学者提出了更一般化的渗透率模型。Parsons(1966)提出了考虑了裂缝的数量和裂缝与流体渗流方向间夹角的双重介质渗透率模型；美国学者Snow(1969)导出了具有可变方向的N条裂缝组成岩石的渗透率张量形式的模型；Oda(1985)引入裂缝张量的概念，导出了球形域的渗透率张量公式；Doolin和Mauldon(1996)建立了含有正交裂缝网络二维岩层的渗透率模型。但无论是渗透率简化模型还是一般化模型，国内外都把研究裂缝性碳酸盐岩储层渗透率模型的重点放到不同裂缝参数下的裂缝渗透率公式的推导，并没有考虑裂缝与岩石基质之间的耦合作用。而通过实验发现，当裂缝开度为10um时，如果不考虑裂缝与岩石基质之间的耦合作用，计算得到的渗透率结果相对于实际渗透率误差高达25.1％。将这些理论模型应用到测井解释中，由于没有考虑裂缝与基质之间耦合作用，会导致储层评价准确度降低。为了得到更符合实际的裂缝性碳酸盐岩储层渗透率模型，提高渗透率测井解释的精度，必须要考虑耦合作用。为此我们专利技术了一种新的裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，解决了以上技术问题。
技术实现思路
本专利技术的目的是提供一种通过研究不同基质孔隙结构参数，不同裂缝-孔隙组合下耦合作用的变化，构建裂缝性碳酸盐岩储层基质与裂缝耦合作用表征方法。本专利技术的目的可通过如下技术措施来实现：裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，该裂缝性碳酸盐岩基质与裂缝耦合作用表征方法包括：步骤1，构建不同裂缝开度的平板裂缝模型；步骤2，将不同裂缝开度的平板裂缝模型施加到CT扫描碳酸盐岩基质数字岩心中构建裂缝性碳酸盐岩三维数字岩心；步骤3，利用格子玻尔兹曼方法模拟三维数字岩心渗流场分布，并计算三维数字岩心模型的渗透率；步骤4，引入渗透率耦合系数，分析渗透率结果并确定渗透率耦合系数表达式。本专利技术的目的还可通过如下技术措施来实现：在步骤1中，构建完全为骨架的三维数据体，该三维数据体的体素个数及分辨率与建立的岩心基质三位数字岩心的参数一致；采用控制变量方法，只改变裂缝开度，固定其它裂缝参数，基于数字岩心叠加方法，构建不同裂缝开度的光滑平板裂缝模型。在步骤2中，基于数字岩心叠加方法，分别将构建的不同裂缝开度平板裂缝模型施加到CT扫描的不同孔隙结构的碳酸盐岩基质三维数字岩心中，构建裂缝性碳酸盐岩三维数字岩心。在步骤3中，采用格子玻尔兹曼方法中的D3Q19模型模拟岩石的渗流特性，D3Q19模型粒子分布函数的演化方程：等式中i的范围是[0,18]，表示19个不同的方向，左边为传播步，右边为碰撞步，fi(x,t)是t时刻在x位置第i个方向上流体粒子的传播函数，fi(x+eiΔt,t+Δt)是t+Δt时刻在x位置第i个方向上流体粒子的传播函数；在D3Q19模型中，τ是弛豫时间常量，fieq(x,t)是t时刻相同位置处流体粒子的平衡态分布函数，表达式为：式中，ρ是流体密度，c＝Δx/Δt.ωi是权重因子，ω0＝1/3,ωi＝1,2…6＝1/18,ωi＝7,8…18＝1/36,ei是速度失量，i表示每个不同的方向，u是根据动量守恒定律和质量守恒得到的流体线性宏观速度，单位是lu/ts，u可以写成：给系统施加一定的宏观压力梯度▽P使其开始演化，碰撞步和传播步重复执行直到粒子分布函数达到稳定态，局部流量可以通过下式计算得到：由达西定律得：式中，<q>表示沿着x方向通过岩心的流量，单位是cm3/s；dP/dx表示x方向压力梯度，单位是MPa/cm；A表示垂直于x方向岩心的横截面积，单位是cm2；η表示流体的动力粘度系数，单位是MPa；K表示在x方向岩心渗透率，单位是um2；通过计算每个格点处的宏观流速，得到三维空间渗流场分布，通过求解达西定律各个参数，进而得到某个三维数字岩心的渗透率。在步骤3中，利用格子玻尔兹曼方法，对构建的三维数字岩心进行渗流特性模拟，分别对构建的三维数字岩心施加相同的宏观压力梯度开始进行渗流模拟，直到粒子分布函数达到稳定态，得到渗流稳定态的三维空间渗流场分布，并根据出口端流量计算每个数字岩心的渗透率。在步骤4中，裂缝与基质之间的渗流耦合作用受裂缝开度的影响最大，而且与岩石基质物性相关，渗流耦合系数kx模型的形式为：kx＝pHf-q式中，Hf为裂缝开度，单位μm；参数p，q反映岩心基质物性，通过数字岩心方法进行拟合求取。本专利技术中的裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，研究裂缝性碳酸盐岩储层基质与裂缝耦合作用，填补目前裂缝性碳酸盐岩储层基质与裂缝耦合作用的空白，提出裂缝性碳酸盐岩基质与裂缝耦合作用表征方法。以往的裂缝性碳酸盐岩储层渗流研究都只单纯针对基质或者裂缝，很少考虑渗流在基质与裂缝同时存在时耦合作用，而基质与裂缝之间的耦合作用的存在一定程度上影响了流体通过基质、裂缝时的流速，改变了储层渗透率大小，通过耦合作用定量表征可以更真实的反映碳酸盐岩裂缝性储层流体流动能力。该方法能够定量表征裂缝性碳酸盐岩基质与裂缝渗流耦合作用，填补目前碳酸盐岩裂缝性储层基质与裂缝渗流耦合作用的空白，解决裂缝性碳酸盐岩渗流耦合作用对总渗透率造成的影响，更真实的反映碳酸盐岩裂缝性储层流体流动能力，为测井评价储层以及油田开发提供帮助。附图说明图1为本专利技术的裂缝性碳酸盐岩基质与裂缝耦合作用表征方法的一具体实施例的流程图；图2为本专利技术的一具体实施例中数字岩心基质渗流场分布图；图3为本专利技术的一具体实施例中10μm裂缝开度的平板裂缝模型渗流场分布图；图4为本专利技术的一具体实施例中10μm裂缝开度下裂缝性碳酸盐岩三维数字岩心渗流场分布图；图5为本专利技术的一具体实施例中数字岩心施加裂缝开度与耦合渗透率系数之间的关系图；图6为本专利技术的一具体实施例本文档来自技高网...

【技术保护点】
裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，其特征在于，该裂缝性碳酸盐岩基质与裂缝耦合作用表征方法包括：步骤1，构建不同裂缝开度的平板裂缝模型；步骤2，将不同裂缝开度的平板裂缝模型施加到CT扫描碳酸盐岩基质数字岩心中构建裂缝性碳酸盐岩三维数字岩心；步骤3，利用格子玻尔兹曼方法模拟三维数字岩心渗流场分布，并计算三维数字岩心模型的渗透率；步骤4，引入渗透率耦合系数，分析渗透率结果并确定渗透率耦合系数表达式。

【技术特征摘要】
1.裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，其特征在于，该裂缝性碳酸盐岩基质与裂缝耦合作用表征方法包括：步骤1，构建不同裂缝开度的平板裂缝模型；步骤2，将不同裂缝开度的平板裂缝模型施加到CT扫描碳酸盐岩基质数字岩心中构建裂缝性碳酸盐岩三维数字岩心；步骤3，利用格子玻尔兹曼方法模拟三维数字岩心渗流场分布，并计算三维数字岩心模型的渗透率；步骤4，引入渗透率耦合系数，分析渗透率结果并确定渗透率耦合系数表达式。2.根据权利要求1所述的裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，其特征在于，在步骤1中，构建完全为骨架的三维数据体，该三维数据体的体素个数及分辨率与建立的岩心基质三位数字岩心的参数一致；采用控制变量方法，只改变裂缝开度，固定其它裂缝参数，基于数字岩心叠加方法，构建不同裂缝开度的光滑平板裂缝模型。3.根据权利要求1所述的裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，其特征在于，在步骤2中，基于数字岩心叠加方法，分别将构建的不同裂缝开度平板裂缝模型施加到CT扫描的不同孔隙结构的碳酸盐岩基质三维数字岩心中，构建裂缝性碳酸盐岩三维数字岩心。4.根据权利要求1所述的裂缝性碳酸盐岩基质与裂缝耦合作用表征方法，其特征在于，在步骤3中，采用格子玻尔兹曼方法中的D3Q19模型模拟岩石的渗流特性，D3Q19模型粒子分布函数的演化方程：等式中i的范围是[0,18]，表示19个不同的方向，左边为传播步，右边为碰撞步，fi(x,t)是t时刻在x位置第i个方向上流体粒子的传播函数，fi(x+eiΔt,t+Δt)是t+Δt时刻在x位置第i个方向上流体粒子的传播函数；在D3Q19模型中，τ是弛豫时间常量，fieq(x,t)是t时刻相同位置处流体粒子的平衡态分布函数，表达式为：

【专利技术属性】
技术研发人员：乌洪翠，张丽艳，孙建孟，闫伟超，隋淑玲，杜玉山，赵国欣，王筱文，王伟，巴志明，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司勘探开发研究院，
类型：发明
国别省市：山东,37