Wellbore multiphase pipe flow calculation method and device based on well trajectory point thinning

Abstract

The application discloses a kind of wellbore multiphase pipe flow calculation method and device based on well trajectory point thinning. The method comprises: obtaining original well trajectory data points, retaining the first and last data points, preliminarily screening the intermediate data points, and secondarily screening the preliminarily screened intermediate data points, which includes: calculating the difference between each intermediate data point and the adjacent point, comparing with the preset first filtering threshold, removing the points with a difference lower than the first filtering threshold, determining whether the number of the first data point, the intermediate data points after secondary screening and the last data point is less than or equal to the preset optimized number, and if not, adjusting the first filtering threshold and re-screening until the number meets the preset optimized number, and outputting the first data point, the intermediate data points after secondary screening and the last data point. The application effectively reduces the number of data points while ensuring the consistency of well trajectory geometric features, significantly reduces the calculation time in the multiphase flow simulation process and reduces the consumption of computing resources.

Description

Technical Field

[0001] This invention relates to the field of oil and gas extraction, and in particular to a method and apparatus for calculating multiphase flow in wellbore based on well trajectory point thinning. Background Technology

[0002] In recent years, with the continuous growth of global energy demand and the development of artificial intelligence, the development of oil and gas resources has become increasingly reliant on efficient production technologies and optimization schemes. Wellbore multiphase flow refers to the flow of fluids in oil wells, gas wells, and other wells where two or more phases (such as oil, gas, and water) coexist. This flow pattern is very common in oil and gas extraction. Fluids of different phases possess different physical properties, such as density and viscosity, and their flow patterns within the wellbore are complex and influenced by various factors. Wellbore multiphase flow simulation calculations mainly involve simulating temperature, pressure, and flow regime along the wellbore, and are a crucial step in oil and gas field development and production. The flow characteristics within the wellbore directly affect production efficiency and development safety assessment. However, due to the nonlinearity, complexity, and uncertainty of multiphase flow problems, there is often a trade-off between the accuracy and computational efficiency of wellbore multiphase flow calculations.

[0003] Currently, the traditional wellbore multiphase flow calculation process takes significantly longer to calculate when dealing with complex well trajectories, resulting in slow calculation speed and low efficiency, making it difficult to meet the real-time decision analysis needs of actual production processes. Summary of the Invention

[0004] In view of the above problems, the present invention is proposed to provide a method and apparatus for calculating multiphase flow in wellbore based on well trajectory point thinning to overcome or at least partially solve the above problems.

[0005] In a first aspect, embodiments of the present invention provide a method for thinning well trajectory points, comprising:

[0006] Obtain the original well trajectory data points;

[0007] The first and last data points of the original well trajectory data points are retained, and the intermediate data points other than the first and last data points are initially screened.

[0008] For each intermediate data point obtained after the initial screening, a secondary screening is performed; the secondary screening includes: calculating the difference between each intermediate data point and the previous intermediate data point of the adjacent well trajectory, comparing the calculated difference with a preset first filtering threshold, and removing intermediate data points whose absolute value of the difference is lower than the first filtering threshold.

[0009] Determine whether the number of the first data point, the intermediate data points after secondary filtering, and the tail data points is less than or equal to the preset number of optimized well trajectory data points;

[0010] If the requirement cannot be met, the first filtering threshold of the secondary screening is adjusted, and each intermediate data point is re-screened until the number of each intermediate data point, first data point, and last data point after the secondary screening is less than or equal to the preset number of optimized well trajectory data points.

[0011] Output the first data point, each intermediate data point after secondary filtering, and the last data point.

[0012] In one embodiment, the data for each well trajectory data point includes: the well inclination angle at the well trajectory measurement location.

[0013] In one embodiment, the initial screening of all intermediate data points other than the first and last data points includes:

[0014] The following calculations are performed using the well inclination angle of the intermediate data points:

[0015] Compare the well inclination angle of each intermediate data point with the second filter threshold of the preset well inclination angle;

[0016] If the well inclination angle is less than or equal to the second filtering threshold of the well inclination angle, then the current intermediate data point is removed; otherwise, the intermediate data point is retained.

[0017] In one embodiment, during the secondary filtering, the difference between each intermediate data point and the previous intermediate data point of the adjacent well trajectory is calculated. The calculated difference is compared with a preset filtering threshold, and intermediate data points whose absolute value of the difference is lower than the first filtering threshold are removed, including:

[0018] The inclination angle of each intermediate data point is compared with the inclination angle of the previous intermediate data point in the adjacent well trajectory. The absolute value of the calculated inclination angle difference is compared with a preset first filtering threshold for inclination angle. If the absolute value of the inclination angle difference corresponding to the intermediate data point is less than or equal to the first filtering threshold, the intermediate data point is removed.

[0019] In one embodiment, adjusting the first filtering threshold of the secondary screening includes:

[0020] The first filtering threshold is adjusted by calculating α = α + η;

[0021] Where η is a constant representing the magnitude of change of the first filtering threshold α;

[0022] When n is greater than m, η takes a negative value; conversely, when n is less than m, η takes a positive value; n is the total number of the first data point, the intermediate data point and the tail data point obtained by the second round of screening, and m is the number of the preset optimized well trajectory data points.

[0023] Secondly, embodiments of the present invention provide a method for calculating multiphase pipe flow in a wellbore based on well trajectory point thinning, including:

[0024] This is used to thin out the collected raw well trajectory data points using the well trajectory point thinning method described above.

[0025] Using the data from the well trajectory points after thinning, the calculation of multiphase flow in the wellbore is performed.

[0026] Thirdly, embodiments of the present invention provide a dilution device for well trajectory points, comprising:

[0027] The acquisition module is used to acquire raw well trajectory data points;

[0028] The initial screening module is used to retain the first and last data points of the original well trajectory data points and to perform initial screening on all intermediate data points other than the first and last data points.

[0029] A secondary filtering module is used to further filter each intermediate data point obtained after the initial screening. The secondary filtering includes: calculating the difference between each intermediate data point and the previous intermediate data point of the adjacent well trajectory; comparing the calculated difference with a preset first filtering threshold; removing intermediate data points whose absolute value of the difference is lower than the first filtering threshold; determining whether the number of the first data point, the intermediate data points after secondary filtering, and the tail data point is less than or equal to the preset number of optimized well trajectory data points; if not, adjusting the first filtering threshold of the secondary filtering and re-filtering each intermediate data point until the number of intermediate data points, the first data point, and the tail data point after secondary filtering meets the requirement.

[0030] The output module is used to output the first data point, each intermediate data point after secondary filtering, and the last data point.

[0031] Fourthly, embodiments of the present invention provide a wellbore multiphase pipe flow calculation device based on well trajectory point thinning, comprising:

[0032] The preprocessing module is used to thin out the collected raw well trajectory data points using the well trajectory point thinning method described above.

[0033] The wellbore multiphase flow calculation module is used to perform wellbore multiphase flow calculations using data from well trajectory points after thinning.

[0034] Fifthly, embodiments of the present invention provide a computing device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the program executed by the processor implements either the well trajectory point thinning method described above, or the wellbore multiphase pipe flow calculation method described above.

[0035] In a sixth aspect, embodiments of the present invention provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements either the well trajectory point thinning method described above or the wellbore multiphase pipe flow calculation method described above.

[0036] The beneficial effects of the above-described technical solutions provided in the embodiments of the present invention include at least the following:

[0037] The above-mentioned wellbore multiphase flow calculation method and apparatus based on well trajectory point thinning provided in this invention calculates the difference between adjacent well trajectory data points and uses the change in adjacent data points as a filtering reference index to perform thinning calculation on the collected original well trajectory data points, filtering out data points that have not undergone significant changes. Thus, while ensuring that the geometric characteristics of the well trajectory are basically consistent, the number of data points is effectively reduced. This significantly reduces the calculation time in the multiphase flow simulation process, improves calculation efficiency, reduces the consumption of computing resources and calculation costs, and achieves efficient simulation of temperature, pressure and flow regime along the wellbore while ensuring the accuracy of the multiphase flow simulation results.

[0038] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.

[0039] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0040] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings:

[0041] Figure 1 A flowchart of the well trajectory point thinning method provided in an embodiment of the present invention;

[0042] Figure 2 Distribution map of the original well trajectory data provided in the embodiments of the present invention;

[0043] Figure 3This is a schematic diagram showing the distribution of well trajectory data points obtained by the well trajectory point thinning method provided in an embodiment of the present invention;

[0044] Figure 4 A comparison of computation time before and after implementing the well trajectory point thinning method provided in this embodiment of the invention for a complex well trajectory;

[0045] Figure 5 A comparison diagram of the geometric feature preservation degree before and after implementing the well trajectory point thinning method provided in this embodiment of the invention;

[0046] Figure 6 A line graph showing the simulation results of well trajectory pressure before and after implementing the well trajectory point thinning method provided in this embodiment of the invention;

[0047] Figure 7 A line graph showing the simulation results of well trajectory temperature before and after the implementation of the well trajectory point thinning method provided in this embodiment of the invention;

[0048] Figure 8 A flowchart of a wellbore multiphase pipe flow calculation method provided in an embodiment of the present invention;

[0049] Figure 9 This is a structural block diagram of the well trajectory point dilution device provided in an embodiment of the present invention;

[0050] Figure 10 This is a structural block diagram of a wellbore multiphase pipe flow calculation device provided in an embodiment of the present invention. Detailed Implementation

[0051] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0052] The inventors of this invention discovered that, in order to solve the problem that the traditional wellbore multiphase flow calculation process in the prior art increases the calculation time significantly when facing complex well trajectories, resulting in slow calculation speed and low efficiency, and making it difficult to meet the real-time decision analysis in actual production processes, many petroleum development workers have proposed various well trajectory thinning methods based on data point density, curvature, etc., to improve the multiphase flow calculation speed. However, in practical applications, it is still difficult to balance accuracy and efficiency.

[0053] To find a balance between calculation speed and accuracy for multiphase flow in wellbore pipes, this invention provides a method for thinning well trajectory points, referring to... Figure 1 As shown, it includes the following steps:

[0054] S11. Obtain the original well trajectory data points;

[0055] S12. Retain the first and last data points of the original well trajectory data points, and perform preliminary screening on each intermediate data point other than the first and last data points.

[0056] The initial screening of intermediate data points (excluding the first and last data points) in step S12 can be implemented through the following steps:

[0057] The following calculations are performed using the well inclination angle of the intermediate data points:

[0058] Compare the well inclination angle of each intermediate data point with the second filter threshold of the preset well inclination angle;

[0059] If the well inclination angle is less than or equal to the second filtering threshold of the well inclination angle, then remove the current intermediate data point; otherwise, retain the intermediate data point.

[0060] The reason for retaining the first and last data points of the well trajectory while filtering the intermediate data points is that the first and last data points can ensure that the total length of the original well trajectory points remains unchanged, thus avoiding huge errors after filtering them out.

[0061] S13. For each intermediate data point obtained after the initial screening, a second screening is performed. The second screening is as follows: the difference between each intermediate data point obtained after the initial screening and the previous intermediate data point of the adjacent well trajectory is calculated, the calculated difference is compared with a preset first filtering threshold, and intermediate data points whose absolute value of the difference is lower than the first filtering threshold are removed.

[0062] S14. Determine whether the number of the first data point, the intermediate data points after secondary filtering, and the last data point is less than or equal to the preset number of optimized well trajectory data points; if not, proceed to step S15 below.

[0063] S15. Adjust the first filtering threshold of the secondary screening and re-execute each step of the secondary screening of each intermediate data point until the number of each intermediate data point, first data point and last data point after secondary screening can meet the requirement of being less than or equal to the preset number of optimized well trajectory data points, then execute S16.

[0064] S16. Output the first data point, each intermediate data point after secondary filtering, and the last data point.

[0065] In this embodiment of the invention, the well trajectory refers to the trajectory formed from the start of drilling to the end of drilling (reaching the bottom of the well). The trajectory consists of a series of trajectory points. The first point of the trajectory corresponds to the relevant data of the position where the drill bit starts drilling, i.e., the aforementioned "first data point". This position is usually located at the starting position of drilling on the drilling platform. The last point of the trajectory represents the relevant data of the position of the drill bit at the end of drilling, i.e., the aforementioned "tail data point". The data points located between the "first data point" and the "tail data point" are called "intermediate data points".

[0066] The well trajectory point thinning method provided in this embodiment calculates the difference between adjacent well trajectory data points and uses the change in adjacent data points as a filtering reference index to perform thinning calculation on the collected original well trajectory data points, filtering out data points that have not undergone significant changes. Thus, while ensuring that the geometric characteristics of the well trajectory are basically consistent, the number of data points is effectively reduced. This significantly reduces the calculation time in the multiphase flow simulation process, improves calculation efficiency, reduces the consumption of computing resources and calculation costs, and achieves efficient simulation of temperature, pressure and flow regime along the wellbore while ensuring the accuracy of the multiphase flow simulation results.

[0067] Furthermore, in the above-mentioned well trajectory point thinning method, the data for each well trajectory data point may include: the measured well depth, vertical depth, well inclination angle, and azimuth angle of the well trajectory measurement location.

[0068] The wellbore inclination angle is the angle between the wellbore axis and the plumb line. Simply put, it measures the degree to which the wellbore deviates from its vertical direction. If the wellbore is perfectly vertical downwards, the inclination angle is 0°; as the wellbore begins to tilt, the inclination angle becomes greater than 0°. For example, during the drilling of deviated or directional wells, the inclination angle gradually increases as the drilling direction changes.

[0069] Well depth refers to the actual length from the wellhead to the bottom of the well, measured along the centerline of the wellbore. It is a comprehensive indicator that takes into account the length of the wellbore in both vertical and horizontal directions. Regardless of whether the wellbore is vertical, inclined, or horizontal, the well depth reflects the length of the path traveled by the drill bit throughout the drilling process. For example, in a complex well with a horizontal section, the well depth would include the sum of the lengths of the vertical, inclined, and horizontal sections.

[0070] Vertical depth refers to the vertical distance from the wellhead to the bottom of the well, which is the depth in the direction of gravity. It only considers the depth of the wellbore in the vertical direction and does not take into account any deviation of the wellbore in the horizontal direction. For example, in a vertical well, the well depth and vertical depth are equal; but in an inclined or horizontal well, the vertical depth is less than the well depth.

[0071] Well trajectory points are discrete locations of the wellbore trajectory in three-dimensional space. These points record the specific spatial location of the well at different depths, including information such as vertical depth, horizontal displacement, and orientation. Each trajectory point is like a "landmark" on the wellbore path, and multiple such points can be used to depict the shape of the entire well trajectory.

[0072] In this embodiment of the invention, the aforementioned angle value, such as the well inclination angle, can be used as reference data for initial screening and secondary screening of data points.

[0073] In this embodiment of the invention, the filtering threshold used in the primary screening and the secondary screening are referred to as the second filtering threshold only to distinguish them. The filtering threshold used in the primary screening is referred to as the first filtering threshold. The terms "first" and "second" are used only for distinction and do not limit their importance or priority.

[0074] In one embodiment, during the secondary screening, the difference between the well inclination angle of each intermediate data point and the well inclination angle of the previous intermediate data point of the adjacent well trajectory is calculated. The absolute value of the calculated well inclination angle difference is compared with a preset first filtering threshold for well inclination angle. If the absolute value of the well inclination angle difference corresponding to the intermediate data point is less than or equal to the first filtering threshold, the intermediate data point is removed.

[0075] For example:

[0076] For every two adjacent well trajectory data points, the difference in well inclination angle between them is calculated using equation (1).

[0077] Δθ i =θ i+1 -θ i ;

[0078] In the above formula, Δθ i θ represents the difference in wellbore inclination angle between the i-th data point and the (i+1)-th data point; i θ represents the well inclination angle of the i-th data point. i+1 It represents the well inclination angle of the (i+1)th data point.

[0079] Assuming the first filtering threshold for the well inclination angle is set to α, when judging |Δθ i When |>α, retain θ i+1 Data points are selected; otherwise, θ is filtered out. i+1 Data points.

[0080] The principle for secondary screening of intermediate data points is to retain intermediate data points with significant changes in well inclination angle and to filter out those intermediate data points with insignificant changes in well inclination angle.

[0081] This invention, through calculating the difference in well inclination angle between adjacent data points and using the change in well inclination angle (or orientation angle) as a filtering reference index, can significantly reduce the number of well trajectory data points while maintaining the basic geometric characteristics of the well trajectory, thereby improving the efficiency of multiphase pipe flow calculation. Using this method, computation time and resource consumption can be significantly reduced while preserving the geometric characteristics of the well trajectory, achieving efficient simulation of temperature, pressure, and flow regime along the wellbore.

[0082] The method provided in this invention is applicable to well trajectory thinning requirements of varying complexity, exhibiting wide applicability and good versatility. By reducing the number of well trajectory data points, it significantly reduces computational costs while maintaining accuracy, thereby improving the accuracy and reliability of multiphase flow calculations.

[0083] In one embodiment, the first filtering threshold of the secondary screening is adjusted in step S15. Specifically, the first filtering threshold can be adjusted by calculating α = α + η; α is the first filtering threshold of the well inclination angle, and η is a constant of the change range of the angle threshold α.

[0084] Let n be the total number of the first data point, the intermediate data points, and the last data point obtained from the second round of screening, and m be the preset number of optimized well trajectory data points. When n is greater than m, it indicates that the current step size setting may be too large, resulting in too many points after thinning, and the step size needs to be reduced. In this case, η takes a negative value. Conversely, when n is less than m, it indicates that the step size may be too small, and the step size needs to be increased, and η takes a positive value. The specific value of η can be adjusted according to the actual situation and experience to ensure that the final number of well trajectory points after thinning meets the requirements.

[0085] In one embodiment, the number of the aforementioned preset optimized well trajectory data points can be determined in advance based on computation time and resource consumption.

[0086] Let's illustrate this with an example.

[0087] The large number of well trajectory data points leads to low computational efficiency in multiphase pipe flow calculations. Existing thinning methods, when handling complex well trajectories, have long computation times and are prone to losing key geometric features, resulting in low accuracy. For example, a complex single-well trajectory in a certain block contains 241 data points, and the multiphase pipe flow calculation in this case takes 80 seconds. (Refer to...) Figure 2 The distribution map of the original well trajectory data shown is as follows. Figure 2 The original data point distribution of a complex well trajectory is shown. A total of 241 data points were collected in the original well trajectory, which shows the characteristics of the high-density data point distribution of the well trajectory.

[0088] By calculating the difference in well inclination angle between well trajectory data points and using the change in well inclination angle as a filtering reference index for thinning calculations, the number of well trajectory data points in the embodiment is reduced to 37 after processing, while maintaining the original well trajectory shape to the greatest extent possible. This significantly improves the calculation efficiency of multiphase pipe flow. Actual tests show that the algorithm of this invention reduces the number of data points by approximately 85%, shortens the calculation time to 10 seconds, maintains a basically consistent geometric appearance, and the multiphase flow calculation simulation results after thinning are basically consistent with the calculation results before thinning, with an error of less than 1%. By implementing the embodiment of this invention, the number of well trajectory data points can be significantly reduced, improving the efficiency of multiphase pipe flow calculations and providing an efficient and reliable solution for wellbore multiphase flow calculation simulations.

[0089] The distribution of well trajectory data points obtained by the well trajectory point thinning method provided in this embodiment of the invention can be referred to... Figure 3 As shown, the number of data points was reduced to 37, the main geometric features of the well trajectory were preserved, and the thinning effect was significant.

[0090] Figure 4 The invention demonstrates a comparison of computation time before and after implementing the well trajectory point thinning method provided in this embodiment for a complex well trajectory. After implementing the well trajectory point thinning method provided in this embodiment, the computation time was reduced from 80 seconds to 10 seconds, and the computation efficiency was improved by 8 times.

[0091] Figure 5 The geometric feature retention is compared before and after implementing the well trajectory point thinning method provided in the embodiments of the present invention, as shown in the form of a line graph, proving that the geometric feature retention is basically consistent after implementing the well trajectory point thinning method provided in the embodiments of the present invention.

[0092] Figure 6 The simulation results of well trajectory pressure before and after implementation are presented in the form of a line graph, proving that the simulation results are basically consistent after implementing the present invention.

[0093] Figure 7 The simulation results of well trajectory temperature before and after implementation are shown in the form of a line graph, proving that the simulation results are basically consistent after implementing the present invention.

[0094] The embodiments of the present invention are applicable to well trajectory thinning requirements of varying complexity and can be stably applied under various wellbore conditions, further improving the speed and accuracy of multiphase flow calculation.

[0095] This invention also provides a method for calculating multiphase pipe flow in wellbore based on well trajectory point thinning, referring to... Figure 8 As shown, the method includes the following steps:

[0096] S81, used to thin out the collected original well trajectory data points using the well trajectory point thinning method as described in the foregoing embodiments;

[0097] The specific implementation process of this step can be referred to the aforementioned embodiments, and will not be repeated here.

[0098] S82. Using the data from the well trajectory data points after thinning, perform the calculation of multiphase flow in the wellbore.

[0099] The calculation of multiphase flow in wellbore can be carried out using various existing calculation methods. The general steps include: determining the basic data and parameters of fluid and wellbore, selecting a suitable calculation model, calculating friction, gravity, and acceleration pressure gradient, determining the flow pattern, iteratively calculating the pressure distribution, calculating other parameters, and finally analyzing and verifying the results. The specific implementation process will not be detailed here.

[0100] Based on the same inventive concept, embodiments of the present invention also provide a well trajectory point thinning device and a wellbore multiphase pipe flow calculation device based on well trajectory point thinning. Since the principle of solving the problem by these devices is similar to the aforementioned well trajectory point thinning method and wellbore multiphase pipe flow calculation method, the implementation of the device can refer to the implementation of the aforementioned method, and the repeated parts will not be described again.

[0101] An embodiment of the present invention provides a dilution device for well trajectory points, referring to... Figure 9 As shown, it includes:

[0102] Module 91 is used to acquire raw well trajectory data points;

[0103] The initial screening module 92 is used to retain the first and last data points of the original well trajectory data points and to perform initial screening on each intermediate data point other than the first and last data points.

[0104] The secondary screening module 93 is used to further screen each intermediate data point obtained after the initial screening. The secondary screening includes: calculating the difference between each intermediate data point and the previous intermediate data point of the adjacent well trajectory; comparing the calculated difference with a preset first filtering threshold; removing intermediate data points whose absolute value of the difference is lower than the first filtering threshold; determining whether the number of the first data point, the intermediate data points after secondary screening, and the tail data points is less than or equal to the preset number of optimized well trajectory data points; if not, adjusting the first filtering threshold of the secondary screening and re-screening each intermediate data point until the number of intermediate data points, the first data point, and the tail data points after secondary screening meets the requirement.

[0105] Output module 94 is used to output the first data point, each intermediate data point after secondary filtering, and the last data point.

[0106] This invention provides a wellbore multiphase pipe flow calculation device based on well trajectory point thinning, referring to... Figure 10 As shown, it includes:

[0107] The preprocessing module 101 is used to thin out the collected original well trajectory data points using the well trajectory point thinning method described in the foregoing embodiments.

[0108] The calculation module 102 for multiphase flow in the wellbore is used to perform calculations of multiphase flow in the wellbore using data from well trajectory data points after thinning.

[0109] This invention also provides a computing device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor. The program executed by the processor implements either the well trajectory point thinning method described above or the wellbore multiphase pipe flow calculation method described above.

[0110] This invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements either the well trajectory point thinning method described above or the wellbore multiphase pipe flow calculation method described above.

[0111] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage and optical storage) containing computer-usable program code.

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

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

[0114] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the functions specified in one or more boxes. Obviously, those skilled in the art can make various modifications and variations to this invention without departing from the spirit and scope of the invention. Therefore, if these modifications and variations of the invention fall within the scope of the claims of the invention and their equivalents, the invention is also intended to include these modifications and variations.

Claims

1. A method for thinning well trajectory points, characterized in that, include: Obtain the original well trajectory data points; The first and last data points of the original well trajectory data points are retained, and the intermediate data points other than the first and last data points are initially screened. Each intermediate data point obtained after the initial screening is further screened a second time. The secondary screening includes: calculating the difference between each intermediate data point and the previous intermediate data point of the adjacent well trajectory; comparing the calculated difference with a preset first filtering threshold; and removing intermediate data points whose absolute value of the difference is lower than the first filtering threshold. Determine whether the number of the first data point, the intermediate data points after secondary filtering, and the tail data points is less than or equal to the preset number of optimized well trajectory data points; If the requirement cannot be met, the first filtering threshold of the secondary screening is adjusted, and each intermediate data point is re-screened until the number of each intermediate data point, first data point, and last data point after the secondary screening is less than or equal to the preset number of optimized well trajectory data points. Output the first data point, each intermediate data point after secondary filtering, and the last data point.

2. The method as described in claim 1, characterized in that, The data for each well trajectory data point includes: the well inclination angle at the well trajectory measurement location.

3. The method as described in claim 2, characterized in that, The initial screening of all intermediate data points other than the first and last data points includes: The following calculations are performed using the inclination angle of the intermediate data points: Compare the well inclination angle of each intermediate data point with the second filter threshold of the preset well inclination angle; If the well inclination angle is less than or equal to the second filtering threshold of the well inclination angle, then the current intermediate data point is removed; otherwise, the intermediate data point is retained.

4. The method as described in claim 2, characterized in that, In the secondary screening, the difference between each intermediate data point and the previous intermediate data point of the adjacent well trajectory is calculated. The calculated difference is compared with a preset first filtering threshold. Intermediate data points whose absolute difference is lower than the first filtering threshold are removed, including: The inclination angle of each intermediate data point is compared with the inclination angle of the previous intermediate data point in the adjacent well trajectory. The absolute value of the calculated inclination angle difference is compared with a preset first filtering threshold for inclination angle. If the absolute value of the inclination angle difference corresponding to the intermediate data point is less than or equal to the first filtering threshold, the intermediate data point is removed.

5. The method according to any one of claims 1-4, characterized in that, Adjusting the first filtering threshold of the secondary screening includes: The first filtering threshold is adjusted by calculating α = α + η; Where η is a constant representing the magnitude of change of the first filtering threshold α; When n is greater than m, η takes a negative value; conversely, when n is less than m, η takes a positive value; n is the total number of the first data point, the intermediate data point and the tail data point obtained by the second round of screening, and m is the number of the preset optimized well trajectory data points.

6. A method for calculating multiphase pipe flow in a wellbore based on well trajectory point thinning, characterized in that, include: Used to thin out the original well trajectory data points acquired using the well trajectory point thinning method as described in any one of claims 1-5; Using the data from the well trajectory points after thinning, the calculation of multiphase flow in the wellbore is performed.

7. A dilution device for well trajectory points, characterized in that, include: The acquisition module is used to acquire raw well trajectory data points; The initial screening module is used to retain the first and last data points of the original well trajectory data points and to perform initial screening on all intermediate data points other than the first and last data points. The secondary filtering module is used to further filter each intermediate data point obtained after the initial screening. The secondary filtering includes: calculating the difference between each intermediate data point and the previous intermediate data point of the adjacent well trajectory; comparing the calculated difference with a preset first filtering threshold; removing intermediate data points whose absolute difference is lower than the first filtering threshold; determining whether the number of the first data point, the intermediate data points after secondary filtering, and the tail data point is less than or equal to the preset number of optimized well trajectory data points; if not, adjusting the first filtering threshold of the secondary filtering and re-filtering each intermediate data point until the number of each intermediate data point, the first data point, and the tail data point after secondary filtering meets the requirement; The output module is used to output the first data point, each intermediate data point after secondary filtering, and the last data point.

8. A wellbore multiphase pipe flow calculation device based on well trajectory point thinning, characterized in that, include: The preprocessing module is used to thin out the collected original well trajectory data points using the well trajectory point thinning method as described in any one of claims 1-5. The wellbore multiphase flow calculation module is used to perform wellbore multiphase flow calculations using data from well trajectory points after thinning.

9. A computing device, characterized in that, include: The memory, the processor, and the computer program stored in the memory and executable on the processor, wherein the program executed by the processor implements the well trajectory point thinning method as described in any one of claims 1-5, or implements the wellbore multiphase pipe flow calculation method as described in claim 6.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the well trajectory point thinning method as described in any one of claims 1-5, or the wellbore multiphase pipe flow calculation method as described in claim 6.

Images