# Patent

## Title
System and process for determining and analysing surface property parameters of substance based on kinetic method

## Abstract

A system and process for determining and analyzing surface property parameters of a substance based on kinetic method is provided. The system comprises a sample processing system and a detection system. The sample processing system includes a reactor (3), a collector for liquid to be tested (5), and a container for liquid to be tested (6). The detection system includes a detecting electrode (13), a concentration and activity operator, a kinetic data processor, a surface property operation module, and a result output module. The process comprises: having the substance to be tested to be treated with an electrolyte solution, measuring activity of liquid to be tested upon reaction at a pre-set time interval, and processing with the kinetic data processor and the surface property operation module, so as to obtain surface property parameters of the substance to be tested. The present invention adopts kinetic method of ion exchange to overcome the issues associated with long reaction equilibrium time, not easy to determine equilibrium and not easy to accurately determine the value of m. Five parameters of the substance surface property can be calculated with only intercept and gradient of the linear regression equation obtained from the kinetic data, and the total amount of surface charge can also be directly determined.

## Background
This application is a National Stage Application of PCT/CN2012/070635, filed 20 Jan. 2012, which claims benefit of Serial No. 201110378979.2, filed 24 Nov. 2011 in China and Serial No. 201110378966.5, filed 24 Nov. 2011 and which applications are incorporated herein by reference. To the extent appropriate, a claim of priority is made to each of the above disclosed applications.
TECHNICAL FIELD
The present invention relates to fields of colloid chemistry, interface chemistry, chemistry and chemical engineering, environmental protection, and material science and technology, and specially relates to a system and process for determining and analyzing surface property parameters of a substance based on kinetic method.
BACKGROUND
The surface properties of a substance, such as surface potential, total amount of electric charge on the surface, surface charge density, electric field intensity on the surface and specific surface area, etc., not only are widely used in science researches in the fields of colloid and interface chemistry, material science, life science, soil science, ecological and environmental science and the like, but also are widely used in chemical engineering fields such as papermaking, cement, ceramics, chemical mechanical polishing, coal slurry, coating, cosmetics, food industry, mixing and dispersion system. Thus, determination on the surface property parameters of a substance appears particularly important.
Currently, for determining the total amount of electric charge on the surface, ion adsorption indicator and potentiometric titration method are commonly used. When using ion adsorption, one first needs to know the amount contributed to electrostatic adsorption from the total adsorption amount at H+ or OH−. However, the adsorption amount of involved in electrostatic adsorption cannot be predicted, due to H+ and OH− also involve adsorption of chemical bond. Thus, this method cannot determine the total amount of surface electric charge of a system containing variable charge under any pH value, any electrolyte concentration and any temperature. Further, the potentiometric titration method not only is not suitable for determining the total charge amount of a system containing permanent charge, but also is always questioned on its reliability, even for a variable charge system. Thus, so far there is no common method suited for determining total amount of substance surface charge in different conditions and different systems.
Currently, one method for determining a surface charge density of a substance is based on the following formula:




  
    
      σ
      0
    
    =
    
      
        T
        C
      
      S
    
  
  ,




Wherein, σ0 is surface charge density, Tc is total amount of surface charge, and S is specific surface area.
Since the parameter of total amount of surface charge is required, the problem with determination of total amount of surface charge must also exist in the determination of surface charge density. Further, in the above formula, the specific surface area is also required. However, the result of the specific surface area may greatly vary, if different methods may be used to determine the specific surface area. Thus, for a method that determines surface charge density based on the parameter of specific surface area, the reliability of its result is hard to expect.
Conventionally, a second method for determining surface charge density is: upon obtaining the surface potential of the substance, correlation formula of Gouy-Chapman is used to indirectly obtain the surface charge density. However, up until now no widely applied method is available to accurately determine the surface potential. Thus, there are still difficulties to apply such method to determine surface charge density.
Currently, electric field intensity can be determined based on the below formula:




  
    
      E
      0
    
    =
    
      
        
          4
          ⁢
          π
        
        ɛ
      
      ⁢
      
        σ
        0
      
    
  
  ,




Wherein, E0 is the surface electric field intensity, ∈ is the media dielectric constant, e.g., water has ∈=8.9×10−10 C2/Jdm. Due to depending on the surface charge density, same problem in determining surface charge density exists in determining surface electric field intensity.
In the prior art, there are a number of method for determining a substance specific surface area, such as commonly used inert gas adsorption method, ion negative adsorption method, glycol ethylene ether adsorption method, or glycerin adsorption method, etc. However, the results vary substantially by using these different methods for same substance.
In the present invention, substance surface potential is defined as the potential on the original surface of diffusion layer or OHP (outer Helmholtz surface) in the double electric layer. In the prior art, the method for determining substance surface potential includes: charge density method, negative adsorption method, positive adsorption method, dual-stage resonance (generation) method, pH indicator molecule method, fluorescence method, atomic force microscope method and Zeta potential method, and the like. All these methods have their own limitations. The charge density method, negative adsorption method, positive adsorption method and dual-stage resonance method are only suitable for determining surface potential of a constant charge sample of a single electrolyte system under neutral condition. The pH indicator molecule method, fluorescence method and atomic force microscope method may damage the condition of the substance surface per se, and thus the reliability of the result is difficult to say. Zeta potential method does not measure the surface potential. Rather, it measures the potential on the shear plane (or sliding surface) during electrophoresis, and the shear plane is often remote from the defined surface. The Zeta potential can be measured under different pH, electrolytes and temperature conditions. Since there is no method that can be widely used and accurately determine the surface potential under various conditions, Zeta potential has to be used as a substitute for surface potential. However, numerous studies over the years have shown that using Zeta potential to determine surface potential often is qualitative. In addition, Zeta potential method demands a very severe condition on the object to be tested, i.e., it requires that the particulate density of colloid suspension cannot be too high, while the size of particulate cannot be too large. Even for the new Zeta potentiometer (model Zetaprobe) by Colloidal Dynamics Inc., USA, the highest particulate density is only 60% (volume density). Thus, it is impossible to obtain a system with higher density, or “original state” measurement of solid particulate material.
Currently, Li Hang, et al. proposes a new method, i.e., with ion exchange equilibrium experiment and electrode method to realize combined measurement of five parameters including surface potential, specific surface area, surface charge density, total amount of surface charge and surface electric field intensity. This method marks a breakthrough in the field of determining surface property parameters of a substance. However, this method also has some weakness: (1) This method may need quite long time to determine surface property, due to the method is based on ion exchange equilibrium experiment, and the ion exchange equilibrium in actual system often requires considerable time. (2) Due to different materials have different surface charge amount and charge density, the time required for ion exchange equilibrium may vary significantly. Thus, in practice, it may be difficult to control the equilibrium state with this method. (3) In this method, it is required to calibrate three parameters, i.e., βA, βB and m with standard sample, in which βA, βB may be easy to be accurately calibrated, but accurate calibration of m may have difficulty. (4) This method may first need to use strong acid to treat the sample to make the tested sample to be H+ saturated, and thus may bring certain change to the substance surface property. In addition, the electrostatic bond between H+ ion and surface is far greater than that of ions such as Ca2+ and Na+, which may cause difficulty to reach exchange equilibrium.
Thus, there is a need for a system and method for determining and analyzing surface property parameters of a substance, which can overcome the shortcomings such as long equilibrium time, not easy to determine equilibrium and not easy to accurately determine m.
SUMMARY
In view of above, the present invention provides a system and a method for determining and analyzing surface property parameters of a substance, which can overcome the shortcomings such as long equilibrium time, not easy to determine equilibrium and not easy to accurately determine the value of m. With kinetic measurement of ion exchange, the kinetic data obtained within short time is used to accurately predict equilibrium state, to overcome the difficulty of long equilibrium time and not easy to determine equilibrium, to directly determine the total amount of surface charge. Thus, in determining the specific surface area, the parameter of m would not be needed. This can overcome the difficulty in determining m. In addition, there is no need to make the sample to be H+ saturated, to overcome the relevant issues brought by H+ saturation.
One objective of the present invention is to provide a system for determining and analyzing surface property parameters of a substance based on a kinetic method. A further objective of the present invention is to provide a method for determining and analyzing surface property parameters of a substance based on a kinetic method.
One objective of the present invention can be achieved with the following technical solutions:
A system for determining and analyzing surface property parameters of a substance based on a kinetic method, comprising: a sample processing system and a detection system. The sample processing system includes a reactor, a collector for liquid to be tested, and a container for liquid to be tested. The reactor is used to receive substance to be tested and reacting liquid. The reactor and the collector for liquid to be tested are connected via a tube. The collector for liquid to be tested and the container for liquid to be tested are connected via a tube. The reacting liquid in the reactor flows into the collector for liquid to be tested. The liquid in the collector for liquid to be tested flows into the container for liquid to be tested.
The detection system includes a detecting electrode, a concentration and activity operator, a kinetic data processor, a surface property operation module, and a result output module. The detecting electrode is used to detect ion activity of solution in a sample containing unit. The concentration and activity operator is used to calculate concentration of the solution based on the ion activity obtained from the detecting electrode. The kinetic data processor is used to process and store data obtained from detecting and establish a linear regression equation, so as to obtain intercept and gradient of the line. The surface property operation module is used to calculate a surface property parameter of the substance based on the detected data. The result output module is used to output the surface property parameter of the substance. The detecting electrode, the concentration and activity operator, the kinetic data processor, the surface property operation module and the result output module are connected in order.
Yet, the surface property operation module includes a surface charge total amount operator, a surface potential operator, a surface charge density operator, a surface electric field intensity operator, and a specific surface area operator. The kinetic data processor is connected with the surface charge total amount operator and the surface potential operator. The surface potential operator is connected with the surface charge density operator. The surface charge density operator is connected with the surface electric field intensity operator. The surface charge density operator and the surface charge total amount operator are connected with the specific surface area operator. The result output module is connected with the surface charge total amount operator, the surface potential operator, the surface charge density operator, the surface electric field intensity operator, and the specific surface area operator, respectively.
Yet, the kinetic data processor establishes the linear regression equation through the following steps, to obtain the intercept and gradient of the line:
S1: y-coordinate: y=[NB(ti+1)−NB(ti)]/(ti+1−ti),

    
    
        x-coordinate: x=NB(ti)+0.5[NB(ti+1)−NB(ti)];
    
    


S2: drawing a linear image based on positions of the x-coordinate and y-coordinate in a rectangular coordinate system;
S3: obtaining the intercept p and gradient q of the line;
wherein, NB(ti) indicates accumulated adsorption amount of ion Bm+, NB(ti) is the accumulated adsorption amount of the sample to ion Bm+, upon reaction for t=ti minutes, in which t0 is a pre-set interval, and ti is reaction time.
Yet, the system includes: a millivoltmeter, a constant flow pump, a liquid container, a constant temperature chamber and a thermostat (controller for constant temperature). The input of the millivoltmeter is connected with the detecting electrode, and the output of the millivoltmeter is connected with the concentration and activity operator. The constant flow pump is positioned between the reactor and the collector for liquid to be tested, for controlling the liquid flow velocity from the reactor to the collector for liquid to be tested. The system further includes: a first control valve, a second control valve, a third control valve, an automatic controller for valve, and a negative pressure suction pump. The first control valve is positioned between the collector for liquid to be tested and the container for liquid to be tested. The container for liquid to be tested is provided with the second control valve, which is connected with the negative pressure pump. The container for liquid to be tested is further provided with the third control valve that is connected with the negative pressure pump. The third control valve and the negative pressure pump are further provided with a waste liquid discharge outlet, which is used for discharging solution in the container for liquid to be tested. The automatic controller for valve is connected with the first control valve, the second control valve and the third control valve, respectively. The thermostat is connected with the constant flow pump. The liquid container is connected with the reactor. The sample processing system is positioned within the constant temperature chamber.
Yet, the surface charge total amount operator calculates the surface charge total amount of the substance to be tested by the following equation:




  
    
      S
      ⁢
      
        
      
      ⁢
      C
      ⁢
      
        
      
      ⁢
      N
    
    =
    
      
        -
        m
      
      ⁢
      
        
          p
          0
        
        
          q
          0
        
      
    
  
  ,




Wherein, SCN is the total amount of surface charge, m is the valence of ion Bm+, p0 is the intercept of the line, and q0 is the gradient of the line.
The surface potential operator calculates the potential on the surface of the substance to be tested by the following equation:




  
    
      φ
      0
    
    =
    
      
        -
        
          RT
          
            
              (
              
                
                  m
                  ⁢
                  
                    
                  
                  ⁢
                  
                    β
                    B
                  
                
                -
                
                  n
                  ⁢
                  
                    
                  
                  ⁢
                  
                    β
                    A
                  
                
              
              )
            
            ⁢
            F
          
        
      
      ⁢
      
        
          ln
          [
          
            
              
                a
                A
                0
              
              
                a
                B
                0
              
            
            ⁢
            
              p
              
                (
                
                  
                    q
                    ⁢
                    
                      
                        
                        
                          S
                          ⁢
                          
                            
                          
                          ⁢
                          C
                          ⁢
                          
                            
                          
                          ⁢
                          N
                        
                        
                      
                      n
                    
                  
                  +
                  
                    
                      m
                      n
                    
                    ⁢
                    p
                  
                
                )
              
            
          
          ]
        
        2
      
    
  
  ,




Wherein, φ0 is the surface potential of the substance to be tested, βB and βA are effective charge coefficient of ion Bm+ and ion An+, respectively, F indicates Faraday constant, R indicates gas constant, T indicates temperature, m indicates valence of ion Bm+, n indicates valence of An+, p indicates intercept, q indicates gradient, aA0 and aB0 are respective activities of ion An+ and ion Bm+ in the bulk solution when reaction reaches equilibrium.
The surface charge density operator calculates the surface charge density of the substance to be tested by the following equation:




  
    
      σ
      0
    
    ≈
    
      
        sign
        ⁡
        
          (
          
            φ
            0
          
          )
        
      
      ⁢
      
        
          
            
              ɛ
              ⁢
              
                
              
              ⁢
              RT
            
            
              2
              ⁢
              π
              ⁢
              
                
              
              ⁢
              
                F
                2
              
            
          
          [
          
            
              
                a
                B
                0
              
              (
              
                
                  e
                  
                    -
                    
                      
                        m
                        ⁢
                        
                          
                        
                        ⁢
                        
                          β
                          B
                        
                        ⁢
                        F
                        ⁢
                        
                          
                        
                        ⁢
                        
                          φ
                          0
                        
                      
                      RT
                    
                  
                
                -
                1
              
              )
            
            +
            
              
                a
                A
                0
              
              (
              
                
                  e
                  
                    -
                    
                      
                        n
                        ⁢
                        
                          
                        
                        ⁢
                        
                          β
                          A
                        
                        ⁢
                        F
                        ⁢
                        
                          
                        
                        ⁢
                        
                          φ
                          0
                        
                      
                      RT
                    
                  
                
                -
                1
              
              )
            
          
          ]
        
      
    
  
  ,




Wherein, σ0 is the surface charge density, φ0 is the surface potential, and c is the dielectric constant of water.
The surface electric field intensity operator calculates the surface electric field intensity of the substance to be tested by the following equation:




  
    
      E
      0
    
    =
    
      
        
          4
          ⁢
          π
          ⁢
          
            
          
          ⁢
          F
        
        ɛ
      
      ⁢
      
        σ
        0
      
    
  
  ,




Wherein, E0, is the surface electric field intensity of the substance to be tested.
The specific surface area operator calculates the specific surface area of the substance to be tested by the following equation:




  
    S
    =
    
      
        S
        ⁢
        
          
        
        ⁢
        C
        ⁢
        
          
        
        ⁢
        N
      
      
        σ
        0
      
    
  
  ,




Wherein, S is the specific surface area of the substance to be tested.
Another objective of the present invention can be achieved by the following technical solutions:
A method for determining and analyzing surface property parameters of a substance based on a kinetic method, comprising the following steps:
S1: conducting An+ saturation to a surface of the substance to be tested with AXn electrolyte solution having known concentration;
wherein, the AXn electrolyte solution represents an electrolyte solution of metallic cation with charge amount of n, and ion An+ represents metallic cation with charge amount of n;
S2: having a solution of BXm electrolyte with known concentration to flow through the surface of the substance to be tested, to have ion An+ and ion Bm+ on the surface of the substance to be tested to perform fusion replacement;
wherein, solution of BXm electrolyte represents an electrolyte solution of metallic cation with charge amount of m, and ion Bm+ represents metallic cation with charge amount of m;
S3: collecting flow-out liquid that flows through the surface of the substance to be tested at pre-set interval;
S4: measuring activity of ion Bm+ in the flow-out liquid, and calculating concentration of ion Bm+;
S5: calculating accumulated adsorption amount of ion B+ according to the following equation by using the concentration of ion B+:




  
    
      
        N
        B
      
      ⁡
      
        (
        
          t
          i
        
        )
      
    
    =
    
      
        1
        m
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          
            
              t
              i
            
            /
            
              t
              0
            
          
        
        ⁢
        
          V
          ⁡
          
            [
            
              
                f
                
                  B
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
              -
              
                
                  f
                  B
                
                ⁡
                
                  (
                  
                    t
                    i
                  
                  )
                
              
            
            ]
          
        
      
    
  
  ,




wherein, NB(ti) is the accumulated adsorption amount of the sample to ion Bm+ upon reaction for t=ti minutes, m is sample mass, V is the volume of solution collected within the pre-set interval, fB0 is the concentration of solution before ion Bm+ flows into the sample, fB(ti) is the concentration of ion Bm+ in solution collected by ith time, t0 is the pre-set interval, and ti is reaction time;
S6: using NB(ti) (accumulated adsorption amount of ion Bm+), with y-coordinate as y=[NB(ti+1)−NB(ti)]/(ti+1−ti) and corresponding x-coordinate as x=NB(ti)+0.5[NB(ti+1)−NB(ti)], to draw a line image based on positions of the x-coordinate and y-coordinate in a rectangular coordinate system, and to obtain intercept of the line as p0 and gradient of the line as q0;
S7: calculating the surface charge total amount of the sample with intercept p0 and gradient q0 of the line according to the below equation:




  
    
      S
      ⁢
      
        
      
      ⁢
      C
      ⁢
      
        
      
      ⁢
      N
    
    =
    
      
        -
        m
      
      ⁢
      
        
          p
          0
        
        
          q
          0
        
      
    
  
  ,




wherein, SCN is the total amount of the surface charge, and m is valence of ion Bm+.
Yet, the step S4 comprises the following sub-steps to calculate concentration of ion Bm+ by using iterative algorithm:
S41: having the solution of BXm electrolyte with known concentration to flow through the surface of the substance to be tested, to have ion An+ and ion Bm+ on the surface of the substance to be tested to perform fusion replacement: mSoil−A+nBXm=nSoil−B+mAXn; obtaining a solution having An+, Bm+ and X−, to calculate concentration of ion An+ and concentration of ion Bm+:




  
    
      C
      A
    
    =
    
      
        (
        
          m
          n
        
        )
      
      ⁢
      
        C
        B
      
    
  
  ,




CX=mnCB;

wherein, CA represents concentration of ion An+, CB represents concentration of ion Bm+, CX represents concentration of ion X−, Soil−A represents soil saturated with ion An+, and Soil−B represents soil saturated with ion Bm+;
S42: measuring activity of ion Bm+ in aqueous solution, with a measurement of ion Bm+ activity as aB, using the ion activity as an initial concentration of ion Bm+ in the first iterative operation, i.e., CB(0)=aB. Likewise, the initial concentration of ion An+ is:




  
    C
    A
    
      (
      0
      )
    
  
  =
  
    
      (
      
        m
        n
      
      )
    
    ⁢
    
      
        a
        B
      
      .
    
  




S43: calculating concentration of ion X− in the first iterative operation to be CX(0)=mnaB, and ionic strength I(1) in the first iterative operation to be:




  
    
      I
      
        (
        1
        )
      
    
    =
    
      
        
          1
          2
        
        ⁢
        
          
            ∑
            i
          
          ⁢
          
            
              c
              i
              
                (
                0
                )
              
            
            ⁢
            
              Z
              i
              2
            
          
        
      
      =
      
        
          
            1
            2
          
          ⁢
          
            (
            
              
                
                  m
                  2
                
                ⁢
                
                  a
                  B
                
              
              +
              
                mna
                B
              
              +
              
                mna
                B
              
            
            )
          
        
        =
        
          
            1
            2
          
          ⁢
          
            (
            
              
                m
                2
              
              +
              
                2
                ⁢
                mn
              
            
            )
          
          ⁢
          
            a
            B
          
        
      
    
  
  ,




wherein, I(1) represents ion strength coefficient in the first iterative operation, Zi represents valence of ion i, ci(0) represents initial concentration of ion i in the iterative operation.
S44: obtaining initial activity coefficient of ion Bm+ in the system from Debye-Hückel limiting formula:




  
    
      γ
      B
      
        (
        1
        )
      
    
    =
    
      exp
      (
      
        -
        
          
            6030.2
            ×
            
              T
              
                -
                
                  3
                  2
                
              
            
            ⁢
            
              
              
                Z
                +
                2
              
              
            
            ⁢
            
              
                I
                
                  (
                  1
                  )
                
              
            
          
          
            1
            +
            
              
                I
                
                  (
                  1
                  )
                
              
            
          
        
      
      )
    
  
  ;




wherein, γB(1) represents initial activity coefficient of ion Bm+, Z+ represents charge amount of positive ion, I(1) represents ion strength of the first iterative operation, and T represents temperature;
S45: calculating concentration cB(1) of ion Bm+ after the first iteration:




  
    
      c
      B
      
        (
        1
        )
      
    
    =
    
      
        a
        B
      
      
        γ
        B
        
          (
          1
          )
        
      
    
  
  ;




S46: calculating concentration CX(1) of ion X after the first iteration:

CX(1)=mncB(1);

S47: respectively calculating ion strength and activity coefficients of the second iteration as:



  
    
      
        I
        
          (
          2
          )
        
      
      =
      
        
          
            1
            2
          
          ⁢
          
            
              ∑
              i
            
            ⁢
            
              
                c
                i
                
                  (
                  1
                  )
                
              
              ⁢
              
                Z
                i
                2
              
            
          
        
        =
        
          
            1
            2
          
          ⁢
          
            (
            
              
                m
                2
              
              +
              
                2
                ⁢
                mn
              
            
            )
          
          ⁢
          
            c
            B
            
              (
              1
              )
            
          
        
      
    
    ;
  




  
    
      
        γ
        B
        
          (
          2
          )
        
      
      =
      
        exp
        (
        
          -
          
            
              6030.2
              ×
              
                T
                
                  -
                  
                    3
                    2
                  
                
              
              ⁢
              
                
                
                  Z
                  +
                  2
                
                
              
              ⁢
              
                
                  I
                  
                    (
                    2
                    )
                  
                
              
            
            
              1
              +
              
                
                  I
                  
                    (
                    2
                    )
                  
                
              
            
          
        
        )
      
    
    ;
  



wherein, γB(2) represents ion activity coefficient of the second iteration, cB(1) represents concentration of ion Bm+ after the first iteration, Zi represents valence of ion i, and I(2) represents ion strength of the second iteration;
S48: calculating concentration cB(2) of ion Bm+ after the second iteration:




  
    
      c
      B
      
        (
        2
        )
      
    
    =
    
      
        a
        B
      
      
        γ
        B
        
          (
          2
          )
        
      
    
  
  ;




S49: concentration of ion X after the second iterative operation shall be:

CX(2)=mncB(2);

S410: repeatedly iterating n times of the above steps S41-S49, and calculating according to the following equation:



  
    
      
        I
        
          (
          n
          )
        
      
      =
      
        
          
            1
            2
          
          ⁢
          
            
              ∑
              i
            
            ⁢
            
              
                c
                i
                
                  (
                  
                    n
                    -
                    1
                  
                  )
                
              
              ⁢
              
                Z
                i
                2
              
            
          
        
        =
        
          
            1
            2
          
          ⁢
          
            (
            
              
                m
                2
              
              +
              
                2
                ⁢
                mn
              
            
            )
          
          ⁢
          
            c
            B
            
              (
              
                n
                -
                1
              
              )
            
          
        
      
    
    ;
  




  
    
      
        γ
        B
        
          (
          n
          )
        
      
      =
      
        exp
        (
        
          -
          
            
              6030.2
              ×
              
                T
                
                  -
                  
                    3
                    2
                  
                
              
              ⁢
              
                
                
                  Z
                  +
                  2
                
                
              
              ⁢
              
                
                  I
                  
                    (
                    n
                    )
                  
                
              
            
            
              1
              +
              
                
                  I
                  
                    (
                    n
                    )
                  
                
              
            
          
        
        )
      
    
    ;
  



when satisfying (I(n)−I(n−1))/I(n)<0.001, iteration is stopped, and the resulting concentration cB(n) is:




  
    
      c
      B
      
        (
        n
        )
      
    
    =
    
      
        a
        B
      
      
        γ
        B
        
          (
          n
          )
        
      
    
  
  ;




wherein, I(n) represents ion strength of nth iteration, I(n−1) represents ion strength of (n−1)th iteration, cB(n) represents concentration of ion Bm+ after nth iteration, and γB(n) represents ion activity coefficient of nth iteration.
Yet, the below steps may be included following step S7:
S8: calculating ion strength I of the flow solution of mixed electrolyte having ion An+ and ion Bm+ according to the following equation:

I=1⁄2[(n+n2)fA0+(m+m2)fB0]

wherein, fA0 is the solution concentration before ion An+ flows into the surface of the substance to be tested, and fB0 is the solution concentration before ion B+ flows into the surface of the substance to be tested;
S9: Substituting the ion strength of the flow liquid in Davies equation, to calculate ion An+ activity coefficient γA0 and ion Bm+ activity coefficient γB0;
S10: calculating activity of ion An+ and ion Bm+ according to the following equation:

aA0=fA0γA0 aB0=fB0γB0 

wherein, γA0 and γB0 are respective activity coefficient of ion An+ and ion Bm+ in bulk solution when reaction reaches equilibrium; aA0 and aB0 are respective activity of ion An+ and ion Bm+ in bulk solution when reaction reaches equilibrium.
The method may include the following steps after step S10:
S11: under a given pH condition, having a mixed electrolyte solution by AXn electrolyte solution with known concentration fA0 and BXm electrolyte solution with known concentration fB0 to flow through the surface of the substance to be tested, to have ion An+ and ion Bm+ on the surface of the substance to be tested to perform fusion replacement, and collecting the flow-out solution that flows through the surface of the substance to be tested in pre-set interval.
S12: repeating the steps of S1 to S6;
S13: using NB(ti) (accumulated adsorption amount of ion Bm+), with y-coordinate as y=[NB(ti+1)−NB(ti)]/(ti+1−ti) and corresponding x-coordinate as x=NB(ti)+0.5[NB(ti+1)−NB(ti)], to draw a line image based on positions of x-coordinate and y-coordinate in a rectangular coordinate system, and to obtain intercept p and gradient q of the line.
The method further comprises the following steps after step S13:
S14: using the total charge amount obtained from step S7, the activity obtained from step S10, and the intercept p and gradient q of the line obtained from step S13, to substitute into the below equation to calculate the electric potential on the surface of the substance to be tested:




  
    
      φ
      0
    
    =
    
      
        -
        
          RT
          
            
              (
              
                
                  m
                  ⁢
                  
                    
                  
                  ⁢
                  
                    β
                    B
                  
                
                -
                
                  n
                  ⁢
                  
                    
                  
                  ⁢
                  
                    β
                    A
                  
                
              
              )
            
            ⁢
            F
          
        
      
      ⁢
      
        
          ln
          [
          
            
              
                a
                A
                0
              
              
                a
                B
                0
              
            
            ⁢
            
              p
              
                (
                
                  
                    q
                    ⁢
                    
                      
                        
                        SCN
                        
                      
                      n
                    
                  
                  +
                  
                    
                      m
                      n
                    
                    ⁢
                    p
                  
                
                )
              
            
          
          ]
        
        2
      
    
  
  ,




wherein, φ0 is the surface potential of the substance to be tested, βB and βA are effective charge coefficients of ion Bm+ and ion An+, F represents Faraday constant, R represents gas constant, and T represents temperature;
S15: substituting the electric potential on the surface of the substance to be tested in the following equation to calculate the surface charge density of the substance to be tested:




  
    
      σ
      0
    
    ≈
    
      
        sign
        ⁡
        
          (
          
            φ
            0
          
          )
        
      
      ⁢
      
        
          
            
              ɛ
              ⁢
              
                
              
              ⁢
              RT
            
            
              2
              ⁢
              π
              ⁢
              
                
              
              ⁢
              
                F
                2
              
            
          
          ⁡
          
            [
            
              
                
                  a
                  B
                  0
                
                ⁡
                
                  (
                  
                    
                      e
                      
                        -
                        
                          
                            m
                            ⁢
                            
                              
                            
                            ⁢
                            
                              β
                              B
                            
                            ⁢
                            F
                            ⁢
                            
                              
                            
                            ⁢
                            
                              φ
                              0
                            
                          
                          RT
                        
                      
                    
                    -
                    1
                  
                  )
                
              
              +
              
                
                  a
                  A
                  0
                
                ⁡
                
                  (
                  
                    
                      e
                      
                        -
                        
                          
                            n
                            ⁢
                            
                              
                            
                            ⁢
                            
                              β
                              A
                            
                            ⁢
                            F
                            ⁢
                            
                              
                            
                            ⁢
                            
                              φ
                              0
                            
                          
                          RT
                        
                      
                    
                    -
                    1
                  
                  )
                
              
            
            ]
          
        
      
    
  
  ,




wherein, σ0 is the surface charge density (mol/dm2), ∈ is the dielectric constant of medium, in which water has ∈=8.9×10−10 C2/(J·dm);
S16: substituting the surface charge density of the substance to be tested in the following equation to calculate the surface electric field intensity of the substance to be tested:




  
    
      E
      0
    
    =
    
      
        
          4
          ⁢
          π
          ⁢
          
            
          
          ⁢
          F
        
        ɛ
      
      ⁢
      
        σ
        0
      
    
  
  ,




wherein, E0 represents the surface electric field intensity of the substance to be tested;
S17: substituting the total charge amount and the surface charge density of the substance to be tested in the following equation to calculate the specific surface area of the sample:




  
    S
    =
    
      SCN
      
        σ
        0
      
    
  
  ,




wherein, S represents the specific surface area of the substance to be tested;
S18: outputting the measurement of the surface parameters of the substance to be tested, the output measurement of surface parameters of the substance to be tested includes: measurements of the surface potential of the substance to be tested, the surface charge total amount of the substance to be tested, the surface charge density of the substance to be tested, the surface electric field intensity of the substance to be tested, and the specific surface area of the substance to be tested.
Advantages of the present invention include: the present invention adopts an ion exchange kinetic method to overcome the difficulties such as long reaction equilibrium time, not easy to determine equilibrium and not easy to accurately determine the value of m, when measuring/determining surface property parameters of a substance. With kinetic data obtained in short time to accurately calculate, the above shortcomings can be overcome. Five parameters of the substance surface property can be calculated with only intercept and gradient of the linear regression equation obtained from the kinetic data. Further, with kinetic method, the total amount of surface charge can be directly determined, and thus the parameter of m no longer is needed in determining the specific surface area, which solves the issue occurred in determining m. In addition, this method does not need the sample to be H+ saturated, which solves the relevant issues brought by H+ saturation, thereby solving the issues associated with change in substance surface property brought by treating the sample with strong acid and hard to reach ion exchange equilibrium in the Fr saturated sample. Further, when the activity of ion Bm+ in the equilibrium liquid is much greater than that of H+ in the solution, the electrostatic adsorption of H+ can be negligible. This ensures that the present invention can be a common method suitable for any pH and any substance type.
Other advantages, objectives and features will be described in the following specification, and to some extent, which will be obvious to people skilled in the art based on studies, or which can be taught by implementing the present invention. The objectives and other advantages of the present invention can be achieved and obtained from the description, claims and drawings.



BRIEF DESCRIPTION OF DRAWINGS
To make the objectives, technical solutions and advantages of the present invention more clear, detailed description will be made to the present invention in combination with drawings, wherein:
FIG. 1 is structural diagram showing a system for analyzing surface property parameters of a substance based on kinetic method, according to an embodiment of the present invention.
FIG. 2 is structural diagram showing a detection system in the system for analyzing surface property parameters of a substance based on kinetic method, according to an embodiment of the present invention.
FIG. 3 is flow chart of a system for analyzing surface property parameters of a substance based on kinetic method, according to an embodiment of the present invention.



In the drawings, 1—constant temperature chamber, 2—liquid container, 3—reactor, 4—constant flow pump, 5—collector for liquid to be tested, 6—container for liquid to be tested, 7—first control valve, 8—second control valve, 9—third control valve, 10—thermostat (controller for constant temperature), 11—automatic controller for valve, 12—negative pressure pump, 13—detecting electrode, 14—detection system, 15—waste liquid discharge outlet.
DETAILED DESCRIPTION
Preferred embodiments will be described in detail in combination with the accompanying drawings. It shall be understood that the preferred embodiments are for illustrative only, and not for limiting the protection scope of the present invention.
FIG. 1 is structural diagram showing a system for analyzing surface property parameters of a substance based on kinetic method in accordance with an embodiment of the present invention. As shown in the drawing, a system for analyzing surface property parameters of a substance based on kinetic method is provided. The system comprises a sample processing system and a detection system. The sample processing system includes a reactor 3, a collector 5 for liquid to be tested, a container 6 for liquid to be tested, a millivoltmeter, a constant flow pump 4, a liquid container 2, a constant temperature chamber 1 and a thermostat (constant temperature controller) 10. The reactor 3 is used for containing the substance to be tested and the reaction liquid. The reactor 3 is connected with the collector 5 for liquid to be tested via a conduit. The collector 5 for liquid to be tested is connected with the container 6 for liquid to be tested via a conduit. The reaction liquid in the reactor 3 flows into the collector 5 for liquid to be tested. The liquid in the collector 5 for liquid to be tested flows into the container 6 for liquid to be tested. The input port of the millivoltmeter is connected with the detecting electrode 13, and the output port of the millivoltmeter is connected with the concentration and activity operator. The constant flow pump 4 is positioned between the reactor 3 and the collector 5 for liquid to be tested, for controlling liquid flow velocity from the reactor 3 to the collector 5 for liquid to be tested. The system further includes a first control valve 7, a second control valve 8, a third control valve 9, an automatic controller 11 for valve, and a negative pressure pump 12. The first control valve 7 is positioned between the collector 5 for liquid to be tested and the container 6 for liquid to be tested. The container 6 for liquid to be tested is provided with the second control valve 8 which is connected with the negative pressure pump 12. The container 6 for liquid to be tested is further provided with the third control valve 9 that is connected with the negative pressure pump 12. The third control valve 9 and the negative pressure pump 12 are provided with a waste liquid discharge outlet 15 that is used for draining solution in the container 6 for liquid to be tested. The automatic controller 11 for valve is connected with the first control valve 7, the second control valve 8 and the third control valve 9, respectively. The thermostat (controller) 10 is connected with the constant flow pump 4. The liquid container 2 is connected with the reactor 3. The sample processing system is disposed within the constant temperature chamber 1.
FIG. 2 is structural diagram showing a detection system in the system for analyzing surface property parameters of a substance based on kinetic method in accordance with one embodiment of the present invention. As shown in FIG. 2, a detection system 14 for is provided. The detection system 14 comprises a detecting electrode 13, a concentration and activity operator, a kinetic data processor, a surface property operation module, and a result output module. The detecting electrode is used to detect ion activity of the solution in the sample containing unit. The concentration and activity operator is used to calculate concentration of the solution based on the ion activity obtained from the detecting electrode. The kinetic data processor is used to process and store the detected data and establish a linear regression equation, so as to obtain intercept and gradient of the line. The surface property operation module is used to calculate the substance surface property parameters based on the detected data. The result output module is used to output the substance surface property parameters. The detecting electrode, the concentration and activity operator, the kinetic data processor, the surface property operation module, and the result output module are connected in order.
The surface property operation module comprises a surface charge total amount operator, a surface potential operator, a surface charge density operator, a surface electric field intensity operator, and a specific surface area operator. The kinetic data processor is connected with the surface charge total amount operator and the surface potential operator. The surface potential operator is connected with the surface charge density operator. The surface charge density operator is connected with the surface electric field intensity operator. The surface charge density operator and the surface charge total amount operator are connected with the specific surface area operator. The result output module is connected with the surface charge total amount operator, the surface potential operator, the surface charge density operator, the surface electric field intensity operator, and the specific surface area operator, respectively.
FIG. 3 is a flow chart of a system for analyzing surface property parameters of a substance based on kinetic method in accordance with one embodiment of the present invention. As shown in FIG. 3, a kinetic method for determining surface property parameters of a substance is provided. The analyzing system operates as follows.
An+ saturation is performed to the substance to be tested with AXn electrolyte solution having known concentration. The An+ saturated substance to be tested then is disposed in the reactor. The constant flow pump is set to control the flowing liquid to flow at a desired velocity (e.g., 1 ml/min). Under a given pH condition, the electrolyte solution of another positive ion Bm+ having known concentration of charge amount m evenly flows through the sample in the reactor at a pre-set flow velocity, where fusion replacement occurs between An+ and Bm+. When the first drop of liquid flows out the reactor outlet, the automatic controller for valve starts the clock (at this moment, the first control valve 7, the second control valve 8 and the third control valve 9 are in automatic close state). Based on the set interval (e.g., 5 minute), the first control valve 7, the third control valve 9 and the negative pressure pump simultaneously open/turn on, to transfer the collected liquid from the collector for liquid to be tested to the container for liquid to be tested. Then, the first control valve 7 and the third control valve 9 are closed. The activity of Bm+ in the container for liquid to be tested is measured with a detector. Then, concentration of Bm+ is calculated with the concentration and activity operator in the detector by using iterative operation. The obtained concentration is stored in the “kinetic data processor” in the detector. Then, the second control valve 8 and the negative pressure pump will open/turn on, to drain the solution in the liquid container. The system automatically repeat collection of liquid from the reactor outlet, so as to obtain and store concentration of Bm+ in the collected liquid at each time interval (e.g., 5 minute). With concentration of Bm+ in the collected liquid obtained at each time interval, accumulated adsorption amounts of Bm+ at different time are calculated and stored in the kinetic data processor in the detector. By using the data of accumulated adsorption amounts obtained, data is automatically processed and stored in the kinetic data processor. All the operation results of individual operators in the detector are output, including the values of the surface potential of the substance to be tested, the surface charge total amount of the substance to be tested, the surface charge density of the substance to be tested, the surface electric field intensity of the substance to be tested, and the specific surface area of the substance to be tested.
A kinetic method for determining surface property parameters of a substance is provided, comprising the steps of:
S1: performing An+ saturation to a surface of the substance to be tested with AXn electrolyte solution having known concentration;
wherein, AXn electrolyte solution represents an electrolyte solution of metallic cation with a charge amount of n, and ion An+ represents metallic cation with a charge amount of n;
S2: having BXm electrolyte solution with known concentration to flow through the surface of the substance to be tested, to have fusion replacement occurred between An+ and Bm+ on the surface of the substance to be tested;
wherein, BXm electrolyte solution represents an electrolyte solution of metallic cation with a charge amount of m, and ion Bm+ represents metallic cation with a charge amount of m;
S3: collecting liquid flowing out from the surface of the substance to be tested at a pre-set time interval;
S4: determining activity of ion Bm+ in the flow-out liquid, and calculating concentration of ion Bm+ by using iterative operation;
S5: calculating accumulated adsorption amount of ion Bm+ with the concentration of ion Bm+ according to the following iterative equation:




  
    
      
        
          
            N
            B
          
          ⁡
          
            (
            
              t
              i
            
            )
          
        
        =
        
          
            1
            m
          
          ⁢
          
            
              ∑
              
                i
                =
                1
              
              
                
                  t
                  i
                
                /
                
                  t
                  0
                
              
            
            ⁢
            
              
            
            ⁢
            
              V
              ⁡
              
                [
                
                  
                    f
                    
                      B
                      ⁢
                      
                        
                      
                      ⁢
                      0
                    
                  
                  -
                  
                    
                      f
                      B
                    
                    ⁡
                    
                      (
                      
                        t
                        i
                      
                      )
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




wherein, NB(ti) is the accumulated adsorption amount of ion Bm+ by the sample upon reaction for t=ti minutes, with unit as mol/g; m is sample mass, with unit as g; V is the volume of solution collected within pre-set interval, with unit as liter, i.e., 1; fB0 is the solution concentration before ion Bm+ flows into the sample, with unit as mol/l; fB(ti) is concentration of ion Bm+ in the solution collected at ith time, with unit as mol/l; t0 is pre-set time interval, and ti is reaction time, with unit as min;
S6: using the accumulated adsorption amount NB(ti) of ion Bm+, with y-coordinate as y=[NB(ti+1)−NB(ti)]/(ti+1−ti) and corresponding x-coordinate as x=NB(ti)+0.5[NB(ti+1)−NB(ti)], to draw a line image based on positions of x-coordinate and y-coordinate in a rectangular coordinate system, to obtain intercept p0 and gradient q0 of the line;
S7: calculating the surface charge total amount of the sample based on the intercept p0 and gradient q0 of the line according to the following equation:




  
    
      
        SCN
        =
        
          
            -
            m
          
          ⁢
          
            
              p
              0
            
            
              q
              0
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




wherein, SCN is the total amount of surface charge, with a unit as mol/g; and, m is valence of ion B+.
S8: calculating ion strength of the mixed electrolyte flow liquid having ion An+ and ion Bm+ according to the following equation:

I=1⁄2[(n+n2)fA0+(m+m2)fB0]  (3)

wherein, fA0 is the solution concentration before ion An+ flows into the surface of the substance to be tested, and, fB0 is the solution concentration before ion Bm+ flows into the surface of the substance to be tested;
S9: Substituting the ion strength of the flow liquid in Davies equation, to calculate ion An+ activity coefficient γA0 and ion Bm+ activity coefficient γB0;
S10: calculating activity of ion An+ and ion Bm+ according to the following equation:

aA0=f0γA0 aB0=fB0γB0  (4)

wherein, γA0 and γB0 are respective activity coefficient of ion An+ and ion Bm+ in bulk solution when reaction reaches equilibrium; and, aA0 and aB0 are respective activity of ion An+ and ion Bm+ in bulk solution when reaction reaches equilibrium.
S11: under a given pH condition, having a mixed electrolyte solution of AXn electrolyte solution with a known concentration fA0 and BXm electrolyte solution with a known concentration fB0 to flow through the surface of the substance to be tested, to have fusion replacement occurred between ion An+ and ion Bm+ on the surface of the substance to be tested, and, collecting liquid flowing out from the surface of the substance to be tested at the pre-set time interval.
S12: repeating the steps of S1-S6;
S13: using the accumulated adsorption amount NB(ti) of ion B+, with y-coordinate as y=[NB(ti+1)−NB(ti)]/(ti+1−ti) and corresponding x-coordinate as x=NB(ti)+0.5[NB(ti+1)−NB(ti)], to draw a line image based on positions of x-coordinate and y-coordinate in a rectangular coordinate system, to obtain intercept p and gradient q of the line.
S14: substituting the charge total amount obtained from step S7, the activity obtained from step S10, and the intercept p and gradient q of the line obtained from step S13 in the following equation, to calculate electric potential on the surface of the substance to be tested:




  
    
      
        
          φ
          0
        
        =
        
          
            -
            
              RT
              
                
                  (
                  
                    
                      m
                      ⁢
                      
                        
                      
                      ⁢
                      
                        β
                        B
                      
                    
                    -
                    
                      n
                      ⁢
                      
                        
                      
                      ⁢
                      
                        β
                        A
                      
                    
                  
                  )
                
                ⁢
                F
              
            
          
          ⁢
          
            
              ln
              [
              
                
                  
                    a
                    A
                    0
                  
                  
                    a
                    B
                    0
                  
                
                ⁢
                
                  p
                  
                    (
                    
                      
                        q
                        ⁢
                        
                          
                            
                            SCN
                            
                          
                          n
                        
                      
                      +
                      
                        
                          m
                          n
                        
                        ⁢
                        p
                      
                    
                    )
                  
                
              
              ]
            
            2
          
        
      
    
    
      
        (
        5
        )
      
    
  




wherein, φ0 is the surface electric potential of the substance to be tested, βB and βA are effective charge coefficients of ion Bm+ and ion An+, F represents Faraday constant, R represents gas constant, and T represents temperature.
S15: substituting the surface electric potential of the substance to be tested in the following equation to calculate surface charge density of the substance to be tested:




  
    
      
        
          σ
          0
        
        ≈
        
          
            sign
            ⁡
            
              (
              
                φ
                0
              
              )
            
          
          ⁢
          
            
              
                
                  ɛ
                  ⁢
                  
                    
                  
                  ⁢
                  RT
                
                
                  2
                  ⁢
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  
                    F
                    2
                  
                
              
              ⁡
              
                [
                
                  
                    
                      a
                      B
                      0
                    
                    ⁡
                    
                      (
                      
                        
                          e
                          
                            -
                            
                              
                                m
                                ⁢
                                
                                  
                                
                                ⁢
                                
                                  β
                                  B
                                
                                ⁢
                                F
                                ⁢
                                
                                  
                                
                                ⁢
                                
                                  φ
                                  0
                                
                              
                              RT
                            
                          
                        
                        -
                        1
                      
                      )
                    
                  
                  +
                  
                    
                      a
                      A
                      0
                    
                    ⁡
                    
                      (
                      
                        
                          e
                          
                            -
                            
                              
                                n
                                ⁢
                                
                                  
                                
                                ⁢
                                
                                  β
                                  A
                                
                                ⁢
                                F
                                ⁢
                                
                                  
                                
                                ⁢
                                
                                  φ
                                  0
                                
                              
                              RT
                            
                          
                        
                        -
                        1
                      
                      )
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




wherein, σ0 is surface charge density (mol/dm2), s is medium dielectric constant, in which water has ∈=8.9×10−10 C2/J·dm.

S16: substituting the surface charge density of the substance to be tested in the following equation to calculate surface electric field intensity of the substance to be tested:




  
    
      
        
          E
          0
        
        =
        
          
            
              4
              ⁢
              π
              ⁢
              
                
              
              ⁢
              F
            
            ɛ
          
          ⁢
          
            σ
            0
          
        
      
    
    
      
        (
        7
        )
      
    
  




wherein, E0 represents the surface electric field intensity (V/dm) of the substance to be tested;
S17: substituting the surface charge total amount and the surface charge density of the substance to be tested in the following equation, to obtain specific surface area of the sample:




  
    
      
        S
        =
        
          SCN
          
            σ
            0
          
        
      
    
    
      
        (
        8
        )
      
    
  




wherein, S represents specific surface area (dm2/g) of the substance to be tested.
S18: outputting the value of the surface potential of the substance to be tested, the surface charge total amount of the substance to be tested, the surface charge density of the substance to be tested, the surface electric field intensity of the substance to be tested, and the specific surface area of the substance to be tested.
The above embodiments are only the preferred embodiments of the present invention, and shall not be used to limit the present invention. Apparently, people skilled in the art can modify the present invention, without departing from the spirit and scope of the present invention. As such, if such modifications fall in the scope of the appended claims and/or equivalent thereto, they are intended to be included in the present invention.

## Claims

1. A system for determining and analyzing surface property parameters of a substance based on kinetic method, comprising a sample processing system and a detection system;
the sample processing system including a reactor, a collector for liquid to be tested and a container for liquid to be tested; the reactor is used for containing the substance to be tested and reaction liquid; the reactor is connected with the collector for liquid to be tested via a conduit; the collector for liquid to be tested is connected with the container for liquid to be tested via a conduit; the reaction liquid in the reactor flows into the collector for liquid to be tested; the liquid in the collector for liquid to be tested flows into the container for liquid to be tested;
the detection system including a detecting electrode, a concentration and activity operator, a kinetic data processor, a surface property operation module, and a result output module; the detecting electrode is used for detecting an ion activity of solution in a sample container; the concentration and activity operator is used for calculating concentration of the solution based on the ion activity obtained from the detecting electrode; the kinetic data processor is used for processing and storing data obtained from detecting and establishes a linear regression equation to obtain intercept and gradient of a line; the surface property operation module is used for calculating the surface property parameters of the substance based on detected data; the result output module is used for outputting the surface property parameters of the substance; the detecting electrode, the operator concentration and activity, the kinetic data processor, the surface property operation module, and the result output module are connected in order.

2. The system for determining and analyzing surface property parameters of a substance based on kinetic method according to claim 1, wherein: the surface property operation module comprises a surface charge total amount operator, a surface potential operator, a surface charge density operator, a surface electric field intensity operator, and a specific surface area operator; the kinetic data processor is connected with the surface charge total amount operator and the surface potential operator; the surface potential operator is connected with the surface charge density operator; the surface charge density operator is connected with the surface electric field intensity operator; the surface charge density operator and the surface charge total amount operator are connected with the specific surface area operator; the result output module is respectively connected with the surface charge total amount operator, the surface potential operator, the surface charge density operator, the surface electric field intensity operator, and the specific surface area operator.

3. The system for determining and analyzing surface property parameters of a substance based on kinetic method according to claim 1, wherein: the kinetic data processor establishes the linear regression equation through the following steps, to obtain the intercept and gradient of the line:
S1: y-coordinate is y=[NB(ti+1)−NB(ti)]/(ti+1−ti),
x-coordinate is x=NB(ti)+0.5[NB(ti+1)−NB(ti)];

S2: drawing a line image based on x-coordinate and y-coordinate in a rectangular coordinate system;
S3: obtaining the intercept p and gradient q of the line;
wherein, NB(ti) represents an accumulated adsorption amount of ion Bm+, NB(ti) is the accumulated adsorption amount of the sample to ion Bm+ upon reaction for t=ti minutes, t0 is a pre-set time interval, and ti is a reaction time.

4. The system for determining and analyzing surface property parameters of a substance based on kinetic method according to claim 3, further comprising: a millivoltmeter, a constant flow pump, a liquid container, a constant temperature chamber and a thermostat; an input port of the millivoltmeter is connected with the detecting electrode; an output port of the millivoltmeter is connected with the concentration and activity operator; the constant flow pump is positioned between the reactor and the collector for liquid to be tested, for controlling liquid flow velocity from the reactor to the collector for liquid to be tested;
further comprising: a first control valve, a second control valve, a third control valve, an automatic controller for valve and a negative pressure pump; the first control valve is positioned between the collector for liquid to be tested and the container for liquid to be tested; the container for liquid to be tested is provided with the second control valve that is connected with the negative pressure pump; the container for liquid to be tested is further provided with the third control valve that is connected with the negative pressure pump; the third control valve and the negative pressure pump are provided with a waste liquid discharge outlet for draining liquid from the container for liquid to be tested; the automatic controller for valve is connected with the first control valve, the second control valve and the third control valve, respectively; the thermostat is connected with the constant flow pump; the liquid container is connected with the reactor; and, the sample processing system is positioned within the constant temperature chamber.

5. The system for determining and analyzing surface property parameters of a substance based on kinetic method according to claim 2, wherein: the surface charge total amount operator calculates the total amount of surface charge of the substance to be tested with the following equation:




  
    SCN
    =
    
      
        -
        m
      
      ⁢
      
        
          p
          0
        
        
          q
          0
        
      
    
  
  ,



wherein, SCN is the total amount of surface charge, m is valence of ion Bm+, p0 is intercept of the line, and q0 is gradient of the line.

6. The system for determining and analyzing surface property parameters of a substance based on kinetic method according to claim 2, wherein: the surface potential operator calculates the surface potential of the substance to be tested with the following equation:




  
    
      φ
      0
    
    =
    
      
        -
        
          RT
          
            
              (
              
                
                  m
                  ⁢
                  
                    
                  
                  ⁢
                  
                    β
                    B
                  
                
                -
                
                  n
                  ⁢
                  
                    
                  
                  ⁢
                  
                    β
                    A
                  
                
              
              )
            
            ⁢
            F
          
        
      
      ⁢
      
        
          ln
          [
          
            
              
                a
                A
                0
              
              
                a
                B
                0
              
            
            ⁢
            
              p
              
                (
                
                  
                    q
                    ⁢
                    
                      
                        
                        SCN
                        
                      
                      n
                    
                  
                  +
                  
                    
                      m
                      n
                    
                    ⁢
                    p
                  
                
                )
              
            
          
          ]
        
        2
      
    
  
  ,



wherein, φ0 is the surface potential of the substance to be tested, βB and βA are effective charge coefficients of ion Bm+ and ion An+, F represents Faraday constant, R represents gas constant, T represents temperature, m represents valence of ion Bm+, n represents valence of ion An+, p represents intercept, q represents gradient, aA0 and aB0 are respective activity of ion An+ and ion Bm+ in bulk solution when reaction reaches equilibrium;
the surface charge density operator calculates the surface charge density of the substance to be tested with the following equation:





  
    
      σ
      0
    
    ≈
    
      
        sign
        ⁡
        
          (
          
            φ
            0
          
          )
        
      
      ⁢
      
        
          
            
              ɛ
              ⁢
              
                
              
              ⁢
              RT
            
            
              2
              ⁢
              π
              ⁢
              
                
              
              ⁢
              
                F
                2
              
            
          
          ⁡
          
            [
            
              
                
                  a
                  B
                  0
                
                ⁡
                
                  (
                  
                    
                      e
                      
                        -
                        
                          
                            m
                            ⁢
                            
                              
                            
                            ⁢
                            
                              β
                              B
                            
                            ⁢
                            F
                            ⁢
                            
                              
                            
                            ⁢
                            
                              φ
                              0
                            
                          
                          RT
                        
                      
                    
                    -
                    1
                  
                  )
                
              
              +
              
                
                  a
                  A
                  0
                
                ⁡
                
                  (
                  
                    
                      e
                      
                        -
                        
                          
                            n
                            ⁢
                            
                              
                            
                            ⁢
                            
                              β
                              A
                            
                            ⁢
                            F
                            ⁢
                            
                              
                            
                            ⁢
                            
                              φ
                              0
                            
                          
                          RT
                        
                      
                    
                    -
                    1
                  
                  )
                
              
            
            ]
          
        
      
    
  
  ,



wherein, σ0 is the surface charge density, φ0 represents the surface potential, and ∈ is the medium dielectric constant of water;
the surface electric field intensity operator calculates the surface electric field intensity of the substance to be tested with the following equation:





  
    
      E
      0
    
    =
    
      
        
          4
          ⁢
          π
          ⁢
          
            
          
          ⁢
          F
        
        ɛ
      
      ⁢
      
        σ
        0
      
    
  
  ,



wherein, E0 represents the surface electric field intensity of the substance to be tested;
the specific surface area operator calculates the specific surface area of the substance to be tested with the following equation:





  
    S
    =
    
      SCN
      
        σ
        0
      
    
  
  ,



wherein, S represents the specific surface area of the substance to be tested.

7. A process according to the system for determining and analyzing surface property parameters of a substance based on kinetic method of claim 1, comprising the following steps:
S1: performing An+ saturation to the surface of the substance to be tested with AXn electrolyte solution having a known concentration;
wherein, the AXn electrolyte solution represents an electrolyte solution of metallic cation with a charge amount of n, and ion An+ represents a metallic cation with a charge amount of n;
S2: having a BXm electrolyte solution with a known concentration to flow through the surface of the substance to be tested, to have fusion replacement occurred between ion An+ and ion Bm+ on the surface of the substance to be tested;
wherein, BXm electrolyte solution represents an electrolyte solution of a metallic cation with a charge amount of m, and ion Bm+ represents a metallic cation with a charge amount of m;
S3: collecting liquid that flows out from the surface of the substance to be tested at the pre-set time interval;
S4: determining an activity of ion Bm+ in the flow-out liquid, and calculating an concentration of ion Bm+;
S5: using the concentration of ion Bm+ to calculate an accumulated adsorption amount of ion Bm+ according to the following equation:





  
    
      
        N
        B
      
      ⁡
      
        (
        
          t
          i
        
        )
      
    
    =
    
      
        1
        m
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          
            
              t
              i
            
            /
            
              t
              0
            
          
        
        ⁢
        
          V
          ⁡
          
            [
            
              
                f
                
                  B
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
              -
              
                
                  f
                  B
                
                ⁡
                
                  (
                  
                    t
                    i
                  
                  )
                
              
            
            ]
          
        
      
    
  
  ,



wherein, NB(ti) is the accumulated adsorption amount of the sample to ion Bm+ upon reaction for t=ti minutes, m is sample mass, V is an volume of liquid collected within the pre-set time interval, fB0 is a solution concentration before ion Bm+ flows into the sample, fB(ti) is concentration of ion Bm+ in liquid collected at ith time, t0 is the pre-set time interval, and ti is reaction time;
S6: using the accumulated adsorption amount NB(ti) of ion Bm+, with y-coordinate as y=[NB(ti+1)−NB(ti)]/(ti+1−ti) and corresponding x-coordinate as x=NB(ti)+0.5[NB(ti+1)−NB(ti)], to draw a line image based on positions of x-coordinate and y-coordinate in a rectangular coordinate system, to obtain intercept p0 and gradient q0 of the line;
S7: calculating the surface charge total amount of the sample based on intercept p0 and gradient q0 of the line according to the following equation:





  
    SCN
    =
    
      
        -
        m
      
      ⁢
      
        
          p
          0
        
        
          q
          0
        
      
    
  
  ,



wherein, SCN is the total amount of surface charge, and m is valence of ion Bm+.

8. The process of claim 7, wherein: step S4 further comprises the following steps:
S41: having BXm electrolyte solution with a known concentration to flow through the surface of the substance to be tested, to have fusion replacement occurred between ion An+ and ion Bm+ on the surface of the substance to be tested: mSoil−A+nBXm=nSoil−B+mAXn; obtaining a solution having An+, Bm+ and X−, and calculating concentration of ion An+ and concentration of ion Bm+:





  
    
      C
      A
    
    =
    
      
        (
        
          m
          n
        
        )
      
      ⁢
      
        C
        B
      
    
  
  ,



 CX=mnCB;
wherein, CA represents concentration of ion An+, CB represents concentration of ion Bm+, CX represents concentration of X−, Soil−A represents soil saturated by ion An+, and Soil−B represents soil saturated by ion Bm+;
S42: determining activity of ion Bm+ in aqueous solution, with a measurement of ion Bm+ activity as aB; using the ion activity as an initial concentration of ion Bm+ for a first iterative operation CB(0)=aB, and an initial concentration of ion An+ is:





  
    
      C
      A
      
        (
        0
        )
      
    
    =
    
      
        (
        
          m
          n
        
        )
      
      ⁢
      
        a
        B
      
    
  
  ;



S43: calculating concentration of ion X in the first iterative operation to be CX(0)=mnaB; and ion strength I(1) in the first iterative operation is:





  
    
      I
      
        (
        1
        )
      
    
    =
    
      
        
          1
          2
        
        ⁢
        
          
            ∑
            i
          
          ⁢
          
            
              c
              i
              
                (
                0
                )
              
            
            ⁢
            
              Z
              i
              2
            
          
        
      
      =
      
        
          
            1
            2
          
          ⁢
          
            (
            
              
                
                  m
                  2
                
                ⁢
                
                  a
                  B
                
              
              +
              
                mna
                B
              
              +
              
                mna
                B
              
            
            )
          
        
        =
        
          
            1
            2
          
          ⁢
          
            (
            
              
                m
                2
              
              +
              
                2
                ⁢
                mn
              
            
            )
          
          ⁢
          
            a
            B
          
        
      
    
  
  ,



wherein, I(1) represents ion strength coefficient in the first iterative operation, Zi is valence of ion i, ci(0) is the initial concentration of ion i in iterative operation;
S44: with Debye-Hückel limiting formula to obtain an initial activity coefficient of ion Bm+ as:





  
    
      γ
      B
      
        (
        1
        )
      
    
    =
    
      exp
      (
      
        -
        
          
            6030.2
            ×
            
              T
              
                -
                
                  3
                  2
                
              
            
            ⁢
            
              
              
                Z
                +
                2
              
              
            
            ⁢
            
              
                I
                
                  (
                  1
                  )
                
              
            
          
          
            1
            +
            
              
                I
                
                  (
                  1
                  )
                
              
            
          
        
      
      )
    
  
  ;



wherein, γB(1) represents initial activity coefficient of ion Bm+, Z+ represents charge amount positive ion, I(1) represents ion strength of the first iterative operation, and T represents temperature;
S45: calculating concentration cB(1) of ion Bm+ after the first iteration operation:





  
    
      c
      B
      
        (
        1
        )
      
    
    =
    
      
        a
        B
      
      
        γ
        B
        
          (
          1
          )
        
      
    
  
  ;



S46: calculating concentration of ion X after the first iterated operation as CX(1):


CX(1)=mncB(1);


S47: respectively calculating ion strength and activity coefficient of a second iteration as:




  
    
      
        I
        
          (
          2
          )
        
      
      =
      
        
          
            1
            2
          
          ⁢
          
            
              ∑
              i
            
            ⁢
            
              
                c
                i
                
                  (
                  1
                  )
                
              
              ⁢
              
                Z
                i
                2
              
            
          
        
        =
        
          
            1
            2
          
          ⁢
          
            (
            
              
                m
                2
              
              +
              
                2
                ⁢
                mn
              
            
            )
          
          ⁢
          
            c
            B
            
              (
              1
              )
            
          
        
      
    
    ;
  




  
    
      
        γ
        B
        
          (
          2
          )
        
      
      =
      
        exp
        (
        
          -
          
            
              6030.2
              ×
              
                T
                
                  -
                  
                    3
                    2
                  
                
              
              ⁢
              
                
                
                  Z
                  +
                  2
                
                
              
              ⁢
              
                
                  I
                  
                    (
                    2
                    )
                  
                
              
            
            
              1
              +
              
                
                  I
                  
                    (
                    2
                    )
                  
                
              
            
          
        
        )
      
    
    ;
  


wherein, γB(2) represents ion activity coefficient of the second iteration, cB(1) represents concentration of ion Bm+ after the first iteration, represents valence of ion i, and I(2) represents ion strength of the second iteration;
S48: calculating concentration cB(2) of ion Bm+ after the second iteration:





  
    
      c
      B
      
        (
        2
        )
      
    
    =
    
      
        a
        B
      
      
        γ
        B
        
          (
          2
          )
        
      
    
  
  ;



S49: concentration of ion X after the second iteration is:


CX(2)=mncB(2);


S410: repeating iteration for n times according to the steps S41-S49, and calculating in accordance with the following equation:




  
    
      
        I
        
          (
          n
          )
        
      
      =
      
        
          
            1
            2
          
          ⁢
          
            
              ∑
              i
            
            ⁢
            
              
                c
                i
                
                  (
                  
                    n
                    -
                    1
                  
                  )
                
              
              ⁢
              
                Z
                i
                2
              
            
          
        
        =
        
          
            1
            2
          
          ⁢
          
            (
            
              
                m
                2
              
              +
              
                2
                ⁢
                mn
              
            
            )
          
          ⁢
          
            c
            B
            
              (
              
                n
                -
                1
              
              )
            
          
        
      
    
    ;
  




  
    
      
        γ
        B
        
          (
          n
          )
        
      
      =
      
        exp
        (
        
          -
          
            
              6030.2
              ×
              
                T
                
                  -
                  
                    3
                    2
                  
                
              
              ⁢
              
                
                
                  Z
                  +
                  2
                
                
              
              ⁢
              
                
                  I
                  
                    (
                    n
                    )
                  
                
              
            
            
              1
              +
              
                
                  I
                  
                    (
                    n
                    )
                  
                
              
            
          
        
        )
      
    
    ;
  


when satisfying (I(n)−I(n−1))/I(n)<0.001, iteration is stopped, and a resulting concentration cB(n) is:





  
    
      c
      B
      
        (
        n
        )
      
    
    =
    
      
        a
        B
      
      
        γ
        B
        
          (
          n
          )
        
      
    
  
  ;



wherein, I(n) represents ion strength of nth iteration, I(n−1) represents ion strength of (n−1)th iteration, cB(n) represents concentration of ion Bm+ after nth iteration, and γB(n) represents ion activity coefficient after nth iteration.

9. The process of claim 7, further comprising the below steps following step S7:
S8: calculating ion strength I of the mixed electrolyte flow liquid containing ion An+ and ion Bm+ according to the following equation:


I=1⁄2[(n+n2)fA0+(m+m2)fB0]


wherein, fA0 is a solution concentration before ion An+ flows into the surface of the substance to be tested, and fB0 is a solution concentration before ion Bm+ flows into the surface of the substance to be tested;
S9: substituting the ion strength of the flow liquid in Davies equation, to calculate ion An+ activity coefficient γA0 and ion Bm+ activity coefficient γB0;
S10: respectively calculating activity of ion An+ and ion Bm+ according to the following equation:


aA0=fA0γA0 aB0=fB0γB0 


wherein, γA0 and γB0 are respective activity coefficient of ion An+ and ion Bm+ in bulk solution when reaction reaches equilibrium; and, aA0 and aB0 are respective activity of ion An+ and ion Bm+ in bulk solution when reaction reaches equilibrium.

10. The process of claim 9, further comprising the below steps following step S10:
S11: under a given pH condition, having a mixed electrolyte solution composed by AXn electrolyte solution with a known concentration fA0 and BXm electrolyte solution with a known concentration fB0 to flow through the surface of the substance to be tested, to have fusion replacement occurred between ion An+ and ion Bm+ on the surface of the substance to be tested, and collecting liquid that flows out from the surface of the substance to be tested at the pre-set time interval;
S12: repeating steps S1 to S6;
S13: using the accumulated adsorption amount NB(ti) of ion Bm+, with y-coordinate as y=[NB(ti+1)−NB(ti)]/(ti+1−ti) and corresponding x-coordinate as x=NB(ti)+0.5[NB(ti+1)−NB(ti)], to draw a line image based on positions of x-coordinate and y-coordinate in a rectangular coordinate system, to obtain intercept p and gradient q of the line.

11. The process of claim 10, further comprising the below steps following step S13:
S14: substituting the charge total amount obtained from step S7, the activity obtained from step S10, and the intercept p and gradient q of the line obtained from step S13 in the following equation to calculate the surface potential of the substance to be tested:





  
    
      φ
      0
    
    =
    
      
        -
        
          RT
          
            
              (
              
                
                  m
                  ⁢
                  
                    
                  
                  ⁢
                  
                    β
                    B
                  
                
                -
                
                  n
                  ⁢
                  
                    
                  
                  ⁢
                  
                    β
                    A
                  
                
              
              )
            
            ⁢
            F
          
        
      
      ⁢
      
        
          ln
          [
          
            
              
                a
                A
                0
              
              
                a
                B
                0
              
            
            ⁢
            
              p
              
                (
                
                  
                    q
                    ⁢
                    
                      
                        
                        SCN
                        
                      
                      n
                    
                  
                  +
                  
                    
                      m
                      n
                    
                    ⁢
                    p
                  
                
                )
              
            
          
          ]
        
        2
      
    
  
  ,



wherein, φ0 is the surface potential of the substance to be tested, βB and βA are effective charge coefficients of ion Bm+ and ion An+, F represents Faraday constant, R represents gas constant, and T represents temperature.

12. The process of claim 11, further comprising the below steps following step S14:
S15: substituting the surface potential of the substance to be tested in the following equation to calculate the surface charge density of the substance to be tested:





  
    
      σ
      0
    
    ≈
    
      
        sign
        ⁡
        
          (
          
            φ
            0
          
          )
        
      
      ⁢
      
        
          
            
              ɛ
              ⁢
              
                
              
              ⁢
              RT
            
            
              2
              ⁢
              π
              ⁢
              
                
              
              ⁢
              
                F
                2
              
            
          
          [
          
            
              
                a
                B
                0
              
              (
              
                
                  e
                  
                    -
                    
                      
                        m
                        ⁢
                        
                          
                        
                        ⁢
                        
                          β
                          B
                        
                        ⁢
                        F
                        ⁢
                        
                          
                        
                        ⁢
                        
                          φ
                          0
                        
                      
                      RT
                    
                  
                
                -
                1
              
              )
            
            +
            
              
                a
                A
                0
              
              (
              
                
                  e
                  
                    -
                    
                      
                        n
                        ⁢
                        
                          
                        
                        ⁢
                        
                          β
                          A
                        
                        ⁢
                        F
                        ⁢
                        
                          
                        
                        ⁢
                        
                          φ
                          0
                        
                      
                      RT
                    
                  
                
                -
                1
              
              )
            
          
          ]
        
      
    
  
  ,



wherein, σ0 is the surface charge density (mol/dm2), and ∈ is the medium dielectric constant, in which water has ∈=8.9×10−10 C2/(J·dm).

13. The process of claim 12, further comprising the below steps following step S15:
S16: substituting the surface charge density of the substance to be tested in the following equation to calculate the surface electric field intensity of the substance to be tested:





  
    
      E
      0
    
    =
    
      
        
          4
          ⁢
          π
          ⁢
          
            
          
          ⁢
          F
        
        ɛ
      
      ⁢
      
        σ
        0
      
    
  
  ,



wherein, E0 represents the surface electric field intensity of the substance to be tested.

14. The process of claim 13, further comprising the below steps following step S16:
S17: substituting the surface charge total amount and the surface charge density of the substance to be tested in the following equation to obtain the specific surface area of the sample:





  
    S
    =
    
      SCN
      
        σ
        0
      
    
  
  ,



wherein, S represents the specific surface area of the substance to be tested.

15. The process of claim 14, further comprising the below steps following step S17:
S18: outputting surface parameters of the substance to be tested, the surface parameters includes: surface potential of the substance to be tested, surface charge total amount of the substance to be tested, surface charge density of the substance to be tested, surface electric field intensity of the substance to be tested, and specific surface area of the substance to be tested.