一种深水测试井筒压力、温度场模拟及水合物预测方法
一种深水测试井筒压力、温度场模拟及水合物预测方法
【技术领域】
[0001] 本发明涉及一种深水测试井筒压力、温度场及水合物预测方法，特别是关于一种 在海水温度影响下井筒稳态传热和非稳态传热下的深水测试井筒压力、温度场模拟及水合 物预测方法。
【背景技术】
[0002] 目前我国能源开发形式日益严峻，海洋深水区域逐渐成为勘探开发的焦点。但对 于深水油气田来说，普遍具有埋藏深，储层高温、高压、渗透性好等特征，即正常情况下，深 海油气田中的生产井往往是高产井，但是伴随着高产的条件下，高流量气体结合海水的降 温作用以及针对高压力的多级节流，很容易在井筒中出现水合物，堵塞测试管柱和生产管 柱。但是目前还没有一套完整的深水测试井筒压力、温度场模拟及水合物预测的技术方法。
[0003] 井筒压力剖面计算极大地依赖于井筒温度剖面的预测，目前常规方法是假设井筒 温度剖面线性化，井筒流体温度与时间无关，这种处理方式对处于稳定生产状态的油气井 压力计算是合适的，但是在油气井不稳定测试过程中流量、压力、温度可能都处于不稳态过 程，在此情况下沿用常规方法估计井筒温度误差较大，导致压力计算偏差。针对井筒动态压 力和温度场模拟技术及水合物预测的技术，国外已经形成有较为成熟商业软件-PIPESIM。 该软件中的黑油模型可以对油、气、水三相，气液两相，以及单相流体进行计算井筒温度场、 压力剖面模拟。但是通过PIPESIM软件进行井筒温度场、压力剖面模拟预测，其存在以下缺 点：1、井筒压力剖面预测是基于假设井筒温度剖面为线性变化的计算上进行的，对于油气 井不稳定测试过程中，计算结果偏差大。2、没有考虑试油测试短期过程中井筒与地层的非 稳态传热。3、没有考虑海水温度变化。
【发明内容】
[0004] 针对上述问题，本发明的目的是提供一种深水测试井筒压力、温度场模拟及水合 物预测方法，该方法对试采阶段的长期稳定生产过程，采用解析方法建立井筒稳态传热模 型，预测非线性井温剖面；针对试油测试的短期过程，建立井筒非稳态传热模型，预测不同 测试制度和时间下的井温剖面；并能计算井筒温度，压力结果可靠，准确度高。
[0005] 本发明的另一目的是在温度剖面预测基础上，进行压力剖面和水合物预测。
[0006] 为实现上述目的，本发明采取以下技术方案：一种深水测试井筒压力、温度场模拟 及水合物预测方法，其包括以下步骤：1)准备参数：气藏参数、流体物性参数、井身结构参 数和海水温度剖面等；2)对井筒进行离散化，并对井筒温度场和压力场进行初始化，初始 化井筒温度等于环境温度；3)利用气藏渗流数学模型求解地层流量和井底压力分布信息； 4)将地层流量和井底压力分布信息输入井筒流动模型内进行迭代，进而求解出井筒压力剖 面，并判断井筒压力是否收敛，若收敛则进入步骤5)，反之继续判断井筒压力是否收敛；其 中，井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流的井 筒压力计算管流模型；5)计算测试井内流体和环空流体在当前温度压力条件下的热物性 参数，并计算地层段和海水段的总传热系数；6)将井筒压力和总传热系数代入井筒温度场 模型，从井底往井口递推计算各节点井筒温度剖面，并修正地层段和海水段的环空温度差， 更新流体热物性参数；各节点是井筒每离散段的节点；7)判断井筒温度剖面是否收敛，若 不收敛，则返回步骤6)进行迭代，直至井筒温度剖面收敛；若收敛，则修正井筒的压力和温 度；8)判断tn是否小于tn+1，若小于则输出井温-井深曲线以及井筒压力对应的水合物生 成温度-井深曲线，并根据曲线判断是否有水合生成；反之若不小于，则返回步骤3)，预测 下一时刻的井筒压力、温度场模拟及水合物情况。
[0007] 所述步骤3)中，所述气藏渗流数学模型为：3.1)无限大均质油气藏中心直井 Laplace空间压力相应解为：
[0008]
[0009] 式中，为井底Laplace空间压力；u=sf(s)，其中s为Laplace空间时间变量， f(s) = 1 &分别为修正的零阶和一阶第二类Bessel函数；S为表皮系数；CD为无因次 井筒储集系数；3. 2)无限大油藏斜井在Laplace空间压力相应解：
[0010]
[0011] 式中，hfD为无因次射孔井段长度；s为Laplace空间时间变量；&为无因次半径； n为级数系数；hD为无因次储层厚度；zD为井筒深度z的无因次值、为井段中心井筒深度z的无因次值。
[0012] 所述步骤4)中，所述井筒流动模型包括海水温度场效应的稳态和非稳态下的井 筒温度场模型、多相流的井筒压力计算管流模型。
[0013] 所述稳态井筒温度场模型为：
[0014] 地层段的流体温度乙为：
[0015] n [_
F J 乂、
[0016] 式中，z为井筒深度；步为倾斜校正系数；
[0017]
[0018] 海水段的流体温度Tf表达式：
[0019]
[0020] 式中，F。为焦耳-汤普森效应与动能作用系数；Tf(l为井筒入口温度；Tei(lS地层温 度；Z为井筒深度；LK为松弛距离参数；9为井筒与水平面的夹角；g为重力加速度；It3为 组合参数；定义组合参数t为：
[0021] 所述非稳态井筒温度场模型为： /-* -nt\
[0022]
[0023] 式中，a为温度随时间变化的衰减系数；1^为任意深度处的原始地层温度；Cp为油 管流体比热，单位为J/kg-K;LK为松弛距离参数；L为垂直井筒总长度，单位为m;ge为地温 梯度，单位为K/m;gse为海水段温度梯度，单位为K/m;t为时间，单位为s。
[0024] 所述多相流的井筒压力计算管流模型建立方法如下：多相管流模型的基本方程 为：
[0025]
, (dP\ (dP\
[0026] 式中，为位差压力梯度；7为摩檫压力梯度；为加速度 dZ止 JmwV必y加速度 Ap Ap 压力梯度；为总压力梯度；其中，总压力梯度的计算步骤如下：a)计算弗鲁德aZaZ 数、无滑脱持液率^和LpL2,L3,L4:
[0027]
[0028]
[0029]L1= 316E^ 302 ；
[0030] L2= 0. 0009252EL^ 4682；
[0031] L3= 0? 1EL-14516;
[0032] L4= 0? 5E,.738;
[0033] 式中，匕为无滑脱持液率；为液体体积流量；（2)根据步骤⑴中计算的基本 参数，采用Beggs-Brill法判断流型：分离流、间歇流或分散流：当匕< 0. 01时Nft<L 0.01时Nft<L2Jj流型为离流；当0.01，L2<K1^时，流型为过渡流；当 0? 01 彡EL< 0? 4 时L3<NFr<L丨，EL彡 0? 4 时L3<NFr<L4，流型为间歇流；当EL< 0? 4 时Lp彡0. 4时NL4，流型为分散流；（3)持液率计算：倾角为0的气液两相流 动的持液率扎（0 )为：扎（0 )=扎(0)步；其中，扎(0)为同样流动参数下，水平流动时的持 液率；
[0034]
[0035]式中，步为倾斜校正系数，!，= 1 + (.1>丨11(1.86))-^如_5(1.册)]，对于垂直管：1])= 1+0. 3C，C为系数，与无滑脱持液率、弗鲁德系数和液相速度有关；b、c均为流型常数；(4)计 算液相密度PyPPJt+Pg(l-H) ;(5)确定流动阻力系数入：
[0040] 其中，h为在P、T下的液体粘度；yg为在P、T下的气体粘度；T为温度；(6)计 算井筒压力梯度：将步骤（1)~步骤（5)中计算的参数代入多相管流模型的基本方程中计 算得到井筒压力梯度。
[0041] 所述步骤5)中，所述总传热系数采用热阻串联法计算。
[0042] 本发明由于采取以上技术方案，其具有以下优点：1、本发明采用针对油气井不稳 定测试过程中流量、压力、温度可能都处于的不稳态过程，建立了井筒非稳态传热模型，计 算井筒温度、压力结果可靠，准确度高。2、本发明为了克服传统井筒温度、压力场预测方法 存在的缺陷，同时解决深水产层高流量气体结合海水的降温作用以及针对高压力的多级节 流，很容易在井筒中出现水合物的问题，本发明对试采阶段的长期稳定生产过程，采用解析 方法建立井筒稳态传热模型，预测的非线性井温剖面；针对试油测试的短期过程，建立井筒 非稳态传热模型，预测不同测试制度和时间下的井温剖面。进一步在温度剖面预测基础上， 开展压力剖面和水合物预测。本发明可以广泛在石油勘探过程中应用。
【附图说明】
[0043] 图1是本发明的整体流程示意图。
【具体实施方式】
[0044] 下面结合附图和实施例对本发明进行详细的描述。
[0045] 如图1所示，本发明提供一种深水测试井筒压力、温度场模拟及水合物预测方法， 利用气藏渗流数学模型和井筒流动模型，将两个模型耦合求解地层和井筒的瞬态温度、压 力剖面，进而实现对井筒压力、温度场模拟及水合物的预测，其包括以下步骤：
[0046] 1)准备参数：气藏参数、流体物性参数、井身结构参数和海水温度剖面等；
[0047] 2)对井筒进行离散化，并对井筒温度场和压力场进行初始化，初始化井筒温度等 于环境温度；
[0048] 3)利用气藏渗流数学模型求解地层流量和井底压力分布信息；其中，气藏渗流数 学模型为：
[0049] 利用渗流力学原理和质量守恒定律建立气藏不同完井方式（裸眼、射孔、筛管 等）、边界条件（封闭、定压、无限大）和不同井型（直井、斜井、水平井）下的渗流数学模 型。其典型模型如下：
[0050] 3. 1)无限大均质油气藏中心直井Laplace空间压力相应解为：
[0051](1) **-1v* / '一n
-uv ~ 1 /j^
[0052] 式中，之为井底Laplace空间压力；u=sf(s)，其中s为Laplace空间时间变量， f(s) = 1 &分别为修正的零阶和一阶第二类Bessel函数；S为表皮系数；CD为无因次 井筒储集系数。
[0053] 3. 2)无限大油藏斜井在Laplace空间压力相应解：
[0054]
[0055] 式中，为无因次射孔井段长度；s为Laplace空间时间变量；&为无因次半径；n 为级数系数；hD为无因次储层厚度；zD为井筒深度z的无因次值、为井段中心井筒深度 z的无因次值。
[0056] 4)将地层流量和井底压力分布信息输入井筒流动模型内进行迭代，进而求解出井 筒压力剖面，并判断井筒压力是否收敛，若收敛则进入步骤5)，反之继续判断井筒压力是否 收敛；
[0057] 其中，井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、 多相流的井筒压力计算管流模型；
[0058] 5)计算测试井内流体和环空流体在当前温度压力条件下的热物性参数，并计算地 层段和海水段的总传热系数；其中，总传热系数采用热阻串联法计算；
[0059] 6)将井筒压力和总传热系数代入井筒温度场模型，从井底往井口递推计算各节点 井筒温度剖面，并修正地层段和海水段的环空温度差，更新流体热物性参数；其中，各节点 是井筒每离散段的节点；
[0060] 7)判断井筒温度剖面是否收敛，若不收敛，则返回步骤6)进行迭代，直至井筒温 度剖面收敛，一般迭代3~5次；若收敛，则修正井筒的压力和温度；
[0061] 8)判断tn是否小于丨1"1，若小于则输出井温-井深曲线以及井筒压力对应的水合 物生成温度-井深曲线，并根据曲线判断是否有水合生成；反之若不小于，则返回步骤3)， 预测下一时刻的井筒压力、温度场模拟及水合物情况；
[0062] 其中，在计算第n+1个时间步长时，将用第n个时间步长的井筒温度场和压力场作 为初始条件。
[0063] 上述步骤4)中，稳态和非稳态下的井筒温度场模型建立方法如下：
[0064] 由于本发明针对试采阶段的长期稳定生产过程，因此采用解析方法建立稳态井筒 温度场模型，预测非线性井温剖面；针对试油测试的短期过程，建立非稳态井筒温度场模 型，预测不同测试制度和时间下的井温剖面。根据能量守恒定律、焦耳-汤普森效应、井筒 径向传热以及海上温度场影响，得到稳态井筒温度场模型和非稳态井筒温度场模型。
[0065] 4. 1)稳态井筒温度场模型
[0066] 地层段的流体温度乙为：
[0067]
[0068] 式中，z为井筒深度；步为倾斜校正系数。
[0069]
⑷
[0070] 海水段的流体温度Tf表达式：
[0071]
[0072] 式中，F。为焦耳-汤普森效应与动能作用系数；Tf(l为井筒入口温度；Tei(lS地层温 度；z为井筒深度；LK为松弛距离参数；9为井筒与水平面的夹角；g为重力加速度；步3为 组合参数；
[0073] 定义组合参数$s为：
[0074]
(6)
[0075] 4. 2)非稳态井筒温度场模型为：
[0076]
(7)
[0077] 式⑶~（7)中，a为温度随时间变化的衰减系数；Tei为任意深度处的原始地层 温度；Cp为油管流体比热，单位为J/kg-K;LK为松弛距离参数；L为垂直井筒总长度，单位为 m;ge为地温梯度，单位为K/m;gse为海水段温度梯度，单位为K/m;t为时间，单位为s。
[0078] 上述步骤4)中，多相流的井筒压力计算管流模型建立方法如下：
[0079]由于海上油气田一般是大斜度井或水平井，同时结合多种管流计算模型的适应性 条件，本发明优选Beggs和Brill方法作为井筒压力计算模型。本发明中Beggs和Brill 根据在长为15m、直径为25. 4mm和38. 1mm的管子中，用空气和水进行实验的基础上提出了 处理所有范围的多相流及任何角度的管流的普适化相关式。其中，Beggs和Brill多相管 流模型的基本方程、流型判别条件以及主要参数确定如下：
[0080] 根据机械能量守恒定律，气液两相气液两相管流的压力降消耗于三个方面：位差、 摩擦和加速度，即多相管流模型的基本方程为：
[0081]
[0082] 式中各项参数的解释如表1所示。 (8)
[0083] 表1式⑶中各项参数的解释表
[0084]
[0085] 其中，总压力梯度的计算步骤如下： dZ
[0086] (1)计算弗鲁德数NFr、无滑脱持液率匕和L^L2,L3,L4:
[0087]
[0088]
[0089] L:= 316E^ 302 ；
[0090] L2= 0. 0009252EL^ 4682；
[0091] L3= 0? 1El_14516;
[0092] L4= 0? 5EJ6.738;
[0093] 式中，无滑脱持液率；Q液体体积流量；
[0094] (2)根据步骤（1)中计算的基本参数，采用Beggs-Brill法判断流型：分离流、间 歇流或分散流。
[0095]当EL< 0.01时NFr<LpEl彡0.01时NFr<L2，则流型为离流；
[0096] 当El彡0? 01，L2<NFr<L3时，流型为过渡流；
[0097]当 0? 01 彡EL< 0? 4 时L3<NFr<L丨，EL彡 0? 4 时L3<NFr<L4，流型为间歇流；
[0098] 当EL<0?4时NFr彡L:，EL彡0?4时NFr>L4,流型为分散流。
[0099] (3)持液率 计算
[0100] 倾角为0的气液两相流动的持液率扎（9 )为：
[0101] HL(e) =Hl(0)It (9)
[0102] 式（9)中，扎(0)为同样流动参数下，水平流动时的持液率。
[0103]
(1〇)
[0104] 式（9)中，$为倾斜校正系1
对于垂直管：步 =1+0. 3C，C为系数，与无滑脱持液率、弗鲁德系数和液相速度有关；b、c均为流型常数；
[0105] (4)计算液相密度Pl:
[0106] pL=pLHL+Pg(l-HL) (11)
[0107] (5)确定流动阻力系数入：
[0108]
"2)
[0109] 式中：
[0110]
[0111]
[0112] 其中，h为在P、T下的液体粘度；yg为在P、T下的气体粘度，T为温度。
[0113] (6)计算井筒压力梯度：将步骤（1)~步骤（5)中计算的参数代入式（8)中计算 得到井筒压力梯度。
[0114] 上述各实施例仅用于说明本发明，各步骤都是可以有所变化的，在本发明技术方 案的基础上，凡根据本发明原理对个别步骤进行的改进和等同变换，均不应排除在本发明 的保护范围之外。
【主权项】
1. 一种深水测试井筒压力、温度场模拟及水合物预测方法，其包括以下步骤： 1) 准备参数：气藏参数、流体物性参数、井身结构参数和海水温度剖面等； 2) 对井筒进行离散化，并对井筒温度场和压力场进行初始化，初始化井筒温度等于环 境温度； 3) 利用气藏渗流数学模型求解地层流量和井底压力分布信息； 4) 将地层流量和井底压力分布信息输入井筒流动模型内进行迭代，进而求解出井筒 压力剖面，并判断井筒压力是否收敛，若收敛则进入步骤5)，反之继续判断井筒压力是否收 敛；其中，井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流 的井筒压力计算管流模型； 5) 计算测试井内流体和环空流体在当前温度压力条件下的热物性参数，并计算地层段 和海水段的总传热系数； 6) 将井筒压力和总传热系数代入井筒温度场模型，从井底往井口递推计算各节点井筒 温度剖面，并修正地层段和海水段的环空温度差，更新流体热物性参数；各节点是井筒每离 散段的节点； 7) 判断井筒温度剖面是否收敛，若不收敛，则返回步骤6)进行迭代，直至井筒温度剖 面收敛；若收敛，则修正井筒的压力和温度； 8) 判断tn是否小于t n+1，若小于则输出井温-井深曲线以及井筒压力对应的水合物生 成温度-井深曲线，并根据曲线判断是否有水合生成；反之若不小于，则返回步骤3)，预测 下一时刻的井筒压力、温度场模拟及水合物情况。2. 如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述步骤3)中，所述气藏渗流数学模型为： 3. 1)无限大均质油气藏中心直井Laplace空间压力相应解为：式中，九ο为井底Laplace空间压力；u = sf(s)，其中s为Laplace空间时间变量，f (s) =1 ;&、K1分别为修正的零阶和一阶第二类Bessel函数；S为表皮系数；C D为无因次井筒 储集系数； 3. 2)无限大油藏斜井在Laplace空间压力相应解：式中，hfD为无因次射孔井段长度；s为Laplace空间时间变量；&为无因次半径；η为 级数系数；hD为无因次储层厚度；ζ D为井筒深度ζ的无因次值、为井段中心井筒深度ζ的无因次值。3.如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述步骤4)中，所述井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒 温度场模型、多相流的井筒压力计算管流模型。4. 如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述稳态井筒温度场模型为： 地层段的流体温度乙为：式中，z为井筒深度；Φ为倾斜校正系数；海水段的流体温度Tf表达式：式中，F。为焦耳-汤普森效应与动能作用系数；Tftl为井筒入口温度；Teitl为地层温度；z 为井筒深度；LK为松弛距离参数；Θ为井筒与水平面的夹角；g为重力加速度；φ 3为组合参 数； 定义组合参数Φ 3为：5. 如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述非稳态井筒温度场模型为：式中，a为温度随时间变化的衰减系数；Tei为任意深度处的原始地层温度；CpS油管流 体比热，单位为J/kg-K ;Lk为松弛距离参数；L为垂直井筒总长度，单位为m ;ge为地温梯度， 单位为K/m ;gse为海水段温度梯度，单位为K/m ;t为时间，单位为s。6. 如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述多相流的井筒压力计算管流模型建立方法如下： 多相管流模型的基本方程为：/ (?Ρλ ((IP λ 式中，（t)为位差压力梯度；为摩擦压力梯度；^ 为加速度压力梯 dZ Jmm V似夕加速度 AP 度；为总压力梯度； dZ dP 其中，总压力梯度的计算步骤如下： CiZ (1) 计算弗鲁德数Nft、无滑脱持液率EjP L i，L2, L3, L4:式中，无滑脱持液率；Ql为液体体积流量； (2) 根据步骤（1)中计算的基本参数，采用Beggs-Brill法判断流型：分离流、间歇流 或分散流： 当EL< 0· Ol时N Fr< L p El彡0· Ol时N Fr< L 2，则流型为离流； 当0. 01，L 2< N L 3时，流型为过渡流； 当 0· 01 彡 EL< 0· 4 时 L 3< N Fr< L p El彡 0· 4 时 L 3< N Fr< L 4，流型为间歇流； 当EL< 0· 4时N Fr彡L i，0· 4时N Fr> L 4，流型为分散流； (3) 持液率计算 倾角为Θ的气液两相流动的持液率扎（0)为： HL( Θ )= Hl(O) Φ， 其中，扎(0)为同样流动参数下，水平流动时的持液率；式中，Φ为倾斜校正系数，对于垂直管：Φ = 1+0. 3C，C为系数，与无滑脱持液率、弗鲁德系数和液相速度有关；b、c均为流型常数； (4) 计算液相密度P ^l= f LhL+Ps(I-Hl); (5) 确定流动阻力系数λ :其中，μ ^3在P、T下的液体粘度；μ g为在P、T下的气体粘度；T为温度； (6) 计算井筒压力梯度：将步骤（1)~步骤（5)中计算的参数代入多相管流模型的基 本方程中计算得到井筒压力梯度。7.如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述步骤5)中，所述总传热系数采用热阻串联法计算。
【专利摘要】本发明涉及一种深水测试井筒压力、温度场模拟及水合物预测方法，其步骤：准备参数；对井筒温度场和压力场进行初始化；利用气藏渗流数学模型求解地层流量和井底压力分布信息；将地层流量和井底压力分布信息输入井筒流动模型内进行迭代，求解出井筒压力剖面并判断井筒压力是否收敛；计算热物性参数，并计算地层段和海水段的总传热系数；将井筒压力和总传热系数代入井筒温度场模型，计算各节点井筒温度剖面；判断井筒温度剖面是否收敛；判断tn是否小于tn+1，若小于则输出井温-井深曲线以及井筒压力对应的水合物生成温度-井深曲线，并根据曲线判断是否有水合生成；反之则返回预测下一时刻的井筒压力、温度场模拟及水合物情况。本发明可以广泛在石油勘探过程中应用。
【IPC分类】E21B49/00, E21B47/07, E21B47/06
【公开号】CN104895560
【申请号】CN201510333445
【发明人】何玉发, 蒋世全, 段永刚, 周建良, 李迅科, 杨秀夫, 陈伟, 方全堂
【申请人】中国海洋石油总公司, 中海油研究总院
【公开日】2015年9月9日
【申请日】2015年6月16日
【技术领域】
[0001] 本发明涉及一种深水测试井筒压力、温度场及水合物预测方法，特别是关于一种 在海水温度影响下井筒稳态传热和非稳态传热下的深水测试井筒压力、温度场模拟及水合 物预测方法。
【背景技术】
[0002] 目前我国能源开发形式日益严峻，海洋深水区域逐渐成为勘探开发的焦点。但对 于深水油气田来说，普遍具有埋藏深，储层高温、高压、渗透性好等特征，即正常情况下，深 海油气田中的生产井往往是高产井，但是伴随着高产的条件下，高流量气体结合海水的降 温作用以及针对高压力的多级节流，很容易在井筒中出现水合物，堵塞测试管柱和生产管 柱。但是目前还没有一套完整的深水测试井筒压力、温度场模拟及水合物预测的技术方法。
[0003] 井筒压力剖面计算极大地依赖于井筒温度剖面的预测，目前常规方法是假设井筒 温度剖面线性化，井筒流体温度与时间无关，这种处理方式对处于稳定生产状态的油气井 压力计算是合适的，但是在油气井不稳定测试过程中流量、压力、温度可能都处于不稳态过 程，在此情况下沿用常规方法估计井筒温度误差较大，导致压力计算偏差。针对井筒动态压 力和温度场模拟技术及水合物预测的技术，国外已经形成有较为成熟商业软件-PIPESIM。 该软件中的黑油模型可以对油、气、水三相，气液两相，以及单相流体进行计算井筒温度场、 压力剖面模拟。但是通过PIPESIM软件进行井筒温度场、压力剖面模拟预测，其存在以下缺 点：1、井筒压力剖面预测是基于假设井筒温度剖面为线性变化的计算上进行的，对于油气 井不稳定测试过程中，计算结果偏差大。2、没有考虑试油测试短期过程中井筒与地层的非 稳态传热。3、没有考虑海水温度变化。
【发明内容】
[0004] 针对上述问题，本发明的目的是提供一种深水测试井筒压力、温度场模拟及水合 物预测方法，该方法对试采阶段的长期稳定生产过程，采用解析方法建立井筒稳态传热模 型，预测非线性井温剖面；针对试油测试的短期过程，建立井筒非稳态传热模型，预测不同 测试制度和时间下的井温剖面；并能计算井筒温度，压力结果可靠，准确度高。
[0005] 本发明的另一目的是在温度剖面预测基础上，进行压力剖面和水合物预测。
[0006] 为实现上述目的，本发明采取以下技术方案：一种深水测试井筒压力、温度场模拟 及水合物预测方法，其包括以下步骤：1)准备参数：气藏参数、流体物性参数、井身结构参 数和海水温度剖面等；2)对井筒进行离散化，并对井筒温度场和压力场进行初始化，初始 化井筒温度等于环境温度；3)利用气藏渗流数学模型求解地层流量和井底压力分布信息； 4)将地层流量和井底压力分布信息输入井筒流动模型内进行迭代，进而求解出井筒压力剖 面，并判断井筒压力是否收敛，若收敛则进入步骤5)，反之继续判断井筒压力是否收敛；其 中，井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流的井 筒压力计算管流模型；5)计算测试井内流体和环空流体在当前温度压力条件下的热物性 参数，并计算地层段和海水段的总传热系数；6)将井筒压力和总传热系数代入井筒温度场 模型，从井底往井口递推计算各节点井筒温度剖面，并修正地层段和海水段的环空温度差， 更新流体热物性参数；各节点是井筒每离散段的节点；7)判断井筒温度剖面是否收敛，若 不收敛，则返回步骤6)进行迭代，直至井筒温度剖面收敛；若收敛，则修正井筒的压力和温 度；8)判断tn是否小于tn+1，若小于则输出井温-井深曲线以及井筒压力对应的水合物生 成温度-井深曲线，并根据曲线判断是否有水合生成；反之若不小于，则返回步骤3)，预测 下一时刻的井筒压力、温度场模拟及水合物情况。
[0007] 所述步骤3)中，所述气藏渗流数学模型为：3.1)无限大均质油气藏中心直井 Laplace空间压力相应解为：
[0008]
[0009] 式中，为井底Laplace空间压力；u=sf(s)，其中s为Laplace空间时间变量， f(s) = 1 &分别为修正的零阶和一阶第二类Bessel函数；S为表皮系数；CD为无因次 井筒储集系数；3. 2)无限大油藏斜井在Laplace空间压力相应解：
[0010]
[0011] 式中，hfD为无因次射孔井段长度；s为Laplace空间时间变量；&为无因次半径； n为级数系数；hD为无因次储层厚度；zD为井筒深度z的无因次值、为井段中心井筒深度z的无因次值。
[0012] 所述步骤4)中，所述井筒流动模型包括海水温度场效应的稳态和非稳态下的井 筒温度场模型、多相流的井筒压力计算管流模型。
[0013] 所述稳态井筒温度场模型为：
[0014] 地层段的流体温度乙为：
[0015] n [_
F J 乂、
[0016] 式中，z为井筒深度；步为倾斜校正系数；
[0017]
[0018] 海水段的流体温度Tf表达式：
[0019]
[0020] 式中，F。为焦耳-汤普森效应与动能作用系数；Tf(l为井筒入口温度；Tei(lS地层温 度；Z为井筒深度；LK为松弛距离参数；9为井筒与水平面的夹角；g为重力加速度；It3为 组合参数；定义组合参数t为：
[0021] 所述非稳态井筒温度场模型为： /-* -nt\
[0022]
[0023] 式中，a为温度随时间变化的衰减系数；1^为任意深度处的原始地层温度；Cp为油 管流体比热，单位为J/kg-K;LK为松弛距离参数；L为垂直井筒总长度，单位为m;ge为地温 梯度，单位为K/m;gse为海水段温度梯度，单位为K/m;t为时间，单位为s。
[0024] 所述多相流的井筒压力计算管流模型建立方法如下：多相管流模型的基本方程 为：
[0025]
, (dP\ (dP\
[0026] 式中，为位差压力梯度；7为摩檫压力梯度；为加速度 dZ止 JmwV必y加速度 Ap Ap 压力梯度；为总压力梯度；其中，总压力梯度的计算步骤如下：a)计算弗鲁德aZaZ 数、无滑脱持液率^和LpL2,L3,L4:
[0027]
[0028]
[0029]L1= 316E^ 302 ；
[0030] L2= 0. 0009252EL^ 4682；
[0031] L3= 0? 1EL-14516;
[0032] L4= 0? 5E,.738;
[0033] 式中，匕为无滑脱持液率；为液体体积流量；（2)根据步骤⑴中计算的基本 参数，采用Beggs-Brill法判断流型：分离流、间歇流或分散流：当匕< 0. 01时Nft<L 0.01时Nft<L2Jj流型为离流；当0.01，L2<K1^时，流型为过渡流；当 0? 01 彡EL< 0? 4 时L3<NFr<L丨，EL彡 0? 4 时L3<NFr<L4，流型为间歇流；当EL< 0? 4 时Lp彡0. 4时NL4，流型为分散流；（3)持液率计算：倾角为0的气液两相流 动的持液率扎（0 )为：扎（0 )=扎(0)步；其中，扎(0)为同样流动参数下，水平流动时的持 液率；
[0034]
[0035]式中，步为倾斜校正系数，!，= 1 + (.1>丨11(1.86))-^如_5(1.册)]，对于垂直管：1])= 1+0. 3C，C为系数，与无滑脱持液率、弗鲁德系数和液相速度有关；b、c均为流型常数；(4)计 算液相密度PyPPJt+Pg(l-H) ;(5)确定流动阻力系数入：
[0040] 其中，h为在P、T下的液体粘度；yg为在P、T下的气体粘度；T为温度；(6)计 算井筒压力梯度：将步骤（1)~步骤（5)中计算的参数代入多相管流模型的基本方程中计 算得到井筒压力梯度。
[0041] 所述步骤5)中，所述总传热系数采用热阻串联法计算。
[0042] 本发明由于采取以上技术方案，其具有以下优点：1、本发明采用针对油气井不稳 定测试过程中流量、压力、温度可能都处于的不稳态过程，建立了井筒非稳态传热模型，计 算井筒温度、压力结果可靠，准确度高。2、本发明为了克服传统井筒温度、压力场预测方法 存在的缺陷，同时解决深水产层高流量气体结合海水的降温作用以及针对高压力的多级节 流，很容易在井筒中出现水合物的问题，本发明对试采阶段的长期稳定生产过程，采用解析 方法建立井筒稳态传热模型，预测的非线性井温剖面；针对试油测试的短期过程，建立井筒 非稳态传热模型，预测不同测试制度和时间下的井温剖面。进一步在温度剖面预测基础上， 开展压力剖面和水合物预测。本发明可以广泛在石油勘探过程中应用。
【附图说明】
[0043] 图1是本发明的整体流程示意图。
【具体实施方式】
[0044] 下面结合附图和实施例对本发明进行详细的描述。
[0045] 如图1所示，本发明提供一种深水测试井筒压力、温度场模拟及水合物预测方法， 利用气藏渗流数学模型和井筒流动模型，将两个模型耦合求解地层和井筒的瞬态温度、压 力剖面，进而实现对井筒压力、温度场模拟及水合物的预测，其包括以下步骤：
[0046] 1)准备参数：气藏参数、流体物性参数、井身结构参数和海水温度剖面等；
[0047] 2)对井筒进行离散化，并对井筒温度场和压力场进行初始化，初始化井筒温度等 于环境温度；
[0048] 3)利用气藏渗流数学模型求解地层流量和井底压力分布信息；其中，气藏渗流数 学模型为：
[0049] 利用渗流力学原理和质量守恒定律建立气藏不同完井方式（裸眼、射孔、筛管 等）、边界条件（封闭、定压、无限大）和不同井型（直井、斜井、水平井）下的渗流数学模 型。其典型模型如下：
[0050] 3. 1)无限大均质油气藏中心直井Laplace空间压力相应解为：
[0051](1) **-1v* / '一n
-uv ~ 1 /j^
[0052] 式中，之为井底Laplace空间压力；u=sf(s)，其中s为Laplace空间时间变量， f(s) = 1 &分别为修正的零阶和一阶第二类Bessel函数；S为表皮系数；CD为无因次 井筒储集系数。
[0053] 3. 2)无限大油藏斜井在Laplace空间压力相应解：
[0054]
[0055] 式中，为无因次射孔井段长度；s为Laplace空间时间变量；&为无因次半径；n 为级数系数；hD为无因次储层厚度；zD为井筒深度z的无因次值、为井段中心井筒深度 z的无因次值。
[0056] 4)将地层流量和井底压力分布信息输入井筒流动模型内进行迭代，进而求解出井 筒压力剖面，并判断井筒压力是否收敛，若收敛则进入步骤5)，反之继续判断井筒压力是否 收敛；
[0057] 其中，井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、 多相流的井筒压力计算管流模型；
[0058] 5)计算测试井内流体和环空流体在当前温度压力条件下的热物性参数，并计算地 层段和海水段的总传热系数；其中，总传热系数采用热阻串联法计算；
[0059] 6)将井筒压力和总传热系数代入井筒温度场模型，从井底往井口递推计算各节点 井筒温度剖面，并修正地层段和海水段的环空温度差，更新流体热物性参数；其中，各节点 是井筒每离散段的节点；
[0060] 7)判断井筒温度剖面是否收敛，若不收敛，则返回步骤6)进行迭代，直至井筒温 度剖面收敛，一般迭代3~5次；若收敛，则修正井筒的压力和温度；
[0061] 8)判断tn是否小于丨1"1，若小于则输出井温-井深曲线以及井筒压力对应的水合 物生成温度-井深曲线，并根据曲线判断是否有水合生成；反之若不小于，则返回步骤3)， 预测下一时刻的井筒压力、温度场模拟及水合物情况；
[0062] 其中，在计算第n+1个时间步长时，将用第n个时间步长的井筒温度场和压力场作 为初始条件。
[0063] 上述步骤4)中，稳态和非稳态下的井筒温度场模型建立方法如下：
[0064] 由于本发明针对试采阶段的长期稳定生产过程，因此采用解析方法建立稳态井筒 温度场模型，预测非线性井温剖面；针对试油测试的短期过程，建立非稳态井筒温度场模 型，预测不同测试制度和时间下的井温剖面。根据能量守恒定律、焦耳-汤普森效应、井筒 径向传热以及海上温度场影响，得到稳态井筒温度场模型和非稳态井筒温度场模型。
[0065] 4. 1)稳态井筒温度场模型
[0066] 地层段的流体温度乙为：
[0067]
[0068] 式中，z为井筒深度；步为倾斜校正系数。
[0069]
⑷
[0070] 海水段的流体温度Tf表达式：
[0071]
[0072] 式中，F。为焦耳-汤普森效应与动能作用系数；Tf(l为井筒入口温度；Tei(lS地层温 度；z为井筒深度；LK为松弛距离参数；9为井筒与水平面的夹角；g为重力加速度；步3为 组合参数；
[0073] 定义组合参数$s为：
[0074]
(6)
[0075] 4. 2)非稳态井筒温度场模型为：
[0076]
(7)
[0077] 式⑶~（7)中，a为温度随时间变化的衰减系数；Tei为任意深度处的原始地层 温度；Cp为油管流体比热，单位为J/kg-K;LK为松弛距离参数；L为垂直井筒总长度，单位为 m;ge为地温梯度，单位为K/m;gse为海水段温度梯度，单位为K/m;t为时间，单位为s。
[0078] 上述步骤4)中，多相流的井筒压力计算管流模型建立方法如下：
[0079]由于海上油气田一般是大斜度井或水平井，同时结合多种管流计算模型的适应性 条件，本发明优选Beggs和Brill方法作为井筒压力计算模型。本发明中Beggs和Brill 根据在长为15m、直径为25. 4mm和38. 1mm的管子中，用空气和水进行实验的基础上提出了 处理所有范围的多相流及任何角度的管流的普适化相关式。其中，Beggs和Brill多相管 流模型的基本方程、流型判别条件以及主要参数确定如下：
[0080] 根据机械能量守恒定律，气液两相气液两相管流的压力降消耗于三个方面：位差、 摩擦和加速度，即多相管流模型的基本方程为：
[0081]
[0082] 式中各项参数的解释如表1所示。 (8)
[0083] 表1式⑶中各项参数的解释表
[0084]
[0085] 其中，总压力梯度的计算步骤如下： dZ
[0086] (1)计算弗鲁德数NFr、无滑脱持液率匕和L^L2,L3,L4:
[0087]
[0088]
[0089] L:= 316E^ 302 ；
[0090] L2= 0. 0009252EL^ 4682；
[0091] L3= 0? 1El_14516;
[0092] L4= 0? 5EJ6.738;
[0093] 式中，无滑脱持液率；Q液体体积流量；
[0094] (2)根据步骤（1)中计算的基本参数，采用Beggs-Brill法判断流型：分离流、间 歇流或分散流。
[0095]当EL< 0.01时NFr<LpEl彡0.01时NFr<L2，则流型为离流；
[0096] 当El彡0? 01，L2<NFr<L3时，流型为过渡流；
[0097]当 0? 01 彡EL< 0? 4 时L3<NFr<L丨，EL彡 0? 4 时L3<NFr<L4，流型为间歇流；
[0098] 当EL<0?4时NFr彡L:，EL彡0?4时NFr>L4,流型为分散流。
[0099] (3)持液率 计算
[0100] 倾角为0的气液两相流动的持液率扎（9 )为：
[0101] HL(e) =Hl(0)It (9)
[0102] 式（9)中，扎(0)为同样流动参数下，水平流动时的持液率。
[0103]
(1〇)
[0104] 式（9)中，$为倾斜校正系1
对于垂直管：步 =1+0. 3C，C为系数，与无滑脱持液率、弗鲁德系数和液相速度有关；b、c均为流型常数；
[0105] (4)计算液相密度Pl:
[0106] pL=pLHL+Pg(l-HL) (11)
[0107] (5)确定流动阻力系数入：
[0108]
"2)
[0109] 式中：
[0110]
[0111]
[0112] 其中，h为在P、T下的液体粘度；yg为在P、T下的气体粘度，T为温度。
[0113] (6)计算井筒压力梯度：将步骤（1)~步骤（5)中计算的参数代入式（8)中计算 得到井筒压力梯度。
[0114] 上述各实施例仅用于说明本发明，各步骤都是可以有所变化的，在本发明技术方 案的基础上，凡根据本发明原理对个别步骤进行的改进和等同变换，均不应排除在本发明 的保护范围之外。
【主权项】
1. 一种深水测试井筒压力、温度场模拟及水合物预测方法，其包括以下步骤： 1) 准备参数：气藏参数、流体物性参数、井身结构参数和海水温度剖面等； 2) 对井筒进行离散化，并对井筒温度场和压力场进行初始化，初始化井筒温度等于环 境温度； 3) 利用气藏渗流数学模型求解地层流量和井底压力分布信息； 4) 将地层流量和井底压力分布信息输入井筒流动模型内进行迭代，进而求解出井筒 压力剖面，并判断井筒压力是否收敛，若收敛则进入步骤5)，反之继续判断井筒压力是否收 敛；其中，井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒温度场模型、多相流 的井筒压力计算管流模型； 5) 计算测试井内流体和环空流体在当前温度压力条件下的热物性参数，并计算地层段 和海水段的总传热系数； 6) 将井筒压力和总传热系数代入井筒温度场模型，从井底往井口递推计算各节点井筒 温度剖面，并修正地层段和海水段的环空温度差，更新流体热物性参数；各节点是井筒每离 散段的节点； 7) 判断井筒温度剖面是否收敛，若不收敛，则返回步骤6)进行迭代，直至井筒温度剖 面收敛；若收敛，则修正井筒的压力和温度； 8) 判断tn是否小于t n+1，若小于则输出井温-井深曲线以及井筒压力对应的水合物生 成温度-井深曲线，并根据曲线判断是否有水合生成；反之若不小于，则返回步骤3)，预测 下一时刻的井筒压力、温度场模拟及水合物情况。2. 如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述步骤3)中，所述气藏渗流数学模型为： 3. 1)无限大均质油气藏中心直井Laplace空间压力相应解为：式中，九ο为井底Laplace空间压力；u = sf(s)，其中s为Laplace空间时间变量，f (s) =1 ;&、K1分别为修正的零阶和一阶第二类Bessel函数；S为表皮系数；C D为无因次井筒 储集系数； 3. 2)无限大油藏斜井在Laplace空间压力相应解：式中，hfD为无因次射孔井段长度；s为Laplace空间时间变量；&为无因次半径；η为 级数系数；hD为无因次储层厚度；ζ D为井筒深度ζ的无因次值、为井段中心井筒深度ζ的无因次值。3.如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述步骤4)中，所述井筒流动模型包括海水温度场效应的稳态和非稳态下的井筒 温度场模型、多相流的井筒压力计算管流模型。4. 如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述稳态井筒温度场模型为： 地层段的流体温度乙为：式中，z为井筒深度；Φ为倾斜校正系数；海水段的流体温度Tf表达式：式中，F。为焦耳-汤普森效应与动能作用系数；Tftl为井筒入口温度；Teitl为地层温度；z 为井筒深度；LK为松弛距离参数；Θ为井筒与水平面的夹角；g为重力加速度；φ 3为组合参 数； 定义组合参数Φ 3为：5. 如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述非稳态井筒温度场模型为：式中，a为温度随时间变化的衰减系数；Tei为任意深度处的原始地层温度；CpS油管流 体比热，单位为J/kg-K ;Lk为松弛距离参数；L为垂直井筒总长度，单位为m ;ge为地温梯度， 单位为K/m ;gse为海水段温度梯度，单位为K/m ;t为时间，单位为s。6. 如权利要求3所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述多相流的井筒压力计算管流模型建立方法如下： 多相管流模型的基本方程为：/ (?Ρλ ((IP λ 式中，（t)为位差压力梯度；为摩擦压力梯度；^ 为加速度压力梯 dZ Jmm V似夕加速度 AP 度；为总压力梯度； dZ dP 其中，总压力梯度的计算步骤如下： CiZ (1) 计算弗鲁德数Nft、无滑脱持液率EjP L i，L2, L3, L4:式中，无滑脱持液率；Ql为液体体积流量； (2) 根据步骤（1)中计算的基本参数，采用Beggs-Brill法判断流型：分离流、间歇流 或分散流： 当EL< 0· Ol时N Fr< L p El彡0· Ol时N Fr< L 2，则流型为离流； 当0. 01，L 2< N L 3时，流型为过渡流； 当 0· 01 彡 EL< 0· 4 时 L 3< N Fr< L p El彡 0· 4 时 L 3< N Fr< L 4，流型为间歇流； 当EL< 0· 4时N Fr彡L i，0· 4时N Fr> L 4，流型为分散流； (3) 持液率计算 倾角为Θ的气液两相流动的持液率扎（0)为： HL( Θ )= Hl(O) Φ， 其中，扎(0)为同样流动参数下，水平流动时的持液率；式中，Φ为倾斜校正系数，对于垂直管：Φ = 1+0. 3C，C为系数，与无滑脱持液率、弗鲁德系数和液相速度有关；b、c均为流型常数； (4) 计算液相密度P ^l= f LhL+Ps(I-Hl); (5) 确定流动阻力系数λ :其中，μ ^3在P、T下的液体粘度；μ g为在P、T下的气体粘度；T为温度； (6) 计算井筒压力梯度：将步骤（1)~步骤（5)中计算的参数代入多相管流模型的基 本方程中计算得到井筒压力梯度。7.如权利要求1所述的一种深水测试井筒压力、温度场模拟及水合物预测方法，其特 征在于：所述步骤5)中，所述总传热系数采用热阻串联法计算。
【专利摘要】本发明涉及一种深水测试井筒压力、温度场模拟及水合物预测方法，其步骤：准备参数；对井筒温度场和压力场进行初始化；利用气藏渗流数学模型求解地层流量和井底压力分布信息；将地层流量和井底压力分布信息输入井筒流动模型内进行迭代，求解出井筒压力剖面并判断井筒压力是否收敛；计算热物性参数，并计算地层段和海水段的总传热系数；将井筒压力和总传热系数代入井筒温度场模型，计算各节点井筒温度剖面；判断井筒温度剖面是否收敛；判断tn是否小于tn+1，若小于则输出井温-井深曲线以及井筒压力对应的水合物生成温度-井深曲线，并根据曲线判断是否有水合生成；反之则返回预测下一时刻的井筒压力、温度场模拟及水合物情况。本发明可以广泛在石油勘探过程中应用。
【IPC分类】E21B49/00, E21B47/07, E21B47/06
【公开号】CN104895560
【申请号】CN201510333445
【发明人】何玉发, 蒋世全, 段永刚, 周建良, 李迅科, 杨秀夫, 陈伟, 方全堂
【申请人】中国海洋石油总公司, 中海油研究总院
【公开日】2015年9月9日
【申请日】2015年6月16日
最新回复(0)