基于经验模态分解和能量熵判别的碳酸盐岩储集体测井识别方法
基于经验模态分解和能量熵判别的碳酸盐岩储集体测井识别方法
【技术领域】
[0001] 本发明设及一种测井方法，具体地说是设及一种基于经验模态分解和能量赌判别 的碳酸盐岩储集体测井识别方法。
【背景技术】
[0002] 碳酸盐岩油气藏在世界范围内广泛发育，油气产量已占全球油气总产量的60%。 但碳酸盐岩油气藏非均质性强，储集空间类型多样，且不同地区、不同层系差异较大，除不 整合面和古岩溶发育的溶桐带、断裂及其伴生的裂缝发育带是前人油气勘探和开发的优选 区带外，受原生、次生成因控制的微孔隙微裂缝含油也是近几年致密油研究的热点。如何准 确的识别不同尺度的有效储集体类型，是后期勘探开发的基础。目前国内外学者主要利用 常规测井资料通过多重分形聚类、测井最优化处理、自适应有效元法、曲边单位造型技术、 多重分形理论和二维有效元数值模拟、赌函数动态模糊判别等算法对碳酸盐岩的储层识别 和定量表征进行了研究，但是常规测井资料精度有限，单层计算误差较大。随着探测方法的 不断改进，利用成像测井、斯通利波测井、远探测测井、CT扫描Ξ维成像、核磁共振测井等对 储层及其流体进行可视化识别标定成为测井研究的新方向，但新方法测井毕竟成本较高， 普及较困难，实际运用较为困难。此外学者们还从缝桐储层的响应特征模拟入手对测井资 料进行标定，但数值模拟受模型条件约束较为理想化，无法满足复杂的实际地质需求。
【发明内容】
[0003] 基于上述技术问题，本发明提供一种基于经验模态分解和能量赌判别的碳酸盐岩 储集体测井识别方法。
[0004] 本发明所采用的技术解决方案是：
[0005] -种基于经验模态分解和能量赌判别的碳酸盐岩储集体测井识别方法，包括W下 步骤：
[0006] a储集体划分及基本特征
[0007] 将碳酸盐岩储集体划分为溶桐型储集体、缝桐复合型储集体和微孔隙型储集体， 研究不同储集体的测井响应特征；
[000引b储集体测井识别标准的建立
[0009] bl测井信号的经验模态分解：
[0010] bll首先将地层进行分类，确定储集层;对地层从上到下分别处理，确定需要进行 分析的井段后，找出该井段原始测井曲线s(t)的极大值和极小值点，拟合原始包络线并计 算其均值，记为mi(t)，得到第一个本征模态函数为：
[0011] hi(t) = s(t)-mi(t) (1)
[001 ^ b 12对第一个本征模态函数hi (t)重复进行上述过程k次，直到hik( t)的平均值为0; 运样，就得到了第1个本征模态函数分量ci(t)，它代表原始信号s(t)中最高频率的分量
[OOU] hi(k-i)(t)-mik(t)=hik(t) (2)
[0014] ci(t)=hik(t)
[0015] bl3将ci(t)从s(t)中分离出来，即得到一个去掉高频分量的差值信号ri(t)，可得 到残余函数
[0016] ri(t) = s(t)-ci(t) (3)
[0017] bl4将ri(t)作为原始信号，重复步骤(1)、（2)和(3)，得到第2个本征模态函数分量 C2(t)，重复η次，得到η个本征模态函数分量cn(t)，运样就有
[001 引 C4)
[0019] blf5当Cn(t)或rn(t)满足终止条件时，循环结束，由式(3)和式(4)可得到
[0020]
(5)
[0021] 其中，rn(t)为残余函数，而各个本征模态函数分量Cl(t)，C2(t)，…，Cn(t)分别包 含了信号不同时间特征尺度大小的成分；
[0022] b2测井参数选取W及能量赌的计算：
[0023] b21首先选取能有效反应储层变化的敏感测井曲线:深侧向电阻率(RD )、浅侧向电 阻率(RS)、自然伽马(GR)、密度化EN)、中子孔隙度(CNL)、声波时差(AC);对各测井信号分别 进行经验模态分解后可得到一系列本征模态函数，每个本征模态函数包括一个残余函数rn (t)和η个本征模态函数分量ci(t)，C2(t)，…，Cn(t);本征模态函数的能量Ei可由式(6)计 算，
[0024]
[0025] b22通过本征模态函数计算的能量Ei可W构成向量E= [Ei，E2，L，En]，对其进行归 一化处理可得：
[0026] E'=[pi，p2，L，pn] (7)
[0027] 其中
91 = 61作4 = 1，2，^11;该测井信号的能量赌可写为：
[0028]
(8)
[0029] b23选取测井数据进行计算，得到自然伽马(GR)、中子孔隙度(CNL)、深侧向电阻率 (RD)、浅侧向电阻率(RS)、声波时差(AC)和密度(DEN)的能量赌，分别表示为Hengr、Hencnl、 也NRD、也NRS、也NAC和也NDEN ;
[0030] c基于F i sher判别的储集体类型划分
[0031 ] 义用Fisher方法对能量赌矩阵X=[化NGR,化NCNL,化NRD,HeNRS,化NAC,化NDEN]进行判别分 析，得到每类储集体的判别公式：
[0032] y激剛綱$= 29.285 X Henrs+12.201 X He 隱巧0.288 X Hem;r+27.426 X Henden (9)
[0033] +12.309 X Hencnl+15.66 X Henac-9.829
[0038] 在实际操作过程中，对每个储层段分别处理，得到能量赌并应用式(9)至式（11)计 算到y激激廳绍日ytMi旌最大值所对应的类型就是所判别的储集体类型。
[0039] 优选的，步骤bl5中:所述终止条件是使rn(t)成为1个单调函数。
[0040] 本发明的有益技术效果是：
[0041] 与常规测井识别方法相比，本发明擬弃了单层取典型测井值的不确定性，有效提 高了解释精度，识别效果较好，整体符合率高。
【附图说明】
[0042] 图1为基于经验模态分解和能量赌判别分析的储集体类型识别流程图。
[0043] 图2a-图2f示出经验模态分解法能量赌计算结果。其中，图2a为伽马和中子交会 图；图化为深浅侧向电阻率交会图；图2c为密度和声波交会图；图2d为伽马和中子经验模态 能量赌交会图；图2e为深浅侧向电阻率经验模态能量赌交会图；图2f为密度和声波经验模 态能量赌交会图。
【具体实施方式】
[0044] 下面W塔河油田为例，对本发明作详细说明。
[0045] 塔河油田位于塔里木盆地，是我国最大的海相油气田，其寒武一奥陶系碳酸盐岩 受多期构造作用改造，具有沉积类型多样、埋藏深度大、埋藏-成岩历史漫长而复杂等特点， 且不同尺度的碳酸盐岩储集体垂向上呈多层楼房式叠置发育，储集空间差异较大，油气富 集规律不同，造成不同井组不同层位开发规律认识不清，产量变化较大，含水上升快，水淹 严重，稳产较困难。
[0046] 本发明针对塔河油田碳酸盐岩储集体类型复杂、常规测井识别难度大等问题，首 先通过岩屯、、薄片等资料将储集体进行精细刻画与分类，应用成像和阵列声波测井资料对 单井储集体进行标定，提取典型储集体的测井响应特征。引入经验模态分解算法将常规测 井资料转换为多个频带的本征模态函数集。计算每类测井参数的经验模态函数的能量赌， 研究不同储集体的能量赌特征，应用Fisher判别算法实现碳酸盐岩储集体的精确识别。
[0047] -种基于经验模态分解和能量赌判别的碳酸盐岩储集体测井识别方法，包括W下 步骤：
[004引a储集体划分依据
[0049] 将碳酸盐岩储集体划分为溶桐型储集体、缝桐复合型储集体和微孔隙型储集体， 研究不同储集体的测井响应特征。
[0050] 塔河油田孔-缝-桐等多种类型的储集体呈多成因、多期次、多层楼房式叠置发育， 且不同类型储集体侧向上相互交替，垂向上相互叠加。塔河油田常见的储集体类型主要为 溶蚀孔桐、裂缝和基质孔隙。在前人研究的基础上，通过对35口取屯、井1759.6m岩屯、观察、 2145张普通薄片、24則长显微巧光和1357张铸体薄片，并结合测井、地震和录井资料的研究， 将该区优势储集体类型依次划分为溶桐型储集体、岩溶缝桐复合型储集体和微孔隙型储集 体Ξ类。
[0051 ] b储集体测井识别标准的建立
[0052] bl测井信号的经验模态分解:经验模态分解法化MD)是一种新型自适应信号时频 处理方法，依据数据自身的时间尺度特征将复杂信号分解成一系列表征时间尺度的分量， 得到具有明显物理意义的瞬时频率。该方法与小波变换、傅里叶变换具有本质差别，在处理 非平稳随机信号上具有独特优势，已广泛应用于地震、机械、航空遥感等学科。由于不同储 集体的常规测井响应差异小，单层选取测井值时误差较大，因此常用的交会图法、判别分析 法难W精细区分储集体类型。本发明引入经验模态分解算法对测井曲线进行处理，W期得 到对储集类型敏感的信号。
[0053] bll首先找出原始信号s(t)所有的极大值和极小值点，用Ξ次样条函数拟合出原 始信号的包络线，并计算包络线的均值，记为mi(t)，那么原始信号的第1个本征模态函数 (IMF)可写为：
[0054] hi(t) = s(t)-mi(t) (1)
[0化5] bl2理论上，hi(t)是1个IMF,-般不满足IMF分量条件。为此，对hi(t)重复进行上述 过程k次，直到hi(t)符合本征模态函数的定义要求，所得到的均值趋于零为止，运样就得到 了第1个本征模态函数分量ci(t)，它代表原始信号s(t)中最高频率的分量
[0056] hi(k-i)(t)-mik(t)=hik(t) (2)
[0057] ci(t)=hik(t)
[0058] bl3将ci(t)从s(t)中分离出来，即得到一个去掉高频分量的差值信号ri(t)，即有
[0059] ri(t) = s(t)-ci(t) (3) [0060] bl4将ri(t)作为原始信号，重复步骤(1)、（2)和(3)，得到第2个本征模态函数分量 C2 (t )，重复η次，得到η个本征模态函数分量，运样就有
[0061] (4)
[0062] bl5当Cn(t)或rn(t)满足终止条件(通常使rn(t)成为1个单调函数）时，循环结束， 由式(3)和式(4)可得到
[0063]
(日)
[0064] 其中，rn(t)为残余函数，代表信号的平均趋势，而各个本征模态函数分量ci(t)，C2 (t)，…，Cn(t)分别包含了信号不同时间特征尺度大小的成分，其尺度依次由小到大。因此， 各分量也就相应地包含了从高到低不同频率段的成分，每一个频率段所包含的频率成分都 是不同的，且随信号本身的变化而变化。
[0065] 构建EMD分解流程图如图1所示，可W看出模态函数必须满足极值点数目和过0点 数目相等或相差1且在任意点处最大最小包络线的均值为0。该条件既保证了数据局部最大 值为正，最小值为负，同时也剔除了瞬时频率中波形不对称引起的震荡。
[0066] b2测井参数选取W及能量赌的计算：
[0067] b21首先选取能有效反应储层变化的敏感测井曲线:畑、RS、GR、DEN、CNL、AC;对各 测井信号s(t)分别进行上述经验模态分解后，可得到一系列的固有模态函数(IMF),且每个 IMF包括一个残余函数rnU)和η个IMF分量ci(t)，C2(t)，···，cn(t)。那么本征模态函数分量 的能量El则可W由式(6)计算得到，
[006引 Ei = J I ci(t) 12化，1 = 1,2,L,n (6)
[0069] b22由于经验模态分解法分解的正交性，因此在忽略残余量rn(t)的情况下，η个本 征模态函数的分量能量之和等于原始信号s(t)的总能量。η个本征模态函数分量分别包含 了不同的频率成分，即本征模态函数的能量向量Ε=化，E2，L，En]是单一测井信号能量在频 率域的一种自动划分。将上面的能量向量做归一化处理，可W得到：
[0070] E'=[pi，p2，L，pn] (7)
[0071] 其牛
Pi = Ei/E，i = l，2，L，n，pi反应了每一个本征模态函数在总能 量中所占的权重。
[0072] 即可W将该测井信号的模态能量赌定义为：
[0073]
(8)
[0074] b23在岩屯、刻度测井的基础上，选取工区158个典型储层段(溶桐型储集体样本数 为45，岩溶缝桐复合型储集体样本数为60，微孔隙型储集体样本数为53)的测井数据进行计 算，得到自然伽马、中子、深侧向电阻率、浅侧向电阻率、声波和密度的模态能量赌，分别表 巧为：Hengr、Henwl、Henrd、Henrs、Henac 和 Henden 〇
[0075] 图2a-图2f示出经验模态分解法能量赌计算结果。图2a、图化、图2c是Ξ种储集体 的测井响应交会图，图2d、图2e、图2f是它们的对应能量赌交会图。具体的：图2a为伽马和中 子交会图；图化为深浅侧向电阻率交会图；图2c为密度和声波交会图；图2d为伽马和中子经 验模态能量赌交会图；图2e为深浅侧向电阻率经验模态能量赌交会图；图2f为密度和声波 经验模态能量赌交会图。对比可知，不同储集体在能量赌交会图中的分布特征较为明显，说 明所提取的能量赌对储集体类型具有较好的表征能力。不同类型的储集体的能量赌分布如 表1所示。
[0076] 表 1
[0077]
[0078] C基于Fisher判别的储集体类型划分
[0079] 采用均等的先验概率，构建Ξ类储集体的判别公式为：
[0080] y激剛綱$= 29.285 X Henrs+1 2.201 X He 隱巧0.288 X Hem;r+27.426 X Henden (13) [0081 ] +12.309 X HEN(m+15.66 X Henac-9.829
[0086] 将参与建模的样本数据代入判别公式进行交叉检验回判，其符合率为98.1% ;将 未参与建模的另外55个样本点进行判别分析，仅有2个样本点判断错误，整体符合率为 96.4%;将所有样本的判别结果进行统计，整体符合率为97.7%。研究结果表明，应用能量 赌判别分析的方法来识别碳酸盐岩储集体类型效果较好。
[0087] 综上所述，本发明在前人研究的基础上，从露头和岩屯、资料出发将碳酸盐岩储集 体类型划分为溶桐型储集体、岩溶缝桐复合型储集体和微孔隙型储集体Ξ类，明确了不同 储集体的微宏观特征。用岩屯、资料、成像测井和阵列声波测井资料对单井储集体进行了标 定，定性描述了不同储集体的典型测井响应特征。提出基于测井资料经验模态分解和能量 赌判别分析的储集体自动识别方法，有效地提高了储集体识别的符合率。该方法擬弃了 W 往大都依靠单层取值的不确定性，更加细致地提取了层内测井响应的变化规律，有效地解 决了碳酸盐岩油藏开发中储集体类型难W识别的难题。
【主权项】
1. 一种基于经验模态分解和能量熵判别的碳酸盐岩储集体测井识别方法，其特征在于 包括以下步骤： a储集体划分及基本特征 将碳酸盐岩储集体划分为溶洞型储集体、缝洞复合型储集体和微孔隙型储集体，研究 不同储集体的测井响应特征； b储集体测井识别标准的建立 bl测井信号的经验模态分解： bl 1首先将地层进行分类，确定储集层;对地层从上到下分别处理，确定需要进行分析 的井段后，找出该井段原始测井曲线s(t)的极大值和极小值点，拟合原始包络线并计算其 均值，记为mi (t )，得到第一个本征模态函数为： hi(t) = s(t)-mi(t) (I) bl2对第一个本征模态函数匕^)重复进行上述过程k次，直到hlk(t)的平均值为0;这 样，就得到了第1个本征模态函数分量C1(t)，它代表原始信号s(t)中最高频率的分量 hi(k-i)(t)_mik(t) =hik(t) (2) ci(t) =hik(t) 1313将(^(〇从s(t)中分离出来，即得到一个去掉高频分量的差值信号n(t)，可得到残 余函数 ri(t) = s(t)-ci(t) (3) bl##ri(t)作为原始信号，重复步骤（I)、（2)和（3)，得到第2个本征模态函数分量C2 (t)，重复η次，得到η个本征模态函数分量cn(t)，这样就有bl5icn⑴或rn⑴满足终止条件时，循环结束，由式⑶和式⑷可得到其中，rn(t)为残余函数，而各个本征模态函数分量C1(t)，C2(t)，…，c n(t)分别包含了信 号不同时间特征尺度大小的成分； b2测井参数选取以及能量熵的计算： b21首先选取能有效反应储层变化的敏感测井曲线：深侧向电阻率、浅侧向电阻率、自 然伽马、密度、中子孔隙度、声波时差;对各测井信号分别进行经验模态分解后可得到一系 列本征模态函数，每个本征模态函数包括一个残余函数r n(t)和η个本征模态函数分量C1 (t)，C2(t)，…，cn(t);本征模态函数的能量Ei可由式(6)计算， Ei = J |ci(t) |2dt，i = l，2，L，n (6) b22通过本征模态函数计算的能量Ei可以构成向量E= [Ei,E2,L，En]，对其进行归一化处 理可得： E7 =[pi,P2,L,pn] (7) 其中i = l，2，L，n;该测井信号的能量熵可写为：b23选取测井数据进行计算，得到自然伽马、中子孔隙度、深侧向电阻率、浅侧向电阻 率、声波时差和密度的能量熵，分别表亦为Hengr、Hencnl、Henrd、Henrs、Henac和Henden; c基于Fisher 判别的储集体类型划分 米用Fisher方法对能量熵矩阵X= [Hengr，Hencnl，Henrd，Henrs，Henac ,Henden ]进行判别分析， 得到每类储集体的判别公式： y溶剛^$=29 · 285 XHenrs+12 · 201 XHenrd+20 · 288 XHengr+27.426 XHenden (9) +12.309 XHencnl+15.66 X Henac_9.829 領!顧$=77 · 226 XHenrs+30 · 819 XHem?d+36 · 891 XHengr+21 · 766 XHenden (10) +28.818 XHencnl+26.115XHenac_31.909 5W酬f顧114 · 053 X Henrs+47 · 13 X Hem?d+47 · 938 X Hengr+48 · 093 X Henden (11) +36.612 XHencnl+45.675 X Henac_72.256 在实际操作过程中，对每个储层段分别处理，得到能量熵并应用式(9)至式（11)计算到 乃容醒(_、y湖麼廳??，最大值所对应的类型就是所判别的储集体类型。2.根据权利要求1所述的一种基于经验模态分解和能量熵判别的碳酸盐岩储集体测井 识别方法，其特征在于，步骤bl5中：所述终止条件是使化(〇成为1个单调函数。
【专利摘要】本发明公开了一种基于经验模态分解和能量熵判别的碳酸盐岩储集体测井识别方法，该方法从露头和岩心资料出发将碳酸盐岩储集体类型划分为溶洞型储集体、岩溶缝洞复合型储集体和微孔隙型储集体三类，明确了不同储集体的微宏观特征。用岩心资料、成像测井和阵列声波测井资料对单井储集体进行了标定，定性描述了不同储集体的典型测井响应特征。提出基于测井资料经验模态分解和能量熵判别分析的储集体自动识别方法，有效地提高了储集体识别的符合率。该方法摒弃了以往大都依靠单层取值的不确定性，更加细致地提取了层内测井响应的变化规律，有效地解决了碳酸盐岩油藏开发中储集体类型难以识别的难题。
【IPC分类】G01V11/00
【公开号】CN105487136
【申请号】CN201511029583
【发明人】葛新民, 范宜仁, 邓少贵, 范卓颖, 邢东辉, 韩玉娇, 李润泽
【申请人】中国石油大学(华东)
【公开日】2016年4月13日
【申请日】2015年12月31日
【技术领域】
[0001] 本发明设及一种测井方法，具体地说是设及一种基于经验模态分解和能量赌判别 的碳酸盐岩储集体测井识别方法。
【背景技术】
[0002] 碳酸盐岩油气藏在世界范围内广泛发育，油气产量已占全球油气总产量的60%。 但碳酸盐岩油气藏非均质性强，储集空间类型多样，且不同地区、不同层系差异较大，除不 整合面和古岩溶发育的溶桐带、断裂及其伴生的裂缝发育带是前人油气勘探和开发的优选 区带外，受原生、次生成因控制的微孔隙微裂缝含油也是近几年致密油研究的热点。如何准 确的识别不同尺度的有效储集体类型，是后期勘探开发的基础。目前国内外学者主要利用 常规测井资料通过多重分形聚类、测井最优化处理、自适应有效元法、曲边单位造型技术、 多重分形理论和二维有效元数值模拟、赌函数动态模糊判别等算法对碳酸盐岩的储层识别 和定量表征进行了研究，但是常规测井资料精度有限，单层计算误差较大。随着探测方法的 不断改进，利用成像测井、斯通利波测井、远探测测井、CT扫描Ξ维成像、核磁共振测井等对 储层及其流体进行可视化识别标定成为测井研究的新方向，但新方法测井毕竟成本较高， 普及较困难，实际运用较为困难。此外学者们还从缝桐储层的响应特征模拟入手对测井资 料进行标定，但数值模拟受模型条件约束较为理想化，无法满足复杂的实际地质需求。
【发明内容】
[0003] 基于上述技术问题，本发明提供一种基于经验模态分解和能量赌判别的碳酸盐岩 储集体测井识别方法。
[0004] 本发明所采用的技术解决方案是：
[0005] -种基于经验模态分解和能量赌判别的碳酸盐岩储集体测井识别方法，包括W下 步骤：
[0006] a储集体划分及基本特征
[0007] 将碳酸盐岩储集体划分为溶桐型储集体、缝桐复合型储集体和微孔隙型储集体， 研究不同储集体的测井响应特征；
[000引b储集体测井识别标准的建立
[0009] bl测井信号的经验模态分解：
[0010] bll首先将地层进行分类，确定储集层;对地层从上到下分别处理，确定需要进行 分析的井段后，找出该井段原始测井曲线s(t)的极大值和极小值点，拟合原始包络线并计 算其均值，记为mi(t)，得到第一个本征模态函数为：
[0011] hi(t) = s(t)-mi(t) (1)
[001 ^ b 12对第一个本征模态函数hi (t)重复进行上述过程k次，直到hik( t)的平均值为0; 运样，就得到了第1个本征模态函数分量ci(t)，它代表原始信号s(t)中最高频率的分量
[OOU] hi(k-i)(t)-mik(t)=hik(t) (2)
[0014] ci(t)=hik(t)
[0015] bl3将ci(t)从s(t)中分离出来，即得到一个去掉高频分量的差值信号ri(t)，可得 到残余函数
[0016] ri(t) = s(t)-ci(t) (3)
[0017] bl4将ri(t)作为原始信号，重复步骤(1)、（2)和(3)，得到第2个本征模态函数分量 C2(t)，重复η次，得到η个本征模态函数分量cn(t)，运样就有
[001 引 C4)
[0019] blf5当Cn(t)或rn(t)满足终止条件时，循环结束，由式(3)和式(4)可得到
[0020]
(5)
[0021] 其中，rn(t)为残余函数，而各个本征模态函数分量Cl(t)，C2(t)，…，Cn(t)分别包 含了信号不同时间特征尺度大小的成分；
[0022] b2测井参数选取W及能量赌的计算：
[0023] b21首先选取能有效反应储层变化的敏感测井曲线:深侧向电阻率(RD )、浅侧向电 阻率(RS)、自然伽马(GR)、密度化EN)、中子孔隙度(CNL)、声波时差(AC);对各测井信号分别 进行经验模态分解后可得到一系列本征模态函数，每个本征模态函数包括一个残余函数rn (t)和η个本征模态函数分量ci(t)，C2(t)，…，Cn(t);本征模态函数的能量Ei可由式(6)计 算，
[0024]
[0025] b22通过本征模态函数计算的能量Ei可W构成向量E= [Ei，E2，L，En]，对其进行归 一化处理可得：
[0026] E'=[pi，p2，L，pn] (7)
[0027] 其中
91 = 61作4 = 1，2，^11;该测井信号的能量赌可写为：
[0028]
(8)
[0029] b23选取测井数据进行计算，得到自然伽马(GR)、中子孔隙度(CNL)、深侧向电阻率 (RD)、浅侧向电阻率(RS)、声波时差(AC)和密度(DEN)的能量赌，分别表示为Hengr、Hencnl、 也NRD、也NRS、也NAC和也NDEN ;
[0030] c基于F i sher判别的储集体类型划分
[0031 ] 义用Fisher方法对能量赌矩阵X=[化NGR,化NCNL,化NRD,HeNRS,化NAC,化NDEN]进行判别分 析，得到每类储集体的判别公式：
[0032] y激剛綱$= 29.285 X Henrs+12.201 X He 隱巧0.288 X Hem;r+27.426 X Henden (9)
[0033] +12.309 X Hencnl+15.66 X Henac-9.829
[0038] 在实际操作过程中，对每个储层段分别处理，得到能量赌并应用式(9)至式（11)计 算到y激激廳绍日ytMi旌最大值所对应的类型就是所判别的储集体类型。
[0039] 优选的，步骤bl5中:所述终止条件是使rn(t)成为1个单调函数。
[0040] 本发明的有益技术效果是：
[0041] 与常规测井识别方法相比，本发明擬弃了单层取典型测井值的不确定性，有效提 高了解释精度，识别效果较好，整体符合率高。
【附图说明】
[0042] 图1为基于经验模态分解和能量赌判别分析的储集体类型识别流程图。
[0043] 图2a-图2f示出经验模态分解法能量赌计算结果。其中，图2a为伽马和中子交会 图；图化为深浅侧向电阻率交会图；图2c为密度和声波交会图；图2d为伽马和中子经验模态 能量赌交会图；图2e为深浅侧向电阻率经验模态能量赌交会图；图2f为密度和声波经验模 态能量赌交会图。
【具体实施方式】
[0044] 下面W塔河油田为例，对本发明作详细说明。
[0045] 塔河油田位于塔里木盆地，是我国最大的海相油气田，其寒武一奥陶系碳酸盐岩 受多期构造作用改造，具有沉积类型多样、埋藏深度大、埋藏-成岩历史漫长而复杂等特点， 且不同尺度的碳酸盐岩储集体垂向上呈多层楼房式叠置发育，储集空间差异较大，油气富 集规律不同，造成不同井组不同层位开发规律认识不清，产量变化较大，含水上升快，水淹 严重，稳产较困难。
[0046] 本发明针对塔河油田碳酸盐岩储集体类型复杂、常规测井识别难度大等问题，首 先通过岩屯、、薄片等资料将储集体进行精细刻画与分类，应用成像和阵列声波测井资料对 单井储集体进行标定，提取典型储集体的测井响应特征。引入经验模态分解算法将常规测 井资料转换为多个频带的本征模态函数集。计算每类测井参数的经验模态函数的能量赌， 研究不同储集体的能量赌特征，应用Fisher判别算法实现碳酸盐岩储集体的精确识别。
[0047] -种基于经验模态分解和能量赌判别的碳酸盐岩储集体测井识别方法，包括W下 步骤：
[004引a储集体划分依据
[0049] 将碳酸盐岩储集体划分为溶桐型储集体、缝桐复合型储集体和微孔隙型储集体， 研究不同储集体的测井响应特征。
[0050] 塔河油田孔-缝-桐等多种类型的储集体呈多成因、多期次、多层楼房式叠置发育， 且不同类型储集体侧向上相互交替，垂向上相互叠加。塔河油田常见的储集体类型主要为 溶蚀孔桐、裂缝和基质孔隙。在前人研究的基础上，通过对35口取屯、井1759.6m岩屯、观察、 2145张普通薄片、24則长显微巧光和1357张铸体薄片，并结合测井、地震和录井资料的研究， 将该区优势储集体类型依次划分为溶桐型储集体、岩溶缝桐复合型储集体和微孔隙型储集 体Ξ类。
[0051 ] b储集体测井识别标准的建立
[0052] bl测井信号的经验模态分解:经验模态分解法化MD)是一种新型自适应信号时频 处理方法，依据数据自身的时间尺度特征将复杂信号分解成一系列表征时间尺度的分量， 得到具有明显物理意义的瞬时频率。该方法与小波变换、傅里叶变换具有本质差别，在处理 非平稳随机信号上具有独特优势，已广泛应用于地震、机械、航空遥感等学科。由于不同储 集体的常规测井响应差异小，单层选取测井值时误差较大，因此常用的交会图法、判别分析 法难W精细区分储集体类型。本发明引入经验模态分解算法对测井曲线进行处理，W期得 到对储集类型敏感的信号。
[0053] bll首先找出原始信号s(t)所有的极大值和极小值点，用Ξ次样条函数拟合出原 始信号的包络线，并计算包络线的均值，记为mi(t)，那么原始信号的第1个本征模态函数 (IMF)可写为：
[0054] hi(t) = s(t)-mi(t) (1)
[0化5] bl2理论上，hi(t)是1个IMF,-般不满足IMF分量条件。为此，对hi(t)重复进行上述 过程k次，直到hi(t)符合本征模态函数的定义要求，所得到的均值趋于零为止，运样就得到 了第1个本征模态函数分量ci(t)，它代表原始信号s(t)中最高频率的分量
[0056] hi(k-i)(t)-mik(t)=hik(t) (2)
[0057] ci(t)=hik(t)
[0058] bl3将ci(t)从s(t)中分离出来，即得到一个去掉高频分量的差值信号ri(t)，即有
[0059] ri(t) = s(t)-ci(t) (3) [0060] bl4将ri(t)作为原始信号，重复步骤(1)、（2)和(3)，得到第2个本征模态函数分量 C2 (t )，重复η次，得到η个本征模态函数分量，运样就有
[0061] (4)
[0062] bl5当Cn(t)或rn(t)满足终止条件(通常使rn(t)成为1个单调函数）时，循环结束， 由式(3)和式(4)可得到
[0063]
(日)
[0064] 其中，rn(t)为残余函数，代表信号的平均趋势，而各个本征模态函数分量ci(t)，C2 (t)，…，Cn(t)分别包含了信号不同时间特征尺度大小的成分，其尺度依次由小到大。因此， 各分量也就相应地包含了从高到低不同频率段的成分，每一个频率段所包含的频率成分都 是不同的，且随信号本身的变化而变化。
[0065] 构建EMD分解流程图如图1所示，可W看出模态函数必须满足极值点数目和过0点 数目相等或相差1且在任意点处最大最小包络线的均值为0。该条件既保证了数据局部最大 值为正，最小值为负，同时也剔除了瞬时频率中波形不对称引起的震荡。
[0066] b2测井参数选取W及能量赌的计算：
[0067] b21首先选取能有效反应储层变化的敏感测井曲线:畑、RS、GR、DEN、CNL、AC;对各 测井信号s(t)分别进行上述经验模态分解后，可得到一系列的固有模态函数(IMF),且每个 IMF包括一个残余函数rnU)和η个IMF分量ci(t)，C2(t)，···，cn(t)。那么本征模态函数分量 的能量El则可W由式(6)计算得到，
[006引 Ei = J I ci(t) 12化，1 = 1,2,L,n (6)
[0069] b22由于经验模态分解法分解的正交性，因此在忽略残余量rn(t)的情况下，η个本 征模态函数的分量能量之和等于原始信号s(t)的总能量。η个本征模态函数分量分别包含 了不同的频率成分，即本征模态函数的能量向量Ε=化，E2，L，En]是单一测井信号能量在频 率域的一种自动划分。将上面的能量向量做归一化处理，可W得到：
[0070] E'=[pi，p2，L，pn] (7)
[0071] 其牛
Pi = Ei/E，i = l，2，L，n，pi反应了每一个本征模态函数在总能 量中所占的权重。
[0072] 即可W将该测井信号的模态能量赌定义为：
[0073]
(8)
[0074] b23在岩屯、刻度测井的基础上，选取工区158个典型储层段(溶桐型储集体样本数 为45，岩溶缝桐复合型储集体样本数为60，微孔隙型储集体样本数为53)的测井数据进行计 算，得到自然伽马、中子、深侧向电阻率、浅侧向电阻率、声波和密度的模态能量赌，分别表 巧为：Hengr、Henwl、Henrd、Henrs、Henac 和 Henden 〇
[0075] 图2a-图2f示出经验模态分解法能量赌计算结果。图2a、图化、图2c是Ξ种储集体 的测井响应交会图，图2d、图2e、图2f是它们的对应能量赌交会图。具体的：图2a为伽马和中 子交会图；图化为深浅侧向电阻率交会图；图2c为密度和声波交会图；图2d为伽马和中子经 验模态能量赌交会图；图2e为深浅侧向电阻率经验模态能量赌交会图；图2f为密度和声波 经验模态能量赌交会图。对比可知，不同储集体在能量赌交会图中的分布特征较为明显，说 明所提取的能量赌对储集体类型具有较好的表征能力。不同类型的储集体的能量赌分布如 表1所示。
[0076] 表 1
[0077]
[0078] C基于Fisher判别的储集体类型划分
[0079] 采用均等的先验概率，构建Ξ类储集体的判别公式为：
[0080] y激剛綱$= 29.285 X Henrs+1 2.201 X He 隱巧0.288 X Hem;r+27.426 X Henden (13) [0081 ] +12.309 X HEN(m+15.66 X Henac-9.829
[0086] 将参与建模的样本数据代入判别公式进行交叉检验回判，其符合率为98.1% ;将 未参与建模的另外55个样本点进行判别分析，仅有2个样本点判断错误，整体符合率为 96.4%;将所有样本的判别结果进行统计，整体符合率为97.7%。研究结果表明，应用能量 赌判别分析的方法来识别碳酸盐岩储集体类型效果较好。
[0087] 综上所述，本发明在前人研究的基础上，从露头和岩屯、资料出发将碳酸盐岩储集 体类型划分为溶桐型储集体、岩溶缝桐复合型储集体和微孔隙型储集体Ξ类，明确了不同 储集体的微宏观特征。用岩屯、资料、成像测井和阵列声波测井资料对单井储集体进行了标 定，定性描述了不同储集体的典型测井响应特征。提出基于测井资料经验模态分解和能量 赌判别分析的储集体自动识别方法，有效地提高了储集体识别的符合率。该方法擬弃了 W 往大都依靠单层取值的不确定性，更加细致地提取了层内测井响应的变化规律，有效地解 决了碳酸盐岩油藏开发中储集体类型难W识别的难题。
【主权项】
1. 一种基于经验模态分解和能量熵判别的碳酸盐岩储集体测井识别方法，其特征在于 包括以下步骤： a储集体划分及基本特征 将碳酸盐岩储集体划分为溶洞型储集体、缝洞复合型储集体和微孔隙型储集体，研究 不同储集体的测井响应特征； b储集体测井识别标准的建立 bl测井信号的经验模态分解： bl 1首先将地层进行分类，确定储集层;对地层从上到下分别处理，确定需要进行分析 的井段后，找出该井段原始测井曲线s(t)的极大值和极小值点，拟合原始包络线并计算其 均值，记为mi (t )，得到第一个本征模态函数为： hi(t) = s(t)-mi(t) (I) bl2对第一个本征模态函数匕^)重复进行上述过程k次，直到hlk(t)的平均值为0;这 样，就得到了第1个本征模态函数分量C1(t)，它代表原始信号s(t)中最高频率的分量 hi(k-i)(t)_mik(t) =hik(t) (2) ci(t) =hik(t) 1313将(^(〇从s(t)中分离出来，即得到一个去掉高频分量的差值信号n(t)，可得到残 余函数 ri(t) = s(t)-ci(t) (3) bl##ri(t)作为原始信号，重复步骤（I)、（2)和（3)，得到第2个本征模态函数分量C2 (t)，重复η次，得到η个本征模态函数分量cn(t)，这样就有bl5icn⑴或rn⑴满足终止条件时，循环结束，由式⑶和式⑷可得到其中，rn(t)为残余函数，而各个本征模态函数分量C1(t)，C2(t)，…，c n(t)分别包含了信 号不同时间特征尺度大小的成分； b2测井参数选取以及能量熵的计算： b21首先选取能有效反应储层变化的敏感测井曲线：深侧向电阻率、浅侧向电阻率、自 然伽马、密度、中子孔隙度、声波时差;对各测井信号分别进行经验模态分解后可得到一系 列本征模态函数，每个本征模态函数包括一个残余函数r n(t)和η个本征模态函数分量C1 (t)，C2(t)，…，cn(t);本征模态函数的能量Ei可由式(6)计算， Ei = J |ci(t) |2dt，i = l，2，L，n (6) b22通过本征模态函数计算的能量Ei可以构成向量E= [Ei,E2,L，En]，对其进行归一化处 理可得： E7 =[pi,P2,L,pn] (7) 其中i = l，2，L，n;该测井信号的能量熵可写为：b23选取测井数据进行计算，得到自然伽马、中子孔隙度、深侧向电阻率、浅侧向电阻 率、声波时差和密度的能量熵，分别表亦为Hengr、Hencnl、Henrd、Henrs、Henac和Henden; c基于Fisher 判别的储集体类型划分 米用Fisher方法对能量熵矩阵X= [Hengr，Hencnl，Henrd，Henrs，Henac ,Henden ]进行判别分析， 得到每类储集体的判别公式： y溶剛^$=29 · 285 XHenrs+12 · 201 XHenrd+20 · 288 XHengr+27.426 XHenden (9) +12.309 XHencnl+15.66 X Henac_9.829 領!顧$=77 · 226 XHenrs+30 · 819 XHem?d+36 · 891 XHengr+21 · 766 XHenden (10) +28.818 XHencnl+26.115XHenac_31.909 5W酬f顧114 · 053 X Henrs+47 · 13 X Hem?d+47 · 938 X Hengr+48 · 093 X Henden (11) +36.612 XHencnl+45.675 X Henac_72.256 在实际操作过程中，对每个储层段分别处理，得到能量熵并应用式(9)至式（11)计算到 乃容醒(_、y湖麼廳??，最大值所对应的类型就是所判别的储集体类型。2.根据权利要求1所述的一种基于经验模态分解和能量熵判别的碳酸盐岩储集体测井 识别方法，其特征在于，步骤bl5中：所述终止条件是使化(〇成为1个单调函数。
【专利摘要】本发明公开了一种基于经验模态分解和能量熵判别的碳酸盐岩储集体测井识别方法，该方法从露头和岩心资料出发将碳酸盐岩储集体类型划分为溶洞型储集体、岩溶缝洞复合型储集体和微孔隙型储集体三类，明确了不同储集体的微宏观特征。用岩心资料、成像测井和阵列声波测井资料对单井储集体进行了标定，定性描述了不同储集体的典型测井响应特征。提出基于测井资料经验模态分解和能量熵判别分析的储集体自动识别方法，有效地提高了储集体识别的符合率。该方法摒弃了以往大都依靠单层取值的不确定性，更加细致地提取了层内测井响应的变化规律，有效地解决了碳酸盐岩油藏开发中储集体类型难以识别的难题。
【IPC分类】G01V11/00
【公开号】CN105487136
【申请号】CN201511029583
【发明人】葛新民, 范宜仁, 邓少贵, 范卓颖, 邢东辉, 韩玉娇, 李润泽
【申请人】中国石油大学(华东)
【公开日】2016年4月13日
【申请日】2015年12月31日
最新回复(0)