CN104751002A - Method for determining effective sandstone reservoir - Google Patents

Method for determining effective sandstone reservoir Download PDF

Info

Publication number
CN104751002A
CN104751002A CN201510171073.1A CN201510171073A CN104751002A CN 104751002 A CN104751002 A CN 104751002A CN 201510171073 A CN201510171073 A CN 201510171073A CN 104751002 A CN104751002 A CN 104751002A
Authority
CN
China
Prior art keywords
sandstone
porosity
buried depth
potential
interface
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.)
Pending
Application number
CN201510171073.1A
Other languages
Chinese (zh)
Inventor
庞雄奇
姜航
陈冬霞
姜福杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201510171073.1A priority Critical patent/CN104751002A/en
Publication of CN104751002A publication Critical patent/CN104751002A/en
Pending legal-status Critical Current

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The embodiment of the invention provides a method for determining effective sandstone reservoir. The method comprises the steps of building relationship between sandstone porosity and buried depth; building relationship between mudstone porosity and buried depth; building a calculation model of sandstone and mudstone pore throat radius; calculating interfacial potentials of sandstone and mudstone; calculating ratio of interfacial potential of mudstone to that of sandstone; building mudstone and sandstone interfacial potential ratio based buried depth distribution model; determining the maximum and minimum limits of the ratio distribution; calculating interfacial potential ratio of the mudstone and the sandstone a hydrocarbon reservoir and a dry layer; determining the critical value of interfacial potential ratio of mudstone to sandstone in the effective sandstone layer according to the interfacial potential ratio of mudstone to sandstone in the hydrocarbon reservoir and the dry layer. The method can overcome the technical problem of inaccurate determining of the effective sandstone reservoir caused by the nondeterminacy of subjectivity and determining standard in the effective sandstone reservoir determining method in the prior art, and reaches the technical effect of accurately and effectively determining the effective sandstone reservoir.

Description

Differentiate the method for effective sandstone reservoir
Technical field
The present invention relates to oil-gas exploration technical field, particularly a kind of method differentiating effective sandstone reservoir.
Background technology
Effective Reservoirs refers to and can preserve and seepage flow fluid (hydro carbons or local water), and extraction can have the reservoir of commercial value Liquid output (hydro carbons or hydro carbons mix with water) under technical conditions instantly.Effective Reservoirs is not equal to effective oil reservoir, and in Effective Reservoirs, the fluid of extraction both can be hydro carbons, and also can be water, therefore Effective Reservoirs comprises effective oil reservoir.Dried layer refers to that reservoir properties is poor, Liquid output (hydro carbons or water) lower than the rock stratum of dried layer output standard, be non-effective reservoir.The validity of reservoir is the concentrated expression of reservoir energy collecting power and filtration flow ability of immigrants, therefore differentiates that Effective Reservoirs is mainly by assessing the storage ability of oozing of reservoir, states with Critical porosity usually.At present, mainly effective sandstone reservoir is determined according to rock core Physical Property Analysis, formation testing and production test data etc.Mainly contain following several method:
1) empirical coefficient method, is mainly multiplied by the permeability limits as this oil field after 5% using field-wide mean permeability value, then adds up and the relational expression of matching porosity and permeability, thus asks for the critical pore angle value of Effective Reservoirs.The experience factor adopted in this method is too simple, determines that the foundation of Effective Reservoirs criterion is abundant not.
2) physical property formation testing method, mainly adopts the factor of porosity of formation testing conclusion data plotting oil reservoir, water layer and dried layer and permeability to cross figure, thus finds out the separatrix of oil-water-layer and dried layer factor of porosity, be defined as Effective Reservoirs critical pore angle value.This method be based upon there is a large amount of formation testing findings data and reservoir porosity, permeability data basis on, if sample point is less, the material point especially near the angle value of reservoir critical pore is less, so can produce larger error.
3) mud filtrate invasion method, mainly cross average for the reservoir of water-base drilling fluid core hole analysis original water saturation with the factor of porosity of respective layer, permeability figure respectively, the reservoir that drilling fluid can not invade substantially is non-effective reservoir, the reservoir of mud filtrate invasion is had to be Effective Reservoirs, Effective Reservoirs and non-effective reservoir can be determined according to this standard, this method be based upon there is water-base drilling fluid coring data basis on, if study area is this data not, then the method cannot be adopted.
4) tail-flick method, mainly does its cumulative frequency on factor of porosity energy F-histogram, and porosity value corresponding when frequency increases rapidly is Effective Reservoirs Critical porosity.This method is fuzzyyer for the restriction of " frequency increases rapidly " this index, therefore accurately cannot determine Effective Reservoirs.
5) shale index method, mainly has no commercial value according to the reservoir when shale index in sandstone is greater than 0.4, thus determines that sandstone porosity is now Effective Reservoirs Critical porosity.The key of this method is to ask accurate shale index critical point, and in sandstone, shale index is the percentage that dispersed shale volume accounts for total porosity, and it and net porosity all can be obtained by well logging interpretation, but practical operation difficulty is comparatively large, and uncertainty is also larger.
6) minimum effective pore throat method is the critical pore throat minimum value first determining Effective Reservoirs, then according to the figure that crosses of pore throat intermediate value and factor of porosity, limits porosity cutoff, thus determine Effective Reservoirs Critical porosity by minimum effective pore throat value.The shortcoming of this method is that asking for of minimum effective pore throat also exists larger difficulty and uncertainty.
7) factor of porosity-permeability intersection, that factor of porosity and permeability X plot are divided into three line segments: First Line sector hole porosity increases and permeability increase is very micro-, second line segment permeability increases with factor of porosity and obviously increases, 3rd line segment is that factor of porosity increase is very micro-and permeability sharply increases, and then factor of porosity corresponding for the turning point between the first line segment with the second line segment is defined as Effective Reservoirs Critical porosity.The shortcoming of this method does not have clear and definite method and principle to determine the first line segment and the second line segment, is therefore difficult to the turning point determining the first line segment and the second line segment exactly.
8) Purcell method, mainly utilize capillary pressure curve data, first the reservoir penetrating power between different aperture radial regions and accumulative penetrating power is calculated, pore radius corresponding when then reaching 99.99% by accumulative penetrating power is as effective pore throat radius lower limit, and the right factor of porosity of this pore throat radius is Effective Reservoirs critical pore angle value.This method needs to obtain the capillary pressure curve characterizing whole oil reservoir, but the capillary pressure curve of experimental determination only represents a bit in payzone, and therefore data point exists one-sidedness.
9) Permeability sensitivity method, mainly according to the susceptibility of permeability formation pressure, statistics ground permeability and the relation of zone permeability, when ground permeability is less, zone permeability relative to the change of ground infiltration rate be significantly less than ground infiltration rate larger time.Determine the Critical permeability value of Effective Reservoirs according to variation tendency turning point between the two, then according to the relation between factor of porosity and permeability, ask for corresponding Effective Reservoirs critical pore angle value.In practical operation, due to shale in the reservoir that permeability is very little or landwaste content high, multiplely educate flat or tabular venturi, confining pressure increases, small venturi can be caused to close, thus make the Effective Reservoirs Critical permeability determined by the method inaccurate, and there is comparatively big error in the variation relation fitting formula asking for ground infiltration rate and zone permeability in process, finally causes the Effective Reservoirs Critical porosity that obtains inaccurate.
10) throat distribution method is the number percent accounting for total pore throat according to the pore throat of different radii, determines Effective Reservoirs pore throat radius critical value, and then determines Effective Reservoirs critical pore angle value according to the relation of factor of porosity and pore throat radius.Owing to not having quantitative index to pore throat radius critical value, the intercepting of critical value has certain variability, therefore also has certain inexactness to asking for of Effective Reservoirs Critical porosity.
11) permeability-ratio curve and capillary pressure curve overlay technique, according to water saturation and relative permeability variation relation that is oily and water, the critical value of water saturation is determined according to the bottom flex point of the rel ative permeability curve to oil, on capillary pressure curve, find the pore throat radius value that this water saturation is corresponding again, the factor of porosity that pore throat radius is corresponding is therewith Effective Reservoirs Critical porosity.The theoretical property of this method is comparatively strong, but there is the shortcoming similar to Purcell method, and capillary pressure curve data only represent a point of underground, therefore there is the problem of data point one-sidedness.
Can be found by above-mentioned analysis, existingly determine that effective sandstone reservoir method mainly differentiates the Critical porosity of Effective Reservoirs.These methods all also exist certain defect, and what utilize these methods above-mentioned to obtain is all the physical property lower limit of Effective Reservoirs now, but not Pool-forming time Effective Reservoirs physical property boundary.Filling process is not considered in the research of reservoir oil-containing physical property lower limit now, and reservoir oil-containing physical property lower limit is that reservoir critical property a series of geological process after Pool-forming time runs up to result now now.Some reservoir has good physical property condition at Pool-forming time and becomes the useful space preserving oil gas, should be classified as in Effective Reservoirs category.Due to reservoir generation densification after Pool-forming time, physical property condition is deteriorated, and when adding up the physical property condition of these reservoirs now under condition, the critical physical property condition of Effective Reservoirs can be reduced.Therefore, it is obviously inappropriate for dividing Effective Reservoirs by statistics Effective Reservoirs physical property lower limit now.
Summary of the invention
Embodiments provide a kind of method differentiating effective sandstone reservoir, to solve the inaccurate technical matters of effective sandstone reservoir determined in prior art, the method comprises:
According to the sandstone porosity φ of the study area different buried depth obtained from the logging trace data of study area sand mudstone porosity φ m, draw sandstone porosity φ sloose point with buried depth change crosses figure and mudstone porosity φ mloose point with buried depth change crosses figure, and simulates mudstone porosity φ mwith the first relational expression of buried depth change;
According to the sandstone porosity φ obtained from the pressure mercury data of study area s1, permeability K s1with pore throat radius r s1, simulate the second relational expression between factor of porosity φ and permeability K, and pore throat radius r is with the 3rd relational expression of factor of porosity φ and permeability K change;
According to described second relational expression, calculate the sandstone porosity φ of different buried depth scorresponding permeability K swith mudstone porosity φ mcorresponding permeability K m;
According to described 3rd relational expression, calculate the sandstone porosity φ of different buried depth swith permeability K scorresponding pore throat radius r s, and mudstone porosity φ mwith permeability K mcorresponding pore throat radius r m;
According to the sandstone pore throat radius r of different buried depth swith mud stone pore throat radius r m, calculate the interface-potential of sandstone and the interface-potential of mud stone respectively, draw out the mud stone interface-potential of different buried depth and the ratio of the sandstone interface-potential distributed model with buried depth;
According to the well logging interpretation conclusion data of study area, obtain the sandstone porosity φ of hydrocarbon zone and dried layer s2with permeability K s2, then the sandstone pore throat radius r of hydrocarbon zone and dried layer is calculated according to described 3rd relational expression s2;
The mudstone porosity φ of hydrocarbon zone and dried layer place buried depth point is calculated according to described first relational expression m2, the mud stone permeability K of hydrocarbon zone and dried layer place buried depth point is calculated according to described second relational expression m2, and the mud stone pore throat radius r of hydrocarbon zone and dried layer place buried depth point is calculated according to described 3rd relational expression m2;
According to the sandstone pore throat radius r of described hydrocarbon zone and dried layer place buried depth point s2with mud stone pore throat radius r m2, calculate mud stone interface-potential and the sandstone interface-potential ratio of hydrocarbon zone and dried layer;
The mud stone interface-potential of the hydrocarbon zone calculated and dried layer and sandstone interface-potential ratio are projected in described distributed model, obtains the critical condition of separatrix as effective sandstone reservoir distribution of hydrocarbon zone and dried layer
In embodiments of the present invention, according to sandstone reservoir and the interface-potential ratio of mud stone around it, draw the Changing Pattern of effective sandstone reservoir properties lower limit with buried depth, thus determine the critical interfaces Potence ratio of each depth of burial effective sandstone reservoir, make the effective sandstone reservoir that obtains more reasonable accurately.And the method is for determining that the effective sandstone reservoir in petroclastic rock basin has general applicability, and easily operate, adopt the method quantitatively calculated simultaneously, eliminate the uncertainty in the subjectivity differentiated in effective sandstone reservoir method and discrimination standard, further, data selected by the method easily obtain, and the effective sandstone reservoir not having in geologic information stratum in order to prediction, therefore for scientific and reasonable planning with organize oil-gas exploration to have important directive significance.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:
Fig. 1 is the method flow diagram of the determination effective sandstone reservoir according to the embodiment of the present invention;
Fig. 2 is that storehouse car down warping region sandstone according to the embodiment of the present invention, mudstone porosity are with buried depth distributed model schematic diagram;
Fig. 3 is that storehouse car down warping region sandstone according to the embodiment of the present invention, mud stone pore throat radius are with buried depth distributed model schematic diagram;
Fig. 4 is with buried depth distributed model schematic diagram according to the storehouse car down warping region mud stone interface-potential of the embodiment of the present invention and sandstone interface-potential ratio;
Fig. 5 is storehouse car down warping region effective sandstone reservoir mud stone interface-potential according to the embodiment of the present invention and sandstone interface-potential ratio distribution schematic diagram.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with embodiment and accompanying drawing, the present invention is described in further details.At this, exemplary embodiment of the present invention and illustrating for explaining the present invention, but not as a limitation of the invention.
Inventor finds, during different buried depth, the critical physical property condition of effective sandstone reservoir is different, itself and buried depth have certain correlativity, on the basis of statistical study, think that effective sandstone reservoir porosity, permeability limits are not definite value, but with oil property and depth of burial, there is certain correlativity.Propose a kind of method determining effective sandstone reservoir in this example, the method is mainly considered can according to sandstone reservoir and the interface-potential ratio of mud stone around it, draw the Changing Pattern of effective sandstone reservoir properties lower limit with buried depth, thus determine the critical interfaces Potence ratio of each depth of burial effective sandstone reservoir, concrete, as shown in Figure 1, can comprise the following steps:
Step 101: according to the sandstone porosity φ of the study area different buried depth obtained from the logging trace data of study area sand mudstone porosity φ m, draw sandstone porosity φ sloose point with buried depth change crosses figure and mudstone porosity φ mloose point with buried depth change crosses figure, and simulates mudstone porosity φ mwith the first relational expression of buried depth change;
In this step, sandstone can be the sandstone of siltstone, packsand, middle sandstone and sandstone grit four kinds of particle diameters.
Concrete, the mudstone porosity φ that the matching that matching obtains obtains mthe first relational expression with buried depth change can be expressed as:
φ m = a × e - ( Z + c b ) 2 + d × e - ( Z + f e ) 2
Wherein, φ mrepresent mudstone porosity, unit is that %, Z represent buried depth, and unit is m, a, b, c, d, e, f represent fitting parameter, i.e. the process of matching is exactly determine the process of fitting parameter a, b, c, d, e, f.
Step 102: according to the sandstone porosity φ obtained from the pressure mercury data of study area s1, permeability K s1with pore throat radius r s1, simulate the second relational expression between factor of porosity φ and permeability K, and pore throat radius r is with the 3rd relational expression of factor of porosity φ and permeability K change;
Concrete, in this step, the second relational expression between the factor of porosity φ simulated and permeability K can be expressed as:
K＝ge hφ
Wherein, K represents permeability, and unit is 10 -3um 2, φ represents factor of porosity, and unit is %, g, h represents fitting parameter.
The pore throat radius r simulated can be expressed as with the 3rd relational expression of factor of porosity φ and permeability K change:
lgr＝mlgφ+nlgK-p
Wherein, r represents pore throat radius, and unit is that um, φ represent factor of porosity, and unit is that %, K represent permeability, and unit is 10 -3um 2, m, n, p represent fitting parameter.
It is exactly the process determining fitting parameter g, h and m, n, p in the process of above-mentioned matching.
Although it should be noted that in the process of matching it is according to sandstone porosity φ s1, permeability K s1with pore throat radius r s1the second relational expression that matching obtains and the 3rd relational expression, but this relational expression can simultaneous adaptation in sandstone and mud stone, therefore, just can calculate the sandstone porosity φ of different buried depth according to this second relational expression and the 3rd relational expression scorresponding permeability K swith mudstone porosity φ mcorresponding permeability K m, and the sandstone porosity φ of different buried depth swith permeability K scorresponding pore throat radius r s, and mudstone porosity φ mwith permeability K mcorresponding pore throat radius r m.
Step 103: according to described second relational expression, calculates the sandstone porosity φ of different buried depth scorresponding permeability K swith mudstone porosity φ mcorresponding permeability K m;
Step 104: according to described 3rd relational expression, calculate the sandstone porosity φ of different buried depth swith permeability K scorresponding pore throat radius r s, and mudstone porosity φ mwith permeability K mcorresponding pore throat radius r m;
Step 105: according to the sandstone pore throat radius r of different buried depth swith mud stone pore throat radius r m, calculate the mud stone interface-potential of different buried depth and the ratio of the sandstone interface-potential distributed model with buried depth;
Because mud stone interface-potential and sandstone interface-potential ratio are numerically equivalent to the ratio of sandstone pore throat radius and mud stone pore throat radius, the mud stone interface-potential of different buried depth and the ratio of sandstone interface-potential can be obtained according to following formulae discovery in this step:
Φ ratio = 2 σ cos θ r m / 2 σ cos θ r s = r s r m
Wherein, Φ ratiorepresent mud stone interface-potential and sandstone interface-potential ratio, dimensionless, σ represents interfacial tension, and unit is that N/m, θ represent wetting angle, and unit is °, r mrepresent mud stone pore throat radius, unit is um, r srepresent sandstone pore throat radius, unit is um.
The mud stone interface-potential of the different buried depth obtained in this step and the ratio of sandstone interface-potential have maximum boundary and minimum border with the distributed model of buried depth, this mainly because: at drafting sandstone porosity φ sloose point with buried depth change crosses figure and mudstone porosity φ mloose point with buried depth change crosses in the process of figure, be using sandstone in the factor of porosity minimum value of each buried depth and maximal value as the minimum border of sandstone porosity and maximum boundary, using the mudstone porosity value of the factor of porosity mean value of each for mud stone buried depth as this buried depth point, to draw out sandstone porosity φ sloose point with buried depth change crosses figure and mudstone porosity φ mloose point with buried depth change crosses figure, then according to described 3rd relational expression, calculates the sandstone porosity φ of different buried depth swith permeability K scorresponding pore throat radius r s, and mudstone porosity φ mwith permeability K mcorresponding pore throat radius r mprocess in, sandstone constitutes the minimum border of sandstone pore throat radius in the minimum value of the sandstone pore throat radius of each buried depth, sandstone constitutes the maximum boundary of sandstone pore throat radius in the maximal value of the sandstone pore throat radius of each buried depth, and mud stone only exists a mud stone pore throat radius value at each buried depth.Therefore, the mud stone interface-potential of the different buried depth obtained and the ratio of sandstone interface-potential is finally made to be have maximum boundary and minimum border with the distributed model of buried depth.
Step 106: according to the well logging interpretation conclusion data of study area, obtains the sandstone porosity φ of hydrocarbon zone and dried layer s2with permeability K s2, then the sandstone pore throat radius r of hydrocarbon zone and dried layer is calculated according to described 3rd relational expression s2;
So-called hydrocarbon zone can comprise with lower floor: oil reservoir, gas-bearing formation, water layer, oil-containing water layer, moisture oil reservoir, gassiness water layer, moisture gas-bearing formation, oil-water common-layer, the same layer of air water, the same layer of oil gas, difference hydrocarbon zone, difference gas-bearing formation, difference oil reservoir, hydrocarbon zone and dried layer all belong to sandstone reservoir.
Step 107: the mudstone porosity φ calculating hydrocarbon zone and dried layer place buried depth point according to described first relational expression m2, the mud stone permeability K of hydrocarbon zone and dried layer place buried depth point is calculated according to described second relational expression m2, and the mud stone pore throat radius r of hydrocarbon zone and dried layer place buried depth point is calculated according to described 3rd relational expression m2;
Step 108: according to the sandstone pore throat radius r of described hydrocarbon zone and dried layer place buried depth point s2with mud stone pore throat radius r m2, calculate mud stone interface-potential and the sandstone interface-potential ratio of hydrocarbon zone and dried layer;
Step 109: the mud stone interface-potential of the hydrocarbon zone calculated and dried layer and sandstone interface-potential ratio are projected in described distributed model, obtains the critical condition of separatrix as effective sandstone reservoir distribution of hydrocarbon zone and dried layer.
In this step, concrete, the determination of critical condition can be carried out: the mud stone interface-potential of hydrocarbon zone and dried layer and sandstone interface-potential ratio have separatrix with the distribution range of buried depth according to following rule, wherein, be mud stone interface-potential and the sandstone interface-potential ratio distributive province of hydrocarbon zone on described marginal right side, be mud stone interface-potential and the sandstone interface-potential ratio distributive province of dried layer in described marginal left side, the reservoir that the mud stone interface-potential of each buried depth point and sandstone interface-potential ratio are greater than this buried depth point critical value is effective sandstone reservoir.
Concrete, above-mentioned sandstone can be the sandstone in petroclastic rock basin, and above-mentioned mud stone is exactly the compactness rock stratum of playing capping oil-gas reactivation around the sandstone reservoir of petroclastic rock basin.
In the above-described embodiments, according to sandstone reservoir and the interface-potential ratio of mud stone around it, draw the Changing Pattern of effective sandstone reservoir properties lower limit with buried depth, thus determine the critical interfaces Potence ratio of each depth of burial effective sandstone reservoir, thus make the effective sandstone reservoir that obtains more reasonable accurately.And the method is for determining that the sandstone effective sandstone reservoir in petroclastic rock basin has general applicability, and easily operate, adopt the method quantitatively calculated simultaneously, eliminate the uncertainty in the subjectivity differentiated in effective sandstone reservoir method and discrimination standard, further, data selected by the method easily obtain, and the effective sandstone reservoir not having in geologic information stratum in order to prediction, therefore for scientific and reasonable planning with organize oil-gas exploration to have important directive significance.
In order to better the present invention is described, determine that the method for effective sandstone reservoir is specifically described below in conjunction with a specific embodiment to above-mentioned, but it should be noted that this specific embodiment carries out, in order to better the present invention is described, not forming inappropriate limitation of the present invention.
In this example, for Kuche Depression of Talimu Basin to how determining that the critical interfaces Potence ratio of storehouse car down warping region effective sandstone reservoir is specifically described, can comprise the following steps:
Step 1: set up storehouse car down warping region sandstone, model that mudstone porosity changes with buried depth, specifically comprises:
S1: the relation that the factor of porosity adding up sandstone changes with buried depth, generates scatter diagram.Concrete foundation is: the sandstone porosity of the same degree of depth has certain interval distribution range in the horizontal, the factor of porosity minimum value of all buried depth points constitutes the minimum border of sandstone porosity, and the factor of porosity maximal value of all buried depth points constitutes the maximum boundary of sandstone porosity;
S2: the relation that statistics mudstone porosity changes with buried depth, generate scatter diagram, concrete, can get the mudstone porosity value of factor of porosity mean value as this buried depth point of each buried depth point, matching obtains the relational expression that the mudstone porosity as formula 1 changes with buried depth:
φ m = 4551 × e - ( Z + 5979 2409 ) 2 + 8.363 × 10 17 × e - ( Z + 36660 58160 ) 2 (formula 1)
Wherein, φ mrepresent mudstone porosity, unit is that %, Z represent buried depth, and unit is m.
S3: utilize above-mentioned formula 1 to calculate the mudstone porosity at buried depth point place, each sandstone place, and the sandstone porosity corresponding to each buried depth, mudstone porosity are generated the mode chart that storehouse car down warping region sandstone, mudstone porosity as shown in Figure 2 change with buried depth.
Step 2: the permeability of statistics sandstone and porosity data, carries out matching to both sides relation and obtains, if the permeability of the storehouse car down warping region sandstone of formula 2 is with the distributed model of factor of porosity, just can calculating permeability corresponding to mudstone porosity according to formula 2.
K s = 0.0047 e 0.5292 φ s (formula 2)
R 2＝0.9999
Wherein, Φ srepresent sandstone porosity, unit is %, K srepresent Permeability of Sandstone, unit is 10 -3um 2, R 2be the index of Trendline fitting degree, its numerical values recited can reflect the fitting degree between the estimated value of Trendline and corresponding real data, and fitting degree is higher, and the reliability of Trendline is higher.
Step 3: determine the pore throat radius of sandstone, the pore throat radius of mud stone.
Concrete, can comprise:
S1: add up the factor of porosity of car area, storehouse sandstone, permeability and pore throat radius data, and matching obtains the governing equation of sandstone pore throat radius as shown in Equation 3 with factor of porosity, permeability variation, just can calculate the pore throat radius of mud stone according to formula 3:
Lgr s=0.171lg φ s+ 0.4123lgK s-0.4176 (formula 3)
R 2＝0.8343
Wherein, r srepresent sandstone pore throat radius, unit is um, φ srepresent sandstone porosity, unit is %, K srepresent Permeability of Sandstone, unit is 10 -3um 2.
Although it should be noted that formula 2 and 3 obtains according to sandstone data fitting, the characteristic simultaneous adaptation of formula, in sandstone and mud stone, therefore can calculate the parameter of mud stone according to above-mentioned formula 2 and 3.
S2: utilize formula 3 to calculate the pore throat radius of sandstone, utilizes formula 2 to calculate the mud stone permeability at buried depth point place, each sandstone place, and then calculates the mud stone pore throat radius of each buried depth point according to formula 3.Be illustrated in figure 3 the pore throat radius schematic diagram of mud stone and sandstone, the sandstone pore throat radius minimum value of each buried depth constitutes the minimum border of sandstone pore throat radius jointly, the sandstone pore throat radius maximal value of each buried depth constitutes the maximum boundary of sandstone pore throat radius jointly, two borders constitute the distributed model of sandstone pore throat radius with buried depth jointly, obtain each buried depth point according to formula 1 and only have mudstone porosity data, the mud stone therefore obtaining each buried depth point place only has a pore throat radius value.
Step 4: calculate mud stone interface-potential and sandstone interface-potential ratio.
Concrete, mud stone interface-potential and sandstone interface-potential ratio can be calculated according to following formula 4, this ratio is numerically equivalent to the ratio of sandstone pore throat radius and mud stone pore throat radius, according to sandstone pore throat radius, the mud stone pore throat radius of each buried depth point that step 3 is asked for, make its ratio, the mud stone interface-potential of different buried depth as shown in Figure 4 and the ratio of sandstone interface-potential can be obtained:
Φ ratio = 2 σ cos θ r m / 2 σ cos θ r s = r s r m (formula 4)
Wherein, Φ ratiorepresent the ratio of mud stone interface-potential and sandstone interface-potential, dimensionless, σ represents interfacial tension, and unit is that N/m, θ represent wetting angle, and unit is °, r nrepresent mud stone pore throat radius, unit is um, r srepresent sandstone pore throat radius, unit is um.
Step 5: statistics storehouse car down warping region Region of Oil-gas hides the porosity and permeability of each depth point, and the oil gas corresponding to each depth point explains conclusion.Concrete, this area's oil gas explains that conclusion mainly can comprise: oil reservoir, gas-bearing formation, water layer, dried layer, oil-containing water layer, moisture oil reservoir, gassiness water layer, moisture gas-bearing formation, oil-water common-layer, the same layer of air water, the same layer of oil gas, difference hydrocarbon zone, difference gas-bearing formation, difference oil reservoir, dried layer can be classified as a class, explain that conclusion is classified as another kind of by all the other 13 kinds, be referred to as hydrocarbon zone.
Step 6: the model utilizing step 3 to set up calculates the pore throat radius of hydrocarbon zone and dried layer, recycle the distributed model of mudstone porosity with buried depth of step 1 foundation, calculate the mudstone porosity of hydrocarbon zone and dried layer place buried depth point, then calculate the pore throat radius of mud stone according to step 3, then calculate mud stone interface-potential and the sandstone interface-potential ratio of hydrocarbon zone and dried layer.
Step 7: the mud stone interface-potential of hydrocarbon zone and dried layer and sandstone interface-potential ratio data point are thrown on mode chart that storehouse car area mud stone interface-potential and sandstone interface-potential ratio changes with buried depth, in the mud stone interface-potential that all data points (comprising hydrocarbon zone and dried layer) are all determined in step 4 and the scope that the minimum and maximum border of sandstone interface-potential ratio is drawn a circle to approve.Meanwhile, be distributed in the data point set of dried layer near the minimum border of mud stone interface-potential and sandstone interface-potential ratio, in the data point set of hydrocarbon zone, be distributed in mud stone interface-potential and the larger region of sandstone interface-potential ratio.Have obvious separatrix as shown in Figure 5 between the mud stone interface-potential of hydrocarbon zone and dried layer and sandstone interface-potential ratio data, this separatrix is the critical condition of effective sandstone reservoir distribution.The reservoir that each depth of burial point mud stone interface-potential and sandstone interface-potential ratio are greater than this buried depth point critical value is effective sandstone reservoir.
The mud stone of the storehouse car down warping region effective sandstone reservoir finally obtained and sandstone interface-potential ratio critical value can be expressed as with the computing formula of depth of burial:
Φ ratio'=1.93e 0.00012Z(formula 5)
Wherein, Φ ratio' representing effective sandstone reservoir mud stone and sandstone interface-potential ratio critical value, dimensionless, Z represents depth of burial, and unit is m.
In another embodiment, additionally provide a kind of software, this software is for performing the technical scheme described in above-described embodiment and preferred implementation.In another embodiment, additionally provide a kind of storage medium, store above-mentioned software in this storage medium, this storage medium includes but not limited to: CD, floppy disk, hard disk, scratch pad memory etc.
From above description, can find out, the embodiment of the present invention achieves following technique effect: according to sandstone reservoir and the interface-potential ratio of mud stone around it, draw the Changing Pattern of effective sandstone reservoir properties lower limit with buried depth, thus determine the critical interfaces Potence ratio of each depth of burial effective sandstone reservoir, make the effective sandstone reservoir that obtains more reasonable accurately.And the method is for determining that the sandstone effective sandstone reservoir in petroclastic rock basin has general applicability, and easily operate, adopt the method quantitatively calculated simultaneously, eliminate the uncertainty in the subjectivity differentiated in effective sandstone reservoir method and discrimination standard, further, data selected by the method easily obtain, and the effective sandstone reservoir not having in geologic information stratum in order to prediction, therefore for scientific and reasonable planning with organize oil-gas exploration to have important directive significance.
Obviously, those skilled in the art should be understood that, each module of the above-mentioned embodiment of the present invention or each step can realize with general calculation element, they can concentrate on single calculation element, or be distributed on network that multiple calculation element forms, alternatively, they can realize with the executable program code of calculation element, thus, they can be stored and be performed by calculation element in the storage device, and in some cases, step shown or described by can performing with the order be different from herein, or they are made into each integrated circuit modules respectively, or the multiple module in them or step are made into single integrated circuit module to realize.Like this, the embodiment of the present invention is not restricted to any specific hardware and software combination.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (11)

1. differentiate a method for effective sandstone reservoir, it is characterized in that, comprising:
According to the sandstone porosity φ of the study area different buried depth obtained from the logging trace data of study area sand mudstone porosity φ m, draw sandstone porosity φ sloose point with buried depth change crosses figure and mudstone porosity φ mloose point with buried depth change crosses figure, and simulates mudstone porosity φ mwith the first relational expression of buried depth change;
According to the sandstone porosity φ obtained from the pressure mercury data of study area s1, permeability K s1with pore throat radius r s1, simulate the second relational expression between factor of porosity φ and permeability K, and pore throat radius r is with the 3rd relational expression of factor of porosity φ and permeability K change;
According to described second relational expression, calculate the sandstone porosity φ of different buried depth scorresponding permeability K swith mudstone porosity φ mcorresponding permeability K m;
According to described 3rd relational expression, calculate the sandstone porosity φ of different buried depth swith permeability K scorresponding pore throat radius r s, and mudstone porosity φ mwith permeability K mcorresponding pore throat radius r m;
According to the sandstone pore throat radius r of different buried depth swith mud stone pore throat radius r m, calculate the interface-potential of sandstone and the interface-potential of mud stone respectively, draw out the mud stone interface-potential of different buried depth and the ratio of the sandstone interface-potential distributed model with buried depth;
According to the well logging interpretation conclusion data of study area, obtain the sandstone porosity φ of hydrocarbon zone and dried layer s2with permeability K s2, then the sandstone pore throat radius r of hydrocarbon zone and dried layer is calculated according to described 3rd relational expression s2;
The mudstone porosity φ of hydrocarbon zone and dried layer place buried depth point is calculated according to described first relational expression m2, the mud stone permeability K of hydrocarbon zone and dried layer place buried depth point is calculated according to described second relational expression m2, and the mud stone pore throat radius r of hydrocarbon zone and dried layer place buried depth point is calculated according to described 3rd relational expression m2;
According to the sandstone pore throat radius r of described hydrocarbon zone and dried layer place buried depth point s2with mud stone pore throat radius r m2, calculate mud stone interface-potential and the sandstone interface-potential ratio of hydrocarbon zone and dried layer;
The mud stone interface-potential of the hydrocarbon zone calculated and dried layer and sandstone interface-potential ratio are projected in described distributed model, obtains the critical condition of separatrix as effective sandstone reservoir distribution of hydrocarbon zone and dried layer.
2. method according to claim 1, is characterized in that, described mud stone is the compactness rock stratum of playing capping oil-gas reactivation around the sandstone reservoir of petroclastic rock basin.
3. method according to claim 1, it is characterized in that, described hydrocarbon zone comprise following one of at least: oil reservoir, gas-bearing formation, water layer, oil-containing water layer, moisture oil reservoir, gassiness water layer, moisture gas-bearing formation, oil-water common-layer, the same layer of air water, the same layer of oil gas, difference hydrocarbon zone, difference gas-bearing formation, difference oil reservoir.
4. method according to claim 1, is characterized in that, draws sandstone porosity φ sloose point with buried depth change crosses figure and mudstone porosity φ mloose point with buried depth change crosses figure, comprising:
Using sandstone in the factor of porosity minimum value of each buried depth and maximal value as the minimum border of sandstone porosity and maximum boundary, using the mudstone porosity value of the factor of porosity mean value of each for mud stone buried depth as this buried depth point, to draw out sandstone porosity φ sloose point with buried depth change crosses figure and mudstone porosity φ mloose point with buried depth change crosses figure.
5. method according to claim 4, is characterized in that, according to described 3rd relational expression, calculates the sandstone porosity φ of different buried depth swith permeability K scorresponding pore throat radius r s, and mudstone porosity φ mwith permeability K mcorresponding pore throat radius r mprocess in, sandstone constitutes the minimum border of sandstone pore throat radius in the minimum value of the sandstone pore throat radius of each buried depth, sandstone constitutes the maximum boundary of sandstone pore throat radius in the maximal value of the sandstone pore throat radius of each buried depth, and mud stone only exists a mud stone pore throat radius value at each buried depth.
6. method according to claim 5, is characterized in that, the mud stone interface-potential of the different buried depth obtained and the ratio of sandstone interface-potential have maximum boundary and minimum border with the distributed model of buried depth.
7. method according to claim 1, is characterized in that, according to the sandstone pore throat radius r of different buried depth swith mud stone pore throat radius r m, the ratio of the mud stone interface-potential and sandstone interface-potential of drawing out different buried depth, with the distributed model of buried depth, comprising:
The mud stone interface-potential of different buried depth and the ratio of sandstone interface-potential is obtained according to following formulae discovery:
Φ ratio = 2 σ cos θ r m / 2 σ cos θ r s = r s r m
Wherein, Φ ratiorepresent mud stone interface-potential and sandstone interface-potential ratio, dimensionless, σ represents interfacial tension, and unit is that N/m, θ represent wetting angle, and unit is °, r mrepresent mud stone pore throat radius, unit is um, r srepresent sandstone pore throat radius, unit is um.
8. method according to claim 1, it is characterized in that, the mud stone interface-potential of the hydrocarbon zone calculated and dried layer and sandstone interface-potential ratio are projected in described distributed model, the separatrix obtaining hydrocarbon zone and dried layer, as the critical condition of effective sandstone reservoir distribution, comprising:
The critical condition of effective sandstone reservoir distribution is determined: the mud stone interface-potential of hydrocarbon zone and dried layer and sandstone interface-potential ratio have separatrix with the distribution range of buried depth according to following rule, wherein, be mud stone interface-potential and the sandstone interface-potential ratio distributive province of hydrocarbon zone on described marginal right side, be mud stone interface-potential and the sandstone interface-potential ratio distributive province of dried layer in described marginal left side, wherein, the mud stone interface-potential that described separatrix is corresponding and sandstone interface-potential ratio are the standards differentiating sandstone reservoir, the reservoir that the mud stone interface-potential of each buried depth point and sandstone interface-potential ratio are greater than this buried depth point critical value is effective sandstone reservoir.
9. method according to any one of claim 1 to 8, is characterized in that, the mudstone porosity φ that matching obtains mthe first relational expression with buried depth change is:
φ m = a × e - ( Z + c b ) 2 + d × e - ( Z + f e ) 2
Wherein, φ mrepresent mudstone porosity, unit is that %, Z represent buried depth, and unit is m, a, b, c, d, e, f represent fitting parameter.
10. method according to any one of claim 1 to 8, is characterized in that, the second relational expression between the factor of porosity φ simulated and permeability K is:
K＝ge hφ
Wherein, K represents permeability, and unit is 10 -3um 2, φ represents factor of porosity, and unit is %, g, h represents fitting parameter.
11. methods according to any one of claim 1 to 8, is characterized in that, the pore throat radius r simulated with the 3rd relational expression of factor of porosity φ and permeability K change is:
lg r＝m lg φ+n lg K-p
Wherein, r represents pore throat radius, and unit is that um, φ represent factor of porosity, and unit is that %, K represent permeability, and unit is 10 -3um 2, m, n, p represent fitting parameter.
CN201510171073.1A 2015-04-10 2015-04-10 Method for determining effective sandstone reservoir Pending CN104751002A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510171073.1A CN104751002A (en) 2015-04-10 2015-04-10 Method for determining effective sandstone reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510171073.1A CN104751002A (en) 2015-04-10 2015-04-10 Method for determining effective sandstone reservoir

Publications (1)

Publication Number Publication Date
CN104751002A true CN104751002A (en) 2015-07-01

Family

ID=53590677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510171073.1A Pending CN104751002A (en) 2015-04-10 2015-04-10 Method for determining effective sandstone reservoir

Country Status (1)

Country Link
CN (1) CN104751002A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089659A (en) * 2015-07-07 2015-11-25 中国石油天然气股份有限公司 Conglomerate oil reservoir seepage unit recognizing method
CN105386757A (en) * 2015-12-21 2016-03-09 中国石油大学(华东) Clastic rock reservoir porosity normal compaction trend line determination method
CN106368694A (en) * 2016-11-16 2017-02-01 中国海洋石油总公司 Porosity evolution recovery and physical property prediction method for complex area reservoir
CN110107283A (en) * 2019-04-26 2019-08-09 中国石油化工股份有限公司 The method and system of Accumulation zone prediction are carried out according to potential energy value and physics value
CN113887078A (en) * 2021-10-26 2022-01-04 中国科学院地质与地球物理研究所 Deep-ultra-deep sea-facies carbonate reservoir occurrence depth lower limit prediction method
CN114396257A (en) * 2021-12-31 2022-04-26 中国石油大学(北京) Three-dimensional stratum pressure prediction method and device based on seismic data

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103645518A (en) * 2013-12-19 2014-03-19 中国石油大学(北京) Method for confirming sand lens gathering oil gas critical conditions
CN103678861A (en) * 2013-05-15 2014-03-26 中国石油大学(北京) Method for determining reservoir forming termination depth and reservoir forming range of tight sandstone gas reservoir
CN103926632A (en) * 2013-01-10 2014-07-16 中国石油大学(北京) Method for determining hydrocarbon accumulation threshold of hydrocarbon containing basin fragmentary rock

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103926632A (en) * 2013-01-10 2014-07-16 中国石油大学(北京) Method for determining hydrocarbon accumulation threshold of hydrocarbon containing basin fragmentary rock
CN103678861A (en) * 2013-05-15 2014-03-26 中国石油大学(北京) Method for determining reservoir forming termination depth and reservoir forming range of tight sandstone gas reservoir
CN103645518A (en) * 2013-12-19 2014-03-19 中国石油大学(北京) Method for confirming sand lens gathering oil gas critical conditions

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
姜航等: "基于毛细管力的有效储层物性下限判别", 《地质论评》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089659A (en) * 2015-07-07 2015-11-25 中国石油天然气股份有限公司 Conglomerate oil reservoir seepage unit recognizing method
CN105089659B (en) * 2015-07-07 2017-10-17 中国石油天然气股份有限公司 A kind of Conglomerate Reservoir permeable unit recognition methods
CN105386757A (en) * 2015-12-21 2016-03-09 中国石油大学(华东) Clastic rock reservoir porosity normal compaction trend line determination method
CN106368694A (en) * 2016-11-16 2017-02-01 中国海洋石油总公司 Porosity evolution recovery and physical property prediction method for complex area reservoir
CN106368694B (en) * 2016-11-16 2019-05-21 中国海洋石油总公司 One kind develops for complicated field reservoir pore space and restores and physical property prediction technique
CN110107283A (en) * 2019-04-26 2019-08-09 中国石油化工股份有限公司 The method and system of Accumulation zone prediction are carried out according to potential energy value and physics value
CN113887078A (en) * 2021-10-26 2022-01-04 中国科学院地质与地球物理研究所 Deep-ultra-deep sea-facies carbonate reservoir occurrence depth lower limit prediction method
CN113887078B (en) * 2021-10-26 2022-03-01 中国科学院地质与地球物理研究所 Deep-ultra-deep sea-facies carbonate reservoir occurrence depth lower limit prediction method
US11435494B1 (en) 2021-10-26 2022-09-06 Institute Of Geology And Geophysics, Chinese Academy Of Sciences Method for predicting oil accumulation depth limit of deep and ultra-deep marine carbonate reservoirs
CN114396257A (en) * 2021-12-31 2022-04-26 中国石油大学(北京) Three-dimensional stratum pressure prediction method and device based on seismic data
CN114396257B (en) * 2021-12-31 2023-10-31 中国石油大学(北京) Three-dimensional stratum pressure prediction method and device based on seismic data

Similar Documents

Publication Publication Date Title
CN104751002A (en) Method for determining effective sandstone reservoir
Lieu et al. The robustness of Sr/Y and La/Yb as proxies for crust thickness in modern arcs
Nordbotten et al. Uncertainties in practical simulation of CO2 storage
CN104712329B (en) A kind of computation model of mud shale hydrocarbon saturation
US10641923B2 (en) Method for characterizing and exploiting a subterranean formation comprising a network of fractures
CN103774988B (en) Real-time while-drilling gamma forward modeling geo-steering drilling method
CN105701319B (en) A kind of microfacies modeling method under horizontal well rule well pattern
CN103114840A (en) Calculating method and device of organic carbon content of high-over high mature shale
Thiele et al. Polymer-flood modeling using streamlines
CN102749647B (en) Quantitative method for space of holes in reservoir
CN106503295B (en) A kind of method and device for explaining oil field Water Flooding Layer using state-space model
CN104632202A (en) Method and device for determining dry clay three-porosity logging parameter values
Reeves et al. Development and Calibration of Dual-Permeability Flow Models with DiscontinuousFault Networks
CN106355571A (en) Determination method and device for dolostone reservoir quality
Al-Fatlawi et al. A new practical method for predicting equivalent drainage area of well in tight gas reservoirs
CN111177946B (en) Simulation calculation method for total skin and compacted zone skin coefficients of high-pressure dense gas reservoir perforation
CN111155980B (en) Water flow dominant channel identification method and device
Suarez et al. Fracturing-to-Production Simulation Approach for Completion Optimization in the Vaca Muerta Shale
CN107766689B (en) Method for establishing reservoir permeability time-varying model for developing dynamic constraint
CN113960288B (en) Quantitative evaluation method for source-sink system
Jambayev Discrete fracture network modeling for a carbonate reservoir
CN110954942A (en) Method and system for making chart with clay content identification function
Olson et al. Reservoir characterization of the giant Hugoton gas field, Kansas
Garamhegyi et al. Revision of archive recovery tests using analytical and numerical methods on thermal water wells in sandstone and fractured carbonate aquifers in the vicinity of Budapest, Hungary
Nielsen et al. Pilot phase of the Aguada Federal Block, black-oil window

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150701