可变岩电耦合指数含水饱和度计算方法
【专利摘要】本发明公开了一种可变岩电耦合指数含水饱和度计算方法，包括如下步骤：a、建立新的含水饱和度计算模型；b、建立饱和度岩心数据库；c、对饱和度岩心数据库的信息进行挖掘，建立岩电耦合指数计算图版；d、对新输入的测井资料进行处理，得出新输入测井资料的岩电耦合指数；e、将步骤d得到的岩电耦合指数用于步骤a中的计算模型，计算得出新输入测井资料的含水饱和度。本发明首次提出了岩电耦合指数概念，并找到了相应的计算方法，达到了提高计算精度的目的。
【专利说明】可变岩电耦合指数含水饱和度计算方法
【技术领域】
[0001]本发明涉及一种可变岩电耦合指数含水饱和度计算方法，用于石油地质测井资料处理。
【背景技术】
[0002]石油地质测井资料处理有一项重要的内容，就是含水饱和度的计算。目前，计算模型多达30多种，这些计算模型均从阿尔奇公式演变而来，针对不同的地区，不同的地质体，这些计算模型演变的思路总体上离不开以下三点:①对阿尔奇公式做泥质校正对阿尔奇公式做矿化度校正对阿尔奇公式做温度校正。例如，2007年05期《测井技术》公开了含水饱和度指数η值计算方法探讨，研究了不同样品在不同实验条件下的含水饱和度指数η值的计算方法。根据岩电实验中含水饱和度Sw与电阻增大率I交会曲线的特征，分析了样品孔隙结构及含水饱和度对含水饱和度指数η值的影响，使用了 1-Sw分段回归法、n-Sw相关分析法以更精确地计算样品的含水饱和度指数η值。实验分析了不同的油藏参数与含水饱和度指数η值的关系，认为储层性质和流体性质是影响含水饱和度指数η值的关键参数，对实际油藏应通过多因素拟合的方法求解含水饱和度指数η值。
[0003]但现有的方法还存在如下不足:每做一项校正，都会给阿尔奇公式增加新的项或新的待定参数，其结果会导致计算模型的适应性变窄。另外，阿尔奇公式自身有四个待参数，它们都需要岩心实验数据来标定，如果增加了更多待定参数，会使计算模型计算结果的稳定性变差。理论和实验研究均表明，即使同种岩性，阿尔奇公式胶结指数(m，无量纲)和饱和度指数(n，无量纲)并不固定，而目前计算含水饱和度均采用固定的m，η值，这是阿尔奇公式计算精度不高的主要原因。
【发明内容】
[0004]本发明的目的在于克服现有得到含水饱和度的方法存在的上述问题，提供一种可变岩电耦合指数含水饱和度计算方法。本发明首次提出了岩电耦合指数概念，并得出相应的计算方法，最大限度地发挥含水饱和岩心实验数据库的作用，简化传统含水饱和度计算模型，同时提高了含水饱和度的精度。
[0005]为实现上述目的，本发明采用的技术方案如下:
一种可变岩电耦合指数含水饱和度计算方法，其特征在于，包括如下步骤:
a、建立新的含水饱和度计算模型；
b、建立饱和度岩心数据库；
c、对饱和度岩心数据库的信息进行挖掘，建立岩电耦合指数计算图版；
d、对新输入的测井资料进行处理，得出新输入测井资料的岩电耦合指数；
e、将步骤d得到的岩电耦合指数用于步骤a中的计算模型，计算得出新输入测井资料的含水饱和度。
[0006]所述步骤a中，建立新的含水饱和度计算模型过程如下:将胶结指数》、饱和度指数Il对含水饱和度&计算结果的影响耦合成一个指数，称之为岩电耦合指数，等价关系表示为:
【权利要求】
1.一种可变岩电耦合指数含水饱和度计算方法，其特征在于，包括如下步骤: a、建立新的含水饱和度计算模型； b、建立饱和度岩心数据库； c、对饱和度岩心数据库的信息进行挖掘，建立岩电耦合指数计算图版； d、对新输入的测井资料进行处理，得出新输入测井资料的岩电耦合指数； e、将步骤d得到的岩电耦合指数用于步骤a中的计算模型，计算得出新输入测井资料的含水饱和度。
2.根据权利要求1所述的可变岩电耦合指数含水饱和度计算方法，其特征在于:所述步骤a中，建立新的含水饱和度计算模型过程如下: 将胶结指数3B、饱和度指数?对含水饱和度尽汁算结果的影响耦合成一个指数，称之为岩电耦合指数，等价关系表示为:
3.根据权利要求1或2所述的可变岩电耦合指数含水饱和度计算方法，其特征在于:所述步骤b建立饱和度岩心数据库中，饱和度岩心数据库字段主要包括岩心实验孔隙度，岩心实验得到的含水饱和度，经过测井深度校正后得到的原状地层电阻率，自然伽玛值，补偿中子值，补偿声波以及补偿密度值。
4.根据权利要求1或2所述的可变岩电耦合指数含水饱和度计算方法，其特征在于:所述步骤c中，对饱和度岩心数据库的信息进行挖掘的方法如下: 先利用步骤a建立的计算模型，反演《 ;对于任意一个岩心样品，已知岩心实验含水饱和度Sb ,岩心实验孔隙度#,地层水电阻率Λ?，岩石孔隙结构弯曲度β ,然后经过测井深度校正找到对应的原状地层电阻率馬，通过数值反演，求解s值，计算式力:
5.根据权利要求4所述的可变岩电耦合指数含水饱和度计算方法，其特征在于:所述对数据点进行分类的步骤如下: al、构建由测井项目测量值形成的5维数学空间； a2、采用最短距离法对5维数学空间里的数据点进行聚类分析，实现对数据点的分类； a3、得到每一类的数据点的集合，通过公式3反演得到的岩电耦合指数? a4、得到每一类的重心，重心是指:该类所有数据点各个测井项目测量值的均值。
6.根据权利要求1、2或5所述的可变岩电耦合指数含水饱和度计算方法，其特征在于:所述步骤c中，建立每一类? = /(?9计算图版的步骤如下: Cl、利用每一类数据集合里的补偿声波测量值，根据威利平均公式，求出声波孔隙度Φα (单位:小数)； c2、根据声波孔隙度#?与》具有一一对应关系，建立该类的? = /(的函数关系式； c3、所有类的? = /--函数关系式就形成了计算图版。
7.根据权利要求6所述的可变岩电耦合指数含水饱和度计算方法，其特征在于:所述步骤d按下述步骤对新输入测井资料进行处理，步骤如下: dl、根据新输入测井资料的测井值归并到某一类； d2、然后按该类已经建立的? = /W)函数关系式计算出岩电耦合指数β。
8.根据权利要求7所述的可变岩电耦合指数含水饱和度计算方法，其特征在于:所述步骤dl中，新输入的测井资料是5个测井项目构成的5维空间当中的一个点，求出该点与上述每一类的重心之间的距离，该点与某一类的重心距离最短，则归并到那一类。
9.根据权利要求1、2、5、7或8所述的可变岩电耦合指数含水饱和度计算方法，其特征在于:所述步骤e中，原状地层电阻率、孔隙度均为新输入的测井资料，岩石孔隙结构弯曲度β，地层水电阻率Ab通过岩心实验分析的地层水样分析均已知，将步骤d得到的岩电耦合指数代入步骤a建立的含水饱和度计算模型即式，求出新输入测井资料类ωΨ的含水饱和度。
【文档编号】E21B47/00GK103527172SQ201310483685
【公开日】2014年1月22日 申请日期:2013年10月16日 优先权日:2013年10月16日
【发明者】王勇军, 罗利, 张树东, 甘秀娥, 陈邦定 申请人:中国石油集团川庆钻探工程有限公司