CN106094053A - Hydrocarbon source rock hydrocarbon generation and discharge strength calculation method and device - Google Patents
Hydrocarbon source rock hydrocarbon generation and discharge strength calculation method and device Download PDFInfo
- Publication number
- CN106094053A CN106094053A CN201610395936.8A CN201610395936A CN106094053A CN 106094053 A CN106094053 A CN 106094053A CN 201610395936 A CN201610395936 A CN 201610395936A CN 106094053 A CN106094053 A CN 106094053A
- Authority
- CN
- China
- Prior art keywords
- hydrocarbon
- resistivity
- natural gamma
- lithology
- transit time
- 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
Links
- 229930195733 hydrocarbon Natural products 0.000 title claims abstract description 185
- 150000002430 hydrocarbons Chemical class 0.000 title claims abstract description 185
- 239000004215 Carbon black (E152) Substances 0.000 title claims abstract description 184
- 239000011435 rock Substances 0.000 title claims abstract description 90
- 238000004364 calculation method Methods 0.000 title claims abstract description 38
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 8
- 229910052799 carbon Inorganic materials 0.000 claims description 8
- 238000000034 method Methods 0.000 abstract description 6
- 238000011156 evaluation Methods 0.000 description 11
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 8
- 241000208340 Araliaceae Species 0.000 description 5
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 5
- 235000003140 Panax quinquefolius Nutrition 0.000 description 5
- 235000008434 ginseng Nutrition 0.000 description 5
- 239000004575 stone Substances 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method and a device for calculating hydrocarbon generation and discharge strength of a hydrocarbon source rock, wherein the calculating method comprises the following steps: acquiring logging data of a single well, wherein the logging data at least comprises: acoustic time difference, natural gamma, resistivity and depth scale; identifying lithology in the single well according to the acoustic time difference, the natural gamma, the resistivity, the depth scale and respective parameter distinguishing intervals; substituting the acoustic time difference, the natural gamma and the resistivity into a hydrocarbon generation intensity calculation formula corresponding to the lithology to calculate the hydrocarbon generation intensity of the lithology; and substituting the hydrocarbon generation strength into the hydrocarbon generation and expulsion strength relational expression to calculate the hydrocarbon expulsion strength of the lithology. By implementing the method, various parameter factors are considered, the obtained calculation result is more accurate, the actual geological environment of the logging is better met, and more accurate data support is provided for subsequent logging work.
Description
Technical field
The present invention is about oil reservoir logging technology, specifically, be about a kind of hydrocarbon source rock life residence strength calculation method and
Device.
Background technology
Along with the development of petroleum exploration in China technology, Exploration Domain constantly expands, by conventional clastic oil-gas pools to
Carbonate rock hydrocarbon reservoir develops, and wherein Carbonate Source Rocks evaluation becomes an important content.In Carbonate Source Rocks,
There are different kinds of rocks, they have different hydrocarbon generation capacities, though and same lithology, time environment during owing to being formed
Difference, the hydrocarbon potentiality being had are the most different.So exact evaluation hydrocarbon source rock becomes in the emphasis of carbonate rock RESERVE EVALUATION
Hold.
At petroleum exploration field, Evaluation of source rocks is the primary work evaluating a block with or without Exploration Potential, and raw row
Hydrocarbon intensity evaluation is the core of Evaluation of source rocks.Multiple oil fields block degree of prospecting now of China is the most relatively low, to some
The hydrocarbon source rock position of new block, life residence evaluation is the primary work of exploration, and the old block higher to some degrees of prospecting
For new hydrocarbon source rock position, carry out life residence evaluation, find that new hydrocarbon layer position is also to increase the important channel of old filed yield.
For hydrocarbon source rock position, the evaluation of generation and expulsion intension is the most important thing, and the evaluation calculation of hydrocarbon intensity can clear and definite exploration strata
The hydrocarbon potentiality of position, row's hydrocarbon intensity is then to judging the discharge of hydro carbons and becoming Tibetan to have great importance.But, existing Logging Evaluation
Generation and expulsion intension can not be given evaluate accurately in conjunction with the multiple different factor in calibrating by valency means.
Summary of the invention
The main purpose of the embodiment of the present invention is to provide a kind of hydrocarbon source rock life residence strength calculation method and device, with more
Add the generation and expulsion intension calculating well logging exactly, and result of calculation more conforms to the actual environment of well logging.
To achieve these goals, the embodiment of the present invention provides a kind of hydrocarbon source rock life residence strength calculation method, described
Hydrocarbon source rock life residence strength calculation method includes: obtaining the log data of individual well, described log data at least includes: during sound wave
Difference, natural gamma, resistivity and deep calibration；According to described interval transit time, natural gamma, resistivity and deep calibration with respective
Parameter differentiate that interval identifies the lithology in described individual well；By described interval transit time normalizing parameter, natural gamma standardization ginseng
Number, resistivity normalizing parameter substitute into hydrocarbon strength calculation formula corresponding to described lithology and calculate the hydrocarbon intensity of described lithology；
Described hydrocarbon intensity is substituted into generation and expulsion intension relational expression and calculates row's hydrocarbon intensity of described lithology.
In one embodiment, above-mentioned according to described interval transit time, natural gamma, resistivity and deep calibration with respective
Parameter differentiates that interval identifies the lithology in described individual well, specifically includes: interval transit time, natural gamma and the resistivity described in statistics
Fall the frequency in respective demarcation interval, obtains the probability distribution of described interval transit time, natural gamma and resistivity respectively；Root
Interval transit time standardization is determined with the probability distribution of interval transit time, natural gamma and resistivity in key hole according to described probability distribution
Parameter, natural gamma normalizing parameter and resistivity normalizing parameter；According to described interval transit time normalizing parameter, natural gamma
Normalizing parameter, resistivity normalizing parameter and described deep calibration and respective parameter differentiate that interval identifies in described individual well
Lithology.
In one embodiment, above-mentioned according to described probability distribution and interval transit time, natural gamma and resistance in key hole
The probability distribution of rate determines interval transit time normalizing parameter, natural gamma normalizing parameter and resistivity normalizing parameter, specifically
Including: respectively by the interval transit time of described individual well, natural gamma, resistivity probability distribution in numerical value section residing for peak value, with
The interval transit time of described key hole, natural gamma, resistivity probability distribution in the numerical value section residing for peak value difference definition
For interval transit time, natural gamma, the meter constant of resistivity；By the interval transit time of described individual well, natural gamma, resistivity respectively
Plus respective meter constant, obtain described interval transit time normalizing parameter, natural gamma normalizing parameter and resistivity mark
Standardization parameter.
In one embodiment, above-mentioned hydrocarbon strength calculation formula is:Its
In, H is the thickness of hydrocarbon source rock；ρ is the rock density in described lithology；KGasFor the hydrocarbon rate in described lithology；TOCiFor each
The content of organic carbon in unit mass rock in lithology described in the TOC demarcation interval preset；QGas iAccording to each TOCiRespectively
Carry out calculating the hydrocarbon intensity obtained；I=1,2 ..., j, j are positive integer.
In one embodiment, above-mentioned generation and expulsion intension relational expression is: row's hydrocarbon intensity=hydrocarbon intensity × row's hydrocarbon rate.
The embodiment of the present invention also provides for a kind of hydrocarbon source rock life residence Strength co-mputation device, described hydrocarbon source rock generation and expulsion intension
Calculating device includes: log data acquiring unit, and for obtaining the log data of individual well, described log data at least includes: sound
The ripple time difference, natural gamma, resistivity and deep calibration；Lithology Discrimination unit, for according to described interval transit time, natural gamma,
With respective parameter, resistivity and deep calibration differentiate that interval identifies the lithology in described individual well；Hydrocarbon Strength co-mputation unit, uses
In described interval transit time normalizing parameter, natural gamma normalizing parameter, resistivity normalizing parameter are substituted into described lithology pair
The hydrocarbon strength calculation formula answered calculates the hydrocarbon intensity of described lithology；Row's hydrocarbon Strength co-mputation unit, for by strong for described hydrocarbon
Degree substitutes into generation and expulsion intension relational expression and calculates row's hydrocarbon intensity of described lithology.
In one embodiment, above-mentioned Lithology Discrimination unit specifically for: interval transit time described in statistics, natural gamma and
Resistivity falls the frequency in respective demarcation interval, obtains the probability of described interval transit time, natural gamma and resistivity respectively
Distribution；According to when in described probability distribution and key hole, the probability distribution of interval transit time, natural gamma and resistivity determines sound wave
Difference normalizing parameter, natural gamma normalizing parameter and resistivity normalizing parameter；According to described interval transit time normalizing parameter,
It is described that natural gamma normalizing parameter, resistivity normalizing parameter and described deep calibration and respective parameter differentiate that interval identifies
Lithology in individual well.
In one embodiment, above-mentioned according to described probability distribution and interval transit time, natural gamma and resistance in key hole
The probability distribution of rate determines interval transit time normalizing parameter, natural gamma normalizing parameter and resistivity normalizing parameter, specifically
Including: respectively by the interval transit time of described individual well, natural gamma, resistivity probability distribution in numerical value section residing for peak value, with
The interval transit time of described key hole, natural gamma, resistivity probability distribution in the numerical value section residing for peak value difference definition
For interval transit time, natural gamma, the meter constant of resistivity；By the interval transit time of described individual well, natural gamma, resistivity respectively
Plus respective meter constant, obtain described interval transit time normalizing parameter, natural gamma normalizing parameter and resistivity mark
Standardization parameter.
In one embodiment, above-mentioned hydrocarbon strength calculation formula is:Its
In, H is the thickness of hydrocarbon source rock；ρ is the rock density in described lithology；KGasFor the hydrocarbon rate in described lithology；TOCiFor each
The content of organic carbon in unit mass rock in lithology described in the TOC demarcation interval preset；QGas iAccording to each TOCiRespectively
Carry out calculating the hydrocarbon intensity obtained；I=1,2 ..., j, j are positive integer.
In one embodiment, above-mentioned generation and expulsion intension relational expression is: row's hydrocarbon intensity=hydrocarbon intensity × row's hydrocarbon rate.
Having the beneficial effects that of the embodiment of the present invention, it is possible to study area hydrocarbon source rock is carried out Lithology Discrimination accurately, passes through
To TOC and the statistical computation of hydrocarbon source rock thickness, obtain hydrocarbon primary rock producing hydrocarbon intensity, row's hydrocarbon in conjunction with parameters such as lithology and density logging, hydrocarbon rates
Intensity.Because considering many kinds of parameters factor, calculated results is more accurate, and more conforms to residing actual geology ring of logging well
Border, provides more accurate data support for follow-up well logging.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, embodiment will be described below
The accompanying drawing used required in is briefly described, it should be apparent that, the accompanying drawing in describing below is only some of the present invention
Embodiment, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to according to these
Accompanying drawing obtains other accompanying drawing.
Fig. 1 is the flow chart of the hydrocarbon source rock life residence strength calculation method according to the embodiment of the present invention；
Fig. 2 A to Fig. 2 C is the probability distribution curve figure of the log data according to the embodiment of the present invention；
Fig. 3 is the structural representation of the hydrocarbon source rock life residence Strength co-mputation device according to the embodiment of the present invention.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Describe, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments wholely.Based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under not making creative work premise
Embodiment, broadly falls into the scope of protection of the invention.
The embodiment of the present invention provides a kind of hydrocarbon source rock life residence strength calculation method and device.Below in conjunction with accompanying drawing to this
Bright it is described in detail.
The embodiment of the present invention provides a kind of hydrocarbon source rock life residence strength calculation method, as it is shown in figure 1, this hydrocarbon source rock life residence
Strength calculation method specifically includes that
Step S101: obtaining the log data of individual well, log data at least includes: interval transit time, natural gamma, resistivity
And deep calibration；
Step S102: differentiate that interval is known with respective parameter according to interval transit time, natural gamma, resistivity, deep calibration
Lithology in other individual well；
Step S103: interval transit time, natural gamma, resistivity are substituted into hydrocarbon strength calculation formula corresponding to lithology and calculates
The hydrocarbon intensity of lithology；
Step S104: hydrocarbon intensity is substituted into generation and expulsion intension relational expression and calculates row's hydrocarbon intensity of lithology.
By above-mentioned step S101 to step S104, the hydrocarbon source rock life residence strength calculation method of the embodiment of the present invention,
By oil field conventional logging data are analyzed, select to carry out interval division by parameter, differentiate boring chance lithology, logical
Cross the statistical computation to TOC and hydrocarbon source rock thickness combine the parameter such as lithology and density logging, hydrocarbon rate obtain hydrocarbon primary rock producing hydrocarbon intensity, row hydrocarbon
Intensity, the result of calculation obtained is more accurate, and owing to considering multiple logging environment factor, acquired results also more conforms to survey
Actual environment residing for well, provides significantly more efficient data support for follow-up well logging of carrying out.
Below in conjunction with concrete example, to each step in the hydrocarbon source rock life residence strength calculation method of the embodiment of the present invention
Elaborate.
Above-mentioned step S101, obtains the log data of individual well, and log data at least includes: interval transit time (AC), hole diameter
(CAL), natural gamma (GR), resistivity (Rt) and deep calibration (DEPT).
Above-mentioned steps S102, differentiates district according to interval transit time, natural gamma, resistivity, deep calibration with respective parameter
Between identify the lithology in individual well.
Specifically, first add up interval transit time, natural gamma and resistivity and fall the frequency in respective demarcation interval, point
Do not obtain the probability distribution of interval transit time, natural gamma and resistivity.Wherein, the maximum magnitude of the demarcation interval of parameters is
The numerical range of corresponding every a bite this parameter value of well, and minizone graduate is according between this parameter maximum, minima
The ratio of difference and total frequency determine.Then filter out the interval that frequency occurs, and add up its interior frequency value, generate general
Rate scattergram, as shown in Fig. 2 A to Fig. 2 C.
After the probability distribution generating the interval transit time of well logging, natural gamma and resistivity, can according to probability distribution with
In key hole, the probability distribution of interval transit time, natural gamma and resistivity determines interval transit time normalizing parameter, natural gamma mark
Standardization parameter and resistivity normalizing parameter.
Key hole is the well being positioned at hydrocarbon primary rock producing hydrocarbon down warping region, and the parameter of this key hole can obtain from oil field prospecting institute.School
Normal number is obtained by the probability distribution curve contrast of well logging with key hole.Due to the particularity of Marine Carbonate Rocks hydrocarbon source rock,
The affinity relevant to shale content makes selected interval transit time, natural gamma and resistivity these three parameter divide at peak value
Should have concordance, i.e. log parameter peak value on cloth and should occur in same numerical value section, therefore, by the probability of well logging with key hole
In distribution curve, the difference of the numerical value section at two peak value places is defined as meter constant, wherein, when different peak value occurs, be with
Maximum is as the peak value of final foundation.And give, by this meter constant, full well section of logging well, by its all well logging numerical value plus being somebody's turn to do
Meter constant so that two well peak values are in same numerical value section, i.e. completes standardization correction, obtains interval transit time standardization ginseng
Number, natural gamma normalizing parameter and resistivity normalizing parameter.
The differentiation of lithology is according to interval transit time (AC), hole diameter (CAL), natural gamma (GR), resistivity (Rt) and the degree of depth
Four parameters of scale (DEPT) are carried out, such as when 100 < AC < 210 and 2 < GR < 15,130 < RT < when 350, it is determined that this lithology is
Mud stone, the rest may be inferred, and the lithology such as carbonate rock and transition lithology has corresponding parameter to differentiate interval, according to these several rocks respectively
Property respective parameter differentiate interval carry out lithology discriminant.
Above-mentioned step S103, substitutes into, by interval transit time, natural gamma, resistivity, the hydrocarbon Strength co-mputation public affairs that lithology is corresponding
Formula calculates the hydrocarbon intensity of lithology.
Specifically, it is that above-mentioned interval transit time, natural gamma, resistivity are substituted into formula
In, calculate the hydrocarbon intensity of each lithology respectively, wherein, H is the thickness of hydrocarbon source rock；ρ is the rock density in described lithology；KGasFor
Hydrocarbon rate in described lithology；TOC is the content of organic carbon in unit mass rock in described lithology.TOCiPreset for each
The content of organic carbon in unit mass rock in lithology described in TOC demarcation interval；QGas iAccording to each TOCiCount respectively
Calculate the hydrocarbon intensity obtained；I=1,2 ..., j, j are positive integer.Such as, the TOC=0.001316GR+0.00228AC-of mud stone
0.056148lgRT-0.30755, correspondingly, various lithology have each self-corresponding formula to calculate its TOC, the ginseng used
Number is identical, adjusts the coefficient before each parameter according to different lithology, numerous to list herein.
TOC is partitioned into different TOC demarcation intervals (e.g. according to 0～0.3,0.3 by the value of calculation according to described TOC
～0.6,0.6～1, > 1 four TOC demarcation intervals divide) calculate described hydrocarbon intensity respectively, then each is divided
Interval corresponding hydrocarbon intensity sums up calculating.In actual applications, follow-up can be corresponding according to above-mentioned TOC demarcation interval
Hydrocarbon intensity carries out area-method calculating, obtains the stock number of target zone, compares and does not divides the money that more can accurately calculate target zone
Source is measured.
After calculating the hydrocarbon intensity obtaining each lithology, above-mentioned steps S106 can be passed through, hydrocarbon intensity is updated to raw
Row's hydrocarbon strength relationships calculates row's hydrocarbon intensity of each lithology.Specifically, this generation and expulsion intension relational expression is: row's hydrocarbon intensity=life
Hydrocarbon intensity × row's hydrocarbon rate, this row's hydrocarbon rate can obtain row's hydrocarbon rate of various lithology by experiment.
The hydrocarbon source rock life residence strength calculation method of the embodiment of the present invention, it is possible to study area hydrocarbon source rock is carried out rock accurately
Property identification, by TOC and the statistical computation of hydrocarbon source rock thickness, obtains primary rock producing hydrocarbon in conjunction with parameters such as lithology and density logging, hydrocarbon rates
Hydrocarbon intensity, row's hydrocarbon intensity.Because considering many kinds of parameters factor, calculated results is more accurate, and it is residing to more conform to well logging
Actual geological environment, provide more accurate data support for follow-up well logging.
The embodiment of the present invention provides a kind of hydrocarbon source rock life residence Strength co-mputation device, as it is shown on figure 3, this hydrocarbon source rock life residence
Strength co-mputation device specifically includes that log data acquiring unit 1, Lithology Discrimination unit 2, hydrocarbon Strength co-mputation unit 3 and row's hydrocarbon
Strength co-mputation unit 4.
Wherein, above-mentioned log data acquiring unit 1 is for obtaining the log data of individual well, and log data at least includes:
Interval transit time, natural gamma, resistivity and deep calibration；Lithology Discrimination unit 2 is for according to interval transit time, natural gamma, electricity
Resistance rate, deep calibration and respective parameter differentiate that interval identifies the lithology in individual well；Hydrocarbon Strength co-mputation unit 3 is for by sound wave
The time difference, natural gamma, resistivity substitute into hydrocarbon strength calculation formula corresponding to lithology and calculate the hydrocarbon intensity of lithology；Row's hydrocarbon intensity
Computing unit 4 calculates row's hydrocarbon intensity of lithology for hydrocarbon intensity is substituted into generation and expulsion intension relational expression.
By each above-mentioned ingredient, the hydrocarbon source rock life residence Strength co-mputation device of the embodiment of the present invention, by right
Oil field conventional logging data are analyzed, and select to carry out interval division by parameter, differentiate, by right boring chance lithology
The statistical computation of TOC and hydrocarbon source rock thickness combine the parameter such as lithology and density logging, hydrocarbon rate obtain hydrocarbon primary rock producing hydrocarbon intensity, row hydrocarbon strong
Degree, the result of calculation obtained is more accurate, and owing to considering multiple logging environment factor, acquired results also more conforms to well logging
Residing actual environment, provides significantly more efficient data support for follow-up well logging of carrying out.
Below in conjunction with concrete example, to each composition in the hydrocarbon source rock life residence Strength co-mputation device of the embodiment of the present invention
Part elaborates.
Above-mentioned log data acquiring unit 1, for obtaining the log data of individual well, log data at least includes: sound wave
The time difference (AC), hole diameter (CAL), natural gamma (GR), resistivity (Rt) and deep calibration (DEPT).
Above-mentioned Lithology Discrimination unit 2, is used for according to interval transit time, natural gamma, resistivity, deep calibration with respective
Parameter differentiates that interval identifies the lithology in individual well.
Specifically, first this Lithology Discrimination unit 2 is added up interval transit time, natural gamma and resistivity and is fallen in respective division
Frequency in interval, obtains the probability distribution of interval transit time, natural gamma and resistivity respectively.Wherein, the division of parameters
Interval maximum magnitude is the numerical range of corresponding every a bite this parameter value of well, and minizone graduate is according to this parameter
The ratio of big difference between value, minima and total frequency determines.Then filter out the interval that frequency occurs, and add up in it
Frequency value, generating probability scattergram, as shown in Fig. 2 A to Fig. 2 C.
After the probability distribution generating the interval transit time of well logging, natural gamma and resistivity, can according to probability distribution with
In key hole, the probability distribution of interval transit time, natural gamma and resistivity determines interval transit time normalizing parameter, natural gamma mark
Standardization parameter and resistivity normalizing parameter.
Key hole is the well being positioned at hydrocarbon primary rock producing hydrocarbon down warping region, and the parameter of this key hole can obtain from oil field prospecting institute.School
Normal number is obtained by the probability distribution curve contrast of well logging with key hole.Due to the particularity of Marine Carbonate Rocks hydrocarbon source rock,
The affinity relevant to shale content makes selected interval transit time, natural gamma and resistivity these three parameter divide at peak value
Should have concordance, i.e. log parameter peak value on cloth and should occur in same numerical value section, therefore, by the probability of well logging with key hole
In distribution curve, the difference of the numerical value section at two peak value places is defined as meter constant, wherein, when different peak value occurs, be with
Maximum is as the peak value of final foundation.And give, by this meter constant, full well section of logging well, by its all well logging numerical value plus being somebody's turn to do
Meter constant so that two well peak values are in same numerical value section, i.e. completes standardization correction, obtains interval transit time standardization ginseng
Number, natural gamma normalizing parameter and resistivity normalizing parameter.
The differentiation of lithology is according to interval transit time (AC), hole diameter (CAL), natural gamma (GR), resistivity (Rt) and the degree of depth
Four parameters of scale (DEPT) are carried out, such as when 100 < AC < 210 and 2 < GR < 15,130 < RT < when 350, it is determined that this lithology is
Mud stone, the rest may be inferred, and the lithology such as carbonate rock and transition lithology has corresponding parameter to differentiate interval, according to these several rocks respectively
Property respective parameter differentiate interval carry out lithology discriminant.
Above-mentioned hydrocarbon Strength co-mputation unit 3, for by interval transit time, natural gamma, resistivity substitute into lithology corresponding
Hydrocarbon strength calculation formula calculates the hydrocarbon intensity of lithology.
Specifically, it is that above-mentioned interval transit time, natural gamma, resistivity are substituted into formula
In, calculate the hydrocarbon intensity of each lithology respectively, wherein, H is the thickness of hydrocarbon source rock；ρ is the rock density in described lithology；KGasFor
Hydrocarbon rate in described lithology；TOC is the content of organic carbon in unit mass rock in described lithology.TOCiPreset for each
The content of organic carbon in unit mass rock in lithology described in TOC demarcation interval；QGas iAccording to each TOCiCount respectively
Calculate the hydrocarbon intensity obtained；I=1,2 ..., j, j are positive integer.Such as, the TOC=0.001316GR+0.00228AC-of mud stone
0.056148lgRT-0.30755, correspondingly, various lithology have each self-corresponding formula to calculate its TOC, the ginseng used
Number is identical, adjusts the coefficient before each parameter according to different lithology, numerous to list herein.
TOC is partitioned into different TOC demarcation intervals (e.g. according to 0～0.3,0.3 by the value of calculation according to described TOC
～0.6,0.6～1, > 1 four TOC demarcation intervals divide) calculate described hydrocarbon intensity respectively, then each is divided
Interval corresponding hydrocarbon intensity sums up calculating.In actual applications, follow-up can be corresponding according to above-mentioned TOC demarcation interval
Hydrocarbon intensity carries out area-method calculating, obtains the stock number of target zone, compares and does not divides the money that more can accurately calculate target zone
Source is measured.
After calculating the hydrocarbon intensity obtaining each lithology, can be strong by hydrocarbon by above-mentioned row's hydrocarbon Strength co-mputation unit 4
Degree is updated to calculate in generation and expulsion intension relational expression row's hydrocarbon intensity of each lithology.Specifically, this generation and expulsion intension relational expression is: row
Hydrocarbon intensity=hydrocarbon intensity × row's hydrocarbon rate, this row's hydrocarbon rate can obtain row's hydrocarbon rate of various lithology by experiment.
The hydrocarbon source rock life residence Strength co-mputation device of the embodiment of the present invention, it is possible to study area hydrocarbon source rock is carried out rock accurately
Property identification, by TOC and the statistical computation of hydrocarbon source rock thickness, obtains primary rock producing hydrocarbon in conjunction with parameters such as lithology and density logging, hydrocarbon rates
Hydrocarbon intensity, row's hydrocarbon intensity.Because considering many kinds of parameters factor, calculated results is more accurate, and it is residing to more conform to well logging
Actual geological environment, provide more accurate data support for follow-up well logging.
One of ordinary skill in the art will appreciate that all or part of step realizing in above-described embodiment method can be led to
Program of crossing completes to instruct relevant hardware, and this program can be stored in a computer read/write memory medium, such as
ROM/RAM, magnetic disc, CD etc..
Particular embodiments described above, has been carried out the purpose of the present invention, technical scheme and beneficial effect the most in detail
Describe in detail bright, be it should be understood that the specific embodiment that the foregoing is only the present invention, the guarantor being not intended to limit the present invention
Protect scope, all within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. done, should be included in this
Within the protection domain of invention.
Claims (10)
1. a hydrocarbon source rock life residence strength calculation method, it is characterised in that described hydrocarbon source rock life residence strength calculation method
Including:
Obtaining the log data of individual well, described log data at least includes: interval transit time, natural gamma, resistivity and the degree of depth are carved
Degree；
Differentiate that interval identifies described list according to described interval transit time, natural gamma, resistivity and deep calibration and respective parameter
Lithology in well；
Described interval transit time, natural gamma, resistivity are substituted into hydrocarbon strength calculation formula corresponding to described lithology and calculates described
The hydrocarbon intensity of lithology；
Described hydrocarbon intensity is substituted into generation and expulsion intension relational expression and calculates row's hydrocarbon intensity of described lithology.
Hydrocarbon source rock life residence strength calculation method the most according to claim 1, it is characterised in that during according to described sound wave
Difference, natural gamma, resistivity and deep calibration and respective parameter differentiate that interval identifies the lithology in described individual well, specifically wrap
Include:
Interval transit time, natural gamma and resistivity described in statistics falls the frequency in respective demarcation interval, obtains institute respectively
State the probability distribution of interval transit time, natural gamma and resistivity；
Interval transit time is determined with the probability distribution of interval transit time, natural gamma and resistivity in key hole according to described probability distribution
Normalizing parameter, natural gamma normalizing parameter and resistivity normalizing parameter；
According to described interval transit time normalizing parameter, natural gamma normalizing parameter, resistivity normalizing parameter and the described degree of depth
Scale and respective parameter differentiate that interval identifies the lithology in described individual well.
Hydrocarbon source rock life residence strength calculation method the most according to claim 2, it is characterised in that according to described probability distribution
Interval transit time normalizing parameter, natural gamma is determined with the probability distribution of interval transit time, natural gamma and resistivity in key hole
Normalizing parameter and resistivity normalizing parameter, specifically include:
Respectively by the interval transit time of described individual well, natural gamma, resistivity probability distribution in numerical value section residing for peak value, with institute
The difference stating the numerical value section residing for the peak value in the probability distribution of the interval transit time of key hole, natural gamma, resistivity is defined as
Interval transit time, natural gamma, the meter constant of resistivity；
By the interval transit time of described individual well, natural gamma, resistivity respectively plus respective meter constant, obtain described sound wave
Time difference normalizing parameter, natural gamma normalizing parameter and resistivity normalizing parameter.
Hydrocarbon source rock life residence strength calculation method the most according to claim 1, it is characterised in that described hydrocarbon intensitometer
Calculation formula is:
Wherein, H is the thickness of hydrocarbon source rock；ρ is the rock density in described lithology；KGasFor the hydrocarbon rate in described lithology；TOCiFor
The content of organic carbon in unit mass rock in lithology described in each TOC demarcation interval preset；QGas iAccording to each TOCi
It is respectively calculated the hydrocarbon intensity of acquisition；I=1,2 ..., j, j are positive integer.
Hydrocarbon source rock life residence strength calculation method the most according to claim 1, it is characterised in that described generation and expulsion intension
Relational expression is: row's hydrocarbon intensity=hydrocarbon intensity × row's hydrocarbon rate.
6. a hydrocarbon source rock life residence Strength co-mputation device, it is characterised in that described hydrocarbon source rock life residence Strength co-mputation device
Including:
Log data acquiring unit, for obtaining the log data of individual well, described log data at least includes: interval transit time, from
So gamma, resistivity and deep calibration；
Lithology Discrimination unit, for sentencing with respective parameter according to described interval transit time, natural gamma, resistivity and deep calibration
The not interval lithology identified in described individual well；
Hydrocarbon Strength co-mputation unit, for substituting into, by described interval transit time, natural gamma, resistivity, the hydrocarbon that described lithology is corresponding
Strength calculation formula calculates the hydrocarbon intensity of described lithology；
Row's hydrocarbon Strength co-mputation unit, calculates row's hydrocarbon of described lithology for described hydrocarbon intensity is substituted into generation and expulsion intension relational expression
Intensity.
Hydrocarbon source rock life residence Strength co-mputation device the most according to claim 5, it is characterised in that described Lithology Discrimination list
Unit specifically for:
Interval transit time, natural gamma and resistivity described in statistics falls the frequency in respective demarcation interval, obtains institute respectively
State the probability distribution of interval transit time, natural gamma and resistivity；
Interval transit time is determined with the probability distribution of interval transit time, natural gamma and resistivity in key hole according to described probability distribution
Normalizing parameter, natural gamma normalizing parameter and resistivity normalizing parameter；
According to described interval transit time normalizing parameter, natural gamma normalizing parameter, resistivity normalizing parameter and the described degree of depth
Scale and respective parameter differentiate that interval identifies the lithology in described individual well.
Hydrocarbon source rock life residence Strength co-mputation device the most according to claim 7, it is characterised in that according to described probability distribution
Interval transit time normalizing parameter, natural gamma is determined with the probability distribution of interval transit time, natural gamma and resistivity in key hole
Normalizing parameter and resistivity normalizing parameter, specifically include:
Respectively by the interval transit time of described individual well, natural gamma, resistivity probability distribution in numerical value section residing for peak value, with institute
The difference stating the numerical value section residing for the peak value in the probability distribution of the interval transit time of key hole, natural gamma, resistivity is defined as
Interval transit time, natural gamma, the meter constant of resistivity；
By the interval transit time of described individual well, natural gamma, resistivity respectively plus respective meter constant, obtain described sound wave
Time difference normalizing parameter, natural gamma normalizing parameter and resistivity normalizing parameter.
Hydrocarbon source rock life residence Strength co-mputation device the most according to claim 6, it is characterised in that described hydrocarbon intensitometer
Calculation formula is:
Wherein, H is the thickness of hydrocarbon source rock；ρ is the rock density in described lithology；KGasFor the hydrocarbon rate in described lithology；TOCiFor
The content of organic carbon in unit mass rock in lithology described in each TOC demarcation interval preset；QGas iAccording to each TOCi
It is respectively calculated the hydrocarbon intensity of acquisition；I=1,2 ..., j, j are positive integer.
Hydrocarbon source rock life residence Strength co-mputation device the most according to claim 6, it is characterised in that described life residence is strong
Degree relational expression is: row's hydrocarbon intensity=hydrocarbon intensity × row's hydrocarbon rate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610395936.8A CN106094053A (en) | 2016-06-07 | 2016-06-07 | Hydrocarbon source rock hydrocarbon generation and discharge strength calculation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610395936.8A CN106094053A (en) | 2016-06-07 | 2016-06-07 | Hydrocarbon source rock hydrocarbon generation and discharge strength calculation method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106094053A true CN106094053A (en) | 2016-11-09 |
Family
ID=57448410
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610395936.8A Pending CN106094053A (en) | 2016-06-07 | 2016-06-07 | Hydrocarbon source rock hydrocarbon generation and discharge strength calculation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106094053A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106886058A (en) * | 2017-04-17 | 2017-06-23 | 中国石油大学(北京) | The determination method and apparatus of compact oil reservoir Effective Reservoirs physical property lower limit |
| CN111624676A (en) * | 2019-02-27 | 2020-09-04 | 中国石油天然气股份有限公司 | Hydrocarbon-generating potential prediction method for hydrocarbon source rock |
| CN112832738A (en) * | 2021-01-05 | 2021-05-25 | 中国石油天然气股份有限公司 | Clastic rock cumulative hydrocarbon generation strength determination method and dessert layer identification and evaluation method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2747426A1 (en) * | 2009-01-09 | 2010-07-15 | Exxonmobil Upstream Research Company | Hydrocarbon detection with passive seismic data |
| CN104111483A (en) * | 2014-03-26 | 2014-10-22 | 中国石油化工股份有限公司 | Exploration method of faulted basin slope oil reservoir distribution |
| US20150233846A1 (en) * | 2014-02-14 | 2015-08-20 | Chevron U.S.A. Inc. | System and method for identifying hydrocarbon potential in a rock formation using x-ray fluorescence |
| CN105257285A (en) * | 2015-10-28 | 2016-01-20 | 中国石油大学(华东) | Evaluation method of oil gas filling capacity from source rock layer to karstic reservoir |
-
2016
- 2016-06-07 CN CN201610395936.8A patent/CN106094053A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2747426A1 (en) * | 2009-01-09 | 2010-07-15 | Exxonmobil Upstream Research Company | Hydrocarbon detection with passive seismic data |
| US20150233846A1 (en) * | 2014-02-14 | 2015-08-20 | Chevron U.S.A. Inc. | System and method for identifying hydrocarbon potential in a rock formation using x-ray fluorescence |
| CN104111483A (en) * | 2014-03-26 | 2014-10-22 | 中国石油化工股份有限公司 | Exploration method of faulted basin slope oil reservoir distribution |
| CN105257285A (en) * | 2015-10-28 | 2016-01-20 | 中国石油大学(华东) | Evaluation method of oil gas filling capacity from source rock layer to karstic reservoir |
Non-Patent Citations (2)
| Title |
|---|
| 刘海英: "呼和湖凹陷烃源岩评价", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 魏之焜: "鄂尔多斯盆地城壕地区长7段致密油成藏条件研究", 《中国优秀硕士论文全文数据库 基础科学辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106886058A (en) * | 2017-04-17 | 2017-06-23 | 中国石油大学(北京) | The determination method and apparatus of compact oil reservoir Effective Reservoirs physical property lower limit |
| CN106886058B (en) * | 2017-04-17 | 2019-04-23 | 中国石油大学(北京) | The determination method and apparatus of compact oil reservoir effective reservoir physical property lower limit value |
| CN111624676A (en) * | 2019-02-27 | 2020-09-04 | 中国石油天然气股份有限公司 | Hydrocarbon-generating potential prediction method for hydrocarbon source rock |
| CN112832738A (en) * | 2021-01-05 | 2021-05-25 | 中国石油天然气股份有限公司 | Clastic rock cumulative hydrocarbon generation strength determination method and dessert layer identification and evaluation method |
| CN112832738B (en) * | 2021-01-05 | 2023-08-22 | 中国石油天然气股份有限公司 | Determination method for strength of hydrocarbon accumulated in clastic rock and identification and evaluation method for dessert layer |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105134191B (en) | Method for evaluating reserves of tight oil well | |
| CN103993871B (en) | Method and device for processing well logging information of thin interbed stratums in standardization mode | |
| CN111665560B (en) | Hydrocarbon reservoir identification method, apparatus, computer device and readable storage medium | |
| CN103437760B (en) | Method for rapidly evaluating oil-water layer by using array induction data | |
| CN104278991A (en) | Multivariate well logging computing method for total organic carbon and hydrocarbon generating potential of hydrocarbon source rocks in salt-lake facies | |
| CN117648590B (en) | Omnibearing gamma logging data optimization processing method | |
| CN106094053A (en) | Hydrocarbon source rock hydrocarbon generation and discharge strength calculation method and device | |
| CN103678778B (en) | Method for radioactive geophysical and geochemical exploration information integration | |
| CN105240006B (en) | A kind of oil-water layer identification method suitable for volcanic rock reservoir | |
| CN107977480A (en) | A kind of shale gas reservoir aerogenesis fast appraisement method | |
| EA016499B1 (en) | Automated borehole image interpretation | |
| CN110727035A (en) | Low-permeability strong heterogeneous gas reservoir gas-water layer identification method | |
| CN113503156B (en) | Quantitative determination method, system, medium, equipment and terminal for vertical heterogeneity of reservoir | |
| CN112576244B (en) | Porosity calculation method and device | |
| CN105257285A (en) | Evaluation method of oil gas filling capacity from source rock layer to karstic reservoir | |
| CN116168224A (en) | Machine learning lithology automatic identification method based on imaging gravel content | |
| Barach et al. | Development and identification of petrophysical rock types for effective reservoir characterization: Case study of the Kristine Field, Offshore Sabah | |
| CN106199754B (en) | Oil gas drilling target integrates optimizing evaluation method | |
| CN104463686A (en) | Method for identifying shale gas reservoir while drilling by using discriminant analysis method | |
| CN112989257B (en) | Gas production amount measuring method for sea shale oil-gas reservoir | |
| CN114047547B (en) | Quantitative evaluation method for sealing capacity of gypsum-salt rock type cover layer | |
| CN112541523B (en) | Method and device for calculating mud content | |
| Deng et al. | Machine learning based stereoscopic triple sweet spot evaluation method for shale reservoirs | |
| CN112069444B (en) | Method and computer for calculating reservoir well testing permeability by using well logging data | |
| CN109598049B (en) | Method for drilling rock fracture development degree and regional rock fracture development rule |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161109 |