A method and apparatus for correcting well axis centering for ultrasonic three-dimensional imaging

Abstract

The application provides a well shaft centering degree correction method and device for ultrasonic three-dimensional imaging, and the method comprises the following steps: capturing information of a plurality of target data points on a target casing by an ultrasonic detector along a circumferential scan; calculating a plane rectangular coordinate of each target data point in the plurality of target data points through the information of the target data points; obtaining an arc length value of a position of each target data point in the plurality of target data points according to the plane rectangular coordinate of each target data point in the plurality of target data points; determining an offset of well shaft centering of the ultrasonic detector according to the arc length values of the positions of the plurality of target data points; and correcting the plane rectangular coordinate of each target data point in the plurality of target data points according to the offset. The above scheme solves the problem of low three-dimensional imaging accuracy caused by the fact that the well shaft centering degree cannot be checked, and effectively improves the three-dimensional imaging accuracy.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of geological exploration, and particularly relates to a method and device for correcting the well shaft centering degree of ultrasonic three-dimensional imaging. BACKGROUND

[0002] Casing completion is the most important completion method in well construction. The working state of the casing is affected by factors such as cementing quality, formation dislocation, wellbore service life, hydraulic fracturing stimulation operation, injection-production measures, alternating injection-production of gas storage, and fluid properties, thereby causing deformation of the casing and posing a serious threat to the safety of oil and gas production. From the current production practice, the prevention, detection and treatment of casing deformation are very important for safety production.

[0003] Among them, casing deformation detection is the premise of casing deformation evaluation, repair and other production modification measures. At present, casing deformation detection technologies mainly include multi-arm caliper detection, electromagnetic flaw detection, and ultrasonic three-dimensional imaging detection. Among these detection technologies, the ultrasonic three-dimensional imaging detection technology has the characteristics of good directivity, strong penetration ability, and high measurement accuracy, and can realize accurate measurement and description of the deformation form of the casing, effectively solving the problems of low wellbore coverage, and difficult to guarantee measurement accuracy and resolution. The principle of casing deformation ultrasonic three-dimensional imaging detection technology is that the ultrasonic probe is installed at the lowermost end of the detection instrument. During the downhole detection process, the instrument ascends or descends at a predetermined speed, the ultrasonic probe rotates at a high speed while ascending or descending with the detection instrument, and emits ultrasonic pulse signals towards the wellbore direction. The ultrasonic pulse signals propagate to the inner wall of the casing through the wellbore working fluid, are reflected back, and the final echo signal is received by the ultrasonic probe. In this way, the data of a sampling point is captured through self-emission and self-reception. Based on the running speed of the detection instrument, the rotation speed of the ultrasonic probe, and the time required for the ultrasonic signal propagation, the time interval of the ultrasonic pulse signal excitation is optimized, and a large amount of data of sampling points can be captured by continuously exciting pulse signals according to the rules. According to the time difference of the echo signal of each sampling point, the distance between the instrument and the inner wall of the casing can be calculated, and the position coordinates of each measurement point can be obtained by combining the rotation angle of the ultrasonic probe and the azimuth correction of the gyroscope. The three-dimensional characteristics of the casing deformation form can be depicted by presenting these measurement point coordinates in three-dimensional space.

[0004] However, in actual detection conditions, especially in long horizontal sections, the centralizer cannot guarantee the stable centering state of the detection instrument, and the instrument will deviate irregularly to the well circumference direction, which changes the well axis reference at different well depths, and the relative position relationship cannot be accurately judged. In view of how to realize the casing deformation ultrasonic three-dimensional imaging well axis centering degree checking, in order to improve the accuracy of three-dimensional imaging, and then better present the real shape of casing deformation, at present, no effective solution has been put forward. SUMMARY

[0005] The present application aims to provide a kind of well axis centering degree correction method and device of ultrasonic three-dimensional imaging, can be corrected to the deviation generated by the deviation of well axis direction, so as to improve the accuracy of three-dimensional imaging.

[0006] The present application provides a kind of well axis centering degree correction method and device of ultrasonic three-dimensional imaging is realized as follows:

[0007] A kind of well axis centering degree correction method of ultrasonic three-dimensional imaging, the method comprises:

[0008] acquire the information of multiple target data points on the target casing by ultrasonic detector along the circumferential scanning;

[0009] by the information of the target data point, the plane rectangular coordinates of each target data point in the multiple target data points are calculated;

[0010] according to the plane rectangular coordinates of each target data point in the multiple target data points, the arc length value of the position of each target data point in the multiple target data points is obtained;

[0011] according to the arc length value of the position of the multiple target data points, the offset of the well axis centering of the ultrasonic detector is determined;

[0012] according to the offset, the plane rectangular coordinates of each target data point in the multiple target data points are corrected.

[0013] In one embodiment, the information of the multiple target data points includes: the azimuth angle and distance of the multiple target data points;

[0014] correspondingly, by the information of the target data point, the plane rectangular coordinates of each target data point in the multiple target data points are calculated, including:

[0015] the polar coordinates represented by the azimuth angle and distance of the multiple target data points are converted into horizontal coordinates and vertical coordinates represented plane rectangular coordinates according to the following formula;

[0016] x i =r ix cos(θ i )

[0017] y i = r i x sin(θ i )

[0018] wherein, x i represents the horizontal coordinate of the i i represents the vertical coordinate of the i i represents the distance of the i i represents the azimuth angle of the i

[0019] In one embodiment, the arc length value of the position of each of the plurality of target data points is calculated according to the plane rectangular coordinates of each of the plurality of target data points, comprising:

[0020] The arc length value of the position of each of the plurality of target data points is calculated according to the following formula:

[0021]

[0022] wherein, i takes values [1, n-1], s i represents the arc length value of the position of the i i represents the horizontal coordinate of the i i represents the vertical coordinate of the i i+1 represents the horizontal coordinate of the i+1 i+1 represents the vertical coordinate of the i+1

[0023]

[0024] wherein, s n represents the arc length value of the position of the n n represents the horizontal coordinate of the n n represents the vertical coordinate of the n

[0025] In one embodiment, the offset of the well shaft centering of the ultrasonic detector is determined according to the arc length values of the positions of the plurality of target data points, comprising:

[0026] The minimum arc length value is selected from the arc length values of the positions of the plurality of target data points as a reference value;

[0027] a ratio between the arc length value of each target data point and the reference value is calculated, and the calculation result is used as a weighting factor of each target data point;

[0028] a geometric center coordinate of the circumferential curve is calculated according to the weighting factor of each target data point;

[0029] a coordinate of a position where the ultrasonic detector is located is obtained as a reference origin coordinate;

[0030] the geometric center coordinate is compared with the reference origin coordinate to determine an offset of the ultrasonic detector from the well axis.

[0031] In one embodiment, the geometric center coordinate of the circumferential curve is calculated according to the weighting factor of each target data point, including:

[0032] the geometric center coordinate of the circumferential curve is calculated according to the following formula:

[0033]

[0034]

[0035] wherein X represents an abscissa of the geometric center coordinate, Y represents an ordinate of the geometric center coordinate, w i represents the weighting factor of the i-th target data point, w sum represents an accumulated value of the weighting factors of the target data points, x i represents the abscissa of the i-th target data point, y i represents the ordinate of the i-th target data point.

[0036] In one embodiment, the ultrasonic detector is applied in detection of casing deformation.

[0037] A device for correcting a well axis centering degree of ultrasonic three-dimensional imaging, including:

[0038] a scanning module for capturing information of a plurality of target data points on a target casing by an ultrasonic detector along a circumferential scanning;

[0039] a calculation module for calculating a plane rectangular coordinate of each target data point in the plurality of target data points through the information of the target data points;

[0040] a calculation module for calculating a plane rectangular coordinate of each target data point in the plurality of target data points through the information of the target data points;

[0041] determining a deviation of the shaft center of the ultrasonic detector according to the arc length values of the positions of the target data points;

[0042] correcting the plane rectangular coordinates of each of the target data points according to the deviation.

[0043] In an embodiment, the determining module comprises:

[0044] selecting a minimum arc length value from the arc length values of the positions of the target data points as a reference value;

[0045] calculating a ratio between the arc length value of each of the target data points and the reference value, and taking the calculation result as a weighting factor of each of the target data points;

[0046] calculating the geometric center coordinates of the circumferential curve according to the weighting factor of each of the target data points;

[0047] acquiring the position coordinates of the ultrasonic detector as reference origin coordinates;

[0048] comparing the geometric center coordinates with the reference origin coordinates to determine the deviation of the shaft center of the ultrasonic detector.

[0049] An electronic device comprising a processor and a memory for storing processor-executable instructions, the processor implementing the steps of the above method when executing the instructions.

[0050] A computer-readable storage medium having stored thereon a computer program / instructions, the computer program / instructions being executed by a processor to implement the steps of the above method.

[0051] The well axis centering degree correction method for ultrasonic three-dimensional imaging provided in the application comprises the following steps: an ultrasonic detector captures information of a plurality of target data points on a target casing through circumferential scanning; the planar rectangular coordinates of each target data point in the plurality of target data points are calculated based on the information of the target data points; the arc length values of the positions of each target data point in the plurality of target data points are obtained based on the planar rectangular coordinates of each target data point in the plurality of target data points; the offset of the well axis centering of the ultrasonic detector is determined based on the arc length values of the positions of the plurality of target data points; and the planar rectangular coordinates of each target data point in the plurality of target data points are corrected based on the offset. That is, the eccentricity of the detector is determined by using the arc length difference between the data points, and the eccentricity obtained by calculation can be used to check the well axis centering degree of the casing deformation ultrasonic three-dimensional imaging, so as to solve the problem of low three-dimensional imaging accuracy caused by the fact that the well axis centering degree cannot be checked, and achieve the technical effect of effectively improving the three-dimensional imaging accuracy. BRIEF DESCRIPTION OF DRAWINGS

[0052] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed in the embodiments or the prior art description will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments described in the present application, and other drawings can also be obtained by those skilled in the art without creative labor.

[0053] Figure 1 is a geometric relationship schematic diagram of the casing deformation ultrasonic three-dimensional imaging detector in an eccentric position in the working process of the casing deformation ultrasonic three-dimensional imaging detector provided in the present application;

[0054] Figure 2 is a standard circular casing horizontal and vertical coordinate offset distance schematic diagram provided in the present application;

[0055] Figure 3 is a standard circular casing horizontal and vertical coordinate offset distance correction calculation value change relationship schematic diagram with the number of data points provided in the present application;

[0056] Figure 4 is a method flowchart of one embodiment of the well axis centering degree correction method for ultrasonic three-dimensional imaging provided in the present application;

[0057] Figure 5 is a position relationship schematic diagram of a standard circular casing when scanned and detected according to 32 data points provided in the present application;

[0058] Figure 6 is a position and distance distribution curve of different data points when a standard circular casing is scanned and detected according to 32 data points provided in the present application;

[0059] Figure 7 is a schematic diagram of the relative position relationship of the standard circular sleeve before and after the centering correction provided in the present application;

[0060] Figure 8 is a schematic diagram of the position relationship when the elliptical sleeve is scanned and detected according to 32 data points provided in the present application;

[0061] Figure 9 is a curve of the azimuth and distance distribution of different data points when the elliptical sleeve is scanned and detected according to 32 data points provided in the present application;

[0062] Figure 10 is a schematic diagram of the relative position relationship of the elliptical sleeve before and after the centering correction provided in the present application;

[0063] Figure 11 is a schematic diagram of the position relationship when the irregularly shaped sleeve is scanned and detected according to 32 data points provided in the present application;

[0064] Figure 12 is a curve of the azimuth and distance distribution of different data points when the irregularly shaped sleeve is scanned and detected according to 32 data points provided in the present application;

[0065] Figure 13 is a schematic diagram of the relative position relationship of the irregularly shaped sleeve before and after the centering correction provided in the present application;

[0066] Figure 14 is a hardware structure block diagram of an electronic device of a well axis centering correction method of ultrasonic three-dimensional imaging provided in the present application;

[0067] Figure 15 is a schematic diagram of the module structure of an embodiment of the well axis centering correction device of ultrasonic three-dimensional imaging provided in the present application. DETAILED DESCRIPTION

[0068] In order to enable persons skilled in the art to better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be described clearly and completely in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by persons skilled in the art without creative labor should be within the scope of protection of the present application.

[0069] In order to solve the technical problems that the casing deformation ultrasonic three-dimensional imaging detector is not centered, the relative position relationship of the casing imaging at different well depths cannot be accurately judged, and the real shape distortion of the casing deformation is distorted, in this case, a casing deformation ultrasonic three-dimensional imaging well axis centering degree checking method is provided. Each data point captured by the ultrasonic three-dimensional imaging detector contains the azimuth of the casing at the point and the distance from the ultrasonic probe. The data points around the circle are simplified as a plane, that is, the coordinates of the point in the plane coordinate system can be calculated according to the angle and distance of the data point. Because the ultrasonic probe rotates at a constant speed, the difference between the azimuth angles of two adjacent data points is equal. As shown in Figure 1 If the casing is a standard circle and the detector is at the center position at the same time, the distribution of the data points obtained by scanning on the circumference is uniform, that is, the arc length of the adjacent data points is equal. If the casing is a standard circle, but the detector is at the eccentric position, the distribution of the data points obtained by scanning on the circumference is not uniform, the data points near the detector position are more dense, and the data points far from the detector position are more sparse. The density and sparsity of the data point distribution is directly reflected in the arc length between adjacent data points, that is, when the data point distribution is dense, the arc length between adjacent data points is short, and when the data point distribution is sparse, the arc length between adjacent data points is long.

[0070] Based on this, the arc length difference between the data points can be used to assign a weighted value to each data point, and the weighted average of the horizontal and vertical coordinates of all data points is calculated according to the weighted value of each data point. In the case that the number of detection data is sufficient, the weighted average can approximately represent the geometric center coordinates of the circumference curve within the error range, and the difference between the calculated circumference curve geometric center coordinates and the origin coordinates is the eccentricity of the detector. The eccentricity obtained by calculation can check the casing deformation ultrasonic three-dimensional imaging well axis centering degree.

[0071] That is, the plane coordinates of each data point can be calculated according to the data point information captured by the instrument, for example, one week calculation, 360° data point coordinates are obtained, and the horizontal and vertical coordinates of the i th data point are respectively denoted as x i and y i The i th data point information captured by the instrument includes the azimuth angle and the distance, which are respectively denoted as θ i and r i There are n data points in 360°, and the polar coordinates are converted into rectangular coordinates, that is:

[0072] x i =r i ×cos(θ i ) (Formula 1)

[0073] yi = r i x sin(θ i ) (Formula 2)

[0074] Where, i = 1, 2…n, the distribution of data points captured by the detector in one scan is uneven on the circumference curve, and the geometric center (representing the central point position) of the circumference curve cannot be directly calculated by the average of the horizontal and vertical coordinates of each data point. The density and sparsity of data points are directly reflected in the arc length between adjacent data points, that is, when the data points are densely distributed, the arc length between adjacent data points is short, and when the data points are sparsely distributed, the arc length between adjacent data points is long. Therefore, the arc length difference between data points can be used to assign a weighted value to each data point, and the weighted average of the horizontal and vertical coordinates of each data point is used to solve the weighted average of all data points. In the case of sufficient detection data, the weighted average can approximately represent the geometric center coordinates of the circumference curve within the error range, and by comparing the calculated geometric center coordinates with the origin coordinates of the polar coordinate to Cartesian coordinate conversion reference, the offset distance of the current three-dimensional imaging curve detector can be obtained.

[0075] First, the arc length of each data point is calculated, and the distance between a data point with a certain serial number and the next data point is used to approximate the arc length of the data point with the serial number, so that:

[0076]

[0077] Where, i = 1, 2…n-1.

[0078] The arc length of the nth data point can be calculated according to the following formula:

[0079]

[0080] The smallest value is selected from all arc lengths, and the weighted factor is calculated according to the ratio of the arc length of each data point to the minimum value as the reference. The weighted average of the horizontal and vertical coordinates of the data points is calculated, and the minimum value in the arc length can be represented as:

[0081] s min = min(s) (Formula 5)

[0082] The weighted factor of each data point is:

[0083]

[0084] The sum of the weighted factors of all data points is:

[0085]

[0086] The geometric center coordinates of the circumferential curve are:

[0087]

[0088]

[0089] Because, in processing the original data point information, the position of the detector is the origin when the polar coordinates are directly converted into rectangular coordinates, therefore, the difference between the weighted average coordinates calculated by formula 8 and formula 9 and the origin is the offset distance, so the offset distance is:

[0090] X 偏移距 = X - 0 (formula 10)

[0091] Y 偏移距 = Y - 0 (formula 11)

[0092] The data points are checked from the eccentric position to the centered position, and are converted as follows:

[0093] x i校正值 = x i - X 偏移距 (formula 12)

[0094] y i校正值 = y i - Y 偏移距 (formula 13)

[0095] wherein formula 12 and formula 13 are the corrected data point coordinates.

[0096] Numerical verification is performed on the above method, as shown in Figure 2 is a standard circular sleeve of a plane, the centered point is as shown in Figure 2 (1), located at the center, the eccentric point is as shown in Figure 2 (2), located in the third quadrant of the coordinate system, when the instrument is located at the eccentric point, the path of the emitted signal is as shown in Figure 2 (3), these paths form a straight line system, which is composed of straight lines passing through the eccentric point and having different slopes, the intersection coordinates of the straight line and the circle can be calculated from the equations of the straight line and the circle, the center coordinates of the circle are (0, 0), the radius is r, the angle between the straight line and the x-axis is α, that is, the slope k = tan(α), the eccentric point coordinates passed by the straight line are (-dx, -dy), dx is the offset distance of the eccentric point in the x-axis direction, dy is the offset distance of the eccentric point in the y-axis direction, the intercept b of the straight line on the y-axis is b = k × dx - dy, and then the intersection coordinates of the straight line and the circle are calculated.

[0097] First, several calculation parameters are defined:

[0098] d a= 1 + k 2 ( Equation 14 )

[0099] d b = 2k x b ( Equation 15 )

[0100] d c = b 2 - r 2 ( Equation 16 )

[0101] The straight line and the circle have two intersection points, and the intersection point coordinates are respectively:

[0102]

[0103] Y1 = -k x X1 + b ( Equation 18 )

[0104]

[0105] Y2 = k x X2 + b ( Equation 20 )

[0106] Suppose there are n data points (n is an even number), that is, n straight lines intersect the circle, according to the symmetry principle, only half of the circle, that is, the rotation of the straight line within the range of 180 degrees can calculate all the data points on the circle, so all the data points on the circle are calculated according to the horizontal and vertical coordinates of the formula 17 to formula 20. According to the horizontal and vertical coordinates of each data point, the distance from the data point to the eccentric point can be solved, so the distance from each data point to the detector and the rotation position can be solved, and then substituted into the above formula 1 to formula 13, the offset distance can be calculated.

[0107] The above formula can generate simulation detection data points of a standard circular sleeve, and verify the established calibration model. As shown in Figure 1 , it is assumed that the standard circle radius is 1 and the eccentric point coordinates are (-0.5, -0.5), that is, the horizontal coordinate offset distance of the detector is 0.5 and the vertical coordinate offset distance is 0.5. Then according to the above model, select different number of data points to calculate the offset distance of the detector, and compare with the actual value of the offset distance. As shown in Figure 3 , the relationship between the calibration calculation value of the horizontal and vertical coordinate offset distance and the number of data points is shown. With the increase of the number of data points, the calibration calculation value of the horizontal and vertical coordinate offset distance gradually approaches the actual value. When the number of data points reaches 202, the calibration calculation error of the horizontal coordinate offset distance is 0.64%, and the calibration calculation error of the vertical coordinate offset distance is 0.66%, which reaches a high calculation precision.

[0108] The above method can be applied in the ultrasonic detection of the deformation of the casing of a shale gas well. Based on the feature that the data points captured by the ultrasonic detector located at the eccentric position are unevenly distributed in the circumferential scanning, different weight proportions are assigned according to the arc length difference between the data points, the weighted average of the horizontal and vertical coordinates of all the data points is calculated, and then the geometric center coordinates of the circumferential curve and the offset distance of the detector are calculated according to the weighted average. Finally, the data point coordinates on the casing around the well are transformed in combination with the geometric change relationship, so that the casing deformation ultrasonic three-dimensional imaging well axis centralization degree data checking can be realized. The above scheme checks the centralization degree of the well axis, so that the offset distance of the detector can be quickly and accurately calculated, the relative position relationship of the casing imaging at different well depths can be accurately judged, the casing deformation ultrasonic three-dimensional imaging well axis centralization degree checking is realized, and the three-dimensional imaging accuracy is improved, so that the real shape of the casing deformation is better presented.

[0109] Figure 4 is a method flowchart of an embodiment of the well axis centralization degree correction method of the ultrasonic three-dimensional imaging provided by the present application. Although the present application provides the method operation steps or device structures as described in the following embodiments or drawings, more or fewer operation steps or module units can be included in the method or device based on conventional or non-creative labor. In steps or structures that do not have necessary causal relationships in logic, the execution order of these steps or the module structure of the device is not limited to the execution order or module structure shown in the embodiments described and the drawings of the present application. When the method or module structure is applied in actual device or terminal product, it can be sequentially executed or executed in parallel (for example, parallel processor or multi-thread processing environment, even distributed processing environment) according to the method or module structure shown in the embodiments or drawings.

[0110] Specifically, as shown in Figure 4 the above-mentioned well axis centralization degree correction method of the ultrasonic three-dimensional imaging can include the following steps:

[0111] Step 401: capturing information of a plurality of target data points on a target casing by an ultrasonic detector in circumferential scanning;

[0112] Step 402: calculating the plane rectangular coordinates of each target data point in the plurality of target data points according to the information of the target data points;

[0113] The information of the plurality of target data points can include the azimuth angle and distance of the plurality of target data points.

[0114] Correspondingly, the polar coordinates represented by the azimuth angle and distance of the plurality of target data points can be converted into the plane rectangular coordinates represented by the horizontal and vertical coordinates according to the following formula:

[0115] x i = r i cos(θ i )

[0116] y i = r i sin(θ i )

[0117] Wherein, x i represents the horizontal coordinate of the i-th target data point, y i represents the vertical coordinate of the i-th target data point, r i represents the distance of the i-th target data point, and θ i represents the azimuth angle of the i-th target data point.

[0118] Step 403: According to the plane rectangular coordinates of each target data point in the plurality of target data points, the arc length value of the position of each target data point in the plurality of target data points is calculated.

[0119] Specifically, the arc length value of the position of each target data point in the plurality of target data points is calculated according to the following formula:

[0120]

[0121] Wherein, i takes the value of [1, n-1], s i represents the arc length value of the position of the i-th target data point, x i represents the horizontal coordinate of the i-th target data point, y i represents the vertical coordinate of the i-th target data point, x i+1 represents the horizontal coordinate of the i+1-th target data point, and y i+1 represents the vertical coordinate of the i+1-th target data point.

[0122]

[0123] Wherein, s n represents the arc length value of the position of the n-th target data point, x1 represents the horizontal coordinate of the 1st target data point, y1 represents the vertical coordinate of the 1st target data point, x n represents the horizontal coordinate of the n-th target data point, and y n represents the vertical coordinate of the n-th target data point.

[0124] That is, it is divided into the 1st to n-1st data points, and the n-th data point, and the two kinds of data points are calculated by different calculation objects, because the last data point is adjacent to the first data point.

[0125] Step 404: determining the offset of the well shaft centering of the ultrasonic detector according to the arc length values of the positions of the plurality of target data points;

[0126] Step 405: correcting the plane rectangular coordinates of each target data point in the plurality of target data points according to the offset.

[0127] In the above step of determining the offset of the well shaft centering of the ultrasonic detector according to the arc length values of the positions of the plurality of target data points, the following steps can be included:

[0128] S1: selecting the minimum arc length value from the arc length values of the positions of the plurality of target data points as a reference value;

[0129] S2: calculating the ratio between the arc length value of the position of each target data point in the plurality of target data points and the reference value, and taking the calculation result as the weighting factor of each target data point;

[0130] S3: calculating the geometric center coordinates of the circumferential curve according to the weighting factor of each target data point;

[0131] Specifically, the geometric center coordinates of the circumferential curve can be calculated according to the following formula:

[0132]

[0133]

[0134] wherein X represents the horizontal coordinate of the geometric center coordinates, Y represents the vertical coordinate of the geometric center coordinates, w i represents the weighting factor of the i-th target data point, w sum represents the cumulative value of the weighting factors of the target data points, x i represents the horizontal coordinate of the i-th target data point, y i represents the vertical coordinate of the i-th target data point.

[0135] S4: obtaining the position coordinates of the ultrasonic detector as the reference origin coordinates;

[0136] S5: comparing the geometric center coordinates with the reference origin coordinates to determine the offset of the well shaft centering of the ultrasonic detector.

[0137] The ultrasonic detector described above can be applied in the detection of casing deformation. By correcting the offset of the well shaft centering of the ultrasonic detector, the accuracy of the determination of casing deformation can be effectively improved.

[0138] The following calculates the centering correction for different shapes of the sleeve, and uses a simulated calculation to detect the relative geometric relationship of the data points to obtain the coordinate information of the data points, including the azimuth angle of the data points and the distance of the data points from the instrument. In the calculation, the number of data points is fixed at 32, i.e. 32 data points are selected per week, and the included angle between the data points is 11.25°.

[0139] 1) The shape of the sleeve is a standard circle:

[0140] Assuming that the radius of the circle is 1, the position of the detection instrument is offset by 0.5 in the negative direction of the horizontal axis and by 0.5 in the negative direction of the vertical axis compared to the center of the centering point. As shown in FIG. 1, in the process of three-dimensional imaging of ultrasonic detection, the detection instrument rotates around itself as the center of rotation and scans along the circumferential direction to capture the information of each detection data point, which includes the azimuth angle of the data point and the distance. According to the geometric relationship of the rotation scanning of the detector in FIG. 1, the information of each data point is solved, and the detection instrument rotates at a uniform speed. Figure 5 Figure 5

[0141] Figure 6 The information of 32 data points is described, i.e. the distribution curve of the rotation angle of different data points and the distance from the detector. The rotation angle and the distance of the data points in FIG. 1 are converted into the rectangular coordinate system with the detection instrument as the coordinate origin, and the plane rectangular coordinates of each data point can be obtained. According to the above formula 1 to formula 13, the offset distance is calculated, and the results are as follows: the horizontal axis offset distance is 0.516, and the vertical axis offset distance is 0.473. Compared with the actual horizontal axis offset distance of 0.5 and the actual vertical axis offset distance of 0.5, the correction error of the horizontal axis offset distance is 3.3%, and the correction error of the vertical axis offset distance is 5.4%. Figure 6 The relative position relationship of the sleeve shape before and after correction is shown, where the origin coordinate (0, 0) point is the position of the instrument. Figure 7

[0142] 2) The shape of the sleeve is a standard ellipse:

[0143] Assuming that the semi-major axis of the ellipse is 1 and the semi-minor axis is 0.5, the position of the detection instrument is offset by 0.25 in the negative direction of the horizontal axis and by 0.25 in the negative direction of the vertical axis compared to the centering point. As shown in FIG. 2, in the process of three-dimensional imaging of ultrasonic detection, the detection instrument rotates around itself as the center of rotation and scans along the circumferential direction to capture the information of each detection data point, which includes the azimuth angle of the data point and the distance. Similarly, according to the geometric relationship of the rotation scanning of the detector in FIG. 2, the information of each data point is solved, and the detection instrument rotates at a uniform speed. Figure 8 Figure 8

[0144] Figure 9 ​​​​​Information on 32 data points is described, namely the rotation angle of different data points and the distance distribution curve of the distance detector. Figure 9 The rotation angles and distances of the data points are transformed into a Cartesian coordinate system with the testing instrument as the origin, yielding the Cartesian coordinates of each data point. Using formulas 1 to 13, the offset is calculated, with the following results: the horizontal offset is 0.245, and the vertical offset is 0.236. Compared to the actual horizontal and vertical offsets of 0.25, the correction error for the horizontal offset is 2.0%, and the correction error for the vertical offset is 5.5%. Figure 10 The relative positional relationship of the sleeve shape before and after calibration is shown, where the origin coordinates (0, 0) represent the location of the instrument.

[0145] 3) The sleeve has an irregular shape:

[0146] To facilitate comparison with the actual offset, the irregular shape is set to be symmetrical about the origin, so its geometric center is the origin. The intercept of this shape with the horizontal axis is 1, and the intercept with the vertical axis is also 1. Compared to the center point, the position of the detection instrument is offset by 0.5 in the negative direction of the horizontal axis and 0.5 in the negative direction of the vertical axis. Figure 11 As shown, in the ultrasonic 3D imaging process, the instrument rotates circumferentially around itself as the center of rotation, capturing information from each data point. This information includes the azimuth angle and distance of the data points located on the irregularly shaped circumference. Similarly, according to... Figure 11 The geometric relationship of the rotating scan of the detector is used to solve the information of each data point. The detector rotates at a constant speed.

[0147] Figure 12 Information on 32 data points is described, namely the rotation angle of different data points and the distance distribution curve of the distance detector. Figure 12 The rotation angles and distances of the data points are transformed into a Cartesian coordinate system with the testing instrument as the origin, yielding the Cartesian coordinates of each data point. The offset is calculated using formulas 1 to 13, with the following results: the horizontal offset is 0.509, and the vertical offset is 0.473. Compared to the actual horizontal and vertical offsets of 0.5, the correction error for the horizontal offset is 1.8%, and the correction error for the vertical offset is 5.4%. Figure 13 The relative positional relationship of the sleeve shape before and after calibration is shown, where the origin coordinates (0, 0) represent the location of the instrument.

[0148] The three-dimensional imaging centering degree of the standard circle is analyzed. In the three-dimensional checking calculation, the parameters to be set include the number of data points of a circle, the number of turns of the spiral line, the offset amount when the detection instrument is offset, the longitudinal position difference between adjacent data points, etc. It is assumed that the casing is a standard circle. In the simulated well section, the two ends each account for 30% of the total well length. It is assumed that the detection instrument is always in the center position of the well axis in the two sections. The middle section accounts for 40% of the total well length. In the middle section, the detection instrument is eccentric. The data point information is simulated by setting the eccentricity parameter. It is assumed that the standard circle radius is 1, the maximum eccentricity of the middle section is 0.4, the number of data points in one turn is 100, a total of 150 turns of the spiral line (i.e., the spiral line rotates 150 turns), and the longitudinal coordinate difference between adjacent data points is 0.00024. After the checking, the data points on the casing circle in the eccentric well section are still unevenly distributed, but the three-dimensional shape formed by the data points is already in the center position, and can better reflect the true shape of the casing. As can be seen from the three-dimensional graph and the top view, the data checking result at the junction of the eccentric and the center has a certain error. According to the quantitative calculation result, the maximum checking calculation error is 3% when 100 data points are selected in one turn, which can meet the engineering needs.

[0149] The three-dimensional imaging centering degree of the standard circle is analyzed. In the three-dimensional checking calculation, the parameters to be set include the number of data points of a circle, the number of turns of the spiral line, the offset amount when the detection instrument is offset, the longitudinal position difference between adjacent data points, etc. It is assumed that the casing is a standard circle. In the simulated well section, the two ends each account for 30% of the total well length. It is assumed that the detection instrument is always in the center position of the well axis in the two sections. The middle section accounts for 40% of the total well length. In the middle section, the detection instrument is eccentric. The data point information is simulated by setting the eccentricity parameter. It is assumed that the standard circle radius is 1, the maximum eccentricity of the middle section is 0.4, the number of data points in one turn is 100, a total of 150 turns of the spiral line (i.e., the spiral line rotates 150 turns), and the longitudinal coordinate difference between adjacent data points is 0.00024. After the checking, the data points on the casing circle in the eccentric well section are still unevenly distributed, but the three-dimensional shape formed by the data points is already in the center position, and can better reflect the true shape of the casing. As can be seen from the three-dimensional graph and the top view, the data checking result at the junction of the eccentric and the center has a certain error. According to the quantitative calculation result, the maximum checking calculation error is 3% when 100 data points are selected in one turn, which can meet the engineering needs.

[0150] The method provided in the above embodiments of the application can be executed in a mobile terminal, a computer terminal or a similar computing device. Taking the running on an electronic device as an example, Figure 14 is a hardware structure block diagram of an electronic device for a centering degree correction method of an ultrasonic three-dimensional imaging well axis provided by the application. As Figure 14 shown, the electronic device 10 can include one or more (only one is shown in the figure) processors 02 (the processor 02 can include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA, etc.), a memory 04 for storing data, and a transmission module 06 for communication function. Those skilled in the art can understand, Figure 14The structure shown is only schematic and does not limit the structure of the electronic device. For example, the electronic device 10 can further include more or fewer components than those shown, or have a different configuration of components than those shown. Figure 14 The structure shown is only schematic and does not limit the structure of the electronic device. For example, the electronic device 10 can further include more or fewer components than those shown, or have a different configuration of components than those shown. Figure 14 The structure shown is only schematic and does not limit the structure of the electronic device. For example, the electronic device 10 can further include more or fewer components than those shown, or have a different configuration of components than those shown.

[0151] The memory 04 can be used to store software programs and modules of application software, such as the program instructions / modules of the well axis centricity correction method for ultrasonic three-dimensional imaging in the embodiments of the present application. The processor 02 executes various functional applications and data processing by running the software programs and modules stored in the memory 04, i.e., implements the well axis centricity correction method for ultrasonic three-dimensional imaging of the application program. The memory 04 can include a high-speed random access memory, and can further include a non-volatile memory, such as one or more magnetic storage devices, flash memories, or other non-volatile solid-state memories. In some examples, the memory 04 can further include a memory remotely arranged with respect to the processor 02, which can be connected to the electronic device 10 through a network. Examples of the network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and a combination thereof.

[0152] The transmission module 06 is used to receive or send data via a network. Specific examples of the network can include a wireless network provided by a communication provider of the electronic device 10. In one example, the transmission module 06 includes a network adapter (NIC), which can be connected to other network devices through a base station so as to communicate with the Internet. In one example, the transmission module 06 can be a radio frequency (RF) module, which is used to communicate with the Internet in a wireless manner.

[0153] In the software layer, the device can include, as shown in Figure 15

[0154] The scanning module 1501 is used to capture information of a plurality of target data points on a target casing by an ultrasonic detector along a circumferential scanning;

[0155] The calculation module 1502 is used to calculate a planar rectangular coordinate of each target data point in the plurality of target data points according to the information of the target data points;

[0156] The calculation module 1502 is used to calculate a planar rectangular coordinate of each target data point in the plurality of target data points according to the information of the target data points;

[0157] ​The determining module 1504 is configured to determine an offset of the well axis of the ultrasonic detector from the center according to the arc length value of the position of the target data points.

[0158] The correcting module 1505 is configured to correct the plane rectangular coordinates of each target data point according to the offset.

[0159] In an embodiment, the information of the target data points can include the azimuth angle and the distance of the target data points, and the calculating module 1502 can convert the polar coordinates represented by the azimuth angle and the distance of the target data points into the plane rectangular coordinates represented by the horizontal coordinate and the vertical coordinate according to the following formula:

[0160] x i = r i × cos (θ i )

[0161] y i = r i × sin (θ i )

[0162] wherein x i represents the horizontal coordinate of the i-th target data point, y i+1 represents the vertical coordinate of the i-th target data point, r i+1 represents the distance of the i-th target data point, and θ n represents the azimuth angle of the i-th target data point.

[0163] In an embodiment, the obtaining module 1503 can calculate the arc length value of the position of each target data point according to the following formula:

[0164]

[0165] wherein i is an integer between 1 and n-1, s n represents the arc length value of the position of the i-th target data point, x n represents the horizontal coordinate of the i-th target data point, y i represents the vertical coordinate of the i-th target data point, x sum represents the horizontal coordinate of the i+1-th target data point, and y i represents the vertical coordinate of the i+1-th target data point.

[0166]

[0167] wherein s iLet x1 represent the arc length of the nth target data point, x1 represent the x-coordinate of the 1st target data point, and y1 represent the y-coordinate of the 1st target data point. n The x-coordinate of the nth target data point is represented by y. n This represents the ordinate of the nth target data point.

[0168] In one embodiment, the determining module 1504 may include: a selection unit, configured to select the minimum arc length value from the arc length values ​​of the positions of the plurality of target data points as a reference value; a first calculation unit, configured to calculate the ratio between the arc length value of each target data point and the reference value, and use the calculation result as a weighting factor for each target data point; a second calculation unit, configured to calculate the geometric center coordinates of the circumferential curve based on the weighting factors of each target data point; an acquisition module, configured to acquire the position coordinates of the ultrasonic detector as the reference origin coordinates; and a determining unit, configured to compare the geometric center coordinates with the reference origin coordinates to determine the offset of the ultrasonic detector's well axis centering.

[0169] In one implementation, calculating the geometric center coordinates of the circular curve based on the weighting factors of each target data point may include:

[0170] Calculate the coordinates of the geometric center of the circular curve using the following formula:

[0171]

[0172]

[0173] Where X represents the x-coordinate of the geometric center, Y represents the y-coordinate of the geometric center, and w i w represents the weighting factor for the i-th target data point. sum x represents the cumulative value of the weighting factors for each target data point. i The x-coordinate of the i-th target data point is represented by y. i This represents the ordinate of the i-th target data point.

[0174] This application also provides a specific implementation of an electronic device capable of implementing all steps of the ultrasonic three-dimensional imaging well-axis centering correction method in the above embodiments. The electronic device specifically includes: a processor, a memory, a communication interface, and a bus; wherein the processor, memory, and communication interface communicate with each other via the bus; the processor is used to call a computer program in the memory, and when the processor executes the computer program, it implements all steps of the ultrasonic three-dimensional imaging well-axis centering correction method in the above embodiments. For example, when the processor executes the computer program, it implements the following steps:

[0175] Step 1: Capture information from multiple target data points on the target sleeve by scanning circumferentially with an ultrasonic detector;

[0176] Step 2: Using the information of the target data points, calculate the Cartesian coordinates of each target data point among the plurality of target data points;

[0177] Step 3: Based on the Cartesian coordinates of each target data point among the multiple target data points, calculate the arc length of the position of each target data point among the multiple target data points;

[0178] Step 4: Determine the offset of the ultrasonic detector's well axis centering based on the arc length values ​​of the multiple target data points.

[0179] Step 5: Correct the Cartesian coordinates of each target data point among the plurality of target data points according to the offset.

[0180] Embodiments of this application also provide a computer-readable storage medium capable of implementing all steps of the ultrasonic three-dimensional imaging well-axis centering correction method in the above embodiments. The computer-readable storage medium stores a computer program that, when executed by a processor, implements all steps of the ultrasonic three-dimensional imaging well-axis centering correction method in the above embodiments. For example, when the processor executes the computer program, it implements the following steps:

[0181] Step 1: Capture information from multiple target data points on the target sleeve by scanning circumferentially with an ultrasonic detector;

[0182] Step 2: Using the information of the target data points, calculate the Cartesian coordinates of each target data point among the plurality of target data points;

[0183] Step 3: Based on the Cartesian coordinates of each target data point among the multiple target data points, calculate the arc length of the position of each target data point among the multiple target data points;

[0184] Step 4: Determine the offset of the ultrasonic detector's well axis centering based on the arc length values ​​of the multiple target data points.

[0185] Step 5: Correct the Cartesian coordinates of each target data point among the plurality of target data points according to the offset.

[0186] As described above, this embodiment of the application captures information of multiple target data points on the target casing by scanning circumferentially with an ultrasonic detector; calculates the Cartesian coordinates of each target data point based on the information of the target data points; calculates the arc length of the position of each target data point based on the Cartesian coordinates of the target data points; determines the offset of the ultrasonic detector's well axis centering based on the arc length of the position of the target data points; and corrects the Cartesian coordinates of each target data point based on the offset. That is, the eccentricity of the detection instrument is determined by using the arc length difference between data points, and the calculated eccentricity can be used to verify the well axis centering of the ultrasonic three-dimensional imaging of casing deformation, thereby solving the problem of low three-dimensional imaging accuracy caused by the inability to verify the well axis centering in existing methods, and achieving a significant improvement in the technical effect of three-dimensional imaging accuracy.

[0187] The user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties. Furthermore, the collection, use and processing of the relevant data must comply with the relevant laws, regulations and standards of the relevant countries and regions, and corresponding operation entry points are provided for users to choose to authorize or refuse.

[0188] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to interchangeably. Each embodiment focuses on its differences from other embodiments. In particular, hardware + program embodiments are relatively simple in description because they are fundamentally similar to method embodiments; relevant parts can be referred to the descriptions in the method embodiments.

[0189] The foregoing has described specific embodiments of this specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims may be performed in a different order than that shown in the embodiments and may still achieve the desired result. Furthermore, the processes depicted in the drawings do not necessarily require the specific or sequential order shown to achieve the desired result. In some embodiments, multitasking and parallel processing are possible or may be advantageous.

[0190] While this application provides the method operation steps as described in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-inventive labor. The order of steps listed in the embodiments is merely one possible execution order among many and does not represent the only execution order. In actual device or client product execution, the methods shown in the embodiments or drawings can be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment).

[0191] The systems, devices, modules, or units described in the above embodiments can be implemented by computer chips or entities, or by products with certain functions. A typical implementation device is a computer. Specifically, a computer can be, for example, a personal computer, a laptop computer, an in-vehicle human-machine interaction device, a cellular phone, a camera phone, a smartphone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or any combination of these devices.

[0192] While this specification provides method operation steps as described in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-inventive means. The order of steps listed in the embodiments is merely one possible execution order among many and does not represent the only execution order. In actual device or end product execution, the methods shown in the embodiments or drawings may be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment, or even a distributed data processing environment). The terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, product, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, product, or apparatus. Without further limitations, the presence of other identical or equivalent elements in the process, method, product, or apparatus that includes said elements is not excluded.

[0193] For ease of description, the above devices are described in terms of function, divided into various modules. Of course, in implementing the embodiments of this specification, the functions of each module can be implemented in one or more software and / or hardware components, or a module that performs the same function can be implemented by a combination of multiple sub-modules or sub-units. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division; in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces, or indirect coupling or communication connection between devices or units, and may be electrical, mechanical, or other forms.

[0194] Those skilled in the art will also know that, besides implementing the controller using purely computer-readable program code, the same functions can be achieved by logically programming the method steps, making the controller function as logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers (PLCs), and embedded microcontrollers. Therefore, such a controller can be considered a hardware component, and the devices within it used to implement various functions can also be considered structures within that hardware component. Alternatively, the devices used to implement various functions can be considered as both software modules implementing the method and structures within a hardware component.

[0195] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0196] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0197] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0198] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0199] Memory may include non-persistent storage in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.

[0200] Computer-readable media includes both permanent and non-permanent, removable and non-removable media that can store information using any method or technology. Information can be computer-readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, disk storage or other magnetic storage devices, or any other non-transferable medium that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transient computer-readable media, such as modulated data signals and carrier waves.

[0201] Those skilled in the art will understand that the embodiments of this specification can be provided as methods, systems, or computer program products. Therefore, the embodiments of this specification can take the form of entirely hardware embodiments, entirely software embodiments, or embodiments combining software and hardware aspects. Furthermore, the embodiments of this specification can take the form of computer program products implemented on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0202] The embodiments described in this specification can be described in the general context of computer-executable instructions, such as program modules, that are executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform a specific task or implement a specific abstract data type. The embodiments of this specification can also be practiced in distributed computing environments where tasks are performed by remote processing devices connected via a communication network. In distributed computing environments, program modules can reside in local and remote computer storage media, including storage devices.

[0203] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments. In particular, system embodiments are basically similar to method embodiments, so the description is relatively simple; relevant parts can be referred to the descriptions in the method embodiments. In the description of this specification, the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of the embodiments in this specification. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described can be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0204] The above description is merely an embodiment of the present specification and is not intended to limit the embodiments of the present specification. For those skilled in the art, various modifications and variations can be made to the embodiments of the present specification. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principle of the embodiments of the present specification should be included within the scope of the claims of the embodiments of the present specification.

Claims

1. A method for correcting the centering of a well axis in ultrasonic three-dimensional imaging, characterized in that, The method includes: Information on multiple target data points on the target sleeve is captured by scanning along the circumferential direction using an ultrasonic detector. Using the information of the target data points, the Cartesian coordinates of each target data point among the plurality of target data points are calculated; Based on the Cartesian coordinates of each target data point among the plurality of target data points, the arc length value of the position of each target data point among the plurality of target data points is obtained; Based on the arc length values ​​of the locations of the multiple target data points, the offset of the ultrasonic detector's well axis is determined. Based on the offset, the Cartesian coordinates of each of the plurality of target data points are corrected. The method for determining the offset of the ultrasonic detector's well axis centering based on the arc length values ​​of the multiple target data points includes: The minimum arc length value is selected from the arc length values ​​of the multiple target data points as the reference value; Calculate the ratio between the arc length of each target data point and the reference value, and use the calculation result as a weighting factor for each target data point; Calculate the geometric center coordinates of the circular curve based on the weighting factors of each target data point; The coordinates of the ultrasonic detector's location are used as the reference origin coordinates. The coordinates of the geometric center are compared with the coordinates of the reference origin to determine the offset of the ultrasonic detector's well axis centering.

2. The method according to claim 1, characterized in that, The information of the plurality of target data points includes: the azimuth angle and distance of the plurality of target data points; Accordingly, using the information of the target data points, the Cartesian coordinates of each target data point among the plurality of target data points are calculated, including: The polar coordinates, represented by the azimuth angles and distances of the multiple target data points, are converted into Cartesian coordinates, represented by the horizontal and vertical coordinates, according to the following formula. x i =r i ×cos(θ i ) y i =r i ×sin(θ i ) Where, x i The x-coordinate of the i-th target data point is represented by y. i Let r represent the ordinate of the i-th target data point. i θ represents the distance to the i-th target data point. i This represents the azimuth angle of the i-th target data point.

3. The method according to claim 1, characterized in that, Based on the Cartesian coordinates of each target data point among the plurality of target data points, the arc length value of the position of each target data point among the plurality of target data points is calculated, including: Calculate the arc length of each target data point in a set of multiple target data points using the following formula: Where i takes the value [1, n-1], s i x represents the arc length of the location of the i-th target data point. i The x-coordinate of the i-th target data point is represented by y. i Let x represent the ordinate of the i-th target data point. i+1 The x-coordinate of the (i+1)th target data point is represented by y. i+1 Represents the ordinate of the (i+1)th target data point; Among them, s n Let x1 represent the arc length of the nth target data point, x1 represent the x-coordinate of the 1st target data point, and y1 represent the y-coordinate of the 1st target data point. n The x-coordinate of the nth target data point is represented by y. n This represents the ordinate of the nth target data point.

4. The method according to claim 1, characterized in that, Based on the weighting factors for each target data point, calculate the coordinates of the geometric center of the circular curve, including: Calculate the coordinates of the geometric center of the circular curve using the following formula: Where X represents the x-coordinate of the geometric center, Y represents the y-coordinate of the geometric center, and w i w represents the weighting factor for the i-th target data point. sum x represents the cumulative value of the weighting factors for each target data point. i The x-coordinate of the i-th target data point is represented by y. i This represents the ordinate of the i-th target data point.

5. The method according to any one of claims 1 to 4, characterized in that, The ultrasonic testing instrument is used for detecting sleeve deformation.

6. A device for correcting the centering of a well axis in ultrasonic three-dimensional imaging, characterized in that, include: The scanning module is used to capture information of multiple target data points on the target sleeve by scanning along the circumferential direction with an ultrasonic detector; The calculation module is used to calculate the Cartesian coordinates of each target data point among the plurality of target data points using the information of the target data points; The calculation module is used to calculate the arc length value of the position of each target data point in the plurality of target data points based on the Cartesian coordinates of each target data point in the plurality of target data points; The determination module is used to determine the offset of the ultrasonic detector's well axis centering based on the arc length values ​​of the positions of the multiple target data points; The correction module is used to correct the Cartesian coordinates of each of the plurality of target data points according to the offset. The determining module includes: The selection unit is used to select the minimum arc length value from the arc length values ​​of the positions of the plurality of target data points as a reference value; The first calculation unit is used to calculate the ratio between the arc length value of the position of each target data point among the plurality of target data points and the reference value, and use the calculation result as a weighting factor for each target data point; The second calculation unit is used to calculate the geometric center coordinates of the circular curve based on the weighting factors of each target data point. The acquisition module is used to acquire the coordinates of the ultrasonic detector's location as the reference origin coordinates; A determining unit is used to compare the coordinates of the geometric center with the coordinates of the reference origin to determine the offset of the well axis of the ultrasonic detector.

7. An electronic device comprising a processor and a memory for storing processor-executable instructions, characterized in that, When the processor executes the instructions, it implements the steps of the method according to any one of claims 1 to 5.

8. A computer-readable storage medium having a computer program / instructions stored thereon, characterized in that, When the computer program / instructions are executed by the processor, they implement the steps of the method according to any one of claims 1 to 5.

Images