CN111177950A - 阵列感应测井曲线环境影响规律的信息处理系统及方法 - Google Patents
阵列感应测井曲线环境影响规律的信息处理系统及方法 Download PDFInfo
- Publication number
- CN111177950A CN111177950A CN202010054506.6A CN202010054506A CN111177950A CN 111177950 A CN111177950 A CN 111177950A CN 202010054506 A CN202010054506 A CN 202010054506A CN 111177950 A CN111177950 A CN 111177950A
- Authority
- CN
- China
- Prior art keywords
- curve
- stratum
- invasion
- borehole
- drilling fluid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000006698 induction Effects 0.000 title claims abstract description 126
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000010365 information processing Effects 0.000 title claims abstract description 29
- 238000005553 drilling Methods 0.000 claims abstract description 119
- 239000012530 fluid Substances 0.000 claims abstract description 118
- 230000009545 invasion Effects 0.000 claims abstract description 115
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 80
- 230000004044 response Effects 0.000 claims abstract description 62
- 239000011435 rock Substances 0.000 claims abstract description 58
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 23
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 239000004576 sand Substances 0.000 claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 9
- 238000011156 evaluation Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims description 39
- 238000001514 detection method Methods 0.000 claims description 36
- 230000007613 environmental effect Effects 0.000 claims description 18
- 238000012545 processing Methods 0.000 claims description 18
- 238000011160 research Methods 0.000 claims description 17
- 239000013505 freshwater Substances 0.000 claims description 15
- 238000003672 processing method Methods 0.000 claims description 11
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims description 10
- 239000011780 sodium chloride Substances 0.000 claims description 10
- 238000012790 confirmation Methods 0.000 claims description 6
- 230000005684 electric field Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000011010 flushing procedure Methods 0.000 claims description 2
- 238000005755 formation reaction Methods 0.000 description 54
- 238000005259 measurement Methods 0.000 description 15
- 239000000243 solution Substances 0.000 description 15
- 238000010586 diagram Methods 0.000 description 12
- 238000004088 simulation Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 11
- 230000002500 effect on skin Effects 0.000 description 10
- 238000003491 array Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000012937 correction Methods 0.000 description 6
- 239000000523 sample Substances 0.000 description 6
- 239000012267 brine Substances 0.000 description 5
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 230000035515 penetration Effects 0.000 description 4
- 238000005316 response function Methods 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 230000033558 biomineral tissue development Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 238000001308 synthesis method Methods 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于信息处理技术领域，公开了一种阵列感应测井曲线环境影响规律的信息处理系统及方法，建立地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围；典型地层的模型建立与响应计算，针对均匀地层、有井眼地层、无侵入地层、有侵入地层以及不同侵入类型的地层中阵列感应测井原始曲线特征研究；砂泥岩剖面中钻井液电阻率和井眼尺寸影响的曲线特征分析；砂泥岩剖面中侵入、围岩影响的曲线特征分析；选取油水识别工作提供依据的合成曲线。本发明可提高阵列感应测井曲线的解释成果的准确性，为地层评价及油气识别提供理论依据。
Description
技术领域
本发明属于信息处理技术领域，尤其涉及一种阵列感应测井曲线环境影响规律的信息处理系统及方法。
背景技术
阵列感应测井是在电子和计算机技术高速发展的基础上实现的感应类电阻率测井方法。早在1957年，Poupon就提出了阵列感应测井和“软件聚焦”的思想，但是，直到二十世纪九十年代才出现商用阵列感应测井仪器。目前，阵列感应测井仪器已经比较成熟，已经有AIT、HDIL、HRAI、HIL和MIT等仪器。上述各种阵列感应测井仪器均已经在中国国内应用，由于它们的分辨率较高和对电阻率剖面的良好反映，已经得到较普遍的认可，因此应用越来越多。但是，通过总结中国国内各油田阵列感应测井资料及应用情况发现，不同研究工区仍然有许多问题需要深入确定和解决，涉及的内容包括盐水井眼问题、钻井液侵入问题和薄层问题等，这些问题严重影响和制约了阵列感应测井的定量应用，需要根据具体情况进行分析或解决。目前，国内外学者针对阵列感应测井响应受到的环境影响做了多方位的研究，江国明等(2005年)，Kruger(2007年)，仵杰等(2016年)完成了井斜对阵列感应测井响应特征的影响分析；夏宏泉等(2017年)完成了鄂尔多斯盆地阵列感应测井响应的钻井液侵入影响分析及校正工作；白彦等(2018年)研究了垂直井中井眼和各向异性地层的阵列感应测井响应特征；仵杰等(2018年)针对井眼、地层、仪器偏心、倾角因素对阵列感应测井原始信号进行了分析。
综合来看，虽然已有学者通过多种环境条件来研究阵列感应测井响应特征，也在研究中取得了较好的效果，但目前尚未有系统性地在各种环境因素条件下综合研究原始曲线和合成曲线特征并结合实际数据处理的分析方法。
目前的阵列感应测井仪器有多种，如AIT(Array Induction Imager Tool)仪器、HRAI(High Resolution Array Induction)仪器和HDIL(High-Definition InductionLog)仪器等。目前，HDIL在中国陆地油田和海洋石油测井中应用较多，效果较好，测井业内普遍认为，HDIL仪器是目前国内应用最好的阵列感应测井仪器。本发明以HDIL仪器为仪器模型进行阵列感应测井方法受环境影响变化规律的说明。HDIL是多道全数字频谱感应测井系统，每个子阵列均采用三线圈系结构(一个发射单元，两个基本接收单元)，HDIL共有7组接收阵列，共用一个发射线圈。使用八种频率(10、30、50、70、90、110、130、150kHz)同时工作，井下仪器共测量112个原始实部分量和虚部分量信号，经计算机处理，实现数字聚焦和井眼校正，得到三种纵向分辨率(1、2、4ft)，六种探测深度(10、20、30、60、90、120in.)的阵列感应测井电阻率曲线。假定地层旋转对称，因而地层电导率仅是纵向和径向深度的函数。地层二维电导率信号经过多个子阵列(响应函数)得到多个输出(各子阵列测量信号)，近似反映二维地层信息。它是一个单输入多输出系统。由于受井眼、侵入、围岩等环境影响和趋肤效应影响以及响应函数有限的分辨率，测量信号不能有效地反映原状地层信息。感应测井信号处理的目的就是消除测量信号中不必要的影响，求出地层真电阻(导)率。因此，在HDIL的信号处理中，环境因素对曲线的处理结果影响很大，若不重视环境因素的影响则会导致阵列感应测井曲线的异常，降低阵列感应测井曲线质量，无法准确识别油气层，进行有效的油气勘探开发工作。阵列感应测井是一种重要的电法测井方法，是电磁波理论应用于在地球物理测井的一种具体技术体现。电磁波运动规律遵从Maxwell方程，研究在普通非均匀介质中的方程解就是感应测井的正演问题。
目前，数值模式匹配法(NMM：Numerical Mode Matching Method)是一种成熟的感应测井响应数值模拟法，已经得到普遍的认可和广泛的应用。其核心是将二维的数值计算转化为一维的解析递推和另一维的数值计算：在轴向上用解析解，在径向上用有限元素法(FEM)求解。
综上所述，现有技术存在的问题是：由于受井眼、侵入、围岩等环境影响和趋肤效应影响以及响应函数有限的分辨率，测量信号不能有效地反映原状地层信息；导致阵列感应测井曲线的异常，阵列感应测井曲线质量降低，无法准确识别油气层，进行有效的油气勘探开发工作。
目前，阵列感应测井的响应特征研究问题主要面临两个难点。第一，合理的地层模型设计和电磁场推导及求解方法的建立；第二，质量可靠的原始曲线选取及合理的合成方法进行分辨率匹配及曲线合成工作。
基于上述难点分析，提出一种阵列感应测井曲线环境影响规律的信息处理系统及方法，可以为选取质量可靠的原始曲线和分辨率匹配、曲线合成工作提供理论依据，从而完成更有效的储层识别，为油气勘探开发奠定基础。
发明内容
针对现有技术存在的问题，本发明提供了一种阵列感应测井曲线环境影响规律的信息处理系统及方法。
本发明是这样实现的，一种阵列感应测井曲线环境影响规律的信息处理方法，所述阵列感应测井曲线环境影响规律的信息处理方法包括以下步骤：
第一步，建立地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围；
第二步，典型地层的模型建立与响应计算，针对均匀地层、有井眼地层、无侵入地层、有侵入地层以及不同侵入类型的地层中阵列感应测井原始曲线特征研究，为后续合成处理提供质量可靠的原始曲线；
第三步，砂泥岩剖面中钻井液电阻率和井眼尺寸影响的曲线特征分析；针对不同地层厚度和不同钻井液电阻率的砂泥岩地层模型中阵列感应测井合成曲线特征研究，评估深浅探测曲线质量，为地层评价工作提供可靠曲线；
第四步，砂泥岩剖面中侵入、围岩影响的曲线特征分析；针对有井眼、无侵入地层模型的围岩影响的阵列感应测井曲线特征研究和有井眼、有侵入地层模型的围岩影响的阵列感应测井合成曲线特征，选取为油水识别工作提供依据的合成曲线。
进一步，所述第一步建立地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围包括：建立存在M个水平界面的地层模型，每个层都是径向不均匀的，由井、冲洗带、侵入带和原状地层组成。采用柱坐标
z为井轴，发射线圈位于层面zS上，用子午面上的一个点(a,zS)代表，a是其半径，发射线圈产生一个轴对称电磁场；发射线圈通过交变电流IT＝I0eiωt，电磁场含有时间因子eiωt，基于经典微分形式的Maxwell方程组可得：
其中，JT为发射电流密度，k2＝ω2μ(ε-iσ/ω)；根据轴对称性可知，E与
无关，对无源区，在第m层内
则：
波数
其中，ω是角频率，μm(r),εm(r),σm(r)分别是第m层的磁导率、介电常数和电导率，径向坐标r的函数；公式
采用分离变量法求解得：
其中，
是N阶对角阵；
建立基函数应满足的边界条件后，对公式
转置变换，公式
中的fm(r)只有数值解，需要通过求解广义复特征值问题来解决；公式
中有解析表达式，得电场分布，计算电导率；
将均匀无限大地层模型电阻率设置为从0.1Ω.m变化到1000Ω.m，观察不同子阵列得到的原始曲线特征，得到阵列感应测井仪器的适用条件。
进一步，所述第二步具体包括：
(1)建立考虑井眼的无限厚地层模型，在不同钻井液条件下对比不同频率的多个子阵列的原始响应信号；
(2)建立考虑井眼、围岩的多层地层模型，在不同层厚围岩条件下对比不同频率的多个子阵列的原始响应信号；
(3)建立考虑井眼、围岩、侵入的多层地层模型，在高侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号；
(4)建立考虑井眼、围岩、侵入的多层地层模型，在低侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号。
进一步，所述第三步具体包括：
(1)建立多层地层模型，设置井眼半径为常用尺寸，钻井液类型从盐水钻井液变化到淡水钻井液，地层厚度从薄到厚变换，确定井眼钻井液影响对阵列感应合成曲线特征；
(2)建立均匀无限大地层模型，在不同钻井液条件下确定，当井眼变化时相同分辨率的不同探测深度曲线特征。
进一步，所述第四步具体包括：
(1)建立盐水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征；
(2)建立淡水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征。
本发明的另一目的在于提供一种实施所述阵列感应测井曲线环境影响规律的信息处理方法的阵列感应测井曲线环境影响规律的信息处理系统，所述阵列感应测井曲线环境影响规律的信息处理系统包括：
地层模型分析模块，用于建立简单的地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围；
模型建立与响应计算模块，用于实现针对典型不同地层模型进行正演模拟；
第一曲线特征分析模块，用于实现砂泥岩剖面中钻井液电阻率和井眼尺寸影响的曲线特征分析；
第二曲线特征分析模块，用于实现砂泥岩剖面中侵入、层厚和围岩影响的曲线特征分析。
进一步，所述模型建立与响应计算模块包括：
第一原始响应信号对比单元，用于建立考虑井眼的无限厚地层模型，在不同钻井液条件下对比不同频率的多个子阵列的原始响应信号；
第二原始响应信号对比单元，用于建立考虑井眼、围岩的多层地层模型，在不同层厚围岩条件下对比不同频率的多个子阵列的原始响应信号；
第三原始响应信号对比单元，用于建立考虑井眼、围岩、侵入的多层地层模型，在高侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号；
第四原始响应信号对比单元，用于建立考虑井眼、围岩、侵入的多层地层模型，在低侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号。
进一步，所述第一曲线特征分析模块包括：
第一阵列感应合成曲线特征确认单元，用于建立多层地层模型，设置井眼半径为常用尺寸，钻井液类型从盐水钻井液变化到淡水钻井液，地层厚度从薄到厚变换，确定井眼钻井液影响对阵列感应合成曲线特征；
探测深度曲线特征确认单元，用于建立均匀无限大地层模型，在不同钻井液条件下确定，当井眼变化时相同分辨率的不同探测深度曲线特征；
第二阵列感应合成曲线特征确认单元，用于建立盐水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征；
第三阵列感应合成曲线特征阵列感应合成曲线特征建立淡水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征。
本发明的另一目的在于提供一种实现所述阵列感应测井曲线环境影响规律的信息处理方法的信息数据处理终端。
本发明的另一目的在于提供一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行所述的阵列感应测井曲线环境影响规律的信息处理方法。
综上所述，本发明的优点及积极效果为：本发明采用数值模式匹配法，建立二维地层模型，进行电磁场分解、波场递推、磁分量求解及电导率计算，对以下几点影响阵列感应测井响应的主要环境因素进行系统分析：1)地层电阻率；2)钻井液电阻率；3)钻井液侵入；4)井眼尺寸；5)地层厚度。建立均匀无限大地层模型，完成不同导电性质地层中的阵列感应测井响应特征分析，确定趋肤效应影响的地层条件，规范仪器适用范围。建立有井眼、无限厚地层模型，将钻井液导电性质从好到差，目的层电阻率从低到高变化，与均匀无限大厚层响应对比，得到井眼影响和钻井液电阻率影响的曲线变化特征；同时对比仪器多组发射频率下的原始信号，分析趋肤效应影响的曲线变化特征，确定质量好的原始曲线类型。建立钻井液侵入半径从浅至深变化的地层模型，钻井液侵入关系可分为高侵、低侵两种状态，分析合成曲线的变化特征，划分质量好的合成曲线，为侵入带反演奠定基础。建立井眼尺寸不同的地层模型，分析阵列感应测井仪器的合成曲线特征，评估深浅探测曲线质量，为地层评价工作提供可靠曲线。建立不同层厚的地层模型，考虑钻井液侵入分为高侵、低侵两种情况，分析各子阵列的原始曲线质量，避免出现不同源距的曲线交叉现象。
本发明利用数值模式匹配法对影响阵列感应测井曲线的主要环境因素进行系统分析的方法，针对复杂地层和井眼条件设计不同的地层模型并进行数值模拟计算，从电磁场分解、波场递推、磁分量求解到电导率计算、分辨率匹配等得到不同地层模型下的阵列感应测井响应曲线，并分析总结在钻井液电阻率、钻井液侵入、井眼尺寸、地层厚度以及地层电阻率中某一因素单独影响下的曲线特征规律。本发明可提高阵列感应测井曲线的解释成果的准确性，为地层评价及油气识别提供理论依据。
附图说明
图1是本发明实施例提供的阵列感应测井曲线环境影响规律的信息处理系统的结构示意图；
图中：1、地层模型分析模块；2、模型建立与响应计算模块；3、第一曲线特征分析模块；4、第二曲线特征分析模块。
图2是本发明实施例提供的阵列感应测井曲线环境影响规律的信息处理方法流程图。
图3是本发明实施例提供的阵列感应测井曲线环境影响分析系统框图。
图4是本发明实施例提供的均匀地层模型(I)HDIL仪器原始信号模拟结果示意图。
图5是本发明实施例提供的有井眼地层模型(II)HDIL仪器原始信号模拟结果(Rm＝10Ω.m)示意图。
图6是本发明实施例提供的有井眼地层模型(II)HDIL仪器原始信号模拟结果(Rm＝100Ω.m)示意图。
图7是本发明实施例提供的无侵入地层模型(III)HDIL响应(对比度10，Rt/Rs＝100/10)示意图。
图8是本发明实施例提供的无侵入地层模型(III)HDIL响应(对比度100，Rt/Rs＝1000/10)示意图。
图9是本发明实施例提供的高侵地层模型(IV)HDIL响应(侵入半径0.25m，Rxo/Rt＝300/80)示意图。
图10是本发明实施例提供的高侵地层模型(IV)HDIL响应(侵入半径1.0m，Rxo/Rt＝300/80)示意图。
图11是本发明实施例提供的低侵地层模型(V)HDIL响应(侵入半径0.25m，Rxo/Rt＝5/20)
图12是本发明实施例提供的低侵地层模型(V)HDIL响应(侵入半径0.5m，Rxo/Rt＝5/20)示意图。
图13是本发明实施例提供的有井眼不同厚度地层模型及HDIL合成曲线(Rm＝0.005Ω.m，Rs＝3.0Ω.m，Rt/Rs＝20/3)示意图。
图14是本发明实施例提供的有井眼不同厚度地层模型及HDIL合成曲线(Rm＝1.0Ω.m，Rs＝3.0Ω.m，Rt/Rs＝20/3)示意图。
图15是本发明实施例提供的井眼尺寸(12in.)对120in.探测深度曲线M1RX的影响示意图。
图16是本发明实施例提供的井眼尺寸(16in.)对120in.探测深度曲线M1RX的影响示意图。
图17是本发明实施例提供的井眼尺寸(12in.)对10in.探测深度曲线M1R1影响示意图。
图18是本发明实施例提供的井眼尺寸(16in.)对10in.探测深度曲线M1R1影响示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
针对现有技术存在的问题，本发明提供了一种阵列感应测井曲线环境影响规律的信息处理系统及方法，下面结合附图对本发明作详细的描述。
如图1所示，本发明实施例提供的阵列感应测井曲线环境影响规律的信息处理系统包括：
地层模型分析模块1，用于建立简单的地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围。
模型建立与响应计算模块2，用于实现针对典型不同地层模型进行正演模拟。
第一曲线特征分析模块3，用于实现砂泥岩剖面中钻井液电阻率和井眼尺寸影响的曲线特征分析。
第二曲线特征分析模块4，用于实现砂泥岩剖面中侵入、围岩影响的曲线特征分析。
如图2所示，本发明实施例提供的阵列感应测井曲线环境影响规律的信息处理方法包括以下步骤：
S201：建立简单的地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围；
S202：典型地层的模型建立与响应计算，针对典型不同地层模型进行正演模拟，主要是针对均匀地层、有井眼地层、无侵入地层、有侵入地层以及不同侵入类型的地层中阵列感应测井原始曲线特征研究，为后续合成处理提供质量可靠的原始曲线；
S203：砂泥岩剖面中钻井液电阻率和井眼尺寸影响的曲线特征分析。针对不同地层厚度和不同钻井液电阻率的砂泥岩地层模型中阵列感应测井合成曲线特征研究，评估深浅探测曲线质量，为地层评价工作提供可靠曲线；
S204：砂泥岩剖面中侵入、和围岩影响的曲线特征分析。针对有井眼、无侵入地层模型的围岩影响的阵列感应测井曲线特征研究和有井眼、有侵入地层模型的围岩影响的阵列感应测井合成曲线特征研究，选取可为油水识别工作提供依据的合成曲线。
本发明实施例提供的阵列感应测井曲线环境影响规律的信息处理方法具体包括以下步骤：
第一步，建立简单的地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围；
考虑轴对称的二维地层模型情况，建立存在M个水平界面的地层模型，每个层都是径向不均匀的，由井、冲洗带、侵入带和原状地层组成。采用柱坐标
z为井轴，发射线圈位于层面zS上，用子午面上的一个点(a,zS)代表，a是其半径，发射线圈产生一个轴对称电磁场。设发射线圈通过交变电流IT＝I0eiωt，电磁场含有时间因子eiωt，基于经典微分形式的Maxwell方程组可得：
其中，JT为发射电流密度，k2＝ω2μ(ε-iσ/ω)。根据轴对称性可知，E与
无关，对无源区，设在第m层内
则：
波数
其中，ω是角频率，μm(r),εm(r),σm(r)分别是第m层的磁导率、介电常数和电导率，它们都是径向坐标r的函数。公式(2)可采用分离变量法求解可得：
其中，
是N阶对角阵。
建立基函数应满足的边界条件后，对公式(3)转置变换，公式(3)中的fm(r)只有数值解，需要通过求解广义复特征值问题来解决。而公式(4)中有解析表达式，可得电场分布，从而计算电导率。
将均匀无限大地层模型电阻率设置为从0.1Ω.m变化到1000Ω.m，观察不同子阵列得到的原始曲线特征，从而得到阵列感应测井仪器的适用条件。
第二步，典型地层的模型建立与响应计算，针对典型不同地层模型进行正演模拟，主要是针对均匀地层、有井眼地层、无侵入地层、有侵入地层以及不同侵入类型的地层中阵列感应测井原始曲线特征研究，为后续合成处理提供质量可靠的原始曲线。
(1)建立考虑井眼的无限厚地层模型，在不同钻井液条件下对比不同频率的多个子阵列的原始响应信号；
(2)建立考虑井眼、围岩的多层地层模型，在不同层厚围岩条件下对比不同频率的多个子阵列的原始响应信号；
(3)建立考虑井眼、围岩、侵入的多层地层模型，在高侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号；
(4)建立考虑井眼、围岩、侵入的多层地层模型，在低侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号。
第三步，砂泥岩剖面中钻井液电阻率和井眼尺寸影响的曲线特征分析。针对不同地层厚度和不同钻井液电阻率的砂泥岩地层模型中阵列感应测井合成曲线特征研究，评估深浅探测曲线质量，为地层评价工作提供可靠曲线。
(1)建立多层地层模型，设置井眼半径为常用尺寸，钻井液类型从盐水钻井液变化到淡水钻井液，地层厚度从薄到厚变换，确定井眼钻井液影响对阵列感应合成曲线特征；
(2)建立均匀无限大地层模型，在不同钻井液条件下确定，当井眼变化时相同分辨率的不同探测深度曲线特征。
第四步，砂泥岩剖面中侵入、围岩影响的曲线特征分析。针对有井眼、无侵入地层模型的围岩影响的阵列感应测井曲线特征研究和有井眼、有侵入地层模型的围岩影响的阵列感应测井合成曲线特征研究，选取可为油水识别工作提供依据的合成曲线。
(1)建立盐水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征；
(2)建立淡水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征。
下面结合附图对本发明的技术方案作进一步的描述。
本发明采用NMM正演方法计算阵列感应HDIL的测井响应，针对不同的影响因素，设计针对性的地层模型，对响应规律进行总结和分析；具体步骤如下：
1、典型地层阵列感应测井原始曲线特征：
图4～图12给出HDIL仪器在多种地层模型条件下的正演模拟结果，主要地层参数已经在对应的图3中给出。针对每一个具体的地层模型，可以把同一源距的不同频率的实部信号模拟结果放在同一道显示，也可以把同一频率的不同源距的实部信号模拟结果放在同一道显示，以满足不同分析的需要。
模拟的原始曲线中包含7组子阵列、八种频率的结果，共有56条实部视电阻率曲线。在图4～图12中，原始曲线名称如下：
A0R1～A0R15：对应子阵列A0，八种频率10kHz、30kHz、…、150kHz的8条实部视电阻率曲线；
其它类推：A1R1～A1R15、A2R1～A2R15、A3R1～A3R15、A4R1～A4R15、A5R1～A5R15和A6R1～A6R15：分别对应子阵列A1～A6，八种频率10kHz、30kHz、…、150kHz的8条实部视电阻率曲线。
此外，需要说明的是，图4～图12中各种模拟响应曲线均为原始响应信号，即未进行井眼影响校正、趋肤效应校正等信号处理过程。
下面结合具体地层模型进行说明分析，表1为设计地层模型的主要信息汇总。
表1设计地层模型的主要信息
(1)图4为均匀地层的HDIL仪器正演模拟结果；可以看出，地层高阻时，地层电阻率与视电阻率基本一致；低阻条件下，视电阻率与地层真实电阻率差别明显，传播效应(趋肤效应)影响明显。
(2)图5和图6为有井眼、无限厚地层的HDIL仪器正演模拟结果，分别为对应Rm＝1.0和100.0Ω.m时两种频率(10kHz和130kHz)条件下的原始响应信号(实部视电阻率)模拟结果。可以看出，在钻井液电阻率较低的条件下，无论地层高阻或低阻，对应的视电阻率与无井眼时的结果严重偏离，同时由于趋肤效应影响，很难进行井眼影响校正。频率越高，趋肤效应影响严重；而钻井液电阻率越大，有井眼地层的阵列感应测井的原始信号与无井眼地层的原始信号图4更为接近，说明钻井液电阻率越高，井眼影响越小，井眼影响校正比较容易，阵列感应测井合成聚焦曲线将更为合理。
(3)图7和图8为无侵地层正演模拟结果。图7中，高阻层电阻率(100.0Ω.m)与围岩电阻率(10.0Ω.m)的对比度为10:1；图8中，高阻层电阻率(1000.0Ω.m)与围岩电阻率(10.0Ω.m)对比度为100:1。
如图7、图8所示，地层越薄，测量结果与实际数值差别越大，说明地层厚度对测量影响明显；无钻井液侵入时，短线圈子阵列能够比较好地反映地层电阻率和地层界面，而长线圈子阵列的测量数值与实际地层差别较大，主要原因在于围岩影响和趋肤效应明显。高对比度(Rt:Rs＝100:1)的地层，围岩影响特别明显，如图8，各子阵列在高阻层与地层模型对应都很差。而且，在该模型的薄层段，长线圈子阵列数值与地层模型对应关系很差。
(5)图9和附图10为高侵地层模型正演模拟结果。高阻层电阻率为80.0Ω.m，侵入带电阻率为300.0Ω.m；低阻层电阻率为80.0Ω.m，无侵入。图9中，高阻层的侵入半径为0.25m；附图11中，高阻层的侵入半径为1.0m。
如图9、图10所示，地层厚度对测量结果影响明显；侵入半径越大，则侵入影响越明显；在钻井液高侵的条件下，各阵列原始曲线的一致性关系良好，在厚层，未出现不同源距曲线的交叉现象；短线圈子阵列反映侵入带电阻率明显；而长线圈子阵列的测量结果受侵入带影响不明显，在厚层可以较好反映原状地层特性。
(6)图11和图12为低侵地层模型正演模拟结果。高阻层电阻率为20.0Ω.m，侵入带电阻率为5.0Ω.m；低阻层电阻率为1.0Ω.m，无侵入。图11中，高阻层的侵入半径为0.25m；图12中，高阻层的侵入半径为0.50m。
如图11、图12所示，地层厚度对测量结果影响明显；侵入半径越大，则侵入影响越明显；而且，在钻井液低侵条件下，即使在厚层，也出现不同源距曲线的交叉现象；短线圈子阵列受侵入带影响明显，主要反映侵入带电阻率；而长线圈子阵列的测量结果同样明显受侵入带影响，使数据处理和分析的难度增大。
2、钻井液电阻率和井眼尺寸影响合成曲线特征
阵列感应测井的主要特点是软件合成聚焦，在合成聚焦之前通常需要对原始信号进行必要的预处理，保证为后续处理提供良好的测量信号。井眼校正是其中非常重要的一个处理环节，如果井眼影响过大导致原始信号严重扭曲，则很难进行井眼校正，从而影响最终的合成聚焦结果。井眼尺寸的影响通常与钻井液电阻率影响结合起来确定。
步骤1)钻井液电阻率影响分析
该组地层模型是井眼半径为4in.，钻井液电阻率从0.005Ω.m变化到1.0Ω.m，地层厚度从0.3m变化到3m；本组模型主要确定井眼钻井液影响。图13和图14为针对不同地层模型的HDIL合成曲线结果。
图13、图14中，第1道：地层模型电阻率参数，Rm为钻井液电阻率，Rt为地层模型的真电阻率；第2道：深度道；第3道至第5道：对应HDIL的1ft、2ft和4ft三种分辨率曲线：M1R1～M1RX、M2R1～M2RX和M4R1～M4RX，探测深度为10in.、20in.、30in.、60in.、90in.和120in.。
当钻井液为盐水钻井液时，如Rm＝0.005Ω.m，阵列感应曲线质量比较差，10in.、20in.即M2R1和M2R2曲线受井眼钻井液影响明显，严重偏离地层模型真实值。30in.、60in.、90in.和120in.即M2R3～M2RX曲线能反应部分地层特性，即由于受钻井液影响，M2R3～M2RX数值与地层真实电阻率数值有明显的差别。
当钻井液为淡水钻井液时，例如，当钻井液电阻率Rm＝1.0Ω.m时，阵列感应曲线质量比较好，各曲线基本与模型值基本重合，一致性较好。
根据数值模拟及合成处理结果可以得到：
(1)在比较淡的钻井液或油基钻井液中，不同探测深度的阵列感应测井曲线均能较好地反映地层的真实信息。因此，在淡水或油基钻井液钻井条件下，应用阵列感应测井是有效的，这与以往的结论完全一致。
(2)随着钻井液电阻率的降低(即钻井液矿化度增高)，测量结果与地层真实数值差别越来越大，即钻井液电阻率越小(钻井液越咸，矿化度越大)，对HDIL合成测井曲线质量影响越大，阵列感应测井反映地层的能力减弱。主要原因在于：井眼钻井液影响明显，趋肤效应明显，探测深度降低，地层真实信息贡献降低，从而影响合成曲线质量。
(3)钻井液电阻率较低或地层与钻井液电阻率对比度较大时，即使没有钻井液侵入，但是不同探测深度的曲线仍然会出现差异，这属于“假侵入”现象，需要引起重视。
(4)阵列感应测井深探测曲线对原状地层的反映能力通常要好于浅探测曲线，这与一般的理论和理解是吻合的。但是，根据图13的可以看出，即使钻井液电阻率很小(Rm分别为0.005Ω.m)，深探测曲线(如M2R9和M2RX)仍然能够较大程度地反映地层真实信息。因此，在盐水钻井液条件下，可以谨慎和有条件地使用阵列感应测井资料，尤其是其深探测曲线仍然可有效反映地层真实信息。
(5)在薄层处，阵列感应测井除受井眼钻井液影响外，还要受围岩影响，因此，合成曲线反映地层的能力变弱；同时，在Rt与Rm接近时，即在低对比度条件下，合成曲线的质量较好，尤其是低阻地层，数据更为合理。
步骤2)井眼尺寸影响分析
井眼尺寸也是影响测井曲线质量的重要因素。下面给出当井眼变化时的地层模型及其合成处理结果。图15、图16和图17、图18给出12in.和16in.井眼对120in.探测深度曲线(M1RX)和10in.探测深度曲线(M1R1)的影响特征。图15-图18中的横坐标为地层真电阻率Rt，纵坐标为合成的视电阻率与真电阻率的比值Ra/Rt。
结合前面的处理结果，根据图15、图16和图17、图18，可以有如下结论：
(1)井眼越大，对测量结果越不利，尤其是盐水钻井液条件的大井眼，对HDIL的测量极为不利，使得HDIL仪器的应用范围明显变小；
(2)Rt/Rm的大小(即地层与钻井液电阻率对比度)对测量结果影响很大，若井眼增大，则影响更为明显；
(3)浅探测曲线较深探测曲线受井眼影响更为明显，因此，深探测曲线有更大的应用范围。
3、砂泥岩剖面中侵入、围岩影响的合成曲线特征。
围岩及钻井液侵入是影响阵列感应测井重要因素，通过设计不同钻井液侵入地层模型的模拟算例，分析围岩和钻井液侵入对HDIL响应的影响，综合分析HDIL的综合反映地层的能力，研究HDIL仪器的适用性原则和地层真实参数的求解方法。
阵列感应测井信号要合成多个探测深度曲线，需要各子阵列原始信号探测深度要与合成的理想探测深度存在合理的关系。HDIL的合成曲线探测深度均在线圈系原始径向探测深度之间，两相邻子阵列的探测深度之间仅有一个合成探测深度。
步骤1)盐水钻井液条件下考虑井眼、侵入地层模型，给出有井眼、侵入和围岩影响的地层模型以及HDIL合成曲线结果，主要考虑钻井液侵入对阵列感应测井的影响。
HDIL合成曲线结果对应盐水钻井液(Rm＝0.06Ω.m)和不同侵入半径，侵入半径依次对应为0.1m到3m。可以发现，当钻井液侵入比较浅时，10in.探测深度曲线受钻井液影响比较严重，其它探测深度曲线能够反应真实地层情况；随着钻井液侵入半径的增加，HDIL合成曲线受钻井液侵入影响明显，反应冲洗带情况的能力增强，与Rt数值偏离明显。
步骤2)淡水钻井液条件下考虑井眼、有侵入地层模型
HDIL合成曲线结果对应淡水钻井液(Rm＝0.4Ω.m)和不同侵入半径，侵入半径依次对应为0.1m到3m。可以发现，当没有钻井液侵入时，所有曲线与地层模型基本重合，说明淡水钻井液情况下，HDIL合成曲线能很好地探测到真实地层情况，这与前面的结论一致。随着钻井液侵入半径的增加，探测深度为10in.(M2R1)和20in.(M2R2)的浅探测曲线主要反应钻井液和侵入带信息，其它中、深探测曲线也受钻井液侵入影响。
在钻井液侵入比较浅时，合成曲线能够反应地层真实情况；当钻井液侵入较深时，合成曲线要反应地层真实情况比较困难，需对原始曲线进行针对性校正，采用合理的处理合成方法，获取可靠的阵列感应测井曲线，完成精确的地层评价及油水识别工作。
下面结合实验对本发明的技术效果作详细的描述。
HDIL合成曲线结果对应盐水泥浆侵入水层模型，钻井液电阻率值为0.05Ω.m，侵入带深度为2.0m，冲洗带电阻率为2.0Ω.m，目的层电阻率为1.0Ω.m，该情况较为不理想测量环境条件下的探测结果。依据上述总结的原始信号特征，选取合适的合成方法得到质量可靠的合成曲线后，选取浅探测曲线用作侵入带信息确定，分辨率小的数据作薄层信息处理，最终处理得到的目的层电阻率及侵入带电阻率信息，其精度可达80％以上。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种阵列感应测井曲线环境影响规律的信息处理方法，其特征在于，所述阵列感应测井曲线环境影响规律的信息处理方法包括以下步骤：
第一步，建立地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围；
第二步，典型地层的模型建立与响应计算，针对均匀地层、有井眼地层、无侵入地层、有侵入地层以及不同侵入类型的地层中阵列感应测井原始曲线特征研究，为后续合成处理提供原始曲线；
第三步，砂泥岩剖面中钻井液电阻率和井眼尺寸影响的曲线特征分析；针对不同地层厚度和不同钻井液电阻率的砂泥岩地层模型中阵列感应测井合成曲线特征研究，评估深浅探测曲线质量，为地层评价工作提供可靠曲线；
第四步，砂泥岩剖面中侵入、围岩影响的曲线特征分析；针对有井眼、无侵入地层模型的围岩影响的阵列感应测井曲线特征研究和有井眼、有侵入地层模型的围岩影响的阵列感应测井合成曲线特征，选取为油水识别工作提供依据的合成曲线。
2.如权利要求1所述的阵列感应测井曲线环境影响规律的信息处理方法，其特征在于，所述第一步建立地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围包括：建立存在M个水平界面的地层模型，每个层都是径向不均匀的，由井、冲洗带、侵入带和原状地层组成；采用柱坐标
z为井轴，发射线圈位于层面zS上，用子午面上的一个点(a,zS)代表，a是其半径，发射线圈产生一个轴对称电磁场；发射线圈通过交变电流IT＝I0eiωt，电磁场含有时间因子eiωt，基于经典微分形式的Maxwell方程组可得：
其中，JT为发射电流密度，k2＝ω2μ(ε-iσ/ω)；根据轴对称性可知，E与
无关，对无源区，在第m层内
则：
波数
其中，ω是角频率，μm(r),εm(r),σm(r)分别是第m层的磁导率、介电常数和电导率，径向坐标r的函数；公式
采用分离变量法求解得：
其中，
是N阶对角阵；
建立基函数应满足的边界条件后，对公式
转置变换，公式
中的fm(r)只有数值解，需要通过求解广义复特征值问题来解决；公式
中有解析表达式，得电场分布，计算电导率；
将均匀无限大地层模型电阻率设置为从0.1Ω.m变化到1000Ω.m，观察不同子阵列得到的原始曲线特征，得到阵列感应测井仪器的适用条件。
3.如权利要求1所述的阵列感应测井曲线环境影响规律的信息处理方法，其特征在于，所述第二步具体包括：
(1)建立考虑井眼的无限厚地层模型，在不同钻井液条件下对比不同频率的多个子阵列的原始响应信号；
(2)建立考虑井眼、围岩的多层地层模型，在不同层厚围岩条件下对比不同频率的多个子阵列的原始响应信号；
(3)建立考虑井眼、围岩、侵入的多层地层模型，在高侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号；
(4)建立考虑井眼、围岩、侵入的多层地层模型，在低侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号。
4.如权利要求1所述的阵列感应测井曲线环境影响规律的信息处理方法，其特征在于，所述第三步具体包括：
(1)建立多层地层模型，设置井眼半径为常用尺寸，钻井液类型从盐水钻井液变化到淡水钻井液，地层厚度从薄到厚变换，确定井眼钻井液影响对阵列感应合成曲线特征；
(2)建立均匀无限大地层模型，在不同钻井液条件下确定，当井眼变化时相同分辨率的不同探测深度曲线特征。
5.如权利要求1所述的阵列感应测井曲线环境影响规律的信息处理方法，其特征在于，所述第四步具体包括：
(1)建立盐水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征；
(2)建立淡水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征。
6.一种实施权利要求1～5任意一项所述阵列感应测井曲线环境影响规律的信息处理方法的阵列感应测井曲线环境影响规律的信息处理系统，其特征在于，所述阵列感应测井曲线环境影响规律的信息处理系统包括：
地层模型分析模块，用于建立简单的地层模型，基于数值模式匹配法，完成进行电磁场分解、波场递推、磁分量求解及电导率计算、曲线合成，分析仪器的适用范围；
模型建立与响应计算模块，用于实现针对典型不同地层模型进行正演模拟；
第一曲线特征分析模块，用于实现砂泥岩剖面中钻井液电阻率和井眼尺寸影响的曲线特征分析；
第二曲线特征分析模块，用于实现砂泥岩剖面中侵入、围岩影响的曲线特征分析。
7.如权利要求6所述的阵列感应测井曲线环境影响规律的信息处理系统，其特征在于，所述模型建立与响应计算模块包括：
第一原始响应信号对比单元，用于建立考虑井眼的无限厚地层模型，在不同钻井液条件下对比不同频率的多个子阵列的原始响应信号；
第二原始响应信号对比单元，用于建立考虑井眼、围岩的多层地层模型，在不同层厚围岩条件下对比不同频率的多个子阵列的原始响应信号；
第三原始响应信号对比单元，用于建立考虑井眼、围岩、侵入的多层地层模型，在高侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号；
第四原始响应信号对比单元，用于建立考虑井眼、围岩、侵入的多层地层模型，在低侵条件下对比不同径向侵入深度时的不同频率的多个子阵列的原始响应信号。
8.如权利要求6所述的阵列感应测井曲线环境影响规律的信息处理系统，其特征在于，所述第一曲线特征分析模块包括：
第一阵列感应合成曲线特征确认单元，用于建立多层地层模型，设置井眼半径为常用尺寸，钻井液类型从盐水钻井液变化到淡水钻井液，地层厚度从薄到厚变换，确定井眼钻井液影响对阵列感应合成曲线特征；
探测深度曲线特征确认单元，用于建立均匀无限大地层模型，在不同钻井液条件下确定，当井眼变化时相同分辨率的不同探测深度曲线特征；
第二阵列感应合成曲线特征确认单元，用于建立盐水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征；
第三阵列感应合成曲线特征阵列感应合成曲线特征建立淡水泥浆条件下，设置井眼半径为常用尺寸，在低阻围岩条件下考虑侵入的三层地层模型中确定井眼钻井液影响对阵列感应合成曲线特征。
9.一种实现权利要求1～5任意一项所述阵列感应测井曲线环境影响规律的信息处理方法的信息数据处理终端。
10.一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如权利要求1-5任意一项所述的阵列感应测井曲线环境影响规律的信息处理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010054506.6A CN111177950B (zh) | 2020-01-17 | 2020-01-17 | 阵列感应测井曲线环境影响规律的信息处理系统及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010054506.6A CN111177950B (zh) | 2020-01-17 | 2020-01-17 | 阵列感应测井曲线环境影响规律的信息处理系统及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111177950A true CN111177950A (zh) | 2020-05-19 |
| CN111177950B CN111177950B (zh) | 2022-06-21 |
Family
ID=70646667
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010054506.6A Active CN111177950B (zh) | 2020-01-17 | 2020-01-17 | 阵列感应测井曲线环境影响规律的信息处理系统及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111177950B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112034522A (zh) * | 2020-08-21 | 2020-12-04 | 中石化石油工程技术服务有限公司 | 一种六子阵列线圈系测量地层电阻率的方法 |
| CN113312796A (zh) * | 2021-06-22 | 2021-08-27 | 成都理工大学 | 考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 |
| CN114647912A (zh) * | 2022-02-28 | 2022-06-21 | 电子科技大学 | 一种过套管水泥侵入模型建立方法及过套管地层三孔隙度测量方法 |
| CN116084912A (zh) * | 2021-11-08 | 2023-05-09 | 中国石油天然气集团有限公司 | 多深度测井仪器信号标定装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20010039477A1 (en) * | 1999-03-08 | 2001-11-08 | Jaiqi Xiao | Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well |
| CN102789003A (zh) * | 2012-08-15 | 2012-11-21 | 中国石油天然气股份有限公司 | 利用阵列感应测井资料进行地层参数反演的方法及装置 |
| CN103437760A (zh) * | 2013-08-30 | 2013-12-11 | 中国石油天然气股份有限公司 | 一种使用阵列感应数据快速评价油水层的方法 |
| CN107451313A (zh) * | 2016-06-01 | 2017-12-08 | 中国石油化工股份有限公司 | 一种基于自适应模型的地层泥浆侵入判断方法及系统 |
-
2020
- 2020-01-17 CN CN202010054506.6A patent/CN111177950B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20010039477A1 (en) * | 1999-03-08 | 2001-11-08 | Jaiqi Xiao | Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well |
| CN102789003A (zh) * | 2012-08-15 | 2012-11-21 | 中国石油天然气股份有限公司 | 利用阵列感应测井资料进行地层参数反演的方法及装置 |
| CN103437760A (zh) * | 2013-08-30 | 2013-12-11 | 中国石油天然气股份有限公司 | 一种使用阵列感应数据快速评价油水层的方法 |
| CN107451313A (zh) * | 2016-06-01 | 2017-12-08 | 中国石油化工股份有限公司 | 一种基于自适应模型的地层泥浆侵入判断方法及系统 |
Non-Patent Citations (5)
| Title |
|---|
| GAO JIE等: "Inversion of array induction logs and its application", 《PETROLEUM SCIENCE》 * |
| 仵杰等: "阵列侧向测井仪器的正演响应分析", 《西安石油大学学报(自然科学版)》 * |
| 王熊等: "阵列感应测井的影响因素分析", 《国外测井技术》 * |
| 范宜仁等: "高分辨率阵列侧向测井响应数值模拟", 《测井技术》 * |
| 高杰等: "大斜度井侧向测井三维正演数值模拟及曲线快速校正方法研究", 《石油勘探与开发》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112034522A (zh) * | 2020-08-21 | 2020-12-04 | 中石化石油工程技术服务有限公司 | 一种六子阵列线圈系测量地层电阻率的方法 |
| CN112034522B (zh) * | 2020-08-21 | 2023-09-19 | 中石化石油工程技术服务有限公司 | 一种六子阵列线圈系测量地层电阻率的方法 |
| CN113312796A (zh) * | 2021-06-22 | 2021-08-27 | 成都理工大学 | 考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 |
| CN113312796B (zh) * | 2021-06-22 | 2022-06-03 | 成都理工大学 | 考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 |
| CN116084912A (zh) * | 2021-11-08 | 2023-05-09 | 中国石油天然气集团有限公司 | 多深度测井仪器信号标定装置 |
| CN114647912A (zh) * | 2022-02-28 | 2022-06-21 | 电子科技大学 | 一种过套管水泥侵入模型建立方法及过套管地层三孔隙度测量方法 |
| CN114647912B (zh) * | 2022-02-28 | 2023-03-28 | 电子科技大学 | 一种过套管水泥侵入模型建立方法及过套管地层三孔隙度测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111177950B (zh) | 2022-06-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111177950B (zh) | 阵列感应测井曲线环境影响规律的信息处理系统及方法 | |
| Wu et al. | A new ultra-deep azimuthal electromagnetic LWD sensor for reservoir insight | |
| Li et al. | New directional electromagnetic tool for proactive geosteering and accurate formation evaluation while drilling | |
| US7359800B2 (en) | Determination of fracture orientation and length using multi-component and multi-array induction data | |
| Rosthal et al. | Field test results of an experimental fully-triaxial induction tool | |
| CN102042009B (zh) | 一种测量地层电阻率的阵列感应测井线圈系 | |
| Barber et al. | Determining formation resistivity anisotropy in the presence of invasion | |
| US10451763B2 (en) | Evaluation of formation utilizing wideband electromagnetic measurements | |
| US6381542B1 (en) | Generic, accurate, and real time borehole correction for resistivity tools | |
| Thiel et al. | 2D reservoir imaging using deep directional resistivity measurements | |
| da Rocha et al. | Permeability estimation and analysis of fracture networks using resistivity logs in an offshore Aptian carbonate reservoir pre-salt, in the Southeastern Santos Basin | |
| Noh et al. | Real-time 2.5 D inversion of LWD resistivity measurements using deep learning for geosteering applications across faulted formations | |
| Rabinovich et al. | Effect of fractures on multi-component and multiarray induction logs | |
| Nan et al. | Numerical simulation, response analysis, and physical experiment of induction logging in an inclined fractured formation | |
| Mikhaylov et al. | Joint inversion of induction and galvanic logging data in axisymmetric geological models | |
| Hou et al. | Characterization of formation fractures with multicomponent induction logging based on biaxial anisotropy models: Method and case studies | |
| Hou et al. | Real-time borehole correction for a new multicomponent array induction logging tool in OBM wells | |
| CN107939385B (zh) | 定量计算极化值及应用的方法 | |
| Gao et al. | SVR-based weighted processing method for electrical imaging logging in oil-based mud | |
| Qin et al. | Application of conventional propagation resistivity logging for formation boundary identification in geosteering | |
| Zhdanov et al. | Fast 3D imaging from a single borehole using tensor induction logging data | |
| Shetty et al. | New multiphysics, multiscale inversion for imaging petrophysical properties in anisotropic, laminated formations | |
| Yang et al. | Bed-boundary effect removal to aid formation resistivity interpretation from LWD propagation measurements at all dip angles | |
| Nan et al. | Influential factors of inclined fracture induction logging based on numerical simulation and physical experiment | |
| Hou et al. | Correction for the borehole effect of multi-component array induction log data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |