CN102678102B - 基于阵列电成像测井的储层油水识别方法及识别系统 - Google Patents
基于阵列电成像测井的储层油水识别方法及识别系统 Download PDFInfo
- Publication number
- CN102678102B CN102678102B CN201210090672.7A CN201210090672A CN102678102B CN 102678102 B CN102678102 B CN 102678102B CN 201210090672 A CN201210090672 A CN 201210090672A CN 102678102 B CN102678102 B CN 102678102B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- formation
- module
- logging
- parameter
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Developing Agents For Electrophotography (AREA)
Abstract
本发明涉及一种基于阵列电成像测井的储层油水识别方法及识别系统,针对泥浆侵入地层电阻率分布特点,建立“五参数地层模型”,适应了模拟不同泥浆类型、不同侵入地层电阻率分布特点要求,在阵列电成像测井正反演模型和处理方法方面具有开创性;基于阵列电成像测井“五参数反演”方法,进行泥浆侵入地层电阻率分布剖面重构,反映了泥浆侵入油层和水层电阻率分布特征或变化规律差异的实际不同,该方法具有较强的实用性,在储层流体性质定性识别和定量计算的实际应用方面具有开创性。
Description
技术领域
本发明是油气资源地球物理勘查技术领域,主要用于泥浆侵入地层电阻率分布重构,进而进行孔隙流体识别的识别方法和识别系统。
背景技术
测井学科是地球科学的一个重要领域,是油气田勘探开发的重要工程技术之一。在油气田勘探开发中,利用测井资料进行储层流体识别、储层参数计算以及产能评价具有十分重要的意义,其中电测井是目前油气评价最重要的核心技术。但实际电测井提供的视电阻率信息和地层真实情况差异很大,特别是淡水泥浆侵入油层的电阻率分布特征和规律与水层存在实际不同,但往往不能被测井实测视电阻率直观显示,因此,进行电测井反演,进行相关影响因素校正,恢复地层真实情况,对储层油水识别具有实际应用意义,但现有的“图版校正”、“三参数反演”和“线性四参数反演”和实际地层情况差异较大,致使反演效果存在很多不确定性。
1、基于视电阻率高低、深浅电阻率幅度差异定性识别储层流体性质存在不确定性,因为复杂油气层的电阻率信息直观差异较小,甚至出现油、水层电性特征倒挂现象;
2、基于常规“图版校正”方法,地层模型简单,手工工作量大,处理解释经验性要求高;
3、基于常规双感应/侧向的“三参数反演”方法,地层模型为台阶状二段式地层模型,和实际泥浆侵入地层电阻率连续分布差异太大;
4、基于阵列电测井的“四参数反演”方法,地层模型虽为三段式地层模型,但侵入带部分为线性单调变化模型,不满足淡水泥浆侵入油层电阻率非线性、非单调分布模拟的要求;
发明内容
本发明的目的在于提供一套阵列电成像测井的储层油水识别系统及识别方法,用于常规储层/复杂储层油气识别和评价。
本发明的技术方案为:一种基于阵列电成像测井的储层油水识别系统,所述的识别系统包括地层选取模块、地层参数值预设模块、电阻率分布模块、多层阵列感应测井/阵列侧向测 井响应模块,测井实测值模块、判断模块、修正模块、地层参数输出模块和油水识别模块;
首先,地层选取模块对处理井段多地层精细划分选取,然后,测井实测值模块利用测井仪对地层电阻率信息进行检测,实测阵列感应/阵列侧向测井曲线特征值选取,得到测井实测值;
同时,地层参数值预设模块对地层选取模块选取的地层进行地层参数预设:冲洗带半径rxo、过渡带半径ri、冲洗带电阻率Rxo、过渡带中部地层电阻率Ri0.5以及原状地层电阻率Rt;
电阻率分布模块根据地层参数值预设模块对预设的地层参数形成三段式地层电阻率分布模型,多层阵列感应测井/阵列侧向测井响应模块利用三段式地层电阻率分布模型和预设的地层参数计算多层阵列感应测井/阵列侧向测井响应,得到计算测井响应;
判断模块判断计算测井响应与测井实测值是否一致,若一致,输出模块输出真实的地层参数,油水识别模块利用真实的地层参数对三段式泥浆侵入剖面重构,进行油水识别;若不一致,修正模块对预设的地层参数进行修正,直至计算测井响应与测井实测值一致得到真实的地层参数,输出模块输出真实的地层参数,油水识别模块利用真实的地层参数对三段式泥浆侵入剖面重构,进行油水识别。
一种基于阵列电成像测井的储层油水识别方法,包括如下步骤:
①处理井段多地层精细划分选取,确定径向深度r;
②对地层参数初值预设:冲洗带半径rxo、过渡带半径ri、冲洗带电阻率Rxo、过渡带中部地层电阻率Ri0.5以及原状地层电阻率Rt;
形成参数向量形式:x=(x1,x2,x3,x4,x5)T (1)
③建立三段式地层电阻率分布模型:
冲洗带:Rxo,r<rxo
过渡带:Ri=ar2+br+c,rxo<r<ri
原状地层:Rt,r>ri
其中:过渡带电阻率系数确定如下式:
ri0.5代表侵入带一半;Ri0.5代表该点处电阻率
④根据地层电阻率分布,计算阵列电阻率测井响应
fi=fi(x),i=1,......,5 (2)
⑤确定5×5的Frechet导数阵
⑥确定实际测井值yi与地层参数向量(x)计算的测井值f(x)之差构成的差值向量
B=(y1-f1,y2-f2,y3-f3,y4-f4,y5-f5)T (4)
⑦求解线性方程组
AΔx=B (5)
⑧修正地层参数
x=x+Δx (6)
⑨确定测井值与计算值平方差
将(6)式回代,重复步骤③~⑨,直到目标函数Φ为最小为止,这时计算的电阻率和测井的电阻率基本一致,得到真实的地层参数,从而利用真实的地层参数对三段式泥浆侵入剖面重构,进行油水识别。
本发明的有益效果为:针对泥浆侵入地层电阻率分布特点,建立“五参数地层模型”,适应了模拟不同泥浆类型、不同侵入地层电阻率分布特点要求,在阵列电成像测井正反演模型和处理方法方面具有开创性;基于阵列电成像测井“五参数反演”方法,进行泥浆侵入地层电阻率分布剖面重构,反映了泥浆侵入油层和水层电阻率分布特征或变化规律差异的实际不同,该方法具有较强的实用性,在储层流体性质定性识别和定量计算的实际应用方面具有开创性。
附图说明
图1本发明基于阵列电成像测井的储层油水识别系统及识别方法流程图。
图2为本发明对阵列感应资料自动分层与特征值读取示意图。
图3为本发明反演重构阵列感应曲线与实测曲线对比图。
图4为本发明地层电阻率分布剖面。
图5为本发明对阵列侧向资料处理成果示意图。
具体实施方式
下面结合附图说明本发明的具体实施方式:
附图1为基于阵列电成像测井的储层油水识别系统及识别方法流程图。其主要由1测井层段精细划分;2测井曲线特征值读取;3地层参数初值确定;4模拟地层电测井响应;5地层参数调整;6地层五参数输出;7地层电阻率剖面重构组成。测井层段精细划分主要实现基于阵列电测井资料自动分层,确定层界面与层厚;测井曲线特征值读取实现对已分层阵列电测井资料特征值读取;地层参数初值确定主要依据层厚与地层电阻率特征值初步估计地层参数;模拟地层阵列电测井响应采用几何因子法与有限元法实现对感应测井与侧向测井的快速正演;地层参数调整实现基于实测曲线与模拟曲线对比产生迭代步长;地层五参数输出实现对满足精度的地层参数反演结果输出;地层电阻率剖面重构实现基于五参数的地层径向电阻率剖面重构显示。
附图2为本发明对阵列感应资料进行自动分层与特征值自动提取,其中8为自动分层曲线。本发明可选择针对不同分辨率(1ft、2ft、4ft)阵列感应测井曲线进行分层,并输出自动分层曲线;曲线特征值自动提取,可实现对阵列电测井资料响应特征值逐层提取,并生成特征值表格。同时,根据测井特征值产生地层参数预测初值:根据不同探测深度曲线产生地层侵入带、过渡带与原状地层电阻率,并根据曲线差异与特征生成冲洗带半径、过渡带半径初值。
附图3为本发明基于五参数模型重构阵列感应曲线与实测曲线对比图。由五参数初值,采用如下方法确定地层电阻率:
冲洗带:Rxo,r<rxo
过渡带:Ri=ar2+br+c,rxo<r<ri
原状地层:Rt,r>ri
并采用几何因子理论与有限元理论对阵列感应、阵列侧向曲线进行正演重构,并将重构曲线与实测曲线进行对比,判断反演结果是否满足精度,若满足精度则进入6进行参数成果输出与地层电阻率剖面重构,若不满足,则采用公式(3)-(6)式确定参数迭代步长,重新迭代求解。
附图4为五参数反演参数输出。本发明和输出某层电阻率径向分布剖面,其中9为由反 演后五参数确定的地层电阻率分布曲线;10为反演重构该层阵列感应(侧向)响应;11为实测响应。基于五参数重构的地层电阻率分布剖面能直观的显示储层电阻率分布的非线性变化,五参数模型能很好的兼容“三参数”、“四参数”模型,如图4(a)显示了储层“低阻环”特征,而图4(b)中过渡带电阻率为线性变化。同时,成果输出中实测特征值与重构特征值的对比可以直观显示反演结果的精度。
附图5为本发明对阵列侧向处理结果。图5(a)为实测曲线与自动分层数据;图5(b)为五参数成果输出与储层径向重构电阻率剖面。
Claims (2)
1.一种基于阵列电成像测井的储层油水识别方法,其特征在于:包括如下步骤:
(1)处理井段多地层精细划分选取,确定径向深度r;
(2)对地层参数初值预设:冲洗带半径rxo、过渡带半径ri、冲洗带电阻率Rxo、过渡带中部地层电阻率Ri0.5以及原状地层电阻率Rt;
形成参数向量形式:x=(x1,x2,x3,x4,x5)T (1)
(3)建立三段式地层电阻率分布模型:
冲洗带:Rxo,r<rxo
过渡带:Ri=ar2+br+c,rxo<r<ri
原状地层:Rt,r>ri
其中:过渡带电阻率系数确定如下式:
(4)根据地层电阻率分布,计算阵列电阻率测井响应
fi=fi(x),i=1,......,5 (2)
(5)确定5×5的Frechet导数阵
(6)确定实际测井值yi与地层参数向量(x)计算的测井值f(x)之差构成的差值向量
B=(y1-f1,y2-f2,y3-f3,y4-f4,y5-f5)T (3)
(7)求解线性方程组
AΔx=B (3)
(8)修正地层参数
x=x+Δx (4)
(9)确定测井值与计算值平方差
将(4)式回代,重复步骤(3)~(9),直到目标函数Φ为最小为止,这时计算的电阻率和测井的电阻率基本一致,得到真实的地层参数,从而利用真实的地层参数对三段式泥浆侵入剖面重构,进行油水识别。
2.一种基于阵列电成像测井的储层油水识别系统,其特征在于:所述的识别系统包括地层选取模块、地层参数值预设模块、电阻率分布模块、多层阵列感应测井/阵列侧向测井响应模块,测井实测值模块、判断模块、修正模块、地层参数输出模块和油水识别模块;
首先,地层选取模块对处理井段多地层精细划分选取,然后,测井实测值模块利用测井仪对地层电阻率信息进行检测,实测阵列感应/阵列侧向测井曲线特征值选取,得到测井实测值;
同时,地层参数值预设模块对地层选取模块选取的地层进行地层参数预设:冲洗带半径rxo、过渡带半径ri、冲洗带电阻率Rxo、过渡带中部地层电阻率Ri0.5以及原状地层电阻率Rt;
电阻率分布模块根据地层参数值预设模块对预设的地层参数形成三段式地层电阻率分布模型,多层阵列感应测井/阵列侧向测井响应模块利用三段式地层电阻率分布模型和预设的地层参数计算多层阵列感应测井/阵列侧向测井响应,得到计算测井响应;
判断模块判断计算测井响应与测井实测值是否一致,若一致,输出模块输出真实的地层参数,油水识别模块利用真实的地层参数对三段式泥浆侵入剖面重构,进行油水识别;若不一致,修正模块对预设的地层参数进行修正,直至计算测井响应与测井实测值一致得到真实的地层参数,输出模块输出真实的地层参数,油水识别模块利用真实的地层参数对三段式泥浆侵入剖面重构,进行油水识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210090672.7A CN102678102B (zh) | 2012-03-31 | 2012-03-31 | 基于阵列电成像测井的储层油水识别方法及识别系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210090672.7A CN102678102B (zh) | 2012-03-31 | 2012-03-31 | 基于阵列电成像测井的储层油水识别方法及识别系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102678102A CN102678102A (zh) | 2012-09-19 |
CN102678102B true CN102678102B (zh) | 2015-02-25 |
Family
ID=46810606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210090672.7A Active CN102678102B (zh) | 2012-03-31 | 2012-03-31 | 基于阵列电成像测井的储层油水识别方法及识别系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102678102B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9746574B2 (en) | 2016-01-21 | 2017-08-29 | Baker Hughes Incorporated | Resistivity imager for conductive and non-conductive mud |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103821506B (zh) * | 2014-03-04 | 2016-04-20 | 中国石油大学(华东) | 井周介质电阻率三维阵列成像测量方法 |
CN107451313B (zh) * | 2016-06-01 | 2020-10-02 | 中国石油化工股份有限公司 | 一种基于自适应模型的地层泥浆侵入判断方法及系统 |
CN108021726A (zh) * | 2016-10-28 | 2018-05-11 | 中国石油天然气股份有限公司 | 储层有效性的评价方法和装置 |
CN108252707B (zh) * | 2017-11-30 | 2024-04-12 | 杭州迅美科技有限公司 | 一种电成像测井图像增强显示处理方法 |
CN109614584B (zh) * | 2018-11-16 | 2020-10-27 | 中国科学院计算技术研究所 | 一种资源测井数据的重建方法 |
CN109815636B (zh) * | 2019-03-01 | 2023-07-21 | 同济大学 | 一种可用于回转粘度仪的水泥净浆流变性迭代算法 |
CN110488367B (zh) * | 2019-08-23 | 2021-04-13 | 中海石油(中国)有限公司深圳分公司 | 一种基于阵列侧向测井资料的电阻率反演初值选取方法 |
CN111042808B (zh) * | 2019-12-27 | 2021-03-16 | 中国石油天然气集团有限公司 | 一种用于层状砂砾岩储层的电成像图像刻度方法及系统 |
CN112415614B (zh) * | 2020-11-13 | 2021-08-06 | 重庆科技学院 | 一种测井探测尺度三维精细地层模型构建方法 |
CN113740928B (zh) * | 2021-08-16 | 2023-09-15 | 中海油田服务股份有限公司 | 一种致密砂岩储层的流体识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1209499A (zh) * | 1997-03-19 | 1999-03-03 | 江汉石油管理局测井研究所 | 阵列电极系成象测井方法 |
US5889729A (en) * | 1996-09-30 | 1999-03-30 | Western Atlas International, Inc. | Well logging data interpretation systems and methods |
US6594584B1 (en) * | 1999-10-21 | 2003-07-15 | Schlumberger Technology Corporation | Method for calculating a distance between a well logging instrument and a formation boundary by inversion processing measurements from the logging instrument |
CA2380382C (en) * | 1999-08-05 | 2010-06-01 | Schlumberger Canada Limited | A method of determining parameters of formations through which a borehole passes |
CN101775981A (zh) * | 2009-01-09 | 2010-07-14 | 中国石油天然气集团公司 | 一种确定地层真电阻率的方法 |
CN101899973A (zh) * | 2010-05-10 | 2010-12-01 | 中国石油天然气集团公司 | 一种地层水电阻率的测量方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7755361B2 (en) * | 2004-07-14 | 2010-07-13 | Schlumberger Technology Corporation | Apparatus and system for well placement and reservoir characterization |
-
2012
- 2012-03-31 CN CN201210090672.7A patent/CN102678102B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5889729A (en) * | 1996-09-30 | 1999-03-30 | Western Atlas International, Inc. | Well logging data interpretation systems and methods |
CN1209499A (zh) * | 1997-03-19 | 1999-03-03 | 江汉石油管理局测井研究所 | 阵列电极系成象测井方法 |
CA2380382C (en) * | 1999-08-05 | 2010-06-01 | Schlumberger Canada Limited | A method of determining parameters of formations through which a borehole passes |
US6594584B1 (en) * | 1999-10-21 | 2003-07-15 | Schlumberger Technology Corporation | Method for calculating a distance between a well logging instrument and a formation boundary by inversion processing measurements from the logging instrument |
CN101775981A (zh) * | 2009-01-09 | 2010-07-14 | 中国石油天然气集团公司 | 一种确定地层真电阻率的方法 |
CN101899973A (zh) * | 2010-05-10 | 2010-12-01 | 中国石油天然气集团公司 | 一种地层水电阻率的测量方法及装置 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9746574B2 (en) | 2016-01-21 | 2017-08-29 | Baker Hughes Incorporated | Resistivity imager for conductive and non-conductive mud |
Also Published As
Publication number | Publication date |
---|---|
CN102678102A (zh) | 2012-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102678102B (zh) | 基于阵列电成像测井的储层油水识别方法及识别系统 | |
Cardiff et al. | A potential‐based inversion of unconfined steady‐state hydraulic tomography | |
CN103149589B (zh) | 火成岩油气勘探方法及装置 | |
CN100429531C (zh) | 目标最小化的三维电磁快速反演方法 | |
CN102042011B (zh) | 利用常规测井资料构建伪核磁t2谱的方法 | |
BR112019018396B1 (pt) | Método e sistema para executar uma operação de perfuração em uma formação de terra | |
CN104514552A (zh) | 一种煤层气储层识别与丰度预测的方法 | |
Iravani et al. | Application of Seismic Attribute Technique to estimate the 3D model of Hydraulic Flow Units: A case study of a gas field in Iran | |
CN104727813A (zh) | 一种源储共生型地层的孔隙度测量方法 | |
CN109885927B (zh) | 一种地层径向电阻率连续反演方法 | |
CN109667576B (zh) | 一种高矿化度成因低阻油层测井识别方法 | |
CN106154342B (zh) | 一种确定洞穴充填物电阻率的方法 | |
CN1164651A (zh) | 地震-测井多信息储层参数反演方法 | |
Danilovskiy et al. | Deep-learning-based noniterative 2D-inversion of unfocused lateral logs | |
Wen et al. | Simultaneous integration of pressure, water cut, 1 and 4-D seismic data in geostatistical reservoir modeling | |
CN113775326A (zh) | 可动水饱和度的评价方法、装置、电子设备及介质 | |
Sunmonu et al. | Reservoir characterization and by-passed pay analysis of Philus Field in Niger Delta, Nigeria | |
Rastegarnia et al. | Estimation of flow zone indicator distribution by using seismic data: a case study from a central Iranian oilfield | |
Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis | |
Hajizadeh et al. | Designing a pragmatic solution for complex numerical modeling problem in thinly laminated reservoirs | |
Buddo et al. | Joint interpretation of seismic and TEM data from the Kovykta gas-condensate field, East Siberia | |
CN116146170A (zh) | 一种三维感应各向异性泥页岩储层饱和度定量计算方法和装置 | |
Feng et al. | Evaluation of reservoir permeability using array induction logging | |
Fajana et al. | Comparison of modified Waxman-Smith algorithms and Archie models in prospectivity analysis of saturations in shaly-sand reservoirs. A case study of Pennay field, Niger-Delta | |
Dell’Aversana et al. | Quantitative estimation of oil saturation from marine CSEM data: a case history |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |