NL2028480A - Stratum dynamic correction method of guidance while drilling - Google Patents

Stratum dynamic correction method of guidance while drilling Download PDF

Info

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
Application number
NL2028480A
Other languages
Dutch (nl)
Other versions
NL2028480B1 (en
Inventor
Wang Zhonghao
Liu Shaohua
Original Assignee
Univ Yangtze
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 Univ Yangtze filed Critical Univ Yangtze
Publication of NL2028480A publication Critical patent/NL2028480A/en
Application granted granted Critical
Publication of NL2028480B1 publication Critical patent/NL2028480B1/en

Links

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

STRATUM DYNAMIC CORRECTION METHOD OF GUIDANCE WHILE DRILLING TECHNICAL FIELD
[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.
BACKGROUND ART
[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.
SUMMARY
[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.
BRIEFT DESCRIPTION OF THE DRAWINGS
[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.
DETAILED DESCRIPTION OF THE EMBODIMENTS
[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)

CONCLUSIESCONCLUSIONS 1. Een methode voor het boren van stratum bestaande uit het toepassen van een stratumdynamische correctie van de begeleiding tijdens het boren, waarbij: de methode bestaat uit de volgende stappen: het verkrijgen van een expertmodel en gelaagde gegevenspunten van stratuminterpretatie; het toevoegen van controlepunten volgens de realtime evaluatiegegevens van het expertmodel; het volgens de gelaagde gegevenspunten en de controlepunten berekenen van hoogtes van rasters die door alle controlepunten worden beïnvloed; en het corrigeren van het expertmodel aan de hand van de hoogtes van de rasters die door alle controlepunten worden beïnvloed; en het aanpassen van het boorgattraject van de boorapparatuur volgens het gecorrigeerde model.A method for drilling stratum comprising applying a stratum dynamic correction of the guidance during drilling, wherein: the method consists of the following steps: obtaining an expert model and layered data points of stratum interpretation; adding checkpoints according to the real-time evaluation data of the expert model; calculating heights of grids affected by all control points according to the layered data points and the control points; and correcting the expert model based on the heights of the grids affected by all control points; and adjusting the borehole trajectory of the drilling equipment according to the corrected model. 2. De methode volgens conclusie 1, waarbij de stap voor het berekenen van de hoogte als volgt verloopt: het berekenen van de afstand van het huidig rasterpunt tot het gelaagd gegevenspunt en het controlepunt in de zone die door het controlepunt wordt beïnvloed; het volgens de afstand tussen het huidige rasterpunt, het gelaagde gegevenspunt en het controlepunt, bepalen van de hoogtegewichten van het gelaagd gegevenspunt en het controlepunt; het volgens de hoogte van het gelaagde gegevenspunt en het controlepunt en hun respectievelijke hoogtegewichten, verkrijgen van de hoogte van het huidig rasterpunt.The method of claim 1, wherein the step of calculating the height is as follows: calculating the distance from the current grid point to the layered data point and the checkpoint in the zone affected by the checkpoint; determining the height weights of the layered data point and the control point according to the distance between the current grid point, the layered data point and the control point; obtaining the height of the current grid point according to the height of the layered data point and the control point and their respective height weights. 3. De methode volgens conclusie 2, waarbij de afstand ten minste één van de euclidische afstand, Manhattan-afstand, Chebyshev-afstand, Hammingafstand en Minkowski-afstand is.The method of claim 2, wherein the distance is at least one of Euclidean distance, Manhattan distance, Chebyshev distance, Hamming distance and Minkowski distance. 4. De methode volgens conclusie 3, waarbij de afstand euclidische afstand is.The method of claim 3, wherein the distance is Euclidean distance. 5. De methode volgens conclusie 2 waarbij het hoogtegewicht omgekeerd evenredig is met de afstand van het huidige rasterpunt tot het gelaagde gegevenspunt en het controlepunt.The method of claim 2 wherein the height weight is inversely proportional to the distance from the current grid point to the layered data point and the control point. 6. De methode volgens conclusie 4, waarbij de berekeningsmethode van het hoogtegewicht als volgt verloopt: de Td waarbij pi het gewicht is van het i-de controlepunt tot het huidige rasterpunt, rastefrowjcot, di de afstand is van het i-de controlepunt tot het huidige rasterpunt, rasteTiowjcot, en dj de afstand is van het j-de controlepunt tot het huidige rasterpunt, rasterirow Jeol.The method of claim 4, wherein the height weight calculation method is as follows: the Td where pi is the weight from the ith control point to the current grid point, rastefrowjcot, ie the distance from the ith control point to the current grid point, rasteTiowjcot, and dj is the distance from the jth control point to the current grid point, rasterirow Jeol. 7. De methode volgens conclusie 6, waarbij de berekeningsmethode van de hoogtewaarde van het huidige rasterpunt als volgt verloopt: UR, Siro foal = > PES; js waarbij irow en jcol het rijnummer en het kolomnummer van het vorige raster zijn, Ziow jor de hoogte is van het huidige rasterpunt, rasteriowjco, px het overeenkomstige hoogtegewicht is van het controlepunt of het gelaagde gegevenspunt, en z; de hoogte is van het controlepunt of het gelaagde gegevenspunt.The method of claim 6, wherein the calculation method of the height value of the current grid point is as follows: UR, Sirofoal = > PES; js where irow and jcol are the row number and column number of the previous grid, Ziow jor is the height of the current grid point, rasteriowjco, px is the corresponding height weight of the control point or layered data point, and z; is the height of the control point or layered data point. 8. De methode volgens conclusie 2, waarbij de stap van het berekenen van de hoogte verder bestaat uit het stoppen van de berekening wanneer de hoogteverandering van het huidige rasterpunt kleiner is dan een drempelwaarde.The method of claim 2, wherein the step of calculating the height further comprises stopping the calculation when the height change of the current grid point is less than a threshold value. 9. De methode volgens conclusie 1 of 2, waarbij de drempelwaarde niet hoger is dan 0,01 m.The method of claim 1 or 2, wherein the threshold value is not greater than 0.01 m. 10. De methode volgens conclusie 1 of 2, waarbij de realtime evaluatiegegevens allemaal de hoogteverandering bevatten die is bepaald door feitelijk gemeten boorgatgegevens, profielgegevens, vlakgegevens en resultaatgegevens die door deskundigen zijn herkend.The method of claim 1 or 2, wherein the real-time evaluation data all include the elevation change determined by actually measured borehole data, profile data, face data, and result data recognized by those skilled in the art.
NL2028480A 2020-06-22 2021-06-18 Stratum dynamic correction method of guidance while drilling NL2028480B1 (en)

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)

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

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

Patent Citations (3)

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