Well track determination method and device of horizontal well and electronic equipment

Abstract

The invention discloses a borehole trajectory determination method and device of a horizontal well and electronic equipment. The method comprises the steps that wellhead coordinates and target coordinates corresponding to a well track of a horizontal well are obtained, and the target coordinates are used for indicating coordinates corresponding to well drilling points in a reservoir; determining a value range of borehole trajectory parameters of the borehole trajectory based on the wellhead coordinates and the target point coordinates; establishing a borehole trajectory model based on the value range of the borehole trajectory parameters; on the basis of the borehole trajectory model, borehole trajectory information of the horizontal well is obtained, and the borehole trajectory information comprises a borehole trajectory and the trajectory length and the friction resistance value corresponding to the borehole trajectory. According to the method, the technical problems of low borehole trajectory determination efficiency and high cost of a borehole trajectory determination method of a horizontal well in the related technology are solved.

Description

Technical Field

[0001] The present invention relates to the field of drilling trajectory design, and in particular to a method, device and electronic equipment for determining a wellbore trajectory of a horizontal well. Background Art

[0002] Drilling engineering is a key link in establishing a channel connecting the ground and the reservoir and developing oil and gas resources. With the continuous improvement of the development of conventional oil and gas resources, unconventional resources such as shale oil and tight oil and gas have gradually become important replacement resources and the main development direction in the future. At present, unconventional resources such as shale oil and tight oil and gas are mainly developed by cluster well groups and horizontal wells. Due to the restrictions of ground conditions such as pastoral areas, residential areas, and environmental protection areas, the offset distance and horizontal section length of horizontal wells are constantly increasing, and the difficulty and cost of drilling engineering have increased significantly. The drilling cost can account for more than 50% of the well construction cost. For horizontal wells, excessive friction during drilling is a key factor that affects the extension of the horizontal section and hinders the smooth lowering of the casing in the later stage.

[0003] The horizontal wellbore trajectory determination method in related technologies does not comprehensively consider the target location and engineering condition restrictions, and the trajectory designed using conventional software cannot directly and accurately measure its quality, and requires continuous trajectory comparison and adjustment, resulting in low design efficiency. The designers have a heavy workload and have to pay high software usage and maintenance fees to software companies every year, resulting in high design costs.

[0004] To address the above-mentioned problems, no effective solution has been proposed yet. Summary of the invention

[0005] The embodiments of the present invention provide a method, device and electronic equipment for determining a wellbore trajectory of a horizontal well, so as to at least solve the technical problems of low wellbore trajectory determination efficiency and high cost existing in the wellbore trajectory determination method of a horizontal well in the related art.

[0006] According to one aspect of an embodiment of the present invention, a method for determining a wellbore trajectory of a horizontal well is provided, comprising: obtaining wellhead coordinates and target coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target coordinates are used to indicate coordinates corresponding to a drilling point in a reservoir; determining a value range of wellbore trajectory parameters of the wellbore trajectory based on the wellhead coordinates and the target coordinates; establishing a wellbore trajectory model based on the value range of the wellbore trajectory parameters; and obtaining wellbore trajectory information of the horizontal well based on the wellbore trajectory model, wherein the wellbore trajectory information includes the wellbore trajectory, and a trajectory length and a friction value corresponding to the wellbore trajectory.

[0007] According to another aspect of an embodiment of the present invention, a device for determining a wellbore trajectory of a horizontal well is provided, comprising: a first acquisition module, used to acquire wellhead coordinates and target point coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target point coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir; a first determination module, used to determine a value range of wellbore trajectory parameters of the wellbore trajectory based on the wellhead coordinates and the target point coordinates; a first establishment module, used to establish a wellbore trajectory model based on the value range of the wellbore trajectory parameters; and a second determination module, used to obtain the wellbore trajectory information of the horizontal well based on the wellbore trajectory model, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory.

[0008] According to another aspect of an embodiment of the present invention, there is also provided an electronic device, comprising one or more processors and a memory, wherein the memory is used to store one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement any one of the methods for determining the wellbore trajectory of a horizontal well.

[0009] In an embodiment of the present invention, by obtaining the wellhead coordinates and target coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir; based on the wellhead coordinates and the target coordinates, the value range of the wellbore trajectory parameters of the wellbore trajectory is determined; based on the value range of the wellbore trajectory parameters, a wellbore trajectory model is established; based on the wellbore trajectory model, the wellbore trajectory information of the horizontal well is obtained, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory, the purpose of establishing a wellbore trajectory model based on the value range of the parameters corresponding to the wellbore trajectory is achieved, so as to efficiently obtain the wellbore trajectory determination result of the horizontal well, thereby achieving the technical effect of improving the wellbore trajectory determination efficiency of the horizontal well and reducing the wellbore trajectory determination cost of the horizontal well, thereby solving the technical problems of low wellbore trajectory determination efficiency and high cost in the wellbore trajectory determination method of the horizontal well in the related art. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] The drawings described herein are used to provide a further understanding of the present invention and constitute a part of this application. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0011] Figure 1 is a schematic diagram of a method for determining a wellbore trajectory of a horizontal well according to an embodiment of the present invention;

[0012] Figure 2 is a schematic diagram of an optional wellbore trajectory design according to an embodiment of the present invention;

[0013] Figure 3 is a schematic diagram of an optional inclined arc interpolation method according to an embodiment of the present invention;

[0014] Figure 4 is a schematic diagram of an optional coordinate system according to an embodiment of the present invention;

[0015] Figure 5 is a schematic diagram of an optional force analysis according to an embodiment of the present invention;

[0016] Figure 6 is a flow chart of an optional wellbore trajectory model solution according to an embodiment of the present invention;

[0017] Figure 7 is a schematic diagram of a wellbore trajectory determination device for a horizontal well according to an embodiment of the present invention;

[0018] Figure 8 is a schematic diagram of an electronic device according to an embodiment of the present invention. DETAILED DESCRIPTION

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

[0020] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0021] According to an embodiment of the present invention, a method embodiment for determining the wellbore trajectory of a horizontal well is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0022] Figure 1 is a flow chart of a method for determining a wellbore trajectory of a horizontal well according to an embodiment of the present invention. Figure 1 As shown, the method comprises the following steps:

[0023] Step S102, obtaining the wellhead coordinates and target point coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target point coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir.

[0024] Optional, Figure 2 is a schematic diagram of an optional wellbore trajectory design according to an embodiment of the present invention, such as Figure 2 As shown, the wellhead O coordinate (N O ,E O ,H O ), where N0 is the north-south displacement of point O, E0 is the east-west displacement of point O, and H0 is the vertical depth of point O. The coordinates of the first target point T0 can be determined using geological exploration data or geological models. in, is the north-south displacement of point T0, is the east-west displacement of point T0, is the vertical depth of point T0, and the coordinates of the final target point T1: in, is the north-south displacement of point T1, is the east-west displacement of point T1, is the vertical depth of point T1. By obtaining the wellhead coordinates and target coordinates, the geographical location of the horizontal well can be determined, providing accurate positioning information to facilitate operations and management on the ground; the shape and path of the wellbore trajectory can be determined for engineering design, drilling operations and production management; it can also be used to evaluate the quality of the wellbore, for example, by comparing the difference between the actual wellbore trajectory and the designed wellbore trajectory, the degree of deviation and quality of the wellbore can be evaluated; it can also guide drilling operations, determine drilling plans and operating methods, and ensure the smooth progress of drilling operations.

[0025] Step S104, determining the value range of the wellbore trajectory parameters of the wellbore trajectory based on the wellhead coordinates and the target point coordinates.

[0026] Optionally, based on the wellhead coordinates, target coordinates, and taking into account the influence of formation stability, collision avoidance, construction capacity, engineering conditions and other factors, the value range of the wellbore trajectory parameters can be reasonably determined to ensure the stability of the wellbore trajectory and the accuracy of the drilling target during the drilling process.

[0027] In an optional embodiment, based on the wellhead coordinates and the target point coordinates, the value range of the wellbore trajectory parameters of the wellbore trajectory is determined, including: segmenting the wellbore trajectory to obtain multiple wellbore trajectory segments; based on the wellhead coordinates and the target point coordinates, determining the value ranges of the wellbore trajectory parameters corresponding to the multiple wellbore trajectory segments.

[0028] Optionally, the borehole trajectory can be divided into 5 sections including: straight-increasing-stable-increasing and twisting-horizontal section, or the borehole trajectory can be divided into 7 sections including: straight well section, first deflection section, first stable deflection section, second deflection section, second stable deflection section, third deflection section, horizontal section. Since each section of the borehole trajectory can be optimized according to different geological conditions and engineering requirements, dividing the borehole trajectory into 7 sections can more accurately describe the shape and characteristics of the borehole, so as to accurately determine the borehole trajectory with the shortest length. Based on the wellhead coordinates, target coordinates, and taking into account the influence of formation stability, collision avoidance, construction capacity, engineering conditions and other factors, the value ranges of the borehole trajectory parameters corresponding to multiple borehole trajectory sections can be determined, so as to accurately determine the borehole trajectory with the shortest length and the least friction, and to improve the precision and accuracy of drilling.

[0029] In an optional embodiment, the method further includes: a plurality of wellbore trajectory segments at least including: a straight well segment, a first deflection segment, a first stable deflection segment, a second deflection segment, a second stable deflection segment, a third deflection segment, and a horizontal segment, wherein the straight well segment is a straight line segment with a slope equal to a preset first threshold, the first deflection segment is a curved segment with a slope greater than or equal to a preset second threshold and less than a preset first threshold, the first stable deflection segment is a straight line segment with a slope greater than or equal to a preset third threshold and less than a preset second threshold, the second deflection segment is a curved segment with a slope greater than or equal to a preset fourth threshold and less than a preset third threshold, the second stable deflection segment is a straight line segment with a slope greater than or equal to a preset fifth threshold and less than a preset fourth threshold, and the third deflection segment is a curved segment with a slope greater than or equal to a preset fifth threshold and less than a preset fourth threshold. The straight line segment whose slope is greater than the preset sixth threshold and less than the preset fifth threshold, the horizontal segment is a straight line segment whose slope is equal to the preset sixth threshold, and the order from the preset first threshold to the preset sixth threshold from large to small is: the preset first threshold, the preset second threshold, the preset third threshold, the preset fourth threshold, the preset fifth threshold, and the preset sixth threshold; the wellbore trajectory parameters corresponding to the multiple wellbore trajectory segments include: the inclination point, the full-angle change rate, the terminal well inclination, and the terminal azimuth of the first inclination segment, the length of the first stable inclination segment, the full-angle change rate, the terminal well inclination, and the terminal azimuth of the second inclination segment, the length of the second stable inclination segment, the full-angle change rate, the terminal well inclination, the terminal azimuth of the third inclination segment, and the target accuracy.

[0030] Optional, still as Figure 2 As shown, the borehole trajectory can be segmented to obtain 7 borehole trajectory segments, including: the vertical well segment OA, where the slope of the vertical well segment is 1, which is a vertical well segment perpendicular to the reference plane (ground plane), the first deflection segment AB, i.e., the inclined arc 1, the first stable deflection segment BC, i.e., the inclined straight segment 1, the second deflection segment CD, i.e., the inclined arc 2, the second stable deflection segment DE, i.e., the inclined straight segment 2, the third deflection segment ET0, i.e., the inclined arc 3, and the horizontal segment T0T1, where the slope of the horizontal segment is 0, which is a horizontal well segment parallel to the reference plane. The borehole trajectory parameters include: H A is the inclination point of the first inclination section, K AB is the full angle change rate of the first deflection section, α B is the end well inclination angle of the first deflection section, is the end azimuth of the first deflection section, L BC is the length of the first stable slope section, K CD is the full angle change rate of the second deflection section, α D is the end well inclination angle of the second inclination section, is the end azimuth of the second deflection section, L DE is the length of the second stable slope section, is the full angle change rate of the third deflection section, is the end well inclination angle of the third inclination section, is the terminal azimuth of the third deflection section, and the target accuracy. Among them, the target accuracy is the sum of the squares of the difference between the given geological target point coordinates and the designed trajectory target point coordinates. The range of the target accuracy can be expressed as:

[0031]

[0032] Where, e is the target accuracy, The north-south displacement of the first target point T0, is the east-west displacement of point T0, is the vertical depth of point T0, i is the designed track entry point of any section among the vertical well section, the first deflection section, the first stable deflection section, the second deflection section, the second stable deflection section, and the third deflection section, ΔN i is the north-south displacement of the i-th target entry point, ΔE i is the east-west displacement of the i-th target point, ΔH i is the vertical depth of the i-th target point.

[0033] Step S106, establishing a wellbore trajectory model based on the value range of the wellbore trajectory parameters of the wellbore trajectory.

[0034] Optionally, by establishing a wellbore trajectory model through the value range of the wellbore trajectory parameters, the well trajectory can be automatically optimized, the workload of designers can be reduced, the friction of drilling construction can be reduced, and useless footage can be reduced, so as to reduce the difficulty and cost of drilling projects, and help the efficient, safe and beneficial development of unconventional resources such as shale oil and tight oil and gas.

[0035] In an optional embodiment, a wellbore trajectory model is established based on the value range of the wellbore trajectory parameters, including: determining a first objective function, wherein the first objective function is used to indicate that the trajectory length corresponding to the wellbore trajectory is the shortest; determining a second objective function, wherein the second objective function is used to indicate that the friction value corresponding to the wellbore trajectory during the drilling process is the smallest; and establishing the wellbore trajectory model based on the value ranges of the wellbore trajectory parameters corresponding to multiple wellbore trajectory segments, the first objective function, and the second objective function.

[0036] Optionally, the wellbore trajectory model is established by:

[0037]

[0038]

[0039]

[0040] Among them, minF1 is the first objective function, minF2 is the second objective function, and L OA , L AB , L BC , LCD , L DE , are the length parameters corresponding to the wellbore trajectory segments OA, AB, BC, CD, DE, ET0, and T0T1, respectively. i is the depth of the well at insertion point i, α i is the well inclination angle of insertion point i, is the azimuth of the insertion point i, is the minimum value of the deflection point of the first deflection section, is the maximum value of the deflection point of the first deflection section, is the minimum value of the full angle change rate of the first deflection section, is the maximum value of the full angle change rate of the first deflection section, is the minimum value of the end well inclination angle of the first deflection section, is the maximum value of the end well inclination angle of the first deflection section, is the minimum value of the terminal azimuth of the first deflection section, is the minimum value of the terminal azimuth of the first deflection section, is the minimum length of the first stable slope segment, is the maximum length of the first stable slope segment, is the minimum value of the full angle change rate of the second deflection section, is the maximum value of the full angle change rate of the second deflection section, is the minimum value of the end well inclination angle of the second deflection section, is the maximum value of the end well inclination angle of the second deflection section, is the minimum value of the terminal azimuth of the second deflection section, is the maximum value of the terminal azimuth of the second deflection section, is the minimum length of the second stable slope segment, is the maximum value of the length of the second stable slope segment, is the minimum value of the full angle change rate of the third deflection section, is the maximum value of the full angle change rate of the third deflection section, is the minimum value of the end well inclination angle of the third inclination section, is the maximum value of the end well inclination angle of the third inclination section, is the minimum value of the terminal azimuth of the third deflection section, It is the minimum value of the terminal azimuth of the third deflection section.

[0041] In an optional embodiment, determining the first objective function includes: obtaining length parameters corresponding to a plurality of wellbore trajectory segments respectively; and determining the first objective function based on the length parameters corresponding to the plurality of wellbore trajectory segments respectively.

[0042] Optionally, in order to minimize the length of the wellbore trajectory, obtain the length parameters L corresponding to the seven wellbore trajectory segments respectively. OA , L AB , L BC , L CD , L DE , Based on the length parameters corresponding to multiple wellbore trajectory segments, the first objective function is determined as follows: Among them, F1 is the total length of the wellbore trajectory.

[0043] In an optional embodiment, determining the second objective function includes: determining an axial force model and a torque model, wherein the axial force model is used to indicate changes in friction resistance experienced by the wellbore trajectory during linear drilling, and the torque model is used to indicate changes in friction resistance experienced by the wellbore trajectory during rotary drilling; determining the second objective function based on the axial force model and the torque model.

[0044] Optionally, the friction resistance of the wellbore trajectory during the drilling process includes axial force and torque, so the axial force model and the torque model can be determined by the following steps:

[0045] Step S1061, calculate the wellbore direction of the horizontal section by the following method and

[0046]

[0047]

[0048] Step S1062, calculate the coordinate increment parameters of any deflection section in the following manner:

[0049]

[0050]

[0051]

[0052]

[0053]

[0054] Among them, L is the length of the deflection section, ΔN is the north-south displacement increment of the deflection section, ΔE is the east-west displacement increment of the deflection section, ΔH is the vertical depth increment of the deflection section, α1 is the well inclination angle at the beginning of the deflection section, and α2 is the well inclination angle at the end of the deflection section. is the azimuth of the head end of the build section, is the azimuth of the end of the build section, γ is the dogleg angle, and λ M is the intermediate parameter of the calculation process.

[0055] Step S1063, calculating the coordinate increment parameters of any stable slope segment in the following manner:

[0056]

[0057]

[0058] ΔH=Lcosα2

[0059] Among them, L is the length of the stable inclination section, ΔN is the north-south displacement increment of the stable inclination section, ΔE is the east-west displacement increment of the BC section of the stable inclination section, ΔH is the vertical depth increment of the stable inclination section, and α2 is the well inclination angle at the end of the stable inclination section. is the azimuth at the end of the stable slope segment.

[0060] Step S1064: The measurement points of the curve in space can be calculated by interpolation using the space arc method. Figure 3 is a schematic diagram of an optional slope arc interpolation method according to an embodiment of the present invention, such as Figure 3 As shown in the figure, a measurement section 1-2 can be assumed to be a circular arc curve on a space inclined plane, and the position parameters of each point of the space arc can be solved. i And other given conditions, the depth length ΔL of the insertion point from the upper end point can be obtained i The well inclination angle α at the insertion point i can be calculated by the following formula: i , azimuth Vertical depth D i , north-south displacement N i , east-west displacement E i :

[0061]

[0062]

[0063]

[0064]

[0065] D i =D1+λ M (cosα1+cosα i )

[0066]

[0067]

[0068] Among them, α1 is the well inclination angle of measuring point 1, and α2 is the well inclination angle of measuring point 2. is the azimuth of measuring point 1, is the azimuth of measuring point 2, γ i is the dogleg length of insertion point i, △L is the distance between measuring point 2 and measuring point 1, △L i is the distance from insertion point i to measurement point 1, α i is the well inclination angle of insertion point i, is the azimuth of insertion point i, D i is the vertical depth of insertion point i, D1 is the vertical depth of measuring point 1, N i is the north-south displacement of insertion point i, N1 is the north-south displacement of measurement point 1, E i is the east-west displacement of insertion point i, and E1 is the east-west displacement of measuring point 1.

[0069] Step S1065, Figure 4 is a schematic diagram of an optional coordinate system according to an embodiment of the present invention, such as Figure 4 As shown, the N, E, and H axes point to the geographic north, geographic east, and center of gravity, respectively. t 、e n 、e b They are the tangent direction, principal normal direction and secondary normal direction of the wellbore axis respectively.

[0070] Step S1066, discretize the entire drill string into micro-element segments, and randomly select a micro-element segment with an arc length of ds on the wellbore trajectory curve. Figure 5 is a schematic diagram of an optional force analysis according to an embodiment of the present invention, such as Figure 5 The figure shows the force analysis of the microelement segment. The differential equation of the overall force of the drill string can be established in the following way:

[0071]

[0072]

[0073] in, In the figure, - stands for pulling, + stands for running, F is the axial tension on the drill string, M is T is the torque of the pipe string, s is the well depth, q is the weight of the drill string per unit length, α is the well inclination angle, EI is the bending stiffness of the drill string, and n t is the contact force between the drill string and the well wall, μ1 is the axial friction coefficient of the drill string, k b is the curvature of the wellbore axis.

[0074] Step S1067, calculate the wellbore curvature k by the following method b :

[0075]

[0076] in, is the azimuth.

[0077] Step S1068, calculate the contact force n on the drill string by the following method t :

[0078]

[0079] in:

[0080]

[0081]

[0082]

[0083] Among them, K s is the curvature of the helix, k n is the torsion of the borehole axis, F is the load, and F sin is the critical load under sinusoidal buckling.

[0084] In step S1069, the axial force model may be determined using a finite difference numerical solution method:

[0085]

[0086] The torque model can be determined using the finite difference numerical solution method:

[0087]

[0088] Among them, F i and F i+1 are the axial forces of the drill string near the ground and near the drill bit at the i-th and i+1-th sections, respectively, M Ti and M Ti+1 are the torques at both ends of the drill string at the i-th section and the i+1-th section respectively, and The curvature of the wellbore at both ends of the i-th section and the i+1-th section, q i is the linear weight of the drill string section i, EI i is the bending stiffness of the drill string section i, n ti is the contact force of the drill string wall in the i-th section, Δs i is the length of the unit end of the i-th section of the drill string wall, D bi is the outer diameter of the unit end of the wellbore wall of the i-th section of the drill string.

[0089] Step S10610, the axial force model and the torque model are modified in combination with the actual drilling conditions. Given the boundary of the differential equation and the friction coefficient value, the second objective function can be determined as:

[0090]

[0091] Among them, stiff is the calculation of the friction model, L i is the depth of the well at insertion point i, α i is the well inclination angle of insertion point i, is the azimuth of the insertion point i.

[0092] Step S108, based on the wellbore trajectory model, obtaining the wellbore trajectory information of the horizontal well, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory.

[0093] Optionally, the fast non-dominated sorting genetic NSGA-Ⅱ algorithm with elite strategy is usually used in the field of multi-objective optimization. It has the characteristics of fast convergence speed and good robustness. Therefore, the fast non-dominated sorting genetic NSGA-Ⅱ algorithm with elite strategy can be used to solve the wellbore trajectory model. In order to speed up the problem of long trajectory solution time, the map function in the scalable parallel operation library SCOOP in Python can be used to perform distributed calculations on the two objective functions in the wellbore trajectory model, which can greatly improve the performance and efficiency of the calculation and obtain the wellbore trajectory efficiently and accurately.

[0094] In an optional embodiment, based on the wellbore trajectory model, a horizontal wellbore trajectory determination result is obtained, including: taking the initial parent population as the current parent population, the initial child population as the current child population, and looping the following operations until a preset number of iterations is reached, wherein the individuals included in the initial parent population and the initial child population are both feasible solutions to the wellbore trajectory model and the number is N: merging the current parent population and the current child population to obtain a merged population, wherein the merged population includes 2N individuals; based on the wellbore trajectory model, determining the function values ​​corresponding to the 2N individuals respectively; based on the function values ​​corresponding to the 2N individuals respectively, performing non-dominated sorting on the 2N individuals respectively to obtain dominated sorting results corresponding to the 2N individuals respectively, wherein the non-dominated sorting is used to sort each of the 2N individuals An individual is divided into multiple levels; based on the dominance sorting results corresponding to the 2N individuals, the crowding degree calculation processing is performed on the 2N individuals to obtain the crowding degree results corresponding to the 2N individuals, wherein the crowding degree is used to indicate the density of the individual in the solution space; based on the dominance sorting results corresponding to the 2N individuals and the crowding degree results corresponding to the 2N individuals, a new parent population is determined; a first processing is performed on the new parent population to determine a new child population, wherein the first processing at least includes: selection processing, crossover processing, and mutation processing; the new parent population is used as the current parent population, and the new child population is used as the current child population; when a preset number of iterations is met, the horizontal wellbore trajectory determination result is determined based on the new child population, wherein the new child population is used to indicate the optimal solution set of the wellbore trajectory model.

[0095] Optionally, the fast non-dominated sorting genetic NSGA-Ⅱ algorithm with elite strategy can be programmed in Python, the environment is Python3.7, the compiler is anaconda, the genetic algorithm uses the "RI" encoding format, the population size can be 100, the genetic algorithm generation number (Generations) is 100, the crossover factor (Crossover Fraction) is 0.8, and the larger crossover probability can ensure that new individuals are generated faster and the algorithm's global search capability is stronger. Figure 6 is a flowchart of an optional wellbore trajectory model solution according to an embodiment of the present invention, such as Figure 6 As shown in the figure, the algorithm can be used to solve the wellbore trajectory model through the following steps:

[0096] Step S1080, using the initial parent population as the current parent population and the initial child population as the current child population, wherein the individuals included in the initial parent population and the initial child population are all feasible solutions of the wellbore trajectory model and the number of individuals is N;

[0097] Step S1081, merging the current parent population and the current child population to obtain a merged population, wherein the merged population includes 2N individuals;

[0098] Step S1082, determine whether to generate a new parent population, if yes, execute step S1087, if no, execute step S1083;

[0099] Step S1083, based on the wellbore trajectory model, determining the function values ​​corresponding to the 2N individuals respectively;

[0100] Step S1084, based on the function values ​​corresponding to the 2N individuals, respectively, non-dominated sorting processing is performed on the 2N individuals to obtain dominated sorting results corresponding to the 2N individuals, wherein the non-dominated sorting processing is used to divide each of the 2N individuals into multiple levels;

[0101] Step S1085, based on the dominance ranking results corresponding to the 2N individuals, respectively, the crowding degree calculation processing is performed on the 2N individuals to obtain the crowding degree results corresponding to the 2N individuals, wherein the crowding degree is used to indicate the density of the individuals in the solution space;

[0102] Step S1086, determining a new parent population based on the dominance ranking results corresponding to the 2N individuals and the crowding degree results corresponding to the 2N individuals;

[0103] Step S1087, performing a first process on the new parent population to determine a new child population, wherein the first process at least includes: selection process, crossover process, and mutation process;

[0104] Step S1088, determine whether the preset number of iterations is met, if yes, execute step S1089, if no, add 1 to the current number of iterations, use the new parent population as the current parent population, use the new child population as the current child population, and execute step S1081;

[0105] Step S1089, determining the horizontal wellbore trajectory determination result based on the new offspring population, wherein the new offspring population is used to indicate the optimal solution set of the wellbore trajectory model.

[0106] In an optional embodiment, the initial parent population is used as the current parent population, and the initial child population is used as the current child population, including: randomly generating an initial parent population, wherein the initial parent population includes N individuals; determining function values ​​corresponding to the N individuals based on the wellbore trajectory model; performing non-dominated sorting processing on the N individuals based on the function values ​​corresponding to the N individuals to obtain non-dominated sorting results corresponding to the N individuals; performing a first processing on the N individuals based on the non-dominated sorting results corresponding to the N individuals to determine the initial child population.

[0107] Optional, still as Figure 6 As shown, the initial parent population can be used as the current parent population and the initial child population as the current child population through the following steps:

[0108] Step S10801, randomly generating an initial parent population, wherein the initial parent population includes N individuals;

[0109] Step S10802, determine whether to generate an initial offspring population, if yes, execute step S10805, if no, execute step S10803;

[0110] Step S10803, based on the wellbore trajectory model, determining the function values ​​corresponding to the N individuals respectively;

[0111] Step S10804, based on the function values ​​corresponding to the N individuals, respectively perform non-dominated sorting processing on the N individuals to obtain non-dominated sorting results corresponding to the N individuals;

[0112] Step S10805, based on the non-dominated sorting results corresponding to the N individuals, perform a first process on the N individuals to determine an initial offspring population;

[0113] Step S10805, the current number of iterations is set to 2.

[0114] Through the above steps S102 to S108, the purpose of establishing a wellbore trajectory model based on the value range of the corresponding parameters of the wellbore trajectory is achieved, so as to efficiently obtain the wellbore trajectory determination result of the horizontal well, thereby achieving the technical effect of improving the wellbore trajectory determination efficiency of the horizontal well and reducing the wellbore trajectory determination cost of the horizontal well, and further solving the technical problems of low wellbore trajectory determination efficiency and high cost existing in the wellbore trajectory determination method of the horizontal well in the related technology.

[0115] Based on the above embodiments and optional embodiments, the present invention proposes a method for determining a wellbore trajectory of an optional horizontal well, the method comprising:

[0116] Step S1, obtaining the wellhead coordinates and target point coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target point coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir.

[0117] Step S2, segmenting the wellbore trajectory to obtain a plurality of wellbore trajectory segments.

[0118] Step S3, based on the wellhead coordinates and the target point coordinates, determine the value ranges of the wellbore trajectory parameters corresponding to the multiple wellbore trajectory segments, wherein the wellbore trajectory parameters include: the inclination point, full-angle change rate, terminal well inclination angle, and terminal azimuth of the first inclination segment, the length of the first stable inclination segment, the full-angle change rate, terminal well inclination angle, and terminal azimuth of the second inclination segment, the full-angle change rate, terminal well inclination angle, and terminal azimuth of the second stable inclination segment and the third inclination segment, as well as the target accuracy.

[0119] Step S4, obtaining length parameters corresponding to a plurality of wellbore trajectory segments respectively.

[0120] Step S5, determining a first objective function based on length parameters corresponding to a plurality of wellbore trajectory segments, wherein the first objective function is used to indicate that the trajectory length corresponding to the wellbore trajectory is the shortest.

[0121] Step S6, determining an axial force model and a torque model, wherein the axial force model is used to indicate the change in friction resistance of the wellbore trajectory during linear drilling, and the torque model is used to indicate the change in friction resistance of the wellbore trajectory during rotary drilling.

[0122] Step S7, based on the axial force model and the torque model, determine the second objective function, wherein the second objective function is used to indicate that the friction value corresponding to the wellbore trajectory during the drilling process is minimum.

[0123] Step S8, establishing a wellbore trajectory model based on the value ranges of the wellbore trajectory parameters corresponding to the plurality of wellbore trajectory segments, the first objective function, and the second objective function.

[0124] Step S9, using a fast non-dominated sorting genetic algorithm NSGA-II with an elite strategy, using Python programming and the map function in the scalable parallel operation library SCOOP in Python to solve the wellbore trajectory model.

[0125] Step S10, obtaining a horizontal wellbore trajectory determination result, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory.

[0126] Step S11, based on the horizontal wellbore trajectory determination result, draw a three-dimensional wellbore trajectory diagram, a vertical profile diagram, a horizontal well projection diagram, key point locations and engineering parameters.

[0127] In this embodiment, a device for determining the borehole trajectory of a horizontal well is also provided, and the device is used to implement the above-mentioned embodiments and preferred implementation modes, and the descriptions that have been made will not be repeated. As used below, the terms "module" and "device" can implement a combination of software and / or hardware for a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0128] According to an embodiment of the present invention, there is also provided an embodiment of a device for implementing the above-mentioned method for determining a wellbore trajectory of a horizontal well. Figure 7 is a schematic diagram of the structure of a wellbore trajectory determination device for a horizontal well according to an embodiment of the present invention. Figure 7 As shown, the wellbore trajectory determination device of the horizontal well includes: a first acquisition module 702, a first determination module 704, a first establishment module 706, and a second determination module 708, wherein:

[0129] The first acquisition module 702 is used to acquire the wellhead coordinates and target point coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target point coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir;

[0130] A first determination module 704, connected to the first acquisition module 702, is used to determine a value range of a wellbore trajectory parameter of the wellbore trajectory based on the wellhead coordinates and the target coordinates;

[0131] A first establishing module 706, connected to the first determining module 704, is used to establish a wellbore trajectory model based on the value range of the wellbore trajectory parameter;

[0132] The second determination module 708 is connected to the first establishment module 706, and is used to obtain the wellbore trajectory information of the horizontal well based on the wellbore trajectory model, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory.

[0133] By setting a first acquisition module 702, it is used to obtain the wellhead coordinates and target coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir; a first determination module 704 is used to determine the value range of the wellbore trajectory parameters of the wellbore trajectory based on the wellhead coordinates and the target coordinates; a first establishment module 706 is used to establish a wellbore trajectory model based on the value range of the wellbore trajectory parameters; a second determination module 708 is used to obtain the wellbore trajectory information of the horizontal well based on the wellbore trajectory model, wherein the wellbore trajectory information includes the wellbore trajectory, and the track length and friction value corresponding to the wellbore trajectory. The purpose of establishing a wellbore trajectory model based on the value range of the parameters corresponding to the wellbore trajectory, thereby efficiently obtaining the wellbore trajectory determination result of the horizontal well, thereby achieving the technical effect of improving the wellbore trajectory determination efficiency of the horizontal well and reducing the wellbore trajectory determination cost of the horizontal well, thereby solving the technical problems of low wellbore trajectory determination efficiency and high cost in the wellbore trajectory determination method of the horizontal well in the related art.

[0134] It should be noted that the above modules can be implemented by software or hardware. For example, for the latter, it can be implemented in the following ways: the above modules can be located in the same processor; or the above modules can be located in different processors in any combination.

[0135] It should be noted that the first acquisition module 702, the first determination module 704, the first establishment module 706, and the second determination module 708 correspond to steps S102 to S108 in the embodiment, and the examples and application scenarios implemented by the modules and the corresponding steps are the same, but are not limited to the contents disclosed in the above embodiments. It should be noted that the modules as part of the device can be run in a computer terminal.

[0136] It should be noted that the optional or preferred implementation of this embodiment can refer to the relevant description in the embodiment, which will not be repeated here.

[0137] The above-mentioned horizontal wellbore trajectory determination device may also include a processor and a memory. The above-mentioned first acquisition module 702, first determination module 704, first establishment module 706, second determination module 708, etc. are all stored in the memory as program modules, and the processor executes the above-mentioned program modules stored in the memory to realize corresponding functions.

[0138] The processor includes a kernel, which retrieves the corresponding program module from the memory. The kernel may be one or more. The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. The memory includes at least one memory chip.

[0139] According to an embodiment of the present application, an embodiment of a non-volatile storage medium is also provided. Optionally, in this embodiment, the non-volatile storage medium includes a stored program, wherein when the program is executed, the device where the non-volatile storage medium is located is controlled to execute any of the above-mentioned methods for determining the wellbore trajectory of a horizontal well.

[0140] Optionally, in this embodiment, the non-volatile storage medium may be located in any computer terminal in a computer terminal group in a computer network, or in any mobile terminal in a mobile terminal group, and the non-volatile storage medium includes a stored program.

[0141] Optionally, when the program is running, the device where the non-volatile storage medium is located is controlled to perform the following functions: obtain the wellhead coordinates and target point coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target point coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir; determine the value range of the wellbore trajectory parameters of the wellbore trajectory based on the wellhead coordinates and the target point coordinates; establish a wellbore trajectory model based on the value range of the wellbore trajectory parameters; obtain the wellbore trajectory information of the horizontal well based on the wellbore trajectory model, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory.

[0142] According to an embodiment of the present application, an embodiment of a processor is also provided. Optionally, in this embodiment, the processor is used to run a program, wherein when the program is run, any of the above-mentioned methods for determining a wellbore trajectory of a horizontal well is executed.

[0143] According to an embodiment of the present application, an embodiment of a computer program product is also provided. When executed on a data processing device, the program is suitable for executing the steps of initializing the wellbore trajectory determination method of any one of the horizontal wells described above.

[0144] Optionally, the above-mentioned computer program product, when executed on a data processing device, is suitable for executing an initialization program having the following method steps: obtaining the wellhead coordinates and target point coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target point coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir; determining the value range of the wellbore trajectory parameters of the wellbore trajectory based on the wellhead coordinates and the target point coordinates; establishing a wellbore trajectory model based on the value range of the wellbore trajectory parameters; and obtaining the wellbore trajectory information of the horizontal well based on the wellbore trajectory model, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory.

[0145] Optionally, the computer program product is also suitable for executing a program initialized with the following method steps: segmenting the wellbore trajectory to obtain multiple wellbore trajectory segments; based on the wellhead coordinates and the target point coordinates, determining the value ranges of the wellbore trajectory parameters corresponding to the multiple wellbore trajectory segments.

[0146] Optionally, the above-mentioned computer program product is also suitable for executing an initialized program having the following method steps: determining a first objective function, wherein the first objective function is used to indicate that the trajectory length corresponding to the wellbore trajectory is the shortest; determining a second objective function, wherein the second objective function is used to indicate that the friction value corresponding to the wellbore trajectory during the drilling process is the smallest; and establishing a wellbore trajectory model based on the value ranges of wellbore trajectory parameters corresponding to multiple wellbore trajectory segments, the first objective function, and the second objective function.

[0147] Optionally, the computer program product is further suitable for executing a program initialized with the following method steps: obtaining length parameters corresponding to a plurality of wellbore trajectory segments respectively; and determining a first objective function based on the length parameters corresponding to the plurality of wellbore trajectory segments respectively.

[0148] Optionally, the above-mentioned computer program product is also suitable for executing an initialization program having the following method steps: determining an axial force model and a torque model, wherein the axial force model is used to indicate the change in friction resistance experienced by the wellbore trajectory during linear drilling, and the torque model is used to indicate the change in friction resistance experienced by the wellbore trajectory during rotary drilling; based on the axial force model and the torque model, determining the second objective function.

[0149] Optionally, the above-mentioned computer program product is also suitable for executing a program that is initialized with the following method steps: using the initial parent population as the current parent population, the initial child population as the current child population, and looping to perform the following operations until a preset number of iterations is reached, wherein the individuals included in the initial parent population and the initial child population are all feasible solutions to the wellbore trajectory model and the number is N: merging the current parent population and the current child population to obtain a merged population, wherein the merged population includes 2N individuals; based on the wellbore trajectory model, determining the function values ​​corresponding to the 2N individuals respectively; based on the function values ​​corresponding to the 2N individuals respectively, performing non-dominated sorting processing on the 2N individuals respectively to obtain dominated sorting results corresponding to the 2N individuals respectively, wherein the non-dominated sorting processing is used to sort each of the 2N individuals. The invention discloses a method for dividing a population into multiple levels; based on the dominance sorting results corresponding to the 2N individuals, respectively, performing crowding calculation processing on the 2N individuals to obtain crowding results corresponding to the 2N individuals, wherein the crowding is used to indicate the density of the individual in the solution space; based on the dominance sorting results corresponding to the 2N individuals and the crowding results corresponding to the 2N individuals, respectively, determining a new parent population; performing a first processing on the new parent population to determine a new child population, wherein the first processing at least includes: selection processing, crossover processing, and mutation processing; using the new parent population as the current parent population and the new child population as the current child population; and determining the horizontal wellbore trajectory determination result based on the new child population when a preset number of iterations is met, wherein the new child population is used to indicate the optimal solution set of the wellbore trajectory model.

[0150] Optionally, the above-mentioned computer program product is also suitable for executing a program initialized with the following method steps: randomly generating an initial parent population, wherein the initial parent population includes N individuals; based on the wellbore trajectory model, determining the function values ​​corresponding to the N individuals respectively; based on the function values ​​corresponding to the N individuals respectively, performing non-dominated sorting processing on the N individuals respectively to obtain non-dominated sorting results corresponding to the N individuals respectively; based on the non-dominated sorting results corresponding to the N individuals respectively, performing a first processing on the N individuals to determine the initial offspring population.

[0151] An embodiment of the present invention provides an electronic device, such as Figure 8As shown, the electronic device 10 includes a processor, a memory, and a program stored in the memory and executable on the processor. When the processor executes the program, the following steps are implemented: obtaining the wellhead coordinates and target coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir; determining the value range of the wellbore trajectory parameters of the wellbore trajectory based on the wellhead coordinates and the target coordinates; establishing a wellbore trajectory model based on the value range of the wellbore trajectory parameters; and obtaining the wellbore trajectory information of the horizontal well based on the wellbore trajectory model, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory.

[0152] In the above embodiments of the present invention, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0153] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. Among them, the device embodiments described above are only schematic. For example, the division of the above modules can be a logical function division. There may be other division methods in actual implementation, such as multiple modules 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 through some interfaces, modules or indirect coupling or communication connection of modules, which can be electrical or other forms.

[0154] The modules described above as separate components may or may not be physically separated, and the components shown as modules may or may not be physical modules, that is, they may be located in one place or distributed on multiple modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the present embodiment.

[0155] In addition, each functional module in each embodiment of the present invention may be integrated into one processing module, or each module may exist physically separately, or two or more modules may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of software functional modules.

[0156] If the above-mentioned integrated module is implemented in the form of a software function module and sold or used as an independent product, it can be stored in a computer-readable non-volatile storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a non-volatile storage medium, including a number of instructions to enable a computer device (which can be a personal computer, server or network device, etc.) to perform all or part of the steps of the various embodiments of the present invention. The aforementioned non-volatile storage medium includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes.

[0157] The above are only preferred embodiments of the present invention. It should be pointed out that, for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A method for determining a wellbore trajectory of a horizontal well, characterized in that: include: Obtaining the wellhead coordinates and target point coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target point coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir; Determining a value range of a wellbore trajectory parameter of the wellbore trajectory based on the wellhead coordinates and the target point coordinates; Based on the value range of the wellbore trajectory parameter, establishing a wellbore trajectory model; Based on the wellbore trajectory model, wellbore trajectory information of the horizontal well is obtained, wherein the wellbore trajectory information includes the wellbore trajectory, and a trajectory length and a friction resistance value corresponding to the wellbore trajectory.

2. The method according to claim 1, characterized in that The determining of the value range of the wellbore trajectory parameter of the wellbore trajectory based on the wellhead coordinates and the target point coordinates includes: Segmenting the wellbore trajectory to obtain a plurality of wellbore trajectory segments; Based on the wellhead coordinates and the target point coordinates, the value ranges of the wellbore trajectory parameters respectively corresponding to the multiple wellbore trajectory segments are determined.

3. The method according to claim 2, characterized in that The step of establishing a wellbore trajectory model based on the value range of the wellbore trajectory parameter comprises: Determining a first objective function, wherein the first objective function is used to indicate that the trajectory length corresponding to the wellbore trajectory is the shortest; Determining a second objective function, wherein the second objective function is used to indicate that the friction value corresponding to the wellbore trajectory during the drilling process is minimal; The wellbore trajectory model is established based on the value ranges of the wellbore trajectory parameters respectively corresponding to the multiple wellbore trajectory segments, the first objective function, and the second objective function.

4. The method according to claim 3, characterized in that: The determining of the first objective function comprises: Obtaining length parameters corresponding to the plurality of wellbore trajectory segments respectively; The first objective function is determined based on the length parameters respectively corresponding to the plurality of wellbore trajectory segments.

5. The method according to claim 3, characterized in that: The determining of the second objective function comprises: Determine an axial force model and a torque model, wherein the axial force model is used to indicate the change of friction resistance of the wellbore trajectory during linear drilling, and the torque model is used to indicate the change of friction resistance of the wellbore trajectory during rotary drilling; The second objective function is determined based on the axial force model and the torque model.

6. The method according to any one of claims 2 to 5, characterized in that The method further comprises: The multiple wellbore trajectory segments at least include: a straight well segment, a first deflection segment, a first stable deflection segment, a second deflection segment, a second stable deflection segment, a third deflection segment, and a horizontal segment, wherein the straight well segment is a straight line segment with a slope equal to a preset first threshold value, the first deflection segment is a curve segment with a slope greater than or equal to a preset second threshold value and less than the preset first threshold value, the first stable deflection segment is a straight line segment with a slope greater than or equal to a preset third threshold value and less than the preset second threshold value, the second deflection segment is a curve segment with a slope greater than or equal to a preset fourth threshold value and less than the preset third threshold value, the second stable deflection segment is a straight line segment with a slope greater than or equal to a preset fifth threshold value and less than the preset fourth threshold value, the third deflection segment is a straight line segment with a slope greater than a preset sixth threshold value and less than the preset fifth threshold value, and the horizontal segment is a straight line segment with a slope equal to the preset sixth threshold value, and the order from the preset first threshold value to the preset sixth threshold value from large to small is: the preset first threshold value, the preset second threshold value, the preset third threshold value, the preset fourth threshold value, the preset fifth threshold value, and the preset sixth threshold value; The wellbore trajectory parameters corresponding to the multiple wellbore trajectory segments include: the inclination point, full-angle change rate, terminal well inclination angle, and terminal azimuth of the first inclination segment, the length of the first stable inclination segment, the full-angle change rate, terminal well inclination angle, and terminal azimuth of the second inclination segment, the length of the second stable inclination segment, the full-angle change rate, terminal well inclination angle, terminal azimuth, and target accuracy of the third inclination segment.

7. The method according to claim 1, characterized in that The step of obtaining a horizontal wellbore trajectory determination result based on the wellbore trajectory model includes: The initial parent population is used as the current parent population, the initial child population is used as the current child population, and the following operations are performed cyclically until a preset number of iterations is reached, wherein the individuals included in the initial parent population and the initial child population are both feasible solutions of the wellbore trajectory model and the number of individuals is N: Merging the current parent population and the current child population to obtain a merged population, wherein the merged population includes 2N individuals; Based on the wellbore trajectory model, determining function values ​​corresponding to the 2N individuals respectively; Based on the function values ​​corresponding to the 2N individuals, respectively, non-dominated sorting processing is performed on the 2N individuals to obtain dominated sorting results corresponding to the 2N individuals, wherein the non-dominated sorting processing is used to divide each of the 2N individuals into multiple levels; Based on the dominance sorting results corresponding to the 2N individuals, respectively, the 2N individuals are respectively subjected to crowding calculation processing to obtain crowding results corresponding to the 2N individuals, wherein the crowding is used to indicate the density of the individuals in the solution space; Determine a new parent population based on the dominance sorting results corresponding to the 2N individuals and the crowding degree results corresponding to the 2N individuals; Performing a first process on the new parent population to determine a new child population, wherein the first process at least includes: selection process, crossover process, and mutation process; Using the new parent population as the current parent population and the new child population as the current child population; When the preset number of iterations is met, the horizontal wellbore trajectory determination result is determined based on the new offspring population, wherein the new offspring population is used to indicate an optimal solution set of the wellbore trajectory model.

8. The method according to claim 7, characterized in that The method of taking the initial parent population as the current parent population and the initial child population as the current child population includes: Randomly generate an initial parent population, wherein the initial parent population includes N individuals; Based on the wellbore trajectory model, determining function values ​​corresponding to the N individuals respectively; Based on the function values ​​corresponding to the N individuals, respectively, performing the non-dominated sorting process on the N individuals, to obtain the non-dominated sorting results corresponding to the N individuals; Based on the non-dominated sorting results respectively corresponding to the N individuals, the first processing is performed on the N individuals to determine the initial offspring population.

9. A device for determining a wellbore trajectory of a horizontal well, characterized in that: include: A first acquisition module is used to acquire the wellhead coordinates and target point coordinates corresponding to the wellbore trajectory of the horizontal well, wherein the target point coordinates are used to indicate the coordinates corresponding to the drilling point in the reservoir; A first determination module, configured to determine a value range of a wellbore trajectory parameter of the wellbore trajectory based on the wellhead coordinates and the target point coordinates; A first establishing module, used for establishing a wellbore trajectory model based on the value range of the wellbore trajectory parameter; The second determination module is used to obtain the wellbore trajectory information of the horizontal well based on the wellbore trajectory model, wherein the wellbore trajectory information includes the wellbore trajectory, and the trajectory length and friction value corresponding to the wellbore trajectory.

10. An electronic device, characterized in that: It includes one or more processors and a memory, wherein the memory is used to store one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the method for determining the wellbore trajectory of a horizontal well as described in any one of claims 1 to 7.

Images