CN103454685B - Study of The Impedence Inversion Restrained By Well Log is utilized to predict the method and apparatus of sand thickness - Google Patents
Study of The Impedence Inversion Restrained By Well Log is utilized to predict the method and apparatus of sand thickness Download PDFInfo
- Publication number
- CN103454685B CN103454685B CN201310346850.2A CN201310346850A CN103454685B CN 103454685 B CN103454685 B CN 103454685B CN 201310346850 A CN201310346850 A CN 201310346850A CN 103454685 B CN103454685 B CN 103454685B
- Authority
- CN
- China
- Prior art keywords
- sand
- wave impedance
- thickness
- group
- total wave
- 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 provides a kind of method and apparatus utilizing Study of The Impedence Inversion Restrained By Well Log to predict sand thickness, the method comprises the following steps: obtain corresponding geological data and geologic logging data carry out seismic prospecting drilling well sampling in target exploration district after; With described geologic logging data for constraint condition, according to described geological data, Study of The Impedence Inversion Restrained By Well Log is carried out to described target exploration district, obtain the Study of The Impedence Inversion Restrained By Well Log model in described target exploration district; On the inverting section of described Study of The Impedence Inversion Restrained By Well Log model, each the sand group for described target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group; The sand thickness of corresponding sand group is obtained according to described relation curve.Invention increases for complex geologic conditions, reservoir thickness differs, physical property changes the sand-body distribution on stratum greatly and the reliability of thickness prediction.
Description
Technical field
The present invention relates to oil-gas exploration technical field, particularly relate to a kind of method and apparatus utilizing Study of The Impedence Inversion Restrained By Well Log to predict sand thickness.
Background technology
In geologic prospecting, sand thickness technology conventional at present comprises seismic attribute and Study of The Impedence Inversion Restrained By Well Log technology etc.
Seismic attribute abstraction refers to the concrete measurement of the geometric shape of the relevant seismic event that geological data is derived through mathematic(al) manipulation, kinematics character, dynamic characteristic and statistics feature.It is a complicated process that seismic event is propagated in the earth formation, is a kind of concentrated expression to subsurface formations feature, and the feature of seismic signal is directly caused by rock stratum physical property and variation thereof.The spatial variations of subsurface formations character, must cause the change of seismic reflection wave characteristic, and then affects the change of seismic properties.When subsurface formations oily, seismic response can change, and same seismic properties also changes thereupon.Seismic event often changes than physical properties of rock the attribute abnormal caused by the seismic properties ANOMALOUS VARIATIONS produced during hydrocarbon-bearing formation to be given prominence to.Therefore, on the one hand, seismic properties carries subsurface formations information, also certainly exists the inner link of certain form on the other hand, the theoretical foundation of Discussion of Earthquake Attribute Technology predicting oil/gas that Here it is between seismic properties and the oil-gas possibility of reservoir.
In general in seismic properties, time generic performance provides the information of structure; Amplitude generic attribute can be used for identifying following geological information: the reefs, unconformity surface, SEQUENCE STRATIGRAPHIC change, tuning effect etc. of gas-bearing sandstone, lithology, pore formation, river course and delta sandbody, some type, wherein instantaneous amplitude may give prominence to change and the details of storage and collection performance obviously, and in window, average amplitude and net amplitude can reflect regularity and the overall picture of reservoir better; Contained by the change of frequency generic attribute general formation structure and stratum, fluid reflection is more responsive, and geologic geophysical meaning is also relatively clearer and more definite, can disclose Fractured Zone and gassiness uptake zone.When extracting multiple seismic properties, amplitude generic attribute and the general correlation of accumulation property, frequency generic attribute is then because absorption and the attenuation of hole and internal flow thereof is stronger and show negative correlation than nonreservoir.This two generic attribute generally also has good correlativity and stability on reflection Reservoir and oil-gas possibility ability.
Study of The Impedence Inversion Restrained By Well Log normally utilizes actual three dimensional seismic data, with geological drilling and well logging information for constraint condition, to the process that tectonic structure and reservoir properties solve.
The basic ideas of Study of The Impedence Inversion Restrained By Well Log first set up an impedance initial value model, then model carries out Seismic forward thus, try to achieve theogram, by theogram compared with real seismic record, according to comparative result, the speed of amendment underground surge impedance model, density, one-tenth-value thickness 1/10 and wavelet, just drilling again and asking for theogram, after comparing with real seismic record, continue amendment surge impedance model, so repeatedly repeatedly, thus constantly pass through iterative modifications, until residual error little to a certain extent (i.e. theogram and real seismic record closest), termination of iterations, obtain final surge impedance model (i.e. Study of The Impedence Inversion Restrained By Well Log model).Model inversion can be formulated as:
Wherein, M is the model upgraded; M0 is initial model; G is sensitivity matrix, or claims Jacobi's operator, the matrix that it is made up of a series of partial derivative; Cn is noise covariance matrix; Cm is model covariance matrix; S is geological data; D is the geological data calculated; S-D claims residual deviation or residual error.Common high-resolution seismic exploration section can not differentiate thin reservoir, and the high-frequency information that enriches with well-log information of Study of The Impedence Inversion Restrained By Well Log technology and complete low-frequency component supplement the band-limited deficiency of earthquake, comprehensive geology, well logging information, as constraint condition, obtain high-precision wave impedance data.
But along with the development of oil exploration technology, Exploration Domain constantly expands.For complex geologic conditions, reservoir thickness differs, and physical property changes stratum greatly, get its reservoir characteristic and the regularity of distribution thereof clear, only uses single conventional method to have difficulties.
Summary of the invention
The object of the embodiment of the present invention is to provide a kind of method and apparatus utilizing Study of The Impedence Inversion Restrained By Well Log to predict sand thickness, and to improve for complex geologic conditions, reservoir thickness differs, physical property changes the sand-body distribution on stratum greatly and the reliability of thickness prediction.
For achieving the above object, the invention provides a kind of method utilizing Study of The Impedence Inversion Restrained By Well Log to predict sand thickness, comprise the following steps:
Corresponding geological data and geologic logging data are obtained carry out seismic prospecting drilling well sampling in target exploration district after;
With described geologic logging data for constraint condition, according to described geological data, Study of The Impedence Inversion Restrained By Well Log is carried out to described target exploration district, obtain the Study of The Impedence Inversion Restrained By Well Log model in described target exploration district;
On the inverting section of described Study of The Impedence Inversion Restrained By Well Log model, each the sand group for described target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group;
The sand thickness of corresponding sand group is obtained according to described relation curve.
For achieving the above object, the present invention also provides a kind of device utilizing Study of The Impedence Inversion Restrained By Well Log to predict sand thickness, comprising:
Data acquisition module, obtains corresponding geological data and geologic logging data for carry out seismic prospecting drilling well sampling in target exploration district after;
Wave impedance inversion module, for described geologic logging data for constraint condition, according to described geological data, Study of The Impedence Inversion Restrained By Well Log is carried out to described target exploration district, obtains the Study of The Impedence Inversion Restrained By Well Log model in described target exploration district;
Sand thickness prediction module, for the inverting section at described Study of The Impedence Inversion Restrained By Well Log model, each sand group for described target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group, obtains the sand thickness of corresponding sand group according to described relation curve.
The present invention, on the basis that Study of The Impedence Inversion Restrained By Well Log is predicted, adopts sampled point statistical method prediction sand thickness.Namely, on the basis of Study of The Impedence Inversion Restrained By Well Log prediction, each the sand group for target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group; Then the sand thickness of corresponding sand group is obtained according to relation curve.Because the lithological change of Study of The Impedence Inversion Restrained By Well Log model reflection is abundant in content, longitudinal layering is good, lateral resolution is higher, and sampling number according to statistics rule derive sand body distribution and thickness thereof according to some actual samples point data research, thus the reliability of the identifiability of reservoir and inversion prediction result is significantly improved, thus it can improve for complex geologic conditions, reservoir thickness differs, physical property changes the sand-body distribution on stratum greatly and the reliability of thickness prediction.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:
Fig. 1 is that the Study of The Impedence Inversion Restrained By Well Log that utilizes of the embodiment of the present invention predicts the process flow diagram of the method for sand thickness;
Fig. 2 is that the Study of The Impedence Inversion Restrained By Well Log that utilizes of the embodiment of the present invention predicts the structural representation of the device of sand thickness;
Fig. 3 is the 3rd of A sloped region in the embodiment of the present invention is palaeogeologic map before deposition;
Fig. 4 is 1-field of razor clam sea, the field of razor clam sea 10 well Permian System comparison diagrams of A sloped region in the embodiment of the present invention;
Fig. 5 a is the Permian system sand body wave impedance probability statistics figure of A sloped region in the embodiment of the present invention;
Fig. 5 b is the Permian system wave impedance of mudstone probability statistics figure of A sloped region in the embodiment of the present invention;
Fig. 6 is the Permian system Study of The Impedence Inversion Restrained By Well Log sectional view of A sloped region in the embodiment of the present invention;
Fig. 7 is total wave impedance value planimetric map of the Permian system one sand group of A sloped region in the embodiment of the present invention;
Fig. 8 is the total wave impedance value of sand body of the Permian system one sand group of A sloped region in the embodiment of the present invention and the graph of relation of sand thickness;
Fig. 9 is the sand thickness figure of the Permian system one sand group of A sloped region in the embodiment of the present invention.
Embodiment
For making the object of the embodiment of the present invention, technical scheme and advantage clearly understand, below in conjunction with embodiment and accompanying drawing, the embodiment of the present invention is described in further details.At this, schematic description and description of the present invention is for explaining the present invention, but not as a limitation of the invention.
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is described in further detail.
Refer to shown in Fig. 1, Fig. 1 is that the Study of The Impedence Inversion Restrained By Well Log that utilizes of the embodiment of the present invention predicts the method for sand thickness, and the method comprises the following steps:
Step 101, the geological data carrying out the rear acquisition correspondence of seismic prospecting drilling well sampling in target exploration district and geologic logging data.
Step 102, with geologic logging data for constraint condition, according to geological data, Study of The Impedence Inversion Restrained By Well Log is carried out to target exploration district, obtains the Study of The Impedence Inversion Restrained By Well Log model in target exploration district.For certain A sloped region, its Study of The Impedence Inversion Restrained By Well Log process is as follows:
1), fine structural interpretation.
At work area, A sloped region area 340km2, be all the full overlay area of 3-D seismics, three-dimensional bin is 25 × 25m, and seismic data frequency is mainly distributed between 5 ~ 80Hz, and wherein zone of interest dominant frequency is 30Hz.
According to the result of fine reservoir prediction, carried out tracing of horizons comparative interpretation to 5 seismic reflection layers of local area, being respectively the 3rd is (E) end circle, the secondary group (Mz) end circle, the Carboniferous system (C) end circle, the Permian system (P) end circle and the Ordovician system (O) end circle.According to geology log data, the secondary group (Mz) stratum is based on Jurassic systerm (J), the Triassic system (T), and Cretaceous System (K) degrades seriously, and the whole district lacks.
Through the Fine structural interpretation to A sloped region, completing this district the 3rd is palaeogeologic map (as shown in Figure 3) before deposition.As can be seen from Figure 3 the secondary group, A sloped region, the Permian system, the Carboniferous system, Ordovician strata are upper super to high-order bit, slope (southeastern direction) successively, are successively degraded.Therefore, Permian System is in south, work area disappearance.
2), meticulous individual well Reservoir Analysis.
The prospect pit drilling the Permian system in study area is less, only field of razor clam Hai1Jing and field of razor clam Hai10Jing.Geologic logging data according to ancient 1 well of field of razor clam Hai1Jing, field of razor clam Hai10Jing and sea, adjacent area are known, and local area Permian system lithology is light gray, aubergine, white sand body, folder celadon sand body and mud stone, press from both sides 3 ~ 4 cover coal seams in middle and lower part brownish red, grey mud stone.In conjunction with seismic data analysis, geologic logging data think that south, work area Permian System seriously degraded, lost strata.As shown in Figure 4, by these two mouthfuls of wells of comparative analysis, Permian system sand body is divided into a sand group, two sand groups, three sand groups from top to bottom.Viewed from the 1 well data of field of razor clam sea, sand body gross thickness and thickness in monolayer are thinning all gradually from top to bottom for three sand groups, grow four cover thin seams in three sand groups.On this basis, the sand body in the 10 well Permian Systems of field of razor clam sea is all divided into three sand groups, and this well point place lacks a sand group and two sand groups.
3), Study of The Impedence Inversion Restrained By Well Log prediction sand-body distribution.
There are geologic logging data according to oneself, added up the sand of Permian system field of razor clam Hai1Jing, field of razor clam Hai10Jing and village Hai4x1Jing (not drilling the Permian system), wave impedance of mudstone value, respectively as shown in figure 5 a and 5b.
Permian system sand body wave impedance scope is at 9000 ~ 10500m/s.g/cm
3, wave impedance of mudstone scope is at 7500 ~ 9500m/s.g/cm
3, sand, wave impedance of mudstone difference are obvious, use Study of The Impedence Inversion Restrained By Well Log can predict Sandbody Reservoirs preferably.Utilize the geologic logging data of field of razor clam Hai1Jing, field of razor clam Hai10Jing and village Hai4x1Jing as constraint condition, in conjunction with geological data, Study of The Impedence Inversion Restrained By Well Log has been done to Permian System in work area.Fig. 6 was the Permian system Study of The Impedence Inversion Restrained By Well Log section of field of razor clam Hai1Jing, according to the position, layer position of three sand groups that drilling data divides, field of razor clam Hai1Jing and field of razor clam sea 10 well well point place accurate calibrations three sand groups (field of razor clam Hai10Jing only has three sand groups) on Study of The Impedence Inversion Restrained By Well Log model, as can be seen from Figure 6, the lithological change of this inverting section reflection is abundant in content, longitudinal layering is good, and lateral resolution is higher.Final by field of razor clam Hai1Jing and field of razor clam Hai10Jing, successively follow the trail of the top bottom boundary of three sand groups, three sand groups discrimination power on inverting section is higher, thus improves the accuracy of the prediction Permian system three sand group distribution ranges.
In general, on Study of The Impedence Inversion Restrained By Well Log model, large few content substantially reflecting sand body of wave impedance value in uniform thickness situation, wave impedance value is higher, and sand body content is more, and content of mudstone is less.According to this conclusion, Study of The Impedence Inversion Restrained By Well Log model is extracted its total wave impedance value to reflect the overall distribution of sand body respectively to three sand groups, Fig. 7 is the Permian system one sand group total wave impedance value planimetric map, as can be seen from Figure 7, the northern total wave impedance value in work area is comparatively large, reduces gradually southwards.Illustrate that the northern sand body in work area is thicker, be thinned to sand body pinching gradually southwards.
Step 103, on the inverting section of Study of The Impedence Inversion Restrained By Well Log model, each the sand group for target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group.For above-mentioned A sloped region, its detailed process is as follows:
1), the inverting section of Study of The Impedence Inversion Restrained By Well Log model reads the zone thickness of a certain sample point one sand group, then to set wave impedance value (such as 10000m/s.g/cm
3) read the sand thickness of sample point one sand group according to the sand added up in geology log data, wave impedance of mudstone value for boundary.In addition, sandy ground can also be calculated accordingly than the rationality checking setting wave impedance value value in order to the later stage possibly.
2) on total wave impedance value planimetric map of a sand group of Study of The Impedence Inversion Restrained By Well Log model, read total wave impedance value of sample point one sand group.
3), by total wave impedance value planimetric map of a sand group and the comparative analysis of geologic logging data, the sand body of sample point one sand group and total wave impedance boundary value (such as 100000m/s.g/cm of mud stone is determined
3), and with sand body pinching place of total wave impedance boundary value for boundary's delimitation sample point one sand group.
4), deduct by total wave impedance value of sample point one sand group the total wave impedance value of sand body that namely total wave impedance boundary value obtains sample point one sand group.
5) sand thickness and the total wave impedance value of sand body that above step obtains other sample point one sand groups in 20 sampled points respectively, is repeated.Obtain a Permian system one sand group sampled point data statistic as shown in table 1 below.
Table 1 Permian system one sand group sampled point data statistic
6), by the sand thickness of a sand group that reads at 20 sampled points and sand body total wave impedance value the total wave impedance value of sand body of a sand group and the relation curve of its sand thickness is simulated, as shown in Figure 8.During for ensureing that sand thickness is zero, the total wave impedance value of sand body is also zero, and the relation curve simulated must cross initial point.
7), repeat above step and obtain the total wave impedance value of sand body of other sand groups in target exploration district and the relation curve of its sand thickness respectively.
Step 104, the sand thickness of corresponding sand group can be obtained according to relation curve.As shown in Figure 9, Permian system sand body affects by thing source, northeast and comparatively grows, and by the region comparatively enrichment between village Hai4x1Jing, field of razor clam Hai5Jing, field of razor clam Hai2Jing and sea, the village 12 wells of sand body inside East and West direction faults control, maximum gauge is at about 180m, all in all northern sand body is thicker, thinning to south gradually.One sand group area is minimum, and maximum gauge is at about 50m, and northern thick south is thin, by southern faults control around field of razor clam Hai1Jing, forms a little sand body enrichment region; Two sand group sand bodies have two parts in work area, and area major part is comparatively smooth, and variation in thickness is little, and thickness is about about 45m, and area smaller portions variation in thickness is comparatively large, and thinning gradually by north orientation south, maximum gauge is about 55m; Three sand group areas are maximum, and average thickness is maximum, has two parts equally in work area, and area major part maximum gauge is at about 70m, and area smaller portions maximum gauge is at about 95m, and sand body variation tendency and a sand group, two sand groups are identical.
In embodiments of the present invention, after obtaining corresponding geological data and geologic logging data, seismic properties (such as amplitude generic attribute and frequency generic attribute etc.) that can also be corresponding according to seismic data acquisition, and the Sand-body Prediction distribution obtaining described target exploration district according to seismic properties.Due to amplitude attribute, to comprise the sand-body distribution situation that RMS amplitude attribute, average energy attribute and gross energy attribute show with it close, and arc length attribute difference in frequency attribute is comparatively far away, when therefore predicting the distribution of the Permian system three sand groups can amplitude generic attribute as a reference.And to obtain the sand thickness of corresponding sand group according to relation curve in step 104 after, for making result of calculation more meet geologic rule, exceptional value and the anomaly trend of the sand thickness of corresponding sand group can be revised according to Sand-body Prediction distribution.
The embodiment of the present invention utilize Study of The Impedence Inversion Restrained By Well Log to predict the method for sand thickness is on the basis of Study of The Impedence Inversion Restrained By Well Log prediction, adopt sampled point statistical method prediction sand thickness.Namely, on the basis of Study of The Impedence Inversion Restrained By Well Log prediction, each the sand group for target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group; Then the sand thickness of corresponding sand group is obtained according to relation curve.Because the lithological change of Study of The Impedence Inversion Restrained By Well Log model reflection is abundant in content, longitudinal layering is good, lateral resolution is higher, and sampling number according to statistics rule derive sand body distribution and thickness thereof according to some actual samples point data research, thus the reliability of the identifiability of reservoir and inversion prediction result is significantly improved, thus it can improve for complex geologic conditions, reservoir thickness differs, physical property changes the sand-body distribution on stratum greatly and the reliability of thickness prediction.
Shown in composition graphs 2, the device of Study of The Impedence Inversion Restrained By Well Log prediction sand thickness that utilizes of the embodiment of the present invention comprises: data acquisition module 21, wave impedance inversion module 22, sand thickness prediction module 23, sand-body distribution prediction module 24(are optional) and sand 25 body thickness correcting module (optional).Wherein: data acquisition module 21 obtains corresponding geological data and geologic logging data for carry out seismic prospecting drilling well sampling in target exploration district after. wave impedance inversion module 22 for geologic logging data for constraint condition, according to geological data, Study of The Impedence Inversion Restrained By Well Log is carried out to target exploration district, obtain the Study of The Impedence Inversion Restrained By Well Log model in target exploration district.Sand thickness prediction module 23 is for the inverting section at Study of The Impedence Inversion Restrained By Well Log model, each sand group for target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group, obtains the sand thickness of corresponding sand group according to relation curve.Sand-body distribution prediction module 24 is for after obtaining corresponding geological data and geologic logging data, and the seismic properties corresponding according to seismic data acquisition, obtains the Sand-body Prediction distribution in target exploration district according to seismic properties.Due to amplitude attribute, to comprise the sand-body distribution situation that RMS amplitude attribute, average energy attribute and gross energy attribute show with it close, and arc length attribute difference in frequency attribute is comparatively far away, when therefore predicting the distribution of the Permian system three sand groups can amplitude generic attribute as a reference.Sand thickness correcting module 25, for after the sand thickness obtaining corresponding sand group according to relation curve, for making result of calculation more meet geologic rule, can revise exceptional value and the anomaly trend of the sand thickness of corresponding sand group according to Sand-body Prediction distribution.
As can be seen from Figure 2, sand thickness prediction module 23 specifically comprises again sampled point sand thickness acquisition submodule 231, the total wave impedance of sand group obtains submodule 232, sand group total wave impedance boundary value determination submodule 233, sand body total wave impedance value calculating sub module 234, first control submodule 235, relation curve matching submodule 236, sand thickness obtain submodule 237 and second and control submodule 238.Wherein, sampled point sand thickness obtains submodule 231 for reading the zone thickness of a certain sample point one sand group on the inverting section of Study of The Impedence Inversion Restrained By Well Log model, then to set wave impedance value (such as 10000m/s.g/cm
3) read the sand thickness of sample point one sand group according to the sand added up in geology log data, wave impedance of mudstone value for boundary.The total wave impedance of sand group obtains total wave impedance value that submodule 232 reads sample point one sand group on total wave impedance value planimetric map of the sand group at Study of The Impedence Inversion Restrained By Well Log model.Sand group total wave impedance boundary value determination submodule 233, for by total wave impedance value planimetric map of a sand group and the comparative analysis of geologic logging data, determines the sand body of sample point one sand group and total wave impedance boundary value (100000m/s.g/cm of mud stone
3), and with sand body pinching place of total wave impedance boundary value for boundary's delimitation sample point one sand group.Sand body total wave impedance value calculating sub module 234 is for deducting the total wave impedance value of sand body that namely total wave impedance boundary value obtains sample point one sand group by total wave impedance value of sample point one sand group.First controls submodule 235, and for controlling, sampled point sand thickness obtains submodule 231, the total wave impedance of sand group obtains submodule 232, sand group total wave impedance boundary value determination submodule 233 and sand body total wave impedance value calculating sub module 234 repeated work to obtain sand thickness and the total wave impedance value of sand body of other sample point one sand groups in several sampled points respectively.Relation curve matching submodule 236 simulates the total wave impedance value of sand body of a sand group and the relation curve of its sand thickness for the sand thickness of the sand group by reading at several sampled points and sand body total wave impedance value, and wherein relation curve crosses initial point.Sand thickness obtains submodule 237 for obtaining the sand thickness of corresponding sand group according to relation curve.Second controls submodule 238 obtains submodule 237 repeated work with the relation curve of the total wave impedance value of the sand body obtaining other sand groups in target exploration district respectively and its sand thickness for controlling the first control submodule 235, relation curve matching submodule 236 and sand thickness.
The Study of The Impedence Inversion Restrained By Well Log that utilizes of the embodiment of the present invention is predicted in the device of sand thickness, on the basis of Study of The Impedence Inversion Restrained By Well Log prediction, adopts sampled point statistical method prediction sand thickness.Namely, on the basis of Study of The Impedence Inversion Restrained By Well Log prediction, each the sand group for target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group; Then the sand thickness of corresponding sand group is obtained according to relation curve.Because the lithological change of Study of The Impedence Inversion Restrained By Well Log model reflection is abundant in content, longitudinal layering is good, lateral resolution is higher, and sampling number according to statistics rule derive sand body distribution and thickness thereof according to some actual samples point data research, thus the reliability of the identifiability of reservoir and inversion prediction result is significantly improved, thus it can improve for complex geologic conditions, reservoir thickness differs, physical property changes the sand-body distribution on stratum greatly and the reliability of thickness prediction.
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any correction done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (8)
1. utilize Study of The Impedence Inversion Restrained By Well Log to predict a method for sand thickness, it is characterized in that, comprise the following steps:
Corresponding geological data and geologic logging data are obtained carry out seismic prospecting drilling well sampling in target exploration district after;
With described geologic logging data for constraint condition, according to described geological data, Study of The Impedence Inversion Restrained By Well Log is carried out to described target exploration district, obtain the Study of The Impedence Inversion Restrained By Well Log model in described target exploration district;
On the inverting section of described Study of The Impedence Inversion Restrained By Well Log model, each the sand group for described target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group;
The sand thickness of corresponding sand group is obtained according to described relation curve; Wherein:
Described on the inverting section of Study of The Impedence Inversion Restrained By Well Log model, each sand group for described target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group, specifically comprises the steps:
The inverting section of described Study of The Impedence Inversion Restrained By Well Log model reads the zone thickness of a certain sample point one sand group, then reads the sand thickness of described sample point one sand group according to the sand added up in described geologic logging data, wave impedance of mudstone value to set wave impedance value as boundary;
Total wave impedance value planimetric map of a described sand group of described Study of The Impedence Inversion Restrained By Well Log model reads total wave impedance value of described sample point one sand group;
By to total wave impedance value planimetric map of a described sand group and the comparative analysis of described geologic logging data, determine the described sand body of sample point one sand group and total wave impedance boundary value of mud stone, and with sand body pinching place of described total wave impedance boundary value described sample point one sand group for boundary delimit;
The total wave impedance value of sand body that namely described total wave impedance boundary value obtains described sample point one sand group is deducted by total wave impedance value of described sample point one sand group;
Repeat sand thickness and the total wave impedance value of sand body that above step obtains other sample point one sand groups in several sampled points described respectively;
Simulate the total wave impedance value of sand body of a described sand group and the relation curve of its sand thickness by the sand thickness of a sand group that reads at described several sampled points and sand body total wave impedance value, wherein said relation curve crosses initial point;
Repeat above step and obtain the total wave impedance value of sand body of other sand groups in described target exploration district and the relation curve of its sand thickness respectively.
2. method according to claim 1, is characterized in that, described obtain corresponding geological data and geologic logging data after, also comprise:
The seismic properties corresponding according to described seismic data acquisition, obtains the Sand-body Prediction distribution in described target exploration district according to described seismic properties;
Described obtain the sand thickness of corresponding sand group according to relation curve after, revise exceptional value and the anomaly trend of the sand thickness of described corresponding sand group according to the distribution of described Sand-body Prediction.
3. method according to claim 2, is characterized in that, described seismic properties comprises amplitude generic attribute and frequency generic attribute; Described amplitude generic attribute comprises RMS amplitude attribute, average energy attribute and gross energy attribute; Described frequency generic attribute comprises arc length attribute.
4. method according to claim 1, is characterized in that, described setting wave impedance value comprises 10000m/s.g/cm
3, described total wave impedance boundary value comprises 100000m/s.g/cm
3.
5. utilize Study of The Impedence Inversion Restrained By Well Log to predict a device for sand thickness, it is characterized in that, comprising:
Data acquisition module, obtains corresponding geological data and geologic logging data for carry out seismic prospecting drilling well sampling in target exploration district after;
Wave impedance inversion module, for described geologic logging data for constraint condition, according to described geological data, Study of The Impedence Inversion Restrained By Well Log is carried out to described target exploration district, obtains the Study of The Impedence Inversion Restrained By Well Log model in described target exploration district;
Sand thickness prediction module, for the inverting section at described Study of The Impedence Inversion Restrained By Well Log model, each sand group for described target exploration district evenly chooses several sampled points respectively with the relation curve of the total wave impedance value of sand body and its sand thickness that simulate corresponding sand group, obtains the sand thickness of corresponding sand group according to described relation curve; Wherein:
Described sand thickness prediction module, specifically comprises:
Sampled point sand thickness obtains submodule, on the inverting section at described Study of The Impedence Inversion Restrained By Well Log model, read the zone thickness of a certain sample point one sand group, then read the sand thickness of described sample point one sand group according to the sand added up in described geologic logging data, wave impedance of mudstone value to set wave impedance value as boundary;
The total wave impedance of sand group obtains submodule, for reading total wave impedance value of described sample point one sand group on total wave impedance value planimetric map of a described sand group of described Study of The Impedence Inversion Restrained By Well Log model;
Sand group total wave impedance boundary value determination submodule, for passing through total wave impedance value planimetric map of a described sand group and the comparative analysis of described geologic logging data, determine the described sand body of sample point one sand group and total wave impedance boundary value of mud stone, and with sand body pinching place of described total wave impedance boundary value described sample point one sand group for boundary delimit;
Sand body total wave impedance value calculating sub module, for deducting by total wave impedance value of described sample point one sand group the total wave impedance value of sand body that namely described total wave impedance boundary value obtains described sample point one sand group;
First controls submodule, and for controlling, described sampled point sand thickness obtains submodule, the total wave impedance of described sand group obtains submodule, described sand group total wave impedance boundary value determination submodule and described sand body total wave impedance value calculating sub module repeated work to obtain sand thickness and the total wave impedance value of sand body of other sample point one sand groups in several sampled points described respectively;
Relation curve matching submodule, simulate the total wave impedance value of sand body of a described sand group and the relation curve of its sand thickness for the sand thickness of the sand group by reading at described several sampled points and sand body total wave impedance value, wherein said relation curve crosses initial point;
Sand thickness obtains submodule, for obtaining the sand thickness of corresponding sand group according to described relation curve;
Second controls submodule, controls submodule, described relation curve matching submodule and described sand thickness acquisition submodule repeated work with the relation curve of the total wave impedance value of sand body and its sand thickness that obtain other sand groups in described target exploration district respectively for controlling described first.
6. device according to claim 5, is characterized in that, also comprises:
Sand-body distribution prediction module, for after the geological data and geologic logging data of described acquisition correspondence, the seismic properties corresponding according to described seismic data acquisition, obtains the Sand-body Prediction distribution in described target exploration district according to described seismic properties;
Sand thickness correcting module, for described obtain the sand thickness of corresponding sand group according to relation curve after, revise exceptional value and the anomaly trend of the sand thickness of described corresponding sand group according to the distribution of described Sand-body Prediction.
7. device according to claim 6, is characterized in that, described seismic properties comprises amplitude generic attribute and frequency generic attribute; Described amplitude generic attribute comprises RMS amplitude attribute, average energy attribute and gross energy attribute; Described frequency generic attribute comprises arc length attribute.
8. device according to claim 5, is characterized in that, described setting wave impedance value comprises 10000m/s.g/cm
3, described total wave impedance boundary value comprises 100000m/s.g/cm
3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310346850.2A CN103454685B (en) | 2013-08-09 | 2013-08-09 | Study of The Impedence Inversion Restrained By Well Log is utilized to predict the method and apparatus of sand thickness |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310346850.2A CN103454685B (en) | 2013-08-09 | 2013-08-09 | Study of The Impedence Inversion Restrained By Well Log is utilized to predict the method and apparatus of sand thickness |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103454685A CN103454685A (en) | 2013-12-18 |
| CN103454685B true CN103454685B (en) | 2016-04-06 |
Family
ID=49737253
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310346850.2A Active CN103454685B (en) | 2013-08-09 | 2013-08-09 | Study of The Impedence Inversion Restrained By Well Log is utilized to predict the method and apparatus of sand thickness |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103454685B (en) |
Families Citing this family (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103954995B (en) * | 2014-04-22 | 2016-12-07 | 核工业北京地质研究院 | A kind of sand body recognition methods in exploration of sandstone type uranium deposits |
| CN104090303B (en) * | 2014-07-10 | 2016-11-02 | 中国海洋石油总公司 | A kind of seismic inversion method and device |
| CN104133250B (en) * | 2014-07-23 | 2017-06-13 | 中国石油集团川庆钻探工程有限公司 | Gypsum-salt bed geology card layer method |
| CN105353411B (en) * | 2015-10-27 | 2017-10-13 | 中国石油大学(华东) | Sandstone reservoir is portrayed and description method on ultra deep coal |
| CN106772615B (en) * | 2016-11-09 | 2018-11-13 | 中国石油大学(华东) | A kind of broad-band teleseismic multiple domain joint AVO inversion methods |
| CN106842317B (en) * | 2016-12-22 | 2019-02-15 | 中国石油天然气集团公司 | A kind of method and device of prediction oil sand body distribution |
| CN106772614A (en) * | 2017-02-28 | 2017-05-31 | 中国石油天然气股份有限公司 | The Forecasting Methodology of High-quality Reservoir in a kind of Gravity-flow Channel Sandbody |
| CN107121699B (en) * | 2017-02-28 | 2019-06-21 | 中国石油化工股份有限公司 | A kind of sedimentary facies identification method under earthquake phase control |
| CN107422374A (en) * | 2017-04-21 | 2017-12-01 | 中国石油大港油田勘探开发研究院 | A kind of sags rich in oil and gas slope belt fine granularing scalability method |
| CN108931813B (en) * | 2017-05-26 | 2020-04-10 | 中国石油化工股份有限公司 | Method for identifying reef cover, reef core and reef base |
| CN107272081B (en) * | 2017-07-18 | 2018-12-28 | 中国石油大学(华东) | A kind of submountain region extension of sedimentary facies belts horizon prediction method |
| CN107664776A (en) * | 2017-10-25 | 2018-02-06 | 中国石油天然气股份有限公司 | A kind of Weathering Crusts of Carbonate Rocks karst reservoir thickness earthquake prediction method and device |
| CN109752760B (en) * | 2017-11-01 | 2021-01-01 | 中国石油天然气股份有限公司 | Method and device for determining reliability of seismic attribute |
| CN107942405B (en) * | 2017-11-15 | 2019-09-24 | 中国石油化工股份有限公司 | The method for predicting thin sand-mud interbed sand body cumulative thickness |
| CN110542922B (en) * | 2019-03-05 | 2020-12-15 | 中国石油化工股份有限公司 | Thin mutual reservoir description method taking reconstructed sand-to-ground ratio as seismic description target |
| CN111101935B (en) * | 2019-12-25 | 2023-01-03 | 中海石油(中国)有限公司 | Oil shale prediction method under few-well condition |
| CN113703042B (en) * | 2020-05-21 | 2023-09-29 | 中国石油天然气股份有限公司 | Method and device for determining accumulation thickness of interlayer clastic rock sand |
| CN112037306A (en) * | 2020-07-14 | 2020-12-04 | 中国石油天然气集团有限公司 | Sand body thickness map drawing method and device, computer equipment and storage medium |
| CN111665562B (en) * | 2020-07-20 | 2022-03-01 | 西南石油大学 | High-precision seismic sequence dividing method |
| CN111983721A (en) * | 2020-08-26 | 2020-11-24 | 核工业北京地质研究院 | Sandstone uranium ore mud-sand-mud geological structure identification method and system |
| CN113406703A (en) * | 2021-05-31 | 2021-09-17 | 中国石油天然气集团有限公司 | Method and device for predicting thickness of stratum in carboniferous work area |
| CN113296150B (en) * | 2021-07-02 | 2022-06-10 | 清华大学 | High-dimensional closed-loop network seismic inversion method under logging constraint |
| CN113703054B (en) * | 2021-09-03 | 2022-02-11 | 中国科学院地质与地球物理研究所 | Sedimentary facies evolution quantitative characterization method and system based on geology-geophysical |
| CN114859418B (en) * | 2022-06-15 | 2022-12-20 | 北京阳光吉澳能源技术有限公司 | River channel sand reservoir identification method and device and electronic equipment |
| CN115903026B (en) * | 2023-01-09 | 2023-05-09 | 东北石油大学三亚海洋油气研究院 | Method, equipment and medium for analyzing composite sand body configuration |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6502037B1 (en) * | 1999-04-02 | 2002-12-31 | Conoco Inc. | Method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
| CN101980053A (en) * | 2010-10-18 | 2011-02-23 | 中国石油化工股份有限公司 | Complicated reef flat reservoir predicting method |
| CN102109613A (en) * | 2009-12-23 | 2011-06-29 | 中国石油天然气股份有限公司 | Method for defining effective thickness of target reservoir bed under complex geological conditions |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8379483B2 (en) * | 2009-06-04 | 2013-02-19 | Baker Hughes Incorporated | Method and apparatus for determining radial shear velocity variation from dipole acoustic logging |
-
2013
- 2013-08-09 CN CN201310346850.2A patent/CN103454685B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6502037B1 (en) * | 1999-04-02 | 2002-12-31 | Conoco Inc. | Method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
| CN102109613A (en) * | 2009-12-23 | 2011-06-29 | 中国石油天然气股份有限公司 | Method for defining effective thickness of target reservoir bed under complex geological conditions |
| CN101980053A (en) * | 2010-10-18 | 2011-02-23 | 中国石油化工股份有限公司 | Complicated reef flat reservoir predicting method |
Non-Patent Citations (4)
| Title |
|---|
| 基于地震相分析的砂体储层厚度描述;蒋炼 等;《断块油气田》;20110531;第18卷(第3期);第274页第2 栏第2段第1-2行 * |
| 多井约束反演在煤层厚度定量预测中的应用;师素珍 等;《采矿与安全工程学报》;20110630;第28卷(第2期);第329页第1栏第3-5行 * |
| 曲线重构反演在储层横向预测中的应用;熊冉 等;《西南石油大学学报（自然科学版）》;20120229;第34卷(第1期);第86页第1栏倒数第1段，第2栏第1-4行 * |
| 陈家庄南区馆陶组下段砂体识别方法研究;赵长勋;《胜利油田职工大学学报》;20061031;第20卷(第5期);第35页第1栏第4段 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103454685A (en) | 2013-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103454685B (en) | Study of The Impedence Inversion Restrained By Well Log is utilized to predict the method and apparatus of sand thickness | |
| CN104977618B (en) | A kind of method evaluated shale gas reservoir and find dessert area | |
| CN104047597B (en) | Oily mud shale stratum log standardized method | |
| Sharma | Environmental and engineering geophysics | |
| CN101158724B (en) | Reservoir thickness prediction method based on dipolar wavelet | |
| CN103135135B (en) | A kind of method and apparatus carrying out hydro carbons quantitative forecast based on understressed sand model | |
| CN109425896A (en) | Dolomite oil and gas reservoir distribution forecasting method and device | |
| CN105044770A (en) | Compact glutenite gas reservoir quantificational prediction method | |
| CN104853822A (en) | Method for evaluating shale gas reservoir and searching sweet spot region | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| CN103245971A (en) | Well breaking point guided earthquake minor fault interpretation method and device | |
| CN112746837A (en) | Shale oil reservoir exploration data acquisition system and method based on distributed optical fiber sensing | |
| CN104516018A (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
| CN104570067A (en) | Phase-controlled earthquake inversion method in geophysical exploration | |
| CN214576965U (en) | Shale oil reservoir exploration data acquisition system based on distributed optical fiber sensing | |
| CN104049283A (en) | Method for identification and prediction of fore-land basin extremely-thick conglomerate body | |
| Han et al. | Seismic imaging of the metamorphism of young sediment into new crystalline crust in the actively rifting I mperial V alley, C alifornia | |
| Khalil et al. | Engineering applications of seismic refraction method: A field example at Wadi Wardan, Northeast Gulf of Suez, Sinai, Egypt | |
| Li et al. | An Integrated quantitative modeling approach for fault-related fractures in tight sandstone reservoirs | |
| CN102096098A (en) | Seismic attribute extracting method | |
| CN110187389A (en) | A kind of AVA inversion method based on thin bed reflection theory | |
| Khan et al. | Characterizing seismo-stratigraphic and structural framework of late cretaceous-recent succession of offshore Indus Pakistan | |
| CN103278852B (en) | Utilize the method for seismic data volume waveform configuration characteristic model predicting oil/gas | |
| Deng et al. | Hydrocarbon accumulation conditions and key exploration and development technologies for PL 19–3 oilfield | |
| Li et al. | Seismic Survey in Lesser Himalayan Thrust Belt, Western Nepal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |