CN108121009B - Development late stage complex fault block velocity-variant depth mapping method - Google Patents
Development late stage complex fault block velocity-variant depth mapping method Download PDFInfo
- Publication number
- CN108121009B CN108121009B CN201611079816.3A CN201611079816A CN108121009B CN 108121009 B CN108121009 B CN 108121009B CN 201611079816 A CN201611079816 A CN 201611079816A CN 108121009 B CN108121009 B CN 108121009B
- Authority
- CN
- China
- Prior art keywords
- well
- velocity
- interpolation
- interval velocity
- point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000013507 mapping Methods 0.000 title claims abstract description 21
- 238000011161 development Methods 0.000 title claims abstract description 18
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 8
- 238000012546 transfer Methods 0.000 claims abstract description 7
- 238000003325 tomography Methods 0.000 claims description 11
- 239000000284 extract Substances 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- 238000005215 recombination Methods 0.000 abstract 1
- 230000006798 recombination Effects 0.000 abstract 1
- 244000264242 Descurainia sophia Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of development late stage complex fault block velocity-variant depth mapping method, comprising: step 1, the interval velocity of individual well is obtained by seismic interpretation layer position and geological layering；Step 2, being superimposed along layer seismic properties, and attributed graph with T0 figure is waited for variation of lithological can be reflected by extracting；Step 3, Wu Jing is encrypted based on the consistent principle of formation velocity at identical lithology and same time and well is unevenly distributed the speed point in region；Step 4, using seismic interpretation layer position be trend surface and fault boundary is that interpolation boundary carries out interpolation to individual well interval velocity, establishes two-dimensional layer velocity field；Step 5, it is multiplied using the two-dimentional interval velocity grid of generation with equal T0 grid and does time and depth transfer, obtain final structural map.This method provides accurate basic guide for Study on residual oil distribution and the well pattern recombination of complex fault block meticulous pool description, has directive function to the later development of such oil reservoir, popularization and application foreground is wide, remarkable economic and social benefits.
Description
Technical field
The present invention relates to complex fault block meticulous pool description fields, especially relate to a kind of development late stage complex fault block change
Fast tectonic ore-forming method.
Background technique
Complex fault block by lithology, terrain configuration and fracture system due to being influenced, and in more regions, formation velocity is longitudinal, horizontal
It is larger to velocity variations, and development late stage is higher to tectonic ore-forming precision, is fitted using individual well average speed or more wells public
The error that formula carries out on the structural map and well of time and depth transfer is larger, is not able to satisfy required precision, needs using velocity-variant depth mapping.
It is relatively more to the research of variable speed mapping method at present, mainly have from well point speed interpolation velocity field and seeks three from stack velocity
Tie up average velocity field both.Wherein divide uneven or well spacing density diluter from the method for individual well speed interpolation velocity field in well
The velocity variations rule of underground cannot accurately be showed in thin work area；From stack velocity field computation three-dimensional average velocity in Wu Jing
In the case of more can truly show the changing rule of subsurface velocities field.However, this calculate three-dimensional velocity from stack velocity
The method of field is more troublesome, and structure interpretation personnel are that do not have stack velocity data under normal circumstances, and what can be done is exactly from well
Spot speed calculating speed field.Although the fault block oil reservoir well of development late stage is more, the method for well point speed interpolation velocity field is in well
Less, well is unevenly distributed, there are the local interpolation precision of tomography is inadequate；Therefore, this in the place that craven fault and stratum change
The speed that method calculates cannot really show subsurface velocities variation characteristic, structural configuration can be caused to be distorted after time and depth transfer.For
This solves the above technical problem we have invented a kind of new development late stage complex fault block velocity-variant depth mapping method.
Summary of the invention
The object of the present invention is to provide it is a kind of using with the encryption of the well interval velocity of fixed well same time and lithology without well and
Well is unevenly distributed the speed point in region, and well point speed difference velocity field under seismic horizon and fault boundary constraint is opened
Send out later period complex fault block velocity-variant depth mapping method.
The purpose of the present invention can be achieved by the following technical measures: development late stage complex fault block velocity-variant depth mapping side
Method, which includes: step 1, is obtained by seismic interpretation layer position and geological layering
Obtain the interval velocity of individual well；Step 2, being superimposed along layer seismic properties, and attributed graph with T0 figure is waited for variation of lithological can be reflected by extracting；
Step 3, Wu Jing is encrypted based on the consistent principle of formation velocity at identical lithology and same time and well is unevenly distributed region
Speed point；Step 4, using seismic interpretation layer position be trend surface and fault boundary to be that interpolation boundary carries out individual well interval velocity slotting
Value, establishes two-dimensional layer velocity field；Step 5, it is multiplied using the two-dimentional interval velocity grid of generation with equal T0 grid and does time and depth transfer, obtained
To final structural map.
The purpose of the present invention can be also achieved by the following technical measures:
In step 1, on the basis of completing whole district's fine reservoir prediction, layer position and tomography closure explanation, reading layer position
The time value of the geologic horizon corresponding seismic interpretation layer position for the individual well that calibration is accurate and hole deviation is errorless, and utilize the layering on well
Depth and the time value of reading calculate individual well interval velocity.
In step 2, to research layer variation of lithological analyze, and extract can reflect variation of lithological along layer earthquake category
Property, while the T0 figure such as generation, attributed graph and equal T0 figure overlapping.
In step 3, local enciphering rate point unevenly distributed to no well or well, which sought in the whole district
Look for and need interval velocity of the interval velocity of the lithology and seismic horizon time all identical fixed well at encryption as the point.
In step 4, it is trend surface using equal T0 figure, using fault line as interpolation boundary, two dimension is carried out to individual well point interval velocity
Gridding generates two-dimentional interval velocity；Using minimum trend gridding interpolation method during interpolation, using well point speed as seed
Point, using seismic interpretation layer position as trend between two o'clock, when encountering fault line, interpolation terminates.
Development late stage complex fault block velocity-variant depth mapping method in the present invention, utilizes variation of lithological and seismic interpretation layer position
With the variable speed mapping method of the well point interpolation speed field under tomography co- controlling, depth relationship and the completion whole district when implementing the whole district
It on the basis of structure interpretation, is unevenly distributed for the lateral speed change of complex fault block, well, more than tomography and without stack velocity,
It is mainly the variation due to caused by lithology, structural configuration and tomography from analysis complex fault block subsurface velocities, utilizes identical rock
The speed of the fixed well of property and same time is unevenly distributed the speed point of block to encrypt Wu Jing and well, and with seismic interpretation layer
Position and fault boundary are constraint condition interpolation individual well interval velocity, to carry out variable speed mapping.This method solve due to well quantity
It is limited, be unevenly distributed, tomography it is more, the problem of velocity field inaccuracy is sought using individual well interval velocity, is not being folded by this method
The cross directional variations rule that more can truly show underground horizon velocity in the case where acceleration, is especially distributed not in well
, the difference of speed can be highlighted in the case that tomography is more, formation variation is fast.The inventive method can effectively improve complex fault block speed change
Mapping precision, and it is simple and practical.The technological invention is the Study on residual oil distribution of complex fault block meticulous pool description and opens
It sends out well pattern reformation and accurate basic guide is provided, there is directive function, popularization and application foreground to the later development of such oil reservoir
It is wide, remarkable economic and social benefits.
Detailed description of the invention
Fig. 1 is the flow chart of a specific embodiment of development late stage complex fault block velocity-variant depth mapping method of the invention；
Deep friendship remittance analysis chart when Fig. 2 is the 84 mouthfuls of wells in work area in a specific embodiment of the invention；
Fig. 3 is the RMS amplitude attribute and interpretation horizon that 10MS is extracted above and below layer in a specific embodiment of the invention
Equal T0 overlapping figure；
Fig. 4 is enciphering rate point position and corresponding fixed well signal and overlapping (attribute in a specific embodiment of the invention
With equal T0) figure；
Fig. 5 is two-dimentional interval velocity grid chart in a specific embodiment of the invention；
Fig. 6 is variable speed structure figure in a specific embodiment of the invention.
Specific embodiment
To enable above and other objects, features and advantages of the invention to be clearer and more comprehensible, preferably implementation is cited below particularly out
Example, and cooperate shown in attached drawing, it is described in detail below.
As shown in FIG. 1, FIG. 1 is the flow charts of development late stage complex fault block velocity-variant depth mapping method of the invention.
Step 101, on the basis of completing whole district's fine reservoir prediction, layer position and tomography closure explanation, reading layer position mark
The time value of the geologic horizon corresponding seismic interpretation layer position of true and errorless hole deviation individual well is fixed, and deep using the layering on well
Degree and the time value read calculate individual well interval velocity, and step enters 102.
Step 102, research layer variation of lithological is analyzed, and extract can reflect variation of lithological along layer seismic properties,
The T0 such as generation figure simultaneously, attributed graph and equal T0 figure overlapping, step enters 103.
Step 103, local enciphering rate point unevenly distributed to no well or well, the encryption method are found in the whole district
With need interval velocity of the interval velocity of lithology and seismic horizon time all identical fixed well at encryption as the point, step into
Enter 104.
Step 104, be trend surface using equal T0 figure, fault boundary is that interpolation boundary carries out interpolation to the interval velocity of individual well point,
Interval velocity grid is generated, step enters 105.
Step 105, interval velocity grid is multiplied with equal T0 grid and does time and depth transfer, process terminates.
The technology of the present invention is on the basis of completing whole district's fine reservoir prediction and whole district's closed structure is explained, after exploitation
Phase complex fault-block structure is at figure problems faced and needs, and the velocity field obtained in view of conventional individual well interpolation cannot really show
The variation of complex fault block underground horizon velocity, place unevenly distributed to Wu Jing and well is using known to identical lithology and time
Well point speed is encrypted as the spot speed, while with the variation tendency of seismic interpretation layer position being constraint, with fault boundary being
Interpolation boundary carries out interpolation to whole district's speed point；It ensures that there is speed point control in the whole district in this way, is establishing interval velocity field
When using well point speed as seed point, interpolation is carried out with the variation tendency of seismic interpretation layer position between two well points, it is disconnected when encountering
Interpolation terminates when layer line.The present invention utilizes the formation velocity one of the same seismic reflection layer position of identical lithology and same time
Then the principle of cause, place unevenly distributed to Wu Jing and well recycle seismic interpretation with fixed well speed enciphering rate point
The time trend of layer position constrains individual well interval velocity interpolation, and well point interpolation grid is enabled more truly to show underground
The variation of interval velocity field, and it is simple and practical.
In an application specific embodiment of the invention, example 1 is positioned at the northern Complex Fault Block Oil Reservoir of Gudong field, quilt
The fracture of two groups of the Near East north-south Xi Hejin is mutually cut, and multiple small fault blocks, research range 27Km2,84 mouthfuls of wells are formd.Study layer
Section SHAHEJIE FORMATION overlap of beds is coated on the buried hill Gu Dong, and formation thickness changes greatly, and Oil Reservoir Types are predominantly by tomography control
The fault block of system and disconnected nose oil reservoir.Presently, there are the problem of mainly have and verify oil domain and current construction mismatch, fault block is drawn
Point not fine, the nose that breaks, micro-structure development but fulfilling degree are low etc..
(1), deep when being done to all wells of the whole district's horizon calibration on the basis of whole district's horizon calibration and fine structural interpretation
Cross figure (Fig. 2) and browsed well profile analysis, check whether there is or not calibration inaccuracy and there are the wells of hole deviation.Then reading layer position is marked
Time value of the geological layering of true and errorless hole deviation individual well on seismic interpretation layer position is fixed, and utilizes the depth of seam division on well
And the time value read calculates individual well interval velocity.
(2), the lithologic character of earthquake interpretation horizon is analyzed, and extracts and can reflects that this layer of lithology sea level changes are special
The RMS amplitude attribute of sign, while the T0 isograms such as interpretation horizon generation and attributed graph are overlapped (Fig. 3).
(3) to no well area and well region unevenly distributed, lithology identical as from this and identical is found out from overlapping figure
Speed enciphering rate point (Fig. 4, table 1) of the interval velocity of the fixed well of time as the point.
Enciphering rate point and fixed well information table in 1 example 1 of table
(4) it is trend surface using equal T0 figure, using fault line as interpolation boundary, two-dimensional mesh is carried out to individual well point interval velocity and is formatted,
Generate two-dimentional interval velocity (Fig. 5).Using minimum trend gridding interpolation method during interpolation, using well point speed as seed
Point, using seismic interpretation layer position as trend between two o'clock, when encountering fault line, interpolation terminates.
(5) it is multiplied using the two-dimentional interval velocity grid of generation with equal T0 grid and does time and depth transfer, obtain final structural map
(Fig. 6).
Wu Jing is encrypted by the fixed well spot speed using identical lithology and time and well is unevenly distributed the speed in region
Point, and more really show using seismic interpretation layer position and tomography as the individual well speed field energy of constraint condition interpolation the variation of speed
Rule, the method is coincide preferably with achieving in variable speed mapping in example 1 preferably with effect with well point depth, convex
The difference in height of micro-structure and craven fault upper lower burrs is shown.By being compared to whole district's well depth of seam division and variable speed structure figure depth
It analyzes (table 2).
Whole district's well depth of seam division and variable speed mapping depth correlation analytical table in 2 example 1 of table
84 mouthfuls of wells of the whole district, error are up to 5.71 meters, minimal error 0；Wherein error is greater than 5 meters of Jing5Kou, accounts for
5.9%, Jing46Kou of the error between 5-1 meters, account for 54.8%, and Jing33Kou of the error less than 1 meter accounts for 39.3%.Fitting formula
The worst error of depth is 25 meters on tectonic ore-forming and well, and error is more greater than 5 meters of well, is mentioned using the variable speed mapping method
High structural map precision.Development late stage complex fault-block structure is 5 meters at most interval is schemed, and error is less than in the invention method
5 meters of probability is 94.1%, and precision is higher, can satisfy the needs of current complex fault-block structure mapping precision, and correct amount is smaller,
And this drawing methods speed is fast, easy to operate.
Development late stage complex fault block velocity-variant depth mapping method in the present invention, is by seismic interpretation layer position and ground first
Matter layering obtains the interval velocity of individual well；Then to explanation of seismic layer position extract can reflect variation of lithological along layer seismic properties, and
Attributed graph is superimposed with equal T0 figure；Again to the equally distributed area of no well, found and lithology and time from pass point from the whole district
Speed of the interval velocity of all identical fixed well as the point；It is finally trend surface and fault boundary is using seismic interpretation layer position
Interpolation boundary carries out interpolation to individual well interval velocity, two-dimensional layer velocity field is established, to reach complex fault-block structure variable speed mapping
Purpose.
Claims (1)
1. development late stage complex fault block velocity-variant depth mapping method, which is characterized in that the development late stage complex fault block variable speed structure
Drawing methods include:
Step 1, the interval velocity of individual well is obtained by seismic interpretation layer position and geological layering；
Step 2, being superimposed along layer seismic properties, and attributed graph with T0 figure is waited for variation of lithological can be reflected by extracting；
Step 3, Wu Jing is encrypted based on the consistent principle of formation velocity at identical lithology and same time and well is unevenly distributed
The speed point in region；
Step 4, using seismic interpretation layer position be trend surface and fault boundary is that interpolation boundary carries out interpolation to individual well interval velocity, is built
Vertical two-dimensional layer velocity field；
Step 5, it is multiplied using the two-dimentional interval velocity grid of generation with equal T0 grid and does time and depth transfer, obtain final structural map；
In step 1, on the basis of completing whole district's fine reservoir prediction, layer position and tomography closure explanation, horizon calibration is read
The time value of the geologic horizon corresponding seismic interpretation layer position of accurate and errorless hole deviation individual well, and utilize the depth of seam division on well
And the time value read calculates individual well interval velocity；
In step 2, to research layer variation of lithological analyze, and extract can reflect variation of lithological along layer seismic properties, together
The T0 such as Shi Shengcheng figure attributed graph and waits T0 figure to overlap；
In step 3, local enciphering rate point unevenly distributed to no well or well, the encryption method be the whole district find with
Interval velocity of the interval velocity of the lithology and seismic horizon time all identical fixed well that need to encrypt place as the point；
In step 4, it is trend surface using equal T0 figure, using fault line as interpolation boundary, two-dimensional grid is carried out to individual well point interval velocity
Change, generates two-dimentional interval velocity；Using minimum trend gridding interpolation method during interpolation, using well point speed as seed point,
Using seismic interpretation layer position as trend between two o'clock, when encountering fault line, interpolation terminates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611079816.3A CN108121009B (en) | 2016-11-30 | 2016-11-30 | Development late stage complex fault block velocity-variant depth mapping method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611079816.3A CN108121009B (en) | 2016-11-30 | 2016-11-30 | Development late stage complex fault block velocity-variant depth mapping method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108121009A CN108121009A (en) | 2018-06-05 |
| CN108121009B true CN108121009B (en) | 2019-07-12 |
Family
ID=62227107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611079816.3A Active CN108121009B (en) | 2016-11-30 | 2016-11-30 | Development late stage complex fault block velocity-variant depth mapping method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108121009B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110888163B (en) * | 2018-09-10 | 2022-01-04 | 中国石油天然气股份有限公司 | Method and device for determining development position of micro-fracture |
| CN110673209B (en) * | 2019-10-13 | 2021-06-04 | 东北石油大学 | Well-seismic calibration method |
| CN111308549B (en) * | 2020-03-09 | 2022-01-11 | 中国石油化工股份有限公司 | Variable-speed mapping method based on model inversion |
| CN112068199B (en) * | 2020-09-07 | 2021-04-09 | 北京中恒利华石油技术研究所 | Plane section rapid fault interpretation method |
| CN113325474B (en) * | 2021-06-03 | 2022-05-13 | 西南石油大学 | Method for discriminating biological reef |
| CN114355443A (en) * | 2021-11-25 | 2022-04-15 | 中国石油天然气股份有限公司 | Method and device for predicting residual oil under close well pattern condition |
| CN116381782A (en) * | 2023-02-17 | 2023-07-04 | 中海石油(中国)有限公司深圳分公司 | Low-amplitude structured graph forming method, device and storage device |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103885083A (en) * | 2014-01-29 | 2014-06-25 | 上海石油天然气有限公司 | Paleotectonic map mapping method and device thereof |
| CN104111476A (en) * | 2014-07-31 | 2014-10-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method and device for building formation velocity field |
| CN104678440A (en) * | 2015-02-15 | 2015-06-03 | 山东科技大学 | Well-constrained two-dimensional seismic variable velocity field nonlinear error correction method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2755059B1 (en) * | 2013-01-15 | 2021-11-03 | CGG Services SAS | Seismic data processing including data-constrained surface-consistent correction |
-
2016
- 2016-11-30 CN CN201611079816.3A patent/CN108121009B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103885083A (en) * | 2014-01-29 | 2014-06-25 | 上海石油天然气有限公司 | Paleotectonic map mapping method and device thereof |
| CN104111476A (en) * | 2014-07-31 | 2014-10-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method and device for building formation velocity field |
| CN104678440A (en) * | 2015-02-15 | 2015-06-03 | 山东科技大学 | Well-constrained two-dimensional seismic variable velocity field nonlinear error correction method |
Non-Patent Citations (1)
| Title |
|---|
| Depth Express变速成图在段六拨油田的应用;姜玲玲;《中国石油和化工标准与质量》;20161031(第18期);第84-85页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108121009A (en) | 2018-06-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108121009B (en) | Development late stage complex fault block velocity-variant depth mapping method | |
| CN104698497B (en) | A kind of fine calibration of thin sand-mud interbed layer position | |
| Hennings et al. | Stability of the fault systems that host‐induced earthquakes in the Delaware Basin of West Texas and Southeast New Mexico | |
| CN103454686A (en) | Method and system for reservoir prediction based on small-scale sedimentary facies of stratigraphic slice | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| CN109541685B (en) | River channel sand body identification method | |
| CN104375172B (en) | Volcanic underlayer structural configuration correct imaging method | |
| Abdel-Fattah et al. | 3D geometric modeling of the Abu Madi reservoirs and its implication on the gas development in Baltim area (Offshore Nile Delta, Egypt) | |
| Xia et al. | Application of 3D fine seismic interpretation technique in Dawangzhuang area, Bohai Bay Basin, Northeast China | |
| CN106226841A (en) | A kind of fluvial facies three-dimensional sedimentary facies model Decided modelling method | |
| Horne et al. | Structural characteristics of shallow faults in the Delaware Basin | |
| CN107269269A (en) | The quick look method of unconformable strata layer position | |
| CN105116447A (en) | A geological river direction discrimination method based on curvature-abnormal stripes | |
| Qi et al. | Hierarchy and subsurface correlation of muddy baffles in lacustrine delta fronts: a case study in the X Oilfield, Subei Basin, China | |
| Kim et al. | Lithofacies modeling by multipoint statistics and economic evaluation by NPV volume for the early Cretaceous Wabiskaw Member in Athabasca oilsands area, Canada | |
| CN109324171A (en) | A kind of sedimentary facies quantitative identification method based on lithology statistics | |
| Louie et al. | Advanced seismic imaging for geothermal development | |
| Singleton | Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas | |
| Zheng et al. | High-resolution sequence stratigraphic division and distribution of tidal deposits in the Zhuhai Formation, Huizhou Sag, Pearl River Mouth Basin, South China Sea | |
| Anderson | Structural relationships between surface lineaments and basement faulting in the northeastern Williston Basin | |
| Zhao et al. | Research on meandering river reservoir deposition architecture and 3D modeling of the Gudao Oil field in the Bohai Bay Basin | |
| Hagstrom et al. | Comparison of the morphology, facies, and reservoir quality of valley fills in the southern Athabasca Oil Sands Region, Alberta, Canada | |
| Hu et al. | Seismic interpretation of sandstone-type uranium deposits in the Songliao Basin, Northeast China | |
| Hohn et al. | Geostatistical analysis of oil production and potential using indicator kriging | |
| Duvbiama et al. | 3D Static Modelling of an offshore field in the Niger-delta |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |