CN108427143B - 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
- CN108427143B CN108427143B CN201810445569.7A CN201810445569A CN108427143B CN 108427143 B CN108427143 B CN 108427143B CN 201810445569 A CN201810445569 A CN 201810445569A CN 108427143 B CN108427143 B CN 108427143B
- 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.)
- Active
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 the quantitatively characterizing methods and device of a kind of rock stratum internal fissure development characteristics, this method comprises: carrying out lithologic log explanation based on log and electric imaging logging image, obtain 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 be realized the quantitatively characterizing of rock stratum internal fissure development characteristics.
Description
Technical field
The present invention relates to reservoir assessment technology field, in particular to a kind of quantitatively characterizing side of rock stratum internal fissure development characteristics
Method and device.
Background technique
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, when from geology angle research Causes of Cracking, FRACTURE CHARACTERISTICS, fracture stage, the microcosmic development characteristics in crack are laid particular emphasis on
Qualitative description.What it is due to crack is vertically and horizontally in inhomogeneous distribution feature, causes the uncertainty of macroscopical Study on regularity larger,
Across comparison scarce capacity.
2, the crack identification method and technology research for angle of logging well and parameter calculate, and can obtain accurate fracture parameters,
It is subsequent it is fixed to carry out development degree of micro cracks in oil based on fracture spacing, fracture porosity including fracture spacing, fracture porosity etc.
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 is difficult to obtain
Obtain fracture parameters.
4, the crack modeling and numerical simulation based on geomechanics etc., carries out mainly for the development mechanism of build joint
Research, combines 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.
Summary of the invention
The object of the present invention is to provide the quantitatively characterizing methods and device of a kind of rock stratum internal fissure development characteristics, can overcome
Defect in the prior art can be realized the quantitatively characterizing of rock stratum internal fissure development characteristics.
Following technical proposal can be used to realize in above-mentioned purpose of the 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 are as follows:
In above formula: FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；∑LiIt is single
Crack total number in lithologic body, unit are item.
In one preferred embodiment, the characteristics of fracture development includes distribution position of the crack in single rock stratum
It sets, which comprises quantitatively characterizing, tool are carried out to rock stratum internal fissure distributing position based on the crack frequency accumulation curve
Body are as follows:
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 greater 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 greater 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 and rock stratum internal fissure are sent out
The relationship educated between feature includes:
The first kind: fracture layers density is high and the crack frequency accumulation curve has the same slope greater 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 close to the corresponding slope in top be greater than close to bottom 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 greater 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 obtains module, for being directed to specific rock type, carries out rock 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 obtains 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；∑LiIt is single
Crack total number in lithologic body, unit are 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 being carried out based on the crack frequency accumulation curve
Quantitatively characterizing, specifically:
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 in petrology about layer internal fissure has been carried out standard by the development degree and distribution situation at the rock stratum position Nei Ge
It really portrays, and is described with computer language, the quantitative assessment for later period reservoir provides reliable basis.
Referring to following description and accompanying drawings, specific implementations of the present application are disclosed in detail, specify 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 benefit 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 uses in a other embodiment, is combined with the feature in other embodiment, or the feature in substitution other embodiment.
It should be emphasized that term "comprises/comprising" refers to the presence of feature, one integral piece, step or component when using herein, but simultaneously
It is not excluded for the presence or additional of one or more other features, one integral piece, step or component.
Detailed description of the invention
Fig. 1 is a kind of step process of the quantitatively characterizing method of rock stratum internal fissure development characteristics in the application embodiment
Figure；
Fig. 2 a, which is 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 is 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 is 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 is a kind of andesite Number of Fractures cumulative frequency under depth crack frequency corresponding with single formation depth
Rate cumulative curve chart；
Fig. 3 b is andesite Number of Fractures cumulative frequency under another depth crack frequency 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 embodiment
Below in conjunction with the drawings and specific embodiments, elaborate to technical solution 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 the present invention has been read, this field skill
Art personnel each fall in the application range as defined in the appended claims the modification of various equivalent forms of the invention.
Unless otherwise defined, all technical and scientific terms used herein and the technical field for belonging to 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 it can directly on the other element when element is referred to as " being set to " another element
Or there may also be elements placed in the middle.When an element is considered as " connection " another element, it, which can be, is 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 meant to be the only embodiment.
The present invention provides the quantitatively characterizing method and device of a kind 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: being directed to specific rock type, carries out lithologic log solution 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 determine the rock type for needing to carry out lithologic interpretation for research work area target zone geological condition.
Specifically, can be merged according to Research Requirements by lithology statistical result, simplify 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 determines, 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 one layer.
Specific the step of merging upwards are as follows: in lithologic log interpretation process, the rock type of some continuous depth if it exists
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 feature in oblique crack is 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 are as follows:
In above formula: FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；∑LiIt is single
Crack total number in lithologic body, unit are 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 into single rock along depth direction
Property body fractue spacing histogram and Number of Fractures generate corresponding with the single formation depth crack of Number of Fractures cumulative frequency
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
Characterization, specifically: 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 biggish depth segment of rate.
In general, having at least one slope in the crack frequency accumulation curve.When the slope is greater 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 greater 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 greater 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, then crack is evenly distributed in layer when the crack frequency accumulation curve is one-part form；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
When for three-stage, 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 below with reference to specific crack frequency accumulation curve graph.
Fig. 2 a is 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 is 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 is two-part, i.e., when entire accumulation curve is equivalent to tool there are two types of slope, wherein slope compared with
Close the top of (relatively steep) greatly, smaller (more gentle) 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 greater 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 (heterogeneity 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 heterogeneity by the slope of curve.Slope difference is big, and heterogeneity is strong, and slope difference is small,
It is distributed relatively uniform, is finally unistage type.Certainly, according to the size of the slope, fractue spacing position can be determined, slope is larger
, the fractue spacing for corresponding to depth segment is more, and slope is lesser, and the fractue spacing for corresponding to depth segment is less.
Fig. 2 c is 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 is three-stage, i.e., when entire accumulation curve is equivalent to tool there are three types of slope, wherein slope compared with
Close the top and bottom of (relatively steep) greatly, slope is smaller (more gentle) 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 as follows: according to crack interpretation results, be spaced constant depth (generally 2m, 5m etc.) and count each interval deeply
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
Divide ratio, percentage summation is 100%.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.
It finally can be 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.
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 heterogeneity in rock stratum.For this purpose, proposing that the method for the present invention solves this
Problem is determined by the way that tracing pattern and fracture layers density size are comprehensive, 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 greater 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 close to the corresponding slope in top be greater than close to bottom 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 greater 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 in other stretchings, which can also have other slope characteristics,
The application not reinflated narration herein.
It is subsequent 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；∑LiIt is single
Crack total number in lithologic body, unit are item, determine 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:
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 the actual situation
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 development degree and distribution of the crack at the rock stratum position Nei Ge are explicitly pointed out using above-mentioned density and curve
Main distribution pattern in petrology about layer internal fissure accurately portray, and has been retouched with computer language by situation
It states, the quantitative assessment for later period 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
There is the deficiency of technology, provides reliable basis for fracture-type reservoir efficiency evaluation and reservoir reconstruction, have a extensive future.
Please refer to Fig. 4, for 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 obtains module 10, for being directed to specific rock type, is carried out 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 obtains 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；∑LiIt is single
Crack total number in lithologic body, unit are 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 is performed the steps of when being executed by the processor
Step S10: being directed to specific rock type, carries out lithologic log solution 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, usually 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: to store the device of information, such as RAM, ROM 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；Using the device of optical mode storage information, such as CD or DVD.
Certainly, there are also memories 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 (micro-) processor execute computer readable program code (such as
Software or firmware) computer-readable medium, logic gate, switch, specific integrated circuit (Application Specific
Integrated Circuit, ASIC), programmable logic controller (PLC) and the form etc. for being embedded in microcontroller.
In one embodiment, there is at least one slope in the crack frequency accumulation curve；When the 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 limit 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 an 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 all possibility for the numerical value enumerated between minimum and peak
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 cover that " 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 combined term " substantially by ... constitute " 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 to
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 embodiments 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 scope of patent protection of the invention, still must be with appended claims
Subject to the range that book is defined.
Claims (8)
1. a kind of quantitatively characterizing method of rock stratum internal fissure development characteristics characterized by comprising
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, the characteristics of fracture development include distributing position of the crack in single rock stratum, comprising: are based on the crack frequency
Accumulation curve carries out quantitatively characterizing to rock stratum internal fissure distributing position, specifically: it obtains in the crack frequency accumulation curve
Slope, distributing position of the crack in single rock stratum is determined based on the slope, wherein the crack frequency accumulation is bent
There is at least one slope in line；The fractue spacing is in the biggish depth segment of slope.
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 are as follows:
In above formula: FdFor fracture layers density, unit is item/m；H is single lithology body thickness, unit m；∑LiFor single lithology
Crack total number in body, unit are 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 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 slope of the crack frequency accumulation curve from top to bottom depth direction with same greater than predetermined value,
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 greater 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.
6. the quantitatively characterizing method of internal fissure development characteristics in rock stratum as claimed in claim 5, 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 greater 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 greater 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 greater 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.
7. a kind of quantitatively characterizing device of rock stratum internal fissure development characteristics, characterized in that it comprises:
Thickness obtains module, for being directed to specific rock type, carries out lithology survey 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 obtains 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；The characteristics of fracture development includes distributing position of the crack in single rock stratum, described
Quantitatively characterizing module is configured as: carrying out quantitative table to rock stratum internal fissure distributing position based on the crack frequency accumulation curve
Sign, specifically: the slope in the crack frequency accumulation curve is obtained, determines the crack in single rock based on the slope
Distributing position in layer, wherein there is at least one slope in the crack frequency accumulation curve；The fractue spacing is in slope
Biggish depth segment.
8. the quantitatively characterizing device of internal fissure development characteristics in rock stratum as claimed in claim 7, 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 lithology
Crack total number in body, unit are item.
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 CN108427143A (en) | 2018-08-21 |
| CN108427143B true 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 (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN108008464A (en) * | 2017-11-29 | 2018-05-08 | 中国科学院地质与地球物理研究所兰州油气资源研究中心 | Crack anisotropism quantitatively characterizing method and its system |
| CN108008467A (en) * | 2017-11-29 | 2018-05-08 | 中国科学院地质与地球物理研究所兰州油气资源研究中心 | Fractue spacing quantitatively characterizing method and its system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10809400B2 (en) * | 2015-10-27 | 2020-10-20 | Schlumberger Technology Corporation | Determining shear slowness based on a higher order formation flexural acoustic mode |
-
2018
- 2018-05-11 CN CN201810445569.7A patent/CN108427143B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN108008464A (en) * | 2017-11-29 | 2018-05-08 | 中国科学院地质与地球物理研究所兰州油气资源研究中心 | Crack anisotropism quantitatively characterizing method and its system |
| CN108008467A (en) * | 2017-11-29 | 2018-05-08 | 中国科学院地质与地球物理研究所兰州油气资源研究中心 | Fractue spacing quantitatively characterizing method and its system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108427143A (en) | 2018-08-21 |
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 | |
| US20160370499A1 (en) | Determining Geomechanics Completion Quality | |
| EP2113796A1 (en) | Method and Apparatus for Analyzing Three-Dimensional Data | |
| CN111458767B (en) | Method and system for identifying lithology based on intersection graph method | |
| CN111665560B (en) | Hydrocarbon reservoir identification method, apparatus, computer device and readable storage medium | |
| CN102749647B (en) | Quantitative method for space of holes in reservoir | |
| Zhang et al. | 3D porosity modeling of a carbonate reservoir using continuous multiple-point statistics simulation | |
| CN105019892B (en) | A method of simulation fracture hole type Reservoir Body electric logging response | |
| 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 | |
| Fernández-Ibáñez et al. | Quantitative assessment of karst pore volume in carbonate reservoirs | |
| CN112861890A (en) | Reservoir evaluation model construction method and reservoir identification method | |
| 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 | |
| CN106285661B (en) | Low-resistance oil layer identification method and device based on judgment index | |
| CN114153002A (en) | Three-dimensional geological modeling method and device for natural fracture of reservoir, electronic equipment and medium | |
| US10114135B2 (en) | System and method for optimizing seismic data analysis | |
| de Oliveira et al. | Ensemble of heterogeneous classifiers applied to lithofacies classification using logs from different wells | |
| CN113625359B (en) | Method and device for calculating oil and gas containing probability of lithology trap of tight sandstone | |
| Yu et al. | Machine Learning Empowered Well Placement in a Large Unconventional Tight Gas Field in China | |
| Ateeq et al. | An integrated study to characterize and model natural fracture networks of gas condensate carbonate reservoirs, onshore Abu Dhabi, UAE | |
| CN113589370B (en) | Structure identification method and device under high-speed stratum coverage | |
| Boisvert | Mineral deposit modeling with pseudo-genetically constructed training images |
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 |