Nothing Special   »   [go: up one dir, main page]

CN114332391A - Three-dimensional geologic body probability model modeling method - Google Patents

Three-dimensional geologic body probability model modeling method Download PDF

Info

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
Application number
CN202111634129.4A
Other languages
Chinese (zh)
Inventor
欧阳进武
周翠英
刘镇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202111634129.4A priority Critical patent/CN114332391A/en
Publication of CN114332391A publication Critical patent/CN114332391A/en
Pending legal-status Critical Current

Links

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

Three-dimensional geologic body probability model modeling method
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
Figure BDA0003427192630000021
Figure BDA0003427192630000022
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):
Figure BDA0003427192630000031
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).
Figure BDA0003427192630000032
Figure BDA0003427192630000033
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;
Figure BDA0003427192630000034
a theta step probability transfer matrix for the virtual borehole;
Figure BDA0003427192630000035
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.
Figure BDA0003427192630000041
(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.
CN202111634129.4A 2021-12-22 2021-12-22 Three-dimensional geologic body probability model modeling method Pending CN114332391A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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