测井盲反褶积及在储层构型界面识别中的应用

《测井盲反褶积及在储层构型界面识别中的应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、322j【JJ西南石油大学学报(自然科学版)Vo】．32N0．22【)l0旬4JJJournalolSouthweslPetroleumIJllive~’sity(Science＆I1echnologyEdilion)Apt’．2010文章编号：1674—5086(2010)02—0063—06测井盲反褶积及在储层构型界面识别中的应用张宪国，林承焰，张涛(中国石油大学(华东)地球资源与信息学院，山东东营257061)摘要：老油田早期测井曲线分辨率低、对薄层沉积界面识别能力差，难以满足储层建筑结构研究的需要，针对这一问题，将信号反褶积原理应用于测井曲线处理。

2、通过将反褶积问题转化为期望输出与实际输出误差的极值问题构造方程组，利用测井曲线自相关间接求取线性方程组中的褶积因子自相关矩阵，求解方程组得到反褶积因子，从而实现测井曲线的盲反褶积处理该处理消除了围岩的影响，增强了测井信号频谱中高频段能量，对测井曲线的低幅度细微变化产生“放大”效应，对曲线斜率的微小变化响应敏感，提高了测井曲线对沉积界面的识别能力。通过取芯井资料研究，证明测井曲线反褶积处理结果中高频成分具有沉积学意义，将该方法应用于秘鲁D油田非取芯井储层建筑结构研究，在I)6井V层识别出5个五级界面，3个四级界面和23个三级界面。关键词：测井；盲反褶积；储层

3、建筑结构；薄层；分辨率中图分类号：TEl32；P631文献标识码：ADOI：10．3863／i．issn．1674—5086．2010．02．012目的。从信号原理的角度分析，可以将测井曲线引言抽象为仪器探测范围内地层的综合响应叠加信号，以储层构型(也称为储层建筑结构)是Miall结合自然伽马曲线为例，某一深度点的测井曲线值并不是Allen关于沉积界面的观点提出的～种有力的露头该点的真实测量响应，而是测点及其附近上下测量值地质研究方法¨，国内外沉积学者应用该方法结的加权平均值]，其测井响应可以简单表示为合油田地下资料、野外露头及现代沉积做了大量的l，=H：

4、l：X+(1)研究工作。。，建立了不同类型储层的构型模式中式，探讨了该方法应用于地下剩余油分布预测y一--测井信号值；研究的可能性引。随着油田开发程度的提高，对一某一深度点的测井参数真实值；老油田储层构型的研究成为复杂油藏剩余油挖潜研日一褶积因子，理想状况下这一褶积因子为仪究的重要手段，但由于受到围岩、测井速率等因素的器的响应函数，在实际测井中则包含了围岩影响、仪影响，测井资料垂向分辨率难以满足三级和四器参数、测井速率等因素；级层次界面识别的需要，影响了复杂储层构型研究一噪音。的开展。针对这一问题，存秘鲁D油田储层构型研由式(1)可以将测井响应的真实值表示

5、为究中，将盲反褶积方法应用于测井曲线处理，显著提x=(Y—)h(2)高了测井资料对三级和四级储层构型界面的识别能式中力，取得了良好效果。一反褶积因子。通过式(2)可以发现，只要对测井曲线进行反1测井曲线盲反褶积的原理及实现褶积处理就可以得到地层测井响应的真实值，从而消除围岩的影响，同时还能够消除测井速率、采样间反褶积是地震资料处理中常用的一种方法，它通隔等因素的影响，提高测井曲线对地层沉积界面过消除子波的影响达到提高地震资料垂向分辨率的的识别能力。收稿日期：2009—02—11基金项目：国家科技重大专项课题(2009ZX05009—003)；国家自然科学基

6、金(40872094)；山东省自然科学基金(Z2008E01)。作者简介：张宪目(1982一)，男(汉族)，山东泰安人，博士研究生，主要从事油藏描述、储层地质学及地震沉积学研究。20101．1确定性反褶积处理的原理也可以将其写成矩阵的形式从式(2)可以看出，上述反褶积处理的关键问rhh(0)rI1b(0)⋯r}Ih(m)0(0)题是求取反褶积因子，理论上讲反褶积因子是无限b(1)rhb(0)⋯rhb(m一1)a(1)长的，本研究利用最小乘方反褶积方法，设计有●；：限长度反褶积算子进行处理。rhb(m)b(m一1)⋯r(0)a(m)假设输入信号为b(t)=(

7、b(o)，6(1)，⋯，6(n))，rdb(0)反褶积因子为a(t)=(a(O)，a(1)，⋯，a(m))，则rdb(1)实际输出(9)●：()=0()6()=∑n(r)b(t—tr)=rf1h(m)(y(o)，(1)，⋯，z(M))(3)通过求解方程(9)即可求出反褶积因子a(t)。其中，M=m+n。求解方程时，考虑一种特殊情况，当d(t)=(t)时输出误差e(t)=d(t)一(t)r曲(s)=∑d(f)6(t—s)=式中，d(￡)一期望输出，且d(t)：(d(o)，d(1)，⋯，0⋯=d(M))。t=uLU．U误差能量这时方程(9)可以简化为Q：∑e(

8、)=∑(y()一d())(4)(0)(0)⋯rh(m)a(O)期望