大斜度井中利用交叉分量快速反演井斜角

利用三维感应测井交叉分量之差响应曲线在边界附近的响应行为,建立了求解大斜度井中井斜角快速反演方法.响应曲线随真实垂直深度(TVD)变化的对数衰减率与随测量距离(MD)变化的对数衰减率与之比是井斜角的余弦函数,通过简单三角函数计算可快速反演井斜角.响应曲线形态规则,在层状各向异性地层中与垂直电导率无关,主要反映水平电导率反差造成的地层边界信息.因此本文的井斜角快速反演方法既适用于各向同性地层也适用于各向异性地层.我们考察了测井环境及仪器参数对响应曲线衰减率的影响以及衰减率估算偏差对井斜角反演精度的影响.数值模拟结果验证了该反演方法在大斜度井中的有效性和实用性.     

水平井和大斜度井中的感应测井响应计算

感应测井受井斜影响十分严重，在水平井和大斜度井中，感应测井曲线几乎面目全非．水平井和大斜度井中感应测井响应的正演计算是正确认识这种复杂环境中感应测井响应特征的重要手段，是进行井斜环境影响校正的必要准备，是研究用感应测井曲线探测水平井井眼离最近界面距离的基础，同时也是感应测井理论的新发展．本文将论述水平井和大斜度井中感应测井响应正演计算的理论和方法；分析水平井、大斜度井中感应测井的响应特征     

张量感应测井资料处理中若干问题的研究

在横向各向同性（TI）介质中,张量感应测井响应与地层的水平和垂直电导率σ_h和σ_v、井眼斜角α、仪器方位角φ同时有关,其资料处理和解释相当复杂,一般采用先定出φ,然后在井眼坐标系中同时反演σ_h、σ_v和α的资料处理方式.我们发现将测井响应变换到井眼坐标系时实际上无法确定其交叉分量的符号而只能确定其绝对值,这个问题尚未引起重视.本文给出了无需确定φ就可由测井响应直接计算井眼坐标系磁感应张量的公式,此方法简单,所得结果不受φ提取精度的影响.在此基础上本文给出了利用井眼坐标系的磁感应张量的交叉分量曲线提取地层纵向边界和判别相邻地层电导率高低的一种特征识别方法.本文还给出了大斜度井中TI介质的电导率视值的一种简化近似公式.     

用于固体矿床多分量感应测井响应模拟的矢量有限元法

本文开发了基于非结构化四面体网格的三维矢量有限元法,实现了固体矿床井眼中多分量感应测井响应的数值模拟,并分析了多分量感应测井仪器在复杂矿床模型中的响应特征.本文通过采用几何因子背景场,有效地避免了源的奇异性问题;同时,在井眼边界采用非均质网格并用Gauss-Legendre积分计算四面体单元的等效电导率.利用LU分解求解线性方程组,实现了一次网格划分多点的数值计算,提高了计算效率,从而实现快速连续的多分量感应测井模拟.非结构化的四面体网格确保了该方法可以模拟实际问题中所能遇到的复杂的矿体模型.基于水平三层分层和径向分层模型,验证了算法在各向同性和各向异性介质中的可靠性.我们还以三个不同的矿床模型为例,研究了多分量感应测井仪的不同分量的探测特性,结果表明,结合九个分量的信息,可以探测矿体的深度,也可以识别矿体的方位和走向,为精确地描述矿体的三维分布特征打下了基础.     

随钻方位电磁波测井反演方法研究及在水平井中的应用

针对随钻方位电磁波测井仪器测井响应特征进行考察,提出了水平井和大斜度井多层界面三维电阻率反演方法,获取准确的水平井/大斜度井地层电阻率剖面和井眼地层位置关系,为储层评价提供可靠依据.以哈里伯顿ADR方位电阻率测井数据为例,利用基于Marquardt方法的三维电阻率反演方法计算地层电阻率和多层边界距离等参数,数值模拟及实际处理井结果显示,反演方法精度较高,反演结果准确可靠.多层界面三维电阻率反演对储层精确评价和储量复算具有重要意义,在实时地质导向中具有广阔应用前景.     

地震波数值模拟的非规则网格PML吸收边界

以格子法为基础,以声波方程为例研究非规则网格PML(Perfectly Matched Layer)方法.本方法的核心是建立局部坐标系下的分裂方程和基于积分近似的微分方程弱形式.该非规则网格模拟方法允许在计算域内设置任意形状的人工边界.对于二维半空间问题,与采用矩形人工边界相比,采用半圆形人工边界可减少计算量20%以上.采用光滑的曲边界,不仅可减少计算区域,还可避免常规的PML吸收边界在吸收带角点区域的特殊处理.本方法事先计算和存储边界单元的局部几何参数,在计算的每一时间步查表调用这些参数,与常规的直边界PML方法相比,不增加任何计算量.     

横向同性地层中井间电磁响应的三维频域耦合势正演模拟研究

本文针对电导率横向同性地层中三轴发射-接收线圈系的井间电磁响应,开展三维正演算法及模拟研究.首先将发射线圈简化为三轴正交磁偶极子源.为克服低频电磁产生的数值迭代收敛缓慢问题,利用低感应数预处理方法将频域Maxwell方程组转化为基于矢势与标势的Helmholtz方程.采用非等间距的Yee氏交错网格、电导率的体积加权平均技术以及积分形式的有限差分格式实现对控制方程的离散化.借助于不完全LU分解预处理的稳定双共轭梯度法计算井间电磁响应的数值解,并通过与有限元软件COMSOL的数值结果对比验证本文模型与算法的有效性.数值计算结果表明：在横向同性地层的垂直井模型中,测井响应横向分量xx和yy能反映出地层电导率各向异性信息,但横向探测能力较差.轴向分量zz虽只反映地层水平电阻率信息,但具有较好的横向探测能力,且对地层中异常体非常敏感,而交叉分量xz、zx则具备较强的层边界识别能力.     

含环状天线槽正交方位随钻电磁波测井仪器三维有限体积法数值模拟

针对由多环状天线槽构成的方位随钻电磁波测井仪器真实结构,基于柱坐标系下耦合势Helmholtz方程和三维有限体积法研究建立一套相应的三维电磁响应数值模拟算法.首先,采用柱坐标系下三维不规则Yee氏交错网格对计算区域进行剖分,保证了钻铤表面附近的剖分网格与非光滑刻槽钻铤表面同形;并利用电场延拓边界条件逼近钻铤表面的理想电导体边界条件,保证非光滑刻槽钻铤表面耦合势方程离散精度.在此基础上,应用标准均质化技术确定非均质网格的等效电导率和磁导率,并借助三维有限体积法对柱坐标系中耦合势方程以及环状天线槽内的发射天线进行离散,得到一个大型非对称稀疏代数方程,通过PARDISO直接求解器求解离散方程.最后,通过与柱状均匀地层中模式匹配算法所得的数值结果的对比验证该三维算法的有效性,并进一步通过数值结果研究考察金属钻铤和天线槽对方位随钻电磁响应的影响,以及层状TI地层模型中该仪器的响应特征.     

层状各向异性地层中三维感应测井响应快速计算及资料处理

本文采用广义反射系数法推导了水平层状各向异性地层中电磁场的积分解析解,并利用快速汉克尔变换技术实现了三维感应仪器测井响应的快速计算.三维感应测井响应与地层水平电导率、垂直电导率和井斜角及仪器方位角同时有关,单一分量的测井曲线不能满足资料解释的需要.通过对仪器测量分量响应特征的考察,本文提出了一种基于组合量测井曲线的资料直观解释方法.数值模拟显示,交叉分量相关组合量可准确划分地层纵向边界,并可直观识别各向异性层;与单独分量相比,主分量相关组合量提高了纵向分辨率、减弱了与地层电导率参数的非线性关系.     

基于耦合标势与矢势的有限体积法模拟非均匀各向异性地层中多分量感应测井三维响应

为快速有效地研究、考察各向异性地层条件下多分量感应测井的响应特征,本文利用电场标势与矢势的有限体积法研制出三维频率域电磁场响应的数值模拟算法,克服由低频发射或高阻地层产生的低感应数问题,有效提高了三维电磁数值模拟算法的应用范围和计算效率.首先利用电场的标势与矢势将Maxwell方程转化为满足库仑规范条件的耦合势Helmholtz方程,以Yee氏交错非均匀网格中不同位置上的节点为中心建立四种控制体积单元,通过对控制体积单元中电磁场与电导率的积分平均实现耦合势方程和磁偶极子旋度的离散,并得到一个对角占优的大型稀疏复线性代数方程组,然后,通过不完全LU分解预处理和稳定双共轭梯度法快速求解离散方程.数值结果证明了该算法的有效性,并进一步考查了仪器偏心、倾斜井、垂直裂缝等复杂条件下多分量感应的响应特征.     

随钻核磁共振测井的地层界面响应特征

随钻核磁共振测井的地层界面响应特征对地质导向和原状地层评价具有重要意义,但钻井轨迹的复杂性决定其响应特征是多重因素综合作用的结果,且较难直接给出统一显式表达式.本文根据仪器在地层内的运动特征,建立了基于敏感区域体积元的随钻核磁共振测井响应方程和离散化计算方法;通过正演、反演结合的数值模拟,研究了单、双界面水平层状地层中,不同井斜角度、天线长度和目标地层厚度等条件下的随钻核磁共振测井响应特征,对关键结果给出了成立条件和定量关系.随钻核磁共振测井在斜井段中的T₂分布和孔隙度响应与直井中差异明显;井斜角越大,视地层厚度越大,地层界面在T₂分布和孔隙度曲线上的过渡段越长,T₂分布过渡起始位置已不能确定界面深度;仪器的最高纵向分辨率与天线长度、井眼条件和仪器参数有关;受围岩作用影响,较薄目的层的测井响应特征上可能出现异常"夹层".基于分析结果,对随钻核磁共振测井资料解释提出了建议.     

计算全三维电磁测井响应的有限元程序(英文)

同时考虑电导率和位移电流的影响 ,导出了基于A -Φ的一个新泛函表达式 ,并开发了能在宽频段内模拟电磁测井响应的三维有限元程序。计算了水平磁偶线圈和垂直磁偶线圈的场分布 ,也给出了倾斜井中感应测井的响应 ,最后分析了线圈半径和各向异性的影响     

Study on the simulation of acoustic logging measurements in horizontal and deviated wells

The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.     

用磁流源并矢Green函数的递推矩阵方法计算层状各向异性地层中多分量感应测井响应

采用递推矩阵方法计算各向异性介质的磁流源并矢Green函数,并利用上述并矢Green函数对层状各向异性倾斜地层中多分量感应测井的响应进行数值模拟,分析了线圈距、层厚、倾角和围岩对多分量感应测井响应的影响.计算发现,共面视电导率比共轴视电导率变化规律复杂,不能反映地层电导率的真实情况.各向异性地层共轴视电导率随井眼相对倾角的增加而减小,而共面视电导率随井眼相对倾角的增加而增加.仪器的垂向分辨率、围岩各向异性对目的层响应的影响程度取决于线圈距,围岩各向异性对高电导率目的层中测井响应的影响大于对高电阻率目的层中测井响应的影响.     

反插值法实现地球物理数据快速网格化

本文借鉴反演理论，采用反插值法实现地球物理数据的快速网格化.首先由已知点与未知网格点的反距离拓扑关系建立反演方程，用已知点值直接计算出其所在网格的未知网格点值，再利用Laplacian算子滤波，使得模型光滑，且能量最小化.利用预条件共轭梯度法求解网格化方程组，并结合螺旋坐标系思想和Wilson\|Burg谱分解法，将二维数据的滤波处理转换到一维空间进行处理，实现二维数据快速滤波.另外，引入了阻尼约束，保证求解稳定、迭代收敛.最后，应用该方法对合成数据和实际资料进行了试验.     

Analysis and application of the response characteristics of DLL and LWD resistivity in horizontal well

There exist different response characteristics in the resistivity measurements of dual laterolog (DLL) and logging while drilling (LWD) electromagnetic wave propagation logging in highly deviated and horizontal wells due to the difference in their measuring principles. In this study, we first use the integral equation method simulated the response characteristics of LWD resistivity and use the three dimensional finite element method (3D-FEM) simulated the response characteristics of DLL resistivity in horizontal wells, and then analyzed the response differences between the DLL and LWD resistivity. The comparative analysis indicated that the response differences may be caused by different factors such as differences in the angle of instrument inclination, anisotropy, formation interface, and mud intrusion. In the interface, the curves of the LWD resistivity become sharp with increases in the deviation while those of the DLL resistivity gradually become smooth. Both curves are affected by the anisotropy although the effect on DLL resistivity is lower than the LWD resistivity. These differences aid in providing a reasonable explanation in the horizontal well. However, this can also simultaneously lead to false results. At the end of the study, we explain the effects of the differences in the interpretation of the horizontal well based on the results and actual data analysis.     

用并矢Green函数的矢量本征函数展开式模拟随钻感应测井仪器的响应

采用圆柱坐标系下完纯导电劈模型的磁流源并矢Green函数的矢量本征函数展开式计算钻铤V型槽内随钻感应测井仪器的响应.计算发现,V型槽张角的变化对电磁场不同分量的影响规律不同.若只采用传统的感应测井仪器,V型槽的张角可尽量小一些以最大限度地增加信号强度.若采用多分量感应测井仪器,需对V型槽张角的取值进行折衷考虑.由于在具...     

井间地震波动方程正演数值模拟及其应用

井间地震正演模拟是在给定地层参数及边界条件的情况下求波动方程数值解的过程.本文详细推导了井间二维声波方程四阶差分方程;根据假设的包含楔形体和丘状体的理论地层模型,对其波场进行了模拟计算分析;据胜利油田垦71区块两口井的实际声波测井资料建立了精细速度模型,对其正演模拟的合成记录波场进行分析,并把合成记录反射波成像与实际的过井三维地震剖面进行了对比,验证了所设计观测系统的合理性.