CN107741605A - 基于时间推移的微元化导电模型求水淹层相关参数的方法 - Google Patents
基于时间推移的微元化导电模型求水淹层相关参数的方法 Download PDFInfo
- Publication number
- CN107741605A CN107741605A CN201710704339.3A CN201710704339A CN107741605A CN 107741605 A CN107741605 A CN 107741605A CN 201710704339 A CN201710704339 A CN 201710704339A CN 107741605 A CN107741605 A CN 107741605A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mfrac
- water
- swb
- 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.)
- Withdrawn
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000002474 experimental method Methods 0.000 claims description 7
- 239000008398 formation water Substances 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000008859 change Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000008346 aqueous phase Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002269 spontaneous effect Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 1
Classifications
-
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Health & Medical Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Dispersion Chemistry (AREA)
- Remote Sensing (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明提供了一种基于时间推移的微元化导电模型求水淹层相关参数的方法,经循环一次性求解出混合水电阻率Rwz、含水饱和度Sw以及产水率Fw等水淹层关键参数。解决了水淹层含水饱和度难以求准的问题,提高了水淹油田开发经济效益。
Description
技术领域
本发明涉及地球物理测井技术领域,尤其涉及基于时间推移的微元化导电模型求水淹层相关参数的方法。
背景技术
目前水淹层的水淹级别主要是根据产水率来划分的,要想求准产水率,准确计算出水淹层含水饱和度是关键。油层水淹后由于注入水矿化度与原生水矿化度不同,导致混合水电阻率不断变化,使得含水饱和度难以求准。目前水淹层混合水电阻率求取方法主要有根据自然电位曲线计算以及通过理论推导或者数值模拟方法建模两种,但是这些方法尚存在如下技术问题:
1、受限于测井系列的影响,大量未测量自然电位曲线的水淹井无法应用自然电位法计算混合水电阻率;
2、绝大部分的理论推导或者数值模拟方法都是以未水淹时为起始点,而未考虑到水淹其实是一个长期的过程,各相关参数在水淹过程中会发生实时的变化。
发明内容
本发明提供了一种基于时间推移的微元化导电模型求水淹层相关参数的方法,首先通过岩心分析相关实验建立束缚水、残余油、相对渗透率及产水率相关模型,然后采用基于时间推移的微元化导电模型迭代联合求解出水淹层各关键参数,提高了水淹层解释的精度。
为解决上述技术问题,本发明提供了一种基于时间推移的微元化导电模型求水淹层相关参数的方法,包括以下步骤:
S1:通过岩电实验数据,计算出区块岩电参数a、b、m、n;
S2:通过毛管压力试验、离心机实验及油水相渗实验建立区块束缚水、残余油和油水相对渗透率计算模型;
S3:通过岩性曲线(GR、SP曲线)计算出泥质含量Vsh,通过物性曲线:DEN、CNL和AC曲线计算出孔隙度φ,通过电性曲线读出泥岩电阻率Rsh;
S4:令含水饱和度由Swi到1-Sor变化,步长ΔSw=0.1%;
S5:含水饱和度取某一个定值时,令未参与混合的原生水饱和度Swb由Swi到0变化,步长ΔSwb=0.1%;
S6:计算出此时的产水率Fw,再根据产水率计算出产出水矿化度为:
根据Pz2计算出混合水电阻率Rw2,然后再根据并联导电原理,将混合水电阻率与未参与混合的原生地层水并联,得到视地层水电阻率:
再由印尼公式:
计算出新的Sw2’,算出含水饱和度误差Δ=|Sw2'-Sw2|。
S7:循环Swb、Sw,采用时间推移微元化思想,以上一时刻的混合液电阻率和含水饱和度作为起点,计算新的混合液电阻率和含水饱和度,找到令Δ最小的Sw以及Swb作为真实值,同时还可以得到混合水电阻率Rwz以及产水率Fw。
本申请实施例中提供的一个或多个技术方案,至少具有如下技术效果或优点:
本发明通过时间推移的思想,建立了能更好的反映水淹层客观规律的混合液矿化度变化规律模型,经循环一次性求解出混合水电阻率Rwz、含水饱和度Sw以及产水率Fw等多个水淹层关键参数。解决了水淹层含水饱和度难以求准的问题,提高了水淹油田开发经济效益。
附图说明
图1方法处理效果图。
具体实施方式
为了更好的理解上述技术方案,下面以具体的实施方式对上述技术方案进行详细的说明。
本实施例所述的一种基于时间推移的微元化导电模型求水淹层相关参数的方法,S1:通过岩电实验数据,计算出区块岩电参数a、b、m、n;
S2:通过毛管压力试验、离心机实验及油水相渗实验建立区块束缚水、残余油和油水相对渗透率计算模型;
S3:通过岩性曲线(GR、SP曲线)计算出泥质含量Vsh,通过物性曲线:DEN、CNL和AC曲线计算出孔隙度φ,通过电性曲线读出泥岩电阻率Rsh;
S4:令含水饱和度由Swi到1-Sor变化,步长ΔSw=0.1%;
S5:含水饱和度取某一个定值时,令未参与混合的原生水饱和度Swb由Swi到0变化,步长ΔSwb=0.1%;
S6:计算出此时的产水率Fw,再根据产水率计算出产出水矿化度为:
根据Pz2计算出混合水电阻率Rw2,然后再根据并联导电原理,将混合水电阻率与未参与混合的原生地层水并联,得到视地层水电阻率:
再由印尼公式:
计算出新的Sw2’,算出含水饱和度误差Δ=|Sw2'-Sw2|。
S7:循环Swb、Sw,采用时间推移微元化思想,以上一时刻的混合液电阻率和含水饱和度作为起点,计算新的混合液电阻率和含水饱和度,找到令Δ最小的Sw以及Swb作为真实值,同时还可以得到混合水电阻率Rwz以及产水率Fw。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (1)
1.基于时间推移的微元化导电模型求水淹层相关参数的方法,其特征在于,包括如下步骤:
S1:通过岩电实验数据,计算出区块岩电参数a、b、m、n;
S2:通过毛管压力试验、离心机实验及油水相渗实验建立区块束缚水、残余油和油水相对渗透率计算模型;
S3:通过岩性曲线(GR、SP曲线)计算出泥质含量Vsh,通过物性曲线:DEN、CNL和AC曲线计算出孔隙度φ,通过电性曲线读出泥岩电阻率Rsh;
S4:令含水饱和度由Swi到1-Sor变化,步长ΔSw=0.1%;
S5:含水饱和度取某一个定值时,令未参与混合的原生水饱和度Swb由Swi到0变化,步长ΔSwb=0.1%;
S6:计算出此时的产水率Fw,再根据产水率计算出产出水矿化度为:
<mrow>
<msub>
<mi>Pz</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>F</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>Sw</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Sw</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Sw</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Swb</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>Pz</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Swb</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Swb</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>P</mi>
<mi>i</mi>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>F</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>Sw</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Sw</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<msub>
<mi>Sw</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Swb</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
根据Pz2计算出混合水电阻率Rw2,然后再根据并联导电原理,将混合水电阻率与未参与混合的原生地层水并联,得到视地层水电阻率:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>w</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Sw</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>Swb</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mi>R</mi>
<mi>w</mi>
<mi>i</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>Sw</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Swb</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>Rw</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mfrac>
</mrow>
再由印尼公式:
<mrow>
<msup>
<msub>
<mi>S</mi>
<mrow>
<mi>w</mi>
<mi>z</mi>
</mrow>
</msub>
<mfrac>
<mi>n</mi>
<mn>2</mn>
</mfrac>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
</msqrt>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>0.5</mn>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
<msqrt>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
</msqrt>
</mfrac>
<mo>+</mo>
<msqrt>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&phi;</mi>
<mi>e</mi>
</msub>
<mi>m</mi>
</msup>
</mrow>
<mrow>
<msub>
<mi>aR</mi>
<mrow>
<mi>w</mi>
<mi>z</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
计算出新的Sw2’,算出含水饱和度误差Δ=|Sw2'-Sw2|。
S7:循环Swb、Sw,采用时间推移微元化思想,以上一时刻的混合液电阻率和含水饱和度作为起点,计算新的混合液电阻率和含水饱和度,找到令Δ最小的Sw以及Swb作为真实值,同时还可以得到混合水电阻率Rwz以及产水率Fw。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710704339.3A CN107741605A (zh) | 2017-08-17 | 2017-08-17 | 基于时间推移的微元化导电模型求水淹层相关参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710704339.3A CN107741605A (zh) | 2017-08-17 | 2017-08-17 | 基于时间推移的微元化导电模型求水淹层相关参数的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107741605A true CN107741605A (zh) | 2018-02-27 |
Family
ID=61235527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710704339.3A Withdrawn CN107741605A (zh) | 2017-08-17 | 2017-08-17 | 基于时间推移的微元化导电模型求水淹层相关参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107741605A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111025409A (zh) * | 2019-12-23 | 2020-04-17 | 中国石油大学(北京) | 一种水淹层评价方法、装置及存储介质 |
CN111594115A (zh) * | 2020-05-22 | 2020-08-28 | 中国石油天然气集团有限公司 | 一种基于水淹层饱和度变化计算累积注水量的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5663499A (en) * | 1995-10-20 | 1997-09-02 | Semmelbeck; Mark E. | Method for estimating permeability from multi-array induction logs |
CN104847340A (zh) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | 一种水淹层测井定量评价方法 |
CN105888631A (zh) * | 2016-05-10 | 2016-08-24 | 中国海洋石油总公司 | 一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法 |
CN107038266A (zh) * | 2016-02-04 | 2017-08-11 | 中石化石油工程技术服务有限公司 | 基于数字岩心的水淹层新饱和度模型建立方法 |
-
2017
- 2017-08-17 CN CN201710704339.3A patent/CN107741605A/zh not_active Withdrawn
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5663499A (en) * | 1995-10-20 | 1997-09-02 | Semmelbeck; Mark E. | Method for estimating permeability from multi-array induction logs |
CN104847340A (zh) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | 一种水淹层测井定量评价方法 |
CN107038266A (zh) * | 2016-02-04 | 2017-08-11 | 中石化石油工程技术服务有限公司 | 基于数字岩心的水淹层新饱和度模型建立方法 |
CN105888631A (zh) * | 2016-05-10 | 2016-08-24 | 中国海洋石油总公司 | 一种裂缝性稠油油藏热水驱物理模拟相似准数的确定方法 |
Non-Patent Citations (3)
Title |
---|
王向公,等: "《地球物理测井数字处理方法》", 31 December 2013 * |
范翔宇: "《钻井储层污染损害精细评价方法》", 31 August 2012 * |
谭伟,等: "注水开发油田水淹层混合液电阻率求取新方法", 《石油钻探技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111025409A (zh) * | 2019-12-23 | 2020-04-17 | 中国石油大学(北京) | 一种水淹层评价方法、装置及存储介质 |
CN111594115A (zh) * | 2020-05-22 | 2020-08-28 | 中国石油天然气集团有限公司 | 一种基于水淹层饱和度变化计算累积注水量的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2583121B1 (en) | Production simulator for simulating a mature hydrocarbon field | |
US7983886B2 (en) | Methods for performing simulation of surfactant flooding of a hydrocarbon reservoir | |
CN105740563B (zh) | 一种成熟油田二次开发之优势通道识别方法 | |
CN104806231B (zh) | 稠油蒸汽吞吐水淹层的定量评价方法 | |
CN108603402A (zh) | 对矿物质沉淀和溶解造成的多孔介质中毛细管压力和相对渗透率的变化进行建模和预测 | |
CN106501144A (zh) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 | |
CN103485770B (zh) | 一种基于人工神经网络求取含油饱和度的方法及系统 | |
CN104246484A (zh) | 用于从 nmr 确定润湿性的方法 | |
CN107622139A (zh) | 裂缝渗透率的计算方法 | |
McPhee et al. | Best practice in coring and core analysis | |
Andersen et al. | Simulation interpretation of capillary pressure and relative permeability from laboratory waterflooding experiments in preferentially oil-wet porous media | |
CN106932324A (zh) | 一种确定高含水砂岩油藏储层渗透率变化规律的方法 | |
CN106503295B (zh) | 一种利用状态空间模型解释油田水淹层的方法及装置 | |
CN107741605A (zh) | 基于时间推移的微元化导电模型求水淹层相关参数的方法 | |
US20230160304A1 (en) | Method and system for predicting relative permeability curve based on machine learning | |
Mahzari et al. | Co-history matching: a way forward for estimating representative saturation functions | |
CN104500051A (zh) | 一种确定水淹层混合地层水宏观俘获截面的方法 | |
Schwalbert et al. | Skin-Factor Equations for Anisotropic Wormhole Networks and Limited-Entry Completions | |
Aghabozorgi et al. | An investigation of polymer adsorption in porous media using pore network modelling | |
CA3041614A1 (en) | Global surface paleo-temperature modeling tool | |
CN111827991B (zh) | 一种计算水淹级别水淹层对应剩余油饱和度的方法及系统 | |
CN111079260B (zh) | 一种非线性渗流数值模拟方法 | |
CN115618750B (zh) | 一种基于耦合神经网络的地下油水渗流代理模型 | |
Triadis et al. | The Green–Ampt limit with reference to infiltration coefficients | |
Yaralidarani et al. | Co-estimation of saturation functions (kr and P c) from unsteady-state core-flood experiment in tight carbonate rocks |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180227 |
|
WW01 | Invention patent application withdrawn after publication |