CN114332391A - Three-dimensional geologic body probability model modeling method - Google Patents
Three-dimensional geologic body probability model modeling method Download PDFInfo
- Publication number
- CN114332391A CN114332391A CN202111634129.4A CN202111634129A CN114332391A CN 114332391 A CN114332391 A CN 114332391A CN 202111634129 A CN202111634129 A CN 202111634129A CN 114332391 A CN114332391 A CN 114332391A
- Authority
- CN
- China
- Prior art keywords
- stratum
- probability
- borehole
- probability model
- geologic body
- 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
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000005553 drilling Methods 0.000 claims abstract description 19
- 230000008859 change Effects 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000011160 research Methods 0.000 claims description 10
- 230000007704 transition Effects 0.000 claims description 6
- 238000004088 simulation Methods 0.000 abstract description 13
- 238000012795 verification Methods 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 2
- 230000015572 biosynthetic process Effects 0.000 description 29
- 238000005755 formation reaction Methods 0.000 description 29
- 238000012546 transfer Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000019771 cognition Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000002948 stochastic simulation Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000010205 computational analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a three-dimensional geologic body probability model modeling method, which is characterized by comprising the following steps: collecting real geological drilling stratum information, and performing interval division on a drilling elevation range; acquiring a stratum sequence state sequence of the drilled hole to obtain a stratum change rule of a real drilled hole; establishing a probability model of the stratum thickness, and generating a virtual borehole through random simulation; establishing a total probability model of the three-dimensional geologic body; determining a stratum demarcation point and a stratum section line probability according to the verification drill hole; and obtaining a virtual borehole according to the most probable stratum section line to generate a three-dimensional geologic body maximum probability model. The three-dimensional geologic body probability model constructed based on the probability analysis and random simulation method has the advantages that the randomness and various possibilities of determining the stratum section lines in the modeling process are considered, the accuracy of three-dimensional geologic body modeling can be improved, and the reliability is high.
Description
Technical Field
The invention belongs to the technical field of three-dimensional geological modeling, and particularly relates to a three-dimensional geological model probability modeling method.
Background
The stratigraphic layers of a geologic body are the result of the combined action of multiple factors on the earth during long-term historical evolution. It is not only affected by the laws of regional deposition, but also by external effects such as crustal movement, climate change, human activities, etc. The distribution of the formations is complex and random. In the range from the bedrock surface to the earth surface, a certain stratum can only appear once, and can also appear or not appear for many times, such as special geological phenomena of inversion, repetition, deletion, fault and the like; meanwhile, the thickness of the stratum can be greatly different in scale according to different research areas, such as pinch-out geological phenomena. And because different stratums have great difference in physical and chemical properties, the method has important influence on site selection, design, construction and the like of engineering construction, and therefore, the determination of the stratum in the research area has important significance on underground space analysis and decision. With the development of computer technology and computer graphics technology, it has become a current research hotspot to establish a three-dimensional geological model through limited geological drilling to visually display the stratum condition.
The three-dimensional geological model established by the traditional method needs to be combined with an interpolation method to determine the positions of the demarcation points of different stratums of the virtual drilling hole, and then the demarcation points of the real drilling hole and the virtual drilling hole which belong to the same stratum are sequentially connected to obtain a geological section line; and then adjusting and determining the stratum section line of the model according to the guidance of expert experience. Due to deviation in measurement, difference of interpolation methods, uncertainty and complexity of formation attribute spatial distribution and difference of formation rule cognition, connection lines of formation profiles have randomness and multiple possibilities, and effective evaluation on the formation cannot be achieved. The probability method has inherent advantages in the aspect of processing the problems of randomness and probability, and can make up for the defects of the traditional method to a great extent. But at present, the research on the three-dimensional geological probability model is not much at home and abroad. And the method is a good idea for researching the stratigraphic section line by performing three-dimensional geological modeling from the angle of probability. Therefore, in view of the current situation, it is urgently needed to establish a three-dimensional geologic probability modeling method to meet the needs of practical use.
Disclosure of Invention
The invention aims to provide a three-dimensional geological model probability modeling method to solve the problems in the background technology.
In order to achieve the purpose, the invention provides the following technical scheme:
(1) collecting real geological drilling information of a research area, and performing interval division on drilling elevations to obtain a stratum state sequence of drilling holes;
(2) according to the Markov chain principle, solving a probability transition matrix of the real borehole, and analyzing a stratum change rule of the real borehole;
(3) establishing a stratum probability model, solving a probability transition matrix of the virtual drilling, and analyzing a stratum change rule of the virtual drilling position;
(4) randomly simulating and obtaining the stratum state of the virtual borehole;
(5) establishing an overall probability model according to a randomly simulated virtual borehole stratum state result;
(6) and determining the probability of the stratum demarcation point and the probability of the stratum section line according to the verified borehole analysis to obtain a maximum probability model which accords with the reality.
The invention has the following advantages:
the three-dimensional geologic body probability model established by the invention integrates random simulation and probability into the modeling process, can reduce the difference selected by a measurement and interpolation method, the uncertainty of stratum attribute space distribution and the deviation caused by stratum rule cognition difference, considers the randomness and various possibilities of the connecting line of a stratum section and carries out probability analysis on the positions of a stratum dividing point and a stratum sectioning line; and the maximum probability model which accords with the research region is determined through probability analysis, so that the rationality of the three-dimensional geological model can be improved.
Drawings
FIG. 1 is a flow chart of the three-dimensional geologic probability model modeling of the present invention.
FIG. 2 is a schematic diagram of a sequence of formation states in accordance with the present invention.
FIG. 3 is a flow chart of a virtual borehole formation sequence generation stochastic simulation of the present invention.
FIG. 4 is a schematic representation of the formation property probability calculation of the present invention.
FIG. 5 is a flow chart of an iterative calculation of virtual borehole depth given formation property probabilities
FIG. 6 is a plot of formation sectioning probability determination in accordance with the present invention.
Detailed Description
The technical scheme of the invention is further explained by combining the attached drawings 1-4 as follows:
(1) and collecting real geological drilling information of the research area, including information of stratum thickness, earth surface elevation, stratum depth elevation, drilling names, coordinates, a bottom layer and the like, and determining a uniform stratum sequence of the research area.
(2) Traverse all boreholes of the investigation region, obtainTaking the maximum elevation Z of the depth range of the borehole in the investigation regionmaxAnd minimum value Zmin. Then the depth interval [ Zmin,Zmax]Is equally divided into n segments with the breaking distance of Then from top to bottom we will divide the interval [ Z ] separatelymin,Zmin+d]、[Zmin,Zmin+2d]、[Zmin,Zmin+3d]、……、[Zmax-d,Zmax]Marked as Q1、Q2、Q3、……、Qn. According to the stratum sequence of each borehole, the stratum state of each borehole in different depth intervals, namely the stratum state sequence, can be obtained. Drilling in the interval Q1、Q2、Q3、……、QnThe corresponding stratum state is marked as I1、I2、I3、……、InIn which In(n-1, 2,3, … …) is the formation, as shown in fig. 2. For a borehole without formation information corresponding to a depth interval, we specify that the formation corresponding to the depth interval is "empty" and marked as "0". For example, if the starting depth h of a borehole is1Greater than ZminIn the depth interval [ Z ]min,h1]If the soil does not exist, marking the stratum state at the moment as 0; same if the end depth h of a bore hole2Less than ZmaxIn the depth interval [ h ]1,Zmax]And the data of the soil state is lacked, so that the stratum state at the moment is marked as 0 for the convenience of computational analysis of the model.
(3) According to the Markov chain principle, solving a probability transition matrix of the real borehole, and analyzing a stratum change rule of the real borehole; probability transition matrix PuThe calculation formula is determined according to equations (1) and (2):
Pij=Mij/Mi (2)
where θ is 1,2,3, … …, representing the step size; pijRepresenting the probability of formation i mutating to formation j; miThe number of samples of the stratum at a certain depth of the borehole in the state i; mijAnd (4) transferring the formation state i to the formation state j through theta steps for a certain depth of the borehole.
(4) Counting formation information, fitting probability distribution of the formation information, establishing a formation probability model, and determining specific parameters, such as mean value and variance, of the corresponding probability distribution model;
(5) generating t (t >0) virtual drill holes VB [ i ] (i is 1,2, … …, t), determining a probability transfer matrix at the positions of the virtual drill holes according to the weight coefficients of the drill holes, and acquiring the stratum change probability of the virtual drill holes; the determination of the weighting coefficients and the probability transfer matrix for the virtual borehole can be represented by equations (3) and (4).
Wherein d isiThe distance between the virtual geological borehole and the real geological borehole point; m is the number of real drilling points; β i is a weight coefficient; u is the power of the power and can be 1 or 2;a theta step probability transfer matrix for the virtual borehole;a theta step probability transfer matrix of the mth real drilling hole; and m is the number of real drilled holes.
(6) The 1 st stochastic simulation was performed according to the monte carlo principle, as shown in fig. 3. According to the unified stratum sequence of the research area, starting from the stratum 1, judging the stratum with the maximum virtual drilling probability, and then generating a thickness random number according to a thickness probability distribution function corresponding to the stratum with the maximum probability; next, judging the stratum with the highest probability of the next layer, and if the stratum with the highest probability of the next layer is not empty, generating a thickness random number according to a thickness probability distribution function corresponding to the stratum; … …, respectively; when the next stratum is empty, the random simulation is finished, and the stratum state of the 1 st randomly simulated virtual borehole can be obtained;
(7) the 2 nd, 3 rd and … … th random simulations were performed in this order until the predetermined number of random simulations was reached, and the simulation was terminated. At the moment, the stratum state of multiple random simulation of virtual drilling can be achieved.
(8) And (5) repeating the steps (5) to (7) to obtain the formation states of the rest virtual boreholes through multiple random simulation.
(9) Generating a plurality of three-dimensional geologic body models according to the stratum state of the virtual borehole obtained by each random simulation, wherein the three-dimensional geologic body models generated by the random simulation jointly form an overall probability model of the three-dimensional geologic body;
(10) and according to different stratum demarcation points of the real drilling holes or the adjacent real drilling holes (such as excavation of new geological drilling holes and generation of new drilling holes according to a geological profile), counting and solving stratum attribute probability values at the stratum demarcation points of the virtual drilling holes obtained by random simulation. The stratum attribute probability refers to the probability value of a point in space belonging to a stratum, and is closely related to the spatial uncertainty of the stratum boundary position. Assume that a region is composed of strata 1 and 2, where AD and BC are possible geological profiles obtained by random modeling, and point F is the probability of a certain geological point to be determined, as shown in FIG. 4. If the section line is AD, F belongs to the stratum 1; if the section line is BC, F belongs to the formation 2. Assuming that the number of random simulations is d, we can obtain d possible geological section lines. We can solve the probability P that the point F belongs to the stratum 1 by only calculating the number of times that the point F belongs to the stratum 1, assuming γ, as shown in equation (5). The probability that the point F belongs to the formation 2 can also be derived in the same way.
(11) For the probability determination of the stratum section line, a virtual borehole or an adjacent real borehole is selected as a verification borehole Z1, and the stratum property probability P of the verification borehole at the stratum demarcation point is solvedz1The probability of the formation profile line is estimated. Then, the probability of the virtual borehole stratum attribute on the section is solved in an iterative mode to be Pz1The position of the point in time (as shown in fig. 5) and connecting the point in time, and the obtained connecting line can be used as the section line of the stratum. The verification borehole Z2 may then continue to be selected (assuming Z2 is located to the right of Z2), as shown in fig. 6. Verifying formation property probability P of borehole Z2 at formation demarcation point by solvingz2The probability of the zone verifying the formation profile to the right of borehole Z2 is estimated. And the probability of the formation profile line of the virtual drill between the verification borehole Z1 and the verification borehole Z2 is taken as Pz1And Pz2Half of the sum of (a). Similar methods can be used for determining the formation section lines of 3 or more verification boreholes, and the description is omitted.
(12) And obtaining stratum information of a plurality of virtual drill holes according to the determined stratum section line, performing three-dimensional geological modeling according to the real geological drill holes and the virtual drill holes, and rendering a map through three-dimensional software to finally obtain a three-dimensional geological body maximum probability model which accords with the reality.
Claims (1)
1. A three-dimensional geologic body probability model modeling method comprises the following steps:
(1) collecting real geological drilling information of a research area, and performing interval division on drilling elevations to obtain a stratum state sequence of drilling holes;
(2) according to the Markov chain principle, solving a probability transition matrix of the real borehole, and analyzing a stratum change rule of the real borehole;
(3) establishing a stratum probability model, solving a probability transition matrix of the virtual drilling, and analyzing a stratum change rule of the virtual drilling position;
(4) randomly simulating and obtaining the stratum state of the virtual borehole;
(5) establishing an overall probability model according to a randomly simulated virtual borehole stratum state result;
(6) and determining the probability of the stratum demarcation point and the probability of the stratum section line according to the verified borehole analysis to obtain a maximum probability model which accords with the reality.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111634129.4A CN114332391A (en) | 2021-12-22 | 2021-12-22 | Three-dimensional geologic body probability model modeling method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111634129.4A CN114332391A (en) | 2021-12-22 | 2021-12-22 | Three-dimensional geologic body probability model modeling method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114332391A true CN114332391A (en) | 2022-04-12 |
Family
ID=81016469
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111634129.4A Pending CN114332391A (en) | 2021-12-22 | 2021-12-22 | Three-dimensional geologic body probability model modeling method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114332391A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116958470A (en) * | 2023-07-25 | 2023-10-27 | 中山大学 | Geological modeling method and device integrating Markov chain and multipoint statistics |
CN117934767A (en) * | 2024-03-22 | 2024-04-26 | 华中科技大学 | Three-dimensional geological profile generation method and device and geological information modeling system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160024232A (en) * | 2014-08-25 | 2016-03-04 | 서울대학교산학협력단 | 3-Dimensional Space Modeling Method based on the Geotechnical Information |
CN106194181A (en) * | 2016-08-08 | 2016-12-07 | 西安科技大学 | Intelligent work surface coal-rock interface identification method based on geologic data |
CN110176070A (en) * | 2019-05-07 | 2019-08-27 | 中山大学 | A kind of bed succession analogy method carrying out machine learning based on borehole data |
CN111127630A (en) * | 2019-12-17 | 2020-05-08 | 核工业北京地质研究院 | Random modeling method for hyperspectral alteration information of drill core in three-dimensional space |
-
2021
- 2021-12-22 CN CN202111634129.4A patent/CN114332391A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160024232A (en) * | 2014-08-25 | 2016-03-04 | 서울대학교산학협력단 | 3-Dimensional Space Modeling Method based on the Geotechnical Information |
CN106194181A (en) * | 2016-08-08 | 2016-12-07 | 西安科技大学 | Intelligent work surface coal-rock interface identification method based on geologic data |
CN110176070A (en) * | 2019-05-07 | 2019-08-27 | 中山大学 | A kind of bed succession analogy method carrying out machine learning based on borehole data |
CN111127630A (en) * | 2019-12-17 | 2020-05-08 | 核工业北京地质研究院 | Random modeling method for hyperspectral alteration information of drill core in three-dimensional space |
Non-Patent Citations (1)
Title |
---|
JINWU OUYANG ET AL: "Triangulated irregular network-based probabilistic 3D geological modelling using Markov Chain and Monte Carlo simulation", 《ENGINEERING GEOLOGY》, 2 May 2023 (2023-05-02) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116958470A (en) * | 2023-07-25 | 2023-10-27 | 中山大学 | Geological modeling method and device integrating Markov chain and multipoint statistics |
CN116958470B (en) * | 2023-07-25 | 2024-05-07 | 中山大学 | Geological modeling method and device integrating Markov chain and multipoint statistics |
CN117934767A (en) * | 2024-03-22 | 2024-04-26 | 华中科技大学 | Three-dimensional geological profile generation method and device and geological information modeling system |
CN117934767B (en) * | 2024-03-22 | 2024-06-11 | 华中科技大学 | Three-dimensional geological profile generation method and device and geological information modeling system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105651676B (en) | A kind of reservoir heterogeneity characterizing method under horizontal well rule well pattern | |
CN109116428B (en) | Fracture-cavity carbonate reservoir uncertainty modeling method and device | |
Che et al. | Three-dimensional geological modeling of coal seams using weighted Kriging method and multi-source data | |
US10853533B2 (en) | Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof | |
CN103077558A (en) | Method for establishing large-sized karst cave reservoir body distribution model of fractured-vuggy carbonate rock reservoir | |
CN110244021B (en) | Stratum layering method based on anisotropic interpolation | |
CN114332391A (en) | Three-dimensional geologic body probability model modeling method | |
CN106940450A (en) | Geosteering method for establishing model based on time-domain seismic volume | |
EP3455459A2 (en) | Three-dimensional fracture abundance evaluation of subsurface formations | |
CN117852416B (en) | Multimode grouting precontrolled analysis method and system based on digital geological model | |
SA519410882B1 (en) | A new porosity independent methodology for permeability prediction based on micro-resistivity images and laterolog resistivities | |
CN105242328A (en) | Paleo-thermal lithosphere thickness determination method and device | |
US10650107B2 (en) | Three-dimensional subsurface formation evaluation using projection-based area operations | |
CN109584365B (en) | Voronoi diagram geological layered model construction method for earth volume manufacturing cost accounting of ocean reclamation land-making engineering | |
CN109472867B (en) | Method for quantifying influence range of drilling position information | |
CN110309597A (en) | Method, apparatus and storage medium are determined based on the choked flow band model at configuration interface | |
CN113742962B (en) | Shale reservoir attribute three-dimensional modeling method based on horizontal well | |
CN113738343B (en) | VR underground drilling machine state detection method, system, device and equipment | |
Lyu et al. | Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China | |
CN111815769B (en) | Modeling method, computing device and storage medium for thrust covered zone construction | |
Erzeybek Balan | Characterization and modeling of paleokarst reservoirs using multiple-point statistics on a non-gridded basis | |
CN110687591A (en) | Method for determining physical property parameters of coal bed and surrounding rock based on waveform matching of prior data | |
CN115469361B (en) | Clastic rock stratum three-dimensional geological modeling method | |
Witty et al. | Application of Geostatistical Sequential Simulation Methods for Probabilistic 3D Subsoil Modeling and Uncertainty Quantification Concept and Examples | |
CN109541683B (en) | Construction method and system of underground point dam inclined stratum model grid |
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 |