[发明专利]一种核磁共振测井资料定量计算储层产液剖面方法

权利要求书

1.一种核磁共振测井资料定量计算储层产液剖面方法，其特征在于，包括以下步骤：

(1)利用已知的岩心核磁资料得到T2截止值、束缚水饱和度、T2几何平均值，建立本区域T2截止值与T2几何平均值的模型；

(2)利用已知的核磁共振测井资料得到储层逐深度点的T2谱及几何平均值，结合步骤(1)建立的模型计算不同深度的变T2截止值以及束缚水饱和度，所述的变T2截止值束缚水饱和度是将实际核磁共振测井资料的T2几何平均值代入岩心建立的T2截止值与T2几何平均值模型得到逐深度点的变T2截止值，然后得到束缚水饱和度，束缚水饱和度公式如(1)式所表示：

式中BFV为束缚水孔隙，单位为％；φ_t为总孔隙，单位为％；S_wirr为束缚水饱和度，单位为小数；

(3)利用轻烃的密度与核磁共振测井计算含烃校正总孔隙度和冲洗带含烃饱和度，后者可被定为残余烃饱和度，所述的密度与核磁共振测井资料计算冲洗带含烃饱和度是气或者轻烃的密度、核磁两个不同响应方程计算的结果，如公式(2)、(3)所示：

ρ_b＝ρ_ma(1-φ_t)+ρ_fφ_t(1-S_gxo)+ρ_gφ_tS_gxo (2)

TCMR＝φ_tS_gxo(HI)_gP_g+φ_t(1-S_gxo)(HI)_f (3)

式中ρ_b是指测井密度，单位为g/cm³；ρ_ma是骨架密度，单位为g/cm³；ρ_f是流体密度，单位为g/cm³；TCMR是核磁共振总孔隙度，单位为％；f是地层水缩写，无量纲；g是气体缩写，无量纲；

式中S_gxo为地层中的含烃饱和度，单位为小数；ρ_g是烃密度，单位为g/cm³；P_g是极化因子，无量纲；HI是含氢指数，无量纲；

(4)利用随钻深电阻率测井数据可计算储层原始含水饱和度，所述的原始含水饱和度是阿尔奇公式计算的原状地层总含水饱和度，公式如(4)所示：

式中a、b、m、n为阿尔奇指数，无量纲；R_w是地层水电阻率，单位Ωm；R_t是原状地层电阻率，单位Ωm；S_w是原状地层总含水饱和度，单位％；

(5)利用已知的岩心毛管压力与核磁T2谱建立非线性幂函数拟合关系，应用到实际核磁共振测井T2谱便可转换得到逐深度点的毛管压力，所述的岩心毛管压力与核磁T2谱幂函数拟合关系是指：使T2谱的总幅度之和为100％，将T2谱从大孔向小孔进行反向累加，得到一条在物理意义和形态上与压汞毛管压力曲线相似的T2谱累积曲线，不再用简单线性公式进行拟合而是在对数坐标下用非线性幂函数进行拟合，如式(5)所示：

其中P_c为压汞过程中施加的毛管压力，单位为MPa；T2为核磁共振测井横向弛豫时间，单位为ms；m、n均为岩心分析参数，无量纲；

(6)建立Burdine相对渗透率方程，输入步骤(2)中的束缚水饱和度，输入步骤(3)中的残余烃饱和度，输入步骤(5)中逐深度点的毛管压力，得到逐深度点的相对渗透率曲线，所述的Burdine相对渗透率是利用毛管压力曲线、束缚水、残余烃饱和度计算相渗曲线的一种方法，其为公式(6)、(7)所示：

其中，ds代表饱和度Sw的积分方式，K_rwt、K_rnwt分别是湿相与非湿相的相渗曲线，单位小数；τ_rwt为孔隙介质湿相迂曲度，τ_nrwt为孔隙介质非湿相迂曲度，单位无量纲；

(7)根据步骤(6)中得到的相对渗透率曲线，结合PVT样品的流体粘度，利用产水率方程计算产水率，所述产水率方程为公式(8)所示：

其中，F_w是任何已知饱和度时的分流量，单位小数， 是该饱和度时水相和烃相相对渗透率的函数，K_o是烃的相对渗透率，K_w是水的相对渗透率，单位小数；μ₀是烃的粘度，μ_w是水的粘度，单位mPa.s；

(8)根据步骤(7)中的产水率，把它作为求积函数，在各深度把储层渗透率下限以上的渗透率用产水率刻度，然后把这些相对的水、烃贡献在测试层段深度范围内累加得到l00％，即可得到流量剖面，所述的储层渗透率下限为实际试油或者测压取样资料确定的储层产出渗透率下限，被积分函数如(9)式所示：

其中，btm是指深度的底深、top是指深度的顶深；F_wi是某深度点的分流量，单位小数；K_i是某深度点渗透率值，单位μm²。

2.根据权利要求1所述的一种核磁共振测井资料定量计算储层产液剖面方法，其特征在于：步骤(1)所述的T2截止值、束缚水饱和度、T2几何平均值，分别是岩心核磁饱和水T2谱与离心T2谱确定的T2截止值、岩心T2截止值计算岩心束缚水饱和度、饱和水岩心核磁T2谱的几何均值。