seilog反演在薄储层预测中的应用

《seilog反演在薄储层预测中的应用》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、Seilog反演在缝洞型储层预测中的应用——以阆中大安寨段为例孙华（西南油气分公司勘探开发研究院川东北油气开发所）摘要：Seislog反演是以声波测井为基础，将地震记录波反演成与井旁道目的层测井声波相似的速度曲线，然后根据这种拟层速度的变化来分析岩性、储层及流体变化特征。其反演结果忠实于地震记录，具有较高的确定性，能够明显地反映岩相、岩性的空间变化，在岩性相对稳定的条件下，能较好地反映储层的物性变化。关键词：Seislog反演大安寨薄储层介屑滩介屑灰岩1引言阆中～南部区块行政区划隶属于四川省阆中市境内，勘探面积641.64km2，矿权属于中石化西南分公司，区域构造位

2、置属川中隆起北部斜坡带（图1）。区内自1979年川石30井在大安寨段获得高产油气流以来，相继在该地区完钻了38口井，获得了一批中高产油气井（如川石43井、石龙2井、石龙13井、石龙15井、石龙16井等），已累计采出原油19.3×104t、天然气6.21×108m3。大安寨段油气藏属于构造～岩性复合、薄层状、低丰度、无边（底）水、裂缝性未饱和油气藏，油气藏的重要控制因素是沉积相、储层发育程度。图1阆中-南部区块位置图2储集层特征2.1岩石学特征储层岩石类型为残余介屑灰岩、粉-粗晶介屑灰岩及亮晶介屑灰岩。残余介屑灰岩：介屑经重结晶作用后仅存残余碎片，多形成于成岩晚期。介

3、屑含量77%(均值)，泥质平均值含量1.86%，矿物成分方解石含量98%-100%，微裂缝、缝合线发育。粉－粗晶介屑灰岩：岩石基质经重结晶作用，使基质晶体增大。介屑含量75%-80%，泥质含量2%-5%，方解石含量98%-100%。亮晶介屑灰岩：在水动力较强的环境中形成的。介屑含量85%，方解石含量95%，泥质含量5%，基质为针柱状的方解石。2.2物性及孔隙结构特征岩芯样品物性分析统计，基质孔隙度0.19%～3.22%，平均孔隙度1.082%，＜2%的样品占91.34%；渗透率介于＜0.1～159.5×10-3μm2，其中：＜1.0×10-3μm2的样品占87.5%

4、，＞10×10-3μm2者仅占3.36%，属于特低孔低渗致密储层（图2）。N=543N=543图2大安寨段物性分布直方图储层基质物性致密，主要的储集空间是次生孔、洞及缝。从岩石薄片观察到得储层储集空间主要有：溶孔、粒间孔、粒内孔（图3），孔径在0.01～0.68mm之间，面孔率大都小于1%，但是裂缝、缝合线比较发育，缝宽大多在在0.005～0.7mm之间，局部缝宽可达2mm。岩芯中裂缝较发育，裂缝类型主要有水平裂缝、低斜裂缝、高斜裂缝和直立裂缝以及网状裂缝等类型（图4），以低角度缝及水平缝为主，次为高斜缝与直立缝；溶蚀孔洞在岩芯中也较发育，洞径大小不等，一般在1～5

5、mm，局部大的溶蚀洞径15mm，还有少量不规则溶蚀洞65×15mm，这些溶蚀孔洞大部分与裂缝互相连通，但也有部分溶蚀洞呈孤立存在。大型裂缝、溶洞在岩芯及薄片中均无法观察，但可从钻井放空和钻井液漏失得到证实，石龙2井2830.6～2832.4m漏失密度1.56～1.57g/cm3泥浆共444.4m3；石龙11井大一亚段2843～2847m，漏失密度1.46g/cm3泥浆共354.4m3；石龙10井千二段2751～2755.50m放空0.25m。图3石龙2井，大一，发育的成岩压裂缝和孔，图4石龙4井，大一段，井深2825.51～2825.69，为方解石充填，单偏光，×2

6、50垂直缝与低角度缝相交3Seislog反演方法原理Seilog反演是以地震资料为基础，以实际声波测井曲线做约束进行的地震资料处理技术，把测井资料纵向分辨率高和地震资料可横向密集追踪的优点有机结合起来，从而达到对储集层进行描述的目的，其方法是以声波测井曲线进行带通滤波和低通滤波后求得的地层速度作为标准，与三维地震资料合成的声波速度进行对比．通过不断地修改有关参数，使二者达到最佳吻合，从而将界面型地震剖面转换成能与钻井、测井资料直接对比的岩层型声波剖面，最终实现反演过程。地震记录的模型为：S(t)=r(t)*W(t)+n(t)其中：S(t)为地震记录，r(t)为地层反

7、射系数，W(t)为地震子波，当地震记录的信噪比较高时，可勿略N(t)。在速度和密度有差异的界面，则任何界面上的反射吸收系数R为：Ri=(ρi+1Vi+1-ρiVi)/(ρi+1Vi+1+ρiVi)(1)i=1，2，3，4，……在大多数情况下，密度的作用很小，可将其看成速度的线性函数，因此，公式（1）可合理近似为：Ri=(Vi+1-Vi)/(Vi+1+Vi)(2)将声波测井积分到以很小间隔采样的垂直刻度，则可以在每一对相邻样点之间的面上计算反射系数，从公式(2)得到：Vi+1=Vi(1+Ri)/(1-Ri)(3)通过(3)式可以实现反演过程，从地震反射数据计算合成