ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ ПРОЦЕССОМ КАТАЛИТИЧЕСКОГО РИФОРМИНГА БЕНЗИНОВЫХ ФРАКЦИЙ
Авторы:
1.
Волго-Каспийский морской рыбопромышленный колледж ФГБОУ ВО «Астраханский государственный технический университет»; Астраханский государственный университет имени В. Н. Татищева
(кандидат технических наук; преподаватель механического отделения; магистрант факультета физики, математики и инженерных технологий)
2.
Астраханский государственный технический университет
Тип:
Статья
Страницы:
с 34
по 42
Статус:
Опубликован
Получено:
05.04.2018
Одобрено:
25.04.2018
Опубликовано:
25.04.2018
Классификаторы:
ВАК 05.12.2013 Системы, сети и устройства телекоммуникаций
ВАК 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
ВАК 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
ВАК 05.13.10 Управление в социальных и экономических системах
ВАК 05.13.18 Математическое моделирование, численные методы и комплексы программ
ВАК 05.13.19 Методы и системы защиты информации, информационная безопасность
УДК 681.5
ГРНТИ 20.01 Общие вопросы информатики
ГРНТИ 28.01 Общие вопросы кибернетики
ГРНТИ 49.01 Общие вопросы связи
ГРНТИ 50.01 Общие вопросы автоматики и вычислительной техники
ГРНТИ 82.01 Общие вопросы организации и управления
ВАК 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
ВАК 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
ВАК 05.13.10 Управление в социальных и экономических системах
ВАК 05.13.18 Математическое моделирование, численные методы и комплексы программ
ВАК 05.13.19 Методы и системы защиты информации, информационная безопасность
УДК 681.5
ГРНТИ 20.01 Общие вопросы информатики
ГРНТИ 28.01 Общие вопросы кибернетики
ГРНТИ 49.01 Общие вопросы связи
ГРНТИ 50.01 Общие вопросы автоматики и вычислительной техники
ГРНТИ 82.01 Общие вопросы организации и управления
Язык материала:
английский
Ключевые слова:
каталитический риформинг, обобщенный критерий оптимальности, Парето-оптимальное решение, гибридная модель, классический ПИД-регулятор, нечеткий ПИД-регулятор, система оптимального управления
Аннотация (русский):
Рассматривается актуальная задача разработки системы оптимального управления процессом каталитического риформинга, обеспечивающей достижение оптимальных соотношений октанового числа бензина и производственных затрат. Выполнена постановка задачи оптимального управления процессом каталитического риформинга, отличающаяся использованием обобщенного критерия оптимальности в качестве целевого функционала. Разработана гибридная математическая модель процесса каталитического риформинга, которая учитывает влияние параметров, характеризующих производственные затраты и октановое число бензина, на эффективность процесса и позволяет вычислять значения обобщенного критерия оптимальности в зависимости от значений входных переменных процесса. На основе модели для установки ЛЧ-35-11/1000 выполен расчет максимальных годовых затрат (4,05 млрд руб.) и минимального октанового числа бензина 92,83. Выполнен многокритериальный выбор системы регулирования температуры на выходе из печи. Синтезирован алгоритм определения настроек регулятора температуры на выходе из печи в различных режимах процесса каталитического риформинга. Разработана методика оптимального управления процессом каталитического риформинга, учитывающая экспертную информацию при формализации нечетких целей и ограничений при управлении процессом и позволяющая рассчитывать оптимальные управляющие воздействия по схеме Беллмана-Заде в зависимости от значений входных переменных процесса. На основе данной методики определен вектор управлений, соответствующий минимуму обобщенного критерия оптимальности ( I = 0,964). Выполнен поиск Парето-оптимальных решений при управлении процессом каталитического риформинга. Получены положительные эффекты при использовании разработанной системы оптимального управления процессом каталитического риформинга: снижение средних затрат на 0,33 млн руб.; повышение среднего октанового числа на 0,53; снижение среднего значения обобщенного критерия оптимальности на 0,025.
Рассматривается актуальная задача разработки системы оптимального управления процессом каталитического риформинга, обеспечивающей достижение оптимальных соотношений октанового числа бензина и производственных затрат. Выполнена постановка задачи оптимального управления процессом каталитического риформинга, отличающаяся использованием обобщенного критерия оптимальности в качестве целевого функционала. Разработана гибридная математическая модель процесса каталитического риформинга, которая учитывает влияние параметров, характеризующих производственные затраты и октановое число бензина, на эффективность процесса и позволяет вычислять значения обобщенного критерия оптимальности в зависимости от значений входных переменных процесса. На основе модели для установки ЛЧ-35-11/1000 выполен расчет максимальных годовых затрат (4,05 млрд руб.) и минимального октанового числа бензина 92,83. Выполнен многокритериальный выбор системы регулирования температуры на выходе из печи. Синтезирован алгоритм определения настроек регулятора температуры на выходе из печи в различных режимах процесса каталитического риформинга. Разработана методика оптимального управления процессом каталитического риформинга, учитывающая экспертную информацию при формализации нечетких целей и ограничений при управлении процессом и позволяющая рассчитывать оптимальные управляющие воздействия по схеме Беллмана-Заде в зависимости от значений входных переменных процесса. На основе данной методики определен вектор управлений, соответствующий минимуму обобщенного критерия оптимальности ( I = 0,964). Выполнен поиск Парето-оптимальных решений при управлении процессом каталитического риформинга. Получены положительные эффекты при использовании разработанной системы оптимального управления процессом каталитического риформинга: снижение средних затрат на 0,33 млн руб.; повышение среднего октанового числа на 0,53; снижение среднего значения обобщенного критерия оптимальности на 0,025.
Ключевые слова:
каталитический риформинг, обобщенный критерий оптимальности, Парето-оптимальное решение, гибридная модель, классический ПИД-регулятор, нечеткий ПИД-регулятор, система оптимального управления
каталитический риформинг, обобщенный критерий оптимальности, Парето-оптимальное решение, гибридная модель, классический ПИД-регулятор, нечеткий ПИД-регулятор, система оптимального управления
Introduction Catalytic reforming (CR) is the most important technological process of modern oil refining and petrochemistry. The CR process is used to produce high-octane gasoline, aromatic hydrocarbons and hydrogen. The increasing urgency of improving the process of CR is determined by the growing demand for high-octane motor fuel. As of 2016 the share of reforming processes in the total number of petroleum refining processes in Eastern Europe was 13.4% [1]. At present, the process capacity of CR in Russia is 9.3% from the total capacity of primary oil refining. As the component composition of gasoline contains 54.1% of reformate, the process of CR takes a leading place at the domestic refineries [2]. In existing studies on the management and modeling of the process of CR, the task of improving the efficiency of this process is solved by developing automatic control systems of technological parameters, as well as optimal control systems by basic technical and economic indicators (income, octane number, etc.) [3]. Since the efficiency of gasoline production depends not only on increasing the octane number of produced fuel, but also on reducing its cost [3], it is necessary to build a system for managing the process of CR ensuring the achievement of optimal ratios of the octane number of gasoline and production costs. In the present study, as the target functional for the optimal control of the process of CR, it was proposed to use the generalized optimality criterion (GOC), which includes indicators characterizing the octane number of gasoline and production costs. A great opportunity for taking into account the influence of these indicators on the efficiency of the process of CR is the construction of a hybrid mathematical model (MM) of this process [4]. It is necessary to develop a system of optimal control of the process of CR using a hybrid MM. In order to determine the optimal ratios of particular control criteria of the process of CR, such as the octane number of gasoline and production costs, this paper suggests the use of the well-known search method of Pareto-optimal solutions [5]. Thus, the development of the system of optimal control of the process of CR, ensuring the achievement of Pareto-optimal ratios of the octane number of gasoline and production costs, is an actual scientific and practical task. Formulation of the problem The aim of the study is to increase the effectiveness of the management of the process of CR by developing a system of the optimal management of this process, ensuring the achievement of Pareto-optimal ratios of the octane number of gasoline and production costs. In accordance with the purpose of the research, the formulation of the problem of optimal control of the CR process has completed, which is distinguished by the use of GOC as the objective functional: , (1) where ON0 - minimum value of octane number ON; Z0 - maximum cost value Z; I1 - a criterion inverse to the criterion of maximization of the average normalized value ON0 of the octane number ON at the time interval [0, T]; I2 - a minimization criterion of the average normalized value Z0 of value of costs Z at the time interval [0, T]; k1, k2 - weight coefficients that regulate the relative importance of the criteria I1, I2 with regard to the adoption of Pareto-optimal solutions, k1 + k2 = 1; 0 < k1 < 1; 0 < k2 < 1; T - the operating time of the installation of the CR, including inter-regeneration and inter-repair intervals. In a general form, the problem of optimal control of the CR process is formulated. With the given input variables Xin find the control actions U that ensure the minimum of the GOC (1) for superimposed connections in the form of a hybrid MM process of the CR and limitations where V - set of values of control actions U; W - set of values of variables of the state of the process A. To achieve this goal, it is necessary to solve the following tasks: 1. Development of a hybrid mathematical model of the process of CR. 2. Multi-criteria selection of automatic control systems for technological parameters of the process. 3. Synthesis of the algorithm for determining the temperature controller settings at the exit from the reforming furnace in various modes. 4. Development of methods for optimal control of the process of CR. 5. Search for Pareto-optimal solutions in the control of the CR process. 6. Development of the optimal management system of the CR process. Hybrid mathematical model of the catalytic reforming process A hybrid MM of the CR process was developed [6], which is generally represented (2) where Xout - a vector of output variables, including the octane number of gasoline ON and production costs Z; F - hybrid connection function between the variables of the hybrid MM of the CR process. The scientific novelty of the hybrid MM process of the CR is to take into account the influence of the parameters characterizing the production costs and the octane number of gasoline on the efficiency of the process and the calculation of the values of the generalized optimality criterion, depending on the current values of the input variables of the process. Fig. 1 shows the graphs of the change in the normalized criterion for minimizing the production costs I1(t), the normalized criterion I2(t), the inverse of criterion of minimization of the octane number and GOC I(t) respectively. Fig. 1. Change in the normalized criteria for the optimality of the catalytic reforming process According to (1) the definition of GOC I(t) for a hybrid MM (2) is the finding a weighted sum of normalized partial criteria I1(t), I2(t). The graphs are made on the basis of a hybrid MM for a typical set of input variables, characterizing the normal operation of the process at the LCH-35-11/1000 CR installation for one month [7]. Verification of the adequacy of the hybrid MM was performed using the Fisher criterion. The maximum error in calculations for the hybrid MM relative to the results of the experiments does not exceed 5%. Multi-criteria selection of automatic control systems of technological parameters of the catalytic reforming process Based on the method of fuzzy pair comparisons, mathematical expressions for the multi-criteria choice of automatic control systems (ACS) of the temperature at the exit from the furnace of 5 alternatives according to 7 criteria were developed [8]. A fuzzy solution on the scheme of Bellman-Zadeh is determined by crossing the particular criteria , where - the function of belonging (FB) of the element ki to the fuzzy set X; «~» - designation of fuzzy dimension. As a result of calculation by expressions, it was found that ACS with fuzzy controller of PID type MT20-V-E FOTEK + PID Fuzzy, corresponding to the maximum FB of fuzzy solution equal to 0.954, meets all the criteria most fully. The adequacy of these expressions was checked using the Fisher test. Software has been developed for multi-criteria selection of CR process control tools in the C # language, which allows choosing the best ACS temperature at the furnace outlet from the reformer from 5 alternatives according to 7 criteria. Using the software and the database, ACS temperature was chosen at the furnace output with a fuzzy PID controller of the type MT20-R-E MT-20E FOTEK + PID Fuzzy. Algorithm for determining the temperature regulator settings at the furnace outlet from the reformer in various modes The scientific novelty of the algorithm is to take into account the regime parameters in the transfer function of the disturbance quantitatively and in the determination of the temperature regulator settings at the outlet from the reforming furnace, depending on the current regime parameters. The heating process in the reforming furnace is considered as a control object (CO). The determination of the temperature controller settings at the furnace outlet is carried out under the action of CO disturbances - temperature difference in the furnace ΔT [9]. To determine the settings of the fuzzy temperature PID controller at the furnace outlet, a mathematical description of CO and perturbations in the form of transfer functions has been performed, and the PID controller settings in the MATLAB system have been determined. As a result of comparing the quality of transients in the ACS with the classical and fuzzy PID regulators (Fig. 2), positive effects were obtained from using a fuzzy PID controller: the reduction of overshoot by 10.4 times; the decrease of the quadratic integral criterion by 4.3 times. Fig. 2. Transients in ACS temperature at the outlet of the furnace: 1 - with a classic PID controller; 2 - with fuzzy PID controller Based on the algorithm for determining the settings of the temperature regulator at the furnace outlet in different modes of the CR process, implemented in the MATLAB system, the settings of the fuzzy PID controller MT20-R-E MT-20E FOTEK + PID Fuzzy are defined: FB type gaussmf; Mamdani's fuzzy inference system; P = 0.345, I = 0.017, D = 1.725, ensuring the achievement of the best indicators of the quality of the transient: overshoot is 0.005°С and the quadratic integral criterion is 1216. The method of optimal control of the catalytic reforming process The solution to the problem of optimizing the CR process using the method of uniform-dimensional search was substantiated. To do this, the dependence of the GOC I on two variables - the temperature of the feed mixture at the entrance to the furnace Tin and the temperature of the feed mixture at the furnace outlet Tout is considered [10]. The optimization algorithm is implemented in the form of software C #. As a result of the calculation by the algorithm, the following optimum values of the parameters of the CR process were obtained: the temperature at the entrance to the furnace Tinopt = 180°С; the temperature at the furnace outlet Toutop t= 210°С; GOC Iopt = 0.964. Based on the algorithm of optimizing the process of the CR, the method for the optimal control of this process has been developed. Scientific novelty of the method of optimal control of the process of the CR is in taking into account the expert information when formalizing fuzzy goals and constraints in the process control and determining the optimal control actions under the scheme of Bellman-Zadeh, depending on the current values of the input variables of the process. This technique is intended for optimal control of the CR process in the case of an undefined goal and fuzzy constraints that are specified by the operator in the form of verbal utterances (formulations) and is an extension (addition) of the algorithm for optimizing the CR process in the case of using verbal formulations in determining the purpose and constraints on management. A fuzzy decision in the management of the process of the CR is determined by the Bellman-Zadeh scheme [11] as the intersection of a fuzzy goal and fuzzy constraints : . (3) On the basis of (3) the alternative x = 212°С is defined, corresponding to the maximum FB of the fuzzy solution. Using the hybrid MM for given input variables, we obtain optimal control actions at x = 212°С for setting the LCH-35-11/1000. On the basis of the technique, the vector of controls corresponding to the minimum of the GOC I = 0.964 is determined. The vector of Uz controls that affect production costs is allocated: raw material consumption Gr 1000 t, electricity consumption Ge 82 MJ, fuel gas consumption Gfg 25 m3, catalyst consumption Gc 0.02 kg, reagent consumption Grg 0.001 kg. The vector of Uon controls influencing the octane number of gasoline is also determined: the productivity of the centrifugal compressor is Qcc 210 000 nm3/h, the volume flow of raw materials is Qr 130 m3/h, the amount of the discharged hydrogen-containing gas is Qhg 90 000 nm3, the volume flow of fuel gas to the first furnace is Qfg1 750 m3/h, the volume flow of fuel gas to the second furnace is Qfg2 850 m3/h, the volume flow of fuel gas to the third furnace is Qfg3 950 m3/h. The program implementation of the optimal control methodology for the process of the CR in the C # language is made, which allows determining the vector of controls corresponding to the minimum of the GOC. For the obtained control vector, the efficiency indicators of the process of the CR were determined: production costs are Z = 0.825 billion rubles; octane number of gasoline is ON = 92.83; GOC is I = 0.964. Search for Pareto-optimal solutions for controlling the catalytic reforming process The use of normalized values of the criteria is associated with the possibility of combining them on one criteria plane and searching for Pareto-optimal solutions (Fig. 3). Fig. 3. Search for the optimal solution on the Pareto set The values of the partial criteria and GOC are determined using the hybrid MM process of the CR [6]. The optimization algorithm is implemented in the form of software C #. As a result of the calculation by the algorithm, the following optimum values of the parameters of the CR process were obtained: the temperature at the entrance to the furnace Tinopt = 180°С; the temperature at the furnace outlet Toutopt = 210°С; GOC Iopt = 0.964. The Pareto set (Fig. 3) is constructed on the basis of the method of optimal control of the CR process for a typical set of input variables characterizing the normal operation of the process at the LCH-35-11/1000 installation for one month [7]. Based on Fig. 3, the best solution corresponding to the minimum of GOC is determined: , (4) where - the values of the partial criteria for the 1st (best) solution. The optimal solution (4) corresponds to the values of the criteria: I1 = 0.959; I2 = 0.967; I = 0.964 at weighting factors k1 = k2 = 0.5. Development of the system for optimal control of the catalytic reforming process The system of optimal control of the process of the CR is developed, the structure of which is presented in Fig. 4. Fig. 4. Structure of the optimal control system for the catalytic reforming process The control object (CO) - 1 is the CR process, which receives the input variables XIN and control actions U. At the output of CO the values of the octane number of gasoline ON and production costs Z are obtained. The optimal control system - 2 serves to control the CR process with a fuzzy goal and the fuzzy constraints specified by operator - 3 when implementing the optimal control technique - 4. To determine GOC - 5, the operator - 3 sets the values of the fuzzy input variables X*IN (catalyst activity, state of the reforming furnace, etc.) to the input of the hybrid MM process - 6. On the basis of the algorithm for optimizing the process - 7, with the given explicit constraints, the alternative "x" is determined - the temperature at the exit from the furnace, at which GOC I(x) is minimal. Based on the method of optimal control process - 4 using the Bellman-Zadeh scheme, an alternative to "x’" is defined, different from "x", at which the minimum of GOC I(x) is reached. For the obtained alternative "x’" and for alternatives from the range, when all components are equal the GOC is minimal with respect to two other local minima, using the search procedure for the Pareto optimal solution - 8 and the hybrid MM - 6, the control vector U - 9 corresponding to the minimum GOC is determined. Based on the multi-criteria selection of ACS - 10 and the algorithm for determining the settings of the CR process controllers in various modes - 11 the best model of ACS is determined and regulator settings for the three parameters of the CR process are the temperature at the outlet from the reforming furnace (ACS1); pressure in the reforming reactor (ACS2); level in the re-forming separator (ACS3). To determine the initial information in the selection of ACS, the database of tools for managing the CR process was developed. Positive effects were obtained using the developed system of optimal control of the CR process: a decrease in average costs by 0.33 million rubles; an increase in the average octane number by 0.53; decrease in the mean value of the GOC by 0.025. Conclusion 1. The problem of optimal control of the CR process has been formulated. 2. The hybrid MM of the CR process has been developed, based on which the maximum annual costs of 4.05 billion rubles and the minimum octane number of gasoline is 92.83 were obtained for the LCH-35-11/1000 CR installation. 3. The algorithm for determining the settings of the temperature regulator at the outlet from the furnace in different modes of the CR process has been synthesized. Based on the algorithm for the given mode of the CR process, the parameters of the fuzzy PID controller MT20-R-E MT-20E FOTEK + PID Fuzzy have been determined, which ensure the achievement of the best parameters of the transient: overshoot 0.005°С and integral quality index 1216. 4. A method for optimal control of the CR process has been developed, on the basis of which a control vector corresponding to the minimum of the GOC I = 0.964 has been determined. 5. A set of Pareto-optimal solutions has been constructed for a typical set of input variables characterizing the normal operation of the process at the LCH-35-11/1000 CR installation for one month. On the Pareto set, the best solution corresponding to the values of the criteria I1 = 0.959; I2 = 0.967; I = 0.964 is defined at weighting factors k1 = k2 = 0.5. 6. The structure of the system of optimal control of the CR process has been developed. 7. The assessment of improving the management effectiveness of the CD process has been made. Positive effects have been determined while using the developed system of optimal control of the CR process: the decrease in average costs by 0.33 million rubles; the increase in the average octane number by 0.53 points; the decrease in the mean value of the GOC by 0.025.
1. Imashev U. B., Tiurin A. A., Udalova E. A. Osobennosti razvitiia protsessa kataliticheskogo riforminga v Rossii [Peculiarities of the development of the catalytic reforming process in Russia]. Bashkirskii khimicheskii zhurnal, 2009, vol. 16, no. 4, pp. 184-186.
2. Sharipov R. A., Sidorov G. M., Zinnatullin R. R., Dmitriev Iu. K. Rol' protsessa kataliticheskogo krekinga v proizvodstve vysokooktanovykh avtomobil'nykh benzinov [The role of the process of catalytic cracking in the production of high-octane gasolines]. Sovremennye problemy nauki i obrazovaniia, 2015, no. 1-1, p. 134.
3. Ivanchina E. D., Deriglazov V. V., Zanin I. K. Povyshenie tekhniko-ekonomicheskoi effektivnosti kataliticheskogo riforminga s ispol'zovaniem komp'iuternoi modeliruiushchei sistemy [Increase of technical and economic efficiency of catalytic reforming using a computer modeling system]. Izvestiia Tomskogo politekhnicheskogo universiteta. Inzhiniring georesursov, 2011, vol. 319, no. 3, pp. 105-109.
4. Protalinskii O. M. Proverka dostovernosti pervichnoi informatsii v ASU TP s ispol'zovaniem nechetkikh mnozhestv [Verification of the reliability of primary information in the process control system using fuzzy sets]. Izvestiia vysshikh uchebnykh zavedenii. Severo-Kavkazskii region. Seriia: Tekhnicheskie nauki, 2003, no. 3, p. 60.
5. Litvinskaia O. S. Formalizatsiia priniatiia reshenii na osnove tselevogo funktsionala [Formalization of decision-making on the basis of the target functional]. XXI vek: itogi proshlogo i problemy nastoiashchego plius, 2012, no. 1 (5), pp. 74-79.
6. Dzhambekov A. M. Upravlenie protsessom kataliticheskogo riforminga na osnove ekspertnoi informatsii [Controlling the process of catalytic reforming on the basis of expert information]. Novye materialy i tekhnologii: sostoianie voprosa i perspektivy razvitiia: sbornik materialov Vserossiiskoi molodezhnoi nauchnoi konferentsii (Saratov, 24-26 iiunia 2014 g.). Saratov, OOO «Izdat. tsentr «Nauka», 2014. Pp. 382-387.
7. Vlasov S. S., Shumikhin A. G. Modelirovanie protsessa otbenzinivaniia nefti pri prognozirovanii pokazatelei kachestva benzina [Modeling the process of oil stripping while predicting the quality of gasoline]. Vestnik Saratovskogo gosudarstvennogo tekhnicheskogo universiteta, 2012, vol. 1, no. 1 (63), pp. 90-94.
8. Dzhambekov A. M. Mnogokriterial'nyi vybor sistem avtomaticheskogo regulirovaniia tekhnologicheskikh parametrov v usloviiakh neopredelennosti [Multicriteria choice of systems of automatic control of technological parameters under conditions of uncertainty]. Nauchnye issledovaniia i perspektivnye proekty - 2016: sbornik trudov I Nauchno-prakticheskoi konferentsii aspirantov, prepodavatelei, uchenykh (Ufa, 26-29 aprelia 2016 g.). Ufa, OOO «Klaud Kepital» Publ., 2016. Pp. 95-105.
9. Dzhambekov A. M. Nechetkaia sistema upravleniia protsessom kataliticheskogo riforminga [Fuzzy control system of the catalytic reforming process]. Prikaspiiskii zhurnal: upravlenie i vysokie tekhnologii, 2015, no. 4 (32), pp. 268-280.
10. Vorobovich N. P. Primenenie teorii optimal'nogo upravleniia dlia resheniia mnogokriterial'noi zadachi setevogo planirovaniia s nechetkimi ogranicheniiami po resursam [Application of the optimal control theory for the solution of the multicriteria network planning problem with fuzzy resource constraints]. Vestnik Krasnoiarskogo gosudarstvennogo agrarnogo universiteta, 2009, no. 8, pp. 15-18.
11. Shumikhin A. G., Vialykh A. G. Formirovanie funktsii prinadlezhnosti dlia algoritma nechetkogo upravleniia tekhnologicheskim protsessom kataliticheskogo krekinga [Formation of membership functions for the algorithm of fuzzy control of the technological process of catalytic cracking]. Izvestiia Tomskogo politekhnicheskogo universiteta. Inzhiniring georesursov, 2010, vol. 316, no. 5, pp. 132-136.
