Next Article in Journal
Hegemony and Colonialization in the Water Management Sector: Issues and Lessons for IWRM
Previous Article in Journal
Research on Glacier Changes and Their Influencing Factors in the Yigong Zangbo River Basin of the Tibetan Plateau, China, Based on ICESat-2 Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization Method for Digital Scheduling of Oilfield Sewage System

1
School of Petroleum Engineering, Northeast Petroleum University, Daqing 163318, China
2
Postdoctoral Programme of Daqing Oilfield, Daqing 163318, China
3
Daging Oilfield Design Institute Co., Ltd., Daqing 163318, China
4
Key Laboratory of Continental Shale Hydrocarbon Accumulation and Efficient Development, Northeast Petroleum University, Daqing 163318, China
5
School of Computer and Information Technology, Northeast Petroleum University, Daqing 163318, China
*
Author to whom correspondence should be addressed.
Water 2024, 16(18), 2623; https://doi.org/10.3390/w16182623
Submission received: 6 August 2024 / Revised: 12 September 2024 / Accepted: 12 September 2024 / Published: 15 September 2024
(This article belongs to the Section Urban Water Management)

Abstract

:
Oilfield Sewage System Scheduling is a complicated, large-scale, nonlinear system problem with multiple variables. The complexity of the sewage system pipeline network connection grows along with the ongoing building of oilfield stations, and the shortcomings of the sewage system water quantity scheduling program based on human experience decision-making become increasingly apparent. The key to solving this problem is to realize the digital and intelligent scheduling of sewage systems. Taking the sewage system of an oil production plant in Daqing oilfield as the research object, the water scheduling model of the sewage system is established in this paper. Aiming at the complex nonlinear characteristics of the model, the Levy flight speed updating operator, the adaptive stochastic offset operator, and the Brownian motion selection optimization operator are established by taking advantage of the particle swarm optimization (PSO) and the cuckoo search (CS) algorithm. Based on these operators, a hybrid PSO-CS algorithm is proposed, which jumps out of the local optimum and has a strong global search capability. Comparing PSO-CS with other algorithms on the CEC2022 test set, it was found that the PSO-CS algorithm ranked first in all 12 test functions, proving the excellent solving performance of the PSO-CS algorithm. Finally, the PSO-CS is applied to solve a water scheduling model for the sewage system of an oil production plant in Daqing Oilfield. It is found that the scheduling plan optimized by PSO-CS has a 100% water supply rate to the downstream water injection station, and the total energy consumption of the scheduling plan on the same day is reduced from 879.95 × 106 m5/d to 712.84 × 106 m5/d, which is a 19% reduction in energy consumption. The number of water balance stations in the sewage station increased by 7, which effectively improved the water resource utilization rate of the sewage station.

1. Introduction

Oilfield sewage is a by-product of the oilfield extraction process, consisting of oilfield-produced water, drilling sewage, and other types of oily sewage in the sewage station [1]. It is the primary issue with the water quality in oilfields and presents a major risk to human health and the environment [2]. The sustainable growth and use of water resources are realized through the effective functioning of the oilfield sewage system, which also lowers environmental contamination and encourages clean oilfield production [3,4]. Though labor-intensive and time-consuming, the current oilfield sewage system water scheduling scheme primarily depends on human experience, which increases the likelihood of operational and technical issues. Moreover, it is challenging to guarantee the water scheduling scheme’s quality in real time [5,6]. Digital scheduling is based on the digitization of production demand, and dynamic water allocation is achieved daily through remote regulation. However, there is no scientific decision-making method for digital scheduling. Therefore, it is of great significance to carry out an optimization method for the digital scheduling of oilfield wastewater systems in order to improve the operational efficiency of oilfield wastewater systems, the utilization of water resources, and the protection of the environment.
Relatively few studies have been performed recently on the subject of oilfield sewage system water scheduling optimization. The research field of water resources scheduling includes the study of water scheduling in oilfield sewage systems, and many academics have worked extensively in this area to better understand how water resources are scheduled. Jia et al. [7] propose a multi-stage reservoir scheduling and water allocation model based on multi-stage coupling. The model can select more reasonable water transfer methods according to different development needs so that the water shortage rate of each user in the receiving area is kept below 8% and the reservoir idle rate is controlled at about 10%, which effectively improves the water resource allocation capacity. Yang et al. [8] establish a multi-objective optimal scheduling model and splits the multi-objective optimal scheduling model using a two-stage method in order to optimize the water transfer engineering scheme. Wu et al. [9]. proposed a multi-objective reservoir optimal scheduling model as a water transfer rule for inter-basin water transfer with the objectives of economic efficiency, power generation, and ecological water transfer. Wang et al. [10] developed a mixed-integer linear programming model considering water demand, pumping station transmission capacity, and pipeline operation constraints. The model minimizes the total cost of the oilfield water delivery system and improves water resource utilization. However, with the increasing complexity of the water resources scheduling operation model, traditional algorithms such as linear programming, nonlinear programming, and dynamic programming can no longer meet the development trend of the water allocation model [11].
Heuristic optimization methods, which may quickly converge to the global optimal solution in a complicated solution space, have become widely employed in water resource allocation research due to advancements in computer technology and artificial intelligence [12]. Under difficult restrictions, it offers an effective solution approach for the water resource allocation model [13]. Ezzeldin et al. [14] proposed a novel dual hybrid optimization algorithm model of Genetic Algorithm-Particle Swarm Algorithm to solve the problem of optimal design of urban water distribution networks to support urban water allocation. Deng et al. [15] proposed a deoxyribonucleic acid (DNA) algorithm as a biological algorithm for solving the rural postman’s problem, which solves the water allocation problem. Turci et al. [16] proposed an improved genetic algorithm to solve the pump scheduling model. Dong et al. [17] proposed a data-driven water pumping station unit optimization scheduling model based on the fusion of IoT data and PSO using a Long Short-Term Memory (LSTM) network for optimal scheduling of water pumping stations using a predictive control framework. The particle swarm algorithm and its improvement methods have become one of the most widely used optimization algorithms in water scheduling engineering due to its advantages, such as easy implementation, algorithmic simplicity, and high efficiency [13]. Eryiğit et al. [18] proposed a modified immune algorithm to solve an urban water supply scheduling model in the Iraqi desert for optimal allocation of urban water management. Ma et al. [19] developed a parallel particle swarm algorithm to solve a multi-objective scheduling model of a reservoir cluster in the Heihe River Basin. An adaptive stochastic inertia-weighted ARIW-PSO was proposed by Chen et al. [13] to address a multi-objective reservoir cluster scheduling model. PSO-GA was proposed by Wu et al. [20] as a solution to the Drinking River Basin multi-objective ecological dispatching model.
In this paper, an improved particle swarm optimization algorithm with a hybrid particle swarm-cuckoo (PSO-CS) algorithm is proposed in combination with a cuckoo intelligent optimization algorithm. The algorithm is used to compensate for the shortcomings of the standard PSO algorithm, which is easy to fall into the local optimum and has fast convergence speed. In order to accurately assess the reliability of the improved PSO, two non-parametric statistical tests were used: the Friedman test and the Bonferroni–Dunn test. The Friedman test is suitable for dealing with repeated measurements and is able to provide robust statistical results by effectively comparing the differences between multiple samples of interest without relying on the assumption of normal distribution. The Bonferroni–Dunn test, on the other hand, enhances the accuracy and reliability of the results when making multiple comparisons by adjusting the significance level to control the rate of false positives and reduce the risk of false findings. The combined use of these two tests ensures that our assessment of the performance of the improved particle swarm algorithm is more comprehensive and rigorous. Taking the production and operation situation of a sewage plant in Daqing oilfield as a case study, considering the oilfield sewage system process flow and water transfer rules, constructing a sewage system water quantity scheduling optimization model, and solving the sewage system water quantity scheduling model by using the improved PSO-CS algorithm. This study innovatively combines an improved particle swarm algorithm with a water scheduling model for oilfield wastewater systems. The optimization efficiency of the scheduling scheme is improved by introducing a swarm intelligent optimization algorithm. Compared with the traditional method, we used the improved particle swarm algorithm to enhance its adaptability and performance in complex wastewater treatment systems. The objective of this study is to provide an optimal water transfer scheme for the production and operation of oilfield wastewater systems through the improved particle swarm algorithm.

2. Materials and Methods

2.1. Overview of the Sewage System in the Study Area

This paper uses the Daqing oilfield in China as its study area. The study area’s sewage system is primarily made up of water-driven sewage treatment stations, poly-containing and ternary sewage treatment stations, deep sewage treatment stations, water injection stations, and a network of sewage pipelines. It consists of eleven water-driven sewage treatment plants, seven comprising poly and ternary treatment plants, nine depth treatment plants, and twenty-two water injection plants. To maintain the regular operation of the water injection system, the upstream sewage treatment station uses pipelines to deliver treated sewage containing oil to the downstream water injection station. Figure 1 depicts the sewage system’s overall layout. Water-driven, poly, and ternary sewage treatment stations with depth and water-driven capabilities to handle 5 mg/L of water quality and 20 mg/L of water quality, respectively. Using the plant’s daily sewage system production statistics from 2024 as an example, the system received 270,698 m3 of raw water each day, of which 164,964 m3 had a high concentration of poly sewage (carrying more than 150 mg/L of polymer). 105,734 m3 of low-containing poly sewage (with a poly content of less than 150 mg/L). The system discharged 109,646.3 m3, of which 63,264.95 m3 were ordinary water and 46,381.35 m3 were water at a depth. The scheduling scheme for this sewage system currently relies on manual experience, which makes it challenging to guarantee the system’s stable and efficient operation. It also faces issues with declining well production, the oilfield sewage system’s high energy consumption, high production costs, etc. Researching the improvement of the sewage system’s operating plan is necessary to boost the system’s operational efficiency and attain steady, uninterrupted oilfield operations.

2.2. Mathematical Model for Optimization of Water Dispatch in Sewage System

2.2.1. Objective Function

The goal of the oilfield sewage system schedule optimization is to identify the best scheme out of those available to satisfy the various needs of the system. The energy consumption of the system during the transfer process is used to weigh the benefits and drawbacks of each scheme. This paper uses the pipeline length, the product of the pipeline’s flow rate, and its pressure drop along the route to describe the energy consumption of the system, and this simplifies the model. The sewage system’s scheduling optimization process objective function may be stated as Formula (1).
min F ( Q ) = i = 1 m j = 1 n Q i , j × L i , j × Δ H i , j
Δ H i , j = β Q i , j 2 k L i , j V k / D i , j 5 k
In Formulas (1) and (2),
Q i , j —Volume of water in the pipeline from the sewage station i to the water injection station j in the system, m 3 / d ;
L i , j —Length of pipeline from sewage station i to water injection station j in the system, m ;
Δ H i , j —Pressure drop along the pipeline from the sewage station i to the water injection station j in the system, m ;
D i , j —Pipeline diameter from sewage station i to water injection station j in the system, m ;
V —Kinematic viscosity of the fluid in the pipe, m 2 / s ;
m —Number of sewer stations in the system, se;
n —Number of water injection stations in the system, se;
β —Pressure drop calculation factor;
k —Fluid state index, values are shown in Table 1;

2.2.2. Constraint Conditions

The restrictions between the transfer pipeline network are established starting with the injection and sewage systems combined. Determining the volume of water to be moved from each sewage treatment station reservoir to the downstream injection station within the intricate network structure is the goal of the sewage system dispatch optimization strategy. Combining the needs and operational features of the water injection system and the sewage system results in the following constraints.
  • Supply and demand constraints
The sum of the water volume of the upstream sewage treatment station of the oilfield sewage system is greater than or equal to the sum of the water required by each downstream water injection station.
i = 1 m a i j = 1 n b j
In Formula (3),
a i —Volume of water from sewage station i , m 3 / d ;
b j —Water requirements for water injection station j , m 3 / d ;
2.
Water quantity constraints
The total water output of any sewage station is equal to the sum of the water flow out of the station, which is equal to the sum of the flow in the pipeline in the direction of its downstream conveyance water, and the total water demand of any water injection station is equal to the sum of the water flow into the station, which is equal to the sum of the flow in the pipeline in the direction of its upstream incoming water.
j = 1 n Q i , j = a i , i = 1 , 2 , 3 , ... , m
j = 1 m Q i , j = b i , j = 1 , 2 , 3 , ... , n
Q i , j 0 , i = 1 , 2 , 3 , ... , m , j = 1 , 2 , 3 , ... , n
In Formulas (4)–(6),
Q i , j —Volume of water in the pipeline from the sewage station i to the water injection station j in the system, m 3 / d ;
a i —Volume of water from sewage station i , m 3 / d ;
b j —Water requirements for water injection station j , m 3 / d ;
m —Number of sewage stations in the system, se;
n —Number of water injection stations in the system, se;
3.
Injection station effluent concentration constraints
One of the indications used to assess whether the water transfer plan satisfies production demand is the poly concentration of sewage, which is a crucial parameter in the sewage system dispatch optimization process. Each downstream injection station in the sewage system must have incoming sewage concentrations that are either less than or equal to the injection concentration requirements of that station.
i = 1 m Q i , j × η i i = 1 m Q i , j η j , max , j = 1 , 2 , 3 , ... , n
In the Formula (7):
η i —Concentration at sewage station i , mg / L ;
η j , max —Maximum effluent concentration required at the jth injection station, mg / L ;
Q i , j —Pipeline flow from sewage station i to water injection station j , m 3 / d ;
m —Number of sewage stations in the system, se;
n —Number of water injection stations in the system, se;
4.
Water injection station pressure constraints
The inlet pressure of each pipe at the injection station downstream of the sewage system is greater than or equal to the required pressure at the injection station.
P j i n P j n
In Formula (8),
P j i n —Inlet pressure at the injection station j in the system, Mpa;
P j n —Required injection pressure at the j th injection station in the system, Mpa.

2.3. PSO-CS Algorithm

2.3.1. PSO Algorithm

PSO algorithm is a method to optimize continuous nonlinear function based on iterative updating of random population, which originates from the study of bird predation behavior [21]. The main process is that, in each iteration cycle, the potential position of each particle is randomly distributed between the individual historical optimal position and the population historical optimal position, and PSO uses the velocity–position update formula to make each candidate solution follow the current optimal solution, and ultimately find the global optimal solution after many iterations [22]. The mathematical model of the standard PSO is as follows.
The symbols x and v represent the position of the particle and the search velocity of the particle, which correspond to the solution of the problem and the search velocity of the solution, respectively. Suppose that in a D-dimensional target search space, there are N particles forming a population, the solution set is X = ( x 1 , x 2 , , x N ) T , where x i = ( x i 1 , x i 2 , , x i D ) is the i-th solution and v i = ( v i 1 , v i 2 , , v i D ) is the search velocity of the current solution. p b e s t i = ( p b e s t i 1 , p b e s t i 2 , , p b e s t i D ) and g b e s t i = ( g b e s t i 1 , g b e s t i 2 , , g b e s t i D ) are individual historical optimal solutions and global historical optimal solutions, respectively. The velocity iteration update formula is based on Equation (9), and the solution iteration update formula is based on Equation (10).
v i ( t + 1 ) = ω ( t ) × v i ( t ) + c 1 × r 1 ( t ) ( p b e s t i ( t ) x i ( t ) ) + c 2 × r 2 ( t ) × ( g b e s t i ( t ) x i ( t ) )
x i ( t + 1 ) = x i ( t ) + v i ( t + 1 )
where, ω ( t ) is the inertia coefficient of the solution x i , and the value ranges from 0.4 to 0.9 [23]. c 1 and c 2 are individual learning factors and social learning factors respectively, and both of which are positive numbers.
r 1 ( t ) and r 2 ( t ) are uniform random numbers between 0 and 1.
The steps of the PSO optimization algorithm are described as follows:
Step 1: Randomly initialize the position x i ( 0 ) and velocity v i ( 0 ) (water flow rate of each pipeline from the sewage station to the injection station) of each particle, calculate the fitness value of x i ( 0 ) (the objective function value of the mathematical model for optimization of the sewage system scheduling), and determine the values of p b e s t i and g b e s t i .
Step 2: Update the positions and velocities of the particles according to Formulas (9) and (10), and then calculate the water dispatch energy consumption according to Formulas (1) and (2).
Step 3: Evaluate the fitness value of x i ( t + 1 ) . If the value is desirable (smaller value of the sewage system scheduling optimization objective function), then replace p b e s t i ( t ) and g b e s t i ( t ) with the new p b e s t i ( t + 1 ) and g b e s t i ( t + 1 ) .
Step 4: If t < T , then t = t + 1 and return to step 2; otherwise, optimization stops.

2.3.2. CS Algorithm

The cuckoo search algorithm (CS) is a meta-heuristic algorithm proposed to mimic the natural phenomenon of cuckoos laying eggs on borrowed nests. In the classical CS algorithm, the cuckoo nest location update is based on Levi’s flights [24]. In order to facilitate the simulation of the cuckoo’s reproduction strategy, the cuckoo algorithm will be assumed to be in the following three ideal states [25].
(1)
Each cuckoo lays one egg at a time and randomly selects a location for the nest to incubate;
(2)
In a randomly selected group of nests, the high-quality nests will be retained for the next generation;
(3)
The number of available nests is fixed, and the probability that the nest owner discovers an alien cuckoo egg is Pa ∈ [0, 1]. When a nest host discovers an alien cuckoo egg, it will either discard the cuckoo egg or re-establish a new nest.
Based on the above three rational states, the updated formula for the cuckoo’s location and path to find the host nest is as follows.
x i ( t + 1 ) = x i ( t ) + α L e v y ( β )
In Formula (11), α is the step control factor, when α > 0 , then the typical value α = 1 ; L e v y ( β ) denotes the Lévy random search path, where β is the control parameter of the path, and the specific formula for L e v y ( β ) is shown below.
L e v y ( β ) = ϕ u v 1 β
In the formula, ϕ = Γ ( 1 + β ) × sin ( π × β / 2 ) Γ [ ( 1 + β ) / 2 ] × β × 2 ( β 1 ) / 2 1 / β ; u , v is a random number obeying a normal distribution.
The number of available nests is fixed, and the probability that the nest owner discovers an alien bird’s egg is Pa ∈ [0, 1]. After updating the location, a random number r is compared with Pa. If r > pa, the nest location is updated randomly once; otherwise, the nest location remains unchanged. The updated formula for the location of the discovered nest is as follows.
x i ( t + 1 ) = x i ( t ) + r ( x j ( t ) x k ( t ) )
In Formula (13), r is the scaling factor; x j ( t ) , x k ( t ) is the random position of the nest in the t th generation.

2.3.3. PSO-CS Hybrid Algorithm

1.
Feasibility Analysis
The optimization of individual search behavior is the key to improving the PSO-CS intelligent optimization algorithm. To demonstrate intuitively that the optimization performance will be effectively enhanced by mixing the two algorithms, Figure 2 is plotted to discuss the optimization mechanism of the mixing of the two algorithms. Particle i moves to point 2 for the traditional PSO algorithm, as illustrated in Figure 2. The particle is then guided to point 3 by the historical optimal position and moves to point 4 from point 3 under the combined influence of social cognition and the current optimal position. That is, if no improvement occurs, the particle can move from point 1 to point 4. A large step movement can be performed by the individual with a certain probability when they move to point 1 of the PSO-CS algorithm. This step movement incorporates the Levy flight’s step movement and is drawn to the historically optimal position to point 5. Then, the individual performs an adaptive stochastic offset operator, taking into account the effect of the current group densities, using the offset operation to move point 6 from point 5. Ultimately, the PSO-CS algorithm’s addition of the Brownian motion selection operators causes the current optimal individual to perform a random Brownian selection operation within a specific neighborhood, effectively enhancing both the optimal position and the current optimal individual. As a result, the individual moves from point 6 to point 7. That means after adding the three optimization operators, the individual moves from point 1 to point 7. In comparison to the traditional PSO algorithm, the hybrid algorithm more effectively incorporates the CS algorithm’s stochastic search capability. Additionally, it increases the individual’s global optimization searchability and the probability of jumping out of the local optimal solution under the influence of the random offset and the current optimal individual selection optimization. This can also cause the algorithm to jump out of the local optimum, as demonstrated by the hybrid algorithm’s individual jumping out of the local optimal region (represented by the blue region in the figure). The aforementioned study leads to the conclusion that combining the PSO and CS algorithms to create a superior algorithm works well.
2.
Levy Flight Speed Update Operators
Levy flight generates a class of random motion trajectories distinct from Brownian motion and is a random walk with a heavy-tailed probability distribution of step lengths. Its movement pattern matches the pattern of birds that forage in one place for food but travel to other regions when food becomes scarce. It also moves with a high step length with a given probability. Levy flight is incorporated into particle updating based on the PSO algorithm’s process architecture. This results in the formation of a new Levy flight velocity updating operator by converting the conventional constant learning factor into Levy flight coefficients.
For the PSO algorithm, the speed update of particles is mainly performed concerning the historical optimal individuals and the current global optimal individuals. Although it accelerates the speed of the population to discover high-quality solutions to a certain extent, it causes the population to gather around excellent individuals and fall into the local optimum. To overcome the shortcomings of the PSO algorithm that the diversity decreases in the late iteration and makes the individuals jump out of the local optimum, combined with the advantage of Levy flight that can search for large step size, the following formula for the speed update of the individuals in the hybrid algorithm is obtained.
v i ( t + 1 ) = w v i ( t ) + L v r 1 [ p b e s t i ( t ) x i ( t ) ] + c 2 r 2 [ g b e s t ( t ) x i ( t ) ]
In this formula, L v is the Levy flight coefficient, which is calculated as follows:
L v = μ p Γ ( 1 + β ) × sin ( π × β / 2 ) Γ [ ( 1 + β ) / 2 ] × β × 2 ( β 1 ) / 2 1 / β q 1 β
In Formula (15), β is the control parameter of the Levy flight path, μ is the Levy flight control factor, and p , q is a random number obeying a normal distribution.
3.
Adaptive Stochastic Offset Operators
By using the concept of randomly discarding individuals from the CS algorithm, the adaptive random offset operator selects a certain number of individuals from the population—excluding the currently optimal individuals—for random offset after each generation of the hybrid algorithm. This process can significantly increase the population’s capacity to explore all possible domains. The selection of these individuals is not according to the level of fitness value but according to the random selection. A different offset group can be created by randomly choosing a specific number of people. The alternative offset group can then use the change rule of group diversity to adjust the corresponding offset probability and assess whether or not to perform the offset operation. The effect of group density is taken into account when calculating the probability of offset individuals by measuring the size of the group diversity. For each individual in the alternative group, it is determined one by one if their offset probability crosses the threshold; if it does, the offset operation is carried out, as indicated by the discriminant formula of the offset operation of the adaptive optimization process.
P A x i ( t ) = e κ / min x j ( t ) x k ( t ) > λ A , min x j ( t ) x k ( t ) 0
In this formula, P A ( ) is the probability of the i th individual’s bias; κ is the bias control parameter; min ( ) denotes the minimum value of the formula taken over the inter-vectorial formula of individuals in parentheses; x j ( t ) , x k ( t ) is the randomly selected individual in the population of the hybrid algorithm; and λ A is the bias threshold.
For the individuals whose offset probability is greater than the threshold, the offset is performed based on the following formula to obtain new individuals, and the new individuals are mixed into the original population to continue performing the subsequent iterations.
x i ( t ) = x i ( t ) + r A ( t ) x j ( t ) x k ( t )
In the formula above, x i ( t ) is the vector form of the individual performing the offset operation through the offset discrimination; x i ( t ) is the vector form of the individual after performing the offset operator; x j ( t ) , x k ( t ) is the vector of the two randomly selected individuals in the discrimination formula; and r A ( t ) is the scaling factor, which is calculated as shown in the following Formula (18):
r A ( t ) = r min + ( r max r min ) π t 2 I max
In this formula, r A ( t ) is the scaling factor of the t th iteration, r min is the minimum value of the scaled factor, r max is the maximum value of the scaling factor, and I max is the maximum number of iterations.
4.
Brownian Motion Selection Operators
Brownian motion random search is used to optimize the entire population after improving the current optimal individual, hence expanding the algorithm’s global search capability. The selection space of the current optimal individual is generated by the Brownian motion using a random step. Each state in the space is a possible state of the current optimal individual, and the selection space is denoted by the { g b 1 , g b 2 , , g b m B } . The current optimal individual and the step of the Brownian motion combine to form each state in the selection space, with g b j representing the j th state and m B denoting the number of states. Each state in the selection space consists of a combination of the current optimal individual and the Brownian motion step, which is formulated as follows:
g b j = g b ( t ) + χ B s ( t )
In this formula, g b ( t ) is the current optimal individual of the t th iteration, χ is the Brownian motion step size control factor, and B s is the approximate Brownian motion step size, where each dimension of the step size is calculated as shown in the following formula:
B s , j ( t ) = k = 0 N a k 2 k + 1 2 sin ( k + 1 2 ) π t
B s , j ( t ) is the value of the j th dimension of the approximate Brownian motion step; a is a random number from a standardized nontrivial distribution; and N is a positive integer corresponding to the Brownian motion.
Based on the above selection space, the best state adapted by the current optimal individual is selected from this space, and the specific formula for the optimal selection has been given.
g b ( t ) = min f i t ( g b j ) g b j { g b 1 , g b 2 , , g b m B }
In this, g b ( t ) is the current optimal individual after executing the Brownian motion selection operator; f i t ( ) is the fitness function to find each state in the selection space.
5.
Hybrid Algorithm Solution Process
The PSO-CS algorithm is a novel kind of swarm intelligent optimization algorithm that combines the PSO and CS algorithms. It uses the PSO algorithm as its foundation, and its primary parameters and initial values are consistent with those of the classical algorithm. The PSO-CS algorithm’s main steps are stated as follows:
(1)
Set the initial values for the inertia weights and learning factors in the classical particle swarm and cuckoo algorithms; set the initial values for the Lévy flight speed updating operator, the adaptive stochastic offset operator, and the Brownian motion selection operator. These include the Lévy flight speed updating operator’s control parameter β and control factor μ , the adaptive stochastic offset operator’s offset control parameter κ and the offset threshold λ A ; the Brownian motion step length, population size, and termination conditions; read and store the objective function and constraints;
(2)
Generate an initial population of hybrid algorithms, compute the fitness function value, and store the historical optimal individual x p b ( 0 ) and the current global optimal individual x g b ( 0 ) ;
(3)
Calculate the Levy flight coefficients according to Formula (15), which in turn updates the velocities and positions of all individuals based on the Levy flight velocity update operator;
(4)
Judge whether the updated individual satisfies the constraints; if yes, go to step (6); otherwise, go to step (5);
(5)
Adjust the individuals that do not meet the constraints;
(6)
Calculate the fitness function values of all individuals in the hybrid algorithm population and update the historical optimal individual x p b ( t ) and the current global optimal individual x g b ( t ) ;
(7)
Based on the current global optimal individual, combine the Brownian motion step formula to form the potential states of the current optimal individual, construct the state space, calculate the fitness function value corresponding to each state, and update the current global optimal individual optimally;
(8)
Calculate the fitness function value of the hybrid algorithm population to update the historical optimal individual x p b ( t ) and the current global optimal individual x g b ( t ) ;
(9)
Judge whether the termination condition is satisfied; if yes, go to step (12); if not, go to step (3);
(10)
Output the global optimal solution. The PSO-CS algorithm flow is shown in Figure 3.

2.4. PSO-CS Algorithm Tests Based on CEC2022 Test Function Sets

The optimization results of the PSO-CS algorithm are compared with those of the PSO, CS, WOA, SA, COA, and HO algorithms based on the 12 test functions of the CEC2022 optimization function test set to confirm the effectiveness of the PSO-CS method. To ensure fairness in the comparison, each method has a population size of 100, a maximum number of iterations of 1000, and is run separately 20 times. The final optimization result is determined by averaging the 20 optimization outcomes. The accuracy and stability of each method’s optimization are reflected by the optimal value mean and standard deviation, respectively, and the solutions of each algorithm are ranked. Table 2 displays the outcomes of the comparison.
Table 2 shows that, when it comes to the mean and standard deviation of the optimal values that are acquired after 1000 iterations of solving, the PSO-CS method performs better than the other six algorithms. The PSO-CS algorithm works particularly well on the functions F 8 , F 10 and F 12 . Furthermore, the PSO-CS algorithm receives the highest score among all the functions that are examined, indicating that it is more successful at enhancing the PSO algorithm’s optimization.
The convergence curves of PSO-CS compared to other algorithms on the CEC2022 test set are presented as indicated in Figure 4 to visualize the convergence performance of the PSO-CS algorithm. When compared to other methods, PSO-CS exhibits faster and more accurate convergence in the majority of the test functions. This suggests that the optimization technique used in this research can successfully increase PSO’s convergence speed and accuracy, leading to superior optimization outcomes. At the same time, other algorithms easily fall into local optimization, which leads to slower convergence or even stagnation.
In addition to analyzing the numerical results obtained by the PSO-CS algorithm and other algorithms, the superiority of the PSO-CS algorithm is demonstrated by two non-parametric statistical tests, the Friedman test and the Bonferroni–Dunn test.
A multiple comparison test called the Friedman test is used to see whether there are any appreciable variations between algorithms. The Friedman test’s algorithm ranking principle states that the algorithm with the highest performance is ranked the smallest, while the algorithm with the lowest performance is ranked the largest. The optimal values and standard deviations produced by various methods in Table 2 were subjected to the Friedman test; the computed outcomes are displayed in Table 3.
Friedman’s test results in Table 3 show that all p-values are below the significance level of α = 0.01, suggesting that these seven algorithms differ significantly from one another. The PSO-CS algorithm performs the best in terms of optimizing the solution, followed by the CS, WOA, and other algorithms. The test results of the mean and standard deviation of the optimal values also show this. Therefore, it can be said that the PSO-CS method has greater accuracy than the other algorithms. Nevertheless, the Bonferroni–Dunn test is used to assess the optimization performance of the PS-FW algorithm because Friedman’s test can only determine whether there is a significant difference between all algorithms from an overall perspective and cannot specifically compare the performance difference between the PSO-CS algorithm and each other algorithm.
Finding significant differences between two or more algorithms may be achieved quite simply with the help of the Bonferroni–Dunn test. The following equation determines the critical difference, which must be exceeded by the performance rankings of the two algorithms for there to be a significant difference between them for the Bonferroni–Dunn test [26].
C D α = q α N i ( N i + 1 ) 6 N f
In this formula, N i and N f are the number of algorithms and benchmark functions, respectively, q α is the critical value at the significance level α , which is as follows at different significance levels:
q 0.05 = 2.77 ,   q 0.1 = 2.54
Combining Formulas (22) and (23) yields the following critical differences at different significance levels:
C D 0 . 05 = 2 . 44 ,   C D 0 . 1 = 2 . 24
The Bonferroni–Dunn test was conducted on the optimal value mean and optimal value standard deviation, utilizing the ranks acquired from Friedman’s test. The histograms based on the optimal value mean and optimal value standard deviation are plotted as shown in Figure 5. For more intuitively displaying the results of the Bonferroni-Dunn test combined with the critical difference. The pink horizontal line represents the ranking of the best algorithms, and the two black and red horizontal lines, respectively, correspond to the threshold level values under the significance level of α = 0 . 05 and α = 0 . 1 , which are equal to the sum of the values for the critical difference and the lowest ranked value. An algorithm is shown to perform worse than the optimal ranking algorithm if the height of its histogram is greater than the threshold level’s horizontal line. The PSO-CS algorithm’s histogram is the lowest of all the algorithms, as shown in Figure 5a. The PSO-CS algorithm outperforms the PSO, WOA, SA, COA, and HO algorithms at both the 90% and 95% significance levels. In addition, PSO outperforms CS but non-significantly. This indicates that the PSO-CS algorithm incorporates the correctness of the three optimization operators and performs well in terms of global optimization solving power. Figure 5b shows that the PSO-CS algorithm has relatively better robustness compared to other improved swarm intelligence optimization algorithms. It outperforms PSO, WOA, SA, COA, and HO algorithms for the PSO-CS algorithm at 90% and 95% confidence levels. It shows that the PSO-CS algorithm is significantly better than PSO, WOA, SA, COA, and HO algorithms in terms of the stability performance of the algorithm in solving the optimization problem. CS algorithm is also better in terms of robustness, so the PSO-CS algorithm does not significantly outperform the NRBO algorithm in the Bonferroni–Dunn test.

3. Results and Discussions

3.1. Scheduling Optimization Results for Sewage System on a Certain Day

3.1.1. Scheduling Optimization Scheme

The sewage and injection stations for this sewage system, as well as the water pipeline paths and pipeline distances between them, are shown in Table 4.
The sewage system operating model is solved using the hybrid PSO-CS optimization technique presented in this research, taking into account the production data of the sewage system on that particular day. The sewage system’s pre- and post-optimization water scheduling schemes are depicted in Figure 6 and Figure 7, respectively. It is clear from these that the optimized water transfer scheme has undergone significant modification in comparison to the pre-optimization one. The improved water transfer strategy is more focused on the injection station’s water supply, whereas the pre-optimization water transfer scheme has scattered distribution and a lower water volume. Using the sewage station “NIII-1” as an example, water is sent to injection stations downstream in four directions: “ZN8”, “ZN7”, “ZN19”, and “ZN10”. The water supply plan before optimization indicates that “NIII-1” exclusively provided water to “ZN8”, “ZN7”, and “ZN19”. Following optimization, “NIII-1” solely provides water to the “ZN10” injection station; the volume of water supplied is 6559 m3, greater than that of “ZN19”. The water amount is 6559 m3, which is more than it was before optimization, and “NIII-1” solely provides water to the “ZN10” injection station after optimization. This demonstrates how the optimized scheduling technique successfully raises operating efficiency, lowers energy consumption, and enhances the sewage system’s water scheduling strategy.

3.1.2. Comparison and Result Analysis before and after Optimization

Based on the water scheduling results for the sewage system, compute the energy consumption of the sewage system scheduling before and after optimization, as illustrated in Figure 8. The total operational energy consumption of the optimized system is 712.84 × 106 m5/d, and the operational energy consumption of the optimized sewage system is reduced by 19% when compared to the system operation before optimization. The total operational energy consumption generated by the scheduling plan on that day is 879.95 × 106 m5/d.
Table 5 displays the volume of water needed for the production and operation of each injection station on that particular day. The water supply for each injection station can be provided by the scheduling plan before and during optimization, as seen in Figure 9. For example, “ZN19” and “ZNII-1” require 5930 m3 and 8352 m3, respectively, and the water supply from the upstream sewage station before the optimization is 3802 m3 and 4640 m3, respectively. These amounts represent 64.1% and 55.6% of the required water quantity. Based on the optimized scheduling plan, the water supply from the injection stations “ZN19” and “ZNII-1” is 5930 m3 and 8352 m3, respectively, which accounts for 100% of the injection stations’ water demand. This indicates that the optimized scheduling plan can effectively meet the injection stations’ water demand and can support the water injection system’s efficient operation.
There were 27 sewage stations produced and running on that day. The water quantity entering and exiting the sewage stations was determined based on the day’s production statistics and optimization outcomes, as indicated in Table 6. The water quantity in and out of the sewage stations was not adequately balanced by the current scheduling plan, which prevented some of the inlet water from being effectively transferred to the downstream injection station. This led to the hoarding of water in the sewage stations. Before optimization, the water quantity out of the sewage stations was typically lower than that of the inlet water. The amount of water in and out of 7 sewage stations reaches a balanced state after optimization, and the amount of water out of the sewage station rises significantly, with the majority of the water out of the sewage station being close to or equal to the amount of inlet water. This indicates that the water transfer scheme effectively improves the utilization rate of water resources in the sewage stations, and the oilfield sewage system’s operation efficiency has been greatly enhanced.

3.2. Optimization Results of Sewage System Scheduling for a Particular Month

The production data of the plant in June is utilized as the dataset to optimize the water transfer scheme for 30 production days of the month to confirm the stability of the PSO-CS optimization algorithm. Figure 10 displays the injection stations’ water supply both before and after optimization. There are few injection stations that meet the water requirement before optimization. In particular, on Days 9, 10, and 18, the water transfer scheme only provides for the water requirements of 15 injection stations, or 68.18% of the total number of injection stations. This has an impact on the downstream injection system’s ability to operate efficiently. Furthermore, the optimized water transfer scheme demonstrates that the number of water injection stations meeting the water requirement is 22 per day, which can effectively ensure the normal operation of the downstream water injection system. The above analysis demonstrates that the water scheduling scheme using human experience cannot adapt to the changes in the system operating conditions, resulting in the inability to effectively meet the water requirements of the injection stations. Using the intervals from Day 17 to Day 18 and Day 10 to Day 11 as examples, the water transfer scheme also demonstrates that the number of injection stations meeting the water quantity requirement fluctuates greatly during these two periods, indicating that the stability of the water transfer scheme is low. The ideal water supply scheme can be modified in response to modifications in system operating conditions, guaranteeing the stability of the scheduling decision. This is because the scheduling scheme produced by PSO-CS optimization does not depend on subjective considerations.
The total energy used by the sewage system for each day of operation, both before and after the optimization in that particular month, is displayed in Figure 11. The sewage system’s energy consumption ranges from 800 × 106 m5/d to 950 × 106 m5/d before optimization. Following optimization, the energy consumption remains within the same range, ranging from 650 × 106 m5/d to 750 × 106 m5/d, indicating a notable reduction in comparison to the pre-optimization levels. This demonstrates that the sewage system that has been adjusted may effectively reduce its energy usage while maintaining stable operation.

4. Conclusions

This paper establishes a sewage system scheduling optimization model by taking the actual production and operation of the sewage system of an oil extraction plant in the Daqing oilfield as an example. Aiming at the nonlinear, multi-constraint, and non-convex characteristics of the oilfield wastewater system scheduling model. The improved particle swarm optimization algorithm PSO-CS, which integrates the cuckoo optimization algorithm and particle swarm optimization algorithm, is proposed to solve this complex problem, and the optimal water quantity scheduling strategy of the sewage system is obtained. The main findings of this study are as follows:
(1)
The optimized sewerage system water supply direction and quantity were effectively adjusted. This indicates that the optimized water scheduling strategy can effectively improve the operation of the system;
(2)
The operational energy consumption of the sewage system was effectively reduced after optimization. The total operational energy consumption generated by the scheduling program on the same day was reduced from 879.95 × 106 m5/d to 712.84 × 106 m5/d, and the operational energy consumption of the sewage system after optimization was 19% lower than that of the system before optimization;
(3)
Significant increase in water output from sewage stations after optimization, with the number of sewage stations in and out of the water balance rising to 7 compared to 0 before optimization and the rate of meeting the water requirements of downstream water injection stations reaching 100%;
(4)
Monthly scheduling optimization of the sewage system has a significant reduction in energy consumption. Before optimization, the energy consumption of the sewage system in June was 800 × 106 m5/d–950 × 106 m5/d. After optimization, the energy consumption was 650 × 106 m5/d–750 × 106 m5/d, and the energy consumption of the sewage system operation after optimization was significantly reduced.
Although this research significantly improves the optimization of water scheduling in oilfield sewage systems by improving the PSO algorithm, there are still some shortcomings: the computational complexity of the improved algorithm may be higher when dealing with very large-scale systems, and it needs to be further optimized to improve the computational efficiency. Future research can further explore and develop more efficient optimization algorithms to cope with larger-scale and more complex scheduling problems.

Author Contributions

Conceptualization, S.C. and S.Z.; methodology, S.C.; software, S.Z.; validation, S.C., S.Z. and M.J.; formal analysis, Y.L.; investigation, S.Z.; resources, S.C.; data curation, J.X.; writing—original draft preparation, S.C.; writing—review and editing, S.Z.; visualization, B.G.; supervision, M.J.; project administration, J.X.; funding acquisition, S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China (Nos. 52104065, 52074090), China Postdoctoral Science Foundation (Nos. 2022T150089, 2020M681064), Natural Science Foundation of Heilongjiang Province (No. LH2021E019), Heilongjiang Postdoctoral Foundation (No. LBH-Z20101), Scientific Research Personnel Training Foundation of Northeast Petroleum University (No. XNYXLY202103) and Northeast Petroleum University Scientific Research Foundation (No. 2019KQ54).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

Author Yuchun Li was employed by the company Daging Oilfield Design Institute Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Qaderi, F. Efficiency Optimization of Petroleum Wastewater Treatment by Using of Serial Moving Bed Biofilm Reactors. J. Clean. Prod. 2018, 192, 665–677. [Google Scholar] [CrossRef]
  2. Farhadian, M.; Duchez, D.; Vachelard, C.; Larroche, C. Monoaromatics Removal from Polluted Water through Bioreactors A Review. Water Res. 2008, 42, 1325–1341. [Google Scholar] [CrossRef] [PubMed]
  3. Roudneshin, M.; Azadeh, A. A Novel Multi-objective Fuzzy Model for Optimization of Oil Sludge Management by Considering Health, Safety and Environment (HSE) and Resiliency Indicators in a Gas Refinery. J. Clean. Prod. 2018, 206, 559–571. [Google Scholar] [CrossRef]
  4. Zhou, X.; Zhang, H.; Qiu, R.; Lv, M.; Xiang, C.; Long, Y.; Liang, Y. A Two-stage Stochastic Programming Model for the Optimal Planning of a Coal-to-liquids Supply Chain Under Demand Uncertainty. J. Clean. Prod. 2018, 22, 10–28. [Google Scholar] [CrossRef]
  5. Smith, K.; Guo, S.; Zhu, Q.; Dong, X.; Liu, S. An Evaluation of the Environmental Benefit and Energy Footprint of China’s Stricter Wastewater Standards: Can Benefit be Increased? J. Clean. Prod. 2019, 219, 723–733. [Google Scholar] [CrossRef]
  6. Zhou, X.; Zhang, H.; Qiu, R.; Liang, Y.; Wu, G.; Xiang, C.; Yan, X. A hybrid time MILP model for the pump scheduling of multi-product pipelines based on the rigorous description of the pipeline hydraulic loss changes. Comput. Chem. Eng. 2019, 121, 174–199. [Google Scholar] [CrossRef]
  7. Jia, D.; Zhang, T.; Wu, L.; Su, X.; Bai, T.; Huang, Q. Multi-objective cooperative optimization of reservoir scheduling and water resources allocation for inter-basin water transfer project based on multi-stage coupling model. J. Hydrol. 2024, 630, 130673. [Google Scholar] [CrossRef]
  8. Yang, Y.; Chen, S.J.; Zhou, Y.R.; Ma, G.W.; Huang, W.B.; Zhu, Y.M. Method for quantitatively assessing the impact of an inter-basin water transfer project on ecological environment-power generation in a water supply region. J. Hydrol. 2023, 618, 129250. [Google Scholar] [CrossRef]
  9. Wu, L.Z.; Bai, T.; Huang, Q. Tradeoff analysis between economic and ecological benefits of the inter basin water transfer project under changing environment and its operation rules. J. Clean. Prod. 2020, 248, 119294. [Google Scholar] [CrossRef]
  10. Wang, B.; Liang, Y.; Yuan, M. Water transport system optimisation inoilfields: Environmental and economic benefits. J. Clean. Prod. 2019, 237, 117768. [Google Scholar] [CrossRef]
  11. Wang, Z.; Tian, J.; Feng, K. Optimal Allocation of Regional Water Resources Based on Simulated Annealing Particle Swarm Optimization Algorithm. Energy Rep. 2022, 8, 9119–9126. [Google Scholar] [CrossRef]
  12. Al-Abaji, M.A. Cuckoo Search Algorithm: Review and its Application. Tikrit J. Pure Sci. 2022, 26, 137–144. [Google Scholar] [CrossRef]
  13. Chen, H.; Wang, W.; Chen, X.; Qiu, L. Multi-objective Reservoir Operation Using Particle Swarm Optimization with Adaptive Random Inertia Weights. Water Sci. Eng. 2020, 13, 136–144. [Google Scholar] [CrossRef]
  14. Ezzeldin, R.; Zelenakova, M.; Abd-Elhamid, H.F.; Pietrucha-Urbanik, K.; Elabd, S. Hybrid Optimization Algorithms of Firefly with GA and PSO for the Optimal Design of Water Distribution Networks. Water 2023, 15, 1906. [Google Scholar] [CrossRef]
  15. Deng, A.N.; Wang, Z.C.; Liu, H.F.; Wu, T.H. A bio-inspired algorithm for a classical water resources allocation problem based on Adleman–Lipton model. Desalin. Water Treat. 2020, 185, 168–174. [Google Scholar] [CrossRef]
  16. de Oliveira Turci, L.; Sun, H.; Bai, M.; Wang, J.; Hu, P. Water Pump Station Scheduling Optimization Using an Improved Genetic Algorithm Approach. In Proceedings of the 2019 IEEE Congress on Evolutionary Computation (CEC), Wellington, New Zealand, 10–13 June 2019. [Google Scholar]
  17. Dong, W.; Yang, Q. Data-Driven Solution for Optimal Pumping Units Scheduling of Smart Water Conservancy. IEEE Internet Things J. 2019, 7, 1919–1926. [Google Scholar] [CrossRef]
  18. Eryiğit, M.; Sulaiman, S.O.; Najm, A.B.A.; Mhedi, N.M. Optimal management of multiple water resources by a heuristic optimization for a water supply in the desert cities of Western Iraq. Desalination Water Treat. 2023, 281, 7–14. [Google Scholar] [CrossRef]
  19. Ma, H.; Shen, S.; Yu, M.; Yang, Z.; Fei, M.; Zhou, H. Multi-population techniques in nature inspired optimization algorithms: A comprehensive survey. Swarm Evol. Comput. 2019, 44, 365–387. [Google Scholar] [CrossRef]
  20. Wu, X.; Shen, X.; Wei, C.; Xie, X.; Li, J. A Hybrid Particle Swarm Optimization-Genetic Algorithm for Multiobjective Reservoir Ecological Dispatching. Water Resour. Manag. 2024, 38, 2229–2249. [Google Scholar] [CrossRef]
  21. Kennedy, J.; Eberhart, R. Particle Swarm Optimization. In Proceedings of the ICNN’95-International Conference on Neural Networks, Perth, WA, USA, 27 November–1 December 1995; pp. 1942–1948. [Google Scholar]
  22. Kim, Y.-G.; Sun, B.-Q.; Kim, P.; Jo, M.-B.; Ri, T.-H.; Pak, G.-H. A Study on Optimal Operation of Gate-controlled Reservoir System for Flood Control Based on PSO Algorithm Combined with Rearrangement Method of Partial Solution Groups. J. Hydrol. 2021, 593, 125783. [Google Scholar] [CrossRef]
  23. Shi, Y.H.; Eberhart, R. Population Diversity of Particle Swarm. In Proceedings of the 2008 IEEE Congress on Evolutionary Computation (IEEE World Congress on Computational Intelligence), Hong Kong, China, 1–6 June 2008. [Google Scholar]
  24. Cuong-Le, T.; Minh, H.-L.; Khatir, S.; Wahab, M.A.; Tran, M.T.; Mirjalili, S. A Novel Version of Cuckoo Search Algorithm for Solving Optimization Problems. Expert Syst. Appl. 2021, 186, 115669. [Google Scholar] [CrossRef]
  25. Yang, X.-S.; Deb, S. Engineering Optimisation by Cuckoo Search. Math. Model. Numer. Optim. 2010, 1, 330–343. [Google Scholar] [CrossRef]
  26. Dunn, O.J. Multiple Comparisons among Means. J. Am. Stat. Assoc. 1961, 56, 52–64. [Google Scholar] [CrossRef]
Figure 1. Summary map of the sewage system in the study area.
Figure 1. Summary map of the sewage system in the study area.
Water 16 02623 g001
Figure 2. Optimization mechanism after fusion of PSO algorithm and CS algorithm.
Figure 2. Optimization mechanism after fusion of PSO algorithm and CS algorithm.
Water 16 02623 g002
Figure 3. Flowchart of hybrid PSO algorithm and CS algorithm.
Figure 3. Flowchart of hybrid PSO algorithm and CS algorithm.
Water 16 02623 g003
Figure 4. Convergence curves of the seven algorithms on the CEC2022 test function (20 dimensions).
Figure 4. Convergence curves of the seven algorithms on the CEC2022 test function (20 dimensions).
Water 16 02623 g004
Figure 5. Histogram of Bonferroni–Dunn detection results for PSO-CS algorithm and other algorithms based on the mean and standard deviation of the optimal values in Table 3. (a) Mean, (b) standard deviation.
Figure 5. Histogram of Bonferroni–Dunn detection results for PSO-CS algorithm and other algorithms based on the mean and standard deviation of the optimal values in Table 3. (a) Mean, (b) standard deviation.
Water 16 02623 g005
Figure 6. Water scheduling scheme of sewage system before optimization.
Figure 6. Water scheduling scheme of sewage system before optimization.
Water 16 02623 g006
Figure 7. Optimized water scheduling scheme for the sewage system.
Figure 7. Optimized water scheduling scheme for the sewage system.
Water 16 02623 g007
Figure 8. Comparison of system energy consumption before and after optimization.
Figure 8. Comparison of system energy consumption before and after optimization.
Water 16 02623 g008
Figure 9. Water supply of each injection station before and after optimization.
Figure 9. Water supply of each injection station before and after optimization.
Water 16 02623 g009
Figure 10. Number of injection stations meeting water requirements before and after optimization for the month.
Figure 10. Number of injection stations meeting water requirements before and after optimization for the month.
Water 16 02623 g010
Figure 11. System operation energy consumption before and after optimization.
Figure 11. System operation energy consumption before and after optimization.
Water 16 02623 g011
Table 1. Fluid state index value.
Table 1. Fluid state index value.
Flow Pattern k
Laminar flow1
Turbulent flowHydraulically smooth0.25
Mixed friction0.123
Completely rough0
Table 2. Performance comparison between PSO-CS and other algorithms on the CEC2022 test set.
Table 2. Performance comparison between PSO-CS and other algorithms on the CEC2022 test set.
FunctionDimValuesPSOPSO-CSWOASACSCOAHO
F 1 10Mean 4.5945 × 1032.9713 × 1026.0795 × 1032.7305 × 1043.0000 × 1023.0000 × 1029.0867 × 103
Std1.6852 × 1030.0000 × 1005.6259 × 1021.9669 × 1040.0000 × 1000.0000 × 1002.2343 × 103
Rank3146225
20Mean 2.4409 × 1043.0001 × 1021.3581 × 1038.6333 × 1043.0469 × 1021.5951 × 1054.9359 × 104
Std5.2925 × 1032.4612 × 10−24.8853 × 1021.8075 × 1041.3221 × 1003.0678 × 10−111.3501 × 104
Rank4136275
F 2 10Mean 4.7014 × 10 × 1024.0004 × 1024.3611 × 1025.1237 × 1024.0607 × 1025.0172 × 1047.2524 × 102
Std6.0995 × 103.2406 × 10−23.8456 × 101.7181 × 102.8569 × 1002.4964 × 1021.9308 × 102
Rank4135276
20Mean 7.9187 × 1024.0000 × 1024.5007 × 1029.6509 × 1024.1323 × 1027.2730 × 103 1.6328 × 103
Std7.7825 × 102.5139 × 103.3160 × 102.4981 × 1021.5138 × 100.0000 × 1003.9209 × 102
Rank4135276
F 3 10Mean 6.2905 × 1026.0313 × 1026.2463 × 102 6.7429 × 102 6.0417 × 102 6.9074 × 102 6.3250 × 102
Std1.7356 × 1008.7382 × 10−11.7667 × 102.6685 × 1000.0000 × 1000.0000 × 1002.0709 × 100
Rank4136275
20Mean 6.4835 × 1026.0321 × 1026.5840 × 102 7.1149 × 1026.1815 × 102 7.4688 × 1026.7073 × 102
Std5.3247 × 1004.8070 × 1001.4006 × 109.1742 × 1003.8591 × 1001.1984 × 10−137.8182 × 100
Rank3146275
F 4 10Mean 8.4823 × 1028.1559 × 108.4278 × 1029.1932 × 1028.1940 × 1029.0608 × 102 8.3458 × 102
Std7.9324 × 1002.7722 × 1001.9700 × 108.8567 × 1002.1106 × 1000.0000 × 1001.4224 × 101
Rank5147263
20Mean 9.5769 × 102 8.6786 × 1029.1931 × 1021.0924 × 1038.9122 × 1021.1106 × 103 9.5849 × 102
Std9.8838 × 1002.2900 × 1013.4061 × 1012.1452 × 1017.6080 × 1002.3967 × 10−131.7615 × 101
Rank4136275
F 5 10Mean 1.1892 × 103 9.0531 × 1021.3677 × 1031.9351 × 103 9.0638 × 1023.8096 × 103 1.1502 × 103
Std5.8726 × 1010.0000 × 1003.7566 × 1015.3324 × 1022.8555 × 1000.0000 × 1002.2494 × 102
Rank4156273
20Mean 4.3709 × 103 1.4389 × 1033.5909 × 103 8.3292 × 1031.4948 × 103 5.5721 × 103 3.0879 × 103
Std5.8027 × 1023.3996 × 1021.2000 × 1031.4561 × 1036.9299 × 1019.5869 × 10−135.1781 × 102
Rank5147263
F 6 10Mean 9.0851 × 105 1.8008 × 103 5.0653 × 103 6.6556 × 103 6.9061 × 103 2.1947 × 108 2.2076 × 106
Std4.0209 × 105 3.3110 × 10−19.7850 × 1021.9873 × 1037.3874 × 1010.0000 × 1003.1119 × 106
Rank5123476
20Mean 1.1266 × 108 1.8440 × 103 8.1254 × 103 1.1053 × 105 1.0345 × 104 1.0197 × 1010 6.2254 × 108
Std8.1138 × 10 74.0885 × 1007.9809 × 1032.7820 × 1059.5156 × 1032.0105 × 1062.9250 × 108
Rank6124357
F 7 10Mean 2.0510 × 103 2.0272 × 103 2.0684 × 1032.1363 × 103 2.0298 × 103 2.2848 × 103 2.0770 × 103
Std4.2911 × 1001.3434 × 1001.2496 × 1011.3888 × 1011.4291 × 1002.0497 × 1014.0105 × 101
Rank3145267
20Mean 2.1344 × 1032.0698 × 103 2.1733 × 103 2.3343 × 103 2.0885 × 103 2.5103 × 103 2.1858 × 103
Std1.2132 × 1012.9131 × 1014.3172 × 1014.2579 × 1015.4936 × 1004.7935 × 1014.3149 × 101
Rank3146275
F 8 10Mean 2.2295 × 103 2.2114 × 103 2.2293 × 103 2.2459 × 103 2.2209 × 103 2.9437 × 103 2.2351 × 103
Std1.8604 × 1001.0935 × 1005.1993 × 1006.3473 × 1011.2369 × 1001.3207 × 1028.0561 × 100
Rank4136275
20Mean 2.3684 × 103 2.2339 × 103 2.3328 × 103 2.3181 × 103 2.2519 × 1033.3954 × 104 2.4012 × 103
Std1.5466 × 1011.1876 × 1001.2793 × 1021.9100 × 1017.8982 × 10−12.5153 × 1041.7064 × 102
Rank6154372
F 9 10Mean 2.6231 × 1032.5273 × 103 2.5294 × 103 2.7251 × 103 2.5289 × 103 4.0536 × 1032.6748 × 103
Std1.1439 × 1021.1962 × 10−11.3342 × 10−21.0071 × 1027.4448 × 10−22.2546 × 1024.7994 × 101
Rank4126375
20Mean 2.5549 × 103 2.4808 × 103 2.4829 × 103 3.0348 × 103 2.4809 × 103 5.1532 × 103 2.9109 × 103
Std4.6963 × 1011.0868 × 10−32.2575 × 1001.6038 × 1022.3199 × 10−26.9633 × 1025.9396 × 101
Rank4136275
F 10 10Mean 2.5078 × 103 2.5004 × 103 2.5688 × 103 2.6107 × 1032.5009 × 103 4.9524 × 103 2.5165 × 103
Std1.1964 × 1001.0658 × 10−29.6564 × 10 1.0047 × 1027.2031 × 10−23.9130 × 10 1.1755 × 101
Rank3156274
20Mean 5.5576 × 1032.5008 × 103 5.2629 × 1036.2856 × 103 2.5009 × 103 8.7814 × 103 5.1262 × 103
Std2.5691 × 1021.0372 × 10−15.1585 × 1021.1019 × 1033.7636 × 10 5.9494 × 10 1.7837 × 103
Rank6134275
F 11 10Mean 2.7946 × 103 2.6000 × 103 2.8753 × 103 3.1918 × 103 2.7504 × 103 5.1052 × 1032.8180 × 103
Std1.4312 × 10 2.8584 × 10−91.7655 × 1024.9260 × 1023.5493 × 10−90.0000 × 1004.2883 × 10
Rank3156274
20Mean 5.1297 × 103 2.9251 × 103 3.0096 × 103 9.5241 × 103 2.9500 × 103 1.0783 × 104 6.6099 × 103
Std1.6095 × 1031.2504 × 1021.0273 × 1021.1613 × 1037.0711 × 10 0.0000 × 1004.1787 × 102
Rank4136275
F 12 10Mean 2.8651 × 103 2.8589 × 1032.8696 × 1032.9014 × 103 2.8626 × 103 2.9907 × 103 2.9887 × 103
Std1.6179 × 10−11.6050 × 10−18.0333 × 1004.9578 × 1005.2891 × 10−14.3252 × 10 4.1039 × 10
Rank3145276
20Mean 2.9702 × 103 2.9474 × 103 3.0309 × 103 3.2915 × 103 2.9510 × 103 3.5523 × 1033.3798 × 103
Std2.0820 × 10 7.9699 × 10−11.5233 × 10 9.7613 × 10 7.6760 × 1002.3911 × 10 2.5929 × 102
Rank3146275
Table 3. Friedman’s test results for PSO-CS and other algorithms based on solving for the mean and standard deviation of the optimal values in Table 2 (best rankings are in bold).
Table 3. Friedman’s test results for PSO-CS and other algorithms based on solving for the mean and standard deviation of the optimal values in Table 2 (best rankings are in bold).
Mean of Optimal ValuesSD
N1212
Test resultsChi-square57.35859.964
p-value<0.01<0.01
PSO3.834.33
PSO-CS1.001.00
WOA3.583.42
Friedman test valueSA5.715.50
CS2.292.17
COA6.426.75
HO5.174.83
Table 4. Sewage system scheduling optimization scheme.
Table 4. Sewage system scheduling optimization scheme.
Name of Sewage StationName of Water Injection StationPipeline PathLength (m)
NIII;-1DZN3NIII;-1D→ZN32600
NIII;-1DZN11NIII;-1D→ZN113100
NIII;-1DZN16-PNIII;-1D→ZN16-P1100
NIII;-1DZN17-PNIII;-1D→ZN17-P1200
NIII;-1ZN8NIII;-1→ZN81207
NIII;-1ZN7NIII;-1→ZN72250
NIII;-1ZN19NIII;-1→ZN19500
NIII;-1ZN10NIII;-1→ZN101382
NII-2-GZNII-4NII-2-G→ZNII-45150
NII-2-GZN22NII-2-G→ZN221350
NII-2ZNII-4NII-2→ZNII-45240
NII-2ZN8NII-2→ZN8650
NII-2ZN22NII-2→ZN22340
NII-2ZN11NII-2→ZN112390
N16ZNII-4N16→ZNII-42550
N16ZN22N16→ZN22920
N16ZN7N16→ZN7800
N16ZN11N16→ZN111100
N16ZN16-PN16→ZN16-P70
N11-DZN11N11-D→ZN1150
N11-DZN16-PN11-D→ZN16-P820
NII-7 TZNII-4NII-7 T→ZNII-42750
NII-7 TZN22NII-7 T→ZN221925
NII-7 TZN16-PNII-7 T→ZN16-P1100
NII-7 DZN16-PNII-7 D→ZN16-P1000
NII-7 DZN17-PNII-7 D→ZN17-P980
NIII;-2DZNIII;-2NIII;-2D→ZNIII;-260
NIII;-2DZN16-PNIII;-2D→ZN16-P3500
SN2801ZN23SN2801→ZN232780
SN2801ZN9SN2801→ZN9600
N20-DZN20N20-D→ZN201378
N4ZN20N4→ZN20850
N4ZN7N4→ZN7230
N17ZNIII;-2N17→ZNIII;-21750
N17ZN8N17→ZN8800
N17ZN7N17→ZN7450
N17ZN10N17→ZN101200
N17ZN17-PN17→ZN17-P1350
NIII;-7TZN8NIII;-7T→ZN81350
NIII;-7TZN7NIII;-7T→ZN72110
NIII;-7TZN17-PNIII;-7T→ZN17-P800
N5-PZN5N5-P→ZN5780
N5ZN7N5→ZN72556
N5ZN10N5→ZN104200
N5ZN5N5→ZN51378
NII-1ZNII-1NII-1→ZNII-1190
NII-1ZNII-4NII-1→ZNII-43125
NII-1ZN23NII-1→ZN23436
SN25-DZNIII;-2SN25-D→ZNIII;-2450
SN25-DZN20SN25-D→ZN205000
SN25-DZN17-PSN25-D→ZN17-P500
N12-DZN12N12-D→ZN12140
N12-DZN13N12-D→ZN133360
N12-DZN12-PN12-D→ZN12-P1340
N15ZNII-1N15→ZNII-11632
N15ZN8N15→ZN8896
N15ZN9N15→ZN91750
N15ZN3N15→ZN33840
N15ZN14N15→ZN142000
N15ZN12-PN15→ZN12-P530
NIII;-6TZN8NIII;-6T→ZN81280
NIII;-6TZN3NIII;-6T→ZN31060
N6-PZN6-PN6-P→ZN6-P580
N6-PZN6-WN6-P→ZN6-W780
N6-PZN19N6-P→ZN193136
N3ZN6-WN3→ZN6-W3206
N3ZN9N3→ZN92750
N3ZN3N3→ZN3300
N3ZN10N3→ZN102800
N6-DZN6-WN6-D→ZN6-W2340
N6-DZN19N6-D→ZN19575
N6ZN6-PN6→ZN6-P1720
N6ZN6-WN6→ZN6-W630
N13-PZNII-1N13-P→ZNII-16890
N13-PZNII-4N13-P→ZNII-46920
N13-PZN23N13-P→ZN23579
N13-PZN13N13-P→ZN13950
N13-PZN14N13-P→ZN14920
N13-NDZN12N13-ND→ZN123565
N13-NDZN13N13-ND→ZN131300
N13-DZN23N13-D→ZN232205
N13-DZN13N13-D→ZN13380
Table 5. Water requirements for each injection station for the day.
Table 5. Water requirements for each injection station for the day.
NumberNameWater Demand
(m3/d)
1ZNII-18352
2ZNII-40
3ZNIII;-29997
4ZN814,795
5ZN229795
6ZN239548
7ZN2016,415
8ZN97109
9ZN6-P0
10ZN6-W7908
11ZN714,563
12ZN39008
13ZN127068
14ZN195930
15ZN1316,825
16ZN147628
17ZN119259
18ZN107321
19ZN57558
20ZN12-P11,456
21ZN16-P5276
22ZN17-P4383
Table 6. Volume of water entering and leaving the sewage station.
Table 6. Volume of water entering and leaving the sewage station.
NumberNameThe Volume of Incoming Water for the Day (m3)Water Output for the Day (m3)Optimized Water Output (m3)
1NII-110,83348218352
2NII-211,8475545.559795
3NII-7D711022495582
4NII-7T967300
5NIII-1D15,28700
6NIII-1655954886559
7NIII-2D420138073201
8NIII-6T616700
9NIII-7T16,5813552.0516,581
10N20-D376132262021
11N6-P984400
12N610,37207908
13N6-D778411205930
14N310,92794339008
15N12-D956638567068
16N1722,76812,652.420,911
17N13-D698827916988
18N13-P18,462863518,462
19N13-ND761152785841
20N1514,68814,08013,688
21N11-D85633881.36563
22N414,394753214,394
23N5-P728120007281
24N596472000277
25N16975582422755
26SN25-D244217461442
27SN2801758732137587
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, S.; Zhou, S.; Li, Y.; Jiang, M.; Guan, B.; Xi, J. Optimization Method for Digital Scheduling of Oilfield Sewage System. Water 2024, 16, 2623. https://doi.org/10.3390/w16182623

AMA Style

Chen S, Zhou S, Li Y, Jiang M, Guan B, Xi J. Optimization Method for Digital Scheduling of Oilfield Sewage System. Water. 2024; 16(18):2623. https://doi.org/10.3390/w16182623

Chicago/Turabian Style

Chen, Shuangqing, Shun Zhou, Yuchun Li, Minghu Jiang, Bing Guan, and Jiahao Xi. 2024. "Optimization Method for Digital Scheduling of Oilfield Sewage System" Water 16, no. 18: 2623. https://doi.org/10.3390/w16182623

APA Style

Chen, S., Zhou, S., Li, Y., Jiang, M., Guan, B., & Xi, J. (2024). Optimization Method for Digital Scheduling of Oilfield Sewage System. Water, 16(18), 2623. https://doi.org/10.3390/w16182623

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop