# Patent

## Title
On-line performance management of membrane separation process

## Abstract

Method for on-line prediction of performance of an RO based desalination plant is disclosed. The method includes: (i) a mathematical model of the RO unit; (ii) on-line estimation of membrane physical parameters of the nonlinear mathematical model representing the RO unit; and (iii) analysis of the estimated membrane transport parameter with respect to time. Based on the analysis of these estimated parameters, plant operators can clean the membranes to restore the performance of the RO desalination plant. The method can be implemented in a computer based control system used for data acquisition and control of an RO based desalination plant. The method can help in maintaining the performance of the RO based desalination plants at a desired level and increase membrane life without affecting quality of permeate water produced.

## Background
RELATED APPLICATION
This application claims priority as a continuation application under 35 U.S.C. §120 to PCT/IB2008/000359 filed as an International Application on Feb. 19, 2008 designating the U.S., the entire content of which is hereby incorporated by reference in its entirety


FIELD
The present disclosure relates to on-line performance monitoring, such as monitoring of reverse osmosis/nanofiltration plants by analyzing physical parameters of a membrane using a membrane transport phenomenological model.
BACKGROUND INFORMATION
Reverse Osmosis/Nanofiltration/Ultrafiltration (RO/NF/UF) is a pressure driven membrane separation process used in various industries such as desalination, wastewater treatment and chemical manufacturing. RO/NF/UF is used in plants to produce potable water from sea/brackish water. In an RO/NF/UF process, high pressure is applied on the feed side of the membrane to overcome the osmotic pressure of solute and cause transport of the solvent from the feed side to a permeate side and solute accumulates near the membrane surface. As a result, the concentration of the solute near the membrane surface increases gradually over a period, adversely affecting the performance of membrane. This phenomena is called concentration polarization. The concentration polarization is inversely proportional to the feed velocity across the membrane module. As recovery increases, the flow velocity across the membrane decreases, causing increased concentration polarization. Product recovery depends on other variables like feed concentration, pressure and temperature. In RO/NF/UF plants, the membrane fouling rate due to concentration polarization is influenced by multiple factors such as changes in feed concentration, temperature, and pressure, and it is difficult for the plant operator to determine the root cause for a changing fouling rate in a RO/NF/UF plant. Prediction of the changes in the fouling rate would help the plant operators in taking maintenance actions like cleaning the membrane to restore the performance to a desired level.
In industry, cleaning of the membrane is carried out in at least two ways; either based on a pressure drop between the feed and reject being more than a threshold value, or at predetermined fixed periodic intervals as per a recommendation by a membrane manufacturer. In the first method, the membrane may get damaged due to permanent fouling, and in the second method, membrane cleaning is independent of the actual fouling taking place in the membrane modules. Thus, both these methods of membrane cleaning are not satisfactory since the fouling rate changes with time and is dependent on the feed flow rate, concentration, pressure and temperature.
Different methods have been reported in literature for online cleaning and performance monitoring of a membrane separation process. Ooe Kenji and Okada Shingo [28] reported online method for performance analysis of an RO plant based on an ASTM D-4516 [1] method. The ASTM D-4516 method does not allow for discovering the development of membrane fouling or scaling until it results in significant loss of product quality such as product flow, and salt passage. In addition, this technique is applicable only where the plant is operated as per the design conditions and capacity with recovery being equal to or less than 15%.
Mohamad Amin Saad [16] extended the ASTM method to measure “Fouling Monitor” (FM) to monitor the performance of an RO plant. The FM is defined as a percentage difference between the normalized flux at design conditions and actual flux at the operating conditions of the RO plant. A cleaning scheduling of a membrane is arrived at based on the value of the FM. This method cannot predict the fouling of the membrane based on the operating conditions before normalized flux deviates from a design value. In addition, the method based on normalized flux may not be sufficient to predict fouling of a membrane accurately.
Nalco chemical company [18-27] has developed a method for monitoring the performance of a membrane separation process. As per the method, a tracer is injected in the feed stream and the concentration of tracer in outlet streams was estimated experimentally by using external sensors. The tracer concentrations in the feed and the outlet streams are used to monitor the fouling taking place in the membrane separation processes. This technique involves external sensors and tracer injection systems for implementation.
University Technology Corporation [U.S. Pat. No. 6,161,435] has developed a method and apparatus for monitoring membrane modules by using an ultrasonic sound technique. Due to fouling, the membrane thickness increases from the original value. Cleaning of the membrane is scheduled based on the monitoring of the membrane thickness using an ultrasonic technique. This method involves an individual ultrasonic transducer to monitor fouling at each membrane module.
The methods described above are not based on actual plant operating conditions and do not account for any time varying nature of fouling taking place in the membrane units.
Several mathematical models dealing with solute transfer in a membrane separation process have been reported in literature. Broadly, these membrane transport models may be divided in two categories (i) for neutral (reverse osmosis) membranes and (ii) for charged (nanofiltration and charged reverse osmosis) membranes. The mathematical models like preferential sorption-capillary flow model [2], Solution Diffusion model [3], Irreversible Thermodynamic model (Kedem-Katchalsky model [4] and Spiegler-Kedem model [5]), and Langmuir-type model [6] have been used for neutral membranes. In the case of charged membranes, the Nernst-Planck equation [7], electrostatic and steric hindrance model [8] have been used. Data driven models based on neural networks [9] have also been used to predict both permeate concentration and flux without solving any membrane transport equation.
Models proposed for charged membranes are developed by considering the chemical and physical properties of the solute and membrane such as solute size, solute charge, pore size of membrane and charge of membrane etc. On the other hand models based on the irreversible thermodynamics [4, 5] are developed by considering the membrane as a black box which has fluxes (permeate and solute flux) corresponding to the driving forces (pressure difference and concentration difference) of the transport process. The phenomenological constants are used to correlate flux and driving force, and physical parameters of the membrane are derived from these phenomenological constants. With irreversible thermodynamic models, the physical parameters of the membrane can be estimated for experimental data without knowing properties of membrane and solute. Soltanieh and Gill [10] compared the performance of the SK model and the KK model and observed that at no fouling condition, the membrane physical parameters of the KK model were found to be a function of feed concentration, while SK model parameters were found to be constant with respect to feed concentration. Several authors [11] compared the Solution Diffusion (SD) model with the SK model and concluded that the SK model predicts better than the SD model.
Murthy and Gupta. [12] proposed new a model, namely a Combined Film Spiegler-Kedem (CFSK) model, by including both membrane transport and concentration polarization effects. They concluded that CFSK model predictions are better than other models available in literature. Senthilmurugan et al [13] and Abhijit et al., [14] extended the CFSK model to spiral wound and hollow fiber modules respectively, and validated the models with experimental data with good results.
SUMMARY
A method is disclosed for real time performance management of membrane separation processes, comprising: predicting a state of fouling of a membrane based on an estimation of a time varying physical parameter of the membrane from plant data; and scheduling cleaning of the membrane based on a comparison of an estimated time varying physical parameter with a pre-defined threshold value.
A system is disclosed for real time performance management of a membrane separation process by performing a computer implemented program on a computer to implement a method comprising: predicting a state of fouling based on estimation of a time varying physical parameter of a membrane from plant data; and scheduling cleaning of the membrane based on a comparison of an estimated time varying physical parameter with a pre-defined threshold value.
A system is disclosed for real time estimation of a time varying physical parameter of a membrane separation process from plant data, comprising: means to measure a plant process variable in real-time; means to store a real-time measurement of plant operation data in a computer based control system; means to process the plant operation data stored in the computer based control system to remove noise; means to estimate a physical parameter of a membrane using a mathematical model; and means to store an estimated physical parameter in the computer based control system.



BRIEF DESCRIPTION OF THE DRAWINGS
Other features and advantages will become apparent to those skilled in the art upon reading the description of the preferred exemplary embodiments in conjunction with the accompanying drawings, wherein:
FIG. 1 is schematic of an exemplary RO/NF/UF plant with an associated instrumentation and control system;
FIG. 2 is an exemplary schematic of on-line performance monitoring system of an RO/NF/UF plant; and
FIG. 3 is an exemplary schematic of a mathematical model.



DETAILED DESCRIPTION
Exemplary embodiments can implement an on-line method that can analyze available plant data in terms of fouling of membranes and suggest an appropriate membrane cleaning schedule to plant operators to maintain the performance of the RO plant and also extend the life of the membrane. So far, such on-line performance monitoring methods based on a membrane transport phenomenological model have not been reported for RO/NF/UF plants and the present disclosure is aimed at filling such a gap.
In the present disclosure, an exemplary method for real-time estimation of a state of fouling and cleaning scheduling for RO/NF/UF plant is proposed. An exemplary method includes periodically executing the following steps: (i) using a phenomenological model to calculate the performance of an RO/NF/UF plant; (ii) on-line estimation of a membrane transport parameter of a phenomenological model at periodic intervals; and (iii) analysis of the membrane transport parameter to determine the state of the fouling of the membrane.
The following exemplary physical parameters of the membrane are estimated online:

    
    
        Hydrodynamic permeability of membrane (A)
        Solute permeability i.e. Permeability of solute with respect to membrane (Pm)
        Reflection coefficient of membrane (σ)
    
    


The present disclosure provides for performance monitoring of a membrane unit through online analysis of physical parameters of the phenomenological model of the membrane transport process. This method provides information about the time varying rate of fouling in a membrane unit, which can be used in scheduling of the membrane cleaning.
A proposed on-line performance monitoring method includes:

    
    
        A mathematical model for both hollow fiber and spiral wound membrane modules with listed inputs (model parameters, feed flow rate, conductivity and pressure) and outputs (permeate flow rate, permeate conductivity, and reject pressure) as shown in FIG. 3.
        Online estimation of both physical parameters of the membrane and parameters related to the configuration of a membrane module by minimizing the error between measured and predicted values of permeate conductivity, flow rate, and reject pressure. A non-linear optimization technique can be used for minimizing the error between the predicted and measured values.
        Online validation of a developed RO/NF/UF unit model with plant data by comparing the model predicted values with actual operation data from the plant.
        Online estimation of physical parameters of the membrane in a regular time interval by minimizing the error between the measured and predicted values of permeate conductivity, flow rate, and reject pressure. A non-linear optimization technique can be used for minimizing the error between the predicted and measured values.
        Analysis of the estimated membrane transport parameter by comparing the current estimated parameter values with predefined threshold values. If the values of the current estimates of the parameter values is more than the threshold values, recommend cleaning of the membranes.
    
    


FIG. 1 illustrates a schematic of an RO/NF/UF plant with associated instrumentation and plant control system
The RO/NF/UF based desalination plant has following streams, namely feed, reject and permeate streams. The feed is pretreated 2 before being pumped to RO/NF/UF membrane module through high pressure pump 3. The properties of the feed stream such as conductivity, pressure and flow rate are measured by corresponding sensors 4, 5, 6. The RO Modules network 7 is connected to the sensors 6 and 8, the RO/NF/UF membrane module purifies the feed water and purified water is collected at a permeate end and concentrated water is collected at a reject end. The process variables such as reject flow rate and pressure are measured at corresponding sensors 8, 9. Similarly, other process variables such as permeate flow rate and conductivity are measured by corresponding sensors 10, 11. This measured data from sensors are stored in plant control system 1. These measurements are carried under two conditions such as (i) normal operating conditions and (2) introducing at least one disturbance such as a step change in any one process variable.
FIG. 2 illustrates an exemplary schematic of the online performance monitoring system of an RO/NF/UF plant using mathematical model 22 and analysis of the estimated model parameters, which change with time depending upon the plant operating conditions.
Various exemplary steps involved in online parameter estimation method are:

    
    
        The data stored in Distributed Control System 1 is processed 16 to remove the noise.
        Estimation of model parameters or physical parameters of membrane 17
        Validation of model parameter 18.
        Analysis of model parameter 19.
        Recommendation goes to the operator panel 21 for cleaning, on expiry of the estimated time 21.
    
    


The parameter estimation can be carried out by minimizing the error between the predicted and measured process variable under normal operating conditions. The error minimization can be performed by a non-linear optimization technique. Further estimated model parameters are used to validate the model using measured process variables.
FIG. 3 illustrates an exemplary mathematical model of an RO/NF/UF plant. The mathematical model for the membrane module will change depending upon the configuration of the module used in a plant, namely a Hollow Fiber (HF) module or a Spiral Wound (SW) module or a tubular module 23. A brief description of the mathematical models of both HF and SW modules is outlined below
Model for HF Module
The permeate flow rate and solute concentration obtained from a given HF module can be predicted [14] by solving a set of equations which describe the mass transfer processes in the module. These equations namely, the membrane transport model, concentration polarization model, local solvent and solute mass balances are all applicable at any point within the permeator. The system of coupled differential equations may be solved numerically using the finite difference method.
The following assumptions have been made in the development of our analysis:

    
    
        The bulk stream flows radially outward and there is sufficient axial mixing in the bulk stream. This implies that the bulk flow variables are only dependent on r and it allows for replacing the partial derivative terms that appear in the material balance equations and the pressure drop equation with ordinary derivatives.
        The element chosen for finite difference analysis within the permeator is much larger than the fiber dimensions. Hence, for all practical purposes the shell side of the membrane can be assumed to be a continuous phase.
        Membrane structure is uniform throughout the module. All model parameters within the permeator are constant.
        There is no variation in bulk flow properties of the feed stream.
        Solution contains only one salt and a solvent (binary solution).
        Film theory is applicable within the membrane module.
        Fluid properties and diffusivities remain constant inside the module.
    
    


By combining the membrane transport equation of Spiegler-Kedem [5] model and film theory based concentration polarization model [12] equation, we obtain:




  
    
      
        
          Permeate
          ⁢
          
            
          
          ⁢
          
            flux
            ⁡
            
              (
              
                
                  m
                  3
                
                ⁢
                
                  /
                
                ⁢
                
                  
                    m
                    2
                  
                  ·
                  s
                
              
              )
            
          
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            J
            v
          
        
        =
        
          
            A
            ρ
          
          ⁡
          
            [
            
              
                (
                
                  
                    P
                    b
                  
                  -
                  
                    P
                    p
                  
                
                )
              
              -
              
                σ
                ⁢
                
                  
                    
                      vR
                      G
                    
                    ⁢
                    T
                  
                  
                    M
                    w
                  
                
                ⁢
                φ
                ⁢
                
                  
                
                ⁢
                
                  
                    C
                    b
                  
                  ⁡
                  
                    (
                    
                      
                        
                          1
                          -
                          F
                        
                        φ
                      
                      
                        
                          
                            1
                            -
                            F
                          
                          φ
                        
                        +
                        
                          
                            1
                            -
                            σ
                          
                          σ
                        
                      
                    
                    )
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        1
        )
      
    
  




Where A is membrane hydrodynamic permeability (m3/m2·s·Pa), σ is reflection coefficient of membrane (−), ρ is density of sea water, Pb and Pp are pressures of feed side bulk stream and permeate stream at membrane local point (Pa), ν is vont-hoff factor of solute (−), RG is gas constant (J·kmol−1·° K−1), T is temperature (° K), Mw is molecular weight of solute (kg/kmol), φ is concentration polarization defined by equation (3), Cb is concentration of bulk feed at membrane local point (kg/m3), F is intermediate dummy variable which defined by equation (3),





  
    
      
        
          Permeate
          ⁢
          
            
          
          ⁢
          concentration
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            C
            p
          
        
        =
        
          
            C
            b
          
          
            1
            +
            
              
                σ
                
                  1
                  -
                  σ
                
              
              ·
              
                
                  1
                  -
                  F
                
                φ
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  
  
    
      
        Where
        ,
        
          φ
          =
          
            exp
            ⁡
            
              (
              
                
                  J
                  v
                
                k
              
              )
            
          
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            F
          
          =
          
            exp
            ⁡
            
              (
              
                
                  -
                  
                    J
                    v
                  
                
                ⁢
                
                  
                    1
                    -
                    σ
                  
                  
                    P
                    m
                  
                
              
              )
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




Where, Pm is solute permeability (m/s).

The mass transfer coefficient (k) used in equation (3) can be expressed as a function of the Reynolds and Schmidt numbers.

Sh=aRebSc1/3  (4)

Equations of the same form are used in literature for estimating the mass transfer coefficients. The values of ‘a’ and ‘b’ for a hollow fiber module have been reviewed by Masaaki Sekino [29] for an HFRO module.
The pressure difference across the membrane which is used in equation 1 for obtaining the permeate flux varies throughout the membrane because of friction losses. The pressure drop for the permeate and bulk streams can be estimated using a Hagen-Poiseuille equation and the modified Ergun's [15] equation respectively. These equations are given below




  
    
      
        
          Hagen
          ⁢
          
            -
          
          ⁢
          Poiseuille
          ⁢
          
            
          
          ⁢
          equation
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            d
            
              d
              z
            
          
          ⁢
          
            P
            p
          
        
        =
        
          
            -
            
              
                32
                ⁢
                μ
              
              
                d
                i
                2
              
            
          
          ⁢
          
            v
            p
          
        
      
    
    
      
        (
        5
        )
      
    
  




Where vp is permeate velocity (m/s) at inside the fiber bore, di is an inside diameter of hollow fiber (m), μ is viscosity of water (Pa·s), z is axial coordinate

The modified Ergun [15] equation for pressure drop per length of the packed bed at a turbulent condition can be written as




  
    
      
        
          
            ∂
            
              P
              b
            
          
          
            ∂
            r
          
        
        =
        
          c
          ·
          
            v
            r
            d
          
          ·
          
            J
            v
            e
          
        
      
    
    
      
        (
        6
        )
      
    
  




Where, vr is superficial velocity of feed stream (m/s), c, d, e constants are used in equation (6).

The material balance equations for both solute and solvent streams within the module are given below




  
    
      
        
          
            Permeate
            ⁢
            
              
            
            ⁢
            stream
            ⁢
            
              :
            
            ⁢
            
              
            
            ⁢
            
              d
              
                d
                z
              
            
            ⁢
            
              v
              p
            
          
          =
          
            
              J
              v
            
            ⁢
            
              ζ
              /
              θ
            
            ⁢
            
              
            
            ⁢
            such
            ⁢
            
              
            
            ⁢
            that
            ⁢
            
              
            
            ⁢
            BC
          
        
        ,
        
          

        
        ⁢
        
          
            
              v
              p
            
            ⁢
            
              ❘
              
                z
                =
                0
              
            
          
          =
          0
        
        ,
        
          
            0
            ≤
            z
            ≤
            
              L
              ⁢
              
                

              
              ⁢
              Here
              ⁢
              
                
              
              ⁢
              θ
            
          
          =
          
            
              
                d
                i
                2
              
              ⁢
              N
            
            
              
                D
                o
                2
              
              -
              
                D
                i
                2
              
            
          
        
        ,
        
          ζ
          =
          
            
              
                4
                ⁢
                θ
                ⁢
                
                  
                
                ⁢
                
                  d
                  o
                
              
              
                d
                i
                2
              
            
            ·
            
              
                L
                
                  L
                  m
                
              
              .
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




The length of a hollow fiber is given as, L=√{square root over (Lm2+4(πrW)2)}
Lm is length of module (m).
Bulk stream solute concentration:




  
    
      
        
          
            d
            
              d
              r
            
          
          ⁢
          
            (
            
              rv
              r
            
            )
          
        
        =
        
          
            
              -
              θ
            
            ⁢
            
              
                rv
                p
              
              L
            
          
          ⁢
          
            |
            
              z
              =
              L
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




subject to BC, νr|r=Di/2=νF 

vF is velocity of feed at feed header (m/s)

Likewise for the solute,





  
    
      
        
          
            d
            
              d
              r
            
          
          ⁢
          
            (
            
              
                rv
                r
              
              ⁢
              
                C
                b
              
            
            )
          
        
        =
        
          
            
              -
              θ
            
            ⁢
            
              
                
                  rv
                  p
                
                ⁢
                
                  C
                  p
                
              
              L
            
          
          ⁢
          
            |
            
              z
              =
              L
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




subject to BC, Cb|r=Di/2=CF for Di/2≦r≦Do/2

CF is feed concentration (kg/m3).

Differentiation of equation (3.8) and subsequent substitution into equation (3.11) leads to:





  
    
      
        
          
            
              d
              2
            
            
              d
              
                z
                2
              
            
          
          ⁢
          
            P
            p
          
        
        =
        
          
            
              -
              
                
                  32
                  ⁢
                  μ
                
                
                  d
                  i
                  2
                
              
            
            ·
            
              ζ
              θ
            
            ·
            
              J
              v
            
          
          ⁢
          
            
          
          ⁢
          BC
          ⁢
          
            
          
          ⁢
          
            {
            
              
                
                  
                    
                      
                        
                          d
                          
                            d
                            z
                          
                        
                        ⁢
                        
                          P
                          p
                        
                      
                      ⁢
                      
                        |
                        
                          z
                          =
                          0
                        
                      
                    
                    =
                    0
                  
                
              
              
                
                  
                    
                      P
                      p
                    
                    =
                    
                      
                        
                          P
                          atm
                        
                        -
                        
                          
                            
                              l
                              s
                            
                            ·
                            
                              
                                32
                                ⁢
                                μ
                              
                              
                                d
                                i
                                2
                              
                            
                          
                          ⁢
                          
                            v
                            p
                          
                        
                      
                      ⁢
                      
                        |
                        
                          z
                          =
                          L
                        
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




Where, Is is length of epoxy seal (m), Patm is atmosphere pressure (Pa).

The above equations (1) to (9) are solved numerically by the finite difference method with each of the variables being expressed as a discrete value. Since the permeate flow variables vary only along the z-axis while the bulk flow terms vary along the r-axis, the equations are solved sequentially by proceeding from r=Di/2 to Do/2 while solving all the z-axis dependent difference equations at a particular radial grid location. The bulk flow terms at r=Di/2 are known; Pb=PF, Cb=CF and vr=VF.
Model for SW Module
The following assumptions have been made in the development of our analysis.


    
    
        Membrane structure is uniform throughout the module. All model parameters within the permeator are constant.
        There is no variation in bulk flow properties.
        Solution contains only one salt and a solvent (binary solution).
        Film theory is applicable within the membrane module.
        Fluid properties and diffusivities remain constant inside the module.
    
    


The mass transport equations of membrane will be same for both HF and SW modules. Therefore, the equations 1-4 are solved with the following pressure drop and mass balance equation given below for SW module [13].
The pressure drop in both the channels is based on the assumption that Darcy's law is applicable. This leads to the following expression for the pressure drops:




  
    
      
        
          Feed
          ⁢
          
            
          
          ⁢
          Channel
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            
              d
              
                P
                b
              
            
            
              d
              x
            
          
        
        =
        
          
            k
            fb
          
          ·
          μ
          ·
          
            U
            b
          
        
      
    
    
      
        (
        10
        )
      
    
  
  
    
      
        
          Permeate
          ⁢
          
            
          
          ⁢
          Channel
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            
              d
              
                P
                p
              
            
            
              d
              y
            
          
        
        =
        
          
            k
            fp
          
          ·
          μ
          ·
          
            U
            p
          
        
      
    
    
      
        (
        11
        )
      
    
  




Where: kfb is the friction parameter in the feed channel (1/m2), kfp is the friction parameter in the permeate channel (1/m2), Ub, Up is the velocity of the solution in feed and permeate channels (m/s) and μ is the viscosity of solution (Pa·s). Here both friction parameters are experimentally determined for a given module, and x and y are directions of feed and permeate flow when in a module unwind condition.

The overall material balance for the feed and the permeate sides are given by the following equations:




  
    
      
        
          
            d
            
              U
              b
            
          
          
            d
            x
          
        
        =
        
          
            -
            2
          
          ⁢
          
            
              J
              v
            
            /
            
              h
              b
            
          
        
      
    
    
      
        (
        12
        )
      
    
  
  
    
      
        
          
            d
            
              U
              p
            
          
          
            d
            y
          
        
        =
        
          2
          ⁢
          
            
              J
              v
            
            /
            
              h
              p
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




Where, hb, hp are thickness of feed and permeate side spacer (m).

Similarly the material balance for the solute on the feed side is represented by the following equation:





  
    
      
        
          
            
              d
              
                U
                b
              
            
            ⁢
            
              C
              b
            
          
          
            d
            x
          
        
        =
        
          
            -
            2
          
          ⁢
          
            J
            v
          
          ⁢
          
            
              C
              p
            
            /
            
              h
              b
            
          
        
      
    
    
      
        (
        14
        )
      
    
  




Differentiating equation (10) with respect of “x” and substituting in equation (12), we obtain:





  
    
      
        
          
            
              d
              2
            
            ⁢
            
              P
              b
            
          
          
            d
            
              x
              2
            
          
        
        =
        
          2
          ⁢
          
            k
            fb
          
          ⁢
          μ
          ⁢
          
            
          
          ⁢
          
            
              J
              v
            
            /
            
              h
              b
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




with boundary conditions




  
    
      P
      b
    
    =
    
      
        
          P
          F
        
        ⁢
        
          
            
          
          ⁢
          
            
          
        
        ⁢
        at
        ⁢
        
          
        
        ⁢
        x
      
      =
      0
    
  




  and




  
    
      P
      R
    
    =
    
      
        
          P
          F
        
        +
        
          
            ∫
            0
            L
          
          ⁢
          
            
              (
              
                
                  d
                  
                    P
                    b
                  
                
                
                  d
                  x
                
              
              )
            
            ⁢
            
              d
              x
            
            ⁢
            
              
            
            ⁢
            at
            ⁢
            
              
            
            ⁢
            x
          
        
      
      =
      L
    
  



Similarly differentiating equation (11) with respect to y and substituting in equation (13) we obtain:





  
    
      
        
          
            
              d
              2
            
            ⁢
            
              P
              p
            
          
          
            d
            
              y
              2
            
          
        
        =
        
          
            -
            2
          
          ⁢
          
            k
            fp
          
          ⁢
          μ
          ⁢
          
            
          
          ⁢
          
            
              J
              v
            
            /
            
              h
              p
            
          
        
      
    
    
      
        (
        16
        )
      
    
  




with boundary conditions:




  
    
      P
      p
    
    =
    
      
        
          P
          atm
        
        ⁢
        
          
        
        ⁢
        at
        ⁢
        
          
        
        ⁢
        y
      
      =
      w
    
  




  
    
      
        and
        ⁢
        
          
        
        ⁢
        at
        ⁢
        
          
        
        ⁢
        y
      
      =
      0
    
    ,
    
      

    
    ⁢
    
      
        P
        pw
      
      =
      
        
          P
          atm
        
        -
        
          
            ∫
            0
            w
          
          ⁢
          
            
              (
              
                
                  d
                  
                    P
                    b
                  
                
                
                  d
                  y
                
              
              )
            
            ⁢
            
              d
              y
            
          
        
      
    
  



PR is reject pressure (Pa), L is length of spiral wound module (m), and w=width of module with respect to number of wounds (m).

The above equations are solved using the method of finite differences. The feed flow path (x direction) is divided into m segments while the permeate flow path (y direction) is divided into n segments.
By solving the above model equations of HF and SW modules, the permeate flux, and concentration at local points of the membrane module can be estimated. The overall permeate concentration and flow rate can be estimated by the following equations:




  
    
      
        
          Q
          p
        
        =
        
          
            ∫
            
              x
              =
              0
            
            
              x
              =
              m
            
          
          ⁢
          
            
              ∫
              
                y
                =
                0
              
              
                x
                =
                n
              
            
            ⁢
            
              
                J
                v
              
              ⁢
              
                S
                m
              
              ⁢
              
                d
                y
              
              ⁢
              
                d
                x
              
            
          
        
      
    
    
      
        (
        17
        )
      
    
  
  
    
      
        
          C
          pt
        
        =
        
          
            ∫
            
              x
              =
              0
            
            
              x
              =
              m
            
          
          ⁢
          
            
              ∫
              
                y
                =
                0
              
              
                x
                =
                n
              
            
            ⁢
            
              
                J
                v
              
              ⁢
              
                S
                m
              
              ⁢
              
                C
                p
              
              ⁢
              
                d
                y
              
              ⁢
              
                d
                x
              
            
          
        
      
    
    
      
        (
        18
        )
      
    
  




Where Sm is surface area of a membrane corresponding finite element.

The exemplary lists of physical parameters 24 used in the model are:

    
    
        Membrane hydrodynamic permeability (A)
        Reflection coefficient of membrane (σ)
        Solute permeability (Pm)
        The constants of mass transfer coefficients correlation (a and b)
        The constants of modified Ergun's equation for HFRO module (c, d, e) or Darcy's law constant for feed and permeate channel for spiral wound module (kfb, kfp).
    
    


The conductivity of permeate can be estimated from the permeate concentration.
The above described mathematical models are used in the present method to describe the physical phenomena occurring in membrane separation processes. The models include parameters such as solute permeability, hydrodynamic permeability and membrane reflection coefficient to characterize the fouling phenomena. These model parameters are time varying in nature and are estimated periodically from the RO plant data such as flow rate, temperature, pressure and quality of feed, reject and permeate. Analysis of these estimated parameters will indicate the rate of fouling taking place in the RO plant and the cleaning of the membrane is recommended whenever the values of these parameters exceed a pre-defined threshold value.
It will be appreciated by those skilled in the art that the present invention can be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The presently disclosed embodiments are therefore considered in all respects to be illustrative and not restricted. The scope of the invention is indicated by the appended claims rather than the foregoing description and all changes that come within the meaning and range and equivalence thereof are intended to be embraced therein.
References, all of which are incorporated herein by reference in their entireties and referenced by number in the specification.


    1. ASTM D4516-00 (2006)e1 Standard Practice for Standardizing Reverse Osmosis Performance Data
    2. R. Rangarajan, T. matsuura, E. C. Goodhe, and S. Sourrirjan, Predictability of reverse osmosis performance of porous cellulose acetate membranes for mixed uni-valent electrolytes in aqueous solutions, Ind. Eng. Chem. Prod. Des. Dev, 17 (1978) 46-56
    3. J. G. Wijmans and R. W. Baker, The solution-diffusion model: a review, Journal of Membrane Science 107 (1995) 1-21
    4. O. Kedem and A. Katchalsky, Thermodynamic analysis of the permeability of biological membranes to non electrolytes, Biochim. Bio-Phys. Acta. 27 (1958) 229.
    5. K. S. Spiegler and O. Kedem, Thermodynamics of hyperfiltration (reverse osmosis): criteria for efficient membranes, Desalination 1 (1966) 311-326.
    6. M. Soltanieh, and S. Sahebdelfar, Interaction effects in multi-component separation by reverse osmosis, J. Membr. Sci, 183 (2001) 15-27
    7. M. W. Vonk and J. A. M. Smit, Positive and negative ion retention curves of mixed electrolytes in reverse osmosis with a cellulose acetate membrane. An analysis on the basis of the generalized Nernst-Planck equation, J. of Colloid and Interface Sci., 96 (1983) 121-134
    8. X. Wang, T. Tsuru, M. Togoh S.-I. Nakao and S. Kimura, Transport of organic electrolytes with electrostatic and steric-hindrance effects through Nan filtration membranes, J. Chem. Engg. Japan, 28 (1995) (372-380)
    9. Grishma R. Shetty a, Shankararaman Chellam, Predicting membrane fouling during municipal drinking water nanofiltration using artificial neural networks, J. Membr. Sci., 217 (2003) 69-86
    10. M. Soltanieh and W. N. Gill, Review of reverse osmosis membranes and transport models, Chemical Engg. Comm., 12 (1981) 279
    11. A. Mason, H. K. Lonsdale, Statistical mechanical theory of membrane transport, J. Membr. Sci. 51 (1990) 1.
    12 Z. V. P. Murthy and S. K. Gupta, Thin film composite polyamide membrane parameters estimation for phenol-water system by reverse osmosis, Sep. Sci. Technol., 33(16)(1998) 2541E.
    13 S. Senthilmurugan, Aruj Ahluwalia and Sharad K. Gupta, Modeling of a spiral wound reverse osmosis module and estimation of model parameters using numerical techniques”, Desalination, 173, 269-286, 2005
    14 Abhijit Chatterjee, Aruj Ahluwalia, S. Senthilmurugan and Sharad K. Gupta, Modeling of a Radial flow hollow fiber module and estimation of model parameters using numerical techniques”, Journal of Membrane Science, 236, 1-16, 2004
    15. Senthilmurugan s and Babji B S, Hydrodynamics studies in radial flow hollow fiber reverse osmosis module, International Conference on Modeling and Simulation, Coimbatore, 27-29 Aug. 2007
    16. Mohamad Amin Saad, Early discovery of RO membrane fouling and real-time monitoring of plant performance for optimizing cost of water, Desalination 165 (2004) 183-191
    17. U.S. Pat. No. 6,161,435 Method and apparatus for determining the state of fouling cleaning of membrane modules
    18. U.S. Pat. No. 6,699,684 Method of monitoring Biofouling membrane separation processes
    19. U.S. Pat. No. 6,730,227 Method of monitoring membrane separation processes
    20. U.S. Pat. No. 6,821,428 Method of monitoring membrane separation processes
    21. U.S. Pat. No. 6,838,001 Method of monitoring membrane separation processes
    22. U.S. Pat. No. 6,838,002 Method of monitoring membrane separation processes
    23. U.S. Pat. No. 7,060,136 Method of monitoring membrane cleaning processes
    24. U.S. Pat. No. 7,169,236 Method of monitoring membrane cleaning processes
    25. U.S. Pat. No. 6,475,394 Pseudo-fouling detector and use thereof to control an industrial water process
    26. U.S. Pat. No. 6,017,459 Apparatus and method for the monitoring of membrane deposition
    27. U.S. Pat. No. 7,252,096 Methods of simultaneously cleaning and disinfecting industrial water systems
    28. Ooe Kenji and Okada Shingo, “eCUBE aqua” application portfolio for reverses osmosis membrane diagnosis. Yaokogawa Technical Report English Edition No 38 (2004).
    29. Masaaki Sekino, Mass Transfer Characteristics of Hollow Fiber RO Modules, Journal of Chemical Engineering of Japan, 28 (1995) 843-846

## Claims

1. A method for real time performance management of a membrane in a plant controlled by a computer based control system, wherein the performance management of the membrane includes monitoring of fouling of the membrane taking place during a membrane separation process, and cleaning of the membrane for restoration of performance of the membrane, the method comprising:
measuring a plant process variable in real time during the fouling of the membrane; and
processing the measured plant process variable in the computer based control system by:
storing the real time measurement of the plant process variable;
estimating, at periodic instants in time during the fouling of the membrane, time varying physical parameters of the membrane from the stored real time measurement of plant process variable using a mathematical model, the physical parameters including at least a hydrodynamic permeability of the membrane and a solute permeability with respect to the membrane;
analyzing the estimated time varying physical parameters to determine a time varying rate of fouling;
predicting a state of fouling of the membrane based on the determined time varying rate of fouling and estimating a scheduled time for the cleaning of the membrane based on a comparison of the estimated time varying physical parameters with predefined threshold values at a respective periodic instant of time during the fouling of the membrane; and
scheduling the cleaning of the membrane based on the estimated schedule time for the cleaning of the membrane by the computer based control system;
sending a signal to an operator panel for recommending the cleaning of the membrane based on the estimated schedule time for cleaning; and

cleaning of the membrane at the estimated scheduled time based on the recommendation sent to the operator panel.

2. The method of claim 1, comprising:
on-line or off-line predicting and scheduling.

3. The method of claim 1, wherein the membrane used in a membrane separation process is one of a reverse osmosis, nanofiltration or ultrafiltration.

4. The method of claim 1 wherein real time estimation of the time varying physical parameters of a membrane separation process from the plant data comprises:
processing the data stored in the computer based control system to remove noise; and
storing the estimated parameters in the computer based control system.

5. The method of claim 1 wherein the physical parameters of the membrane include the hydrodynamic permeability of the membrane, the solute permeability with respect to the membrane, and a reflection coefficient of the membrane.

6. The method of claim 4, wherein the plant process variable includes at least one of flow rate, pressure, temperature, and conductivity of feed, flow rate and conductivity of permeate, and pressure of reject stream.

7. The method of claim 4, wherein the computer based control system includes hardware and software of a distributed control system, a programmable logic controller (PLC) or any microprocessor based embedded system.

8. The method of claim 4 wherein the mathematical model represents one of a hollow fiber, spiral wound or tubular type membrane module.

9. The method of claim 1, wherein on-line measurement of the plant process variables is performed under normal operating conditions.

10. The method of claim 1, comprising:
performing on-line measurement of a plant process variable by introducing at least one disturbance in any one process variable.

11. The method of claim 1 for scheduling cleaning of the membrane, comprising:
formulating an empirical model to characterize changes in the physical parameters of the membrane;
analyzing the estimated physical parameters of the membrane by comparing each estimated physical parameter with a predefined threshold value; and
recommending membrane cleaning when the predicted physical parameter of the membrane exceeds the predefined threshold value.

12. The method of claim 1, comprising:
validating the mathematical model based on a comparison of estimated time varying physical parameters of the membrane with the stored real time measurements of the plant process variable.