[发明专利]盐间页岩油储层数值处理方法、装置、设备和存储介质

权利要求书

1.一种盐间页岩油储层数值处理方法，其特征在于，包括：

建立待模拟盐间页岩油储层的数值模拟模型；

获取所述数值模拟模型中的初始特性参数以及注水后每个时间步的过水倍数；

根据每个时间步的过水倍数，采用预设的储层孔渗变化数学表征模型，分别计算盐溶作用下和重结晶作用下的每个时间步的孔渗参数变化倍数；

基于所述初始特性参数、盐溶作用下和重结晶作用下的每个时间步的孔渗参数变化倍数，分别获取盐溶作用下和重结晶作用下每个时间步的孔渗参数。

2.根据权利要求1所述的方法，其特征在于，所述初始特性参数包括初始孔隙体积和初始传导率，所述孔渗参数包括孔隙体积和传导率，所述孔渗参数变化倍数包括孔隙体积变化倍数和传导率变化倍数。

3.根据权利要求2所述的方法，其特征在于，所述预设的储层孔渗变化数学表征模型包括渗透率变化数学表征模型、盐溶作用下的孔隙度变化数学表征模型和重结晶作用下的孔隙度变化数学表征模型；所述根据每个时间步的过水倍数，采用预设的储层孔渗变化数学表征模型，分别计算盐溶作用下和重结晶作用下的每个时间步的孔渗参数变化倍数，包括：

根据每个时间步的过水倍数，基于盐溶作用下的孔隙度变化数学表征模型计算盐溶作用下的孔隙度变化倍数，得到盐溶作用下的孔隙体积变化倍数；

根据每个时间步的过水倍数，基于重结晶作用下的孔隙度变化数学表征模型计算重结晶作用下的孔隙度变化倍数，得到重结晶作用下的孔隙体积变化倍数；

分别根据盐溶作用下的孔隙度变化倍数和重结晶作用下的孔隙度变化倍数，基于所述渗透率变化数学表征模型，计算盐溶作用下的渗透率变化倍数和重结晶作用下的渗透率变化倍数，得到盐溶作用下的传导率变化倍数和重结晶作用下的传导率变化倍数。

4.根据权利要求3所述的方法，其特征在于，所述盐溶作用下的孔隙度变化数学表征模型包括：

其中，φ为当前的孔隙度，无量纲；φ₀为初始的孔隙度，无量纲；为孔隙度变化倍数，无量纲；m_sample为岩样质量，单位g；a_NaCl为岩样中NaCl的质量百分数，无量纲；ρ_NaCl为NaCl的密度，单位g/L；PV_flow为过水倍数，无量纲；为岩样中NaCl全溶解所需的水孔隙倍数，无量纲；b_glauberite为岩样中钙芒硝的质量百分数，小数；ρ_glauberite为钙芒硝的密度，单位g/L；C为钙芒硝重结晶系数，无量纲；为岩样中钙芒硝全溶解所需的水孔隙倍数，无量纲；V_sample为岩样的体积，单位L。

5.根据权利要求3所述的方法，其特征在于，所述重结晶作用下的孔隙度变化数学表征模型包括：

其中，φ为当前的孔隙度，无量纲；φ₀为初始的孔隙度，无量纲；为孔隙度变化倍数，无量纲；PV_flowout为过水倍数，无量纲；α_NaCl为当前的温度、压力下的NaCl结晶析出率，单位g/L；ρ_NaCl为NaCl的密度，单位g/L；β_glauberite为当前的温度、压力下的钙芒硝结晶析出率，单位g/L；ρ_glauberite为钙芒硝的密度，单位g/L；V_sample为岩样的体积，单位L。

6.根据权利要求3所述的方法，其特征在于，所述渗透率变化数学表征模型包括：

K(φ)/K₀＝(φ/φ₀)^ckpower×[(1-φ₀)/(1-φ)]²；

其中，K(φ)为当前的渗透率，单位D；K₀为初始的渗透率，单位D；K(φ)/K₀为渗透率变化倍数，无量纲；φ为当前的孔隙度，无量纲；φ₀为初始的孔隙度，无量纲；φφ₀为孔隙度变化倍数；ckpower为指数，无量纲。