前判别函数的岩性识别正确率；
[0073] 选取子单元243用于在判别函数的岩性识别正确率大于或等于90%时，选取识别 函数为目标判别函数；
[0074] 剔除子单元244用于在判别函数的岩性识别正确率小于90%时，剔除岩性识别正 确率最低的岩性类别。
[0075] 在本发明提供的碳酸盐岩地层录井综合柱状图的校正系统中，测井曲线为能够反 映岩石物理性质的测井曲线，可以包括密度测井曲线、声波时差测井曲线、自然伽马测井曲 线、电阻率系列测井曲线以及中子孔隙度测井曲线等。
[0076] 本发明充分利用了钻井取心资料与测井资料，建立了岩性与测井响应之间的关 系，实现了基于Fisher判别分析的岩性再识别，对原有的精度较低的录井柱状图进行了校 正，提高了岩性剖面精度。
[0077] 为了对上述本发明提供的碳酸盐岩地层录井综合柱状图的校正方法及系统进行 更为清楚的解释，下面结合一个具体的实施例来进行说明，然而值得注意的是该实施例仅 是为了更好地说明本发明，并不构成对本发明不当的限定。
[0078] 塔里木盆地寒武系广泛发育海相碳酸盐岩沉积，并且普遍被白云岩化。目前寒武 系的油气勘探程度很低，钻井数量少。从地层录井综合柱状图看，寒武系碳酸盐岩的岩性基 本都定名为白云岩。如巴楚地区的和4井寒武系丘里塔格组发育730余米的碳酸盐岩沉 积，地层录井剖面上岩性十分单一，除发育50米厚的辉长岩外，其余均为碳酸盐岩，其中包 括35米厚的泥晶灰岩。单一的岩性剖面对于沉积相、层序地层以及储层研究都是不利的。 该井没有对灰岩段进行钻井取心，只对白云岩段进行了取心。从取心看，白云岩的岩性类别 丰富。同时该井进行了多种测井系列的测井，因此，可以利用岩心及测井资料对白云岩段的 岩性的识别，从而对地层录井综合柱状图进行校正。
[0079] 和4井寒武系丘里塔格组共取心4筒次，累计进尺17. 87米，累计取心16. 99米。 以岩心观察的方式对岩心进行详细的观察描述，并对岩心的岩性进行定名并描述岩性的变 化序列。每隔20cm对岩心进行取样，制备岩心薄片。在实验室利用光学显微镜观察薄片， 对岩性进行最终的定名，绘制岩心录井综合柱状图（如图4所示）。通过岩心、薄片观察，白 云岩均为结晶白云岩，可细分为泥晶白云岩、粉晶白云岩、细晶白云岩、中晶白云岩。
[0080] 分析要进行校正的井段的测井资料质量。根据井径测井曲线的变化来判断井壁是 否发生垮塌。当井径值大于基线值50%以上，同时其他测井曲线也在相应的深度值发生值 的突变，则认为该井段的测井资料质量低，不满足岩性识别预测的要求。对于连续扩径井 段厚度小于20m，可认为岩性在这一深度范围内没有发生明显的变化，利用相邻非扩径井段 的测井曲线值对其进行校正。当连续扩径井段厚度大于20m，应将该井段的岩性类别与相 应的测井曲线值剔除，不纳入到岩性识别样本集中，同时不再对该井段进行岩性的识别和 校正。如和4井4470m-4477m的钻头尺寸通过井径测井可知为9. 6in，其中4475. 875m至 4476. 625m井壁发生明显的垮塌，扩径超过50%，同时其他测井曲线也发生率明显的突变， 由此可认为4475. 875m至4476. 625m的测井资料质量低，需要进行校正（如图5所示）。在 操作中，需要对整个待进行录井柱状图校正的井段的测井资料都进行评价，按评价结果校 正或剔除（如图6所示）。
[0081] 和4井共进行了 20余种测井，选择对岩石物理性质敏感的5条测井曲线进行岩性 的识别和校正，包括密度测井（DEN)、声波时差测井（AC)、自然伽马测井（GR)、深侧向电阻 率测井（RD)、中子孔隙度测井（CNL)。建立岩性识别样本集，见表一，其中对于深侧向电阻 率测井值取对数。注意到岩心中中晶白云岩的厚度仅为1米，每种测井曲线仅对应了 8个 测井值。由于该种岩性可用的测井数据量太少，故可以将细晶白云岩和中晶白云岩的岩性 合并统称为细中晶白云岩。在岩性识别样本中，泥晶白云岩每条测井曲线有46个样本点， 粉晶白云岩则有42个样本点，细中晶白云岩共有32个样本点。将岩性用数字代码表不，泥 晶白云岩为1、粉晶白云岩为2、细中晶白云岩为3。
[0082] 表一
[0083]
[0084] 将岩性识别样本集加载到数据统计系统SPSS，在判别分析模块下，将岩性列定义 为分组变量，定义分组变量的最大值与最小值分别为岩性数字代码的最大值与最小值，将 测井值列定义为自变量，见表二。
[0085] 表二
[0086]
[0087] 通过Fisher判别分析，可以得到以下两个判别函数：
[0088] fl = 0? 97GR-0. 329RD+15. 467DEN+5. 587CNL-0. 102AC-36. 258
[0089] f2 = 0? 171GR+1. 107RD+2. 697DEN-2. 560CNL+0. 552AC-41. 082
[0090] 利用上述判别函数对岩性识别样本集进行岩性分类，分类正确率为92. 9%，认为 计算得到的判别函数是有效的，可以对其他井段的岩性进行识别（如图7所示）。其中泥晶 白云岩的分类正确率100%，粉晶白云岩的分类正确率为89. 7%，细中晶白云岩的分类正 确率为85. 2%，见表三。将待岩性识别地层段的测井值输入SPSS系统，利用相同的判别函 数，对岩性进行预测。最后，利用预测的岩性结果对待校正的碳酸盐岩地层录井综合柱状图 进行校正，得到校正后的录井综合柱状图，如图8所示。
[0091] 表三
[0092]
[0093] 本发明通过获取待进行录井综合柱状图校正的碳酸盐岩地层的测井资料以及取 心段碳酸盐岩地层的岩性分类，然后根据取心段碳酸盐岩地层的岩性分类和测井资料，利 用Fisher判别方法生成目标判别函数，最后利用目标判别函数，对待进行录井综合柱状图 校正的碳酸盐岩地层的未取心段的碳酸盐岩地层的岩性进行识别，并根据识别结果校正待 进行录井综合柱状图校正的碳酸盐岩地层的录井综合柱状图。其充分利用了钻井取心资料 与测井资料，建立了岩性与测井响应之间的关系，实现了基于Fisher判别分析的岩性再识 另IJ，对原有的精度较低的录井柱状图进行了校正，提高了岩性剖面精度。
[0094] 以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保 护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本 发明的保护范围之内。
【主权项】
1. 一种碳酸盐岩地层录井综合柱状图的校正方法，其特征在于，包括以下步骤： S100,获取待进行录井综合柱状图校正的碳酸盐岩地层的测井资料以及取心段碳酸盐 岩地层的岩性分类； S200,根据所述取心段碳酸盐岩地层的岩性分类和所述测井资料，利用Fisher判别方 法生成目标判别函数； S300,利用所述目标判别函数，对所述待进行录井综合柱状图校正的碳酸盐岩地层的 未取心段的碳酸盐岩地层的岩性进行识别，并根据识别结果校正所述待进行录井综合柱状 图校正的碳酸盐岩地层的录井综合柱状图。2. 根据权利要求1所述的碳酸盐岩地层录井综合柱状图的校正方法，其特征在于，所 述步骤SlOO包括以下步骤： S110,获取所述待进行录井综合柱状图校正的碳酸盐岩