ARTICLE DETAIL

资讯详情

深耕网站建设与运营推广的一线实战洞察。

水平井DTS温度剖面反演产出剖面:从原理到应用

水平井DTS温度剖面反演产出剖面:从原理到应用 简介这是一份围绕水平井产出剖面解释的博士学位论文资源核心是利用分布式光纤温度测试DTS技术开展理论建模与反演应用适合具有油气藏工程、生产测井背景的技术人员和研究者。整个资源包只有一个PDF文件20.69MB目前已有63人参与学习。内容系统地从油藏渗流和井筒多相流动机理出发综合考虑微热效应、热膨胀、热传导、热对流、粘性耗散以及地层伤害等因素构建了温度预测复合模型在此正演模型基础上又分别建立L-M与MCMC反演模型通过对不同出水形式底水突破、油水同层、不同出水位置及渗透率分布的真实数据进行拟合定量反演温度与压力剖面结果表明MCMC方法的解释精度优于L-M尤其在确定水产出位置方面更为突出。最后针对实际油田DTS测试数据进行了实例验证与传统生产测井结果一致性高可有效减少找水作业时间与成本为水平井动态监控和控水稳油提供可靠的技术手段。 水平井产量下来了井口计量只有一个总数可你想知道水平段哪一段在产、哪一段在出水、哪些压裂段根本没干活。传统生产测井仪器在直井里还能凑合一到水平井就麻烦爬行器推不动、多相流流型混乱、解释结果五花八门。分布式光纤温度测试DTS给了我一个完全不同的思路——把整个水平段变成一根连续的“体温计”再通过温度场的变化反推出各段产出贡献。这篇博文就把我从温度剖面到产出剖面的解释理论和应用过程完整写下来包含可复现的模型思路和参数也分享几个容易踩的坑。1. 水平井产出剖面测试的痛点仪器下不去、测不准、还多解水平井在现在油气田开发里已经是绝对主力特别是页岩油气、致密砂岩和低渗油藏基本是“水平井分段压裂”的组合拳。但是开发效果怎么评价各段产了多少油、多少水、多少气这问题直接关系到后续调剖、堵水、重复压裂的决策。可惜的是水平井的产出剖面测试一直是个老大难。传统生产测井PLT依靠涡轮流量计、持率计、压力计、温度计等组合仪器在直井或大斜度井里比较成熟到水平井段就开始拉胯。首先是力学问题仪器要靠爬行器才能进水平段井筒里如果有结垢、出砂、弯段曲率过大爬行器很容易卡住。其次是流态问题水平井段里气液两相或油水两相会出现严重的重力分异水倾向于走底部气或轻质油倾向于走顶部流型从泡状流到段塞流再到分层流变化很快。PLT在这么大流动条件下测到的局部持率和流速几乎不能直接代表整段流量。更麻烦的是解释多解性。同样的井口产量不同的分段产出组合可能会产生相似的PLT响应尤其在高含水井里涡轮流量计的启动门槛高低产段根本转不起来。很多情况下解释人员只能靠经验“分配”产量精度自然没法保证。于是行业里一直在找一种不干扰生产、能覆盖整个水平段的连续测量手段分布式光纤温度测试就是这么被推到前台的。DTS最大的优势是把离散的点测变成连续的线测。光纤沿水平段全程铺设测出来的温度剖面是整口井热状态的连续快照。更关键的是产出流体的温度效应是“自带标签”的——油、气、水从地层进入井筒时因为焦耳-汤姆逊效应会在对应的射孔段或压裂缝位置产生可识别的温度异常。只要能把这些温度异常定量解释出来就能反推各段产出贡献。这就是全文的核心逻辑温度剖面是载体产出剖面是目标中间靠传热模型和反演算法搭桥。2. DTS测试原理与现场实施怎样把油井变成一根“体温计”2.1 拉曼散射与温度解调DTS的原理本质上不复杂一句话概括就是“用光测量井筒温度”。井下布设一根单模或多模光纤地面仪器向光纤中发射窄激光脉冲光在传播过程中会发生背向散射其中拉曼散射包含斯托克斯光和反斯托克斯光两个分量。关键点在反斯托克斯光对温度非常敏感而斯托克斯光对温度基本不敏感所以通过两个光强信号的比值就能解调出光纤所在位置的温度。这一比值与温度之间有明确的指数关系经过系统标定后温度测量精度能做到0.1℃甚至更高。我实际用过的DTS系统常见的空间分辨率在0.5米到1米之间采样间隔可以做到0.25米到1米测试一条2000米的水平段大概十几秒到几分钟就能完成一次全井剖面的扫描。这个参数组合对产层识别来说很关键空间分辨率决定了你能分辨的最小产出段长度而采样频率决定了能不能捕捉到瞬态温度变化。参数典型值对解释的影响空间分辨率0.5 ~ 1.0 m决定可识别的最小出液段宽度采样间隔0.25 ~ 1.0 m影响深度定位精度温度分辨率0.01 ~ 0.1 ℃影响小产量段识别能力测量时间10 s ~ 几分钟影响瞬态过程捕捉能力最大测量距离几公里覆盖水平段长度这里要提醒一个误区空间分辨率不等于采样间隔光纤上某一测量点的读数实际上是该点周围一定范围内的平均温度。解释的时候不要把局部尖锐异常当成真实裂缝位置得先做空间响应修正我后面会细说。2.2 三种典型的光纤布设方式光纤要想到达水平段施工方式直接决定资料质量。目前工程上主要分三类第一是永久安装在完井时将光纤绑在油管外壁或套管外壁随管柱下入固井或坐封后光纤就和井筒固定在一起。这种方式信号最好光纤与地层、流体接触充分但后期一旦光纤断裂或损伤没法更换。页岩气井和海上井用得比较多因为测完产剖之后还要长期做压裂监测和井完整性监测。第二是连续油管/过钻杆方式把光纤预置在连续油管或专用保护管内作业完成后可以回收复用。这种方式灵活适合老井接产或者临时监测但连续油管本身会占据一定的生产通道影响流量同时光纤与井筒的耦合不如永久安装紧密温度响应会有滞后。第三是毛细管井筒内自由悬挂光纤封装在毛细钢管里靠自重下入水平段适合低产液、无复杂管柱的井。成本低但只能测到井筒流体温度测量不到套管外或地层温度解释精度受限制。我自己的建议是如果目标是水平井长期产出剖面监测优先选永久安装如果是一次性测试或老井评价可以用连续油管方案。布设方式还会直接影响热传导模型里“光纤所在地层/井筒位置”的边界条件解释的时候必须注意。2.3 现场测试质量控制深度校正与光缆蠕变DTS原始数据拿到手第一件事不是解释而是做质量控制。光纤深度的准确标定是所有解释的基础井下光纤经常有接箍、封隔器等磁性标记通过原始信号幅值上对应的接箍位置与工程管柱表核对做一阶线性深度校正才能把温度剖面对到准确的测深位置。光缆蠕变是我在现场吃过大亏的地方。永久安装的光纤在高温高压下护套材料会发生蠕变导致光纤与管柱之间出现局部松弛表现在温度剖面上就是深度错位。特别是注气热采井这个问题尤其严重必须定期用接箍标记重新校正深度。另外还要注意损耗补偿。一段长光纤可能存在几处微弯或熔接点导致后续光强衰减解调温度出现系统性偏差。通常会在光纤末端设置一段已知温度参考段比如井口附近的恒温段用参考段实测值反算损耗系数并做动态补偿。3. 温度剖面和产出剖面之间的物理桥梁焦耳-汤姆逊效应与井筒换热温度剖面本身并不是产出剖面。它只是结果产出的流体和地层、井筒之间的热交换共同决定了最终温度场。想从中提取产出剖面必须先搞清楚温度到底受哪些因素控制。第一个主导因素是地层原始温度梯度。浅层温度低深层温度高水平段一般沿固定深度延伸所以如果没有流体产出光纤测到的温度值会接近地层温度且沿井筒方向变化很小。一旦有流体产出温度剖面就会“偏离”地温线。第二个关键因素就是焦耳-汤姆逊效应Joule-Thomson Effect简称JT效应这是解释的核心物理机制。流体从地层孔隙渗流到井筒压力从地层压力降为井底流压这个节流降压过程中流体的温度会发生显著变化。不同相态的流体JT系数正负和大小差别很大水JT系数在常见压力温度范围内基本为负也就是说产水会导致局部温度下降。水产量越大压降越大温度下降越明显。油JT效应也比较明显同样表现为节流冷却中等产量下温度下降幅度在零点几摄氏度到几摄氏度之间。气JT效应更复杂。在低压范围气体节流通常变冷但在高压条件下很多气体特别是甲烷可能出现节流升温。处理气井时必须结合具体压力和温度状态查JT系数表不能想当然。正因如此不同流体在同一段层位产出时温度剖面的“异常形态”是不一样的。产水段往往出现清晰的低温凹陷产气高压段可能出现温度上升而产油段的响应介于两者之间。再加上多相混合后的流体在井筒中继续与周围地层换热整体温度剖面就会被重塑。第三个因素是井筒流体与地层的持续换热。流体从产出位置沿井筒流到井口沿途每时每刻都在和套管、水泥环及地层之间进行热交换。在关井或低产情况下井筒温度趋近于地温梯度生产时冷流体低于地温沿井筒流动会被逐渐加热热流体高于地温则被冷却。这种环境换热会“抹平”局部的温度异常导致裂缝位置在远端的温度剖面上可能只表现为温差变化率的转折点而不是一个明显“坑”。所以解释界有一个共识指望单看温度剖面上的“冷点”就直接判断产液段太天真了。温度剖面上能看到的是综合热过程的叠加结果要还原出各段产量必须建立考虑流动、压降、传热耦合的正演模型再用反演方法求产量分布。有朋友问过我既然关井温度基本等于地层温度那温度异常越大的地方是不是产液越多这个说法在单井单相产水时勉强成立但一旦存在多相、多段产出、段间流体混合或者井筒换热强的情况温度异常幅度和产量之间就不是单调对应了。我在实际算例里遇到过中等产量气段的温升比大产量油段的温降还明显如果只用温度幅度定性解释会严重误判主产段。4. 从温度场反演产出剖面的数学模型与求解流程4.1 正演模型井筒-地层耦合传热可复现的解释流程必须有扎实的正演模型。我的方案是采用经典的井筒-地层轴对称瞬态传热模型配合沿井筒的流体质量守恒和动量守恒方程形成一套完整的“流动-传热”联合求解过程。先设定简化的单井物理模型水平井段沿深度不变井筒半径为r_w光纤布置在套管外壁。每段射孔或压裂缝段被视作一个“流量汇/源”从该位置向井筒内注入一定质量的产出流体。假设地层径向无限大沿井筒方向的热传导远小于径向地层内热传导按一维径向瞬态导热计算常用的解析近似是Ramey提出的无因次时间函数也可以用有限差分差分求解径向温度场。井筒内的流体温度控制方程本质是能量守恒[ \dot m c_p \frac{dT_f}{dz} q_{form} q_{JT} q_{fric} ]每一项的含义左边是沿井筒方向流体携带的热量变化右边第一项是地层对井筒的径向导热第二项是节流温降第三项是摩擦生热。对水平井来说重力项在水平段内可以忽略但压降和流型对换热系数影响很大不能省略。[ \frac{dp}{dz} - \rho g \sin\theta - f \frac{\rho v^2}{2d} - \frac{1}{2} \rho v \frac{dv}{dz} ]其中摩擦系数f要根据多相流流型选择合适的经验关联式比如Beggs-Brill或Lockhart-Martinelli。油气水三相流动的持率、折算速度和混合温度都要按局部压力状态迭代计算。这个正演模型的输入是一套“假设的产物剖面”——每段流量和含水/气率输出是沿井筒的温度剖面T(z)。正演算得越准反演反推的产量就越可靠。实操中我会先用稳态模型做初始正演再叠加生产时间内的瞬态恢复过程。4.2 反演方法目标函数、正则化与不确定性反演的目标是找出一组各段流量、含水率参数让正演温度剖面尽可能逼近实测DTS温度剖面。最常见的数学形式是最小二乘目标函数[ J(\mathbf{x}) | T_{obs} - T_{sim}(\mathbf{x}) |^2 \lambda | \mathbf{R} \mathbf{x} |^2 ]第一项是数据拟合误差第二项是正则化约束项。这里R矩阵可以是差分矩阵约束相邻段的产量不要剧烈跳动λ是正则化系数要靠L曲线或交叉验证选取。为什么要加正则化因为这就是一个典型的病态反问题。温度剖面是有限分辨率测量空间内连续分布的产量参数却远多于测量点数量不加以约束反演结果会出现剧烈的振荡某一小段产量为负、另一段产量翻倍。加上平滑约束后才能获得合理且稳定的分段产量分布。我目前比较推荐的做法是两阶段反演粗定位先对温度剖面求空间导数识别温度梯度突变位置用这些位置作为候选产液段边界减少未知参数维度精定量在固定段位置的情况下反演每段的液量和含水率。这个子问题未知数少可以用Levenberg-Marquardt或贝叶斯MCMC方法求解顺便得到参数的不确定性区间。贝叶斯方法还有个额外好处它能给出后验分布从而评估解释结果的可信度比如某段产量95%置信区间是20~50方/天还是35~40方/天这对现场决策非常重要。4.3 可复现模拟算例六段水平井的温度反演我设计了一个便于复现的模拟算例以满足“按参数自己就能跑通”的要求。基本参数如下水平段长度1000 m油藏中部垂深1800 m地温梯度每100 m约2.5℃井筒内径0.1 m生产压差1.8 MPa6条压裂缝段分别位于500、600、700、800、900、1000 m测深位置单段长10 m缝间间距约100 m总液量200 m³/d含水率按段变化第一段10%第二段40%第三段5%第四段70%第五段20%第六段90%正演模拟的温度剖面显示高含水段第4、第6段出现显著低温凹陷凹陷幅度与产水量正相关中低产油段温度变化幅度较小但温度剖面的梯度仍有明显转折。从定性的角度看高含水段基本可以一眼识别但产量定量需要依赖反演。我采用上述两阶段反演方法处理后结果如下段位置 (m)设定液量 (m³/d)反演液量 (m³/d)含水率设定 (%)含水率反演 (%)5002021.510126003534.140387001514.2568005052.370669002523.8202210005554.09088液量绝对误差基本在1~2 m³/天以内含水率误差在2~4个百分点。对于工程应用来说这个精度足够定位主产液段和高水段了。建议的复现文件结构是dts_profiling/ data/ raw_dts.csv # 原始温度-深度数据 completion.csv # 管柱和射孔/压裂段表 models/ wellbore_heat.py # 井筒传热正演 inversion.py # 反演主程序 config/ params.json # 所有可调参数集中管理 results/ inverted_profiles.csv figures/我习惯把所有参数集中在params.json里一方面是为了可复现另一方面是现场井参数经常更新改一个参数文件比改代码安全得多。整个算例用Python实现核心求解器用numpy和scipy不需要商业软件普通笔记本跑反演只需要几分钟。5. 应用场景找水、产气剖面和压裂效果评价5.1 分段产液剖面与找水DTS产出剖面解释在水平井上最直接的应用就是找水。前面说过水由于JT效应强在温度剖面上会产生明显的低温异常。实际油田中这种情况特别典型一口水平井产油下降井口含水快速上升需要知道是哪一个压裂段出水。常规测井很难分段测量DTS则可以通过对比不同时间测得的温度剖面捕捉含水段位置的温度降低进而指导堵水措施。我在一个实际油井案例里见过这种情况DTS温度剖面上有两处明显低温凹陷其中一个是原本认为的主产油段另一个是之前射孔后从未参与改造的低渗段。通过温度反演发现低渗段其实在大量产水名副其实的“水层”随即堵掉该段之后井口含水率下降了20个百分点。正因为DTS能提供连续温度资料这种找水定位的准确度要比传统经验判断高很多。5.2 压裂后返排评价水平井分段压裂后各段的返排贡献率直接决定压裂改造效果。常规做法是每段下封隔器或桥塞分别放喷计量费时费力又影响投产进度。DTS可以做到压裂后直接进行连续温度监测压裂液进入地层后温度通常低于地层温度返排期间低温流体从各裂缝段汇入井筒会在每一段位置留下温度异常。通过监测关井后温度恢复速率和返排时的温度剖面再结合反演模型可以估计每段的返排液量贡献。把这套方法和微地震事件密度对比往往能发现“温度显示产出的段”和“微地震显示压开的段”并不完全重合这种情况对压裂设计的价值非常大微地震反映的是裂缝几何形态而DTS反映的是有效产出能力两者叠加才是完整评价。5.3 注水/注汽剖面监测DTS不光用于产出剖面也经常被部署在注入井里做剖面对应监测。注水井中各段吸水量差异大吸水多的层段温度响应明显。特别是注冷水时低温水的推进会在温度剖面上表现出明显的“冷前锋”根据冷前锋推进速度能估算吸水剖面。热采井注汽时DTS能测到蒸汽沿水平段的分配情况。蒸汽沿井筒流动时会产生极其明显的温度响应高干度蒸汽到达某段该段温度会迅速上升。通过监测各段升温时序可以判断蒸汽突破位置优化注汽参数防止汽窜。气体或蒸汽的JT效应在高压井中可能表现为升温解释时一定不要套用水井的解释模板。5.4 与DAS/PLT联合解释DTS虽然好但只有温度一个物理量。近年来分布式声波传感DAS逐步成熟声波信号对流体流动和段间活跃度非常敏感。DTS和DAS联合解释的思路是DAS负责“看谁在动”——识别各段是否有流体流动、流动是否持续DTS负责“看动多少”——定量反演流量。两者一静一动配合起来能大幅降低单一温度反演的多解性。如果井内有条件再下测一次PLT哪怕只测到倾斜段也可以作为“锚点”校正DTS反演结果的绝对产液量。温度剖面反演能给出各段的相对贡献分布但总流量通常依赖井口计量如果井口计量不准绝对产量就会系统偏移。此时一条PLT的实测总液量数据就能非常有效地消除这个偏差。6. 实操经验与避坑清单从资料质量到解释结果验证DTS解释的成败一半在资料质量一半在算法。我在好几口井上栽过跟头总结下来的实操要点都在这张清单里先检查原始数据。如果光纤某一段信号强度骤降那很可能是微弯或熔接损耗对应的温度数据基本不可用不要硬解释。深度校正必须先做用接箍信号确认光纤实际深度与工程测深一致否则后续所有定段都会偏。关注时间窗口。刚开井或刚换工作制度时井筒温度处于瞬态恢复阶段此时测到的温度剖面偏离稳定产出的热状态如果直接拿来做定量反演结果误差非常大。我一般会等生产稳定至少几天后再取“准稳态”温度曲线如果要做瞬态解释必须让正演模型把时间维也考虑进去。反演多解性问题。就算加了正则化产量在线性相关的情况下仍然可能多解。实战中我的做法是同时用多个不同的初始猜测做反演看结果是否收敛到同一组解如果不一致说明约束不够需要添加先验信息比如压裂段位置、射孔段长度、地质认识等。把待反演的段数尽量压缩也是降低多解性的有效方法。验证结果是最后一步。我会把反演出的各段液量相加和井口计量对比相对误差超过10%就要回头找问题。有条件的井用PLT或井下流量计做单点验证没有条件时至少做敏感性分析改变正则化系数、段边界位置看产量结果变化多少以此估算不确定度。还有一个容易被忽视的经验相对产量比绝对产量可靠。温度反演对段间产出差异的辨识能力强对总流量的绝对标定能力弱。所以汇报时不要只给一张“各段产液量”的表要同时给出各段相对贡献占比和置信区间让地质工程师把相对贡献和地质认识去对照这样解释结果才更容易被现场接受。最后分享一个小技巧做DTS解释前先把实测温度剖面减去“无水无流”的地温投影线得到“温度偏差剖面”。这一步看起来简单但能显著放大局部异常帮你在定性阶段快速定位疑似产液段给后续反演提供更好的初始分段。这条经验我用了很多年对新手特别友好。本文还有配套的精品资源点击获取
返回列表