CN103454685A - Method and device for predicating sand body thicknesses through logging constraint wave impedance inversion - Google Patents
Method and device for predicating sand body thicknesses through logging constraint wave impedance inversion Download PDFInfo
- Publication number
- CN103454685A CN103454685A CN2013103468502A CN201310346850A CN103454685A CN 103454685 A CN103454685 A CN 103454685A CN 2013103468502 A CN2013103468502 A CN 2013103468502A CN 201310346850 A CN201310346850 A CN 201310346850A CN 103454685 A CN103454685 A CN 103454685A
- 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.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method and device for predicating sand body thicknesses through logging constraint wave impedance inversion. The method includes the steps of acquiring corresponding seismic data and geological logging data after carrying out seismic exploration drilling sampling in a target exploration area, carrying out logging constraint wave impedance inversion on the target exploration area according to the seismic data with the geological logging data as constraint conditions to obtain a logging constraint wave impedance inversion model of the target exploration area, selecting a plurality of sampling points respectively and evenly for each sand group in the target exploration area on an inversion section of the logging constraint wave impedance inversion model to obtain sand body total wave impedance values of corresponding sand groups and sand body thickness relation curves of the corresponding sand groups through fitting, and obtaining the sand body thicknesses of the corresponding sand groups according to the relation curves. The method and device for predicating the sand body thicknesses through logging constraint wave impedance inversion improve reliability of sand body distribution and thickness prediction in stratums with complicated geological conditions, different reservoir stratum thicknesses and large physical property changes.
Description
Technical field
The present invention relates to the oil-gas exploration technical field, particularly relate to a kind of method and apparatus that utilizes Study of The Impedence Inversion Restrained By Well Log prediction sand thickness.
Background technology
In geologic prospecting, sand thickness technology commonly used comprises seismic attribute and Study of The Impedence Inversion Restrained By Well Log technology etc. at present.
Seismic properties is extracted the concrete measurement of geometric shape, kinematics character, dynamic characteristic and the statistics feature of the relevant seismic event that refers to that geological data is derived through mathematic(al) manipulation.It is a complicated process that seismic event is propagated in stratum, is a kind of concentrated expression to the 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 variation of seismic reflection wave characteristic, and then affect the variation of seismic properties.When the subsurface formations oily, seismic response can change, and same seismic properties also changes thereupon.It is outstanding that the seismic properties ANOMALOUS VARIATIONS that seismic event produces during by hydrocarbon-bearing formation often changes than physical properties of rock the attribute abnormal caused.Therefore, on the one hand, seismic properties carries subsurface formations information, also certainly exists on the other hand the inner link of certain form between the oil-gas possibility of seismic properties and reservoir, the theoretical foundation of Discussion of Earthquake Attribute Technology predicting oil/gas that Here it is.
In seismic properties, time generic performance provides the information of structure in general; The amplitude generic attribute can be used for identifying following geological information: the reefs of gas-bearing sandstone, lithology, pore formation, river course and delta sandbody, some type, unconformity surface, SEQUENCE STRATIGRAPHIC variation, tuning effect etc., wherein instantaneous amplitude may be given prominence to variation and the details of storage and collection performance apparent in viewly, and the interior average amplitude of window and net amplitude can reflect regularity and the overall picture of reservoir better; The frequency generic attribute is generally more responsive to the contained fluid reflection of variation and stratum of stratal configuration, and the geologic geophysical meaning is also clearer and more definite, can disclose Fractured Zone and gassiness uptake zone.When extracting multiple seismic properties, the amplitude generic attribute with preserve the general correlation of characteristic, the frequency generic attribute due to the absorption of hole and internal flow thereof with attenuation than the strong negative correlation that shows of 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, take geological drilling and well logging information as constraint condition, the process that tectonic structure and reservoir properties are solved.
The basic ideas of Study of The Impedence Inversion Restrained By Well Log are first to set up an impedance initial value model, then model carries out Seismic forward thus, try to achieve theogram, theogram is compared with real seismic record, according to comparative result, revise the speed of underground surge impedance model, density, one-tenth-value thickness 1/10 and wavelet, just drilling and asking for theogram again, with real seismic record relatively after, continue to revise surge impedance model, so repeatedly repeatedly, thereby constantly pass through iterative modifications, until residual error little to a certain extent (being that theogram and real seismic record are the most approaching), termination of iterations, obtain final surge impedance model (being the Study of The Impedence Inversion Restrained By Well Log model).Model inversion can be formulated as:
Wherein, the model of M for upgrading; M0 is initial model; G is sensitivity matrix, or claims Jacobi's operator, the matrix that it is comprised of a series of partial derivatives; Cn is noise covariance matrix; Cm is the 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 be differentiated thin reservoir, and Study of The Impedence Inversion Restrained By Well Log technology high-frequency information and the complete low-frequency component abundant with well-log information supplements the band-limited deficiency of earthquake, comprehensive geology, well logging information, as constraint condition, obtain high-precision wave impedance data.
Yet, along with the development of oil exploration technology, Exploration Domain constantly enlarges.For complex geologic conditions, the 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 purpose of the embodiment of the present invention is to provide a kind of method and apparatus that utilizes Study of The Impedence Inversion Restrained By Well Log prediction sand thickness, to improve for the sand body that complex geologic conditions, reservoir thickness differ, physical property changes stratum greatly, distributes and the reliability of thickness prediction.
For achieving the above object, the invention provides a kind of method of utilizing Study of The Impedence Inversion Restrained By Well Log prediction sand thickness, comprise the following steps:
After carrying out seismic prospecting drilling well sampling, the target exploration district obtains corresponding geological data and geologic logging data;
The described geologic logging data of take are constraint condition, according to described geological data, described target exploration district are carried out to Study of The Impedence Inversion Restrained By Well Log, 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, for each sand group in described target exploration district, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness;
Obtain the sand thickness of corresponding sand group according to described relation curve.
For achieving the above object, the present invention also provides a kind of device that utilizes Study of The Impedence Inversion Restrained By Well Log prediction sand thickness, comprising:
Data acquisition module, obtain corresponding geological data and geologic logging data for carrying out in the target exploration district after seismic prospecting drilling well sampling;
The wave impedance inversion module, be constraint condition for take described geologic logging data, according to described geological data, described target exploration district carried out to Study of The Impedence Inversion Restrained By Well Log, obtains the Study of The Impedence Inversion Restrained By Well Log model in described target exploration district;
The sand thickness prediction module, for the inverting section at described Study of The Impedence Inversion Restrained By Well Log model, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness for each sand group in described target exploration district, obtain the sand thickness of corresponding sand group according to described relation curve.
The present invention, on the basis of Study of The Impedence Inversion Restrained By Well Log prediction, adopts sampled point statistical method prediction sand thickness.On the basis of Study of The Impedence Inversion Restrained By Well Log prediction, for each sand group in target exploration district, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness; Then obtain the sand thickness of corresponding sand group according to relation curve.Because the lithological change of Study of The Impedence Inversion Restrained By Well Log model reflection is abundant in content, vertically layering is good, lateral resolution is higher, and sampling number according to statistics rule according to some actual samples point data researches derivation sand body distribution and thickness thereof, thereby the identifiability of reservoir and the reliability of inversion prediction result are significantly improved, thereby it can improve for complex geologic conditions, reservoir thickness differ, physical property changes stratum greatly sand body and distribute and the reliability of thickness prediction.
The accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms the application's a part, does not form limitation of the invention.In the accompanying drawings:
The process flow diagram of the method for utilizing Study of The Impedence Inversion Restrained By Well Log prediction sand thickness that Fig. 1 is the embodiment of the present invention;
The structural representation of the device that utilizes Study of The Impedence Inversion Restrained By Well Log prediction sand thickness that Fig. 2 is the embodiment of the present invention;
Fig. 3 be in the embodiment of the present invention A sloped region the 3rd be the deposition before palaeogeologic map;
1-field of razor clam sea, the field of razor clam sea 10 well Permian System comparison diagrams that Fig. 4 is A sloped region in the embodiment of the present invention;
The Permian system sand body wave impedance probability statistics figure that Fig. 5 a is A sloped region in the embodiment of the present invention;
The Permian system wave impedance of mudstone probability statistics figure that Fig. 5 b is A sloped region in the embodiment of the present invention;
The Permian system Study of The Impedence Inversion Restrained By Well Log sectional view that Fig. 6 is A sloped region in the embodiment of the present invention;
Total wave impedance value planimetric map of the Permian system one sand group that Fig. 7 is A sloped region in the embodiment of the present invention;
The total wave impedance value of sand body of the Permian system one sand group that Fig. 8 is A sloped region in the embodiment of the present invention and the graph of relation of sand thickness;
The sand thickness figure of the Permian system one sand group that Fig. 9 is A sloped region in the embodiment of the present invention.
Embodiment
For the purpose, technical scheme and the advantage that make the embodiment of the present invention is clearer, 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, the method for utilizing Study of The Impedence Inversion Restrained By Well Log prediction sand thickness that Fig. 1 is the embodiment of the present invention, the method comprises the following steps:
1), fine structures is explained.
At A sloped region work area area 340km2, be all the full overlay area of 3-D seismics, three-dimensional bin is 25 * 25m, and the seismic data frequency mainly is distributed between 5～80Hz, and wherein the zone of interest dominant frequency is 30Hz.
Result according to meticulous horizon calibration, 5 seismic reflection layers to local area have carried out the tracing of horizons comparative interpretation, and being respectively the 3rd is (E) circle, the end, circle, the end, the secondary group (Mz), the Carboniferous system (C) circle, the end, the Permian system (P) circle, the end and the Ordovician system (O) circle, the end.Known according to the geology log data, Jurassic systerm (J), the Triassic system (T) be take as main in stratum, the secondary group (Mz), and Cretaceous System (K) degrades seriously, whole district's disappearance.
Through the Fine structural interpretation to the A sloped region, completing Liao Gai district the 3rd is palaeogeologic map (as shown in Figure 3) before deposition.As can be seen from Figure 3 the A sloped region secondary group, 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 that drills 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, and 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.The geologic logging data think that in conjunction with seismic data analysis 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 to a sand group, two sand groups, three sand groups from top to bottom.See three sand groups equal attenuate gradually of sand body gross thickness and thickness in monolayer from top to bottom from field of razor clam sea 1 well data, in three sand groups, grow the quadruplet thin seams.On this basis, the sand body in the 10 well Permian Systems of field of razor clam sea all is divided into to three sand groups, these place, well point disappearance one sand group and two sand groups.
3), Study of The Impedence Inversion Restrained By Well Log prediction sand body distributes.
According to oneself, the geologic logging data are arranged, added up sand, the wave impedance of mudstone value of Permian system field of razor clam Hai1Jing, field of razor clam Hai10Jing and village Hai4x1Jing (not drilling the Permian system), respectively as shown in Fig. 5 a and Fig. 5 b.
Permian system sand body wave impedance scope is at 9000～10500m/s.g/cm
3, the 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 preferably Sandbody Reservoirs.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, Permian System in work area has been done to Study of The Impedence Inversion Restrained By Well Log.The Permian system Study of The Impedence Inversion Restrained By Well Log section that Fig. 6 was field of razor clam Hai1Jing, the position, layer position of three sand groups of dividing according to drilling data, on the Study of The Impedence Inversion Restrained By Well Log model field of razor clam Hai1Jing and field of razor clam sea 10 place, well well points accurate calibrations three sand groups (field of razor clam Hai10Jing only has three sand groups), as can be seen from Figure 6, the lithological change of this inverting section reflection is abundant in content, vertically layering is good, and lateral resolution is higher.Finally by field of razor clam Hai1Jing and field of razor clam Hai10Jing, set out, successively follow the trail of the top bottom boundary of three sand groups, three sand groups discrimination power on the inverting section is higher, thereby has improved the accuracy of three sand component cloth scopes of the prediction Permian system.
In general, on the Study of The Impedence Inversion Restrained By Well Log model, large few content that has substantially reflected sand body of wave impedance value in the uniform thickness situation, the wave impedance value is higher, and sand body content is more, and content of mudstone is less.According to this conclusion, on the Study of The Impedence Inversion Restrained By Well Log model, three sand components have been got to the overall distribution that its total wave impedance value reflects sand body indescribably, Fig. 7 is the total wave impedance value of the Permian system one sand group planimetric map, as can be seen from Figure 7, the northern total wave impedance value in work area is larger, reduces gradually southwards.Illustrate that the northern sand body in work area is thicker, be thinned to gradually southwards the sand body pinching.
1), read 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 for example, to set wave impedance value (10000m/s.g/cm
3) be boundary the sand thickness that reads sample point one sand group according to the sand of adding up in the geology log data, wave impedance of mudstone value.In addition, can also calculate accordingly sandy ground than check possibly the rationality of setting wave impedance value value in order to the later stage.
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, determine the sand body of sample point one sand group and total wave impedance boundary value of mud stone (100000m/s.g/cm for example
3), and take total wave impedance boundary value and delimit sand body pinching place of sample point one sand group as boundary.
4), by total wave impedance value of sample point one sand group, deduct the total wave impedance value of sand body that total wave impedance boundary value obtains sample point one sand group.
5), repeat sand thickness and the total wave impedance value of sand body that above step is obtained respectively other sample point one sand groups in 20 sampled points.Obtain the Permian system one a sand group sampled point data statistic as shown in table 1 below.
Table 1 Permian system one sand group sampled point data statistic
6), the sand thickness by the sand group that reads at 20 sampled points and the total wave impedance value of sand body that the total wave impedance value of sand body simulates a sand group and the relation curve of its sand thickness, as shown in Figure 8.When guaranteeing that sand thickness is zero, the total wave impedance value of sand body is also zero, and the relation curve simulated must be crossed initial point.
7), repeat above step and obtain respectively the total wave impedance value of sand body of other sand groups in the target exploration district and the relation curve of its sand thickness.
In embodiments of the present invention, after obtaining corresponding geological data and geologic logging data, seismic properties that can also be corresponding according to seismic data acquisition (such as amplitude generic attribute and frequency generic attribute etc.), and the Sand-body Prediction distribution of obtaining described target exploration district according to seismic properties.Because amplitude attribute comprises that the sand body distribution situation that RMS amplitude attribute, average energy attribute and gross energy attribute show with it is close, and the arc length attribute in frequency attribute differs far away, therefore can the amplitude generic attribute while predicting the distribution of three sand groups of the Permian system as a reference.And, after step 104 obtains the sand thickness of corresponding sand group according to relation curve, for making result of calculation, more meet geologic rule, can revise according to the Sand-body Prediction distribution exceptional value and the anomaly trend of the sand thickness of corresponding sand group.
The method of utilizing Study of The Impedence Inversion Restrained By Well Log prediction sand thickness of the embodiment of the present invention is on the basis of Study of The Impedence Inversion Restrained By Well Log prediction, adopts sampled point statistical method prediction sand thickness.On the basis of Study of The Impedence Inversion Restrained By Well Log prediction, for each sand group in target exploration district, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness; Then obtain the sand thickness of corresponding sand group according to relation curve.Because the lithological change of Study of The Impedence Inversion Restrained By Well Log model reflection is abundant in content, vertically layering is good, lateral resolution is higher, and sampling number according to statistics rule according to some actual samples point data researches derivation sand body distribution and thickness thereof, thereby the identifiability of reservoir and the reliability of inversion prediction result are significantly improved, thereby it can improve for complex geologic conditions, reservoir thickness differ, physical property changes stratum greatly sand body and distribute and the reliability of thickness prediction.
Shown in Fig. 2, the device that utilizes Study of The Impedence Inversion Restrained By Well Log prediction sand thickness of the embodiment of the present invention comprises: 24(is optional for data acquisition module 21, wave impedance inversion module 22, sand thickness prediction module 23, sand body forecast of distribution module) and sand 25 body thickness correcting modules (optional).Wherein: data acquisition module 21 obtains corresponding geological data and geologic logging data for carrying out in the target exploration district after seismic prospecting drilling well sampling. and wave impedance inversion module 22 is for take the geologic logging data as constraint condition, according to geological data, the target exploration district is carried out to Study of The Impedence Inversion Restrained By Well Log, 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 the Study of The Impedence Inversion Restrained By Well Log model, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness for each sand group in target exploration district, obtain the sand thickness of corresponding sand group according to relation curve.Sand body forecast of distribution module 24 is for after obtaining corresponding geological data and geologic logging data, the corresponding seismic properties according to seismic data acquisition, and the Sand-body Prediction of obtaining the target exploration district according to seismic properties distributes.Because amplitude attribute comprises that the sand body distribution situation that RMS amplitude attribute, average energy attribute and gross energy attribute show with it is close, and the arc length attribute in frequency attribute differs far away, therefore can the amplitude generic attribute while predicting the distribution of three sand groups of the Permian system as a reference.Sand thickness correcting module 25, for after according to relation curve, obtaining the sand thickness of corresponding sand group, meets geologic rule more for making result of calculation, can revise according to the Sand-body Prediction distribution exceptional value and the anomaly trend of the sand thickness of corresponding sand group.
As can be seen from Figure 2, sand thickness prediction module 23 specifically comprises that again the sampled point sand thickness obtains submodule 231, the total wave impedance of sand group and obtain submodule 232, the total wave impedance boundary value of sand group and determine that submodule 233, the total wave impedance value of sand body calculating sub module 234, first control submodule 235, relation curve matching submodule 236, sand thickness and obtain submodule 237 and second and control submodule 238.Wherein, the sampled point sand thickness is obtained submodule 231 for reading the zone thickness of a certain sample point one sand group on the inverting section at the Study of The Impedence Inversion Restrained By Well Log model, then for example, to set wave impedance value (10000m/s.g/cm
3) be boundary the sand thickness that reads sample point one sand group according to the sand of adding up in the geology log data, wave impedance of mudstone value.The total wave impedance of sand group is obtained the total wave impedance value that reads sample point one sand group on total wave impedance value planimetric map of submodule 232 for the sand group at the Study of The Impedence Inversion Restrained By Well Log model.The total wave impedance boundary value of sand group determines that 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 take total wave impedance boundary value and delimit sand body pinching place of sample point one sand group as boundary.The total wave impedance value of sand body calculating sub module 234 deducts for the total wave impedance value by sample point one sand group the total wave impedance value of sand body that total wave impedance boundary value obtains sample point one sand group.First controls submodule 235 obtains for controlling the sampled point sand thickness that submodule 231, the total wave impedance of sand group obtain submodule 232, the total wave impedance boundary value of sand group is determined submodule 233 and sand thickness and the sand body total wave impedance value of the total wave impedance value of sand body calculating sub module 234 repeated works to obtain respectively other sample point one sand groups of several sampled points.The total wave impedance value of sand body that relation curve matching submodule 236 simulates a sand group for sand thickness and the total wave impedance value of sand body of the sand group by reading at several sampled points and the relation curve of its sand thickness, wherein relation curve is crossed initial point.Sand thickness is obtained submodule 237 for obtain the sand thickness of corresponding sand group according to relation curve.Second controls submodule 238 obtains submodule 237 repeated works with the total wave impedance value of sand body of obtaining respectively other sand groups of target exploration district and the relation curve of its sand thickness for controlling the first control submodule 235, relation curve matching submodule 236 and sand thickness.
In the device that utilizes Study of The Impedence Inversion Restrained By Well Log prediction sand thickness of the embodiment of the present invention, on the basis of Study of The Impedence Inversion Restrained By Well Log prediction, adopt sampled point statistical method prediction sand thickness.On the basis of Study of The Impedence Inversion Restrained By Well Log prediction, for each sand group in target exploration district, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness; Then obtain the sand thickness of corresponding sand group according to relation curve.Because the lithological change of Study of The Impedence Inversion Restrained By Well Log model reflection is abundant in content, vertically layering is good, lateral resolution is higher, and sampling number according to statistics rule according to some actual samples point data researches derivation sand body distribution and thickness thereof, thereby the identifiability of reservoir and the reliability of inversion prediction result are significantly improved, thereby it can improve for complex geologic conditions, reservoir thickness differ, physical property changes stratum greatly sand body and distribute and the reliability of thickness prediction.
Above-described specific embodiment; purpose of the present invention, technical scheme and beneficial effect are further described; institute is 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 of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.
Claims (10)
1. a method of utilizing Study of The Impedence Inversion Restrained By Well Log prediction sand thickness, is characterized in that, comprises the following steps:
After carrying out seismic prospecting drilling well sampling, the target exploration district obtains corresponding geological data and geologic logging data;
The described geologic logging data of take are constraint condition, according to described geological data, described target exploration district are carried out to Study of The Impedence Inversion Restrained By Well Log, 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, for each sand group in described target exploration district, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness;
Obtain the sand thickness of corresponding sand group according to described relation curve.
2. method according to claim 1, it is characterized in that, on the described section of the inverting at the Study of The Impedence Inversion Restrained By Well Log model, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness for each sand group in described target exploration district, specifically comprise the steps:
Read the zone thickness of a certain sample point one sand group on the inverting section of described Study of The Impedence Inversion Restrained By Well Log model, then take the sand thickness of setting the wave impedance value as boundary and reading described sample point one sand group according to the sand of adding up in described geologic logging data, wave impedance of mudstone value;
Read total wave impedance value of described sample point one sand group on total wave impedance value planimetric map of the described sand group of described Study of The Impedence Inversion Restrained By Well Log model;
By the total wave impedance value planimetric map to a described sand group and the comparative analysis of described geologic logging data, determine the sand body of described sample point one sand group and total wave impedance boundary value of mud stone, and take described total wave impedance boundary value and delimit sand body pinching place of described sample point one sand group as boundary;
Deduct by total wave impedance value of described sample point one sand group the total wave impedance value of sand body that described total wave impedance boundary value obtains described sample point one sand group;
Repeat above step and obtain respectively sand thickness and the total wave impedance value of sand body of other sample point one sand groups in described several sampled points;
Sand thickness by the sand group that reads at described several sampled points and the total wave impedance value of sand body that the total wave impedance value of sand body simulates a described sand group and the relation curve of its sand thickness, wherein said relation curve is crossed initial point;
Repeat above step and obtain respectively 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.
3. method according to claim 1, is characterized in that, described obtain corresponding geological data and geologic logging data after, also comprise:
The corresponding seismic properties according to described seismic data acquisition, the Sand-body Prediction of obtaining described target exploration district according to described seismic properties distributes;
Described obtain the sand thickness of corresponding sand group according to relation curve after, according to described Sand-body Prediction distribute exceptional value and the anomaly trend of the sand thickness of revising described corresponding sand group.
4. method according to claim 3, 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 the arc length attribute.
5. method according to claim 2, 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.
6. a device that utilizes Study of The Impedence Inversion Restrained By Well Log prediction sand thickness, is characterized in that, comprising:
Data acquisition module, obtain corresponding geological data and geologic logging data for carrying out in the target exploration district after seismic prospecting drilling well sampling;
The wave impedance inversion module, be constraint condition for take described geologic logging data, according to described geological data, described target exploration district carried out to Study of The Impedence Inversion Restrained By Well Log, obtains the Study of The Impedence Inversion Restrained By Well Log model in described target exploration district;
The sand thickness prediction module, for the inverting section at described Study of The Impedence Inversion Restrained By Well Log model, evenly choose respectively the relation curve of several sampled points with the total wave impedance value of the sand body that simulates corresponding sand group and its sand thickness for each sand group in described target exploration district, obtain the sand thickness of corresponding sand group according to described relation curve.
7. device according to claim 6, is characterized in that, described sand thickness prediction module specifically comprises:
The sampled point sand thickness is obtained submodule, for reading the zone thickness of a certain sample point one sand group on the inverting section of described Study of The Impedence Inversion Restrained By Well Log model, then take the sand thickness of setting the wave impedance value as boundary and reading described sample point one sand group according to the sand of adding up in described geologic logging data, wave impedance of mudstone value;
The total wave impedance of sand group is obtained submodule, reads total wave impedance value of described sample point one sand group on the total wave impedance value planimetric map for the described sand group at described Study of The Impedence Inversion Restrained By Well Log model;
The total wave impedance boundary value of sand group is determined 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 sand body of described sample point one sand group and total wave impedance boundary value of mud stone, and take described total wave impedance boundary value and delimit sand body pinching place of described sample point one sand group as boundary;
The total wave impedance value of sand body calculating sub module, deduct for the total wave impedance value by described sample point one sand group the total wave impedance value of sand body that described total wave impedance boundary value obtains described sample point one sand group;
First controls submodule, for controlling described sampled point sand thickness, obtains that submodule, the total wave impedance of described sand group obtain submodule, the total wave impedance boundary value of described sand group is determined submodule and sand thickness and the sand body total wave impedance value of the total wave impedance value of described sand body calculating sub module repeated work to obtain respectively other sample point one sand groups of described several sampled points;
Relation curve matching submodule, for the sand thickness by the sand group that reads at described several sampled points and the total wave impedance value of sand body that the total wave impedance value of sand body simulates a described sand group and the relation curve of its sand thickness, wherein said relation curve is crossed initial point;
Sand thickness is obtained submodule, for obtain the sand thickness of corresponding sand group according to described relation curve;
Second controls submodule, for controlling described first, controls submodule, described relation curve matching submodule and described sand thickness and obtains the submodule repeated work with the total wave impedance value of sand body of obtaining respectively other sand groups of described target exploration district and the relation curve of its sand thickness.
8. device according to claim 6, is characterized in that, also comprises:
Sand body forecast of distribution module, for described obtain corresponding geological data and geologic logging data after, the corresponding seismic properties according to described seismic data acquisition, the Sand-body Prediction of obtaining described target exploration district according to described seismic properties distributes;
The sand thickness correcting module, for described obtain the sand thickness of corresponding sand group according to relation curve after, according to described Sand-body Prediction distribute exceptional value and the anomaly trend of the sand thickness of revising described corresponding sand group.
9. device according to claim 8, 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 the arc length attribute.
10. device according to claim 7, 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 true CN103454685A (en) | 2013-12-18 |
| CN103454685B 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) |
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103954995A (en) * | 2014-04-22 | 2014-07-30 | 核工业北京地质研究院 | Sand body reorganization method in sandstone-type uranium deposit exploration |
| CN104090303A (en) * | 2014-07-10 | 2014-10-08 | 中国海洋石油总公司 | Seismic inversion method and device |
| CN104133250A (en) * | 2014-07-23 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司 | Gypsiferous salt stratum geological layering method |
| CN105353411A (en) * | 2015-10-27 | 2016-02-24 | 中国石油大学(华东) | Description method of reservoir sand bodies on ultra-deep coal layer |
| CN106772614A (en) * | 2017-02-28 | 2017-05-31 | 中国石油天然气股份有限公司 | The Forecasting Methodology of High-quality Reservoir in a kind of Gravity-flow Channel Sandbody |
| CN106772615A (en) * | 2016-11-09 | 2017-05-31 | 中国石油大学(华东) | A kind of broad-band teleseismic multiple domain combines AVO inversion methods |
| CN106842317A (en) * | 2016-12-22 | 2017-06-13 | 中国石油天然气集团公司 | A kind of method and device for predicting oil sand body distribution |
| CN107121699A (en) * | 2017-02-28 | 2017-09-01 | 中国石油化工股份有限公司 | A kind of sedimentary facies identification method under earthquake phase control |
| CN107272081A (en) * | 2017-07-18 | 2017-10-20 | 中国石油大学(华东) | A kind of submountain region extension of sedimentary facies belts horizon prediction method |
| CN107422374A (en) * | 2017-04-21 | 2017-12-01 | 中国石油大港油田勘探开发研究院 | A kind of sags rich in oil and gas slope belt fine granularing scalability 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 |
| CN107942405A (en) * | 2017-11-15 | 2018-04-20 | 中国石油化工股份有限公司 | The method for predicting thin sand-mud interbed sand body cumulative thickness |
| CN108931813A (en) * | 2017-05-26 | 2018-12-04 | 中国石油化工股份有限公司 | A method of identification reef lid, reef core, reef base |
| CN109752760A (en) * | 2017-11-01 | 2019-05-14 | 中国石油天然气股份有限公司 | The determination method and device of seismic properties confidence level |
| CN110542922A (en) * | 2019-03-05 | 2019-12-06 | 中国石油化工股份有限公司 | Thin mutual reservoir description method taking reconstructed sand-to-ground ratio as seismic description target |
| CN111101935A (en) * | 2019-12-25 | 2020-05-05 | 中海石油(中国)有限公司 | Oil shale prediction method under few-well condition |
| CN111665562A (en) * | 2020-07-20 | 2020-09-15 | 西南石油大学 | 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 |
| CN112037306A (en) * | 2020-07-14 | 2020-12-04 | 中国石油天然气集团有限公司 | Sand body thickness map drawing method and device, computer equipment and storage medium |
| CN113296150A (en) * | 2021-07-02 | 2021-08-24 | 清华大学 | High-dimensional closed-loop network seismic inversion method under logging constraint |
| CN113406703A (en) * | 2021-05-31 | 2021-09-17 | 中国石油天然气集团有限公司 | Method and device for predicting thickness of stratum in carboniferous work area |
| CN113703054A (en) * | 2021-09-03 | 2021-11-26 | 中国科学院地质与地球物理研究所 | Sedimentary facies evolution quantitative characterization method and system based on geology-geophysical |
| CN113703042A (en) * | 2020-05-21 | 2021-11-26 | 中国石油天然气股份有限公司 | Method and device for determining accumulated thickness of interlayer clastic rock sand body |
| CN114859418A (en) * | 2022-06-15 | 2022-08-05 | 北京阳光吉澳能源技术有限公司 | River channel sand reservoir identification method and device and electronic equipment |
| CN115903026A (en) * | 2023-01-09 | 2023-04-04 | 东北石油大学三亚海洋油气研究院 | Method, equipment and medium for analyzing configuration of composite sand body |
Citations (4)
| 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 |
| US20100309748A1 (en) * | 2009-06-04 | 2010-12-09 | Baker Hughes Incorporated | Method and Apparatus for Determining Radial Shear Velocity Variation From Dipole Acoustic Logging |
| 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 |
-
2013
- 2013-08-09 CN CN201310346850.2A patent/CN103454685B/en active Active
Patent Citations (4)
| 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 |
| US20100309748A1 (en) * | 2009-06-04 | 2010-12-09 | Baker Hughes Incorporated | Method and Apparatus for Determining Radial Shear Velocity Variation From Dipole Acoustic Logging |
| 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 |
|---|
| 师素珍 等: "多井约束反演在煤层厚度定量预测中的应用", 《采矿与安全工程学报》 * |
| 熊冉 等: "曲线重构反演在储层横向预测中的应用", 《西南石油大学学报（自然科学版）》 * |
| 蒋炼 等: "基于地震相分析的砂体储层厚度描述", 《断块油气田》 * |
| 赵长勋: "陈家庄南区馆陶组下段砂体识别方法研究", 《胜利油田职工大学学报》 * |
Cited By (34)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103954995A (en) * | 2014-04-22 | 2014-07-30 | 核工业北京地质研究院 | Sand body reorganization method in sandstone-type uranium deposit exploration |
| CN104090303A (en) * | 2014-07-10 | 2014-10-08 | 中国海洋石油总公司 | Seismic inversion method and device |
| CN104090303B (en) * | 2014-07-10 | 2016-11-02 | 中国海洋石油总公司 | A kind of seismic inversion method and device |
| CN104133250A (en) * | 2014-07-23 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司 | Gypsiferous salt stratum geological layering method |
| CN104133250B (en) * | 2014-07-23 | 2017-06-13 | 中国石油集团川庆钻探工程有限公司 | Gypsum-salt bed geology card layer method |
| CN105353411A (en) * | 2015-10-27 | 2016-02-24 | 中国石油大学(华东) | Description method of reservoir sand bodies on ultra-deep coal layer |
| CN106772615B (en) * | 2016-11-09 | 2018-11-13 | 中国石油大学(华东) | A kind of broad-band teleseismic multiple domain joint AVO inversion methods |
| CN106772615A (en) * | 2016-11-09 | 2017-05-31 | 中国石油大学(华东) | A kind of broad-band teleseismic multiple domain combines AVO inversion methods |
| CN106842317A (en) * | 2016-12-22 | 2017-06-13 | 中国石油天然气集团公司 | A kind of method and device for predicting oil sand body distribution |
| CN107121699A (en) * | 2017-02-28 | 2017-09-01 | 中国石油化工股份有限公司 | A kind of sedimentary facies identification method under earthquake phase control |
| CN107121699B (en) * | 2017-02-28 | 2019-06-21 | 中国石油化工股份有限公司 | A kind of sedimentary facies identification method under earthquake phase control |
| CN106772614A (en) * | 2017-02-28 | 2017-05-31 | 中国石油天然气股份有限公司 | The Forecasting Methodology of High-quality Reservoir in a kind of Gravity-flow Channel Sandbody |
| CN107422374A (en) * | 2017-04-21 | 2017-12-01 | 中国石油大港油田勘探开发研究院 | A kind of sags rich in oil and gas slope belt fine granularing scalability method |
| CN108931813A (en) * | 2017-05-26 | 2018-12-04 | 中国石油化工股份有限公司 | A method of identification reef lid, reef core, reef base |
| CN107272081A (en) * | 2017-07-18 | 2017-10-20 | 中国石油大学(华东) | 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 |
| CN109752760A (en) * | 2017-11-01 | 2019-05-14 | 中国石油天然气股份有限公司 | The determination method and device of seismic properties confidence level |
| CN109752760B (en) * | 2017-11-01 | 2021-01-01 | 中国石油天然气股份有限公司 | Method and device for determining reliability of seismic attribute |
| CN107942405A (en) * | 2017-11-15 | 2018-04-20 | 中国石油化工股份有限公司 | The method for predicting thin sand-mud interbed sand body cumulative thickness |
| CN110542922A (en) * | 2019-03-05 | 2019-12-06 | 中国石油化工股份有限公司 | Thin mutual reservoir description method taking reconstructed sand-to-ground ratio as seismic description target |
| CN111101935A (en) * | 2019-12-25 | 2020-05-05 | 中海石油(中国)有限公司 | Oil shale prediction method under few-well condition |
| CN111101935B (en) * | 2019-12-25 | 2023-01-03 | 中海石油(中国)有限公司 | Oil shale prediction method under few-well condition |
| CN113703042A (en) * | 2020-05-21 | 2021-11-26 | 中国石油天然气股份有限公司 | Method and device for determining accumulated thickness of interlayer clastic rock sand body |
| 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 |
| CN111665562A (en) * | 2020-07-20 | 2020-09-15 | 西南石油大学 | 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 |
| CN113296150A (en) * | 2021-07-02 | 2021-08-24 | 清华大学 | High-dimensional closed-loop network seismic inversion method under logging constraint |
| CN113703054A (en) * | 2021-09-03 | 2021-11-26 | 中国科学院地质与地球物理研究所 | Sedimentary facies evolution quantitative characterization method and system based on geology-geophysical |
| CN114859418A (en) * | 2022-06-15 | 2022-08-05 | 北京阳光吉澳能源技术有限公司 | River channel sand reservoir identification method and device and electronic equipment |
| CN114859418B (en) * | 2022-06-15 | 2022-12-20 | 北京阳光吉澳能源技术有限公司 | River channel sand reservoir identification method and device and electronic equipment |
| CN115903026A (en) * | 2023-01-09 | 2023-04-04 | 东北石油大学三亚海洋油气研究院 | Method, equipment and medium for analyzing configuration of composite sand body |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103454685B (en) | 2016-04-06 |
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 | |
| CN104047597B (en) | Oily mud shale stratum log standardized method | |
| CN103527184B (en) | A kind of Forecasting Methodology of dolostone reservoirs and system | |
| CN104977618B (en) | A kind of method evaluated shale gas reservoir and find dessert area | |
| Wang et al. | Cenozoic structure and tectonic evolution of the Kuqa fold belt, southern Tianshan, China | |
| CN101158724B (en) | Reservoir thickness prediction method based on dipolar wavelet | |
| CN104516018B (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
| CN109425896A (en) | Dolomite oil and gas reservoir distribution forecasting method and device | |
| CN105044770A (en) | Compact glutenite gas reservoir quantificational prediction method | |
| CN104267429A (en) | Method and device for determining formation pressure | |
| CN103670388A (en) | Method for evaluating organic carbon content of shale | |
| CN101634717A (en) | Fine shear-wave (S-wave) impedance access technology based on logging and prestack channel set seismic data | |
| CN104570067A (en) | Phase-controlled earthquake inversion method in geophysical exploration | |
| CN105005077A (en) | Thin layer thickness prediction method with combination of real drilling wells and virtual wells under rare well condition | |
| CN104049283A (en) | Method for identification and prediction of fore-land basin extremely-thick conglomerate body | |
| Fajana et al. | Seismic interpretation and petrophysical analysis for hydrocarbon resource evaluation of ‘Pennay’field, Niger Delta | |
| Xia et al. | Application of 3D fine seismic interpretation technique in Dawangzhuang area, Bohai Bay Basin, Northeast China | |
| CN104142516A (en) | Method for predicting thickness of thin single sand bed | |
| Khan et al. | 3D structural modeling integrated with seismic attribute and petrophysical evaluation for hydrocarbon prospecting at the Dhulian Oilfield, Pakistan | |
| CN110187389A (en) | A kind of AVA inversion method based on thin bed reflection theory | |
| CN103278852B (en) | Utilize the method for seismic data volume waveform configuration characteristic model predicting oil/gas | |
| Abdelrahman et al. | Seismic vulnerability assessment in the new urban area of Diriyah Governorate, Riyadh, Saudi Arabia | |
| Wylie Jr et al. | Log-curve amplitude slicing: Visualization of log data and depositional trends in the Middle Devonian Traverse Group, Michigan basin, United States | |
| Asmari et al. | An integrated geophysical study for inferring the subsurface setting of water-bearing aquifers: case study at Wadi al-Dawasir area, Saudi Arabia | |
| Gargouri | Multicomponent 3D seismic interpretation of the Marcellus shale Bradford county, Pennsylvania |
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 |