CN106321087B - A method of obtaining lithostratigraphy factor - Google Patents

A method of obtaining lithostratigraphy factor Download PDF

Info

Publication number
CN106321087B
CN106321087B CN201510388742.0A CN201510388742A CN106321087B CN 106321087 B CN106321087 B CN 106321087B CN 201510388742 A CN201510388742 A CN 201510388742A CN 106321087 B CN106321087 B CN 106321087B
Authority
CN
China
Prior art keywords
porosity
factor
lithostratigraphy
rock
improvement type
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510388742.0A
Other languages
Chinese (zh)
Other versions
CN106321087A (en
Inventor
张军
李军
苏俊磊
武清钊
南泽宇
郝士博
王晓畅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201510388742.0A priority Critical patent/CN106321087B/en
Publication of CN106321087A publication Critical patent/CN106321087A/en
Application granted granted Critical
Publication of CN106321087B publication Critical patent/CN106321087B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of methods for obtaining lithostratigraphy factor, the lithostratigraphy factor is obtained according to the corresponding porosity calculation of two kinds of different pore structures of lithostratigraphy factor influence degree, it the described method comprises the following steps: Archie formula improvement type is constructed according to Archie formula, in the Archie formula improvement type, lithostratigraphy factor porosity corresponding to described two pore structures is related;The lithostratigraphy factor is obtained according to the Archie formula improvement type.Compared with prior art, the available lithostratigraphy factor for more meeting practical geological condition of the present invention, to substantially increase the accuracy of geological prospecting result, the oil-gas exploration and development for after provides more accurate strong support.

Description

A method of obtaining lithostratigraphy factor
Technical field
The present invention relates to geological exploration fields, in particular relate to a kind of method for obtaining lithostratigraphy factor.
Background technique
Determine that lithostratigraphy factor is an important content in the evaluation of oil log interpretation saturation degree.In the prior art In, determine lithostratigraphy factor generally use U.S. Well Logging Engineering teacher G.E.Archie 1942 deliver achievement (referred to as Ah The odd formula of that).
G.E.Archie is had found by a large amount of rock-electric test data analysis: under log-log coordinate, " pure " sandstone rock Stone formation factor is in a linear relationship with porosity, is represented by formula (1):
Wherein: F is lithostratigraphy factor, dimensionless;
R0For resistivity of the rock sample 100% when aqueous, unit is ohm meter;
RwFor pore water resistivity, ohm meter;
φ is rock porosity, decimal;
A, m is constant related with rock property.
The research achievement of G.E.Archie is based primarily upon the rock core sample porosity for being derived from United States Gulf area, resistivity is surveyed It measures as a result, rock sample porosity ranges are between 10%~40%.In recent years, with for tight sand (with porosity less than 10% Based on) petrophysics experiment research increase, a large amount of experimental data shows: under log-log coordinate, sandstone lithostratigraphy because Element and porosity are no longer single linear relationships, but piecewise linearity is presented.In some cases, Archie formula can not be very Good is suitable for tight sandstone reservoir well log interpretation.
For this purpose, domestic and foreign scholars explore tight sand lithostratigraphy factor calculation method, proposes and achieve Corresponding achievement, such as:
Mao Zhiqiang (1997), Zhang Minglu (2005), Zheng Qinglin (2006), Wu Zeyun (2009) etc. analyze different basins and cause Close sandstone HYDROCARBON-BEARING REGION porosity-resistivity relation according to porosity piecewise fitting Archie formula parameter, or is built with quadratic function Vertical porosity and lithostratigraphy factor relation, the lithostratigraphy Lai Tigao factor computational accuracy;
Tight sand resistance is considered as the parallel connection of Free water and micropore water portions of resistance by Li Xia etc. (2012), is proposed Consider the basis of dual porosity saturation computation model of pore structure;
Changxi Li et al. (2011) is established comprising porosity, permeability, is averaged by model analysis and core experiment The tight sand lithostratigraphy factor computation model of pore throat radius.
In general, above-mentioned new lithostratigraphy factor calculation method/computation model Consideration is more, and form is compared with A Er Odd formula is complicated, explains in application there is certain limitation in real logging data processing.
Therefore, for the acquisition problem of lithostratigraphy factor in the prior art, need a kind of new method of relative ease with Obtain more accurate lithostratigraphy factor.
Summary of the invention
For the acquisition problem of lithostratigraphy factor in the prior art, the present invention provides a kind of acquisition lithostratigraphy factors Method, obtain the rock according to the corresponding porosity calculation of two kinds of different pore structures of lithostratigraphy factor influence degree Stone formation factor, the described method comprises the following steps:
Archie formula improvement type, in the Archie formula improvement type, the rock are constructed according to Archie formula Formation factor porosity corresponding to described two pore structures is related;
The lithostratigraphy factor is obtained according to the Archie formula improvement type.
In one embodiment, during constructing the Archie formula improvement type, defining specific pore structure is Macro conduction hole, the corresponding porosity of the macro conductive hole is macro conductive porosity, in which:
It can convert respectively according to the different described two pore structures of special ratios as the macro conductive hole;
Specific linear relationship is presented in the macro conductive porosity and the lithostratigraphy factor under log-log coordinate.
In one embodiment, described two holes are divided according to the difference of the inflection point of the lithostratigraphy factor and porosity relationship Gap structure.
In one embodiment, described two pore structures are respectively large aperture hole and small-bore hole, corresponding hole Degree is respectively large aperture porosity and small size pores porosity, in which:
Using the pore throat radius of the corresponding hole of the inflection point as cut off value;
The hole that the pore throat radius is greater than the cut off value is divided into the large aperture hole;
The hole that the pore throat radius is less than or equal to the cut off value is divided into the small-bore hole.
In one embodiment, using porosity 10% as the inflection point.
In one embodiment, graph of pore diameter distribution is drawn to obtain different aperture pair using the pressure mercury measured data of experiment of rock sample The pore throat radius and the porosity answered.
In one embodiment, during obtaining the lithostratigraphy factor according to the Archie formula improvement type:
The parameter value for obtaining the Archie formula improvement type is calculated according to the measured data of experiment of rock sample;
The parameter value is substituted into the Archie formula improvement type and obtains the lithostratigraphy factor to calculate.
In one embodiment, during obtaining the parameter value:
The lithostratigraphy factor of the rock sample is obtained by the rock-electric test measurement data of the rock sample;
The corresponding porosity of described two pore structures is obtained by the pressure mercury measured data of experiment of the rock sample;
Lithostratigraphy factor and described two hole knots based on the Archie formula improvement type according to the rock sample The corresponding porosity fitting of structure obtains the parameter value.
In one embodiment, during calculating the acquisition lithostratigraphy factor:
The corresponding porosity of described two pore structures is obtained respectively by pressure mercury experiment measurement;
Based on the Archie formula improvement type according to the corresponding porosity calculation acquisition of described two pore structures Lithostratigraphy factor.
In one embodiment, the pseudo- macro conductive porosity of definition, the macro conduction porosity of puppet under log-log coordinate with institute It states lithostratigraphy factor and single linear relationship is presented, in the process that the parameter value is substituted into the Archie formula improvement type In:
Described in pressure mercury measured data of experiment acquisition based on the Archie formula improvement type by the relevant rock sample The correlation of the total porosity of rock sample and pseudo- macro conductive porosity;
The total porosity for obtaining rock/rock stratum to be analyzed is calculated according to the correlation obtains the rock to be analyzed The pseudo- macro conductive porosity in stone/rock stratum and then the lithostratigraphy factor for calculating the acquisition rock/rock stratum to be analyzed.
Compared with prior art, the available lithostratigraphy factor for more meeting practical geological condition of the present invention, thus The accuracy for substantially increasing geological prospecting result, the oil-gas exploration and development for after provide more accurate strong support.
Other feature or advantage of the invention will illustrate in the following description.Also, Partial Feature of the invention or Advantage will be become apparent by specification, or be appreciated that by implementing the present invention.The purpose of the present invention and part Advantage can be realized or be obtained by step specifically noted in the specification, claims and drawings.
Detailed description of the invention
Attached drawing is used to provide further understanding of the present invention, and constitutes part of specification, with reality of the invention It applies example and is used together to explain the present invention, be not construed as limiting the invention.In the accompanying drawings:
Fig. 1 is flow chart according to an embodiment of the invention;
Fig. 2 is the pseudo- macro conductive porosity of an embodiment and lithostratigraphy Factor Experiment data record figure according to the present invention;
Fig. 3 is an embodiment rock sample pore throat radius distribution map according to the present invention;
Fig. 4 is the experimental data of an embodiment and the prior art-calculating entity relationship diagram according to the present invention;
Fig. 5 is the macro conductive porosity relational graph of an embodiment porosity-puppet according to the present invention;
Fig. 6 is experimental data according to another embodiment of the present invention-calculating entity relationship diagram.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, implementation personnel of the invention whereby Can fully understand that how the invention applies technical means to solve technical problems, and reach technical effect realization process and according to The present invention is embodied according to above-mentioned realization process.As long as each embodiment it should be noted that do not constitute conflict, in the present invention And each feature in each embodiment can be combined with each other, be formed by technical solution protection scope of the present invention it It is interior.
In the prior art, determine that lithostratigraphy factor generallys use Archie formula (formula 1).Due to Archie formula It is proposed is based primarily upon the rock core sample porosity for being derived from United States Gulf area, measurements resistivity, and rock sample on probation is main It is porosity ranges between 10%~40%.In some cases, Archie formula can not be suitable for tight sand well Reservoir well log interpretation.
Meanwhile other obtain the method for lithostratigraphy factors in the prior art and the different degrees of presence of computation model is examined Worry factor is more, and form explains in application there is certain limitation etc. compared with Archie formula complexity, in real logging data processing Problem.
To solve the above problems, streamline operation on the basis of obtaining more accurate lithostratigraphy factor, this hair A kind of method of the bright acquisition lithostratigraphy factor for proposing relative ease.Method of the invention is according to lithostratigraphy factor shadow The corresponding porosity calculation of two kinds of different pore structures of the degree of sound obtains lithostratigraphy factor.It is utilized in compared to the prior art Archie formula obtains lithostratigraphy factor, and calculated result of the invention more meets practical geological condition.
Next based on the attached drawing specific implementation procedure that the present invention will be described in detail, step shown in the flow chart of attached drawing can To be executed in the computer system comprising such as a group of computer-executable instructions.Although showing each step in flow charts Logical order but in some cases can be with the steps shown or described are performed in an order that is different from the one herein.
In the present embodiment, Archie formula is carried out as shown in Figure 1, being primarily based on the research to practical geological condition To construct Archie formula improvement type, i.e. execution step S100 constructs Archie formula improvement type according to Archie formula for deformation. Archie formula improvement type by more meeting practical geological condition, which calculates, obtains lithostratigraphy factor, thus greatly improve The accuracy of lithostratigraphy factor acquisition result.Archie formula improvement type is specific in following the present invention is described in detail first Construction process.
It can be found that lithostratigraphy factor and porosity relationship usually exist obviously based on the research to practical geological condition Inflection point.Since pore structure is an important factor for another influences electric conduction of rock in addition to porosity.The present invention is directed to rock Stone formation factor and the inflection point of porosity relationship analyze the pore structure of rock.By lithostratigraphy factor and porosity Piecewise linear relationship and Different porosities under the conditions of the pore-size distribution difference of rock sample can speculate: the two of corner position correspondence Kind pore structure has notable difference to rock sample Conductivity.
Based on the studies above as a result, the present invention will have the pore structure of notable difference to distinguish rock sample Conductivity It divides, blowhole is defined as large aperture hole and two kinds of small-bore hole.Large aperture hole and small-bore hole are to rock Sample Conductivity has notable difference.Large aperture hole porosity corresponding with small-bore hole be respectively large aperture porosity and Small size pores porosity.It is constructed by the relationship established between large aperture porosity and small size pores porosity and lithostratigraphy factor Archie formula improvement type.
In Archie formula improvement type, factor porosity corresponding to two kinds of pore structures in lithostratigraphy is related.In this way The inflection point that lithostratigraphy factor and porosity relationship can be embodied in Archie formula improvement type, so that A Erqi is public Formula improvement type is more in line with practical geological condition.
For the relationship established between large aperture porosity and small size pores porosity and lithostratigraphy factor, the present invention is first Construct a kind of specific pore structure of reference that can describe large aperture porosity and small size pores porosity simultaneously.Assuming that There are a kind of pore structure of ideal strong conductive hole (being named as " macro conduction hole "), and corresponding porosity is macro conduction Porosity (is denoted as φs)。
Macro conductive porosity after converting in rock plays a decisive role to electric conduction of rock, with reference to Archie formula (formula 1), Specific linear relationship is presented in macro conduction porosity and lithostratigraphy factor under log-log coordinate.Thus it proposes with obtaining rock The new formula of layer factor:
Wherein: F is lithostratigraphy factor;a0、m0For constant related with rock property.
Large aperture hole and small-bore hole can be converted respectively according to different special ratios as macro conductive hole.It is aobvious So, small-bore hole conversion ratio (is denoted as cn) with large aperture hole conversion ratio (be denoted as cx) different, it is macro after being converted in rock Conductive porosity φsIt may be expressed as:
φs=cxφx+cnφn (3)
Wherein: φxFor large aperture porosity;φnFor small size pores porosity.
Formula (3) are substituted into formula (2), the new formula that can will acquire lithostratigraphy factor is written as:
Wherein:
Formula (4) is the Archie formula improvement type that the present invention constructs.
Next lithostratigraphy factor can be obtained using Archie formula improvement type (formula 4).Due to parameters of formula a1、 m1、c1It can be regarded as constant related with rock property.Therefore it in specifically used Archie formula improvement type, first has to obtain Take the corresponding parameter a in lithostratigraphy1、m1、c1Value.
It is calculated in the present embodiment first with the experimental data of the rock sample of lithostratigraphy to be analyzed and obtains parameters of formula a1、m1、 c1Value.Then the parameters of formula a comprising specific value is utilized1、m1、c1Archie formula improvement type obtain lithostratigraphy Factor.
Specifically, utilizing rock rock electricity-pressure mercury translocation fitting experimental data parameters of formula a1、m1、c1.S111 is first carried out, The case where determining porosity step, obtaining the porosity of rock sample.Determine the corresponding large aperture porosity of rock sample and small-bore Porosity.
In step S111, the large aperture of lithostratigraphy is first determined according to the inflection point of lithostratigraphy factor and porosity relationship The range of definition of hole and small-bore hole.Its detailed process is as follows:
Graph of pore diameter distribution is drawn using the pressure mercury measured data of experiment of rock sample to obtain the corresponding pore throat radius of different aperture And porosity;
Since there are apparent inflection points usually at porosity 10% for lithostratigraphy factor and porosity relationship, at this It with porosity 10% is inflection point in embodiment, with the pore throat radius r of the corresponding hole of inflection point0As cut off value;
Pore throat radius is greater than pore throat radius cut off value r0Hole be divided into large aperture hole, by pore throat radius be less than etc. In pore throat radius cut off value r0Be divided into small-bore hole, corresponding porosity is respectively large aperture porosity and small size pores Porosity.That is:
φx=φ × ∑ Δ Hgi (5)
φn=φ × ∑ Δ Hgj (6)
φ is rock porosity, is obtained by rock sample physical property experiment, Δ HgiIt is greater than r for pore throat radius0The corresponding mercury of pore throat Saturation degree is poor, Δ HgjIt is less than r for pore throat radius0The corresponding mercury saturation degree of pore throat it is poor.
For a certain specific rock sample, the F in formula (4) is lithostratigraphy factor, next executes step S112, is determined Lithostratigraphy factor step.Obtained in step S112 by rock-electric test the lithostratigraphy of the rock sample of lithostratigraphy to be analyzed because Element.
In formula (4), lithostratigraphy factor F is obtained by rock-electric test, φx、φnIt is surveyed by pressure mercury experimental data combination porosity Amount result obtains.Therefore step S113 next can be executed, fitting obtains fitting formula parameter a1、m1、c1
After getting parameters of formula, so that it may bring parameters of formula into further to obtain corresponding lithostratigraphy factor. During obtaining lithostratigraphy factor, need to consider whether rock sample has pressure mercury experimental data.
For there is rock sample/rock stratum of pressure mercury experimental data: first passing through pressure mercury experiment measurement and obtain two kinds of pore structures respectively Corresponding porosity;Again according to the corresponding porosity calculation acquisition rock of two kinds of pore structures based on Archie formula improvement type Layer factor.
For rock sample/rock stratum of no pressure mercury experimental data: then first passing through the pressure mercury experiment measurement number of other related rock samples According to the correlation of the total porosity for obtaining rock sample and pseudo- macro conductive porosity;The total pore space of rock stratum/rock to be analyzed is obtained again Degree;The macro conductive porosity of puppet for obtaining rock stratum/rock sample to be analyzed is finally calculated according to above-mentioned correlation, then public based on A Erqi Formula improvement type calculates lithostratigraphy factor.
For convenience of formula (4) application, the pseudo- macro conductive porosity φ of definitionps, enable
φpsx+c1φn (7)
Then the Archie formula improvement type of formula (4) is represented by
Obviously, φ in formula (8)psIt is not equal with macro conductive porosity described previously, but has certain ratio between the two Relationship.Rock electricity-pressure mercury translocation experimental data shows: under log-log coordinate, pseudo- macro conductive porosity and lithostratigraphy factor are presented Good single linear relationship.
Fig. 2 is pseudo- macro conductive porosity and lithostratigraphy Factor Experiment data record figure, and abscissa is pseudo- macro conductive hole Degree, ordinate are lithostratigraphy factor.As shown in Fig. 2, good single linear is presented in the Experiment Data Records point distribution in Fig. 2 Relationship (is evenly distributed near an oblique line), it was confirmed that the validity of formula (8).Therefore lithostratigraphy factor is next being obtained During, step S114 need to be only first carried out, pseudo- macro conductive porosity step is obtained;Step S115 is executed again, with obtaining rock Layer factor step obtains corresponding lithostratigraphy factor according to pseudo- macro conductive porosity calculation based on formula (8).
In step S114, it can be appreciated that for there is rock sample/rock stratum lithostratigraphy factor of pressure mercury experimental data, first By acquiring rock sample large aperture porosity φ respectively by pressure mercury experiment measurement (rock sample presses mercury experimental data)xAnd small size pores Porosity φn, pseudo- macro conductive porosity is then calculated by formula (7).
It is in the present embodiment, then then first first for rock sample/rock stratum lithostratigraphy factor of no pressure mercury experimental data The total porosity and puppet of rock sample are obtained by the pressure mercury measured data of experiment of other related rock samples based on Archie formula improvement type The correlation (being equivalent to the correlation for obtaining total porosity and lithostratigraphy factor) of macro conduction porosity;It is obtained again wait divide Analyse rock sample/rock stratum total porosity;The macro conductive hole of puppet for obtaining rock stratum/rock sample to be analyzed is finally calculated according to above-mentioned correlation Porosity.
Illustrate specific implementation procedure of the invention with a specific example application.
Lithostratigraphy factor is obtained to certain oil field tight sand rock sample using the present invention.Rock sample is carried out first Mercury experiment is pressed, draws graph of pore diameter distribution, as shown in Figure 3.Abscissa is pore throat radius (unit: micron μm) in Fig. 3, and ordinate is Mercury saturation degree is poor, and the point on dotted line represents hole of the porosity less than 10%, and the point on solid line represents porosity and is greater than 10% Hole.As seen from Figure 3, rock sample pore throat radius of the porosity less than 10% is approximately less than 0.3 μm, and porosity is greater than 10% Rock sample pore throat radius be greater than 0.3 μm based on.It is thus determined that porosity is greater than 10% and aperture of the porosity less than 10% point Dividing value r0=0.3 μm.
Next the large aperture porosity φ of each rock sample is determinedxWith small size pores porosity φn, measured in conjunction with rock-electric test Lithostratigraphy factor F is fitted to obtain formula (4) parameter a1=5.23, m1=0.867, c1=0.280, the area 43 is determined using formula (7) Block has the pseudo- macro conductive porosity of the core sample of pressure mercury experimental data, and formula (8) is recycled to determine lithostratigraphy factor.
The rock sample lithostratigraphy factor that Binding experiment obtains draws experimental data-calculating entity relationship diagram.As shown in figure 4, Abscissa represents the lithostratigraphy factor that experiment measurement obtains, and ordinate represents the lithostratigraphy factor for calculating and obtaining.In Fig. 4 In, the ordinate of solid data points represents the lithostratigraphy factor for calculating acquisition according to the method for the present invention, hollow data points Ordinate, which is represented, calculates the lithostratigraphy factor obtained according to the prior art (Archie formula).
As seen from Figure 4, solid data points are distributed in oblique line nearby (oblique line are 45 ° of lines, and the point transverse and longitudinal on oblique line is sat Mark equal), the lithostratigraphy factor that display calculates acquisition according to the present invention is coincide substantially with experimental rock formation factor, and (calculating obtains The lithostratigraphy factor and experimental rock formation factor taken is essentially identical).
Meanwhile compared to hollow data points, the distribution of solid data points is more close to oblique line (most of solid data points collection In near oblique line, compared to solid data points, the distribution of hollow data points it is more dispersed and apart from oblique line farther out).This explanation Compared with the prior art calculates the lithostratigraphy factor obtained by Archie formula, according to the lithostratigraphy of this calculatings acquisition because Element is more in line with actual conditions.Compared to the prior art, method of the invention shows preferably to calculate effect.
To 72 pieces of the area without pressing the core sample of mercury experimental data, closed first according to porosity and pseudo- macro conductive porosity System determines pseudo- macro conductive porosity.Herein, the macro conductive porosity relational graph of porosity-puppet is drawn first, as shown in figure 5, horizontal Coordinate is porosity, and ordinate is pseudo- macro conductive porosity.Pressure mercury experimental data based on other samples is drawn in relational graph A series of data point determines the relationship between porosity and pseudo- macro conductive porosity according to the distribution of data point.Such as institute in Fig. 5 Show, the oblique line (void in figure of the relationship between description porosity and pseudo- macro conductive porosity is determined according to the distribution of data point Line).Trend (slope) based on the dotted line in figure, which calculates, obtains pseudo- macro conductive porosity.
Finally lithostratigraphy factor is calculated using formula (8).In conjunction with the lithostratigraphy factor and calculating of rock sample experiment measurement The lithostratigraphy factor of acquisition draws experimental data-calculating entity relationship diagram.As shown in fig. 6, abscissa is that experiment measurement obtains Lithostratigraphy factor, ordinate be according to the present invention calculate obtain lithostratigraphy factor.It is attached that data point is distributed in oblique line Closely (oblique line is 45 ° of lines, and the point transverse and longitudinal coordinate on oblique line is equal), indicate to calculate according to the present invention the lithostratigraphy factor of acquisition with Experimental rock formation factor coincide substantially and (calculates the lithostratigraphy factor obtained and experimental rock formation factor is essentially identical).
To sum up, compared with prior art, the available lithostratigraphy factor for more meeting practical geological condition of the present invention, To substantially increase the accuracy of geological prospecting result, the oil-gas exploration and development for after provides more accurate strong branch It holds.
While it is disclosed that embodiment content as above but described only to facilitate understanding the present invention and adopting Embodiment is not intended to limit the invention.Method of the present invention can also have other various embodiments.Without departing substantially from In the case where essence of the present invention, those skilled in the art make various corresponding changes or change in accordance with the present invention Shape, but these corresponding changes or deformation all should belong to scope of protection of the claims of the invention.

Claims (14)

1. a kind of method for obtaining lithostratigraphy factor, which is characterized in that according to different to lithostratigraphy factor influence degree The corresponding porosity calculation of two kinds of pore structures obtains the lithostratigraphy factor, the described method comprises the following steps:
Archie formula improvement type, in the Archie formula improvement type, the lithostratigraphy are constructed according to Archie formula Factor porosity corresponding to described two pore structures is related;
The lithostratigraphy factor is obtained according to the Archie formula improvement type;
Wherein, during constructing the Archie formula improvement type, defmacro conduction hole, the macro conductive hole is corresponded to Porosity be macro conductive porosity, in which:
Described two pore structures are converted respectively as the macro conductive hole, following expression is based on, to obtain described macro lead Electric porosity,
φs=cxφx+cnφn
Wherein,Indicate macro conductive porosity, φx、φnRespectively correspond the porosity for indicating two kinds of pore structures, cx、cnIt is right respectively It should indicate the conversion ratio of two kinds of pore structures;
Single linear relationship is presented in the macro conductive porosity and the lithostratigraphy factor under log-log coordinate.
2. the method according to claim 1, wherein turning according to the lithostratigraphy factor and porosity relationship Point difference divides described two pore structures.
3. according to the method described in claim 2, it is characterized in that, described two pore structures are respectively large aperture hole and small Aperture hole, corresponding porosity are respectively large aperture porosity and small size pores porosity, in which:
Using the pore throat radius of the corresponding hole of the inflection point as cut off value;
The hole that the pore throat radius is greater than the cut off value is divided into the large aperture hole;
The hole that the pore throat radius is less than or equal to the cut off value is divided into the small-bore hole.
4. according to the method described in claim 3, it is characterized in that, the pressure mercury measured data of experiment using rock sample draws aperture point Butut is to obtain the corresponding pore throat radius of different aperture and the porosity.
5. the method according to any one of claim 2-4, which is characterized in that using porosity 10% as the inflection point.
6. method according to any of claims 1-4, which is characterized in that according to the Archie formula improvement type During obtaining the lithostratigraphy factor:
The parameter value for obtaining the Archie formula improvement type is calculated according to the measured data of experiment of rock sample;
The parameter value is substituted into the Archie formula improvement type and obtains the lithostratigraphy factor to calculate.
7. according to the method described in claim 5, it is characterized in that, obtaining the rock according to the Archie formula improvement type During stone formation factor:
The parameter value for obtaining the Archie formula improvement type is calculated according to the measured data of experiment of rock sample;
The parameter value is substituted into the Archie formula improvement type and obtains the lithostratigraphy factor to calculate.
8. according to the method described in claim 6, it is characterized in that, during obtaining the parameter value:
The lithostratigraphy factor of the rock sample is obtained by the rock-electric test measurement data of the rock sample;
The corresponding porosity of described two pore structures is obtained by the pressure mercury measured data of experiment of the rock sample;
Lithostratigraphy factor and described two pore structures pair based on the Archie formula improvement type according to the rock sample The porosity fitting answered obtains the parameter value.
9. the method according to the description of claim 7 is characterized in that during obtaining the parameter value:
The lithostratigraphy factor of the rock sample is obtained by the rock-electric test measurement data of the rock sample;
The corresponding porosity of described two pore structures is obtained by the pressure mercury measured data of experiment of the rock sample;
Lithostratigraphy factor and described two pore structures pair based on the Archie formula improvement type according to the rock sample The porosity fitting answered obtains the parameter value.
10. according to the method described in claim 6, it is characterized in that, during calculating the acquisition lithostratigraphy factor:
The corresponding porosity of described two pore structures is obtained respectively by pressure mercury experiment measurement;
The rock is obtained according to the corresponding porosity calculation of described two pore structures based on the Archie formula improvement type Formation factor.
11. the method according to any one of claim 7-9, which is characterized in that calculate obtain the lithostratigraphy because In the process of element:
The corresponding porosity of described two pore structures is obtained respectively by pressure mercury experiment measurement;
The rock is obtained according to the corresponding porosity calculation of described two pore structures based on the Archie formula improvement type Formation factor.
12. according to the method described in claim 6, it is characterized in that, the pseudo- macro conductive porosity of definition, the macro conductive hole of the puppet Single linear relationship is presented with the lithostratigraphy factor under log-log coordinate in degree, and the parameter value is being substituted into the A Er During odd formula improvement type:
The rock sample is obtained by the pressure mercury measured data of experiment of the relevant rock sample based on the Archie formula improvement type Total porosity and pseudo- macro conductive porosity correlation;
Obtain the total porosity of rock/rock stratum to be analyzed calculated according to the correlation obtain the rock to be analyzed/ The macro conductive porosity of the puppet of rock stratum and then the lithostratigraphy factor for calculating the acquisition rock/rock stratum to be analyzed.
13. according to the method for claim 11, which is characterized in that the pseudo- macro conductive porosity of definition, the macro conductive hole of puppet Single linear relationship is presented with the lithostratigraphy factor under log-log coordinate in porosity, by the parameter value substitute into Ah During your odd formula improvement type:
The rock sample is obtained by the pressure mercury measured data of experiment of the relevant rock sample based on the Archie formula improvement type Total porosity and pseudo- macro conductive porosity correlation;
Obtain the total porosity of rock/rock stratum to be analyzed calculated according to the correlation obtain the rock to be analyzed/ The macro conductive porosity of the puppet of rock stratum and then the lithostratigraphy factor for calculating the acquisition rock/rock stratum to be analyzed.
14. the method according to any one of claim 7-10, which is characterized in that the pseudo- macro conductive porosity of definition, it is described Single linear relationship is presented with the lithostratigraphy factor under log-log coordinate in pseudo- macro conductive porosity, by the parameter value During substituting into the Archie formula improvement type:
The rock sample is obtained by the pressure mercury measured data of experiment of the relevant rock sample based on the Archie formula improvement type Total porosity and pseudo- macro conductive porosity correlation;
Obtain the total porosity of rock/rock stratum to be analyzed calculated according to the correlation obtain the rock to be analyzed/ The macro conductive porosity of the puppet of rock stratum and then the lithostratigraphy factor for calculating the acquisition rock/rock stratum to be analyzed.
CN201510388742.0A 2015-07-06 2015-07-06 A method of obtaining lithostratigraphy factor Active CN106321087B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510388742.0A CN106321087B (en) 2015-07-06 2015-07-06 A method of obtaining lithostratigraphy factor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510388742.0A CN106321087B (en) 2015-07-06 2015-07-06 A method of obtaining lithostratigraphy factor

Publications (2)

Publication Number Publication Date
CN106321087A CN106321087A (en) 2017-01-11
CN106321087B true CN106321087B (en) 2019-05-07

Family

ID=57727224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510388742.0A Active CN106321087B (en) 2015-07-06 2015-07-06 A method of obtaining lithostratigraphy factor

Country Status (1)

Country Link
CN (1) CN106321087B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107366539B (en) * 2017-07-18 2020-05-22 中国石油大学(华东) Method for predicting porosity of sandstone reservoir in mountain front area
CN111595752B (en) * 2020-06-22 2021-07-06 中国科学技术大学 Method for determining effective porosity of rock
CN112539057B (en) * 2021-01-21 2023-11-28 陕西延长石油(集团)有限责任公司 Template and method for identifying fluid properties of tight sandstone reservoir

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102979517A (en) * 2012-12-04 2013-03-20 中国海洋石油总公司 Method for quantitatively evaluating saturation of complex oil and gas reservoir
CN103225506A (en) * 2013-04-19 2013-07-31 中国石油大学(华东) Method for establishing three-part automatic parallel-series electric conduction saturation degree model
CN103912271A (en) * 2014-04-25 2014-07-09 中国石油大学(华东) Method for grading evaluation of tight sandstone gas resource

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8061442B2 (en) * 2008-07-07 2011-11-22 Bp Corporation North America Inc. Method to detect formation pore pressure from resistivity measurements ahead of the bit during drilling of a well
WO2015017190A2 (en) * 2013-08-01 2015-02-05 Conocophillips Company Dynamic in-situ measurement of reservoir wettability

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102979517A (en) * 2012-12-04 2013-03-20 中国海洋石油总公司 Method for quantitatively evaluating saturation of complex oil and gas reservoir
CN103225506A (en) * 2013-04-19 2013-07-31 中国石油大学(华东) Method for establishing three-part automatic parallel-series electric conduction saturation degree model
CN103912271A (en) * 2014-04-25 2014-07-09 中国石油大学(华东) Method for grading evaluation of tight sandstone gas resource

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
完全含水多孔岩石电学性质及其孔隙结构实验研究;毛志强等;《石油学报》;19970805;第18卷(第3期);第51-55页

Also Published As

Publication number Publication date
CN106321087A (en) 2017-01-11

Similar Documents

Publication Publication Date Title
RU2523776C2 (en) Method for quantitative calculation of saturation of fractured reservoir with hydrocarbons
CN103422853B (en) A kind of method and device determining stratum water saturation
CN101725344B (en) Method for determining litho-electric parameters
CN103353462B (en) A kind of rock nonuniformity method for quantitatively evaluating based on Magnetic resonance imaging
Clerke Permeability, Relative Permeability, Microscopic Displacement Efficiency and Pore Geometry of M_1 Bimodal Pore Systems in Arab D Limestone
CN102175832A (en) Method for determining optimal saturation calculation model of typical reservoir
CN103884633A (en) Method and device for determining rock permeability
Wang et al. The gravity matching area selection criteria for underwater gravity-aided navigation application based on the comprehensive characteristic parameter
CN105223116A (en) A kind of method calculating irreducible water saturation based on nuclear magnetic resoance spectrum Y-factor method Y
CN106321087B (en) A method of obtaining lithostratigraphy factor
CN108049866A (en) Two-dimensional NMR Logging tight gas reservoir quantitative evaluation method
Zou et al. Porosity type analysis and permeability model for micro-trans-pores, meso-macro-pores and cleats of coal samples
CN109164033B (en) Method for acquiring hydraulic conductivity of unsaturated soil based on resistivity method
CN104863574A (en) Fluid identification method applicable to compact sandstone reservoir
CN105676309A (en) Fractured reservoir reverse imbibition recovery prediction method
CN109387468A (en) Shale reservoir nanoaperture structural characteristic parameter method for testing and analyzing and system
Liu et al. Permeability prediction from mercury injection capillary pressure curves by partial least squares regression method in tight sandstone reservoirs
CN109580689A (en) A kind of node-by-node algorithm method of nuclear magnetic resonance log T2 cutoff value
Zhou et al. Self-constrained Euler deconvolution using potential field data of different altitudes
CN104374679B (en) A kind of method for determining reservoir permeability contribution margin
CN109856688B (en) Fluid property identification method based on nuclear magnetic logging double TW polarization enhancement method
CN113640890B (en) Reservoir gas-water relative permeability curve construction method and system based on nuclear magnetic resonance logging
Khalili et al. Formation factor for heterogeneous carbonate rocks using multi-scale Xray-CT images
CN108444895B (en) Efficient method for obtaining unsaturated permeability parameters of loess body
CN104514553A (en) Recognition method of gas reservoir in lithic sandstone layer and application thereof

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant