本发明涉及一种石油天然气储集层内流体性质的评价方法，具体为一种基于灰色关联分析识别储层流体性质的方法。
背景技术：
目前，录井油气层解释评价一般以气测解释为基础，采用岩屑录井、定量荧光目、分析化验等多参数综合应用，分区块、分层位建立评价模型或评价图版的方法进行储层流体性质判识。
常规解释方法一般依据单项技术的全烃值高低、峰型饱满程度、组分等特征进行单参数储层流体性质评价，但由于地质条件的复杂性，不同层位、不同深度、不同的油气水性质的差异，随着油田地质勘探目标的隐蔽性、复杂性增加，仅依靠单项技术给油气层解释评价带来了的难度。有必要多项录井资料综合起来求取权重，继而实现综合评价储层流体性质。
灰色关联分析是根据序列的曲线几何形状的相似程度来判断其联系是否紧密，曲线越接近，相应序列之间关联度越大，反之越小。这种方法对样本量多寡和样本有无规律都同样适应，且计算量小、应用方便，更不会出现量化结果与定性分析结果不符的情况。
权系数是某一评价因素在决定总体特征时所占有的重要性程度，而不同地区不同层位各参数所占权重也不一样，因此，权重确定的合理性直接关系到评价结果的可靠性，这是油气层解释评价的一个无法回避的重要难题。
在储集层流体性质评价中，各参数对解释评价结果的影响主要体现在两个方面：一是评价者对各指标的重视程度不同，二是指标本身在评价过程中传给被评判体的信息量不同。各指标的权重大小可以反映用该指标单独评价结果的准确程度，并能良好反映出所选指标的实用性及其合理性。因此，合理确定评价指标的权重系数，客观反映它在综合评价中的重要性，是直接提高解释符合率的有效途径。
灰色关联分析的基本思想是根据序列曲线几何开关的相似程度来判断其联系是否紧密。曲线越接近，相应序列之间关联度就越大，该子因素对母因素的影响就越大，因此关联度与权重在基本意义上是相通的。
目前尚未见到将灰色关联分析方法综合运用到多参数储层流体性质评价中。
技术实现要素：
本发明的目的是：提供一种基于灰色关联分析的，适用于多参数综合评价储层流体性质的方法，弥补单参数评价储层油气水性质的局限性和多解性，提高录井油气水层解释符合率，为试油选层提供可靠支持。
为此，本发明的技术方案如下：
一种基于灰色关联分析识别储层流体性质的方法，包括如下步骤：
1)收集待分析地层及与其属于同一油气藏的已有试油结论层的气测、热解、定量荧光、孔隙度、渗透率和声波时差六项参数值；
优选，所述已有试油结论层包括含油水层、干层、油水同层和油层；
2)分别将所述六项参数值进行无量纲处理，具体包括以下两步：
①计算各项参数的极差，即各项参数的最大值与最小值之差；
②将各层的六项参数分别减去各项参数的最小值，再除以该项参数的极差，得到处理后的六项参数值；
3)求关联度；
基于步骤2)得到的处理后的六项参数值，以气测值为母因素，以热解和定量荧光为子因素，分别求各层各子因素与母因素的关联系数；
基于步骤2)得到的处理后的六项参数值，以孔隙度为母因素，以渗透率和声波时差为子因素，分别求各层各子因素与母因素的关联系数；
优选，求取各子因素与母因素的关联度的方法如下：
a)求取各层子因素与母因素的关联系数的公式如下：
其中：将与待分析地层属于同一油气藏的已有试油结论层分别标记为1至n层；p代表1至n中的任意一层，即p＝1，2，…，n；
式中，δi(p)＝▏x0ˊ(p)-xiˊ(p)▏；
x0ˊ(p)---第p层母因素值；
xiˊ(p)---第p层xi子因素值，i＝1代表第一个子因素，i＝2代表第二个子因素；
δi(p)---第p层xi子因素与母因素绝对差值；
m＝maxδi(p)；m＝minδi(p)；
m---xi子因素与母因素绝对差值最大值；
m---xi子因素与母因素绝对差值最小值；
ri(p)---第p层xi个参数的关联系数；
b)各子因素与母因素的关联度的计算方法如下：
其中：ri---xi子因素与母因素之间的关联度；
同时，母因素与母因素的关联度计为1；
4)求权重系数
计母因素、x1子因素、x2子因素的权重系数为a、b、c；
a＝1/(1+r1+r2)；b＝r1/(1+r1+r2)；c＝r2/(1+r1+r2)
求取气测、热解、定量荧光、孔隙度、渗透率和声波的六项参数的权重系数；
5)求油气丰度指数和物性指数；
p层的油气丰度指数为步骤2)处理后的气测、热解、定量荧光三项参数分别乘以步骤4)得到的各项参数的权重后求和得到的值；
p层的物性指数为步骤2)处理后的孔隙度、渗透率、声波时差三项参数分别乘以步骤4)得到的各项参数的权重后求和得到的值；
6)以p层的油气丰度指数为横轴，以p层的物性指数为纵轴，或者以p层的物性指数为横轴，以p层油气丰度指数的为纵轴；结合已有试油结论层的性质，构建二维灰色关联识别储层流体性质的解释图版，分别标记不同区域，具体可以为含油水层、干层、油水同层和油层的区域；
7)利用步骤2)得到的待分析地层处理后的气测、热解、定量荧光三项参数分别乘以步骤4)得到的各参数的权重后求和得到的待分析地层的油气丰度指数；
利用步骤2)得到的待分析地层处理后的孔隙度、渗透率、声波时差三项参数分别乘以步骤4)得到的各参数的权重后求和得到的待分析地层的物性指数；
将得到的待分析地层的油气丰度指数和物性指数在步骤6)构建的解释图版上描点，该点对应的区域即代表储层流体性质。
与现有技术相比，该技术方案具有如下优点：
这种基于灰色关联分析的识别储层流体性质的方法评价油气层评价界线明确、应用简单，结果标准、客观、科学、唯一，避免了单参数解释评价中的多解性、差异性、交叉性，且解决了长期困扰解释评价人员的各参数间的数据归一问题。有效地解决了单参数评价储层油气水性质的不足，客观的反映了储层流体性质，提高了录井油气水层解释符合率。
附图说明
图1为具体实施方式中构建的二维灰色关联识别储层流体性质的解释图版。
具体实施方式
以下结合附图和具体实施方式对本发明的技术方案进行详细描述。
一种基于灰色关联分析识别储层流体性质的方法，包括如下步骤：
1)收集待分析地层及与其属于同一油气藏的已有试油结论层的气测、热解、定量荧光、孔隙度、渗透率和声波时差六项参数值；具体数据如表1所示，表中全烃代表气测数据，后文同此。
表1wa斜坡油气显示数据表
2)分别将所述六项参数值进行无量纲处理，具体包括以下两步：
①计算各项参数的极差，即各项参数的最大值与最小值之差；
②将各层的六项参数分别减去各项参数的最小值，再除以该项参数的极差，得到处理后的六项参数值；得到数据见表2；
表2wa斜坡油气显示无量纲化数据表
3)求关联度；
基于步骤2)得到的处理后的六项参数值，以气测值为母因素，以热解和定量荧光为子因素，分别求各层各子因素与母因素的关联系数；
基于步骤2)得到的处理后的六项参数值，以孔隙度为母因素，以渗透率和声波时差为子因素，分别求各层各子因素与母因素的关联系数；
在本申请的一个具体实施方式中，求取各子因素与母因素的关联度的方法如下：
a)求取各层子因素与母因素的关联系数的公式如下：
其中：将与待分析地层属于同一油气藏的已有试油结论层分别标记为1至n层；p代表1至n中的任意一层，即p＝1，2，…，n；
式中，δi(p)＝▏x0ˊ(p)-xiˊ(p)▏；
x0ˊ(p)---第p层母因素值；
xiˊ(p)---第p层xi子因素值，i＝1代表第一个子因素，i＝2代表第二个子因素；
δi(p)---第p层xi子因素与母因素绝对差值；
m＝maxδi(p)；m＝minδi(p)；
m---xi子因素与母因素绝对差值最大值；
m---xi子因素与母因素绝对差值最小值；
ri(p)---第p层xi个参数的关联系数；
b)各子因素与母因素的关联度的计算方法如下：
其中：ri---xi子因素与母因素之间的关联度；
同时，母因素与母因素的关联度计为1；
利用以上方法得到的气测、热解、定量荧光的关联度分别为：1，0.6997，0.882；
求得的孔隙度、渗透率和声波时差的关联度分别为1，0.6316，0.5035。
4)求权重系数
计母因素、x1子因素、x2子因素的权重系数为a、b、c；
a＝1/(1+r1+r2)；b＝r1/(1+r1+r2)；c＝r2/(1+r1+r2)
得到气测、热解、定量荧光三项参数的权重系数依次为0.3873，0.2710，0.3416；
得到孔隙度、渗透率和声波三项参数的权重系数依次为0.4684，0.2958，0.2358；
5)求油气丰度指数和物性指数；
p层的油气丰度指数为步骤2)处理后的气测、热解、定量荧光三项参数分别乘以步骤4)得到的各项参数的权重后求和得到的值；
p层的物性指数为步骤2)处理后的孔隙度、渗透率、声波时差三项参数分别乘以步骤4)得到的各项参数的权重后求和得到的值；具体如表3所示；
表3wa斜坡灰色关联解释数据表
6)以p层的油气丰度指数为横轴，以p层的物性指数为纵轴，结合已有试油结论层的性质，构建二维灰色关联识别储层流体性质的解释图版，分别标记为含油水层、干层、油水同层和油层的区域，如图1所示；在具体绘图过程中，横坐标、纵坐标可以互换；
7)利用步骤2)得到的待分析地层处理后的气测、热解、定量荧光三项参数分别乘以步骤4)得到的各参数的权重后求和得到的待分析地层的油气丰度指数；
利用步骤2)得到的待分析地层处理后的孔隙度、渗透率、声波时差三项参数分别乘以步骤4)得到的各参数的权重后求和得到的待分析地层的物性指数；
将得到的待分析地层的油气丰度指数和物性指数在步骤6)构建的解释图版上描点，即图1中圆点，该点对应的区域为油水同层区，即代表储层流体性质为油水同层。
待测层试油情况：
待测层：解释井段2265-2269m，解释结论油水同层；试油井段2265-2267m，抽汲日产油3.17t，水1.19m3，试油结论为油水同层；与前期推测结论一致。