CN110017136B - 一种基于视水层电阻率的水淹层识别与产水率预测方法 - Google Patents
一种基于视水层电阻率的水淹层识别与产水率预测方法 Download PDFInfo
- Publication number
- CN110017136B CN110017136B CN201910194063.8A CN201910194063A CN110017136B CN 110017136 B CN110017136 B CN 110017136B CN 201910194063 A CN201910194063 A CN 201910194063A CN 110017136 B CN110017136 B CN 110017136B
- Authority
- CN
- China
- Prior art keywords
- water
- layer
- resistivity
- apparent
- curve
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 153
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000012216 screening Methods 0.000 claims abstract description 6
- 239000008398 formation water Substances 0.000 claims abstract description 4
- 239000000523 sample Substances 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 230000000704 physical effect Effects 0.000 claims description 2
- 238000005094 computer simulation Methods 0.000 claims 1
- 238000011158 quantitative evaluation Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 8
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 206010051015 Radiation hepatitis Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000035699 permeability 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)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
一种基于视水层电阻率的水淹层识别与产水率预测方法,包括以下步骤:通过筛选测井曲线、油气测试及生产数据,建立视水层电阻率模型;根据步骤1的视水层电阻率模型,通过常规测井曲线计算目标井的地层水电阻率曲线;水淹层定性识别:在测井曲线图的对数刻度道上,将目标井的视水层电阻率曲线与深探测电阻率曲线进行重叠,分析两条曲线的幅度差异,实现水淹层的定性识别;构建产水率模型,针对水淹层的定性识别结果计算产水率曲线,实现产水率预测。采用视水层电阻率曲线与深探测电阻率曲线进行重叠,分析两条曲线的幅度差异,实现水淹层的定性识别;并且计算得到产水率曲线,实现水淹层的定量评价。
Description
技术领域
本发明属于油藏开发技术领域,特别涉及一种基于视水层电阻率的水淹层识别与产水率预测方法。
背景技术
水淹层识别与产水率预测方法经常运用于一些开采已经进入高含水阶段的油田,其主要依据是通过观察常规测井系列测得的测井曲线响应特征,对储层进行定性解释。对于不同的油田,其地质概况不同,所采用的开发方案也不同,使得各油田的测井曲线响应非常复杂,没有统一的规律。因此,目前难以找到一种通用的方法来进行水淹层的识别与产水率的预测。
目前常用的水淹层识别方法有自然电位基线偏移法、自然电位与电阻率曲线对应分析法、冲洗带电阻率法、测井曲线交会图法等。然而这些方法在实际应用时,往往有很多局限性。例如,对于自然电位基线偏移法,当油层大段被水淹时,难以确定原始自然电位基线,其使用效果较差。目前常用的产水率预测方法主要依赖于含水饱和度、束缚水饱和度的准确计算及相渗曲线的获取,然而在有些油田,这些数据并不容易准确计算或获得。因此研发一套易获得(基于常规测井曲线)且适应性强的水淹层识别与产水率预测方法具有重要的意义及广阔的前景。
发明内容
本发明的目的在于提供一种基于视水层电阻率的水淹层识别与产水率预测方法,以解决上述问题。
为实现上述目的,本发明采用以下技术方案:
一种基于视水层电阻率的水淹层识别与产水率预测方法,包括以下步骤:
步骤1,通过筛选测井曲线、油气测试及生产数据,建立视水层电阻率模型;
步骤2,根据步骤1的视水层电阻率模型,通过常规测井曲线计算目标井的地层水电阻率曲线;
步骤3,水淹层定性识别:在测井曲线图的对数刻度道上,将目标井的视水层电阻率曲线与深探测电阻率曲线进行重叠,分析两条曲线的幅度差异,实现水淹层的定性识别;
步骤4,构建产水率模型,针对水淹层的定性识别结果计算产水率曲线,实现产水率预测。
进一步的,步骤1中,筛选数据包括:选取能够反映地层泥值含量、物性特征的测井曲线,选取已被油气测试及生产数据证实的纯水层,读取深探测电阻率曲线的数值,优选并读取常规测井曲线的数值,得到与深探测电阻率相关性相对较好的测井曲线。
进一步的,步骤1中,选择一批水层,读取其常规曲线的值采用多元回归的方法,拟合得到视水层电阻率曲线的计算模型:
log10(Ro)=a1×C1+a2×C2+…+an×CN+b (1)
式中,Ro为视水层电阻率,C1~Cn为测井曲线值,a1~an、b为拟合系数。
进一步的,步骤3中,具体包括:
1)若深探测电阻率值Rt>视水层电阻率值Ro,且幅度差异大(Rt>4Ro),该层为油层或气层,未被水淹;
2)若深探测电阻率值Rt>视水层电阻率值Ro,且幅度差异小(Ro<Rt<4Ro),该层为弱~中水淹层;
3)若深探测电阻率值Rt<=视水层电阻率值Ro,该层为强水淹层。
进一步的,步骤4中,具体包括:
(1)产水率样本数据准备
根据油气水生产数据,通过式(2)计算样本储层的含水率,为产水率模型构建提供基础数据;
Fw=qw/(qo+qw) (2)
式中,Fw为产水率,qo为日产油量,qw为日产水量;
(2)产水率模型构建
首先根据式(3)及式(4)对深探测电阻率与视水层电阻率进行对数归一化,并利用式(5)计算两者之间的差值K,再通过拟合方法得到样本储层的产水率Fw与K的关系式6;
ΔlogRt=(log10(Rt)-log10(Rmin))/(log10(Rmax)-log10(Rmin)) (3)
ΔlogR0=(log10(R0)-log10(Rmin))/(log10(Rmax)-log10(Rmin)) (4)
k=ΔlogRt-ΔlogR0 (5)
Fw=m×kn (6)
式中,ΔlogRt与ΔlogR0分别为深探测电阻率Rt与视水层电阻率R0的对数归一化值,
K为深探测电阻率Rt和视水层电阻率R0的对数归一化差值,m、n为拟合系数;
(3)产水率预测
对于水淹层,利用公式(3)、(4)、(5)、(6)计算产水率曲线,实现产水率的定量预测。
进一步的,当出现Rt<R0,会使得k<=0,此时令k=0.01;当出现计算得到的Fw>1,此时令Fw=1。
与现有技术相比,本发明有以下技术效果:
本发明的实现首先要计算视水层电阻率RO,其计算方法所需要的原始数据非常容易获得,计算公式比较容易实现。在计算得到视水层电阻率RO以后,采用视水层电阻率曲线与深探测电阻率曲线进行重叠,分析两条曲线的幅度差异,该操作简单、实用,实现了水淹层的定性识别。并且还可以进一步计算得到产水率曲线,能够利用常见的测井曲线及数据实现水淹层的定量评价,方法适用性强。本发明的评价结果能够准确识别水淹层,优化投产层位,最终助力产能建设。
附图说明
图1为本发明流程图;
图2产水率预测强水淹层示意图;
图3a和图3b为深探测电阻率对数值与GR、DEN曲线的相关性分析图;
图4为视水层电阻率曲线的计算;
图5深探测电阻率与视水层电阻率重叠识别油气层;
图6深探测电阻率与视水层电阻率重叠识别弱~中水淹层
图7深探测电阻率与视水层电阻率重叠识别强水淹层
图8线性拟合方法构建产水率模型
图9产水率预测(油气层)
图10产水率预测(弱~中水淹层)
图11产水率预测(强水淹层)。
具体实施方式
以下结合附图对本发明进一步说明:
以XX油田为例,说明本发明的具体实施方式。
1、视水层电阻率模型建立
(1)测井曲线、油气测试及生产数据等资料的筛选
首先,选取能够反映地层泥值含量、物性特征的测井曲线,例如GR、U、TH、K、PE、SP、AC、DEN、CNL等曲线。然后,选取已被油气测试及生产数据证实的纯水层,读取深探测电阻率曲线(例如RILD或RLLD)的数值,优选并读取常规测井曲线的数值。例如,XX区块通过深探测电阻率对数值与常规测井曲线的相关性分析,优选得到与深探测电阻率相关性相对较好的GR、DEN测井曲线(图3)。
(2)视水层电阻率模型建立
选择一批水层,读取其常规曲线(例如GR、SP、AC、CNL、DEN、RT等)的值,然后,通过多元回归的方法,建立视水层电阻率计算模型(式7):
log10(R0)=0.004358021×GR+1.597117845×DEN-4.016456738 (7)
2、视地层水电阻率计算
利用式7,通过GR和DEN曲线计算目标井的视水层电阻率曲线(R0)(图4)。
3、水淹层的定性识别
利用视水层电阻率曲线(R0)与深探测电阻率曲线(Rt)重叠后的幅度差异,可以有效的实现水淹层的定性识别。图5、6、7分别为深探测电阻率与视水层电阻率重叠法识别出的油气层、弱~中水淹层及强水淹层。
4、产水率样本数据准备
选择有试油或生产数据的井,统计其射孔层段深探测电阻率Rt和视水层电阻率R0,并按照式(3)计算得到视水层电阻率与深探测电阻率对数增大系数K(表1)。
表1产水率样本数据示例
5、产水率模型构建
利用表1数据及式(6),通过线性拟合方法(图8)构建产水率模型,得到式(8):
Fw=0.0432k-1.072 (8)
6、产水率预测
利用式(8),计算得到产水率曲线,实现目的层段的产水率预测(图9~图11)。图9中1799.2-1809米为油气层,预测产水率为0.1左右;图10中2150-2155.3米为弱~中水淹层。预测产水率为0.3~0.7之间;图11中1730-1737米为强水淹层,产水率为1。
Claims (4)
1.一种基于视水层电阻率的水淹层识别与产水率预测方法,其特征在于,包括以下步骤:
步骤1,通过筛选测井曲线、油气测试及生产数据,建立视水层电阻率模型;
步骤2,根据步骤1的视水层电阻率模型,通过常规测井曲线计算目标井的地层水电阻率曲线;
步骤3,水淹层定性识别:在测井曲线图的对数刻度道上,将目标井的视水层电阻率曲线与深探测电阻率曲线进行重叠,分析两条曲线的幅度差异,实现水淹层的定性识别;
步骤4,构建产水率模型,针对水淹层的定性识别结果计算产水率曲线,实现产水率预测;
步骤3中,具体包括:
1)若深探测电阻率值Rt>视水层电阻率值Ro,且幅度差异大:Rt>4Ro,该层为油层或气层,未被水淹;
2)若深探测电阻率值Rt>视水层电阻率值Ro,且幅度差异小:Ro<Rt<4Ro,该层为弱~中水淹层;
3)若深探测电阻率值Rt<=视水层电阻率值Ro,该层为强水淹层;
步骤4中,具体包括:
(1)产水率样本数据准备
根据油气水生产数据,通过式(2)计算样本储层的产水率,为产水率模型构建提供基础数据;
Fw=qw/(qo+qw) (2)
式中,Fw为产水率,qo为日产油量,qw为日产水量;
(2)产水率模型构建
首先根据式(3)及式(4)对深探测电阻率与视水层电阻率进行对数归一化,并利用式(5)计算两者之间的差值K,再通过拟合方法得到样本储层的产水率Fw与K的关系式6;
ΔlogRt=(log10(Rt)-log10(Rmin))/(log10(Rmax)-log10(Rmin)) (3)
ΔlogR0=(log10(R0)-log10(Rmin))/(log10(Rmax)-log10(Rmin)) (4)
k=ΔlogRt-ΔlogR0 (5)
Fw=m×kn (6)
式中,ΔlogRt与ΔlogR0分别为深探测电阻率Rt与视水层电阻率R0的对数归一化值,K为深探测电阻率Rt和视水层电阻率R0的对数归一化差值,m、n为拟合系数;
(3)产水率预测
对于水淹层,利用公式(3)、(4)、(5)、(6)计算产水率曲线,实现产水率的定量预测。
2.根据权利要求1所述的一种基于视水层电阻率的水淹层识别与产水率预测方法,其特征在于,步骤1中,筛选数据包括:选取能够反映地层泥值含量、物性特征的测井曲线,选取已被油气测试及生产数据证实的纯水层,读取深探测电阻率曲线的数值,选出并读取常规测井曲线的数值,得到与深探测电阻率相关性较好的测井曲线。
3.根据权利要求2所述的一种基于视水层电阻率的水淹层识别与产水率预测方法,其特征在于,步骤1中,选择一批水层,读取其常规曲线的值采用多元回归的方法,拟合得到视水层电阻率曲线的计算模型:
log10(Ro)=a1×C1+a2×C2+…+an×Cn+b (1)
式中,Ro为视水层电阻率,C1~Cn为测井曲线值,a1~an、b为拟合系数。
4.根据权利要求1所述的一种基于视水层电阻率的水淹层识别与产水率预测方法,其特征在于,当出现Rt<R0,会使得k<=0,此时令k=0.01;当出现计算得到的Fw>1,此时令Fw=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910194063.8A CN110017136B (zh) | 2019-03-14 | 2019-03-14 | 一种基于视水层电阻率的水淹层识别与产水率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910194063.8A CN110017136B (zh) | 2019-03-14 | 2019-03-14 | 一种基于视水层电阻率的水淹层识别与产水率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110017136A CN110017136A (zh) | 2019-07-16 |
CN110017136B true CN110017136B (zh) | 2023-01-10 |
Family
ID=67189507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910194063.8A Active CN110017136B (zh) | 2019-03-14 | 2019-03-14 | 一种基于视水层电阻率的水淹层识别与产水率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110017136B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111025409B (zh) * | 2019-12-23 | 2021-03-02 | 中国石油大学(北京) | 一种水淹层评价方法、装置及存储介质 |
CN117266843B (zh) * | 2023-09-27 | 2024-03-26 | 广东海洋大学 | 油藏水淹层识别方法、系统、装置及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1040245A (zh) * | 1988-07-25 | 1990-03-07 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
CN103670390A (zh) * | 2013-12-20 | 2014-03-26 | 中国石油天然气集团公司 | 一种水淹层测井评价方法及系统 |
CN104806231A (zh) * | 2014-01-27 | 2015-07-29 | 中国石油化工股份有限公司 | 稠油蒸汽吞吐水淹层的定量评价方法 |
CN104847340A (zh) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | 一种水淹层测井定量评价方法 |
CN105447762A (zh) * | 2015-12-08 | 2016-03-30 | 中国石油天然气集团公司 | 一种流体替换的低渗透油藏水淹信息的计算方法 |
CN106503295A (zh) * | 2016-09-22 | 2017-03-15 | 中国石油天然气股份有限公司 | 一种利用状态空间模型解释油田水淹层的方法及装置 |
CN106869913A (zh) * | 2017-03-03 | 2017-06-20 | 中国地质大学(北京) | 一种利用井地电位技术探测油田注水井水驱前缘的方法 |
CN107725044A (zh) * | 2017-09-27 | 2018-02-23 | 中国石油天然气集团公司 | 基于阵列感应、侧向测井的砂岩含气储层产水率预测的方法 |
CN107832574A (zh) * | 2017-09-19 | 2018-03-23 | 中国石油天然气集团公司 | 一种基于随钻测井的水平井水淹层解释方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101749012B (zh) * | 2008-12-08 | 2012-12-12 | 中国石油天然气集团公司 | 一种油层开采程度的确定方法 |
US10386531B2 (en) * | 2013-03-08 | 2019-08-20 | Schlumberger Technology Corporation | Geological model analysis incorporating cross-well electromagnetic measurements |
US20160154133A1 (en) * | 2013-05-07 | 2016-06-02 | Halliburton Energy Services, Inc. | Systems and methods of providing compensated geological measurements |
CN104863574B (zh) * | 2014-02-20 | 2018-03-13 | 中国石油化工股份有限公司 | 一种适用于致密砂岩储层的流体识别方法 |
WO2017086956A1 (en) * | 2015-11-18 | 2017-05-26 | Halliburton Energy Services, Inc. | Monitoring water floods using potentials between casing-mounted electrodes |
-
2019
- 2019-03-14 CN CN201910194063.8A patent/CN110017136B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1040245A (zh) * | 1988-07-25 | 1990-03-07 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
CN103670390A (zh) * | 2013-12-20 | 2014-03-26 | 中国石油天然气集团公司 | 一种水淹层测井评价方法及系统 |
CN104806231A (zh) * | 2014-01-27 | 2015-07-29 | 中国石油化工股份有限公司 | 稠油蒸汽吞吐水淹层的定量评价方法 |
CN104847340A (zh) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | 一种水淹层测井定量评价方法 |
CN105447762A (zh) * | 2015-12-08 | 2016-03-30 | 中国石油天然气集团公司 | 一种流体替换的低渗透油藏水淹信息的计算方法 |
CN106503295A (zh) * | 2016-09-22 | 2017-03-15 | 中国石油天然气股份有限公司 | 一种利用状态空间模型解释油田水淹层的方法及装置 |
CN106869913A (zh) * | 2017-03-03 | 2017-06-20 | 中国地质大学(北京) | 一种利用井地电位技术探测油田注水井水驱前缘的方法 |
CN107832574A (zh) * | 2017-09-19 | 2018-03-23 | 中国石油天然气集团公司 | 一种基于随钻测井的水平井水淹层解释方法 |
CN107725044A (zh) * | 2017-09-27 | 2018-02-23 | 中国石油天然气集团公司 | 基于阵列感应、侧向测井的砂岩含气储层产水率预测的方法 |
Non-Patent Citations (3)
Title |
---|
LC砂砾岩油藏出水类型及控水对策研究;姚军朋等;《长江大学学报(自科版)》;20160215(第05期);62-66 * |
Osteopontin promotes cancer cell drug resistance, invasion, and lactate production and is associated with poor outcome of patients with advanced non-small-cell lung cancer;Xiaoping Ouyang;《OncoTargets and Therapy》;20181231;5933-5941 * |
新站油田水淹层识别技术研究;刘江南;《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》;20150315(第3期);5-41 * |
Also Published As
Publication number | Publication date |
---|---|
CN110017136A (zh) | 2019-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106468172B (zh) | 一种超低渗砂岩油藏低阻储层测井解释方法 | |
CN109838230B (zh) | 油藏水淹层的定量评价方法 | |
CN101892837B (zh) | 地层因数确定方法及含油饱和度确定方法 | |
RU2315339C2 (ru) | Система петрофизической оценки в реальном времени | |
CN110847901B (zh) | 一种变矿化度地层水下的致密砂岩储层流体识别方法 | |
CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
CN108049866B (zh) | 二维核磁共振测井致密气藏定量评价方法 | |
CN108694264B (zh) | 一种确定页岩气储层渗透率的方法 | |
CN106154343B (zh) | 计算致密油储层的含油饱和度的方法 | |
CN110043254B (zh) | 一种基于电缆地层测试资料地层有效渗透率的获取方法 | |
CN107590550A (zh) | 有关超低渗透储层油田产能评价及预测的方法 | |
CN104863574A (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
Qin et al. | Fast prediction method of Archie’s cementation exponent | |
CN112363242B (zh) | 基于测录井融合的储层流体识别方法与装置 | |
CN112835124B (zh) | 基于成像测井和阵列声波测井资料的裂缝有效性评价方法 | |
CN110017136B (zh) | 一种基于视水层电阻率的水淹层识别与产水率预测方法 | |
CN110905493B (zh) | 一种测量井下地层流体污染率的方法 | |
Feng et al. | Accurate determination of water saturation in tight sandstone gas reservoirs based on optimized Gaussian process regression | |
CN106777649B (zh) | 一种裂缝型储集层孔隙结构定量评价方法 | |
CN113216945B (zh) | 一种致密砂岩储层渗透性定量评价方法 | |
CN104632206B (zh) | 一种石油测井标准层确定方法及装置 | |
CN111625750A (zh) | 一种针对致密砂岩储层气层与气水同层的识别方法 | |
CN114428365A (zh) | 一种基于含气特征的致密砂岩气层测井识别方法及装置 | |
CN108412488B (zh) | 快速确定页岩气储层有机孔隙度的测井方法 | |
CN103197348A (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 |