套管非匀质性对过套管电阻率测井影响的数值模拟与分析

套管非匀质性对于过套管电阻率测井响应会产生显著影响.首次采用数值模式匹配法模拟过套管电阻率测井响应,在此基础上分析了各种套管局部变化对于测井响应的影响.模拟结果表明,套管异常区域及其附近的测井响应出现异常变化,异常变化程度与套管参数局部变化程度呈现正相关性;通过模拟不同类型电极系还发现,双注入电极系对套管非匀质性具有很好的补偿作用.文中首次模拟了井温梯度对过套管电阻率测井响应的影响.计算结果表明,该数值方法能够很好地解决包含套管非匀质性因素的过套管电阻率测井数值模拟问题.     

水平层状介质中双侧向资料的全参数正则化迭代反演与应用

水平层状介质中的模型向量通常由地层原状电阻率、侵人带电阻率、侵入半径以及地层界面深度等参数组成,整个测井记录实质上是模型向量的非线性函数.本文将研究双侧向仪器Born近似表达式并建立用双侧向测井资料同时反演所有模型的全参数正则化迭代算法(RWPII).双侧向电极系上的电流和电位必须满足特定的聚焦条件,其正演模拟与Born近似表达式与非聚集型仪器完全不同.应用摄动原理,首先研究双侧向仪器响应的Born近似,给出电极系上的电流密度、电流和电位微小变化与电导率摄动之间的积分关系.然后,根据层状地层中电阻率分布特征以及半解析正演模拟技术,建立Fréchet导数矩阵的快速算法,该Fréchet导数矩阵包含了视电阻率对所有模型参数的偏导数.在此基础上,研究同时反演所有模型参数的正则化迭代算法,并给出正则化因子的具体计算方法.最后,应用理论模型和实际资料反演结果验证该反演算法的有效性.     

基于阵列感应与自然电位联合反演地层水电阻率

原状地层水电阻率是重要的储层参数,也是进行精细储层评价的基础.基于泥浆侵入数值模拟与侵入过程中井周岩石物理特征分析,确定了薄膜电位的产生位置,针对储层高、低侵等不同侵入特征,提出了可适用于包括存在"低阻环"等不同侵入特征时储层电阻率分布的"五参数"地层模型,基于几何因子理论与有限元方法,建立了阵列感应与自然电位测井联合反演算法,实现了地层电阻率参数反演,重构了地层径向电阻率剖面,进而精确求取了地层水电阻率.通过对实际资料处理表明:反演算法稳定可靠,"五参数"模型能很好地表征储层电阻率分布形态,重构储层电阻率剖面,确定薄膜电位产生位置;基于阵列感应与自然电位的联合反演,能精确计算原状地层水电阻率,为储层评价与流体性质识别提供依据.     

基于非结构网格有限元的三维井地电阻率法任意各向异性正演模拟

地层介质的电各向异性增加了井地电阻率法响应的复杂性,开展基于电各向异性介质模型的井地电阻率法响应规律研究对于正确解释各向异性显著地区的观测数据特征至关重要.针对垂直线源井地电阻率法的任意各向异性响应模拟问题,本文提出了一种基于非结构网格有限元三维正演算法,通过引入3×3的对称正定张量来表征任意各向异性的电导率,采用非结构四面体网格有限元方法来离散电位的边值问题,通过将垂直线源等效为一系列点源问题,进而实现了任意各向异性介质中井地电阻率法的高效数值计算.通过与三个地电模型的解的对比,验证了本文数值解算法的精度和有效性.针对线源远离和垂直穿过异常体的两类模型,分别考察了当围岩或异常体为各向异性介质时的井地视电阻率响应特征.结果表明,对于各向异性地层,围岩和异常体的主轴电阻率值和旋转角均会对井地视电阻率的幅值及分布产生显著的影响.研究结果对于提高井地电阻率法的认识和资料解释水平具有重要的理论和实际意义.     

双侧向型复电阻率测井响应的数值模拟

在二维地层模型条件下,应用有限元方法模拟了双侧向模式的复电阻率测井仪器在不同频率下的响应,结果表明当频率很低时,复电阻率测井的响应接近双侧向测井的响应;随着频率的提高,复电阻率测井的响应幅度降低。当地层具有一定的频散效应时,可以利用不同频率下的复电阻率测井响应直观地区分油水层     

方位梯度电极系的正反演

针对环井周的非均质地层模型,利用三维有限元方法模拟了方位梯度电极系的测井响应．结果表明,在非均质程度比较小的情况下,各个方位上的响应之间差别不大,与真实的电阻率分布差别也不大;反之,则差别很大．其特点是所测得的响应远远低于较高的真电阻率值,而接近较低的真电阻率值．因此,有必要对测井响应进行校正．考虑到传统的校正方法的弊端,本文利用非线性反演方法求取真实的电阻率分布,数值模拟的结果表明了此方法的可行性．     

水平井感应测井响应影响因素数值模拟分析

地层电阻率是影响油气或煤层气评价的重要参数,水平井中常利用感应测井确定地层电阻率,由于水平井与直井测量环境有较大差异,响应方式势必与直井不同.水平井测井响应的环境影响复杂,不能使用基于几何因子的信号处理,必须通过三维正演数值模拟方法研究响应特征.本文利用基于有限元方法的COMSOL Multiphysics软件对三层对称地层模型中的水平井感应测井响应进行了数值模拟分析.通过分析比较源距为0.3m、0.6m、0.9m、1.2m的四个水平线圈系的响应特征发现,不同源距的线圈系测得的电导率值受仪器位置、上下围岩电导率、目的层厚度和仪器偏心的影响规律不同.数值模拟结果对水平井的解释工作及仪器研发具有指导性意义.     

双侧向测井响应计算洞穴充填物电阻率方法

碳酸盐岩缝洞储集体具有丰富的油气资源,洞穴是其最主要的储集空间,电性特征是评价洞穴充填情况的重要指标,利用双侧向测井资料表征洞穴充填物真实的电性特征对洞穴充填情况的准确评价具有重要的作用.应用三维有限元法对等效洞穴地层模型进行双侧向测井响应数值模拟,分析洞穴尺寸、洞穴充填物和基岩对双侧向测井响应的影响.分析结果表明:当洞穴较大时,浅侧向测井响应受基岩和洞穴尺寸影响很小,可直接将浅侧向测井响应作为洞穴充填物电阻率;当洞穴较小时,双侧向测井响应无法直接反映洞穴充填物电阻率,因此应用洞穴尺寸、基岩电阻率和浅侧向测井响应建立洞穴充填物电阻率的计算模型.该模型计算得到的洞穴充填物电阻率与给定的洞穴充填物电阻率具有很好的相似性,在西部某油田实际井资料处理解释中亦取得了良好的应用效果.洞穴充填物电阻率的准确计算为缝洞储集体的测井识别和解释评价提供了有力支撑.     

自适应hp-FEM在随钻电阻率测井仪器响应数值模拟中的应用

随钻电阻率测井仪器响应数值模拟在电法测井仪器研制和资料处理与解释中具有重要的意义.本文利用自适应hp有限元算法(hp-FEM)模拟和分析了地层各向异性以及相对倾斜对仪器接收信号幅值比和相位差的影响,并对仪器响应特性进行了评价;同时,利用HERMES有限元软件包与VC++开发平台相结合,可以对随钻电阻率测井仪器响应数值模...     

用正演数值计算方法开展双侧向测井对裂缝的响应研究

采用有限元分析软件ANSYS对裂缝的双侧向测井视电阻率与裂缝孔隙度、泥浆电阻率、裂缝倾角和基岩电阻率的关系进行了计算. 在大量正演数据的基础上得出双侧向测井响应反演公式和裂缝孔隙度计算公式,提出更为精细的裂缝倾角的弹性划分模型,用双侧向视电阻率值近似估算裂缝倾角的方法,提高了利用双侧向测井求裂缝产状与裂缝孔隙度（裂缝宽度）的精度.     

盐水泥浆条件下定向井双侧向测井环境校正方法研究

在定向井中,双侧向测井数据由于受井眼、围岩、井斜、径向侵入等环境因素的影响,难以准确地估算储层的电阻率和更准确地识别储层流体.本文拟提出图版法和反演法相结合的思路来研究地层视电阻率的校正问题.首先,针对研究区的测井环境和钻井液电阻率,计算了不同尺寸井眼的双侧向测井响应,绘制了井眼校正图版;然后,针对斜度井地层模型,采用三维有限元方法计算了不同井斜、不同厚度地层的测井响应,构建了井斜-围岩/层厚图版,利用该图版实现了井眼校正和井斜-围岩/层厚的电阻率快速校正.最后,针对钻井液侵入的影响,采用反演方法计算了侵入半径和地层真电阻率,最终实现了双侧向测井的环境校正.利用上述理论与方法对海洋中定向井的双侧向测井资料进行了环境校正,提高了储层流体定性识别和定量评价的精度.     

斜井泥浆侵入仿真及其阵列侧向测井响应数值模拟

为了研究斜井泥浆侵入储层的井周岩石物理特征,采用广义曲线坐标系统模拟斜井泥浆侵入特性,进而研究泥浆侵入过程的阵列侧向测井响应.数值仿真表明,采用广义曲线坐标系消除了笛卡儿坐标系统在斜井边界处网格划分的锯齿现象;储层非均质造成泥浆侵入不均匀,储层渗透性越好,泥浆滤液侵入越深;泥浆滤液侵入使得井周地层饱和度、孔隙水矿化度的分布发生变化,进而造成斜井井周岩石电阻率分布剖面复杂变化;泥浆侵入过程中阵列侧向视电阻率变化以及深浅视电阻率的幅度差异,有效反映了斜井泥浆侵入特性;文中采用了基于Marquardt方法的阵列侧向测井四参数反演,有效恢复了斜井储层侵入剖面,得到了原状地层电阻率.斜井泥浆侵入特征及其阵列侧向测井响应模拟分析,有助于正确认识斜井井周岩石物理特征和利用阵列侧向测井进行斜井储层评价.     

3D laterolog array sonde design and response simulation

A new three-dimensional laterolog array sonde (3D-LS) is presented. The 3DLS is based on existing high-resolution laterolog array and azimuthal resistivity imaging sondes with radial, longitudinal, and circumferential detection abilities. Six investigation modes are designed using the 3D finite-element method and different investigation depths are simulated based on the pseudo-geometrical factor of the six modes. The invasion profile is described using multi-array radial logs. From the analysis of the pseudo-geometrical factor, the investigation depth of the 3D-LS is about 1.5 m for conductive invasion, which is close to that of the dual laterolog tool but greater than that of the highly integrated azimuthal laterolog sonde. The vertical and azimuthal resolution is also analyzed with the same method. The 3DLS can detect low-resistivity anomalies of 0.5 m thickness and 15° around the borehole for infinitely thick formations. This study lays the foundation for more work on 3D laterolog array sonde for evaluating low-resistivity anomalies.     

低阻油层双侧向测井的反演研究

对于双侧向测井反演，本文将广义逆反演方法用于消除侵入影响的三参数(侵入半径、侵入带电阻率、原状地层电阻率)反演中。通过实际多口井连续曲线的反演处理表明，本文提出的反演方法令理而且有效。同时，广义逆方法提高了反演计算的稳定性和精度。对低阻油层实际资料的处理结果与试油结论一致。     

多辐射场源地空瞬变电磁法多分量全域视电阻率定义

本文针对多辐射场源地空瞬变电磁法理论,建立了相应的多分量全域视电阻率定义方法.分析了利用磁场强度进行多辐射场源地空系统全域视电阻率定义的优点,针对磁场强度的各个分量,提出了各自的全域视电阻率算法,实现了多分量、全时域、全空域视电阻率计算,并分析了偏移距对全域视电阻率的影响.通过调整源的相对位置及电流方向等参数,多辐射源瞬变电磁地空系统不仅可以加强不同分量信号强度,削弱随机干扰,还可更好地分辨地下异常体的位置.通过对所设计模型的处理,证实了多辐射场源地空系统多分量全域视电阻率算法的有效性,也验证了多辐射场源地空系统的优势.     

各向异性层状介质中视电阻率与磁场响应研究

针对任意各向异性地层,利用极向型和环向型标量位函数,导出相应的直流电视电阻率和磁电阻率的磁场响应关系.计算了各向异性地层的直流电视电阻率和磁电阻率响应特征,重点分析了电阻率测深方法与磁电阻率探测方法对地下各向异性介质的探测能力.文中采用状态矩阵的分析方法,首先采用极向型和环向型标量位构造了各向异性层状介质电场与磁场的通解,利用各层界面电场、磁场的连续性及地面激励源的耦合条件,推导了不同层之间电磁场的状态矩阵,建立了空间电场与磁场的递归计算关系.其次,针对递归计算中指数项数值计算的不稳定性,借用状态矩阵的性质,导出了将不稳定指数计算项转化为稳定的指数项的转换关系.针对横向各向同性(TI)介质中极向位与环向位解耦的特点,导出了电磁场的直接积分解.最后,采用解析解验证了算法的正确性,给出了多层各向异性地层模型的视电阻率和磁场响应曲线,分析了直流电法探测裂缝性地层、估计裂缝分布性状的可能性.     

地-空瞬变电磁法电阻率成像研究与应用

地-空瞬变电磁法(Semi-airborne transient electromagnetic,SATEM)凭借适应能力强、探测深度大、实时性强等特点,适用于湖泊、沼泽、山区等地形复杂地区观测.SATEM接收数据量大、精度要求高,传统成像方法基于地面视电阻率计算未考虑飞行高度,开展SATEM的高精度快速成像研究对实际...     

Dual laterolog borehole correction based on dynamic tool constants

The positive and negative difference of deep and shallow resisitivity in formation without invasion is caused mainly by the fixed two tool constants from numerical simulations. A dynamic calibration method for tool constants is proposed based on the effects of the mud and formation resisitivity ratio on the two constants calculated using the finite element method (FEM). Finally, four specific examples are given to validate the dynamic calibration method. It is an automatic borehole correction method and can give more accurate formation resistivity. The method is useful for dual laterolog logging.     

点源两层大地三维地电模型视电阻率边界元解

本文用两层均匀大地人工点源作用下的空间电位分布作为基本解，导出了下伏基岩中陷含三维地质体情况下，关于电位的边界积分方程组，用边界单元法解积分方程求得地表电位，从而求得了一定装置下的视电阻率。本文方法的独特优点是积分只须在异常体边界面上进行了测线可以沿任意方向。用联剖法计算了均匀半空间下含球形异常体的ρｓ曲线，所得结果与其他学者用解析方法所求得的结果一致。用中间梯度法计算下层基岩中含有一长方体，测线