1.一种基于流压测试的多层合采气井分层地层压力计算方法,其特征在于,包括以下步骤,第一步:收集多层合采气井的井底流压测试数据、小层物性和小层产量数据;
第二步:根据流压和流压梯度测试值计算小层出气口对应的井底流压;
第三步:根据气井小层的原始含水饱和度和目前产气量、产水量,计算小层目前的含水饱和度;
第四步:根据不同泥质含量岩心的相渗实验数据,回归确定不同泥质含量岩心的相渗曲线方程;
第五步:根据小层泥质含量匹配相渗曲线方程,计算小层气相相对渗透率;
第六步:根据气井产能方程计算小层的地层压力。
2.如权利要求1所述的一种基于流压测试的多层合采气井分层地层压力计算方法,其特征在于,所述第一步包括,通过井底流压测试,获得井底流压测试值 产层中部的流压梯度值 通过测井解释,获得小层的渗透率、有效厚度、原始含水饱和度、泥质含量等数据;通过生产数据现场录取,获得测试时期的气井产气量、产水量。
3.如权利要求1所述的一种基于流压测试的多层合采气井分层地层压力计算方法,其特征在于,所述第二步包括,根据以下公式计算小层出气口对应的井底流压:式中, 表示气井井底流压,单位是MPa;
表示井底流压梯度,单位是MPa/m;
表示流压测试点对应的深度,单位是m;
Hi表示第i层的深度,单位是m;
i表示多层合采小层的序号,i=1…N。
4.如权利要求1所述的一种基于流压测试的多层合采气井分层地层压力计算方法,其特征在于,所述第三步包括,根据以下公式计算目前小层的含水饱和度:式中Sw表示任意时刻的含水饱和度,无因次;
Swr表示束缚水饱和度,无因次;
qw表示日产水量,单位是m3/d;
Bw表示地层水体积系数,无因次;
Bg表示天然气体积系数,无因次;
i表示多层合采小层的序号,i=1…N。
5.如权利要求1所述的一种基于流压测试的多层合采气井分层地层压力计算方法,其特征在于,所述第四步包括,收集不同泥质含量岩心的相渗实验数据,利用指数相渗曲线方程拟合实验数据,确定以下不同泥质含量小层的相渗曲线方程,式中,a、b、m、n表示待定拟合系数,无因次,由指数式相渗曲线方程回归拟合得到。
Sw表示目前小层的含水饱和度,无因次。
6.如权利要求1所述的一种基于流压测试的多层合采气井分层地层压力计算方法,其特征在于,所述第五步包括,根据小层泥质含量多少,匹配泥质含量最接近的相渗曲线方程,并将第三步计算的目前含水饱和度代入第四步获得的小层对应的相渗曲线方程中计算小层气相渗透率Krgi。
7.如权利要求1所述的一种基于流压测试的多层合采气井分层地层压力计算方法,其特征在于,所述第五步包括,根据以下公式分别计算小层的地层压力,式中pr表示小层地层压力,单位是MPa;为待求值;
pwf表示井底流压,单位是MPa;在第二步中计算得出;
K表示渗透率,单位是mD;测井解释得到;
Krg表示小层气相相对渗透率,无因次;在第五步中计算得出;
h表示小层厚度,单位是m;测井解释得到;
qsc表示小层产量,单位是m3/d;产量劈分或者产气剖面测试得到;
re表示泄流半径,单位是m;由不稳定典型曲线拟合得到;
rw表示井半径,单位是m;完井参数得到;
S表示表皮系数,无因次;试井解释得到;
T表示地层温度,单位是K;由测试得到μ表示气体粘度,mPa.s;高压物性分析得到;
Z表示气体偏差系数,无因次;高压物性分析得到;
下标i表示多层合采小层的序号,i=1…N。