Method and system for determining the through capacity of a casing deformation while drilling tool

Abstract

The application provides a method and system for determining the passing capacity of a casing deformation post-entry tool, and the method comprises the following steps: obtaining the center of a casing section according to a well diameter two-dimensional space coordinate; calculating the three-dimensional space coordinates of the center of the casing section and each arm end point according to the obtained well trajectory, depth and well diameter data, and establishing a three-dimensional wellbore model after casing deformation; and calculating the maximum passing diameter of the tool in the casing deformation well section by using a projection method. The calculation system can realize the method for determining the passing capacity of the casing deformation post-entry tool. The method for determining the passing capacity of the casing deformation post-entry tool provided by the application overcomes the distortion drawbacks caused by the previous geometric algorithm based on the wellbore curvature change, and improves the calculation accuracy of the tool passing capacity.

Description

Technical Field

[0001] This invention relates to the fields of geophysical logging and oil and gas development engineering, specifically to a method for determining the passability of a well-entry tool after casing deformation, a system for determining the passability of a well-entry tool after casing deformation, a computer device, and a computer-readable storage medium storing a computer program. Background Technology

[0002] Shale oil and gas has gradually become the main battleground supporting the vigorous development of my country's energy industry, with its production share increasing year by year. During shale oil and gas exploration and development, casing well deformation events occur frequently, seriously affecting the efficiency of oil and gas exploration and development and the timeliness of construction operations. After casing deformation, the calculation of the well tool string's throughput capacity and the optimization of the tubing string are crucial to the success or failure of subsequent construction operations. However, the complex geometric spatial morphology and difficulty in characterizing deformed casing severely restrict the calculation of the well tool string's throughput capacity. How to accurately characterize the three-dimensional wellbore spatial characteristics after casing deformation and accurately calculate the well tool string's throughput capacity has become a technical challenge that urgently needs to be overcome.

[0003] Current research on downhole tool passability mainly employs geometric and mechanical analysis methods. Most studies use drill strings and drill pipes as examples, simplifying the tool string into a uniform cross-section beam model subjected to the lateral component of net weight and axial pressure. Furthermore, most studies assume the tool string gets stuck at the maximum dogleg. Ran Jing et al. (1990) proposed a method for determining downhole tool passability under rigid conditions. Zhao Junping et al. (1993) improved the rigid passability model and established a mechanical model for solving downhole tool passability under flexible conditions using the longitudinal and transverse bending method. Di Qinfeng et al. (1996) considered the influence of centralizers on drill string passability, but neglected the influence of axial stress. He Shiming et al. (1997) established a friction calculation model for casing in horizontal wells through mechanical analysis. Later, Sui Mancang et al. (1999) considered the influence of rotational running in casing friction analysis. Wang Yanhong (2008) discussed the mechanical analysis method for downhole tool passability. Zhu Xiuxing et al. (2013) established the governing equation for the length of the joint perforation string under rigid conditions, but did not consider the deformation of the string. Feng Ding et al. (2016) conducted a study on the passability of multi-layer perforation strings under the assumption that the location of the maximum curvature of the wellbore trajectory is the jamming point of the sub-injection string. Liu Jun et al. (2021) established a wellbore passability analysis model for cable-pumped cluster perforation strings based on a comprehensive consideration of factors such as friction between downhole tools and the wellbore, wellbore geometric constraints, pump thrust, axial tension, variable cross-section of the string, cable head tension, and tool elastic deformation. They then used geometric analysis and longitudinal and transverse bending methods to establish and solve the complex coefficient equations of the model.

[0004] Chinese patents CN111241684B and CN111101931B disclose two methods for determining the throughput capacity of clustered perforated tubing strings in wells with cable pumps. Both methods assume the tubing string throughput capacity is determined under conditions of undeformed cylindrical casing. These methods yield good calculation results for high-rigidity tool strings or undeformed cylindrical casing. However, after casing deformation, the wellbore is no longer a regular cylinder, and the wellbore axis changes. Using existing geometric or mechanical analysis methods for determination will lead to erroneous calculation results and may even cause more complex accidents. Therefore, a new method for determining the throughput capacity of tool strings in wells with deformed casing is needed to meet the needs of field production. Summary of the Invention

[0005] In view of the shortcomings of the prior art, the purpose of this invention is to solve one or more of the problems existing in the prior art. For example, one objective of this invention is to provide a method and system for determining the passability of a well-entry tool after casing deformation.

[0006] To achieve the above objectives, the present invention provides a method for determining the passage capability of a well tool after casing deformation. The method may include the following steps: determining the center of the casing cross-section based on the two-dimensional spatial coordinates of the wellbore diameter; calculating the three-dimensional spatial coordinates of the center of the casing cross-section and the endpoints of each arm based on the acquired wellbore trajectory, depth, and wellbore diameter data, and establishing a three-dimensional wellbore model after casing deformation; and calculating the maximum passage diameter of the tool in the casing deformation section using the projection method.

[0007] According to one or more exemplary embodiments of one aspect of the present invention, the method for determining the center of the casing cross-section may include at least one of the centroid method, the least squares method, and the iterative search method.

[0008] According to one or more exemplary embodiments of one aspect of the present invention, the casing cross section may include at least one of a casing interface of a specific depth and a casing deformation cross section.

[0009] According to one or more exemplary embodiments of the present invention, the method for calculating the three-dimensional spatial coordinates of the center of the casing cross-section and the endpoints of each arm may include a translational rotation replacement method, wherein the translational rotation replacement method is used to complete well trajectory reconstruction and wellbore imaging, and the well trajectory reconstruction process may include:

[0010] Set parameters, set the first The three-dimensional coordinates of the center of the casing section at each depth point Translation vector That is, the first Original trajectory point coordinates at each depth point; rotation vector ,in, For the first Tangent of the trajectory point at each depth point for The normal to the plane, , and x The axis is parallel; the rotation angle For the first There are several depth points, with depth valued as Depth, Depth= h 0+ d k , k =0, 1, ..., m-1, d The depth sampling interval for the data conversion and merging unit;

[0011] Translate the coordinates to move the point According to the translation vector Translate to point The calculation formula includes Equation 1: ;

[0012] Rotate the coordinates to make the point Along the axis of rotation Rotation angle After that, I got points. The calculation formula includes Equation 2: ;

[0013] Replace the coordinates with the central 3D spatial coordinates obtained after translation and rotation. Replace the Original trajectory coordinates of each depth point ;

[0014] The wellbore imaging process includes: using the above-mentioned translation and rotation method to image the first... Coordinates of a depth point Perform rotation, initial number The three-dimensional coordinate vectors of the endpoints are: Its rotation axis along a rotation axis Rotation angle Then, the rotated vector is obtained. Its calculation formula includes Equation 3: ,Will Substitute them into the equation to obtain the rotated point set. Next, a displacement operation is performed, resulting in the set of points after displacement. for: .

[0015] According to one or more exemplary embodiments of one aspect of the present invention, establishing a three-dimensional wellbore model after casing deformation may include reconstructing the well trajectory and performing three-dimensional wellbore modeling.

[0016] According to one or more exemplary embodiments of one aspect of the present invention, the method for establishing a three-dimensional wellbore model after casing deformation may include data sample interpolation using Lagrange interpolation.

[0017] According to one or more exemplary embodiments of the present invention, the calculation of the maximum through diameter of the tool in the casing deformation section using the projection method may include the following steps: Step 1: Set the projection plane, and project the first The normal plane of the well trajectory at each depth point is set as the target projection plane, which is also of length [missing information]. L The bottom cut plane of the three-dimensional deformable wellbore, from the first From the depth point to the The coordinates of all caliper endpoints at a given depth point are set as the points to be projected, where... Step 2: Set the depth sampling interval for data conversion and merging units; Step 3: Set the projection direction, mesh the target projection plane, and obtain the set of mesh nodes. Then the first Well trajectory coordinates at each depth point Pointing to a certain grid node Projection direction The expression includes Equation 4: Step 3: Calculate the projected coordinates and set them. A point on the well wall to be projected Coordinates are The normal to the target projection plane is , The projection point on the target projection plane is , where represents , The system of equations is as follows: Step 4: Obtain the inner boundary of the projection area, and set it from the first... From the depth point to the The wellbore at the depth point to be projected is at the... The projection point of the normal plane containing the well trajectory at depth n and the nth depth point The set consisting of wellhead endpoints at depth points is Get point set Step 5: Calculate the largest inscribed circle of the inner boundary point set D of the enclosed region, denoted as . Step 6: Determine the length as L The tool with the longest diameter can pass through; Step 7: Repeat steps 2 to 6 to calculate the maximum inscribed circle diameter in each projection direction, and determine the maximum value of all diameters as the length. L The maximum diameter of the meter tool .

[0018] According to one or more exemplary embodiments of one aspect of the present invention, the method for calculating the maximum inscribed circle of the inner boundary set D may include an iterative search method.

[0019] Another aspect of the present invention provides a system for determining the passability of a wellbore tool after casing deformation. The system can implement the method described above for determining the passability of a wellbore tool after casing deformation. The system may include: a data conversion and merging unit configured to integrate the original two-dimensional coordinate data of the well trajectory and the two-dimensional spatial coordinates of the well diameter after casing deformation; a well diameter data correction unit configured to determine the center of the casing cross-section and compare the obtained centers; a three-dimensional wellbore modeling unit configured to calculate the three-dimensional spatial coordinates of the center and each arm endpoint, reconstruct the well trajectory, and establish a three-dimensional wellbore model after casing deformation; and a passability calculation unit configured to calculate the maximum passability diameter of the tool in the casing deformation section.

[0020] According to one or more exemplary embodiments of another aspect of the present invention, the data conversion and merging unit may include at least one of well trajectory data, well caliper logging data, casing inner diameter and depth sampling interval input modules; the well caliper data correction unit may include at least one of raw data calculation and storage, geometric center method data calculation and storage, least squares method data calculation and storage, iterative search method data calculation and storage, and calculation center two-dimensional drawing module; the three-dimensional wellbore modeling unit may include at least one of well section setting, well trajectory three-dimensional space display, and well wall three-dimensional space modeling module; the throughput capacity calculation unit may include at least one of tool length input, calculation step length setting, and projection method calculation module.

[0021] Another aspect of the present invention provides a computer device, which may include: a processor; and a memory storing a computer program, which, when executed by the processor, enables the method described above for determining the passability of the well tool after casing deformation.

[0022] Another aspect of the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, enables the method described above for determining the passability of a well tool after casing deformation.

[0023] Compared with the prior art, the beneficial effects of the present invention include at least one of the following:

[0024] (1) The method for determining the well passability of the casing deformation tool provided by the present invention is based on the three-dimensional features of the wellbore of the casing deformation well. It uses the translation and rotation replacement method to reconstruct the well trajectory and truly restore the three-dimensional wellbore spatial features of the casing deformation well.

[0025] (2) The method for determining the passage capacity of the casing after deformation provided by the present invention extracts the well trajectory coordinates and well wall endpoint coordinates based on the spatial real shape characteristics of the three-dimensional wellbore after deformation, and uses the interval projection method to determine the maximum passage diameter of the cylindrical tool with a length of L meters. This method overcomes the distortion drawbacks of the previous geometric algorithm based on the curvature change of the wellbore and improves the accuracy of the tool passage capacity calculation. Attached Figure Description

[0026] The above and other objects and / or features of the present invention will become clearer from the following description taken in conjunction with the accompanying drawings, in which:

[0027] Figure 1 A flowchart illustrating a method for determining the wellbore penetration capability after casing deformation in an exemplary embodiment of the present invention is shown.

[0028] Figure 2 A schematic diagram illustrating the use of a translation-rotation substitution method to translate coordinates in an exemplary embodiment of the present invention is shown;

[0029] Figure 3 A schematic diagram illustrating the use of a translation-rotation substitution method to rotate coordinates in an exemplary embodiment of the present invention is shown;

[0030] Figure 4 A schematic diagram illustrating the use of a translation and rotation substitution method to replace coordinates in an exemplary embodiment of the present invention is shown;

[0031] Figure 5 A schematic diagram illustrating the use of a projection method to set the projection plane in an exemplary embodiment of the present invention is shown;

[0032] Figure 6 A schematic diagram illustrating the use of a projection method to set the projection direction in an exemplary embodiment of the present invention is shown;

[0033] Figure 7 A schematic diagram illustrating the calculation of projected coordinates using a projection method in an exemplary embodiment of the present invention is shown;

[0034] Figure 8 A schematic diagram illustrating the use of a projection method to obtain the inner boundary of a projection area in an exemplary embodiment of the present invention is shown.

[0035] Figure 9 A schematic diagram illustrating the calculation of the maximum inscribed circle of the inner boundary using the projection method in an exemplary embodiment of the present invention is shown.

[0036] Figure 10 A schematic diagram of the composition of a system for calculating the wellbore penetration capability after casing deformation is shown in an exemplary embodiment of the present invention.

[0037] Explanation of reference numerals in the attached figures:

[0038] 100 - Data conversion and merging unit; 200 - Wellbore data correction unit; 300 - 3D wellbore modeling unit; 400 - Through capacity calculation unit. Detailed Implementation

[0039] In the following sections, a method and system for determining the passability of a casing-deformed well tool will be described in detail with reference to the accompanying drawings and exemplary embodiments.

[0040] The method for calculating the passage capacity of the well tool after casing deformation of the present invention utilizes wellbore trajectory data such as measured well inclination angle and well inclination azimuth, as well as well diameter data reflecting the well wall or casing wall measured after casing deformation. The method uses translation and rotation replacement to establish the true three-dimensional wellbore shape after casing deformation, and uses projection method to calculate the passage capacity of a cylindrical tool of preset specifications.

[0041] Exemplary Example 1

[0042] Figure 1 A flowchart illustrating a method for determining the wellbore penetration capability after casing deformation in an exemplary embodiment of the present invention is shown. Figure 2 A schematic diagram illustrating the use of a translation-rotation substitution method to translate coordinates in an exemplary embodiment of the present invention is shown; Figure 3 A schematic diagram illustrating the use of a translation-rotation substitution method to rotate coordinates in an exemplary embodiment of the present invention is shown; Figure 4 A schematic diagram illustrating the use of a translation and rotation substitution method to replace coordinates in an exemplary embodiment of the present invention is shown; Figure 5 A schematic diagram illustrating the use of a projection method to set the projection plane in an exemplary embodiment of the present invention is shown; Figure 6 A schematic diagram illustrating the use of a projection method to set the projection direction in an exemplary embodiment of the present invention is shown; Figure 7 A schematic diagram illustrating the calculation of projected coordinates using a projection method in an exemplary embodiment of the present invention is shown; Figure 8 A schematic diagram illustrating the use of a projection method to obtain the inner boundary of a projection area in an exemplary embodiment of the present invention is shown. Figure 9 This diagram illustrates a method for calculating the maximum inscribed circle of the inner boundary using a projection method in an exemplary embodiment of the present invention.

[0043] This exemplary embodiment provides a method for determining the passability of a well-entry tool after casing deformation.

[0044] like Figure 1 As shown, the method for determining the wellbore penetration capability after casing deformation mainly includes the following steps:

[0045] S1. Obtain the two-dimensional spatial coordinates of the wellbore after casing deformation, and determine the center of the casing deformation section at a specific depth. Here, by obtaining the wellbore diameter data reflecting the well wall or casing wall after casing deformation, the two-dimensional spatial coordinates of the wellbore diameter are obtained, and the center of the casing section (including the casing deformation section) at a specific depth is calculated.

[0046] S2. Calculate the three-dimensional spatial coordinates of the center and each arm endpoint of the casing cross-section at a specific depth, reconstruct the well trajectory, and establish a three-dimensional wellbore model after casing deformation. Here, by acquiring wellbore trajectory data such as well inclination angle and well inclination azimuth, as well as depth data, the three-dimensional wellbore trajectory spatial coordinates are obtained, and the three-dimensional spatial coordinates of the center and each arm endpoint of the casing cross-section (including the casing deformation cross-section) at a specific depth are calculated to reconstruct the well trajectory, perform wellbore imaging, and establish a three-dimensional wellbore model after casing deformation.

[0047] S3. Calculate the maximum through diameter of the rigid cylindrical tool in the casing section using the projection method.

[0048] In this exemplary embodiment, the method for determining the center of the casing cross-section in step S1 may include at least one of the centroid method, the least squares method, and the iterative search method.

[0049] In this exemplary embodiment, in step S2, a translational rotation replacement method can be used to complete well trajectory reconstruction and wellbore imaging. The well trajectory reconstruction process may include:

[0050] Set parameters, set the first The three-dimensional coordinates of the center of the casing section at each depth point Translation vector That is, the first Original trajectory point coordinates at each depth point; rotation vector ,in, For the first Tangent of the trajectory point at each depth point for The normal to the plane, Its and x The axis is parallel; the rotation angle For the first There are several depth points, with depth valued as Depth, Depth= h 0+ d k , k =0, 1, ..., m-1, d The depth sampling interval for the data conversion and merging unit;

[0051] Translation coordinates, such as Figure 2 As shown, the point According to the translation vector Translate to point The calculation formula includes Equation 1: ;

[0052] Rotate coordinates, such as Figures 3 to 4 As shown, the point Along the axis of rotation Rotation angle After that, I got points. The calculation formula includes Equation 2: ;

[0053] Replace coordinates, such as Figure 4 As shown, the central three-dimensional spatial coordinates are obtained after translation and rotation. Replace the Original trajectory coordinates of each depth point ;

[0054] The wellbore imaging process may include: using the above-described translation and rotation method to image the first... Coordinates of a depth point Perform rotation, initial number The three-dimensional coordinate vectors of the endpoints are: Its rotation axis along a rotation axis Rotation angle Then, the rotated vector is obtained. Its calculation formula includes Equation 3: ,Will Substitute them into the equation to obtain the rotated point set. Next, a displacement operation is performed, resulting in the set of points after displacement. for: .

[0055] In this exemplary embodiment, establishing a three-dimensional wellbore model after casing deformation may include reconstructing the well trajectory and performing three-dimensional wellbore modeling. The method for establishing a three-dimensional wellbore model after casing deformation may include interpolating data samples using Lagrange interpolation to establish a more accurate three-dimensional wellbore model.

[0056] In this exemplary embodiment, calculating the maximum passing diameter of the tool in the deformed section of the casing using the projection method may include the following steps:

[0057] Step 1: Set the projection plane, such as Figure 5 As shown, the first The normal plane of the well trajectory at each depth point is set as the target projection plane, which is also of length [missing information]. L The bottom cut plane of the three-dimensional deformable wellbore, from the first From the depth point to the The coordinates of all caliper endpoints at a given depth point are set as the points to be projected, where... This refers to the depth sampling interval of the data conversion and merging unit.

[0058] Step 2: Set the projection direction, such as... Figure 6 As shown, the target projection plane is meshed to obtain a set of mesh nodes. Then the first Well trajectory coordinates at each depth point Pointing to a certain grid node Projection direction The expression includes Equation 4: .

[0059] Step 3: Calculate the projected coordinates, such as Figure 7 As shown, setting A point on the well wall to be projected Coordinates are The normal to the target projection plane is , The projection point on the target projection plane is , where represents , The system of equations is as follows: .

[0060] Step 4: Obtain the inner boundary of the projection area, such as... Figure 8 As shown, the setting starts from the first From the depth point to the The wellbore at the depth point to be projected is at the... The projection point of the normal plane containing the well trajectory at depth n and the nth depth point The set consisting of wellhead endpoints at depth points is Get point set The set of inner boundary points D of the enclosed region.

[0061] Step 5: Calculate the largest incircle of the inner boundary set D, such as... Figure 9 As shown, the diameter is denoted as Here, methods for calculating the maximum inscribed circle of the inner boundary set D may include iterative search methods.

[0062] Step 6: Determine the length as L The tool with the largest diameter that can pass through is the meter.

[0063] Step 7: Traverse all feasible projection directions, repeating steps 2 to 6, calculate the maximum inscribed circle diameter for each projection direction, and determine the maximum value of all diameters as the length. L The maximum diameter of the cylindrical tool in meters .

[0064] Exemplary Example 2

[0065] Figure 10A schematic diagram of the composition of a system for calculating the wellbore penetration capability after casing deformation is shown in an exemplary embodiment of the present invention.

[0066] This exemplary embodiment provides a calculation system for the passability of a well-entry tool after casing deformation. This calculation system can implement the method for determining the passability of a well-entry tool after casing deformation as described in Exemplary Embodiment 1.

[0067] like Figure 10 As shown, the system for calculating the throughput capacity of the casing-deformed well tool may include: a data conversion and merging unit 100, configured to integrate the original well trajectory two-dimensional coordinate data and the two-dimensional spatial coordinates of the well diameter after casing deformation; a well diameter data correction unit 200, configured to determine the center of the casing section and compare the obtained center; a three-dimensional wellbore modeling unit 300, configured to calculate the three-dimensional spatial coordinates of the center and each arm endpoint, reconstruct the well trajectory, and establish a three-dimensional wellbore model after casing deformation; and a throughput capacity calculation unit 400, configured to calculate the maximum throughput diameter of the tool in the casing-deformed section.

[0068] In this exemplary embodiment, the data conversion and merging unit may include at least one of the following: well trajectory data, well caliper logging data, casing inner diameter and depth sampling interval input modules; the well caliper data correction unit may include at least one of the following: raw data calculation and storage, geometric center method data calculation and storage, least squares method data calculation and storage, iterative search method data calculation and storage, and calculation center two-dimensional drawing module; the three-dimensional wellbore modeling unit may include at least one of the following: well section setting, well trajectory three-dimensional space display, and well wall three-dimensional space modeling module; and the capability calculation unit may include at least one of the following: tool length input, calculation step length setting, and projection method calculation module.

[0069] In this exemplary embodiment, the method / process of using the system for calculating the throughput capacity of the casing-deformed tool may include: calling the data conversion and merging unit 100 to associate the original well trajectory coordinate data (including parameters such as depth, inclination angle, and azimuth angle) and the casing-deformed wellbore logging data (including parameters such as depth, wellbore length vector, and wellbore endpoint azimuth vector) with depth as the index; using the casing-deformed wellbore logging data as input, calling the wellbore data correction unit 200 to obtain the center of the two-dimensional space of the casing cross-section and outputting the obtained center and endpoint coordinates; using the original well trajectory coordinate data and the center and endpoint coordinates of the two-dimensional space of the casing cross-section as input, calling the three-dimensional wellbore modeling unit 300 to calculate the three-dimensional spatial coordinates of the center and each arm endpoint, reconstructing the well trajectory, imaging the wellbore wall, and establishing a three-dimensional wellbore model after casing deformation; and calling the throughput capacity calculation unit 400 to calculate the maximum throughput diameter of the tool in the casing-deformed section using the three-dimensional wellbore model after casing deformation.

[0070] Exemplary Example 3

[0071] This exemplary embodiment provides a computer device. The computer device includes a processor and a memory. The memory stores a computer program. The computer program is executed by the processor, causing the processor to perform a method for determining the passability of a casing-deformed well tool according to the present invention.

[0072] Exemplary Example 4

[0073] This exemplary embodiment provides a computer-readable storage medium storing a computer program. The computer-readable storage medium stores a computer program that, when executed by a processor, causes the processor to perform a method for determining the passability of a casing-deformed well tool according to the present invention. The computer-readable recording medium is any data storage device capable of storing data readable by a computer system. Examples of computer-readable recording media include: read-only memory, random access memory, read-only optical disc, magnetic tape, floppy disk, optical data storage device, and carrier waves (such as data transmission via the Internet through wired or wireless transmission paths).

[0074] To better understand the exemplary embodiments described above, further explanation will be provided below with reference to specific examples.

[0075] Example 1

[0076] Calculations show that the casing deformation in well XX is most severe at a depth of 3331.75m. The casing deformation parameters are shown in Table 1.

[0077] Table 1. Statistical results of casing deformation parameters in the 3300-3350m section of Well XX.

[0078]

[0079] Where, the maximum deformation = max{the first k The depth point i The maximum deformation rate is calculated as follows: (Inner diameter of each endpoint - standard inner diameter) / (Maximum deformation amount / standard inner diameter).

[0080] In this example, the average maximum passable outer diameter of the 5m tool string is 93.06mm, the average maximum passable outer diameter of the 10m tool string is 93.0mm, the average maximum passable outer diameter of the 15m tool string is 92.45mm, the average maximum passable outer diameter of the 20m tool string is 85.74mm, and the average maximum passable outer diameter of the 25m tool string is 77.72mm. The specific results are shown in Table 2.

[0081] Table 2. Results of tool strings of different lengths at depths of 3300-3350m in Well XX, based on diameter.

[0082]

[0083] As shown in Table 2, the results match the actual situation by 90%, which is 50% more accurate than the conventional geometric algorithm based on wellbore curvature variation. This provides a reliable reference for handling complex on-site engineering projects, reduces operational risks, and improves operational efficiency.

[0084] In summary, the advantages proposed by this invention include at least the following:

[0085] (1) The method for determining the passage capacity of the well entry tool after casing deformation provided by the present invention effectively solves the problem of accurately calculating the passage capacity of the well entry string after casing deformation, and provides a strong basis for the selection decision of the well entry tool in casing deformation engineering.

[0086] (2) The method for determining the wellbore passing capacity after casing deformation provided by the present invention establishes a real and reliable three-dimensional wellbore with casing deformation, and the calculated wellbore passing capacity of the casing deformation tool string is accurate.

[0087] Although a method and system for determining the passability of a casing-deformed well tool has been described above in conjunction with exemplary embodiments, those skilled in the art will understand that various modifications and changes can be made to the exemplary embodiments of the present invention without departing from the spirit and scope defined by the claims.

Claims

1. A method for determining the wellbore penetration capability after casing deformation, characterized in that, The method includes the following steps: The center of the casing section is determined based on the two-dimensional spatial coordinates of the well diameter. Based on the acquired wellbore trajectory, depth, and caliber data, the three-dimensional spatial coordinates of the casing cross-section center and the endpoints of each arm are calculated, and a three-dimensional wellbore model after casing deformation is established; and The maximum through diameter of the tool in the casing deformation section was calculated using the projection method. The maximum through diameter of the tool in the deformed section of the casing is calculated using a three-dimensional wellbore model after casing deformation. The calculation of the maximum passing diameter of the tool in the deformed section of the casing using the projection method includes the following steps: Step 1: Set the projection plane, and place the first... The normal plane of the well trajectory at each depth point is set as the target projection plane, which is also of length [missing information]. L The bottom cut plane of the three-dimensional deformable wellbore, from the first From the depth point to the The coordinates of all caliper endpoints at a given depth point are set as the points to be projected, where... The depth sampling interval for the data conversion and merging unit; Step 2: Set the projection direction, and mesh the target projection plane to obtain a set of mesh nodes. Then the first Well trajectory coordinates at each depth point Pointing to a certain grid node Projection direction The expression includes Equation 4: ; Step 3: Calculate the projected coordinates and set them. A point on the well wall to be projected Coordinates are The normal to the target projection plane is , The projection point on the target projection plane is , where represents , The system of equations is as follows: ; Step 4: Obtain the inner boundary of the projection area, and set it from the first... From the depth point to the The wellbore at the depth point to be projected is at the... The projection point of the normal plane containing the well trajectory at depth n and the nth depth point The set consisting of wellhead endpoints at depth points is Get point set The set of inner boundary points D of the enclosed region; Step 5: Calculate the largest incircle of the inner boundary set D, denoted by its diameter. ; Step 6: Determine the length as L The maximum diameter through which a meter-sized tool can pass; Step 7: Repeat steps 2 to 6 to calculate the maximum inscribed circle diameter in the feasible projection direction, and determine the maximum value of all diameters as the length. L The maximum diameter of the meter tool .

2. The method for determining the wellbore penetration capability after casing deformation according to claim 1, characterized in that, The methods for determining the center of the casing cross-section include at least one of the centroid method, the least squares method, and the iterative search method.

3. The method for determining the wellbore penetration capability after casing deformation according to claim 1, characterized in that, The casing cross-section includes at least one of the casing interface at a specific depth and the casing deformation cross-section.

4. The method for determining the wellbore penetration capability after casing deformation according to claim 1, characterized in that, The method for calculating the three-dimensional spatial coordinates of the center of the casing cross-section and the endpoints of each arm includes a translational rotation replacement method. This method is used to reconstruct the well trajectory and image the wellbore. The well trajectory reconstruction process includes: Set parameters, set the first The three-dimensional coordinates of the center of the casing section at each depth point Translation vector That is, the first Original trajectory point coordinates at each depth point; rotation vector ,in, For the first Tangent of the trajectory point at each depth point for The normal to the plane, , and x The axis is parallel; the rotation angle For the first There are several depth points, with depth valued as Depth, Depth= h 0+ d k , k =0, 1, ..., m-1, d The depth sampling interval for the data conversion and merging unit; Translate the coordinates to move the point According to the translation vector Translate to point The calculation formula includes Equation 1: ; Rotate the coordinates to make the point Along the axis of rotation Rotation angle After that, I got points. The calculation formula includes Equation 2: ; Replace the coordinates with the central 3D spatial coordinates obtained after translation and rotation. Replace the Original trajectory coordinates of each depth point ; The wellbore imaging process includes: using the above-mentioned translation and rotation method to image the first... Coordinates of a depth point Perform rotation, initial number The three-dimensional coordinate vectors of the endpoints are: Its rotation axis along a rotation axis Rotation angle Then, the rotated vector is obtained. Its calculation formula includes Equation 3: ,Will Substitute them into the equation to obtain the rotated point set. Next, a displacement operation is performed, resulting in the set of points after displacement. for: .

5. The method for determining the wellbore penetration capability after casing deformation according to claim 1, characterized in that, The establishment of the three-dimensional wellbore model after casing deformation includes reconstructing the well trajectory and performing three-dimensional wellbore modeling.

6. The method for determining the wellbore penetration capability after casing deformation according to claim 1 or 5, characterized in that, The method for establishing a three-dimensional wellbore model after casing deformation includes data sample interpolation using the Lagrange interpolation method.

7. The method for determining the wellbore penetration capability after casing deformation according to claim 1, characterized in that, The method for calculating the maximum inscribed circle of the inner boundary set D includes an iterative search method.

8. A system for determining the passability of a wellbore insertion tool after casing deformation, characterized in that, The system enables the determination of the wellbore penetration capability after casing deformation as described in any one of claims 1 to 7, and the system comprises: The data conversion and merging unit is configured to integrate the original well trajectory two-dimensional coordinate data and the two-dimensional spatial coordinates of the well diameter after casing deformation; The wellbore data correction unit is configured to determine the center of the casing cross-section and compare the obtained centers; The three-dimensional wellbore modeling unit is configured as the three-dimensional spatial coordinates of the calculation center and each arm endpoint to reconstruct the well trajectory and establish a three-dimensional wellbore model after casing deformation. The capability calculation unit is configured to calculate the maximum through diameter of the tool in the casing deformation section.

9. The system for determining the wellbore penetration capability after casing deformation according to claim 8, characterized in that, The data conversion and merging unit includes at least one of the following: well trajectory data, well diameter logging data, casing inner diameter, and depth sampling interval input modules; The wellbore data correction unit includes at least one of the following: raw data calculation and storage, geometric center method data calculation and storage, least squares method data calculation and storage, iterative search method data calculation and storage, and a calculation center two-dimensional drawing module; The three-dimensional wellbore modeling unit includes at least one of the following modules: well section setting, well trajectory three-dimensional space display, and well wall three-dimensional space modeling module; The throughput calculation unit includes at least one of the following: tool length input, calculation step setting, and projection method calculation module.

10. A computer device, characterized in that, The device includes: processor; and The memory stores a computer program that, when executed by a processor, implements the method for determining the wellbore penetration capability after casing deformation as described in any one of claims 1 to 7.

11. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for determining the well-entry tool's passability after casing deformation as described in any one of claims 1 to 7.

Images