NL2028480A - Stratum dynamic correction method of guidance while drilling - Google Patents
Stratum dynamic correction method of guidance while drilling Download PDFInfo
- Publication number
- NL2028480A NL2028480A NL2028480A NL2028480A NL2028480A NL 2028480 A NL2028480 A NL 2028480A NL 2028480 A NL2028480 A NL 2028480A NL 2028480 A NL2028480 A NL 2028480A NL 2028480 A NL2028480 A NL 2028480A
- Authority
- NL
- Netherlands
- Prior art keywords
- point
- distance
- height
- current grid
- data
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005553 drilling Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- 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
-
- 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/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
- G01V2200/10—Miscellaneous details
- G01V2200/16—Measure-while-drilling or logging-while-drilling
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/54—Borehole-related corrections
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/641—Continuity of geobodies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
- G01V2210/665—Subsurface modeling using geostatistical modeling
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present disclosure relates to a stratum dynamic correction method of guidance while drilling. The method comprises the steps of obtaining an expert model and layered data points of stratum interpretation; adding control points according to the real-time elevation data of the expert model; according to the layered data points and the control points, calculating elevations of grids influenced by all the control points, and correcting the expert model according to the elevations of the grids influenced by all the control points. According to the method, the models are selectively updated according to the grid points of the nearby influence areas where the control points are located, elevation values of all the grid models do not need to be reconstructed, the calculated amount is reduced, the updating speed of the models is greatly increased, and the real-time performance and precision of expert models are improved.
Description
[01] The present disclosure relates to the field of geological exploration, relates to a correction method of a geological model, in particular to a stratum dynamic correction method of guidance while drilling.
[02] In the interpretation process of real-time well logging while drilling, interpretation experts can obtain the expert model of stratum interpretation through the interpretation of the relationship between a borehole trajectory and a stratum, and make fine adjustments to the current stratum interface through this model. At this time, it is necessary to locally correct the initial structural model established before, so that the modeled stratum model is consistent with the actual stratum, providing a more accurate stratum interface for subsequent geological guidance and providing technical support for guidance.
[03] Elevation refers to the distance from a certain point to absolute datum in the direction of the plumb line, which is referred to as absolute elevation and is simplify referred to as elevation. The distance from a certain point to the datum of an assumed level is referred to as the assumed elevation.
[04] Horizon refers to a specific position in the stratigraphic succession. There are many stratum horizons, such as lithologic horizons with special lithologic characteristics, fossil horizons with special fossils, chronological horizons with specific era, seismic horizons, electrical logging horizons, etc. Therefore, the stratum horizon can be the boundary line of stratum units, and can also be a marker layer belonging to a specific era.
[05] Aiming at the technical problems existing in the prior art, the present disclosure provides a method for drilling stratum comprising applying a dynamic correction guidance while drilling, so as to improve the precision and real-time performance of an expert model.
[06] According to the present disclosure, the technical scheme for solving the technical problems is as follows. A method for drilling stratum comprising applying a stratum dynamic correction of guidance while drilling comprises the steps of: S1, obtaining an expert model and layered data points of stratum interpretation; S2, adding control points according to the real-time elevation data of the expert model, S3, according to the layered data points and the control points, calculating elevations of grids influenced by all the control points; and S4, correcting the expert model according to the elevations of the grids influenced by all the control points.
[07] In some embodiments of the present disclosure, the step of calculating the elevation in S3 is as follows: calculating the distance from the current grid point to the layered data point and the control point in the area influenced by the control point; according to the distance between the current grid point, the layered data point and the control point, determining the elevation weights of the layered data point and the control point; according to the elevation of the layered data point and the control point and their respective elevation weights, obtaining the elevation of the current grid point.
[08] Common distances include at least one of Euclidean distance, Manhattan distance, Chebyshev distance, Hamming distance and Minkowski distance. Further, the distance is Euclidean distance.
[09] In some embodiments of the present disclosure, the distance between the grid point and the layered data point and the control point calculates the elevation weight, such as Euclidean distance. The elevation weight is inversely proportional to the distance from the current grid point to the layered data point and the control point.
[10] In some of the above embodiments, the calculation method of the elevation weight is as follows: bos | dy jl a;
[11] :
[12] where pi is the weight of the i-th control point to the current grid point gridirow jet, di 18 the distance from the i-th control point to the current grid point gridiow jcot, and dj is the distance from the j-th control point to the current grid point gridiow jeol.
[13] In some embodiments of the determined elevation weight calculation method described above, the calculation method of the elevation value of the current grid point is as follows: Spore goad > > Pry
[14] Fd
[15] where irow and jcol are the row number and column number of the previous grid, Ziow jcot ís the elevation of the current grid point gridicow jeot, pris the corresponding elevation weight of the control point or the layered data point, and z; is the elevation of the control point or the layered data point.
[16] In some embodiments of the present disclosure, the step of calculating the elevation further comprises stopping the calculation when the elevation change of the current grid point is less than a threshold value. Preferably, the threshold value is not more than 0.01m.
[17] In the above embodiments, the real-time elevation data all include the elevation change determined by actually measured borehole data, profile data, plane data and result data recognized by experts.
[18] The method has the beneficial effects that the models are selectively updated according to the grid points of the nearby influence areas where the control points are located, elevation values of all the grid models do not need to be reconstructed, and the updating speed of the models is greatly increased.
[19] FIG. 1 is a basic flow chart of a stratum dynamic correction method of a geological model according to some embodiments of the present disclosure;
[20] FIG. 2 is a schematic diagram of a three-dimensional model of a stratum according to some embodiments of the present disclosure;
[21] FIG. 3 is a schematic diagram of an adjusted three-dimensional model of a stratum according to some embodiments of the present disclosure.
[22] The principles and features of the present disclosure will be described hereinafter with reference to the accompanying drawings. The examples are only used to explain the present disclosure, rather than limit the scope of the present disclosure.
[23] The present disclosure relates to a stratum dynamic correction method of guidance while drilling, wherein: the method comprises the steps of: obtaining an expert model and layered data points of stratum interpretation; adding control points according to the real-time elevation data of the expert model; according to the layered data points and the control points, calculating elevations of grids influenced by all the control points; and correcting the expert model according to the elevations of the grids influenced by all the control points.
[24] In some embodiments of the present disclosure, the step of calculating the elevation is as follows: calculating the distance from the current grid point to the layered data point and the control point in the area influenced by the control point; according to the distance between the current grid point, the layered data point and the control point, determining the elevation weights of the layered data point and the control point; according to the elevation of the layered data point and the control point and their respective elevation weights, obtaining the elevation of the current grid point.
[25] Common distances include Euclidean distance, Manhattan distance, Chebyshev distance, Hamming distance and Minkowski distance. In order to improve the model precision, one or more of the distances can be selected for combined calculation. Preferably, the distance from the current grid point to the layered data point and the control point is Euclidean distance.
[26] In some embodiments of the present disclosure, the elevation weight is inversely proportional to the distance from the current grid point to the layered data point and the control point. Specifically, the calculation method of the elevation weight is as follows: oS]
[27] © formula (1)
[28] where pi is the weight of the i-th control point to the current grid point gridirowjcot, di is the distance from the i-th control point to the current grid point 5 gridiowjcot, and dj is the distance from the j-th control point to the current grid point gridisow col.
[29] In some embodiments of the determined elevation weight calculation method described above, the calculation method of the elevation value of the current grid point is as follows: 3 i. Sirew geel > > Pry 0 [30] Fed formula (2)
[31] where irow and jcol are the row number and column number of the previous grid, Zirow jcot Is the elevation of the current grid point gridicow,jcot, Px is the corresponding elevation weight of the control point or the layered data point, and z; is the elevation of the control point or the layered data point.
[32] In some embodiments of the present disclosure, a stratum dynamic correction method of guidance while drilling comprises the following specific steps:
[33] (1) adding three-dimensional control points in the target horizon;
[34] (2) calculating the distance between the current grid point gridiow jcot in the area influenced by the control point and the layered data point of the target layer;
[35] (3) calculating the distance between the current grid point gridiow jc: and the control point;
[36] (4) calculating the weight according to the above formula 1;
[37] (5) calculating the elevation value of the current grid point gridiow jest according to the above formula 2;
[38] (0) circulating all grids in the area influenced by control points, and repeating the above steps 2-5 to complete the adjustment of the model.
[39] In some embodiments of the present disclosure, the step of calculating the elevation further comprises stopping the calculation when the elevation change of the current grid point is less than a threshold value. Preferably, the threshold value is not more than 0.01m.
[40] Some embodiments of the present disclosure will be explained with reference to the attached drawings. The red dots P1 to P8 in FIG. 2 are layered data control points. interpretation experts can obtain an expert model of stratum interpretation by interpreting the relationship between a borehole trajectory and a stratum, and add a new control point P9 from elevation data. When the point pointed by the arrow in FIG. 3 is the layered data point which is the added stratum control point, the model needs to be adjusted, but it does not need to be re-modeled, and only the elevation of the grid points near P9 needs to be recalculated.
[41] According to the coordinates of point P9, the row and column number where it is located are calculated, diffusing around the grid, and increasing outward step by step at a time. The elevation of the current grid point gridiowjcot is calculated. First, the distance between the grid point gridiow jot and points P1 to P9 is calculated, the weight pi is calculated by formula (1), and then the elevation of the current grid point gridirow jcot is calculated by formula (2). After outward diffusion, the change (absolute value of the difference) between the calculated elevation of the grid point and the original elevation 1s less than a threshold value, and the threshold value is not more than 0.01m. The outward diffusion is stopped. At this time, the model adjustment of the area influenced by the new control point P9 ends, and the whole model completes updating. The updated model is shown in FIG. 3.
[42] In the above embodiments, the real-time elevation data all include the elevation change determined by actually measured borehole data, profile data, plane data and result data recognized by experts.
[43] The above is only a preferred embodiment of the present disclosure, and is not intended to limit the present disclosure. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present disclosure shall be included in the scope of protection of the present disclosure.
Claims (10)
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010576436.0A CN111784832B (en) | 2020-06-22 | 2020-06-22 | Stratum dynamic correction method of geological model |
Publications (2)
Publication Number | Publication Date |
---|---|
NL2028480A true NL2028480A (en) | 2021-09-16 |
NL2028480B1 NL2028480B1 (en) | 2022-06-03 |
Family
ID=72757667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NL2028480A NL2028480B1 (en) | 2020-06-22 | 2021-06-18 | Stratum dynamic correction method of guidance while drilling |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN111784832B (en) |
NL (1) | NL2028480B1 (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150066460A1 (en) * | 2013-08-30 | 2015-03-05 | Jimmy Klinger | Stratigraphic function |
CN105317375B (en) * | 2015-10-15 | 2017-10-17 | 中国石油天然气集团公司 | Horizontal well is inducted into Target process and device |
CN110685600A (en) * | 2018-06-20 | 2020-01-14 | 中国石油化工股份有限公司 | Drill bit adjustment prediction method for geosteering |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5838634A (en) * | 1996-04-04 | 1998-11-17 | Exxon Production Research Company | Method of generating 3-D geologic models incorporating geologic and geophysical constraints |
NO336002B1 (en) * | 2013-09-17 | 2015-04-20 | Roxar Software Solutions As | Method for providing a geological model based on measured geological data |
CN106894383A (en) * | 2017-03-10 | 2017-06-27 | 中国电建集团成都勘测设计研究院有限公司 | A kind of Three-Dimensional Dynamic method for designing based on building stones exploitation and sub-area utilization |
CN107180452A (en) * | 2017-05-02 | 2017-09-19 | 武汉中地数码科技有限公司 | The stratified geological objects modeling method of the multi-source data coupling on the stratum containing standard |
CN111311746B (en) * | 2019-12-09 | 2023-09-08 | 中交广州航道局有限公司 | Intelligent three-dimensional geological modeling method based on drilling data |
-
2020
- 2020-06-22 CN CN202010576436.0A patent/CN111784832B/en active Active
-
2021
- 2021-06-18 NL NL2028480A patent/NL2028480B1/en active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150066460A1 (en) * | 2013-08-30 | 2015-03-05 | Jimmy Klinger | Stratigraphic function |
CN105317375B (en) * | 2015-10-15 | 2017-10-17 | 中国石油天然气集团公司 | Horizontal well is inducted into Target process and device |
CN110685600A (en) * | 2018-06-20 | 2020-01-14 | 中国石油化工股份有限公司 | Drill bit adjustment prediction method for geosteering |
Also Published As
Publication number | Publication date |
---|---|
CN111784832A (en) | 2020-10-16 |
NL2028480B1 (en) | 2022-06-03 |
CN111784832B (en) | 2024-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105317375B (en) | Horizontal well is inducted into Target process and device | |
CN106194181B (en) | Intelligent working face coal-rock interface recognition method based on geological data | |
CN105807316B (en) | Ground observation microseism velocity model corrections method based on amplitude superposition | |
CN107132578B (en) | A kind of microseism ground monitoring velocity model corrections algorithm | |
CN107843928A (en) | A kind of shot point method for correcting error based on grid-search algorithms | |
CN110073246B (en) | Improved method relating to quality control | |
CN105093319B (en) | Ground micro-seismic static correcting method based on 3D seismic data | |
CN105631753A (en) | Attitude of stratum-based horizontal well oil reservoir profile modeling method | |
KR102155101B1 (en) | System and method for predicting ground layer information, and a recording medium having computer readable program for executing the method | |
CN108333627A (en) | Igneous rock area is broken recognition methods and the device of the true and false | |
CN105911592B (en) | A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint | |
CN106054252B (en) | A kind of method and device of pre-stack time migration | |
NL2028480B1 (en) | Stratum dynamic correction method of guidance while drilling | |
CN114707253A (en) | Feature-guidance-based three-dimensional measurement robot scanning track optimization method | |
CN113236249B (en) | Self-adaptive cutting method for fusing geological baselines | |
WO2014149561A1 (en) | Migration velocity analysis method for vertical seismic profile data | |
CN107818217B (en) | Seismic guiding method based on horizontal well exploration type visual analysis | |
CN111999768B (en) | Construction elevation determination method based on construction model in drilling layered design | |
Yilmaz et al. | Application of artificial neural networks to height transformation | |
KR102232616B1 (en) | Method of tracking position of mobile object in space and apparatus of tracking using the method | |
WO2021260563A3 (en) | System and method for subsurface structural interpretation | |
GB2615216A (en) | Utilizing micro-services for optimization workflows of borehole operations | |
CN111577250A (en) | Horizontal well geosteering method and device | |
CN117420569B (en) | Inversion method of non-uniform wind field based on Doppler laser wind finding radar | |
CN107862100A (en) | A kind of shot point method for correcting error based on particle cluster algorithm |