CN105178943A - 一种实时校正井筒压力的方法 - Google Patents

Abstract

本发明提供了一种实时校正井筒压力的方法，属于石油钻井工程技术领域。所述方法包括：根据钻井井筒压力基本参数和控压钻井边界条件计算获得待测井的校前环空压耗值，并通过计算获得待测井的静液柱压力值；通过压力随钻监测系统PWD获取待测井的实测井底压力值，并基于压力守恒原理对井筒压力中的任意一种进行实时校正，所述井筒压力包括井底压力、环空压耗和井口回压。本发明提供的方法基于压力守恒原理对计算获得和通过PWD采集的井筒压力数据对井筒压力进行实时校正，具有校正井筒压力简便快速且精确的较高的特点，并且不需要复杂的网格划分及控制方程组的离散求解，对计算机性能要求较低，更适用于钻井施工的现场工况。

Description

一种实时校正井筒压力的方法

技术领域

本发明涉及一种利用井下压力测量工具实时校正井筒压力的方法，属于石油钻井工程技术领域。

背景技术

随着石油勘探向深部复杂地层发展，窄密度窗口、喷漏同层、喷漏同存等井下复杂情况已成为影响和制约石油钻井施工的技术瓶颈。窄密度安全窗口地层孔隙压力和破裂压力、漏失压力之间的安全钻井窗口较小，在常规钻井方式下环空压力波动范围将超过此安全钻井窗口，易造成井漏、井涌、井壁失稳、卡钻等复杂事故，导致生产时效低、钻井周期长、成本高等问题。因此，对石油、天然气钻井井筒动态压力的精确、快速计算和精细控制提出了更高的要求。

目前，现有的石油钻井井筒气液两相流模拟计算方法主要基于连续气柱理论和气液两相流理论。气液两相流理论通过划分不同的流型，建立气液两相连续性方程、动量方程来模拟流动状态，不同的计算方法误差均较大，精度难以满足压力敏感地层精细控压钻井井筒动态压力的精确、快速计算需求。

具体的，在地层压力敏感的窄密度窗口地层，井底压力如果控制不当就会形成溢流、井口等复杂事故。控压钻井能够很好的控制井筒压力，减少事故的发生，而控压钻井气液两相流动井筒动态压力精确快速计算是控压钻井实现更加精确的井筒压力控制问题的核心前提。现有技术提出了一种通过建立井筒环空内的多相流动控制方程组的方法，将多相流动计算的时间和空间域进行网格划分，然后计算得到压力参数，但是该方法不仅基础参数获取麻烦且计算较为繁琐，计算结构的误差也较大。

发明内容

本发明为解决现有的井筒压力计算方法存在的计算过程较复杂且结果误差较大、计算精度难以满足压力敏感地层精细控压钻井井筒压力计算需求的问题，进而提出了一种利用井下压力测量工具实时校正井筒压力的方法，具体包括如下的技术方案：

一种实时校正井筒压力的方法，包括：

根据钻井井筒压力基本参数和控压钻井边界条件计算获得待测井的校前环空压耗值，并通过计算获得待测井的静液柱压力值；

通过压力随钻监测系统PWD获取待测井的实测井底压力值，并基于压力守恒原理对井筒压力中的任意一种进行实时校正，所述井筒压力包括井底压力、环空压耗和井口回压。

在本发明所述的实时校正井筒压力的方法中，所述方法还包括：

若无法通过PWD获取待测井的实测井底压力值，则通过实时校正后的井底压力对待测井的环空压耗和井口回压进行校正。

在本发明所述的实时校正井筒压力的方法中，所述井底压力值通过以下公式计算获得：

P_b(t)＝P_h(t)+P_f(t)+P_w(t)

其中，P_b(t)表示t时刻的井底压力；P_h(t)表示t时刻的静液柱压力；P_f(t)表示t时刻的环空压耗；P_w(t)表示t时刻的井口回压。

在本发明所述的实时校正井筒压力的方法中，在对所述井筒压力中的任意一种进行实时校正的过程中，对所述环压空耗的校正包括：

将环空摩阻校核系数带入所述环空压耗的计算公式，对所述环压空耗进行校正，所述环空摩阻校核系数通过以下公式计算获得：

ΔP(t)＝P_b(t)-P_pwd(t)

P`_f＝P_f(t)-ΔP(t)

$f^{`}=\frac{{p^{`}}_f}{p_f}f$

其中，P_pwd(t)表示PWD在t时刻获取的实测井底压力值；ΔP(t)表示计算获得的井底压力与PWD在t时刻获取的实测井底压力值之间的差值；f`表示环空摩阻校核系数，f表示环空摩阻系数。

在本发明所述的实时校正井筒压力的方法中，所述环空压耗通过以下公式计算获得：

$P_f\left(t\right)=f\frac{{\mathrm{\ρ}}_{mix}{\mathrm{Lv}}_{mix}^2\left(t\right)}{2D_{\mathrm{\α}}}$

$v_{mix}\left(t\right)=\frac{Q_{mix}\left(t\right)}{A}$

$D_{\mathrm{\α}}=\frac{D_o-D_i}{2}$

其中，ρ_mix表示井筒内混合流体的密度，D_α表示水力直径，D₀表示井筒直径，D_i表示钻具外径；Q_mix(t)表示质量流量计的测量值，A表示环空过流面积，L表示井筒环空长度，v_mix表示混合物流速。

在本发明所述的实时校正井筒压力的方法中，所述环空摩阻系数通过以下公式计算获得：

$\frac{I}{\sqrt{f}}=-1.08{\log }_{10}\[\frac{6.9}{\mathrm{Re}}+K{\left(\frac{\mathrm{\ϵ}/D_{\mathrm{\α}}}{3.7}\right)}^{1.11}\]$

$\mathrm{Re}=\frac{{\mathrm{\ρ}}_{mix}{v}_{mix}\left(t\right)D_{\mathrm{\α}}}{\mathrm{\μ}}$

其中，ε/D_α表示相对粗糙度，Re表示雷诺数。

在本发明所述的实时校正井筒压力的方法中，在对所述井筒压力中的任意一种进行实时校正的过程中，对所述井底压力的校正包括：

将t+1时刻的井底压力带入所述井底压力的计算公式，对所述井底压力进行校正，所述t+1时刻的井底压力通过以下公式计算获得：

$\mathrm{\α}=\frac{P_{pwd}\left(t\right)}{P_b\left(t\right)}$

P`_b＝αP_b(t)

其中，α表示PWD在t时刻获取的实测井底压力值与t时刻的井底压力之间的比值；P`_b表示t+1时刻的井底压力。

在本发明所述的实时校正井筒压力的方法中，所述井口回压通过以下公式计算获得：

P_w(t)＝P_w0-ΔP_h(t)+ΔP_safe

${\mathrm{\ΔP}}_h\left(t\right)=-\frac{({\mathrm{\ρ}}_l-{\mathrm{\ρ}}_g)V_g}{v}gH$

其中，P_w0表示未发生溢流时的井口回压值，ΔP_h(t)表示静液柱压力的变化值，ΔP_safe表示安全附加压力值，ρ_l表示环空液体密度，ρ_g表示平均高压力下的气体密度，V_g表示井下溢流量，g表示重力加速度，H表示待测井的井深。

在本发明所述的实时校正井筒压力的方法中，所述静液柱压力通过以下公式计算获得：

P_h(t)＝ρ_mixgH(t)

${\mathrm{\ρ}}_{mix}=\frac{m_g+m_l}{V\left(t\right)}$

m_g＝ρ_gV_g

其中，ρ_mix表示井筒内混合流体的密度；m_g表示井筒环空中气体质量且，ρ_g表示平均高压力下的气体密度，V_g表示井下溢流量，m_l表示环空液体质量；V(t)表示t时刻的环空体积并由井身结构、裸眼段直径及入井钻柱体积求得。

一种井筒压力测量装置，包括：

PWD，用于测量井底压力；

压力测量计，用于测量井口回压；

质量流量计，用于测量井筒内混合流体的质量；

液位计，用于测量井下溢流量。

本发明的有益效果是：基于压力守恒原理对计算获得和通过PWD采集的井筒压力数据对井筒压力进行实时校正，具有校正井筒压力简便快速且精确的较高的特点，并且不需要复杂的网格划分及控制方程组的离散求解，对计算机性能要求较低，更适用于钻井施工的现场工况。

附图说明

图1是以示例的方式示出了实时校正井筒压力的方法的流程图。

图2是以示例的方式示出了井筒压力测量装置的结构图。

图3为一优选实施例提供的井筒压力计算方法流程图。

具体实施方式

在本领域的现有技术中，由于多相流的压力波和密度波的不稳定性造成其压力传播更加复杂，因此不易进行控压钻井。而在控压钻井作业过程中，目前只能处理含气量不超过15％的多相流体控压钻井安全钻进，这也是精细控压钻井系统采用的自动节流阀节流进行压力控制技术本质决定的：另外，节流阀不能在气相含量过高的条件下保证节流压力的稳定性。因此可认为控压钻井的目的主要是进行微过平衡的压力控制方式钻井，所以本发明根据以井筒气液两相流动理论，由质量及压力守恒基本原理，建立适用于控压钻井的井筒压力水力计算模拟的快速、简便方法，并采用PWD(PressureWhileDrilling，压力随钻监测系统)采集的环空压力数据对水力计算模型进行实时校正，从而大幅度优化提高钻井井筒动态压力流动方程的计算精度，优化后的水力计算模型可用于各种工况下的控压钻井井筒动态的水力参数实时计算。

其中，利用PWD实时校正的井筒压力精确、快速计算方法是由质量及压力守恒基本原理，依据井筒内及钻井液循环罐总流体质量变化与井筒内各压力变化的关系，建立钻井井筒动态压力流动方程，其核心两个基础公式为：井筒环空体积+泥浆罐体积＝常量，井底压力＝静液柱压力+环空摩阻+井口回压。

下面通过具体的实施例对本发明提出的实时校正井筒压力的方法进行详细说明，结合图1所示，所述方法包括：

步骤11，根据钻井井筒压力基本参数和控压钻井边界条件计算获得待测井的校前环空压耗值，并通过计算获得待测井的静液柱压力值。

由于地层出水或液体的情况与循环钻井液密度相差不大，使井筒压力的压力变化较为缓慢，控压钻井相对较为容易，所以本实施例不考虑出水或其它流体的情况，仅考虑地层是在出气的情况下进行控压钻井的井筒压力计算。本实施例以井筒气液两相流动理论为基础，由质量及压力守恒原理，建立适用于控压钻井的钻井井筒动态压力流动方程。

图2所示的是井筒压力测量装置，包括有PWD21、压力测量计22、质量流量计23和液位计24。首先，泥浆泵25将井筒内混合液体抽送至PWD21，由PWD21侧量t时刻的井底压力，位于钻井液地面的压力测量计22用于测量t时刻的井口回压P_w(t)，质量流量计23用于测量t时刻的井筒内混合流体的质量Q_mix(t)，超声波液位计24用于测量t时刻的井下溢流量V_g(t)。

由于总的钻井液体积＝钻具水眼体积+井筒环空体积+泥浆罐体积＝常量，并且钻具在一定时间段内可认为是不变的，则钻具水眼体积相对不变，即可认为：井筒环空体积+泥浆罐体积＝常量。

进一步地，不考虑由于流体加速运动，则井底压力可以表示为：

P_b(t)＝P_h(t)+P_f(t)+P_w(t)(1)

式(1)中的P_b(t)表示t时刻的井底压力；P_h(t)表示t时刻的静液柱压力；P_f(t)表示t时刻的环空压耗；P_w(t)表示t时刻的井口回压。

由于地层中的气体注入井底后会继续沿环空上返，因此需要考虑气体压缩性，静液柱压力变化也是由于混合物密度的变化。所以式(1)中的P_b(t)可通过用于控压钻井的井筒压力计算模拟的方程计算获得、预测或由PWD直接测量获得，而式(1)中的t时刻的静液柱压力P_h(t)也可实时测量。所述t时刻的静液柱压力P_h(t)可通过以下公式计算获得：

P_h(t)＝ρ_mixgH(t)(2)

${\mathrm{\ρ}}_{mix}=\frac{m_g+m_l}{V\left(t\right)}---\left(3\right)$

式(2)中的ρ_mix表示井筒内混合流体的密度；式(3)中的m_g表示井筒环空中气体质量且m_g＝ρ_gV_g，其中的ρ_g表示平均高压力((P_b+P_w)/2)下的气体密度；m_l表示环空液体质量；V(t)表示t时刻的环空体积并由井身结构及裸眼段直径及入井钻柱体积求得，g表示重力加速度，H表示待测井的井深。其中，所述环空液体质量通过以下公式计算获得：

m_l＝ρ_l(V(t)-V_g(t))(4)

式(4)中的ρ_l表示环空液体体积；V_g(t)表示t时刻的井下溢流量。其中，所述t时刻的井下溢流量V_g(t)通过以下公式计算获得：

$V_g\left(t\right)={\∫}_0^t{q}_g\left(t\right)dt---\left(5\right)$

式(5)中的q_g(t)表示t时刻的溢流速度并通过泥浆罐液面测量得到。对一定深度条件下，公式(2)对时间t进行求导，即可计算获得单位体积下t时刻的静液柱压力：

$\frac{{\mathrm{dP}}_h\left(t\right)}{dt}=-\frac{({\mathrm{\ρ}}_l-{\mathrm{\ρ}}_g)q_g\left(t\right)}{v}gH---\left(6\right)$

进一步地，t时刻的环空压耗P_f(t)可通过以下公式计算获得：

$P_f\left(t\right)=f\frac{{\mathrm{\ρ}}_{mix}{\mathrm{LV}}_{mix}^2\left(t\right)}{2D_{\mathrm{\α}}}---\left(7\right)$

$v_{mix}\left(t\right)=\frac{Q_{mix}\left(t\right)}{A}---\left(8\right)$

式(7)中的f表示摩阻系数，ρ_mix表示井筒内混合流体的密度，D_α表示水力直径且D₀表示井筒直径，D_i表示钻具外径；Q_mix(t)表示质量流量计测量获得的井筒内混合流体的质量，A表示环空过流面积，L表示井筒环空长度，v_mix表示混合物流速。其中，所述摩阻系数通过以下公式计算获得：

$\frac{I}{\sqrt{f}}=-1.08{\log }_{10}\[\frac{6.9}{\mathrm{Re}}+K{\left(\frac{\mathrm{\ϵ}/D_{\mathrm{\α}}}{3.7}\right)}^{1.11}\]---\left(9\right)$

$\mathrm{Re}=\frac{{\mathrm{\ρ}}_{mix}{v}_{mix}\left(t\right)D_{\mathrm{\α}}}{\mathrm{\μ}}---\left(10\right)$

式(9)中的ε/D_α表示相对粗糙度，Re表示雷诺数，用于衡量流动状态的物理量。

由公式(6)可判断钻井过程中静液柱压力变化，由于井下溢流主要是钻井液静液柱压力损失导致的井底压力急剧下降，并引起恶性溢流甚至井喷，因此为了保持井底压力恒定，静液柱压力的变化值ΔP_h(t)表可由以下公式计算：

${\mathrm{\ΔP}}_h\left(t\right)=-\frac{({\mathrm{\ρ}}_l-{\mathrm{\ρ}}_g)V_g}{v}gH---\left(11\right)$

并且，t时刻的井口回压P_w(t)通过以下公式计算获得：

P_w(t)＝P_w0-ΔP_h(t)+ΔP_safe(12)

式(12)中的ΔP_safe表示安全附加压力值；P_w0表示未发生溢流时井口回压。

步骤12，通过压力随钻监测系统PWD获取待测井的实测井底压力值，并基于压力守恒原理对井筒压力中的任意一种进行实时校正，所述井筒压力包括井底压力、环空压耗和井口回压。

本实施例通过采用PWD井下压力随钻测量工具采集的环空压力的方程(即公式(1))进行实时校正，能够大幅度优化提高钻井井筒动态压力流动方程的计算精度，优化后的钻井井筒动态压力流动方程可用于各种工况下的控压钻井井筒动态的水力参数实时计算。PWD快速校正方法有两种，为了两种不同的目的服务。第一种是校正环空压耗的摩阻系数，可迅速计算出准确的环空压耗。第二种是为了快速预测井底压力，计算PWD实测值与井底压力计算值的比值，通过该比值可进行井底压力快速、准确预测分析。下面以对环空压耗、井底压力和井口回压进行实时校正为例进行说明。

校正方式一：校正环空压耗

ΔP(t)＝P_b(t)-P_pwd(t)(13)

P`_f＝P_f(t)-ΔP(t)(14)

则环空摩阻系数校核系数为：

$f^{`}=\frac{{p^{`}}_f}{p_f}f---\left(15\right)$

其中，P_pwd(t)表示PWD测得的井底压力值；ΔP(t)表示计算获得的井底压力与PWD实测值的差值；f`表示环空摩阻校核系数，f表示环空摩阻系数。

由此，代入公式(7)即可计算校正后的环空压耗。

校正方式二：校正井底压力和井口回压

$\mathrm{\α}=\frac{P_{pwd}\left(t\right)}{P_b\left(t\right)}---\left(16\right)$

P`_b＝αP_b(t)(17)

式(16)中的α表示PWD实测压力值与t时刻的井底压力计算值的比值；式(17)中的P`_b表示t+1时刻的井底压力值。根据α值即可计算下一时刻的井底压力值，然后根据公式(1)，即可计算该时刻校正后的井口回压值，通过该井口回压值对节流阀进行控制。

结合图3所示，本实施例提出的实时校正井筒压力的方法的流程包括：首先获取钻井井筒压力计算的基本参数；其次由控压钻井边界条件，即包括已知井身结构、钻具组合及尺寸、钻井液密度、性能等非实时测量参数，以及井口回压、钻井液流量、钻井液循环罐体积变化等实时测量参数；然后根据建立的钻井井筒动态压力流动方程P_b(t)＝P_h(t)+P_f(t)+P_w(t)，计算井底压力和环空压耗；再后通过PWD实测井底压力，利用公式(13)和(15)校正环空摩阻系数并计算获得校正后的环空压耗，并进一步校正井底压力，以便快速进行下部井底压力计算，进而根据公式(1)、(13)、(15)确定合适的井口回压；最后通过实时对比及动态实时获取井底压力、井口回压、钻井液流量、钻井液循环罐体积变化等实时测量参数，从而实现在窄密度窗口地层井筒动态压力的准确计算、实时校正与控制。

进一步优选的，本实施例提出的实时校正井筒压力的方法还可以包括：若无法通过PWD获取待测井的实测井底压力值，则通过实时校正后的井底压力对待测井的环空压耗或井口回压进行校正。

通过采用PWD获取的实测井底压力对井筒压力进行实时校正，能够大幅度优化钻井井筒动态压力流动方程的计算精度，校正后的钻井井筒动态压力流动方程可用于各种工况下的控压钻井井筒动态的水力参数实时计算。但是在某些工况下，例如当井筒出气较多时，导致PWD信号无法传输，可直接通过校正后的钻井井筒动态压力流动方程计算所需井底压力、井口回压等关键参数。下面以计算井口回压为例进行说明：

校正后的钻井井筒动态压力流动方程为：P_b(t)＝P_h(t)+P_f(t)+P_w(t)

首先根据钻井井筒压力基本参数和控压钻井边界条件计算获得t时刻的井底压力P_b(t)，然后计算获得t时刻的静液柱压力P_h(t)，若通过本实施例提出的实时校正井筒压力的方法校正的是t时刻的环空压耗P_f(t)，则可根据校正后的钻井井筒动态压力流动方程计算获得t时刻的井口回压P_w(t)；若通过本实施例提出的实时校正井筒压力的方法校正的是t时刻的井口回压P_w(t)，则可根据校正后的钻井井筒动态压力流动方程计算获得环空压耗P_f(t)。

本具体实施方式提供的实时校正井筒压力的方法具有以下优点：

1、利用PWD实时计算及校正井筒动态压力，具有计算过程快速及计算结果精确的特点，能够计算常规钻井及控压钻井过程中不同井口条件下的井筒动态压力；

2、不需要复杂的网格划分及控制方程组的离散求解，对计算机性能要求较低，更适用于钻井施工的现场工况；

3、利用PWD可以准确、实时的测量井底环空压力等关键数据，利用该测量数据实时校正建立的钻井井筒动态压力流动方程，使得建立的水力计算模型具有较高的精度，并且当井筒出气较多导致PWD信号无法传输时，可直接通过钻井井筒动态压力流动方程计算所需井底及井口回压等关键参数，从而更加精确的控制井筒压力，能够更好的适应井下多变的复杂工况。

本具体实施方式是对本发明的技术方案进行清楚、完整地描述，其中的实施例仅仅是本发明的一部分实施例，而并不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有经过创造性劳动的前提下所获得的所有其它实施方式都属于本发明的保护范围。

Claims (10)

1.一种实时校正井筒压力的方法，其特征在于，包括：

根据钻井井筒压力基本参数和控压钻井边界条件计算获得待测井的校前环空压耗值，并通过计算获得待测井的静液柱压力值；

通过压力随钻监测系统PWD获取待测井的实测井底压力值，并基于压力守恒原理对井筒压力中的任意一种进行实时校正，所述井筒压力包括井底压力、环空压耗和井口回压。

2.如权利要求1所述的方法，其特征在于，所述方法还包括：

若无法通过PWD获取待测井的实测井底压力值，则通过实时校正后的井底压力对待测井的环空压耗或井口回压进行校正。

3.如权利要求1所述的方法，其特征在于，所述井底压力值通过以下公式计算获得：

P_b(t)＝P_h(t)+P_f(t)+P_w(t)

其中，P_b(t)表示t时刻的井底压力；P_h(t)表示t时刻的静液柱压力；P_f(t)表示t时刻的环空压耗；P_w(t)表示t时刻的井口回压。

4.如权利要求3所述的方法，其特征在于，在对所述井筒压力中的任意一种进行实时校正的过程中，对所述环压空耗的校正包括：

将环空摩阻校核系数带入所述环空压耗的计算公式，对所述环压空耗进行校正，所述环空摩阻校核系数通过以下公式计算获得：

ΔP(t)＝P_b(t)-P_pwd(t)

P`_f＝P_f(t)-ΔP(t)

$f^{`}=\frac{{p^{`}}_f}{p_f}f$

其中，P_pwd(t)表示PWD在t时刻获取的实测井底压力值；ΔP(t)表示计算获得的井底压力与PWD在t时刻获取的实测井底压力值之间的差值；f`表示环空摩阻校核系数，f表示环空摩阻系数。

5.如权利要求4所述的方法，其特征在于，所述环空压耗通过以下公式计算获得：

$P_f\left(t\right)=f\frac{{\mathrm{\ρ}}_{mix}{\mathrm{Lv}}_{mix}^2(tj}{2D_{\mathrm{\α}}}$

$v_{mix}\left(t\right)=\frac{Q_{mix}\left(t\right)}{A}$

$D_{\mathrm{\α}}=\frac{D_o-D_i}{2}$

其中，ρ_mix表示井筒内混合流体的密度，D_α表示水力直径，D₀表示井筒直径，D_i表示钻具外径；Q_mix(t)表示质量流量计的测量值，A表示环空过流面积，L表示井筒环空长度，v_mix表示混合物流速。

6.如权利要求5所述的方法，其特征在于，所述环空摩阻系数通过以下公式计算获得：

$\frac{I}{\sqrt{f}}=-1.08{\log }_{10}\[\frac{6.9}{\mathrm{Re}}+{\left(\frac{\mathrm{\ϵ}/D_{\mathrm{\α}}}{3.7}\right)}^{1.11}\]$

$\mathrm{Re}=\frac{{\mathrm{\ρ}}_{mix}{v}_{mix}\left(t\right)D_{\mathrm{\α}}}{\mathrm{\μ}}$

其中，ε/D_α表示相对粗糙度，Re表示雷诺数。

7.如权利要求3所述的方法，其特征在于，在对所述井筒压力中的任意一种进行实时校正的过程中，对所述井底压力的校正包括：

将t+1时刻的井底压力带入所述井底压力的计算公式，对所述井底压力进行校正，所述t+1时刻的井底压力通过以下公式计算获得：

$\mathrm{\α}=\frac{P_{pwd}\left(t\right)}{P_b\left(t\right)}$

P`_b＝αP_b(t)

其中，α表示PWD在t时刻获取的实测井底压力值与t时刻的井底压力之间的比值；P`_b表示t+1时刻的井底压力。

8.如权利要求7所述的方法，其特征在于，所述井口回压通过以下公式计算获得：

P_w(t)＝P_w0-ΔP_h(t)+ΔP_safe

${\mathrm{\ΔP}}_h\left(t\right)=-\frac{({\mathrm{\ρ}}_l-{\mathrm{\ρ}}_g)V_g}{v}gH$

其中，P_w0表示未发生溢流时的井口回压值，ΔP_h(t)表示静液柱压力的变化值，ΔP_safe表示安全附加压力值，ρ_l表示环空液体密度，ρ_g表示平均高压力下的气体密度，V_g表示井下溢流量；g表示重力加速度，H表示待测井的井深。

9.如权利要求3所述的方法，其特征在于，所述静液柱压力通过以下公式计算获得：

P_h(t)＝ρ_mixgH(t)

${\mathrm{\ρ}}_{mix}=\frac{m_g+m_l}{V\left(t\right)}$

m_g＝ρ_gV_g

其中，ρ_mix表示井筒内混合流体的密度；m_g表示井筒环空中气体质量且，ρ_g表示平均高压力下的气体密度，V_g表示井下溢流量，m_l表示环空液体质量；V(t)表示t时刻的环空体积并由井身结构、裸眼段直径及入井钻柱体积求得。

10.一种基于权利要求1至9任意一项所述的实时校正井筒压力的方法的井筒压力测量装置，其特征在于，包括：

PWD，用于测量井底压力；

压力测量计，用于测量井口回压；

质量流量计，用于测量井筒内混合流体的质量；

液位计，用于测量井下溢流量。