Method and apparatus for correcting well axis centering degree in ultrasonic three-dimensional imaging

Abstract

Provided in the present application are a method and apparatus for correcting a well axis centering degree in ultrasonic three-dimensional imaging. The method comprises: an ultrasonic detector performing scanning in a circumferential direction, so as to capture information of a plurality of target data points on a target casing; performing calculation by means of the information of the target data points, so as to obtain rectangular plane coordinates of each target data point among the plurality of target data points; on the basis of the rectangular plane coordinates of each target data point among the plurality of target data points, solving for an arc length value of each position where each target data point among the plurality of target data points is located; on the basis of the arc length values of the positions where the plurality of target data points are located, determining an offset of the position, where the ultrasonic detector is located, with respect to well axis centering; and on the basis of the offset, correcting the rectangular plane coordinates of each target data point among the plurality of target data points. By means of the solution, the existing problem of three-dimensional imaging having relatively low accuracy due to the fact that a well axis centering degree cannot be checked is solved, thereby achieving the technical effect of effectively improving the accuracy of three-dimensional imaging.

Description

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

[0001] Related applications

[0002] This application claims priority to Chinese patent application No. 202311125616.7 filed on September 1, 2023, and cites the contents disclosed in the above patent application as part of this application. Technical Field

[0003] The present application belongs to the field of geological exploration technology, and in particular relates to a method and device for correcting the well axis centering of ultrasonic three-dimensional imaging. Background Art

[0004] Casing completion is the primary completion method during well construction. Factors such as cementing quality, formation movement, wellbore age, hydraulic fracturing stimulation, injection-production measures, alternating injection-production cycles in gas storage facilities, and fluid properties can all affect the working condition of the casing, leading to casing deformation and posing a serious threat to oil and gas production safety. Current production practices demonstrate that the prevention, detection, and management of casing deformation are crucial for safe production.

[0005] Casing deformation detection is a prerequisite for evaluating and repairing casing deformation, as well as for other production improvement measures. Currently, the main detection technologies for casing deformation include multi-arm caliper testing, electromagnetic flaw detection, and ultrasonic 3D imaging. Among these detection technologies, ultrasonic 3D imaging boasts excellent directionality, strong penetration, and high measurement accuracy. It enables precise measurement and characterization of casing deformation, effectively addressing issues such as low wellbore coverage and difficulties in ensuring measurement accuracy and resolution. The principle of ultrasonic 3D imaging for casing deformation detection is as follows: an ultrasonic probe is mounted at the bottom of the detection instrument. During the downhole inspection process, the instrument ascends or descends at a predetermined speed. The ultrasonic probe rotates at high speed as the instrument ascends or descends, emitting ultrasonic pulse signals toward the wellbore. These pulse signals propagate through the wellbore fluid to the inner wall of the casing, where they are reflected back. The final echo signal is then received by the ultrasonic probe, capturing data at each sampling point through a self-transmitting and self-receiving process. Based on the operating speed of the detection instrument, the rotation speed of the ultrasonic probe, and the time required for the ultrasonic signal to propagate, the time interval for the excitation of the ultrasonic pulse signal is optimized, and then the pulse signal can be continuously excited according to a regular pattern, thereby capturing data from a large number of sampling points. Based on the time difference of the echo signal at each sampling point, the distance between the instrument and the inner wall of the casing can be calculated. Combined with the rotation angle of the ultrasonic probe and the azimuth correction of the gyroscope, the position coordinates of each measuring point can be obtained. By presenting these measuring point coordinates in three-dimensional space, the three-dimensional characteristics of the casing deformation morphology can be depicted. Under ideal conditions, the detection instrument will be stably located at the center of the wellbore axis with the support of the stabilizer. In this way, the data captured by the detection will have a unified center position at different well depths, and the relative position relationships presented can all be referenced with the unified center position.

[0006] However, in actual inspection conditions, especially in long horizontal sections, the centralizer cannot ensure the stable centering of the inspection instrument. The instrument will erratically deviate circumferentially, causing the well axis reference to vary at different well depths, making it difficult to accurately determine relative position. Currently, no effective solution has been proposed for verifying the well axis centering in ultrasonic 3D imaging of casing deformation to improve the accuracy of 3D imaging and better visualize the true form of casing deformation.

[0007] Summary of the Invention

[0008] The purpose of this application is to provide a method and device for correcting the well axis centering of ultrasonic three-dimensional imaging, which can correct the deviation caused by the deviation of the well axis direction, thereby improving the accuracy of three-dimensional imaging.

[0009] The present application provides a method and apparatus for correcting well axis centering in ultrasonic three-dimensional imaging, which is implemented as follows:

[0010] A method for correcting well axis centering in ultrasonic three-dimensional imaging, comprising:

[0011] The ultrasonic detector is used to scan along the circumference to capture information of multiple target data points on the target casing;

[0012] Calculating the plane rectangular coordinates of each target data point among the multiple target data points based on the information of the target data point;

[0013] Calculating the arc length value of each target data point in the plurality of target data points according to the plane rectangular coordinates of each target data point; the arc length value refers to the length of the arc path formed by the data point during the rotation scanning process of the detector;

[0014] Determining the offset of the ultrasonic detector from the center of the well axis according to the arc length values ​​of the positions of the multiple target data points, that is, the offset of the position of the ultrasonic detector relative to the center of the well axis;

[0015] The plane rectangular coordinates of each target data point among the multiple target data points are corrected according to the offset.

[0016] In one embodiment, the information of the plurality of target data points includes: azimuth angles and distances of the plurality of target data points;

[0017] Accordingly, the plane rectangular coordinates of each target data point among the multiple target data points are calculated based on the information of the target data point, including:

[0018] The polar coordinates represented by the azimuth angles and distances of multiple target data points are converted into plane rectangular coordinates represented by horizontal and vertical coordinates according to the following formula; i =r i ×cos(θ i ) y i =r i ×sin(θ i )

[0019] Among them, x i Indicates the horizontal coordinate of the i-th target data point, y i Represents the ordinate of the i-th target data point, r i represents the distance to the i-th target data point, θ i Represents the azimuth angle of the i-th target data point.

[0020] In one embodiment, obtaining the arc length value of each target data point among the multiple target data points according to the plane rectangular coordinates of each target data point includes:

[0021] Calculate the arc length of each target data point among multiple target data points according to the following formula:

[0022] Among them, i is [1, n-1], s i Indicates the arc length value of the location of the i-th target data point, x i Indicates the horizontal coordinate of the i-th target data point, y i Indicates the vertical coordinate of the i-th target data point, x i+1 Indicates the horizontal coordinate of the i+1th target data point, y i+1 Indicates the ordinate of the i+1th target data point;

[0023] Among them, s n Indicates the arc length value of the location of the nth target data point, x1 indicates the horizontal coordinate of the first target data point, y1 indicates the vertical coordinate of the first target data point, x n Indicates the horizontal coordinate of the nth target data point, y n Indicates the vertical coordinate of the nth target data point.

[0024] In one embodiment, determining the offset of the ultrasonic detector relative to the center of the well axis based on the arc length values ​​of the locations of the plurality of target data points includes:

[0025] Select the minimum arc length value from the arc length values ​​of the locations of multiple target data points as the reference value;

[0026] Calculating the ratio between the arc length value of each target data point and the reference value, and using the calculation result as a weighting factor of each target data point;

[0027] Calculate the geometric center coordinates of the circular curve according to the weighting factors of each target data point;

[0028] Obtain the coordinates of the ultrasonic detector's location as the reference origin coordinates;

[0029] The geometric center coordinates are compared with the reference origin coordinates to determine the offset of the ultrasonic detector position relative to the center of the well axis.

[0030] In one embodiment, calculating the geometric center coordinates of the circular curve according to the weighting factors of each target data point includes:

[0031] The geometric center coordinates of the circular curve are calculated according to the following formula:

[0032] Among them, X represents the horizontal coordinate of the geometric center coordinate, Y represents the vertical coordinate of the geometric center coordinate, and wi represents the weighting factor of the i-th target data point, w sum Represents the cumulative value of the weighting factor of each target data point, x i Indicates the horizontal coordinate of the i-th target data point, y i Represents the vertical coordinate of the i-th target data point.

[0033] In one embodiment, an ultrasonic detector is used to detect casing deformation.

[0034] A device for correcting the well axis centering of ultrasonic three-dimensional imaging, comprising:

[0035] A scanning module, used for capturing information of multiple target data points on the target casing by scanning along the circumferential direction using an ultrasonic detector;

[0036] A calculation module, configured to calculate the plane rectangular coordinates of each target data point among the plurality of target data points based on the information of the target data point;

[0037] an obtaining module, for obtaining an arc length value of a position of each target data point among the multiple target data points according to the plane rectangular coordinates of each target data point among the multiple target data points;

[0038] a determination module for determining an offset of the ultrasonic detector relative to the center of the well axis based on arc length values ​​of the locations of the plurality of target data points;

[0039] The correction module is used to correct the plane rectangular coordinates of each target data point among the multiple target data points according to the offset.

[0040] In one embodiment, the determining module includes:

[0041] A selection unit is used to select a minimum arc length value from arc length values ​​of the positions of multiple target data points as a reference value;

[0042] a first calculation unit, configured to calculate a ratio between an arc length value of each target data point among a plurality of target data points and a reference value, and use the calculation result as a weighting factor for each target data point;

[0043] A second calculation unit is used to calculate the geometric center coordinates of the circular curve according to the weighting factors of each target data point;

[0044] An acquisition module is used to obtain the coordinates of the position of the ultrasonic detector as the reference origin coordinates;

[0045] The determination unit is used to compare the geometric center coordinates with the reference origin coordinates to determine the offset of the ultrasonic detector position relative to the center of the well axis.

[0046] An electronic device includes a processor and a memory for storing instructions executable by the processor, and the steps of the above method are implemented when the processor executes the instructions.

[0047] A computer-readable storage medium stores a computer program / instruction thereon, which implements the steps of the above method when executed by a processor.

[0048] The present application provides a method for correcting the well axis centering of ultrasonic three-dimensional imaging. The method comprises: capturing information of multiple target data points on a target casing by scanning the target casing along a circumferential direction with an ultrasonic detector; calculating the plane rectangular coordinates of each target data point in the multiple target data points based on the information of the target data points; determining the arc length value of each target data point in the multiple target data points based on the plane rectangular coordinates of each target data point; determining the offset of the ultrasonic detector position relative to the well axis centering based on the arc length values ​​of the multiple target data points; and correcting the plane rectangular coordinates of each target data point in the multiple target data points based on the offsets. Specifically, the arc length difference between the data points is used to determine the eccentricity of the detector. The calculated eccentricity can be used to verify the well axis centering of ultrasonic three-dimensional imaging of casing deformation, thereby resolving the problem of low three-dimensional imaging accuracy caused by the inability to verify the well axis centering in existing methods, thereby achieving the technical effect of effectively improving the accuracy of three-dimensional imaging. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments recorded in this application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative labor.

[0050] FIG1 is a schematic diagram of the geometric relationship of the casing deformation ultrasonic three-dimensional imaging detection instrument provided by the present application when it is in an eccentric position during operation;

[0051] FIG2 is a schematic diagram of the horizontal and vertical coordinate offset of a standard circular casing provided by the present application;

[0052] FIG3 is a schematic diagram showing the relationship between the calculated values ​​of the horizontal and vertical coordinate offsets of a standard circular casing provided by the present application and the number of data points;

[0053] FIG4 is a flow chart of an embodiment of a method for correcting well axis centering in ultrasonic three-dimensional imaging provided by the present application;

[0054] FIG5 is a schematic diagram of the positional relationship of a standard circular sleeve provided by the present application when scanning and detecting 32 data points;

[0055] FIG6 is a distribution curve of the orientation and distance of different data points when scanning and testing the standard circular casing provided by the present application according to 32 data points;

[0056] FIG7 is a schematic diagram of the relative position relationship of the standard circular sleeve provided by the present application before and after the centering correction;

[0057] FIG8 is a schematic diagram of the positional relationship of the elliptical sleeve provided by the present application when scanning and detecting 32 data points;

[0058] FIG9 is a distribution curve of the orientation and distance of different data points when scanning and testing the elliptical sleeve provided by the present application at 32 data points;

[0059] FIG10 is a schematic diagram of the relative position relationship of the elliptical sleeve provided by the present application before and after the centering correction;

[0060] FIG11 is a schematic diagram of the positional relationship of an irregularly shaped casing provided by the present application when scanned and tested according to 32 data points;

[0061] FIG12 is a distribution curve of the position and distance of different data points when scanning and testing an irregular-shaped casing according to 32 data points provided by the present application;

[0062] FIG13 is a schematic diagram of the relative position relationship of the irregular-shaped sleeve provided by the present application before and after the centering correction;

[0063] FIG14 is a hardware block diagram of an electronic device for a method for correcting well axis centering in ultrasonic three-dimensional imaging provided by the present application;

[0064] FIG15 is a schematic diagram of the module structure of an embodiment of a device for correcting the well axis centering of ultrasonic three-dimensional imaging provided in the present application. DETAILED DESCRIPTION

[0065] In order to enable those skilled in the art to better understand the technical solutions in this application, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of this application.

[0066] To address the technical issues of misalignment of the ultrasonic three-dimensional imaging instrument for casing deformation detection, resulting in inaccurate determination of the relative position of casing images at different well depths and distortion of the true shape of casing deformation, this paper proposes a method for verifying the well axis centering of ultrasonic three-dimensional imaging of casing deformation. Each data point captured by the ultrasonic three-dimensional imaging instrument contains the casing's orientation and distance from the ultrasonic probe at that point. By simplifying the data points around a circle into a plane, the coordinates of that point in the plane coordinate system can be calculated based on the angle and distance of the data point's orientation. Because the ultrasonic probe rotates at a constant speed, the difference in azimuth angle between two adjacent data points is equal. As shown in Figure 1, if the casing is a perfect circle and the instrument is at the center of the circle, the distribution of the scanned data points along the circumference is uniform, meaning that the arc lengths between adjacent data points are equal. If the casing is a perfect circle but the instrument is off-center, the distribution of the scanned data points along the circumference is uneven, with data points closer to the instrument being more densely distributed and data points farther away being more sparsely distributed. The density and sparseness of the data point distribution 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 shorter, and when the data points are sparsely distributed, the arc length between adjacent data points is longer.

[0067] Based on this, the difference in arc length between data points can be used to assign a weighted value to each data point. Based on the weighted value of each data point, the weighted average of the horizontal and vertical coordinates of all data points can be calculated. With a sufficient amount of test data, this weighted average can approximate the geometric center coordinates of the circular curve within the error range. The difference between the calculated geometric center coordinates of the circular curve and the coordinates of the origin is the eccentricity of the testing instrument. The calculated eccentricity can be used to verify the well axis centering of the casing deformation ultrasonic 3D imaging.

[0068] That is, we can first calculate the plane coordinates of each data point based on the data point information captured by the instrument. Taking a circle as an example, we can get the coordinates of the data points for a circle of 360°. The horizontal and vertical coordinates of the i-th data point are x and y respectively. i and y i The information of the i-th data point captured by the instrument includes: azimuth angle (the angle of the i-th data point relative to the instrument) and distance, respectively denoted as θ i and r i , there are n data points in a circle of 360°, so converting polar coordinates into rectangular coordinates, we can get: i =r i ×cos(θ i ) (Formula 1) yi =r i ×sin(θ i ) (Formula 2)

[0069] Among them, i = 1, 2...n. The data points captured by the detector during a scan are unevenly distributed on the circular curve, and the geometric center of the circular curve (representing the center point position) cannot be directly calculated using the average value of the horizontal and vertical coordinates of each data point. The density and sparseness of the data point distribution are directly reflected in the length of the arc between adjacent data points. That is, when the data points are densely distributed, the arc length between adjacent data points is shorter, and when the data points are sparsely distributed, the arc length between adjacent data points is longer. Therefore, the difference in arc length between data points can be used to assign a weighted value to each data point, and the weighted value of each data point can be used to solve the weighted average of the horizontal and vertical coordinates of all data points. When the number of detection data is large enough, the weighted average can approximately represent the geometric center coordinates of the circular curve within the error range. The calculated geometric center coordinates are compared with the origin coordinates referenced by the polar coordinate to rectangular coordinate conversion of Formula 1 and Formula 2 to obtain the offset distance of the current three-dimensional imaging curve detection instrument.

[0070] First, calculate the arc length of each data point. Use the distance between a certain data point and the next data point to approximate the arc length of the data point. Then we can get:

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

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

[0073] The smallest value is selected from all arc lengths. Based on this minimum value, the weighted average of the horizontal and vertical coordinates of the data points is calculated according to the ratio of the arc length of each data point to the minimum value. The minimum value in the arc length can be expressed as: s min =min(s) (Formula 5)

[0074] The weighting factor for each data point is:

[0075] The sum of the weighting factors of all data points is:

[0076] Then the geometric center coordinates of the circular curve are:

[0077] Because when processing the original data point information, the polar coordinates are directly converted to rectangular coordinates, so the position of the detector during the conversion process is the origin. Therefore, the difference between the weighted average coordinates calculated by formula 8 and formula 9 and the origin is the offset, so the offset is: X 偏移距 =X-0 (Formula 10) Y 偏移距 =Y-0 (Formula 11)

[0078] To calibrate each data point from the eccentric position to the central position, the conversion is as follows: i校正值 =x i -X 偏移距 (Formula 12) y i校正值 =y i -Y 偏移距 (Formula 13)

[0079] Among them, Formula 12 and Formula 13 are the corrected data point coordinates.

[0080] The above method is numerically verified. As shown in FIG2 , a standard circular sleeve is shown in FIG2 . The center point is shown in FIG2 (1) and is located at the center of the circle. The eccentric point is shown in FIG2 (2) and is located in the third quadrant of the coordinate system. When the instrument is located at the eccentric point, the path of the transmitted signal is shown in FIG2 (3). These paths constitute a straight line system. This straight line system is composed of straight lines passing through the eccentric point and having different slopes. The coordinates of the intersection of the two can be calculated by the equations of the straight line and the circle. The coordinates of the center 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 coordinates of the eccentric point through which the straight line passes 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, and the intercept of the straight line on the y-axis is b = k×dx-dy. Next, the coordinates of the intersection of the straight line and the circle are calculated.

[0081] First, define several calculation parameters: d a =1+k 2 (Formula 14) d b=2k×b (Formula 15) d c =b 2 -r 2 (Formula 16)

[0082] The straight line and the circle have two intersection points, and the coordinates of the intersection points are: Y1=-k×X1+b (Formula 18) Y2=k×X2+b (Formula 20)

[0083] Assuming there are n data points (n is an even number), that is, n straight lines intersecting the circle, according to the principle of symmetry, only half of the circumference, that is, a rotation line within a range of 180 degrees, is needed to calculate all the data points on the circumference. Then, the horizontal and vertical coordinates of all data points on the circumference are calculated according to Formulas 17 to 20. The distance from the data point to the eccentric point can be solved based on the horizontal and vertical coordinates of each data point. Therefore, the distance from each data point to the detector and the rotational orientation can be calculated. Substituting these into Formulas 1 to 13 above, the offset can be calculated.

[0084] The above formula can be used to generate simulated detection data points of standard circular casing to verify the established calibration model. As shown in Figure 1, assuming that the radius of the standard circle is 1 and the coordinates of the eccentric point are (-0.5, -0.5), the horizontal coordinate offset of the detection instrument is 0.5 and the vertical coordinate offset is 0.5. Then, according to the above-mentioned established model, different numbers of data points are selected to calculate the offset of the detection instrument, and then compared with the actual value of the offset. As shown in Figure 3, the calibration calculation value of the horizontal and vertical coordinate offset changes with the number of data points. As the number of data points increases, the calibration calculation value of the horizontal and vertical coordinate offset gradually approaches the actual value. When the number of data points reaches 202, the calibration calculation error of the horizontal coordinate offset is 0.64%, and the calibration calculation error of the vertical coordinate offset is 0.66%, achieving a high calculation accuracy.

[0085] The above method can be applied to ultrasonic detection of shale gas well casing deformation. Based on the uneven distribution of data points captured along the circumferential scan when the ultrasonic detector is located at an eccentric position, different weight ratios are assigned according to the arc length differences between the data points, and the weighted average of the horizontal and vertical coordinates of all data points is calculated. Based on this weighted average, the geometric center coordinates of the circular curve and the offset of the detection instrument can be calculated. Finally, combined with the geometric change relationship, the coordinates of the data points on the casing around the well are transformed, thereby realizing the data verification of the well axis centering of the ultrasonic three-dimensional imaging of casing deformation. By verifying the centering of the well axis through the above scheme, the offset of the detection instrument can be quickly and accurately calculated, and the relative position relationship of the casing imaging at different well depths can be accurately determined to realize the verification of the well axis centering of the ultrasonic three-dimensional imaging of casing deformation, thereby improving the accuracy of the three-dimensional imaging and better presenting the true shape of the casing deformation.

[0086] FIG4 is a method flow chart of an embodiment of a method for correcting the well axis centering of ultrasonic three-dimensional imaging provided by the present application. Although the present application provides method operation steps or device structures as shown in the following embodiments or drawings, more or fewer operation steps or module units may be included in the method or device based on routine or no creative labor. In the steps or structures where there is no necessary causal relationship logically, the execution order of these steps or the module structure of the device is not limited to the execution order or module structure described in the embodiments of the present application and shown in the drawings. When the method or module structure is applied to an actual device or terminal product, it can be connected according to the method or module structure shown in the embodiments or drawings for sequential execution or parallel execution (for example, a parallel processor or multi-threaded processing environment, or even a distributed processing environment).

[0087] Specifically, as shown in FIG4 , the method for correcting the well axis centering of ultrasonic three-dimensional imaging may include the following steps:

[0088] Step 401: Capturing information of multiple target data points on a target casing by scanning along a circumferential direction using an ultrasonic detector;

[0089] Step 402: Calculate the plane rectangular coordinates of each target data point among the multiple target data points based on the information of the target data point;

[0090] The information of the plurality of target data points may include: azimuth angles and distances of the plurality of target data points;

[0091] Accordingly, the polar coordinates represented by the azimuth angles and distances of multiple target data points can be converted into plane rectangular coordinates represented by the horizontal coordinate and the vertical coordinate according to the following formula; i =r i ×cos(θ i) y i =r i ×sin(θ i )

[0092] Among them, x i Indicates the horizontal coordinate of the i-th target data point, y i Represents the ordinate of the i-th target data point, r i represents the distance to the i-th target data point, θ i Represents the azimuth angle of the i-th target data point.

[0093] Step 403: Calculating the arc length of each target data point in the plurality of target data points based on the plane rectangular coordinates of each target data point; the arc length refers to the length of the arc path formed by the data point during the rotation scanning process of the detector;

[0094] Specifically, the arc length value of each target data point among multiple target data points is calculated according to the following formula:

[0095] Among them, i is [1, n-1], s i Indicates the arc length value of the location of the i-th target data point, x i Indicates the horizontal coordinate of the i-th target data point, y i Indicates the vertical coordinate of the i-th target data point, x i+1 Indicates the horizontal coordinate of the i+1th target data point, y i+1 Indicates the ordinate of the i+1th target data point;

[0096] Among them, s n Indicates the arc length value of the location of the nth target data point, x1 indicates the horizontal coordinate of the first target data point, y1 indicates the vertical coordinate of the first target data point, x n Indicates the horizontal coordinate of the nth target data point, y n Indicates the vertical coordinate of the nth target data point.

[0097] That is, the data points are divided into the 1st to n-1th data points and the nth data point. The two data points use different calculation objects because the last data point is adjacent to the first data point.

[0098] Step 404: determining an offset of the ultrasonic detector relative to the center of the well axis based on the arc length values ​​of the locations of the multiple target data points;

[0099] Step 405: Correcting the plane rectangular coordinates of each target data point among the multiple target data points according to the offset.

[0100] When determining the offset of the ultrasonic detector relative to the center of the well axis based on the arc length values ​​of the positions of the multiple target data points, the following steps may be included:

[0101] S1: Select the minimum arc length value from the arc length values ​​of multiple target data points as the reference value;

[0102] S2: Calculate the ratio between the arc length value of each target data point and the reference value, and use the calculation result as the weighting factor of each target data point;

[0103] S3: Calculate the geometric center coordinates of the circular curve according to the weighting factors of each target data point;

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

[0105] Among them, X represents the horizontal coordinate of the geometric center coordinate, Y represents the vertical coordinate of the geometric center coordinate, and w i represents the weighting factor of the i-th target data point, w sum Represents the cumulative value of the weighting factor of each target data point, x i Indicates the horizontal coordinate of the i-th target data point, y i Represents the vertical coordinate of the i-th target data point.

[0106] S4: Obtain the coordinates of the position of the ultrasonic detector as the reference origin coordinates;

[0107] S5: Compare the geometric center coordinates with the reference origin coordinates to determine the offset of the ultrasonic detector position relative to the center of the well axis.

[0108] The above-mentioned ultrasonic detector can be used in the detection of casing deformation. By correcting the offset of the position of the ultrasonic detector relative to the well axis, the accuracy of the casing deformation determination can be effectively improved.

[0109] The following calculations perform centering correction for casing of varying shapes. Specifically, they determine whether the detector is centered on the well axis during scanning. Simulations are used to calculate the relative geometric relationships of the test data points to obtain their coordinate information, including their azimuth angle and distance from the instrument. The calculations are performed with a fixed number of 32 data points, i.e., 32 points are selected per circle, with an angle of 11.25° between each point.

[0110] 1) The shape of the casing is a standard circle:

[0111] Assuming the radius of the circle is 1, the position of the testing instrument is offset by 0.5 in the negative horizontal direction and 0.5 in the negative vertical direction relative to the center of the circle. As shown in Figure 5, during the ultrasonic 3D imaging process, the testing instrument rotates and scans around itself, capturing information about each test data point, including its azimuth angle and distance. Based on the geometric relationship of the detector's rotational scanning in Figure 5, the information of each data point is solved, and the testing instrument rotates at a constant speed.

[0112] Figure 6 describes the information of 32 data points, namely the distribution curves of the rotation angles and distances of different data points from the detector. By converting the rotation angles and distances of the data points in Figure 6 into a rectangular coordinate system with the detection instrument as the coordinate origin, the plane rectangular coordinates of each data point can be obtained. The offset distance is calculated according to the above formulas 1 to 13, and the calculation 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 7 shows the relative position relationship of the casing shape before and after correction, where the origin coordinate point (0, 0) is the location of the instrument.

[0113] 2) The shape of the casing is a standard ellipse:

[0114] Assuming the semi-major axis of the ellipse is 1 and the semi-minor axis is 0.5, the position of the testing instrument is offset by 0.25 in the negative direction of the horizontal axis and 0.25 in the negative direction of the vertical axis relative to the center point. As shown in Figure 8, during the ultrasonic three-dimensional imaging process, the testing instrument rotates and scans along the circumference with itself as the center of rotation, capturing information for each test data point. This information includes the azimuth angle and distance of the data point located on the circumference of the ellipse. Similarly, based on the geometric relationship of the detector's rotational scanning in Figure 8, the information of each data point is solved, and the testing instrument rotates at a constant speed.

[0115] Figure 9 describes the information of 32 data points, namely the distribution curves of the rotation angles and distances of different data points from the detector. By converting the rotation angles and distances of the data points in Figure 9 into a rectangular coordinate system with the detection instrument as the coordinate origin, the plane rectangular coordinates of each data point can be obtained. According to the above formulas 1 to 13, the offset is calculated, and the calculation results are as follows: the horizontal axis offset is 0.245, and the vertical axis offset is 0.236. Compared with the actual horizontal axis offset of 0.25 and the actual vertical axis offset of 0.25, the correction error of the horizontal axis offset is 2.0%, and the correction error of the vertical axis offset is 5.5%. Figure 10 shows the relative position relationship of the casing morphology before and after correction, where the origin coordinate point (0, 0) is the location of the instrument.

[0116] 3) The shape of the casing is irregular:

[0117] 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 the shape with the horizontal axis is 1, and the intercept with the vertical axis is also 1. Compared with 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. As shown in Figure 11, during the three-dimensional imaging process of ultrasonic detection, the detection instrument uses itself as the center of rotation and scans along the circumference to capture information on each detection data point. This information includes the azimuth angle and distance of the data point on the circumference of the irregular shape. Similarly, based on the geometric relationship of the detection instrument's rotation scan in Figure 11, the information of each data point is solved, and the detection instrument rotates at a constant speed.

[0118] Figure 12 describes the information of 32 data points, namely the distribution curve of the rotation angle and distance of different data points from the detector. By converting the rotation angle and distance of the data points in Figure 12 into a rectangular coordinate system with the detection instrument as the coordinate origin, the plane rectangular coordinates of each data point can be obtained. According to the above formulas 1 to 13, the offset is calculated, and the calculation results are as follows: the horizontal axis offset is 0.509, and the vertical axis offset is 0.473. Compared with the actual horizontal axis offset of 0.5 and the actual vertical axis offset of 0.5, the correction error of the horizontal axis offset is 1.8%, and the correction error of the vertical axis offset is 5.4%. Figure 13 shows the relative position relationship of the casing shape before and after correction, where the origin coordinate (0, 0) is the position of the instrument.

[0119] The 3D imaging centering calibration will be analyzed based on a standard circle. Parameters required for the 3D calibration calculation include the number of data points per circumference, the number of turns of the helix, the offset when the instrument deviates, and the longitudinal position difference between adjacent data points. The casing is assumed to be a standard circle. In the simulated well section, the two ends each account for 30% of the total well length. The instrument is assumed to remain centered along the well axis throughout these two sections. The middle section, which accounts for 40% of the total well length, experiences instrument eccentricity. Data point information is simulated by setting eccentricity parameters. Assuming a standard circle radius of 1, a maximum eccentricity of 0.4 in the middle section, and 100 data points per circumference, a total of 150 circles (i.e., 150 turns of the helix) are calculated. The longitudinal coordinate difference between adjacent data points is 0.00024. After calibration, while the data points on the casing circumference in the eccentric section remain uneven, the resulting 3D shape is now centered and effectively reflects the true shape of the casing. It can be seen from the three-dimensional diagram and the top view that there is a certain error in the data verification results at the junction of eccentricity and centering. From the quantitative calculation results, the maximum verification calculation error is 3% when 100 data points are selected in a week, which can meet the engineering needs.

[0120] The above 1) to 3) respectively calibrate and calculate the well axis centering of the ultrasonic three-dimensional imaging of the casing deformation forms of standard circles, ellipses and irregular shapes. The number of data points selected in a week is 32. The calculation results are all good, and the errors are within the controllable range. The maximum calculation error is the longitudinal axis offset of the ellipse, which is 5.5%. The above three-dimensional imaging centering calibration is based on the standard circle analysis. From the quantitative calculation results, the maximum calibration calculation error is 3% when 100 data points are selected in a week, which can meet the engineering needs. According to the simulation calculation results of this method, when the number of data points continues to increase, the centering calibration calculation error will further decrease. Therefore, in actual engineering applications, an appropriate number of data points can be selected for detection and capture as needed to achieve the purpose of reducing the calculation error.

[0121] The method embodiments provided in the above embodiments of the present application can be executed in a mobile terminal, a computer terminal or a similar computing device. Taking operation on an electronic device as an example, FIG14 is a hardware structure block diagram of an electronic device for a method for correcting the well axis centering of ultrasonic three-dimensional imaging provided by the present application. As shown in FIG14 , the electronic device 10 may include one or more (only one is shown in the figure) processors 02 (the processor 02 may include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA), a memory 04 for storing data, and a transmission module 06 for communication functions. It can be understood by those skilled in the art that the structure shown in FIG14 is only for illustration, and it does not limit the structure of the above-mentioned electronic device. For example, the electronic device 10 may also include more or fewer components than those shown in FIG14 , or have a configuration different from that shown in FIG14 .

[0122] Memory 04 can be used to store software programs and modules for application software, such as the program instructions / modules corresponding to the method for correcting the well axis centering for ultrasonic three-dimensional imaging in the embodiments of the present application. Processor 02 executes the software programs and modules stored in memory 04 to perform various functional applications and data processing, thereby implementing the method for correcting the well axis centering for ultrasonic three-dimensional imaging in the aforementioned application. Memory 04 can include high-speed random access memory (RAM) and can also include non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some instances, memory 04 can further include memory remotely located relative to processor 02, which can be connected to electronic device 10 via a network. Examples of such networks include, but are not limited to, the Internet, corporate intranets, local area networks, mobile communication networks, and combinations thereof.

[0123] The transmission module 06 is used to receive or send data via a network. Specific examples of the aforementioned network may include a wireless network provided by the communication provider of the electronic device 10. In one embodiment, the transmission module 06 includes a network interface controller (NIC), which can be connected to other network devices via a base station to enable communication with the Internet. In another embodiment, the transmission module 06 may be a radio frequency (RF) module, which is used to communicate with the Internet wirelessly.

[0124] At the software level, the above device may be as shown in FIG15 and may include:

[0125] The scanning module 1501 is used to capture information of multiple target data points on the target casing by scanning along the circumferential direction using an ultrasonic detector;

[0126] A calculation module 1502 is configured to calculate the plane rectangular coordinates of each target data point among the plurality of target data points based on the information of the target data point;

[0127] An obtaining module 1503 is configured to obtain an arc length value of a position of each target data point among the multiple target data points based on the plane rectangular coordinates of each target data point among the multiple target data points;

[0128] A determination module 1504 is configured to determine an offset of the ultrasonic detector relative to the center of the well axis based on arc length values ​​of the locations of the plurality of target data points;

[0129] The correction module 1505 is configured to correct the plane rectangular coordinates of each target data point among the multiple target data points according to the offset.

[0130] In one embodiment, the information of the plurality of target data points may include: the azimuth angles and distances of the plurality of target data points; the calculation module 1502 may specifically convert the polar coordinates represented by the azimuth angles and distances of the plurality of target data points into plane rectangular coordinates represented by the horizontal coordinate and the vertical coordinate according to the following formula: i =r i ×cos(θ i ) y i =r i ×sin(θ i )

[0131] Among them, x i Indicates the horizontal coordinate of the i-th target data point, y i Represents the ordinate of the i-th target data point, r i represents the distance to the i-th target data point, θ iRepresents the azimuth angle of the i-th target data point.

[0132] In one embodiment, the obtaining module 1503 may specifically calculate the arc length value of each target data point in the plurality of target data points according to the following formula:

[0133] Among them, i is [1, n-1], s i Indicates the arc length value of the location of the i-th target data point, x i Indicates the horizontal coordinate of the i-th target data point, y i Indicates the vertical coordinate of the i-th target data point, x i+1 Indicates the horizontal coordinate of the i+1th target data point, y i+1 Indicates the ordinate of the i+1th target data point;

[0134] Among them, s n Indicates the arc length value of the location of the nth target data point, x1 indicates the horizontal coordinate of the first target data point, y1 indicates the vertical coordinate of the first target data point, x n Indicates the horizontal coordinate of the nth target data point, y n Indicates the vertical coordinate of the nth target data point.

[0135] In one embodiment, the above-mentioned determination module 1504 may include: a selection unit, used to select the minimum arc length value from the arc length values ​​of the positions of multiple target data points as a reference value; a first calculation unit, used to calculate the ratio between the arc length value of the position of each target data point among the multiple target data points and the reference value, and use the calculation result as the weighting factor of each target data point; a second calculation unit, used to calculate the geometric center coordinates of the circular curve according to the weighting factor of each target data point; an acquisition module, used to obtain the coordinates of the position of the ultrasonic detector as the reference origin coordinates; a determination unit, used to compare the geometric center coordinates with the reference origin coordinates to determine the offset of the position of the ultrasonic detector relative to the center of the well axis.

[0136] In one embodiment, calculating the geometric center coordinates of the circular curve according to the weighting factors of each target data point may include:

[0137] The geometric center coordinates of the circular curve are calculated according to the following formula:

[0138] Among them, X represents the horizontal coordinate of the geometric center coordinate, Y represents the vertical coordinate of the geometric center coordinate, and w i represents the weighting factor of the i-th target data point, w sumRepresents the cumulative value of the weighting factor of each target data point, x i Indicates the horizontal coordinate of the i-th target data point, y i Represents the vertical coordinate of the i-th target data point.

[0139] The embodiments of the present application also provide a specific implementation of an electronic device capable of implementing all steps of the method for correcting well axis centering for ultrasonic three-dimensional imaging in the above-mentioned embodiment. The electronic device specifically includes the following: a processor, a memory, a communication interface, and a bus; wherein the processor, the memory, and the communication interface communicate with each other via the bus; the processor is configured to call a computer program in the memory, and when the processor executes the computer program, all steps of the method for correcting well axis centering for ultrasonic three-dimensional imaging in the above-mentioned embodiment are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0140] Step 1: Use an ultrasonic detector to scan along the circumference to capture information of multiple target data points on the target casing;

[0141] Step 2: Calculate the plane rectangular coordinates of each target data point among the multiple target data points based on the information of the target data point;

[0142] Step 3: Calculate the arc length value of each target data point among the multiple target data points based on the plane rectangular coordinates of each target data point;

[0143] Step 4: Determine the offset of the ultrasonic detector relative to the center of the well axis based on the arc length values ​​of the locations of the multiple target data points;

[0144] Step 5: Correct the plane rectangular coordinates of each target data point among the multiple target data points according to the offset.

[0145] Embodiments of the present application also provide a computer-readable storage medium capable of implementing all steps of the method for correcting well axis centering for ultrasonic three-dimensional imaging in the above-mentioned embodiment. The computer-readable storage medium stores a computer program. When the computer program is executed by a processor, the computer program implements all steps of the method for correcting well axis centering for ultrasonic three-dimensional imaging in the above-mentioned embodiment. For example, when the processor executes the computer program, the following steps are implemented:

[0146] Step 1: Use an ultrasonic detector to scan along the circumference to capture information of multiple target data points on the target casing;

[0147] Step 2: Calculate the plane rectangular coordinates of each target data point among the multiple target data points based on the information of the target data point;

[0148] Step 3: Calculate the arc length value of each target data point among the multiple target data points based on the plane rectangular coordinates of each target data point;

[0149] Step 4: Determine the offset of the ultrasonic detector relative to the center of the well axis based on the arc length values ​​of the locations of the multiple target data points;

[0150] Step 5: Correct the plane rectangular coordinates of each target data point among the multiple target data points according to the offset.

[0151] As can be seen from the above description, the embodiment of the present application captures information of multiple target data points on the target casing by scanning along the circumferential direction with an ultrasonic detector; calculates the plane rectangular coordinates of each target data point in the multiple target data points based on the information of the target data points; obtains the arc length value of each target data point in the multiple target data points based on the plane rectangular coordinates of each target data point in the multiple target data points; determines the offset of the ultrasonic detector position relative to the well axis centering based on the arc length value of the multiple target data points; and corrects the plane rectangular coordinates of each target data point in the multiple target data points based on the offset. That is, the arc length difference between the data points is used to determine the eccentricity of the detection instrument, and the calculated eccentricity can be used to check the well axis centering of the casing deformation ultrasonic three-dimensional imaging, thereby solving the problem of low three-dimensional imaging accuracy caused by the inability to check the well axis centering, and achieving the technical effect of effectively improving the three-dimensional imaging accuracy.

[0152] 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, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with the relevant laws, regulations and standards of relevant countries and regions, and provide corresponding operation entrances for users to choose to authorize or refuse.

[0153] The various embodiments in this specification are described in a progressive manner. Similar portions between the various embodiments can be referenced to each other. Each embodiment focuses on the differences between the other embodiments. In particular, the hardware + program embodiments are generally similar to the method embodiments, so their description is relatively simple. For relevant portions, refer to the description of the method embodiments.

[0154] The foregoing description of this specification describes specific embodiments. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims can be performed in an order different from that described in the embodiments and still achieve the desired results. Furthermore, the processes depicted in the accompanying drawings do not necessarily require the specific order shown or the sequential order to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0155] Although the present application provides method operation steps as described in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-creative work. The order of steps listed in the embodiments is only one way of executing the steps among many steps and does not represent the only execution order. When the actual device or client product is executed, it can be executed in sequence or in parallel according to the method shown in the embodiments or the drawings (for example, in a parallel processor or multi-threaded processing environment).

[0156] The systems, devices, modules, or units described in the above embodiments may be implemented by computer chips or entities, or by products having certain functions. A typical implementation device is a computer. Specifically, the computer may be, for example, a personal computer, a laptop computer, an in-vehicle human-computer 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 a combination of any of these devices.

[0157] Although the present specification embodiment provides the method operation steps as described in the embodiment or flow chart, more or less operation steps may be included based on conventional or non-creative means. The order of steps listed in the embodiment is only one way in the order of execution of many steps and does not represent a unique execution order. When the device or terminal product in practice is executed, it can be performed in sequence or in parallel according to the method shown in the embodiment or the accompanying drawings (such as a parallel processor or a multi-threaded processing environment, or even a distributed data processing environment). The term "comprise", "comprising" or any other variant thereof is intended to cover non-exclusive inclusion, so that the process, method, product or equipment including a series of elements not only include those elements, but also include other elements not clearly listed, or also include elements inherent to such process, method, product or equipment. In the absence of more restrictions, it is not excluded that there are other identical or equivalent elements in the process, method, product or equipment including the elements.

[0158] For the convenience of description, the above devices are described in terms of functions divided into various modules. Of course, when implementing the embodiments of this specification, the functions of each module can be implemented in the same or multiple software and / or hardware, or the modules that implement the same function can be implemented by a combination of multiple sub-modules or sub-units, etc. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, which can be electrical, mechanical or other forms.

[0159] Those skilled in the art will also appreciate that, in addition to implementing the controller in pure computer-readable program code, it is entirely possible to implement the same functionality by logically programming the method steps in the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, embedded microcontrollers, and the like. Therefore, such a controller can be considered a hardware component, and the devices included therein for implementing various functions can also be considered structures within the hardware component. Alternatively, the devices for implementing various functions can be considered both software modules implementing the method and structures within the hardware component.

[0160] The present application is described with reference to the flow chart and / or block diagram of the method, device (system), and computer program product according to the embodiment of the present application. It should be understood that each flow process and / or box in the flow chart and / or block diagram and the combination of the flow process and / or box in the flow chart and / or block diagram can be realized by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processing machine or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for realizing the function specified in one flow chart flow or multiple flows and / or one box or multiple boxes of the block diagram.

[0161] These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a product including an instruction device that implements the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

[0162] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

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

[0164] Memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. Memory is an example of a computer-readable medium.

[0165] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. The information can be computer-readable instructions, data structures, program modules, 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 technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic disk storage or other magnetic storage devices, or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transitory media such as modulated data signals and carrier waves.

[0166] Those skilled in the art will appreciate that the embodiments of this specification may be provided as methods, systems, or computer program products. Therefore, the embodiments of this specification may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the embodiments of this specification may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0167] Embodiments of this specification may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, and the like that perform specific tasks or implement specific abstract data types. Embodiments of this specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected via a communications network. In distributed computing environments, program modules may be located in local and remote computer storage media, including storage devices.

[0168] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between the various embodiments can be referenced across them. Each embodiment focuses on the differences from the other embodiments. In particular, since the system embodiments are generally similar to the method embodiments, their description is relatively simple. For relevant parts, reference can be made to the description of the method embodiments. Throughout this specification, reference to the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with that embodiment or example are included in at least one embodiment or example of the embodiments in this specification. In this specification, the schematic representations of these terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, those skilled in the art may combine and integrate the different embodiments or examples, and features of different embodiments or examples, described in this specification, without conflict.

[0169] The above description is merely an example of the embodiments of this specification and is not intended to limit the embodiments of this specification. For those skilled in the art, various modifications and variations of the embodiments of this specification are possible. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of this specification shall be included within the scope of the claims of the embodiments of this specification.

Claims

1. A method for correcting the well axis centering of ultrasonic three-dimensional imaging, characterized in that: The method comprises: The ultrasonic detector is used to scan along the circumferential direction to capture information of multiple target data points on the target casing; Calculating the plane rectangular coordinates of each target data point among the multiple target data points through the information of the target data point; Calculating the arc length value of each target data point among the multiple target data points according to the plane rectangular coordinates of each target data point among the multiple target data points; Determining the offset of the ultrasonic detector relative to the center of the well axis according to the arc length values ​​of the locations of the multiple target data points; According to the offset, the plane rectangular coordinates of each target data point among the multiple target data points are corrected.

2. The method according to claim 1, characterized in that The information of the plurality of target data points includes: the azimuth angles and distances of the plurality of target data points; Correspondingly, the plane rectangular coordinates of each target data point in the plurality of target data points are calculated based on the information of the target data point, including: The polar coordinates represented by the azimuth angles and distances of the plurality of target data points are converted into plane rectangular coordinates represented by abscissas and ordinates according to the following formula; x i =r i ×cos(θ i ) y i =r i ×sin(θ i ) Among them, x i Represents the horizontal coordinate of the i-th target data point, y i represents the ordinate of the i-th target data point, r i represents the distance of the i-th target data point, θ i Represents the azimuth angle of the i-th target data point.

3. The method according to claim 1, characterized in that Obtaining the arc length value of each target data point among the multiple target data points according to the plane rectangular coordinates of each target data point, including: The arc length value of each target data point among multiple target data points is calculated according to the following formula: Among them, i is [1, n-1], s i Indicates the arc length value of the location of the i-th target data point, x i Represents the horizontal coordinate of the i-th target data point, y i Represents the ordinate of the i-th target data point, x i+1 Indicates the horizontal coordinate of the i+1th target data point, y i+1 Indicates the ordinate of the i+1th target data point; Among them, s n represents the arc length value of the location of the nth target data point, x1 represents the horizontal coordinate of the first target data point, y1 represents the vertical coordinate of the first target data point, and x n Indicates the horizontal coordinate of the nth target data point, y n Indicates the ordinate of the nth target data point.

4. The method according to claim 1, characterized in that: Determining the offset of the ultrasonic detector relative to the center of the well axis according to the arc length values ​​of the locations of the multiple target data points includes: Selecting a minimum arc length value from the arc length values ​​of the locations of the plurality of target data points as a reference value; Calculating the ratio between the arc length value of each target data point in the plurality of target data points and the reference value, and using the calculation result as a weighting factor of each target data point; Calculate the geometric center coordinates of the circular curve according to the weighting factors of each target data point; Obtaining the coordinates of the position of the ultrasonic detector as the reference origin coordinates; The geometric center coordinates are compared with the reference origin coordinates to determine the offset of the ultrasonic detector relative to the center of the well axis.

5. The method according to claim 4, characterized in that According to the weighting factors of each target data point, the geometric center coordinates of the circular curve are calculated, including: The geometric center coordinates of the circular curve are calculated according to the following formula: Among them, X represents the horizontal coordinate of the geometric center coordinate, Y represents the vertical coordinate of the geometric center coordinate, and w i represents the weighting factor of the i-th target data point, w sum Represents the cumulative value of the weighted factors of each target data point, x i Represents the horizontal coordinate of the i-th target data point, y i Represents the ordinate of the i-th target data point.

6. The method according to any one of claims 1 to 5, characterized in that The ultrasonic detector is used in detecting casing deformation.

7. A device for correcting the well axis centering of ultrasonic three-dimensional imaging, characterized in that: include: A scanning module, used for capturing information of multiple target data points on the target casing by scanning along the circumferential direction through an ultrasonic detector; A calculation module, used for calculating the plane rectangular coordinates of each target data point among the multiple target data points through information of the target data point; A obtaining module is used to obtain the plane rectangular coordinates of each target data point in the plurality of target data points. The arc length value of each target data point among multiple target data points; A determination module, used to determine the offset of the position of the ultrasonic detector relative to the center of the well axis according to the arc length values ​​of the positions of the multiple target data points; A correction module is used to correct the plane rectangular coordinates of each target data point among the multiple target data points according to the offset.

8. The device according to claim 7, characterized in that The determination module comprises: A selection unit, configured to select a minimum arc length value from the arc length values ​​of the locations of the plurality of target data points as a reference value; A first calculation unit, used for calculating the ratio between the arc length value of each target data point among the plurality of target data points and the reference value, and using the calculation result as a weighting factor of each target data point; A second calculation unit, used to calculate the geometric center coordinates of the circular curve according to the weighting factors of each target data point; An acquisition module, used to acquire the coordinates of the position of the ultrasonic detector as the reference origin coordinates; The determination unit is used to compare the geometric center coordinates with the reference origin coordinates to determine the offset of the position of the ultrasonic detector relative to the center of the well axis.

9. An electronic device comprising a processor and a memory for storing instructions executable by the processor, characterized in that: When the processor executes the instructions, the steps of the method according to any one of claims 1 to 6 are implemented.

10. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instructions are executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images