Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2015, Article ID 706897, 9 pages
http://dx.doi.org/10.1155/2015/706897
Research Article

Developing a Robust Surrogate Model of Chemical Flooding Based on the Artificial Neural Network for Enhanced Oil Recovery Implications

Department of Petroleum Engineering, Ahwaz Faculty of Petroleum Engineering, Petroleum University of Technology, P.O. Box 63431, Ahwaz 6198144471, Iran

Received 10 June 2014; Revised 5 October 2014; Accepted 13 October 2014

Academic Editor: Shifei Ding

Copyright © 2015 Mohammad Ali Ahmadi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Application of chemical flooding in petroleum reservoirs turns into hot topic of the recent researches. Development strategies of the aforementioned technique are more robust and precise when we consider both economical points of view (net present value, NPV) and technical points of view (recovery factor, RF). In current study many attempts have been made to propose predictive model for estimation of efficiency of chemical flooding in oil reservoirs. To gain this end, a couple of swarm intelligence and artificial neural network (ANN) is employed. Also, lucrative and high precise chemical flooding data banks reported in previous attentions are utilized to test and validate proposed intelligent model. According to the mean square error (MSE), correlation coefficient, and average absolute relative deviation, the suggested swarm approach has acceptable reliability, integrity and robustness. Thus, the proposed intelligent model can be considered as an alternative model to predict the efficiency of chemical flooding in oil reservoir when the required experimental data are not available or accessible.

1. Introduction

The oil and gas upstream industries are recently encountered with the difficulties and challenges of dealing with hydrocarbon resources whose productions with conventional technologies are following an upward trend of technical limitations. It is because of achieving the stage of decline phase by most of oilfields around the world. Therefore, how to postpone the abandonment of reservoirs has tuned into the priority of researchers in the worldwide. Their researches normally highlight the concept of great necessities for inventions of new techniques, normally classified as tertiary oil recovery methods, having abilities of maintaining the economic production rate [13].

Chemical enhanced oil recovery approaches as one of the most effective subsets of tertiary methods are known as a key to unlock the exploitation of referred resources. Different methods for this process have been developed, such as polymer, surfactant/polymer (SP), and alkaline/surfactant/polymer (ASP) flooding. These methods are applied to increase the rate of oil production through focusing on both lowering the interfacial tension and reducing the water mobility. In more details, it has enormously been declared in previous literatures that in order to design, manage, and run a chemical enhanced oil recovery operation it is highly required to set very expensive and time-consuming but precise experimental procedures which their generated results must be gained to plan effectively the process of injecting chemical materials [49].

The laboratorial generated outputs are then used to conclude two parameters, recovery factor (RF) and net present value (NPV), which are used to evaluate the performance of the chemical flooding which is one of the most popular methods of chemical enhanced oil recovery. Having knowledge about these two parameters is essentially vital to make decisions if it is beneficial to run the referred operation. Unfortunately, there are no global methods to interpret simultaneously both aforementioned factors although there are numerous numbers of different software and numerical or analytical methods which are capable of making very precise quantitative decisions about the amount of one of the RF or NPV [1012].

Hence, there is a great need in oilfield for having access to a solution or model which can predict the amount of these two parameters at the same time. The major aim of current study is to execute new kind of artificial intelligence approaches to suggest robust and accurate predictive method to forecast efficiency of the chemical flooding through petroleum reservoirs. To gain successfully this referred goal, hybridization of artificial neural network and particle swarm optimization (PSO) was executed on the previous literature data bases. The integrity and performance of the proposed predictive approaches in estimating recovery factor (RF) and net present value (NPV) from the literature are described in details.

2. Data Gathering

The data utilized throughout this research have been gathered from previous attentions [9] in which chemical flooding had been simulated in Benoist sand reservoir, by executing UTCHEM simulator. That reservoir has been produced under primary and secondary processes over fifty years. The original dataset contained 202 data. Each data had 7 inputs: surfactant slug size, surfactant concentration in surfactant slug, polymer concentration in surfactant slug, polymer drive size, and polymer concentration in polymer drive, ratio, and salinity of polymer drive. In addition, the outputs were RF and NPV. The ranges of implemented data banks are reported in Table 1 [9].

Table 1: Statistical analysis of the implemented chemical flooding data samples [9].

3. Artificial Neural Network and Particle Swarm Optimization

Artificial neural network (ANN) includes simple nodes, named as neurons, which are bonded to each other to construct a network model. Indeed, the biological nervous systems can be simulated with the ANN system, somehow. Characterization of an ANN model is normally performed through three ways including (a) certain patterns between various layers, (b) connection between input and output via activation function, and (c) updating the interconnection weights through training process [1324].

In fact, the main purpose of an ANN model is to determine target function through internal computation during the training phase if the values of input variables are provided. The most common type of ANN is the multilayer feed forward neural network which is made up of group of interconnected neurons organized in the form of layers: input layer, hidden layer(s), and output layer where each layer comprises a group of neurons as presented in Figure 1. This network is strictly an acyclic type since signals propagate only in a forward direction from the input neurons to the output neurons and no signals are allowed to be fed-back among the neurons. The number of neurons in the input and output layers is decided by the number of input and output variables that are planned for the predictive tool. However, the optimal number of neurons in hidden layer(s) is a strong function of nonlinearity and dimensionality of the problem under study [1324].

Figure 1: Architecture of multilayer feed forward neural network. The symbol denotes the synaptic weight between the output of the th neuron in the hidden layer and the input of the th neuron in output layer. The symbol denotes the bias of the th neuron in hidden layer. The superscript stands for output layer.

The artificial neuron is the fundamental part of the neural networks. Each artificial neuron—excluding neurons at the input layer—takes and processes inputs gathered from other neurons. Given further information, each artificial neuron is a mathematical information-processing unit. The processed information is presented at the output end of the neuron. Figure 2 addresses the procedure in which an artificial neuron treats the data and information entered in the model. Each input signal () is primarily multiplied by the corresponding weight value () and the resultant products are summed up to generate a total weight in the form of . The sum of the weighted inputs and the bias () forms the input to the activation function, . An activation function processes this sum and gives out the output, . Indeed, the resulting sum is processed by a neuron activation function to obtain the ultimate output of the neuron as follows [1326]: This output will be the input signal for the neurons in the following layer. The linear (purelin) transfer, tan-sigmoid (tansig) activation, and log-sigmoid (logsig) activation functions are mostly employed in the practical cases with applications in science and engineering disciplines. The corresponding relationships for these functions are defined, respectively, by (2)–(4), as given below [1327]: The weight factors are generally considered as the adaptive parameters in the network to obtain the strength of the input signals. A bias is characterized with a weight which is not responsible for connecting an input of two neurons to an output. A particular level of a neuron output signal is represented by a set of bias that does not depend on the input signals. The weight factors and biases are tuned during the course of training phase such that the network is able to forecast the accurate target parameter for a given set of inputs. There are a number of training algorithms with different methodologies in the context of intelligence system. A variety of optimization tools such as particle swarm optimization (PSO) [15, 18, 19], genetic algorithm (GA) [21], hybrid genetic algorithm and particle swarm optimization (HGAPSO) [13, 16], unified particle swarm optimization (UPSO) [14], and imperialist competitive algorithm (ICA) [17, 20, 23] for weight training of neural networks have been used.

Figure 2: Information processing by an artificial neuron.

Kennedy [27] introduced the PSO as a strong stochastic optimization technique which simulates the social manners of birds within a group, based on population concept. It searches for an optimum solution by iteratively updating a swarm of particles.

The model originally includes a group of random particles (solutions). A random velocity is attributed to each candidate particle which flies within the problem space. The solutions consist of memory and try to attain the best position or/and fitness. This parameter is symbolized by “” that is linked only to a specific particle. The model also retains the best fitness, known as “,” which is found among the entire solutions (particles) in the swarm. The candidate particle that obtains this fitness is the global best in the population [2528]. In the current study, a particle’s fitness is calculated through determination of the network output for every point in the training part and then computing the sum of squares of the resultant errors (MSE) for performance evaluation. The basic PSO theory involves variation of each particle velocity toward its and locations at each time interval. The particles’ new velocity and position are updated according to the following equations [1328]: where and are velocities of particle at iterations and ; and are positions of particle at iterations and ; represents the inertia weight that directs the exploitation and exploration of the search space as it continuously updates velocity; and are termed as cognition and social components, respectively. They are considered as the acceleration constants which alter the velocity of a solution in the direction of and [1328]; and and refer to the two random variables uniformly distributed in the interval of .

Herein, PSO algorithm has been used in evolving weights of multilayer feed forward neural network. In this case, a particle’s position at any iteration is described as a particle whose coordinates are connection weights. The vectors of weights for each particle will be called . Throughout the training process the above equations (equations (5) and (6)) will customize the network weights until a criterion is met. In this case, a lower MSE, as a sufficiently good fitness, is achieved; nevertheless, a maximum number of iterations are used to terminate the iterative search process if no improvement is observed over a number of consecutive generations in an appropriate time. The flowchart of the PSO-based training algorithm for the ANN is shown in Figure 3.

Figure 3: PSO-based algorithm flowchart in optimization of the weights of ANN.

The PSO utilizes a random procedure in the search space of the problem such that particles in the population are directed toward optimum positions but not in or between optimal areas [27]. Thus, PSO can be used to train neural networks with nondifferentiable (even discontinuous) neurons activation functions. It can be also implemented in cases where gradient or error information is not accessible. PSO is easy to implement and there are few parameters to be adjusted. However, the uniqueness of the algorithm lies in the dynamic interactions among the particles that turn it into a social-psychological model of knowledge management [27].

4. Results and Discussion

According to the study accomplished by Cybenko [29], a network that consists of only one single hidden layer has the ability to approximate nearly any kind of nonlinear function. However, determination of the ideal number of neurons in the hidden layer is a challenging task; few neurons will not give adequate precision and too many hidden neurons may lead to overfitting. It means that the training data might be fitted adequately; however considerable oscillations between the points are noticed in the fitting curve, resulting in poor interpolation and extrapolation. The network performance is evaluated as demonstrated in Figure 4, when different number of neurons is tested. A smart model with one hidden layer (including just one neuron) was primarily built in the current study to predict the recovery factor and net present value (NPV) of chemical flooding in oil reservoirs. Prediction accuracy was further analyzed by an increase in the number of neurons to 10 to decide on the most precise technique. As clear from the results demonstrated in Figure 4, a 3-6-1 architecture (6 neurons in the hidden layer, 3 neurons in the input layer, and one neuron in output layer) offers the best model for recovery factor and net present value (NPV) prediction in terms of and , since the optimum structure achieved by the trial and error procedure has a very low mean squared error of and a satisfactory coefficient of determination of , on the basis of comparison between the predicted and real data.

Figure 4: Effect of number of hidden neuron on PSO-ANN accuracy of (a) recovery and (b) NPV predictions in terms of MSE and -squared.

The generated results of the proposed intelligent approach are depicted through Figures 5 to 10. The existing contrasts between suggested intelligent approach and related recovery factor (RF) of the chemical flooding in oil reservoir in the regression plot have been depicted in Figure 5. As shown in Figure 5 which is a graphical and scatter presentation of the PSO-ANN results versus corresponding determined recovery factor (RF) data, the PSO-ANN outputs lie over the line , the fact that indicates the identity of outputs gained from suggested PSO-ANN model and relevant recovery factor data samples. To serve better understanding about generated results of the proposed PSO-ANN model, the comparison between gained recovery factor from the addressed model and real recovery factor data versus relevant data index has been illustrated in Figure 6. As illustrated in Figure 6, the obtained results of proposed model are as close as possible to real recovery factor (RF) data samples. To put it another way, the outputs of the PSO-ANN approach have the same behaviour as actual data do. The high considerable level of efficiency and accuracy related to the PSO-ANN approach in prediction of the recovery factor dataset of chemical flooding has once again been certified in Figure 6. Moreover, the robustness of the PSO-ANN has been demonstrated in terms of the relative deviations of PSO-ANN model outputs from corresponding determined recovery factor data in Figure 7. As could be observed in Figure 7, the highest deviations of the suggested approach results are subjected to the early boundary of recovery factor data samples. 5% is the maximum degree of relative deviation shown in Figure 7.

Figure 5: Performance plot of the suggested network model for determining recovery factor of chemical flooding owing to correlation coefficient (): (a) training phase and (b) testing phase.
Figure 6: Comparison between suggested network model and recovery factor versus relevant data index: (a) training phase and (b) testing phase.
Figure 7: Relative error distribution of the proposed approach versus actual recovery factor (RF).

The draw parallel between our proposed intelligent PSO-ANN model results and related net present value (NPV) of the chemical flooding in oil reservoir in the regression plot has been shown in Figure 8. As shown in Figure 8 which is a graphical and scatter presentation of the PSO-ANN results versus corresponding determined net present value (NPV) data, the PSO-ANN outputs lie over the line , the fact that indicates the identity of outputs gained from suggested PSO-ANN model and relevant net present value (NPV) data samples. The comparison between generated net present value (NPV) from the addressed approach and real net present value (NPV) data versus relevant data index has been shown in Figure 9. As illustrated in Figure 9, the obtained results of proposed model are as close as possible to net present value (NPV) data samples. To put it another way, the outputs of the PSO-ANN approach have the same behaviour as actual data do. Furthermore, the effectiveness of the proposed intelligent model has been depicted in terms of the relative deviations of PSO-ANN model outputs from corresponding indicated net present value (NPV) data in Figure 10. As can be seen from Figure 10, the highest deviations of the suggested approach results are subjected to the early boundary of net present value (NPV) data. 6% is the maximum degree of relative deviation depicted in Figure 10.

Figure 8: Performance plot of the suggested network model for determining net present value (NPV) of chemical flooding owing to correlation coefficient (): (a) training phase and (b) testing phase.
Figure 9: Comparison between suggested network model and net present value (NPV) versus relevant data index: (a) training phase and (b) testing phase.
Figure 10: Relative error distribution of the proposed approach versus actual net present value (NPV).

The performance efficiency of the selected network is assessed using the various error analysis parameters. Table 2 tabulates the PSO-ANN accuracy in terms of correlation coefficient (), coefficient of determination (), mean absolute error (MAE), and mean squared error (MSE), which are defined as follows: in which represents the total number of data points including either training, testing, or whole data set (input and output pairs), refers to the actual value at the sampling point , is the th output of the model, and and stand for the average magnitudes of the actual and predicted data, respectively.

Table 2: Statistical parameters of the proposed approaches in prediction of efficiency of chemical flooding in oil reservoirs.

5. Conclusions

Owing to the gained results of this contribution the following major conclusions can be drawn.(1)Adequate agreement between gain dew point pressure from the developed intelligent model and corresponding real recovery factor/net present value (NPV) values is observed. In other words, the conventional approaches fail to monitor real recovery factor/net present value (NPV) of chemical flooding dedicated to the gained statistical criteria such as mean square error (MSE) and correlation coefficient.(2)The evolved intelligent network model for monitoring real recovery factor/net present value (NPV) of chemical flooding is user friendly, fast, and cheap for implementation. Moreover, it is very useful and user friendly for evolving the accuracy and robustness of the commercial reservoir simulators like ECLIPSE and computer modelling group (CMG) software for enhanced oil recovery (EOR) from oil reservoirs.

Nomenclature

Abbreviations
:Artificial neural network
:Back propagation
:Genetic algorithm
:Hybrid genetic algorithm and particle swarm optimization
:Imperialist competitive algorithm
:Mean absolute error (MAE)
:Mean squared error (MSE)
:Particle swarm optimization
:Particle swarm optimization
:Coefficient of determination
:Unified particle swarm optimization.
Variables
:The average of the predicted data
:The average of the actual data
:Bias
:Cognition component
:Social components
:The total number of data points
: Output
: Two random numbers
: Sum of interconnection weights
: Velocity of particle
: Interconnection weights in network model
: Position of particle
: The th output of the model
: The actual at the sampling point .
Greek Letters
:Absolute relative error
:Activation function
:The inertia weight.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

  1. B. Shaker Shiran and A. Skauge, “Enhanced oil recovery (EOR) by combined low salinity water/polymer flooding,” Energy and Fuels, vol. 27, no. 3, pp. 1223–1235, 2013. View at Publisher · View at Google Scholar · View at Scopus
  2. S. Strand, T. Puntervold, and T. Austad, “Effect of temperature on enhanced oil recovery from mixed-wet chalk cores by spontaneous imbibition and forced displacement using seawater,” Energy and Fuels, vol. 22, no. 5, pp. 3222–3225, 2008. View at Publisher · View at Google Scholar · View at Scopus
  3. H. Zhang, M. Dong, and S. Zhao, “Which one is more important in chemical flooding for enhanced court heavy oil recovery, lowering interfacial tension or reducing water mobility?” Energy & Fuels, vol. 24, no. 3, pp. 1829–1836, 2010. View at Publisher · View at Google Scholar · View at Scopus
  4. A. Z. Abidin, T. Puspasari, and W. A. Nugroho, “Polymers for enhanced oil recovery technology,” Procedia Chemistry, vol. 4, pp. 11–16, 2012. View at Google Scholar
  5. M. A. Ahmadi, Y. Arabsahebi, S. R. Shadizadeh, and S. Shokrollahzadeh Behbahani, “Preliminary evaluation of mulberry leaf-derived surfactant on interfacial tension in an oil-aqueous system: EOR application,” Fuel, vol. 117, pp. 749–755, 2014. View at Publisher · View at Google Scholar · View at Scopus
  6. P. Luo, Y. Zhang, and S. Huang, “A promising chemical-augmented WAG process for enhanced heavy oil recovery,” Fuel, vol. 104, pp. 333–341, 2013. View at Publisher · View at Google Scholar · View at Scopus
  7. A. Mandal, A. Samanta, A. Bera, and K. Ojha, “Characterization of oil-water emulsion and its use in enhanced oil recovery,” Industrial and Engineering Chemistry Research, vol. 49, no. 24, pp. 12756–12761, 2010. View at Publisher · View at Google Scholar · View at Scopus
  8. A. Sabhapondit, A. Borthakur, and I. Haque, “Water soluble acrylamidomethyl propane sulfonate (AMPS) copolymer as an enhanced oil recovery chemical,” Energy and Fuels, vol. 17, no. 3, pp. 683–688, 2003. View at Publisher · View at Google Scholar · View at Scopus
  9. M. S. Karambeigi, R. Zabihi, and Z. Hekmat, “Neuro-simulation modeling of chemical flooding,” Journal of Petroleum Science and Engineering, vol. 78, no. 2, pp. 208–219, 2011. View at Publisher · View at Google Scholar · View at Scopus
  10. M. A. Ahmadi and S. R. Shadizadeh, “Implementation of a high-performance surfactant for enhanced oil recovery from carbonate reservoirs,” Journal of Petroleum Science and Engineering, vol. 110, pp. 66–73, 2013. View at Publisher · View at Google Scholar · View at Scopus
  11. G. Chen, X. Wang, Z. Liang et al., “Tontiwachwuthikul P: simulation of CO2-oil Minimum Miscibility Pressure (MMP) for CO2 Enhanced Oil Recovery (EOR) using neural networks,” Energy Procedia, vol. 37, pp. 6877–6884, 2013. View at Publisher · View at Google Scholar
  12. N. Loahardjo, X. Xie, and N. R. Morrow, “Oil recovery by sequential waterflooding of mixed-wet sandstone and limestone,” Energy & Fuels, vol. 24, no. 9, pp. 5073–5080, 2010. View at Publisher · View at Google Scholar · View at Scopus
  13. M. Ali Ahmadi and M. Golshadi, “Neural network based swarm concept for prediction asphaltene precipitation due to natural depletion,” Journal of Petroleum Science and Engineering, vol. 98-99, pp. 40–49, 2012. View at Publisher · View at Google Scholar · View at Scopus
  14. M. A. Ahmadi, “Neural network based unified particle swarm optimization for prediction of asphaltene precipitation,” Fluid Phase Equilibria, vol. 314, pp. 46–51, 2012. View at Publisher · View at Google Scholar · View at Scopus
  15. S. Zendehboudi, M. A. Ahmadi, L. James, and I. Chatzis, “Prediction of condensate-to-gas ratio for retrograde gas condensate reservoirs using artificial neural network with particle swarm optimization,” Energy & Fuels, vol. 26, no. 6, pp. 3432–3447, 2012. View at Publisher · View at Google Scholar · View at Scopus
  16. M. Ali Ahmadi, S. Zendehboudi, A. Lohi, A. Elkamel, and I. Chatzis, “Reservoir permeability prediction by neural networks combined with hybrid genetic algorithm and particle swarm optimization,” Geophysical Prospecting, vol. 61, no. 3, pp. 582–598, 2013. View at Publisher · View at Google Scholar · View at Scopus
  17. S. Zendehboudi, M. A. Ahmadi, O. Mohammadzadeh, A. Bahadori, and I. Chatzis, “Thermodynamic investigation of asphaltene precipitation during primary oil production: laboratory and smart technique,” Industrial and Engineering Chemistry Research, vol. 52, no. 17, pp. 6009–6031, 2013. View at Publisher · View at Google Scholar · View at Scopus
  18. S. Zendehboudi, M. A. Ahmadi, A. Bahadori, A. Shafiei, and T. Babadagli, “A developed smart technique to predict minimum miscible pressure—EOR implications,” Canadian Journal of Chemical Engineering, vol. 91, no. 7, pp. 1325–1337, 2013. View at Publisher · View at Google Scholar · View at Scopus
  19. M. A. Ahmadi and S. R. Shadizadeh, “New approach for prediction of asphaltene precipitation due to natural depletion by using evolutionary algorithm concept,” Fuel, vol. 102, pp. 716–723, 2012. View at Publisher · View at Google Scholar · View at Scopus
  20. M. A. Ahmadi, “Prediction of asphaltene precipitation using artificial neural network optimized by imperialist competitive algorithm,” Journal of Petroleum Exploration and Production Technology, vol. 1, no. 2–4, pp. 99–106, 2011. View at Publisher · View at Google Scholar · View at Scopus
  21. S. Zendehboudi, A. R. Rajabzadeh, A. Bahadori et al., “Connectionist model to estimate performance of steam-assisted gravity drainage in fractured and unfractured petroleum reservoirs: enhanced oil recovery implications,” Industrial and Engineering Chemistry Research, vol. 53, no. 4, pp. 1645–1662, 2014. View at Publisher · View at Google Scholar · View at Scopus
  22. M. A. Ahmadi, M. Ebadi, and S. M. Hosseini, “Prediction breakthrough time of water coning in the fractured reservoirs by implementing low parameter support vector machine approach,” Fuel, vol. 117, part A, pp. 579–589, 2014. View at Publisher · View at Google Scholar · View at Scopus
  23. M. A. Ahmadi, M. Ebadi, A. Shokrollahi, and S. M. J. Majidi, “Evolving artificial neural network and imperialist competitive algorithm for prediction oil flow rate of the reservoir,” Applied Soft Computing Journal, vol. 13, no. 2, pp. 1085–1098, 2013. View at Publisher · View at Google Scholar · View at Scopus
  24. M. A. Ahmadi and M. Ebadi, “Evolving smart approach for determination dew point pressure through condensate gas reservoirs,” Fuel, vol. 117, part B, pp. 1074–1084, 2014. View at Publisher · View at Google Scholar · View at Scopus
  25. M. A. Ahmadi, R. Soleimani, and A. Bahadori, “A computational intelligence scheme for prediction equilibrium water dew point of natural gas in TEG dehydration systems,” Fuel, vol. 137, pp. 145–154, 2014. View at Google Scholar
  26. M. A. Ahmadi, M. Ebadi, and A. Yazdanpanah, “Robust intelligent tool for estimating dew point pressure in retrograded condensate gas reservoirs: application of particle swarm optimization,” Journal of Petroleum Science and Engineering, 2014. View at Publisher · View at Google Scholar
  27. J. Kennedy, “Particle swarm: social adaptation of knowledge,” in Proceedings of the IEEE International Conference on Evolutionary Computation (ICEC '97), pp. 303–308, Indianapolis, Ind, USA, April 1997. View at Scopus
  28. M.-A. Ahmadi, M. Masumi, R. Kharrat, and A. H. Mohammadi, “Gas analysis by in situ combustion in heavy-oil recovery process: experimental and modeling studies,” Chemical Engineering & Technology, vol. 37, no. 3, pp. 409–418, 2014. View at Publisher · View at Google Scholar · View at Scopus
  29. G. V. Cybenko, “Approximation by superpositions of a sigmoidal function,” Mathematics of Control, Signals, and Systems, vol. 2, no. 4, pp. 303–314, 1989. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus