CN104977617B - Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods - Google Patents
Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods Download PDFInfo
- Publication number
- CN104977617B CN104977617B CN201410130972.2A CN201410130972A CN104977617B CN 104977617 B CN104977617 B CN 104977617B CN 201410130972 A CN201410130972 A CN 201410130972A CN 104977617 B CN104977617 B CN 104977617B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- mrow
- analyzed
- interval
- curve
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of Reservoir Fracture recognition methods and a kind of imaging logging Reservoir Fracture recognition methods, wherein the Reservoir Fracture recognition methods includes:Porosity curve determines step, based on the lithological information of reservoir to be analyzed, calculates acoustic porosity curve, the density porosity curve of reservoir to be analyzed respectively according to default porosity curve model;Porosity curve combination determines step, and in the dried layer of reservoir to be analyzed, by acoustic porosity curve and density porosity curve co-insides, the porosity curve as reservoir to be analyzed combines;Types of fractures identification step, in other intervals of reservoir to be analyzed, compare the acoustic porosity curve and density porosity curve in porosity curve combination, determine the types of fractures of corresponding interval according to default fractured model based on comparative result.The present invention can intuitively, quickly and reliably carry out the identification of Reservoir Fracture type, can effectively solve the problem that multiresolution issue present in existing method and leakage interpretation problems.
Description
Technical field
The present invention relates to oil-gas exploration and development technical field, specifically, be related to a kind of Reservoir Fracture recognition methods and into
As logging reservoir crack identification method.
Background technology
With the utilization of continually developing of world oil natural gas resource, the direction of oil-gas exploration and development is gradually by conventional gas and oil
Hide and turn to special oil and gas pools.Wherein, the exploration and development of crack elimination oil and gas reserves is particularly important, and Reservoir Fracture type
Identification be then fracture-type reservoir well log interpretation evaluation basis and emphasis.
Existing Reservoir Fracture type identification technology is typically to be split according to dual laterolog or tri-porosity logging
Seam identification.However, when carrying out crack identification using dual laterolog, the factor such as fluidic response easily causes bilaterally to produce difference
It is different, so as to cause more solution risks be present using dual laterolog identification crack.And carry out crack identification using tri-porosity logging
When, the change of the resolution ratio, lithology of instrument, the change of properties of fluid in bearing stratum may all cause three porosity to increase so that only from
The response characteristic of porosity curve, which is set out, to be difficult to accurately identify crack.
Imaging Logging System comes out from early 1990s, belongs to the logging method for being capable of direct detection crack attribute.
Image Logging Data can intuitively, image, clearly show the geologic feature of borehole wall two-dimensional space.But imaging logging exists
Lack in there is also experimental data, pattern base information is not abundant enough, geologic feature identification has the problems such as multi-solution.Example
Such as, when carrying out crack identification using imaging logging, half filling seam is often explained because having disguise in response to unobvious by leakage;It is low
Angle crack is not easy to be identified by imaging logging again because close with stratification response.
In addition, the relatively advanced crack identification technology such as under-mine TV, stratigraphic dip is have also appeared in recent years, but it is higher
Cost limits widely using for they.
Based on the above situation, a kind of inexpensive, method that can effectively identify various Reservoir Fracture types is needed badly.
The content of the invention
To solve the above problems, the invention provides a kind of Reservoir Fracture recognition methods, methods described includes:
Porosity curve determines step, based on the lithological information of reservoir to be analyzed, according to default porosity curve model point
Acoustic porosity curve, the density porosity curve of the reservoir to be analyzed are not calculated;
Porosity curve combination determines step, in the dried layer of the reservoir to be analyzed, by acoustic porosity curve and density
Porosity curve overlaps, and the porosity curve as the reservoir to be analyzed combines;
Types of fractures identification step, in other intervals of the reservoir to be analyzed, the porosity curve combination
The acoustic porosity curve and density porosity curve, corresponding interval is determined according to default fractured model based on comparative result
Types of fractures.
According to one embodiment of present invention, before the porosity curve determines step, methods described also includes:
Lithological information obtaining step, according to the Core information of reservoir to be analyzed, geological analysis information and log information,
Determine the lithological information of reservoir to be analyzed.
According to one embodiment of present invention, in the default porosity curve model:
According to the lithological information of the reservoir to be analyzed, the rock matrix interval transit time and rock of the reservoir to be analyzed are determined
Stone skeletal density;
Based on rock matrix interval transit time and matrix density, the interval transit time curve and volume of reservoir to be analyzed are utilized
Density curve, the sound of the reservoir to be analyzed is calculated according to acoustic porosity response model and density porosity response model respectively
Porosity curve and density porosity curve.
According to one embodiment of present invention, the acoustic porosity response model includes:
Wherein, φACThe acoustic porosity of interval is analyzed in expression, and Δ t represents to analyze the interval transit time of interval, Δ tma
The rock matrix interval transit time of interval, Δ t are analyzed in expressionfThe interval transit time of the blowhole fluid of interval is analyzed in expression,
VshThe shale content of interval, Δ t are analyzed in expressionshShale interval transit time is represented, k represents compaction correction coefficient.
According to one embodiment of present invention, the density porosity response model includes:
Wherein, φDENThe density porosity of interval is analyzed in expression, and ρ represents to analyze the rock volume density of interval, ρma
The matrix density of interval, ρ are analyzed in expressionfThe blowhole fluid density of interval is analyzed in expression, and n represents Shale Correction
Coefficient, VshThe shale content of interval, ρ are analyzed in expressionshThe shale bulk density of interval is analyzed in expression.
According to one embodiment of present invention,
Porosity curve determines that step also includes:
Based on the neutron curve of the reservoir to be analyzed, the neutron porosity curve of the reservoir to be analyzed is determined;
In the porosity curve combines determination step:
In the dried layer of the reservoir to be analyzed, using the neutron porosity curve as reference curve, by acoustic porosity
Curve, density porosity curve and the neutron porosity curve co-insides, the porosity curve group as the reservoir to be analyzed
Close.
According to one embodiment of present invention, before the porosity curve combines determination step, methods described is also wrapped
Include dried layer and determine step, the dried layer determines that step includes:
The stratum micro resistor of reservoir to be analyzed is obtained, the region that bright white is shown as in imaging is defined as
The dried layer of the reservoir to be analyzed;And/or
The bilaterally curve of reservoir to be analyzed is obtained, the region that high resistant is shown as in the bilaterally curve and is overlapped is true
It is set to the dried layer of the reservoir to be analyzed.
According to one embodiment of present invention, in the types of fractures identification step:
When the density porosity curve is higher than the acoustic porosity curve, judge the interval development for high angle
Crack;
When the acoustic porosity curve is higher than the density porosity curve, judge the interval development for low angle
Crack;
When the acoustic porosity curve and the density porosity curve intersect, judge interval development is
Chicken-wire cracking.
According to one embodiment of present invention, after the types of fractures identification step, methods described also includes checking
Step, the verification step include:
Depth playback is carried out to the rock core of the reservoir to be analyzed of acquisition, judges each of the rock core after depth playback
Whether the types of fractures of corresponding interval of the types of fractures of individual interval to identifying to obtain in the crack identification step matches;With/
Or,
The types of fractures of each interval of reservoir to be analyzed is identified using imaging logging, judges the imaging logging identification
To each interval types of fractures and the types of fractures of corresponding interval that identifies to obtain in the crack identification step whether
Match somebody with somebody.
Present invention also offers a kind of imaging logging Reservoir Fracture recognition methods, the described method comprises the following steps:
Imaging logging identification step, the crack in reservoir to be identified is identified using imaging logging, obtains described treat
Identify the type and its geometric parameter in crack in reservoir;
Conventional logging identification step, crack in the reservoir to be identified is identified using method as described above,
Obtain the type in crack in the reservoir to be identified;
Correct step, according to the type in crack in the reservoir to be identified obtained in the conventional logging identification step with
The type in the crack obtained in the imaging logging identification step, judge in the imaging logging identification step with the presence or absence of leakage solution
The crack released, to correct the recognition result of the imaging logging identification step.
Crack identification method provided by the invention uses conventional logging, is recalculated by conventional tri-porosity logging curve
Arrangement, to carry out the identification of Reservoir Fracture type.Compared to existing imaging logging, by the present invention in that with conventional logging,
While the accuracy of recognition result is effectively ensured, additionally it is possible to effectively reduce cost, also make it that implementation process is more directly perceived, fast
Speed, operability are stronger.Meanwhile rearranging by three porosity curve of the invention, efficiently solve sharp in existing method
With the three porosity response characteristic and bilaterally multiresolution issue in positive and negative difference identification crack.In addition, will be provided by the present invention
Crack identification method is combined with imaging logging, additionally it is possible to effectively overcomes existing imaging logging because half filling seam response is not clear
The problem of explaining is leaked in aobvious, low angle crack and stratification the response crack caused by, and this is established extensively for imaging logging
Basis.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
Obtain it is clear that or being understood by implementing the present invention.The purpose of the present invention and other advantages can be by specification, rights
Specifically noted structure is realized and obtained in claim and accompanying drawing.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the accompanying drawing required in technology description to do simple introduction:
Fig. 1 is the flow chart of Reservoir Fracture recognition methods according to an embodiment of the invention;
Fig. 2 is the flow chart of Reservoir Fracture recognition methods in accordance with another embodiment of the present invention;
Fig. 3 is certain fine and close Fractured Gas Reservoir X1 wells low angle frature identification figure according to an embodiment of the invention;
Fig. 4 is certain fine and close Fractured Gas Reservoir X2 wells high angle frature identification figure according to an embodiment of the invention;
Fig. 5 is certain fine and close Fractured Gas Reservoir X3 well patterns shape seam identification figure according to an embodiment of the invention;
Fig. 6 is the flow chart of the Reservoir Fracture recognition methods according to an embodiment of the invention based on imaging logging;
Fig. 7 is certain fine and close filling irrigating sealant identification figure of Fractured Gas Reservoir X4 wells half according to an embodiment of the invention.
Embodiment
Embodiments of the present invention are described in detail below with reference to drawings and Examples, and how the present invention is applied whereby
Technological means solves technical problem, and the implementation process for reaching technique effect can fully understand and implement according to this.Need to illustrate
As long as not forming conflict, each embodiment in the present invention and each feature in each embodiment can be combined with each other,
The technical scheme formed is within protection scope of the present invention.
Embodiment one:
For dual laterolog, in fracture development section, due to the intrusion of drilling fluid and drilling fluid filtrate, the resistance at crack
Rate shows as the relatively low resistivity value in high resistant background than the obvious reduction of matrix resistivity of densification.The heterogeneity in crack
Property make resistivity curve form be in height between mutual, up-and-down more needle patterns.
Generally, the al-lateral resistivity value at crack decreases, and its reduction amplitude and development degree of micro cracks in oil, splits
The fluid properties stitched in aperture and crack are relevant.In the fracture development section of oily, al-lateral resistivity difference becomes big, curve
In double track phenomenon.
In high angle fracture, resistivity curve is serrated, and deep lateral resistivity change is little, shallow lateral resistivity by
In being influenceed substantially to reduce by crack, the lateral resistivity of the depth forms larger positive separation(I.e. deep lateral resistivity is more than
Shallow lateral resistivity).When low angle crack in stratum be present, the lateral resistivity of the depth substantially reduces, and shows therebetween
For indifference XOR negative variance, i.e., shallow lateral resistivity is more than or equal to deep lateral resistivity.
For tri-porosity logging, the occurrence in crack around fracture development interval, interval transit time curve change and pit shaft
It is relevant with development degree.Sound wave can bypass crack as far as possible by shortest time selection sound path in propagation.Therefore, sound wave is to high angle
Crack reflection is insensitive.
Stitched for horizontal fracture and oblique, acoustic wave propagation path is orthogonal to that, therefore interval transit time can reflect that level is split
Seam and oblique seam, curve is in spiculation.When running into big horizontal joint, interval transit time significantly increases, it is possible to cycle occurs
Jump.Typically in fracture development interval, for low angle frature crack, density value reduces, neutron porosity increase;For the angle of elevation
For degree seam, density curve value substantially reduces, and density porosity increase, amplitude difference, neutron porosity occurs with neutron porosity
Increase.
It is dual laterolog and tri-porosity logging as described above mostly using Conventional Logs in the prior art
In crack, the response characteristic of section carries out the crack identification of reservoir, and its particular technique includes:
(1)Identification and evaluation fracture reservoir are combined with Image Logging Data using Conventional Logs.The main root of the technology
Logged well according to al-lateral resistivity, acoustic logging, density log, natural gamma ray spectrometry log and dipmeter logging be in chicken-wire cracking, height
Different response characteristics on angle crack and low angle crack, in conjunction with stratum micro resistor(Formaiton
MicroScanner Image, referred to as FMI), to carry out crack identification.
(2)Crack identification index is defined, using imaging and conventional data, establishes the one of different types of fractures and conventional data
One corresponding data sample, establish the cross plot such as fracture index and resistivity, three porosity curve, gamma-ray spectrometry.
(3)Research is identified to research area's fracture reservoir development degree using combined chance method.The technology is ground first
Study carefully the conventional method that fracture reservoir is identified according to well-log information, including it is logging anomalies ratio method, rock pore structure index method, full
With degree ratio method and porosity ratio method.The subsequent technology sees the various features parameter curve that above method is calculated with rock core
The data such as survey, geological logging and formation testing are contrasted, and analyze the ability in each characteristic parameter reflection crack, and determine its weighting system
Number, draws the overall target for reflecting whether crack and its development degree be present.
(4)The identification and evaluation in crack is carried out based on Conventional Logs based on BP neural network research.This technology is for grinding
Study carefully the concrete condition in area, selectively choose interval transit time, compensation density, compensated neutron, microballoon, the depth of reflection FRACTURE CHARACTERISTICS
It is shallow laterally to wait curve, sample mode is established, curve abnormality value is extracted, is normalized, establishes neural network recognization model.
Finally the identification of Reservoir Fracture is carried out using the neural network recognization model established.
In crack identification method described above, the response characteristic or three holes of triangle porosity curve fracture are utilized
Degree difference/ratio is to carry out the most commonly used step.But due to causing the factor that three porosity increases a lot, only from porosity
The response characteristic of curve set out be difficult to standard go identify crack.
Meanwhile if using the identification that Reservoir Fracture is carried out the methods of under-mine TV, stratigraphic dip, imaging logging, deposit again
Cost is too high, experimental data lacks the shortcomings of.
The present invention is according to the defects of above-mentioned crack identification method, it is proposed that a kind of Reservoir Fracture based on Logging Curves
Recognition methods, Fig. 1 show the flow chart of this method.
As shown in figure 1, the lithological information of reservoir to be analyzed is determined in step S101 first.In the present embodiment, based on normal
Core information, geological analysis information and the log information for the reservoir to be analyzed that rule well logging is got, by combining three kinds of letters
Cease to obtain the lithological information of reservoir to be analyzed, wherein log includes interval transit time curve, bulk density curve and neutron
Curve.It should be noted that in other embodiments of the invention, reservoir to be analyzed can also be obtained by other rational methods
Lithological information, the invention is not restricted to this.
Then perform porosity curve and determine step, it is bent according to default porosity based on the lithological information of reservoir to be analyzed
Line model calculates the acoustic porosity curve of the reservoir to be analyzed, density porosity curve respectively.
As shown in figure 1, in the present embodiment, default porosity curve model includes step S102 and step S103.In step
In S102, according to the lithological information of reservoir to be analyzed, by search rock essential mineral sound wave, matrix density method come
It is determined that the rock matrix interval transit time corresponding with the lithology of reservoir to be analyzed and matrix density.
In step s 103, based on the rock matrix interval transit time and matrix density obtained in step S102, difference
Using the interval transit time curve and bulk density curve of reservoir to be analyzed, according to acoustic porosity response model and density porosity
Response model, calculate the acoustic porosity curve and density porosity curve of reservoir to be analyzed.So far it is true to complete porosity curve
Determine step.
In the present embodiment, acoustic porosity response model can use formula(1)Represent:
Wherein, φACThe acoustic porosity of interval is analyzed in expression, and Δ t represents to analyze the interval transit time of interval, Δ tma
The rock matrix interval transit time of interval, Δ t are analyzed in expressionfThe interval transit time of the blowhole fluid of interval is analyzed in expression,
VshThe shale content of interval, Δ t are analyzed in expressionshShale interval transit time is represented, k represents compaction correction coefficient.
Density porosity response model can use formula(2)Represent:
Wherein, φDENThe density porosity of interval is analyzed in expression, and ρ represents to analyze the rock volume density of interval, ρma
The matrix density of interval, ρ are analyzed in expressionfThe blowhole fluid density of interval is analyzed in expression, and n represents Shale Correction
Coefficient, VshThe shale content of interval, ρ are analyzed in expressionshThe shale bulk density of interval is analyzed in expression.
It should be noted that in other embodiments of the invention, acoustic porosity response model and density porosity ring
Model is answered all to use other reasonable function statements, the invention is not restricted to this.
In order to solve to respond form and bilaterally positive and negative difference hardly possible using three porosity curve in existing crack identification method
The problem of to accurately identify Reservoir Fracture type, the present invention according to conventional logging three porosity curve by rearranging combination
To overcome disadvantages mentioned above.In the present embodiment, the combination that rearranges of porosity curve refers to determine porosity curve in step
Obtained acoustic porosity curve and density porosity curve overlaps in the dried layer of reservoir to be analyzed.
Again as shown in figure 1, determining the dried layer of reservoir to be analyzed in step S104.In the present embodiment, based on to be analyzed
The FMI of reservoir is imaged to determine dried layer.In the method, the FMI imagings of reservoir to be analyzed are obtained first, are then imaged FMI
In show as the region of bright white and be defined as the dried layer of reservoir to be analyzed.
It should be noted that in other embodiments of the invention, can also determine to treat point using other rational methods
The dried layer of reservoir is analysed, such as the dried layer of the reservoir is determined according to the bilaterally curve of reservoir to be analyzed.In the method, first
The bilaterally curve of reservoir to be analyzed is obtained, the region that high resistant is shown as in bilaterally curve and is overlapped then is defined as the storage
The dried layer of layer.
After the dried layer of reservoir to be analyzed is obtained, perform porosity curve combination and determine step S105, in reservoir to be analyzed
Dried layer, by acoustic porosity curve and density porosity curve co-insides, and using the curve combination after coincidence as storage to be analyzed
The porosity curve combination of layer.This makes it possible to same interval acoustic porosity and density porosity are intuitively presented for user
State, provided convenience for the further analysis of data.
After porosity curve is tortuous to be determined, types of fractures identification step S106, the porosity overlapped with dried layer are performed
Curve combination is standard, in other intervals of reservoir to be analyzed, compare in porosity curve combination acoustic porosity curve with it is close
Porosity curve is spent, and is based on comparative result, the crack class of corresponding interval in reservoir to be analyzed is determined according to default fractured model
Type.
In this implementation, default fractured model includes:When density porosity curve is higher than acoustic porosity curve, judging should
Interval development for high angle fracture;When acoustic porosity curve is higher than density porosity curve, interval development is judged
For low angle crack;When acoustic porosity curve and the density porosity curve height each other, interval development is judged
For chicken-wire cracking.
By comparing acoustic porosity curve corresponding to a certain interval and close in porosity curve combination in step s 106
Spend the height situation of porosity curve, you can identify the type in crack in the interval.
In order to further determine that the recognition effect of Reservoir Fracture recognition methods that the present embodiment provided, in the present embodiment,
It further comprises verification step S107.The rock core of reservoir to be analyzed is obtained first, by the rock core for observing reservoir to be analyzed, you can obtain
Obtain the types of fractures of each interval of reservoir to be analyzed.Depth playback then is carried out to the rock core got, makes reservoir to be analyzed
Each interval is corresponding with each interval in step S106 in rock core, just intentional so when being compared to the two
Justice.Finally, the types of fractures and step of each interval of the rock core of the reservoir to be analyzed after being playbacked by the depth for judging to observe
Whether the types of fractures of the corresponding interval recognized in rapid S106 is matched, and the knowledge to method fracture provided by the invention is carried out with this
Sorrow of separation condition is verified.
It should be noted that in other embodiments in accordance with the invention, can also be using other rational methods come to this
The recognition methods that invention provides is verified.Such as verified using imaging logging, treated first with imaging logging identification
The types of fractures of each interval of reservoir is analyzed, subsequently determines whether the types of fractures and step of each interval that imaging logging recognizes
Whether the types of fractures for the corresponding interval for identifying to obtain in S106 matches.
The present invention is calculated by conventional logging three porosity curve, rearranges combination, is formd a kind of identification reservoir and is split
The method for stitching type, solve and respond form and bilaterally positive and negative difference is difficult to accurately identify storage with three porosity curve merely
The problem of layer types of fractures.
This method successively practices in Sichuan Basin gas field, Ordos Basin oil gas field, above-mentioned two oil gas
The most of drilled well in field only measures Logging Curves, and only a small amount of well has rock core or Image Logging Data.Rock is respectively adopted
Crack evidence known to the heart or imaging logging can the FRACTURE CHARACTERISTICS of real example method provided by the invention is verified, it was demonstrated that this
Method is not only accurately and reliably, and the half filling seam for being difficult to accurate judgement for imaging logging can also be accomplished to accurately identify.Therefore,
Method provided by the present invention accomplishes that only just energy is convenient, fast and height is recognized accurately with most basic Conventional Logs
Angle seam, low angle frature and netted seam etc., judge to provide important evidence for the production capacity of fine and close fracture reservoir.Using the hair
It is bright, fine and close slit formation oilfield prospecting developing cost can be reduced, improves the exploration and development efficiency of oil gas field, improves oil gas harvesting
Rate.
Embodiment two:
When containing special mineral in reservoir to be analyzed, the different of acoustic porosity and density porosity curve may be caused
Often mutation.For the problem, neutron porosity curve is introduced in the present embodiment, by the comprehensive descision of three porosity curve,
To avoid merely because density porosity curve is higher than acoustic porosity curve, or acoustic porosity curve higher than density porosity song
Line, and cause the erroneous judgement of types of fractures.
Fig. 2 shows the flow chart for the crack identification method that neutron porosity curve is introduced in the present embodiment.
As shown in Fig. 2 in the present embodiment, the lithological information of reservoir to be analyzed is determined in step s 201 first, itself and step
Rapid S101 principle is identical, will not be repeated here
In step S202, according to the lithological information of reservoir to be analyzed, by searching rock essential mineral sound wave, lithosome
Product density and the method for compensated neutron determine the rock matrix interval transit time corresponding with the lithology of reservoir to be analyzed and rock
Skeletal density.
In step S203, based on the rock matrix interval transit time and matrix density obtained in step S202, utilize
Interval transit time curve, bulk density curve and the neutron curve of reservoir to be analyzed, respectively according to corresponding porosity response model,
Calculate the acoustic porosity curve, density porosity curve and neutron porosity curve of reservoir to be analyzed.In the present embodiment, sound wave
The Computing Principle of porosity and density porosity is identical with the step S102 in embodiment one, will not be repeated here.For neutron
Porosity curve, its neutron curve are neutron porosity curve.
The dried layer of reservoir to be analyzed, its principle and the step S104 phases in embodiment one are then determined in step S204
Together, will not be repeated here.In step S205, in the dried layer of reservoir to be analyzed, using neutron porosity curve as reference curve,
By acoustic porosity curve, density porosity curve and neutron porosity curve co-insides, using the curve after coincidence as to be analyzed
The porosity curve combination of reservoir.
The identification of Reservoir Fracture type to be analyzed is then carried out in step S206, in the present embodiment in step S207
Recognition result is verified that its principle is identical with the step S106 in embodiment one and step S107, will not be repeated here.
Crack identification, figure are carried out using the crack identification method provided in the present embodiment Fractured Gas Reservoir X1 wells fine and close to one
3 show the identification figure of gas reservoir X1 wells.
As shown in figure 3, the in log the 1st shows gas reservoir X1 natural potential(Spontaneous
Poential, referred to as SP)Curve and natural gamma(Natural Gamma-ray, referred to as GR), the 2nd shows sound wave
The time difference(Acoustic, referred to as AC)Curve, rock volume density(Density, referred to as DEN)Curve and compensated neutron
(CNL)Curve.3rd shows that porosity curve combines, and sound wave curve and density curve are correspondingly converted into sound wave hole by it
Porosity curve and density porosity curve, and this two curves and neutron curve, i.e. neutron porosity curve are placed on together together
In one curve road.
From Fig. 3 in log the 4th in correlation curve can be seen that resistance near 4090 meters of depth
The resistivity that rate curve RXO and resistivity curve RT are characterized is very high, imaging logging FMI display colors near the depth
For highlighted white, it is nearby dried layer to confirm the depth.By acoustic porosity curve, density porosity curve and neutron porosity
Curve is in the dried layer determined(I.e. near 4090 meters of depth)Overlap, obtain the porosity curve combination of gas reservoir X1 wells.
As shown in figure 3, porosity curve combination, in 4085~4089 meters of depth boundses, acoustic porosity curve is higher than close
Porosity curve is spent, by crack identification method provided by the invention, judges that the depth segment has low angle frature development.From Fig. 3
In imaging logging map in as can be seen that 4086~4089 meters of depth segments low angle frature development really be present, with the present invention provide
Method crack identification result it is identical.
Equally near 4091 meters of depth segments, acoustic porosity curve also apparently higher than density porosity curve, passes through this
Invent the crack identification method provided, it can be determined that going out the depth, nearby there is also low angle frature development.Surveyed from the imaging in Fig. 3
Nearby low angle frature development really be present in well it can be seen from the figure that, 4091 meters of depth segments.
It can thus be seen that identification of the crack identification method provided by the invention for low angle frature is reliable, accurate.
Meanwhile also carry out crack using the crack identification method provided in the present embodiment Fractured Gas Reservoir X2 wells fine and close to one
Identification, Fig. 4 show the identification figure of gas reservoir X2 wells.
From fig. 4, it can be seen that near 3959 meters of depth, the resistivity curve RXO in the 4th in log
And although the resistivity that resistivity curve RT is characterized is not very high, three articles of porosity curves in the 3rd are at this
Depth point is all very small nearby, and this again shows that the depth is nearby dried layer.So by acoustic porosity curve, density porosity
Curve and neutron porosity curve it is determined that dried layer(I.e. near 3959 meters of depth)Place overlaps, and obtains the porosity of gas reservoir X2 wells
Curve combination.
As shown in figure 4, porosity curve combination, in 3954~3957 meters of depth boundses, density porosity curve is higher than sound
Porosity curve.According to crack identification method provided by the invention, judge that the depth segment has high angle frature development.Further according to
Imaging logging map comprehensive analysis can be seen that 3954~3957 meters of depth segments and high angle frature development really be present, be carried with the present invention
The crack identification result of the method for confession is identical.
It can thus be seen that the high angle in reservoir can be also recognized accurately in crack identification method provided by the invention
Seam.
Meanwhile in order to verify that the crack identification method of the present embodiment offer to the recognition effect of netted seam, now utilizes the party
Method Fractured Gas Reservoir X3 wells fine and close to one carry out crack identification, and Fig. 5 shows the identification figure of gas reservoir X3 wells.
From fig. 5, it can be seen that near 4025 meters of depth, the resistivity curve RXO in the 4th in log
And the resistivity ratio that resistivity curve RT is characterized is higher, and three articles of porosity curves in the 3rd are attached in the depth point
Near all very small, it is nearby dried layer to show the depth.So by acoustic porosity curve, density porosity curve and neutron hole
Write music line it is determined that dried layer(I.e. near 4025 meters of depth)Place overlaps, and obtains the porosity curve combination of gas reservoir X3 wells.
As shown in figure 5, porosity curve is combined in 4035~4039 meters of depth boundses, density porosity curve and sound wave
Porosity curve intersects, i.e. two curves just, according to crack identification method provided by the invention, judge the depth each other
Spend the netted seam development of section.4035~4039 meters of depth segments, which are can be seen that, further according to imaging logging map comprehensive analysis net really be present
Shape seam development, it is identical with the crack identification result of method provided by the invention.
Equally in 4030~4031 meters of depth segments, density porosity curve is handed over acoustic porosity curve there is also mutual
The situation of fork, pass through crack identification method provided by the invention, it can be determined that going out the depth, nearby there is also netted seam development.From
As can be seen that really there are netted seam in 4030~4031 meters of depth segments in imaging logging map.
Thus see to find out, there is also have in terms of the identification of netted seam is carried out for crack identification method provided by the invention
Effect.
Embodiment three:
Imaging logging has the characteristics of can identifying to obtain Reservoir Fracture type and fracture geometry parameter, and it is needing to obtain
The occasion of fracture geometry parameter has extensive use.But because half filling seam response unobvious, part low angle frature and stratification are rung
Should approach etc. reason, existing imaging logging to above-mentioned crack exist leakage explain the problem of, this largely constrains imaging
The application of well logging.For the problem, present invention also offers a kind of Reservoir Fracture recognition methods based on imaging logging.
It was found from from foregoing description, the identification of the Reservoir Fracture based on conventional logging described in embodiment one and embodiment two
Method can have for the various types of fractures in reservoir efficiently and accurately to be identified.So in the present embodiment, routinely surveyed using this
Whether well method is modified to the imaging logging map that imaging logging obtains, determine to have in imaging logging map the leakage explain with this
The half filling crack and low angle crack similar to stratification response, so as to correct the recognition result of imaging logging.
Fig. 6 shows the flow chart of the Reservoir Fracture recognition methods based on imaging logging in the present embodiment.
As shown in fig. 6, in imaging logging identification step S601, the crack in reservoir to be identified is entered using imaging logging
Row identification, obtains the type in the crack of each interval in reservoir to be identified.
In conventional logging identification step S602, using being provided in such as embodiment one or embodiment two based on conventional logging
Crack in reservoir to be identified is identified for the Reservoir Fracture recognition methods of data, obtains each interval in reservoir to be identified
The type in crack.
In the positive step S603 of last repairing, crack in the reservoir to be identified obtained according to S602 in conventional logging identification step
Type and imaging logging identification step S601 in the obtained types of fractures of corresponding interval, judge imaging logging identification step
The crack explained in S601 with the presence or absence of leakage.If there is leakage explain crack, then the types of fractures in step S603 come
The crack of the corresponding interval of additional explanation.
Fig. 7 shows to carry out crack knowledge to fine and close Fractured Gas Reservoir X4 wells using the crack identification method provided in the present embodiment
Other identification figure.
From figure 7 it can be seen that near 3798 meters of depth, the resistivity curve RXO in the 4th in log
And the resistivity ratio that resistivity curve RT is characterized is higher, and three articles of porosity curves in the 3rd are attached in the depth point
Near all very small, it is nearby dried layer to show the depth.So by acoustic porosity curve, density porosity curve and neutron hole
Write music line it is determined that dried layer(I.e. near 3798 meters of depth)Place overlaps, and obtains the porosity curve combination of gas reservoir X3 wells.
As shown in fig. 7, porosity curve combination, in 3784~3786 meters of depth boundses, density porosity curve is higher than sound
Porosity curve, the crack identification method provided according to embodiment one or embodiment two, judges height in the depth segment be present
Angle seam development.If the high angle frature fully opened, can be easier to identify from imaging logging map.If split
Seam is half filling, then more difficult from imaging logging map to identify.And in the imaging logging map obtained using imaging logging
In the 2nd in be difficult to identify that in the depth segment there is high angle frature development.So it is combined with embodiment one or embodiment two
The high angle fracture that the crack identification method of offer identifies, so as to identify leakage explanation in the 2nd in imaging logging map
Half filling high angle frature.With reference in the 1st in imaging logging map as can be seen that figure in show half filling high angle frature with this
The half filling high angle frature that method identifies is consistent, so the method recognition result that the present embodiment is provided is accurate.
In addition, the Reservoir Fracture recognition methods provided by conjunction with the embodiments one or embodiment two, using imaging logging also
The low angle frature for leaking and explaining because similar to stratification response, its principle and the above-mentioned knowledge for half-and-half filling high angle frature can be identified
It is not identical, it will not be repeated here.
Thus see to find out, the present embodiment provide the crack identification method based on imaging logging effectively overcome it is existing into
As well logging leaks the problem of explaining to the low angle frature for responding unconspicuous half filling seam and being close with stratification response so that imaging
The reliability of well logging has obtained further raising, contributes to promoting the use of for imaging logging.
Although disclosed herein embodiment as above, described content only to facilitate understand the present invention and adopt
Embodiment, it is not limited to the present invention.Any those skilled in the art to which this invention pertains, this is not being departed from
On the premise of the disclosed spirit and scope of invention, any modification and change can be made in the implementing form and in details,
But the scope of patent protection of the present invention, still should be subject to the scope of the claims as defined in the appended claims.
Claims (22)
1. a kind of Reservoir Fracture recognition methods, it is characterised in that methods described includes:
Porosity curve determines step, based on the lithological information of reservoir to be analyzed, is counted respectively according to default porosity curve model
Calculate acoustic porosity curve, the density porosity curve of the reservoir to be analyzed;
Porosity curve combination determines step, in the dried layer of the reservoir to be analyzed, by acoustic porosity curve and density hole
Curve co-insides are spent, the porosity curve as the reservoir to be analyzed combines;
Types of fractures identification step, the institute in other intervals of the reservoir to be analyzed, the porosity curve combination
Acoustic porosity curve and density porosity curve are stated, splitting for corresponding interval is determined according to default fractured model based on comparative result
Stitch type.
2. the method as described in claim 1, it is characterised in that before the porosity curve determines step, methods described
Also include:
Lithological information obtaining step, according to the Core information of reservoir to be analyzed, geological analysis information and log information, it is determined that
The lithological information of reservoir to be analyzed.
3. the method as described in claim 1, it is characterised in that in the default porosity curve model:
According to the lithological information of the reservoir to be analyzed, the rock matrix interval transit time and rock bone of the reservoir to be analyzed are determined
Frame density;
Based on rock matrix interval transit time and matrix density, the interval transit time curve and bulk density of reservoir to be analyzed are utilized
Curve, the sound wave hole of the reservoir to be analyzed is calculated according to acoustic porosity response model and density porosity response model respectively
Porosity curve and density porosity curve.
4. method as claimed in claim 2, it is characterised in that in the default porosity curve model:
According to the lithological information of the reservoir to be analyzed, the rock matrix interval transit time and rock bone of the reservoir to be analyzed are determined
Frame density;
Based on rock matrix interval transit time and matrix density, the interval transit time curve and bulk density of reservoir to be analyzed are utilized
Curve, the sound wave hole of the reservoir to be analyzed is calculated according to acoustic porosity response model and density porosity response model respectively
Porosity curve and density porosity curve.
5. method as claimed in claim 3, it is characterised in that the acoustic porosity response model includes:
<mrow>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
Wherein, φACThe acoustic porosity of interval is analyzed in expression, and Δ t represents to analyze the interval transit time of interval, Δ tmaRepresent institute
Analyze the rock matrix interval transit time of interval, Δ tfThe interval transit time of the blowhole fluid of interval, V are analyzed in expressionshRepresent
The shale content of analyzed interval, Δ tshShale interval transit time is represented, k represents compaction correction coefficient.
6. method as claimed in claim 4, it is characterised in that the acoustic porosity response model includes:
<mrow>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
Wherein, φACThe acoustic porosity of interval is analyzed in expression, and Δ t represents to analyze the interval transit time of interval, Δ tmaRepresent institute
Analyze the rock matrix interval transit time of interval, Δ tfThe interval transit time of the blowhole fluid of interval, V are analyzed in expressionshRepresent
The shale content of analyzed interval, Δ tshShale interval transit time is represented, k represents compaction correction coefficient.
7. the method as any one of claim 3~6, it is characterised in that the density porosity response model includes:
<mrow>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>N</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&rho;</mi>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
Wherein, φDENThe density porosity of interval is analyzed in expression, and ρ represents to analyze the rock volume density of interval, ρmaRepresent
The matrix density of analyzed interval, ρfThe blowhole fluid density of interval is analyzed in expression, and n represents Shale Correction system
Number, VshThe shale content of interval, ρ are analyzed in expressionshThe shale bulk density of interval is analyzed in expression.
8. such as method according to any one of claims 1 to 6, it is characterised in that
Porosity curve determines that step also includes:
Based on the neutron curve of the reservoir to be analyzed, the neutron porosity curve of the reservoir to be analyzed is determined;
In the porosity curve combines determination step:
In the dried layer of the reservoir to be analyzed, using the neutron porosity curve as reference curve, by acoustic porosity curve,
Density porosity curve and the neutron porosity curve co-insides, the porosity curve as the reservoir to be analyzed combine.
9. method as claimed in claim 7, it is characterised in that
Porosity curve determines that step also includes:
Based on the neutron curve of the reservoir to be analyzed, the neutron porosity curve of the reservoir to be analyzed is determined;
In the porosity curve combines determination step:
In the dried layer of the reservoir to be analyzed, using the neutron porosity curve as reference curve, by acoustic porosity curve,
Density porosity curve and the neutron porosity curve co-insides, the porosity curve as the reservoir to be analyzed combine.
10. such as the method any one of claim 1~6,9, it is characterised in that combine and determine in the porosity curve
Before step, methods described also determines step including dried layer, and the dried layer determines that step includes:
The stratum micro resistor of reservoir to be analyzed is obtained, the region that bright white is shown as in imaging is defined as described
The dried layer of reservoir to be analyzed;And/or
The bilaterally curve of reservoir to be analyzed is obtained, the region that high resistant is shown as in the bilaterally curve and is overlapped is defined as
The dried layer of the reservoir to be analyzed.
11. method as claimed in claim 7, it is characterised in that described before the porosity curve combines determination step
Method also determines step including dried layer, and the dried layer determines that step includes:
The stratum micro resistor of reservoir to be analyzed is obtained, the region that bright white is shown as in imaging is defined as described
The dried layer of reservoir to be analyzed;And/or
The bilaterally curve of reservoir to be analyzed is obtained, the region that high resistant is shown as in the bilaterally curve and is overlapped is defined as
The dried layer of the reservoir to be analyzed.
12. method as claimed in claim 8, it is characterised in that described before the porosity curve combines determination step
Method also determines step including dried layer, and the dried layer determines that step includes:
The stratum micro resistor of reservoir to be analyzed is obtained, the region that bright white is shown as in imaging is defined as described
The dried layer of reservoir to be analyzed;And/or
The bilaterally curve of reservoir to be analyzed is obtained, the region that high resistant is shown as in the bilaterally curve and is overlapped is defined as
The dried layer of the reservoir to be analyzed.
13. such as the method any one of claim 1~6,9,11,12, it is characterised in that identified in the types of fractures
In step:
When the density porosity curve is higher than the acoustic porosity curve, being split for high angle for interval development is judged
Seam;
When the acoustic porosity curve is higher than the density porosity curve, being split for low angle for interval development is judged
Seam;
When the acoustic porosity curve and the density porosity curve intersect, judge interval development to be netted
Crack.
14. method as claimed in claim 7, it is characterised in that in the types of fractures identification step:
When the density porosity curve is higher than the acoustic porosity curve, being split for high angle for interval development is judged
Seam;
When the acoustic porosity curve is higher than the density porosity curve, being split for low angle for interval development is judged
Seam;
When the acoustic porosity curve and the density porosity curve intersect, judge interval development to be netted
Crack.
15. method as claimed in claim 8, it is characterised in that in the types of fractures identification step:
When the density porosity curve is higher than the acoustic porosity curve, being split for high angle for interval development is judged
Seam;
When the acoustic porosity curve is higher than the density porosity curve, being split for low angle for interval development is judged
Seam;
When the acoustic porosity curve and the density porosity curve intersect, judge interval development to be netted
Crack.
16. method as claimed in claim 10, it is characterised in that in the types of fractures identification step:
When the density porosity curve is higher than the acoustic porosity curve, being split for high angle for interval development is judged
Seam;
When the acoustic porosity curve is higher than the density porosity curve, being split for low angle for interval development is judged
Seam;
When the acoustic porosity curve and the density porosity curve intersect, judge interval development to be netted
Crack.
17. such as the method any one of claim 1~6,9,11,12,14-16, it is characterised in that in the crack class
After type identification step, methods described also includes verification step, and the verification step includes:
Depth playback is carried out to the rock core of the reservoir to be analyzed of acquisition, judges each layer of the rock core after depth playback
Whether the types of fractures of corresponding interval of the types of fractures of section to identifying to obtain in the crack identification step matches;And/or
The types of fractures of each interval of reservoir to be analyzed is identified using imaging logging, judges what the imaging logging recognized
Whether the types of fractures of corresponding interval of the types of fractures of each interval to identifying to obtain in the crack identification step matches.
18. method as claimed in claim 7, it is characterised in that after the types of fractures identification step, methods described is also
Including verification step, the verification step includes:
Depth playback is carried out to the rock core of the reservoir to be analyzed of acquisition, judges each layer of the rock core after depth playback
Whether the types of fractures of corresponding interval of the types of fractures of section to identifying to obtain in the crack identification step matches;And/or
The types of fractures of each interval of reservoir to be analyzed is identified using imaging logging, judges what the imaging logging recognized
Whether the types of fractures of corresponding interval of the types of fractures of each interval to identifying to obtain in the crack identification step matches.
19. method as claimed in claim 8, it is characterised in that after the types of fractures identification step, methods described is also
Including verification step, the verification step includes:
Depth playback is carried out to the rock core of the reservoir to be analyzed of acquisition, judges each layer of the rock core after depth playback
Whether the types of fractures of corresponding interval of the types of fractures of section to identifying to obtain in the crack identification step matches;And/or
The types of fractures of each interval of reservoir to be analyzed is identified using imaging logging, judges what the imaging logging recognized
Whether the types of fractures of corresponding interval of the types of fractures of each interval to identifying to obtain in the crack identification step matches.
20. method as claimed in claim 10, it is characterised in that after the types of fractures identification step, methods described
Also include verification step, the verification step includes:
Depth playback is carried out to the rock core of the reservoir to be analyzed of acquisition, judges each layer of the rock core after depth playback
Whether the types of fractures of corresponding interval of the types of fractures of section to identifying to obtain in the crack identification step matches;And/or
The types of fractures of each interval of reservoir to be analyzed is identified using imaging logging, judges what the imaging logging recognized
Whether the types of fractures of corresponding interval of the types of fractures of each interval to identifying to obtain in the crack identification step matches.
21. method as claimed in claim 13, it is characterised in that after the types of fractures identification step, methods described
Also include verification step, the verification step includes:
Depth playback is carried out to the rock core of the reservoir to be analyzed of acquisition, judges each layer of the rock core after depth playback
Whether the types of fractures of corresponding interval of the types of fractures of section to identifying to obtain in the crack identification step matches;And/or
The types of fractures of each interval of reservoir to be analyzed is identified using imaging logging, judges what the imaging logging recognized
Whether the types of fractures of corresponding interval of the types of fractures of each interval to identifying to obtain in the crack identification step matches.
22. a kind of imaging logging Reservoir Fracture recognition methods, it is characterised in that the described method comprises the following steps:
Imaging logging identification step, the crack in reservoir to be identified is identified using imaging logging, obtained described to be identified
The type in crack in reservoir;
Conventional logging identification step, using the method as any one of claim 1~21 in the reservoir to be identified
Crack be identified, obtain the type in crack in the reservoir to be identified;
Correct step, according to the type in crack in the reservoir to be identified obtained in the conventional logging identification step with it is described
The type in the crack obtained in imaging logging identification step, judge what is explained in the imaging logging identification step with the presence or absence of leakage
Crack, to correct the recognition result of the imaging logging identification step.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410130972.2A CN104977617B (en) | 2014-04-02 | 2014-04-02 | Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410130972.2A CN104977617B (en) | 2014-04-02 | 2014-04-02 | Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104977617A CN104977617A (en) | 2015-10-14 |
CN104977617B true CN104977617B (en) | 2018-03-13 |
Family
ID=54274292
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410130972.2A Active CN104977617B (en) | 2014-04-02 | 2014-04-02 | Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104977617B (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105422089B (en) * | 2015-11-30 | 2018-07-27 | 长江大学 | A method of it generating porosity using well week ultrasonoscopy and composes |
CN107301255B (en) * | 2016-04-15 | 2021-01-05 | 中国石油化工股份有限公司 | Method for calculating reservoir fracture porosity |
CN108802192B (en) * | 2017-05-03 | 2021-05-14 | 中国石油化工股份有限公司 | Calcium debris sandstone reservoir pore type identification method |
CN109100793B (en) * | 2017-06-20 | 2020-06-23 | 中国石油化工股份有限公司 | Method for quantitatively analyzing influence of fracture factors on reservoir |
CN110895704B (en) * | 2018-09-12 | 2023-04-18 | 中国石油天然气股份有限公司 | Microorganism dune complex reservoir type identification method and device and storage medium |
CN111239849A (en) * | 2020-02-10 | 2020-06-05 | 东北石油大学 | Reservoir fracture prediction method and prediction system based on stress release |
CN113818870B (en) * | 2020-06-02 | 2024-03-26 | 中国石油天然气股份有限公司 | Reservoir fluid property identification method and system |
CN111764895B (en) * | 2020-08-12 | 2023-01-31 | 重庆地质矿产研究院 | Logging evaluation method suitable for shale gas reservoir geological model |
CN112682034B (en) * | 2020-12-04 | 2022-12-09 | 中国地质大学(北京) | Method and device for crack identification and dip angle characterization based on tight sandstone reservoir |
CN112946743B (en) * | 2021-03-19 | 2024-02-27 | 成都捷科思石油天然气技术发展有限公司 | Method for distinguishing reservoir types |
CN113484215B (en) * | 2021-07-06 | 2022-09-23 | 西安石油大学 | Sandstone reservoir microscopic pore analysis system based on geological cause |
CN116010789B (en) * | 2023-03-21 | 2023-06-23 | 中国石油天然气股份有限公司 | Carbonate reservoir type identification method, device, equipment and application |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6473696B1 (en) * | 2001-03-13 | 2002-10-29 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
CN102454398B (en) * | 2010-10-28 | 2015-06-17 | 中国石油化工股份有限公司 | Gas layer and water layer recognition method for low-porosity and low-permeability reservoir |
CN102749651B (en) * | 2012-07-12 | 2015-07-08 | 恒泰艾普石油天然气技术服务股份有限公司 | Calibration method for fracture logging of reservoir |
CN103592690B (en) * | 2013-10-24 | 2016-09-14 | 长江大学 | Method based on electric imaging logging porosity spectrum information automatic Reservoir Recognition crack |
-
2014
- 2014-04-02 CN CN201410130972.2A patent/CN104977617B/en active Active
Non-Patent Citations (2)
Title |
---|
储层裂缝识别与预测;王时林 等;《断块油气田》;20090930;第16卷(第5期);第31-33页 * |
镇泾地区延长组8 段致密储层裂缝识别与预测;张世懋 等;《测井技术》;20110228;第35卷(第1期);第36-40页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104977617A (en) | 2015-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104977617B (en) | Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods | |
Ameen et al. | Diverse fracture properties and their impact on performance in conventional and tight-gas reservoirs, Saudi Arabia: The Unayzah, South Haradh case study | |
Kundert et al. | Proper evaluation of shale gas reservoirs leads to a more effective hydraulic-fracture stimulation | |
CN104818982B (en) | Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir | |
CN106951660A (en) | Sea facies clastic rock horizontal well reservoir logging interpretation method and device | |
Dance et al. | Illuminating the geology: Post-injection reservoir characterisation of the CO2CRC Otway site | |
Belaidi et al. | The Lancaster Field: progress in opening the UK's fractured basement play | |
Ejofodomi et al. | Integrating all available data to improve production in the Marcellus Shale | |
CN109085663A (en) | A kind of tight sandstone reservoir stratification seam recognition methods | |
KR101591430B1 (en) | Method of measuring subsurface structure and method of drilling shale gas using the same | |
CN106338765B (en) | Method and system for determining the composite elastic parameter oil-gas recognition factor | |
Neuhaus et al. | Integrated microseismic monitoring for field optimization in the Marcellus Shale-A case study | |
CN105986813A (en) | Quasi tight reservoir rapid evaluation method and quasi tight reservoir multi-index evaluation method | |
Ouenes et al. | Integrated characterization and simulation of the fractured Tensleep Reservoir at Teapot Dome for CO2 injection design | |
Ejofodomi et al. | Development of an optimized completion strategy in the Vaca Muerta Shale with an anisotropic geomechanical model | |
CN107229076A (en) | A kind of method that temperature-responsive signature analysis is carried out based on well-log information | |
CN109752761A (en) | Method and device for evaluating reservoir characteristics of dissolved hydrocarbon reservoir | |
CN105785445B (en) | Crack pattern recognition method and device | |
CN106199690A (en) | The Forecasting Methodology in mud shale crack | |
Burns et al. | The challenges of shale gas exploration and appraisal in Europe and North Africa | |
Rauch-Davies et al. | Diffraction imaging applied to pre-existing 3D seismic data to map fracture corridors in an unconventional play | |
Altman et al. | Coupling production data with the DFN to unravel natural fracture networks in a tight naturally fractured jurassic Reservoir, Kuwait | |
Bello et al. | Multiuse of seismic and attribute mapping for field appraisal and development in the Eagle Ford Shale: Mapping TOC, porosity and seal integrity | |
Pavlyukov et al. | TIV-Anisotropy in Geomechanical Modeling for Planning of Hydraulic Fracturing at the Kharampurskoye Field | |
CN107766690A (en) | Shale gas reservoir geologic assessment modeling for single well method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |