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 PDF

Info

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
Application number
CN201410130972.2A
Other languages
Chinese (zh)
Other versions
CN104977617A (en
Inventor
魏修平
李�浩
王丹丹
冯琼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201410130972.2A priority Critical patent/CN104977617B/en
Publication of CN104977617A publication Critical patent/CN104977617A/en
Application granted granted Critical
Publication of CN104977617B publication Critical patent/CN104977617B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of 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

Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods
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>&amp;phi;</mi> <mrow> <mi>A</mi> <mi>C</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>-</mo> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;Delta;t</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> </mfrac> <mo>&amp;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>&amp;times;</mo> <mfrac> <mrow> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>s</mi> <mi>h</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;Delta;t</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>&amp;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>&amp;phi;</mi> <mrow> <mi>A</mi> <mi>C</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>-</mo> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;Delta;t</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> </mfrac> <mo>&amp;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>&amp;times;</mo> <mfrac> <mrow> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>s</mi> <mi>h</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;Delta;t</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>&amp;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>&amp;phi;</mi> <mrow> <mi>D</mi> <mi>E</mi> <mi>N</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;rho;</mi> <mo>-</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;rho;</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> </mfrac> <mo>&amp;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>&amp;times;</mo> <mfrac> <mrow> <msub> <mi>&amp;rho;</mi> <mrow> <mi>s</mi> <mi>h</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>m</mi> <mi>a</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;rho;</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>&amp;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.
CN201410130972.2A 2014-04-02 2014-04-02 Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods Active CN104977617B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
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