CN113124808B - Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir - Google Patents

Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir Download PDF

Info

Publication number
CN113124808B
CN113124808B CN202110389388.9A CN202110389388A CN113124808B CN 113124808 B CN113124808 B CN 113124808B CN 202110389388 A CN202110389388 A CN 202110389388A CN 113124808 B CN113124808 B CN 113124808B
Authority
CN
China
Prior art keywords
radius
microsphere
pore throat
frequency distribution
distribution 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
Application number
CN202110389388.9A
Other languages
Chinese (zh)
Other versions
CN113124808A (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.)
Shaanxi Yanchang Petroleum Group Co Ltd
Original Assignee
Shaanxi Yanchang Petroleum Group Co Ltd
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 Shaanxi Yanchang Petroleum Group Co Ltd filed Critical Shaanxi Yanchang Petroleum Group Co Ltd
Priority to CN202110389388.9A priority Critical patent/CN113124808B/en
Publication of CN113124808A publication Critical patent/CN113124808A/en
Application granted granted Critical
Publication of CN113124808B publication Critical patent/CN113124808B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/10—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring diameters
    • 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
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Abstract

The invention relates to a method for evaluating the matching of an initial radius of a microsphere and a pore throat structure of a reservoir, which comprises the following steps: (1) respectively measuring the pore throat radius and the initial microsphere radius of the reservoir core to obtain a reservoir core pore throat radius frequency distribution curve, a microsphere radius frequency distribution curve, a reservoir core pore throat radius accumulated frequency distribution curve and a microsphere radius accumulated frequency distribution curve; (2) sequentially determining r on the cumulative frequency distribution curve of the pore throat radius of the reservoir rock core and the cumulative frequency distribution curve of the microsphere radius min 、r 25 、r 50 、r 75 、r max And R min 、R 25 、R 50 、R 75 、R max (ii) a (3) Fitting Q on frequency distribution curve of pore throat radius of reservoir core 1 (R)、Q 2 (R)、Q 3 (R) and Q 4 (R); fitting f on the microsphere radius frequency distribution curve 1 (r)、f 2 (r)、f 3 (r) and f 4 (r); (4) and evaluating the matching of the initial radius of the microsphere and the pore throat structure of the reservoir. The method provided by the invention avoids long-period microsphere expansibility experiments, is simple, convenient, efficient and quick, and can quickly screen out a microsphere system which meets physical property conditions of a target reservoir.

Description

Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir
Technical Field
The invention belongs to the technical field of oil exploitation, and relates to a method for quantitatively and quickly evaluating the matching of an initial radius of a microsphere and a pore throat structure of a reservoir.
Background
As oil field development progresses, water flooding has become the most widely used technology at present. Because most of the Ordos basins are fractured low-permeability oil reservoirs, natural micro fractures of the reservoirs and artificial fractures generated by fracturing modification commonly exist, the self-heterogeneity of the oil reservoirs is severe, a dominant channeling channel is easy to form in the reservoirs due to long-term water injection development, and once the channeling channel is formed, water flooding wave and efficiency are rapidly reduced.
The foam, microsphere and gel profile control plugging is the most common technical method for treating the water channeling problem, compared with other plugging agents, the microsphere has the characteristics of swelling, elastic deformation and the like, so that the channeling channels with different sizes can be plugged in a self-adjusting mode, the swept volume of a residual oil enrichment area is further improved, and the technology is favored by various oil fields in recent years. The research on the matching of the radius of the microsphere and the pore throat structure of the target reservoir stratum is a necessary step before the injection test of the microsphere mine field and is also a key link influencing the later effect, however, the evaluation on the matching of the microsphere and the pore throat structure of the reservoir stratum at the present stage generally firstly tests the radius of the microsphere after expansion stability through a long-period experiment, and then simply compares the average diameters of the microsphere and the pore throat structure to determine the radius, and cannot directly carry out quick judgment through the initial radius parameters of the microsphere sample.
Disclosure of Invention
The invention aims to solve the problems and provides a method for evaluating the matching between the initial radius of a microsphere and the pore throat structure of a reservoir.
The technical scheme of the invention is as follows:
a method for evaluating the matching of the initial radius of a microsphere and a pore throat structure of a reservoir comprises the following steps:
(1) respectively measuring the pore throat radius and the initial microsphere radius of the reservoir core to obtain a reservoir core pore throat radius frequency distribution curve, a microsphere radius frequency distribution curve, a reservoir core pore throat radius accumulated frequency distribution curve and a microsphere radius accumulated frequency distribution curve;
(2) sequentially determining r on the cumulative frequency distribution curve of the pore throat radius of the reservoir rock core and the cumulative frequency distribution curve of the microsphere radius min 、r 25 、r 50 、r 75 、r max And R min 、R 25 、R 50 、R 75 、R max ；
(3) Fitting Q on frequency distribution curve of pore throat radius of reservoir core 1 (R)、Q 2 (R)、Q 3 (R) and Q 4 (R); fitting f on the microsphere radius frequency distribution curve 1 (r)、f 2 (r)、f 3 (r) and f 4 (r)；
(4) Evaluating the matching of the initial radius of the microsphere and the pore throat structure of the reservoir by using a formula (1), a formula (2), a formula (3) and a formula (4);
Figure GDA0003642384200000021
Figure GDA0003642384200000022
Figure GDA0003642384200000023
Figure GDA0003642384200000024
when the formula (1), the formula (2), the formula (3) and the formula (4) are all established, the initial radius of the microsphere is better matched with the pore throat structure of the target stratum; otherwise, the two are not matched;
wherein r is 25 、r 50 、r 75 And r max Respectively representing the radiuses of 25%, 50%, 75% and 100% of the corresponding accumulated frequency on the distribution curve of the accumulated frequency of the radius of the microsphere; r 25 、R 50 、R 75 And R max Respectively representing the radiuses of 25%, 50%, 75% and 100% of corresponding accumulated frequency on the cumulative frequency distribution curve of the core pore throat radius; r is min And R min The minimum radius values of the microspheres and pore throats, respectively; f. of 1 (r)、f 2 (r)、f 3 (r) and f 4 (r) represents a radius of the microsphere, r min ～r 25 、r 25 ～r 50 、r 50 ～r 75 And r 75 ～r max Corresponding to the radius of the microsphere, accumulating the frequency distribution fitting function; q 1 (R)、Q 2 (R)、Q 3 (R) and Q 4 (R) respectively represents that the radius of the pore throat of the reservoir core is R min ～R 25 、R 25 ～R 50 、R 50 ～R 75 And R 75 ～R max And fitting a function of cumulative frequency distribution of pore throat radius of the reservoir core corresponding to the interval.
The invention has the technical effects that:
the method provided by the invention is used for evaluating the matching of the initial radius of the microsphere and the pore throat structure of the reservoir, so that a long-period microsphere expansibility experiment is avoided, the method is simple, convenient, efficient and rapid, and a microsphere system meeting the physical property condition of a target reservoir can be rapidly screened out.
Drawings
FIG. 1 is a cumulative frequency distribution curve of pore throat radius and a cumulative frequency distribution curve of microsphere radius for a reservoir core according to the present invention.
FIG. 2 is a frequency distribution curve of pore throat radius of a reservoir core.
FIG. 3 is a frequency distribution curve of the radius of the microsphere.
FIG. 4 is a graph showing the change in radius of microspheres at different time conditions.
Detailed Description
A method for evaluating the matching of the initial radius of a microsphere and a pore throat structure of a reservoir comprises the following steps:
(1) respectively measuring the pore throat radius and the initial microsphere radius of the reservoir core to obtain a reservoir core pore throat radius frequency distribution curve (figure 2), a microsphere radius frequency distribution curve (figure 3), a reservoir core pore throat radius accumulated frequency distribution curve and a microsphere radius accumulated frequency distribution curve (figure 1);
(2) sequentially determining r on the cumulative frequency distribution curve of the pore throat radius of the reservoir rock core and the cumulative frequency distribution curve of the microsphere radius min 、r 25 、r 50 、r 75 、r max And R min 、R 25 、R 50 、R 75 、R max ；r min ＝1.41μm、r 25 ＝4.49μm、r 50 ＝8.20μm、r 75 ＝15.64μm、r max ＝115.79μm；R min ＝1.26μm、R 25 ＝5.66μm、R 50 ＝10.15μm、R 75 ＝20.09μm、R max ＝123.55μm；
(3) Fitting Q on frequency distribution curve of pore throat radius of reservoir core 1 (R)、Q 2 (R)、Q 3 (R) and Q 4 (R); fitting f on the microsphere radius frequency distribution curve 1 (r)、f 2 (r)、f 3 (r) and f 4 (r); the method comprises the following specific steps:
Q 1 (R)＝-0.071R 3 +0.603R 2 -0.075R-0.817
Q 2 (R)＝0.013R 3 -0.398R 2 +3.974R-6.979
Q 3 (R)＝0.003R 3 -0.147R 2 +1.961R-2.03
Q 4 (R)＝681.8R -1.64
f 1 (r)＝-0.135r 2 +2.424r-3.171
f 2 (r)＝0.025r 3 -0.647r 2 +5.463r-8.807
f 3 (r)＝-0.005r 2 -0.058r+7.842
f 4 (r)＝503.5r -1.53
(4) r determined above min 、r 25 、r 50 、r 75 、r max And R min 、R 25 、R 50 、R 75 、R max And fitted Q 1 (R)、Q 2 (R)、Q 3 (R) and Q 4 (R) and f 1 (r)、f 2 (r)、f 3 (r) and f 4 (r) substituting the formula (1), the formula (2), the formula (3) and the formula (4) to evaluate the matching of the initial radius of the microsphere and the pore throat structure of the reservoir;
(5) respectively calculate
Figure GDA0003642384200000031
Figure GDA0003642384200000041
Figure GDA0003642384200000042
Figure GDA0003642384200000043
The constraint conditions of formula (1), formula (2), formula (3) and formula (4) are respectively satisfied. The evaluation results show that: the microsphere system has good matching property with a target reservoir pore throat structure.
Example verification currently, the matching relationship between the initial radius of the microsphere and the pore-throat structure of the reservoir is evaluated by adopting a formula (5):
Figure GDA0003642384200000044
wherein δ is a radius matching coefficient; d 50 The average radius of the microspheres after expansion stabilization; d 50 The average pore throat diameter of the reservoir core is shown. When the radii are matchedWhen the coefficient delta is 1.4 to 1.5, the matching property between the two is considered to be good.
The microspheres remained stable in radius after 3 days of expansion (FIG. 4), with a radius of about 29.10 μm. The average pore throat radius of the target reservoir is 10.15 mu m, and the matching calculation of the microsphere radius and the pore throat scale of the target reservoir is carried out by using the formula (5), so that the calculation result is 1.43. Also shows that the radius is better matched with the pore throat structure of the target reservoir.

Claims (1)

1. A method for evaluating the matching of the initial radius of a microsphere and a pore throat structure of a reservoir is characterized by comprising the following steps: the method comprises the following steps:
(1) respectively measuring the pore throat radius and the initial microsphere radius of the reservoir core to obtain a reservoir core pore throat radius frequency distribution curve, a microsphere radius frequency distribution curve, a reservoir core pore throat radius accumulated frequency distribution curve and a microsphere radius accumulated frequency distribution curve;
(2) sequentially determining r on the cumulative frequency distribution curve of the pore throat radius of the reservoir rock core and the cumulative frequency distribution curve of the microsphere radius min 、r 25 、r 50 、r 75 、r max And R min 、R 25 、R 50 、R 75 、R max ；
(3) Fitting Q on frequency distribution curve of pore throat radius of reservoir core 1 (R)、Q 2 (R)、Q 3 (R) and Q 4 (R); fitting f on the microsphere radius frequency distribution curve 1 (r)、f 2 (r)、f 3 (r) and f 4 (r)；
(4) Evaluating the matching of the initial radius of the microsphere and the pore throat structure of the reservoir by using a formula (1), a formula (2), a formula (3) and a formula (4);
Figure FDA0003642384190000011
Figure FDA0003642384190000012
Figure FDA0003642384190000013
Figure FDA0003642384190000014
when the formula (1), the formula (2), the formula (3) and the formula (4) are all established, the initial radius of the microsphere is better matched with the pore throat structure of the target stratum; otherwise, the two are not matched;
wherein r is 25 、r 50 、r 75 And r max Respectively representing the radiuses of 25%, 50%, 75% and 100% of the corresponding accumulated frequency on the distribution curve of the accumulated frequency of the radius of the microsphere; r 25 、R 50 、R 75 And R max Respectively representing the radiuses of 25%, 50%, 75% and 100% of corresponding accumulated frequency on the cumulative frequency distribution curve of the core pore throat radius; r is min And R min The minimum radius values of the microspheres and pore throats, respectively; f. of 1 (r)、f 2 (r)、f 3 (r) and f 4 (r) represents a radius of the microsphere, r min ～r 25 、r 25 ～r 50 、r 50 ～r 75 And r 75 ～r max Corresponding to the radius of the microsphere, accumulating the frequency distribution fitting function; q 1 (R)、Q 2 (R)、Q 3 (R) and Q 4 (R) respectively represents that the radius of the pore throat of the reservoir core is R min ～R 25 、R 25 ～R 50 、R 50 ～R 75 And R 75 ～R max And fitting a function of cumulative frequency distribution of pore throat radius of the reservoir core corresponding to the interval.
CN202110389388.9A 2021-04-12 2021-04-12 Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir Active CN113124808B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110389388.9A CN113124808B (en) 2021-04-12 2021-04-12 Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110389388.9A CN113124808B (en) 2021-04-12 2021-04-12 Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir

Publications (2)

Publication Number Publication Date
CN113124808A CN113124808A (en) 2021-07-16
CN113124808B true CN113124808B (en) 2022-08-12

Family

ID=76776006

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110389388.9A Active CN113124808B (en) 2021-04-12 2021-04-12 Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir

Country Status (1)

Country Link
CN (1) CN113124808B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115267141B (en) * 2022-09-26 2023-01-24 中国石油天然气股份有限公司 Reservoir comprehensive evaluation method and device

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6008645A (en) * 1997-03-11 1999-12-28 Conoco Inc. Prediction of permeability from capillary pressure curves derived from nuclear magnetic resonance pore size distributions
JP2004184721A (en) * 2002-12-04 2004-07-02 Canon Inc Manufacturing method of dry type toner, and dry type toner
CN104034645A (en) * 2014-06-07 2014-09-10 向丹 Detection algorithm for pore structure of shale gas reservoir
CN105334149A (en) * 2015-11-24 2016-02-17 中国石油大学(北京) Micro-pore structure evaluation and reservoir classification method for tight reservoirs
CN106383078A (en) * 2016-09-20 2017-02-08 中国石油天然气股份有限公司 Method and device for determining waterflood efficiency of rock
CN108507920A (en) * 2018-03-02 2018-09-07 中国石油天然气股份有限公司 A kind of analogy method of water-oil phase relative permeability
CN109856019A (en) * 2019-03-04 2019-06-07 东南大学 A kind of measurement method of particle system particle diameter distribution

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA1250644A (en) * 1984-12-31 1989-02-28 Robert Ehrlich Method of making rock-pore micromodel involving generation of an accurate and reliable template image of an actual reservoir rock pore system

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6008645A (en) * 1997-03-11 1999-12-28 Conoco Inc. Prediction of permeability from capillary pressure curves derived from nuclear magnetic resonance pore size distributions
JP2004184721A (en) * 2002-12-04 2004-07-02 Canon Inc Manufacturing method of dry type toner, and dry type toner
CN104034645A (en) * 2014-06-07 2014-09-10 向丹 Detection algorithm for pore structure of shale gas reservoir
CN105334149A (en) * 2015-11-24 2016-02-17 中国石油大学(北京) Micro-pore structure evaluation and reservoir classification method for tight reservoirs
CN106383078A (en) * 2016-09-20 2017-02-08 中国石油天然气股份有限公司 Method and device for determining waterflood efficiency of rock
CN108507920A (en) * 2018-03-02 2018-09-07 中国石油天然气股份有限公司 A kind of analogy method of water-oil phase relative permeability
CN109856019A (en) * 2019-03-04 2019-06-07 东南大学 A kind of measurement method of particle system particle diameter distribution

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Shuai Zhao等.A comprehensive investigation of polymer microspheres (PMs) migration in porous media: EOR implication.《Fuel》.2018,第249–258页. *
梁守成等.聚合物微球粒径与岩芯孔喉的匹配关系研究.《西南石油大学（自然科学版）》.2016,第38卷(第1期),第140-145页. *

Also Published As

Publication number Publication date
CN113124808A (en) 2021-07-16

Similar Documents

Publication Publication Date Title
US7971644B2 (en) Fluid loss additive and breaking agent
US8838427B2 (en) Method for determining the closure pressure of a hydraulic fracture
CN113124808B (en) Method for evaluating matching of initial radius of microsphere and pore throat structure of reservoir
CN110984949B (en) Shale continuous sand-adding fracturing process
Chuanliang et al. A study on mechanism of oil displacement by imbibition
CN110094201A (en) A kind of polymer microballoon deep profile correction injection partial size selection method
CN111058824B (en) Temporary plugging and shunting acidification method for strong heterogeneous sandstone reservoir through sieve tube
CN114152554B (en) Hydraulic shear stimulation hot dry rock reservoir permeability increasing simulation experiment system and experiment method
CN107476776B (en) Temporary plugging agent plugging performance experimental method for fracturing
CN113513295B (en) Method for improving uniform extension and support of multi-cluster cracks in section
Ramos et al. Sand production in vertical and horizontal wells in a friable sandstone formation, North Sea
CN110397436B (en) Oil reservoir plugging adjustment feasibility analysis method and system
CN111810108A (en) Dynamic adjusting system and method for back-flow oil nozzle after shale gas horizontal well pressure
Liu et al. Fracturing fluid leakoff behavior in natural fractures: effects of fluid rheology, natural fracture deformability, and hydraulic fracture propagation
CN108952677A (en) A kind of small scale fracturing test method being suitable for shale gas well
CN107288595B (en) Evaluation method for water injection utilization rate
CN114075974B (en) Temporary plugging steering fracturing method, device, control equipment and storage medium
CN110617040B (en) Periodic water injection method for fractured reservoir well group
CN113051698A (en) Method and system for determining closing time of intermittent gas well
Yinpeng et al. Fracturing Technologies with Profile Control and Water Shutoff for Medium and High Water-Cut Wells in Ultra-Low Permeability Reservoirs of Changqing Oilfield
CN115828639B (en) Profile control and profile control combined process scheme adjustment method
CN113338882B (en) Method and system for determining injection quantity of medicament for controlling bottom water of fracture-cavity type oil reservoir separator
CN114201927A (en) Calculation method for shale oil volume fracturing reconstruction space permeability
CN116629165B (en) Reservoir fracturing reconstruction area and non-reconstruction area parameter inversion method, system and equipment
CN112818513B (en) Evaluation method and evaluation model for permeability of fractured shale and construction method thereof

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