各向异性地层随钻方位电磁波测井探测特性分析（英文）

随钻方位电磁波测井仪器已广泛应用于大斜度井/水平井钻井,但地层各向异性使得地质导向及储层评价更为困难。为定量分析地层各向异性对随钻方位电磁波测井响应及资料处理影响,本文借助正反演技术,详细考察了测井响应对各向异性的敏感性及测井数据反演精度。首先,定义了敏感因子,定量分析了不同磁场分量、合成信号对各向异性的敏感性;然后,模拟了层状各向异性地层中的随钻方位电磁波测井仪器响应,对比分析了对称设计和非对称设计仪器的响应特征及对地层参数的敏感性;最后,探讨了反演模型的选取对各向异性地层界面反演精度的影响。数值模拟及反演结果表明:随钻方位电磁波测井响应受各向异性影响严重,通过对称补偿设计可降低地质信号对各向异性的敏感性,进而增强对地层界面的探测能力;各向异性地层中,视电阻率曲线出现分离及"犄角"等现象,界面附近方位信号则可能出现复杂非线性变化;反演中若不考虑各向异性的影响,难以获取准确的地层界面和电阻率信息,进而可能导致地质导向决策的失误。     

超深随钻方位电磁波测井探测特性及参数敏感性分析

本文给出超深随钻方位电磁波测井的仪器天线设计及探测模式信号定义;分析各探测模式的探测性能及对地层电阻率、倾角及各向异性的表征能力;定义方位角、井斜角及地层电阻率的敏感性函数,定性分析不同模式对地层参数的敏感性.研究结果表明,超深随钻方位电磁波测井各探测模式具有不同探测特性,探测模式Ⅰ利用对称设计识别地层边界,消除井斜角...     

水平井随钻电磁波测井与双侧向测井响应差异及其解释应用（英文）

大斜度井/水平井中,随钻电磁波电阻率测井与双侧向测井由于测量原理不同,引起的测井响应特征有着较大差别。首先依据积分方程法模拟了随钻电磁波电阻率测井在大斜度井/水平井中的测井响应特征,采用三维有限元素法模拟了双侧向测井在大斜度井/水平井中的测井响应特征,然后对比分析了二者在水平地层中的测井响应差异和产生的原因。结果表明,水平井中随钻电磁波电阻率测井和双侧向测井响应差异的主要影响因素有:不同井斜角、地层各向异性、地层界面以及泥浆侵入等。随着相对井斜角的增大,随钻电磁波电阻率测井在地层界面处的测井响应异常大,而双侧向测井在地层界面附近受围岩影响逐渐变得平滑;随着相对井斜角的增加,各向异性的影响逐渐增大,但各向异性对双侧向测井的影响程度要小于对随钻电磁波电阻率的影响。这些差异如不能正确认识将导致水平井解释过程中出现偏差。最后根据实际资料,分析如何这些差异进行水平井测井解释。     

随钻方位电磁波测井多参数快速反演

大斜度井/水平井随钻测井技术已被广泛的应用于各种复杂油气藏中,但传统随钻电磁波测井仅提供仪器轴向磁场分量,井斜和各向异性对测井响应影响耦合在一起,难以剥离,无法基于反演同时获取多个参数,进行相应的地层评价,且不能提供方位信息,难以实现地质导向功能.随钻方位电磁波测井在提供仪器轴向磁场分量的基础上,还提供了仪器横向磁场分量,其提取的地质导向信号,对地层界面非常敏感,受井斜、地层电阻率对比度和各向异性影响弱,可以反演得到准确的界面位置.此外,额外的磁场分量信息,有效降低了井斜与各向异性的耦合程度,可基于反演算法实现多参数联合反演,满足电阻率测井精细解释评价需求.

     

定向随钻电磁波电阻率测井仪器地质导向方法初探

随着石油勘探开发需要的转变,随钻测井技术的地质导向功能正受到越来越多的使用.本文利用各向异性水平层状地层中的磁流源并矢Green函数对定向随钻电磁波电阻率测井仪器的地质导向方法进行了研究.利用由接收线圈的电压定义的导向向量对低阻围岩的指向性并结合定向地质导向信号的变化,可以确定仪器与地层界面的相对位置.对定向数据进行方位成像可以得到地层方位和井斜角的相关信息,并且不受地层各向异性的影响.数值模拟结果表明:线圈距越长,仪器的探测范围越远且成像的对比度越大;频率越高,目的层与围岩层的电阻率对比度及井斜角越大时,成像的对比度也越大,但井斜角越大,定向信号发生明显变化的深度区间会越窄.利用单界面模型的定向地质导向信号制作出的交会图版可以使用插值法估计仪器到界面距离.所得结果为定向随钻电磁波电阻率测井仪器的地质导向方法提供了参考.     

水平井随钻电阻率与双侧向响应差异影响因素分析

在大斜度井和水平井中,不同电阻率测井系列测量的结果并非完全一致,测量结果的差异性分析对于正确认识和评价地层非常重要.依据双侧向测井和随钻电磁波测井原理,采用三维有限元素法模拟双侧向在水平井中的测井响应特征,采用积分方程法模拟随钻电磁波电阻率在水平井中的测井响应特征,通过对比分析值模拟结果,造成大斜度井和水平井随钻电阻率与双侧向电阻率测井响应差异的主要因素有:井斜、各向异性、地层界面以及泥浆侵入.随钻电阻率一般受泥浆侵入影响小,但是它在水平井中容易产生极化现象,有助于地层界面识别,同时受各向异性影响明显;而双侧向电阻率在水平地层中对地层界面不敏感,且受各向异性影响相对较弱,但会受泥浆侵入影响严重;正确认识这些差异,可以有效指导水平井的合理解释.     

定向电磁波测井的倾斜天线系结构研究

采用水平层状各向异性介质中的磁流源并矢Green函数数值模拟了倾斜天线系随钻电磁波电阻率测井仪的电磁波传播特性.比较了轴向天线和倾斜天线在地层界面处的电流特性,从理论上揭示了倾斜天线系电磁波传播随钻测量探测地层界面位置及方位的原理.分析了采用倾斜天线系时电压的幅值衰减和相位移,理论上给出天线系倾角的选择方法.详细阐述了发射频率和源距参数的选择方法,给出了不同发射频率适用的源距范围.研究结果与现有两种国外定向随钻电磁波电阻率测井仪的参数一致,为国内定向测井仪器的研发及测井数据解释奠定理论基础.     

电性各向异性地层频率域电磁响应模拟

在现有计算机内存和运算速度条件下,电性各向异性地层中频率域电磁响应的数值模拟是一个难点. 为探索各向异性介质模型电磁响应的有效算法,同时为进一步研究各向异性地层测量响应的校正方法奠定基础,从Maxwell电磁响应方程出发,使用由三个彼此垂直的发射线圈和三个彼此垂直的接收线圈组成的线圈结构模型,应用交错网格有限差分法推导了三维频率域电磁响应的差分计算格式,采用预条件双共轭梯度法对所形成的复系数矩阵进行了迭代求解,从而得到电性各向异性地层模型下的线圈测量响应. 应用这些方法进行了多种条件下电磁响应的计算,计算结果表明,所形成的计算方法是准确可靠的,预条件双共轭梯度法的运算效率也较高;电性各向异性地层的电磁响应与地层倾角（或井斜角）、仪器测量方位、地层的水平电阻率和垂直电阻率有关;根据水平方向磁场强度分量可以确定地层界面;在电磁响应计算过程中,背景电导率的取值要尽量接近整个计算区域的平均电导率.     

随钻方位电磁波测井响应快速正演方法与地质导向应用

随钻方位电磁波测井解析快速正演是实现仪器响应特征分析、结构优化设计以及地层参数精确反演的基础.本文基于赫兹势函数法求解水平层状各向异性介质中的电磁场, 通过引入地层界面位置, 改造解势函数中的指数增大项, 提取溢出因子, 消除了势函数系数递推法中的数值溢出, 实现了频率-波数域电磁场收敛计算.将算法应用于新型随钻方位电磁波测井技术, 分析自主研发的Directional Wave Propagation Resistivity (DWPR)仪器原始响应信号特征, 提取地质信号、电阻率信号.结果表明: DWPR新型仪器结构对地层界面及电阻率敏感, 可提供准确的实时地质导向信息; 水槽实验和实际地质导向应用结果进一步验证了本文算法的实用性和高效性.

     

用传播矩阵法研究层状正交各向异性地层中多分量感应测井响应

本文采用传播矩阵技术研究并建立了层状正交各向异性地层中多分量感应测井响应的有效算法.首先通过Fourier变换将频率空间域中的Maxwell方程组求解问题转化为频率波数域中关于电磁场水平分量常微分方程组的定解问题.利用该方程组系数矩阵的本征值和归一化本征向量将电磁场分解成上行波和下行波模式的组合,推导出均匀正交各向异性介质中由任意方向磁偶极子产生的电磁波模式解析表达式;在此基础上,利用叠加原理和边界条件研究了电磁波在层状正交各向异性地层中的反射和透射,给出各个界面上的广义反射系数和不同地层中电磁波振幅的递推公式,进而得到电磁波模式的解析解.为了有效确定频率空间域中的电磁场,采用二维Patterson自适应求积算法结合有限连分式展开技术计算傅氏逆变换.最后通过数值模拟结果证明了该算法的有效性,考察了不同各向异性系数、不同井眼倾角以及仪器长度和工作频率变化等情况下的多分量感应测井响应特征.     

基于交联线圈的随钻方位电阻率仪的地层电阻率成像仿真

本文针对具有交联接收天线的随钻方位电阻率仪器在穿过不同地层模型时,对地层模型电阻率成像进行仿真.交联接收天线由极化方向正交的两个线圈组成.交联接收天线在仪器旋转一周内的电压响应可以分离出反映地层背景电阻率信息的响应分量和反映地层电阻率在不同方位上相对变化的响应分量.两种分量进行线性迭加后可以实现对地层电阻率的全方位成像.本文通过对具有交联接收线圈的随钻方位电阻率仪器在两层以及两层以上的复杂地层模型中的电阻率成像仿真,提供了实际仪器在地质导向应用的参考模板,对于现场解释具有重要指导作用.     

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

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

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

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

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

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

层状各向异性地层中含环状天线槽随钻方位电磁波测井几何因子算法

本文基于感应测井中的几何因子与一阶Born近似理论,通过柱坐标系下耦合势Helmholtz方程三维有限体积法研究建立了各向异性地层中含环状天线槽的随钻方位电磁波测井几何因子.首先,通过引入环状电流源与径向磁偶极子源电磁场Green函数,并结合Green第二积分公式,推导出随钻方位电磁波测井仪器中环形电场和横向磁场微小变化与各向异性地层中水平和垂直电导率相对摄动量之间的关系,得到轴向和横向感应电动势微小变化量计算方法与柱坐标系中轴向和横向分量空间灵敏度函数（三维几何因子）表达式.在此基础上,进一步给出径向和纵向微分几何因子计算公式.然后,应用耦合势三维有限体积法确定发射线圈产生的感应电磁场以及环状电流源与径向磁偶极子源电磁场Green函数的数值解.最后,给出不同频率、不同倾角、不同各向异性系数等多种情况下该仪器的空间灵敏度函数与微分几何因子的数值结果,用于分析考查仪器的响应特征与空间探测能力.

     

基于方位电磁波的煤岩层界面及低阻异常体探测方法（英文）

以煤岩介质和煤岩界面作为电磁波传播响应研究对象,构建与煤岩电磁参数相关的地质模型。针对煤矿井下随钻方位电磁波探测技术在煤岩地质导向中的响应特性进行深入分析,研究煤岩界面探测识别方法及其适用性。模拟仪器在煤岩层响应特性规律,建立煤层顶底板界面的方位、仪器到煤层界面距离的计算方法。利用拟牛顿法对随钻方位电磁波电阻率仪器响应理论模型进行反演计算,并对模拟数据数据进行实时反演,降低了Jacobian矩阵的计算次数,提高了反演速度并降低了对初值的依赖性。基于正反演数值仿真模拟,获得低阻含水层、异常地质体（陷落柱、夹矸层等）方位电磁波探测的信号响应规律。通过正交线圈的定向电动势模拟成像方法研究,验证了高阻煤层中方位电磁波电阻率应用的可行性、能够为煤层界面、地质异常体探测及地质导向应用提供依据。     

用有限差分法计算各向异性介质中多分量感应测井的响应

感应测井电阻率是砂泥岩储层中定量评价含油气性的主要参数之一,然而由于储层电各向异性的存在,传统的感应测井仪得到的电阻率主要响应地层的水平电阻率,比地层的垂直电阻率小,由此计算的含水饱和度偏高．多分量感应测井则可以通过不同方向的发射和接收,测量多个分量,弥补上述缺陷,为储层的含油气性评价提供准确信息．本文用交错网格有限差分法导出了主轴各向异性介质中Maxwell方程的离散化关系式,计算了三维各向异性介质中多分量感应测井的电磁响应．利用三层模型的数值模拟,分析了不同磁场分量对各向异性介质的响应特性,用多层模型和倾斜各向异性地层的响应计算检验了所编制的交错网格有限差分法程序．     

用模式匹配算法研究层状各向异性倾斜地层中多分量感应测井响应

利用数值模式匹配技术,研究并建立层状各向异性倾斜地层中多分量感应测井响应的快速数值模拟算法.首先将位于井轴上的三个相互正交磁偶极子转化成关于极角θ的三个谐变分量的叠加,以便将三个正交磁偶极子电磁场的正演问题完全简化成三个谐变分量电磁场的轴对称定解问题,并给出电磁场各个谐变分量在井轴上满足的边界条件,保证电磁场在井轴附近仍然可解.然后利用数值模式匹配技术建立电磁场各个谐变分量的正演过程,得到层状各向异性地层中磁流源并矢Green函数的半解析表达式,给出计算层状各向异性倾斜地层中多分量感应测井响应的具体方法,最后通过数值计算结果证明该算法的有效性并考察几种不同情况下多分量感应测井响应特征.     

用模式匹配算法研究层状各向异性倾斜

利用数值模式匹配技术,研究并建立层状各向异性倾斜地层中多分量感应测井响应的快速数值模拟算法.首先将位于井轴上的三个相互正交磁偶极子转化成关于极角θ的三个谐变分量的叠加,以便将三个正交磁偶极子电磁场的正演问题完全简化成三个谐变分量电磁场的轴对称定解问题,并给出电磁场各个谐变分量在井轴上满足的边界条件,保证电磁场在井轴附近仍然可解.然后利用数值模式匹配技术建立电磁场各个谐变分量的正演过程,得到层状各向异性地层中磁流源并矢Green函数的半解析表达式,给出计算层状各向异性倾斜地层中多分量感应测井响应的具体方法,最后通过数值计算结果证明该算法的有效性并考察几种不同情况下多分量感应测井响应特征.     

3-D rheologic model of earthquake preparation (Ⅱ):Strain field and its applications

On the basis of the three-dimensional elastic inclusion model, the analytic expression of viscoelastic strain field is derived, i.e., the analytic expression of viscoelastic strain at an arbitrary point (x, y, z) in x-axis, y-axis and z-axis produced by three-dimension inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model, namely the normal strains εxx(r, t), εyy(r, t) and εzz(r, t), the shear strains εxy(r, t) and εyx(r, t), εyz(r, t)and εzy(r, t), εxz(r, t) and εzx(r, t), and the bulk-strain θ (r, t). By computing the spatial-temporal variation of bulk strain on the ground produced by a spherical rheologic inclusion in a semi-infinite rheologic medium, we obtained some significant results that the bulk-strain variation with time produced by a hard inclusion has three stages (α, β,γ) with different characteristics, which are similar to those of most geodetic deformation curves, but not the case for those by a soft inclusion. It is meaningful that these theoretical results have been applied to explain preliminarily the characteristics of stage variation of spatial-temporal evolution, the pattern and quadrant distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to found the physical model of earthquake precursors and a reference to predict physically the earthquakes.