CN107909649A - A kind of plate scale new 3 D geological model construction method - Google Patents
A kind of plate scale new 3 D geological model construction method Download PDFInfo
- Publication number
- CN107909649A CN107909649A CN201711474429.4A CN201711474429A CN107909649A CN 107909649 A CN107909649 A CN 107909649A CN 201711474429 A CN201711474429 A CN 201711474429A CN 107909649 A CN107909649 A CN 107909649A
- Authority
- CN
- China
- Prior art keywords
- earth
- crust
- dimensional
- plate
- rock
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Abstract
The invention discloses a kind of plate scale new 3 D geological model construction method, consider earth surface arc geometric shape under plate scale, internal height fluctuating layering interfaces and three-dimensional rock mechanics parameters non-the features such as changing, provide fine geologic basis and reference, the accuracy of research precision and its result is improved, a new thinking and approach is opened for plate dimensional configurations dynamic analysis.
Description
Technical field
The present invention relates to belong to microcomputer modelling technical field, and in particular, to a kind of plate scale three-dimensional geological model
Construction method.
Background technology
Plate be epigeosphere not of uniform size, the thickness very little that is constructed that active belt is divided into, area it is very big, by the earth
The curved block of surface profile.Plate dimensional configurations dynamic analysis is the importance of architectonics research, due to research
Scope is big, complex geologic conditions are changeable, carries out correlative study using numerical simulation means more, and the order of accuarcy of analog result depends on
In the fine degree of three-dimensional geological model.
In existing plate scale three-dimensional geological model, be mostly flat plate model or rock mechanics uniform soft soil base, exist with
Lower both sides defect：(1) the objective of plate scale arc lower surface configuration and height fluctuating interior laminate layer interface is not accounted for
It is true；(2) not the characteristics of not accounting for the three dimensions heterogeneous change of rock mechanics parameters.Thus, the plate based on numerical simulation
Block scale dynamic analysis result is relative to the situation of necessary being often with there is larger deviation.
In consideration of it, be necessary for conventional art there are the problem of, there is provided a kind of structure of plate scale three-dimensional geological model
Construction method, considers earth surface arc geometric shape under plate scale, internal height fluctuating layering interfaces and three-dimensional rock
Mechanics parameter non-the features such as changing, a new thinking and approach is opened for plate dimensional configurations dynamic analysis, there is provided
Fine geologic basis and reference, improve the accuracy of research precision and its result.
The content of the invention
It is three-dimensional it is an object of the invention in view of the above-mentioned problems, overcome the shortage of prior art, disclose a kind of plate scale
Geological model construction method, provides a new thinking and approach for plate dimensional configurations dynamic analysis, improves research
The accuracy of precision and its result.
To achieve the above object, the present invention uses following technical scheme：
A kind of plate scale new 3 D geological model construction method, it is characterised in that comprise the following steps：
Step 1, determine research area's three dimensions scope, specifically includes in the longitude and latitude and longitudinal extent in lateral extent
The earth's crust/earth mantle boundary line depth, wherein agreement longitude westwards gradually becomes -180 °, gradually become eastwards using 0 ° of warp as boundary
180 °, be warp-wise positive direction by west to east；Latitude is again -90 °~0 °~90 ° changes to the arctic by the South Pole to equator, by the South Pole
It is broadwise positive direction to the arctic；Longitudinal direction is positive direction by earth's surface to the earth's core
In step 2, lateral extent, will research area according to graticules precision transverse grid, be denoted as XiYj (i, j=1,2,
3……)；In longitudinal extent, research area is divided into sedimentary, upper crust, the middle earth's crust and lower crust from top to bottom, is denoted as Zk (k
=1,2,3 ...), wherein Xi is warp-wise, and Yj is broadwise, and Zk is longitudinal direction；
Step 3, establish basic database, including each depth of seam division (Hxyz) of the earth's surface elevation of each grid and longitudinal direction, horizontal stroke
Longitudinal wave velocity (Pxyz/Sxyz), rock density (Dxyz) data；
Step 4, follow earth geometric shape, establishes ball system coordinate screen work；
Step 5, the layering interfaces data (Hxyz) of all grids in area are studied in extraction from basic database, as basic
Control point, under ball system coordinate screen work, establishes the fundamental space curved surface in research area；
Step 6, under ball system coordinate screen work, establish what is limited by fundamental space curved surface, by sedimentary, upper crust, middle
The three-dimensional entity model of shell and lower crust composition；
Step 7, according to specific trellis precision by three-dimensional entity model mesh generation into limited a element entity (Vxyz)；
Step 8, from master database extraction research area in all grids horizontal longitudinal wave velocity data (Pxyz/Sxyz) and
Rock density (Dxyz), calculates the rock mechanics parameters such as Young's modulus (E), the Poisson's ratio (υ) of each grid；
Step 9, each element entity sequence number of matched indicia, will be each single with calculating the rock mechanics parameters sequence number obtained
The rock mechanics parameters calculated in first entity are loaded into the corresponding unit entity of three-dimensional entity model, and structure is based on ball system coordinate
Screen work and the refined model of the three-dimensional heterogeneous change of rock mechanics parameters.
Preferably, the fundamental space curved surface in the step 5 include earth surface, sedimentary bottom surface, upper crust bottom surface, in
Earth's crust bottom surface and lower crust bottom surface.
Beneficial effects of the present invention are：The present invention considers earth surface arc geometric shape, inside under plate scale
Height fluctuating layering interfaces and rock mechanics parameters three dimensions non-the characteristics of changing, it is proposed that a kind of fine plate scale
New 3 D geological model construction method.
A new thinking and approach is provided for plate dimensional configurations dynamic analysis, improves research precision and its knot
The accuracy of fruit.
Brief description of the drawings
Fig. 1 is scale three-dimensional geological model flow chart of the present invention；
Fig. 2 is African Plate three-dimensional entity model schematic diagram of the present invention；
Fig. 3 is African Plate three-dimensional geological model schematic diagram of the present invention.
Embodiment
The present invention is further illustrated below in conjunction with the accompanying drawings.
The present invention's simultaneously establishes database from establishment research area's scope and gridding, collection research area related data, imports
Spatial three-dimensional dispersion model is established after ball system coordinate screen work and mechanics parameter, builds plate scale three-dimensional geological model.
As described in Figure 1, the present invention specifically includes following aspect：
Step 1, determine research area's three dimensions scope, specifically include in the longitude and latitude and longitudinal extent in lateral extent
The earth's crust/earth mantle boundary line depth.For convenience, arrange：Longitude westwards gradually becomes -180 °, eastwards gradually using 0 ° of warp as boundary
It is changed into 180 °, is warp-wise positive direction by west to east；Latitude is again -90 °~0 °~90 ° changes to the arctic by the South Pole to equator, by
The South Pole to the arctic be broadwise positive direction；Longitudinal direction is positive direction by earth's surface to the earth's core.
Such as：African Plate lateral extent after simplification draws a circle to approve (form by following 25 longitudes and latitudes points：Sequence number, longitude/X,
Latitude/Y)：
3.70 °, 11.90 ° of Isosorbide-5-Nitrae
2,34.20 °, 27.10 °
3,35.50 °, 34.00 °
4,23.70 °, 34.40 °
5,20.00 °, 38.20 °
6, -10.00 °, 36.00 °
7, -16.00 °, 37.80 °
8, -24.10 °, 36.82 °
9, -29.10 °, 39.30 °
10, -44.90 °, 23.70 °
11, -45.00 °, 15.30 °
12, -24.60 °, -1.20 °
13, -15.80 °, -1.40 °
14, -13.30 °, -32.30 °
15, -16.50 °, -41.70 °
16,0.80 °, -54.90 °
17,15.10 °, -52.20 °
18,27.80 °, -52.80 °
19,69.30 °, -25.80 °
20,66.40 °, -12.88 °
21,68.80 °, -5.24 °
22,66.80 °, 2.30 °
23,57.20 °, 9.90 °
24,58.30 °, 12.70 °
25,54.50 °, 15.2 °
Step 2, in lateral extent, African Plate will be studied into area's transverse grid according to 1 ° × 1 ° graticules precision；It is vertical
In the range of, African Plate is divided into sedimentary, upper crust, the middle earth's crust and lower crust from top to bottom, be denoted as Zk (k=1,2,
3,4).Thus, some grid can be labeled as XiYj Zk in African Plate (Xi is warp-wise, and Yj is broadwise, and Zk is longitudinal direction).
Step 3, the ground of each grid determined by step 2 in African Plate basic database, including the plate is established
The data such as each depth of seam division (Hxyz) of table elevation and longitudinal direction, horizontal longitudinal wave velocity (Pxyz/Sxyz), rock density (Dxyz).
Step 4, on the basis of earth mean radius 6371km, it then follows spherical geometric shape, establishes ball system coordinate screen work.
Step 5, layering circle of all grids in African Plate is extracted from the African Plate basic database that step 3 is established
Face data (Hxyz), as basic control point, under the ball system coordinate screen work that step 4 is established, establish the substantially empty of African Plate
Between curved surface, including：Earth surface, sedimentary bottom surface, upper crust bottom surface, middle earth's crust bottom surface and lower crust bottom surface.
Step 6, under the ball system coordinate screen work that step 4 is established, it is bent to establish the African Plate fundamental space established by step 5
Face limitation, the African Plate three-dimensional entity model that the sedimentary that is determined by step 2, upper crust, the middle earth's crust and lower crust form.
Step 7, according to step 2 specific trellis precision by African Plate three-dimensional entity model mesh generation into limited a entity
Unit (Vxyz).
Step 8, the transverse and longitudinal ripple of all grids in African Plate is extracted from the African Plate master database that step 3 is established
Velocity of wave data (Pxyz/Sxyz) and rock density (Dxyz), calculate Young's modulus (E), Poisson's ratio (υ) of each element entity etc.
Rock mechanics parameters.The rock mechanics parameters for calculating each element entity obtained are labeled as：Young's modulus (Exyz), Poisson's ratio
(υxyz)。
The calculating of rock mechanics parameters, with reference to following equation：
In formula：E is rock Young's modulus, and υ is Rock Poisson Ratio Using, and ρ is rock density, VpFor rock longitudinal wave velocity, VsFor rock
Stone transverse wave speed.
Step 9, each element entity sequence number that matching step 7 marks calculates the rock mechanics parameters sequence obtained with step 8
Number, the rock mechanics parameters of each calculating are loaded into the corresponding unit entity of three-dimensional entity model, structure is sat based on ball system
Style frame and the refined model of the three-dimensional heterogeneous change of rock mechanics parameters.
Fig. 2 and Fig. 3 is the African Plate three-dimensional entity model schematic diagram and African Plate three-dimensional geological after present invention modeling
Model schematic, compared with prior art, the present invention advantage is：Consider earth geometric shape and three-dimensional rock mechanics is non-
The objective fact of homogeneous, the present invention provides a kind of plate scale new 3 D geological model construction method, with existing correlation technique
Compare, overcome tablet and the drawback such as homogeneous in plate scale geological model structure, analyzed for tectonodynamics and provide one
Bar new thinking and approach, improve the accuracy of research precision and result.
Above-described embodiment, simply presently preferred embodiments of the present invention, is not used for limiting the scope of the present invention, therefore all with this
The equivalence changes that content described in invention claim is done, should all be included within scope of the invention as claimed.
Claims (2)
1. a kind of plate scale new 3 D geological model construction method, it is characterised in that comprise the following steps：
Step 1, determine research area's three dimensions scope, specifically includes the ground in the longitude and latitude and longitudinal extent in lateral extent
Shell/earth mantle boundary line depth, wherein agreement longitude westwards gradually becomes -180 °, gradually becomes 180 ° eastwards using 0 ° of warp as boundary,
It is warp-wise positive direction by west to east；Latitude is again -90 °~0 °~90 ° changes to the arctic by the South Pole to equator, by the South Pole to the arctic
For broadwise positive direction；Longitudinal direction is positive direction by earth's surface to the earth's core；
In step 2, lateral extent, will research area according to graticules precision transverse grid, be denoted as XiYj (i, j=1,2,
3……)；In longitudinal extent, research area is divided into sedimentary, upper crust, the middle earth's crust and lower crust from top to bottom, is denoted as Zk (k
=1,2,3 ...), wherein Xi is warp-wise, and Yj is broadwise, and Zk is longitudinal direction；
Step 3, establish basic database, including each depth of seam division (Hxyz) of the earth's surface elevation of each grid and longitudinal direction, transverse and longitudinal ripple
Velocity of wave (Pxyz/Sxyz), rock density (Dxyz) data；
Step 4, follow earth geometric shape, establishes ball system coordinate screen work.；
Step 5, the layering interfaces data (Hxyz) of all grids in area are studied in extraction from basic database, as basic control
Point, under ball system coordinate screen work, establishes the fundamental space curved surface in research area；
Step 6, under ball system coordinate screen work, establish what is limited by fundamental space curved surface, by sedimentary, upper crust, the middle earth's crust and
The three-dimensional entity model of lower crust composition；
Step 7, according to specific trellis precision by three-dimensional entity model mesh generation into limited a element entity (Vxyz)；
Step 8, the horizontal longitudinal wave velocity data (Pxyz/Sxyz) of all grids and rock in area are studied in extraction from master database
Density (Dxyz), calculates the rock mechanics parameters such as Young's modulus (E), the Poisson's ratio (υ) of each grid；
Step 9, each element entity sequence number of matched indicia are real by each unit with calculating the rock mechanics parameters sequence number obtained
The rock mechanics parameters calculated in vivo are loaded into the corresponding unit entity of three-dimensional entity model, and structure is based on ball system coordinate screen work
With the refined model of the three-dimensional heterogeneous change of rock mechanics parameters.
A kind of 2. plate scale new 3 D geological model construction method according to claim 1, it is characterised in that the step
Fundamental space curved surface in 5 includes earth surface, sedimentary bottom surface, upper crust bottom surface, middle earth's crust bottom surface and lower crust bottom surface.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711474429.4A CN107909649A (en) | 2017-12-29 | 2017-12-29 | A kind of plate scale new 3 D geological model construction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711474429.4A CN107909649A (en) | 2017-12-29 | 2017-12-29 | A kind of plate scale new 3 D geological model construction method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107909649A true CN107909649A (en) | 2018-04-13 |
Family
ID=61871966
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711474429.4A Pending CN107909649A (en) | 2017-12-29 | 2017-12-29 | A kind of plate scale new 3 D geological model construction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107909649A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109509397A (en) * | 2018-11-13 | 2019-03-22 | 浙江大学 | A kind of data processing method for the experiment of hypergravity geologic structure physical modelling |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102945570A (en) * | 2012-11-23 | 2013-02-27 | 华东师范大学 | Method for constructing full-space three-dimensional digital earth model |
| CN103345774A (en) * | 2013-07-17 | 2013-10-09 | 中国人民解放军第三军医大学 | Method for building three-dimensional multi-scale vectorization model |
| CN103838936A (en) * | 2014-03-25 | 2014-06-04 | 中国石油大学(华东) | High-precision tectonic stress field simulation method applicable to turbidite sand low-permeability reservoirs |
| CN106680104A (en) * | 2016-11-29 | 2017-05-17 | 西南石油大学 | Dynamic and static rock mechanical parameter correction method and device |
-
2017
- 2017-12-29 CN CN201711474429.4A patent/CN107909649A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102945570A (en) * | 2012-11-23 | 2013-02-27 | 华东师范大学 | Method for constructing full-space three-dimensional digital earth model |
| CN103345774A (en) * | 2013-07-17 | 2013-10-09 | 中国人民解放军第三军医大学 | Method for building three-dimensional multi-scale vectorization model |
| CN103838936A (en) * | 2014-03-25 | 2014-06-04 | 中国石油大学(华东) | High-precision tectonic stress field simulation method applicable to turbidite sand low-permeability reservoirs |
| CN106680104A (en) * | 2016-11-29 | 2017-05-17 | 西南石油大学 | Dynamic and static rock mechanical parameter correction method and device |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109509397A (en) * | 2018-11-13 | 2019-03-22 | 浙江大学 | A kind of data processing method for the experiment of hypergravity geologic structure physical modelling |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109753707B (en) | Method for extracting stratum boundary line by utilizing exploration line profile to develop three-dimensional modeling | |
| CN103454686B (en) | Small scale sedimentary facies based on strata slicing carries out the method and system of reservoir prediction | |
| CN104866653B (en) | Method for acquiring underground three-dimensional density structure | |
| CN103838936B (en) | It is applicable to the high accuracy tectonic stress field analogy method of turbid long-pending sand low permeability reservoir | |
| CN107944086B (en) | Rapid modeling method based on drilling data | |
| CN105974479B (en) | The chromatography 2D/3D anisotropy Depth Domain velocity modeling method of GPU space lattice body | |
| CN105386756B (en) | A method of brittle formation porosity is calculated using dependent variable | |
| CN109061728B (en) | A kind of ach-bar sandbody detailed predicting method | |
| CN110286421A (en) | A kind of tight sandstone reservoir intrinsic fracture modeling method | |
| CN103514631A (en) | Gridding type parallel construction method for three-dimensional geological model based on massive borehole data | |
| CN110244352A (en) | A kind of CRUSTAL THICKNESS gravitational inversion method based on variable density | |
| CN104574513B (en) | The accurate characterizing method for portraying interlayer distributed in three dimensions based on geometric shape | |
| CN107766978A (en) | The intelligent optimization method of skewed pattern | |
| CN103942842A (en) | Embedded type meandering river sand modeling method | |
| CN111754623A (en) | Discrete fracture modeling method based on multi-point geostatistics | |
| CN105652341B (en) | A kind of quick quantification palaeogeomorphology for considering paleostructure influence | |
| CN105022856B (en) | Predict the reservoir modeling methodologies of high camber meandering channel reservoir internal structure | |
| CN107939372A (en) | For the optimal well location dispositions method and device of craven fault block reservoirs | |
| CN107909649A (en) | A kind of plate scale new 3 D geological model construction method | |
| CN110863818A (en) | Description method and device for residual oil/gas distribution | |
| CN112346139A (en) | Gravity data multilayer equivalent source continuation and data conversion method | |
| CN104462649B (en) | A kind of automatic update method of ore body block segment model reserves | |
| CN110428497A (en) | Braided stream training image generation method | |
| CN110287560A (en) | A kind of complexity form earth material field excavated volume calculation method | |
| CN115392032A (en) | GIS-MPM seamless integrated dynamic three-dimensional geological model construction method |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180413 |
|
| RJ01 | Rejection of invention patent application after publication |