CN104535391B - A kind of Geophysical Data Processing method based on layered geology model - Google Patents
A kind of Geophysical Data Processing method based on layered geology model Download PDFInfo
- Publication number
- CN104535391B CN104535391B CN201410844376.0A CN201410844376A CN104535391B CN 104535391 B CN104535391 B CN 104535391B CN 201410844376 A CN201410844376 A CN 201410844376A CN 104535391 B CN104535391 B CN 104535391B
- Authority
- CN
- China
- Prior art keywords
- model
- layer
- geophysical
- layered
- subdivision
- 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
Abstract
The present invention is a kind of method setting up underground medium model for geophysical computing.This method is assumed based on underground medium stratified model, with boundary at the bottom of each stratum as the end, the multiple layer combination model with overlapping region is set up with surface layer for top, every layer model can set up spatial model each via any subdivision of vertical polygons prism, sets up real laminar formation character by multilayered model physical property attribute simulation.The method can fully realize the fine description of Layered medium over the ground and improve the computational efficiency of geophysical effect.This method can be used for the applications such as the FORWARD AND INVERSE PROBLEMS calculating of geophysics Gravity-Magnetic Survey.
Description
Technical field
The present invention relates to geophysical probing technique, at a kind of geophysical data based on layered geology model
Reason method, the energy complex layered model of fine description also improves geophysical computing efficiency.
Background technology
Geophysics modeling is by the basis of geophysics forward modelling.Modeling method determines underground Media Description
Mode and geophysical computing is had a direct impact.Underground medium is generally of Lamellar character, therefore at geophysical research
In be modeled being the most general method based on layered medium.When possessing multilayer dielectricity, generally by interlayer is situated between
The model facetization of matter describes underground medium character between layers.And due to different layers medium complexity journey in Practical Calculation
Degree difference, in the bottom of model meshes subdivision upper interlayer dielectric and top layer partition patterns and differ, generally this can cause complexity
Geophysical computing, affect geophysics and just drilling and the efficiency of Inversion Calculation.As a example by gravitational inversion, with prism subdivision
Mode calculated gravity anomaly is the simplest effectively, and for adjacent inter-level dielectric, intermediate interface is the end circle of last layer medium
And as the top circle of next layer, if the partition patterns difference of two-layer certainly will cause repeatedly subdivision and the calculating of intermediate layer model,
This have impact on the fine description to intermediate layer character, is also unfavorable for improving computational efficiency.Therefore, it is to avoid to shared layer position repeatedly
Subdivision then can realize the fine description to layer position, and improves computational efficiency, and just drilling geophysics and Inversion Calculation etc. has weight
Want meaning.
Summary of the invention
Based on problems of the prior art and reality need, the present invention provides a kind of, based on layered geology model
Geophysical Data Processing method, its objective is to set up multiple layer combination geophysical model based on layered geology model theory, can
Realize the fine description of multi-layer underground medium, and geophysical computing efficiency can be improved, can be directly used for gravity and magnetic anomalies
Forward modelling and geophysical exploration research.
First technical scheme of the foundation present invention, a kind of Geophysical Data Processing method based on layered geology model,
It, for building geophysical model based on layered medium, comprises the following steps:
(1), initial data obtain, input Bottom surfaces of strata distribution and depth information;
(2), by plane space scope, this layer model being carried out subdivision, subdivision trellis depth value can be inserted by input depth information
It is worth to;
(3), with subdivision grid as cross section, vertically extending to earth's surface, form the triangulated space of stratified model;
(4) this layer model subdivision grid Geophysical Properties parameter, is given;
(5), set up other layer models by (1)-(4) step, set up to all layer models;
(6), interlayer model Geophysical Properties parameter is calculated.The inter-level dielectric geophysics ginseng that boundary at the bottom of adjacent two layers determines
Number is the Geophysical Properties parameter value sum of stratified model below boundary at the bottom of this inter-level dielectric and this layer.
Second technical scheme of the foundation present invention, a kind of Geophysical Data Processing method based on layered geology model,
Its for build geophysical model based on layered medium, it is characterised in that this Geophysical Data Processing method include with
Lower step:
The first step, initial data obtain, and obtain the subsurface formations degree of depth according to geophysical exploration (earthquake, drilling well etc.);
Second step, depth of stratum Data Data based on acquisition in first step structure geophysics mould based on layered medium
Type, concrete methods of realizing includes:
Underground medium model is divided into N shell from earth's surface to underground, from earth's surface following ground floor to bottommost layer numbering be labeled as 1 to
N;I-th (i=1,2 ..., N) layer by layer shape model representation be i-th layer of bottom interface to the underground space on earth's surface, the equivalent layer model earth
Physical parameter be labeled as Pi (i=1,2 ..., N);, N is positive integer;
I-th layer model sets up process: according to the i-th Ceng Di circle degree of depth, by its planar distribution coordinate, (x, y) subdivision is some
Arbitrary polygon grid, and set up the vertical prism to surface with subdivision grid for cross section, and give each subdivision grid
Corresponding Geophysical Properties value Pi;In the range of i-th Ceng Di circle to i-th-1 Ceng Di circle, medium Geophysical Properties may be characterized as public affairs
Formula
3rd step, above step are set up has specific Geophysical Properties (such as density, susceptibility, resistivity, speed
Deng) layered geology model, can be used for geophysical anomaly forward modelling.
Wherein, each layer model is to carry out mesh generation with this layer of bottom interface form and vertically split this stratified model sky
Between.
Preferably, surface layer is the interface, public top of each stratified model.Each layer model carries out subdivision net according to bottom interface
The shape of lattice can be arbitrary polygon.Each layer model each layer subdivision grid can be according to practical situation self-adaptative adjustment, each layer net
Lattice subdivision form is not required for identical.Each layer model has corresponding geophysical parameters.
Further, each layer model can have transformable geophysical parameters according to subdivision grid.Geophysics is joined
Number is characterized as density or the susceptibility of model.Geophysics by this interlayer model of interlayer model determined by boundary at the bottom of adjacent layer
Character can be added by the Geophysical Properties of all layer models including this interlayer model part and obtain.
Each Ceng Di circle and the determined stratified model of earth's surface observation layer are combined obtaining by the modeling method that the present invention uses
Round body underground medium distributed model, simulate underground medium practical situation by the distribution of Geophysical Properties parameter.
Present invention, avoiding the subdivision of interlayer model common interface, modeling method has more adaptivity, is suitable to stratiform and is situated between
The fine description of matter, is conducive to improving computational efficiency.
Accompanying drawing explanation
Fig. 1 is monolayer geophysical model subdivision cross-sectional view;
Fig. 2 is monolayer geophysical model subdivision longitudinal section;
Fig. 3 is two-layer geophysical model subdivision built-up pattern profilograph.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Ground describes, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments.Based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under not making creative work premise
Embodiment, broadly falls into the scope of protection of the invention.Additionally, described embodiment is only that the present invention is expanded on further,
Rather than limitation of the present invention.
First technical scheme of the foundation present invention, a kind of Geophysical Data Processing method based on layered geology model,
It, for building geophysical model based on layered medium, comprises the following steps:
(1), initial data obtain, input Bottom surfaces of strata distribution and depth information;
(2), by plane space scope, this layer model being carried out subdivision, subdivision trellis depth value can be inserted by input depth information
It is worth to;
(3), with subdivision grid as cross section, vertically extending to earth's surface, form the triangulated space of stratified model;
(4) this layer model subdivision grid Geophysical Properties parameter, is given;
(5), set up other layer models by (1)-(4) step, set up to all layer models;
(6), interlayer model Geophysical Properties parameter is calculated.The inter-level dielectric geophysics ginseng that boundary at the bottom of adjacent two layers determines
Number is the Geophysical Properties parameter value sum of stratified model below boundary at the bottom of this inter-level dielectric and this layer.
Second technical scheme of the foundation present invention, a kind of Geophysical Data Processing method based on layered geology model,
Its for build geophysical model based on layered medium, it is characterised in that this Geophysical Data Processing method include with
Lower step:
The first step, initial data obtain, and obtain the subsurface formations degree of depth according to geophysical exploration (earthquake, drilling well etc.);
Second step, depth of stratum Data Data based on acquisition in first step structure geophysics mould based on layered medium
Type, concrete methods of realizing includes:
Underground medium model is divided into N shell from earth's surface to underground, from earth's surface following ground floor to bottommost layer numbering be labeled as 1 to
N;I-th (i=1,2 ..., N) layer by layer shape model representation be i-th layer of bottom interface to the underground space on earth's surface, the equivalent layer model earth
Physical parameter be labeled as Pi (i=1,2 ..., N);, N is positive integer;
I-th layer model sets up process: according to the i-th Ceng Di circle degree of depth, by its planar distribution coordinate, (x, y) subdivision is some
Arbitrary polygon grid, and set up the vertical prism to surface with subdivision grid for cross section, and give each subdivision grid
Corresponding Geophysical Properties value Pi;In the range of i-th Ceng Di circle to i-th-1 Ceng Di circle, medium Geophysical Properties may be characterized as public affairs
Formula
3rd step, above step are set up has specific Geophysical Properties (such as density, susceptibility, resistivity, speed
Deng) layered geology model, can be used for geophysical anomaly forward modelling.
Wherein, each layer model is to carry out mesh generation with this layer of bottom interface form and vertically split this stratified model sky
Between.
Preferably, surface layer is the interface, public top of each stratified model.Each layer model carries out subdivision net according to bottom interface
The shape of lattice can be arbitrary polygon.Each layer model each layer subdivision grid can be according to practical situation self-adaptative adjustment, each layer net
Lattice subdivision form is not required for identical.Each layer model has corresponding geophysical parameters.
Further, each layer model can have transformable geophysical parameters according to subdivision grid.Geophysics is joined
Number is characterized as density or the susceptibility of model.Geophysics by this interlayer model of interlayer model determined by boundary at the bottom of adjacent layer
Character can be added by the Geophysical Properties of all layer models including this interlayer model part and obtain.
Below in conjunction with specific embodiment, the present invention is further explained.
Embodiment 1: carry out implementing underground medium density modeling as a example by gravity anomaly calculating by the modeling of two-layer underground medium.Press
Of the present invention, modeling specifically comprises the following steps that
(1) input initial data, in this example, stratum is two-layer, for upper crust and lower crust.Input data bit upper crust
The bottom interface degree of depth and the lower crust bottom interface degree of depth;
(2) lower crust layer layer model accordingly is set up.According to lower crust distribution of strata data, to this layer of distribution along xy
Plane carries out rectangular mesh subdivision, and this layer depth that each grid represents is determined (see accompanying drawing 1) by input value interpolation;
(3) with Rectangular Partition grid shown in step (2) be cross section set up upright this layer of bottom interface of hexahedron subdivision to ground
Table space scope (accompanying drawing 2 is corresponding vertical section);
(4) density value giving this layer model subdivision grid is 3300kg/m3;
(5) upper crust layer model accordingly is set up.According to upper crust distribution of strata data, this layer of distribution is put down along xy
Face carries out holding mesh generation, and this layer depth that each grid represents is determined by input value interpolation, upper crust subdivision network style
Identical with step (2), layer data distribution in mesh generation big I base area changes, and in this example, upper crust model is corresponding
Model meshes is smaller in size than lower crust model correspondence model meshes size, and (accompanying drawing 1 mid portion size of mesh opening represents upper crust net
Lattice subdivision size);
(6) it is that cross section sets up upright this layer of bottom interface of hexahedron subdivision with rectangular mesh shown in step (5) empty to earth's surface
Between scope (accompanying drawing 3 upper formation model be corresponding vertical section);Setting up ground floor dielectric model by (1)-(4), this layer model is close
Angle value is-600kg/m3
(7) collectively constituted by lower crust model (accompanying drawing 2) and upper crust model there is the geological model of double-deck crustal structure
(accompanying drawing 3), the Media density parameter of lower crust model corresponding with upper crust is respectively second layer 3300kg/m3, ground floor is
2700kg/m3。
Obviously above-described embodiment is only the specific implementation process clearly describing the present invention.The present embodiment only illustrates this
The citing done of invention, and not restriction to embodiment.For those of ordinary skill in the field, state upper
On the basis of bright, it is also possible to make change or the variation of other multi-forms, here without also all embodiments being given
With exhaustive.That is thus amplified obviously changes or changes among still in the protection domain of the invention.
Claims (3)
1. a Geophysical Data Processing method based on layered geology model, for building earth thing based on layered medium
Reason model, it is characterised in that this Geophysical Data Processing method comprises the following steps:
The first step, initial data obtain, and obtain the subsurface formations degree of depth according to geophysical exploration;
Second step, based in the first step obtain depth of stratum data construct geophysical model based on layered medium, specifically
Implementation method includes: underground medium model is divided into N shell from earth's surface to underground, and from earth's surface, following ground floor is to bottommost layer numbering labelling
Be 1 to N;I-th (i=1,2 ..., N) layer by layer shape model representation be i-th layer of bottom interface to the underground space on earth's surface, equivalent layer mould
Type geophysical parameters be labeled as Pi (i=1,2 ..., N), N is positive integer;I-th layer model sets up process: according to i-th layer
By its planar distribution coordinate, (x, y) subdivision is some arbitrary polygon grids to the end circle degree of depth, and builds with subdivision grid for cross section
Found the vertical prism to surface, and give each subdivision grid corresponding Geophysical Properties value Pi;I-th Ceng Di circle to i-th-
In the range of 1 Ceng Di circle, medium Geophysical Properties may be characterized as formula
3rd step, above step set up the layered geology model with specific Geophysical Properties, can be used for geophysics different
Often forward modelling.
2. Geophysical Data Processing method based on layered geology model as claimed in claim 1, it is characterised in that each
Layer model has transformable geophysical parameters according to subdivision grid.
3. Geophysical Data Processing method based on layered geology model as claimed in claim 1, it is characterised in that the earth
Physical parameter is characterized as density of earth formations or the susceptibility of model.
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 CN104535391A (en) | 2015-04-22 |
CN104535391B true 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) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104866653B (en) * | 2015-04-29 | 2017-10-24 | 中国地质科学院矿产资源研究所 | Method for acquiring underground three-dimensional density structure |
CN105549106B (en) * | 2016-01-07 | 2018-06-08 | 中国科学院地质与地球物理研究所 | A kind of gravity multiple solutions inversion method |
CN113917561B (en) * | 2021-09-17 | 2023-06-20 | 东华理工大学 | Method and system for determining target geologic body in three-dimensional geologic modeling of mining area |
Citations (3)
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 |
CN104181585A (en) * | 2014-06-27 | 2014-12-03 | 中国石油化工股份有限公司 | Shear wave estimation method and system in geophysical exploration |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2987903B1 (en) * | 2012-03-09 | 2014-05-09 | Schlumberger Services Petrol | GEOLOGICAL FAILURE STRUCTURES CONTAINING NONCONFORMITIES. |
-
2014
- 2014-12-30 CN CN201410844376.0A patent/CN104535391B/en active Active
Patent Citations (3)
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 |
CN104181585A (en) * | 2014-06-27 | 2014-12-03 | 中国石油化工股份有限公司 | Shear wave estimation method and system in geophysical exploration |
Non-Patent Citations (3)
Title |
---|
Construction of regional and local seismic anisotropic structures from wide-angle seismic data:crustal deformation in the southeast of China;Z.Zhang et al.;《J Seismol》;20080904;第13卷;241-252 * |
Regional gravity anomaly separation using wavelet transform and spectrum analysis;Ya Xu et al.;《J.Geophys.Eng》;20090710;第6卷;279-287 * |
渤海残留盆地分布综合地球物理研究;徐亚 等;《地球物理学报》;20070531;第50卷(第3期);868-881 * |
Also Published As
Publication number | Publication date |
---|---|
CN104535391A (en) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10838108B2 (en) | Method and system for interpolating discontinuous functions in a subsurface model | |
CN102222365B (en) | Method for reconstructing curved surface of complex space | |
CN110286421B (en) | Method for modeling natural fractures of tight sandstone reservoir | |
CN104574511B (en) | A kind of quick progressive three-dimensional geological modeling method | |
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 | |
EP2686714B1 (en) | Grid from depositional space | |
US20140136171A1 (en) | Unstructured Grids For Modeling Reservoirs | |
US20130218539A1 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
CN106687827B (en) | Stratum modeling method for fault | |
CN106981093A (en) | A kind of three-dimensional formation parallel modeling method of subregion constraint coupling | |
CN104504754B (en) | A kind of method and device of oil and gas reservoir multi-point statistic modeling | |
EP2631685A2 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
US11041976B2 (en) | Method and system for creating and using a subsurface model in hydrocarbon operations | |
CN105513131A (en) | Three-dimensional geological structure model angular point grid dissection method considering folds | |
CN104535391B (en) | A kind of Geophysical Data Processing method based on layered geology model | |
CN112562078A (en) | Three-dimensional geological analysis prediction model construction method | |
CN104297792B (en) | The phased inversion method of water channel reservoir is stacked on a kind of fan | |
CN109063324A (en) | Finite element crustal stress analogue technique method based on Corner-point Grids | |
AU2015414307B2 (en) | Target object simulation using undulating surfaces | |
CN114820959A (en) | Method for digitally modeling multiple geological elements of rock burst mine | |
CN106294900A (en) | A kind of structure and methods for using them of the digital rock mass for drilling well | |
CN112780252B (en) | Geosteering model and forward modeling method suitable for well logging while drilling forward modeling | |
US20220171091A1 (en) | Hydrocarbon flow simulation | |
CN106324668B (en) | A kind of thin seismic reservoir the Forward Modelings based on double change Review of geologic model building techniques | |
WO2017082871A1 (en) | Target object simulation using orbit propagation |
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 |