CN108729912A - Production split method suitable for reservoir numerical simulation - Google Patents
Production split method suitable for reservoir numerical simulation Download PDFInfo
- Publication number
- CN108729912A CN108729912A CN201810356826.XA CN201810356826A CN108729912A CN 108729912 A CN108729912 A CN 108729912A CN 201810356826 A CN201810356826 A CN 201810356826A CN 108729912 A CN108729912 A CN 108729912A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- production
- numerical simulation
- perforation
- payzone
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 title claims abstract description 34
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical compound C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 claims abstract description 28
- 239000012530 fluid Substances 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 12
- 239000007789 gas Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The present invention provides a kind of production split method suitable for reservoir numerical simulation, should include suitable for the production split method of reservoir numerical simulation：Step 1, the output coefficient of payzone is calculated in the active porosity volume of the oil saturation of each perforation unit, oil phase fluid density and each layer according to well；Step 2, full well production coefficient when stratified reservoir is developed simultaneously is calculated；Step 3, a point coefficient is split using what the output coefficient of payzone and full well production coefficient calculated payzone.The production split method suitable for reservoir numerical simulation establishes the output situation for capableing of accurate description oil reservoir development and each payzone reservoir numerical simulation, is of great significance to oil-gas reservoir fine development, lays the foundation to improve Research Numerical Simulation Techique.
Description
Technical field
The present invention relates to oil reservoir development technical field, especially relates to a kind of yield suitable for reservoir numerical simulation and split
Divide method.
Background technology
Conventional KH production split methods are mainly the influence of the permeability and core intersection that consider substratum, by producing well
Production split is to each substratum.But this method cannot identify reservoir property.Phase is introduced on the basis of conventional KH splits point method
A point method is split in influence to permeability, referred to as KHK, and this method improves the calculating of formation capacity, but for storage
The heterogeneous stronger thin interbed sandstone oil reservoir applicability of layer is poor.KNK production split methods are oily on the basis of KHK methods
Well substratum controlled geological reserves are modified the model.The above method is on the basis for considering substratum fluid ability
Upper foundation, do not reflect the degree of saturation of near wellbore zone fluid.For this purpose, we have invented one kind being suitable for reservoir numerical simulation
Production split method, the reserves scale of concentrated expression payzone, the production status of near wellbore zone and fluid degree of saturation.Thus
We have invented a kind of new production split methods suitable for reservoir numerical simulation, solve the above technical problem.
Invention content
The object of the present invention is to provide a kind of reserves scale of concentrated expression payzone, the production status and fluid of near wellbore zone
Degree of saturation, to improve the production split side suitable for reservoir numerical simulation that Research Numerical Simulation Techique provides technical support
Method.
The purpose of the present invention can be achieved by the following technical measures：Production split side suitable for reservoir numerical simulation
Method, the production split method that should be suitable for reservoir numerical simulation include：Step 1, the oil-containing according to well in each perforation unit is saturated
The active porosity volume of degree, oil phase fluid density and each layer, calculates the output coefficient of payzone；Step 2, stratified reservoir is calculated simultaneously
Full well production coefficient when exploitation；Step 3, calculating splitting for payzone using the output coefficient of payzone and full well production coefficient point is
Number.
The purpose of the present invention can be also achieved by the following technical measures：
In step 1, the oil saturation of each perforation grid cell where well, the density and grid of oil phase fluid are obtained
The active porosity volume of unit calculates the output coefficient M of each perforation grid cell using formula (1)i, calculation formula is：
Mi=ρo,i·so,i·Vi (1)
In formula, MiFor the output coefficient of perforation node unit i, kg；ρo,iFor perforation node unit i oil phase fluid it is close
Degree, kg/m3；so,iFor the oil saturation of perforation node unit i；ViFor the active porosity body of payzone where perforation node unit i
Product, m3。
In step 2, according to producing well all perforation units output coefficient Mi, producing well is calculated using formula (2)
Full well production coefficient, calculation formula is：
In formula, M is the full well production coefficient of producing well, kg；N is the perforation payzone sum of producing well, integer.
In step 3, according to producing well perforation node unit output coefficient MiWith full well production coefficient M, using public affairs
Formula (3) calculates producing well and splits a point coefficient in each payzone, and calculation formula is：
In formula, ciFor producing well a point coefficient, decimal are split in perforation node unit i.
In step 3, in reservoir numerical simulation, set a well a certain given time walk under full well production as Q,
A point coefficient c is split using each perforation node uniti, the production Q i of each analogue perforation node unit is obtained according to formula (4),
Calculation formula is：
Qi=ci×Q (4)。
The production split method suitable for reservoir numerical simulation in the present invention is Reservoir Development technology, research oil
With underground fluids important method of seepage flow characteristics and numerical simulation for oil-gas reservoir technology in reservoir such as natural gas.Production split is more
The important method of layer oil reservoir fine development, each payzone reserves exploitation situation of accurate description and recovery percent of reserves, is reservoir numerical simulation
Whether technology can be with each payzone of accurate description to the important method of contribution to yield degree.This method is directed in numerical simulation for oil-gas reservoir
The problem of production split method, answers on the basis of each payzone of acquisition is averaged oil phase fluid density, oil saturation and effective volume
With production split model, the production split coefficient of each payzone is calculated.This method establishes being capable of accurate description oil reservoir development and oil
The output situation for hiding each payzone in numerical simulation, is of great significance to oil-gas reservoir fine development, to improve numerical reservoir mould
Quasi- technology lays the foundation.
Description of the drawings
Fig. 1 is the flow chart of a specific embodiment of the production split method suitable for reservoir numerical simulation of the present invention.
Specific implementation mode
For enable the present invention above and other objects, features and advantages be clearer and more comprehensible, it is cited below particularly go out preferable implementation
Example, and coordinate shown in attached drawing, it is described in detail below.
As shown in FIG. 1, FIG. 1 is the flow charts of the production split method suitable for reservoir numerical simulation of the present invention.
In a step 101, the oil saturation of each perforation grid cell where well, the density and net of oil phase fluid are obtained
The active porosity volume of lattice unit calculates the output coefficient M of each perforation grid cell using formula (1)i, calculation formula is：
Mi=ρo,i·so,i·Vi (1)
In formula, MiFor the output coefficient of perforation node unit i, kg；ρo,iFor perforation node unit i oil phase fluid it is close
Degree, kg/m3；so,iFor the oil saturation of perforation node unit i；ViFor the active porosity body of payzone where perforation node unit i
Product, m3。
In a step 102, according to producing well obtained above all perforation units output coefficient Mi, using formula
(2) the full well production coefficient of producing well is calculated, calculation formula is：
In formula, M is the full well production coefficient of producing well, kg；N is the perforation payzone sum of producing well, integer.
In step 103, the output coefficient according to producing well obtained above in perforation node unit and full well production system
Number calculates producing well using formula (3) and splits a point coefficient in each payzone, and calculation formula is：
In formula, ciFor producing well a point coefficient, decimal are split in perforation node unit i.
In reservoir numerical simulation, full well production of a well in the case where a certain given time walks is set as Q, utilizes step
The each perforation node unit being calculated in 103 splits a point coefficient ci, each analogue perforation node list is obtained according to formula (4)
The production Q i of member, calculation formula are：
Qi=ci×Q (4)
Flow terminates.
The production split method suitable for reservoir numerical simulation in the present invention, the concentrated expression reserves scale of payzone,
The production status of near wellbore zone and the degree of saturation of fluid, are of great significance to oil-gas reservoir fine development.Meanwhile in oil-gas reservoir
In terms of numerical simulation, suitable production split method can accurately describe the life of well with the output situation of each payzone of accurate description
Occurrence condition improves numerical simulation for oil-gas reservoir accuracy and confidence, lays the foundation to improve numerical simulation for oil-gas reservoir technology.
Claims (5)
1. the production split method suitable for reservoir numerical simulation, which is characterized in that should be suitable for the yield of reservoir numerical simulation
Splitting a point method includes：
Step 1, according to well in the active porosity volume of the oil saturation of each perforation unit, oil phase fluid density and each layer, meter
Calculate the output coefficient of payzone；
Step 2, full well production coefficient when stratified reservoir is developed simultaneously is calculated；
Step 3, a point coefficient is split using what the output coefficient of payzone and full well production coefficient calculated payzone.
2. the production split method according to claim 1 suitable for reservoir numerical simulation, which is characterized in that in step 1
In, obtain oil saturation, the density of oil phase fluid and the active porosity of grid cell of each perforation grid cell where well
Volume calculates the output coefficient M of each perforation grid cell using formula (1)i, calculation formula is：
Mi=ρo,i·so,i·Vi (1)
In formula, MiFor the output coefficient of perforation node unit i, kg；ρo,iFor the density of the oil phase fluid of perforation node unit i, kg/
m3；so,iFor the oil saturation of perforation node unit i；ViFor the active porosity volume of payzone where perforation node unit i, m3。
3. the production split method according to claim 1 suitable for reservoir numerical simulation, which is characterized in that in step 2
In, according to producing well all perforation units output coefficient Mi, the full well production coefficient of producing well is calculated using formula (2),
Calculation formula is：
In formula, M is the full well production coefficient of producing well, kg；N is the perforation payzone sum of producing well, integer.
4. the production split method according to claim 1 suitable for reservoir numerical simulation, which is characterized in that in step 3
In, according to producing well perforation node unit output coefficient MiWith full well production coefficient M, producing well is calculated using formula (3)
A point coefficient is split in each payzone, calculation formula is：
In formula, ciFor producing well a point coefficient, decimal are split in perforation node unit i.
5. the production split method according to claim 4 suitable for reservoir numerical simulation, which is characterized in that in step 3
In, in reservoir numerical simulation, full well production of a well in the case where a certain given time walks is set as Q, utilizes each perforation section
Dot element splits a point coefficient ci, the production Q i of each analogue perforation node unit is obtained according to formula (4), calculation formula is：
Qi=ci×Q (4)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810356826.XA CN108729912B (en) | 2018-04-19 | 2018-04-19 | Yield splitting method suitable for numerical reservoir simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810356826.XA CN108729912B (en) | 2018-04-19 | 2018-04-19 | Yield splitting method suitable for numerical reservoir simulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108729912A true CN108729912A (en) | 2018-11-02 |
| CN108729912B CN108729912B (en) | 2022-04-19 |
Family
ID=63939689
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810356826.XA Active CN108729912B (en) | 2018-04-19 | 2018-04-19 | Yield splitting method suitable for numerical reservoir simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108729912B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109577944A (en) * | 2018-11-28 | 2019-04-05 | 王国亭 | Each layer dynamic control reserve measuring method of multilayer system Low permeability and competent sand gas well |
| CN109630102A (en) * | 2018-11-28 | 2019-04-16 | 郭建林 | Each layer gas-bearing area measuring method of multilayer system Low permeability and competent sand gas well |
| CN110984976A (en) * | 2019-12-20 | 2020-04-10 | 陕西延长石油(集团)有限责任公司研究院 | Yield splitting method for multilayer commingled production gas well |
| CN111382523A (en) * | 2020-03-17 | 2020-07-07 | 中国石油化工股份有限公司 | New method for splitting oil well yield |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104933246A (en) * | 2015-06-15 | 2015-09-23 | 中国海洋石油总公司 | Method for quantitatively representing contribution of gas field reserves parameters to reserves change |
| US20160130916A1 (en) * | 2014-11-06 | 2016-05-12 | Schlumberger Technology Corporation | Local layer geometry engine with work zone generated from buffer defined relative to a wellbore trajectory |
| CN107339087A (en) * | 2017-08-10 | 2017-11-10 | 中国石油天然气股份有限公司 | A kind of water injection rate splits a point method and device |
| CN107795319A (en) * | 2016-08-29 | 2018-03-13 | 中国石油天然气股份有限公司 | Reservoir geology reserves acquisition methods and device with gravitational differentiation phenomenon |
| CN107869348A (en) * | 2017-10-27 | 2018-04-03 | 西北大学 | A kind of method of thick-layer sandstone oil reservoir producing well production split |
-
2018
- 2018-04-19 CN CN201810356826.XA patent/CN108729912B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160130916A1 (en) * | 2014-11-06 | 2016-05-12 | Schlumberger Technology Corporation | Local layer geometry engine with work zone generated from buffer defined relative to a wellbore trajectory |
| CN104933246A (en) * | 2015-06-15 | 2015-09-23 | 中国海洋石油总公司 | Method for quantitatively representing contribution of gas field reserves parameters to reserves change |
| CN107795319A (en) * | 2016-08-29 | 2018-03-13 | 中国石油天然气股份有限公司 | Reservoir geology reserves acquisition methods and device with gravitational differentiation phenomenon |
| CN107339087A (en) * | 2017-08-10 | 2017-11-10 | 中国石油天然气股份有限公司 | A kind of water injection rate splits a point method and device |
| CN107869348A (en) * | 2017-10-27 | 2018-04-03 | 西北大学 | A kind of method of thick-layer sandstone oil reservoir producing well production split |
Non-Patent Citations (1)
| Title |
|---|
| 曹肖萌 等: "薄互储层多层合采条件下的产量劈分新方法研究", 《科学技术与工程》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109577944A (en) * | 2018-11-28 | 2019-04-05 | 王国亭 | Each layer dynamic control reserve measuring method of multilayer system Low permeability and competent sand gas well |
| CN109630102A (en) * | 2018-11-28 | 2019-04-16 | 郭建林 | Each layer gas-bearing area measuring method of multilayer system Low permeability and competent sand gas well |
| CN109630102B (en) * | 2018-11-28 | 2022-11-11 | 郭建林 | Method for measuring gas-containing area of each layer of multilayer low-permeability tight sandstone gas well |
| CN110984976A (en) * | 2019-12-20 | 2020-04-10 | 陕西延长石油(集团)有限责任公司研究院 | Yield splitting method for multilayer commingled production gas well |
| CN110984976B (en) * | 2019-12-20 | 2023-02-28 | 陕西延长石油(集团)有限责任公司研究院 | Yield splitting method for multilayer commingled production gas well |
| CN111382523A (en) * | 2020-03-17 | 2020-07-07 | 中国石油化工股份有限公司 | New method for splitting oil well yield |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108729912B (en) | 2022-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108729912A (en) | Production split method suitable for reservoir numerical simulation | |
| CN107829718B (en) | Oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory | |
| CN105224704B (en) | A kind of evaluation method of seabed erosion | |
| CN102041995B (en) | System for monitoring complicated oil deposit flooding conditions | |
| CN104632200B (en) | Fluvial facies multilayer and thick-layer oil reservoir fine modeling method | |
| CN103485770B (en) | A kind of method and system asking for oil saturation based on artificial neural network | |
| CN108089227B (en) | Novel stratum pore pressure prediction method based on three-dimensional seismic data | |
| CN109670220A (en) | A kind of horizontal well air water two-phase method for numerical simulation based on unstrctured grid | |
| CN105808311B (en) | A kind of reservoir modeling Fast Fitting method based on Dimension-reduced Strategy | |
| CN110056346B (en) | Oil reservoir three-dimensional original water saturation simulation method based on trend change function | |
| CN110838175B (en) | Geological model building method for gas injection development oil reservoir | |
| CN107991714A (en) | The quantitative approach recovered based on lake basin paleotopography | |
| CN106150477A (en) | A kind of method determining single well controlled reserves | |
| CN110439515A (en) | Note adopts parameter optimization method and device | |
| CN109684685B (en) | method for analyzing productivity and reservoir stability under pressure-reducing exploitation condition of hydrate in multilateral well | |
| CN105095642B (en) | History matching adjusting parameter determines method and device in reservoir numerical simulation | |
| CN105205273B (en) | The analogy method and its device that oil gas relay is flowed in compact reservoir multi-dielectric | |
| CN110263434A (en) | A kind of flow unit method for numerical simulation based on multiple dimensioned mixed finite element | |
| CN108518219A (en) | It is a kind of that the evaluation method of well is selected based on the natural energy scale extract that develops oil fields | |
| CN109958413A (en) | A kind of ultra-high water cut stage Reservoir behavior flow unit division methods | |
| CN106321057B (en) | Fractured horizontal well crack well pattern design methods | |
| CN108457630A (en) | The method for improving oil recovery using polygamy point small pressure difference water filling | |
| CN104156617B (en) | For the six stage modeling methods that multilayer sandstone reservoirs gas-bearing formation attribute classification is characterized | |
| CN104183018B (en) | For the six stage modeling methods characterized with aqueous carbonate rock salt gas reservoir water-gas distribution | |
| Zhao et al. | INSIM-FPT-3D: A data-driven model for history matching, water-breakthrough prediction and well-connectivity characterization in three-dimensional reservoirs |
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 |