CN110162878B - Method for determining lower limit of crude oil filling pore throat of tight reservoir - Google Patents
Method for determining lower limit of crude oil filling pore throat of tight reservoir Download PDFInfo
- Publication number
- CN110162878B CN110162878B CN201910422497.9A CN201910422497A CN110162878B CN 110162878 B CN110162878 B CN 110162878B CN 201910422497 A CN201910422497 A CN 201910422497A CN 110162878 B CN110162878 B CN 110162878B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- filling
- pore throat
- crude oil
- lower limit
- 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
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 54
- 239000010779 crude oil Substances 0.000 title claims abstract description 46
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 47
- 239000003921 oil Substances 0.000 claims abstract description 36
- 238000005429 filling process Methods 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 19
- 238000004088 simulation Methods 0.000 claims description 18
- 238000009933 burial Methods 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 16
- 239000012530 fluid Substances 0.000 claims description 12
- 239000008398 formation water Substances 0.000 claims description 8
- 230000002706 hydrostatic effect Effects 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- 238000009736 wetting Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 abstract description 7
- 239000003208 petroleum Substances 0.000 abstract description 7
- 238000011160 research Methods 0.000 abstract description 5
- 238000009826 distribution Methods 0.000 abstract description 3
- 239000013589 supplement Substances 0.000 abstract description 3
- 230000000875 corresponding effect Effects 0.000 description 7
- 239000007789 gas Substances 0.000 description 7
- 238000011161 development Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008278 dynamic mechanism Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000010426 asphalt Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Strategic Management (AREA)
- Mining & Mineral Resources (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Geology (AREA)
- Entrepreneurship & Innovation (AREA)
- Agronomy & Crop Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Primary Health Care (AREA)
- Animal Husbandry (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Computer Hardware Design (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a calculation method in the field of petroleum geology research, in particular to a method for determining the lower limit of a crude oil filling pore throat of a tight reservoir, which comprises the following steps: (1) Based on the relation between the power and the resistance in the filling process, when the residual pressure P is equal to the capillary force Pc, the lower filling limit is reached; (2) And calculating and determining the lower limit of the filling pore throat of the crude oil according to a balance relation of power and resistance in the filling process. The invention deepens the understanding of the mechanism of filling the compact oil, supplements the geological theory of the compact oil reservoir formation, and provides an effective means for predicting the distribution range of the compact oil, thereby having great guiding value and wide application prospect in the current and future unconventional oil and gas exploration fields.
Description
Technical Field
The invention relates to a calculation method in the field of petroleum geology research, in particular to a method for determining a compact reservoir interface of an oil-gas-containing basin and a lower limit of a filling pore throat of crude oil in a reservoir by means of a basin numerical simulation method.
Background
Dense oil is an important exploration field of current and future oil and gas resources, and the accumulation mechanism of the dense oil, such as a filling mechanism, a seepage mechanism, an occurrence mechanism and an aggregation mechanism, is not clearly known. The lower limit of the filling pore throat of the dense oil and the change of the lower limit of the filling pore throat in the process of accumulation are the primary problems to be solved. The pore throat filling lower limit is the minimum physical boundary of the oil and gas to overcome migration resistance under specific geological conditions so as to fill the reservoir in the oil and gas reservoir forming process. Thus, different reservoir conditions have different pore throat filling lower limits that vary during the reservoir evolution even under specific reservoir conditions.
The research institute Zou Cai of petroleum exploration and development can add the thickness of the bound water film and the maximum diameter (4 nm) of the asphalt molecules to calculate that the pore throat lower limit of the compact oil reservoir is 54nm, and a calculation model is shown in figure 1, wherein B in figure 1 is a local section of A in figure 1, (a) is the critical pore throat radius of 54nm, (B) is the pore throat diameter of 200nm, and (c) is the average thickness of the bound water film of 50nm. The estimation method adopting the summation of the thickness of the bound water film and the diameter of the oil molecule is simple, convenient and quick, and can obtain an approximate lower limit of pore throat. However, this method is limited to static descriptions, departs to some extent from actual geological conditions and hydrocarbon reservoir conditions, and fails to account for changes in the lower fill throat limit during crude oil filling. The lower filling pore throat limit and the change of the filling pore throat limit in the crude oil filling process are determined by the balance condition of power and resistance, and are comprehensively influenced by a plurality of factors such as the intrinsic characteristics of a reservoir stratum, the characteristics of source rocks, the properties of the crude oil, the burial depth and the burial history.
According to the statistical relationship between the residual oil content and the average pore diameter, the medium petroleum exploration and development institute Cui Jingwei and the like think that the filling lower limit of the seven-segment compact oil of the extended group length of the Ordos basin is about 15nm. This example statistics demonstrates that the lower dense oil filling pore throat limit can be smaller, unlike the theoretical value that Zou Cai can yield, etc. However, the method is completely based on the statistical relationship between the aperture and the residual quantity, on one hand, the statistical data points are too few to reflect the actual geological condition, and on the other hand, the limit values of the oil-containing layer and the non-oil-containing layer are determined based on the mutation points on the statistical relationship and lack of geological basis.
The physical simulation experiment of dense oil reservoir formation such as China university of Petroleum (Beijing) Wu Kangjun and the like obtains that the lower limit of the filling pore throat of dense oil to form a large-scale oil reservoir is 0.12 mu m. The physical simulation experiment shows that the lower fill orifice throat limit is related to the source reservoir pressure differential. However, the experimental conditions of the physical simulation experiment are far from the actual geological conditions, the selected samples are not continuous, and a relationship formula with a dynamic connotation in a general sense is established on the basis of data and knowledge obtained by the physical simulation experiment, so that a more accurate lower filling limit of the pore throat can be deduced.
The relation between the water film thickness and the throat radius is established by Wang Weiming and the like of China university of Petroleum (east China) at different formation pressures, and the chart is shown in figure 2, wherein a straight line C represents that the water film thickness is equal to the pore throat diameter, an area A above the straight line C represents that the throat radius is smaller than the water film thickness, and an area B represents that the throat radius is larger than the water film thickness. And obtaining a filling porosity lower limit according to the relation between the water film thickness and the porosity. Analysis shows that when the radius of the throat is smaller than the thickness of the water film, the corresponding pore throat and the controlled micro pores are saturated by bound water, and only when the radius of the throat is larger than the thickness of the water film, the throat can become an effective filling channel for dense oil. The focus of this understanding is, however, the relationship of the water film thickness to the formation pressure, and whether the throat is fully occupied by a water film, the emphasis being on "whether there is residual space that can be filled with crude oil". In fact, there is such a residual space that the crude oil does not necessarily fill in.
Based on a filling mechanical equilibrium relationship, the research institute Zhang Hong for petroleum exploration and development obtains filling pore throats of compact oil source storage interfaces of an Ordorsi basin extension group, a Jurassic system in a Sichuan basin and a Bakken group in a West basin in the United states, wherein the lower limits of the filling pores near the storage interfaces are respectively 15.74nm, 29.06nm and 14.22nm, and the lower limits of the filling pores in a storage layer are respectively 39.45nm, 37.20nm and 52.32nm. The basic principle to which this method applies is also a charge kinetic equilibrium relationship, recognizing that the pore-throat charge lower limits at the source-reservoir interface and the reservoir interior are different, refining the knowledge of the charge mechanism. However, in this method, the conversion of both the source reservoir pressure differential and the residual pressure within the reservoir due to hydrocarbon production pressurization to maximum depth of burial is simplistic, especially in the presence of lift-and-subsidence (over time), the changes in source reservoir pressure differential and residual pressure are not simply directly correlated with depth. It is questionable that the application of formation fracture pressure without fluid pressure in the event of a formation fracture may be questionable because even if the formation fractures, the fluid pressure in the formation is present and not fully released. In addition, the selection of key parameters in the method can not obtain the change of the pore throat filling lower limit in the occlusion process, thereby causing the undersea understanding of the filling process.
Disclosure of Invention
Aiming at the defects of the prior art and the research, the invention provides a calculation method for determining the lower limit of the filling pore throat of crude oil in a reservoir interface and the reservoir based on a basin numerical simulation method and an oil-gas filling dynamic balance principle.
The invention provides a method for determining the lower limit of a crude oil filling pore throat of a tight reservoir, which comprises the following steps:
(1) Based on the relation between the power and the resistance in the filling process, when the residual pressure P is equal to the capillary force Pc, the lower filling limit is reached;
(2) And calculating and determining the lower limit of the filling pore throat of the crude oil according to a balance relation of power and resistance in the filling process.
Wherein, the balance relation of the power and the resistance in the filling process in the step (2) is as follows:
wherein, P is charging power and MPa; pc is filling resistance, namely capillary force, MPa; sigma is the gas-water interfacial tension, N/m; theta is the wetting angle, °; r is the critical pore throat radius of filling, m.
When the filling pore throat is positioned on the reservoir interface, the evolution process of the residual pressure of the reservoir interface is obtained by a basin simulation method and quantified as a function of time; when the filling pore throat is positioned in the reservoir, the residual pressure in the reservoir is calculated by a fluid state equation, and a formation temperature evolution history is obtained by using a basin simulation method, so that a quantitative expression of the residual pressure evolution in the reservoir can be obtained.
Wherein the relation of the residual pressure of the reservoir interface changing along with the time is as follows:
P S-R =f(t)
wherein, P S-R The residual pressure at the reservoir interface, MPa; t is time, ma.
Wherein, the lower limit of the crude oil filling pore throat of the compact reservoir interface is as follows:
wherein d is S-R Filling the reservoir interface with crude oil down the throatLimit, m, typically converted to nm; r is a radical of hydrogen S-R The crude oil filling pore throat radius, m, for the reservoir interface is typically converted to nm.
Wherein the relation of the residual pressure in the reservoir varying with the formation temperature is as follows:
wherein, P R The residual pressure inside the reservoir, MPa; p F Is the fluid pressure, MPa, inside the reservoir; p w Is the hydrostatic pressure, MPa, inside the reservoir; r = 8.3145J/(mol · K); t is the formation temperature, k; v m Is the molar volume of the crude oil, ml/mol; a =3.0346V w -5.43;B=0.1975A 2.0679 ;α=3.1725A 3.0456 ;β=1.539ω+2.071;V w Is the molar volume of formation water, V w =150mL/mol; ω =0.68; rho is the density of the formation water, kg/m3; g is the gravity acceleration, and the value is 9.8N/kg; h is the formation burial depth, m.
The residual pressure inside the reservoir is mainly related to the formation temperature, the molar volume and the formation burial depth, the change of the molar volume is mainly caused by the change of components in the oil production process, a quantitative expression of the evolution of the residual pressure inside the reservoir can be obtained by utilizing the history of the evolution of the formation temperature along with time and the history of the evolution of the formation burial depth along with time, and the relational expression of the fluid pressure inside the reservoir along with the change of the formation temperature is as follows:
P F =f(T)=f'(t)
wherein T is the formation temperature, k; t is time, ma;
the formation burial depth can be determined as a function of time by a basin simulation method, so the evolution process of the hydrostatic pressure is as follows:
P W =f”(t);
the residual pressure inside the reservoir as a function of time is:
P R =P F -P W =f'(t)-f”(t)。
wherein, the lower limit of the filling pore throat of the crude oil in the compact reservoir is as follows:
wherein d is R Filling the lower pore throat limit, m, for the crude oil inside the reservoir, typically converted to nm; r is R The pore throat radius, m, is filled with the crude oil inside the reservoir, typically converted to nm.
The invention has the beneficial effects that: the invention discloses a method for determining a lower limit of a crude oil filling pore throat of a tight reservoir, (1) the method indicates that the lower limit of the crude oil filling pore throat is changed in the reservoir forming process and is not a fixed value; (2) the crude oil filling process is illustrated as a residual pressure decay process, and the corresponding pore throat filling lower limit is gradually increased; (3) the change characteristics of the lower pore throat filling limit are explained from the aspect of the dynamic mechanism, and a calculation formula of the lower pore throat filling limit is established. The invention deepens the understanding of the mechanism of filling the compact oil, supplements the geological theory of the compact oil reservoir formation, and provides an effective means for predicting the distribution range of the compact oil, thereby having great guiding value and wide application prospect in the current and future unconventional oil and gas exploration fields.
Drawings
FIG. 1 is a graphical representation of a derivation of the lower dense pore throat fill threshold based on adsorbed water film thickness;
FIG. 2 is a graphical depiction of the relationship of water film thickness to throat radius under different formation pressure conditions;
FIG. 3 is a graphical representation of the residual pressure at the reservoir interface as a function of time for an implementation of the present invention;
FIG. 4 is a graphical representation of formation temperature as a function of time for an implementation of the present invention;
FIG. 5 is a graphical representation of the formation burial depth as a function of time for an implementation of the present invention, wherein FIGS. 5A and 5B are corresponding schematic and graphical representations;
fig. 6 is a graph of the change in residual pressure during inflation and the corresponding lower throat fill limit for an implementation of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings and examples.
The method for determining the lower limit of the filling pore throat of crude oil at the compact reservoir interface and in the reservoir comprises the following steps:
(1) Based on the relation between the power and the resistance in the filling process, when the residual pressure (P) is equal to the capillary force (Pc), the lower filling limit is reached;
(2) The evolution process of the residual pressure of the reservoir interface is obtained by a basin simulation method and quantified as a function of time;
the functional relationship between the residual pressure evolution of the reservoir interface obtained by the basin simulation method and the time is as follows:
P S-R =f(t);
wherein, P S-R The residual pressure at the reservoir interface, MPa; t is time, ma.
In the invention, the evolution process of the residual pressure in the reservoir is obtained by a fluid state equation;
wherein, P R The residual pressure inside the reservoir, MPa; r = 8.3145J/(mol · K); t is the formation temperature, k; v m Is the molar volume of the crude oil, ml/mol; a =3.0346V w -5.43;B=0.1975A 2.0679 ;α=3.1725A 3.0456 ;β=1.539ω+2.071;V w Is the molar volume of formation water, V w =150mL/mol; ω =0.68; rho is the density of the formation water, kg/m3; g is the gravity acceleration, and the value is 9.8N/kg; h is the formation burial depth, m.
The residual pressure in the reservoir is mainly related to the temperature of the stratum, the molar volume and the buried depth of the stratum, and the change of the molar volume is mainly caused by the change of components in the crude oil process, so the crude oil component change is not considered for the moment. Therefore, the evolution quantitative expression of the residual pressure in the reservoir can be obtained by utilizing the evolution history of the formation temperature along with time and the evolution history of the formation burial depth along with time.
As shown in fig. 4, the formation temperature may be determined by basin simulation method as the evolution process:
T=f'(t);
wherein T is the formation temperature, k; t is time, ma.
As shown in fig. 5, the stratum burial depth can be determined by the basin simulation method as follows:
h=f”(t);
the function relationship between the residual pressure in the reservoir and the time is as follows:
P R =f(T)-ρgh=f'(t)-f”(t);
Wherein, sigma is the gas-water interfacial tension, N/m; theta is the wetting angle, °; r is the critical pore throat radius, m;
the pore throat fill lower limit of the reservoir interface of one of the present invention is:
wherein d is S-R The lower limit of the reservoir interface crude oil filling pore throat, m, is usually converted to nm; r is S-R The crude oil filling pore throat radius, m, for the reservoir interface is typically converted to nm.
The pore throat fill lower limit inside the reservoir of one aspect of the invention is:
wherein, d is R Filling the lower pore throat limit, m, for the crude oil inside the reservoir, typically converted to nm; r is S-R The pore throat radius, m, is filled with the crude oil inside the reservoir, typically converted to nm.
The invention has the beneficial effects that: the lower throat filling limit is illustrated as varying during occlusion, and is not a fixed value; the crude oil filling process is illustrated as a residual pressure decay process, and the corresponding pore throat filling lower limit is gradually increased; the change characteristics of the lower pore throat filling limit are explained from the aspect of the dynamic mechanism, and a calculation formula of the lower pore throat filling limit is established. The invention deepens the understanding of the mechanism of filling the compact oil, supplements the geological theory of the compact oil reservoir formation, and provides an effective means for predicting the distribution range of the compact oil, thereby having great guiding value and wide application prospect in the current and future unconventional oil and gas exploration fields.
FIG. 3 shows a formation pressure evolution process obtained by a basin numerical simulation method; as shown in fig. 3, the evolution process of the residual pressure at the reservoir interface is fitted as a function of time:
P S-R =f(t)=-0.00008t 5 +0.004t 4 -0.101t 3 +0.951t 2 -5.299t+32.11;
wherein, P S-R The residual pressure at the reservoir interface, MPa; t is time, ma.
The capillary force calculation formula is as follows:
wherein, sigma is the gas-water interfacial tension, N/m; theta is the wetting angle, °; r is the critical pore throat radius, m;
according to the dynamic equilibrium relation at the pore throat filling lower limit, the pore throat filling lower limit of the reservoir interface can be obtained as follows:
wherein σ =0.172 × ρ oil -113.55;ρ oil =916kg/m 3 (ii) a θ =0 °; t is between 20Ma and 0 Ma.
As shown in fig. 6, the evolution process of the residual pressure at the reservoir interface between 20Ma and 0Ma and the corresponding variation process of the pore throat filling lower limit can be obtained.
And calculating according to a fluid state equation to obtain the residual pressure inside the reservoir as follows:
P R the residual pressure inside the reservoir, MPa; p F Is the fluid pressure, MPa, inside the reservoir; p w Hydrostatic pressure, MPa, inside the reservoir; r = 8.3145J/(mol · K); t is the formation temperature, k; v m Is the molar volume of the crude oil, ml/mol; a =3.0346V w -5.43;B=0.1975A 2.0679 ;α=3.1725A 3.0456 ;β=1.539ω+2.071;V w Is the molar volume of formation water, V w =150mL/mol; ω =0.68; rho is the density of the formation water, kg/m3; g is the gravity acceleration, and the value is 9.8N/kg; h is the formation burial depth, m.
The fluid pressure inside the reservoir may be expressed as a function of formation temperature:
P F =f(T)=f'(t)=-0.009t 2 -0.086t+59.9
wherein T is the formation temperature, k; t is time, ma.
The formation burial depth can be determined as a function of time by a basin simulation method, so the evolution process of the hydrostatic pressure is as follows:
P W =f”(t)=-0.003t 3 +0.094t 2 -1.841t+38.63
the residual pressure inside the reservoir as a function of time is:
P R =P F -P W =f'(t)-f”(t)=0.003t 3 -0.103t 2 +1.755t+21.27
according to the dynamic balance relation at the pore throat filling lower limit, the pore throat filling lower limit in the reservoir can be obtained as follows:
wherein σ =0.172 × ρ oil -113.55;ρ oil =916kg/m 3 (ii) a θ =0 °; t is between 20Ma and 0 Ma.
As shown in fig. 6, the evolution process of the residual pressure inside the reservoir between 20Ma and 0Ma and the corresponding variation process of the pore throat filling lower limit can be obtained.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only illustrative of the present invention and are not intended to limit the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (5)
1. A method for determining a lower limit of a crude oil filling pore throat of a tight reservoir, comprising the steps of:
(1) Based on the relation between the power and the resistance in the filling process, when the residual pressure P is equal to the capillary force Pc, the lower filling limit is reached;
(2) Calculating and determining the lower limit of the filling pore throat of the crude oil according to a balance relation of the power and the resistance in the filling process, wherein the balance relation of the power and the resistance in the filling process in the step (2) is as follows:
wherein, P is filling power, MPa; pc is filling resistance, namely capillary force, MPa; sigma is gas-water interfacial tension, N/m; theta is the wetting angle, °; r is the radius of the filling critical pore throat, m;
when the filling pore throat is positioned on the reservoir interface, the evolution process of the residual pressure of the reservoir interface is obtained by a basin simulation method and quantified as a function of time; when the filling pore throat is positioned in the reservoir, the residual pressure in the reservoir is calculated by a fluid state equation, and a basin simulation method is utilized to obtain a formation temperature evolution history so as to obtain a residual pressure evolution quantitative expression in the reservoir;
the relation of the residual pressure in the reservoir along with the change of the formation temperature is as follows:
wherein, P R The residual pressure inside the reservoir, MPa; p is F Is the fluid pressure, MPa, inside the reservoir; p w Is the hydrostatic pressure, MPa, inside the reservoir; r = 8.3145J/(mol · K); t is the formation temperature, k; v m Is the molar volume of the crude oil, ml/mol; a =3.0346V w -5.43;B=0.1975A 2.0679 ;α=3.1725A 3.0456 ;β=1.539ω+2.071;V w Is the molar volume of formation water, V w =150mL/mol; ω =0.68; rho is the density of the formation water, kg/m3; g is the gravity acceleration, and the value is 9.8N/kg; h is the formation burial depth, m.
2. The method of determining a lower crude oil charge pore throat limit for a tight reservoir of claim 1, wherein the residual pressure at the reservoir interface over time is related by:
P S-R =f(t)
wherein, P S-R Is the residual pressure at the reservoir interface, MPa; t is time, ma.
3. The method of determining a tight reservoir crude oil charge pore throat lower limit of claim 2, wherein the tight reservoir interface crude oil charge pore throat lower limit is:
wherein d is S-R The lower limit of the reservoir interface crude oil filling pore throat, m, is usually converted to nm; r is S-R The crude oil filling pore throat radius, m, for the reservoir interface is typically converted to nm.
4. The method for determining the lower limit of the crude oil filling pore throat of the tight reservoir as claimed in claim 1, wherein the residual pressure inside the reservoir is mainly related to the formation temperature, the molar volume and the formation burial depth, the change of the molar volume is mainly caused by the change of components in the oil production process, a quantitative expression of the evolution of the residual pressure inside the reservoir can be obtained by using the history of the evolution of the formation temperature with time and the history of the evolution of the formation burial depth with time, and the relation of the fluid pressure inside the reservoir with the change of the formation temperature is as follows:
P F =f(T)=f'(t)
wherein T is the formation temperature, k; t is time, ma;
the formation burial depth can be determined as a function of time by a basin simulation method, so the evolution process of the hydrostatic pressure is as follows:
P W =f”(t);
the residual pressure inside the reservoir as a function of time is:
P R =P F -P W =f'(t)-f”(t)。
5. the method of determining a lower limit of a tight reservoir crude oil filling pore throat according to claim 4, wherein the lower limit of the tight reservoir crude oil filling pore throat is:
wherein d is R Filling the lower pore throat limit, m, for the crude oil inside the reservoir, typically converted to nm; r is R The pore throat radius, m, is filled with the crude oil inside the reservoir, typically converted to nm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910422497.9A CN110162878B (en) | 2019-05-21 | 2019-05-21 | Method for determining lower limit of crude oil filling pore throat of tight reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910422497.9A CN110162878B (en) | 2019-05-21 | 2019-05-21 | Method for determining lower limit of crude oil filling pore throat of tight reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110162878A CN110162878A (en) | 2019-08-23 |
CN110162878B true CN110162878B (en) | 2023-04-18 |
Family
ID=67631625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910422497.9A Active CN110162878B (en) | 2019-05-21 | 2019-05-21 | Method for determining lower limit of crude oil filling pore throat of tight reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110162878B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111075443B (en) * | 2019-12-31 | 2021-08-27 | 成都理工大学 | Natural gas filling semi-quantitative measuring system and method suitable for low-abundance gas reservoir |
CN113759099B (en) * | 2021-09-07 | 2023-07-21 | 重庆科技学院 | Quantitative evaluation method for oil-gas filling capacity of source-storage side-connected oil-gas reservoir |
CN117826247B (en) * | 2024-01-06 | 2024-07-05 | 中国地质科学院地质力学研究所 | Carbonate rock oil-gas reservoir geological process reconstruction method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675236A (en) * | 2013-12-19 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and device for measuring critical tight oil filling pore throat radius threshold |
CN107228934A (en) * | 2017-06-27 | 2017-10-03 | 中国石油大学(华东) | The determination method of tight sand hydrocarbon charge pore throat radius lower limit |
-
2019
- 2019-05-21 CN CN201910422497.9A patent/CN110162878B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675236A (en) * | 2013-12-19 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and device for measuring critical tight oil filling pore throat radius threshold |
CN107228934A (en) * | 2017-06-27 | 2017-10-03 | 中国石油大学(华东) | The determination method of tight sand hydrocarbon charge pore throat radius lower limit |
Non-Patent Citations (1)
Title |
---|
酒泉盆地青西凹陷下沟组混积层系致密油成藏机理与富集影响因素;郭迎春等;《石油与天然气地质》;20180802;第39卷(第04期);第1节第1段至第5节最后一段 * |
Also Published As
Publication number | Publication date |
---|---|
CN110162878A (en) | 2019-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110162878B (en) | Method for determining lower limit of crude oil filling pore throat of tight reservoir | |
Daigle | Relative permeability to water or gas in the presence of hydrates in porous media from critical path analysis | |
Lei et al. | Laboratory strategies for hydrate formation in fine‐grained sediments | |
Konno et al. | Hydraulic fracturing in methane-hydrate-bearing sand | |
Kneafsey et al. | Permeability of laboratory-formed methane-hydrate-bearing sand: measurements and observations using X-ray computed tomography | |
Huang et al. | Evaluation on the gas production potential of different lithological hydrate accumulations in marine environment | |
Hyodo et al. | Mechanical behavior of gas‐saturated methane hydrate‐bearing sediments | |
Delli et al. | Prediction performance of permeability models in gas-hydrate-bearing sands | |
Su et al. | Evaluation on gas production potential from laminar hydrate deposits in Shenhu Area of South China Sea through depressurization using vertical wells | |
Ning et al. | Invasion of drilling mud into gas-hydrate-bearing sediments. Part I: effect of drilling mud properties | |
Mingjie et al. | Coupling relationship between sandstone reservoir densification and hydrocarbon accumulation: A case from the Yanchang Formation of the Xifeng and Ansai areas, Ordos Basin | |
Su et al. | Sedimentary evolution of the central canyon system in Q iongdongnan Basin, northern South China Sea | |
Ruppel | Tapping methane hydrates for unconventional natural gas | |
CN110646332B (en) | Method for determining movable water saturation of gas-water interbed gas reservoir under high-temperature and high-pressure conditions | |
Zhao et al. | Gas permeability characteristics of marine sediments with and without methane hydrates in a core holder | |
Yu et al. | Distribution and controlling factors of tight sandstone oil in Fuyu oil layers of Da'an area, Songliao Basin, NE China | |
Jian et al. | Application of mercury injection and rate-controlled mercury penetration in quantitative characterization of microscopic pore structure of tight reservoirs: A case study of the Chang7 reservoir in Huachi-Heshui area, the Ordos Basin | |
Zang et al. | Gas hydrate formation in fine sand | |
CN104200105A (en) | Method for determining filling property lower limit of tight sandstone gas | |
Shi et al. | Research into polymer injection timing for Bohai heavy oil reservoirs | |
CN111577264A (en) | Method and device for predicting capacity of fractured-pore oil reservoir horizontal well | |
Cheng et al. | Gas–liquid–solid migration characteristics of gas hydrate sediments in depressurization combined with thermal stimulation dissociation | |
CN110242261B (en) | Method and system for predicting microcosmic seepage rule of oil, gas and water in vertical gas injection pore | |
CN111400854B (en) | Gas injection breakthrough time prediction method for gas injection miscible oil displacement reservoir | |
Minagawa et al. | Water permeability of porous media containing methane hydrate as controlled by the methane-hydrate growth process |
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 |