CN110488367B - 一种基于阵列侧向测井资料的电阻率反演初值选取方法 - Google Patents
一种基于阵列侧向测井资料的电阻率反演初值选取方法 Download PDFInfo
- Publication number
- CN110488367B CN110488367B CN201910786024.7A CN201910786024A CN110488367B CN 110488367 B CN110488367 B CN 110488367B CN 201910786024 A CN201910786024 A CN 201910786024A CN 110488367 B CN110488367 B CN 110488367B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- initial value
- mlr4
- mlr1
- logging
- 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.)
- Active
Links
- 238000010187 selection method Methods 0.000 title claims abstract description 14
- 238000011010 flushing procedure Methods 0.000 claims abstract description 21
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 238000012417 linear regression Methods 0.000 claims abstract description 4
- 238000012937 correction Methods 0.000 claims description 7
- 239000011435 rock Substances 0.000 claims description 7
- 101000605076 Homo sapiens Ligand-dependent nuclear receptor corepressor-like protein Proteins 0.000 abstract description 4
- 102100038259 Ligand-dependent nuclear receptor corepressor-like protein Human genes 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 230000001502 supplementing effect Effects 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 description 14
- 230000009545 invasion Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000009977 dual effect Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000006698 induction Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000004005 microsphere Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 238000004088 simulation Methods 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Developing Agents For Electrophotography (AREA)
Abstract
本发明公开一种基于阵列侧向测井资料的电阻率反演初值选取方法,包括以下步骤:基于径向阶跃地层模型进行正演计算,定义电阻率幅度差系数S1、S2、S3、S4、S5、S6,获取径向阶跃地层模型针对的阵列侧向测井中MLR1a、MLR2a、MLR3a、MLR4a四条基础测井曲线的视电阻率,建立电阻率幅度差系数随冲洗带半径的变化图版;基于电阻率幅度差系数随冲洗带半径的变化图版,通过线性回归建立冲洗带半径初值rxocz与电阻率幅度差系数的定量关系;将阵列侧向测井提供的针对基础测井曲线MLR4a的视电阻率作为地层电阻率初值Rtcz,然后基于伪几何因子定义得到冲洗带电阻率初值Rxocz,该方法无需补充其他测井资料即可有效减小初始值与真实值间的偏差,有助于提高电阻率测井反演精度和速度。
Description
技术领域
本发明涉及油气资源勘探开发中的测井评价技术领域,尤其涉及一种基于阵列侧向测井资料的电阻率反演初值选取方法。
背景技术
储层电阻率资料在油层评价方面具有重要用途,电阻率不仅可以用来定性识别油气层(一般油层显现出高阻特征),也可以用以定量计算含油饱和度(利用阿尔奇或印度尼西亚等公式)。多种电阻率测井方法均在油田现场得到了广泛应用,如双侧向测井、双感应测井、阵列侧向测井、阵列感应测井、电成像测井等,随着油气藏评价的日益精细化,阵列型电阻率测井仪器相较传统双侧向侧向、双感应测井仪器可以提供多条纵向分辨率更高的视电阻率曲线,从而越来越多地被油田现场所应用,相应的阵列型电阻率测井资料的处理工作也日益成为亟待解决的重点问题。
侧向测井仪器响应受泥浆侵入、井眼、围岩等因素的影响使得测得的视电阻率曲线不能真实反映储层电阻率,需要对其进行反演处理,恢复泥侵剖面电阻率信息。目前在油田现场对阵列侧向测井资料的处理主要分为两种。一种是直接将视电阻率当作储层真实电阻率来使用,这种方法在储层泥浆侵入不深时评价误差在可接受范围内,但在泥侵比较严重的层位会造成严重误差。另一种处理方法是利用微球或微侧向测井近似估计侵入带电阻率,然后利用泥浆侵入校正图版对视电阻率进行校正,得到侵入带半径和地层电阻率,这种方法需要进行大量正演计算构建数据库,然后从数据库中查找校正图版,工作量一般很大,同时校正结果直接取决于正演计算的精度。
随着计算机技术的发展,人们将反演思想引入电测井反演领域。电测井反演的思想就是先构建电阻率测井反演模型,通过相应的反演算法进行反复迭代,当正演算子和实际值的误差小于预先设定的阈值,则将此时对应的正演模型参数当做反演参数输出。通过对现场实时测量的测井数据进行快速反演处理来消除泥浆侵入影响,可以初步恢复地层原始特征,为后期精细解释提供基础数据,同时也是测井人员判断测井资料好坏的重要依据。马奎特算法是一种油田现场运用最为广泛的电阻率反演方法。该反演算法收敛速度极快,利于测井资料的实时反演,但在实际运用过程中这类算法的弊端也暴露了出来,那就是反演结果受初值影响极大,若初值设置不合理(偏离真实值较远),则很难准确反演出地层真实电阻率值,严重影响反演精度。为解决马奎特算法存在的问题,有学者将模拟退火法、遗传算法、人工神经网络算法、免疫算法、差分进化算法、粒子群等全局寻优算法引入了电测井资料的反演处理中,此类算法具备较强的全局搜索能力,可以有效提高反演精度,但收敛速度普遍较慢,不适用于电测井资料的快速处理,使用受限,马奎特算法仍是现场阵列侧向测井资料实时反演的首要选择。
发明内容
本发明要解决的技术问题在于,提供一种基于阵列侧向测井资料的电阻率反演初值选取方法。
本发明解决其技术问题所采用的技术方案是:构造一种基于阵列侧向测井资料的电阻率反演初值选取方法,包括以下步骤:
步骤S100:基于径向阶跃地层模型进行正演计算,定义电阻率幅度差系数S1、S2、S3、S4、S5、S6,获取所述径向阶跃地层模型针对的阵列侧向测井中MLR1a、MLR2a、MLR3a、MLR4a四条基础测井曲线的视电阻率,并进行井眼、围岩校正,建立电阻率幅度差系数随冲洗带半径的变化图版;
步骤S200:基于所述电阻率幅度差系数随冲洗带半径的变化图版,通过线性回归建立冲洗带半径初值rxocz与电阻率幅度差系数的定量关系;
步骤S300:将阵列侧向测井提供的针对基础测井曲线MLR4a的视电阻率作为地层电阻率初值Rtcz,然后基于伪几何因子定义得到冲洗带电阻率初值Rxocz。
优选地,在所述步骤S100中,S1、S2、S3、S4、S5、S6分别定义如下:
S1=2·(MLR4a-MLR1a)/(MLR4a+MLR1a) (1)
S2=2·(MLR4a-MLR2a)/(MLR4a+MLR2a) (2)
S3=2·(MLR4a-MLR3a)/(MLR4a+MLR3a) (3)
S4=2·(MLR3a-MLR1a)/(MLR3a+MLR1a) (4)
S5=2·(MLR3a-MLR2a)/(MLR3a+MLR2a) (5)
S6=2·(MLR2a-MLR1a)/(MLR2a+MLR1a) (6)
其中,MLR1a、MLR2a、MLR3a、MLR4a代表阵列侧向测井提供的针对四条基础测井曲线的视电阻率,单位为Ω·m;S1、S2、S3、S4、S5、S6无量纲。
优选地,在步骤S200中,建立的冲洗带半径初值rxocz与电阻率幅度差系数的定量关系为:
rxocz=925.14·S1+1004.04·S2-1835.46·S3-1295.69·S4-221.40·S5+70.02·S6+1.39 (7)。
优选地,在步骤S300中,将阵列侧向测井提供的针对基础测井曲线MLR4a的视电阻率作为地层电阻率初值Rtcz,然后基于伪几何因子定义得到冲洗带电阻率初值Rxocz:
其中,K1为MLR1a对应的仪器常数;r0为井眼半径;h1为MLR1a对应的主电流厚度;K1、r0、h1均为已知值。
优选地,还包括
步骤S400:将所述冲洗带电阻率初值Rxocz与实际地层参数对比。
实施本发明具有以下有益效果:本发明的基于阵列侧向测井资料的电阻率反演初值选取方法无需补充其他测井资料即可有效减小初始值与真实值间的偏差,有助于提高电阻率测井反演精度和速度。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明基于阵列侧向测井资料的电阻率反演初值选取方法的示意图;
图2是本发明径向阶跃地层模型的示意图;
图3是本发明电阻率幅度差系数随冲洗带半径的变化图版的示意图。
具体实施方式
如图1所示,是本发明的基于阵列侧向测井资料的电阻率反演初值选取方法,解决目前利用马奎特算法进行阵列侧向测井反演初始值选取不合适的不足,有利于提高反演精度和算法收敛速度。
包括以下步骤:
步骤S100:基于径向阶跃地层模型进行正演计算,定义电阻率幅度差系数S1、S2、S3、S4、S5、S6,获取径向阶跃地层模型针对的阵列侧向测井中MLR1a、MLR2a、MLR3a、MLR4a四条基础测井曲线的视电阻率,并进行井眼、围岩校正,建立电阻率幅度差系数随冲洗带半径的变化图版。其中,如图2所示,为径向阶跃地层模型。如图3所示,为电阻率幅度差系数随冲洗带半径的变化图版。
在该步骤S100中,S1、S2、S3、S4、S5、S6分别定义如下:
S1=2·(MLR4a-MLR1a)/(MLR4a+MLR1a) (1)
S2=2·(MLR4a-MLR2a)/(MLR4a+MLR2a) (2)
S3=2·(MLR4a-MLR3a)/(MLR4a+MLR3a) (3)
S4=2·(MLR3a-MLR1a)/(MLR3a+MLR1a) (4)
S5=2·(MLR3a-MLR2a)/(MLR3a+MLR2a) (5)
S6=2·(MLR2a-MLR1a)/(MLR2a+MLR1a) (6)
其中,MLR1a、MLR2a、MLR3a、MLR4a代表阵列侧向测井提供的针对四条基础测井曲线的视电阻率,单位为Ω·m;S1、S2、S3、S4、S5、S6无量纲。
进一步的,获取径向阶跃地层模型(理论模型)针对的MLR1a、MLR2a、MLR3a、MLR4a基础测井的视电阻率,并进行井眼、围岩校正,如表1。
表1理论模型阵列侧向测井响应
针对阵列侧向测井响应MLR1a、MLR2a、MLR3a、MLR4a,优选相应的井眼、围岩校正图版,将图版读值拟合成校正公式的基础上进行数值模拟计算机自动校正,消除井眼、围岩对阵列侧向测井响应的影响,以免其他因素对电阻率反演初值选取精确程度产生影响。
步骤S200:基于电阻率幅度差系数随冲洗带半径的变化图版,通过线性回归建立冲洗带半径初值rxocz与电阻率幅度差系数的定量关系。
在该步骤S200中,建立的冲洗带半径初值rxocz与电阻率幅度差系数的定量关系为:
rxocz=925.14·S1+1004.04·S2-1835.46·S3-1295.69·S4-221.40·S5+70.02·S6+1.39 (7)。
其中,rxocz为冲洗带半径初值,in。
步骤S300:将阵列侧向测井提供的针对基础测井曲线MLR4a的视电阻率作为地层电阻率初值Rtcz,然后基于伪几何因子定义得到冲洗带电阻率初值Rxocz。
在该步骤S300中,将阵列侧向测井提供的针对基础测井曲线MLR4a的视电阻率作为地层电阻率初值Rtcz,然后基于伪几何因子定义得到冲洗带电阻率初值Rxocz:
其中,K1为MLR1a对应的仪器常数;r0为井眼半径;h1为MLR1a对应的主电流厚度;K1、r0、h1均为已知值。
在本实施例中,还包括
步骤S400:将冲洗带电阻率初值Rxocz与实际地层参数对比。可以理解的,反演初值与实际地层参数对比,验证基于阵列侧向测井资料的电阻率反演初值选取方法的有效性和优越性,对比结果如表2。
表2反演初值与实际地层参数对比
经两种反演初值选取方法的对比,基于阵列侧向测井资料的电阻率测井初始值选取方法效果更好,可推广应用到现场。
本发明的基于阵列侧向测井资料的电阻率反演初值选取方法无需补充其他测井资料即可有效减小初始值与真实值间的偏差,有助于提高电阻率测井反演精度和速度。
可以理解的,以上实施例仅表达了本发明的优选实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制;应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,可以对上述技术特点进行自由组合,还可以做出若干变形和改进,这些都属于本发明的保护范围;因此,凡跟本发明权利要求范围所做的等同变换与修饰,均应属于本发明权利要求的涵盖范围。
Claims (5)
1.一种基于阵列侧向测井资料的电阻率反演初值选取方法,其特征在于,包括以下步骤:
步骤S100:基于径向阶跃地层模型进行正演计算,定义电阻率幅度差系数S1、S2、S3、S4、S5、S6,获取所述径向阶跃地层模型针对的阵列侧向测井中MLR1a、MLR2a、MLR3a、MLR4a四条基础测井曲线的视电阻率,并进行井眼、围岩校正,建立电阻率幅度差系数随冲洗带半径的变化图版;
步骤S200:基于所述电阻率幅度差系数随冲洗带半径的变化图版,通过线性回归建立冲洗带半径初值rxocz与电阻率幅度差系数的定量关系;
步骤S300:将阵列侧向测井提供的针对基础测井曲线MLR4a的视电阻率作为地层电阻率初值Rtcz,然后基于伪几何因子定义得到冲洗带电阻率初值Rxocz。
2.根据权利要求1所述的基于阵列侧向测井资料的电阻率反演初值选取方法,其特征在于,在所述步骤S100中,S1、S2、S3、S4、S5、S6分别定义如下:
S1=2·(MLR4a-MLR1a)/(MLR4a+MLR1a) (1)
S2=2·(MLR4a-MLR2a)/(MLR4a+MLR2a) (2)
S3=2·(MLR4a-MLR3a)/(MLR4a+MLR3a) (3)
S4=2·(MLR3a-MLR1a)/(MLR3a+MLR1a) (4)
S5=2·(MLR3a-MLR2a)/(MLR3a+MLR2a) (5)
S6=2·(MLR2a-MLR1a)/(MLR2a+MLR1a) (6)
其中,MLR1a、MLR2a、MLR3a、MLR4a代表阵列侧向测井提供的针对四条基础测井曲线的视电阻率,单位为Ω·m;S1、S2、S3、S4、S5、S6无量纲。
3.根据权利要求2所述的基于阵列侧向测井资料的电阻率反演初值选取方法,其特征在于,在步骤S200中,建立的冲洗带半径初值rxocz与电阻率幅度差系数的定量关系为:
rxocz=925.14·S1+1004.04·S2-1835.46·S3-1295.69·S4-221.40·S5+70.02·S6+1.39 (7)。
5.根据权利要求1所述的基于阵列侧向测井资料的电阻率反演初值选取方法,其特征在于,还包括
步骤S400:将所述冲洗带电阻率初值Rxocz与实际地层参数对比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910786024.7A CN110488367B (zh) | 2019-08-23 | 2019-08-23 | 一种基于阵列侧向测井资料的电阻率反演初值选取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910786024.7A CN110488367B (zh) | 2019-08-23 | 2019-08-23 | 一种基于阵列侧向测井资料的电阻率反演初值选取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110488367A CN110488367A (zh) | 2019-11-22 |
CN110488367B true CN110488367B (zh) | 2021-04-13 |
Family
ID=68553495
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910786024.7A Active CN110488367B (zh) | 2019-08-23 | 2019-08-23 | 一种基于阵列侧向测井资料的电阻率反演初值选取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110488367B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111305834B (zh) * | 2020-02-25 | 2022-09-23 | 北京工业大学 | 基于多探测模式电阻率测井的三维反演初始模型构建方法 |
CN112253090B (zh) * | 2020-10-14 | 2022-10-18 | 中海油田服务股份有限公司 | 一种多频电成像的数据参数反演方法和装置 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1188690A1 (ru) * | 1984-05-14 | 1985-10-30 | Всесоюзный научно-исследовательский институт нефтепромысловой геофизики | Устройство дл бокового каротажа скважин |
RU2402047C1 (ru) * | 2009-09-08 | 2010-10-20 | Владимир Алексеевич Королев | Способ измерения при боковом каротаже скважин |
RU2478223C1 (ru) * | 2011-09-28 | 2013-03-27 | Закрытое акционерное общество Научно-производственная фирма "ГИТАС" (ЗАО НПФ "ГИТАС") | Способ оценки удельного электрического сопротивления пласта при проведении исследований скважин, обсаженных металлической колонной |
CN102678102B (zh) * | 2012-03-31 | 2015-02-25 | 中国石油大学(华东) | 基于阵列电成像测井的储层油水识别方法及识别系统 |
CN103821506B (zh) * | 2014-03-04 | 2016-04-20 | 中国石油大学(华东) | 井周介质电阻率三维阵列成像测量方法 |
CN108729898B (zh) * | 2017-04-17 | 2023-11-28 | 中国石油集团长城钻探工程有限公司 | 侧向测井仪电极系以及侧向测井仪 |
CN109577963A (zh) * | 2018-10-19 | 2019-04-05 | 中国石油天然气股份有限公司 | 模拟各向异性地层阵列侧向测井响应的装置及其构建方法 |
CN109885927B (zh) * | 2019-02-18 | 2022-08-02 | 中国石油天然气集团有限公司 | 一种地层径向电阻率连续反演方法 |
-
2019
- 2019-08-23 CN CN201910786024.7A patent/CN110488367B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110488367A (zh) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11118441B2 (en) | Optimized geosteering using real-time geological models | |
CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
EP2810101B1 (en) | Improving efficiency of pixel-based inversion algorithms | |
CN108345719B (zh) | 一种大斜度井/水平井随钻电磁波测井资料钻后反演方法 | |
CN106154322B (zh) | 测井曲线校正方法和装置 | |
CN108005646B (zh) | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 | |
CN110488367B (zh) | 一种基于阵列侧向测井资料的电阻率反演初值选取方法 | |
CN107218033B (zh) | 一种识别地层产状的方法和地层参数的反演方法 | |
DE112016005364T5 (de) | Verfahren zum Auswählen eines Erdmodells von einer Vielzahl von Erdmodellen | |
US11841476B2 (en) | Methods and systems of determining parameters characterizing porous media from data gathered by a plurality of different tools | |
Zhou et al. | Interpretation of azimuthal propagation resistivity measurements: Modeling, inversion, application and discussion | |
CN109885927B (zh) | 一种地层径向电阻率连续反演方法 | |
CN108665545B (zh) | 录井参数三维地质模型建立方法 | |
CN111965712B (zh) | 一种校正可控源音频大地电磁法静态效应的方法 | |
CN112114374B (zh) | 复杂地质体的多密度界面反演方法 | |
CN112379462A (zh) | 电磁地震数据联合处理方法及装置 | |
Wu et al. | Enhanced reservoir geosteering and geomapping from refined models of ultra-deep LWD resistivity inversions using machine-learning algorithms | |
CN115263287A (zh) | 一种高精度地热资源综合勘查方法 | |
CN110513105A (zh) | 一种三维阵列感应测井仪器数据处理系统及方法 | |
CN115267927B (zh) | 一种基于蚁群-梯度串联算法的多边界幕式地质导向方法 | |
CN117950065B (zh) | 一种水平井阵列电阻率测井资料层厚校正方法 | |
CN113468727B (zh) | 基于先验结构和已知点双重约束的层状密度建模方法 | |
CN114109349B (zh) | 一种致密砂岩储层孔隙度指数/饱和度指数的确定方法 | |
Boesing et al. | High-Angle Formation Response Modeling for More Accurate Petrophysics–Challenges and Experiences from Multiple Cases | |
CN118011509A (zh) | 地层裂缝有效性评价方法及系统 |
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 |