Study of Kinetic coefficients of a Membrane Bioreactor (MBR) for municipal wastewater treatment

Similar documents
Application of MBR Technology in Municipal Wastewater Treatment

OPTIMIZATION OF AN INTERMITTENTLY AERATED AND FED SUBMERGED MEMBRANE BIOREACTOR

Characteristics of Nutrient Removal in Vertical Membrane Bioreactors

Module 19 : Aerobic Secondary Treatment Of Wastewater. Lecture 26 : Aerobic Secondary Treatment Of Wastewater (Contd.)

operation of continuous and batch reactors. Contrary to what happens in the batch reactor, the substrate (BOD) of the wastewater in the continuous rea

Effect of organic load of refinery wastewater in municipal wastewater treatment using bench scale extended aeration activated sludge system

Application of MBR for the treatment of textile wastewater

A Study on Sludge Reduction in Sewage using Microbial Catalysts

Environmental Engineering I Jagadish Torlapati, PhD Fall 2017 MODULE 3 WASTEWATER TREATMENT CONTROL PARAMETERS QS = VX F M

Impact of Petroleum Refinery Wastewater on Activated Sludge

1/11/2016. Types and Characteristics of Microorganisms. Topic VI: Biological Treatment Processes. Learning Objectives:

Effect of the start-up length on the biological nutrient removal process

Permeate Characteristics and Operation Conditions of a Membrane Bioreactor with Dispersing Modified Hollow Fibers

membrane bioreactors MBR vs. Conventional Activated Sludge Systems CAS

Iranica Journal of Energy & Environment

COMPARISON STUDY BETWEEN INTEGRATED FIXED FILM ACTIVATED SLUDGE (IFAS), MEMBRANE BIOREACTOR (MBR) AND CONVENTIONAL ACTIVATED SLUDGE (AS) PROCESSES

CSR Process Simulations Can Help Municipalities Meet Stringent Nutrient Removal Requirements

ASSESSMENT OF AN MBR SYSTEM FOR SEGREGATED HOUSEHOLD WASTEWATER BY USING SIMULATION MATHEMATICAL MODEL

ANOXIC BIOREACTOR SIZING

The optimum operational condition of membrane bioreactor (MBR): cost estimation of aeration and sludge treatment

Designing Single-Sludge Bionutrient Removal Systems

Biochemical Kinetics of Cross flow Membrane Bioreactor Processes in the Treatment of Refinery Wastewater

HOW TO SELECT SRT. SRT is selected based on the following treatment requirements/objectives:

Mathematical Modelling of Wastewater Treatment Plant of Žiar nad Hronom

Wastewater Tools: Activated Sludge and Energy Use Analysis

et al., 2005). Therefore, in the MBR systems, the DO concentration can easily fluctuate above 4 mg O 2

Comnletely Mixed Activated SIud2e (CMAS) Bioreactor Desi2n Eguations

Wastewater Treatment by Anaerobic Digestion Coupled with Membrane Processing

SWIM and Horizon 2020 Support Mechanism

1. Introduction. 2. Material and methods. Soheil Farajzadehha 1+, Jalal Shayegan 2, S.A.Mirbagheri 1, Soroush Farajzadehha 4

Membrane BioReactor: Technology for Waste Water Reclamation

Removal of High C and N Contents in Synthetic Wastewater Using Internal Circulation of Anaerobic and Anoxic/Oxic Activated Sludge Processes

EFFECT OF DYE CONCENTRATION ON SEQUENCING BATCH REACTOR PERFORMANCE

DEGRADATION OF AMMONIA IN AN INTEGRATED ANOXIC-AEROBIC CLARIFIER SYSTEM

Kinetic Modelling For Rotating Biological Contactor for Treatment of Greywastewater

Secondary Wastewater Treatment

Use of Sequencing Batch Reactors (SBRs) in Treatment of Wood Fiber Wastewater

Bio-Tiger Model Biokinetic Model that Simulates Activated Sludge Processes

PERFORMANCE EVALUATION OF BENCH SCALE MEMBRANE BIOREACTOR (MBR) SYSTEMS

Pilot Studies on Performance of Membrane Bio-Reactor in Treating Hong Kong Freshwater and Saline Sewage and Its Virus Rejection Ability and Mechanism

SIMULTANEOUS NITRIFICATION DENITRIFICATION IN AN AEROBIC MBR TREATING HIGH STRENGTH PET FOOD WASTEWATER Chandresh Acharya 1, G. Nakhla 2, A.

19. AEROBIC SECONDARY TREATMENT OF WASTEWATER

Performance of a Novel Hybrid Membrane Bioreactor: Effect of Bacterial Floc Size on Fouling

ENHANCED BIOLOGICAL PHOSPHORUS REMOVAL WITHIN MEMBRANE BIOREACTORS. 255 Consumers Road Toronto, ON, Canada, M2J 5B6

New Developments in BioWin 5.3

Optimization of a Combined UASB and Continuous-flow SBR System for Sludge Reduction and Biogas Production

PHENOLIC WASTEWATER TREATMENT

PERFORMANCE OF AERATED SUBMERGED FIXED-FILM BIOREACTOR FOR TREATMENT OF ACRYLONITRILE-CONTAINING WASTEWATER

Determination of Kinetic Coefficients for Secondary Aerobic Treatment of Digested Spentwash based on COD Analysis

a new membrane-aerated biofilm reactor for low energy wastewater treatment: pilot results

MODELING OF BIOOLOGICAL PHOSPHORUS REMOVAL USING MEMBRANE BIOREACTOR: PART II

Figure Trickling Filter

Probabilistic Modeling of Two-Stage Biological Nitrogen Removal Process: Formulation of Control Strategy for Enhanced Process Certainty

A mini-review of modeling studies on membrane bioreactor (MBR) treatment for municipal wastewaters

Key Points. The Importance of MCRT/SRT for Activated Sludge Control. Other (Confusing) Definitions. Definitions of SRT

Innovating Water Treatment Technologies in Managing Water-Energy Demand

Modelling of Wastewater Treatment Plants

Troubleshooting Activated Sludge Processes. PNCWA - Southeast Idaho Operators Section Pocatello, ID February 11, 2016 Jim Goodley, P.E.

BIOLOGICAL PHOSPHOROUS REMOVAL AN OPERATOR S GUIDE

Advances in wastewater treatment by combined microbial fuel cell-membrane bioreactor

REMOVAL OF THE HIGH AMMONIA NITROGEN CONCENTRATION AT THE DIFFERENT SLUDGE AGES IN THE MEMBRANE BIOREACTOR

COMPARING THE EFFICIENCY OF UAFF AND UASB WITH HYBRID REACTOR IN TREATING WOOD FIBER WASTEWATER

Department of Civil Engineering-I.I.T. Delhi CVL723 Problem Set_2_Feb6_15

New Developments in BioWin 4.0

Petroleum Refinery Effluent Biodegradation in Sequencing Batch Reactor

Energy Efficient Wastewater Treatment. Emily Gonthier & Jennifer Lawrence April 30, 2018

UASB. Performance of the UASB Pilot Scale Reactor in Treating MSW Leachate in the Hot Climate City of Ahvaz

Energy and chemical efficient nitrogen removal at a full-scale MBR water reuse facility

Online at International Journal of Current Innovation Research Issue, 6(A), pp , June, 2018.

EFFECT OF FERROUS SULPHATE ON THE SIMULTANEOUS ORGANIC MATTER AND NUTRIENT REMOVAL PERFORMANCE OF SEQUENCING BATCH REACTOR

Journal of Separation Science and Engineering Vol. 3, No. 1, pp (MBR)

Asian Institute of Technology School of Environment, Resources and Development Bangkok, Thailand August, 1996

Study on nitrification and denitrification of high nitrogen and COD load wastewater in moving bed biofilm reactor

USE OF A ROTATING BIOLOGICAL CONTACTOR FOR APPROPRIATE TECHNOLOGY WASTEWATER TREATMENT

Investigation of Excess Biological Sludge Reduction in Sequencing Bach Reactor

HOSPITAL WASTEWATER TREATMENT USING MIXED MEDIA BIOFILTER

SIMULATING THE EFFECTS OF INTRODUCTION OF MEMBRANE FILTRATION AT A MUNICIPAL BNR PLANT

Evaluation of Conventional Activated Sludge Compared to Membrane Bioreactors

ENVE 302 Environmental Engineering Unit Processes DENITRIFICATION

DEVELOPMENT OF THE. Ken Mikkelson, Ph.D. Ed Lang Lloyd Johnson, P.E. Aqua Aerobic Systems, Inc.

USING NUMERICAL SIMULATION SOFTWARE FOR IMPROVING WASTEWATER TREATMENT EFFICIENCY

Treatment of Swine Wastewater using Sequencing Batch Reactor*

Aqua-Aerobic MBR Topics

COD Reduction by Moving Bed Biofilm Reactor

Ammonia Removal using Sequencing Batch Reactor: The Effects of Organic Loading Rate

Chapter 9 Nitrification

Operation of a small scale MBR system for wastewater reuse

Submerged Membrane Technology for MBR

Water and Wastewater Engineering Dr. Ligy Philip Department of Civil Engineering Indian Institute of Technology, Madras

Design of Anaerobic Treatment/Digestion Processes

THE INFLUENCE OF CARRIER SIZE AND SHAPE IN THE MOVING BED BIOFILM PROCESS

ENHANCING THE PERFORMANCE OF OXIDATION DITCHES. Larry W. Moore, Ph.D., P.E., DEE Professor of Environmental Engineering The University of Memphis

Pradyut Kundu 1 Lecturer, Department of Food Processing Technology Acharyya Prafulla Chandra Ray Polytechnic Jadavpur, Kolkata, India

Factors Affecting Filtration Characteristics in Submerged Membrane Bioreactor for Wastewater Treatment

CEE 370 Fall Homework #8

Aerobic digestion efficiency for mixed sludge (primary and secondary) of West Ahvaz wastewater Treatment plant in Batch reactor

Appendix D JWPCP Background and NDN

Membrane Bio-Reactors (MBRs) The Future of Wastewater Technology, Science and Economy Aspects

NEW BIOLOGICAL PHOSPHORUS REMOVAL CONCEPT SUCCESSFULLY APPLIED IN A T-DITCH PROCESS WASTEWATER TREATMENT PLANT

Activated Sludge Processes

Transcription:

RESEARCH ARTICLE Arch Hyg Sci 213;2(3):18-113 Journal Homepage: http://jhygiene.muq.ac.ir Study of Kinetic coefficients of a Membrane Bioreactor (MBR) for municipal wastewater treatment Ali Naghizadeh a *, Elham Derakhshani b, Ahmad Reza Yari c a Department of Environmental Health Engineering, School of Public Health, Birjand University of Medical Sciences, Birjand, Iran. b Department of Environmental Health Engineering, School of Public Health, Shahid Beheshti University of Medical Sciences, Tehran, Iran. c Research Center for Environmental Pollutants, Qom University of Medical Sciences, Qom, Iran. * Correspondence should be addressed to Dr. Ali Naghizadeh; Email: al.naghizadeh@yahoo.com A-R-T-I-C-L-E I-N-F-O Article Notes: Received: Oct 3, 213 Received in revised form: November 12, 213 Accepted: November 16, 213 Available Online: November 22, 213 Keywords: Biokinetic coefficient Membrane bioreactor municipal wastewater A-B-S-T-R-A-C-T Background & Aims of the Study: In order to design membrane bioreactors (MBR) properly, it is essential to comprehend the behavior of microorganisms in such wastewater treatment processes. Materials & Methods: In this study, a lab-scale MBR process was operated to determine the biokinetic coefficients of the MBR system under different MLSS concentrations of 68, 7, 74, and 78 mg/l and organic loading rates of.5 kg /m 3 /day. Results: The results of this study showed that the yield of microorganisms (Y), the endogenous decay coefficient (k d ), the maximum specific growth rate (μ max ) and the saturation constant (K s ) were in the range of.67 g VSS/g,.56 d 1, 1.86 d 1 and 6.65 mg /l, respectively. Conclusions: The kinetic coefficients in this study can be used to improve the operation and design the MBR system in full scale. Please cite this article as: Ali Naghizadeh*, Elham Derakhshani, Ahmad Reza Yari, Study of Kinetic coefficients of a Membrane Bioreactor (MBR) for municipal wastewater treatment. Arch Hyg Sci 213;2(3):18-113. Background The membrane bioreactor (MBR) is a system that combines biological treatment with membrane filtration into a single process. The first application of MBR technology was reported in 1969, when an ultrafiltration membrane was used to separate activated sludge from the final effluent of a biological wastewater treatment system and the sludge was recycled back into the aeration tank (1). Submerged membrane bioreactors are increasingly used for domestic wastewater treatment, and there is a growing need to provide valuable modeling tools for their design and operation. In the years ago, the biological wastewater treatment processes design was based on the parameters such as hydraulic loading, organic loading, and detention time. Today, the design utilizes empirical and rational parameters based on biological kinetic equations. These equations describe growth of biological solids, substrate utilization rates, food-to-microorganisms ratio, and the mean cell residence time. Reactor volume, substrate utilization, biomass growth, and effluent quality can be calculated from those equations. Biokinetic coefficients of 18

Study of kinetic coefficient of a membrane... Naghizadeh A, et al./ Arch Hyg Sci 213;2(3):18-113 different wastewater treatment processes were evaluated by several investigators (2,3). However, kinetic coefficients used in the design of activated sludge processes include specific growth rate (μ), maximum rate of substrate utilization per unit mass of microorganisms (k), half-velocity constant, or substrate concentration at one-half the maximum specific growth rate (K s ), maximum cell yield (Y), and endogenous decay coefficient (k d ). Typical values of kinetic coefficients for activated sludge are shown in Table 1 (4). Table 1) Typical values of kinetic coefficients for activated sludge process Coefficient Basis Value Range Typical k (k=μ max /Y) day 1 2 8 4 k d day 1.3.7.5 K s mg/l, BOD 5 4 12 8 mg/l, 2 8 4 Y VSS/BOD 5.3.7.5 VSS/.2.5.4 The initial activated sludge system has been widely studied previously, leading to complete kinetic knowledge for modeling the main heterotrophic and autotrophic biological processes (5). Many of the studies are focused on the operational stability and usefulness of MBR for treating various wastewaters under different operating conditions (6, 7). Little has been found reported on the kinetic properties of the MBR process which are important in understanding the mechanism of the MBR process and in designing and controlling the system. Aims of the study: The purpose of this paper was to study the kinetic properties of the MBR system when it is applied to urban wastewater treatment. Materials & Methods Composition of Synthetic Wastewater (SWW) and operational conditions: A synthetic substrate was prepared based on the : N: P ratio 1: 5:1. Glucose, NH 4 Cl and KH 2 PO 4 were used as sources of carbon, nitrogen and phosphorus, respectively. The synthetic wastewater used in this study simulates municipal wastewater. It has been designed to provide all the inorganics and micronutrients, as well as nitrogen, phosphorous for the development of the biomass. The influent substrate concentration was fixed to mg/l for the biokinetic studies. The list of operational conditions is presented in Table 2. The MBR system was operated with a flow rate of 96 L/d. The hydraulic retention times (HRTs) of anoxic tank and oxic tank were 2.6 and 12 hours, respectively. Membrane Cleaning: Membrane fabric cleaning was achieved by combination of various methods including intermittent pumping of the suction pump, continuous air scouring by air from stone aerators in the reactor, back washing with water, back washing with air and mechanical cleaning of the membrane fabric with a brush. Technical data for membrane is presented in Table 3. Analytical Procedures: For the continuous reactor experiments, the sampling from the reactor and permeate were carried out periodically and analyzed for the following Physio-chemical parameters, by the methods described in the Standard Methods for the examination of wastewater (8). Model description: The purpose of studying the kinetic coefficients was to obtain information on the rate of cell growth and consumption of substrate. 19

Naghizadeh A, et al./ Arch Hyg Sci 213;2(3):18-113 Study of kinetic coefficient of a membrane... Table 2) Operating conditions of the submerged MBR Parameters UNIT Run 1 Run 2 Run 3 Run 4 Influent flow L/d 9 12 2 288 HRT h 12.8 9.6 5.8 4 SRT day 17 1 6 3 MLSS g/l 7.8 7.4 7.9 6.8 MLVSS g/l 5.8 5.6 5.4 5.3 Waste sludge L/d 4.8 2.4 1.2.6 OLR Kg /m3.d.5.5.5.5 Mixed liquor recycle % 2.5Q 2.5Q 2.5Q 2.5Q *HRT hydraulic retention time * OLR- Organic loading rate *The amount of MLSS and MLVSS are based on arithmetic mean in each run Table 3) Specifications of hollow fiber membrane Physical property Specifications Raw material polypropylene Inside diameter (μm) 32 Pore size (μm) Pore density (%) 4-5 This enabled the required volume of the reactor to be calculated and simulation of the system can be used for process control. The kinetic coefficients of a biological system have generally been determined experimentally using either continuous flow, completely mixed or batch lab-scale reactors. For soluble substrate, the substrate utilization rate in biological systems can be modeled with the following expression (9): The biomass growth rate is proportional to the substrate utilization rate by the yield coefficient, and the biomass decay is proportional to the amount of biomass present. When the substrate is being consumed at its maximum rate, the bacterial growth rate is also at its maximum (1). By substituting (k=µ max /Y) in Eq. (1), we will have: And also, Dividing both sides of Eq. (4) by the biomass concentration X, we obtain the specific growth rate as: The specific biomass growth rate (µ) can be defined as inverse of the solid retention time (SRT). Thus, Eq. (4) can be rearranged as follows: Results A summary of the results at the steady state conditions as well as kinetic coefficients obtained from the continuous system studies is present in Table 4 and 5 respectively. The plot for determination of non biodegradable soluble is shown in Figure 1. Fig 2 and 3 show the curves of the biokinetic determination. Taking into consideration influent and effluent substrate concentration: 11

1/U (day) U(day-1) Study of kinetic coefficient of a membrane... Naghizadeh A, et al./ Arch Hyg Sci 213;2(3):18-113 Run number Run 1 Run 2 Run 3 Run 4 Influent Effluent 5 8 1.4 14 Table4) Experimental parameters values in 4 runs Non biodegradable Biodegradable 3.8 7 9.5 13 1/S e (l/mg).25 4.7 1/SRT (day -1 ).59 7.34 Q ( S S ) e VX 6.22.38.55 Q VX ( S S e 6.31 4.57 2.64 1.82 ) Table 5) Summary of biokinetic coefficient for membrane bio reactor Type of study This study CF-MBR Municipal wastewater Industrial wastewater Y (mg/mg).67.276.4-.8.3-.72 k d ( day) -1.56.7.25-.48.45 µ m (day) -1 1.86.653 2-5.18.77 K s 6.65 62.3 11-372 5.8.5.4.3.2 y =.38x -.31 R² =.9633 5 1 15 2 1/SRT (day) -1.35.3.25.2 5.5 y =.674x -.557 R² =.9651 3.23.33.43.53.63 Fig. 1) determination of non biodegradable soluble 7 6 5 4 3 2 1 S Fig. 2) Determination of K and K s biokinetic coefficients y = 23.68x +.365 R² =.9227.5 5.2.25.3 1/Se (day) -1 Fig. 3) Determination of Y and K d biokinetic Discussion U (day) -1 coefficients In continuous flow completely mixed reactor, the determination of the kinetic coefficients is usually achieved by collecting data from lab-scale or pilot-plant experiments. Operating the system at various hydraulic retention times (HRT) and/or at various sludge retention times, and by allowing (at each adapted stage or HRT or SRT), a steady state condition to prevail. Accurate measurements of the biomass and permeate substrate concentration are then recorded. The parameters such as K s, μ, Y and k d can be determined through linearization of equation 5 as follows: 111

Naghizadeh A, et al./ Arch Hyg Sci 213;2(3):18-113 Study of kinetic coefficient of a membrane... To determine the kinetic coefficients, μm and Ks, equation 6 can be re arranged as: If equation 8 is plotted as Q(S -S)/VX versus 1/SRT, then from the slope and the intercept, it is possible to determine, the kinetic coefficients, k d and Y. Substituting the obtained value of k d in equation 9 and plotting SRT/[1+(SRT k d )] versus 1/S, then from the slope and the intercept it is possible to determine, the kinetic coefficients, μm and K s. A summary of the results at the steady state conditions as well as kinetic coefficients obtained from the continuous system studies is presented in table 4 and 5 respectively. It appears that the kinetic coefficients vary significantly with the change in MLSS concentration. However, this variability does not follow any definite pattern, and is not straight forward to draw a firm conclusion. This variability might be attributed to the nature of the system itself, since the system could be a selective process and the kinetic coefficients obtained might represent a different species. These finding are in accordance with Henze et al. (11). The plot for determination of non biodegradable soluble is shown in figure 1. As can be seen in this figure non biodegradable which is the interval between interception of the curve and the X axis is mg/l. this result is similar with the results of Rahman and Al-Malack (12). Fig 2 and 3 show the curves for determination of the biokientic determination. According to Fig 2 and from the slope and the intercept it is possible to determine, the kinetic coefficients, K and K s. table 2 show amount of these parameters. From the slope and the intercept of the curve presented in Fig 3, the biokinetic of Y and k d can be determined. Generally, the values of the kinetic coefficients which are presented in table 2 are (8) (9) within the normal range for the activated sludge process, but the values of K s, especially for MLSS of 74 mg/l are much higher than those reported in the literature. The estimation of the K s value is affected by an estimation of the decay rate, k d, thus any uncertainty in estimating k d will be reflected in K s. Y values were within the range of reported values treating glucose based synthetic wastewaters. However, Y values were lower than upper rage that for industrial plants treating the industrial wastewaters. This clearly shows that type of substrate has a role to play in the kinetic coefficients. These results are agreement with finding of Wisniewski et al. (13). Conclusion Kinetic analysis of the MBR system was successfully performed using different MLSS concentration. As a result of this study the biokinetic coefficients of MBR were similar to activated sludge process. According to this study Y, k d, μ max, and K s coefficients were found to be.67 g VSS/g,.56 d 1, 1.86 d 1 and 6.65 mg /l, respectively. The kinetic parameters determined in this study can be used to improve the design and operation of the MBR system on full scale. Footnotes Conflict of Interest: The authors declare no conflict of interest. References 1.Aileen N.L. Ng, A.S.K., A mini-review of modeling studies on membrane bioreactor (MBR) treatment for municipal wastewaters. Desalination 27; 212(1): 261-281. 2. Mardani, Sh., Mirbagheri, A., Amin, M. M., Ghasemian, M., determination of biokinetic coefficients for activated sludge processes on municipal wastewater, Iran. J. Environ. Health. Sci. Eng. 211; 8(1): 25-34. 3. Al-Malack, M., Determination of biokinetic coefficients of an immersed membrane bioreactor. Journal of Membrane Science 26; 271(1): 47 58. 112

Study of kinetic coefficient of a membrane... Naghizadeh A, et al./ Arch Hyg Sci 213;2(3):18-113 4. Qasim, S.R., Wastewater Treatment Plants: Planning, Design, and Operation. Second ed. 1999, Pennsylvania, USA: Technomic Publishing Company Inc, Lancaster. 5. M. Spérandio, M.C.E., Modelling an aerobic submerged membrane bioreactor with ASM models on a large range of sludge retention time. Desalination 28; 231(1): 82-9. 6. Naghizadeh, A., Mahvi, A.H., Vaezi, F., Naddafi, K., Evaluation of hollow fiber membrane bioreactor for municipal wastewater treatment. Iranian journal of environmentall health and engineering 28. 5(4): 257-268. 7. Naghizadeh, A., Mahvi, A.H, Mesdaghinia, A. R., Alimohammadi, M., Application of MBR Technology in Municipal Wastewater Treatment. Arab J Sci Eng 211; 36(1): 3-1. 8. APHA, WPCF, AWWA, (1999).Standard Methods for the Examination of Water and Wastewater, 2th ed., American Public Health Association (APHA), Washington, DC. 9. Haydar, S., Aziz, A., Kinetic Coefficients for the Biological Treatment of Tannery Wastewater Using Activated Sludge Process, Pak. J. Engg. & Appl. Sci 29; 5: 39-43. 1. Jin, R.C.; and Zheng, P. Kinetics of nitrogen removal in high rate anammox upflow filter, J. Hazard. Mater 29;17(2): 652 656. 11. Henze, M., Loosdrecht, M., Ekama, G., Brdjanovic, D., (28), Biological WastewaterTreatment, principles, modeling and design, IWA Publishing, New York. 12. Rahman, M. M., Al-Malack, M. H. Biochemical Kinetics of Cross flow Membrane Bioreactor Processes in the Treatment of Refinery Wastewater, Int. J. Environ. Res 212; 6(1):285-296. 13. Wisniewski, C., Leon Cruz, A., Grasmick, A., Kinetics of organic carbon removal by a mixed culture in a membrane bioreactor, Biochemical Engineering Journal 1999;3 (1): 61-69. 113