CN106646608B - A kind of layer model method for building up for observation system parameter demonstration - Google Patents
A kind of layer model method for building up for observation system parameter demonstration Download PDFInfo
- Publication number
- CN106646608B CN106646608B CN201611240398.1A CN201611240398A CN106646608B CN 106646608 B CN106646608 B CN 106646608B CN 201611240398 A CN201611240398 A CN 201611240398A CN 106646608 B CN106646608 B CN 106646608B
- Authority
- CN
- China
- Prior art keywords
- layer
- root mean
- depth
- grid
- sequare velocity
- 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 25
- 238000004364 calculation method Methods 0.000 claims abstract description 3
- 238000005070 sampling Methods 0.000 claims description 17
- 238000009499 grossing Methods 0.000 claims description 8
- 238000005096 rolling process Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 7
- 238000007689 inspection Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 4
- 235000013399 edible fruits Nutrition 0.000 claims description 3
- 238000004513 sizing Methods 0.000 claims description 3
- 238000003892 spreading Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 238000013213 extrapolation Methods 0.000 abstract 1
- 238000009776 industrial production Methods 0.000 abstract 1
- 239000009671 shengli Substances 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 4
- 230000007812 deficiency Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000035800 maturation Effects 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- G01V20/00—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Abstract
The invention discloses it is a kind of for observation system parameter demonstration layer model method for building up, including：According to the existing formation information in work area and earthquake root mean sequare velocity, gridding is carried out to the layer position to be handled, depth value and inclination angle by all layers of position grid of interpolation calculation, the earthquake root mean sequare velocity of all discrete points is right to being converted into " depth speed " by " time speed ", reflection interval and root mean sequare velocity value are obtained according to the layer bit depth of the discrete point, all layers of position grid are traveled through to obtain the layer position grid of no reflection interval and root mean sequare velocity assignment using the equidistant inverse ratio weighted interpolation of progressively extrapolation based on grid index, final establish has depth, inclination angle, reflection interval and the unified layer model of root mean sequare velocity information.This method provides reliable basic data for the seismic observation system design in ripe oil-gas exploration area, available in the industrial production of seismic prospecting.
Description
Technical field
The present invention relates to oil-gas seismic exploration Data Acquisition Design fields, and middle level is proved more particularly to observation system parameter
The method that bit model is established.
Background technology
With deepening continuously for In Oil Field Exploration And Development, the geologic objective faced becomes increasingly complex, to seismic data into
As precision, resolution capability propose higher requirement.Seismic data acquisition is the basis of seismic prospecting, and observation system
Design plays the role of the quality of seismic data vital.In oil-gas exploration maturation zone, seismic prospecting has been introduced into
Secondary, even four acquisition phases three times, how according to the existing geological knowledge in work area and exploration results to new acquisition
It is the effective ways for improving exploration success ratio that observation system, which is designed,.
In parameter demonstration in observation system design, it is desirable to provide depth, inclination angle, reflection interval and the root mean square on stratum
The information such as speed, and do not include the depth information on stratum in earthquake root mean sequare velocity file, while in stratigraphic horizon file also not
Root mean sequare velocity information including stratum, therefore can not convenient, effectively establish the unification for observation system parameter demonstration
Layer model causes current observation system to choose the parameter demonstration that indivedual " proving a little " is observed system in designing mostly, and
The non-parameter demonstration established layer model and carry out entire purpose of horizon.On how to establish with depth, inclination angle, reflection interval and
The method of the unified layer model of root mean sequare velocity, has no that pertinent literature is delivered.
The content of the invention
In order to overcome the above-mentioned deficiencies of the prior art, join the present invention provides a kind of convenience, reliably for observation system
The unified layer model method for building up of number theory card.
The technical solution adopted in the present invention is as follows.
A kind of layer model method for building up for observation system parameter demonstration includes the following steps.
Step 1：It determines the stratigraphic distribution scope of target zone, obtains the layer position data of target zone and the earthquake root mean square in the work area
Speed data.
Step 2：The sizing grid designed according to exploration deployment sets the sizing grid of demonstration and position carries out layer by layer to purpose
Gridding is handled.
Step 3：According to existing purpose, position data carry out assignment to the layer position trellis depth with the sampling of layer position layer by layer.
Step 4:The layer position of layer position sampling is not had using the calculating of equidistant inverse ratio weight interpolation algorithm according to existing depth value
Trellis depth value.
Step 5：Using depth value as colour code, target zone depth model is drawn, checks the purpose effect of bit depth modeling layer by layer,
If depth model value needs to be smoothed depth model there are apparent zigzag.
Step 6：According to the inclination angle of all layers of position grid of each mesh coordinate and depth calculation.
Step 7：Using inclination value as colour code, purpose inclination layer model, inspection layer position inclination angle modeling effect, if inclined are drawn
Angle mould offset then needs to be smoothed inclination angle model there are apparent zigzag.
Step 8：By the data of earthquake root mean sequare velocity discrete sampling point by " when m- speed " to being converted into " depth-speed
Degree " is right.
Step 9：The depth of layer position grid obtains the reflection of the grid according to where earthquake root mean sequare velocity discrete sampling point
Time and root mean sequare velocity.
Step 10：All layers of position grid are traveled through, to the grid that no reflection interval and root mean sequare velocity sample using equidistantly
From reflection interval and the root mean sequare velocity that this layer of position grid is calculated in inverse ratio weight interpolation algorithm.
Step 11：Using reflection interval as colour code, target zone reflection interval model, the modeling of inspection layer position reflection interval are drawn
Effect, if reflection interval model value needs to be smoothed reflection interval model there are apparent zigzag.
Step 12：Using root mean sequare velocity as colour code, target zone root mean sequare velocity model, inspection layer position root mean square speed are drawn
Degree modeling effect, if root mean sequare velocity model value needs to put down root mean sequare velocity model there are apparent zigzag
Sliding processing.
Step 13：The depth of each grid of target zone, inclination angle, reflection interval and root mean sequare velocity is defeated according to grid index
Go out, that is, complete the foundation that the target zone unifies layer model.
Further, in the step 1, the layer position data of acquisition be derived from work area three dimensional seismic data geologic interpretation into
Fruit is exported after spatial spreading samples including sample point coordinate, the text file of layer locating depth angle value.
Further, in the step 1, the earthquake root mean square data of acquisition are derived from work area D seismic data processing
Achievement, be exported after spatial spreading samples include sample point coordinate, " when m- speed " to text file.
Further, in the step 4, without the grid of layer position data sampling, depth value is added by equidistant inverse ratio
Power interpolation obtains, and the search strategy in Interpolation Process is to carry out horizontal and vertical search according to grid index, and according to default
Control points the problem of progressively extrapolating, avoiding the inefficiency that global search is brought；The wherein default several evidences in control point
Layer position sampled point density degree determines.
Further, in the step 5, when smoothing processing algorithm uses square neighborhood rolling average algorithm to save search
Between.
Further, in the step 6, the scope at inclination angle is defined as [0,90 °], and grid phase is used and calculated when specifically calculating
The plane of the central point composition of three layer position grids of distribution adjacent and triangular in shape, each layer position grid and this plane plane with
The angle of horizontal plane is the inclination angle for calculating grid.
Further, in the step 7, when smoothing processing algorithm uses square neighborhood rolling average algorithm to save search
Between.
Further, in the step 8, " when m- speed " to " depth-speed " to conversion be directly with the time with
Root speed is multiplied to obtain depth, and root mean sequare velocity remains unchanged.
Further, in the step 9, the speed of layer position grid is to respective value according to its depth in " depth-speed "
What former and later two sample value linear interpolations calculated.
Further, in the step 10, without reflection interval and the grid of root mean sequare velocity data sampling, reflection interval
It is obtained with root mean sequare velocity value by equidistant inverse ratio weighted interpolation, the search strategy in Interpolation Process is according to grid rope
The horizontal and vertical search of row is introduced, and is progressively extrapolated according to default control points, it is low to avoid the efficiency that global search is brought
The problem of lower.
Further, in the step 11, smoothing processing algorithm uses square neighborhood rolling average algorithm to save search
Time.
Further, in the step 12, smoothing processing algorithm uses square neighborhood rolling average algorithm to save search
Time.
Further, in the step 12, the file for indexing output is binary file.
Compared with prior art, the beneficial effects of the invention are as follows：For the unification of observation system parameter demonstration in the present invention
Layer model method for building up, convenient to carry out, reliable results, the whole needed for the model integrated built observation system parameter demonstration
Information compensates for the deficiency that the conventional parameter with " point " generation " face " is proved, can be calculated, adapted to for entire target zone
High precision seismic exploration developing direction.
Description of the drawings
Fig. 1 is one work area target zone depth model plan view of Shengli Oil Field Dongying City.
Fig. 2 is the inclination layer model plan view in one work area of Shengli Oil Field Dongying City.
Fig. 3 is the layer reflection interval model plan view in one work area of Shengli Oil Field Dongying City.
Fig. 4 is the layer root mean sequare velocity model plan view in one work area of Shengli Oil Field Dongying City.
Fig. 5 is one work area target zone depth model plan view of Shengli Oil Field Xinjiang exploratory area.
Fig. 6 is the inclination layer model plan view in one work area of Shengli Oil Field Xinjiang exploratory area.
Fig. 7 is the layer reflection interval model plan view in one work area of Shengli Oil Field Xinjiang exploratory area.
Fig. 8 is the layer root mean sequare velocity model plan view in one work area of Shengli Oil Field Xinjiang exploratory area.
Specific embodiment
For enable the present invention above and other objects, features and advantages be clearer and more comprehensible, it is cited below particularly go out two preferably
Embodiment, and coordinate institute's accompanying drawings, it is described in detail below.
Embodiment 1.From one three-dimensional seismic acquisition observation system design object of sinopec Shengli Oil Field, work area position
In in Shandong Province's Dongying City, specific embodiment is：
(1) determine work area scope, obtain the layer position instrument of interpretation of target zone and earthquake root mean sequare velocity file；
(2) required according to observation system Parameter analysis, setting analysis grid length and width is 50m, and net is carried out to purpose of horizon
It formats processing；
(3) according to existing purpose, data are explained in position layer by layer, establish layer bit depth model (attached drawing 1) and inclination angle model is (attached
Fig. 2)；
(4) by the data of earthquake root mean sequare velocity discrete sampling point by " when m- speed " to being converted into " depth-speed "
It is right.
(5) reflection interval of layer position and root mean sequare velocity model are established according to the result of (3) and (4), as a result such as 3 He of attached drawing
Shown in attached drawing 4.
(6) according to the layer model established, it is as follows that result is proved to T1 target zone observation system major parameters：Face element ruler
It is very little to be less than 15m, maximum offset 1041m~1450m, it receives line-spacing and is less than 150m.
Embodiment 2.From one, sinopec Shengli Oil Field Xinjiang exploratory area three-dimensional seismic acquisition observation system design object,
The work area is located in the uighur's culture autonomous region of Xinjiang, and specific embodiment is：
(1) determine work area scope, obtain the layer position instrument of interpretation of target zone and earthquake root mean sequare velocity file；
(2) required according to observation system Parameter analysis, setting analysis grid length and width is 50m, and net is carried out to purpose of horizon
It formats processing；
(3) according to existing purpose, data are explained in position layer by layer, establish layer bit depth model (attached drawing 5) and inclination angle model is (attached
Fig. 6)；
(4) by the data of earthquake root mean sequare velocity discrete sampling point by " when m- speed " to being converted into " depth-speed "
It is right.
(5) reflection interval of layer position and root mean sequare velocity model are established according to the result of (3) and (4), as a result such as 7 He of attached drawing
Shown in attached drawing 8.
(6) according to the layer model established, it is as follows that result is proved to T2 target zone observation system major parameters：Face element ruler
It is very little to be less than 25m, maximum offset 1481m~2150m, it receives line-spacing and is less than 300m.
Claims (12)
1. a kind of layer model method for building up for observation system parameter demonstration, which is characterized in that include the following steps:
Step 1：It determines the stratigraphic distribution scope of target zone, obtains the layer position data of target zone and the earthquake root mean sequare velocity in work area
Data；
Step 2：The grid designed according to exploration deployment sets the sizing grid of demonstration and position is carried out at gridding layer by layer to purpose
Reason；
Step 3：According to existing purpose, position data carry out assignment to the layer position trellis depth with the sampling of layer position layer by layer；
Step 4:The layer position grid of layer position sampling is not had using the calculating of equidistant inverse ratio weight interpolation algorithm according to existing depth value
Depth value；
Step 5：Using depth value as colour code, target zone depth model is drawn, checks the purpose effect of bit depth modeling layer by layer, if
Depth model value then needs to be smoothed depth model there are apparent zigzag；
Step 6：According to the inclination angle of all layers of position grid of each mesh coordinate and depth calculation；
Step 7：Using inclination value as colour code, purpose inclination layer model, inspection layer position inclination angle modeling effect, if inclination angle mould are drawn
Type then needs to be smoothed inclination angle model there are apparent zigzag；
Step 8：The data of earthquake root mean sequare velocity discrete sampling point are right to being converted into " depth-speed " by " when m- speed "；
Step 9：The depth of layer position grid obtains the reflection interval of the grid according to where earthquake root mean sequare velocity discrete sampling point
And root mean sequare velocity；
Step 10：All layers of position grid are traveled through, the grid that no reflection interval and root mean sequare velocity sample are utilized equidistant anti-
Than reflection interval and the root mean sequare velocity that this layer of position grid is calculated in weight interpolation algorithm；
Step 11：Using reflection interval as colour code, target zone reflection interval model, inspection layer position reflection interval modeling effect are drawn
Fruit, if reflection interval model needs to be smoothed reflection interval model there are apparent zigzag；
Step 12：Using root mean sequare velocity as colour code, target zone root mean sequare velocity model is drawn, inspection layer position root mean sequare velocity is built
Fruit is imitated, if root mean sequare velocity model needs to be smoothed root mean sequare velocity model there are apparent zigzag；
Step 13：The depth of each grid of target zone, inclination angle, reflection interval and root mean sequare velocity are exported according to grid index, i.e.,
Complete the foundation that the target zone unifies layer model.
2. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 1, the layer position data of acquisition are derived from work area three dimensional seismic data geologic interpretation achievement, are by space
What is exported after discrete sampling includes sample point coordinate, the text file of layer locating depth angle value.
3. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 1, the earthquake root mean square data of acquisition are derived from the achievement of work area D seismic data processing, be by
Exported after spatial spreading sampling include sample point coordinate, " when m- speed " to text file.
4. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 4, without the grid of layer position data sampling, depth value is obtained by equidistant inverse ratio weighted interpolation,
Search strategy in Interpolation Process is to carry out horizontal and vertical search according to grid index, and according to default control points progressively
The problem of extrapolating, avoiding the inefficiency that global search is brought；Wherein default control point is several according to layer position sampled point density
Degree determines.
5. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 5, smoothing processing algorithm uses square neighborhood rolling average algorithm to save search time.
6. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 6, the scope at inclination angle is defined as [0,90o], it is used when specifically calculating adjacent and triangular in shape with calculating grid
The angle of the plane of the central point composition of three layer position grids of distribution, the plane and horizontal plane is the inclination angle for calculating grid.
7. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 7, smoothing processing algorithm uses square neighborhood rolling average algorithm to save search time.
8. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 8, " when m- speed " to " depth-speed " to conversion be directly with time and root mean sequare velocity phase
It is multiplied to arrive depth, and root mean sequare velocity remains unchanged.
9. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 9, the speed of layer position grid is former and later two sampling points according to its depth in " depth-speed " to respective value
It is worth what linear interpolation calculated.
10. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 10, without reflection interval and the grid of root mean sequare velocity data sampling, reflection interval and root mean sequare velocity
Value is obtained by equidistant inverse ratio weighted interpolation, the search strategy in Interpolation Process be according to grid index carry out laterally and
Longitudinal searching, and according to default the problem of points is controlled progressively to extrapolate, avoid the inefficiency that global search is brought.
11. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 11, smoothing processing algorithm uses square neighborhood rolling average algorithm to save search time.
12. a kind of layer model method for building up for observation system parameter demonstration according to claim 1, feature exist
In：In the step 12, smoothing processing algorithm uses square neighborhood rolling average algorithm to save search time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611240398.1A CN106646608B (en) | 2016-12-28 | 2016-12-28 | A kind of layer model method for building up for observation system parameter demonstration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611240398.1A CN106646608B (en) | 2016-12-28 | 2016-12-28 | A kind of layer model method for building up for observation system parameter demonstration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106646608A CN106646608A (en) | 2017-05-10 |
| CN106646608B true CN106646608B (en) | 2018-05-18 |
Family
ID=58836395
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611240398.1A Active CN106646608B (en) | 2016-12-28 | 2016-12-28 | A kind of layer model method for building up for observation system parameter demonstration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106646608B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110244383B (en) * | 2019-06-27 | 2021-06-22 | 中国石油化工股份有限公司 | Geological lithology comprehensive model establishing method based on near-surface data |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1412575A (en) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | Oil-gas detection method based on multiphase medium theory |
| CN104834009A (en) * | 2014-12-02 | 2015-08-12 | 夏正元 | Non-hyperbolic velocity scanning and dynamic correction method |
| CN105223630A (en) * | 2014-06-26 | 2016-01-06 | 中石化石油工程地球物理有限公司胜利分公司 | Based on the omnibearing observation systematic parameter Demonstration Method of geologic model |
| CN105549082A (en) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | Establishing method and system of three-dimensional geomechanical field of extra-deep carbonate reservoir |
| KR20160095690A (en) * | 2015-02-03 | 2016-08-12 | 한국해양대학교 산학협력단 | Seismic waveform inversion method and apparatus for underground velocity information derived from the single-channel seismic data |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2479172B (en) * | 2010-03-31 | 2016-02-10 | Statoil Petroleum As | Estimating interval velocities |
-
2016
- 2016-12-28 CN CN201611240398.1A patent/CN106646608B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1412575A (en) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | Oil-gas detection method based on multiphase medium theory |
| CN105223630A (en) * | 2014-06-26 | 2016-01-06 | 中石化石油工程地球物理有限公司胜利分公司 | Based on the omnibearing observation systematic parameter Demonstration Method of geologic model |
| CN105549082A (en) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | Establishing method and system of three-dimensional geomechanical field of extra-deep carbonate reservoir |
| CN104834009A (en) * | 2014-12-02 | 2015-08-12 | 夏正元 | Non-hyperbolic velocity scanning and dynamic correction method |
| KR20160095690A (en) * | 2015-02-03 | 2016-08-12 | 한국해양대학교 산학협력단 | Seismic waveform inversion method and apparatus for underground velocity information derived from the single-channel seismic data |
Non-Patent Citations (3)
| Title |
|---|
| 三维速度模型及其应用;吴国明等;《石油地球物理勘探》;19880630;第23卷(第3期);第263-276页 * |
| 叠前偏移技术在准噶尔盆地逆冲断裂带地震资料处理中的应用;蔡春江等;《石油勘探与开发》;20040831;第31卷(第4期);第55-57页 * |
| 层控网格层析速度建模技术在陆上盐丘区的应用;王兆旗等;《天然气地球科学》;20161130;第27卷(第11期);第2070-2076页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106646608A (en) | 2017-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102222365B (en) | Method for reconstructing curved surface of complex space | |
| CN103454686B (en) | Small scale sedimentary facies based on strata slicing carries out the method and system of reservoir prediction | |
| CN102867330B (en) | Region-division-based spatial complex horizon reconstruction method | |
| CN103529474B (en) | Adopt lithology segmentation to realize the method for petrofacies fine description | |
| CN106097445B (en) | A kind of method for drafting of three-dimensional formation curved surface | |
| CN105651676A (en) | Reservoir heterogeneity characterization method under regular development well pattern of horizontal wells | |
| CN109388817A (en) | A kind of Reservoir Fracture three-dimensional modeling method | |
| CN104459768B (en) | One kind is based on visual Three dimensional Targets geologic body method for tracing | |
| CN103901478A (en) | Method for determining deposition characteristics and distribution of reservoirs by combining logging and seismic information | |
| CN102385067B (en) | Drawing method for isoline containing reverse fault | |
| CN109884710B (en) | Micro-logging tomography method aiming at excitation well depth design | |
| CN110097637A (en) | A kind of three-dimensional geological attribute model temporal-spatial interpolating method and system | |
| CN106204312A (en) | Rapid economic evaluation method for ion adsorption type rare earth ore | |
| CN103824329A (en) | Geological exploration three-dimensional visual reserve estimation method | |
| CN110824563A (en) | Reservoir lithology prediction method based on Xgboost algorithm | |
| CN112150582B (en) | Multi-modal data-oriented geological profile approximate expression method | |
| CN103226847A (en) | Meandering river point bar 3D data generating method and device | |
| CN107766978A (en) | The intelligent optimization method of skewed pattern | |
| CN102867332A (en) | Multi-level subdivided mesh surface fitting method based on complicated boundary constraint | |
| CN103376463A (en) | Inversion modeling method based on fault control | |
| CN109143337A (en) | A kind of lithologic sand body reservoir modeling methodologies | |
| CN105022856B (en) | Predict the reservoir modeling methodologies of high camber meandering channel reservoir internal structure | |
| CN106646608B (en) | A kind of layer model method for building up for observation system parameter demonstration | |
| CN102830430B (en) | A kind of horizon velocity modeling method | |
| CN105589096B (en) | A kind of sedimentary facies belt division methods based on D-S evidence theory |
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 |