Numerical Modeling and Experimental Study of MFL Technique for Detection of Sub-surface Pittings in Pipelines a

Similar documents
NUMERICAL SIMULATION AND EXPERIMENTAL RESEARCH ON CRACK MAGNETIC FLUX LEAKAGE FIELD

STRESS CHARACTERIZATION BY LOCAL MAGNETIC MEASUREMENTS

Capabilities of MFL Inspection in DUPLEX Steel Pipelines

Fast Corrosion Screening Technique Saturation Low Frequency Eddy Current

CHARACTERIZATION OF MECHANICAL DAMAGE IN GAS TRANSMISSION. V.Kamat, C.Yeoh, P.lvanov, D.Kim, Y.Sun, S.Mandayam, L.Udpa, S.Udpa, and W.

Detection of hidden defects in thin steel plates using GMR sensor arrays

(MFL) 1 Mining University, St. Petersburg, Russia. Contact

Effect of Crack Orientation in Aluminium Plate using Tone Burst Eddy Current Thermography

Detecting metal loss and pipeline wall thinning with a combination of Ultrasonic Testing (UT) and Magnetic Flux Leakage (MFL)

Techniques for Detection of Surface Defects in Tubing and Pipe

Recent Trends in Electromagnetic NDE Techniques and Future Directions

Linear Magnetic Stray Flux Array based on GMR- Gradiometers

Feasibility Study of Measuring the Stress in Prestressing Tendons Using Induced Magnetic Field

Research Collection. Detection of breaks in reinforcing bars with the magnetic flux leakage method. Conference Paper. ETH Library

Advances In ILI Allow Assessing Unpiggable Pipelines. By Lisa Barkdull and Ian Smith, Quest Integrity Group, Houston

Department of Chemical & Polymer Engineering University of Engineering & Technology Lahore (FSD Campus)

APPLYING ULTRASOUND FOR IN-LINE INSPECTION: FACTS AND ISSUES. By Dr. Michael Beller, NDT Systems & Services AG, Stutensee, Germany

IN-LINE INSPECTION OF PIPES USING CORROSION RESISTANT ALLOYS (CRA) M.Sc. Johannes Keuter PPSA Seminar Aberdeen 19-Nov-2014

NEW APPROACHES FOR FLOW-ACCELERATED CORROSION

On the modelling of a Barkhausen sensor

SPE MS. Copyright 2015, Society of Petroleum Engineers

Comparison of MEC and MFL Technique for Tank Floor Inspection

Introduction and Overview of Condition Assessment. UKSTT Condition Assessment Masterclass

A comparison of DCVG measurements with in-line inspection results on a pipeline under AC influence

Integrative In-situ Detection of Defects in Aeroengine Blades Combining Borescope and Eddy Current Techniques

Finite Element Analysis of Workpiece Temperature during Surface Grinding of Inconel 718 Alloy

Detection of Sensitization for 600 Alloy and Austenitic Stainless Steel by Magnetic Field Sensor

Causes for residual magnetism on parts. Residual magnetism of steel products:

Fatigue Damage Evaluation of SUS304 Steel Using Magnetism Change in Fatigue Process

Abstract. 1. Introduction. 2. Literature Review

Magnetic Barkhausen Noise for hardness checking on steel

In-Line Inspection using Ultrasonic Technology

Non Destructive Testing with GMR Magnetic Sensor Arrays

Magnetic Sensor for Nondestructive Evaluation of Deteriorated Prestressing Strand Phase II. Final Report

Reading assignment. Nondestructive evaluation (NDE) Nondestructive testing (NDT) Penetrant. Conventional NDE mthods. Topic 7

Modeling of Residual Stresses in a Butt-welded Joint with Experimental Validation

LTSH boiler tube inspection in power plant by in-house developed Robotic inspection system Rakesh Kumar Chakraborty and Kishore Aggarwal

FATIGUE LIFE ASSESSMENT FOR COMPOSITE MATERIALS

Application Information Guidelines for Designing a Concentrator for High-Current Sensing Applications with an Allegro Hall-Effect Sensor IC

An Indian Journal FULL PAPER ABSTRACT KEYWORDS. Trade Science Inc. Visualization method research of butt weld magnetic flux leakage

EDDY CURRENT INSPECTION OF CLADDING MATERIAL ON VVER 440 REACTOR PRESSURE VESSELS

Feeler Pig A New Tool for Multi-Size Pipeline Inspection

Thickness Detection of Corroded Steel Plate by Low-Frequency Eddy Current Testing

ALLOY 600 REACTOR VESSEL HEAD PENETRATION ISSUES. Dudenstrasse 44, Mannheim, Germany. Post Office Box 355 Pittsburgh, PA ABSTRACT

Burst Pressure Prediction of Cylindrical Shell Intersection

Evaluating effect of magnetic flux leakage signals on fatigue crack growth of mild steel

When an axial load is applied to a bar, normal stresses are produced on a cross section perpendicular to the axis of the bar.

Bonded Neo Magnetization Guide

Finite Element Analysis of Concrete Filled Steel Tube Flexural Model

Numerical study of nugget formation in resistance spot welding

Buried Pipe NDE Overview

Dates of Innovation Development: January, 2014 to July 2015 Web site: Summary Description: Cathodic Protection

Corrosion related electromagnetic signatures measurements and modelling on a 1:40 th scaled model

Analysis and Experimental Study on Acoustic Emission Source Characteristics of Corrosion Process in Atmospheric Vertical Tank Bottom

EDDY CURRENT ARRAY TECHNOLOGY ON CRA LINER PIPES.

Determination of Metallic Fatigue in Nitrided Steel using a MDK Magnetic Nondestructive Tester

EWGAE 2010 Vienna, 8th to 10th September

Near field array for inspection of ferritic/ferromagnetic stainless steel (NFA FS)

Abstract. Nomenclature. A Porosity function for momentum equations L Latent heat of melting (J/Kg) c Specific heat (J/kg-K) s Liquid fraction

Finite Element Simulation of Residual Stresses in Butt Welding of Two AISI 304 Stainless Steel Plates

Crack detection in composite cantilever beam by Vibration analysis and Numerical method

SCIENCE & TECHNOLOGY

FLEXURAL BEHAVIOUR OF COLD FORM STEEL SECTIONS SUBJECTED TO STATIC LOADING

Development and Validation of an Automated Ultrasonic System for the Non- Destructive Evaluation of Welded Joints in Thermoplastic Storage Tanks

Other Major Component Inspection I. NDE for Preventing Radioactive Nuclides into Groundwater N. Muthu, EPRI, USA

Non-destructive Evaluation of Fatigue Damage in the Stainless Steel by Using Electromagnetic Method

Numerical and Experimental Modeling of Producer Gas Carburettor

Magnetic Flux Leakage Nondestructive Testing System and Services

Magnetic Sensor for Nondestructive Evaluation of Deteriorated Prestressing Strand. Phase I Final Report

Vibration-Induced Fatigue Failure of Discharge Bottle Fitting Connections A Case Study

Advancements in Imaging Corrosion Under Insulation (CUI) for Piping and Vessels

A Study on the Critical Thickness of the Inner Tube for Magnetic Pulse Welding Using FEM and BEM

TUBE INSPECTION SOLUTIONS TUBE INSPECTION SOLUTIONS. Selection of heat exchanger tube inspection techniques

Computation and analysis of temperature distribution in the crosssection

FUEL TUBE SPACER-PAD SPOT-WELD QUALITY ESTIMATION USING GUIDED ULTRASONIC WAVES

STATIC STRUCTURAL ANALYSIS OF 3 AXIS CNC MACHINE TABLE USING FINITE ELEMENT ANALYSIS

NEW TECHNOLOGIES OF SPETSNEFTEGAZ MAGNETIC IN-LINE INSPECTION

Glossary of Magnetic Terms

NEW SOUTH WALES DEPARTMENT OF EDUCATION AND TRAINING. Manufacturing and Engineering ESD. Sample Examination EA607 MAGNETIC PARTICLE TESTING

Development and validation of an automated ultrasonic system for the non-destructive evaluation of welded joints in thermoplastic storage tanks

INTRODUCTION TO ULTRASONIC IN-LINE INSPECTION OF CRA PIPELINES. Herbert Willems, Gerhard Kopp NDT Global, Germany

Evaluation of the Material Degradation of Austenitic Stainless Steel under Pulsating Tension Stress Using Magnetic Method

Preventing and Mitigating Corrosion to Advance Pipeline Integrity Pipeline Safety Trust Conference November 2, 2017 New Orleans, LA

Fitness-for-Service Evaluation of a Pipe with an Incomplete Penetration Flaw

Arch. Metall. Mater. 62 (2017), 3,

IMPACT ECHO TESTING OF IN-SITU PRECAST CONCRETE CYLINDER PIPE

Automated ultrasonic tube-to-tube sheet weld scanning

Ultrasonic Monitoring of Corrosion with Permanently Installed Sensors (PIMS)

Evaluation of eddy current and magnetic techniques for inspecting rebars in bridge barrier rails

Research Article Feasibility Study of Stress Measurement in Prestressing Tendons Using Villari Effect and Induced Magnetic Field

FINITE ELEMENT METHOD ANALYSIS OF STRESS INTENSITY FACTOR IN I CHANNEL SECTION

Weld defects analysis of 60 mm thick SS316L mock-ups of TIG and EB welds by Ultrasonic inspection for fusion reactor vacuum vessel applications

Overview of Non-Destructive Testing and Inspection Course Syllabus

Study on the Heat Transfer Characteristics Performed in the Infrared Thermography Detection of Welded Structure

Corrosion Imaging through Insulation using MR-MWM-Arrays

Statistical Analysis of TIG Arc Weldment Characteristics

NON-DESTRUCTIVE PRACTICAL EVALUATION OF PROD BASED MAGNETIC PARTICLE INDICATIONS OF HEAT AFFECTED ZONE IN AUSTENITIC GRADE RANGE MATERIAL.

Use of Piezoelectric Films for NDT Applications

ADVANCED POSSIBILITIES FOR CORROSION INSPECTION OF GAS PIPELINES USING EMAT TECHNOLOGY

University of Huddersfield Repository

Transcription:

Proceedings of the 5th Iranian International NDT Conference November 4-5, 2018, Olympic Hotel, Tehran Numerical Modeling and Experimental Study of MFL Technique for Detection of Sub-surface Pittings in Pipelines a Turaj Azizzadeh 1, Mir Saeed Safizadeh 2 More info about this article: http://www.ndt.net/?id=24158 1 PhD student, Iran university of science and technology; azizzadeh@mecheng.iust.ac.ir 2 Associate professor, Iran university of science and technology; safizadeh@iust.ac.ir Abstract Magnetic flux leakage (MFL) technique is currently the most commonly used method for corrosion monitoring in oil and gas pipelines. The present work employs 3D finite element modeling (FEM) to investigate the performance of MFL technique in detecting sub-surface pittings in pipelines. In order to simulate sub-surface pitting corrosion of oil and gas pipelines, rectangular pittings with different depths are modeled on a 10mm-thick steel plate. Then, the axial and radial components of MFL signals are obtained to investigate the influence of defect depth on detectability of sub-surface pittings by the MFL technique. Also, the influence of lift-off distance on MFL signals has been analyzed for enhancing the reliability of detection of defects. The results indicate that the designed magnetizer assembly renders detectability of pittings located at 9mm below surface in a 10mm-thick steel plate, which can meet the industry demand for in-service inspection. The MFL signals have also been obtained experimentally using an experimental set-up at the laboratory facilities on a steel plate containing sub-surface pittings. A good agreement between FEM and experimental results confirms that 3D FEM is an effective analysis method for MFL technique in the pipeline inspection. Keywords: Pipeline, Corrosion, Magnetic flux leakage. 1- Introduction Pipelines are a major means for transporting the gas and oil products over long distances. The pipelines exist in a corrosive environment. In such conditions, corrosion and erosion cause the thickness of the steel pipelines to decrease and fail and can even lead to leakage and rupture which cause fires, explosions and pollution of the environment. Therefore, periodic corrosion monitoring of the pipelines is one of the most important tasks for pipelines. In particular, the detection and characterization of sub-surface defects has always been of great industrial interest. Over the years, several non-destructive testing (NDT) methods have been investigated for corrosion monitoring of pipelines. Pipelines are often made of ferromagnetic materials with good permeability, so that the use of magnetic flux leakage (MFL) method is quite suitable for inspecting the corrosions of the pipelines and locating defects. MFL is one type of electromagnetic methods which is widely used for non-destructive inspection of steel pipelines. The basic principle of this method is that a powerful permanent magnet such as neodymium magnets is used to magnetize to saturation the steel under inspection. At areas where there is corrosion or metal loss, magnetic flux leaks from the steel and this flux is detected by a set of magnetic sensors. Then, the sensors signals will be processed to identify damaged areas and to estimate width, length and depth of metal loss [1-7]. Finite element method (FEM) is a convenient tool to model process of corrosion detection by means of MFL method. Two-dimensional (2D) simulation has been used to study the MFL signals from different defect sizes, materials, magnetizing levels and so on and it has also shown to be an effective technique. However, in 2D simulations, defects are considered as 2D geometry instead of really 3D geometry. Also, the 2D simulation falls short off estimating accurately the required magnetic saturation level in the test piece. Some of the research works on MFL technique using FEM have been conducted worldwide [8-11]. In this study, COMSOL-based 3D FEM model is used to simulate MFL signals from sub-surface pitting corrosions of different depths in the 10 mm-thick ferromagnetic plate. Also, the influence of lift-off distance on the MFL signals is studied. The designed magnetizer assembly renders detectability of pitting defects located at 9mm below surface in a 10mm-thick steel plate, which can meet the industry demand for in-service inspection. Then, the designed magnetizer assembly and test piece are manufactured. Finally, MFL tests are carried out and the numerical simulation results are validated with the experimental results 2- Finite element modeling (FEM) Maxwell s equations are used to model the MFL technique with the magneto static assumption. These equations are presented below [12-13]. H = 0 (1). B = 0 (2) Where H [A/m] is the intensity of the magnetic field, B [T] is the magnetic flux density. In addition, the following equations have to be used.

0H B = µ 0( H + M ) µ (3) Where µ 0 is the permeability of free space (4π 10-7 Tm/A) and M is the magnetization of the ferromagnet and permanent magnet. Vector potential A can be defined according to the below equation. Then, the system equations follow as: B = A (4) 1 ( A M ) = 0 µ 0 1 ( A) = 0 µ 0 The equation (5) is applied to ferromagnetic materials and permanent magnet. The equation (6) is applied to free space. Boundary conditions i.e. the value of A (vector potential) is considered zero at infinity and continuous at interfaces. The equations (5) and (6) are solved subject to defined boundary conditions to obtain the magnetic flux density detected by magnetic sensors. 3D FEM is a powerful tool for the analysis of the magnetizer assembly used in the MFL technique. It allows solving the non-linear equations governing the physical behavior of the system, and allows estimating the value of the magnetic field components at every point in the system. In this paper, COMSOL software was used to set up the 3D MFL model in the magneto static mode. 3D view of the MFL system together with the magnetizer arrangement is shown in Figure 1. The specimen length was 400mm, width was 140mm, and thickness was 10mm. The magnetizer assembly had two magnets, two couplings and one yoke. Yoke was meant to complete the magnetic circuit. Two neodymium magnets of thickness 40mm were used to magnetize the specimen to saturation. In this study, Magnetization of the permanent magnets was considered equal to 10e6 A/m. (5) (6) Figure 1. 3D model of MFL system The material of the specimen, couplings and yoke was steel X52. Figure 2 shows the B-H curve of X52 material. According to Figure 2, the strength of the magnetic field required to magnetically saturate the specimen is about 2Tesla. Figure 2. BH curve of X52 material Figure 3 schematically shows the specimen containing sub-surface pitting and lift-off distance. The rectangular subsurface pitting was assumed to be present at the center of the specimen. The width and length of the pitting were both

20mm. These dimensions were kept constant, while the depth of the pitting was varied from 10% (1mm) to 90% (9mm) of the thickness of the specimen. Figure 3. Schematic illustration of sub-surface pitting and lift-off distance Figure 4 shows the 3D model discretized by tetrahedral elements. The elements near defects were refined to obtain precise results. Figure 4. Finite element mesh After defining geometry, materials, boundary conditions and mesh, the solution by the finite element method can be calculated. Figure 5 shows the saturation levels inside the specimen obtained from 3D model and how the field lines travel through each element of the magnetizer assembly. As shown in Figure 5, the flux lines originate at the north pole of the right magnet, pass through a ferromagnetic coupling, enter the plate wall, pass through the other ferromagnetic coupling, and then enter the south pole of the left magnet. Inside the ferromagnetic yoke, the flux lines flow from the left to the right, thus completing the magnetic circuit. The appropriate magnetization level of the specimen is an important parameter of MFL technique for pipeline inspection. Figure 5. Magnetic flux behavior in the magnetizer assembly COMSOL allows visualizing of the amplitude of the magnetic flux components onto the scanning area as a C-scan image. Figure 6 shows the C-scan image of the axial and radial MFL components for the sub-surface pitting with a depth of 1mm. Similarly, Figure 7 shows the C-scan image of the axial and radial MFL components for the sub-surface pitting with a depth of 9mm. The solid squares in the C-scans represent the real position of the pitting. The lift-off was considered 2mm. The intensities of both components of the leakage fields for the pitting are found to increase with the increase in the pitting depth. Also, the length and width of the pitting can be easily predicted from the span of MFL signals (both radial and axial).

Figure 6. MFL C-scans for the sub-surface pitting with a depth of 1mm: (a) axial and (b) radial MFL component, respectively. Figure 7. MFL C-scans for the sub-surface pitting with a depth of 9mm: (a) axial and (b) radial MFL component, respectively. The axial and radial MFL B-scans for a range of pitting depths from 1mm to 9mm are shown in Figures 8(a) and (b), respectively. It can be seen that the intensity of MFL signals decreases with the decreasing depths of the sub-surface pitting. However, sub-surface pittings with different depths can be discriminated easily. Figure 8. MFL B-scans for a range of pitting depths from 1mm to 9mm: (a) axial and (b) radial MFL component, respectively. Pitting depth is the most important parameter for the pipeline integrity. As shown from the Figure 8, the designed magnetizer assembly renders detectability of pitting located at 9mm below surface in a 10mm-thick steel plate, which can meet the industry demand for in-service inspection. The peak-to-peak value of radial flux resulting from the sub-surface pitting corrosion of 10%t (1mm) depth is about 100gauss, which can easily be detected by magnetic sensors. Signal amplitudes for sub-surface pittings of different depths have been determined and are plotted in Figures 9(a) and (b) for both axial and radial MFL components. It can be seen that the relation between MFL amplitudes and pitting depths is somewhat linear.

Figure 9. The relation between MFL amplitudes and pitting depths: (a) axial and (b) radial MFL component, respectively. Figures 10(a) and (b) show the relation between MFL amplitudes and pitting depths at different lift-off distances for both axial and radial MFL components. It can be seen that the MFL amplitudes decrease with the increase in lift-off distance. However, the effect of lift-off distance on the radial MFL amplitudes is more pronounced than axial MFL amplitudes. Figure 10. The relation between MFL amplitudes and pitting depths at different lift-off distances: (a) axial and (b) radial MFL component, respectively. 3- Experimental tests Following the FEM, an experimental MFL system was established as shown in Figure 11(a). The measurement system consisted of a magnetizer assembly, a three-axis scanner, a data acquisition system and other associated units. The axial component of the MFL signals was measured with an mg-3002 gauss meter from Lutron. The distance between the specimen surface and the sensor was kept constant (2mm), as the sensor was moving by means of a scanner. The measurements were made on a ferromagnetic steel plate with a thickness of 10mm. The photo of the steel plate is shown in Figure 11(b). Rectangular pitting defects with different depths were machined at one side of the specimen. Figure 11(c) shows the details of the machined pittings. During the MFL measurements, the sensor was placed on the defect-free side of the specimen to detect the sub-surface pittings. The step resolution of the scanner was set to 4mm in axial direction and 5mm in circumferential direction. At every x-y position in the scanning area, the MFL signals were precisely measured and recorded on a computer for the further processing. Figure 11. (a) The experimental MFL system, (b) the photo of the steel plate, and (c) the geometrical details of the machined sub-surface pittings.

Figure 12 shows the result of the experimental test as a C-scan image. The positions of the pittings are given by the black squares. Figure 12 shows that magnetic flux spreads in the circumferential direction beyond the defects and overestimate the corrosion areas. This effect often called blooming effect. Figure 12. The experimental C-scan image Figure 13 shows the MFL B-scans of the experiment and FEM in the vicinity of the sub-surface pittings with depths of 4mm and 6mm. The FEM signal patterns are very similar to the experimental signals. The mean error is 9%. The reason for this relative error is mainly due to the small inaccuracies in the magnetic properties of the material used in the FEM with respect to the magnetic properties of the real material used in the experiment. Figure 13. Comparison between the experimental and FEM MFL signal in the vicinity of sub-surface pittings of 4mm and 6mm depth. 4- Conclusions This paper has investigated the performance of MFL technique for the detection of localized sub-surface pittings in ferromagnetic pipelines. The study was carried out through extensive finite element analysis in 3D using COMSOL software. In order to simulate sub-surface pitting corrosion of oil and gas pipelines, rectangular pittings with different depths were modeled on a 10mm-thick steel plate. Then, the axial and radial components of MFL signals were obtained to investigate the influence of defect depth on detectability of sub-surface pittings by the MFL technique. Also, the influence of lift-off distance on MFL signals was analyzed for enhancing the reliability of detection of defects. The results indicated that the designed magnetizer assembly renders detectability of pittings located at 9mm below surface in a 10mm-thick steel plate, which can meet the industry demand for in-service inspection. Finally, an experimental setup was established for evaluating the FEM results, which were in good agreement with experimental measurements. References 1. Singh, S., Rao, B., Vaidyanathan, S., Jayakumar, T., and Raj, B., Detection of Leakage Magnetic Flux from Near side and Far side Defects in Carbon Steel Plates using Giant-magneto resistive Sensor, 2008, Measurement Science and Technology, Vol. 19, pp. 15702-15710. 2. Gloria, N.B.S., Areiza, M.C.L., Miranda, I.V.J., and Rebello, J.M.A., Development of a magnetic sensor for detection and sizing of internal pipeline corrosion defects, 2009, NDT&E International, Vol. 42, pp. 669 677. 3. LiYing, S., YiBo, L., LiBo, S., and LingGe, L., 2012, Comparison of Magnetic Flux Leakage (MFL) and Acoustic Emission (AE) Techniques in corrosion Inspection for Pressure Pipelines, Proceedings of the 31 st Chinese Control Conference, Hefei, China. 4. Kim, D.K., Yoo, H.R., Yoo, J.S., Kim, D.K., Cho, S.H., Koo, S.J., Woo, R.Y., and Jung, H.K., 2013, Development of MFL System for In-Pipe Robot for Unpiggable Natural Gas Pipelines, 10 th International Conference on Ubiquitous Robots and Ambient Intelligence (URAI), Jeju, Korea. 5. Bubenic, T., Electromagnetic methods for detecting corrosion in underground pipelines: magnetic flux leakage (MFL), Woodhead Publishing Limited, Det Norske Veritas (U.S.A.) Inc., USA, 2014.

6. Mao, B., Lu, Y., Wu, P., Mao, B., and Li, P., Signal processing and defect analysis of pipeline inspection applying magnetic flux leakage methods, 2014, Intel. Serv. Robotics, Vol. 7, pp. 203 209. 7. Vanaei, H.R., Eslami, A., and Egbewande, A., A review on pipeline corrosion, in-line inspection (ILI), and corrosion growth rate models, 2017, International Journal of Pressure Vessels and Piping, Vol. 149, pp. 43-54. 8. Mitsuaki, K., Noritaka, M., and Kazumasa, N., Modeling of the yoke-magnetization in MFL-testing by finite elements, 2003, NDT&E International, Vol. 36, pp. 479-486. 9. Fengzhu, J., Changlong, W., Shiyu, S., and Weiguo, W., Application of 3-D FEM in the simulation analysis for MFL signals, 2009, Insight, Vol. 32, pp. 1784-1794. 10. Nemitz, O., and Schmitte, T., 2010, Simulation of Flaw Signals in a Magnetic Flux Leakage Inspection Procedure, Proceedings of the COMSOL Conference, Paris, France. 11. Singh, W., Thirunavukkarasu, S., Mahadevan, S., Rao, B., Mukhopadhyay, C., and Jayakumar, T., 2010, Three- Dimensional Finite Element Modeling of Magnetic Flux Leakage Technique for Detection of Defects in Carbon Steel Plates, Proceedings of the COMSOL Conference, India. 12. COMSOL Inc., 2016, URL: http://www.comsol.com. 13. Bray, E., Stanley, R.K., Nondestructive Evaluation - A tool in Design, Manufacturing, and Service, CRC Press Inc., 1997.