Next Article in Journal
Impact-Damage Equivalency for Twisted Composite Blades with Symmetrical Configurations
Next Article in Special Issue
Statistical Criteria of Nanofluid Flows over a Stretching Sheet with the Effects of Magnetic Field and Viscous Dissipation
Previous Article in Journal
An Optimized Abstractive Text Summarization Model Using Peephole Convolutional LSTM
Previous Article in Special Issue
On Energies of Charged Particles with Magnetic Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Pigeon-Inspired Optimisation Algorithm and Its Application in Parameter Inversion

1
School of Computer Science, Wuhan Institute of Ship Building Technology, Wuhan 430050, China
2
School of Computer Science, China University of Geosciences, Wuhan 430074, China
3
Hubei Key Laboratary of Intelligent Geo-Information Processing, Wuhan 430074, China
4
Faculty of Computer Science and Engineering, Wuhan Institute of Technology, Wuhan 430205, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(10), 1291; https://doi.org/10.3390/sym11101291
Submission received: 20 August 2019 / Revised: 6 September 2019 / Accepted: 6 September 2019 / Published: 15 October 2019

Abstract

:
Pre-stack amplitude variation with offset (AVO) elastic parameter inversion is a nonlinear, multi-solution optimisation problem. The techniques that combine intelligent optimisation algorithms and AVO inversion provide an effective identification method for oil and gas exploration. However, these techniques also have shortcomings in solving nonlinear geophysical inversion problems. The evolutionary optimisation algorithms have recognised disadvantages, such as the tendency of convergence to a local optimum resulting in poor local optimisation performance when dealing with multimodal search problems, decreasing diversity and leading to the prematurity of the population as the number of evolutionary iterations increases. The pre-stack AVO elastic parameter inversion is nonlinear with slow convergence, while the pigeon-inspired optimisation (PIO) algorithm has the advantage of fast convergence and better optimisation characteristics. In this study, based on the characteristics of the pre-stack AVO elastic parameter inversion problem, an improved PIO algorithm (IPIO) is proposed by introducing the particle swarm optimisation (PSO) algorithm, an inverse factor, and a Gaussian factor into the PIO algorithm. The experimental comparisons indicate that the proposed IPIO algorithm can achieve better inversion results.

1. Introduction

First proposed in the early 1990s, the swarm intelligence algorithm is a stochastic optimisation algorithm that is largely inspired by the phenomena of biological swarm intelligence in nature and mimics the behaviour of social animals [1]. Swarm intelligence means that a swarm constituted by a large number of simple individuals without intelligence shows intelligent behaviour through the simple cooperation of the individuals with each other. Under the guidance of appropriate evolutionary mechanisms, swarm intelligence takes an explicit form through the interaction among individuals. The swarm intelligence algorithm simulates the search and optimisation process of the individuals’ evolution or foraging process. In these algorithms, the individuals in nature are represented by the points in the search space. The objective function of the problem at issue is measured as the individuals’ ability to adapt to the environment, while the individuals’ survival of the fittest process or foraging process is analogous to the iterative process. In the process, a poor feasible solution is replaced by a good feasible solution, which gives rise to an iterative search algorithm characterised with the “generation + test” feature and represents a self-adaptive artificial intelligence technique for solving the extremum problems [2]. Based on swarm intelligence research, a variety of swarm intelligence optimisation algorithms have been proposed.
In 1991, Dorigo et al. [3] proposed the ant colony optimisation (ACO) algorithm by simulating the foraging behaviour of ant colonies in nature. The ACO algorithm enables the reorganisation of its own database according to changes in the environment to achieve the evolution of the algorithm’s solving ability. However, the ACO algorithm requires a long search time, convergences slowly, and is easily trapped in local optima.
In 1995, Kennedy and Eberhart [4] proposed the particle swarm optimisation (PSO) algorithm through the study of swarm phenomena in nature. The PSO algorithm manifests strong robustness and fast convergence when processing complex multimodal optimisation problems, but it still has the inherent shortcomings of traditional bionic intelligent algorithms and is easily trapped in local optima.
Delice et al. [5] proposed a new PSO algorithm with negative knowledge to address the mixed-model two-sided assembly line balancing problem, showing an enhanced solving ability. Chatterjee [6] proposed a PSO-based approach to train the neural network (NN-PSO) to address the issues of local optimum trapping and premature convergence, common to traditional learning methods. The NN-PSO can solve the failure problem of multi-story reinforced concrete structures by detecting the failure probability of the structures. This detection is of great significance in the prevention of enormous losses due to structural failure.
In long-range homing, pigeons need the help of both the geomagnetic field and the landform, using the former to determine the general direction and the latter to amend the actual direction, thereby achieving the purpose of accurately locking the position. Inspired by the homing behaviour of pigeons in nature, in 2014, Duan et al. [7] proposed the pigeon-inspired optimisation (PIO) algorithm, a pigeon homing behaviour-based swarm intelligence optimisation algorithm that has made remarkable achievements in recent years in various fields, such as unmanned aerial vehicle (UAV) formation, parameter change control, and image processing.
Duan et al. [8] proposed a cooperative control method for UAV tight formation based on the optimisation of the predatory escape behaviour of pigeons. Following the basic idea of PIO, the structure is adjusted, and the predatory escape mechanism is introduced to improve the overall performance and to address the issue of easily trapping the PIO in local optima. The three-dimensional route planning for unmanned combat aerial vehicles (UCAVs) is a complex optimisation problem that focuses on the routing optimisation of aircraft with different constraints in complex, dynamic environments [9]. To address this problem, Zhang and Duan proposed a new predator-prey PIO (PPPIO) algorithm, which improves the global optimisation performance and convergence rate by adopting the concept of predator and prey. Dou et al. [10] proposed a new method to obtain the optimal parameters for the carrier-based aircraft controller, wherein the PIO is used in the optimisation of the prediction control parameters of the model to quickly obtain the global optimal solution while enabling easier and more effective choice of parameters through the PIO algorithm. Duan et al. [11] used the PIO to optimise the parameters of the echo state neural network to restore blurred and noisy images. In addition, the PIO improved using the orthogonal strategy which exhibits higher adaptability in performance than the existing bionic optimisation algorithms. Saranya et al. [12] proposed a PIO-based method to obtain the optimal threshold for automatic detection and extraction of the features of follicular cysts from ovary images. They demonstrated through contrast experiments that the proposed method can find better solutions with faster convergence rates. Lei et al. [13] successfully applied the PIO method to protein complex detection, which is key to understanding the principles of cellular organisation. By optimising the parameters of ε and MinPts in density-based spatial clustering of applications with noise (DBSCAN) using PIO, Lei et al. developed a new method, the pigeon-inspired optimisation DBSCAN (P_DBSCAN) to detect protein complexes in the dynamic protein-protein interaction network. In recent years, the PIO algorithms have been extensively applied.
The main contributions of this paper are as follows:
  • An improved PIO algorithm (IPIO) is proposed by introducing the particle swarm optimisation (PSO) algorithm, an inverse factor, and a Gaussian factor into the PIO algorithm. In order to verify the effectiveness of the improvement, the benchmark functions are used to do the experiments and compare the experimental results with other algorithms.
  • The improved algorithm for the pre-stack AVO elastic parameter inversion is used. The experimental comparison results with other optimisation algorithms indicate that the proposed IPIO algorithm can achieve better inversion results.

2. PIO Algorithm

2.1. PIO Algorithm and Its Improvement

The main factors affecting the pigeons’ homing behaviour are the sun, the Earth’s magnetic field, and landmarks. The pigeons use different navigation tools at different stages of the journey. Studies on the pigeons’ navigation behaviour have revealed that in the initial stage of homing, the pigeons use the compass-like navigational tool most of the time. In the middle part of the journey, they switch to the landmark navigation tool and, at the same time, re-evaluate their route and make necessary corrections.
According to the pigeons’ homing behaviour and by imitating the mechanism that pigeons use different navigation tools at different stages, two different operator models are proposed for the PIO algorithm [7] as follows:
(1) The map and compass operator. The pigeons perceive the Earth’s magnetic field via magnetic objects and then form a map in their mind. They use the height of the sun as a compass to adjust their flight direction, and as they approach the destination, their dependence on the sun and magnetic objects is reduced.
Similar to the PSO algorithm, in the PIO algorithm, pigeons have their own position and velocity information, with the following updated Formulas (1) and (2) for each pigeon.
V i N c = V i N c 1 e R × N c + r a n d ( X g b e s t X i N c 1 )
X i N c = X i N c 1 + V i N c
where R is the map and compass factor, with a range of 0 to 1; rand is a random number in the range of 0 to 1; Nc is the current number of iterations; and Xgbest is the global optimal position obtained by comparing the positions of all the pigeons after Nc-1 iteration cycles. Once the required number of iterations is reached, the operation of the map and compass operator is stopped and the landmark operator is activated.
(2) The landmark operator. The landmark operator is used to simulate the influence of landmarks on the pigeons’ navigation. When pigeons fly near the destination, they rely more on nearby landmarks. If they are familiar with the landmarks, they fly directly to the destination. Otherwise, they follow the pigeons that are familiar with the landmarks. In each iteration, the number of pigeons is halved by discarding the pigeons that are far from the destination since they are unfamiliar with the landmarks and no longer have the ability to distinguish the route. The central position of the remaining pigeons is Xcenter, which is used as a landmark, i.e., as a reference flight direction. Accordingly, the position of the ith pigeon (Xi) is updated using the Formulas (3)–(5).
X c e n t e r N c 1 = i = 1 N N c 1 X i N c 1 F ( X i N c 1 ) N N c 1 i = 1 N N c 1 F ( X i N c 1 )
N N c = N N c 1 2
X i = X i N c 1 + r a n d ( X c e n t e r N c 1 X i N c 1 )
where, F ( X i N c 1 ) = { 1 f i t n e s s ( X i N c 1 ) + ε , F o r   m i n i m i s a t i o n   p r o b l e m s f i t n e s s ( X i N c 1 ) , F o r   m a x i m i s a t i o n   p r o b l e m s . Similarly, once the maximum number of iterations is reached, the operation of the landmark operator is stopped.
The PIO algorithm has a fast convergence rate but can be very easily trapped in local optima. Moreover, the population at every step shows an aggregation state, making it impossible to perform a large number of iterations. In general, after the iterations of the first operator, all the individuals are essentially overlapped, leading to the situation where the second operator has no chance to have a turn. To address the above shortcomings of the PIO algorithm, some modifications have been made.
(1) Hybrid particle swarm algorithm
In the map and compass operator of the PIO algorithm, all the individuals move towards the optimal or sub-optimal individual. When the initial position is not near the true global optimal solution, it can cause other individuals to search in the wrong space along the route, and ultimately, be unable to find the optimal solution. Moreover, the population exhibits an aggregation state at each step, decreasing the diversity.
The PIO algorithm is similar to the PSO algorithm in that both simulate the flight characteristics of living creatures, both use the position of the current optimal individual as a reference, and both have the advantages of fast convergence and simple implementation. The differences between the two are that in the PSO algorithm, individuals not only use the position of the current optimal individual as a reference, but also take their own optimal positions as references, using the updated Formulas (6) and (7).
V i N c = ω V i N c 1 + c 1 r a n d ( X p b e s t X i X c 1 ) + c 2 r a n d ( X g b e s t X i X c 1 )
X i N c = X i X c 1 + V i N c
Compared with the map and compass operator of the PIO algorithm, the speed update formula of the PSO algorithm introduces the individuals’ own optimal positions, which can avoid the aggregation state of all the particles, and thus increase the diversity of the swarm. Therefore, for the improved PIO algorithm (IPIO) proposed in this study, the original map and compass operator are replaced by the PSO algorithm.
(2) Reverse learning mechanism
After the operation of the map and compass operator, if all the individuals are basically overlapped, the landmark operator is likely to be invalid, causing the algorithm to be trapped in local optima, and thus seriously affecting its performance. To avoid such a situation, reverse learning is introduced to the algorithm and combined with particle learning of the particle swarm to effectively balance the global search and local mining capabilities.
The concept of reverse learning was proposed by Tizhoosh in 2005, who argued that the efficiency of the algorithm can be enhanced by simultaneously considering the current solution and its reverse solution in the search process, and then selecting the better one of the two [14]. In the proposed algorithm, the elitist selection strategy based on the current solution and the reverse solution is used to accelerate the convergence of the algorithm. To avoid local optima and improve algorithm performance, a reverse learning strategy is introduced to perform general reverse learning on the individuals’ position attributes and the next-generation population is selected through a comparison. In other words, based on the original population, a reverse population is generated, in which the individuals correspond pair-wisely to those in the original population. When making selections, the fitness value of each individual of the two populations is calculated, and if the reverse individual of an individual is superior to the individual itself in terms of the fitness value, the reverse individual is retained, and, otherwise, the original individual is retained. The next-generation population thus selected remains unchanged in size. The reverse solution is obtained using Formula (8).
X i = a i + b i X i
where a i and b i are the upper and lower thresholds, respectively, of the solution. A comparison of the current solution with its reverse solution causes the solution to be more complete and avoids entrapment in local optima while ensuring that the second operator (i.e., the landmark operator) plays a better role.
(3) Gaussian factor.
The Gaussian factor has been introduced into the PIO algorithm [15], and it was noted that its use of random numbers with a uniform distribution results in, to some extent, a full search capability. A Gaussian distribution, also known as a normal distribution, has been widely used in natural and social sciences. It has two parameters, i.e., μ and σ2. The parameter μ is the mean or expected value (which is also the mid-value and the mode) of the distribution and the parameter σ is the standard deviation (and its variance is σ2).
In many cases, the optimisation algorithm should have the ability to ensure that the destination is correct when carrying out a focused search. However, this requirement cannot be satisfied using a method for obtaining random numbers with a uniform distribution. The search equation in the landmark operator of the PIO satisfies the latent premise of a Gaussian distribution and can be improved to achieve global optima. Therefore, the Gaussian random numbers are introduced into the landmark operator to enhance their ability to search for the optima.
The position update formula of the improved landmark operator is as Formulas (9) and (10).
{ X i ( T ) = X i ( T 1 ) + 2 ( R 1 0.5 ) · ( x c ( T ) X i ( T 1 ) ) · m n i f ( R 2 > p ) X i ( T ) = X i ( T 1 ) + 2 ( R 1 0.5 ) · ( x c ( T ) X i ( T 1 ) ) · 2 n i f ( R 2 p )
{ m = | R n | n = 0.5 0.25 T T 2
where p is a flexible parameter that is used to balance a Gaussian distribution and a uniform distribution, R1 and R2 are two random numbers between 0 and 1, Rn is a Gaussian random number between 0 and 1, and T2 is the maximum number of iterations.
From what has been discussed above, the steps of the proposed IPIO algorithm are as follows and the flow chart of the improved algorithm is shown in Figure 1.
  • Initialise the parameters, initialise the population, calculate the fitness value of each individual in the population, and select the optimal position X g b e s t of the population.
  • Update the position and speed of each individual according to the PSO algorithm, and calculate the position of the individual’s corresponding reverse individual.
  • Calculate the fitness value, compare the fitness values of the individual and its reverse individual, retain the better one of the two, remove the poor one, and update the global optimum and the historical optimum of each individual.
  • Determine whether the maximum number of iterations of the particle swarm operator is reached. If yes, proceed to the next step, otherwise return to Step 2.
  • Calculate the centre position of the population using the landmark operator of the PIO algorithm.
  • Update the position of each individual based on the improved landmark operator.
  • Calculate the fitness value, and update the global optimum.
  • Determine whether the maximum number of iterations of the landmark operator is reached. If yes, terminate the operation, otherwise return to Step 6 to operate until the iteration is stopped.

2.2. Experimental Results and Analysis

To verify the effectiveness of the proposed IPIO algorithm and compare it with the basic genetic algorithm (GA), the differential evolution (DE) algorithm, the PSO algorithm, and the PIO algorithm, tests were performed using a variety of single-modal and multi-modal test functions, as shown in Formulas (11)–(19).
f 1 ( x ) = i = 1 D t i + 200 D t i = { 160 + y i 2 i f ( y i < 0 ) 160 15 ( y i 15 ) i f ( 0 y i 15 ) 200 5 ( 15 y i ) i f ( 15 y i 20 ) 200 + ( y i 20 ) 2 i f ( y i > 20 ) y = x + 20
f 2 ( x ) = i = 1 D t i + D t i = { y i 2 i f ( y i < 0 o r y i > 1 ) exp [ 2 log ( 2 ) ( y i 0.1 0.8 ) 2 ] sin 6 ( 5 π y i ) i f ( 0 y i 1 ) i = 1 , 2 , , D y = x + 0.1
f 3 ( x ) = i = 1 D [ x i 2 10 cos ( 2 π x i ) + 10 ]
f 4 ( x ) = x 1 2 + 10 6 i = 2 D x i 2
f 5 ( x ) = i = 1 n 1 [ ( x i 2 ) x i + 1 2 + 1 + ( x i + 1 2 ) x i 2 + 1 ] x i [ 1 , 4 ]
f 6 ( x ) = i = 1 n [ ( x i 1 ) 2 + ( x 1 x i 2 ) 2 ] x [ 0 , 10 ]
f 7 ( x ) = i = 1 D t i + 200 D t i = { 200 + x i 2 i f ( x i < 0 ) 80 ( 2.5 x i ) i f ( 0 x i < 2.5 ) 64 ( x i 2.5 ) i f ( 2.5 x i < 5 ) 64 ( 7.5 x i ) i f ( 5 x i < 7.5 ) 28 ( x i 7.5 ) i f ( 7.5 x i < 12.5 ) 28 ( 17.5 x i ) i f ( 12.5 x i < 17.5 ) 32 ( x i 17.5 ) i f ( 17.5 x i < 22.5 ) 32 ( 27.5 x i ) i f ( 22.5 x i < 27.5 ) 80 ( x i 27.5 ) i f ( 27.5 x i 30 ) 200 + ( x i 30 ) 2 i f ( x i > 30 )
f 8 ( x ) = i = 1 D t i + D t i = { y i 2 i f ( y i < 0 o r y i > 1 ) sin 6 ( 5 π y i ) i f ( 0 y i 1 ) i = 1 , 2 , , D y = x + 0.1
f 9 ( x ) = i = 1 , 3 , 5 , D 1 [ ( y i 2 + y i + 1 11 ) 2 + ( y i + y i + 1 2 7 ) 2 ] y i = { x i + 3 i f   i   i s   o d d   n u m b e r x i + 2 i f   i   i s   e v e n   n u m b e r i = 1 , 2 , , D D   m u s t   b e   a n   e v e n   n u m b e r
The independent variable range, dimension, optimal value, and solution category of each test function are shown in Table 1.
The experimental results of the test functions are shown in Table 2.
The above experimental comparison indicates that for test Functions 1, 3, 5, 6 and 7, the IPIO algorithm significantly outperforms other methods (GA, DE, PSO, and PIO), and specifically, the results for test Functions 5 and 6 using the IPIO algorithm are very near to the optimal values. For test Functions 2, 8, and 9, the IPIO algorithm showed a slightly higher performance than the GA, while the DE, the PSO, and the PIO algorithms performed poorly. For test Function 4, all algorithms performed poorly during optimisation. In summary, the measurements of nine multimodal functions indicate that the proposed IPIO algorithm has a high optimisation ability and can be used to solve the pre-stack amplitude variation with offset (AVO) elastic parameter inversion problem.

3. Pre-stack AVO Elastic Parameter Inversion Problem

In seismic exploration, seismic information is used for oil exploration. This information may reflect the trend of change in reservoir parameters. Therefore, it is also used to predict reservoir parameters. Pre-stack seismic data contain more fluid information than post-stack seismic data, and the pre-stack inversion method presents significant advantages, such as a stable result, high resolution, and strong controllability. The AVO technique uses the elastic wave theory as the basis to study and analyse the characteristics of the variation of the seismic reflection amplitude with offset using the pre-stack common depth point to understand the relation between the reflection coefficient and the incident angle. This is then used to analyse the lithological characteristics and physical parameters above and below the reflection interface, and further predict and judge the fluid properties and lithology of the reservoir [16,17,18,19]. Pre-stack seismic data contains numerous useful information that can be used to predict underground oil and gas conditions, of which three elastic parameters, i.e., P-wave velocity V p , S-wave velocity V s , and density ρ , are key parameters that indirectly reflect the saturation state of underground oil and gas [20,21,22,23,24,25]. By using the AVO information to solve the approximate formula of the Zoeppritz equation, the pre-stack inversion obtains the elastic parameters that reflect the characteristics of underground rock directly, i.e., P-wave velocity, S-wave velocity and density. The inversion of pre-stack AVO elastic parameters requires the construction and then the optimisation of a suitable objective function, which is generally nonlinear. In theory, multiplicity is present in inversion. In the inversion process, multiple sets of solutions that can fit the observation data to the same extent can be obtained. However, it is difficult to confirm which set of solutions is closer to the real model. Therefore, the inversion problem exhibits multiplicity and non-uniqueness. Relative to linear inversion methods, intelligent optimisation algorithms do not depend on the initial model and are easy to implement. Therefore, they are suitable for the inversion of multi-parameter and multi-extremum geophysical problems. In the mid-1980s, nonlinear global intelligent optimisation inversion techniques began to receive the attention of researchers in the geophysical field. New ideas and methods in other areas have been continuously introduced since that time and have many successful applications [22,23,24,25,26,27,28,29,30,31,32,33,34,35,36].
The basic steps for solving the pre-stack AVO elastic parameter inversion problem are as follows. First, the three parameters are obtained through a series of operations, next the reflection coefficient is calculated by substituting the three parameters into the approximate equation, then the seismic data are obtained by convolving the reflection coefficient and the seismic wavelet, and finally, the obtained seismic data are compared with the actual seismic data to determine the inversion accuracy, i.e., when the two sets of data are close and the three calculated parameters are similar to the three actual parameters, the inversion accuracy is considered high. Since the three calculated parameters are within the range calculated based on well-logging data, and three smooth curves can be generated by connecting the data of all the layers corresponding to the three parameters, the pre-stack AVO elastic parameter inversion can also be viewed as a global numerical optimisation problem with continuous variables.

3.1. Inversion Model

One of the main steps in solving the pre-stack AVO elastic parameter inversion problem is to establish the inversion convolution model. The main steps of establishing the convolution model are as following: First, calculate the reflection coefficient R pp . In this study, the approximation formula of Aki and Richard [37] is used to calculate R pp as Formula (20).
R pp ( θ ) = 1 2 ( 1 + tan 2 θ ) Δ V p V p ¯ ( 4 γ 2 sin 2 θ ) Δ V s V s ¯ + 1 2 ( 1 4 γ 2 sin 2 θ ) Δ ρ ρ ¯
where Δ V p , Δ V s and Δ ρ are the differences of V p , V s and ρ , respectively, between the upper and lower layers; V p ¯ , V s ¯ and ρ ¯ are the mean values of V p , V s and ρ , respectively, of the upper and lower layers; θ is the angle; and γ = V s V p , which is calculated based on the actual measurements. Therefore, the R pp calculated using this formula can be used as a component of the seismic record convolution operation.
Second, obtain the seismic wavelet. The seismic wavelet is another component of the seismic record convolution model. Seismic record data can be obtained by convolving the seismic wavelet and the reflection coefficient, which is suitable for establishing the forward model and synthesising seismic data gathers. In this study, the Ricker wavelet, a zero-phase seismic wavelet, is adopted, using Formula (21).
f ( t ) = ( 1 2 π 2 V m 2 t 2 ) e π 2 V m 2 t 2
where V m is the main frequency, and t is the time, which can be set manually.
Third, perform the convolution operation on the Ricker wavelet and the reflection coefficient, using Formula (22).
s ( θ ) = R pp ( θ ) f ( t ) + n ( t )
where R pp ( θ ) is the reflection coefficient function, f ( t ) is the seismic wavelet, and n ( t ) is the noise, which is neglected in this study. The calculated s ( θ ) is used to construct the objective function.

3.2. Inversion Results Evaluation

The simulation-optimisation method converts the pre-stack AVO elastic parameter inversion problem to an optimisation problem, which is then solved using the optimisation algorithm. From the perspective of optimisation, when the difference between the inversion seismic data generated by the optimised elastic parameter and the actual seismic record data is 0 or less than a specified threshold, the elastic parameter is considered satisfactory. As the optimisation algorithm judges the performance of each individual based on the fitness function converted from the objective function, the performance of the objective function constructed specifically for the inversion problem is the main factor affecting the effectiveness of the pre-stack AVO elastic parameter inversion. In this study, the approximation equation of Aki and Richard is first used to obtain R p p , i.e., the reflection coefficient of the p-wave, then the synthetic seismic record data is obtained by convolving R p p and the wavelet. Assume the number of sampling points is n, and for each point, m different angles are required, giving rise to n × m entries of seismic record data. Lastly, the sum of squares of the differences between the m sets of seismic record data obtained through the optimisation and the actual record data of each sampling point is calculated, which is then divided by m, and the sum of squares of difference data of n sampling points is then divided by n. Taking the square root of the final result yields the solution. Based on the above formulas, the inversion objective function is Formula (23).
f ( x ) = i = 1 n j = 1 m ( s ( θ i , j ) s ( θ i , j ) ) 2 n × m
where s ( θ i , j ) is the forward seismic record and s ( θ i , j ) is the inverse seismic record.
Since the reflection coefficient R pp is calculated through the elastic parameters V p , V s and ρ , and the same set of inverse seismic record data obtained through the calculation are generated from infinite combinations of V p , V s and ρ , there are cases in which the inverse seismic record data are calculated with all three parameters having errors identical to those calculated with one parameter having errors and two parameters not having errors. To better evaluate the performance of the optimisation algorithm in the pre-stack AVO elastic parameter inversion, the Pearson product-moment correlation coefficient (PPMCC or PCCs) is adopted to measure the correlation between the inversed three parameters and the actual three parameters. The higher the absolute value of the correlation coefficient is, the stronger the correlation; the closer the correlation coefficient is to 1 or −1, the stronger the correlation; and the closer the correlation coefficient is to 0, the weaker the correlation. The correlation coefficient function is Formula (24).
r = i = 1 n ( X i X ¯ ) ( Y i Y ¯ ) i = 1 n ( X i X ¯ ) 2 i = 1 n ( Y i Y ¯ ) 2
where X i is the standard value of a certain parameter of the three parameters, and Y i is the corresponding inversion value, and X ¯ , Y ¯ are the means of a set of values.
From the complexity of the solution process of seismic data, a small objective function does not necessarily mean a high correlation coefficient of the three parameters, and vice versa. However, when the objective function value reaches its theoretical optimal value of 0, the correlation coefficients of V p , V s and ρ also reach their theoretical optimal value of 1. Therefore, in this study, the objective function and the correlation coefficient are combined to evaluate the performance of the inversion. The ultimate goal is to have a small objective function through the inversion and simultaneously, the correlation coefficients for the three parameters are large.

4. Experimental Simulation and Analysis

4.1. Parameter Setting

The algorithm’s parameter setting and experimental environment parameters are shown in Table 3 and Table 4.

4.2. Simulation Experiment

The well-logging curve data in Dataset 1 are collected from 241 sampling points, including the p-wave velocity V p , the S-wave velocity V s , and the density ρ . At each sampling point, the data from eight different angles (0°, 6°, 11°, 17°, 23°, 29°, 34°, and 40°) were acquired and included in each dataset. The formula of Aki and Richard is used to forward the theoretical model of the logging curve, which is then used to calculate the reflection coefficient that is subsequently further convolved with the wavelet. When generating the seismic record data, it is necessary to use the relation between the upper and lower groups of sampling points. Therefore, the seismic record data contain 240 × 8 entries. The well-logging data thus generated are shown in Figure 2 and Figure 3.
The pre-stack AVO elastic parameter inversion was performed, respectively, using the proposed IPIO algorithm as well as the GA, PSO, DE and PIO algorithms. The experiments were conducted according to the algorithm parameter settings. The experimental results are shown in Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8.
The mean correlation coefficients of the three parameters ( V p , V s and ρ ) obtained from the experiments on the five algorithms (i.e., GA, PSO, DE, PIO and IPIO) are compared. The results are shown in Figure 9.
The comparison results shown in Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 and Table 5 indicate that the proposed IPIO algorithm significantly outperforms the GA, PSO, DE and PIO algorithms, with a drastically increased correlation coefficient. Form Table 5, it can be seen for the GA algorithm, the correlation coefficients for the three elastic parameters V p , V s and ρ are: 0.556364, 0.650805 and 0.462346. For the PSO algorithm, the correlation coefficients for the three elastic parameters V p , V s and ρ are: 0.765481, 0.832.94 and 0.649089. For the DE algorithm, the correlation coefficients for the three elastic parameters V p , V s and ρ are: 0.701792, 0.787196 and 0.606702. For the PIO algorithm, the correlation coefficients for the three elastic parameters V p , V s and ρ are: 0.980645, 0.894879 and 0.537096, but for the proposed algorithm, the correlation coefficients for the three elastic parameters V p , V s and ρ are: 0.937951, 0.908213 and 0.925746. The comparison results shown in Figure 9 indicate that the proposed IPIO algorithm effectively improves the correlation coefficients for the three elastic parameters V p , V s and ρ , indicating that the inversed seismic data are closer to the actual seismic data.
The objective function values from the inversion are summarised in Figure 10. The original PIO algorithm has the worst performance in solving the inversion problem. The performance of the DE algorithm in solving the inversion problem is not satisfactory. GA outperforms the PIO and DE algorithms, with an objective function value of approximately 0.0035. The objective function value of the PSO algorithm is approximately 0.0018, which is better than the GA. The objective function value of the IPIO algorithm reaches 0.00094, which is an order of magnitude better than the other algorithms. The proposed IPIO algorithm greatly reduces the objective function value of the inversion. The comparison results in Figure 10 indicate that the proposed IPIO algorithm shows a rapid convergence rate and has a higher optimisation effectiveness than the GA, PSO and DE algorithms.

5. Conclusions

This study proposed an improved PIO algorithm that is more suitable for solving the inversion problem of pre-stack AVO elastic parameters. The proposed algorithm combines the PSO algorithm, the reverse particle, and the Gaussian factor, along with an improved initialisation strategy, to optimise the inversion, thus improving inversion accuracy. Through a large number of experiments, the authors confirm that the proposed IPIO algorithm is able to improve the accuracy of the inversion, making the inversed seismic data more consistent with the actual seismic data.
Future research directions can be divided into two levels: (1) At the algorithm level, more improvement strategies can be proposed for the PIO algorithm, so as to enhance the optimization ability of the algorithm; (2) At the problem level, the elastic parameters of the pre-stack seismic data inversion problem are complex. The actual cases and the geological structure are complex, from drilling to data acquisition. Each step has considerable errors. This article uses the data as an ideal data without error, and the actual situation, access to data, logging data and seismic data have some noise error. The inversion further increases the difficulty, and how to develop with noise error of pre-stack seismic data inversion of elastic parameters. Therefore, further research is needed.

Author Contributions

H.L. and Q.W. designed the research; X.Y. and H.L. conducted the experiments; H.L. and Q.W. analyzed the data; H.L. wrote the manuscript; Q.W. promoted the manuscript.

Funding

This research was funded by Natural Science Foundation of China (No. 61673354), the Fundamental Research Funds for the Central Universities, China University of Geosciences (Wuhan) (CUGGC03), the State Key Lab of Digital Manufacturing Equipment & Technology, Huazhong University of Science & Technology (DMETKF2018020) and Open Research Project of the Hubei Key Laboratory of Intelligent Geo-Information Processing (China University of Geosciences (Wuhan)).

Conflicts of Interest

The authors declare there is no conflicts of interest.

References

  1. Kennedy, J. Swarm Intelligence. Swarm Intelligence; Morgan Kaufmann Publishers Inc.: San Francisco, CA, USA, 2001. [Google Scholar]
  2. Krause, J.; Ruxton, G.D.; Krause, S. Swarm intelligence in animals and humans. Trends Ecol. Evol. 2010, 25, 28–34. [Google Scholar] [CrossRef] [PubMed]
  3. Dorigo, M.; Stützle, T. Ant Colony Optimization: Overview and Recent Advances. In Handbook of Metaheuristics; Springer: Cham, Switzerland, 2010. [Google Scholar]
  4. Eberhart, R.; Kennedy, J. A New Optimizer Using Particle Swarm Theory. In Proceedings of the Sixth International Symposium, Nagoya, Japan, 4–6 October 1995; pp. 39–43. [Google Scholar]
  5. Delice, Y.; Aydoğan, E.K.; Özcan, U.; İlkay, M.S. A Modified Particle Swarm Optimization Algorithm to Mixed-Model Two-Sided Assembly Line Balancing. J. Intell. Manuf. 2017, 28, 23–36. [Google Scholar] [CrossRef]
  6. Chatterjee, S.; Sarkar, S.; Hore, S.; Dey, N.; Ashour, A.S.; Balas, V.E. Particle Swarm Optimization Trained Neural Network for Structural Failure Prediction of Multistoried RC Buildings. Neural Comput. Appl. 2017, 28, 2005–2016. [Google Scholar] [CrossRef]
  7. Duan, H.; Qiao, P. Pigeon-Inspired Optimization: A New Swarm Intelligence Optimizer for Air Robot Path Planning. Int. J. Intell. Comput. Cybern. 2014, 7, 24–37. [Google Scholar] [CrossRef]
  8. Duan, H.B.; Qiu, H.X.; Fan, Y.M. Optimal Combat Coordination Control for Unmanned Aerial Vehicle Using Predator-Escaped Pigeon. Chin. Sci. Technol. Sci. 2015, 45, 559–572. [Google Scholar]
  9. Zhang, B.; Duan, H. Three-Dimensional Path Planning for Uninhabited Combat Aerial Vehicle Based on Predator-Prey Pigeon-Inspired Optimization in Dynamic Environment; IEEE Computer Society Press: Los Alamitos, CA, USA, 2017. [Google Scholar]
  10. Dou, R.; Duan, H. Pigeon Inspired Optimization Approach to Model Prediction Control for Unmanned Air Vehicles. Aircr. Eng. Aerosp. Technol. Int. J. 2016, 88, 108–116. [Google Scholar] [CrossRef]
  11. Duan, H.; Wang, X. Echo State Networks with Orthogonal Pigeon-Inspired Optimization for Image Restoration. IEEE Trans. Neural Netw. Learn. Syst. 2015, 27, 2413–2425. [Google Scholar] [CrossRef] [PubMed]
  12. Sankareswaran, U.M.; Rajendran, S. A Novel Pigeon Inspired Optimization in Ovarian Cyst Detection. Curr. Med. Imaging Rev. 2016, 12, 43–49. [Google Scholar]
  13. Lei, X.; Ding, Y.; Wu, F.X. Detecting Protein Complexes from DPINs by Density Based Clustering with Pigeon-Inspired Optimization Algorithm. Sci. China Inf. Sci. 2016, 59, 070103. [Google Scholar] [CrossRef]
  14. Tizhoosh, H.R. Opposition-Based Learning: A New Scheme for Machine Intelligence. In International Conference on Computational Intelligence for Modelling, Control and Automation and International Conference on Intelligent Agents, Web Technologies and Internet Commerce (CIMCA-IAWTIC’06); IEEE: New York, NY, USA, 2005; pp. 695–701. [Google Scholar]
  15. Zhang, S.; Duan, H. Gaussian Pigeon-Inspired Optimization Approach to Orbital Spacecraft Formation Reconfiguration. Chin. J. Aeronaut. 2015, 28, 200–205. [Google Scholar] [CrossRef]
  16. Rodriguez-Zalapa, O.; Huerta-Ruelas, J.A.; Rangel-Miranda, D.; Morales-Sanchez, E.; Hernandez-Zavala, A. Csimfs: An Algorithm to Tune Fuzzy Logic Controllers. J. Intell. Fuzzy Syst. 2017, 33, 679–691. [Google Scholar] [CrossRef]
  17. Alfonso, M.J.F.I.; Milán, D.P.S.; Córdoba, J.V.A.I.; Colmenero, N.P. Some Improvements on Relativistic Positioning Systems. Appl. Math. Nonlinear Sci. 2018, 3, 161–166. [Google Scholar] [CrossRef] [Green Version]
  18. Plata, S.A.; Sáez, S.B. After Notes on Chebyshev’s Iterative Method. Appl. Math. Nonlinear Sci. 2017, 2, 1–12. [Google Scholar]
  19. Neidell, N.S. Amplitude Variation with Offset. Lead. Edge 1986, 5, 47–51. [Google Scholar] [CrossRef]
  20. Gao, W.; Zhu, L.; Guo, Y.; Wang, K. Ontology Learning Algorithm for Similarity Measuring and Ontology Mapping Using Linear Programming. J. Intell. Fuzzy Syst. 2017, 33, 3153–3163. [Google Scholar] [CrossRef]
  21. Jin, J.; Mi, W. An Aimms-Based Decision-Making Model for Optimizing the Intelligent Stowage of Export Containers in a Single Bay. Discret. Contin. Dyn. Syst. Ser. S 2019, 12, 1101–1115. [Google Scholar] [CrossRef]
  22. Wu, Q.; Zhu, Z.; Yan, X. Research on the Parameter Inversion Problem of Prestack Seismic Data Based on Improved Differential Evolution Algorithm. Clust. Comput. 2017, 20, 2881–2890. [Google Scholar] [CrossRef]
  23. Yan, X.; Zhu, Z.; Wu, Q.; Gong, W.; Wang, L. Elastic Parameter Inversion Problem Based on Brain Storm Optimization Algorithm. Memet. Comput. 2019, 11, 143–153. [Google Scholar] [CrossRef]
  24. Wu, Q.; Zhu, Z.; Yan, X.; Gong, W. An Improved Particle Swarm Optimization Algorithm for AVO Elastic Parameter Inversion Problem. Concurr. Comput. Pract. Exp. 2019, 31, 1–16. [Google Scholar] [CrossRef]
  25. Yan, X.; Zhu, Z.; Hu, C.; Gong, W.; Wu, Q. Spark-Based Intelligent Parameter Inversion Method for Prestack Seismic Data. Neural Comput. Appl. 2018. [Google Scholar] [CrossRef]
  26. Mallick, S. Model-based Inversion of Amplitude-variations-with-offset Data Using a Genetic Algorithm. Geophysics 1995, 60, 939–954. [Google Scholar] [CrossRef]
  27. Priezzhev, I.; Shmaryan, L.; Bejarano, G. Nonlinear Multitrace Seismic Inversion Using Neural Network and Genetic Algorithm-Genetic Inversion. In Proceedings of the 3rd EAGE St. Petersburg International Conference and Exhibition on Geosciences-Geosciences: From New Ideas to New Discoveries 2008, Saint Petersburg, Russia, 2008. [Google Scholar]
  28. Soupios, P.; Akca, I.; Mpogiatzis, P.; Basokur, A.T.; Papazachos, C. Applications of Hybrid Genetic Algorithms in Seismic Tomography. J. Appl. Geophys. 2011, 75, 479–489. [Google Scholar] [CrossRef]
  29. Agarwal, A.; Sain, K.; Shalivahan, S. Traveltime and Constrained AVO Inversion Using FDR PSO. In SEG Technical Program Expanded Abstracts 2016; Society of Exploration Geophysicists: Tulsa, OK, USA, 2016; pp. 577–581. [Google Scholar]
  30. Sun, S.Z.; Chen, L.; Bai, Y.; Hu, L.G. PSO Non-linear Pre-stack Inversion Method and the Application in Reservoir Prediction. In SEG Technical Program Expanded Abstracts 2012; Society of Exploration Geophysicists: Tulsa, OK, USA, 2012; pp. 1–5. [Google Scholar]
  31. Sun, S.Z.; Liu, L. A Numerical Study on Non-linear AVO Inversion Using Chaotic Quantum Particle Swarm Optimization. J. Seism. Explor. 2014, 23, 379–392. [Google Scholar]
  32. Zhou, Y.; Nie, Z.; Jia, Z. An Improved Differential Evolution Algorithm for Nonlinear Inversion of Earthquake Dislocation. Geod. Geodyn. 2014, 5, 49–56. [Google Scholar]
  33. Song, X.; Li, L.; Zhang, X.; Huang, J.; Shi, X.; Jin, S.; Bai, Y. Differential Evolution Algorithm for Nonlinear Inversion of High-frequency Rayleigh Wave Dispersion Curves. J. Appl. Geophys. 2014, 109, 47–61. [Google Scholar] [CrossRef]
  34. Lu, Y.F.; Bao, T.F.; Li, J.M.; Wang, T. Inversion of Geotechnical Mechanical Parameters Based on Improved Differential Evolution Algorithm Online Support Vector Regression and ABAQUS. J. Yangtze River Sci. Res. Inst. 2017, 34, 81–87. [Google Scholar]
  35. Gao, Z.; Pan, Z.; Gao, J. Multimutation Differential Evolution Algorithm and Its Application to Seismic Inversion. IEEE Trans. Geosci. Remote Sens. 2016, 54, 3626–3636. [Google Scholar] [CrossRef]
  36. Yan, X.; Zhu, Z.; Wu, Q. Intelligent Inversion Method for Pre-stack Seismic Big Data Based on MapReduce. Comput. Geosci. 2018, 110, 81–89. [Google Scholar] [CrossRef]
  37. Aid, K.; Richards, P.G. Quantitative Seismology: Theory and Methods; W. H. Freeman and Co.: San Francisco, CA, USA, 1980. [Google Scholar]
Figure 1. The improved pigeon-inspired optimization (IPIO) algorithm flow chart.
Figure 1. The improved pigeon-inspired optimization (IPIO) algorithm flow chart.
Symmetry 11 01291 g001
Figure 2. Original logging curves.
Figure 2. Original logging curves.
Symmetry 11 01291 g002
Figure 3. Seismic records of theory model.
Figure 3. Seismic records of theory model.
Symmetry 11 01291 g003
Figure 4. Logging curves and seismic records of inversion generated by basic genetic algorithm (GA).
Figure 4. Logging curves and seismic records of inversion generated by basic genetic algorithm (GA).
Symmetry 11 01291 g004
Figure 5. Logging curves and seismic records of inversion generated by particle swarm optimisation (PSO).
Figure 5. Logging curves and seismic records of inversion generated by particle swarm optimisation (PSO).
Symmetry 11 01291 g005
Figure 6. Logging curves and seismic records of inversion generated by differential evolution (DE).
Figure 6. Logging curves and seismic records of inversion generated by differential evolution (DE).
Symmetry 11 01291 g006
Figure 7. Logging curves and seismic records of inversion generated by pigeon-inspired optimisation (PIO).
Figure 7. Logging curves and seismic records of inversion generated by pigeon-inspired optimisation (PIO).
Symmetry 11 01291 g007
Figure 8. Logging curves and seismic records of inversion generated by IPIO.
Figure 8. Logging curves and seismic records of inversion generated by IPIO.
Symmetry 11 01291 g008
Figure 9. Comparison of the mean correlation coefficients of the three parameters for the GA, PSO, DE, PIO and IPIO algorithms.
Figure 9. Comparison of the mean correlation coefficients of the three parameters for the GA, PSO, DE, PIO and IPIO algorithms.
Symmetry 11 01291 g009
Figure 10. Convergence curve of the algorithm.
Figure 10. Convergence curve of the algorithm.
Symmetry 11 01291 g010
Table 1. Test function information.
Table 1. Test function information.
FunctionRange of Independent VariableDimensionOptimal ValueCategory
f(1)[−100,100]50Single solution
f(2)[−100,100]50Single solution
f(3)[−100,100]60Single solution
f(4)[−100,100]50Single solution
f(5)[−1,4]50Single solution
f(6)[0,10]50Single solution
f(7)[−100,100]50Multiple solutions
f(8)[−100,100]50Multiple solutions
f(9)[−100,100]60Multiple solutions
Table 2. Experimental results.
Table 2. Experimental results.
FunctionAlgorithmMinimum ValueMaximum ValueMean
f(1)GA2.56 × 10−35.14 × 10−11.01 × 10−1
DE3.94 × 10−245.0859927.59336
PSO3.60 × 10−14.2241611.934762
PIO1.32 × 10−73.55 × 10−28.44 × 10−3
IPIO5.69 × 10−83.71 × 10−55.42 × 10−6
f(2)GA1.19 × 10−11.6045135.17 × 10−1
DE4.63 × 10−711.10021.684983
PSO2.2138674.6072233.56288
PIO8.84 × 10−53.3769681.703053
IPIO3.82 × 10−142.15 × 10−11.69 × 10−2
f(3)GA5.81 × 10−22.7067441.060794
DE2.44 × 10−855.000467.248065
PSO7.60982520.1941415.15247
PIO8.81 × 10−98.3284671.403083
IPIO1.07 × 10−141.15 × 10−15.73 × 10−3
f(4)GA1.90 × 10−24659.428594.2779
DE9.68 × 10−222078.023841.383
PSO11827.933459811362612.9
PIO5.26E + 02114849.423962.32
IPIO3.88 × 10−8318704.31.63E+04
f(5)GA1.67 × 10−119.68 × 10−61.81 × 10−6
DE2.29 × 10−64.11 × 10−13.09 × 10−2
PSO1.21 × 10−48.53 × 10−44.84 × 10−4
PIO1.62 × 10−72.32 × 10−33.23 × 10−4
IPIO3.87 × 10−245.45 × 10−197.82 × 10−20
f(6)GA6.56 × 10−47.93 × 10−33.40 × 10−3
DE3.35 × 10−630.110477.679355
PSO4.68 × 10−41.679668.87 × 10−2
PIO4.02 × 10−13.062161.59
IPIO9.86 × 10−326.99 × 10−156.79 × 10−16
f(7)GA2.99 × 10−32.68 × 10−21.09 × 10−2
DE7.31 × 10−25.1633481.671903
PSO3.07 × 10−14.4848611.868947
PIO0.0086.5113333.82977
IPIO7.45 × 10−96.22 × 10−54.70 × 10−6
f(8)GA2.82 × 10−42.18 × 10−13.22 × 10−2
DE8.06 × 10−68.3538499.78 × 10−1
PSO1.7720183.9917543.012487
PIO3.02 × 10−53.41 × 10−21.66 × 10−4
IPIO3.55 × 10−153.73 × 10−32.26 × 10−4
f(9)GA3.30 × 10−41.9081914.45 × 10−1
DE4.88 × 10−3295.510421.82682
PSO9.76199947.8248224.92847
PIO3.06 × 10−42.2360641.41 × 10−1
IPIO5.50 × 10−120.0585372.93 × 10−3
Table 3. Parameter settings for the IPIO algorithm.
Table 3. Parameter settings for the IPIO algorithm.
NωC1C2pMax_Iteration1Max_Iteration2
400.5220.340001000
Table 4. Experimental environment parameters.
Table 4. Experimental environment parameters.
Experimental EnvironmentParameter Description
JAVA version1.8.0_111-b14
Compiler environmentEclipse-jee-luna-SR1a-win32-x86_64
ProcessorIntel(R) Core (TM) i5-6500 CPU @ 3.10 GHZ
Installed Memory (RAM)8.00 GB
System type64-bit operating system
Table 5. Comparison of the mean correlation coefficients of the three parameters for the five algorithms.
Table 5. Comparison of the mean correlation coefficients of the three parameters for the five algorithms.
GAPSODEPIOIPIO
V p 0.5563640.7654810.7017920.9806450.937951
V s 0.6508050.8320940.7871960.8948790.908213
ρ 0.4623460.6490890.6067020.5370960.925746

Share and Cite

MDPI and ACS Style

Liu, H.; Yan, X.; Wu, Q. An Improved Pigeon-Inspired Optimisation Algorithm and Its Application in Parameter Inversion. Symmetry 2019, 11, 1291. https://doi.org/10.3390/sym11101291

AMA Style

Liu H, Yan X, Wu Q. An Improved Pigeon-Inspired Optimisation Algorithm and Its Application in Parameter Inversion. Symmetry. 2019; 11(10):1291. https://doi.org/10.3390/sym11101291

Chicago/Turabian Style

Liu, Hanmin, Xuesong Yan, and Qinghua Wu. 2019. "An Improved Pigeon-Inspired Optimisation Algorithm and Its Application in Parameter Inversion" Symmetry 11, no. 10: 1291. https://doi.org/10.3390/sym11101291

APA Style

Liu, H., Yan, X., & Wu, Q. (2019). An Improved Pigeon-Inspired Optimisation Algorithm and Its Application in Parameter Inversion. Symmetry, 11(10), 1291. https://doi.org/10.3390/sym11101291

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