Temperature Field Prediction Model for Multi-Layer Commingled Production Wellbore in Intelligent Wells and It's Application
-
摘要:
智能井多层合采过程中为了优选温度监测设备和确定测点位置，需要准确预测井筒温度剖面。根据智能井多层合采过程中的井筒内流体流动特征，考虑流体经过流量控制阀时，节流效应对井筒内流体流动参数的影响，建立了含流量控制阀的单油管多层合采井筒温度预测模型，并结合生产井的工况进行了数值模拟。模型预测结果表明，井筒温度随产层产出液性质、产液量、产层厚度、产层配比和地层温度梯度的变化呈规律性变化；与各产层单独开采相比，合采时的井筒温度高于各产层单独开采时的平均温度，且合采时的温度梯度最低。为了有效应用多层合采井筒温度场预测模型，基于流量控制阀处温度测量误差最小的原则，提出了温度传感器指标及测点的优选方法；基于井筒温度、温度梯度及流量控制阀处温降变化规律，提出了产层温度异常的解释方法。多层合采智能井井筒温度场预测模型为多层合采智能井温度测量装置的优选和温度变化规律的解释提供了理论依据。
Abstract:In order to optimize the temperature monitoring equipment and properly determine the positions of measuring points for the multi-layer commingled production of intelligent wells, it is necessary to accurately predict the whole wellbore temperature profile. According to the characteristics of fluid flow in the commingled production wellbore of intelligent wells, a prediction model for single-tubing multi-layer commingled production wellbore temperature with the flow control valve was established. Taking the throttling effect of flow control valve on fluids flow parameters in the system into account, a numerical simulation was carried out in combination with the working conditions of production wells. The model prediction results showed that the wellbore temperature changed regularly with several variables, including the produced fluids properties, liquid production rate, thickness of payzone, production allocation of each layer and geothermal gradient. Compared with the individual production layer, the wellbore temperature of commingled production wells is higher than the average temperature of wells with individually produced layer, and the temperature gradient of commingled production wellbore is the lowest. Based on the principle of minimum temperature measurement error at the flow control valve, an optimization method for temperature sensor index and measuring point was proposed. Taking into consideration the wellbore temperature, temperature gradient and the laws of temperature drop at flow control valve, it was possible to obtain an interpretation method of production layer temperature anomaly. The temperature field prediction model of multi-layer commingled production wellbore in intelligent wells provided a theoretical basis for optimizing the temperature measuring equipment and temperature data interpretation in multi-layer commingled production.
-
随着全球油气资源需求快速增长，难采难动用储量逐渐成为勘探开发目标，非均质、多产层油藏等复杂油藏的开发技术成为研究热点。智能井技术可以在单根油管中实现多层合采，生产前同时射开同井筒的多个储层，各层之间用封隔器封隔，并在各产层安装远程流量控制阀[1–3]，实现生产过程中各层段的开启和关闭而不停产，从而以顺序开采或者多层同采的方式开发多产层油藏[4–5]。
井下的长期监测数据在智能井获取油藏参数、监控油藏状况、制定开采计划和预测油井及油藏生产动态等过程中起到非常关键的作用[6–8]。国内多采用井下永久传感器、点式传感器和分布式温度传感器监测智能井的井筒温度[9–10]。H. J. Ramey Jr等人[11–12]在井筒流体温度预测方面进行了开创性的研究工作。20世纪90年代以来，C. S. Kabir等人[13–14]将流体温度看作是深度和循环时间的函数，提出了用于预测井内循环流体温度分布的方法。对于有外界流体流入的井筒内流体流动的问题，现有模型大多近似认为外来流体混合前为等温流动，由混合点处流体质量守恒和热量守恒得到混合后的温度[15]。但是，实际上外来流体在混合前通常要经过流量控制阀等节流装置，受节流效应的影响，通过节流装置时温度、压力会发生变化，尤其是含气量较高的油气井，节流效应十分显著，会导致温度、压力的计算误差较大[16]。
一般来说，多层合采油井的井筒温度剖面受各种因素的影响更为复杂，温度数据解释难度更大，目前这方面的研究仍然较少。因此，笔者建立了多层合采井筒温度场预测模型，分析了各种因素对井筒温度场的影响规律，并针对现场应用，提出了优选多层合采智能井温度测量装置和解释温度变化规律的方法。
1. 单油管多层合采井筒温度场模型
1.1 物理模型的建立
智能井多层合采是通过每个产层处的流量控制阀[17]控制该层流体的流入量，从而使各产层流体在同一根油管内混合并且采出的过程。为了提高多层合采智能井的井筒温度预测精度，在基本传热理论[18]和A.R.Hasan等人[19]研究成果的基础上，考虑流量控制阀的节流效应，将流动由等温流动改变为等熵流动，建立了预测井筒温度的物理模型（见图1）。
![]()
图 1 含流量控制阀的单油管多层合采井筒流动模型①为环空与环境传热；②为井筒与环空传热；Ta为环空内流体的温度，K；ΔTva为流量控制阀处的温度损失，K；Tt,u和Tt,d分别为混合点前、后油管内流体的温度，K。Figure 1. The flow model for single-tubing multi-layer commingled production wellbore with flow control valve流体从地层流到井口的过程分为3个阶段：1）流体从地层进入环空，在环空上下压差作用下流向流量控制阀入口；2）流体在节流压差作用下经流量控制阀流入井筒；3）流体在井筒上下压差作用下自下而上流向井口。
1.2 数学模型的建立
为了简化计算，首先进行以下假设：1）流体在井筒、环空内的流动为一维稳态流动，且同一截面上各点的温度、压力相等；2）油管到水泥环外缘间的传热为径向稳态传热，忽略油管壁和套管壁的热阻；3）井筒周围地层中的传热为非稳态，且不考虑沿井身方向的纵向传热；4）流体以紊流状态流经流量控制阀、无旋涡，流体充满控制阀并连续通过，流量不随时间变化；5）井筒及环空压力剖面、相组分已知。
1.2.1 流体从环空流向流量控制阀入口
流体在环空内为一维稳态流动，沿井深方向取环空流体微元，通过对流流入微元体的热量等于环空与地层和井筒交换的热量，其可以表述为[19]：
dqadz=qta−qF (1) 式中：qa为单位长度环空流体放出的热量，J/（m·s）；qta为环空与井筒交换的热量，J/（m2·s）；qF为环空与外界环境交换的热量，J/（m2·s）；z为井深，m。
利用稳态传热关系表示热流量qta和qF，则有：
(wCp)adTadz=2πrtUt(Ta−Tt)−2πrcUatDKe+rcUatD(Tei−Ta) (2) 式中：w为流体质量流量，kg/s；Cp为流体比热，J/（kg·K）；rt为油管外半径，m；Ut为井筒与环空的热交换系数，W/（m2·K）；Tt为井筒温度，K；rc为套管外半径，m；tD为无因次时间函数；Ua为井眼与外界环境的热交换系数，W/（m2·K）；Ke为环境导热系数，W/（m·K）；Tei为外界环境温度，K。
1.2.2 流体从流量控制阀流入井筒
在一定压力差下，通过改变流量控制阀节流口的流动阻力来控制节流口的流量，可以用喷嘴模型[20]模拟。流体在从流量控制阀入口流向出口的过程中没有热量输入，摩擦和耗散引起的热量损失可以忽略不计，流动过程可以看作等熵流动，流动过程中没有压力损失。
假设喷嘴上游、下游的流体密度相同，利用喷嘴流量系数CD对真实流体状态进行修正，可得：
CD=wacwis≈vacvis (3) 式中：wac为通过喷嘴的真实质量流量，kg/s；wis为通过喷嘴的等熵质量流量，kg/s；vac为通过喷嘴的真实流速，m/s；vis为通过喷嘴的等熵流速，m/s。
式（3）中的流速之比为流速系数，流速系数的平方与喷嘴效率系数ηN相等。
(vacvis)2=ηN=h0u−hd,ach0u−hd,is≈hu−hd,achu−hd,is (4) 式中：ηN为喷嘴效率系数；h0u为喷嘴上游流体的滞止焓，J/kg；hd,ac为喷嘴下游截面上流体的真实焓，J/kg；hd,is为喷嘴下游截面上流体的等熵焓，J/kg；hu为喷嘴上游截面上流体的焓，J/kg。
忽略动能项，滞止焓等于截面上流体的焓，于是得到：
$${C_{\rm{D}}}^2 \approx \frac{{{h_{\rm{u}}} - {h_{{\rm{d,ac}}}}}}{{{h_{\rm{u}}} - {h_{{\rm{d,is}}}}}} = \frac{{{C_{\rm{p}}}\left( {{T_{\rm{u}}} - {T_{{\rm{d,ac}}}}} \right)}}{{{C_{\rm{p}}}\left( {{T_{\rm{u}}} - {T_{{\rm{d,is}}}}} \right)}} = \frac{{{T_{\rm{u}}} - {T_{{\rm{d,ac}}}}}}{{{T_{\rm{u}}} - {T_{{\rm{d,is}}}}}}$ 式中： T u为喷嘴上游截面上流体的温度，K；T d,ac为喷嘴下游截面上流体的真实温度，K；T d,is为喷嘴下游截面上流体的等熵温度，K。$(5) 等熵温度变化的计算公式为：
(wCp)u(Tu−Td,is)=(pu−pd)(wCpKJT)u (6) 式中：pu和pd分别为喷嘴上游和下游流体的压力，MPa；KJT为Joule-Thomson系数。
联立式（5）和式（6），得到喷嘴上游和下游真实流体的温度差ΔTva：
ΔTva=Tu−Td,ac=CD2(Tu−Td,is)=CD2(pu−pd)(wCpKJT)u(wCp)u (7) 环空内流体经过流量控制阀后与井筒内已有流体混合，混合后控制阀下游流体的温度为：
Tt,d=(wCp)t,uTt,u+(wCp)a,u(Ta,u−ΔTva)(wCp)t,d (8) 式中：Ta,u为环空内喷嘴上游流体的温度，K；Tt,d为环空内喷嘴下游流体的温度，K；下标a和t分别表示流体在环空和井筒内的状态；下标u和d分别表示流体在喷嘴上下游的状态。
1.2.3 流体在井筒内自下而上流向井口
井筒内流动也可看作一维稳态流动，沿井深方向取井筒流体微元，通过对流流入微元体的热量等于井筒与环空交换的热量，能量守恒方程为：
dqtdz=−qta (9) 式中：qt为单位长度井筒内流体放出的热量，J/(m·s)。
利用稳态传热关系表示传热过程：
(wCp)tdTtdz=−2πrtUt(Ta−Tt) (10) 1.2.4 方程求解
联立式（2）与式（10），消去Ta，得：
d2Ttdz2=[2πrtUt(wCp)t+2πrtUt(wCp)a+1(wCp)a]dTtdz+2πrtUt(wCp)a(wCp)t2πrcUaTDKe+rcUaTD(Tei−Tt) (11) 假设环境温度函数为Tei=f(z)，井底的热量交换为0，则边界条件为：
Tt|z=zmax=Tei|z=zmax (12) dTtdz|z=zmax=0 (13) 式中：zmax为井底的深度，m。
考虑流量控制阀处的温度变化，连续条件为：
Tt|z=z+va=(wCp)t|z=z−vaTt|z=z−va+(wCp)a|z=z−va(Ta|z=z−va−ΔTva)(wCp)t|z=z+va (14) 式中：zva为流量控制阀所在处井深，m；ΔTva由式（7）计算；上标–、+分别表示流体混合前和混合后。
若环境温度函数为线性函数，则式（11）可以化简为二阶线性微分方程，采用数值方法和解析方法均可求解；若环境温度函数复杂，可以采用数值计算方法求解。模型求解时，取：
Tei=Te0+gGz (15) 采用差分法构造四阶精度的差分格式[21]，采用迭代法求解。
1.3 模型的验证
利用所建模型模拟文献[22]中多层合采智能井的井筒温度场，并与该井实测井筒温度进行对比，结果如图2所示。由图2可以看出，利用所建模型计算出的井筒温度与实测井筒温度的误差小于0.8 K，相对误差小于1.5%，表明所建模型的计算精度满足工程要求，可以用来预测多层合采智能井的井筒温度场。文献[22]中多层合采智能井的基本参数见表1（压力11.03～12.31 MPa）。
![]()
2. 井筒温度影响因素分析
多产层油藏合采时，合采井段的井筒温度变化复杂，更能体现各种因素对井筒温度的影响。首先，同一时间各产层的产液情况不同，不同时间相同产层的产液情况也不同，因此每个产层的产液状态都会影响井筒温度分布，如产出液性质、产液量、产层长度和产层配比等。此外，井筒温度还会受到环境温度、井身结构等因素的影响。为了分析各种因素对井筒温度的影响程度，以典型多层合采井M–X井为例，利用多层合采井筒温度预测模型，模拟不同参数下的井筒温度。
M–X井为三层合采直井，基本参数为：井深3 000.00 m；油管内径75.9 mm，外径88.9 mm；套管内径230.7 mm，外径244.5 mm；水泥环直径311.1 mm；下部产层处于2 950.00～3 000.00 m井段，中间产层处于2 900.00～2 950.00 m井段，上部产层处于2 850.00～2 900.00 m井段；油井总产液量150 m3/d，上中下各产层的产量配比为1∶1∶1，各产层产出状态相同，产量均为50 m3/d，产出液密度0.90 kg/L，产出液比热容2 200 J/（K·kg），产出液导热系数0.14 W/（m·K）；地温梯度5.0 K/100m，地层密度2.50 g/cm3，地层导热系数2.7 W/（m·K），地层比热容830 J/（K·kg）。
2.1 产层产出液性质
分别改变M–X 井中间产层产出液的密度、比热容和导热系数等流体性质参数（见表2），模拟不同产出液性质下的井筒温度剖面，结果见图3。
表 2 产出液的基本性质参数Table 2. Basic properties of output liquid序号 密度/（kg·m–3） 比热容/（J·K–1·kg–1） 导热系数/（W·m–1·K–1） A 900 2 200 0.14 B 1 000 2 200 0.14 C 900 2 640 0.14 D 900 2 200 0.17 ![]()
图 3 产出液的密度、比热容和导热系数对井筒温度的影响曲线Figure 3. Influence curves of the density, specific heat capacity and thermal conductivity of produced fluids on wellbore temperature由图3可知：1）中间产层产出液的性质发生改变，下部产层的井筒温度不变；2）中间产层产出液的密度和比热容增大，控制中间产层流量控制阀（井深2 950.00 m）处的温降变小；中间产层产出液的导热系数增大，控制中间产层流量控制阀处的温降变化不显著；3）中间产层产出液的密度和比热容增大，中间产层井段的井筒温度升高，温度梯度降低；中间产层产出液的导热系数增大，中间产层井段的井筒温度降低，温度梯度升高；4）中间产层产出液的密度和比热容增大，控制上部产层流量控制阀（井深2 900.00 m）处的温降变小；中间产层产出液的导热系数增大，控制上部产层流量控制阀（井深2 900.00 m）处的温降变大；5）中间产层产出液的性质改变，上部产层井段的井筒温度有所变化，但温度梯度无明显变化；6）产出液的比热容增大20%与导热系数增大20%相比，比热容增大对井筒温度的影响更为显著。
某产层的产出液性质改变，导致流量控制阀处温降变化，并影响该产层及以上井段的井筒温度。这是因为该产层的产出液经流量控制阀进入井筒，2种流体在短时间内快速混合产生温降；在相同地层温度、井身结构和产液量条件下，虽然导热传递的热量会因导热系数变化而改变，但总体来讲，相同环境条件和井身结构时，定量液体与地层交换的总热量趋于稳定；交换相同热量后，井筒温降随产出液性质变化而变化。随着产出液向上流动，比热容越大，放出相同热量后，温度降低越慢；而对于固定产层，产出液密度增大往往伴随着含水量增大或含气量减小，含水增大或含气减小都会导致产出液的比热容增大，传递相同热量的温度差变小；此外，密度增大引起的摩阻增大对井筒温度也有一定程度的影响。随着产出液继续流动，新产层的产出液进入井筒混合，这种变化逐渐变得不明显。
2.2 产层产液量
保持上部产层和下部产层的产液量不变，且不改 M–X 井的其他参数，只改变中间产层的产液量，模拟中间产层不同产液量下的井筒温度剖面，结果如图4所示。
![]()
图 4 中间产层产液量对井筒温度的影响Figure 4. Effect of Intermediate production layer liquid production rate on wellbore temperature由图4可知，中间产层产液量由40 m3/d增大至70 m3/d，流量控制阀处温降变小，该产层及以上井段的井筒温度都升高，温度梯度降低明显。这是因为在油管尺寸相同的条件下，产液量增大引起流量控制阀入口和油管内液体流速加快，热量交换时间变短，沿液体流动方向的温度变化较慢；新产层产出液进入井筒后，这种趋势仍然存在。
2.3 产层厚度
在井深、上部产层和下部产层的厚度与M–X 井相同的条件下，调整中间产层上边界的深度，从而改变中间产层厚度，模拟不同中间层厚度下的井筒温度剖面，结果如图5所示。
![]()
图 5 中间产层厚度对井筒温度的影响Figure 5. Effect of Intermediate production layer payzone thickness on wellbore temperature从图5可以看出，当总产液量不变而中间产层厚度增大后，流量控制阀处温降基本不变，但由于单位长度的产液量变小，油管内流体流速变慢，热交换更充分，该产层相同井深的井筒温度降低得更快，这种趋势一直延续到上部产层，整个合采井段井筒上下的温度差变大。
2.4 产层配比
保持油井总产液量150 m3/d不变，调整产层产量配比依次为1∶0∶0、1∶1∶0、1∶1∶1、1∶0∶1、0∶1∶0、0∶1∶1和0∶0∶1（1代表该层开采，0代表该层不开采），模拟各自的井筒温度剖面，结果如图6所示。
从图6可以看出：当产层配比依次为1∶0∶0、1∶1∶0、1∶1∶1、1∶0∶1、0∶1∶0、0∶1∶1和0∶0∶1时，井筒温度差呈由小到大的变化趋势；单独开采上、中、下产层时，井深2 850.00 m处的温度分别为390.50，390.97和392.90 K，平均为391.46 K，而3个产层合采时的该处温度为391.90 K，比单独开采时的平均温度高0.4%；单独开采上、中、下产层时，井深2 850.00 m处的井温梯度分别为0.015，0.020和0.017 K/m，3个产层合采时该井深处的温度梯度为0.012 K/m，比单独开采时低。较高的井内温度有利于抑制蜡质凝结等流动障碍的形成，因此定产量生产时采用产层合采的制度，可以简化油井生产措施，更有利于油井长期安全高效生产。
3. 井筒温度场预测模型的应用
3.1 智能井温度测量装置优选
智能井温度监测多采用光纤分布式温度传感器（DTS），与传统点式传感器相比，DTS无需在检测区域内来回移动，不会影响井内的温度平衡状态。但由于受传感器温度分辨率和空间分辨率的限制，温度测量存在一定误差，尤其是流量控制阀处的误差较大。因此，为了保证合采井段温度的监测精度，必须合理选择温度传感器类型、温度分辨率、空间分辨率及测点间距等。
3.1.1 流量控制阀处的温度测量误差分析
流量控制阀处存在温度突变，受分布式温度传感器测量特点影响，该处的测量误差较大。从测量方法上看，由于温度传感器空间分辨率有限，测点之间有间隔，只能以距流量控制阀最近的两相邻测点的温度差作为流量控制阀处的温降，这样就会存在一定的误差。从仪器精度上看，传感器自身的温度分辨能力有限，会使测量值偏离测点的真实值，再次加大了测量误差。
流量控制阀处的温度测量误差原因分析如图7所示。图7中，A(h1，t1)、B(h2，t2)为真实温度点，测量点A′(h1，t1′)、B′(h2，t2′)沿蓝色虚线随机分布；M(H，tM)、N(H，tN)为流量控制阀上、下游的真实温度点；δ为传感器的温度分辨率，d为传感器的采样间隔，Δh为测点A′与流量控制阀的间距；Δt为流量控制阀处的真实温降，Δt′为测点A′、B′的测量温度差。以测点温度差为流量控制阀处的温降，则流量控制阀处温降的相对测量误差为：
![]()
图 7 流量控制阀处的温度测量误差示意Figure 7. Schematic diagram of temperature measurement error at the flow control valveεr≈Δt′−ΔtΔt×100%=(t′2−t′1)−(tN−tM)tN−tM×100% (16) 某测点的最大绝对偏差为±δ，故|εr|最大的可能极限情况为：1）t1′=t1+δ，t2′=t2–δ，此时Δt′取到最小值，Δt′=t2–t1–2δ；2）t1′=t1-δ，t2′=t2+δ，此时Δt′取到最大值，Δt′=t2–t1+2δ。比较2种极限情况下的|εr|，其最大值即为该测点位置的最大相对误差|εr|max。
在算例井M–X井的温度剖面预测结果基础上，计算选取不同测点位置时流量控制阀处的测量误差。计算中选用的DTS–SR传感器空间分辨率为1.00 m，采样间隔0.50 m，温度分辨率0.01 K，假设预测得到的井筒温度等于真实井筒温度，计算结果见表3。
表 3 DTS–SR系统监测算例井温度实测值与所建立模型计算值的相对误差Table 3. The relative error between the measured temperature of example well with DTS-SR monitoring system and the calculated value by the established model测点间距Δh/m H=2 900 m相对误差εr，% H=2 950 m相对误差εr，% 平均|εr|max，% t1′=t1+0.01K
t2′=t2–0.01Kt1′=t1–0.01K
t2′=t2+0.01K|εr|max t1′=t1+0.01K
t2′=t2–0.01Kt1′=t1–0.01K
t2′=t2+0.01K|εr|max 0.10 –1.58 8.98 8.98 –1.60 9.96 9.96 9.47 0.20 –2.63 7.93 7.93 –2.78 8.78 8.78 8.35 0.30 –3.68 6.88 6.88 –3.96 7.60 7.60 7.24 0.40 –5.26 5.30 5.30 –5.73 5.83 5.83 5.56 3.1.2 基于误差分析的温度测量装置优选
以M–X井为例，对各种分布式温度测量系统进行优选，其过程如下：
1）预测M–X井的井筒温度剖面。
2）针对每种分布式温度测量系统，计算选取不同测点位置时每个流量控制阀处的εr及|εr|max，然后计算得到各个流量控制阀处|εr|max的平均值。
3）对比不同分布式温度测量系统在不同测点位置|εr|max的平均值，其最小值所对应的|εr|max即为该测量系统的最大测量误差，而对应的测点位置即为该测量系统的最优测点间距。如表3所示，当Δh=0.40 m时，M–X井流量控制阀处|εr|max的平均值最小，此时DTS–SR传感器的最大测量误差为5.56%。
4）比较各种测量系统的最大测量误差，优选误差较小的测量系统。
4种DTS传感器的技术性能指标及最大相对测量误差见表4。由表4可知，DTS–SR传感器的测量误差最小，因此选择DTS–SR传感器监测M–X井的井筒温度，其空间分辨率1.00 m，采样间隔0.50 m，温度分辨率0.01 K，最优测点应取为2 899.60，2 900.10，…，2 949.60，2 950.10 m，…。
表 4 各种DTS传感器的技术性能指标及监测M-X井时的最大相对测量误差Table 4. Technical performance indicators of various DTS systems and the maximum relative measurement error when monitoring Well M-X传感器 空间分辨率/m 采样间隔/m 温度分辨率/K 最大相对测量误差，% 最优测点间距Δh/m DTS-SR 1.00 0.50 0.01 5.83 0.40 DTS-FR 0.30 0.30 0.05 26.92 0.20 Sensa DTS 1.00 1.00 0.05 22.17 0.30 WFT-E10 2.00 1.00 0.10 33.21 0.90 3.2 温度测量数据解释
智能井井筒温度测量数据的解释是一个复杂的过滤、融合与反演过程，这需要一个庞大的数据解释处理平台和软件系统支持。但在数据处理之前，需要对井筒温度剖面进行定性分析。根据井筒温度剖面可以确定产出状态发生变化的位置，分析井筒温度预测结果（见表5）。
表 5 某井段的井筒温度剖面变化与产层状态/环境变化的关系Table 5. Relationship between wellbore temperature profile variation and production layer status/environmental changes in a well section产层状态/环境变化 温度是否变化 温度梯度是否变化 控制该井段的流量控制阀处温降是否变化 井段以上的产层状态改变 否 否 否 井段所在的产层状态改变 是 是 是/否 井段以下的产层状态改变 是 是/否 是/否 井段以上的环境改变 是 否 否 井段所在层的环境改变 是 是 是 井段以下的环境改变 是 否 是 若某井段井筒温度有变化，首先根据表5判断环境温度梯度是否发生变化，排除其干扰；然后自下而上判断其他产层的产出状态是否发生了变化，即可确定变化产层的位置。一般来说，某产层的产出状态发生变化后，该井段的井筒温度变化，并且井筒温度梯度变化明显；而某井段只有温度变化显著但温度梯度基本不变时，该井段产层的产出状态可能没有变化，其温度变化可能是该井段以下的其他产层的产出状态发生变化导致的。
根据各种因素对多层合采温度场分布的影响规律，可以对井筒温度剖面作出定性解释，分析油井产出状态的变化，从而发现水侵、气侵、砂堵等的早期特征，及时调整或关闭产层、采取防砂措施等。井筒温度剖面的定性解释结果见表6，若已经确定为某产层状态变化，可参照表6分析该产层所在井段的温度剖面，从而确定可能变化的参数及可能出现的油井异常情况。例如，某井段的温度剖面显示温度升高、温度梯度变小且该层流量控制阀处温降变小，已确定地层温度梯度等环境条件未改变，则优先考虑该层或该层以下产层的产量增大，再考虑该产层是否发生水侵。
表 6 多层合采井筒温度剖面的定性解释Table 6. Qualitative interpretation on the wellbore temperature profile in multi-layer commingled production wells温度变化情况 可能的参数变化情况及解释结果 优先等级 温度升高 产出液密度增大，考虑水侵 2 产出液比热容增大，考虑水侵 2 产出液导热系数减小，考虑气侵 4 产层厚度减小，考虑砂堵 3 该层或该层以下的产层产量增大 1 地层温度梯度增大 1 温度梯度减小 产出液密度增大，考虑水侵 2 产出液比热容增大，考虑水侵 2 产出液导热系数减小，考虑气侵 4 产层厚度减小，考虑砂堵 3 该层或该层以下的产层产量增大 1 地层温度梯度减小 1 流量控制阀处温差增大 产出液密度减小，考虑气侵 2 产出液比热容减小，考虑气侵 2 该层或该层以下的产层产量减小 1 地层温度梯度增大 1 注：优先等级1、2、3和4代表优先次序由高到低。 4. 结 论
1）考虑外来流体混合前通过节流装置时，节流效应对流体流动参数的影响，将控制阀处的流动看作等熵节流流动，建立了多层合采智能井井筒温度剖面预测模型。
2）根据流动特征，将含流量控制阀的单油管多层合采流体的流动划分为环空内流动、流量控制阀节流和井筒内的流动等3个阶段。所建立模型的预测结果与文献中实测数据的绝对误差在0.8 K以内，相对误差小于1.5%。
3）模拟分析了产出液性质、产液量、产层厚度、产层配比和地层温度梯度等参数对井筒温度分布的影响，结果表明，三层合采时的井筒温度比各层分别配产时高0.4%左右，且三层合采时的井筒温度梯度最低，比单独开采时的温度梯度低40%。
4）根据模型模拟结果，以流量控制阀处的温度测量误差最小为原则，提出了温度传感器指标及测点的优选方法；以井筒温度变化规律为基础，提出了异常产层位置和异常参数变化的确定方法和温度解释方法。
-
![]()
图 1 含流量控制阀的单油管多层合采井筒流动模型
①为环空与环境传热；②为井筒与环空传热；Ta为环空内流体的温度，K；ΔTva为流量控制阀处的温度损失，K；Tt,u和Tt,d分别为混合点前、后油管内流体的温度，K。
Figure 1. The flow model for single-tubing multi-layer commingled production wellbore with flow control valve
![]()
![]()
图 3 产出液的密度、比热容和导热系数对井筒温度的影响曲线
Figure 3. Influence curves of the density, specific heat capacity and thermal conductivity of produced fluids on wellbore temperature
![]()
图 4 中间产层产液量对井筒温度的影响
Figure 4. Effect of Intermediate production layer liquid production rate on wellbore temperature
![]()
图 5 中间产层厚度对井筒温度的影响
Figure 5. Effect of Intermediate production layer payzone thickness on wellbore temperature
![]()
图 7 流量控制阀处的温度测量误差示意
Figure 7. Schematic diagram of temperature measurement error at the flow control valve
表 1 文献[22]多层合采智能井基本参数
Table 1 Basic parameters of the multi-layer commingled production intelligent wells in the reference[22]
层段 深度/m 产液量/（m3·d–1） 上部油层 1 646.00～1 737.00 91.0 中部油层 1 737.00～1 859.00 161.7 下部油层 1 859.00～1 890.00 163.7 
下载: 导出CSV
表 2 产出液的基本性质参数
Table 2 Basic properties of output liquid
序号 密度/（kg·m–3） 比热容/（J·K–1·kg–1） 导热系数/（W·m–1·K–1） A 900 2 200 0.14 B 1 000 2 200 0.14 C 900 2 640 0.14 D 900 2 200 0.17
下载: 导出CSV
表 3 DTS–SR系统监测算例井温度实测值与所建立模型计算值的相对误差
Table 3 The relative error between the measured temperature of example well with DTS-SR monitoring system and the calculated value by the established model
测点间距Δh/m H=2 900 m相对误差εr，% H=2 950 m相对误差εr，% 平均|εr|max，% t1′=t1+0.01K
t2′=t2–0.01Kt1′=t1–0.01K
t2′=t2+0.01K|εr|max t1′=t1+0.01K
t2′=t2–0.01Kt1′=t1–0.01K
t2′=t2+0.01K|εr|max 0.10 –1.58 8.98 8.98 –1.60 9.96 9.96 9.47 0.20 –2.63 7.93 7.93 –2.78 8.78 8.78 8.35 0.30 –3.68 6.88 6.88 –3.96 7.60 7.60 7.24 0.40 –5.26 5.30 5.30 –5.73 5.83 5.83 5.56
下载: 导出CSV
表 4 各种DTS传感器的技术性能指标及监测M-X井时的最大相对测量误差
Table 4 Technical performance indicators of various DTS systems and the maximum relative measurement error when monitoring Well M-X
传感器 空间分辨率/m 采样间隔/m 温度分辨率/K 最大相对测量误差，% 最优测点间距Δh/m DTS-SR 1.00 0.50 0.01 5.83 0.40 DTS-FR 0.30 0.30 0.05 26.92 0.20 Sensa DTS 1.00 1.00 0.05 22.17 0.30 WFT-E10 2.00 1.00 0.10 33.21 0.90
下载: 导出CSV
表 5 某井段的井筒温度剖面变化与产层状态/环境变化的关系
Table 5 Relationship between wellbore temperature profile variation and production layer status/environmental changes in a well section
产层状态/环境变化 温度是否变化 温度梯度是否变化 控制该井段的流量控制阀处温降是否变化 井段以上的产层状态改变 否 否 否 井段所在的产层状态改变 是 是 是/否 井段以下的产层状态改变 是 是/否 是/否 井段以上的环境改变 是 否 否 井段所在层的环境改变 是 是 是 井段以下的环境改变 是 否 是
下载: 导出CSV
表 6 多层合采井筒温度剖面的定性解释
Table 6 Qualitative interpretation on the wellbore temperature profile in multi-layer commingled production wells
温度变化情况 可能的参数变化情况及解释结果 优先等级 温度升高 产出液密度增大，考虑水侵 2 产出液比热容增大，考虑水侵 2 产出液导热系数减小，考虑气侵 4 产层厚度减小，考虑砂堵 3 该层或该层以下的产层产量增大 1 地层温度梯度增大 1 温度梯度减小 产出液密度增大，考虑水侵 2 产出液比热容增大，考虑水侵 2 产出液导热系数减小，考虑气侵 4 产层厚度减小，考虑砂堵 3 该层或该层以下的产层产量增大 1 地层温度梯度减小 1 流量控制阀处温差增大 产出液密度减小，考虑气侵 2 产出液比热容减小，考虑气侵 2 该层或该层以下的产层产量减小 1 地层温度梯度增大 1 注：优先等级1、2、3和4代表优先次序由高到低。
下载: 导出CSV
-
[1] 单彦魁, 韦红术, 张俊斌, 等. 南海深水油田智能井修井挑战及对策[J]. 钻采工艺, 2017, 40(6): 63–65. doi: 10.3969/J.ISSN.1006-768X.2017.06.19 SHAN Yankui, WEI Hongshu, ZHANG Junbin, et al. Challenges and countermeasures for workover on deepwater smart well in South China Sea[J]. Drilling & Production Technology, 2017, 40(6): 63–65. doi: 10.3969/J.ISSN.1006-768X.2017.06.19
[2] 石军太, 李骞, 张磊, 等. 多层合采气井产能指示曲线异常的原因与校正方法[J]. 天然气工业, 2018, 38(3): 50–59. doi: 10.3787/j.issn.1000-0976.2018.03.006 SHI Juntai, LI Qian, ZHANG Lei, et al. An abnormality of productivity indicative curves for multi-layer gas wells: reason analysis and a correction method[J]. Natural Gas Industry, 2018, 38(3): 50–59. doi: 10.3787/j.issn.1000-0976.2018.03.006
[3] 何祖清, 梁承春, 彭汉修, 等. 鄂尔多斯盆地南部致密油藏水平井智能分采技术研究与试验[J]. 石油钻探技术, 2017, 45(3): 88–94. HE Zuqing, LIANG Chengchun, PENG Hanxiu, et al. Research and tests on horizontal well smart layering exploiting technology in tight oil reservoirs in Southern Ordos Basin[J]. Petroleum Drilling Tech-niques, 2017, 45(3): 88–94.
[4] 石艺. 国际石油2015年十大科技进展(四):全电动智能井系统取得重大进展[J]. 石油钻采工艺, 2016, 38(1): 92. SHI Yi. Top ten scientific and technological progress of inter-national petroleum in 2015(4): major progress of all-electric smart well system[J]. Oil Drilling & Production Technology, 2016, 38(1): 92.
[5] 王浩. 智能完井生产优化方法研究[D]. 西安: 西安石油大学, 2017. WANG Hao. Research on production optimization method of intelligent completion[D]. Xi’an: Xi’an Shiyou University, 2017.
[6] 姚军, 刘钧荣, 张凯. 国外智能井技术[M]. 北京: 石油工业出版社, 2011: 21-55. YAO Jun, LIU Junrong, ZHANG Kai. Intelligent well technology abroad[M]. Beijing: Petroleum Industry Press, 2011: 21-55.
[7] TIWARI A, BROWN W P. System and method for simultaneous visualization of fluid flow within well completions and a reservoir: US13811826[P]. 2013-07-11.
[8] 刘涛, 刘生平. 智能分层控采技术在牛263井的应用[J]. 化工管理, 2016, 29(16): 157. doi: 10.3969/j.issn.1008-4800.2016.16.084 LIU Tao, LIU Shengping. Application of intelligent stratification and production control technology in Niu 263 Well[J]. Chemical Enterprise Management, 2016, 29(16): 157. doi: 10.3969/j.issn.1008-4800.2016.16.084
[9] 张振奇, 童茂松, 曹庆芳, 等. 光纤传感器及其在石油钻探中的应用[J]. 石油仪器, 2007, 21(4): 1–6. doi: 10.3969/j.issn.1004-9134.2007.04.001 ZHANG Zhenqi, TONG Maosong, CAO Qingfang, et al. Optic fiber sensor and its application in drilling and exploration[J]. Petroleum Instruments, 2007, 21(4): 1–6. doi: 10.3969/j.issn.1004-9134.2007.04.001
[10] DA SILVA M F, MURADOV K M, DAVIES D R. Review, analysis and comparison of intelligent well monitoring systems[R]. SPE 150195, 2012.
[11] RAMEY H J Jr. Wellbore heat transmission[J]. Journal of Petroleum Technology, 1962, 14(4): 427–435. doi: 10.2118/96-PA
[12] WILLHITE G P. Over-all heat transfer coefficients in steam and hot water injection wells[J]. Journal of Petroleum Technology, 1967, 19(5): 607–615. doi: 10.2118/1449-PA
[13] KABIR C S, HASAN A R, JORDAN D L, et al. A well-bore/reservoir simulator for testing gas wells in high-temperature reservoirs[J]. SPE Formation Evaluation, 1996, 11(2): 128–134. doi: 10.2118/28402-PA
[14] HASAN A R, KABIR C S, AMEEN M M, et al. A fluid circulating temperature model for workover operations[J]. SPE Journal, 1996, 1(2): 133–144. doi: 10.2118/27848-PA
[15] YIN Bangtang, LI Xiangfang, LIU Gang. A mechanistic model of heat transfer for gas-liquid flow in vertical wellbore annuli[J]. Petroleum Science, 2018, 15(1): 135–145. doi: 10.1007/s12182-017-0193-y
[16] 廖成龙, 黄鹏, 李明, 等. 智能完井用井下液控多级流量控制阀研究[J]. 石油机械, 2016, 44(12): 32–37. LIAO Chenglong, HUANG Peng, LI Ming, et al. Downhole hydraulic operated multistage flow control valve for intelligent well completion system[J]. China Petroleum Machinery, 2016, 44(12): 32–37.
[17] AJAYI A A, KONOPCZYNSKI M R. Theory and application of probabilistic method of designing customized interval control valves choke trim for multizone intelligent well systems[R]. SPE 110600, 2007.
[18] HOLMAN J P. Heat transfer[M]. New York: McGraw-Hill Book Company, 1976.
[19] HASAN A R, KABIR C S. Fluid flow and heat transfer in wellbores[M]. 2nd ed. Richardson: Society of Petroleum Engineers, 2018.
[20] CENGEL Y A, BOLES M A. Thermodynamics: an engineering approach[M]. 3rd ed. Hightstown: WCB/McGraw-Hill, 1998.
[21] 李维国, 同登科. 数值计算方法[M]. 2版. 东营: 中国石油大学出版社, 2009: 243-249. LI Weiguo, TONG Dengke. Numerical computation method[M]. 2nd ed. Dongying: China University of Petroleum Press, 2009: 243-249.
[22] MURADOV K M, DAVIES D R. Zonal rate allocation in intelligent wells[R]. SPE 121055, 2009.
-
期刊类型引用(4)
1. 冯高城，尹彦君，马良帅，张亮，王伟. 海上油田智能井技术发展应用及探讨. 西南石油大学学报(自然科学版). 2022(04): 153-164 . 
百度学术
2. 赵向阳，赵聪，王鹏，梁晓阳，杨谋. 超深井井筒温度数值模型与解析模型计算精度对比研究. 石油钻探技术. 2022(04): 69-75 . 本站查看
3. 卫胜利，王亚强，高致斌，张刚刚. 低渗砂岩油藏多层合采层间干扰因素分析. 中外企业家. 2020(20): 230 . 百度学术
4. 毛新军，曹植纲，陈超峰，胡广文，封猛，咸玉席. 高产井地层特征的热流耦合温度反演新方法. 石油钻探技术. 2020(04): 118-123 . 本站查看
其他类型引用(4)
计量
- 文章访问数: 1056
- HTML全文浏览量: 693
- PDF下载量: 77
- 被引次数: 8
