CN108717498A - 一种油藏相渗曲线模型及相渗曲线计算方法 - Google Patents
一种油藏相渗曲线模型及相渗曲线计算方法 Download PDFInfo
- Publication number
- CN108717498A CN108717498A CN201810505945.7A CN201810505945A CN108717498A CN 108717498 A CN108717498 A CN 108717498A CN 201810505945 A CN201810505945 A CN 201810505945A CN 108717498 A CN108717498 A CN 108717498A
- Authority
- CN
- China
- Prior art keywords
- model
- percolation curve
- models
- oil reservoir
- phase percolation
- 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.)
- Granted
Links
- 238000005325 percolation Methods 0.000 title claims abstract description 45
- 238000000205 computational method Methods 0.000 title claims abstract description 11
- 238000002474 experimental method Methods 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims abstract description 5
- 230000035699 permeability Effects 0.000 claims description 22
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 17
- 238000000034 method Methods 0.000 abstract 1
- 238000011161 development Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种油藏相渗曲线模型及相渗曲线计算方法，包括一种油藏相渗曲线模型和油藏相渗曲线计算方法。模型包括括GML模型和GBL模型，模型获取方法依次包括以下步骤：(1)将实验测试所得相渗曲线和GML模型、GBL模型进行拟合，得到多组β、γ和m的参数值。(2)确定模型中各参数的平均值。将步骤1拟合得到的多组β、γ等参数值求取平均值，得到通用模型参数，将通用模型参数代入模型中，得到通用模型并验证通用模型的拟合度。(3)获取油藏相渗曲线模型。将步骤2得到的β、γ等参数的通用模型参数值代入GML模型和GBL模型中，计算得到油藏相渗曲线模型。(4)根据步骤(3)所得通用模型计算得到油藏相渗曲线。
Description
技术领域
本发明属于油气田开发室内实验技术领域，具体地涉及一种油藏相渗曲线模型及相渗曲线计算方法。
背景技术
相对渗透率曲线(相渗曲线)是描述油藏中流体渗流时，岩石和流体相互作用的重要参数，是油藏开发预测和评价中最重要的基础资料，其准确与否直接决定了油藏开发效果预测和评价的精度。自Puecell(Purcell W R.Capillary pressures-theirmeasurement using mercury and the calculation of permeability therefrom[J].Journal of Petroleum Technology，1949，1(2)：39-48)通过毛管束模型推导出由毛管压力计算相对渗透率的模型后，出现了很多改进方案或复杂网络模型，如考虑了毛管对润湿相束缚作用的Burdine模型(Burdine N T.Relative permeability calculations frompore size distribution data[J].Journal of Petroleum Technology，1953，5(3)：71-78)，和考虑孔隙非均质性的Mualem模型(Mualem Y.A new model for predicting thehydraulic conductivity of unsaturated porous media[J].Water ResourcesResearch，1976，12(3)：513-522)。但是，这些模型没有充分考虑动态毛管效应等因素对相对渗透率的影响，实验数据拟合度较低，因此需要准确度更高、拟合度更高的相渗曲线模型，并计算得到准确的油藏相渗曲线。
发明内容
本发明的目的在于提供一种油藏相渗曲线模型及相渗曲线计算方法，得到新的油藏相渗曲线模型，克服现有油藏相渗曲线模型拟合度不高的不足，并计算得到准确的油藏相渗曲线。该油藏相渗曲线模型充分考虑动态毛管效应对相对渗透率的影响，能够较好地拟合油藏注水过程中的相对渗透率实验数据，油藏相渗曲线的计算方法较为简便和准确，可为评价油藏注水开发效果提供准确的数据支持。本发明的目的通过以下技术方案来实现：
1.一种油藏相渗曲线模型，其特征在于：模型包括GML模型和GBL模型，所述GML模型为：
式中，Krw为水相相对渗透率；
Kro为油相相对渗透率；
为有效含水饱和度，无量纲；
Sw为含水饱和度，无量纲；
Swi为初始含水饱和度，无量纲；
β、γ和m为参数，可以通过实验拟合得到，其中0<m<1。
所述GBL模型为：
式中，0<m<1。
2.一种油藏相渗曲线计算方法，依次包括以下步骤：
(1)将实验测试所得相渗曲线和GML模型、GBL模型进行拟合，得到多组β、γ和m的参数值。
(2)确定模型中各参数的平均值。将步骤1拟合得到的多组β、γ等参数值求取平均值，得到通用模型参数，并验证通用模型参数的拟合度。
(3)得到油藏相渗曲线模型。将步骤2得到的β、γ等参数的通用模型参数值代入GML模型和GBL模型中，计算得到油藏相渗曲线模型。
(4)根据步骤(3)所得模型计算得到油藏相渗曲线
本发明具有以下优点：
在计算相对渗透率时，充分考虑了毛管压力、动态毛管效应、孔隙大小、润湿性和孔隙非均质性等因素对相对渗透率的影响，因此实验得到的油藏相对渗透率曲线与模型拟合效果好，计算得到的油藏相渗曲线可为评价油藏开发效果提供准确的数据支持。
附图说明
图1为本一种油藏相渗曲线计算方法的流程图。
图2和图3为实施例中油藏相渗曲线模型代表性拟合结果。
图4为实施例中相渗曲线的代表性计算结果。
具体实施方式
为使本发明实施例的目的，技术方案和优点更加清楚，下面将结合本发明实施例中的附图和实验数据处理方法应用实例进一步说明本发明。显然，所描述的实施例是本发明一部分施例，不是全部实施例。对于本领域普通技术人员来说，在不付出创造性劳动的前提下，可以根据这些附图获得其他的附图。
具体如下：
如图1所示，为该油藏相渗曲线计算方法的具体流程图，步骤如下：
1.将实验测试所得相渗曲线和GML模型、GBL模型进行拟合，所述GML模型为：
式中，Krw为水相相对渗透率；
Kro为油相相对渗透率；
为有效含水饱和度，无量纲；
Sw为含水饱和度，无量纲；
Swi为初始含水饱和度，无量纲；
β、γ和m为参数，通过实验拟合得到，其中0<m<1。
所述GBL模型为：
式中，0<m<1。
代表性拟合结果如图2所示，拟合得到的β、γ和m如表所示：
拟合结果表明，实验水相相对渗透率数据与GML模型拟合较好，油相相对渗透率数据与GBL模型拟合较好。因此计算油藏相对渗透率曲线时，采用GML模型计算水相相对渗透率，GBL模型计算油相相对渗透率。
2.将步骤1拟合得到的多组β、γ等参数求取平均值，得到通用模型参数，并验证通用模型参数的拟合度，通用模型参数值和拟合度如下表所示：
3.获取油藏相渗曲线模型。将步骤2得到的β、γ等参数的通用模型参数值代入GML模型和GBL模型中，得到该油藏相渗曲线模型，计算油藏水相相对渗透率的GML模型为：
计算该油藏油相相对渗透率的GBL模型为：
4.根据步骤3所得模型计算得到油藏相渗曲线，代表性结果如图4所示。
Claims (2)
1.一种油藏相渗曲线模型，其特征在于：模型包括GML模型和GBL模型，所述GML模型为：
式中，Krw为水相相对渗透率；
Kro为油相相对渗透率；
为有效含水饱和度，无量纲；
Sw为含水饱和度，无量纲；
Swi为初始含水饱和度，无量纲；
β、γ和m为参数，可以通过实验拟合得到，其中0<m<1；
所述GBL模型为：
式中，0<m<1。
2.一种油藏相渗曲线计算方法，依次包括以下步骤：
(1)将实验测试所得相渗曲线和GML模型、GBL模型进行拟合，得到多组β、γ和m的参数值；
(2)确定模型中各参数的平均值，将步骤1拟合得到的多组β、γ等参数值求取平均值，得到通用模型参数，将通用模型参数代入模型中得到通用模型，并验证通用模型参数的拟合度；
(3)获取油藏相渗曲线模型，将步骤2得到的β、γ等参数的通用模型参数值代入GML模型和GBL模型中，得到油藏相渗曲线模型；
(4)根据步骤(3)所得模型计算得到油藏相渗曲线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810505945.7A CN108717498B (zh) | 2018-05-19 | 2018-05-19 | 一种油藏相渗曲线模型及相渗曲线计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810505945.7A CN108717498B (zh) | 2018-05-19 | 2018-05-19 | 一种油藏相渗曲线模型及相渗曲线计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108717498A true CN108717498A (zh) | 2018-10-30 |
| CN108717498B CN108717498B (zh) | 2022-02-15 |
Family
ID=63900478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810505945.7A Active CN108717498B (zh) | 2018-05-19 | 2018-05-19 | 一种油藏相渗曲线模型及相渗曲线计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108717498B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020258515A1 (zh) * | 2019-06-28 | 2020-12-30 | 浙江大学 | 多柱连续流层析设计及分析的方法 |
| CN113806998A (zh) * | 2020-06-16 | 2021-12-17 | 中国石油化工股份有限公司 | 一种储层相渗曲线仿真方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104912550A (zh) * | 2015-05-05 | 2015-09-16 | 中国海洋石油总公司 | 一种核磁共振测井资料定量计算储层产液剖面方法 |
| CN106204304A (zh) * | 2016-07-12 | 2016-12-07 | 中国石油新疆油田分公司勘探开发研究院 | 一种砾岩油藏聚驱相对渗透率曲线的确定方法 |
| US20170169139A1 (en) * | 2015-12-10 | 2017-06-15 | Saudi Arabian Oil Company | Predicting and modeling changes in capillary pressure and relative permeabilities in a porous medium due to mineral precipitation and dissolution |
| CN107340219A (zh) * | 2017-07-07 | 2017-11-10 | 西南石油大学 | 一种油藏动态毛管效应评价实验数据处理方法 |
-
2018
- 2018-05-19 CN CN201810505945.7A patent/CN108717498B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104912550A (zh) * | 2015-05-05 | 2015-09-16 | 中国海洋石油总公司 | 一种核磁共振测井资料定量计算储层产液剖面方法 |
| US20170169139A1 (en) * | 2015-12-10 | 2017-06-15 | Saudi Arabian Oil Company | Predicting and modeling changes in capillary pressure and relative permeabilities in a porous medium due to mineral precipitation and dissolution |
| CN106204304A (zh) * | 2016-07-12 | 2016-12-07 | 中国石油新疆油田分公司勘探开发研究院 | 一种砾岩油藏聚驱相对渗透率曲线的确定方法 |
| CN107340219A (zh) * | 2017-07-07 | 2017-11-10 | 西南石油大学 | 一种油藏动态毛管效应评价实验数据处理方法 |
Non-Patent Citations (3)
| Title |
|---|
| YING LI ET AL: "Capillarity characters measurement and effects analysis in different permeability formations during waterflooding", 《FUEL》 * |
| YONGZHONG LIU ET AL.: "Effects of capillary pressure–fluid saturation–relative permeability relationships on predicting carbon dioxide migration during injection into saline aquifers", 《ENERGY PROCEDIA》 * |
| 闫国亮: "基于数字岩心储层渗透率模型研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020258515A1 (zh) * | 2019-06-28 | 2020-12-30 | 浙江大学 | 多柱连续流层析设计及分析的方法 |
| CN113806998A (zh) * | 2020-06-16 | 2021-12-17 | 中国石油化工股份有限公司 | 一种储层相渗曲线仿真方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108717498B (zh) | 2022-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106599449B (zh) | 一种溶洞体积计算的试井解释方法 | |
| CN107525746B (zh) | 一种表征致密油藏低速非达西渗流特征的方法及装置 | |
| CN105569641B (zh) | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 | |
| CN107340219B (zh) | 一种油藏动态毛管效应评价实验数据处理方法 | |
| CN108487904B (zh) | 一种基于图版的消除末端效应的相渗曲线校正方法 | |
| CN110346258A (zh) | 一种致密岩石油相相对渗透率测定的方法 | |
| CN107462936B (zh) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 | |
| CN105781539A (zh) | 一种致密油气储层饱和度测井计算方法 | |
| RU2020141650A (ru) | Оценка сетей потоков | |
| CN107063967A (zh) | 一种基于微观孔隙通道的饱和∕非饱和渗透系数预测方法 | |
| CN110598167A (zh) | 低渗透油藏油水相对渗透率实验数据的处理方法 | |
| CN109948219A (zh) | 一种预测不同孔隙比条件下非饱和绝对渗透系数的方法 | |
| CN112986097B (zh) | 一种确定致密储层稳态法相对渗透率曲线的实验测量方法 | |
| CN108717498A (zh) | 一种油藏相渗曲线模型及相渗曲线计算方法 | |
| CN102915406A (zh) | 径向流条件下油水相对渗透率曲线的计算方法 | |
| CN107063935B (zh) | 一种测量co2在油水两相间扩散过程及扩散系数的方法 | |
| CN104569350A (zh) | 密闭取芯不规则全直径岩样原始含水饱和度的测试方法 | |
| CN110161071A (zh) | 一种评价致密储层可动流体最小喉道半径的方法 | |
| CN207073409U (zh) | 一种气藏非均匀平面水侵模拟实验装置 | |
| CN102052074A (zh) | 一种岩心数字化确定岩电参数的方法 | |
| CN106771071B (zh) | 一种基于油水相渗的密闭取心饱和度校正方法 | |
| CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
| CN105844011B (zh) | 一种基于毛管模型的渗透率计算方法 | |
| CN108489878B (zh) | 一种基于数值模拟迭代消除末端效应的相渗曲线校正方法 | |
| CN105938084B (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 |