[发明专利]基于核磁视自由水孔隙度反演的致密储层高阻水层识别方法

说明书

基于核磁视自由水孔隙度反演的致密储层高阻水层识别方法。本发明利用最小二乘法迭代自动搜索最小面积差及最大相关性系数对应的转换系数计算伪毛管压力曲线，基于伪毛管压力曲线与岩心毛管压力曲线绝对误差的大小优选反演伪毛管压力曲线的方法，通过统计N个岩样的绝对误差，取绝对误差最小的W组转换系数的平均值作为区域最优转换系数计算最优的伪毛管压力曲线，再由最优的伪毛管压力曲线计算核磁孔隙度和孔隙结构参数。然后利用试气资料、最大孔喉半径、核磁孔隙度和粘土束缚水孔隙度综合分析确定毛管束缚水孔隙度和视自由水孔隙度，最终利用视自由水孔隙度、深侧向电阻率和试气资料做交会图，确定气水划分标准，完成气层和高阻水层识别。

技术领域

本发明属于地质勘探领域，具体涉及基于核磁视自由水孔隙度反演的致密储层高阻水层识别方法。

背景技术

致密储层的石英、高岭石、伊利石(水云母)等矿物具有强吸水特征，储层微观孔隙发育、孔隙结构复杂导致束缚水饱和度增高，加之物性差、储层非均质性强和气藏充注程度不足等和低地层水矿化度等因素造成储层出水，并出现“高阻水层”异常现象。导致测井交会图法、基于岩性和水性的流体性质识别法、纵横波速度比法、相关系数法等常规方法识别流体难度非常大。

发明内容

本发明要解决的技术问题是：(1)分段面积法或相似对比法计算的伪毛管压力曲线不能准确反演孔径参数，(2)常规测井解释方法不能准确识别致密储层高阻水层。

本发明提供基于核磁视自由水孔隙度反演的致密储层高阻水层识别方法,包括以下步骤，

S1基于核磁测井求解T2谱，包括，

岩石核磁共振所测得的总磁化强度信号为M(t)，t＝1，…，n，采用以下公式求解出各类孔隙T₂弛豫时间以及孔隙在总孔隙中所占的份额P_i，

…

t(j)＝j.T_E j＝1,2,3,…,n (4)

式中，P_i为第i种孔隙在总孔隙中所占的份额，T_2i为第i种孔隙的T₂弛豫时间，T_E为回波间隔，ε(j)为随机噪声，j＝1，…，n，t表示时间脉冲信号，n表示观测到的回波数，t(n)表示第n个回波的时间脉冲信号，m表示驰豫组分种类数,i表示孔隙的种类；

S2计算T2谱反向累计谱，包括，

将累计孔隙度直接转换成毛管压力曲线的进汞饱和度S_Hg，将T₂谱从大孔隙部分向小孔隙部分进行反向累加，得到压汞毛管压力曲线相似的T₂谱积分曲线，其累积饱和度为压汞毛管压力曲线的进汞饱和度，获取进汞饱和度；

S3设置大小孔径分界初始值，包括，设置大小孔径分界初始值，对不同孔径尺度进行分段拟合，大小孔径分界初始值为最大驱动压力/T₂截止值，T₂截止值通过岩心离心实验获取，从离心前的岩样T₂分布幅度之和中找到某一点，使得所述某一点左侧各点的幅度和等于离心后的T₂分布幅度和，把所述某一点作为T₂截止值T_2cuoff；

S4计算伪毛管压力曲线，包括，

已知T₂分布的条件下，根据式(5)得到伪毛细管压力曲线,基于第S3步得到的大小孔径分界初始值，分别计算大小孔径尺度分界初始值两侧的不同孔径尺度下实测压汞曲线和伪毛细管压力曲线包络面积比值,该比值分别为大、小孔径尺度部分的纵向刻度转换系数；