面向碳酸盐岩复杂储层数字岩心的各组分提取方法
【技术领域】
[0001]本发明涉及测井技术数据处理领域，具体涉及一种面向碳酸盐岩复杂储层数字岩心的各组分提取方法。
【背景技术】
[0002]电成像测井资料纵向分辨率高，能够清晰地反映井眼所穿过地层的各种储层特征，通过研究基于成像测井资料的储层各组分提取，结合岩心铸体薄片资料，进行三维数字岩心构建，使数字岩心模型最大程度地与实际的岩石相同，为精确地预测储层的岩石物理特性与分析不同的岩石物理响应之间的内在关系提供研究基础，碳酸盐岩储层各组分(孔隙、连通部分和骨架)提取是三维数字岩心构建的基础，对三维数字岩心的建模精度影响很大。目前，主要利用电成像测井资料，采用图像分割方法提取储层各组分。其得到的结果精度较低，数字岩心模型与实际的岩石具有较大差异，因此，急待研究一种精度高，且得到的数字岩心模型最大程度地与实际的岩石相同的碳酸盐岩储层各组分提取方法。
【发明内容】
[0003]有鉴于此，有必要提供一种精度高，且得到的数字岩心模型最大程度地与实际的岩石相同的碳酸盐岩储层各组分提取方法。
[0004]—种面向碳酸盐岩复杂储层数字岩心的各组分提取方法，其包括以下步骤:
[0005]S1、加载电成像测井资料，将电成像测井资料进行预处理，得到电成像测井静态图像；
[0006]S2、采用图像修复方法生成全井眼电成像图像；
[0007]S3、将全井眼电成像图像转换为灰度图像，在灰度图像上采用双类间最大方差法进行碳酸盐岩储层的组分自动提取；
[0008]S4、基于Visual c++6.0开发平台，设置阈值调节滑动条，采用人机交互方式手动局部调节分割双阈值，对碳酸盐岩三类储层各组分交互精细提取。
[0009]本发明所述面向碳酸盐岩复杂储层数字岩心的各组分提取方法，其利用电成像测井资料，采用双类间最大方差法进行碳酸盐岩储层的组分自动提取，实现准确对储层的孔隙部分、连通部分和骨架部分进行提取。该方法具有精度高，可靠性强的特点。解决了碳酸盐岩三维数字岩心构建，使数字岩心模型最大程度地与实际的岩石相同。
【附图说明】
[0010]图1为本发明所述面向碳酸盐岩复杂储层数字岩心的各组分提取方法的流程框图；
[0011]图2为图1中步骤S3的子流程框图；
[0012]图3为图1中步骤S3的另一子流程框图；
[0013]图4为全井眼电成像图像与碳酸盐岩三类储层各组分自动提取结果图；其中，左边为原始电成像图像，右边为双类间最大方差法碳酸盐岩三类储层组分自动提取结果；
[0014]图5为碳酸盐岩三类储层各组分交互精细提取结果图；其中，左边为双阈值交互界面，右边为碳酸盐岩三类储层组分交互精细提取结果。
【具体实施方式】
[0015]为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明，应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
[0016]本发明提供一种面向碳酸盐岩复杂储层数字岩心的各组分提取方法，如图1所示，其包括以下步骤:
[0017]S1、加载电成像测井资料，将电成像测井资料进行预处理，得到电成像测井静态图像，具体的，所述预处理包括极板对齐、均衡化、产生动态与静态图像；
[0018]S2、采用图像修复方法生成全井眼电成像图像；
[0019]S3、将全井眼电成像图像转换为灰度图像，在灰度图像上采用双类间最大方差法进行碳酸盐岩储层的组分自动提取；
[0020]S4、基于Visual c++6.0开发平台，设置阈值调节滑动条，采用人机交互方式手动局部调节分割双阈值，对碳酸盐岩三类储层各组分交互精细提取。
[0021]其中，如图2所示，所述步骤S3包括以下子步骤:
[0022]S31、使用类间最大方差法，根据前景和背景的之间方差最大得到碳酸盐岩复杂储层数字岩心分类的最佳阈值，通过最佳阈值将图像分割成前景图像和背景图像两类；
[0023]S32、在分割后的两类中，再次使用类间最大方差法分别计算前景图像和背景图像的子类中的最佳分类，并求得两个最佳分类的类间最大方差；
[0024]S33、分别选取两个类间最大方差的对应阈值作为碳酸盐岩三类储层的分割阈值。
[0025]所述碳酸盐岩三类储层包括孔隙部分、连通部分和骨架部分，对碳酸盐岩储层各组分等三类分类问题，需要两个阈值进行分类，因此，提出双类间最大方差法，所述双类间最大方差法为对灰度图像进行两次最大类间方差法分隔，属于基于类间最大方差法的一种面向三类问题分割算法。
[0026]具体的，在双类间最大方差法中，设t为前景图像与背景图像的分割阈值，前景图像点数占图像比例为wo，平均灰度为UQ，背景图像点数占图像比例为W1，平均灰度为m;
[0027]则所述图像的总平均灰度值:
[0028]u=wo*uo+wi*ui (I)
[0029]所述前景图象和背景图象的方差:
[0030]g=wo* (Uo-U )*( Uo-U)+wi (u1-u )*( u1-u) (2)
[0031 ] =wo*wi*(uo-ui)*(uo-ui)
[0032]其中，当前景和背景类间方差最大时，前景图像和背景图像差异最大，所对应的灰度t是两类分类的最佳阈值。
[0033]如图3所示，所述采用双类间最大方差法进行碳酸盐岩储层的组分自动提取包括以下分步骤:
[0034]S31、根据最佳阈值t将碳酸盐岩电成像静态灰度图像P分割成pi和P2两部分；
[0035]在分割后的两部分中，对P1部分采用类间最大方差法确定阈值^，对?2部分采用类间最大方差法确定阈值t3，并根据阈值t2将Pl分割为Pll和P12，根据阈值t3将P2分割为P21和
P22；
[0036]从t2和t3两个阈值中选取两者中类间方差最大的对应阈值作为最佳阈值，由和