CN115880482B - Logging image card identification and inclination angle correction method and device and computing equipment - Google Patents
Logging image card identification and inclination angle correction method and device and computing equipment Download PDFInfo
- Publication number
- CN115880482B CN115880482B CN202310130957.7A CN202310130957A CN115880482B CN 115880482 B CN115880482 B CN 115880482B CN 202310130957 A CN202310130957 A CN 202310130957A CN 115880482 B CN115880482 B CN 115880482B
- Authority
- CN
- China
- Prior art keywords
- card
- imaging
- encountering
- depth
- target
- 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.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000003384 imaging method Methods 0.000 claims abstract description 185
- 238000012545 processing Methods 0.000 claims abstract description 18
- 238000004891 communication Methods 0.000 claims description 15
- 238000013519 translation Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 10
- 238000003702 image correction Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method, a device and a computing device for identifying a well logging image when encountering a card and correcting an inclination angle, according to the technical scheme provided by the invention, ultrasonic imaging well logging is carried out in a target depth interval in a well to obtain an original amplitude imaging and an original travel time imaging; processing the original travel time imaging to obtain a card-encountering indication curve of the target depth interval; determining the length of a card encountering depth section in a card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain a corrected amplitude imaging of a target depth section; and performing inclination angle identification and correction on the correction amplitude imaging to obtain a target inclination angle value. According to the invention, the indication curve of the encountering card can be obtained based on the analysis of the original logging data, so that the stretching filling of the original amplitude imaging is completed, the angle and the depth of the crack in the original amplitude imaging are corrected, the real dip angle value and the real dip depth in the geological structure are obtained, and the detection performance of the instrument is effectively improved.
Description
Technical Field
The invention relates to the field of oil well exploration, in particular to a method, a device, computing equipment and a storage medium for identifying a well logging image when encountering a card and correcting an inclination angle.
Background
In the oil and gas exploration and development process, the ultrasonic imaging logging instrument can visually give out geological structure information of stratum on the inner wall of the well through amplitude and time imaging. The cable ultrasonic imaging logging instrument performs complex non-uniform motion underground, serious distortion can occur in certain stratum images, and difficulty is brought to identification of geological structures such as cracks, layer interfaces and the like. This is mainly due to the fact that the instrument is not centered in highly deviated wells or to the damping effect of irregular wellbores and mud cakes etc. on the instrument. Furthermore, the instrument in the blocking section always collects the geological structure information of the current depth, and the geological structure information of other depth points is not collected. Therefore, only the geological structure information of the instrument when the instrument is released from the clamp can be measured on the ultrasonic imaging curve, and the geological structure information of the interval when the instrument is in contact with the clamp is missing.
In the prior art, the method for solving the problem of ultrasonic imaging image deletion is a common Kalman filtering method, and the method uses the variance and the absolute average value of an acceleration curve within a certain window length as an indication of encountering a card. However, when the acceleration value is not changed greatly, the method is difficult to distinguish whether the card is encountered or not.
Disclosure of Invention
The present invention has been made in view of the above problems, and provides a log image encounter card identification and inclination correction method, apparatus, computing device, and computer storage medium that overcome or at least partially solve the above problems.
According to one aspect of the invention, there is provided a method for identifying a log image and correcting an inclination angle, comprising:
performing ultrasonic imaging logging in a target depth interval in a well to obtain an original amplitude imaging and an original travel time imaging;
processing the original travel time imaging to obtain a card indication curve of the target depth interval;
determining the length of a card encountering depth section in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain a corrected amplitude imaging of the target depth section;
and performing inclination angle identification and correction on the correction amplitude imaging to obtain a target inclination angle value.
In the above scheme, the processing the original travel time imaging to obtain the card indication curve of the target depth interval further includes:
performing eccentric correction on the original travel time imaging to obtain a corrected travel time imaging when the logging instrument is positioned at the central position in the well;
calculating a corresponding reflected wave travel time average value of one circle in a scanning well of the logging instrument at each depth point position in the target depth interval, and obtaining a travel time average value curve of the target depth interval;
and deriving the travel time average value curve to obtain a card-contacting indication curve of the target depth interval.
In the above scheme, the performing the eccentric correction on the original travel time imaging to obtain the corrected travel time imaging when the logging instrument is located at the central position in the well further includes:
and carrying out eccentric correction on the original travel-time imaging by adopting an ellipse fitting method to obtain the corrected travel-time imaging.
In the above scheme, determining the length of the card encountering depth section in the card encountering indication curve according to the preset threshold, and stretching the original amplitude imaging according to the length of the card encountering depth section, so as to obtain the corrected amplitude imaging of the target depth section further comprises:
determining a depth point position of which the derivative value is smaller than the preset threshold value in the card encountering indication curve as a card encountering position, determining depth section lengths corresponding to a plurality of continuous card encountering positions as card encountering depth section lengths, and determining other depth section lengths except the card encountering depth section lengths in the card encountering indication curve as non-card encountering depth section lengths;
and stretching the image corresponding to the length of the non-card-encountering depth section in the original amplitude imaging to fill the image corresponding to the length of the card-encountering depth section in the original amplitude imaging, so as to obtain the corrected amplitude imaging.
In the above solution, the stretching the image corresponding to the length of the non-card-encountering depth segment in the original amplitude imaging to fill the image corresponding to the length of the card-encountering depth segment in the original amplitude imaging, and obtaining the corrected amplitude imaging further includes:
for each card encountering depth section length, determining the length of the non-card encountering depth section with depth before and adjacent to the card encountering depth section length as the target non-card encountering depth section length;
and stretching the corresponding image of the target non-card-encountering depth section length in the original amplitude imaging in the depth direction so as to fill the corresponding image of the card-encountering depth section length in the original amplitude imaging.
In the above solution, the performing inclination angle identification and correction on the corrected amplitude imaging to obtain a target inclination angle value further includes:
performing inclination angle identification on the corrected amplitude imaging according to the geological structure characteristics in the corrected amplitude imaging to obtain a first inclination angle value;
correcting the first inclination angle value according to a preset angle correction formula to obtain a second inclination angle value;
and carrying out translation on the second inclination angle value in the depth direction to obtain the target inclination angle value.
In the above scheme, after obtaining the target inclination angle value, the method further includes:
and obtaining a target stratum structure dip angle curve according to the target dip angle value.
According to another aspect of the present invention, there is provided a log image encounter card recognition and inclination correction apparatus, the apparatus comprising: the device comprises an acquisition module, an imaging processing module, an imaging correction module and an inclination correction module; wherein,,
the acquisition module is suitable for carrying out ultrasonic imaging logging in a target depth interval in the well to obtain an original amplitude imaging and an original travel time imaging;
the imaging processing module is suitable for processing the original travel time imaging to obtain a card indication curve of the target depth interval;
the imaging correction module is suitable for determining the length of a card encountering depth section in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain corrected amplitude imaging of the target depth section;
and the inclination angle correction module is suitable for carrying out inclination angle identification and correction on the correction amplitude imaging to obtain a target inclination angle value.
According to another aspect of the present invention, there is provided a computing device comprising: the device comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete communication with each other through the communication bus;
the memory is used for storing at least one executable instruction, and the executable instruction enables the processor to execute the operation corresponding to the inclination angle correction method and the log image card identification.
According to another aspect of the present invention, there is provided a computer storage medium having stored therein at least one executable instruction for causing a processor to perform operations corresponding to the method for dip correction in a log image as described above.
According to the technical scheme provided by the invention, ultrasonic imaging logging is carried out in a target depth interval in a well to obtain an original amplitude imaging and an original travel time imaging; processing the original travel time imaging to obtain a card indication curve of the target depth interval; determining the length of a card encountering depth section in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain a corrected amplitude imaging of the target depth section; and performing inclination angle identification and correction on the correction amplitude imaging to obtain a target inclination angle value. The scheme provides an effectual card recognition scheme, in having solved prior art from this, the distortion appears in ultrasonic imaging, the logging instrument can't acquire the geological structure information of this degree of depth, and then lead to the geological structure information of meeting the card interval to miss, and be difficult to distinguish whether meet the problem of card, get through original time imaging and meet the card instruction curve, and then confirm the position and the length of meeting the card, in carrying out tensile correction to original range imaging based on this position and length, and carry out inclination discernment and correction in the range imaging of correction, obtain target inclination value and target geological structure inclination curve. Therefore, the influence caused by distortion in amplitude imaging is basically eliminated, the inclination angle information in the geological structure can be accurately obtained, the real inclination angle value and the depth thereof in the geological structure can be accurately obtained through correction, and the accuracy of well logging and the detection performance of a well logging instrument are effectively improved.
The foregoing description is only an overview of the present invention, and is intended to be implemented in accordance with the teachings of the present invention in order that the same may be more clearly understood and to make the same and other objects, features and advantages of the present invention more readily apparent.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to designate like parts throughout the figures. In the drawings:
FIG. 1 is a flow chart of a method for identifying a log image and correcting an inclination angle according to an embodiment of the invention;
FIG. 2 is a flow chart of a method for determining a card indicator curve according to another embodiment of the invention;
FIG. 3 is a schematic view showing the evolution of a curve during acquisition of a card indicating curve according to another embodiment of the present invention;
FIG. 4 shows a flow diagram of a method of amplitude image correction according to one embodiment of the invention;
FIG. 5 is a flow chart of a tilt angle identification and correction method according to one embodiment of the invention;
FIG. 6 shows a schematic diagram of the evolution of a curve during correction of an amplitude image according to one embodiment of the invention;
FIG. 7 illustrates a schematic diagram of image changes during tilt angle identification and correction according to one embodiment of the present invention;
FIG. 8 shows a block diagram of a log image encounter card identification and dip correction apparatus in accordance with one embodiment of the invention;
FIG. 9 illustrates a schematic diagram of a computing device, according to an embodiment of the invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
FIG. 1 is a flow chart of a method for identifying a log image and correcting an inclination angle according to an embodiment of the invention; as shown in fig. 1, the method comprises the steps of:
step S101, ultrasonic imaging logging is carried out in a target depth interval in the well, and original amplitude imaging and original travel time imaging are obtained.
Specifically, performing ultrasonic imaging logging in the target depth interval to obtain imaging data of original amplitude and reflected wave travel time; and obtaining the original travel-time imaging according to the reflected wave travel-time imaging data.
Step S102, processing the original travel time imaging to obtain a card indication curve of the target depth interval.
Step S103, determining the length of a card encountering depth section in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain the corrected amplitude imaging of the target depth section.
Step S104, inclination angle identification and correction are carried out on the correction amplitude imaging, and a target inclination angle value is obtained.
According to the method for identifying the well logging image by the aid of the card and correcting the inclination angle, ultrasonic imaging well logging is conducted in a target depth zone in a well, and an original amplitude imaging and an original travel time imaging are obtained; processing the original travel time imaging to obtain a card indication curve of the target depth interval; determining the length of a card encountering depth section in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain a corrected amplitude imaging of the target depth section; and performing inclination angle identification and correction on the correction amplitude imaging to obtain a target inclination angle value. Obtaining a card indication curve through original travel imaging, further determining the position and the length of the card, carrying out stretching correction on the original amplitude imaging based on the position and the length, and carrying out inclination angle identification and correction from the corrected amplitude imaging to obtain a target inclination angle value and a target geological structure inclination angle curve. Therefore, the influence caused by distortion in amplitude imaging is basically eliminated, the inclination angle information in the geological structure can be accurately obtained, the real inclination angle value and the depth thereof in the geological structure can be accurately obtained through correction, and the accuracy of well logging and the detection performance of a well logging instrument are effectively improved.
FIG. 2 is a flow chart of a method for determining a card indicator curve according to another embodiment of the invention; as shown in fig. 2, the method comprises the steps of:
step S201, performing eccentricity correction on the original travel time imaging, to obtain a corrected travel time imaging when the logging instrument is located at a central position in the well.
Specifically, the performing the eccentric correction on the original travel time imaging to obtain a corrected travel time imaging when the logging instrument is located at a central position in the well further includes:
and carrying out eccentric correction on the original travel-time imaging by adopting an ellipse fitting method to obtain the corrected travel-time imaging.
Step S202, calculating a travel time average value of reflected waves corresponding to one circle in a scanning well at each depth point position in the target depth interval of the logging instrument, and obtaining a travel time average value curve of the target depth interval.
Preferably, the formula for calculating the running average of the reflected wave is:
wherein N is the azimuth point number acquired by the ultrasonic logging instrument in one-circle scanning at one depth point, and Time (i) is the reflected wave travel Time of each azimuth point of the ultrasonic logging instrument, AVET m And scanning an average value of one circle of reflected wave travel time at the m depth point for the ultrasonic logging instrument.
And calculating the travel time average value of the reflected waves of all depth points in the target depth interval, and forming a travel time average value curve of the target depth interval according to the travel time average value.
And step S203, deriving the travel time average value curve to obtain a card-meeting indication curve of the target depth interval.
Preferably, the derivative formula of the running average curve is:
DZ m =（ABS(AVET m+1 -AVET m )+ABS(AVET m - AVET m-1 )）/（SMOTH-1）
wherein DZ is m A derivative value representing the reflected wave travel time average value of an m depth point; ABS is the absolute sign; SMOTH is the number of smoothing points in the depth direction, and may be 1,3,5,7, … …; AVET m+1 Scanning an average value of one-circle reflected wave travel time at an m+1 depth point for the ultrasonic logging instrument; AVET m-1 And scanning the average value of one-circle reflected wave travel time of the ultrasonic logging instrument at the m-1 depth point.
Preferably, after determining the derivative value of the reflected wave travel time average value corresponding to one depth point at a time, determining whether the determination operation of the derivative value of the reflected wave travel time average value has been performed on all depth points in the target depth interval, and if not, performing the operation on the next depth point; if so, finally obtaining a travel time average curve of the target depth interval according to the derivative value of the travel time average value of the reflected waves corresponding to each depth point in the target depth interval, and taking the travel time average curve of the target depth interval as a card-contacting indication curve in the depth interval.
Preferably, FIG. 3 is a schematic view showing the evolution of the curve during acquisition of the indication curve of a card according to another embodiment of the present invention. As shown in fig. 3, each image has a depth on the ordinate. The first image is the original amplitude image, and the image is unfolded in the clockwise direction of north-east-south-west-north, so that amplitude values in all angles in the horizontal direction are obtained; the second path is the imaging of the original travel time and represents the reflected wave travel time of each azimuth; the third time is the corrected travel time imaging obtained after correction by an ellipse fitting method; the fourth path is a travel time average value curve generated according to the travel time average value obtained by scanning each depth point for one week; the fifth path is the card-encountering indication curve obtained by deriving according to the travel time average value curve;
the depth section corresponding to the value in the card indication curve approaching zero is the card layer section when the logging instrument measures.
Fig. 4 shows a flow diagram of a method of amplitude image correction according to one embodiment of the invention, as shown in fig. 4, wherein:
and S401, determining the depth point position of which the derivative value is smaller than the preset threshold value in the card encountering indication curve as a card encountering position, determining the depth section lengths corresponding to a plurality of continuous card encountering positions as card encountering depth section lengths, and determining the other depth section lengths except the card encountering depth section lengths in the card encountering indication curve as non-card encountering depth section lengths.
Preferably, the preset threshold may be δ=0.2; screening the card-encountering indication curve according to the preset threshold, and determining the depth point position corresponding to the derivative value as the card-encountering position if the derivative value in the card-encountering indication curve is smaller than 0.2; otherwise, the card is not in the card position.
Further, the lengths of the card encountering depth sections corresponding to the plurality of card encountering positions are continuous; and in the target depth interval, the length of the depth section without encountering a card is the length of the depth section without encountering a card except the depth section with a card.
Step S402, performing stretching processing on an image corresponding to the length of the non-card-encountering depth segment in the original amplitude imaging, so as to fill the image corresponding to the length of the card-encountering depth segment in the original amplitude imaging, thereby obtaining the corrected amplitude imaging.
Specifically, for each card encountering depth section length, determining the length of the non-card encountering depth section with depth positioned before and adjacent to the card encountering depth section length as a target non-card encountering depth section length;
and stretching the corresponding image of the target non-card-encountering depth section length in the original amplitude imaging in the depth direction so as to fill the corresponding image of the card-encountering depth section length in the original amplitude imaging.
Preferably, for one card encountering depth segment d_down, taking an uncongested depth segment d_up with depth before the card encountering depth segment d_down as a target uncongested depth segment length;
and stretching the length of the target non-card-encountering depth section, namely the image corresponding to the non-card-encountering depth section D_UP in the original amplitude imaging in the depth direction, so as to fill the imaging corresponding to the card-encountering depth section D_DOWN in the original amplitude imaging.
Preferably, after the above-mentioned stretching correction operation is performed on the length of one card-encountering depth section at a time, judging whether the stretching correction operation has been performed on all the lengths of the card-encountering depth sections in the target depth section, if not, performing the above-mentioned operation on the length of the next card-encountering depth section; if yes, generating the corrected amplitude imaging for the whole target depth interval according to the result of the stretching correction operation performed on the lengths of the depth sections of the contact card in the target depth interval.
Fig. 5 shows a flow chart of a tilt angle identification and correction method according to an embodiment of the invention, as shown in fig. 5, wherein:
step S501, performing inclination angle identification on the corrected amplitude imaging according to the geologic structure feature in the corrected amplitude imaging, and obtaining a first inclination angle value.
Specifically, in the corrected amplitude imaging after the stretching correction, a plurality of points are set in an image according to the geologic structure characteristics to perform sinusoidal curve fitting, and an inclination angle value of a geologic structure (crack) is calculated, wherein the inclination angle calculation formula is as follows:
θ=arctgH/D
wherein θ is the first inclination value, and represents an inclination of a crack in the correction amplitude imaging; h is the depth corresponding to the wave crest and the wave trough of the sinusoidal curve; d is the wellbore diameter.
Step S502, correcting the first inclination angle value according to a preset angle correction formula to obtain a second inclination angle value.
Specifically, the preset angle correction formula is:
θ c =tgθ-L/D
wherein θ c The second inclination angle value is the length of the depth section of the reflected wave amplitude meeting the clamp, and the diameter of the borehole is D.
Step S503, performing a translation on the second inclination angle value in the depth direction, to obtain the target inclination angle value.
Preferably, the second inclination angle value determined in step S502 is translated in the depth direction, where the translation distance may be L/2, and L is the length of the corresponding depth segment of the card.
Further, in step S501, step S502 and step S503, in addition to determining the first inclination angle value, the second inclination angle value and the target inclination angle value, the inclination angle value is the included angle between the point on the crack axis and the plumb line, and the inclination angle is the horizontal included angle between the point on the crack axis and the target direction line from the direction line in the clockwise direction.
Preferably, after the target dip angle value is obtained, a target formation dip angle curve is obtained according to the target dip angle value.
Preferably, the multiple images of the implementation of the method shown in fig. 4 and fig. 5 are shown in fig. 6 and fig. 7:
FIG. 6 shows a schematic diagram of the evolution of a curve in an amplitude image correction process according to one embodiment of the invention, wherein:
the first image is a well deviation and Azimuth curve, wherein the well deviation curve is generated according to well deviation angles (DEVs) corresponding to the depth points, and the Azimuth curve is generated according to Azimuth angles (AZ, azimuths) corresponding to the depth points; the well inclination angle represents an included angle between a tangent line at a certain point on the axis of the well bore and a plumb line; the azimuth angle represents a horizontal angle between a point on the axis of the borehole from the north-pointing direction line, clockwise, to the target direction line.
The second image is the original amplitude imaging, wherein D_DOWN is the length of a depth section of a card, and D_UP is the length of a depth section of a target non-card corresponding to D_DOWN; D_DOWN+ is the length of another card encountering depth segment, and D_UP+ is the length of a target non-card encountering depth segment corresponding to D_DOWN+.
The third image is an original inclination angle value of the crack, wherein the inclination angle value also comprises a tendency value of the crack; the inclination value is the horizontal included angle between the point on the crack axis and the target direction line from the north direction line, and the clockwise direction is the horizontal included angle between the point on the crack axis and the target direction line; as shown in the figure, the inclination angle value is 63 degrees, and the inclination value is 279 degrees.
And the fourth image is the corrected amplitude image obtained by performing the stretching correction operation on the original amplitude image.
A fifth image is the second inclination angle value of the corrected crack obtained through imaging calculation according to the correction amplitude; the dip angle value also comprises a tendency value of the crack, wherein the dip angle value is 72 degrees, and the tendency value is 279 degrees.
FIG. 7 is a schematic diagram showing image changes during tilt angle identification and correction according to one embodiment of the invention, wherein:
the first image is the card-encountering indication curve obtained by the card-encountering indication curve determining method; a second image is the original amplitude imaging; a third image is the corrected amplitude imaging; the fourth image is the second inclination angle value, wherein the four groups of inclination angle values and inclination values are 79 degrees/185 degrees, 74 degrees/183 degrees, 81 degrees/221 degrees and 75 degrees/214 degrees respectively; the fifth image is the target inclination angle value, wherein four groups of target inclination angle values corresponding to the fourth image are 78 degrees/185 degrees, 70 degrees/183 degrees, 80 degrees/221 degrees and 71 degrees/214 degrees respectively.
As can be seen from the figure, the inclination angle in the fifth image is adjusted in the depth direction compared with the inclination angle corresponding to the fourth image.
The method utilizes the characteristic that the ultrasonic reflection echo travel time curve is unchanged in time when a card is encountered, and evaluates the card-encountering sign in the image by the derivative thereof through the methods shown in fig. 2, 4 and 5; filling the reflection wave amplitude deformity section with the length of the non-blocking depth section in the target depth section by adopting a stretching correction method; and further performing angle and depth correction on the inclination angle of the geological structure (fracture) obtained from the corrected amplitude curve. Therefore, the length of the clamping position and the clamping depth section is accurately determined in the well logging process, meanwhile, correction of various numerical values of the crack is completed, and the capacity of accurately acquiring geological structure information of the well logging instrument is greatly improved.
FIG. 8 shows a block diagram of a log image encounter card identification and dip correction apparatus, as shown in FIG. 8, according to one embodiment of the invention, comprising: an acquisition module 801, an imaging processing module 802, an imaging correction module 803, and an inclination correction module 804. Wherein,,
the acquiring module 801 performs ultrasonic imaging logging in a target depth interval in the well to obtain an original amplitude imaging and an original travel time imaging.
The imaging processing module 802 is adapted to process the original travel time imaging to obtain a card indication curve of the target depth interval.
Specifically, the imaging processing module 802 is further adapted to: performing eccentric correction on the original travel time imaging to obtain a corrected travel time imaging when the logging instrument is positioned at the central position in the well; calculating a corresponding reflected wave travel time average value of one circle in a scanning well of the logging instrument at each depth point position in the target depth interval, and obtaining a travel time average value curve of the target depth interval; and deriving the travel time average value curve to obtain a card-contacting indication curve of the target depth interval.
Specifically, the imaging processing module 802 is further adapted to: and carrying out eccentric correction on the original travel-time imaging by adopting an ellipse fitting method to obtain the corrected travel-time imaging.
The imaging correction module 803 is adapted to determine a length of a card encountering depth segment in the card encountering indication curve according to a preset threshold, and stretch the original amplitude imaging according to the length of the card encountering depth segment, so as to obtain a corrected amplitude imaging of the target depth interval.
Specifically, the imaging correction module 803 is further adapted to: determining a depth point position of which the derivative value is smaller than the preset threshold value in the card encountering indication curve as a card encountering position, determining depth section lengths corresponding to a plurality of continuous card encountering positions as card encountering depth section lengths, and determining other depth section lengths except the card encountering depth section lengths in the card encountering indication curve as non-card encountering depth section lengths; and stretching the image corresponding to the length of the non-card-encountering depth section in the original amplitude imaging to fill the image corresponding to the length of the card-encountering depth section in the original amplitude imaging, so as to obtain the corrected amplitude imaging.
Specifically, the imaging correction module 803 is further adapted to: for each card encountering depth section length, determining the length of the non-card encountering depth section with depth before and adjacent to the card encountering depth section length as the target non-card encountering depth section length; and stretching the corresponding image of the target non-card-encountering depth section length in the original amplitude imaging in the depth direction so as to fill the corresponding image of the card-encountering depth section length in the original amplitude imaging.
The inclination angle correction module 804 is adapted to perform inclination angle recognition and correction on the corrected amplitude imaging to obtain a target inclination angle value.
Specifically, the tilt correction module 804 is further adapted to: performing inclination angle identification on the corrected amplitude imaging according to the geological structure characteristics in the corrected amplitude imaging to obtain a first inclination angle value; correcting the first inclination angle value according to a preset angle correction formula to obtain a second inclination angle value; and carrying out translation on the second inclination angle value in the depth direction to obtain the target inclination angle value.
Specifically, after obtaining the target tilt angle value, the tilt angle correction module 804 is further adapted to: and obtaining a target stratum structure dip angle curve according to the target dip angle value.
According to the device for identifying the well logging image and correcting the inclination angle when the well logging image meets the card, ultrasonic imaging well logging is carried out in a target depth zone in a well, and an original amplitude imaging and an original travel time imaging are obtained; processing the original travel time imaging to obtain a card indication curve of the target depth interval; determining the length of a card encountering depth section in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain a corrected amplitude imaging of the target depth section; and performing inclination angle identification and correction on the correction amplitude imaging to obtain a target inclination angle value. Obtaining a card indication curve through original travel imaging, further determining the position and the length of the card, carrying out stretching correction on the original amplitude imaging based on the position and the length, and carrying out inclination angle identification and correction from the corrected amplitude imaging to obtain a target inclination angle value and a target geological structure inclination angle curve. Therefore, the influence caused by distortion in amplitude imaging is basically eliminated, the inclination angle information in the geological structure can be accurately obtained, the real inclination angle value and the depth thereof in the geological structure can be accurately obtained through correction, and the accuracy of well logging and the detection performance of a well logging instrument are effectively improved.
The invention also provides a nonvolatile computer storage medium, and the computer storage medium stores at least one executable instruction which can execute the method for identifying the well logging image and correcting the inclination angle in the embodiment of any method.
FIG. 9 illustrates a schematic diagram of a computing device, according to an embodiment of the invention, the particular embodiment of the invention not being limited to a particular implementation of the computing device.
As shown in fig. 9, the computing device may include: a processor 902, a communication interface 904, a memory 906, and a communication bus 908.
Wherein:
processor 902, communication interface 904, and memory 906 communicate with each other via a communication bus 908.
A communication interface 904 for communicating with network elements of other devices, such as clients or other servers.
The processor 902 is configured to execute the program 910, and specifically may execute the steps related to the above-mentioned embodiment of the method for identifying a log image and correcting an inclination angle.
In particular, the program 910 may include program code including computer-operating instructions.
The processor 902 may be a central processing unit, CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement embodiments of the present invention. The one or more processors included by the computing device may be the same type of processor, such as one or more CPUs; but may also be different types of processors such as one or more CPUs and one or more ASICs.
A memory 906 for storing a program 910. Memory 906 may comprise high-speed RAM memory or may also include non-volatile memory (non-volatile memory), such as at least one disk memory.
The algorithms and displays presented herein are not inherently related to any particular computer, virtual system, or other apparatus. Various general-purpose systems may also be used with the teachings herein. The required structure for a construction of such a system is apparent from the description above. In addition, the present invention is not directed to any particular programming language. It will be appreciated that the teachings of the present invention described herein may be implemented in a variety of programming languages, and the above description of specific languages is provided for disclosure of enablement and best mode of the present invention.
In the description provided herein, numerous specific details are set forth. However, it is understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Similarly, it should be appreciated that in the foregoing description of exemplary embodiments of the invention, various features of the invention are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects. However, the disclosed method should not be construed as reflecting the intention that: i.e., the claimed invention requires more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive aspects lie in less than all features of a single foregoing disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate embodiment of this invention.
Those skilled in the art will appreciate that the modules in the apparatus of the embodiments may be adaptively changed and disposed in one or more apparatuses different from the embodiments. The modules or units or components of the embodiments may be combined into one module or unit or component and, furthermore, they may be divided into a plurality of sub-modules or sub-units or sub-components. Any combination of all features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or units of any method or apparatus so disclosed, may be used in combination, except insofar as at least some of such features and/or processes or units are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings), may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise.
Furthermore, those skilled in the art will appreciate that while some embodiments described herein include some features but not others included in other embodiments, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments. For example, in the claims, any of the claimed embodiments may be used in any combination.
Various component embodiments of the invention may be implemented in hardware, or in software modules running on one or more processors, or in a combination thereof. Those skilled in the art will appreciate that some or all of the functions of some or all of the components in accordance with embodiments of the present invention may be implemented in practice using a microprocessor or Digital Signal Processor (DSP). The present invention can also be implemented as an apparatus or device program (e.g., a computer program and a computer program product) for performing a portion or all of the methods described herein. Such a program embodying the present invention may be stored on a computer readable medium, or may have the form of one or more signals. Such signals may be downloaded from an internet website, provided on a carrier signal, or provided in any other form.
It should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and that those skilled in the art will be able to design alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The invention may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the unit claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The use of the words first, second, third, etc. do not denote any order. These words may be interpreted as names.
Claims (9)
1. A method for identifying a well logging image and correcting an inclination angle is characterized by comprising the following steps:
performing ultrasonic imaging logging in a target depth interval in a well to obtain an original amplitude imaging and an original travel time imaging;
performing eccentric correction on the original travel time imaging to obtain a corrected travel time imaging when the logging instrument is positioned at the central position in the well; calculating a corresponding reflected wave travel time average value of one circle in a scanning well of the logging instrument at each depth point position in the target depth interval, and obtaining a travel time average value curve of the target depth interval; deriving the travel time average value curve to obtain a card-contacting indication curve of the target depth interval;
determining the length of a card encountering depth section in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain a corrected amplitude imaging of the target depth section;
and performing inclination angle identification and correction on the correction amplitude imaging to obtain a target inclination angle value.
2. The method of claim 1, wherein the decentering correction of the raw travel-time imaging to obtain a corrected travel-time imaging with the logging instrument in a centered position in the well further comprises:
and carrying out eccentric correction on the original travel-time imaging by adopting an ellipse fitting method to obtain the corrected travel-time imaging.
3. The method of claim 1, wherein determining a length of a card encountering depth segment in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth segment, and obtaining the corrected amplitude imaging of the target depth interval further comprises:
determining a depth point position of which the derivative value is smaller than the preset threshold value in the card encountering indication curve as a card encountering position, determining depth section lengths corresponding to a plurality of continuous card encountering positions as card encountering depth section lengths, and determining other depth section lengths except the card encountering depth section lengths in the card encountering indication curve as non-card encountering depth section lengths;
and stretching the image corresponding to the length of the non-card-encountering depth section in the original amplitude imaging to fill the image corresponding to the length of the card-encountering depth section in the original amplitude imaging, so as to obtain the corrected amplitude imaging.
4. A method according to claim 3, wherein stretching the image corresponding to the length of the non-carded depth segment in the original amplitude imaging to fill the image corresponding to the length of the carded depth segment in the original amplitude imaging, the obtaining the corrected amplitude imaging further comprises:
for each card encountering depth section length, determining the length of the non-card encountering depth section with depth before and adjacent to the card encountering depth section length as the target non-card encountering depth section length;
and stretching the corresponding image of the target non-card-encountering depth section length in the original amplitude imaging in the depth direction so as to fill the corresponding image of the card-encountering depth section length in the original amplitude imaging.
5. The method of claim 1, wherein said tilt identifying and correcting said corrected magnitude imaging to obtain a target tilt value further comprises:
performing inclination angle identification on the corrected amplitude imaging according to the geological structure characteristics in the corrected amplitude imaging to obtain a first inclination angle value;
correcting the first inclination angle value according to a preset angle correction formula to obtain a second inclination angle value;
and carrying out translation on the second inclination angle value in the depth direction to obtain the target inclination angle value.
6. The method according to any one of claims 1-5, wherein after obtaining the target tilt angle value, the method further comprises:
and obtaining a target stratum structure dip angle curve according to the target dip angle value.
7. A log image encounter card identification and dip correction device, the device comprising: the device comprises an acquisition module, an imaging processing module, an imaging correction module and an inclination correction module; wherein,,
the acquisition module is suitable for carrying out ultrasonic imaging logging in a target depth interval in the well to obtain an original amplitude imaging and an original travel time imaging;
the imaging processing module is suitable for carrying out eccentric correction on the original travel time imaging to obtain corrected travel time imaging when the logging instrument is positioned at the central position in the well; calculating a corresponding reflected wave travel time average value of one circle in a scanning well of the logging instrument at each depth point position in the target depth interval, and obtaining a travel time average value curve of the target depth interval; deriving the travel time average value curve to obtain a card-contacting indication curve of the target depth interval;
the imaging correction module is suitable for determining the length of a card encountering depth section in the card encountering indication curve according to a preset threshold value, and stretching the original amplitude imaging according to the length of the card encountering depth section to obtain corrected amplitude imaging of the target depth section;
and the inclination angle correction module is suitable for carrying out inclination angle identification and correction on the correction amplitude imaging to obtain a target inclination angle value.
8. A computing device, comprising: the device comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete communication with each other through the communication bus;
the memory is configured to store at least one executable instruction that causes the processor to perform the operations of the log image encounter card identification method of any one of claims 1-6, corresponding to the dip correction method.
9. A computer storage medium having stored therein at least one executable instruction for causing a processor to perform operations corresponding to the borehole image encounter card identification and dip angle correction method of any one of claims 1-6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310130957.7A CN115880482B (en) | 2023-02-17 | 2023-02-17 | Logging image card identification and inclination angle correction method and device and computing equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310130957.7A CN115880482B (en) | 2023-02-17 | 2023-02-17 | Logging image card identification and inclination angle correction method and device and computing equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115880482A CN115880482A (en) | 2023-03-31 |
| CN115880482B true CN115880482B (en) | 2023-05-09 |
Family
ID=85761241
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310130957.7A Active CN115880482B (en) | 2023-02-17 | 2023-02-17 | Logging image card identification and inclination angle correction method and device and computing equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115880482B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103913776A (en) * | 2014-04-23 | 2014-07-09 | 深圳市市政工程总公司 | Method for detecting spherical boulders in shield construction |
| CN104278990A (en) * | 2013-07-02 | 2015-01-14 | 中国石油天然气集团公司 | Logging data quality recovery method and device |
| CN114233274A (en) * | 2022-02-25 | 2022-03-25 | 中海油田服务股份有限公司 | Image generation method and device based on while-drilling electrical imaging |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6704655B2 (en) * | 2000-10-12 | 2004-03-09 | Schlumberger Technology Corporation | Method and apparatus for correcting the depth index for well-log data |
| CN103806899B (en) * | 2012-11-14 | 2018-04-27 | 中国石油集团长城钻探工程有限公司 | A kind of method and apparatus for determining depth of the instrument in well |
| WO2016108906A1 (en) * | 2014-12-31 | 2016-07-07 | Halliburton Energy Services, Inc. | Pulse reflection travel time analysis to track position of a downhole object |
| CN111677494B (en) * | 2020-02-06 | 2022-03-11 | 西南石油大学 | Method, system, terminal and storage medium for judging and calculating differential pressure of logging instrument in three-dimensional curved well |
| CN114016998A (en) * | 2020-07-16 | 2022-02-08 | 中国石油天然气集团有限公司 | Prediction method and device for logging encountering block |
-
2023
- 2023-02-17 CN CN202310130957.7A patent/CN115880482B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104278990A (en) * | 2013-07-02 | 2015-01-14 | 中国石油天然气集团公司 | Logging data quality recovery method and device |
| CN103913776A (en) * | 2014-04-23 | 2014-07-09 | 深圳市市政工程总公司 | Method for detecting spherical boulders in shield construction |
| CN114233274A (en) * | 2022-02-25 | 2022-03-25 | 中海油田服务股份有限公司 | Image generation method and device based on while-drilling electrical imaging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115880482A (en) | 2023-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11808143B2 (en) | Methods and apparatus to measure formation features | |
| US10121261B2 (en) | Automatic dip picking in borehole images | |
| US10970814B2 (en) | Subsurface formation imaging | |
| CN110424945B (en) | Shale oil horizontal well deployment method and device | |
| CN115144906B (en) | Amplitude eccentricity correction method and device for ultrasonic logging | |
| CN105116448A (en) | Conversion wave azimuthal anisotropy correction method and device thereof | |
| WO2019200242A1 (en) | Evaluating casing cement using automated detection of clinging compression wave (p) arrivals | |
| CN111551927A (en) | Underground pipeline diameter measuring method based on three-dimensional ground penetrating radar | |
| CN115992692B (en) | Cement ring thickness measuring method and device, electronic equipment and storage medium | |
| CN113640808A (en) | Shallow water submarine cable buried depth detection method and device | |
| CN115880482B (en) | Logging image card identification and inclination angle correction method and device and computing equipment | |
| CN116556925A (en) | Method, device, computing equipment and storage medium for computing borehole ovality | |
| CN107780926A (en) | Borehole section shape measuring method, device and terminal | |
| CN114961703B (en) | Method and device for positioning cement plug well, electronic equipment and storage medium | |
| CN114233274B (en) | Image generation method and device based on while-drilling electrical imaging | |
| CN114876454B (en) | Method and system for identifying complex lithology flooding layer of highly-deviated well | |
| CN114439464A (en) | Horizontal well output profile flow calculation method and device based on multi-trip combination | |
| CN1696742A (en) | Method for processing data of well logging to calculate stratum attitude in drilling well | |
| CN112904348A (en) | Three-dimensional detection method, device, equipment and storage medium | |
| CN111691876A (en) | Method and device for imaging adjacent well by using acoustic logging and storage medium | |
| CN116201526B (en) | Micro annular gap detection method, micro annular gap detection device, computing equipment and storage medium | |
| CN112394415B (en) | Method and device for determining horizontal azimuth angle of submarine node | |
| CN113187457B (en) | Method and device for determining fracturing effect of reservoir and storage medium | |
| NL2028501B1 (en) | PROCEDURE FOR DETECTION OF LOW LEVEL BASED ON FULL HOLE VIEW OF ELECTRICAL IMAGE RECORDING | |
| CN115977617B (en) | Method and device for identifying string groove outside casing, computing equipment and storage medium |
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 |