CN107290259B - The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs - Google Patents

The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs Download PDF

Info

Publication number
CN107290259B
CN107290259B CN201610191329.XA CN201610191329A CN107290259B CN 107290259 B CN107290259 B CN 107290259B CN 201610191329 A CN201610191329 A CN 201610191329A CN 107290259 B CN107290259 B CN 107290259B
Authority
CN
China
Prior art keywords
pore throat
radius
throat radius
effective
seepage flow
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
Application number
CN201610191329.XA
Other languages
Chinese (zh)
Other versions
CN107290259A (en
Inventor
殷艳玲
孙志刚
王军
房会春
荣毅
孙强
李培伦
李新
罗惠平
杨伟宇
陈亚宁
张玉利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201610191329.XA priority Critical patent/CN107290259B/en
Publication of CN107290259A publication Critical patent/CN107290259A/en
Application granted granted Critical
Publication of CN107290259B publication Critical patent/CN107290259B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials

Abstract

The present invention provides a kind of calculation method of effective seepage flow pore throat radius of low permeability sandstone reservoirs, and the calculation method of the effective seepage flow pore throat radius of the low permeability sandstone reservoirs includes: step 1, and the capillary pressure curve of testing rock core obtains reservoir micropore structure distributed data;Step 2, single pore throat is calculated to the contribution rate of permeability according to the distributed data of pore structure;Step 3, each pore throat radius is subtracted into pore throat radius lower limit value, obtains the distributed data of effective radius;Step 4, according to the distribution of effective pore throat data, effective seepage flow pore throat radius is calculated.The calculation method of the effective seepage flow pore throat radius of the low permeability sandstone reservoirs can quantitatively calculate the effective radius of the practical seepage flow of low permeability reservoir, and for the percolation ability of rationally assessment LOW PERMEABILITY RESERVOIR, effective exploitation lays the foundation.

Description

The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs
Technical field
The present invention relates to porous media features to the theoretical research technology neck of the underground fluids seepage effect such as petroleum, natural gas Domain especially relates to a kind of calculation method of effective seepage flow pore throat radius of low permeability sandstone reservoirs.
Background technique
When flowing of the Study of Fluid in porous media, it is necessary to consider the influence of porous media.In porous media spy In sign, it is most important that its pore structure characteristic, the i.e. shape, size of pore throat, distribution and its interconnected relationship.Wherein pore throat Size is the most common characterization parameter of pore structure.Since the pore structure of low permeability reservoir belongs to the thin larynx type of aperture, permeability Contribution margin is mainly provided by larger pore throat, and therefore, in low permeability reservoir there are a part of pore throats does not contribute permeability.Cause This, mean radius the problem of cannot sufficiently reflecting reservoir fluid flow radius of currently used characterization Pore throat size.For this purpose, I From low permeability reservoir pore structure characteristic, consider that liquid acts on admittedly, invented a kind of meter of low permeability sandstone reservoirs effective radius Calculation method solves the above technical problem.
Summary of the invention
The object of the present invention is to provide a kind of effective radiuses that can quantitatively calculate the practical seepage flow of low permeability reservoir, rationally to comment Estimate the percolation ability of LOW PERMEABILITY RESERVOIR, the calculating side of the effective seepage flow pore throat radius of the low permeability sandstone reservoirs that effective exploitation lays the foundation Method.
The purpose of the present invention can be achieved by the following technical measures: the meter of the effective seepage flow pore throat radius of low permeability sandstone reservoirs Calculation method, the calculation method of the effective seepage flow pore throat radius of the low permeability sandstone reservoirs include: step 1, the capillary pressure of testing rock core Curve obtains reservoir micropore structure distributed data;Step 2, single pore throat pair is calculated according to the distributed data of pore structure The contribution rate of permeability;Step 3, each pore throat radius is subtracted into pore throat radius lower limit value, obtains the distributed data of effective radius; Step 4, according to the distribution of effective pore throat data, effective seepage flow pore throat radius is calculated.
The purpose of the present invention can be also achieved by the following technical measures:
In step 1, formula is utilizedCalculate pore throat radius, in formula, riFor rock core pore throat radius, μm;PciFor Into mercury pressure, MPa;Utilize formulaCalculate throat distribution frequency, in formula, fiFor distribution frequency, decimal;ΔsiFor phase The difference of adjacent mercury saturation degree, percentage;smaxIt is maximum into mercury saturation degree, percentage.
In step 2, by formulaSingle pore throat is calculated to the contribution rate of permeability, in formula, KJIt is a certain The permeability contribution that pore throat radius can be provided, decimal;riFor rock core pore throat radius, μm;fiFor distribution frequency, decimal.
It in step 2, adds up to all pore throat contribution rates, pore throat radius is done into figure with corresponding accumulation contribution rate, is used in combination Binomial relationship is fitted, and obtains the relational expression K of accumulation contribution rate and pore throat radiusJ=ar2+ br+c, in formula, KJFor The permeability contribution that a certain pore throat radius can be provided, decimal;R is rock core pore throat radius, μm;A, b, c are constant, utilize this Relational expression seeks corresponding pore throat radius lower limit value r when permeability accumulation contribution rate is 98%98
In step 4, formula R is utilizedy=∑ ryi·fiCalculate effective seepage flow pore throat radius, in formula, RyFor effective seepage flow pore throat Radius, μm;ryiFor rock core pore throat effective radius, μm;fiFor distribution frequency, decimal.
The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs in the present invention is related to microscopic seepage experiment point Analysis is the important method for studying underground fluids percolation laws in porous media such as petroleum gas.With Exploitation Level of Oilfield Continuous improvement, further increase the recovery percent of reserves of low permeability reservoir, be guarantee oil field sustainable and stable development important directions.The party Method considers that liquid acts on admittedly from low permeability reservoir pore structure characteristic, not for currently used mainstream radius and mean radius The problem of capable of sufficiently reflecting reservoir Test Liquid Permeability of Core, establishes the calculation method for being suitable for the effective seepage flow pore throat radius of LOW PERMEABILITY RESERVOIR, Preferably to characterize the effective percolation ability of low permeability reservoir.The effective seepage flow pore throat radius calculation method of the low permeability sandstone reservoirs can determine Scale sign environmental liquids enter actual effective percolation ability after reservoir, lay the foundation for the effective exploitation of LOW PERMEABILITY RESERVOIR.
Detailed description of the invention
Fig. 1 is the stream of a specific embodiment of the calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs of the invention Cheng Tu;
Fig. 2 is the relational graph of cumulative contribution rate and pore throat radius in a specific embodiment of the invention.
Specific embodiment
To enable above and other objects, features and advantages of the invention to be clearer and more comprehensible, preferably implementation is cited below particularly out Example, and cooperate shown in attached drawing, it is described in detail below.
As shown in FIG. 1, FIG. 1 is the processes of the calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs of the invention Figure.
Step 101, the capillary pressure curve of testing rock core obtains reservoir micropore structure distributed data, utilizes formulaCalculate pore throat radius, in formula, riFor rock core pore throat radius, μm;PciFor into mercury pressure, MPa;Utilize formulaCalculate throat distribution frequency, in formula, fiFor distribution frequency, decimal;ΔsiFor the difference of adjacent mercury saturation degree, percentage Number;smaxIt is maximum into mercury saturation degree, percentage.
Step 102, formula is pressed according to the distributed data of pore structureSingle pore throat is calculated to permeability Contribution rate, in formula, KJFor the permeability contribution that a certain pore throat radius can be provided, decimal;riFor rock core pore throat radius, μm;fi For distribution frequency, decimal.It adds up to all pore throat contribution rates, pore throat radius is done into figure with corresponding accumulation contribution rate, is used in combination Binomial relationship is fitted, and obtains the relational expression K of accumulation contribution rate and pore throat radiusJ=ar2+ br+c, in formula, KJFor The permeability contribution that a certain pore throat radius can be provided, decimal;R is rock core pore throat radius, μm;A, b, c are constant.Utilize this Relational expression seeks corresponding pore throat radius lower limit value r when permeability accumulation contribution rate is 98%98
Step 103, by each pore throat radius riSubtract pore throat radius lower limit value r98, obtain effective radius ryiDistribution number According to;
Step 104, according to the distribution of effective pore throat data, formula R is utilizedy=∑ ryi·fiCalculate effective seepage flow pore throat half Diameter, in formula, RyFor effective seepage flow pore throat radius, μm;ryiFor rock core pore throat effective radius, μm;fiFor distribution frequency, decimal.
Illustrate the specific technical solution of the invention by taking the calculating of Shengli Oil Field somewhere block rock sample effective radius as an example:
Step 1, the capillary pressure curve of testing rock core obtains micropore structure distributed data (as shown in table 1);Pore throat The calculation formula of radius are as follows:
Wherein: riFor rock core pore throat radius, μm;PciFor into mercury pressure, MPa.
The calculation formula of distribution frequency are as follows:
Wherein: fiFor distribution frequency, decimal;ΔsiFor the difference of adjacent mercury saturation degree, percentage;smaxIt is full into mercury for maximum And degree, percentage.
1 capillary pressure test data of table
Step 2, single pore throat is calculated to the contribution rate of permeability and right by formula (3) according to the distributed data of pore structure All pore throat contribution rates add up, and calculating process is as shown in table 2;
In formula: KJFor the permeability contribution that a certain pore throat radius can be provided, decimal.
Contribution rate calculating process of 2 pore throat of table to permeability
Cumulative contribution rate and pore throat radius are done into figure, such as Fig. 2, obtain the relation formula (4) of the two, is obtained by formula (4) KJCorresponding pore throat radius lower limit value is 0.369 μm when=0.98;
KJ=-0.5239r2+0.1562·r+0.992 (4)
Step 3, pore throat radius each in rock core micropore structure is subtracted to the contributive pore throat radius of permeability 0.369 μm of lower limit value, the distributed data for obtaining effective pore throat is as shown in table 3;
The effective throat distribution of table 3
Step 4, according to the distribution of effective pore throat data, effective seepage flow pore throat radius is calculated.
Ry=0.3691 × 0.3980+0.1247 × 0.1354=0.164 μm.

Claims (5)

1. the calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs, which is characterized in that the low permeability sandstone reservoirs are effectively seeped Stream pore throat radius calculation method include:
Step 1, the capillary pressure curve of testing rock core obtains reservoir micropore structure distributed data;
Step 2, single pore throat is calculated to the contribution rate of permeability according to the distributed data of pore structure;
Step 3, each pore throat radius is subtracted into pore throat radius lower limit value, obtains the distributed data of effective radius;
Step 4, according to the distribution of effective pore throat data, effective seepage flow pore throat radius is calculated.
2. the calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs according to claim 1, which is characterized in that Step 1, formula is utilizedCalculate pore throat radius, in formula, riFor rock core pore throat radius, μm;PciFor into mercury pressure, MPa;Utilize formulaCalculate throat distribution frequency, in formula, fiFor distribution frequency, decimal;ΔsiFor adjacent mercury saturation degree Difference, percentage;smaxIt is maximum into mercury saturation degree, percentage.
3. the calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs according to claim 1, which is characterized in that Step 2, by formulaSingle pore throat is calculated to the contribution rate of permeability, in formula, KJFor a certain pore throat radius institute The permeability contribution that can be provided, decimal;riFor rock core pore throat radius, μm;fiFor distribution frequency, decimal.
4. the calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs according to claim 3, which is characterized in that Step 2, it adds up to all pore throat contribution rates, pore throat radius and corresponding accumulation contribution rate is done into figure, and with binomial relationship It is fitted, obtains the relational expression K of accumulation contribution rate and pore throat radiusJ=ar2+ br+c, in formula, KJFor a certain pore throat half The permeability contribution that diameter can be provided, decimal;R is rock core pore throat radius, μm;A, b, c are constant, are sought using this relational expression Permeability accumulates corresponding pore throat radius lower limit value r when contribution rate is 98%98
5. the calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs according to claim 1, which is characterized in that Step 4, formula R is utilizedy=∑ ryi·fiCalculate effective seepage flow pore throat radius, in formula, RyFor effective seepage flow pore throat radius, μm; ryiFor rock core pore throat effective radius, μm;fiFor distribution frequency, decimal.
CN201610191329.XA 2016-03-30 2016-03-30 The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs Active CN107290259B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610191329.XA CN107290259B (en) 2016-03-30 2016-03-30 The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610191329.XA CN107290259B (en) 2016-03-30 2016-03-30 The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs

Publications (2)

Publication Number Publication Date
CN107290259A CN107290259A (en) 2017-10-24
CN107290259B true CN107290259B (en) 2019-09-10

Family

ID=60087114

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610191329.XA Active CN107290259B (en) 2016-03-30 2016-03-30 The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs

Country Status (1)

Country Link
CN (1) CN107290259B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108303348B (en) * 2017-12-20 2020-07-14 中国石油化工股份有限公司 Device and system for simulating tensile degradation characteristics of polymer in oil reservoir pore throat
CN108959723B (en) * 2018-06-11 2023-05-30 东北石油大学 Method for discussing permeability main control factors of tight sandstone reservoir
CN109100278B (en) * 2018-07-18 2020-02-14 西南石油大学 Apparent permeability calculation method considering shale pore size distribution characteristics
CN110186843B (en) * 2019-04-24 2020-03-20 中国矿业大学 Device and method for measuring effective influence radius of high-pressure medium injected into drill hole
CN112326523A (en) * 2020-09-18 2021-02-05 中国石油天然气股份有限公司 Method and device for acquiring pore throat characteristic parameters, computer equipment and storage medium
CN112966365B (en) * 2021-02-04 2023-09-12 中海石油(中国)有限公司 Method for evaluating reverse condensation injury of ultralow condensation gas reservoir

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706913A (en) * 2012-06-19 2012-10-03 中国石油天然气股份有限公司 Method and device for determining radius lower limit of reservoir formation pore throat of compact sandstone reservoir
CN103675236A (en) * 2013-12-19 2014-03-26 中国石油天然气股份有限公司 Method and device for testing critical tight oil filling pore throat radius threshold
CN104453882A (en) * 2014-12-10 2015-03-25 中国石油天然气股份有限公司 Method and device for determining MOIP (movable oil in place) of crude oil of tight reservoir
CN104777181A (en) * 2014-12-31 2015-07-15 中国石油天然气股份有限公司 Method for determining nuclear magnetic resonance T2 cutoff value and fluid saturation of tight oil, and apparatus thereof

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706913A (en) * 2012-06-19 2012-10-03 中国石油天然气股份有限公司 Method and device for determining radius lower limit of reservoir formation pore throat of compact sandstone reservoir
CN103675236A (en) * 2013-12-19 2014-03-26 中国石油天然气股份有限公司 Method and device for testing critical tight oil filling pore throat radius threshold
CN104453882A (en) * 2014-12-10 2015-03-25 中国石油天然气股份有限公司 Method and device for determining MOIP (movable oil in place) of crude oil of tight reservoir
CN104777181A (en) * 2014-12-31 2015-07-15 中国石油天然气股份有限公司 Method for determining nuclear magnetic resonance T2 cutoff value and fluid saturation of tight oil, and apparatus thereof

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Pore-throat radius and tortuosity estimation from formation resistivity data for tight-gas sandstone reservoirs";Ali S. Ziarani 等;《Journal of Applied Geophysics》;20120322;第83卷;第65-73页
"低渗砂岩与微裂缝泥灰岩储集层毛细管压力曲线特征对比";臧士宾 等;《新疆石油地质》;20111031;第32卷(第5期);第504-507页
"孔隙结构研究在低渗透储层有效性评价中的应用";耿斌 等;《断块油气田》;20110331;第18卷(第2期);第187-190页

Also Published As

Publication number Publication date
CN107290259A (en) 2017-10-24

Similar Documents

Publication Publication Date Title
CN107290259B (en) The calculation method of the effective seepage flow pore throat radius of low permeability sandstone reservoirs
Guo et al. Study on gas flow through nano pores of shale gas reservoirs
CN106599449B (en) Well test interpretation method for karst cave volume calculation
CN105574320B (en) The evaluation method of the effective percolation ability of low permeability sandstone reservoirs
CN110598167B (en) Processing method of oil-water relative permeability experimental data of low-permeability reservoir
CN108798634A (en) A kind of fractured-vuggy reservoir list hole-mean value of single layer Well Test Data Analysis Method
Guo et al. Study on gas permeability in nano pores of shale gas reservoirs
Yanze et al. The role of diffusion for nonequilibrium gas injection into a fractured reservoir
Xu et al. An analysis of the radial flow in the heterogeneous porous media based on fractal and constructal tree networks
Li et al. An improved method to study CO 2–oil relative permeability under miscible conditions
You et al. Pore structure and limit pressure of gas slippage effect in tight sandstone
CN106547930A (en) Consider the gas drainage radius computational methods of tight gas reservoir seepage flow mechanism
CN104200105A (en) Method for determining filling property lower limit of tight sandstone gas
Meng et al. Effect of fluid viscosity on correlation of oil recovery by linear counter-current spontaneous imbibition
Liu et al. Improving production in gas/condensate reservoirs by wettability alteration to gas wetness
Wei et al. Mechanism model for shale gas transport considering diffusion, adsorption/desorption and Darcy flow
CN106769745B (en) The characterizing method that monophasic fluid passes through the non linear fluid flow through porous medium feature of flow in low permeability core
Wang et al. Performance of drainage experiments with Orinoco Belt heavy oil in a long laboratory core in simulated reservoir conditions
Zhang et al. Analytical method for performance evaluation of fractured horizontal wells in tight reservoirs
Ibarra et al. Experimental study of a poor boy downhole gas separator under continuous gas-liquid flow
Yang et al. Pore shape factors in shale: calculation and impact evaluation on fluid imbibition
Jamiolahmady et al. Gas/Condensate Relative Permeabilities in Propped-Fracture Porous Media: Coupling vs. Inertia
Poe et al. Novel methodology to estimate reservoir pressure and productivity index in unconventional and conventional reservoirs using production data
CN203772691U (en) Ultra-low permeability core gas permeability testing device
Chenggang et al. Tubing Pressure Characteristics and Fracturing Response Analysis of Horizontal Well in Tight Reservoir An example from three horizontal wells in Changji oilfield of eastern Junggar basin

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