简介面向油气田开发工程领域研究人员与现场工程师该博士论文系统阐述了基于分布式光纤温度测试DTS的水平井产出剖面解释理论与应用。论文从油藏渗流和井筒流动机理出发考虑微热效应、热膨胀、热传导、热对流及粘性耗散和地层伤害建立多相流水平井温度预测复合模型并基于L-M与MCMC反演构建正反演框架实现对温度/压力剖面最佳拟合。通过底水突破、油水同层等不同出水形式和位置及渗透率变化案例验证表明MCMC反演在定量解析产水位置方面更优。资源为PDF单文件约20.69MB适合从事水平井生产测井、找水堵水的技术人员参考目前已有63人学习可用于指导实际水平井出水位置定量解释与生产动态分析。 做水平井产出剖面解释的人十有八九都被传统生产测井折腾过水平段那么长仪器送不下去即便靠爬行器进去了多相流流型乱得一塌糊涂低产层段谁在出液、出多少根本说不清楚。我后来把分布式光纤温度测试DTS引进解释流程才真正把这个问题从“凭经验猜”变成了“用数据反演”。DTS能沿井筒连续测出温度剖面配合水平井产出剖面解释理论就能把各段产液贡献反演出来。这篇文章我把整套方法、建模逻辑、实操流程和踩过的坑一次说透适合生产测井解释、油藏动态监测、水平井管理的工程师参考。1. 为什么水平井产出剖面偏偏要用DTS1.1 常规测井在水平井面前有多无力传统产液剖面测试主要靠生产测井仪器PLT下井测量涡轮流量、持水率、密度、压力等参数。这套方法在直井里成熟可靠一到水平井就卡壳。首先是仪器输送问题水平段无法靠电缆自重下放必须用爬行器爬行器本身有遇卡风险作业周期长一趟测井成本动辄几十万还经常拿不到合格数据。其次是测量环境问题。水平段倾斜角大气液两相常常出现层状流、段塞流、搅拌流交替流型极不稳定。涡轮流量计在这种流型下转速忽高忽低持水率计测出来的数值也很难和产量建立稳定关系。更麻烦的是很多水平井产液量不高单段产液量很低时涡轮可能根本转不动仪器响应进入“死区”剖面上某个层段是否有流体产出只能靠猜。再者PLT本质是点测仪器从井底往上提拉测到的是某几个深度点的快照数据。水平段动辄几百上千米点测根本无法交代清楚每一段的产出贡献。这就倒逼我们换一种能覆盖全井筒、不干扰流场、还能实时监测的技术路线DTS恰好满足了这些需求。1.2 DTS到底测了什么怎么变成温度剖面DTS的原理不复杂向井下光纤发射激光脉冲光纤分子与光子发生非弹性碰撞产生拉曼散射光。其中反斯托克斯光对温度敏感斯托克斯光对温度相对不敏感通过计算两者光强比值就能解调出光纤上每个位置的温度。光纤本身就是传感器不用通电、不用动井筒结构空间分辨率可以做到1米左右温度分辨率能做到0.1℃甚至更高单次采集就能得到一条沿井深的连续温度剖面。在实际应用中DTS光纤有两种典型布设方式。一种是永久安装在固井时把光纤绑在套管外壁或油管外壁随井寿命长期监测另一种是临时下入用连续油管携带光纤下到目标井段测完再起出适合探井或短期监测。布设方式不同后续解释时热边界条件不一样——光纤贴在套管外和悬在油管内感受到的温度延迟和幅值差异是不同的解释模型必须区分处理。DTS真正厉害的地方在于它产出的是温度-深度-时间的三维数据体。地层流体在哪一段进入井筒会把井筒温度拉偏多少这个信号随产量变化、随生产时间演化。我们要做的就是从这张“温度指纹”里反推各段产量这就是产出剖面解释理论的核心任务。2. 解释理论温度剖面到产量剖面的换算逻辑2.1 三个物理机制决定了温度剖面长什么样温度剖面不是在地层里量体温它是多种物理过程叠加的结果。要想从温度反演产量先得搞明白温度为什么会被产量改变。第一个机制是焦耳-汤姆逊效应J-T效应。流体从储层流向井筒压力降低、节流膨胀温度随之变化。不同流体J-T系数差异很大天然气J-T系数典型为正值且较大节流后明显降温可以降几度甚至十几度水J-T系数接近零低温时可能轻微降温、高温时轻微升温轻质油通常降温重油有时反而升温。因此产气段温度会出现明显低温异常产水段温度变化微弱这个差异是判断流体性质的重要线索。第二个机制是井筒-地层换热。流体进入井筒后沿水平段流动不断和周围地层交换热量。热传导方向取决于流体温度和地温的温差地层温度随深度递增水平段各点的围岩温度本底不同产量大的层段新鲜流体刚进入井筒、温度还未充分平衡和围岩温差大换热强烈产量小的层段流体在井筒里滞留时间长温度早被周围地层“同化”。所以温度剖面上的“冷点”或“热点”位置、幅度和产量分布有直接定量关系。第三个机制是地温梯度本身。水平井在造斜段和水平段穿过的地层深度不同即便没有流体产出静止井筒内温度也会沿深度变化。解释时必须先建立一条“零产出基线”把所有温度异常都和这条基线做差才能把产量引起的温度信号剥离出来。我见过不少解释翻车的案例都是没把基线做准把地温梯度当成产层信号结果反演出一堆根本不存在的“高产层”。2.2 从正演到反演先会算温度才能反推产量产出剖面解释在数学上是个典型的反问题。解法也不是硬刚而是先走正演再走反演。正演很简单假设我知道各段产液量那我就能用传热模型算出井筒温度剖面。反演就是反过来已知实测温度剖面倒推各段产量。正演模型通常把水平段划分成若干小段每一段赋予一个产液量q_i流体从这一段进入井筒后与原有流体混合沿井筒向跟端流动。井筒流体温度变化由能量守恒方程控制流体携带的热量、与地层的径向导热、流体节流引起的焦耳-汤姆逊温变三者的平衡决定了下一段的温度。地层的径向导热可以采用类似Ramey模型的瞬态传热解考虑水泥环、套管和地层的热阻。所有这些方程联立给定各段q_i就能算出理论温度剖面T_cal(z)。反演过程则是优化的思路。我先给各段q_i一个初始猜测正演出T_cal(z)和实测T_obs(z)对比计算残差。然后不断调整q_i让残差最小化。这里有个关键问题温度信号对产量分布不是一一对应的两个完全不同的产量分布可能算出几乎相同的温度剖面也就是反演高度不适定。所以目标函数不能只有拟合误差还得加正则项比如要求相邻段产量变化平滑或者产量值必须在某个合理区间内。正则化系数的选取很考验经验太小结果震荡太大结果被抹平。反演算法我习惯先用粒子群或差分进化这类全局优化算法粗扫一遍找到较好初值再用梯度类算法精调。纯梯度法很容易陷入局部极小尤其是水平段分段多、参数维度高的时候纯全局算法收敛又太慢。两者结合效率和稳定性都能兼顾。3. 实操过程一口典型水平井的完整解释流程3.1 数据准备与预处理我完整走过一遍从原始DTS数据到最终产出剖面的流程这里以一口水平段长度约800米、6段压裂的水平井为例。实测数据是光纤在套管内下放、连续监测得到的T(z,t)数据体生产数据是井口油压、日产液量、含水率。拿到原始数据后第一件事不是建模而是预处理。DTS原始温度存在几个常见问题激光脉冲衰减引起的深度越深信号越弱、部分深度点数据跳变、光纤接头位置温度异常。我的预处理顺序固定为坏点剔除、数据平滑、深度对齐、温度标定。坏点用滑动窗口中位数滤波处理阈值超过3倍标准差就剔除。平滑用窗口宽度5米的中值滤波既能去噪又不至于把1米量级的局部热点磨掉。深度对齐很重要因为光纤下放过程中可能伸长同一口井不同时期测得的深度坐标要对齐到统一的参考深度否则时间维度上的对比全乱套。温度标定是最容易忽略的一步。DTS解调出的温度是相对值不同通道、不同时间的基线会漂移。我的做法是找一个稳定的参考段比如造斜点以上的无产出井段把该段实测温度校正到已知的静止地温可由井温测试或地温梯度推算再用这个偏移量校正全井段。这一步不做后面所有温度异常的绝对值都是失真的J-T效应的幅度判断也会出问题。3.2 模型参数与解释段划分预处理完成后进入参数准备阶段。产出剖面解释需要的参数分四类井筒与完井参数、地层热物性参数、流体物性参数、生产边界参数。我把常用的参数列个表方便大家对照准备参数类别关键参数获取方式备注井筒与完井井筒内径、套管外径、水泥环厚度、完井方式完井设计文件分段压裂井需拿到射孔/滑套位置地层热物性地温梯度、地层导热系数、地层热扩散系数区域地温资料、岩心热物性测试导热系数影响换热强度敏感度高流体物性油气水J-T系数、比热容、密度、黏度PVT报告、相态软件计算多相流时按体积加权生产边界井口油压、日产液量、含水率、井底流压生产日报、压力计数据用来约束反演总量解释段划分的原则是“先按完井段再按信号特征细分”。以这口井为例6段压裂滑套位置明确我先把水平段按滑套位置分成6大段。然后看温度曲线上是否有明显的局部低温尖峰——如果有说明该位置附近流体进入井筒后J-T效应显著会在这个尖峰附近细化一个或多个子段。最终我把水平段划成10个解释段每段长度60到100米不等。3.3 正演-反演迭代与结果校验建模时我先算“零产出基线”假设整个水平段不产液根据地温梯度和井筒静止传热算出理论温度剖面作为所有温度异常的参考面。然后以日产液量75立方米为总量约束初始化各解释段产量。第一次正演算完拟合残差很大这是正常的重点是看各段对温度的贡献方向是否合理。迭代过程中我盯三个指标残差均方根误差RMSE、各段产量是否在合理区间、相邻段产量是否突变。正常两三小时能收敛。这口井最终反演结果和压裂设计规模高度吻合——第2段和第5段压裂加砂量大、裂缝半长短反演产液贡献分别占22%和26%最低的是第4段只有7%和微地震监测显示的改造弱区对应上了。为了验证可靠性我们还和井口分层流量计数据做了总量对比误差在5%以内。结果校验是解释报告里必须做的一件事。光有反演结果没有验证手段我自己心里都没底。我常用的校验方式有三种总量约束校验所有段产量加起来必须等于井口实测产液量点测对照校验如果有PLT点测或示踪剂测试数据直接按深度段对比生产时间演化校验DTS是在线监测的解释同一口井不同时间的温度剖面产量分配随时间的变化趋势要符合生产动态逻辑比如某段地层压力亏空后贡献逐渐下降这是很好的物理一致性检查。4. 常见问题与排查技巧实录4.1 温度数据本身带来的坑所有解释的“锅”最后都能甩到数据质量上但现实是DTS数据往往没你想的那么干净。最常见的坑是温度漂移。光纤在井下长期工作受氢黑效应、辐射衰减影响反斯托克斯信号会衰减解调温度可能整体偏移而且不同深度的漂移量还不一样。排查方法是看长时间稳定监测曲线若温度基线缓慢整体平移就是漂移用我前面说的参考段标定法重新校正。第二个坑是空间分辨率限制。DTS的空间分辨率通常1米这意味着极短产层段的温度异常会被“抹平”。我遇到过一口井的单段压裂簇间距只有2米温度剖面上两个簇的低温异常连成一个宽平台分都分不开。解决办法是引入完井信息做先验约束既然知道滑套位置就强制每段产量只在滑套位置附近才能非零不要指望纯靠数据把子簇分出来。第三个坑是光纤局部失效。高温、流体冲刷、光纤弯折都可能导致某段信号异常跳动或彻底丢失。我的排查技巧是画彩图看时间演化——单条曲线不好看出问题但T(z,t)彩图上断点的位置通常表现为固定深度上的条带状异常一眼就能识别果断把这段从解释域里剔除。4.2 解释模型参数里隐藏的坑模型参数里最要命的是地温梯度不准确。DTS解释对地温梯度极其敏感地温梯度差0.5℃/100米反演的产量分布可能面目全非。水平井造斜段长、水平段横穿的地层深度变化大用区域平均地温梯度往往不够我建议优先用实测静温数据生产前停井足够长时间后的DTS静态温度剖面就是最真实的地温基线比任何经验公式都可靠。J-T系数用错也会引发系统性偏差。油、气、水的J-T系数差异很大如果井筒里是气液两相流用单一流体参数计算就会把产气段的降温幅度低估或高估。我踩过的坑是某井油压下降后部分层段流体发生相态变化反演结果从“高产液”突然变成“高含气”后来查了PVT报告才发现J-T系数没按相态切换。现在我的标准做法是先根据井口压力和流体相态模拟软件算出各段压力下的J-T系数场而不是用地面条件下的常数。还有井筒积液问题。水平井低产期井筒下部容易积液积液段的温度动态和产液段完全不同会把解释结果带偏。判断方法很简单——无产出阶段温度剖面是否随时间抬升如果某段温度持续向地层温度逼近说明液体在静止换热这段很可能积液或封隔不佳解释时要么单独处理要么剔除。4.3 反演结果不唯一怎么办反演最容易被挑战的弱点是多解性。不同产量分布得到几乎相同温度剖面这在数学上是常态。我的处理原则是“框架约束不确定性量化”。框架约束包括总量守恒约束各段产量加和必须等于井口实测值空间连续性正则化相邻段产量不会突变物性界限约束每段产液量不能为负、不能超过单段最大产能。不确定性量化的手段是蒙特卡洛。我对每个解释段设置合理的参数扰动范围比如地温梯度扰动±0.3℃/100米、导热系数扰动±10%、J-T系数扰动±5%批量跑几百次反演得到每个解释段产量的概率分布而不是单一数值。解释报告里我现在不只写“第3段产液贡献22%”而会写“P50贡献22%P10-P90区间为17%-28%”。这种表达在工程决策里更有价值也更能应对质疑。5. 应用场景可以做哪些扩展这套方法并不局限于油井水平井产出剖面。页岩气压裂水平井是目前应用最热的方向之一压裂段数多、每段产量差异大DTS温度解释能评估各压裂段的有效改造长度和产气贡献。我在实际项目里拿DTS反演的段产气量和分布式声波传感DAS的产气剖面做交叉验证两者趋势一致性很好但DTS对低产气段的响应更稳定DAS则对流型变化更敏感结合使用效果明显好于单一手段。注水井的吸水剖面也可以借鉴这套框架。把“产出流体进入井筒”换成“注入水从井筒进入地层”J-T效应变成了注入水与地层的温差主导同样是用温度剖面反演各段吸水量。注水剖面最怕水窜DTS能实时捕捉水窜通道的温度异常变化比定期下流量计及时得多。此外DTS温度数据体还能做生产动态的连续监测。同一口井不同时间的温度剖面解释结果能反映各段产量的递减规律和地层压力变化为后续调层、堵水、补孔提供依据。我目前正在做的方向是把深度学习的代理模型和正演模型结合起来用大量正演样本训练神经网络替代耗时的数值反演让解释从小时级降到分钟级实现边采集边解释。这个方向还在完善中成熟后我会把可复现的案例再整理出来。最后分享一点个人体会DTS这套解释方法难的不是数学也不是软件的某一步而是每个环节里“看似无关的参数”之间强烈的耦合关系。地温梯度偏一点J-T系数偏一点解释结果可能南辕北辙。所以我现在养成一个习惯每次解释第一版结果出来先不急着调反演算法而是回过头重新核数据标定和基线计算往往问题就出在那里。如果大家在复现中遇到问题欢迎直接联系我交流尤其是参数标定和正则化系数选择这两块多聊几次比看十遍文献都管用。本文还有配套的精品资源点击获取