CN106814392B - The velocity modeling method in three-dimensional secondary closely plast fusion face - Google Patents
The velocity modeling method in three-dimensional secondary closely plast fusion face Download PDFInfo
- Publication number
- CN106814392B CN106814392B CN201510853950.3A CN201510853950A CN106814392B CN 106814392 B CN106814392 B CN 106814392B CN 201510853950 A CN201510853950 A CN 201510853950A CN 106814392 B CN106814392 B CN 106814392B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- fusion
- modeling method
- subzone
- data
- 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
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/282—Application of seismic models, synthetic seismograms
Abstract
The invention discloses the velocity modeling methods in three-dimensional secondary closely plast fusion face, the basic thought of this method combined structure modeling, closing geologic body secondary one by one is regarded as closed subzone, the intersection in search face, determine that subzone is decomposed into the form on one or more Interface of superior top or bottom circle by the level of intersection according to the complexity of intersection, fusion and the calculating of space interpolation are carried out, a closed three-dimensional subzone body is eventually formed.Geologic horizon is considered as the curved surface that can decompose fusion by this thought namely from the viewpoint of iconology, forms closed space one by one.It is not limited to specific physical significance, can easier handle the closing of the special geobodies such as tomography, reversed fault.This is a kind of completely new three-dimensional tectonic model method for building up, effectively compensates for the deficiency of traditional modeling method.
Description
Technical field
The invention belongs to geophysics treatment technologies, are a kind of velocity modeling sides in three-dimensional secondary closely plast fusion face
Method.
Background technology
As exploration engineering constantly extends to exploitation, exploration object is increasingly complicated, and two-dimentional forward simulation has been unable to meet reality
Border is applied, and D integral pin-fin tube simulation is more in line with the design of field layout, more really reflects underground propagation and anti-
Penetrate feature.Therefore the foundation of complex geologic body three-dimensional geological model is very important.The modeling software of mainstream currently on the market,
As the modeling function of Petrel softwares and RMS softwares is powerful and abundant, but modeling procedure is complicated, and it is long to expend the time.
Invention content
It is an object of the present invention in view of the deficiencies of the prior art, provide a kind of three-dimensional secondary closely plast fusion face
Velocity modeling method.
The technical solution adopted by the present invention is as follows.
The velocity modeling method in three-dimensional secondary closely plast fusion face, includes the following steps：
Step 1：It determines structure interpretation scheme, obtains layer data, layer position data and well location data；
Step 2：Loaded floor by floor data format to data progress space interpolation and three dimensional network；
Step 3：Three-dimensional visualization monitors difference and gridding process, generation level simultaneously repair loophole in time；
Step 4：Closing geologic body secondary one by one is regarded as closed subzone, the secondary closing geologic body is
Intersect the closing geologic body formed with layer position, layer position and tomography by layer position, secondary is decomposed closely using subzone blending algorithms
Plastid determines the top bottom interface of subzone that is, according to layer position and layer position, layer position and tomography overlapping relation；
Step 5：Fusion and the calculating of space interpolation are carried out, closed three-dimensional geologic is eventually formed；To all tomographies and
Layer position carries out decomposition fusion, forms closed space one by one；
Step 6：Three-dimensional prolongs the point-by-point speed change assignment of layer and output speed model.
Further, in step 1, it when determining structure interpretation scheme, in conjunction with three dimensional seismic data, is shaken by well and combines determination
The rational structure interpretation scheme in work area.
Further, in step 2, it when space interpolation and three dimensional network are formatted, is carried out according to fault strike and dip direction empty
Interpolation and three dimensional network are formatted.
Further, in step 2, it is constraint item with well location data when carrying out space interpolation and gridding to input data
Part.
Further, in step 4, when decomposing secondary closing geologic body, according to the complexity of intersection, secondary is closed
The top interface decomposition of geologic body is multiple Interface of superior top.
Further, in steps of 5, when fusion forms closed three-dimensional geologic, successively by the Interface of superior top of several intersections
Fusion forms top interface, will finally push up the bottom interface that interface is intersected with it and merges.
Further, in steps of 5, it when carrying out fusion and the calculating of space interpolation, is merged by curved surface and determines subzone,
Using three-dimensional inverse distance-weighting, gram in golden spatial interpolation methods, form closed three-dimensional geologic.
Further, step 6 beginning before, different direction section by the three-dimensional tectonic model preliminarily formed with practically
Shake data, well location data etc. are compared, if repeating step 4 to 5 when the two error is larger until error is met the requirements.
Further, in step 6, it refers to that migrated seismic section data volume is combined to carry out three that three-dimensional, which prolongs the point-by-point speed change assignment of layer,
Dimension is along layer speed change assignment.
There is fault block not when complicated structure band area carries out multi-fault mixing intersection in traditional three-dimensional construction modeling technique
The problems such as closure, tectonic model distortion and deviation actual configuration pattern.The velocity modeling in three-dimensional secondary closely plast fusion face
Method combines the basic thought of construction modeling, and closing geologic body secondary one by one is regarded as closed subzone, searches for face
Intersection, determine intersection level subzone is decomposed by one or more Shang Ding circle according to the complexity of intersection
The form in face or bottom circle carries out fusion and the calculating of space interpolation, eventually forms a closed three-dimensional subzone body.This
Geologic horizon is considered as the curved surface that can decompose fusion by kind thought namely from the viewpoint of iconology, and formation is sealed one by one
The space closed.It is not limited to specific physical significance, can easier handle the envelope of the special geobodies such as tomography, reversed fault
It closes.
This is a kind of completely new three-dimensional tectonic model method for building up, effectively compensates for the deficiency of traditional modeling method.
The three-dimensional modeling method in face is merged based on subzone, has modeling efficiency high, model accuracy is high, can realize complicated geological
The characteristics of body three-dimensional geological model is established, has very important significance.
In short, the three-dimensional modeling method based on three-dimensional secondary closely plast fusion face is more in the progress of complicated structure band area
When tomography mixing intersection, it is not easy to the problems such as fault block is not closed, tectonic model distorts and deviates actual configuration pattern occur, have
Modeling efficiency is high, and model accuracy is high, can easier handle the closing of the special geobodies such as tomography, reversed fault.
Description of the drawings
Fig. 1 is that subzone merges face dimensional Modeling Technology flow chart in embodiment 1.
The load of well location data and three-dimensional visualization design sketch in Fig. 2 embodiments 1.
The design sketch of the 3 D geological horizon of gridding in Fig. 3 embodiments 1.
The design sketch of three-dimensional geological level in Fig. 4 embodiments 1 after leak repairing.
The schematic diagram of the decomposition of subzone and fusion process in Fig. 5 embodiments 1.
Subzone decomposes the effect of visualization figure of the three-dimensional velocity structure after fusion in Fig. 6 embodiments 1.
Specific implementation mode
In the following, the present invention will be further described in conjunction with the accompanying drawings and embodiments.
Embodiment 1.A kind of velocity modeling method in three-dimensional secondary closely plast fusion face, includes the following steps.
Step 1：It determines structure interpretation scheme, obtains layer data, layer position data and well location data；
Step 2：Loaded floor by floor data, process of formatting to data progress space interpolation and three dimensional network, it is simultaneously timely to generate level
Repair loophole；
Step 4：Secondary closing geologic body is decomposed using subzone blending algorithms, i.e., by closing geology secondary one by one
It is stereoscopic to make closed subzone, subzone is decomposed into the form of one or more top interface or bottom interface, searches for it
Intersection, cross surface；
Step 5：Fusion and the calculating of space interpolation are carried out, closed three-dimensional geologic is eventually formed；To all tomographies and
Layer position carries out decomposition fusion, forms closed space one by one；
Step 6：Three-dimensional prolongs the point-by-point speed change assignment of layer and output speed model.
In the following, this modeling method is described with reference to the drawings.
As shown in Figure 1, the velocity modeling method in three-dimensional secondary closely plast fusion face, includes the following steps.
The first step：In conjunction with three dimensional seismic data, is shaken by well and combine the rational structure interpretation scheme in work area that determines, acquisition disconnected
Layer data, layer position data and well location data, well location data are as shown in Figure 2.
Second step：Loaded floor by floor data carry out space interpolation and three dimensional network to data according to fault strike and dip direction
It formats, the three-dimension layer position of gridding is as shown in Figure 3.
Third walks：The realization process of second step is monitored by three-dimensional visualization, generation level simultaneously repairs loophole, loophole in time
Three-dimensional geological level after repairing is as shown in Figure 4.
4th step：Subzone algorithms decompose geologic body, determine the tomography of each subzone, the quantity of layer position cross surface, melt
Close the subzone bodies for forming three dimensional closure.As shown in figure 5, first, h4 is merged with f2, it is assumed that merged again with f3 for layer position 1,1
To the top interface of subzone, i.e. white portion in Fig. 5 (c) finally merges to obtain a subzone geologic body with f4.Figure
(a) subzone schematic diagrames are independent, the title that letter in (b) is layer position and tomography is schemed, scheme what number representative in (d) was formed
Subzone serial numbers.
5th step：Different direction section compares three-dimensional tectonic model and actual seismic data, and the two goodness of fit is low,
The 4th step is repeated, migrated seismic section data volume is otherwise combined to carry out three-dimensional along layer speed change assignment, output construction model such as Fig. 6 institutes
Show.
Using the layer position data of 3-D interpretation achievement, layer data and known well location data, loaded floor by floor and with well location
Data are that constraints carries out gridding and space interpolation to input data, and three-dimensional visualization monitors interpolation form and gridding
Journey, and loophole is repaired, it is directed to closed area using the fusion face algorithm of subzone, the top bottom interface of geologic body is decomposed, determines phase
The quantity of cross surface, fusion form closed subzone bodies.According to the step for, all tomographies, layer position are decomposed, merged,
Establish subzone bodies.It finally combines migrated seismic section data volume to carry out three-dimensional and prolongs the point-by-point speed change assignment of layer, export 3D solid
Geological model lays the foundation for subsequent analysis seismic wave filed character, forward simulation and imaging method research.
The basic thought that modeling method combined structure described in the present embodiment models regards closing geologic body secondary one by one
Make closed subzone, search for the intersection in face, determine the level of intersection, according to the complexity of intersection, by subzone points
Solution is the form on one or more Interface of superior top or bottom circle, carries out fusion and the calculating of space interpolation, eventually forms one
Closed three-dimensional subzone bodies.Geologic horizon is considered as to decompose and melt by this thought namely from the viewpoint of iconology
The curved surface of conjunction forms closed space one by one.It is not limited to specific physical significance, can easier handle tomography, inverse
The closing of the special geobodies such as tomography.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirits and principles of the present invention, made by
Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.
Claims (9)
1. the velocity modeling method in three-dimensional secondary closely plast fusion face, which is characterized in that include the following steps：
Step 1：It determines structure interpretation scheme, obtains layer data, layer position data and well location data；
Step 2：Loaded floor by floor data format to data progress space interpolation and three dimensional network；
Step 3：Three-dimensional visualization monitors interpolation and gridding process, generation level simultaneously repair loophole in time；
Step 4：Closing geologic body secondary one by one is regarded as closed subzone, the secondary closing geologic body is by layer
Position intersects the closing geologic body of composition with layer position, layer position with tomography, and secondary closing geology is decomposed using subzone blending algorithms
Body determines the top bottom interface of subzone that is, according to layer position and layer position, layer position and tomography overlapping relation；
Step 5：Fusion and the calculating of space interpolation are carried out, closed three-dimensional geologic is eventually formed；To all tomographies and layer position
Decomposition fusion is carried out, closed space one by one is formed；
Step 6：Three-dimensional is along the point-by-point speed change assignment of layer and output speed model.
2. the velocity modeling method in three-dimensional secondary closely plast fusion face as described in claim 1, it is characterised in that：In step
In rapid 1, when determining structure interpretation scheme, in conjunction with three dimensional seismic data, is shaken by well and combine the determining rational structure interpretation in work area
Scheme.
3. the velocity modeling method in three-dimensional secondary closely plast fusion face as described in claim 1, it is characterised in that：In step
In rapid 2, when space interpolation and three dimensional network are formatted, carries out space interpolation according to fault strike and dip direction and three dimensional network is formatted.
4. the velocity modeling method in three-dimensional secondary closely plast fusion face as described in claim 1, it is characterised in that：In step
In rapid 2, when carrying out space interpolation and gridding to input data, using well location data as constraints.
5. the velocity modeling method in three-dimensional secondary closely plast fusion face as described in claim 1, it is characterised in that：In step
In rapid 4, when decomposing secondary closing geologic body, according to the complexity of intersection, by the top interface decomposition of secondary closing geologic body
For multiple Interface of superior top.
6. the velocity modeling method in three-dimensional secondary closely plast fusion face as claimed in claim 5, it is characterised in that：In step
In rapid 5, when fusion forms closed three-dimensional geologic, the Interface of superior top of several intersections is merged to formation top interface successively, finally
The bottom interface fusion intersected with it at interface will be pushed up.
7. the velocity modeling method in three-dimensional secondary closely plast fusion face as described in claim 1, it is characterised in that：In step
In rapid 5, carry out fusion and space interpolation calculating when, by curved surface merge determine subzone, using three-dimensional inverse distance-weighting,
Golden spatial interpolation methods, form closed three-dimensional geologic in gram.
8. the velocity modeling method in three-dimensional secondary closely plast fusion face as described in claim 1, it is characterised in that：In step
Before rapid 6 start, the three-dimensional tectonic model preliminarily formed and actual seismic data, well location data are carried out in different direction section
Comparison, if repeating step 4 to 5 when the two error is larger until error is met the requirements.
9. the velocity modeling method in three-dimensional secondary closely plast fusion face as described in claim 1, it is characterised in that：In step
In rapid 6, three-dimensional refers to that migrated seismic section data volume is combined to carry out three-dimensional along layer speed change assignment along the point-by-point speed change assignment of layer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510853950.3A CN106814392B (en) | 2015-11-27 | 2015-11-27 | The velocity modeling method in three-dimensional secondary closely plast fusion face |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510853950.3A CN106814392B (en) | 2015-11-27 | 2015-11-27 | The velocity modeling method in three-dimensional secondary closely plast fusion face |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106814392A CN106814392A (en) | 2017-06-09 |
| CN106814392B true CN106814392B (en) | 2018-10-30 |
Family
ID=59156157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510853950.3A Active CN106814392B (en) | 2015-11-27 | 2015-11-27 | The velocity modeling method in three-dimensional secondary closely plast fusion face |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106814392B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110941029B (en) * | 2018-09-25 | 2022-08-16 | 中国石油化工股份有限公司 | Speed modeling method related to geological capping |
| CN109658505B (en) * | 2018-11-05 | 2023-05-09 | 武汉中地数码科技有限公司 | Three-dimensional space conformation method and system with fault plane |
| CN113379907A (en) * | 2020-03-10 | 2021-09-10 | 中国石油天然气股份有限公司 | Method and device for constructing fault block geological model |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7330791B2 (en) * | 2002-10-18 | 2008-02-12 | Exxonmobil Upstream Research Co. | Method for rapid fault interpretation of fault surfaces generated to fit three-dimensional seismic discontinuity data |
| CN102054292A (en) * | 2009-10-30 | 2011-05-11 | 吴立新 | TIN coupling-based overground/underground seamless integrated three-dimensional space modeling method |
| RU2451951C2 (en) * | 2010-08-09 | 2012-05-27 | ОБЩЕСТВО С ОГРАНИЧЕННОЙ ОТВЕТСТВЕННОСТЬЮ "Антел-нефть" | Method of searching for hydrocarbon deposits confined to fissured-cavernous collectors |
| CN102819039A (en) * | 2012-08-22 | 2012-12-12 | 电子科技大学 | Interval velocity model building method under complicated geological conditions |
| CN103646423A (en) * | 2013-12-24 | 2014-03-19 | 中国科学院地质与地球物理研究所 | Three-dimensional geological modeling method and device |
| 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 |
| CN104574511A (en) * | 2014-12-17 | 2015-04-29 | 中国地质大学（武汉） | Rapidly progressive three-dimensional geologic modeling method |
-
2015
- 2015-11-27 CN CN201510853950.3A patent/CN106814392B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7330791B2 (en) * | 2002-10-18 | 2008-02-12 | Exxonmobil Upstream Research Co. | Method for rapid fault interpretation of fault surfaces generated to fit three-dimensional seismic discontinuity data |
| CN102054292A (en) * | 2009-10-30 | 2011-05-11 | 吴立新 | TIN coupling-based overground/underground seamless integrated three-dimensional space modeling method |
| RU2451951C2 (en) * | 2010-08-09 | 2012-05-27 | ОБЩЕСТВО С ОГРАНИЧЕННОЙ ОТВЕТСТВЕННОСТЬЮ "Антел-нефть" | Method of searching for hydrocarbon deposits confined to fissured-cavernous collectors |
| CN102819039A (en) * | 2012-08-22 | 2012-12-12 | 电子科技大学 | Interval velocity model building method under complicated geological conditions |
| AU2014227491A1 (en) * | 2013-09-19 | 2015-04-02 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
| CN103646423A (en) * | 2013-12-24 | 2014-03-19 | 中国科学院地质与地球物理研究所 | Three-dimensional geological modeling method and device |
| CN104330825A (en) * | 2014-07-31 | 2015-02-04 | 中国石油化工股份有限公司 | Fault interpretation processing method by slicing along structural trend surface |
| CN104574511A (en) * | 2014-12-17 | 2015-04-29 | 中国地质大学（武汉） | Rapidly progressive three-dimensional geologic modeling method |
Non-Patent Citations (5)
| Title |
|---|
| Three-dimensional modeling and application in geological exploration engineering;Jianya Gong et al.;《Computers & Geosciences》;20041231;第30卷;第391-404页 * |
| Topological relations embodied in a generalized tri-prism (GTP) model for a 3D geoscience modeling system;Wu Lixin;《Computers & Geoscience》;20041231;第30卷;第405-418页 * |
| 基于Kriging的三维地质体建模方法研究;朱燕玲等;《计算机时代》;20101231(第1期);第26-27、30页 * |
| 建立盆地三维构造-地层格架的插值方法;田宜平等;《地球科学-中国地质大学学报》;20000331;第25卷(第2期);第191-194页 * |
| 约束Delaunay四面体剖分在三维地质建模中的应用;余淑娟等;《地理与地理信息科学》;20130131;第29卷(第1期);第41-44页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106814392A (en) | 2017-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10838108B2 (en) | Method and system for interpolating discontinuous functions in a subsurface model | |
| CN104574511B (en) | A kind of quick progressive three-dimensional geological modeling method | |
| US6907392B2 (en) | Method of generating a hybrid grid allowing modelling of a heterogeneous formation crossed by one or more wells | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| US20180348401A1 (en) | Constructing Structural Models of the Subsurface | |
| Sprague et al. | Interpretive tools for 3-D structural geological modelling part II: Surface design from sparse spatial data | |
| CN102609982B (en) | Topology discovery method of space geological data based on unstructured mode | |
| CN101906965A (en) | Technology for identifying and analyzing stratums in VR (Virtual Reality) geological environment | |
| CN106097445B (en) | A kind of method for drafting of three-dimensional formation curved surface | |
| CN106023313A (en) | 3D geologic digital reconnaissance report generating method | |
| CN106338761B (en) | Three-dimensional inserted body modeling method | |
| CN102867332B (en) | Based on the multistage subdivided meshes curved surface fitting method of complex boundary constraint | |
| CN103218850B (en) | One very three-dimensional mining blasting unit model method for building up | |
| CN110244021B (en) | Stratum layering method based on anisotropic interpolation | |
| CN106814392B (en) | The velocity modeling method in three-dimensional secondary closely plast fusion face | |
| CN108898670B (en) | Three-dimensional geological modeling method based on entity and profile | |
| Tonini et al. | 3D geo-mapping based on surface data for preliminary study of underground works: A case study in Val Topina (Central Italy) | |
| CN105719342B (en) | A kind of the three-dimensional modeling method for visualizing and device of ground fissure geologic body | |
| CN110349262A (en) | A kind of tunnel geological forecast three-dimensional modeling method and system based on BIM | |
| CN111859687A (en) | Mixed geological modeling method and system for depicting geological structure of uranium-bearing sand layer | |
| CN106920275A (en) | A kind of complex properties border trivector alternative manner and application system | |
| CN115035258A (en) | Efficient urban three-dimensional geological modeling method based on CAD (computer-aided design) drilling histogram | |
| Wu et al. | An automatic three-dimensional geological engineering modeling method based on tri-prism | |
| Liu et al. | Study on a computing technique suitable for true 3D modeling of complex geologic bodies | |
| Guo et al. | An automatic method for generating curvilinear geological section with stratal pinch-out |
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 |