CN106338761B - Three-dimensional inserted body modeling method - Google Patents
Three-dimensional inserted body modeling method Download PDFInfo
- Publication number
- CN106338761B CN106338761B CN201510394355.8A CN201510394355A CN106338761B CN 106338761 B CN106338761 B CN 106338761B CN 201510394355 A CN201510394355 A CN 201510394355A CN 106338761 B CN106338761 B CN 106338761B
- Authority
- CN
- China
- Prior art keywords
- model
- dimensional
- tectonic
- layer position
- modeling method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of three-dimensional inserted body modeling method, and this method includes：Step 1, structure interpretation scheme is determined, obtains tomography and layer position data；Step 2, according to section or construction trend, gridding direction and boundary condition are determined；Step 3, fault block model is generated, and judges whether the fault block model of generation coincide with structure interpretation scheme；Step 4, when fault block model and structure interpretation scheme are coincide, according to tomography and layer position contact relation, layer position high density scatterplot generation aspect；Step 5, Tectonic Framework model is established, and judges whether Tectonic Framework model meets with original earthquake imaging section；Step 6, when Tectonic Framework model and original earthquake imaging section meet, complete three-dimensional tectonic model is established.The inserted body modeling method of the three-dimensional overcomes the defects of traditional modeling pattern is present, and can more really reduce subsurface geology feature and structural configuration, more effectively to instruct work of exploration and development to establish geological foundation.
Description
Technical field
The present invention relates to seism processing application field, especially relates to a kind of three-dimensional inserted solid modelling side
Method.
Background technology
With the increasingly raising of exploration and development degree, complicated structure band is increasingly becoming the key areas of oil-gas exploration.Construction
Modeling technique has very important as cognition subsurface geology situation and the important means of construction feature in oil-gas exploration and development
Effect.Three-dimensional inserted body modeling method is that solve a kind of effective ways that complicated structure regional model is established.Tradition three
Dimension construction modeling technique is to establish closed model using the constraint connection between tomography and tomography, tomography and aspect, and its defect exists
Fault block is not closed when fault block modeling is carried out using the tomography that tendency is different and mutually cuts, and ridge fault block breakthrough stratum,
Stratigraphic model can not be established correctly when angular unconformity contacts between stratum, cause tectonic model to distort.Traditional three-dimensional construction
Modeling technique carried out in complicated structure band area multi-fault mixing it is intersecting when, occur fault block do not close, tectonic model distortion and partially
The problems such as from actual configuration pattern.For this we have invented a kind of new three-dimensional inserted body modeling method, solves the above
Technical problem.
The content of the invention
It is an object of the invention to provide one kind to utilize building for " high order level faults control screen work, the covering of foreword level fault mosaic "
Mould mode, overcome the three-dimensional inserted body modeling method for the defects of traditional modeling pattern is present.
The purpose of the present invention can be achieved by the following technical measures：Three-dimensional inserted body modeling method, three-dimensional edge
Embedded body modeling method includes：Step 1, structure interpretation scheme is determined, obtains tomography and layer position data；Step 2, according to section
Or construction, determine gridding direction and boundary condition；Step 3, generate fault block model, and judge generation fault block model whether with
Structure interpretation scheme is coincide；Step 4, when fault block model and structure interpretation scheme are coincide, foundation tomography and layer position contact relation,
Layer position high density scatterplot generation aspect；Step 5, establish Tectonic Framework model, and judge Tectonic Framework model whether with primitively
Shake imaging section meets；Step 6, when Tectonic Framework model and original earthquake imaging section meet, complete three-dimensional structure is established
Modeling type.
The purpose of the present invention can be also achieved by the following technical measures：
In step 1, analyzed area tectonic movement, with reference to three dimensional seismic data, combination is shaken by well and determines that work area is reasonable
Structure interpretation scheme, output construction interpretation results simultaneously carry out time and depth transfer.
In step 2, tomography net boundary is determined, intersecting section is connected, determines high order level fault boundary, walked according to tomography
Three-dimensional grid direction and grid trend are determined to dip direction.
In step 3, using Ke Lijin, inverse distance-weighting algorithm and digitlization filtering method, suitable grid is set to join
It is several that tomography high density scatterplot data are carried out with gridding generation fault block model.
In step 3, when fault block model and structure interpretation scheme misfit, repeat step 2.
In step 4, it is determined that layer position net boundary, according to the contact relation between tomography and layer position, set distance value and
Mesh parameter, layer position high density scatterplot data are carried out with gridding generation aspect.
In steps of 5, the top bottom surface on different the sedimentary age stratum is determined using index bed aspect, establishes Tectonic Framework mould
Type, index bed refer to that distribution is stable, have one layer or one group of the obvious specific characteristic rock stratum that can be followed the trail of or contrast as stratum.
In steps of 5, when Tectonic Framework model and original earthquake imaging section are not inconsistent, repeat step 4.
In step 6, the fault block model and special geology body Model of low order fault composition are established, is mounted to construction bone
In frame model, complete three-dimensional tectonic model is established.
The inserted body modeling method of the three-dimensional also includes, and after step 6, different azimuth section is by three-dimensional tectonic model
Contrasted with actual seismic data, when both goodnesses of fit are low, repeat step 6；
When both goodnesses of fit are high, carried out with reference to migrated seismic section data volume it is three-dimensional along the horizontal point-by-point speed change assignment in layer position,
Export three-dimensional modeling data.
Three-dimensional inserted body modeling method in the present invention, on the basis of effective structural analysis, utilizes 3-D seismics
Data and well shake joint determine rational structure interpretation scheme, and interpretation results are carried out into time and depth transfer；Utilize Ke Lijin, anti-distance
Weighting algorithm and digitlization filtering method carry out grid to fault interpretation achievement and establish section, and intersecting high order level section is carried out
Connection, plane fault block border is determined, set fault block grid direction and grid trend, grid to establish full work area three-dimension disclocation model；
The contact relation of tomography and aspect is determined, net is carried out to layer position interpretation results using inverse distance-weighting algorithm and digitlization filter method
Lattice establish aspect；Three-dimensional Tectonic Framework model is obtained after defining stratigraphic model and grid using index bed；Establish small sequence tomography
The tectonic model of control, inlay to whole district's three-dimensional tectonic framework model, obtain complete construction model；Finally cutd open with reference to seismic migration
Face data volume carries out point-by-point speed change assignment in three-dimension layer, exports 3D solid geological model, be subsequent analysis seismic wave filed character,
Forward simulation and imaging method research lay the foundation.
Brief description of the drawings
Fig. 1 is the flow chart of a specific embodiment of the three-dimensional inserted body modeling method of the present invention；
Fig. 2 is the schematic diagram of three-dimensional Tectonic Framework model in one embodiment of the invention；
Fig. 3 is the schematic diagram of Complete three-dimensional tectonic model in one embodiment of the invention；
Fig. 4 is the schematic diagram of three-dimensional tectonic model section in one embodiment of the invention.
Embodiment
For enable the present invention above and other objects, features and advantages become apparent, it is cited below particularly go out preferable implementation
Example, and coordinate institute's accompanying drawings, it is described in detail below.
As shown in figure 1, Fig. 1 is the flow chart of the three-dimensional inserted body modeling method of the present invention.
In step 101, structure interpretation scheme is determined, obtains tomography and layer position data.Analyzed area tectonic movement, with reference to three
Seismic data is tieed up, the rational structure interpretation scheme in work area that determines that combines is shaken by well, output construction interpretation results are simultaneously deep when carrying out
Conversion.Flow enters step 102.
In step 102, according to section or construction trend, gridding direction and boundary condition are determined.Determine tomography Grid Edge
Boundary, connects intersecting section, determines high order level fault boundary, and three-dimensional grid direction and net are determined according to fault strike and dip direction
Lattice trend.Flow enters step 103.
In step 103, using Ke Lijin, inverse distance-weighting algorithm and digitlization filtering method, suitable grid is set to join
It is several that tomography high density scatterplot data are carried out with gridding generation fault block model.Fault block model is misfitted with structure interpretation scheme, weight
Multiple step 102.If coincideing, flow enters step 104.
In step 104, according to tomography and layer position contact relation, layer position high density scatterplot generates aspect.It is determined that layer position grid
Border, according to the contact relation between tomography and layer position, correct distance value and suitable mesh parameter are set, it is highly dense to layer position
Spend scatterplot data and carry out gridding generation aspect.Flow enters step 105.
In step 105, the top bottom surface on different the sedimentary age stratum is determined using index bed aspect, establishes Tectonic Framework model
Such as Fig. 2.Index bed refers to that distribution is stable, has one layer or one group of the obvious specific characteristic rock that can be followed the trail of or contrast as stratum
Layer.Tectonic Framework model and original earthquake imaging section are contrasted, both are not inconsistent repeat step 104.If meeting, flow
Enter step 106.
In step 106, the fault block model and special geology body Model of low order fault composition are established, is mounted to construction bone
It is final to establish complete three-dimensional tectonic model such as Fig. 3 in frame model.Flow enters step 107.
In step 107, different azimuth section is contrasted three-dimensional tectonic model and actual seismic data, both goodnesses of fit
Low, repeat step 106, otherwise, flow enters step 108.
In step 108, when both goodnesses of fit are high, it is laterally point-by-point along layer position to carry out three-dimensional with reference to migrated seismic section data volume
Speed change assignment, output three-dimensional modeling data such as Fig. 4.Flow terminates.
Three-dimensional inserted body modeling method in the present invention, combines geological knowledge, seismic imaging section and well location
Multi-source information, gridding interpolation is carried out using Ke Lijin, inverse distance-weighting algorithm and digitlization filtering method, with reference to special geobody
It is a kind of brand-new three-dimensional tectonic model method for building up with the inserted processing of low order fault, Complicated Fault Zone is three-dimensional to be constructed
Modeling technique, for fracture complicated band area by the way of small facies tract is inlayed, it is ensured that the closure on fault block and stratum and
The integrality of construction, subsurface geology feature and structural configuration can be more really reduced, effectively compensate for traditional modeling
Deficiency, for more really cognition subsurface geology situation, and then can be three-dimensional construction forward modeling and oilfield prospecting developing subsequently
Model media is provided.
Claims (7)
1. a kind of three-dimensional inserted body modeling method, it is characterised in that this method includes：
Step 1, structure interpretation scheme is determined, obtains tomography and layer position data；
Step 2, according to section or construction, gridding direction and boundary condition are determined；
Step 3, fault block model is generated, and judges whether the fault block model of generation coincide with structure interpretation scheme；
Step 4, when fault block model and structure interpretation scheme are coincide, according to tomography and layer position contact relation, layer position high density dissipates
Point generation aspect；
Step 5, Tectonic Framework model is established, and judges whether Tectonic Framework model meets with original earthquake imaging section；
Step 6, when Tectonic Framework model and original earthquake imaging section meet, complete three-dimensional tectonic model is established；
In step 3, when fault block model and structure interpretation scheme misfit, repeat step 2；
In steps of 5, when Tectonic Framework model and original earthquake imaging section are not inconsistent, repeat step 4；
After step 6, different azimuth section is contrasted three-dimensional tectonic model and actual seismic data, when both goodnesses of fit
When low, repeat step 6；When both goodnesses of fit are high, carry out three-dimensional with reference to migrated seismic section data volume becomes along layer position is laterally point-by-point
Fast assignment, export three-dimensional modeling data.
A kind of 2. three-dimensional inserted body modeling method according to claim 1, it is characterised in that in step 1, analysis
Regional structure campaign, with reference to three dimensional seismic data, shaken by well and combine the determination rational structure interpretation scheme in work area, output construction
Interpretation results simultaneously carry out time and depth transfer.
A kind of 3. three-dimensional inserted body modeling method according to claim 1, it is characterised in that in step 2, it is determined that
Tomography net boundary, connects intersecting section, determines high order level fault boundary, and three dimensional network is determined according to fault strike and dip direction
Lattice direction and grid trend.
4. a kind of three-dimensional inserted body modeling method according to claim 1, it is characterised in that in step 3, utilize
Ke Lijin, inverse distance-weighting algorithm and digitlization filtering method, set suitable mesh parameter to tomography high density scatterplot data
Carry out gridding generation fault block model.
A kind of 5. three-dimensional inserted body modeling method according to claim 1, it is characterised in that in step 4, it is determined that
Layer position net boundary, according to the contact relation between tomography and layer position, sets distance value and mesh parameter, layer position high density is dissipated
Point data carries out gridding generation aspect.
6. a kind of three-dimensional inserted body modeling method according to claim 1, it is characterised in that in steps of 5, utilize
Index bed aspect determines the top bottom surface on different the sedimentary age stratum, establishes Tectonic Framework model, and index bed refers to that distribution is stable, has
There is the one layer or one group rock stratum that can be followed the trail of or contrast as stratum of obvious specific characteristic.
7. a kind of three-dimensional inserted body modeling method according to claim 1, it is characterised in that in step 6, establish
The fault block model and special geology body Model of low order fault composition, are mounted in Tectonic Framework model, establish complete three
Tie up tectonic model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510394355.8A CN106338761B (en) | 2015-07-07 | 2015-07-07 | Three-dimensional inserted body modeling method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510394355.8A CN106338761B (en) | 2015-07-07 | 2015-07-07 | Three-dimensional inserted body modeling method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106338761A CN106338761A (en) | 2017-01-18 |
| CN106338761B true CN106338761B (en) | 2018-03-02 |
Family
ID=57827115
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510394355.8A Active CN106338761B (en) | 2015-07-07 | 2015-07-07 | Three-dimensional inserted body modeling method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106338761B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107703542B (en) * | 2017-08-31 | 2020-04-10 | 中国石油天然气集团公司 | Method and device for determining seismic horizon of reverse fault area |
| CN107728216B (en) * | 2017-09-29 | 2019-08-30 | 中国石油化工股份有限公司 | A kind of D integral pin-fin tube method for establishing model based on flume experiment |
| CN109655898A (en) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | The attribute modeling method and system of the embedded filling of complex geologic body |
| CN109035411B (en) * | 2018-07-30 | 2021-07-20 | 中国石油化工股份有限公司 | Three-dimensional stratum corner grid voxel modeling method based on forward modeling |
| CN112562075A (en) * | 2019-09-26 | 2021-03-26 | 中国石油天然气集团有限公司 | Three-dimensional geological modeling method and device for multi-slip-layer structure |
| CN113379907A (en) * | 2020-03-10 | 2021-09-10 | 中国石油天然气股份有限公司 | Method and device for constructing fault block geological model |
| CN113589367B (en) * | 2020-04-30 | 2023-04-28 | 中国石油化工股份有限公司 | Method for correcting nearby structure trend based on area conservation and large fracture |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303414A (en) * | 2008-05-22 | 2008-11-12 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
| RU2451951C2 (en) * | 2010-08-09 | 2012-05-27 | ОБЩЕСТВО С ОГРАНИЧЕННОЙ ОТВЕТСТВЕННОСТЬЮ "Антел-нефть" | Method of searching for hydrocarbon deposits confined to fissured-cavernous collectors |
| CN104330825A (en) * | 2014-07-31 | 2015-02-04 | 中国石油化工股份有限公司 | Fault interpretation processing method by slicing along structural trend surface |
| AU2014227491A1 (en) * | 2013-09-19 | 2015-04-02 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2003275021A1 (en) * | 2002-10-18 | 2004-05-13 | Exxonmobil Upstream Research Company | A method for rapid fault interpretation of fault surfaces generated to fit three-dimensional seismic discontinuity data |
-
2015
- 2015-07-07 CN CN201510394355.8A patent/CN106338761B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303414A (en) * | 2008-05-22 | 2008-11-12 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
| RU2451951C2 (en) * | 2010-08-09 | 2012-05-27 | ОБЩЕСТВО С ОГРАНИЧЕННОЙ ОТВЕТСТВЕННОСТЬЮ "Антел-нефть" | Method of searching for hydrocarbon deposits confined to fissured-cavernous collectors |
| AU2014227491A1 (en) * | 2013-09-19 | 2015-04-02 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
| CN104330825A (en) * | 2014-07-31 | 2015-02-04 | 中国石油化工股份有限公司 | Fault interpretation processing method by slicing along structural trend surface |
Non-Patent Citations (2)
| Title |
|---|
| 卫2块油气藏三维构造建模研究;刘龙龙 等;《中国化工贸易》;20150423(第11期);第1、4页 * |
| 建立盆地三维构造-地层格架的插值方法;田益平 等;《地球科学》;20000331;第25卷(第2期);第191-194页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106338761A (en) | 2017-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106338761B (en) | Three-dimensional inserted body modeling method | |
| CN107817535B (en) | The determination method and apparatus of short lap | |
| CA2579138C (en) | Geologic models of subsurface sedimentary volumes | |
| Molinari et al. | Development and testing of a 3D seismic velocity model of the Po Plain sedimentary basin, Italy | |
| Purvis et al. | Complex reservoir geometries in a deep water clastic sequence, Gryphon Field, UKCS: injection structures, geological modelling and reservoir simulation | |
| NO347906B1 (en) | Stratigraphic function | |
| CN105467464B (en) | The effective connected component of Analyze of Nearshore Subaqueous Fans glutenite is divided and control methods | |
| Ampomah et al. | An integrated approach for characterizing a sandstone reservoir in the Anadarko Basin | |
| Tang et al. | A three dimensional high-resolution reservoir model of Napo Formation in Oriente Basin, Ecuador, integrating sediment dynamic simulation and geostatistics | |
| Zhang et al. | New Insights Into the Rift‐To‐Drift Process of the Northern South China Sea Margin Constrained by a Three‐Dimensional Wide‐Angle Seismic Velocity Model | |
| CN104459798A (en) | Velocity modeling method based on RTM imaging | |
| CN106814392B (en) | The velocity modeling method in three-dimensional secondary closely plast fusion face | |
| CN111859687A (en) | Mixed geological modeling method and system for depicting geological structure of uranium-bearing sand layer | |
| Laribi et al. | Integrated fractured reservoir characterization and simulation: Application to Sidi el Kilani Field, Tunisia | |
| Norden et al. | From pilot site knowledge via integrated reservoir characterization to utilization perspectives of a deep geothermal reservoir: 3D geological model at the research platform Groß Schönebeck in the Northeast German Basin | |
| Abdelnabi et al. | Integrated geologic, geophysical, and petrophysical data to construct full field geologic model of Cambrian-Ordovician and Upper Cretaceous reservoir formations, Central Western Sirte Basin, Libya | |
| Lozano et al. | Integrated Geocellular Static Model for Geomechanical and Dynamic Simulations in the Vaca Muerta Formation | |
| Wang et al. | Modeling dolomitized carbonate-ramp reservoirs: A case study of the Seminole San Andres unit—Part II, Seismic modeling, reservoir geostatistics, and reservoir simulation | |
| Obafemi et al. | 3D facies and reservoir property prediction of deepwater turbidite sands; case study of an offshore Niger delta field | |
| Bu et al. | Geological modeling of combined incised channel system in deep-water gravity flow sedimentation: the case of OML130 in Niger delta basin, West Africa | |
| Zeynalov et al. | Petro physical Modeling of the Tertiary Reservoirs in the Small Giant Bai Hassan Oil Field, Northern Iraq | |
| Prianda et al. | Enhancing Geostatistical Simulation using Additional Constrained Trend from Seismic-Derived Physical Parametere for a Quantitative and Geologically Reasonable Sand Facies Distribution Model at Fluvial Environment | |
| Jiang et al. | Paleogene sedimentary system and sedimentary dynamics of windfield-source-basin system in the Dongying Sag | |
| Ali et al. | Prospectivity and Resource Potential of the Malaysia-Thailand Joint Development Area JDA: Results from an Integrated Play Based Exploration Project | |
| MacMahon | Understanding the tectonic evolution of the West Orphan Basin, offshore Canada, and the conjugate Rockall Basin, offshore Ireland, using a seismic megatransect |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |