# Parameter Identification of Multistage Fracturing Horizontal Well Based on PSO-RBF Neural Network

Authors: Rongwang Yin; Qingyu Li; Peichao Li; Detang Lu
Journal: Scientific Programming
 (2020)

---

## Abstract

In order to more accurately identify multistage fracturing horizontal well (MFHW) parameters and address the heterogeneity of reservoirs and the randomness of well-production data, a new method based on the PSO-RBF neural network model is proposed. First, the GPU parallel program is used to calculate the bottomhole pressure of a multistage fracturing horizontal well. Second, most of the above pressure data are imported into the RBF neural network model for training. In the training process, the optimization function of the global optimal solution of the PSO algorithm is employed to optimize the parameters of the RBF neural network, and eventually, the required PSO-RBF neural network model is established. Third, the resulting neural network is tested using the remaining data. Finally, a field case of a multistage fracturing horizontal well is studied by using the presented PSO-RBF neural network model. The results show that in most cases, the proposed model performs better than other models, with the highest correlation coefficient, the lowest mean, and absolute error. This proves that the PSO-RBF neural network model can be applied effectively to horizontal well parameter identification. The proposed model has great potential to improve the prediction accuracy of reservoir physical parameters.


---

## Body

## 1. Introduction

Reservoir description is an important means of estimating the physical properties of a reservoir in the petroleum industry. These reservoir physical parameters require better methods to improve the prediction accuracy to further enhance the subsequent success rate of exploration and development. Multiphase flow in oil and gas fields usually refers to the simultaneous flow of multiple fluids in a reservoir, which generally occurs during oil and gas production [1, 2]. For the old oilfields, oil and gas development has experienced stages such as self-spraying, secondary, and tertiary oil recovery, which makes the process of the oil and gas production complex and challenging. In order to deal with this complexity, a lot of studies have been done at home and abroad, and the reservoir numerical simulations have been introduced to estimate the multiphase flow characteristics of oil, natural gas, and water [3, 4].In the petroleum industry, well-logging data have been commonly used to predict formation parameters such as permeability, formation boundaries, and fracture length. Permeability is defined as a measure of the ability of a porous medium to allow fluid to pass through it. The concept of permeability is important for determining accurate reservoir description, simulation, and management. Therefore, prior to any modeling or calculation, the permeability of porous media must be determined. The earliest method of permeability prediction is the empirical correlation between permeability and other petrophysical properties such as porosity and water saturation. These correlations have made some successes in sandstone reservoirs. However, it is not applicable to heterogeneous formations.In recent years, the artificial neural network (ANN) has been applied to the reservoir field, solving many highly complex nonlinear problems [5]. The ANN is considered to be a nonlinear tool that can predict complex nonlinear system behavior well. It is also used to solve many different types of problems in petroleum engineering, such as reservoir characterization, permeability prediction, vertical multiphase flow downhole pressure prediction, and dissolved gas driven flow prediction [6–8].The ANN is considered to be a method that can effectively replace traditional techniques of predicting permeability from well-logging data. The disadvantage of the ANN application is the requirement to identify the structure of the network. In order to overcome the shortcoming, in this study, an MFHW parameter identification model based on the PSO-RBF (particle swarm optimization, PSO; radial basis function, RBF) neural network is proposed. The main purpose is to seek a better prediction of reservoir properties. The parameters of the radial basis neural network are designed by using the particle swarm algorithm to obtain higher prediction accuracy. The parameter identification model based on the PSO-RBF neural network has high learning ability and flexibility when encountering complex problems, which makes the proposed model applicable to most engineering practices.


## 2. Review of Related Works

In the past few decades, researchers have done a lot of work to study reservoir parameter estimation problems, such as permeability, porosity, and fracture length. Taking the estimation of permeability as an example, there are currently three main methods: core analysis, well test analysis, and artificial intelligence methods.The most accurate method of permeability prediction is to perform core analysis in the laboratory by applying Darcy’s law. Among the different techniques, the permeability obtained from the core analysis is more efficient than that obtained by other methods, and the permeability from the core analysis can be used to verify other estimation models. However, the core analysis process is time consuming and difficult to use widely.Another method of obtaining reservoir permeability is well test analysis. The analytical solution is solved by using the porous flow equation, and some rock properties such as permeability and porosity can be then obtained inversely from the pressure curve. The data available after well testing help petroleum engineers calculate formation permeability, skin factor, and wellbore storage [9]. Due to the cost of testing, well testing is not a good solution for any reservoirs that have been developed.ANNs are increasingly being employed to predict reservoir properties based on geophysical well-logging data [10, 11]. Mohaghegh et al. pointed out that neural networks are powerful tools for identifying the relationship between permeability and well-logging data [12]. In addition, Aminian et al.’s research shows that artificial neural networks can be used to predict formation permeability, even in highly heterogeneous reservoirs [13, 14]. Based on the artificial neural network and particle swarm optimization algorithm, Ahmadi et al. proposed a method for predicting horizontal well productivity in a pseudosteady state, and the experimental results achieved a high fitting accuracy (<0.82%) [15]. Zhang et al. proposed a hybrid algorithm that combines a particle swarm optimization algorithm with a back-propagation algorithm to train the weights of a feed-forward neural network. The results show that the hybrid algorithm is superior to the adaptive particle swarm optimization algorithm and backpropagation neural network (BPNN) in both convergence speed and convergence accuracy [16]. Based on the RBF and BPNN, Chen developed two water treatment models. The model has been trained and checked separately from actual data from the water plant. Compared with the model based on the BPNN, the model based on the RBF neural network has the characteristics of better approximation and faster convergence [17]. Based on the models of the BPNN and RBFNN (the radial basis function neural network, RBFNN), Yang et al. proposed a model for cement decomposition and calcination process. The results show that the model based on the RBFNN can achieve very high fitting results [18].Saemi et al. improved the predictive performance of neural networks by using the GA to optimize the network parameters of the ANN [19]. In addition, Kaydani et al.’s findings indicate that the GA and ANN network structures designed by subregion can predict the permeability of a heterogeneous reservoir in Iran [20]. Tahmasebi et al. proposed four methods for predicting permeability of different neural network structures and statistically compared the results obtained. Finally, a new method of the modular neural network (MNN) was obtained [21]. By combining cuckoo optimization algorithm (COA), particle swarm optimization algorithm (PSO), imperialist competitive algorithm (ICA), and Levenberg–Marquardt algorithm (LM), Kaydan et al. proposed a new method to estimate the permeability [22].Rough set theory can be successfully applied to the prediction of permeability of porous media [23–28]. By establishing a rough set model, Ilkhchi et al. successfully predicted the permeability of an offshore gas field in Iran [29].Support vector regression (SVR) based on the principle of structural risk minimization is a promising machine learning method. Through experimental comparisons, Gholami et al. found that the SVR method has faster speed and higher accuracy in predicting reservoir permeability [30].


## 3. Problem Formulation

### 3.1. Optimized Design

After fracturing tight sandstone reservoirs, the physical properties of the reservoirs are characterized by extremely low porosity, a wide range of permeability changes, and complex pore-permeability relationships. Therefore, conventional production prediction models often fail to meet the requirements in terms of prediction accuracy.Effective algorithms for complex nonlinear problems are essentially optimization problems. The so-called optimization refers to the problem of seeking the maximum or minimum value of the given objective function according to the change of design variables under the condition of meeting some certain constraints. It can be described by the following equation:(1)minfxormaxfx,s.t.gix=0,i=1,2,…,hjx≥0orhjx≤0,j=1,2,…,where x∈S and S is the design parameter space that meets the limiting conditions, called the solution space. fx=fx1,x2,…,xn is the objective function to be optimized. x=x1,x2,…,xnT, x1,x2,…,xn is the optimization design variable, and minfxormaxfx is the minimum or maximum value of the objective function. gix=0,i=1,2,…, are the equality constraints, and hjx≥0orhjx≤0 denotes the inequality constraints.Taking MFHW as the studied object, the PSO-RBF neural network was used to predict permeability and fracture half-length. The optimization objective function is the minimum RMSE (root mean square error) value of permeability and half-length of fracture. The input variables of the network are the formation pressure, the number of fractures, stratum thickness, the well storage constant, the well storage skin, and so on, and the output variables are the RMSE values.


### 3.2. Data Set

The data set describes the characteristics and behavior of the reservoir’s input and output parameters and is the training data set for the reservoir model. It is obtained using static and dynamic data such as porosity, permeability, pressure, and reservoir production values. For a given source unit, according to the line source superposition principle, the pressure distribution at any point (xD, yD, and zD) in the formation space shall be(2)pDijxD,yD,zD,tD=∫0tDqDijτGxDxD,τDGyDyD,τDGzDzD,τDdτ.The coefficient matrix is expressed as follows:(3)M11M12⋯M1NMi−1M21M22⋯M2NMi−1⋮⋮⋮⋮⋮MNMi1MNMi2⋯MNMiNMi⋮11⋯10⋅qD11qD12⋮qDNMipDwf=00⋮01.As for the definitions of various symbols appearing in equations (2) and (3), the reader can refer to our previous paper [31]. Note that, the formation pressure data of MFHW are obtained by the GPU parallel calculation method proposed in the reference [31]. All training data sets are divided into two parts before the PSO-RBF training begins. One is the input data set, and the other is the output data set. The input and output data sets are normalized within a specific range. Before the training begins, a standard normalization function is used to limit the range of the input and output data sets between −1 and 1, and the mathematical function is given in the following equation:(4)y=e2x−1e2x+1.During PSO-RBF training, the training data set is divided into three parts: training, validation, and testing. Training data are used in the training, and validation data are also used during training. However, it should be mentioned that validation data are used to check the network learning instead of training the network.


## 3.1. Optimized Design

After fracturing tight sandstone reservoirs, the physical properties of the reservoirs are characterized by extremely low porosity, a wide range of permeability changes, and complex pore-permeability relationships. Therefore, conventional production prediction models often fail to meet the requirements in terms of prediction accuracy.Effective algorithms for complex nonlinear problems are essentially optimization problems. The so-called optimization refers to the problem of seeking the maximum or minimum value of the given objective function according to the change of design variables under the condition of meeting some certain constraints. It can be described by the following equation:(1)minfxormaxfx,s.t.gix=0,i=1,2,…,hjx≥0orhjx≤0,j=1,2,…,where x∈S and S is the design parameter space that meets the limiting conditions, called the solution space. fx=fx1,x2,…,xn is the objective function to be optimized. x=x1,x2,…,xnT, x1,x2,…,xn is the optimization design variable, and minfxormaxfx is the minimum or maximum value of the objective function. gix=0,i=1,2,…, are the equality constraints, and hjx≥0orhjx≤0 denotes the inequality constraints.Taking MFHW as the studied object, the PSO-RBF neural network was used to predict permeability and fracture half-length. The optimization objective function is the minimum RMSE (root mean square error) value of permeability and half-length of fracture. The input variables of the network are the formation pressure, the number of fractures, stratum thickness, the well storage constant, the well storage skin, and so on, and the output variables are the RMSE values.


## 3.2. Data Set

The data set describes the characteristics and behavior of the reservoir’s input and output parameters and is the training data set for the reservoir model. It is obtained using static and dynamic data such as porosity, permeability, pressure, and reservoir production values. For a given source unit, according to the line source superposition principle, the pressure distribution at any point (xD, yD, and zD) in the formation space shall be(2)pDijxD,yD,zD,tD=∫0tDqDijτGxDxD,τDGyDyD,τDGzDzD,τDdτ.The coefficient matrix is expressed as follows:(3)M11M12⋯M1NMi−1M21M22⋯M2NMi−1⋮⋮⋮⋮⋮MNMi1MNMi2⋯MNMiNMi⋮11⋯10⋅qD11qD12⋮qDNMipDwf=00⋮01.As for the definitions of various symbols appearing in equations (2) and (3), the reader can refer to our previous paper [31]. Note that, the formation pressure data of MFHW are obtained by the GPU parallel calculation method proposed in the reference [31]. All training data sets are divided into two parts before the PSO-RBF training begins. One is the input data set, and the other is the output data set. The input and output data sets are normalized within a specific range. Before the training begins, a standard normalization function is used to limit the range of the input and output data sets between −1 and 1, and the mathematical function is given in the following equation:(4)y=e2x−1e2x+1.During PSO-RBF training, the training data set is divided into three parts: training, validation, and testing. Training data are used in the training, and validation data are also used during training. However, it should be mentioned that validation data are used to check the network learning instead of training the network.


## 4. PSO-RBF Model

The advantage of ANNs over other conventional techniques is the ability to perform complex and highly nonlinear tasks accurately and swiftly. In most previous works related to reservoir models, the researchers used the BPNN to construct reservoir models. However, the BPNN has a problem of falling into local minimum during training time. In most cases, the network does not reach the global minimum to find the minimum error value. The global optimization ability of the PSO algorithm makes the radial basis neural network model based on the PSO algorithm have no local minimum problem.


### 4.1. PSO Algorithm

Inspired by the movement of the birds, the PSO algorithm was developed by Kennedy [32]. In this approach, every possible solution is considered as a particle. Each particle is characterized by its position and velocity. The position of the particle is defined in the hyperspace, and its dimension is equal to the number of nonoptimized parameters. Several particles are initially defined in hyperspace, which iteratively change their positions to determine the best position. The fitness of a particle is determined by a fitness function such as RMSE. This algorithm is similar to how birds search for food in a wide area.During the iterative execution of the algorithm, the values ofpbest and gbest are constantly updated. pbest is defined as the optimal position of the particle in the hyperspace, determined by the fitness value. gbest is the overall best position for all particles. In each iteration step, the speed is updated first and then the location is updated. The particles are accelerated to pbest and gbest by updating the speed:(5)vk+1i=wi⋅vki+c1⋅r1⋅pbesti−pki+c2⋅r2⋅gbest−pki,where vk+1i is the speed of the next iteration, wi is the inertia weight, vki is the current speed, r1 and r2 are random numbers, and pki is the current position of the particle. Update the position of the particle using the following equation:(6)pk+1i=pki+vk+1i⋅Δt.The initial position and velocity of each particle are randomly distributed. After initializing the position and velocity of all particles, the fitness is calculated. In the subsequent steps, the position and velocity are iteratively updated by the local best parameters and the global best parameters, as shown in Figure1.Figure 1
Flowchart of PSO algorithm.The entire flowchart can be divided into four parts, namely initialization, fitness calculation, status check, and updating of speed and position.


### 4.2. Radial Basis Function Neural Network

RBFNN has a small number of hidden layer neurons and are used to develop networks with good generalization capabilities [33]. RBFNN is considered to be a special type of artificial neural network because its architecture requires only a hidden layer, which allows the input space to be represented in new spaces of different hidden layer neurons. The structure is shown in Figure 2. In the training process of RBFNN, it appears as a linear model because all hidden neuron centers and calculations are fixed. The RBFNN hidden layer neurons perform a nonlinear transformation and map all inputs into the new input space. The output layer is considered as a linear converter and is applied to the new input space. The performance of RBFNN can be determined by adjusting the center (width) of the hidden layer neurons, and there is no specific formula to choose the width of the radial basis function [34]. RBFNN has been widely used in system prediction, pattern recognition, speech recognition, and adaptive control. It has also been used to solve the problems of oil and gas fields, e.g., oil and gas ratios of reservoirs, electromagnetics, resistivity, log data inversion, log data prediction, seismic properties, and nonlinear relationships between reservoir properties and seismic properties [35–37].Figure 2
RBFNN structure.The growth and overall structure of RBFNN is affected by the RBF. The RBFNN input is directly connected to each of the basic functions and produces an output. That is,(7)ϕi=exp−x−uTx−uσ2,where x represents the input data of the network, u is the center of the radial basis function, and σ is the width of the RBF σ>0. First, the hidden layer neurons are calculated based on the radius of the radial basis function and then transmitted to the output layer. Here, the sum of the products between the hidden layer neurons and the weight vectors is calculated to obtain the final network output yn; that is,(8)yn=∑m=1mwi⋅ϕi.


### 4.3. PSO-RBF Neural Network Method

The prediction accuracy of the RBFNN mainly depends on the central vectoru of the radial basis kernel function, the radius σ of the radial basis, and the connection weight wi between the output layer and the hidden layer. The traditional RBF neural network uses local information based on parameter space to set parameters, which results in the values of u, σ, and wi being local optimal solutions rather than global optimal solutions. In view of these defects in the RBFNN, this work employs the PSO algorithm to optimize the traditional RBFNN when identifying MFHW parameters. The RBFNN parameters optimized by the PSO algorithm are global optimal parameters, thus avoiding the problem of low reliability of RBFNN learning. The specific flowchart of the PSO-RBF algorithm is shown in Figure 3.Figure 3
PSO-RBF-based algorithm flowchart.PSO optimization RBFNN parameters are divided into two steps: the first step is to determine the center value and width of the basis function, and the second step is to determine the connection weight between the hidden layer and the input layer. In the optimization process, the data obtained by GPU parallel computing are used to train and verify the network.The algorithm flow of the first step is as follows.(1) Collection of reservoir model data samples; (2) cluster analysis of samples to generate centeru and width σ of basis functions; (3) initialization of the particle swarm algorithm using u and σ; (4) calculation of each particle root mean square error (RMSE); (5) updating local optimal solution and global optimal solution; (6) updating particle position and velocity; (7) repeating steps (4–6) until the accuracy requirement or number of iterations is reached; and (8) obtaining the center u and width σ of the basis function.The algorithm flow of the second step is as follows.(1) Calculation of the output of the hidden layer; (2) initialization of the weightwi and reinitialization of the particle group; (3) calculating the cumulative error of each particle; (4) updating the local optimal solution and the global optimal solution; (5) adjusting the position and velocity of the particles; (6) repeating steps (3–5); and (7) obtaining the parameter wi of the RBF neural network.


## 4.1. PSO Algorithm

Inspired by the movement of the birds, the PSO algorithm was developed by Kennedy [32]. In this approach, every possible solution is considered as a particle. Each particle is characterized by its position and velocity. The position of the particle is defined in the hyperspace, and its dimension is equal to the number of nonoptimized parameters. Several particles are initially defined in hyperspace, which iteratively change their positions to determine the best position. The fitness of a particle is determined by a fitness function such as RMSE. This algorithm is similar to how birds search for food in a wide area.During the iterative execution of the algorithm, the values ofpbest and gbest are constantly updated. pbest is defined as the optimal position of the particle in the hyperspace, determined by the fitness value. gbest is the overall best position for all particles. In each iteration step, the speed is updated first and then the location is updated. The particles are accelerated to pbest and gbest by updating the speed:(5)vk+1i=wi⋅vki+c1⋅r1⋅pbesti−pki+c2⋅r2⋅gbest−pki,where vk+1i is the speed of the next iteration, wi is the inertia weight, vki is the current speed, r1 and r2 are random numbers, and pki is the current position of the particle. Update the position of the particle using the following equation:(6)pk+1i=pki+vk+1i⋅Δt.The initial position and velocity of each particle are randomly distributed. After initializing the position and velocity of all particles, the fitness is calculated. In the subsequent steps, the position and velocity are iteratively updated by the local best parameters and the global best parameters, as shown in Figure1.Figure 1
Flowchart of PSO algorithm.The entire flowchart can be divided into four parts, namely initialization, fitness calculation, status check, and updating of speed and position.


## 4.2. Radial Basis Function Neural Network

RBFNN has a small number of hidden layer neurons and are used to develop networks with good generalization capabilities [33]. RBFNN is considered to be a special type of artificial neural network because its architecture requires only a hidden layer, which allows the input space to be represented in new spaces of different hidden layer neurons. The structure is shown in Figure 2. In the training process of RBFNN, it appears as a linear model because all hidden neuron centers and calculations are fixed. The RBFNN hidden layer neurons perform a nonlinear transformation and map all inputs into the new input space. The output layer is considered as a linear converter and is applied to the new input space. The performance of RBFNN can be determined by adjusting the center (width) of the hidden layer neurons, and there is no specific formula to choose the width of the radial basis function [34]. RBFNN has been widely used in system prediction, pattern recognition, speech recognition, and adaptive control. It has also been used to solve the problems of oil and gas fields, e.g., oil and gas ratios of reservoirs, electromagnetics, resistivity, log data inversion, log data prediction, seismic properties, and nonlinear relationships between reservoir properties and seismic properties [35–37].Figure 2
RBFNN structure.The growth and overall structure of RBFNN is affected by the RBF. The RBFNN input is directly connected to each of the basic functions and produces an output. That is,(7)ϕi=exp−x−uTx−uσ2,where x represents the input data of the network, u is the center of the radial basis function, and σ is the width of the RBF σ>0. First, the hidden layer neurons are calculated based on the radius of the radial basis function and then transmitted to the output layer. Here, the sum of the products between the hidden layer neurons and the weight vectors is calculated to obtain the final network output yn; that is,(8)yn=∑m=1mwi⋅ϕi.


## 4.3. PSO-RBF Neural Network Method

The prediction accuracy of the RBFNN mainly depends on the central vectoru of the radial basis kernel function, the radius σ of the radial basis, and the connection weight wi between the output layer and the hidden layer. The traditional RBF neural network uses local information based on parameter space to set parameters, which results in the values of u, σ, and wi being local optimal solutions rather than global optimal solutions. In view of these defects in the RBFNN, this work employs the PSO algorithm to optimize the traditional RBFNN when identifying MFHW parameters. The RBFNN parameters optimized by the PSO algorithm are global optimal parameters, thus avoiding the problem of low reliability of RBFNN learning. The specific flowchart of the PSO-RBF algorithm is shown in Figure 3.Figure 3
PSO-RBF-based algorithm flowchart.PSO optimization RBFNN parameters are divided into two steps: the first step is to determine the center value and width of the basis function, and the second step is to determine the connection weight between the hidden layer and the input layer. In the optimization process, the data obtained by GPU parallel computing are used to train and verify the network.The algorithm flow of the first step is as follows.(1) Collection of reservoir model data samples; (2) cluster analysis of samples to generate centeru and width σ of basis functions; (3) initialization of the particle swarm algorithm using u and σ; (4) calculation of each particle root mean square error (RMSE); (5) updating local optimal solution and global optimal solution; (6) updating particle position and velocity; (7) repeating steps (4–6) until the accuracy requirement or number of iterations is reached; and (8) obtaining the center u and width σ of the basis function.The algorithm flow of the second step is as follows.(1) Calculation of the output of the hidden layer; (2) initialization of the weightwi and reinitialization of the particle group; (3) calculating the cumulative error of each particle; (4) updating the local optimal solution and the global optimal solution; (5) adjusting the position and velocity of the particles; (6) repeating steps (3–5); and (7) obtaining the parameter wi of the RBF neural network.


## 5. A Field Case Study

### 5.1. Computation

The bottomhole flow pressure is analyzed and calculated by using the measured data of multistage fracturing horizontal wells in an oilfield in western China. The three-dimensional PEBI grid of fracturing of underground horizontal wells is shown in Figure4. The reservoir stratum is homogeneous, the fluid flow satisfies Darcy’s law, the horizontal and vertical permeability of the fracture are different, and the horizontal well is fractured in a closed rectangular reservoir. Using the GPU-based MFHW bottomhole pressure calculation method proposed in our previous work [31], the bottomhole flow pressure of the horizontal well is calculated, and the data set is obtained.Figure 4
Schematic diagram of horizontal well fracturing PEBI grid.


### 5.2. Training Model

The various parameter settings in the particle swarm algorithm are described in detail in the sensitivity study section. The RBFNN consists of 5 input nodes and 1 output node. The PSO optimization steps presented in the previous section are used to optimize the training of hidden layer neurons 4, 8, 12, 16, 20, 24, and 28, respectively. The test results are shown in Figure5. As can be seen from the figure, when the number of hidden layer neurons is 16, the error is the lowest. Therefore, an RBN network structure of 5-16-1 is obtained.Figure 5
RMSE of different hidden layer neurons.


### 5.3. Prediction

After obtaining the RBF neural network parameters through the PSO algorithm, the network needs to be tested. To this end, the remaining data are imported into the optimized RBF neural network to obtain the predicted values of the bottomhole permeability and the fracture length of MFHW.In order to further detect the performance of the network, the PSO-RBF network structure proposed in the paper is compared with the PSO-BPNN and the SVR algorithm. The result is shown in Figure6. The permeability prediction results and errors of the three algorithms are shown in Table 1. As can be obtained from Table 1, the RMSE of the PSO-BPNN is 0.332, the RMSE of the SVR algorithm is 0.308, and the RMSE of the PSO-RBF is 0.178. The fracture length is predicted by three algorithms, and the result is shown in Figure 7. The fracture length prediction results and errors of the three algorithms are shown in Table 2. As can be obtained from Table 2, the RMSE of the PSO-BPNN is 3.776, the RMSE of the SVR algorithm is 3.319, and the RMSE of the PSO-RBF is 2.250. It can be seen from the comparison results that the RMSE of the PSO-RBF algorithm is much lower than that of the other two algorithms, showing that the algorithm has higher prediction accuracy and better performance.Figure 6
Prediction of permeability of different algorithms.Table 1
Permeability prediction results of different algorithms.
Fracture numberMeasured valuePSO-BPNN algorithmSVR algorithmPSO-RBN algorithm11.5491.4631.2941.66320.0110.0850.0940.05632.8372.6682.7622.79344.6715.1724.8544.70151.4371.6601.6871.53261.4521.8871.8541.56471.8182.0241.6871.90081.7562.0342.1141.902919.50019.06018.96519.2251017.53917.83717.36517.3851112.63213.09212.39413.0651211.96611.67412.32112.2111313.38913.15713.25613.395145.2505.6385.0545.341155.6335.1596.2545.406Figure 7
Prediction of fracture length in different algorithms.Table 2
Prediction results of fracture lengths for different algorithms.
Fracture numberMeasured valuePSO-BPNN algorithmSVR algorithmPSO-RBN algorithm123423823623223783693823763286293288289430630430230452762802722796233236235230725725426125382682672712699254252251252101671691691661114915114714812241239237243132392412422411424825025124615305307301303


## 5.1. Computation

The bottomhole flow pressure is analyzed and calculated by using the measured data of multistage fracturing horizontal wells in an oilfield in western China. The three-dimensional PEBI grid of fracturing of underground horizontal wells is shown in Figure4. The reservoir stratum is homogeneous, the fluid flow satisfies Darcy’s law, the horizontal and vertical permeability of the fracture are different, and the horizontal well is fractured in a closed rectangular reservoir. Using the GPU-based MFHW bottomhole pressure calculation method proposed in our previous work [31], the bottomhole flow pressure of the horizontal well is calculated, and the data set is obtained.Figure 4
Schematic diagram of horizontal well fracturing PEBI grid.


## 5.2. Training Model

The various parameter settings in the particle swarm algorithm are described in detail in the sensitivity study section. The RBFNN consists of 5 input nodes and 1 output node. The PSO optimization steps presented in the previous section are used to optimize the training of hidden layer neurons 4, 8, 12, 16, 20, 24, and 28, respectively. The test results are shown in Figure5. As can be seen from the figure, when the number of hidden layer neurons is 16, the error is the lowest. Therefore, an RBN network structure of 5-16-1 is obtained.Figure 5
RMSE of different hidden layer neurons.


## 5.3. Prediction

After obtaining the RBF neural network parameters through the PSO algorithm, the network needs to be tested. To this end, the remaining data are imported into the optimized RBF neural network to obtain the predicted values of the bottomhole permeability and the fracture length of MFHW.In order to further detect the performance of the network, the PSO-RBF network structure proposed in the paper is compared with the PSO-BPNN and the SVR algorithm. The result is shown in Figure6. The permeability prediction results and errors of the three algorithms are shown in Table 1. As can be obtained from Table 1, the RMSE of the PSO-BPNN is 0.332, the RMSE of the SVR algorithm is 0.308, and the RMSE of the PSO-RBF is 0.178. The fracture length is predicted by three algorithms, and the result is shown in Figure 7. The fracture length prediction results and errors of the three algorithms are shown in Table 2. As can be obtained from Table 2, the RMSE of the PSO-BPNN is 3.776, the RMSE of the SVR algorithm is 3.319, and the RMSE of the PSO-RBF is 2.250. It can be seen from the comparison results that the RMSE of the PSO-RBF algorithm is much lower than that of the other two algorithms, showing that the algorithm has higher prediction accuracy and better performance.Figure 6
Prediction of permeability of different algorithms.Table 1
Permeability prediction results of different algorithms.
Fracture numberMeasured valuePSO-BPNN algorithmSVR algorithmPSO-RBN algorithm11.5491.4631.2941.66320.0110.0850.0940.05632.8372.6682.7622.79344.6715.1724.8544.70151.4371.6601.6871.53261.4521.8871.8541.56471.8182.0241.6871.90081.7562.0342.1141.902919.50019.06018.96519.2251017.53917.83717.36517.3851112.63213.09212.39413.0651211.96611.67412.32112.2111313.38913.15713.25613.395145.2505.6385.0545.341155.6335.1596.2545.406Figure 7
Prediction of fracture length in different algorithms.Table 2
Prediction results of fracture lengths for different algorithms.
Fracture numberMeasured valuePSO-BPNN algorithmSVR algorithmPSO-RBN algorithm123423823623223783693823763286293288289430630430230452762802722796233236235230725725426125382682672712699254252251252101671691691661114915114714812241239237243132392412422411424825025124615305307301303


## 6. Sensitivity Study

In the particle swarm optimization algorithm, the population number, particle speed, inertia weight, learning factor, and random number all have certain effects on the performance of the algorithm. According to the studies at home and abroad, the improvement of the PSO algorithm mainly focuses on the inertia weights and learning factors [38, 39]. The inertia weight plays a major role in the convergence of the particle swarm algorithm. The larger the value of wi, the stronger the global search ability and, conversely, the stronger the local search capability. c1 and c2 are the maximum step sizes of the particles to adjust to the individual optimal or group optimal direction, respectively. When the learning factor is small, the particles may linger away from the target area, while the particles can quickly move toward the target area or even exceed the target area when the learning factor is large.The values of the inertia weight and learning factor in the particle swarm algorithm can change during the optimization process. During the search process, the inertia weight is set to change with the number of iterations. The formula of inertia weight is expressed by the following equation:(9)ωi=ωmax−i×ωmax−ωminimax,where ωmax is the maximum value of inertia weight, ωmin is the minimum value of inertia weight, i is the current number of iterations, and imax is the maximum number of iterations. The calculation formula of the learning factor is shown in the following equation:(10)c1=c1,start+i×c1,end−c1,startimax,c2=c2,start+i×c2,end−c2,startimax,where c1,startandc1,end are the initial and final values of the learning factor c1; c2,startandc2,end are the initial and final values of the learning factor c2; i is the current number of iterations; and imax is the maximum number of iterations.In order to illustrate the effect of parameter settings of inertial weight and learning factor on the prediction results of the PSO-RBF model, a classic particle swarm algorithm with fixed parameter values and an improved particle swarm algorithm with parameter value changing were used to predict the permeability and the half-length of the fracture, respectively. The basic parameters of the two algorithms are shown in Table3. The permeability prediction data and fracture length prediction data of the two algorithms are shown in Tables 4 and 5, respectively. The prediction results of the two algorithms are shown in Figures 8–11, respectively. From Table 4, it can be concluded that when making permeability predictions, the RMSE of the classic algorithm is 0.294 and the RMSE of the improved algorithm is 0.178. From Table 5, it can be concluded that when predicting the fracture length, the RMSE of the classic algorithm is 3.823 and the RMSE of the improved algorithm is 2.25. Therefore, the improved particle swarm algorithm is applied to the neural network training process in this paper.Table 3
Particle swarm algorithm parameters.
AlgorithmParameter nameValueClassic PSOInertia weightωmax=1.2ωmin=0.4Learning factorc1=c2=2Position and speed range−1,1Maximum number of iterations and population sizeimax=800N=50Improved PSOInertia weightωi=ωmax−i×ωmax−ωmin/imaxωmax=1.2ωmin=0.4Learning factorc1=c1,start+i×c1,end−c1,start/imaxc2=c2,start+i×c2,end−c2,start/imaxc1,start=c2,start=0.5c1,end=c2,end=2.5Position and speed range−1,1Maximum number of iterations and population sizeimax=800N=50Table 4
Permeability prediction data of two algorithms.
Fracture numberMeasured valueClassic algorithmImproved algorithm11.5491.7121.66320.0110.0610.05632.8372.3152.79344.6714.4234.70151.4371.6001.53261.4521.5821.56471.8181.6421.90081.7561.9231.902919.50019.74219.2251017.53917.16217.3851112.63212.22313.0651211.96611.57612.2111313.38913.62313.395145.2505.3615.341155.6335.1425.406Table 5
Fracture length prediction data of two algorithms.
Fracture numberMeasured valueClassic algorithmImproved algorithm123423023223783723763286291289430630230452762732796233236230725726125382682652699254250252101671701661114914514812241238243132392362411424825224615305302303Figure 8
Permeability prediction of classic algorithm.Figure 9
Permeability prediction of improved algorithm.Figure 10
Fracture length prediction of classic algorithm.Figure 11
Fracture length prediction of improved algorithm.


## 7. Conclusions

In this study, a neural network model of PSO-RBN was proposed to predict the permeability and fracture length in a MFHW in western China. A comparison of the performance of the PSO-RBN algorithm and core permeability measurements shows that the presented model can adequately estimate reservoir permeability and fracture length despite the highly nonlinear relationship between reservoir parameters. In addition, the PSO-RBN algorithm and other algorithms are statistically compared in the prediction of permeability and fracture length. The results depict that the presented model has some advantages over the mentioned algorithms. Therefore, the proposed PSO-RBN model can provide more accurate and efficient predictions to reservoir physical parameters.



---

*Source: 6810903-2020-07-03.xml*