一种分析储层孔隙度的方法
【技术领域】
[0001] 本发明属于石油地球物理勘探技术领域，涉及一种分析储层孔隙度的方法，该方 法用于测井资料处理解释和储量分析中的孔隙度分析。
【背景技术】
[0002] 孔隙度分析是测井解释和储量计算中的重要环节，通常利用测井参数与岩心分析 孔隙度建立数理统计关系，并由此进行孔隙度分析。由于测井资料质量受井眼状况、测量仪 器工作状态等因素的影响，同时岩心分析实验也存在一定的测量误差，因此建立测井参数 与岩心分析孔隙度数理统计关系时常存在一些异常点，这些异常点不符合测井参数与岩心 分析孔隙度的整体关系趋势，如何准确地识别并放弃这些异常点，从而得到能代表测井参 数与岩心分析孔隙度整体关系趋势的特征数据，是孔隙度建模面临的一个重要问题。当前 测井解释人员常用的方法是先利用全部数据建模，然后凭借肉眼一次删除一个或几个在游 离于整体关系趋势之外的异常点，之后再利用剩下的数据点来建立测井参数与岩心分析孔 隙度数理统计关系，再次删除异常点，如此多次，直到得到自己觉得满意的孔隙度模型。这 个过程较繁复，又依赖于测井解释人员的判断。
[0003] 因此目前需要研究开发一种方法，能快速自动识别异常点，挑选出能代表测井参 数与岩心分析孔隙度整体趋势的特征数据来准确建立测井参数与岩心分析孔隙度数理统 计关系，从而能够由测井参数准确的对孔隙度进行分析。
【发明内容】
[0004] 本发明所要解决的技术问题是针对现有技术的不足，提供一种分析储层孔隙度的 方法。该方法在建立储层的岩心分析孔隙度和测井参数的对应关系过程中快速识别并排除 异常点，挑选能代表测井参数与岩心分析孔隙度整体趋势的特征数据快速建立较为准确的 建立储层的岩心分析孔隙度和测井参数的对应关系，并利用该对应关系通过待解释储层的 测井参数来准确分析待解释储层的孔隙度。
[0005] 为此，本发明提供了一种分析储层孔隙度的方法，其包括：
[0006] 步骤A，通过测井资料获取储层的测井参数；
[0007] 步骤B，测量储层中n个深度的岩心分析孔隙度；
[0008] 步骤C，建立储层的岩心分析孔隙度与测井参数的对应关系；
[0009] 步骤D，利用步骤C中所建立的储层的岩心分析孔隙度与测井参数的对应关系通 过测井参数来确定待解释储层的孔隙度。
[0010] 本发明中，所述测井参数包括测井实测值和/或测井计算值。
[0011] 所述测井实测值包括声波时差测井值、密度测井值、中子测井值中的至少一个。所 述测井计算值由声波时差测井值、密度测井值、中子测井值中的至少一个计算得到。例如， 测井计算值包括但不限于由声波时差测井值计算得到的声波孔隙度、密度测井值计算得到 的密度孔隙度、中子测井值计算得到的中子孔隙度、中子测井值与密度测井值(或声波时差 测井值）计算得到的平均孔隙度等。
[0012] 本发明中，n为正整数。
[0013] 根据本发明，步骤C包括：
[0014] 步骤K，以测井参数为横坐标，以岩心分析孔隙度为纵坐标，建立平面直角坐标系， 并将n对岩心分析孔隙度和相同深度的测井参数所对应的n个数据点绘制于坐标系内，并 且所述数据点的坐标为（ Xj，yP，
[0015] 其中，Xj为测井参数值；
[0016] y」为岩心分析孔隙度；
[0017] j=l、2、…、n_i ;
[0018] i=0、l、2、...n-1 ;
[0019] 步骤L，对n-i对岩心分析孔隙度和相同深度的测井参数所对应的n-i个数据点进 行线性拟合获得如式（I )所示的拟合直线；
[0020] y=ai ? x+bi 式（I )
[0021] 其中，ai为所述拟合直线的斜率，bi为所述拟合直线的截距，ai和bi均为步骤L 直线拟合过程中直接获得；
[0022] 步骤M，通过式（II)计算各数据点与步骤L中的拟合直线的距离；
【主权项】
1. 一种分析储层孔隙度的方法，其包括： 步骤A，通过测井资料获取储层的测井参数； 步骤B，测量储层中η个深度的岩心分析孔隙度； 步骤C，建立储层的岩心分析孔隙度与测井参数的对应关系； 步骤D，利用步骤C中所建立的储层的岩心分析孔隙度与测井参数的对应关系通过测 井参数来确定待解释储层的孔隙度。
2. 根据权利要求1所述的方法，其特征在于，所述测井参数包括测井实测值和/或测井 计算值。
3. 根据权利要求2所述的方法，其特征在于，所述测井实测值包括声波时差测井值、密 度测井值、中子测井值中的至少一个。
4. 根据权利要求2所述的方法，其特征在于，所述测井计算值由声波时差测井值、密度 测井值、中子测井值中的至少一个计算得到。
5. 根据权利要求1-4中任意一项所述的方法，其特征在于，步骤C包括： 步骤Κ，以测井参数为横坐标，以岩心分析孔隙度为纵坐标，建立平面直角坐标系，并将 η对岩心分析孔隙度和相同深度的测井参数所对应的η个数据点绘制于坐标系内，并且所 述数据点的坐标为（Xj，yj)， 其中，&为测井参数值； yj为岩心分析孔隙度； j=l、2、…、n_i ; i=0、1、2、…n_l ; 步骤L，对n-i对岩心分析孔隙度和相同深度的测井参数所对应的n-i个数据点进行线 性拟合获得如式（I )所示的拟合直线； Y=SLi · x+bj 式（I ) 其中，ai为所述拟合直线的斜率，bi为所述拟合直线的截距，ai和bi均为步骤L直线 拟合过程中直接获得； 步骤M，通过式（II)计算各数据点与步骤L中的拟合直线的距离；
步骤N，依据各数据点与步骤L中的拟合直线的距离判断是否存在异常点，并确定储层 的岩心分析孔隙度与测井参数的对应关系。
6. 根据权利要求5所述的方法，其特征在于，在步骤N中，依据各数据点与步骤L中的 拟合直线的距离对n-i个数据点进行排序，选出距离排序最大的数据点和距离排序中部的 数据点，通过式(III)计算距离排序最大的数据点和距离排序中部的数据点的距离倍数，然 后将其与距离倍数界限值进行比较，并且，当所计算的距离倍数大于距离倍数界限值时，则 判断该距离排序最大的数据点为异常点，剔除该数据点，i值加1，重复步骤L到N，直到所 计算的距离倍数大于距离倍数界限值，则判断为无异常数据点，步骤L中所拟合的直线即 为储层的岩心分析孔隙度与测井参数的对应关系； l=lenmax/len50 式(III) 其中，1为距离最大的数据点的距离倍数； Ienmax为距离最大的数据点的距离； len50为排序在中部的数据点的距离,该数据点的排序为（n-i) /2。
7.根据权利要求6所述的方法，其特征在于，所述距离倍数界限值为2-5。
【专利摘要】本发明涉及一种分析孔隙度的方法。该方法在建立已知储层的测井参数和岩心分析孔隙度的对应关系过程中快速识别并排除异常点，挑选能代表测井参数与岩心分析孔隙度整体趋势的特征数据快速建立较为准确的建立已知储层的测井参数和岩心分析孔隙度的对应关系，并利用该关系，通过待解释储层的测井参数来准确分析待解释储层的孔隙度。由于本发明方法每次只删除一个异常数据，不影响测井参数与岩心分析孔隙度的整体关系趋势，避免了误删数据，提高了工作效率。
【IPC分类】E21B47-08, E21B49-00
【公开号】CN104863572
【申请号】CN201410058847
【发明人】张军, 李军, 武清钊, 苏俊磊, 南泽宇, 王晓畅
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油勘探开发研究院
【公开日】2015年8月26日
【申请日】2014年2月20日