[发明专利]渗透率确定方法及系统

摘要

本发明提供一种渗透率确定方法及系统。该渗透率确定方法包括：创建初始微粒运移质量守恒模型和初始颗粒捕集动态模型；根据孔隙介质表面瞬间脱落的颗粒百分数及无量纲数组对初始微粒运移质量守恒模型进行优化处理，得到微粒运移质量守恒模型；根据无量纲数组对初始颗粒捕集动态模型进行优化处理，得到颗粒捕集动态模型；对微粒运移质量守恒模型进行变形，得到变形模型；对变形模型进行拉普拉斯变换处理和拉普拉斯逆变换处理，得到流体悬浮颗粒浓度；根据颗粒捕集动态模型和流体悬浮颗粒浓度得到捕集参数；根据捕集参数得到渗透率。本发明可以得到更加准确的渗透率随时间的变化趋势。

主权项

1.一种渗透率确定方法，其特征在于，包括：创建初始微粒运移质量守恒模型和初始颗粒捕集动态模型；根据孔隙介质表面瞬间脱落的颗粒百分数及无量纲数组对所述初始微粒运移质量守恒模型进行优化处理，得到微粒运移质量守恒模型；根据所述无量纲数组对所述初始颗粒捕集动态模型进行优化处理，得到颗粒捕集动态模型；对所述微粒运移质量守恒模型进行变形，得到变形模型；对所述变形模型进行拉普拉斯变换处理和拉普拉斯逆变换处理，得到流体悬浮颗粒浓度；根据所述颗粒捕集动态模型和所述流体悬浮颗粒浓度得到捕集参数；根据所述捕集参数得到渗透率；所述初始微粒运移质量守恒模型如下：所述初始颗粒捕集动态模型如下：通过如下公式根据孔隙介质表面瞬间脱落的颗粒百分数及无量纲数组对所述初始微粒运移质量守恒模型进行优化处理：t＝0时，所述无量纲数组如下：所述微粒运移质量守恒模型如下：所述颗粒捕集动态模型如下：其中，x为空间坐标，t为时间坐标，φ为孔隙度，c为流体悬浮颗粒浓度，σ_s为被孔喉捕集颗粒的浓度，σ_s(x,t)为t时刻x位置的被孔喉捕集颗粒的浓度，σ_a为壁面吸附颗粒的浓度，σ_a(x,t)为t时刻x位置的壁面吸附颗粒的浓度，U_s为颗粒运移速度，Δσ_a为高矿化度环境所对应的最大吸附浓度与低矿化度环境所对应的最大吸附浓度之差，σ_a0为原始壁面吸附颗粒的浓度，σ_a,max为壁面吸附颗粒的最大浓度，X为空间距离，T为时间，S_S(X,T)为X距离T时刻的捕集参数，L为岩心长度，U为流体运移速度，λ为渗率参数，α为颗粒运移滞后系数，f(T)如下：其中，b为颗粒脱附速率与吸附颗粒浓度的比例系数，Δσ_a为高矿化度环境所对应的最大吸附浓度与低矿化度环境所对应的最大吸附浓度之差，f₁为孔隙介质表面瞬间脱落的颗粒百分数；第一捕集参数如下：第二捕集参数如下：其中，S_s1(T)为T时刻的第一捕集参数，S_s2(T)为T时刻的第二捕集参数，X为空间距离，T为时间，α为颗粒运移滞后系数，U_s为颗粒运移速度，Λ＝λL，L为岩心长度，λ为渗率参数，b为颗粒脱附速率与吸附颗粒浓度的比例系数，Δσ_a为高矿化度环境所对应的最大吸附浓度与低矿化度环境所对应的最大吸附浓度之差，f₁为孔隙介质表面瞬间脱落的颗粒百分数，U为流体运移速度，φ为孔隙度，第三捕集参数如下：第四捕集参数如下：其中，S_s3(X,T)为X距离T时刻的第三捕集参数，S_s4(X,T)为X距离T时刻的第四捕集参数，X为空间距离，T为时间，α为颗粒运移滞后系数，U_s为颗粒运移速度，Λ＝λL，L为岩心长度，λ为渗率参数，b为颗粒脱附速率与吸附颗粒浓度的比例系数，Δσ_a为高矿化度环境所对应的最大吸附浓度与低矿化度环境所对应的最大吸附浓度之差，f₁为孔隙介质表面瞬间脱落的颗粒百分数，U为流体运移速度，φ为孔隙度，通过如下公式得到第一渗透率：其中，或，通过如下公式得到第二渗透率：其中，其中，k₁(T)为T时刻的第一渗透率，k₂(T)为T时刻的第二渗透率，k₀为初始渗透率，β为储层损害系数，S_s1(T)为T时刻的第一捕集参数，S_s2(T)为T时刻的第二捕集参数，S_s3(X,T)为X距离T时刻的第三捕集参数，S_s4(X,T)为X距离T时刻的第四捕集参数，α为颗粒运移滞后系数，w为测点距离与岩心距离之比。