[发明专利]一种确定地层中洞体积的方法

权利要求书

1.一种确定地层中洞体积的方法，其特征在于，所述方法包括：

步骤一、建立缝洞型油藏试井模型，并将地层中的溶洞视作点源，所述缝洞型油藏试井模型中包括至少两个溶洞，第一溶洞与井筒连接；

步骤二、根据所述缝洞型油藏试井模型，采用叠加原理获得地层中多个溶洞同时存在时的真实空间井底压力解，所述步骤二包括：

步骤a、对所述缝洞型油藏试井模型进行拉普拉斯变换，得到第一拉普拉斯空间井底压力解函数；

步骤b、采用叠加原理，确定源汇在井筒处的无因次压力函数；

步骤c、根据所述源汇在井筒处的无因次压力函数进行拉普拉斯变换，得到第二拉普拉斯空间井底压力解函数；

步骤d、根据所述第一拉普拉斯空间井底压力解函数和第二拉普拉斯空间井底压力解函数确定考虑点源函数后的拉普拉斯空间上的井底压力解函数；

步骤e、根据所述考虑点源函数后的拉普拉斯空间上的井底压力解函数反演得到地层中多个溶洞同时存在时的真实空间井底压力解；

步骤三、利用多个溶洞同时存在时的真实空间井底压力解与实测的井底压力数据进行拟合，根据拟合结果确定出第二溶洞到第一溶洞的无因次距离；

步骤四、根据所述第二溶洞到第一溶洞的无因次距离确定第二溶洞的无因次储集系数，并根据所述第二溶洞的无因次储集系数确定所述第二溶洞的体积，其中，根据如下表达式计算所述第二溶洞的体积：

其中，V₂表示第二溶洞的体积，φ表示外部介质孔隙度，C_t表示外部介质压缩系数，h₁表示井筒相连地层的厚度，h₂表示溶洞相连地层的厚度，r_w表示井筒半径，C_2vD表示第二溶洞的无因次储集系数，C_f表示流体压缩系数，其中包括：

步骤f、根据所述第二溶洞到第一溶洞的无因次距离和无因次生产时间确定所述第一溶洞在所述第二溶洞到第一溶洞的无因次距离处的无因次地层压力，在所述步骤f中，先利用洞的渗流方程，确定拉普拉斯空间下存在第一溶洞时的地层压力分布，再根据拉普拉斯空间下存在第一溶洞时的地层压力分布反演确定所述第一溶洞在所述第二溶洞到第一溶洞的无因次距离处的无因次地层压力；

步骤g、根据所述无因次生产时间和第一溶洞在所述第二溶洞到第一溶洞的无因次距离处的无因次地层压力计算所述第二溶洞的无因次储集系数，其中，根据如下表达式计算所述第二溶洞的无因次储集系数：

其中，C_2vD表示第二溶洞的无因次储集系数，P_D表示无因次地层压力，L_D表示第二溶洞到第一溶洞的无因次距离，t_pD表示无因次生产时间。

2.如权利要求1所述的方法，其特征在于，在所述步骤a中，所述第一拉普拉斯空间井底压力解函数表示为：

其中，表示第一拉普拉斯空间井底压力解函数，u表示拉普拉斯变量，M₁表示中间参数，B表示体积系数。

3.如权利要求2所述的方法，其特征在于，根据如下表达式计算所述体积系数B：

其中，M₁、M₂、M₃、M₄、T表示中间参数，K₀表示零阶第二类贝塞尔函数，s_w表示井筒的表皮系数，s_v表示溶洞的表皮系数，r_vD表示无因次溶洞半径，K₁表示一阶第二类贝塞尔函数，C_wD表示无因次井筒储集系数，λ表示无因次高度，C_vD表示无因次溶洞储集系数，α表示修正系数，β表示波动系数，γ表示阻尼系数。