CN110909466A - Method for improving calculation precision of differential pore sorting reservoir average capillary pressure curve - Google Patents
Method for improving calculation precision of differential pore sorting reservoir average capillary pressure curve Download PDFInfo
- Publication number
- CN110909466A CN110909466A CN201911147846.7A CN201911147846A CN110909466A CN 110909466 A CN110909466 A CN 110909466A CN 201911147846 A CN201911147846 A CN 201911147846A CN 110909466 A CN110909466 A CN 110909466A
- Authority
- CN
- China
- Prior art keywords
- formula
- regression
- capillary pressure
- reservoir
- pressure 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000011148 porous material Substances 0.000 title abstract description 36
- 238000004364 calculation method Methods 0.000 title abstract description 19
- 238000009736 wetting Methods 0.000 claims abstract description 21
- 238000002474 experimental method Methods 0.000 claims abstract description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 42
- 238000012545 processing Methods 0.000 claims description 7
- 241000764238 Isis Species 0.000 claims description 4
- 230000035699 permeability Effects 0.000 abstract description 12
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 16
- 229910052753 mercury Inorganic materials 0.000 description 16
- 235000019738 Limestone Nutrition 0.000 description 14
- 239000006028 limestone Substances 0.000 description 14
- 230000008569 process Effects 0.000 description 7
- 239000011435 rock Substances 0.000 description 7
- 238000009795 derivation Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 235000020681 well water Nutrition 0.000 description 3
- 239000002349 well water Substances 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 241000209117 Panicum Species 0.000 description 1
- 235000006443 Panicum miliaceum subsp. miliaceum Nutrition 0.000 description 1
- 235000009037 Panicum miliaceum subsp. ruderale Nutrition 0.000 description 1
- 206010048669 Terminal state Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- AHLBNYSZXLDEJQ-FWEHEUNISA-N orlistat Chemical compound CCCCCCCCCCC[C@H](OC(=O)[C@H](CC(C)C)NC=O)C[C@@H]1OC(=O)[C@H]1CCCCCC AHLBNYSZXLDEJQ-FWEHEUNISA-N 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- 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
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a regression method for effectively improving the calculation accuracy of an average capillary pressure curve of a differential pore sorting reservoir, which comprises the following steps: selecting a plurality of core samples on the same reservoir stratum to carry out capillary pressure experiments, recording corresponding numerical values of capillary pressure and wetting phase saturation, carrying out J function dimensionless treatment on all capillary pressure curves, recording corresponding numerical values of J values and wetting phase saturation, and respectively carrying out multiple nonlinear regression on data by using fractal J function regression formulas 1 and 2 to obtain two sets of regression parameters and a determination coefficient R2Selecting the regression result to be within the data point range, and determining the coefficient R2The formula closer to 1 is used as a basic formula of the average capillary pressure curve of the reservoir, and the average capillary pressure and the saturation of the wetting phase are calculated according to the J function definition formula by utilizing the regression parameters, the interfacial tension, the wetting angle, the average porosity and the average absolute permeability of the reservoirI.e. the mean capillary pressure curve.
Description
Technical Field
The invention relates to the field of oilfield development, and discloses an application method for effectively improving the calculation accuracy of an average capillary pressure curve of a differential pore sorting reservoir.
Background
The method mainly comprises a J function method for calculating an average capillary pressure curve representing the averaging property of pore structures of the same type of reservoirs in the field of oil field development at present, wherein the method introduces a capillary pressure curve for testing a single core sample of a reservoir into the capillary pressure curveAnd performing dimensionless operation, performing power function regression on the J function values and the water saturation of a plurality of cores of the reservoir, and obtaining an average capillary pressure curve of the reservoir together with the porosity and the absolute permeability averaged by the reservoir. In theoretical and practical application, the method mainly has the following two problems: (1) measured pore sorting differenceThe J function curve of the reservoir is not in an empirical power function form, namely the regression precision of an empirical formula is low, and the applicability is limited; (2) the empirical formula of the J function lacks a theoretical basis, and regression parameters have no physical significance.
The average capillary pressure curve of the reservoir has important significance in representing the averaging property of complex pore structures containing microcracks, karst caves and the like, is used as a middle parameter for calculating the initial water saturation of an oil reservoir to influence the reserve calculation result and also used as a direct input parameter for predicting numerical simulation development indexes to influence the yield prediction result, and for the reservoir with poor pore sorting, the regression precision of the current empirical power function on the dimensionless capillary pressure curve is low, so that research needs to be carried out aiming at the core of the technology, namely the dimensionless J function regression method, so as to improve the regression precision of the reservoir and provide theoretical explanation.
Disclosure of Invention
In view of the above problems, the present invention provides a regression analysis formula for effectively improving the regression accuracy of the J function of the poor reservoir.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method of calculating an average capillary pressure curve for a reservoir, comprising the steps of:
1) selecting a plurality of core samples on the same reservoir, performing a group of capillary pressure experiments on each core sample, and recording capillary pressure and wetting phase saturation data;
2) for each core sample, using the data recorded in step 1, respectively using the formulaCarrying out dimensionless J function processing, and establishing a numerical value corresponding relation between the J function and the water saturation;
3) respectively using the data recorded in the step 2 to respectively use the formula for all core samples of the same reservoir Performing regression solution, wherein the parameter variation range limits A2>0,C2>0,B2Is greater than 0, and regression parameter A is obtained2、B2、C2And a first determination coefficient
Respectively using the data recorded in the step 2 to respectively use the formula for all core samples of the same reservoir Performing regression solution, wherein the parameter variation range limits A4>0,C4>0,0<B4< 1, obtaining a regression parameter A4、B4、C4And a second determination coefficient
4) And (3) plotting the regression result and the data recorded in the step 2, selecting a formula with the regression result positioned in the range of the water saturation and the J function in the step 2, determining a formula with a coefficient closer to 1 as a basic formula of the average capillary pressure curve, and taking the regression parameters of all core samples of the same type of reservoir in the step 3) corresponding to the basic formula as the regression parameters corresponding to the average capillary pressure curve of the reservoir.
Wherein, the unit conversion coefficient a is different when processing because the capillary pressure units obtained by different testing methods are different.
Wherein the first and second determining coefficients are of the same type, i.e., the first and second determining coefficients are both of the same typeWherein SSR (regression sum of squares) is regression sum of squares, SSE (error sum of squares) is residual sum of squares, and SST (total sum of squares) is total sum of squares.
The selection method for selecting the optimal formula comprises the following steps:
3-1) taking the J function and water saturation relation data points of a plurality of dimensionless core samples of the same reservoir as true values, calculating the J function regression value corresponding to the water saturation by using the formula (2) and the parameters obtained by regression of the true values and the formula (2) as regression values 1, and recording the determination coefficient 1 as the regression value 1Calculating a J function regression value corresponding to the water saturation by using the formula (3) and a parameter obtained by the regression of the true value as the formula (3), taking the J function regression value as a regression value 2, and recording a determination coefficient 2 asSelecting a corresponding formula with a regression value within the range of the true value as an optimal formula;
3-2) for the case that all 2 formulas are within the range of true values, the magnitude of the coefficient is determined by comparison:
if it isIs greater than or equal toThe basic formula of the average capillary pressure curve of the reservoir to be detected is a formula (2);
if it isIs greater than or equal toThe basic formula of the average capillary pressure curve of the reservoir to be detected is formula (3).
In the step 3-2), when the basic formula of the average capillary pressure curve of the reservoir to be tested is the formula (2), the average capillary pressure curve of the reservoir to be tested isWherein A is2>0,C2>0,B2>0;
When the basic formula of the average capillary pressure curve of the reservoir to be tested is formula (3), the average capillary pressure curve of the reservoir to be tested isWherein A is4>0,C4>0,0<B4<1。
Wherein the number of the cores subjected to the capillary pressure experiment in the step 3-1) is more than or equal to 1.
Wherein, each core sample is subjected to 1 group of capillary pressure experiments; and recording the corresponding relation between the water saturation and the capillary pressure value, the wetting angle, the two-phase interfacial tension, and the porosity and the absolute permeability of the rock core in the capillary pressure experiment.
Wherein the regression parameter A4、A2、B4、B2Related to fractal dimension, inlet pressure, maximum capillary pressure, C4And C2Related to the fractal dimension.
Wherein, A is2、B2、C2Can be represented by formulaWherein A is2>0,C2>0,B2Is greater than 0; orWherein A is2>0,C2>0,B2More than 1, obtained by multivariate nonlinear regression, saturation of wetting phase SwNon-wetting phase saturation SnwHas a relationship of Sw+Snw=1;
A is described4、B4、C4Can be represented by formulaWherein A is4>0,C4>0,0<B4< 1 orWherein A is4>0,C4>0,B4Greater than 0, obtained by multivariate nonlinear regression, wet phase saturation SwNon-wetting phase saturation SnwHas a relationship of Sw+Snw=1。
The invention deduces the fractal expression of J function regression formula based on fractal geometry theory(applicable conditions D < 3 and rmin→0)、J=A(Sw+B)C(applicable conditions D < 3 and rmin→ epsilon > 0) and J ═ A (1-BS)w)C(applicable conditions 3 < D < 5 and rminGreater than 0), establishing the relation between the regression parameter of the original empirical power-law formula and the fractal dimension theoretically, and indicating that the formula is better for separating pores, the fractal dimension D is less than 3, and the minimum pore radius rmin→ 0, on the other hand, the function form of the original regression equation is generalized, the separation of pores is poor, the fractal dimension D is less than 3, rmin→ epsilon > 0 and 3 < D < 5, rminFor reservoirs larger than 0, two wide-fit fractal regression methods J ═ A (S) are providedw+B)CAnd J ═ A [1-BS ]w]CThrough the verification of actual oil field data, the two methods can effectively improve the J function regression precision (determining coefficient R) of the reservoir2Closer to 1), the initial water saturation obtained by the method is closer to the interpretation result of the saturation logging curve, thereby indirectly proving that the average capillary pressure curve of the reservoir obtained by the method has higher precision and is closer to the underground practical situation, and laterally proving the more universal applicability of the invention.
Due to the adoption of the technical scheme, the invention has the following advantages: 1. the method is characterized in that a relation is established between regression parameters and fractal dimension, maximum capillary pressure and minimum capillary pressure based on mathematical derivation of a fractal theory, the physical significance is clear, a solid theoretical basis is achieved, and a theoretical basis is provided for a general empirical formula (power function) at present; 2. the saturation does not need to be normalized, so that the data processing link is simplified; 3. the method has the advantages that the J function regression precision, the initial water saturation calculation precision and the capillary pressure curve calculation precision of the reservoir with poor pore sorting are effectively improved, and the application range is wider.
Drawings
FIG. 1 is a graph of regression of J-function for a particular actual limestone reservoir type as compared to the current method in example 1;
FIG. 2 is a graph of the regression of the J function of the present invention versus the current method for a practical limestone reservoir type II reservoir in example 1;
FIG. 3 is a graph of regression of J function for a practical limestone reservoir of example 1 in accordance with the present invention and the present method;
FIG. 4 is a graph of regression of J function for the present invention versus the current method for four types of reservoirs of an actual limestone reservoir in example 1;
FIG. 5 is a comparison of the average capillary pressure curve calculation results of a practical limestone reservoir of example 1, compared with the average capillary pressure curve calculation results of the prior art;
FIG. 6 is a comparison of the average capillary pressure curve calculation results of a practical limestone reservoir of example 1 in accordance with the present invention and the prior art;
FIG. 7 is a comparison of the average capillary pressure curve calculation results of a practical limestone reservoir of example 1 in accordance with the present invention and the prior art;
FIG. 8 is a comparison of the average capillary pressure curve calculation results of the present invention compared to the prior art for four types of actual limestone reservoirs in example 1;
FIG. 9 is the initial water saturation field for one to four reservoir types of an actual limestone reservoir determined using the present method in example 1;
FIG. 10 is the initial water saturation field for one to four reservoir types of an actual limestone reservoir determined using the method of the present invention in example 1;
FIG. 11 is a histogram statistics of the initial water saturation field of one to four reservoir types of an actual limestone reservoir determined using the present method in example 1;
FIG. 12 is a histogram statistics of the initial water saturation field of one to four types of reservoir of an actual limestone reservoir determined using the method of the present invention in example 1;
FIG. 13 is a comparison of the single well initial water saturation and logging interpretation saturation for four reservoir types of actual limestone reservoir, as determined by the present invention and the present method; from left to right, the five figures are respectively: BU-28 well, BU-29 well, BU-30 well, BU-31 well and BU-41 well, wherein the black line is a water saturation logging curve, the green line is a saturation curve calculated by the existing method, and the red line is a saturation curve calculated by the method.
Detailed Description
The invention is described in detail below with reference to the following formula derivations, the figures and the examples.
The expression of the dimensionless J function of the capillary pressure isIts existing regression formulaBroad-fit fractal regression formula 1:and broad-fit fractal regression formula 2:and the applicable conditions of each formula, and the derivation process of the parameter variation range is as follows.
According to the fractal geometry theory, if the pore size distribution of the reservoir pores conforms to the fractal structure, the fractal structure can be filled with a certain self-similar basic unit body, and according to a non-equal-diameter equal-length straight capillary bundle seepage model, the unit body can be a straight capillary with the radius and the length of r. In the mercury pressing process, in a certain mercury saturation increasing interval, the power function relationship exists between the number N (> r) of capillaries filled with mercury and the radius r of the cluster of capillaries just filled with mercury in the interval:
in the formula (I), the compound is shown in the specification,rmaxthe maximum capillary radius, P (r) is the pore size distribution density function, D is the fractal dimension, and a is the proportionality constant.
The expression of the pore size distribution density function P (r) can be obtained by deriving r according to the formula (1):
P(r)=-Da·r-D-1formula (2)
Substituting the formula (2) into the integral of the formula to obtain the accumulated pore volume of the interval with the capillary radius less than r and without being filled with mercury
In the formula (I), the compound is shown in the specification,and obtaining a total pore volume expression of the core in the same way:
wherein V (< r) represents the pore volume occupied by the wetting phase and V (r) represents the total pore volume.
The wetting phase air saturation in the interval is easy to obtain:
the relationship between capillary pressure and pore radius is:
where σ is the mercury-air interfacial tension, typically 480mN/m, and θ is the mercury-air wetting angle, typically 140 °.
Formula (6) is substituted for formula (5) to give:
if using SHg=1-SAirReplacement ofWhen the mercury saturation is satisfied, the formula (7) becomes the formula (8):
in the formulae (7) and (8), PmaxMaximum capillary pressure, P, corresponding to minimum capillary radiuseThe displacement pressure is the minimum capillary pressure corresponding to the maximum capillary radius.
The fractal dimension D of the medium and high porosity permeability reservoir stratum with the separated pores is less than 3, rmin<<rmaxThe pore radius is mainly concentrated above rminIf r is in the range of mercury intrusion to residual air saturationmin→ 0, known by the nature of power functionsCorresponds to Pmax→+∞,The formulas (7) and (8) can be simplified into the formulas (9) and (10):
the fractal function of the mercury intrusion into the heredity ancestor (9) or (10) is not normalized in saturation compared with the Brooks-Corey function, and because the fractal dimension D is less than 3, the point with the normalized saturation of 0 exceeds the definition domain when the power exponent of the power function is negative, the fractal function of the mercury intrusion into the heredity ancestor (9 or 10) is more reasonable.
Poor pore sorting reservoir with large fractal dimension, rmin<rmaxThe concentration of the small pore radius becomes high, rminNot to be ignored, two ranges of fractal dimensions are discussed:
(1) d < 3, poor separation, mercury intrusion to residual air saturation, ① if rmin→ 0 for Pmax→ infinity, also by power function natureStill obtain the formula (9), the formula (10); ② if rmin→ ε > 0 corresponds to PmaxThe maximum value, but not infinity, is derived from equation (7):
from formula (8):
(2) d is more than 3 and less than 5, sorting is extremely poor, the concentration degree of small pore radius is further increased, and the pore radius r corresponding to the saturation degree of residual airmin> 0, by power function propertiesCorresponds to Pmax>0,From formula (7):
from formula (8):
the formula (11) or the formula (12) is a mercury intrusion fractal function regression formula with D < 3, rminThe form in the case of → epsilon > 0, the form in which the formula (13) or the formula (14) is a fractal function regression of mercury intrusion in the case of 3 < D < 5, pore sorting is extremely poor, and the saturation is likewise not normalized as compared with the Li function.
And simplifying the seepage model with equal-diameter capillary bundles into rock pore space in a seepage state. From the law of panicum pratense, pressure differenceUnder the action of delta p, mercury with viscosity of mu flows through the radius of riThe flow of a single bundle of bristles of length L is:
n (r) root capillary flow rate is added to obtain the total flow rate of the porous medium:
the single capillary volume is expressed in terms of mercury saturation as:
in the formula (18), A is the core sectional area, L is the core length, and SHgiIs the mercury saturation of a single capillary tube. Substituting formula (18) for formula (17) yields:
alli darcy's law:
obtaining an expression for absolute permeability, and considering the continuous distribution of pores, yields:
the J function is defined as dimensionless capillary force and eliminates the effects of porosity and permeability:
for a conventional medium and high porosity and pore-sorted reservoir, formula (6), formula (10) and formula (21) are substituted into formula (22), and the formula can be obtained:
in practical applications, the equation (23) is simplified in form as follows, without being tied to the intrinsic correlation of each parameter:
the formula (24) is a J function power form empirical regression formula widely applied at present, and the application conditions of the formula (24) are revealed to be that the fractal dimension D is less than 3 and r is under the residual air saturation through derivation based on the fractal theorymin→ 0, the formula relates empirical parameter A0、B0Has a relationship with fractal dimension, wherein A0>0,B0>0。
When formula (6), formula (12) or formula (21) is substituted for formula (22), it is possible to obtain:
if formula (11) is substituted instead of formula (12), then:
equation 25 (equation 26) is a generalization of the empirical regression equation for the power-law type of J function, formally simplified as:
the method is applicable under the conditions that the fractal dimension D is less than 3 and r is under the saturation of residual airmin→ ε > 0, if rmin→ 0, corresponding to Pmax→+∞,In formula 25In formula 26The two equations are returned to the power law equation 24. Regression parameter A1、A2、B1、B2Related to fractal dimension, inlet pressure, maximum capillary pressure, C1、C2Only related to the fractal dimension. In the derivation process of the formula, r caused by the concentration of the pore diameter of a poor sorting reservoir in the direction of the small pore is mainly consideredminGreater than 0, but because the fractal dimension D is less than 3, the rock pore sorting cannot be too poor, rmin→ ε, so the application conditions of this formula are D < 3 and r at residual air saturationmin→ epsilon > 0, wherein A1>0,A2>0,C1>0,C2>0,B1>1,B2>0。
For the worse low-porosity and low-permeability reservoir, formula (6), formula (14) and formula (21) are substituted into formula (22):
if formula (11) is substituted instead of formula (12), then:
equations (29) and (30) may be simplified in form to:
formula (31) or formula (32), i.e., J-function fractal regression, at 3 < D < 5 and residual air saturation rminFunctional form in the case > 0, corresponding to the physical case where the worse-sorted reservoir pore size is further concentrated towards the small pore direction, resulting in rminAnd the fractal dimension 3 < D < 5 is the applicable condition of the formula. Regression parameter A3、A4、B3、B4Related to fractal dimension, inlet pressure, maximum capillary pressure, C3、C4Only related to the fractal dimension. Wherein A is3>0,A4>0,C3>0,C4>0,B3>0,0<B4<1。
It should be noted that the wetting angle and interfacial tension related to the fluid properties are eliminated in the derivation process, i.e., J-function fractal regression equations (27), (28), equations (31) and (32) are not only applicable to mercury intrusion process, but also applicable to any two-phase fluid displacement process, i.e., SAirCorresponding to the saturation of the wetting phase Sw,SHgCorresponding to the non-wetting phase saturation SnwThe relationship between the two is Sw+Snw=1。
In summary, the following steps: in order to improve the calculation precision of the average capillary pressure curve of the differentially selected reservoir, a formula (28) And type (32)Respectively in J function wide-fit fractal expression form, and the applicable conditions are D < 3 and rmin> 0 and 3 < D < 5, rmin>0。
1) First use formula (22)The method comprises the steps of performing J function dimensionless processing on actually measured capillary pressure curves (function relation between capillary pressure and water saturation) of a plurality of cores of the same reservoir, establishing numerical value corresponding relation between the J function and the water saturation, and paying attention to that unit conversion coefficients a are different during processing due to different capillary pressure units obtained by different testing methods, and the detailed table is shown in a table 1, and common interfacial tension and wetting angle data are shown in a table 2.
TABLE 1 tubing pressure curve J function dimensionless unit conversion coefficient table obtained by different testing methods
TABLE 2 data sheet for interfacial tension and wetting angle in general
2) By means of existing regression modelsRegression of J function and water saturation data of the same type of reservoir to obtain regression result and determination coefficient R2(ii) a Using a wide-fitting fractal function 1And a broad fractal function 2Regression is carried out on J functions of the same type of reservoir, the regression result is selected to be located in the J function and water saturation data point range of the reservoir, and a coefficient R is determined2The fractal function closer to 1 serves as an optimal regression method.
3) The parameter result can be obtained by solving a multivariate nonlinear regression program, namely a wide-fit fractal function 1 Has a parameter variation range of A2>0,C2>0,B2Is greater than 0; wide-fit fractal function 2 Has a parameter variation range of A4>0,C4>0,0<B4<1。
4) According to the average porosity and absolute permeability of the same type of reservoir, according to S in the 4 th stepwAnd PcThe average capillary pressure curve of the reservoir from initial to terminal states is solved for subsequent simulation of the development scheme.
In order to prove that the average capillary pressure curve of the reservoir obtained by the regression method of the invention is more consistent with the underground actual condition, the initial water saturation curve of the reservoir and the water saturation curve of the single well logging obtained by calculation under the control of the regression method of the invention and the existing method can be compared.
5) Analytical relationship between capillary pressure and J function in equation (22)Establishing an analytical relational expression of capillary pressure and water saturation, and for the same type of reservoir:
Since the calculator function in the mainstream modeling digital-analog integration software Petrel (version 2018.2) does not have a function with a power exponent of any power, three equations can be equivalently rewritten with an exponent and a logarithmic function with a base 10 as:
6) and establishing the attribute h (m) of each grid distance free water interface height difference in the model. According to the relation P of the initial capillary pressure and the height differencecAnd (4) solving the initial capillary pressure value of each grid of the model as delta rho gh.
7) According to S established in the step 4wAnd PcSubstituting the parameter value regressed in the step into the same type of reservoir:
③ if fractal function 2 is used, i.e.Then substitute for A4、B4、C4. Substituting the porosity and permeability of each grid to obtain the initial water saturation field S of the modelwiIt can be demonstrated by way of example that the single well water saturation curve obtained by the present invention will be closer to the single well logging saturation curve than that obtained by the prior art method, i.e. that the regression accuracy (determination coefficient R) is higher by using the present invention2The J function fractal regression equation which is closer to 1) is adopted, the solved average capillary pressure curve is more in line with the actual pore structure of the underground reservoir, and the oil field production dynamics and the recovery ratio simulated by the average capillary pressure curve are more in line with the underground reality.
Example 1
Taking a limestone reservoir of a middle east B oil field as an example, calculating the average capillary pressure curve of the reservoir according to the method
1. The B field reservoirs were classified into four categories, one category, two category, three category and four category according to the method for reservoir classification in "description of carbonate rock reservoir based on rock type" (pit school, liberty, authored by zhang jun, china oil university press, eastern 2016, 6 months, 1 st edition). The method comprises the following steps that a first type reservoir comprises 4 core samples, a second type reservoir comprises 19 core samples, a third type reservoir comprises 16 core samples, a fourth type reservoir comprises 4 core samples and 43 core samples, a group of capillary pressure experiments are carried out on each core sample, and capillary pressure and wetting phase saturation data are recorded;
2. for 43 capillary pressure curves corresponding to 43 rock samples respectively, using the formula (22)Dimensionless J function processing is carried out, the numerical value corresponding relation between the J function and the water saturation is established, and relevant parameters are shown in tables 1 and 2. And respectively aggregating the J function and the Sw data according to 4 types of reservoirs to obtain 4 groups of J function and Sw data point sets to be regressed, wherein the data points are shown in the data scatter points of the figures 1-4.
3. By means of existing regression modelsSequentially regressing the relation of J functions Sw of the same type of reservoir, and determining the coefficient R2The reservoir type 0.9482, the reservoir type 0.7551, the reservoir type 0.7101 and the reservoir type 0.7693 are respectively two fractal regression forms by using multivariate nonlinear regression software, wherein the fractal regression forms are shown in a formula (28)And type (32) Regression solution, in which the parameter variation limits A2>0,C2>0,0<B2<+∞,A4>0,C4>0,0<B4< 1, obtaining a regression parameter A2、B2、C2、A4、B4、C4And determining the coefficients
4. The regression curve is projected on the data scatter points of the graphs 1 to 4, the regression result of the formula (32) is obviously not in the range of data points for one type of reservoir, the regression results of the formula (28) and the formula (32) for the second, third and fourth types of reservoirs are in the range of data points, and coefficients are determined by comparisonAndlast-class reservoirThe fractal form (28) is selected for the three-type reservoir and the four-type reservoir, the fractal form (32) is selected for the two-type reservoir, and after the fractal regression formula is adopted, the regression results (solid lines in figures 1-4) are all located in the range of actually measured data points, the regression precision is improved compared with the existing regression formula, and the coefficient R is determined2In aspects, one type of reservoir was raised to 0.9978, two type of reservoir was raised to 0.9714, three type of reservoir was raised to 0.9978, and four type of reservoir was raised to 0.8129, as detailed in table 3.
TABLE 3 comparison of regression results for two fractal J-functions
5. According to the average porosity and absolute permeability of various reservoirs (see table 4 in detail), substituting the analytical relationship between capillary pressure and J function in formula (22)The average capillary pressure curve (the relation between capillary pressure and water saturation) in the whole process of oil field development can be solved and directly input parameters as a model for numerical simulation development index prediction, wherein four types of reservoirs respectively correspond to the parameters in the table 3, and the comparison result of the method with the prior method is shown in the graph 5-8. The results of the oil quantity, the water quantity, the pressure and the like of the single well obtained by inputting the calculation result of the improved formula into a numerical model for simulation are more consistent with the underground actual conditions. To demonstrate this conclusion, the following is a comparison of the single well water saturation calculation under regression control of the J-function of the present invention with the single well water saturation curve.
TABLE 4 mean porosity and Absolute Permeability of middle east B oilfield reservoir
Reservoir classification | Porosity, decimal fraction | Absolute permeability, mD |
Reservoir stratum | 0.141 | 228.0 |
Reservoir of the second kind | 0.159 | 12.5 |
Three types of reservoirs | 0.165 | 3.3 |
Four types of reservoirs | 0.161 | 1.1 |
6. Based on the original regression formulaThe established relationship between capillary pressure and water saturation is as follows:
the relationship between capillary pressure and water saturation established based on the improved fractal regression formula is as follows:
7. establishing the attribute h (m) of the height difference of each grid distance free water interface in a geological model, and according to the relation P between the initial capillary pressure and the height differencecThe initial capillary pressure of each grid is worked out as delta rho gh, the initial water saturation fields of the four types of reservoirs are worked out by four improved formulas in the last step, the initial water saturation fields are compared with the initial saturation fields worked out by the original formulas and are shown in figures 9-10, the statistical histograms are shown in figures 11-12, and as is obvious from figures 11-12, the water saturation peak calculated by the original method is all 0.2, while the method is about 0.35. This shows that the original water saturation field found by the present invention directly affects the model reserve calculation.
8. The comparison of the initial saturation curve of the single well and the saturation of the logging calculated by the improved formula and the original formula is shown in fig. 13, and it can be seen from the trend that compared with the calculation result of the existing method (solid line and hollow square in fig. 13), the water saturation (solid line and solid circle in fig. 13) under the control of the improved formula of the invention is obviously closer to the interpreted saturation curve of the single well actual logging (solid line in fig. 13), and has more obvious effect on poor second-class reservoirs, third-class reservoirs and fourth-class reservoirs (I, II marked by Facies is corresponding to the first-class reservoirs, III and IV are corresponding to the second-class reservoirs, V is corresponding to the third-class reservoirs, and VI is corresponding to the fourth-class reservoirs). The improved regression method improves the regression accuracy, and means that the calculation result of the original geological reserve model between wells controlled by the method is closer to the single-well logging result than the existing method, namely better accords with the actual underground condition.
The above examples are only used to illustrate the present invention, wherein the capillary pressure testing method, the rock core lithology, the porosity, the permeability, etc. may be varied, and all equivalent changes and modifications based on the technical solution of the present invention should not be excluded from the protection scope of the present invention.
Claims (8)
1. A method of calculating an average capillary pressure curve for a reservoir, comprising the steps of:
1) selecting a plurality of core samples on the same type of reservoir, performing a group of capillary pressure experiments on each core sample, and recording capillary pressure and wetting phase saturation data;
2) for each core sample, using the data recorded in step 1) with the formula respectivelyCarrying out dimensionless J function processing, and establishing a numerical value corresponding relation between the J function and the water saturation;
3) respectively using the data recorded in the step 2) to respectively use the formula for all core samples of the same reservoir Performing regression solution, wherein the parameter variation range limits A2>0,C2>0,B2Is greater than 0, and regression parameter A is obtained2、B2、C2And a first determination coefficientRespectively using the data recorded in the step 2) to respectively use the formula for all core samples of the same reservoirPerforming regression solution, wherein the parameter variation range limits A4>0,C4>0,0<B4<1, obtaining a regression parameter A4、B4、C4And a second determination coefficient
4) And (3) plotting the regression result and the data recorded in the step 2), selecting a formula with the regression result positioned in the range of the water saturation and the J function in the step 2), determining a formula with a coefficient closer to 1 as a basic formula of the average capillary pressure curve, and taking the regression parameters of all core samples of the same type of reservoir in the step 3) corresponding to the basic formula as the regression parameters corresponding to the average capillary pressure curve of the reservoir.
2. A method of calculating a reservoir average capillary pressure curve as claimed in claim 1, wherein said first and second determining coefficients are of the same type, i.e. both said first and second determining coefficients are of the same typeWherein SSR is regression sum of squares, SSE is residual sum of squares, SST is total sum of squares.
3. The method of calculating a reservoir average capillary pressure curve of claim 1, wherein the optimal formula is selected by:
3-1) taking the J function and water saturation relation data points of a plurality of dimensionless core samples of the same reservoir as true values, calculating the J function regression value corresponding to the water saturation by using the formula (2) and the parameters obtained by regression of the true values and the formula (2) as regression values 1, and recording the determination coefficient 1 as the regression value 1The formula (3) and the true valueThe regression value of the J function corresponding to the water saturation is calculated as a regression value 2 by using the regression formula (3) as the regression parameter obtained by regression, and the determination coefficient 2 is recorded asSelecting a corresponding formula with a regression value within the range of the true value as an optimal formula;
3-2) for the case that all 2 formulas are within the range of true values, the magnitude of the coefficient is determined by comparison:
if it isIs greater than or equal toThe basic formula of the average capillary pressure curve of the reservoir to be detected is a formula (2);
4. The method for calculating an average capillary pressure curve of a reservoir according to claim 3, wherein in the step 2), when the basic formula of the average capillary pressure curve of the reservoir to be calculated is formula (2), the average capillary pressure curve of the reservoir to be calculated isWherein A is2>0,C2>0,B2>0;
5. A method for calculating an average capillary pressure curve of a reservoir according to any one of claims 1 to 4, wherein the number of cores subjected to a capillary pressure experiment in the step 1) is more than or equal to 1.
6. A method of calculating an average capillary pressure curve for a reservoir according to any one of claims 1 to 4, wherein 1 set of capillary pressure experiments is performed on each core sample; and recording the water saturation, capillary pressure, wetting angle and two-phase interfacial tension in the capillary pressure experiment.
7. A method of calculating a reservoir average capillary pressure curve as claimed in any one of claims 1 to 6, wherein said regression parameter A4、A2、B4、B2Related to fractal dimension, inlet pressure, maximum capillary pressure, C4And C2Related to the fractal dimension.
8. The method of calculating a reservoir average capillary pressure curve of claim 7, wherein A is2、B2、C2Can be represented by formulaWherein A is2>0,C2>0,B2Is greater than 0; orWherein A is2>0,C2>0,B2More than 1, obtained by multivariate nonlinear regression; a is described4、B4、C4Can be represented by formulaWherein A is4>0,C4>0,0<B4<1 orWherein A is4>0,C4>0,B4Greater than 0, obtained by multivariate nonlinear regression, wet phase saturation SwNon-wetting phase saturation SnwHas a relationship of Sw+Snw=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911147846.7A CN110909466B (en) | 2019-11-21 | 2019-11-21 | Method for improving calculation accuracy of average capillary pressure curve of differential pore sorting reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911147846.7A CN110909466B (en) | 2019-11-21 | 2019-11-21 | Method for improving calculation accuracy of average capillary pressure curve of differential pore sorting reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110909466A true CN110909466A (en) | 2020-03-24 |
CN110909466B CN110909466B (en) | 2023-11-03 |
Family
ID=69818410
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911147846.7A Active CN110909466B (en) | 2019-11-21 | 2019-11-21 | Method for improving calculation accuracy of average capillary pressure curve of differential pore sorting reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110909466B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113092337A (en) * | 2021-04-08 | 2021-07-09 | 西南石油大学 | Method for establishing initial water saturation of compact rock core under in-situ condition |
CN113433035A (en) * | 2021-06-28 | 2021-09-24 | 西南石油大学 | Multi-capillary force curve averaging method based on multi-sample overall virtual measurement |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2634104A1 (en) * | 2007-06-08 | 2008-12-08 | University Of New Brunswick | Methods suitable for measuring capillary pressure and relative permeability curves of porous rocks |
CN105353419A (en) * | 2015-11-25 | 2016-02-24 | 中国地质大学(北京) | Nuclear magnetism capillary pressure curve construction method based on reservoir stratum classification |
CN108763648A (en) * | 2018-04-26 | 2018-11-06 | 中国石油大学(北京) | Method and apparatus based on nuclear magnetic resonance T2 distributed acquisition capillary pressure curves |
CN108983295A (en) * | 2018-07-20 | 2018-12-11 | 中国海洋石油集团有限公司 | A kind of reservoir critical property acquisition methods based on pore structure |
-
2019
- 2019-11-21 CN CN201911147846.7A patent/CN110909466B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2634104A1 (en) * | 2007-06-08 | 2008-12-08 | University Of New Brunswick | Methods suitable for measuring capillary pressure and relative permeability curves of porous rocks |
CN105353419A (en) * | 2015-11-25 | 2016-02-24 | 中国地质大学(北京) | Nuclear magnetism capillary pressure curve construction method based on reservoir stratum classification |
CN108763648A (en) * | 2018-04-26 | 2018-11-06 | 中国石油大学(北京) | Method and apparatus based on nuclear magnetic resonance T2 distributed acquisition capillary pressure curves |
CN108983295A (en) * | 2018-07-20 | 2018-12-11 | 中国海洋石油集团有限公司 | A kind of reservoir critical property acquisition methods based on pore structure |
Non-Patent Citations (5)
Title |
---|
W.S. XU: "A Prediction Model of the Capillary Pressure J-Function", PLOS ONE, pages 1 - 9 * |
代百祥,方小宇,郇金来,代云娇,甘永年: "利用J函数拟合法建立海上油田饱和度模型", 河北地质大学学报, pages 1 - 5 * |
刘义坤等: "毛管压力曲线和分形理论在储层分类中的应用", 岩性油气藏, vol. 26, no. 3, pages 89 - 92 * |
王庆勇;: "利用毛细管压力计算低渗储层可动水饱和度", 断块油气田, vol. 20, no. 01, pages 85 - 88 * |
贺承祖,华明琪: "储层孔隙结构的分形几何描述", 石油与天然气地质, vol. 19, no. 01, pages 15 - 23 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113092337A (en) * | 2021-04-08 | 2021-07-09 | 西南石油大学 | Method for establishing initial water saturation of compact rock core under in-situ condition |
CN113092337B (en) * | 2021-04-08 | 2022-01-28 | 西南石油大学 | Method for establishing initial water saturation of compact rock core under in-situ condition |
CN113433035A (en) * | 2021-06-28 | 2021-09-24 | 西南石油大学 | Multi-capillary force curve averaging method based on multi-sample overall virtual measurement |
CN113433035B (en) * | 2021-06-28 | 2022-04-15 | 西南石油大学 | Multi-capillary force curve averaging method based on multi-sample overall virtual measurement |
Also Published As
Publication number | Publication date |
---|---|
CN110909466B (en) | 2023-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105781539B (en) | A kind of densification oil and gas reservoir saturation degree well logging calculation method | |
CN107387051B (en) | Repeated fracturing well selection method for multi-stage fractured horizontal well with low-permeability heterogeneous oil reservoir | |
CN104134101A (en) | Low-seepage reservoir natural gas productivity prediction method | |
CN104806232B (en) | A kind of method for determining porosity lower limit of fracture | |
CN107291667B (en) | Method and system for determining communication degree between wells | |
CN105608334B (en) | A method of evaluation cigarette auxiliary material parameter influence degree | |
CN108982320B (en) | Method for calculating permeability of reservoir with complex pore structure by using particle size parameters | |
CN110909466A (en) | Method for improving calculation precision of differential pore sorting reservoir average capillary pressure curve | |
CN104278989A (en) | Method for obtaining saturability index of low porosity and low permeability reservoir | |
CN106897531B (en) | Quantitative evaluation method for permeability of low-permeability limestone reservoir | |
CN107132171A (en) | A kind of method for determining compact reservoir pore-size distribution based on pressure mercury N2 adsorption translocation data | |
CN110245686A (en) | A kind of lithology method for quickly identifying calculating quartzy percentage contents | |
CN105447298A (en) | Reservoir stratum analysis method and application thereof | |
CN106204302B (en) | A kind of method and application calculating original water saturation | |
CN106600440A (en) | Method for profile-control water-shutoff dynamic index well selection for low-permeability reservoir | |
CN105628584B (en) | The quantitatively characterizing method of the oily reservoir quality of densification based on combination weighting coefficient | |
CN107153733A (en) | A kind of capillary pressure curve approximating method of Complicated Pore Structures rock | |
CN106771071B (en) | A kind of sealing core drilling saturation correction method mutually seeped based on grease | |
CN114370269B (en) | Comprehensive determination method for physical property lower limit of effective reservoir of deep carbonate reservoir | |
CN106481315B (en) | Land sandstone oil reservoir individual well recoverable reserves quickly determines model and method for building up | |
CN112145165A (en) | Dynamic and static permeability conversion method for micro-crack-pore type reservoir | |
CN104374679A (en) | Method for determining permeability contribution value of reservoir | |
CN108427654B (en) | Rapid calculation method for silted storage capacity of medium-sized or over-sized check dam | |
CN109655394B (en) | Nuclear magnetic resonance T2 spectrum permeability calculation method under constraint of throat structure parameters | |
CN110322363A (en) | Shale gas reservoir reconstruction volume calculation method and system |
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 |