CN108427143A - The quantitatively characterizing method and device of rock stratum internal fissure development characteristics - Google Patents
The quantitatively characterizing method and device of rock stratum internal fissure development characteristics Download PDFInfo
- Publication number
- CN108427143A CN108427143A CN201810445569.7A CN201810445569A CN108427143A CN 108427143 A CN108427143 A CN 108427143A CN 201810445569 A CN201810445569 A CN 201810445569A CN 108427143 A CN108427143 A CN 108427143A
- Authority
- CN
- China
- Prior art keywords
- crack
- rock stratum
- slope
- accumulation curve
- internal fissure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Abstract
The invention discloses a kind of quantitatively characterizing method and device of rock stratum internal fissure development characteristics, this method includes：Lithologic log explanation is carried out based on log and electric imaging logging image, is obtained along the lithologic interpretation result on depth direction；Single layer lithologic body will be formed along lithology point consecutive identical on depth direction, obtains single lithology body thickness；Crack is released according to the electric imaging logging image characteristic solution, and counts the Number of Fractures in each single lithologic body；Fracture layers density is determined based on Number of Fractures in the single lithology body thickness and single lithologic body；Number of Fractures cumulative frequency crack frequency accumulation curve corresponding with single formation depth is generated based on the Number of Fractures in the single lithologic body；Quantitatively characterizing is carried out to rock stratum internal fissure development characteristics based on the fracture layers density and crack frequency accumulation curve.The present invention can realize the quantitatively characterizing of rock stratum internal fissure development characteristics.
Description
Technical field
The present invention relates to reservoir assessment technology field, more particularly to a kind of quantitatively characterizing side of rock stratum internal fissure development characteristics
Method and device.
Background technology
Crack is not only the important reservoir space of volcanic reservoirs, and is the main thoroughfare of oil and gas flow, can
Control oil-gas migration and aggregation.Therefore, fractue spacing properties study determines exploration success or failure, influences economic benefit.
For many years, researcher is illustrated below for crack developing numerous studies.
1, from geology angle research Causes of Cracking, FRACTURE CHARACTERISTICS, fracture stage when, lay particular emphasis on the microcosmic development characteristics in crack
Qualitative description.What it is due to crack is in vertically and horizontally inhomogeneous distribution feature, causes the uncertainty of macroscopical Study on regularity larger,
Across comparison scarce capacity.
2, the crack identification method and technology research of well logging angle and parameter calculate, and can obtain accurate fracture parameters,
Including fracture spacing, fracture porosity etc., subsequently it is fixed that development degree of micro cracks in oil can be carried out based on fracture spacing, fracture porosity
Amount evaluation, but the regularity of distribution is then difficult to describe.
3, from the FRACTURE PREDICTION technical research of earthquake angle, micro-fracture, fracture developing zone can be predicted, but be difficult to obtain
Obtain fracture parameters.
4, the crack modeling based on geomechanics etc. and numerical simulation, carry out mainly for the development mechanism of build joint
Research, is combined with stress field, can predict fractue spacing to a certain extent.But since model and parameter can not reflect reality
Border strata condition causes error larger.
As it can be seen that fracture development mechanism, feature and study on the efficiency are more deep at present, but single lithologic body internal fissure is sent out
The depth of the characterization and Study on regularity of educating degree is inadequate, still based on geologic description, lacks quantitatively characterizing method.
Invention content
The object of the present invention is to provide a kind of quantitatively characterizing methods and device of rock stratum internal fissure development characteristics, can overcome
Defect in the prior art can realize the quantitatively characterizing of rock stratum internal fissure development characteristics.
Following technical proposal can be used to realize in the above-mentioned purpose of the present invention：
A kind of quantitatively characterizing method of rock stratum internal fissure development characteristics comprising：
For specific rock type, lithologic log explanation is carried out based on log and electric imaging logging image, obtains edge
The lithologic interpretation result on depth direction；Based on the lithologic interpretation as a result, by along rock consecutive identical on depth direction
Property point form single layer lithologic body, obtain single lithology body thickness；
Crack is released according to the electric imaging logging image characteristic solution, and counts the crack in each single lithologic body
Item number；
Fracture layers density is determined based on Number of Fractures in the single lithology body thickness and single lithologic body；
It is opposite with single formation depth that Number of Fractures cumulative frequency is generated based on the Number of Fractures in the single lithologic body
The crack frequency accumulation curve answered；
Based at least one of the fracture layers density and crack frequency accumulation curve to rock stratum internal fissure development characteristics
Carry out quantitatively characterizing.
In one preferred embodiment, the single lithology body thickness is not less than 0.5 meter.
In one preferred embodiment, on along depth direction, when the single lithology body thickness is less than 0.5 meter
When, which is grouped into a upper rock stratum.
In one preferred embodiment, the calculation formula of the fracture layers density is：
In above formula：FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；LiFor single rock
Crack total number in property body, unit is item.
In one preferred embodiment, the characteristics of fracture development includes distribution position of the crack in single rock stratum
It sets, the method includes：Quantitatively characterizing, tool are carried out to rock stratum internal fissure distributing position based on the crack frequency accumulation curve
Body is：
The slope in the crack frequency accumulation curve is obtained, determines the crack in single rock stratum based on the slope
Distributing position, wherein there is at least one slope in the crack frequency accumulation curve；The fractue spacing is larger in slope
Depth segment.
In one preferred embodiment, the slope obtained in the crack frequency accumulation curve, based on described
Slope determines that distributing position of the crack in single rock stratum includes：
When the crack frequency accumulation curve is oblique more than predetermined value with same from top to bottom depth direction
When rate, then crack is evenly distributed in layer；
When the crack frequency accumulation curve has from top to bottom depth direction, there are two types of slopes, and right close to top
The slope answered is more than the slope close to bottom, then crack is mainly distributed on top, and bottom is less；
When the crack frequency accumulation curve has from top to bottom depth direction, there are three types of slopes, and right close to middle part
The slope answered is less than the slope close to bottom and top, then crack is mainly distributed at the top and bottom of rock stratum, and middle part is relatively fewer.
In one preferred embodiment, the fracture layers density and crack frequency accumulation curve are sent out with rock stratum internal fissure
The relationship educated between feature includes：
The first kind：Fracture layers density is high and the crack frequency accumulation curve has the same slope more than predetermined value,
Show that single rock stratum internal fissure development degree is high and is evenly distributed；
Second class：Fracture layers density is higher and the crack frequency accumulation curve has from top to bottom depth direction
Three kinds of slopes, and it is less than the slope close to bottom and top close to the corresponding slope in middle part, show single rock stratum internal fissure development
Degree is higher, and it is less that crack is mainly distributed on rock stratum top/bottom part, middle part；
Third class：Fracture layers density is higher and the crack frequency accumulation curve has from top to bottom depth direction
Two kinds of slopes, and be more than the slope close to bottom close to top corresponding slope, show single rock stratum internal fissure development degree compared with
Height, crack are mainly distributed at the top of rock stratum, and bottom is less；
4th class：Fracture layers density is low and the crack frequency accumulation curve has the same slope less than predetermined value,
Show that single rock stratum internal fissure is owed to educate, there is a small amount of crack at each position in rock stratum；
5th class：Fracture layers density is low and the crack frequency accumulation curve has two from top to bottom depth direction
Kind slope, and it is more than the slope close to bottom close to the corresponding slope in top, show that single rock stratum internal fissure is owed and educates, top hair
Crack is given birth to, and development degree is higher than the 4th class；
6th class：Fracture layers density is low and the crack frequency accumulation curve has three from top to bottom depth direction
Kind slope, shows that single rock stratum internal fissure is owed and educates, and top/bottom part development in rock stratum has crack, and development degree is higher than the 5th class,
But it is poorer than preceding three classes.
A kind of quantitatively characterizing device of rock stratum internal fissure development characteristics comprising：
Thickness acquisition module carries out rock for being directed to specific rock type based on log and electric imaging logging image
Property well log interpretation, obtain along the lithologic interpretation result on depth direction；As a result, it will be along depth side based on the lithologic interpretation
Consecutive identical lithology point forms single layer lithologic body upwards, obtains single lithology body thickness；
Number of Fractures acquisition module for releasing crack according to the electric imaging logging image characteristic solution, and counts each
Number of Fractures in the single lithologic body；
Fracture layers density determining module, for true based on Number of Fractures in the single lithology body thickness and single lithologic body
Determine fracture layers density；
Crack frequency accumulation curve acquisition module, for generating crack item based on the Number of Fractures in the single lithologic body
Number cumulative frequency crack frequency accumulation curve corresponding with single formation depth；
Quantitatively characterizing module, for being based at least one of the fracture layers density and crack frequency accumulation curve to rock
Layer internal fissure development characteristics carry out quantitatively characterizing.
In one preferred embodiment, the fracture layers density determining module is configured as：
In above formula：FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；LiFor single rock
Crack total number in property body, unit is item.
In one preferred embodiment, the characteristics of fracture development includes distribution position of the crack in single rock stratum
It sets, the quantitatively characterizing module is configured as：Rock stratum internal fissure distributing position is carried out based on the crack frequency accumulation curve
Quantitatively characterizing is specially：
The slope in the crack frequency accumulation curve is obtained, determines the crack in single rock stratum based on the slope
Distributing position, wherein there is at least one slope in the crack frequency accumulation curve；The fractue spacing is larger in slope
Depth segment.
The features and advantages of the invention are：The quantitatively characterizing method of rock stratum internal fissure development characteristics provided herein and
Device obtains fracture layers density by Number of Fractures in single lithology body thickness and single lithologic body, and close using fracture layers
It is more accurate to spend quantitatively characterizing development degree of the crack in single lithologic body, geologic description method than before；Further
The form of single rock stratum internal fissure frequency accumulation probability curve is obtained, and explicitly points out crack using above-mentioned density and curve and exists
Main distribution pattern about layer internal fissure in petrology has been carried out standard by the development degree and distribution situation at the rock stratum positions Nei Ge
It really portrays, computer language is used in combination to be described, the quantitative assessment for later stage reservoir provides reliable basis.
With reference to following description and accompanying drawings, the particular implementation of the application is disclosed in detail, specifies the original of the application
Reason can be in a manner of adopted.It should be understood that presently filed embodiment is not so limited in range.In appended power
In the range of the spirit and terms that profit requires, presently filed embodiment includes many changes, modifications and is equal.
The feature for describing and/or showing for a kind of embodiment can be in a manner of same or similar one or more
It is used in a other embodiment, it is combined with the feature in other embodiment, or substitute the feature in other embodiment.
It should be emphasized that term "comprises/comprising" refers to the presence of feature, one integral piece, step or component when being used herein, but simultaneously
It is not excluded for the presence or additional of one or more other features, one integral piece, step or component.
Description of the drawings
Fig. 1 is a kind of step flow of the quantitatively characterizing method of rock stratum internal fissure development characteristics in the application embodiment
Figure；
Fig. 2 a, which are that the first Number of Fractures cumulative frequency is corresponding with single formation depth in the application embodiment, to be split
Stitch frequency accumulation curve graph；
Fig. 2 b, which are that second of Number of Fractures cumulative frequency is corresponding with single formation depth in the application embodiment, to be split
Stitch frequency accumulation curve graph；
Fig. 2 c, which are that the third Number of Fractures cumulative frequency is corresponding with single formation depth in the application embodiment, to be split
Stitch frequency accumulation curve graph；
Fig. 3 a are a kind of andesite Number of Fractures cumulative frequency under depth crack frequencies corresponding with single formation depth
Rate cumulative curve chart；
Fig. 3 b are andesite Number of Fractures cumulative frequency under another depth crack frequencies corresponding with single formation depth
Rate cumulative curve chart；
Fig. 4 is a kind of module signal of the quantitatively characterizing device of rock stratum internal fissure development characteristics in the application embodiment
Figure.
Specific implementation mode
Below in conjunction with the drawings and specific embodiments, elaborate to technical scheme of the present invention, it should be understood that these
Embodiment is only illustrative of the invention and is not intended to limit the scope of the invention, after having read the present invention, this field skill
Art personnel each fall in the application range as defined in the appended claims the modification of the various equivalent forms of the present invention.
Unless otherwise defined, all of technologies and scientific terms used here by the article and belong to the technical field of the application
The normally understood meaning of technical staff is identical.The term used in the description of the present application is intended merely to description tool herein
The purpose of the embodiment of body, it is not intended that in limitation the application.
It should be noted that when element is referred to as " being set to " another element, it can be directly on another element
Or there may also be elements placed in the middle.When an element is considered as " connection " another element, it can be directly connected to
To another element or it may be simultaneously present centering elements.Term as used herein " vertical ", " horizontal ", "upper",
"lower", "left", "right" and similar statement for illustrative purposes only, are not offered as being unique embodiment.
The present invention provides a kind of quantitatively characterizing method and device of rock stratum internal fissure development characteristics, can overcome the prior art
The defects of, realize the quantitatively characterizing of rock stratum internal fissure development characteristics.
Referring to Fig. 1, a kind of quantitatively characterizing method of rock stratum internal fissure development characteristics is provided in the application embodiment,
It may include steps of.
Step S10：For specific rock type, lithologic log solution is carried out based on log and electric imaging logging image
It releases, obtains along the lithologic interpretation result on depth direction；As a result, it will be along continuous on depth direction based on the lithologic interpretation
Identical lithology point forms single layer lithologic body, obtains single lithology body thickness；
Step S12：Crack is released according to electric imaging logging image characteristic solution, and is counted in each single lithologic body
Number of Fractures；
Step S14：Fracture layers density is determined based on Number of Fractures in the single lithology body thickness and single lithologic body；
Step S16：Number of Fractures cumulative frequency and single rock stratum are generated based on the Number of Fractures in the single lithologic body
The corresponding crack frequency accumulation curve of depth；
Step S18：Based at least one of the fracture layers density and crack frequency accumulation curve to rock stratum internal fissure
Development characteristics carry out quantitatively characterizing.
It is possible, firstly, to for research work area target zone geological condition, the rock type for needing to carry out lithologic interpretation is determined.
Specifically, being merged according to Research Requirements by lithology statistical result, simplifying similar lithology, determine rock type.
That is the quantitatively characterizing method of rock stratum internal fissure development characteristics provided herein be rock type explicitly before
It puts, is carried out for specific rock type.Specifically, the rock stratum internal fissure development characteristics may include rock stratum implosion
The distribution characteristics of seam and development degree etc..
After rock type determination, lithologic log explanation is carried out according to log and electric imaging logging image, obtains edge
The lithologic interpretation result on depth direction.Specifically, can be in conjunction with Oil Field coring analysis and geologic description, based on conventional
Well logging and electric imaging logging image carry out depth point lithologic interpretation one by one using intersection drawing method and characteristics of image.
Longitudinal upper lithology point continuous, identical (on depth direction) is formed into single layer lithologic body, it is thick to obtain single lithologic body
Degree.Wherein, single lithology body thickness is not less than 0.5m.When the thickness of single lithologic body is less than 0.5m, then upward merging, by the thickness
Degree is grouped into last layer.
Specifically upward the step of merging, is：In lithologic log interpretation process, if there are the rock types of some continuous depth
Type thickness is less than 0.5m, then the thickness is grouped into upper layer, and rock type is identical as upper layer.
Then crack can be released according to electric imaging logging image characteristic solution, and counted in each single lithologic body
Number of Fractures.Wherein, the method in electric logging image interpretation crack is highly developed, the application not reinflated elaboration herein.One
As in the case of, the electric imaging logging image in oblique crack is characterized as black sine curve.
It is then possible to determine fracture layers density based on Number of Fractures in the single lithology body thickness and single lithologic body.
Wherein, the calculation formula of the fracture layers density is：
In above formula：FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；LiFor single rock
Crack total number in property body, unit is item.
Followed by Number of Fractures cumulative frequency and single rock can be generated based on the Number of Fractures in the single lithologic body
The corresponding crack frequency accumulation curve of layer depth.
In the present embodiment, the Number of Fractures in the single lithologic body can be generated single rock along depth direction
Property body fractue spacing histogram and Number of Fractures generate Number of Fractures cumulative frequency crack corresponding with single formation depth
Frequency accumulation curve.Using the crack, frequency accumulation curve describes crack inhomogeneous distribution.Especially it is available with crack frequency
Slope in rate accumulation curve determines distributing position feature of the crack in single rock stratum.
In one embodiment, rock stratum internal fissure distributing position is quantified based on the crack frequency accumulation curve
It characterizes, is specially：The slope in the crack frequency accumulation curve is obtained, determines the crack single based on the slope
Distributing position in rock stratum, wherein there is at least one slope in the crack frequency accumulation curve；The fractue spacing is oblique
The larger depth segment of rate.
In general, having at least one slope in the crack frequency accumulation curve.When the slope is more than predetermined value, table
Bright rock is more in the correspondence depth fractue spacing.Wherein the predetermined value can be according to actual rock type and other parameters
Comprehensive determination is carried out, the application does not make specific limit herein.
In a specific embodiment, the slope obtained in the crack frequency accumulation curve, based on described
Slope determines that distributing position of the crack in single rock stratum may include：
When the crack frequency accumulation curve is oblique more than predetermined value with same from top to bottom depth direction
When rate, then crack is evenly distributed in layer；
When the crack frequency accumulation curve has from top to bottom depth direction, there are two types of slopes, and right close to top
The slope answered is more than the slope close to bottom, then crack is mainly distributed on top, and bottom is less；
When the crack frequency accumulation curve has from top to bottom depth direction, there are three types of slopes, and right close to middle part
The slope answered is less than the slope close to bottom and top, then crack is mainly distributed at the top and bottom of rock stratum, and middle part is relatively fewer.
Specifically, when the crack frequency accumulation curve is one-part form, then crack is evenly distributed in layer；It is split when described
When seam frequency accumulation curve is two-part, then crack is mainly distributed on top, and bottom is less；When the crack frequency accumulation curve
For three-stage when, then crack is mainly distributed at the top and bottom of rock stratum, and middle part is relatively fewer.
Above-mentioned three kinds of situations are described in detail with reference to specific crack frequency accumulation curve graph.
Fig. 2 a are please referred to, abscissa is single formation depth in figure, and unit is m (rice)；Extremely from the top of rock stratum
The depth direction of bottom gradually extends；Ordinate in figure indicates Number of Fractures cumulative frequency, unit %.
(1) if accumulation curve is one-part form, i.e., when entire accumulation curve is equivalent to a kind of slope, then it is assumed that crack exists
It is evenly distributed in layer.
Fig. 2 b are please referred to, abscissa is single formation depth in figure, and unit is m (rice)；Extremely from the top of rock stratum
The depth direction of bottom gradually extends；Ordinate in figure indicates Number of Fractures cumulative frequency, unit %, can be according to splitting
Distribution histogram is stitched to obtain.
(2) if accumulation curve be two-part, i.e., entire accumulation curve be equivalent to tool there are two types of slope when, wherein slope compared with
Close the top of (relatively steep) greatly, smaller (shallower) the close bottom of slope, then it represents that at the top of crack is mainly distributed on, bottom compared with
It is few.
If accumulation curve is two-part, each section of the slope of curve is different, and is recognized in conjunction with the generality that fractue spacing is studied
Know, the slope of first segment is more than second segment certainly, shows that crack major developmental is distributed in first segment, i.e. crack is mainly distributed on top
Portion, bottom are less.
As for the degree of strength (anisotropism distributed degrees) of this segmentation distribution characteristics, quantitative table can not be carried out at present
Sign is based on this point consideration, and the present invention is characterized point on the basis of the frequency statistics of crack by this accumulation curve form
Cloth feature further characterizes the power of anisotropism by the slope of curve.Slope difference is big, and anisotropism is strong, and slope difference is small,
It is distributed relatively uniform, is finally unistage type.Certainly, according to the size of the slope, it may be determined that fractue spacing position, slope are larger
, the fractue spacing for corresponding to depth segment is more, and slope is smaller, and the fractue spacing for corresponding to depth segment is less.
Fig. 2 c are please referred to, abscissa is single formation depth in figure, and unit is m (rice)；Extremely from the top of rock stratum
The depth direction of bottom gradually extends；Ordinate in figure indicates Number of Fractures cumulative frequency, unit %, can be according to splitting
Distribution histogram is stitched to obtain.
(3) if accumulation curve be three-stage, i.e., entire accumulation curve be equivalent to tool there are three types of slope when, wherein slope compared with
Close the top and bottom of (relatively steep) greatly, slope is smaller (shallower) to be located at middle part, then it represents that at the top of crack is mainly distributed on
And bottom, middle part are less.
If tracing pattern be three-stage, from the slope of accumulation curve can find out crack top bottom distribution it is more,
Middle part is less, also consistent with the generality of existing research understanding, while having confirmed this method and can describe crack in rock stratum
In homogenieity distribution.
The intension of ordinate is the distribution characteristics in order to express crack in single rock stratum in above-mentioned Fig. 2 a to Fig. 2 c.This is vertical
The acquisition modes of coordinate are：According to crack interpretation results, it is deep that interval constant depth (generally 2m, 5m etc.) counts each interval
Number of Fractures in degree, and the Number of Fractures for calculating each gap depth accounts for hundred of total Number of Fractures in entire single lithologic body
It is 100% to divide ratio, percentage summation.Ordinate in figure is the percentage accumulation curve of each gap depth, i.e. the latter interval
The cumulative frequency value (frequency values, that is, percentage) of depth is the summation of each gap depth cumulative frequency in front.
At least one of the fracture layers density and crack frequency accumulation curve can be finally based on to rock stratum internal fissure
Development characteristics carry out quantitatively characterizing.
Specifically, can be in summary in step for crack frequency accumulation curve and fracture layers density, to single rock stratum
Interior characteristics of fracture development carries out quantitatively characterizing.
Petrological study has shown at present：There are a variety of situations for distribution of the crack in single rock stratum, are split using current
It is strong and weak that seam parameter cannot evaluate crack distribution mode and anisotropism in rock stratum.For this purpose, proposing that the method for the present invention solves this
Problem is determined by tracing pattern and fracture layers density size synthesis, unified not to the utmost as the criteria for classifying of development degree power,
But can be with digital representation using fracture layers density, strong and weak comparison of being more convenient for.
Wherein, the relationship packet between the fracture layers density and crack frequency accumulation curve and rock stratum internal fissure development characteristics
It includes：
The first kind：Fracture layers density is high and the crack frequency accumulation curve has the same slope more than predetermined value,
Show that single rock stratum internal fissure development degree is high and is evenly distributed；
Second class：Fracture layers density is higher and the crack frequency accumulation curve has from top to bottom depth direction
Three kinds of slopes, and it is less than the slope close to bottom and top close to the corresponding slope in middle part, show single rock stratum internal fissure development
Degree is higher, and it is less that crack is mainly distributed on rock stratum top/bottom part, middle part；
Third class：Fracture layers density is higher and the crack frequency accumulation curve has from top to bottom depth direction
Two kinds of slopes, and be more than the slope close to bottom close to top corresponding slope, show single rock stratum internal fissure development degree compared with
Height, crack are mainly distributed at the top of rock stratum, and bottom is less；
4th class：Fracture layers density is low and the crack frequency accumulation curve has the same slope less than predetermined value,
Show that single rock stratum internal fissure is owed to educate, there is a small amount of crack at each position in rock stratum；
5th class：Fracture layers density is low and the crack frequency accumulation curve has two from top to bottom depth direction
Kind slope, and it is more than the slope close to bottom close to the corresponding slope in top, show that single rock stratum internal fissure is owed and educates, top hair
Crack is given birth to, and development degree is higher than the 4th class；
6th class：Fracture layers density is low and the crack frequency accumulation curve has three from top to bottom depth direction
Kind slope, shows that single rock stratum internal fissure is owed and educates, and top/bottom part development in rock stratum has crack, and development degree is higher than the 5th class,
But it is poorer than preceding three classes.
Certainly, in this embodiment, the classification that the applicant carries out just for existing common several situations, needle
To different rock type, different field conditions, those skilled in the art can also be according to fracture layers density and crack frequency
The slope characteristics of accumulation curve carry out categorical rating, and the application does not make unique limit herein.
Incorporated by reference to Fig. 3 a and Fig. 3 b, under a specific application scenarios, such as the andesite of predetermined depth,
It has different slopes in different depth, for example, between 3723 meters to 3737 meters, with same slope,
There are three types of slopes for tool between 3615 meters to 3650 meters.Certainly it is stretched out at other, which can also have other slope characteristics,
The application not reinflated narration herein.
It subsequently can be in conjunction with the calculation formula of the fracture layers density：
In above formula：FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；LiFor single rock
Crack total number in property body, unit is item, determines the density of fracture layers；
It is finally referred between the fracture layers density and crack frequency accumulation curve and rock stratum internal fissure development degree
Relationship determine each rock stratum internal fissure development characteristics of the andesite.
Relationship between above-mentioned fracture layers density and crack frequency accumulation curve and rock stratum internal fissure development degree includes：The
It is a kind of：Fracture layers density is high and the crack frequency accumulation curve is one-part form, shows that single rock stratum internal fissure development degree is high
And it is evenly distributed；
Second class：Fracture layers density is higher and the crack frequency accumulation curve is three-stage feature, shows single rock stratum
Internal fissure development degree is higher, and it is less that crack is mainly distributed on rock stratum top/bottom part, middle part；
Third class：Fracture layers density is higher and the crack frequency accumulation curve is two-part feature, shows single rock stratum
Internal fissure development degree is higher, and crack is mainly distributed at the top of rock stratum, and bottom is less；
4th class：Fracture layers density is low and the crack frequency accumulation curve is unistage type feature, shows in single rock stratum
Crack, which is owed, educates, and there is a small amount of crack at each position in rock stratum；
5th class：Fracture layers density is low and the crack frequency accumulation curve is two-part feature, shows in single rock stratum
Crack, which is owed, educates, and a small amount of crack is developed at top, and development degree is higher than the 4th class；
6th class：Fracture layers density is low and the crack frequency accumulation curve is three-stage feature, shows in single rock stratum
Crack, which is owed, educates, and rock stratum top/bottom part develops a small amount of crack, and development degree is higher than the 5th class, but poorer than preceding three classes.
In the present embodiment, the high or low of the fracture layers density is an opposite concept.For different oil fields, no
There is various criterion in same layer position.But it is directed to particular studies area, a unified standard can be set, in the system according to actual conditions
One standard judges the opposite height of fracture layers density.Specifically, the application not reinflated description herein.
The quantitatively characterizing method of rock stratum internal fissure development characteristics provided herein, passes through single lithology body thickness and list
Number of Fractures obtains fracture layers density in one lithologic body, and using fracture layers density quantitatively characterizing crack in single lithologic body
In development degree, geologic description method than before is more accurate；Further obtain single rock stratum internal fissure frequency accumulation
The form of probability curve, and explicitly point out development degree and distribution of the crack at the rock stratum positions Nei Ge using above-mentioned density and curve
Main distribution pattern about layer internal fissure in petrology accurately portray, computer language is used in combination to be retouched by situation
It states, the quantitative assessment for later stage reservoir provides reliable basis.
The present invention is suitable for the fracture development distribution characteristics research in the fracture-type reservoirs such as volcanic rock, carbonate rock.In rock
On the basis of stone type and single lithology body thickness determine, can preferable Prediction of fracture degree and distribution characteristics, make up existing
The deficiency for having technology provides reliable basis for fracture-type reservoir efficiency evaluation and reservoir reconstruction, has a extensive future.
Please refer to Fig. 4, for the above embodiment provide rock stratum internal fissure development characteristics quantitatively characterizing method,
A kind of quantitatively characterizing device of rock stratum internal fissure development characteristics is correspondingly provided in the application embodiment, which can wrap
It includes：
Thickness acquisition module 10 is carried out for being directed to specific rock type based on log and electric imaging logging image
Lithologic log is explained, is obtained along the lithologic interpretation result on depth direction；As a result, it will be along depth based on the lithologic interpretation
Consecutive identical lithology point forms single layer lithologic body on direction, obtains single lithology body thickness；
Number of Fractures acquisition module 12 for releasing crack according to the electric imaging logging image characteristic solution, and counts every
Number of Fractures in a single lithologic body；
Fracture layers density determining module 14, for based on Number of Fractures in the single lithology body thickness and single lithologic body
Determine fracture layers density；
Crack frequency accumulation curve acquisition module 16, for generating crack based on the Number of Fractures in the single lithologic body
Item number cumulative frequency crack frequency accumulation curve corresponding with single formation depth；
Quantitatively characterizing module 18, for being based at least one of the fracture layers density and crack frequency accumulation curve pair
Rock stratum internal fissure development characteristics carry out quantitatively characterizing.
Wherein, the fracture layers density determining module 14 is configured as：
In above formula：FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；LiFor single rock
Crack total number in property body, unit is item.
In the present embodiment, modules are realized in the quantitatively characterizing device of the rock stratum internal fissure development characteristics
Concrete function corresponds to each step in the quantitatively characterizing method of rock stratum internal fissure development characteristics.Specifically, in the rock stratum
The quantitatively characterizing device of characteristics of fracture development may include memory and processor, and computer program is stored in memory, described
Computer program realizes following steps when being executed by the processor：
Step S10：For specific rock type, lithologic log solution is carried out based on log and electric imaging logging image
It releases, obtains along the lithologic interpretation result on depth direction；As a result, it will be along continuous on depth direction based on the lithologic interpretation
Identical lithology point forms single layer lithologic body, obtains single lithology body thickness；
Step S12：Crack is released according to electric imaging logging image characteristic solution, and is counted in each single lithologic body
Number of Fractures；
Step S14：Fracture layers density is determined based on Number of Fractures in the single lithology body thickness and single lithologic body；
Step S16：Number of Fractures cumulative frequency and single rock stratum are generated based on the Number of Fractures in the single lithologic body
The corresponding crack frequency accumulation curve of depth；
Step S18：Rock stratum internal fissure development characteristics are carried out based on the fracture layers density and crack frequency accumulation curve
Quantitatively characterizing.
In the present embodiment, the memory may include the physical unit for storing information, typically by information
It is stored again with the media using the methods of electricity, magnetic or optics after digitlization.Memory described in present embodiment again may be used
To include：The device of information, such as RAM, ROM are stored in the way of electric energy；The device of information is stored in the way of magnetic energy, it is such as hard
Disk, floppy disk, tape, core memory, magnetic bubble memory, USB flash disk；The device that information is stored using optical mode, such as CD or DVD.
Certainly, also memory of other modes, such as quantum memory, graphene memory etc..
In the present embodiment, the processor can be implemented in any suitable manner.For example, the processor can be with
Take such as microprocessor or processor and storage can by computer readable program code that (micro-) processor executes (such as
Software or firmware) computer-readable medium, logic gate, switch, application-specific integrated circuit (Application Specific
Integrated Circuit, ASIC), the form etc. of programmable logic controller (PLC) and embedded microcontroller.
In one embodiment, there is at least one slope in the crack frequency accumulation curve；When crack frequency
When rate accumulation curve is one-part form, then crack is evenly distributed in layer；When the crack frequency accumulation curve is two-part, then
Crack is mainly distributed on top, and bottom is less；When the crack frequency accumulation curve is three-stage, then crack is mainly distributed on
At the top and bottom of rock stratum, middle part is relatively fewer.
The quantitatively characterizing method and device of internal fissure development characteristics in rock stratum provided by the present invention can be adapted for volcanic rock,
Fracture development distribution characteristics research in the fracture-type reservoirs such as carbonate rock.Base is determined in rock type and single lithology body thickness
On plinth, can preferable Prediction of fracture degree and distribution characteristics, make up the deficiencies in the prior art, be that fracture-type reservoir is effective
Property evaluation and reservoir reconstruction provide reliable basis, have a extensive future.
Herein cited any digital value all include between lower limiting value to upper limit value with the lower value of an incremented and
The all values of upper value, there are the intervals of at least two units between any lower value and any much higher value.For example, such as
Fruit elaborates that the quantity of component or the value of process variable (such as temperature, pressure, time etc.) are from 1 to 90, preferably from 20
To 80, more preferably from 30 to 70, then purpose is arrived in order to illustrate also clearly listing such as 15 to 85,22 in the specification
68,43 to 51,30 to 32 is equivalent.For the value less than 1, suitably think that a unit is 0.0001,0.001,0.01,0.1.
These are only intended to the example clearly expressed, it is believed that the numerical value enumerated between minimum and peak is possible to
Combination is all expressly set forth in the specification in a similar manner.
Unless otherwise indicated, all ranges all include all numbers between endpoint and endpoint.It is used together with range
" about " or " approximation " be suitable for two endpoints of the range.Thus, " about 20 to 30 " are intended to covering " about 20 to about
30 ", including at least the endpoint indicated.
All articles and reference disclosed herein, including patent application and publication, for various purposes by helping
Draw and is hereby incorporated by.Describe combination term " substantially by ... constitutes " should including identified element, ingredient, component or step with
And other elements, ingredient, component or step essentially without the basic novel feature for influencing the combination.Use term "comprising"
Or " comprising " describes the combination of element here, ingredient, component or step it is also contemplated that substantially by these elements, ingredient, portion
The embodiment that part or step are constituted.Here by using term " can with ", it is intended to illustrate described that " can with " includes
What attribute is all optional.
Multiple element, ingredient, component or step can be provided by single integrated component, ingredient, component or step.Optionally
Ground, single integrated component, ingredient, component or step can be divided into multiple element, ingredient, component or the step of separation.It is used for
The open "a" or "an" for describing element, ingredient, component or step is not said to exclude other elements, ingredient, component
Or step.
Above-mentioned each embodiment in this specification is described in a progressive manner, identical between each embodiment
Similar portion is cross-referenced, and what each embodiment stressed is and other embodiment difference.
The foregoing is merely several embodiments of the invention, although disclosed herein embodiment it is as above, institute
Content is stated only to facilitate the embodiment for understanding the present invention and using, is not intended to limit the present invention.Any institute of the present invention
Belong to those skilled in the art, do not depart from disclosed herein spirit and scope under the premise of, can be in embodiment
Formal and details on make any modification and variation, but the scope of patent protection of the present invention still must be with appended claims
Subject to the range that book is defined.
Claims (10)
1. a kind of quantitatively characterizing method of rock stratum internal fissure development characteristics, which is characterized in that including：
For specific rock type, lithologic log explanation is carried out based on log and electric imaging logging image, is obtained along depth
Spend the lithologic interpretation result on direction；Based on the lithologic interpretation as a result, by along lithology point consecutive identical on depth direction
Single layer lithologic body is formed, single lithology body thickness is obtained；
Crack is released according to the electric imaging logging image characteristic solution, and counts the crack item in each single lithologic body
Number；
Fracture layers density is determined based on Number of Fractures in the single lithology body thickness and single lithologic body；
It is corresponding with single formation depth that Number of Fractures cumulative frequency is generated based on the Number of Fractures in the single lithologic body
Crack frequency accumulation curve；
Rock stratum internal fissure development characteristics are carried out based at least one of the fracture layers density and crack frequency accumulation curve
Quantitatively characterizing.
2. the quantitatively characterizing method of internal fissure development characteristics in rock stratum as described in claim 1, which is characterized in that the single rock
Property body thickness be not less than 0.5 meter.
3. the quantitatively characterizing method of internal fissure development characteristics in rock stratum as claimed in claim 2, which is characterized in that along depth
On direction, when the single lithology body thickness is less than 0.5 meter, which is grouped into a upper rock stratum.
4. the quantitatively characterizing method of internal fissure development characteristics in rock stratum as described in claim 1, which is characterized in that the fracture layers
The calculation formula of density is：
In above formula：FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；LiFor single lithologic body
In crack total number, unit is item.
5. the quantitatively characterizing method of internal fissure development characteristics in rock stratum as described in claim 1, which is characterized in that the crack hair
It includes distributing position of the crack in single rock stratum to educate feature, the method includes：Based on the crack frequency accumulation curve pair
Rock stratum internal fissure distributing position carries out quantitatively characterizing, is specially：
The slope in the crack frequency accumulation curve is obtained, point of the crack in single rock stratum is determined based on the slope
Cloth position, wherein there is at least one slope in the crack frequency accumulation curve；The fractue spacing depth larger in slope
Spend section.
6. the quantitatively characterizing method of internal fissure development characteristics in rock stratum as claimed in claim 5, which is characterized in that the acquisition institute
The slope in the frequency accumulation curve of crack is stated, distributing position packet of the crack in single rock stratum is determined based on the slope
It includes：
When the crack frequency accumulation curve has same slope more than predetermined value from top to bottom depth direction,
Then crack is evenly distributed in layer；
When the crack frequency accumulation curve has from top to bottom depth direction, there are two types of slopes, and corresponding close to top
Slope is more than the slope close to bottom, then crack is mainly distributed on top, and bottom is less；
When the crack frequency accumulation curve has from top to bottom depth direction, there are three types of slopes, and corresponding close to middle part
Slope is less than the slope close to bottom and top, then crack is mainly distributed at the top and bottom of rock stratum, and middle part is relatively fewer.
7. the quantitatively characterizing method of internal fissure development characteristics in rock stratum as claimed in claim 6, which is characterized in that the fracture layers
Relationship between density and crack frequency accumulation curve and rock stratum internal fissure development characteristics includes：
The first kind：Fracture layers density is high and the crack frequency accumulation curve has the same slope more than predetermined value, shows
Single rock stratum internal fissure development degree is high and is evenly distributed；
Second class：There are three types of fracture layers density is higher and the crack frequency accumulation curve has from top to bottom depth direction
Slope, and it is less than the slope close to bottom and top close to the corresponding slope in middle part, show single rock stratum internal fissure development degree
Higher, it is less that crack is mainly distributed on rock stratum top/bottom part, middle part；
Third class：There are two types of fracture layers density is higher and the crack frequency accumulation curve has from top to bottom depth direction
Slope, and it is more than the slope close to bottom close to the corresponding slope in top, show that single rock stratum internal fissure development degree is higher, splits
Seam is mainly distributed at the top of rock stratum, and bottom is less；
4th class：Fracture layers density is low and the crack frequency accumulation curve has the same slope less than predetermined value, shows
Single rock stratum internal fissure, which is owed, educates, and there is a small amount of crack at each position in rock stratum；
5th class：Fracture layers density is low and the crack frequency accumulation curve has that there are two types of tiltedly from top to bottom depth direction
Rate, and it is more than the slope close to bottom close to the corresponding slope in top, show that single rock stratum internal fissure is owed and educate, top development has
Crack, and development degree is higher than the 4th class；
6th class：Fracture layers density is low and the crack frequency accumulation curve has that there are three types of tiltedly from top to bottom depth direction
Rate shows that single rock stratum internal fissure is owed and educates that top/bottom part development in rock stratum has crack, and development degree is higher than the 5th class, but compares
Preceding three classes are poor.
8. a kind of quantitatively characterizing device of rock stratum internal fissure development characteristics, which is characterized in that it includes：
Thickness acquisition module carries out lithology survey for being directed to specific rock type based on log and electric imaging logging image
Well is explained, is obtained along the lithologic interpretation result on depth direction；As a result, it will be along on depth direction based on the lithologic interpretation
Consecutive identical lithology point forms single layer lithologic body, obtains single lithology body thickness；
Number of Fractures acquisition module for releasing crack according to the electric imaging logging image characteristic solution, and counts each described
Number of Fractures in single lithologic body；
Fracture layers density determining module, for being split based on Number of Fractures determination in the single lithology body thickness and single lithologic body
Stitch layer density；
Crack frequency accumulation curve acquisition module, it is tired for generating Number of Fractures based on the Number of Fractures in the single lithologic body
Product frequency crack frequency accumulation curve corresponding with single formation depth；
Quantitatively characterizing module, for being based at least one of the fracture layers density and crack frequency accumulation curve in rock stratum
Characteristics of fracture development carries out quantitatively characterizing.
9. the quantitatively characterizing device of internal fissure development characteristics in rock stratum as claimed in claim 8, which is characterized in that the fracture layers
Density determining module is configured as：
In above formula：FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；LiFor single lithologic body
In crack total number, unit is item.
10. the quantitatively characterizing device of internal fissure development characteristics in rock stratum as claimed in claim 8, which is characterized in that the crack
Development characteristics include distributing position of the crack in single rock stratum, and the quantitatively characterizing module is configured as：Based on the crack
Frequency accumulation curve carries out quantitatively characterizing to rock stratum internal fissure distributing position, is specially：
The slope in the crack frequency accumulation curve is obtained, point of the crack in single rock stratum is determined based on the slope
Cloth position, wherein there is at least one slope in the crack frequency accumulation curve；The fractue spacing depth larger in slope
Spend section.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810445569.7A CN108427143B (en) | 2018-05-11 | 2018-05-11 | The quantitatively characterizing method and device of rock stratum internal fissure development characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810445569.7A CN108427143B (en) | 2018-05-11 | 2018-05-11 | The quantitatively characterizing method and device of rock stratum internal fissure development characteristics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108427143A true CN108427143A (en) | 2018-08-21 |
| CN108427143B CN108427143B (en) | 2019-10-11 |
Family
ID=63162638
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810445569.7A Active CN108427143B (en) | 2018-05-11 | 2018-05-11 | The quantitatively characterizing method and device of rock stratum internal fissure development characteristics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108427143B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170115414A1 (en) * | 2015-10-27 | 2017-04-27 | Schlumberger Technology Corporation | Determining shear slowness based on a higher order formation flexural acoustic mode |
| CN107917865A (en) * | 2016-10-11 | 2018-04-17 | 中国石油化工股份有限公司 | A kind of tight sandstone reservoir multi-parameter Permeability Prediction method |
| CN107956472A (en) * | 2017-11-10 | 2018-04-24 | 中国石油天然气股份有限公司 | Fissured carbonate reservoir completion mode determines method, apparatus and medium |
| CN108008467A (en) * | 2017-11-29 | 2018-05-08 | 中国科学院地质与地球物理研究所兰州油气资源研究中心 | Fractue spacing quantitatively characterizing method and its system |
| CN108008464A (en) * | 2017-11-29 | 2018-05-08 | 中国科学院地质与地球物理研究所兰州油气资源研究中心 | Crack anisotropism quantitatively characterizing method and its system |
-
2018
- 2018-05-11 CN CN201810445569.7A patent/CN108427143B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170115414A1 (en) * | 2015-10-27 | 2017-04-27 | Schlumberger Technology Corporation | Determining shear slowness based on a higher order formation flexural acoustic mode |
| CN107917865A (en) * | 2016-10-11 | 2018-04-17 | 中国石油化工股份有限公司 | A kind of tight sandstone reservoir multi-parameter Permeability Prediction method |
| CN107956472A (en) * | 2017-11-10 | 2018-04-24 | 中国石油天然气股份有限公司 | Fissured carbonate reservoir completion mode determines method, apparatus and medium |
| CN108008467A (en) * | 2017-11-29 | 2018-05-08 | 中国科学院地质与地球物理研究所兰州油气资源研究中心 | Fractue spacing quantitatively characterizing method and its system |
| CN108008464A (en) * | 2017-11-29 | 2018-05-08 | 中国科学院地质与地球物理研究所兰州油气资源研究中心 | Crack anisotropism quantitatively characterizing method and its system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108427143B (en) | 2019-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8370122B2 (en) | Method of predicting connectivity between parts of a potential hydrocarbon reservoir and analyzing 3D data in a subsurface region | |
| Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
| KR101148835B1 (en) | Prediction system and method for subsurface lithology in oil sands reservoir using statistical analysis of well logging data | |
| CN111665560B (en) | Hydrocarbon reservoir identification method, apparatus, computer device and readable storage medium | |
| AU2014287371A1 (en) | Determining geomechanics completion quality | |
| EP2113796A1 (en) | Method and Apparatus for Analyzing Three-Dimensional Data | |
| CN102749647B (en) | Quantitative method for space of holes in reservoir | |
| Yang et al. | Origin and architecture of fractured-cavernous carbonate reservoirs and their influences on seismic amplitudes | |
| Zhang et al. | 3D porosity modeling of a carbonate reservoir using continuous multiple-point statistics simulation | |
| Euzen et al. | Well log cluster analysis: an innovative tool for unconventional exploration | |
| CN103376468A (en) | Reservoir parameter quantitative characterization method based on neural network function approximation algorithm | |
| CN109298464A (en) | Tight sandstone reservoir Diagenetic Facies Logging Identification Method and device | |
| CN108427143B (en) | The quantitatively characterizing method and device of rock stratum internal fissure development characteristics | |
| CN110566194B (en) | Comprehensive quantitative evaluation method and device for multilayer oil-containing system reservoir | |
| US9880321B2 (en) | Defining non-linear petrofacies for a reservoir simulation model | |
| CN116796231A (en) | Method and device for automatically dividing lithology and comparing lithology based on logging curve | |
| CN112282751B (en) | Geological engineering three-dimensional coupling compact oil gas horizontal well exploitation detection method | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| Palmer | Using AI and Machine Learning to Indicate Shale Anisotropy and Assist in Completions Design | |
| CN114153002A (en) | Three-dimensional geological modeling method and device for natural fracture of reservoir, electronic equipment and medium | |
| CN106285661A (en) | A kind of based on the Low Resistivity Reservoir Identification Methods and the device that judge index | |
| Zambrano et al. | Influence of geological features on well test behavior | |
| Sahin et al. | Stochastic modeling of porosity distribution in a multi-zonal carbonate reservoir | |
| Ateeq et al. | An integrated study to characterize and model natural fracture networks of gas condensate carbonate reservoirs, onshore Abu Dhabi, UAE | |
| Mehranfar et al. | A comprehensive workflow for reservoir characterization and history matching applicable to cretaceous naturally fractured reservoirs in Mexico |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |