1.一种利用液面测试数据进行井底压力计算的方法,其特征在于,包括:获取井筒液面测试数据;将井筒沿着管柱液面到油层中深分隔成多个微元段;确定各个微元段的井斜角;根据所述井斜角确定各个微元段的温度;根据所述温度确定各个微元段的地层水的物性参数;根据所述温度确定各个微元段的油的物性参数;根据所述地层水的物性参数和所述油的物性参数确定各个微元段油水混合物的密度;根据所述油水混合物的密度以及所述井斜角计算井底压力。
2.如权利要求1所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述液面测试数据包括:井筒液面高度、井斜角、套压、油层段顶深、油层段底深、地表温度、温度梯度、地层水矿化度、地层水类型、天然气相对密度、原油相对密度、原油体积系数、质量含水率。
3.如权利要求2所述的利用液面测试数据进行井底压力计算的方法,其特征在于,在所述将井筒沿着管柱液面到油层中深分隔成多个微元段步骤中,包括:计算油层中深;其计算公式为:计算微元段长度;其计算公式为:上式中,H u 为油层顶深,H d 为油层底深,H为油层中深,h 0 为实测液面高度,h i 为沿着管柱液面至油层中深间各个微元段长度,单位均为米;n为微元段个数。
4.如权利要求3所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述根据所述井斜角确定各个微元段的温度包括:计算液面位置温度;其计算公式为:T 0 =T d +T t h 0 cosθ 0 /100计算各个微元段的温度:其计算公式为:上式中,T i 为各个微元段的温度,T 0 为液面位置温度,T d 为地表温度,单位均为摄氏度;T t 为地温梯度,摄氏度/米;θ i 为各个微元段的井斜角,度。
5.如权利要求4所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述根据所述温度确定各个微元段的地层水的物性参数包括:计算各个微元段纯水的密度;计算各个微元段地层水的密度;计算各个微元段纯水的体积系数;计算各个微元段地层水的体积系数。
6.如权利要求5所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段纯水的密度的公式为:ρ iw ′=0.9998+6.4993×10 -5 T i -7.9987×10 -6 T i 2 +4.451×10 -8 T i 3 -1.2357×10 -10 T i 4上式中,ρ iw ′为各个微元段纯水密度,克/立方厘米。
7.如权利要求6所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段地层水的密度的公式为:上式中,m′为地层水的矿化度,毫克/升;ρ iwt 为各个微元段地层温度下地层水的密度,克/立方厘米。
8.如权利要求7所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段纯水的体积系数的公式为:B iw =A 1 +A 2 (145.03P i )+A 3 (145.03P i ) 2A j =a 1 +a 2 M i +a 3 M i 2 j=1,2,3M i =1.8T i +32上式中,B iw 为各个微元段纯水的体积系数,无因次;P i 为各个微元段压力,兆帕。
9.如权利要求8所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段地层水的体积系数的公式为:B iwb =B iw {a 1 (145.03P i )+[a 2 +a 3 (145.03P i )](M i -60)+[a 4 +a 5 (145.03P i )](M i -60) 2 }S i +B iw a<sub>1</sub> a<sub>2</sub> a<sub>3</sub> a<sub>4</sub> a<sub>5</sub> 5.1×10<sup>-8</sup> 5.47×10<sup>-6</sup> -1.95×10<sup>-10</sup> -3.23×10<sup>-8</sup> 8.5×10<sup>-13</sup>上式中,B iwb 为各个微元段地层水的体积系数,无因次。
10.如权利要求9所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述根据所述温度确定各个微元段的油的物性参数包括:计算各个微元段的原油溶解气油比;计算各个微元段的原油泡点压力;计算各个微元段的原油密度;计算各个微元段原油的体积系数。
11.如权利要求10所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段的原油溶解气油比的公式为: 系数 r<sub>o</sub>≥0.8762 r<sub>o</sub><0.8762 c<sub>1</sub> 0.0362 0.0178 c<sub>2</sub> 1.0937 1.187 c<sub>3</sub> 25.724 23.931上式中,R is 为各个微元段原油溶解气油比,立方米/立方米;r g 为天然气的相对密度,无因次;r o 为原油的相对密度,无因次。
12.如权利要求11所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段的原油泡点压力的公式为:α i =0.00091M i -0.0125D上式中,P ib 为各个微元段原油泡点压力,兆帕。
13.如权利要求12所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段的原油密度包括:当P i >P ib 时,各个微元段的原油为饱和原油;各个微元段饱和原油密度为:上式中,ρ iob 为各个微元段饱和原油的密度,千克/立方米。
14.如权利要求12所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段的原油密度包括:当P i ≤P ib 时,各个微元段的原油为不饱和原油;各个微元段不饱和原油密度为:F i =10 -5 (-1433+28.075+17.2M i -1180r g +12.61D)上式中,ρ io 为各个微元段不饱和原油的密度,千克/立方米。
15.如权利要求13所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段原油的体积系数的公式为:B io =B iob exp[-c io (P i -P ib )]a 1 =-2540.8 a 2 =28.07 a 3 =30.96 a 4 =-1180.0 a 5 =1784.3 a 6 =100000B iob =0.972+1.1213×10 -2 F i ′ 1.175上式中,B io 为各个微元段原油的体积系数,无因次。
16.如权利要求14所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段原油的体积系数的公式为: 系数 r<sub>o</sub>≥0.8762 r<sub>o</sub><0.8762 e<sub>1</sub> 2.620×10<sup>-3</sup> 2.622×10<sup>-3</sup> e<sub>2</sub> 1.751×10<sup>-5</sup> 1.100×10<sup>-5</sup> e<sub>3</sub> -1.062×10<sup>-7</sup> 7.507×10<sup>-9</sup>上式中,B io 为各个微元段原油的体积系数,无因次。
17.如权利要求15或16所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述根据所述地层水的物性参数和所述油的物性参数确定各个微元段油水混合物的密度包括:计算地面原油体积含水率;计算各个微元段原油体积含水率;计算各个微元段油水混合物的密度。
18.如权利要求17所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算地面原油体积含水率的公式为:上式中,f w ′为地面原油体积含水率,无因次;f w 为质量含水率,无因次。
19.如权利要求18所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述计算各个微元段原油体积含水率的公式为:上式中,f iw ′为各个微元段体积含水率。
20.如权利要求19所述的利用液面测试数据进行井底压力计算的方法,其特征在于,计算各个微元段油水混合物的密度的公式为:ρ ic =ρ io (1-f iw ′)+ρ iwt f iw ′上式中,f iw ′为各个微元段油水混合物的密度,千克/立方米。
21.如权利要求20所述的利用液面测试数据进行井底压力计算的方法,其特征在于,所述根据所述油水混合物的密度以及所述井斜角计算井底压力的公式为:P H =P n =P t +ρ c1 gh 1 cosθ 1 +ρ c2 gh 2 cosθ 2 +…+ρ cn gh n cosθ n上式中,P H 为井底压力,兆帕;P t 为套压,兆帕。
22.一种利用液面测试数据进行井底压力计算的装置,其特征在于,包括:获取数据模块,用于获取井筒液面测试数据;井筒分割模块,用于将井筒沿着管柱液面到油层中深分隔成多个微元段;井斜角确定模块,用于确定各个微元段的井斜角;温度计算模块,用于根据所述井斜角确定各个微元段的温度;第一计算模块,用于根据所述温度确定各个微元段的地层水的物性参数;第二计算模块,用于根据所述温度确定各个微元段的油的物性参数;第三计算模块,用于根据所述地层水的物性参数和所述油的物性参数确定各个微元段油水混合物的密度;井底压力计算模块,用于根据所述油水混合物的密度以及所述井斜角计算井底压力。