CN113482596A - LWD-based while-drilling instrument real-time three-dimensional gamma imaging data processing method - Google Patents
LWD-based while-drilling instrument real-time three-dimensional gamma imaging data processing method Download PDFInfo
- Publication number
- CN113482596A CN113482596A CN202110919136.2A CN202110919136A CN113482596A CN 113482596 A CN113482596 A CN 113482596A CN 202110919136 A CN202110919136 A CN 202110919136A CN 113482596 A CN113482596 A CN 113482596A
- Authority
- CN
- China
- Prior art keywords
- depthmwd
- azm
- gamma
- tvd
- depthgr
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 29
- 238000005553 drilling Methods 0.000 title claims abstract description 24
- 238000003672 processing method Methods 0.000 title claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 16
- 238000009877 rendering Methods 0.000 claims abstract description 8
- 238000005259 measurement Methods 0.000 claims description 41
- 238000000034 method Methods 0.000 claims description 13
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 7
- 210000002445 nipple Anatomy 0.000 claims description 7
- 230000005251 gamma ray Effects 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 238000005192 partition Methods 0.000 claims description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 4
- 238000011161 development Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 239000003345 natural gas Substances 0.000 abstract description 2
- 238000005755 formation reaction Methods 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 6
- 230000009471 action Effects 0.000 description 3
- 239000007789 gas Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- 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
- E21B47/00—Survey of boreholes or wells
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/02—Determining slope or direction
- E21B47/022—Determining slope or direction of the borehole, e.g. using geomagnetism
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/04—Measuring depth or liquid level
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Image Processing (AREA)
Abstract
The invention provides a real-time three-dimensional gamma imaging data processing method based on an LWD (logging while drilling) instrument, which relates to the technical field of oil and natural gas exploration and development and comprises the following steps of firstly carrying out depth completion on well deviation and azimuth data (depthmwd, inc, azm) and gamma data (depthgr, GR) of all measuring points, secondly carrying out equidistant D interpolation on the depthmwd, and then calculating each depthmwd point (depthmwd)d_m，tvdd_m，ewd_m，nsd_m) Finally, calculating the updated data and rendering images to finish real-time three-dimensional imaging, wherein the real-time three-dimensional gamma imaging is based on the LWD while-drilling instrumentThe data processing method can enable the formation gamma characteristics measured by the orientation logging while drilling instrument in the LWD to be represented in a three-dimensional mode, and is beneficial to a field geosteering engineer to observe the formation gamma characteristics more intuitively and make an accurate decision for geosteering.
Description
Technical Field
The invention relates to the technical field of petroleum and natural gas exploration and development, in particular to a real-time three-dimensional gamma imaging data processing method based on an LWD (logging while drilling) instrument.
Background
In recent years, with the arrival of the times of offshore and unconventional oil and gas resource development, the technology of large-pitch and horizontal wells is also developed vigorously, wherein the Logging While Drilling technology (LWD Logging While Drilling) is not available. LWD is a Measurement While Drilling (MWD) technology developed on the basis of MWD measurement While Drilling (MWD measurement While Drilling) technology, and the LWD technology has more complete functions and more complex structure. The method mainly comprises the items of MWD, azimuth natural gamma logging, neutron density logging and the like.
The azimuth gamma logging is one of essential items of a logging-while-drilling technology in petroleum and gas exploration, and real-time three-dimensional gamma imaging can be realized by combining drill bit track information measured by MWD. The real-time three-dimensional gamma imaging can visually observe the change of the well wall and the conditions of the stratum encountered by drilling when the well is drilled around the well body, thereby effectively controlling the landing and the trend of the well track. Historically, the drilling trajectory, represented in the form of a "line", has not modeled three-dimensionally the "wellbore", and has been less able to achieve a roaming view of the formation being drilled in the wellbore, and therefore, has been far from the actual requirements of geosteering while drilling.
Currently, among LWD (gamma-containing) instruments used in various large oil fields, mwd (measurement While drilling) instruments measure information on the orientation of a drill bit, and azimuthal gamma logging instruments (nipples) measure gamma characteristics of formations. Generally, the azimuth gamma measuring nipple is closer to the drill bit than the well deviation and the azimuth MWD nipple of the measuring drill bit, and the difference between the azimuth gamma measuring nipple and the measurement drill bit is about 8 meters, namely, a blind area exists. In addition, because the frequency of the deviation and the azimuth of the LWD measurement is lower, an inflection point can appear when the real-time drawn borehole trajectory is drawn, namely the borehole is not slippery.
Disclosure of Invention
Technical problem to be solved
Aiming at the defects of the prior art, the invention provides a real-time three-dimensional gamma imaging data processing method based on an LWD (logging while drilling) instrument, and solves the problems that an azimuth gamma measuring pup joint is closer to a drill bit than a well deviation of a measuring drill bit and an azimuth MWD pup joint, the difference between the azimuth gamma measuring pup joint and the measurement drill bit is about 8 meters, a blind area exists, an inflection point appears when a real-time drawn well track, and a well is not smooth.
(II) technical scheme
In order to achieve the purpose, the invention is realized by the following technical scheme: in LWD actual measurement, the data of real-time three-dimensional gamma imaging comprises well depth (md), well deviation (inc), azimuth (azm) and GRi(GR denotes the formation gamma characteristics measured by the azimuth gamma nipple) — 1, 2,. n, n denotes the number of sectors of the azimuth gamma logging maximum measurement, n is typically 2, 4, 8, 16;
where md denotes the depth of the well, and when the tool is operated downhole, the depth at which the tool is located can be considered md; inc represents the angle of the well, azm represents the orientation of the well, typically denoted (depthmwd, inc, azm), where depthmwd-md-zero length (well deviation, orientation); GR is expressed as the formation gamma characteristic measured for azimuthal gamma sub, where the data is expressed as (depthgr, GR)i) Wherein depthgr is md-zero length (GR); the LWD measurement well deviation and azimuth frequency is much less than the formation gamma characteristic frequency;
before three-dimensional gamma imaging is achieved, the field engineer will give joint points, mainly including depthmwdti、incti、azmti、tvdti、nstiChinese and ancient ewtiSix parameters, among which depthmwdtiWell depth, inc, denoted as the junction pointtiDenoted as the junction point well deviation, azmtiIndicates the orientation of the joint points, tvdtiDenotes the vertical depth of the joint point, nstiIndicates the north-south displacement of the joint pointtiRepresenting east-west displacement of joint points;
the method comprises the following steps:
step 1, in LWD instrument, the sub for azimuth gamma measurement is closer to the drill bit than the sub for well deviation and azimuth measurement, so that the current measurement (depthgr, GR)i) Greater than (depthmwd, inc, azm), i.e. depthgr is greater than depthmwd; in order to measure (depthgr-depthmwd, GRi) The gamma value of the range is matched with well deviation and direction and needs to be carried outLine (depthgr-depthmwd, inc, azm) completion, well deviation and azimuth data format (depthmwd) for all measurement pointsti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgr，inck，amzk) (ii) a In the formula (depthmwd)ti，incti，azmti) This junction point, (depthgr)last，inck，amzk) Gamma-ray fill-in point, depthgrrlastThe depth of the last measurement gamma is k, and the number of the measurement point is k;
step 2, performing equidistant D interpolation on depthmwd;
step 2.1, according to (depthmwd)ti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgrrlast，inck，amzk) Using a natural parameter interpolation model (the variation of the well deviation change rate and the azimuth change rate of interpolation is kept constant), performing equidistant D interpolation on depthmwd, wherein D is 0.1m, and performing depthmwdtiAnd depthgrlastOne decimal fraction is reserved, and the partition interval is several equal divisions of 0.1m, denoted depthmwdd_1，depthmwdd_2，……，depthmwdd_mWherein m represents the number of sequence numbers;
step 2.2, calculate depthmwdd_1，depthmwdd_2，……，depthmwdd_mInc and azm, in particular depthmwdd_mIn (depthmwd)ti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgrrlast，inck，amzk) Two adjacent depth points are (depthmwd)n，incn，azmn)，(depthmwdn+1，inc n+1，azm n+1) Then depthmwdd_mWell deviation inc of pointsd_m，azmd_mIn order to realize the purpose,
wherein depthmwdd_mWell deviation inc of pointsd_m，azmd_mFormula, all points (depthmwd) can be obtainedd_1，incd_1，azmd_1)，(depthmwdd_2，incd_2，azmd_2)，……，(depthmwdd_m，incd_m，azmd_m)；
Step 3, calculating each measuring point (depthmwd)d_m，tvdd_m，nsd_m，ewd_m)；
Step 3.1, data (depthmwd) calculated in step 2d_1，incd_1，azmd_1)，(depthmwdd_2，incd_2，azmd_2)，……，(depthmwdd_m，incd_m，azmd_m) Tvd, ns, ew for each point can be calculated, and the currently measured point tvd2, ns2 and ew2 can be calculated from the (depthmwd2, inc2, azm2) measured by the instrument at the current point, as follows,
wherein:
γ＝arccos[cos(inc1)cos(inc2)+sin(inc1)sin(inc2）cos(azm2-azm1)]
ΔL＝depthmwd2-depthmwd1
in the above formula, depthmwd1, inc1, azm1, tvd1, ns1 and ew1 are respectively expressed as the depth, inclination, azimuth, dip, north-south displacement and east-west displacement of the last measured point;
step 3.2, the following data (depthmwd) can be obtained through cyclic calculationti，tvdti，nsti，ewti)，(depthmwdd_1，tvdd_1，nsd_1，ewd_1)，(depthmwdd_2，tvdd_2，nsd_2，ewd_2)，......，(depthmwdd_m，tvdd_m，nsd_m，ewd_m)；
Step 4, rendering and calculating the updated data and images to finish real-time imaging;
step 4.1, in order to make the three-dimensional gamma imaging have a real-time function, the data volume of one-time rendering should not be too large, so the calculated (depthmwd, tvd, ns, ew) and (depthgr, GR) need to be comparedi) Matching is carried out, and only updated data are rendered each time;
step 4.2, respectively taking out the latest measured depth depthmwd of the well deviation azimuthkAnd depth depthgr of the latest measurement point of the formation gamma GRlastAnd is in (depthmwd)ti，tvdti，nsti，ewti)，(depthmwdd_1，tvdd_1，nsd_1，ewd_1)，(depthmwdd_2，tvdd_2，nsd_2，ewd_2)，......，(depthmwdd_m，tvdd_m，nsd_m，ewd_m) And (depthgr, GR)i) Find out the depth of [ depthmwdk，depthgrlast]Data within the interval, the portion of data being used to render the image.
(III) advantageous effects
The invention provides a real-time three-dimensional gamma imaging data processing method based on an LWD (logging while drilling) instrument. The method has the following beneficial effects:
the real-time three-dimensional gamma imaging data processing method based on the LWD while-drilling instrument solves the problem of a distance blind zone between an MWD (measurement while drilling) and a gamma short section in LWD measured data, and compared with the previous method, a field engineer can make accurate judgment on the position of the blind zone. In addition, the drawn real-time three-dimensional borehole trajectory has no inflection point, and the actual shape of the borehole is more truly restored.
Drawings
FIG. 1 is a three-dimensional gamma borehole imaging of measured data after processing of raw data according to the present invention.
FIG. 2 is a three-dimensional gamma borehole imaging plot of raw survey data.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
An embodiment of the present invention provides a real-time three-dimensional gamma imaging data processing method based on an LWD while-drilling instrument, as shown in fig. 1-2, in an LWD actual measurement, data of real-time three-dimensional gamma imaging includes well depth (md), well deviation (inc), azimuth (azm), and GRi(GR denotes the formation gamma characteristics measured by the azimuth gamma nipple) — 1, 2,. n, n denotes the number of sectors of the azimuth gamma logging maximum measurement, n is typically 2, 4, 8, 16;
where md denotes the depth of the well, and when the tool is operated downhole, the depth at which the tool is located can be considered md; inc represents the angle of the well, azm represents the orientation of the well, typically denoted (depthmwd, inc, azm), where depthmwd-md-zero length (well deviation, orientation); GR is expressed as the formation gamma characteristic measured for azimuthal gamma sub, where the data is expressed as (depthgr, GR)i) Whereindepthgr-md-zero length (GR); the LWD measurement well deviation and azimuth frequency is much less than the formation gamma characteristic frequency;
before three-dimensional gamma imaging is achieved, the field engineer will give joint points, mainly including depthmwdti、incti、azmti、tvdti、nstiChinese and ancient ewtiSix parameters, among which depthmwdtiWell depth, inc, denoted as the junction pointtiDenoted as the junction point well deviation, azmtiIndicates the orientation of the joint points, tvdtiDenotes the vertical depth of the joint point, nstiIndicates the north-south displacement of the joint pointtiRepresenting east-west displacement of joint points;
the method comprises the following steps:
step 1, in LWD instrument, the sub for azimuth gamma measurement is closer to the drill bit than the sub for well deviation and azimuth measurement, so that the current measurement (depthgr, GR)i) Greater than (depthmwd, inc, azm), i.e. depthgr is greater than depthmwd; in order to measure (depthgr-depthmwd, GRi) The gamma values of the range are matched with well deviation and azimuth, the (depthgr-depthmwd, inc, azm) is required to be supplemented, and the data format (depthmwd) of the well deviation and the azimuth of all the measuring points is adoptedti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgr，inck，amzk) (ii) a In the formula (depthmwd)ti，incti，azmti) This junction point, (depthgr)last，inck，amzk) Gamma-ray fill-in point, depthgrrlastThe depth of the last measurement gamma is k, and the number of the measurement point is k;
step 2, performing equidistant D interpolation on depthmwd;
step 2.1, according to (depthmwd)ti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgrrlast，inck，amzk) Using a natural parameter interpolation model (the variation of the well deviation change rate and the azimuth change rate of interpolation is kept constant), performing equidistant D interpolation on depthmwd, wherein D is 0.1m, and performing depthmwdtiAnd depthgrlastOne decimal fraction is reserved, and the partition interval is several equal divisions of 0.1m, denoted depthmwdd_1，depthmwdd_2，……，depthmwdd_mWherein m represents the number of sequence numbers;
step 2.2, calculate depthmwdd_1，depthmwdd_2，……，depthmwdd_mInc and azm, in particular depthmwdd_mIn (depthmwd)ti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgrrlast，inck，amzk) Two adjacent depth points are (depthmwd)n，incn，azmn)，(depthmwdn+1，inc n+1，azm n+1) Then depthmwdd_mWell deviation inc of pointsd_m，azmd_mIn order to realize the purpose,
wherein depthmwdd_mWell deviation inc of pointsd_m，azmd_mFormula, all points (depthmwd) can be obtainedd_1，incd_1，azmd_1)，(depthmwdd_2，incd_2，azmd_2)，……，(depthmwdd_m，incd_m，azmd_m)；
Step 3, calculating each measuring point (dept)hmwdd_m，tvdd_m，nsd_m，ewd_m)；
Step 3.1, data (depthmwd) calculated in step 2d_1，incd_1，azmd_1)，(depthmwdd_2，incd_2，azmd_2)，……，(depthmwdd_m，incd_m，azmd_m) Tvd, ns, ew for each point can be calculated, and the currently measured point tvd2, ns2 and ew2 can be calculated from the (depthmwd2, inc2, azm2) measured by the instrument at the current point, as follows,
wherein:
r＝arccos[cos(inc1)cos(inc2)+sin(inc1)sin(inc2)cos(azm2-azm1)]
ΔL＝depthmwd2-depthmwd1
in the above formula, depthmwd1, inc1, azm1, tvd1, ns1 and ew1 are respectively expressed as the depth, inclination, azimuth, dip, north-south displacement and east-west displacement of the last measured point;
step 3.2, the following data (depthmwd) can be obtained through cyclic calculationti，tvdti，nsti，ewti)，(depthmwdd_1，tvdd_1，nsd_1，ewd_1)，(depthmwdd_2，tvdd_2，nsd_2，ewd_2)，......，(depthmwdd_m，tvdd_m，nsd_m，ewd_m)；
Step 4, rendering and calculating the updated data and images to finish real-time imaging;
step 4.1, in order to make the three-dimensional gamma imaging have a real-time function, the data volume of one-time rendering should not be too large, so the calculated (depthmwd, tvd, ns, ew) and (depthgr, GR) need to be comparedi) Matching is carried out, and only updated data are rendered each time;
step 4.2, respectively taking out the latest measured depth depthmwd of the well deviation azimuthkAnd depth depthgr of the latest measurement point of the formation gamma GRlastAnd is in (depthmwd)ti，tvdti，nsti，ewti)，(depthmwdd_1，tvdd_1，nsd_1，ewd_1)，(depthmwdd_2，tvdd_2，nsd_2，ewd_2)，......，(depthmwdd_m，tvdd_m，nsd_m，ewd_m) And (depthgr, GR)i) Find out the depth of [ depthmwdk，depthgrlast]Data within the interval, the portion of data being used to render the image.
According to the method, the LWD measurement while drilling real-time data is calculated and processed, and finally real-time imaging is carried out. The data processing method of the invention effectively eliminates the data blind area and also improves the smoothness of drawing the well track curve.
TABLE 1 survey points well depth, well deviation, and azimuth data
| Serial number | Depth of field | Well deviation | Orientation | Serial number | Depth of field | Well deviation | Orientation |
| 1 | 3143.99 | 80.45 | 12.08 | 24 | 3794.01 | 82.02 | 7.6 |
| 2 | 3172.69 | 80.6 | 12.08 | 24 | 3818.9 | 82.77 | 7.69 |
| 3 | 3199.82 | 80.7 | 13.3 | 24 | 3851.62 | 83.47 | 8.13 |
| 4 | 3228.86 | 79.7 | 13.5 | 24 | 3880.65 | 84 | 8.48 |
| 5 | 3257.47 | 81.7 | 10.6 | 24 | 3908.8 | 82.86 | 8.83 |
| 6 | 3287.4 | 79.3 | 10.8 | 24 | 3936.45 | 82.46 | 8.48 |
| 7 | 3316 | 75.69 | 10.8 | 24 | 3967.62 | 81.67 | 8.74 |
| 8 | 3344.41 | 76.5 | 10.8 | 24 | 3996.3 | 81.49 | 8.92 |
| 9 | 3373.28 | 78.9 | 9.5 | 24 | 4024.91 | 82.33 | 9.45 |
| 10 | 3402.05 | 79.2 | 10.3 | 24 | 4052.08 | 80.97 | 9.27 |
| 11 | 3402.05 | 79.2 | 10.3 | 24 | 4081.62 | 80.35 | 8.74 |
| 12 | 3431.78 | 80.18 | 11.13 | 24 | 4110.41 | 80.84 | 10.06 |
| 13 | 3459.01 | 84.09 | 10.69 | 24 | 4138.5 | 82.11 | 10.06 |
| 14 | 3488.46 | 82.46 | 10.6 | 24 | 4165.91 | 84.53 | 9.45 |
| 15 | 3517.01 | 82.2 | 9.46 | 24 | 4195.54 | 85.05 | 9.97 |
| 16 | 3544.84 | 82.24 | 9.28 | 24 | 4225.39 | 82.99 | 9.8 |
| 17 | 3573.63 | 82.02 | 7.52 | 24 | 4253.85 | 83.16 | 10.32 |
| 18 | 3602.28 | 83.25 | 9.87 | 24 | 4282.42 | 83.74 | 10.59 |
| 19 | 3632.18 | 82.99 | 9.52 | 24 | 4310.46 | 84.53 | 10.94 |
| 20 | 3660.46 | 81.63 | 7.69 | 24 | 4338.86 | 83.82 | 11.03 |
| 21 | 3689.92 | 81.1 | 7.25 | 24 | 4366.96 | 82.33 | 10.85 |
| 22 | 3736.26 | 81.93 | 7.6 | 24 | 4425.36 | 81.32 | 10.76 |
| 23 | 3765.48 | 82.11 | 7.69 |
TABLE 2 measurement of azimuthal gamma data
Data of tables 1 and 2 are from data of LWD at a certain time in actual operation, from the raw measurements of Table 1 above
| Serial number | Depth of field | Upper gamma ray | Lower gamma | Left gamma | Right gamma |
| 1 | 3158.2 | 77.94 | 76.81 | 72.65 | 73.97 |
| 2 | 3158.3 | 76.77 | 75.39 | 76.58 | 72.06 |
| 3 | 3158.4 | 79.69 | 74.04 | 75.81 | 72.75 |
| 4 | 3158.5 | 84.36 | 72.73 | 73.01 | 74.55 |
| 5 | 3158.6 | 89.03 | 71.41 | 70.22 | 76.35 |
| 6 | 3158.7 | 93.7 | 70.09 | 67.42 | 78.15 |
| 7 | 3158.8 | 96.7 | 72.59 | 68.37 | 76.35 |
| 8 | 3158.9 | 96.54 | 77.78 | 70.4 | 73.59 |
| 9 | 3159 | 95.33 | 83.43 | 72.21 | 71.07 |
| 10 | 3159.1 | 94.11 | 89.08 | 74.01 | 68.55 |
| …… | …… | …… | …… | …… | …… |
| 12731 | 4431.3 | 78.8 | 82.12 | 81.89 | 76.17 |
| 12731 | 4431.4 | 78.88 | 82.18 | 82.27 | 76.4 |
| 12731 | 4431.5 | 78.96 | 82.24 | 82.66 | 76.63 |
| 12731 | 4431.6 | 79.03 | 82.3 | 83.04 | 76.86 |
| 12731 | 4431.7 | 79.11 | 82.36 | 83.42 | 77.09 |
| 12731 | 4431.8 | 79.19 | 82.43 | 83.8 | 77.32 |
| 12731 | 4431.9 | 79.26 | 82.49 | 84.18 | 77.55 |
| 12731 | 4432 | 79.34 | 82.55 | 84.56 | 77.78 |
| 12731 | 4432.1 | 79.42 | 82.61 | 84.94 | 78.01 |
| 12731 | 4432.2 | 79.49 | 82.67 | 85.32 | 78.24 |
| 44323 | 4432.3 | 79.57 | 82.73 | 85.7 | 78.48 |
| 44324 | 4432.4 | 79.65 | 82.8 | 86.08 | 78.71 |
The data shows that the interval of depth sampling points of the original data is large, the frequency of well deviation and azimuth measured by LWD is low, so that the phenomenon that an inflection point appears when a well track is drawn in real time, namely the well track is not smooth, and due to the structure of an instrument, a dead zone exists in the measurement.
As can be seen from comparison between FIG. 1 and FIG. 2, the method for processing the LWD measurement-while-drilling raw data in real time effectively eliminates the inflection point of the well trajectory, smoothes the well trajectory, and covers the blind zone existing during measurement. Therefore, the method provided by the invention can meet the requirement of effective application in engineering.
The working principle is as follows: when in use, the data (depthgr-depthmwd, inc, azm) of the well deviation and the azimuth of all the measuring points are firstly supplemented, then the equal-interval D interpolation is carried out on the depthmwd, and then each measuring point (depthmwd) is calculatedd_m，tvdd_m，nsd_m，ewd_m) And finally, calculating the updated data and rendering images to finish real-time three-dimensional imaging.
It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (1)
1. The real-time three-dimensional gamma imaging data processing method based on the LWD while-drilling instrument is characterized in that in LWD actual measurement, the real-time three-dimensional gamma imaging data comprise well depth (md), well deviation (inc), azimuth (azm) and GRi(GR denotes the formation gamma characteristics measured by the azimuth gamma nipple) — 1, 2,. n, n denotes the number of sectors of the azimuth gamma logging maximum measurement, n is typically 2, 4, 8, 16;
where md denotes the depth of the well, the tool may be considered to be located when it is operated downholeThe depth of (a) is md; inc represents the angle of the well, azm represents the orientation of the well, typically denoted (depthmwd, inc, azm), where depthmwd-md-zero length (well deviation, orientation); GR is expressed as the formation gamma characteristic measured for azimuthal gamma sub, where the data is expressed as (depthgr, GR)i) Wherein depthgr is md-zero length (GR); the LWD measurement well deviation and azimuth frequency is much less than the formation gamma characteristic frequency;
before three-dimensional gamma imaging is achieved, the field engineer will give joint points, mainly including depthmwdti、incti、azmti、tvdti、nstiChinese and ancient ewtiSix parameters, among which depthmwdtiWell depth, inc, denoted as the junction pointtiDenoted as the junction point well deviation, azmtiIndicates the orientation of the joint points, tvdtiDenotes the vertical depth of the joint point, nstiIndicates the north-south displacement of the joint pointtiRepresenting east-west displacement of joint points;
the method comprises the following steps:
step 1, in LWD instrument, the sub for azimuth gamma measurement is closer to the drill bit than the sub for well deviation and azimuth measurement, so that the current measurement (depthgr, GR)i) Greater than (depthmwd, inc, azm), i.e. depthgr is greater than depthmwd; in order to measure (depthgr, GRi) The gamma values of the range are matched with well deviation and azimuth, the (depthmwd, inc, azm) completion is needed, and the data format (depthmwd) of the well deviation and the azimuth of all the measuring points is adoptedti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgr，inck，amzk) (ii) a In the formula (depthmwd)ti，incti，azmti) This junction point, (depthgr)last，inck，amzk) Gamma-ray fill-in point, depthgrrlastThe depth of the last measurement gamma is k, and the number of the measurement point is k;
step 2, performing equidistant D interpolation on depthmwd;
step 2.1, according to (depthmwd)ti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgrrlast，inck，amzk) Using a natural parameter interpolation model (the variation of the well deviation change rate and the azimuth change rate of interpolation is kept constant), performing equidistant D interpolation on depthmwd, wherein D is 0.1m, and performing depthmwdtiAnd depthgrlastOne decimal fraction is reserved, and the partition interval is several equal divisions of 0.1m, denoted depthmwdd_1，depthmwdd_2，……，depthmwdd_mWherein m represents the number of sequence numbers;
step 2.2, calculate depthmwdd_1，depthmwdd_2，……，depthmwdd_mInc and azm, in particular depthmwdd_mIn (depthmwd)ti，incti，azmti)，(depthmwd1，inc1，azm1)，(depthmwd2，inc2，azm2)，……(depthmwdk，inck，azmk)，(depthgrrlast，inck，amzk) Two adjacent depth points are (depthmwd)n，incn，azmn)，(depthmwdn+1，inc n+1，azm n+1) Then depthmwdd_mWell deviation inc of pointsd_m，azmd_mIn order to realize the purpose,
wherein depthmwdd_mWell deviation inc of pointsd_m，azmd_mFormula, all points (depthmwd) can be obtainedd_1，incd_1，azmd_1)，(depthmwdd_2，incd_2，azmd_2)，……，(depthmwdd_m，incd_m，azmd_m)；
Step 3, calculating each measuring point (depthmwd)d_m，tvdd_m，nsd_m，ewd_m)；
Step 3.1, data (depthmwd) calculated in step 2d_1，incd_1，azmd_1)，(depthmwdd_2，incd_2，azmd_2)，……，(depthmwdd_m，incd_m，azmd_m) Tvd, ns, ew for each point can be calculated, and the currently measured point tvd2, ns2 and ew2 can be calculated from the (depthmwd2, inc2, azm2) measured by the instrument at the current point, as follows,
wherein:
γ＝arccos[cos(inc1)cos(inc2)+sin(iHc1)sin(inc2)cos(azm2-azm1)]
ΔL＝depthmwd2-depthmwd1
in the above formula, depthmwd1, inc1, azm1, tvd1, ns1 and ew1 are respectively expressed as the depth, inclination, azimuth, dip, north-south displacement and east-west displacement of the last measured point;
step 3.2, the following data (depthmwd) can be obtained through cyclic calculationti，tvdti，nsti，ewti)，(depthmwdd_1，tvdd_1，nsd_1，ewd_1)，(depthmwdd_2，tvdd_2，nsd_2，ewd_2)，……，(depthmwdd_m，tvdd_m，nsd_m，ewd_m)；
Step 4, rendering and calculating the updated data and images to finish real-time imaging;
step 4.1, in order to make the three-dimensional gamma imaging have a real-time function, the data volume of one-time rendering should not be too large, so the calculated (depthmwd, tvd, ns, ew) and (depthgr, GR) need to be comparedi) Matching is carried out, and only updated data are rendered each time;
step 4.2, respectively taking out the latest measured depth depthmwd of the well deviation azimuthkAnd depth depthgr of the latest measurement point of the formation gamma GRlastAnd is in (depthmwd)ti，tvdti，nsti，ewti)，(depthmwdd_1，tvdd_1，nsd_1，ewd_1)，(depthmwdd_2，tvdd_2，nsd_2，ewd_2)，……，(depthmwdd_m，tvdd_m，nsd_m，ewd_m) And (depthgr, GR)i) Find out the depth of [ depthmwdk,depthgrlast]Data within the interval, the portion of data being used to render the image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110919136.2A CN113482596B (en) | 2021-08-11 | 2021-08-11 | Real-time three-dimensional gamma imaging data processing method based on LWD while-drilling instrument |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110919136.2A CN113482596B (en) | 2021-08-11 | 2021-08-11 | Real-time three-dimensional gamma imaging data processing method based on LWD while-drilling instrument |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113482596A true CN113482596A (en) | 2021-10-08 |
| CN113482596B CN113482596B (en) | 2023-09-15 |
Family
ID=77946276
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110919136.2A Active CN113482596B (en) | 2021-08-11 | 2021-08-11 | Real-time three-dimensional gamma imaging data processing method based on LWD while-drilling instrument |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113482596B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0208767D0 (en) * | 2001-05-02 | 2002-05-29 | Schlumberger Holdings | Welllbore caliper measurement method using measurements from a gamma-gamma density instrument |
| CN101839132A (en) * | 2010-04-21 | 2010-09-22 | 中国石油化工集团 | Electromagnetic induction MWD (Measurement While Drilling) data transmission system |
| CN102140889A (en) * | 2011-03-09 | 2011-08-03 | 邱世军 | Horizontal well logging while drilling geo-steering method |
| CN203742581U (en) * | 2014-03-12 | 2014-07-30 | 北京捷威思特科技有限公司 | Pi logging while drilling system |
| CN105484732A (en) * | 2015-12-28 | 2016-04-13 | 上海神开石油设备有限公司 | Method for processing well depth in drilling geosteering construction process of horizontal well |
| CN105545284A (en) * | 2015-12-14 | 2016-05-04 | 中国石油天然气集团公司 | While-drilling gamma imaging data processing method |
| CN111308568A (en) * | 2020-03-26 | 2020-06-19 | 中国石油天然气集团有限公司 | Formation dip angle automatic pickup method based on while-drilling gamma imaging logging |
-
2021
- 2021-08-11 CN CN202110919136.2A patent/CN113482596B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0208767D0 (en) * | 2001-05-02 | 2002-05-29 | Schlumberger Holdings | Welllbore caliper measurement method using measurements from a gamma-gamma density instrument |
| CN101839132A (en) * | 2010-04-21 | 2010-09-22 | 中国石油化工集团 | Electromagnetic induction MWD (Measurement While Drilling) data transmission system |
| CN102140889A (en) * | 2011-03-09 | 2011-08-03 | 邱世军 | Horizontal well logging while drilling geo-steering method |
| CN203742581U (en) * | 2014-03-12 | 2014-07-30 | 北京捷威思特科技有限公司 | Pi logging while drilling system |
| CN105545284A (en) * | 2015-12-14 | 2016-05-04 | 中国石油天然气集团公司 | While-drilling gamma imaging data processing method |
| CN105484732A (en) * | 2015-12-28 | 2016-04-13 | 上海神开石油设备有限公司 | Method for processing well depth in drilling geosteering construction process of horizontal well |
| CN111308568A (en) * | 2020-03-26 | 2020-06-19 | 中国石油天然气集团有限公司 | Formation dip angle automatic pickup method based on while-drilling gamma imaging logging |
Non-Patent Citations (1)
| Title |
|---|
| 刘春生;陈晓林;侯岳;刘聃;何楠;: "方位伽马随钻测量技术在土耳其天然碱溶采对接井中的应用", 探矿工程(岩土钻掘工程), no. 08 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113482596B (en) | 2023-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105899749B (en) | The stratum and interpretation of structure of deviated wellbore and horizontal hole | |
| CN106597533B (en) | A kind of Depth Domain velocity modeling method for mountain front seism processing | |
| US7925483B2 (en) | Methods for visualizing distances between wellbore and formation boundaries | |
| CN107783187B (en) | Method for establishing three-dimensional velocity field by combining logging velocity and seismic velocity | |
| US8473212B2 (en) | Log processing in highly deviated wellbores | |
| CN108005646B (en) | Stratum anisotropic resistivity extraction method based on electromagnetic wave logging while drilling data | |
| CN109061764A (en) | A kind of frequency dividing fusion Optimum Impedance Inversion Method | |
| CN107861917B (en) | Method of the multivariate data combined calculation wellbore to stratigraphic boundary distance in horizontal well | |
| CN110107282B (en) | Method for controlling while-drilling track of thin-layer shale oil horizontal well | |
| US11435497B2 (en) | Three dimensional visualization from point-by-point one dimensional inversion with bed azimuth | |
| US20140222345A1 (en) | Methods and systems of incorporating pseudo-surface pick locations in seismic velocity models | |
| CN104809266A (en) | Spline based accurate predicating method for face coal seam occurrence condition | |
| US11525353B2 (en) | System and methods for evaluating a formation using pixelated solutions of formation data | |
| DE112016005364T5 (en) | A method of selecting an earth model from a plurality of earth models | |
| CN113482596A (en) | LWD-based while-drilling instrument real-time three-dimensional gamma imaging data processing method | |
| US11703610B2 (en) | Enhanced inversion imaging | |
| CN115434688B (en) | Drilling curve control method for logging while drilling of horizontal well | |
| Shao et al. | Study of real-time LWD data visual interpretation and geo-steering technology | |
| CN106504319B (en) | Reservoir Three Dimensional Contrast map generalization method and device between well | |
| CN110761780A (en) | Three-dimensional geosteering method based on well-seismic combination | |
| CN116953797B (en) | Lake-phase thin-layer sandstone reservoir prediction method based on typical high-frequency gyratory isochrone contrast | |
| CN105484732B (en) | Processing method for horizontal well drilling geosteering work progress well depth | |
| CN112901070B (en) | Method for tracking horizontal well while drilling and control method for controlling drilling direction of drill bit | |
| CN108734781A (en) | A kind of stratigraphic model construction method | |
| CN113530527A (en) | While-drilling azimuth gamma imaging data processing method based on single detector |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |