CN104535391A - Physical geography data processing method based on layered geography models - Google Patents
Physical geography data processing method based on layered geography models Download PDFInfo
- Publication number
- CN104535391A CN104535391A CN201410844376.0A CN201410844376A CN104535391A CN 104535391 A CN104535391 A CN 104535391A CN 201410844376 A CN201410844376 A CN 201410844376A CN 104535391 A CN104535391 A CN 104535391A
- Authority
- CN
- China
- Prior art keywords
- model
- layer
- geophysical
- layered
- processing 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method used for establishing underground medium models through physical geography calculation. On the basis of the hypothesis of underground medium layered models, with bottom boundaries of all stratums being bottoms and the surface layer being the top, a multi-layer combined model with coincidence regions is established, each layer of the model can be dissected freely through a vertical polygonal prism to establish a space model, and through performing attribute simulation on physical properties of the multi-layer model, real layered stratum properties are established. Through the method, underground layered media can be sufficiently and fine described, and the calculation efficiency of the physical geography effect is improved. The method can be used for the application field of forward-inversion calculation of physical geography gravity and magnetic prospecting and others.
Description
Technical field
The present invention relates to geophysical probing technique, specifically a kind of Geophysical Data Processing method based on layered geology model, the energy complex layered model of fine description also improves geophysical computing efficiency.
Background technology
Geophysics modeling is the basis of carrying out geophysics forward modelling.Modeling method determines the mode of underground Media Description and has a direct impact geophysical computing.Underground medium has Lamellar character usually, and therefore in geophysical research, carrying out modeling based on layered medium is comparatively general method.When possessing multilayered medium, usually by describing underground medium character between layers to the model facetization of inter-level dielectric.And because different layers medium complexity is different in actual computation, not identical with top layer partition patterns at the bottom of model meshes subdivision upper interlayer dielectric, this can cause complicated geophysical computing usually, affect geophysics and just drilling and the efficiency of Inversion Calculation.For gravitational inversion, effectively the simplest with the mode calculated gravity anomaly of prism subdivision, for adjacent inter-level dielectric, intermediate interface be last layer medium end circle and as the top circle of lower one deck, if two-layer partition patterns difference certainly will cause repeatedly subdivision and the calculating of middle layer model, this have impact on the fine description to middle layer character, is also unfavorable for improving counting yield.Therefore, avoid the fine description that then can realize layer position to the repeatedly subdivision of shared layer position, and improve counting yield, geophysics is just being drilled and Inversion Calculation etc. significant.
Summary of the invention
Based on problems of the prior art and reality need, the invention provides one, based on the Geophysical Data Processing method of layered geology model, its objective is and set up multiple layer combination geophysical model based on layered geology model theory, the fine description of multi-layer underground medium can be realized, and geophysical computing efficiency can be improved, forward modelling and the geophysical exploration research of gravity and magnetic anomalies can be directly used in.
According to the first technical scheme of the present invention, a kind of Geophysical Data Processing method based on layered geology model, it, for building the geophysical model based on layered medium, comprises the following steps:
(1), raw data obtains, input Bottom surfaces of strata distribution range and depth information;
(2), by this layer model carry out subdivision by plane space scope, subdivision trellis depth value can be obtained by input depth information interpolation;
(3), with subdivision grid for xsect, vertically extend to earth's surface, form the triangulated space of stratified model;
(4) this layer model subdivision grid Geophysical Properties parameter, is given;
(5), by (1)-(4) step set up other layer models, set up to all layer models;
(6) interlayer model Geophysical Properties parameter, is calculated.The inter-level dielectric geophysical parameters that boundary at the bottom of adjacent two layers determines is the Geophysical Properties parameter value sum of stratified model below boundary at the bottom of this inter-level dielectric and this layer.
According to the second technical scheme of the present invention, a kind of Geophysical Data Processing method based on layered geology model, it, for building the geophysical model based on layered medium, is characterized in that, this Geophysical Data Processing method comprises the following steps:
The first step, raw data obtain, and obtain the subsurface formations degree of depth according to geophysical exploration (earthquake, drilling well etc.);
Second step, build the geophysical model based on layered medium based on the depth of stratum Data Data that obtains in the first step, concrete methods of realizing comprises:
Underground medium model is divided into N layer from earth's surface to underground, and from earth's surface, following ground floor is labeled as 1 to N to bottommost layer numbering; I-th (i=1,2 ..., N) layer by layer shape model representation be the underground space of i-th layer of bottom boundary to earth's surface, equivalent layer model geophysical parameters be labeled as Pi (i=1,2 ..., N); , N is positive integer;
I-th layer model process of establishing is: according to the i-th Ceng Di circle degree of depth by its planar distribution coordinate (x, y) subdivision is some arbitrary polygon grids, and with subdivision grid for xsect is set up to surperficial vertical prism, and give each subdivision grid corresponding Geophysical Properties value Pi; In i-th Ceng Di circle to the i-th-1 Ceng Di circle scope, medium Geophysical Properties can be characterized by formula
3rd step, above step set up the layered geology model with specific Geophysical Properties (as density, magnetic susceptibility, resistivity, speed etc.), can be used for geophysical anomaly forward modelling.
Wherein, each layer model carries out mesh generation with this layer of bottom boundary form and vertically splits this stratified model space.
Preferably, surface layer is the interface, public top of each stratified model.The shape that each layer model carries out subdivision grid according to bottom boundary can be arbitrary polygon.Each layer model each layer subdivision grid can according to actual conditions self-adaptative adjustment, and each layer mesh generation form does not require identical.Each layer model has corresponding geophysical parameters.
Further, each layer model can have transformable geophysical parameters according to subdivision grid.Geophysical parameters is characterized by density or the magnetic susceptibility of model.Can be added by the Geophysical Properties of all layer models comprising this interlayer model part by the Geophysical Properties of this interlayer model of the determined interlayer model in boundary at the bottom of adjacent layer and obtain.
Boundary at the bottom of every one deck and earth's surface stratified model that observation layer is determined are carried out combination acquisition and integrally descend dielectric distribution model, by the distribution simulate formation medium actual conditions of Geophysical Properties parameter by the modeling method that the present invention adopts.
Present invention, avoiding the subdivision of interlayer model common interface, modeling method has more adaptivity, is suitable for the fine description of layered medium, is conducive to improving counting yield.
Accompanying drawing explanation
Fig. 1 is individual layer geophysical model subdivision cross-sectional view;
Fig. 2 is individual layer geophysical model subdivision longitudinal section;
Fig. 3 is two-layer geophysical model subdivision built-up pattern longitudinal diagram.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only a part of embodiment of the present invention, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.Additionally, described embodiment is only to further elaboration of the present invention, but not limitation of the present invention.
According to the first technical scheme of the present invention, a kind of Geophysical Data Processing method based on layered geology model, it, for building the geophysical model based on layered medium, comprises the following steps:
(1), raw data obtains, input Bottom surfaces of strata distribution range and depth information;
(2), by this layer model carry out subdivision by plane space scope, subdivision trellis depth value can be obtained by input depth information interpolation;
(3), with subdivision grid for xsect, vertically extend to earth's surface, form the triangulated space of stratified model;
(4) this layer model subdivision grid Geophysical Properties parameter, is given;
(5), by (1)-(4) step set up other layer models, set up to all layer models;
(6) interlayer model Geophysical Properties parameter, is calculated.The inter-level dielectric geophysical parameters that boundary at the bottom of adjacent two layers determines is the Geophysical Properties parameter value sum of stratified model below boundary at the bottom of this inter-level dielectric and this layer.
According to the second technical scheme of the present invention, a kind of Geophysical Data Processing method based on layered geology model, it, for building the geophysical model based on layered medium, is characterized in that, this Geophysical Data Processing method comprises the following steps:
The first step, raw data obtain, and obtain the subsurface formations degree of depth according to geophysical exploration (earthquake, drilling well etc.);
Second step, build the geophysical model based on layered medium based on the depth of stratum Data Data that obtains in the first step, concrete methods of realizing comprises:
Underground medium model is divided into N layer from earth's surface to underground, and from earth's surface, following ground floor is labeled as 1 to N to bottommost layer numbering; I-th (i=1,2 ..., N) layer by layer shape model representation be the underground space of i-th layer of bottom boundary to earth's surface, equivalent layer model geophysical parameters be labeled as Pi (i=1,2 ..., N); , N is positive integer;
I-th layer model process of establishing is: according to the i-th Ceng Di circle degree of depth by its planar distribution coordinate (x, y) subdivision is some arbitrary polygon grids, and with subdivision grid for xsect is set up to surperficial vertical prism, and give each subdivision grid corresponding Geophysical Properties value Pi; In i-th Ceng Di circle to the i-th-1 Ceng Di circle scope, medium Geophysical Properties can be characterized by formula
3rd step, above step set up the layered geology model with specific Geophysical Properties (as density, magnetic susceptibility, resistivity, speed etc.), can be used for geophysical anomaly forward modelling.
Wherein, each layer model carries out mesh generation with this layer of bottom boundary form and vertically splits this stratified model space.
Preferably, surface layer is the interface, public top of each stratified model.The shape that each layer model carries out subdivision grid according to bottom boundary can be arbitrary polygon.Each layer model each layer subdivision grid can according to actual conditions self-adaptative adjustment, and each layer mesh generation form does not require identical.Each layer model has corresponding geophysical parameters.
Further, each layer model can have transformable geophysical parameters according to subdivision grid.Geophysical parameters is characterized by density or the magnetic susceptibility of model.Can be added by the Geophysical Properties of all layer models comprising this interlayer model part by the Geophysical Properties of this interlayer model of the determined interlayer model in boundary at the bottom of adjacent layer and obtain.
Below in conjunction with specific embodiment, give to explain explanation further to the present invention.
Embodiment 1: carry out gravity anomaly with two-layer underground medium modeling and be calculated as example enforcement underground medium density modeling.By of the present invention, modeling concrete steps are as follows:
(1) input raw data, in this example, stratum is two-layer, is upper crust and lower crust.The input data bit upper crust bottom boundary degree of depth and the lower crust bottom boundary degree of depth;
(2) lower crust layer layer model is accordingly set up.According to lower crust stratigraphic distribution data, carry out rectangular node subdivision to this layer of distribution range along xy plane, this layer depth of each grid representative is determined (see accompanying drawing 1) by input value interpolation;
(3) be that xsect sets up upright hexahedron subdivision this layer of bottom boundary to spatial surface scope (accompanying drawing 2 for corresponding longitudinal profile) with Rectangular Partition grid step (2) Suo Shi;
(4) density value giving this layer model subdivision grid is 3300kg/m
3;
(5) upper crust layer model is accordingly set up.According to upper crust stratigraphic distribution data, along xy plane, mesh generation is held to this layer of distribution range, this layer depth of each grid representative is determined by input value interpolation, upper crust subdivision network style is identical with step (2), the layer data distribution of mesh generation large I base area changes, in this example, the model meshes size that upper crust model is corresponding is less than corresponding model meshes size (accompanying drawing 1 center section size of mesh opening represents upper crust mesh generation size) of lower crust model;
(6) be that xsect sets up upright hexahedron subdivision this layer of bottom boundary to spatial surface scope (accompanying drawing 3 upper formation model is corresponding longitudinal profile) with rectangular node step (5) Suo Shi; Set up ground floor dielectric model by (1)-(4), this layer model density value is-600kg/m
3
(7) be jointly made up of the geologic model (accompanying drawing 3) with double-deck crustal structure lower crust model (accompanying drawing 2) and upper crust model, the Media density parameter of lower crust and the corresponding model of upper crust is respectively second layer 3300kg/m
3, ground floor is 2700kg/m
3.
Obvious above-described embodiment is only and clearly describes specific embodiment of the invention process.The present embodiment is only the citing illustrating that the present invention does, and the restriction not to embodiment.For those of ordinary skill in the field, on the basis of the above description, other multi-form change or variations can also be made, here without the need to also giving exhaustive to all embodiments.Apparent change amplified thus or variation are still among the protection domain of the invention.
Claims (10)
1., based on a Geophysical Data Processing method for layered geology model, it, for building the geophysical model based on layered medium, comprises the following steps:
(1), raw data obtains, input Bottom surfaces of strata distribution range and depth information;
(2), by this layer model carry out subdivision by plane space scope, subdivision trellis depth value can be obtained by input depth information interpolation;
(3), with subdivision grid for xsect, vertically extend to earth's surface, form the triangulated space of stratified model;
(4) this layer model subdivision grid Geophysical Properties parameter, is given;
(5), by (1)-(4) step set up other layer models, set up to all layer models;
(6) interlayer model Geophysical Properties parameter, is calculated.The inter-level dielectric geophysical parameters that boundary at the bottom of adjacent two layers determines is the Geophysical Properties parameter value sum of stratified model below boundary at the bottom of this inter-level dielectric and this layer.
2. based on a Geophysical Data Processing method for layered geology model, it, for building the geophysical model based on layered medium, is characterized in that, this Geophysical Data Processing method comprises the following steps:
The first step, raw data obtain, and obtain the subsurface formations degree of depth according to geophysical exploration (earthquake, drilling well etc.);
Second step, build the geophysical model based on layered medium based on the depth of stratum Data Data that obtains in the first step, concrete methods of realizing comprises:
Underground medium model is divided into N layer from earth's surface to underground, and from earth's surface, following ground floor is labeled as 1 to N to bottommost layer numbering; I-th (i=1,2 ..., N) layer by layer shape model representation be the underground space of i-th layer of bottom boundary to earth's surface, equivalent layer model geophysical parameters be labeled as Pi (i=1,2 ..., N); , N is positive integer;
I-th layer model process of establishing is: according to the i-th Ceng Di circle degree of depth by its planar distribution coordinate (x, y) subdivision is some arbitrary polygon grids, and with subdivision grid for xsect is set up to surperficial vertical prism, and give each subdivision grid corresponding Geophysical Properties value Pi; In i-th Ceng Di circle to the i-th-1 Ceng Di circle scope, medium Geophysical Properties can be characterized by formula
3rd step, above step set up the layered geology model with specific Geophysical Properties (as density, magnetic susceptibility, resistivity, speed etc.), can be used for geophysical anomaly forward modelling.
3. as claimed in claim 1 or 2 based on the Geophysical Data Processing method of layered geology model, it is characterized in that, each layer model carries out mesh generation with this layer of bottom boundary form and vertically splits this stratified model space.
4. as claimed in claim 1 or 2 based on the Geophysical Data Processing method of layered geology model, it is characterized in that, surface layer is the interface, public top of each stratified model.
5. as Claims 2 or 3 based on the Geophysical Data Processing method of layered geology model, it is characterized in that, the shape that each layer model carries out subdivision grid according to bottom boundary can be arbitrary polygon.
6. as described in Claims 2 or 3 based on the Geophysical Data Processing method of layered geology model, it is characterized in that, each layer model each layer subdivision grid can according to actual conditions self-adaptative adjustment, and each layer mesh generation form does not require identical.
7., as claimed in claim 1 or 2 based on the Geophysical Data Processing method of layered geology model, it is characterized in that, each layer model has corresponding geophysical parameters.
8., as claimed in claims 6 or 7 based on the Geophysical Data Processing method of layered geology model, it is characterized in that, each layer model can have transformable geophysical parameters according to subdivision grid.
9. as claimed in claims 6 or 7 based on the Geophysical Data Processing method of layered geology model, it is characterized in that, geophysical parameters is characterized by density or the magnetic susceptibility of model.
10. as described in claim 1 or claim 7 or 8 based on the Geophysical Data Processing method of layered geology model, it is characterized in that, can be added by the Geophysical Properties of all layer models comprising this interlayer model part by the Geophysical Properties of this interlayer model of the determined interlayer model in boundary at the bottom of adjacent layer and obtain.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410844376.0A CN104535391B (en) | 2014-12-30 | 2014-12-30 | A kind of Geophysical Data Processing method based on layered geology model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410844376.0A CN104535391B (en) | 2014-12-30 | 2014-12-30 | A kind of Geophysical Data Processing method based on layered geology model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104535391A true CN104535391A (en) | 2015-04-22 |
CN104535391B CN104535391B (en) | 2017-01-04 |
Family
ID=52850965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410844376.0A Active CN104535391B (en) | 2014-12-30 | 2014-12-30 | A kind of Geophysical Data Processing method based on layered geology model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104535391B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104866653A (en) * | 2015-04-29 | 2015-08-26 | 中国地质科学院矿产资源研究所 | Method for acquiring underground three-dimensional density structure |
CN105549106A (en) * | 2016-01-07 | 2016-05-04 | 中国科学院地质与地球物理研究所 | Gravity multi-interface inversion method |
CN113917561A (en) * | 2021-09-17 | 2022-01-11 | 东华理工大学 | Method and system for determining target geologic body in three-dimensional geological modeling of mining area |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101661115A (en) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | Standard framework-based method for quick three-dimensional gravity and magnetic inversion for physical properties |
CN103207410A (en) * | 2013-04-17 | 2013-07-17 | 中国海洋石油总公司 | Rugged seabed aimed hybrid grid model building method |
US20130238297A1 (en) * | 2012-03-09 | 2013-09-12 | Schlumberger Technology Corporation | Faulted geological structures containing unconformities |
CN104181585A (en) * | 2014-06-27 | 2014-12-03 | 中国石油化工股份有限公司 | Shear wave estimation method and system in geophysical exploration |
-
2014
- 2014-12-30 CN CN201410844376.0A patent/CN104535391B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101661115A (en) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | Standard framework-based method for quick three-dimensional gravity and magnetic inversion for physical properties |
US20130238297A1 (en) * | 2012-03-09 | 2013-09-12 | Schlumberger Technology Corporation | Faulted geological structures containing unconformities |
CN103207410A (en) * | 2013-04-17 | 2013-07-17 | 中国海洋石油总公司 | Rugged seabed aimed hybrid grid model building method |
CN104181585A (en) * | 2014-06-27 | 2014-12-03 | 中国石油化工股份有限公司 | Shear wave estimation method and system in geophysical exploration |
Non-Patent Citations (3)
Title |
---|
YA XU ET AL.: "Regional gravity anomaly separation using wavelet transform and spectrum analysis", 《J.GEOPHYS.ENG》 * |
Z.ZHANG ET AL.: "Construction of regional and local seismic anisotropic structures from wide-angle seismic data:crustal deformation in the southeast of China", 《J SEISMOL》 * |
徐亚 等: "渤海残留盆地分布综合地球物理研究", 《地球物理学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104866653A (en) * | 2015-04-29 | 2015-08-26 | 中国地质科学院矿产资源研究所 | Method for acquiring underground three-dimensional density structure |
CN104866653B (en) * | 2015-04-29 | 2017-10-24 | 中国地质科学院矿产资源研究所 | Method for acquiring underground three-dimensional density structure |
CN105549106A (en) * | 2016-01-07 | 2016-05-04 | 中国科学院地质与地球物理研究所 | Gravity multi-interface inversion method |
CN113917561A (en) * | 2021-09-17 | 2022-01-11 | 东华理工大学 | Method and system for determining target geologic body in three-dimensional geological modeling of mining area |
CN113917561B (en) * | 2021-09-17 | 2023-06-20 | 东华理工大学 | Method and system for determining target geologic body in three-dimensional geologic modeling of mining area |
Also Published As
Publication number | Publication date |
---|---|
CN104535391B (en) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110286421B (en) | Method for modeling natural fractures of tight sandstone reservoir | |
US10838108B2 (en) | Method and system for interpolating discontinuous functions in a subsurface model | |
CN102222365B (en) | Method for reconstructing curved surface of complex space | |
CN103279986B (en) | Three-dimensional horizontal geologic profile figure preparation method and application thereof | |
CN104574511B (en) | A kind of quick progressive three-dimensional geological modeling method | |
US20160376885A1 (en) | Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs | |
US20140136171A1 (en) | Unstructured Grids For Modeling Reservoirs | |
CN106687827B (en) | Stratum modeling method for fault | |
WO2016070587A1 (en) | Method for looking for oil and gas reservoir on the basis of trap-3d software | |
US10571601B2 (en) | Geologic fault representation for flow simulation | |
CN105701319B (en) | A kind of microfacies modeling method under horizontal well rule well pattern | |
CN104504754B (en) | Oil and gas reservoir multipoint statistical modeling method and device | |
CN104850682A (en) | Multiple-point geostatistics modeling method based on position | |
CN110244021B (en) | Stratum layering method based on anisotropic interpolation | |
CA2543128A1 (en) | Shoulder bed effects removal | |
CN104535391B (en) | A kind of Geophysical Data Processing method based on layered geology model | |
WO2017082870A1 (en) | Fracture network triangle mesh adjustment | |
WO2017151100A1 (en) | Hybrid 3d geocellular representation of selected natural fracture network subsets | |
CN104297792A (en) | Phase control inversion method for waterway reservoir stratums stacked on alluvial fan | |
AU2015414307B2 (en) | Target object simulation using undulating surfaces | |
CN116778096A (en) | Three-dimensional geologic body modeling method based on BIM three-dimensional cloud data | |
CN103473428B (en) | Method and a device for analyzing a space of a mine | |
CN112507615B (en) | Intelligent identification and visualization method for lithofacies of continental tight reservoir | |
CN114820959A (en) | Method for digitally modeling multiple geological elements of rock burst mine | |
CN112052591A (en) | Interlayer fine depicting and embedded modeling method under reservoir configuration constraint |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |