阅读说明：本技术 一种丛式井平台位置优选方法、系统及存储介质 (Cluster well platform position optimization method, system and storage medium ) 是由 袁俊亮 李中 范白涛 幸雪松 谢仁军 何英明 孙翀 吴怡 于 2021-07-21 设计创作，主要内容包括：本发明公开了一种丛式井平台位置优选方法、系统及存储介质。所述方法包括如下步骤：S1、确定目标油气田的井数n、靶点坐标、靶点垂深以及受限区的边界范围；S2、利用粒子群算法对平台位置进行T轮迭代,T为最大迭代次数,依次记录每轮迭代的整个粒子群的全局最优平台位置g和最小总进尺gbest,T轮迭代后的g和gbest就是考虑受限区影响的钻井总进尺最小的平台位置及相应总进尺数。本发明能够考虑受限区对海上油气田开发过程中平台位置优选的影响,突破了现有方法未考虑受限区的缺陷,为海上油气田开发降低钻井进尺和成本提供了保障。(The invention discloses a cluster well platform position optimization method, a cluster well platform position optimization system and a storage medium. The method comprises the following steps: s1, determining the well number n, the target point coordinates, the target point vertical depth and the boundary range of the limited area of the target oil-gas field; and S2, performing T-round iteration on the platform position by using a particle swarm algorithm, wherein T is the maximum iteration number, and sequentially recording the global optimal platform position g and the minimum total footage gbest of the whole particle swarm of each round of iteration, wherein the g and the gbest after the T-round iteration are the platform position with the minimum total footage of the drilling well and the corresponding total footage number which are influenced by the restricted area. The method can consider the influence of the restricted area on the optimization of the platform position in the offshore oil and gas field development process, overcomes the defect that the restricted area is not considered in the conventional method, and provides guarantee for reducing the drilling footage and the cost in the offshore oil and gas field development.)

1. A cluster well platform position optimization method comprises the following steps:

s1, determining the well number n, the target point coordinates, the target point vertical depth and the boundary range of the limited area of the target oil-gas field;

and S2, performing T-round iteration on the platform position by using a particle swarm algorithm, wherein T represents the maximum iteration times, and sequentially recording the global optimal platform position g and the minimum total footage gbest of the whole particle swarm of each round of iteration, and the global optimal platform position g and the final global minimum total footage gbest after the T-round iteration are the determined global optimal platform position and the global minimum total footage.

2. The preferred method of claim 1, wherein step S2 includes the steps of:

s21, presetting the number N of particles, the maximum iteration number T, learning factors C1 and C2, an inertia weight maximum value Wmax and an inertia weight minimum value Wmin, and a speed maximum value Vmax and a speed minimum value Vmin;

s22, initializing a position x (i, j) of each particle, wherein i represents the ith particle, j is 1 or 2, x (i,1) represents the east-west coordinates of the particle, x (i,2) represents the north-south coordinates of the particle, and in the initial stage, the positions of the particles are randomly distributed in the oil field expansion range and outside the restricted area;

s23, initializing the speeds v (i, j) of the N particles, wherein the initial speed is between Vmin and Vmax, the speed is positive to increase x (i, j), and the speed is negative to decrease x (i, j);

s24, calculating an initial optimal position p (i, j) and an initial minimum total footage pbest (i) of each particle, wherein the initial optimal position p (i, j) is the initial position x (i, j) in the step S22, the initial minimum total footage pbest (i) is the total drilling footage corresponding to each position, and N positions correspond to N total drilling footages;

s25, calculating the global optimal platform position g and the global minimum total footage of the whole particle swarm in the initial stagegbest：

Traversing from the 1 st particle to the Nth particle, selecting the minimum total footage pbest (i) of all the particles as an initial global minimum total footage gbest, and taking the corresponding platform position p (i, j) as an initial global optimal position g;

s26, carrying out a first iteration, updating the positions x (i, j) of N particles to obtain new x (i, j), recalculating the total footage pbest (i) of each particle, selecting the minimum value as the global minimum total footage gbest after the first iteration, and newly using the corresponding particle position x (i, j) as the global optimal platform position g after the first iteration;

and S27, performing 2 nd to T th iterations according to the iteration method of the step S26, and sequentially recording the global optimal platform position g and the global minimum total footage gbest of each iteration, wherein the global optimal platform position g and the final global minimum total footage gbest after the T iterations are the determined global optimal platform position and the determined global minimum total footage.

3. The preferred method of claim 2, wherein: in step S21, the number N of particles is 2 × m, the maximum number of iterations T is 10 × m, the learning factor C1 is C2 is 1.5, the maximum recommended inertia weight Wmax is 0.9, the minimum Wmin is 0.4, the maximum speed Vmax is 1, and the minimum Vmin is-1;

where m represents the area of the oilfield region, km²。

4. A preferred method according to claim 2 or 3, characterized in that: in step S24, the initial minimum total footage pbest (i) of the ith particle is obtained as follows:

wherein n represents the number of wells, i is 1 to n; dep(s) represents the footage of the s-th well, including straight well sections, deflecting sections, steady deflecting sections, and the like.

5. The preferred method according to any one of claims 2-4, characterized in that: in step S26, the method for updating the position x (i, j) includes the following steps:

first, calculating a dynamic inertia weight w:

w＝Wmax-(Wmax-Wmin)×k/T

wherein k represents the kth iteration;

the velocity v (i, j) of the particle is updated as follows:

v (i, j) new w × v (i, j) + C1 × rand x [ p (i, j) -x (i, j) ] + C2 × rand x [ g-x (i, j) ]

Wherein, C1 and C2 are learning factors, v (i, j) ═ rand x (Vmax-Vmin) + Vmin;

if the obtained v (i, j) is not between Vmin and Vmax, the regularization processing is carried out:

the position x (i, j) of the particle is then updated, as follows:

x (i, j) new ═ x (i, j) + vx (i, j)

If the obtained x (i, j) is not within the expanded range and is outside the restricted area, the regularization processing is carried out:

wherein Xmin, Xman, Ymin and Ymax are the boundaries of the west, east, south and north coordinates of the expanded range respectively.

6. The preferred method according to any one of claims 2-5, characterized in that: in step S27, after the particle swarm updates the position, if the global minimum total footage gbest is smaller than the global minimum total footage gbest of the previous round, the corresponding platform position g is recorded, and the final global optimal platform position g and the final global minimum total footage gbest are obtained through T-round iteration.

7. A system for optimizing cluster well platform locations comprising a processor and a memory storing a computer program; the processor is configured to execute the computer program to implement the cluster well platform location-preferred method of any of claims 1-6.

8. A computer storage medium, characterized in that: the computer storage medium having stored thereon a computer program that, when executed by a processor, implements the cluster well platform location-preferred method of any of claims 1-6.

Technical Field

The invention relates to a cluster well platform position optimization method, a cluster well platform position optimization system and a storage medium, and belongs to the technical field of petroleum engineering artificial intelligence application.

Background

In the process of offshore oil and gas field development, sea contradiction is prominent to face a big problem, Bohai oil field is taken as an example, the Bohai oil field is influenced by various factors such as national defense, navigation channel, fishery, environmental protection and the like, the platform position location is greatly limited, the platform position cannot be deployed on the sea surface of a 'limited area' such as the navigation channel, military area and the like, and the position can only be optimized in the sea area outside the limited area. Therefore, how to preferably select the platform position with the minimum drilling footage outside the restricted area is an urgent problem to be solved.

Depending on the starting point, a preferred method of cluster well platform location includes: a minimum total footage method, a minimum total horizontal displacement method, a difficulty weight method and the like. Firstly, the minimum total footage method is preferably used for obtaining the position with the minimum total footage sum of all target wells, but the platform position is mostly positioned in the center of the target area, the drilling difficulty is high, and the target area is probably positioned in a limited area. Secondly, the minimum total horizontal displacement method preferably selects the position with the minimum sum of the horizontal displacements from each target point on the plane, but can only represent the minimum sum of the plane distances of the target points, and cannot display the influence of the depth. And thirdly, the difficulty weighting method balances the drilling difficulty on the basis of considering the total footage, adopts a mode of subjectively selecting the position of the platform by hand to select, and cannot ensure the minimum total footage. The method does not consider the constraint influence of the restricted area, and the selected platform position is probably positioned in the restricted area, so that the actual drilling operation cannot be implemented.

Disclosure of Invention

The invention aims to provide a cluster well platform position optimization method, which is suitable for solving the problem of platform position optimization with the minimum drilling total footage in the process of offshore oil and gas field development under the constraint of a limited area.

The invention provides a cluster well platform position optimization method, which comprises the following steps:

s1, determining the well number n, the target point coordinates, the target point vertical depth and the boundary range of the limited area of the target oil-gas field;

and S2, performing T-round iteration on the platform position by using a particle swarm algorithm, wherein T represents the maximum iteration times, and sequentially recording the global optimal platform position g and the minimum total footage gbest of the whole particle swarm of each round of iteration, and the global optimal platform position g and the final global minimum total footage gbest after the T-round iteration are the determined global optimal platform position and the global minimum total footage.

Specifically, step S2 includes the steps of:

s21, presetting the number N of particles (representing N platform positions), the maximum iteration number T, learning factors C1 and C2, the maximum value Wmax of inertia weight, the minimum value Wmin of inertia weight, the maximum value Vmax of speed and the minimum value Vmin of speed;

s22, initializing a position x (i, j) of each particle, wherein i represents the ith particle, j is 1 or 2, x (i,1) represents the east-west coordinates of the particle, x (i,2) represents the north-south coordinates of the particle, and in the initial stage, the positions of the particles are randomly distributed in the oil field expansion range and outside the restricted area;

s23, initializing the speeds v (i, j) of the N particles, wherein the initial speed is between Vmin and Vmax, the speed is positive to increase x (i, j), and the speed is negative to decrease x (i, j);

s24, calculating an initial optimal position p (i, j) and an initial minimum total footage pbest (i) of each particle, wherein the initial optimal position p (i, j) is the initial position x (i, j) in the step S22, the initial minimum total footage pbest (i) is the total drilling footage corresponding to each position, and N positions correspond to N total drilling footages;

s25, calculating the global optimal platform position g and the global minimum total footage of the whole particle swarm in the initial stagegbest：

Traversing from the 1 st particle to the Nth particle, selecting the minimum total footage pbest (i) of all the particles as an initial global minimum total footage gbest, and taking the corresponding platform position p (i, j) as an initial global optimal position g;

s26, carrying out a first iteration, updating the positions x (i, j) of N particles to obtain new x (i, j), recalculating the total footage pbest (i) of each particle, selecting the minimum value as the global minimum total footage gbest after the first iteration, and newly using the corresponding particle position x (i, j) as the global optimal platform position g after the first iteration;

and S27, performing 2 nd to T th iterations according to the iteration method of the step S26, and sequentially recording the global optimal platform position g and the global minimum total footage gbest of each iteration, wherein the global optimal platform position g and the final global minimum total footage gbest after the T iterations are the determined global optimal platform position and the determined global minimum total footage.

Specifically, in step S21, the number N of particles is 2 × m, the maximum number of iterations T is 10 × m, the learning factor C1 is C2 is 1.5, the maximum recommended inertia weight Wmax is 0.9, the minimum value Wmin is 0.4, the maximum speed Vmax is 1, and the minimum value Vmin is-1;

where m represents the area of the oilfield region, km²。

Specifically, in step S24, the initial minimum total footage pbest (i) of the ith particle is obtained as follows:

wherein n represents the number of wells, i is 1 to n; dep(s) represents the footage of the s-th well, including straight well sections, deflecting sections, steady deflecting sections, and the like.

Specifically, in step S22, the initial position of each particle is randomly distributed within a range of 10 times the field area enlargement (referred to as "enlargement range"), and the number of times can be adjusted according to actual conditions.

Specifically, in step S26, the method for updating the position x (i, j) includes the following steps:

first, calculating a dynamic inertia weight w:

w＝Wmax-(Wmax-Wmin)×k/T

wherein k represents the kth iteration;

the velocity v (i, j) of the particle is updated as follows:

v (i, j) new w × v (i, j) + C1 × rand x [ p (i, j) -x (i, j) ] + C2 × rand x [ g-x (i, j) ]

Wherein, C1 and C2 are learning factors, v (i, j) ═ rand x (Vmax-Vmin) + Vmin;

if the obtained v (i, j) is not between Vmin and Vmax, the regularization processing is carried out:

the position x (i, j) of the particle is then updated, as follows:

x (i, j) new ═ x (i, j) + vx (i, j)

If the obtained x (i, j) is not within the expanded range and is outside the restricted area, the regularization processing is carried out:

wherein Xmin, Xman, Ymin and Ymax are the boundaries of the west, east, south and north coordinates of the expanded range respectively.

Specifically, in step S27, after the particle swarm updates the position, if the global minimum total footage gbest is smaller than the global minimum total footage gbest of the previous round, the corresponding platform position g is recorded, and the final global optimal platform position g and the final global minimum total footage gbest are obtained through T-round iteration.

The final global optimal platform position g determined by the method represents that after the influence of the restricted area is considered, the platform is selected at the position, and the total drilling footage is minimum.

The invention also provides a system for optimizing cluster well platform locations, comprising a processor and a memory storing a computer program; the processor is configured to execute the computer program to implement the cluster well platform location preferred method of the present invention.

The present invention still further provides a computer storage medium having a computer program stored thereon which, when executed by a processor, implements the cluster well platform location optimization method of the present invention.

The method considers the influence of the restricted area on the optimization of the platform position in the offshore oil and gas field development process, and based on the particle swarm algorithm, the method for optimizing the minimum total footage platform position outside the restricted area breaks through the defect that the existing method for optimizing the minimum total footage platform position does not consider the restricted area (the platform position calculated by the existing method is probably in the restricted area), and provides guarantee for reducing the drilling footage and the cost for the offshore oil and gas field development.

Drawings

FIG. 1 is a schematic (two-dimensional, initial) diagram of a preferred method of cluster well platform placement in accordance with one embodiment of the present invention.

FIG. 2 is a schematic representation (three-dimensional, results) of a preferred method of cluster well platform placement in accordance with one embodiment of the present invention.

Detailed Description

The experimental procedures used in the following examples are all conventional procedures unless otherwise specified.

Materials, reagents and the like used in the following examples are commercially available unless otherwise specified.

The invention considers the influence of the restricted area on the optimization of the platform position in the offshore oil and gas field development process, provides the optimization method of the cluster well platform position under the restriction of the restricted area based on the particle swarm algorithm, breaks through the defect that the existing optimization method of the minimum total footage platform position does not consider the restricted area (the platform position calculated by the existing method is probably in the restricted area), and provides guarantee for reducing the drilling footage and the cost for the offshore oil and gas field development.

The method is used for simulating a bird group foraging flight route based on a particle swarm algorithm, an individual experience and a group experience are learned and are cooperated with each other to iteratively find an optimal solution, the algorithm is easy to program and realize, and after target point coordinates and a limited area boundary are input, a program is operated, so that the platform position with the minimum total drilling footage under the constraint of the limited area can be obtained.

In order to achieve the purpose, the invention adopts the following technical scheme:

1) and determining the number n of wells, the coordinates of the target points, the vertical depth of the target points and the boundary range of the limited area of the target oil-gas field.

2) The method comprises the steps of presetting the number N of particles (representing N platform positions), the maximum iteration number T, learning factors C1 and C2, the maximum value Wmax and the minimum value Wmin of inertia weight, and the maximum value Vmax and the minimum value Vmin of speed.

3) Initializing the positions x (i, j) of N particles, wherein i represents the ith particle, j is 1 or 2, x (i,1) represents the east-west coordinates of the particle, x (i,2) represents the north-south coordinates of the particle, and the positions of the particles are randomly distributed in the expanded range of the oil field and outside the limited area in the initial stage;

the velocities v (i, j) of the N particles are initialized, the initial velocity is between Vmin and Vmax, a positive velocity indicates an increase in x (i, j), and a negative velocity indicates a decrease in x (i, j).

4) Calculating an initial optimal position p (i, j) and an initial minimum total footage pbest (i) of each particle, wherein the initial optimal position p (i, j) of each particle is an initial position x (i, j) of the 3 rd step, the initial minimum total footage pbest (i) is a drilling total footage corresponding to each position, and the pbest (i) is obtained after the calculation formula.

5) Calculating the global optimal platform position g and the global minimum total footage gbest of the whole particle swarm in the initial stage:

traversing from the 1 st particle to the Nth particle, selecting the minimum total footage pbest (i) of all the particles as the initial global minimum total footage gbest, and taking the corresponding platform position p (i, j) as the initial global optimal position g.

6) And (4) carrying out a first iteration (T rounds in total), updating the positions x (i, j) of all N particles to obtain new x (i, j), recalculating the total footage pbest (i) of each particle by using the step 4), selecting the minimum value as the global minimum total footage gbest after the first round of iteration, newly using the corresponding particle position x (i, j) as the global optimal platform position g after the first round of iteration, and carrying out the updating steps of the positions x (i, j).

7) And then performing iteration from 2 to T, and sequentially recording the global optimal platform position g and the global minimum total footage gbest of each iteration, wherein the g and the gbest after the T iteration are the final global optimal platform position and the global minimum total footage.

In the above method, in step 2), the number N of particles is set to 2 × m, the maximum number of iterations T is set to 10 × m, the learning factor C1 is set to C2 is set to 1.5, the maximum recommended inertial weight Wmax is set to 0.9 and the minimum recommended inertial weight Wmin is set to 0.4, and the maximum speed Vmax is set to 1 and the minimum speed Vmin is set to-1.

Where m represents the area of the oilfield region, km²。

The parameter setting in the above steps is composed of a fixed value and a function value, wherein learning factors C1 and C2, inertia weight maximum value Wmax and minimum value Wmin, and speed maximum value Vmax and minimum value Vmin are fixed values; the number N of particles and the maximum iteration time T are not fixed values but functions changing along with the area of the oil field, so that the value can ensure the operation of the algorithm and the accuracy of the result.

In the conventional particle swarm algorithm, the values of the number N of particles and the maximum iteration number T are fixed values, but the value taking method can cause that the number of the particles is not enough to cover the oil field area and the iteration number is not enough, so that the algorithm cannot be converged, and therefore, a function value mode is adopted in the invention to ensure that a program can normally run. In addition, when the maximum value and the minimum value of the speed are the fixed values, the calculation convergence speed is high, and the accuracy is high.

In the method, the initial position of each particle in the step 3-4) is within a range of 10 times of enlargement of the oil field area (called 'enlargement range') and is randomly distributed outside the limited area, and the times can be adjusted according to actual conditions.

In the above method, in step 4), the corresponding drilling total advancing length pbest (i) of the i-th particle is as follows:

wherein n is the number of wells; i is 1 to n; dep(s) represents the footage of the s-th well, including straight well sections, deflecting sections, steady deflecting sections and the like.

In the above method, in step 6), the updating step of the particle position x (i, j) is as follows:

first, calculating a dynamic inertia weight w:

w＝Wmax-(Wmax-Wmin)×k/T

wherein k represents the kth iteration;

the velocity v (i, j) of the particle is updated as follows:

v (i, j) new w × v (i, j) + C1 × rand x [ p (i, j) -x (i, j) ] + C2 × rand x [ g-x (i, j) ]

Wherein, C1 and C2 are learning factors, v (i, j) ═ rand x (Vmax-Vmin) + Vmin;

if the obtained v (i, j) is not between Vmin and Vmax, the regularization processing is carried out:

the position x (i, j) of the particle is then updated, as follows:

x (i, j) new ═ x (i, j) + vx (i, j)

If the obtained x (i, j) is not within the expanded range and is outside the restricted area, the regularization processing is carried out:

wherein Xmin, Xman, Ymin and Ymax are the boundaries of the west, east, south and north coordinates of the expanded range respectively.

And if the global minimum total footage gbest of the particle swarm after the position updating is smaller than the global minimum total footage gbest of the previous round, recording the global optimal platform position g corresponding to the global minimum total footage gbest. And after T-round iteration, obtaining a final global optimal platform position g and a final global minimum total footage gbest.

In step 6), the velocity v (i, j) is updated by calculating the dynamic inertia weight w, such as v (i, j) new w × v (i, j) + C1 × rand x [ p (i, j) -x (i, j) ] + C2 × rand x [ g-x (i, j) ], where the first term is the historical inertia, the second term is the individual experience, and the third term is the group experience. In the step, individual experience and group experience are comprehensively learned to obtain new speed, so that the fault tolerance is stronger. The inertia weight w is a dynamic value which changes along with the iteration turns, and can provide a proper inertia weight in each iteration turn, thereby preventing the phenomenon that the program can not be converged due to a static value.

And 6-7), further calculating x (i, j) new after v (i, j) new is obtained, calculating the global minimum total footage by using the particle swarm after the position is updated, and recording the corresponding global optimal platform position g if the global minimum total footage gbest is smaller than the global minimum total footage gbest of the previous round. And after T-round iteration, obtaining a final global optimal platform position g and a final global minimum total footage gbest.

And the final global optimal platform position g represents that after the influence of the restricted area is considered, the platform is selected at the position, the total drilling footage is minimum, and the total footage is gbest. Through practical verification, the setting of the parameters (including fixed values and function values) can ensure stable program operation and high convergence rate.

The following describes a preferred process of the method of the present invention, taking a certain oil and gas field as an example:

1) the number of the wells of the target oil-gas field is determined to be 16 directional wells, and the coordinates and the vertical depth of the target point are shown in table 1.

In this example, the type of the restricted area is a channel area, and the channel boundary is the restricted area boundary. The relative positions are shown in fig. 1.

TABLE 1 target point coordinates and target point vertical depth of directional well

2) Area of about 30km in oil field area²Therefore, the number of particles N is set to 60, the maximum number of iterations T is set to 300, the learning factor C1 is set to C2 is set to 1.5, the maximum value Wmax of the inertia weight is set to 0.9 and the minimum value Wmin is set to 0.4, and the maximum value Vmax is set to 1 and the minimum value Vmin is set to-1.

3) The positions x (i, j) of 60 particles are initialized, as shown in table 2, x (i,1) represents the east-west coordinates of the ith particle, x (i,2) represents the north-south coordinates of the ith particle, and the positions of the particles are randomly distributed within the oil field expansion range and outside the restricted area in the initial stage. As shown in fig. 1.

Table 2 particle position initialization

Particle number i	East-west coordinate x (i,1)	North-south coordinate x (i,2)
			1	387125	4352260
2	385082	4345976
			3	394274	4353569
…	…	…
			60	392595	4353149

The velocities v (i, j) of the 60 particles are initialized, the initial velocity being between Vmin and Vmax.

4) The initial optimal platform position p (i, j) and the initial minimum total footage pbest (i) for each particle are calculated.

The initial optimal platform position p (i, j) of each particle is the initial platform position x (i, j) in the step 3, and the footage of all 16 wells is accumulated to obtain the initial minimum total footage pbest (i) corresponding to each particle, as shown in table 3.

TABLE 3 initial optimal plateau position and initial minimum Total footage for each particle

5) And calculating the global optimal platform position g and the global minimum total footage gbest of the whole particle swarm in the initial stage.

Traversing from the 1 st particle to the 60 th particle, selecting the minimum pbest (i) of the 60 particles as an initial global minimum total footage gbest, and taking the corresponding particle position x (i, j) as an initial global optimal platform position g.

In this example, the initial global minimum total footage gbest is 91633m, and the corresponding initial global optimal platform position is shown in table 4.

TABLE 4 initial global optimal platform position and initial global minimum gross footage for a population of particles

Particle number i	East-west coordinate	North-south coordinates	Initial global minimum total footage/m
				29	393522	4355331	91633

6) A first iteration (300 total iterations) is performed in which the particle swarm locations are updated. If the particle position is out of the "expanded range" or enters the restricted area, then a corresponding adjustment is made according to step 6) in the embodiment to obtain new positions x (i, j) of the particle group that are both within the expanded range and outside the restricted area. Based on the new positions of 60 particles, the pbest (i) of each particle is recalculated using step 4), as shown in table 5.

TABLE 5 after iteration 1, new position and minimum Total footage for each particle

Particle number i	East-west coordinate	North-south coordinates	Minimum total footage/m
				1	392954	4350239	107121
2	392785	4350047	157172
				3	393452	4350363	93654
…	…	…
				60	392624	4353156	134339

And selecting the minimum value as the global minimum total footage gbest after the first iteration, and newly taking the corresponding particle position x (i, j) as the global optimal platform position g after the first iteration.

In this example, the global minimum total footage gbest after the first iteration is 84509/m. The corresponding global optimal platform position is shown in table 6.

TABLE 6 after 1 st iteration, global optimal platform position and global minimum total footage for particle swarm

Particle number i	East-west coordinate	North-south coordinates	Global minimum total footage/m
				29	393499	4355672	84509

7) And then performing 2 nd to 300 th iteration, and sequentially recording the global optimal platform position g and the global minimum total footage gbest of each iteration, wherein the g and the gbest after the 300 th iteration are the final global optimal platform position and the global minimum total footage. The final global optimal platform position and minimum overall footage in this example are shown in table 7.

TABLE 7 after 300 th iteration, global optimal platform position and global minimum gross footage for particle swarm

The results show that: under the restriction of the navigation channel, the platform position is selected to be (393348, 4350337), and the total drilling footage is minimum and is 71596 m. As shown in fig. 1 and 2.

The particle swarm algorithm based on the artificial intelligence bionic algorithm is an artificial intelligence bionic algorithm, through simulating foraging behavior of a bird swarm, comprehensively considering individual experience and swarm experience, iteratively searching for an optimal solution, and is high in searching speed, stable in operation after algorithm programming and reliable in result.

The above description is only an exemplary embodiment of the present invention, and should not be taken as limiting the scope of the invention, and any person skilled in the art should understand that they can make equivalent changes and modifications without departing from the concept and principle of the present invention. It should be noted that the components of the present invention are not limited to the above-mentioned whole application, and various technical features described in the present specification can be selected to be used alone or in combination according to actual needs, so that the present invention naturally covers other combinations and specific applications related to the present invention.