CN107784159B - 一种储层电阻率各向异性系数的确定方法 - Google Patents
一种储层电阻率各向异性系数的确定方法 Download PDFInfo
- Publication number
- CN107784159B CN107784159B CN201710848308.5A CN201710848308A CN107784159B CN 107784159 B CN107784159 B CN 107784159B CN 201710848308 A CN201710848308 A CN 201710848308A CN 107784159 B CN107784159 B CN 107784159B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- reservoir
- coefficient
- value
- natural gamma
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Mining & Mineral Resources (AREA)
- General Engineering & Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种储层电阻率各向异性系数的确定方法，包括：首先选取研究区具有自然伽马、电阻率、地层倾角测井资料的井，获得纯泥岩段电阻率值与测量井段的地层倾角；选取目的层不同粒径岩心开展筛析粒度实验，结合自然伽马测井资料建立泥质含量计算模型；结合泥质含量、水平电阻率、垂直电阻率、各向异性系数、地层视电阻率模型公式构建方程组；结合电阻率、地层倾角资料，令各向异性系数初值为1.0，代入方程组迭代求解储层各向异性系数。本发明在实际应用中，实现了在没有岩心电阻率各向异性实验的基础上，利用测井资料与理论模型，逐点迭代求解储层电阻率各向异性系数，为储层电阻率准确校正提供基础参数。
Description
技术领域
本发明属于石油勘探中的复杂构造油气藏测井资料校正领域，涉及一种储层电阻率各向异性系数的确定方法。
背景技术
目前，储层电阻率各向异性系数的确定方法有两种，一种是通过岩心实验测量得到某一固定深度点的各向异性系数，另一种是通过三轴感应电阻率测井仪器测量计算得到储层各向异性系数。岩心实验方法是利用目的层的取心资料进行岩心取样，平行于岩心层理方向测量岩样的电阻率为水平电阻率，垂直于岩心层理方向测量岩样的电阻率为垂直电阻率，通过下式计算各向异性系数：
式中：Rh为储层水平电阻率；Rv为储层垂直电阻率；λ为储层各向异性系数。三轴感应电阻率测井仪器(以Rtscanner为例)包括一个三轴发射器，三个用于井筒校正的短间距单轴接收器，六个三轴接受器，利用高级反演技术从测量到的张量数据中提取到电阻率各向异性、层界位置和相对倾角，进而确定储层的垂直与水平电阻率。
以上两种技术方法存在以下问题，通过岩心实验测量只能得到某一深度点的各向异性系数，无法得到储层连续变化的各向异性系数及未取心井段的各向异性系数，该方法受到限制较多；三轴感应电阻率测井可以准确测量储层连续变化的各向异性系数、垂直电阻率与水平电阻率，但该方法测井费用较高，在一些老井中无法应用，不能充分利用老井中已测量的电阻率资料进行各向异性系数求取。
发明内容
为解决现有技术中存在的上述缺陷，本发明的目的在于提供一种储层电阻率各向异性系数的确定方法，实现了在没有岩心电阻率各向异性实验的基础上，利用测井资料与理论模型，逐点迭代求解储层电阻率各向异性系数，为储层电阻率准确校正提供基础参数。
本发明是通过下述技术方案来实现的。
步骤101，选取研究区具有自然伽马、电阻率和地层倾角测井资料的井，获得纯泥岩段电阻率值与测量井段的地层倾角；
步骤102，选取目的层不同粒径岩心开展筛析粒度实验，结合自然伽马测井资料建立泥质含量计算模型；
步骤103，结合泥质含量、水平电阻率、垂直电阻率、地层视电阻率和各向异性系数模型公式构建方程组；
步骤104，结合电阻率和地层倾角测井资料，令各向异性系数初值为1.0，代入方程组迭代求解储层各向异性系数。
进一步地，所述步骤101，选取研究区具有自然伽马、电阻率、地层倾角测井资料的井，通过测井资料获得纯泥岩段电阻率值与测量井段的地层倾角，确定求解方程组中的定量参数。
进一步地，所述步骤102，选取目的层不同粒径岩心开展筛析粒度实验，结合自然伽马测井资料建立泥质含量计算模型如下：
Vsh＝10(a×ΔGR+b) (1)
式中：Vsh为储层泥质含量；ΔGR为自然伽马相对值；GR为自然伽马测井值；GRmin为纯砂岩自然伽马测井值；GRmax为纯泥岩段自然伽马测井值；a、b为模型公式中的系数，通过筛析粒度实验分析数据与测井数据，采用最小二乘法拟合获得。
进一步地，所述的筛析粒度实验按照《岩心分析方法SY/T 5434-2009》标准规定的流程进行。
进一步地，所述步骤103，水平电阻率、垂直电阻率、地层视电阻率、各向异性系数模型公式如下：
式中：Rh为储层水平电阻率；Rsd为储层砂岩电阻率；Rsh为纯泥岩段电阻率；Rv为储层垂直电阻率；λ为储层各向异性系数；Ra为地层视电阻率，即电阻率测井值；θ为储层的地层倾角。
进一步地，所述步骤104，将纯泥岩电阻率、地层倾角、测井电阻率代入式(1)、式(3)、式(4)、式(5)、式(6)组成的方程组，各向异性系数初值为1.0，在迭代过程中逐渐增加各向异性系数值，通过迭代算法求解方程组的解，确定储层电阻率各向异性系数。
本发明的有益效果在于，实现了在没有岩心电阻率各向异性实验的基础上，利用自然伽马、电阻率、地层倾角测井资料，在构建的理论模型方程组的基础上，逐点迭代求解方程组确定储层电阻率各向异性系数，能够更加快速准确确定储层电阻率各向异性系数，为储层电阻率准确校正提供基础参数。
附图说明
图1为本发明提供的一种储层电阻率各向异性系数的确定方法流程图；
图2为本发明实施例中提供的泥质含量计算模型图；
图3为本发明实施例中提供的某口井各向异性系数计算成果图。
具体实施方式
下面通过对本实施例的具体实施情况做进一步详细说明，以支持本发明所要解决的技术问题，但并不作为对发明做任何限制的依据。
如图1所示，本发明实施例提供的一种基于随钻测井的水平井水淹层解释方法，包括如下步骤：
步骤101：选取研究区具有自然伽马、电阻率和地层倾角测井资料的井，获得纯泥岩段电阻率值与测量井段的地层倾角；
步骤102：选取目的层不同粒径岩心开展筛析粒度实验，结合自然伽马测井资料建立泥质含量计算模型；
步骤103：结合泥质含量、水平电阻率、垂直电阻率、地层视电阻率和各向异性系数模型公式构建方程组；
步骤104：结合电阻率、地层倾角测井资料，令各向异性系数初值为1.0，代入方程组迭代求解储层各向异性系数。
下面，通过对本实施例的具体实施情况做进一步详细说明，以支持本发明所要解决的技术问题，按照以下步骤进行操作：
步骤一，选取研究区具有自然伽马、电阻率、地层倾角测井资料的XX207井，得到纯泥岩段电阻率值为2.4Ω·m与测量井段的地层倾角为40度。
步骤二，选取目的层不同粒径岩心开展筛析粒度实验，结合自然伽马测井资料通过最小二乘法拟合获得泥质含量计算模型(如图2所示)：
Vsh＝10(1.2094×ΔGR+0.6435) (1)
式中：Vsh为储层泥质含量；ΔGR为自然伽马相对值；GR为自然伽马测井值；GRmin为纯砂岩自然伽马测井值；GRmax为纯泥岩段自然伽马测井值。
步骤三，水平电阻率、垂直电阻率、地层视电阻率、各向异性系数模型公式如下：
式中：Rh为储层水平电阻率；Rsd为储层砂岩电阻率；Rsh为纯泥岩段电阻率；Rv为储层垂直电阻率；λ为储层各向异性系数；Ra为地层视电阻率，即电阻率测井值；θ为储层的地层倾角。
步骤四，将XX207井的纯泥岩电阻率、地层倾角、测井电阻率代入式(1)、式(3)、式(4)、式(5)、式(6)组成的方程组，各向异性系数初值为1.0，在迭代过程中逐渐增加各向异性系数值，通过迭代算法求解方程组的解，确定储层电阻率各向异性系数。
将上述方法进行编程，实现方法的处理模块化，具体处理结果如图3所示。从图中可以看出，利用本方法通过阵列感应测井M2RX确定的储层电阻率各向异性系数与斯伦贝谢RTscanner测井确定的各向异性系数及岩石物理实验分析的各向异性系数基本一致，验证了本方法的可靠性。
最后所应说明的是，以上具体实施方式仅用以说明本发明的技术方案而非限制，尽管参照实例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1.一种储层电阻率各向异性系数的确定方法，其特征在于，包括以下步骤：
步骤101，选取研究区具有自然伽马、电阻率和地层倾角测井资料的井，获得纯泥岩段电阻率值为2.4Ω·m与测量井段的地层倾角为40度；
步骤102，选取目的层不同粒径岩心开展筛析粒度实验，结合自然伽马测井资料建立泥质含量计算模型；
选取目的层不同粒径岩心开展筛析粒度实验，结合自然伽马测井资料建立泥质含量计算模型如下：
Vsh＝10(a×ΔGR+b) (1)
式中：Vsh为储层泥质含量；ΔGR为自然伽马相对值；GR为自然伽马测井值；GRmin为纯砂岩自然伽马测井值；GRmax为纯泥岩段自然伽马测井值；a、b为模型公式中的系数，通过筛析粒度实验分析数据与测井数据，采用最小二乘法拟合获得；
步骤103，结合泥质含量、水平电阻率、垂直电阻率、地层视电阻率和各向异性系数模型公式构建方程组；
水平电阻率、垂直电阻率、地层视电阻率、各向异性系数模型公式如下：
式中：Rh为储层水平电阻率；Vsh为储层泥质含量；Rsd为储层砂岩电阻率；Rsh为纯泥岩段电阻率；Rv为储层垂直电阻率；λ为储层各向异性系数；Ra为地层视电阻率，即电阻率测井值；θ为储层的地层倾角；
步骤104，结合电阻率和地层倾角测井资料，令各向异性系数初值为1.0，代入方程组迭代求解储层各向异性系数。
2.如权利要求1所述的一种储层电阻率各向异性系数的确定方法，其特征在于，所述步骤101，选取研究区具有自然伽马、电阻率、地层倾角测井资料的井，通过测井资料获得纯泥岩段电阻率值与测量井段的地层倾角，确定求解方程组中的定量参数。
3.如权利要求1所述的一种储层电阻率各向异性系数的确定方法，其特征在于，所述的筛析粒度实验按照《岩心分析方法SY/T 5434-2009》标准规定的流程进行。
4.如权利要求1所述的一种储层电阻率各向异性系数的确定方法，其特征在于，所述步骤104，将纯泥岩电阻率、地层倾角、测井电阻率代入式(1)、式(3)、式(4)、式(5)、式(6)组成的方程组，各向异性系数初值为1.0，在迭代过程中逐渐增加各向异性系数值，通过迭代算法求解方程组的解，确定储层电阻率各向异性系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710848308.5A CN107784159B (zh) | 2017-09-19 | 2017-09-19 | 一种储层电阻率各向异性系数的确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710848308.5A CN107784159B (zh) | 2017-09-19 | 2017-09-19 | 一种储层电阻率各向异性系数的确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107784159A CN107784159A (zh) | 2018-03-09 |
| CN107784159B true CN107784159B (zh) | 2021-09-28 |
Family
ID=61437701
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710848308.5A Active CN107784159B (zh) | 2017-09-19 | 2017-09-19 | 一种储层电阻率各向异性系数的确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107784159B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108643895B (zh) * | 2018-04-03 | 2022-03-01 | 中国石油天然气股份有限公司 | 一种渐变地层的电阻率测井响应计算方法及装置 |
| CN109236284B (zh) * | 2018-09-06 | 2022-05-10 | 中国石油天然气集团有限公司 | 一种感应与侧向联测确定高陡地层真电阻率的方法 |
| CN112696195B (zh) * | 2019-10-23 | 2023-11-28 | 中国石油天然气股份有限公司 | 地层电阻率方位各向异性确定方法及装置 |
| CN111794738B (zh) * | 2020-07-07 | 2023-07-25 | 中国石油天然气集团有限公司 | 一种基于阵列感应与各向异性联合反演地层倾角的方法 |
| CN112115592B (zh) * | 2020-08-31 | 2024-11-05 | 中国石油大学(华东) | 一种基于岩石物理实验的大斜度/水平井电阻率校正方法 |
| CN115898368A (zh) * | 2022-11-10 | 2023-04-04 | 青岛理工大学 | 一种水平井页岩油储层各向异性电阻率提取方法 |
| CN116910422B (zh) * | 2023-09-12 | 2023-12-26 | 广州海洋地质调查局三亚南海地质研究所 | 一种天然气水合物储层泥质含量计算方法、系统及装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101957457A (zh) * | 2009-07-15 | 2011-01-26 | Pgs地球物理公司 | 用于使用电磁测量来确定电阻率各向异性的方法 |
| CN105089663A (zh) * | 2015-07-31 | 2015-11-25 | 中国石油天然气集团公司 | 一种大斜度/水平井地层电阻率各向异性校正方法 |
| CN105507891A (zh) * | 2015-11-30 | 2016-04-20 | 中国石油天然气股份有限公司 | 获取大斜度井地层电阻率各向异性系数的方法及装置 |
| CN106324689A (zh) * | 2016-06-24 | 2017-01-11 | 杭州迅美科技有限公司 | 一种水平井地层环境下电阻率各向异性识别方法 |
| CN107066749A (zh) * | 2017-04-25 | 2017-08-18 | 西安石油大学 | 一种定量评价煤层顶底板封盖性能的方法 |
-
2017
- 2017-09-19 CN CN201710848308.5A patent/CN107784159B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101957457A (zh) * | 2009-07-15 | 2011-01-26 | Pgs地球物理公司 | 用于使用电磁测量来确定电阻率各向异性的方法 |
| CN105089663A (zh) * | 2015-07-31 | 2015-11-25 | 中国石油天然气集团公司 | 一种大斜度/水平井地层电阻率各向异性校正方法 |
| CN105507891A (zh) * | 2015-11-30 | 2016-04-20 | 中国石油天然气股份有限公司 | 获取大斜度井地层电阻率各向异性系数的方法及装置 |
| CN106324689A (zh) * | 2016-06-24 | 2017-01-11 | 杭州迅美科技有限公司 | 一种水平井地层环境下电阻率各向异性识别方法 |
| CN107066749A (zh) * | 2017-04-25 | 2017-08-18 | 西安石油大学 | 一种定量评价煤层顶底板封盖性能的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107784159A (zh) | 2018-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107784159B (zh) | 一种储层电阻率各向异性系数的确定方法 | |
| US11048012B2 (en) | Formation characterization system | |
| Lai et al. | Geophysical well-log evaluation in the era of unconventional hydrocarbon resources: A review on current status and prospects | |
| CN106951660B (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
| Buller et al. | Petrophysical evaluation for enhancing hydraulic stimulation in horizontal shale gas wells | |
| US10254432B2 (en) | Multi-electrode electric field downhole logging tool | |
| US10400590B1 (en) | Method and system for determining a distribution of rock types in geological cells around a wellbore | |
| CN106154322B (zh) | 测井曲线校正方法和装置 | |
| NO20101136L (no) | Karakterisering av bruddlengder og formasjonsresistivitet ut ifra matrise induksjonsdata | |
| CN103790579B (zh) | 随钻地质导向中确定钻头与地层界面距离的方法及装置 | |
| CN110927794A (zh) | 一种致密储层裂缝识别与孔隙度定量计算的方法 | |
| Cluff et al. | Petrophysics of the Lance sandstone reservoirs in Jonah field, Sublette County, Wyoming | |
| US20210332692A1 (en) | Dip correction for non-circular borehole and off-center logging | |
| Lai et al. | Typical misinterpretations and scientific concepts in well-logging geologic studies | |
| CN109657258B (zh) | 地层岩石力学特征剖面的处理方法和装置 | |
| Vij et al. | LWD as the absolute formation evaluation technology: present-day capabilities, limitations, and future developments of LWD technology | |
| US11391860B2 (en) | Anisotropic formation modeling | |
| Christensen et al. | Formation evaluation challenges in Tamar field, offshore Israel | |
| US10295696B2 (en) | Multi-component induction logging data processing in non-circular boreholes | |
| US9720125B2 (en) | Subterranean formation oil mobility quicklook | |
| CN109236284B (zh) | 一种感应与侧向联测确定高陡地层真电阻率的方法 | |
| CA2967266C (en) | Attirbute-indexed multi-instrument logging of drill cuttings | |
| CN105484732B (zh) | 用于水平井钻井地质导向施工过程井深的处理方法 | |
| CN119414467A (zh) | 基于孔内测井联合波形相控技术的岩性精细表征方法 | |
| Wu et al. | Reservoir Characterization and Fluid Typing of Tight Carbonate Reservoir Assisted by High-Resolution Electrical Image Logs: A Case Study of Central Sichuan Basin, China |
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 |